OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge pull request #777 from joakim-hove/well2-prediction-mode

Browse Source 
Add method bool Well2::predictionMode()
This commit is contained in:
Joakim Hove
2019-05-21 06:52:59 +02:00
committed by GitHub
parent 272d0fc9a1 8231b84a11
commit 16dabd5dc4
3 changed files with 73 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
opm/parser/eclipse/EclipseState/Schedule/Well/Well2.hpp
View File

@@ -68,6 +68,7 @@ public:
double getGuideRateScalingFactor() const;
std::size_t firstTimeStep() const;
bool predictionMode() const;
bool canOpen() const;
bool isProducer() const;
size_t seqIndex() const;
@@ -106,6 +107,7 @@ public:
*/
std::map<int, std::vector<Connection>> getCompletions() const;
bool updatePrediction(bool prediction_mode);
bool updateAutoShutin(bool auto_shutin);
bool updateCrossFlow(bool allow_cross_flow);
bool updateHead(int I, int J);
@@ -153,6 +155,7 @@ private:
WellGuideRate guide_rate;
double efficiency_factor;
double solvent_fraction;
bool prediction_mode = true;
std::shared_ptr<const WellEconProductionLimits> econ_limits;
std::shared_ptr<const WellPolymerProperties> polymer_properties;

62
src/opm/parser/eclipse/EclipseState/Schedule/Schedule.cpp
View File

@@ -713,6 +713,7 @@ namespace Opm {
auto well2 = std::make_shared<Well2>(*dynamic_state[currentStep]);
bool switching_from_injector = !well2->isProducer();
auto properties = std::make_shared<WellProductionProperties>(well2->getProductionProperties());
bool update_well = false;
properties->handleWCONHIST(record);
if (switching_from_injector) {
@@ -720,11 +721,21 @@ namespace Opm {
well2->updateProducer(true);
}
if (well2->updateProduction(properties) || switching_from_injector) {
if (well2->updateProducer(true))
update_well = true;
if (well2->updateProduction(properties))
update_well = true;
if (well2->updatePrediction(false))
update_well = true;
if (update_well) {
m_events.addEvent( ScheduleEvents::PRODUCTION_UPDATE , currentStep);
this->addWellEvent( well2->name(), ScheduleEvents::PRODUCTION_UPDATE, currentStep);
this->updateWell(well2, currentStep);
}
if ( !well2->getAllowCrossFlow() && (properties->OilRate + properties->WaterRate + properties->GasRate) == 0 ) {
std::string msg =
"Well " + well2->name() + " is a history matched well with zero rate where crossflow is banned. " +
@@ -779,9 +790,9 @@ namespace Opm {
auto well2 = std::make_shared<Well2>(*dynamic_state[currentStep]);
bool switching_from_injector = !well2->isProducer();
auto properties = std::make_shared<WellProductionProperties>(well2->getProductionProperties());
bool update_well = switching_from_injector;
properties->clearControls();
if (well2->isAvailableForGroupControl())
properties->addProductionControl(WellProducer::GRUP);
@@ -790,7 +801,16 @@ namespace Opm {
if (switching_from_injector)
properties->resetDefaultBHPLimit();
if (well2->updateProduction(properties)) {
if (well2->updateProducer(true))
update_well = true;
if (well2->updateProduction(properties))
update_well = true;
if (well2->updatePrediction(true))
update_well = true;
if (update_well) {
m_events.addEvent( ScheduleEvents::PRODUCTION_UPDATE , currentStep);
this->addWellEvent( well2->name(), ScheduleEvents::PRODUCTION_UPDATE, currentStep);
this->updateWell(well2, currentStep);
@@ -898,8 +918,15 @@ namespace Opm {
auto injection = std::make_shared<WellInjectionProperties>(well2->getInjectionProperties());
injection->handleWCONINJE(record, well2->isAvailableForGroupControl(), well_name, section.unitSystem());
update_well = well2->updateProducer(false);
update_well |= well2->updateInjection(injection);
if (well2->updateProducer(false))
update_well = true;
if (well2->updateInjection(injection))
update_well = true;
if (well2->updatePrediction(true))
update_well = true;
if (update_well) {
this->updateWell(well2, currentStep);
m_events.addEvent( ScheduleEvents::INJECTION_UPDATE , currentStep );
@@ -957,8 +984,15 @@ namespace Opm {
auto injection = std::make_shared<WellInjectionProperties>(well2->getInjectionProperties());
injection->handleWCONINJH(record, well2->isProducer(), well_name, section.unitSystem());
update_well = well2->updateProducer(false);
update_well |= well2->updateInjection(injection);
if (well2->updateProducer(false))
update_well = true;
if (well2->updateInjection(injection))
update_well = true;
if (well2->updatePrediction(false))
update_well = true;
if (update_well) {
this->updateWell(well2, currentStep);
m_events.addEvent( ScheduleEvents::INJECTION_UPDATE , currentStep );
@@ -2731,7 +2765,6 @@ namespace Opm {
const auto& well2 = *(dynamic_state[last_step]);
bool equal = true;
if (well2.getWellConnectionOrdering() != well.getWellConnectionOrdering())
equal = false;
@@ -2763,6 +2796,19 @@ namespace Opm {
for (std::size_t step = first_step; step < this->size(); step++) {
const auto& well2 = *(dynamic_state[step]);
{
bool predictionMode;
if (well.isProducer(step))
predictionMode = well.getProductionProperties(step).predictionMode;
else
predictionMode = well.getInjectionProperties(step).predictionMode;
if (well2.predictionMode() != predictionMode) {
std::string msg = "Predictionmode error for well:" + well2.name() + " at " + std::to_string(step);
printf("Well:%d Well2:%d \n", predictionMode, well2.predictionMode());
parseContext.handleError(ParseContext::SCHEDULE_WELL_ERROR, msg, errors);
}
}
{
const auto& conn1 = well.getConnections(step);
const auto& conn2 = well2.getConnections();

16
src/opm/parser/eclipse/EclipseState/Schedule/Well/Well2.cpp
View File

@@ -599,6 +599,22 @@ bool Well2::canOpen() const {
}
}
bool Well2::predictionMode() const {
return this->prediction_mode;
}
bool Well2::updatePrediction(bool prediction_mode) {
if (this->prediction_mode != prediction_mode) {
this->prediction_mode = prediction_mode;
return true;
}
return false;
}
WellCompletion::CompletionOrderEnum Well2::getWellConnectionOrdering() const {
return this->ordering;
}