静校正方法.ppt

工区地形立体图 工区低速带厚度分布图 一、概述 静校正是在复杂地表地区进行地震勘探的一个重要问题，它影响着地震资料的处理结果和勘探效果。在复杂地表工区，一般情况下，地震资料的信噪比都很低，它对地震资料处理的各个环节要求都较高，尤其是对静校正和速度分析的精度要求较高。 未做折射静校的叠加剖面 折射静校后的叠加剖面 折射静校和剩余静校及二次速度分析后的叠加剖面 折射静校和二次剩余静校及二次速度分析后的叠加剖面 3炮原始地震记录 交互静校正处理的结果 交互静校正处理的剖面 用模拟退火算法 计算静校正量 一、引 言 获取准确的静校正，取决于准确地计算道与道之间的时间延迟。常规的延迟计算方法是由叠加模型道与所在的同一CDP的其它未叠加道做互相关，假定模型道与未叠加道之间的差别只是时移和不相干的噪声，那么，互相关函数的最大值就被当作真正时间延迟的最佳估计。 初至波的旅行时间C是近地表模型参数pi(i=1,2,…，n)的函数。我们把拾取的值C0与计算的值Cc进行比较 △C=C0-Cc` 把△C按泰勒级数展开 忽略高次项，写成矩阵形式为： 可以简写成 A.ΔP=ΔC 式中，A是(m.n)矩阵，可称为灵敏度矩阵；△P是n维的列向量，可称为模型参数向量；△C是m维的列向量，可称为观测误差向量。如果我们把△P视为近地表模型参数（深度、速度）的初始值的修正量，△C为计算出的初至时间与观测的初至时间之差，则上式就是确定地表模型的线性方程组。根据线性代数的理论，矩阵Amn可以分解成， Amn=UmmDmnVtnn 式中，U和V分别是（m·m）、(n·n)的正交矩阵；D是一个由奇异值组成的对角矩阵。令矩阵A的广义逆为 A+=V·D+UT ? 则地表模型参数的修正量矩阵为 △P=A+△C 为了得到准确的近地表模型，需要进行迭代运算，迭代过程直到满足收敛条件（小于给定的误差范围或大于指定的迭代次数）为止。 ? 一、处理流程 理论试算与实际数据实例 理论模型Ⅰ的计算 理论模型Ⅱ的计算 (a) 理论模型； (a) 理论模型； (b) 用于计算的初始模型； (b) 用于计算的初始模型； (c) 3次迭代后的计算结果 (c) 6次迭代后的计算结果 ? 有干扰时的计算结果 常规CMP叠加剖面 用近地表模拟实现静校正后的CMP叠加剖面 用近地表模型实现静校正的方法，适应性较强，只要使用适当，长短波长静态影响均可消除。广义线性反演参数修正量不依赖于未知量的初值，所以应该得到比较准确的结果。但由于广义矩阵求逆存在非唯一性问题，特征值病态分布会严重影响计算精度的提高，因此最终效果与数据的质量和参数的选择均有着很大的关系，还需要进一步地进行研究。另外，由于采用了迭代运算，计算量较大。 交互迭代折射波静校正 然而对于地表很复杂的情况，现有的静校正方法都存在一定的局限性。这里提出了一种交互迭代的折射静校正方法，它适合于解决一些地表复杂工区的静校正问题。 静校正是陆地地震资料常规处理流程中必不可少的一环。在我国西北地区，地表条件比较复杂，静校正问题尤为严重。目前地震勘探的重点主要在我国的西部，在这些地区，静校正问题严重制约着地震勘探的效果，解决好静校正问题具有重要的理论意义和实际意义。 我们在推导反射波时距曲线方程时，假设观测面是一个水平面，地下传播介质是均匀的。但实际情况并非如此，观测面不是一个水平面，通常是起伏不平的，地下传播介质通常也不是均匀的，其表层还存在着低降速带的横向变化。因此野外观测得到的反射波到达时间，不满足教科书中给出的双曲线方程，而是一条畸变了的双曲线。静校正就是研究由于地形起伏、地表低降速带横向变化对地震波传播时间的影响，并进行校正。著名地球物理学家迪克斯教授生前曾说