趙靜 古斌

摘要：為研究多重流動(dòng)機(jī)制各向異性對(duì)頁(yè)巖氣運(yùn)移的影響，基于Biot線彈性孔隙介質(zhì)模型，考慮黏性流動(dòng)、氣體流動(dòng)動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散等流動(dòng)機(jī)制的各向異性，采用各向同性應(yīng)力依賴模型，建立流-固耦合的視滲透率模型，對(duì)應(yīng)力約束邊界條件下的頁(yè)巖氣運(yùn)移過(guò)程進(jìn)行數(shù)值模擬，定量分析流動(dòng)各向異性對(duì)頁(yè)巖氣運(yùn)移的影響。結(jié)果表明各種流動(dòng)機(jī)制的各向異性均對(duì)頁(yè)巖氣運(yùn)移有顯著影響：固有滲透率各向異性可明顯改變頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和流通量大小，其影響隨頁(yè)巖儲(chǔ)層滲透率增加而放大；當(dāng)固有滲透率較小時(shí)，動(dòng)態(tài)效應(yīng)各向異性會(huì)削弱氣體流動(dòng)的動(dòng)態(tài)效應(yīng)，對(duì)頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和流通量的影響十分顯著，但在固有滲透率較大時(shí)其影響可以忽略；表面吸附擴(kuò)散各向異性的影響與固有滲透率各向異性相似。準(zhǔn)確評(píng)估頁(yè)巖氣的運(yùn)移能力和產(chǎn)能應(yīng)考慮運(yùn)移過(guò)程中的流動(dòng)各向異性。

關(guān)鍵詞：頁(yè)巖氣運(yùn)移多重流動(dòng)機(jī)制流動(dòng)各向異性視滲透率模型數(shù)值模擬

中圖分類號(hào)：TE3文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1671-8755（2023）01-0040-07

Abstract：BasedonBiotsmodeloflinearelasticporousmedium，takingintoaccountfortheanisotropyofflowmechanismssuchasviscousflow，dynamiceffectsofgasflowandsurfaceadsorptivediffusion，andadoptingtheisotropicstressdependencemodel，afluid-solidcouplingapparentpermeabilitymodelwasestablishedtoinvestigatetheeffectofanisotropyofmultipleflowmechanismsonshalegastransport.Numericalsimulationontheshalegasmigrationprocessunderthestressconstraintboundaryconditionswasthenconductedtoquantitativelyanalyzetheinfluenceofflowanisotropyonshalegastransfer.Theresultsshowthattheanisotropyofeachflowmechanismhassignificantimpactonshalegastransport.Theanisotropyofintrinsicpermeabilitycangiverisetoremarkablechangesonthebeginningofshalegastransferandthegasflux，andsuchinfluenceisenhancedbyincreasingtheshalepermeability.Theanisotropyofdynamiceffectscandiminishthedynamiceffectsofgasflowandhasasignificantinfluenceonthebeginningofshalegastransferandthegasfluxwhentheintrinsicpermeabilityislow.Whiletheinfluencecanbeneglectedinthecaseoflargeintrinsicpermeability.Theeffectofanisotropyofsurfaceadsorptivediffusionissimilartothatofintrinsicpermeabilityanisotropy.Therefore，itisnecessarytoconsidertheflowanisotropyinthemigrationprocesswhenaccuratelyevaluatingthemigrationabilityandproductivityofshalegas.

Keywords：Shalegasmigration;Multipleflowmechanisms;Flowanisotropy;Apparentpermeabilitymodel;Numericalsimulation

頁(yè)巖儲(chǔ)層構(gòu)造復(fù)雜、孔隙類型和大小多樣、承受的地應(yīng)力高、呈現(xiàn)超低孔隙率和滲透率、各向異性和非均質(zhì)等特性。頁(yè)巖氣是自生自儲(chǔ)在頁(yè)巖中的非常規(guī)天然氣，其在儲(chǔ)層中的傳輸是一個(gè)多尺度和多場(chǎng)耦合的復(fù)雜過(guò)程［1-2］。目前，利用連續(xù)介質(zhì)力學(xué)理論建立視滲透率模型是研究頁(yè)巖氣運(yùn)移規(guī)律的主要手段之一。現(xiàn)有的視滲透率模型主要通過(guò)不同程度地考慮黏性流動(dòng)、氣體分子與孔壁碰撞引起的氣體流動(dòng)動(dòng)態(tài)效應(yīng)、表面吸附擴(kuò)散等多種流動(dòng)機(jī)制和固體變形（即應(yīng)力依賴）的綜合影響，對(duì)頁(yè)巖氣運(yùn)移的機(jī)制進(jìn)行闡述和分析［3-8］。研究發(fā)現(xiàn)頁(yè)巖氣運(yùn)移能力的改變?nèi)Q于3個(gè)因素的相互競(jìng)爭(zhēng)：表征黏性流動(dòng)的固有滲透率、動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散。本課題組的前期研究［7］建立了一個(gè)流-固耦合的視滲透率模型，全面考慮固體變形和上述多種流動(dòng)機(jī)制，分析了常應(yīng)力假設(shè)、初始變形狀態(tài)以及自由氣與吸附氣權(quán)重對(duì)頁(yè)巖氣運(yùn)移能力的影響。在此基礎(chǔ)上，進(jìn)一步討論了吸附效應(yīng)引起的孔隙體積模量改變以及該變化對(duì)頁(yè)巖氣運(yùn)移的影響［8］。但該模型與其他大多數(shù)視滲透率模型類似，假設(shè)與頁(yè)巖運(yùn)移相關(guān)的力學(xué)性質(zhì)均為各向同性。

眾多的實(shí)驗(yàn)研究表明，頁(yè)巖的力學(xué)性質(zhì)和頁(yè)巖氣流動(dòng)，如滲透率、彈性模量、氣體擴(kuò)散等，都存在明顯的各向異性。對(duì)不同的頁(yè)巖，平行和垂直于層理方向的滲透率比值可以達(dá)到幾至幾百［9-10］；對(duì)近似橫觀各向同性的頁(yè)巖，平行和垂直于各向同性面方向的彈性模量比值最大可達(dá)2.7［11］；氣體擴(kuò)散系數(shù)在平行和垂直于層理方向上表現(xiàn)出很大的差別［12-13］。數(shù)值模擬結(jié)果也揭示氣體滲透率和擴(kuò)散系數(shù)具有顯著的各向異性［14-15］。現(xiàn)有的視滲透率模型中只考慮了黏性流動(dòng)的各向異性，即采用各向異性固有滲透率。大多數(shù)研究關(guān)注由固體變形引起的固有滲透率各向異性?；凇盎鸩癜簟钡壤硐牖碾x散物理介質(zhì)模型，通過(guò)建立孔隙率和固有滲透率與應(yīng)變的關(guān)系，多種適用于裂縫型孔隙的各向異性應(yīng)力依賴模型已被提出和應(yīng)用［16-21］。在這些應(yīng)力依賴模型中，儲(chǔ)層、基體和孔隙可處理為具有各向同性或各向異性彈性參數(shù)的介質(zhì)，但這些研究極少關(guān)注頁(yè)巖氣運(yùn)移中的多重流動(dòng)機(jī)制。也有學(xué)者［22-23］利用不同的滲透率主值和孔隙取向來(lái)描述固有滲透率各向異性，結(jié)合經(jīng)驗(yàn)性的各向同性應(yīng)力依賴模型［22］或忽略應(yīng)力依賴性［23］，對(duì)頁(yè)巖氣的運(yùn)移過(guò)程進(jìn)行分析。盡管現(xiàn)有的視滲透率模型存在一定的不足，相應(yīng)的研究結(jié)果均發(fā)現(xiàn)流動(dòng)各向異性對(duì)非常規(guī)天然氣的產(chǎn)氣速率和產(chǎn)量有不可忽略的影響。鑒于固有滲透率各向異性的研究存在不足，頁(yè)巖氣運(yùn)移中動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散各向異性影響的研究至今尚未見(jiàn)報(bào)道，本工作在前期研究的基礎(chǔ)上，利用固有滲透率、Knudsen數(shù)和表面吸附擴(kuò)散系數(shù)的各向異性分別表征黏性流動(dòng)、動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散等流動(dòng)機(jī)制的各向異性特性，建立流-固耦合的視滲透率模型，對(duì)應(yīng)力約束邊界條件下的頁(yè)巖氣運(yùn)移過(guò)程進(jìn)行數(shù)值模擬，定量分析不同流動(dòng)各向異性對(duì)產(chǎn)能的影響。

1視滲透率模型和氣體運(yùn)移控制方程

假設(shè)平面應(yīng)變狀態(tài)下的頁(yè)巖氣儲(chǔ)層流動(dòng)相關(guān)的性質(zhì)為各向異性，固體變形相關(guān)的性質(zhì)為各向同性?？紤]黏性流動(dòng)、動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散3種流動(dòng)機(jī)制和固體變形的影響，前期建立的各向同性流-固耦合視滲透率模型［7］可擴(kuò)展為以下各向異性形式：

數(shù)值計(jì)算中輸入?yún)?shù)的取值見(jiàn)表1。為了表征固有滲透率和表面擴(kuò)散系數(shù)主值的各向異性，引入?yún)?shù)tk=kx∞0/ky∞0，tad=Dxad0/Dyad0。當(dāng)tk=tad=1時(shí)，退化為各向同性情況。為了討論取向角的影響，本文考慮所有取向角相同的狀態(tài)（固有滲透率和表面擴(kuò)散系數(shù)的取向角可不一樣），即Pk（θ）=δ（θ-θ0），Pad（θ）=δ（θ-θ′0）。

2.1固有滲透率各向異性的影響

為了研究固有滲透率各向異性對(duì)流通量的獨(dú)立影響，假設(shè)吸附氣流動(dòng)相關(guān)的參數(shù)保持不變，相應(yīng)的表面擴(kuò)散系數(shù)取值為：Dxad0=10-6m2/s，tad=2，Pad（θ）=δ（θ-π/4）。同時(shí)設(shè)動(dòng)態(tài)效應(yīng)為各向異性，為了降低動(dòng)態(tài)效應(yīng)的影響，固有滲透率主值取較大值。圖2和圖3分別表示當(dāng)kx∞0=5×10-19m2，tk=0.5時(shí)，頁(yè)巖氣運(yùn)移過(guò)程中x方向和y方向的流通量在不同固有滲透率取向角情況下隨時(shí)間的變化規(guī)律。

此外，圖中還給出了當(dāng)取向角均勻分布（即Pk（θ）=1/π）時(shí)流通量的變化。由圖可知，頁(yè)巖氣的運(yùn)移主要發(fā)生在壓力變化的某段時(shí)間內(nèi)，流通量的大小隨

時(shí)間先增大至峰值，然后逐漸減小。運(yùn)移開(kāi)始時(shí)間隨取向角變化，而且流通量的峰值大小也與取向角明顯相關(guān)。在所研究的例子中，x方向的流通量在取向角為0時(shí)最大，取向角為π/2時(shí)最小。y方向的流通量在取向角為0時(shí)最小，取向角為-π/4時(shí)最大。同時(shí)注意到，隨著取向角的變化，y方向的流通量發(fā)生流入和流出孔隙的轉(zhuǎn)變。由于模型中假設(shè)右邊界為井筒壓力作用邊界，因此x方向?yàn)轫?yè)巖氣的主要運(yùn)移方向，y方向流通量的大小預(yù)期要遠(yuǎn)小于x方向流通量，這與數(shù)值結(jié)果相符。在后續(xù)討論中，我們重點(diǎn)討論主要運(yùn)移方向上頁(yè)巖氣的產(chǎn)能，即x方向的流通量。圖4表示當(dāng)kx∞0=5×10-19m2，Pk（θ）=δ（θ-π/6）時(shí)，不同初始固有滲透率主值比對(duì)x方向流通量的影響。隨著主值比的增加，頁(yè)巖氣運(yùn)移開(kāi)始越早，x方向流通量峰值越高。而且與各向同性情況（tk=1）相比，頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和x方向流通量大小的差異都隨著初始固有滲透率主值各向異性程度增加而增加。綜合圖2和圖4可知，固有滲透率各向異性對(duì)頁(yè)巖氣運(yùn)移的產(chǎn)能有著重要的影響，而且各向異性的影響隨著整體固有滲透率的增大而放大，如圖5所示。圖5給出了4種不同的固有滲透率各向異性程度當(dāng)x方向初始固有滲透率主值由kx∞0=10-14m2增加為kx∞0=10-18m2時(shí)流通量的變化。結(jié)果顯示實(shí)心標(biāo)記曲線間的差異（kx∞0=10-19m2）遠(yuǎn)小于空心標(biāo)記曲線間（kx∞0=10-18m2）的差異。

2.2表面吸附擴(kuò)散各向異性的影響

為了研究表面吸附擴(kuò)散各向異性對(duì)流通量的影響，假設(shè)動(dòng)態(tài)效應(yīng)為各向異性且自由氣流動(dòng)相關(guān)的參數(shù)保持不變，相應(yīng)的固有滲透率等參數(shù)取值為

kx∞0=10-19m2，tk=0.5，Pk（θ）=δ（θ-π/6）。圖6和圖7分別給出了表面擴(kuò)散系數(shù)取向角（Dxad0=6×10-6m2/s，tad=2保持不變）和表面擴(kuò)散系數(shù)主值比（Dxad0=6×10-6m2/s，Pad（θ）=δ（θ-π/4）保持不變）對(duì)x方向流通量的影響規(guī)律。由圖可知，頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和流通量的大小均受表面吸附擴(kuò)散各向異性的影響。由于此時(shí)y方向的表面擴(kuò)散系數(shù)主值大于x方向的主值，表面擴(kuò)散系數(shù)取向角對(duì)流通量的影響規(guī)律與固有滲透率各向異性情況相反，如流通量峰值在取向角為0時(shí)最小，取向角為π/2時(shí)最大。當(dāng)tad<1時(shí)，兩個(gè)參數(shù)取向角的影響規(guī)律相似（相應(yīng)圖形省略）。表面擴(kuò)散系數(shù)主值比的影響規(guī)律與固有滲透率各向異性情況相似，如主值比越高，流通量越大。比較圖6與圖2或圖7與圖4可以發(fā)現(xiàn)，表面吸附擴(kuò)散各向異性引起的流通

量變化要小于固有滲透率各向異性，這是因?yàn)樵谟?jì)算中參數(shù)取值條件下，吸附氣流動(dòng)的貢獻(xiàn)要小于自由氣流動(dòng)。此外，頁(yè)巖儲(chǔ)層表面擴(kuò)散系數(shù)的增大也將放大表面吸附擴(kuò)散各向異性的影響（相關(guān)圖形略）。

2.3動(dòng)態(tài)效應(yīng)各向異性的影響

通過(guò)比較動(dòng)態(tài)效應(yīng)為各向異性和各向同性時(shí)x方向流通量的差別可分析動(dòng)態(tài)效應(yīng)各向異性的影響。由于表征動(dòng)態(tài)效應(yīng)的Knudsen數(shù)與固有滲透率大小相關(guān)，固有滲透率越小，動(dòng)態(tài)效應(yīng)越強(qiáng)。本文保持x方向初始固有滲透率主值和取向角不變（kx∞0=10-19m2，Pk（θ）=δ（θ-π/6），調(diào)整y方向初始固有滲透率主值的大小，考慮固有滲透率主值比tk=0.01，0.5，5等3種情況。此外，在數(shù)值計(jì)算中以下參數(shù)取值不變：Dxad0=6×10-6m2/s，tad=2，Pad（θ）=δ（θ-π/4）。圖8顯示了對(duì)不同程度的固有滲透率各向異性、動(dòng)態(tài)效應(yīng)分別處理為各向異性和各向同性時(shí)x方向流通量在頁(yè)巖氣運(yùn)移過(guò)程中的演化。從圖8可以看出，當(dāng)固有滲透率主值比較大（tk=5）時(shí)，動(dòng)態(tài)效應(yīng)各向異性的結(jié)果與各向同性時(shí)幾乎一致，表明此時(shí)動(dòng)態(tài)效應(yīng)各向異性的影響可以忽略。但是當(dāng)固有滲透率主值比較?。╰k=0.01）時(shí)，動(dòng)態(tài)效應(yīng)各向異性時(shí)的流通量遠(yuǎn)小于動(dòng)態(tài)效應(yīng)各向同性的情況，而且頁(yè)巖氣運(yùn)移也更晚開(kāi)始。同時(shí)注意到，當(dāng)動(dòng)態(tài)效應(yīng)為各向同性時(shí)，流通量峰值大小隨著tk取值的增加先減小后增加，這是由于動(dòng)態(tài)效應(yīng)與固有滲透率的影響相互競(jìng)爭(zhēng)的結(jié)果。當(dāng)動(dòng)態(tài)效應(yīng)為各向異性時(shí)，流通量峰值大小隨著tk取值的增加而單調(diào)增加，說(shuō)明固有滲透率較小時(shí)的動(dòng)態(tài)效應(yīng)受到抑制，從而改變了固有滲透率主值比的影響規(guī)律。以上結(jié)果表明動(dòng)態(tài)效應(yīng)各向異性對(duì)頁(yè)巖氣運(yùn)移產(chǎn)能有明顯影響。

頁(yè)巖氣氣藏產(chǎn)能的現(xiàn)場(chǎng)數(shù)據(jù)由于多種原因很難獲取，而文獻(xiàn)中也只提供儲(chǔ)層的各向異性力學(xué)性質(zhì)參數(shù)［9-13］，同時(shí)基于各向同性視滲透率模型的研究未能提供相應(yīng)的儲(chǔ)層各向異性參數(shù)，因此本工作未通過(guò)實(shí)例驗(yàn)證提出的考慮多重流動(dòng)機(jī)制各向異性的視滲透率模型的有效性。但本文建立的視滲透率模型可退化為實(shí)例驗(yàn)證的各向同性模型［7］，而且分析中采用的各向異性參數(shù)均依據(jù)頁(yè)巖儲(chǔ)層的實(shí)驗(yàn)數(shù)據(jù)，所以本工作中的理論模型和定量分析結(jié)果具有合理性和有效性。

3結(jié)論

（1）基于Biot線彈性孔隙連續(xù)介質(zhì)模型，利用張量形式的固有滲透率、Knudsen數(shù)和表面擴(kuò)散系數(shù)分別表征黏性流動(dòng)、動(dòng)態(tài)效應(yīng)和表面吸附擴(kuò)散3種流動(dòng)機(jī)制的各向異性，同時(shí)考慮固有滲透率的應(yīng)力依賴性，發(fā)展了流-固耦合的視滲透率模型。利用數(shù)值方法，該新模型可對(duì)頁(yè)巖氣運(yùn)移過(guò)程中的運(yùn)移能力和產(chǎn)能進(jìn)行分析。

（2）固有滲透率各向異性可明顯改變頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和流通量大小，其影響隨頁(yè)巖儲(chǔ)層滲透率增加而放大。表面吸附擴(kuò)散各向異性的影響與固有滲透率各向異性相似。

（3）當(dāng)固有滲透率較小時(shí)，動(dòng)態(tài)效應(yīng)各向異性會(huì)削弱氣體流動(dòng)的動(dòng)態(tài)效應(yīng)，對(duì)頁(yè)巖氣運(yùn)移開(kāi)始時(shí)間和流通量的影響十分顯著，但在固有滲透率較大時(shí)，其影響可以忽略。

參考文獻(xiàn)

［1］CURTISJB.Fracturedshalegassystems［J］.AAPGBulletin，2002，86（11）：1921-1938.

［2］GENSTERBLUMY，GHANIZADEHA，ROBERTCJ，etal.Gastransportandstoragecapacityinshalegasreservoirs—Areview.PartA：Transportprocesses［J］.JournalofUnconventionalOilandGasResources，2015，12：87-122.

［3］BESKOKA，KARNIADAKISGE.Amodelforflowsinchannels，pipes，andductsatmicroandnanoscales［J］.MicroscaleThermophysicalEngineering，1999，3（1）：43-77.

［4］JAVADPOURF.Nanoporesandapparentpermeabilityofgasflowinmudrocks（shalesandsiltstone）［J］.JournalofCanadianPetroleumTechnology，2009，48（8）：16-21.

［5］WUKL，LIXF，GUOCH，etal.Aunifiedmodelforgastransferinnanoporesofshalegasreservoirs：couplingporediffusionandsurfacediffusion［J］.SPEJournal，2016，21（5）：1583-1611.

［6］CAOP，LIUJS，LEONGY-K.Combinedimpactofflowregimesandeffectivestressontheevolutionofshaleapparentpermeability［J］.JournalofUnconventionalOilandGasResources，2016，14：32-43.

［7］YANGX，ZHANGHB，WUWB，etal.Gasmigrationinthereservoirsofultralowporosityandpermeabilitybasedonanimprovedapparentpermeabilitymodel［J］.JournalofPetroleumScienceandEngineering，2020，185：106614.

［8］YINGQ，ZHANGHB，ZHAOJ，etal.InfluenceofconsideringthesorptioneffectintheBetti-Maxwellreciprocaltheoremongastransportcapacityinunconventionalreservoirs［J］.TransportinPorousMedia，2021，137（3）：451-469.

［9］BHANDARIAR，F(xiàn)LEMINGSPB，POLITOPJ，etal.AnisotropyandstressdependenceofpermeabilityintheBarnettshale［J］.TransportinPorousMedia，2015，108（2）：393-411.

［10］PANZJ，MAY，CONNELLLD，etal.Measuringanisotropicpermeabilityusingacubicshalesampleinatriaxialcell［J］.JournalofNaturalGasScience&Engineering，2015，26：336-344.

［11］CHOJW，KIMH，JEONS，etal.DeformationandstrengthanisotropyofAsangneiss，BoryeongshaleandYeoncheonschist［J］.InternationalJournalofRockMechanics&MiningSciences，2012，50：158-169.

［12］WANY，PANZJ，TANGSH，etal.AnexperimentalinvestigationofdiffusivityandporosityanisotropyofaChinesegasshale［J］.JournalofNaturalGasScienceandEngineering，2015，23：70-79.

［13］李赟.C7頁(yè)巖氣藏吸附擴(kuò)散微觀機(jī)理研究［D］.成都：西南石油大學(xué)，2017.

［14］YINY，QUZG，ZHANGT，etal.Threedimensionalporescalestudyofmethanegasmassdiffusioninshalewithspatiallyheterogeneousandanisotropicfeatures［J］.Fuel，2020，273：117750.

［15］李滔，李閩，荊雪琪，等.孔隙尺度各向異性與孔隙分布非均質(zhì)性對(duì)多孔介質(zhì)滲透率的影響機(jī)理［J］.石油勘探與開(kāi)發(fā)，2019，46（3）：569-579.

［16］WANGGX，MASSAROTTOP，RUDOLPHV.AnimprovedpermeabilitymodelofcoalforcoalbedmethanerecoveryandCO2geosequestration［J］.InternationalJournalofCoalGeology，2009，77（1/2）：127-136.

［17］PANZJ，CONNELLLD.Modellingofanisotropiccoalswellinganditsimpactonpermeabilitybehaviourforprimaryandenhancedcoalbedmethanerecovery［J］.InternationalJournalofCoalGeology，2011，85（3/4）：257-267.

［18］GUFG，CHALATURNYKR.Permeabilityandporositymodelsconsideringanisotropyanddiscontinuityofcoalbedsandapplicationincoupledsimulation［J］.JournalofPetroleumScienceandEngineering，2010，74：113-131.

［19］WUY，LIUJS，DEREKE，etal.Developmentofanisotropicpermeabilityduringcoalbedmethaneproduction［J］.JournalofNaturalGasScienceandEngineering，2010，2（4）：197-210.

［20］WANGJG，LIUJS，KABIRA.Combinedeffectsofdirectionalcompaction，nonDarcyflowandanisotropicswellingoncoalseamgasextraction［J］.InternationalJournalofCoalGeology，2013，109：1-14.

［21］WANGK，ZANGJ，WANGGD，etal.Anisotropicpermeabilityevolutionofcoalwitheffectivestressvariationandgassorption：Modeldevelopmentandanalysis［J］.InternationalJournalofCoalGeology，2014，130：53-65.

［22］HUANGT，TAOZW，LIE，etal.Effectofpermeabilityanisotropyontheproductionofmultiscaleshalegasreservoirs［J］.Energies，2017，10（10）：1549.

［23］單嫻，姚軍.基于滲透率張量的各向異性油藏兩相滲流數(shù)值模擬［J］.中國(guó)石油大學(xué)學(xué)報(bào)，2011，35（2）：101-106.

［24］AL-YOUSEFHY.Permeabilityanisotropymeasurementonwholecores—analyticalsolutionandapplication［C］∥14thSPEMiddleEastOil&GasShowandConference.Bahrain，2005.