OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#1901 Fix error in peaceman radius calculation. Used int values 1/2 and 1/4 as exponent, which ends in 0.

Browse Source
This commit is contained in:
Jacob Støren 2017-09-20 11:02:19 +02:00
parent 8334c36eff
commit 50a88dff4a
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ApplicationCode/ReservoirDataModel/RigTransmissibilityEquations.cpp
View File

@ -33,11 +33,11 @@ double RigTransmissibilityEquations::peacemanRadius(double permeabilityNormalDir
double cellSizeNormalDirection2)
{
double numerator = cvf::Math::sqrt(
pow(cellSizeNormalDirection2, 2) * pow(permeabilityNormalDirection1 / permeabilityNormalDirection2, 1 / 2)
+ pow(cellSizeNormalDirection1, 2) * pow(permeabilityNormalDirection2 / permeabilityNormalDirection1, 1 / 2));
pow(cellSizeNormalDirection2, 2.0) * pow(permeabilityNormalDirection1 / permeabilityNormalDirection2, 0.5)
+ pow(cellSizeNormalDirection1, 2.0) * pow(permeabilityNormalDirection2 / permeabilityNormalDirection1, 0.5) );
double denominator = pow((permeabilityNormalDirection1 / permeabilityNormalDirection2), 1 / 4)
+ pow((permeabilityNormalDirection2 / permeabilityNormalDirection1), 1 / 4);
double denominator = pow((permeabilityNormalDirection1 / permeabilityNormalDirection2), 0.25 )
+ pow((permeabilityNormalDirection2 / permeabilityNormalDirection1), 0.25 );
double r0 = 0.28 * numerator / denominator;
@ -77,6 +77,6 @@ double RigTransmissibilityEquations::wellBoreTransmissibilityComponent(double ce
double RigTransmissibilityEquations::totalConnectionFactor(double transX, double transY, double transZ)
{
return cvf::Math::sqrt(
pow(transX, 2.0) + pow(transY, 2.0) + pow(transZ, 2));
pow(transX, 2.0) + pow(transY, 2.0) + pow(transZ, 2.0));
}