本申请属于地球物理勘探技术领域，尤其涉及一种确定地质资源量的数据处理方法及装置。
背景技术：
随着油气勘探工作的逐步深入，准确的评价和计算地质资源量和指明剩余资源量的分布，对于勘探潜力的评价和下一步的勘探方向具有重要的指导作用，是勘探决策的重要依据。
现有技术中，应用较广泛的是采用静态容积法计算地质资源量，该方法是将整个烃源岩层作为有效烃源岩来计算烃源岩的体积，按照单井点单位体积烃源岩生烃量计算地质资源量。利用该方法计算地质资源量的方法考虑的因素比较单一，没有考虑不同区带的差异性，计算结果往往偏大，不够精确，适用于盆地早期勘探阶段，对其他地区尤其是成熟探区的资源评价和勘探方向的决策往往不适用。因此，业内亟需一种能够提高地质资源量计算准确性的实施方案。
技术实现要素：
本申请目的在于提供一种确定地质资源量的数据处理方法及装置，提高了地质资源量计算结果的准确性，为油气的勘探、开发提供了准确的理论基础。
一方面本申请提供了一种确定地质资源量的数据处理方法，包括：
基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；
根据所述地温分布、所述成熟度分布，划分生排烃区带；
根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；
利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
进一步地，所述方法的另一个实施例中，所述根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率，包括：
构建所述成熟度和所述生烃有效率、所述排烃有效率之间的成熟度效率对应关系；
根据所述成熟度效率对应关系和所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
所述成熟度效率关系包括：
上式中，ge表示所述生烃有效率，ee表示所述排烃有效率，ro表示所述成熟度值，ge1表示第一生烃有效率值，ee1表示第一排烃有效率值，ro1表示第一成熟度阈值，ro2表示第二成熟度阈值，ge2表示第二生烃有效率值，ee2表示第二排烃有效率值，ro3表示第三成熟度阈值，ge3表示第三生烃有效率值，ee3表示第三排烃有效率值，ro1＜ro2＜ro3。
进一步地，所述方法的另一个实施例中，所述根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量，包括：
根据所述地温分布，构建所述地温和所述单位体积烃源岩的排烃量之间的温度生烃关系，所述单位体积烃源岩的排烃量包括：单位体积烃源岩的排油量、单位体积烃源岩的排气量；
根据所述温度生烃关系和所述生排烃区带对应的地温，确定所述生排烃区带对应的所述单位体积烃源岩的排气量、所述单位体积烃源岩的排油量；
所述温度生烃关系包括：
上述公式中，voil表示单位体积烃源岩的排油量，v1、v2、v3、v4表示排量系数，t表示所述地温，t1表示第一温度阈值，t2表示第二温度阈值，t3表示第三温度阈值，t1＜t2＜t3。
进一步地，所述方法的另一个实施例中，所述有效烃源岩的厚度的获取方法包括：
利用有机碳含量曲线和地震波阻抗反演技术，计算所述有效烃源岩的厚度。
进一步地，所述方法的另一个实施例中，所述成熟度分布的确定方法包括：
根据公式ro＝r1×d-r2确定所述烃源岩的成熟度与埋深之间的关系，确定所述成熟度分布；
上式中，ro表示所述烃源岩的成熟度，r1、r2表示成熟度系数，d表示所述烃源岩埋深。
进一步地，所述方法的另一个实施例中，所述计算所述生排烃区带对应的地质资源量，包括：
分别计算所述生排烃区带对应的油资源量和气资源量，将所述油资源量和气资源量到的总和作为所述地质资源量；
所述油资源量包括所述生排烃区带的烃源岩的体积、单位体积烃源岩的排油量、生烃有效率、排烃有效率之间的乘积；
所述气资源量包括所述生排烃区带的烃源岩的体积、单位体积烃源岩的排气量、生烃有效率、排烃有效率之间的乘积。
进一步地，所述方法的另一个实施例中，所述方法还包括：
根据所述生排烃区带对应的地质资源量，绘制资源量分布图。
另一方面，本申请提供了一种确定地质资源量的数据处理装置，包括：
温度成熟度确定模块，用于基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；
区带划分模块，用于根据所述地温分布、所述成熟度分布，划分生排烃区带；
有效率确定模块，用于根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
生烃系数确定模块，用于根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；
资源量确定模块，用于利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
又一方面，本申请还提供了一种确定地质资源量的数据处理装置，包括：包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；
根据所述地温分布、所述成熟度分布，划分生排烃区带；
根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；
利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
再一方面，本申请还提供了一种确定地质资源量的数据处理系统，包括：处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现上述确定地质资源量的数据处理方法。
本申请提供的确定地质资源量的数据处理方法及装置，根据待分析地区的烃源岩的埋深，获得待分析地区的地温分布、成熟度分布，根据地温分布以及成熟度分布，对待分析地区进行生排烃区带划分。根据地温与单位体积烃源岩的排烃量之间的关系，可以获得生排烃区带对应的单位体积烃源岩的排烃量，根据成熟度与烃源岩的生排烃有效率之间的关系，可以确定生排烃区带的排烃有效率和生烃有效率。结合生排烃区带的烃源岩的体积、生排烃有效率、排烃量等，可以确定出各个生排烃区带的地质资源量。利用了地温、成熟度对烃源岩的生排烃有效率以及生烃量的影响，针对不同成熟度的烃源岩进行地质资源量的计算，提高了地质资源量的计算结果的准确性，为油气的勘探和开发提供了准确的理论依据。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本申请提供的一种确定地质资源量的数据处理方法一个实施例的方法流程示意图；
图2是本申请实施例中根据钻井数据统计获得的烃源岩的埋深和成熟度之间的关系图；
图3是本申请一个实施例中将烃源岩埋深图转换的烃源岩地温图；
图4是本申请一个实施例中将烃源岩埋深图转换的烃源岩成熟度图；
图5是本申请一个实施例中地温与单位体积烃源岩累积排烃量的关系示意图；
图6是本申请一个实施例中分区带确定的资源量分布示意图；
图7是本申请一个实施例中转换获得的烃源岩排油体积系数示意图；
图8是本申请一个实施例中转换获得的烃源岩排气体积系数示意图；
图9是本申请实施例中toc测井解释模板和g1井toc测井解释成果对比示意图；
图10是本申请一个实施例中烃源岩层段稀疏脉冲反演剖面和转换获得的toc剖面示意图；
图11是本申请实施例中通过地震波阻抗反演技术获得的有效烃源岩的厚度示意图；
图12是本申请提供的确定地质资源量的数据处理装置一个实施例的模块结构示意图；
图13是本申请提供的一种确定地质资源量的数据处理系统实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
地质资源量可以包括可采的油气储量，根据待开发地区或已经开发的地区进行地质资源量的计算，可以为勘探、开发潜力的评价提供一定的理论依据，也可以为勘探或者开发方案提供一定的理论指导。
本申请实施例通过对待分析地区的温度、成熟度进行分析，根据不同区带内的温度以及成熟度之间的差异，确定待分析地区不同区带对应的地质资源量，进一步可以获得整个待分析地区的地质资源量。本申请实施例提供了一种对不同成熟度的烃源岩生排烃有效率进行评价计算地质资源量的方法，能精确计算地质资源量，为勘探、开发决策提供准确的理论依据。
具体地，图1是本申请提供的一种确定地质资源量的数据处理方法一个实施例的方法流程示意图，本申请提供的确定地质资源量的数据处理方法包括：
s1、基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布。
具体地，烃源岩也可以叫生油岩，是一种能够产生或已经产生可移动烃类的岩石。可以根据地质资料、地震解释等数据获得待分析地区的烃源岩的埋深分布。烃源岩的成熟度通常是相对碎屑岩而言的，可以分为结构成熟度和成分成熟度两种。成分成熟度是指砂岩中碎屑组分在风化、搬运、沉积作用的改造下接近最稳定的终极产物的程度。结构成熟度又称物理成熟度，是指碎屑沉积物在其风化，搬运和沉积作用的改造下接近终极结构特征的程度。可以通过分析烃源岩的埋深与烃源岩的地温之间的关系，根据烃源岩的埋深确定烃源岩的地温分布。同样的，可以通过分析烃源岩的埋深与烃源岩的成熟度之间的关系，根据获得的烃源岩的埋深确定烃源岩的成熟度分布。
烃源岩的埋深和地温之间的关系、以及埋深和成熟度之间的关系可以通过分析历史数据或实验数据，利用函数拟合、表格分析等获得。本申请一个实施例中，通过数据分析，烃源岩的埋深和地温之间的关系可以采用如下公式(1)：
t＝d×g+t0(1)
上式中，t可以表示烃源岩段的地温，d可以表示烃源岩埋深，g可以表示盆地地温梯度，t0可以表示地表温度。
在获得待分析地区烃源岩的埋深d后，可以根据上述公式(1)，确定出烃源岩的地温分布t。
图2是本申请实施例中根据钻井数据统计获得的烃源岩的埋深和成熟度之间的关系图，如图2所示，每一条折线图可以表示一口井的烃源岩的埋深与成熟度之间的关系。通过统计分析，可以大致确定烃源岩的埋深和成熟度之间呈直线关系，本申请一个实施例中，通过数据分析，可以采用如下公式(2)计算不同烃源岩的埋深对应的成熟度。
ro＝r1×d-r2(2)
上式中，ro可以表示烃源岩的成熟度，d可以表示烃源岩埋深，r1、r2可以表示成熟度系数，可以根据具体情况进行设置。
当然，不同地区的数据，可能统计获得不同的计算公式，根据实际需要，可以调整烃源岩的埋深与成熟度之间的函数关系式、以及埋深和地温之间的函数关系式，本申请实施例不作具体限定。
本申请一个实施例中，可以将r1的值设置为0.0002，将r2的值设置为0.5903，即烃源岩的埋深和成熟度之间的关系可以采用如下公式(3)进行计算：
ro＝0.0002×d-0.5903(3)
上式中，ro可以表示烃源岩的成熟度，d可以表示烃源岩埋深。式中具体的成熟度系数：0.0002、0.5903，可以根据实际需要进行调整，本申请实施例不做具体限定。
利用公式(1)可以获得不同烃源岩的埋深对应的地温，利用公式(2)或(3)可以计算不同烃源岩的埋深对应的成熟度。本申请实施例中可以根据地质资料、地震解释、测井解释等数据，获得烃源岩的埋深图，根据烃源岩的地温与埋深之间的关系，可以将烃源岩层段的埋深图转换成地温图。图3是本申请一个实施例中将烃源岩埋深图转换的烃源岩地温图，如图3所示，可以在烃源岩埋深图中标示出不同区域对应的地温。还可以根据烃源岩的埋深与成熟度之间的关系，将烃源岩埋深图转换成成熟度图，图4是本申请一个实施例中将烃源岩埋深图转换的烃源岩成熟度图，如图4所示，可以在图中标示出不同区域对应的成熟度。
s2、根据所述地温分布、所述成熟度分布，划分生排烃区带。
具体地，确定出待分析地区的烃源岩的地温分布和成熟度分布后，可以根据实际需要将待分析地区划分出不同的生排烃区带，可以定义地温和成熟度在一定范围内的属于一个生排烃区带。图5是本申请一个实施例中地温与单位体积烃源岩累积排烃量的关系示意图，如：可以将地温小于120℃，成熟度值小于0.62的区域划分为一个生排烃区带，在本申请实施例中图2和图5中，当地温小于120℃时，烃源岩中的有机质不能生成和排出烃，因此生烃有效率和排烃有效率均为0；当地温大于120℃而小于150℃，而成熟度值大于0.62小于1.1时划分为第二个生排烃区带。具体划分的方式可以根据实际需要进行划分，本申请实施例不作具体限定。
s3、根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率。
具体地，生烃可以表示烃源岩生成烃类物质的过程，排烃可以表示烃源岩生烃后，经过初次运移，将烃类物质排驱到输导层(或储集层)中的过程。生烃有效率、排烃有效率可以表示烃源岩中的有机质具体能够生成烃或排出烃所占总有机质的比例。可以预先分析烃源岩的成熟度与生烃有效率、排烃有效率之间的规律，如根据历史数据获得实验数据，通过函数拟合、表格分析等，确定烃源岩的成熟度与生烃有效率、排烃有效率之间的函数关系。根据确定的烃源岩的成熟度与生烃有效率、排烃有效率之间的函数关系，将生排烃区带对应的成熟度值代入确定的函数关系中，可以获得各个生排烃区带对应的生烃有效率、排烃有效率。
s4、根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量。
具体地，可以预先分析烃源岩的地温与单位体积烃源岩的排烃量之间的规律，单位体积烃源岩的排烃量可以表示单位体积烃源岩排出烃类物质的含量，可以包括单位体积烃源岩的排油量和单位体积烃源岩的排气量。如：根据历史数据获得实验数据，通过函数拟合、表格分析等，确定烃源岩的地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间的函数关系。根据确定的烃源岩的地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间的函数关系，将生排烃区带对应的地温代入确定的函数关系中，可以获得各个生排烃区带对应的单位体积烃源岩的排气量、单位体积烃源岩的排油量。
s5、利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
根据划分的生排烃区带的面积以及生排烃区带对应的有效烃源岩的厚度，可以获得生排烃区带对应的体积，结合获得的生排烃区带对应的单位体积烃源岩的排烃量、生烃有效率以及排烃有效率，可以确定各个生排烃区带的地质资源量。
本申请一个实施例中，所述计算所述生排烃区带对应的地质资源量，可以包括：
分别计算所述生排烃区带对应的油资源量和气资源量，将所述油资源量和气资源量到的总和作为所述地质资源量；
所述油资源量包括所述生排烃区带的烃源岩的体积、单位体积烃源岩的排油量、生烃有效率、排烃有效率之间的乘积；
所述气资源量包括所述生排烃区带的烃源岩的体积、单位体积烃源岩的排气量、生烃有效率、排烃有效率之间的乘积。
具体地，可以根据获得的生排烃区带对应的数据，获得生排烃区带的烃源岩的体积以及厚度，计算出生排烃区带的烃源岩的体积，进一步计算各个生排烃区带对应的油资源量和气资源量，将油资源量和气资源量的总和作为生排烃区带的地质资源量。其中，油资源量的计算方式可以采用如下公式(4)：
reso＝h×s×voil×ge×ee(4)
生排烃区带的气资源量的计算方式可以采用如下公式(5)：
resg＝h×s×vgas×ge×ee(5)
各个生排烃区带的地质资源量：res＝reso+resg
上述各式中：reso可以表示油资源量，h可以表示生排烃区带对应的有效烃源岩的厚度，s可以表示生排烃区带对应的面积，voil可以表示生排烃区带对应的单位体积烃源岩的排油量，ge可以表示生排烃区带对应的生烃有效率，ee可以表示生排烃区带对应的排烃有效率，resg可以表示气资源量，vgas可以表示生排烃区带对应的单位体积烃源岩的排气量，res可以表示生排烃区带对应的地质资源量。
本申请一个实施例中，还可以根据计算获得的各个生排烃区带对应的地质资源量，绘制资源量分布图，以便更加直观的对待分析地区的地质资源量进行分析和统计。图6是本申请一个实施例中分区带确定的资源量分布示意图，如图6所示，本申请实施例中可以将待分析地区划分为4个区带，即4个生排烃区带，在每个区带中标注出对应的地质资源量，可以分别将对应的油资源量和气资源量标注出来。图6中g-1，g-3，t-1，t-2，t-3均表示井号，例如：图6中区带3中标注的油：1304，可以表示区带3的油资源量为1304百万桶，气：823，可以表示区带3的气资源量为823百万桶，若没有标注油或气的则可以表示该区带不含油资源量或气资源量。
本申请实施例提供的确定地质资源量的数据处理方法，根据待分析地区的烃源岩的埋深，获得待分析地区的地温分布、成熟度分布，根据地温分布以及成熟度分布，对待分析地区进行生排烃区带划分。根据地温与单位体积烃源岩的排烃量之间的关系，可以获得生排烃区带对应的单位体积烃源岩的排烃量，根据成熟度与烃源岩的生排烃有效率之间的关系，确定生排烃区带的排烃有效率和生烃有效率。结合生排烃区带的烃源岩的体积、生排烃有效率、排烃量等，可以确定出各个生排烃区带的地质资源量。利用了地温、成熟度对烃源岩的生排烃有效率以及生烃量的影响，针对不同成熟度的烃源岩进行地质资源量的计算，提高了地质资源量的计算结果的准确性，为油气的勘探和开发提供了准确的理论依据。
在上述实施例的基础上，本申请一个实施例中，所述根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率，可以包括：
构建所述成熟度和所述生烃有效率、所述排烃有效率之间的成熟度效率对应关系；
根据所述成熟度效率对应关系和所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
所述成熟度效率关系包括：
上式中，ge可以表示所述生烃有效率，ee可以表示所述排烃有效率，ro可以表示所述成熟度值，ge1可以表示第一生烃有效率值，ee1可以表示第一排烃有效率值，ro1可以表示第一成熟度阈值，ro2可以表示第二成熟度阈值，ge2可以表示第二生烃有效率值，ee2可以表示第二排烃有效率值，ro3可以表示第三成熟度阈值，ge3可以表示第三生烃有效率值，ee3可以表示第三排烃有效率值，ro1＜ro2＜ro3。
具体地，在确定生排烃区带的生烃有效率和排烃有效率时，可以预先构建成熟度和生排烃有效率之间的函数关系，如：可以预先根据成熟度与生烃有效率、排烃有效率之间的关系，利用历史数据或实验数据，通过函数拟合、表格分析等数据处理方法，构建出成熟度与生烃有效率、排烃有效率之间的函数关系(即成熟度效率关系)。将生排烃区带对应的成熟度值代入预先构建的成熟度效率关系式中，即可以获得生排烃区带对应的生烃有效率和拍烃有效率。
成熟度与生烃有效率、排烃有效率之间的成熟度效率关系可以表示为上述公式(6)，公式(6)中的第一生烃有效率值ge1、第一排烃有效率值ee1、第一成熟度阈值ro1、第二成熟度阈值ro2、第二生烃有效率值ge2、第二排烃有效率值ee2、第三成熟度阈值ro3、第三生烃有效率值ge3、第三排烃有效率值ee3，可以通过历史数据或实验数据分析或专家经验进行确定。本申请一个示例中，可以将各个参数的值设置为：ge1＝20％，ee1＝10％，ro1＝0.62，ro2＝1.1，ge2＝90％，ee2＝50％，ro3＝1.9，ge3＝100％，ee3＝70％。即上述公式(6)可以表示为：
当成熟度值等于第一成熟度阈值ro1或第二成熟度阈值ro2或第三成熟度阈值ro3时，可以根据实际需要取与阈值相关的任意一个区间内的值，如：若成熟度值等于第一成熟度阈值ro1，则可以将生烃有效率、排烃有效率均确定为0，也可以将生烃有效率确定为第一生烃有效率值ge1、排烃有效率确定为第一排烃有效率值ee1，可以根据实际需要进行选择。
本申请提供的确定地质资源量的数据处理方法，通过对烃源岩的成熟度与生烃有效率、排烃有效率之间关系进行分析，确定了成熟度与生烃有效率、排烃有效率之间的成熟度效率关系。利用成熟度效率关系可以确定各个生排烃区带对应的生烃有效率、排烃有效率，考虑了不同成熟度对烃源岩的生烃、排烃有效率的影响，进一步提高了地质资源量计算结果的准确性。
在上述实施例的基础上，本申请一个实施例中，所述根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量，可以包括：
根据所述地温分布，构建所述地温和所述单位体积烃源岩的排烃量之间的温度生烃关系，所述单位体积烃源岩的排烃量包括：单位体积烃源岩的排油量、单位体积烃源岩的排气量；
根据所述温度生烃关系和所述生排烃区带对应的地温，确定所述生排烃区带对应的所述单位体积烃源岩的排气量、所述单位体积烃源岩的排油量；
所述温度生烃关系包括：
上述公式中，voil可以表示单位体积烃源岩的排油量，v1、v2、v3、v4可以表示排量系数，t可以表示所述地温，t1可以表示第一温度阈值，t2可以表示第二温度阈值，t3可以表示第三温度阈值，t1＜t2＜t3。
具体地，在确定生排烃区带对应的单位体积烃源岩的排气量、单位体积烃源岩的排油量时，可以预先根据地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间的关系，获得温度生烃关系。如：可以通过对历史数据、实验数据等的分析，利用函数拟合、表格分析等，确定出地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间的温度生烃关系。再将生排烃区带的地温代入到温度生烃关系中，可以确定出各个生排烃区带的单位体积烃源岩的排气量、单位体积烃源岩的排油量。
如图5所示，根据历史数据或实验数据进行统计分析，可以确定地温与单位体积烃源岩的排油量和排气量的关系。本申请一个实施例中，通过数据分析等，根据图5可以确定出的温度生烃关系可以参考上述公式(8)，公式(8)中的排量系数v1、v2、v3、v4，第一温度阈值t1，第二温度阈值t2，第三温度阈值t3可以预先根据历史数据、实验数据等分析确定。本申请一个示例中，上述各个参数可以确定为：v1＝0.649、v2＝73.335、v3＝0.8356、v4＝119.11、t1＝120℃、t2＝150℃、t3＝230℃。即上述公式(8)可以表示为：
当地温等于第一温度阈值t1，第二温度阈值t2，第三温度阈值t3时，可以根据实际需要取与阈值相关的任意一个区间内的值，如：若地温等于第一温度阈值t1，则可以将单位体积烃源岩的排油量确定为0，也可以将单位体积烃源岩的排油量确定为voil＝v1×t-v2，可以根据实际需要进行选择，相应地，单位体积烃源岩的排气量可以确定为0。
本申请一个实施例中，可以根据地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间关系，将根据烃源岩埋深获得的烃源岩地温图，转换成烃源岩排油体积系数图和排气体积系数图。图7是本申请一个实施例中转换获得的烃源岩排油体积系数示意图，图8是本申请一个实施例中转换获得的烃源岩排气体积系数示意图。
本申请提供的确定地质资源量的数据处理方法，通过对烃源岩的地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间关系进行分析，确定了地温与单位体积烃源岩的排气量、单位体积烃源岩的排油量之间的温度生烃关系，利用温度生烃关系可以确定各个生排烃区带对应的单位体积烃源岩的排气量、单位体积烃源岩的排油量。考虑了不同地温对烃源岩的排油量、排气量的影响，进一步提高了地质资源量计算结果的准确性。
在上述实施例的基础上，本申请一个实施例中，所述生排烃区带的烃源岩的体积的获取方法可以包括：
利用有机碳含量曲线和地震波阻抗反演技术，计算所述有效烃源岩的厚度；
将所述生排烃区带的面积和所述厚度的乘积作为所述生排烃区带的烃源岩的体积。
具体地，可以利用测井解释，获得声波时差曲线、电阻率曲线、伽玛曲线等，利用passey博士提出的公式确定烃源岩总有机碳含量toc与声波时差和电阻率的关系式，获得toc曲线。图9是本申请实施例中toc测井解释模板和g1井toc测井解释成果对比示意图，图9中左侧是toc测井解释模板，右侧是本申请实施例中g1井toc测井解释成果，图9中gr可以表示伽马曲线，api可以表示伽马曲线值的单位。利用测井解释资料和passey公式，可以获得总有机碳含量toc曲线。本申请一个实施例中依据g1井烃源岩热解地化分析数据确定出的passey公式可以表示为：
toc＝4.12×﹛log﹙r/r基线值﹚+0.02﹙δt-δt基线值﹚+0.25﹜(10)
上式中，r可以表示电阻率值，△t可以表示声波时差值，参数r基线值可以取7.7661，δt基线值可以取1.4216，当然，根据实际需要，参数r基线值和δt基线值也可以取其他数值，本申请实施例不作具体限定。
通过地震波阻抗反演技术如：地震稀疏脉冲阻抗反演，可以得到阻抗数据体。利用反演获得的阻抗数据体，利用阻抗与上述公式(10)中电阻率值r、以及声波时差值△t的关系，可以确定出电阻率值r、以及声波时差值△t，进而转换获得烃源岩有机碳数据体toc。根据反演获得的烃源岩有机碳数据体，利用数学统计等数据分析方法，设定有机碳含量门限值，可以在烃源岩有机碳数据体中确定出有效烃源岩的厚度。本申请一个实施例中，有机碳含量门限值可以设置为toc>1％。图10是本申请一个实施例中烃源岩层段稀疏脉冲反演剖面和转换获得的toc剖面示意图，图10中最左侧的纵坐标可以表示色度，可以表示有机碳含量的变化，另一个纵坐标可以表示烃源岩的深度。图中黑色虚线可以表示有效烃源岩层段(即toc>有机碳含量门限值对应的烃源岩，如toc>1％的烃源岩层段)。通过统计大于有机碳含量门限值的烃源岩的厚度，可以确定出有效烃源岩的厚度。
图11是本申请实施例中通过地震波阻抗反演技术获得的有效烃源岩的厚度示意图，利用确定出的有效烃源岩的厚度、生排烃区带的面积(即生排烃区带的烃源岩的面积)，可以确定出生排烃区带有效烃源岩的体积，再结合生烃有效率、排烃有效率、单位体积烃源岩的排油量、单位体积烃源岩的排气量，可以准确的计算出各个生排烃区带的地质资源量。
下面结合具体的示例介绍本申请实施例的方案：
步骤一、利用测井资料的声波时差曲线、电阻率曲线、伽马曲线依据passey公式确定烃源岩总有机碳含量toc与声波时差和电阻率的关系式。本申请一个示例中，依据g1井烃源岩热解地化分析数据确定passey公式可以表示为上述公式(10)。
步骤二、通过地震稀疏脉冲阻抗反演得到阻抗数据体，利用passey公式转换为烃源岩有机碳数据体toc。
步骤三、按照有机碳含量门限值(toc>1％)在烃源岩有机碳数据体toc中确定有效烃源岩的厚度h。
步骤四、按照地温梯度与烃源岩埋深之间的关系，获得烃源岩的地温分布。本申请一个示例中，地温梯度与烃源岩的埋深之间的关系可以表示为：
t＝d×3.85℃/100m+26.6℃
上式中，t可以表示烃源岩的地温，d可以表示烃源岩的埋深，26.6℃可以表示地表温度。
步骤五、依据烃源岩的埋深d与成熟ro关系，确定出烃源岩的成熟度分布。烃源岩的埋深和成熟度之间的关系式可以表示为：
ro＝0.0002×d-0.5903
步骤六、依据烃源岩的地温t与单位体积烃源岩的排油量和排气量的关系，确定出单位体积烃源岩的排油量和排气量，具体可以参考上述公式(9)。
步骤七、依据不同成熟度值和地温划分生排烃区带，并统计各个生排烃区带的面积s。计算各个生排烃区带的生烃有效率ge和排烃有效率ee，具体可以参考上述公式(7)计算各生排烃区带对应的生烃有效率ge和排烃有效率ee。
步骤八、计算盆地资源量res。可以依据上述公式(4)-(5)分别计算出生排烃区带对应的油资源量和气资源量，可以将油资源量和气资源量的和作为生排烃区带的地质资源量。
需要说明的是，上述各实施例中的公式仅仅是一种示意性的表达公式，可以根据实际需要对各个公式进行调整和变形，本申请实施例不作具体的限定。
本申请提供的确定地质资源量的数据处理方法，综合考虑了地温和成熟度的差异对烃源岩的生排烃效率、以及生烃量的影响，提拱了一种准确计算待分析地区的地质资源量的方法。提高了地质资源量计算结果的准确性，进一步为后续油气的勘探、开发提供了准确的理论依据。
基于上述所述的确定地质资源量的数据处理方法，本说明书一个或多个实施例还提供一种确定地质资源量的数据处理装置。所述的装置可以包括使用了本说明书实施例所述方法的系统(包括分布式系统)、软件(应用)、模块、组件、服务器、客户端等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
具体地，图12是本申请提供的确定地质资源量的数据处理装置一个实施例的模块结构示意图，如图12所示，本申请中提供的确定地质资源量的数据处理装置包括：温度成熟度确定模块121，区带划分模块122，有效率确定模块123，生烃系数确定模块124、资源量确定模块125。
温度成熟度确定模块121，可以用于基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；
区带划分模块122，可以用于根据所述地温分布、所述成熟度分布，划分生排烃区带；
有效率确定模块123，可以用于根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
生烃系数确定模块124，可以用于根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；
资源量确定模块125，可以用于利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
本申请提供的确定地质资源量的数据处理装置，根据待分析地区的烃源岩的埋深，获得待分析地区的地温分布、成熟度分布，根据地温分布以及成熟度分布，对待分析地区进行生排烃区带划分。根据地温与单位体积烃源岩的排烃量之间的关系，可以获得生排烃区带对应的单位体积烃源岩的排烃量，根据成熟度与烃源岩的生排烃有效率之间的关系，确定生排烃区带的排烃有效率和生烃有效率。结合生排烃区带的烃源岩的体积，可以确定出各个生排烃区带的地质资源量。利用了地温、成熟度对烃源岩的生排烃有效率以及生烃量的影响，针对不同成熟度的烃源岩进行地质资源量的计算，提高了地质资源量的计算结果的准确性，为油气的勘探和开发提供了准确的理论依据。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种确定地质资源量的数据处理装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；
根据所述地温分布、所述成熟度分布，划分生排烃区带；
根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；
根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；
利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如ram、rom等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、u盘；利用光学方式存储信息的装置如，cd或dvd。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书实施例提供的上述确定地质资源量的数据处理方法或装置可以在计算机中由处理器执行相应的程序指令来实现，如使用windows操作系统的c++语言在pc端实现、linux系统实现，或其他例如使用android、ios系统程序设计语言在智能终端实现，以及基于量子计算机的处理逻辑实现等。本说明书提供的一种确定地质资源量的数据处理系统的一个实施例中，图13是本申请提供的一种确定地质资源量的数据处理系统实施例的模块结构示意图，如图13所示，本申请另一实施例提供的确定地质资源量的数据处理装置可以包括处理器131以及用于存储处理器可执行指令的存储器132，
处理器131和存储器132通过总线133完成相互间的通信；
所述处理器131用于调用所述存储器132中的程序指令，以执行上述各确定地质资源量的数据处理方法实施例所提供的方法，例如包括：基于烃源岩的埋深与所述烃源岩的地温、成熟度之间的关系，分别确定所述烃源岩的地温分布、成熟度分布；根据所述地温分布、所述成熟度分布，划分生排烃区带；根据所述生排烃区带对应的成熟度值，确定所述生排烃区带对应的生烃有效率、排烃有效率；根据所述生排烃区带对应的地温，确定所述生排烃区带对应的单位体积烃源岩的排烃量；利用所述生烃有效率、所述排烃有效率、所述单位体积烃源岩的排烃量、所述生排烃区带的烃源岩的体积，计算所述生排烃区带对应的地质资源量。
需要说明的是说明书上述所述的装置根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本申请中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
本说明书实施例并不局限于必须是符合行业通信标准、标准计算机数据处理和数据存储规则或本说明书一个或多个实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书实施例的可选实施方案范围之内。
在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(programmablelogicdevice，pld)(例如现场可编程门阵列(fieldprogrammablegatearray，fpga))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片pld上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logiccompiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(hardwaredescriptionlanguage，hdl)，而hdl也并非仅有一种，而是有许多种，如abel(advancedbooleanexpressionlanguage)、ahdl(alterahardwaredescriptionlanguage)、confluence、cupl(cornelluniversityprogramminglanguage)、hdcal、jhdl(javahardwaredescriptionlanguage)、lava、lola、myhdl、palasm、rhdl(rubyhardwaredescriptionlanguage)等，目前最普遍使用的是vhdl(very-high-speedintegratedcircuithardwaredescriptionlanguage)与verilog。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
控制器可以按任何适当的方式实现，例如，控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(applicationspecificintegratedcircuit，asic)、可编程逻辑控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：arc625d、atmelat91sam、microchippic18f26k20以及siliconelabsc8051f320，存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
虽然本说明书一个或多个实施例提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或终端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。第一，第二等词语用来表示名称，而并不表示任何特定的顺序。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本发明是参照根据本发明实施例的方法、装置(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flashram)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器(cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储、石墨烯存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitorymedia)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书一个或多个实施例的实施例而已，并不用于限制本本说明书一个或多个实施例。对于本领域技术人员来说，本说明书一个或多个实施例可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在权利要求范围之内。