Support for loading only active cell geometry

ApplicationLibCode/Application/RiaPreferencesGrid.cpp

@@ -157,7 +157,8 @@ RifReaderSettings RiaPreferencesGrid::gridOnlyReaderSettings()
         false, // useResultIndexFile
         true, // skipWellData
         false, // import summary data
-        "" // include prefix
+        "", // include prefix,
+        false // only active cells
     };
     return rs;
 }
@@ -174,7 +175,8 @@ RifReaderSettings RiaPreferencesGrid::readerSettings()
                           m_useResultIndexFile,
                           m_skipWellData,
                           true, // import summary data
-                          m_includeFileAbsolutePathPrefix };
+                          m_includeFileAbsolutePathPrefix,
+                          m_onlyLoadActiveCells };
     return rs;
 }
 

ApplicationLibCode/FileInterface/CMakeLists_files.cmake

@@ -86,6 +86,7 @@ set(SOURCE_GROUP_HEADER_FILES
     ${CMAKE_CURRENT_LIST_DIR}/RifRevealCsvSectionSummaryReader.h
     ${CMAKE_CURRENT_LIST_DIR}/RifStimPlanCsvSummaryReader.h
     ${CMAKE_CURRENT_LIST_DIR}/RifReaderOpmCommon.h
+    ${CMAKE_CURRENT_LIST_DIR}/RifReaderOpmCommonActive.h
     ${CMAKE_CURRENT_LIST_DIR}/RifEclipseReportKeywords.h
     ${CMAKE_CURRENT_LIST_DIR}/RifInpExportTools.h
     ${CMAKE_CURRENT_LIST_DIR}/RifFaultReactivationModelExporter.h
@@ -188,6 +189,7 @@ set(SOURCE_GROUP_SOURCE_FILES
     ${CMAKE_CURRENT_LIST_DIR}/RifRevealCsvSectionSummaryReader.cpp
     ${CMAKE_CURRENT_LIST_DIR}/RifStimPlanCsvSummaryReader.cpp
     ${CMAKE_CURRENT_LIST_DIR}/RifReaderOpmCommon.cpp
+    ${CMAKE_CURRENT_LIST_DIR}/RifReaderOpmCommonActive.cpp
     ${CMAKE_CURRENT_LIST_DIR}/RifInpExportTools.cpp
     ${CMAKE_CURRENT_LIST_DIR}/RifFaultReactivationModelExporter.cpp
     ${CMAKE_CURRENT_LIST_DIR}/RifThermalToStimPlanFractureXmlOutput.cpp

ApplicationLibCode/FileInterface/RifReaderEclipseWell.cpp

@@ -132,6 +132,11 @@ size_t RifReaderEclipseWell::localGridCellIndexFromErtConnection( const RigGridB
         return cvf::UNDEFINED_SIZE_T;
     }
 
+    if ( ( cellI < 0 ) || ( cellJ < 0 ) )
+    {
+        return cvf::UNDEFINED_SIZE_T;
+    }
+
     return grid->cellIndexFromIJK( cellI, cellJ, cellK );
 }
 
@@ -158,6 +163,17 @@ RigWellResultPoint RifReaderEclipseWell::createWellResultPoint( const RigEclipse
 
     RigWellResultPoint resultPoint;
 
+    if ( ( grid->cellCount() == 0 ) || ( gridCellIndex > grid->cellCount() - 1 ) )
+    {
+        return resultPoint;
+    }
+
+    const RigCell& c = grid->cell( gridCellIndex );
+    if ( c.isInvalid() )
+    {
+        return resultPoint;
+    }
+
     if ( gridCellIndex != cvf::UNDEFINED_SIZE_T )
     {
         int branchId = -1, segmentId = -1, outletBranchId = -1, outletSegmentId = -1;
@@ -857,11 +873,17 @@ void RifReaderEclipseWell::readWellCells( RifEclipseRestartDataAccess* restartDa
                                 {
                                     prevResPoint = wellResFrame.wellHead();
                                 }
-                                else
+                                else if ( rpIdx > 0 )
                                 {
                                     prevResPoint = wellResultBranch.branchResultPoints()[rpIdx - 1];
                                 }
 
+                                if ( !prevResPoint.isCell() )
+                                {
+                                    // When importing only active cells, this situation can occur if the well head is a inactive cell.
+                                    continue;
+                                }
+
                                 cvf::Vec3d lastConnectionPos = grids[prevResPoint.gridIndex()]->cell( prevResPoint.cellIndex() ).center();
 
                                 SegmentPositionContribution

ApplicationLibCode/FileInterface/RifReaderInterface.cpp

@@ -80,6 +80,14 @@ const QString RifReaderInterface::faultIncludeFileAbsolutePathPrefix() const
     return m_readerSettings.includeFileAbsolutePathPrefix;
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+bool RifReaderInterface::onlyLoadActiveCells() const
+{
+    return m_readerSettings.onlyLoadActiveCells;
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------

ApplicationLibCode/FileInterface/RifReaderInterface.h

@@ -57,6 +57,7 @@ class RifReaderInterface : public cvf::Object
     bool          includeInactiveCellsInFaultGeometry() const;
     bool          loadWellDataEnabled() const;
     const QString faultIncludeFileAbsolutePathPrefix() const;
+    bool          onlyLoadActiveCells() const;
 
     void setReaderSettings( RifReaderSettings readerSettings );
 

ApplicationLibCode/FileInterface/RifReaderOpmCommon.cpp

@@ -29,6 +29,7 @@
 #include "RifOpmRadialGridTools.h"
 #include "RifReaderEclipseWell.h"
 
+#include "RigActiveCellGrid.h"
 #include "RigActiveCellInfo.h"
 #include "RigCaseCellResultsData.h"
 #include "RigEclipseCaseData.h"
@@ -87,18 +88,16 @@ bool RifReaderOpmCommon::open( const QString& fileName, RigEclipseCaseData* ecli
             return false;
         }
 
+        if ( isFaultImportEnabled() )
         {
             auto task = progress.task( "Reading faults", 25 );
 
-            if ( isFaultImportEnabled() )
-            {
-                cvf::Collection<RigFault> faults;
+            cvf::Collection<RigFault> faults;
 
-                importFaults( fileSet, &faults );
+            importFaults( fileSet, &faults );
 
-                RigMainGrid* mainGrid = eclipseCaseData->mainGrid();
-                mainGrid->setFaults( faults );
-            }
+            RigMainGrid* mainGrid = eclipseCaseData->mainGrid();
+            mainGrid->setFaults( faults );
         }
 
         m_eclipseCaseData = eclipseCaseData;
@@ -108,9 +107,10 @@ bool RifReaderOpmCommon::open( const QString& fileName, RigEclipseCaseData* ecli
             buildMetaData( eclipseCaseData, progress );
         }
 
-        auto task = progress.task( "Handling NCC Result data", 25 );
         if ( isNNCsEnabled() )
         {
+            auto task = progress.task( "Handling NCC Result data", 25 );
+
             caf::ProgressInfo nncProgress( 10, "" );
             RigMainGrid*      mainGrid = eclipseCaseData->mainGrid();
 
@@ -278,9 +278,11 @@ bool RifReaderOpmCommon::importGrid( RigMainGrid* mainGrid, RigEclipseCaseData*
             mapAxes[i] = opmMapAxes[i];
         }
 
+        double norm_denominator = mapAxes[2] * mapAxes[5] - mapAxes[4] * mapAxes[3];
+
         // Set the map axes transformation matrix on the main grid
         mainGrid->setMapAxes( mapAxes );
-        mainGrid->setUseMapAxes( true );
+        mainGrid->setUseMapAxes( norm_denominator != 0.0 );
 
         auto transform = mainGrid->mapAxisTransform();
 
@@ -414,6 +416,7 @@ void RifReaderOpmCommon::transferGeometry( Opm::EclIO::EGrid&  opmMainGrid,
 
     RigCell defaultCell;
     defaultCell.setHostGrid( localGrid );
+
     mainGrid->globalCellArray().resize( cellStartIndex + cellCount, defaultCell );
 
     mainGrid->nodes().resize( nodeStartIndex + cellCount * 8, cvf::Vec3d( 0, 0, 0 ) );
@@ -724,12 +727,26 @@ void RifReaderOpmCommon::setupInitAndRestartAccess()
 {
     if ( ( m_initFile == nullptr ) && !m_initFileName.empty() )
     {
-        m_initFile = std::make_unique<EclIO::EInit>( m_initFileName );
+        try
+        {
+            m_initFile = std::make_unique<EclIO::EInit>( m_initFileName );
+        }
+        catch ( ... )
+        {
+            m_initFile = nullptr;
+        }
     }
 
     if ( ( m_restartFile == nullptr ) && !m_restartFileName.empty() )
     {
-        m_restartFile = std::make_unique<EclIO::ERst>( m_restartFileName );
+        try
+        {
+            m_restartFile = std::make_unique<EclIO::ERst>( m_restartFileName );
+        }
+        catch ( ... )
+        {
+            m_restartFile = nullptr;
+        }
     }
 }
 

ApplicationLibCode/FileInterface/RifReaderOpmCommon.h

@@ -32,6 +32,7 @@ class EGrid;
 } // namespace Opm::EclIO
 
 class RigMainGrid;
+class RigActiveCellGrid;
 class RigGridBase;
 class RigEclipseCaseData;
 class RigEclipseTimeStepInfo;
@@ -58,35 +59,37 @@ class RifReaderOpmCommon : public RifReaderInterface
 
     std::vector<QDateTime> timeStepsOnFile( QString gridFileName );
 
-private:
-    void buildMetaData( RigEclipseCaseData* caseData, caf::ProgressInfo& progress );
-    bool importGrid( RigMainGrid* mainGrid, RigEclipseCaseData* caseData );
-    void transferGeometry( Opm::EclIO::EGrid&  opmMainGrid,
-                           Opm::EclIO::EGrid&  opmGrid,
-                           RigMainGrid*        riMainGrid,
-                           RigGridBase*        riGrid,
-                           RigEclipseCaseData* caseData );
+protected:
+    virtual bool importGrid( RigMainGrid* mainGrid, RigEclipseCaseData* caseData );
 
     void transferActiveCells( Opm::EclIO::EGrid&  opmGrid,
                               size_t              cellStartIndex,
                               RigEclipseCaseData* eclipseCaseData,
                               size_t              matrixActiveStartIndex,
                               size_t              fractureActiveStartIndex );
-
     void transferStaticNNCData( Opm::EclIO::EGrid& opmMainGrid, std::vector<Opm::EclIO::EGrid>& lgrGrids, RigMainGrid* mainGrid );
-    void transferDynamicNNCData( RigMainGrid* mainGrid );
 
-    void locateInitAndRestartFiles( QString gridFileName );
-    void setupInitAndRestartAccess();
+    bool verifyActiveCellInfo( int activeSizeMat, int activeSizeFrac );
+    void updateActiveCellInfo( RigEclipseCaseData*             eclipseCaseData,
+                               Opm::EclIO::EGrid&              opmGrid,
+                               std::vector<Opm::EclIO::EGrid>& lgrGrids,
+                               RigMainGrid*                    mainGrid );
+
+private:
+    void buildMetaData( RigEclipseCaseData* caseData, caf::ProgressInfo& progress );
 
     std::vector<RigEclipseTimeStepInfo> createFilteredTimeStepInfos();
+    std::vector<std::vector<int>>       readActiveCellInfoFromPorv( RigEclipseCaseData* eclipseCaseData, bool isDualPorosity );
 
-    bool                          verifyActiveCellInfo( int activeSizeMat, int activeSizeFrac );
-    std::vector<std::vector<int>> readActiveCellInfoFromPorv( RigEclipseCaseData* eclipseCaseData, bool isDualPorosity );
-    void                          updateActiveCellInfo( RigEclipseCaseData*             eclipseCaseData,
-                                                        Opm::EclIO::EGrid&              opmGrid,
-                                                        std::vector<Opm::EclIO::EGrid>& lgrGrids,
-                                                        RigMainGrid*                    mainGrid );
+    void transferGeometry( Opm::EclIO::EGrid&  opmMainGrid,
+                           Opm::EclIO::EGrid&  opmGrid,
+                           RigMainGrid*        riMainGrid,
+                           RigGridBase*        riGrid,
+                           RigEclipseCaseData* caseData );
+    void transferDynamicNNCData( RigMainGrid* mainGrid );
+
+    void locateInitAndRestartFiles( QString gridFileName );
+    void setupInitAndRestartAccess();
 
     struct TimeDataFile
     {
@@ -99,22 +102,22 @@ class RifReaderOpmCommon : public RifReaderInterface
 
     std::vector<TimeDataFile> readTimeSteps();
 
-private:
+protected:
     enum class ActiveType
     {
         ACTIVE_MATRIX_VALUE   = 1,
         ACTIVE_FRACTURE_VALUE = 2
     };
 
-    std::string m_gridFileName;
-    std::string m_initFileName;
-    std::string m_restartFileName;
-    int         m_gridUnit;
+    std::string              m_gridFileName;
+    int                      m_gridUnit;
+    std::vector<std::string> m_gridNames;
 
     RigEclipseCaseData* m_eclipseCaseData;
 
+private:
+    std::string                        m_initFileName;
+    std::string                        m_restartFileName;
     std::unique_ptr<Opm::EclIO::ERst>  m_restartFile;
     std::unique_ptr<Opm::EclIO::EInit> m_initFile;
-
-    std::vector<std::string> m_gridNames;
 };