Added helper class for measured depth

ApplicationCode/ReservoirDataModel/CMakeLists_files.cmake

@@ -49,6 +49,7 @@ ${CEE_CURRENT_LIST_DIR}RigTimeHistoryResultAccessor.h
 ${CEE_CURRENT_LIST_DIR}RigCurveDataTools.h
 ${CEE_CURRENT_LIST_DIR}RigSummaryCaseData.h
 ${CEE_CURRENT_LIST_DIR}RigLasFileExporter.h
+${CEE_CURRENT_LIST_DIR}RigSimulationWellCoordsAndMD.h
 )
 
 set (SOURCE_GROUP_SOURCE_FILES
@@ -92,6 +93,7 @@ ${CEE_CURRENT_LIST_DIR}RigTimeHistoryResultAccessor.cpp
 ${CEE_CURRENT_LIST_DIR}RigCurveDataTools.cpp    
 ${CEE_CURRENT_LIST_DIR}RigSummaryCaseData.cpp
 ${CEE_CURRENT_LIST_DIR}RigLasFileExporter.cpp
+${CEE_CURRENT_LIST_DIR}RigSimulationWellCoordsAndMD.cpp
 )
 
 list(APPEND CODE_HEADER_FILES

ApplicationCode/ReservoirDataModel/RigSimulationWellCoordsAndMD.cpp

@@ -0,0 +1,104 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  Copyright (C) 2017     Statoil 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 "RigSimulationWellCoordsAndMD.h"
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+RigSimulationWellCoordsAndMD::RigSimulationWellCoordsAndMD(const std::vector<cvf::Vec3d>& wellPathPoints)
+{
+    m_wellPathPoints = wellPathPoints;
+
+    computeMeasuredDepths();
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+const std::vector<cvf::Vec3d>& RigSimulationWellCoordsAndMD::wellPathPoints() const
+{
+    return m_wellPathPoints;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+const std::vector<double>& RigSimulationWellCoordsAndMD::measuredDepths() const
+{
+    return m_measuredDepths;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+cvf::Vec3d RigSimulationWellCoordsAndMD::interpolatedPointAlongWellPath(double measuredDepth) const
+{
+    cvf::Vec3d wellPathPoint = cvf::Vec3d::ZERO;
+
+    size_t i = 0;
+    while (i < m_measuredDepths.size() && m_measuredDepths.at(i) < measuredDepth)
+    {
+        i++;
+    }
+
+    if (m_measuredDepths.size() > i)
+    {
+        if (i == 0)
+        {
+            //For measuredDepth same or lower than first point, use this first point
+            wellPathPoint = m_wellPathPoints.at(0);
+        }
+        else
+        {
+            //Do interpolation
+            double stepsize = (measuredDepth - m_measuredDepths.at(i - 1)) /
+                (m_measuredDepths.at(i) - m_measuredDepths.at(i - 1));
+            wellPathPoint = m_wellPathPoints.at(i - 1) + stepsize * (m_wellPathPoints.at(i) - m_wellPathPoints.at(i - 1));
+        }
+    }
+    else
+    {
+        //Use endpoint if measuredDepth same or higher than last point
+        wellPathPoint = m_wellPathPoints.at(i - 1);
+    }
+
+    return wellPathPoint;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+void RigSimulationWellCoordsAndMD::computeMeasuredDepths()
+{
+    cvf::Vec3d prev = cvf::Vec3d::UNDEFINED;
+
+    double accumulatedMD = 0;
+
+    for (const auto& point : m_wellPathPoints)
+    {
+        if (!prev.isUndefined())
+        {
+            accumulatedMD += point.pointDistance(prev);
+        }
+
+        m_measuredDepths.push_back(accumulatedMD);
+
+        prev = point;
+    }
+}

ApplicationCode/ReservoirDataModel/RigSimulationWellCoordsAndMD.h

@@ -0,0 +1,49 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  Copyright (C) 2017     Statoil 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 "cvfBase.h"
+#include "cvfVector3.h"
+
+#include <vector>
+
+
+
+
+//==================================================================================================
+///  
+///  
+//==================================================================================================
+class RigSimulationWellCoordsAndMD
+{
+public:
+    RigSimulationWellCoordsAndMD(const std::vector<cvf::Vec3d>& wellPathPoints);
+
+    const std::vector<cvf::Vec3d>&  wellPathPoints() const;
+    const std::vector<double>&      measuredDepths() const;
+
+    cvf::Vec3d                      interpolatedPointAlongWellPath(double measuredDepth) const;
+
+private:
+    void                            computeMeasuredDepths();
+
+private:
+    std::vector<cvf::Vec3d> m_wellPathPoints;
+    std::vector<double>     m_measuredDepths;
+};