技术领域
[0001] 本申请涉及石油与天然气勘探技术领域，特别地涉及一种有效烃源岩分布的预测方法、装置、电子设备及存储介质。
相关背景技术
[0002] 烃源岩是指富含有机质、大量生成油气与排出油气的岩石。对于盆地的油气资源评价，首先要开展的研究工作就是盆地烃源岩的评价，同时，盆地烃源岩的评价结果也是衡量盆地油气资源潜力的关键指标。评价烃源岩的参数通常有三个，即有机质丰度、有机质类型和有机质成熟度。其中，有机质类型和有机质成熟度均首先对野外露头或已钻井岩心的烃源岩进行样品采集，再对烃源岩样品开展针对性的实验室分析测试，最终得到对应的有机质类型及有机质成熟度参数。对于烃源岩有机质丰度的评价，按研究范围可以分为两大类：对于单点的研究，如野外露头或某个已钻井烃源岩有机质丰度的评价，需要先开展烃源岩的样品采集，再对烃源岩样品进行实验室分析测试，得到总有机碳含量、烃源岩热解参数中的可溶烃、热解烃、氯仿沥青和总烃含量等表征有机质丰度的参数；对于某个盆地的区域范围的有机质丰度评价，一般以有效烃源岩分布的预测为主。具体的，盆地有效烃源岩分布的预测通常包含两方面内容，第一，烃源岩总有机碳含量的计算，第二，有效烃源岩厚度的计算，其中有效烃源岩是指总有机碳含量大于0.5％的烃源岩。常规的盆地有效烃源岩的预测研究，一方面是利用实验室样品分析测试技术，得到烃源岩的有机质总有机碳含量；另一方面，通过对野外露头烃源岩厚度的观测，估算地下岩层中有效烃源岩的厚度。但是，野外露头通常在盆地边缘，并且在勘探初期部署论证风险探井时，钻井岩心资料相对较少，这种情况下，没有办法实现对有效烃源岩分布的预测。
具体实施方式
[0053] 为了使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请作进一步地详细描述，所描述的实施例不应视为对本申请的限制，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本申请保护的范围。
[0054] 在以下的描述中，涉及到“一些实施例”，其描述了所有可能实施例的子集，但是可以理解，“一些实施例”可以是所有可能实施例的相同子集或不同子集，并且可以在不冲突的情况下相互结合。
[0055] 如果申请文件中出现“第一\第二\第三”的类似描述则增加以下的说明，在以下的描述中，所涉及的术语“第一\第二\第三”仅仅是是区别类似的对象，不代表针对对象的特定排序，可以理解地，“第一\第二\第三”在允许的情况下可以互换特定的顺序或先后次序，以使这里描述的本申请实施例能够以除了在这里图示或描述的以外的顺序实施。
[0056] 除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述本申请实施例的目的，不是旨在限制本申请。
[0057] 在介绍本申请实施例之前，对相关技术进行简单介绍：
[0058] 目前，盆地有效烃源岩分布的预测方法专利较少、文献较多，具体的技术方法主要分为三类：
[0059] 第一类：利用区域内已钻井的烃源岩分析测试数据，来概括性的评价有效烃源岩的分布。这种方法采用地质统计和实验室分析测试相结合的手段，建立评价和预测优质烃源岩的方法体系，预测盆地或某个区域范围的有效烃源岩分布。这种方法存在两个问题：(1)需要足够多的已钻井以及对应的分析测试数据，来构建评价和预测优质烃源岩的方法体系，而勘探初期往往没有足够多的分析测试数据；(2)利用这种方法得到的有效烃源岩分布特征，在无井处没有足够的数据约束，其预测精度受钻井密度的影响较大。
[0060] 第二类：利用测井曲线拟合的方法来预测有效烃源岩的分布，该方法通过测井曲线、分析测试数据以及某研究区小范围的三维地震数据来预测有效烃源岩的分布特征，但是，该类方法在解决盆地这种大范围的有效烃源岩分布预测问题时，该方法不适用。
[0061] 第三类：通过地质手段确定烃源岩有利沉积相的分布特征，用以表征有效烃源岩的分布。该方法通过分析测试手段构建沉积速率等沉积相参数与总有机碳含量的关系，并通过沉积相的分布来确定有利烃源岩的分布特征；但是，这种方法没有考虑沉积相的变化，而在盆地内的不同区域，经常会出现同一地质时期发育不同沉积相，即等时异向的现象，因此，这种方法也无法有效解决盆地这种大范围有效烃源岩的分布预测问题。
[0062] 从相关技术的方法来看，在勘探初期基础资料较为缺乏的前提下，尚无有效解决盆地有效烃源岩分布预测问题的技术方法。
[0063] 基于相关技术中存在的问题，本申请实施例提供一种有效烃源岩分布的预测方法，所述方法应用于电子设备，例如计算机、移动终端等。本申请实施例提供的有效烃源岩分布的预测方法所实现的功能可以通过电子设备的处理器调用程序代码来实现，其中，程序代码可以保存在计算机存储介质中。
[0064] 实施例一
[0065] 本申请实施例提供一种有效烃源岩分布的预测方法，图1为本申请实施例提供的一种有效烃源岩分布的预测方法的实现流程示意图，如图1所示，包括：
[0066] 步骤S1，获取目标区域中已钻井的总有机碳含量数据、钻井分层数据、测井曲线以及所述目标区域的二维地震测线数据。
[0067] 本申请实施例中，所述测井曲线可以包括：声波测井曲线、密度测井曲线、自然伽马曲线、孔隙度曲线、电阻率曲线等。
[0068] 本申请实施例中，钻井分层数据可以表征所述目标区域的地层的分层。
[0069] 本申请实施例中，可以通过测试设备、存储设备等的输入来目标区域中已钻井的总有机碳含量数据、钻井分层数据、测井曲线以及所述目标区域的二维地震测线数据。所述测试设备可以是总有机碳含量数据的分析设备、测井设备及地震数据测试设备等。
[0070] 在一些实施例中，电子设备可以和数据库通信连接，从数据库中获取目标区域中已钻井的总有机碳含量数据、钻井分层数据、测井曲线以及所述目标区域的二维地震测线数据。
[0071] 本申请实施例中，所述目标区域可以是盆地。所述盆地可以为低勘探程度区，其钻井基础资料较少。
[0072] 本申请实施例中，各个测井曲线的振幅数据在同一值域范围内。
[0073] 本申请实施例中，二维地震测线数据的振幅范围相同，且通过了闭合差处理。
[0074] 在一些实施例中，在步骤S1之前，所述方法还包括：
[0075] 步骤S11，获取所述目标区域中已钻井的初始测井曲线。
[0076] 步骤S12，对所述初始测井曲线进行标准化处理，得到各个测井曲线，其中，各个测井曲线的振幅数据在同一值域范围内。
[0077] 本申请实施例中，由于初始测井曲线通常是在不同测量时期，由不同测井仪器测量得到，因此需要消除不同测量时期、不同测井仪器对测井曲线的影响，可以对各个初始测井曲线进行标准化处理，使这些初始测井曲线的振幅数据处于同一值域范围内。
[0078] 在一些实施例中，在步骤S1之前，所述方法还包括：
[0079] 步骤S13，获取所述目标区域中不同区域的初始二维地震测线数据。
[0080] 本申请实施例中，所述不同区域可以包括：不同级构造单元或构造单元内的具体研究区。
[0081] 步骤S14，基于所述目标区域中至少两个标志性地层界面确定各个初始二维地震测线数据上所述两个标志性地层界面对应的时间深度值。
[0082] 本申请实施例中，由于这些不同级构造单元或构造单元内具体研究区的地表情况不同，造成不同级构造单元或构造单元内具体研究区的二维地震测线数据存在一定的时间差，即二维地震测线数据的闭合差。在利用这些二维地震测线数据开展盆地的烃源岩评价研究时，首先要消除这些闭合差，以提升烃源岩评价的精度。因此，需要对初始二维地震测线数据进行消除闭合差处理。
[0083] 本申请实施例中，可以通过目标区域中的地质认识来确定出盆地内至少两个标志性地层界面，示例性地，所述至少两个标志性地层界面可以包括：寒武系筇竹寺组顶界面以及志留系底界面作为标志层。
[0084] 本申请实施例中，至少两个标志性地层界面在地震剖面上较容易被识别，因此，可以在二维地震测线数据上读取标志性地层界面对应的时间深度值。
[0085] 步骤S15，基于各个时间深度值确定各个初始二维地震测线数据对应的闭合差。
[0086] 本申请实施例中，在确定了各个时间深度值后可求差，从而确定各个初始二维地震测线数据对应的闭合差。
[0087] 步骤S16，基于所述各个初始二维地震测线数据对应的闭合差对各个初始地震测线数据进行移动，以消除各个初始二维地震测线数据的闭合差，得到中间二维地震测线数据。
[0088] 本申请实施例中，可以基于闭合差来对各个初始地震数据进行整体的时间移动，从而实现不同区域的二维地震测线数据的闭合差消除。
[0089] 步骤S17，对所述中间二维地震测线数据进行振幅一致性处理，得到所述二维地震测线数据。
[0090] 本申请实施例中，在对盆地内不同区域的二维地震测线数据进行地震资料采集和处理时，由于采集时间不同、采集的仪器设备不同、使用的地震资料处理系统不同、处理的技术方法不同，都会造成不同级构造单元或构造单元内具体研究区的二维地震测线数据的振幅值值域不同，这会极大影响利用二维地震测线数据进行有效烃源岩评价的预测精度，因此需要对中间二维地震测线数据进行振幅一致性处理，得到所述二维地震测线数据，以使二维地震测线数据的振幅值在标准值域范围内。
[0091] 本申请实施例中，步骤S17可以通过以下步骤实现：
[0092] 步骤S171，统计所有中间二维地震测线数据的振幅值的值域范围。
[0093] 步骤S172，基于所述值域范围确定振幅值的标准值域范围。
[0094] 本申请实施例中，标准值域范围可以进行设定，可以是值域范围中的任意一个值域范围确定为标准值域范围，标准值域范围也可以是值域范围的中间值，也可以将最大的值域范围确定为标准值域范围。
[0095] 步骤S173，基于所述标准值域范围确定各个中间二维地震测线数据的幅度的校正系数。
[0096] 本申请实施例中，所述校正系数可以认为是放大或缩小的比例。
[0097] 步骤S174，基于所述校正系数对对应的二维地震测线数据的振幅值进行校正以得到所述二维地震测线数据，其中，各个二维地震测线数据的振幅值的值域在标准值域范围内。
[0098] 本申请实施例中，可以将各个二维地震测线数据的振幅值乘以校正系数，从而实现对二维地震测线数据的振幅值进行校正。
[0099] 本申请实施例中，在校正后，实现对不同区域内二维地震测线数据的振幅一致性处理。
[0100] 步骤S2，对所述总有机碳含量数据和所述测井曲线进行交会分析确定目标测井曲线，并确定所述目标测井曲线对应的关键参数与总有机碳含量之间的关系式。
[0101] 示例性地，可以开展烃源岩总有机碳含量与声波、密度、自然伽玛等测井曲线的交汇分析，制作交会分析图，从而可以确定与烃源岩总有机碳含量相关性最高的测井曲线或几条测井曲线的运算组合确定为目标测井曲线，则目标测井曲线对应的参数则为关键参数，示例性地，如果声波与烃源岩总有机碳含量之间的相关性最大，则关键参数则为声波。
[0102] 本申请实施例中，通过交会分析也可以确定烃源岩识别的关键参数与总有机碳含量的关系式。
[0103] 步骤S3，基于所述钻井分层数据和所述目标区域的地层剖面上的标志性地层界面特征对所述目标测井曲线和所述二维地震测线数据进行井震标定以确定时深关系。
[0104] 本申请实施例中，测井曲线都为深度域数据，而地震数据既有时间域也有深度域，二维地震测线数据通常是时间域数据，在进行地震‑关键参数反演之前，要充分利用盆地内少数已钻井资料进行井震标定，以建立准确的时间‑深度关系(即时深关系)。
[0105] 本申请实施例中，可以基于地质认识确定出若干个标志性地层界面，然后，结合钻井分层数据以及地震剖面上的标志性地层界面特征进行已钻井井震标定，从而确定出准确的时深关系。
[0106] 步骤S4，基于所述二维地震测线数据和所述时深关系进行层位构造解释得到构造解释层位，并基于所述构造解释层位确定所述目标区域的烃源岩地层框架模型。
[0107] 本申请实施例中，可以在时深关系的基础上基于二维地震测线数据开展烃源岩地层的层位构造解释，从而确定构造解释层位，利用二维地震测线数据的构造解释层位，构建出盆地的烃源岩地层框架模型。本申请实施例中，烃源岩地层框架模型可以标准目标区域的地层结构，二维地震测线数据可以对应于各个烃源岩地层框架模型中的位置。
[0108] 步骤S5，基于所述二维地震测线数据、烃源岩地层框架模型对所述关键参数进行反演计算，得到所述关键参数在所述烃源岩地层框架模型中的分布特征。
[0109] 本申请实施例中，由于二维地震测线数据和关键参数存在相关性，可以基于二维地震测线数据、烃源岩地层框架模型对所述关键参数进行反演计算，得到所述关键参数在所述烃源岩地层框架模型中的分布特征。
[0110] 步骤S6，基于所述关系式、所述分布特征确定所述目标区域中有效烃源岩分布特征。
[0111] 本申请实施例中，所述有效烃源岩分布特征包括：分布范围和厚度。
[0112] 本申请实施例中，步骤S6，可以通过以下步骤实现：
[0113] 步骤S61，确定有效烃源岩的总有机碳含量的门槛值。
[0114] 本申请实施例中，有效烃源岩是指总有机碳含量大于0.5％的烃源岩，则门槛值为总有机碳含量为0.5％。
[0115] 步骤S62，基于所述关系式和所述门槛值确定关键参数的门槛值。
[0116] 本申请实施例中，可以通过门槛值和关系式计算得到关键参数的门槛值。即求得烃源岩总有机碳含量为0.5％时关键参数的门槛值。
[0117] 步骤S63，基于所述关键参数的门槛值和所述分布特征，确定目标关键参数分布特征。
[0118] 本申请实施例中，可以基于关键参数的门槛值来剔除分布特征中的数据，从而确定目标关键参数分布特征。
[0119] 可以将大于总有机碳含量为0.5％对应的关键参数确定为目标关键参数分布特征。
[0120] 步骤S64，基于所述目标关键参数分布特征和所述关系式确定在所述目标区域中总有机碳含量的分布范围。
[0121] 本申请实施例中，可以将目标关键参数分布特征输入至关系式中，从而确定目标区域中总有机碳含量的分布范围。
[0122] 步骤S65，基于所述总有机碳含量分布范围确定有效烃源岩分布特征。
[0123] 可以通过烃源岩总有机碳含量的分布范围计出有效烃源岩的厚度，从而实现对盆地有效烃源岩分布特征的预测。
[0124] 本申请实施例提供的有效烃源岩分布的预测方法，通过获取目标区域中已钻井的总有机碳含量数据、钻井分层数据、测井曲线以及所述目标区域的二维地震测线数据；对所述总有机碳含量数据和所述测井曲线进行交会分析确定目标测井曲线，并确定所述目标测井曲线对应的关键参数与总有机碳含量之间的关系式；基于所述钻井分层数据和所述目标区域的地层剖面上的标志性地层界面特征对所述目标测井曲线和所述二维地震测线数据进行井震标定以确定时深关系；基于所述二维地震测线数据和所述时深关系进行层位构造解释得到构造解释层位，并基于所述构造解释层位确定所述目标区域的烃源岩地层框架模型；基于所述二维地震测线数据、烃源岩地层框架模型对所述关键参数进行反演计算，得到所述关键参数在所述烃源岩地层框架模型中的分布特征；基于所述关系式、所述分布特征确定所述目标区域中有效烃源岩分布特征，能够实现在基础资料(如岩心资料)较为缺乏的情况下，对有效烃源岩分布的预测。
[0125] 实施例二
[0126] 基于前述的各个实施例，本申请实施例再提供一种有效烃源岩分布的预测方法，结合钻井、测井及地震数据，提出了一种盆地有效烃源岩分布的预测方法。该方法主要如图2所示的步骤，包括烃源岩数据分析、二维地震测线数据校正、地震关键参数反演、烃源岩分布特征预测等四个步骤：
[0127] 步骤S101，烃源岩数据分析。
[0128] 本申请实施例中，烃源岩数据分析主要包括烃源岩数据收集整理和烃源岩识别关键参数确定两个步骤。
[0129] 步骤S1011，烃源岩数据收集整理。
[0130] 在低勘探程度区，往往钻井基础资料较少，因此，需收集盆地内以及周边地区已钻井的烃源岩分析测试数据、钻井及测井数据，包括烃源岩的总有机碳含量分析测试数据(同上述实施例中的总有机碳含量数据)、钻井分层数据以及声波、密度、自然伽玛等测井曲线数据。考虑到这些测井曲线数据是在不同测量时期，由不同测井仪器测量得到，因此，需要对收集到的声波、密度、自然伽玛等测井曲线进行标准化处理，消除不同测量时期、不同测井仪器对测井曲线数据的影响，使这些测井曲线数据处于同一值域范围内。
[0131] 步骤S1012，烃源岩识别关键参数确定。
[0132] 以已钻井的烃源岩分析测试数据以及标准化后的声波、密度、自然伽玛等测井曲线为基础，开展烃源岩总有机碳含量与声波、密度、自然伽玛等测井曲线的交汇分析，确定与烃源岩总有机碳含量相关性最高的测井曲线或几条测井曲线的运算组合，作为烃源岩识别的关键参数；同时，确定出烃源岩识别关键参数与总有机碳含量的关系式以及有效烃源岩识别关键参数的门槛值，即烃源岩总有机碳含量大于0.5％时关键参数的门槛值。
[0133] 步骤S102，二维地震测线数据校正。
[0134] 本申请实施例中，二维地震测线数据校正主要包括二维地震测线闭合差消除和振幅一致性处理两个步骤，其中，
[0135] 步骤S1021，二维地震测线数据闭合差消除步骤。
[0136] 在勘探初期，在盆地内进行的二维地震资料采集与处理通常是针对盆地内具体的不同级构造单元或构造单元内的具体研究区。由于这些不同级构造单元或构造单元内具体研究区的地表情况不同，造成不同级构造单元或构造单元内具体研究区的二维地震测线数据存在一定的时间差，即二维地震测线数据的闭合差。在利用这些二维地震测线数据开展盆地的烃源岩评价研究时，首先要消除这些闭合差。具体的，首先根据地质认识确定出盆地内的至少两个标志性地层界面，这种标志性地层界面通常在地震剖面上易于识别；然后，分别在不同级构造单元或构造单元内具体研究区的二维地震测线数据上读取标志性地层界面对应的时间深度值，并计算出闭合差的具体数值；最后，通过对二维地震测线数据进行整体的时间移动，实现对不同级构造单元或构造单元内具体研究区的二维地震测线数据的闭合差消除。
[0137] 步骤S1022，振幅一致性处理步骤。
[0138] 在对盆地内不同级构造单元或构造单元内具体研究区的二维地震测线数据进行地震资料采集和处理时，由于采集时间不同、采集的仪器设备不同、使用的地震资料处理系统不同、处理的技术方法不同，都会造成不同级构造单元或构造单元内具体研究区的二维地震测线数据的振幅值值域不同，这会极大影响利用二维地震测线数据进行有效烃源岩评价的预测精度。因此，需对不同级构造单元或构造单元内具体研究区的二维地震测线数据开展振幅一致性处理。具体的，首先，分别统计研究中使用的所有二维地震测线数据的振幅值值域范围，然后，根据所有二维地震测线数据的振幅值值域范围，确定出一个振幅值的标准值域范围；最后，通过振幅系数校正处理的方法，分别将不同二维地震测线数据的振幅值值域处理到同一个标准值域范围内，从而实现对不同级构造单元或构造单元内具体研究区二维地震测线数据的振幅一致性处理。
[0139] 步骤S103，地震关键参数反演。
[0140] 本申请实施例中，地震关键参数反演包括：已钻井井震标定、盆地烃源岩地层框架模型建立和二维地震烃源岩识别关键参数反演三个步骤，其中，各步骤介绍如下：
[0141] 步骤S1031，已钻井井震标定。
[0142] 测井曲线均为深度域数据，而地震数据既有时间域也有深度域。二维地震测线数据通常是时间域数据，那么，在进行地震关键参数反演之前，要充分利用盆地内少数已钻井资料进行井震标定，以建立准确的时间‑深度关系。具体的，首先根据地质认识确定出若干个标志性地层界面，然后，结合钻井分层数据以及地震剖面上的标志性地层界面特征进行已钻井井震标定，确定出准确的时深关系。
[0143] 步骤S1032，盆地烃源岩地层框架模型建立。
[0144] 首先，在建立的时深关系基础上，基于二维地震测线数据开展烃源岩地层的层位构造解释，然后，利用二维地震测线数据的构造解释层位，构建出盆地的烃源岩地层框架模型。
[0145] 步骤S1033，二维地震烃源岩识别关键参数反演。
[0146] 基于盆地内二维地震测线数据，结合得到的烃源岩识别的关键参数以及得到的盆地内的烃源岩地层框架模型，开展二维地震烃源岩识别关键参数反演，得到烃源岩识别关键参数在盆地内的分布特征。
[0147] 步骤S104，烃源岩分布特征预测。
[0148] 本申请实施例中，利用得到的烃源岩识别关键参数与总有机碳含量的关系式、有效烃源岩识别关键参数门槛值以及得到的烃源岩识别关键参数在盆地内的分布特征，可以换算出盆地烃源岩总有机碳含量的分布范围，并统计出有效烃源岩的厚度，从而实现对盆地有效烃源岩分布特征的预测。
[0149] 本申请实施例提供的方法可以在勘探初期基础资料较为匮乏的情况下，实现对盆地有效烃源岩分布特征的预测，将有利于低勘探程度区的油气勘探突破，并提高油气勘探的经济效益。
[0150] 基于前述的各个实施例，本申请实施例提供一种有效烃源岩分布的预测方法的应用示例，某盆地钻至寒武系筇竹寺组的钻井相对较少，在未使用该技术前，某盆地寒武系筇竹寺组的有效烃源岩分布特征不明确。具体应用中实施了如下几个步骤：
[0151] (1)收集整理盆地内已钻井的烃源岩总有机碳含量等分析测试数据、钻井分层数据以及测井曲线数据，并对已钻井的测井曲线进行标准化处理；然后，开展烃源岩总有机碳含量与声波、密度、自然伽玛等测井曲线的交汇分析，最终确定出纵波阻抗为烃源岩识别的关键参数，并得到总有机碳含量与纵波阻抗的关系式，以及有效烃源岩的纵波阻抗门槛值。
[0152] (2)选取寒武系筇竹寺组顶界面以及志留系底界面作为标志层，对盆地的二维地震测线数据进行闭合差校正，图3为本申请实施例提供的一种校准前后对比图，如图3所示，图3中的上图为闭合差校正前的地震剖面示意图，图3中的下图为闭合差校正后的地震剖面示意图。然后，统计出所有二维地震测线数据的振幅值值域范围；之后，利用振幅系数校正处理的方法，将所有二维地震测线数据的振幅值处理到同一振幅值值域范围内。
[0153] (3)开展已钻井井震标定，并在二维地震测线上解释出寒武系筇竹寺组的顶底界面，结合已钻井数据开展二维地震纵波阻抗反演，得到纵波阻抗在盆地内的分布特征。
[0154] (4)利用总有机碳含量与纵波阻抗的关系式、有效烃源岩的纵波阻抗门槛值以及纵波阻抗在盆地内的分布特征，计算出盆地烃源岩总有机碳含量的分布范围，并统计出有效烃源岩的厚度，从而实现对某盆地寒武系筇竹寺组有效烃源岩分布特征的预测。
[0155] 本申请实施例提供的有效烃源岩分布的预测方法，可以在勘探初期基础资料较为匮乏的情况下，实现对盆地有效烃源岩分布的预测。将有利于低勘探程度区的油气勘探突破，并提高油气勘探的经济效益。该方法包括：对盆地内烃源岩数据进行整理分析，并对已钻井的测井曲线进行标准化，结合烃源岩样品的分析测试数据确定烃源岩识别的关键参数、关键参数与总有机碳含量的关系式以及有效烃源岩的关键参数门槛值；利用盆地内的二维地震资料开展二维地震测线数据校正；之后，开展盆地范围的二维地震烃源岩识别关键参数反演，得到烃源岩识别关键参数在盆地内的分布特征；最终，基于烃源岩识别关键参数与总有机碳含量的关系式、有效烃源岩的关键参数门槛值以及关键参数在盆地内的分布特征，统计出盆地有效烃源岩的分布特征和厚度，从而实现在低勘探程度区对盆地有效烃源岩分布的准确预测。
[0156] 实施例三
[0157] 基于前述的实施例，本申请实施例提供一种有效烃源岩分布的预测装置，该装置包括的各模块、以及各模块包括的各单元，可以通过计算机设备中的处理器来实现；当然也可通过具体的逻辑电路实现；在实施的过程中，处理器可以为中央处理器(CPU，Central Processing Unit)、微处理器(MPU，Microprocessor Unit)、数字信号处理器(DSP，Digital Signal Processing)或现场可编程门阵列(FPGA，Field Programmable Gate Array)等。
[0158] 本申请实施例提供一种有效烃源岩分布的预测装置，有效烃源岩分布的预测装置包括：
[0159] 获取模块，用于获取目标区域中已钻井的总有机碳含量数据、钻井分层数据、测井曲线以及所述目标区域的二维地震测线数据；
[0160] 第一确定模块，用于对所述总有机碳含量数据和所述测井曲线进行交会分析确定目标测井曲线，并确定所述目标测井曲线对应的关键参数与总有机碳含量之间的关系式；
[0161] 第二确定模块，用于基于所述钻井分层数据和所述目标区域的地层剖面上的标志性地层界面特征对所述目标测井曲线和所述二维地震测线数据进行井震标定以确定时深关系；
[0162] 第三确定模块，用于基于所述二维地震测线数据和所述时深关系进行层位构造解释得到构造解释层位，并基于所述构造解释层位确定所述目标区域的烃源岩地层框架模型；
[0163] 反演模块，用于基于所述二维地震测线数据、烃源岩地层框架模型对所述关键参数进行反演计算，得到所述关键参数在所述烃源岩地层框架模型中的分布特征；
[0164] 第四确定模块，用于基于所述关系式、所述分布特征确定所述目标区域中有效烃源岩分布特征。
[0165] 在一些实施例中，所述有效烃源岩分布的预测装置还用于：
[0166] 获取所述目标区域中已钻井的初始测井曲线；
[0167] 对所述初始测井曲线进行标准化处理，得到各个测井曲线，其中，各个测井曲线的振幅数据在同一值域范围内。
[0168] 在一些实施例中，所述有效烃源岩分布的预测装置还用于：
[0169] 获取所述目标区域中不同区域的初始二维地震测线数据；
[0170] 基于所述目标区域中至少两个标志性地层界面确定各个初始二维地震测线数据上所述两个标志性地层界面对应的时间深度值；
[0171] 基于各个时间深度值确定各个初始二维地震测线数据对应的闭合差；
[0172] 基于所述各个初始二维地震测线数据对应的闭合差对各个初始地震测线数据进行移动，以消除各个初始二维地震测线数据的闭合差，得到中间二维地震测线数据；
[0173] 对所述中间二维地震测线数据进行振幅一致性处理，得到所述二维地震测线数据。
[0174] 在一些实施例中，所述对所述中间二维地震测线数据进行振幅一致性处理，得到所述二维地震测线数据，包括：
[0175] 统计所有中间二维地震测线数据的振幅值的值域范围；
[0176] 基于所述值域范围确定振幅值的标准值域范围；
[0177] 基于所述标准值域范围确定各个中间二维地震测线数据的幅度的校正系数；
[0178] 基于所述校正系数对对应的二维地震测线数据的振幅值进行校正以得到所述二维地震测线数据，其中，各个二维地震测线数据的振幅值的值域在标准值域范围内。
[0179] 在一些实施例中，所述基于所述总有机碳含量数据和所述测井曲线进行交会分析确定目标测井曲线，包括：
[0180] 将所述总有机碳含量数据和所述测井曲线进行交会分析确定总有机碳含量与各个测井曲线之间的相关性；
[0181] 基于所述相关性确定目标测井曲线。
[0182] 在一些实施例中，所述基于所述关系式、所述分布特征确定目标区域中有效烃源岩分布特征，包括：
[0183] 确定有效烃源岩的总有机碳含量的门槛值；
[0184] 基于所述关系式和所述门槛值确定关键参数的门槛值；
[0185] 基于所述关键参数的门槛值和所述分布特征，确定目标关键参数分布特征；
[0186] 基于所述目标关键参数分布特征和所述关系式确定在所述目标区域中总有机碳含量的分布范围；
[0187] 基于所述总有机碳含量分布范围确定有效烃源岩分布特征。
[0188] 需要说明的是，本申请实施例中，如果以软件功能模块的形式实现上述的有效烃源岩分布的预测方法，并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。这样，本申请实施例不限制于任何特定的硬件和软件结合。
[0189] 相应地，本申请实施例提供一种存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现上述实施例中提供的有效烃源岩分布的预测方法中的步骤。
[0190] 实施例四
[0191] 本申请实施例提供一种电子设备；图4为本申请实施例提供的电子设备的组成结构示意图，如图4所示，所述电子设备700包括：一个处理器701、至少一个通信总线702、用户接口703、至少一个外部通信接口704、存储器705。其中，通信总线702配置为实现这些组件之间的连接通信。其中，用户接口703可以包括显示屏，外部通信接口704可以包括标准的有线接口和无线接口。所述处理器701配置为执行存储器中存储的有效烃源岩分布的预测方法的程序，以实现以上述实施例提供的有效烃源岩分布的预测方法中的步骤。
[0192] 以上电子设备和存储介质实施例的描述，与上述方法实施例的描述是类似的，具有同方法实施例相似的有益效果。对于本申请计算机设备和存储介质实施例中未披露的技术细节，请参照本申请方法实施例的描述而理解。
[0193] 应理解，说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本申请的至少一个实施例中。因此，在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外，这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解，在本申请的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。上述本申请实施例序号仅仅为了描述，不代表实施例的优劣。
[0194] 需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
[0195] 在本申请所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。
[0196] 上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
[0197] 另外，在本申请各实施例中的各功能单元可以全部集成在一个处理单元中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
[0198] 本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器(ROM，Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0199] 或者，本申请上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台控制器执行本申请各个实施例所述方法的全部或部分。而前述的存储介质包括：移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。
[0200] 以上所述，仅为本申请的实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。
