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

#6651: Allow wild card in formation name when importing elastic properties.

Browse Source
This commit is contained in:
Kristian Bendiksen 2020-10-07 12:41:10 +02:00
parent 5e31df53bc
commit 6c8dda3869
2 changed files with 86 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

103
ApplicationCode/Commands/RicElasticPropertiesImportTools.cpp
View File

@ -20,15 +20,20 @@
#include "RiaLogging.h"
#include "RimElasticProperties.h"
#include "RimFractureModelTemplate.h"
#include "RifCsvUserDataParser.h"
#include "RifElasticPropertiesReader.h"
#include "RifFileParseTools.h"
#include "RigElasticProperties.h"
#include "RigFormationNames.h"
#include "RimElasticProperties.h"
#include "RimFormationNames.h"
#include "RimFormationNamesCollection.h"
#include "RimFractureModelTemplate.h"
#include "RimOilField.h"
#include "RimProject.h"
#include <set>
#include <vector>
@ -36,7 +41,8 @@
///
//--------------------------------------------------------------------------------------------------
void RicElasticPropertiesImportTools::importElasticPropertiesFromFile( const QString& filePath,
RimFractureModelTemplate* fractureModelTemplate )
RimFractureModelTemplate* fractureModelTemplate,
const QString& formationWildCard )
{
RifCsvUserDataFileParser csvParser( filePath );
QString separator = csvParser.tryDetermineCellSeparator();
@ -57,12 +63,26 @@ void RicElasticPropertiesImportTools::importElasticPropertiesFromFile( const QSt
return;
}
std::vector<QString> formationNames = getFormationNames();
// Find the unique facies keys (combination of field, formation and facies names)
std::set<FaciesKey> faciesKeys;
for ( RifElasticProperties item : rifElasticProperties )
{
FaciesKey faciesKey = std::make_tuple( item.fieldName, item.formationName, item.faciesName );
faciesKeys.insert( faciesKey );
if ( item.formationName == formationWildCard )
{
// Found wildcard: generate keys for all formations
for ( QString formationName : formationNames )
{
FaciesKey faciesKey = std::make_tuple( item.fieldName, formationName, item.faciesName );
faciesKeys.insert( faciesKey );
}
}
else
{
FaciesKey faciesKey = std::make_tuple( item.fieldName, item.formationName, item.faciesName );
faciesKeys.insert( faciesKey );
}
}
RimElasticProperties* rimElasticProperties = new RimElasticProperties;
@ -77,7 +97,9 @@ void RicElasticPropertiesImportTools::importElasticPropertiesFromFile( const QSt
// Group the items with a given facies key
for ( RifElasticProperties item : rifElasticProperties )
{
if ( item.fieldName == fieldName && item.formationName == formationName && item.faciesName == faciesName )
if ( item.fieldName == fieldName &&
( item.formationName == formationName || item.formationName == formationWildCard ) &&
item.faciesName == faciesName )
{
matchingFacies.push_back( item );
}
@ -88,28 +110,63 @@ void RicElasticPropertiesImportTools::importElasticPropertiesFromFile( const QSt
matchingFacies.end(),
[]( const RifElasticProperties& a, const RifElasticProperties& b ) { return a.porosity < b.porosity; } );
// Finally add the values
RigElasticProperties rigElasticProperties( fieldName, formationName, faciesName );
for ( RifElasticProperties item : matchingFacies )
std::vector<QString> matchingFormations;
if ( formationName == formationWildCard )
{
rigElasticProperties.appendValues( item.porosity,
item.youngsModulus,
item.poissonsRatio,
item.K_Ic,
item.proppantEmbedment,
item.biotCoefficient,
item.k0,
item.fluidLossCoefficient,
item.spurtLoss,
item.immobileFluidSaturation );
// Duplicate values for all formations when encountering a wild card
matchingFormations = formationNames;
}
else
{
matchingFormations.push_back( formationName );
}
// Avoid using the field name in the match for now
FaciesKey noFieldKey = std::make_tuple( "", formationName, faciesName );
rimElasticProperties->setPropertiesForFacies( noFieldKey, rigElasticProperties );
for ( QString matchingFormationName : matchingFormations )
{
// Finally add the values
RigElasticProperties rigElasticProperties( fieldName, matchingFormationName, faciesName );
for ( RifElasticProperties item : matchingFacies )
{
rigElasticProperties.appendValues( item.porosity,
item.youngsModulus,
item.poissonsRatio,
item.K_Ic,
item.proppantEmbedment,
item.biotCoefficient,
item.k0,
item.fluidLossCoefficient,
item.spurtLoss,
item.immobileFluidSaturation );
}
// Avoid using the field name in the match for now
FaciesKey noFieldKey = std::make_tuple( "", matchingFormationName, faciesName );
rimElasticProperties->setPropertiesForFacies( noFieldKey, rigElasticProperties );
}
}
rimElasticProperties->setFilePath( filePath );
fractureModelTemplate->setElasticProperties( rimElasticProperties );
fractureModelTemplate->updateConnectedEditors();
}
//--------------------------------------------------------------------------------------------------
/// Finds first formation names.
//--------------------------------------------------------------------------------------------------
std::vector<QString> RicElasticPropertiesImportTools::getFormationNames()
{
RimProject* project = RimProject::current();
RimOilField* oilField = project->activeOilField();
RimFormationNamesCollection* formationNamesCollection = oilField->formationNamesCollection();
for ( RimFormationNames* formationNames : formationNamesCollection->formationNamesList().childObjects() )
{
if ( formationNames && formationNames->formationNamesData() )
{
RigFormationNames* rigFormationNames = formationNames->formationNamesData();
if ( rigFormationNames ) return rigFormationNames->formationNames();
}
}
return std::vector<QString>();
}

7
ApplicationCode/Commands/RicElasticPropertiesImportTools.h
View File

@ -19,6 +19,7 @@
#pragma once
#include <QString>
#include <vector>
class RimFractureModelTemplate;
@ -28,9 +29,13 @@ class RimFractureModelTemplate;
class RicElasticPropertiesImportTools
{
public:
static void importElasticPropertiesFromFile( const QString& filePath, RimFractureModelTemplate* fractureModelTemplate );
static void importElasticPropertiesFromFile( const QString& filePath,
RimFractureModelTemplate* fractureModelTemplate,
const QString& formationWildCard = "1*" );
private:
// Hidden to avoid instantiation
RicElasticPropertiesImportTools();
static std::vector<QString> getFormationNames();
};