CN106405644A - 裂缝确定方法和装置 - Google Patents
裂缝确定方法和装置 Download PDFInfo
- Publication number
- CN106405644A CN106405644A CN201610772174.9A CN201610772174A CN106405644A CN 106405644 A CN106405644 A CN 106405644A CN 201610772174 A CN201610772174 A CN 201610772174A CN 106405644 A CN106405644 A CN 106405644A
- Authority
- CN
- China
- Prior art keywords
- point
- cdp
- research area
- fracture
- channel set
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 112
- 238000011160 research Methods 0.000 claims abstract description 162
- 238000004458 analytical method Methods 0.000 claims abstract description 35
- 230000008569 process Effects 0.000 claims description 36
- 238000012937 correction Methods 0.000 claims description 34
- 238000009826 distribution Methods 0.000 claims description 28
- 230000008859 change Effects 0.000 claims description 21
- 230000001174 ascending effect Effects 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 10
- 238000004141 dimensional analysis Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 4
- 235000013399 edible fruits Nutrition 0.000 claims description 4
- 239000003973 paint Substances 0.000 claims description 4
- 238000003325 tomography Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 238000010304 firing Methods 0.000 claims description 2
- 108091006146 Channels Proteins 0.000 description 106
- 238000012545 processing Methods 0.000 description 10
- 235000012489 doughnuts Nutrition 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 101150118300 cos gene Proteins 0.000 description 3
- 238000000151 deposition Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000012141 concentrate Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002017 high-resolution X-ray diffraction Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 1
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 1
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000000819 phase cycle Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
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
本发明实施例提供了一种裂缝确定方法和装置，其中，该方法包括：生成研究区内各个CDP点的角道集的方位超道集和角道集的方位超道集的数据库文件；根据上述数据，通过单点分析，确定裂缝预测参数；根据裂缝预测参数、各个CDP点的角道集的方位超道集和角道集的方位超道集的数据库文件，得到研究区的裂缝预测结果数据体；根据实测结果，确定裂缝预测结果数据体是否满足预设要求，如果不满足预设要求，则校正裂缝预测参数，直至裂缝预测结果数据体满足预设要求，将此时的裂缝预测结果数据体确定为该研究区的实际裂缝数据体。由于该方案利用超道集通过方位梯度HOAGVAA法确定裂缝，因而解决了现有技术预测精度差、误差大的问题，提高了裂缝预测的准确性。
Description
技术领域
本发明涉及地质勘探技术领域，特别涉及一种裂缝确定方法和装置。
背景技术
在地质勘探领域中，含油气储层一般大多为具有一定裂缝的地层。其中，储层中往往存在着大量发育的小级别裂缝及微裂缝，这些裂缝可以改造储层的储集空间，使储层孔隙度大大提高，从而使得储层更好地存储石油，这样形成的微裂缝集中发育区的油气藏就是裂缝型油气藏。因此，裂缝预测是准确进行含油储层预测的一项重要内容。
具体实施时，目前，裂缝确定方法主要有以下三种：第一种，利用横波、转换波及多分量检测裂缝；第二种，利用多方位VSP(Vertical Seismic Profiling，垂直地震剖面)检测裂缝；第三种，利用纵波分析检测裂缝。
在实际实施过程中，第一种方法由于使用的横波信噪比低、频率低，得到的裂缝结果往往精度比较低，另外该方法使用的多波采集和处理的成本比较高，转换波的处理过程也很复杂；第二种方法的采集和处理过程相对比较复杂，并且成本也比较高，而且由于方法自身的局限性，其作用范围有限；第三种方法由于一般多是直接利用原始OVT道集和原始CRP(Common Reflection Point，共反射点)道集进行裂缝的确定，而所使用的原始OVT道集和原始CRP道集大都信噪比较低、分辨率较低，且常常含有较多噪音，因此通过第三种方法得到的裂缝预测结果往往存在精度低、误差大的技术问题。
针对现有的裂缝确定方法所存在的处理过程复杂、精度低的问题，目前尚未提出有效的解决方案。
发明内容
针对目前裂缝预测方法存在的技术问题，本发明实施例提供了一种裂缝确定方法和装置，以达到兼顾计算复杂度和裂缝识别精度的目的。
本发明实施例提供了一种裂缝确定方法，包括：
根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件；
根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
根据所述裂缝预测参数、所述各个CDP点的符合规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，求解得到所述研究区的叠前裂缝预测结果数据体；
根据所述研究区的叠后断层预测结果、所述研究区的地震剖面和所述研究区的成像测井的实测结果，确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求；
如果所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，则将所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
在一个实施方式中，在确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求之后，所述方法还包括：
如果所述研究区的叠前裂缝预测结果数据体不满足所述预设精度要求，则对所述裂缝预测参数进行校正，直至根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，并将根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
在一个实施方式中，根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件，包括：
将所述多个CDP点的单方位上的所有主测线上的CRP道集文件合并，得到所述多个CDP点的单方位上的CRP道集数据体；
获取所述多个CDP点的单方位上的CRP道集数据体的低频信息；
根据所述多个CDP点的单方位上的CRP道集数据体和所述多个CDP点的单方位上的CRP道集数据体的低频信息，通过高分辨率层序地层处理，得到所述多个CDP点的单方位上的高分辨率的CRP道集数据体；
利用层速度和所述多个CDP点的单方位上的高分辨率的CRP道集数据体，生成所述多个CDP点的单方位上的角道集；
根据所述研究区的目的层，截短所述多个CDP点的单方位上的角道集，得到所述多个CDP点的单方位上截短后的角道集；
通过格式转换，将所述多个CDP点的单方位上截短后的角道集转换为所述多个CDP点的单方位上预设格式的角道集；
根据所述多个CDP点的单方位上预设格式的角道集，生成所述多个CDP点的单方位上的角道集的数据库文件；
合并所述多个CDP点的各个单方位上预设格式的角道集，得到所述多个CDP点中各个CDP点的角道集的方位超道集；
修改所述多个CDP点中各个CDP点的角道集的方位超道集的SEG-Y道头，得到所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集；
根据所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件。
在一个实施方式中，所述层速度，包括：HANDVEL格式层速度或RMS均方根速度。
在一个实施方式中，获取所述多个CDP点的单方位上的CRP道集数据体的低频信息，包括：
对所述多个CDP点的单方位上的CRP道集数据体的多次波、声波、绕射波和断面波进行低频处理，获取所述多个CDP点的单方位上的CRP道集数据体的低频信息。
在一个实施方式中，根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数，包括：
根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的所述研究区的CDP点分布范围；
从所述用于裂缝预测的所述研究区的CDP点分布范围中选取一点作为进行所述单点分析的单点，并将所述单点的位置作为所述单点的空间参考位置，其中，所述单点的预设区域范围内存在用于获得参考结果的地震资料；
对所述研究区的目的层进行层位处理，得到所述目的层的层位地震解释数据，以所述目的层的层位地震解释数据作为目标约束文件；
根据所述单点的空间参考位置和所述目标约束文件，将所述目的层中用于叠后预测的断层的时间位置作为所述单点的时间参考位置；
根据所述单点的空间参考位置、所述单点的时间参考位置、所述入射角范围和所述方位角范围，对所述单点进行裂缝预测，得到所述单点的裂缝预测结果；
将所述单点的裂缝预测结果与所述单点的参考结果进行比较，得到比较误差，如果所述比较误差小于等于预设阈值，则确定所述入射角范围、所述方位角范围和所述预测时窗为裂缝预测参数；
如果所述比较误差大于所述预设阈值，则调整所述入射角范围、所述方位角范围和所述预测时窗，直到根据调整后的所述入射角范围、所述方位角范围和所述预测时窗得到的所述单点的比较误差小于等于所述预设阈值，则将调整后的所述入射角范围、调整后的所述方位角范围和调整后的所述预测时窗作为所述裂缝预测参数。
在一个实施方式中，根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的所述研究区的CDP点分布范围，包括：
根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取所述各个CDP点的入射角、方位角和覆盖次数；
求取所述入射角的最大值和最小值，求取所述方位角的最大值和最小值；
根据所述入射角的最大值和最小值与所述方位角的最大值和最小值，通过所述入射角和所述方位角的交绘分析，确定所述入射角的范围和所述方位角的范围；
根据所述覆盖次数的平面分析，确定用于裂缝预测的所述研究区的CDP点分布范围；
根据所述多个CDP点的单方位上预设格式的角道集的同相轴拉平程度，确定所述预测时窗。
在一个实施方式中，根据所述单点的空间参考位置、所述单点的时间参考位置、所述入射角范围和所述方位角范围，对所述单点进行裂缝预测，得到所述单点的裂缝预测结果，包括：
根据所述方位角范围，按照所述方位角由小到大的顺序，将所述单点的符合预设规则的角道集的方位超道集划分为多组方位超道集，其中，所述多组方位超道集中的各组方位超道集的覆盖次数是相同的；
根据所述入射角范围，按照所述入射角由小到大的顺序，将所述多组方位超道集中的各组方位超道集分别划分为多个子单元，其中，所述多个子单元中各个子单元的方位超道集的覆盖次数是相同的；
按照所述入射角由小到大的顺序，逐次确定所述多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，并根据所述多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，确定所述多个子单元中各个子单元的裂缝预测结果；
根据所述多个子单元中各个子单元的裂缝预测结果，通过中值滤波，求解得到所述单点的裂缝预测结果。
在一个实施方式中，所述裂缝预测结果包括：裂缝密度和裂缝方向。
在一个实施方式中，所述用于获得参考结果的地震资料包括以下至少之一：测井资料，和/或，用于叠后预测的地震资料。
在一个实施方式中，对所述研究区的目的层进行层位处理，包括：
根据所述研究区的目的层的合成记录标定出所述研究区的目的层的层位位置；
根据所述目的层的层位位置，通过地震资料解释得到所述目的层的层位的地震解释层位数据；
将所述目的层的层位的地震解释层位数据按照预设格式进行转换，将转换得到的结果作为目标约束文件。
在一个实施方式中，对所述裂缝预测参数进行校正，包括：
根据所述研究区的叠后断层预测结果、所述研究区地震剖面和所述研究区的成像测井的实测结果，调整所述入射角范围、所述方位角范围和所述预测时窗，根据校正后的所述入射角范围、校正后的所述方位角范围和校正后的所述预测时窗，通过单点分析，确定校正后的裂缝预测参数。
本发明实施例还提供了一种裂缝确定装置，包括：
预处理单元，用于根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件；
确定裂缝预测参数单元，用于根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
求解裂缝预测结果数据体单元，用于根据所述裂缝预测参数、所述各个CDP点的符合规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，求解得到所述研究区的叠前裂缝预测结果数据体；
校正单元，用于根据所述研究区的叠后断层预测结果、所述研究区的地震剖面和所述研究区的成像测井的实测结果，确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求；如果所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，则将所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体；如果所述研究区的叠前裂缝预测结果数据体不满足所述预设精度要求，则对所述裂缝预测参数进行校正，直至根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，并将根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
本发明实施例利用超道集通过方位梯度HOAGVAA法确定裂缝，有效解决了现有技术中由于使用的原始数据信噪比低、分辨率低且数据含噪音多导致的预测精度低、误差大的技术问题，提高了裂缝预测的准确性。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是根据本发明实施例的裂缝确定方法的处理流程图；
图2是根据本发明实施例的裂缝确定方法中的生成方位超道集的处理流程图；
图3是根据本发明实施例的裂缝确定方法中的对单点应用高精度方位梯度HOAGVAA法确定裂缝的具体处理流程图；
图4；是根据本发明实施例的裂缝确定方法中的结合圆环具体使用高精度方位梯度HOAGVAA法的处理过程示意图；
图5是根据本发明实施例的裂缝确定装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
考虑到目前主要采用的裂缝确定方法，大都直接利用原始OVT道集和原始CRP道集进行裂缝预测。而实际实施时，由于原始OVT道集和原始CRP道集信噪比低、分辨率低，且含有的噪音多，导致直接利用原始OVT道集和原始CRP道集通过常规方法预测裂缝，得到的裂缝结果时往往存在精度低、误差大的技术问题。针对产生上述技术问题的根本原因，考虑可以先对原始数据进行预处理以提高用于确定裂缝的相关数据的信噪比和分辨率，降低所使用数据中含有的噪音，进一步，再通过利用高精度的方位梯度HOAGVAA(High OrderAmplitude Gradient versus AzimuthalAngle，高阶振幅梯度随方位角变化)法确定裂缝，从而有针对性地解决现有方法中存在的裂缝确定精度差、误差大的技术问题，有效地提高了裂缝确定的精度和准确性。
基于上述思路，本发明实施例提供了一种裂缝确定方法，如图1所述，该方法具体可以包括：
步骤101：根据研究区内多个CDP点的单方位CRP道集文件，生成该研究区内多个CDP点的单方位角道集，进而生成多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成各个CDP点的角道集的方位超道集的数据库文件；
上述实施方式中的CDP点(Common Depth Point，共深度点)和CRP点(CommonReflecrtion Point，共反射点)都是地震资料具体的采集过程中所涉及的点。具体地，在地震资料采集过程中，当反射界面水平时，在测线上不同的共炮点道集中，总能找到不同的道，它们都来自地下界面上的某个共同点，该点称为共深度点(CDP)或共反射点(CRP)，这些具有共同深度或反射点的相应各记录道组成共深度点或共反射点道集称为CDP道集或CRP道集。
考虑到目前主要使用的裂缝确定方法，大都由于直接使用未经处理的原始OVT道集和CRT道集数据，这类原始数据由于未经相应处理，往往精度低、分辨率低，且含有较多的噪音，这些特点都会影响最后得到的裂缝预测结果，导致裂缝预测的精度低、误差大。针对上述情况，本发明实施例考虑可以将研究区内多个CDP点单方位的初始CRP道集数据通过一系列的预处理，转换成相对精度更高、分辨率更高的CDP点的角道集的方位超道集，再利用CDP点的角道集的方位超道集对裂缝进行预测和确定。本实施例实际上在确定裂缝之前，对研究区内多个CDP点的相关数据都进行了处理，即，将多个CDP点的单方位的初始CRP道集数据转换生成CDP点的角道集的方位超道集的预处理，具体处理过程可以如图2所示，包括：
步骤201：将研究区多个CDP点的单方位上的所有主测线上的CRP道集文件合并，得到多个CDP点的单方位上的CRP道集数据体；
在实际实施过程中，由于采集得到的数据量比较大，研究区内所有CDP点的每个方位上的CRP道集数据都被分成了若干个文件来记盘。在根据本发明实施例进行裂缝确定工作之前，首先要将各个单方位上的所有CRP道集文件合并起来，形成该方位一个完整的CRP道集数据体。在一个具体的实施例中，可以先将360°全方位CRP道集均匀地等分成6等份，例如，可以等分为以下6个方位：0°～30°-180°～210°、30°～60°-210°～240°、60°～90°-240°～270°、90°～120°-270°～300°、120°～150°-300°～330°、150°～180°-330°～380°，并分别建立相应的6个子目录：0-30_180-210、30-60_210-240、60-90_240-270、90-120_270-300、120-150_300-330、150-180_330-380用来存放对应范围的数据。例如，其中一个方位子目录“0-30_180-210”具体用于保存相应方位“0°～30°-180°～210°”的CRP道集数据。另外，在实际操作时，例如165°角度是150°～180°的中心方位，通常可以直接利用这个中心方位来代表该方位范围。这样记录，可以达到记录简练、表达简单的效果。
步骤202：获取多个CDP点的单方位上的CRP道集数据体的低频信息；
在一个具体实施方式中，为了通过后续步骤203的高分辨层序地层处理，得到进行裂缝预测时所要用的精度更好、效果更佳的CDP点的单方位上的高分辨率的CRP道集数据体之前首先要获取CDP点的单方位上的CRP道集数据体的低频信息。因此，根据步骤202对步骤201得到的多个CDP点的单方位上的CRP道集数据体进行处理。即，可以通过步骤202进行低频处理，得到多个CDP点单方位上的CRP道集数据体的低频信息，其中，步骤202涉及的具体去噪处理具体可以包括：对多个CDP点的单方位上的CRP道集数据体的多次波、声波、绕射波和断面波进行低频处理，获取多个CDP点的单方位上的CRP道集数据体的低频信息。
步骤203：根据多个CDP点的单方位上的CRP道集数据体和多个CDP点的单方位上的CRP道集数据体的低频信息，通过高分辨率层序地层处理，得到多个CDP点的单方位上的高分辨率的CRP道集数据体；
上述实施方式中的多个CDP点的单方位上的高分辨率的CRP道集数据体实际上也是一种单方位上的CRP道集数据体的低频信息。
上述实施方式中的高分辨率层序地层处理一般是指应用高分辨层序地层学理论，对目标地层进行处理，从而预测相关地层分布型式和具体相类型的分布规律。其中，上述高分辨层序地层学是对地层记录中反映的基准面变化旋回的时间地层单元进行二元划分，该理论及技术应用的关键是识别地层中多级次的基准面旋回和等时对比，核心内容为基准面旋回变化过程中，对应A/S比值的变化在相同沉积相体系域或相域中发生沉积物体积分配作用和相分异作用，以及其所属导致的沉积物保存程度、堆积样式、相序、相类型及岩石结构的变化，而上述变化是相关地层在基准面旋回中所处位置的函数，因而地层分布型式和相类型的分布规律是可以预测的。另外，通过高精度分辨率的基准面旋回等时对比，还可以将钻井的一维信息转换为三维信息。
步骤204：利用层速度和多个CDP点的单方位上的高分辨率的CRP道集数据体，生成多个CDP点的单方位上的角道集；
上述实施例中的层速度一般指在层状地层中地震波传播的速度。它直接反映地层的岩性，能用来划分地层。一般是用地震测井或声波测井测得，并且指的是纵波的速度。也可以利用反射记录计算得到。在地震勘探中，一般把层速度低于1400米/秒的地层叫做低速层，把高于3500米/秒的地层视为高速层。但是，习惯上在折射波法中高速层是相对的概念。具体在本发明实施例中，所涉及的层速度一般可以包括：HANDVEL格式层速度或RMS均方根速度。其中，HANDVEL是一种文本格式速度记录的标识符，RMS是具体应用时对相关特定公式进行处理得到的均方根速度。然而，值得注意的是，上述所列举的层速度仅是一种示意性描述，是为了更好地说明本发明，在实际执行的过程中，还可以包括其它可用于辅助生成CDP点各个单方位上的角道集的层速度，具体可以根据实际需要进行选取，本申请对此不作限定。
步骤205：根据研究区的目的层，截短多个CDP点的单方位上的角道集，得到多个CDP点的单方位上截短后的角道集；
实际实施时，由于步骤204获得的多个CDP点的各个单方位上的角道集含有丰富的资料数据，而本发明实施例提出的裂缝确定方法具体实施时是根据具体要求仅对某个目的层进行裂缝的确定，所以，多个CDP点的各个单方位上的角道集所包含的其他层段的信息在本发明实施例中并没有被直接用处。结合实际情况，为了减少存储空间的占用，以免无用信息浪费存储资源，并提高具体实施过程的效率，本发明实施例通过步骤205根据研究区的目的层，对多个CDP点的各个单方位上的角道集进行截短处理，删除其所包含的其他层段的信息，仅保留与所研究的目的层有关的信息，作为后续步骤使用的多个CDP点的各个单方位上截短后的角道集。
步骤206：通过格式转换，将多个CDP点的单方位上截短后的角道集转换为多个CDP点的单方位上预设格式的角道集；
需要说明的是，在具体实施时，由于后续步骤无法直接对多个CDP点的单方位上截短后的初始格式的角道集进行处理，因此可以通过步骤206对步骤205得到的数据进行格式转换，将初始格式转换成后续步骤能处理的格式。具体的格式转换过程可以包括：将多个CDP点的单方位上截短后的角道集以SGY IBM 32float格式输出，得到上述的多个CDP点的单方位上预设格式的角道集。
步骤207：根据多个CDP点的单方位上预设格式的角道集，生成多个CDP点的角道集的数据库文件；
具体实施时，将多个CDP点的各个单方位上预设格式的角道集里具体的各种数据根据规定格式存储在上述多个CDP点的角道集数据库文件中，并由该CDP点的角道集数据库文件进行统一、集中管理。这种处理不仅可以准确地存储根据本发明实施例的裂缝确定方法所要使用的具体的各种数据，即，可以依次存放多个CDP点的各个单方位上预设格式的角道集中的InLine线号、XLine线号、角道集起始道序号、角道集终了道序号和覆盖次数等数据；还能有效地管理、使用上述相应的具体数据。从而，提高应用本发明实施例的来进行裂缝确定的实施效率。
步骤208：合并多个CDP点的各个单方位上预设格式的角道集，得到多个CDP点中各个CDP点的角道集的方位超道集；
在一个具体实施方式，步骤208可以包括：将同InLine线号和同XLine线号的同位置的所有方位角道集合并在一起，从而得到多个CDP点中各个CDP点的角道集的方位超道集。其中，上述实施方式中的超道集一般是指把相邻的好几个CMP/CRP点的地震道全部合在一起，构成1个新的道集，这个新的道集称之为超道集。其中，使用超道集的优点可以包括：1、通过提高覆盖次数提高了信噪比。2、缩小了空间采样间隔，防止空间假频。正因为这两个优点，超道集在实际施工的处理、解释中都有广泛应用：在处理中，可以提高速度分析精度、提高多次波压制效果、提高模型道信噪比；在解释中，因为提高了空间采样率，所以在AVO、AVA等属性分析、角道集抽取等方面都有一定优势。
步骤209：修改多个CDP点中各个CDP点的角道集的方位超道集的SEG-Y道头，得到多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集；
在具体实施过程中，通过步骤208得到的各个CDP点的角道集的方位超道集的格式还不能完全符合后续处理所要求的格式，还需要对各个CDP点的角道集的方位超道集进行进一步格式处理，具体可以根据步骤209修改多个CDP点中各个CDP点的角道集的方位超道集的SEG-Y道头内一条或多条相关信息。在本发明实施例中，该具体修改步骤一般可以是增添方位角道头，使得各个CDP点的角道集的方位超道集的SEG-Y道头符合后续处理要求的预设规则，以便后续进行裂缝确定时可正确使用本数据；再将修改后的各个CDP点的角道集的方位超道集作为多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集。
步骤210：根据多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，生成多个CDP点中各个CDP点的角道集的方位超道集的数据库文件。
类似于前面生成各个CDP点的角道集的数据库文件的过程，具体实施时，将多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集里具体的各类数据根据规定格式存储在各个CDP点的角道集的方位超道集的数据库文件中，并由各个CDP点的角道集的方位超道集的数据库文件统一、集中管理。这种处理不仅可以准确地存储根据本发明实施例的进行裂缝确定时所要使用的具体数据，还能有效地管理、使用相应的数据，从而提高了根据本发明实施例进行裂缝确定的实施效率。
步骤102：根据各个CDP点的符合预设规则的角道集的方位超道集和各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
为了提高对研究区内所有CPD点进行裂缝预测的效率和精度，可以考虑在研究区中先选取某个具有测井或具有可用于进行叠后预测的地震资料的CDP点作为实验单点进行单点分析。具体可以先对该点应用高精度方位梯度HOAGVAA法进行裂缝确定，得到裂缝预测结果，将得到的裂缝预测结果与根据该CDP点附近的测井的测井资料得到的参考结果或者根据该CDP点附近的地震资料进行叠后预测得到的参考结果进行比较，根据比较情况有针对性地调整裂缝预测时具体涉及的裂缝预测参数，直到基于调整后的裂缝预测参数得到的裂缝预测结果与根据测井资料得到的参考结果或者根据该CDP点附近的地震资料通过叠后预测得到的参考结果的误差在可接受的范围内，将这时的裂缝预测参数确定为后续步骤使用的裂缝预测参数，并将该裂缝预测参数推广应用到该研究区内所有CDP点的裂缝确定的过程中，从而避免了对研究区内所有CDP点每个CDP点分别选择、调试裂缝预测参数，以确定最优的裂缝预测参数的繁琐过程，提高了裂缝确定过程的效率和准确度。基于上述思路，步骤102具体可以包括：
S1：根据研究区内各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的该研究区的CDP点分布范围；
具体实施时，在后续步骤通过单点分析确定裂缝预测参数之前，本发明实施例要先根据研究区内CDP点的符合预设规则的角道集的方位超道集和CDP点的角道集的方位超道集的数据库文件，通过步骤S1获取用于单点分析的具体参数，这些具体参数一般可以包括：覆盖次数、预测时窗、入射角范围和方位角范围，具体实施时，S1可以包括：
根据各个CDP点的符合预设规则的角道集的方位超道集，获取各个CDP点的入射角、方位角和覆盖次数；
求取入射角的最大值和最小值，求取方位角的最大值和最小值；
根据入射角的最大值和最小值与方位角的最大值和最小值，通过入射角和方位角的交绘分析，确定入射角的范围和方位角的范围；
根据覆盖次数的平面分析，确定用于裂缝预测的研究区的CDP点分布范围；
根据多个CDP点的单方位上预设格式的角道集的同相轴拉平程度，确定所述预测时窗。
需要说明的是，在上述操作中，通过对入射角与方位角的交绘分析，可以分析得到数据的入射角方位展布特征，根据得到具体的展布特征，可以为后续的裂缝确定选择合适、恰当的且具有较为宽广方位分布的入射角范围，从而提高根据本发明实施例进行裂缝确定的精度和准确性。
需要说明的是上述实施方式中所获取的覆盖次数一般可以反映CDP点的采集数据量，一般情况下，一个CDP点的覆盖次数越多，意味着该CDP点的采集数据量就越多，相应地，利用采集数据对该CDP点进行裂缝预测得到的裂缝预测结果也就越准确。在具体实施时，一般由于获取采集数据的实测方法自身的原因，研究区边界附近的CDP点的覆盖次数相对会比较少，对应的，这些位置上的点的采集数据量也相对比较少，在正式对该研究区的进行裂缝确定时，本发明实施例为了提高该研究区的裂缝预测结果数据整体的准确度，一般会根据CDP点的覆盖次数，删除例如研究区边界附近的覆盖次数比较少的CDP点，将覆盖次数比较多的CDP点作为用于实际裂缝预测的该研究区的CDP点分布范围。在正式对该研究区进行裂缝确定时，一般只根据上述用于裂缝预测的研究区的CDP点分布范围中各个CDP点的裂缝预测结果，整理得到该研究区的裂缝预测结果。
S2：从用于裂缝预测的所述研究区的CDP点的分布范围中选取一点作为进行单点分析的单点，并将该单点的位置作为单点的空间参考位置，其中，该单点的预设区域范围内存在用于获得参考结果的地震资料；
上述实施方式中的用于获得参考结果的地震资料一般可以包括以下至少之一：单点的预设区域范围内的测井资料，和/或，单点的预设区域范围内的用于叠后预测的地震资料。然而，值得注意的是，上述所列举的地震资料仅是一种示意性描述，是为了更好地说明本发明，在实际执行的过程中，还可以包括其它相关的可以用于得到后续步骤与预测结果进行比较的参考结果的地震资料，可以根据实际需要进行选取，本申请对此不作限定。根据具体已有的用于获得参考结果的地震资料，具体实施时，在研究区的CDP点分布范围中一般选择CDP点附近有测井资料或者有可用于进行叠后预测的地震资料的CDP点作为单点分析的单点。相应的，如果选择的单点附近有测井资料，则步骤S6中用于与步骤S5获得的单点裂缝预测结果进行的比较的参考结果就是该CDP点附近的测井探测结果；如果选择的单点附近有用于叠后预测的地震资料，则步骤S6中用于与步骤S5获得的单点裂缝预测结果进行的比较的参考结果就是根据该CDP点附近的地震资料通过叠后预测得到的叠后断层预测结果。
S3：对研究区的目的层进行层位处理，得到目的层的层位地震解释数据，以目的层的层位地震解释数据作为目标约束文件；
考虑到本发明实施例只针对某个待研究的目的层进行裂缝确定，在实际实施时，可以根据具体要研究的目的层，即对要研究的目的层进行层位处理，得到用于后续裂缝确定过程的目标约束文件。因此，步骤S3对研究区的目的层进行层位处理具体可以包括：根据研究区的目的层的合成记录标定出该研究区的目的层的层位位置；根据该目的层的层位位置，通过地震资料解释得到该目的层的层位的地震解释层位数据；再将目的层的层位的地震解释层位数据按照预设格式进行转换，将转换得到的结果作为目标约束文件。其中，具体实施时，研究区目的层地震合成记录的标定可以用于识别和定位目标层段顶界和底界面在地震剖面上的纵向时间位置；根据合成记录标定出的层位位置，在所有主测线和联络测线地震剖面上解释出该层位；实际执行时，该层位从解释软件中输出时都是按“主测线—联络测线—顶(底)解释时间T0”格式来存放；因此要将这种“主测线—联络测线—顶(底)界面时间T0”格式转换成本案层位格式即“主测线—联络测线—顶界面时间T0—底界面间T0”，这就是上述实施方式中的目标约束文件。
S4：根据单点的空间参考位置和目标约束文件，将目的层中用于叠后预测的断层的时间位置作为该单点的时间参考位置；
S5：根据单点的空间参考位置、单点的时间参考位置、入射角范围和方位角范围，应用高精度方位梯度HOAGVAA法对单点进行裂缝预测，得到单点的裂缝预测结果；
具体实施时过程中，为了根据该研究区的各个CDP点的符合预设规则的角道集的方位超道集得到具体应用高精度方位梯度HOAGVAA法进行裂缝确定时所需要的裂缝预测参数，在一个实施例中，具体可以按照下面的步骤执行：根据研究区各个CDP点的符合预设规则的角道集的方位超道集，获取各个CDP点的入射角、方位角和覆盖次数；求取入射角的最大值和最小值，求取方位角的最大值和最小值；根据入射角的最大值和最小值与方位角的最大值和最小值，通过入射角和方位角的交绘分析，确定入射角的范围和方位角的范围；根据覆盖次数的平面分析，确定用于裂缝预测的研究区的CDP点分布范围；根据多个CDP点的单方位上预设格式的角道集的同相轴拉平程度，确定预测时窗。其中，通过上述步骤得到的入射角范围、方位角范围和预测视窗等参数都是对该研究区的所有CDP点应用高精度方位梯度HOAGVAA法进行裂缝确定所要使用的裂缝预测参数。然而，值得注意的是，上述所列举的裂缝预测参数仅是一种示意性描述，是为了更好地说明本发明，在实际执行的过程中，还可以包括其它相关的应用高精度方位梯度HOAGVAA法预测裂缝所涉及的裂缝预测参数，具体可以根据实际需要进行选取，本申请对此不作限定。
考虑到现有裂缝确定方法本身的局限性，常规方法本身存在裂缝预测过程繁琐，且裂缝预测结果精度不高的技术问题，本发明实施例对裂缝确定的方法进行改进，具体通过采用高精度的数据，并应用方位梯度HOAGVAA法对裂缝进行预测，从而提高了裂缝预测的效率和精度。上述实施方式中，对单点应用的高精度方位梯度HOAGVAA法确定裂缝预测结果，具体可以如图3所示，包括：
步骤301：根据方位角范围，按照方位角由小到大的顺序，将该单点的符合预设规则的角道集的方位超道集划分为多组方位超道集，其中，多组方位超道集中的各组方位超道集的覆盖次数是相同的；
步骤302：根据入射角范围，按照入射角由小到大的顺序，将多组方位超道集中的各组方位超道集分别划分为多个子单元，其中，多个子单元中各个子单元的方位超道集的覆盖次数是相同的；
步骤303：按照入射角由小到大的顺序，逐次确定多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，并根据多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，确定多个子单元中各个子单元的裂缝预测结果；
步骤304：根据多个子单元中各个子单元的裂缝预测结果，通过中值滤波，求解得到该单点的裂缝预测结果。
在一个具体实施例中，结合图4的处理过程示意图，利用圆环作为辅助，具体说明一下应用上述高精度方位梯度HOAGVAA法具体确定裂缝的具体流程：
1)某个CDP点的所有角道集就是图4的整个圆，将该点所有方位超道集按照方位角由小到达进行等分，等分成多组方位超道集，每一组方位超道集在图4对应一个扇形区，其中各组方位超道集的覆盖次数是相同的；
2)将每组方位超道集即图4中每个窄扇形区内的数据沿入射角方向等分成若干个子单元即分成若干个同心圆环，每一个同心圆环对应一个子单元，其中，每个同心圆环的覆盖次数一样即地震道数是相同的；
3)按入射角大小由小到大顺序依此求解每个同心圆环的不同方位的振幅随角度变化而变化的地震响应强度变化曲线的曲率或梯度(方位梯度)，即每个子单元的方位超道集的地震响应强度变化曲线的曲率或梯度；
4)按圆环逐环做叠前裂缝预测，即用第1个同心圆环上方位A梯度来进行叠前裂缝预测，求得第1个同心圆环地震数据所预测得到的裂缝参数(裂缝密度、裂缝方向)；用第2个同心圆环方位梯度来进行叠前裂缝预测，求得第2个同心圆环地震数据所预测得到的裂缝参数(裂缝密度、裂缝方向)……如此可以求得所有同心圆环地震数据所预测得到的裂缝参数(裂缝密度、裂缝方向)，即求解得到了每个子单元的裂缝预测结果；
5)最终整理所有叠前裂缝预测结果，对各分环叠前裂缝预测的裂缝参数进行中值滤波等处理，即通过整理各个子单元的裂缝预测结果得到最终叠前裂缝预测结果即最终裂缝密度和裂缝方向。
需要说明的是，之所以能够应用上述的高精度方位梯度HOAGVAA法进行较为准确的裂缝确定，其具体原理如下所述：
具体实施时，纵波方位属性如振幅等不仅与入射角有关，而且还与方位角有关，一般可表示为：
F(θ,φ)＝A(θ)+B(θ)cos2φ+C(θ)cos22φ (1)
其中，F表示访问属性如P波的反射时间或反射振幅或速度等；θ是入射角；Φ是相对裂缝方向的方位角；A(θ)、B(θ)、C(θ)是与方位角无关的系数。
在某一固定方位，则有：
R(θ)＝A+Bsin2θ+Ctan2θsin2θ (2)
其中
对(2)变形有
R(θ)＝A+Bsin2θ+Csin4θ/cos2θ
令a＝A,b＝B,c＝C/cos2θ,x＝sin2θ,y＝R(θ)则近似的有
y＝a+bx+cx2 (3)
对式(3)求导有
y′＝b+cx (4)
令g＝y′,其中，g为梯度函数。
(4)表示振幅变化是线性分段变化的，即振幅在纵向上按偏移距或入射角可划分为若干段，每段的振幅变化是线性的，即：
这样，就可以对某个方位上的随入射角或偏移距的变化的地震响应强度变化曲线(AVA)进行剖分，即：对整个AVA变化曲线进行剖分，使每一剖分段都趋线性化变化，进而求解每一剖分段的AVA曲线梯度，使得求解出的梯度更精确，更能准确地反映振幅的变化规律，从而增加裂缝预测的精度和可靠性。
S6：将单点的裂缝预测结果与单点的参考结果进行比较，得到比较误差，如果比较误差小于等于预设阈值，则确定此时的入射角范围、方位角范围和预测时窗为裂缝预测参数；
S7：如果比较误差大于预设阈值，则调整入射角范围、方位角范围和预测时窗，直到根据调整后的入射角范围、方位角范围和预测时窗得到的单点的比较误差小于等于预设阈值，则将调整后的入射角范围、调整后的方位角范围和调整后的预测时窗作为裂缝预测参数。
上述实施方式中，需要进一步说明的是得到的到的裂缝预测结果一般可以包括裂缝密度数据和裂缝方向数据。其中，通过裂缝密度数据可以确定该CDP点的裂缝密度，通过裂缝方向数据可以确定该CDP点的裂缝位置。因此，通过这上述个裂缝预测结果数据可以基本确定该CDP点的裂缝。
步骤103：根据裂缝预测参数、各个CDP点的符合预设规则的角道集的方位超道集和各个CDP点的角道集的方位超道集的数据库文件，求解得到该研究区的叠前裂缝预测结果数据体；
类似于单点分析中对所单点的裂缝确定，同样，对于用于裂缝预测的该研究区CPD点分布范围内的各个点，根据步骤102单点分析确定的裂缝预测参数的处理过程，通过步骤103具体应用高精度方位梯度HOAGVVA法对用于裂缝预测的该研究区CDP点范围内各个CDP点进行裂缝确定，再根据该研究区CDP点范围内各个CDP点的裂缝预测结果数据，整理得到该研究区的叠前裂缝预测结果数据体。
步骤104：根据该研究区的叠后断层预测结果、该研究区的地震剖面和该研究区的成像测井的实测结果，确定该研究区的叠前裂缝预测结果数据体是否满足预设精度要求；
步骤105：如果该研究区的叠前裂缝预测结果数据体满足预设精度要求，则将该研究区的叠前裂缝预测结果数据体确定为该研究区的实际裂缝数据体。
需要说明的是，利用上述步骤103求解得到的研究区的叠前裂缝预测结果数据体进行步骤104与该研究区的叠后断层预测结果、该研究区的地震剖面和该研究区的成像测井的实测结果进行比较时，可能出现不满足预设精度要求的情况，这时得到的该地区的叠前裂缝预测结果与实际的裂缝预测结果还存在较大的误差，不能直接将其通过步骤105确定为该研究区的实际裂缝数据体，因此还要根据比较结果，对裂缝预测参数进行一次或多次的校正，直到基于校正后的裂缝预测参数得到裂缝预测结果数据体接近或等于实际裂缝数据体，再将这时的裂缝预测结果数据体确定为实际的裂缝数据体。具体处理，可以包括：在确定该研究区的叠前裂缝预测结果数据体是否满足预设精度要求之后，如果该研究区的叠前裂缝预测结果数据体不满足预设精度要求，则对裂缝预测参数进行校正，即，根据与实测结果的比较的情况针对性地调整入射角范围、方位角范围和预测时窗，直至根据校正后的裂缝预测参数得到的研究区的叠前裂缝预测结果数据体满足预设精度要求，并将根据校正后的裂缝预测参数得到的研究区的叠前裂缝预测结果数据体确定为该研究区的实际裂缝数据体。
上述实施方式中的校正处理，一般可以包括：根据该研究区的叠后断层预测结果、该研究区地震剖面和该研究区的成像测井的实测结果，调整入射角范围、方位角范围和预测时窗，根据校正后的入射角范围、校正后的方位角范围和校正后的预测时窗，通过单点分析，确定校正后的裂缝预测参数。
本发明实施例利用超道集通过方位梯度HOAGVAA法确定裂缝，有效解决了现有技术中由于直接使用的原始数据的信噪比低、分辨率低且数据含噪音多导致的预测精度低、误差大的技术问题，提高了裂缝预测的准确性。
基于同一发明构思，本发明实施例中还提供了一种裂缝确定装置，如下面的实施例所述。由于裂缝确定装置解决问题的原理与裂缝确定方法相似，因此裂缝确定装置的实施可以参见裂缝确定方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图5是本发明实施例的裂缝确定装置的一种结构框图，如图5所示，可以包括：预处理模块501、确定裂缝预测参数模块502、求解裂缝预测结果数据体模块503和校正模块504，下面对该结构进行说明。
预处理模块501，用于根据研究区内多个CDP点的单方位CRP道集文件，生成多个CDP点的单方位角道集，进而生成多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成多个CDP点中各个CDP点的角道集的方位超道集的数据库文件；
确定裂缝预测参数模块502，用于根据各个CDP点的符合预设规则的角道集的方位超道集和各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
求解裂缝预测结果数据体模块503，用于根据裂缝预测参数、各个CDP点的符合预设规则的角道集的方位超道集和各个CDP点的角道集的方位超道集的数据库文件，求解得到该研究区的叠前裂缝预测结果数据体；
校正模块504，用于根据研究区的叠后断层预测结果、研究区的地震剖面和研究区的成像测井的实测结果，确定研究区的叠前裂缝预测结果数据体是否满足预设精度要求；如果研究区的叠前裂缝预测结果数据体满足预设精度要求，则将研究区的叠前裂缝预测结果数据体确定为研究区的实际裂缝数据体；如果研究区的叠前裂缝预测结果数据体不满足预设精度要求，则对裂缝预测参数进行校正，直至根据校正后的裂缝预测参数得到的研究区的叠前裂缝预测结果数据体满足预设精度要求，并将根据校正后的裂缝预测参数得到的该研究区的叠前裂缝预测结果数据体确定为该研究区的实际裂缝数据体。
在一个具体实施例中，预处理模块501对该研究区内多个CDP点所执行的操作具体可以包括：
将多个CDP点的单方位上的所有主测线上的CRP道集文件合并，得到多个CDP点的单方位上的CRP道集数据体；
获取多个CDP点的单方位上的CRP道集数据体的低频信息；
根据多个CDP点的单方位上的CRP道集数据体和多个CDP点的单方位上的CRP道集数据体的低频信息，通过高分辨率层序地层处理，得到多个CDP点的单方位上的高分辨率的CRP道集数据体；
利用层速度和多个CDP点的单方位上的高分辨率的CRP道集数据体，生成多个CDP点的单方位上的角道集；
根据研究区的目的层，截短多个CDP点的单方位上的角道集，得到多个CDP点的单方位上截短后的角道集；
通过格式转换，将多个CDP点的单方位上截短后的角道集转换为多个CDP点的单方位上预设格式的角道集；
根据多个CDP点的单方位上预设格式的角道集，生成多个CDP点的角道集的数据库文件；
合并多个CDP点的各个单方位上预设格式的角道集，得到多个CDP点中各个CDP点的角道集的方位超道集；
修改多个CDP点中各个CDP点的角道集的方位超道集的SEG-Y道头，得到各个CDP点的符合预设规则的角道集的方位超道集；
根据各个CDP点的符合预设规则的角道集的方位超道集，生成各个CDP点的角道集的方位超道集的数据库文件。
在一个具体实施例中，确定裂缝预测参数模块502一般通过单点分析确定裂缝预测参数，具体处理可以包括：
根据各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的该研究区的CDP点分布范围；
从用于裂缝预测的研究区的CDP点分布范围中选取一点作为进行单点分析的单点，并将单点的位置作为单点的空间参考位置，其中，单点的预设区域范围内存在用于获得参考结果的地震资料；
对研究区的目的层进行层位处理，得到目的层的层位地震解释数据，并以目的层的层位地震解释数据作为目标约束文件；
根据单点的空间参考位置和目标约束文件，将目的层中用于叠后预测的断层的时间位置作为单点的时间参考位置；
根据单点的空间参考位置、单点的时间参考位置、入射角范围和方位角范围，对单点进行裂缝预测，得到单点的裂缝预测结果；
将单点的裂缝预测结果与单点的参考结果进行比较，得到比较误差，如果比较误差小于等于预设阈值，则确定入射角范围、方位角范围和预测时窗为裂缝预测参数；
如果比较误差大于预设阈值，则调整入射角范围、方位角范围和预测时窗，直到根据调整后的入射角范围、方位角范围和预测时窗得到的单点的比较误差小于等于预设阈值，则将调整后的入射角范围、调整后的方位角范围和调整后的预测时窗作为裂缝预测参数。
需要说明的是，上述用于获得参考结果的地震资料包括以下至少之一：测井资料，和/或，用于叠后预测的地震资料。
其中，上述实施例中利用确定裂缝预测参数模块502获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的研究区的CDP点分布范围的过程，可以包括：
根据各个CDP点的符合预设规则的角道集的方位超道集，获取各个CDP点的入射角、方位角和覆盖次数；
求取入射角的最大值和最小值，求取方位角的最大值和最小值；
根据入射角的最大值和最小值与方位角的最大值和最小值，通过入射角和方位角的交绘分析，确定入射角的范围和方位角的范围；
根据覆盖次数的平面分析，确定用于裂缝预测的研究区的CDP点分布范围；
根据多个CDP点的单方位上预设格式的角道集的同相轴拉平程度，确定预测时窗。
在一个实施例中，确定裂缝预测参数模块502进行单点分析对单点进行裂缝确定的方法和求解裂缝预测结果数据体模块503对用于裂缝预测的该研究的CDP点范围内的各个CDP点进行裂缝预测得到各个CDP的裂缝预测结果的方法是相同的，一般都是对一个CDP点应用高精度方位梯度HOAGVAA法确定CDP的裂缝。其中，确定裂缝预测参数模块502和求解裂缝预测结果数据体模块503利用该方法确定裂缝的具体过程可以包括：
根据方位角范围，按照方位角由小到大的顺序，将该单点的符合预设规则的角道集的方位超道集划分为多组方位超道集，其中，多组方位超道集中的各组方位超道集的覆盖次数是相同的；
根据入射角范围，按照入射角由小到大的顺序，将多组方位超道集中的各组方位超道集分别划分为多个子单元，其中，多个子单元中各个子单元的方位超道集的覆盖次数是相同的；
按照入射角由小到大的顺序，逐次确定多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，并根据多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，确定多个子单元中各个子单元的裂缝预测结果；
根据多个子单元中各个子单元的裂缝预测结果，通过中值滤波，求解得到该单点的裂缝预测结果。
在一个具体实施例中，校正模块504一般具体利用下面的方法进行校正处理：根据研究区的叠后断层预测结果、研究区地震剖面和研究区的成像测井的实测结果，调整入射角范围、方位角范围和预测时窗，根据校正后的入射角范围、校正后的方位角范围和校正后的预测时窗，通过单点分析，确定校正后的裂缝预测参数。
从以上的描述中，可以看出，本发明实施例实现了如下技术效果：通过对原始的CRP道集进行预处理，得到并使用分辨率更高、精度更高的角道集的超道集进行裂缝确定，有效地解决了现有方法中由于直接使用分辨率低、精度低且含噪音多的原始OVT道集和原始CRP道集导致得到的裂缝预测结果精度低、误差大的技术问题；通过先进行单点分析确定裂缝预测参数，再将上述得到的裂缝预测参数推广应用到研究区内所有CDP点的裂缝确定中，提高了对该研究区所有CDP点进行裂缝确定的效率；通过采用高精度方位梯度HOAGVAA法确定裂缝，充分并有效地应用了高精度数据，并且避免了常规方法存在的局限，减小了裂缝预测的误差，有效提高了裂缝预测的精度和准确性。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (13)
1.一种裂缝确定方法，其特征在于，包括：
根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件；
根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
根据所述裂缝预测参数、所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，求解得到所述研究区的叠前裂缝预测结果数据体；
根据所述研究区的叠后断层预测结果、所述研究区的地震剖面和所述研究区的成像测井的实测结果，确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求；
如果所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，则将所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
2.根据权利要求1所述的方法，其特征在于，在确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求之后，所述方法还包括：
如果所述研究区的叠前裂缝预测结果数据体不满足所述预设精度要求，则对所述裂缝预测参数进行校正，直至根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，并将根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
3.根据权利要求1所述的方法，其特征在于，根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件，包括：
将所述多个CDP点的单方位上的所有主测线上的CRP道集文件合并，得到所述多个CDP点的单方位上的CRP道集数据体；
获取所述多个CDP点的单方位上的CRP道集数据体的低频信息；
根据所述多个CDP点的单方位上的CRP道集数据体和所述多个CDP点的单方位上的CRP道集数据体的低频信息，通过高分辨率层序地层处理，得到所述多个CDP点的单方位上的高分辨率的CRP道集数据体；
利用层速度和所述多个CDP点的单方位上的高分辨率的CRP道集数据体，生成所述多个CDP点的单方位上的角道集；
根据所述研究区的目的层，截短所述多个CDP点的单方位上的角道集，得到所述多个CDP点的单方位上截短后的角道集；
通过格式转换，将所述多个CDP点的单方位上截短后的角道集转换为所述多个CDP点的单方位上预设格式的角道集；
根据所述多个CDP点的单方位上预设格式的角道集，生成所述多个CDP点的单方位上的角道集的数据库文件；
合并所述多个CDP点的各个单方位上预设格式的角道集，得到所述多个CDP点中各个CDP点的角道集的方位超道集；
修改所述多个CDP点中各个CDP点的角道集的方位超道集的SEG-Y道头，得到所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集；
根据所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件。
4.根据权利要求3所述的方法，其特征在于，所述层速度，包括：HANDVEL格式层速度或RMS均方根速度。
5.根据权利要求3所述的方法，其特征在于，获取所述多个CDP点的单方位上的CRP道集数据体的低频信息，包括：
对所述多个CDP点的单方位上的CRP道集数据体的多次波、声波、绕射波和断面波进行低频处理，获取所述多个CDP点的单方位上的CRP道集数据体的低频信息。
6.根据权利要求1所述的方法，其特征在于，根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数，包括：
根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的所述研究区的CDP点分布范围；
从所述用于裂缝预测的所述研究区的CDP点分布范围中选取一点作为进行所述单点分析的单点，并将所述单点的位置作为所述单点的空间参考位置，其中，所述单点的预设区域范围内存在用于获得参考结果的地震资料；
对所述研究区的目的层进行层位处理，得到所述目的层的层位地震解释数据，根据所述目的层的层位地震解释数据作为目标约束文件；
根据所述单点的空间参考位置和所述目标约束文件，将所述目的层中用于叠后预测的断层的时间位置作为所述单点的时间参考位置；
根据所述单点的空间参考位置、所述单点的时间参考位置、所述入射角范围和所述方位角范围，对所述单点进行裂缝预测，得到所述单点的裂缝预测结果；
将所述单点的裂缝预测结果与所述单点的参考结果进行比较，得到比较误差，如果所述比较误差小于等于预设阈值，则确定所述入射角范围、所述方位角范围和所述预测时窗为裂缝预测参数；
如果所述比较误差大于所述预设阈值，则调整所述入射角范围、所述方位角范围和所述预测时窗，直到根据调整后的所述入射角范围、所述方位角范围和所述预测时窗得到的所述单点的比较误差小于等于所述预设阈值，则将调整后的所述入射角范围、调整后的所述方位角范围和调整后的所述预测时窗作为所述裂缝预测参数。
7.根据权利要求6所述的方法，其特征在于，根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取覆盖次数、预测时窗、入射角范围和方位角范围，并确定用于裂缝预测的所述研究区的CDP点分布范围，包括：
根据所述各个CDP点的符合预设规则的角道集的方位超道集，获取所述各个CDP点的入射角、方位角和覆盖次数；
求取所述入射角的最大值和最小值，求取所述方位角的最大值和最小值；
根据所述入射角的最大值和最小值、所述方位角的最大值和最小值，通过所述入射角和所述方位角的交绘分析，确定所述入射角的范围和所述方位角的范围；
根据所述覆盖次数的平面分析，确定用于裂缝预测的所述研究区的CDP点分布范围；
根据所述多个CDP点的单方位上预设格式的角道集的同相轴拉平程度，确定所述预测时窗。
8.根据权利要求6所述的方法，其特征在于，根据所述单点的空间参考位置、所述单点的时间参考位置、所述入射角范围和所述方位角范围，对所述单点进行裂缝预测，得到所述单点的裂缝预测结果，包括：
根据所述方位角范围，按照所述方位角由小到大的顺序，将所述单点的符合预设规则的角道集的方位超道集划分为多组方位超道集，其中，所述多组方位超道集中的各组方位超道集的覆盖次数是相同的；
根据所述入射角范围，按照所述入射角由小到大的顺序，将所述多组方位超道集中的各组方位超道集分别划分为多个子单元，其中，所述多个子单元中各个子单元的方位超道集的覆盖次数是相同的；
按照所述入射角由小到大的顺序，逐次确定所述多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，并根据所述多个子单元中各个子单元的方位超道集的地震响应强度随入射角度变化曲线的梯度，确定所述多个子单元中各个子单元的裂缝预测结果；
根据所述多个子单元中各个子单元的裂缝预测结果，通过中值滤波，求解得到所述单点的裂缝预测结果。
9.根据权利要求8所述的方法，其特征在于，所述裂缝预测结果包括：裂缝密度和裂缝方向。
10.根据权利要求6所述的方法，其特征在于，所述用于获得参考结果的地震资料包括以下至少之一：测井资料，和/或，用于叠后预测的地震资料。
11.根据权利要求6所述的方法，其特征在于，对所述研究区的目的层进行层位处理，得到所述目的层的层位地震解释数据，根据所述目的层的层位地震解释数据作为目标约束文件，包括：
根据所述研究区的目的层的合成记录标定出所述研究区的目的层的层位位置；
根据所述目的层的层位位置，通过地震资料解释得到所述目的层的层位的地震解释层位数据；
将所述目的层的层位的地震解释层位数据按照预设格式进行转换，将转换得到的结果作为目标约束文件。
12.根据权利要求1所述的方法，其特征在于，对所述裂缝预测参数进行校正，包括：
根据所述研究区的叠后断层预测结果、所述研究区地震剖面和所述研究区的成像测井的实测结果，调整所述入射角范围、所述方位角范围和所述预测时窗，根据校正后的所述入射角范围、校正后的所述方位角范围和校正后的所述预测时窗，通过单点分析，确定校正后的裂缝预测参数。
13.一种裂缝确定装置，其特征在于，包括：
预处理单元，用于根据研究区内多个CDP点的单方位CRP道集文件，生成所述多个CDP点的单方位角道集，进而生成所述多个CDP点中各个CDP点的符合预设规则的角道集的方位超道集，并生成所述多个CDP点中各个CDP点的角道集的方位超道集的数据库文件；
确定裂缝预测参数单元，用于根据所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，通过单点分析，确定裂缝预测参数；
求解裂缝预测结果数据体单元，用于根据所述裂缝预测参数、所述各个CDP点的符合预设规则的角道集的方位超道集和所述各个CDP点的角道集的方位超道集的数据库文件，求解得到所述研究区的叠前裂缝预测结果数据体；
校正单元，用于根据所述研究区的叠后断层预测结果、所述研究区的地震剖面和所述研究区的成像测井的实测结果，确定所述研究区的叠前裂缝预测结果数据体是否满足预设精度要求；如果所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，则将所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体；如果所述研究区的叠前裂缝预测结果数据体不满足所述预设精度要求，则对所述裂缝预测参数进行校正，直至根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体满足所述预设精度要求，并将根据校正后的裂缝预测参数得到的所述研究区的叠前裂缝预测结果数据体确定为所述研究区的实际裂缝数据体。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610772174.9A CN106405644B (zh) | 2016-08-30 | 2016-08-30 | 裂缝确定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610772174.9A CN106405644B (zh) | 2016-08-30 | 2016-08-30 | 裂缝确定方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106405644A true CN106405644A (zh) | 2017-02-15 |
| CN106405644B CN106405644B (zh) | 2018-10-16 |
Family
ID=58004018
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610772174.9A Active CN106405644B (zh) | 2016-08-30 | 2016-08-30 | 裂缝确定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106405644B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107153725A (zh) * | 2017-04-20 | 2017-09-12 | 中国地质大学（武汉） | 一种关于陆相层序地层格架构建的成图系列及技术方法 |
| CN109469477A (zh) * | 2018-10-18 | 2019-03-15 | 中国海洋石油集团有限公司 | 一种人工裂缝延伸方向的预测方法和装置 |
| CN109655902A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 用于速度分析的地震数据超道集生成方法及系统 |
| CN110286410A (zh) * | 2019-08-06 | 2019-09-27 | 中国矿业大学(北京) | 基于绕射波能量的裂缝反演方法和装置 |
| CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
| US20020053430A1 (en) * | 2000-06-19 | 2002-05-09 | Curtis Michael P. | Apparatus and methods for applying time lapse VSP to monitor a reservoir |
| CN101907725A (zh) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | 裂缝预测方法和装置 |
| CN102879800A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种横波分裂裂缝检测的方法 |
| CN104316965A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种裂缝方位和强度的预测方法及系统 |
| CN104375182A (zh) * | 2014-11-10 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种裂缝型储层流体的识别方法及装置 |
-
2016
- 2016-08-30 CN CN201610772174.9A patent/CN106405644B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
| US20020053430A1 (en) * | 2000-06-19 | 2002-05-09 | Curtis Michael P. | Apparatus and methods for applying time lapse VSP to monitor a reservoir |
| CN101907725A (zh) * | 2010-06-12 | 2010-12-08 | 中国石油天然气股份有限公司 | 裂缝预测方法和装置 |
| CN102879800A (zh) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | 一种横波分裂裂缝检测的方法 |
| CN104316965A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种裂缝方位和强度的预测方法及系统 |
| CN104375182A (zh) * | 2014-11-10 | 2015-02-25 | 中国石油天然气股份有限公司 | 一种裂缝型储层流体的识别方法及装置 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107153725A (zh) * | 2017-04-20 | 2017-09-12 | 中国地质大学（武汉） | 一种关于陆相层序地层格架构建的成图系列及技术方法 |
| CN107153725B (zh) * | 2017-04-20 | 2020-07-03 | 中国地质大学（武汉） | 一种关于陆相层序地层格架构建的成图系列及技术方法 |
| CN109655902A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | 用于速度分析的地震数据超道集生成方法及系统 |
| CN109469477A (zh) * | 2018-10-18 | 2019-03-15 | 中国海洋石油集团有限公司 | 一种人工裂缝延伸方向的预测方法和装置 |
| CN109469477B (zh) * | 2018-10-18 | 2022-08-02 | 中国海洋石油集团有限公司 | 一种人工裂缝延伸方向的预测方法和装置 |
| CN112305594A (zh) * | 2019-07-23 | 2021-02-02 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
| CN112305594B (zh) * | 2019-07-23 | 2023-10-27 | 中国石油天然气股份有限公司 | 非均质储层的油气分布确定方法及系统 |
| CN110286410A (zh) * | 2019-08-06 | 2019-09-27 | 中国矿业大学(北京) | 基于绕射波能量的裂缝反演方法和装置 |
| CN110286410B (zh) * | 2019-08-06 | 2020-06-09 | 中国矿业大学(北京) | 基于绕射波能量的裂缝反演方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106405644B (zh) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106405644B (zh) | 裂缝确定方法和装置 | |
| US20190302298A1 (en) | Selection of pick values for automated near surface analysis | |
| EP3253946B1 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| US10371842B2 (en) | Anisotropy analysis using direct and reflected arrivals in seismic survey data | |
| CN102176052B (zh) | 一种面向三维层面网格生成的层序分析方法 | |
| CN105259581B (zh) | 一种地震资料时深转换方法 | |
| US11221425B1 (en) | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations | |
| CN107956465A (zh) | 基于关联井的全区多井测井曲线标准化方法及装置 | |
| CN105549084A (zh) | 一种三维高精度速度建模方法和系统 | |
| US20220221604A1 (en) | Method of analysing seismic data | |
| CN104375178A (zh) | 一种碳酸盐岩缝洞储层预测方法和装置 | |
| US20160377752A1 (en) | Method of Digitally Identifying Structural Traps | |
| CN106353807B (zh) | 裂缝识别方法和装置 | |
| US20120099396A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
| CN104111476B (zh) | 构建地层速度场的方法及装置 | |
| CN107991705B (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
| CN113419274B (zh) | 基于高精度层序格架模型的三维地震切片属性体提取方法 | |
| CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 | |
| CN109991663A (zh) | 工区地震速度体校正方法和装置 | |
| CN104199107B (zh) | 基于垂直地震多波数据的钻前深度预测方法以及系统 | |
| Li et al. | Novel strategies for complex foothills seismic imaging—Part 1: Mega-near-surface velocity estimation | |
| Jiang et al. | Quantitative identification of pseudofaults underlying igneous rocks through wave-equation-based velocity modeling: A case study of the Shuntuoguole uplift, Tarim Basin, China | |
| CN114859414B (zh) | 从地震数据中自动提取地层倾角信息的方法和装置 | |
| US11768303B2 (en) | Automatic data enhancement for full waveform inversion in the midpoint-offset domain | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |