Seismic Amplitude Research Articles

Volcanic-controlled basin architecture variation and dynamic sediment filling in the South Lufeng Sag, South China Sea

Volcanic-controlled changes of basin architectures and sediments have given rise to diverse and rapidly evolving “dynamic geomorphology,” which poses challenges to successful petroleum development. Nonetheless, the accompanying volcanic deposits and landforms often produce rich and distinct seismic responses within conventional sedimentary layers, offering vital support in addressing this critical issue. Based on this, we compared the volcanic activity stages, evolution of geomorphology, and sediment dispersal processes to determine the process of basinal evolution by analyzing the episodes and distribution of volcanic-related deposits, sub-stages, seismic facies, and attributes of deposits of the South Lufeng Sag in the Pearl River Mouth Basin, South China Sea. The widespread development of strong amplitude attributes and chaotic reflection variance cube attribute spatial perspectives reveal the evolution of both intrusive and effusive facies in the lower Wenchang Formation across three stages of volcanic activity. The seismic amplitude difference of the andesite and andesitic tuff unveils the transition of the depositional pathway from north to south in the west depression with change of depocenters, lays the foundation for the dynamic geomorphic evolution of the “transitional basin” in the region. In the middle Huizhou-Lufeng uplift, the chaotic reflection with weak amplitude and obvious diapir structures revealed by andesite constrain the distribution of sedimentary space and the transformation of sedimentary pathways within the basin, and the evolution process of the “restricted basin” is summarised in terms of sedimentary filling processes. The seismic amplitude differences of basalt and andesitic tuff - related sediments in the east depression record the process of transformation from a single sedimentary channel to multi-pathway transport, which is summarised as the evolutionary model of the “diffusive basin”. This study offers a meticulous exploration into the evolutionary trajectory of dynamic geomorphology and sedimentary accumulation in a basin regulated by volcanic activity. The findings pave the way for the extensive application of this methodology in petroleum development across various geographic regions.

Enhancement of CO2 monitoring in the sleipner field (north sea) using seismic inversion based on simulated annealing of time-lapse seismic data

The primary aim of this research is to enhance seismic data interpretation and CO2 monitoring by utilizing seismic inversion techniques based on the simulated annealing method. Simulated annealing is a global optimization technique employed for inverting seismic data and provides better results as compared with local optimization-based inversion. This methodology is implemented in the Utsira Formation, located at a depth of 1000 m within the Sleipner Field, Norway. The study encompasses the analysis of three sets of time-lapse seismic data, first from 1994 (pre-injection), followed by surveys in 1999 and 2001, corresponding to the injection of 2.35 million tonnes and 4.26 million tonnes of CO2, respectively. Firstly, synthetic data is used to check the reliability of the algorithm followed by real data application. This process starts by performing the inversion analysis on the synthetic data which shows a decrease in the impedance values observed at the injection site whereas the seismic amplitude increases. The qualitative as well as quantitative analysis depicts that the algorithm works satisfactorily. The same process is applied to the real data from the Sleipner field. Acoustic impedances are calculated using a simulated annealing-based inversion scheme for the pre-injection case in 1994 and post-injection scenarios in 1999 and 2001. Because of the presence of injected CO2 in the years 1999 and 2001, a low impedance zone that ranged from 2000 m/s*g/cc to 2400 m/s*g/cc appeared at the time interval of 0.85–1.10sec. The interpretation of the inverted impedance section and seismic attribute analysis show no signature of CO2 leakage. The results indicated that the inverted section which is derived from the SA optimization technique shows very clear CO2 information offering a more realistic representation with enhanced resolution of the CO2 plume and its migratory paths.

Application of multi-component seismic data in identifying dolomite reservoirs in the Sichuan Basin

Abstract A comprehensive drilling of wells has been conducted in the Permian Qixia Formation in the central Sichuan Basin, revealing a significant number of dolomite reservoirs. High- and medium-porosity dolomite reservoirs are the main gas-producing reservoirs in the Qixia Formation. Seismic PP-wave data show a “bright spot” for high-porosity dolomite reservoir formations but weak responses for medium-porosity dolomite reservoir formations, which is attributed to the inability of P waves to distinguish between medium-porosity reservoirs and limestone. However, medium-porosity dolomite and limestone have different S-wave velocities. Therefore, in this study, the identification of different-porosity dolomite reservoirs using multi-component seismic data was investigated. A comprehensive analysis of the elastic waves by forward modeling shows that the PS-wave amplitude is more sensitive to medium-porosity dolomite than the PP-wave amplitude. Therefore, medium-porosity dolomite reservoirs can be predicted using the amplitude attributes of the PS-wave, and high-porosity dolomite reservoirs can be characterized using the PP-wave. Meanwhile, the elastic parameter λρ (the product of Lame constant λ and density ρ), which is highly correlated with the dolomite content, can be used as an indicator of dolomite formations. Furthermore, compared to the results of PP-wave inversion, the elastic parameters derived from the joint inversion of PP- and PS-waves exhibited a better correspondence with the well-logging results. The comprehensive use of the seismic amplitude responses of PP- and PS-waves and multi-component seismic joint inversion can effectively predict high- and medium-porosity dolomite reservoirs. The predicted results can support the exploration and development of the Qixia Formation.

Examining 22 Years of Ambient Seismic Wavefield at Mount St. Helens

Abstract An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

A fast amplitude preserving three-parameter 3D parabolic Radon transform and its application on multiple attenuation

Seismic wavefields propagate through three-dimensional (3D) space, and their precise characterization is crucial for understanding subsurface structures. Traditional 2D algorithms, due to their limitations, are insufficient to fully represent three-dimensional wavefields. The classic 3D Radon transform algorithm assumes that the wavefield's propagation characteristics are consistent in all directions, which often does not hold true in complex underground media. To address this issue, we present an improved 3D three-parameter Radon algorithm that considers the wavefield variation with azimuth and provides a more accurate wavefield description. However, introducing new parameters to describe the azimuthal variation also poses computational challenges. The new Radon transform operator involves five variables and cannot be simply decomposed into small matrices for efficient computation; instead, it requires large matrix multiplication and inversion operations, significantly increasing the computational load. To overcome this challenge, we have integrated the curvature and frequency parameters, simplifying all frequency operators to the same, thereby significantly improving computation efficiency. Furthermore, existing transform algorithms neglect the lateral variation of seismic amplitudes, leading to discrepancies between the estimated multiples and those in the data. To enhance the amplitude preservation of the algorithm, we employ orthogonal polynomial fitting to capture the amplitude spatial variation in 3D seismic data. Combining these improvements, we propose a fast, amplitude-preserving, 3D three-parameter Radon transform algorithm. This algorithm not only enhances computational efficiency while maintaining the original wavefield characteristics, but also improves the representation of seismic data by increasing amplitude fidelity. We validated the algorithm in multiple attenuation using both synthetic and real seismic data. The results demonstrate that the new algorithm significantly improves both accuracy and computational efficiency, providing an effective tool for analyzing seismic wavefields in complex subsurface structures.

Gas prediction in tight sandstones based on the rock-physics-derived seismic amplitude variation versus offset method

Estimating gas enrichments is a key objective in exploring sweet spots within tight sandstone gas reservoirs. However, the low sensitivity of elastic parameters to gas saturations in such formations makes it a significant challenge to reliably estimate gas enrichments using seismic methods. Through rock physical modeling and reservoir parameter analyses conducted in this study, a more suitable indicator for estimating gas enrichment, termed the gas content indicator, has been proposed. This indicator is formulated based on effective fluid bulk modulus and shear modulus and demonstrates a clear positive correlation with gas content in tight sandstones. Moreover, a new seismic amplitude variation versus offset (AVO) equation is derived to directly extract reservoir properties, such as the gas content indicator and porosity, from prestack seismic data. The accuracy of this proposed AVO equation is validated through comparison with the exact solutions provided by the Zoeppritz equation. To ensure reliable estimations of reservoir properties from partial angle-stacked seismic data, the proposed AVO equation is reformulated within the elastic impedance inversion framework. The estimated gas content indicator and porosity exhibit favorable agreement with logging data, suggesting that the obtained results are suitable for reliable predictions of tight sandstones with high gas enrichments. Furthermore, the proposed methods have the potential to stimulate the advancement of other suitable inversion techniques for directly estimating reservoir properties from seismic data across various petroleum resources.

Experiments on Landquakes Generated by Free‐Falling Granular Masses: Implications for Rockfall Impact Dynamics

AbstractThe properties of rockfalls, such as the volume and geometry of the detached rock mass as well as the number and diameter of rock fragments, greatly affect their propagated behavior, and different strategies are required to mitigate rockfall hazards. Landquakes generated by rockfalls are regarded as good proxies for those properties. To explore the free‐fall impact dynamics of rockfalls, a series of free‐fall experiments on granular masses were designed and conducted with the characteristics of the generated landquakes versus initial variables being analyzed. The results show that the particle diameter is a dominant factor affecting landquakes, with the most significant effect on the mean frequency. The impact area is another key factor, with which both the maximum seismic amplitude and radiated seismic energy have strong positive correlations, while the mean frequency of landquakes shows a weak negative correlation. The maximum seismic amplitude appears to have a poor correlation with the numbers of layers of the free‐falling granular masses, and it was proven that the maximum seismic amplitude is determined by only the lower part of the granular masses. Quantitative relationships between the number of particles and seismic signal characteristics indicate that disintegration needs to be accounted for when landquakes are used to investigate the characteristics of a rockfall. Moreover, the frequency of landquakes decreases with a decrease in the vertical velocity and increase in the horizontal velocity during granular mass collapse.

A High-Resolution subtle fault mapping in Marlin oilfield, Brazil

Discriminating geological discontinuities is important for reservoir management because they influence storage capacity and fluid flow. Structural heterogeneities, such as faults and fractures, significantly impact the compartmentalization of the reservoir’s facies and affect the flow of hydrocarbons. These structures can act as the pore space for fluid storage, flow-control conduits, or barriers that allow the accumulation of exploitable volumes of oil and gas. Understanding these structures’ spatial and temporal development can significantly impact hydrocarbon exploration. Reservoir characterization is a critical process in the oil and gas industry that aims to comprehensively understand reservoir rocks and their fluid content and distribution. It is a multidisciplinary approach that combines geological, geophysical, and engineering datasets through statistical and deterministic mathematical models. Seismic datasets play a crucial role in reservoir characterization. Interpreters can use seismic data processing and imaging to interpret geological horizons and map faults, building complex geological models. Seismic inversion is also employed to estimate acoustic impedance and other petrophysical properties. Faults and fractures mapping are highly relevant themes in reservoir characterization. Identifying subtle faults and fractures is essential in specifying hydrocarbon migration paths and identifying bypassed oil accumulations. However, traditional seismic amplitude data interpretation should consider these structures more. In this study, we apply a high-resolution fault mapping workflow to map the subtle faults system at the Marlim field in the offshore portion of the Campos Basin. As an outcome, we imaged two separate fault trends inside and outside the Marlim reservoir.

Rapid Finite-Fault Models for the 2023 Mw 7.8 Kahramanmaraş, Türkiye, Earthquake Sequence

Abstract In the immediate aftermath of devastating earthquakes such as in the 6 February 2023 Kahramanmaraş sequence in southcentral Türkiye, key stakeholders and the public demand timely and accurate earthquake information. Especially for large events, finite-fault models provide important insights into the rupture process and enable interpretation of the observed ground shaking, which can improve situational awareness and facilitate rapid assessment of future hazards. Using strong-motion waveforms recorded during the Kahramanmaraş sequence, we simulate a real-time playback and calculate how a finite-source model computed with the Finite-fault rupture Detector (FinDer) algorithm would evolve for the Mw 7.8 Pazarcık, Mw 7.6 Elbistan, and Mw 6.4 Yayladağı earthquakes. Using template matching FinDer compares observed and predicted ground-motion acceleration amplitudes to determine the orientation and spatial extent of fault rupture. We test both generic crustal and fault-specific templates from ground-motion models and rupture geometries of the east Anatolian and Çardak–Sürgü faults. In the second step, we estimate the seismic slip along the source models from the backprojection of the seismic displacement amplitudes. The algorithms achieve excellent performance for all three earthquakes, and the final source models and slip profiles available within tens of seconds of the rupture nucleation match well with models computed days to weeks after the events occurred. The temporal evolution of the source models for the Pazarcık and Elbistan earthquakes suggests that FinDer can provide insight into the rupture kinematics of large earthquakes. Cascading instrument failures as well as power and data telemetry interruptions during the Pazarcık earthquake led to an early termination of signals at a significant number of near-source stations. We show that FinDer is robust enough to cope with this type of degradation in network performance that can occur in large earthquakes, in general.

Seismic and radon signatures: A multiparametric approach to monitor surface dynamics of a hazardous 2021 rock–ice avalanche, Chamoli Himalaya

AbstractThe observation of precursory signals of the 2021 Chamoli rock–ice avalanche provides an opportunity to investigate the multidisciplinary analysis approach of rock failure. On 7 February 2021, a huge rock–ice mass detached from the Raunthi peak at Chamoli district in Uttarakhand, India. The tragic catastrophe resulted in more than 200 deaths and significant economic losses. Here, we analyse radon concentration and seismic signals to characterise the potential precursory anomalies prior to the detachment. Continuous peaks of radon anomalies were observed from the afternoon of 5 to 7 February and decreased suddenly after the event, while a cumulative number of seismic tremors and amplitude variations are more intensified ~2.30 h before the main event, indicating a static to dynamic phase change within the weak zone. This study not only characterises abnormal signals but also models the rock failure mechanisms. The analysis unveils three time‐dependent nucleation phases, physical mechanisms of signal generation and a complete scenario of physical factors that affected the degree of criticality of slope failure. The results of this study suggest gradual progression of rock cracks/joints, subsequent material creep and slip advancement acceleration preceded the final failure. Furthermore, the study highlights the importance of an early warning system to mitigate the impact of events like the 2021 Chamoli rock–ice avalanche.

Amplitude‐Dependent Energy Dissipation Inferred From Spectral Analyses of Intraslab Earthquakes

AbstractWhile the Earth (rocky planet) can be approximated as an elastic body, it also possesses anelasticity, Q−1, which causes energy dissipation (intrinsic attenuation) during deformation. Rock deformation experiments have demonstrated that Q−1 shows a marked amplitude dependency for strain of >∼10−6, but seismic‐waves‐induced strain is generally much smaller than the threshold value. Therefore, all seismological analyses have assumed amplitude‐independent attenuation. However, theoretical models of dislocation‐induced attenuation have predicted amplitude‐dependent attenuation under high temperatures and pressures equivalent to the crustal and uppermost mantle conditions. This study investigates whether attenuation is actually independent of seismic‐wave amplitudes by a systematic analysis of a large number of spectral amplitudes of co‐located earthquake pairs with different observed amplitudes. We assume the ω2 source and circular crack models to correct for the source parameters of earthquakes. The obtained results suggest that amplitude‐dependent attenuation is required to explain the observations, which contradicts a long‐standing recognition that amplitude‐dependent attenuation is negligible for the propagation of seismic waves. Our model calculation further suggests that Q−1 is proportional to the seismic amplitude, A, whereby Q−1 ∝ A0.1–0.2, which introduces attenuation variations along a given raypath, with significantly enhanced attenuation near the hypocenter. Our result highlights the need to revisit and potentially rethink current models of the physical mechanisms influencing intrinsic attenuation.

Classification and Prediction Method for Dolomite Reservoirs in the Permian Maokou Formation in the Sichuan Basin

The Lower Permian Maokou Formation in the Hechuan area is one of the key horizons for natural gas exploration in the Sichuan Basin, showing a low degree of exploration and great exploration potential. The Lower Permian gas reservoir is mainly controlled by the reservoir, therefore it is subject to the development of the Maokou Formation reservoir. However, due to the strong heterogeneity and significant variations in physical properties of the Maokou Formation dolomite reservoir in the Hechuan area, there are significant differences in seismic response characteristics of reservoirs with different porosities. Through comprehensive analysis of gas testing wells in the Maokou Formation, criteria for reservoir classification were established based on lithology, average porosity, plane porosity, and thickness of dolomite reservoirs. The effective reservoirs in the Maokou Formation were divided into two types: high-porosity reservoirs and medium-porosity reservoirs. On the basis of establishing a physical quantity model for the Lower Permian carbonate rock in the research area, the rock physical characteristics and seismic response characteristics of reservoirs with different pores are determined through model forward modeling, and appropriate methods are selected for prediction. By using the optimization processing method of gathers and wavelet decomposition to remove sidelobes, the issue of seismic abnormal response caused by lateral velocity changes in the underlying Maokou Formation has been solved. The relative strength relationship of seismic response amplitude can effectively characterize the true situation of reservoir development. Then, by using post-stack inversion and pre-stack simultaneous inversion methods, the high-porosity and medium-porosity dolomite reservoirs in the Maokou Formation are classified and predicted, thereby to effectively improve the prediction accuracy of reservoirs, ultimately providing a reliable basis for the exploration and development of the Maokou Formation in the Hechuan area in the later stage.

Basin‐scale prediction of S‐wave Sonic Logs using Machine Learning techniques from conventional logs

AbstractS‐wave velocity plays a crucial role in various applications but often remains unavailable in vintage wells. To address this practical challenge, we propose a machine learning framework utilizing an enhanced bidirectional long short‐term memory algorithm for estimating S‐wave sonic logs from conventional logs, including P‐wave sonic, gamma ray, total porosity, and bulk density. These input logs are selected based on traditional rock physics models, integrating geological and geophysical relations existing in the data. Our study, encompassing 34 wells across diverse formations in the Delaware Basin, Texas, demonstrates the superiority of machine learning models over traditional methods like Greenberg–Castagna equations, without prior geological and geophysical information. Among these machine learning models, the enhanced bidirectional long short‐term memory model with self‐attention yields the highest performance, achieving an R‐squared value of 0.81. Blind tests on five wells without prior geologic information validate the reliability of our approach. The estimated S‐wave velocity values enable the creation of a basin‐scale S‐wave velocity model through interpolation and extrapolation of these prediction models. Additionally, the bidirectional long short‐term memory model excels not only in predicting S‐wave velocity but also in estimating S‐wave reflectivity for seismic amplitude variation with offset applications in exploration seismology. In conclusion, these S‐wave velocity estimates facilitate the prediction of further elastic properties, aiding in the comprehension of petrophysical and geomechanical property variations within the basin and enhancing earthquake hypocentral depth estimation.

1300-m-High gas plume from pockmarks in the north Nansha waters, South China sea

The northeastern region of Nansha Waters, located in the southern margin of the South China Sea, is abundant in gas resources. However, it remains unclear if the central and northern areas of Nansha Waters possess gas resources and if there are any associated leaks. The Jiuzhang Basin is located in the northern Nansha Waters. Pre-stack depth migration, analysis of seismic attributes and multibeam data processing were conducted to image the seabed topography and backscatter intensity, as well as the seismic structure of the water and sedimentary layers in the Jiuzhang Basin. The result reveals a presence of large plumes, approximately 1500 m in width and up to 1300 m in height, from a pockmark. The plume is characterized by a strong seismic amplitude featuring chaotic reflections, a high instantaneous frequency (100–150 Hz), and a strong instantaneous amplitude. These seismic features commonly appear in hydrocarbon-bearing seawater, indicating that the plume is composed of gas. The relationship between the volume fraction of bubbles in the gas plumes and the amplitude of multi-channel seismic waves suggests that the volume fraction of bubbles within the gas plume can reach up to 5.0 × 10−3. Meanwhile, a new pockmark with a length of 248 m, a width of 90 m, and a depth of 20.4 m, and an abrupt 22 m drop in water depth of an existing pockmark are observed in the gas plume area. These observations, combined with the gas plume revealed by the multi-channel seismic collected in 2013 and multibeam collected in 2020, suggest a very strong and recently active gas leaks in the pockmark area of the northern Jiuzhang Basin. Furthermore, in the pockmark area, there are gas-bearing formations with weak reflections, and mud diapirs charactered by low seismic velocity, low amplitude chaotic reflections, and pulled up reflections on both sides, and faults extend to the seabed. These observations suggest that gas plume originate from deep strata of the northern Jiuzhang Basin. Large gas plumes and mud diapirs are often indicative of the existence of deep gas reservoirs. Therefore, the presence of such active gas plumes and mud diapirs in the northern region of the Jiuzhang Basin signifies the presence of abundant gas resources.

Comparing the Properties of Well Logs to Seismic Amplitude Analysis to Characterize the Reservoir of the KD Oilfield in Niger Delta, Nigeria

Well logs and data are integrated in most reservoir exploration to provide a reservoir description, the depth and width of the holding subsurface structure are two reservoir parameters that must be understood in order to quantify viable hydrocarbon reservoirs, these metrics are important because they offer suitable contributions for measuring reservoir volume, the objective is to cross-plot well logs with seismic amplitude to provide reservoir characteristics, map reservoir pay zones and determine reservoir depth and thickness, the method applied in this research were collection of set of data consisting of ten wells information cutting across the Agbada formation and 3D post-stack seismic data acquired from the Niger Delta, set of well log suites and software (Word, Excel and Petrel 2014 version) used for analysis, the result analyzed after various cross plots comparing the amplitudes with well logs reveals several trends were discovered, this indicates that higher amplitudes correspond to higher concentrations of hydrocarbons or water based on peak amplitude versus neutron log cross-plot for well KR-1. This study examines the relationship between well logs and amplitude peaks with troughs, as well as potential changes in log signal resulting from variations in seismic amplitude when assessing reservoir sand.

Reservoir Characterization of KD Field Located Onshore of Akwa-Ibom State, Niger Delta, Nigeria Using Well Logs Property Comparison to Seismic Amplitude Analysis

Characterization of reservoirs defines how well they can generate and store hydrocarbon, reservoir parameters are used to determine the behavior of reservoir fluids under various conditions and to identify the best production practices that can optimize hydrocarbon production, the objective is to identify and map reservoir sand zones using basic logs, provide information about the reservoir depth and thickness using well log data, interpret faults and sealing system using the seismic data and generate amplitude maps and model basic reservoir facies, the method applied was sectioned into three distinct parts that encompasses qualitative well log analysis, seismic interpretation and cross plotting of basic well logs versus amplitude for reservoir characterization, the well data analyzed shows well number 5, 6 and 8 have incomplete GR logs resulting in missing lithology while well 10 have some missing logs, such as, density (DEN) log, neutron porosity (NEU) log and porosity (POR) log, five logs needed for this study were available in the remaining complete wells. These logs include gamma ray (GR) log, resistivity (RES) log, density (DEN) log, neutron porosity (NEU) log, and porosity (POR) log. This study looks at a possible change in log signal due to an increase or decrease in seismic amplitude in characterizing a reservoir sand and how logs are influence with Amplitude peaks and trough.

Use of resistivity and density borehole image logs to identify and distribute facies in the Pikka unit — A case study from the Nanushuk Formation, North Slope, Alaska

Over the past few decades, the North Slope of Alaska has yielded several major hydrocarbon discoveries in the deltaic topsets of the Brookian Nanushuk Formation. Together, the Nanushuk topsets and genetically related foreset and bottomset beds of the Torok Formation comprise part of a giant clinothem system that prograded across the Colville Foreland Basin during the lower Cretaceous (the Aptian through the Cenomanian). The Nanushuk Topset Play contains stratigraphically trapped hydrocarbons within multiple fairways trending roughly north to south along the basin’s extent. Inside the Pikka Unit, the Nanushuk Formation represents a shelf-edge delta and attached shoreface system, defined primarily from sedimentological interpretations of conventional cores integrated with 3D seismic reflection data and well-log motifs. In this study, we build a structural and stratigraphic framework for the Nanushuk reservoir across the Pikka Unit through integrated observations of the borehole image logs, seismic data, and core. In the described workflow, the boundaries delimiting dip zones interpreted in the wells are correlated to higher order stratigraphic cycles within the Nanushuk clinothem zone. These boundaries are observed in seismic profiles and provide the basis of the structural framework. In addition, core facies interpretations are calibrated to CT-scan data to establish each facies character and vertical distribution within the wells. In the absence of core data, core-calibrated borehole images are instrumental in guiding sedimentological and structural interpretations. This extrapolation is dependent on the image log quality and resolution. The calibrated facies interpretations are applied to wells without core data and guided by seismic amplitude interpretation to extrapolate reservoir presence and distribution across the area.

Sedimentary architecture of thin-layer beach bar sand bodies in the G oilfield, Niger

The G oilfield in the southeastern Termit Basin of Niger is characterized by thin-bed beach bar deposits exhibiting strong reservoir heterogeneity and suboptimal production efficiency, necessitating internal structural dissection of the beach bar sand bodies. Employing a well-seismic integration approach, we systematically dissect the architecture of these sand bodies layer by layer to determine their spatial distribution. A classification scheme for beach bar architecture is proposed, with core and log data analysis identifying essential architectural elements and their logging responses. Seismic amplitudes, thin bed delineation, frequency decomposition inversion attributes, and attribute fusion technology delineate the architectural boundaries. Integrating five indicators from four-level architectural recognition at wellbores—shallow lake mudstone appearance, bar margin/beach microfacies occurrence, logging curve morphology differences, beach bar thickness variations, and elevation differences between adjacent bars—enables detailed dissection of the beach bar architecture, corroborated by connectivity analysis. In the study area, beach bar distribution primarily develops in two modes: vertical stacking (accumulation of multiple main bars from different episodes) and isolated (stable mudstone interlayers between main bar sand bodies appearing relatively isolated). This research provides a basis for dissecting beach bar architecture reservoirs under sparse well conditions.

Forecasting the onset of volcanic eruptions using the increase in seismicity during magma ascent

The Failure Forecast Method (FFM) has been used to forecast the onset of volcanic eruptions with varying degrees of success. The method involves fitting its empirical equation to precursory observables, e.g., seismic data. Current models explaining the empirical equation assume that the seismic observables used (e.g., seismic event-rate, seismic moment-rate or seismic energy-rate) are related to accelerated fracture growth. In these models, such fracture growth results in an apparent increase in seismic activity. In this study, however, we propose an alternative explanation for the increase in seismicity culminating in volcanic eruptions. We suggest that the increase in seismicity, particularly in the seismic amplitudes, can result from ascending seismic sources towards the surface and that the increase can be used for eruption forecasts using the FFM. Our argument is based on the need for magma to migrate towards the surface to produce an eruption.To test our argument, we used synthetic and field/real data and fitted them with the FFM model. The synthetic data is generated using the body wave attenuation equation, and the real data is from the Piton de la Fournaise eruption on January 2nd, 2010. We used seismic amplitudes and seismic amplitude ratios as the precursory observables and used a combination of Simulated Annealing and Bayesian Inversion algorithms to obtain the best-fit parameters of the FFM model and their posterior distributions. For the synthetic cases, we found that the method gave reasonable estimates of failure or eruption time, suggesting that the seismic amplitudes and the seismic amplitude ratios can be used for forecasting using the FFM. For the field/real cases, the results gave information about the timing of magma chamber failure, magma propagation, and the eruption onset. Our results highlight an example of modelling volcano monitoring data with the FFM that gives reasonable estimates of eruption time even if the data may not fully represent failure processes. We suggest that the result may still be reasonable owing to the empirical and asymptotic nature of the FFM combined with the underlying physics of magma and seismic migrations.

Seismic amplitude inversion based on a new PP-wave reflection coefficient approximation equation for vertical transversely isotropic media

We develop a method for the simultaneous amplitude-variation-with-offset or -angle inversion of anisotropic parameters for transversely isotropic media with vertical axis of symmetry (VTI media). First, we introduce a nonlinear PP-wave reflection coefficient approximation equation in terms of only P- and S-wave impedances for isotropic elastic media. Then, by replacing the isotropic part of Rüger’s equation with this equation, we obtain a new PP-wave reflection coefficient approximation equation called the ASI Rüger equation for VTI media. To invert the parameters for VTI media based on the ASI Rüger equation, we adopt the Bayesian generalized linear inversion (GLI) method, a combination of GLI and Bayesian linear inversion, in which the noise and model perturbation are assumed to conform to the zero-mean Gaussian distribution. Compared with Rüger’s equation, the ASI Rüger equation decreases the trade-off between the parameters and reduces the ill-posedness of the inverse problem. The synthetic and field data tests demonstrate the feasibility of our method for inverting VTI media parameters (the vertical P-wave impedance, the vertical S-wave impedance, and Thomsen’s parameters [Formula: see text] and [Formula: see text]).