CN102562048A - 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 - Google Patents

Abstract

本发明涉及一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，属油气勘探开发技术领域，其特征在于通过对影响孔隙度的地质作用过程分析、模拟和定量评价，建立起地质过程响应关系，确定储层原始孔隙保留和次生孔隙增减情况，进而定量评价处于高成岩阶段的低孔低渗储层目前的孔隙度，预测有效储层发育，为油藏评价提供依据；本发明预测性强，能够准确预测处于高成岩各阶段各类碎屑岩有效储层的分布，可对研究目标区段有效储层分布进行精细预测，预测单元平面可达到百米级；理论依据扎实，避免了传统方法中对采样点数据的过份依赖；定量性好，适用性广，资料要求相对较低；为不同勘探开发阶段具有不同资料基础的评价提供了方便。

Description

一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法

技术领域：

本发明涉及一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，属油气勘探开发技术领域。

背景技术：

处于高成岩阶段的低孔、低渗储层作为油气勘探开发的新领域受到了广泛关注。低孔、低渗储层在我国各大油区都有发现，如长庆、大庆、辽河、吉林、克拉玛依、吐哈、中原、青海、四川等油气田，都有大量的处于高成岩阶段的低孔、低渗油气藏，特别是近年来随着开发工艺的提高，低渗油气藏的产量不断提高，大量的低孔低渗储量得到有效开发，使得低渗储层的重要性日渐突出，而近年来新发现的油气储量中，有近一半来自低孔、渗的碎屑岩储层，使得低孔渗储层的地位更为重要。

有效储层预测是低孔渗储层油藏评价的难题和关键。对受成岩作用影响较大的低孔、低渗储层的预测，目前多采用分析化验资料进行统计预测、测井资料进行井间插值和地震资料反演等方法进行预测。

基于分析化验资料所提供的直接资料进行有利储层预测有两种方法，一是通过对储层直接的分析物性或微观孔隙特征的研究，统计面孔率或测量其储、渗参数，进而确定其有效性。这种方法对于有钻井，特别是有取心的井点可以完成；而对于没有钻井控制的区域，或没有直接取心不具有分析化验资料的区域则无法进行相应的评价。而且由于其成本较高，往往取心有限，因而难于全区控制和预测。另一种方法则是通过成岩相研究，通过对成岩作用现象分析，确定其所处的成岩阶段，明确其主要孔隙类型，按成岩阶段大致确定储层的有效性。该方法预测过于模式化，而且没有考虑不同岩石类型由于沉积过程差异所造就的原始物性的差异，虽处于同一成岩阶段但所处的成岩环境的差异所造成的成岩作用效果的差异，因而对储层的预测不够准确，难于识别有效储层发育带，并且不具有定量性。

测井预测是最常用的储层预测方法，但其不足也是显而易见的。测井预测主要基于岩电关系建立的测井解释模版，对已钻井未进行取样的层段进行预测，确定其储油物性。若解释模版可靠，该方法能够较好地给出测井段储层的物性参数。受测井资料的影响导致的解释模版的可靠性对储层物性解释的影响较大，常使得解释结果与实际情况有所偏差。而该方法本身也是一种基于井点数据的方法，对于井间参数分布无法给出，常采用统计分析方法，对井间进行插值给出井间储层参数，确定井间储层分布。但基于井点的井间插值本身多只考虑了测井点的数据分布所形成的地质场的特征，而这一特征可能受取样点的限制而不能很好地反映地下真正的地质参数场的分布，从而不能很好地体现不同的地质作用过程造成的储层特征的不均一性，使得预测结果的可靠性受到极大的影响。特别重要的是在勘探初期，井点数据非常有限，常不能进行较好的平面控制，因而难于对储层进行有效的预测。

地震资料具有较好的平面连续性，因而对储层平面上的特征可具有较连续的反映。利用地震资料预测储层参数主要是依据地震反射波在不同地层中传播特征的差异分析来实现的。在精细井震标定基础上，得到储层的地震响应特征，按地震响应特征预测储层发育及其属性特征，因而井震响应模型便成为预测的关键。然而由于地震属性受多种因素影响，同一种响应可能是不同属性作用的结果，难于从中筛选出特别好的属性关系，特别是对于深层高成岩阶段低渗储层，常由于埋深大而地震分辨率低导致标定困难或标定不准，使得难于形成较好的回归模型，预测结果往往很不理想，常常难于识别出储层，或难于区分出质量较好的储层。

基于这一背景，利用间接的储层测量数据已不能较好地反映出储层的发育，特别是有效储层的发育分布状况，影响了处于高成岩阶段低孔、低渗储层的油气勘探开发，必须寻找出一种不同于常规的评价方法，更好地预测储层分布，为油气勘探开发提供依据。

过程响应原理是基本的科学原理，有什么样的作用过程，就有什么样的响应特征。储层是地质历史中沉积、成岩和构造作用的共同结果，考虑到不同时期储层发生的各种地质作用，特别是不同时期地质作用发生时条件的差异，以及这些差异化地质条件作用的结果，可以为处于高成岩阶段的储层预测提供很好的手段。

发明内容：

为了克服现有技术的不足，本发明的目的在于提供一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，通过对影响孔隙度的地质作用过程的分析、模拟和定量评价，建立起地质过程响应关系，确定储层原始孔隙保留和次生孔隙增减情况，进而定量评价处于高成岩阶段的低孔低渗储层目前的孔隙度，预测有效储层发育，为油藏评价提供依据。

本发明是通过如下技术方案来实现上述目的的。

1、地质资料的收集。

收集的地质资料包括以下内容：目标地区的地层序列及各地质年代持续时间；目标地区钻井分层数据、目标层段沉积相研究成果图件、不同地层界面构造图；目标层段的热成熟度分析数据、胶结物含量、溶蚀孔隙度、残余孔隙度等数据。

2、原始资料数据库的建立。

原始资料数据库的建立包括以下几个方面：

(1)热成熟度指标TTI值与成熟度Ro值及成岩阶段之间关系的理论关系数据库；

(2)地质年代符号与地质年代持续时间数据库；

(3)沉积相类型、符号序列数据库；

(4)结构相与原始孔隙度关系数据库；

(5)用于约束模型的单井模拟数据库，包括单井分层数据、单井沉积相分析成果数据、单井储层结构参数、单井成熟度分析数据、单井成油成气时代分析数据、单井孔隙度分析数据、单井孔隙类型及孔隙构成数据；

(6)地质模型数据库，包括研究区钻井分层数据、不同地层界面构造图、目标层段沉积相平面图等。

3、以井点资料为基础，模拟约束井的沉积充填史、热源演化史，确定不同时期的地温梯度。

(1)以井点分层数据为基础，回朔井点的沉积过程，恢复沉积史。

地层埋藏深度恢复即用已知某地质时间的地层埋深减去在此地质时间后的地层厚度作为该地质时间的地层埋藏深度。同时对地层的压实作用造成的地层厚度的变化进行恢复，减除地层压实造成的埋深的影响。而构造恢复则主要考虑构造造成的地层的分析目标层段埋藏后几何形态，确定其所处的构造形态特征，明确其流体流动指向，推断其能否接受外部流体或与外界流体有交流，确定其与外界的连通状况，为分析胶结、溶解作用提供封闭环境判定依据。

(2)以井点分层热成熟度分析数据为基础，模拟井点热演化史，确定井点地温梯度。

即在埋藏史分析基础上，给定原始的地温场模拟目标层段的热演化，得到目前的热演化程度，与实测的热演化程度进行对比，按照差异对温度场进行调整，最终使得模拟井点热演化程度与井点测值达到一致，设定此时的地温场为此井点的真实地温场。

4、以井点资料为基础，模拟约束点的各沉积相的结构相类型，确定对不同沉积相原始孔隙度的影响参数，建立不同类型储层原始孔隙度的预测模型。

通过对储层的岩石类型、岩石结构、沉积相等方面研究，按照理论指导，确定其原始孔隙度。

(1)读入模拟井点的沉积相数参数、储层结构参数、原始孔隙度分析参数；

(2)由沉积相参数与储层结构参数判定其结构相；

(3)回归结构相、原始孔隙度关系，确定不同结构相的孔隙影响系数。

5、以井点资料为基础，模拟约束点的压实演化史，确定压实减孔参数，固化压实减孔模型。

(1)压实演化史模拟确定压实减孔系数。

压实作用导致储层颗粒排列紧密，从而减少储层孔隙。而压实作用又与储层本身岩石结构、颗粒成份、颗粒抗压能力、早期胶结、孔隙流体外排情况、埋藏过程等相关。根据压实减孔量与深度、结构相之间的关系，可以确定不同埋深下不同沉积砂体的减孔量。

6、恢复钙质胶结史、硅质胶结史、粘土矿物胶结史，确定钙质胶结指数、硅质胶结指数、粘土胶结指数，固化胶结减孔模型。

胶结作用是储层成岩作用中的重要作用，通过胶结作用，使原有的相对松散的沉积粘结一起而转化为岩石。在胶结作用过程中，在岩石系统内形成了新物质，而这些物质只能占据原有的孔隙空间，从而使孔隙减少。对于不同的岩石类型、不同的地质温度场和化学场，特别是孔隙介质的酸碱度，会大大影响储层的胶结作用。同时胶结物的来源也是影响胶结作用的重要因素。通过对不同环境下影响胶结物形成因素的分析，建立起不同类型胶结物(粘土、钙质、硅质、铁质)在不同的酸碱度、温度、压力、外部连通情况、岩石组构下的析出模型，通过计算得到胶结物的生成量，即胶结物减孔量。将获得的胶结物量与测点处实测数据进行对比，找出差异，调整理论模型中相应参数，再进行正演获取新胶结物量，再对比，直至两者一致，固定模型参数。

7、恢复长石等组构溶解史，确定各组构溶解系数，固化各组构溶解模型。

溶蚀孔隙是高成岩阶段储层的主要孔隙类型，溶蚀孔隙发育取决于储层中原有物质或后期生成物质的溶解作用。影响这些物质溶解的因素一是这些溶解物质的含量，二是溶解时的物化环境，三是溶解时的流体系统的封闭性。

通过建立的不同环境下的溶解模型，可对不同物质成份的溶解状况，主要是钙质、长石、岩屑等易溶物质的溶解量进行定量计算，确定新增的溶解孔隙量。

利用建立的溶解模型，设定原始模拟参数，求取测点减孔量，将之与测量值对比，找出差异，调整理论模型中相应参数，再进行正演获取新溶解量，再对比，直至两者一致，固定模型参数，对所有网格进行计算，获得真实的溶解物总量，获取溶解增孔量。

8、读入分层数据、沉积相数据及构造模型，建立地质模型。

利用成熟的建模软件(如PETREL)，建立研究区的地质模型。

①读入地层分层数据，建立不同地质单元的构造模型，体现不同地质单元的地层的顶底面的形态及其变化；

②读入模拟目标层段的沉积相研究成果，建立目标层段沉积相模型。

9、模拟目标区的沉积充填史。

由地质模型中的各层顶底面深度数据，逐点恢复目标地区的地层充填埋藏史。

10、模拟目标区的热演化史，确定目标地层不同地质时间时的温度、热成熟度、所处成岩阶段。

(1)由地质模型和井点的地温梯度计算各点各地质时期的地温梯度。

(2)依据各地质时期的地温梯度计算目标层各点的热成熟度指标TTI值。

(3)将TTI值转化Ro值。

(4)由Ro值确定各点所处的成岩作用阶段。

(5)由成岩阶段数据库和成岩阶段与孔隙类型关系判断目标层段在各地质时期的储集空间类型。

11、依据地质模型，模拟目标层段的原始孔隙度。

(1)读取模型中目标层段各点的沉积相类型。

(2)由沉积相数据计算各类储层的结构相。

(3)按结构相类型判断其原始孔隙度因子。

(4)由原始孔隙度因子与其区域储层孔隙度的值计算各网格点的原始孔隙度。

12、模拟目标区的压实减孔史，确定不同时期的累积压实减孔量。

(1)逐个读取目标层段在不同时期的地层埋藏深度值和地温梯度值。

(2)按压实减孔模型，由老到新逐期、逐点计算目标层段不同时代的地层压实减孔量。

(3)由老到新计算不同地质时代累计压实减孔量。

13、模拟目标区的钙质胶结史、硅质胶结史、粘土胶结史，确定不同时期的累积胶结减孔量。

(1)逐个读取目标层段在不同时期的地层埋藏深度值和地温梯度值。

(2)按胶结减孔模型，由老到新逐期、逐点计算目标层段不同时代的地层胶结减孔量。

(3)由老到新计算不同地质时代累计钙、硅和粘土胶结减孔量。

14、模拟目标区的矿物溶解史，确定不同时期的溶解增孔量。

(1)逐网格读取目标层段在不同时期的地层埋藏深度值和地温梯度值。

(2)按溶解增孔模型，由老到新逐期、逐点计算目标层段不同时代的地层溶解增孔量。

(3)由老到新计算不同地质时代累积溶解增孔量。

15、计算不同地质时期目标层段的总孔隙度。

高成岩阶段的总孔隙是原生孔隙经改造后的剩余孔隙，是原生孔隙中残留孔隙与新生的溶解孔隙的总和，依据计算的原生孔隙，加上新增的溶解孔隙，减去压实、胶结减少的孔隙便可得到高成岩储层目前的孔隙度即：

其中：

为某地质时刻储层总的孔隙度

为某储层原始孔隙度

为第i种地质作用增减的孔隙度

16、按总孔隙度对储层进行评价，确定有效储层的发育区间。

读入储层判定标准，按储层的孔隙度大小对各点的储层有效性进行判断，确定其有效性，进而由Petrel给出有效储层的发育范围。

本发明与现有的技术相比，具有如下有益效果：

1、预测性强，能够准确预测处于高成岩各阶段各类碎屑岩有效储层的分布；通过成岩阶段的细分、成岩事件的详细恢复、模拟网格的细化，可对研究目标区段有效储层分布进行精细预测，预测单元可达到百米级。

2、理论依据扎实，避免了传统方法中对采样点数据的过份依赖，依据储层本身的演化过程，从过程—响应原理进行模拟，充分体现了过程对结果的控制。

3、定量性好，预测过程具有定量性，可对储层的发育情况进行定量评价。

4、适用性广，对不同地质背景具有较好的适用性；由于在目标实现过程中，充分考虑了不同地质环境的影响，基本上涵盖了目前所知的各种有效碎屑岩储层的形成过程。

5、资料要求相对较低，由于是基于过程进行预测的，而地质过程本身具有延展性，即具有一定的延展范围，地质场中的大量因素可通过相应的与埋深等相关的因素得以获取，因而在本发明中只需要地层格架、沉积相研究成果及个别井钻井成果，数据要求较少，为不同勘探开发阶段具有不同资料基础的评价提供了方便。

具体实施方式：

以下通过实施例对本发明作进一步的说明。

鄂尔多斯盆地山西组为一成岩作用很强的低孔、低渗储层，分为两段，不同地层段沉积环境和岩性也不相同，目前所处成岩作用阶段、发生的成岩作用类型也不相同，有效储层分布也不相同。应用本发明的技术方案对该地区的有效储层分布进行预测。

1、地质资料收集。

收集了本区以下地质资料：

(1)鄂尔多斯盆地上古生界及以上的地层序列及各地质年代持续时间；

(2)鄂尔多斯盆地200口钻遇上古生界的探井的分层数据；

(3)鄂尔多斯盆地200口上古生界探井山西组和石盆子组各段地层的沉积相类型、各段的沉积相分布图；

(4)鄂尔多斯盆地5口上古生界探井山西组和石盆子组各段的热成熟度分析数据、各类胶结物含量、溶蚀孔隙度、残余孔隙度数据；

(5)鄂尔多斯盆地上古生界储层判别标准。

2、原始资料数据库的建立。

按照Petrel格式准备建模数据，同时在EXCEL数据管理器中建立相关的反演数据库，具体数据库包括：

(1)TTI值与成熟度Ro及成岩阶段之间关系理论关系数据库；

(2)上古生界山西组以上地层的地质年代符号与地质年代持续时间数据库；

(3)山西组上段、下段各井的沉积相类型数据库；

(4)山西组沉积相类型及其粒度、成份数据库；

(5)5口典型探井的单井模拟数据库，包括单井分层数据、单井沉积相分析成果数据、单井储层结构参数、单井成熟度分析数据、单井成气时代分析数据、单井孔隙度分析数据、单井孔隙类型及孔隙构成数据。

(6)地质模型数据库，包括200口探井的山西组以上各层的分层数据。

3、以5口探井资料为基础，模拟各井的沉积充填史、热演化史，确定不同时期的地温梯度。

(1)以井点分层数据为基础，回朔井点的沉积过程，恢复沉积史，具体过程如下：

①由分层数据获得各地质单元的顶底深，进而求得各层的地层厚度；

②目标层的地层厚度即为目标层沉积完成时的目标层底面埋深；

③目标层上部相邻地层单元沉积后目标层底面埋深为目标层的地层厚度加上上部相邻层的地层厚度；

④依次求取各地层单元沉积后目标层的底面埋深，直到目前底面埋深，即可恢复地层的埋藏史。

(2)以5口井分层热成熟度分析数据为基础，模拟井点热演化史，确定井点地温梯度，具体过程如下：

①给定初始平均地温梯度3℃/100M；

②计算自沉积始各地质时期的地层温度；

③按10度的温度区间，计算井中有测试数据的层段的热成熟度TTI值；

④将TTI值转化为Ro值；

⑤将计算的Ro值与实验测得的Ro值进行对比，确定两者的差异，调整地温梯度；

⑥重复②～⑤直至计算的Ro与实测的Ro一致，此时的地温梯度即为该井真实地温梯度。

依次对5口探井进行模拟，得到5口井不同地质年代的地温梯度介于2.8-3.5之间。

4、以5口探井资料为基础，模拟井点各沉积相结构相类型，确定不同沉积相的原始孔隙度影响参数，具体过程如下：

(1)读入井点沉积相参数、储层结构参数、原始孔隙度分析参数；

(2)由沉积相参数与储层结构参数判定其结构相；

(3)回归结构相、原始孔隙度关系，确定不同结构相的孔隙影响系数。

按照5口井的模拟结果，最终确定各沉积相的结构相及其影响参数。

5、以5口探井资料为基础，模拟井点的压实演化史，确定压实减孔参数，固化压实减孔模型，具体过程如下：

(1)压实演化史模拟确定压实减孔系数。

按照建立的压实模型，给定压实减孔初始参数fd1、fd2、深度校正系数D1，对模型进行正演运算，得到不同目标层段压实减孔量；将压实减孔量与实测的压实减孔量进行对比，确定压实减孔量的差异；调整压实减孔参数fd1、fd2及深度较正系数D1，直到计算的压实减孔量与实测的压实减孔量吻合；确定此时的压实减孔参数fd1、fd2及深度较正系数D1，实现拟合模型的固化。

6、恢复钙质胶结史、硅质胶结史、粘土矿物胶结史，确定钙质胶结指数、硅质胶结指数、粘土胶结指数，固化胶结减孔模型，具体过程如下：

(1)钙质胶结史恢复与钙质胶结系数确定。

按照建立的模型，给定各胶结高峰初始胶结系数fc_i，从分析井中读入各胶结高峰标准深度，对分析井的钙质胶结作用过程进行模拟，得到不同层段的钙质胶结物的含量，计算出目前钙质胶结减孔量，将之与实测值进行对比，得到钙质胶结减孔量差异；调整胶结系数，直到模拟结果与井中实测结果一致，确定胶结系数，固定模型，即为该区钙质胶结模拟模型。

(2)硅质胶结史恢复与硅质胶结系数确定。

给定热模拟计算的石英胶结系数a、b值，读取该地区的储层的砂岩粒度分布数据、原始孔隙度数据、石英含量数据，埋藏史、地温史计算结果，石英胶结减孔量数据；

按埋藏、热演化史，以1百万年为单位计算从老到新计算石英的析出量，算出目前目标层段的总减孔量；

将这一减孔量结果与实测数值进行对比，确定两者的差别，据此进行石英胶结系数a、b值的调整；

重复以上步骤直至模拟结果与实测结果一致，此时的石英胶结系数a、b值即为该区的石英胶结系数，则石英胶结模型可据此确定。

7、恢复长石等组构溶解史，确定各组构溶解系数，固化各组构溶解模型，具体过程如下：

①读入测试井不同层段的三个高峰溶蚀埋藏深度、标准地温梯度、溶蚀孔隙度值；

②给定初始标准成熟度、标准埋藏深度；

③计算各测试层段的溶蚀量；

④将计算的溶蚀量与测试的溶蚀量进行对比；

⑤根据对比结果调整标准成熟度和标准埋藏深度；

⑥重复③～⑤直到模拟结果与实测结果一致，此时固定相应的参数，即为长石溶蚀模拟具体的模拟模型。

8、读入分层数据、沉积相数据及构造模型，建立地质模型，利用成熟的建模软件(如PETREL)，建立研究区的地质模型，具体过程如下：

①读入地层分层数据，建立不同地质单元的构造模型，体现不同地质单元的地层的顶底面的形态及其变化；

②读入模拟目标层段的沉积相研究成果，建立目标层段沉积相模型。

9、模拟目标区的沉积充填史。

由地质模型中的各层顶底面深度数据，由老到新逐点恢复目标地区地层充填埋藏史。

10、模拟目标区的热演化史，确定目标地层不同地质时期的温度、热成熟度、所处的成岩阶段，具体过程如下：

(1)由地质模型和井点的地温梯度计算各点各地质时期的地温梯度；

(2)依据各地质时期的地温梯度计算目标层各点的热成熟度指标TTI值；

(3)将TTI值转化为Ro值；

(4)由Ro值确定各点所处的成岩作用阶段；

(5)由成岩阶段数据库和成岩阶段与孔隙类型关系判断目标层段在各地质时期的储集空间类型。

11、依据地质模型，模拟目标层段的原始孔隙度，具体过程如下：

(1)读取模型中目标层段各点的沉积相类型；

(2)由沉积相数据计算各类储层的结构相；

(3)按结构相类型判断其原始孔隙度因子；

(4)由原始孔隙度因子与其区域储层孔隙度的值计算各网格点的原始孔隙度；表明目标层段原始孔隙度介于36％-42.2％，平均39％。

12、模拟目标区的压实减孔史，确定不同时期的累积压实减孔量，具体过程如下：

(1)逐个读取目标层段在不同时期的地层埋藏深度值和地温梯度值；

(2)按压实减孔模型，由老到新逐期、逐点计算目标层段不同时代的地层压实减孔量；

(3)由老到新计算不同地质时代累计压实减孔量；表明目前目标层段压实介于18％-22％，平均20.6％。

13、模拟目标区的钙质胶结史、硅质胶结史、粘土胶结史，确定不同时期的累积胶结减孔量，具体过程如下：

(1)逐个读取目标层段在不同时代的地层埋藏深度值和地温梯度值；

(2)按胶结减孔模型(包括钙质、硅质和粘土矿物胶结)由老到新逐期、逐点计算目标层段不同时代的地层胶结减孔量；

(3)由老到新计算不同地质时代累计胶结减孔量，表明目前目标层段胶结减孔量介于12％～18.2％，平均15.6％。

14、模拟目标区的矿物溶解史，确定不同时期的溶解增孔量，具体过程如下：

(1)逐个读取目标层段在不同时期的地层埋藏深度值和地温梯度值；

(2)按溶解增孔模型，由老到新逐期、逐点计算目标层段不同时代的地层溶解增孔量；

(3)由老到新计算不同地质时代累计溶解增孔量，表明目前地层的总溶蚀增孔量介于3％～9％之间，平均为4.5％。

15、计算不同地质时期目标层段的总孔隙度。

依据计算的原生孔隙，加上新增的溶解孔隙量，减去压实、胶结减少的孔隙便可得到高成岩储层目前的总孔隙度，表明目前目标层段地层的总孔隙度介于4％～13％之间，平均为7.9％。

16、按总孔隙度对储层进行评价，确定有效储层的发育区间。

读入储层判定标准，按储层总孔隙度的大小对各点的储层有效性进行判断，确定其有效性，进而由Petrel给出有效储层的发育范围。

Claims (4)

1.一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，其特征在于，具体步骤如下：

(1)地质资料的收集；

(2)原始资料数据库的建立；

(3)以井点资料为基础，模拟约束井的沉积充填史、热源演化史，确定不同时期的地温梯度；

(4)以井点资料为基础，模拟约束点的各沉积相的结构相类型，确定对不同沉积相原始孔隙度的影响参数，建立不同类型储层原始孔隙度的预测模型；

(5)以井点资料为基础，模拟约束点的压实演化史，确定压实减孔参数，固化压实减孔模型；

(6)恢复钙质胶结史、硅质胶结史、粘土矿物胶结史，确定钙质胶结指数、硅质胶结指数、粘土胶结指数，固化胶结减孔模型；

(7)恢复长石等组构溶解史，确定各组构溶解系数，固化各组构溶解模型；

(8)读入分层数据、沉积相数据及构造模型，建立地质模型；

(9)模拟目标区的沉积充填史；

(10)模拟目标区的热演化史，确定目标地层不同地质时间时的温度、热成熟度、所处成岩阶段；

(11)依据地质模型，模拟目标层段的原始孔隙度；

(12)模拟目标区的压实减孔史，确定不同时期的累积压实减孔量；

(13)模拟目标区的钙质胶结史、硅质胶结史、粘土胶结史，确定不同时期累积胶结减孔量；

(14)模拟目标区的溶解史，确定不同时期溶解增孔量；

(15)计算不同地质时期目标层段的总孔隙度；

(16)按总孔隙度对储层进行评价，确定有效储层的发育区间。

2.根据权利要求1所述的一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，其特征在于所述的地质资料的收集包括以下内容：目标地区的地层序列及各地质年代持续时间；目标地区钻井分层数据、目标层段沉积相研究成果图件、不同地层界面构造图；目标层段的热成熟度分析数据、胶结物含量、溶蚀孔隙度、残余孔隙度等数据。

3.根据权利要求1所述的一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，其特征在于所述的原始资料数据库的建立包括以下几个方面：

(1)热成熟度指标TTI值与成熟度Ro值及成岩阶段之间关系的理论关系数据库；

(2)地质年代符号与地质年代持续时间数据库；

(3)沉积相类型、符号序列库数据库；

(4)结构相与原始孔隙度关系数据库；

(5)用于约束模型的单井模拟数据库，包括单井分层数据、单井沉积相分析成果数据、单井储层结构参数、单井成熟度分析数据、单井成油成气时代分析数据、单井孔隙度分析数据、单井孔隙类型及孔隙构成数据；

(6)地质模型数据库，包括研究区钻井分层数据、不同地层界面构造图、目标层段沉积相平面图等。

4.根据权利要求1所述的一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法，其特征在于所述的总孔隙度是原生孔隙中残留孔隙与新生的溶解孔隙的总和，即原生孔隙，加上新增的溶解孔隙量，减去压实、胶结减少的孔隙便可得到高成岩储层目前的孔隙度即：

其中：

为某地质时刻储层总的孔隙度；

为某储层原始孔隙度；

为第i种地质作用增减的孔隙度。