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

Roff: import properties from Roff grid files.

Browse Source
This commit is contained in:
Kristian Bendiksen
2022-12-19 15:56:50 +01:00
parent 492f2dc4f6
commit ca5064e7da
4 changed files with 130 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

116
ApplicationLibCode/FileInterface/RifRoffFileTools.cpp
View File

@@ -486,6 +486,34 @@ void RifRoffFileTools::convertToReservoirIndexOrder( int nx
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RifRoffFileTools::convertToReservoirIndexOrder( int nx,
int ny,
int nz,
const std::vector<float>& in,
std::vector<float>& out )
{
CAF_ASSERT( static_cast<size_t>( nx ) * ny * nz == in.size() );
out.resize( in.size(), 0.0 );
int outIdx = 0;
for ( int k = 0; k < nz; k++ )
{
for ( int j = 0; j < ny; j++ )
{
for ( int i = 0; i < nx; i++ )
{
int inIdx = i * ny * nz + j * nz + ( nz - k - 1 );
out[outIdx] = in[inIdx];
outIdx++;
}
}
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
@@ -508,3 +536,91 @@ size_t RifRoffFileTools::computeActiveCellMatrixIndex( std::vector<int>& activeC
return activeMatrixIndex;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
bool RifRoffFileTools::createInputProperties( const QString& fileName, RigEclipseCaseData* eclipseCaseData )
{
RiaLogging::info( QString( "Opening roff file: %1" ).arg( fileName ) );
std::string filename = fileName.toStdString();
std::ifstream stream( filename, std::ios::binary );
if ( !stream.good() )
{
RiaLogging::error( "Unable to open roff file" );
return false;
}
try
{
Reader reader( stream );
reader.parse();
std::vector<std::pair<std::string, Token::Kind>> arrayTypes = reader.getNamedArrayTypes();
for ( auto [keyword, kind] : arrayTypes )
{
size_t keywordLength = reader.getArrayLength( keyword );
RiaLogging::info( QString( "Array found: %1 . Type: %2 Size: %3" )
.arg( QString::fromStdString( keyword ) )
.arg( QString::fromStdString( Token::kindToString( kind ) ) )
.arg( keywordLength ) );
if ( eclipseCaseData->mainGrid()->cellCount() == keywordLength )
{
QString newResultName = eclipseCaseData->results( RiaDefines::PorosityModelType::MATRIX_MODEL )
->makeResultNameUnique( QString::fromStdString( keyword ) );
if ( !appendNewInputPropertyResult( eclipseCaseData, newResultName, keyword, reader ) )
{
RiaLogging::error( QString( "Unable to import result '%1' from %2" )
.arg( QString::fromStdString( keyword ) )
.arg( fileName ) );
return false;
}
}
}
}
catch ( std::runtime_error& err )
{
RiaLogging::error( QString( "Roff property file import failed: %1" ).arg( err.what() ) );
return false;
}
return true;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
bool RifRoffFileTools::appendNewInputPropertyResult( RigEclipseCaseData* caseData,
const QString& resultName,
const std::string& keyword,
Reader& reader )
{
CVF_ASSERT( caseData );
std::vector<float> values = reader.getFloatArray( keyword );
CAF_ASSERT( values.size() == caseData->mainGrid()->cellCount() );
int nx = caseData->mainGrid()->cellCountI();
int ny = caseData->mainGrid()->cellCountJ();
int nz = caseData->mainGrid()->cellCountK();
std::vector<float> valuesReservoirOrder;
convertToReservoirIndexOrder( nx, ny, nz, values, valuesReservoirOrder );
RigEclipseResultAddress resAddr( RiaDefines::ResultCatType::INPUT_PROPERTY, resultName );
caseData->results( RiaDefines::PorosityModelType::MATRIX_MODEL )->createResultEntry( resAddr, false );
auto newPropertyData =
caseData->results( RiaDefines::PorosityModelType::MATRIX_MODEL )->modifiableCellScalarResultTimesteps( resAddr );
std::vector<double> doubleVals;
doubleVals.insert( doubleVals.begin(), valuesReservoirOrder.begin(), valuesReservoirOrder.end() );
newPropertyData->push_back( doubleVals );
return true;
}

11
ApplicationLibCode/FileInterface/RifRoffFileTools.h
View File

@@ -27,6 +27,7 @@
class RigEclipseCaseData;
class RigMainGrid;
class Reader;
//==================================================================================================
//
@@ -41,6 +42,8 @@ public:
static bool openGridFile( const QString& fileName, RigEclipseCaseData* eclipseCase, QString* errorMessages );
static bool createInputProperties( const QString& fileName, RigEclipseCaseData* eclipseCase );
private:
static void interpretSplitenzData( int nz,
float zoffset,
@@ -52,6 +55,9 @@ private:
static void
convertToReservoirIndexOrder( int nx, int ny, int nz, const std::vector<char>& activeIn, std::vector<int>& activeOut );
static void
convertToReservoirIndexOrder( int nx, int ny, int nz, const std::vector<float>& in, std::vector<float>& out );
static size_t computeActiveCellMatrixIndex( std::vector<int>& activeCells );
static cvf::Vec3d getCorner( const RigMainGrid& grid,
@@ -63,4 +69,9 @@ private:
const cvf::Vec3d& scale );
static double interpolate( const cvf::Vec3d& top, const cvf::Vec3d& bottom, double z, int idx );
static bool appendNewInputPropertyResult( RigEclipseCaseData* caseData,
const QString& resultName,
const std::string& keyword,
Reader& reader );
};