mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
18c29d5fc7
The "absolute" filter cake adjustment should only be used for cells when filter cake is present (i.e filtrate thickness is not zero).
85 lines
4.0 KiB
C++
85 lines
4.0 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) 2022 Statoil 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 "RigEclipseToThermalCellTransmissibilityCalculator.h"
|
|
|
|
#include "RimEclipseCase.h"
|
|
#include "RimFracture.h"
|
|
#include "RimThermalFractureTemplate.h"
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RigEclipseToThermalCellTransmissibilityCalculator::RigEclipseToThermalCellTransmissibilityCalculator(
|
|
const RimEclipseCase* caseToApply,
|
|
cvf::Mat4d fractureTransform,
|
|
double skinFactor,
|
|
double cDarcy,
|
|
const RigFractureCell& stimPlanCell,
|
|
const RimFracture* fracture,
|
|
RimThermalFractureTemplate::FilterCakePressureDrop filterCakePressureDrop,
|
|
double injectivityFactor,
|
|
double filterCakeMobility,
|
|
double filtrateThickness,
|
|
double viscosity,
|
|
double relativePermeability )
|
|
: RigEclipseToStimPlanCellTransmissibilityCalculator( caseToApply, fractureTransform, skinFactor, cDarcy, stimPlanCell, fracture )
|
|
, m_filterCakePressureDrop( filterCakePressureDrop )
|
|
, m_injectivityFactor( injectivityFactor )
|
|
, m_filterCakeMobility( filterCakeMobility )
|
|
, m_filtrateThickness( filtrateThickness )
|
|
, m_viscosity( viscosity )
|
|
, m_relativePermeability( relativePermeability )
|
|
{
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
double RigEclipseToThermalCellTransmissibilityCalculator::calculateTransmissibility( const cvf::Vec3d& transmissibilityVector,
|
|
double fractureArea )
|
|
{
|
|
double fractureMatrixTransimissibility = transmissibilityVector.length();
|
|
|
|
if ( m_filterCakePressureDrop == RimThermalFractureTemplate::FilterCakePressureDrop::RELATIVE )
|
|
{
|
|
return m_injectivityFactor * fractureMatrixTransimissibility;
|
|
}
|
|
else if ( m_filterCakePressureDrop == RimThermalFractureTemplate::FilterCakePressureDrop::ABSOLUTE )
|
|
{
|
|
// Take filter cake transmissibility into account only when filter cake exists (i.e. has thickness).
|
|
if ( m_filtrateThickness > 0.0 )
|
|
{
|
|
double filterCakeTransmissibility = ( m_viscosity / m_relativePermeability ) * fractureArea * m_filterCakeMobility;
|
|
|
|
// Harmonic mean
|
|
return ( fractureMatrixTransimissibility * filterCakeTransmissibility ) /
|
|
( fractureMatrixTransimissibility + filterCakeTransmissibility );
|
|
}
|
|
else
|
|
{
|
|
return fractureMatrixTransimissibility;
|
|
}
|
|
}
|
|
else
|
|
{
|
|
CAF_ASSERT( m_filterCakePressureDrop == RimThermalFractureTemplate::FilterCakePressureDrop::NONE );
|
|
return fractureMatrixTransimissibility;
|
|
}
|
|
}
|