WellBhpThpCalculator.hpp

Go to the documentation of this file.

/*
  Copyright 2017 SINTEF Digital, Mathematics and Cybernetics.
  Copyright 2017 Statoil ASA.
  Copyright 2017 IRIS
  Copyright 2019 Norce
 
  This file is part of the Open Porous Media project (OPM).
 
  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
 
 
#ifndef OPM_WELL_BPH_THP_CALCULATOR_HEADER_INCLUDED
#define OPM_WELL_BPH_THP_CALCULATOR_HEADER_INCLUDED
 
#include <functional>
#include <optional>
#include <string>
#include <vector>
 
namespace Opm
{
 
class DeferredLogger;
class SummaryState;
class Well;
class WellInterfaceGeneric;
template<class Scalar> class WellState;
 
class WellBhpThpCalculator {
public:
    WellBhpThpCalculator(const WellInterfaceGeneric& well) : well_(well) {}
 
    bool wellHasTHPConstraints(const SummaryState& summaryState) const;
 
    double getTHPConstraint(const SummaryState& summaryState) const;
 
    double mostStrictBhpFromBhpLimits(const SummaryState& summaryState) const;
 
    double calculateThpFromBhp(const std::vector<double>& rates,
                               const double bhp,
                               const double rho,
                               const std::optional<double>& alq,
                               const double thp_limit,
                               DeferredLogger& deferred_logger) const;
 
    std::optional<double>
    computeBhpAtThpLimitProd(const std::function<std::vector<double>(const double)>& frates,
                             const SummaryState& summary_state,
                             const double maxPerfPress,
                             const double rho,
                             const double alq_value,
                             const double thp_limit,
                             DeferredLogger& deferred_logger) const;
 
    std::optional<double>
    computeBhpAtThpLimitInj(const std::function<std::vector<double>(const double)>& frates,
                            const SummaryState& summary_state,
                            const double rho,
                            const double flo_rel_tol,
                            const int max_iteration,
                            const bool throwOnError,
                            DeferredLogger& deferred_logger) const;
 
    void updateThp(const double rho,
                   const bool stop_or_zero_rate_target,
                   const std::function<double()>& alq_value,
                   const std::array<unsigned,3>& active,
                   WellState<double>& well_state,
                   const SummaryState& summary_state,
                   DeferredLogger& deferred_logger) const;
 
    template<class EvalWell>
    EvalWell calculateBhpFromThp(const WellState<double>& well_state,
                                 const std::vector<EvalWell>& rates,
                                 const Well& well,
                                 const SummaryState& summaryState,
                                 const double rho,
                                 DeferredLogger& deferred_logger) const;
 
  double calculateMinimumBhpFromThp(const WellState<double>& well_state,
                                    const Well& well,
                                    const SummaryState& summaryState,
                                    const double rho) const;
 
  bool isStableSolution(const WellState<double>& well_state,
                        const Well& well,
                        const std::vector<double>& rates,
                        const SummaryState& summaryState) const;   
 
  std::optional<double>
  estimateStableBhp (const WellState<double>& well_state,
                    const Well& well,
                    const std::vector<double>& rates,
                    const double rho,
                    const SummaryState& summaryState) const;    
 
  std::pair<double, double>
  getFloIPR(const WellState<double>& well_state,
            const Well& well, 
            const SummaryState& summary_state) const;                                                                                              
 
private:
    template<class ErrorPolicy>
    std::optional<double>
    computeBhpAtThpLimitInjImpl(const std::function<std::vector<double>(const double)>& frates,
                                const SummaryState& summary_state,
                                const double rho,
                                const double flo_rel_tol,
                                const int max_iteration,
                                DeferredLogger& deferred_logger) const;
 
    std::optional<double>
    bhpMax(const std::function<double(const double)>& fflo,
           const double bhp_limit,
           const double maxPerfPress,
           const double vfp_flo_front,
           DeferredLogger& deferred_logger) const;
 
    std::optional<double>
    computeBhpAtThpLimit(const std::function<std::vector<double>(const double)>& frates,
                         const std::function<double(const std::vector<double>)>& fbhp,
                         const std::array<double, 2>& range,
                         DeferredLogger& deferred_logger) const;
 
    double getVfpBhpAdjustment(const double bph_tab, const double thp_limit) const;
 
    bool bisectBracket(const std::function<double(const double)>& eq,
                       const std::array<double, 2>& range,
                       double& low, double& high,
                       std::optional<double>& approximate_solution,
                       DeferredLogger& deferred_logger) const;
 
    static bool bruteForceBracket(const std::function<double(const double)>& eq,
                                  const std::array<double, 2>& range,
                                  double& low, double& high,
                                  DeferredLogger& deferred_logger);
 
    double findThpFromBhpIteratively(const std::function<double(const double, const double)>& thp_func,
                                     const double bhp,
                                     const double thp_limit,
                                     const double dp,
                                     DeferredLogger& deferred_logger) const;
 
    const WellInterfaceGeneric& well_; 
};
 
}
 
#endif // OPM_WELL_BHP_THP_CALCULATOR_HEADER_INCLUDED