OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-16 20:24:48 -06:00
Code Issues Packages Projects Releases Wiki Activity

incoporate keyword WINJCLN

Browse Source 
still in draft state.
This commit is contained in:
Kai Bao 2023-01-03 11:39:33 +01:00
parent 43196800b6
commit 4cb3574e0f
3 changed files with 17 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@ -768,6 +768,8 @@ namespace Opm
// Pressure drawdown (also used to determine direction of flow)
const Value drawdown = cell_press_at_perf - perf_press;
const double effectiveTw = this->isInjector() ? this->inj_fc_multiplier_[perf] * Tw : Tw;
// producing perforations
if (drawdown > 0.0) {
// Do nothing if crossflow is not allowed
@ -777,7 +779,7 @@ namespace Opm
// compute component volumetric rates at standard conditions
for (int comp_idx = 0; comp_idx < this->numComponents(); ++comp_idx) {
const Value cq_p = - Tw * (mob_perfcells[comp_idx] * drawdown);
const Value cq_p = - effectiveTw * (mob_perfcells[comp_idx] * drawdown);
cq_s[comp_idx] = b_perfcells[comp_idx] * cq_p;
}
@ -802,7 +804,7 @@ namespace Opm
}
// injection perforations total volume rates
const Value cqt_i = - Tw * (total_mob * drawdown);
const Value cqt_i = - effectiveTw * (total_mob * drawdown);
// compute volume ratio between connection and at standard conditions
Value volume_ratio = 0.0;

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

@ -240,9 +240,6 @@ namespace Opm
drawdown += skin_pressure;
}
const double effectiveTw = this->isInjector() ? this->inj_fc_multiplier_[perf] * Tw : Tw;
if (this->isInjector()) {
// std::cout << " well " << this->name() << " perf " << perf << " Tw " << Tw << " scaling " << this->inj_fc_multiplier_[perf] << " effectiveTw " << effectiveTw << std::endl;
}
// producing perforations
if (drawdown > 0) {

31
opm/simulators/wells/WellInterface_impl.hpp
View File

@ -1410,25 +1410,27 @@ namespace Opm
const auto connection = this->well_ecl_.getConnections()[perf_ecl_index];
const auto& filter_cake = this->well_ecl_.getConnections()[perf_ecl_index].getFilterCake();
if (filter_cake.active()) {
const double poro = filter_cake.poro;
const double perm = filter_cake.perm;
const double conc = this->well_ecl_.getFilterConc();
const double rw = connection.getFilterCakeRadius();
const auto cr0 = connection.r0();
const auto crw = connection.rw();
const auto cskinfactor = connection.skinFactor();
const double K = connection.Kh() / connection.connectionLength();
const double factor = filter_cake.sf_multiplier;
// we do the work here, the main thing here is to compute a multiplier for the transmissibility
if (filter_cake.geometry == Connection::FilterCakeGeometry::LINEAR) {
const double poro = filter_cake.poro;
const double perm = filter_cake.perm;
// TODO: do we want to use this rw?
const double rw = connection.getFilterCakeRadius();
const double area = connection.getFilterCakeArea();
const double conc = this->well_ecl_.getFilterConc();
const double thickness = water_inj_volume[perf] * conc / (area*(1.-poro));
std::cout << " perf " << perf << " water_injection_volume " << water_inj_volume[perf] << " conc " << conc
<< " area " << area << " poro " << poro << " thickness " << thickness << std::endl;
// TODO: this formulation might not apply for different situation
// but we are using this form just for first prototype
const double K = connection.Kh() / connection.connectionLength();
const double skin_factor = thickness / rw * K / perm;
std::cout << " sf_multiplier " << factor;
const double skin_factor = thickness / rw * K / perm * factor;
std::cout << " K " << K << " skin_factor " << skin_factor << std::endl;
const auto cr0 = connection.r0();
const auto crw = connection.rw();
const auto cskinfactor = connection.skinFactor();
const auto denom = std::log(cr0 / std::min(crw, cr0)) + cskinfactor;
const auto denom2 = std::log(cr0 / std::min(crw, cr0)) + cskinfactor + skin_factor;
const auto scaling = denom / denom2;
@ -1438,19 +1440,12 @@ namespace Opm
// TODO: basically, rescale the well connectivity index with the following formulation
// CF = angle * Kh / (std::log(r0 / std::min(rw, r0)) + skin_factor);
} else if (filter_cake.geometry == Connection::FilterCakeGeometry::RADIAL) {
const double poro = filter_cake.poro;
const double perm = filter_cake.perm;
const double conc = this->well_ecl_.getFilterConc();
const double rw = connection.getFilterCakeRadius();
const double rc = std::sqrt(rw * rw - conc * water_inj_volume[perf]/(3.1415926*(1-poro)));
std::cout << " perf " << perf << " rw " << rw << " rc " << rc;
const double K = connection.Kh() / connection.connectionLength();
std::cout << " K " << K << " perm " << perm;
const double skin_factor = K / perm * std::log(rw/rc);
std::cout << " sf_multiplier " << factor;
const double skin_factor = K / perm * std::log(rw/rc) * factor;
std::cout << " skin_factor " << skin_factor << std::endl;
const auto cr0 = connection.r0();
const auto crw = connection.rw();
const auto cskinfactor = connection.skinFactor();
const auto denom = std::log(cr0 / std::min(crw, cr0)) + cskinfactor;
const auto denom2 = std::log(cr0 / std::min(crw, cr0)) + cskinfactor + skin_factor;
const auto scaling = denom / denom2;