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

Allow multi-segment wells with no grid connections - warning is already issued.

Browse Source
This commit is contained in:
Vegard Kippe
2023-09-25 13:28:08 +02:00
parent 53301b78d4
commit 776597e00b
1 changed files with 19 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

26
src/opm/input/eclipse/Schedule/Schedule.cpp
View File

@@ -529,7 +529,9 @@ namespace
/// \brief Check whether each MS well has COMPSEGS entry andissue error if not.
/// \param welsegs All wells with a WELSEGS entry together with the location.
/// \param compegs All wells with a COMPSEGS entry
void check_compsegs_consistency(::Opm::Schedule::WelSegsSet& welsegs, std::set<std::string>& compsegs)
void check_compsegs_consistency(::Opm::Schedule::WelSegsSet& welsegs,
std::set<std::string>& compsegs,
const std::vector<::Opm::Well>& wells)
{
std::vector<std::pair<std::string,::Opm::KeywordLocation>> difference;
difference.reserve(welsegs.size());
@@ -537,17 +539,27 @@ void check_compsegs_consistency(::Opm::Schedule::WelSegsSet& welsegs, std::set<s
compsegs.begin(), compsegs.end(),
std::back_inserter(difference),
::Opm::Schedule::PairComp());
// Ignore wells without connections
const auto empty_conn = [wells](const std::pair<std::string,::Opm::KeywordLocation> &x) -> bool {
const auto& wname = x.first;
const auto it = std::find_if(wells.begin(), wells.end(),
[wname](const ::Opm::Well& well) -> bool {
return well.name()==wname && well.getConnections().empty();});
return it != wells.end();
};
difference.erase(std::remove_if(difference.begin(), difference.end(), empty_conn), difference.end());
if (difference.size()) {
std::string wells = "well";
std::string well_str = "well";
if (difference.size()>1) {
wells.append("s");
well_str.append("s");
}
wells.append(":");
well_str.append(":");
for(const auto& [name, location] : difference) {
wells.append(fmt::format("\n {} in {} at line {}", name, location.filename, location.lineno));
well_str.append(fmt::format("\n {} in {} at line {}", name, location.filename, location.lineno));
}
auto msg = fmt::format("Missing COMPSEGS keyword for the following multisegment {}.", wells);
auto msg = fmt::format("Missing COMPSEGS keyword for the following multisegment {}.", well_str);
throw Opm::OpmInputError(msg, std::get<1>(difference[0]));
}
}
@@ -697,7 +709,7 @@ void Schedule::iterateScheduleSection(std::size_t load_start, std::size_t load_e
keyword_index++;
}
check_compsegs_consistency(welsegs_wells, compsegs_wells);
check_compsegs_consistency(welsegs_wells, compsegs_wells, this->getWells(report_step));
this->applyGlobalWPIMULT(wpimult_global_factor);
this->end_report(report_step);