高 英，朱維耀，岳 明，李?lèi)?ài)山，張燎原，宋洪慶

（1.北京科技大學(xué) 土木與環(huán)境工程學(xué)院，北京 100083； 2.中石化勝利油田分公司 采油工藝研究院，山東 東營(yíng)257000）

0 引言

世界石油工業(yè)正向非常規(guī)油氣跨越，最早認(rèn)知的非常規(guī)油氣為致密砂巖氣、深盆氣、煤層氣、頁(yè)巖氣等致密儲(chǔ)層天然氣.隨著頁(yè)巖氣理論研究的深入及勘探開(kāi)發(fā)的迅速發(fā)展，頁(yè)巖油也成為油氣開(kāi)發(fā)的另一新熱點(diǎn)，成為油氣供給的重要來(lái)源［1-6］.頁(yè)巖油不以浮力作用為聚集動(dòng)力，具有源儲(chǔ)一體、納米級(jí)孔喉、儲(chǔ)集物性致密、裂縫系統(tǒng)發(fā)育、儲(chǔ)集層脆性指數(shù)較高等典型非常規(guī)油氣特點(diǎn).頁(yè)巖油儲(chǔ)層屬于超低滲透致密儲(chǔ)層，一般無(wú)自然產(chǎn)能或低產(chǎn)，鄒才能等［7］認(rèn)為需要通過(guò)“人造滲透率”，即體積壓裂、重復(fù)壓裂技術(shù)手段改造儲(chǔ)層，實(shí)現(xiàn)頁(yè)巖油有效開(kāi)發(fā).

由于受頁(yè)巖油儲(chǔ)層的脆性特征及天然裂縫發(fā)育等影響，水力壓裂裂縫不再是傳統(tǒng)的雙翼對(duì)稱(chēng)裂縫容易形成錯(cuò)綜復(fù)雜的裂縫網(wǎng)絡(luò)系統(tǒng)［8］.人們對(duì)低滲透油藏、頁(yè)巖氣體積壓裂已有研究［9-15］，但對(duì)頁(yè)巖油體積壓裂井的滲流理論和數(shù)學(xué)模型研究尚處于起步階段，多數(shù)是采用數(shù)值模擬軟件研究低滲透油藏體積壓裂井產(chǎn)能［16-17］，并不能很好體現(xiàn)頁(yè)巖油儲(chǔ)層的滲流特性.筆者結(jié)合質(zhì)量和動(dòng)量守恒方程，根據(jù)頁(yè)巖油在井周?chē)黧w的不同流動(dòng)特征［18-19］，將裂縫周?chē)牧黧w流動(dòng)劃分為2個(gè)區(qū)域，建立體積壓裂成縫網(wǎng)的頁(yè)巖油儲(chǔ)層直井二區(qū)耦合滲流數(shù)學(xué)模型，推導(dǎo)解析解，得到頁(yè)巖油體積壓裂直井產(chǎn)能方程，能夠直觀評(píng)價(jià)體積壓裂縫網(wǎng)效果和預(yù)測(cè)直井產(chǎn)能.

1 納微米孔喉滲流特征

頁(yè)巖油儲(chǔ)集層中廣泛發(fā)育納米級(jí)孔喉［5，7］，孔徑主要為50～300nm，孔隙度一般小于10%，滲透率小于0.01×10-3μm2.儲(chǔ)層巖石的微觀孔喉結(jié)構(gòu)直接影響儲(chǔ)層的儲(chǔ)集滲流能力，并最終決定油氣藏產(chǎn)能分布的差異.利用巖心流動(dòng)實(shí)驗(yàn)，分別測(cè)定流體在不同流量條件下通過(guò)巖心時(shí)的壓力梯度（見(jiàn)圖1），并繪制啟動(dòng)壓力梯度與滲透率關(guān)系曲線（見(jiàn)圖2）.

由圖1可知，頁(yè)巖油開(kāi)發(fā)過(guò)程中，直線段的延伸與壓力梯度軸的交點(diǎn)不經(jīng)過(guò)坐標(biāo)原點(diǎn)，曲線具有明顯的非線性特征，流體滲流不再遵循達(dá)西定律，出現(xiàn)低速非達(dá)西滲流，具有擬啟動(dòng)壓力梯度.低滲巖心的啟動(dòng)壓力梯度隨流量的增加而增加，巖心的滲透率越小，啟動(dòng)壓力梯度越大.

由圖2可知，流體在頁(yè)巖油儲(chǔ)層中流動(dòng)時(shí)，啟動(dòng)壓力梯度與地層的滲透率密切相關(guān)，隨著滲透率的降低，地層流體的啟動(dòng)壓力梯度急劇增加，啟動(dòng)壓力梯度與滲透率之間呈冪函數(shù)關(guān)系，即

式中：G為啟動(dòng)壓力梯度；K為滲透率.

頁(yè)巖油儲(chǔ)層滲透率小于0.01×10-3μm2，啟動(dòng)壓力梯度達(dá)到5MPa／m，遠(yuǎn)遠(yuǎn)大于低滲透油藏的0.02 MPa／m的啟動(dòng)壓力梯度［20］，存在明顯的強(qiáng)非線性滲流特征.因此，頁(yè)巖油儲(chǔ)層納微米孔喉的滲流運(yùn)動(dòng)方程可用廣義公式表示為

式中：v為滲流速度；μ為黏度；▽p為壓力梯度.

2 表征模型

2.1 裂縫形態(tài)

頁(yè)巖油儲(chǔ)層中發(fā)育微裂縫，主要有2種類(lèi)型：一是發(fā)育在顆粒內(nèi)部；二是發(fā)育在碎屑顆粒邊緣.顆粒內(nèi)部的微裂縫一般比較平直，曲折度較小，少有膠結(jié)物充填.顆粒間的微裂縫呈鋸齒狀彎曲.微裂縫長(zhǎng)度為5.5～12.0μm，裂縫間距可達(dá)50nm以上，但延伸很短.存在微裂縫的區(qū)域，巖石脆性指數(shù)較高，易形成微裂縫網(wǎng)絡(luò)，從而成為頁(yè)巖中微觀尺度上油氣滲流的主要通道.

體積壓裂過(guò)程中，存在剪切、錯(cuò)斷和滑移等復(fù)雜的力學(xué)行為，在張性主裂縫存在的同時(shí)還存在大量的次生裂縫，并與微裂縫形成錯(cuò)綜復(fù)雜的網(wǎng)絡(luò)系統(tǒng)（見(jiàn)圖3），增大改造體積，溝通流場(chǎng).斷裂和裂縫的分布和幾何形體具有明顯的分形結(jié)構(gòu)，斷裂系統(tǒng)可以使用分維度、統(tǒng)計(jì)自相似性和冪函數(shù)等研究，以分維數(shù)度量［21-23］.因此，采用分形理論描述體積壓裂改造后裂縫的分布特征，并表征體積壓裂“人造滲透率”.

2.2 “人造滲透率”

為了研究滲流速度，建立理想模型：地層厚度為h、半徑為r、井徑為rw、井底流壓為pw、流量為Q的圓柱型徑向滲流區(qū)，見(jiàn)圖4.取dr的微元，設(shè)通過(guò)單位橫截面A0的裂縫總數(shù)為Mc，裂縫的最大寬度和最小寬度分別為λa和λi，裂縫寬度?大于λ的裂縫數(shù)量Mc滿(mǎn)足

式中：D為裂縫寬度分形維數(shù)，量綱為一；-d Mc＞0為單位截面A0中所通過(guò)的寬度在λ與λ＋dλ之間裂縫折合數(shù).

考慮實(shí)際地層中裂縫的彎曲特性，實(shí)際長(zhǎng)度Le應(yīng)遵循分形的基本公式，即

式中：δ為迂曲分維，量綱為一；L為裂縫外觀長(zhǎng)度.

根據(jù)裂縫面密度定義（觀測(cè)面積內(nèi)全部裂縫長(zhǎng)度的總和與觀測(cè)面積之比），裂縫寬度為λ的單條裂縫長(zhǎng)度為dLe，圖4中圓管面積為2πrh／A0，圓筒區(qū)域的裂縫面密度Df為

由流量定義可得到通過(guò)單位橫截面流量為單條裂縫流量Q1（λ）與裂縫條數(shù)的乘積，即

式中：η0為裂縫深度；V為單位橫截面的流量；p為壓力；dp／dLe為壓力梯度.

將式（4）代入式（6）中，整理可得

根據(jù)一般徑向流達(dá)西公式，得

對(duì)式（7-8）從（rw，pw）到（re，pe）積分，可以得到裂縫的有效滲透率 Kn為

式中：re為供給半徑；pe為供給壓力.

裂縫體積系數(shù)fn為圓管區(qū)域的裂縫體積Vf與圓管總體積Vr之比，即

則“人造滲透率”，即體積壓裂改造區(qū)的等效滲透率Ke為

式中：Km為基質(zhì)滲透率.

3 產(chǎn)能數(shù)學(xué)模型

假設(shè)壓裂改造體積為沿主裂縫對(duì)稱(chēng)的2×xf×b的橢圓柱體［9，13］，根據(jù)頁(yè)巖油儲(chǔ)層體積壓裂改造后的流動(dòng)特點(diǎn)，儲(chǔ)層流體滲流分為2個(gè)區(qū)域（見(jiàn)圖5）：一是體積壓裂主裂縫內(nèi)的線性達(dá)西流動(dòng)；二是體積壓裂改造形成的橢圓縫網(wǎng)區(qū)域的低速非達(dá)西流動(dòng).

3.1 基本假設(shè)

主要根據(jù)頁(yè)巖油儲(chǔ)層體積壓裂井生產(chǎn)時(shí)流體的流動(dòng)特征，考慮頁(yè)巖油儲(chǔ)層非線性滲流特征，建立頁(yè)巖油體積壓裂改造儲(chǔ)層直井產(chǎn)能預(yù)測(cè)模型.假設(shè)：（1）頁(yè)巖油儲(chǔ)層為上下封閉且無(wú)限大地層；（2）對(duì)直井進(jìn)行體積壓裂，儲(chǔ)層體積壓裂改造后形成橢圓形的縫網(wǎng)，橢圓形體積改造區(qū)域短半軸長(zhǎng)為b，焦距為主裂縫半長(zhǎng)，見(jiàn)圖5；（3）油藏和裂縫內(nèi)流體為單相流體，不可壓縮，滲流為等溫穩(wěn)定滲流，不考慮重力影響；（4）滲流過(guò)程中考慮啟動(dòng)壓力梯度的影響.

3.2 主裂縫的線性流動(dòng)

主干縫內(nèi)流體的流動(dòng)服從達(dá)西定律，屬于線性流動(dòng)，其運(yùn)動(dòng)方程為

其穩(wěn)態(tài)流動(dòng)表達(dá)式為

式中：Kf為主干縫滲透率.

對(duì)式（13）進(jìn)行分離變量，并從（xf，pxf）到（rw，pw）積分，可以得到裂縫內(nèi)流體的流量和壓差之間的關(guān)系表達(dá)式，即

式中：pxf為主裂縫兩端的壓力；wf為主裂縫寬度；xf為主裂縫半長(zhǎng).

3.3 改造區(qū)的橢圓滲流

改造后儲(chǔ)層橢圓滲流區(qū)直角坐標(biāo)系（x，y）和橢圓坐標(biāo)系（η，ξ）的關(guān)系為

式中：a、b分別為橢圓的長(zhǎng)軸和短軸.

對(duì)于改造后的頁(yè)巖油儲(chǔ)層，其等效滲透率Ke即“人造滲透率”依然不高，該區(qū)域流動(dòng)為非線性滲流，存在啟動(dòng)壓力梯度，因此橢圓區(qū)穩(wěn)態(tài)滲流的流量表達(dá)式為

式中：B為體積壓縮因子.

結(jié)合式（17），對(duì)式（16）從（ξw，pxf）到（ξi，pe）進(jìn)行積分，可以得到橢圓滲流區(qū)的流量與壓差公式為

式中：ξw為井筒附近橢圓坐標(biāo)；ξi為泄油區(qū)橢圓坐標(biāo).

3.4 二區(qū)耦合后產(chǎn)能解析式

根據(jù)兩區(qū)交界處壓力相等，對(duì)式（14）和式（18）聯(lián)立，可以得出頁(yè)巖油儲(chǔ)層體積壓裂改造后直井的產(chǎn)能解析解表達(dá)式，即

4 單井產(chǎn)能影響因素分析

以勝利油田某實(shí)驗(yàn)區(qū)塊一口頁(yè)巖油儲(chǔ)層體積壓裂直井為例，儲(chǔ)層厚度為61m，原始地層壓力為45.75 MPa，生產(chǎn)壓差為20MPa，原油黏度為6.82mPa·s，原油密度為816.1kg·m-3，基質(zhì)滲透率為0.001×10-3μm2，孔隙度為0.024，壓裂裂縫半長(zhǎng)為200m，縫寬為4mm，主裂縫滲透率為10μm2.初始分形參數(shù)：縫寬分形維數(shù)為1.1，迂曲分維為1.1.啟動(dòng)壓力梯度由式（1）求得.將地質(zhì)參數(shù)代入模型進(jìn)行計(jì)算并分析參數(shù)影響因素.模型計(jì)算生產(chǎn)井初期產(chǎn)能為4.0t·d-1，而實(shí)際產(chǎn)量為3.7t·d-1，誤差較小，說(shuō)明本模型可靠性強(qiáng).

4.1 體積壓裂與常規(guī)壓裂

體積壓裂和常規(guī)壓裂的直井產(chǎn)能見(jiàn)圖6.由圖6可知，由于頁(yè)巖油儲(chǔ)層具有致密性，常規(guī)壓裂后基本沒(méi)有產(chǎn)能，只能通過(guò)體積壓裂改造儲(chǔ)層溝通天然裂縫，形成縫網(wǎng)才能獲得較高產(chǎn)能.

4.2 縫寬分形維數(shù)與裂縫面密度、基質(zhì)—裂縫等效滲透率

不同迂曲分維時(shí)縫寬分形維數(shù)與裂縫面密度和基質(zhì)—裂縫等效滲透率的關(guān)系見(jiàn)圖7.由圖7（a）可知，裂縫的縫寬分維和迂曲分維對(duì)裂縫面密度有很大影響，相同迂曲分維下，裂縫的面密度越大，裂縫分形維數(shù)越大，當(dāng)迂曲分維為1.1，分形維數(shù)從1.0到1.8變化時(shí)，裂縫面密度從3m-1增長(zhǎng)到20m-1；迂曲分維的增加即裂縫彎曲程度的增加，變相增加裂縫長(zhǎng)度，也增加裂縫面密度.

由圖7（b）可知，裂縫的分形維數(shù)和迂曲分維對(duì)等效后的基質(zhì)—裂縫等效滲透率有很大影響，相同迂曲分維下，裂縫分形維數(shù)越大，裂縫條數(shù)越多，等效滲透率越大；迂曲分維的增加相當(dāng)于延長(zhǎng)流體流動(dòng)的距離，盡管裂縫面密度增加，但等效滲透率減小.

4.3 產(chǎn)能影響

4.3.1 裂縫面密度、體積壓裂寬度

不同裂縫面密度時(shí)體積壓裂改造寬度與直井產(chǎn)能的關(guān)系見(jiàn)圖8.由圖8可知，在相同的裂縫面密度條件下，體積壓裂改造儲(chǔ)層越寬，產(chǎn)能越高；當(dāng)裂縫面密度較小時(shí)，儲(chǔ)層改造寬度的增加對(duì)提升單井產(chǎn)能影響較小；當(dāng)裂縫面密度較大時(shí)，儲(chǔ)層改造寬度的增加對(duì)提升單井產(chǎn)能影響巨大.

4.3.2 裂縫導(dǎo)流能力

不同次生縫導(dǎo)流能力時(shí)裂縫導(dǎo)流能力對(duì)產(chǎn)能的影響見(jiàn)圖9.由圖9可知，在相同次生縫導(dǎo)流能力條件下，主干縫導(dǎo)流能力越強(qiáng)，產(chǎn)能越高.在次生縫導(dǎo)流能力較低的情況下，主干縫導(dǎo)流能力的增加對(duì)產(chǎn)能提升較??；在次生縫導(dǎo)流能力很強(qiáng)的情況下，主干縫導(dǎo)流能力的增加對(duì)產(chǎn)能影響很大.

5 結(jié)論

（1）頁(yè)巖油儲(chǔ)層流體滲流存在啟動(dòng)壓力梯度，具有強(qiáng)非線性特征，啟動(dòng)壓力梯度大小與滲透率之間呈冪函數(shù)關(guān)系.

（2）采用分形理論，推導(dǎo)出裂縫儲(chǔ)層等效滲透率、裂縫體積系數(shù)和裂縫面密度參數(shù)表達(dá)式，分形維數(shù)越大，裂縫面密度越大，等效滲透率越大.

（3）當(dāng)裂縫面密度較大時(shí)，增加儲(chǔ)層改造寬度對(duì)提升產(chǎn)能效果顯著，當(dāng)裂縫面密度較小時(shí)，增加儲(chǔ)層改造寬度對(duì)產(chǎn)能提升效果很??；當(dāng)次生縫導(dǎo)流能力較高時(shí)，增加主干縫導(dǎo)流能力對(duì)產(chǎn)能提升效果較強(qiáng).因此現(xiàn)場(chǎng)應(yīng)盡可能增加裂縫面密度和次生縫的導(dǎo)流能力提高產(chǎn)能.

（References）：

［1］周慶凡，楊國(guó)豐.致密油與頁(yè)巖油的概念與應(yīng)用［J］.石油與天然氣地質(zhì)，2012，33（4）：541-544.Zhou Qingfan，Yang Guofeng.Definition and application of tight oil and shale oil terms［J］.Oil ＆ Gas Geology，2012，33（4）：541-544.

［2］車(chē)長(zhǎng)波，朱杰，李富兵，等.全球油氣資源形勢(shì)［J］.天然氣工業(yè)，2010，30（1）：1-4.Che Changbo，Zhu Jie，Li Fubing，et al.The status quo and developing trend of global hydrocarbon resources［J］.Natural Gas Industry，2010，30（1）：1-4.

［3］楊華，李士祥，劉顯陽(yáng).鄂爾多斯盆地致密油、頁(yè)巖油特征及資源潛力［J］.石油學(xué)報(bào)，2013，34（1）：1-11.Yang Hua，Li Shixiang，Liu Xianyang.Characteristics and resource prospects of tight oil and shale oil in Ordos basin［J］.Acta Petrolei Sinica，2013，34（1）：1-11.

［4］Bowker K A.Barnett shale gas production，fort worth basin：Issues and discussion［J］.AAPG Bulletin，2007，91（4）：523-533.

［5］Gale J F，Reed R M，Holder J，et al.Natural fracture in the barnett shale and their importance for hydraulic fracture treatments［J］.AAPG Bulletin，2007，91（4）：603-622.

［6］久凱，丁文龍，黃文輝，等.渤海灣盆地濟(jì)陽(yáng)拗陷沙三段頁(yè)巖氣地質(zhì)條件分析［J］.大慶石油學(xué)院學(xué)報(bào)，2012，36（2）：65-70.Jiu Kai，Ding Wenlong，Huang Wenhui，et al.Analysis of geological condition of the formation of Shahejie shale gas in Jiyang depression in Bohai bay basin［J］.Journal of Daqing petroleum institute，2012，36（2）：65-70.

［7］鄒才能，楊智，崔景偉，等.頁(yè)巖油形成機(jī)制、地質(zhì)特征及發(fā)展對(duì)策［J］.石油勘探與開(kāi)發(fā)，2013，40（1）：14-26.Zou Caineng，Yang Zhi，Cui Jingwei，et al.Formation mechanism，geological characteristics and development strategy of nonmarine shale oil in China［J］.Petroleum Exploration and Development，2013，40（1）：14-26.

［8］趙金洲，任嵐，胡永全.頁(yè)巖油儲(chǔ)層壓裂縫成網(wǎng)延伸的受控因素分析［J］.西南石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2013，35（1）：2-6.Zhao Jinzhou，Ren Lan，Hu Yongquan.Controlling factors of hydraulic fractures extending into network in shale formations［J］.Journal of Southwest Petroleum University：Science ＆ Technology Edition，2013，35（1）：2-6.

［9］雷群，胥云，蔣廷學(xué)，等.用于提高低—特低滲透油氣藏改造效果的縫網(wǎng)壓裂技術(shù)［J］.石油學(xué)報(bào)，2009，30（2）：237-241.Lei Qun，Xu Yun，Jiang Tingxue，et al.'Fracture network'fracturing technique for improving post-fracturing performance of low and ultra-low permeability reservoirs［J］.Acta Petrolei Sinica，2009，30（2）：237-241.

［10］陳作，薛承瑾，蔣廷學(xué)，等.頁(yè)巖氣井體積壓裂技術(shù)在我國(guó)的應(yīng)用建議［J］.天然氣工業(yè)，2010，30（10）：30-32.Chen Zuo，Xue Chengjin，Jiang Tingxue，et al.Proposals for the application of fracturing by stimulated reservoir volume（SRV）in shale gas wells in China［J］.Natural Gas Industry，2010，30（10）：30-32.

［11］鄧佳，朱維耀，劉錦霞，等.考慮應(yīng)力敏感性的頁(yè)巖氣產(chǎn)能預(yù)測(cè)模型［J］.天然氣地球科學(xué)，2013，24（3）：456-460.Deng Jia，Zhu Weiyao，Liu Jinxia，et al.Productivity prediction model of shale gas considering stress sensitivity［J］.Natural Gas Geoscience，2013，24（3）：456.

［12］程遠(yuǎn)方，李友志，時(shí)賢，等.頁(yè)巖氣體積壓裂縫網(wǎng)模型分析及應(yīng)用［J］.天然氣工業(yè)，2013，33（9）：53-59.Cheng Yuanfang，Li Youzhi，Shi Xian，et al.Analysis and application of fracture network models of volume fracturing in shale gas reservoirs［J］.Natural Gas Industry，2013，33（9）：53-59.

［13］劉雄，田昌炳，姜龍燕，等.致密油藏直井體積壓裂穩(wěn)態(tài)產(chǎn)能評(píng)價(jià)模型［J］.東北石油大學(xué)學(xué)報(bào)，2014，38（1）：91-95.Liu Xiong，Tian Changbing，Jiang Longyan，et al.Steady deliverability evaluation model of fracture network reconstruction vertical well in tight oil reservoir［J］.Journal of northeast petroleum university，2014，38（1）：91-95.

［14］Zheng Q，Xu J，Yang B，et al.Research on the effective gas diffusion coefficient in dry porous media embedded with a fractal-like tree network［J］.Physica A，2013（392）：1557-1566.

［15］Zhu W Y，Yue M，Ma D X，et al.The micro-structure and seepage characteristics of shale reservoir［J］.Journal of Chemical and Pharmaceutical Research，2014，6（1）：312-315.

［16］Tran T，Sinurat P，Wattenbarger R A.Production characteristics of the bakken shale oil［R］.SPE145684，2011.

［17］Giger F M.Horizontal wells production techniques in heterogeneous reservoirs［R］.SPE13710，1988，307.

［18］朱維耀，岳明，宋洪慶，等.考慮解吸擴(kuò)散過(guò)程的煤層氣壓裂直井產(chǎn)能模型研究［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2013，21（5）：951-960.Zhu Weiyao，Yue Ming，Song Hongqing，et al.Productivity model of gas flow in CBM fractured reservoirs considering desorption and diffusion［J］.Journal of Basic Science and Engineering，2013，21（5）：951-960.

［19］王志平，朱維耀，岳明，等.低、特低滲透油藏壓裂水平井產(chǎn)能計(jì)算方法［J］.北京科技大學(xué)學(xué)報(bào)，2012，34（7）：750-754.Wang Zhiping，Zhu Weiyao，Yue Ming，et al.A method to predict the production of fractured horizontal wells in low／ultralow permeability reservoirs［J］.Journal of University of Science and Technology Beijing，2012，34（7）：750-754.

［20］郝斐，程林松，李春蘭，等.特低滲透油藏啟動(dòng)壓力梯度研究［J］.西南石油學(xué)院學(xué)報(bào)，2006，28（6）：29-32.Hao Fei，Cheng Linsong，Li Chunlan，et al，Study on threshold pressure gradient in ultea-low permeability reservoir［J］.Journal of Southwest Petroleum Institute，2006，28（6）：29-32.

［21］馮陣東，戴俊生，鄧航，等.利用分形幾何定量評(píng)價(jià)克拉2氣田裂縫［J］.石油與天然氣地質(zhì)，2011，32（54）：928-933.Feng Zhendong，Dai Junsheng，Deng Hang，et al.Quantitative evaluation of fractures with fractal geometry in Kela-2gas field［J］.Oil ＆ Gas Geology，2011，32（54）：928-933.

［22］Yu B M，Cheng P.A fractal permeability model for bi-dispersed porous media［J］.International Journal of Heat and Mass Transfer，2002（45）：2983-2993.

［23］張吉昌，田國(guó)清，劉建中.儲(chǔ)層構(gòu)造裂縫的分形分析［J］.石油勘探與開(kāi)發(fā)，1996，23（4）：65.Zhang Jichang，Tian Guoqing，Liu Jianzhong.A fractal analysis on structural fractures of reservoirs［J］.Petroleum Exploration and Development，1996，23（4）：65.