油页岩反射能量确定方法、系统、存储介质及电子设备与流程

1.本发明属于地质勘探技术领域，尤其涉及一种油页岩反射能量确定方法、系统、存储介质及电子设备。


背景技术：

2.由于油页岩、煤层、火成岩等地质体在地震剖面上出现的强反射响应特征，会严重影响弱反射目标储层的识别和研究。在现有技术中，一般是利用匹配追踪等方法，将强反射能量从储层能量中进行分离，以更好地研究弱反射目标储层的反射能量。但是在油页岩反射能量分析方面，目前并没有相关的研究。在油页岩反射能量分析中使用的大部分方法一般是通过参数测试或者人工经验来确定油页岩的能量关系。但是在依靠人工经验来进行参数选取会出现较大的误差，从而导致油页岩能量分析缺乏相应依据和可靠性，影响了油页岩能量分离的效果和对后续储层预测的精度。因此，需要提出一种新的油页岩反射能量定量分析方法，从而提高能量分析的可靠性。


技术实现要素：

3.本发明正是基于现有技术中通过参数测试或者人工经验来确定油页岩的能量关系，会导致油页岩能量分析缺乏相应依据和可靠性，影响了油页岩能量分离的效果和对后续储层预测的精度的技术问题，提出了一种油页岩反射能量确定方法、系统、存储介质及电子设备。
4.第一方面，本发明实施例提供了一种油页岩反射能量确定方法，包括：
5.获取待研究工区的测井资料；
6.基于所述测井资料确定所述待研究工区的全反射系数序列；
7.从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列；
8.分别对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形；
9.根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比；
10.根据所述百分比以及所述待研究工区的储层总能量的估计值分析所述待研究工区的油页岩反射能量。
11.可选地，从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列，包括：
12.基于油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度；
13.基于所述油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度，从所述全反射系数序列中提取在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层
段的能量反射强度的反射系数，获得油页岩反射系数序列。
14.可选地，所述测井资料包括声波测井数据以及密度测井数据。
15.可选地，基于所述测井资料确定所述待研究工区的全反射系数序列，包括：
16.根据所述声波测井数据，计算地层声波速度，获得速度场数据；
17.根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列；
18.基于时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得所述待研究工区的全反射系数序列。
19.可选地，根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列，包括：
20.根据所述速度场数据以及所述密度测井数据，利用预设计算式计算得到初始反射系数序列；其中，所述预设计算式为：
[0021][0022]
其中，r(i)表示初始反射系数序列中第i个采样点的反射系数，ρ(i)表示密度测井数据中的第i个采样点的密度数据，ρ(i 1)表示密度测井数据中的第i 1个采样点的密度数据，v(i)表示速度场数据中的第i个采样点的速度场数据，v(i 1)表示速度场数据中的第i 1个采样点的速度场数据。
[0023]
可选地，根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，包括：
[0024]
将所述油页岩反射波形与所述全反射波形进行叠合；
[0025]
计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0026]
第二方面，本发明实施例还提供了一种油页岩反射能量确定系统，包括：
[0027]
获取模块，用于获取待研究工区的测井资料；
[0028]
全反射系数序列计算模块，用于基于所述测井资料确定所述待研究工区的全反射系数序列；
[0029]
油页岩反射系数序列提取模块，用于从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列；
[0030]
地震正演模拟模块，用于分别对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形；
[0031]
油页岩反射能量百分比确定模块，用于根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比；
[0032]
分析模块，用于根据所述百分比以及所述待研究工区的储层总能量的估计值分析所述待研究工区的油页岩反射能量。
[0033]
可选地，所述油页岩反射系数序列提取模块，包括：
[0034]
时窗确定单元，用于基于油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度；
[0035]
提取单元，用于基于所述油页岩发育层段所在的时窗以及油页岩发育层段的能量
反射强度，从所述全反射系数序列中提取在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，获得油页岩反射系数序列。
[0036]
可选地，所述测井资料包括声波测井数据以及密度测井数据。
[0037]
可选地，所述全反射系数序列计算模块包括：
[0038]
速度场数据获取单元，用于根据所述声波测井数据，计算地层声波速度，获得速度场数据；
[0039]
初始反射系数序列计算单元，用于根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列；
[0040]
时深转换单元，基于时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得所述待研究工区的全反射系数序列。
[0041]
可选地，所述初始反射系数序列计算单元具体用于：
[0042]
根据所述速度场数据以及所述密度测井数据，利用预设计算式计算得到初始反射系数序列；其中，所述预设计算式为：
[0043][0044]
其中，r(i)表示初始反射系数序列中第i个采样点的反射系数，ρ(i)表示密度测井数据中的第i个采样点的密度数据，ρ(i 1)表示密度测井数据中的第i 1个采样点的密度数据，v(i)表示速度场数据中的第i个采样点的速度场数据，v(i 1)表示速度场数据中的第i 1个采样点的速度场数据。
[0045]
可选地，所述油页岩反射能量百分比确定模块包括：
[0046]
叠合单元，用于将所述油页岩反射波形与所述全反射波形进行叠合；
[0047]
油页岩反射能量计算单元，用于计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0048]
第三方面，本发明实施例还提供了一种存储介质，所述存储介质上存储有程序代码，所述程序代码被处理器执行时，实现如上述实施例中任一项所述的油页岩反射能量确定方法。
[0049]
第四方面，本发明实施例还提供了一种电子设备，所述电子设备包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的程序代码，所述程序代码被所述处理器执行时，实现如上述实施例中任一项所述的油页岩反射能量确定方法。
[0050]
在本发明实施例提供的一种油页岩反射能量确定方法、系统、存储介质及电子设备，通过利用待研究工区的测井资料计算该待研究工区的全反射系数序列，进而从该全反射系数序列中提取出与油页岩特征相关的反射系数，可以获得在油页岩发育时窗内的油页岩反射系数序列，从而对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形，进而能够根据油页岩反射波形以及全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，实现对所述待研究工区的储层能量进行分析。可见，本发明实施例提供的油页岩反射能量确定方法，能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强
反射特征分离提供有效依据，同时也具有很强的适用性。
附图说明
[0051]
通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：
[0052]
图1示出了本发明实施例一提出的一种油页岩反射能量确定方法的流程示意图；
[0053]
图2示出了本发明实施例提出的全反射系数序列的示意图；
[0054]
图3示出了本发明实施例提出的提取油页岩反射系数序列的示意图；
[0055]
图4示出了本发明实施例提出的地震正演模拟的示意图；
[0056]
图5示出了本发明实施例提出的油页岩能量分析的示意图。
具体实施方式
[0057]
为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。
[0058]
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
[0059]
实施例一
[0060]
根据本发明的实施例，提供了一种油页岩反射能量确定方法，图1示出了本发明实施例一提出的一种油页岩反射能量确定方法的流程示意图，如图1所示，该油页岩反射能量确定方法可以包括：步骤110至步骤160。
[0061]
在步骤110中，获取待研究工区的测井资料。
[0062]
这里，测井是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法，属于应用地球物理方法之一。在石油钻井时，在钻到设计井深深度后都必须进行测井，又称完井电测，以获得各种石油地质及工程技术资料。其中，测井资料包括电阻率测井资料、孔隙度测井资料、自然电位测井资料以及自然伽马测井资料等，通过对测井资料进行分析，可以了解待研究工区储层评价情况以及综合地质情况等信息。
[0063]
值得说明的是，待研究工区的测井资料可以是利用在该待研究工区中的一口钻井获取到的测井资料，也可以是利用该待研究工区中的多口钻井获取到的测井资料。
[0064]
在步骤120中，基于所述测井资料确定所述待研究工区的全反射系数序列。
[0065]
这里，通过对测井资料进行计算，即可获得待研究工区的全反射系数序列。其中，全反射系数序列实质上是一个数据集，其包括该从测井资料中获得的多个反射系数。
[0066]
值得说明的是，反射系数的获取是建立在速度资料和密度资料的解释的基础上的，由此可以得出波阻抗曲线，最后利用计算式计算出反射系数。其中，速度资料可以由速度测井资料获取，密度资料可以从密度测井中获得。如果研究区没有密度资料，则可以假定密度不变。
[0067]
图2示出了本发明实施例提出的全反射系数序列的示意图，如图2所示，通过对测
井资料进行计算，即可获得如图2所示的全反射系数序列。
[0068]
在步骤130中，从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列。
[0069]
这里，全反射系数序列是通过测井获得的该待研究地区的完整反射系数序列，其包含了不同储层段的反射系数。通过从所述全反射系数序列中提取与油页岩特征相关的反射系数，可以获得该待研究地区的油页岩反射系数序列。
[0070]
其中，油页岩对不同测井曲线响应的敏感性研究表明，油页岩在测井曲线上一般表现为高电阻率、高自然伽马、高声波时差和低密度的特征。因此，该油页岩反射系数序列实质上是包含了在油页岩发育时窗内的反射系数，其可以用来对油页岩发育层段进行分析。
[0071]
在步骤140中，分别对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形。
[0072]
这里，通过对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形。
[0073]
其中，地震正演模拟是根据地质体的形状、产状和物性数据，通过构造数学模型计算得到其理论值，或通过构造实体模型来观测模型所产生的地球物理效应的数值。通过地震正演模拟获得的油页岩反射波形体现了在该待研究工区中油页岩能量的情况，而通过地震正演模拟获得的全反射波形体现了在该待研究工区中所有能量的情况，其不仅包含了具有强反射特征的油页岩能量，也包含其他低反射特征的储层能量。
[0074]
在步骤150中，根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比。
[0075]
在步骤160中，根据所述百分比以及所述待研究工区的储层总能量的估计值分析所述待研究工区的油页岩反射能量。
[0076]
这里，通过计算所述油页岩反射波形在所述全反射波形中的百分比，就可以确定到在该待研究工区中油页岩反射能量所占的比重，而根据该比重就可以在给后续储层能量分析提供依据。
[0077]
在本实施例中，通过利用待研究工区的测井资料计算该待研究工区的全反射系数序列，进而从该全反射系数序列中提取出与油页岩特征相关的反射系数，可以获得在油页岩发育时窗内的油页岩反射系数序列，从而对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形，进而能够根据油页岩反射波形以及全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，实现对所述待研究工区的储层能量进行分析。能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强反射特征分离提供有效依据，同时也具有很强的适用性。
[0078]
实施例二
[0079]
在上述实施例的基础上，本发明的实施例二还可以提供一种油页岩反射能量确定方法。该油页岩反射能量确定方法可以包括：步骤210至步骤260。
[0080]
在步骤210中，获取待研究工区的测井资料。
[0081]
这里，测井是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法，属于应用地球物理方法之一。在石油钻井时，在钻到设计井深深度后都必须进行测井，又称完井电测，以获得各种石油地质及工程技术资料。其中，测井资料包括电阻率测井资料、孔隙度测井资料、自然电位测井资料以及自然伽马测井资料等，通过对测井资料进行分析，可以了解待研究工区储层评价情况以及综合地质情况等信息。
[0082]
值得说明的是，待研究工区的测井资料可以是利用在该待研究工区中的一口钻井获取到的测井资料，也可以是利用该待研究工区中的多口钻井获取到的测井资料。
[0083]
其中，所述测井资料包括声波测井数据以及密度测井数据。
[0084]
这里，声波测井数据则是利用声波测井获得的声波测井数据，根据岩石的声学物理特性来测量地层声波速度。其主要用来确定岩石孔隙度、划分地层、判断岩性以及判断气层等。密度测井数据是利用密度测井获得密度数据，密度测井的目的是测量地层的体积密度，包括岩石骨架密度和流体密度。其可以用来确定地层孔隙度、计算矿物含量、曲线重叠确定岩性等。
[0085]
在步骤220中，基于所述测井资料确定所述待研究工区的全反射系数序列。
[0086]
这里，通过对测井资料进行计算，即可获得待研究工区的全反射系数序列。其中，全反射系数序列实质上是一个数据集，其包括该从测井资料中获得的多个反射系数。
[0087]
值得说明的是，反射系数的获取是建立在速度资料和密度资料的解释的基础上的，由此可以得出波阻抗曲线，最后利用计算式计算出反射系数。其中，速度资料可以由速度测井资料获取，密度资料可以从密度测井中获得。如果研究区没有密度资料，则可以假定密度不变。
[0088]
在一个可选的实施方式中，步骤220中，基于所述测井资料确定所述待研究工区的全反射系数序列，包括步骤221至步骤223。
[0089]
步骤221，根据所述声波测井数据，计算地层声波速度，获得速度场数据；
[0090]
步骤222，根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列；
[0091]
步骤223，基于时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得所述待研究工区的全反射系数序列。
[0092]
这里，时深转换是将地震数据从时间域向深度域转换的过程，它依靠速度，向人们展示了地下的构造情况。时深转换就是对地震数据处理、速度分析、测井资料研究几个过程的不断迭代往复，以优化转换结果。声波测井数据、校验炮观测和垂直地震剖面都可以帮助优化转换结果，还可以促进测井数据、钻井数据和地震数据的关联。其中，时深转换关系可以根据测井数据，利用反射率法得到时移地震差异反演中的正演算子；然后根据时移地震差异反演中的正演算子，利用测井数据来正演模拟角度道集进行井震标定，与井旁地震道进行相关对比，确定时深转换关系。
[0093]
步骤221至步骤223实际过程是：首先利用声波测井数据得到声波在岩石中的传播速度，即获得速度场数据，再根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列。该初始反射系数序列是属于深度域的反射系数序列，因此，需要根据时深转换
关系将所述初始反射系数序列从深度域转换至时间域，获得处于时间域的全反射系数序列。
[0094]
其中，步骤221至步骤223的过程实际上是利用密度测井数据代表的密度曲线以及声波测井数据代表的声波时差曲线来计算反射系数。
[0095]
在一个可选的实施方式中，步骤222中，根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列，包括：
[0096]
根据所述速度场数据以及所述密度测井数据，利用预设计算式计算得到初始反射系数序列；其中，所述预设计算式为：
[0097][0098]
其中，r(i)表示初始反射系数序列中第i个采样点的反射系数，ρ(i)表示密度测井数据中的第i个采样点的密度数据，ρ(i 1)表示密度测井数据中的第i 1个采样点的密度数据，v(i)表示速度场数据中的第i个采样点的速度场数据，v(i 1)表示速度场数据中的第i 1个采样点的速度场数据。
[0099]
这里，通过该预设计算式对依次采样的采样点计算，即可获得该待研究工区的完整的反射系数，从而获得全反射系数序列。
[0100]
在步骤230中，从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列。
[0101]
这里，全反射系数序列是通过测井获得的该待研究地区的完整反射系数序列，其包含了不同储层段的反射系数。通过从所述全反射系数序列中提取与油页岩特征相关的反射系数，可以获得该待研究地区的油页岩反射系数序列。
[0102]
其中，油页岩对不同测井曲线响应的敏感性研究表明，油页岩在测井曲线上一般表现为高电阻率、高自然伽马、高声波时差和低密度的特征。因此，该油页岩反射系数序列实质上是包含了在油页岩发育时窗内的反射系数，其可以用来对油页岩发育层段进行分析。
[0103]
在一个可选的实施方式中，步骤230中，从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列，可以包括步骤231至步骤232。
[0104]
在步骤231中，基于油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度；
[0105]
在步骤232中，基于所述油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度，从所述全反射系数序列中提取在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，获得油页岩反射系数序列。
[0106]
这里，由于在全反射系数序列中不仅包含了具有强反射特征的油页岩反射能量，也包含了其他储层的低反射特征的能量。因此，先通过油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度，然后再从全反射系数序列中分离出在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，从而将分离出的反射系数作为油页岩反射系数序列。
[0107]
图3示出了本发明实施例提出的提取油页岩反射系数序列的示意图，如图3所示，trace1是全反射系数序列，trace2是提取的油页岩反射系数序列，trace3是非油页岩反射
系数系列。
[0108]
其中，非油页岩反射系数系列是全反射系数序列减去油页岩反射系数序列，该非油页岩反射系数系列中并不含有油页岩放射系数，其代表的是储层段的反射系数序列。
[0109]
在步骤240中，分别对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形。
[0110]
这里，通过对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形。
[0111]
其中，地震正演模拟是根据地质体的形状、产状和物性数据，通过构造数学模型计算得到其理论值，或通过构造实体模型来观测模型所产生的地球物理效应的数值。通过地震正演模拟获得的油页岩反射波形体现了在该待研究工区中油页岩能量的情况，而通过地震正演模拟获得的全反射波形体现了在该待研究工区中所有能量的情况，其不仅包含了具有强反射特征的油页岩能量，也包含其他低反射特征的储层能量。
[0112]
图4示出了本发明实施例提出的地震正演模拟的示意图，如4所示，trace1是对全反射系数序列进行地震正演模拟获得的全反射波形，trace2是由油页岩反射系数序列进行地震正演模拟获得的油页岩反射波形，trace3是对非油页岩反射系数系列进行正演模拟获得的非油页岩反射波形。
[0113]
在步骤250中，根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比。
[0114]
这里，通过计算所述油页岩反射波形在所述全反射波形中的百分比，就可以确定到在该待研究工区中油页岩反射能量所占的比重，而根据该比重就可以在给后续储层能量分析提供依据。
[0115]
在一个实施方式中，步骤250中，根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，可以包括步骤251至步骤252。
[0116]
步骤251，将所述油页岩反射波形与所述全反射波形进行叠合；
[0117]
步骤252，计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0118]
这里，图5示出了本发明实施例提出的油页岩能量分析的示意图，如图5所示，将图4中的全反射波形与油页岩反射波形进行叠合得到如图5所示的油页岩能量分析示意图，然后计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0119]
在步骤260中，根据所述百分比以及所述待研究工区的储层总能量的估计值分析所述待研究工区的油页岩反射能量。
[0120]
从图5可知，该油页岩反射波形对应的油页岩反射能量占该全反射波形在该油页岩发育时窗内的总能量的80％左右，则非油页岩反射波形代表的反射能量占该全反射波形在该油页岩发育时窗内的总能量的20％左右，这个对油页岩反射能量的定量化认识为后续储层研究奠定基础。
[0121]
在本实施例中，通过利用待研究工区的测井资料计算该待研究工区的全反射系数
序列，进而从该全反射系数序列中提取出与油页岩特征相关的反射系数，可以获得在油页岩发育时窗内的油页岩反射系数序列，从而对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形，进而能够根据油页岩反射波形以及全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，实现对所述待研究工区的储层能量进行分析。能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强反射特征分离提供有效依据，同时也具有很强的适用性。
[0122]
实施例三
[0123]
在上述实施例的基础上，本实施例通过具体应用案例来对上述实施例中所述的油页岩反射能量确定方法进行说明。
[0124]
本实施例对鄂尔多斯盆地某探区中的油页岩反射能量进行了确定。
[0125]
在鄂尔多斯盆地某探区，其储层段上部广泛发育油页岩，在地震剖面上显示为强反射特征，该强反射特征对下伏储层的识别带来了很大的影响。因此，如何确定油页岩反射能量的强弱，成为后续储层预测的一个关键环节。
[0126]
首先，利用鄂尔多斯盆地某探区中的待研究工区的一口钻井，获得该待研究工区的测井资料。其中，测井资料包括声波测井数据以及密度测井数据。
[0127]
进而通过测井数据获得该工区的全反射系数序列，具体步骤是：
[0128]
通过根据所述声波测井数据，计算地层声波速度，获得速度场数据；
[0129]
根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列；
[0130]
基于时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得所述待研究工区的全反射系数序列。
[0131]
其中，该全反射系数序列如图2所示。
[0132]
全反射系数序列的具体计算过程是：首先利用声波测井数据得到声波在岩石中的传播速度，即获得速度场数据。再根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列。该初始反射系数序列是属于深度域的反射系数序列，因此，需要根据时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得处于时间域的全反射系数序列。
[0133]
其中，根据所述速度场数据以及所述密度测井数据，利用预设计算式计算得到初始反射系数序列；其中，所述预设计算式为：
[0134][0135]
其中，r(i)表示初始反射系数序列中第i个采样点的反射系数，ρ(i)表示密度测井数据中的第i个采样点的密度数据，ρ(i 1)表示密度测井数据中的第i 1个采样点的密度数据，v(i)表示速度场数据中的第i个采样点的速度场数据，v(i 1)表示速度场数据中的第i 1个采样点的速度场数据。
[0136]
这里，通过该预设计算式对依次采样的采样点计算，即可获得该待研究工区的完整的反射系数，从而获得全反射系数序列。
[0137]
由于油页岩段反射系数较强，通过确定确定油页岩发育层段所在的时窗以及油页
岩发育层段的能量反射强度，从全反射系数序列中提取出在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，获得油页岩反射系数序列。
[0138]
其中，该油页岩反射系数序列如图3所示。
[0139]
这里，全反射系数序列是通过测井获得的该待研究地区的完整反射系数序列，其包含了不同储层段的反射系数。通过从所述全反射系数序列中提取与油页岩特征相关的反射系数，可以获得该待研究地区的油页岩反射系数序列。
[0140]
其中，油页岩对不同测井曲线响应的敏感性研究表明，油页岩在测井曲线上一般表现为高电阻率、高自然伽马、高声波时差和低密度的特征。因此，该油页岩反射系数序列实质上是包含了在油页岩发育时窗内的反射系数，其可以用来对油页岩发育层段进行分析。
[0141]
其次，利用地震正演模拟技术对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形。
[0142]
这里，还可以对非油页岩反射波形进行正演模拟，获得如图4所示的3道地震反射波形。
[0143]
由图4可知，油页岩反射能量较强，非油页岩(储层)反射特征较弱。通过波形叠合开展能量对比，可以分析油页岩反射能量与原始地震反射波形的相互关系，具体是：
[0144]
将所述油页岩反射波形与所述全反射波形进行叠合；
[0145]
计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0146]
如图5所示，该油页岩反射波形对应的油页岩反射能量占该全反射波形在该油页岩发育时窗内的总能量的80％左右，则非油页岩反射波形代表的反射能量占该全反射波形在该油页岩发育时窗内的总能量的20％左右，这个对油页岩反射能量的定量化认识为后续储层研究奠定基础。
[0147]
由此，本发明实施例提供的油页岩反射能量确定方法，能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强反射特征分离提供有效依据，同时也具有很强的适用性。
[0148]
实施例四
[0149]
根据本发明的实施例，还提供了一种油页岩反射能量确定系统，包括：
[0150]
获取模块，用于获取待研究工区的测井资料；
[0151]
全反射系数序列计算模块，用于基于所述测井资料确定所述待研究工区的全反射系数序列；
[0152]
油页岩反射系数序列提取模块，用于从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列；
[0153]
地震正演模拟模块，用于分别对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形；
[0154]
油页岩反射能量百分比确定模块，用于根据所述油页岩反射波形以及所述全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比；
[0155]
分析模块，用于根据所述百分比以及所述待研究工区的储层总能量的估计值分析
所述待研究工区的油页岩反射能量。
[0156]
这里，测井是利用岩层的电化学特性、导电特性、声学特性、放射性等地球物理特性，测量地球物理参数的方法，属于应用地球物理方法之一。在石油钻井时，在钻到设计井深深度后都必须进行测井，又称完井电测，以获得各种石油地质及工程技术资料。其中，测井资料包括电阻率测井资料、孔隙度测井资料、自然电位测井资料以及自然伽马测井资料等，通过对测井资料进行分析，可以了解待研究工区储层评价情况以及综合地质情况等信息。
[0157]
其中，通过全反射系数序列计算模块对获取模块采集到的测井资料，计算待研究工区的全反射系数序列。进而通过油页岩反射系数序列提取模块从所述全反射系数序列中提取与油页岩特征相关的反射系数，获得油页岩反射系数序列，然后地震正演模拟模块对所述油页岩反射系数序列以及所述全反射系数序列进行地震正演模拟，获得油页岩反射波形以及全反射波形。
[0158]
这里，通过对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形。
[0159]
其中，地震正演模拟是根据地质体的形状、产状和物性数据，通过构造数学模型计算得到其理论值，或通过构造实体模型来观测模型所产生的地球物理效应的数值。通过地震正演模拟获得的油页岩反射波形体现了在该待研究工区中油页岩能量的情况，而通过地震正演模拟获得的全反射波形体现了在该待研究工区中所有能量的情况，其不仅包含了具有强反射特征的油页岩能量，也包含其他低反射特征的储层能量。
[0160]
因此，通过油页岩反射能量百分比确定模块确定所述油页岩反射波形在所述全反射波形中的百分比。其中，通过计算所述油页岩反射波形在所述全反射波形中的百分比，就可以确定到在该待研究工区中油页岩反射能量所占的比重，而根据该比重就可以在给后续储层能量分析提供依据。
[0161]
可选地，所述油页岩反射系数序列提取模块，包括：
[0162]
时窗确定单元，用于基于油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度；
[0163]
提取单元，用于基于所述油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度，从所述全反射系数序列中提取在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，获得油页岩反射系数序列。
[0164]
这里，由于在全反射系数序列中不仅包含了具有强反射特征的油页岩反射能量，也包含了其他储层的低反射特征的能量。因此，先通过油页岩特征确定油页岩发育层段所在的时窗以及油页岩发育层段的能量反射强度，然后再从全反射系数序列中分离出在所述油页岩发育层段所在的时窗内且符合所述油页岩发育层段的能量反射强度的反射系数，从而将分离出的反射系数作为油页岩反射系数序列。
[0165]
图3示出了本发明实施例提出的提取油页岩反射系数序列的示意图，如图3所示，trace1是全反射系数序列，trace2是提取的油页岩反射系数序列，trace3是非油页岩反射系数系列。
[0166]
其中，非油页岩反射系数系列是全反射系数序列减去油页岩反射系数序列，该非
油页岩反射系数系列中并不含有油页岩放射系数，其代表的是储层段的反射系数序列。
[0167]
可选地，所述测井资料包括声波测井数据以及密度测井数据。
[0168]
可选地，所述全反射系数序列计算模块包括：
[0169]
速度场数据获取单元，用于根据所述声波测井数据，计算地层声波速度，获得速度场数据；
[0170]
初始反射系数序列计算单元，用于根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列；
[0171]
时深转换单元，基于时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得所述待研究工区的全反射系数序列。
[0172]
这里，首先利用声波测井数据得到声波在岩石中的传播速度，即获得速度场数据，再根据所述速度场数据以及所述密度测井数据，计算得到初始反射系数序列。该初始反射系数序列是属于深度域的反射系数序列，因此，需要根据时深转换关系将所述初始反射系数序列从深度域转换至时间域，获得处于时间域的全反射系数序列。
[0173]
可选地，所述初始反射系数序列计算单元具体用于：
[0174]
根据所述速度场数据以及所述密度测井数据，利用预设计算式计算得到初始反射系数序列；其中，所述预设计算式为：
[0175][0176]
其中，r(i)表示初始反射系数序列中第i个采样点的反射系数，ρ(i)表示密度测井数据中的第i个采样点的密度数据，ρ(i 1)表示密度测井数据中的第i 1个采样点的密度数据，v(i)表示速度场数据中的第i个采样点的速度场数据，v(i 1)表示速度场数据中的第i 1个采样点的速度场数据。
[0177]
这里，通过该预设计算式对依次采样的采样点计算，即可获得该待研究工区的完整的反射系数，从而获得全反射系数序列。
[0178]
可选地，所述油页岩反射能量百分比确定模块包括：
[0179]
叠合单元，用于将所述油页岩反射波形与所述全反射波形进行叠合；
[0180]
油页岩反射能量计算单元，用于计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0181]
这里，图5示出了本发明实施例提出的油页岩能量分析的示意图，如图5所示，将图4中的全反射波形与油页岩反射波形进行叠合得到如图5所示的油页岩能量分析示意图，然后计算所述油页岩反射波形在所述全反射波形处于油页岩发育层段所在的时窗内的反射波形中所占的百分比，从而确定所述油页岩反射波形在所述全反射波形中的百分比。
[0182]
从图5可知，该油页岩反射波形对应的油页岩反射能量占该全反射波形在该油页岩发育时窗内的总能量的80％左右，则非油页岩反射波形代表的反射能量占该全反射波形在该油页岩发育时窗内的总能量的20％左右，这个对油页岩反射能量的定量化认识为后续储层研究奠定基础。
[0183]
在本实施例中，通过利用待研究工区的测井资料计算该待研究工区的全反射系数序列，进而从该全反射系数序列中提取出与油页岩特征相关的反射系数，可以获得在油页
岩发育时窗内的油页岩反射系数序列，从而对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形，进而能够根据油页岩反射波形以及全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，实现对所述待研究工区的储层能量进行分析。能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强反射特征分离提供有效依据，同时也具有很强的适用性。
[0184]
实施例五
[0185]
根据本发明的实施例，还提供了一种存储介质，所述存储介质上存储有程序代码，所述程序代码被处理器执行时，实现如上述实施例任一项所述的油页岩反射能量确定方法。
[0186]
其中，存储介质可以是如闪存、硬盘、多媒体卡、卡型存储器(例如，sd或dx存储器等)、随机访问存储器(ram)、静态随机访问存储器(sram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、可编程只读存储器(prom)、磁性存储器、磁盘、光盘、服务器、app应用商城等等。
[0187]
实施例六
[0188]
根据本发明的实施例，还提供了一种电子设备，所述电子设备包括存储器、处理器，所述存储器上存储有可在所述处理器上运行的程序代码，所述程序代码被所述处理器执行时，实现如上述实施例任一项所述的油页岩反射能量确定方法。
[0189]
可以理解，电子设备还可以包括，输入/输出(i/o)接口，以及通信组件。
[0190]
其中，处理器用于执行如述实施例任一项所述的暗河厚度预测方法中的全部或部分步骤。存储器用于存储各种类型的数据，这些数据例如可以包括电子设备中的任何应用程序或方法的指令，以及应用程序相关的数据。
[0191]
所述处理器可以是专用集成电路(application specific integrated circuit，简称asic)、数字信号处理器(digital signal processor，简称dsp)、数字信号处理设备(digital signal processing device，简称dspd)、可编程逻辑器件(programmable logic device，简称pld)、现场可编程门阵列(field programmable gate array，简称fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行述实施例任一项所述的暗河厚度预测方法。
[0192]
所述存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，例如静态随机存取存储器(static random access memory，简称sram)，电可擦除可编程只读存储器(electrically erasable programmable read-only memory，简称eeprom)，可擦除可编程只读存储器(erasable programmable read-only memory，简称eprom)，可编程只读存储器(programmable read-only memory，简称prom)，只读存储器(read-only memory，简称rom)，磁存储器，快闪存储器，磁盘或光盘。
[0193]
以上结合附图详细说明了本发明的技术方案，考虑到相关技术中，通过参数测试或者人工经验来确定油页岩的能量关系，会导致油页岩能量分析缺乏相应依据和可靠性，影响了油页岩能量分离的效果和对后续储层预测的精度。本发明提供一种油页岩反射能量确定方法、系统、存储介质及电子设备，通过利用待研究工区的测井资料计算该待研究工区
的全反射系数序列，进而从该全反射系数序列中提取出与油页岩特征相关的反射系数，可以获得在油页岩发育时窗内的油页岩反射系数序列，从而对该油页岩反射系数序列以及该全反射系数序列进行地震正演模拟，获得油页岩反射系数序列对应的油页岩反射波形以及该全反射系数序列对应的全反射波形，进而能够根据油页岩反射波形以及全反射波形，确定所述油页岩反射波形在所述全反射波形中的百分比，实现对所述待研究工区的储层能量进行分析。可见，本发明实施例提供的油页岩反射能量确定方法，能够根据待研究工区的测井数据来进行油页岩反射能量的分析工作，充分尊重了原始的测井资料的特征，使得能量分析具有较高的可靠性，为后续强反射特征分离提供有效依据，同时也具有很强的适用性。
[0194]
在本技术所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
[0195]
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
[0196]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0197]
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台电子设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-only memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0198]
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。