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 pull request #3879 from totto82/rsrvinj2

Browse Source 
WIP implement rsrv item 10 WCONINJE
This commit is contained in:
Tor Harald Sandve
2022-06-16 08:30:31 +02:00
committed by GitHub
parent 7512c0c247 1f876d1800
commit 60d0fe0f87
5 changed files with 32 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@@ -76,6 +76,11 @@ namespace Opm
if constexpr (Base::has_watVapor) { if constexpr (Base::has_watVapor) {
OPM_THROW(std::runtime_error, "water evaporation is not supported by multisegment well yet"); OPM_THROW(std::runtime_error, "water evaporation is not supported by multisegment well yet");
} }
if(this->rsRvInj() > 0) {
OPM_THROW(std::runtime_error, "dissolved gas/ vapporized oil in injected oil/gas not supported by multisegment well yet."
<< " \n See (WCONINJE item 10 / WCONHIST item 8)");
}
} }

15
opm/simulators/wells/StandardWellEval.cpp
View File

@@ -195,12 +195,19 @@ getQs(const int comp_idx) const
if (Indices::enableSolvent && comp_idx == Indices::contiSolventEqIdx) { // solvent if (Indices::enableSolvent && comp_idx == Indices::contiSolventEqIdx) { // solvent
inj_frac = baseif_.wsolvent(); inj_frac = baseif_.wsolvent();
} else if (comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::gasCompIdx))) { } else if (comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::gasCompIdx))) {
inj_frac = Indices::enableSolvent ? 1.0 - baseif_.wsolvent() : 1.0; inj_frac = 1.0 - baseif_.rsRvInj();
if (Indices::enableSolvent) {
inj_frac -= baseif_.wsolvent();
}
} else if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) && comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::oilCompIdx))) {
inj_frac = baseif_.rsRvInj();
} }
break; break;
case InjectorType::OIL: case InjectorType::OIL:
if (comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::oilCompIdx))) { if (comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::oilCompIdx))) {
inj_frac = 1.0; inj_frac = 1.0 - baseif_.rsRvInj();
} else if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx) && comp_idx == int(Indices::canonicalToActiveComponentIndex(FluidSystem::gasCompIdx))) {
inj_frac = baseif_.rsRvInj();
} }
break; break;
case InjectorType::MULTI: case InjectorType::MULTI:
@@ -316,9 +323,9 @@ updatePrimaryVariables(const WellState& well_state, DeferredLogger& deferred_log
if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) { if (FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
auto phase = baseif_.wellEcl().getInjectionProperties().injectorType; auto phase = baseif_.wellEcl().getInjectionProperties().injectorType;
if (phase == InjectorType::GAS) { if (phase == InjectorType::GAS) {
primary_variables_[GFrac] = 1.0; primary_variables_[GFrac] = (1.0 - baseif_.rsRvInj());
if constexpr (Indices::enableSolvent) { if constexpr (Indices::enableSolvent) {
primary_variables_[GFrac] = 1.0 - baseif_.wsolvent(); primary_variables_[GFrac] = 1.0 - baseif_.rsRvInj() - baseif_.wsolvent();
primary_variables_[SFrac] = baseif_.wsolvent(); primary_variables_[SFrac] = baseif_.wsolvent();
} }
} else { } else {

5
opm/simulators/wells/WellInterfaceGeneric.cpp
View File

@@ -173,6 +173,11 @@ double WellInterfaceGeneric::wsolvent() const
return wsolvent_; return wsolvent_;
} }
double WellInterfaceGeneric::rsRvInj() const
{
return well_ecl_.getInjectionProperties().rsRvInj;
}
bool WellInterfaceGeneric::wellHasTHPConstraints(const SummaryState& summaryState) const bool WellInterfaceGeneric::wellHasTHPConstraints(const SummaryState& summaryState) const
{ {
if (dynamic_thp_limit_) { if (dynamic_thp_limit_) {

1
opm/simulators/wells/WellInterfaceGeneric.hpp
View File

@@ -170,6 +170,7 @@ public:
double getTHPConstraint(const SummaryState& summaryState) const; double getTHPConstraint(const SummaryState& summaryState) const;
double getALQ(const WellState& well_state) const; double getALQ(const WellState& well_state) const;
double wsolvent() const; double wsolvent() const;
double rsRvInj() const;
// whether a well is specified with a non-zero and valid VFP table number // whether a well is specified with a non-zero and valid VFP table number
bool isVFPActive(DeferredLogger& deferred_logger) const; bool isVFPActive(DeferredLogger& deferred_logger) const;

11
opm/simulators/wells/WellInterface_impl.hpp
View File

@@ -718,7 +718,16 @@ namespace Opm
switch(current) { switch(current) {
case Well::InjectorCMode::RATE: case Well::InjectorCMode::RATE:
{ {
ws.surface_rates[phasePos] = controls.surface_rate; ws.surface_rates[phasePos] = (1.0 - this->rsRvInj()) * controls.surface_rate;
if(this->rsRvInj() > 0) {
if (injectorType == InjectorType::OIL && FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
ws.surface_rates[pu.phase_pos[BlackoilPhases::Vapour]] = controls.surface_rate * this->rsRvInj();
} else if (injectorType == InjectorType::GAS && FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx)) {
ws.surface_rates[pu.phase_pos[BlackoilPhases::Liquid]] = controls.surface_rate * this->rsRvInj();
} else {
OPM_DEFLOG_THROW(std::runtime_error, "Expected OIL or GAS as type for injectors when RS/RV (item 10) is non-zero " + this->name(), deferred_logger );
}
}
break; break;
} }