Externally coupled workflows that rely on exchanging Inflow Performance Relationships (IPRs) at the coupling points, such as those between reservoir and surface network simulators, may exhibit instabilities due to the fact that the IPR calculated at the beginning of a timestep is not representative of the IPR at the end of the timestep. One solution is to have an implicitly coupled reservoir-surface system. This is often impractical because reservoir and surface networks are modeled using different applications, but also in cases where the resulting coupled system becomes too large and complex to solve implicitly due to processing time and convergence issues. We propose the calculation of multipoint, nonlinear IPRs by solving near-well subdomains for the subsequent timestep. A flexible reservoir simulation architecture enables the dynamic creation and simulation of near-well subdomains at run time. Subdomains are automatically created within the vicinity of the well, or defined dynamically from the pressure gradient. These near-well subdomain simulations are embedded within the full-field simulation, and extract all the required model properties (PVT, rock, etc.) from the full-field model. The most recent fluxes from the global solution are used as boundary conditions of the near-well subdomains. The IPRs created through subdomain solutions agree with those generated by carrying out the full-field simulation multiple times. In this paper, the subdomain IPRs are used within reservoir-network coupling workflows in which regular IPRs result in oscillations and high errors. Sensitivity analysis is carried out on the extent of the subdomains and the size of the coupling timestep. It is shown on a real field case that the subdomain IPRs result in stable coupling for the otherwise unstable model and also help in reducing the balancing errors.

You can access this article if you purchase or spend a download.