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

refactoring function updateWellProductionTargets()

Browse Source 
To handle different types of control mode.
This commit is contained in:
Kai Bao 2016-10-13 12:37:07 +02:00
parent be5f1e5c58
commit a67bff245b
3 changed files with 49 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
opm/core/wells/ProductionSpecification.cpp
View File

@ -82,6 +82,7 @@ namespace Opm
case GuideRateType::OIL : return "OIL" ; case GuideRateType::OIL : return "OIL" ;
case GuideRateType::GAS : return "GAS" ; case GuideRateType::GAS : return "GAS" ;
case GuideRateType::WATER : return "WATER" ; case GuideRateType::WATER : return "WATER" ;
case GuideRateType::LIQ : return "LIQ" ;
case GuideRateType::NONE_GRT: return "NONE_GRT"; case GuideRateType::NONE_GRT: return "NONE_GRT";
} }
OPM_THROW(std::domain_error, "Unknown guide rate type " << type << " encountered in production specification"); OPM_THROW(std::domain_error, "Unknown guide rate type " << type << " encountered in production specification");

78
opm/core/wells/WellsGroup.cpp
View File

@ -680,23 +680,39 @@ namespace Opm
// While there will be somre more complication invloved for sure. // While there will be somre more complication invloved for sure.
// Basically, we need to update the target rates for the wells still under group control. // Basically, we need to update the target rates for the wells still under group control.
ProductionSpecification::ControlMode control_mode = prodSpec().control_mode_; ProductionSpecification::ControlMode prod_mode = prodSpec().control_mode_;
double target_rate = prodSpec().liquid_max_rate_; double target_rate = -1.0;
if (ProductionSpecification::toString(control_mode) == "FLD") { switch(prod_mode) {
auto* parent_node = getParent(); case ProductionSpecification::FLD :
control_mode = parent_node->prodSpec().control_mode_; {
target_rate = parent_node->prodSpec().liquid_max_rate_; auto* parent_node = getParent();
prod_mode = parent_node->prodSpec().control_mode_;
target_rate = parent_node->getTarget(prod_mode);
break;
}
case ProductionSpecification::LRAT :
case ProductionSpecification::ORAT :
case ProductionSpecification::GRAT :
case ProductionSpecification::WRAT :
target_rate = getTarget(prod_mode);
break;
default:
OPM_THROW(std::runtime_error, "Not supporting type " << ProductionSpecification::toString(prod_mode) <<
" when updating well targets ");
} }
// the rates contributed from wells under individual control due to their own limits.
// TODO: will handle wells specified not to join group control later.
double rate_individual_control = 0.; double rate_individual_control = 0.;
for (size_t i = 0; i < children_.size(); ++i) { for (size_t i = 0; i < children_.size(); ++i) {
if (children_[i]->individualControl() && children_[i]->isProducer()) { if (children_[i]->individualControl() && children_[i]->isProducer()) {
rate_individual_control += std::abs(children_[i]->getLiquidProductionRate(well_rates)); rate_individual_control += std::abs(children_[i]->getProductionRate(well_rates, prod_mode));
} }
} }
// the rates left for the wells under group control to split
const double rate_for_group_control = target_rate - rate_individual_control; const double rate_for_group_control = target_rate - rate_individual_control;
const double my_guide_rate = productionGuideRate(true); const double my_guide_rate = productionGuideRate(true);
@ -704,7 +720,7 @@ namespace Opm
for (size_t i = 0; i < children_.size(); ++i) { for (size_t i = 0; i < children_.size(); ++i) {
if (!children_[i]->individualControl() && children_[i]->isProducer()) { if (!children_[i]->individualControl() && children_[i]->isProducer()) {
const double children_guide_rate = children_[i]->productionGuideRate(true); const double children_guide_rate = children_[i]->productionGuideRate(true);
children_[i]->applyProdGroupControl(control_mode, (children_guide_rate/my_guide_rate) * rate_for_group_control, true); children_[i]->applyProdGroupControl(prod_mode, (children_guide_rate/my_guide_rate) * rate_for_group_control, true);
children_[i]->setShouldUpdateWellTargets(false); children_[i]->setShouldUpdateWellTargets(false);
} }
} }
@ -732,25 +748,13 @@ namespace Opm
return false; return false;
} }
double WellsGroup::getLiquidProductionRate(const std::vector<double>& /*well_rates*/) const double WellsGroup::getProductionRate(const std::vector<double>& /* well_rates */,
const ProductionSpecification::ControlMode /* prod_mode */) const
{ {
// TODO: to be implemented // TODO: to be implemented
return -1.e98; return -1.e98;
} }
double WellsGroup::getOilProductionRate(const std::vector<double>& /*well_rates*/) const
{
// TODO: to be implemented
return -1.e98;
}
double WellsGroup::getWaterProductionRate(const std::vector<double>& /*well_rates*/) const
{
// TODO: to be implemented
return -1.e98;
}
// ============== WellNode members ============ // ============== WellNode members ============
@ -1133,21 +1137,23 @@ namespace Opm
double WellNode::getLiquidProductionRate(const std::vector<double>& well_rates) const double WellNode::getProductionRate(const std::vector<double>& well_rates,
const ProductionSpecification::ControlMode prod_mode) const
{ {
return ( getTotalProductionFlow(well_rates, BlackoilPhases::Liquid) + switch(prod_mode) {
getTotalProductionFlow(well_rates, BlackoilPhases::Aqua) ); case ProductionSpecification::LRAT :
return ( getTotalProductionFlow(well_rates, BlackoilPhases::Liquid) +
} getTotalProductionFlow(well_rates, BlackoilPhases::Aqua) );
case ProductionSpecification::ORAT :
double WellNode::getOilProductionRate(const std::vector<double>& well_rates) const return getTotalProductionFlow(well_rates, BlackoilPhases::Liquid);
{ case ProductionSpecification::WRAT :
return getTotalProductionFlow(well_rates, BlackoilPhases::Liquid); return getTotalProductionFlow(well_rates, BlackoilPhases::Aqua);
} case ProductionSpecification::GRAT :
return getTotalProductionFlow(well_rates, BlackoilPhases::Vapour);
double WellNode::getWaterProductionRate(const std::vector<double>& well_rates) const default:
{ OPM_THROW(std::runtime_error, "Not supporting type " << ProductionSpecification::toString(prod_mode) <<
return getTotalProductionFlow(well_rates, BlackoilPhases::Aqua); " for production rate calculation ");
}
} }
void WellNode::updateWellProductionTargets(const std::vector<double>& /*well_rates*/) void WellNode::updateWellProductionTargets(const std::vector<double>& /*well_rates*/)

21
opm/core/wells/WellsGroup.hpp
View File

@ -225,11 +225,8 @@ namespace Opm
/// Should only appy for WellNode /// Should only appy for WellNode
virtual bool isInjector() const = 0; virtual bool isInjector() const = 0;
virtual double getLiquidProductionRate(const std::vector<double>& well_rates) const = 0; virtual double getProductionRate(const std::vector<double>& well_rates,
const ProductionSpecification::ControlMode prod_mode) const = 0;
virtual double getOilProductionRate(const std::vector<double>& well_rates) const = 0;
virtual double getWaterProductionRate(const std::vector<double>& well_rates) const = 0;
virtual void updateWellProductionTargets(const std::vector<double>& well_rates) = 0; virtual void updateWellProductionTargets(const std::vector<double>& well_rates) = 0;
@ -356,11 +353,8 @@ namespace Opm
/// Should only appy for WellNode /// Should only appy for WellNode
virtual bool isInjector() const; virtual bool isInjector() const;
virtual double getLiquidProductionRate(const std::vector<double>& well_rates) const; virtual double getProductionRate(const std::vector<double>& well_rates,
const ProductionSpecification::ControlMode prod_mode) const;
virtual double getOilProductionRate(const std::vector<double>& well_rates) const;
virtual double getWaterProductionRate(const std::vector<double>& well_rates) const;
private: private:
std::vector<std::shared_ptr<WellsGroupInterface> > children_; std::vector<std::shared_ptr<WellsGroupInterface> > children_;
@ -459,11 +453,8 @@ namespace Opm
virtual bool isInjector() const; virtual bool isInjector() const;
virtual double getLiquidProductionRate(const std::vector<double>& well_rates) const; virtual double getProductionRate(const std::vector<double>& well_rates,
const ProductionSpecification::ControlMode prod_mode) const;
virtual double getOilProductionRate(const std::vector<double>& well_rates) const;
virtual double getWaterProductionRate(const std::vector<double>& well_rates) const;
virtual void updateWellProductionTargets(const std::vector<double>& well_rates); virtual void updateWellProductionTargets(const std::vector<double>& well_rates);