CN112394408A - 各向异性介质裂缝预测方法及装置 - Google Patents
各向异性介质裂缝预测方法及装置 Download PDFInfo
- Publication number
- CN112394408A CN112394408A CN202011189478.5A CN202011189478A CN112394408A CN 112394408 A CN112394408 A CN 112394408A CN 202011189478 A CN202011189478 A CN 202011189478A CN 112394408 A CN112394408 A CN 112394408A
- Authority
- CN
- China
- Prior art keywords
- crack
- anisotropic medium
- azimuth
- determining
- anisotropic
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000001228 spectrum Methods 0.000 claims abstract description 78
- 238000004590 computer program Methods 0.000 claims description 15
- 230000003595 spectral effect Effects 0.000 claims description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 21
- 230000008569 process Effects 0.000 abstract description 9
- 230000000875 corresponding effect Effects 0.000 description 65
- 230000005012 migration Effects 0.000 description 18
- 238000013508 migration Methods 0.000 description 18
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 239000004215 Carbon black (E152) Substances 0.000 description 10
- 229930195733 hydrocarbon Natural products 0.000 description 10
- 150000002430 hydrocarbons Chemical class 0.000 description 10
- 238000004891 communication Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/626—Physical property of subsurface with anisotropy
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
本申请实施例提供一种各向异性介质裂缝预测方法及装置，方法包括：基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。本申请能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性。
Description
技术领域
本申请涉及地球物理勘探地震资料计算机识别技术领域，属于地球物理勘探数据处理和解释范畴，具体涉及各向异性介质裂缝预测方法及装置。
背景技术
在石油地球物理勘探过程中，为了更好地对裂缝型油气藏进行解释，往往需要根据地震勘探资料对各向异性介质裂缝体进行解释和预测，以便于解释地质目标，更好地预测油气藏。目前对各向异性介质进行预测的方法主要有基于OVT宽方位道集的振幅各向异性检测、频率各向异性检测、速度各向异性检测等，但在一些地区由于资料性噪比等问题往往通过这些特征无法检测各向异性介质，无法对裂缝进行有效地预测。
发明内容
针对现有技术中的问题，本申请提供一种各向异性介质裂缝预测方法及装置，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性。
为解决上述技术问题，本申请提供以下技术方案：
第一方面，本申请提供一种各向异性介质裂缝预测方法，包括：
基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；
根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；
将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
进一步地，所述基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道，包括：
获取目标区域的OVT宽方位道集资料；
基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位；
将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
进一步地，所述根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度，包括：
对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱；
在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱；
基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
进一步地，所述将各个所述方位分别对应的频谱梯度进行椭圆拟合，包括：
将各个所述方位分别对应的频谱梯度以极坐标进行显示；
根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；
其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
进一步地，所述基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果，包括：
将所述椭圆的长轴方位作为各向异性介质的裂缝方向；
根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度；
输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
进一步地，所述输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果，包括：
判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征；
输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
第二方面，本申请提供一种各向异性介质裂缝预测装置，包括：
叠后地震道确定模块，用于基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；
频率梯度确定模块，用于根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；
裂缝预测模块，用于将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
进一步地，所述叠后地震道确定模块用于执行下述内容
获取目标区域的OVT宽方位道集资料；
基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位；
将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
进一步地，所述频率梯度确定模块用于执行下述内容：
对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱；
在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱；
基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
进一步地，所述裂缝预测模块，包括：
椭圆拟合单元，用于执行下述内容：
将各个所述方位分别对应的频谱梯度以极坐标进行显示；
根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；
其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
进一步地，所述裂缝预测模块，还包括：
方向确定单元，用于将所述椭圆的长轴方位作为各向异性介质的裂缝方向；
强度确定单元，用于根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度；
结果输出单元，用于输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
进一步地，所述结果输出单元用于执行下述内容：
判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征；
输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
第三方面，本申请提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的各向异性介质裂缝预测方法。
第四方面，本申请提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的各向异性介质裂缝预测方法。
由上述技术方案可知，本申请提供的一种各向异性介质裂缝预测方法及装置，方法包括：基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性，进而能够应用该各向异性介质裂缝预测结果来提高裂缝型油气藏的解释精度和效率，并能够有效提高应用裂缝型油气藏的解释结果进行油藏开采等处理的准确性及可靠性，并有效节省人力、物力成本。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例中的各向异性介质裂缝预测方法的流程示意图。
图2是本申请实施例中的各向异性介质裂缝预测方法中步骤100的流程示意图。
图3是本申请实施例中的各向异性介质裂缝预测方法中步骤200的流程示意图。
图4是本申请实施例中的各向异性介质裂缝预测方法中步骤300的第一种具体流程示意图。
图5是本申请实施例中的各向异性介质裂缝预测方法中步骤300的第二种具体流程示意图。
图6是本申请实施例中的各向异性介质裂缝预测方法中步骤350的具体流程示意图。
图7是本申请实施例中的各向异性介质裂缝预测装置的结构示意图。
图8是本申请应用实例中的方位角划分示意图。
图9是本申请应用实例中的频谱梯度计算示意图。
图10是本申请应用实例中的椭圆拟合效果图。
图11是本申请应用实例中的裂缝显示效果示意图。
图12是本申请实施例中的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本申请目地在于提供一种使用OVT宽方位道集资料进行分方位频谱梯度各向异性分析进而预测裂缝方法。在该方法中，首先对宽方位OVT道集资料进行方位角划分，不同方位角的道集资料进行叠加。然后使用傅里叶变化对不同方位的叠后地震资料进行分频变换，计算每个方位角特定频率范围内的频谱变化率。最后将不同方位角内的频谱变化率进行椭圆拟合，从而通过椭圆的长轴方向来预测裂缝方向，通过椭圆的长短轴比来预测裂缝强度。解释人员可以通过裂缝的方向和强度来进行油气藏预测，进而指导下一步的油田开发工作。
基于上述内容，本申请提供一种各向异性介质裂缝预测方法的实施例，参见图1，所述各向异性介质裂缝预测方法具体包含有如下内容：
步骤100：基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道。
可以理解的是，偏移距矢量片OVT技术是一种面向宽方位地震资料的处理方法。
步骤200：根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度。
可以理解的是，叠加偏移剖面是在水平叠加剖面的基础上进行偏移，最初发展的目的是解决水平叠加剖面中倾斜反射层的空间错位和绕射波发散引起的横向分辨率低的问题。根据叠加和偏移进行的顺序，偏移可以分为叠前偏移和叠后偏移，叠加偏移剖面属于叠后偏移。
步骤300：将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
从上述描述可知，本申请实施例提供的各向异性介质裂缝预测方法，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性，进而能够应用该各向异性介质裂缝预测结果来提高裂缝型油气藏的解释精度和效率，并能够有效提高应用裂缝型油气藏的解释结果进行油藏开采等处理的准确性及可靠性，并有效节省人力、物力成本。
在一种各向异性介质裂缝预测方法的实施例中，参见图2，所述各向异性介质裂缝预测方法中的步骤100具体包含有如下内容：
步骤110：获取目标区域的OVT宽方位道集资料。
步骤120：基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位。
步骤130：将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
具体来说，将OVT宽方位道集进行方位角划分，由于地震资料各向异性特征的极坐标对称性，将360度全方位道集按180度方位角平均划分为固定的方位角范围。通常为6个方位。将道集资料按照划分的方位进行叠加，形成分方位叠后地震数据，对于地震道中的每一个采样点，叠加方法为通过所有地震道求平均值(n为地震道个数)。
在一种各向异性介质裂缝预测方法的实施例中，参见图3，所述各向异性介质裂缝预测方法中的步骤200还可以具体包含有如下内容：
步骤210：对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱。
步骤220：在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱。
步骤230：基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
具体来说，通过指定的窗口长度，对每个方位的叠后地震数据序列进行短时窗傅里叶变换。统计频谱值，找到最大频谱位置对应频率，以该频率为起始频率，频谱值接近0对应频率为结束频率。从这些频谱对中取出频谱值计算其频谱梯度。
在一种各向异性介质裂缝预测方法的实施例中，参见图4，所述各向异性介质裂缝预测方法中的步骤300具体包含有如下内容：
步骤310：将各个所述方位分别对应的频谱梯度以极坐标进行显示。
步骤320：根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
具体来说，将不同方位的频谱梯度，根据其对应的方位变换为直角坐标系下的离散点，以便于进行椭圆拟合，将拟合的椭圆在二维视图上显示出来。
在一种各向异性介质裂缝预测方法的实施例中，参见图5，所述各向异性介质裂缝预测方法中的步骤300还具体包含有如下内容：
步骤330：将所述椭圆的长轴方位作为各向异性介质的裂缝方向。
步骤340：根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度。
步骤350：输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
在一种各向异性介质裂缝预测方法的实施例中，参见图6，所述各向异性介质裂缝预测方法中的步骤350还可以具体包含有如下内容：
步骤351：判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征。
步骤352：输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
具体来说，将椭圆长轴方向作为裂缝方向，椭圆长短轴的比值作为裂缝强度进行输出，裂缝强度小于0.8的认为没有裂缝特征。
综上所述，本申请实现了一种使用宽方位OVT道集进行裂缝预测的方法，该方法通过不同方位道集资料的频谱梯度的各向异性来预测裂缝，该方法适用于道集资料的振幅、频率、速度各向异性特征不明显时，无法使用常规方法来进行裂缝预测，使用该方法来进行裂缝预测，从而更有效地对裂缝油气藏进行解释和开发。
从软件层面来说，本申请提供一种用于实现所述各向异性介质裂缝预测方法中全部或部分内容的各向异性介质裂缝预测装置的实施例，参见图7，所述各向异性介质裂缝预测装置具体包含有如下内容：
叠后地震道确定模块10，用于基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道。
频率梯度确定模块20，用于根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度。
裂缝预测模块30，用于将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
本申请提供的各向异性介质裂缝预测装置的实施例具体可以用于执行上述实施例中的各向异性介质裂缝预测方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
从上述描述可知，本申请实施例提供的各向异性介质裂缝预测装置，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性，进而能够应用该各向异性介质裂缝预测结果来提高裂缝型油气藏的解释精度和效率，并能够有效提高应用裂缝型油气藏的解释结果进行油藏开采等处理的准确性及可靠性，并有效节省人力、物力成本。
在一种各向异性介质裂缝预测装置的实施例中，所述各向异性介质裂缝预测装置中的叠后地震道确定模块10用于执行下述内容：
步骤110：获取目标区域的OVT宽方位道集资料。
步骤120：基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位。
步骤130：将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
在一种各向异性介质裂缝预测装置的实施例中，所述各向异性介质裂缝预测装置中频率梯度确定模块20用于执行下述内容：
步骤210：对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱。
步骤220：在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱。
步骤230：基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
在一种各向异性介质裂缝预测装置的实施例中，所述各向异性介质裂缝预测装置中的所述裂缝预测模块30具体用于包含有如下内容：
椭圆拟合单元，用于执行下述内容：
将各个所述方位分别对应的频谱梯度以极坐标进行显示；
根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；
其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
在一种各向异性介质裂缝预测装置的实施例中，所述各向异性介质裂缝预测装置中的所述裂缝预测模块30还具体用于包含有如下内容：
子波提取机器学习模型建立模块，用于执行下述内容：
方向确定单元，用于将所述椭圆的长轴方位作为各向异性介质的裂缝方向；
强度确定单元，用于根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度；
结果输出单元，用于输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
在一种各向异性介质裂缝预测装置的实施例中，所述各向异性介质裂缝预测装置中的所述结果输出单元还具体用于执行下述内容：
步骤351：判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征。
步骤352：输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
为进一步说明本方案，本申请还提供一种应用各向异性介质裂缝预测装置实现各向异性介质裂缝预测方法的具体应用实例，首先对宽方位OVT道集资料进行方位角划分，不同方位角的道集资料进行叠加。然后使用傅里叶变化对不同方位的叠后地震资料进行分频变换，计算每个方位角特定频率范围内的频谱变化率。最后将不同方位角内的频谱变化率进行椭圆拟合，从而通过椭圆的长轴方向来预测裂缝方向，通过椭圆的长短轴比来预测裂缝强度。解释人员可以通过裂缝的方向和强度来进行油气藏预测，进而指导下一步的油田开发工作。
具体实施步骤如下：
1)将OVT宽方位道集进行方位角划分，由于地震资料各向异性特征的极坐标对称性，将360度全方位道集按180度方位角平均划分为固定的方位角范围，通常为6个方位。
2)将道集资料按照划分的方位进行叠加，形成分方位叠后地震数据，每一个方位为一个叠后地震道。
3)对每个方位的叠后地震道进行傅里叶变换，分解为不同频率的频谱。
4)取高频部分频谱，计算每个方位内的特定频率范围内的频率梯度。
5)将不同方位的频谱梯度在极坐标展示出来并进行椭圆拟合,计算椭圆的长轴和短轴长度以及长轴方位角。
6)将椭圆长轴方向作为裂缝方向，椭圆长短轴的比值作为裂缝强度进行输出，裂缝强度小于0.8的认为没有裂缝特征。
基于上述内容，本申请应用实体提供的各向异性介质裂缝预测方法的具体实例如下：
1)将OVT宽方位道集进行方位角划分，由于地震资料各向异性特征的极坐标对称性，将360度全方位道集按180度方位角平均划分为固定的方位角范围。通常为6个方位，如图8所示。
2)将道集资料按照划分的方位进行叠加，形成分方位叠后地震数据，对于地震道中的每一个采样点，叠加方法为通过所有地震道求平均值(n为地震道个数)：
3)通过指定的窗口长度γ，对每个方位的叠后地震数据序列
进行短时窗傅里叶变换：
F[n]为输入一维单道地震数据，k为时窗长度,w特定频率，t为该频率对应频谱信息，变换后叠后地震道分解为不同频率的频谱(w，s)。
4)统计频谱值，找到最大频谱位置对应频率，以该频率为起始频率Fstart，频谱值接近0对应频率为结束频率Fend。从这些频谱对中取出频谱值计算其频谱梯度：
拟合方法对应效果如图9所示。
5)将不同方位的频谱梯度G，根据其对应的方位θ变换为直角坐标系下的离散点，以便于进行椭圆拟合：
Pi(x,y)＝(Gi*cos(θ),Gi*sin(θ))
因为拟合的椭圆中心点位于(0，0)，所以设椭圆方程为：
Ax2+Bxy+Cy2+Dx+Ey＝1
根据最小二乘法计算对应的A，B，C，D，E：
椭圆对应的长轴a，短轴b和旋转角度β分别为：
最后根据a，b，β将拟合的椭圆在二维视图上显示出来，如图10所示。
6)将椭圆长轴方向β作为裂缝方向，椭圆长短轴的比值c作为裂缝强度：
c＝a/b
裂缝强度小于0.8的认为没有裂缝特征,以线条形式将裂缝显示在二维平面图上，线条长度为裂缝强度c，线条朝向为裂缝的方位β，最终效果如图11所示。
本申请应用实例提供的各向异性介质裂缝预测方法，使用OVT宽方位道集资料进行裂缝预测，首先将OVT宽方位道集资料进行方位角划分，划分为不同的方位角范围，将不同方位角范围内的地震道进行叠加，将叠后地震道使用傅里叶变化进行频谱分解，然后计算特定频率范围内频谱的变化梯度，最后在极坐标下将频谱梯度进行椭圆拟合，最后根据椭圆的方向和长短轴比值作为裂缝的方向和强度，进而实现裂缝预测。
从硬件层面来说，本申请提供一种用于实现所述各向异性介质裂缝预测方法中的全部或部分内容的电子设备的实施例，所述电子设备具体包含有如下内容：
图12为本申请实施例的电子设备9600的装置构成的示意框图。如图12所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图12是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
在一实施例中，各向异性介质裂缝预测功能可以被集成到中央处理器中。其中，中央处理器可以被配置为进行如下控制：
步骤100：基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道。
可以理解的是，偏移距矢量片OVT技术是一种面向宽方位地震资料的处理方法。
步骤200：根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度。
可以理解的是，叠加偏移剖面是在水平叠加剖面的基础上进行偏移，最初发展的目的是解决水平叠加剖面中倾斜反射层的空间错位和绕射波发散引起的横向分辨率低的问题。根据叠加和偏移进行的顺序，偏移可以分为叠前偏移和叠后偏移，叠加偏移剖面属于叠后偏移。
步骤300：将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
从上述描述可知，本申请实施例提供的电子设备，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性，进而能够应用该各向异性介质裂缝预测结果来提高裂缝型油气藏的解释精度和效率，并能够有效提高应用裂缝型油气藏的解释结果进行油藏开采等处理的准确性及可靠性，并有效节省人力、物力成本。
在另一个实施方式中，各向异性介质裂缝预测装置可以与中央处理器9100分开配置，例如可以将各向异性介质裂缝预测装置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现各向异性介质裂缝预测功能。
如图12所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图12中所示的所有部件；此外，电子设备9600还可以包括图12中没有示出的部件，可以参考现有技术。
如图12所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
本申请的实施例还提供能够实现上述实施例中的各向异性介质裂缝预测方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的执行主体为服务器或客户端的各向异性介质裂缝预测方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤100：基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道。
可以理解的是，偏移距矢量片OVT技术是一种面向宽方位地震资料的处理方法。
步骤200：根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度。
可以理解的是，叠加偏移剖面是在水平叠加剖面的基础上进行偏移，最初发展的目的是解决水平叠加剖面中倾斜反射层的空间错位和绕射波发散引起的横向分辨率低的问题。根据叠加和偏移进行的顺序，偏移可以分为叠前偏移和叠后偏移，叠加偏移剖面属于叠后偏移。
步骤300：将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
从上述描述可知，本申请实施例提供的计算机可读存储介质，能够对各向异性介质裂缝进行有效地预测，并能够有效提高各向异性介质裂缝预测过程的处理效率、自动化程度及准确性，进而能够应用该各向异性介质裂缝预测结果来提高裂缝型油气藏的解释精度和效率，并能够有效提高应用裂缝型油气藏的解释结果进行油藏开采等处理的准确性及可靠性，并有效节省人力、物力成本。
本领域内的技术人员应明白，本发明的实施例可提供为方法、装置、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(装置)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (14)
1.一种各向异性介质裂缝预测方法，其特征在于，包括：
基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；
根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；
将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
2.根据权利要求1所述的各向异性介质裂缝预测方法，其特征在于，所述基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道，包括：
获取目标区域的OVT宽方位道集资料；
基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位；
将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
3.根据权利要求1所述的各向异性介质裂缝预测方法，其特征在于，所述根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度，包括：
对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱；
在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱；
基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
4.根据权利要求1所述的各向异性介质裂缝预测方法，其特征在于，所述将各个所述方位分别对应的频谱梯度进行椭圆拟合，包括：
将各个所述方位分别对应的频谱梯度以极坐标进行显示；
根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；
其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
5.根据权利要求4所述的各向异性介质裂缝预测方法，其特征在于，所述基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果，包括：
将所述椭圆的长轴方位作为各向异性介质的裂缝方向；
根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度；
输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
6.根据权利要求5所述的各向异性介质裂缝预测方法，其特征在于，所述输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果，包括：
判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征；
输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
7.一种各向异性介质裂缝预测装置，其特征在于，包括：
叠后地震道确定模块，用于基于目标区域的OVT宽方位道集资料确定各个方位分别对应的叠后地震道；
频率梯度确定模块，用于根据各个所述方位分别对应的叠后地震道，分别确定各个方位内的预设频率范围内的频率梯度；
裂缝预测模块，用于将各个所述方位分别对应的频谱梯度进行椭圆拟合，并基于对应的拟合结果确定目标区域的各向异性介质裂缝预测结果。
8.根据权利要求7所述的各向异性介质裂缝预测装置，其特征在于，所述叠后地震道确定模块用于执行下述内容：
获取目标区域的OVT宽方位道集资料；
基于各向异性特征的极坐标对称性，对所述OVT宽方位道集资料进行方位角划分，得到对应的多个方位；
将所述OVT宽方位道集资料按照划分的各个所述方位进行叠加，形成由各个所述方位分别对应的叠后地震道组成的分方位叠后地震数据。
9.根据权利要求7所述的各向异性介质裂缝预测装置，其特征在于，所述频率梯度确定模块用于执行下述内容：
对各个所述方位分别对应的叠后地震道进行傅里叶变换，分解为各个不同频率的频谱；
在分解得到的各个频谱中选取频谱值大于频率阈值的频谱作为目标频谱；
基于所述目标频谱确定各个所述方位内的预设频率范围内的频率梯度。
10.根据权利要求7所述的各向异性介质裂缝预测装置，其特征在于，所述裂缝预测模块，包括：
椭圆拟合单元，用于执行下述内容：
将各个所述方位分别对应的频谱梯度以极坐标进行显示；
根据极坐标显示结果对各个所述方位分别对应的频谱梯度进行椭圆拟合，得到对应的拟合结果；
其中，所述拟合结果包括椭圆的长轴长度、短轴长度和长轴方位。
11.根据权利要求10所述的各向异性介质裂缝预测装置，其特征在于，所述裂缝预测模块，还包括：
方向确定单元，用于将所述椭圆的长轴方位作为各向异性介质的裂缝方向；
强度确定单元，用于根据所述椭圆的长轴长度和短轴长度的比值确定所述各向异性介质的裂缝强度；
结果输出单元，用于输出包含有所述裂缝方向和裂缝强度的各向异性介质裂缝预测结果。
12.根据权利要求11所述的各向异性介质裂缝预测装置，其特征在于，所述结果输出单元用于执行下述内容：
判断所述裂缝强度是否等于或大于强度阈值，若是，则确定各向异性介质存在裂缝特征；
输出所述各向异性介质裂缝预测结果，其中，该各向异性介质裂缝预测结果中包含有裂缝方向、裂缝强度和用于表示各向异性介质存在裂缝特征的信息。
13.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至6任一项所述的各向异性介质裂缝预测方法。
14.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至6任一项所述的各向异性介质裂缝预测方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011189478.5A CN112394408A (zh) | 2020-10-30 | 2020-10-30 | 各向异性介质裂缝预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011189478.5A CN112394408A (zh) | 2020-10-30 | 2020-10-30 | 各向异性介质裂缝预测方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112394408A true CN112394408A (zh) | 2021-02-23 |
Family
ID=74599037
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011189478.5A Pending CN112394408A (zh) | 2020-10-30 | 2020-10-30 | 各向异性介质裂缝预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112394408A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114019565A (zh) * | 2021-11-02 | 2022-02-08 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6473696B1 (en) * | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN102788994A (zh) * | 2012-07-12 | 2012-11-21 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种储层裂缝的确定方法 |
| US20160025880A1 (en) * | 2014-07-24 | 2016-01-28 | Conocophillips Company | Target-oriented process for estimating fracture attributes from seismic data |
| CN105403917A (zh) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | 一种裂缝检测的方法及装置 |
| CN107526103A (zh) * | 2017-08-04 | 2017-12-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地震资料处理方法及其阙值和有效信号频率的求取方法 |
| CN108562936A (zh) * | 2018-03-28 | 2018-09-21 | 北京博达瑞恒科技有限公司 | 一种裂缝预测方法、系统、存储介质及终端 |
| CN109212592A (zh) * | 2017-07-01 | 2019-01-15 | 中国石油化工股份有限公司 | 一种叠前方位p波各向异性裂缝检测方法及装置 |
| WO2020140803A1 (zh) * | 2019-01-02 | 2020-07-09 | 中国石油天然气股份有限公司 | 高角度裂缝预测方法、计算机设备及计算机可读存储介质 |
-
2020
- 2020-10-30 CN CN202011189478.5A patent/CN112394408A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6473696B1 (en) * | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| CN102788994A (zh) * | 2012-07-12 | 2012-11-21 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种储层裂缝的确定方法 |
| US20160025880A1 (en) * | 2014-07-24 | 2016-01-28 | Conocophillips Company | Target-oriented process for estimating fracture attributes from seismic data |
| CN105403917A (zh) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | 一种裂缝检测的方法及装置 |
| CN109212592A (zh) * | 2017-07-01 | 2019-01-15 | 中国石油化工股份有限公司 | 一种叠前方位p波各向异性裂缝检测方法及装置 |
| CN107526103A (zh) * | 2017-08-04 | 2017-12-29 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地震资料处理方法及其阙值和有效信号频率的求取方法 |
| CN108562936A (zh) * | 2018-03-28 | 2018-09-21 | 北京博达瑞恒科技有限公司 | 一种裂缝预测方法、系统、存储介质及终端 |
| WO2020140803A1 (zh) * | 2019-01-02 | 2020-07-09 | 中国石油天然气股份有限公司 | 高角度裂缝预测方法、计算机设备及计算机可读存储介质 |
Non-Patent Citations (2)
| Title |
|---|
| XIANGWEN LI ET AL.: "Pre-stack Anisotropy Fractures Detection Based on OVT Domain Gather Data and Its Application in Ultra-deep Burial Carbonate Rocks", EARTH SCIENCES, 22 October 2020 (2020-10-22) * |
| 孙成田 等: "多属性分析及断裂识别技术在酒泉盆地油气勘探中的应用", 石油地球物理勘探, 31 December 2017 (2017-12-31) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114019565A (zh) * | 2021-11-02 | 2022-02-08 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
| CN114019565B (zh) * | 2021-11-02 | 2023-10-13 | 中海石油(中国)有限公司 | 快速叠前裂缝预测方法、装置、计算机存储介质及设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112415591A (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
| CN112394408A (zh) | 各向异性介质裂缝预测方法及装置 | |
| CN104481526B (zh) | 一种双交叉偶极子测井方法 | |
| CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
| CN109597964A (zh) | 理论油藏渗流参数值确定方法、系统、设备、可读介质 | |
| US20160320502A1 (en) | Electroformed nickel-chromium alloy | |
| CN112147688B (zh) | 一种储层油气识别方法和装置 | |
| CN111983679B (zh) | 基于沉积参数的相控型白云岩储层地震预测方法和装置 | |
| CN112255679B (zh) | 地震资料绕射深度偏移处理方法及装置 | |
| US20150153466A1 (en) | Source Start Time Determination | |
| CN107942373B (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
| US20150117150A1 (en) | In-phase precursor seismic sweep | |
| CN112817042B (zh) | 沉积储层中油气识别方法及装置 | |
| CA2693247A1 (en) | Seismic attributes for structural analysis | |
| CN110632650B (zh) | 全波形反演数据的分时窗提取方法及装置 | |
| CN114114392A (zh) | 层速度模型建立方法及装置 | |
| CN112130209A (zh) | 岩溶储层预测方法及装置 | |
| CN112731519A (zh) | 确定颤动时间区间的方法和装置 | |
| CN114114391B (zh) | 地震数据负面样本识别方法及装置 | |
| CN110824555A (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
| CN117849875B (zh) | 一种地震信号分析方法、系统、装置及存储介质 | |
| CN112596102B (zh) | 油藏渗流特征参数的地球物理预测方法、装置及介质 | |
| CN113075747A (zh) | 储层裂缝发育区域的预测方法及装置 | |
| CN112859171B (zh) | 利用有效裂缝进行储层连通性分析的方法和装置 | |
| CN107060741B (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 |