[发明专利]压实减孔量的定量预测方法

权利要求书

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；

F(Sedimentary facies)为沉积相属性的函数，用以表达不同空间位置沉积相属性；

(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包括Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ、Ⅷ、Ⅸ、Ⅹ总共十种模型；

Ⅰ:f₁＝1.9467*Ln(Dep)-0.8689；

Ⅱ:f₂＝2.9200*Ln(Dep)-1.3034；

Ⅲ:f₃＝4.3800*Ln(Dep)-1.9551；

Ⅳ:f₄＝7.0080*Ln(Dep)-3.1282；

Ⅴ:f₅＝8.7600*Ln(Dep)-3.9103；

Ⅵ:f₆＝11.6800*Ln(Dep)-5.2137；

Ⅶ:f₇＝14.0160*Ln(Dep)-6.2564；

Ⅷ:f₈＝17.5200*Ln(Dep)-7.8205；

Ⅸ:f₉＝20.6118*Ln(Dep)-9.2006；

Ⅹ:f₁₀＝23.3600*Ln(Dep)-10.4273；

其中，f_i为压实率，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)；其中，Wm(xm,ym,zm)为每个网格化的储层的空间函数，xm是横坐标，ym是纵坐标，zm是垂直坐标；

(2)确定每个网格化的储层沉积相

根据上述建立的沉积相类型数据库，分析每个网格化的储层的沉积相类型；

W_m(x_m，y_m，z_m)＝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)确定演化期次

根据研究区埋藏史数据，确定每个网格化的储层即为网格化的目标储层在地质历史时期演化期次，演化期次以t j表示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)确定网格化的目标储层不同演化阶段压实模型及孔隙度大小

g_j＝f_i；

其中，g_j为演化阶段j的压实率改变量，f_i为相应的压实率模型，为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.根据研究区地质资料确定每个网格的演化阶段t_j，j＝1，2，3……；

c.根据原始组构确定每个网格初始压实模型IM；

d.根据成岩作用以及原始组构确定每个网格不同演化阶段t_j地质指数；

e.在步骤c和步骤d的基础上确定每个网格不同演化阶段t_j压实模型g_j；

f.根据初始孔隙度初始压实模型IM，不同演化阶段压实模型g_j计算不同演化阶段t_j时机械压实作用下储层孔隙减少量。