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

updating perforation pressure and rates for MSW well

Browse Source
This commit is contained in:
Kai Bao 2019-04-05 10:01:16 +02:00
parent 32feeefae3
commit e2b8edd968
2 changed files with 31 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/autodiff/MultisegmentWell.hpp
View File

@ -295,14 +295,15 @@ namespace Opm
// basically Q_p / \sigma_p Q_p // basically Q_p / \sigma_p Q_p
EvalWell surfaceVolumeFraction(const int seg, const int comp_idx) const; EvalWell surfaceVolumeFraction(const int seg, const int comp_idx) const;
void computePerfRate(const IntensiveQuantities& int_quants, void computePerfRatePressure(const IntensiveQuantities& int_quants,
const std::vector<EvalWell>& mob_perfcells, const std::vector<EvalWell>& mob_perfcells,
const int seg, const int seg,
const int perf, const int perf,
const EvalWell& segment_pressure, const EvalWell& segment_pressure,
const bool& allow_cf, const bool& allow_cf,
std::vector<EvalWell>& cq_s, std::vector<EvalWell>& cq_s,
Opm::DeferredLogger& deferred_logger) const; EvalWell& perf_press,
Opm::DeferredLogger& deferred_logger) const;
// convert a Eval from reservoir to contain the derivative related to wells // convert a Eval from reservoir to contain the derivative related to wells
EvalWell extendEval(const Eval& in) const; EvalWell extendEval(const Eval& in) const;

34
opm/autodiff/MultisegmentWell_impl.hpp
View File

@ -971,14 +971,15 @@ namespace Opm
template <typename TypeTag> template <typename TypeTag>
void void
MultisegmentWell<TypeTag>:: MultisegmentWell<TypeTag>::
computePerfRate(const IntensiveQuantities& int_quants, computePerfRatePressure(const IntensiveQuantities& int_quants,
const std::vector<EvalWell>& mob_perfcells, const std::vector<EvalWell>& mob_perfcells,
const int seg, const int seg,
const int perf, const int perf,
const EvalWell& segment_pressure, const EvalWell& segment_pressure,
const bool& allow_cf, const bool& allow_cf,
std::vector<EvalWell>& cq_s, std::vector<EvalWell>& cq_s,
Opm::DeferredLogger& deferred_logger) const EvalWell& perf_press,
Opm::DeferredLogger& deferred_logger) const
{ {
std::vector<EvalWell> cmix_s(num_components_, 0.0); std::vector<EvalWell> cmix_s(num_components_, 0.0);
@ -1011,9 +1012,10 @@ namespace Opm
// pressure difference between the perforation and the grid cell // pressure difference between the perforation and the grid cell
const double cell_perf_press_diff = cell_perforation_pressure_diffs_[perf]; const double cell_perf_press_diff = cell_perforation_pressure_diffs_[perf];
perf_press = pressure_cell + cell_perf_press_diff;
// Pressure drawdown (also used to determine direction of flow) // Pressure drawdown (also used to determine direction of flow)
// TODO: not 100% sure about the sign of the seg_perf_press_diff // TODO: not 100% sure about the sign of the seg_perf_press_diff
const EvalWell drawdown = (pressure_cell + cell_perf_press_diff) - (segment_pressure + perf_seg_press_diff); const EvalWell drawdown = perf_press - (segment_pressure + perf_seg_press_diff);
// producing perforations // producing perforations
if ( drawdown > 0.0) { if ( drawdown > 0.0) {
@ -1861,7 +1863,15 @@ namespace Opm
std::vector<EvalWell> mob(num_components_, 0.0); std::vector<EvalWell> mob(num_components_, 0.0);
getMobility(ebosSimulator, perf, mob); getMobility(ebosSimulator, perf, mob);
std::vector<EvalWell> cq_s(num_components_, 0.0); std::vector<EvalWell> cq_s(num_components_, 0.0);
computePerfRate(int_quants, mob, seg, perf, seg_pressure, allow_cf, cq_s, deferred_logger); EvalWell perf_press;
computePerfRatePressure(int_quants, mob, seg, perf, seg_pressure, allow_cf, cq_s, perf_press, deferred_logger);
// store the perf pressure and rates
const int rate_start_offset = (first_perf_ + perf) * number_of_phases_;
for (int comp_idx = 0; comp_idx < num_components_; ++comp_idx) {
well_state.perfPhaseRates()[rate_start_offset + ebosCompIdxToFlowCompIdx(comp_idx)] = cq_s[comp_idx].value();
}
well_state.perfPress()[first_perf_ + perf] = perf_press.value();
for (int comp_idx = 0; comp_idx < num_components_; ++comp_idx) { for (int comp_idx = 0; comp_idx < num_components_; ++comp_idx) {
// the cq_s entering mass balance equations need to consider the efficiency factors. // the cq_s entering mass balance equations need to consider the efficiency factors.
@ -1877,14 +1887,14 @@ namespace Opm
// also need to consider the efficiency factor when manipulating the jacobians. // also need to consider the efficiency factor when manipulating the jacobians.
duneC_[seg][cell_idx][pv_idx][comp_idx] -= cq_s_effective.derivative(pv_idx + numEq); // intput in transformed matrix duneC_[seg][cell_idx][pv_idx][comp_idx] -= cq_s_effective.derivative(pv_idx + numEq); // intput in transformed matrix
// the index name for the D should be eq_idx / pv_idx // the index name for the D should be eq_idx / pv_idx
duneD_[seg][seg][comp_idx][pv_idx] -= cq_s_effective.derivative(pv_idx + numEq); duneD_[seg][seg][comp_idx][pv_idx] -= cq_s_effective.derivative(pv_idx + numEq);
} }
for (int pv_idx = 0; pv_idx < numEq; ++pv_idx) { for (int pv_idx = 0; pv_idx < numEq; ++pv_idx) {
// also need to consider the efficiency factor when manipulating the jacobians. // also need to consider the efficiency factor when manipulating the jacobians.
duneB_[seg][cell_idx][comp_idx][pv_idx] -= cq_s_effective.derivative(pv_idx); duneB_[seg][cell_idx][comp_idx][pv_idx] -= cq_s_effective.derivative(pv_idx);
} }
} }
// TODO: we should save the perforation pressure and preforation rates? // TODO: we should save the perforation pressure and preforation rates?