[发明专利]基于佐布里兹方程约束的地震数据规则化处理方法

权利要求书

1.基于佐布里兹方程约束的地震数据规则化处理方法，其特征在于，该基于佐布里兹方程约束的地震数据规则化处理方法包括：

步骤1，对输入的炮集记录进行道头扫描，确定正常炮点间距及炮点缺失位置；

步骤2，根据临近炮点炮检关系，创建缺失位置炮点炮检关系，并为新添加炮点创建数据存储位置；

步骤3，对数据进行道集抽取，抽取共中心点道集，对新增加的炮点创建空道数据；

步骤4，根据已有数据进行振幅高次函数拟合，获得高次函数权值系数；

步骤5，使用FK域凸集投影方法进行数据规则化处理，并进行收敛判断；

步骤6，所有共中心点道集规则化处理结束后，输出规则化处理结果；

在步骤3中，对原有及新补数据进行道集抽取，抽取共中心点道集，并对共中心点数据集按照方位角和偏移距进行排序，构建出三维数据体，如果数据是二维的，则共中心点仅按照偏移距排序，构建的数据体为二维；

抽取共中心点道集，并为新增加炮点创建空道；创建新炮点时，同时创建了数据道头和空道数据，并在炮表中标注当前新增加的炮点为特殊标志，用以后期数据规则化处理的时候使用；

在入射角不是特别大的时候，通过高次函数基本可以拟合振幅变化规律，采用高次函数来替代佐布里兹方程对振幅能量进行约束；

在步骤5中，使用FK域凸集投影算法对共中心点道集数据进行规则化处理，处理采用阈值迭代算法进行处理；进行处理之前首先使用已有数据道振幅对数据振幅变化进行高次函数拟合，拟合公式如下：

a(t)、b(t)、c(t)是随时间变化的权系数，x_i为当前数据对应偏移距；

通过已知数据拟合，计算出各时间位置权值，待插值数据的振幅根据公式(1)进行预估，阈值迭代中，实际插值数据与期望振幅数据误差最小时，即如公式(2)所示情况下，输出对应的a(t)、b(t)、c(t)，并带入(1)式获得实际输出值，

I_i(t)为规则化处理后某偏移距处插值振幅结果，Q_min(t)是指Q取得的最小值；

在步骤4中，在每个共中心点道集内，每个时间位置均首先进行振幅高次函数拟合，拟合使用公式(1)进行；具体高次项项数根据处理精度和实际数据分布情况而定。

2.根据权利要求1所述的基于佐布里兹方程约束的地震数据规则化处理方法，其特征在于，在步骤1中，进行道头扫描，扫描获得炮点X,Y坐标，然后需要根据观测系统分布方位进行整体旋转，重置炮点X,Y坐标。

3.根据权利要求1所述的基于佐布里兹方程约束的地震数据规则化处理方法，其特征在于，在步骤2中，对缺失位置增加炮点，炮点对应的检波点，采用临近炮点对应炮检关系，或根据观测系统模板重新设置炮检关系。

4.根据权利要求1所述的基于佐布里兹方程约束的地震数据规则化处理方法，其特征在于，在步骤5中，在数据规则化处理过程中，加入内部数据交叉验证算法，即使用已知数据对数据规则化参数进行内部验证，保证数据插值参数的可靠性。

5.根据权利要求1所述的基于佐布里兹方程约束的地震数据规则化处理方法，其特征在于，在步骤6中，在共中心点域进行数据规则化处理后，将数据重新抽取到原始炮集记录，炮内道序按照其炮检关系排列，并将数据重新保存。