OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

WIP: create occupancy grid

Browse Source
This commit is contained in:
Kristian Bendiksen 2024-10-04 13:24:53 +02:00
parent 3351cacd33
commit 20db9141f2
1 changed files with 11 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
ApplicationLibCode/ReservoirDataModel/Well/RigWellTargetCandidatesGenerator.cpp
View File

@ -614,6 +614,8 @@ void RigWellTargetCandidatesGenerator::generateEnsembleCandidates( RimEclipseCas
VolumeResultType volumeResultType, VolumeResultType volumeResultType,
const ClusteringLimits& limits ) const ClusteringLimits& limits )
{ {
RiaLogging::debug( "Generating ensemble statistics" );
for ( auto eclipseCase : ensemble.cases() ) for ( auto eclipseCase : ensemble.cases() )
{ {
generateCandidates( eclipseCase, timeStepIdx, volumeType, volumesType, volumeResultType, limits ); generateCandidates( eclipseCase, timeStepIdx, volumeType, volumesType, volumeResultType, limits );
@ -624,32 +626,32 @@ void RigWellTargetCandidatesGenerator::generateEnsembleCandidates( RimEclipseCas
const size_t targetNumReservoirCells = targetActiveCellInfo->reservoirCellCount(); const size_t targetNumReservoirCells = targetActiveCellInfo->reservoirCellCount();
const size_t targetNumActiveCells = targetActiveCellInfo->reservoirActiveCellCount(); const size_t targetNumActiveCells = targetActiveCellInfo->reservoirActiveCellCount();
const int numGrids = ensemble.cases().size();
std::vector<int> occupancy( targetNumActiveCells, 0 ); std::vector<int> occupancy( targetNumActiveCells, 0 );
std::vector<double> pressure( numGrids, std::numeric_limits<double>::infinity() );
for ( auto eclipseCase : ensemble.cases() ) for ( auto eclipseCase : ensemble.cases() )
{ {
RigCaseCellResultsData* resultsData = eclipseCase->results( RiaDefines::PorosityModelType::MATRIX_MODEL ); RigCaseCellResultsData* resultsData = eclipseCase->results( RiaDefines::PorosityModelType::MATRIX_MODEL );
const RigMainGrid* mainGrid = eclipseCase->mainGrid(); const RigMainGrid* mainGrid = eclipseCase->mainGrid();
const RigActiveCellInfo* activeCellInfo = resultsData->activeCellInfo(); const RigActiveCellInfo* activeCellInfo = resultsData->activeCellInfo();
RigEclipseResultAddress pressureAddress( RiaDefines::ResultCatType::GENERATED, "CLUSTER_NUM" ); RigEclipseResultAddress pressureAddress( RiaDefines::ResultCatType::GENERATED, "CLUSTERS_NUM" );
resultsData->ensureKnownResultLoaded( pressureAddress ); resultsData->ensureKnownResultLoaded( pressureAddress );
const std::vector<double>& clusterNum = resultsData->cellScalarResults( pressureAddress, timeStepIdx ); const std::vector<double>& clusterNum = resultsData->cellScalarResults( pressureAddress, 0 );
for ( size_t targetCellIdx = 0; targetCellIdx < targetNumReservoirCells; targetCellIdx++ ) for ( size_t targetCellIdx = 0; targetCellIdx < targetNumReservoirCells; targetCellIdx++ )
{ {
const RigCell& nativeCell = targetCase.mainGrid()->globalCellArray()[targetCellIdx]; const RigCell& nativeCell = targetCase.mainGrid()->globalCellArray()[targetCellIdx];
cvf::Vec3d cellCenter = nativeCell.center(); cvf::Vec3d cellCenter = nativeCell.center();
size_t targetResultIndex = targetActiveCellInfo->cellResultIndex( targetCellIdx );
size_t cellIdx = mainGrid->findReservoirCellIndexFromPoint( cellCenter ); size_t cellIdx = mainGrid->findReservoirCellIndexFromPoint( cellCenter );
if ( cellIdx != cvf::UNDEFINED_SIZE_T && activeCellInfo->isActive( cellIdx ) ) if ( cellIdx != cvf::UNDEFINED_SIZE_T && activeCellInfo->isActive( cellIdx ) && targetResultIndex != cvf::UNDEFINED_SIZE_T )
{ {
size_t resultIndex = targetResultsData->activeCellInfo()->cellResultIndex( cellIdx ); size_t resultIndex = resultsData->activeCellInfo()->cellResultIndex( cellIdx );
if ( clusterNum[resultIndex] > 0 ) if ( !std::isinf( clusterNum[resultIndex] ) && clusterNum[resultIndex] > 0 )
{ {
occupancy[resultIndex]++; occupancy[targetResultIndex]++;
} }
} }
} }