层析速度反演静校正.pdfVIP

层析速度反演静校正 李家康 余云英 梁 宏 齐 莉 (石油勘撵开发科学研究院地球物理所) 摘 要 复杂地袭条件区的近地袭速度变化剧烈，严重影响着地震成泉质量，搞清近地表速度结构 是提高地震处理质量的首要任务。本戈采用程函方程有限差舟击计算通过速度模型的初至旅 行时射线路径，鼎后使用约束最小二秉正麦分解法CLsOR进行走时反演．从而获得近地表速度 结构。这种层析反演的方法以块体为基本单元，可以模拟和反演复杂速度结构，它不曲识别波 形，计算速度快。本方法应用于高陡山区实际散据．获得的浅近地袁速度横向变化很大。利用 近地表速度进行静校正所得的成泉质量比其它方法清晰。 引 言 近地表区的速度不均匀极大地影响着深层成象与低幅度构造形态。常规的野外静校正 和剩余静校正不能在复杂表层结构区取得满意效果，而折射静校正在速度变化剧烈地带因 模型建立困难也难取得好的效果。在复杂地表条件和低信噪比地区，反射信号往往不清楚， 折射渡初至的信噪比是最高的。折射波走时包含着地下浅层的速度信息，是求取表层速度 的有效途径。本文采用程函方程初至旅行时层析反演法求取表层速度分布。它可以处理直 达波、折射波、透射波。有了明确的速度分布，就可在此基础上进行依据速度场的静校正、波 动方程基准面校正和叠前深度偏移。 一、基本原理 地震渡旅行时与速度的关系由程函方程描述 (爱)2+(笺)2=(南)2 ㈨ 根据程函方程，利用速度模型可以计算出时间场分布，由此可以进一步追踪出射线路径P。 依据走时方程 t=Ij女d1 (2) Jlv、jz， 进行层析反演即可以得到地下的速度分布。网格剖分可有多种形式，最简便的方式是采用 矩形网格，整个区域的速度分布函数由网格结点处的速度确定。最便于处理的情形是每个 单元块内速度为常数。这种方式不仅方便而且十分适于复杂速度结构的模拟与反演。 变化的复杂情况。 对走时方程离散化得到线性方程组 I：P．s (3) 可以求出速度场(”=l／s)的分布。但走时l依赖于速度-(#。：)和射线路径P，现在的办法 是固定p求解口。如果p不正确．则解得的”也会不正确，尽管 乱：一I笃gdl (4) JI％、，， 用初始模型的射线路径代替真实地层中的射线路径在时差计算上所引起的误差仅仅是二阶 小量，但并不保证时差的二阶小量误差不会在解中被放大成一个很大的误差。解决的根本 办法是采用约束。约柬条件包括：已知的先验信息、对解的光滑度要求、正则化条件。这3 类条件统一于约束方程 C．∞：h (5) 之中，与走时方程联台求解，并利用迭代法， △‘=G．△j (6) 不断修改慢度场直至满意为止．最终可得近地表区的速度结构。 这里所涉及的方程组是个特大型非良态非线性方程组，阶数达数千乃至数万阶且射线 分布不均匀，常规求解方法通常失效。为此，我们发展了一种稳定且具有高效率的约束最小 二乘cJ．sQR算法。 层析反演整个过程分7步进行：①人工或自动交互拾取初至；⑦由初至时间给出初始速 度模型；③根据初始模型速度计算走时与射线路径；④由先验信息构制约束方程；⑤计算观 测走时