[发明专利]一种地震数据处理VIA双参数分析方法

权利要求书

1.一种地震数据处理VIA双参数分析方法，其特征在于，包括如下步骤：

步骤S1，布置地震数据采集观测系统，包括设置炮点、检波点的位置以及二者的排列关系，利用地震数据采集观测系统采集地震数据，并设置地震波传播相关参数，其中，所述地震波传播相关参数，包括：炮点S坐标为x_s，检波点R坐标为x_r，l为炮点S与检波点R之间的距离，即偏移距为l＝x_r-x_s，v为均匀各向同性介质中的传播速度，l₀表示反射点O处反射界面的法线在地面出射点C与炮点的距离，地震波入射角和反射角相等，用β表示，θ₁为反射点O处反射界面的倾角；

步骤S2，通过炮点S、检波点R和反射点O构建一个圆，对于圆上每一个反射点，入射线和反射线之间的夹角固定不变，在每个反射点处做角平分线，该角平分线即为反射点处反射界面的法线，所有法线相交于一点，该交点为极点；

令t表示地震波沿入射波射线的旅行时t₁和反射波射线的旅行时t₂的总旅行时；t₀表示法向射线OC段的双程旅行时；令表示法向射线延长到P点即OP段对应的双程旅行时，v表示地震波速度；l表示炮点S与检波点R之间的距离，l₀表示反射点O处反射界面的法线在SR上的出射点C与炮点S的距离；计算的VIA双参数分析方法的变换方程为：

所述计算的VIA双参数分析方法的变换方程，包括如下步骤:

已知如下数学关系：

由圆的性质，有如下关系式：

把(2)中的关系代入(1),得到:

通过计算可以求取圆的直径d：

d＝l/(2sin2β)(4)

进而可以推出如下关系式，

利用公式(1)-(4)，可以推出

结合椭圆展开变换方程，有

公式(6)和(7)联立,可以推出

进一步可以推出v和β双参数关系式：

步骤S3，通过上述变换方程，计算生成速度分析界面，速度分析界面包括速度谱和叠加能量剖面，通过速度分析获得速度值；(1)生成速度谱；速度谱用于速度分析估算速度参数值；(2)将t转换为t₀，实现地震数据的展开等时线的相应振幅的绝对值叠加，获得叠加能量剖面，该剖面可以确定主要地层的层数和大致位置，用于辅助进行速度分析；

步骤S4，根据步骤S3得到的速度，通过步骤S2的变换方程进行入射角扫描地震数据成像，生成入射角扫描道集；

步骤S5，根据步骤S4得到的入射角扫描道集，分析估算不同位置、不同成像时间参与计算的入射角参数值范围。

2.如权利要求1所述的地震数据处理VIA双参数分析方法，其特征在于，在所述步骤S2中，进一步采用参数展开共反射点叠加成像方式，所有展开等时线，当双参数正确时，等时线具有公切点，形成同相叠加；其它等时线相互干涉形成变换噪音。

3.如权利要求1所述的地震数据处理VIA双参数分析方法，其特征在于，设θ₂为法向射线与垂向线的夹角，其中，θ₂＝θ₁。