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

#1560 Use radius, not diameter, when computing fishbones transmissibilities

Browse Source
This commit is contained in:
Bjørnar Grip Fjær
2017-06-08 11:56:06 +02:00
parent 0f188a2200
commit af95dacf03
1 changed files with 9 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
ApplicationCode/Commands/CompletionCommands/RicWellPathExportCompletionDataFeature.cpp
View File

@@ -189,6 +189,8 @@ void RicWellPathExportCompletionDataFeature::exportCompletions(const std::vector
{
generateWpimultTable(formatter, completions);
}
RiaLogging::info(QString("Successfully exported completion data to %1").arg(exportSettings.fileName()));
}
//--------------------------------------------------------------------------------------------------
@@ -337,7 +339,7 @@ std::vector<RigCompletionData> RicWellPathExportCompletionDataFeature::generateF
wellPath,
intersection.lengthsInCell,
location.fishbonesSubs->skinFactor(),
location.fishbonesSubs->holeDiameter(),
diameter / 2,
intersection.cellIndex);
completion.setFromFishbone(transmissibility, location.fishbonesSubs->skinFactor());
}
@@ -360,12 +362,16 @@ std::vector<RigCompletionData> RicWellPathExportCompletionDataFeature::generateF
{
std::vector<RigCompletionData> completionData;
std::vector<size_t> wellPathCells = findIntersectingCells(settings.caseToApply()->eclipseCaseData(), wellPath->wellPathGeometry()->m_wellPathPoints);
double diameter = wellPath->fishbonesCollection()->wellPathCollection()->holeDiameter() / 1000;
for (const RimFishboneWellPath* fishbonesPath : wellPath->fishbonesCollection()->wellPathCollection()->wellPaths())
{
std::vector<WellPathCellIntersectionInfo> intersectedCells = RigWellPathIntersectionTools::findCellsIntersectedByPath(settings.caseToApply->eclipseCaseData(), fishbonesPath->coordinates());
for (auto& cell : intersectedCells)
{
if (std::find(wellPathCells.begin(), wellPathCells.end(), cell.cellIndex) != wellPathCells.end()) continue;
size_t i, j, k;
settings.caseToApply->eclipseCaseData()->mainGrid()->ijkFromCellIndex(cell.cellIndex, &i, &j, &k);
RigCompletionData completion(wellPath->name(), IJKCellIndex(i, j, k));
@@ -373,11 +379,11 @@ std::vector<RigCompletionData> RicWellPathExportCompletionDataFeature::generateF
if (settings.computeTransmissibility())
{
double skinFactor = wellPath->fishbonesCollection()->wellPathCollection()->skinFactor();
double transmissibility = calculateTransmissibility(settings.caseToApply,
double transmissibility = calculateTransmissibility(settings.caseToApply(),
wellPath,
cell.internalCellLengths,
skinFactor,
wellPath->fishbonesCollection()->wellPathCollection()->holeDiameter(),
diameter / 2,
cell.cellIndex);
completion.setFromFishbone(transmissibility, skinFactor);
}