CN113009577B - 一种低幅构造校正方法、装置、计算机设备及存储介质 - Google Patents
一种低幅构造校正方法、装置、计算机设备及存储介质 Download PDFInfo
- Publication number
- CN113009577B CN113009577B CN202110309456.6A CN202110309456A CN113009577B CN 113009577 B CN113009577 B CN 113009577B CN 202110309456 A CN202110309456 A CN 202110309456A CN 113009577 B CN113009577 B CN 113009577B
- Authority
- CN
- China
- Prior art keywords
- amplitude
- low
- abnormal body
- speed abnormal
- mean square
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012937 correction Methods 0.000 title claims abstract description 45
- 230000002159 abnormal effect Effects 0.000 claims abstract description 113
- 238000011161 development Methods 0.000 claims abstract description 16
- 238000010586 diagram Methods 0.000 claims description 34
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 239000002609 medium Substances 0.000 description 20
- 235000019738 Limestone Nutrition 0.000 description 14
- 239000006028 limestone Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 239000012120 mounting media Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 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
-
- 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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明实施例公开了一种低幅构造校正方法、装置、计算机设备及存储介质。该方法包括：获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；基于真实地层岩性组合，通过岩性替换得到高速异常体层段的厚度渐变的岩性组合模型；基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；提取高速异常体层段的均方根振幅平面图，并根据均方根振幅平面图和定量关系确定高速异常体层段的平面厚度图；根据平面厚度图确定平面时间误差网格图；根据平面时间误差网格图对低幅构造进行校正。本发明实施例所提供的技术方案，实现了高效精确的对低幅构造进行校正，从而降低了上覆发育薄层高速异常岩性体对低幅构造的影响。
Description
技术领域
本发明实施例涉及地震勘测技术领域，尤其涉及一种低幅构造校正方法、装置、计算机设备及存储介质。
背景技术
低幅构造一般是指构造幅度小于20m的构造。由于其本身构造幅度低，当其上覆发育高速异常岩性体并且横向厚度存在变化时，会造成目的层低幅构造形态发生变化，从而影响目的层资源量的评估和油藏储量规模的落实。如果高速异常体厚度较大，可以对其顶底地震反射作精细解释，刻画出高速异常体平面厚度，消除其对低幅构造的影响。但是如果高速异常体厚度较小且横向存在变化，在常规地震剖面上难以直接刻画，则需要借助其他方法来刻画薄层高速异常体，以校正其对低幅构造的影响。
当前针对受上覆高速异常岩性影响的低幅构造的校正方法包括：网格层析速度建模、地震反演和地震正演等。其中，网格层析速度建模方法为在存在薄层高速异常体影响的地层加密网格分布，从而实现根据地质条件划分网格的策略，进而通过解析矩阵方程实现速度模型的迭代更新，以精细刻画出高速异常体的速度，再用迭代更新后的速度场进行变速成图，以校正低幅构造，但是网格划分精度很难达到薄层的要求，且计算量较大。地震反演方法为通过宽频反演、地质统计学反演以及相控反演等不同的反演方法，反演得到精细的地层速度体，最终用反演得到的速度进行变速成图，以校正低幅构造，但是反演的结果具有多解性，会受到低频速度模型的制约，且反演的精度受地震资料分辨率的限制，难以精细刻画出薄层速度异常。地震正演方法为基于理论模型正演，定性或半定量分析高速异常体对下覆低幅构造的影响，从而得到经验公式，以校正低幅构造，但是由于是基于理论模型，通常与实际地层岩性组合存在差异，得到的经验公式校正低幅构造的误差较大。
发明内容
本发明实施例提供一种低幅构造校正方法、装置、计算机设备及存储介质，以高效精确的对低幅构造进行校正，从而降低上覆发育薄层高速异常岩性体对低幅构造的影响。
第一方面，本发明实施例提供了一种低幅构造校正方法，该方法包括：
获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
基于所述真实地层岩性组合，通过岩性替换得到所述高速异常体层段的厚度渐变的岩性组合模型；
基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
提取所述高速异常体层段的均方根振幅平面图，并根据所述均方根振幅平面图和所述定量关系确定所述高速异常体层段的平面厚度图；
根据所述平面厚度图确定平面时间误差网格图；
根据所述平面时间误差网格图对所述低幅构造进行校正。
第二方面，本发明实施例还提供了一种低幅构造校正装置，该装置包括：
岩性组合获取模块，用于获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
岩性替换模块，用于基于所述真实地层岩性组合，通过岩性替换得到所述高速异常体层段的厚度渐变的岩性组合模型；
定量关系确定模块，用于基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
厚度图确定模块，用于提取所述高速异常体层段的均方根振幅平面图，并根据所述均方根振幅平面图和所述定量关系确定所述高速异常体层段的平面厚度图；
时间误差网格图确定模块，用于根据所述平面厚度图确定平面时间误差网格图；
校正模块，用于根据所述平面时间误差网格图对所述低幅构造进行校正。
第三方面，本发明实施例还提供了一种计算机设备，该计算机设备包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现本发明任意实施例所提供的低幅构造校正方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现本发明任意实施例所提供的低幅构造校正方法。
本发明实施例提供了一种低幅构造校正方法，首先获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，通过岩性变换，逐渐改变该高速异常体层段的厚度，得到厚度渐变的岩性组合模型，然后基于该岩性组合模型进行地震正演，以得到均方根振幅与高速异常体厚度之间的定量关系，接着提取该高度异常体层段的均方根振幅平面图，从而对应该定量关系确定该高速异常体层段的平面厚度图，再根据该平面厚度图确定平面时间误差网格图，从而根据该平面时间误差网格图对低幅构造进行校正。本发明实施例所提供的低幅构造校正方法，通过基于真实地层岩性组合进行岩性替换后作地震正演，并确定均方根振幅与高速异常体厚度之间的定量关系，进而确定平面时间误差网格对低幅构造进行校正，使得对低幅构造的校正更加高效精确，从而降低了其上覆发育薄层高速异常岩性体对低幅构造的影响，进而提高了目的层资源量评估和油藏储量规模落实的准确性。
附图说明
图1为本发明实施例一提供的低幅构造校正方法的流程图；
图2为本发明实施例二提供的低幅构造校正装置的结构示意图；
图3为本发明实施例三提供的计算机设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
在更加详细地讨论示例性实施例之前应当提到的是，一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各步骤描述成顺序的处理，但是其中的许多步骤可以被并行地、并发地或者同时实施。此外，各步骤的顺序可以被重新安排。当其操作完成时所述处理可以被终止，但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
实施例一
图1为本发明实施例一提供的低幅构造校正方法的流程图。本实施例可适用于对受上覆薄层高速异常体影响的低幅构造进行校正的情况，该方法可以由本发明实施例所提供的低幅构造校正装置来执行，该装置可以由硬件和/或软件的方式来实现，一般可集成于计算机设备中。如图1所示，具体包括如下步骤：
S11、获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合。
具体的，首先获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，可选的，获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，包括：根据已钻井的测井资料，获取真实地层岩性组合，测井资料包括泥岩、砂岩和高速异常体层段的密度和声波速度。其中，该已钻井位于待校正的低幅构造的周边。根据实际的测井资料，即可绘制出高速异常体层段以及其附近层段的真实地层岩性组合。
S12、基于真实地层岩性组合，通过岩性替换得到高速异常体层段的厚度渐变的岩性组合模型。
具体的，基于真实地层岩性组合，通过岩性替换，逐渐改变薄层高速异常体的厚度，从而得到高速异常体层段厚度渐变的岩性组合模型，具体可以是在被砂岩相隔的上下两层相邻的高速异常体层段中，从左至右的逐渐向上增加下层的高速异常体厚度，使得下层的高速异常体逐渐靠近上层的高速异常体，并且下层的高速异常体厚度可以呈近线性增长。
S13、基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系。
在获得岩性组合模型之后，可以基于该岩性组合模型进行地震正演，以区别于传统在理论模型上进行的地震正演。通过地震正演，即可获得均方根振幅属性与高速异常体厚度属性之间的定量关系。可选的，在基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系之前，还包括：根据井旁道地震数据，获取高速异常体层段的地震主频；相应的，基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系，包括：基于岩性组合模型和主频为地震主频的雷克子波进行地震正演，得到岩性组合模型的合成地震反射记录；从合成地震反射记录中提取高速异常体层段的反射均方根振幅；根据反射均方根振幅确定定量关系。具体的，通过分析井旁道地震数据，可以获取高速异常体层段的地震主频，然后可以基于得到的岩性组合模型和主频为该地震主频的雷克子波进行地震正演。通过地震正演可以得到岩性组合模型的合成地震反射记录，然后可以在该合成地震反射记录上提取高速异常体层段上各厚度对应的反射均方根振幅，即得到了均方根振幅与高速异常体厚度之间的对应关系，接着可以根据该对应关系进行拟合，从而获得均方根振幅与高速异常体厚度之间的函数关系式，即所需的定量关系。
S14、提取高速异常体层段的均方根振幅平面图，并根据均方根振幅平面图和定量关系确定高速异常体层段的平面厚度图。
然后即可通过提取高速异常体层段的均方根振幅平面图并代入上述确定的定量关系来确定对应的平面厚度图。可选的，提取高速异常体层段的均方根振幅平面图，包括：在实际地震数据体上对高速异常体层段的反射同相轴进行精细解释，以提取均方根振幅平面图，并对均方根振幅平面图的属性值进行尺度标定。具体的，在提取到均方根振幅平面图之后，可以将均方根振幅平面图的属性值的范围标定到合成地震反射记录上的均方根振幅属性值的范围内，以便于对应确定厚度。然后即可基于定量关系和标定后的均方根振幅平面图计算得到高速异常体层段的平面厚度图。
S15、根据平面厚度图确定平面时间误差网格图。
具体的，在获得平面厚度图之后，可以根据该平面厚度图计算高速异常体层段引起的时间误差，从而得到平面时间误差网格图。
S16、根据平面时间误差网格图对低幅构造进行校正。
可选的，根据平面时间误差网格图对低幅构造进行校正，包括：根据平面时间误差网格图对低幅构造的时间网格进行校正，以得到校正后的时间构造图；根据已钻井的时深关系，对校正后的时间构造图进行单井时深转换，以得到校正后的深度构造图。具体的，可以使用平面时间误差网格图对低幅构造的时间网格进行精细校正，再通过单井时深转换即可得到校正后低幅构造的深度构造图。
示例性的，选定珠江口盆地某构造为应用实例区进行说明。该构造的主要目的层H1层构造为典型的低幅构造，上覆发育薄层高速灰岩影响其构造落实。通过地震频谱分析，发育薄层灰岩层段的地震主频约为40赫兹。在该构造的北高点钻探一口井，通过已钻井得到灰岩发育层段的真实地层岩性组合，并获得其中的砂岩的声波速度为2800米每秒，密度为2.06千克每立方米，泥岩的声波速度为3000米每秒，密度为2.29千克每立方米，灰岩的声波速度为4500米每秒，密度为2.7千克每立方米。然后通过岩性替换，逐渐改变灰岩的厚度，得到岩性替换后的岩性组合模型。基于得到的岩性组合模型和主频为40赫兹的雷克子波，通过地震正演，得到该岩性组合模型的合成地震反射记录。在该合成地震反射记录上，提取灰岩的反射均方根振幅，从而得到均方根振幅与灰岩厚度之间的对应关系，然后通过拟合可得到二者之间的函数关系式为y＝0.0181*x-0.0697，其中，y表示均方根振幅大小，x表示灰岩厚度。在实际地震数据体上对薄层灰岩的反射同相轴进行精细解释，提取其均方根振幅，并对提取的均方根振幅进行尺度标定。基于拟合的函数关系式和标定后的均方根振幅平面图，计算薄层灰岩的平面厚度图。根据已钻井的真实地层岩性组合和薄层灰岩的平面厚度图，计算灰岩引起的时间误差，得到平面时间误差网格图。使用平面时间误差网格对目的层低幅构造的时间网格进行校正，得到校正后的时间构造图。最后根据已钻井的时深关系，对校正后的时间构造图进行单井时深转换，即可得到校正后的低幅构造深度构造图，从而完成对低幅构造的校正过程，使用该深度构造图进行目的层资源量的评估和油藏储量规模的落实等将获得更加准确的结果。
本发明实施例所提供的技术方案，首先获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，通过岩性变换，逐渐改变该高速异常体层段的厚度，得到厚度渐变的岩性组合模型，然后基于该岩性组合模型进行地震正演，以得到均方根振幅与高速异常体厚度之间的定量关系，接着提取该高度异常体层段的均方根振幅平面图，从而对应该定量关系确定该高速异常体层段的平面厚度图，再根据该平面厚度图确定平面时间误差网格图，从而根据该平面时间误差网格图对低幅构造进行校正。通过基于真实地层岩性组合进行岩性替换后作地震正演，并确定均方根振幅与高速异常体厚度之间的定量关系，进而确定平面时间误差网格对低幅构造进行校正，使得对低幅构造的校正更加高效精确，从而降低了其上覆发育薄层高速异常岩性体对低幅构造的影响，进而提高了目的层资源量评估和油藏储量规模落实的准确性。
实施例二
图2为本发明实施例二提供的低幅构造校正装置的结构示意图，该装置可以由硬件和/或软件的方式来实现，一般可集成于计算机设备中，用于执行本发明任意实施例所提供的低幅构造校正方法。如图2所示，该装置包括：
岩性组合获取模块21，用于获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
岩性替换模块22，用于基于真实地层岩性组合，通过岩性替换得到高速异常体层段的厚度渐变的岩性组合模型；
定量关系确定模块23，用于基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
厚度图确定模块24，用于提取高速异常体层段的均方根振幅平面图，并根据均方根振幅平面图和定量关系确定高速异常体层段的平面厚度图；
时间误差网格图确定模块25，用于根据平面厚度图确定平面时间误差网格图；
校正模块26，用于根据平面时间误差网格图对低幅构造进行校正。
本发明实施例所提供的技术方案，首先获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，通过岩性变换，逐渐改变该高速异常体层段的厚度，得到厚度渐变的岩性组合模型，然后基于该岩性组合模型进行地震正演，以得到均方根振幅与高速异常体厚度之间的定量关系，接着提取该高度异常体层段的均方根振幅平面图，从而对应该定量关系确定该高速异常体层段的平面厚度图，再根据该平面厚度图确定平面时间误差网格图，从而根据该平面时间误差网格图对低幅构造进行校正。通过基于真实地层岩性组合进行岩性替换后作地震正演，并确定均方根振幅与高速异常体厚度之间的定量关系，进而确定平面时间误差网格对低幅构造进行校正，使得对低幅构造的校正更加高效精确，从而降低了其上覆发育薄层高速异常岩性体对低幅构造的影响，进而提高了目的层资源量评估和油藏储量规模落实的准确性。
在上述技术方案的基础上，可选的，该低幅构造校正装置，还包括：
地震主频获取模块，用于在基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系之前，根据井旁道地震数据，获取高速异常体层段的地震主频；
相应的，定量关系确定模块，包括：
地震正演单元，用于基于岩性组合模型和主频为地震主频的雷克子波进行地震正演，得到岩性组合模型的合成地震反射记录；
反射均方根振幅提取单元，用于从合成地震反射记录中提取高速异常体层段的反射均方根振幅；
定量关系确定单元，用于根据反射均方根振幅确定定量关系。
在上述技术方案的基础上，可选的，厚度图确定模块24具体用于：
在实际地震数据体上对高速异常体层段的反射同相轴进行精细解释，以提取均方根振幅平面图，并对均方根振幅平面图的属性值进行尺度标定。
在上述技术方案的基础上，可选的，校正模块26，包括：
时间构造图确定单元，用于根据平面时间误差网格图对低幅构造的时间网格进行校正，以得到校正后的时间构造图；
深度构造图确定单元，用于根据已钻井的时深关系，对校正后的时间构造图进行单井时深转换，以得到校正后的深度构造图。
在上述技术方案的基础上，可选的，岩性组合获取模块21具体用于：
根据已钻井的测井资料，获取真实地层岩性组合，测井资料包括泥岩、砂岩和高速异常体层段的密度和声波速度。
本发明实施例所提供的低幅构造校正装置可执行本发明任意实施例所提供的低幅构造校正方法，具备执行方法相应的功能模块和有益效果。
值得注意的是，在上述低幅构造校正装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
实施例三
图3为本发明实施例三提供的计算机设备的结构示意图，示出了适于用来实现本发明实施方式的示例性计算机设备的框图。图3显示的计算机设备仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。如图3所示，该计算机设备包括处理器31、存储器32、输入装置33及输出装置34；计算机设备中处理器31的数量可以是一个或多个，图3中以一个处理器31为例，计算机设备中的处理器31、存储器32、输入装置33及输出装置34可以通过总线或其他方式连接，图3中以通过总线连接为例。
存储器32作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的低幅构造校正方法对应的程序指令/模块(例如，低幅构造校正装置中的岩性组合获取模块21、岩性替换模块22、定量关系确定模块23、厚度图确定模块24、时间误差网格图确定模块25及校正模块26)。处理器31通过运行存储在存储器32中的软件程序、指令以及模块，从而执行计算机设备的各种功能应用以及数据处理，即实现上述的低幅构造校正方法。
存储器32可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据计算机设备的使用所创建的数据等。此外，存储器32可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器32可进一步包括相对于处理器31远程设置的存储器，这些远程存储器可以通过网络连接至计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置33可用于获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，以及产生与计算机设备的用户设置和功能控制有关的键信号输入等。输出装置34可包括显示屏等设备，可用于向用户展示校正过程中所使用的各类数据以及校正后的数据等等。
实施例四
本发明实施例四还提供一种包含计算机可执行指令的存储介质，该计算机可执行指令在由计算机处理器执行时用于执行一种低幅构造校正方法，该方法包括：
获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
基于真实地层岩性组合，通过岩性替换得到高速异常体层段的厚度渐变的岩性组合模型；
基于岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
提取高速异常体层段的均方根振幅平面图，并根据均方根振幅平面图和定量关系确定高速异常体层段的平面厚度图；
根据平面厚度图确定平面时间误差网格图；
根据平面时间误差网格图对低幅构造进行校正。
存储介质可以是任何的各种类型的存储器设备或存储设备。术语“存储介质”旨在包括：安装介质，例如CD-ROM、软盘或磁带装置；计算机系统存储器或随机存取存储器，诸如DRAM、DDR RAM、SRAM、EDO RAM，兰巴斯(Rambus)RAM等；非易失性存储器，诸如闪存、磁介质(例如硬盘或光存储)；寄存器或其它相似类型的存储器元件等。存储介质可以还包括其它类型的存储器或其组合。另外，存储介质可以位于程序在其中被执行的计算机系统中，或者可以位于不同的第二计算机系统中，第二计算机系统通过网络(诸如因特网)连接到计算机系统。第二计算机系统可以提供程序指令给计算机用于执行。术语“存储介质”可以包括可以驻留在不同位置中(例如在通过网络连接的不同计算机系统中)的两个或更多存储介质。存储介质可以存储可由一个或多个处理器执行的程序指令(例如具体实现为计算机程序)。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的低幅构造校正方法中的相关操作。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、电线、光缆、RF等等，或者上述的任意合适的组合。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (8)
1.一种低幅构造校正方法，其特征在于，包括：
获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
基于所述真实地层岩性组合，通过岩性替换得到所述高速异常体层段的厚度渐变的岩性组合模型；
基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
提取所述高速异常体层段的均方根振幅平面图，并根据所述均方根振幅平面图和所述定量关系确定所述高速异常体层段的平面厚度图；
根据所述平面厚度图确定平面时间误差网格图；
根据所述平面时间误差网格图对所述低幅构造进行校正；
在所述基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系之前，还包括：
根据井旁道地震数据，获取所述高速异常体层段的地震主频；
相应的，所述基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系，包括：
基于所述岩性组合模型和主频为所述地震主频的雷克子波进行地震正演，得到所述岩性组合模型的合成地震反射记录；
从所述合成地震反射记录中提取所述高速异常体层段的反射均方根振幅；
根据所述反射均方根振幅确定所述定量关系。
2.根据权利要求1所述的低幅构造校正方法，其特征在于，所述提取所述高速异常体层段的均方根振幅平面图，包括：
在实际地震数据体上对所述高速异常体层段的反射同相轴进行精细解释，以提取所述均方根振幅平面图，并对所述均方根振幅平面图的属性值进行尺度标定。
3.根据权利要求1所述的低幅构造校正方法，其特征在于，所述根据所述平面时间误差网格图对所述低幅构造进行校正，包括：
根据所述平面时间误差网格图对所述低幅构造的时间网格进行校正，以得到校正后的时间构造图；
根据已钻井的时深关系，对所述校正后的时间构造图进行单井时深转换，以得到校正后的深度构造图。
4.根据权利要求1所述的低幅构造校正方法，其特征在于，所述获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合，包括：
根据已钻井的测井资料，获取所述真实地层岩性组合，所述测井资料包括泥岩、砂岩和所述高速异常体层段的密度和声波速度。
5.一种低幅构造校正装置，其特征在于，包括：
岩性组合获取模块，用于获取低幅构造的上覆发育薄层高速异常体层段的真实地层岩性组合；
岩性替换模块，用于基于所述真实地层岩性组合，通过岩性替换得到所述高速异常体层段的厚度渐变的岩性组合模型；
定量关系确定模块，用于基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系；
厚度图确定模块，用于提取所述高速异常体层段的均方根振幅平面图，并根据所述均方根振幅平面图和所述定量关系确定所述高速异常体层段的平面厚度图；
时间误差网格图确定模块，用于根据所述平面厚度图确定平面时间误差网格图；
校正模块，用于根据所述平面时间误差网格图对所述低幅构造进行校正；
所述装置还包括：
地震主频获取模块，用于在所述基于所述岩性组合模型进行地震正演，得到均方根振幅与高速异常体厚度之间的定量关系之前，根据井旁道地震数据，获取所述高速异常体层段的地震主频；
相应的，所述定量关系确定模块，包括：
地震正演单元，用于基于所述岩性组合模型和主频为所述地震主频的雷克子波进行地震正演，得到所述岩性组合模型的合成地震反射记录；
反射均方根振幅提取单元，用于从所述合成地震反射记录中提取所述高速异常体层段的反射均方根振幅；
定量关系确定单元，用于根据所述反射均方根振幅确定所述定量关系。
6.根据权利要求5所述的低幅构造校正装置，其特征在于，所述校正模块，包括：
时间构造图确定单元，用于根据所述平面时间误差网格图对所述低幅构造的时间网格进行校正，以得到校正后的时间构造图；
深度构造图确定单元，用于根据已钻井的时深关系，对所述校正后的时间构造图进行单井时深转换，以得到校正后的深度构造图。
7.一种计算机设备，其特征在于，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-4中任一所述的低幅构造校正方法。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-4中任一所述的低幅构造校正方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110309456.6A CN113009577B (zh) | 2021-03-23 | 2021-03-23 | 一种低幅构造校正方法、装置、计算机设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110309456.6A CN113009577B (zh) | 2021-03-23 | 2021-03-23 | 一种低幅构造校正方法、装置、计算机设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113009577A CN113009577A (zh) | 2021-06-22 |
| CN113009577B true CN113009577B (zh) | 2023-07-21 |
Family
ID=76405566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110309456.6A Active CN113009577B (zh) | 2021-03-23 | 2021-03-23 | 一种低幅构造校正方法、装置、计算机设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113009577B (zh) |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2191414C1 (ru) * | 2001-12-28 | 2002-10-20 | Копилевич Ефим Абрамович | Способ геофизической разведки для выявления малоамплитудных тектонических нарушений нефтегазопродуктивных горных пород |
| CN102129084B (zh) * | 2010-12-17 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
| CA2904008C (en) * | 2013-03-15 | 2020-10-27 | Schlumberger Canada Limited | Methods of characterizing earth formations using physiochemical model |
| CN105911592B (zh) * | 2016-06-30 | 2017-11-10 | 中国石油化工股份有限公司 | 一种实钻约束的三维地震资料长波长静校正方法 |
| CN107688200B (zh) * | 2017-04-12 | 2019-09-03 | 恒泰艾普集团股份有限公司 | 针对低幅度构造畸变的速度校正方法 |
| CN111208558B (zh) * | 2018-11-22 | 2022-05-06 | 中国石油天然气股份有限公司 | 超深低幅度三维地质构造的建立方法及装置 |
| CN111913218B (zh) * | 2019-05-09 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种基于多尺度混合反演的薄层反演方法、设备以及系统 |
-
2021
- 2021-03-23 CN CN202110309456.6A patent/CN113009577B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113009577A (zh) | 2021-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Holcombe et al. | Three-dimensional terrain corrections in resistivity surveys | |
| US20210003731A1 (en) | Method for determining favorable time window of infill well in unconventional oil and gas reservoir | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| CN105158796B (zh) | 确定toc含量的方法和装置 | |
| US9366771B2 (en) | System and method for migration velocity modeling | |
| CN107587871A (zh) | 确定水平裂缝宽度的方法及装置 | |
| CN112965103A (zh) | 一种多重孔隙储层叠前地震概率化多道反演方法 | |
| CN105275459A (zh) | 一种确定页岩地层可动水体积含量的方法 | |
| CN113009577B (zh) | 一种低幅构造校正方法、装置、计算机设备及存储介质 | |
| CN111239822B (zh) | 一种消除井震闭合差的断层下速度建模方法及处理终端 | |
| CN113945978A (zh) | 一种深度域地震数据优化方法及装置 | |
| CN109188522B (zh) | 速度场构建方法及装置 | |
| CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| CN107991705B (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
| CN111257969B (zh) | 一种断层控制下高精度速度建模方法及处理终端 | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
| CN113820745A (zh) | 地震速度建模方法、装置、电子设备及介质 | |
| CN113296157B (zh) | 一种利用广义高斯分布进行储层的预测方法及装置 | |
| CN115903026B (zh) | 复合砂体构型解析方法、设备及介质 | |
| CN113219533B (zh) | 一种叠前时间偏移速度建模方法、装置、介质及电子设备 | |
| CN113219534B (zh) | 一种叠前深度偏移速度质控方法、装置、介质及电子设备 | |
| CN114706126B (zh) | 一种地质认识约束下的地震地质体定量雕刻方法及系统 | |
| CN118330739B (zh) | 基于井震联合约束的深度域数据体获取方法、系统及存储介质 | |
| US20240151868A1 (en) | Prediction of Permeability Distribution of Deep Buried Sandstones Using Machine Learning | |
| CN111198400B (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 |