OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-26 00:06:49 -06:00
Code Issues Packages Projects Releases Wiki Activity

#3018 Fracture Truncation : Evaluate fault throw as diff between average z value for two base edge vertices

Browse Source
This commit is contained in:
Magne Sjaastad 2018-06-07 22:10:50 +02:00
parent dae054cf8f
commit 10d0c66262
1 changed files with 52 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

69
ApplicationCode/ProjectDataModel/Completions/RimFractureContainmentTools.cpp
View File

@ -56,6 +56,39 @@ void RimFractureContainmentTools::appendNeighborCellForFace(const std::set<size_
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double computeAverageZFromTwoDeepestZ(const RigMainGrid* mainGrid,
size_t globalReservoirCellIndex,
cvf::StructGridInterface::FaceType face)
{
const RigCell& currentCell = mainGrid->globalCellArray()[globalReservoirCellIndex];
cvf::Vec3d hexCorners[8];
mainGrid->cellCornerVertices(globalReservoirCellIndex, hexCorners);
double avgZ = 0.0;
cvf::ubyte faceVertexIndices[4];
cvf::StructGridInterface::cellFaceVertexIndices(face, faceVertexIndices);
for (const auto& faceIdx : faceVertexIndices)
{
// Face indices 0-3 are defined to have deepest Z
// See void StructGridInterface::cellFaceVertexIndices(FaceType face, cvf::ubyte vertexIndices[4])
if (faceIdx < 4)
{
avgZ += hexCorners[faceIdx].z();
}
}
avgZ /= 2.0;
return avgZ;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
@ -77,29 +110,31 @@ void RimFractureContainmentTools::checkFaultAndAppendNeighborCell(const std::set
{
// See RigMainGrid::calculateFaults()
// TODO: Remove when we know if LGR can have faults
size_t neighborGlobalReservoirCellIndex = 0;
{
const RigGridBase* hostGrid = nullptr;
size_t gridLocalCellIndex;
hostGrid = mainGrid->gridAndGridLocalIdxFromGlobalCellIdx(globalReservoirCellIndex, &gridLocalCellIndex);
CVF_ASSERT(hostGrid);
const RigGridBase* hostGrid = nullptr;
size_t i, j, k;
hostGrid->ijkFromCellIndex(gridLocalCellIndex, &i, &j, &k);
size_t gridLocalCellIndex;
hostGrid = mainGrid->gridAndGridLocalIdxFromGlobalCellIdx(globalReservoirCellIndex, &gridLocalCellIndex);
size_t neighborGridLocalCellIndex;
size_t i, j, k;
hostGrid->ijkFromCellIndex(gridLocalCellIndex, &i, &j, &k);
bool foundCell = hostGrid->cellIJKNeighbor(i, j, k, face, &neighborGridLocalCellIndex);
CVF_ASSERT(foundCell);
size_t neighborGridLocalCellIndex;
neighborGlobalReservoirCellIndex =
mainGrid->reservoirCellIndexByGridAndGridLocalCellIndex(hostGrid->gridIndex(), neighborGridLocalCellIndex);
}
bool foundCell = hostGrid->cellIJKNeighbor(i, j, k, face, &neighborGridLocalCellIndex);
CVF_ASSERT(foundCell);
double currentCellAvgZ = computeAverageZFromTwoDeepestZ(mainGrid, globalReservoirCellIndex, face);
double neighborCellAvgZ = computeAverageZFromTwoDeepestZ(
mainGrid, neighborGlobalReservoirCellIndex, cvf::StructGridInterface::oppositeFace(face));
size_t neighborGlobalReservoirCellIndex = hostGrid->reservoirCellIndex(neighborGridLocalCellIndex);
const RigCell& currentCell = mainGrid->globalCellArray()[globalReservoirCellIndex];
const RigCell& neightborCell = mainGrid->globalCellArray()[neighborGlobalReservoirCellIndex];
auto diffBetweenFaceCenters =
currentCell.faceCenter(face) - neightborCell.faceCenter(cvf::StructGridInterface::oppositeFace(face));
if (diffBetweenFaceCenters.length() > maximumFaultThrow)
double faultThrow = fabs(currentCellAvgZ - neighborCellAvgZ);
if (faultThrow > maximumFaultThrow)
{
return;
}