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

#3089 Non-Darcy longitudinal : Use weighted average calculator

Browse Source
This commit is contained in:
Magne Sjaastad 2018-08-15 15:52:00 +02:00
parent 2e78bddaf6
commit fa25f448a1
2 changed files with 39 additions and 32 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ApplicationCode/Application/Tools/RiaWeightedAverageCalculator.cpp
View File

@ -18,6 +18,8 @@
#include "RiaWeightedAverageCalculator.h"
#include "cvfAssert.h"
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
@ -32,6 +34,8 @@ RiaWeightedAverageCalculator::RiaWeightedAverageCalculator()
//--------------------------------------------------------------------------------------------------
void RiaWeightedAverageCalculator::addValueAndWeight(double value, double weight)
{
CVF_ASSERT(weight >= 0.0);
m_aggregatedValue += value * weight;
m_aggregatedWeight += weight;
}

67
ApplicationCode/ProjectDataModel/Completions/RimStimPlanFractureTemplate.cpp
View File

@ -21,6 +21,7 @@
#include "RiaApplication.h"
#include "RiaFractureDefines.h"
#include "RiaLogging.h"
#include "RiaWeightedAverageCalculator.h"
#include "RifStimPlanXmlReader.h"
@ -457,51 +458,53 @@ FractureWidthAndConductivity
if (rimWellPath && rimWellPath->wellPathGeometry())
{
RigWellPathStimplanIntersector intersector(rimWellPath->wellPathGeometry(), fractureInstance);
double totalLength = 0.0;
for (const auto& v : intersector.intersections())
{
totalLength += v.second.computeLength();
}
std::vector<double> widthResultValues;
{
auto nameUnit = widthParameterNameAndUnit();
widthResultValues = fractureGridResultsForUnitSystem(
nameUnit.first, nameUnit.second, m_activeTimeStepIndex, fractureTemplateUnit());
}
std::vector<double> conductivityResultValues;
{
auto nameUnit = conductivityParameterNameAndUnit();
conductivityResultValues = fractureGridResultsForUnitSystem(
nameUnit.first, nameUnit.second, m_activeTimeStepIndex, fractureTemplateUnit());
}
double totalLength = 0.0;
double weightedConductivity = 0.0;
double weightedWidth = 0.0;
for (const auto& v : intersector.intersections())
{
size_t fractureGlobalCellIndex = v.first;
if (fractureGlobalCellIndex < widthResultValues.size())
std::vector<double> widthResultValues;
{
weightedWidth += widthResultValues[fractureGlobalCellIndex] * v.second.computeLength();
auto nameUnit = widthParameterNameAndUnit();
widthResultValues = fractureGridResultsForUnitSystem(
nameUnit.first, nameUnit.second, m_activeTimeStepIndex, fractureTemplateUnit());
}
if (fractureGlobalCellIndex < conductivityResultValues.size())
std::vector<double> conductivityResultValues;
{
weightedConductivity += conductivityResultValues[fractureGlobalCellIndex] * v.second.computeLength();
auto nameUnit = conductivityParameterNameAndUnit();
conductivityResultValues = fractureGridResultsForUnitSystem(
nameUnit.first, nameUnit.second, m_activeTimeStepIndex, fractureTemplateUnit());
}
RiaWeightedAverageCalculator widthCalc;
RiaWeightedAverageCalculator conductivityCalc;
RigWellPathStimplanIntersector intersector(rimWellPath->wellPathGeometry(), fractureInstance);
for (const auto& v : intersector.intersections())
{
size_t fractureGlobalCellIndex = v.first;
double intersectionLength = v.second.computeLength();
if (fractureGlobalCellIndex < widthResultValues.size())
{
widthCalc.addValueAndWeight(widthResultValues[fractureGlobalCellIndex], intersectionLength);
}
if (fractureGlobalCellIndex < conductivityResultValues.size())
{
conductivityCalc.addValueAndWeight(conductivityResultValues[fractureGlobalCellIndex],
intersectionLength);
}
}
weightedConductivity = conductivityCalc.weightedAverage();
weightedWidth = widthCalc.weightedAverage();
totalLength = widthCalc.aggregatedWeight();
}
if (totalLength > 1e-7)
{
weightedWidth /= totalLength;
weightedConductivity /= totalLength;
values.m_width = weightedWidth;
values.m_conductivity = weightedConductivity;
}