1.一种压实减孔量的定量预测方法，其特征在于：包括以下步骤：
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(Sedimentary facies)＝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为地质参数，其大小根据压实级别G来定，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(Sedimentary facies)
(3)确定每个网格化的储层的原始组构
根据研究区中每个网格化的储层的沉积相数据确定每个空间网格原始组构数据，并根据研究区已知井“沉积相-原始组构数据”对未知区域“沉积相-原始组构”数据F(Sedimentary facies)＝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.根据权利要求1所述压实减孔量的定量预测方法，其特征在于：所述步骤1)中，收集研究区的地质资料包括：
①研究区目标层沉积相研究成果图；
②目标层界面构造等值线图；
③目标层上覆地层界面构造等值线图；
④研究区储层物性资料；
⑤研究区储层成岩相研究成果图；
⑥研究区成岩演化序列；
⑦研究区埋藏史、热演化史。
3.根据权利要求1所述压实减孔量的定量预测方法，其特征在于：所述步骤2)的第(2)小步中，
a.建立粒度D_a数据库
按照碎屑岩粒度大小，建立粒度数据库；
表碎屑岩粒度D_a类型及其分类
粒级划分 数据库代码 颗粒直径(mm) 压实级别(G_1)
巨砾 D_1 ＞1000 1
粗砾 D_2 ＞100～1000 2
中砾 D_3 ＞10～100 3
细砾 D_4 ＞1～10 4
粗砂 D_5 ＞0.5～1 5
中砂 D_6 ＞0.25～0.5 6
细砂 D_7 ＞0.1～0.25 7
粗粉砂 D_8 ＞0.05～0.1 8
细粉砂 D_9 ＞0.01～0.05 9
粘土 D_10 ≤0.01 10
b.建立分选性S_b数据库
按照碎屑岩颗粒大小的标准偏差即均一性，建立碎屑岩分选性数据库；
表碎屑岩分选性S_b类型及其分类
c.建立磨圆度R_c数据库
按照碎屑岩磨圆情况，建立磨圆度数据库；
表碎屑岩磨圆度R_c类型及其分类
磨圆度 数据库代码 压实级别(G_3)
棱角状 R_1 1
次棱角状 R_2 2
次圆状 R_3 3
圆状 R_4 4
d.建立刚性颗粒V_d数据库
按照碎屑岩中刚性颗粒含量，建立刚性颗粒数据库；
表碎屑岩刚性组分V_d及其分类
刚性颗粒 数据库代码 压实级别(G_4)
＞75 V_1 1
＞60～75 V_2 2
＞45～60 V_3 3
＞30～45 V_4 4
≤30 V_5 5
e.建立矿物组分M_e数据库
按照碎屑岩中矿物组分，建立矿物组分数据库；
表碎屑岩矿物组分M_e及其分类
矿物组分 数据库代码 压实级别(G_5)
石英、长石、方解石、白云石、云母、黏土矿物 M_1 1
石英、方解石、白云石 M_2 2
石英、方解石、白云石、黏土矿物 M_3 3
石英、长石 M_4 4
长石、方解石、白云石 M_5 5
石英、方解石 M_6 6
f.建立矿物比分O_f数据库
按照碎屑岩中长石含量建立矿物比分数据库；
表碎屑岩矿物比分Mn及其分类
长石百分比(％) 数据库代码 压实级别(G_6)
＞75 O_1 1
＞60～75 O_2 2
＞45～60 O_3 3
＞30～45 O_4 4
≤30 O_5 5
g.建立胶结类型C_g数据库
按照碎屑岩胶结类型，建立胶结类型数据库；
表胶结类型C_g数据库
胶结类型 数据库代码 压实级别(G_7)
钙质 C_1 1
硅质 C_2 2
铁质 C_3 3
黏土 C_4 4
h.建立胶结强度Q_h数据库
按照碎屑岩的胶结强度，建立胶结强度Q_h数据库；
表胶结类型Q_h数据库
i.建立胶结时间T_i数据库
按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；
表胶结时间T_i数据库
胶结时间 数据库代码 压实级别(G_9)
ⅠA T_1 5
ⅠB T_2 4
ⅡA T_3 3
ⅡB T_4 2
Ⅲ T_5 1
j.建立溶蚀时间F_j数据库
按照碎屑岩成岩阶段将碎屑岩胶结时间进行分类，建立储层演化过程中胶结时间数据库；
表溶蚀时间F_j数据库
溶蚀时间 数据库代码 压实级别(G_10)
ⅠA F_1 1
ⅠB F_2 2
ⅡA F_3 3
ⅡB F_4 4
Ⅲ F_5 5
k.建立泥质含量N_k数据库
按照碎屑岩中的泥质含量，建立泥质含量数据库；
表11泥质含量N_k分类
泥质含量(％) 数据库代码 压实级别(G_11)
＜2 N_1 1
≥2～5 N_2 2
＞5～8 N_3 3
＞8 N_4 4
l.建立岩屑含量X_l数据库
按照碎屑岩中岩屑含量，建立岩屑含量数据库；
表岩屑含量X_l分类
岩屑含量(％) 数据库代码 压实级别(G_12)
≤2 X_1 1
＞2～6 X_2 2
＞6～10 X_3 3
＞10 X_4 4
m.建立沉积速率(E_m)数据库
按照碎屑岩的沉积速率将碎屑岩储层进行分类，建立沉积速率数据库；
表沉积速率E_n分类
沉积速率(m/Ma) 数据库代码 压实级别(G_13)
＜100 E_1 1
≥100～500 E_2 2
＞500～1000 E_3 3
＞1000 E_4 4
n.超压(P_n)数据库
按照压力系数大小将碎屑岩储层演化过程中压力大小进行分类，建立超压数据库；
压力系数αp＝Pp/Ph
其中，Pp为实测地层压力，Ph为静水压力；
表超压P_n分类