#5472 Add indication of water depth to formation and casing track

* Also add lighter formation lines as an option and make this default for WBS plots

ApplicationCode/Application/Tools/RiaColorTables.cpp

@@ -484,6 +484,20 @@ const caf::ColorTable& RiaColorTables::wellPathsPaletteColors()
     return colorTable;
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+const caf::ColorTable& RiaColorTables::waterAndRockPaletteColors()
+{
+    static std::vector<cvf::Color3ub> colors{
+        cvf::Color3ub( 127, 205, 255 ), // Sea Blue
+        cvf::Color3ub( 220, 220, 220 ) // Light Gray
+    };
+    static caf::ColorTable colorTable = caf::ColorTable( colors );
+
+    return colorTable;
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------

ApplicationCode/Application/Tools/RiaColorTables.h

@@ -59,6 +59,7 @@ public:
     static const caf::ColorTable& timestepsPaletteColors();
     static const caf::ColorTable& editableWellPathsPaletteColors();
     static const caf::ColorTable& wellPathsPaletteColors();
+    static const caf::ColorTable& waterAndRockPaletteColors();
 
     static cvf::Color3f undefinedCellColor();
 

ApplicationCode/Commands/WellLogCommands/RicNewWellBoreStabilityPlotFeature.cpp

@@ -311,7 +311,8 @@ void RicNewWellBoreStabilityPlotFeature::createStabilityCurvesTrack( RimWellBore
     stabilityCurvesTrack->setXAxisGridVisibility( RimWellLogPlot::AXIS_GRID_MAJOR_AND_MINOR );
     stabilityCurvesTrack->setFormationWellPath( wellPath );
     stabilityCurvesTrack->setFormationCase( geoMechCase );
-    stabilityCurvesTrack->setAnnotationType( RiuPlotAnnotationTool::NO_ANNOTATIONS );
+    stabilityCurvesTrack->setAnnotationType( RiuPlotAnnotationTool::FORMATION_ANNOTATIONS );
+    stabilityCurvesTrack->setAnnotationDisplay( RiuPlotAnnotationTool::LIGHT_LINES );
     stabilityCurvesTrack->setShowRegionLabels( true );
 
     std::vector<QString> resultNames = RiaDefines::wbsDerivedResultNames();
@@ -414,6 +415,7 @@ void RicNewWellBoreStabilityPlotFeature::createAnglesTrack( RimWellBoreStability
     wellPathAnglesTrack->setXAxisGridVisibility( RimWellLogPlot::AXIS_GRID_MAJOR_AND_MINOR );
     wellPathAnglesTrack->setFormationWellPath( wellPath );
     wellPathAnglesTrack->setFormationCase( geoMechCase );
-    wellPathAnglesTrack->setAnnotationType( RiuPlotAnnotationTool::NO_ANNOTATIONS );
+    wellPathAnglesTrack->setAnnotationType( RiuPlotAnnotationTool::FORMATION_ANNOTATIONS );
+    wellPathAnglesTrack->setAnnotationDisplay( RiuPlotAnnotationTool::LIGHT_LINES );
     wellPathAnglesTrack->setShowRegionLabels( false );
 }

ApplicationCode/ProjectDataModel/RimWellLogTrack.cpp

@@ -139,6 +139,7 @@ template <>
 void AppEnum<RiuPlotAnnotationTool::RegionDisplay>::setUp()
 {
     addItem( RiuPlotAnnotationTool::DARK_LINES, "DARK_LINES", "Dark Lines" );
+    addItem( RiuPlotAnnotationTool::LIGHT_LINES, "LIGHT_LINES", "Light Lines" );
     addItem( RiuPlotAnnotationTool::COLORED_LINES, "COLORED_LINES", "Colored Lines" );
     addItem( RiuPlotAnnotationTool::COLOR_SHADING, "COLOR_SHADING", "Color Shading" );
     addItem( RiuPlotAnnotationTool::COLOR_SHADING_AND_LINES, "SHADING_AND_LINES", "Color Shading and Lines" );
@@ -1487,6 +1488,14 @@ void RimWellLogTrack::setAnnotationDisplay( RiuPlotAnnotationTool::RegionDisplay
     m_regionAnnotationDisplay = annotationDisplay;
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+void RimWellLogTrack::setAnnotationTransparency( int percent )
+{
+    m_colorShadingTransparency = percent;
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
@@ -1535,6 +1544,14 @@ void RimWellLogTrack::setShowWellPathAttributes( bool on )
     m_showWellPathAttributes = on;
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+void RimWellLogTrack::setShowBothSidesOfWell( bool on )
+{
+    m_showWellPathComponentsBothSides = on;
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
@@ -1840,6 +1857,40 @@ void RimWellLogTrack::handleWheelEvent( QWheelEvent* event )
     }
 }
 
+//--------------------------------------------------------------------------------------------------
+///
+//--------------------------------------------------------------------------------------------------
+std::vector<std::pair<double, double>> RimWellLogTrack::waterAndRockRegions( RiaDefines::DepthTypeEnum  depthType,
+                                                                             const RigWellLogExtractor* extractor ) const
+{
+    if ( depthType == RiaDefines::MEASURED_DEPTH )
+    {
+        double waterStartMD = 0.0;
+        if ( extractor->wellPathData()->rkbDiff() != std::numeric_limits<double>::infinity() )
+        {
+            waterStartMD += extractor->wellPathData()->rkbDiff();
+        }
+        double waterEndMD = extractor->cellIntersectionMDs().front();
+        double rockEndMD  = extractor->cellIntersectionMDs().back();
+        return {{waterStartMD, waterEndMD}, {waterEndMD, rockEndMD}};
+    }
+    else if ( depthType == RiaDefines::TRUE_VERTICAL_DEPTH )
+    {
+        double waterStartTVD = 0.0;
+        double waterEndTVD   = extractor->cellIntersectionTVDs().front();
+        double rockEndTVD    = extractor->cellIntersectionTVDs().back();
+        return {{waterStartTVD, waterEndTVD}, {waterEndTVD, rockEndTVD}};
+    }
+    else if ( depthType == RiaDefines::TRUE_VERTICAL_DEPTH_RKB )
+    {
+        double waterStartTVDRKB = extractor->wellPathData()->rkbDiff();
+        double waterEndTVDRKB   = extractor->cellIntersectionTVDs().front() + extractor->wellPathData()->rkbDiff();
+        double rockEndTVDRKB    = extractor->cellIntersectionTVDs().back() + extractor->wellPathData()->rkbDiff();
+        return {{waterStartTVDRKB, waterEndTVDRKB}, {waterEndTVDRKB, rockEndTVDRKB}};
+    }
+    return {};
+}
+
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
@@ -2174,17 +2225,43 @@ void RimWellLogTrack::updateRegionAnnotationsOnPlot()
 //--------------------------------------------------------------------------------------------------
 void RimWellLogTrack::updateFormationNamesOnPlot()
 {
-    std::vector<QString> formationNamesToPlot;
-
     RimWellLogPlot* plot = nullptr;
     firstAncestorOrThisOfTypeAsserted( plot );
 
-    if ( m_formationSource == CASE )
+    if ( m_formationSource() == WELL_PICK_FILTER )
     {
-        if ( ( m_formationSimWellName == QString( "None" ) && m_formationWellPathForSourceCase == nullptr ) ||
+        if ( m_formationWellPathForSourceWellPath == nullptr ) return;
-             m_formationCase == nullptr )
-            return;
 
+        if ( !( plot->depthType() == RiaDefines::MEASURED_DEPTH || plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH ||
+                plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH_RKB ) )
+        {
+            return;
+        }
+
+        std::vector<double> yValues;
+
+        const RigWellPathFormations* formations = m_formationWellPathForSourceWellPath->formationsGeometry();
+        if ( !formations ) return;
+
+        std::vector<QString> formationNamesToPlot;
+        formations->depthAndFormationNamesUpToLevel( m_formationLevel(),
+                                                     &formationNamesToPlot,
+                                                     &yValues,
+                                                     m_showformationFluids(),
+                                                     plot->depthType() );
+
+        if ( plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH_RKB )
+        {
+            for ( double& depthValue : yValues )
+            {
+                depthValue += m_formationWellPathForSourceWellPath->wellPathGeometry()->rkbDiff();
+            }
+        }
+
+        m_annotationTool->attachWellPicks( m_plotWidget, formationNamesToPlot, yValues );
+    }
+    else
+    {
         RimMainPlotCollection* mainPlotCollection;
         this->firstAncestorOrThisOfTypeAsserted( mainPlotCollection );
 
@@ -2194,6 +2271,7 @@ void RimWellLogTrack::updateFormationNamesOnPlot()
 
         RigEclipseWellLogExtractor* eclWellLogExtractor     = nullptr;
         RigGeoMechWellLogExtractor* geoMechWellLogExtractor = nullptr;
+        RigWellLogExtractor*        extractor               = nullptr;
 
         if ( m_formationTrajectoryType == SIMULATION_WELL )
         {
@@ -2222,6 +2300,7 @@ void RimWellLogTrack::updateFormationNamesOnPlot()
                                                                   RiaDefines::activeFormationNamesResultName() ) );
 
             curveData = RimWellLogTrack::curveSamplingPointData( eclWellLogExtractor, resultAccessor.p() );
+            extractor = eclWellLogExtractor;
         }
         else
         {
@@ -2235,60 +2314,52 @@ void RimWellLogTrack::updateFormationNamesOnPlot()
                                                                  RigFemResultAddress( RIG_FORMATION_NAMES,
                                                                                       activeFormationNamesResultName,
                                                                                       "" ) );
+            extractor = geoMechWellLogExtractor;
         }
 
-        std::vector<std::pair<double, double>> yValues;
+        // Attach water and rock base formations
-        std::vector<QString> formationNamesVector = RimWellLogTrack::formationNamesVector( m_formationCase );
+        const std::pair<double, double> xRange = std::make_pair( m_visibleXRangeMin(), m_visibleXRangeMax() );
-
-        RimWellLogTrack::findRegionNamesToPlot( curveData,
-                                                formationNamesVector,
-                                                plot->depthType(),
-                                                &formationNamesToPlot,
-                                                &yValues );
-
-        std::pair<double, double> xRange = std::make_pair( m_visibleXRangeMin(), m_visibleXRangeMax() );
-
-        caf::ColorTable colorTable( RimRegularLegendConfig::colorArrayFromColorType( m_colorShadingPalette() ) );
 
+        const caf::ColorTable                        waterAndRockColors = RiaColorTables::waterAndRockPaletteColors();
+        const std::vector<std::pair<double, double>> waterAndRockIntervals    = waterAndRockRegions( plot->depthType(),
+                                                                                                  extractor );
+        int                                          waterAndRockTransparency = m_colorShadingTransparency / 2;
         m_annotationTool->attachNamedRegions( m_plotWidget,
-                                              formationNamesToPlot,
+                                              {"Water", ""},
                                               xRange,
-                                              yValues,
+                                              waterAndRockIntervals,
                                               m_regionAnnotationDisplay(),
-                                              colorTable,
+                                              waterAndRockColors,
-                                              ( ( 100 - m_colorShadingTransparency ) * 255 ) / 100,
+                                              ( ( 100 - waterAndRockTransparency ) * 255 ) / 100,
                                               m_showRegionLabels() );
-    }
-    else if ( m_formationSource() == WELL_PICK_FILTER )
-    {
-        if ( m_formationWellPathForSourceWellPath == nullptr ) return;
 
-        if ( !( plot->depthType() == RiaDefines::MEASURED_DEPTH || plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH ||
+        if ( m_formationSource == CASE )
-                plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH_RKB ) )
         {
-            return;
+            if ( ( m_formationSimWellName == QString( "None" ) && m_formationWellPathForSourceCase == nullptr ) ||
+                 m_formationCase == nullptr )
+                return;
+
+            std::vector<std::pair<double, double>> yValues;
+            std::vector<QString> formationNamesVector = RimWellLogTrack::formationNamesVector( m_formationCase );
+
+            std::vector<QString> formationNamesToPlot;
+            RimWellLogTrack::findRegionNamesToPlot( curveData,
+                                                    formationNamesVector,
+                                                    plot->depthType(),
+                                                    &formationNamesToPlot,
+                                                    &yValues );
+
+            caf::ColorTable colorTable( RimRegularLegendConfig::colorArrayFromColorType( m_colorShadingPalette() ) );
+
+            m_annotationTool->attachNamedRegions( m_plotWidget,
+                                                  formationNamesToPlot,
+                                                  xRange,
+                                                  yValues,
+                                                  m_regionAnnotationDisplay(),
+                                                  colorTable,
+                                                  ( ( 100 - m_colorShadingTransparency ) * 255 ) / 100,
+                                                  m_showRegionLabels() );
         }
-
-        std::vector<double> yValues;
-
-        const RigWellPathFormations* formations = m_formationWellPathForSourceWellPath->formationsGeometry();
-        if ( !formations ) return;
-
-        formations->depthAndFormationNamesUpToLevel( m_formationLevel(),
-                                                     &formationNamesToPlot,
-                                                     &yValues,
-                                                     m_showformationFluids(),
-                                                     plot->depthType() );
-
-        if ( plot->depthType() == RiaDefines::TRUE_VERTICAL_DEPTH_RKB )
-        {
-            for ( double& depthValue : yValues )
-            {
-                depthValue += m_formationWellPathForSourceWellPath->wellPathGeometry()->rkbDiff();
-            }
-        }
-
-        m_annotationTool->attachWellPicks( m_plotWidget, formationNamesToPlot, yValues );
     }
 }
 

ApplicationCode/ProjectDataModel/RimWellLogTrack.h

@@ -53,6 +53,7 @@ class RimWellLogPlotCollection;
 class RigGeoMechWellLogExtractor;
 class RigResultAccessor;
 class RigFemResultAddress;
+class RigWellLogExtractor;
 
 class QwtPlotCurve;
 
@@ -163,6 +164,7 @@ public:
 
     void setAnnotationType( RiuPlotAnnotationTool::RegionAnnotationType annotationType );
     void setAnnotationDisplay( RiuPlotAnnotationTool::RegionDisplay annotationDisplay );
+    void setAnnotationTransparency( int percent );
 
     RiuPlotAnnotationTool::RegionAnnotationType annotationType() const;
     RiuPlotAnnotationTool::RegionDisplay        annotationDisplay() const;
@@ -173,6 +175,7 @@ public:
 
     bool showWellPathAttributes() const;
     void setShowWellPathAttributes( bool on );
+    void setShowBothSidesOfWell( bool on );
     void setWellPathAttributesSource( RimWellPath* wellPath );
 
     RimWellPath* wellPathAttributeSource() const;
@@ -279,6 +282,9 @@ private:
 
     void handleWheelEvent( QWheelEvent* event ) override;
 
+    std::vector<std::pair<double, double>> waterAndRockRegions( RiaDefines::DepthTypeEnum  depthType,
+                                                                const RigWellLogExtractor* extractor ) const;
+
 private:
     QString m_xAxisTitle;
 

ApplicationCode/ReservoirDataModel/RigWellLogExtractor.cpp

@@ -39,7 +39,7 @@ RigWellLogExtractor::~RigWellLogExtractor() {}
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
-const std::vector<double>& RigWellLogExtractor::cellIntersectionMDs()
+const std::vector<double>& RigWellLogExtractor::cellIntersectionMDs() const
 {
     return m_intersectionMeasuredDepths;
 }
@@ -47,7 +47,7 @@ const std::vector<double>& RigWellLogExtractor::cellIntersectionMDs()
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
-const std::vector<double>& RigWellLogExtractor::cellIntersectionTVDs()
+const std::vector<double>& RigWellLogExtractor::cellIntersectionTVDs() const
 {
     return m_intersectionTVDs;
 }
@@ -88,7 +88,7 @@ std::vector<WellPathCellIntersectionInfo> RigWellLogExtractor::cellIntersectionI
 //--------------------------------------------------------------------------------------------------
 ///
 //--------------------------------------------------------------------------------------------------
-const std::vector<size_t>& RigWellLogExtractor::intersectedCellsGlobIdx()
+const std::vector<size_t>& RigWellLogExtractor::intersectedCellsGlobIdx() const
 {
     return m_intersectedCellsGlobIdx;
 }

ApplicationCode/ReservoirDataModel/RigWellLogExtractor.h

@@ -58,9 +58,9 @@ public:
     RigWellLogExtractor( const RigWellPath* wellpath, const std::string& wellCaseErrorMsgName );
     ~RigWellLogExtractor() override;
 
-    const std::vector<double>& cellIntersectionMDs();
+    const std::vector<double>& cellIntersectionMDs() const;
-    const std::vector<double>& cellIntersectionTVDs();
+    const std::vector<double>& cellIntersectionTVDs() const;
-    const std::vector<size_t>& intersectedCellsGlobIdx();
+    const std::vector<size_t>& intersectedCellsGlobIdx() const;
 
     const RigWellPath* wellPathData() const;
 

ApplicationCode/UserInterface/RiuPlotAnnotationTool.cpp

@@ -97,12 +97,23 @@ void RiuPlotAnnotationTool::attachNamedRegions( QwtPlot*
 
         QColor lineColor( 0, 0, 0, 0 );
         QColor textColor( 0, 0, 0, 255 );
-        if ( regionDisplay & DARK_LINES || regionDisplay & COLORED_LINES )
+        if ( regionDisplay & DARK_LINES || regionDisplay & COLORED_LINES || regionDisplay & LIGHT_LINES )
         {
             cvf::Color3ub cvfColor = catMapper.mapToColor( static_cast<double>( i ) );
             QColor        cycledColor( cvfColor.r(), cvfColor.g(), cvfColor.b() );
 
-            lineColor = regionDisplay & DARK_LINES ? QColor( 0, 0, 100 ) : cycledColor;
+            if ( regionDisplay & DARK_LINES )
+            {
+                lineColor = QColor( 50, 50, 100 );
+            }
+            else if ( regionDisplay & LIGHT_LINES )
+            {
+                lineColor = QColor( 200, 200, 200 );
+            }
+            else
+            {
+                lineColor = cycledColor;
+            }
             textColor = lineColor;
         }
         Qt::Alignment horizontalAlignment = trackTextAlignment( trackSpan );

ApplicationCode/UserInterface/RiuPlotAnnotationTool.h

@@ -44,7 +44,8 @@ public:
         DARK_LINES              = 0x01,
         COLORED_LINES           = 0x02,
         COLOR_SHADING           = 0x04,
-        COLOR_SHADING_AND_LINES = 0x05
+        COLOR_SHADING_AND_LINES = 0x05,
+        LIGHT_LINES             = 0x08,
     };
     enum TrackSpan
     {