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

Use PerfData for polymer rates

Browse Source
This commit is contained in:
Joakim Hove 2021-06-07 08:26:14 +02:00
parent 3010b96611
commit 49e9b51269
3 changed files with 4 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -407,7 +407,7 @@ namespace Opm
cq_s_poly *= this->extendEval(intQuants.polymerConcentration() * intQuants.polymerViscosityCorrection());
}
// Note. Efficiency factor is handled in the output layer
auto& perf_rate_polymer = well_state.perfRatePolymer(this->index_of_well_);
auto& perf_rate_polymer = perf_data.polymer_rates;
perf_rate_polymer[perf] = cq_s_poly.value();
cq_s_poly *= well_efficiency_factor_;

7
opm/simulators/wells/WellState.cpp
View File

@ -45,7 +45,6 @@ void WellState::base_init(const std::vector<double>& cellPressures,
this->perf_water_throughput_.clear();
this->perf_water_velocity_.clear();
this->perfRateSolvent_.clear();
this->perfRatePolymer_.clear();
this->status_.clear();
this->well_perf_data_.clear();
this->parallel_well_info_.clear();
@ -110,7 +109,6 @@ void WellState::initSingleWell(const std::vector<double>& cellPressures,
this->perf_skin_pressure_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perf_water_velocity_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perf_water_throughput_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perfRatePolymer_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perfRateSolvent_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perfRateBrine_.add(well.name(), std::vector<double>(num_perf_this_well, 0));
this->perfdata.add(well.name(), PerfData{static_cast<std::size_t>(num_perf_this_well), this->phase_usage_});
@ -737,7 +735,7 @@ void WellState::reportConnections(data::Well& well,
comp.rates.set( pi [ i ], connPI[i] );
}
if ( pu.has_polymer ) {
const auto& perf_polymer_rate = this->perfRatePolymer(well_index);
const auto& perf_polymer_rate = perf_data.polymer_rates;
comp.rates.set( rt::polymer, perf_polymer_rate[local_comp_index]);
}
if ( pu.has_brine ) {
@ -915,7 +913,8 @@ double WellState::solventWellRate(const int w) const
double WellState::polymerWellRate(const int w) const
{
const auto& perf_rates_polymer = this->perfRatePolymer(w);
const auto& perf_data = this->perfData(w);
const auto& perf_rates_polymer = perf_data.polymer_rates;
return parallel_well_info_[w]->sumPerfValues(perf_rates_polymer.begin(), perf_rates_polymer.end());
}

12
opm/simulators/wells/WellState.hpp
View File

@ -156,15 +156,6 @@ public:
/// One rate pr well
double solventWellRate(const int w) const;
/// One rate pr well connection.
std::vector<double>& perfRatePolymer(std::size_t well_index) {
return this->perfRatePolymer_[well_index];
}
const std::vector<double>& perfRatePolymer(std::size_t well_index) const {
return this->perfRatePolymer_[well_index];
}
/// One rate pr well
double polymerWellRate(const int w) const;
@ -432,9 +423,6 @@ private:
WellContainer<std::vector<double>> perfRateSolvent_;
// only for output
WellContainer<std::vector<double>> perfRatePolymer_;
// it is the throughput of water flow through the perforations
// it is used as a measure of formation damage around well-bore due to particle deposition
// it will only be used for injectors to check the injectivity