Fault reactivation: Use same pore pressure gradient for POR and stress.

OPM · Dec 20, 2023 · 9dc6f65 · 9dc6f65
1 parent 29890ee
commit 9dc6f65
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp b/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp
@@ -48,8 +48,8 @@ RimFaultReactivationDataAccess::RimFaultReactivationDataAccess( RimEclipseCase*
                                                                 const std::vector<size_t>& timeSteps )
     : m_timeSteps( timeSteps )
 {
-    // TODO: correct default pore pressure gradient?
-    m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorPorePressure>( thecase, 1.0 ) );
+    double porePressureGradient = 1.0;
+    m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorPorePressure>( thecase, porePressureGradient ) );
     m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorVoidRatio>( thecase, 0.0001 ) );
     m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorTemperature>( thecase ) );
     if ( geoMechCase )
@@ -71,7 +71,7 @@ RimFaultReactivationDataAccess::RimFaultReactivationDataAccess( RimEclipseCase*
                                                                          RimFaultReactivation::Property::LateralStressComponentY };
         for ( auto property : stressProperties )
         {
-            m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorStress>( geoMechCase, property, 0.003 ) );
+            m_accessors.push_back( std::make_shared<RimFaultReactivationDataAccessorStress>( geoMechCase, property, porePressureGradient ) );
         }
     }
 }

diff --git a/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStress.cpp b/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStress.cpp
@@ -381,7 +381,7 @@ void RimFaultReactivationDataAccessorStress::fillInMissingValues( const std::vec
     auto [lastOverburdenIndex, firstUnderburdenIndex] = findOverburdenAndUnderburdenIndex( values );
 
     // Fill in overburden values using gradient
-    double topPorePressure = calculatePorePressure( std::abs( intersections[0].z() ), 1.0 );
+    double topPorePressure = calculatePorePressure( std::abs( intersections[0].z() ), gradient );
     double overburdenGradient =
         computeGradient( intersections[0].z(), topPorePressure, intersections[lastOverburdenIndex].z(), values[lastOverburdenIndex] );
 
@@ -392,7 +392,7 @@ void RimFaultReactivationDataAccessorStress::fillInMissingValues( const std::vec
 
     // Fill in underburden values using gradient
     int    lastElementIndex    = static_cast<int>( values.size() ) - 1;
-    double bottomPorePressure  = calculatePorePressure( std::abs( intersections[lastElementIndex].z() ), 1.0 );
+    double bottomPorePressure  = calculatePorePressure( std::abs( intersections[lastElementIndex].z() ), gradient );
     double underburdenGradient = computeGradient( intersections[firstUnderburdenIndex].z(),
                                                   values[firstUnderburdenIndex],
                                                   intersections[lastElementIndex].z(),