[發(fā)明專利]一種基于變徑毛細(xì)管束模型的滲吸驅(qū)替表征方法

在審

• 申請?zhí)枺?/td>: 202211078060.6
• 申請日: 2022-09-05
• 公開（公告）號(hào): CN115436257A
• 公開（公告）日: 2022-12-06
• 發(fā)明（設(shè)計(jì)）人: 高濤;王錳;王超;郭睿;付晨陽;倪軍;王維波;劉文強(qiáng);楊亞潔;高怡文;薛文卓;管雅倩;孫淑娟;田偉偉;李紅
• 申請（專利權(quán)）人: 陜西延長石油（集團(tuán)）有限責(zé)任公司
• 主分類號(hào): G01N15/08
• 分類號(hào): G01N15/08
• 代理公司: 西安億諾專利代理有限公司 61220
• 代理人: 何春江

1.一種基于變徑毛細(xì)管束模型的滲吸驅(qū)替表征方法，其特征在于，包括以下步驟：

(1)根據(jù)真實(shí)巖心實(shí)驗(yàn)確定變徑毛細(xì)管生成參數(shù)，包括毛細(xì)管長度L、半徑分布數(shù)據(jù)集r以及迂曲度范圍τ；

(2)根據(jù)步驟(1)中的參數(shù)，采用隨機(jī)算法建立巖心尺度隨機(jī)變徑毛細(xì)管束；

(3)根據(jù)步驟(2)中得到的隨機(jī)變徑毛細(xì)管束，分別計(jì)算毛細(xì)管束中各個(gè)毛細(xì)管各段的空間坐標(biāo)；

(4)根據(jù)步驟(2)、步驟(3)得到的隨機(jī)變徑毛細(xì)管結(jié)構(gòu)參數(shù)對(duì)滲吸速度進(jìn)行求解；

(5)根據(jù)步驟(4)的計(jì)算結(jié)果，求解隨機(jī)變徑毛細(xì)管束采出程度、含水率。

2.如權(quán)利要求1所述一種基于變徑毛細(xì)管束模型的滲吸驅(qū)替表征方法，其特征在于，所述步驟(1)需確定的變徑毛細(xì)管生成參數(shù)如下：

根據(jù)真實(shí)巖心長度確定毛細(xì)管長度L；

毛細(xì)管段數(shù)范圍k，且k＝(L，5L)；

根據(jù)真實(shí)巖心壓汞或核磁數(shù)據(jù)，確定半徑分布數(shù)據(jù)集r，該數(shù)據(jù)集r生成頻率之和為1；

根據(jù)真實(shí)巖心實(shí)驗(yàn)數(shù)據(jù)確定迂曲度范圍τ；

并采用下式計(jì)算毛細(xì)管真實(shí)流動(dòng)路徑L_e；

3.如權(quán)利要求2所述一種基于變徑毛細(xì)管模型的滲吸驅(qū)替表征方法，其特征在于，所述步驟(2)根據(jù)以下方式建立巖心尺度變徑毛細(xì)管束：

(a)根據(jù)毛細(xì)管束n及毛細(xì)管段數(shù)范圍k，生成段數(shù)列表K_L＝{k₁,k₂,k₃…k_n}，其中k_i∈k，則該毛細(xì)管束共具有n_k＝∑k_i個(gè)微管段；

(b)根據(jù)毛細(xì)管束n及迂曲度范圍τ生成迂曲度列表τ_L＝{τ₁,τ₂,τ₃…τ_n}，其中τ_i∈τ；

(c)根據(jù)半徑分布數(shù)據(jù)集r和nk值，生成半徑列表R_L＝{r₁,r₂,r₃…r_nk}，其中r_i∈r；

(d)對(duì)任意一個(gè)毛細(xì)管i，從迂曲度列表τ_L中為其指定一個(gè)τ_i，從段數(shù)列表K_L中為其指定一個(gè)k_i，并從半徑列表R_L取k_i個(gè)值賦予各段，即可得到一個(gè)具有空間坐標(biāo)且段數(shù)及半徑隨機(jī)分布的隨機(jī)變徑毛細(xì)管束模型。

4.如權(quán)利要求3所述一種基于變徑毛細(xì)管模型的滲吸驅(qū)替表征方法，其特征在于，所述步驟(3)根據(jù)以下方式確定隨機(jī)變徑毛細(xì)管束模型中每一根隨機(jī)變徑毛細(xì)管各段空間坐標(biāo)：對(duì)所有變徑毛細(xì)管進(jìn)行以下計(jì)算，將毛細(xì)管長度L分成k份得到單段前進(jìn)長度L_k，將當(dāng)前毛細(xì)管真實(shí)流動(dòng)路徑L_e分成k份得到單段長度L_e，k，采用下式計(jì)算各段對(duì)應(yīng)縱向偏轉(zhuǎn)量L_p，k

令變徑毛細(xì)管在x軸方向的單段增量Δx_i＝L_k，變徑毛細(xì)管在y軸方向的單段增量Δy_i＝L_p，k，變徑管各段結(jié)尾坐標(biāo)表示如下：