[發明專利]裂縫滲透率的計算方法

有效

• 申請號: 201610563073.0
• 申請日: 2016-07-15
• 公開（公告）號: CN107622139B
• 公開（公告）日: 2020-08-07
• 發明（設計）人: 肖承文;江同文;王謙;肖香嬌;周波;信毅;韓闖;袁龍;孫雄偉
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;E21B49/00
• 代理公司: 北京同立鈞成知識產權代理有限公司 11205
• 代理人: 劉丹;黃健

1.一種裂縫滲透率的計算方法，其特征在于，包括：

根據試井滲透率、基質滲透率以及裂縫孔隙度得到裂縫平均滲透率，其中，所述試井滲透率根據無阻流量獲得；

建立所述裂縫平均滲透率、裂縫平均孔隙度、平均流體移動指數之間的第一計算模型；

根據所述試井滲透率對所述第一計算模型進行標定，建立裂縫滲透率、裂縫孔隙度、流體移動指數間的第二計算模型；

根據所述第二計算模型得到所述裂縫滲透率；

其中，所述根據試井滲透率、基質滲透率以及裂縫孔隙度得到裂縫平均滲透率，包括：根據式1得到所述裂縫平均滲透率：

其中，K_s為試井滲透率，md；K_d(m)對應深度點處的基質滲透率，md；K_l(m)為對應深度點處的裂縫滲透率，md；FVPA(m)為對應深度點處的裂縫孔隙度，％；dep1為試井段頂深，m；dep2為試井段底深，m；dep3為裂縫段頂深，m；dep4為裂縫段底深，m；

所述試井滲透率根據無阻流量獲得，包括：根據式2得到所述試井滲透率：

其中，Q為ΔP壓差下測試段的無阻流量，萬方/天；K_s為所述試井滲透率，md；H為測試段厚度，m；A_m為單位井段的滲流面積，m²；μ為流體粘度，mPa·s；L為測試段供氣半徑，m；ΔP為測試段地層壓力，MPa；

所述建立裂縫平均滲透率、裂縫平均孔隙度、平均流體移動指數之間的第一計算模型，包括：建立所述裂縫平均滲透率與裂縫平均孔隙度的第一關系模型，建立所述裂縫平均滲透率與平均流體移動指數的第二關系模型；根據所述第一關系模型和第二關系模型，利用多元回歸法，建立所述第一計算模型；所述第一計算模型為式3：

其中，為裂縫平均滲透率，md；為裂縫平均孔隙度，％；為裂縫段平均流體移動指數，us/ft；a₄、b₄、c₄為系數，通過第一關系模型和第二關系模型采用多元回歸擬合獲得；

所述根據試井滲透率對所述第一計算模型進行標定，包括：根據裂縫孔隙度、流體移動指數以及所述第一計算模型，得到第一參照裂縫滲透率；根據所述第一參照裂縫滲透率、裂縫孔隙度、流體移動指數、基質滲透率，通過式1，得到參照試井滲透率；比較所述參照試井滲透率和試井滲透率，線性調整所述第一計算模型中的模型系數a₄為M，使所述參照試井滲透率和試井滲透率相等，得到標定計算模型；所述標定計算模型為式4：

建立裂縫滲透率、裂縫孔隙度、流體移動指數間的第二計算模型，包括：根據所述裂縫孔隙度、流體移動指數以及標定計算模型，得到第二參照裂縫滲透率；建立所述第二參照裂縫滲透率與裂縫孔隙度的第三關系模型，建立所述第二參照裂縫滲透率與流體移動指數的第四關系模型；根據所述第三關系模型和第四關系模型，利用多元回歸法，建立所述第二計算模型；所述第二計算模型為式5：

K_l＝A*FVPA^B*QFM^C 式5

其中，K_l為裂縫滲透率，md；FVPA為裂縫孔隙度，％；QFM為流體移動指數，us/ft；A、B、C為系數，通過第三關系模型和第四關系模型采用多元回歸擬合獲得。

2.根據權利要求1所述的裂縫滲透率的計算方法，其特征在于，通過電成像測井得到所述裂縫孔隙度，通過偶極聲波成像測井得到所述流體移動指數。

3.根據權利要求2所述的裂縫滲透率的計算方法，其特征在于，所述通過偶極聲波成像測井得到所述流體移動指數，包括：

根據式6-式7得到所述流體移動指數；

QFM＝S-S_e 式6

其中，DEN為密度測井曲線，g/cm³；ρ_f為泥漿密度，g/cm³；Δt_f為泥漿時差，us/ft；Δt_S為地層橫波時差，us/ft；S_e為理論斯通利波時差，us/ft；S為實測斯通利波時差，us/ft；QFM為流體移動指數，us/ft。

4.根據權利要求1所述的裂縫滲透率的計算方法，其特征在于，所述基質滲透率通過式8得到：

其中，K_d為基質滲透率，md；為基質孔隙度，％；a₂、b₂為系數，通過巖石物理物性實驗，采用最小二乘法擬合獲得。