mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
1abd0afe6e
Improvements: * Well Target Candidates: Add filtering by binary result * Well Target Candidates: make result grid cell count configurable. * Well Target Candidates: Improve ui ordering. * Well Target Candidates: compute probability, and simplify vector names. * Optimization: Avoid expensive unintended duplication of data
215 lines
11 KiB
C++
215 lines
11 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 <map>
|
|
#include <optional>
|
|
|
|
class RigActiveCellInfo;
|
|
class RigCaseCellResultsData;
|
|
class RimEclipseCase;
|
|
class RimEclipseCaseEnsemble;
|
|
class RimRegularGridCase;
|
|
|
|
//==================================================================================================
|
|
///
|
|
///
|
|
//==================================================================================================
|
|
class RigWellTargetCandidatesGenerator
|
|
{
|
|
public:
|
|
enum class VolumeType
|
|
{
|
|
OIL,
|
|
GAS,
|
|
HYDROCARBON
|
|
};
|
|
|
|
enum class VolumeResultType
|
|
{
|
|
MOBILE,
|
|
TOTAL
|
|
};
|
|
|
|
enum class VolumesType
|
|
{
|
|
RESERVOIR_VOLUMES,
|
|
SURFACE_VOLUMES,
|
|
COMPUTED_VOLUMES
|
|
};
|
|
|
|
struct ClusteringLimits
|
|
{
|
|
double volume;
|
|
double permeability;
|
|
double pressure;
|
|
double transmissibility;
|
|
int maxClusters;
|
|
int maxIterations;
|
|
RigEclipseResultAddress filterAddress;
|
|
};
|
|
|
|
static void generateCandidates( RimEclipseCase* eclipseCase,
|
|
size_t timeStepIdx,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
const ClusteringLimits& limits );
|
|
|
|
static std::vector<double> getVolumeVector( RigCaseCellResultsData& resultsData,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
size_t timeStepIdx );
|
|
|
|
static RimRegularGridCase* generateEnsembleCandidates( RimEclipseCaseEnsemble& ensemble,
|
|
size_t timeStepIdx,
|
|
const cvf::Vec3st& resultGridCellCount,
|
|
VolumeType volumeType,
|
|
VolumesType volumesType,
|
|
VolumeResultType volumeResultType,
|
|
const ClusteringLimits& limits );
|
|
|
|
class ClusterStatistics
|
|
{
|
|
public:
|
|
ClusterStatistics()
|
|
: id( -1 )
|
|
, numCells( 0 )
|
|
, totalPorvSoil( 0.0 )
|
|
, totalPorvSgas( 0.0 )
|
|
, totalPorvSoilAndSgas( 0.0 )
|
|
, totalFipOil( 0.0 )
|
|
, totalFipGas( 0.0 )
|
|
, permeability( 0.0 )
|
|
, pressure( 0.0 )
|
|
{
|
|
}
|
|
|
|
int id;
|
|
size_t numCells;
|
|
double totalPorvSoil;
|
|
double totalPorvSgas;
|
|
double totalPorvSoilAndSgas;
|
|
double totalFipOil;
|
|
double totalFipGas;
|
|
double permeability;
|
|
double pressure;
|
|
};
|
|
|
|
private:
|
|
static std::optional<caf::VecIjk> findStartCell( RimEclipseCase* eclipseCase,
|
|
size_t timeStepIdx,
|
|
const ClusteringLimits& limits,
|
|
const std::vector<double>& volume,
|
|
const std::vector<double>& pressure,
|
|
const std::vector<double>& permeabilityX,
|
|
const std::vector<double>& permeabilityY,
|
|
const std::vector<double>& permeabilityZ,
|
|
const std::vector<double>& transmissibilityX,
|
|
const std::vector<double>& transmissibilityY,
|
|
const std::vector<double>& transmissibilityZ,
|
|
const std::vector<double>& filterVector,
|
|
const std::vector<int>& clusters );
|
|
|
|
static void growCluster( RimEclipseCase* eclipseCase,
|
|
const caf::VecIjk& startCell,
|
|
const ClusteringLimits& limits,
|
|
const std::vector<double>& volume,
|
|
const std::vector<double>& pressure,
|
|
const std::vector<double>& permeabilityX,
|
|
const std::vector<double>& permeabilityY,
|
|
const std::vector<double>& permeabilityZ,
|
|
const std::vector<double>& transmissibilityX,
|
|
const std::vector<double>& transmissibilityY,
|
|
const std::vector<double>& transmissibilityZ,
|
|
const std::vector<double>& filterVector,
|
|
std::vector<int>& clusters,
|
|
int clusterId,
|
|
size_t timeStepIdx,
|
|
int maxIterations );
|
|
|
|
static std::vector<size_t> findCandidates( const RimEclipseCase& eclipseCase,
|
|
const std::vector<size_t>& previousCells,
|
|
const ClusteringLimits& limits,
|
|
const std::vector<double>& volume,
|
|
const std::vector<double>& pressure,
|
|
const std::vector<double>& permeabilityX,
|
|
const std::vector<double>& permeabilityY,
|
|
const std::vector<double>& permeabilityZ,
|
|
const std::vector<double>& transmissibilityX,
|
|
const std::vector<double>& transmissibilityY,
|
|
const std::vector<double>& transmissibilityZ,
|
|
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 double getValueForFace( const std::vector<double>& x,
|
|
const std::vector<double>& y,
|
|
const std::vector<double>& z,
|
|
cvf::StructGridInterface::FaceType face,
|
|
size_t resultIndex );
|
|
|
|
static double getTransmissibilityValueForFace( const std::vector<double>& x,
|
|
const std::vector<double>& y,
|
|
const std::vector<double>& z,
|
|
cvf::StructGridInterface::FaceType face,
|
|
size_t resultIndex,
|
|
size_t neighborResultIndex );
|
|
|
|
static std::vector<RigWellTargetCandidatesGenerator::ClusterStatistics> generateStatistics( RimEclipseCase* eclipseCase,
|
|
const std::vector<double>& pressure,
|
|
const std::vector<double>& permeabilityX,
|
|
const std::vector<double>& permeabilityY,
|
|
const std::vector<double>& permeabilityZ,
|
|
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 );
|
|
};
|