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) {};
 
private:
 
    string name; 
    vector<USI> wId; 
    vector<USI> wIdINJ; 
    vector<USI> wIdPROD; 
    
    // for Reinjection
    bool reInj{false}; 
    USI injPhase; 
    USI prodGroup{0}; 
    OCP_DBL saleRate;
    OCP_DBL xi;
    vector<OCP_DBL> zi;
    OCP_DBL factor; 
 
    OCP_DBL GOPR{0}; 
    OCP_DBL GGPR{0}; 
    OCP_DBL GWPR{0}; 
    OCP_DBL GOPT{0}; 
    OCP_DBL GGPT{0}; 
    OCP_DBL GWPT{0}; 
    OCP_DBL GGIR{0}; 
    OCP_DBL GWIR{0}; 
    OCP_DBL GGIT{0}; 
    OCP_DBL GWIT{0}; 
};
 
 
class AllWells
{
 
    // temp
    friend class Reservoir;
    friend class DetailInfo;
 
public:
    AllWells() = default;
 
    // General
 
public:
    void InputParam(const ParamWell& paramWell);
    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 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 CalProdWeight(const Bulk& myBulk);
    void CaldG(const Bulk& myBulk);
    void CalIPRT(const Bulk& myBulk, OCP_DBL dt);
    void CalReInjFluid(const Bulk& myBulk);
    void AllocateMat(LinearSystem& myLS, const USI& bulknum) const;
    void UpdateLastBHP() { for (auto& w : wells) w.lBHP = w.BHP; }
    void ResetBHP();
    void UpdateLastDg()
    {
        for (auto& w : wells) w.ldG = w.dG;
    }
    void ResetDg()
    {
        for (auto& w : wells) w.dG = w.ldG;
    }
    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 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 GetFGIT() const { return FGIT; }
    OCP_DBL GetFWIR() const { return FWIR; }
    OCP_DBL GetFWIT() const { return FWIT; }
    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; }
    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 GetWBHP(const USI& w) const { 
        if (wells[w].WellState())
            return wells[w].BHP;
        else
            return 0;
    }
    OCP_DBL GetWellDg(const USI& w, const USI& p) const { return wells[w].dG[p]; }
    OCP_DBL CalWellQT();
    void ShowWellStatus(const Bulk& myBulk) { for (USI w = 0; w < numWell; w++) wells[w].ShowPerfStatus(myBulk); }
    bool GetWellChange() const { return wellChange; }
 
private:
    USI          numWell;   
    vector<Well> wells; 
    USI          numGroup; 
    vector<WellGroup> wellGroup; 
 
    bool         wellChange; 
    vector<SolventINJ> solvents; 
    OCP_DBL            dPmax{0}; 
 
    vector<Mixture*> flashCal; 
 
    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}; 
 
    OCP_DBL QT{0}; 
 
    // IMPEC
 
public:
    void CalCFL(const Bulk& myBulk, const OCP_DBL& dt) const;
    void MassConserveIMPEC(Bulk& myBulk, OCP_DBL dt);
    void AssemblaMatIMPEC(LinearSystem& myLS, const Bulk& myBulk,
                          const OCP_DBL& dt) const;
    void GetSolIMPEC(const vector<OCP_DBL>& u, const OCP_USI& bId);
 
    // FIM
 
public:
    void AssemblaMatFIM(LinearSystem& myLS, const Bulk& myBulk,
                        const OCP_DBL& dt) const;
    void GetSolFIM(const vector<OCP_DBL>& u, const OCP_USI& bId, const USI& len);
    void GetSol01FIM(const vector<OCP_DBL>& u, const OCP_USI& bId, const USI& len,
                     const OCP_DBL& alpha);
    void CalResFIM(ResFIM& resFIM, const Bulk& myBulk, const OCP_DBL& dt) const;
    void ShowRes(const vector<OCP_DBL>& res, const Bulk& myBulk) const;
 
 
    // FIM
    void AssemblaMatFIM_new(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
    void AssemblaMatFIM_new_n(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
 
    // AIMt, AIMs
 
public:
    void AssemblaMatAIMt(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) const;
    void CalResAIMt(ResFIM& resFIM, const Bulk& myBulk, const OCP_DBL& dt) const;
    void GetSolAIMt(const vector<OCP_DBL>& u, const OCP_USI& bId, const USI& len);
 
    void AssemblaMatAIMs(LinearSystem& myLS, const Bulk& myBulk,
        const OCP_DBL& dt) 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                   */
/*----------------------------------------------------------------------------*/