牟其海，李 帥，張 爽

(1.沈陽(yáng)市市政工程設(shè)計(jì)研究院，遼寧沈陽(yáng)110819；2.東北大學(xué)資源與土木工程學(xué)院，遼寧沈陽(yáng)110819)

0 引言

公路工程建設(shè)中，一般采用疏降地下水來(lái)確保路基和區(qū)間施工順利進(jìn)行，尤其是對(duì)于廣西這樣年平均降水量為1 530.1 mm[1]的地區(qū)。大規(guī)模降水一方面增加了公路建設(shè)成本，同時(shí)降水過(guò)度引起地表不均勻沉降，影響附近地表建筑物、管網(wǎng)等設(shè)施的安全。所以，開(kāi)展含高富水土層公路路基的施工方案優(yōu)化及理論研究，對(duì)于提高公路建設(shè)的經(jīng)濟(jì)效益，保證施工安全，縮短施工周期都具有重大的實(shí)際意義。

二灰土在公路路基建設(shè)過(guò)程中有很多優(yōu)勢(shì)特性，國(guó)內(nèi)外許多學(xué)者開(kāi)展了大量的理論和實(shí)驗(yàn)研究。崔小軍[2]在“二灰改良高含水量土用作路基填料的試驗(yàn)研究”中提到，灰改良高含水量土的強(qiáng)度、壓實(shí)度、水穩(wěn)定性都能滿足規(guī)范要求，在公路工程中可以應(yīng)用；并且給出生石灰粉∶粉煤灰=6∶18時(shí)CBR值達(dá)到峰值，隨二灰劑量繼續(xù)增加，CBR值會(huì)逐漸下降。此外，路面基層材料中鹽分聚集能降低路面強(qiáng)度，從而導(dǎo)致天然路面的不規(guī)則變形。陳陽(yáng)亭[3]等在“濱海鹽漬土路基改良試驗(yàn)研究”中提到，二灰改良濱海鹽漬土中，二灰土的最佳配比為石灰∶粉煤灰∶土=6%∶18%∶76%，并且通過(guò)實(shí)驗(yàn)證明二灰改良鹽漬土均能滿足高速公路和一級(jí)公路強(qiáng)度的要求。陸明忠[4]在“二灰土施工中對(duì)高液塑限弱膨脹土的處理應(yīng)用”中提到，二灰改良高液塑限弱膨脹土?xí)r，要根據(jù)WL、IP、FS的不同，合理選用消石灰與磨細(xì)生石灰粉砂化，可有效降低WL、IP、FS，加快土顆粒粉碎與含水量的降低，增強(qiáng)二灰土的整體強(qiáng)度，并建議：一般摻3%消石灰砂化7 d最佳含水量為30%，摻2%磨細(xì)生石灰粉砂化3 d最佳砂化含水量為40%～45%。軟土的含水量大，壓縮性強(qiáng)，具有軟塑性。楊福增[5]在“二灰法軟土路基處理”指出，采用二灰法處理軟土路基，一般工后沉降不超過(guò)5 cm，或基本上不沉降；通過(guò)實(shí)驗(yàn)得出石灰用量從8%增加到14%，其強(qiáng)度不斷增加，早期強(qiáng)度尤為明顯，但到14%以后，其強(qiáng)度并不明顯；粉煤灰用量從15%～45%增加含量其強(qiáng)度也不斷增加，特別是后期強(qiáng)度（1個(gè)月、2個(gè)月）比較明顯，但45%以后，強(qiáng)度反而下降。郁萬(wàn)彬采用ABAQUS軟件，建立等級(jí)公路層間接觸三維有限元模型，觸非線性力學(xué)特性進(jìn)行了分析，研究了不良層間接觸面位置和層間接觸狀態(tài)變化對(duì)等級(jí)公路結(jié)構(gòu)體系與疲勞壽命的影響[6]。劉俊新以FLAC2D為基礎(chǔ)建立考慮水氣兩相非飽和滲流的正演有限元模型，利用解析解對(duì)坡面產(chǎn)流數(shù)學(xué)模型的數(shù)值解進(jìn)行了驗(yàn)證[7]。薛新華以軟土路基有限元簡(jiǎn)化模型為例，利用自主編制的單標(biāo)量損傷變量流一固耦合有限元程序EPDAP_S、雙標(biāo)量損傷變量有限元程序EDAP_D以及有限變形損傷程序EDAP_L，對(duì)計(jì)算模型的沉降、孔壓、應(yīng)力、應(yīng)變以及損傷和彈性能釋放率等進(jìn)行分析討論，并以典型節(jié)點(diǎn)為例，計(jì)算了模型沉降、孔壓值等的變化曲線[8]。付振東基于多尺度有限元法(MSFEM)的基本思想，構(gòu)建了飽和多孔介質(zhì)流體相多尺度基函數(shù)，考慮非均質(zhì)飽和多孔介質(zhì)各非均勻參數(shù)場(chǎng)如滲透率、彈性模量、密度、孔隙率等[9]。

鑒于以上觀點(diǎn),采用多物理場(chǎng)耦合分析軟件Comsol Multiphysics，基于Boit固結(jié)理論，建立了考慮孔隙率受土骨架變形影響，及流體壓力和路基受車(chē)輛荷載共同作用均質(zhì)飽和多孔介質(zhì)動(dòng)態(tài)變形過(guò)程的二維滲流固結(jié)耦合模型，以實(shí)際工程為背景，對(duì)經(jīng)過(guò)二灰土改良后的高富水砂層公路路基與一般沒(méi)有處理的高富水砂層公路路基進(jìn)行對(duì)比；采用沉降位移，水壓力分布等地表變形指標(biāo)來(lái)綜合評(píng)價(jià)二灰土改良公路路基的有效性。

1 理論模型

1.1 應(yīng)力-滲流耦合理論模型

1.1.1 固體力學(xué)平衡方程

假設(shè)多孔介質(zhì)為各向同性的線彈性介質(zhì)，滿足如下以總應(yīng)力σij，應(yīng)變?chǔ)舏j，孔隙流體壓力的改變P表示的本構(gòu)方程[10]：

式（1）中：G 為剪切模量；ν為泊松比；δij為 Kronecker符號(hào)；α 是 Biot’s系數(shù)。

1.1.2 滲流連續(xù)性方程

非穩(wěn)定滲流連續(xù)性方程，考慮土體骨架變形對(duì)滲流的影響[11]，即

式（2）中：k為土體的滲透系數(shù)；ρ為土體的密度，kg/m3；D 為土體厚度，m；g是重力加速度，m2/s；φ 為水勢(shì)，m；I0為匯源項(xiàng)，kg·m-3·s-1；Ss為比貯水系數(shù)，定義為：

式（3）中：a為多孔介質(zhì)的體積壓縮系數(shù)；n是土體的孔隙率；β為土體的孔隙水壓縮系數(shù)。

式（1）為固體平衡的控制方程，體現(xiàn)了流體對(duì)固體體力平衡的影響；另外，土層的孔隙率又與所處的應(yīng)力狀態(tài)有關(guān)，這種關(guān)系可以表示為[12]：

式（4）中：φ0為零應(yīng)力狀態(tài)時(shí)的孔隙率；αφ為孔隙率的應(yīng)力敏感系數(shù)，其可以取值為5.0×10-8Pa-1；φr為高壓縮應(yīng)力狀態(tài)下的孔隙率；σ1、σ2、σ3分別為 3個(gè)主應(yīng)力，Pa；α 為 Biot’s系數(shù)。

式（2）包含了固體變形對(duì)滲流的影響。二者聯(lián)立就可解決地下水滲流的流固耦合問(wèn)題。

1.1.3 土體的線性本構(gòu)模型

廣西地區(qū)以強(qiáng)透水土層為主，土體孔隙率大，滲透性強(qiáng)，所以在流固耦合計(jì)算過(guò)程中，將土體骨架變形簡(jiǎn)化為線彈性變形，因此不考慮土體變形而引起的滲透性的耦合響應(yīng)，假設(shè)土體的滲透系數(shù)為一個(gè)常量。另外，基于Boit固結(jié)理論，土體骨架變形是微小的、線彈性的，因此不考慮損傷對(duì)流固參數(shù)的影響。

2 應(yīng)用實(shí)例

2.1 工程概況

廣西新建鐵路客運(yùn)專線玉鐵線的公路道路改移中心點(diǎn)里程為DK116+100，路基寬5m，長(zhǎng)30m，中心填高3.05m，含水層為中-粗砂、礫砂層，屬孔隙潛水。地層滲透系數(shù)為27.6m/d，地下水主要由河流的側(cè)向滲流及大氣降水為補(bǔ)給來(lái)源。不良地質(zhì)水文表現(xiàn)為地下水滲透，逕向流動(dòng)作用下對(duì)附近地層產(chǎn)生擾動(dòng)，影響路基基底的穩(wěn)定和產(chǎn)生路面沉降。路堤橫斷面如圖1所示。

圖1 路堤橫斷面示意圖（單位：m）

2.2 數(shù)值模型

2.2.1 模型介紹

根據(jù)施工資料，對(duì)路堤區(qū)域的土層進(jìn)行概化離散，建立動(dòng)態(tài)沉降數(shù)值模型如圖2所示，模型尺寸(5.0+17.0)×4.05×0.5（m2）,劃分為 2160個(gè)單元。

2.2.2 邊界條件

滲流場(chǎng)：模型內(nèi)部定為初始的定水頭邊界，砂土層內(nèi)部初始水壓力0.5MPa，左、右及下表面均采用零通量；二灰土層內(nèi)部初始水壓力0.0MPa，左、右表面采用零通量；水泥穩(wěn)定碎石層內(nèi)部初始水壓力0.0MPa，左、右表面采用零通量；瀝青混凝土層的上表面采用零通量，見(jiàn)圖2。

應(yīng)力場(chǎng)：模型砂土層側(cè)面約束x方向的位移，底部約束x方向的位移，上表面自由；二灰土層和水泥穩(wěn)定碎石層的左右及上下面為自由；瀝青混凝土層上表面受軸距為1.8m的對(duì)稱車(chē)輛荷載[13]，為100kN，側(cè)面及下表面自由見(jiàn)圖2。

圖2 動(dòng)態(tài)路堤沉降模型及邊界條件

2.3 模型計(jì)算結(jié)果分析

二灰土應(yīng)用于路堤基底土最終目的是降低基底土層的滲透性，進(jìn)而減小水在土層中流動(dòng)，進(jìn)而提高路堤的強(qiáng)度和減小路堤路面沉降，從而起到預(yù)防路堤發(fā)生失穩(wěn)破壞的作用。但采用該方法進(jìn)行路堤地基土處理，尚存在需進(jìn)行探索和研究的問(wèn)題，即二灰土改良土減透效果的影響參數(shù)有哪些，非飽和二灰土改良土對(duì)減透效果影響如何?；诼返掏翆痈黜?xiàng)異性滲透的表征和定量化，來(lái)研究實(shí)施二灰土改良后各項(xiàng)異性滲透對(duì)于路堤非飽和土層中水和氣體的運(yùn)移規(guī)律的影響。表1為土體的物理力學(xué)以及水文地質(zhì)參數(shù)。

表1 土體的物理力學(xué)以及水文地質(zhì)參數(shù)

2.3.1 實(shí)施二灰土改良模擬結(jié)果

圖3明顯看出在車(chē)輛車(chē)輪荷載作用時(shí)首先在路面周?chē)霈F(xiàn)（藍(lán)色）拉應(yīng)力，0.7MPa；隨后在路堤自重力作用下在坡腳產(chǎn)生（紅色）壓應(yīng)力。由于路堤土層自重和地基土地應(yīng)力的共同作用，靠近路堤基底附近的砂土也呈現(xiàn)出較強(qiáng)的拉應(yīng)力，藍(lán)色部分表示土層出現(xiàn)拉應(yīng)力區(qū)域，紅色部分表示土層出現(xiàn)壓應(yīng)力區(qū)域。

圖3 二灰土改良路基正應(yīng)力σy分布

路堤土層含水壓力的高低，是判斷失穩(wěn)破壞潛能大小的主要標(biāo)志，也是選擇防失技術(shù)措施的重要依據(jù)。同時(shí)，防失技術(shù)措施執(zhí)行效果好壞，關(guān)鍵表現(xiàn)在水壓力的降低程度上。為了能夠明顯體現(xiàn)高富水土層實(shí)施二灰土改良措施后，高滲透性土層水壓力的提高減緩，本節(jié)采用二灰土改良路基底層后水壓力變化時(shí)間歷程模型。圖4為模擬得到的實(shí)施二灰土改良路基基底土層后（1y,10y,50y）水壓力變化規(guī)律。假設(shè)改良前砂土層中水自由流動(dòng)，其初始水壓力為0.50 MPa左右。結(jié)果表明，在T=1y時(shí)，由于路堤水泥穩(wěn)定碎石土層中的孔隙水壓很低，水流從砂土層中向路堤土層內(nèi)不斷滲流，經(jīng)過(guò)50y的土體滲流，由于基底二灰土層的低滲透性，導(dǎo)致路堤水泥穩(wěn)定碎石土層水壓力沒(méi)有變化，圖4中的紅色箭頭代表水流滲流方向。

圖4 二灰土改良路基動(dòng)態(tài)水壓力分布

圖5 二灰土改良路基動(dòng)態(tài)水壓力曲線

圖6 二灰土改良路基動(dòng)態(tài)土體沉降量曲線

由于二灰土墊層的存在，路堤土水壓力沒(méi)有升高很快，如圖5所示，在路堤中心填埋高度1.5～4 m左右范圍內(nèi)水壓很小。在路堤使用年限達(dá)到1 a時(shí)，其水壓在14.9 Pa；在達(dá)到10 a時(shí)，其水壓在25.0 Pa;在達(dá)到50 a時(shí)，其水壓在70.0 Pa。在二灰土填埋高度1～1.5 m處，其土層內(nèi)的水壓力變化很大，越靠近高富水砂土層其壓力值越高，但是在總體上，二灰土層的水壓力隨時(shí)間變化不大。例如，路基填埋高度在1.3 m范圍內(nèi)，路堤使用年限達(dá)到1 a時(shí)，二灰土層內(nèi)的水壓力為4.6 kPa；使用年限達(dá)到10 a時(shí)，土層內(nèi)的水壓力為6.0 kPa；使用年限達(dá)到50 a時(shí)，土層內(nèi)的水壓力為55.6 kPa?？梢?jiàn)二灰土改良高富水砂土效果明顯，提高路堤使用年限。

路堤土體的沉降量是評(píng)價(jià)路堤穩(wěn)定的一個(gè)重要指標(biāo)。根據(jù)圖6所示，在路堤中心填埋高度1.5～4 m范圍內(nèi)沉降量很小，變化幅度不大。在路堤使用年限達(dá)到1 a時(shí)，路堤水泥穩(wěn)定碎石土層的沉降值穩(wěn)定在6.0 mm左右；在路堤使用年限達(dá)到10 a時(shí)，其路堤土層的沉降值穩(wěn)定在6.2 mm左右;在路堤使用年限達(dá)到50 a時(shí)，其路堤土層的沉降值穩(wěn)定在6.6 mm左右。另外，路堤沉降的動(dòng)態(tài)變化很小，經(jīng)過(guò)50 a后沉降值增加0.4 mm。由此可見(jiàn)，由于二灰土改良路堤基底土體，減緩高富砂土層的水流的流動(dòng)，沒(méi)有弱化路堤本體強(qiáng)度。隨著時(shí)間的推移，路堤土體的沉降量增加不大，為路堤施工期間及投入使用期的安全提供依據(jù)。

圖7為實(shí)施二灰土改良路基基底土層后在T=1y,10y,50y時(shí)路基沉降的動(dòng)態(tài)變化規(guī)律。在路堤使用年限達(dá)到1 a時(shí)，路堤回填土層的最大沉降值為6.294 mm，在路堤坡腳兩側(cè)的土層及靠近路基基底土層的沉降很小;在路堤使用年限達(dá)到10 a時(shí)，路堤回填土層的最大沉降值為6.394 mm；在路堤使用年限達(dá)到50 a時(shí)，路堤回填土層的最大沉降值為6.837 mm。路堤坡腳兩側(cè)的土層及靠近路基基底土層的沉降和路堤使用年限達(dá)到1 a時(shí)的沉降相比，變化不大。

2.3.2 沒(méi)有實(shí)施二灰土改良模擬結(jié)果

在車(chē)輛荷載和土體自重共同作用下，首先在路面周?chē)霈F(xiàn)拉應(yīng)力區(qū)，而在車(chē)輪作用的地方拉應(yīng)力更加明顯，其值為0.7 MPa；隨后在路堤坡腳產(chǎn)生(紅色)壓應(yīng)力。由于路堤土層自重和地基土地應(yīng)力的共同作用，靠近路堤基底附近的砂土也呈現(xiàn)出較強(qiáng)的拉應(yīng)力。與圖3經(jīng)過(guò)二灰土改良過(guò)的路堤相比，由于水的滲透，對(duì)路基的侵蝕，使得壓應(yīng)力和應(yīng)力區(qū)域擴(kuò)大。藍(lán)色部分表示土層出現(xiàn)拉應(yīng)力區(qū)域，紅色部分表示土層出現(xiàn)壓應(yīng)力區(qū)域。圖8為一般路基正應(yīng)力σy分布。

圖7 二灰土改良路基動(dòng)態(tài)沉降值

圖8 一般路基正應(yīng)力σy分布

一般含有高富水層的軟土路基，由于水勢(shì)差的存在，使得水流從壓力相對(duì)高的砂土層向壓力相對(duì)低的路堤土層滲流。如圖9所示，在路堤中心填埋高度1.0～4 m范圍內(nèi)，使用年限達(dá)到1 a時(shí)，水壓304.6～276.5 kPa；在使用年限達(dá)到10 a時(shí)，其水壓在330.0～307.1 kPa；在使用年限達(dá)到50 a時(shí)，其水壓在442.8 kPa左右；在路堤填埋高度1～1.5 m處，高富水的砂土層內(nèi)的水壓力隨時(shí)間變化不大，基底越靠近高富水砂土層其壓力值越高。例如，路基填埋高度在0.5 m范圍內(nèi)，路堤使用年限達(dá)到1 a時(shí)，砂土層內(nèi)的水壓力為470.9 kPa；使用年限達(dá)到10 a時(shí)，砂土層內(nèi)的水壓力為467.0 kPa；使用年限達(dá)到50 a時(shí)，砂土層內(nèi)的水壓力為460.7 kPa。另外，砂土層高度從0～1.0 m處，水壓493.9～305.5 kPa。

圖10為模擬得到的一般路基（1y,10y,50y）時(shí)水壓力變化規(guī)律。假設(shè)高富水砂土層中水自由流動(dòng)，其初始水壓力為0.50 MPa左右。由于壓差的作用，使得水流向路堤本體穩(wěn)定滲流。結(jié)果表明，在T=1y時(shí)，由于路堤水泥穩(wěn)定碎石土層中的孔隙水壓很低，水流從砂土層中向路堤土層內(nèi)不斷滲流，此時(shí)路堤土層內(nèi)的水壓最小值為0.277 MPa，高富水砂土層的水壓值為0.499 MPa；經(jīng)過(guò)10y的土體滲流，路堤土層內(nèi)的水壓最小值為0.307 MPa，高富水砂土層的水壓值為0.496 MPa；經(jīng)過(guò)50y的土體滲流，路堤土層內(nèi)的水壓最小值為0.443 MPa，高富水砂土層的水壓值為0.491 MPa。可見(jiàn)，和圖5相比，沒(méi)有經(jīng)過(guò)二灰土改良的路基土層內(nèi)的水流滲流很快，在T=50y路堤土層基本達(dá)到飽和，這對(duì)于路堤的穩(wěn)定是極為不利的，圖中紅色代表水流滲流方向。

圖9 一般路基動(dòng)態(tài)水壓力曲線

圖11為含有高富水砂土層的路基在1y,10y,50y時(shí)路基沉降的動(dòng)態(tài)變化規(guī)律。在路堤使用年限達(dá)到1 a時(shí)，路堤回填土層的最大沉降值為19.5 mm，在路堤坡腳兩側(cè)的土層及靠近路基基底土層的沉降很??；在路堤使用年限達(dá)到10 a時(shí)，路堤回填土層的最大沉降值為20.6 mm；在路堤使用年限達(dá)到50 a時(shí)，路堤回填土層的最大沉降值為25.9 mm。路堤坡腳兩側(cè)的土層及靠近路基基底土層的沉降和路堤使用年限達(dá)到1 a時(shí)的沉降相比，變化不大。圖11中紅色箭頭代表路基沉降方向。

圖10 一般路基動(dòng)態(tài)水壓力分布

圖11 一般路基動(dòng)態(tài)沉降值

在路基基底含有高富水砂土層的路堤中，由于水的滲流作用，使得路堤土體的強(qiáng)度降低，引起路面沉降。根據(jù)圖12所示，在路堤中心填埋高度1.5～4 m左右范圍內(nèi)，沉降量很小,變化幅度不大。在路堤使用年限達(dá)到1 a時(shí)，路堤水泥穩(wěn)定碎石土層的沉降值穩(wěn)定在19.0 mm左右；在路堤使用年限達(dá)到10 a時(shí)，其路堤土層的沉降值穩(wěn)定在20.5 mm左右;在路堤使用年限達(dá)到50 a時(shí)，其路堤土層的沉降值穩(wěn)定在25.8 mm左右。另外，路堤沉降的動(dòng)態(tài)變化很大，經(jīng)過(guò)50 a后，沉降值增加6.8 mm。

圖12 一般路基動(dòng)態(tài)土體沉降量曲線

3 結(jié)論

基于Biot固結(jié)理論對(duì)高富水砂層公路路基和經(jīng)過(guò)二灰土改良后的高富水砂層公路路基進(jìn)行數(shù)值模擬對(duì)比，得出了以下結(jié)論。

（1）經(jīng)過(guò)二灰土改良后的高富水砂層公路路基比一般沒(méi)有處理的高富水砂層公路路基的路面沉降顯著減少。使用年限1 a時(shí)，減少13.21 mm；使用年限10 a時(shí)，減少14.21 mm;使用年限50 a時(shí)，減少19.06 mm。另外，隨著使用年限的增加，經(jīng)過(guò)二灰土路基改良后的路面沉降值變化很小，從6.294～6.837 mm，而沒(méi)有改良的路面沉降變化很大，從19.5～25.9 mm。

（2）有二灰土基底路堤，由于其滲透率低，整體性及密實(shí)度高等特點(diǎn)，減緩水流滲透作用，從而使得路堤土層內(nèi)的水壓上升緩慢，經(jīng)過(guò)10 a后，水壓上升10.1 Pa；而沒(méi)有經(jīng)過(guò)處理的路堤，經(jīng)過(guò)10 a后，水壓上升25.4 kPa；

（3）車(chē)輛荷載、土體自重及水力梯度共同作用下，對(duì)比路堤基底二灰土改良和沒(méi)有進(jìn)行改良的正應(yīng)力σy分布圖，靠近路堤基底附近的砂土層及車(chē)輛荷載作用路面處產(chǎn)生的拉應(yīng)力區(qū)域要顯著增加，另外，在路堤坡腳處的壓應(yīng)力區(qū)域也明顯擴(kuò)大。

[1]涂方旭.廣西近80年來(lái)的降水量序列及變化特征[J].廣西水利水電，1998(4):21-25.

[2]崔小軍.二灰改良高含水量土用作路基填料的試驗(yàn)研究[J].路基工程，2009(2):188-189.

[3]陳陽(yáng)亭,廖陳林,鄭慧振.濱海鹽漬土路基改良試驗(yàn)研究[J].公路與汽車(chē),2009(6):71-73.

[4]陸明忠,宋曄,柴如飛.二灰土施工中對(duì)高液塑限弱膨脹土的處理應(yīng)用[J].浙江建筑,2008,25(8):43-46.

[5]楊福增.二灰法軟土路基處理[J].石家莊鐵路職業(yè)技術(shù)學(xué)院學(xué)報(bào),2006,5(增刊):6-9.

[6]郁萬(wàn)彬.黃泛區(qū)二灰土底基層道路結(jié)構(gòu)界面設(shè)計(jì)研究[M].濟(jì)南：山東大學(xué),2010.

[7]劉俊新.非飽和滲流條件下紅層路堤穩(wěn)定性研究[D].成都：西南交通大學(xué),2007.

[8]薛新華.巖土介質(zhì)流一固耦合非線性損傷力學(xué)理論與數(shù)值分析[D].杭州：浙江大學(xué),2008.

[9]付振東.非均質(zhì)飽和多孔介質(zhì)準(zhǔn)靜態(tài)行為分析的耦合多尺度及耦合升尺度有限元法[D].大連：大連理工大學(xué),2010.

[10]章夢(mèng)濤,潘一山,梁冰,等.煤巖流體力學(xué)[D].北京:科學(xué)出版社,1995.

[11]魏晨慧.巖體開(kāi)挖損傷區(qū)的熱-流-力耦合模型研究 [M].沈陽(yáng)：東北大學(xué),2008.

[12]Rutqvist J,Tsang C-F.A study of caprock hydromechanical changes associated with CO2-injection into a brine formation[J].Enviromental Geology,2002,42:296-305.

[13]JTG B01-2003，公路工程技術(shù)標(biāo)準(zhǔn)[S].