双复杂构造区叠前深度偏移速度建模方法及装置

1.本技术涉及地震勘探技术领域，特别涉及双复杂构造区叠前深度偏移速度建模方法及装置。


背景技术：

2.双复杂构造具体是指地表构造和地下构造双重复杂的地区，主要分布于山地地震采集区域，如前陆冲断带发育逆演推覆构造，伴随着逆冲叠瓦状构造样式，形成断块、断鼻、挤压背斜等构造。双复杂构造区地震波场十分复杂，基于水平层状介质假设的水平叠加及叠后偏移技术存在一定的问题。叠前时间偏移是建立在双平方根方程的非零偏移距成像理论，它沿非零炮检距的绕射曲线旅行轨迹对振幅求和，对每个共炮检距剖面单独成像，把反射波能量转移到真实的地下位置处，保证分析点速度代表地下成像点位置处的真实均方根速度，在一定程度上解决了不同倾角地层速度多解性问题，但还是无法满足双复杂构造成像精度的要求。叠前深度偏移技术，尤其是基于双程波旅行时的波动方程逆时偏移技术，可适应地下地震波场横向变化剧烈、高陡构造成像困难、高速地质体导致垂向速度反转等问题，有效提高双复杂构造区偏移成像质量，为勘探开发提供优质精准的基础资料。但是，叠前深度偏移对速度建模的精度要求非常高。因此，如何提高双复杂构造区叠前深度偏移速度建模精度，历来是地球物理学术界和工业界研究的热点和难点，本发明正是为了解决该问题。


技术实现要素：

3.本技术的目的是提供双复杂构造区叠前深度偏移速度建模方法及装置，以提高双复杂构造区叠前深度偏移速度建模精度。
4.为解决上述技术问题，本说明书第一方面提供一种双复杂构造区叠前深度偏移速度建模方法，包括：根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度；根据各目标直井的测井数据，确定各目标直井上的多个层速度值，并将所述层速度值作为层速度值集合中的元素；对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面；重复执行以下方法，直至达到预定条件后得到所述第二速度层的第二目标速度体：在相邻两层速度趋势面之间做插值，得到新的速度趋势面；采用层析反演不断迭代修正所述新的速度趋势面。
5.在一些实施例中，在对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面之前，包括：根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；确定第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；将所述交点的层速度值作为所述层速度值集合的元素。
6.在一些实施例中，根据各目标直井的测井数据，确定各目标直井上的多个层速度
值，包括采用以下方法确定各目标直井上的多个层速度值：根据当前目标直井的gr曲线值区分砂泥岩标志层；对当前目标直井进行岩心归位和岩性标定，并确定各标志层的分界面；确定当前目标直井与各标志层分界面的交点；根据当前目标直井的声波测井曲线数据，确定各交点处的层速度值。
7.在一些实施例中，所述预定条件包括所述第二速度层的速度趋势面的准确度达到预定准确度阈值；其中，所述第二速度层的速度趋势面的准确度通过以下方法计算得到：根据目标直井的测井数据，确定深度域初始cip道集；通过cip道集动校正后的同相轴拉平程度，确定所述第二速度层的速度趋势面的准确度。
8.在一些实施例中，在确定所述第二速度层的第二目标速度体之后，还包括：计算所述层速度值集合中各层速度值的实际值、所述实际值对应的拟合值之间的差值的平方和，作为剩余值；其中所述拟合值是指第二目标速度体中各速度趋势面上对应的层速度值；将所述剩余值作为所述第二速度层的第二目标速度体的一部分。
9.在一些实施例中，所述方法还包括：通过微测井数据和大炮初至数据，建立所述第一速度层的第一目标速度体。
10.在一些实施例中，在基于所述层速度值集合，确定所述第二速度层的速度趋势面之前，包括：根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；将所述初始速度体中第一速度层的部分替换为所述第一速度层的第一目标速度体；确定替换后的第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；将所述交点的层速度值作为所述层速度值集合的元素。
11.本说明书第二方面提供一种双复杂构造区叠前深度偏移速度建模装置，包括：第一确定单元，用于根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度；第二确定单元，用于根据各目标直井的测井数据，确定各目标直井上的多个层速度值，并将所述层速度值作为层速度值集合中的元素；拟合单元，用于对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面；迭代单元，用于重复执行以下方法，直至达到预定条件后得到所述第二速度层的第二目标速度体：在相邻两层速度趋势面之间做插值，得到新的速度趋势面；采用层析反演不断迭代修正所述新的速度趋势面。
12.在一些实施例中，所述装置还包括：第三确定单元，用于根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；第四确定单元，用于确定第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；第五确定单元，用于将所述交点的层速度值作为所述层速度值集合的元素。
13.在一些实施例中，所述第二确定单元包括：第一确定子单元，用于根据当前目标直井的gr曲线值区分砂泥岩标志层；第二确定子单元，用于对当前目标直井进行岩心归位和岩性标定，并确定各标志层的分界面；第三确定子单元，用于确定当前目标直井与各标志层分界面的交点；第四确定子单元，用于根据当前目标直井的声波测井曲线数据，确定各交点处的层速度值。
14.在一些实施例中，所述预定条件包括所述第二速度层的速度趋势面的准确度达到预定准确度阈值；其中，所述第二速度层的速度趋势面的准确度通过以下装置计算得到：根据目标直井的测井数据，确定深度域初始cip道集；通过cip道集动校正后的同相轴拉平程
度，确定所述第二速度层的速度趋势面的准确度。
15.在一些实施例中，所述装置还包括：计算单元，用于计算所述层速度值集合中各层速度值的实际值、所述实际值对应的拟合值之间的差值的平方和，作为剩余值；其中所述拟合值是指第二目标速度体中各速度趋势面上对应的层速度值；第六确定单元，用于将所述剩余值作为所述第二速度层的第二目标速度体的一部分。
16.在一些实施例中，所述装置还包括：第七确定单元，用于通过微测井数据和大炮初至数据，建立所述第一速度层的第一目标速度体。
17.在一些实施例中，在对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面之前，包括：第八确定单元，用于根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；替换单元，用于将所述初始速度体中第一速度层的部分替换为所述第一速度层的第一目标速度体；第九确定单元，用于确定替换后的第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；第十确定单元，用于将所述交点的层速度值作为所述层速度值集合的元素。
18.本说明书第三方面提供一种电子设备，包括：存储器和处理器，所述处理器和所述存储器之间互相通信连接，所述存储器中存储有计算机指令，所述处理器通过执行所述计算机指令，从而实现第一方面任一项所述方法的步骤。
19.本说明书第四方面提供一种计算机存储介质，所述计算机存储介质存储有计算机程序指令，所述计算机程序指令被执行时实现第一方面任一项所述方法的步骤。
20.本说明书提供了双复杂构造区叠前深度偏移速度建模方法及装置，本方案对于近地表浅层使用微测井和大炮初至数据，采用层间相关系数法建立表层速度模型。对于中深层的地层，基于目标直井的测井数据和地震处理速度谱数据联合确定种子点层速度值集合，然后对种子点层速度值集合中的层速度值进行多项式拟合以确定较深地层的速度趋势面；多次迭代，每次迭代时，在相邻两层速度趋势面之间做插值、分小层，并采用层析反演不断修正新的速度趋势面。由此可见，该方法通过控制多项式拟合参数，调整速度分界面与实际地质情况的拟合程度，能够有效控制宏观大幅度构造形态。同时，通过表层模型替换精确速度，外加中深层速度趋势面约束的层析反演建模，获得更加精确的初始深度偏移速度模型。并在后续通过多次迭代修正速度模型过程中，不断调整和增加速度趋势面的数量和形态、细化速度趋势面约束尺度和范围，逐步刻画微观小幅度构造形态，直到最终获得高精度叠前深度偏移速度模型，满足更复杂偏移算法对偏移成像的需求。
附图说明
21.为了更清楚地说明本技术实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
22.图1示出了本说明书提供的一种双复杂构造区叠前深度偏移速度建模方法的流程图；
23.图2示出了根据本说明书所提供的双复杂构造区叠前深度偏移速度建模方法所确定的模型的示意图；
24.图3示出了本说明书提供的另一种双复杂构造区叠前深度偏移速度建模方法的流程图；
25.图4示出了本说明书提供的又一种双复杂构造区叠前深度偏移速度建模方法的流程图；
26.图5示出了微测井的示意图；
27.图6示出了微测井数据的示意图；
28.图7示出了采用空变的层间形态相关系数进行插值操作的示意图；
29.图8示出了本说明书提供的一种双复杂构造区叠前深度偏移速度建模装置的原理框图；
30.图9示出了本说明书提供的一种电子设备的结构示意图。
具体实施方式
31.为了使本技术领域的人员更好地理解本技术中的技术方案，下面将结合本技术实施方式中的附图，对本技术实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本技术一部分实施方式，而不是全部的实施方式。基于本技术中的实施方式，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施方式，都应当属于本技术保护的范围。
32.为突破中西部含油气盆地双复杂构造区地震资料品质差、偏移成像精度低的技术难题，最大限度地发挥叠前深度偏移技术在双复杂构造区的成像优势，迫切需要建立一套系统的、科学的、有效的叠前深度偏移速度建模方法，提高该区地震资料的偏移成像精度，进一步提升地震资料品质，满足勘探开发需求，助力油气增储上产。
33.本说明书提供一种双复杂构造区叠前深度偏移速度建模方法，如图1所示，该方法包括如下步骤：
34.s10：根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度。
35.微测井，是指通过井中激发、地面接收，或者地面激发、井中接收，或者井中激发、井中接收等方式采集地震波信息，并进一步求取近地表地球物理参数的方法。微测井方法中所涉及到的数据即为微测井数据。如图2所示，包括微测井1和微测井2。
36.如图2所示，第一速度层可以包括图2中的低速层和降速层，第二速度层可以是图2中的高速层。微测井的深度可以达到第一速度层与第二速度层的分界面，也即微测井的深度可以达到降速层与高速层的分界面。从而，根据微测井数据能够较为准确地确定第一速度层与第二速度层的分界面，能够较为准确地确定图2中低速层与降速层的分界面、降速层与高速层的分界面。
37.层速度指在层状地层中地震波传播的速度。层速度直接反映地层的岩性，能用来划分地层。一般是用地震测井或声波测井测得，并且指的是纵波的速度；也可以利用反射记录计算得到。在地震勘探中，一般把层速度低于1400米/秒的地层叫做低速层，把高于3500米/秒的地层视为高速层，层速度位于1400-3500米/秒之间的地层视为降速层。
38.将地震采集获得的微测井数据进行岩性标定并重新解释，精细划分每口微测井处的低速层、降速层、高速层的分界线。具体地，不同层的声波传播速度不同，也就是测井曲线
的斜率不同，因此可以通过测井曲线的斜率来区分不同层。由于微测井数量较少，以往建模的方法多采用厚度线性内插法，这种方法受人为因素影响较大，所解释的结果带有不确定性。本说明书所提供的建模方法，在微测井无法确定层速度的空间范围内，使用大炮初至拾取，计算层速度，进而划分低速层、降速层、高速层。
39.s20：根据各目标直井的测井数据，确定各目标直井上的多个层速度值，并将所述层速度值作为层速度值集合中的元素。
40.在本说明书中，层速度值集合中的元素也称为种子点层速度。种子点是指层速度值集合中各层速度值对应的空间位置点。
41.直井，与定向井相对应，是指设计轨道是一条铅垂线的井。在直井的设计轨迹中，理论上轨迹上所有点的井斜角都为零，但在实际钻井中是有误差的，只要误差在合格井的允许范围内，不管实际井眼的井斜角是多少，它仍叫直井。
42.为保证叠前深度偏移速度建模精度，原则上目标直井的井点坐标位置到目标速度体的纵向目标线垂直距离必须小于100米。
43.微测井与直井可以选址在双复杂构造区的地表凸点、地表凹点。
44.如图2所示，直井的深度可以达到高速层中，因此，可以通过直井的测井数据进行高速层的速度建模。
45.直井的测井数据通常在深度维度上分布较为密集，为提高建模效率，在进行高速层的速度建模时无需采用如此密集的数据，为此，首先需要从密集的直井测井数据中筛选出建模所需的数据集，具体地，可以筛选出目标速度体的横向目标线所在的深度位置的直井测井数据形成建模所需的数据集。
46.对此，在一些实施例中，s20包括采用以下方法确定各目标直井上的多个层速度值：
47.s21：根据当前目标直井的gr曲线值区分砂泥岩标志层。
48.s22：对当前目标直井进行岩心归位和岩性标定，并确定各标志层的分界面。
49.s23：确定当前目标直井与各标志层分界面的交点。
50.s24：根据当前目标直井的声波测井曲线数据，确定各交点处的层速度值。
51.标志层，是指已知地质年代并具有明显特征(如岩石的性质、颜色、厚度及所含的化石、杂质等特征)、易于辨认的岩层，可作为地层对比标志。
52.如图2所示，通过上述步骤s21至s24，可以确定目标直井a上的a1、a2、a3、a4点处的层速度值，可以确定目标直井b上的b1、b2、b3、b4、b5点处的层速度值，可以确定目标直井c上的c1、c2、c3、c4点处的层速度值。
53.由于双复杂构造区地表形态较为复杂，不便于设置很多直井，因此仅根据目标直井上多个点处的层速度值所建立的速度趋势面精度较低。为提高所建立的速度趋势面的精度，可以通过叠前时间偏移速度谱确定叠前深度偏移的初始速度体，并基于该初始速度体扩充层速度值集合中的层速度值的数量，通过扩充后的层速度值集合确定第二速度层的速度趋势面。扩充的层速度值(即初始速度体中的层速度值)可能不是很精准，但是其含有变化趋势信息，能够辅助提高速度趋势面的准确度。相应地，如图3所示，在步骤s30之前，还包括如下步骤s70、s80和s90。
54.s70：根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度
体。
55.将叠前时间偏移速度体，经dix公式转换得到叠前深度偏移速度体。
56.本说明书中的“速度体”是指在空间中分布的速度值的集合，即速度体中包含了多个速度值，各速度值具有对应的空间位置。
57.s80：确定第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值。
58.通过智能化层位追踪拾取技术可以快速获得解释精度较高的、具有较强地震反射特征的多套地质层位，如图2所示，主要包括：1、高速层顶界(i类/ii类标志层)；2、区域可追踪且层位起伏变化平缓(i类标志层)，如水平沉积、单斜构造等；3、区域可追踪且层位起伏变化剧烈(ii类标志层)，如前积构造、逆冲构造等；4、局部可追踪的形态各异的地质体反射界面(iii类标志层)，如潜山构造、盐丘构造等。
59.目标速度体的纵横目标线，包括横向目标线和纵向目标线(即与直井几乎平行的目标线)，是缩减模型的网格线，相当于空间中的虚拟网格线。横向目标线位于垂直于纵向目标线平面上。
60.步骤s80可以为，在确定第二目标速度体的纵横目标线与层位分界线的交点坐标之后，可以确定在初始速度体中该交点坐标对应的速度值，并将初始速度体中该交点坐标的速度值作为第二目标速度体中该交点坐标的速度值。
61.s90：将所述交点的层速度值作为所述层速度值集合的元素。
62.综合步骤s20和s80-s90，层速度值集合中的元素既包括了目标直井上的层速度值，又包括了根据初始速度体得到的层速度值。从而，步骤s30所确定的速度趋势面，既可以利用到目标直井上的层速度确定速度趋势面的关键点，建立速度趋势面骨架，又可以利用初始速度体所含有的变化趋势信息确定速度趋势面的变化趋势。
63.上述步骤s70、s80和s90所确定的层速度值可以仅用于确定第二速度层的第二速度体。第一速度层为近地表层，便于设置多个微测井，因此，可以通过微测井数据和大炮初至数据，为第一速度层建立较为准确的第一速度体。
64.在一些实施例中，如图4所示，在步骤s80之前，还可以包括如下步骤s100：将所述初始速度体中第一速度层的部分替换为所述第一速度层的第一目标速度体。由于第一目标速度体是根据微测井数据和大炮初至数据确定的，因此，相对于初始速度体中的第一速度层部分，第一目标速度体更为准确。将初始速度体中第一速度层的部分替换为第一目标速度体，可以便于步骤s80直接根据替换后的初始速度体确定层速度值集合，而无需再筛选出初始速度体中的第二速度层部分，以根据第二速度层部分确定第二速度体的速度趋势面。
65.在一些实施例中，在确定第一速度体时，可以基于微测井数据确定各观测点的深度和接收到地震波的时间，基于该深度和该时间生成时距图，然后进行解释检验和重新解释。图5示出了微测井的示意图，其中s所对应的井中设置多个地震源，r所对应的井中设置有多个观测点，x0表示井s和井r之间的距离，q1、q2……qn
表示观测点。图6示出了微测井数据的示意图，其中，time表示观测点接收到地震波的时间，trace表示观测点的编号，可以根据各观测点的编号及深度确定观测点相对于激发点的距离。时距图，是指地震波从震源出发，传播到观测线上各观测点的传播时间与观测点相对于激发点(即坐标原点)距离之间的关系。直达波的时距曲线是一条直线。
66.根据各观测点对应的深度和时间，做时深曲线。各深度点对应的时间可以表示为tj(h)＝ti+h/vi，其中tj表示目标深度点对应的时间，hi表示目标深度点所在层的上一层的地层总深度，h表示目标深度点在其所在层的深度，vi表示目标深度点所在层的层速度，j表示目标深度点的编号，i表示地层的编号。例如，低速层的深度为10m，目标深度点位于降速层，且在降速层的深度为6m，那么目标深度点对应的时间为tj(h)＝10+6/v
降速层
。也就是说，目标深度点所对应的时间可以根据关联层的层速度及厚度重新求取。
67.依据时深曲线上曲线的斜率划分多个斜率段，根据划分的斜率段对应划分地层，对各层的深度点做直线回归拟合，求出回归系数a、b和平均标准差s。在深度点数大于3的情况下，可剔除个别与拟合线偏差大于平均标准差2.5倍的点，然后再作最后的直线回归拟合，得到最终的a、b值。
68.系数a也即直线的斜率，系数b也即直线的截距。对应于公式：tj(h)＝ti+h/vi，1/vi是直线的斜率，对应系数a，速度单位为m/s，若将速度单位转换成m/ms，则速度值为1/a
×
1000；hi是直线的截距，对应系数b。
69.第i层的时深曲线为：t
i+1
(h)＝ti+h/vi。根据两个时深曲线在相交处的深度相等，可以得到：ti+h/vi＝t
i+1
+h/v
i+1
，左侧为第i层的深度，右侧为第i+1层的深度。那么可以得到：
70.第i层的厚度(注：t的单位ms)。通过该方法可以求取各层的厚度。
71.为得到微测井数据所设置的观测点的数量通常是有限的，因此仅通过微测井数据只能够确定地层分层的轮廓，该地层分层轮廓虽然比第二速度体准确，但是对于一些分层细节依然无法确定。本说明书提出采用空变的层间形态相关系数进行插值操作，以进一步优化建模结果。
72.如图7所示，对于一个层界面(图7中的加粗黑色曲线表示层界面)，可以通过上述方法求得层界面上的两个控制点hm、h
m+1
，对于这两个控制点之间的各点可以采用插值方法确定。
73.具体地，对于各层界面，控制点m、m+1之间的控制点层界面高程whj的公式为：
74.wh
0j
＝fh
0j
+dh
0j
×
k0ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(第一层的底界)
75.wh
1j
＝fh
1j
+dh
1j
×
k1ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(第二层的底界)
76.wh
2j
＝fh
2j
+dh
2j
×
k2ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(第三层的底界)
77.……
78.其中，wh表示实测为测井的物理点，fh表示内插得到的点，dh是插值点所对应的上一层的地层深度，该点处的最终深度如何，需要调节k值确定。
79.每一道地震道的地表海拔高度。各地震道均有高程，当然肯定有某一地震道的高层是最大的，该道理论上应该为最接近工区海拔最高点的位置。
80.k0为层间相关系数，即低速层底界面与地表的相关系数。
81.k值大小与层间相关度成正比，k＝0，不相关，k＝1，强相关。
82.k值的确定方法为：
83.a)根据经验，建模前定义，批量应用；
84.b)根据一组或一个局部区域内表层控制点上形态的层间相关度统计求取出空间变化的相关系数；
85.c)对特殊地段采用交互建模，即通过交互改变相关系数动态生成模型，直到满意为止。具体地，可以为：主测线方向计算一次，同时联络测线方向再计算一次，两次计算深度误差值足够小，或者小于某个门槛值，比如半个采样间隔，就证明达到建模精度要求。
86.s30：对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面。
87.如图2所示，基于各目标直井上各点处的层速度值，可以确定第二速度层的速度趋势面，图2中加粗黑线即代表速度趋势面。
88.根据层速度值集合，按照纵横向深度域速度建模目标线网格间距，构建工区层速度和主要标志层位的栅状速度剖面格架，并输出当前全部离散点数据。
89.在此基础上，加载当前速度值、层位数据、井点坐标等全部离散点数据进入速度建模软件，并设置网格化参数(方位角、平滑半径、网格化方法、网格边长和网格边界缩进距离等)后；将离散点按照方位角和网格边界缩进距离转换坐标系，以正方形网格划分整个区域从而得到网格点x、y坐标，落在网格点平滑半径内的离散点按照网格化方法(距离加权平均法)计算出该网格点z值，如果网格点z值不允许为空值，需要逐级扩大平滑半径直到有离散点落在平滑半径内；如需要对网格数据进行平滑和加密，可对其进行平滑滤波(box平滑)和加密插值，使网格点z值能够充分平滑；以色块图的形式检查网格数据是否满足需要，如满足需要则保存到网格文件并与种子点层速度数据库建立关联，以色块图的形式检查相邻最近的层速度种子点的坐标是否能够投影到不同的网格之中，如果不能则证明网格边长偏大，需要设置更密的网格，重新进行网格化。其中，种子点是指层速度值集合中各速度值对应的空间位置点。
90.因以往地震解释方案带有一定的主观性和多解性，所以在新采集地震资料处理的深度偏移速度建模过程中，需要处理人员带着批判的态度，有选择性的使用有关数据。如果完全按照先前解释的构造格架，很可能会误导新一轮的地震资料处理，以至于造成对地下构造成像的“将错就错”。所以，本说明书提供一种基于多项式拟合建立速度趋势面的方法，在宏观尺度范围内使用相对可靠的地震解释方案，通过控制多项式拟合参数，调整速度分界面与实际地质情况的拟合程度，能够有效控制宏观大幅度构造形态。同时，通过表层模型替换精确速度，外加中深层速度趋势面约束的层析反演建模，获得更加精确的初始深度偏移速度模型。并在后续通过多次迭代修正速度模型过程中，不断调整和增加速度趋势面的数量和形态、细化速度趋势面约束尺度和范围，逐步刻画微观小幅度构造形态，直到最终获得高精度叠前深度偏移速度模型，满足更复杂偏移算法对偏移成像的需求。随着多次迭代速度模型精度的提高，每轮迭代速度建模过程中，都根据需要修改原有速度趋势面，用来更好地拟合地下宏观趋势，同时增加新的更小尺度的速度趋势面，用来刻画比宏观构造尺度小一级别的微幅构造。然后在下一轮速度模型更新迭代中，再去增加速度趋势面去刻画更小尺度的微幅构造，直到建模精度符合叠前深度偏移精度要求为止。该方法把偏移成像的主动权交还给更先进的偏移算法，而不是人为设置边界和障碍，以上述网格化为基础，通过多项式拟合速度趋势面约束速度建模，经过多次迭代获得精度更高的最终层速度模型。
91.多项式拟合速度趋势面的原理：从实际数据中分解出“趋势值”和“剩余值”，从而
揭示地理要素空间分布的趋势与规律。
92.a)建立速度趋势面模型
93.以二元函数为例，以z表示某一点初始层速度值，即以zi(xi,yi)，(1，2，...，n)表示，用表示各个点的层速度趋势拟合值，则是平面点坐标(x,y)的函数，即：
[0094][0095][0096]
式中：εi即为速度剩余值(残差值)。
[0097]
显然，当(xi,yi)在空间上变动时，上式(2)就刻画速度要素的实际分布曲面、趋势面和剩余面之间的互动关系。
[0098]
速度趋势面分析的核心是从实际地下观测值出发推算趋势面，一般采用回归分析方法，使得残差平方和趋于最小，即：
[0099][0100]
这就是在最小二乘法意义下的趋势曲面拟合。
[0101]
用来计算趋势面的数学方程式有多项式函数和傅立叶级数，其中最为常用的是多项式函数形式。因为任何一个函数都可以在一个适当的范围内用多项式来逼近，而且调整多项式的次数，可使所求的回归方程适合实际问题的需要。
[0102]
多项式趋势面的形式：
[0103]
(1)一次趋势面模型：(平面：适用于近水平或单斜标志层)
[0104]
z＝a0+a1x+a2y
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0105]
(2)二次趋势面模型：(双曲面：适用于起伏变化较小标志层)
[0106]
z＝a0+a1x+a2y+a3x2+a4xy+a5y2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0107]
(3)三次趋势面模型：(三次曲面：适用于起伏变化剧烈标志层)
[0108]
z＝a0+a1x+a2y+a3x2+a4xy+a5y2+a6x3+a7x2y+a8xy2+a9y3ꢀꢀ
(6)
[0109]
b)趋势面的极值数值位置
[0110]
拟合函数f(x,y)在(x0,y0)点达到极值的必要条件是：
[0111][0112]
以二次为例，设f(x,y)＝2a0+a
1 x+a
2 y+a
3 x2+a
4 xy+a
5 y2，则极值点的地理坐标应该满足以下方程组：
[0113][0114]
上述方程组(8)为一个二元一次方程组，有一组解：
[0115][0116]
即二维二次曲面有一个极值点，由上式(9)所决定。
[0117]
二维三次曲面的回归方程为：
[0118]
z＝a0+a1x+a2y+a3x2+a4xy+a5y2+a6x3+a7x2y+a8xy2+a9y3ꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0119]
其极值点的位置应满足如下的方程组：
[0120][0121]
这是一个二元二次的方程组，应该有4组解，因此，最多有4个极值点。同理，二维四次趋势面最多应有9个极值点。对二维趋势面来说，设趋势面次数最高为i，则极值点个数最大值为(i-1)2。
[0122]
c)趋势面模型的参数估计
[0123]
如果趋势面函数为任意多项式。即：
[0124][0125]
上式中：
[0126]
——为速度趋势值；
[0127]
l——为趋势面的次数；
[0128][0129]
所谓趋势面模型的参数估计就是根据观测值zi(xi,yi)，(1，2，...，n)来确定多项式系数a0,a1,...,a
p
，使残差平方和q为最小。因而，可以将多项式回归问题转化为多元线性回归问题来解决。
[0130]
记：
[0131]
x＝x1,y＝x2,x2＝x3,xy＝x4,y2＝x5,
…ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0132]
则(公式1-12)可以写为
[0133][0134]
这里，残差可以写为：
[0135][0136]
由式(1-14)可以看出，q是a0,a1,...,a
p
的函数。
[0137]
根据最小二乘法的原则，选择使q达到极小这样的系数a0,a1,...,a
p
，(p《n)。为此，需要求q对a0,a1,...,a
p
的偏导数：
[0138][0139]
令即得：
[0140][0141]
在本说明书的公式中，大写的“z”与小写的“z”表示同一含义。(公式1-16)包括p+1个线性方程组，解此线性代数方程组即得p+1个系数a0,a1,...,a
p
。(公式1-16)称为正规方程组，写成矩阵形式，如记为：
[0142][0143]
则正规方程组(16)可以写为：
[0144]
x'xa＝x'z
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0145]
其中，x'为x的转置矩阵：
[0146][0147]
理论上，到这里剩下的工作只是一个计算的问题，也就是说，只是一个如何解正规方程组的问题了。可以证明，a的系数矩阵
[0148]
a＝x'x
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(20)
[0149]
是非奇异的，所以解
[0150]
a＝[x'x]-1
x'z
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(21)
[0151]
是唯一确定的。
[0152]
为了求解a0,a1,...,a
p
，我们首先改变正规方程组(16)的写法，记：
[0153]
[0154]
并引进符号：
[0155][0156][0157][0158]
其中，——某点层速度提取值的平均值；
[0159]
——第j个变量xj的平均值；
[0160]
——z的总方差。
[0161]
在方程组(16)中，由第一个方程，可以得到
[0162][0163]
用上面的符号，可以写成
[0164][0165]
将(23)代入(14)式得：
[0166][0167]
由(24)式出发，求q对aj的偏导数，并使可得方程组：
[0168][0169]
方程组(13)和(16)是等价的。
[0170]
方程组(13)包括p个方程，由此解出a0,a1,...,a
p
，在代入(23)式即得a0，这样就求出了全部系数a0,a1,...,a
p
。
[0171]
d)趋势面模型的拟合精度
[0172]
根据提取值zi及其坐标(xi，yi)进行拟合的某次趋势面，在一定程度上指出了坐标(xi，yi)与观测值zi的关系，显然，不同次数的趋势面对原始数据的逼近程度是有差别的。所以，我们需要对趋势面原始数据的逼近情况进行研究。
[0173]
趋势面的拟合精度是一个极其重要的参数，它是指趋势面方程所代表的数学曲面与实际观测数据面的逼近程度。它随趋势面次数的不同而变化，一般而言，趋势面的次数越高拟合精度越大。但是，趋势面的拟合精度具有特殊性，它并不要求很高的拟合精度，相反过高的拟合精度会导致数学曲面过于逼近实测数据的分布面而难以反映分布的主体规律，达不到描述空间总体趋势的目的。在采用多项式趋势面分析中可通过改变多项式的次数来
控制拟合精度，以达到满意的分析效果。实际应用中，人们一般将趋势面的拟合精度c用下式表示：
[0174][0175]
其中，
[0176]
zi：层速度的提取值；
[0177]
层速度的趋势值；
[0178]
层速度的平均值；
[0179]
为zi的残差平方和；即s
偏
；
[0180]
为层速度数据zi与其平均值之差的平方和，即s
总
。
[0181]
c表示由速度趋势面反映出来的波动平方和占层速度采样数据的总波动平方和中的百分比。c百分比越高，则拟合程度越好，但过高的拟合精度将导致趋势剩余的消失，从而漏掉有价值的异常带。
[0182]
e)趋势面模型的剩余分析
[0183]
趋势面分析的目的包括两方面：一方面要确定某些元素的区域性变化规律，所以要求绘出趋势面图形；另一方面是要确定局部异常带，找到不易发现的低缓异常，由此就需要解出剩余值，作剩余分析。根据剩余值的涵义，它是某点的观测值和该点趋势值作差。将所有的剩余值解出之后也就能够画出剩余值等值线图，可得剩余(残差)值为：
[0184][0185]
由于剩余值的引进，原始数据就被分解为趋势值和剩余值两个部分。一般情况下，趋势值反映区域性大范围内的总体变化情况，剩余值则反映局部性小范围变化特点。两方面综合起来可以有助于更深入地开展层速度方面的研究工作。一次趋势面剩余值，把二次曲面以上的复杂变化都看成局部异常，二次趋势面剩余值则把三次及其以上的曲面变化视为局部异常。因此，由低到高不同次数的趋势面剩余值反映了由大到小不同规模的局部异常特征，利用这种方法可以发现有意义的低缓异常带。
[0186]
基于此，在一些实施例中，在确定所述第二速度层的第二目标速度体之后，还可以包括如下步骤sa1和sa2。
[0187]
sa1：计算所述层速度值集合中各层速度值的实际值、所述实际值对应的拟合值之间的差值的平方和，作为剩余值；其中所述拟合值是指第二目标速度体中各速度趋势面上对应的层速度值。
[0188]
sa2：将所述剩余值作为所述第二速度层的第二目标速度体的一部分。
[0189]
s40：在相邻两层速度趋势面之间做插值，得到新的速度趋势面。
[0190]
具体地，步骤s40可以是对于各纵向目标线执行如下操作：获取当前纵向目标线分别与第一速度趋势面的第一交点、与第二速度趋势面的第二交点，确定第一交点位置的速
度值和第二交点位置的速度值的平均值，确定第一交点与第二交点的连线的中点，将该平均值作为该中点位置的速度值。这只是插值的一个示例，在实际实施的过程中，也可以采用其他插值方法。
[0191]
s50：采用层析反演不断迭代修正所述新的速度趋势面。
[0192]
层析与反演方法，是地震勘探领域的常用方法。通过层析反演，可以对误差数据进行校正。
[0193]
s60：判断是否达到预定条件。若是，则得到第二速度层的第二目标速度体；否则，跳转至步骤s40继续执行。
[0194]
首先按照相邻两层速度趋势面做层速度线性差值，通过分小层速度控制速度渐变快慢，最后通过层析反演不断迭代修正。在未达到预定条件的情况下，重复执行步骤s40和s50，也就是在相邻速度趋势面之间不断分小层。
[0195]
在一些实施例中，所述预定条件包括所述第二速度层的速度趋势面的准确度达到预定准确度阈值。所述第二速度层的速度趋势面的准确度可以通过以下方法计算得到：根据目标直井的测井数据，确定深度域初始cip道集；通过cip道集动校正后的同相轴拉平程度，确定所述第二速度层的速度趋势面的准确度。可以直接将同相轴拉平程度作为第二速度层的速度趋势面的准确度，也可以对同相轴拉平程度做一定的数据处理，得到第二速度层的速度趋势面的准确度。也可以采用上式(26)计算所述第二速度层的速度趋势面的准确度，即将拟合精度c作为第二速度层的速度趋势面的准确度。
[0196]
在速度建模过程中，目标线网格从疏到密，例如通过网格800*800m对的目标线网格建立速度趋势面，可以修正400*400m网格处的速度线处的层速度残差值，逐渐加密迭代修正，最终获得高精度，符合地下地质构造规律的层速度模型。
[0197]
上述建模方法，通过建立速度趋势面，来替代传统的层位约束参数叠前深度偏移层析成像反演速度建模。该建模方法，可以人为通过控制多项式拟合参数，调整速度分界面与实际地质情况的拟合程度，控制宏观大幅度构造形态，滤除微观小幅度构造形态，通过表层模型替换精确速度，外加中深层速度趋势面约束的层析反演建模，得到最终叠前深度偏移速度模型，满足更复杂偏移算法对偏移成像的需求。
[0198]
本说明书所提供的双复杂构造区叠前深度偏移速度建模方法，对于较深的地层，基于目标直井的测井数据确定种子点层速度值集合，然后对种子点层速度值集合中的层速度值进行多项式拟合以确定较深地层的速度趋势面；多次迭代，每次迭代时，在相邻两层速度趋势面之间做插值、分小层，并采用层析反演不断修正新的速度趋势面。由此可见，该方法通过控制多项式拟合参数，调整速度分界面与实际地质情况的拟合程度，能够控制宏观大幅度构造形态，滤除微观小幅度构造形态，通过表层模型替换精确速度，外加中深层速度趋势面约束的层析反演建模，得到最终叠前深度偏移速度模型，满足更复杂偏移算法对偏移成像的需求。
[0199]
本说明书提供的一种双复杂构造区叠前深度偏移速度建模装置，如图8所示，该装置包括第一确定单元10、第二确定单元20、拟合单元30和迭代单元40。
[0200]
第一确定单元10用于根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度。
[0201]
第二确定单元20用于根据各目标直井的测井数据，确定各目标直井上的多个层速
度值，并将所述层速度值作为层速度值集合中的元素。
[0202]
拟合单元30用于对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面。
[0203]
迭代单元40用于重复执行以下方法，直至达到预定条件后得到所述第二速度层的第二目标速度体：在相邻两层速度趋势面之间做插值，得到新的速度趋势面；采用层析反演不断迭代修正所述新的速度趋势面。
[0204]
在一些实施例中，所述装置还包括：第三确定单元，用于根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；第四确定单元，用于确定第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；第五确定单元，用于将所述交点的层速度值作为所述层速度值集合的元素。
[0205]
在一些实施例中，所述第二确定单元包括：第一确定子单元，用于根据当前目标直井的gr曲线值区分砂泥岩标志层；第二确定子单元，用于对当前目标直井进行岩心归位和岩性标定，并确定各标志层的分界面；第三确定子单元，用于确定当前目标直井与各标志层分界面的交点；第四确定子单元，用于根据当前目标直井的声波测井曲线数据，确定各交点处的层速度值。
[0206]
在一些实施例中，所述预定条件包括所述第二速度层的速度趋势面的准确度达到预定准确度阈值；其中，所述第二速度层的速度趋势面的准确度通过以下装置计算得到：根据目标直井的测井数据，确定深度域初始cip道集；通过cip道集动校正后的同相轴拉平程度，确定所述第二速度层的速度趋势面的准确度。
[0207]
在一些实施例中，所述装置还包括：计算单元，用于计算所述层速度值集合中各层速度值的实际值、所述实际值对应的拟合值之间的差值的平方和，作为剩余值；其中所述拟合值是指第二目标速度体中各速度趋势面上对应的层速度值；第六确定单元，用于将所述剩余值作为所述第二速度层的第二目标速度体的一部分。
[0208]
在一些实施例中，所述装置还包括：第七确定单元，用于通过微测井数据和大炮初至数据，建立所述第一速度层的第一目标速度体。
[0209]
在一些实施例中，在对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面之前，包括：第八确定单元，用于根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；替换单元，用于将所述初始速度体中第一速度层的部分替换为所述第一速度层的第一目标速度体；第九确定单元，用于确定替换后的第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；第十确定单元，用于将所述交点的层速度值作为所述层速度值集合的元素。
[0210]
上述装置的具体描述及有益效果请参阅对应的方法实施例，不再赘述。
[0211]
本发明实施例还提供了一种电子设备，如图9所示，该电子设备可以包括处理器901和存储器902，其中处理器901和存储器902可以通过总线或者其他方式连接，图9中以通过总线连接为例。
[0212]
处理器901可以为中央处理器(central processing unit，cpu)。处理器901还可以为其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、
分立硬件组件等芯片，或者上述各类芯片的组合。
[0213]
存储器902作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如本发明实施例中的双复杂构造区叠前深度偏移速度建模方法对应的程序指令/模块(例如，图8所示的迭代单元40)。处理器901通过运行存储在存储器902中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据分类，即实现上述方法实施例中的双复杂构造区叠前深度偏移速度建模方法。
[0214]
存储器902可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器901所创建的数据等。此外，存储器902可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器902可选包括相对于处理器901远程设置的存储器，这些远程存储器可以通过网络连接至处理器901。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
[0215]
所述一个或者多个模块存储在所述存储器902中，当被所述处理器901执行时，执行上述双复杂构造区叠前深度偏移速度建模方法。
[0216]
上述电子设备具体细节可以参阅对应实施例中的相关描述和效果进行理解，此处不再赘述。
[0217]
上述实施方式阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
[0218]
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本技术时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
[0219]
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
[0220]
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
[0221]
以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。

技术特征：
1.一种双复杂构造区叠前深度偏移速度建模方法，其特征在于，包括：根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度；根据各目标直井的测井数据，确定各目标直井上的多个层速度值，并将所述层速度值作为层速度值集合中的元素；对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面；重复执行以下方法，直至达到预定条件后得到所述第二速度层的第二目标速度体：在相邻两层速度趋势面之间做插值，得到新的速度趋势面；采用层析反演不断迭代修正所述新的速度趋势面。2.根据权利要求1所述的方法，其特征在于，在对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面之前，包括：根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；确定第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；将所述交点的层速度值作为所述层速度值集合的元素。3.根据权利要求1所述的方法，其特征在于，根据各目标直井的测井数据，确定各目标直井上的多个层速度值，包括采用以下方法确定各目标直井上的多个层速度值：根据当前目标直井的gr曲线值区分砂泥岩标志层；对当前目标直井进行岩心归位和岩性标定，并确定各标志层的分界面；确定当前目标直井与各标志层分界面的交点；根据当前目标直井的声波测井曲线数据，确定各交点处的层速度值。4.根据权利要求1所述的方法，其特征在于，所述预定条件包括所述第二速度层的速度趋势面的准确度达到预定准确度阈值；其中，所述第二速度层的速度趋势面的准确度通过以下方法计算得到：根据目标直井的测井数据，确定深度域初始cip道集；通过cip道集动校正后的同相轴拉平程度，确定所述第二速度层的速度趋势面的准确度。5.根据权利要求1所述的方法，其特征在于，在确定所述第二速度层的第二目标速度体之后，还包括：计算所述层速度值集合中各层速度值的实际值、所述实际值对应的拟合值之间的差值的平方和，作为剩余值；其中所述拟合值是指第二目标速度体中各速度趋势面上对应的层速度值；将所述剩余值作为所述第二速度层的第二目标速度体的一部分。6.根据权利要求1所述的方法，其特征在于，所述方法还包括：通过微测井数据和大炮初至数据，建立所述第一速度层的第一目标速度体。7.根据权利要求6所述的方法，其特征在于，在基于所述层速度值集合，确定所述第二速度层的速度趋势面之前，包括：根据双复杂构造区的叠前时间偏移速度谱，确定叠前深度偏移的初始速度体；
将所述初始速度体中第一速度层的部分替换为所述第一速度层的第一目标速度体；确定替换后的第二目标速度体的纵横目标线与层位分界面的交点，并根据所述初始速度体确定所述交点的层速度值；将所述交点的层速度值作为所述层速度值集合的元素。8.一种双复杂构造区叠前深度偏移速度建模装置，其特征在于，包括：第一确定单元，用于根据双复杂构造区的微测井数据，确定各微测井上的第一速度层、第二速度层的分界面；所述第二速度层的深度大于第一速度层的深度；第二确定单元，用于根据各目标直井的测井数据，确定各目标直井上的多个层速度值，并将所述层速度值作为层速度值集合中的元素；拟合单元，用于对所述层速度值集合中的层速度值进行多项式拟合，以确定所述第二速度层的速度趋势面；迭代单元，用于重复执行以下方法，直至达到预定条件后得到所述第二速度层的第二目标速度体：在相邻两层速度趋势面之间做插值，得到新的速度趋势面；采用层析反演不断迭代修正所述新的速度趋势面。9.一种电子设备，其特征在于，包括：存储器和处理器，所述处理器和所述存储器之间互相通信连接，所述存储器中存储有计算机指令，所述处理器通过执行所述计算机指令，从而实现权利要求1至7任一项所述方法的步骤。10.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机程序指令，所述计算机程序指令被执行时实现权利要求1至7任一项所述方法的步骤。

技术总结
本说明书提供了双复杂构造区叠前深度偏移速度建模方法及装置，对于中深层的地层，基于目标直井的测井数据和地震处理速度谱数据联合确定种子点层速度值集合，然后对种子点层速度值集合中的层速度值进行多项式拟合以确定较深地层的速度趋势面；多次迭代，每次迭代时，在相邻两层速度趋势面之间做插值、分小层，并采用层析反演不断修正新的速度趋势面。由此可见，该方法通过控制多项式拟合参数，调整速度分界面与实际地质情况的拟合程度，能够控制宏观大幅度构造形态，并在后续通过多次迭代修正速度模型过程中，不断调整和增加速度趋势面，逐步刻画微观小幅度构造形态，最终获得高精度叠前深度偏移速度模型，满足更复杂偏移算法对偏移成像的需求。法对偏移成像的需求。法对偏移成像的需求。


技术研发人员：杜磊 李国发 饶莹 詹仕凡 杨娜霞
受保护的技术使用者：中国石油大学（北京）
技术研发日：2023.03.27
技术公布日：2023/7/18