Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study fo...Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.展开更多
Forward modeling of seismic wave propagation is crucial for the realization of reverse time migration(RTM) and full waveform inversion(FWI) in attenuating transversely isotropic media. To describe the attenuation and ...Forward modeling of seismic wave propagation is crucial for the realization of reverse time migration(RTM) and full waveform inversion(FWI) in attenuating transversely isotropic media. To describe the attenuation and anisotropy properties of subsurface media, the pure-viscoacoustic anisotropic wave equations are established for wavefield simulations, because they can provide clear and stable wavefields. However, due to the use of several approximations in deriving the wave equation and the introduction of a fractional Laplacian approximation in solving the derived equation, the wavefields simulated by the previous pure-viscoacoustic tilted transversely isotropic(TTI) wave equations has low accuracy. To accurately simulate wavefields in media with velocity anisotropy and attenuation anisotropy, we first derive a new pure-viscoacoustic TTI wave equation from the exact complex-valued dispersion formula in viscoelastic vertical transversely isotropic(VTI) media. Then, we present the hybrid finite-difference and low-rank decomposition(HFDLRD) method to accurately solve our proposed pure-viscoacoustic TTI wave equation. Theoretical analysis and numerical examples suggest that our pure-viscoacoustic TTI wave equation has higher accuracy than previous pure-viscoacoustic TTI wave equations in describing q P-wave kinematic and attenuation characteristics. Additionally, the numerical experiment in a simple two-layer model shows that the HFDLRD technique outperforms the hybrid finite-difference and pseudo-spectral(HFDPS) method in terms of accuracy of wavefield modeling.展开更多
The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and...The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.展开更多
基金financially supported by the Natural Science Foundation of Sichuan Province(Grant Nos.2023NSFSC0767 and2024NSFSC0809)the China Postdoctoral Science Foundation(Grant No.2024MF750281)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20230326)。
文摘Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+1 种基金the Funds of Creative Research Groups of China(No.41821002)National Natural Science Foundation of China Outstanding Youth Science Fund Project(Overseas)(No.ZX20230152)。
文摘Forward modeling of seismic wave propagation is crucial for the realization of reverse time migration(RTM) and full waveform inversion(FWI) in attenuating transversely isotropic media. To describe the attenuation and anisotropy properties of subsurface media, the pure-viscoacoustic anisotropic wave equations are established for wavefield simulations, because they can provide clear and stable wavefields. However, due to the use of several approximations in deriving the wave equation and the introduction of a fractional Laplacian approximation in solving the derived equation, the wavefields simulated by the previous pure-viscoacoustic tilted transversely isotropic(TTI) wave equations has low accuracy. To accurately simulate wavefields in media with velocity anisotropy and attenuation anisotropy, we first derive a new pure-viscoacoustic TTI wave equation from the exact complex-valued dispersion formula in viscoelastic vertical transversely isotropic(VTI) media. Then, we present the hybrid finite-difference and low-rank decomposition(HFDLRD) method to accurately solve our proposed pure-viscoacoustic TTI wave equation. Theoretical analysis and numerical examples suggest that our pure-viscoacoustic TTI wave equation has higher accuracy than previous pure-viscoacoustic TTI wave equations in describing q P-wave kinematic and attenuation characteristics. Additionally, the numerical experiment in a simple two-layer model shows that the HFDLRD technique outperforms the hybrid finite-difference and pseudo-spectral(HFDPS) method in terms of accuracy of wavefield modeling.
基金financially supported by the Sichuan Science and Technology Program (Grant No. 2023ZYD0158)the National Natural Science Foundation of China (Grant Nos. 42304147 and 42304076)。
文摘The main objective is to optimize the development of shale gas-rich areas by predicting seismic sweet spot parameters in shale reservoirs. We systematically assessed the fracture development, fracture gas content, and rock brittleness in fractured gas-bearing shale reservoirs. To better characterize gas-bearing shale reservoirs with tilted fractures, we optimized the petrophysical modeling based on the equivalent medium theory. Based on the advantages of shale petrophysical modeling, we not only considered the brittle mineral fraction but also the combined effect of shale porosity, gas saturation, and total organic carbon(TOC) when optimizing the brittleness index. Due to fractures generally functioning as essential channels for fluid storage and movement, fracture density and fracture fluid identification factors are critical geophysical parameters for fractured reservoir prediction. We defined a new fracture gas indication factor(GFI) to detect fracture-effective gas content. A new linear PP-wave reflection coefficient equation for a tilted transversely isotropic(TTI) medium was rederived, realizing the direct prediction of anisotropic fracture parameters and the isotropic elasticity parameters from offset vector tile(OVT)-domain seismic data. Synthetic seismic data experiments demonstrated that the inversion algorithm based on the L_P quasinorm sparsity constraint and the split-component inversion strategy exhibits high stability and noise resistance. Finally, we applied our new prediction method to evaluate fractured gas-bearing shale reservoirs in the Sichuan Basin of China, demonstrating its effectiveness.
