OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #4559 from totto82/rockfraction

Browse Source 
compute rock fraction
This commit is contained in:
Bård Skaflestad 2023-03-29 11:44:01 +02:00 committed by GitHub
commit 5107cd242d
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
2 changed files with 25 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
ebos/eclgenericproblem.cc
View File

@ -347,6 +347,20 @@ porosity(unsigned globalSpaceIdx, unsigned timeIdx) const
return this->referencePorosity_[timeIdx][globalSpaceIdx];
}
template<class GridView, class FluidSystem, class Scalar>
Scalar EclGenericProblem<GridView,FluidSystem,Scalar>::
rockFraction(unsigned elementIdx, unsigned timeIdx) const
{
const auto& fp = eclState_.fieldProps();
const std::vector<double>& poro = fp.get_double("PORO");
// the reference porosity is defined as the accumulated pore volume divided by the
// geometric volume of the element. Note that it can
// be larger than 1.0 if porevolume multipliers are used
// to for instance implement larger boundary cells
Scalar porosity = poro[elementIdx];
return referencePorosity(elementIdx, timeIdx) / porosity * (1 - porosity);
}
template<class GridView, class FluidSystem, class Scalar>
template<class T>
void EclGenericProblem<GridView,FluidSystem,Scalar>::

11
ebos/eclgenericproblem.hh
View File

@ -138,6 +138,17 @@ public:
Scalar referencePorosity(unsigned elementIdx, unsigned timeIdx) const
{ return referencePorosity_[timeIdx][elementIdx]; }
/*!
* \brief Returns the rockFraction of an element
*
* Usually 1 - porosity, but if pvmult is used to modify porosity
* we will apply the same multiplier to the rock fraction
* i.e. pvmult*(1 - porosity) and thus interpret multpv as a volume
* multiplier. This is to avoid negative rock volume for pvmult*porosity > 1
*/
Scalar rockFraction(unsigned elementIdx, unsigned timeIdx) const;
/*!
* \brief Sets the porosity of an element
*