From 9dc6f651058ffc8a1413791159e42a2c09ca07ee Mon Sep 17 00:00:00 2001 From: Kristian Bendiksen Date: Wed, 20 Dec 2023 14:13:12 +0100 Subject: [PATCH] Fault reactivation: Use same pore pressure gradient for POR and stress. --- .../Faults/RimFaultReactivationDataAccess.cpp | 6 +++--- .../Faults/RimFaultReactivationDataAccessorStress.cpp | 4 ++-- 2 files changed, 5 insertions(+), 5 deletions(-) diff --git a/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp b/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp index f2e7691c0f..b095b21731 100644 --- a/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp +++ b/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccess.cpp @@ -48,8 +48,8 @@ RimFaultReactivationDataAccess::RimFaultReactivationDataAccess( RimEclipseCase* const std::vector& timeSteps ) : m_timeSteps( timeSteps ) { - // TODO: correct default pore pressure gradient? - m_accessors.push_back( std::make_shared( thecase, 1.0 ) ); + double porePressureGradient = 1.0; + m_accessors.push_back( std::make_shared( thecase, porePressureGradient ) ); m_accessors.push_back( std::make_shared( thecase, 0.0001 ) ); m_accessors.push_back( std::make_shared( thecase ) ); if ( geoMechCase ) @@ -71,7 +71,7 @@ RimFaultReactivationDataAccess::RimFaultReactivationDataAccess( RimEclipseCase* RimFaultReactivation::Property::LateralStressComponentY }; for ( auto property : stressProperties ) { - m_accessors.push_back( std::make_shared( geoMechCase, property, 0.003 ) ); + m_accessors.push_back( std::make_shared( geoMechCase, property, porePressureGradient ) ); } } } diff --git a/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStress.cpp b/ApplicationLibCode/ProjectDataModel/Faults/RimFaultReactivationDataAccessorStress.cpp index 67559300bf..5c4a633e60 100644 --- 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( 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(),