OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-28 03:53:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

The water/gas pressures are only calculated if woc and goc lies within

Browse Source 
the reservoar
The water/gas pressure is set to -inf when woc and goc is above or below
the reservoar.
This commit is contained in:
Tor Harald Sandve 2014-08-15 08:56:35 +02:00 committed by Andreas Lauser
parent 480f5337b6
commit b1d20aa83d
1 changed files with 49 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -394,22 +394,30 @@ namespace Opm
assign(G, opress, z0, cells, press); assign(G, opress, z0, cells, press);
// Default pressure at WOC and GOC to -inf
po_woc = -std::numeric_limits<double>::max(); po_woc = -std::numeric_limits<double>::max();
po_goc = -std::numeric_limits<double>::max(); po_goc = -std::numeric_limits<double>::max();
const double woc = reg.zwoc(); const double woc = reg.zwoc();
// Compute Oil pressure at WOC // Compute Oil pressure at WOC,
// if WOC is within the reservoir.
if ( ( woc > span[0] ) & ( woc < span[1] ) ){
if (z0 > woc) { po_woc = opress[0](woc); } // WOC above datum if (z0 > woc) { po_woc = opress[0](woc); } // WOC above datum
else if (z0 < woc) { po_woc = opress[1](woc); } // WOC below datum else if (z0 < woc) { po_woc = opress[1](woc); } // WOC below datum
else { po_woc = p0; } // WOC *at* datum else { po_woc = p0; } // WOC *at* datum
}
const double goc = reg.zgoc(); const double goc = reg.zgoc();
// Compute Oil pressure at GOC // Compute Oil pressure at GOC,
// if GOC is within the reservoir.
if ( ( goc > span[0] ) & ( goc < span[1] ) ){
if (z0 > goc) { po_goc = opress[0](goc); } // GOC above datum if (z0 > goc) { po_goc = opress[0](goc); } // GOC above datum
else if (z0 < goc) { po_goc = opress[1](goc); } // GOC below datum else if (z0 < goc) { po_goc = opress[1](goc); } // GOC below datum
else { po_goc = p0; } // GOC *at* datum else { po_goc = p0; } // GOC *at* datum
} }
}
template <class Region, template <class Region,
class CellRange> class CellRange>
void void
@ -474,14 +482,43 @@ namespace Opm
if (PhaseUsed::water(pu)) { if (PhaseUsed::water(pu)) {
const int wix = PhaseIndex::water(pu); const int wix = PhaseIndex::water(pu);
// If woc is above or below the reservoar,
// po_woc is -inf and the water pressure
// is set to -inf.
if ( po_woc > 0 ){
PhasePressure::water(G, reg, span, grav, po_woc, PhasePressure::water(G, reg, span, grav, po_woc,
cells, press[ wix ]); cells, press[ wix ]);
} else {
std::vector<double>::size_type local_index = 0;
for (typename CellRange::const_iterator ci = cells.begin(); ci != cells.end(); ++ci, ++local_index) {
const int cell = *ci;
press[wix][cell] = po_woc;
}
}
} }
if (PhaseUsed::gas(pu)) { if (PhaseUsed::gas(pu)) {
const int gix = PhaseIndex::gas(pu); const int gix = PhaseIndex::gas(pu);
// If woc is above or below the reservoar,
// po_woc is -inf and the water pressure
// is set to -inf.
if (po_goc > 0){
PhasePressure::gas(G, reg, span, grav, po_goc, PhasePressure::gas(G, reg, span, grav, po_goc,
cells, press[ gix ]); cells, press[ gix ]);
} else {
std::vector<double>::size_type local_index = 0;
for (typename CellRange::const_iterator ci = cells.begin(); ci != cells.end(); ++ci, ++local_index) {
const int cell = *ci;
press[gix][cell] = po_goc;
}
}
} }
} }
} // namespace Details } // namespace Details