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

Thermal Fracture: compute filter cake mobility.

Browse Source
This commit is contained in:
Kristian Bendiksen 2022-09-23 15:48:27 +02:00
parent ce65d3738a
commit 7d774d733c
1 changed files with 26 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
ApplicationLibCode/ProjectDataModel/Completions/RimThermalFractureTemplate.cpp
View File

@ -172,7 +172,33 @@ void RimThermalFractureTemplate::loadDataAndUpdate()
} }
}; };
auto addFilterCakeMobility = []( std::shared_ptr<RigThermalFractureDefinition> def ) {
int leakoffPressureDropIndex = def->getPropertyIndex( "LeakoffPressureDrop" );
int filtratePressureDropIndex = def->getPropertyIndex( "FiltratePressureDrop" );
int leakoffMobilityIndex = def->getPropertyIndex( "LeakoffMobility" );
QString filterCakeMobilityValueTag = "FilterCakeMobility";
def->addProperty( filterCakeMobilityValueTag, "m/day/bar" );
int filterCakeMobilityDeclineIndex = def->getPropertyIndex( filterCakeMobilityValueTag );
for ( size_t nodeIndex = 0; nodeIndex < def->numNodes(); nodeIndex++ )
{
for ( size_t timeStepIndex = 0; timeStepIndex < def->numTimeSteps(); timeStepIndex++ )
{
double leakoffPressureDrop =
def->getPropertyValue( leakoffPressureDropIndex, nodeIndex, timeStepIndex );
double filtratePressureDrop =
def->getPropertyValue( filtratePressureDropIndex, nodeIndex, timeStepIndex );
double leakoffMobility = def->getPropertyValue( leakoffMobilityIndex, nodeIndex, timeStepIndex );
double filterCakeMobilityValue = leakoffMobility * ( leakoffPressureDrop / filtratePressureDrop );
def->appendPropertyValue( filterCakeMobilityDeclineIndex, nodeIndex, filterCakeMobilityValue );
}
}
};
addInjectivityDecline( m_fractureDefinitionData ); addInjectivityDecline( m_fractureDefinitionData );
addFilterCakeMobility( m_fractureDefinitionData );
setDefaultConductivityResultIfEmpty(); setDefaultConductivityResultIfEmpty();