论文部分内容阅读
地震层析成像是通过对观测到的地震波各种震相的运动学(走时、射线路径)和动力学(波形、振幅、相位、频率)资料的分析,进而反演由大量射线覆盖的地下介质的结构、速度分布及其弹性参数等重要信息的一种地球物理方法。该方法最初是从医学CT理论发展而来的,随着在地学上的不断完善和发展,已经取得了丰硕的成果,因而也被越来越多的人接受和推广。本论文简述了地震层析成像的发展历史,地震层析成像的基本原理以及模型构建、射线追踪及反演方法;对比分析了阻尼最小二乘法和遗传算法两种不同的反演算法的理论异同,并将以上两种反演算法应用于祁连山中东段地区。利用该区近十年的S ? P数字地震记录数据,联合使用Pg、S g、Pm、S m、Pn、S n等多种震相,基于相同的初始速度模型,进行了震源参数与速度结构的联合反演,获得了不同深度的速度平面分布图像、不同位置的地震深度与速度剖面图像以及莫霍面的深度分布特征。对比两种反演算法理论以及结果,结合地质构造资料和前人研究结果,可以得到如下结论:1、两种不同反演算法具有理论差异:阻尼最小二乘法属于线性解法,阻尼因子λ2既起着促进平衡的作用,同时也影响着收敛速度的快慢,可直接给出解估计的分辨矩阵和协方差矩阵;遗传算法是全局优化方法,采用随机转移规则进行最优解搜索,适合于各种类型的优化问题,但无法直接求解模型参数的分辨率、后验方差等信息,只有用线性化的办法获得协方差矩阵。2、从两种算法的反演结果看:阻尼最小二乘法所得结果的走时残差要大于遗传算法所得结果,另外其计算效率较低;对于边缘地震参数定位较差;所得速度结构高低速异常较大;整体上显示两种反演算法所得震中分布相对更加集中,速度平面分布及莫霍面深度分布具有相似的特点,速度剖面起伏较为一致。3、研究区内断裂多从高低速过渡带穿过,呈北西西向走向,地震多发生在断裂及其边缘,震源深度优势分布于上、中地壳;通过分析,认为本文研究的结果与地质构造资料及相关区域的其它地震层析成像结果在许多方面可得到印证。