OilfieldToolsNet/ResInsight
4
0 Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2026-10-30 09:00:20 -05:00
Code Issues Packages Projects Releases Wiki Activity

#12302 Well Target Mapping: use P90 as default + "Reset to P90"-button

Browse Source
This commit is contained in:
Kristian Bendiksen authored and Magne Sjaastad committed 2025-03-27 13:28:10 +01:00
1 parent 59f9104352
commit fba8abbc9f
2 files changed
+57 -12
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

ApplicationLibCode/ProjectDataModel/RimWellTargetMapping.cpp
+51 -12
View File
@@ -27,6 +27,7 @@
#include "RigCaseCellResultsData.h"
#include "RigEclipseResultAddress.h"
#include "RigStatisticsMath.h"
#include "Well/RigWellTargetMapping.h"
#include "RimEclipseCase.h"
@@ -102,6 +103,10 @@ RimWellTargetMapping::RimWellTargetMapping()
CAF_PDM_InitField( &m_transmissibility, "Transmissibility", 0.0, "Transmissibility" );
m_transmissibility.uiCapability()->setUiEditorTypeName( caf::PdmUiDoubleSliderEditor::uiEditorTypeName() );
CAF_PDM_InitField( &m_resetDefaultButton, "ResetDefaultButton", true, "Reset to P90" );
caf::PdmUiPushButtonEditor::configureEditorLabelHidden( &m_resetDefaultButton );
m_resetDefaultButton.xmlCapability()->disableIO();
CAF_PDM_InitField( &m_maxIterations, "Iterations", 100000, "Max Iterations" );
CAF_PDM_InitField( &m_maxNumTargets, "MaxNumTargets", 5, "Maximum Number of Well Targets" );
@@ -124,12 +129,15 @@ RimWellTargetMapping::RimWellTargetMapping()
m_minimumPressure = cvf::UNDEFINED_DOUBLE;
m_maximumPressure = cvf::UNDEFINED_DOUBLE;
m_defaultPressure = cvf::UNDEFINED_DOUBLE;
m_minimumPermeability = cvf::UNDEFINED_DOUBLE;
m_maximumPermeability = cvf::UNDEFINED_DOUBLE;
m_defaultPermeability = cvf::UNDEFINED_DOUBLE;
m_minimumTransmissibility = cvf::UNDEFINED_DOUBLE;
m_maximumTransmissibility = cvf::UNDEFINED_DOUBLE;
m_defaultTransmissibility = cvf::UNDEFINED_DOUBLE;
setDeletable( true );
}
@@ -169,6 +177,10 @@ void RimWellTargetMapping::fieldChangedByUi( const caf::PdmFieldHandle* changedF
}
}
}
else if ( changedField == &m_resetDefaultButton )
{
resetMinimumCellValuesToDefault();
}
}
//--------------------------------------------------------------------------------------------------
@@ -212,10 +224,11 @@ void RimWellTargetMapping::updateAllBoundaries()
const int timeStepIdx = m_timeStep();
auto updateBoundaryValues =
[]( auto resultsData, const std::vector<RigEclipseResultAddress>& addresses, size_t timeStepIdx ) -> std::pair<double, double>
[]( auto resultsData, const std::vector<RigEclipseResultAddress>& addresses, size_t timeStepIdx ) -> std::tuple<double, double, double>
{
double globalMin = std::numeric_limits<double>::max();
double globalMax = -std::numeric_limits<double>::max();
double globalMin = std::numeric_limits<double>::max();
double globalMax = -std::numeric_limits<double>::max();
std::vector<double> allValues;
for ( auto address : addresses )
{
double currentMinimum;
@@ -223,27 +236,31 @@ void RimWellTargetMapping::updateAllBoundaries()
if ( resultsData->ensureKnownResultLoaded( address ) )
{
resultsData->minMaxCellScalarValues( address, timeStepIdx, currentMinimum, currentMaximum );
globalMin = std::min( globalMin, currentMinimum );
globalMax = std::max( globalMax, currentMaximum );
globalMin = std::min( globalMin, currentMinimum );
globalMax = std::max( globalMax, currentMaximum );
const std::vector<double>& values = resultsData->cellScalarResults( address, timeStepIdx );
allValues.insert( allValues.end(), values.begin(), values.end() );
}
}
return { globalMin, globalMax };
double p10, p50, p90, mean;
RigStatisticsMath::calculateStatisticsCurves( allValues, &p10, &p50, &p90, &mean, RigStatisticsMath::PercentileStyle::SWITCHED );
return { globalMin, globalMax, p90 };
};
std::tie( m_minimumPressure, m_maximumPressure ) =
std::tie( m_minimumPressure, m_maximumPressure, m_defaultPressure ) =
updateBoundaryValues( resultsData, { RigEclipseResultAddress( RiaDefines::ResultCatType::DYNAMIC_NATIVE, "PRESSURE" ) }, timeStepIdx );
std::vector<double> volume =
RigWellTargetMapping::getVolumeVector( *resultsData, m_volumeType(), m_volumesType(), m_volumeResultType(), timeStepIdx );
std::tie( m_minimumPermeability, m_maximumPermeability ) =
std::tie( m_minimumPermeability, m_maximumPermeability, m_defaultPermeability ) =
updateBoundaryValues( resultsData,
{ RigEclipseResultAddress( RiaDefines::ResultCatType::STATIC_NATIVE, "PERMX" ),
RigEclipseResultAddress( RiaDefines::ResultCatType::STATIC_NATIVE, "PERMY" ),
RigEclipseResultAddress( RiaDefines::ResultCatType::STATIC_NATIVE, "PERMZ" ) },
0 );
std::tie( m_minimumTransmissibility, m_maximumTransmissibility ) =
std::tie( m_minimumTransmissibility, m_maximumTransmissibility, m_defaultTransmissibility ) =
updateBoundaryValues( resultsData,
{ RigEclipseResultAddress( RiaDefines::ResultCatType::STATIC_NATIVE, "TRANX" ),
RigEclipseResultAddress( RiaDefines::ResultCatType::STATIC_NATIVE, "TRANY" ),
@@ -294,6 +311,14 @@ void RimWellTargetMapping::defineEditorAttribute( const caf::PdmFieldHandle* fie
attrib->m_buttonText = "Generate";
}
}
if ( field == &m_resetDefaultButton )
{
if ( auto attrib = dynamic_cast<caf::PdmUiPushButtonEditorAttribute*>( attribute ) )
{
attrib->m_buttonText = "Reset to P90";
}
}
}
//--------------------------------------------------------------------------------------------------
@@ -365,6 +390,7 @@ void RimWellTargetMapping::defineUiOrdering( QString uiConfigName, caf::PdmUiOrd
minimumCellValuesGroup->add( &m_pressure );
minimumCellValuesGroup->add( &m_permeability );
minimumCellValuesGroup->add( &m_transmissibility );
minimumCellValuesGroup->add( &m_resetDefaultButton );
auto hasEnsembleParent = firstAncestorOrThisOfType<RimEclipseCaseEnsemble>() != nullptr;
@@ -468,9 +494,22 @@ void RimWellTargetMapping::setDefaults()
RigWellTargetMapping::VolumesType::COMPUTED_VOLUMES };
std::vector<RigWellTargetMapping::VolumesType> availableVolumesTypes = findAvailableVolumesTypes( eclipseCase );
m_volumesType = findFirstByPriority( volumesTypesByPriority, availableVolumesTypes );
updateAllBoundaries();
resetMinimumCellValuesToDefault();
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RimWellTargetMapping::resetMinimumCellValuesToDefault()
{
m_pressure = std::clamp( m_defaultPressure, m_minimumPressure, m_maximumPressure );
m_permeability = std::clamp( m_defaultPermeability, m_minimumPermeability, m_maximumPermeability );
m_transmissibility = std::clamp( m_defaultTransmissibility, std::max( m_minimumTransmissibility, 0.1 ), m_maximumTransmissibility );
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
ApplicationLibCode/ProjectDataModel/RimWellTargetMapping.h
+6
View File
@@ -65,6 +65,8 @@ private:
RigWellTargetMapping::ClusteringLimits getClusteringLimits() const;
void resetMinimumCellValuesToDefault();
caf::PdmField<int> m_timeStep;
caf::PdmField<caf::AppEnum<RigWellTargetMapping::VolumeType>> m_volumeType;
@@ -74,6 +76,7 @@ private:
caf::PdmField<double> m_pressure;
caf::PdmField<double> m_permeability;
caf::PdmField<double> m_transmissibility;
caf::PdmField<bool> m_resetDefaultButton;
caf::PdmField<int> m_maxIterations;
caf::PdmField<int> m_maxNumTargets;
@@ -90,10 +93,13 @@ private:
double m_minimumPressure;
double m_maximumPressure;
double m_defaultPressure;
double m_minimumPermeability;
double m_maximumPermeability;
double m_defaultPermeability;
double m_minimumTransmissibility;
double m_maximumTransmissibility;
double m_defaultTransmissibility;
};