Tight reservoirs around the world contain a significant volume of hydrocarbons; however, the heterogeneity of these reservoirs limits the recovery of the original oil in place to less than 20%. Accurate characterization is therefore needed to understand variations in reservoir properties and their effects on production. Water saturation (Sw) has always been challenging to estimate in ultra-tight reservoirs such as the Bakken Formation due to the inaccuracy of resistivity-based methods. While machine learning (ML) has proven to be a powerful tool for predicting rock properties in many tight formations, few studies have been conducted in reservoirs of similar complexity to the Bakken Formation, which is an ultra-tight, multimineral, low-resistivity reservoir. This study presents a workflow for Sw prediction using well logs, core data, and ML algorithms. Logs and core data were gathered from 29 wells drilled in the Bakken Formation. Due to the inaccuracy and lack of robustness of the tried and tested regression models (e.g., linear regression, random forest regression) in predicting Sw as a continuous variable, the problem was reformulated as a classification task. Instead of exact values, the Sw predictions were made in intervals of 10% increments representing 10 classes from 0% to 100%. Gradient boosting and random forest classifiers scored the best classification accuracy, and these two models were used to construct a voting classifier that achieved the best accuracy of 85.53%. The ML model achieved much better accuracy than conventional resistivity-based methods. By conducting this study, we aim to develop a new workflow to improve the prediction of Sw in reservoirs where conventional methods have poor performance.

Multilateral well drilling technology has recently assisted the drilling industry in improving borehole contact area and reducing operation time, while maintaining a competitive cost. The most advanced multilateral well drilling method is Fishbone drilling (FbD). This method has been utilized in several hydrocarbon fields worldwide, resulting in high recovery enhancement and reduced carbon emissions from drilling. FbD involves drilling several branches from laterals and can be considered as an alternative method to hydraulic fracturing to increase the stimulated reservoir volume. However, the expected productivity of applying a Fishbone well from one field to another can vary due to various challenges such as Fishbone well design, reservoir lithology, and accessibility. Another challenge is the lack of existing analytical models and the effect of each Fishbone parameter on the cumulative production, as well as the interaction between them. In this paper, analytical and empirical productivity models were modified for FbD in a dry gas reservoir. The modified analytical model showed a higher accuracy with respect to the existing model. It was also compared with the modified empirical model, which proved its higher accuracy. Finally, machine learning algorithms were developed to predict FbD productivity, which showed close results with both analytical and empirical models.

Salt caverns serve as underground storage for crude oil, natural gas, compressed air, carbon dioxide, and hydrogen. Key stages of cavern development for storage purposes include design, construction, storage, and abandonment. The design phase addresses optimal cavern shape, size, pillar dimensions, number of caverns, the impact of interbeds, and cyclic loading while considering the creep behavior of salt and the mechanical behavior of surrounding layers. During this phase, geological factors such as depth, thickness, and the quality of salt are considered. For construction, two main methods—direct leaching and reverse leaching—are chosen based on design specifications. The storage stage includes the injection and withdrawal of gases in a cyclic manner with specific injection rates and pressures. After 30 to 50 years, the caverns are plugged and abandoned. The geological limitation of salt domes makes it essential to look for more bedded evaporites. This study provides a comprehensive review of bedded evaporites, including their origin and depositional environment. The stability of caverns in all these stages heavily relies on geomechanical analysis. Factors affecting the geomechanics of bedded salts such as mineralogy, physical properties, and mechanical properties are reviewed. A list of bedded evaporites in the U.S. and Canada, including their depth, thickness, and existing caverns, is provided. Additionally, this study discusses the main geomechanical considerations influencing design, solution mining, cyclic loading, and abandonment of caverns in bedded salt caverns.

Accurate prediction of the fracture geometry before the operation of a hydraulic fracture (HF) job is important for the treatment design. Simplified planar fracture models, which may be applicable to predict the fracture geometry in homogeneous and continuous formations, fail in case of fractured reservoirs and laminated formations such as shales. To gain a better understanding of the fracture propagation mechanism in laminated formations and their vertical geometry to be specific, a series of numerical models were run using XSite, a lattice-based simulator. The results were studied to understand the impact of the mechanical properties of caprock and injection parameters on HF propagation. The tensile and shear stimulated areas were used to determine the ability of HF to propagate vertically and horizontally. The results indicated that larger caprock Young’s modulus increases the stimulated area (SA) in both vertical and horizontal directions, whereas it reduces the fracture aperture. Also, larger vertical stress anisotropy and tensile strength of caprock and natural interfaces inhibit the horizontal fracture propagation with an inconsiderable effect in vertical propagation, which collectively reduces the total SA. It was also observed that an increased fluid injection rate suppresses vertical fracture propagation with an insignificant effect on horizontal propagation. The dimensionless parameters defined in this study were used to characterize the transition of HF propagation behavior between horizontal and vertical HFs.

Correct estimation of hydraulic fracture height growth is a critical step in the design of Hydraulic Fracturing (HF) treatment, as it maximizes the reservoir stimulation and returns on investment. The height of the fractures is governed by several in situ conditions, especially stress variation with depth. The common workflow to estimate stress is by building the mechanical earth model (MEM) and calibrating it using the Diagnostic Fracture Injection Test (DFIT). However, DFIT interpretation is a complex task, and depending on the method used, different results may be obtained that will consequently affect the predicted hydraulic fracture height. This work used the tangent and compliance methods for DFIT interpretation, along with isotropic and anisotropic stress profiles, to estimate the HF height growth using numerical modeling in a 3D planar HF simulator. Data from two wells in the Anadarko Basin were used in this study. The predicted height was compared with microseismic data. The results showed that even though the tangent method fits better to the isotropic stress profile, HF did not match with the microseismic data. On the contrary, the anisotropic stress profile showed a good match between the compliance DFIT model and the microseismic events. Based on the discussions presented in this study, the validity of the DFIT interpretation is debatable, and when the formations are anisotropic, the isotropic model fails to correctly estimate the minimum stress profile, which is the main input for the estimation of the fracture height. This is in addition to the fact that some researchers have questioned the use of the tangent method in low-permeability formations.

Multi-stage hydraulic fracturing (MsHF) is the main technology to improve hydrocarbon recovery from shale plays. Associated with their rich organic contents and laminated depositional environments, shales exhibit transverse isotropic (TI) characteristics. In several cases, the lamination planes are horizontal in shale formations with a symmetric axis that are vertical to the bedding plane; hence, shale formations are known as transverse isotropic vertical (TIV) rocks. Ignoring the TIV nature of shale formations leads to erroneous estimates of in situ stresses and consequently to inefficient designs of fractural geometry, which negatively affects the ultimate recovery. The goal of this study is to investigate the effects of TIV medium characteristics on fractural geometry, spacing, and stress shadow development in the Jurassic Tuwaiq Mountain formation (TMF) in the Jafurah basin, which is a potential unconventional world-class play. This formation is the main source for prolific Jurassic oil reservoirs in Saudi Arabia. On the basis of a petrophysical evaluation in the Jafurah basin, TMF exhibited exceptional unconventional gas characteristics, such as high total organic content (TOC) and low clay content, and it was in the proper maturity window for oil and gas generation. The unconventional Jafurah field covers a large area that is comparable to the size of the Eagle Ford shale play in South Texas, and it is planned for development through multi-stage hydraulic fracturing technology. In this study, analytical modeling was performed to estimate the fractural geometry and in situ stresses in the anisotropic medium. The results show that the Young’s modulus anisotropy had a noticeable impact on fractural width, whereas the impact of Poisson’s ratio was minimal. Moreover, we investigated the impact of stress anisotropy and other rock properties on the stress shadow, and found that a large stress anisotropy could result in fractures being positioned close to one another or theoretically without minimal fractural spacing concerns. Additionally, we estimated the fractural aspect ratio in different propagation regimes and observed that the highest aspect ratio had occurred in the fractural toughness-dominated regime. This study also compares the elastic properties and confirms that TMF exhibited greater anisotropic properties than those of Eagle Ford. These findings have practical implications for field operations, particularly with regard to the fractural geometry and proppant placement.