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

scaling the gas rates a little bit for MS well.

Browse Source 
It can possibly benefit the StandardWell, while it is something needs a
little bit more investigation before using it.
This commit is contained in:
Kai Bao 2017-09-25 14:03:10 +02:00
parent f7e1623b7a
commit 3255626b2e
3 changed files with 18 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/autodiff/MultisegmentWell_impl.hpp
View File

@ -550,9 +550,11 @@ namespace Opm
// check convergence for flux residuals
for ( int comp_idx = 0; comp_idx < numComponents(); ++comp_idx)
{
report.converged = report.converged && (maximum_residual[comp_idx] < param.tolerance_wells_ * 10.);
// report.converged = report.converged && (maximum_residual[comp_idx] < param.tolerance_wells_ * 10.);
report.converged = report.converged && (maximum_residual[comp_idx] < param.tolerance_wells_);
}
// TODO: it is not good to use a hard-coded value.
report.converged = report.converged && (maximum_residual[SPres] < 100.0);
} else { // abnormal values found and no need to check the convergence
report.converged = false;

2
opm/autodiff/SimulatorFullyImplicitBlackoilEbos.hpp
View File

@ -301,7 +301,7 @@ public:
// handling MS well related
for (const auto& well : wells_ecl) {
if (well->isMultiSegment(timer.currentStepNum())) { // there is one well is MS well
well_state.initWellStateMSWell(wells, wells_ecl, timer.currentStepNum(), prev_well_state);
well_state.initWellStateMSWell(wells, wells_ecl, timer.currentStepNum(), phaseUsage_, prev_well_state);
break;
}
}

15
opm/autodiff/WellStateFullyImplicitBlackoilDense.hpp
View File

@ -114,7 +114,7 @@ namespace Opm
/// init the MS well related.
template <typename PrevWellState>
void initWellStateMSWell(const Wells* wells, const std::vector<const Well*>& wells_ecl,
const int time_step, const PrevWellState& prev_well_state)
const int time_step, const PhaseUsage& pu, const PrevWellState& prev_well_state)
{
// still using the order in wells
const int nw = wells->number_of_wells;
@ -178,12 +178,25 @@ namespace Opm
}
}
// for the segrates_, now it becomes a recursive solution procedure.
{
const int np = numPhases();
const int start_perf = wells->well_connpos[w];
const int start_perf_next_well = wells->well_connpos[w + 1];
assert(nperf == (start_perf_next_well - start_perf)); // make sure the information from wells_ecl consistent with wells
if (pu.phase_used[Gas]) {
const int gaspos = pu.phase_pos[Gas];
// scale the phase rates for Gas to avoid too bad initial guess for gas fraction
// it will probably benefit the standard well too, while it needs to be justified
// TODO: to see if this strategy can benefit StandardWell too
// TODO: it might cause big problem for gas rate control or if there is gas rate limit
for (int perf = 0; perf < nperf; perf++) {
const int perf_pos = start_perf + perf;
perfPhaseRates()[np * perf_pos + gaspos] *= 100.;
}
}
const std::vector<double> perforation_rates(perfPhaseRates().begin() + np * start_perf,
perfPhaseRates().begin() + np * start_perf_next_well); // the perforation rates for this well
std::vector<double> segment_rates;