[發明專利]地質歷史時期砂巖儲層滲透率演化恢復方法

有效

• 申請號: 201210235345.6
• 申請日: 2012-07-10
• 公開（公告）號: CN102778421A
• 公開（公告）日: 2012-11-14
• 發明（設計）人: 王艷忠;操應長;葸克來;張善文;王永詩;劉惠民;陳林;馬奔奔
• 申請（專利權）人: 中國石油大學(華東)
• 主分類號: G01N15/08
• 分類號: G01N15/08
• 代理公司: 暫無信息
• 代理人: 暫無信息

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/Φ的關系，結合孔隙結構恢復與分類結果，求取地質歷史時期砂巖儲層滲透率；

第四步；建立地質歷史時期砂巖儲層滲透率演化曲線。