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

#1091 - pre-proto - Moving transform of polygons to fracture plane out of function for calculating area.

Browse Source
This commit is contained in:
astridkbjorke
2017-01-23 13:48:19 +01:00
parent 86636de274
commit eb423dd879
1 changed files with 16 additions and 21 deletions
Download Patch File Download Diff File Expand all files Collapse all files

37
ApplicationCode/ProjectDataModel/RimFracture.cpp
View File

@@ -271,6 +271,17 @@ void RimFracture::computeTransmissibility()
bool isPlanIntersected = planeCellIntersectionPolygons(fracCell, planeCellPolygons);
if (!isPlanIntersected || planeCellPolygons.size()==0) continue;
//Transform planCell polygon(s) to x/y coordinate system (where fracturePolygon already is located)
cvf::Mat4f invertedTransMatrix = transformMatrix().getInverted();
for (std::vector<cvf::Vec3d> planeCellPolygon : planeCellPolygons)
{
for (cvf::Vec3d& v : planeCellPolygon)
{
v.transformPoint(static_cast<cvf::Mat4d>(invertedTransMatrix));
}
}
RigFractureData fracData;
fracData.reservoirCellIndex = fracCell;
@@ -280,6 +291,7 @@ void RimFracture::computeTransmissibility()
{
double areaOfCellPlaneFractureOverlap = 0.0;
std::vector<cvf::Vec3f> fracPolygon = attachedFractureDefinition()->fracturePolygon();
calculateFracturePlaneCellPolygonOverlapArea(planeCellPolygons, fracPolygon, areaOfCellPlaneFractureOverlap);
//TODO: get correct input values...
@@ -288,8 +300,7 @@ void RimFracture::computeTransmissibility()
double c = 0.008527; // TODO: Get value with units, is defined in RimReservoirCellResultsStorage
//transmissibility = areaOfCellPlaneFractureOverlap; //TODO: Only for debug - write area to file instead of transmissibility!
//TODO: Les perm fra data-strukturer...
//TODO: read permeability from file (should use matrix permeability and not fracture perm)...
transmissibility = 8 * c * attachedFractureDefinition()->permeability * areaOfCellPlaneFractureOverlap /
( flowLength + (attachedFractureDefinition()->skinFactor * fractureLength) / cvf::PI_D);
}
@@ -466,21 +477,12 @@ void RimFracture::createPolygonFromLineSegments(std::list<std::pair<cvf::Vec3d,
void RimFracture::calculateFracturePlaneCellPolygonOverlapArea(std::vector<std::vector<cvf::Vec3d> > planeCellPolygons,
std::vector<cvf::Vec3f> fracturePolygon, double & areaOfCellPlaneFractureOverlap)
{
cvf::Mat4f invertedTransMatrix = transformMatrix().getInverted();
int polygonScaleFactor = 10000; //For transform to clipper int
int xInt, yInt;
cvf::Vec3d areaVector = cvf::Vec3d::ZERO;
for (std::vector<cvf::Vec3d> planeCellPolygon : planeCellPolygons)
{
//Transform planCell polygon(s) to x/y coordinate system (where fracturePolygon already is located)
for (cvf::Vec3d& v : planeCellPolygon)
{
v.transformPoint(static_cast<cvf::Mat4d>(invertedTransMatrix));
}
//Convert to int for clipper library and store as clipper "path"
ClipperLib::Path planeCellPath;
for (cvf::Vec3d& v : planeCellPolygon)
@@ -498,12 +500,10 @@ void RimFracture::calculateFracturePlaneCellPolygonOverlapArea(std::vector<std::
fracturePath.push_back(ClipperLib::IntPoint(xInt, yInt));
}
ClipperLib::Clipper clpr;
clpr.AddPath(fracturePath, ClipperLib::ptSubject, true);
clpr.AddPath(planeCellPath, ClipperLib::ptClip, true);
ClipperLib::Paths solution;
clpr.Execute(ClipperLib::ctIntersection, solution, ClipperLib::pftEvenOdd, ClipperLib::pftEvenOdd);
@@ -524,15 +524,10 @@ void RimFracture::calculateFracturePlaneCellPolygonOverlapArea(std::vector<std::
for (std::vector<cvf::Vec3d> areaPolygon : clippedPolygons)
{
areaVector += cvf::GeometryTools::polygonAreaNormal3D(areaPolygon);
areaVector = cvf::GeometryTools::polygonAreaNormal3D(areaPolygon);
areaOfCellPlaneFractureOverlap += areaVector.length();
}
}
areaOfCellPlaneFractureOverlap = areaVector.length();
}
}
//--------------------------------------------------------------------------------------------------