CN103217711B - 获取共偏移距共方位角域道集的方法及系统 - Google Patents

Abstract

本发明涉及获取共偏移距共方位角域道集的方法及系统，将地震数据分选成n*m个交叉排列地震数据道集；确定每个交叉排列的中心点，地震数据坐标系平移获取新的坐标系，新的坐标系原点为中心点；将每个平行四边形的两平行边旋转至与新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；对矩形反射点数据域划分共偏移距共方位角单元，并进行编号；收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集。

Description

获取共偏移距共方位角域道集的方法及系统

技术领域

本发明涉及地震数据处理领域，特别涉及获取共偏移距共方位角域道集的方法及系统。

背景技术

基于共偏移距共方位角域的宽方位地震数据处理技术是目前工业界研究和应用的热点，主要用于宽方位地震数据的叠前成像、分方位各向异性处理和裂缝预测。它不同于传统的共偏移距域地震数据，传统的共偏移距域地震数据偏移后方位角信息丢失。共偏移距共方位角域地震数据偏移是在限定偏移距和方位角的前提下进行偏移，更能提高偏移结果的精度，提高横向分辨率。偏移后道集具有方位角信息，有利于开展叠前反演，也使得方位AVO（AZAVO）分析成为可能，具有非常广阔的应用前景。

形成共偏移距共方位角域道集是这一技术的关键步骤。目前的技术只适用于在正交观测系统下形成共偏移距共方位角域道集，不能解决广泛使用的斜交观测系统和其他更复杂的观测系统中共偏移距共方位角域道集生成问题。

发明内容

本发明的目的是针对现有技术不适应斜交观测系统的局限性，提出一种在斜交观测系统中获取共偏移距共方位角域道集的方法及系统，解决了在斜交观测系统下共偏移距共方位角域宽方位地震数据处理技术应用难题。

为实现上述目的，本发明提供了在斜交观测系统中获取共偏移距共方位角域道集的方法，包括：

将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集。

可选的，在本发明一实施例中，所述对所述矩形反射点数据域划分共偏移距共方位角单元的步骤包括：

根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

可选的，在本发明一实施例中，所述将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点的步骤包括：

根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角。

可选的，在本发明一实施例中，所述将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点的步骤包括：

根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角。

可选的，在本发明一实施例中，所述获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域的步骤包括：

获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

可选的，在本发明一实施例中，所述获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域的步骤包括：

获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

为实现上述目的，本发明还提供了在斜交观测系统中获取共偏移距共方位角域道集的系统，包括：

地震数据处理单元，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

坐标系平移变换单元，用于确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

矩形反射点数据域获取单元，包括旋转模块和映射模块，所述旋转模块用于将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；所述映射模块用于获取经所述旋转模块处理后与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

划分编号单元，用于对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

共偏移距共方位角域道集获取单元，用于收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集。

可选的，在本发明一实施例中，所述划分编号单元根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

可选的，在本发明一实施例中，所述旋转模块根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角。

可选的，在本发明一实施例中，所述旋转模块根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角。

可选的，在本发明一实施例中，所述映射模块获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

可选的，在本发明一实施例中，所述映射模块获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

上述技术方案具有如下有益效果：本申请的技术方案解决了在斜交观测系统中获取共偏移距共方位角域道集的难题，使得斜交观测系统可以成功应用共偏移距共方位角域宽方位处理技术。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提出的获取共偏移距共方位角域道集的方法流程图；

图2为本发明提出的获取共偏移距共方位角域道集的系统方框图；

图3为实施例中交叉排列地震数据道集反射点分布图；

图4为实施例中一交叉排列经过原点平移后的反射点分布图；

图5为实施例中经过平移的平行四边形经坐标旋转后的反射点分布图；

图6为实施例中经过平移和旋转后的平行四边形映射后的反射点分布图；

图7为实施例中进行共偏移距共方位角域数据单元划分以及编号后的示意图；

图8为实施例中在平行四边形上共偏移距共方位角域数据单元划分以及编号在地震数据坐标系的示意图；

图9为现有技术的共偏移距共方位角域数据单元划分以及编号后在地震数据坐标系的示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，为本发明提出的获取共偏移距共方位角域道集的方法流程图。该方法包括：

步骤101）：将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

步骤102）：确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

步骤103）：将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

步骤104）：对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

步骤105）：收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集。

在图1中，步骤105）具体包括：根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

步骤103）包括两种情况，第一种情况为：根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角。在此情况下，获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

第二种情况为：根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角。在此情况下，获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

如图2所示，为本发明提出的获取共偏移距共方位角域道集的系统方框图。包括：

地震数据处理单元201，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

坐标系平移变换单元202，用于确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

矩形反射点数据域获取单元203，包括旋转模块2031和映射模块2032，所述旋转模块2031用于将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；所述映射模块2032用于获取经所述旋转模块处理后与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

划分编号单元204，用于对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

共偏移距共方位角域道集获取单元205，用于收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集。

优选地，所述划分编号单元204根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

在图2中，旋转模块2031一种情况下根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角。此种情况下，所述映射模块2032获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

旋转模块2031另一种情况下根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角。此种情况下，所述映射模块获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。实施例：

首先，地震数据准备，包括原始地震数据读入、地震数据编辑、观测系统加载和预处理等；

然后，将三维地震数据分选成交叉排列地震数据道集；

在本实施例中，实际地震资料观测系统共有36条炮线，67条检波线，所以交叉排列地震数据数目为36×67个，即2412个；如图3所示，为实施例中交叉排列地震数据道集反射点分布图。从图3可知，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形，在本实施例中，有2412个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形。

确定每个交叉排列的中心点坐标，将新坐标原点移到交叉排列的中心点位置；

在本实施例中，有36条炮线，67条检波线，每个交叉排列的中心点为炮线与检波线的交点，例如，2412个交叉排列的中心点中的一个中心点坐标在地震数据坐标系中为（X₀,Y₀），将对应地交叉排列中每个反射点数据的原有坐标减去中心点坐标（X₀,Y₀）即可以实现将交叉排列中心点平移至新的坐标系原点位置；如图4所示，为实施例中交叉排列经过原点平移后的反射点分布图。从图4可知，新的坐标系原点位于平行四边形的中心。

将平行四边形的两平行边旋转到与坐标系中一坐标轴至平行状态；在本实施例中，如图5所示，为实施例中经过平移的平行四边形经坐标旋转后的反射点分布图。根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角。量出图5中的平行四边形的α角，获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，在图5中，A点、B点是水平坐标轴与平行四边形的另外两平行边的交点，保持A、B和坐标原点不动，应用映射公式 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right),$ 将平行四边形反射点数据域映射为矩形反射点数据域，如图6所示。

还有一种情况，根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角。在此情况下，获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

如图7所示，为实施例中进行共偏移距共方位角域数据单元划分以及编号后的示意图。根据炮线距、检波线距和覆盖次数等信息，在矩形反射点数据域上进行共偏移距共方位角域数据单元划分以及编号。

收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即形成共偏移距共方位角域道集。

如图8所示，为实施例中在平行四边形上共偏移距共方位角域数据单元划分以及编号在地震数据坐标系的示意图；如图9所示，为现有技术的共偏移距共方位角域数据单元划分以及编号后在地震数据坐标系的示意图。图8和图9对比可知，本申请提出的技术方案划分结果正确无误。现有技术的划分结果在边界处单元形状不规则，说明现有技术划分不合理。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.在斜交观测系统中获取共偏移距共方位角域道集的方法，其特征在于，包括：

将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集；

其中，所述将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点的步骤包括：

根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角；

所述获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域的步骤包括：

获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

2.根据权利要求1所述的方法，其特征在于，所述对所述矩形反射点数据域划分共偏移距共方位角单元的步骤包括：

根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

3.在斜交观测系统中获取共偏移距共方位角域道集的方法，其特征在于，包括：

将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集；

其中，所述将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点的步骤包括：

根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角；

所述获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域的步骤包括：

获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

4.根据权利要求3所述的方法，其特征在于，所述对所述矩形反射点数据域划分共偏移距共方位角单元的步骤包括：

根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

5.在斜交观测系统中获取共偏移距共方位角域道集的系统，其特征在于，包括：

地震数据处理单元，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

坐标系平移变换单元，用于确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

矩形反射点数据域获取单元，包括旋转模块和映射模块，所述旋转模块用于将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；所述映射模块用于获取经所述旋转模块处理后与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

划分编号单元，用于对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

共偏移距共方位角域道集获取单元，用于收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集；

其中，所述旋转模块根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& \sin \mathrm{\θ}\\ -\sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的水平坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的底边与水平坐标轴的夹角；

所述映射模块获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}1& -\mathrm{ctg\α}\\ 0& \csc \mathrm{\α}\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

6.根据权利要求5所述的系统，其特征在于，所述划分编号单元根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。

7.在斜交观测系统中获取共偏移距共方位角域道集的系统，其特征在于，包括：

地震数据处理单元，用于将地震数据分选成n*m个交叉排列地震数据道集；其中，n为在斜交观测系统中炮线的数目，m为在斜交观测系统中检波线的数目，每个交叉排列地震数据道集反射点分布在地震数据坐标系中均表现为平行四边形；

坐标系平移变换单元，用于确定每个交叉排列的中心点，所述地震数据坐标系平移获取新的坐标系，新的坐标系原点为所述中心点；其中，所述中心点为炮线与检波线的交点，也为平行四边形的中心点，中心点数目为n*m；

矩形反射点数据域获取单元，包括旋转模块和映射模块，所述旋转模块用于将每个平行四边形的两平行边旋转至与所述新的坐标系的坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；所述映射模块用于获取经所述旋转模块处理后与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，将平行四边形反射点数据域映射为矩形反射点数据域；

划分编号单元，用于对所述矩形反射点数据域划分共偏移距共方位角单元，并进行编号；

共偏移距共方位角域道集获取单元，用于收集每个矩形反射点数据域的矩形网格中同一象限相同编号的共偏移距共方位角单元，即获取共偏移距共方位角域道集；

其中，所述旋转模块根据旋转公式 $\left(\begin{array}{cc}\cos \mathrm{\θ}& -\sin \mathrm{\θ}\\ \sin \mathrm{\θ}& \cos \mathrm{\θ}\end{array}\right)$ 将每个平行四边形的两平行边旋转至与所述新的坐标系的纵坐标轴平行，旋转后的平行四边形的中心点位所述新的坐标系的原点；其中，θ为平行四边形的侧边与纵坐标轴的夹角；

所述映射模块获取与平行四边形的两平行边平行的坐标轴同平行四边形的另外两平行边的交点，保持交点和坐标原点位置不变，根据映射关系 $\left(\begin{array}{cc}\csc \mathrm{\α}& 0\\ -\mathrm{ctg\α}& 1\end{array}\right)$ 将平行四边形反射点数据域映射为矩形反射点数据域；其中，α为平行四边形为锐角的内角角度。

8.根据权利要求7所述的系统，其特征在于，所述划分编号单元根据炮线距、检波线距和覆盖次数对矩形反射点数据域划分共偏移距共方位角单元以及对每个单元进行编号。