[發(fā)明專利]一種頁巖氣藏壓裂水平井產(chǎn)出剖面反演解釋方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202210226114.2
• 申請(qǐng)日: 2022-03-09
• 公開（公告）號(hào): CN114595504B
• 公開（公告）日: 2023-03-10
• 發(fā)明（設(shè)計(jì)）人: 羅紅文;向雨行;李海濤;于皓;張琴;劉暢;李穎;蔣貝貝;高素娟;聶松;馬欣;劉子民
• 申請(qǐng)（專利權(quán)）人: 西南石油大學(xué)
• 主分類號(hào): G06F30/13
• 分類號(hào): G06F30/13;E21B43/26;E21B47/00;G06F111/10;G06F119/08
• 代理公司: 四川北新律師事務(wù)所 51366
• 代理人: 李琳

1.一種頁巖氣藏壓裂水平井產(chǎn)出剖面反演解釋方法，其特征在于，包括以下步驟：

步驟S1：通過實(shí)測井筒參數(shù)，判斷有效人工裂縫位置以及數(shù)量N，再通過溫度剖面在人工裂縫的溫度變化，對(duì)各級(jí)裂縫半長進(jìn)行初步的估算，取值范圍的最小向量以及中值向量以及最大值向量確定誤差精度ε_T；

步驟S2：預(yù)估裂縫半長取值，第i條裂縫的裂縫半長取值范圍為[x_fa,i,x_fb,i],其中值為x_fc,i＝(x_fa,i+x_fb,i)/2,將代入溫度正演預(yù)測模型，計(jì)算出溫度剖面代入評(píng)價(jià)方程計(jì)算出擬合評(píng)價(jià)函數(shù)

步驟S3：判斷擬合評(píng)價(jià)函數(shù)是否滿足不滿足條件則進(jìn)行步驟S4，否則轉(zhuǎn)到步驟S8；

步驟S4：從趾端到跟端，逐級(jí)進(jìn)行反演，第i條裂縫半長解當(dāng)前所在區(qū)間為[x_fa,i,x_fb,i]，若x_fb,i-x_fa,i＜ε_T則轉(zhuǎn)步驟S6，不滿足則分別將向量[x_fc,1,x_fc,2,…,x_fc,i-1,x_fa,i,x_fc,i+1,…x_fc,N]、[x_fc,1,x_fc,2,…,x_fc,i-1,x_fb,i,x_fc,i+1,…x_fc,N]和代入溫度正演預(yù)測模型計(jì)算溫度剖面和計(jì)算出第i條裂縫處的溫度分別為和

步驟S5：根據(jù)實(shí)測溫度剖面上的i條裂縫處溫度分別計(jì)算第i條裂縫溫度反演誤差函數(shù)值和進(jìn)行判斷：

①若F_εT(x_fc,i)＝0，則進(jìn)行步驟S6；

②若F_εT(x_fa,i)·F_εT(x_fc,i)＜0，則對(duì)區(qū)間[x_fa,i,x_fc,i]進(jìn)行線性插值計(jì)算；

i如果則令進(jìn)行步驟S6；

ii如果則令x_fb,i＝x_fc,i，進(jìn)行步驟S4；

iii如果則令然后轉(zhuǎn)步驟S4；

③若F_εT(x_fa,i)·F_εT(x_fc,i)＞0，則對(duì)區(qū)間[x_fc,i,x_fb,i]進(jìn)行線性插值計(jì)算；

i如果則令進(jìn)行步驟S6；

ii如果則令x_fa,i＝x_fc,i、進(jìn)行步驟S4；

iii如果則令然后轉(zhuǎn)步驟S4；

步驟S6：第i條裂縫半長反演結(jié)果為x_fc,i，進(jìn)行第i+1條裂縫半長反演，用i+1替換i，轉(zhuǎn)至步驟S4；

步驟S7：重復(fù)步驟S4至步驟S6，直至完成所有裂縫半長反演，當(dāng)前反演結(jié)果為將代入溫度正演預(yù)測模型，然后代入評(píng)價(jià)方程計(jì)算目標(biāo)函數(shù)如滿足則轉(zhuǎn)步驟S8，否則重新賦值各級(jí)裂縫半長初值，然后轉(zhuǎn)步驟S1；

步驟S8：反演結(jié)束，輸出裂縫半長反演結(jié)果為然后將代入溫度正演預(yù)測模型，輸出產(chǎn)出剖面解釋。

2.如權(quán)利要求1所述的一種頁巖氣藏壓裂水平井產(chǎn)出剖面反演解釋方法，其特征在于，所述步驟S2和步驟S7中所述評(píng)價(jià)方程如下：

式中，為反演目標(biāo)參數(shù)向量，即裂縫半長，為反演溫度剖面計(jì)算值，為溫度剖面的實(shí)測值。

3.如權(quán)利要求1所述的一種頁巖氣藏壓裂水平井產(chǎn)出剖面反演解釋方法，其特征在于，所述步驟S5中②所述線性插值計(jì)算如下：

式中，為i點(diǎn)處的線性插值，x_fa,i與x_fc,i分別為i點(diǎn)處解的左區(qū)間點(diǎn)和區(qū)間中值點(diǎn)，F(xiàn)_εT為反演誤差函數(shù)。