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

Tracer summary curves, collect well rates.

Browse Source
This commit is contained in:
Ove Sævareid 2021-05-11 13:01:45 +02:00
parent 969fc20154
commit 6d76f6c557
2 changed files with 49 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

53
ebos/ecltracermodel.hh
View File

@ -103,6 +103,28 @@ public:
gasPhaseIdx, oilPhaseIdx, waterPhaseIdx);
}
/*!
* \brief Collect well tracer rates for reporting
*/
void addTracerRatesToWells(Opm::data::WellRates& wellDatas) const
{
if (numTracers()==0)
return; // no tracers today
for (const auto& wTR : wellTracerRate_) {
std::string wellName = wTR.first.first;
std::string tracerName = tracerNames_[wTR.first.second];
double rate = wTR.second;
if (!wellDatas.count(wellName)) {
Opm::data::Well wellData;
wellDatas.emplace(std::make_pair(wellName, wellData));
}
Opm::data::Well& wellData = wellDatas.at(wellName);
wellData.rates.set(Opm::data::Rates::opt::tracer, rate, tracerName);
}
}
void beginTimeStep()
{
if (this->numTracers()==0)
@ -285,6 +307,8 @@ protected:
const auto& wells = simulator_.vanguard().schedule().getWells(episodeIdx);
for (const auto& well : wells) {
wellTracerRate_[std::make_pair(well.name(),tracerIdx)] = 0.0;
if (well.getStatus() == Well::Status::SHUT)
continue;
@ -299,17 +323,34 @@ protected:
cartesianCoordinate[1] = connection.getJ();
cartesianCoordinate[2] = connection.getK();
const size_t cartIdx = simulator_.vanguard().cartesianIndex(cartesianCoordinate);
const int I = this->cartToGlobal_[cartIdx];
Scalar rate = simulator_.problem().wellModel().well(well.name())->volumetricSurfaceRateForConnection(I, this->tracerPhaseIdx_[tracerIdx]);
if (rate > 0)
this->tracerResidual_[I][0] -= rate*wtracer;
else if (rate < 0)
this->tracerResidual_[I][0] -= rate*this->tracerConcentration_[tracerIdx][I];
const int I = cartToGlobal_[cartIdx];
Scalar rate = simulator_.problem().wellModel().well(well.name())->volumetricSurfaceRateForConnection(I, tracerPhaseIdx_[tracerIdx]);
if (rate > 0) {
tracerResidual_[I][0] -= rate*wtracer;
wellTracerRate_.at(std::make_pair(well.name(),tracerIdx)) += rate*wtracer;
}
else if (rate < 0) {
tracerResidual_[I][0] -= rate*tracerConcentration_[tracerIdx][I];
wellTracerRate_.at(std::make_pair(well.name(),tracerIdx)) += rate*tracerConcentration_[tracerIdx][I];
}
}
}
}
Simulator& simulator_;
std::vector<std::string> tracerNames_;
std::vector<int> tracerPhaseIdx_;
std::vector<Dune::BlockVector<Dune::FieldVector<Scalar, 1>>> tracerConcentration_;
std::vector<Dune::BlockVector<Dune::FieldVector<Scalar, 1>>> tracerConcentrationInitial_;
TracerMatrix *tracerMatrix_;
TracerVector tracerResidual_;
std::vector<int> cartToGlobal_;
std::vector<Dune::BlockVector<Dune::FieldVector<Scalar, 1>>> storageOfTimeIndex1_;
// <wellName, tracerIdx> -> wellRate
std::map<std::pair<std::string, int>, double> wellTracerRate_;
};
} // namespace Opm

2
ebos/eclwriter.hh
View File

@ -170,6 +170,8 @@ public:
simulator_.vanguard().externalSetupTime();
const auto localWellData = simulator_.problem().wellModel().wellData();
simulator_.problem().tracerModel().addTracerRatesToWells(const_cast<Opm::data::WellRates&>(localWellData));
const auto localGroupAndNetworkData = simulator_.problem().wellModel()
.groupAndNetworkData(reportStepNum);