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

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

?

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

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

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

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

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

0前言

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

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

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

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

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

按照上述算法,得到的多孔介質(zhì)模型見(jiàn)圖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)計(jì)意義上的自相似兩種情況。長(zhǎng)期以來(lái),儲(chǔ)層孔隙被認(rèn)為是具有統(tǒng)計(jì)意義上的分形體[13-17],本文通過(guò)具體的三維孔隙結(jié)構(gòu)計(jì)算其分形維數(shù)進(jìn)行驗(yàn)證,并且進(jìn)一步探討孔隙度和分形維數(shù)的內(nèi)在關(guān)聯(lián)。

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

(1)

根據(jù)式(1)計(jì)算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ù)也隨之增大,其增長(zhǎng)模式接近線形,擬合程度達(dá)到0.986,進(jìn)一步做出孔隙度對(duì)數(shù)值和分形維數(shù)對(duì)數(shù)值的相關(guān)曲線(圖3)。

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

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

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

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

多重分形指一個(gè)分形集由若干個(gè)子集組成,不同的標(biāo)度情況下,各個(gè)子集的分形維數(shù)也不同。多重分形的盒維數(shù)計(jì)數(shù)方法就是在考慮盒子內(nèi)像素?cái)?shù)的差別,歸一化后得到一個(gè)概率分布集,再用一個(gè)多重分形譜進(jìn)行描述。對(duì)于孔隙結(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)系,說(shuō)明所建立的孔隙模型滿足分形標(biāo)度不變性質(zhì),是多重分形體,可以求解其多重分形譜。出于計(jì)算穩(wěn)定性的要求,本文在尺度r=1/32下,計(jì)算各個(gè)孔隙模型的多重分形譜線,見(jiàn)圖5。

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

圖5中紅色曲線為孔隙度為0.15的多重分形譜線,藍(lán)色曲線為孔隙度為0.35的多重分形譜線,可以看出孔隙度從0.15至0.35的所有模型多重分形譜線都在由紅藍(lán)線構(gòu)成的包絡(luò)線中,因此,不同孔隙度模型的多重分形譜線具有相似性。以包絡(luò)線為代表進(jìn)行具體分析,見(jià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)造等作用,天然的儲(chǔ)層孔隙結(jié)構(gòu)顯然與本文所建立的隨機(jī)孔隙結(jié)構(gòu)有明顯不同,孔隙的生長(zhǎng)在各個(gè)方向并不是等概率生長(zhǎng),本文建立了孔隙度為0.25,存在優(yōu)勢(shì)生長(zhǎng)方向的非均質(zhì)孔隙模型。通過(guò)提高優(yōu)勢(shì)生長(zhǎng)方向的生長(zhǎng)概率來(lái)實(shí)現(xiàn)優(yōu)勢(shì)生長(zhǎng),其優(yōu)勢(shì)方向,沿某一方向或沿某一平面生長(zhǎng),共建立非均質(zhì)模型5個(gè),與等概率隨機(jī)生長(zhǎng)模型對(duì)比,結(jié)果見(jiàn)表2。

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

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

從表2可以看出,存在優(yōu)勢(shì)方向的非均質(zhì)孔隙結(jié)構(gòu),其盒維數(shù)均小于等概率生長(zhǎng)的隨機(jī)孔隙模型,這是由于存在優(yōu)勢(shì)方向,即孔隙生長(zhǎng)會(huì)更有向優(yōu)勢(shì)方向生長(zhǎng)的趨勢(shì),從而減弱了整體占滿整個(gè)構(gòu)造空間的趨勢(shì),因此盒維數(shù)會(huì)有所有下降。由于盒維數(shù)與孔隙度有很強(qiáng)的相關(guān)性,因此盒維數(shù)并不能獨(dú)立反映出孔隙結(jié)構(gòu)的非均質(zhì)性,然而Δα的值相對(duì)于孔隙度是獨(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)系,在雙對(duì)數(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] 郝樂(lè)偉,王琪,唐俊.儲(chǔ)層巖石微觀孔隙結(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] 賴錦,王貴文,陳敏,等.基于巖石物理相的儲(chǔ)集層孔隙結(jié)構(gòu)分類評(píng)價(jià)—以鄂爾多斯盆地姬塬地區(qū)長(zhǎng)8油層組為例[J].石油勘探與開(kāi)發(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] 李海燕,岳大力,張秀娟.蘇里格氣田低滲透儲(chǔ)層微觀孔隙結(jié)構(gòu)特征及其分類評(píng)價(jià)方法[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ī)生長(zhǎng)四參數(shù)生成法構(gòu)造黏土微觀結(jié)構(gòu)[J].浙江大學(xué)學(xué)報(bào)：(工學(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].合肥：中國(guó)科學(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] 李云開(kāi),楊培嶺,任樹(shù)梅,等.土壤水分與溶質(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] 王國(guó)梁,周生路,趙其國(guó).土壤顆粒的體積分形維數(shù)及其在土地利用中的應(yīng)用[J].土壤學(xué)報(bào),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é)報(bào),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)熱與流動(dòng)數(shù)值模擬研究[D].北京：中國(guó)科學(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é)通報(bào),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é)報(bào),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] 馬新仿,張士誠(chéng),郎兆新.儲(chǔ)層巖石孔隙結(jié)構(gòu)的分形研究[J].中國(guó)礦業(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] 張憲國(guó),張濤,林承焰.基于孔隙分形特征的低滲透儲(chǔ)層孔隙結(jié)構(gòu)評(píng)價(jià)[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é)報(bào),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ì)算方法[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),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] 管孝艷,楊培嶺,任樹(shù)梅,等.基于多重分形理論的壤土粒徑分布非均勻性分析[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),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é)報(bào),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

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

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

DOI:10.3969/j.issn.1006-5539.2016.03.013