Overburden characterization with formation pore pressure and anisotropic stress field estimation in the Athabasca Basin, Canada

2019 ◽  
Vol 7 (4) ◽  
pp. T761-T771
Author(s):  
Carmen C. Dumitrescu ◽  
Draga A. Talinga

One of the challenges encountered during the life cycle of an oil-sand thermal-production reservoir is the prediction of the formation pore pressure and in situ stress regime during the assessment phase of the reservoir development and, more importantly, during the development phase. We have investigated the state of formation pore pressure and stress in the overburden — represented by the Clearwater Formation, Grand Rapids Formation, and Colorado Group — of a preproduction oil-sands reservoir situated in the Athabasca Basin of Alberta, Canada. Our methodology integrates pressure data from piezometers, stress data from mini-frac (MF), dipole sonic logs, and elastic properties obtained from multicomponent 3D seismic inversion data. It combines the Terzaghi effective stresses with the Schoenberg and Sayers elastic stiffness matrix for horizontal transversely isotropic fractured materials. The total principal stresses (vertical, minimum, and maximum horizontal stresses) are expressed as functions of the normal fracture weakness (anisotropic correction factor), formation pore pressure, seismic data (Lamé constants), and the Biot-Willis coefficient. The effective principal stresses are estimated from the equivalent total principal stresses and the formation pore pressure multiplied by the Biot-Willis coefficient. On all three overburden intervals analysed, the relations between principal stresses indicate a normal stress regime. The estimated total minimum horizontal stress matches the MF values within 10%. The formation pore pressure, along with the 3D seismically derived estimates of the total and effective principal stresses, allows for better assessment of the caprock integrity and for operational savings based on a reduced number of MF tests. It can also be used for stress estimation within the formations hosting aquifers, which is so important for thermal production. Understanding the subsurface on the reservoir area is important for efficient production, but knowing the subsurface of the overburden is equally important for reducing potential issues due to production.

2019 ◽  
Vol 7 (1) ◽  
pp. T221-T229 ◽  
Author(s):  
Draga A. Talinga ◽  
Carmen C. Dumitrescu

In the Athabasca Basin of Alberta, Canada, the Lower Cretaceous McMurray Formation reservoir is overlain by the Clearwater Formation, a regionally continuous layer composed predominantly of shales with interbedded mudstones. The shales form the reservoir caprock and have the role of blocking the vertical migration of the steam from thermal oil production by confining the stresses and the deformations. We have developed a new method for the 3D seismic reconstruction of the anisotropic stress field, which accounts for the formation pore pressure and the effective stress. We integrated anisotropy estimated from dipole sonic logs with formation pore pressure data from piezometers and elastic properties obtained from multicomponent seismic inversion. The method combines the Terzaghi effective stress with the Schoenberg and Sayers elastic stiffness matrix for horizontal transversely isotropic (HTI) fractured materials. The key points of this method are the estimation of the formation pore pressure in the abnormal regime of the Clearwater Formation, the normal fracture weakness parameter (based on constraints on the compressional velocity of the intact rock under the HTI assumption), and the 3D seismic anisotropic stress field. We expressed the total vertical stress as the weight of the overlying formations, and the total minimum and maximum horizontal stresses as a combination of the total vertical stress, normal fracture weakness, formation pore pressure, Biot-Willis coefficient, and Lamé elastic constants. The effective principal stresses are estimated from the equivalent total principal stresses and the formation pore pressure multiplied by the Biot-Willis coefficient. We observed excellent consistency between the calculated total minimum horizontal stress and mini-frac values. This new method for the 3D seismic reconstruction of the anisotropic stress field allows for the assessment of the caprock integrity and for operational savings based on a reduced number of mini-frac measurements, and it can be used for time-lapse stress estimation within the thermal production reservoir.


2021 ◽  
pp. 1-19
Author(s):  
Aymen Al-Ameri

Summary Sand production is a serious problem in oil and gas wells, and one of the main concerns of production engineers. This problem can damage downhole equipment and surface production facilities. This study presents a sand production case and quantifies sanding risks for an oil field in Iraq. The study applies an integrated workflow of constructing 1D Mechanical Earth Modeling (MEM) and predicting the sand production with multiple criteria such as shear failure during drilling, B index, and critical bottomhole pressure (CBHP) or critical drawdown pressure (CDDP). Wireline log data were used to estimate the mechanical properties of the formations in the field. The predicted sand production propensity was validated based on the sand production history in the field. The interpretation results of some wells anticipated in this study showed that when a shear failure occurs during drilling, the B index is around 2 × 104 MPa or less and the CBHP is equal to the formation pore pressure. For this case, sand control shall be carried out in the initial stage of production. On the other hand, when the shear failure does not exist, the B index is always greater than 2 × 104 MPa, and the CBHP is mostly less than the formation pore pressure. In this case, implementing sand control methods could be postponed as the reservoir pressure undergoes depletion. However, for the anticipated field, sand control is recommended to be carried out in the initial stage of well production even when the CBHP is less than the formation pore pressure since sanding will be inevitable when the reservoir pressure depletes to values close to the initial reservoir pressure. The tentative evaluation of the stress regime showed that a normal fault could be the stress regime for the formations. For a normal fault stress regime, the study explained that when the reservoir permeability is isotropic, an openhole vertical wellbore has less propensity for sand production than a horizontal wellbore. Moreover, when the wellbore azimuth is in the direction of the minimum horizontal stress, the CBHP will be lower than in any other azimuth, and sanding will take place at higher wellbore inclination angles. For the anticipated field, because of the casedhole well completion and the anisotropic reservoir permeability, a horizontal well drilled in the direction of minimum horizontal stress with oriented perforation in the direction of maximum horizontal stress is an alternative method for controlling sand production.


Author(s):  
Mohammad Farsi ◽  
Nima Mohamadian ◽  
Hamzeh Ghorbani ◽  
David A. Wood ◽  
Shadfar Davoodi ◽  
...  

2021 ◽  
Vol 8 ◽  
pp. 55-79
Author(s):  
E. Bakhshi ◽  
A. Shahrabadi ◽  
N. Golsanami ◽  
Sh. Seyedsajadi ◽  
X. Liu ◽  
...  

The more comprehensive information on the reservoir properties will help to better plan drilling and design production. Herein, diagenetic processes and geomechanical properties are notable parameters that determine reservoir quality. Recognizing the geomechanical properties of the reservoir as well as building a mechanical earth model play a strong role in the hydrocarbon reservoir life cycle and are key factors in analyzing wellbore instability, drilling operation optimization, and hydraulic fracturing designing operation. Therefore, the present study focuses on selecting the candidate zone for hydraulic fracturing through a novel approach that simultaneously considers the diagenetic, petrophysical, and geomechanical properties. The diagenetic processes were analyzed to determine the porosity types in the reservoir. After that, based on the laboratory test results for estimating reservoir petrophysical parameters, the zones with suitable reservoir properties were selected. Moreover, based on the reservoir geomechanical parameters and the constructed mechanical earth model, the best zones were selected for hydraulic fracturing operation in one of the Iranian fractured carbonate reservoirs. Finally, a new empirical equation for estimating pore pressure in nine zones of the studied well was developed. This equation provides a more precise estimation of stress profiles and thus leads to more accurate decision-making for candidate zone selection. Based on the results, vuggy porosity was the best porosity type, and zones C2, E2 and G2, having suitable values of porosity, permeability, and water saturation, showed good reservoir properties. Therefore, zone E2 and G2 were chosen as the candidate for hydraulic fracturing simulation based on their E (Young’s modulus) and ν (Poisson’s ratio) values. Based on the mechanical earth model and changes in the acoustic data versus depth, a new equation is introduced for calculating the pore pressure in the studied reservoir. According to the new equation, the dominant stress regime in the whole well, especially in the candidate zones, is SigHmax>SigV>Sighmin, while according to the pore pressure equation presented in the literature, the dominant stress regime in the studied well turns out to be SigHmax>Sighmin>SigV.  


2019 ◽  
Vol 125 ◽  
pp. 15001
Author(s):  
Benny Abraham Bungasalu ◽  
M. Syamsu Rosid ◽  
Don S. Basuki

The subsurface pressure analysis is used to detect the overpressure and problems in the well that will be drilled based on exploration well data. Various problems were found while drilling operations carried out on A and B wells, namely, Kick and Pipe sticking which cause a high Non-Productive Time (NPT). This research is conducted to identify the mechanism of overpressure formation in Tight Sand Gas and Shale Gas in the Jambi Sub-Basin. Furthermore, to predict pore pressure using the Drilling Efficiency and Mechanical Specific Energy (DEMSE) and Bowers method. The final result will be a 3D pore pressure cube in the area based on quantitative analysis of post-stack seismic inversion. The results of the pore pressure analysis from the wells and the 3D pore pressure model indicate that top of overpressure occurs in the Gumai Formation, then it is decreasing gradually approaching the hydrostatic pressure on the Basement. The mechanisms of overpressure are caused by under compaction, fluid expansion (kerogen maturation). The Gumai Formation and Talang Akar Formation are shale rocks so the type of mud weight that is well used is oil based mud (OBM).


SPE Journal ◽  
2016 ◽  
Vol 21 (05) ◽  
pp. 1721-1742 ◽  
Author(s):  
Mazda Irani ◽  
Ian Gates

Summary Li et al. (2004) described three zones at the edge of steam chambers on the basis of drainage conditions: drained, partially drained, and undrained. In the drained zone, the pore pressure is controlled by injection pressure, and fluid mobility within this region is sufficient to drain additional pore pressures because of shear dilation and pore-fluid thermal expansion. The undrained zone lies beyond the partially drained zone and extends to virgin reservoir far beyond the chamber. In this zone shearing behaves under undrained conditions; by this, Li et al. (2004) mean no volume change occurs but shear lead to changes in pore pressure. Li et al. (2004) proposed that the boundaries of these zones are dependent on bitumen viscosity, which relates to the temperature distribution beyond the steam interface. Because drained/undrained conditions affect the geomechanics at the edge of the chamber, we investigate whether the assumption of Li et al. (2004) that there is no volume change within the sheared zone is correct and is supported by field data. Here, we establish the physics associated with the undrained zone at the edge of steam-assisted gravity-drainage steam chamber and explore the pressure front vs. temperature front of different oil-sand field projects. The results reveal that the drained zone governed by pressure-front advancement is greater in extent than the sheared zone. The thermodynamics of the undrained zone are discussed to derive a new theory for mechanothermal phenomena at the edge of the chamber. The results from the theory show that the drained zone extends beyond the temperature front and thus, from a geomechanical point of view, the system solely consists of the drained and partially drained zones.


2006 ◽  
Vol 46 (1) ◽  
pp. 307
Author(s):  
B.A. Camac ◽  
S.P. Hunt ◽  
P.J. Boult ◽  
M. Dillon

In distinct element (DEM) numerical stress modelling, the principal stress magnitudes and orientations are applied to the boundary of the 3D model. Due to data restrictions and typical depths of investigation, it is possible to have much uncertainty in the conventional methodologies used to constrain the regional principal stress magnitudes and orientations.A case study from the Kupe field in the Taranaki Basin, New Zealand is presented where the uncertainty in the input data made it difficult to determine which stress regime—a transitional normal/strike-slip or reverse/thrust—is active at reservoir depth (approximately 3,000 m). The magnitudes and orientation of the principal stresses were constrained using published techniques. A sensitivity analysis was applied to account for the uncertainty in the input data. A model of the Kupe field incorporating 18 major faults was subsequently loaded under both derived stressed regimes, using the calculated magnitudes.Borehole breakout analysis was used to acquire interpreted orientations of the maximum principal stress (Shmax). The work presented herein describes a different or unconventional approach to the general petroleum geomechanics methodology. Typically, the breakout data is averaged to get one data point per well location. Here, all breakout data is retained and displayed vertically. The data is actively used and the variations with depth can be seen to show how faults can generate local perturbations of the regional stress trajectory. These data are then used to compare the observed or field indications of the breakouts along the borehole with the modelled Shmax predicted by both end point DEM stress models. This comparison has provided additional confidence in the derived stress regime and the derived stress models for the Kupe field. The stress models are used to predict areas of enhanced hydrocarbon pooling and low seal integrity.


2016 ◽  
Vol 61 (1) ◽  
pp. 69-92
Author(s):  
Yuliy A. Dashevsky ◽  
Semen Petrov ◽  
Alexandr N. Vasilevskiy ◽  
Oleg B. Bocharov ◽  
Gleb V. Dyatlov

2020 ◽  
Author(s):  
Rebecca O. Salvage ◽  
David W. Eaton

<p>On 30 November 2018, three felt earthquakes occurred in quick succession close to the city of Fort St. John, British Columbia, likely as a direct response to a hydraulic fracturing operation in the area. Events appear tightly clustered spatially within the upper 10 km of the crust. Hypocenters locate at the confluence between a large scale reverse faulting regime (in the north-west, probably due to the influence of the Rocky Mountain fold and thrust belt) and an oblique strike slip faulting regime (in the south-east, probably due to the influence of the Fort St. John Graben), resulting in a variety of focal mechanisms and a very complex local stress regime. Further analysis of the principal stresses suggests that σ<sub>1</sub> is well constrained and close to horizontal, whereas σ<sub>2</sub> and σ<sub>3</sub> are poorly constrained, and can alternate between the horizontal and the vertical plane. Here, we present an overview of the temporal and spatial evolution of this seismic sequence and its relationship to hydraulic fracturing operations in the area, and examine the influence of large-scale regional tectonic structures on the generation of seismicity on this occasion.</p><p> </p>


Sign in / Sign up

Export Citation Format

Share Document