mirror of
https://github.com/OPM/opm-simulators.git
synced 2025-02-25 18:55:30 -06:00
fbb24e2a5a
Check group limits in gas lift stage 1 to avoid adding too much ALQ which must anyway later be removed in stage 2. This should make the optimization more efficient for small ALQ increment values. Also adds MPI support.
196 lines
6.9 KiB
C++
196 lines
6.9 KiB
C++
// -*- 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 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/>.
|
|
|
|
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.
|
|
*/
|
|
#include "config.h"
|
|
|
|
#define BOOST_TEST_MODULE Glift1
|
|
|
|
#include <opm/models/utils/propertysystem.hh>
|
|
#include <opm/models/utils/parametersystem.hh>
|
|
#include <ebos/equil/equilibrationhelpers.hh>
|
|
#include <ebos/eclproblem.hh>
|
|
#include <ebos/ebos.hh>
|
|
#include <opm/models/utils/start.hh>
|
|
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Schedule.hpp>
|
|
#include <opm/parser/eclipse/EclipseState/Schedule/Well/Well.hpp>
|
|
#include <opm/simulators/utils/DeferredLogger.hpp>
|
|
#include <opm/simulators/flow/BlackoilModelEbos.hpp>
|
|
#include <opm/simulators/wells/BlackoilWellModel.hpp>
|
|
#include <opm/simulators/wells/StandardWell.hpp>
|
|
#include <opm/simulators/wells/GasLiftSingleWell.hpp>
|
|
#include <opm/simulators/wells/GasLiftSingleWellGeneric.hpp>
|
|
#include <opm/simulators/wells/GasLiftGroupInfo.hpp>
|
|
//#include <opm/simulators/flow/SimulatorFullyImplicitBlackoilEbos.hpp>
|
|
//#include <flow/flow_ebos_blackoil.hpp>
|
|
#include <opm/simulators/wells/WellState.hpp>
|
|
/// #include <opm/simulators/flow/Main.hpp>
|
|
|
|
#if HAVE_DUNE_FEM
|
|
#include <dune/fem/misc/mpimanager.hh>
|
|
#else
|
|
#include <dune/common/parallel/mpihelper.hh>
|
|
#endif
|
|
|
|
#include <exception>
|
|
#include <iostream>
|
|
#include <memory>
|
|
#include <optional>
|
|
#include <string>
|
|
#include <vector>
|
|
|
|
#include <boost/test/unit_test.hpp>
|
|
#include <boost/version.hpp>
|
|
#if BOOST_VERSION / 100000 == 1 && BOOST_VERSION / 100 % 1000 < 71
|
|
#include <boost/test/floating_point_comparison.hpp>
|
|
#else
|
|
#include <boost/test/tools/floating_point_comparison.hpp>
|
|
#endif
|
|
|
|
namespace Opm::Properties {
|
|
namespace TTag {
|
|
struct TestGliftTypeTag {
|
|
using InheritsFrom = std::tuple<EbosTypeTag>;
|
|
};
|
|
}
|
|
}
|
|
|
|
template <class TypeTag>
|
|
std::unique_ptr<Opm::GetPropType<TypeTag, Opm::Properties::Simulator>>
|
|
initSimulator(const char *filename)
|
|
{
|
|
using Simulator = Opm::GetPropType<TypeTag, Opm::Properties::Simulator>;
|
|
|
|
std::string filename_arg = "--ecl-deck-file-name=";
|
|
filename_arg += filename;
|
|
|
|
const char* argv[] = {
|
|
"test_equil",
|
|
filename_arg.c_str()
|
|
};
|
|
|
|
Opm::setupParameters_<TypeTag>(/*argc=*/sizeof(argv)/sizeof(argv[0]), argv, /*registerParams=*/true);
|
|
|
|
return std::unique_ptr<Simulator>(new Simulator);
|
|
}
|
|
|
|
|
|
namespace {
|
|
|
|
struct GliftFixture {
|
|
GliftFixture() {
|
|
int argc = boost::unit_test::framework::master_test_suite().argc;
|
|
char** argv = boost::unit_test::framework::master_test_suite().argv;
|
|
#if HAVE_DUNE_FEM
|
|
Dune::Fem::MPIManager::initialize(argc, argv);
|
|
#else
|
|
Dune::MPIHelper::instance(argc, argv);
|
|
#endif
|
|
using TypeTag = Opm::Properties::TTag::EclFlowProblem;
|
|
Opm::registerAllParameters_<TypeTag>();
|
|
}
|
|
};
|
|
|
|
}
|
|
|
|
BOOST_GLOBAL_FIXTURE(GliftFixture);
|
|
|
|
BOOST_AUTO_TEST_CASE(G1)
|
|
{
|
|
//using TypeTag = Opm::Properties::TTag::EclFlowProblem;
|
|
using TypeTag = Opm::Properties::TTag::TestGliftTypeTag;
|
|
//using EclProblem = Opm::EclProblem<TypeTag>;
|
|
//using EclWellModel = typename EclProblem::EclWellModel;
|
|
using WellModel = Opm::BlackoilWellModel<TypeTag>;
|
|
using WellState = Opm::WellState;
|
|
using StdWell = Opm::StandardWell<TypeTag>;
|
|
using GasLiftSingleWell = Opm::GasLiftSingleWell<TypeTag>;
|
|
using GasLiftGroupInfo = Opm::GasLiftGroupInfo;
|
|
using GasLiftSingleWellGeneric = Opm::GasLiftSingleWellGeneric;
|
|
using GLiftEclWells = typename GasLiftGroupInfo::GLiftEclWells;
|
|
const std::string filename = "GLIFT1.DATA";
|
|
using GLiftSyncGroups = typename GasLiftSingleWellGeneric::GLiftSyncGroups;
|
|
|
|
auto simulator = initSimulator<TypeTag>(filename.data());
|
|
|
|
simulator->model().applyInitialSolution();
|
|
simulator->setEpisodeIndex(-1);
|
|
simulator->setEpisodeLength(0.0);
|
|
simulator->startNextEpisode(/*episodeStartTime=*/0.0, /*episodeLength=*/1e30);
|
|
simulator->setTimeStepSize(43200); // 12 hours
|
|
simulator->model().newtonMethod().setIterationIndex(0);
|
|
WellModel& well_model = simulator->problem().wellModel();
|
|
int report_step_idx = 0;
|
|
well_model.beginReportStep(report_step_idx);
|
|
well_model.beginTimeStep();
|
|
well_model.updatePerforationIntensiveQuantities();
|
|
Opm::DeferredLogger deferred_logger;
|
|
well_model.calculateExplicitQuantities(deferred_logger);
|
|
well_model.prepareTimeStep(deferred_logger);
|
|
well_model.updateWellControls(deferred_logger, /* check group controls */ true);
|
|
well_model.initPrimaryVariablesEvaluation();
|
|
Opm::WellInterface<TypeTag> *well_ptr = well_model.getWell("B-1H").get();
|
|
StdWell *std_well = dynamic_cast<StdWell *>(well_ptr);
|
|
|
|
const auto& schedule = simulator->vanguard().schedule();
|
|
auto wells_ecl = schedule.getWells(report_step_idx);
|
|
std::optional<std::size_t> idx;
|
|
int num_producers = 0;
|
|
for(std::size_t i = 0; i < wells_ecl.size(); ++i) {
|
|
const auto &well = wells_ecl[i];
|
|
if (well.isProducer()) {
|
|
idx = i;
|
|
num_producers++;
|
|
}
|
|
}
|
|
BOOST_CHECK_EQUAL( num_producers, 1);
|
|
const auto &well = wells_ecl[*idx];
|
|
BOOST_CHECK_EQUAL( well.name(), "B-1H");
|
|
const auto& summary_state = simulator->vanguard().summaryState();
|
|
WellState &well_state = well_model.wellState();
|
|
GLiftEclWells ecl_well_map;
|
|
well_model.initGliftEclWellMap(ecl_well_map);
|
|
const int iteration_idx = simulator->model().newtonMethod().numIterations();
|
|
GasLiftGroupInfo group_info {
|
|
ecl_well_map,
|
|
schedule,
|
|
summary_state,
|
|
simulator->episodeIndex(),
|
|
iteration_idx,
|
|
well_model.phaseUsage(),
|
|
deferred_logger,
|
|
well_state
|
|
};
|
|
GLiftSyncGroups sync_groups;
|
|
GasLiftSingleWell glift {*std_well, *(simulator.get()), summary_state,
|
|
deferred_logger, well_state, group_info, sync_groups};
|
|
auto state = glift.runOptimize(iteration_idx);
|
|
BOOST_CHECK_CLOSE(state->oilRate(), 0.01736111111111111, 1e-8);
|
|
BOOST_CHECK_CLOSE(state->gasRate(), 1.6464646999768586, 1e-8);
|
|
BOOST_CHECK(state->oilIsLimited());
|
|
BOOST_CHECK(!state->gasIsLimited());
|
|
BOOST_CHECK(!state->alqIsLimited());
|
|
BOOST_CHECK_CLOSE(state->alq(), 0.0, 1e-8);
|
|
BOOST_CHECK(!state->increase().has_value());
|
|
}
|
|
|