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

#3657 Temporary LGR : Set transmissibility to zero for temporary grid cells in non-active cells

Browse Source
This commit is contained in:
Magne Sjaastad 2018-11-09 14:12:17 +01:00
parent 71701cd79b
commit ef17a7c057
1 changed files with 15 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
ApplicationCode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp
View File

@ -264,9 +264,23 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM
double fractureAreaWeightedlength = totalAreaXLength / fractureArea; double fractureAreaWeightedlength = totalAreaXLength / fractureArea;
// Transmissibility for inactive cells is set to zero // Transmissibility for inactive cells is set to zero
// Inactive cells must be include in order to compute the fractured area correctly
double transmissibility = 0.0; double transmissibility = 0.0;
if (activeCellInfo->isActive(reservoirCellIndex)) bool isActive = true;
{
// Use main grid cell to evaluate if a cell is active or not.
// All cells in temporary grids are active
const RigCell& cell = mainGrid->globalCellArray()[reservoirCellIndex];
size_t mainGridReservoirIndex = cell.mainGridCellIndex();
if (!activeCellInfo->isActive(mainGridReservoirIndex))
{
isActive = false;
}
}
if (isActive)
{ {
double permX = dataAccessObjectPermX->cellScalarGlobIdx(reservoirCellIndex); double permX = dataAccessObjectPermX->cellScalarGlobIdx(reservoirCellIndex);
double permY = dataAccessObjectPermY->cellScalarGlobIdx(reservoirCellIndex); double permY = dataAccessObjectPermY->cellScalarGlobIdx(reservoirCellIndex);