CN112859167B - 畸变地质体速度场的校正方法及装置 - Google Patents
畸变地质体速度场的校正方法及装置 Download PDFInfo
- Publication number
- CN112859167B CN112859167B CN201911190151.7A CN201911190151A CN112859167B CN 112859167 B CN112859167 B CN 112859167B CN 201911190151 A CN201911190151 A CN 201911190151A CN 112859167 B CN112859167 B CN 112859167B
- Authority
- CN
- China
- Prior art keywords
- geologic body
- corrected
- control layer
- geologic
- distorted
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 120
- 238000000034 method Methods 0.000 title claims abstract description 91
- 238000009499 grossing Methods 0.000 claims abstract description 100
- 239000010410 layer Substances 0.000 claims description 161
- 238000004590 computer program Methods 0.000 claims description 17
- 239000011229 interlayer Substances 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 abstract description 67
- 230000005012 migration Effects 0.000 abstract description 32
- 238000013508 migration Methods 0.000 abstract description 32
- 238000012545 processing Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 238000003908 quality control method Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 238000012512 characterization method Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种畸变地质体速度场的校正方法及装置，该方法包括：根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上下地质体；对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。本方法可以同时解决复杂地质体下伏地层深度偏移成像所有畸变问题，提高复杂地质条件深度偏移成像效率。
Description
技术领域
本发明涉及石油勘探领域，尤其是地球物理勘探领域，具体涉及一种畸变地质体速度场的校正方法及装置。
背景技术
在深度偏移处理中，本领域技术人员经常会被某些复杂地质体引起的下伏地层成像畸变所困扰。复杂地质体是指成像速度刻画意义上复杂的地质体，其速度往往在横向变化剧烈并且边界难以刻画，造成其下伏平缓地层发生一致性成像扭曲假象，即成像畸变，比如火成岩造成的其下平层成像的扭曲或错断(如图1所示)，盐丘造成的其下地层成像形态的扭曲等。这些扭曲靠现有的层析迭代技术很难解决。常规的解决方法是根据畸变地质体，拾取期望得到的目标层位，然后根据目标层位反算出消除畸变的地质体的速度场。常规的方法能较好解决成像畸变，但当工区范围特别大，成像畸变特别多时，对单个地质体拾取工作量很大，并且经常还会漏掉一些需要校正的地方。
发明内容
针对现有技术中的问题，本发明的畸变地质体速度场的校正方法和装置，可以提供一种可一次解决复杂地质体下伏地层深度偏移成像中发生多个畸变的问题，从而提高复杂地质条件深度偏移成像效率。
为解决上述技术问题，本发明提供以下技术方案：
第一方面，本发明提供一种畸变地质体速度场的校正方法，包括：
根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上、下地质体；
对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
一实施例中，所述根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正，包括：
根据平滑结果计算平滑后的待校正地质体与其上的地质体之间的距离以及平滑后的待校正地质体与其下的地质体之间的距离；
根据待校正地质体与其上的地质体之间的距离、待校正地质体与其下的地质体之间的距离、所述平滑后的待校正地质体与其上的地质体之间的距离以及所述平滑后的待校正地质体与其下的地质体之间的距离生成速度场畸变校正因子；
根据所述速度场畸变校正因子对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
一实施例中，所述对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果，包括：
利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数；
根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
一实施例中，畸变地质体速度场的校正方法还包括：根据地质体速度在横向上的变化情况以及所述地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
第二方面，本发明提供一种畸变地质体速度场的校正装置，该装置包括：
地质体拾取单元，用于根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上、下地质体；
地质体平滑单元，用于对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
速度场校正单元，用于根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
一实施例中，所述速度场校正单元包括：
距离计算模块，用于根据平滑结果计算平滑后的待校正地质体与其上的地质体之间的距离以及平滑后的待校正地质体与其下的地质体之间的距离；
因子生成模块，用于根据待校正地质体与其上的地质体之间的距离、待校正地质体与其下的地质体之间的距离、所述平滑后的待校正地质体与其上的地质体之间的距离以及所述平滑后的待校正地质体与其下的地质体之间的距离生成速度场畸变校正因子；
速度场校正模块，用于根据所述速度场畸变校正因子对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
一实施例中，所述地质体平滑单元包括：
参数计算模块，用于利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数；
地质体平滑模块，用于根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
一实施例中，畸变地质体速度场的校正装置还包括畸变判断单元，用于根据地质体速度在横向上的变化情况以及所述地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现畸变地质体速度场的校正方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现畸变地质体速度场的校正方法的步骤。
本发明实施例提供的畸变地质体速度场的校正方法及装置，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的实施例中的畸变地质体速度场的校正方法的畸变校正模型的区域识别示意图；
图2为本发明的实施例中的畸变地质体速度场的校正方法的流程示意图；
图3为本发明的实施例中步骤300的流程示意图；
图4为本发明的实施例中步骤200的流程示意图；
图5为本发明的实施例中的畸变地质体速度场的校正方法的畸变校正模型的另一流程示意图；
图6为本发明的具体应用实例中畸变地质体速度场的校正方法的流程示意图；
图7为本发明的具体应用实例中成像畸变判断示意图；
图8为本发明的具体应用实例中畸变校正模型的层位设置示意图；
图9为本发明的具体应用实例中畸变校正模型层位拾取方案图；
图10为本发明的具体应用实例中畸变校正模型的数学抽象示意图；
图11为本发明的具体应用实例中地质体控层组H12深度以及速度校正前后对比图；
图12为本发明的具体应用实例中地质体下控层H13深度以及速度校正前后对比图；
图13A为本发明的具体应用实例中过Z InLine线速度场校正前速度场示意图；
图13B为本发明的具体应用实例中过Z InLine线速度场校正后速度场示意图；
图14A为本发明的具体应用实例中过Z InLine线深度成像剖面及延迟谱图；
图14B为本发明的具体应用实例中过Z InLine线深度成像剖面及延迟谱校正效果图；
图15为本发明的实施例中畸变地质体速度场的校正装置的结构框图；
图16为本发明的实施例中速度场校正单元结构框图；
图17为本发明的实施例中地质体平滑单元结构框图；
图18为本发明的实施例中畸变地质体速度场的校正装置的另一结构框图；
图19为本发明的实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
针对现有技术速度刻画与成像过程中，复杂地质体常常造成其下伏平缓地层发生一致性成像扭曲假象，即成像畸变。基于此，本发明的实施例提供一种畸变地质体速度场的校正方法的具体实施方式，参见图2，该方法具体包括如下内容：
步骤100：根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上、下地质体。
可以理解的是，振幅反映反射的强度，是波阻抗的间接表示，是能量在传播过程中质点振动的幅度(偏离平衡位置)。同相轴是地震记录上各道振动相位相同的极值(俗称波峰或波谷)的连线。波形是指地震波的形态，包括质点的振动方向、传播空间、传播路径以及各个时刻各指点偏离平衡位置等等。
拾取层位(待校正地质体以及其上、下地质体之间的界限)是以振幅、同相轴连续性以及波形这三个基本准则为基础所进行的，优选地，在二维或者三维地震剖面上利用波形相似性对地层连续反射同相轴进行手动或者自动追踪得到层位线(地层界面)，再对所有的层位线进行插值从而形成层位面，最终得到各自的地质体。
步骤200：对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果。
可以理解的是，地下地质体由于岩性或者其他性质发生变化会发生畸变，如：火成岩造成的其下平层成像的扭曲或错断，盐丘造成的其下地层成像形态的扭曲等。如果想得到准确的速度场，需要对畸变的部分进行校正，可以理解的是，对待校正地质体需要进行较小尺度的平滑，对其下的地质体需要进行较大尺度的平滑。
步骤300：根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
由物理学原理可知，地震波在地层中垂向传播时存在如下关系：
h＝v·t (1)
其中，h为地层的厚度，v为地层的垂向速度，t为地层垂向旅行时间。由公式(1)可以看出，在准确得出了地质体距离的基础上，才可以得出准确的速度场。
本发明实施例提供畸变地质体速度场的校正方法，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。
一实施例中，参见图3所示的步骤300的具体流程，步骤300具体包括：
步骤301：根据平滑结果计算平滑后的待校正地质体与其上的地质体之间的距离以及平滑后的待校正地质体与其下的地质体之间的距离。
步骤302：根据待校正地质体与其上的地质体之间的距离、待校正地质体与其下的地质体之间的距离、所述平滑后的待校正地质体与其上的地质体之间的距离以及所述平滑后的待校正地质体与其下的地质体之间的距离生成速度场畸变校正因子。
步骤303：根据所述速度场畸变校正因子对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
在步骤301至步骤303中，在深度偏移速度刻画与成像过程中，因为地层满足等t(即，地层垂向旅行时间)的原则，所以对于待校正地质体上控层(待校正地质体之上的地质体)与待校正地质体间以及待校正地质体与待校正地质体下控层(待校正地质体之下的地质体)间两套地层分别存在如下关系：
其中，hA为待校正地质体到地质体上控层的距离，hAS为待校正地质体平滑后对应的距离，vA为待校正地质体与地质体上控层层间速度，vAS为消除待校正地质体畸变的层间速度，hB为地质体下控层到地质体上控层的距离，hBS为地质体下控层平滑后对应的距离，vB为地质体下控层与待校正地质体层间速度，vBS消除地质体下控层畸变的层间速度。对式(2)整理，可得：
其中，αA为待校正地质体的畸变校正因子，αB为地质体下控层的畸变校正因子。至此，推导出待校正地质体设置一个层位时各层位的畸变校正因子，进而求得各层位消除畸变的速度场。
一实施例中，参见图4，步骤200具体包括：
步骤201：利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数。
步骤202：根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
需要指出的是，使用消除地质体畸变的层间速度进行层位约束速度建模更新原来的速度模型，接着使用新产生的速度模型进行叠前深度偏移成像，最后质控新偏移结果是否满足畸变校正要求。当不满足时，需要重新修改各层位平滑参数，然后进行后续处理，直到满足畸变校正要求。
一实施例中，参见图5，畸变地质体速度场的校正方法还包括：
步骤400：根据地质体速度在横向上的变化情况以及所述地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
可以理解的是，在步骤100之前，需要进行复杂地质体下伏地层深度偏移成像畸变判断。具体地，成像畸变涉及复杂地质体和下伏地层，从地质上，了解本区复杂地质体岩性，分析下伏地层构造力学特征；从处理上，了解复杂地质体速度刻画，分析下伏地层偏移成像特征。
本发明实施例提供的畸变地质体速度场的校正方法，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
为进一步地说明本方案，以下以某油田A区块为例，结合图6来描述畸变地质体速度场的校正方法的具体应用实例，该具体应用实例具体包括如下内容，参见图6。
S0：判断所述地质体是否为畸变地质体。
成像畸变涉及复杂地质体和下伏地层，从地质上，了解本区复杂地质体岩性，分析下伏地层构造力学特征；从处理上，了解复杂地质体速度刻画，分析下伏地层偏移成像特征。具体地，如图7所示，本区复杂反射层发育火成岩，其速度在横向变化剧烈并且边界难以刻画，且下伏地层等幅度一致性扭曲难以从构造力学解释，故判定该剖面发生成像畸变。本发明实施例中畸变校正模型涉及复杂地质体和下伏平缓地层，该畸变符合此模型。
S1：根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上下地质体。
当深度偏移叠加数据体中存在成像畸变时，对其按照畸变校正模型设置层位(如图1及图8所示)。首先，识别畸变校正模型，畸变校正模型围绕复杂地质体建立，复杂地质体是指成像速度刻画意义上复杂的地质体，比如火成岩、盐丘等，复杂地质体之上为未被干扰的上覆地层，之下为成像一致性扭曲区域。然后，设置畸变校正模型的上控层、下控层等，复杂地质体上边界，即复杂地质体与上覆地层间的界面，称为地质体上控层，该层位未发生成像畸变；复杂地质体的下边界，即成像一致性扭曲区域的上界面，称为地质体的下控层，该层位与其下地层发生等幅度的成像畸变，与其上上控层间即复杂地质体发生过渡畸变；在复杂地质体内部可以追踪地质体控层组，用于畸变校正过程中对复杂地质体内幕构造形态精细控制；与地质体的下控层平行，成像一致性扭曲区域中可以追踪出一系列地层畸变一致性层，这些层位均与地质体下控层层间等厚，用于分析检查该区域成像扭曲及校正情况。
因为复杂地质体内幕成像规律性差并且难以解释，地质体控层组设置比较灵活。当不设置时，复杂地质体内幕构造形态均由地质体下控层控制；当设置一个层位时，复杂地质体就被分为上下两部分，可以通过分别控制该层位和地质体下控层实现对复杂地质体内幕构造形态精细控制；当设置更多层位时，依此类推。优选地，设置一个地质体控层组层位。
确定畸变校正模型后，需对地质体上控层、下控层等层位进行拾取。该拾取需在发生成像畸变的深度偏移叠加数据体上进行，通常可以使用自动拾取方法，效率较高。相比之下，常规的方法除了上述拾取外还需要针对单个地质体手工拾取期望得到的目标层位，效率较低且容易发生遗漏。
设置畸变校正模型控制层(待校正地质体)，自动拾取层位。进行畸变校正模型的区域识别、层位设置以及层位拾取，最终确定畸变校正模型层位拾取结果如图9所示，其中H11为地质体上控层，H12为地质体控层组，H13为地质体下控层，H14为地层畸变一致性层。
S2：对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果。
图10为本具体应用实例中畸变校正模型的数学抽象示意图，具体地，对畸变校正模型进行数学抽象，当地质体控层组设置一个层位时，畸变校正因子推导方法可以用图10所示表达。常规的方法需在每个畸变位置(如图10中A和B位置)手工拾取期望得到的目标层位，效率较低且容易发生遗漏，畸变模型校正方法通过对地质体控层组和地质体下控层进行数学运算，同时自动校正层位所有畸变，可达到与手工拾取相当的效果。这里使用的数学运算为一种常用的层位平滑的方法，依据各层位发生畸变的幅度选择合适的平滑参数。通常，对地质体控层组进行较小尺度的平滑，对地质体下控层进行较大尺度的平滑。
S3：根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
当地质体控层组设置一个层位时，分别推导地质体控层组和地质体下控层畸变校正因子。由物理学原理可知，地震波在地层中垂向传播时存在如下关系：
h＝v·t
其中，h为地层的厚度，v为地层的垂向速度，t为地层垂向旅行时。在深度偏移速度刻画与成像过程中，因为地层满足等t的原则，所以对于地质体上控层与地质体控层组间以及地质体控层组与地质体下控层间两套地层分别存在如下关系：
其中，hA为地质体控层组到地质体上控层的距离，hAS为地质体控层组平滑后对应的距离，vA为地质体控层组与地质体上控层层间速度，vAS为消除地质体控层组畸变的层间速度，hB为地质体下控层到地质体上控层的距离，hBS为地质体下控层平滑后对应的距离，vB为地质体下控层与地质体控层组层间速度，vBS消除地质体下控层畸变的层间速度。对上式整理，可得：
其中，αA为地质体控层组的畸变校正因子，αB为地质体下控层的畸变校正因子。至此，推导出地质体控层组设置一个层位时各层位的畸变校正因子，进而求得各层位消除畸变的速度场。接着，使用平滑后的地质体控层组、消除地质体控层组畸变的层间速度、平滑后的地质体下控层、消除地质体下控层畸变的层间速度进行层位约束速度建模更新原来的速度模型，接着使用新产生的速度模型进行叠前深度偏移成像，最后质控新偏移结果是否满足畸变校正要求。当不满足时，重新修改各层位平滑参数，然后进行后续处理，直到满足畸变校正要求。
针对本具体应用实例，对地质体控层组H12进行较小尺度的平滑，然后求取该层位消除畸变的速度场，结果如图11所示，对地质体下控层H13进行较大尺度的平滑，然后求取该层位消除畸变的速度场，结果如图12所示。从图11和图12可见，层位畸变校正后补充了速度中的高波数成分，同时深度层位畸变得到消除。
S4：进行层位约束速度建模。
使用平滑后的地质体控层组、消除地质体控层组畸变的层间速度、平滑后的地质体下控层、消除地质体下控层畸变的层间速度以及与地质体下控层层间厚度和速度均未发生改变的地层畸变一致性层进行层位约束速度建模更新原来的速度模型，结果如图13A及图13B所示。
S5:进行成像及质控。
进行叠前深度偏移成像，并质控新偏移结果。结果如图14A及图14B所示，图中可见，新的成像结果很好消除了下伏平缓地层扭曲畸变，同时保留了复杂地质体内幕构造细节特征，成像更加准确，解释更加合理。
本发明实施例提供的畸变地质体速度场的校正方法，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
基于相似的发明构思，本申请实施例还提供了畸变地质体速度场的校正装置，优选地用于实现上述实施例所描述的方法。由于畸变地质体速度场的校正装置解决问题的原理与畸变地质体速度场的校正方法相似，因此畸变地质体速度场的校正装置的实施可以参见畸变地质体速度场的校正方法实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的系统较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
本发明的实施例提供一种能够实现畸变地质体速度场的校正方法的畸变地质体速度场的校正装置的具体实施方式，参见图15，畸变地质体速度场的校正装置具体包括如下内容：
地质体拾取单元10，用于根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上下地质体；
地质体平滑单元20，用于对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
速度场校正单元30，用于根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
本发明实施例提供畸变地质体速度场的校正装置，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。
一实施例中，参见图16，所述速度场校正单元30包括：
距离计算模块301，用于根据平滑结果计算平滑后的待校正地质体与其上的地质体之间的距离以及平滑后的待校正地质体与其下的地质体之间的距离；
因子生成模块302，用于根据待校正地质体与其上的地质体之间的距离、待校正地质体与其下的地质体之间的距离、所述平滑后的待校正地质体与其上的地质体之间的距离以及所述平滑后的待校正地质体与其下的地质体之间的距离生成速度场畸变校正因子；
速度场校正模块303，用于根据所述速度场畸变校正因子对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
一实施例中，参见图17，所述地质体平滑单元20包括：
参数计算模块201，用于利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数；
地质体平滑模块202，用于根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
一实施例中，参见图18，畸变地质体速度场的校正装置还包括畸变判断单元40，用于根据地质体速度在横向上的变化情况以及所述地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
本发明实施例提供的畸变地质体速度场的校正装置，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
本申请的实施例还提供能够实现上述实施例中的畸变地质体速度场的校正方法中全部步骤的一种电子设备的具体实施方式，参见图19，电子设备具体包括如下内容：
处理器(processor)1201、存储器(memory)1202、通信接口(CommunicationsInterface)1203和总线1204；
其中，处理器1201、存储器1202、通信接口1203通过总线1204完成相互间的通信；通信接口1203用于实现服务器端设备、校正设备以及用户端设备等相关设备之间的信息传输。
处理器1201用于调用存储器1202中的计算机程序，处理器执行计算机程序时实现上述实施例中的畸变地质体速度场的校正方法中的全部步骤，例如，处理器执行计算机程序时实现下述步骤：
步骤100：根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上下地质体；
步骤200：对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
步骤300：根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
本申请实施例中的电子设备，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
本申请的实施例还提供能够实现上述实施例中的畸变地质体速度场的校正方法中全部步骤的一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的畸变地质体速度场的校正方法的全部步骤，例如，处理器执行计算机程序时实现下述步骤：
步骤100：根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上下地质体；
步骤200：对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
步骤300：根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正。
本申请实施例中的计算机可读存储介质，首先在准确拾取待校正地质体以及其上、下地质体的基础上，对待校正地质体以及其下伏的地质体发生畸变的部位进行平滑，并基于平滑结果结合地震波物理学原理自动、快捷的对发生畸变的地质体的速度场进行校正。本方法在从地质和处理角度判断成像畸变的基础上，通过识别发生畸变的地质体，并对畸变层位进行平滑，进一步推导畸变校正因子，求取消除畸变的地质层位的速度场，接着进行速度建模与叠前深度偏移成像，结果满足质控要求后实现对复杂地质体下伏地层深度偏移成像畸变的校正。本发明中的畸变地质体速度场的校正方法，通过对控制层进行数学运算，同时自动校正层位所有畸变，避免了手工拾取期望得到的目标层位的大量工作，效率得到提高。并且可以灵活设置复杂地质体内幕控制层位，实现对其构造形态精细控制，结果更加准确。通过本发明实施例中的畸变地质体速度场的校正方法，取得的畸变校正后的深度偏移成像，可以作为下一步油田勘探研究的重要基础资料，并且更加高效可靠。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (8)
1.一种畸变地质体速度场的校正方法，其特征在于，所述方法包括：
根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上、下地质体，其中，待校正地质体中发生畸变的部分设置为地质体控层组，待校正地质体与待校正地质体其上的地质体的界面称为地质体上控层，待校正地质体与待校正地质体其下的地质体的界面称为地质体下控层；
对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正；
所述根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正，包括：
根据平滑结果计算平滑后的待校正地质体中地质体控层组与地质体上控层之间的距离以及平滑后的待校正地质体中地质体控层组与地质体下控层之间的距离；
根据待校正地质体中地质体控层组与地质体上控层之间的距离、待校正地质体中地质体控层组与地质体下控层之间的距离、所述平滑后的待校正地质体中地质体控层组与地质体上控层之间的距离以及所述平滑后的待校正地质体中地质体控层组与地质体下控层之间的距离生成速度场畸变校正因子；
地层满足等t的原则，t为地层垂向旅行时间，对于待校正地质体上控层与待校正地质体间以及待校正地质体与待校正地质体下控层间两套地层分别存在如下关系：
，其中，hA为待校正地质体中地质体控层组到地质体上控层的距离，hAS为地质体控层组平滑后到地质体上控层的距离，vA为地质体控层组与地质体上控层层间速度，vAS为消除地质体控层组畸变后地质体控层组与地质体上控层的层间速度，hB为地质体下控层到地质体上控层的距离，hBS为地质体下控层平滑后到地质体上控层的距离，vB为地质体下控层与地质体控层组的层间速度，vBS消除地质体下控层畸变后地质体下控层与地质体控层组的层间速度；
对上述关系进行整理，得到：，其中，，αA为地质体控层组的速度场畸变校正因子，αB为地质体下控层的速度场畸变校正因子；
根据所述地质体控层组的速度场畸变校正因子对所述待校正地质体进行校正以及根据所述地质体下控层的速度场畸变校正因子对所述待校正地质体其下的地质体的速度场进行校正。
2.根据权利要求1所述的校正方法，其特征在于，所述对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果，包括：
利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数；
根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
3.根据权利要求1所述的校正方法，其特征在于，还包括：根据地质体速度在横向上的变化情况以及地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
4.一种畸变地质体速度场的校正装置，其特征在于，所述装置包括：
地质体拾取单元，用于根据地震数据中的振幅参数、同相轴连续性以及波形，分别拾取待校正地质体以及所述待校正地质体上、下地质体，其中，待校正地质体中发生畸变的部分设置为地质体控层组，待校正地质体与待校正地质体其上的地质体的界面称为地质体上控层，待校正地质体与待校正地质体其下的地质体的界面称为地质体下控层；
地质体平滑单元，用于对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑，以生成平滑结果；
速度场校正单元，用于根据所述平滑结果对所述待校正地质体以及所述待校正地质体其下的地质体的速度场进行校正；
所述速度场校正单元包括：
距离计算模块，用于根据平滑结果计算平滑后的待校正地质体中地质体控层组与地质体上控层之间的距离以及平滑后的待校正地质体中地质体控层组与地质体下控层之间的距离；
因子生成模块，用于根据待校正地质体中地质体控层组与地质体上控层之间的距离、待校正地质体中地质体控层组与地质体下控层之间的距离、所述平滑后的待校正地质体中地质体控制组与地质体上控层之间的距离以及所述平滑后的待校正地质体中地质体控层组与地质体下控层之间的距离生成速度场畸变校正因子；
地层满足等t的原则，t为地层垂向旅行时间，对于待校正地质体上控层与待校正地质体间以及待校正地质体与待校正地质体下控层间两套地层分别存在如下关系：
，其中，hA为待校正地质体中地质体控层组到地质体上控层的距离，hAS为地质体控层组平滑后到地质体上控层的距离，vA为地质体控层组与地质体上控层层间速度，vAS为消除地质体控层组畸变后地质体控层组与地质体上控层的层间速度，hB为地质体下控层到地质体上控层的距离，hBS为地质体下控层平滑后到地质体上控层的距离，vB为地质体下控层与地质体控层组的层间速度，vBS消除地质体下控层畸变后地质体下控层与地质体控层组的层间速度；
对上述关系进行整理，得到：
，其中，/>，αA为地质体控层组的速度场畸变校正因子，αB为地质体下控层的速度场畸变校正因子；
速度场校正模块，用于根据所述地质体控层组的速度场畸变校正因子对所述待校正地质体进行校正以及根据所述地质体下控层的速度场畸变校正因子对所述待校正地质体其下的地质体的速度场进行校正。
5.根据权利要求4所述的校正装置，其特征在于，所述地质体平滑单元包括：
参数计算模块，用于利用层位平滑方法，分别计算所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分的平滑参数；
地质体平滑模块，用于根据所述平滑参数对所述待校正地质体发生畸变的部分以及其下的地质体发生畸变的部分进行平滑。
6.根据权利要求4所述的校正装置，其特征在于，所述装置还包括：畸变判断单元，用于根据地质体速度在横向上的变化情况以及地质体下伏地层幅度，判断所述地质体是否为畸变地质体。
7.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至3任一项所述畸变地质体速度场的校正方法的步骤。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至3任一项所述畸变地质体速度场的校正方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911190151.7A CN112859167B (zh) | 2019-11-28 | 2019-11-28 | 畸变地质体速度场的校正方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911190151.7A CN112859167B (zh) | 2019-11-28 | 2019-11-28 | 畸变地质体速度场的校正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112859167A CN112859167A (zh) | 2021-05-28 |
| CN112859167B true CN112859167B (zh) | 2023-08-22 |
Family
ID=75995399
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911190151.7A Active CN112859167B (zh) | 2019-11-28 | 2019-11-28 | 畸变地质体速度场的校正方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112859167B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
| CN107688200A (zh) * | 2017-04-12 | 2018-02-13 | 恒泰艾普集团股份有限公司 | 针对低幅度构造畸变的速度校正方法 |
| CN107817525A (zh) * | 2017-09-29 | 2018-03-20 | 中国石油天然气集团公司 | 一种近地表常速填充方法及装置 |
| CN108107473A (zh) * | 2017-12-12 | 2018-06-01 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9435904B2 (en) * | 2013-12-19 | 2016-09-06 | Chevron U.S.A. Inc. | Method of correcting velocity for complex surface topography |
| US10481287B2 (en) * | 2016-08-05 | 2019-11-19 | Saudi Arabian Oil Company | Surface consistent statics solution and amplification correction |
-
2019
- 2019-11-28 CN CN201911190151.7A patent/CN112859167B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
| CN107688200A (zh) * | 2017-04-12 | 2018-02-13 | 恒泰艾普集团股份有限公司 | 针对低幅度构造畸变的速度校正方法 |
| CN107817525A (zh) * | 2017-09-29 | 2018-03-20 | 中国石油天然气集团公司 | 一种近地表常速填充方法及装置 |
| CN108107473A (zh) * | 2017-12-12 | 2018-06-01 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 温铁民 等.基于地质模型的火山岩体速度建模方法.CPS/SEG北京2018国际地球物理会议暨展览.2018,第293-296页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112859167A (zh) | 2021-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10036829B2 (en) | Fault removal in geological models | |
| CN109001813B (zh) | 一种压制多次波的方法、装置及系统 | |
| CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| CN108614293A (zh) | 砂体预测方法及装置 | |
| US10732311B2 (en) | Method for adjusting attributes of marine seismic survey data | |
| CN109696705A (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
| CN112859167B (zh) | 畸变地质体速度场的校正方法及装置 | |
| US20200049844A1 (en) | Computer implemented method for improving a velocity model for seismic imaging | |
| CN112946746B (zh) | 用于提高薄煤层avo反演准确性的方法及装置 | |
| Toomey et al. | Stress arching and its impact on 4D seismic amplitudes and traveltimes | |
| CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
| US10634804B2 (en) | System and method for dip-guided seismic image stacking | |
| CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
| Yilmaz et al. | Interpretive imaging of seismic data | |
| CN113933897B (zh) | 一种基于二维叠前资料的膏岩尖灭线预测方法及系统 | |
| CN103852786B (zh) | 一种应用于陆上地震数据的逆时偏移成像的方法及系统 | |
| CN113109865B (zh) | 混源激发地震波的处理方法及装置 | |
| CN113933893B (zh) | 二维批量叠前反演方法和白云岩储层预测方法及系统 | |
| CN111596348A (zh) | 基于声波时差测井资料的地层层速度校正方法及装置 | |
| CN114002735B (zh) | 弯曲射线叠前时间偏移速度求取方法及装置 | |
| CN111856557A (zh) | 深度域合成地震记录制作方法及装置 | |
| CN112180446B (zh) | 三维地震数据叠前道集叠加方法及装置 | |
| CN114063165B (zh) | 一种三维地震数据的拼接方法及装置 | |
| CN112824940B (zh) | 反射纵波与反射转换横波的振幅匹配方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |