OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

WellState FIBO: Return Segment Vectors from report()

Browse Source 
This commit extends class WellStateFullyImplictBlackoil to report
segment-related quantities as Opm::data::Segment objects (included
in Opm::data::WellRates objects).  All wells have at least a top
segment in the context of WellState FIBO, so there is a meaningful
value to report for each well.

We put the extraction of segment-related quantities into a new
helper function

    WellStateFullyImplicitBlackoil::reportSegmentResults()

to avoid cluttering up the body of report() more than absolutely
needed.

The primary use-case for this is assigning appropriate values to
items 8 through 11 of restart vector RSEG.  In turn, this will
enable restoring these quantities from a restart file.
This commit is contained in:
Bård Skaflestad
2018-10-09 00:22:18 +02:00
parent 519e6ad174
commit 2418df701f
3 changed files with 335 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

49
opm/autodiff/WellStateFullyImplicitBlackoil.hpp
View File

@@ -579,6 +579,12 @@ namespace Opm
}
}
assert(local_comp_index == this->wells_->well_connpos[ w + 1 ] - this->wells_->well_connpos[ w ]);
const auto nSeg = this->numSegments(w);
for (auto segID = 0*nSeg; segID < nSeg; ++segID) {
well.segments[segID + 1] =
this->reportSegmentResults(pu, w, segID);
}
}
return res;
@@ -935,6 +941,49 @@ namespace Opm
// Well potentials
std::vector<double> well_potentials_;
::Opm::data::Segment
reportSegmentResults(const PhaseUsage& pu,
const int wellID,
const int segmentID) const
{
auto segRes = ::Opm::data::Segment{};
const auto segDoF =
this->topSegmentIndex(wellID) + segmentID;
const auto* rate =
&this->segRates()[segDoF * this->numPhases()];
segRes.pressure = this->segPress()[segDoF];
if (pu.phase_used[Water]) {
segRes.rates.set(data::Rates::opt::wat,
rate[pu.phase_pos[Water]]);
}
if (pu.phase_used[Oil]) {
segRes.rates.set(data::Rates::opt::oil,
rate[pu.phase_pos[Oil]]);
}
if (pu.phase_used[Gas]) {
segRes.rates.set(data::Rates::opt::gas,
rate[pu.phase_pos[Gas]]);
}
segRes.segNumber = segmentID + 1;
return segRes;
}
int numSegments(const int wellID) const
{
const auto topSeg = this->topSegmentIndex(wellID);
return (wellID + 1 == this->numWells()) // Last well?
? (this->numSegment() - topSeg)
: (this->topSegmentIndex(wellID + 1) - topSeg);
}
};
} // namespace Opm