Bagheri, Mohammadreza (Research Centre for Fluid and Complex Systems, Coventry University) | Shariatipour, Seyed M. (Research Centre for Fluid and Complex Systems, Coventry University) | Ganjian, Eshmaiel (School of Energy, Construction and Environment, Built & Natural Environment Research Centre, Coventry University)
The fluid pressure, the stress due to the column of the cement in the annulus of oil and gas wells, and the radial pressure exerted on the cement sheath from the surrounding geological layers all affect the integrity of the cement sheath. This paper studies the impact of CO2-bearing fluids, coupled with the geomechanical alterations within the cement matrix on its integrity. These geochemical and geomechanical alterations within the cement matrix have been coupled to determine the cement lifespan. Two main scenarios including radial cracking and radial compaction, were assumed in order to investigate the behaviour of the cement matrix exposed to CO2-bearing fluids over long periods. If the radial pressure from the surrounding rocks on the cement matrix overcomes the strength of the degraded layers within the cement matrix, cement failure can be postponed, while on the other hand, high vertical stress on the cement matrix in the absence of a proper radial pressure can lead to a reduction in the cement lifespan. The radial cracking process generates local areas of high permeability around the outer face of the cement sheath. Our simulation results show at the shallower depths the cement matrices resist CO2-bearing fluids more and this delays exponentially the travel time of CO2-bearing fluids towards the Earth's surface. This is based on the evolution of CO2 gas from the aqueous phase due to the reduction in the fluid pressure at shallower depths, and consumption of CO2 in the reactions which occur at the deeper locations.
Schumi, Bettina (OMV E&P) | Clemens, Torsten (OMV E&P) | Wegner, Jonas (HOT Microfluidics) | Ganzer, Leonhard (Clausthal University of Technology) | Kaiser, Anton (Clariant) | Hincapie, Rafael E. (OMV E&P) | Leitenmüller, Verena (Montan University Leoben)
Chemical Enhanced Oil Recovery leads to substantial incremental costs over waterflooding of oil reservoirs. Reservoirs containing oil with a high Total Acid Number (TAN) could be produced by injection of alkali. Alkali might lead to generation of soaps and emulsify the oil. However, the generated emulsions are not always stable.
Phase experiments are used to determine the initial amount of emulsions generated and their stability if measured over time. Based on the phase experiments, the minimum concentration of alkali can be determined and the concentration of alkali above which no significant increase in formation of initial emulsions is observed.
Micro-model experiments are performed to investigate the effects on pore scale. For injection of alkali into high TAN number oils, mobilization of residual oil after waterflooding is seen. The oil mobilization is due to breaking-up of oil ganglia or movement of elongated ganglia through the porous medium. As the oil is depleting in surface active components, residual oil saturation is left behind either as isolated ganglia or in down-gradient of grains.
Simultaneous injection of alkali and polymers leads to higher incremental oil production in the micro-models owing to larger pressure drops over the oil ganglia and more effective mobilization accordingly.
Core flood tests confirm the micro-model experiments and additional data are derived from these tests. Alkali co-solvent polymer injection leads to the highest incremental oil recovery of the chemical agents which is difficult to differentiate in micro-model experiments. The polymer adsorption is substantially reduced if alkali is injected with polymers compared with polymer injection only. The reason is the effect of the pH on the polymers. As in the micro-models, the incremental oil recovery is also higher for alkali polymer injection than with alkali injection only.
To evaluate the incremental operating costs of the chemical agents, Equivalent Utility Factors (EqUF) are calculated. The EqUF takes the costs of the various chemicals into account. The lowest EqUF and hence lowest chemical incremental OPEX are incurred by injection of Na2CO3, however, the highest incremental recovery factor is seen with alkali co-solvent polymer injection. It should be noted that the incremental oil recovery owing to macroscopic sweep efficiency improvement by polymer needs to be taken into account to assess the efficiency of the chemical agents.
Recent studies have indicated that Huff-n-Puff (HNP) gas injection has the potential to recover an additional 30-70% oil from multi-fractured horizontal wells in shale reservoirs. Nonetheless, this technique is very sensitive to production constraints and is impacted by uncertainty related to measurement quality (particularly frequency and resolution), and lack of constraining data. In this paper, a Bayesian workflow is provided to optimize the HNP process under uncertainty using a Duvernay shale well as an example.
Compositional simulations are conducted which incorporate a tuned PVT model and a set of measured cyclic injection/compaction pressure-sensitive permeability data. Markov chain Monte Carlo (McMC) is used to estimate the posterior distributions of the model uncertain variables by matching the primary production data. The McMC process is accelerated by employing an accurate proxy model (kriging) which is updated using a highly adaptive sampling algorithm. Gaussian Processes are then used to optimize the HNP control variables by maximizing the lower confidence interval (μ-σ) of cumulative oil production (after 10 years) across a fixed ensemble of uncertain variables sampled from posterior distributions.
The uncertain variable space includes several parameters representing reservoir and fracture properties. The posterior distributions for some parameters, such as primary fracture permeability and effective half-length, are narrower, while wider distributions are obtained for other parameters. The results indicate that the impact of uncertain variables on HNP performance is nonlinear. Some uncertain variables (such as molecular diffusion) that do not show strong sensitivity during the primary production strongly impact gas injection HNP performance. The results of optimization under uncertainty confirm that the lower confidence interval of cumulative oil production can be maximized by an injection time of around 1.5 months, a production time of around 2.5 months, and very short soaking times. In addition, a maximum injection rate and a flowing bottomhole pressure around the bubble point are required to ensure maximum incremental recovery. Analysis of the objective function surface highlights some other sets of production constraints with competitive results. Finally, the optimal set of production constraints, in combination with an ensemble of uncertain variables, results in a median HNP cumulative oil production that is 30% greater than that for primary production.
The application of a Bayesian framework for optimizing the HNP performance in a real shale reservoir is introduced for the first time. This work provides practical guidelines for the efficient application of advanced machine learning techniques for optimization under uncertainty, resulting in better decision making.
Since decades, steam-assisted oil recovery processes have been successfully deployed in heavy oil reservoirs to extract bitumen/heavy oil. Current resource allocation practices mostly involve reservoir model-based open loop optimization at the planning stage and its periodic recurrence. However, such decades-old strategies need a complete overhaul as they ignore dynamic changes in reservoir conditions and surface facilities, ultimately rendering heavy oil production economically unsustainable in the low-oil-price environment. Since steam supply costs account for more than 50% of total operating costs, a data-driven strategy that transforms the data available from various sensors into meaningful steam allocation decisions requires further attention.
In this research, we propose a purely data-driven algorithm that maximizes the economic objective function by allocating an optimal amount of steam to different well pads. The method primarily constitutes two components: forecasting and nonlinear optimization. A dynamic model is used to relate different variables in historical field data that were measured at regular time intervals and can be used to compute economic performance indicators (EPI). The variables in the model are cumulative in nature since they can represent the temporal changes in reservoir conditions. Accurate prediction of EPI is ensured by retraining regression model using the latest available data. Then, predicted EPI is optimized using a nonlinear optimization algorithm subject to amplitude and rate saturation constraints on decision variables i.e., amount of steam allocated to each well pad.
Proposed steam allocation strategy is tested on 2 well pads (each containing 10 wells) of an oil sands reservoir located near Fort McMurray in Alberta, Canada. After exploratory analysis of production history, an output error (OE) model is built between logarithmically transformed cumulative steam injection and cumulative oil production for each well pad. Commonly used net-present-value (NPV) is considered as EPI to be maximized. Optimization of the objective function is subject to distinct operating conditions and realistic constraints. By comparing results with field production history, it can be observed that optimum steam injection profiles for both well pads are significantly different than that of a field. In fact, the proposed algorithm provides smooth and consistent steam injection rates, unlike field injection history. Also, the lower steam-oil ratio is achieved for both well pads, ultimately translating into ~19 % higher NPV when compared with field data.
Inspired from state-of-the-art control techniques, proposed steam allocation algorithm provides a generic data-driven framework that can consider any number of well pads, EPIs, and amount of past data. It is computationally inexpensive as no numerical simulations are required. Overall, it can potentially reduce the energy required to extract heavy oil and increase the revenue while inflicting no additional capital cost and reducing greenhouse gas emissions.
The present study provides a comprehensive set of new analytical expressions to help understand and quantify well interference due to competition for flow space between the hydraulic fractures of parent and child wells. Determination of the optimum fracture spacing is a key factor to improve the economic performance of unconventional oil and gas resources developed with multi-well pads. Analytical and numerical model results are combined in our study to identify, analyze, and visualize the streamline patterns near hydraulic fractures, using physical parameters that control the flow process, such as matrix permeability, hydraulic fracture dimensions and assuming infinite fracture conductivity. The algorithms provided can quantify the effect of changes in fracture spacing on the production performance of both parent and child wells. All results are based on benchmarked analytical methods which allow for fast computation, making use of Excel-based spreadsheets and Matlab-coded scripts. Such practical tools can support petroleum engineers in the planning of field development operations. The theory is presented with examples of its practical application using field data from parent and child wells in the Eagle Ford shale (Brazos County, East Texas). Based on our improved understanding of the mechanism and intensity of production interference, the fracture spacing (this study) and inter-well spacing (companion study) of multifractured horizontal laterals can be optimized to effectively stimulate the reservoir volume to increase the overall recovery factor and improve the economic performance of unconventional oil and gas properties.
Lin, Ran (Southwest Petroleum University) | Ren, Lan (Southwest Petroleum University) | Zhao, Jinzhou (Southwest Petroleum University) | Tao, Yongfu (Exploration and Development Research Institute, Yumen Oilfield Company) | Tan, Xiucheng (Southwest Petroleum University) | Zhao, Jiangyu (Southwest Petroleum University)
Multi-stage & multi-cluster fracturing in horizontal well drilling is the core technology in for commercial exploitation of shale gas resevoir. According to vast field data, there is remarkable positive correlation relationship between stimulated reservoir volume (SRV) and shale gas production. Hence, estimating the SRV is essential for both pre-fracturing design and post-fracturing evaluation. However, the forming process of SRV involves with many complex mechanisms, making it is difficult to be simulated.
In this paper, we establish a mathematical model to estimate the SRV by simulating multiple hydraulic fractures propagate, formation stress change and reservoir pressure rise; consequently, the stress and pressure change might make natural fractures occur tensile failure or shear failure, generating a high-conductivity zone (i.e., SRV) in the shale reservoir.
To solve the model, displacement discontinuity method (DDM) is applied to simulate non-planar propagation of multiple hydraulic fractures and calculate formation stress change. Finite difference method (FDM) is used to compute reservoir pressure rise. The natural fractures failure state is determined by tensor formulae derived from Warpinski's failure theory. This SRV estimation method involves a variety of complex but crucial physical mechanisms during shale fracturing process which include unequal flow-rate distribution in different hydraulic fractures, non-planar hydraulic fractures propagation under stress interference, reservoir permeability increases with SRV expanding, two types of natural fracture failure and so on.
A field case study was performed to show the dynamic processes of hydraulic fractures propagation, reservoir permeability increase, and the SRV expansion during shale gas fracturing. Then we compared the simulation results with analytical solution, published papers and on-site microseismic monitoring data to verify our model. Finally, the influence of geological condition and engineering parameters on SRV was investigated by sensitivity analysis.
The reporting of potential resources is essential to assess the future development plan and profitability of a petroleum discovery, but if the project is under appraised and production data are absent, analysts often use analogs for preliminary estimates of technically recoverable volumes. To address this, a workflow is presented for selecting appropriate analogs for unconventional plays and using them to estimate the target play's potential. The proposed technique is demonstrated with a case study of the as-yet undeveloped Bowland Shale, which is the most prominent of the shale plays in the United Kingdom (UK) and is at the early stage of its assessment. The paper describes the current shale gas activity in the UK, highlighting the enviromental constraints placed on would-be Bowland Shale developers, which impact on drilling and production operations and stem from the geographic proximity of urban developments, infrastructure and nature, which limit the size of well pad footprint in the UK where land use is high. Studies have estimated the play's in-place resources for possible future development, but there are few estimates of its corresponding recoverable volumes due to lack of production history. At the outset, a database is created with published minimum-average-maximum ranges of key parameters such as total organic carbon, maturity level, gas filled porosity, permeability, etc. that play a major role in resources estimation and recovery potential for all unconventional plays. A comparison of triangular distributions, key parameter by key parameter, between the target shale play and the analog database, is then carried out using novel graphical and statistical methods to establish a "confidence factor" relating to the analog's viability. The most appropriate analog for the Bowland Shale is chosen from an exhaustive list of North American shale gas plays. Analytical approaches are then used to transform a model of the published type well performance of the selected analog by exchanging key model parameters with those of the target shale play. The paper shows how UK operational constraints can be statistically incorporated into the workflow and have a marked effect on the estimated recovery from the Bowland Shale.
One of the main uncertainties when designing polymer floods is the polymer injectivity, an important parameter that can affect the economics of the process. Reservoir simulation can be used to forecast injectivity, but the process is not straightforward and can be affected by grid size and other factors. Analytical methods are also available for that purpose, but they are considered too simplistic to deal with realistic reservoir conditions. The aim of this paper is to show that this is not the case and that simple analytical tools can be accurate and of great help to predict or history match polymer injectivity.
The analytical method has been developed by Lake in his classical textbook on Enhanced Oil Recovery, but few applications are documented in the literature. This paper will review the method and corresponding equations before presenting several actual field cases of injectivity in polymer flood pilots or tests from several countries that have been matched analytically.
Although it has not been used very often, the method has been found to give very good results in most of the field cases tested in a variety of situations; these cases will be presented along with recommendations on how to apply the method and a discussion of the results. Sensitivities to the various parameters will also be presented. Once the equations are programmed in a spreadsheet, the matching process takes only a few minutes and it is easy to run various scenarios and sensitivities.
Polymer injectivity remains one of the less understood and less predictable aspects of polymer flood projects. This paper will encourage engineers who are planning such projects to use simple yet accurate analytical tools before embarking in more complex and time-consuming reservoir simulations.
Rate transient analysis using log-log plots of rate-normalized pressure (RNP) and its derivative (RNP') versus material balance time have proven helpful in providing estimates of shale matrix permeability and SRV drainage volumes in multiple transverse fracture wells (MTFW's) (
We have constructed an analytical model of MTFW's that accurately predicts individual fracture flow performance for both constant and variable rate and constant bottom hole pressure inner boundary conditions. Using this model, we can accurately compute the pressure disturbance and rate change seen at the whole well and for individual fractures to quantify the degree of interference between fractures for any number of parallel, equally-spaced, and equally-sized fractures. This model has been validated by simulation using a commercial simulator. With both this analytical model and a series of numerical simulations, we investigated the fundamental mechanisms of flow in MTFW's and how the estimation of telf may be improved.
Previous authors have represented the progression of flow regimes in MTFW's as a linear flow period that transitions to a pseudo steady state (or apparently boundary-dominated) flow regime. We show that the same flow response is exhibited by a fully-infinite linear system, calling into question the nature of the "stimulated reservoir volume" (SRV) as a bounded reservoir system. In addition, we show telf can be detected and interpreted as the beginning of the onset of this fracture interference using the "limit of detectability" concept.
Penghui, Su (PetroChina Research Institute of Petroleum Explorationand and Development) | Zhaohui, Xia (PetroChina Research Institute of Petroleum Explorationand and Development) | Ping, Wang (PetroChina Research Institute of Petroleum Explorationand and Development) | Liangchao, Qu (PetroChina Research Institute of Petroleum Explorationand and Development) | xiangwen, Kong (PetroChina Research Institute of Petroleum Explorationand and Development) | Wenguang, Zhao (PetroChina Research Institute of Petroleum Explorationand and Development)
Interest has spread to potential unconventional shale reservoirs in the last decades, and they have become an increasingly important source of hydrocarbon. Importantly, pore structure of shale has considerable effects on the storage, seepage and output of the fluids in shale reservoirs so that reliable fractal characteristics are essential. To better understand the evolution characteristics of pore structure for a shale gas condensate reservoir and their influence on liquid hydrocarbon occurrences and reservoir physical properties, we conducted high-pressure mercury intrusion tests (HPMIs), field emission scanning electron microscopies (FESEM), total organic carbon (TOC), Rock-Eval pyrolysis and saturation measurements on samples from the Duvernay formation. Furthermore, the fractal theory is applied to calculate the fractal dimension of the capillary pressure curves, and three fractal dimensions D1, D2 and D3 are obtained. The relationships among the characteristics of the Duvernay shale (TOC, organic matter maturity, fluid saturation), the pore structure parameters (permeability, porosity, median pore size), and the fractal dimensions were investigated.
The results show that the fractal dimension D1 ranges from 2.44 to 2.85, D2 ranges from 2.09 to 2.15 and D3 ranges from 2.35 to 2.48. D2 and D3 have a good positive correlation. The pore system studied mainly consists of organic pores and microfractures, with the percentage of micropores being 50.38%. TOC has a positive relationship with porosity and D3 due to the development of organic pores. D3 has a positive correlation with gas saturation. With increased D3, median pore size shows a decreasing trend and an increase in permeability and porosity, demonstrating that D3 has a large effect on pore size distribution and the heterogeneity of pore size. In general, D3 has a better correlation with petrophysical and petrochemical parameters. Fractal theory can be applied to better understand the pore evolution, pore size distribution and fluid storage capacity of shale reservoirs.