隨機(jī)生長四參數(shù)生成孔隙模型及其分形描述

朱維偉　姚約東　鄒　威　王德強(qiáng)

中國石油大學(xué)(北京)石油工程教育部重點(diǎn)實驗室,　北京　昌平　102249

?

隨機(jī)生長四參數(shù)生成孔隙模型及其分形描述

朱維偉姚約東鄒威王德強(qiáng)

中國石油大學(xué)(北京)石油工程教育部重點(diǎn)實驗室,北京昌平102249

摘要：儲層孔隙是油氣儲存和運(yùn)移的主要通道,孔隙結(jié)構(gòu)的研究在油氣田開發(fā)中占有十分重要的位置。通過隨機(jī)生長四參數(shù)算法,建立了21組不同孔隙度、相同概率生長方向和若干組相同孔隙度、具有優(yōu)勢生長方向的孔隙模型,定量計算了各個孔隙模型的分形盒維數(shù)及其多重分形譜,對比研究孔隙結(jié)構(gòu)與分形維數(shù)相關(guān)關(guān)系的結(jié)果表明孔隙介質(zhì)具有很好的分形性質(zhì),孔隙度與盒維數(shù)在雙對數(shù)表示下,具有高度線性相關(guān)性。多重分形譜特征參數(shù)?Δα?與孔隙度相對獨(dú)立,能夠反映孔隙結(jié)構(gòu)的非均質(zhì)性,其值越大,說明孔隙結(jié)構(gòu)的非均質(zhì)性越強(qiáng)。

關(guān)鍵詞：隨機(jī)生長四參數(shù)；孔隙度；分形；盒維數(shù)；多重分形

0前言

孔隙微觀結(jié)構(gòu)的研究對石油與天然氣開發(fā)有重要意義。常規(guī)的孔隙結(jié)構(gòu)研究主要有毛管壓力曲線法,核磁共振法,掃描電鏡法,CT掃描法和測井方法等[1-3],這些實驗或者現(xiàn)場測試方法一般耗時長,費(fèi)用高,過程復(fù)雜。本文采用隨機(jī)生長四參數(shù)法建立孔隙模型,具有操作簡單、結(jié)果直觀、后續(xù)計算方便的優(yōu)勢。儲層孔隙被認(rèn)為是具有統(tǒng)計意義上的分形體,但通常是基于理論推導(dǎo),或二維圖像的計算,缺少三維儲層的分形計算。對于孔隙結(jié)構(gòu)與分形參數(shù)的相關(guān)關(guān)系,也討論甚少,對于局部孔隙結(jié)構(gòu)的表征,也是常規(guī)分形所不能表達(dá)的。本文基于單重和多重分形盒維數(shù)的計算方法,深入研究了孔隙結(jié)構(gòu)和分形維數(shù)的相關(guān)關(guān)系,并且利用多重分形進(jìn)行孔隙結(jié)構(gòu)非均質(zhì)性的表征。

1隨機(jī)生長四參數(shù)法建立孔隙模型

本文采用隨機(jī)生長四參數(shù)算法[4-5],進(jìn)行隨機(jī)孔隙結(jié)構(gòu)建模,建立由孔隙和固體顆粒兩相構(gòu)成的三維多孔介質(zhì)模型,其中孔隙為生長相,固體顆粒為非生長相,初始相全為固體顆粒。具體構(gòu)造步驟如下：

1)本文構(gòu)建立方體區(qū)域作為孔隙模型的生長區(qū)域。在構(gòu)造區(qū)域內(nèi)隨機(jī)分布孔隙節(jié)點(diǎn)生長核,其分布概率為pd,對于構(gòu)造區(qū)域內(nèi)每一個網(wǎng)格節(jié)點(diǎn)在[0,1]區(qū)間內(nèi)生成平均分布隨機(jī)數(shù),隨機(jī)數(shù)不大于pd的節(jié)點(diǎn)即為生長核。

2)對于初始隨機(jī)分布的孔隙節(jié)點(diǎn)進(jìn)行遍歷,將每個孔隙節(jié)點(diǎn)向三維空間進(jìn)行隨機(jī)生長。其生長方向涉及到26個方向,包括6個主方向,12個面對角方向,8個體對角方向。本文算法是對初始孔隙節(jié)點(diǎn)和生長之后的孔隙節(jié)點(diǎn)循環(huán)遍歷,不涉及到單獨(dú)節(jié)點(diǎn)連續(xù)生長的情況,因此為了生成各個方向生長均勻的孔隙模型,各個方向的生長概率ps均相同,以孔隙節(jié)點(diǎn)為生長原點(diǎn),向26個方向上的非孔隙鄰點(diǎn)重新生成隨機(jī)數(shù),當(dāng)生成的隨機(jī)數(shù)小于ps時,該點(diǎn)生長為孔隙節(jié)點(diǎn)。

按照上述算法,得到的多孔介質(zhì)模型見圖1。

圖1　孔隙度為0.15的多孔介質(zhì)模型

2.1孔隙結(jié)構(gòu)及其分形維數(shù)

Mandelbrot提出的分形概念[6]已經(jīng)廣泛運(yùn)用于多學(xué)科的理論研究和工程應(yīng)用[7-12]。分形體的根本性質(zhì)是自相似性,指局部和整體有某種相似的形,包括嚴(yán)格自相似和統(tǒng)計意義上的自相似兩種情況。長期以來,儲層孔隙被認(rèn)為是具有統(tǒng)計意義上的分形體[13-17],本文通過具體的三維孔隙結(jié)構(gòu)計算其分形維數(shù)進(jìn)行驗證,并且進(jìn)一步探討孔隙度和分形維數(shù)的內(nèi)在關(guān)聯(lián)。

分形維數(shù)的定義中,由于盒計數(shù)維數(shù)求取方便,含義簡明直觀,容易理解,因此有較為廣泛的運(yùn)用,在描述孔隙分形中,也通常采用盒計數(shù)維數(shù)進(jìn)行表征[18],其求取公式如下：

(1)

根據(jù)式(1)計算21組不同孔隙度模型相應(yīng)的盒維數(shù),結(jié)果表明隨機(jī)孔隙模型具有自相似性,各模型相關(guān)系數(shù)均高達(dá)0.99以上,不同孔隙度下的盒維數(shù)D相關(guān)關(guān)系曲線如圖2。

圖2　孔隙度與盒維數(shù)相關(guān)曲線

從圖2可以看出,隨著孔隙度的增加,孔隙模型的分形維數(shù)也隨之增大,其增長模式接近線形,擬合程度達(dá)到0.986,進(jìn)一步做出孔隙度對數(shù)值和分形維數(shù)對數(shù)值的相關(guān)曲線(圖3)。

圖3　孔隙度與盒維數(shù)雙對數(shù)相關(guān)曲線

由圖3擬合結(jié)果可以看出,孔隙度和盒維數(shù)的雙對數(shù)值成高度的線性相關(guān)性,相關(guān)系數(shù)達(dá)到0.999,因此儲層孔隙分形維數(shù)對數(shù)值可以嚴(yán)格表征儲層的孔隙度大小。

2.2孔隙結(jié)構(gòu)的多重分形描述

孔隙結(jié)構(gòu)的盒維數(shù)只能反映出儲層孔隙分布的大致情況,不能對孔隙分布的局部信息進(jìn)行表達(dá),因此引入多重分形對孔隙分布的局部信息進(jìn)行描述[19-20]。

多重分形指一個分形集由若干個子集組成,不同的標(biāo)度情況下,各個子集的分形維數(shù)也不同。多重分形的盒維數(shù)計數(shù)方法就是在考慮盒子內(nèi)像素數(shù)的差別,歸一化后得到一個概率分布集,再用一個多重分形譜進(jìn)行描述。對于孔隙結(jié)構(gòu)的多重分形描述,其關(guān)鍵點(diǎn)在于定義孔隙系統(tǒng)的概率分布函數(shù)如式(2)：

(2)

χq(r)=∑Pi(r)q=r(q)

(3)

由式(4)、(5)得到多重分形譜線：

(4)

(5)

2.3多重分形譜的求取

以孔隙度分別為0.15,0.25四組孔隙模型為例,繪制lnχq(r)-ln r曲線(圖4)。

a)孔隙度=0.15

b)孔隙度=0.25圖4　不同孔隙度模型lnχq(r)-ln r曲線

從圖4可以發(fā)現(xiàn),lnχq(r)和ln r成很好的線性關(guān)系,說明所建立的孔隙模型滿足分形標(biāo)度不變性質(zhì),是多重分形體,可以求解其多重分形譜。出于計算穩(wěn)定性的要求,本文在尺度r=1/32下,計算各個孔隙模型的多重分形譜線,見圖5。

圖5　不同孔隙度多重分形譜線

圖5中紅色曲線為孔隙度為0.15的多重分形譜線,藍(lán)色曲線為孔隙度為0.35的多重分形譜線,可以看出孔隙度從0.15至0.35的所有模型多重分形譜線都在由紅藍(lán)線構(gòu)成的包絡(luò)線中,因此,不同孔隙度模型的多重分形譜線具有相似性。以包絡(luò)線為代表進(jìn)行具體分析,見表1。

表1不同孔隙度模型的多重分形參數(shù)比較

孔隙度?D㊣0?Δα㊣?Δf(α)㊣?f(α)㊣max0.151.621.0952.7452.9870.351.521.1732.0912.998

2.4非均質(zhì)孔隙及其分形描述

由于天然沉積、地質(zhì)構(gòu)造等作用,天然的儲層孔隙結(jié)構(gòu)顯然與本文所建立的隨機(jī)孔隙結(jié)構(gòu)有明顯不同,孔隙的生長在各個方向并不是等概率生長,本文建立了孔隙度為0.25,存在優(yōu)勢生長方向的非均質(zhì)孔隙模型。通過提高優(yōu)勢生長方向的生長概率來實現(xiàn)優(yōu)勢生長,其優(yōu)勢方向,沿某一方向或沿某一平面生長,共建立非均質(zhì)模型5個,與等概率隨機(jī)生長模型對比,結(jié)果見表2。

表2存在優(yōu)勢生長方向的隨機(jī)孔隙模型分形參數(shù)對比

模型孔隙度盒維數(shù)?Δα㊣等概率0.252.87871.140219優(yōu)勢生長10.252.87101.160101優(yōu)勢生長20.252.87561.152628優(yōu)勢生長30.252.87301.142911優(yōu)勢生長40.252.87861.143530優(yōu)勢生長50.252.87571.150965

從表2可以看出,存在優(yōu)勢方向的非均質(zhì)孔隙結(jié)構(gòu),其盒維數(shù)均小于等概率生長的隨機(jī)孔隙模型,這是由于存在優(yōu)勢方向,即孔隙生長會更有向優(yōu)勢方向生長的趨勢,從而減弱了整體占滿整個構(gòu)造空間的趨勢,因此盒維數(shù)會有所有下降。由于盒維數(shù)與孔隙度有很強(qiáng)的相關(guān)性,因此盒維數(shù)并不能獨(dú)立反映出孔隙結(jié)構(gòu)的非均質(zhì)性,然而Δα的值相對于孔隙度是獨(dú)立的,能夠獨(dú)立反映出孔隙結(jié)構(gòu)的非均質(zhì)性,從表2可以看出,非均質(zhì)的孔隙結(jié)構(gòu)的Δα值更大。

3結(jié)論

1)孔隙結(jié)構(gòu)是分形體,能夠用分形理論進(jìn)行研究,孔隙結(jié)構(gòu)的分形維數(shù)與孔隙度成較好的線性關(guān)系,在雙對數(shù)坐標(biāo)下,兩者有相當(dāng)高的線性相關(guān)性,相關(guān)系數(shù)達(dá)到0.999以上。

2)分形盒維數(shù)與孔隙度關(guān)聯(lián)性大,但只能反映孔隙結(jié)構(gòu)的總體狀況,不能反映局部分布結(jié)構(gòu)和非均質(zhì)性。

3)孔隙結(jié)構(gòu)具備多重分形的標(biāo)度不變性質(zhì),是多重分形體,其特征參數(shù)Δα與孔隙度相互獨(dú)立,能夠有效描述孔隙結(jié)構(gòu)的局部信息和非均質(zhì)性,其值越大,則表示孔隙結(jié)構(gòu)非均質(zhì)性越強(qiáng)。

參考文獻(xiàn)：

[1] 郝樂偉,王琪,唐俊.儲層巖石微觀孔隙結(jié)構(gòu)研究方法與理論綜述[J].巖性油氣藏,2013,25(5):123-128.

Hao Lewei, Wang Qi, Tang Jun. Research Progress of Reservoir Microscopic Pore Structure [J]. Lithologic Reservoirs, 2013, 25 (5): 123-128.

[2] 賴錦,王貴文,陳敏,等.基于巖石物理相的儲集層孔隙結(jié)構(gòu)分類評價—以鄂爾多斯盆地姬塬地區(qū)長8油層組為例[J].石油勘探與開發(fā),2013,40(5):566-572.

Lai Jin, Wang Guiwen, Chen Min, et al. Pore Structures Evaluation of Low Permeability Clastic Reservoirs Based on Petrophysical Facies: A Case Study on Chang 8 Reservoir in Jiyuan Region, Ordos Basin [J]. Petroleum Exploration and Development, 2013, 40 (5): 566-572.

[3] 李海燕,岳大力,張秀娟.蘇里格氣田低滲透儲層微觀孔隙結(jié)構(gòu)特征及其分類評價方法[J].地學(xué)前緣,2012,19(2):133-139.

Li Haiyan, Yue Dali, Zhang Xiujuan. Characteristics of Pore Structure and Reservoir Evaluation of Low Permeability Reservoir in Sulige Gas Field [J]. Science Frontiers, 2012, 19 (2): 133-139.

[4] 李仁民,劉松玉,方磊,等.采用隨機(jī)生長四參數(shù)生成法構(gòu)造黏土微觀結(jié)構(gòu)[J].浙江大學(xué)學(xué)報：(工學(xué)版),2010,10(10):1897-1901.

Li Renmin, Liu Songyu, Fang Lei, et al. Micro-Structure of Clay Generated by Quartet Structure Generation Set [J]. Journal of Zhejiang University (Engineering Science), 2010, 10 (10): 1897-1901.

[5] Wang M, Pan N. Numerical Analyses of Effective Dielectric Constant of Multiphase Microporous Media [J]. Journal of Applied Physics, 2007, 101 (11): 114102.

[6] 孫霞,吳自勤.分形原理及應(yīng)用[M].合肥：中國科學(xué)技術(shù)大學(xué)出版社,2003：53-74.

Sun Xia, Wu Ziqing. The Principle and Application of Fractal [M]. Hefei: University of Science and Technology of China Press, 2003： 53-74.

[7] 李云開,楊培嶺,任樹梅,等.土壤水分與溶質(zhì)運(yùn)移機(jī)制的分形理論研究進(jìn)展[J].水科學(xué)進(jìn)展,2005,16(6):892-899

Li Yunkai, Yang Peiling, Ren Shumei, et al. Development in Research on the Fractal Theory of Soil Water and Solute Transportation Mechanisms [J]. Advance in Water Science, 2005, 16 (6): 892-899

[8] 王國梁,周生路,趙其國.土壤顆粒的體積分形維數(shù)及其在土地利用中的應(yīng)用[J].土壤學(xué)報,2005,42(4):545-550

Wang Guoliang, Zhou Shenglu, Zhao Qiguo. Volume Fractal Dimension of Soil Particles and Its Applications to Land Use [J]. Acta Pedologica Sinica, 2005, 42 (4): 545-550

[9] 劉勇軍,朱岳明,豐土根,等.隨機(jī)土體結(jié)構(gòu)模型理論及應(yīng)用[J].巖土工程學(xué)報,2001,23(3):354-357.Liu Yongjun, Zhu Yueming, Feng Tugen, et al. Study on Random Soil Frame Model and Its Application [J]. Chinese Journal of Geotechnical Engineering, 2001, 23 (3): 354-357.

[10] 王唯威.分形多孔介質(zhì)內(nèi)導(dǎo)熱與流動數(shù)值模擬研究[D].北京：中國科學(xué)院工程熱物理研究所,2006.

Wang Weiwei. Numerical Study of Heat Conduction in Fractal Porous Media [D]. Beijing: Institute of Engineering Thermophysics, CAS, 2006.

[11] 楊培嶺,羅遠(yuǎn)培,石元春.用粒徑的重量分布表征的土壤分形特征[J].科學(xué)通報,1993,38(20):1896-1899.

Yang Peiling, Luo Yuanpei, Shi Yuanchun. Soil Fractal Characteristics Measured by Mass of Particle-Size Distribution [J]. Chinese Science Bulletin, 1993, 38 (20): 1896-1899.

[12] 黃冠華,詹衛(wèi)華.土壤顆粒的分形特征及其應(yīng)用[J].土壤學(xué)報,2002,39(4):490-497.

Huang Guanhua, Zhan Weihua. Fractal Property of Soil Particle Size Distribution and Its Application [J]. Acta Pedologica Sinica, 2002, 39 (4): 490-497.

[13] 馬新仿,張士誠,郎兆新.儲層巖石孔隙結(jié)構(gòu)的分形研究[J].中國礦業(yè),2003,12(9):46-48.

Ma Xinfang, Zhang Shicheng, Lang Zhaoxin. Fractal Research on Pore Structure in Reservoir Rock [J]. China Mining Magazine, 2003, 12 (9): 46-48.

[14] 張憲國,張濤,林承焰.基于孔隙分形特征的低滲透儲層孔隙結(jié)構(gòu)評價[J].巖性油氣藏2013,25(6):40-45.

Zhang Xianguo, Zhang Tao, Lin Chengyan. Pore Structure Evaluation of Low Permeability Reservoir Based on Pore Fractal Features [J]. Lithologic Reservoirs, 2013, 25 (6): 40-45.

[15] 賈芬淑,沈平平,李克文.砂巖孔隙結(jié)構(gòu)的分形特征及應(yīng)用研究[J].斷塊油氣田,1995,2(1):16-21.Jia Fenshu, Shen Pingping, Li Kewen. Study on the Fracture Characteristics of Sandstone Pore Structure and Its Application [J]. Fault-Block Oil & Gas Field, 1995, 2 (1): 16-21.

[16] 楊建,陳家軍,楊周喜,等.松散砂?？紫督Y(jié)構(gòu)、孔隙分形特征及滲透率研究[J].水文地質(zhì)工程地質(zhì),2008,35(3):93-96.

Yang Jian, Chen Jiajun, Yang Zhouxi et al. A Study of Pore Structure, Pore Fractal Feature and Permeability of Unconsolidated Sand [J]. Hydrogeology and Engineering Geology, 2008, 35 (3): 93-96.

[17] 張東暉,楊浩,施明恒.多孔介質(zhì)分形模型的難點(diǎn)與探索[J].東南大學(xué)學(xué)報,2002,32(5):692-697.

Zhang Donghui, Yang Hao, Shi Mingheng.Important Prob

lems of Fractal Model in Porous Media [J]. Journal ofSoutheast University, 2002, 32 (5): 692-697.

[18] 彭瑞東,謝和平,鞠楊.二維數(shù)字圖像分形維數(shù)的計算方法[J].中國礦業(yè)大學(xué)學(xué)報,2004,33(1):19-24.

Peng Ruidong, Xie Heping, Ju Yang. Calculation Method of Fractal Dimension for 2-D Digital Image [J]. Journal of China University of Mining & Technology, 2004, 33 (1): 19-24.

[19] 管孝艷,楊培嶺,任樹梅,等.基于多重分形理論的壤土粒徑分布非均勻性分析[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報,2009,17(2):196-205.

Guan Xiaoyan, Yang Peiling, Ren Shumei, et al. Heterogeneity Analysis of Particle Size Distribution for Loamy Soil Based on Multifractal Theory [J]. Journal of Basic Science and Engineering, 2009, 17 (2): 196-205.

[20] 周煒星,吳韜,于遵宏.多重分形奇異譜的幾何特性Ⅱ.配分函數(shù)法[J].華東理工大學(xué)學(xué)報,2000,26(4):390-395

Zhou Weixing, Wu Tao, Yu Zundong. Geometrical Characteristics of Singularity Spectra of Multifractals Ⅱ Partition Function Definition [J]. Journal of East China University of Science and Technology, 2000, 26 (4): 390-395.

收稿日期：2015-12-29

基金項目：國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃(973計劃)(2015 CB 250900)；新世紀(jì)優(yōu)秀人才支持計劃(NCET-13-1030)

作者簡介：朱維偉(1993-),男,江蘇南京人,碩士研究生,主要從事油氣田開發(fā)研究工作。

DOI:10.3969/j.issn.1006-5539.2016.03.013