#4048 #4050 3D calculations : Add UI and calculation of case difference

2025-02-25 18:55:39 -06:00 · 2019-02-12 12:48:18 +01:00
parent 93eca21b69
commit 1d4ab1a0b2
7 changed files with 333 additions and 17 deletions
--- a/ApplicationCode/ReservoirDataModel/RigCaseCellResultCalculator.cpp
+++ b/ApplicationCode/ReservoirDataModel/RigCaseCellResultCalculator.cpp
@@ -0,0 +1,162 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  Copyright (C) 2019-    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.
+//
+/////////////////////////////////////////////////////////////////////////////////
+
+#include "RigCaseCellResultCalculator.h"
+
+#include "RiaApplication.h"
+#include "RiaLogging.h"
+
+#include "RigCaseCellResultsData.h"
+#include "RigEclipseCaseData.h"
+#include "RigEclipseResultAddress.h"
+#include "RigGridManager.h"
+#include "RigMainGrid.h"
+#include "RigResultAccessorFactory.h"
+#include "RigResultModifier.h"
+#include "RigResultModifierFactory.h"
+
+#include "RimEclipseCase.h"
+#include "RimProject.h"
+
+#include "cvfAssert.h"
+
+#include <algorithm>
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+bool RigCaseCellResultCalculator::computeDifference(RigEclipseCaseData*            sourceCase,
+                                                    RiaDefines::PorosityModelType  porosityModel,
+                                                    const RigEclipseResultAddress& address)
+{
+    CVF_ASSERT(address.isValid());
+    CVF_ASSERT(address.hasDifferenceCase());
+
+    // Assume at this stage that data for the case is available
+    // It is up to the caller to make sure the case is read from file
+
+    RigEclipseCaseData* otherCase = nullptr;
+
+    {
+        auto eclipseCases = RiaApplication::instance()->project()->eclipseCases();
+        for (RimEclipseCase* c : eclipseCases)
+        {
+            if (c->caseId() == address.m_differenceCaseId)
+            {
+                if (c && c->eclipseCaseData())
+                {
+                    otherCase = c->eclipseCaseData();
+                }
+            }
+        }
+    }
+
+    if (!otherCase || !sourceCase)
+    {
+        RiaLogging::error("Missing input case for difference calculator");
+
+        return false;
+    }
+
+    RigMainGrid* sourceMainGrid = sourceCase->mainGrid();
+    RigMainGrid* otherMainGrid  = otherCase->mainGrid();
+
+    if (!RigGridManager::isEqual(sourceMainGrid, otherMainGrid))
+    {
+        RiaLogging::error("Case difference : Grid cases do not match");
+
+        return false;
+    }
+
+    RigCaseCellResultsData* otherCaseResults  = otherCase->results(porosityModel);
+    RigCaseCellResultsData* sourceCaseResults = sourceCase->results(porosityModel);
+
+    if (!otherCaseResults || !sourceCaseResults)
+    {
+        RiaLogging::error("Missing result data for difference calculator");
+
+        return false;
+    }
+
+    RigEclipseResultAddress nativeAddress(address);
+    nativeAddress.m_differenceCaseId = RigEclipseResultAddress::NO_CASE_DIFF;
+    if (!sourceCaseResults->ensureKnownResultLoaded(nativeAddress))
+    {
+        RiaLogging::error("Failed to load destination diff result");
+
+        return false;
+    }
+
+    if (!otherCaseResults->ensureKnownResultLoaded(nativeAddress))
+    {
+        RiaLogging::error("Failed to load difference result");
+
+        return false;
+    }
+
+    // Initialize difference result with infinity for correct number of time steps and values per time step
+    {
+        const std::vector<std::vector<double>>& srcFrames  = sourceCaseResults->cellScalarResults(nativeAddress);
+        std::vector<std::vector<double>>& diffResultFrames = sourceCaseResults->modifiableCellScalarResultTimesteps(address);
+        diffResultFrames.resize(srcFrames.size());
+        for (size_t fIdx = 0; fIdx < srcFrames.size(); ++fIdx)
+        {
+            const std::vector<double>& srcVals = srcFrames[fIdx];
+            std::vector<double>&       dstVals = diffResultFrames[fIdx];
+
+            dstVals.resize(srcVals.size(), std::numeric_limits<double>::infinity());
+        }
+    }
+
+    size_t otherFrameCount  = otherCaseResults->cellScalarResults(nativeAddress).size();
+    size_t sourceFrameCount = sourceCaseResults->cellScalarResults(nativeAddress).size();
+    size_t maxFrameCount    = std::min(otherFrameCount, sourceFrameCount);
+
+    for (size_t gridIdx = 0; gridIdx < sourceMainGrid->gridCount(); ++gridIdx)
+    {
+        auto                     grid           = sourceMainGrid->gridByIndex(gridIdx);
+        const RigActiveCellInfo* activeCellInfo = sourceCaseResults->activeCellInfo();
+
+        for (size_t fIdx = 0; fIdx < maxFrameCount; ++fIdx)
+        {
+            cvf::ref<RigResultAccessor> sourceResultAccessor =
+                RigResultAccessorFactory::createFromResultAddress(sourceCase, gridIdx, porosityModel, fIdx, nativeAddress);
+
+            cvf::ref<RigResultAccessor> otherResultAccessor =
+                RigResultAccessorFactory::createFromResultAddress(otherCase, gridIdx, porosityModel, fIdx, nativeAddress);
+
+            cvf::ref<RigResultModifier> resultModifier =
+                RigResultModifierFactory::createResultModifier(sourceCase, gridIdx, porosityModel, fIdx, address);
+
+            for (size_t localGridCellIdx = 0; localGridCellIdx < grid->cellCount(); localGridCellIdx++)
+            {
+                size_t reservoirCellIndex = grid->reservoirCellIndex(localGridCellIdx);
+                if (activeCellInfo->isActive(reservoirCellIndex))
+                {
+                    double sourceVal = sourceResultAccessor->cellScalar(localGridCellIdx);
+                    double otherVal  = otherResultAccessor->cellScalar(localGridCellIdx);
+
+                    double difference = otherVal - sourceVal;
+
+                    resultModifier->setCellScalar(localGridCellIdx, difference);
+                }
+            }
+        }
+    }
+
+    return true;
+}