基于多元线性回归分析的页岩气甜点预测方法与流程

本发明涉及页岩气领域，具体而言，涉及一种基于多元线性回归分析的页岩气甜点预测方法。

背景技术：

页岩气资源丰富，有望缓解面临的能源危机。然而，由于页岩地层复杂，其勘探开发难度大、且页岩气钻井成本远高于常规石油钻井。这就需要在对页岩气储层进行开发时，准确预测并识别未来勘探的“甜点”潜力区。

据jarvie等(2007)给出北美页岩甜点的参考指标主要为静态参数，包含三大类型，页岩分布(埋藏深度、页岩厚度与横向分布、断裂分布)；页岩组成(有机质类型和丰度，粘土矿物、脆性矿物)；页岩属性(热成熟度、孔隙度、渗透率)。但我国四川盆地及周缘海相页岩层系热演化程度高、埋藏时间长、经历构造期次多导致其地质条件复杂。其页岩气含气量影响因素不仅受前期静态沉积指标影响，更受后期构造运动导致的保存参数变化影响。可知，现有北美页岩评价方法难以满足对海相页岩甜点的预测。

此外，目前存在的这种甜点地质参数评价与优选甜点区的方法均存在以定性描述为主、简单主观赋值打分、多参数综合评价仅是通过简单叠加、保存条件指标评价单一难以量化等问题。这就导致预测过程中，难以对地质特征参数进行量化，叠加法预测的结果可比较性差；同时，没有考虑各个地质特征参数对含气量影响的权重，叠加法得到的预测平面图也无法判断研究区内的甜点区域或趋势的变化。

技术实现要素：

本发明的目的在于提供一种基于多元线性回归分析的页岩气甜点预测方法，该方法通过多元线性回归分析方法定量的对各甜点地质参数进行评价，得到最终综合评价甜点区域，甜点位置更细化、准确，且甜点顺序面积显示更直观。

为了实现本发明的上述目的，特采用以下技术方案：

一种基于多元线性回归分析的页岩气甜点预测方法，包括：获取储层甜点、保存甜点和压力系数甜点中的各项地质参数，评价与页岩气含气量相关的地质特征参数，并量化成平面分布图，对每个所述地质特征参数与所述含气量进行线性回归拟合，得到每个所述地质参数与所述含气量的回归方程，每个所述地质参数的权重为其对应的所述回归方程的相关系数，根据各地质参数的实测值及其对应的权重，量化叠加后得到所述待勘探区页岩储层的页岩气甜点预测模型。

进一步地，在本发明较佳的实施方式中，在得到回归方程后，还包括对回归方程以及回归系数进行显著性检验。

进一步地，在本发明较佳的实施方式中，上述储层甜点中的地质参数包括：有机碳含量、热演化程度、孔隙度、渗透率及脆性矿物含量；保存甜点中的地质参数包括：埋深和露头距离；压力系数甜点中的地质参数包括压力系数。

进一步地，在本发明较佳的实施方式中，上述页岩气甜点的预测模型如下：

gas＝a·dep+b·toc+c·ro+d·φ+e·k+f·bm+g·p+h·dist

式中，dep为待勘探区页岩储层的埋深；toc为待勘探区页岩储层的有机碳含量；ro为待勘探区页岩储层的热演化程度；φ为待勘探区页岩储层的孔隙度；k为待勘探区页岩储层的渗透率；bm为待勘探区页岩储层中脆性矿物含量；p为待勘探区页岩储层的压力系数；dist为待勘探区页岩储层中页岩露头距离；a、b、c、d、e、f、g以及h均为各相应所述地质参数与所述含气量经线性回归拟合后所得的回归方程中的相关系数。

进一步地，在本发明较佳的实施方式中，上述页岩气甜点的含气量≥2m³/t。

进一步地，在本发明较佳的实施方式中，上述页岩气甜点的有机碳含量≥2.5％。

进一步地，在本发明较佳的实施方式中，上述页岩气甜点的热演化程度为1.1～3.0％。

进一步地，在本发明较佳的实施方式中，上述页岩气甜点的脆性矿物含量≥40％。

有现有技术相比，本发明的有益效果为：

本申请提供的基于多元线性回归分析的页岩气甜点预测方法，通过系统分析储层甜点、保存甜点以及压力系数甜点中各项地质参数对页岩含气量的影响，建立每个地质参数与含气量之间的关系模型，结合多元回归分析，明确不同构造区甜点的地质参数权重并量化预测待勘探区页岩储层的页岩气甜点区的分布。通过这种预测方法所得到的甜点区分布位置更细化、准确，且甜点顺序面积显示更直观，有助于避免因为简单图层叠加造成的主观判断错误，为未来建立成熟探区寻找遗漏的“甜点”。且该方法有助于对低勘探程度下的构造复杂区页岩气甜点进行量化选区。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，以下将对实施例或现有技术描述中所需要使用的附图作简单地介绍。

图1为本发明实施例提供的研究区龙马溪组沉积相分布定量分布图。

图2为本发明实施例提供的研究区含气量与龙马溪组页岩厚度关系图。

图3为本发明实施例提供的研究区龙马溪组富有机质泥岩厚度分布图。

图4为本发明实施例提供的研究区含气量与龙马溪组页岩有机碳含量关系图。

图5为本发明实施例提供的研究区龙马溪组有机碳含量分布图。

图6为本发明实施例提供的研究区含气量与龙马溪组页岩热演化程度关系图。

图7为本发明实施例提供的研究区龙马溪组热演化程度图。

图8为本发明实施例提供的研究区龙马溪组孔隙度分布图。

图9为本发明实施例提供的研究区含气量与孔隙度关系图。

图10为本发明实施例提供的研究区含气量与比表面关系图。

图11为本发明实施例提供的研究区含气量与渗透率关系图。

图12为本发明实施例提供的研究区矿物成分三角关系图。

图13为本发明实施例提供的研究区生物成因硅质含量与有机碳含量关系图。

图14为本发明实施例提供的研究区龙马溪组脆性矿物分布图。

图15为本发明实施例提供的四川盆地盆缘残余向斜龙马溪组距离露头距离图。

图16为本发明实施例提供的研究区含气量与页岩露头出露距离关系图。

图17为本发明实施例提供的研究区龙马溪组埋深等值线图。

图18为本发明实施例提供的研究区龙马溪组地层构造斜率分布图。

图19为本发明实施例提供的研究区龙马溪组埋深等值线图。

图20为本发明实施例提供的各甜点影响因素拟合点图。

图21为本发明实施例提供的多元回归参数与检验结果图。

图22为本发明实施例提供的四川盆地及周缘海相页岩复杂区预测甜点分布图。

具体实施方式

下面将结合实施例对本发明的实施方案进行详细描述，但是本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限制本发明的范围。实施例中未注明具体条件者，按照常规条件或制造商建议的条件进行。所用试剂或仪器未注明生产厂商者，均为可以通过市售购买获得的常规产品。

本实施方式提供一种基于多元线性回归分析的页岩气甜点预测方法。

由于现有技术中对“甜点”的定义多为定型化描述、无法定量，且不同学者对“甜点”的定义不同。在本发明中，“甜点”指油气富集的、在当前经济技术条件下可有效开发的层段和区域，考虑广泛适用性，以含气量>2m³/t作为参考指标。即，在本专利中，页岩气甜点含气量的下限是2m³/t。

此外，选取四川盆地及周缘区域为研究区作为本发明实施方式的方法的应用区，对本发明的方法进行详细说明。换言之，本发明实施方式中的方法也可以应用到其它页岩气有利区的甜点预测。

这种基于多元线性回归分析的页岩气甜点预测方法，包括：

步骤1：分析与待勘探区页岩储层的地质环境相似的多个页岩气钻井的页岩气地质特征参数信息，以获取每个页岩气钻井的含气量、储层甜点、保存甜点以及压力系数甜点中的各项地质参数，并量化成平面分布图。

其中，储层甜点、保存甜点、压力系数甜点的参数为：

1.储层甜点包含优质相带、地球化学甜点、储集物性甜点三方面，其中，优质相带评判参数包括优质页岩厚度及岩相古地理展布；地球化学甜点评价参数的选取主要考虑有机质丰度、成熟度；储集物性甜点包括页岩孔隙度、渗透率、比表面及脆性矿物分布。

2.保存甜点包含构造样式(褶皱样式、断裂样式、褶皱-断裂组合样式、埋藏史类型)、封盖类型(封盖层组合类型、构造封盖类型、顶底板特征)分析。

3.压力系数甜点主要考虑压力系数分布。

较为具体地，影响页岩气甜点的主要地质因素包括：

(1)沉积相带与厚度分布：

针对四川盆地及周缘区域沉积特征、地层分布特征等方面进行研究，表明：海相优质烃源岩的发育受水体中生物生产率、沉积速率、海底深部流体活动、沉积环境和水体环境等多方面因素控制，沉积环境还可以综合反映原始生产力、沉积速率及有机质保存条件等因素。其中，深水陆棚沉积物主要形成缺氧、滞留环境沉积，有利页岩有机质保存，故对页岩气生成有利；向滨岸过度时有机质减少，页岩生成物质基础减弱。因此，在对研究区内沉积相分布及厚度分布等值线图进行数值化处理中，将沉积相按深水陆棚、浅水陆棚及滨岸体系分别赋值>2、>1和<0，如图1所示。

通过统计四川盆地及周缘构造复杂去页岩钻井成果，表明含气量与富有机质泥岩厚度有明显正相关关系(如图2所示)，四川盆地五峰-龙马溪组黑色页岩分布稳定，厚度较大，对页岩甜点区形成提供有利支撑(如图3所示)。

(2)地化甜点：

沉积岩中有机质的丰度和类型是生成油气的物质基础，但有机质只有一定量的有机物质并达到一定的热演化程度才能开始大量生烃。通过对热演化程度与含气量指数关系进行分析，认为最有利于页岩气形成与富集的热演化程度(ro)为1.1～3.0％。

请参照图4和图5，四川盆地内钻井数据表明，toc与含气量的呈正相相关关系。若以含气量2m³/t为甜点下限，toc甜点下限约为2.5％，在区内高值区主要分布于长宁、涪陵-武隆、巫溪一带。

请参照图6和图7，四川盆地盆内与盆外钻井ro与含气量呈先增后减趋势。盆内页岩埋藏深时间长，热演化程度更高，处于干气阶段大量有机质孔隙发育并逐渐开始坍塌，影响储集空间，故含气量有一定程度降低；盆外因抬升早，主体仍处于高-过成熟演化阶段，较盆内热演化程度低，故在热演化较低和适中的区域含气量更高，在高演化区间含气量同样逐步降低。

(3)储集性甜点：

请参照图8，对研究区志留系龙马溪组的孔隙度与含气量进行统计，范围为0.5～7.5％。请参照图9，相关性分析结果表明，含气量随孔隙度增大而线性增加。此外，孔隙度与含气量关系出现分区特征：盆外复杂区：两低(低孔隙度低含气量)；盆内高陡带：两高(高孔隙度高含气量)。

请参照图8，研究区五峰组-龙马溪组孔隙度特征在平面分布特征总体上在威远、宜宾-泸州为高值，盆内向盆外逐步减少。

请参照图10，对研究区内比表面积与含气量进行统计表明，比表面积和含气量呈正相关关系。

请参照图11和12，对志留系龙马溪组含气量和渗透率数据进行相关性分析，四川盆地及周缘龙马溪组页岩含气量与渗透率相关性不甚明显，在盆内外相关系数分别为负相关性和正相关性，说明其控制因素可能存在较大差异。如盆内保存条件好，裂缝相对不发育，渗透率越低含气性经改造后越高，而盆外总体含气量较低，相关性不明显，弱正相关性可能反映出储层孔隙间连通性较差，受裂缝沟通使含气量稍微增加的可能。

脆性矿物含量与含气量密切相关，直接影响页岩气后期压裂。脆性矿物含量高，可压性好。研究区内脆性矿物含量总体大于40％，主要以生物成因硅为主，且与含气量有一定正相关性。由此说明随着生物含量的增加，生气物质基础增加、保存条件变好，进而含气量增加。平面上，硅质矿物含量呈现由盆内西南向盆外东南增加的趋势(如图12、13和14所示)，总体事宜压裂。

(4)保存甜点：

总体上，盆内保存条件较好，仅高陡构造及川南长宁背斜主体附近保存条件中等,主要考虑埋深或者局部深大断裂影响；而盆外褶皱区则需要考虑褶皱宽度、倾角、埋深及断裂分布，产状高陡(或曲率较大)曲率较大部位通常构造变形较强，断裂和裂缝较发育，影响页岩含气性。

请参照图15，盆内构造演化相对简单，但埋深变化大，存在隐伏断裂，故断裂分布和埋深影响更甚；盆外复杂区向斜单元形态各异，构造演化复杂，多期断裂叠加，埋深变化大，故露头距离、构造曲率(倾角)、断裂分布与埋深均需要考虑。

研究钻井离露头出露距离与含气量的关系表明，钻井离露头距离越远含气量越高，呈明显正向相关关系(如图16和表1所示)。

表1.四川盆地页岩气钻井统计表

请参照图17、18和19，盆内川西南与川东高陡褶皱带地层倾角变化相对较大，其余区域产状总体平缓，地层埋藏深度除高陡向斜埋深较浅外其余总体埋深较大；盆外背斜分布区断裂发育，龙马溪组地层被断层、褶皱破坏、地层埋深较小、地层连续连片分布的面积较小，而向斜发育区，地层埋藏深度变化较大、连续连片分布地层面积整体呈现褶皱核部地层较缓、两翼产状陡峭特征。断层本次主要考虑一级、二级断裂，分布主要位于盆缘及盆外部位，距一级断裂10km以上保存条件较好；距二级走滑断裂3-5km以上保存条件较好；距三级断裂1～3km以上保存条件较好。

步骤2：对每个地质参数与含气量进行线性回归拟合，得到每个地质参数与含气量的回归方程，每个地质参数的权重为其对应的回归方程的相关系数，根据各地质参数的实测值及其对应的权重，量化叠加后得到待勘探区页岩储层的页岩气甜点预测模型。

作为优选，在得到回归方程后，还包括对回归方程以及回归系数进行显著性检验。

进一步地，上述页岩气甜点的预测模型如下：

gas＝a·dep+b·toc+c·ro+d·φ+e·k+f·bm+g·p+h·dist

式中，dep为待勘探区页岩储层的埋深；toc为待勘探区页岩储层的有机碳含量；ro为待勘探区页岩储层的热演化程度；φ为待勘探区页岩储层的孔隙度；k为待勘探区页岩储层的渗透率；bm为待勘探区页岩储层中脆性矿物含量；p为待勘探区页岩储层的压力系数；dist为待勘探区页岩储层中页岩露头距离；a、b、c、d、e、f、g以及h均为各相应所述地质参数与所述含气量经线性回归拟合后所得的回归方程中的相关系数。

以下结合实施例对本发明的特征和性能作进一步的详细描述：

实施例

本实施例提供一种基于多元线性回归分析的页岩气甜点预测方法，该方法根据页岩气储层的基本特征，结合四川盆地及周缘地质特点，以地质参数特征为指导，通过对比多元线性回归分析对页岩气“甜点”进行识别预测。

其包括以下步骤：

a.收集四川盆地及周缘页岩气钻井的页岩气地质特征参数信息，得到四川盆地及周缘页岩地质参数，如表2所示，再将各地质参数量化成平面分布图；

b.以含气量为因变量，其他地质参数为自变量，对每个地质参数与含气量进行线性回归拟合，得到每个地质参数与含气量之间的回归方程，每个地质参数的权重为其对应的回归方程的相关系数，如图20所示；

c.对步骤b中所得的回归方程以及回归系数进行显著性检验，结果如图21所示。

d.根据各地质参数的实测值及其对应的权重，量化叠加后得到待勘探区页岩储层的页岩气甜点预测模型：

gas＝0.0019*depth+0.466*toc+0.9594*ro+0.034*φ-0.007*k+0.03*bm-0.63*pressure-0.084*dist.

通过上式所示的甜点预测模型，结合量化参数平面图(选取实测值)，获得四川盆地及周缘海相页岩复杂区甜点分布图(如图22所示)。由图可见，在整个研究区内，靠近四川盆地内灰度较深的部分以及盆内和盆外分裂带灰度变化的部分为页岩气甜点预测位置。并且，能够有效预测到盆地内甜点的位置的变化趋势，有助于在低勘探区进行页岩气的开发。由此可知，通过该方法预测的结果具备较高的可比较性，可直观判断有利区块内甜点区域和非甜点区域。另一方面，通过多元线性回归分析各个影响参数的权重，通过量化叠加最终预测页岩气甜点区域以及优化排序，对未开采的有利区块内甜点区域具有明确指导意义。

表2.四川盆地及周缘钻井参数表

由此可见，本发明提供的这种基于多元线性回归分析的页岩气甜点预测方法，对甜点区分布位置的预测更加细化、准确，且甜点顺序面积显示更直观，有助于避免因为简单图层叠加造成的主观判断错误，为未来建立成熟探区寻找遗漏的“甜点”。且该方法有助于对低勘探程度下的构造复杂区页岩气甜点进行量化选区。

尽管已用具体实施例来说明和描述了本发明，然而应意识到，在不背离本发明的精神和范围的情况下可以作出许多其它的更改和修改。因此，这意味着在所附权利要求中包括属于本发明范围内的所有这些变化和修改。