ParamWell Class Reference

#include <ParamWell.hpp>

Public Member Functions
void	Init ()
	Initialize the inputting the params of wells.
void	InitTime ()
	Initialize the critical time.
void	InputWELSPECS (ifstream &ifs)
void	InputCOMPDAT (ifstream &ifs)
void	InputWCONINJE (ifstream &ifs)
void	InputWCONPROD (ifstream &ifs)
void	InputTSTEP (ifstream &ifs)
void	InputWELTARG (ifstream &ifs)
void	InputWTEMP (ifstream &ifs)
	Input the temperature of injected fluid.
void	InputUNWEIGHT (ifstream &ifs)
	Input injector type – MOBWEIGHT(defaulted) or UNWEIGHT.
void	InputWELLSTRE (ifstream &ifs)
void	InputPSURF (ifstream &ifs)
	Input surface pressure.
void	InputTSURF (ifstream &ifs)
	Input surface temperature.
void	CheckParam (const OCP_BOOL &boModel) const
	Check if wrong params are input.
void	CheckPerf () const
	Check if params of Perforation is wrong.
void	CheckINJFluid () const

Public Attributes
vector< WellParam >	well
	Contains all the information of wells.
vector< OCP_DBL >	criticalTime
	Records the critical time given by users.
vector< Solvent >	solSet
	Sets of Solvent.
OCP_DBL	Psurf {PRESSURE_STD}
	Pressure in surface condition.
OCP_DBL	Tsurf {TEMPERATURE_STD}
	Temperature in surface condition.

Detailed Description

ParamWell is an internal structure used to stores the information of wells from input files. It is an intermediate interface and independent of the main simulator. After all file inputting finishes, the params in it will pass to corresponding modules.

Definition at line 99 of file ParamWell.hpp.

Member Function Documentation

InputCOMPDAT()

void ParamWell::InputCOMPDAT

(

ifstream &

ifs

)

Input the well keyword COMPDAT. COMPDAT contains the information of perforations of wells, for example the location, the trans or directions.

Definition at line 80 of file ParamWell.cpp.

{
    USI            num = well.size();
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        DealDefault(vbuf);
        string            src   = vbuf[0];
        string::size_type pos   = src.find("*");
        OCP_BOOL          match = (pos != string::npos);
        if (match) {
            src.erase(pos);
        }
        OCP_BOOL tmp = OCP_FALSE;
 
        for (USI w = 0; w < num; w++) {
            if (match)
                tmp = (well[w].name.substr(0, pos) == src);
            else
                tmp = (well[w].name == src);
 
            if (tmp) {
 
                USI k1 = stoi(vbuf[3]);
                USI k2 = stoi(vbuf[4]);
 
                for (USI k = k1; k <= k2; k++) {
                    if (vbuf[1] == "DEFAULT" || vbuf[2] == "DEFAULT") {
                        well[w].I_perf.push_back(well[w].I);
                        well[w].J_perf.push_back(well[w].J);
                    } else {
                        well[w].I_perf.push_back(stoi(vbuf[1]));
                        well[w].J_perf.push_back(stoi(vbuf[2]));
                    }
                    well[w].K_perf.push_back(k);
 
                    if (vbuf[5] != "DEFAULT")
                        well[w].WI.push_back(stod(vbuf[5]));
                    else
                        well[w].WI.push_back(-1.0);
 
                    if (vbuf[6] != "DEFAULT")
                        well[w].diameter.push_back(stod(vbuf[6]));
                    else
                        well[w].diameter.push_back(1.0);
 
                    if (vbuf[7] != "DEFAULT")
                        well[w].kh.push_back(stod(vbuf[7]));
                    else
                        well[w].kh.push_back(-1.0);
 
                    if (vbuf[8] != "DEFAULT")
                        well[w].skinFactor.push_back(stod(vbuf[8]));
                    else
                        well[w].skinFactor.push_back(0.0);
 
                    if (vbuf[9] != "DEFAULT")
                        well[w].direction.push_back(vbuf[9]);
                    else
                        well[w].direction.push_back("z");
                }
            }
        }
    }
    // cout << "COMPDAT" << endl;
}

References DealDefault(), ReadLine(), and well.

InputTSTEP()

void ParamWell::InputTSTEP

(

ifstream &

ifs

)

Input the keyword: TSTEP. TSTEP is used to divide the simulation time and ususally the time point is critical, at which for example, the operation mode of well will change. So the params of solving equations could be adjusted correspondingly.

Definition at line 214 of file ParamWell.cpp.

{
    assert(criticalTime.size() > 0);
 
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        DealDefault(vbuf);
        OCP_INT len = vbuf.size();
        for (OCP_INT i = 0; i < len - 1; i++) {
            OCP_DBL t = criticalTime.back() + stod(vbuf[i]);
            criticalTime.push_back(t);
        }
        if (vbuf.back() != "/") {
            OCP_DBL t = criticalTime.back() + stod(vbuf.back());
            criticalTime.push_back(t);
        }
    }
}

References criticalTime, DealDefault(), and ReadLine().

InputWCONINJE()

void ParamWell::InputWCONINJE

(

ifstream &

ifs

)

Input the well keyword WCONINJE. WCONINJE describes the initial operation mode of injection well.

Definition at line 148 of file ParamWell.cpp.

{
    assert(criticalTime.size() > 0);
 
    const USI      d   = criticalTime.size() - 1;
    const USI      num = well.size();
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        DealDefault(vbuf);
        string            src        = vbuf[0];
        string::size_type pos        = src.find("*");
        const OCP_BOOL    fuzzyMatch = (pos != string::npos);
        if (fuzzyMatch) {
            src.erase(pos);
        }
 
        if (fuzzyMatch) {
            for (USI w = 0; w < num; w++) {
                if (well[w].name.find(src) != string::npos) {
                    well[w].optParam.push_back(WellOptPair(d, "INJ", vbuf));
                }
            }
        } else {
            for (USI w = 0; w < num; w++) {
                if (well[w].name == src) {
                    well[w].optParam.push_back(WellOptPair(d, "INJ", vbuf));
                }
            }
        }
    }
    // cout << "WCONINJE" << endl;
}

References criticalTime, DealDefault(), ReadLine(), and well.

InputWCONPROD()

void ParamWell::InputWCONPROD

(

ifstream &

ifs

)

Input the well keyword WCONPROD. WCONPROD describes the initial operation mode of production well.

Definition at line 183 of file ParamWell.cpp.

{
    assert(criticalTime.size() > 0);
 
    const USI      d   = criticalTime.size() - 1;
    const USI      num = well.size();
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        DealDefault(vbuf);
        string            src        = vbuf[0];
        string::size_type pos        = src.find("*");
        const OCP_BOOL    fuzzyMatch = (pos != string::npos);
        if (fuzzyMatch) {
            src.erase(pos);
        }
 
        if (fuzzyMatch) {
            for (USI w = 0; w < num; w++)
                if (well[w].name.find(src) != string::npos)
                    well[w].optParam.push_back(WellOptPair(d, "PROD", vbuf));
        } else {
            for (USI w = 0; w < num; w++)
                if (well[w].name == src)
                    well[w].optParam.push_back(WellOptPair(d, "PROD", vbuf));
        }
    }
    // cout << "WCONPROD" << endl;
}

References criticalTime, DealDefault(), ReadLine(), and well.

InputWELLSTRE()

void ParamWell::InputWELLSTRE

(

ifstream &

ifs

)

Input well keyword: Solvent. It describes the molar fraction of components of fluid injected to reservoir from INJ.

Definition at line 373 of file ParamWell.cpp.

{
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
        solSet.push_back(Solvent(vbuf));
    }
    // cout << "WELLSTRE" << endl;
}

References ReadLine(), and solSet.

InputWELSPECS()

void ParamWell::InputWELSPECS

(

ifstream &

ifs

)

Input the well keyword WELSPECS. WELSPECS defines wells including well name, well location, well depth and son on.

Definition at line 68 of file ParamWell.cpp.

{
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        DealDefault(vbuf);
        well.push_back(WellParam(vbuf));
    }
    // cout << "WELSPECS" << endl;
}

References DealDefault(), ReadLine(), and well.

InputWELTARG()

void ParamWell::InputWELTARG

(

ifstream &

ifs

)

Input the well keyword: WELTARG. WELTARG is used to change the operation mode of well anytime. For example, the oil production rate is changed from 1000 stb/day to 1500 stb/day at the 100th day.

Definition at line 235 of file ParamWell.cpp.

{
    assert(criticalTime.size() > 0);
 
    cout << "\n---------------------" << endl
         << "WELTARG"
         << "\n---------------------" << endl;
 
    const USI      d   = criticalTime.size() - 1;
    const USI      num = well.size();
    vector<string> vbuf;
    while (ReadLine(ifs, vbuf)) {
        if (vbuf[0] == "/") break;
 
        for (const auto& s : vbuf) {
            cout << s << "   ";
        }
        cout << endl;
 
        string            src        = vbuf[0];
        string::size_type pos        = src.find("*");
        const OCP_BOOL    fuzzyMatch = (pos != string::npos);
        if (fuzzyMatch) {
            src.erase(pos);
        }
        OCP_BOOL succMatch = OCP_FALSE;
 
        for (USI w = 0; w < num; w++) {
            if (fuzzyMatch)
                succMatch = (well[w].name.find(src) != string::npos);
            else
                succMatch = (well[w].name == src);
 
            if (succMatch) {
                if (well[w].optParam.size() == 0) {
                    OCP_ABORT("No Well Control Defined in Advance!");
                }
                WellOptPair tar = well[w].optParam.back();
                tar.d           = d;
                tar.opt.optMode = vbuf[1];
                OCP_DBL val     = stod(vbuf[2]);
                if (vbuf[1] == "BHP") {
                    if (tar.opt.type == "INJ")
                        tar.opt.maxBHP = val;
                    else
                        tar.opt.minBHP = val;
                } else {
                    tar.opt.maxRate = val;
                }
                well[w].optParam.push_back(tar);
            }
        }
    }
}

References criticalTime, ReadLine(), and well.

---

The documentation for this class was generated from the following files:

• ParamWell.hpp
• ParamWell.cpp