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

#2232 Well formations: Read new file format

Browse Source
This commit is contained in:
Rebecca Cox 2017-12-07 08:28:20 +01:00
parent 89600393ab
commit 9a4eab2569
1 changed files with 48 additions and 17 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

65
ApplicationCode/FileInterface/RifWellPathFormationReader.cpp
View File

@ -78,40 +78,71 @@ void RifWellPathFormationReader::readFileIntoMap(const QString&
while (header.size() < 3) while (header.size() < 3)
{ {
if (data.atEnd()) return;
QString line = data.readLine().toLower(); QString line = data.readLine().toLower();
removeWhiteSpaces(&line); removeWhiteSpaces(&line);
header = line.split(';', QString::KeepEmptyParts); header = line.split(';', QString::KeepEmptyParts);
} }
const QString wellNameText = "wellname"; static const QString wellNameText = "wellname";
const QString surfaceNameText = "surfacename"; static const QString surfaceNameText = "surfacename";
const QString measuredDepthText = "md"; static const QString measuredDepthText = "md";
int wellNameIndex = header.indexOf(wellNameText); int wellNameIndex = header.indexOf(wellNameText);
int surfaceNameIndex = header.indexOf(surfaceNameText); int surfaceNameIndex = header.indexOf(surfaceNameText);
int measuredDepthIndex = header.indexOf(measuredDepthText); int measuredDepthIndex = header.indexOf(measuredDepthText);
if (wellNameIndex == -1 || surfaceNameIndex == -1 || measuredDepthIndex == -1) if (wellNameIndex != -1 && surfaceNameIndex != -1 && measuredDepthIndex != -1)
{ {
return; do
{
QString line = data.readLine();
QStringList dataLine = line.split(';', QString::KeepEmptyParts);
if (dataLine.size() != header.size()) continue;
bool conversionOk;
double measuredDepth = dataLine[measuredDepthIndex].toDouble(&conversionOk);
if (!conversionOk) continue;
QString wellName = dataLine[wellNameIndex];
QString surfaceName = dataLine[surfaceNameIndex];
(*formations)[wellName].push_back(std::make_pair(measuredDepth, surfaceName));
} while (!data.atEnd());
return;
} }
do static const QString unitNameText = "unitname";
static const QString measuredDepthToptext = "topmd";
static const QString measuredDepthBasetext = "basemd";
int unitNameIndex = header.indexOf(unitNameText);
int measuredDepthTopIndex = header.indexOf(measuredDepthToptext);
int measuredDepthBaseIndex = header.indexOf(measuredDepthBasetext);
if (unitNameIndex != -1 && measuredDepthTopIndex != -1 && measuredDepthBaseIndex != -1)
{ {
QString line = data.readLine(); do
{
QString line = data.readLine();
QStringList dataLine = line.split(';', QString::KeepEmptyParts); QStringList dataLine = line.split(';', QString::KeepEmptyParts);
if (dataLine.size() != header.size()) continue; if (dataLine.size() != header.size()) continue;
bool conversionOk; QString wellName = dataLine[wellNameIndex];
double measuredDepth = dataLine[measuredDepthIndex].toDouble(&conversionOk); QString unitName = dataLine[unitNameIndex];
if (!conversionOk) continue; unitName = unitName.trimmed();
double measuredDepthTop = dataLine[measuredDepthTopIndex].toDouble();
double measuredDepthBase = dataLine[measuredDepthBaseIndex].toDouble();
QString wellName = dataLine[wellNameIndex]; (*formations)[wellName].push_back(std::make_pair(measuredDepthTop, unitName));
QString surfaceName = dataLine[surfaceNameIndex]; (*formations)[wellName].push_back(std::make_pair(measuredDepthBase, unitName));
(*formations)[wellName].push_back(std::make_pair(measuredDepth, surfaceName)); } while (!data.atEnd());
}
} while (!data.atEnd());
} }