1.一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，其特征在于：包括以下步骤：

1）收集研究区地质资料

2）建立原始资料数据库

（1）沉积相类型数据库

对碎屑岩储层沉积相类型进行概括和分类建立沉积相类型数据库；

（2）岩性数据库

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

（3）流体数据库

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

（4）成岩阶段数据库

（5）建立溶蚀作用数据库根据岩石中不同矿物的溶蚀率建立溶蚀作用数据库；

（6）建立溶蚀级别数学模型建立起不同类型溶蚀物的溶蚀级别数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；

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）研究区储层网格化

对研究区储层网格化，研究区储层的每个网格用 （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_j与综合属性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，VCem=∅0*RCem，其中∅0为原始孔隙；

（5）根据每个网格Wi（X，Y）不同演化期次的溶蚀量确定目标储层溶蚀量，即溶蚀增孔量。

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、酸性水环境P_2、碱性水环境P_3；

（4）成岩阶段数据库由以下组成：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；

（5）溶蚀作用数据库中溶蚀作用类型由以下组成：石英强溶蚀作用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%；

（6）溶蚀级别数学模型由以下几类组成：B_1=Ⅰ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/10；

B_2=Ⅱ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/8；

B_3=Ⅲ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6] /6；

B_4=Ⅳ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6] /5；

B_5=Ⅴ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6] /3；

B_6=Ⅵ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/1.6；

B_7=Ⅶ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/1.5；

B_8=Ⅷ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/1.4；

B_9=Ⅸ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/1.2；

B_10=Ⅹ：Y =[‑4.859e‑06* (X)^4+0.001028* (X)^3‑0.0745* (X)^2 +2.862* (X)+

10.6]/1；

其中，X为埋藏时间，单位：Ma，Y为溶蚀率，单位：%。