#1557 Adding equations for peacemans formula, well bore transmissibility for a given direction, and total connection factor based on X, Y and Z trans.

ApplicationCode/ReservoirDataModel/CMakeLists_files.cmake

@@ -55,6 +55,8 @@ ${CEE_CURRENT_LIST_DIR}RigLasFileExporter.h
 ${CEE_CURRENT_LIST_DIR}RigSimulationWellCoordsAndMD.h
 ${CEE_CURRENT_LIST_DIR}RigFishbonesGeometry.h
 ${CEE_CURRENT_LIST_DIR}RigWellPathIntersectionTools.h
+${CEE_CURRENT_LIST_DIR}RigTransmissibilityEquations.h
+
 )
 
 set (SOURCE_GROUP_SOURCE_FILES
@@ -103,6 +105,8 @@ ${CEE_CURRENT_LIST_DIR}RigLasFileExporter.cpp
 ${CEE_CURRENT_LIST_DIR}RigSimulationWellCoordsAndMD.cpp
 ${CEE_CURRENT_LIST_DIR}RigFishbonesGeometry.cpp
 ${CEE_CURRENT_LIST_DIR}RigWellPathIntersectionTools.cpp
+${CEE_CURRENT_LIST_DIR}RigTransmissibilityEquations.cpp
+
 )
 
 list(APPEND CODE_HEADER_FILES

ApplicationCode/ReservoirDataModel/RigTransmissibilityEquations.cpp

@@ -0,0 +1,80 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  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 "RigTransmissibilityEquations.h"
+#include "cvfBase.h"
+#include "cvfMath.h"
+#include <cmath>
+
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+double RigTransmissibilityEquations::peacemanRadius(double permeabilityNormalDirection1, 
+                                                    double permeabilityNormalDirection2, 
+                                                    double cellSizeNormalDirection1, 
+                                                    double cellSizeNormalDirection2)
+{
+    double nominator = cvf::Math::sqrt(
+                       pow(cellSizeNormalDirection2, 2) + pow(permeabilityNormalDirection1 / permeabilityNormalDirection2, 1 / 2)
+                     + pow(cellSizeNormalDirection1, 2) + pow(permeabilityNormalDirection2 / permeabilityNormalDirection1, 1 / 2));
+
+    double denominator = pow((permeabilityNormalDirection1 / permeabilityNormalDirection2), 1 / 4)
+                       + pow((permeabilityNormalDirection2 / permeabilityNormalDirection1), 1 / 4);
+
+    double r0 = 0.28 * nominator / denominator;
+
+    return r0;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+double RigTransmissibilityEquations::wellBoreTransmissibilityComponent(double cellPerforationVectorComponent, 
+                                                                       double permeabilityNormalDirection1, 
+                                                                       double permeabilityNormalDirection2, 
+                                                                       double cellSizeNormalDirection1, 
+                                                                       double cellSizeNormalDirection2, 
+                                                                       double wellRadius,
+                                                                       double skinFactor, 
+                                                                       double cDarcyForRelevantUnit)
+{
+    double K = cvf::Math::sqrt(permeabilityNormalDirection1 * permeabilityNormalDirection2);
+
+    double nominator = cDarcyForRelevantUnit * 2 * cvf::PI_D * K * cellPerforationVectorComponent;
+    
+    double peaceManRad = peacemanRadius(permeabilityNormalDirection1, 
+                                        permeabilityNormalDirection2, 
+                                        cellSizeNormalDirection1, 
+                                        cellSizeNormalDirection2);
+
+    double denominator = log(peaceManRad / wellRadius) + skinFactor;
+
+    double trans = nominator / denominator;
+    return trans;
+}
+
+//--------------------------------------------------------------------------------------------------
+/// 
+//--------------------------------------------------------------------------------------------------
+double RigTransmissibilityEquations::totalConnectionFactor(double transX, double transY, double transZ)
+{
+    return  cvf::Math::sqrt(
+        pow(transX, 2.0) + pow(transY, 2.0) + pow(transZ, 2));
+}
+

ApplicationCode/ReservoirDataModel/RigTransmissibilityEquations.h

@@ -0,0 +1,52 @@
+/////////////////////////////////////////////////////////////////////////////////
+//
+//  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
+
+
+class RigTransmissibilityEquations
+{
+public:
+
+    // Calculations are assuming an orthogonal coordinate system
+
+    //If using wellBoreTransmissibilityComponent to calculate Tx (transmissibility in x direction), 
+    // perforationVectorComponent is the x component (in the cell local coordinate system) of the perforation vector
+    // permeability and cell size for Z and Y are to be specified as "normal directions" 1 and 2
+    // but normal directions 1 and 2 are interchangeable (so Z=1, Y=2 and Z=2, Y=1 gives same result)
+
+    static double               peacemanRadius(double permeabilityNormalDirection1,
+                                               double permeabilityNormalDirection2,
+                                               double cellSizeNormalDirection1,
+                                               double cellSizeNormalDirection2);
+
+    static double               wellBoreTransmissibilityComponent(double cellPerforationVectorComponent, 
+                                                                  double permeabilityNormalDirection1,
+                                                                  double permeabilityNormalDirection2,
+                                                                  double cellSizeNormalDirection1,
+                                                                  double cellSizeNormalDirection2,
+                                                                  double wellRadius, 
+                                                                  double skinFactor, 
+                                                                  double cDarcyForRelevantUnit);
+
+    static double               totalConnectionFactor(double transX,
+                                                      double transY,
+                                                      double transZ);
+
+};
+