Seismic Attribute Analyses and Attenuation Applications for Detecting Gas Hydrate Presence

Identifying gas hydrates in the oceanic subsurface using seismic reflection data supported by the presence of a bottom simulating reflector (BSR) is not an easy task, given the wide range of geophysical methods that have been applied to do so. Though the presence of the BSR is attributed to the attenuation response, as seismic waves transition from hydrate-filled sediment within the gas hydrate stability zone (GHSZ) to free gas-bearing sediment below, few studies have applied a direct attenuation measurement. To improve the detection of gas hydrates and associated features, including the BSR and free gas accumulation beneath the gas hydrates, we apply a recently developed method known as Sparse-Spike Decomposition (SSD) that directly measures attenuation from estimating the quality factor (Q) parameter. In addition to performing attribute analyses using frequency attributes and a spectral decomposition method to improve BSR imaging, using a comprehensive analysis of the three methods, we make several key observations. These include the following: (1) low-frequency shadow zones seem to correlate with large values of attenuation; (2) there is a strong relationship between the amplitude strength of the BSR and the increase of the attenuation response; (3) the resulting interpretation of migration pathways of the free gas using the direct attenuation measurement method; and (4) for the data analyzed, the gas hydrates themselves do not give rise to either impedance or attenuation anomalies that fully differentiate them from nearby non-hydrate zones. From this last observation, we find that, although the SSD method may not directly detect in situ gas hydrates, the same gas hydrates often form an effective seal trapping and deeper free gas accumulation, which can exhibit a large attenuation response, allowing us to infer the likely presence of the overlying hydrates themselves.

3D GPR attributes to generate classified facies models: A case study from a dune island

Ground-penetrating radar (GPR) is a standard geophysical technique to image near-surface structures in sedimentary environments. In such environments, GPR data acquisition and processing are increasingly following 3D strategies. However, the processed GPR data volumes are typically still interpreted using selected 2D slices and manual concepts such as GPR facies analyses. In seismic volume interpretation, the application of (semi-)automated and reproducible approaches such as 3D attribute analyses as well as the production of attribute-based facies models are common practice today. In contrast, the field of 3D GPR attribute analyses and corresponding facies models is largely untapped. We develop and apply a workflow to produce 3D attribute-based GPR facies models comprising the dominant sedimentary reflection patterns in a GPR volume which images complex sandy structures on the dune island of Spiekeroog (Northern Germany). After presenting our field site and details regarding our data acquisition and processing, we calculate and filter 3D texture attributes to generate a database comprising the dominant texture features of our GPR data. Then, we perform a dimensionality reduction of this database to obtain meta texture attributes, which we analyze and integrate using composite imaging and (also considering additional geometric information) fuzzy c-means cluster analysis resulting in a classified GPR facies model. Considering our facies model and a corresponding GPR facies chart, we interpret our GPR data set in terms of near-surface sedimentary units, the corresponding depositional environments, and the recent formation history at our field site. Thus, we demonstrate the potential of the proposed workflow, which represents a novel and clear strategy to perform a more objective and consistent interpretation of 3D GPR data collected across different sedimentary environments.

Application of seismic attribute analysis in fluvial seismic geomorphology

Seismic attributes can be important predictors, either qualitative or quantitative, of reservoir geometries when they are correctly used in reservoir characterization studies. This paper discusses seismic attribute analyses and their usefulness in seismic geomorphology study of Moragot field of Pattani Basin, Gulf of Thailand. Early to Middle Miocene fluvial channel and overbank sands are the reservoirs in Pattani Basin. Due to their limited horizontal and vertical distribution, it is not always possible to predict the geometry and distribution of these sands based on the conventional seismic interpretation. This study utilized various seismic attributes, e.g., RMS amplitude analysis, spectral decomposition, semblance and dip-steered similarity, RGB blending to image the geometry and the spatial distribution of sand bodies in horizon and stratal slices at different stratigraphic intervals. Attribute analyses reveal, at shallow stratigraphic levels, RMS and semblance can successfully identify channel-shaped sand bodies and mud-filled channels associated with channel belts. On the other hand in deeper stratigraphic intervals, sand distribution can be imaged more effectively by using spectral decomposition and dip-steered similarity volumes. High-frequency spectral decomposition slices can image thin sands, and low-frequency slices can image thick sands quite effectively in deeper intervals. RGB blending of different frequency slices is particularly useful in delineating channel systems of various dimensions at deeper intervals. These images show the distribution of sands and mud-filled channels at various stratigraphic levels. The width of channel belts varies from 200 m to 3 km. These channel belts are N–S or NW–SE oriented. From the channel pattern and their dimensions, depositional environments can be predicted. Mud-filled channels identified in the horizon slices will act as a connectivity barrier between sand bodies at either side of the channel. They can also act as lateral and up-dip seal to form stratigraphic traps. The seismic attribute analyses clearly show the geometry and spatial distribution of sand bodies. Hence, this method for predicting sand body geometry might help in field development planning as well as in reducing exploration risk.

Investigating fault continuity associated with geologic carbon storage planning in the Illinois Basin

Because the confinement of [Formula: see text] in a storage reservoir depends on a stratigraphically continuous set of seals to isolate the fluid in the reservoir, the detection of structural anomalies is critical for guiding any assessment of a potential subsurface carbon storage site. Employing a suite of 3D seismic attribute analyses (as opposed to relying upon a single attribute) maximizes the chances of identifying geologic anomalies or discontinuities (e.g., faults) that may affect the integrity of a seal that will confine the stored [Formula: see text] in the reservoir. The Illinois Basin, a major area for potential carbon storage, presents challenges for target assessment because geologic anomalies can be ambiguous and easily misinterpreted when using 2D seismic reflection data, or even 3D data, if only conventional display techniques are used. We procured a small 3D seismic reflection data set in the central part of the basin (Stewardson oil field) to experiment with different strategies for enhancing the appearance of discontinuities by integrating 3D seismic attribute analyses with conventional visualizations. Focusing on zones above and below the target interval of the Cambrian Mt. Simon Sandstone, we computed attribute traveltime slices (combined with vertical views) based on discontinuity computations, crossline-directed amplitude change, azimuth of the dip, shaded relief, and fault likelihood attributes. The results provided instructive examples of how discontinuities (e.g., subseismic scale faults) may be almost “invisible” on conventional displays but become detectable and mappable using an appropriate integration of 3D attributes. Strong discontinuities in underlying Precambrian basement rocks do not necessarily propagate upward into the target carbon storage interval. The origin of these discontinuities is uncertain, but we explored a possible strike-slip role that also explains the localization of a structural embayment developed in Lower Paleozoic strata above the basement discontinuities.