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

Zero initialize wells with no perforations.

Browse Source 
For these wells access its well_cells might read of the bounds
an array if they are the last wells in the struct. Therefore
we cannnot initialiue first_cell and the well control is uninitialized,
to.

With this commit theses wells are now detected and theor bhp, thp, and well_rates
are initialized to zero.
This commit is contained in:
Markus Blatt 2016-09-21 12:23:34 +02:00
parent c89b84eeed
commit 2d115888de
1 changed files with 13 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
opm/core/simulator/WellState.hpp
View File

@ -65,6 +65,7 @@ namespace Opm
for (int w = 0; w < nw; ++w) { for (int w = 0; w < nw; ++w) {
assert((wells->type[w] == INJECTOR) || (wells->type[w] == PRODUCER)); assert((wells->type[w] == INJECTOR) || (wells->type[w] == PRODUCER));
const WellControls* ctrl = wells->ctrls[w]; const WellControls* ctrl = wells->ctrls[w];
const int num_perf_this_well = wells->well_connpos[w + 1] - wells->well_connpos[w];
// setup wellname -> well index mapping // setup wellname -> well index mapping
{ {
@ -75,10 +76,20 @@ namespace Opm
wellMapEntry[ 0 ] = w; wellMapEntry[ 0 ] = w;
wellMapEntry[ 1 ] = wells->well_connpos[w]; wellMapEntry[ 1 ] = wells->well_connpos[w];
// also store the number of perforations in this well // also store the number of perforations in this well
const int num_perf_this_well = wells->well_connpos[w + 1] - wells->well_connpos[w];
wellMapEntry[ 2 ] = num_perf_this_well; wellMapEntry[ 2 ] = num_perf_this_well;
} }
if ( num_perf_this_well == 0 )
{
// No perforations of the well. Initialize to zero.
for (int p = 0; p < np; ++p) {
wellrates_[np*w + p] = 0.0;
}
bhp_[w] = 0;
thp_[w] = 0;
continue;
}
if (well_controls_well_is_stopped(ctrl)) { if (well_controls_well_is_stopped(ctrl)) {
// Stopped well: // Stopped well:
// 1. Rates: assign zero well rates. // 1. Rates: assign zero well rates.
@ -143,7 +154,7 @@ namespace Opm
} else { } else {
thp_[w] = bhp_[w]; thp_[w] = bhp_[w];
} }
} }
} }
// The perforation rates and perforation pressures are // The perforation rates and perforation pressures are