OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

added salt concentration dependency of water-gas ratio

Browse Source
This commit is contained in:
Paul Egberts 2022-10-20 11:22:18 +02:00
parent 6d747b6620
commit d420a78585
1 changed files with 14 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
opm/models/blackoil/blackoilprimaryvariables.hh
View File

@ -450,12 +450,14 @@ public:
return false; return false;
} }
Scalar Sw = 0.0; Scalar Sw = 0.0;
Scalar saltConcentration = 0.0;
if (waterEnabled && primaryVarsMeaning() != Rvw_po_Sg && primaryVarsMeaning() != Rvw_pg_Rv ) if (waterEnabled && primaryVarsMeaning() != Rvw_po_Sg && primaryVarsMeaning() != Rvw_pg_Rv )
Sw = (*this)[Indices::waterSaturationIdx]; Sw = (*this)[Indices::waterSaturationIdx];
if constexpr (enableSaltPrecipitation) { if constexpr (enableSaltPrecipitation) {
Scalar saltSolubility = BrineModule::saltSol(pvtRegionIndex()); Scalar saltSolubility = BrineModule::saltSol(pvtRegionIndex());
if (primaryVarsMeaningBrine() == Sp) { if (primaryVarsMeaningBrine() == Sp) {
saltConcentration = saltSolubility;
Scalar saltSat = (*this)[saltConcentrationIdx]; Scalar saltSat = (*this)[saltConcentrationIdx];
if (saltSat < -eps){ //precipitated salt dissappears if (saltSat < -eps){ //precipitated salt dissappears
setPrimaryVarsMeaningBrine(Cs); setPrimaryVarsMeaningBrine(Cs);
@ -463,8 +465,8 @@ public:
} }
} }
else if (primaryVarsMeaningBrine() == Cs) { else if (primaryVarsMeaningBrine() == Cs) {
Scalar saltConc = (*this)[saltConcentrationIdx]; saltConcentration = (*this)[saltConcentrationIdx];
if (saltConc > saltSolubility + eps){ //salt concentration exceeds solubility limit if (saltConcentration > saltSolubility + eps){ //salt concentration exceeds solubility limit
setPrimaryVarsMeaningBrine(Sp); setPrimaryVarsMeaningBrine(Sp);
(*this)[saltConcentrationIdx] = 0.0; (*this)[saltConcentrationIdx] = 0.0;
} }
@ -519,7 +521,8 @@ public:
} }
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_, Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T, T,
pg); pg,
saltConcentration);
setPrimaryVarsMeaning(Rvw_po_Sg); setPrimaryVarsMeaning(Rvw_po_Sg);
(*this)[Indices::waterSaturationIdx] = RvwSat; //primary variable becomes Rvw (*this)[Indices::waterSaturationIdx] = RvwSat; //primary variable becomes Rvw
@ -535,7 +538,8 @@ public:
Scalar pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]); Scalar pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]);
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_, Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T, T,
pg); pg,
saltConcentration);
setPrimaryVarsMeaning(Rvw_pg_Rv); setPrimaryVarsMeaning(Rvw_pg_Rv);
(*this)[Indices::pressureSwitchIdx] = pg; (*this)[Indices::pressureSwitchIdx] = pg;
@ -683,7 +687,8 @@ public:
Scalar pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]); Scalar pg = po + (pC[gasPhaseIdx] - pC[oilPhaseIdx]);
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_, Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T, T,
pg); pg,
saltConcentration);
Scalar Rvw = (*this)[Indices::waterSaturationIdx]; Scalar Rvw = (*this)[Indices::waterSaturationIdx];
if (Rvw > RvwSat*(1.0 + eps)) { if (Rvw > RvwSat*(1.0 + eps)) {
// water phase appears // water phase appears
@ -722,7 +727,8 @@ public:
Scalar T = asImp_().temperature_(); Scalar T = asImp_().temperature_();
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_, Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T, T,
pg); pg,
saltConcentration);
Scalar Rvw = (*this)[Indices::waterSaturationIdx]; Scalar Rvw = (*this)[Indices::waterSaturationIdx];
if (Rvw > RvwSat*(1.0 + eps)) { if (Rvw > RvwSat*(1.0 + eps)) {
// water phase appears // water phase appears
@ -801,7 +807,8 @@ public:
if(Sw < -eps && FluidSystem::enableVaporizedWater()) { if(Sw < -eps && FluidSystem::enableVaporizedWater()) {
Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_, Scalar RvwSat = FluidSystem::gasPvt().saturatedWaterVaporizationFactor(pvtRegionIdx_,
T, T,
pg); pg,
saltConcentration);
setPrimaryVarsMeaning(Rvw_pg_Rv); setPrimaryVarsMeaning(Rvw_pg_Rv);
(*this)[Indices::waterSaturationIdx] = RvwSat; //primary variable becomes Rvw (*this)[Indices::waterSaturationIdx] = RvwSat; //primary variable becomes Rvw