OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

added water evaporation as function of pressure and salt concentration

Browse Source
This commit is contained in:
Paul Egberts
2022-10-20 11:16:37 +02:00
parent ef7d1a5441
commit 58da92b420
8 changed files with 161 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

3
opm/material/fluidsystems/BlackOilFluidSystem.hpp
View File

@@ -1231,10 +1231,11 @@ public:
const auto& p = decay<LhsEval>(fluidState.pressure(phaseIdx));
const auto& T = decay<LhsEval>(fluidState.temperature(phaseIdx));
const auto& saltConcentration = decay<LhsEval>(fluidState.saltConcentration());
switch (phaseIdx) {
case oilPhaseIdx: return 0.0;
case gasPhaseIdx: return gasPvt_->saturatedWaterVaporizationFactor(regionIdx, T, p);
case gasPhaseIdx: return gasPvt_->saturatedWaterVaporizationFactor(regionIdx, T, p, saltConcentration);
case waterPhaseIdx: return 0.0;
default: throw std::logic_error("Unhandled phase index "+std::to_string(phaseIdx));
}

10
opm/material/fluidsystems/blackoilpvt/Co2GasPvt.hpp
View File

@@ -209,6 +209,16 @@ public:
const Evaluation& /*pressure*/) const
{ return 0.0; /* this is non-humid gas! */ }
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of water saturated gas.
*/
template <class Evaluation = Scalar>
Evaluation saturatedWaterVaporizationFactor(unsigned /*regionIdx*/,
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/,
const Evaluation& /*saltConcentration*/) const
{ return 0.0; }
/*!
* \brief Returns the oil vaporization factor \f$R_v\f$ [m^3/m^3] of the oil phase.
*/

11
opm/material/fluidsystems/blackoilpvt/DryGasPvt.hpp
View File

@@ -279,6 +279,17 @@ public:
const Evaluation& /*pressure*/) const
{ return 0.0; /* this is non-humid gas! */ }
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of water saturated gas.
*/
template <class Evaluation = Scalar>
Evaluation saturatedWaterVaporizationFactor(unsigned /*regionIdx*/,
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/,
const Evaluation& /*saltConcentration*/) const
{ return 0.0; }
/*!
* \brief Returns the oil vaporization factor \f$R_v\f$ [m^3/m^3] of the oil phase.
*/

49
opm/material/fluidsystems/blackoilpvt/DryHumidGasPvt.hpp
View File

@@ -110,6 +110,33 @@ public:
setReferenceDensities(regionIdx, rhoRefO, rhoRefG, rhoRefW);
}
enableRwgSalt_ = !eclState.getTableManager().getRwgSaltTables().empty();
if (enableRwgSalt_)
{
const auto& rwgsaltTables = eclState.getTableManager().getRwgSaltTables();
for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
const auto& rwgsaltTable = rwgsaltTables[regionIdx];
const auto& saturatedTable = rwgsaltTable.getSaturatedTable();
assert(saturatedTable.numRows() > 1);
auto& waterVaporizationFac = saturatedWaterVaporizationSaltFactorTable_[regionIdx];
for (unsigned outerIdx = 0; outerIdx < saturatedTable.numRows(); ++ outerIdx) {
const auto& underSaturatedTable = rwgsaltTable.getUnderSaturatedTable(outerIdx);
Scalar pg = saturatedTable.get("PG" , outerIdx);
waterVaporizationFac.appendXPos(pg);
size_t numRows = underSaturatedTable.numRows();
for (size_t innerIdx = 0; innerIdx < numRows; ++ innerIdx) {
Scalar saltConcentration = underSaturatedTable.get("C_SALT" , innerIdx);
Scalar rvwSat= underSaturatedTable.get("RVW" , innerIdx);
waterVaporizationFac.appendSamplePoint(outerIdx, saltConcentration, rvwSat);
}
}
}
}
for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
const auto& pvtgwTable = pvtgwTables[regionIdx];
@@ -478,6 +505,22 @@ public:
return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure, /*extrapolate=*/true);
}
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of the water phase.
*/
template <class Evaluation>
Evaluation saturatedWaterVaporizationFactor(unsigned regionIdx,
const Evaluation& /*temperature*/,
const Evaluation& pressure,
const Evaluation& saltConcentration) const
{
if (enableRwgSalt_)
return saturatedWaterVaporizationSaltFactorTable_[regionIdx].eval(pressure, saltConcentration, /*extrapolate=*/true);
else {
return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure, /*extrapolate=*/true);
}
}
/*!
* \brief Returns the oil vaporization factor \f$R_v\f$ [m^3/m^3] of the oil phase.
*/
@@ -598,6 +641,10 @@ public:
return saturatedWaterVaporizationFactorTable_;
}
const std::vector<TabulatedTwoDFunction>& saturatedWaterVaporizationSaltFactorTable() const {
return saturatedWaterVaporizationSaltFactorTable_;
}
const std::vector<TabulatedOneDFunction>& saturationPressure() const {
return saturationPressure_;
}
@@ -658,8 +705,10 @@ private:
std::vector<TabulatedTwoDFunction> inverseGasBMu_;
std::vector<TabulatedOneDFunction> inverseSaturatedGasBMu_;
std::vector<TabulatedOneDFunction> saturatedWaterVaporizationFactorTable_;
std::vector<TabulatedTwoDFunction> saturatedWaterVaporizationSaltFactorTable_;
std::vector<TabulatedOneDFunction> saturationPressure_;
bool enableRwgSalt_;
Scalar vapPar1_;
};

10
opm/material/fluidsystems/blackoilpvt/GasPvtMultiplexer.hpp
View File

@@ -302,6 +302,16 @@ public:
const Evaluation& temperature,
const Evaluation& pressure) const
{ OPM_GAS_PVT_MULTIPLEXER_CALL(return pvtImpl.saturatedWaterVaporizationFactor(regionIdx, temperature, pressure)); return 0; }
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of water saturated gas.
*/
template <class Evaluation = Scalar>
Evaluation saturatedWaterVaporizationFactor(unsigned regionIdx,
const Evaluation& temperature,
const Evaluation& pressure,
const Evaluation& saltConcentration) const
{ OPM_GAS_PVT_MULTIPLEXER_CALL(return pvtImpl.saturatedWaterVaporizationFactor(regionIdx, temperature, pressure, saltConcentration)); return 0; }
/*!
* \brief Returns the saturation pressure of the gas phase [Pa]

11
opm/material/fluidsystems/blackoilpvt/GasPvtThermal.hpp
View File

@@ -413,6 +413,17 @@ public:
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/) const
{ return 0.0; }
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of water saturated gas.
*/
template <class Evaluation = Scalar>
Evaluation saturatedWaterVaporizationFactor(unsigned /*regionIdx*/,
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/,
const Evaluation& /*saltConcentration*/) const
{ return 0.0; }
/*!

12
opm/material/fluidsystems/blackoilpvt/WetGasPvt.hpp
View File

@@ -535,7 +535,17 @@ public:
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/) const
{ return 0.0; /* this is non-humid gas! */ }
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of water saturated gas.
*/
template <class Evaluation = Scalar>
Evaluation saturatedWaterVaporizationFactor(unsigned /*regionIdx*/,
const Evaluation& /*temperature*/,
const Evaluation& /*pressure*/,
const Evaluation& /*saltConcentration*/) const
{ return 0.0; }
/*!
* \brief Returns the oil vaporization factor \f$R_v\f$ [m^3/m^3] of the gas phase.
*/

61
opm/material/fluidsystems/blackoilpvt/WetHumidGasPvt.hpp
View File

@@ -73,6 +73,7 @@ public:
const std::vector<TabulatedTwoDFunction>& inverseGasBMuRvSat,
const std::vector<TabulatedOneDFunction>& inverseSaturatedGasBMu,
const std::vector<TabulatedOneDFunction>& saturatedWaterVaporizationFactorTable,
const std::vector<TabulatedTwoDFunction>& saturatedWaterVaporizationSaltFactorTable,
const std::vector<TabulatedOneDFunction>& saturatedOilVaporizationFactorTable,
const std::vector<TabulatedOneDFunction>& saturationPressure,
Scalar vapPar1)
@@ -88,6 +89,7 @@ public:
, inverseGasBMuRvSat_(inverseGasBMuRvSat) // Bg^-1*Mug evaluated at saturated oil-gas ratio (Rv) values; pvtgw
, inverseSaturatedGasBMu_(inverseSaturatedGasBMu) //pvtgw
, saturatedWaterVaporizationFactorTable_(saturatedWaterVaporizationFactorTable) //pvtgw
, saturatedWaterVaporizationSaltFactorTable_(saturatedWaterVaporizationSaltFactorTable) //rwgsalt
, saturatedOilVaporizationFactorTable_(saturatedOilVaporizationFactorTable) //pvtg
, saturationPressure_(saturationPressure)
, vapPar1_(vapPar1)
@@ -120,6 +122,34 @@ public:
setReferenceDensities(regionIdx, rhoRefO, rhoRefG, rhoRefW);
}
enableRwgSalt_ = !eclState.getTableManager().getRwgSaltTables().empty();
if (enableRwgSalt_)
{
const auto& rwgsaltTables = eclState.getTableManager().getRwgSaltTables();
for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
const auto& rwgsaltTable = rwgsaltTables[regionIdx];
const auto& saturatedTable = rwgsaltTable.getSaturatedTable();
assert(saturatedTable.numRows() > 1);
auto& waterVaporizationFac = saturatedWaterVaporizationSaltFactorTable_[regionIdx];
for (unsigned outerIdx = 0; outerIdx < saturatedTable.numRows(); ++ outerIdx) {
const auto& underSaturatedTable = rwgsaltTable.getUnderSaturatedTable(outerIdx);
Scalar pg = saturatedTable.get("PG" , outerIdx);
waterVaporizationFac.appendXPos(pg);
size_t numRows = underSaturatedTable.numRows();
for (size_t innerIdx = 0; innerIdx < numRows; ++ innerIdx) {
Scalar saltConcentration = underSaturatedTable.get("C_SALT" , innerIdx);
Scalar rvwSat= underSaturatedTable.get("RVW" , innerIdx);
waterVaporizationFac.appendSamplePoint(outerIdx, saltConcentration, rvwSat);
}
}
}
}
// Table PVTGW
for (unsigned regionIdx = 0; regionIdx < numRegions; ++ regionIdx) {
const auto& pvtgwTable = pvtgwTables[regionIdx];
@@ -293,13 +323,13 @@ public:
pvtgTable.getUnderSaturatedTable(xIdx),
pvtgTable.getUnderSaturatedTable(masterTableIdx));
}
} //end PVTGW
}
vapPar1_ = 0.0;
const auto& oilVap = schedule[0].oilvap();
if (oilVap.getType() == OilVaporizationProperties::OilVaporization::VAPPARS) {
vapPar1_ = oilVap.vap1();
}
initEnd();
}
@@ -422,6 +452,7 @@ public:
gasMuRvwSat_.resize(numRegions, TabulatedTwoDFunction{TabulatedTwoDFunction::InterpolationPolicy::RightExtreme});
gasMuRvSat_.resize(numRegions, TabulatedTwoDFunction{TabulatedTwoDFunction::InterpolationPolicy::RightExtreme});
saturatedWaterVaporizationFactorTable_.resize(numRegions);
saturatedWaterVaporizationSaltFactorTable_.resize(numRegions, TabulatedTwoDFunction{TabulatedTwoDFunction::InterpolationPolicy::RightExtreme});
saturatedOilVaporizationFactorTable_.resize(numRegions);
saturationPressure_.resize(numRegions);
}
@@ -654,6 +685,23 @@ public:
return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure, /*extrapolate=*/true);
}
/*!
* \brief Returns the water vaporization factor \f$R_vw\f$ [m^3/m^3] of the water phase.
*/
template <class Evaluation>
Evaluation saturatedWaterVaporizationFactor(unsigned regionIdx,
const Evaluation& /*temperature*/,
const Evaluation& pressure,
const Evaluation& saltConcentration) const
{
if (enableRwgSalt_)
return saturatedWaterVaporizationSaltFactorTable_[regionIdx].eval(pressure, saltConcentration, /*extrapolate=*/true);
else {
return saturatedWaterVaporizationFactorTable_[regionIdx].eval(pressure, /*extrapolate=*/true);
}
}
template <class Evaluation>
Evaluation saturatedOilVaporizationFactor(unsigned regionIdx,
const Evaluation& /*temperature*/,
@@ -794,12 +842,15 @@ public:
const std::vector<TabulatedOneDFunction>& saturatedWaterVaporizationFactorTable() const {
return saturatedWaterVaporizationFactorTable_;
}
const std::vector<TabulatedTwoDFunction>& saturatedWaterVaporizationSaltFactorTable() const {
return saturatedWaterVaporizationSaltFactorTable_;
}
const std::vector<TabulatedOneDFunction>& saturatedOilVaporizationFactorTable() const {
return saturatedOilVaporizationFactorTable_;
}
const std::vector<TabulatedOneDFunction>& saturationPressure() const {
return saturationPressure_;
}
@@ -862,9 +913,11 @@ private:
std::vector<TabulatedTwoDFunction> inverseGasBMuRvSat_;
std::vector<TabulatedOneDFunction> inverseSaturatedGasBMu_;
std::vector<TabulatedOneDFunction> saturatedWaterVaporizationFactorTable_;
std::vector<TabulatedTwoDFunction> saturatedWaterVaporizationSaltFactorTable_;
std::vector<TabulatedOneDFunction> saturatedOilVaporizationFactorTable_;
std::vector<TabulatedOneDFunction> saturationPressure_;
bool enableRwgSalt_;
Scalar vapPar1_;
};