同相轴拾取方法及装置与流程

1.本发明涉及油气地震勘探技术领域，尤其涉及同相轴拾取方法及装置。


背景技术：

2.本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
3.层析速度反演在地震勘探领域是一种有效的速度建模手段，可以利用接收到的地震波信息对地下介质的参数进行反演，构建能够反映地下介质低频特征的速度场。层析速度反演方法有很多分支，成像域层析反演因其稳定性和有效性被广泛使用。在成像域进行层析速度反演，要拾取角度域成像点道集的同相轴曲率信息。拾取方法一般包括手动拾取和自动拾取方法。
4.然而，现有常规的同相轴拾取方法不仅拾取效率不高，而且拾取结果也不够精确。


技术实现要素：

5.本发明实施例提供一种同相轴拾取方法，用以提高同相轴拾取效率及拾取精确度，该同相轴拾取方法包括：
6.利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线；
7.根据成像道集内每个同相轴的拟合曲线确定成像道集的拟合曲线；
8.确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。
9.本发明实施例还提供一种同相轴拾取装置，用以提高同相轴拾取效率及拾取精确度，该同相轴拾取装置包括：
10.同相轴拟合模块，用于利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线；
11.道集曲线拟合模块，用于根据成像道集内每个同相轴的拟合曲线确定成像道集的拟合曲线；
12.最终拟合确定模块，用于确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。
13.本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述同相轴拾取方法。
14.本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述同相轴拾取方法的计算机程序。
15.本发明实施例中，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线；根据成像道集内每个同相轴的拟合曲线确定成像道
集的拟合曲线；确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。本发明实施例通过二次函数拟合同相轴，可以不受地震资料分辨率的限制，能够很好的与真实同相轴吻合，提高同相轴拾取精度；同时同相轴拾取无需人工干扰和辅助，能够提高同相轴拾取效率。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
17.图1为本发明实施例提供的同相轴拾取方法的实现流程图；
18.图1-1为本发明实施例提供的角度域成像点道集的扫描示意图；
19.图1-2为本发明实施例提供的成像道集的示意图；
20.图1-3为本发明实施例提供的最终拾取的成像道集的同相轴曲线示意图；
21.图2为本发明实施例提供的同相轴拾取方法的另一实现流程图；
22.图3为本发明实施例提供的同相轴拾取方法中步骤201的实现流程图；
23.图4为本发明实施例提供的同相轴拾取方法的又一实现流程图；
24.图5为本发明实施例提供的同相轴拾取方法中步骤101的实现流程图；
25.图5-1为本发明实施例提供的围绕零角度顶点建立的坐标系示意图；
26.图5-2为本发明实施例提供的角度域共成像点道集中同相轴的初始拟合同相轴道集的示意图；
27.图6为本发明实施例提供的同相轴拾取方法中步骤501的实现流程图；
28.图7为本发明实施例提供的同相轴拾取方法中步骤102的实现流程图；
29.图8为本发明实施例提供的同相轴拾取装置的功能模块图；
30.图9为本发明实施例提供的同相轴拾取装置的另一功能模块图；
31.图10为本发明实施例提供的同相轴拾取装置中预处理模块901的结构框图；
32.图11为本发明实施例提供的同相轴拾取装置的又一功能模块图；
33.图12为本发明实施例提供的同相轴拾取装置中同相轴拟合模块801的结构框图；
34.图13为本发明实施例提供的同相轴拾取装置中初始道集拟合单元1201的结构框图；
35.图14为本发明实施例提供的同相轴拾取装置中道集曲线拟合模块802的结构框图。
具体实施方式
36.为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
37.图1示出了本发明实施例提供的同相轴拾取方法的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
38.如图1所示，同相轴拾取方法，其包括：
39.步骤101，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线；
40.步骤102，根据成像道集内每个同相轴的拟合曲线确定成像道集的拟合曲线；
41.步骤103，确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。
42.在层析速度建模中，成像域层析成像是一种比较成熟的速度建模方法。成像域层析成像根据提取的共成像点道集获得预测模型与真实模型之间的差引起的深度残差，通过一系列方法转化成深度残差从而建立方程反演速度。其中，共成像点道集同相轴的拾取是十分重要的，有较高的精度要求。
43.在初始速度模型和真实速度模型存在差异时，角度域成像点道集是一条弯曲的曲线。图1-1示出了角度域成像点道集的扫描示意，图1-2示出了成像道集的示意。在拾取地震道集的同相轴时，首先获取成像点道集。例如，接收输入的角度域成像点道集。进而，利用二次函数对角度域成像点道集中的同相轴进行拟合。其中，角度域成像点道集包含多个同相轴。通过二次函数拟合的方式，可以获得每个同相轴的拟合曲线。
44.然而，在通过二次函数拟合获得成像道集内每个同相轴的拟合曲线后，这些拟合曲线中包含有无效拟合曲线，需要对其进行过滤，以获得有效的拟合曲线。
45.最后，鉴于成像道集的拟合曲线是无限连续的，在获得成像道集的拟合曲线后，确定拟合曲线的顶底边界，进而将顶底边界范围内的拟合曲线作为最终拾取的曲线。图1-3示出了本发明实施例提供的最终拾取的成像道集的同相轴曲线示意。
46.在本发明实施例中，利用二次函数对每个同相轴进行拟合，获得每个同相轴的拟合曲线，确定拟合曲线的顶底边界，将顶底边界范围内的拟合曲线作为最终拾取的曲线。本发明实施例通过二次函数拟合同相轴，可以不受地震资料分辨率的限制，能够很好的与真实同相轴吻合，提高同相轴拾取精度；同时同相轴拾取无需人工干扰和辅助，能够提高同相轴拾取效率。
47.图2示出了本发明实施例提供的同相轴拾取方法的另一实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
48.在本发明的一实施例中，为了进一步提高同相轴拾取精度，如图2所示，在上述图1所示方法步骤的基础上，同相轴拾取方法，还包括：
49.步骤201，对角度域成像点道集进行预处理，获得预处理后的角度域成像点道集；
50.步骤101，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线，包括：
51.步骤202，利用二次函数对预处理后的角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线。
52.在拾取同相轴道集的过程中，可以首先对获取的角度域成像点道集进行预处理，以获得预处理后的角度域成像点道集。进而，利用二次函数对预处理后的角度域成像点道集中每个同相轴进行拟合，据此可以提高同相轴拾取精度。
53.在本发明实施例中，对角度域成像点道集进行预处理，获得预处理后的角度域成像点道集，进而利用二次函数对预处理后的角度域成像点道集中每个同相轴进行拟合，获
得每个同相轴的拟合曲线，能够进一步提高同相轴拾取精度。
54.图3示出了本发明实施例提供的同相轴拾取方法中步骤201的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
55.在本发明的一实施例中，预处理包括以下一种或多种：频率滤波、炮点切除、自动增益控制。为了提高同相轴拾取精度，如图3所示，步骤201，对角度域成像点道集进行预处理，获得预处理后的角度域成像点道集，包括：
56.步骤301，对角度域成像点道集进行频率滤波，获得频率滤波后的角度域成像点道集；和/或
57.步骤302，对角度域成像点道集进行炮点切除，获得炮点切除后的角度域成像点道集；和/或
58.步骤303，对角度域成像点道集进行自动增益控制，获得自动增益控制后的角度域成像点道集。
59.预处理可以包括频率滤波及炮点切除等，本领域技术人员可以理解的是，还可以包括除上述频率滤波及炮点切除之外的其它预处理方式，例如自动增益控制等，本发明实施例对此不作特别的限制。
60.为了减少人工成本，考虑使用先拾取后过滤的采集策略。首先对整个道集进行去噪、频率滤波、炮点切除及自动增益控制等预处理，使道集能量均衡，更容易识别。因为观测系统的限制，在数据的顶部和边缘可能出现振幅与有效部位振幅差距较大的同相轴，为了排除这些干扰，将上下左右部分窄带区域排除在处理范围之外。
61.其中，可以对角度域成像点道集进行上述频率滤波、炮点切除及自动增益控制等预处理中的一种。例如，对角度域成像点道集仅进行频率滤波，或仅进行炮点切除或自动增益控制等一种预处理。还可以对角度域成像点道集同时进行上述频率滤波、炮点切除及自动增益控制等预处理中的多种。例如，对角度域成像点道集同时进行频率滤波、炮点切除及自动增益控制三种预处理等。
62.在本发明实施例中，对角度域成像点道集进行频率滤波，获得频率滤波后的角度域成像点道集；和对角度域成像点道集进行炮点切除，获得炮点切除后的角度域成像点道集；对角度域成像点道集进行自动增益控制，获得自动增益控制后的角度域成像点道集，能够分别提高同相轴拾取精度。
63.图4示出了本发明实施例提供的同相轴拾取方法的又一实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
64.在本发明的一实施例中，二次函数的二次项系数和常数项分别反映拟合后同相轴道集的弯曲程度及纵向偏移量。为了在保证分辨率的前提下提高同相轴拾取效率，如图4所示，在上述方法步骤的基础上，同相轴拾取方法，还包括：
65.步骤401，通过记录弯曲程度及纵向偏移量存储每个同相轴的拟合曲线。
66.其中，用于拟合同相轴的二次拟合函数可以表示为：
67.y＝ax2+b；
68.其中，y表示拟合同相轴，x表示角度域成像点道集同相轴在坐标系中的角度，a反映拟合同相轴的弯曲程度，b反映拟合同相轴在坐标系中的纵向偏移量。
69.在拾取到拟合同相轴之后还需要对其进行储存，以往的道集存储模式是存储道集
的几个特征点，在使用时对道集进行插值，这种方式不仅分辨率有限，还增加了计算成本。我们从另外的角度出发，假设每一个同相轴都能使用二次函数拟合，在存储时只需要存储道集零点的位置和弯曲程度(参数)的值即可记录整个道集。这种存储方法可以在不损失分辨率的情况下压缩冗余数据，提高同相轴拾取效率。
70.在本发明实施例中，通过记录弯曲程度及纵向偏移量存储每个同相轴的拟合曲线，能够在保证分辨率的前提下提高同相轴拾取效率。
71.图5示出了本发明实施例提供的同相轴拾取方法中步骤101的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
72.在本发明的一实施例中，为了提高拟合同相轴的准确性，进而提高同相轴拾取精度，如图5所示，步骤101，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线，包括：
73.步骤501，围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得与每个同相轴对应的多个初始拟合曲线；
74.步骤502，对每个同相轴对应的每个初始拟合曲线的振幅和绝对值进行比较，将振幅和绝对值最大的拟合曲线作为每个同相轴的拟合曲线。
75.层析反演假设速度的扰动很小，所以角度域共成像点道集同相轴的曲率也很小，同相轴比较平缓，适合二次函数拟合。当层速度估计值相对层速度真实值偏小时，角度域共成像点道集的同相轴(二次函数形态)向上弯，层速度估计值相对层速度真实值偏大时，角度域共成像点道集的同相轴(二次函数形态)向下弯。在零角度处总会存在一个顶点，可以围绕零角度处顶点建立如图5-1所示的坐标系。进而，利用二次函数对角度域成像点道集中每个同相轴进行拟合。其中，在拟合时每个同相轴可以得到对应的多个初始拟合曲线。
76.在预处理之后对角度域共成像点道集进行扫描，选取一个零角度点作为扫面的起始点，从一个负值开始对这个零角度点上的道集形态进行拟合，如图5-2所示。不断增加的值，使拟合函数的曲率变化。通过改变扫描步长可以提高拟合精度，但效率会下降。根据实验，增量值一般取10-3
到10-4
为适宜。
77.当速度值偏小时，零角度是角度域共成像点道集的顶点，通过改变拟合曲线的弯曲程度a的值来改变拟合曲线的弯曲程度，之后和真实的同相轴进行比对，得到和真实同相轴最为适配的拟合曲线。
78.同相轴的形态和子波的形态类似，一般会出现一个波峰和前后两个波谷。如果初始拟合曲线和真实同相轴形态是匹配的，初始拟合曲线经过路径的振幅和相互叠加，得到一个极大值；相反，如果初始拟合曲线的形态和同相轴形态不匹配，拟合同相轴道集经过的路径上的振幅值必定有正有负，叠加之后会有一个比较小的值，通过这个规律，就可以判断出初始拟合曲线的形态是否正确。
79.因此，可以对每个同相轴对应的每个初始拟合曲线路径上的振幅的绝对值之和进行比较，同时将振幅和绝对值最大的初始拟合曲线作为同相轴的拟合曲线。
80.在本发明实施例中，围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得与每个同相轴对应的多个初始拟合曲线，对每个同相轴对应的每个初始拟合曲线的振幅和绝对值进行比较，将
振幅和绝对值最大的拟合曲线作为同相轴的拟合曲线。本发明实施例通过二次函数拟合同相轴，利用振幅和绝对值确定同相轴的拟合曲线，能够提高拟合同相轴的准确性，进而提高同相轴拾取精度。
81.图6示出了本发明实施例提供的同相轴拾取方法中步骤501的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
82.在本发明的一实施例中，为了提高初始拟合曲率的准确性，进而提高同相轴拾取精度，如图6所示，步骤501，围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得与每个同相轴对应的多个初始拟合曲线，包括：
83.步骤601，围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用第一二次函数对角度域成像点道集中每个同相轴角度不大于零的部分进行拟合；利用第二二次函数对角度域成像点道集中每个同相轴角度大于零的部分进行拟合，获得每个同相轴对应的初始拟合曲线。
84.由于地下反射界面倾角等因素的影响，角度域共成像点道集的同相轴具有特殊性，其形态可能发生畸变。考虑用左右两个二次函数进行拟合同相轴道集，左右两边使用不同的弯曲程度(参数)，可以适应同相轴道集畸变的情况。
85.因此，在拟合每个同相轴时，首先围绕角度域成像点道集中每个同相轴道集的零角度顶点建立如图5-1的坐标系，进而利用第一二次函数对角度域成像点道集中每个同相轴角度不大于零的部分进行拟合，确定初始拟合同相轴的左侧部分。然后，利用第二二次函数对角度域成像点道集中每个同相轴角度大于零的部分进行拟合，确定初始拟合同相轴的右侧部分。
86.最后，将初始拟合曲线的左侧部分与初始拟合曲线的右侧部分进行合并，确定每个同相轴对应的初始拟合曲线。
87.其中，初始拟合曲线可以表示如下：
[0088][0089]
其中，y表示拟合曲线，a1表示拟合曲线左侧的弯曲程度，a2表示拟合曲线右侧的弯曲程度，b表示拟合曲线在坐标系上的纵向偏移量，x表示角度域成像点道集同相轴在坐标系中的角度。
[0090]
在本发明实施例中，围绕角度域成像点道集中每个同相轴道集的零角度顶点建立坐标系，利用第一二次函数对角度域成像点道集中每个同相轴角度不大于零的部分进行拟合；利用第二二次函数对角度域成像点道集中每个同相轴角度大于零的部分进行拟合，获得每个同相轴对应的拟合曲线。本发明实施例利用第一二次函数及第二二次函数，对角度域成像点道集中每个同相轴角度不大于零的部分及角度大于零的部分分别进行拟合，得到每个同相轴对应的拟合曲线，能够减轻甚至消除道集畸变，提高拟合曲线的准确性，进而提高同相轴拾取精度。
[0091]
图7示出了本发明实施例提供的同相轴拾取方法中步骤102的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：
[0092]
在本发明的一实施例中，为了进一步提高同相轴拾取精度，如图7所示，步骤102，
根据成像道集内每个同相轴的拟合曲线确定成像道集的拟合曲线，包括：
[0093]
步骤701，将成像道集内同相轴拟合曲线的振幅平均值作为成像道集的拾取过滤阈值；
[0094]
步骤702，将振幅不小于拾取过滤阈值的拟合曲线，作为成像道集的拟合曲线。
[0095]
在确定成像道集的拟合曲线时，扫描整个成像道集获得成像道集内每个同相轴的拟合曲线。然后，计算获得成像道集内每个同相轴的拟合曲线的振幅平均值，将拟合曲线的振幅平均值作为成像道集的过滤阈值。利用该过滤阈值，将振幅不小于成像道集过滤阈值的拟合曲线，作为成像道集的拟合曲线，而将振幅小于成像道集过滤阈值的拟合曲线滤除。
[0096]
在本发明实施例中，将成像道集内每个同相轴的拟合曲线的振幅平均值作为成像道集的过滤阈值，将振幅不小于成像道集过滤阈值的拟合曲线，作为成像道集的拟合曲线，本发明实施例通过将拟合曲线的振幅平均值作为成像道集的过滤阈值，过滤成像道集内的拟合曲线，能够进一步提高同相轴拾取精度。
[0097]
本发明实施例还提供一种同相轴拾取装置，如下面的实施例所述。由于这些装置解决问题的原理与同相轴拾取方法相似，因此这些装置的实施可以参见方法的实施，重复之处不再赘述。
[0098]
图8示出了本发明实施例提供的同相轴拾取装置的功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0099]
参考图8，所述同相轴拾取装置所包含的各个模块用于执行图1对应实施例中的各个步骤，具体请参阅图1以及图1对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述同相轴拾取装置包括同相轴拟合模块801、道集曲线拟合模块802及最终拟合确定模块803。
[0100]
同相轴拟合模块801，用于利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线。
[0101]
道集曲线拟合模块802，用于确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。
[0102]
最终拟合确定模块803，用于确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。
[0103]
在本发明实施例中，同相轴拟合模块801利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线，道集曲线拟合模块802确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线，最终拟合确定模块803，用于确定成像道集拟合曲线的顶底边界，将成像道集剖面顶底边界范围内的拟合曲线作为最终拾取的拟合曲线。本发明实施例通过二次函数拟合同相轴，可以不受地震资料分辨率的限制，能够很好的与真实同相轴吻合，提高同相轴拾取精度；同时同相轴拾取无需人工干扰和辅助，能够提高同相轴拾取效率。
[0104]
图9示出了本发明实施例提供的同相轴拾取装置的另一功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0105]
在本发明的一实施例中，为了进一步提高同相轴拾取精度，参考图9，所述同相轴拾取装置所包含的各个模块用于执行图2对应实施例中的各个步骤，具体请参阅图2以及图2对应实施例中的相关描述，此处不再赘述。本发明实施例中，在上述图8所示模块结构的基
础上，所述同相轴拾取装置还包括预处理模块901。
[0106]
预处理模块901，用于对角度域成像点道集进行预处理，获得预处理后的角度域成像点道集。
[0107]
同相轴拟合模块801，还用于利用二次函数对预处理后的角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线。
[0108]
在本发明实施例中，预处理模块901对角度域成像点道集进行预处理，获得预处理后的角度域成像点道集，进而同相轴拟合模块801利用二次函数对预处理后的角度域成像点道集中每个同相轴进行拟合，获得成像道集内每个同相轴的拟合曲线，能够进一步提高同相轴拾取精度。
[0109]
图10示出了本发明实施例提供的同相轴拾取装置中预处理模块901的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0110]
在本发明的一实施例中，预处理包括以下一种或多种：频率滤波、炮点切除、自动增益控制。为了提高同相轴拾取精度，参考图10，所述预处理模块901所包含的各个单元用于执行图3对应实施例中的各个步骤，具体请参阅图3以及图3对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述预处理模块901包括频率滤波单元1001、炮点切除单元1002及自动增益控制单元1003。
[0111]
频率滤波单元1001，用于对角度域成像点道集进行频率滤波，获得频率滤波后的角度域成像点道集。和/或
[0112]
炮点切除单元1002，用于对角度域成像点道集进行炮点切除，获得炮点切除后的角度域成像点道集。和/或
[0113]
自动增益控制单元1003，用于对角度域成像点道集进行自动增益控制，获得自动增益控制后的角度域成像点道集。
[0114]
在本发明实施例中，频率滤波单元1001对角度域成像点道集进行频率滤波，获得频率滤波后的角度域成像点道集；炮点切除单元1002和对角度域成像点道集进行炮点切除，获得炮点切除后的角度域成像点道集；自动增益控制单元1003对角度域成像点道集进行自动增益控制，获得自动增益控制后的角度域成像点道集，能够分别提高同相轴拾取精度。
[0115]
图11示出了本发明实施例提供的同相轴拾取装置的又一功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0116]
在本发明的一实施例中，二次函数的二次项系数和常数项分别反映拟合后曲线的弯曲程度及纵向偏移量。为了在保证分辨率的前提下提高同相轴拾取效率，参考图11，所述同相轴拾取装置所包含的各个模块用于执行图4对应实施例中的各个步骤，具体请参阅图4以及图4对应实施例中的相关描述，此处不再赘述。本发明实施例中，在上述模块结构的基础上，所述同相轴拾取装置，还包括存储模块1101。
[0117]
存储模块1101，用于通过记录弯曲程度及纵向偏移量存储每个同相轴的拟合曲线。
[0118]
在本发明实施例中，存储模块1101通过记录弯曲程度及纵向偏移量存储每个同相轴的拟合曲线，能够在保证分辨率的前提下提高同相轴拾取效率。
[0119]
图12示出了本发明实施例提供的同相轴拾取装置中同相轴拟合模块801的结构示
意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0120]
在本发明的一实施例中，为了提高拟合同相轴的准确性，参考图12，所述同相轴拟合模块801所包含的各个单元用于执行图5对应实施例中的各个步骤，具体请参阅图5以及图5对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述同相轴拟合模块801包括初始道集拟合单元1201及拟合道集确定单元1202。
[0121]
初始道集拟合单元1201，用于围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得与每个同相轴对应的多个初始拟合曲线。
[0122]
拟合道集确定单元1202，用于对每个同相轴对应的每个初始拟合曲线的振幅和绝对值进行比较，将振幅和绝对值最大的拟合曲线作为每个同相轴的拟合曲线。
[0123]
在本发明实施例中，初始道集拟合单元1201围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用二次函数对角度域成像点道集中每个同相轴进行拟合，获得与每个同相轴对应的多个初始拟合曲线，拟合道集确定单元1202对每个同相轴对应的每个初始拟合曲线的振幅和绝对值进行比较，将振幅和绝对值最大的拟合曲线作为每个同相轴的拟合曲线。本发明实施例通过二次函数拟合同相轴，利用振幅和绝对值确定同相轴的拟合曲线，能够提高拟合同相轴的准确性，进而提高同相轴拾取精度。
[0124]
图13示出了本发明实施例提供的同相轴拾取装置中初始道集拟合单元1201的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0125]
在本发明的一实施例中，为了提高初始拟合同相轴道集的准确性，进而提高同相轴拾取精度，参考图13，所述初始道集拟合单元1201所包含的各个子单元用于执行图6对应实施例中的各个步骤，具体请参阅图6以及图6对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述初始道集拟合单元1201包括拟合子单元1301。
[0126]
拟合子单元1301，用于围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用第一二次函数对角度域成像点道集中每个同相轴角度不大于零的部分进行拟合；利用第二二次函数对角度域成像点道集中每个同相轴角度大于零的部分进行拟合，获得每个同相轴对应的初始拟合曲线。
[0127]
在本发明实施例中，拟合子单元1301围绕角度域成像点道集中每个同相轴的零角度顶点建立坐标系，利用第一二次函数对角度域成像点道集中每个同相轴角度不大于零的部分进行拟合；利用第二二次函数对角度域成像点道集中每个同相轴角度大于零的部分进行拟合，获得每个同相轴对应的初始拟合曲线。本发明实施例利用第一二次函数及第二二次函数，对角度域成像点道集中每个同相轴角度不大于零的部分及角度大于零的部分分别进行拟合，得到每个同相轴道集对应的初始拟合曲线，能够减轻甚至消除道集畸变，提高初始拟合曲线的准确性，进而提高同相轴拾取精度。
[0128]
图14示出了本发明实施例提供的同相轴拾取装置中道集曲线拟合模块802的结构示意，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0129]
在本发明的一实施例中，为了进一步提高同相轴拾取精度，参考图14，所述道集曲线拟合模块802所包含的各个单元用于执行图7对应实施例中的各个步骤，具体请参阅图7以及图7对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述道集曲线拟合模块802包括过滤阈值确定单元1401及道集确定单元1402。
[0130]
过滤阈值确定单元1401，用于将成像道集内同相轴拟合曲线的振幅平均值作为成像道集的拾取过滤阈值。
[0131]
道集确定单元1402，用于将振幅不小于拾取过滤阈值的拟合曲线，作为成像道集的拟合曲线。
[0132]
在本发明实施例中，过滤阈值确定单元1401将成像道集内同相轴拟合曲线的振幅平均值作为成像道集的拾取过滤阈值，道集确定单元1402将振幅不小于拾取过滤阈值的拟合曲线，作为成像道集的拟合曲线，本发明实施例通过将拟合同相轴的振幅平均值作为成像道集的过滤阈值，过滤成像道集内的拟合曲线，能够进一步提高同相轴拾取精度。
[0133]
在成像域层析中，拾取地下层位的目的是为了寻找能够和角度域道集对应的反射点信息，如果拾取的层位没有反射点对应，这个点含有的信息就是无效的。本发明使用拾取的角度域成像点道集作为依据，将角度域共成像点道集的零角度点确定为反射点，符合成像点道集的假设。同时，扫描到的角度域共成像点道集可以根据振幅和等信息进行匹配，寻找到同一层位的反射轴，为层析提供地层倾角数据。在不拾取层位的情况下也能找到反射点和走时残差的关系，避免了各种转换之间的精度损失。
[0134]
通过二次函数拟合角度域共成像点道集同相轴方法可以获得不受地震资料分辨率限制的高精度深度残差信息，并且拥有很高的效率，无需人工干扰，提高了层析速度反演的效率。
[0135]
同常规的同相轴拾取方法相比，二次函数拟合扫描方法考虑了角度域共成像点道集本身特点，从形态和波形两个方面入手，从整个角度域共成像点道集上扫描同相轴。速度不准确时的角度域共成像点道集呈现的双曲线，通过二次函数对其进行拟合，可以比较好的还原同相轴形态。结合拟合路径上的振幅信息，与背景噪音进行对比和区分，除去干扰同相轴，最终达到精确拾取同相轴的目的。
[0136]
该方法具有以下优点：1、有别于其他同相轴拾取方法，本发明充分考虑共成像点道集同相轴特点，利用二次函数进行拟合；2、充分利用角度域共成像点道集的振幅信息，受噪音干扰小；3、利于拟合的方式扫描，存储空间少，拟合分辨率超过角度域共成像点道集分辨率，能很好的与同相轴吻合，拾取出准确的深度残差，为后续反演建立准确的方程组，从而得到准确的速度更新量。
[0137]
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述同相轴拾取方法。
[0138]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述同相轴拾取方法的计算机程序。
[0139]
综上所述，本发明实施例中，利用二次函数对每个同相轴进行拟合，获得每个同相轴的拟合曲线，确定成像道集拟合曲线的顶底边界，将顶底边界范围内的拟合曲线作为最终拾取的曲线。本发明实施例通过二次函数拟合同相轴，可以不受地震资料分辨率的限制，能够很好的与真实同相轴吻合，提高同相轴拾取精度；同时同相轴拾取无需人工干扰和辅助，能够提高同相轴拾取效率。
[0140]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实
施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0141]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0142]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0143]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0144]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。