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

hopefully better initialization of wells

Browse Source
This commit is contained in:
Halvor M Nilsen 2025-02-14 14:37:57 +01:00
parent 3027dc6306
commit 3e3b72d037
3 changed files with 107 additions and 88 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@ -437,9 +437,9 @@ namespace Opm {
// nonzero phase anyway.
for (const auto& well : well_container_) {
//const bool zero_target = well->stoppedOrZeroRateTarget(simulator_, this->wellState(), local_deferredLogger);
if (well->isProducer()){// && !zero_target) {
//if (well->isProducer()){// && !zero_target) {
well->initializeWellState(simulator_, this->groupState(), this->wellState(), local_deferredLogger);
}
//}
}
}

73
opm/simulators/wells/MultisegmentWellPrimaryVariables.cpp
View File

@ -112,50 +112,59 @@ update(const WellState<Scalar>& well_state,
value_[seg][WQTotal] = 0;
}
assert(ws.initializedFromReservoir());
// tot to map old fraction to new perforations for now start from scratch.
bool prim_set = ws.stw_primaryvar.size() == value_.size();
if(prim_set){
//if (std::abs(total_seg_rate) > 0.) {
if (has_wfrac_variable) {
//const int water_pos = pu.phase_pos[Water];
//value_[seg][WFrac] = well_.scalingFactor(water_pos) * segment_rates[well_.numPhases() * seg + water_pos] / total_seg_rate;
value_[seg][WFrac] = ws.multiseg_primaryvar[seg][WFrac];
}
if (has_gfrac_variable) {
//const int gas_pos = pu.phase_pos[Gas];
//value_[seg][GFrac] = well_.scalingFactor(gas_pos) * segment_rates[well_.numPhases() * seg + gas_pos] / total_seg_rate;
value_[seg][GFrac] = ws.multiseg_primaryvar[seg][GFrac];
}
} else {
if (std::abs(total_seg_rate) > 0.) {
if (has_wfrac_variable) {
const int water_pos = pu.phase_pos[Water];
value_[seg][WFrac] = well_.scalingFactor(water_pos) * segment_rates[well_.numPhases() * seg + water_pos] / total_seg_rate;
}
if (has_gfrac_variable) {
const int gas_pos = pu.phase_pos[Gas];
value_[seg][GFrac] = well_.scalingFactor(gas_pos) * segment_rates[well_.numPhases() * seg + gas_pos] / total_seg_rate;
}
// what about water and gas injection?
// } else { // total_seg_rate == 0
// if (well_.isInjector()) {
// // only single phase injection handled
// auto phase = well.getInjectionProperties().injectorType;
} else { // total_seg_rate == 0
if (well_.isInjector()) {
// only single phase injection handled
auto phase = well.getInjectionProperties().injectorType;
// if (has_wfrac_variable) {
// if (phase == InjectorType::WATER) {
// value_[seg][WFrac] = 1.0;
// } else {
// value_[seg][WFrac] = 0.0;
// }
// }
if (has_wfrac_variable) {
if (phase == InjectorType::WATER) {
value_[seg][WFrac] = 1.0;
} else {
value_[seg][WFrac] = 0.0;
}
}
// if (has_gfrac_variable) {
// if (phase == InjectorType::GAS) {
// value_[seg][GFrac] = 1.0;
// } else {
// value_[seg][GFrac] = 0.0;
// }
// }
if (has_gfrac_variable) {
if (phase == InjectorType::GAS) {
value_[seg][GFrac] = 1.0;
} else {
value_[seg][GFrac] = 0.0;
}
}
// } else if (well_.isProducer()) { // producers
// if (has_wfrac_variable) {
// value_[seg][WFrac] = 1.0 / well_.numPhases();
// }
} else if (well_.isProducer()) { // producers
if (has_wfrac_variable) {
value_[seg][WFrac] = 1.0 / well_.numPhases();
}
// if (has_gfrac_variable) {
// value_[seg][GFrac] = 1.0 / well_.numPhases();
// }
// }
// }
if (has_gfrac_variable) {
value_[seg][GFrac] = 1.0 / well_.numPhases();
}
}
}
}
}
}

118
opm/simulators/wells/StandardWellPrimaryVariables.cpp
View File

@ -162,70 +162,80 @@ update(const WellState<Scalar>& well_state,
value_[WQTotal] = 0.;
}
}
assert(ws.initializedFromReservoir());
//if (std::abs(total_well_rate) > 0.) {
assert(ws.initializedFromReservoir());
bool prim_set = ws.stw_primaryvar.size()>0;
if(prim_set){
if constexpr (has_wfrac_variable) {
//value_[WFrac] = well_.scalingFactor(pu.phase_pos[Water]) * ws.surface_rates[pu.phase_pos[Water]] / total_well_rate;
value_[WFrac] = ws.stw_primaryvar[WFrac];
value_[WFrac] = ws.stw_primaryvar[WFrac];
}
if constexpr (has_gfrac_variable) {
// value_[GFrac] = well_.scalingFactor(pu.phase_pos[Gas]) *
// (ws.surface_rates[pu.phase_pos[Gas]] -
// (Indices::enableSolvent ? ws.sum_solvent_rates() : 0.0) ) / total_well_rate ;
value_[GFrac] = ws.stw_primaryvar[GFrac];
value_[GFrac] = ws.stw_primaryvar[GFrac];
}
if constexpr (Indices::enableSolvent) {
value_[SFrac] = ws.stw_primaryvar[SFrac];
//value_[SFrac] = well_.scalingFactor(Indices::contiSolventEqIdx) * ws.sum_solvent_rates() / total_well_rate ;
}
}else{
assert(prim_set==false);
if (std::abs(total_well_rate) > 0.) {
if constexpr (has_wfrac_variable) {
value_[WFrac] = well_.scalingFactor(pu.phase_pos[Water]) * ws.surface_rates[pu.phase_pos[Water]] / total_well_rate;
}
if constexpr (has_gfrac_variable) {
value_[GFrac] = well_.scalingFactor(pu.phase_pos[Gas]) *
(ws.surface_rates[pu.phase_pos[Gas]] -
(Indices::enableSolvent ? ws.sum_solvent_rates() : 0.0) ) / total_well_rate ;
}
if constexpr (Indices::enableSolvent) {
value_[SFrac] = well_.scalingFactor(Indices::contiSolventEqIdx) * ws.sum_solvent_rates() / total_well_rate ;
}
} else { // total_well_rate == 0
if (well_.isInjector()) {
// only single phase injection handled
if constexpr (has_wfrac_variable) {
if (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) {
auto phase = well_.wellEcl().getInjectionProperties().injectorType;
if (phase == InjectorType::WATER) {
value_[WFrac] = 1.0;
} else {
value_[WFrac] = 0.0;
}
}
}
if constexpr (has_gfrac_variable) {
if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
auto phase = well_.wellEcl().getInjectionProperties().injectorType;
if (phase == InjectorType::GAS) {
value_[GFrac] = (1.0 - well_.rsRvInj());
if constexpr (Indices::enableSolvent) {
value_[GFrac] = 1.0 - well_.rsRvInj() - well_.wsolvent();
value_[SFrac] = well_.wsolvent();
}
} else {
value_[GFrac] = 0.0;
}
}
}
// } else { // total_well_rate == 0
// if (well_.isInjector()) {
// // only single phase injection handled
// if constexpr (has_wfrac_variable) {
// if (FluidSystem::phaseIsActive(FluidSystem::waterPhaseIdx)) {
// auto phase = well_.wellEcl().getInjectionProperties().injectorType;
// if (phase == InjectorType::WATER) {
// value_[WFrac] = 1.0;
// } else {
// value_[WFrac] = 0.0;
// }
// }
// }
// if constexpr (has_gfrac_variable) {
// if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
// auto phase = well_.wellEcl().getInjectionProperties().injectorType;
// if (phase == InjectorType::GAS) {
// value_[GFrac] = (1.0 - well_.rsRvInj());
// if constexpr (Indices::enableSolvent) {
// value_[GFrac] = 1.0 - well_.rsRvInj() - well_.wsolvent();
// value_[SFrac] = well_.wsolvent();
// }
// } else {
// value_[GFrac] = 0.0;
// }
// }
// }
// TODO: it is possible to leave injector as a oil well,
// when F_w and F_g both equals to zero, not sure under what kind of circumstance
// this will happen.
} else if (well_.isProducer()) { // producers
// TODO: the following are not addressed for the solvent case yet
if constexpr (has_wfrac_variable) {
value_[WFrac] = 1.0 / np;
}
// // TODO: it is possible to leave injector as a oil well,
// // when F_w and F_g both equals to zero, not sure under what kind of circumstance
// // this will happen.
// } else if (well_.isProducer()) { // producers
// // TODO: the following are not addressed for the solvent case yet
// if constexpr (has_wfrac_variable) {
// value_[WFrac] = 1.0 / np;
// }
// if constexpr (has_gfrac_variable) {
// value_[GFrac] = 1.0 / np;
// }
// } else {
// OPM_DEFLOG_THROW(std::logic_error, "Expected PRODUCER or INJECTOR type of well", deferred_logger);
// }
// }
// BHP
if constexpr (has_gfrac_variable) {
value_[GFrac] = 1.0 / np;
}
} else {
OPM_DEFLOG_THROW(std::logic_error, "Expected PRODUCER or INJECTOR type of well", deferred_logger);
}
}
}
//BHP
value_[Bhp] = ws.bhp;
}