OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-10 23:46:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

#1403 - pre-proto - Adding stimPlan cells transmissibility in fracture to export.

Browse Source
This commit is contained in:
astridkbjorke 2017-04-18 13:34:32 +02:00
parent 2595bf08fe
commit 94d10f547d
6 changed files with 100 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

104
ApplicationCode/FileInterface/RifEclipseExportTools.cpp
View File

@ -103,8 +103,8 @@ bool RifEclipseExportTools::writeFracturesToTextFile(const QString& fileName, c
out << "\n";
//Included for debug / prototyping only
printStimPlanFractureTrans(fractures, out);
printStimPlanCellsMatrixTransContributions(fractures, caseToApply, out, wellPath, simWell, mainGrid);
printBackgroundDataHeaderLine(out);
@ -235,7 +235,36 @@ void RifEclipseExportTools::performStimPlanUpscalingAndPrintResults(const std::v
//--------------------------------------------------------------------------------------------------
void RifEclipseExportTools::printStimPlanCellsMatrixTransContributions(const std::vector<RimFracture *>& fractures, RimEclipseCase* caseToApply, QTextStream &out, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid)
{
out << "StimPlan cells' matrix transmissibility and contributing \n";
out << "StimPlan cells' matrix transmissibility and Eclipse Cell contributions \n";
out << qSetFieldWidth(4);
out << "-- ";
out << qSetFieldWidth(12);
out << "Well name "; // 1. Well name
out << qSetFieldWidth(16);
out << "Fracture name ";
out << qSetFieldWidth(5);
out << "Ec i";
out << "Ec j";
out << "Ec k";
out << qSetFieldWidth(10);
out << "Ecl cell";
out << qSetFieldWidth(5);
out << "SP i";
out << "SP j";
out << qSetFieldWidth(10);
out << "Tm contr";
out << "\n";
for (RimFracture* fracture : fractures)
{
RimStimPlanFractureTemplate* fracTemplateStimPlan;
@ -247,14 +276,14 @@ void RifEclipseExportTools::printStimPlanCellsMatrixTransContributions(const std
RigFractureTransCalc transmissibilityCalculator(caseToApply, fracture);
//TODO: Get these more generally:
std::vector<RigStimPlanCell> stimPlanCells = fracTemplateStimPlan->getStimPlanCells();
for (RigStimPlanCell stimPlanCell : stimPlanCells)
{
if (stimPlanCell.getConductivtyValue() < 1e-7)
{
//Adding empty fracStimPlanCell to get vectors of same length
//If conductivity in stimPlanCell is 0, contributions might not be relevant...
//continue;
}
RigFractureStimPlanCellData fracStimPlanCellData(stimPlanCell.getI(), stimPlanCell.getJ());
@ -281,7 +310,6 @@ void RifEclipseExportTools::printStimPlanCellsMatrixTransContributions(const std
out << qSetFieldWidth(16);
out << fracture->name().left(15) + " ";
out << qSetFieldWidth(5);
size_t ii, jj, kk;
mainGrid->ijkFromCellIndex(stimPlanContributingEclipseCells[i], &ii, &jj, &kk);
@ -303,18 +331,72 @@ void RifEclipseExportTools::printStimPlanCellsMatrixTransContributions(const std
out << QString::number(stimPlanContributingEclipseCellTransmisibilities[i], 'e', 3);
out << "\n";
}
}
}
return;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RifEclipseExportTools::printStimPlanFractureTrans(const std::vector<RimFracture *>& fractures, QTextStream &out)
{
out << "StimPlan cells' fracture transmissibility \n";
out << qSetFieldWidth(4);
out << "-- ";
out << qSetFieldWidth(5);
out << "SP i";
out << "SP j";
out << qSetFieldWidth(10);
out << "Tf_hor";
out << "Tf_vert";
out << "\n";
if (fractures.size() < 1) return;
RimFracture* fracture = fractures[0];
RimStimPlanFractureTemplate* fracTemplateStimPlan;
if (dynamic_cast<RimStimPlanFractureTemplate*>(fracture->attachedFractureDefinition()))
{
fracTemplateStimPlan = dynamic_cast<RimStimPlanFractureTemplate*>(fracture->attachedFractureDefinition());
}
else return;
std::vector<RigStimPlanCell> stimPlanCells = fracTemplateStimPlan->getStimPlanCells();
for (RigStimPlanCell stimPlanCell : stimPlanCells)
{
if (stimPlanCell.getConductivtyValue() < 1e-7)
{
//If conductivity in stimPlanCell is 0, contributions might not be relevant...
continue;
}
RigFractureTransCalc::computeStimPlanCellTransmissibilityInFracture(&stimPlanCell);
out << qSetFieldWidth(5);
size_t spi = stimPlanCell.getI();
size_t spj = stimPlanCell.getJ();
out << spi;
out << spj;
out << qSetFieldWidth(10);
out << QString::number(stimPlanCell.getVerticalTransmissibilityInFracture(), 'e', 3);
out << QString::number(stimPlanCell.getHorizontalTransmissibilityInFracture(), 'e', 3);
out << "\n";
}
return;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------

2
ApplicationCode/FileInterface/RifEclipseExportTools.h
View File

@ -53,6 +53,8 @@ public:
static void performStimPlanUpscalingAndPrintResults(const std::vector<RimFracture *>& fractures, RimEclipseCase* caseToApply, QTextStream &out, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid);
static void printStimPlanCellsMatrixTransContributions(const std::vector<RimFracture *>& fractures, RimEclipseCase* caseToApply, QTextStream &out, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid);
static void printStimPlanFractureTrans(const std::vector<RimFracture *>& fractures, QTextStream &out);
static void printCOMPDATvalues(QTextStream & out, RigFractureData &fracData, RimFracture* fracture, RimWellPath* wellPath, RimEclipseWell* simWell, const RigMainGrid* mainGrid);

2
ApplicationCode/ReservoirDataModel/RigFractureTransCalc.cpp
View File

@ -727,7 +727,7 @@ void RigFractureTransCalc::computeUpscaledPropertyFromStimPlan( QString resultNa
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RigFractureTransCalc::computeFlowInFracture(RigStimPlanCell* stimPlanCell)
void RigFractureTransCalc::computeStimPlanCellTransmissibilityInFracture(RigStimPlanCell* stimPlanCell)
{
std::vector<cvf::Vec3d> polygon = stimPlanCell->getPolygon();

2
ApplicationCode/ReservoirDataModel/RigFractureTransCalc.h
View File

@ -59,7 +59,7 @@ public:
void calculateStimPlanCellsMatrixTransmissibility(RigStimPlanCell* stimPlanCell, RigFractureStimPlanCellData* fracStimPlanCellData);
void computeFlowInFracture(RigStimPlanCell* stimPlanCell);
static void computeStimPlanCellTransmissibilityInFracture(RigStimPlanCell* stimPlanCell);
void computeFlowIntoTransverseWell();

9
ApplicationCode/ReservoirDataModel/RigStimPlanCell.cpp
View File

@ -53,13 +53,4 @@ void RigStimPlanCell::setTransmissibilityInFracture(double valueHorizontal, doub
m_transmissibilityInFractureVertical = valueVertical;
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
// void RigStimPlanCell::addContributingEclipseCell(size_t eclipseCell, double transmissibility)
// {
// contributingEclipseCells.push_back(eclipseCell);
// contributingEclipseCellTransmisibilities.push_back(transmissibility);
//
// }

5
ApplicationCode/ReservoirDataModel/RigStimPlanCell.h
View File

@ -42,13 +42,14 @@ public:
double getDisplayValue() { return m_displayValue; }
size_t getI() { return m_i; }
size_t getJ() { return m_j; }
double getVerticalTransmissibilityInFracture() { return m_transmissibilityInFractureVertical; }
double getHorizontalTransmissibilityInFracture() { return m_transmissibilityInFractureHorizontal; }
void setConductivityValue(double cond) { m_concutivityValue = cond; }
void setDisplayValue(double value) { m_displayValue = value; };
void setTransmissibilityInFracture(double valueHorizontal, double valueVertical);
// void addContributingEclipseCell(size_t fracCell, double transmissibility);
private:
std::vector<cvf::Vec3d> m_polygon;
double m_displayValue;