OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add Transport Protocol for Segment Level Phase Densities

Browse Source 
This commit adds the requisite backing storage and parser support
for capturing and transporting simulator-level calculation of phase
and mixture density value for purpose of summary file output.  To
this end, make 'SegmentQuantity' into a template on a set of defined
items and make SegmentPhaseQuantity into a specialisation of this
template.  Add a new specialisation, SegmentPhaseDensity, which
holds phase densities (oil, gas, water), and fluid mixture densities
with and without flowing fraction exponents.

These will be used to transport the values needed to output segment
level summary vectors

    SDENx -- Phase density of phase 'x' (O, G, W)
    SDENM -- Fluid mixture density without flowing fraction exponents
    SMDEN -- Fluid mixture density with flowing fraction exponents
This commit is contained in:
Bård Skaflestad 2023-07-04 19:04:45 +02:00
parent 784faa2b05
commit e7226d4ca8
4 changed files with 94 additions and 23 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

100
opm/output/data/Wells.hpp
View File

@ -318,15 +318,11 @@ namespace Opm {
}
};
class SegmentPhaseQuantity
template <typename Items>
class SegmentQuantity
{
public:
enum class Item : std::size_t {
Oil, Gas, Water,
// -- Must be last enumerator --
NumItems,
};
using Item = typename Items::Item;
void clear()
{
@ -341,7 +337,7 @@ namespace Opm {
return (i < Size) && this->hasItem(i);
}
bool operator==(const SegmentPhaseQuantity& vec) const
bool operator==(const SegmentQuantity& vec) const
{
return (this->has_ == vec.has_)
&& (this->value_ == vec.value_);
@ -351,21 +347,21 @@ namespace Opm {
{
if (! this->has(p)) {
throw std::invalid_argument {
"Request for Unset Item Value for " + this->itemName(p)
"Request for Unset Item Value for " + Items::itemName(p)
};
}
return this->value_[ this->index(p) ];
}
SegmentPhaseQuantity& set(const Item p, const double value)
SegmentQuantity& set(const Item p, const double value)
{
const auto i = this->index(p);
if (i >= Size) {
throw std::invalid_argument {
"Cannot Assign Item Value for Unsupported Item '"
+ this->itemName(p) + '\''
+ Items::itemName(p) + '\''
};
}
@ -403,12 +399,15 @@ namespace Opm {
serializer(this->value_);
}
static SegmentPhaseQuantity serializationTestObject()
static SegmentQuantity serializationTestObject()
{
return SegmentPhaseQuantity{}
.set(Item::Oil , 1.0)
.set(Item::Gas , 7.0)
.set(Item::Water, 2.9);
auto quant = SegmentQuantity{};
for (const auto& [item, value] : Items::serializationTestItems()) {
quant.set(item, value);
}
return quant;
}
private:
@ -430,8 +429,18 @@ namespace Opm {
{
return (this->has_ & (1 << i)) != 0;
}
};
std::string itemName(const Item p) const
struct PhaseItems
{
enum class Item {
Oil, Gas, Water,
// -- Must be last enumerator --
NumItems,
};
static std::string itemName(const Item p)
{
switch (p) {
case Item::Oil: return "Oil";
@ -442,10 +451,59 @@ namespace Opm {
return "Out of bounds (NumItems)";
}
return "Unknown (" + std::to_string(this->index(p)) + ')';
return "Unknown (" + std::to_string(static_cast<int>(p)) + ')';
}
static auto serializationTestItems()
{
return std::vector {
std::pair { Item::Oil , 1.0 },
std::pair { Item::Gas , 7.0 },
std::pair { Item::Water, 2.9 },
};
}
};
struct DensityItems
{
enum class Item {
Oil, Gas, Water, Mixture, MixtureWithExponents,
// -- Must be last enumerator --
NumItems,
};
static std::string itemName(const Item p)
{
switch (p) {
case Item::Oil: return "Oil";
case Item::Gas: return "Gas";
case Item::Water: return "Water";
case Item::Mixture: return "Mixture";
case Item::MixtureWithExponents: return "MixtureWithExponents";
case Item::NumItems:
return "Out of bounds (NumItems)";
}
return "Unknown (" + std::to_string(static_cast<int>(p)) + ')';
}
static auto serializationTestItems()
{
return std::vector {
std::pair { Item::Oil , 876.54 },
std::pair { Item::Gas , 321.09 },
std::pair { Item::Water , 987.65 },
std::pair { Item::Mixture , 975.31 },
std::pair { Item::MixtureWithExponents, 765.43 },
};
}
};
using SegmentPhaseQuantity = SegmentQuantity<PhaseItems>;
using SegmentPhaseDensity = SegmentQuantity<DensityItems>;
struct Segment
{
Rates rates{};
@ -453,6 +511,7 @@ namespace Opm {
SegmentPhaseQuantity velocity{};
SegmentPhaseQuantity holdup{};
SegmentPhaseQuantity viscosity{};
SegmentPhaseDensity density{};
std::size_t segNumber{};
bool operator==(const Segment& seg2) const
@ -462,6 +521,7 @@ namespace Opm {
&& (velocity == seg2.velocity)
&& (holdup == seg2.holdup)
&& (viscosity == seg2.viscosity)
&& (density == seg2.density)
&& (segNumber == seg2.segNumber);
}
@ -479,6 +539,7 @@ namespace Opm {
serializer(this->velocity);
serializer(this->holdup);
serializer(this->viscosity);
serializer(this->density);
serializer(this->segNumber);
}
@ -490,6 +551,7 @@ namespace Opm {
SegmentPhaseQuantity::serializationTestObject(), // velocity
SegmentPhaseQuantity::serializationTestObject(), // holdup
SegmentPhaseQuantity::serializationTestObject(), // viscosity
SegmentPhaseDensity::serializationTestObject(), // density
10
};
}
@ -1044,6 +1106,7 @@ namespace Opm {
this->velocity.write(buffer);
this->holdup.write(buffer);
this->viscosity.write(buffer);
this->density.write(buffer);
}
template <class MessageBufferType>
@ -1167,6 +1230,7 @@ namespace Opm {
this->velocity.read(buffer);
this->holdup.read(buffer);
this->viscosity.read(buffer);
this->density.read(buffer);
}
template <class MessageBufferType>

1
src/opm/input/eclipse/EclipseState/SummaryConfig/SummaryConfig.cpp
View File

@ -1150,6 +1150,7 @@ inline void keywordMISC( SummaryConfig::keyword_list& list,
"SGFR" , "SGFRF", "SGFRS", "SGFT", "SGFV", "SGHF", "SGVIS",
"SWFR" , "SWFT", "SWFV", "SWHF", "SWVIS",
"SGOR" , "SOGR" , "SWCT" , "SWGR" ,
"SODEN", "SGDEN", "SWDEN", "SMDEN", "SDENM",
"SPR" , "SPRD" , "SPRDH", "SPRDF", "SPRDA",
};

5
src/opm/input/eclipse/share/keywords/000_Eclipse100/S/SEGMENT_PROBE
View File

@ -4,6 +4,9 @@
"SUMMARY"
],
"deck_names": [
"SDENM",
"SMDEN",
"SGDEN",
"SGFR",
"SGFRF",
"SGFRS",
@ -12,6 +15,7 @@
"SGHF",
"SGOR",
"SGVIS",
"SODEN",
"SOFR",
"SOFRF",
"SOFRS",
@ -26,6 +30,7 @@
"SPRDF",
"SPRDH",
"SWCT",
"SWDEN",
"SWFR",
"SWFT",
"SWFV",

11
src/opm/output/eclipse/Summary.cpp
View File

@ -1058,8 +1058,9 @@ inline quantity srate(const fn_args& args)
});
}
template <Opm::data::SegmentPhaseQuantity::Item p>
double segment_phase_quantity_value(const Opm::data::SegmentPhaseQuantity& q)
template <typename Items>
double segment_phase_quantity_value(const Opm::data::SegmentQuantity<Items>& q,
const typename Items::Item p)
{
return q.has(p) ? q.get(p) : 0.0;
}
@ -1071,7 +1072,7 @@ quantity segment_flow_velocity(const fn_args& args)
[](const Opm::data::Segment& segment)
{
// Note: Opposite velocity sign conventions in Flow vs. ECLIPSE.
return - segment_phase_quantity_value<p>(segment.velocity);
return - segment_phase_quantity_value(segment.velocity, p);
});
}
@ -1081,7 +1082,7 @@ quantity segment_holdup_fraction(const fn_args& args)
return segment_quantity(args, measure::identity,
[](const Opm::data::Segment& segment)
{
return segment_phase_quantity_value<p>(segment.holdup);
return segment_phase_quantity_value(segment.holdup, p);
});
}
@ -1091,7 +1092,7 @@ quantity segment_viscosity(const fn_args& args)
return segment_quantity(args, measure::viscosity,
[](const Opm::data::Segment& segment)
{
return segment_phase_quantity_value<p>(segment.viscosity);
return segment_phase_quantity_value(segment.viscosity, p);
});
}