1.一种碎屑岩储层胶结减孔量的定量预测方法，其特征在于：包括以下步骤：

1)收集研究区地质资料

2)建立原始资料数据库，它由以下组成(1)沉积相类型数据库

(2)岩性数据库

(3)酸碱度数据库

(4)沉积旋回数据库

(5)成岩阶段数据库

(6)断层指数数据库

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_j与综合属性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

‑22 (‑1)

其中：a，b为常数，a＝1.98×10 ，b＝0.022℃ ；VQue_Cem为石英胶结作用改变孔隙体3

积百分量，％；M为石英的摩尔质量，g/mol；ρ为石英的密度，g/cm ；T为反应温度，℃；T＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；A0为自生石英开始形成时石英颗粒原始表面积；A0＝6fV/D；D为石英颗粒的粒径，mm；V为单位砂岩体3

积，mm ；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所述碎屑岩储层胶结减孔量的定量预测方法，其特征在于：所述步骤

1)中，收集研究区地质资料由以下组成①研究区目标层沉积相研究成果图；

②目标层界面构造等值线图；

③目标层上覆地层界面构造等值线图；

④研究区储层物性资料；

⑤研究区储层单井成岩作用数据；

⑥研究区成岩演化序列；

⑦研究区埋藏史、热演化史。

3.根据权利要求1所述碎屑岩储层胶结减孔量的定量预测方法，其特征在于：所述步骤

2)中，

(1)沉积相类型数据库

对碎屑岩储层沉积相类型进行概括和分类，其中由以下组成：冲积扇‑旱扇‑扇根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_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_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_1、反旋回C_

2、复合旋回C_3；

(5)成岩阶段数据库

成岩阶段数据库由以下组成：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_

2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；

(6)断层指数数据库

根据研究区断层发育情况确定断层指数Ω，建立断层指数数据库，其中由以下组成：一级指数N_1：Ω＜1；二级指数N_2：1≤Ω≤3；三级指数N_3：3＜Ω≤5；四级指数N_4：5＜Ω。

4.根据权利要求1所述碎屑岩储层胶结减孔量的定量预测方法，其特征在于：所述步骤

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：其孔隙改变量