[发明专利]地质历史时期砂岩储层渗透率演化恢复方法

权利要求书

1.一种地质历史时期砂岩储层渗透率演化恢复方法，具体步骤如下：

(1)地质历史时期砂岩储层孔隙度演化恢复；

(2)地质历史时期砂岩储层孔隙结构演化恢复

第一步：利用反演回剥法，在成岩作用演化序列的约束下，将胶结作用及溶蚀作用成岩事件进行反演回剥；

第二步：对压实损失的孔隙面积进行等比例放大，具体恢复方法如下：

①参量设置：

待恢复状态：

设视域长为L，宽为D

岩石体积：V_2z

孔隙体积：V_2P

视域面积：S_2z(未知)

孔隙面积：S_2P(未知)

面孔率：A₂＝S_2P/S_2z*100％(可以由Φ于A的关系求取，A₂＝F(Φ₂))

孔隙度：Φ₂＝V_2P/V_2z*100％(已知)

已知状态：

设视域长为L+ΔL，宽为D+ΔD

岩石体积：V_1z

孔隙体积：V_1P

视域面积：S_1z(已知)

孔隙面积：S_1P(已知)

面孔率：A₁＝S_1P/S_1z*100％(已知)

孔隙度：Φ1＝V_1P/V_1z*100％(已知)

设由待恢复状态到已知状态压实减小的面积为：S_y，取正直

设由待恢复状态到已知状态胶结减小的面积为：S_j，取正直

设由待恢复状态到已知状态溶解增加的面积为：S_r，取正直

②推导过程：

S_2p-S_y-S_j+S_r＝S_1p    公式1

A₂＝S_2P/S_2z          公式2

S_2z-S_y＝S_1z          公式3

把公式3带入公式2可得，

A₂＝S_2P/S_2z＝S_2P/(S_1z+S_y)

S_2p＝A₂*(S_1z+S_y)            公式4

把公式4带入公式1：

A₂*(S_1z+S_y)-S_y-S_j+S_r＝S_1p

A₂*S_1z+A2*S_y-S_y-S_j+S_r＝S_1p

S_y(1-A₂)＝A₂*S_1z+S_r-S_j-S_1p

S_y＝(A₂*S_1z-S_1p+S_r-S_j)/(1-A₂)    公式5

S_y＝(A₂*S_1z-A₁*S_1z+S_r-S_j)/(1-A2)

S_y＝(S_1z(A₂-A₁)+S_r-S_j)/(1-A2)    公式5’

S_y＝(A₂*S_1z-S_1p-S_j+S_r)/(1-A₂)

设由待恢复状态到已知状态压实作用使总视域长宽成等比例缩小，即

L/D＝(L+ΔL)/(D+ΔD)    公式6

S_1z+S_y＝(L+ΔL)*(D+ΔD)    公式7

将公式6带入公式7可得，

S_1z+S_y＝(L/D)*(D+ΔD)²

D(S_1z+S_y)/L＝(D+ΔD)²

ΔD＝[D(S_1z+S_y)/L]^1/2-D    公式8

ΔL＝L*(D+ΔD)/D-L

ΔL＝L*ΔD/D    公式9

求得ΔD、ΔL后，选择中心颗粒不动，从两端开始按比例移动颗粒，假设某一颗粒横向到纵向中心线的距离为d，则需向边移动的距离Δd＝2ΔL d/(ΔL+L)；

第三步：建立地质历史时期砂岩储层孔隙结构演化图版，并进行各时期孔隙结构分类；

(3)地质历史时期砂岩储层渗透率演化恢复

第一步：利用压汞资料，根据排驱压力(Pd)、毛管压力中值(P50)等参数，建立每一类孔隙结构的K与K/Φ关系，其中：K为渗透率，Φ为孔隙度；

第二步：根据每一类孔隙结构K与K/Φ的关系，结合孔隙结构恢复与分类结果，求取地质历史时期砂岩储层渗透率；

第四步；建立地质历史时期砂岩储层渗透率演化曲线。