王世芳，吳 濤，鄧永菊，戴 偉

（1.湖北第二師范學(xué)院 物理與電子信息學(xué)院，武漢 430205；2.武漢工程大學(xué) 理學(xué)院，武漢 430074）

隨機(jī)分布樹狀分叉網(wǎng)絡(luò)滲流特性的分形研究

王世芳1＊，吳 濤2，鄧永菊1，戴 偉1

（1.湖北第二師范學(xué)院 物理與電子信息學(xué)院，武漢 430205；2.武漢工程大學(xué) 理學(xué)院，武漢 430074）

基于Darcy定律和Hagen－Poiseuille方程，應(yīng)用分形理論研究了母管直徑分布具有分形分布的樹狀分叉網(wǎng)絡(luò)的滲流特性，得到了具有隨機(jī)分布的樹狀分叉網(wǎng)絡(luò)滲透率的解析表達(dá)式.研究結(jié)果表明：通過具有隨機(jī)分布的樹狀分叉網(wǎng)絡(luò)的流量與壓力梯度之間呈線性變化關(guān)系；其滲透率隨分形維數(shù)Df和孔隙率φ的增加而增加，隨長度比的增加而減小.

樹狀分叉網(wǎng)絡(luò)；滲透率；Darcy定律；壓力梯度

樹狀分叉網(wǎng)絡(luò)是自然界一種常見網(wǎng)絡(luò)結(jié)構(gòu)，在生物工程、石油開采工程、材料科學(xué)工程、微電子器件冷卻系統(tǒng)工程等實際應(yīng)用領(lǐng)域得到了廣泛的應(yīng)用，因此樹狀分叉網(wǎng)絡(luò)輸運(yùn)特性的研究引起眾多學(xué)者的關(guān)注［1-5］.

樹狀分叉網(wǎng)絡(luò)的研究始于1926年生理學(xué)家Murray對于血管系統(tǒng)血液流動優(yōu)化管道直徑比為2－1／3的提出［6-7］.Chen等人［8］提出了一個三維樹狀網(wǎng)絡(luò)結(jié)構(gòu)微通道散熱器熱流體動力學(xué)模型，他們的研究結(jié)果表明樹狀網(wǎng)絡(luò)的微通道散熱器在散熱和壓降損失等方面將比傳統(tǒng)結(jié)構(gòu)具有更大的優(yōu)勢.Xu和 Yu［9-10］采用兩種不同的方法研究了“點到線”型和“點到圓”型樹狀分叉網(wǎng)絡(luò)的滲流特性，還討論了迂曲度對網(wǎng)絡(luò)滲流特性的影響，但他們只研究了單個樹狀分叉網(wǎng)絡(luò)滲流特性.Wang和Yu［11］研究了非牛頓冪律流體在單個Y－Shape樹狀分叉網(wǎng)絡(luò)中的滲流特性，分析了網(wǎng)絡(luò)結(jié)構(gòu)參數(shù)對滲透率的影響.Li和Yu［12］研究非牛頓Bingham流體在單個規(guī)則樹狀分叉網(wǎng)絡(luò)中啟動壓力梯度，討論了網(wǎng)絡(luò)微結(jié)構(gòu)參數(shù)對啟動壓力梯度的影響.然而自然界中的分叉網(wǎng)絡(luò)一般都是不規(guī)則的，且大部分是隨機(jī)分布，僅僅用規(guī)則的樹狀分叉網(wǎng)絡(luò)來分析滲流特性往往是不夠的.本文以低滲透油藏中可能存在的裂縫網(wǎng)絡(luò)為研究背景，將母管直徑分布具有分形分布的“點到線”型樹狀分叉網(wǎng)絡(luò)嵌入到多孔介質(zhì)母體材料中，形成復(fù)合材料.假設(shè)這些分叉網(wǎng)絡(luò)在分叉過程中沒有重疊現(xiàn)象，分叉網(wǎng)絡(luò)具有相同的直徑比和長度比.為了研究問題的方便，本文忽略了母體介質(zhì)對滲流的影響，主要研究了母管直徑具有分形分布樹狀分叉網(wǎng)絡(luò)的滲流特性.

1 母管直徑具有分形分布的樹狀分叉網(wǎng)絡(luò)的流量

單個類分形樹狀分叉網(wǎng)絡(luò)一般是由初級到末級組成，設(shè)其最大分叉級數(shù)為m，每個母分支可以分成n個子分叉.圖1是單個“點到線”型樹狀分形分叉網(wǎng)絡(luò)的結(jié)構(gòu)示意圖［13］，第k級分叉的長度和直徑分別為lk和dk，分叉角度為θ，為了描述分叉網(wǎng)絡(luò)的幾何結(jié)構(gòu)，定義尺度因子γk和βk分別表示相鄰兩級分叉子管的長度比和直徑比：

圖1 點到線型樹狀分叉網(wǎng)絡(luò)示意圖［13］Fig.1 Fractal－like tree network model between apoint and a straight line［13］

如果單個樹狀分叉網(wǎng)絡(luò)具有分形特征，必須要求分叉長度比與直徑比保持常數(shù)［1］，即：γk＝γ和βk＝β.第k級分叉管道的總數(shù)為nk，其長度與直徑分別為［1］：

式中，l0和d0分別表示第0級管道的長度和直徑.對于二分叉網(wǎng)絡(luò)結(jié)構(gòu)，分叉數(shù)n＝2.

圖2 隨機(jī)分布的類分形樹狀分叉網(wǎng)絡(luò)嵌入到母體多孔介質(zhì)中（忽略母體滲流的影響）Fig.2 A porous matrix embedded with randomly distributed Y－shaped fractal－like tree networks（The seepage of matrix is ignored）

圖2為母管直徑具有分形分布的樹狀分叉網(wǎng)絡(luò)嵌入到母體多孔介質(zhì)中形成復(fù)合介質(zhì)示意圖，這里忽略了母體的滲流影響.根據(jù)Hagen－Poiseuille方程，通過第k級分叉單管的流量為

式中，Δpk為第k級單管兩端的壓強(qiáng)差，μ為牛頓流體的粘度.由（4）式可以得到流過第k級所有分叉管道的流量：

根據(jù)流體的質(zhì)量守恒定律，第k級分叉網(wǎng)絡(luò)的流量q（d0）也是通過單個樹狀分形分叉網(wǎng)絡(luò)的總流量.由（5）式可得第k級分叉網(wǎng)絡(luò)兩端的壓強(qiáng)降為：

如果忽略網(wǎng)絡(luò)分叉處的局部損失，網(wǎng)絡(luò)的總壓降Δp可近似等于各級壓降Δpk之和：

將方程（3）、（5）和（6）代入（7）式，得到分叉網(wǎng)絡(luò)的總壓降為：

根據(jù)方程（8），則流過母管直徑為d0的單個分叉網(wǎng)絡(luò)的流量重新改寫為：

假設(shè)分叉網(wǎng)絡(luò)的母管直徑滿足分形分布，則母管道直徑分布在d0和d0＋d（d0）區(qū)間內(nèi)的母管數(shù)目為：

根據(jù)Darcy定律：

式中，A和L0分別代表復(fù)合材料的橫截面積和直線距離長度，L0可以由下式計算得［9］

將方程（11）（13）代入（12）式，得到隨機(jī)分布的樹狀分叉網(wǎng)絡(luò)嵌入到母體介質(zhì)形成的復(fù)合介質(zhì)的滲透率的分形解析表達(dá)式：

對于二分叉網(wǎng)絡(luò)（n＝2），在保證分叉網(wǎng)絡(luò)橫截面積守恒的前提下，即每一個母管分叉出兩個相鄰子管道的橫截面積均相等，得到直徑比大小為

類似于多孔介質(zhì)孔隙率的定義［14］，定義母管直徑具有分形分布的樹狀分叉網(wǎng)絡(luò)鑲嵌到多孔介質(zhì)中形成復(fù)合材料的孔隙率：

由上式可以看出復(fù)合材料無量綱滲透率k＋是分叉管道長度比γ、分叉角θ以及孔隙率φ和分形維數(shù)Df的函數(shù)，該表達(dá)式不含任何經(jīng)驗常數(shù)，清楚地揭示了影響裂縫網(wǎng)絡(luò)滲透率的物理機(jī)制.

由此可以看出流量和滲透率對最大母管直徑非常敏感.當(dāng)m＝0或γ＝0時，每個樹狀分形分叉網(wǎng)絡(luò)只有第0級單管，此時裂縫網(wǎng)絡(luò)相當(dāng)于由一束獨立的、不同大小尺寸的平行直圓毛細(xì)管組成的，此時裂縫網(wǎng)絡(luò)流量和滲透率的表達(dá)式與Yu和Cheng［15］給出的結(jié)論一致.

圖3顯示了基于方程（19）在不同最大母管直徑下，通過裂縫網(wǎng)絡(luò)的流量與壓力梯度之間的關(guān)系.在總分叉級數(shù)m＝10、分叉角θ＝、分形維數(shù)Df＝1.1、孔隙率φ＝0.4、粘度μ＝1Pa·s下，通過裂縫網(wǎng)絡(luò)的流量隨壓力梯度的增加而線性增加；當(dāng)長度比γ＝0.6和直徑比β＝時，分叉網(wǎng)絡(luò)的最大母管直徑越大，通過具有分形分布樹狀分叉網(wǎng)絡(luò)的流量也越大.

圖3 在不同最大母管直徑下流量與壓力梯度之間的關(guān)系Fig.3 The flow rate versus the pressure gradient at the different maximum diameter of the 0th branching level

圖4顯示裂縫網(wǎng)絡(luò)無量綱滲透率隨孔隙率的變化趨勢曲線.在m＝10，n＝2，θ＝，μ＝1Pa·s，γ＝0.6，β＝和Df＝1.1時，裂縫網(wǎng)絡(luò)無量綱滲透率隨孔隙率的增加而逐漸增加，這也與物理事實相符合.

圖4 無量綱滲透率與孔隙率的關(guān)系曲線Fig.4 The dimensionless permeability versus the porosity

圖5 在不同的分形維數(shù)下，裂縫網(wǎng)絡(luò)無量綱滲透率隨長度比的變化趨勢圖Fig.5 The dimensionless permeability versus the length ratio

圖5顯示了裂縫網(wǎng)絡(luò)無量綱滲透率隨長度比的變化趨勢圖，其中參數(shù)為m＝10，n＝2，θ＝π／4，β＝從該圖可以看出，無量綱滲透率隨長度比的增加而減小，最后逐漸趨近于0，這是因為流體流動阻力隨著長度的增加而增加.

2 結(jié)論

本文基于分形理論與技術(shù)，將母管直徑具有分形分布的樹狀分叉網(wǎng)絡(luò)嵌入到母體材料中形成復(fù)合材料，在不考慮母體材料的滲透率情況下，得到了通過裂縫網(wǎng)絡(luò)的流量以及滲透率的分形表達(dá)式.結(jié)果表明通過裂縫網(wǎng)絡(luò)的流量與滲透率均對最大母管直徑很敏感.裂縫網(wǎng)絡(luò)滲透率k＋是最大母管直徑d0max、長度比γ以及孔隙率φ和分形維數(shù)Df的函數(shù)，不含任何經(jīng)驗常數(shù)；且滲透率隨分形維數(shù)Df和孔隙率φ的增加而增加，隨長度比的增加而減小.

［1］Mandelbrot B B.The Fractal Geometry of Nature［M］.New York：Freeman，1982.

［2］McCulloh K A，Sperry J S，Adler F R.Water transport in plants obeys Murray's law［J］.Nature，2003，421：939－942.

［3］Mathieu－Potvin F，Gosselin L.Optimal conduction pathways for cooling a heat－generating body：A comparison exercise［J］.International Journal of Heat and Mass Transfer，2007，50：2996－3006.

［4］Lasance C，Simons R E.Advances in high－performance cooling for electronics［J］.Electronics Cooling，2005，11：22－38.

［5］Frey U，Hislop A，Silverman M.Branching properties of the pulmonary arterial tree during pre－and postnatal development［J］.Respiratory Physiology ＆ Neurobiology，2004，139：179－189.

［6］Murray C D.The physiological principle of minimum work.I.The vascular system and the cost of blood volume［J］.Proc Natl Acad Sci U S A，1926，12：207－214.

［7］Murray C D.The physiological principle of minimum work applied to the angle of branching of arteries［J］.J Gen Physiol，1926，9：835－841.

［8］Chen Y P，Chen B Z，Shi M H，et al.Thermal and hydrodynamic characteristics of constructor tree－shaped minichannel heat Sink［J］.AIChE，2010，56：2018－2029.

［9］Xu P，Yu B M，F(xiàn)eng Y J，et al.Analysis of permeability for the fractal－like tree network by parallel and series models［J］.Physica A，2006，369：884－894.

［10］Xu P，Yu B M，F(xiàn)eng Y J，et al.Permeability of the fractal disk－shaped branched network with tortuosity effect［J］.Phys Fluids，2006，18：078103－078104.

［11］Wang Shifang，Yu Boming.Analysis of seepage for powerlaw fluids in the fractal－like tree network［J］.Transp Porous Media，2011，85：191－206.

［12］Li Yan，Yu Boming.Study of the starting pressure gradient in branching network［J］.Sci China Tech Sci，2010，53（9）：2397－2403.

［13］Lorente S，Wechsatol W，Bejan A.Tree－shaped flow structures designed by minimizing path lengths［J］.Int J Heat Mass Transfer，2002，45：3299－3312.

［14］Meijuan Yun，Boming Yu，Jianchao Cai.A fractal model for the starting pressure gradient for Bingham fluids in porous media ［J］.International Journal of Heat and Mass Transfer，2008，51：1402－1408.

［15］Yu Boming.Fractal analysis of permeabilities for porous media［J］.AI Ch E J，2004，50：46－57.

Analysis on seepage of the randomly distributed fractal－like tree network

WANG Shifang1，WU Tao2，DENG Yongju1，DAI Wei1
（1.School Physics and Electronic Information，Hubei University of Education，Wuhan 430205；2.School of Science，Wuhan Institute of Technology，Wuhan 430074）

Based on Darcy's law and Hagen－Poiseuille equation，the flow characteristics of the randomly distributed fractal－like tree network with the diameter d0of the 0th branching level are investigated by use of fractal theory，which follows the fractal scaling law.An analytical expression for permeability of the randomly distributed fractal－like tree network is presented.The results show that the flow rate for the randomly distributed fractal－like tree network scales linearly with pressure gradient.The dimensionless permeability for the randomly distributed fractal－like tree network increases with the increase of fractal dimension and porosity and decreases with the increase of the length ratio.

fractal－like tree network；permeability；Darcy law；pressure gradient

O357

A

1000－1190（2012）04－0406－04

2011－12－05.

湖北省教育廳科學(xué)技術(shù)項目（B20113002）.

＊E－mail：flatime＠sina.com.