摘要：在建立近場動力學非局部多孔介質(zhì)滲透模型的基礎(chǔ)上，引入多個可反映非局部作用程度的核函數(shù)以提高計算收斂精度，推導(dǎo)不同核函數(shù)情況下對應(yīng)的近場動力學滲透系數(shù)，并在二維滲透模型中建立服從威布爾分布的滲透系數(shù)及裂隙網(wǎng)絡(luò)區(qū)域滲流來分別實現(xiàn)多孔介質(zhì)基質(zhì)及裂隙非均勻滲流的模型，彌補了經(jīng)典近場動力學模型無法較好的模擬巖石、土壤等多孔介質(zhì)非均勻滲流的不足。運用不同的核函數(shù)模型對一維滲流問題模擬，分析核函數(shù)對結(jié)果的影響，結(jié)果顯示改進后的模型可以較好的收斂于理論解，并且多項式型核函數(shù)相比其他核函數(shù)收斂精度最高。隨后在二維模型中引入多項式型核函數(shù)并推導(dǎo)對應(yīng)二維滲透系數(shù)，利用本文提出的非均質(zhì)滲流方法對二維不含裂隙及含裂隙滲流過程進行模擬，模擬結(jié)果表明所提出的新模型能夠較好的模擬巖石材料非均勻滲流過程，在多孔介質(zhì)滲流數(shù)值模擬方面有較廣泛的應(yīng)用前景。

　　關(guān)鍵詞：多孔介質(zhì);近場動力學;非局部方法;滲流模擬

　　土壤、巖體等多孔介質(zhì)中賦存大量天然孔隙及構(gòu)造裂隙，為地下水的傳輸、儲存提供了場所[13]。近年來隨著高邊坡及地下隧道等工程增多，由地下水引發(fā)的工程災(zāi)害頻發(fā)，造成重大的經(jīng)濟損失。地下水在多孔介質(zhì)中的滲流特性是土力學、巖石力學等學科中關(guān)注的熱點。多孔介質(zhì)內(nèi)無法定量化的孔隙及裂隙會對滲流特性產(chǎn)生較大影響，使得采取理論手段探究滲流問題時十分困難;并且利用實驗方法探究時存在周期長、成本高及實參數(shù)不易確定等問題。隨著計算機技術(shù)的發(fā)展，越來越多的學者采用數(shù)值模擬的手段探究多孔介質(zhì)滲流特性。

　　力學工程評職知識： 工程力學期刊投稿能評職稱嗎

　　由流量等效原理建立的連續(xù)介質(zhì)模型10將孔隙及裂隙的滲透性均勻分配到連續(xù)的介質(zhì)中，適合模擬宏觀裂隙不顯著的滲流現(xiàn)象。裂隙網(wǎng)絡(luò)模型在考慮裂隙滲流時具有一定優(yōu)勢，該模型將裂隙中的流動視為管道流。綜合兩種理論的雙重介質(zhì)模型采用雙重網(wǎng)格分別考慮裂隙及基質(zhì)的滲透性，較符合非均勻多孔介質(zhì)滲流的真實情況，其中孔隙介質(zhì)儲存水，裂隙介質(zhì)運輸水。近場動力學方法通過構(gòu)建積分方程來來描述物理場特性，可以避免傳統(tǒng)空間中微分方程尖端奇異性等問題，在模擬結(jié)構(gòu)破壞方面已取得廣泛的應(yīng)用。

　　熱傳導(dǎo)、滲流等擴散問題具有非局部特性，同樣可利用近場動力學來進行描述。Gerstle等首次建立了熱、電傳導(dǎo)等擴散問題的近場動力學模型，提出了其表達及數(shù)值解法。Bobaru等20基于Fourier定律和能量守恒建立了非穩(wěn)態(tài)熱擴散方程，并研究含裂隙結(jié)構(gòu)熱傳導(dǎo)問題。黃丹等采用拉格朗日坐標系下的歐拉方程，建立了常規(guī)態(tài)近場動力學熱擴散方程。程站起等利用近場動力學方法研究在溫度荷載作用下功能梯度材料的溫度場的變化表現(xiàn)。對于多孔介質(zhì)滲流問題，Jabakhanji等提出適用于均質(zhì)非飽和土水分擴散的近場動力學模型，并給出了水力傳導(dǎo)系數(shù)與傳統(tǒng)達西方程滲透率之間的關(guān)系及濕度通量表達。

　　周小平等利用近場動力學方法基于經(jīng)典Biot孔隙理論建立了裂隙巖體滲流模型及水壓致裂模型。近場動力學擴散理論中離散點所賦存的物理量在一定距離內(nèi)發(fā)生相互擴散，若超出臨界值則不再相互作用。在非局部理論中作用程度隨物質(zhì)點間距離增加而變?nèi)酰�若沒有考慮作用程度與距離的關(guān)系會導(dǎo)致定量計算精度受到影響。此外無法定量的微裂隙會導(dǎo)致土壤、巖石等多孔介質(zhì)非均質(zhì)性，若采用均質(zhì)的方法進行離散計算則無法體現(xiàn)介質(zhì)的非均勻性擴散。

　　本文在前人基礎(chǔ)上，建立多孔介質(zhì)滲流的近場動力學擴散方程，引入反映非局部作用特性的核函數(shù)提高計算精度，推導(dǎo)不同核函數(shù)下對應(yīng)的滲透系數(shù)，并建立基于威布爾分布理論及裂隙網(wǎng)絡(luò)的多孔介質(zhì)非均勻滲流模型。利用所建立的模型對多孔介質(zhì)滲流過程進行模擬將結(jié)果與理論解對比驗證方法的有效性。同時對非均勻基質(zhì)及裂隙中的滲流過程進行了模擬，結(jié)果表明所建立的方法可以較好的模擬多孔介質(zhì)滲流過程。

　　結(jié)論本文建立了多孔介質(zhì)滲流的非局部近場動力學模型，并利用提出的方法對不同情況的滲流過程進行模擬得到了以下結(jié)論：(1)近場動力學方法在模擬多孔介質(zhì)滲流過程具有一定的效果，由積分形式所建立的控制方程在滲流模擬問題上有較好的應(yīng)用前景。(2)多項式型核函數(shù)與其他核函數(shù)相比可以較大的提高計算精度，在基準問題模擬中結(jié)果收斂于理論解驗證了方法的可行性。(3)基于威布爾分布理論的介質(zhì)滲透系數(shù)可較好的反映異質(zhì)基質(zhì)滲透過程，并且裂隙網(wǎng)格滲透模型能夠模擬水分在裂隙中的流動。

　　參考文獻：

　　[1]CHERTKOVVY.Modellingcrackingstagesofsaturatedsoilsastheydryandshrink[J].EuropeanJournalofSoilScience200253(1):10518.

　　[2]GRECOR.Preferentialflowinmacroporousswellingsoilwithinternalcatchment:modeldevelopmentandapplication[J].JournalofHydrology,2002269(34):150168.

　　[3]JARVISN.Areviewofnonequilibriumwaterflowandsolutetransportinsoilmacropores:principles,controllingfactorsandconsequencesforwaterquality[J].EuropeanJournalofSoilScience,2010,58(3):523546.

　　[4]佘詩剛,林鵬.中國巖石工程若干進展與挑戰(zhàn)[J].巖石力學與工程學報,2014,33(3):433457SHEhigang,LINeng.SomedevelopmentsandchallengingissuesinrockengineeringfieldinChina[J].ChineseJournalofRockMechanicsandEngineering,2014,33(3):433457.

　　作者：馬鵬飛，李樹忱，王修偉，周慧穎，王曼靈，趙一民