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

a draft version for for RADIAL model for WINJDAM

Browse Source
This commit is contained in:
Kai Bao 2023-01-02 14:36:17 +01:00
parent aa5062683d
commit 43196800b6
2 changed files with 21 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -241,7 +241,7 @@ namespace Opm
}
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;
// std::cout << " well " << this->name() << " perf " << perf << " Tw " << Tw << " scaling " << this->inj_fc_multiplier_[perf] << " effectiveTw " << effectiveTw << std::endl;
}
// producing perforations

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

@ -1437,6 +1437,26 @@ namespace Opm
std::cout << " well " << this->name() << " perf " << perf << " inj_fc_scaling " << this->inj_fc_multiplier_[perf] << std::endl;
// 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 << " 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;
std::cout << " scaling will be " << scaling << std::endl;
this->inj_fc_multiplier_[perf] = scaling;
std::cout << " well " << this->name() << " perf " << perf << " inj_fc_scaling " << this->inj_fc_multiplier_[perf] << std::endl;
}
} else {
this->inj_fc_multiplier_[perf] = 1.0;