mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
57a049d39f
Create stored native riTRANX,Y,Z results instead of doing everyting in one go. This makes the statistics calculation work automatically instead of beeing custoized and require a compplete calculation of its own.
71 lines
3.0 KiB
C++
71 lines
3.0 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) Statoil ASA, Ceetron Solutions AS
|
|
//
|
|
// ResInsight is free software: you can redistribute it and/or modify
|
|
// it under the terms of the GNU General Public License as published by
|
|
// the Free Software Foundation, either version 3 of the License, or
|
|
// (at your option) any later version.
|
|
//
|
|
// ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
|
|
// WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
// FITNESS FOR A PARTICULAR PURPOSE.
|
|
//
|
|
// See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
|
|
// for more details.
|
|
//
|
|
/////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#pragma once
|
|
|
|
#include "RigResultAccessor.h"
|
|
|
|
#include "cvfCollection.h"
|
|
|
|
class RigGridBase;
|
|
|
|
|
|
//==================================================================================================
|
|
///
|
|
//==================================================================================================
|
|
class RigCombRiTransResultAccessor : public RigResultAccessor
|
|
{
|
|
public:
|
|
RigCombRiTransResultAccessor(const RigGridBase* grid);
|
|
|
|
void setPermResultAccessors(RigResultAccessor* xPermAccessor,
|
|
RigResultAccessor* yPermAccessor,
|
|
RigResultAccessor* zPermAccessor);
|
|
void setNTGResultAccessor( RigResultAccessor* ntgAccessor);
|
|
|
|
/// CDARCY Darcy's constant
|
|
/// = 0.00852702 (E300); 0.008527 (ECLIPSE 100) (METRIC)
|
|
/// = 0.00112712 (E300); 0.001127 (ECLIPSE 100) (FIELD)
|
|
/// = 3.6 (LAB)
|
|
/// = 0.00864 (PVT-M)
|
|
void setCDARCHY(double cDarchy) { m_cdarchy = cDarchy;}
|
|
|
|
virtual double cellScalar(size_t gridLocalCellIndex) const;
|
|
virtual double cellFaceScalar(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const;
|
|
|
|
private:
|
|
double getPermValue(size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId) const;
|
|
double calculateHalfCellTrans( size_t gridLocalCellIndex, size_t neighborGridCellIdx, cvf::StructGridInterface::FaceType faceId, bool isFaultFace) const;
|
|
|
|
double getNtgValue( size_t gridLocalCellIndex, cvf::StructGridInterface::FaceType faceId ) const;
|
|
static void calculateConnectionGeometry( const RigGridBase* grid, size_t gridLocalCellIndex, size_t neighborGridCellIdx, cvf::StructGridInterface::FaceType faceId, cvf::Vec3d* centerToFace, cvf::Vec3d* faceAreaVec);
|
|
|
|
static double halfCellTransmissibility(double perm, double ntg, const cvf::Vec3d& centerToFace, const cvf::Vec3d& faceAreaVec);
|
|
static double newtran(double cdarchy, double mult, double halfCellTrans, double neighborHalfCellTrans);
|
|
|
|
cvf::ref<RigResultAccessor> m_xPermAccessor;
|
|
cvf::ref<RigResultAccessor> m_yPermAccessor;
|
|
cvf::ref<RigResultAccessor> m_zPermAccessor;
|
|
|
|
cvf::ref<RigResultAccessor> m_ntgAccessor;
|
|
|
|
double m_cdarchy;
|
|
|
|
const RigGridBase* m_grid;
|
|
};
|