1.一种碎屑岩储层埋藏过程中成岩相演化的恢复方法，其特征在于：包括以下步骤：

1)原始资料数据库建立

原始资料数据库包括以下数据库：(1)沉积相类型数据库Fi

沉积相决定储层原始物质分布，按碎屑岩沉积时沉积环境建立沉积相类型数据库Fi，(2)矿物类型数据库Mi

碎屑成岩过程中，以矿物类型分类并建立数据库，得到矿物类型数据库Mi；

(3)成岩相类型数据库Di

成岩相类型数据库Di包括压实相数据库D1、胶结相数据库D2和溶蚀相数据库D3(4)刚性颗粒强度数据库Vi根据岩石的刚性颗粒比大小将岩石的刚性强度分类并建立刚性颗粒强度数据库Vi；

2)建立基础模型

(1)沉积相与碎屑岩原始组构模型基于不同沉积环境下不同碎屑岩原始组成，建立“沉积相—碎屑岩原始物质”模型；

沉积相决定储层原始物质分布，建立“沉积相—原始物质”模型，为地质历史时期成岩相演化提供基础资料；

f(Sedimentary facies)＝f(Vi，Mi)其中：Sedimentary facies为沉积相，指沉积环境及其在该环境中形成的沉积岩特征的综合；Vi为刚性颗粒强度数据库，Mi为矿物类型数据库；

(2)压实相模型

a.根据地质参数R的强弱大小即刚性颗粒强度、胶结作用强度、溶蚀类型强度与压实之间关系，建立压实级别系数Gi模型：强溶蚀作用Y1，Gi＝‑3；中溶蚀作用Y2，Gi＝‑2；弱溶蚀作用Y3，Gi＝‑1；强胶结作用Y4，Gi＝5；中胶结作用Y5，Gi＝3；弱胶结作用Y6，Gi＝1；异常低压Y7，Gi＝1；正常压力Y8，Gi＝3；异常高压Y9，Gi＝5；超强刚性V1，Gi＝‑2；强刚性V2，Gi＝‑

1；中刚性V3，Gi＝‑0.5；弱刚性V4，Gi＝0；

b.根据储层在埋藏过程中，压实作用的强度与地质参数刚性颗粒、胶结作用、溶蚀类型有密切的关系，通过计算压实级别系数总和来判断碎屑岩储层压实作用的强弱；

压实级别系数总和计算公式：

其中，Gi为不同地质参数R的压实级别系数，EI为碎屑岩储层压实级别系数；

当10≥EI＞5，压实相类型为机械弱压实相D1‑3；

当5≥EI＞2，压实相类型为机械中压实相D1‑2；

当2≥EI，压实相类型为机械强压实相D1‑1；

(3)胶结相模型B

根据流体的条件和温度，确定胶结相模型B；

a.淡水‑半咸水条件下

当25＜T≤65℃，

当65＜T≤85℃，

当85＜T≤140℃，Di∈(D12‑2，D13‑2，D14‑2，D15‑2)；

当140＜T≤175℃，

当175＜T≤200℃，

b.酸性水条件下

当25＜T≤65℃，

当65＜T≤85℃，

当85＜T≤140℃，Di∈(D12‑2，D13‑2，D14‑2)；

当140＜T≤175℃，

当175＜T≤200℃，

c.碱性水条件下

当25＜T≤65℃，Di＝D11‑2；

当65＜T≤85℃，Di∈(D11‑2，D13‑3，D14‑3，D15‑1)；

当85＜T≤140℃，Di∈(D12‑3，D13‑2，D14‑2，D15‑2)；

当140＜T≤175℃，Di∈(D12‑2，D13‑2，D14‑2)；

当175＜T≤200℃，Di＝D12‑2；

(4)溶蚀相模型E

根据流体的条件和温度，确定溶蚀相模型E；

a.淡水‑半咸水条件下

当25＜T≤65℃，Di＝D5‑3；

当65＜T≤85℃，Di∈D3‑1，D4‑1，D5‑2，D7，D8，D11‑2，D15‑3)；

当85＜T≤140℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D7，D8，D9，D11‑2，D15‑2)；

当140＜T≤175℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D8，D9，D11‑1，D15‑2)；

当175＜T≤200℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D9，D11‑1，D15‑2)；

b.酸性水条件下

当25＜T≤65℃，Di＝D5‑3；

当65＜T≤85℃，Di∈(D5‑2，D7，D8，D11‑2)；

当85＜T≤140℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D7，D8，D9，D11‑2，D15‑2)；

当140＜T≤175℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D8，D9，D11‑1，D15‑2)；

当175＜T≤200℃，Di∈(D3‑2，D4‑2，D6，D9，D11‑1，D15‑2)；

c.碱性水条件下

当25＜T≤65℃，Di∈(D4‑2，D5‑3，D6，D15‑2)；

当65＜T≤85℃，Di∈(D4‑2，D5‑2，D6，D15‑3)；

当85＜T≤140℃，Di∈(D3‑3，D4‑2，D5‑2，D6，D11‑2，D15‑3)；

当140＜T≤175℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D11‑2)；

当175＜T≤200℃，Di∈(D3‑2，D4‑2，D5‑2，D6，D11‑2)；

其中，成岩相类型数据库Di包括压实相数据库D1、胶结相数据库D2和溶蚀相数据库D3a.压实相数据库

根据碎屑岩在成岩过程中压实强弱，即压实作用下孔隙改变量，建立压实相数据库D1，压实相类型包括：

强机械压实相D1‑1：其孔隙改变量为中机械压实相D1‑2：其孔隙改变量为弱机械压实相D1‑3：其孔隙改变量b.胶结相数据库D2

根据发生胶结作用的矿石类型和胶结强度共同建立胶结相数据库D2，胶结相类型包括：

石英强胶结相D2‑1：其孔隙改变量为石英中胶结相D2‑2：其孔隙改变量为石英弱胶结相D2‑3：其孔隙改变量长石强胶结相D3‑1：孔隙改变量长石中胶结相D3‑2：其孔隙改变量长石弱胶结相D3‑3：其孔隙改变量白云石强胶结相D4‑1：其孔隙改变量白云石中胶结相D4‑2：其孔隙改变量白云石弱胶结相D4‑3：其孔隙改变量方解石强胶结相D5‑1：其孔隙改变量方解石中胶结相D5‑2：其孔隙改变量方解石弱胶结相D5‑3：其孔隙改变量伊利石胶结相D6‑1：其孔隙改变量伊利石胶结相D6‑2：其孔隙改变量伊利石胶结相D6‑3：其孔隙改变量高岭石胶结相D7‑1：其孔隙改变量高岭石胶结相D7‑2：其孔隙改变量高岭石胶结相D7‑3：其孔隙改变量蒙脱石胶结相D8‑1：其孔隙改变量蒙脱石胶结相D8‑2：其孔隙改变量蒙脱石胶结相D8‑3：其孔隙改变量绿泥石胶结相D9‑1：其孔隙改变量绿泥石胶结相D9‑2：其孔隙改变量绿泥石胶结相D9‑3：其孔隙改变量 沸石强胶结相D10‑1：其孔隙改变量沸石中胶结相D10‑2：其孔隙改变量沸石弱胶结相D10‑3：其孔隙改变量c.溶蚀相数据库D3

根据发生溶蚀作用的矿石类型和溶蚀强度共同建立胶结相数据库D3，溶蚀相类型包括：

石英强溶蚀相D11‑1：其孔隙改变量石英中溶蚀相D11‑2、其孔隙改变量石英弱溶蚀相D11‑3、其孔隙改变量长石强溶蚀相D12‑1：其孔隙改变量长石中溶蚀相D12‑2：其孔隙改变量长石弱溶蚀相D12‑3：其孔隙改变量白云石强溶蚀相D13‑1：其孔隙改变量白云石中溶蚀相D13‑2：其孔隙改变量白云石弱溶蚀相D13‑3：其孔隙改变量方解石强溶蚀相D14‑1：其孔隙改变量方解石中溶蚀相D14‑2：其孔隙改变量方解石弱溶蚀相D14‑3：其孔隙改变量沸石强溶蚀相D15‑1：其孔隙改变量沸石中溶蚀相D15‑2：其孔隙改变量沸石弱溶蚀相D15‑3：其孔隙改变量(5)地层埋深模型

(a)确定演化期次

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

2……L2、L1；

(b)确定不同期次埋深

在(a)的基础下，计算目的层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阶段埋深；

(6)地层温度模型

地质体在埋藏的过程中，温度的大小表示为与深度的线性函数关系，通过该模型，计算目的层L在不同时期、不同深度、不同位置地层温度；

温度计算模型公式：

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

(7)压力模型

a.地层压力

地质体在埋藏的过程中，地层压力大小表示为深度的非线性函数关系；通过该模型，计算目的层L在不同时期、不同深度、不同位置的地层压力；

地层压力计算模型公式：

Pf＝a*Dep(Li)^2+b*Dep(Li)+c其中a、b、c为常数，通过井上实测地层压力数据进行数据拟合得到，Dep(Li)为目的层Li阶段埋深；Pf为地层压力b.静水压力

大小与液体的密度和液柱的高度有关，而与液体的形状和大小无关；

PH＝h*PW*g

其中，PH为静水压力，h为静水柱高度，PW为水的密度，g为重力加速度；

c.异常压力

当压力系数小于0.8时，为异常低压Y7；

当压力系数介于0.8～1.2时，为正常压力Y8；当压力系数大于1.2时，为异常高压Y9；

压力系数αp计算公式如下：

其中，Pf为地层压力，PH为静水压力；

3)确定成岩相

(1)储层网格化

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

(2)确定沉积相

根据研究区地质资料建立起每个网格的沉积相属性即为：W(x，y)＝f(Sedimentary facies)；

(3)确定原始组构

根据研究区已知井沉积相与原始组构关系，确定每个空间网格原始组构数据，即：f(Sedimentary facies)＝f(V1，Mi)，W(x，y)＝f(V1，Mi)其中：Sedimentary facies为沉积相，指沉积环境及其在该环境中形成的沉积岩特征的综合；Vi为刚性颗粒强度数据库，Mi为矿物类型数据库；

(4)确定演化期次

根据目的层上覆地层发育情况划分目标层在地质历史时期所经历的演化阶段，如果目的层L上覆地层有i层(，则目的层在地质历史时期演化阶段共有i个，按照时间演化的先后顺序依次计为Li、Li‑1、Li‑2……L2、L1；

(5)确定不同期次埋深

根据确定研究区储层的每个网格的演化期次，确定不同期次的地层埋深；埋深计算公式如下：

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阶段埋深；

(6)确定不同期次温度

将研究区目的层L的每个网格W(x，y)的参数代入温度计算模型公式中：T＝T0+c*(Dep(Ti)‑H0)计算得到目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置地层温度；

(7)确定不同期次压力

将研究区目的层L的每个网格W(x，y)的参数代入地层压力计算模型公式：Pf＝a*Dep(Li)^2+b*Dep(Li)+cPH＝h*PW*g

计算得到目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置异常压力，其中压力包括静水压力、地层压力、异常压力；

(8)计算溶蚀相

将目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置地层温度代入上述溶蚀相模型E中，确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置溶蚀相；

(9)计算胶结相

将目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置地层温度代入上述胶结相模型B中，确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置胶结相；

(10)计算压实相

a.根据步骤3)第(8)小步的溶蚀相确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的溶蚀相；

b.根据步骤3)第(9)小步的胶结相确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的胶结相；

c.确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的压力地质参数；

压力系数αp计算公式如下：

其中，Pf为地层压力，PH为静水压力；

d.确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的刚性地质参数；

刚性地质参数根据研究区“沉积相—原始物质”来确定；

f(Sedimentary facies)＝f(Vi，Mi)其中：Sedimentary facies为沉积相，指沉积环境及其在该环境中形成的沉积岩特征的综合；Vi为刚性颗粒强度数据库，Mi为矿物类型数据库；

e.根据a、b、c、d确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的压实相；

压实级别系数总和计算公式：

其中，Gi为不同地质参数R的压实级别系数，EI为碎屑岩储层压实级别系数；

当10≥EI＞5，Di＝机械弱压实相D1‑3；

当5≥EI＞2，Di＝机械中压实相D1‑2当2≥EI，Di＝机械强压实相D1‑1；

(11)确定成岩相

根据上述计算的压实相、胶结相和溶蚀相，确定目的层L的每个网格W(x，y)在不同时期、不同深度、不同位置的成岩相；

成岩相＝压实相+胶结相+溶蚀相。

2.根据权利要求1所述碎屑岩储层埋藏过程中成岩相演化的恢复方法，其特征在于：所述步骤1)中第(1)小步中，沉积相类型为：冲积扇‑旱扇‑扇根Ⅰ‑A‑a、冲积扇‑旱扇‑扇中Ⅰ‑A‑b、冲积扇‑旱扇‑扇缘Ⅰ‑A‑c、冲积扇‑湿扇‑扇根Ⅰ‑B‑a、冲积扇‑湿扇‑扇中Ⅰ‑B‑b、冲积扇‑湿扇‑扇缘Ⅰ‑B‑c、河流相‑曲流河‑河床亚相‑河床滞留Ⅱ‑A‑a‑1、河流相‑曲流河‑河床亚相‑边滩Ⅱ‑A‑a‑2、河流相‑曲流河‑堤岸亚相‑天然堤Ⅱ‑A‑b‑1、河流相‑曲流河‑堤岸亚相‑决口扇Ⅱ‑A‑b‑2、河流相‑曲流河‑河漫亚相‑河漫滩Ⅱ‑A‑c‑1、河流相‑曲流河‑泛滥盆地‑河漫湖泊Ⅱ‑A‑d‑1、河流相‑曲流河‑河漫沼泽Ⅱ‑A‑e、河流相‑辫状河‑牛轭湖Ⅱ‑B‑a、河流相‑辫状河‑河床滞留Ⅱ‑B‑b、河流相‑辫状河‑心滩、河道Ⅱ‑B‑c、河流相‑辫状河‑泛滥平原Ⅱ‑B‑d、湖泊相‑断陷型‑湖成三角洲Ⅲ‑A‑a、湖泊相‑断陷型‑滨湖Ⅲ‑A‑b、湖泊相‑断陷型‑浅湖Ⅲ‑A‑c、湖泊相‑断陷型‑半深湖Ⅲ‑A‑d、湖泊相‑断陷型‑深湖Ⅲ‑A‑e、湖泊相‑断陷型‑湖湾Ⅲ‑A‑f、湖泊相‑坳陷型‑湖成三角洲Ⅲ‑B‑a、湖泊相‑坳陷型‑滨湖Ⅲ‑B‑b、湖泊相‑坳陷型‑浅湖Ⅲ‑B‑c、湖泊相‑坳陷型‑半深湖Ⅲ‑B‑d、湖泊相‑坳陷型‑深湖Ⅲ‑B‑e、湖泊相‑坳陷型‑湖湾Ⅲ‑B‑f、湖泊相‑前陆型‑湖成三角洲Ⅲ‑C‑a、湖泊相‑前陆型‑滨湖Ⅲ‑C‑b、湖泊相‑前陆型‑浅湖Ⅲ‑C‑c、湖泊相‑前陆型‑半深湖Ⅲ‑C‑d、湖泊相‑前陆型‑深湖Ⅲ‑C‑e、湖泊相‑前陆型‑湖湾Ⅲ‑C‑f、三角洲相‑辫状河三角洲‑三角洲平原‑分支Ⅳ‑A‑a‑1、三角洲相‑辫状河三角洲‑三角洲平原‑河道Ⅳ‑A‑a‑2、三角洲相‑辫状河三角洲‑三角洲平原‑陆上天然堤Ⅳ‑A‑a‑3、三角洲相‑辫状河三角洲‑三角洲平原‑决口扇Ⅳ‑A‑a‑4、三角洲相‑辫状河三角洲‑三角洲平原‑沼泽Ⅳ‑A‑a‑5、三角洲相‑辫状河三角洲‑三角洲平原‑淡水湖泊Ⅳ‑A‑a‑6、三角洲相‑辫状河三角洲‑三角洲前缘‑水下分支河道Ⅳ‑A‑b‑1、三角洲相‑辫状河三角洲‑三角洲前缘‑水下天然堤Ⅳ‑A‑b‑

2、三角洲相‑辫状河三角洲‑三角洲前缘‑支流间湾Ⅳ‑A‑b‑3、三角洲相‑辫状河三角洲‑三角洲前缘‑分支河口砂坝Ⅳ‑A‑b‑4、三角洲相‑辫状河三角洲‑三角洲前缘‑远砂坝Ⅳ‑A‑b‑

5、三角洲相‑辫状河三角洲‑三角洲前缘‑前缘席状砂Ⅳ‑A‑b‑6、三角洲相‑辫状河三角洲‑前三角洲Ⅳ‑A‑c、三角洲相‑扇三角洲‑扇三角洲平原‑泥石流Ⅳ‑B‑a‑1、三角洲相‑扇三角洲‑扇三角洲平原‑河道充填Ⅳ‑B‑a‑2、三角洲相‑扇三角洲‑扇三角洲平原‑漫滩Ⅳ‑B‑a‑3、三角洲相‑扇三角洲‑扇三角洲前缘‑碎屑流Ⅳ‑B‑b‑1、三角洲相‑扇三角洲‑扇三角洲前缘‑水下分流河道Ⅳ‑B‑b‑2、三角洲相‑扇三角洲‑扇三角洲前缘‑支流间湾Ⅳ‑B‑b‑3、三角洲相‑扇三角洲‑扇三角洲前缘‑河口砂坝Ⅳ‑B‑b‑4、三角洲相‑扇三角洲‑扇三角洲前缘‑远砂坝Ⅳ‑B‑b‑5、三角洲相‑扇三角洲‑前扇三角洲Ⅳ‑B‑c。

3.根据权利要求1所述碎屑岩储层埋藏过程中成岩相演化的恢复方法，其特征在于：所述步骤1)中第(2)小步中，矿物类型包括：石英M1、长石M2、白云石M3、方解石M4、伊利石M5、高岭石M6、蒙脱石M7、绿泥石M8、沸石M9、白云母M10、黑云母M11。

4.根据权利要求1所述碎屑岩储层埋藏过程中成岩相演化的恢复方法，其特征在于：所述步骤1)中第(4)小步中，刚性颗粒强度数据库Vi包括超强刚性颗粒V1，刚性颗粒比例＞

75；强刚性颗粒V2，刚性颗粒比例＞60且≤75；中刚性颗粒V3刚性颗粒比例＞45且≤60；弱刚性颗粒V4，刚性颗粒比例≤45。