静校正方法王颜春.ppt

工区地形立体图 工区低速带厚度分布图 一、概述 静校正是在复杂地表地区进行地震勘探的一个重要问题，它影响着地震资料的处理结果和勘探效果。在复杂地表工区，一般情况下，地震资料的信噪比都很低，它对地震资料处理的各个环节要求都较高，尤其是对静校正和速度分析的精度要求较高。 3炮原始地震记录 交互静校正处理的结果 交互静校正处理的剖面 用模拟退火算法 计算静校正量 一、引 言 获取准确的静校正，取决于准确地计算道与道之间的时间延迟。常规的延迟计算方法是由叠加模型道与所在的同一CDP的其它未叠加道做互相关，假定模型道与未叠加道之间的差别只是时移和不相干的噪声，那么，互相关函数的最大值就被当作真正时间延迟的最佳估计。 相似性较好的地震道与相关曲线 模型道与未叠加道具有较好的相似性，互相关函数仅存在一个最大的峰值，对应于8ms。这个峰值无疑是一个可靠的估值。 加入不同随机静校正量后的叠加道 是不同扰动时差的叠加模型道，子波主频25hz.，从左到右，第一道扰动时差为0ms的叠加模型道，后续依次为扰动增量为5ms的叠加模型道，可看出，当扰动小于半个波形周期，虽然波形不断变宽，但仍然保持原有的层位反射特征。而扰动大于半个波形周期时，波形的畸变逐渐变大，以至面目全非 互相关函数出现几个大小近于相等的峰值，最大的峰值对应的时间并不可靠地代表真正的时间延迟。因常规的静校正只是单一地向目标函数减小的方向搜寻，近似的峰值可能以相同的可靠程度被拾取。这就是产生“周期跳跃”的根源。 (相似性较差的地震道与相关曲线) 模拟退火是一个蒙特·卡罗优化过程，它模仿退火的物理过程。将炮点或检波点扰动后的能量，转换成概率曲线，在一定的温度控制下，按概率分布函数随机取值，它即可向目标函数增大的方向搜寻，也可向目标函数减小的方向搜寻，在温度控制下，通过迭代的方法，逐渐向目标函数的全局极值逼近，当温度达到临界点，即得到全局最优解。 1、地表一致性静校正模型 设表示炮点和检波点所对应的实际观测的地震道，表示地表一致性条件下记录的地震道，定义一个算子，它与炮点和检波点对应的静校正量和有关，它把未知的地震道延迟了。那么，对于实际观测的每一地震道，方程 定义了静校正影响前后的地震道之间的关系。像这样的方程对一条地震测线来说，有成千上万个。因为和未知，所以也未知，在无假设条件下，尽管是一个线性算子，但是它的反算子却未知，因此不能用线性反演算出，实际上它是个非线性反演问题。 ? 2、目标函数 所谓求取的全局最优解，也就是这样一组静校正值，在排除动校速度的影响之下，第一不产生周期跳跃、第二使CDP道集相位同向。这时叠加剖面能量达到最大，当在最大值之前加一负号，可以将它转换为最小问题。不同的和反映出叠加剖面能量的变化，叠加剖面的能量是炮点和检波点的多元函数，即： ? 3、静校正反演 那么，如何来反演和呢？。将在和扰动范围内可能的取值当作一个向量，在此向量中的所有元素的排列组合中，必定有一组最佳的静校正，使剖面的叠加能量达到最大。怎样得到这组最佳的静校正值呢？。若彻底搜寻，在现今的计算机运算速度条件下是不可能的，因为静校正问题的可能答案很多。假定每个炮点和检波点的静校正可以取N个值，对于M个炮点和检波点，就存在个可能的答案。 ? 一、处理流程 理论试算与实际数据实例 理论模型Ⅰ的计算 理论模型Ⅱ的计算 (a) 理论模型； (a) 理论模型； (b) 用于计算的初始模型； (b) 用于计算的初始模型； (c) 3次迭代后的计算结果 (c) 6次迭代后的计算结果 ? 有干扰时的计算结果 常规CMP叠加剖面 用近地表模拟实现静校正后的CMP叠加剖面 用近地表模型实现静校正的方法，适应性较强，只要使用适当，长短波长静态影响均可消除。广义线性反演参数修正量不依赖于未知量的初值，所以应该得到比较准确的结果。但由于广义矩阵求逆存在非唯一性问题，特征值病态分布会严重影响计算精度的提高，因此最终效果与数据的质量和参数的选择均有着很大的关系，还需要进一步地进行研究。另外，由于采用了迭代运算，计算量较大。 交互迭代折射波静校正 然而对于地表很复杂的情况，现有的静校正方法都存在一定的局限性。这里提出了一种交互迭代的折射静校正方法，它适合于解决一些地表复杂工区的静校正问