AllWells.hpp

Go to the documentation of this file.

 
#ifndef __WELLGROUP_HEADER__
#define __WELLGROUP_HEADER__
 
// OpenCAEPoro header files
#include "ParamWell.hpp"
#include "Well.hpp"
 
using namespace std;
 
class WellGroup
{
    friend class AllWells;
 
public:
    WellGroup() = default;
    WellGroup(const string& gname)
        : name(gname){};
 
    OCP_BOOL IfReInj() const { return reInj; }
 
private:
    string      name;    
    vector<USI> wId;     
    vector<USI> wIdINJ;  
    vector<USI> wIdPROD; 
 
    // for Reinjection
    OCP_BOOL reInj{OCP_FALSE}; 
    USI      reInjPhase;       
 
    USI     prodGroup{0};
    OCP_DBL saleRate;
    OCP_DBL reInjXi;
    vector<OCP_DBL> reInjZi;
    OCP_DBL
    reInjFactor; 
 
    /*
    OCP_DBL GOPR{0}; ///< Group oil production rate
    OCP_DBL GGPR{0}; ///< Group gas production rate
    OCP_DBL GWPR{0}; ///< Group water production rate
    OCP_DBL GOPT{0}; ///< Group oil total production
    OCP_DBL GGPT{0}; ///< Group gas total production
    OCP_DBL GWPT{0}; ///< Group water total production
    OCP_DBL GGIR{0}; ///< Group gas injection rate
    OCP_DBL GWIR{0}; ///< Group water injection rate
    OCP_DBL GGIT{0}; ///< Group gas total injection
    OCP_DBL GWIT{0}; ///< Group gas total injection
    */
};
 
class AllWells
{
 
    // temp
    friend class Reservoir;
    friend class Out4RPT;
    friend class Out4VTK;
 
    friend class IsoT_FIM;
    friend class IsoT_IMPEC;
    friend class IsoT_AIMc;
    friend class IsoT_FIMn;
    friend class T_FIM;
 
public:
    AllWells() = default;
 
    // Input Param and Setup
 
public:
    void InputParam(const ParamWell& paramWell, const ParamOutput& output_param);
    void Setup(const Grid& myGrid, const Bulk& myBulk);
    void SetupWell(const Grid& myGrid, const Bulk& myBulk);
    void SetupWellGroup(const Bulk& myBulk);
    void SetupMixture(const Bulk& myBulk);
    void SetupWellBulk(Bulk& myBulk) const;
    void SetupConnWell2Bulk(const Bulk& myBulk);
 
    // General
 
public:
    void ApplyControl(const USI& i);
    void InitBHP(const Bulk& myBulk);
    void PrepareWell(const Bulk& myBulk);
    void CalTrans(const Bulk& myBulk);
    void CalFlux(const Bulk& myBulk);
    void CaldG(const Bulk& myBulk);
    void CalIPRT(const Bulk& myBulk, OCP_DBL dt);
    void CalReInjFluid(const Bulk& myBulk);
    void UpdateLastTimeStepBHP()
    {
        for (auto& w : wells) w.lbhp = w.bhp;
    }
    void ResetBHP();
    OCP_INT CheckP(const Bulk& myBulk);
    USI GetWellNum() const { return numWell; }
    string GetWellName(const USI& i) const { return wells[i].name; }
    USI GetIndex(const string& name) const;
    USI GetWellPerfNum(const USI& i) const { return wells[i].numPerf; }
    USI GetWellPerfNum() const;
    USI     GetMaxWellPerNum() const;
    void    CalMaxBHPChange();
    OCP_DBL GetdBHPmax() const { return dPmax; }
 
    OCP_DBL GetWBHP(const USI& w) const
    {
        if (wells[w].IsOpen())
            return wells[w].bhp;
        else
            return 0;
    }
    OCP_DBL GetWellDG(const USI& w, const USI& p) const { return wells[w].dG[p]; }
    void    ShowWellStatus(const Bulk& myBulk)
    {
        for (USI w = 0; w < numWell; w++) wells[w].ShowPerfStatus(myBulk);
    }
    OCP_BOOL    GetWellChange() const { return wellChange; }
    const auto& GetWell2Bulk() const { return well2bulk; }
 
protected:
    USI                     numWell;   
    vector<Well>            wells;     
    USI                     numGroup;  
    vector<WellGroup>       wellGroup; 
    vector<vector<OCP_USI>> well2bulk; 
 
    OCP_BOOL           wellChange; 
    vector<SolventINJ> solvents;   
    OCP_DBL            dPmax{0};   
 
    vector<Mixture*> flashCal;               
    OCP_DBL          Psurf{PRESSURE_STD};    
    OCP_DBL          Tsurf{TEMPERATURE_STD}; 
 
    // Injection/Production Rate
 
public:
    OCP_DBL GetFGIT() const { return FGIT; }
    OCP_DBL GetFWIR() const { return FWIR; }
    OCP_DBL GetFWIT() const { return FWIT; }
    OCP_DBL GetFOPR() const { return FOPR; }
    OCP_DBL GetFOPT() const { return FOPT; }
    OCP_DBL GetFGPR() const { return FGPR; }
    OCP_DBL GetFGPT() const { return FGPt; }
    OCP_DBL GetFWPR() const { return FWPR; }
    OCP_DBL GetFWPT() const { return FWPT; }
    OCP_DBL GetFGIR() const { return FGIR; }
 
    OCP_DBL GetWGIR(const USI& w) const { return wells[w].WGIR; }
    OCP_DBL GetWGIT(const USI& w) const { return wells[w].WGIT; }
    OCP_DBL GetWWIR(const USI& w) const { return wells[w].WWIR; }
    OCP_DBL GetWWIT(const USI& w) const { return wells[w].WWIT; }
    OCP_DBL GetWOPR(const USI& w) const { return wells[w].WOPR; }
    OCP_DBL GetWOPT(const USI& w) const { return wells[w].WOPT; }
    OCP_DBL GetWGPR(const USI& w) const { return wells[w].WGPR; }
    OCP_DBL GetWGPT(const USI& w) const { return wells[w].WGPT; }
    OCP_DBL GetWWPR(const USI& w) const { return wells[w].WWPR; }
    OCP_DBL GetWWPT(const USI& w) const { return wells[w].WWPT; }
 
protected:
    OCP_DBL FGIR{0}; 
    OCP_DBL FGIT{0}; 
    OCP_DBL FWIR{0}; 
    OCP_DBL FWIT{0}; 
    OCP_DBL FOPR{0}; 
    OCP_DBL FOPT{0}; 
    OCP_DBL FGPR{0}; 
    OCP_DBL FGPt{0}; 
    OCP_DBL FWPR{0}; 
    OCP_DBL FWPT{0}; 
 
    // for output
private:
    OCP_BOOL              useVTK{OCP_FALSE};
    vector<OCPpolyhedron> polyhedronWell;
    // When the well is under the BHP control, then give it's BHP
    // When the well is under the RATE control, then give it's RATE
    mutable vector<OCP_DBL> wellVal; 
 
public:
    void SetPolyhedronWell(const Grid& myGrid);
    void SetWellVal() const;
};
 
#endif /* end if __WELLGROUP_HEADER__ */
 
/*----------------------------------------------------------------------------*/
/*  Brief Change History of This File                                         */
/*----------------------------------------------------------------------------*/
/*  Author              Date             Actions                              */
/*----------------------------------------------------------------------------*/
/*  Shizhe Li           Oct/01/2021      Create file                          */
/*  Chensong Zhang      Oct/15/2021      Format file                          */
/*  Shizhe Li           Feb/08/2022      Rename to AllWells                   */
/*----------------------------------------------------------------------------*/