mirror of
https://github.com/OPM/ResInsight.git
synced 2025-02-25 18:55:39 -06:00
289 lines
10 KiB
C++
289 lines
10 KiB
C++
/////////////////////////////////////////////////////////////////////////////////
|
|
//
|
|
// Copyright (C) 2011-2012 Statoil ASA, Ceetron 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.
|
|
//
|
|
/////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#include "RIStdInclude.h"
|
|
#include "RimResultReservoir.h"
|
|
#include "RigEclipseCase.h"
|
|
#include "RifReaderEclipseOutput.h"
|
|
#include "RigReservoirCellResults.h"
|
|
#include "RimReservoirView.h"
|
|
#include "RifReaderMockModel.h"
|
|
#include "RifReaderEclipseInput.h"
|
|
#include "cafProgressInfo.h"
|
|
#include "RimProject.h"
|
|
|
|
|
|
CAF_PDM_SOURCE_INIT(RimResultReservoir, "EclipseCase");
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RimResultReservoir::RimResultReservoir()
|
|
: RimReservoir()
|
|
{
|
|
CAF_PDM_InitObject("Eclipse Case", ":/AppLogo48x48.png", "", "");
|
|
|
|
CAF_PDM_InitField(&caseFileName, "CaseFileName", QString(), "Case file name", "", "" ,"");
|
|
CAF_PDM_InitField(&caseDirectory, "CaseFolder", QString(), "Directory", "", "" ,"");
|
|
}
|
|
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
bool RimResultReservoir::openEclipseGridFile()
|
|
{
|
|
caf::ProgressInfo progInfo(50, "Reading Eclipse Grid File");
|
|
|
|
progInfo.setProgressDescription("Open Grid File");
|
|
progInfo.setNextProgressIncrement(48);
|
|
// Early exit if reservoir data is created
|
|
if (this->reservoirData()) return true;
|
|
|
|
cvf::ref<RifReaderInterface> readerInterface;
|
|
|
|
if (caseName().contains("Result Mock Debug Model"))
|
|
{
|
|
readerInterface = this->createMockModel(this->caseName());
|
|
|
|
}
|
|
else
|
|
{
|
|
QString fname = createAbsoluteFilenameFromCase(caseName);
|
|
if (fname.isEmpty())
|
|
{
|
|
return false;
|
|
}
|
|
|
|
cvf::ref<RigEclipseCase> eclipseCase = new RigEclipseCase;
|
|
readerInterface = new RifReaderEclipseOutput;
|
|
if (!readerInterface->open(fname, eclipseCase.p()))
|
|
{
|
|
return false;
|
|
}
|
|
|
|
this->setReservoirData( eclipseCase.p() );
|
|
}
|
|
|
|
results(RifReaderInterface::MATRIX_RESULTS)->setReaderInterface(readerInterface.p());
|
|
results(RifReaderInterface::FRACTURE_RESULTS)->setReaderInterface(readerInterface.p());
|
|
|
|
progInfo.incrementProgress();
|
|
|
|
CVF_ASSERT(this->reservoirData());
|
|
CVF_ASSERT(readerInterface.notNull());
|
|
|
|
progInfo.setProgressDescription("Computing Case Cache");
|
|
computeCachedData();
|
|
|
|
return true;
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
bool RimResultReservoir::openAndReadActiveCellData(RigEclipseCase* mainEclipseCase)
|
|
{
|
|
cvf::ref<RifReaderInterface> readerInterface;
|
|
|
|
if (caseName().contains("Result Mock Debug Model"))
|
|
{
|
|
readerInterface = this->createMockModel(this->caseName());
|
|
}
|
|
else
|
|
{
|
|
QString fname = createAbsoluteFilenameFromCase(caseName);
|
|
if (fname.isEmpty())
|
|
{
|
|
return false;
|
|
}
|
|
|
|
cvf::ref<RigEclipseCase> eclipseCase = new RigEclipseCase;
|
|
|
|
CVF_ASSERT(mainEclipseCase && mainEclipseCase->mainGrid());
|
|
eclipseCase->setMainGrid(mainEclipseCase->mainGrid());
|
|
|
|
size_t scalarIndexWithMaxTimeStepCount = cvf::UNDEFINED_SIZE_T;
|
|
mainEclipseCase->results(RifReaderInterface::MATRIX_RESULTS)->maxTimeStepCount(&scalarIndexWithMaxTimeStepCount);
|
|
if (scalarIndexWithMaxTimeStepCount == cvf::UNDEFINED_SIZE_T)
|
|
{
|
|
return false;
|
|
}
|
|
|
|
std::vector<QDateTime> timeStepDates = mainEclipseCase->results(RifReaderInterface::MATRIX_RESULTS)->timeStepDates(scalarIndexWithMaxTimeStepCount);
|
|
|
|
cvf::ref<RifReaderEclipseOutput> readerEclipseOutput = new RifReaderEclipseOutput;
|
|
if (!readerEclipseOutput->openAndReadActiveCellData(fname, timeStepDates, eclipseCase.p()))
|
|
{
|
|
return false;
|
|
}
|
|
|
|
readerEclipseOutput->close();
|
|
|
|
this->setReservoirData( eclipseCase.p() );
|
|
|
|
readerInterface = readerEclipseOutput;
|
|
}
|
|
|
|
results(RifReaderInterface::MATRIX_RESULTS)->setReaderInterface(readerInterface.p());
|
|
results(RifReaderInterface::FRACTURE_RESULTS)->setReaderInterface(readerInterface.p());
|
|
|
|
CVF_ASSERT(this->reservoirData());
|
|
CVF_ASSERT(readerInterface.notNull());
|
|
|
|
reservoirData()->computeCachedData();
|
|
|
|
return true;
|
|
}
|
|
|
|
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
cvf::ref<RifReaderInterface> RimResultReservoir::createMockModel(QString modelName)
|
|
{
|
|
cvf::ref<RifReaderMockModel> mockFileInterface = new RifReaderMockModel;
|
|
cvf::ref<RigEclipseCase> reservoir = new RigEclipseCase;
|
|
|
|
if (modelName == "Result Mock Debug Model Simple")
|
|
{
|
|
// Create the mock file interface and and RigSerervoir and set them up.
|
|
mockFileInterface->setWorldCoordinates(cvf::Vec3d(10, 10, 10), cvf::Vec3d(20, 20, 20));
|
|
mockFileInterface->setGridPointDimensions(cvf::Vec3st(4, 5, 6));
|
|
mockFileInterface->addLocalGridRefinement(cvf::Vec3st(0, 2, 2), cvf::Vec3st(0, 2, 2), cvf::Vec3st(3, 3, 3));
|
|
|
|
mockFileInterface->open("", reservoir.p());
|
|
{
|
|
size_t idx = reservoir->mainGrid()->cellIndexFromIJK(1, 3, 4);
|
|
|
|
//TODO: Rewrite active cell info in mock models
|
|
//reservoir->mainGrid()->cell(idx).setActiveIndexInMatrixModel(cvf::UNDEFINED_SIZE_T);
|
|
}
|
|
|
|
{
|
|
size_t idx = reservoir->mainGrid()->cellIndexFromIJK(2, 2, 3);
|
|
|
|
//TODO: Rewrite active cell info in mock models
|
|
//reservoir->mainGrid()->cell(idx).setActiveIndexInMatrixModel(cvf::UNDEFINED_SIZE_T);
|
|
}
|
|
}
|
|
else if (modelName == "Result Mock Debug Model With Results")
|
|
{
|
|
mockFileInterface->setWorldCoordinates(cvf::Vec3d(10, 10, 10), cvf::Vec3d(-20, -20, -20));
|
|
mockFileInterface->setGridPointDimensions(cvf::Vec3st(5, 10, 20));
|
|
mockFileInterface->addLocalGridRefinement(cvf::Vec3st(0, 3, 3), cvf::Vec3st(1, 4, 9), cvf::Vec3st(2, 2, 2));
|
|
mockFileInterface->setResultInfo(3, 10);
|
|
|
|
mockFileInterface->open("", reservoir.p());
|
|
|
|
// Make a fault
|
|
cvf::Vec3d& tmp = reservoir->mainGrid()->nodes()[1];
|
|
tmp += cvf::Vec3d(1, 0, 0);
|
|
}
|
|
else if (modelName =="Result Mock Debug Model Large With Results")
|
|
{
|
|
double startX = 0;
|
|
double startY = 0;
|
|
double startZ = 0;
|
|
|
|
double widthX = 6000;
|
|
double widthY = 12000;
|
|
double widthZ = 500;
|
|
|
|
double offsetX = 0;
|
|
double offsetY = 0;
|
|
double offsetZ = 0;
|
|
|
|
// Test code to simulate UTM coordinates
|
|
offsetX = 400000;
|
|
offsetY = 6000000;
|
|
offsetZ = 0;
|
|
|
|
|
|
mockFileInterface->setWorldCoordinates(cvf::Vec3d(startX + offsetX, startY + offsetY, startZ + offsetZ), cvf::Vec3d(startX + widthX + offsetX, startY + widthY + offsetY, startZ + widthZ + offsetZ));
|
|
mockFileInterface->setGridPointDimensions(cvf::Vec3st(50, 100, 200));
|
|
mockFileInterface->addLocalGridRefinement(cvf::Vec3st(0, 30, 30), cvf::Vec3st(1, 40, 90), cvf::Vec3st(2, 2, 2));
|
|
mockFileInterface->setResultInfo(3, 10);
|
|
|
|
mockFileInterface->open("", reservoir.p());
|
|
|
|
}
|
|
|
|
this->setReservoirData( reservoir.p() );
|
|
|
|
return mockFileInterface.p();
|
|
}
|
|
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
RimResultReservoir::~RimResultReservoir()
|
|
{
|
|
reservoirViews.deleteAllChildObjects();
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
QString RimResultReservoir::locationOnDisc() const
|
|
{
|
|
return caseDirectory;
|
|
}
|
|
|
|
//--------------------------------------------------------------------------------------------------
|
|
///
|
|
//--------------------------------------------------------------------------------------------------
|
|
QString RimResultReservoir::createAbsoluteFilenameFromCase(const QString& caseName)
|
|
{
|
|
QString candidate;
|
|
|
|
candidate = QDir::fromNativeSeparators(caseDirectory.v() + QDir::separator() + caseName + ".EGRID");
|
|
if (QFile::exists(candidate)) return candidate;
|
|
|
|
candidate = QDir::fromNativeSeparators(caseDirectory.v() + QDir::separator() + caseName + ".GRID");
|
|
if (QFile::exists(candidate)) return candidate;
|
|
|
|
std::vector<caf::PdmObject*> parentObjects;
|
|
this->parentObjects(parentObjects);
|
|
|
|
QString projectPath;
|
|
for (size_t i = 0; i < parentObjects.size(); i++)
|
|
{
|
|
caf::PdmObject* obj = parentObjects[i];
|
|
RimProject* proj = dynamic_cast<RimProject*>(obj);
|
|
if (proj)
|
|
{
|
|
QFileInfo fi(proj->fileName);
|
|
projectPath = fi.path();
|
|
}
|
|
}
|
|
|
|
if (!projectPath.isEmpty())
|
|
{
|
|
candidate = QDir::fromNativeSeparators(projectPath + QDir::separator() + caseName + ".EGRID");
|
|
if (QFile::exists(candidate)) return candidate;
|
|
|
|
candidate = QDir::fromNativeSeparators(projectPath + QDir::separator() + caseName + ".GRID");
|
|
if (QFile::exists(candidate)) return candidate;
|
|
}
|
|
|
|
return QString();
|
|
}
|
|
|