OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-18 21:43:27 -06:00
Code Issues Packages Projects Releases Wiki Activity

Update support for VFP tables: Now produces almost identical results to BHP for 'identity' table'

Browse Source
This commit is contained in:
babrodtk 2015-08-06 15:11:11 +02:00
parent 83f6b02941
commit ffac75a873
1 changed files with 23 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

31
opm/core/simulator/WellState.hpp
View File

@ -67,6 +67,7 @@ namespace Opm
bhp_[w] = well_controls_get_current_target( ctrl );
break;
case THP:
assert(false && "Not properly implemented");
thp_[w] = well_controls_get_current_target( ctrl );
break;
default:
@ -97,20 +98,34 @@ namespace Opm
wellrates_[np*w + p] = small_rate * sign;
}
}
// 2. Initialize bhp to be target pressure if
// bhp-controlled well, otherwise set to a
// 2. if we have a thp/bhp control,
// set target (may be overridden later)
// Assumes only one THP / BHP control.
thp_[w] = -1e100;
bhp_[w] = -1e100;
int num_controls = well_controls_get_num( ctrl );
for (int i=0; i<num_controls; ++i) {
switch (well_controls_iget_type( ctrl , i )) {
case BHP:
bhp_[w] = well_controls_iget_target( ctrl , i );
break;
case THP:
thp_[w] = well_controls_iget_target( ctrl , i );
break;
}
}
// 3. Set BHP to a
// little above or below (depending on if
// the well is an injector or producer)
// pressure in first perforation cell.
switch (well_controls_get_current_type(ctrl)) {
case BHP:
bhp_[w] = well_controls_get_current_target( ctrl );
thp_[w] = -1e100;
break;
case THP:
bhp_[w] = -1e100;
thp_[w] = well_controls_get_current_target( ctrl );
//bhp_[w] = thp_[w]; //< TODO: ARB Adding this produces identical results as without THP control for artificial test case
//Already taken care of above in 2.
break;
default:
@ -118,7 +133,7 @@ namespace Opm
const int first_cell = wells->well_cells[wells->well_connpos[w]];
const double safety_factor = (wells->type[w] == INJECTOR) ? 1.01 : 0.99;
bhp_[w] = safety_factor*state.pressure()[first_cell];
thp_[w] = -1e100;
// thp_[w] = -1e100;
}
}
}