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

Use average NTG values for merged cells due to MINPV algorithm

Browse Source 
The volume weighted averaged NTG values are used to get more
reasonable x and y transmissibilities for the cells effected by the
MINPV grid processing algorithm.

This commit starts the process of supporting the MINPV-FILL keyword
where upscaled properties are used for the merged/filled cells.
This commit is contained in:
Tor Harald Sandve 2015-03-19 08:08:07 +01:00
parent 531bc7fd73
commit 65199735e2
1 changed files with 44 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

44
opm/autodiff/GeoProps.hpp
View File

@ -109,6 +109,12 @@ namespace Opm
pvol_[cellIdx] *= eclgrid->getCellVolume(cartesianCellIdx);
}
}
// Use volume weighted arithmetic average of the NTG values for
// the cells effected by the current OPM cpgrid process algorithm
// for MINPV. Note that the change does not effect the pore volume calculations
// as the pore volume is currently defaulted to be comparable to ECLIPSE, but
// only the transmissibility calculations.
minPvFillProps_(grid, eclState,ntg);
// Transmissibility
@ -186,6 +192,11 @@ namespace Opm
const double* perm,
Vector &hTrans);
template <class Grid>
void minPvFillProps_(const Grid &grid,
Opm::EclipseStateConstPtr eclState,
std::vector<double> &ntg);
Vector pvol_ ;
Vector trans_;
Vector gpot_ ;
@ -197,6 +208,39 @@ namespace Opm
};
template <class GridType>
inline void DerivedGeology::minPvFillProps_(const GridType &grid,
Opm::EclipseStateConstPtr eclState,
std::vector<double> &ntg)
{
int numCells = Opm::AutoDiffGrid::numCells(grid);
const int* global_cell = Opm::UgGridHelpers::globalCell(grid);
EclipseGridConstPtr eclgrid = eclState->getEclipseGrid();
std::vector<double> porv = eclState->getDoubleGridProperty("PORV")->getData();
for (int cellIdx = 0; cellIdx < numCells; ++cellIdx) {
const int nx = grid.cartdims[0];
const int ny = grid.cartdims[1];
const int cartesianCellIdx = global_cell[cellIdx];
const int cartesianCellIdxAbove = cartesianCellIdx - nx*ny;
// Average properties when the cell above exist
// and has porv less than the MINPV threshold
if ( cartesianCellIdxAbove > 0 &&
porv[cartesianCellIdxAbove] > 0 &&
porv[cartesianCellIdxAbove] < eclgrid->getMinpvValue() ) {
// Volume weighted arithmetic average of NTG
const double cellVolume = eclgrid->getCellVolume(cartesianCellIdx);
const double cellAboveVolume = eclgrid->getCellVolume(cartesianCellIdxAbove);
const double inv_totalCellVolume = 1/(cellVolume + cellAboveVolume);
ntg[cartesianCellIdx] *= cellAboveVolume;
ntg[cartesianCellIdx] += ntg[cartesianCellIdxAbove]*cellAboveVolume;
ntg[cartesianCellIdx] *= inv_totalCellVolume;
}
}
}
template <class GridType>
inline void DerivedGeology::multiplyHalfIntersections_(const GridType &grid,