EclEpsConfig.hpp

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
/*
  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 2 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/>.
 
  Consult the COPYING file in the top-level source directory of this
  module for the precise wording of the license and the list of
  copyright holders.
*/
#ifndef OPM_ECL_EPS_CONFIG_HPP
#define OPM_ECL_EPS_CONFIG_HPP
 
#if HAVE_ECL_INPUT
#include <opm/input/eclipse/EclipseState/EclipseState.hpp>
#include <stdexcept>
#endif
 
#include <cassert>
#include <string>
 
namespace Opm {
enum EclTwoPhaseSystemType {
    EclGasOilSystem,
    EclOilWaterSystem,
    EclGasWaterSystem
};
 
class EclEpsConfig
{
public:
    EclEpsConfig()
    {
        // by default, do not scale anything
        enableSatScaling_ = false;
        enablePcScaling_ = false;
        enableLeverettScaling_ = false;
        enableKrwScaling_ = false;
        enableKrnScaling_ = false;
        enableThreePointKrSatScaling_ = false;
    }
 
    void setEnableSatScaling(bool yesno)
    { enableSatScaling_ = yesno; }
 
    bool enableSatScaling() const
    { return enableSatScaling_; }
 
    void setEnableThreePointKrSatScaling(bool yesno)
    { enableThreePointKrSatScaling_ = yesno; }
 
    bool enableThreePointKrSatScaling() const
    { return enableThreePointKrSatScaling_; }
 
    void setEnableKrwScaling(bool yesno)
    { enableKrwScaling_ = yesno; }
 
    bool enableKrwScaling() const
    { return enableKrwScaling_; }
 
    void setEnableThreePointKrwScaling(const bool enable)
    { this->enableThreePointKrwScaling_ = enable; }
 
    bool enableThreePointKrwScaling() const
    { return this->enableThreePointKrwScaling_; }
 
    void setEnableThreePointKrnScaling(const bool enable)
    { this->enableThreePointKrnScaling_ = enable; }
 
    bool enableThreePointKrnScaling() const
    { return this->enableThreePointKrnScaling_; }
 
    void setEnableKrnScaling(bool yesno)
    { enableKrnScaling_ = yesno; }
 
    bool enableKrnScaling() const
    { return enableKrnScaling_; }
 
    void setEnablePcScaling(bool yesno)
    { enablePcScaling_ = yesno; }
 
    bool enablePcScaling() const
    { return enablePcScaling_; }
 
    void setEnableLeverettScaling(bool yesno)
    { enableLeverettScaling_ = yesno; }
 
    bool enableLeverettScaling() const
    { return enableLeverettScaling_; }
 
#if HAVE_ECL_INPUT
    void initFromState(const EclipseState& eclState,
                       EclTwoPhaseSystemType twoPhaseSystemType,
                       const std::string& prefix = "",
                       const std::string& suffix = "")
    {
        const auto& endscale = eclState.runspec().endpointScaling();
        // find out if endpoint scaling is used in the first place
        if (!endscale) {
            // it is not used, i.e., just set all enable$Foo attributes to 0 and be done
            // with it.
            enableSatScaling_ = false;
            enableThreePointKrSatScaling_ = false;
            enablePcScaling_ = false;
            enableLeverettScaling_ = false;
            enableKrwScaling_ = false;
            enableKrnScaling_ = false;
            return;
        }
 
        // endpoint scaling is used, i.e., at least saturation scaling needs to be enabled
        enableSatScaling_ = true;
        enableThreePointKrSatScaling_ = endscale.threepoint();
 
        if (eclState.getTableManager().useJFunc()) {
            const auto flag = eclState.getTableManager().getJFunc().flag();
 
            enableLeverettScaling_ = (flag == JFunc::Flag::BOTH)
                || ((twoPhaseSystemType == EclOilWaterSystem) &&
                    (flag == JFunc::Flag::WATER))
                || ((twoPhaseSystemType == EclGasOilSystem) &&
                    (flag == JFunc::Flag::GAS));
        }
 
        const auto& fp = eclState.fieldProps();
        auto hasKR = [&fp, &prefix, &suffix](const std::string& scaling)
        {
            return fp.has_double(prefix + "KR" + scaling + suffix);
        };
        auto hasPC = [&fp, &prefix](const std::string& scaling)
        {
            return fp.has_double(prefix + "PC" + scaling);
        };
 
        // check if we are supposed to scale the Y axis of the capillary pressure
        if (twoPhaseSystemType == EclOilWaterSystem) {
            this->setEnableThreePointKrwScaling(hasKR("WR"));
            this->setEnableThreePointKrnScaling(hasKR("ORW"));
 
            this->enableKrnScaling_ = hasKR("O") || this->enableThreePointKrnScaling();
            this->enableKrwScaling_ = hasKR("W") || this->enableThreePointKrwScaling();
            this->enablePcScaling_  = hasPC("W") || fp.has_double("SWATINIT");
        }
        else if (twoPhaseSystemType == EclGasOilSystem) {
            this->setEnableThreePointKrwScaling(hasKR("ORG"));
            this->setEnableThreePointKrnScaling(hasKR("GR"));
 
            this->enableKrnScaling_ = hasKR("G") || this->enableThreePointKrnScaling();
            this->enableKrwScaling_ = hasKR("O") || this->enableThreePointKrwScaling();
            this->enablePcScaling_  = hasPC("G");
        }
        else {
            assert(twoPhaseSystemType == EclGasWaterSystem);
            //TODO enable endpoint scaling for gaswater system
        }
 
        if (enablePcScaling_ && enableLeverettScaling_) {
            throw std::runtime_error {
                "Capillary pressure scaling and the Leverett scaling function "
                "are mutually exclusive. The deck contains the PCW/PCG property "
                "and the JFUNC keyword applies to the water phase."
            };
        }
    }
#endif
 
private:
    // enable scaling of the input saturations (i.e., rescale the x-Axis)
    bool enableSatScaling_;
 
    // use three (instead of two) points to scale the saturations for the relative
    // permeabilities.
    //
    // This means that two piecewise linear functions are used for saturation scaling
    // instead of a single linear one
    bool enableThreePointKrSatScaling_;
 
    // enable the scaling of the capillary pressure and relative permeability outputs
    // (i.e., rescale the y-Axis)
    bool enablePcScaling_;
    bool enableLeverettScaling_;
    bool enableKrwScaling_;
    bool enableKrnScaling_;
 
    // Employ three-point vertical scaling (e.g., KRWR and KRW).
    bool enableThreePointKrwScaling_{false};
 
    // Employ three-point vertical scaling (e.g., KRORW and KRO).
    bool enableThreePointKrnScaling_{false};
};
 
} // namespace Opm
 
#endif