Numerical simulation of the wavefield in a viscous fluid-saturated two-phase VTI medium based on the constant-Q viscoelastic constitutive relation with a fractional temporal derivative
-
摘要:
分数阶微分算子具有描述历史依赖性和全域相关性的特质,本文利用这种特质描述双相介质固体骨架的黏弹性特征.基于Kjartansson常Q理论将含有分数阶时间导数的黏弹固体骨架各向异性本构关系与双相介质理论有机地结合起来,并引入流变学本构关系描述孔隙流体的黏滞性力学行为,提出一种新的基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI模型.推导了相应的时间域波传播方程,然后对该方程进行了数值模拟.对整数阶导数采用高阶交错网格有限差分算法,对分数阶时间导数采用短时记忆中心差分算法,进行了不同相界、不同品质因子组及双层地质结构情况下该类介质中波场的数值模拟与特征分析.模拟结果表明:将含有分数阶时间导数的常Q黏弹固体骨架各向异性本构关系及孔隙流体的黏滞性本构关系引入双相介质理论是可行的,二者的结合能更好地反映地下介质的黏弹性特征,对于进一步认识波在黏弹各向异性孔隙介质中的传播机理具有重要意义,为反演和重构地下油气储层和结构奠定正演理论基础.
Abstract:Fractional differential operators have the trait of describing historical dependence and global correlation.In this paper, we use this trait to describe the viscoelastic characteristics of a two-phase solid matrix.Based on the Kjartansson constant-Q model theory, the anisotropic constitutive relation of the viscoelastic solid skeleton with fractional temporal derivative is combined with the two-phase medium theory, and the rheological constitutive relation is introduced to describe the viscous mechanical behavior of the pore fluid.A new VTI model of a vicious fluid-saturated two-phase VTI medium based on the constant-Q viscoelastic constitutive relation with a fractional temporal derivative is proposed, and the corresponding wavefield propagation equation is deduced in the time domain.In numerical simulation, the finite difference method on a high-order staggered grid is used for the integer order derivative, and the short time memory center difference algorithm is employed to discretize the fractional order temporal derivative.Following the equations, the wavefield in this medium model is simulated for different phase boundaries, different Q values, and a two-layer structure, and then the wavefield features are analyzed.The results of numerical modeling indicate that it is feasible to introduce the constant-Q viscoelastic solid anisotropic constitutive relation containing a fractional temporal derivative and the viscous constitutive relation of the pore fluid into the two-phase medium theory, and the combination of both the relations can better reflect the viscoelastic characteristics of the underground medium.These results are of importance for further understanding the propagation mechanism in a viscoelastic anisotropic porosity medium, and provide a theoretical basis for inversion and reconstruction of underground oil and gas reservoirs and structures.
-
-
-
Aki K, Richards P G. 1980. Quantitative Seismology:Theory and Methods (Vol. 1 and Vol. 2). New York:Freeman.
Biot M A. 1941. General theory of three-dimensional consolidation. Journal of Applied Physics, 12(2):155-164. doi: 10.1063/1.1712886
Biot M A. 1956a. Theory of propagation of elastic waves in a fluid-saturated porous solid, Ⅰ:Low-frequency range. The Journal of the Acoustical Society of America, 28(1):168-178. http://gji.oxfordjournals.org/external-ref?access_num=10.1121/1.1908239&link_type=DOI
Biot M A. 1956b. Theory of propagation of elastic waves in a fluid-saturated porous solid, Ⅱ:Higher frequency range. The Journal of the Acoustical Society of America, 28(2):179-191. doi: 10.1121/1.1908241
Caputo M, Mainardi F. 1971. Linear models of dissipation in anelastic solids. La Rivista del Nuovo Cimento, 1(2):161-198. doi: 10.1007/BF02820620
Carcione J M. 2009. Theory and modeling of constant-Q P-and S-waves using fractional time derivatives.Geophysics, 74(1):T1-T11. doi: 10.1190/1.3008548
Carcione J M, Cavallini F, Mainardi F, et al. 2002. Time-domain modeling of constant-Q seismic waves using fractional derivatives. //Pšeník I, Červený V eds. Seismic Waves in Laterally Inhomogeneous Media. Basel: Birkhäuser.
Carcione J M, Kosloff R, Kosloff R. 1988. Wave propagation simulation in a linear viscoelastic medium. Geophysical Journal International, 95(3):597-611. doi: 10.1111/j.1365-246X.1988.tb06706.x
Cerjan C, Kosloff D, Kosloff R, et al. 1985. A nonreflecting boundary condition for discrete acoustic and elastic wave equation. Geophysics, 50(4):705-708. doi: 10.1190/1.1441945
Christensen R M. 1982. Theory of Viscoelasticity-An Introduction. 2nd ed. New York:Academic Press Inc.
Diallo M S, Appel E. 2000. Acoustic wave propagation in saturated porous media:Reformulation of the Biot/Squirt flow theory. Journal of Applied Geophysics, 44(4):313-325. doi: 10.1016/S0926-9851(00)00009-4
Dvorkin J, Nur A. 1993. Dynamic poroelasticity:A unified model with the squirt and the Biot mechanisms. Geophysics, 58(4):524-533. doi: 10.1190/1.1443435
Emmerich H, Korn M. 1987. Incorporation of attenuation into time-domain computations of seismic wave fields. Geophysics, 52(9):1252-1264. doi: 10.1190/1.1442386
Guo Z B, Tian K, Li Z C, et al. 2014. Constructing nearly constant-Q viscoelastic model using the nonlinear optimization method. Journal of China University of Petroleum (in Chinese), 38(2):52-58. http://en.cnki.com.cn/Article_en/CJFDTotal-SYDX201402008.htm
Guo Z Q. 2008. Wave field modeling in viscoelastic anisotropic media and reservoir information study [Ph. D. thesis] (in Chinese). Changchun: Jilin University.
Kjartansson E. 1979. Constant Q-wave propagation and attenuation. Journal of Geophysical Research:Solid Earth, 84(B9):4737-4748. doi: 10.1029/JB084iB09p04737
Li H X, Liu C, Tao C H. 2007. Elastic wave high-order staggering grid finite-difference numeric simulation based on transversely isotropic BISQ model. Oil Geophysical Prospecting (in Chinese), 42(6):686-693. http://en.cnki.com.cn/Article_en/CJFDTOTAL-SYDQ200706014.htm
Liu C, Lan H T, Guo Z Q, et al. 2013. Pseudo-spectral modeling and feature analysis of wave propagation in two-phase HTI medium based on reformulated BISQ mechanism. Chinese Journal of Geophysics (in Chinese), 56(10):3461-3473, doi:10.6038/cjg20131021.
Liu H P, Anderson D L, Kanamori H. 1976. Velocity dispersion due to anelasticity:Implications for seismology and mantle composition. Geophysical Journal International, 47(1):41-58. doi: 10.1111/j.1365-246X.1976.tb01261.x
Liu Q R, Katsube N J. 1990. The discovery of a second kind of rotational wave in a fluid-filled porous material. The Journal of the Acoustical Society of America, 88(2):1045-1053. doi: 10.1121/1.399853
Liu Y, Li C C. 2000. Study of elastic wave propagation in two-phase anisotropic media by numerical modeling of pseudospectral method. Acta Seismologica Sinica (in Chinese), 22(2):132-138. https://rd.springer.com/content/pdf/10.1007/s11589-000-0003-1.pdf
Lu M H, Ba J, Yang H Z. 2009. Propagation of elastic waves in a viscous fluid-saturated porous solid. Engineering Mechanics (in Chinese), 26(5):36-40. http://en.cnki.com.cn/Article_en/CJFDTotal-GCLX200905009.htm
Nie J X, Ba J, Yang D H, et al. 2012. BISQ model based on a Kelvin-Voigt viscoelastic frame in a partially saturated porous medium. Applied Geophysics (in Chinese), 9(2):213-222. doi: 10.1007/s11770-012-0332-6
Nie J X, Yang D H, Ba J. 2010. Velocity dispersion and attenuation of waves in low-porosity-permeability anisotropic viscoelastic media with clay. Chinese Journal of Geophysics (in Chinese), 53(2):385-392, doi:10.3969/j.issn.0001-5733.2010.02.016.
Nie J X, Yang D H, Yang H Z. 2004. Inversion of reservoir parameters based on the BISQ model in partially saturated porous media. Chinese Journal of Geophysics (in Chinese), 47(6):1101-1105.
Parra J O. 1997. The transversely isotropic poroelastic wave equation including the Biot and the squirt mechanisms:Theory and application. Geophysics, 62(1):309-318. doi: 10.1190/1.1444132
Parra J O. 2000. Poroelastic model to relate seismic wave attenuation and dispersion to permeability anisotropy. Geophysics, 65(1):202-210. doi: 10.1190/1.1444711
Podlubny I. 1999. Fractional Differential Equations. San Diego:Academic Press.
Qiao Z H, Sun C Y, Wu D S. 2016. Theory and modelling of viscoelastic anisotropic media using fractional time derivative. //78th EAGE Conference and Exhibition 2016. Vienna Austria: EAGE.
Sun C Y, Xiao Y F, Yin X Y, et al. 2010. Study on the stability of finite difference solution of visco-elastic wave equations. Acta Seismologica Sinica (in Chinese), 32(2):147-156. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZXB201002003.htm
Sun C Y, Yin X Y. 2007. Construction of constant-Q viscoelastic model with three parameters. Acta Seismologica Sinica (in Chinese), 29(4):348-357, 448. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxy200704002&dbname=CJFD&dbcode=CJFQ
Yang D H, Zhang Z J. 2002. Poroelastic wave equation including the Biot/squirt mechanism and the solid/fluid coupling anisotropy. Wave Motion, 35(3):223-245. doi: 10.1016/S0165-2125(01)00106-8
Yang D H, Shen Y Q, Lin E R. 2007. Wave propagation and the frequency domain Green's functions in viscoelastic Biot/squirt (BISQ) media. International Journal of Solids and Structures, 44(14-15):4784-4794. doi: 10.1016/j.ijsolstr.2006.12.001
Yang K D, Song G J, Li J S. 2011. FCT compact difference simulation of wave propagation based on the Biot and the squirt-flow coupling interaction. Chinese Journal of Geophysics (in Chinese), 54(5):1348-1357, doi:10.3969/j.issn.0001-5733.2011.05.024.
Yang K D, Yang D H, Wang S Q. 2002. Wave-field simulation based on the Biot-Squirt equation. Chinese Journal of Geophysics (in Chinese), 45(6):853-861. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWX200206013.htm
Zhu T Y, Carcione J M. Harris J M. 2013. Approximating constant-Q seismic propagation in the time domain. Geophysical Prospecting, 61(5):931-940. doi: 10.1111/gpr.2013.61.issue-5
Zhu Y P, Tsvankin I. 2006. Plane-wave propagation in attenuative transversely isotropic media. Geophysics, 71(2):T71-T30. http://www.mendeley.com/research/approximate-dispersion-relations-for-qpqsvwaves-in-transversely-isotropic-media/
郭振波, 田坤, 李振春等. 2014.利用非线性最优化方法构建近似常Q黏弹模型.中国石油大学学报(自然科学版), 38(2):52-58. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sydx201402008&dbname=CJFD&dbcode=CJFQ
郭智奇. 2008. 粘弹各向异性介质波场模拟与储层信息研究[博士论文]. 长春: 吉林大学.
李红星, 刘财, 陶春辉. 2007.基于横向各向同性BISQ模型的弹性波高阶交错网格有限差分数值模拟.石油地球物理勘探, 42(6):686-693. http://mall.cnki.net/magazine/Article/SYDQ200706014.htm
刘财, 兰慧田, 郭智奇等. 2013.基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析.地球物理学报, 56(10):3461-3473, doi:10.6038/cjg20131021. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract9786.shtml
刘洋, 李承楚. 2000.双相各向异性介质中弹性波传播伪谱法数值模拟研究.地震学报, 22(2):132-138. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200002002.htm
卢明辉, 巴晶, 杨慧珠. 2009.含粘滞流体孔隙介质中的弹性波.工程力学, 26(5):36-40. http://www.oalib.com/paper/4189247
聂建新, 巴晶, 杨顶辉等. 2012.基于Kelvin-Voigt黏弹性骨架的含非饱和流体孔隙介质BISQ模型.应用地球物理, 9(2):213-222. http://www.cnki.com.cn/Article/CJFDTotal-DQWJ201501036.htm
聂建新, 杨顶辉, 巴晶. 2010.含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究.地球物理学报, 53(2):385-392, doi:10.3969/j.issn.0001-5733.2010.02.016. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract1278.shtml
聂建新, 杨顶辉, 杨慧珠. 2004.基于非饱和多孔隙介质BISQ模型的储层参数反演.地球物理学报, 47(6):1101-1105. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract563.shtml
孙成禹, 肖云飞, 印兴耀等. 2010.黏弹介质波动方程有限差分解的稳定性研究.地震学报, 32(2):147-156. http://www.cnki.com.cn/Article/CJFDTOTAL-WTHT201206020.htm
孙成禹, 印兴耀. 2007.三参数常Q粘弹性模型构造方法研究.地震学报, 29(4):348-357, 448. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzxb200704002&dbname=CJFD&dbcode=CJFQ
杨宽德, 宋国杰, 李静爽. 2011. Biot流动和喷射流动耦合作用下波传播的FCT紧致差分模拟.地球物理学报, 54(5):1348-1357, doi:10.3969/j.issn.0001-5733.2011.05.024. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract7968.shtml
杨宽德, 杨顶辉, 王书强. 2002.基于Biot-Squirt方程的波场模拟.地球物理学报, 45(6):853-861. http://manu39.magtech.com.cn/Geophy/CN/abstract/abstract3547.shtml
-