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

using a small value drawdown for connection

Browse Source 
when crossflow is banned from the well, while the crossflow is happening
for this connection for mutlisegment well.
This commit is contained in:
Kai Bao
2018-01-30 15:32:07 +01:00
parent 04c35c8916
commit 972fed602e
1 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
opm/autodiff/MultisegmentWell_impl.hpp
View File

@@ -1021,7 +1021,8 @@ namespace Opm
if ( drawdown > 0.0) { if ( drawdown > 0.0) {
// Do nothing is crossflow is not allowed // Do nothing is crossflow is not allowed
if (!allow_cf && well_type_ == INJECTOR) { if (!allow_cf && well_type_ == INJECTOR) {
return; // TODO: actaully, we did not have a case to test if it is the same case for injectors
drawdown = -mini_drawdown;
} }
// compute component volumetric rates at standard conditions // compute component volumetric rates at standard conditions
@@ -1041,7 +1042,7 @@ namespace Opm
} else { // injecting perforations } else { // injecting perforations
// Do nothing if crossflow is not allowed // Do nothing if crossflow is not allowed
if (!allow_cf && well_type_ == PRODUCER) { if (!allow_cf && well_type_ == PRODUCER) {
return; drawdown = mini_drawdown;
} }
// for injecting perforations, we use total mobility // for injecting perforations, we use total mobility
@@ -1896,10 +1897,11 @@ namespace Opm
// 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?
// we do not use it in the simulation for now, while we might need them if
// we handle the pressure in SEG mode.
well_state.perfPhaseRates()[(first_perf_ + perf) * number_of_phases_ + ebosCompIdxToFlowCompIdx(comp_idx)] = cq_s[comp_idx].value();
} }
// TODO: we should save the perforation pressure and preforation rates?
// we do not use it in the simulation for now, while we might need them if
// we handle the pressure in SEG mode.
} }
// the fourth dequation, the pressure drop equation // the fourth dequation, the pressure drop equation