#1792 Avoiding issue with negative transmissibility from fracture to well due to small discretization by forcing peaceman radius r0 > wellRadius.

This commit is contained in:
astridkbjorke
2017-08-24 14:01:07 +02:00
parent c7dcf9f4d0
commit e428bf33f9

View File

@@ -66,9 +66,15 @@ double RigFractureTransmissibilityEquations::fractureCellToWellRadialTrans(doubl
double ro = 0.14 * cvf::Math::sqrt(
pow(fractureCellSizeX, 2.0) + pow(fractureCellSizeZ, 2));
if (ro < (wellRadius * 1.01))
{
ro = wellRadius * 1.01;
}
double Tc = 2 * cvf::PI_D * cDarcyForRelevantUnit * fractureCellConductivity /
(log(ro / wellRadius) + skinFactor );
CVF_TIGHT_ASSERT(Tc > 0);
return Tc;
}