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

Make the computation of the productivity index only when asked for

Browse Source
This commit is contained in:
Tor Harald Sandve 2018-11-08 13:28:52 +01:00
parent 3bc292d168
commit 8dc874dda4
1 changed files with 15 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/autodiff/StandardWell_impl.hpp
View File

@ -437,6 +437,9 @@ namespace Opm
WellState& well_state, WellState& well_state,
bool only_wells) bool only_wells)
{ {
const Opm::SummaryConfig summaryConfig = ebosSimulator.vanguard().summaryConfig();
const int np = number_of_phases_; const int np = number_of_phases_;
// clear all entries // clear all entries
@ -606,13 +609,19 @@ namespace Opm
// Store the perforation pressure for later usage. // Store the perforation pressure for later usage.
well_state.perfPress()[first_perf_ + perf] = well_state.bhp()[index_of_well_] + perf_pressure_diffs_[perf]; well_state.perfPress()[first_perf_ + perf] = well_state.bhp()[index_of_well_] + perf_pressure_diffs_[perf];
// Compute Productivity index // Compute Productivity index if asked for
const auto& pu = phaseUsage();
for (int p = 0; p < np; ++p) { for (int p = 0; p < np; ++p) {
const unsigned int compIdx = flowPhaseToEbosCompIdx(p); if ( (pu.phase_pos[Water] == p && (summaryConfig.hasSummaryKey("WPIW:" + name()) || summaryConfig.hasSummaryKey("WPIL:" + name())))
const double drawdown = well_state.perfPress()[first_perf_ + perf] - intQuants.fluidState().pressure(FluidSystem::oilPhaseIdx).value(); || (pu.phase_pos[Oil] == p && (summaryConfig.hasSummaryKey("WPIO:" + name()) || summaryConfig.hasSummaryKey("WPIL:" + name())))
double productivity_index = cq_s[compIdx].value() / drawdown; || (pu.phase_pos[Gas] == p && summaryConfig.hasSummaryKey("WPIG:" + name()))) {
scaleProductivityIndex(perf, productivity_index);
well_state.productivityIndex()[np*index_of_well_ + p] += productivity_index; const unsigned int compIdx = flowPhaseToEbosCompIdx(p);
const double drawdown = well_state.perfPress()[first_perf_ + perf] - intQuants.fluidState().pressure(FluidSystem::oilPhaseIdx).value();
double productivity_index = cq_s[compIdx].value() / drawdown;
scaleProductivityIndex(perf, productivity_index);
well_state.productivityIndex()[np*index_of_well_ + p] += productivity_index;
}
} }
} }