基于高阶交错网格有限差分的隧道超前探测地震波场模拟

张焕钧; 陈祖斌; 李昊; 杨兴林

doi:10.3973/j.issn.2096-4498.2021.07.011

隧道建设 >

2021 , Vol. 41 >Issue 7: 1172 - 1179

DOI: https://doi.org/10.3973/j.issn.2096-4498.2021.07.011

研究与探索

基于高阶交错网格有限差分的隧道超前探测地震波场模拟

张焕钧 ,
陈祖斌 ,
李昊 ,
杨兴林

展开

（吉林大学仪器科学与电气工程学院, 吉林长春 130026）

张焕钧（1994—），男，新疆乌鲁木齐人，吉林大学仪器仪表工程专业在读硕士，研究方向为地震勘探仪器及数据处理。E-mail： 949987640@qq.com。 *通信作者：陈祖斌， E-mail： czb@jlu.edu.cn。

网络出版日期: 2021-07-29

基金资助

吉林省科技厅重点项目（20160204065GX）

收起

Seismic Wavefield Simulation of Tunnel Advance Detection Based on High Order Staggered Grid Finite Difference Method

ZHANG Huanjun ,
CHEN Zubin ,
LI Hao ,
YANG Xinglin

Expand

（College of Instrumentation & Electrical Engineering, Jilin University, Changchun 130026, Jilin, China)

Online published: 2021-07-29

Fold

摘要

为解决地面勘探中直接采用数值模拟方法产生的运算速度慢、边界吸收效果差、频散严重等一系列问题，提出一种针对隧道超前探测的正演模拟方法，包括：利用高阶交错网格有限差分求解一阶速度-应力弹性波方程来获取任意时刻的波场信息，选择合适的网格精度和差分阶数满足稳定条件，使用一种新衰减函数的PML吸收边界来处理人为边界造成的反射问题，利用交错网格和FCT通量传输矫正法抑制频散。通过建立非均匀介质倾斜层状隧道模型和空洞隧道模型，对上述方法进行仿真模拟，得到清晰的波场和道集数据，证明上述正演方法的可行性。

关键词： 隧道超前探测; 地震波场模拟; 数值计算; PML吸收边界

本文引用格式

张焕钧 , 陈祖斌 , 李昊 , 杨兴林 . 基于高阶交错网格有限差分的隧道超前探测地震波场模拟[J]. 隧道建设, 2021 , 41(7) : 1172 -1179 . DOI: 10.3973/j.issn.2096-4498.2021.07.011

Abstract

The use of numerical simulation methods directly in surface exploration will result in numerous issues, including slow calculation speed, poor boundary absorption effect, and significant dispersion. As a result, a forward simulation method for detecting tunnel in advance is proposed. The firstorder velocitystress elastic wave equation is calculated using a highorder staggeredgrid finite difference method to obtain wavefield information at any time; the appropriate mesh accuracy and difference order are chosen to satisfy the stability conditions, and the perfectly matched layer of a new attenuation function is used to absorb boundary to deal with the reflection problem caused by artificial boundary. Further, the staggered grid and flux corrected transport method is adopted to restrain the dispersion. Finally, in heterogeneous media, an inclined layered tunnel model and a cavity tunnel model are established to simulate the processes and obtain clear wavefields and data, proving the feasibility of the forward modeling method.

Key words： tunnel advance detection; seismic wavefield simulation; numerical calculation; perfectly matched layer absorption boundary

Options

文章导航

〈

〉

模态框（Modal）标题

摘要

本文引用格式

Abstract