地震层析成像的理论研究及其在岩溶发育区的应用.docVIP

PAGE PAGE 1 地震层析成像的理论研究及其在岩溶发育区的应用 　　摘要：基于均匀介质中的波射线传播理论，应用直射线追踪方法进行射线走时正演计算，建立了速度模型。对模型分别进行SIRT算法和LSQR算法反演，对比分析后发现，LSQR算法在精度、效率和稳定性三个方面都能达到较好的平衡。以溶洞发育区某水电站右岸部分平硐勘探为实例，利用地震层析成像法准确探明了岩体速度分布情况和区内K25溶洞的空间分布位置，通过与其它方法结果的对比验证，表明地震层析成像法对于溶洞发育区具有分辨率高、图像直观、易于解释等优点。 　　关键词：地震层析成像；SIRT算法；LSQR算法；溶洞发育区 　　中图分类号：P631文献标识码：A文章编号2013　　地震层析成像是20世纪80年代发展起来的地球物理勘探方法，通过地震波走时和衰减的观测来反演地质体的内部构造并对其成像[1]。地震波穿透岩土介质时，其速度快慢与岩土介质的弹性模量、剪切模量和密度等有关。密度大、强度高的岩石模量大，波速高、走时短；反之，波速低，走时长。因为完整坚硬的岩体对地震波能量的吸收弱，破碎岩体和松散土体的吸收强。因而地震CT波速图像能可靠地反映出各类速度差异明显的岩土体分布界线及其中的破碎岩体的破碎程度和分布[25]。地震CT特别适用于研究各类地层、构造、岩体、土体的分布界线及其力学性状，在工程地质勘查中常被用来探查断裂带、密集节理带、含水带、空洞、风化带等不良地质体的位置、性状及力学强度等[68]。 　　1硐间层析成像原理及离散图像重建原理 　　1.1层析成像原理 　　弹性波层析成像是通过对观测到的弹性波各种震相的运动学（走时、射线路径）和动力学（波形、振幅、相位、频率）资料的分析，进而反演地下介质的结构、速度分布及其弹性参数等重要信息的一种地球物理方法[910]。在实际工作中，需要根据工程岩体平硐之间的地质情况，合理布置观测系统，采用多次单点反射和多点同时接收的扇形穿透，经过逐点激发将在被测区域形成致密的射线交叉网络（见图1）。由激发与接收时间互换原理得知，每条射线声波旅行时间将被唯一地确定。然后根据射线的疏密程度及成像精度，将施测范围划分成若干规划的成像单元[1112]。由此运用适当的反演算法即可精确地获得异常体的展布形态。 　　1.2离散图像重建原理 　　地震波层析成像属于离散图像重建技术。首先通过扇形观测系统获取钻孔或平硐的首波走时数据（ti），然后通过求解大型矩阵方程来获取两孔之间速度剖面图像，根据速度剖面图像可以直观准确地判断异常体大小分布。 　　层析成像中的反演方法可分为线性方法和非线性方法两种，非线性方法主要有：遗传算法、模拟退火法和神经网络法等。在体波层析成像中，使用线性反演方法的较多，如奇异值分解法（SVD）、共轭梯度法（CG）和最小二乘法（LSQR）等（见图2）。 　　典型的联合迭代图像重建法（SIRT）不利于同一迭代的第j条射线的修改值来计算第j+1条射线的修改值，而是把第k轮迭代中由所有射线得到的修改值保存起来，在本轮对射线迭代结束时求某种平均Δ，然后由fk+1=fk+Δ对每个象元的慢度作修改，并留作下一轮迭代使用。 　　最小二乘法（LSQR）求解的是ATAx=ATb，具有计算量小的优点，并且能很容易地利用矩阵的稀疏性简化计算，因而适合求解大型稀疏问题。 　　2数值模拟 　　本文为了比较SIRT算法和LSQR算法求解大型稀疏矩阵方程的能力，建立如图3所示模型进行速度层析成像反演。正演模型一个单元数为9×12，单元边长为30×50的速度传播区域，即横向和纵向步长分别是50m和30m，其中存在一个6m×5m的长方形高速异常区和一个3m×10m长方形低速异常区。左边设置4个激发点，右边12个检波点接收。模型背景速度为3000m/s，异常区域速度分别为5000m/s和2000m/s，迭代收敛精度值e等于01×106。其反演结果见图4。 　　从图4中分析得出，在SIRT算法中，获得的数据量较少，而未知像素较多，得到的系数矩阵通常是大型的稀疏矩阵，因而方程组都是欠定的，甚至是严重欠定的。所以采用SIRT的反演方法，波速异常区域虽然可以大致反映出来，但与原始模型有一定误差，如右图中较小的速度变大，较大的速度变小，此误差源于SIRT方法的局限性；LSQR的反演方法，具有数值稳定、能充分利用矩阵稀疏性减少计算量等特点，所以采用LSQR法取得了较好的结果，波速异常区域在图4（a）中明显反映出来，与原始模型吻合，且波速误差较小。通过两种算法模型的对比分析得出，LSQR算法在精度、效率和稳定性三个方面都能达到较好的平衡，运用LSQR算法进行层析处理能获得较为满意的图像。 　　3工程实例 　　3.1工区概况及勘探目的 　　某水电站坝