A Graph Network Based Approach for Reservoir Modeling

2021 ◽  
Author(s):  
Wenyue Sun ◽  
Sathish Sankaran

Abstract Reservoir management routinely requires assimilating historical data and predicting field performance against multiple production strategies before implementing them in the field. However, traditional numerical methods are often cumbersome to characterize, build and calibrate at a timescale that can be used reliably for such short-term decision cycles such as production forecasting, IOR optimization and production rate control. Simpler analytical models make assumptions and lack the rigor needed to adequately model these systems. Pure data-driven methods may lack physical insights or have limited range of applicability. Model fidelity, speed, interpretability, suitability to automate and ease-of-use are some key modeling traits that are desired for reservoir management purposes. In this work, we propose to use a reservoir graph-network modeling approach (RGNet), based on the concept of diffusive time of flight, to forecast well performance using routinely measured field measurements (e.g. bottomhole pressure and rates). We propose a novel, model order reduction method based on discretized time of flight for multiple wells with interference. It simplifies the 3D reservoir flow problem into a flow network representation that can be solved as a 2D simulation model with any general-purpose reservoir simulator. Parameters in RGNet model cover well productivity index, grid pore volume and transmissibility, which are estimated through a history-matching process. After history matching, multiple posterior RGNet models are generated to quantify subsurface uncertainties. The RGNet modeling approach allows various fluid-flow physics to be modeled within the grids and boundary conditions, and is applicable to a range of conventional and unconventional reservoirs with different flow mechanisms. We applied the proposed approach on a field case reservoir models for multiple wells with interference. By virtue of the reduced complexity, the modeling methodology is highly scalable and still retains physical interpretability. The calibration method produces parsimonious models and provides uncertainty estimates in history matching parameters with range of outcomes. In addition, the RGNet models are much faster to simulate, over 1000x speed up, compared with full-physics models. We then used RGNet models for well-control and flood optimization and achieved significant improvement over field net-present-values. Parameterization of the proposed reservoir graph-network modeling approach provides a unique and sustainable way to reduce model complexity needed for reservoir management purposes. The method is rooted in physical principles and provides an explainable dynamic reservoir model that can be effectively used to understand reservoir behavior and optimize performance. The lightweight model lends itself naturally to fast computation that are required for scenario analysis and optimization.

2021 ◽  
Author(s):  
Tsubasa Onishi ◽  
Hongquan Chen ◽  
Jiang Xie ◽  
Shusei Tanaka ◽  
Dongjae Kam ◽  
...  

Abstract Streamline-based methods have proven to be effective for various subsurface flow and transport modeling problems. However, the applications are limited in dual-porosity and dual-permeability (DPDK) system due to the difficulty in describing interactions between matrix and fracture during streamline tracing. In this work, we present a robust streamline tracing algorithm for DPDK models and apply the new algorithm to rate allocation optimization in a waterflood reservoir. In the proposed method, streamlines are traced in both fracture and matrix domains. The inter-fluxes between fracture and matrix are described by switching streamlines from one domain to another using a probability computed based on the inter-fluxes. The approach is fundamentally similar to the existing streamline tracing technique and can be utilized in streamline-assisted applications, such as flow diagnostics, history matching, and production optimization. The proposed method is benchmarked with a finite-volume based approach where grid-based time-of-flight was obtained by solving the stationary transport equation. We first validated our method using simple examples. Visual time-of-flight comparisons as well as tracer concentration and allocation factors at wells show good agreement. Next, we applied the proposed method to field scale models to demonstrate the robustness. The results show that our method offers reduced numerical artifacts and better represents reservoir heterogeneity and well connectivity with sub-grid resolutions. The proposed method is then used for rate allocation optimization in DPDK models. A streamline-based gradient free algorithm is used to optimize net present value by adjusting both injection and production well rates under operational constraints. The results show that the optimized schedule offers significant improvement in recovery factor, net present value, and sweep efficiency compared to the base scenario using equal rate injection and production. The optimization algorithm is computationally efficient as it requires only a few forward reservoir simulations.


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 71-91 ◽  
Author(s):  
Salam Al-Rbeawi

Summary The objective of this paper is to revisit currently used techniques for analyzing reservoir performance and characterizing the horizontal-well productivity index (PI) in finite-acting oil and gas reservoirs. This paper introduces a new practical and integrated approach for determining the starting time of pseudosteady-state flow and constant-behavior PI. The new approach focuses on the fact that the derivative of PI vanishes to zero when pseudosteady-state flow is developed. At this point, the derivative of transient-state pressure drop and that of pseudosteady-state pressure drop become mathematically identical. This point indicates the starting time of pseudosteady-state flow as well as the constant value of pseudosteady-state PI. The reservoirs of interest in this study are homogeneous and heterogamous, single and dual porous media, undergoing Darcy and non-Darcy flow in the drainage area, and finite-acting, depleted by horizontal wells. The flow in these reservoirs is either single-phase oil flow or single-phase gas flow. Several analytical models are used in this study for describing pressure and pressure-derivative behavior considering different reservoir configurations and wellbore types. These models are developed for heterogeneous and homogeneous formations consisting of single and dual porous media (naturally fractured reservoirs) and experiencing Darcy and non-Darcy flow. Two pressure terms are assembled in these models; the first pressure term represents the time-dependent pressure drop caused by transient-state flow, and the second pressure term represents time-invariant pressure drop controlled by the reservoir boundary. Transient-state PI and pseudosteady-state PI are calculated using the difference between these two pressures assuming constant wellbore flow rate. The analytical models for the pressure derivatives of these two pressure terms are generated. Using the concept that the derivative of constant PI converges to zero, these two pressure derivatives become mathematically equal at a certain production time. This point indicates the starting time of pseudosteady-state flow and the constant behavior of PI. The outcomes of this study are summarized as the following: Understanding pressure, pressure derivative, and PI behavior of bounded reservoirs drained by horizontal wells during transient- and pseudosteady-state production Investigating the effects of different reservoir configurations, wellbore lengths, reservoir homogeneity or heterogeneity, reservoirs as single or dual porous media, and flow pattern in porous media whether it has undergone Darcy or non-Darcy flow Applying the concept of the PI derivative to determine the starting time of pseudosteady-state stabilized PI The novel points in this study are the following: The derivative of the PI can be used to precisely indicate the starting time of pseudosteady-state flow and the constant behavior of PI. The starting time of pseudosteady-state flow determined by the convergence of transient- and pseudosteady-state pressure derivative or by the PI curve is always less than that determined from the curves of total pressure drop and its derivative. Non-Darcy flow may significantly affect the transient-state PI, but pseudosteady-state PI is slightly affected by non-Darcy flow. The starting time of pseudosteady-state flow is not influenced by non-Darcy flow. The convergence of transient- and pseudosteady-state pressure derivatives is affected by reservoir configurations, wellbore lengths, and porous-media characteristics.


2018 ◽  
Author(s):  
Andreas Plach ◽  
Kerim H. Nisancioglu ◽  
Sébastien Le clec’h ◽  
Andreas Born ◽  
Petra M. Langebroek ◽  
...  

Abstract. Understanding the behavior of the Greenland ice sheet in a warmer climate, and particularly its surface mass balance (SMB), is important for assessing Greenland’s potential contribution to future sea level rise. The Eemian interglacial, the most recent warmer-than-present period in Earth’s history approximately 125 000 years ago, provides an analogue for a warm summer climate over Greenland. The Eemian is characterized by a positive Northern Hemisphere summer insolation anomaly, which introduces uncertainties in Eemian SMB when using positive degree day estimates. In this study, we use Eemian global and regional climate simulations in combination with three types of SMB models – a simple positive degree day, an intermediate complexity, and a full surface energy balance model – to evaluate the importance of regional climate and model complexity for estimates of Greenland SMB. We find that all SMB models perform well under the relatively cool pre-industrial and late Eemian. For the relatively warm early Eemian, the differences between SMB models are large which is associated with the representation of insolation in the respective models. For all simulated time slices there is a systematic difference between globally and regionally forced SMB models, due to the different representation of the regional climate over Greenland. We conclude that both the resolution of the simulated climate as well as the method used to estimate the SMB, are important for an accurate simulation of Greenland’s SMB. Whether model resolution or SMB method is most important depends on the climate state and in particular the prevailing insolation pattern. We suggest that future Eemian climate model inter-comparison studies are combined with different SMB models to quantify Eemian SMB uncertainty estimates.


2021 ◽  
Author(s):  
Evgeniy Viktorovich Yudin ◽  
George Aleksandrovich Piotrovskiy ◽  
Maria Vladimirovna Petrova ◽  
Alexey Petrovich Roshchektaev ◽  
Nikita Vladislavovich Shtrobel

Abstract Requirements of targeted optimization are imposed on the hydraulic fracturing operations carried out in the conditions of borderline economic efficiency of fields taking into account geological and technological features. Consequently, the development of new analytical tools foranalyzing and planning the productivity of fractured wells, taking into account the structuralfeatures of the productive reservoir and inhomogeneous distribution of the fracture conductivity, is becoming highly relevant. The paper proposes a new approach of assessing the vertical hydraulic fracture productivityin a rectangular reservoir in a pseudo-steady state, based on reservoir resistivity concept described in the papers of Meyer et al. However, there is a free parameter in the case of modeling the productivity of a hydraulic fracture by the concept. The parameter describes the distribution of the inflow along the plane of the fracture. This paper presents a systematic approach to determining of the parameter. The resulting model allows to conduct an assessment of the influence of various complications in the fracture on the productivity index. During the research a method of determining the free parameter was developed,it was based on the obtained dependence of the inflow distribution on the coordinate along the fracture of finite conductivity. The methodology allowed to refine existent analytical solution of the Meyer et al. model, which, in turn, allowed to assess the influence of different fracture damages in the hydraulic fracture on the productivity index of the well. The work includes the cases of the presence of fracture damages at the beginning and at the end of the fracture. A hydraulic fracture model was built for each of the types of damages, it was based on the developed method, and also the solution of dimensionless productivity ratio was received. The results of the obtained solution were confirmed by comparison with the numerical solutions of commercial simulators and analytical models available in the literature. The advantage of the methodology is the resulting formulas for well productivity are relatively simple, even for exotic cases ofvariable conductivity fractures. The approaches and algorithms described in the paper assume the calculation of the productivity of a hydraulic fracture with variable conductivity and the presence of other complicatingfactors.The methodology of the paper can be used for analysis and diagnosis problems with formation hydraulic fracturing. The efficiency of the calculations allows using the presented methodology to solve inverse problems of determining the efficiency of the hydraulic fracturing operation.


Author(s):  
William T. Cousins ◽  
Milt W. Davis

Modeling compression systems using parallel compressor theory has been used for the analysis of compression system operability since the 1960s. Parallel compressor models have been traditionally designed and used for the analysis of circumferential distortion effects as a means to evaluate the impact of various inlet flow field disturbances on compressor operation. This paper (the first of two) provides a review of the parallel compressor concept and discusses extensions to the original theory. These extensions include the incorporation of dynamic response, application to complex distortions, and the application to inlet swirl. Understanding these effects and the application of parallel compressor theory extensions is required to produce analytical models and computer simulations that can be used to enhance the development testing and the understanding of the response of gas turbine compression systems. Once a computer simulation has been constructed for a particular test article, it can be exercised and results compared against test results where distortion-generator devices (such as distortion screens) have been used, generally with favorable accuracy. The usefulness of the extended parallel compressor model is derived from its ease of use, simplicity, and ability for quick turn-around of results. It is often more desirable to have an analysis capability that is easy and quick to use than to have one that is extremely accurate, especially when understanding basic physics is of primary concern during a test operation. Extreme accuracy may require large amounts of computer resources and take days or weeks to compute a single performance point. While this may be acceptable for design, the limitations of high-fidelity simulations make them impractical to use due to the time constraints imposed by the pace of testing. Applying a timely analysis capability, using a parallel compressor simulation can provide a new physical understanding of the effects of complex distortion during the testing process when comparing the analytical and test results. This concept is presented in two companion papers: the first paper, Part 1, concentrates on the parallel compressor concepts, theory and limitations of the methodology while the second paper, Part 2 [1] presents applications of the approximate methods developed and compares results with experimental data.


Sign in / Sign up

Export Citation Format

Share Document