Grid.hpp

Go to the documentation of this file.

 
#ifndef __GRID_HEADER__
#define __GRID_HEADER__
 
// Standard header files
#include <iostream>
#include <vector>
 
// OpenCAEPoro header files
#include "CornerGrid.hpp"
#include "OCPConst.hpp"
#include "ParamReservoir.hpp"
#include "UtilOutput.hpp"
 
using namespace std;
 
 
 
 
class GPair
{
public:
    GPair() = default;
    GPair(const OCP_USI& Id, const OCP_DBL& Area)
        : id(Id)
        , area(Area){};
    static bool lessG(const GPair& G1, const GPair& G2) { return G1.id < G2.id; }
 
    OCP_USI id;   
    OCP_DBL area; 
};
 
//  Note: GB_Pair contains two variables, which indicates if a grid cell is active or
//  not and its index among active cells if it is active.
class GB_Pair
{
public:
    GB_Pair() = default;
 
    GB_Pair(bool act, OCP_USI i)
        : activity(act)
        , index(i){};
 
    bool IsAct() const { return activity; }
 
    OCP_USI GetId() const { return index; }
 
private:
    bool    activity; 
    OCP_USI index;    
};
 
//  Note: All grid cells are stored here, you can regard it as a database of the
//  reservoir. Considering there exist inactive cells (whose porosity or cell volume is
//  too small) and activeness status might change, the Grid class is necessary. The Grid
//  class is static while simulating, active grids will be stored in bulks, which is
//  "area" for calculating???
class Grid
{
    friend class Bulk;
    friend class BulkConn;
    friend class Well;
 
public:
    Grid() = default;
    void InputParam(const ParamReservoir& rs_param);
    void Setup();
 
    void SetupOrthogonalGrid();
    void SetupNeighborOrthogonalGrid();
    OCP_DBL CalAkdOrthogonalGrid(const OCP_USI& bId, const OCP_USI& eId,
                                const USI& direction);
    void CalDepthVOrthogonalGrid();
 
    void SetupCornerGrid();
    void SetupNeighborCornerGrid(const COORD& CoTmp);
    OCP_DBL CalAkdCornerGrid(const GeneralConnect& conn);
 
    void CalActiveGrid(const OCP_DBL& e1, const OCP_DBL& e2);
    const GB_Pair& MapG2B(const OCP_USI& i) const { return activeMap_G2B[i]; }
    OCP_USI GetGridNx() const { return nx; }
    OCP_USI GetGridNy() const { return ny; }
    OCP_USI GetGridNz() const { return nz; }
    OCP_USI GetGridNum() const { return numGrid; }
    OCP_USI GetConnNum() const { return numConn; }
    OCP_USI GetActiveGridNum() const { return activeGridNum; }
    OCP_USI GetActIndex(const USI& i, const USI& j, const USI& k) const;
    void GetIJKGrid(USI& i, USI& j, USI& k, const OCP_USI& n) const;
    void GetIJKBulk(USI& i, USI& j, USI& k, const OCP_USI& n) const;
    void CalSomeInfo()const;
private:
    USI     nx;      
    USI     ny;      
    USI     nz;      
    OCP_USI numGrid; 
    OCP_USI numConn; 
 
    USI                   gridType;  
    vector<vector<GPair>> gNeighbor; 
 
    // Orthogonal grid
    vector<OCP_DBL> tops;  
    vector<OCP_DBL> depth; 
    vector<OCP_DBL> dx;    
    vector<OCP_DBL> dy;    
    vector<OCP_DBL> dz;    
 
    // Corner-point grid
    vector<OCP_DBL> coord; 
    vector<OCP_DBL> zcorn; 
 
    // Rock properties
    vector<OCP_DBL> v;    
    vector<OCP_DBL> ntg;  
    vector<OCP_DBL> poro; 
    vector<OCP_DBL> kx;   
    vector<OCP_DBL> ky;   
    vector<OCP_DBL> kz;   
 
    // Initial Properties
    vector<OCP_DBL>   SwatInit; 
 
    // Region
    vector<USI> SATNUM; 
    vector<USI> PVTNUM; 
    vector<USI> ACTNUM; 
 
    // Active grid cells
    OCP_USI         activeGridNum; 
    vector<OCP_USI> activeMap_B2G; 
    vector<GB_Pair> activeMap_G2B; 
 
private:
    // Auxiliary variable
    USI             numDigutIJK;  
 
public:
    void CalNumDigutIJK(); 
    USI GetNumDigitIJK() const { return numDigutIJK; }
};
 
#endif /* end if __GRID_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           Nov/18/2021      Add Connections between Grids        */
/*  Chensong Zhang      Jan/16/2022      Finish Doxygen                       */
/*----------------------------------------------------------------------------*/