CN111812712A - 储层中目标层段的油气类型确定方法及系统 - Google Patents
储层中目标层段的油气类型确定方法及系统 Download PDFInfo
- Publication number
- CN111812712A CN111812712A CN201910289789.XA CN201910289789A CN111812712A CN 111812712 A CN111812712 A CN 111812712A CN 201910289789 A CN201910289789 A CN 201910289789A CN 111812712 A CN111812712 A CN 111812712A
- Authority
- CN
- China
- Prior art keywords
- frequency
- interval
- low
- amplitude value
- relative amplitude
- 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 55
- 238000009826 distribution Methods 0.000 claims abstract description 189
- 238000001228 spectrum Methods 0.000 claims abstract description 111
- 238000010586 diagram Methods 0.000 claims abstract description 106
- 230000004044 response Effects 0.000 claims abstract description 97
- 238000001514 detection method Methods 0.000 claims abstract description 41
- 239000004215 Carbon black (E152) Substances 0.000 claims description 112
- 229930195733 hydrocarbon Natural products 0.000 claims description 112
- 150000002430 hydrocarbons Chemical class 0.000 claims description 108
- 230000003595 spectral effect Effects 0.000 claims description 60
- 238000012545 processing Methods 0.000 claims description 41
- 238000010183 spectrum analysis Methods 0.000 claims description 20
- 238000010219 correlation analysis Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 238000012216 screening Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 abstract description 6
- 239000012071 phase Substances 0.000 description 67
- 239000011435 rock Substances 0.000 description 14
- 230000002902 bimodal effect Effects 0.000 description 10
- 238000004891 communication Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 125000001183 hydrocarbyl group Chemical group 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- -1 then Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 239000007790 solid phase Substances 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
-
- 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
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本申请提供一种储层中目标层段的油气类型确定方法及系统，通过进行油气检测获得能量分布图，然后与预设的频谱响应特征标准相进行比对，频谱响应特征标准相是根据同区域内的完井数据得到，这样根据每个区域的特性产生一个适应该区域油气属性的参考标准，后续的油气检测可以直接与该参考标准进行比对，一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
Description
技术领域
本申请涉及石油勘探技术领域，更具体的，涉及一种储层中目标层段的油气类型确定方法及系统。
背景技术
烃类检测技术种类较多，实用频率较高的方法是通过识别低频共振高频衰减现象，对油气层进行预测和识别，它建立在双相介质理论之上，即岩石骨架(即固相)与孔隙流体(即流相)。在应用此技术过程中，不同地区、油气田受构造形态、火成岩干扰、储层结构等诸多因素影响，在利用地震资料开展频谱分析后，仍无法准确确定目的层的含油气性，例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，因此目前的频谱分析技术无法准确标定油气发育区域。
发明内容
为了解决目前基于双相介质理论的通过识别低频共振高频衰减技术无法准确标定油气发育区域的问题，本申请提供一种储层中目标层段的油气类型确定方法及系统。
本申请第一方面实施例提供一种储层中目标层段的油气类型确定方法，包括：
基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
在某些实施例中，还包括：
根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
在某些实施例中，所述根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相，包括：
采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
在某些实施例中，所述设定频率区间包括低频敏感频率区间和高频敏感频率区间，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图和对应高频敏感频率区间的高频能量分布图；
基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型，包括：
对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
在某些实施例中，所述根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型，包括：
计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
在某些实施例中，所述根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型，包括：
计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
在某些实施例中，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
本申请第二方面实施例提供一种储层中目标层段的油气类型确定系统，包括：
油气检测模块，基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
油气类型确定模块，基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
在某些实施例中，还包括：
频谱响应特征标准相建立模块，根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
在某些实施例中，所述频谱响应特征标准相建立模块，包括：
层位标定单元，采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
频谱时窗确定单元，将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
频谱特征获得单元，截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
频谱特征筛选单元，从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
频谱响应特征标准相确定单元，根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
在某些实施例中，所述设定频率区间包括低频敏感频率区间和高频敏感频率区间，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图和对应高频敏感频率区间的高频能量分布图；
所述油气类型确定模块，包括：
商差分布图获取单元，对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
相对振幅值范围确定单元，对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
油气类型确定单元，根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
在某些实施例中，所述油气类型确定单元，包括：
比值差值范围计算单元，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
比对单元，将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
第一求和处理确定单元，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
在某些实施例中，所述油气类型确定单元，包括：
比值差值范围计算单元，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
相关性分析单元，对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
第二求和处理确定单元，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
在某些实施例中，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
本申请第三方面实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的储层中目标层段的油气类型确定方法的步骤。
本申请第四方面实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的储层中目标层段的油气类型确定方法的步骤。
本申请的有益效果如下：
本申请提供一种储层中目标层段的油气类型确定方法及系统，通过进行油气检测获得能量分布图，然后与预设的频谱响应特征标准相进行比对，频谱响应特征标准相是根据同区域内的完井数据得到，这样根据每个区域的特性产生一个适应该区域油气属性的参考标准，后续的油气检测可以直接与该参考标准进行比对，一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本申请实施例中储层中目标层段的油气类型确定方法的流程示意图之一。
图2示出本申请实施例中储层中目标层段的油气类型确定方法的流程示意图之二。
图3示出本申请实施例中图2的步骤S001的具体流程示意图。
图4示出本申请实施例中地震数据示意图。
图5示出本申请实施例中地震剖面上标定出频谱时窗的示意图。
图6示出本申请实施例中频谱响应特征标准相的示意图。
图7a示出本申请实施例中低频能量分布图。
图7b示出本申请实施例中高频能量分布图。
图7c示出本申请实施例中低频高频能量之商分布图。
图7d示出本申请实施例中低频高频能量之差分布图。
图8示出本申请实施例中图1或图2中步骤S200的具体流程示意图。
图9示出本申请实施例中图8中步骤S203的具体流程示意图之一。
图10示出本申请实施例中图8中步骤S203的具体流程示意图之二。
图11示出本申请实施例中储层中目标层段的油气类型确定系统的结构示意图之一。
图12示出本申请实施例中储层中目标层段的油气类型确定系统的结构示意图之二。
图13示出本申请实施例中图12中频谱响应特征标准相建立模块001的具体结构示意图。
图14示出本申请实施例中图11或图12中油气类型确定模块200的具体结构示意图。
图15示出本申请实施例中图14中油气类型确定单元203的具体结构示意图之一。
图16示出本申请实施例中图14中油气类型确定单元203的具体结构示意图之二。
图17示出适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
不同地区、油气田受构造形态、火成岩干扰、储层结构等诸多因素影响，在利用地震资料开展频谱分析后，仍无法准确确定目的层的含油气性，例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，因此目前的频谱分析技术无法准确标定油气发育区域。
基于此，本申请第一方面实施例提供一种储层中目标层段的油气类型确定方法，如图1所示，包括：
S100:基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
S200:基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
本方面首先通过进行油气检测获得能量分布图，然后与预设的频谱响应特征标准相进行比对，频谱响应特征标准相是根据同区域内的完井数据得到，这样根据每个区域的特性产生一个适应该区域油气属性的参考标准，后续的油气检测可以直接与该参考标准进行比对，这样一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
预设频谱响应特征标准相可以是预建立的，即实现建立完成的，也可以包括在本方案中，鉴于此，如图2所示，上述方法还包括：
S001:根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
频谱响应特征标准相时基于多个完井不同油气类型的对应层段的频谱响应特征建立的，因此在建立频谱响应特征标准相时需要进行完整的频谱分析，获得频谱响应特征，然后根据频谱响应特征绘制标准相。
如图3所示，建立频谱响应特征标准相的步骤具体包括：
S001a：采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
S001b：将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
S001c：截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
S001d：从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
S001e：根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
具体的，下面以一实例对上述建立步骤进行详细说明。
首先要收集整理研究区(目标区域)的完钻井的地震资料(如图4所示)和生产数据，确定各井的出油、出气、油水同层等层段的起始深度。通过制作合成记录，将以上层段表示在地震剖面上(如图5所示)，读取层段的顶底界面对应的地震时间，作为油气频谱分析的时窗(图5中的虚线段)。这样就得到以上所说的不同含油气级别层段对应的频谱特征，但同级别含油气层段的频谱特征经统计后，占优势且具有代表性的频谱特征被优选出来，根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相，如图6所示。
本申请基于高频衰减和低频共振的双相理论，因此，本申请的所述设定频率区间包括低频敏感频率区间和高频敏感频率区间。
进一步的，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
在一实施例中，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图(如图7a所示)和对应高频敏感频率区间的高频能量分布图(如图7b所示)。
如图8所示，步骤S200具体包括：
S201:对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
S202:对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
S203:根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
具体的，低频高频能量之商分布图如图7c所示，低频高频能量之差分布图如图7d所示。
步骤S203包括两种不同方案，其一为计算匹配度(比对)，其二为计算相关度(相关性分析)。
对于计算匹配度的技术方案，如图9所示，步骤S203具体包括：
S231a:计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
S232a:将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
S233a:根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
匹配度可以与比对结果具有对应关系，作为一种示例，例如可以根据两个比对范围的重合区域确定匹配度，假设[1-100]，与[2-101]，重合的区域为[2-100]，即两者重合度近乎于99％(近似值)，可以认为匹配度为99％。
此外，对于相关度的技术方案，如图10所示，步骤S203具体包括：
S231b:计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
S232b:对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
S233b:根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
相关度可以根据相关性分析的计算公式计算，例如通过协方差公式计算相关度，本申请不予赘述。
上述步骤可以通过输入油气检测软件实施，只要设定好油气检测软件的参数即可，例如，前期采用人工进行识别和对应的参数设置，后续只要类似输入，即可获得类似的结果，以图7a至图7d为例，人工识别时，图7a至图7d同时具有“三大一小”特征，即低频能量大(图7a)、高频能量小(图7b)、能量之商大(图7c)、能量之差大(图7d)的范围内的构造部位相对较高的层段，可作为油气勘探开发部署的首选目标。
本申请在建立标准相之后，按照标准相中油层、气层的形态将1～10HZ定义为低频敏感段，15～35hz为高频敏感段，其间为鞍部。敏感段确定之后，在油气检测软件分别输入低、高频敏感段，这样就可以得到针对油层和气层而开展的油气检测，同时也避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性。(举个例子：由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，故由标准相所提供的敏感段范围和形态标准，对于识别油气层就会成为一种更加可靠的油气检测手段。取得标准相后，不需要再按照提取已知油气层段频谱提取的方法去进行分析，这样不仅盲目，工作量(标定、地震解释顶底界面、砂体分布未知)巨大且难度较大。故利用油气检测软件，输入上述敏感段、对大套层位或大时间段内的地震资料开展油气检测，就会得到对油气识别性和针对性更强的检测结果，利用其油气检测平面图(即上述的能量分布图)可以在平面上快速确定“低频共振、高频衰减”的优势区域，结合构造图，优先部署于符合“低频共振、高频衰减”规律的构造高部位，在实际应用中取得非常好的效果。
此外，由于油气检测输入的低、高频起始值是基于油层和气层敏感段数值的输入，所以图7a上的低频共振就是油层、气层引起的低频共振，图7b上的高频衰减就是油层、气层引起的高频衰减，同时具有上述这两种特征的层段就是油气部署的有利区带。
实际工作过程中，在已经建立标准相后，对地震上的背斜、楔状体等多种油气可能聚集的某构造部位或地质体开展频谱提取，并将该频谱特征与标准相进行对比，快速确定某构造部位或地质体的频谱特征与标准相中哪一种油气级别地层的频谱特征相似，然后再开展大范围的油气检测，确定构造部位或地质体的平面分布范围，并优选构造高点进行部署。
通过上述实施例的详细描述可以知晓，本方面一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
基于与上述方法相同的发明构思，本申请第二方面实施例提供一种储层中目标层段的油气类型确定系统，如图11所示，包括：
油气检测模块100，基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
油气类型确定模块200，基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
本方面首先通过进行油气检测获得能量分布图，然后与预设的频谱响应特征标准相进行比对，频谱响应特征标准相是根据同区域内的完井数据得到，这样根据每个区域的特性产生一个适应该区域油气属性的参考标准，后续的油气检测可以直接与该参考标准进行比对，这样一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
预设频谱响应特征标准相可以是预建立的，即实现建立完成的，也可以包括在本方案中，鉴于此，如图12所示，上述系统还包括：
频谱响应特征标准相建立模块001，根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
频谱响应特征标准相时基于多个完井不同油气类型的对应层段的频谱响应特征建立的，因此在建立频谱响应特征标准相时需要进行完整的频谱分析，获得频谱响应特征，然后根据频谱响应特征绘制标准相。
如图13所示，所述频谱响应特征标准相建立模块001，包括：
层位标定单元001a，采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
频谱时窗确定单元001b，将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
频谱特征获得单元001c，截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
频谱特征筛选单元001d，从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
频谱响应特征标准相确定单元001e，根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
具体的，下面以一实例对上述建立步骤进行详细说明。
首先要收集整理研究区(目标区域)的完钻井的地震资料(如图4所示)和生产数据，确定各井的出油、出气、油水同层等层段的起始深度。通过制作合成记录，将以上层段表示在地震剖面上(如图5所示)，读取层段的顶底界面对应的地震时间，作为油气频谱分析的时窗(图5中的虚线段)。这样就得到以上所说的不同含油气级别层段对应的频谱特征，但同级别含油气层段的频谱特征经统计后，占优势且具有代表性的频谱特征被优选出来，根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相，如图6所示。
本申请基于高频衰减和低频共振的双相理论，因此，本申请的所述设定频率区间包括低频敏感频率区间和高频敏感频率区间。
进一步的，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
在一实施例中，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图(如图7a所示)和对应高频敏感频率区间的高频能量分布图(如图7b所示)。
如图14所示，所述油气类型确定模块200，包括：
商差分布图获取单元201，对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
相对振幅值范围确定单元202，对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
油气类型确定单元203，根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
具体的，低频高频能量之商分布图如图7c所示，低频高频能量之差分布图如图7d所示。
油气类型确定单元可以通过计算匹配度(比对)确定油气类型，也可以通过计算相关度(相关性分析)确定油气类型。
对于计算匹配度的技术方案，如图15所示，所述油气类型确定单元203，包括：
比值差值范围计算单元231a，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
比对单元232a，将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
第一求和处理确定单元233a，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
匹配度可以与比对结果具有对应关系，作为一种示例，例如可以根据两个比对范围的重合区域确定匹配度，假设[1-100]，与[2-101]，重合的区域为[2-100]，即两者重合度近乎于99％(近似值)，可以认为匹配度为99％。
此外，对于相关度的技术方案，如图16所示，所述油气类型确定单元203，包括：
比值差值范围计算单元231b，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
相关性分析单元232b，对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
第二求和处理确定单元233b，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
油气类型确定模块可以是具有处理能力的处理设备，例如处理器等，其上设置有油气检测软件，例如只要设定好油气检测软件的参数即可，例如，前期采用人工进行识别和对应的参数设置，后续只要类似输入，即可获得类似的结果，以图7a至图7d为例，人工识别时，图7a至图7d同时具有“三大一小”特征，即低频能量大(图7a)、高频能量小(图7b)、能量之商大(图7c)、能量之差大(图7d)的范围内的构造部位相对较高的层段，可作为油气勘探开发部署的首选目标。
本申请在建立标准相之后，按照标准相中油层、气层的形态将1～10HZ定义为低频敏感段，15～35hz为高频敏感段，其间为鞍部。敏感段确定之后，在油气检测软件分别输入低、高频敏感段，这样就可以得到针对油层和气层而开展的油气检测，同时也避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性。(举个例子：由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，故由标准相所提供的敏感段范围和形态标准，对于识别油气层就会成为一种更加可靠的油气检测手段。取得标准相后，不需要再按照提取已知油气层段频谱提取的方法去进行分析，这样不仅盲目，工作量(标定、地震解释顶底界面、砂体分布未知)巨大且难度较大。故利用油气检测软件，输入上述敏感段、对大套层位或大时间段内的地震资料开展油气检测，就会得到对油气识别性和针对性更强的检测结果，利用其油气检测平面图(即上述的能量分布图)可以在平面上快速确定“低频共振、高频衰减”的优势区域，结合构造图，优先部署于符合“低频共振、高频衰减”规律的构造高部位，在实际应用中取得非常好的效果。
此外，由于油气检测输入的低、高频起始值是基于油层和气层敏感段数值的输入，所以图7a上的低频共振就是油层、气层引起的低频共振，图7b上的高频衰减就是油层、气层引起的高频衰减，同时具有上述这两种特征的层段就是油气部署的有利区带。
实际工作过程中，在已经建立标准相后，对地震上的背斜、楔状体等多种油气可能聚集的某构造部位或地质体开展频谱提取，并将该频谱特征与标准相进行对比，快速确定某构造部位或地质体的频谱特征与标准相中哪一种油气级别地层的频谱特征相似，然后再开展大范围的油气检测，确定构造部位或地质体的平面分布范围，并优选构造高点进行部署。
通过上述实施例的详细描述可以知晓，本方面一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种电子设备的具体实施方式，参见图17，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
S200:基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型。
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
从上述描述可知，本申请提供的电子设备，一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
本申请的实施例还提供能够实现上述实施例中的方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
S100:基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
S200:基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型。
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
从上述描述可知，本申请提供的计算机可读存储介质，一方面能够避虑掉因其它非油气层的频率分布对油气层识别的干扰，进而提高准确性(例如由于薄互层干扰、火成岩干扰等影响，即使干层也可能表现出类似油气层双峰特征的频谱分布特征，导致常规频谱特征无法识别出油气层和干层)，另一方面，预设频谱响应特征标准相后，在后续的油气检测时不需要按照提取已知油气层段频谱提取的方法去进行分析，油气检测工作量少，难度低。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。虽然本说明书实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书实施例时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。以上所述仅为本说明书实施例的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (16)
1.一种储层中目标层段的油气类型确定方法，其特征在于，包括：
基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
2.根据权利要求1所述的油气类型确定方法，其特征在于，还包括：
根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
3.根据权利要求2所述的油气类型确定方法，其特征在于，所述根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相，包括：
采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
4.根据权利要求3所述的油气类型确定方法，其特征在于，所述设定频率区间包括低频敏感频率区间和高频敏感频率区间，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图和对应高频敏感频率区间的高频能量分布图；
基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型，包括：
对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
5.根据权利要求4所述的油气类型确定方法，其特征在于，所述根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型，包括：
计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
6.根据权利要求4所述的油气类型确定方法，其特征在于，所述根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型，包括：
计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
7.根据权利要求4所述的油气类型确定方法，其特征在于，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
8.一种储层中目标层段的油气类型确定系统，其特征在于，包括：
油气检测模块，基于目标区域的地震数据，对所述目标区域中的待测钻井开展油气检测，获取所述目标区域中的目标层段在设定频率区间的地震波能量分布图；
油气类型确定模块，基于预设的频谱响应特征标准相，根据所述地震波能量分布图确定所述目标层段的油气类型；
其中，所述频谱响应特征标准相是根据目标区域中多个完井的不同油气类型的对应层段的频谱响应特征得到。
9.根据权利要求8所述的油气类型确定系统，其特征在于，还包括：
频谱响应特征标准相建立模块，根据目标区域中多个完井不同油气类型的对应层段的频谱响应特征建立所述频谱响应特征标准相。
10.根据权利要求9所述的油气类型确定系统，其特征在于，所述频谱响应特征标准相建立模块，包括：
层位标定单元，采集目标区域的地震资料和测井资料，对各完井进行层位标定，确定目标区域的每种油气类型的对应层段的起始深度；
频谱时窗确定单元，将每种油气类型的对应层段的起始深度标记在目标区域的地震剖面上，并从已标记的地震剖面上读取各层段的顶底界面对应的地震时间点，将每个层段的两个地震时间点之间的时间区间作为一个频谱分析时窗；
频谱特征获得单元，截取每个频谱分析时窗中的地震信号，获得每种油气类型的对应层段的频谱特征；
频谱特征筛选单元，从所有所述频谱特征中筛选出特征误差小于设定阈值的频谱特征；
频谱响应特征标准相确定单元，根据筛选出的频谱特征确定每种油气类型的对应层段的相对振幅值随频率的变化形态，并针对每种油气类型的对应层段一一绘制出相对振幅值与频率的对应关系曲线，其中所有油气类型的对应层段的对应关系曲线共同形成所述频谱响应特征标准相。
11.根据权利要求10所述的油气类型确定系统，其特征在于，所述设定频率区间包括低频敏感频率区间和高频敏感频率区间，所述地震波能量分布图包括：对应低频敏感频率区间的低频能量分布图和对应高频敏感频率区间的高频能量分布图；
所述油气类型确定模块，包括：
商差分布图获取单元，对所述低频能量分布图和高频能量分布图进行做商和做差处理，获得低频高频能量之商分布图和低频高频能量之差分布图；
相对振幅值范围确定单元，对所述低频能量分布图、高频能量分布图、低频高频能量之商分布图和低频高频能量之差分布图进行数值化处理，获得各分布图的分布范围最大的能量区间，并将各分布图的能量区间转化为相对振幅值范围；
油气类型确定单元，根据各分布图对应的相对振幅值范围以及所述频谱响应特征标准相中每一个对应关系曲线在高频敏感区间和低频敏感区间上的相对振幅值范围，确定目标层段的油气类型。
12.根据权利要求11所述的油气类型确定系统，其特征在于，所述油气类型确定单元，包括：
比值差值范围计算单元，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
比对单元，将低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行比对，输出各比对的匹配度，并对每个对应关系曲线所对应的所有匹配度进行求和处理；
第一求和处理确定单元，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应匹配度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
13.根据权利要求11所述的油气类型确定系统，其特征在于，所述油气类型确定单元，包括：
比值差值范围计算单元，计算每个对应关系曲线低频敏感区间与高频敏感区间的相对振幅值的比值范围以及差值范围；
相关性分析单元，对低频能量分布图对应的相对振幅值范围与每个对应关系曲线低频敏感区间的相对振幅值范围、高频能量分布图对应的相对振幅值范围与每个对应关系高频敏感区间的相对振幅值范围、低频高频能量之商分布图对应的相对振幅值范围与每个对应关系曲线对应的所述比值范围，以及低频高频能量之差分布图对应的相对振幅值范围与每个对应关系曲线对应的所述差值范围进行相关性分析，输出各相关性分析的相关度，并对每个对应关系曲线所对应的所有相关度进行求和处理；
第二求和处理确定单元，根据所述求和处理的结果确定所述目标层段的油气类型；其中对应相关度总和最高的对应关系曲线所表示的油气类型为所述目标层段的油气类型。
14.根据权利要求11所述的油气类型确定系统，其特征在于，所述低频敏感频率区间的频率范围为1-10HZ，所述高频敏感频率区间的频率范围为20-35HZ。
15.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至7任一项所述的储层中目标层段的油气类型确定方法的步骤。
16.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至7任一项所述的储层中目标层段的油气类型确定方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910289789.XA CN111812712A (zh) | 2019-04-11 | 2019-04-11 | 储层中目标层段的油气类型确定方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910289789.XA CN111812712A (zh) | 2019-04-11 | 2019-04-11 | 储层中目标层段的油气类型确定方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111812712A true CN111812712A (zh) | 2020-10-23 |
Family
ID=72843672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910289789.XA Pending CN111812712A (zh) | 2019-04-11 | 2019-04-11 | 储层中目标层段的油气类型确定方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111812712A (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (zh) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | 一种基于多相介质理论的油气检测方法 |
| WO2012071115A2 (en) * | 2010-11-23 | 2012-05-31 | Apex Spectral Technology, Inc. | Methods of hydrocarbon detection using spectra dominant frequency and measures of energy decay on the low side and high side of spectra dominant frequency |
| CN103645499A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 基于叠后反射波能量统计的地表一致性振幅补偿方法 |
-
2019
- 2019-04-11 CN CN201910289789.XA patent/CN111812712A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (zh) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | 一种基于多相介质理论的油气检测方法 |
| WO2012071115A2 (en) * | 2010-11-23 | 2012-05-31 | Apex Spectral Technology, Inc. | Methods of hydrocarbon detection using spectra dominant frequency and measures of energy decay on the low side and high side of spectra dominant frequency |
| CN103645499A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 基于叠后反射波能量统计的地表一致性振幅补偿方法 |
Non-Patent Citations (1)
| Title |
|---|
| 于亮亮: "利用频谱特征检测油气—以辽河坳陷大平房地区为例", 《石油地球物理勘探》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108122066B (zh) | 储层岩性的确定方法和装置 | |
| CN109709607B (zh) | 一种预测薄层砂岩储层厚度方法及装置 | |
| EP1984759A2 (en) | Methods and apparatus for predicting the hydrocarbon production of a well location | |
| EP1570300B1 (en) | Log permeability model calibration using reservoir fluid flow measurements | |
| CN111399056B (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
| US11879316B2 (en) | Geostatistical analysis of microseismic data in fracture modeling | |
| CN109425900B (zh) | 一种地震储层预测方法 | |
| CN107589469B (zh) | 油水界面的确定方法和装置 | |
| CN107894615A (zh) | 一种定量化评价三维地震属性预测储层参数有效性的方法 | |
| CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
| CN104011566A (zh) | 用于分析地质构造的特性的基于小波变换的系统和方法 | |
| CN112363242A (zh) | 基于测录井融合的储层流体识别方法与装置 | |
| EP2761332B1 (en) | System and method for analyzing seismic data | |
| Pyrcz et al. | Declustering and debiasing | |
| US10746888B2 (en) | Microseismic density mapping | |
| CN111812712A (zh) | 储层中目标层段的油气类型确定方法及系统 | |
| CN108181655B (zh) | 暗河垮塌体系形态的确定方法和装置 | |
| CN107991705B (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
| CN111983671A (zh) | 基于微古地貌恢复的浅水湖盆储层预测方法和装置 | |
| CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 | |
| CN111679318B (zh) | 储层中断裂系统的识别方法及系统 | |
| RU2119583C1 (ru) | Способ контроля за разработкой нефтяных залежей | |
| CN107435537B (zh) | 一种录井目标层自动划分的方法及装置 | |
| CN110795513B (zh) | 河流相源储异位型致密油气甜点区分布的预测方法 | |
| CN110308488B (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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201023 |