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

Rft curve: Fix first depth value

Browse Source
This commit is contained in:
Rebecca Cox
2017-11-01 11:47:47 +01:00
parent e40d41fe68
commit 8d388bbd0e
1 changed files with 22 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

29
ApplicationCode/ProjectDataModel/RimWellLogRftCurve.cpp
View File

@@ -168,7 +168,12 @@ void RimWellLogRftCurve::setDefaultAddress(QString wellName)
}
}
if (!wellNameHasRftData) return;
if (!wellNameHasRftData)
{
m_wellLogChannelName = RifEclipseRftAddress::NONE;
m_timeStep = QDateTime();
return;
}
m_wellLogChannelName = RifEclipseRftAddress::PRESSURE;
@@ -474,12 +479,13 @@ bool RimWellLogRftCurve::createWellPathIdxToRftFileIdxMapping()
if (!eclExtractor) return false;
std::vector<CellIntersectionInfo> intersections = eclExtractor->intersectionInfo();
if (intersections.empty()) return false;
std::map<size_t, size_t> globCellIndicesToIndexInWell;
for (size_t i = 0; i < intersections.size(); i++)
for (size_t idx = 0; idx < intersections.size(); idx++)
{
globCellIndicesToIndexInWell[intersections[i].globCellIndex] = i;
globCellIndicesToIndexInWell[intersections[idx].globCellIndex] = idx;
}
RifEclipseRftAddress depthAddress(m_wellName(), m_timeStep, RifEclipseRftAddress::DEPTH);
@@ -488,12 +494,15 @@ bool RimWellLogRftCurve::createWellPathIdxToRftFileIdxMapping()
const RigMainGrid* mainGrid = eclExtractor->caseData()->mainGrid();
for (size_t i = 0; i < rftIndices.size(); i++)
for (size_t idx = 0; idx < rftIndices.size(); idx++)
{
caf::VecIjk ijkIndex = rftIndices[i];
caf::VecIjk ijkIndex = rftIndices[idx];
size_t globalCellIndex = mainGrid->cellIndexFromIJK(ijkIndex.i(), ijkIndex.j(), ijkIndex.k());
m_idxInWellPathToIdxInRftFile[globCellIndicesToIndexInWell[globalCellIndex]] = i;
if (globCellIndicesToIndexInWell.find(globalCellIndex) != globCellIndicesToIndexInWell.end())
{
m_idxInWellPathToIdxInRftFile[globCellIndicesToIndexInWell[globalCellIndex]] = idx;
}
}
return true;
@@ -619,9 +628,15 @@ std::vector<double> RimWellLogRftCurve::measuredDepthValues()
if (!eclExtractor) return measuredDepthForCells;
std::vector<double> measuredDepthForIntersections = eclExtractor->measuredDepth();
if (measuredDepthForIntersections.empty())
{
return measuredDepthForCells;
}
std::vector< size_t> globCellIndices = eclExtractor->intersectedCellsGlobIdx();
for (size_t i = 0; i < globCellIndices.size() - 1; i = i + 2)
{
double sum = measuredDepthForIntersections[i] + measuredDepthForIntersections[i + 1];