隧道空间有限元地震波场模拟与逆时偏移成像

发布时间：2017-12-28 22:20

  本文关键词：隧道空间有限元地震波场模拟与逆时偏移成像　出处：《中国地质大学》2017年博士论文　论文类型：学位论文

  更多相关文章： 隧道地震超前预报 三维有限元 波场模拟 边界条件 弹性波逆时偏移

【摘要】：隧道工程作为交通建设的咽喉,其安全和高效地施工一直受到人们的重视。但隧道围岩地质情况复杂多变,施工过程中断层破碎带、岩溶等不良地质体极易造成突水突泥、围岩失稳等地质灾害。因此,开展隧道超前预报的工作,准确掌握隧道掌子面前方的地质情况,并提前做好预防工作,是隧道工程安全和高效施工的重要保障。隧道超前预报技术主要是从20世纪70年代开始兴起的,其中地震反射波法超前预报技术以其精度高,探测距离远,而备受工程界青睐。但隧道复杂的围岩地质情况使得地震波场异常复杂,难以识别有效波场和获取较可靠的成像结果。目前,针对隧道超前预报的研究主要是从波场正演模拟与成像两方面进行的,但二者仍然在以下几个方面值得进一步的研究和探讨:1)对于隧道正演模拟来说,通常采用的方法是有限差分,有限差分在处理隧道这种特殊的弧形自由边界时往往比较困难,不利于获得符合实际的隧道空间地震波场响应特征;2)对于隧道超前预报偏移成像方法的研究大多是采用Kirchhoff偏移或是绕射扫描偏移,而对目前公认的成像精度最高的逆时偏移也大多是基于二维空间,缺乏针对隧道地震超前预报三维逆时偏移成像的研究以及有利于成像的观测系统讨论;3)目前针对隧道地震超前预报中获得的多分量地震资料,主要的处理流程是基于声波偏移或是纵横波分离式的标量偏移。由于隧道地震资料的复杂性,往往很难通过对地震记录进行波场分离获得纯纵波或是纯横波资料,在偏移成像时,残留的能量容易对成像造成干扰。针对目前隧道超前预报研究中值得进一步探讨的问题,本文开展了基于有限元法的隧道三维空间弹性波波场模拟和基于弹性波逆时外推的隧道三维纵波逆时偏移成像研究。主要工作及成果如下:(1)基于二阶位移弹性波波动方程,利用虚功原理推导并建立了三维弹性波波场模拟的有限元方程,并介绍了具有二阶精度的显式Newmark差分算法的实现过程;对比分析了行压缩存储(CSR)以及一维变带宽存储对稀疏矩阵的存储效率;在给出了非结构网格下有限元爆炸震源的加载方式下,实现了三维有限元弹性波全波场模拟,为后文开展逆时偏移成像研究提供了基础。(2)基于完美匹配层理论,详细推导并建立了基于二阶位移波动方程的分裂完美匹配层(PML)边界条件和非分裂卷积完美匹配层(CFS-UPML)边界条件的有限元方程,通过与解析解对比验证了本文有限元模拟面波的正确性,对比分析了两种边界条件对面波、掠射角入射波和近地表高泊松比介质下的吸收效果,以及长时间计算的稳定性情况;最后对有限元在三维倾斜自由表面模拟面波的正确性进行了验证。研究结果表明:具有网格剖分灵活以及自动满足自由边界条件等优势的有限元法能够获得较精确的面波记录,并且适应于复杂环境下的波场模拟;CFS-UPML边界条件在对掠射角入射波的吸收效果以及高泊松比介质中的长时间稳定性方面更具优势,但计算效率要低于PML边界。(3)基于隧道施工过程中常见的不良地质体(岩性分界面、溶洞)的结构特点,构建了三维隧道地质地球物理模型,实现了隧道波场的三维弹性波全波场模拟,并进行了隧道波场特征分析,增强了对隧道波场传播特征的系统认识,为隧道地震反射信息的识别提供了理论依据。分析结果显示:炮点位于隧道边墙时,岩性分界面产生的反射横波能量相对来说要强于反射纵波;检波器布置在靠近倾斜界面一侧的隧道侧壁上时,获得的反射纵波能量强于另外一侧,更有利于反射纵波成像;溶洞充填空气时,反射波的能量最强,其次是充填水,而充填软弱土产生的反射波能量最弱。(4)在分析了弹性波逆时偏移成像基本理论的基础上,对逆时偏移成像中的关键步骤(波场分离、成像条件以及成像噪音压制)进行了详细的分析和讨论,并对其在隧道逆时偏移成像中的适用性进行了验证。结果表明:采用基于Helmholtz分解的方法能够有效的对隧道波场进行纵横波解耦,以及边界存储策略的互相关成像条件的实施方式,能够有效解决三维逆时偏移的内存消耗问题;基于Poynting矢量和Laplace滤波的方法能够有效消除隧道逆时偏移成像中的低频噪音干扰。通过上述讨论,最终建立了基于多分量地震资料-弹性波波场外推-纵横波分离-弹性纵波成像的隧道地震超前预报三维弹性纵波逆时偏移成像的流程。(5)在隧道有限元三维弹性波波场模拟以及弹性波逆时偏移成像的研究基础上,实现了基于弹性波逆时外推的隧道三维纵波逆时偏移成像。基于60°倾斜界面模型,对影响隧道成像结果的主要因素(观测系统以及噪音)进行了讨论,结果表明:基于弹性波逆时外推的三维弹性纵波成像能够有效的进行隧道成像;双边观测系统更有利于倾斜界面的成像,理论记录成像结果表明采用一条排列即可满足成像精度要求。噪声讨论结果表明:当噪音比例小于等于50%时,获得的成像剖面均较为可靠;随着噪音比例的增加,成像结果中虚假同相轴的能量不断增加,当噪音比例大于70%时,成像结果才会受到虚假成像同相轴的干扰,可见,本文采用的隧道逆时偏移成像流程具有较好的抗噪性。(6)基于上述三维弹性纵波逆时偏移成像的流程,开展了隧道典型不良地质体的理论模型(与隧道轴线成90°夹角)成像研究。结果表明:本文采用的基于多分量地震资料和弹性波波场外推的三维弹性纵波逆时偏移成像方法能够有效的对掌子面正前方不良地质体理论模型进行成像。通过与倾斜界面以及侧面溶洞的成像结果对比可知,偏移剖面中与隧道轴线成90°夹角的地质体成像能量相对较弱,因此,在高背景噪声情况下,很难实现对与隧道轴线成90°夹角地质体的精确成像。(7)基于本文讨论的观测系统以及三维弹性纵波成像方法进行了实际隧道地震超前预报资料的成像研究,通过对含有噪音记录的成像,说明了本文采用方法具有一定的抗噪性;与二维偏移成像结果的对比,验证了本文采用的成像方法的可行性和有效性。本文的创新性成果主要体现在:(1)基于二维分裂PML边界条件,将其有效的扩展到了三维空间,并实现了基于二阶位移方程的分裂PML边界条件和CFS-UPML边界条件的有限元隧道地震全波场模拟;(2)有效实现了三维隧道纵波逆时偏移成像,并对影响成像精度的关键参数进行了系统讨论与分析,建立了满足成像精度要求的观测系统,形成了基于多分量地震资料-弹性波波场外推-纵横波分离-噪声压制-弹性纵波成像的隧道地震超前预报三维弹性纵波逆时偏移成像的流程。
[Abstract]:As the throat of traffic construction, the safety and efficient construction of tunnel engineering has been paid much attention by people. But the geological condition of tunnel surrounding rock is complex and changeable. During construction process, such unfavorable geological bodies as interruption zone, fracture zone and karst are very likely to cause geological disasters such as water bursting, mud outburst and surrounding rock instability. Therefore, to carry out the tunnel prediction work, accurately grasp the geological conditions of the tunnel face, and advance the prevention work, is an important guarantee for the safe and efficient construction of tunnel engineering. The advanced prediction technology of tunnel is mainly started in 1970s. The advanced prediction technology of seismic reflection wave is favored by engineering circles because of its high accuracy and long detection distance. However, the complex geological condition of the tunnel makes the seismic wave field very complicated, and it is difficult to identify the effective wave field and obtain more reliable imaging results. At present, the research of tunnel prediction is mainly from the wave field forward modeling and imaging of the two aspects, but the two still in the following aspects worthy of further research and discussion: 1) for the simulation of tunnel, is a commonly used method of finite difference, finite difference in the tunnel special arc free boundary is often difficult to get in line with the actual seismic wave field response characteristics of tunnel space; 2) for the prediction of migration method in tunnel mostly uses the Kirchhoff offset or diffraction scan migration, and the currently accepted imaging the highest accuracy of reverse time migration was also based on two-dimensional space in view of the lack of tunnel seismic prediction, 3D reverse time migration imaging research and observation system for imaging the discussion; 3) for the tunnel seismic prediction for multi points The main processing flow of seismic data is based on the wave offset or the scalar offset of the longitudinal and transverse wave separation. Because of the complexity of tunnel seismic data, it is often difficult to get pure longitudinal wave or pure shear wave data by separating wave field from seismic records. When migration imaging, residual energy can easily cause interference to imaging. Aiming at the problems that should be further explored in the study of tunnel advanced prediction, this paper develops three-dimensional elastic wave field simulation based on finite element method and 3D longitudinal wave reverse time migration imaging based on elastic wave inverse time extrapolation. The main work and results are as follows: (1) two order displacement based on the elastic wave equation is derived by using the principle of virtual work, and established a 3D elastic wave field simulation of finite element equations, and introduces explicit Newmark with two order accuracy difference algorithm implementation process; comparative analysis of the compressed storage (CSR) and one dimension variable bandwidth of sparse matrix storage efficiency; finite element is given under unstructured explosive source loading method, the three-dimensional finite element full elastic wave field simulation, for the development of reverse time migration imaging research provides a basis for. (2) based on the theory of perfectly matched layer, is derived and the matching layer two order displacement wave equation based on perfect splitting (PML) boundary conditions and non dividing convolutional perfectly matched layer (CFS-UPML) boundary conditions of finite element equations, by comparison with analytic solutions and verify the correctness of the finite element simulation of surface waves. Comparative analysis of two kinds of boundary conditions of surface wave, the grazing angle of incident wave and near surface absorption of medium high Poisson's ratio, and the stability of long time calculation; finally the finite element simulation of free surface in three-dimensional inclined surface wave are verified. The results show that: with the finite element method and flexible mesh automatically satisfy the free boundary conditions and other advantages can obtain more accurate records of surface wave, wave field simulation and adapt to the complex environment; CFS-UPML boundary conditions on the grazing angle of the incident wave absorption effect and high Poisson's ratio of long time stability of medium the advantages, but the computational efficiency is lower than the PML boundary. (3) common adverse geological tunnel construction process based on (lithologic interface, cave) structure characteristics, build the three-dimensional tunnel geological geophysical model, the 3D elastic wave tunnel wave field full wave field simulation, and analyzed the tunnel wave field characteristics, and enhance the understanding of wave propagation in tunnel the characteristics of the system, provides a theoretical basis for the identification of tunnel seismic reflection information. The analysis results show that the shot point is located in the tunnel wall, shear wave reflection energy generated by lithologic interface relatively stronger than the reflection wave; detector disposed near the tunnel sidewall tilt interface on one side, the other side is stronger than the P-wave reflection energy, more conducive to reflection wave imaging; cave filling the air, the strongest energy reflection wave, followed by water filling, filling and reflection of soft soil produced by wave energy. (4) in the analysis of the basic theory of elastic wave reverse time migration, the key steps in the reverse time migration imaging (wave field separation, imaging conditions and imaging noise suppression) are analyzed and discussed in detail, and the applicability of the offset in the tunnel in reverse time imaging is verified. The results show that the method based on Helmholtz decomposition of effective tunnel wave field wave decoupling, implementation method and boundary storage strategy of cross-correlation imaging condition, can effectively solve the problem of offset memory consumption inverse method; Poynting vector and Laplace filter can effectively eliminate the low-frequency noise migration in tunnel based on the inverse. Based on the above discussions, a three-dimensional elastic longitudinal wave reverse time migration imaging process based on multi-component seismic data elastic wave field extrapolation longitudinal and transverse wave separation elastic longitudinal wave imaging is finally established. (5) on the basis of the study of the three-dimensional elastic wave field simulation and the inverse time migration imaging of the elastic wave in the tunnel finite element, the tunnel based on the extrapolation of the elastic wave inverse time is realized.
【学位授予单位】：中国地质大学
【学位级别】：博士
【学位授予年份】：2017
【分类号】：U452.11;P631.4

【相似文献】

中国期刊全文数据库 前10条

1 薛桂霞,邓世坤,刘秀娟;逆时偏移在探地雷达信号处理中的应用[J];煤田地质与勘探;2004年01期

2 杨勤勇;段心标;;逆时偏移技术发展现状与趋势[J];石油物探;2010年01期

3 张春燕;李振春;孙小东;;逆时偏移方法技术进展综述[J];勘探地球物理进展;2010年05期

4 滕厚华;许卫华;于晓莉;尚新民;贺红蕾;;无反射递推算法叠后逆时偏移的研究与应用[J];油气地球物理;2011年04期

5 许璐;孟小红;刘国峰;;逆时偏移去噪方法研究进展[J];地球物理学进展;2012年04期

6 王童奎;李莹;郭爱华;石文武;杨晓利;孙明;潘欢;;逆时偏移技术在南堡1号构造中的应用研究[J];地球物理学进展;2012年06期

7 康智清;刘恩;司杰戈;陈康;;逆时偏移剖面假象产生机制及其消除方法[J];石油天然气学报;2013年01期

8 胡明顺;潘冬明;董守华;金红娣;赵新;;基于散射理论的逆时偏移随机边界构建策略优化分析[J];地球物理学进展;2013年04期

9 张叔伦;陈静;张关泉;;逆时偏移的方向导数算法[J];石油地球物理勘探;1991年06期

10 Ruey-Chyuan Shih;A.R.Levander;陈炳文;;层剥离逆时偏移[J];国外油气勘探;1995年04期

中国重要会议论文全文数据库 前10条

1 杜启振;秦童;;逆时偏移最佳匹配层边界实现[A];中国地球物理·2009[C];2009年

2 陈康;;逆时偏移剖面的波数域分析及高阶拉普拉斯滤波算法的实现[A];“地球·资源”全国博士生学术论坛会议论文摘要集[C];2011年

3 龙云;韩立国;韩利;邓武斌;;单程波逆时偏移[A];中国地球物理学会第二十七届年会论文集[C];2011年

4 李庆洋;李振春;黄建平;;基于双变算法的精细逆时偏移成像[A];中国地球物理2013——第二十二专题论文集[C];2013年

5 曹晓莉;黄建平;李振春;;最小二乘逆时偏移高精度成像方法研究[A];中国地球物理2013——第二十二专题论文集[C];2013年

6 许璐;孟小红;刘国峰;;成像后压制逆时偏移噪声[A];中国地球物理学会第二十七届年会论文集[C];2011年

7 薛霆哠;张书杰;;逆时偏移中坡印廷向量滤波的改进方法[A];中国地球物理2013——第二十二专题论文集[C];2013年

8 杨仁虎;毛伟健;凌云;;逆时偏移储层成像波形校正研究[A];中国地球物理2013——第二十专题论文集[C];2013年

9 曾同生;李振春;;基于叠前逆时偏移的速度分析方法研究[A];中国地球物理学会第二十七届年会论文集[C];2011年

10 胡昊;刘伊克;王一博;;地震叠前逆时偏移不同边界处理方式研究[A];中国地球物理学会第二十七届年会论文集[C];2011年

中国重要报纸全文数据库 前2条

1 王跃;物探院逆时偏移成像技术跻身国际前列[N];中国石化报;2013年

2 本报记者 刘洪宇;石油HPC走向“平衡计算”[N];中国计算机报;2011年

中国博士学位论文全文数据库 前10条

1 王京;隧道空间有限元地震波场模拟与逆时偏移成像[D];中国地质大学;2017年

2 李静爽;基于近似解析离散化算子的逆时偏移方法及其应用[D];清华大学;2014年

3 程飞;多井井间三维波场数值模拟与逆时偏移成像研究[D];中国地质大学;2015年

4 冯兰兰;保辛几何结构不变的逆时偏移方法研究[D];清华大学;2015年

5 杨佳佳;多分量地震波逆时偏移的关键技术研究[D];中国海洋大学;2015年

6 高新成;基于云计算的逆时偏移数据处理方法研究[D];东北石油大学;2016年

7 Naveed Akram;基于弹性波模拟的逆时偏移以及地震数据解释[D];中国科学技术大学;2017年

8 王娟;TTI介质准纵波逆时偏移方法研究[D];中国石油大学(华东);2015年

9 赵旭;时域声波方程逆时偏移成像及其角道集计算[D];吉林大学;2015年

10 蔡志成;基于Q补偿的黏弹介质逆时偏移成像研究[D];中国地质大学;2017年

中国硕士学位论文全文数据库 前10条

1 陈洪杰;基于声波方程的数值模拟与逆时偏移方法研究[D];吉林大学;2010年

2 魏巍;基于叠前逆时偏移方法的圆管损伤识别[D];哈尔滨工业大学;2015年

3 张敬秋;交错网格有限差分法地震资料逆时偏移[D];成都理工大学;2015年

4 刘诗竹;VSP数据成像方法研究[D];东北石油大学;2015年

5 柯璇;复杂构造叠前逆时偏移成像方法研究[D];东北石油大学;2015年

6 侯维丽;逆时偏移成像[D];东北石油大学;2015年

7 郭雪豹;VSP地震资料逆时偏移方法研究[D];东北石油大学;2015年

8 张千祥;纯P波各向异性介质波场数值模拟及逆时偏移研究[D];吉林大学;2016年

9 于江龙;声波逆时偏移及噪声压制[D];吉林大学;2016年

10 方修政;各向异性介质逆时偏移方法研究[D];东北石油大学;2016年

，

本文编号：1347664