ISSN 1004-4140
CN 11-3017/P

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

频率-空间域CPML和特征分析法组合边界研究

李桂花 王树林 张文波 徐子煜

李桂花, 王树林, 张文波, 等. 频率-空间域CPML和特征分析法组合边界研究[J]. CT理论与应用研究, 2022, 31(3): 269-279. DOI: 10.15953/j.ctta.2021.075
引用本文: 李桂花, 王树林, 张文波, 等. 频率-空间域CPML和特征分析法组合边界研究[J]. CT理论与应用研究, 2022, 31(3): 269-279. DOI: 10.15953/j.ctta.2021.075
LI G H, WANG S L, ZHANG W B, et al. Combined boundary of CPML and feature analysis within frequency-space domain[J]. CT Theory and Applications, 2022, 31(3): 269-279. DOI: 10.15953/j.ctta.2021.075. (in Chinese)
Citation: LI G H, WANG S L, ZHANG W B, et al. Combined boundary of CPML and feature analysis within frequency-space domain[J]. CT Theory and Applications, 2022, 31(3): 269-279. DOI: 10.15953/j.ctta.2021.075. (in Chinese)

频率-空间域CPML和特征分析法组合边界研究

doi: 10.15953/j.ctta.2021.075
基金项目: 山东省自然科学基金面上项目(基于混合模型深度神经网络的多波地震油气藏特征提取与识别(ZR202103050722))
详细信息
    作者简介:

    李桂花:女,山东科技大学地球科学与工程学院副教授、硕士生导师,主要从事地震勘探相关理论与应用教学和研究工作,E-mail:skd993914@sdust.edu.cn

  • 中图分类号: O  242; P  631.4

Combined Boundary of CPML and Feature Analysis within Frequency-Space Domain

  • 摘要: 在地震波场数值模拟过程中,边界反射是影响其模拟结果的一个重要因素。实际地下介质具有各向异性特征,传统的完全匹配层边界(PML)对于小入射角地震波具有良好效果,但该方法并不能有效地吸收低频波和大角度入射波。针对VTI介质边界反射的问题,本文提出在频率-空间域有限差分法数值模拟中采用卷积完全匹配层(CPML)和特征分析法的组合边界条件,并对该组合边界条件进行数值模拟实验和边界反射吸收效果分析,验证所提方法是一种可靠的人工吸收边界条件,能够有效地压制波场模拟过程中产生的边界反射。

     

  • 图  1  CPML边界示意图

    Figure  1.  Boundary diagram of CPML

    图  2  未加边界时的50 Hz单频波快照(a)和时间域150 ms波前快照(b)

    Figure  2.  50 Hz single frequency snapshot (a) and 150 ms wavefront snapshot (b) without boundary attenuation in time domain

    图  3  施加特征分析法边界衰减后的50 Hz单频波快照(a)和150 ms时间域波前快照(b)

    Figure  3.  50 Hz single frequency snapshot (a) and 150 ms wavefront snapshot (b) with the boundary attenuation of characteristic analysis method in time domain

    图  4  施加PML和特征分析组合边界后的50 Hz单频波快照(a)和150 ms时间域波前快照(b)

    Figure  4.  50 Hz single frequency snapshot (a) and 150 ms wavefront snapshot (b) with the combination absorption boundary of PML and characteristic analysis method in time domain

    图  5  施加CPML和特征分析的组合边界后的50 Hz单频波快照(a)和150 ms时间域波前快照(b)

    Figure  5.  50 Hz single frequency snapshot (a) and 150 ms wavefront snapshot (b) with the combination absorption boundary of CPML and characteristic analysis method in time domain

  • [1] 戚艳平. 三维VTI介质qP波正演方法研究[D]. 东营: 中国石油大学(华东), 2008.

    QI Y P. The Study on qP wave forward modeling methods in 3D VTI media[D]. Dongying: China University of Petroleum (East China), 2008. (in Chinese).
    [2] BERENGER J P. A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 1994, 114(2): 185−200. doi: 10.1006/jcph.1994.1159
    [3] BERENGER J P. Application of the CFS-PML to the absorption of evanescent waves in waveguides[J]. IEEE Microwave and Wireless Components Letters, 2002, 12(6): 218−220. doi: 10.1109/LMWC.2002.1010000
    [4] BERENGER J P. Numerical reflection from FDTD-PMLs: A comparison of the split PML with the unsplit and CFS PMLs[J]. IEEE Transactions on Antennas and Propagation, 2002, 50(3): 258−265. doi: 10.1109/8.999615
    [5] 罗玉钦, 刘财. 近似完全匹配层边界条件吸收效果分析及衰减函数的改进[J]. 石油地球物理勘探, 2018,53(5): 903−913,879. DOI: 10.13810/j.cnki.issn.1000-7210.2018.05.003.

    LUO Y Q, LIU C. Absorption effects in nearly perfectly matched layers and damping factor improvement[J]. Oil Geophysical Prospecting, 2018, 53(5): 903−913,879. DOI: 10.13810/j.cnki.issn.1000-7210.2018.05.003. (in Chinese).
    [6] 李青阳, 吴国忱, 梁展源. 基于PML边界的时间高阶伪谱法弹性波场模拟[J]. 地球物理学进展, 2018,33(1): 228−235.

    LI Q Y, WU G C, LIANG Z Y. Time domain high-order pseudo spectral method based on PML boundary for elastic wavefield simulation[J]. Progress in Geophysics, 2018, 33(1): 228−235. (in Chinese).
    [7] 谢志南, 郑永路, 章旭斌, 等. 弱形式时域完美匹配层—滞弹性近场波动数值模拟[J]. 地球物理学报, 2019,62(8): 3140−3154.

    XIE Z N, ZHENG Y L, ZHANG X B, et al. Weak-form time-domain perfectly matched layer for numerical simulation of viscoelastic wave propagation in infinite-domain[J]. Chinese Journal of Geophysics, 2019, 62(8): 3140−3154. (in Chinese).
    [8] 张衡, 刘洪, 李博, 等. VTI介质声波方程非分裂式PML吸收边界条件研究[J]. 石油物探, 2016,55(6): 781−792.

    ZHANG H, LIU H, LI B, et al. The research on unsplit PML absorbing boundary conditions of acoustic equation for VTl media[J]. Geophysical Prospecting for Petroleum, 2016, 55(6): 781−792. (in Chinese).
    [9] 马锐, 邹志辉, 芮拥军,等. 基于SPML和海绵边界的伪谱法弹性波模拟复合吸收边界条件[J]. 石油物探, 2018,57(1): 94−103.

    MA R, ZOU Z H, RUI Y J, et al. A composite absorbing boundary based on the SPML and sponge absorbing boundary for pseudo-spectral elastic wave modeling[J]. Geophysical Prospecting for Petroleum, 2018, 57(1): 94−103. (in Chinese).
    [10] CHEW W C, LIU Q H. Perfectly matched layers for elastodynamics: A new absorbing boundary condition[J]. Journal of Computational Acoustics, 1996, 4(4): 341−359. doi: 10.1142/S0218396X96000118
    [11] WANG T L, TANG X M. Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach[J]. Geophysics, 2003, 68(5): 1749−1755. doi: 10.1190/1.1620648
    [12] 方修政, 钮凤林. 二阶声波方程非分裂式CPML实施新方法[J]. 中国科学: 地球科学, 2021,51(8): 1341−1354. doi: 10.1360/SSTe-2021-0012

    FANG X Z, NIU F L. A new implementation method of non splitting CPML for second-order acoustic equation[J]. Scientia Sinica Terrae, 2021, 51(8): 1341−1354. (in Chinese). doi: 10.1360/SSTe-2021-0012
    [13] KUZUOGLU M, MITTRA R. Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers[J]. IEEE Microwave and Wireless Components Letters, 1996, 6(12): 447−449.
    [14] RODEN J A, GEDNEY S D. Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media[J]. Microwave and Optical Technology Letters, 2000, 27(5): 334−339. doi: 10.1002/1098-2760(20001205)27:5<334::AID-MOP14>3.0.CO;2-A
    [15] DROSSAERT F H, GIANNOPOULOS A. A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves[J]. Geophysics, 2007, 72(2): T9−T17. doi: 10.1190/1.2424888
    [16] DROSSAERT F H, GIANNOPOULOS A. Complex frequency shifted convolution PML for FDTD modelling of elastic waves[J]. Wave Motion, 2007, 44(7): 593−604.
    [17] CHEN H M, ZHOU H, LI Y Q. Application of unsplit convolutional perfectly matched layer for scalar arbitrarily wide-angle wave equation[J]. Geophysics, 2014, 79(6): T313−T321. doi: 10.1190/geo2014-0103.1
    [18] MA X, YANG D H, HUANG X Y, et al. Nonsplit complex-frequency shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 1: Method[J]. Geophysics, 2018, 83(6): T301−T311. doi: 10.1190/geo2017-0603.1
    [19] MA X, YANG D H, HE X J, et al. Nonsplit complex-frequency-shifted perfectly matched layer combined with symplectic methods for solving second-order seismic wave equations — Part 2: Wavefield simulations[J]. Geophysics, 2019: T167−T179.
    [20] MA X, LI Y J, SONG J X. A stable auxiliary differential equation perfectly matched layer condition combined with low-dispersive symplectic methods for solving second-order elastic wave equations[J]. Geophysics, 2019: T193−T206.
    [21] 吴国忱, 梁锴. VTI介质准P波频率空间域组合边界条件研究[J]. 石油物探, 2005,44(4): 301−307, 15.

    WU G C, LIANG K. Combined boundary conditions of quasi-P wave within frequency-space domin in VTI media[J]. Geophysical Prospecting for Petroleum, 2005, 44(4): 301−307, 15. (in Chinese).
    [22] 吴国忱. 各向异性介质地震波传播和成像[M]. 东营: 中国石油大学出版社, 2006.

    WU G C. Seismic wave propagation and imaging in anisotropic media[M]. Dongying: China University of Petroleum Press, 2006. (in Chinese).
    [23] 董良国. 地震波数值模拟与反演中的几个关键问题研究[D]. 上海: 同济大学, 2003.

    DONG L G. Research on several key problems in seismic wave numerical simulation and inversion[D]. Shanghai: Tongji University, 2003. (in Chinese).
    [24] 李桂花, 林年添, 杨思通, 等. 井间地震高分辨率数值模拟方法及波场特征研究[M]. 徐州: 中国矿业大学出版社, 2017.

    LI G H, LIN N T, YANG S T, et al. Study on high resolution numerical simulation method and wave field characteristics of crosswell seismic[M]. Xuzhou: China University of mining and Technology Press, 2017. (in Chinese).
    [25] 张奎涛, 顾汉明, 刘少勇, 等. 基于CPML-RML组合边界条件粘弹TTI介质旋转交错网格有限差分正演模拟[J]. 石油物探, 2019,58(2): 187−198, 218.

    ZHANG K T, GU H M, LIU S Y, et al. Rotated staggered grid finite difference forward modeling for wave propagation in viscoelastic TTI media based on CPML-RML combined boundary condition[J]. Geophysical Prospecting for Petroleum, 2019, 58(2): 187−198, 218. (in Chinese).
  • 加载中
图(5)
计量
  • 文章访问数:  173
  • HTML全文浏览量:  76
  • PDF下载量:  43
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-12-21
  • 录用日期:  2022-01-17
  • 网络出版日期:  2022-01-28
  • 刊出日期:  2022-05-23

目录

    /

    返回文章
    返回