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

making the VREP enforcement into a small method.

Browse Source 
in BlackoilModelBase.
This commit is contained in:
Kai Bao 2016-11-02 14:26:57 +01:00
parent 9a95697fde
commit 9150c217f2
2 changed files with 27 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/autodiff/BlackoilModelBase.hpp
View File

@ -319,6 +319,10 @@ namespace Opm {
std::vector<double>& well_voidage_rates,
std::vector<double>& voidage_conversion_coeffs);
void applyVREPGroupControl(const ReservoirState& reservoir_state,
const WellState& well_state);
protected:
// --------- Types and enums ---------

40
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -1063,11 +1063,8 @@ namespace detail {
converged = getWellConvergence(it);
// Enforce the VREP control
if (it == 0 && asImpl().wellModel().wellCollection()->havingVREPGroups()) {
std::vector<double> well_voidage_rates;
std::vector<double> voidage_conversion_coeffs;
computeWellVoidageRates(reservoir_state, well_state, well_voidage_rates, voidage_conversion_coeffs);
asImpl().wellModel().wellCollection()->applyVREPGroupControls(well_voidage_rates, voidage_conversion_coeffs);
if (it == 0) {
applyVREPGroupControl(reservoir_state, well_state);
}
// When the well targets are just updated or need to be updated, we need at least one more iteration.
@ -1103,12 +1100,7 @@ namespace detail {
asImpl().wellModel().updateWellState(dx.array(), dpMaxRel(), well_state);
// Enforce the VREP control
if (asImpl().wellModel().wellCollection()->havingVREPGroups()) {
std::vector<double> well_voidage_rates;
std::vector<double> voidage_conversion_coeffs;
computeWellVoidageRates(reservoir_state, well_state, well_voidage_rates, voidage_conversion_coeffs);
asImpl().wellModel().wellCollection()->applyVREPGroupControls(well_voidage_rates, voidage_conversion_coeffs);
}
applyVREPGroupControl(reservoir_state, well_state);
}
// We have to update the well controls regardless whether there are local
// wells active or not as parallel logging will take place that needs to
@ -1529,12 +1521,8 @@ namespace detail {
asImpl().wellModel().updateWellState(dwells, dpMaxRel(), well_state);
if (asImpl().wellModel().wellCollection()->havingVREPGroups()) {
std::vector<double> well_voidage_rates;
std::vector<double> voidage_conversion_coeffs;
computeWellVoidageRates(reservoir_state, well_state, well_voidage_rates, voidage_conversion_coeffs);
asImpl().wellModel().wellCollection()->applyVREPGroupControls(well_voidage_rates, voidage_conversion_coeffs);
}
// enforce VREP control when necessary.
applyVREPGroupControl(reservoir_state, well_state);
// Update phase conditions used for property calculations.
updatePhaseCondFromPrimalVariable(reservoir_state);
@ -2677,6 +2665,24 @@ namespace detail {
}
}
template <class Grid, class WellModel, class Implementation>
void
BlackoilModelBase<Grid, WellModel, Implementation>::
applyVREPGroupControl(const ReservoirState& reservoir_state,
const WellState& well_state)
{
if (asImpl().wellModel().wellCollection()->havingVREPGroups()) {
std::vector<double> well_voidage_rates;
std::vector<double> voidage_conversion_coeffs;
computeWellVoidageRates(reservoir_state, well_state, well_voidage_rates, voidage_conversion_coeffs);
asImpl().wellModel().wellCollection()->applyVREPGroupControls(well_voidage_rates, voidage_conversion_coeffs);
}
}
} // namespace Opm
#endif // OPM_BLACKOILMODELBASE_IMPL_HEADER_INCLUDED