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

#1141 - pre-proto - Adding intermediate values from transmissibility calculation to the top of the COMPDAT export file.

Browse Source
This commit is contained in:
astridkbjorke
2017-01-30 14:47:56 +01:00
parent b13ca8bb1c
commit 07f5ea45d2
4 changed files with 189 additions and 119 deletions
Download Patch File Download Diff File Expand all files Collapse all files

275
ApplicationCode/FileInterface/RifEclipseExportTools.cpp
View File

@@ -32,9 +32,9 @@
#include "cafProgressInfo.h"
#include <QTextStream>
#include <QFile>
#include <QString>
#include <QTextStream>
@@ -88,37 +88,10 @@ bool RifEclipseExportTools::writeFracturesToTextFile(const QString& fileName, c
out << "\n";
out << "-- Exported from ResInsight" << "\n\n";
out << "-- Background data for calculation" << "\n\n";
//Write header line
out << qSetFieldWidth(4);
out << "--";
out << qSetFieldWidth(12);
out << "Well name";
out << qSetFieldWidth(16);
out << "Fracture Name";
out << qSetFieldWidth(5);
out << "i";
out << "j";
out << "k";
out << qSetFieldWidth(8);
out << "Ax";
out << " ";
out << "Ay";
out << " ";
out << "Az";
out << " ";
out << "TotalArea";
out << "\n";
printBackgroundDataHeaderLine(out);
//Write background numbers (with -- first to comment out line)
for (RimFracture* fracture : fractures)
{
fracture->computeTransmissibility();
@@ -126,46 +99,12 @@ bool RifEclipseExportTools::writeFracturesToTextFile(const QString& fileName, c
for (RigFractureData fracData : fracDataVector)
{
out << qSetFieldWidth(4);
out << "--";
out << qSetFieldWidth(12);
wellPath, simWell = nullptr;
fracture->firstAncestorOrThisOfType(simWell);
if (simWell) out << simWell->name; // 1. Well name
fracture->firstAncestorOrThisOfType(wellPath);
if (wellPath) out << wellPath->name; // 1. Well name
out << qSetFieldWidth(16);
out << fracture->name();
out << qSetFieldWidth(5);
size_t i, j, k;
mainGrid->ijkFromCellIndex(fracData.reservoirCellIndex, &i, &j, &k);
out << i + 1; // 2. I location grid block, adding 1 to go to eclipse 1-based grid definition
out << j + 1; // 3. J location grid block, adding 1 to go to eclipse 1-based grid definition
out << k + 1; // 4. K location of upper connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << " ";
out << fracData.projectedAreas.x();
out << " ";
out << fracData.projectedAreas.y();
out << " ";
out << fracData.projectedAreas.z();
out << " ";
out << fracData.totalArea;
out << "\n";
printBackgroundData(out, wellPath, simWell, fracture, mainGrid, fracData);
}
}
out << "\n";
//Write COMPDAT keyword and entries
out << "COMPDAT" << "\n" << right << qSetFieldWidth(8);
for (RimFracture* fracture : fractures)
{
@@ -174,53 +113,9 @@ bool RifEclipseExportTools::writeFracturesToTextFile(const QString& fileName, c
for (RigFractureData fracData : fracDataVector)
{
out << qSetFieldWidth(8);
if (fracData.transmissibility == cvf::UNDEFINED_DOUBLE || !(fracture->attachedFractureDefinition())) out << "--"; //Commenting out line in output file
wellPath, simWell = nullptr;
fracture->firstAncestorOrThisOfType(simWell);
if (simWell) out << simWell->name; // 1. Well name
fracture->firstAncestorOrThisOfType(wellPath);
if (wellPath) out << wellPath->name; // 1. Well name
out << qSetFieldWidth(5);
size_t i, j, k;
mainGrid->ijkFromCellIndex(fracData.reservoirCellIndex, &i, &j, &k);
out << i+1; // 2. I location grid block, adding 1 to go to eclipse 1-based grid definition
out << j+1; // 3. J location grid block, adding 1 to go to eclipse 1-based grid definition
out << k+1; // 4. K location of upper connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << k+1; // 5. K location of lower connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << "2*"; // Default value for
//6. Open / Shut flag of connection
// 7. Saturation table number for connection rel perm. Default value
out << qSetFieldWidth(8);
// 8. Transmissibility
if (fracData.transmissibility != cvf::UNDEFINED_DOUBLE) out << fracData.transmissibility;
else out << "UNDEF";
out << qSetFieldWidth(4);
out << "2*"; // Default value for
// 9. Well bore diameter. Set to default
// 10. Effective Kh (perm times width)
if (fracture->attachedFractureDefinition())
{
out << fracture->attachedFractureDefinition()->skinFactor; // 11. Skin factor
}
else //If no attached fracture definition these parameters are set to UNDEF
{
out << "UNDEF";
}
out << "/";
out << " " << fracture->name(); //Fracture name as comment
out << "\n"; // Terminating entry
printCOMPDATvalues(out, fracData, fracture, wellPath, simWell, mainGrid);
}
//TODO: If same cell is used for multiple fractures, the sum of contributions should be added to table.
progress++;
@@ -229,3 +124,165 @@ bool RifEclipseExportTools::writeFracturesToTextFile(const QString& fileName, c
return true;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RifEclipseExportTools::printCOMPDATvalues(QTextStream & out, RigFractureData &fracData, RimFracture* fracture, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid)
{
out << qSetFieldWidth(8);
if (fracData.transmissibility == cvf::UNDEFINED_DOUBLE || !(fracture->attachedFractureDefinition())) out << "--"; //Commenting out line in output file
wellPath, simWell = nullptr;
fracture->firstAncestorOrThisOfType(simWell);
if (simWell) out << simWell->name; // 1. Well name
fracture->firstAncestorOrThisOfType(wellPath);
if (wellPath) out << wellPath->name; // 1. Well name
out << qSetFieldWidth(5);
size_t i, j, k;
mainGrid->ijkFromCellIndex(fracData.reservoirCellIndex, &i, &j, &k);
out << i + 1; // 2. I location grid block, adding 1 to go to eclipse 1-based grid definition
out << j + 1; // 3. J location grid block, adding 1 to go to eclipse 1-based grid definition
out << k + 1; // 4. K location of upper connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << k + 1; // 5. K location of lower connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << "2* "; // Default value for
//6. Open / Shut flag of connection
// 7. Saturation table number for connection rel perm. Default value
out << qSetFieldWidth(8);
// 8. Transmissibility
if (fracData.transmissibility != cvf::UNDEFINED_DOUBLE) out << QString::number(fracData.transmissibility, 'g', 6);
else out << "UNDEF";
out << qSetFieldWidth(4);
out << "2* "; // Default value for
// 9. Well bore diameter. Set to default
// 10. Effective Kh (perm times width)
if (fracture->attachedFractureDefinition())
{
out << fracture->attachedFractureDefinition()->skinFactor; // 11. Skin factor
}
else //If no attached fracture definition these parameters are set to UNDEF
{
out << "UNDEF";
}
out << "/";
out << " " << fracture->name(); //Fracture name as comment
out << "\n"; // Terminating entry
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RifEclipseExportTools::printBackgroundDataHeaderLine(QTextStream & out)
{
out << "-- Background data for calculation" << "\n\n";
//Write header line
out << qSetFieldWidth(4);
out << "--";
out << qSetFieldWidth(12);
out << "Well ";
out << qSetFieldWidth(16);
out << "Fracture ";
out << qSetFieldWidth(5);
out << "i";
out << "j";
out << "k";
out << qSetFieldWidth(12);
out << "Ax";
out << "Ay";
out << "Az";
out << "TotArea";
out << "skinfac";
out << "FracLen";
out << qSetFieldWidth(10);
out << "DX";
out << "DY";
out << "DZ";
out << "PermX";
out << "PermY";
out << "PermZ";
out << "NTG";
out << qSetFieldWidth(12);
out << "T_x";
out << "T_y";
out << "T_z";
out << qSetFieldWidth(15);
out << "Transm";
out << "\n";
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RifEclipseExportTools::printBackgroundData(QTextStream & out, RimWellPath* wellPath, RimEclipseWell* simWell, RimFracture* fracture, const RigMainGrid* mainGrid, RigFractureData &fracData)
{
out << qSetFieldWidth(4);
out << "--";
out << qSetFieldWidth(12);
wellPath, simWell = nullptr;
fracture->firstAncestorOrThisOfType(simWell);
if (simWell) out << simWell->name + " " ; // 1. Well name
fracture->firstAncestorOrThisOfType(wellPath);
if (wellPath) out << wellPath->name + " "; // 1. Well name
out << qSetFieldWidth(16);
out << fracture->name().left(15) + " ";
out << qSetFieldWidth(5);
size_t i, j, k;
mainGrid->ijkFromCellIndex(fracData.reservoirCellIndex, &i, &j, &k);
out << i + 1; // 2. I location grid block, adding 1 to go to eclipse 1-based grid definition
out << j + 1; // 3. J location grid block, adding 1 to go to eclipse 1-based grid definition
out << k + 1; // 4. K location of upper connecting grid block, adding 1 to go to eclipse 1-based grid definition
out << qSetFieldWidth(12);
//Use f for float, e for exponent float and g for best choice of these two.
out << QString::number(fracData.projectedAreas.x(), 'g', 4);
out << QString::number(fracData.projectedAreas.y(), 'g', 4);
out << QString::number(fracData.projectedAreas.z(), 'g', 4);
out << QString::number(fracData.totalArea, 'g', 4);
out << QString::number(fracData.skinFactor, 'f', 2);
out << QString::number(fracData.fractureLenght, 'g', 3);
out << qSetFieldWidth(10);
out << QString::number(fracData.cellSizes.x(), 'f', 2);
out << QString::number(fracData.cellSizes.y(), 'f', 2);
out << QString::number(fracData.cellSizes.z(), 'f', 2);
out << QString::number(fracData.permeabilities.x(), 'f', 4);
out << QString::number(fracData.permeabilities.y(), 'f', 4);
out << QString::number(fracData.permeabilities.z(), 'f', 4);
out << QString::number(fracData.NTG, 'f', 2);
out << qSetFieldWidth(12);
out << QString::number(fracData.transmissibilities.x(), 'e', 3);
out << QString::number(fracData.transmissibilities.y(), 'e', 3);
out << QString::number(fracData.transmissibilities.z(), 'e', 3);
out << qSetFieldWidth(15);
out << QString::number(fracData.transmissibility, 'f', 3);
out << "\n";
}

13
ApplicationCode/FileInterface/RifEclipseExportTools.h
View File

@@ -30,7 +30,12 @@
class RimFracture;
class QFile;
class QTextStream;
class RimWellPath;
class RimEclipseWell;
class RimFracture;
class RigMainGrid;
class RigFractureData;
@@ -47,5 +52,11 @@ public:
static bool writeFracturesToTextFile(const QString& fileName, const std::vector<RimFracture*>& fractures);
static void printCOMPDATvalues(QTextStream & out, RigFractureData &fracData, RimFracture* fracture, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid);
static void printBackgroundDataHeaderLine(QTextStream & out);
static void printBackgroundData(QTextStream & out, RimWellPath* wellPath, RimEclipseWell* simWell, RimFracture* fracture, const RigMainGrid* mainGrid, RigFractureData &fracData);
};

19
ApplicationCode/ProjectDataModel/RimFracture.cpp
View File

@@ -274,7 +274,6 @@ void RimFracture::computeTransmissibility()
size_t i, j, k;
mainGrid->ijkFromCellIndex(fracCell, &i, &j, &k);
caf::PdmObjectHandle* objHandle = dynamic_cast<caf::PdmObjectHandle*>(this);
if (!objHandle) return;
RimEclipseCase* eclipseCase;
@@ -284,8 +283,8 @@ void RimFracture::computeTransmissibility()
RifReaderInterface::PorosityModelResultType porosityModel = RigCaseCellResultsData::convertFromProjectModelPorosityModel(resultColors->porosityModel());
cvf::ref<RigResultAccessor> dataAccessObjectPermX = RigResultAccessorFactory::createFromUiResultName(eclipseCaseData, 0, porosityModel, 0, "PERMX"); //assuming 0 time step and main grid (so grid index =0)
cvf::ref<RigResultAccessor> dataAccessObjectPermY = RigResultAccessorFactory::createFromUiResultName(eclipseCaseData, 0, porosityModel, 0, "PERMX"); //assuming 0 time step and main grid (so grid index =0)
cvf::ref<RigResultAccessor> dataAccessObjectPermZ = RigResultAccessorFactory::createFromUiResultName(eclipseCaseData, 0, porosityModel, 0, "PERMX"); //assuming 0 time step and main grid (so grid index =0)
cvf::ref<RigResultAccessor> dataAccessObjectPermY = RigResultAccessorFactory::createFromUiResultName(eclipseCaseData, 0, porosityModel, 0, "PERMY"); //assuming 0 time step and main grid (so grid index =0)
cvf::ref<RigResultAccessor> dataAccessObjectPermZ = RigResultAccessorFactory::createFromUiResultName(eclipseCaseData, 0, porosityModel, 0, "PERMZ"); //assuming 0 time step and main grid (so grid index =0)
double permX = dataAccessObjectPermX->cellScalarGlobIdx(fracCell);
double permY = dataAccessObjectPermY->cellScalarGlobIdx(fracCell);
double permZ = dataAccessObjectPermZ->cellScalarGlobIdx(fracCell);
@@ -334,7 +333,9 @@ void RimFracture::computeTransmissibility()
double Ay = 0.0;
double Az = 0.0;
double skinfactor = 0.0;
double transmissibility_X = 0.0;
double transmissibility_Y = 0.0;
double transmissibility_Z = 0.0;
if (attachedFractureDefinition())
{
@@ -386,7 +387,7 @@ void RimFracture::computeTransmissibility()
double totalAreaXLength = 0.0;
for (double lengtXarea : lengthXareaOfFractureParts) totalAreaXLength += lengtXarea;
double fractureAreaWeightedlength = totalAreaXLength / fractureArea;
fractureAreaWeightedlength = totalAreaXLength / fractureArea;
double c = 0.008527; // TODO: Get value with units, is defined in RimReservoirCellResultsStorage
@@ -394,10 +395,9 @@ void RimFracture::computeTransmissibility()
double slDivPi = (skinfactor * fractureAreaWeightedlength) / cvf::PI_D;
//TODO: Use permx, permy, permz
double transmissibility_X = 8 * c * ( permY * NTG ) * Ay / (dx + slDivPi);
double transmissibility_Y = 8 * c * ( permX * NTG ) * Ax / (dy + slDivPi);
double transmissibility_Z = 8 * c * permZ * Az / (dz + slDivPi);
transmissibility_X = 8 * c * ( permY * NTG ) * Ay / (dx + slDivPi);
transmissibility_Y = 8 * c * ( permX * NTG ) * Ax / (dy + slDivPi);
transmissibility_Z = 8 * c * permZ * Az / (dz + slDivPi);
transmissibility = sqrt(transmissibility_X * transmissibility_X
+ transmissibility_Y * transmissibility_Y
@@ -410,6 +410,7 @@ void RimFracture::computeTransmissibility()
fracData.transmissibility = transmissibility;
fracData.transmissibilities = cvf::Vec3d(transmissibility_X, transmissibility_Y, transmissibility_Z);
fracData.totalArea = fractureArea;
fracData.projectedAreas = cvf::Vec3d(Ax, Ay, Az);

1
ApplicationCode/ReservoirDataModel/RigFracture.h
View File

@@ -33,6 +33,7 @@ public:
size_t reservoirCellIndex;
double transmissibility;
cvf::Vec3d transmissibilities;
double totalArea;
double fractureLenght;