#2627 Export CF : kh is reported as zero

This commit is contained in:
Magne Sjaastad 2018-03-20 10:48:06 +01:00
parent cb6dbb1bee
commit 39c8834181
3 changed files with 14 additions and 0 deletions

View File

@ -611,6 +611,16 @@ double RimFractureTemplate::dFactor() const
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
double RimFractureTemplate::kh() const double RimFractureTemplate::kh() const
{ {
// kh = permeability * h
// conductivity = permeability * h
auto values = widthAndConductivityAtWellPathIntersection();
if (values.m_conductivity != HUGE_VAL)
{
// If conductivity is found in stim plan file, use this directly
return values.m_conductivity;
}
return effectivePermeability() * fractureWidth(); return effectivePermeability() * fractureWidth();
} }

View File

@ -44,6 +44,7 @@ public:
FractureWidthAndConductivity() FractureWidthAndConductivity()
: m_width(0.0) : m_width(0.0)
, m_permeability(0.0) , m_permeability(0.0)
, m_conductivity(HUGE_VAL)
{ {
} }
@ -60,6 +61,8 @@ public:
// Unit : mD (milliDarcy) // Unit : mD (milliDarcy)
double m_permeability; double m_permeability;
double m_conductivity;
}; };
//================================================================================================== //==================================================================================================

View File

@ -403,6 +403,7 @@ FractureWidthAndConductivity RimStimPlanFractureTemplate::widthAndConductivityAt
const RigFractureCell& wellCell = m_fractureGrid->cellFromIndex(wellCellIndex); const RigFractureCell& wellCell = m_fractureGrid->cellFromIndex(wellCellIndex);
double conductivity = wellCell.getConductivtyValue(); double conductivity = wellCell.getConductivtyValue();
values.m_conductivity = conductivity;
std::vector<std::pair<QString, QString> > propertyNamesUnitsOnFile = m_stimPlanFractureDefinitionData->getStimPlanPropertyNamesUnits(); std::vector<std::pair<QString, QString> > propertyNamesUnitsOnFile = m_stimPlanFractureDefinitionData->getStimPlanPropertyNamesUnits();