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

deleting the old unused computeSegmentDensities()

Browse Source 
which is a first version, whose implementation should not be used
anymore.
This commit is contained in:
Kai Bao 2015-10-13 14:31:25 +02:00
parent 617e4c1786
commit b7b5672c86
2 changed files with 0 additions and 89 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/autodiff/BlackoilMultiSegmentModel.hpp
View File

@ -261,12 +261,6 @@ namespace Opm {
std::vector<ADB>& vars,
SolutionState& state) const;
void
computeSegmentDensities(const SolutionState& state,
const std::vector<ADB>& cq_s,
const std::vector<ADB>& b_perfcells,
const WellState& xw);
void
computeSegmentDensities(const SolutionState& state,
const WellState& xw); //,

83
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -1635,89 +1635,6 @@ namespace Opm {
template <class Grid>
void
BlackoilMultiSegmentModel<Grid>::computeSegmentDensities(const SolutionState& /* state */,
const std::vector<ADB>& cq_s,
const std::vector<ADB>& b_perf,
const WellState& xw)
{
// how to calculate the formation valume factor for the segments is the most tricky part.
// ave_rho = (\sigma(rho_surf * Q_surf)) / (\sigma(rho_surf * B))
// how to calculate the rs and rv to obtain B
// should we relate the segments with different cells
// TODO: as the first solution, we calculate the rs and rv as the average rs and rv from the related perforations.
// Based on the current framework, it has to be done one well by one well
const int nw = xw.numWells();
const int nseg_total = xw.numSegments();
const int np = numPhases();
assert(np == int(b_perf.size()));
well_segment_densities_ = ADB::constant(V::Zero(nseg_total)); // initialize to be zero
int start_segment = 0;
int start_perforation = 0;
// surface density
std::vector<double> surf_dens(fluid_.surfaceDensity(), fluid_.surfaceDensity() + np);
for (int w = 0; w < nw; ++w) {
// TODO: to check to make sure b_perf is B or 1/B
// TODO: now for non-segmented wells, we do not calculate the segment densities.
// If the WELSPECS set the AVG, the density for the wellbore can be same with the segment density.
// TODO: which is one way to unify the implementation
WellMultiSegmentConstPtr well = wellsMultiSegment()[w];
const int nseg = well->numberOfSegments();
const int nperf = well->numberOfPerforations();
if (well->isMultiSegmented()) {
// calculate the the formation volume factor.
// calcuate the rs and rv first
// at the moment, the rs and rv are the averaged rs and rv for the perforation cells related to the segment
std::vector<ADB> b_perf_well(np, ADB::null());
std::vector<ADB> cq_s_well(np, ADB::null());
std::vector<ADB> b_segment(np, ADB::null());
// accumlated phase rates for segments from the perforations
std::vector<ADB> q_ps(np, ADB::null());
std::vector<ADB> b_seg(np, ADB::null());
for (int p = 0; p < np; ++p) {
cq_s_well[p] = subset(cq_s[p], Span(nperf, 1, start_perforation));
b_perf_well[p] = subset(b_perf[p], Span(nperf, 1, start_perforation));
q_ps[p] = well->wellOps().p2s_gather * cq_s_well[p];
b_seg[p] = well->wellOps().p2s_average * b_perf_well[p];
}
ADB density = ADB::constant(V::Zero(nseg));
ADB temp_rho_q = ADB::constant(V::Zero(nseg));
ADB temp_volume_rate = ADB::constant(V::Zero(nseg));
for (int p = 0; p < np; ++p) {
temp_rho_q += q_ps[p] * surf_dens[p];
temp_volume_rate += q_ps[p] * b_seg[p];
}
density = temp_rho_q / temp_volume_rate;
well_segment_densities_ += superset(density, Span(nseg, 1, start_segment), nseg_total);
// formation factor volume B should be calculated based on the segment pressures
// how to handle rs and rv to calculate B
// ADB& seg_press = subset(state.segp, Span(nseg, 1, start_segment));
}
start_segment += nseg;
start_perforation += nperf;
}
}
template <class Grid>
void
BlackoilMultiSegmentModel<Grid>::computeSegmentDensities(const SolutionState& state,