mirror of
https://github.com/OPM/ResInsight.git
synced 2026-10-26 15:10:22 -05:00
Magne Sjaastad
93bd0b9c9d
* Refactor well log track settings into modular classes Major refactoring of RimWellLogTrack to delegate property axis, formation, region annotation, and well path attribute settings to dedicated settings classes. Introduces RimWellLogPropertyAxisSettings, RimWellLogFormationSettings, RimWellLogRegionAnnotationSettings, and RimWellLogWellPathAttributeSettings, each encapsulating related fields, UI logic, and option handling. Old fields are retained for migration and backward compatibility. Updates enums and signatures for consistency.
194 lines
8.3 KiB
C++
194 lines
8.3 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) 2021- Equinor ASA
|
|
//
|
|
// ResInsight is free software: you can redistribute it and/or modify
|
|
// it under the terms of the GNU General Public License as published by
|
|
// the Free Software Foundation, either version 3 of the License, or
|
|
// (at your option) any later version.
|
|
//
|
|
// ResInsight is distributed in the hope that it will be useful, but WITHOUT
|
|
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
// FITNESS FOR A PARTICULAR PURPOSE.
|
|
//
|
|
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
|
|
// for more details.
|
|
//
|
|
/////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#include "RimEnsembleWellLogStatisticsCurve.h"
|
|
|
|
#include "RiaApplication.h"
|
|
#include "RiaDefines.h"
|
|
#include "RiaInterpolationTools.h"
|
|
#include "RiaLogging.h"
|
|
|
|
#include "RimEnsembleWellLogCurveSet.h"
|
|
#include "RimWellLogTrack.h"
|
|
#include "RimWellPath.h"
|
|
|
|
#include "RiuQwtPlotCurve.h"
|
|
#include "RiuQwtPlotWidget.h"
|
|
|
|
#include "cafPdmUiTreeOrdering.h"
|
|
|
|
CAF_PDM_SOURCE_INIT( RimEnsembleWellLogStatisticsCurve, "EnsembleWellLogStatisticsCurve" );
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RimEnsembleWellLogStatisticsCurve::RimEnsembleWellLogStatisticsCurve()
|
|
{
|
|
CAF_PDM_InitObject( "Ensemble Well Log Statistics Curve" );
|
|
|
|
CAF_PDM_InitFieldNoDefault( &m_ensembleWellLogCurveSet, "EnsembleWellLogCurveSet", "Ensemble Well Log Curve Set" );
|
|
m_ensembleWellLogCurveSet.uiCapability()->setUiTreeChildrenHidden( true );
|
|
|
|
CAF_PDM_InitFieldNoDefault( &m_statisticsType, "StatisticsType", "Statistics Type" );
|
|
m_statisticsType.uiCapability()->setUiHidden( true );
|
|
|
|
m_wellPath = nullptr;
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RimEnsembleWellLogStatisticsCurve::~RimEnsembleWellLogStatisticsCurve()
|
|
{
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
void RimEnsembleWellLogStatisticsCurve::setEnsembleWellLogCurveSet( RimEnsembleWellLogCurveSet* ensembleWellLogCurveSet )
|
|
{
|
|
m_ensembleWellLogCurveSet = ensembleWellLogCurveSet;
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
void RimEnsembleWellLogStatisticsCurve::setStatisticsType( RimEnsembleWellLogStatistics::StatisticsType statisticsType )
|
|
{
|
|
m_statisticsType = statisticsType;
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RimEnsembleWellLogStatistics::StatisticsType RimEnsembleWellLogStatisticsCurve::statisticsType() const
|
|
{
|
|
return m_statisticsType();
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
void RimEnsembleWellLogStatisticsCurve::performDataExtraction( bool* isUsingPseudoLength )
|
|
{
|
|
std::vector<double> values;
|
|
std::vector<double> measuredDepthValues;
|
|
std::vector<double> tvDepthValues;
|
|
double rkbDiff = 0.0;
|
|
|
|
*isUsingPseudoLength = false;
|
|
|
|
if ( m_ensembleWellLogCurveSet )
|
|
{
|
|
const RimEnsembleWellLogStatistics* ensembleWellLogStatistics = m_ensembleWellLogCurveSet->ensembleWellLogStatistics();
|
|
|
|
m_ensembleWellLogCurveSet->updateStatistics();
|
|
|
|
RiaDefines::DepthUnitType depthUnit = ensembleWellLogStatistics->depthUnitType();
|
|
QString xUnits = ensembleWellLogStatistics->logChannelUnitString();
|
|
|
|
if ( m_statisticsType == RimEnsembleWellLogStatistics::StatisticsType::MEAN )
|
|
{
|
|
values = ensembleWellLogStatistics->mean();
|
|
measuredDepthValues = ensembleWellLogStatistics->measuredDepths();
|
|
tvDepthValues = ensembleWellLogStatistics->tvDepths();
|
|
}
|
|
else if ( m_statisticsType == RimEnsembleWellLogStatistics::StatisticsType::P10 )
|
|
{
|
|
values = ensembleWellLogStatistics->p10();
|
|
measuredDepthValues = ensembleWellLogStatistics->measuredDepths();
|
|
tvDepthValues = ensembleWellLogStatistics->tvDepths();
|
|
}
|
|
else if ( m_statisticsType == RimEnsembleWellLogStatistics::StatisticsType::P50 )
|
|
{
|
|
values = ensembleWellLogStatistics->p50();
|
|
measuredDepthValues = ensembleWellLogStatistics->measuredDepths();
|
|
tvDepthValues = ensembleWellLogStatistics->tvDepths();
|
|
}
|
|
else if ( m_statisticsType == RimEnsembleWellLogStatistics::StatisticsType::P90 )
|
|
{
|
|
values = ensembleWellLogStatistics->p90();
|
|
measuredDepthValues = ensembleWellLogStatistics->measuredDepths();
|
|
tvDepthValues = ensembleWellLogStatistics->tvDepths();
|
|
}
|
|
|
|
if ( !values.empty() && !measuredDepthValues.empty() && measuredDepthValues.size() == values.size() )
|
|
{
|
|
if ( m_ensembleWellLogCurveSet->depthEqualization() == RimEnsembleWellLogStatistics::DepthEqualization::NONE )
|
|
{
|
|
std::vector<double> dummyValues( values.begin(), values.end() );
|
|
addDatapointsForBottomOfSegment( measuredDepthValues, values );
|
|
// Pass copy of values vector to avoid adding bottom segments twice
|
|
addDatapointsForBottomOfSegment( tvDepthValues, dummyValues );
|
|
}
|
|
|
|
std::map<RiaDefines::DepthType, std::vector<double>> validDepths;
|
|
if ( values.size() == measuredDepthValues.size() )
|
|
{
|
|
validDepths.insert( std::make_pair( RiaDefines::DepthType::MEASURED_DEPTH, measuredDepthValues ) );
|
|
}
|
|
if ( values.size() == tvDepthValues.size() )
|
|
{
|
|
validDepths.insert( std::make_pair( RiaDefines::DepthType::TRUE_VERTICAL_DEPTH, tvDepthValues ) );
|
|
validDepths.insert( std::make_pair( RiaDefines::DepthType::TRUE_VERTICAL_DEPTH_RKB, tvDepthValues ) );
|
|
}
|
|
|
|
bool useLogarithmicScale = false;
|
|
setPropertyValuesAndDepths( values, validDepths, rkbDiff, depthUnit, false, useLogarithmicScale );
|
|
}
|
|
}
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
QString RimEnsembleWellLogStatisticsCurve::createCurveAutoName()
|
|
{
|
|
return caf::AppEnum<RimEnsembleWellLogStatistics::StatisticsType>::uiText( m_statisticsType() );
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
void RimEnsembleWellLogStatisticsCurve::addDatapointsForBottomOfSegment( std::vector<double>& depthValues, std::vector<double>& values )
|
|
{
|
|
std::vector<double> depthValuesWithBottomLayers;
|
|
std::vector<double> valuesWithBottomLayers;
|
|
for ( size_t i = 0; i < values.size(); i++ )
|
|
{
|
|
// Add the data point at top of the layer
|
|
double topLayerDepth = depthValues[i];
|
|
double value = values[i];
|
|
depthValuesWithBottomLayers.push_back( topLayerDepth );
|
|
valuesWithBottomLayers.push_back( value );
|
|
|
|
// Add extra data points for bottom part of the layer
|
|
if ( i < values.size() - 1 )
|
|
{
|
|
double bottomLayerDepth = depthValues[i + 1];
|
|
double bottomValue = value;
|
|
|
|
depthValuesWithBottomLayers.push_back( bottomLayerDepth );
|
|
valuesWithBottomLayers.push_back( bottomValue );
|
|
}
|
|
}
|
|
|
|
values = valuesWithBottomLayers;
|
|
depthValues = depthValuesWithBottomLayers;
|
|
}
|