OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-15 20:01:56 -06:00
Code Issues Packages Projects Releases Wiki Activity

some cleanup

Browse Source
This commit is contained in:
Tor Harald Sandve 2020-12-16 12:58:17 +01:00
parent 30ce605713
commit 273538e1f4
2 changed files with 18 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/simulators/wells/MultisegmentWell.hpp
View File

@ -340,7 +340,8 @@ namespace Opm
// scale the segment rates and pressure based on well rates and bhp // scale the segment rates and pressure based on well rates and bhp
void scaleSegmentRatesAndPressureWithWellRatesAndPressure(WellState& well_state) const; void scaleSegmentRatesWithWellRates(WellState& well_state) const;
void scaleSegmentPressuresWithBhp(WellState& well_state) const;
// computing the accumulation term for later use in well mass equations // computing the accumulation term for later use in well mass equations
void computeInitialSegmentFluids(const Simulator& ebos_simulator); void computeInitialSegmentFluids(const Simulator& ebos_simulator);

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

@ -296,7 +296,6 @@ namespace Opm
const auto& well = well_ecl_; const auto& well = well_ecl_;
const int well_index = index_of_well_; const int well_index = index_of_well_;
const int top_segment_index = well_state.topSegmentIndex(index_of_well_);
const auto& pu = phaseUsage(); const auto& pu = phaseUsage();
const int np = well_state.numPhases(); const int np = well_state.numPhases();
const auto& summaryState = ebos_simulator.vanguard().summaryState(); const auto& summaryState = ebos_simulator.vanguard().summaryState();
@ -512,7 +511,8 @@ namespace Opm
// scale segment rates based on the wellRates // scale segment rates based on the wellRates
// and segment pressure based on bhp // and segment pressure based on bhp
scaleSegmentRatesAndPressureWithWellRatesAndPressure(well_state); scaleSegmentPressuresWithBhp(well_state);
scaleSegmentRatesWithWellRates(well_state);
} }
@ -522,10 +522,9 @@ namespace Opm
template <typename TypeTag> template <typename TypeTag>
void void
MultisegmentWell<TypeTag>:: MultisegmentWell<TypeTag>::
scaleSegmentRatesAndPressureWithWellRatesAndPressure(WellState& well_state) const scaleSegmentRatesWithWellRates(WellState& well_state) const
{ {
const int top_segment_index = well_state.topSegmentIndex(index_of_well_); const int top_segment_index = well_state.topSegmentIndex(index_of_well_);
for (int phase = 0; phase < number_of_phases_; ++phase) { for (int phase = 0; phase < number_of_phases_; ++phase) {
const double well_phase_rate = well_state.wellRates()[number_of_phases_*index_of_well_ + phase]; const double well_phase_rate = well_state.wellRates()[number_of_phases_*index_of_well_ + phase];
const double unscaled_top_seg_rate = well_state.segRates()[number_of_phases_*top_segment_index + phase]; const double unscaled_top_seg_rate = well_state.segRates()[number_of_phases_*top_segment_index + phase];
@ -544,22 +543,27 @@ namespace Opm
} }
std::vector<double> perforation_rates(number_of_phases_ * number_of_perforations_,0.0); std::vector<double> perforation_rates(number_of_phases_ * number_of_perforations_,0.0);
for (int phase = 0; phase < number_of_phases_; ++phase) {
const double perf_phaserate_scaled = well_state.wellRates()[number_of_phases_ * index_of_well_ + phase] / sumTw; const double perf_phaserate_scaled = well_state.wellRates()[number_of_phases_ * index_of_well_ + phase] / sumTw;
for (int perf = 0; perf < number_of_perforations_; ++perf) { for (int perf = 0; perf < number_of_perforations_; ++perf) {
perforation_rates[number_of_phases_ * perf + phase] = well_index_[perf] * perf_phaserate_scaled; perforation_rates[number_of_phases_ * perf + phase] = well_index_[perf] * perf_phaserate_scaled;
} }
}
std::vector<double> segment_rates; std::vector<double> segment_rates;
WellState::calculateSegmentRates(segment_inlets_, segment_perforations_, perforation_rates, number_of_phases_, WellState::calculateSegmentRates(segment_inlets_, segment_perforations_, perforation_rates, number_of_phases_,
0, segment_rates); 0, segment_rates);
const int top_segment_index = well_state.topSegmentIndex(index_of_well_);
std::copy(segment_rates.begin(), segment_rates.end(), std::copy(segment_rates.begin(), segment_rates.end(),
well_state.segRates().begin() + number_of_phases_ * top_segment_index ); well_state.segRates().begin() + number_of_phases_ * top_segment_index );
} }
} }
}
template <typename TypeTag>
void
MultisegmentWell<TypeTag>::
scaleSegmentPressuresWithBhp(WellState& well_state) const
{
//scale segment pressures //scale segment pressures
const int top_segment_index = well_state.topSegmentIndex(index_of_well_);
const double bhp = well_state.bhp()[index_of_well_]; const double bhp = well_state.bhp()[index_of_well_];
const double unscaled_top_seg_pressure = well_state.segPress()[top_segment_index]; const double unscaled_top_seg_pressure = well_state.segPress()[top_segment_index];
for (int seg = 0; seg < numberOfSegments(); ++seg) { for (int seg = 0; seg < numberOfSegments(); ++seg) {
@ -569,7 +573,6 @@ namespace Opm
} }
template <typename TypeTag> template <typename TypeTag>
ConvergenceReport ConvergenceReport
MultisegmentWell<TypeTag>:: MultisegmentWell<TypeTag>::
@ -890,6 +893,7 @@ namespace Opm
well_state_copy.currentProductionControls()[index_of_well_] = Well::ProducerCMode::BHP; well_state_copy.currentProductionControls()[index_of_well_] = Well::ProducerCMode::BHP;
} }
well_state_copy.bhp()[well_copy.index_of_well_] = bhp; well_state_copy.bhp()[well_copy.index_of_well_] = bhp;
well_copy.scaleSegmentPressuresWithBhp(well_state_copy);
// initialized the well rates with the potentials i.e. the well rates based on bhp // initialized the well rates with the potentials i.e. the well rates based on bhp
const int np = number_of_phases_; const int np = number_of_phases_;
@ -898,8 +902,7 @@ namespace Opm
well_state_copy.wellRates()[well_copy.index_of_well_*np + phase] well_state_copy.wellRates()[well_copy.index_of_well_*np + phase]
= sign * well_state_copy.wellPotentials()[well_copy.index_of_well_*np + phase]; = sign * well_state_copy.wellPotentials()[well_copy.index_of_well_*np + phase];
} }
// we also want to adjust the segment rates and pressure well_copy.scaleSegmentRatesWithWellRates(well_state_copy);
well_copy.scaleSegmentRatesAndPressureWithWellRatesAndPressure(well_state_copy);
well_copy.calculateExplicitQuantities(ebosSimulator, well_state_copy, deferred_logger); well_copy.calculateExplicitQuantities(ebosSimulator, well_state_copy, deferred_logger);
const double dt = ebosSimulator.timeStepSize(); const double dt = ebosSimulator.timeStepSize();