OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-08 23:23:01 -06:00
Code Issues Packages Projects Releases Wiki Activity

StimPlan Model Calculator: guard against inf/nan in stress calculation.

Browse Source
This commit is contained in:
Kristian Bendiksen 2021-05-19 15:02:47 +02:00 committed by Magne Sjaastad
parent 3b383f1a57
commit f72890584f
1 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
View File

@ -436,6 +436,11 @@ bool RimStimPlanModelCalculator::calculateStressWithGradients( std::vector<doubl
double Sv_diff = 0.0; double Sv_diff = 0.0;
double deltaHorizontalStress = poissonsRatio / ( 1.0 - poissonsRatio ) * ( Sv_diff - biot * pressureDiff ) + double deltaHorizontalStress = poissonsRatio / ( 1.0 - poissonsRatio ) * ( Sv_diff - biot * pressureDiff ) +
( biot * pressureDiff ); ( biot * pressureDiff );
if ( std::isnan( deltaHorizontalStress ) || std::isinf( deltaHorizontalStress ) )
{
RiaLogging::warning( "Invalid horizontal stress delta calculated. Setting to zero." );
deltaHorizontalStress = 0.0;
}
double depletionStress = Sh_init + deltaHorizontalStress; double depletionStress = Sh_init + deltaHorizontalStress;
stress.push_back( RiaEclipseUnitTools::barToPsi( depletionStress ) ); stress.push_back( RiaEclipseUnitTools::barToPsi( depletionStress ) );