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

Merge pull request #792 from qilicun/fix-polymer-inflow

Browse Source 
look for water injector only, if not find, just continue to next step.
This commit is contained in:
Atgeirr Flø Rasmussen 2016-08-15 17:35:15 +02:00 committed by GitHub
commit c14b19c48b
1 changed files with 18 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
opm/polymer/PolymerInflow.cpp
View File

@ -73,15 +73,27 @@ namespace Opm
{
ScheduleConstPtr schedule = eclipseState->getSchedule();
for (const auto& well : schedule->getWells(currentStep)) {
WellInjectionProperties injection = well->getInjectionProperties(currentStep);
WellPolymerProperties polymer = well->getPolymerProperties(currentStep);
if (well->isInjector(currentStep)) {
WellInjectionProperties injection = well->getInjectionProperties(currentStep);
if (injection.injectorType == WellInjector::WATER) {
WellPolymerProperties polymer = well->getPolymerProperties(currentStep);
wellPolymerRate_.insert(std::make_pair(well->name(), polymer.m_polymerConcentration));
if (well->getStatus(currentStep) != WellCommon::SHUT) {
if (injection.injectorType == WellInjector::WATER) {
wellPolymerRate_.insert(std::make_pair(well->name(), polymer.m_polymerConcentration));
} else {
if (polymer.m_polymerConcentration > 0) {
OpmLog::error("Inject polymer through non-water injector '" + well->name() + "'");
}
}
} else {
OPM_THROW(std::logic_error, "For polymer injector you must have a water injector");
if (polymer.m_polymerConcentration > 0) {
OpmLog::error("Inject polymer through a shut injector '" + well->name() + "'");
}
}
}
if (well->isProducer(currentStep) && polymer.m_polymerConcentration > 0) {
OpmLog::error("Inject polymer through a producer '" + well->name() + "'");
}
}
}
@ -137,7 +149,7 @@ namespace Opm
std::fill(poly_inflow_c.begin(), poly_inflow_c.end(), 0.0);
const int nnz = sparse_inflow_.nonzeroSize();
for (int i = 0; i < nnz; ++i) {
poly_inflow_c[sparse_inflow_.nonzeroIndex(i)] = sparse_inflow_.nonzeroElement(i) ;
poly_inflow_c[sparse_inflow_.nonzeroIndex(i)] = sparse_inflow_.nonzeroElement(i);
}
}