AVO analysis for high amplitude anomalies using 2D pre-stack seismic data

Amplitude variation with offset (AVO) analysis is an 1 efficient tool for hydrocarbon detection and identification of elastic rock properties and fluid types. It has been applied in the present study using reprocessed pre-stack 2D seismic data (1992, Caulerpa) from north-west of the Bonaparte Basin, Australia. The AVO response along the 2D pre-stack seismic data in the Laminaria High NW shelf of Australia was also investigated. Three hypotheses were suggested to investigate the AVO behaviour of the amplitude anomalies in which three different factors; fluid substitution, porosity and thickness (Wedge model) were tested. The AVO models with the synthetic gathers were analysed using log information to find which of these is the controlling parameter on the AVO analysis. AVO cross plots from the real pre-stack seismic data reveal AVO class IV (showing a negative intercept decreasing with offset). This result matches our modelled result of fluid substitution for the seismic synthetics. It is concluded that fluid substitution is the controlling parameter on the AVO analysis and therefore, the high amplitude anomaly on the seabed and the target horizon 9 is the result of changing the fluid content and the lithology along the target horizons. While changing the porosity has little effect on the amplitude variation with offset within the AVO cross plot. Finally, results from the wedge models show that a small change of thickness causes a change in the amplitude; however, this change in thickness gives a different AVO characteristic and a mismatch with the AVO result of the real 2D pre-stack seismic data. Therefore, a constant thin layer with changing fluids is more likely to be the cause of the high amplitude anomalies.

TIME-LAPSE POROSITY AND VELOCITY ANALYSIS USING ROCK PHYSICS MODELS IN NIGER DELTA, NIGERIA

The analysis of 3-D and time-lapse seismic data in Isomu Field has offered the dynamic characterization of the reservoir changes. The changes were analyzed using fluid substitution and seismic velocity models. The results of the initial porosity of the reservoirs was 29.50% with water saturation value of12%.The oil and gas maintained saturation values of 40% and 48% with average compressional and shear wave velocities of 2905m/s and 1634m/s respectfully. However, in fluid substitution modelling, the results reflect a change in fluid properties where average gas and oil saturation assume a new status of 34% and 24% which indicates a decrease by 14% and 16% respectively. The average water saturation increases by 30% with an average value of 42%. The decrease in hydrocarbon saturation and increase in formation water influence the porosity. Thus, porosity decreased by 4.16% which probably arose from the closure of the aspect ratio crack due to pressure increase.

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

Can Fluid-Substitution Models Continue to Ignore the Complex Physico-Chemical Properties of Crude Oil?

The mechanical nature of fluid-substitution models has always been recognized as a major cause of their limited predictive power. Saturants, for instance, are typically treated as simple fluids characterized only by their densities, viscosities, and moduli of elasticity; their chemistry is just ignored, even when that fluid is crude oil. However, crude oil is a complex mixture of several thousands of organic compounds characterized by a variety of molecular weights, polarities, and polarizabilities, and the response of its rheological behavior to acoustic wave propagation is difficult to predict, especially when it resides inthe pore space of rocks. Here, we report ultrasonic-velocity measurements performed on carbonate core plugs saturated with brine and with a light crude oil that are mechanically similar (i.e., having comparable densities, viscosities, and moduli of elasticity) and that show a significant and consistent excess of hardening when the saturant is oil. Dispersion and wettability are excluded as explanations for the data. We hypothesize that asphaltene aggregation and adsorption, as well as paraffin-wax crystallization (and possibly volumetric expansion), combine to cause crude oil to exhibit a dilatant-like behavior within the pore space of carbonates at ultrasonic frequencies. Roughly speaking, the observed effect would be similar to the hardening of oobleck at high deformation rates. This hypothesis could betested in the future by an adequate combination of high-resolution imaging and microfluidic setups. This and similar studies would be beneficial in developing physical fluid-substitution models with a more consistent predictive power.

NMR Fluid Substitution–A New Method of Reconstructing T2 Distributions Under Primary Drainage and Imbibition Conditions

Most nuclear magnetic resonance (NMR)-based petrophysics models, such as pore structure characterization and permeability prediction, were developed using T2 distributions measured at fully water-saturated conditions (i.e., Sw = 1). The downhole implementation of those models across the hydrocarbon zones is disputable due to partial saturation (Sw < 1) conditions; hence, a correction to such effects on T2 distributions is required. This paper provides a critical review of the fluid substitution methods currently available in the industry and presents an improved method for enhanced formation evaluation. In the new method presented, an effective irreducible water saturation model is used to account for the pore structure and capillary pressure effects, which were barely considered by the currently available NMR fluid substitution methods. For water-wet reservoir rocks, the typical NMR T2 distribution at the partial saturation condition displays a clear separation between the wetting and nonwetting phases. The water phase can be classified as irreducible and movable fluid volumes. Then, using a T2 mapping relationship and a total porosity constraint, the T2 distribution of movable water at Sw < 1 is shifted and amplified to determine the T2 distribution of movable water at Sw = 1. To validate the new method, NMR measurements were conducted on sandstone samples at Sw = 1 as well as Sw < 1. The reconstructed T2 distribution at Sw = 1 was compared with the measured T2 distribution at Sw = 1. Results showed that the reconstructed T2 distribution matched very well with the T2 distribution measured at Sw = 1, confirming the robustness of the new technique. Parameters used in the reconstruction methodology are observed to be a good indicator of pore connectivity. During desaturation, the water T2 in large pores shifts to a shorter T2 because of the enhanced surface relaxation as the water volume decreases while the surface area remains constant. Therefore, the amplitude at the short T2 increases. The increased amplitude was remapped to large pores in reconstructing T2 spectra of full saturation.

Rock property modelling and sensitivity analysis for hydrocarbon exploration in OSSY field, Niger Delta Basin

The application of quantitative interpretation techniques for hydrocarbon prospect evaluation from seismic has become so vital. The effective employment of these techniques is dependent on several factors: the quality of the seismic and well data, sparseness of data, the physics of rock, lithological and structural complexity of the field. This study adopts reflection pattern, amplitude versus offset (AVO), Biot–Gassmann fluid substitution and cross-plot models to understand the physics of the reservoir rocks in the field by examining the sensitivity of the basic rock properties; P-wave velocity, S-wave velocity and density, to variation in lithology and fluid types in the pore spaces of reservoirs. This is to ascertain the applicability of quantitative seismic interpretation techniques to explore hydrocarbon prospect in the studied field. The results of reflection pattern and AVO models revealed that the depth of interest is dominated by Class IV AVO sands with a high negative zero offset reflectivity that reduces with offset. The AVO intercept versus gradient plot indicated that both brine and hydrocarbon bearing sands can be discriminated on seismic. Fluid substitution modelling results revealed that the rock properties will favourably respond to variation in oil saturation, but as little as 5% gas presence will result in huge change in the rock properties, which will remain constant upon further increments of gas saturation, thereby making it difficult to differentiate between economical and sub-economical saturations of gas on seismic data. Rock physics cross-plot models revealed separate cluster points typical of shale presence, brine sands and hydrocarbon bearing sands. Thus, the response of the rock properties to the modelling processes adopted favours the application of quantitative interpretation techniques to evaluate hydrocarbon in the field.