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

1)收集研究区地质资料

2)建立原始资料数据库，它包括：

(1)沉积相类型数据库F_m

对碎屑岩储层沉积相类型进行概括和分类，

(2)岩性数据库R_n

根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库；

(3)流体数据库P_o

根据储层孔隙中流体性质，建立流体数据库；

(4)成岩阶段数据库S_q

(5)建立胶结作用数据库D_c

根据岩石中不同矿物Mi的胶结率建立胶结作用数据库；

(6)建立胶结级别数学模型B_k

建立起不同类型胶结物Mi的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；

3)建立成岩阶段预测模型

(1)确定演化期次

根据目的层上覆地层发育情况划分目标层在地质历史时期所经历的演化阶段，目的层L上覆地层有i层，即从上到下依次标记为：L1、L2、L3……Li-1、Li，则目的层在地质历史时期演化阶段共有i个，按照时间演化，目的层L演化期次先后顺序依次计为Li、Li-1、Li-

2……L2、L1；

(2)确定不同期次埋深

在(1)的基础下，计算目的层L不同演化期次的埋深，其计算公式如下：Dep(Li)＝H(L)-H(Li)；

Dep(Li-1)＝H(L)-H(Li-1)；

Dep(Li-2)＝H(L)-H(Li-2)；

……

Dep(L2)＝H(L)-H2；

Dep(L1)＝H(L)-H1；

其中：H(L)为目的层L顶界面，H(Li)为上覆层Li顶界面，Dep(Li)为目的层Li阶段埋深；

(3)确定不同期次地层温度

地质体在埋藏的过程中，温度的大小表示为与深度的线性函数关系，通过该模型，计算目的层L在不同时期、不同深度、不同位置地层温度；其温度计算模型公式：T＝T0+c*(D(Ti)-H0)

其中T0为常温带温度，c为常数，Dep(Ti)为目的层Ti阶段埋深，H0为恒温带埋深，为常数，T为目标层古地温；

(4)确定不同期次地层成岩阶段

当T∈[20～65)，DS为早成岩阶段早期，即为ⅠA；

当T∈[65～85)，DS为早成岩阶段晚期，即为ⅠB；

当T∈[85～140)，DS为中成岩阶段早期，即ⅡA；

当T∈[140～175)，DS为中成岩阶段晚期，即为ⅡB；

当T∈[175～200)，DS为晚成岩阶段，即为Ⅲ；

4)建立不同成岩阶段胶结作用预测模型

(1)研究区储层网格化

将研究区储层进行网格化，研究区储层的每个网格用Wi(X，Y)表示；

(2)确定网格属性

a.根据步骤3)确定网格Wi(X，Y)不同演化期次成岩阶段S_q；

b.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；

c.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；

d.根据研究区流体数据确定网格Wi(X，Y)的流体属性P_o；

e.根据上述步骤a、b、c和d确定每个网格不同演化期次Wi(X，Y)的综合属性G(S_q，F_m，R_n，P_o)，即：Wi(X，Y)＝S_q+F_m+R_n+P_o；

(3)基于研究区地质资料，确定研究区已知井Hj所在网格Wj(X，Y)不同演化期次的成岩阶段S_q_j以及胶结作用D_s_j；

(4)确定已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段的综合属性Gj(F_m_j，R_n_j，P_o_j)，建立已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段胶结作用D_s与综合属性G(F_m_j，R_n_j，P_o_j)的对应关系，即为：Wj(S_q_j，D_s_j)＝F_m_j+R_n_j+P_o_j；

(5)任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)成岩阶段S_q时的综合属性G(F_m，R_n，P_o)，即为：G i＝(F_m_i，R_n_i，P_o_i)；

(6)将网格Wi(X，Y)成岩阶段S_q时的胶结作用综合属性G i与井Hj所在网格Wj(X，Y)成岩阶段S_q的综合属性G j相比较，即为：Gi-Gj＝(F_m_i，R_n_i，P_o_i)-(F_m_j+R_n_j+P_o_j)；

如果满足G i-G j＝0，则未知网格Wi(X，Y)的胶结作用D_s_i与井j具有相同的胶结作用D_s_j；

如果G i-G j≠0，则按不同属性优先级顺序，即为：沉积相F_m一级、岩性R_n二级、流体性质P_o三级对未知网格Wi(X，Y)的胶结作用进行判识，即：ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，P_o_i-P_o_j≠0；

ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，P_o_i-P_o_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)的胶结作用相同；

(7)根据每个网格Wi(X，Y)不同演化期次的胶结作用确定研究区胶结作用演化；

5)计算不同演化期次胶结量

(1)根据成岩序列数据，确定已知井所在网格Wj(X，Y)的成岩阶段S_q所对应的胶结作用下的胶结级别B_k_j，即Wj(S_q_j，D_s_j)＝B_k_j；

(2)将Wi(X，Y)的Wi(S_q_j，D_s_j)与Wj(S_q_j，D_s_j)比对，满足Wi(S_q_j，D_s_j)＝Wj(S_q_j，D_s_j)，则Wi(X，Y)与Wj(S_q_j，D_s_j)具有相同的胶结级别B_k，即B_k_i＝B_k_j；如果出现同一未知网格Wi(X，Y)与多个已知网格Wj(S_q_j，D_s_j)相同，则B_k_i＝1/n∑B_k_j，n为相同网格Wj(S_q_j，D_s_j)的网格数；

(3)根据网格Wi(X，Y)不同演化期次的胶结级别B_k_i计算每个网格Wi(X，Y)不同演化期次的胶结率PCem；

(4)根据每个网格Wi(X，Y)的胶结率RCem计算胶结量VCem， 其中 为原始孔隙；

(5)根据每个网格Wi(X，Y)不同演化期次的胶结量确定目标储层胶结量，即胶结减孔量。

2.根据权利要求1所述基于成岩过程的碎屑岩储层胶结减孔量的定量预测方法，其特征在于：所述步骤1)中，收集研究区地质资料，它包括：①研究区东营组二段胶结作用表；

②研究区平原组构造等值线图；

③研究区明化镇组构造等值线图；

④研究区馆陶组构造等值线图；

⑤研究区东营组二段构造等值线图；

⑥研究区东营组二段沉积相图；

⑦研究区东营组二段岩性分布图；

⑧研究区馆陶组沉积时东营组二段流体分布图；

⑨研究区明化镇组沉积时东营组二段流体分布图；

⑩研究区平原组沉积时东营组二段流体分布图。

3.根据权利要求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、酸性水环境P_2、碱性水环境P_3；

(4)成岩阶段数据库S_q包括：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_

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

(5)建立胶结作用数据库D_c

根据岩石中不同矿物Mi的胶结率建立胶结作用数据库，胶结作用类型包括：石英强胶结作用D1-1：其胶结率RCem＞10％；

石英中胶结作用D1-2：其胶结率5％＜RCem≤10％；

石英弱胶结作用D1-3：其胶结率RCem≤5％；

长石强胶结作用D2-1：其胶结率RCem＞10％；

长石中胶结作用D2-2：其胶结率5％＜RCem≤10％；

长石弱胶结作用D2-3：其胶结率RCem≤5％；

方解石强胶结作用D3-1：其胶结率RCem＞10％；

方解石中胶结作用D3-2：其胶结率5％＜RCem≤10％；

方解石弱胶结作用D3-3：其胶结率RCem≤5％；

白云石强胶结作用D4-1：其胶结率RCem＞10％；

白云石中胶结作用D4-2：其胶结率5％＜RCem≤10％；

白云石弱胶结作用D4-3：其胶结率RCem≤5％；

黏土强胶结作用D5-1：其胶结率RCem＞10％；

黏土中胶结作用D5-2：其胶结率RCem＞10％；

黏土弱胶结作用D5-3：其胶结率RCem≤5％；

(6)建立胶结级别数学模型B_k

建立起不同类型胶结物Mi的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量。

B_1＝Ⅰ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/10；

B_2＝Ⅱ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/8；

B_3＝Ⅲ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/6；

B_4＝Ⅳ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/5；

B_5＝Ⅴ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/3；

B_6＝Ⅵ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/1.6；

B_7＝Ⅶ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/1.5；

B_8＝Ⅷ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/1.4；

B_9＝Ⅸ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/1.2；

B_10＝Ⅹ：Y＝[7.515e-05*(X)^3-0.0158*(X)^2+1.686*(X)+14.38]/1；

其中，X埋藏时间，Ma，Y胶结率，％。