本发明涉及油气勘探开发技术领域，具体地指一种碎屑岩储层胶结减孔量的定量预测方法。
背景技术：
沉积物在埋藏成岩过程中，岩石体积和孔隙空间逐渐缩小，岩石变得越来越致密，其中压实作用是砂岩储集性能的主要控制因素，由压实作用引起储层减孔量占储层总孔隙损失量的绝大部分以上，从而极大地影响了砂体的储集物性。
胶结作用是沉积物在成岩过程中的一种变化，在此过程中，从孔隙溶液里析出晶体能将松散颗粒粘在一起，固结成岩。常见的胶结物有泥质的、铁质的、硅质的和钙质的。在储层演化中，胶结作用是导致砂体内原生孔隙减少的主要原因之一，从而极大地影响了砂体的储集物性。另外，砂岩中很多胶结作用还是重要的保持性成岩作用，是深埋藏条件下砂岩孔隙得以保持的重要机制，这在储层质量预测中具有重要意义。
目前，在研究碎屑岩储层成岩作用的时候，多是关注单矿物方面胶结作用，且多数研究是简单的、定性的描述，而没有开展系统的、深入研究。由于成岩过程中，胶结类型有多有少，胶结时间有先有后，胶结的强度有强有弱，从而使得定量预测碎屑岩储层胶结量存在很大难度。如果通过查清不同地质过程的关键地质因素及其差异以及不同地质过程中各地质因素对储层孔隙度的影响，构建孔隙响应定量方程，依据定量响应方程模拟储层孔隙演化过程，再现沉积物埋藏成岩过程中胶结作用对储层物性变化规律影响，就能对碎屑岩胶结减孔量进行很好定量预测，从而为盆地储层物性研究提供新的思路和依据。
技术实现要素：
本发明针对上述存在的问题，提供了一种碎屑岩储层胶结减孔量的定量预测方法，该方法基于碎屑岩成岩过程中不同演化阶段成岩作用研究，建立起不同成岩作用下胶结模型，确定储层演化过程中原始孔隙减少情况，为油藏评价提供依据。
为实现上述目的，本发明提供的一种碎屑岩储层胶结减孔量的定量预测方法，包括以下步骤：
1)收集研究区地质资料
2)建立原始资料数据库，它包括
(1)沉积相类型数据库f_m
(2)岩性数据库r_n
(3)酸碱度数据库p_o
(4)沉积旋回数据库c_p
(5)成岩阶段数据库s_q
(6)断层指数数据库n_r
3)建立胶结作用数据库d_c
根据岩石中不同矿物mi的孔隙改变量建立胶结作用数据库；
4)建立成岩作用预测模型
(1)研究区储层网格化
对研究区储层网格化，即为将空间上不均匀分布的数据按一定方法(如滑动平均法、克里格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程；研究区储层的每个网格用wi(x，y)表示；
(2)确定网格属性
a.根据研究区沉积相数据确定网格wi(x，y)的沉积相属性f_m；
b.根据研究区岩性数据确定网格wi(x，y)的岩性属性r_n；
c.根据研究区酸碱度数据确定网格wi(x，y)的酸碱度属性p_o；
d.根据研究区沉积旋回数据确定网格wi(x，y)的旋回属性c_p；
e.根据研究区成岩阶段数据确定网格wi(x，y)的成岩阶段属性s_q；
f.根据研究区断层指数数据确定网格wi(x，y)的断层指数属性n_r；
g.根据上述步骤a、b、c、d、e和f确定每个网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)，即：
wi(x，y)＝f_m+r_n+p_o+c_p+s_q+n_r；
(3)基于研究区岩石薄片资料、分析测试资料，确定研究区已知井hj所在网格wj(x，y)的成岩作用d_s_j；
(4)确定已知井hj所在网格wj(x，y)的综合属性gj(f_m_j，r_n_j，p_o_j，c_p_j，s_q_j，n_r_j)，建立已知井hj所在网格wj(x，y)成岩作用d_s与综合属性g(f_m_j，r_n_j，p_o_j，c_p_j，s_q_j，n_r_j)的对应关系，即为：
d_s_j＝f_m_j+r_n_j+p_o_j+c_p_j+s_q_j+n_r_j；
(5)任取一未知网格wi(x，y)，确定未知网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)，即为：
gi＝(f_m_i，r_n_i，p_o_i，c_p_i，s_q_i，n_r_i)；
(6)将网格wi(x，y)的成岩作用综合属性gi与井hj所在网格wj(x，y)的综合属性gj相比较，即为：
gi-gj＝(f_m_i，r_n_i，p_o_i，c_p_i，s_q_i，n_r_i)-(f_m_j+r_n_j+p_o_j+c_p_j+s_q_j+n_r_j)；
如果满足gi-gj＝0，则未知网格wi(x，y)的成岩作用d_s_i与井j具有相同的成岩作用d_s_j；
如果gi-gj≠0，则按不同属性优先级顺序，即：沉积相f_m一级、岩性r_n二级、沉积旋回c_p三级、酸碱度p_o四级、成岩阶段s_q五级、断层指数n_r六级对未知网格wi(x，y)的成岩作用进行判识，即为：
ⅰ：f_m_i-f_m_j＝0，r_n_i-r_n_j≠0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅱ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅲ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅳ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅴ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j＝0，n_r_i-n_r_j≠0；
未知网格wi(x，y)的成岩作用d_s_i满足条件i时，未知网格wi(x，y)与具有相同属性f_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅱ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅲ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅳ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、s_q_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅴ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、p_o_j、s_q_j的网格wj(x，y)的成岩作用相同；
(7)根据每个网格wi(x，y)的成岩作用确定研究区成岩作用横向分布；
5)建立胶结减孔量数学模型
建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
a.石英胶结模型d1
其中：a，b为常数，a＝1.98×10-22，b＝0.022℃(-1)；vque_cem为石英胶结作用改变孔隙体积百分量，％；m为石英的摩尔质量，g/mol；ρ为石英的密度，g/cm3；t为反应温度，℃；t＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；a0为自生石英开始形成时石英颗粒原始表面积；a0＝6fv/d；d为石英颗粒的粒径，mm；v为单位砂岩体积，mm3；f为单位体积砂岩中石英的体积百分含量；c为平均每年地层温度的平均改变量，℃/y；石英强胶结
b.长石胶结模型d2
其中：x埋藏深度(km)；vfel为长石胶结总量，％；长石强胶结
c.方解石胶结模型d3
其中：x埋藏深度(km)；vcal为方解石胶结总量，％；方解石强胶结
d.白云石胶结模型d4
其中：x埋藏深度(km)；vdol为白云石胶结总量，％；白云石强胶结
e.粘土胶结模型d5
其中：x埋藏深度(km)；vcla为胶结总量，％；粘土强胶结
6)胶结量模拟计算
(1)根据研究区埋藏史确定每个网格不同时刻t埋深h；
(2)根据研究区温度史确定每个网格不同时刻t温度t；
(3)根据步骤4)中成岩作用预测模型确定每个网格wi(x，y)胶结减孔量数学模型；
(4)根据步骤(1)中埋深h、步骤(2)中温度t、步骤(3)中胶结减孔量数学模型及其步骤4)中的成岩作用预测模型计算每个网格的胶结量。
进一步地，所述步骤1)中，收集研究区地质资料包括
①研究区目标层沉积相研究成果图；
②目标层界面构造等值线图；
③目标层上覆地层界面构造等值线图；
④研究区储层物性资料；
⑤研究区储层单井成岩作用数据；
⑥研究区成岩演化序列；
⑦研究区埋藏史、热演化史。
再进一步地，所述步骤2)中，(1)沉积相类型数据库f_m
对碎屑岩储层沉积相类型进行概括和分类，其中包括：冲积扇-旱扇-扇根f_1、冲积扇-旱扇-扇中f_2、冲积扇-旱扇-扇缘f_3、冲积扇-湿扇-扇根f_4、冲积扇-湿扇-扇中f_5、冲积扇-湿扇-扇缘f_6、河流相-曲流河-河床亚相-河床滞留f_7、河流相-曲流河-河床亚相-边滩f_8、河流相-曲流河-堤岸亚相-天然堤f_9、河流相-曲流河-堤岸亚相-决口扇f_10、河流相-曲流河-河漫亚相-河漫滩f_11、河流相-曲流河-泛滥盆地-河漫湖泊f_12、河流相-曲流河-河漫沼泽f_13、河流相-辫状河-牛轭湖f_14、河流相-辫状河-河床滞留f_15、河流相-辫状河-心滩f_16、河道f_17、河流相-辫状河-泛滥平原f_18、湖泊相-断陷型-湖成三角洲f_19、湖泊相-断陷型-滨湖f_20、湖泊相-断陷型-浅湖f_21、湖泊相-断陷型-半深湖f_22、湖泊相-断陷型-深湖f_23、湖泊相-断陷型-湖湾f_24、湖泊相-坳陷型-湖成三角洲f_25、湖泊相-坳陷型-滨湖f_26、湖泊相-坳陷型-浅湖f_27、湖泊相-坳陷型-半深湖f_28、湖泊相-坳陷型-深湖f_29、湖泊相-坳陷型-湖湾f_30、湖泊相-前陆型-湖成三角洲f_31、湖泊相-前陆型-滨湖f_32、湖泊相-前陆型-浅湖f_33、湖泊相-前陆型-半深湖f_34、湖泊相-前陆型-深湖f_35、湖泊相-前陆型-湖湾f_36、三角洲相-辫状河三角洲-三角洲平原-分支f_37、三角洲相-辫状河三角洲-三角洲平原-河道f_38、三角洲相-辫状河三角洲-三角洲平原-陆上天然堤f_39、三角洲相-辫状河三角洲-三角洲平原-决口扇f_40、三角洲相-辫状河三角洲-三角洲平原-沼泽f_41、三角洲相-辫状河三角洲-三角洲平原-淡水湖泊f_42、三角洲相-辫状河三角洲-三角洲前缘-水下分支河道f_43、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤f_44、三角洲相-辫状河三角洲-三角洲前缘-支流间湾f_45、三角洲相-辫状河三角洲-三角洲前缘-分支河口砂坝f_46、三角洲相-辫状河三角洲-三角洲前缘-远砂坝f_47、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂f_48、三角洲相-辫状河三角洲-前三角洲f_49、三角洲相-扇三角洲-扇三角洲平原-泥石流f_50、三角洲相-扇三角洲-扇三角洲平原-河道充填f_51、三角洲相-扇三角洲-扇三角洲平原-漫滩f_52、三角洲相-扇三角洲-扇三角洲前缘-碎屑流f_53、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道f_54、三角洲相-扇三角洲-扇三角洲前缘-支流间湾f_55、三角洲相-扇三角洲-扇三角洲前缘-河口砂坝f_56、三角洲相-扇三角洲-扇三角洲前缘-远砂坝f_57、三角洲相-扇三角洲-前扇三角洲f_58；
(2)岩性数据库r_n
根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库，它包括细砾岩r_1、中砾岩r_2、粗砾岩r_3、巨砾岩r_4、石英砂岩r_5、长石质石英砂岩r_6、岩屑质石英砂岩r_7、长石岩屑质石英砂岩r_8、长石砂岩r_9、岩屑质长石砂岩r_10、岩屑长石砂岩r_11、岩屑砂岩r_12、长石质岩屑砂岩r_13、长石岩屑砂岩r_14、黏土r_15、泥岩r_16、页岩r_17；
(3)酸碱度数据库p_o
根据储层孔隙中流体酸碱性大小，建立酸碱度数据库；它包括强酸性p_1：ph＜3；酸性p_2：3≤ph＜5；弱酸性p_3：5≤ph＜7；弱碱性p_4：7≤ph＜9；碱性p_5：9≤ph＜11；强碱性p_6：11≤ph≤14；
(4)沉积旋回数据库c_p
对碎屑岩储层沉积旋回类型进行概括和分类，其中包括：正旋回c_1、反旋回c_2、复合旋回c_3。
(5)成岩阶段数据库s_q
成岩阶段数据库s_q包括：早成岩阶段a期ⅰa，即为：s_1；早成岩阶段b期ⅰb，即为：s_2；中成岩阶段a期ⅱa，即为：s_3；中成岩阶段b期ⅱb，即为s_4；晚成岩阶段ⅲ，即为：s_5；
(6)断层指数数据库n_r
根据研究区断层发育情况确定断层指数ω，建立断层指数数据库，其中包括：一级指数n_1：ω＜1；二级指数n_2：1≤ω≤3；三级指数n_3：3＜ω≤5；四级指数n_4：5＜ω。
再进一步地，所述步骤3)中，胶结作用数据库中胶结作用类型包括：
石英强胶结作用d1-1：其孔隙改变量为
石英中胶结作用d1-2：其孔隙改变量为
石英弱胶结作用d1-3：其孔隙改变量
长石强胶结作用d2-1：孔隙改变量
长石中胶结作用d2-2：其孔隙改变量
长石弱胶结作用d2-3：其孔隙改变量
白云石强胶结作用d3-1：其孔隙改变量
白云石中胶结作用d3-2：其孔隙改变量
白云石弱胶结作用d3-3：其孔隙改变量
方解石强胶结作用d4-1：其孔隙改变量
方解石中胶结作用d4-2：其孔隙改变量
方解石弱胶结作用d4-3：其孔隙改变量
黏土强胶结作用d5-1：其孔隙改变量
黏土中胶结作用d5-2：其孔隙改变量
黏土弱胶结作用d5-3：其孔隙改变量
本发明的有益效果在于：
本发明提供了一种碎屑岩储层胶结减孔量的定量预测方法，该方法基于碎屑岩储层沉积相、岩相、流体性质、沉积旋回、成岩阶段、断层发育情况的研究，建立不同成岩作用下“地质参数-成岩作用-孔隙度”预测模型，能对碎屑岩储层尤其是低孔、低渗致密储层孔隙度大小起到很好预测作用，降低勘探开发成本。
附图说明
图1为研究区埋藏史、温度史；
图2为研究区东营组二段沉积相；
图3为研究区东营组二段等值线图；
图4为研究区东营组二段沉积旋回图；
图5为研究区东营组二段流体酸碱性图；
图6为研究区东营组二段岩性图；
图7为研究区东营组二段成岩阶段图；
图8为研究区东营组二段断层指数图；
图9为研究区东营组二段胶结作用图；
图10为研究区东营组二段孔隙平面分布图。
具体实施方式
为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。
实施例1
胶结减孔量的定量预测方法，包括以下步骤：
1)收集研究区地质资料它包括：
①研究区目标层沉积相研究成果图；
②目标层界面构造等值线图；
③目标层上覆地层界面构造等值线图；
④研究区储层物性资料；
⑤研究区储层单井成岩作用数据；
⑥研究区成岩演化序列；
⑦研究区埋藏史、热演化史；
2)建立原始资料数据库
(1)沉积相类型数据库f_m
对碎屑岩储层沉积相类型进行概括和分类，其中包括：冲积扇-旱扇-扇根f_1、冲积扇-旱扇-扇中f_2、冲积扇-旱扇-扇缘f_3、冲积扇-湿扇-扇根f_4、冲积扇-湿扇-扇中f_5、冲积扇-湿扇-扇缘f_6、河流相-曲流河-河床亚相-河床滞留f_7、河流相-曲流河-河床亚相-边滩f_8、河流相-曲流河-堤岸亚相-天然堤f_9、河流相-曲流河-堤岸亚相-决口扇f_10、河流相-曲流河-河漫亚相-河漫滩f_11、河流相-曲流河-泛滥盆地-河漫湖泊f_12、河流相-曲流河-河漫沼泽f_13、河流相-辫状河-牛轭湖f_14、河流相-辫状河-河床滞留f_15、河流相-辫状河-心滩f_16、河道f_17、河流相-辫状河-泛滥平原f_18、湖泊相-断陷型-湖成三角洲f_19、湖泊相-断陷型-滨湖f_20、湖泊相-断陷型-浅湖f_21、湖泊相-断陷型-半深湖f_22、湖泊相-断陷型-深湖f_23、湖泊相-断陷型-湖湾f_24、湖泊相-坳陷型-湖成三角洲f_25、湖泊相-坳陷型-滨湖f_26、湖泊相-坳陷型-浅湖f_27、湖泊相-坳陷型-半深湖f_28、湖泊相-坳陷型-深湖f_29、湖泊相-坳陷型-湖湾f_30、湖泊相-前陆型-湖成三角洲f_31、湖泊相-前陆型-滨湖f_32、湖泊相-前陆型-浅湖f_33、湖泊相-前陆型-半深湖f_34、湖泊相-前陆型-深湖f_35、湖泊相-前陆型-湖湾f_36、三角洲相-辫状河三角洲-三角洲平原-分支f_37、三角洲相-辫状河三角洲-三角洲平原-河道f_38、三角洲相-辫状河三角洲-三角洲平原-陆上天然堤f_39、三角洲相-辫状河三角洲-三角洲平原-决口扇f_40、三角洲相-辫状河三角洲-三角洲平原-沼泽f_41、三角洲相-辫状河三角洲-三角洲平原-淡水湖泊f_42、三角洲相-辫状河三角洲-三角洲前缘-水下分支河道f_43、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤f_44、三角洲相-辫状河三角洲-三角洲前缘-支流间湾f_45、三角洲相-辫状河三角洲-三角洲前缘-分支河口砂坝f_46、三角洲相-辫状河三角洲-三角洲前缘-远砂坝f_47、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂f_48、三角洲相-辫状河三角洲-前三角洲f_49、三角洲相-扇三角洲-扇三角洲平原-泥石流f_50、三角洲相-扇三角洲-扇三角洲平原-河道充填f_51、三角洲相-扇三角洲-扇三角洲平原-漫滩f_52、三角洲相-扇三角洲-扇三角洲前缘-碎屑流f_53、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道f_54、三角洲相-扇三角洲-扇三角洲前缘-支流间湾f_55、三角洲相-扇三角洲-扇三角洲前缘-河口砂坝f_56、三角洲相-扇三角洲-扇三角洲前缘-远砂坝f_57、三角洲相-扇三角洲-前扇三角洲f_58；
(2)岩性数据库r_n
根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库，它包括细砾岩r_1、中砾岩r_2、粗砾岩r_3、巨砾岩r_4、石英砂岩r_5、长石质石英砂岩r_6、岩屑质石英砂岩r_7、长石岩屑质石英砂岩r_8、长石砂岩r_9、岩屑质长石砂岩r_10、岩屑长石砂岩r_11、岩屑砂岩r_12、长石质岩屑砂岩r_13、长石岩屑砂岩r_14、黏土r_15、泥岩r_16、页岩r_17；
(3)酸碱度数据库p_o
根据储层孔隙中流体酸碱性大小，建立酸碱度数据库；它包括强酸性p_1：ph＜3；酸性p_2：3≤ph＜5；弱酸性p_3：5≤ph＜7；弱碱性p_4：7≤ph＜9；碱性p_5：9≤ph＜11；强碱性p_6：11≤ph≤14；
(4)沉积旋回数据库c_p
对碎屑岩储层沉积旋回类型进行概括和分类，其中包括：正旋回c_1、反旋回c_2、复合旋回c_3。
(5)成岩阶段数据库s_q
成岩阶段数据库s_q包括：早成岩阶段a期ⅰa，即为：s_1；早成岩阶段b期ⅰb，即为：s_2；中成岩阶段a期ⅱa，即为：s_3；中成岩阶段b期ⅱb，即为s_4；晚成岩阶段ⅲ，即为：s_5；
(6)断层指数数据库n_r
根据研究区断层发育情况确定断层指数ω，建立断层指数数据库，其中包括：一级指数n_1：ω＜1；二级指数n_2：1≤ω≤3；三级指数n_3：3＜ω≤5；四级指数n_4：5＜ω；
3)建立胶结作用数据库d_c
根据岩石中不同矿物mi的孔隙改变量建立胶结作用数据库，胶结作用类型包括：
石英强胶结作用d1-1：其孔隙改变量为
石英中胶结作用d1-2：其孔隙改变量为
石英弱胶结作用d1-3：其孔隙改变量
长石强胶结作用d2-1：孔隙改变量
长石中胶结作用d2-2：其孔隙改变量
长石弱胶结作用d2-3：其孔隙改变量
白云石强胶结作用d3-1：其孔隙改变量
白云石中胶结作用d3-2：其孔隙改变量
白云石弱胶结作用d3-3：其孔隙改变量
方解石强胶结作用d4-1：其孔隙改变量
方解石中胶结作用d4-2：其孔隙改变量
方解石弱胶结作用d4-3：其孔隙改变量
黏土强胶结作用d5-1：其孔隙改变量
黏土中胶结作用d5-2：其孔隙改变量
黏土弱胶结作用d5-3：其孔隙改变量
4)建立成岩作用预测模型
(1)研究区储层网格化
对研究区储层网格化，即为将空间上不均匀分布的数据按一定方法(如滑动平均法、克里格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程；研究区储层的每个网格用wi(x，y)表示；
(2)确定网格属性a.根据研究区沉积相数据确定网格wi(x，y)的沉积相属性f_m；
b.根据研究区岩性数据确定网格wi(x，y)的岩性属性r_n；
c.根据研究区酸碱度数据确定网格wi(x，y)的酸碱度属性p_o；
d.根据研究区沉积旋回数据确定网格wi(x，y)的旋回属性c_p；
e.根据研究区成岩阶段数据确定网格wi(x，y)的成岩阶段属性s_q；
f.根据研究区断层指数数据确定网格wi(x，y)的断层指数属性n_r；
g.根据上述步骤a、b、c、d、e和f确定每个网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)，即：
wi(x，y)＝f_m+r_n+p_o+c_p+s_q+n_r；
(3)基于研究区岩石薄片资料、分析测试资料，确定研究区已知井hj所在网格wj(x，y)的成岩作用d_s_j；
(4)确定已知井hj所在网格wj(x，y)的综合属性gj(f_m_j，r_n_j，p_o_j，c_p_j，s_q_j，n_r_j)，建立已知井hj所在网格wj(x，y)成岩作用d_s与综合属性g(f_m_j，r_n_j，p_o_j，c_p_j，s_q_j，n_r_j)的对应关系，即为：
d_s_j＝f_m_j+r_n_j+p_o_j+c_p_j+s_q_j+n_r_j；
(5)任取一未知网格wi(x，y)，确定未知网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)，即为：
gi＝(f_m_i，r_n_i，p_o_i，c_p_i，s_q_i，n_r_i)；
(6)将网格wi(x，y)的成岩作用综合属性gi与井hj所在网格wj(x，y)的综合属性gj相比较，即为：
gi-gj＝(f_m_i，r_n_i，p_o_i，c_p_i，s_q_i，n_r_i)-(f_m_j+r_n_j+p_o_j+c_p_j+s_q_j+n_r_j)；
如果满足gi-gj＝0，则未知网格wi(x，y)的成岩作用d_s_i与井j具有相同的成岩作用d_s_j；
如果gi-gj≠0，则按不同属性优先级顺序，即：沉积相f_m一级、岩性r_n二级、沉积旋回c_p三级、酸碱度p_o四级、成岩阶段s_q五级、断层指数n_r六级对未知网格wi(x，y)的成岩作用进行判识，即为：
ⅰ：f_m_i-f_m_j＝0，r_n_i-r_n_j≠0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅱ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅲ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅳ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅴ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j＝0，n_r_i-n_r_j≠0；
未知网格wi(x，y)的成岩作用d_s_i满足条件i时，未知网格wi(x，y)与具有相同属性f_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅱ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅲ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅳ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、s_q_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅴ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、p_o_j、s_q_j的网格wj(x，y)的成岩作用相同；
(7)根据每个网格wi(x，y)的成岩作用确定研究区成岩作用横向分布；
5)建立胶结减孔量数学模型
建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
a.石英胶结模型d1
其中：a，b为常数，a＝1.98×10-22，b＝0.022℃(-1)；vque_cem为石英胶结作用改变孔隙体积百分量，％；m为石英的摩尔质量，g/mol；ρ为石英的密度，g/cm3；t为反应温度，℃；t＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；a0为自生石英开始形成时石英颗粒原始表面积；a0＝6fv/d；d为石英颗粒的粒径，mm；v为单位砂岩体积，mm3；f为单位体积砂岩中石英的体积百分含量；c为平均每年地层温度的平均改变量，℃/y；石英强胶结
b.长石胶结模型d2
其中：x埋藏深度(km)；vfel为长石胶结总量，％；长石强胶结
c.方解石胶结模型d3
其中：x埋藏深度(km)；vcal为方解石胶结总量，％；方解石强胶结
d.白云石胶结模型d4
其中：x埋藏深度(km)；vdol为白云石胶结总量，％；白云石强胶结
e.粘土胶结模型d5
其中：x埋藏深度(km)；vcla为胶结总量，％；粘土强胶结
6)胶结量模拟计算
(1)根据研究区埋藏史确定每个网格不同时刻t埋深h；
(2)根据研究区温度史确定每个网格不同时刻t温度t；
(3)根据步骤4)中成岩作用预测模型确定每个网格wi(x，y)胶结减孔量数学模型；
(4)根据步骤(1)中埋深h、步骤(2)中温度t、步骤(3)中胶结减孔量数学模型及其步骤4)中的成岩作用预测模型计算每个网格的胶结量。
实施例2
基于上述碎屑岩储层胶结减孔量的定量预测方法预测渤中凹陷东营组二段储层，具体步骤如下：
1、收集研究区地质资料
收集渤中凹陷东营组二段储层的地质资料，并按照原始资料数据库对各种类型数据进行分类。收集的地质资料包括以下内容：东营组二段储层沉积相图、东营组二段储层岩性图、东营组二段储层断裂分布图、东营组二段储层酸碱度图、东营组二段储层成岩阶段图、东营组二段储层沉积旋回图(如图所述)。
2、确定储层属性
(1)网格化工区
将渤中凹陷东营组二段碎屑岩储层进行网格化，每个网格用wi(x，y)表示。
(2)确定网格属性
a.根据研究区沉积相数据确定网格wi(x，y)的沉积相属性f_m；
b.根据研究区岩性数据确定网格wi(x，y)的岩性属性r_n；
c.根据研究区酸碱度数据确定网格wi(x，y)的酸碱度属性p_o；
d.根据研究区沉积旋回数据确定网格wi(x，y)的旋回属性c_p；
e.根据研究区成岩阶段数据确定网格wi(x，y)的成岩阶段属性s_q；
f.根据研究区断层指数数据确定网格wi(x，y)的断层指数属性n_r；
g.根据上述步骤a、b、c、d、e和f确定每个网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)，即为：
wi(x，y)＝f_m+r_n+p_o+c_p+s_q+n_r；
3、确定研究区成岩作用
(1)确定单井成岩作用
收集研究区岩石薄片分析测试数据，根据岩石颗粒接触关系、溶蚀类型、胶结类型、原生孔隙、次生孔隙以及胶结含量，结合成岩作用划分方案，确定不同井东营组二段的成岩作用类型，即已知井hi所在网格wi(x，y)的成岩作用所对应的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)。
(2)确定井间成岩作用
a.任取一未知网格wi(x，y)，确定未知网格wi(x，y)的综合属性g(f_m，r_n，p_o，c_p，s_q，n_r)；即为：
gi＝(f_m_i，r_n_i，p_o_i，c_p_i，s_q_i，n_r_i)
如网格w30＝(f_43，r_6，ⅰb，c_2，p_2，n1)；
如网格w40＝(f_43，r_10，ⅰb，c_3，p_2，n3)；
如网格w50＝(f_26，r_11，ⅱa，c_1，p_4，n2)；
如网格w60＝(f_45，r_6，ⅰb，c_3，p_2，n2)；b.将网格wi(x，y)的成岩作用综合属性gi与井hj(cfd16-1-2d、cfd16-1-1、hz4、cfd17-3-1、hz5、cfd18-1-2d、cfd18-2-2d、cfd18-2e-1、bz13-1-1、bz13-1-2、bz13-1-3、cfd23-1-1、cfd23-3-1、bz19-4-1、cfd24-1-1、cfd16-2-1、cfd16-3-1、cfd18-1n-1、bz19-2-1)所在网格wj(x，y)的综合属性gj相比较，即为：gi-gj＝(f_m_i+r_n_i+p_o_i+c_p_i+s_q_i+n_r_i)-(f_m_j+r_n_j+p_o_j+c_p_j+s_q_j+n_r_j)
如果满足gi-gj＝0，则未知网格wi(x，y)的成岩作用d_s_i与井j具有相同的成岩作用d_s_j；
如果gi-gj≠0，则按不同属性优先级(即为：沉积相一级、岩性二级、旋回三级、酸碱度四级、成岩阶段五级、断层指数六级)对未知网格wi(x，y)的成岩作用进行判识，即：
ⅰ：f_m_i-f_m_j＝0，r_n_i-r_n_j≠0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅱ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j≠0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅲ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j≠0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅳ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j≠0，n_r_i-n_r_j≠0；
ⅴ：f_m_i-f_m_j＝0，r_n_i-r_n_j＝0，c_p_i-c_p_j＝0，p_o_i-p_o_j＝0，s_q_i-s_q_j＝0，n_r_i-n_r_j≠0；
未知网格wi(x，y)的成岩作用d_s_i满足条件i时，未知网格wi(x，y)与具有相同属性f_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅱ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅲ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅳ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、s_q_j的网格wj(x，y)的成岩作用相同；
d_s_i满足条件ⅴ时，未知网格wi(x，y)与具有相同属性f_m_j、r_n_j、c_p_j、s_q_j、p_o_j的网格wj(x，y)的成岩作用相同；
如网格w30＝(f_43，r_6，ⅰb，c_2，p_2，n1)-(f_m_cfd17-3-1+r_n_cfd17-3-1+p_o_cfd17-3-1+c_p_cfd17-3-1+s_q_cfd17-3-1+n_r_cfd17-3-1)＝0，即网格w30与井cfd17-3-1具有相同的成岩作用，也即：w30＝d5-3_d4-1；
如网格w40＝(f_43，r_10，ⅰb，c_3，p_2，n3)-(f_cfd18-2e-1+r_n_cfd18-2e-1+p_o_cfd18-2e-1+c_p_cfd18-2e-1+s_q_cfd18-2e-1+n_r_cfd18-2e-1)＝0，即网格w40与井cfd18-2e-1具有相同的成岩作用，也即：w40＝d5-3_d4-3；
如网格w50＝(f_26，r_11，ⅱa，c_1，p_2，n1)网格属性满足条件ⅳ，即：f_m-f_m_bz-13-1-3＝0，r_n-r_n_bz-13-1-3＝0，c_p-c_p_bz-13-1-3＝0，s_q-s_q_bz-13-1-3＝0，p_o-p_o_bz-13-1-3≠0，n_r-n_r_bz-13-1-3≠0；则w80与井bz-13-1-3具有相同成岩作用，也即：w50＝d5-2_d1-3；
如网格w60＝(f_26，r_11，ⅱa，c_1，p_2，n1)网格属性满足条件ⅴ，即：f_m-f_m_cfd-16-2-1＝0，r_n-r_n_cfd-16-2-1＝0，c_p-c_p_cfd-16-2-1＝0，s_q-s_q_cfd-16-2-1＝0，p_o-p_o_cfd-16-2-1＝0，n_r-n_r_cfd-16-2-1≠0；则w60与井cfd-16-2-1具有相同成岩作用，也即：w60＝d5-2；
c.根据每个网格wi(x，y)的成岩作用确定研究区成岩作用横向分布。
4、确定研究区胶结减孔量数学模型
(1)石英胶结模型d1
其中：a，b为常数，a＝1.98×10-22，b＝0.022℃(-1)；vque_cem为石英胶结作用改变孔隙体积百分量，％；m为石英的摩尔质量，g/mol；ρ为石英的密度，g/cm3；t为反应温度，℃；t＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；a0为自生石英开始形成时石英颗粒原始表面积；a0＝6fv/d；d为石英颗粒的粒径，mm；v为单位砂岩体积，mm3；f为单位体积砂岩中石英的体积百分含量；c为平均每年地层温度的平均改变量，℃/y；石英强胶结
(2)方解石胶结模型d3
其中：x埋藏深度(km)；vcal为方解石胶结总量，％；方解石强胶结
(3)粘土胶结模型d5
其中：x埋藏深度(km)；vcla为胶结总量，％；粘土强胶结
(5)石英胶结模型主要参数选取
a＝1.98×10-22；
b＝0.022℃(-1)；
m＝60g/mol；
ρ＝0.25g/cm3；
c＝0.00000000000002℃/y；
cn＝0.0031℃/m；
dn＝18℃；
d＝0.02mm；
v＝1mm3；
f＝0.65；
t1＝90℃；
t2＝140℃；
6)胶结量模拟计算
(1)根据研究区构造深度数据，确定每个网格埋深h；
(2)根据研究区埋藏史、温度史，确定每个网格不同时刻t的埋深h(t)以及温度t；
(3)根据4)中成岩作用预测模型确定每个网格wi(x，y)胶结减孔量数学模型；
(4)根据步骤(1)埋深、步骤(2)温度、步骤(3)胶结作用及胶结数学模型计算每个网格的胶结量；
如网格w30，发生了黏土弱胶结、方解石强胶结，埋深h＝2300(m)，
网格w30总胶结＝vcla+vcal＝3.5+7.8＝11.3；
如网格w40，发生了黏土弱胶结、方解石弱胶结，埋深h＝2342，
网格w40总胶结＝vcla+vcal＝3.6+4.8＝8.4；
如网格w50，发生了黏土中胶结、石英弱胶结，埋深h＝2150，
网格w50总胶结＝vcla+vcal＝7.1+1.8＝8.9；如网格w60，发生了黏土中胶结，埋深h＝2400，
(5)根据每个网格的胶结量，确定东营组二段储层胶结减孔量。
其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。