mirror of
https://github.com/OPM/opm-simulators.git
synced 2024-12-22 15:33:29 -06:00
749d0e838d
the largest change is that all classes below opm/core/props/pvt take the PVT region index as an argument, the higher-level ones (i.e., BlackoilProps*) take cell indices.
116 lines
4.5 KiB
C++
116 lines
4.5 KiB
C++
/*
|
|
Copyright 2012 SINTEF ICT, Applied Mathematics.
|
|
|
|
This file is part of the Open Porous Media project (OPM).
|
|
|
|
OPM is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OPM is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OPM. If not, see <http://www.gnu.org/licenses/>.
|
|
*/
|
|
|
|
|
|
#include "config.h"
|
|
#include <opm/core/props/pvt/PvtPropertiesIncompFromDeck.hpp>
|
|
#include <opm/core/props/phaseUsageFromDeck.hpp>
|
|
#include <opm/core/utility/Units.hpp>
|
|
#include <opm/core/utility/ErrorMacros.hpp>
|
|
#include <opm/core/props/BlackoilPhases.hpp>
|
|
|
|
|
|
namespace Opm
|
|
{
|
|
|
|
PvtPropertiesIncompFromDeck::PvtPropertiesIncompFromDeck()
|
|
{
|
|
}
|
|
|
|
void PvtPropertiesIncompFromDeck::init(Opm::DeckConstPtr deck)
|
|
{
|
|
// So far, this class only supports a single PVT region. TODO?
|
|
int region_number = 0;
|
|
|
|
PhaseUsage phase_usage = phaseUsageFromDeck(deck);
|
|
if (phase_usage.phase_used[PhaseUsage::Vapour] ||
|
|
!phase_usage.phase_used[PhaseUsage::Aqua] ||
|
|
!phase_usage.phase_used[PhaseUsage::Liquid]) {
|
|
OPM_THROW(std::runtime_error, "PvtPropertiesIncompFromDeck::init() -- must have gas and oil phases (only) in deck input.\n");
|
|
}
|
|
|
|
// Surface densities. Accounting for different orders in eclipse and our code.
|
|
if (deck->hasKeyword("DENSITY")) {
|
|
Opm::DeckRecordConstPtr densityRecord = deck->getKeyword("DENSITY")->getRecord(region_number);
|
|
surface_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] = densityRecord->getItem("OIL")->getSIDouble(0);
|
|
surface_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] = densityRecord->getItem("WATER")->getSIDouble(0);
|
|
} else {
|
|
OPM_THROW(std::runtime_error, "Input is missing DENSITY\n");
|
|
}
|
|
|
|
// Make reservoir densities the same as surface densities initially.
|
|
// We will modify them with formation volume factors if found.
|
|
reservoir_density_ = surface_density_;
|
|
|
|
// Water viscosity.
|
|
if (deck->hasKeyword("PVTW")) {
|
|
Opm::DeckRecordConstPtr pvtwRecord = deck->getKeyword("PVTW")->getRecord(region_number);
|
|
if (pvtwRecord->getItem("WATER_COMPRESSIBILITY")->getSIDouble(0) != 0.0 ||
|
|
pvtwRecord->getItem("WATER_VISCOSIBILITY")->getSIDouble(0) != 0.0) {
|
|
OPM_MESSAGE("Compressibility effects in PVTW are ignored.");
|
|
}
|
|
reservoir_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] /= pvtwRecord->getItem("WATER_VOL_FACTOR")->getSIDouble(0);
|
|
viscosity_[phase_usage.phase_pos[PhaseUsage::Aqua]] = pvtwRecord->getItem("WATER_VISCOSITY")->getSIDouble(0);
|
|
} else {
|
|
// Eclipse 100 default.
|
|
// viscosity_[phase_usage.phase_pos[PhaseUsage::Aqua]] = 0.5*Opm::prefix::centi*Opm::unit::Poise;
|
|
OPM_THROW(std::runtime_error, "Input is missing PVTW\n");
|
|
}
|
|
|
|
// Oil viscosity.
|
|
if (deck->hasKeyword("PVCDO")) {
|
|
Opm::DeckRecordConstPtr pvcdoRecord = deck->getKeyword("PVCDO")->getRecord(region_number);
|
|
|
|
if (pvcdoRecord->getItem("OIL_COMPRESSIBILITY")->getSIDouble(0) != 0.0 ||
|
|
pvcdoRecord->getItem("OIL_VISCOSIBILITY")->getSIDouble(0) != 0.0) {
|
|
OPM_MESSAGE("Compressibility effects in PVCDO are ignored.");
|
|
}
|
|
reservoir_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] /= pvcdoRecord->getItem("OIL_VOL_FACTOR")->getSIDouble(0);
|
|
viscosity_[phase_usage.phase_pos[PhaseUsage::Liquid]] = pvcdoRecord->getItem("OIL_VISCOSITY")->getSIDouble(0);
|
|
} else {
|
|
OPM_THROW(std::runtime_error, "Input is missing PVCDO\n");
|
|
}
|
|
}
|
|
|
|
const double* PvtPropertiesIncompFromDeck::surfaceDensities() const
|
|
{
|
|
return surface_density_.data();
|
|
}
|
|
|
|
|
|
const double* PvtPropertiesIncompFromDeck::reservoirDensities() const
|
|
{
|
|
return reservoir_density_.data();
|
|
}
|
|
|
|
|
|
const double* PvtPropertiesIncompFromDeck::viscosity() const
|
|
{
|
|
return viscosity_.data();
|
|
}
|
|
|
|
|
|
int PvtPropertiesIncompFromDeck::numPhases() const
|
|
{
|
|
return 2;
|
|
}
|
|
|
|
|
|
} // namespace Opm
|