OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#5747 Completion Export: Avoid nan values for zero permeability

Browse Source
This commit is contained in:
Magne Sjaastad 2020-04-03 10:35:09 +02:00 committed by Gaute Lindkvist
parent 66c3587d5b
commit efe79aeacf
1 changed files with 7 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
ApplicationCode/ReservoirDataModel/RigTransmissibilityEquations.cpp
View File

@ -35,6 +35,13 @@ double RigTransmissibilityEquations::wellBoreTransmissibilityComponent( double c
{ {
double K = cvf::Math::sqrt( permeabilityNormalDirection1 * permeabilityNormalDirection2 ); double K = cvf::Math::sqrt( permeabilityNormalDirection1 * permeabilityNormalDirection2 );
const double lowerLimit = 1.0e-9;
if ( std::fabs( permeabilityNormalDirection1 * permeabilityNormalDirection2 ) < lowerLimit )
{
// Guard further computations to avoid nan values
return 0.0;
}
double nominator = cDarcyForRelevantUnit * 2 * cvf::PI_D * K * cellPerforationVectorComponent; double nominator = cDarcyForRelevantUnit * 2 * cvf::PI_D * K * cellPerforationVectorComponent;
double peaceManRad = peacemanRadius( permeabilityNormalDirection1, double peaceManRad = peacemanRadius( permeabilityNormalDirection1,