[发明专利]一种适用于起伏地表的深度域层速度修正方法

摘要

本发明提供了一种适用于起伏地表的深度域层速度修正方法，属于石油地球物理勘探的地震资料处理领域。所述方法首先对偏移距分段，将输入道集按照偏移距进行分组，并进行起伏地表叠前深度偏移，产生共成像点道集；然后进行起伏地表剩余速度分析，计算高阶剩余曲率谱，并根据高阶剩余曲率谱的局部最大值拾取深度‑剩余曲率对，然后求取速度模型修正量；再对速度模型修正量进行空间插值和平滑，并将速度模型修正量跟初始速度模型相加得到更新后的起伏地表深度域层速度模型；最后利用更新后的起伏地表深度域层速度模型代替初始速度模型，重新进行起伏地表叠前深度偏移和起伏地表剩余速度分析，从而不断提高速度模型的精度。

主权项

一种适用于起伏地表的深度域层速度修正方法，其特征在于：所述方法包括：(1)将偏移距分段，并将输入道集按照偏移距进行分组；(2)利用初始速度模型对各段偏移距道集进行起伏地表叠前深度偏移，产生共成像点道集；产生的共成像点道集按照偏移距排列，每个共成像点道集的道数与步骤(1)中的偏移距分段数相同；(3)判断共成像点道集是否拉平，如果已经拉平，则转入步骤(7)，如果没有拉平，则转入步骤(4)；(4)计算起伏地表剩余曲率谱：首先读入每个地表CDP点处偏移产生的共成像点道集，然后利用高阶剩余曲率谱计算公式计算起伏地表剩余曲率谱，并根据局部最大能量值拾取深度‑剩余曲率对；(5)利用剩余速度求取公式对每一个共成像点道集的深度‑剩余曲率对计算深度域层速度修正量；(6)对步骤(5)获得的所有深度域层速度修正量进行空间插值和平滑处理，获得整个空间的速度模型修正量，并将该速度模型修正量与初始速度模型相加，得到更新后的速度模型；转入步骤(2)，同时用更新后的速度模型代替步骤(2)中的初始速度模型；(7)输出偏移剖面；所述高阶剩余曲率谱计算公式为：$Spe(z_m,\mathrm{\γ})=\frac{{\[\underset{i=1}{\overset{Ntr}{\Σ}}U(z_{mi},i)\]}^{2n}}{\underset{i=1}{\overset{Ntr}{\Σ}}{U}^{2n}(z_{mi},i)};$所述起伏地表剩余速度分析公式为：$z_{mi}=({\mathrm{dz}}_{1i}+{\mathrm{dz}}_{2i})\mathrm{/2+}\sqrt{\frac{L_{offseti}^2}{4}({\mathrm{\γ}}^2-1)+\frac{{(2z_i-{\mathrm{dz}}_{1i}-{\mathrm{dz}}_{2i})}^2}{4}{\mathrm{\γ}}^2};$其中，z_m是深度域层速度为v_m时对应的偏移深度，γ＝v_m/v，表示叠前深度偏移所用的深度域层速度与准确层速度之间的比值，即剩余曲率值，U(z_mi，i)表示当前γ值对应的剩余曲率扫描曲线所经过共成像点道集网格(z_mi，i)处的振幅值，i是共成像点道集内的道数，代表了偏移距信息，z_mi是当前γ值对应的剩余曲率扫描曲线在第i道的偏移深度，由起伏地表剩余速度分析公式求得，Ntr表示一个共成像点道集的总道数，n为大于1的自然数；dz_1i表示地表CDP点与炮点的高程差，dz_2i表示地表CDP点与检波点的高程差，z_i表示地表CDP点与地下反射点之间的距离，L_offseti是炮点和检波点之间的偏移距。