mirror of
https://github.com/OPM/ResInsight.git
synced 2026-10-25 14:40:26 -05:00
257 lines
13 KiB
C++
257 lines
13 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) 2024- Equinor ASA
|
|
//
|
|
// ResInsight 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.
|
|
//
|
|
// ResInsight 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 at <http://www.gnu.org/licenses/gpl.html>
|
|
// for more details.
|
|
//
|
|
/////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#pragma once
|
|
|
|
#include "cafVecIjk.h"
|
|
|
|
#include "cvfBoundingBox.h"
|
|
#include "cvfStructGrid.h"
|
|
|
|
#include "RigEclipseResultAddress.h"
|
|
|
|
#include <list>
|
|
#include <map>
|
|
#include <optional>
|
|
#include <utility>
|
|
|
|
class RigActiveCellInfo;
|
|
class RigCaseCellResultsData;
|
|
class RigMainGrid;
|
|
class RimEclipseCase;
|
|
class RimEclipseCaseEnsemble;
|
|
class RimRegularGridCase;
|
|
class RigFloodingSettings;
|
|
|
|
//==================================================================================================
|
|
///
|
|
///
|
|
//==================================================================================================
|
|
class RigWellTargetMapping
|
|
{
|
|
public:
|
|
using CellFaceType = cvf::StructGridInterface::FaceType;
|
|
|
|
enum class VolumeType
|
|
{
|
|
OIL,
|
|
GAS,
|
|
HYDROCARBON
|
|
};
|
|
|
|
enum class VolumeResultType
|
|
{
|
|
MOBILE,
|
|
TOTAL
|
|
};
|
|
|
|
enum class VolumesType
|
|
{
|
|
RESERVOIR_VOLUMES,
|
|
SURFACE_VOLUMES_SFIP,
|
|
SURFACE_VOLUMES_FIP,
|
|
RESERVOIR_VOLUMES_COMPUTED
|
|
};
|
|
|
|
struct ClusteringLimits
|
|
{
|
|
double saturationOil;
|
|
double saturationGas;
|
|
double permeability;
|
|
double pressure;
|
|
double transmissibility;
|
|
int maxNumTargets;
|
|
int maxIterations;
|
|
RigEclipseResultAddress filterAddress;
|
|
std::vector<double> filter;
|
|
};
|
|
|
|
struct DataContainer
|
|
{
|
|
std::vector<double> volume;
|
|
std::vector<double> saturationOil;
|
|
std::vector<double> saturationGas;
|
|
std::vector<double> pressure;
|
|
std::vector<double> permeabilityX;
|
|
std::vector<double> permeabilityNNC;
|
|
std::vector<double> transmissibilityX;
|
|
std::vector<double> transmissibilityY;
|
|
std::vector<double> transmissibilityZ;
|
|
const std::vector<double>* transmissibilityNNC;
|
|
};
|
|
|
|
static void generateCandidates( RimEclipseCase* eclipseCase,
|
|
size_t timeStepIdx,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
const RigFloodingSettings& floodingSettings,
|
|
const ClusteringLimits& limits,
|
|
bool skipUndefinedResults );
|
|
|
|
static std::vector<double> getVolumeVector( RigCaseCellResultsData& resultsData,
|
|
RiaDefines::EclipseUnitSystem unitsType,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
size_t timeStepIdx,
|
|
const RigFloodingSettings& floodingSettings );
|
|
|
|
static RimRegularGridCase* generateEnsembleCandidates( RimEclipseCaseEnsemble& ensemble,
|
|
size_t timeStepIdx,
|
|
const cvf::Vec3st& resultGridCellCount,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
const RigFloodingSettings& floodingSettings,
|
|
const ClusteringLimits& limits );
|
|
|
|
static QString wellTargetResultName();
|
|
|
|
class ClusterStatistics
|
|
{
|
|
public:
|
|
ClusterStatistics()
|
|
: id( -1 )
|
|
, numCells( 0 )
|
|
, totalPorvSoil( 0.0 )
|
|
, totalPorvSgas( 0.0 )
|
|
, totalPorvSoilAndSgas( 0.0 )
|
|
, totalFipOil( 0.0 )
|
|
, totalFipGas( 0.0 )
|
|
, totalRfipOil( 0.0 )
|
|
, totalRfipGas( 0.0 )
|
|
, totalSfipOil( 0.0 )
|
|
, totalSfipGas( 0.0 )
|
|
, permeability( 0.0 )
|
|
, pressure( 0.0 )
|
|
{
|
|
}
|
|
|
|
int id;
|
|
size_t numCells;
|
|
double totalPorvSoil;
|
|
double totalPorvSgas;
|
|
double totalPorvSoilAndSgas;
|
|
double totalFipOil;
|
|
double totalFipGas;
|
|
double totalRfipOil;
|
|
double totalRfipGas;
|
|
double totalSfipOil;
|
|
double totalSfipGas;
|
|
double permeability;
|
|
double pressure;
|
|
};
|
|
|
|
private:
|
|
static std::optional<caf::VecIjk0> findStartCell( RimEclipseCase* eclipseCase,
|
|
size_t timeStepIdx,
|
|
const VolumeType volumeType,
|
|
const ClusteringLimits& limits,
|
|
const DataContainer& data,
|
|
const std::vector<double>& filterVector,
|
|
const std::vector<int>& clusters );
|
|
|
|
static void growCluster( RimEclipseCase* eclipseCase,
|
|
const caf::VecIjk0& startCell,
|
|
const VolumeType volumeType,
|
|
const ClusteringLimits& limits,
|
|
const DataContainer& data,
|
|
const std::vector<double>& filterVector,
|
|
std::vector<int>& clusters,
|
|
int clusterId,
|
|
size_t timeStepIdx,
|
|
int maxIterations );
|
|
|
|
static std::vector<size_t> findCandidates( RimEclipseCase* eclipseCase,
|
|
const std::vector<size_t>& previousCells,
|
|
const VolumeType volumeType,
|
|
const ClusteringLimits& limits,
|
|
const DataContainer& data,
|
|
const std::vector<double>& filterVector,
|
|
std::vector<int>& clusters );
|
|
|
|
static void assignClusterIdToCells( const RigActiveCellInfo& activeCellInfo,
|
|
const std::vector<size_t>& cells,
|
|
std::vector<int>& clusters,
|
|
int clusterId );
|
|
|
|
static std::optional<size_t> getActiveCellCount( RimEclipseCase* eclipseCase );
|
|
|
|
static void createResultVector( RimEclipseCase& eclipseCase, const QString& resultName, const std::vector<int>& clusterIds );
|
|
|
|
static void createResultVector( RimEclipseCase& eclipseCase, const QString& resultName, const std::vector<double>& values );
|
|
|
|
static void createResultVector( RimEclipseCase& eclipseCase, const QString& resultName, const std::vector<int>& clusterIds, double value );
|
|
|
|
static void createResultVectorIfDefined( RimEclipseCase& eclipseCase, const QString& resultName, const std::vector<double>& values );
|
|
|
|
static double getValueForFace( const std::vector<double>& x,
|
|
const std::vector<double>& y,
|
|
const std::vector<double>& z,
|
|
CellFaceType face,
|
|
size_t resultIndex );
|
|
|
|
static double getTransmissibilityValueForFace( const std::vector<double>& x,
|
|
const std::vector<double>& y,
|
|
const std::vector<double>& z,
|
|
CellFaceType face,
|
|
size_t resultIndex,
|
|
size_t neighborResultIndex );
|
|
|
|
static std::vector<RigWellTargetMapping::ClusterStatistics> generateStatistics( RimEclipseCase* eclipseCase,
|
|
const std::vector<double>& pressure,
|
|
const std::vector<double>& permeabilityX,
|
|
int numClustersFound,
|
|
size_t timeStepIdx,
|
|
const QString& clusterResultName );
|
|
|
|
static void computeStatisticsAndCreateVectors( RimEclipseCase& targetCase,
|
|
const QString& resultName,
|
|
const std::vector<std::vector<double>>& vec );
|
|
|
|
static void accumulateResultsForSingleCase( RimEclipseCase& eclipseCase,
|
|
RimEclipseCase& targetCase,
|
|
std::map<QString, std::vector<std::vector<double>>>& resultNamesAndSamples,
|
|
std::vector<int>& occupancy );
|
|
|
|
static cvf::BoundingBox computeBoundingBoxForResult( RimEclipseCase& eclipseCase, const QString& resultName, size_t timeStepIndex );
|
|
|
|
static std::list<std::pair<std::pair<size_t, CellFaceType>, size_t>> nncConnectionCellAndResult( size_t cellIdx, RigMainGrid* mainGrid );
|
|
|
|
static std::vector<double> loadVectorByName( RigCaseCellResultsData& resultsData, const QString& resultName, size_t timeStepIdx );
|
|
|
|
static std::vector<double> loadOilVectorByName( RigCaseCellResultsData& resultsData,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
size_t timeStepIdx,
|
|
const RigFloodingSettings& floodingSettings );
|
|
|
|
static std::vector<double> loadGasVectorByName( RigCaseCellResultsData& resultsData,
|
|
RiaDefines::EclipseUnitSystem unitsType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
size_t timeStepIdx,
|
|
const RigFloodingSettings& floodingSettings );
|
|
|
|
static QString getOilVectorName( VolumesType volumesType );
|
|
static QString getGasVectorName( VolumesType volumesType );
|
|
|
|
static bool isSaturationSufficient( const VolumeType volumeType, const DataContainer& data, const ClusteringLimits& limits, size_t idx );
|
|
};
|