/// Constructor SinglePvtDeadSpline::SinglePvtDeadSpline(const table_t& pvd_table, const int samples) { const int region_number = 0; if (pvd_table.size() != 1) { THROW("More than one PVT-region"); } // Copy data const int sz = pvd_table[region_number][0].size(); std::vector<double> press(sz); std::vector<double> B_inv(sz); std::vector<double> visc(sz); for (int i = 0; i < sz; ++i) { press[i] = pvd_table[region_number][0][i]; B_inv[i] = 1.0 / pvd_table[region_number][1][i]; visc[i] = pvd_table[region_number][2][i]; } buildUniformMonotoneTable(press, B_inv, samples, one_over_B_); buildUniformMonotoneTable(press, visc, samples, viscosity_); // Dumping the created tables. // static int count = 0; // std::ofstream os((std::string("dump-") + boost::lexical_cast<std::string>(count++)).c_str()); // os.precision(15); // os << "1/B\n\n" << one_over_B_ // << "\n\nvisc\n\n" << viscosity_ << std::endl; }
void SatFuncBase<UniformTableLinear<double> >::initializeTableType(UniformTableLinear<double> & table, const std::vector<double>& arg, const std::vector<double>& value, const int samples) { buildUniformMonotoneTable(arg, value, samples, table); }
void SatFuncGwsegUniform::init(const EclipseGridParser& deck, const int table_num, const PhaseUsage phase_usg, const int samples) { phase_usage = phase_usg; double swco = 0.0; double swmax = 1.0; if (phase_usage.phase_used[Aqua]) { const SWOF::table_t& swof_table = deck.getSWOF().swof_; const std::vector<double>& sw = swof_table[table_num][0]; const std::vector<double>& krw = swof_table[table_num][1]; const std::vector<double>& krow = swof_table[table_num][2]; const std::vector<double>& pcow = swof_table[table_num][3]; buildUniformMonotoneTable(sw, krw, samples, krw_); buildUniformMonotoneTable(sw, krow, samples, krow_); buildUniformMonotoneTable(sw, pcow, samples, pcow_); krocw_ = krow[0]; // At connate water -> ecl. SWOF swco = sw[0]; smin_[phase_usage.phase_pos[Aqua]] = sw[0]; swmax = sw.back(); smax_[phase_usage.phase_pos[Aqua]] = sw.back(); } if (phase_usage.phase_used[Vapour]) { const SGOF::table_t& sgof_table = deck.getSGOF().sgof_; const std::vector<double>& sg = sgof_table[table_num][0]; const std::vector<double>& krg = sgof_table[table_num][1]; const std::vector<double>& krog = sgof_table[table_num][2]; const std::vector<double>& pcog = sgof_table[table_num][3]; buildUniformMonotoneTable(sg, krg, samples, krg_); buildUniformMonotoneTable(sg, krog, samples, krog_); buildUniformMonotoneTable(sg, pcog, samples, pcog_); smin_[phase_usage.phase_pos[Vapour]] = sg[0]; if (std::fabs(sg.back() + swco - 1.0) > 1e-3) { THROW("Gas maximum saturation in SGOF table = " << sg.back() << ", should equal (1.0 - connate water sat) = " << (1.0 - swco)); } smax_[phase_usage.phase_pos[Vapour]] = sg.back(); } // These only consider water min/max sats. Consider gas sats? smin_[phase_usage.phase_pos[Liquid]] = 1.0 - swmax; smax_[phase_usage.phase_pos[Liquid]] = 1.0 - swco; }
void SatFuncSimpleUniform::init(const EclipseGridParser& deck, const int table_num, const PhaseUsage phase_usg, const int samples) { phase_usage = phase_usg; double swco = 0.0; double swmax = 1.0; if (phase_usage.phase_used[Aqua]) { const SWOF::table_t& swof_table = deck.getSWOF().swof_; const std::vector<double>& sw = swof_table[table_num][0]; const std::vector<double>& krw = swof_table[table_num][1]; const std::vector<double>& krow = swof_table[table_num][2]; const std::vector<double>& pcow = swof_table[table_num][3]; if (krw.front() != 0.0 || krow.back() != 0.0) { THROW("Error SWOF data - non-zero krw(swco) and/or krow(1-sor)"); } buildUniformMonotoneTable(sw, krw, samples, krw_); buildUniformMonotoneTable(sw, krow, samples, krow_); buildUniformMonotoneTable(sw, pcow, samples, pcow_); krocw_ = krow[0]; // At connate water -> ecl. SWOF swco = sw[0]; smin_[phase_usage.phase_pos[Aqua]] = sw[0]; swmax = sw.back(); smax_[phase_usage.phase_pos[Aqua]] = sw.back(); krwmax_ = krw.back(); kromax_ = krow.front(); swcr_ = swmax; sowcr_ = 1.0 - swco; krwr_ = krw.back(); krorw_ = krow.front(); for (std::vector<double>::size_type i=1; i<sw.size(); ++i) { if (krw[i]> 0.0) { swcr_ = sw[i-1]; krorw_ = krow[i-1]; break; } } for (std::vector<double>::size_type i=sw.size()-1; i>=1; --i) { if (krow[i-1]> 0.0) { sowcr_ = 1.0 - sw[i]; krwr_ = krw[i]; break; } } } if (phase_usage.phase_used[Vapour]) { const SGOF::table_t& sgof_table = deck.getSGOF().sgof_; const std::vector<double>& sg = sgof_table[table_num][0]; const std::vector<double>& krg = sgof_table[table_num][1]; const std::vector<double>& krog = sgof_table[table_num][2]; const std::vector<double>& pcog = sgof_table[table_num][3]; buildUniformMonotoneTable(sg, krg, samples, krg_); buildUniformMonotoneTable(sg, krog, samples, krog_); buildUniformMonotoneTable(sg, pcog, samples, pcog_); smin_[phase_usage.phase_pos[Vapour]] = sg[0]; if (std::fabs(sg.back() + swco - 1.0) > 1e-3) { THROW("Gas maximum saturation in SGOF table = " << sg.back() << ", should equal (1.0 - connate water sat) = " << (1.0 - swco)); } smax_[phase_usage.phase_pos[Vapour]] = sg.back(); } // These only consider water min/max sats. Consider gas sats? smin_[phase_usage.phase_pos[Liquid]] = 1.0 - swmax; smax_[phase_usage.phase_pos[Liquid]] = 1.0 - swco; }
void SatFuncGwsegUniform::init(const EclipseGridParser& deck, const int table_num, const PhaseUsage phase_usg, const int samples) { phase_usage = phase_usg; double swco = 0.0; double swmax = 1.0; if (phase_usage.phase_used[Aqua]) { const SWOF::table_t& swof_table = deck.getSWOF().swof_; const std::vector<double>& sw = swof_table[table_num][0]; const std::vector<double>& krw = swof_table[table_num][1]; const std::vector<double>& krow = swof_table[table_num][2]; const std::vector<double>& pcow = swof_table[table_num][3]; // Extend the tables with constant values such that the // derivatives at the endpoints are zero int n = sw.size(); std::vector<double> sw_ex(n+2); std::vector<double> krw_ex(n+2); std::vector<double> krow_ex(n+2); std::vector<double> pcow_ex(n+2); SatFuncGwsegUniform::ExtendTable(sw,sw_ex,1); SatFuncGwsegUniform::ExtendTable(krw,krw_ex,0); SatFuncGwsegUniform::ExtendTable(krow,krow_ex,0); SatFuncGwsegUniform::ExtendTable(pcow,pcow_ex,0); buildUniformMonotoneTable(sw_ex, krw_ex, samples, krw_); buildUniformMonotoneTable(sw_ex, krow_ex, samples, krow_); buildUniformMonotoneTable(sw_ex, pcow_ex, samples, pcow_); krocw_ = krow[0]; // At connate water -> ecl. SWOF swco = sw[0]; smin_[phase_usage.phase_pos[Aqua]] = sw[0]; swmax = sw.back(); smax_[phase_usage.phase_pos[Aqua]] = sw.back(); } if (phase_usage.phase_used[Vapour]) { const SGOF::table_t& sgof_table = deck.getSGOF().sgof_; const std::vector<double>& sg = sgof_table[table_num][0]; const std::vector<double>& krg = sgof_table[table_num][1]; const std::vector<double>& krog = sgof_table[table_num][2]; const std::vector<double>& pcog = sgof_table[table_num][3]; // Extend the tables with constant values such that the // derivatives at the endpoints are zero int n = sg.size(); std::vector<double> sg_ex(n+2); std::vector<double> krg_ex(n+2); std::vector<double> krog_ex(n+2); std::vector<double> pcog_ex(n+2); SatFuncGwsegUniform::ExtendTable(sg,sg_ex,1); SatFuncGwsegUniform::ExtendTable(krg,krg_ex,0); SatFuncGwsegUniform::ExtendTable(krog,krog_ex,0); SatFuncGwsegUniform::ExtendTable(pcog,pcog_ex,0); buildUniformMonotoneTable(sg_ex, krg_ex, samples, krg_); buildUniformMonotoneTable(sg_ex, krog_ex, samples, krog_); buildUniformMonotoneTable(sg_ex, pcog_ex, samples, pcog_); smin_[phase_usage.phase_pos[Vapour]] = sg[0]; if (std::fabs(sg.back() + swco - 1.0) > 1e-3) { OPM_THROW(std::runtime_error, "Gas maximum saturation in SGOF table = " << sg.back() << ", should equal (1.0 - connate water sat) = " << (1.0 - swco)); } smax_[phase_usage.phase_pos[Vapour]] = sg.back(); } // These only consider water min/max sats. Consider gas sats? smin_[phase_usage.phase_pos[Liquid]] = 1.0 - swmax; smax_[phase_usage.phase_pos[Liquid]] = 1.0 - swco; }