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

Implements initialization for constant capillary pressure functions

Browse Source
This commit is contained in:
Tor Harald Sandve 2014-08-11 11:23:15 +02:00 committed by Andreas Lauser
parent c61fb736fa
commit 70dc519539
2 changed files with 33 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/core/simulator/initStateEquil.hpp
View File

@ -403,7 +403,7 @@ namespace Opm
PVec press = phasePressures(G, eqreg, cells, grav); PVec press = phasePressures(G, eqreg, cells, grav);
const PVec sat = phaseSaturations(eqreg, cells, props, swat_init_, press); const PVec sat = phaseSaturations(G, eqreg, cells, props, swat_init_, press);
const int np = props.numPhases(); const int np = props.numPhases();
for (int p = 0; p < np; ++p) { for (int p = 0; p < np; ++p) {

21
opm/core/simulator/initStateEquil_impl.hpp
View File

@ -578,7 +578,8 @@ namespace Opm
template <class Region, class CellRange> template <class Region, class CellRange>
std::vector< std::vector<double> > std::vector< std::vector<double> >
phaseSaturations(const Region& reg, phaseSaturations(const UnstructuredGrid& G,
const Region& reg,
const CellRange& cells, const CellRange& cells,
BlackoilPropertiesInterface& props, BlackoilPropertiesInterface& props,
const std::vector<double> swat_init, const std::vector<double> swat_init,
@ -598,6 +599,8 @@ namespace Opm
double smin[BlackoilPhases::MaxNumPhases] = { 0.0 }; double smin[BlackoilPhases::MaxNumPhases] = { 0.0 };
double smax[BlackoilPhases::MaxNumPhases] = { 0.0 }; double smax[BlackoilPhases::MaxNumPhases] = { 0.0 };
const double* const cc = & G.cell_centroids[0];
const bool water = reg.phaseUsage().phase_used[BlackoilPhases::Aqua]; const bool water = reg.phaseUsage().phase_used[BlackoilPhases::Aqua];
const bool gas = reg.phaseUsage().phase_used[BlackoilPhases::Vapour]; const bool gas = reg.phaseUsage().phase_used[BlackoilPhases::Vapour];
const int oilpos = reg.phaseUsage().phase_pos[BlackoilPhases::Liquid]; const int oilpos = reg.phaseUsage().phase_pos[BlackoilPhases::Liquid];
@ -611,6 +614,13 @@ namespace Opm
// inverting capillary pressure functions. // inverting capillary pressure functions.
double sw = 0.0; double sw = 0.0;
if (water) { if (water) {
if (isConstPc(props,waterpos,cell)){
const int nd = G.dimensions;
const double cellDepth = cc[nd * cell + nd-1];
sw = satFromDepth(props,cellDepth,zwoc,waterpos,cell,false);
phase_saturations[waterpos][local_index] = sw;
}
else{
const double pcov = phase_pressures[oilpos][local_index] - phase_pressures[waterpos][local_index]; const double pcov = phase_pressures[oilpos][local_index] - phase_pressures[waterpos][local_index];
if (swat_init.empty()) { // Invert Pc to find sw if (swat_init.empty()) { // Invert Pc to find sw
sw = satFromPc(props, waterpos, cell, pcov); sw = satFromPc(props, waterpos, cell, pcov);
@ -621,14 +631,23 @@ namespace Opm
phase_saturations[waterpos][local_index] = sw; phase_saturations[waterpos][local_index] = sw;
} }
} }
}
double sg = 0.0; double sg = 0.0;
if (gas) { if (gas) {
if (isConstPc(props,gaspos,cell)){
const int nd = G.dimensions;
const double cellDepth = cc[nd * cell + nd-1];
sg = satFromDepth(props,cellDepth,zgoc,gaspos,cell,true);
phase_saturations[gaspos][local_index] = sg;
}
else{
// Note that pcog is defined to be (pg - po), not (po - pg). // Note that pcog is defined to be (pg - po), not (po - pg).
const double pcog = phase_pressures[gaspos][local_index] - phase_pressures[oilpos][local_index]; const double pcog = phase_pressures[gaspos][local_index] - phase_pressures[oilpos][local_index];
const double increasing = true; // pcog(sg) expected to be increasing function const double increasing = true; // pcog(sg) expected to be increasing function
sg = satFromPc(props, gaspos, cell, pcog, increasing); sg = satFromPc(props, gaspos, cell, pcog, increasing);
phase_saturations[gaspos][local_index] = sg; phase_saturations[gaspos][local_index] = sg;
} }
}
if (gas && water && (sg + sw > 1.0)) { if (gas && water && (sg + sw > 1.0)) {
// Overlapping gas-oil and oil-water transition // Overlapping gas-oil and oil-water transition
// zones can lead to unphysical saturations when // zones can lead to unphysical saturations when