A Physics-Based Proxy for Surface and Subsurface Coupled Simulation Models

2021 ◽  
Author(s):  
Changdong Yang ◽  
Jincong He ◽  
Song Du ◽  
Zhenzhen Wang ◽  
Tsubasa Onishi ◽  
...  

Abstract Full-physics subsurface simulation models coupled with surface network can be computationally expensive. In this paper, we propose a physics-based subsurface model proxy that significantly reduces the run-time of the coupled model to enable rapid decision-making for reservoir management. In the coupled model the subsurface reservoir simulator generates well inflow performance relationship (IPR) curves which are used by the surface network model to determine well rates that satisfy surface constraints. In the proposed proxy model, the CPU intensive reservoir simulation is replaced with an IPR database constructed from a data pool of one or multiple simulation runs. The IPR database captures well performance that represents subsurface reservoir dynamics. The proxy model can then be used to predict the production performance of new scenarios – for example new drilling sequence – by intelligently looking up the appropriate IPR curves for oil, gas and water phases for each well and solving it with the surface network. All necessary operational events in the surface network and field management logic (such as facility constraints, well conditional shut-in, and group guide rate balancing) for the full-coupled model can be implemented and honored. In the proposed proxy model, while the reservoir simulation component is eliminated for efficiency. The entirety of the surface network model is retained, which offers certain advantages. It is particularly suitable for investigating the impact of different surface operations, such as maintenance schedule and production routing changes, with the aim of minimizing production capacity off-line due to maintenance. Replacing the computationally intensive subsurface simulation with the appropriate IPR significantly improves the run time of the coupled model while preserving the essential physics of the reservoir. The accuracy depends on the difference between the scenarios that the proxy is trained on and the scenarios being evaluated. Initial testing with a complex reservoir with more than 300 wells showed the accuracy of the proxy model to be more than 95%. The computation speedup could be an order of magnitude, depending largely on complexity of the surface network model. Prior work exists in the literature that uses decline curves to replicate subsurface model performance. The use of the multi-phase IPR database and the intelligent lookup mechanism in the proposed method allows it to be more accurate and flexible in handling complexities such as multi-phase flow and interference in the surface network.

2018 ◽  
Vol 140 (12) ◽  
Author(s):  
Jingwen Zheng ◽  
Juliana Y. Leung ◽  
Ronald P. Sawatzky ◽  
Jose M. Alvarez

Artificial intelligence (AI) tools are used to explore the influence of shale barriers on steam-assisted gravity drainage (SAGD) production. The data are derived from synthetic SAGD reservoir simulations based on petrophysical properties and operational constraints gathered from the Suncor's Firebag project, which is representative of Athabasca oil sands reservoirs. The underlying reservoir simulation model is homogeneous and two-dimensional. Reservoir heterogeneities are modeled by superimposing sets of idealized shale barrier configurations on this homogeneous reservoir model. The individual shale barriers are categorized by their location relative to the SAGD well pair and by their geometry. SAGD production for a training set of shale barrier configurations was simulated. A network model based on AI tools was constructed to match the output of the reservoir simulation for this training set of shale barrier configurations, with a focus on the production rate and the steam-oil ratio (SOR). Then the trained AI proxy model was used to predict SAGD production profiles for arbitrary configurations of shale barriers. The predicted results were consistent with the results of the SAGD simulation model with the same shale barrier configurations. The results of this work demonstrate the capability and flexibility of the AI-based network model, and of the parametrization technique for representing the characteristics of the shale barriers, in capturing the effects of complex heterogeneities on SAGD production. It offers the significant potential of providing an indirect method for inferring the presence and distribution of heterogeneous reservoir features from SAGD field production data.


Author(s):  
Stein Krogstad ◽  
Halvor Møll Nilsen

AbstractModel-based optimization of placement and trajectories of wells in petroleum reservoirs by the means of reservoir simulation forecasts is computationally demanding due to the high number of simulations typically required to achieve a local optimum. In this work, we develop an efficient flow-diagnostics proxy for net-present-value (NPV) with adjoint capabilities for efficient computation of well control gradients and approximate sensitivities with respect to placement/trajectory parameters. The suggested flow-diagnostic proxy consists of numerically solving a single pressure equation for the given scenario and the solution of a few inter-well time-of-flight and steady-state tracer equations, typically achieved in a few seconds for a reservoir model of medium size. Although the proxy may not be a particularly good approximation for the full reservoir simulation response, we find that for the cases considered, the correlation is very good and hence the proxy is suitable for use in an optimization loop. The adjoint simulation for the proxy model which provides control gradients and placement sensitivities is of similar computational complexity as the forward proxy model (a few seconds). We employ a version of the generalized reduced gradient for handling individual well constraints (e.g., bottom-hole-pressures and rates). As a result, the individual well constraints are enforced within the flow-diagnostics computations, and hence every parameter update becomes feasible without sacrificing gradient information. We present two numerical experiments illustrating the efficiency and performance of the approach for well placement problems involving trajectories and simulation models of realistic complexity. The suggested placements are evaluated using full simulations. We conclude by discussing limitations and possible enhancements of the methodology.


Author(s):  
Klaus Rollmann ◽  
Aurea Soriano-Vargas ◽  
Forlan Almeida ◽  
Alessandra Davolio ◽  
Denis Jose Schiozer ◽  
...  

2021 ◽  
Author(s):  
Mokhles Mezghani ◽  
Mustafa AlIbrahim ◽  
Majdi Baddourah

Abstract Reservoir simulation is a key tool for predicting the dynamic behavior of the reservoir and optimizing its development. Fine scale CPU demanding simulation grids are necessary to improve the accuracy of the simulation results. We propose a hybrid modeling approach to minimize the weight of the full physics model by dynamically building and updating an artificial intelligence (AI) based model. The AI model can be used to quickly mimic the full physics (FP) model. The methodology that we propose consists of starting with running the FP model, an associated AI model is systematically updated using the newly performed FP runs. Once the mismatch between the two models is below a predefined cutoff the FP model is switch off and only the AI model is used. The FP model is switched on at the end of the exercise either to confirm the AI model decision and stop the study or to reject this decision (high mismatch between FP and AI model) and upgrade the AI model. The proposed workflow was applied to a synthetic reservoir model, where the objective is to match the average reservoir pressure. For this study, to better account for reservoir heterogeneity, fine scale simulation grid (approximately 50 million cells) is necessary to improve the accuracy of the reservoir simulation results. Reservoir simulation using FP model and 1024 CPUs requires approximately 14 hours. During this history matching exercise, six parameters have been selected to be part of the optimization loop. Therefore, a Latin Hypercube Sampling (LHS) using seven FP runs is used to initiate the hybrid approach and build the first AI model. During history matching, only the AI model is used. At the convergence of the optimization loop, a final FP model run is performed either to confirm the convergence for the FP model or to re iterate the same approach starting from the LHS around the converged solution. The following AI model will be updated using all the FP simulations done in the study. This approach allows the achievement of the history matching with very acceptable quality match, however with much less computational resources and CPU time. CPU intensive, multimillion-cell simulation models are commonly utilized in reservoir development. Completing a reservoir study in acceptable timeframe is a real challenge for such a situation. The development of new concepts/techniques is a real need to successfully complete a reservoir study. The hybrid approach that we are proposing is showing very promising results to handle such a challenge.


2021 ◽  
Author(s):  
Obinna Somadina Ezeaneche ◽  
Robinson Osita Madu ◽  
Ishioma Bridget Oshilike ◽  
Orrelo Jerry Athoja ◽  
Mike Obi Onyekonwu

Abstract Proper understanding of reservoir producing mechanism forms a backbone for optimal fluid recovery in any reservoir. Such an understanding is usually fostered by a detailed petrophysical evaluation, structural interpretation, geological description and modelling as well as production performance assessment prior to history matching and reservoir simulation. In this study, gravity drainage mechanism was identified as the primary force for production in reservoir X located in Niger Delta province and this required proper model calibration using variation of vertical anisotropic ratio based on identified facies as against a single value method which does not capture heterogeneity properly. Using structural maps generated from interpretation of seismic data, and other petrophysical parameters from available well logs and core data such as porosity, permeability and facies description based on environment of deposition, a geological model capturing the structural dips, facies distribution and well locations was built. Dynamic modeling was conducted on the base case model and also on the low and high case conceptual models to capture different structural dips of the reservoir. The result from history matching of the base case model reveals that variation of vertical anisotropic ratio (i.e. kv/kh) based on identified facies across the system is more effective in capturing heterogeneity than using a deterministic value that is more popular. In addition, gas segregated fastest in the high case model with the steepest dip compared to the base and low case models. An improved dynamic model saturation match was achieved in line with the geological description and the observed reservoir performance. Quick wins scenarios were identified and this led to an additional reserve yield of over 1MMSTB. Therefore, structural control, facies type, reservoir thickness and nature of oil volatility are key forces driving the gravity drainage mechanism.


Sign in / Sign up

Export Citation Format

Share Document