[發(fā)明專利]基于生產(chǎn)動(dòng)態(tài)自動(dòng)歷史擬合的壓裂裂縫參數(shù)反演方法

在審

• 申請(qǐng)?zhí)枺?/td>: 202011342448.3
• 申請(qǐng)日: 2020-11-25
• 公開（公告）號(hào): CN112539049A
• 公開（公告）日: 2021-03-23
• 發(fā)明（設(shè)計(jì)）人: 姚軍;張凱;李振宇;樊靈;劉凡;翟亮;徐耀東;張黎明;劉均榮;姚傳進(jìn)
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油大學(xué)（華東）
• 主分類號(hào): E21B43/26
• 分類號(hào): E21B43/26;E21B47/00;E21B49/00;G06T17/00
• 代理公司: 北京匯澤知識(shí)產(chǎn)權(quán)代理有限公司 11228
• 代理人: 關(guān)宇辰

1.一種基于生產(chǎn)動(dòng)態(tài)自動(dòng)歷史擬合的壓裂裂縫參數(shù)反演方法，其特征在于，包括以下步驟：

步驟一、應(yīng)用嵌入式離散裂縫模型對(duì)儲(chǔ)層進(jìn)行剖分，并計(jì)算傳遞率耦合裂縫與基質(zhì)間的流動(dòng)

步驟二、建立三維地震監(jiān)測(cè)模型，模擬儲(chǔ)層變化

步驟三、根據(jù)裂縫形態(tài)和離散裂縫反演算法建立復(fù)雜裂縫網(wǎng)絡(luò)生成方法

步驟四、構(gòu)建離散網(wǎng)絡(luò)確定性反演算法。

2.根據(jù)權(quán)利要求1所述的基于生產(chǎn)動(dòng)態(tài)自動(dòng)歷史擬合的壓裂裂縫參數(shù)反演方法，其特征在于，步驟一過程如下：

S1.1：嵌入式離散裂縫模型剖分儲(chǔ)層基質(zhì)和裂縫

將儲(chǔ)層分成裂縫和基質(zhì)兩個(gè)系統(tǒng)并且儲(chǔ)層由結(jié)構(gòu)化網(wǎng)格剖分，根據(jù)裂縫與不同的區(qū)域連接情況將傳遞率分為三個(gè)類型，即裂縫段與儲(chǔ)層基質(zhì)之間連接、不同裂縫段之間連接、同一條裂縫段之間連接，嵌入式離散裂縫模型使用結(jié)構(gòu)化網(wǎng)格剖分儲(chǔ)層基質(zhì)，對(duì)于具有復(fù)雜形態(tài)的裂縫使用了非結(jié)構(gòu)化網(wǎng)格進(jìn)行剖分，二維模型中裂縫的剖分也是使用結(jié)構(gòu)化網(wǎng)格；

S1.2：計(jì)算傳遞率耦合裂縫與基質(zhì)間的流動(dòng)

使用傳遞率大小來近似衡量流體在裂縫中的流動(dòng)，定義如下：

(q_j)_i+1/2＝(T_j)_i+1/2((p_j)_i+1-(p_j)_i) (1)

式中，i+1和i表示相鄰網(wǎng)格屬性的下標(biāo)；i表示相鄰網(wǎng)格的交界面；j表示相的類型；p表示壓力；T表示傳遞率系數(shù)。三種傳遞率類型的具體計(jì)算公式如下：

裂縫網(wǎng)格與基質(zhì)網(wǎng)格連接的傳遞率公式為：

q_mf＝T_mf(p_f-p_m) (2)

式中，mf表示裂縫網(wǎng)格與儲(chǔ)層基質(zhì)網(wǎng)格連接的下標(biāo)；K_mf定義為裂縫網(wǎng)格與其相對(duì)應(yīng)的基質(zhì)網(wǎng)格滲透率的調(diào)和平均數(shù)；A_mf表示裂縫網(wǎng)格與基質(zhì)網(wǎng)格交界面的面積；d_mf表示基質(zhì)網(wǎng)格與裂縫網(wǎng)格的平均距離，該值計(jì)算公式如下：

式中，n表示垂直于裂縫平面的單位法向量；V表示基質(zhì)網(wǎng)格的體積；x表示基質(zhì)網(wǎng)格體積元素到裂縫網(wǎng)格的距離；

不同裂縫網(wǎng)格之間連接的傳遞率為每條裂縫傳遞率指數(shù)的調(diào)和平均值，公式為：

式中，下標(biāo)f和f’表示不同的裂縫網(wǎng)格；k_i表示不同裂縫網(wǎng)格的滲透率；w_i表示不同裂縫的開度；L表示不同裂縫相交線的長(zhǎng)度；d_i表示不同裂縫網(wǎng)格的中心到裂縫相交線的距離；

同一條裂縫網(wǎng)格之間的連接的傳遞率公式：

式中，k_f表示裂縫的滲透率；A_f表示同一條裂縫不同裂縫網(wǎng)格交界面的面積，通常等于裂縫開度與裂縫網(wǎng)格邊界長(zhǎng)度的乘積；d_f表示同一條裂縫不同裂縫網(wǎng)格中心之間的距離。