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

Browse Source
This commit is contained in:
Kjetil Olsen Lye 2012-05-03 15:35:59 +02:00
commit 66b4122be9
4 changed files with 23 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/core/WellsManager.cpp
View File

@ -492,7 +492,7 @@ namespace Opm
control_pos[ProductionControl::ORAT] = w_->ctrls[wix]->num; control_pos[ProductionControl::ORAT] = w_->ctrls[wix]->num;
double distr[3] = { 0.0, 0.0, 0.0 }; double distr[3] = { 0.0, 0.0, 0.0 };
distr[pu.phase_pos[BlackoilPhases::Liquid]] = 1.0; distr[pu.phase_pos[BlackoilPhases::Liquid]] = 1.0;
ok = append_well_controls(SURFACE_RATE, wcp_line.oil_max_rate_, ok = append_well_controls(SURFACE_RATE, -wcp_line.oil_max_rate_,
distr, wix, w_); distr, wix, w_);
} }
if (ok && wcp_line.water_max_rate_ > 0.0) { if (ok && wcp_line.water_max_rate_ > 0.0) {
@ -502,7 +502,7 @@ namespace Opm
control_pos[ProductionControl::WRAT] = w_->ctrls[wix]->num; control_pos[ProductionControl::WRAT] = w_->ctrls[wix]->num;
double distr[3] = { 0.0, 0.0, 0.0 }; double distr[3] = { 0.0, 0.0, 0.0 };
distr[pu.phase_pos[BlackoilPhases::Aqua]] = 1.0; distr[pu.phase_pos[BlackoilPhases::Aqua]] = 1.0;
ok = append_well_controls(SURFACE_RATE, wcp_line.water_max_rate_, ok = append_well_controls(SURFACE_RATE, -wcp_line.water_max_rate_,
distr, wix, w_); distr, wix, w_);
} }
if (ok && wcp_line.gas_max_rate_ > 0.0) { if (ok && wcp_line.gas_max_rate_ > 0.0) {
@ -512,7 +512,7 @@ namespace Opm
control_pos[ProductionControl::GRAT] = w_->ctrls[wix]->num; control_pos[ProductionControl::GRAT] = w_->ctrls[wix]->num;
double distr[3] = { 0.0, 0.0, 0.0 }; double distr[3] = { 0.0, 0.0, 0.0 };
distr[pu.phase_pos[BlackoilPhases::Vapour]] = 1.0; distr[pu.phase_pos[BlackoilPhases::Vapour]] = 1.0;
ok = append_well_controls(SURFACE_RATE, wcp_line.gas_max_rate_, ok = append_well_controls(SURFACE_RATE, -wcp_line.gas_max_rate_,
distr, wix, w_); distr, wix, w_);
} }
if (ok && wcp_line.liquid_max_rate_ > 0.0) { if (ok && wcp_line.liquid_max_rate_ > 0.0) {
@ -526,13 +526,13 @@ namespace Opm
double distr[3] = { 0.0, 0.0, 0.0 }; double distr[3] = { 0.0, 0.0, 0.0 };
distr[pu.phase_pos[BlackoilPhases::Aqua]] = 1.0; distr[pu.phase_pos[BlackoilPhases::Aqua]] = 1.0;
distr[pu.phase_pos[BlackoilPhases::Liquid]] = 1.0; distr[pu.phase_pos[BlackoilPhases::Liquid]] = 1.0;
ok = append_well_controls(SURFACE_RATE, wcp_line.liquid_max_rate_, ok = append_well_controls(SURFACE_RATE, -wcp_line.liquid_max_rate_,
distr, wix, w_); distr, wix, w_);
} }
if (ok && wcp_line.reservoir_flow_max_rate_ > 0.0) { if (ok && wcp_line.reservoir_flow_max_rate_ > 0.0) {
control_pos[ProductionControl::RESV] = w_->ctrls[wix]->num; control_pos[ProductionControl::RESV] = w_->ctrls[wix]->num;
double distr[3] = { 1.0, 1.0, 1.0 }; double distr[3] = { 1.0, 1.0, 1.0 };
ok = append_well_controls(RESERVOIR_RATE, wcp_line.reservoir_flow_max_rate_, ok = append_well_controls(RESERVOIR_RATE, -wcp_line.reservoir_flow_max_rate_,
distr, wix, w_); distr, wix, w_);
} }
if (ok && wcp_line.BHP_limit_ > 0.0) { if (ok && wcp_line.BHP_limit_ > 0.0) {

1
opm/core/pressure/tpfa/ifs_tpfa.c
View File

@ -287,7 +287,6 @@ assemble_well_contrib(int nc ,
/* ---------------------------------------------------------------------- */ /* ---------------------------------------------------------------------- */
{ {
int w, p, np; int w, p, np;
int are_rate;
struct WellControls *ctrls; struct WellControls *ctrls;

12
opm/core/utility/initState_impl.hpp
View File

@ -214,16 +214,16 @@ namespace Opm
// Initialise pressure to hydrostatic state. // Initialise pressure to hydrostatic state.
const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa; const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa;
initHydrostaticPressure(grid, props.density(), woc, gravity, woc, ref_p, state); initHydrostaticPressure(grid, props.density(), woc, gravity, woc, ref_p, state);
} else if (param.has("init_sat")) { } else if (param.has("init_saturation")) {
// Initialise water saturation to init_sat parameter. // Initialise water saturation to init_saturation parameter.
const double init_sat = param.get<double>("init_sat"); const double init_saturation = param.get<double>("init_saturation");
for (int cell = 0; cell < num_cells; ++cell) { for (int cell = 0; cell < num_cells; ++cell) {
state.saturation()[2*cell] = init_sat; state.saturation()[2*cell] = init_saturation;
state.saturation()[2*cell + 1] = 1.0 - init_sat; state.saturation()[2*cell + 1] = 1.0 - init_saturation;
} }
// Initialise pressure to hydrostatic state. // Initialise pressure to hydrostatic state.
const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa; const double ref_p = param.getDefault("ref_pressure", 100)*unit::barsa;
const double rho = props.density()[0]*init_sat + props.density()[1]*(1.0 - init_sat); const double rho = props.density()[0]*init_saturation + props.density()[1]*(1.0 - init_saturation);
const double dens[2] = { rho, rho }; const double dens[2] = { rho, rho };
const double ref_z = grid.cell_centroids[0 + grid.dimensions - 1]; const double ref_z = grid.cell_centroids[0 + grid.dimensions - 1];
initHydrostaticPressure(grid, dens, ref_z, gravity, ref_z, ref_p, state); initHydrostaticPressure(grid, dens, ref_z, gravity, ref_z, ref_p, state);

15
opm/core/utility/miscUtilities.cpp
View File

@ -410,16 +410,25 @@ namespace Opm
} }
src.resize(num_cells); src.resize(num_cells);
for (int w = 0; w < wells.number_of_wells; ++w) { for (int w = 0; w < wells.number_of_wells; ++w) {
const int cur = wells.ctrls[w]->current;
if (wells.ctrls[w]->num != 1) { if (wells.ctrls[w]->num != 1) {
THROW("In wellsToSrc(): well has more than one control."); MESSAGE("In wellsToSrc(): well has more than one control, all but current control will be ignored.");
} }
if (wells.ctrls[w]->type[0] != RESERVOIR_RATE) { if (wells.ctrls[w]->type[cur] != RESERVOIR_RATE) {
THROW("In wellsToSrc(): well is something other than RESERVOIR_RATE."); THROW("In wellsToSrc(): well is something other than RESERVOIR_RATE.");
} }
if (wells.well_connpos[w+1] - wells.well_connpos[w] != 1) { if (wells.well_connpos[w+1] - wells.well_connpos[w] != 1) {
THROW("In wellsToSrc(): well has multiple perforations."); THROW("In wellsToSrc(): well has multiple perforations.");
} }
const double flow = wells.ctrls[w]->target[0] * wells.ctrls[w]->distr[0]; for (int p = 0; p < np; ++p) {
if (wells.ctrls[w]->distr[np*cur + p] != 1.0) {
THROW("In wellsToSrc(): well not controlled on total rate.");
}
}
double flow = wells.ctrls[w]->target[cur];
if (wells.type[w] == INJECTOR) {
flow *= wells.comp_frac[np*w + 0]; // Obtaining water rate for inflow source.
}
const double cell = wells.well_cells[wells.well_connpos[w]]; const double cell = wells.well_cells[wells.well_connpos[w]];
src[cell] = flow; src[cell] = flow;
} }