[发明专利]一种基于复波场的方位角度域偏移成像方法及系统

权利要求书

1.一种基于复波场的方位角度域偏移成像方法，其特征在于：所述方法首先进行炮点波场、检波点波场的复波场的延拓，然后在各个空间窗内实现复波场的局部角度分解，获得炮点波场、检波点波场的局部角度分量，最后对不同方位的角度域道集进行叠加，获得最终的方位角域成像结果。

2.根据权利要求1所述的基于复波场的方位角度域偏移成像方法，其特征在于:所述方法包括：

步骤一：输入数据：震源、震源HT、炮集、炮集HT；

步骤二：进行炮点复波场外推得到炮点复波场，同时进行检波点复波场外推得到检波点复波场；

步骤三：分别将炮点复波场、检波点复波场分解为多个相同的空间窗，然后获得每一个空间窗的所有角度分量；

步骤四，利用所有空间窗的角度分量获得方位角域成像结果。

3.根据权利要求2所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤二的操作包括：

利用下面的公式(1)进行炮点复波场外推得到炮点复波场，利用下面的公式(2)进行检波点复波场外推得到检波点复波场：

其中，Us和Ur分别为炮端波场和检波端波场，v为速度，t为时间，fs为震源子波，fs^h为做完Hilbert变换后的子波，dr和dr^h为炮集和做完Hilbert变换后的炮集，与为炮点复波场与检波点复波场。

4.根据权利要求3所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤三的操作包括：

(1)：将当前空间窗内的复波场变换到波数域；

(2)：计算波数域内每个斜率上的叠加能量；

(3)：筛选出最大叠加能量角度分量；

(4)：将所述最大叠加能量角度分量反变换到空间域，并求得误差泛函；

(5)：判断所述误差泛函的值是否小于设定的阈值，如果否，则转入步骤(6)，如果是，则移动到下一个空间窗，然后返回步骤(1)；

(6)：输出当前角度分量，然后转入步骤(1)。

5.根据权利要求4所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤(1)的操作包括：

利用下面的公式(5)将当前空间窗内的复波场变换到波数域：

其中K_sx＝(k_sx,k_sz)为空间波数矢量，win为空间窗大小，包括横向和纵向；t，x分别为步骤一输入的数据中的时间和空间。

6.根据权利要求5所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤(3)的操作包括：

从步骤(2)计算得到的叠加能量中找到最大的叠加能量，该最大叠加能量对应的斜率即为最大能量角度分量。

7.根据权利要求6所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤(4)中的将所述最大叠加能量角度分量反变换到空间域的操作包括：

利用公式(6)将所述最大叠加能量角度分量进行傅立叶反变换：

其中，其中，代表局部空间波场，

θs表示自变量，θ₀为步骤(3)得到的最大叠加能量角度。

8.根据权利要求7所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤(4)中的求得误差泛函的操作包括：

利用公式(4)计算得到误差泛函：

9.根据权利要求8所述的基于复波场的方位角度域偏移成像方法，其特征在于：所述步骤四的操作包括：

利用下面的公式得到方位角域成像结果：

其中，

10.一种基于复波场的方位角度域偏移成像系统，其特征在于：所述系统包括：

输入单元，用于输入数据：震源、震源HT、炮集、炮集HT；

波场延拓单元，与所述输入单元连接，用于进行炮点复波场外推得到炮点复波场，同时进行检波点复波场外推得到检波点复波场；

角度分量获取单元，与所述波场延拓单元连接，用于分别将炮点复波场、检波点复波场分解为多个相同的空间窗，然后获得每一个空间窗的所有角度分量；

成像单元，与所述角度分量获取单元连接，用于利用所有空间窗的角度分量获得方位角域成像结果。