-
双重孔隙介质模型考虑了岩石非均质性诱发的介观流对弹性波频散和衰减的影响,在非均质储层地震资料定量解释中取得良好的应用效果.基于双重孔隙介质理论模型,利用数值算法开展弹性波模拟工作不仅可以直观显示波的传播特征,同时也为后期地震反演成像工作奠定基础.本文基于Santos-Rayleigh部分饱和双重孔隙介质模型,利用交错网格有限差分算法模拟并分析了双重孔隙介质中的波场快照和波形曲线。采用Zener粘弹性模型近似表征介观流机制,结果表明Zener模型能够较好反映快纵波传播特征,却无法表征慢纵波P3波在低频段的衰减特征.利用时间分裂法解决波动方程的刚性问题,提高计算效率.在利用解析解验证了有限差分算法正确性的基础上,模拟了均匀介质和分层介质中的波场快照和波形曲线,清晰直观的反映了快纵波在地震频段的强衰减特征,与双重孔隙理论模型预测结果一致.模拟结果有助于进一步理解非均匀部分饱和孔隙介质中的弹性波传播特征.Double-porosity poroelastic models, which account for the effect of mesoscopic flow in heterogeneous rocks on wave dispersion and attenuation, are useful for quantitative seismic interpretation. Wavefield simulation based on double-porosity models not only helps visualize the propagation characteristics of the elastic waves but also lays the foundation for seismic imaging. In this paper, we perform wavefield simulation and analysis based on the Santos-Rayleigh model which incorporates mesoscopic and global flow in partially-saturated double-porosity media. Specifically, the mesoscopic flow mechanism is represented with a Zener viscoelastic model. The comparison shows that the Zener model can accurately capture the propagation characteristics of fast P-wave, but fails to represent the attenuation characteristics of slow P3 wave in the low-frequency band. It suggests that Zener viscoelastic model and slow wave modes follow different mechanisms. Then staggered grid finite-difference method is used to simulate wave propagation in double-porosity media, and the stiff problem is solved with a time-splitting algorithm, which can significantly improve computational efficiency. Based on above methods, the correctness of our algorithm is verified with derived analytical solution for a P-wave source in a uniform partially saturated poroelastic media. Analytical and numerical solutions are in good agreement and mean error is 0.33%. We provide some examples of wavefield snapshots and seismograms in homogeneous and layered heterogeneous media at seismic and ultrasonic frequencies. Simulation results demonstrate the strong attenuation of fast P-wave and no change of S-wave in the seismic band due to mesoscopic flow mechanism, which is consistent with the theoretical predictions of double-porosity model. Moreover, energy of fast P-wave is concentrated in solid phase while slow waves are stronger in fluid phases. This work contributes to the understanding of broadband elastic wave propagation in heterogeneous partially saturated porous media and can be applied in the reservoir imaging with broadband geophysical data.
-
Keywords:
- double-porosity media /
- finite difference /
- dispersion and attenuation
-
[1] Biot M A 1956 J. Acoust. Soc. Am. 28179
[2] Santos J E, Corberó J M, Douglas J 1990 J. Acoust. Soc. Am. 871428
[3] Liu L, Zhang X, Wang X 2022 J. Theor. Comp. Acout. 302150002
[4] Berryman J G, Wang H F 2000 Int J Rock Mech Min 3763
[5] Huang J, Yang D, He X, Chang Y 2023 Geophysics 88 T121
[6] Pride S R, Berryman J G, Harris J M 2004 J. Geophys. Res. 109 B01201
[7] Zheng P, Ding B, Sun X 2017 Int J Rock Mech Min 91104
[8] Ba J, Carcione J M, Nie J X 2011 J. Geophys. Res. 116 B06202
[9] Sun W, Ba J, Carcione J M 2016 Geophys. J. Int. 20522
[10] Ba J, Xu W, Fu L Y, Carcione J M, Zhang L 2017 J. Geophys. Res. Solid Earth 1221949
[11] Shi Z Q, He X, Liu L, Chen D H, Wang X M 2023 Acta Phys. Sin. 72069101(in Chinese) [石志奇,何晓,刘琳,陈德华,王秀明2023物理学报72069101]
[12] Shi Z, He X, Chen D, Wang X 2024 Geophys. J. Int. 2361172
[13] Ba J, Nie J X, Cao H, Yang H Z 2008 Geophys. Res. Lett. 35 L04303
[14] Liu X, Greenhalgh S, Zhou B 2009 Geophys. J. Int. 178375
[15] Liu X, Greenhalgh S 2019 Geophysics 84 WA59
[16] Wang E, Carcione J M, Ba J 2019 Geophysics 84 WA11
[17] Carcione J M 2015 Wave fields in real media: wave propagation in anisotropic, anelastic, porous and electromagnetic media (Amsterdam Boston [etc.]: Elsevier)
[18] Jiang Y, Gao Y, Cheng Q, Song Y 2023 Geophys. J. Int. 235970
[19] Virieux J 1986 Geophysics. 51889
[20] Wenzlau F, Müller T M 2009 Geophysics 74 T55
[21] Guan W, Hu H 2011 Commun. Comput. Phys. 10695
[22] Kong L Y, Wang Y B, Yang H Z 2013 Acta Phys. Sin. 62139101(in Chinese) [孔丽云,王一博,杨慧珠2013物理学报62139101]
[23] Liu C, Luo Y Q 2023 Chinese J. Geophys. 663840[刘财,罗玉钦2023地球物理学报663840]
[24] Zhao H, Wang X 2008 Sci. China Ser. G-Phys. Mech. Astron. 51723
[25] Deng W, Fu L Y, Wang Z, Hou W, Han T 2023 Geophys. J. Int. 2351218
[26] Carcione J M, Quiroga-Goode G 1995 J Comput Acoust 3261
[27] Komatitsch D, Martin R 2007 Geophysics 72 SM155
计量
- 文章访问数: 163
- PDF下载量: 6
- 被引次数: 0