[發明專利]一種縫洞型儲層試井確定地層中洞流量方法

有效

• 申請號: 201810745121.7
• 申請日: 2018-07-09
• 公開（公告）號: CN108920849B
• 公開（公告）日: 2022-07-19
• 發明（設計）人: 盧志煒;程誠;李越
• 申請（專利權）人: 合肥辰工科技有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F30/28;G06F113/08;G06F119/14;G06F113/14
• 代理公司: 合肥洪雷知識產權代理事務所(普通合伙) 34164
• 代理人: 張悅

1.一種縫洞型儲層試井確定地層中洞流量方法，其特征在于，包括采用基于能量守恒定律及溶洞波動相結合的縫洞型儲層試井解釋方法，具體步驟如下：

S000：建立縫洞型油藏試井模型，并將地層中的溶洞視作點源；

S001：根據所述試井模型，獲得無限大地層中的真實空間井底壓力解；

S002：根據所述真實空間井底壓力解，采用疊加原理獲得地層中多個洞同時存在時的真實空間井底壓力；

S003：利用所述真實空間井底壓力解與實測的井底壓力數據進行擬合，獲得洞的個數及洞到井的距離；

其中S001具體步驟包括如下過程：

A000：利用拉普拉斯變換對所述試井模型進行求解，獲得拉普拉斯空間井底壓力解函數

A001：利用Stehfest數值反演技術，根據所述拉普拉斯空間井底壓力解函數得到真實空間井底壓力解p_wfD(t_D)；

其中，所述試井模型為：

所述試井模型中的油藏參數包括：原始地層壓力：p_i；裂縫壓力：p_f；溶洞壓力：p_v；井底壓力：p_wf；外部介質滲透率：k；流體粘度：μ；外部介質壓縮系數：C_t；產量：Q；體積系數：B；井筒半徑：r_w；溶洞半徑：r_v；井筒儲集系數：C_w；為溶洞儲集常數：C_v；井筒的表皮系數：s_w；溶洞的表皮系數：s_v；與井筒相連地層的厚度：h₁；與溶洞相連地層的厚度：h₂；時間：t；半徑：r；溶洞無因次半徑：r_vD；

其中，

無因次時間：

無因次半徑：

無因次井筒儲集系數：

無因次溶洞儲集系數：

無因次高度：

無因次波動系數：

無因次阻尼系數：

方程修正系數：α；

數學函數：

其中，k為外部介質滲透率；μ為流體粘度；v₀為流體初始時刻的速度；D為油管直徑；ρ為流體密度；C為管道及流體系統中的波速；

S002包括如下步驟：

B000：確定地層中向試井中提供產量Q1的洞1以及提供產量Q2的洞2；

B001：利用疊加原理得到洞1和洞2同時存在時的真實空間井底壓力解；

其中，利用疊加原理得到洞1和洞2同時存在時的真實空間井底壓力解包括如下過程：

C000：將所述洞1看作點源1、洞2看作點源2；按照源匯理論得到地層中多源匯的連續性方程為：

(1)

其中，q_i為第i個源匯的強度；M₁為點源1在空間上的位置；M₂為點源2在空間上的位置；M為空間上任意點位置；

δ函數表示集中分布的物理量；

C001：采用微可壓縮假設，并將達西定律代入方程(1)得到：

C002：定義無因次量并求解無因次源匯方程，得到源匯在井筒處的無因次壓力：

其中，為無因次壓力；為無因次時間；為洞1的無因次距離；為洞2的無因次距離；為指數積分函數；為洞2提供的流量比；為洞1提供的流量比；

C003：洞1和洞2同時存在時的試井的真實空間井底壓力為：

S003中利用所述真實空間井底壓力解與實測的井底壓力數據進行擬合，獲得洞的個數具體過程如下：

利用所述真實空間井底壓力解p_wfD(t_D)與實測的井底壓力數據進行擬合，獲得無因次流量比Q_iD，(1)由公式Q_i＝Q_iD×Q計算洞提供的流量Q_i；(2)從實測曲線中上翹的段數為洞的個數。