基于TSP观测模型二维正演模拟.docVIP

基于TSP观测模型二维正演模拟 　　摘要：在TSP（隧道超前预报）中，对于前方地质体的探测是工程物探以及隧道工程中的重要的问题。但是在基于隧道模型进行全波场模拟时，隧道的复杂性又使得地震记录变得相对复杂而难以识别。而如果没有对前方地质情况有一定的了解，施工经常就会造成重大的损失。因此对本文利用交错网格有限差分法对基于隧道的模型进行各种不良地质体的模拟，对TSP常出现的地质体进行正演模拟并分析是有必要的。 　　Abstract： While constructing tunnel， its important to know what geological body in front of the tunnel. But seismic of TSP are complex usually because of the TSP observation system. Problems often occurred while we do not know the situation in front of tunnel. So we use staggered grid Finite difference method （FDM） to model varies of geological bodies based on TSP observation system. 　　关键词： TSP；全波场模拟；弹性介质；有限差分；正演模拟 　　Key words： TSP；full wave field simulation；elastic medium；FDM；forward modeling 　　中图分类号：U45 文献标识码：A 文章编号：1006-4311（2015）01-0108-04 　　0 引言 　　隧道属于地下建筑工程，其围岩地质情况复杂多变，可能会遇到断层破碎带、岩溶、松散地层、软土地层等不良地质体受经济条件及技术水平的限制，使得在勘察设计阶段，仅依据地表地质勘察工作，要完全查明上述不良地质体的状态、特性，准确预报隧道开挖过程中可能出现的不良地质体（带）的位置、规模是十分困难的。这给隧道安全施工带来较大的盲目性，因而会出现预料不到的塌方、冒顶、涌水、涌沙等事故。为了保证隧道高效施工，减少经济损失，如何准确地了解掌子面前方不良地质情况成为迫切需要解决的问题。隧道施工地质超前探测技术正是在这种情况下提出的。 　　TSP超前地质预报系统是专门为隧道与地下工程超前地质预报研制开发的设备，其基本原理是在隧道壁上不止一系列爆破点，依次进行激发并接收，产生的地震波信号在隧道周围岩体内传播，当岩石的岩性特征发生变化时，信号的一部分就会被返回并被接经过特殊处理的三分量接收器放大并接收到。经过数据处理软件分析处理以后便可以确定岩体强度变化界面的位置和方位。 　　可见TSP在隧道施工中的重要性，因此对其进行正演，研究地震波在TSP观测系统下的传播规律是有比较的。本文在假定在已知掌子面前方地质模型结构和相应物理参数的情况下，进行了2D地震波数值模拟，来研究地震波在掌子面前方的传播规律，并计算出各观测点所应观测到的数值记录。正演模拟的方法有很多种，考虑到地质状况的复杂性，本文采用有限差分法进行模拟[1]。 　　1 粘弹性介质的有限差分原理 　　有限差分法以其简单，容易实现，运算速度快[2]等优点成为了地震正演模拟的一种主要方法。下面给出二维情况速度应力方程组[3，4]： 　　=B（ + ） 　　=B（ + ） 　　=（λ+2μ） +λ 　　=（λ+2μ） +λ 　　=μ +μ 　　上述方程组中，其中B是疏度，vx（z，x，y）、vz（z，x，y）分别表示x方向和z方向的速度分量。σxx（z，x，y）、σzz（z，x，y）、τxz（z，x，y）分别表示应力分量。λ和μ分别表示拉梅系数。对上述方程进行空间离散化，便可以得到有限差分的公式。这里不详细叙述。正演模型中的边界会对地震记录造成强烈干扰，因此对于模型的边界则采用PML[5]（完美匹配层边界）对地震波能量进行吸收。 　　2 TSP理论计算模型 　　利用编制的软件，建立一个500×200的平面网格模型，网格步长为0.5m。其中设置观测模型如下图，偏移距为15m，炮间距为1.5m建立观测模型。详细观测模型参数见图1。 　　图1 （a）中所示的为正演模型，而不是实际的模型，实际模型会将检波点放在途中炮点处，并且分别在图中检波点分别激发。正演模拟的时候进行了炮检互换，由此可以在不影响结果的情况下提升模拟的效率；（b）则为观测模型的参数。 　　在TSP中，断层与溶洞是经常遇见的两种模型，研究这些地质体地震波响应特征和分析地震波的传播规律，对实际工作是有意义的。