二维叠加速度及均方根速度场闭合差校正方法

摘要

本发明涉及一种二维叠加速度及均方根速度场的两步法校正方法；首先求取T0曲面与T0闭合校正前的每条测线的时间差ΔT，利用ΔT对叠加速度做基准面校正，ΔT＞0时，校正速度用地表充填速度，ΔT＜0时，校正速度用速度谱上的速度；针对在叠加速度转换为均方根速度时由视倾角引起的均方根速度闭合差，首先获取主测线、联络测线叠加剖面的时间斜率，求取不同测线方向的视倾角和均方根速度，计算交点处的加权系数，对交点处的加权系数插值并归一化处理，利用不同测线方向的加权系数对主测线和联络线均方根速度加权叠加得到最终的均方根速度场；本方法科学有效地消除了交叉测线速度闭合差，提高了二维地震工区速度场精度。

主权项

1.一种二维叠加速度及均方根速度场校正方法，其特征在于：(1)、T₀闭合差引起的叠加速度闭合差校正①T₀闭合差的求取首先对整个二维工区的地震解释数据做T₀闭合差校正，然后网格插值，形成整个工区的地震解释等T₀图，用T₀图的曲面文件减去每条地震测线原始的没有做T₀闭合差校正的T₀散点数据，得到沿地震测线分布的ΔT散点数据，当T₀散点数据相对于T₀图上抬时ΔT为负，反之ΔT为正。②T₀闭合差引起的叠加速度闭合差校正用ΔT作为叠加速度谱的校正量，根据ΔT符号的正负，用不同的速度及不同的校正公式对速度谱做校正。当ΔT＞0时，校正速度用地表充填速度v₀第一步：用Dix公式由速度谱上的叠加速度计算第n层的层速度：$v_{\mathrm{cn}}=\sqrt{\frac{v_n^2{t}_n-v_{n-1}^2{t}_{n-1}}{t_n-t_{n-1}}}$第二步：由层速度反算叠加速度：$v_n=\sqrt{\frac{v_0^2\mathrm{\Delta t}+\underset{i=1}{\overset{n}{\Sigma }}{v}_{\mathrm{ci}}^2{\mathrm{\Delta t}}_i}{\mathrm{\Delta t}+\underset{i=1}{\overset{n}{\Sigma }}{\mathrm{\Delta t}}_i}}$当ΔT＜0时，校正速度用速度谱上的速度；第一步：用Dix公式由速度谱上的叠加速度计算第n层的层速度。第二步：由下面公式由层速度反算叠加速度$v_n=\sqrt{\frac{\underset{i=k}{\overset{n}{\Sigma }}{v}_{\mathrm{ci}}^2{\mathrm{\Delta t}}_i}{\underset{i=1}{\overset{n}{\Sigma }}{\mathrm{\Delta t}}_i+\mathrm{\Delta t}}}\begin{array}{c}{t}_{k=1}<\mathrm{\Delta t}<t_k\\ {\mathrm{\Delta t}}_k=t_k-\mathrm{\Delta t}\end{array}$(2)主测线、联络线视倾角不同引起的均方根速度闭合差校正①从主测线、联络测线叠加剖面上，分别拾取出剖面的时间斜率d₁、d₂，②利用地震测线分出主测线、联络测线的叠加速度V_a1、V_a2；③根据如下公式计算主测线、联络测线的视倾角θ₁和θ₂、均方根速度V_R1、V_R2；$V_{R1}=\frac{V_{a1}}{\sqrt{1+{\left(\frac{V_{a1}{d}_1}{2}\right)}^2}},$$V_{R2}=\frac{V_{a2}}{\sqrt{1+{\left(\frac{V_{a2}{d}_2}{2}\right)}^2}}$${\theta }_1={\cos }^{-1}\left[\frac{1}{\sqrt{1+{\left(\frac{V_{a1}{d}_1}{2}\right)}^2}}\right],$${\theta }_2={\cos }^{-1}\left[\frac{1}{\sqrt{1+{\left(\frac{V_{a2}{d}_2}{2}\right)}^2}}\right]$④计算主测线与联络测线交点处的加权系数c₁(i，j)、c₂(i，j)$c_1(i,j)=\frac{{\sin }^2{\theta }_1(i,j)+{\cos }^2{\theta }_2(i,j)}{2}$$c_2(i,j)=\frac{{\cos }^2{\theta }_1(i,j)+{\sin }^2{\theta }_2(i,j)}{2}$⑤对交点处的加权系数c₁(i，j)、c₂(i，j)各自进行插值，得到各自加权系数散点c₁(x，y)c₂(x，y)。⑥归一化处理最终各点的加权系数为：$\begin{array}{cc}\frac{c_2(x,y)}{c_1(x,y)+c_2(x,y)}& \frac{c_1(x,y)}{c_1(x,y)+c_2(x,y)}\end{array}$⑦对主测线和联络线均方根速度插值得到各自均方根速度散点v_R1(x，y)v_R2(x，y)⑧加权叠加得到最终的均方根速度场：$v_R(x,y)=\frac{c_1(x,y)}{c_1(x,y)+c_2(x,y)}{v}_{R1}(x,y)+\frac{c_2(x,y)}{c_1(x,y)+c_2(x,y)}{v}_{R2}(x,y).$