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

#3855 Limit entry and exit lengths in perforation COMPSEGS when perforation doesn't pass through entire cell.

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-12-12 09:44:00 +01:00
parent 762b6ef382
commit 509f77c7f1
4 changed files with 52 additions and 30 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
ApplicationCode/Commands/CompletionExportCommands/RicMswSubSegment.cpp
View File

@ -67,11 +67,11 @@ const cvf::Vec3d& RicMswSubSegmentCellIntersection::lengthsInCell() const
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
RicMswSubSegment::RicMswSubSegment(double startMD, double deltaMD, double startTVD, double deltaTVD) RicMswSubSegment::RicMswSubSegment(double startMD, double endMD, double startTVD, double endTVD)
: m_startMD(startMD) : m_startMD(startMD)
, m_deltaMD(deltaMD) , m_endMD(endMD)
, m_startTVD(startTVD) , m_startTVD(startTVD)
, m_deltaTVD(deltaTVD) , m_endTVD(endTVD)
, m_segmentNumber(-1) , m_segmentNumber(-1)
, m_attachedSegmentNumber(-1) , m_attachedSegmentNumber(-1)
{ {
@ -85,12 +85,20 @@ double RicMswSubSegment::startMD() const
return m_startMD; return m_startMD;
} }
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RicMswSubSegment::endMD() const
{
return m_endMD;
}
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
double RicMswSubSegment::deltaMD() const double RicMswSubSegment::deltaMD() const
{ {
return m_deltaMD; return m_endMD - m_startMD;
} }
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
@ -101,12 +109,20 @@ double RicMswSubSegment::startTVD() const
return m_startTVD; return m_startTVD;
} }
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RicMswSubSegment::endTVD() const
{
return m_endTVD;
}
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
/// ///
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------
double RicMswSubSegment::deltaTVD() const double RicMswSubSegment::deltaTVD() const
{ {
return m_deltaTVD; return m_endTVD - m_startTVD;
} }
//-------------------------------------------------------------------------------------------------- //--------------------------------------------------------------------------------------------------

10
ApplicationCode/Commands/CompletionExportCommands/RicMswSubSegment.h
View File

@ -52,13 +52,15 @@ class RicMswSubSegment
{ {
public: public:
RicMswSubSegment(double startMD, RicMswSubSegment(double startMD,
double deltaMD, double endMD,
double startTVD, double startTVD,
double deltaTVD); double endTVD);
double startMD() const; double startMD() const;
double endMD() const;
double deltaMD() const; double deltaMD() const;
double startTVD() const; double startTVD() const;
double endTVD() const;
double deltaTVD() const; double deltaTVD() const;
int segmentNumber() const; int segmentNumber() const;
@ -74,9 +76,9 @@ public:
private: private:
double m_startMD; double m_startMD;
double m_deltaMD; double m_endMD;
double m_startTVD; double m_startTVD;
double m_deltaTVD; double m_endTVD;
int m_segmentNumber; int m_segmentNumber;
int m_attachedSegmentNumber; int m_attachedSegmentNumber;

38
ApplicationCode/Commands/CompletionExportCommands/RicWellPathExportMswCompletionsImpl.cpp
View File

@ -403,8 +403,8 @@ void RicWellPathExportMswCompletionsImpl::generateWelsegsSegments(
} }
else else
{ {
depth = subSegment->startTVD() + subSegment->deltaTVD(); depth = subSegment->endTVD();
length = subSegment->startMD() + subSegment->deltaMD(); length = subSegment->endMD();
} }
double diameter = segment->effectiveDiameter(); double diameter = segment->effectiveDiameter();
formatter.add(subSegment->segmentNumber()); formatter.add(subSegment->segmentNumber());
@ -520,14 +520,14 @@ void RicWellPathExportMswCompletionsImpl::generateCompsegTable(
generatedHeader = true; generatedHeader = true;
} }
for (std::shared_ptr<RicMswSubSegment> segment : completion->subSegments()) for (const std::shared_ptr<RicMswSubSegment>& subSegment : completion->subSegments())
{ {
if (completion->completionType() == RigCompletionData::FISHBONES_ICD) if (completion->completionType() == RigCompletionData::FISHBONES_ICD)
{ {
startMD = segment->startMD(); startMD = subSegment->startMD();
} }
for (std::shared_ptr<RicMswSubSegmentCellIntersection> intersection : segment->intersections()) for (const std::shared_ptr<RicMswSubSegmentCellIntersection>& intersection : subSegment->intersections())
{ {
bool isSubGridIntersection = !intersection->gridName().isEmpty(); bool isSubGridIntersection = !intersection->gridName().isEmpty();
if (isSubGridIntersection == exportSubGridIntersections) if (isSubGridIntersection == exportSubGridIntersections)
@ -540,14 +540,14 @@ void RicWellPathExportMswCompletionsImpl::generateCompsegTable(
formatter.addOneBasedCellIndex(ijk.x()).addOneBasedCellIndex(ijk.y()).addOneBasedCellIndex(ijk.z()); formatter.addOneBasedCellIndex(ijk.x()).addOneBasedCellIndex(ijk.y()).addOneBasedCellIndex(ijk.z());
formatter.add(completion->branchNumber()); formatter.add(completion->branchNumber());
double startLength = segment->startMD(); double startLength = subSegment->startMD();
if (exportInfo.lengthAndDepthText() == QString("INC") && if (exportInfo.lengthAndDepthText() == QString("INC") &&
completion->completionType() != RigCompletionData::PERFORATION) completion->branchNumber() != 1)
{ {
startLength -= startMD; startLength -= startMD;
} }
formatter.add(startLength); formatter.add(startLength);
formatter.add(startLength + segment->deltaMD()); formatter.add(startLength + subSegment->deltaMD());
formatter.rowCompleted(); formatter.rowCompleted();
} }
@ -848,7 +848,7 @@ RicMswExportInfo RicWellPathExportMswCompletionsImpl::generateFishbonesMswExport
{ {
// Add completion for ICD // Add completion for ICD
std::shared_ptr<RicMswFishbonesICD> icdCompletion(new RicMswFishbonesICD(QString("ICD"))); std::shared_ptr<RicMswFishbonesICD> icdCompletion(new RicMswFishbonesICD(QString("ICD")));
std::shared_ptr<RicMswSubSegment> icdSegment(new RicMswSubSegment(subEndMD, 0.1, subEndTVD, 0.0)); std::shared_ptr<RicMswSubSegment> icdSegment(new RicMswSubSegment(subEndMD, subEndMD + 0.1, subEndTVD, subEndTVD));
icdCompletion->setFlowCoefficient(subs->icdFlowCoefficient()); icdCompletion->setFlowCoefficient(subs->icdFlowCoefficient());
double icdOrificeRadius = subs->icdOrificeDiameter(unitSystem) / 2; double icdOrificeRadius = subs->icdOrificeDiameter(unitSystem) / 2;
icdCompletion->setArea(icdOrificeRadius * icdOrificeRadius * cvf::PI_D * subs->icdCount()); icdCompletion->setArea(icdOrificeRadius * icdOrificeRadius * cvf::PI_D * subs->icdCount());
@ -1095,7 +1095,7 @@ RicWellPathExportMswCompletionsImpl::MainBoreSegments RicWellPathExportMswComple
std::shared_ptr<RicMswCompletion> intervalCompletion(new RicMswPerforation(interval->name())); std::shared_ptr<RicMswCompletion> intervalCompletion(new RicMswPerforation(interval->name()));
std::vector<RigCompletionData> completionData = std::vector<RigCompletionData> completionData =
generatePerforationIntersections(wellPath, interval, timeStep, eclipseCase); generatePerforationIntersections(wellPath, interval, timeStep, eclipseCase);
assignPerforationIntersections(completionData, intervalCompletion, cellIntInfo, foundSubGridIntersections); assignPerforationIntersections(completionData, intervalCompletion, cellIntInfo, overlapStart, overlapEnd, foundSubGridIntersections);
segment->addCompletion(intervalCompletion); segment->addCompletion(intervalCompletion);
} }
} }
@ -1146,7 +1146,7 @@ void RicWellPathExportMswCompletionsImpl::assignSuperValveCompletions(
if (segment->startMD() <= valveMD && valveMD < segment->endMD()) if (segment->startMD() <= valveMD && valveMD < segment->endMD())
{ {
QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2); QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2);
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(valveMD, 0.1, 0.0, 0.0)); std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(valveMD, valveMD + 0.1, 0.0, 0.0));
if (isAicd) if (isAicd)
{ {
@ -1162,14 +1162,14 @@ void RicWellPathExportMswCompletionsImpl::assignSuperValveCompletions(
else if (overlap > 0.0 && (!isAicd && !superICD)) else if (overlap > 0.0 && (!isAicd && !superICD))
{ {
QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2); QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2);
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(overlapStart, 0.1, 0.0, 0.0)); std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(overlapStart, overlapStart + 0.1, 0.0, 0.0));
superICD.reset(new RicMswPerforationICD(valveLabel)); superICD.reset(new RicMswPerforationICD(valveLabel));
superICD->addSubSegment(subSegment); superICD->addSubSegment(subSegment);
} }
else if (overlap > 0.0 && (isAicd && !superAICD)) else if (overlap > 0.0 && (isAicd && !superAICD))
{ {
QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2); QString valveLabel = QString("%1 #%2").arg("Combined Valve for segment").arg(nMainSegment + 2);
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(overlapStart, 0.1, 0.0, 0.0)); std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(overlapStart, overlapStart + 0.1, 0.0, 0.0));
superAICD.reset(new RicMswPerforationAICD(valveLabel)); superAICD.reset(new RicMswPerforationAICD(valveLabel));
superAICD->addSubSegment(subSegment); superAICD->addSubSegment(subSegment);
} }
@ -1393,7 +1393,7 @@ void RicWellPathExportMswCompletionsImpl::assignFishbonesLateralIntersections(co
size_t i = 0u, j = 0u, k = 0u; size_t i = 0u, j = 0u, k = 0u;
localGrid->ijkFromCellIndex(localGridIdx, &i, &j, &k); localGrid->ijkFromCellIndex(localGridIdx, &i, &j, &k);
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment( std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(
previousExitMD, cellIntInfo.endMD - previousExitMD, previousExitTVD, cellIntInfo.endTVD - previousExitTVD)); previousExitMD, cellIntInfo.endMD, previousExitTVD, cellIntInfo.endTVD));
std::shared_ptr<RicMswSubSegmentCellIntersection> intersection(new RicMswSubSegmentCellIntersection( std::shared_ptr<RicMswSubSegmentCellIntersection> intersection(new RicMswSubSegmentCellIntersection(
gridName, cellIntInfo.globCellIndex, cvf::Vec3st(i, j, k), cellIntInfo.intersectionLengthsInCellCS)); gridName, cellIntInfo.globCellIndex, cvf::Vec3st(i, j, k), cellIntInfo.intersectionLengthsInCellCS));
@ -1428,7 +1428,7 @@ void RicWellPathExportMswCompletionsImpl::assignFractureIntersections(const RimE
width = perforationLength; width = perforationLength;
} }
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(position, width, 0.0, 0.0)); std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(position, position + width, 0.0, 0.0));
for (const RigCompletionData& compIntersection : completionData) for (const RigCompletionData& compIntersection : completionData)
{ {
const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell(); const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell();
@ -1487,14 +1487,16 @@ std::vector<RigCompletionData>
void RicWellPathExportMswCompletionsImpl::assignPerforationIntersections(const std::vector<RigCompletionData>& completionData, void RicWellPathExportMswCompletionsImpl::assignPerforationIntersections(const std::vector<RigCompletionData>& completionData,
std::shared_ptr<RicMswCompletion> perforationCompletion, std::shared_ptr<RicMswCompletion> perforationCompletion,
const SubSegmentIntersectionInfo& cellIntInfo, const SubSegmentIntersectionInfo& cellIntInfo,
double overlapStart,
double overlapEnd,
bool* foundSubGridIntersections) bool* foundSubGridIntersections)
{ {
size_t currCellId = cellIntInfo.globCellIndex; size_t currCellId = cellIntInfo.globCellIndex;
std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(cellIntInfo.startMD, std::shared_ptr<RicMswSubSegment> subSegment(new RicMswSubSegment(overlapStart,
cellIntInfo.endMD - cellIntInfo.startMD, overlapEnd,
cellIntInfo.startTVD, cellIntInfo.startTVD,
cellIntInfo.endTVD - cellIntInfo.startTVD)); cellIntInfo.endTVD));
for (const RigCompletionData& compIntersection : completionData) for (const RigCompletionData& compIntersection : completionData)
{ {
const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell(); const RigCompletionDataGridCell& cell = compIntersection.completionDataGridCell();

8
ApplicationCode/Commands/CompletionExportCommands/RicWellPathExportMswCompletionsImpl.h
View File

@ -134,9 +134,11 @@ private:
RimEclipseCase* eclipseCase); RimEclipseCase* eclipseCase);
static void assignPerforationIntersections(const std::vector<RigCompletionData>& completionData, static void assignPerforationIntersections(const std::vector<RigCompletionData>& completionData,
std::shared_ptr<RicMswCompletion> perforationCompletion, std::shared_ptr<RicMswCompletion> perforationCompletion,
const SubSegmentIntersectionInfo& cellIntInfo, const SubSegmentIntersectionInfo& cellIntInfo,
bool* foundSubGridIntersections); double overlapStart,
double overlapEnd,
bool* foundSubGridIntersections);
static void assignBranchAndSegmentNumbers(const RimEclipseCase* caseToApply, static void assignBranchAndSegmentNumbers(const RimEclipseCase* caseToApply,
std::shared_ptr<RicMswSegment> location, std::shared_ptr<RicMswSegment> location,