Sensitivity of seismic attenuation and dispersion to dynamic elastic interactions of connected fractures: Quasi-static finite element modeling study

Prediction of seismic attenuation and dispersion that are inherently sensitive to hydraulic and elastic properties of the medium of interest in the presence of mesoscopic fractures and pores,is of great interest in the characterization of fractured formations.This has been very difficult,however,considering that stress interactions between fractures and pores,related to their spatial distributions,tend to play a crucial role on affecting overall dynamic elastic properties that are largely unexplored.We thus choose to quantitatively investigate frequency-dependent P-wave characteristics in fractured porous rocks at the scale of a representative sample using a numerical scale-up procedure via performing finite element modelling.Based on 2-D numerical quasi-static experiments,effects of fracture and fluid properties on energy dissipation in response to wave-induced fluid flow at the mesoscopic scale are quantified via solving Biot's equations of consolidation.We show that numerical results are sensitive to some key characteristics of probed synthetic rocks containing unconnected and connected fractures,demonstrating that connectivity,aperture and inclination of fractures as well as fracture infills exhibit strong impacts on the two manifestations of WIFF mechanisms in the connected scenario,and on resulting total wave attenuation and phase velocity.This,in turn,illustrates the importance of these two WIFF mechanisms in fractured rocks and thus,a deeper understanding of them may eventually allow for a better characterization of fracture systems using seismic methods.Moreover,this presented work combines rock physics predictions with seismic numerical simulations in frequency domain to illustrate the sensitivity of seismic signatures on the monitoring of an idealized geologic CO_(2) sequestration in fractured reservoirs.The simulation demonstrates that these two WIFF mechanisms can strongly modify seismic records and hence,indicating that incorporating the two energy dissipation mechanisms in the geophysical interpretation can potentially improving the monitoring and surveying of fluid variations in fractured formations.

Modeling the efects of fracture infll on frequency‑dependent anisotropy and AVO response of a fractured porous layer

In a fractured porous hydrocarbon reservoir,wave velocities and refections depend on frequency and incident angle.A proper description of the frequency dependence of amplitude variations with ofset(AVO)signatures should allow efects of fracture inflls and attenuation and dispersion of fractured media.The novelty of this study lies in the introduction of an improved approach for the investigation of incident-angle and frequency variations-associated refection responses.The improved AVO modeling method,using a frequency-domain propagator matrix method,is feasible to accurately consider velocity dispersion predicted from frequency-dependent elasticities from a rock physics modeling.And hence,the method is suitable for use in the case of an anisotropic medium with aligned fractures.Additionally,the proposed modeling approach allows the combined contributions of layer thickness,interbedded structure,impedance contrast and interferences to frequency-dependent refection coefcients and,hence,yielding seismograms of a layered model with a dispersive and attenuative reservoir.Our numerical results show bulk modulus of fracture fuid signifcantly afects anisotropic attenuation,hence causing frequencydependent refection abnormalities.These implications indicate the study of amplitude versus angle and frequency(AVAF)variations provides insights for better interpretation of refection anomalies and hydrocarbon identifcation in a layered reservoir with vertical transverse isotropy(VTI)dispersive media.

Viscoacoustic generalized screen propagator in constant-Q model

When seismic waves propagate through the geological formation,there is a significant loss of energy and a decrease in imaging resolution,because of the viscoacoustic properties of subsurface medium.This profoundly impacts seismic wavefield propagation,imaging and interpretation.To accurately image the true structure of subsurface medium,the consensus among geophysicists is to no longer treat subsurface medium as ideal homogeneous medium,but rather to incorporate the viscoacoustic properties of subsurface medium.Based on the generalized screen propagator using conventional acoustic wave equation(acoustic GSP),our developed method introduces viscoacoustic compensation strategy,and derives a one-way wave generalized screen propagator based on time-fractional viscoacoustic wave equation(viscoacoustic GSP).In numerical experiments,we conducted tests on two-dimensional multi-layer model and the Marmousi model.When comparing with the acoustic GSP using the acoustic data,we found that the imaging results of the viscoacoustic GSP using the viscoacoustic data showed a significant attenuation compensation effect,and achieved imaging results for both algorithms were essentially consistent.However,the imaging results of acoustic GSP using viscoacoustic data showed significant attenuation effects,especially for deep subsurface imaging.This indicates that we have proposed an effective method to compensate the attenuated seismic wavefield.Our application on a set of real seismic data demonstrated that the imaging performance of our proposed method in local areas surpassed that of the conventional acoustic GSP.This suggests that our proposed method holds practical value and can more accurately image real subsurface structures while enhancing imaging resolution compared with the conventional acoustic GSP.Finally,with respect to computational efficiency,we gathered statistics on running time to compare our proposed method with conventional Q-RTM,and it is evident that our method exhibits higher computational efficiency.In summary,our proposed viscoacoustic GSP method takes into account the true properties of the medium,still achieves migration results comparable to conventional acoustic GSP.

The analysis of frequency-dependent characteristics for fluid detection: a physical model experiment

According to the Chapman multi-scale rock physical model, the seismic response characteristics vary for different fluid-saturated reservoirs. For class I AVO reservoirs and gas-saturation, the seismic response is a high-frequency bright spot as the amplitude energy shifts. However, it is a low-frequency shadow for the Class III AVO reservoirs saturated with hydrocarbons. In this paper, we verified the high-frequency bright spot results of Chapman for the Class I AVO response using the frequency-dependent analysis of a physical model dataset. The physical model is designed as inter-bedded thin sand and shale based on real field geology parameters. We observed two datasets using fixed offset and 2D geometry with different fluid- saturated conditions. Spectral and time-frequency analyses methods are applied to the seismic datasets to describe the response characteristics for gas-, water-, and oil-saturation. The results of physical model dataset processing and analysis indicate that reflection wave tuning and fluid-related dispersion are the main seismic response characteristic mechanisms. Additionally, the gas saturation model can be distinguished from water and oil saturation for Class I AVO utilizing the frequency-dependent abnormal characteristic. The frequency-dependent characteristic analysis of the physical model dataset verified the different spectral response characteristics corresponding to the different fluid-saturated models. Therefore, by careful analysis of real field seismic data, we can obtain the abnormal spectral characteristics induced by the fluid variation and implement fluid detection using seismic data directly.