20#ifndef OPM_BLACKOILINITIALIZATION_HEADER_INCLUDED
21#define OPM_BLACKOILINITIALIZATION_HEADER_INCLUDED
23#include <opm/core/utility/parameters/ParameterGroup.hpp>
24#include <opm/common/ErrorMacros.hpp>
25#include <opm/core/utility/Units.hpp>
33 template <
class Simulator>
37 typedef typename Simulator::State
State;
38 typedef typename Simulator::Grid
Grid;
39 typedef typename Simulator::Fluid
Fluid;
41 virtual void init(
const Opm::parameter::ParameterGroup& param,
44 typename Grid::Vector gravity,
51 template <
class Simulator>
55 typedef typename Simulator::State
State;
56 typedef typename Simulator::Grid
Grid;
57 typedef typename Simulator::Fluid
Fluid;
59 virtual void init(
const Opm::parameter::ParameterGroup& param,
62 typename Grid::Vector gravity,
65 typedef typename Fluid::CompVec CompVec;
66 typedef typename Fluid::PhaseVec PhaseVec;
68 if (param.getDefault(
"heterogenous_initial_mix",
false)) {
69 CompVec init_oil(0.0);
70 init_oil[Fluid::Oil] = 1.0;
71 CompVec init_water(0.0);
72 init_water[Fluid::Water] = 1.0;
73 simstate.cell_z_.resize(grid.numCells());
74 std::fill(simstate.cell_z_.begin(), simstate.cell_z_.begin() + simstate.cell_z_.size()/2, init_oil);
75 std::fill(simstate.cell_z_.begin() + simstate.cell_z_.size()/2, simstate.cell_z_.end(), init_water);
76 OPM_MESSAGE(
"******* Assuming zero capillary pressures *******");
77 PhaseVec init_p(100.0*Opm::unit::barsa);
78 simstate.cell_pressure_.resize(grid.numCells(), init_p);
87 }
else if (param.getDefault(
"unstable_initial_mix",
false)) {
88 CompVec init_oil(0.0);
89 init_oil[Fluid::Oil] = 1.0;
90 init_oil[Fluid::Gas] = 0.0;
91 CompVec init_water(0.0);
92 init_water[Fluid::Water] = 1.0;
93 CompVec init_gas(0.0);
94 init_gas[Fluid::Gas] = 150.0;
95 simstate.cell_z_.resize(grid.numCells());
96 std::fill(simstate.cell_z_.begin(),
97 simstate.cell_z_.begin() + simstate.cell_z_.size()/3,
99 std::fill(simstate.cell_z_.begin() + simstate.cell_z_.size()/3,
100 simstate.cell_z_.begin() + 2*(simstate.cell_z_.size()/3),
102 std::fill(simstate.cell_z_.begin() + 2*(simstate.cell_z_.size()/3),
103 simstate.cell_z_.end(),
105 OPM_MESSAGE(
"******* Assuming zero capillary pressures *******");
106 PhaseVec init_p(100.0*Opm::unit::barsa);
107 simstate.cell_pressure_.resize(grid.numCells(), init_p);
109 if (gravity.two_norm() != 0.0) {
111 typename Fluid::FluidState state = fluid.computeState(simstate.cell_pressure_[0], simstate.cell_z_[0]);
112 simstate.cell_z_[0] *= 1.0/state.total_phase_volume_density_;
113 for (
int cell = 1; cell < grid.numCells(); ++cell) {
114 double fluid_vol_dens;
117 double rho = 0.5*((simstate.cell_z_[cell]+simstate.cell_z_[cell-1])*fluid.surfaceDensities());
118 double press = rho*((grid.cellCentroid(cell) - grid.cellCentroid(cell-1))*gravity) + simstate.cell_pressure_[cell-1][0];
119 simstate.cell_pressure_[cell] = PhaseVec(press);
120 typename Fluid::FluidState state = fluid.computeState(simstate.cell_pressure_[cell], simstate.cell_z_[cell]);
121 fluid_vol_dens = state.total_phase_volume_density_;
122 simstate.cell_z_[cell] *= 1.0/fluid_vol_dens;
124 }
while (std::fabs((fluid_vol_dens-1.0)) > 1.0e-8 && cnt < 10);
128 std::cout <<
"---- Exit - BlackoilSimulator.hpp: No gravity, no fun ... ----" << std::endl;
131 }
else if (param.getDefault(
"CO2-injection",
false)) {
132 CompVec init_water(0.0);
135 init_water[Fluid::Oil] = 1.0;
136 simstate.cell_z_.resize(grid.numCells());
137 std::fill(simstate.cell_z_.begin(),simstate.cell_z_.end(),init_water);
139 double datum_pressure_barsa = param.getDefault<
double>(
"datum_pressure", 200.0);
140 double datum_pressure = Opm::unit::convert::from(datum_pressure_barsa, Opm::unit::barsa);
141 PhaseVec init_p(datum_pressure);
142 simstate.cell_pressure_.resize(grid.numCells(), init_p);
145 double zMin = grid.cellCentroid(0)[2];
146 for (
int cell = 1; cell < grid.numCells(); ++cell) {
147 if (grid.cellCentroid(cell)[2] < zMin)
148 zMin = grid.cellCentroid(cell)[2];
151 typename Fluid::FluidState state = fluid.computeState(init_p, init_water);
152 simstate.cell_z_[0] *= 1.0/state.total_phase_volume_density_;
153 double density = (init_water*fluid.surfaceDensities())/state.total_phase_volume_density_;
155 for (
int cell = 0; cell < grid.numCells(); ++cell) {
156 double pressure(datum_pressure + (grid.cellCentroid(cell)[2] - zMin)*gravity[2]*density);
157 simstate.cell_pressure_[cell] = PhaseVec(pressure);
158 state = fluid.computeState(simstate.cell_pressure_[cell], simstate.cell_z_[cell]);
159 simstate.cell_z_[cell] *= 1.0/state.total_phase_volume_density_;
163 double initial_mixture_gas = param.getDefault(
"initial_mixture_gas", 0.0);
164 double initial_mixture_oil = param.getDefault(
"initial_mixture_oil", 1.0);
165 double initial_mixture_water = param.getDefault(
"initial_mixture_water", 0.0);
166 init_z[Fluid::Water] = initial_mixture_water;
167 init_z[Fluid::Gas] = initial_mixture_gas;
168 init_z[Fluid::Oil] = initial_mixture_oil;
170 simstate.cell_z_.resize(grid.numCells(), init_z);
171 OPM_MESSAGE(
"******* Assuming zero capillary pressures *******");
172 PhaseVec init_p(param.getDefault(
"initial_pressure", 100.0*Opm::unit::barsa));
173 simstate.cell_pressure_.resize(grid.numCells(), init_p);
174 if (gravity.two_norm() != 0.0) {
175 double ref_gravpot = grid.cellCentroid(0)*gravity;
176 double rho = init_z*fluid.surfaceDensities();
177 for (
int cell = 1; cell < grid.numCells(); ++cell) {
178 double press = rho*(grid.cellCentroid(cell)*gravity - ref_gravpot) + simstate.cell_pressure_[0][0];
179 simstate.cell_pressure_[cell] = PhaseVec(press);
192 template <
class Simulator>
196 typedef typename Simulator::State
State;
197 typedef typename Simulator::Grid
Grid;
198 typedef typename Simulator::Fluid
Fluid;
200 virtual void init(
const Opm::parameter::ParameterGroup& param,
203 typename Grid::Vector gravity,
206 typedef typename Fluid::CompVec CompVec;
208 double zeroDepth = param.getDefault(
"zero_depth", 2743.2);
209 int nx = param.getDefault<
int>(
"nx", 24);
210 int ny = param.getDefault<
int>(
"ny", 25);
211 int nz = param.getDefault<
int>(
"nz", 15);
214 double datum_pressure_barsa = param.getDefault<
double>(
"datum_pressure", 248.22);
215 double datum_pressure = Opm::unit::convert::from(datum_pressure_barsa, Opm::unit::barsa);
216 double wo_contact_depth = param.getDefault<
double>(
"wo_contact_depth", 3032.76) - zeroDepth;
217 double go_contact_depth = param.getDefault<
double>(
"go_contact_depth", 2682.24) - zeroDepth;
219 double connate_water_saturation = param.getDefault<
double>(
"connate_water_saturation", 0.151090);
220 double residual_oil_saturation = param.getDefault<
double>(
"residual_oil_saturation", 0.118510);
222 double initial_mixture_gas = param.getDefault(
"initial_mixture_gas", 247.43);
223 double initial_mixture_oil = param.getDefault(
"initial_mixture_oil", 1.0);
226 CompVec oil_sample(0.0);
227 oil_sample[Fluid::Oil] = initial_mixture_oil;
228 oil_sample[Fluid::Gas] = initial_mixture_gas;
229 CompVec water_sample(0.0);
230 water_sample[Fluid::Water] = 1.0;
232 simstate.cell_z_.resize(grid.numCells());
233 simstate.cell_pressure_.resize(grid.numCells());
237 simstate.cell_pressure_[0] = datum_pressure;
238 typename Fluid::FluidState state = fluid.computeState(simstate.cell_pressure_[0],oil_sample);
239 simstate.cell_z_[0] = oil_sample;
240 simstate.cell_z_[0] *= (1.0-connate_water_saturation)/state.total_phase_volume_density_;
241 state = fluid.computeState(simstate.cell_pressure_[0],water_sample);
242 simstate.cell_z_[0][Fluid::Water] = water_sample[Fluid::Water];
243 simstate.cell_z_[0][Fluid::Water] *= connate_water_saturation/state.total_phase_volume_density_;
245 for (
int i=0; i<nx; ++i) {
247 for (
int k=0; k<nz; ++k) {
251 kk, kk-nz, wo_contact_depth, go_contact_depth, connate_water_saturation,
252 residual_oil_saturation, simstate);
255 kk, kk-1, wo_contact_depth, go_contact_depth, connate_water_saturation,
256 residual_oil_saturation, simstate);
259 for (
int j=1; j<ny; ++j) {
260 int jj = j*nx*nz + kk;
261 simstate.cell_z_[jj] = simstate.cell_z_[kk];
262 simstate.cell_pressure_[jj] = simstate.cell_pressure_[kk];
272 typename Grid::Vector gravity,
275 double wo_contact_depth,
276 double go_contact_depth,
277 double connate_water_saturation,
278 double residual_oil_saturation,
281 typedef typename Fluid::PhaseVec PhaseVec;
283 const int maxCnt = 30;
284 const double eps = 1.0e-8;
285 simstate.cell_z_[iCell] = simstate.cell_z_[iRef];
287 bool below_wo_contact =
false;
288 if (grid.cellCentroid(iCell)[2] > wo_contact_depth)
289 below_wo_contact =
true;
291 double gZ = (grid.cellCentroid(iCell) - grid.cellCentroid(iRef))*gravity;
292 double fluid_vol_dens;
295 double rho = 0.5*(simstate.cell_z_[iCell]*fluid.surfaceDensities()
296 + simstate.cell_z_[iRef]*fluid.surfaceDensities());
297 double press = rho*gZ + simstate.cell_pressure_[iRef][0];
298 simstate.cell_pressure_[iCell] = PhaseVec(press);
299 typename Fluid::FluidState state = fluid.computeState(simstate.cell_pressure_[iCell], simstate.cell_z_[iCell]);
300 fluid_vol_dens = state.total_phase_volume_density_;
301 double oil_vol_dens = state.phase_volume_density_[Fluid::Liquid]
302 + state.phase_volume_density_[Fluid::Vapour];
303 double wat_vol_dens = state.phase_volume_density_[Fluid::Aqua];
304 if (below_wo_contact) {
305 simstate.cell_z_[iCell][Fluid::Oil] *= residual_oil_saturation/oil_vol_dens;
306 simstate.cell_z_[iCell][Fluid::Gas] *= residual_oil_saturation/oil_vol_dens;
307 simstate.cell_z_[iCell][Fluid::Water] *= (1.0-residual_oil_saturation)/wat_vol_dens;
309 simstate.cell_z_[iCell][Fluid::Oil] *= (1.0-connate_water_saturation)/oil_vol_dens;
310 simstate.cell_z_[iCell][Fluid::Gas] *= (1.0-connate_water_saturation)/oil_vol_dens;
311 simstate.cell_z_[iCell][Fluid::Water] *= connate_water_saturation/wat_vol_dens;
314 }
while (std::fabs(fluid_vol_dens-1.0) > eps && cnt < maxCnt);
317 std::cout <<
"z_cell_[" << iCell <<
"]: " << simstate.cell_z_[iCell]
318 <<
" pressure: " << simstate.cell_pressure_[iCell][Fluid::Liquid]
320 <<
" eps: " << std::fabs(fluid_vol_dens-1.0) << std::endl;
323 return (cnt < maxCnt);
Initialize basic cases.
Definition: BlackoilInitialization.hpp:53
Simulator::State State
Definition: BlackoilInitialization.hpp:55
virtual void init(const Opm::parameter::ParameterGroup ¶m, const Grid &grid, const Fluid &fluid, typename Grid::Vector gravity, State &simstate)
Definition: BlackoilInitialization.hpp:59
Simulator::Fluid Fluid
Definition: BlackoilInitialization.hpp:57
Simulator::Grid Grid
Definition: BlackoilInitialization.hpp:56
Base class for initialization codes.
Definition: BlackoilInitialization.hpp:35
Simulator::State State
Definition: BlackoilInitialization.hpp:37
virtual void init(const Opm::parameter::ParameterGroup ¶m, const Grid &grid, const Fluid &fluid, typename Grid::Vector gravity, State &simstate)=0
Simulator::Fluid Fluid
Definition: BlackoilInitialization.hpp:39
Simulator::Grid Grid
Definition: BlackoilInitialization.hpp:38
Initialize SPE9 type case.
Definition: BlackoilInitialization.hpp:194
virtual void init(const Opm::parameter::ParameterGroup ¶m, const Grid &grid, const Fluid &fluid, typename Grid::Vector gravity, State &simstate)
Definition: BlackoilInitialization.hpp:200
Simulator::State State
Definition: BlackoilInitialization.hpp:196
bool computeCellState(const Grid &grid, const Fluid &fluid, typename Grid::Vector gravity, int iCell, int iRef, double wo_contact_depth, double go_contact_depth, double connate_water_saturation, double residual_oil_saturation, State &simstate)
Definition: BlackoilInitialization.hpp:270
Simulator::Grid Grid
Definition: BlackoilInitialization.hpp:197
Simulator::Fluid Fluid
Definition: BlackoilInitialization.hpp:198
Definition: BlackoilFluid.hpp:32