Estimation of Seismic Wave Attenuation from 3D Seismic Data: A Case Study of OBC Data Acquired in an Offshore Oilfield

Previous studies performed in Abu Dhabi oilfields, United Arab Emirates, revealed the direct link of seismic wave attenuation to petrophysical properties of rocks. However, all those studies were based on zero offset VSP data, which limits the attenuation estimation at one location only. This is due to the difficulty of estimating attenuation from 3D seismic data, especially in carbonate rocks. To overcome this difficulty, we developed a workflow based on the centroid frequency shift method and Gabor transform which is optimized by using VSP data. The workflow was applied on 3D Ocean Bottom Cable seismic data. Distinct attenuation anomalies were observed in highly heterogeneous and saturated zones, such as the reservoirs and aquifers. Scattering shows significant contribution in attenuation anomalies, which is unusual in sandstones. This is due to the complex texture and heterogeneous nature of carbonate rocks. Furthermore, attenuation mechanisms such as frictional relative movement between fluids and solid grains, are most likely other important causes of attenuation anomalies. The slight lateral variation of attenuation reflects the lateral homogeneous stratigraphy of the oilfield. The results demonstrate the potential of seismic wave attenuation for delineating heterogeneous zones with high fluid content, which can substantially help for enhancing oil recovery.

Single-step probabilistic inversion of 3D seismic data of a carbonate reservoir in Southwest Iran

We use a sampling-based Markov chain Monte Carlo method to invert seismic data directly for porosity and quantify its uncertainty distribution in a hard-rock carbonate reservoir in Southwest Iran. The noise that remains on seismic data after the processing flow is correlated with the bandwidth in the range of the seismic wavelet. Hence, to account for the inherent correlated nature of the band-limited seismic noise in the probabilistic inversion of real seismic data, we assume the estimated seismic wavelet as a suitable proxy for capturing the coupling of noise samples. In contrast to the common approach of inserting a delta function on the main diagonal of the covariance matrix, we insert the seismic wavelet on its main diagonal. We also calibrate an empirical and a semi-empirical inclusion-based rock-physics model to characterize the rock-frame elastic moduli via a lithology constrained fitting of the parameters of these models, i.e. the critical porosity and the pore aspect ratio. These calibrated rock-physics models are embedded in the inversion procedure to link petrophysical and elastic properties. In addition, we obtain the pointwise critical porosity and pore aspect ratio, which can potentially facilitate the interpretation of the reservoir for further studies. The inversion results are evaluated by comparing with porosity logs and an existing geological model (porosity model) constructed through a geostatistical simulation approach. We assess the consistency of the geological model through a geomodel-to-seismic modeling approach. The results confirm the performance of the probabilistic inversion in resolving some thin layers and reconstructing the observed seismic data. We present the applicability of the proposed sampling-based approach to invert 3D seismic data for estimating the porosity distribution and its associated uncertainty for four subzones of the reservoir. The porosity time maps and the facies probabilities obtained via porosity cut-offs indicate the relative quality of the reservoir’s subzones.

Use of Attributes on 3D Seismic Data for Studying Mishrif Formation in East Abu-Amoud Field, South-Eastern Iraq

The seismic method depends on the nature of the reflected waves from the interfaces between layers, which in turn depends on the density and velocity of the layer, and this is called acoustic impedance. The seismic sections of the East Abu-Amoud field that is located in Missan Province, south-eastern Iraq, were studied and interpreted for updating the structural picture of the major Mishrif Formation for the reservoir in the field. The Mishrif Formation is rich in petroleum in this area, with an area covering about 820 km2. The horizon was calibrated and defined on the seismic section with well logs data (well tops, check shot, sonic logs, and density logs) in the interpretation process to identify the upper and lower boundaries of the Formation. Seismic attributes were used to study the formation, including instantaneous phase attributes and relative acoustic impedance on time slice of 3D seismic data . Also, relative acoustic impedance was utilized to study the top of the Mishrif Formation. Based on these seismic attributes, karst features of the formation were identified. In addition, the nature of the lithology in the study area and the change in porosity were determined through the relative acoustic impedance The overlap of the top of the Mishrif Formation with the bottom of the Khasib Formation was determined because the Mishrif Formation is considered as an unconformity surface.

Controls of Deep-Seated Faults and Folds on Hydrocarbon Fluid Migration and Accumulation in Sedimentary Basins: A Case Study from the Northwestern Sichuan Basin, China

Deep-seated faults and folds of foreland basin systems have become important exploration targets in the recent years because they are crucial in controlling fluid migration and hydrocarbon accumulation. In this study, we analyzed the characteristics and formation history of these structures in the northwestern Sichuan Basin using recently acquired two-dimensional (2D) and three-dimensional (3D) seismic data. The seismic interpretation revealed that the thrust sheets, tectonic wedges, and foredeep were well developed in the northwestern Sichuan Basin from the mountain to the basin. Forward thrusts, fault-bend folds, and wedges are the main types of structures in the thrust sheets and tectonic wedges. The deep-seated faults and folds were easily recognized in the high-resolution 3D seismic data. The imbricate thrust faults that merged into detachment layers of the Lower Cambrian are the main types of structures in the foredeep, and they show a prominent strike-slip influence in the horizontal direction. The formation of these structures in the foredeep in the northwestern Sichuan Basin mainly endured two stages of thrusting, including those during the Middle-to-Late Triassic and Cenozoic. Based on the tectonic evolution and seismic data, we infer that these deep-seated faults and folds in the foredeep may have formed earlier than the northern Longmen Shan fold-and-thrust belts and they may have been initially active in the late of Early Triassic and reactive during the Cenozoic. Furthermore, evaporites in the Lower and Middle Triassic were crucial in forming these structures. The petroleum exploration data suggested that the deep-seated faults can facilitate hydrocarbon accumulation. The thrust faults in the foredeep were more likely to act as migration pathways for fluids instead of sealing barriers along the horizontal direction. The interconnected reservoirs of deep-seated folds possess a great potential to allow large-scale hydrocarbon accumulation. Our study provides a good example for evaluating the hydrocarbon exploration potential in the deeply buried area in the sedimentary basin.

A Machine-Learning Based Workflow for Predicting Overpressure in a Stiff Dolomitic Formation

Overpressure refers to the abnormally high subsurface pressure that may exceed hydrostatic pressure at a given depth. Its characterization is an important part of subsurface characterization as it allows to complete drilling operations in a safe and optimal way. In dolomitic formations, however, the prediction of such overpressure is especially challenging because of (1) the high degree of lateral variability of the formations, (2) the limited effect of overpressure on tight rocks elastic parameters, and (3) the complexity of physical processes involved to form overpressure. In addition to these factors, existing experimental models generally used to relate elastic parameters to pressure are often not well calibrated to carbonate rocks. The alternative to existing purely physical approaches is a data-driven model that leverages data from offset wells. We show that due to the complexity of the characterization question to be solved, an end-to-end machine learning based approach is deemed to fail. Instead of a fully automated approach, we show a semi-supervised workflow that integrates seismic, geological data, and overpressure observations from previously drilled wells to map overpressure regions. Attribute maps are first extracted from a 3D seismic data set in an overpressured geological formation of interest. An auto-encoder is then used to learn a more compact representation of data, resulting in a reduced number of latent attributes. Then, a hand-tailored semi-supervised approach is applied, which is a combination of clustering method (here based on DBSCAN algorithm) and Bayesian classification to determine overpressure risk degree (no risk, mild, or high risk). The approach described in this study is compared to direct end-to-end models and significantly outperforms them with an error on a blind well prediction of around 25%. The overpressure probability maps resulting from the models can be used later for the optimization of drilling processes and to reduce drilling hazards.

Using AI/ML to Explore & Develop Quickly and Efficiently

See how application of a fully trained Artificial Intelligence (AI) / Machine Learning (ML) technology applied to 3D seismic data volumes delivers an unbiased data driven assessment of entire volumes or corporate seismic data libraries quickly. Whether the analysis is undertaken using onsite hardware or a cloud based mega cluster, this automated approach provides unparalleled insights for the interpretation and prospectivity analysis of any dataset. The Artificial Intelligence (AI) / Machine Learning (ML) technology uses unsupervised genetics algorithms to create families of waveforms, called GeoPopulations, that are used to derive Amplitude, Structure (time or depth depending on the input 3D seismic volume) and the new seismic Fitness attribute. We will show how Fitness is used to interpret paleo geomorphology and facies maps for every peak, trough and zero crossing of the 3D seismic volume. Using the Structure, Amplitude and Fitness attribute maps created for every peak, trough and zero crossing the Exploration and Production (E&P) team can evaluate and mitigate Geological and Geophysical (G&G) risks and uncertainty associated with their petroleum systems quickly using the entire 3D seismic data volume.

Dating submarine landslides using the transient response of gas hydrate stability

Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of slope failure. Herein, we present a novel technique for constraining the age of submarine landslides without sediment core dating. Underneath a submarine landslide in the Orca Basin, Gulf of Mexico, in 3D seismic data we map an irregular bottom simulating reflection (BSR), which mimics the geometry of the pre-slide seafloor rather than the modern bathymetry. Based on the observed BSR, we suggest that the gas hydrate stability zone (GHSZ) is currently adjusting to the post-slide sediment temperature perturbations. We apply transient conductive heat flow modeling to constrain the response of the GHSZ to the slope failure, which yields a most likely age of ~8 ka demonstrating that gas hydrate systems can respond to slope failures even on the millennia timescales. We also provide an analytical approach to rapidly determine the age of submarine slides at any location.

Identification of submarine landslides in the Colombian Caribbean Margin (Southern Sinú Fold Belt) using seismic investigations

Submarine landslides can be triggered by several processes and involve a variety of mechanisms. These phenomena are important sediment transport processes, but they also constitute a significant geohazard. Mapping of the southwestern Caribbean Sea using 3D seismic data has allowed identification of several submarine landslides in the Colombian Margin in the area dominated by the Southern Sinú Fold Belt (SSFB). A poststack depth-migrated seismic cube survey with a 12.5 by 12.5 m bin spacing was used to identify landslides in an area covering 5746 km2. Landslides were interpreted using a seafloor morphologic parameter identification process and the internal deformation of the slope-forming material, as seen from seismic data. A total of 93 landslides were identified and classified based on their movement styles as follows: 52 rotational, 29 translational, and 12 complex landslides. In addition, 12 distinct deformational zones and a zone of mass transport complex (MTC) were identified. Five different ground condition terrains were interpreted based on landslide type and distribution as well as in geologic structures and seismic reflection analysis. Two main processes seem to influence landslides in the study area. First is the folding and faulting involved in the SSFB evolution. This process results in oversteepened slopes that start as deformational zones and then fail as translational or rotational slides. Those individual landslides progressively become complex landslide zones that follow geologic structural orientation. Second is the continental shelf break erosion by debris flows, which fills in intraslope subbasins and continental rise with several MTCs. According to the results, risk of damage by landslides increases in distances shorter than 4 km along structural ridge foothills in the study zone.

Seismic Reservoir Characterisation of Tyumen Formation in Frolov Megadepression

The article discusses the approaches and methods to study the Middle Jurassic deposits of the Tyumen Formation within the Frolov megadepression (West Siberian oil and gas province), which have high hydrocarbon potential. The materials refer to several areas with available 3D seismic data and several dozen oil wells. The problems of seismic interpretation and its application for geological modeling are considered. We also propose several ways to overcome them.