[发明专利]高速层覆盖下基于初至拾取的近地表参数反演方法

权利要求书

1.高速层覆盖下基于初至拾取的近地表参数反演方法，其特征是包括：

步骤101：输入原始地震单炮记录，并对地震记录信噪比进行定性估计，根据信噪比的高、中、低选取合适的阈值ε，根据近地表参数的先验值及透射定律，估算出高速层底界面出射角度θ₁，从而建立初始的“信噪比S/N-阈值ε-高速层底界面出射角度θ₁”量版；

步骤102：基于能量比法引入阈值ε作为约束进行初至拾取；

步骤103：根据拾取到的初至直达波信息做出x²-t²曲线，其中x和t分别是为初至直达波的偏移距及其走时，利用其斜率求取上覆高速层速度v₁；

步骤104：利用“高速帽子”的断点信息以及特定阈值ε对应的高速层出射角度θ₁信息求取上覆高速层厚度h₁、低速带速度v₂；

步骤105：在一定范围内变化v₁，以满足射线时-距关系为准则，通过对初至信息的扫描求取到最合适的v₁,h₁,v₂组合，即为反演出的近地表参数。

2.根据权利要求1所述的高速层覆盖下基于初至拾取的近地表参数反演方法，其特征是，所述步骤102基于能量比法引入阈值ε作为约束进行初至拾取方法：对输入的单炮地震记录数据，在每一道的道头加一个时窗长度为T的零数据，之后根据如下公式在地震记录上滑动时间窗口计算对应每一时刻τ的时窗能量

M τ = Σ t = τ τ + T x 2 ( t ) ]]>

计算相邻两个时刻τ+Δτ和τ的时窗能量比值A_τ,从而得出每道地震记录随时间变化的能量比曲线，基本公式如下：

A τ = M τ + Δ τ M τ = | Σ t = τ + Δ τ τ + Δ τ + T x 2 ( t ) Σ t = τ τ + T x 2 ( t ) | ]]>

其中，A_τ为对应时间点τ的前后时窗内振幅平方和的比值，x(t)为地震记录值，Δτ为滑动时间长度，T为道头增加的零值时间窗口长度；

定义变量ξ为地震记录值A和炮点附近道初至直达波振幅平均值A₀的比值，沿A_τ变化曲线依次寻找极值，舍弃奇异点，直至找到第一个满足ξ≥ε的极值点，在该极值对应时间点的基础上向前推移半个波长取为“有效初至”，即便在该极值点以前存在其他的极值点也都认为是“无效初至”，最终实现用含有阈值ε约束的能量比法完成“高速帽子”的初至拾取。