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

(#675) (#666) Extrapolating intersection by user amount, along well-end directions

Browse Source 
Making vertical-only well paths more robust also.
This commit is contained in:
Jacob Støren 2015-11-26 15:59:02 +01:00
parent 3c4142fc52
commit 7ebd6514cb
2 changed files with 96 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

107
ApplicationCode/ProjectDataModel/RimCrossSection.cpp
View File

@ -67,12 +67,13 @@ RimCrossSection::RimCrossSection()
CAF_PDM_InitField(&isActive, "Active", true, "Active", "", "", "");
isActive.uiCapability()->setUiHidden(true);
CAF_PDM_InitFieldNoDefault(&type, "Type", "Type", "", "", "");
CAF_PDM_InitFieldNoDefault(&direction, "Direction", "Direction", "", "", "");
CAF_PDM_InitFieldNoDefault(&wellPath, "WellPath", "Well Path", "", "", "");
CAF_PDM_InitFieldNoDefault(&simulationWell, "SimulationWell", "Simulation Well", "", "", "");
CAF_PDM_InitField (&branchIndex, "Branch", -1, "Branch", "", "", "");
CAF_PDM_InitFieldNoDefault(&type, "Type", "Type", "", "", "");
CAF_PDM_InitFieldNoDefault(&direction, "Direction", "Direction", "", "", "");
CAF_PDM_InitFieldNoDefault(&wellPath, "WellPath", "Well Path", "", "", "");
CAF_PDM_InitFieldNoDefault(&simulationWell, "SimulationWell", "Simulation Well", "", "", "");
CAF_PDM_InitField (&m_branchIndex, "Branch", -1, "Branch", "", "", "");
CAF_PDM_InitField (&m_extentLength, "ExtentLength", 200.0, "Extent length", "", "", "");
uiCapability()->setUiChildrenHidden(true);
}
@ -86,7 +87,8 @@ void RimCrossSection::fieldChangedByUi(const caf::PdmFieldHandle* changedField,
changedField == &direction ||
changedField == &wellPath ||
changedField == &simulationWell ||
changedField == &branchIndex)
changedField == &m_branchIndex ||
changedField == &m_extentLength)
{
m_crossSectionPartMgr = NULL;
@ -104,7 +106,15 @@ void RimCrossSection::fieldChangedByUi(const caf::PdmFieldHandle* changedField,
{
m_wellBranchCenterlines.clear();
updateWellCenterline();
branchIndex = -1;
if (m_wellBranchCenterlines.size())
{
cvf::Vec3d wellLength = m_wellBranchCenterlines[0].front() - m_wellBranchCenterlines[0].back();
if (wellLength.length() < 200)
{
}
}
m_branchIndex = -1;
}
}
@ -120,6 +130,7 @@ void RimCrossSection::defineUiOrdering(QString uiConfigName, caf::PdmUiOrdering&
if (type == CS_WELL_PATH)
{
uiOrdering.add(&wellPath);
uiOrdering.add(&m_extentLength);
}
else if (type == CS_SIMULATION_WELL)
{
@ -127,14 +138,16 @@ void RimCrossSection::defineUiOrdering(QString uiConfigName, caf::PdmUiOrdering&
updateWellCenterline();
if (simulationWell() && m_wellBranchCenterlines.size() > 1)
{
uiOrdering.add(&branchIndex);
uiOrdering.add(&m_branchIndex);
}
uiOrdering.add(&m_extentLength);
}
else
{
// User defined poly line
}
uiOrdering.setForgetRemainingFields(true);
}
@ -184,7 +197,7 @@ QList<caf::PdmOptionItemInfo> RimCrossSection::calculateValueOptions(const caf::
options.push_front(caf::PdmOptionItemInfo("None", QVariant::fromValue(caf::PdmPointer<caf::PdmObjectHandle>(NULL))));
}
}
else if (fieldNeedingOptions == &branchIndex)
else if (fieldNeedingOptions == &m_branchIndex)
{
updateWellCenterline();
@ -251,12 +264,12 @@ std::vector< std::vector <cvf::Vec3d> > RimCrossSection::polyLines() const
{
updateWellCenterline();
if (0 <= branchIndex && branchIndex < m_wellBranchCenterlines.size())
if (0 <= m_branchIndex && m_branchIndex < m_wellBranchCenterlines.size())
{
lines.push_back(m_wellBranchCenterlines[branchIndex]);
lines.push_back(m_wellBranchCenterlines[m_branchIndex]);
}
if (branchIndex == -1)
if (m_branchIndex == -1)
{
lines = m_wellBranchCenterlines;
}
@ -271,14 +284,10 @@ std::vector< std::vector <cvf::Vec3d> > RimCrossSection::polyLines() const
{
for (int lIdx = 0; lIdx < lines.size(); ++lIdx)
{
cvf::Vec3d startToEnd = (lines[lIdx].back() - lines[lIdx].front());
startToEnd[2] = 0.0;
cvf::Vec3d newStart = lines[lIdx].front() - startToEnd * 0.1;
cvf::Vec3d newEnd = lines[lIdx].back() + startToEnd * 0.1;
std::vector<cvf::Vec3d>& polyLine = lines[lIdx];
addExtents(polyLine);
lines[lIdx].insert(lines[lIdx].begin(), newStart);
lines[lIdx].push_back(newEnd);
}
}
@ -315,3 +324,61 @@ void RimCrossSection::updateWellCenterline() const
}
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
void RimCrossSection::addExtents(std::vector<cvf::Vec3d> &polyLine) const
{
size_t lineVxCount = polyLine.size();
if (lineVxCount == 0) return;
// Add extent at end of well
{
size_t endIdxOffset = lineVxCount > 3 ? 3: lineVxCount;
cvf::Vec3d endDirection = (polyLine[lineVxCount-1] - polyLine[lineVxCount-endIdxOffset]);
endDirection[2] = 0; // Remove z. make extent lenght be horizontally
if (endDirection.length() < 1e-2)
{
endDirection = polyLine.back() - polyLine.front();
endDirection[2] = 0;
if (endDirection.length() < 1e-2)
{
endDirection = cvf::Vec3d::X_AXIS;
}
}
endDirection.normalize();
cvf::Vec3d newEnd = polyLine.back() + endDirection * m_extentLength();
polyLine.push_back(newEnd);
}
// Add extent at start
{
size_t endIdxOffset = lineVxCount > 3 ? 3: lineVxCount-1;
cvf::Vec3d startDirection = (polyLine[0] - polyLine[endIdxOffset]);
startDirection[2] = 0; // Remove z. make extent lenght be horizontally
if (startDirection.length() < 1e-2)
{
startDirection = polyLine.front() - polyLine.back();
startDirection[2] = 0;
if (startDirection.length() < 1e-2)
{
startDirection = -cvf::Vec3d::X_AXIS;
}
}
startDirection.normalize();
cvf::Vec3d newStart = polyLine.front() + startDirection * m_extentLength();
polyLine.insert(polyLine.begin(), newStart);
}
}

12
ApplicationCode/ProjectDataModel/RimCrossSection.h
View File

@ -67,22 +67,28 @@ public:
caf::PdmPtrField<RimWellPath*> wellPath;
caf::PdmPtrField<RimEclipseWell*> simulationWell;
caf::PdmField<int> branchIndex;
std::vector< std::vector <cvf::Vec3d> > polyLines() const;
virtual caf::PdmFieldHandle* userDescriptionField();
virtual caf::PdmFieldHandle* objectToggleField();
void addExtents(std::vector<cvf::Vec3d> &polyLine) const;
RivCrossSectionPartMgr* crossSectionPartMgr();
protected:
virtual caf::PdmFieldHandle* userDescriptionField();
virtual caf::PdmFieldHandle* objectToggleField();
virtual void fieldChangedByUi(const caf::PdmFieldHandle* changedField, const QVariant& oldValue, const QVariant& newValue);
virtual void defineUiOrdering(QString uiConfigName, caf::PdmUiOrdering& uiOrdering);
virtual QList<caf::PdmOptionItemInfo> calculateValueOptions(const caf::PdmFieldHandle* fieldNeedingOptions, bool * useOptionsOnly);
private:
caf::PdmField<int> m_branchIndex;
caf::PdmField<double> m_extentLength;
RimEclipseWellCollection* simulationWellCollection();
void updateWellCenterline() const;