本发明涉及油气勘探开发
技术领域：
，具体地指一种压实减孔量的定量预测方法。
背景技术：
：沉积物在埋藏成岩过程中，岩石体积和孔隙空间逐渐缩小，岩石变得越来越致密，其中机械压实作用是砂岩储集性能的主要控制因素，由机械压实作用引起储层减孔量占储层总孔隙损失量的绝大部分以上，从而极大地影响了砂体的储集物性。影响碎屑岩储层机械压实作用的地质因素除与砂体物质成分、颗粒结构有关，在埋藏过程中成岩环境也是影响机械压实作用的一个重要影响因素。随着成岩过程的进行，各种成岩地质因素互相影响，影响程度表现不一。在沉积物埋藏较浅的成岩早期，影响机械压实作用的主要地质因素为岩石粒度、分选、刚性颗粒、磨圆、泥质含量等，到了成岩中期，影响机械压实作用的主要地质因素除了与组成碎屑岩石颗粒的原始组构之外，成岩过程中胶结类型、胶结强度、胶结时间、异常高压等地质因素也是影响压实强度的重要因素。到了成岩的晚期，岩石已固结成岩，机械压实作用对原生孔隙改造几乎不起作用。从影响机械压实作用的因素出发，通过物理模拟实验来探讨机械压实作用过程中砂体储层物性的变化规律，由于受实验条件限制，模拟参数选择多为温度、压力、岩石组分以及粒度大小，而不能综合考虑多种因素作用下机械压实作用过程中砂体储层物性的变化规律。由于不能系统开展不同岩石成分、不同粒级大小和不同分选程度的单体或混合矿物的压实模拟试验，因此就无法进一步揭示压实过程中碎屑岩储层物性的变化规律，也就无法准确、定量预测压实作用下原生孔隙减少量。通过计算机模拟，可以综合考虑埋深、地质年代、热成熟度、温度、沉积速率、石英含量、胶结程度、构造旋回等地质因素对机械压实作用的影响，并通过因子分析、主次因素分析和多项式拟合，建立砂岩压实减孔量与埋深、颗粒密度、分选性、原始孔隙度、温度等成岩要素之间模型公式，对砂岩储层物性预测能起到良好的指导作用，然而其科学性、合理性、真实性、准确性需依赖合理数学地质模型的建立。技术实现要素：本发明针对上述存在的问题，提供了一种压实减孔量的定量预测方法，该方法综合考虑了沉积相、粒度大小、岩屑颗粒、刚性颗粒、磨圆度、分选性、泥质含量、胶结强度、胶结类型、胶结时间、埋藏深度、异常高压等多种地质因素对压实作用影响，分别在时空域内进行压实作用数值模拟，构建了一个能够反映沉积物在埋藏成岩过程中孔隙演化与压实作用之间定量化模型。该模型科学、合理、准确，为压实减孔量的定量预测提供了一种新的模拟方法。为了实现上述目的，本发明提供了一种压实减孔量的定量预测方法，包括以下步骤：1)收集研究区地质资料2)建立原始资料数据库(1)建立沉积相类型数据库研究不同沉积相类型控制下碎屑岩储层原始组构以及原生孔隙发育情况，即，“沉积相-原始组构”模型和“原始组构-孔隙度”模型(沉积相-原始组构”模型和“原始组构-孔隙度”模型为地质历史时期原生孔隙演化提供基础资料)，建立沉积相类型数据库；(2)建立地质参数qi数据库a.建立粒度(d_a)数据库，b.建立分选性(s_b)数据库，c.建立磨圆度(r_c)数据库，d.建立刚性颗粒(v_d)数据库，e.建立矿物组分(m_e)数据库，f.建立矿物比分(o_f)数据库，g.建立胶结类型(c_g)数据库，h.建立胶结强度(q_h)数据库，i.建立胶结时间(t_i)数据库，j.建立溶蚀时间(f_j)数据库，k.建立泥质含量(n_k)数据库，l.建立岩屑含量(x_l)数据库，m.建立沉积速率(e_m)数据库，n.超压(p_n)数据库；3)基于上述原始资料数据库建立基础模型(1)建立沉积相与碎屑岩原始组构模型f(sedimentaryfacies)＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c)其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；(2)建立原始组构与沉积物原始孔隙度模型；其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；(3)建立压实指数与成岩参数(沉积相、粒度大小、刚性颗粒、磨圆度、分选性、泥质含量、胶结强度、胶结类型、胶结时间、沉积速率、异常高压)模型；pi＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c，c_g，q_h，t_i，f_j，e_m，p_n，dep)其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c，胶结类型c_g，胶结强度q_h，胶结时间t_i，溶蚀时间f_j，沉积速率e_m，超压p_n；(4)建立压实模型cm建立压实模型cm，其中压实模型cm包括ⅰ、ⅱ、ⅲ、ⅳ、ⅴ、ⅵ、ⅶ、ⅷ、ⅸ、ⅹ总共十种模型；ⅰ:f1＝1.9467*ln(dep)-0.8689；ⅱ:f2＝2.9200*ln(dep)-1.3034；ⅲ:f3＝4.3800*ln(dep)-1.9551；ⅳ:f4＝7.0080*ln(dep)-3.1282；ⅴ:f5＝8.7600*ln(dep)-3.9103；ⅵ:f6＝11.6800*ln(dep)-5.2137；ⅶ:f7＝14.0160*ln(dep)-6.2564；ⅷ:f8＝17.5200*ln(dep)-7.8205；ⅸ:f9＝20.6118*ln(dep)-9.2006；ⅹ:f10＝23.3600*ln(dep)-10.4273；其中，fi为压实率，dep为埋深；在储层成岩的过程中，可根据压实指数pi确定储层演化属于的压实模型；pi压实指数，n为地质参数的个数，qi为地质参数，其大小根据压实级别来定，pi为第i个压实指标的权值；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；4)影响参数选取(1)研究区储层网格化将研究区的储层进行网格化，每个网格化的储层的表示为wm(xm，ym，zm)；(2)确定每个网格化的储层沉积相根据上述建立的沉积相类型数据库，分析每个网格化的储层的沉积相类型。wm(xm，ym，zm)＝f(sedimentaryfacies)(3)确定每个网格化的储层的原始组构根据研究区中每个网格化的储层的沉积相数据确定每个空间网格原始组构数据，并根据研究区已知井“沉积相-原始组构数据”对未知区域“沉积相-原始组构”数据f(sedimentaryfacies)＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c)进行校正；(4)确定每个网格化的储层的原始孔隙a.确定每个网格原始组构其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；b.在步骤a基础上，确定每个网格的压实级别g当15≥g≥7，当20≥g≥16，当30≥g≥21，当35≥g≥31，当40≥g≥36，(5)确定演化期次根据研究区埋藏史数据，确定每个网格化的储层即为网格化的目标储层在地质历史时期演化期次，演化期次以tj表示j＝1，2，3……；(6)确定网格化的目标储层的胶结类型和胶结强度、溶蚀时间；根据研究区成岩相数据以及单井成岩序列确定网格化的目标储层每个期次胶结类型和胶结强度、溶蚀时间；(7)确定网格化的目标储层的压实指数以及每个网格初始压实模型im在储层埋藏之前，决定储层初始压实模型im的为分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m地质参数，根据以上地质参数qi确定每个网格初始压实模型im，初始压实模型的确定可用以下公式：计算初始压实模型的地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；(8)确定网格化的目标储层不同演化阶段地质指数根据地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m、胶结类型c_g、胶结强度q_h、胶结时间t_i、溶蚀时间f_j、超压p_n确定成岩过程中压实指数pi。压实指数计算公式如下：计算压实指数pi的地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m、胶结类型c_g、胶结强度q_h、胶结时间t_i、溶蚀时间f_j、超压p_n；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；(9)确定网格化的目标储层不同演化阶段压实模型及孔隙度大小gj＝fi；其中，gj为演化阶段j的压实率改变量，fi为相应的压实率模型，为j-1阶段孔隙度，为第j阶段孔隙度，i＝1，2，3，……10，为原始孔隙度；(10)计算孔隙度随时间演化分布a.根据原始组构确定每个网格的压实级别g当15≥g≥7，当20≥g≥16，当30≥g≥21，当35≥g≥31，当40≥g≥36，b.根据研究区地质资料确定每个网格的演化阶段tj，j＝1，2，3……；c.根据原始组构确定每个网格初始压实模型im；d.根据成岩作用以及原始组构确定每个网格不同演化阶段tj地质指数；e.在步骤c和步骤d的基础上确定每个网格不同演化阶段tj压实模型gj；f.根据初始孔隙度初始压实模型im，不同演化阶段压实模型gj计算不同演化阶段tj时机械压实作用下储层孔隙减少量。进一步地，所述步骤1)中，收集研究区的地质资料包括：①研究区目标层沉积相研究成果图；②目标层界面构造等值线图；③目标层上覆地层界面构造等值线图；④研究区储层物性资料；⑤研究区储层成岩相研究成果图；⑥研究区成岩演化序列；⑦研究区埋藏史、热演化史。再进一步地，所述步骤2)的第(2)小步中，a.建立粒度(d_a)数据库按照碎屑岩粒度大小，建立粒度数据库；表碎屑岩粒度(d_a)类型及其分类粒级划分数据库代码颗粒直径(mm)压实级别(g_1)巨砾d_1＞10001粗砾d_2＞100～10002中砾d_3＞10～1003细砾d_4＞1～104粗砂d_5＞0.5～15中砂d_6＞0.25～0.56细砂d_7＞0.1～0.257粗粉砂d_8＞0.05～0.18细粉砂d_9＞0.01～0.059粘土d_10≤0.0110b.建立分选性(s_b)数据库按照碎屑岩颗粒大小的标准偏差即均一性，建立碎屑岩分选性数据库；表碎屑岩分选性(s_b)类型及其分类分选性级别数据库代码标准偏差压实级别(g_2)分选极好s_1≤0.351分选好s_2＞0.35～0.502分选较好s_3＞0.50～0.713分选中等s_4＞0.71～1.004分选较差s_5＞1.00～2.005分选差s_6＞2.00～4.006分选极差s_7＞4.007c.建立磨圆度(r_c)数据库按照碎屑岩磨圆情况，建立磨圆度数据库；表碎屑岩磨圆度(r_c)类型及其分类磨圆度数据库代码压实级别(g_3)棱角状r_11次棱角状r_22次圆状r_33圆状r_44d.建立刚性颗粒(v_d)数据库按照碎屑岩中刚性颗粒含量，建立刚性颗粒数据库；表碎屑岩刚性组分(v_d)及其分类刚性颗粒数据库代码压实级别(g_4)＞75v_11＞60～75v_22＞45～60v_33＞30～45v_44≤30v_55e.建立矿物组分(m_e)数据库按照碎屑岩中矿物组分，建立矿物组分数据库。表碎屑岩矿物组分(m_e)及其分类矿物组分数据库代码压实级别(g_5)石英、长石、方解石、白云石、云母、黏土矿物m_11石英、方解石、白云石m_22石英、方解石、白云石、黏土矿物m_33石英、长石m_44长石、方解石、白云石m_55石英、方解石m_66f.建立矿物比分(o_f)数据库按照碎屑岩中长石含量建立矿物比分数据库；表碎屑岩矿物比分(mn)及其分类长石百分比(％)数据库代码压实级别(g_6)＞75o_11＞60～75o_22＞45～60o_33＞30～45o_44≤30o_55g.建立胶结类型(c_g)数据库按照碎屑岩胶结类型，建立胶结类型数据库；表胶结类型(c_g)数据库胶结类型数据库代码压实级别(g_7)钙质c_11硅质c_22铁质c_33黏土c_44h.建立胶结强度(q_h)数据库按照碎屑岩的胶结强度，建立胶结强度(q_h)数据库；表8胶结类型(q_h)数据库胶结强度数据库代码压实级别(g_8)弱(胶结率≤2)q_11中(7≥胶结率＞2)q_22强(胶结率＞7)q_33i.建立胶结时间(t_i)数据库按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；表胶结时间(t_i)数据库胶结时间数据库代码压实级别(g_9)ⅰat_15ⅰbt_24ⅱat_33ⅱbt_42ⅲt_51j.建立溶蚀时间(f_j)数据库按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；表溶蚀时间(f_j)数据库溶蚀时间数据库代码压实级别(g_10)ⅰaf_11ⅰbf_22ⅱaf_33ⅱbf_44ⅲf_55k.建立泥质含量(n_k)数据库按照碎屑岩中的泥质含量，建立泥质含量数据库；表泥质含量n_k分类泥质含量(％)数据库代码压实级别(g_11)＜2n_11≥2～5n_22＞5～8n_33＞8n_44l.建立岩屑含量(x_l)数据库按照碎屑岩中岩屑含量，建立岩屑含量数据库。表岩屑含量(x_l)分类岩屑含量(％)数据库代码压实级别(g_12)≤2x_11＞2～6x_22＞6～10x_33＞10x_44m.建立沉积速率(e_m)数据库按照碎屑岩的沉积速率将碎屑岩储层进行分类，建立沉积速率数据库；表沉积速率(e_n)分类沉积速率(m/ma)数据库代码压实级别(g_13)＜100e_11≥100～500e_22＞500～1000e_33＞1000e_44n.超压(p_n)数据库按照压力系数大小将碎屑岩储层演化过程中压力大小进行分类，建立超压数据库；压力系数αp＝pp/ph其中，pp为实测地层压力，ph为静水压力；表超压(p_n)分类超压程度αp数据库代码压实级别(g_14)弱超压＜1.2p_11中超压1.2～1.8p_22强压＞1.8p_33本发明的有益效果在于：(1)本发明在压实减孔量的定量预测上，综合考虑了多种地质参数，其中包括沉积相、粒度大小、岩屑颗粒、刚性颗粒、磨圆度、分选性、泥质含量、胶结强度、胶结类型、胶结时间、溶蚀时间、埋藏深度、异常高压、沉积速率，弥补了计算过程中参数选择不足；(2)本发明在压实模型的建立上，使用了十种理想压实曲线对成岩过程中压实作用进行动态刻画，能准确地恢复整个压实过程；(3)本发明在地质参数的选取上，以动态形式进行综合判识。附图说明图1为研究区埋藏史、温度史、成岩演化序列；具体实施方式为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。实施例1压实减孔量的定量预测方法，包括以下步骤：1)收集研究区地质资料，它包括：①研究区目标层沉积相研究成果图；②目标层界面构造等值线图；③目标层上覆地层界面构造等值线图；④研究区储层物性资料；⑤研究区储层成岩相研究成果图；⑥研究区成岩演化序列；⑦研究区埋藏史、热演化史；2)建立原始资料数据库(1)建立沉积相类型数据库研究不同沉积相类型控制下碎屑岩储层原始组构以及原生孔隙发育情况，即，“沉积相-原始组构”模型和“原始组构-孔隙度”模型(沉积相-原始组构”模型和“原始组构-孔隙度”模型为地质历史时期原生孔隙演化提供基础资料)，建立沉积相类型数据库；(2)建立地质参数qi数据库a.建立粒度(d_a)数据库按照碎屑岩粒度大小，建立粒度数据库；表1碎屑岩粒度(d_a)类型及其分类粒级划分数据库代码颗粒直径(mm)压实级别(g_1)巨砾d_1＞10001粗砾d_2＞100～10002中砾d_3＞10～1003细砾d_4＞1～104粗砂d_5＞0.5～15中砂d_6＞0.25～0.56细砂d_7＞0.1～0.257粗粉砂d_8＞0.05～0.18细粉砂d_9＞0.01～0.059粘土d_10≤0.0110b.建立分选性(s_b)数据库按照碎屑岩颗粒大小的标准偏差即均一性，建立碎屑岩分选性数据库；表2碎屑岩分选性(s_b)类型及其分类分选性级别数据库代码标准偏差压实级别(g_2)分选极好s_1≤0.351分选好s_2＞0.35～0.502分选较好s_3＞0.50～0.713分选中等s_4＞0.71～1.004分选较差s_5＞1.00～2.005分选差s_6＞2.00～4.006分选极差s_7＞4.007c.建立磨圆度(r_c)数据库按照碎屑岩磨圆情况，建立磨圆度数据库；表3碎屑岩磨圆度(r_c)类型及其分类磨圆度数据库代码压实级别(g_3)棱角状r_11次棱角状r_22次圆状r_33圆状r_44d.建立刚性颗粒(v_d)数据库按照碎屑岩中刚性颗粒含量，建立刚性颗粒数据库；表4碎屑岩刚性组分(v_d)及其分类刚性颗粒数据库代码压实级别(g_4)＞75v_11＞60～75v_22＞45～60v_33＞30～45v_44≤30v_55e.建立矿物组分(m_e)数据库按照碎屑岩中矿物组分，建立矿物组分数据库。表5碎屑岩矿物组分(m_e)及其分类矿物组分数据库代码压实级别(g_5)石英、长石、方解石、白云石、云母、黏土矿物m_11石英、方解石、白云石m_22石英、方解石、白云石、黏土矿物m_33石英、长石m_44长石、方解石、白云石m_55石英、方解石m_66f.建立矿物比分(o_f)数据库按照碎屑岩中长石含量建立矿物比分数据库；表6碎屑岩矿物比分(mn)及其分类长石百分比(％)数据库代码压实级别(g_6)＞75o_11＞60～75o_22＞45～60o_33＞30～45o_44≤30o_55g.建立胶结类型(c_g)数据库按照碎屑岩胶结类型，建立胶结类型数据库；表7胶结类型(c_g)数据库h.建立胶结强度(q_h)数据库按照碎屑岩的胶结强度，建立胶结强度(q_h)数据库；表8胶结类型(q_h)数据库胶结强度数据库代码压实级别(g_8)弱(胶结率≤2)q_11中(7≥胶结率＞2)q_22强(胶结率＞7)q_33i.建立胶结时间(t_i)数据库按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；表9胶结时间(t_i)数据库胶结时间数据库代码压实级别(g_9)ⅰat_15ⅰbt_24ⅱat_33ⅱbt_42ⅲt_51j.建立溶蚀时间(f_j)数据库按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；表10溶蚀时间(f_j)数据库溶蚀时间数据库代码压实级别(g_10)ⅰaf_11ⅰbf_22ⅱaf_33ⅱbf_44ⅲf_55k.建立泥质含量(n_k)数据库按照碎屑岩中的泥质含量，建立泥质含量数据库；表11泥质含量n_k分类泥质含量(％)数据库代码压实级别(g_11)＜2n_11≥2～5n_22＞5～8n_33＞8n_44l.建立岩屑含量(x_l)数据库按照碎屑岩中岩屑含量，建立岩屑含量数据库。表12岩屑含量(x_l)分类岩屑含量(％)数据库代码压实级别(g_12)≤2x_11＞2～6x_22＞6～10x_33＞10x_44m.建立沉积速率(e_m)数据库按照碎屑岩的沉积速率将碎屑岩储层进行分类，建立沉积速率数据库；表13沉积速率(e_n)分类沉积速率(m/ma)数据库代码压实级别(g_13)＜100e_11≥100～500e_22＞500～1000e_33＞1000e_44n.超压(p_n)数据库按照压力系数大小将碎屑岩储层演化过程中压力大小进行分类，建立超压数据库；压力系数αp＝pp/ph其中，pp为实测地层压力，ph为静水压力；表14超压(p_n)分类超压程度αp数据库代码压实级别(g_14)弱超压＜1.2p_11中超压1.2～1.8p_22强压＞1.8p_333)基于上述原始资料数据库建立基础模型(1)建立沉积相与碎屑岩原始组构模型f(sedimentaryfacies)＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c)其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；(2)建立原始组构与沉积物原始孔隙度模型；其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；(3)建立压实指数与成岩参数(沉积相、粒度大小、刚性颗粒、磨圆度、分选性、泥质含量、胶结强度、胶结类型、胶结时间、沉积速率、异常高压)模型；pi＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c，c_g，q_h，t_i，f_j，e_m，p_n，dep)其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c，胶结类型c_g，胶结强度q_h，胶结时间t_i，溶蚀时间f_j，沉积速率e_m，超压p_n；(4)建立压实模型cm建立压实模型cm，其中压实模型cm包括ⅰ、ⅱ、ⅲ、ⅳ、ⅴ、ⅵ、ⅶ、ⅷ、ⅸ、ⅹ总共十种模型；ⅰ:f1＝1.9467*ln(dep)-0.8689；ⅱ:f2＝2.9200*ln(dep)-1.3034；ⅲ:f3＝4.3800*ln(dep)-1.9551；ⅳ:f4＝7.0080*ln(dep)-3.1282；ⅴ:f5＝8.7600*ln(dep)-3.9103；ⅵ:f6＝11.6800*ln(dep)-5.2137；ⅶ:f7＝14.0160*ln(dep)-6.2564；ⅷ:f8＝17.5200*ln(dep)-7.8205；ⅸ:f9＝20.6118*ln(dep)-9.2006；ⅹ:f10＝23.3600*ln(dep)-10.4273；其中，fi为压实率，dep为埋深；在储层成岩的过程中，可根据压实指数pi确定储层演化属于的压实模型；pi压实指数，n为地质参数的个数，qi为地质参数，其大小根据压实级别来定，pi为第i个压实指标的权值；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；4)影响参数选取(1)研究区储层网格化将研究区的储层进行网格化，每个网格化的储层的表示为wm(xm，ym，zm)；(2)确定每个网格化的储层沉积相根据上述建立的沉积相类型数据库，分析每个网格化的储层的沉积相类型。wm(xm，ym，zm)＝f(sedimentaryfacies)(3)确定每个网格化的储层的原始组构根据研究区中每个网格化的储层的沉积相数据确定每个空间网格原始组构数据，并根据研究区已知井“沉积相-原始组构数据”对未知区域“沉积相-原始组构”数据f(sedimentaryfacies)＝f(s_b，d_a，x_l，n_k，v_d，m_e，r_c)进行校正；(4)确定每个网格化的储层的原始孔隙a.确定每个网格原始组构其中：分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c；b.在步骤a基础上，确定每个网格的压实级别g当15≥g≥7，当20≥g≥16，当30≥g≥21，当35≥g≥31，当40≥g≥36，(5)确定演化期次根据研究区埋藏史数据，确定每个网格化的储层即为网格化的目标储层在地质历史时期演化期次，演化期次以tj表示j＝1，2，3……；(6)确定网格化的目标储层的胶结类型和胶结强度、溶蚀时间；根据研究区成岩相数据以及单井成岩序列确定网格化的目标储层每个期次胶结类型和胶结强度、溶蚀时间；(7)确定网格化的目标储层的压实指数以及每个网格初始压实模型im在储层埋藏之前，决定储层初始压实模型im的为分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m地质参数，根据以上地质参数qi确定每个网格初始压实模型im，初始压实模型的确定可用以下公式：计算初始压实模型的地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；(8)确定网格化的目标储层不同演化阶段地质指数根据地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m、胶结类型c_g、胶结强度q_h、胶结时间t_i、溶蚀时间f_j、超压p_n确定成岩过程中压实指数pi。压实指数计算公式如下：计算压实指数pi的地质参数qi包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m、胶结类型c_g、胶结强度q_h、胶结时间t_i、溶蚀时间f_j、超压p_n；pi∈[18～12)，cm＝ⅰ；pi∈[12～8)，cm＝ⅱ；pi∈[8～5)，cm＝ⅲ；pi∈[5～4)，cm＝ⅳ；pi∈[4～3)，cm＝ⅴ；pi∈[3～2.5)，cm＝ⅵ；pi∈[2.5～2)，cm＝ⅶ；pi∈[2～1.7)，cm＝ⅷ；pi∈[1.7～1.5)，cm＝ⅸ；pi∈[1.5～0)，cm＝ⅹ；(9)确定网格化的目标储层不同演化阶段压实模型及孔隙度大小gj＝fi；其中，gj为演化阶段j的压实率改变量，fi为相应的压实率模型，为j-1阶段孔隙度，为第j阶段孔隙度，i＝1，2，3，……10，为原始孔隙度；(10)计算孔隙度随时间演化分布a.根据原始组构确定每个网格的压实级别g当15≥g≥7，当20≥g≥16，当30≥g≥21，当35≥g≥31，当40≥g≥36，b.根据研究区地质资料确定每个网格的演化阶段tj，j＝1，2，3……；c.根据原始组构确定每个网格初始压实模型im；d.根据成岩作用以及原始组构确定每个网格不同演化阶段tj地质指数；e.在步骤c和步骤d的基础上确定每个网格不同演化阶段tj压实模型gj；f.根据初始孔隙度初始压实模型im，不同演化阶段压实模型gj计算不同演化阶段tj时机械压实作用下储层孔隙减少量。实施例2基于上述压实减孔量的定量预测方法对渤中凹陷bz19-2-1进行预测：分析bz19-2-1井东二段d2为滨浅湖相，岩性为中粒石英质长石砂岩，岩石分选性较好，粒度次圆状，刚性颗粒平均含量68％，矿物组分以长石、石英为主，碎屑岩中长石含量平均为53％。东二段d2在整个地质历史时期经历了三个演化阶段，分别是：ⅰa、ⅰb、ⅱa。在ⅰa演化阶段，胶结类型为泥质胶结和钙质胶结，胶结强度中等；在ⅰb演化阶段，胶结类型为泥质胶结和钙质胶结，胶结强度中等；在ⅱa演化阶段，胶结类型为泥质胶结，胶结强度弱，且长石发生了溶蚀作用，溶蚀强度中等。整个地质演化过程中，地层无异常高压现象；具体预测方法如下：1、收集研究区地质资料(1)沉积相；(2)研究区埋藏史、热演化史、成岩演化序列；2、建立基础地质模型(1)确定原始组构根据沉积相-原始组构模型，结合研究区沉积相资料、原始组构资料，确定井bz19-2-1东营组二段储层滨湖相下碎屑颗粒粒度d_a、分选性s_b、磨圆度r_c、刚性颗粒v_d、矿物组分m_e、矿物比分o_f、泥质含量n_k、岩屑含量x_l的分布，即：f(ⅲ-a-b)＝f(d_a，s_b，r_c，v_d，m_e，o_f，n_k，x_l)其中：d_a＝d_6，s_b＝s_3，r_c＝r_3，v_d＝v_2，m_e＝m_4，o_f＝o_3，n_k＝n_1，x_l＝x_1；(2)确定压实级别根据原始组构-原始孔隙模型，选择地质参数，结合压实级别计算公式，确定井bz19-2-1东营组二段储层的原始孔隙度；其中，d_a＝d_6＝6，s_b＝s_3＝3，r_c＝r_3＝3，v_d＝v_2＝2，m_e＝m_4＝4，o_f＝o_3＝3，n_k＝n_1＝1，x_l＝x_1＝1；即g＝23，压实级别计算公式(3)确定演化期次根据研究区埋藏史数据，确定bz19-2-1东营组二段储层在地质历史时期演化期次，其由老到新的演化顺序依次是：t1＝ⅰa，t2＝ⅰb，t3＝ⅱa；(4)确定初始压实模型选择参数，其中包括分选性s_b、粒度d_a、岩屑含量x_l、泥质含量n_k、刚性颗粒v_d、磨圆度r_c、矿物比分o_f、矿物组分m_e、泥质含量n_k、沉积速率e_m地质参数，进而确定压实指数。压实指数计算公式：根据压实指数确定bz19-2-1东营组二段储层压实初始模型im。其中，d_a＝d_6＝6，s_b＝s_3＝3，r_c＝r_3＝3，v_d＝v_2＝2，m_e＝m_4＝4，o_f＝o_3＝3，n_k＝n_1＝1，x_l＝x_1＝1，e_m＝3，pi＝0.15，即pi＝3.9；im:g0＝f5＝8.7600*ln(dep)-3.9103(5)确定不同阶段压实指数根据地质参数，其中包括分选性s_b，粒度d_a，岩屑含量x_l，泥质含量n_k，刚性颗粒v_d，矿物组分m_e，磨圆度r_c，胶结类型c_g，胶结强度q_h，胶结时间t_i，溶蚀时间f_j，沉积速率e_m，超压p_n，进而确定不同阶段储层压实指数。压实指数计算公式：a.ⅰa阶段压实指数以及压实模型在ⅰa演化阶段，d_a＝d_6＝6，s_b＝s_3＝3，r_c＝r_3＝3，v_d＝v_2＝2，m_e＝m_4＝4，o_f＝o_3＝3，n_k＝n_1＝1，x_l＝x_1＝1，e_m＝3，胶结类型c_g＝2，胶结强度q_h＝2，胶结时间t_i＝5，溶蚀时间f_j＝0，超压p_n＝0，pi＝0.15，即pi＝5.25；im:g0＝f3＝4.3800*ln(dep)-1.9551；b.ⅰb阶段压实指数以及压实模型在ⅰb演化阶段，d_a＝d_6＝6，s_b＝s_3＝3，r_c＝r_3＝3，v_d＝v_2＝2，m_e＝m_4＝4，o_f＝o_3＝3，n_k＝n_1＝1，x_l＝x_1＝1，e_m＝3，胶结类型c_g＝2，胶结强度q_h＝2，胶结时间t_i＝5，溶蚀时间f_j＝0，超压p_n＝0，pi＝0.20，即pi＝7.0；im:g0＝f3＝4.3800*ln(dep)-1.9551；c.ⅱa阶段压实指数以及压实模型在ⅱa演化阶段，d_a＝d_6＝6，s_b＝s_3＝3，r_c＝r_3＝3，v_d＝v_2＝2，m_e＝m_4＝4，o_f＝o_3＝3，n_k＝n_1＝1，x_l＝x_1＝1，e_m＝3，胶结类型c_g＝4，胶结强度q_h＝1，胶结时间t_i＝3，溶蚀时间f_j＝3，超压p_n＝0，pi＝0.20，即pi＝8.14；im:g0＝f2＝2.9200*ln(dep)-1.3034；(6)确定不同演化阶段机械压实作用下减孔量及孔隙度大小a.目标层初始压实模型压实模型：im:g0＝f5＝8.7600*ln(dep)-3.9103，dep∈[0，150)；孔隙度：孔隙度改变量：b.目标层ⅰa演化阶段压实模型：g1＝f3＝4.3800*ln(dep)-1.9551，dep∈[150，1021)；孔隙度：孔隙度改变量：c.目标层ⅰb演化阶段压实模型：g2＝f3＝4.3800*ln(dep)-1.9551，dep∈[1021，1760)；孔隙度：孔隙度改变量：d.目标层ⅱa演化阶段压实模型：g3＝f2＝2.9200*ln(dep)-1.3034，dep∈[1760，3982)；孔隙度：孔隙度改变量：其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。当前第1页12