吉林大学学报(地球科学版) ›› 2018, Vol. 48 ›› Issue (3): 900-908.doi: 10.13278/j.cnki.jjuese.20170282
胡宁, 刘财
Hu Ning, Liu Cai
摘要: 相对于整数阶导数,分数阶微分算子可以更简洁地描述具有历史依赖性和空间全域相关性的复杂力学和物理过程。但是对分数阶波动方程进行数值模拟,计算量和存储量均较大,尤其对长时间或大计算域的模拟更是如此。文中给出了3种计算方法:全局记忆法、短时记忆法、自适应记忆法,并将这3种方法应用于含黏滞流体黏弹双相VTI (横向各向同性)介质分数阶波传播方程正演。通过对比3种方法的模拟精度、计算时间及占用内存发现:虽然短时记忆法可以通过设置短时记忆长度来调整计算时间与所占内存,但是短时记忆长度越短,精度越差;而自适应记忆法在保证精度的前提下,是短时记忆法与全局记忆法在计算时间与占用内存两方面的折衷。最后对各方法的利弊进行总结,为后续正演模拟及新的分数阶数值算法开发提供方法上的参考。在正演过程中,不仅要使所建模型更贴近实际地下介质,还需对选取的数值算法在计算时间、计算存储量和精度之间进行利弊权衡,以得到一个比较合理的数值算法。
中图分类号:
| [1] Biot M A. General Theory of Three-Dimensional Con-solidation[J]. J Appl Phys, 1941, 12:155-164. [2] Dvorkin J, Nur A. Dynamic Poroelasticity:A Unified Model with the Squirt and the Biot Mechanisms[J]. Geophysics, 1993,58:524-533. [3] Diallo M S, Appel E. Acoustic Wave Propagation in Saturated Porous Media:Reformulation of the Biot/Squirt Flow Theory[J]. J Appl Geophys, 2000, 44:313-325. [4] Liu Q R, Katsube N J. The Discovery of a Second Kind of Rotational Wave in Fluid Filled Porous Material[J]. The Journal of the Acoustical Society of America, 1990, 88(2):1045-1053. [5] Kjartansson E. Constant-Q Wave Propagation and At-tenuation[J]. Journal of Geophysical Research, 1979, 84:4737-4748. [6] Carcione J M. Time-Domain Modeling of Constant-Q Seismic Waves Using Fractional Derivatives[J]. Pure and Applied Geophys, 2002, 159:1719-1736. [7] Podlubny. Fractional Differential Equations[Z]. San Diego:Academic Press, 1999. [8] Sprouse B P. Computational Efficiency of Fractional Diffusion Using Adaptive Time Step Memory[Z]. Dissertations & Theses-Gradworks, 2010. [9] Caputo M, Mainardi F. Linear Models of Dissipation in Anelastic Solids[J]. La Rivista del Nuovo Cimento, 1971, 1(2):161-198. [10] 林孔容. 关于分数阶导数的几种不同定义的分析与比较[J]. 闽江学院学报, 2003, 24(5):3-6. Lin Kongrong, Analysis and Comparision of Different Definition About Fractional Integrals and Derivatives[J]. Journal of Minjiang University, 2003,24(5):3-6. [11] 卢明辉,巴晶,杨慧珠. 含粘滞流体孔隙介质中的弹性波[J].工程力学,2009,26(5):36-40. Lu Minghui,Ba Jing,Yang Huizhu. Propagation of Elastic Waves in a Viscous Fluid-Saturated Prorous Solid[J]. Engineering Mechanics, 2009, 26(5):36-40. [12] Yang D H, Zhang Z J. Poroelastic Wave Equation Including the Biot/Squirt Mechanism and the Solid/Fluid Coupling Anisotropy[J]. Wave Motion, 2002, 35(3):223-245. [13] 刘财,兰慧田,郭智奇,等. 基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析[J].地球物理学报, 2013, 56(10):3461-3473. Liu Cai, Lan Huitian, Guo Zhiqi, et al. Pseudo-Spectral Modeling and Feature Analysis of Wave Propagation in Two-Phase HTI Medium Based on Reformulated BISQ Mechanism[J]. Chinese Journal of Geophysics, 2013, 56(10):3461-3473. [14] Qiao Z H. Theory and Modelling of Viscoelastic Ani-sotropic Media Using Fractional Time Derivative[C]//78th EAGE Conference & Exhibition 2016.[S. l.]:EAGE, 2016. [15] 董良国,马在田,曹景忠. 一阶弹性波动方程交错网格高阶差分解法稳定性研究[J]. 地球物理学报,2000,43(6):856-864. Dong Liangguo, Ma Zaitian, Cao Jingzhong. A Study on Stability of the Staggered-grid High-Order Difference Method of First-Order Elastic Wave Equation[J]. Chinese Journal of Geophysics,2000, 43(6):856-864. [16] Cerjan C, Kosloff D, Kosloff R, et al. A Nonref-lecting Boundary Condition for Discrete Acoustic and Elastic Wave Equation[J]. Geophysics, 1985, 50(4):705-708. |
| [1] | 黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385. |
| [2] | 岳崇旺, 王祝文. Kelvin-Voigt黏弹性介质井孔声场有限差分数值模拟与波形特征[J]. 吉林大学学报(地球科学版), 2021, 51(4): 1268-1275. |
|