mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
#3625 Fracture header : Add computation of permeability
This commit is contained in:
@@ -22,6 +22,7 @@
|
||||
#include "cvfMath.h"
|
||||
|
||||
#include <cmath>
|
||||
#include <limits>
|
||||
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
///
|
||||
@@ -70,6 +71,25 @@ double RigTransmissibilityEquations::totalPermeability(const double permx, const
|
||||
return totalPerm;
|
||||
}
|
||||
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
///
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
double RigTransmissibilityEquations::permeability(const double conductivity, const double width)
|
||||
{
|
||||
double threshold = 1e-7;
|
||||
|
||||
if (std::fabs(width) > threshold)
|
||||
{
|
||||
double perm = conductivity / width;
|
||||
|
||||
return perm;
|
||||
}
|
||||
else
|
||||
{
|
||||
return 0.0;
|
||||
}
|
||||
}
|
||||
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
///
|
||||
//--------------------------------------------------------------------------------------------------
|
||||
|
||||
@@ -36,6 +36,8 @@ public:
|
||||
|
||||
static double totalPermeability(const double permx, const double permy, const double permz);
|
||||
|
||||
static double permeability(const double conductivity, const double width);
|
||||
|
||||
private:
|
||||
// If using wellBoreTransmissibilityComponent to calculate Tx (transmissibility in x direction),
|
||||
// perforationVectorComponent is the x component (in the cell local coordinate system) of the perforation vector
|
||||
|
||||
Reference in New Issue
Block a user