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

Fetch connpos and number of perforations from WellInfo for MSW wells

Browse Source
This commit is contained in:
Joakim Hove
2020-12-14 08:45:03 +01:00
parent 38eb0802d6
commit aa345a36f1
1 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
opm/simulators/wells/WellStateFullyImplicitBlackoil.hpp
View File

@@ -722,12 +722,15 @@ namespace Opm
seg_number_.clear(); seg_number_.clear();
nseg_ = 0; nseg_ = 0;
int connpos = 0;
// in the init function, the well rates and perforation rates have been initialized or copied from prevState // in the init function, the well rates and perforation rates have been initialized or copied from prevState
// what we do here, is to set the segment rates and perforation rates // what we do here, is to set the segment rates and perforation rates
for (int w = 0; w < nw; ++w) { for (int w = 0; w < nw; ++w) {
const auto& well_ecl = wells_ecl[w]; const auto& well_ecl = wells_ecl[w];
int num_perf_this_well = well_perf_data_[w].size(); const auto& wname = wells_ecl[w].name();
const auto& well_info = this->wellMap().at(wname);
const int connpos = well_info[1];
const int num_perf_this_well = well_info[2];
top_segment_index_.push_back(nseg_); top_segment_index_.push_back(nseg_);
if ( !well_ecl.isMultiSegment() ) { // not multi-segment well if ( !well_ecl.isMultiSegment() ) { // not multi-segment well
nseg_ += 1; nseg_ += 1;
@@ -748,6 +751,7 @@ namespace Opm
for (auto segID = 0*well_nseg; segID < well_nseg; ++segID) { for (auto segID = 0*well_nseg; segID < well_nseg; ++segID) {
this->seg_number_.push_back(segment_set[segID].segmentNumber()); this->seg_number_.push_back(segment_set[segID].segmentNumber());
} }
// we need to know for each segment, how many perforation it has and how many segments using it as outlet_segment // we need to know for each segment, how many perforation it has and how many segments using it as outlet_segment
// that is why I think we should use a well model to initialize the WellState here // that is why I think we should use a well model to initialize the WellState here
std::vector<std::vector<int>> segment_perforations(well_nseg); std::vector<std::vector<int>> segment_perforations(well_nseg);
@@ -826,7 +830,6 @@ namespace Opm
} }
} }
} }
connpos += num_perf_this_well;
} }
assert(int(seg_press_.size()) == nseg_); assert(int(seg_press_.size()) == nseg_);
assert(int(seg_rates_.size()) == nseg_ * numPhases() ); assert(int(seg_rates_.size()) == nseg_ * numPhases() );