林中川

（重慶交通大學土木工程學院，重慶 400041）

隨著人類活動范圍的擴展及規(guī)模的擴大，為方便出行，交通行業(yè)的發(fā)展尤為重要。經(jīng)過幾十年的發(fā)展，一些簡單易攻克的問題已經(jīng)得以解決，但是一些復雜的問題還在繼續(xù)研究中，巖土工程就是其中重要的一部分。隨著巖體工程技術的快速發(fā)展，巖體工程的規(guī)模和深度都在不斷增長，需要更加重視多場耦合問題的求解研究。巖體存在于地質條件中，在地質條件中存在著初始地應力、地下水的作用，巖石中化學物質對圍巖的作用等，這些條件的相互作用，改變了原有的地層條件。在地下水豐富的地區(qū)，由于地下水的作用，在巖石中會形成溶洞、溶縫、漏斗等，且存在地下暗河的可能，因此地下水的作用不可忽略。若工程設計在地下水豐富的地區(qū)，就需要采用降低圍巖等級的辦法使隧道結構達到安全可靠，但是現(xiàn)階段這一方式主要還是依賴于一些工程經(jīng)驗，缺少理論的分析。何川[1]根據(jù)重慶主城的過江排水隧道，采用流固耦合的計算方法，運用有限差分理論，研究過江隧道在應力場與滲流場的耦合作用下裂隙巖體的變化規(guī)律；黃鑫、李術才[2]等通過屬性數(shù)學理論建立了隱伏溶洞與隧道巖層間安全厚度風險評估的屬性識別模型，用于開展隱伏溶洞與隧道巖層間安全厚度危險性分級，確保隧道工程施工安全。

1 理論的選用

選用的數(shù)值模擬軟件為FALC3D，其理論基礎為有限差分，將描述滲流問題的微分方程用各階導數(shù)的差商來代替，邊界條件也轉換為差分的形式，這樣就把微分的求解就轉換為代數(shù)方程組的求解。

2 數(shù)值計算與分析

2.1 依托工程概況

本隧道位于重慶市彭水縣高古鎮(zhèn)境內，最大埋深約350m，隧道全長5408m，為時速120Km/h 電化鐵路單線軌道，隧區(qū)主要地層有土層、頁巖、灰?guī)r夾頁巖、角礫巖等。隧址區(qū)可溶巖分布廣泛，節(jié)理裂隙發(fā)育、地表溶蝕較為嚴重，全隧道地下水較為發(fā)育。巖溶形態(tài)主要有沿構造節(jié)理裂隙發(fā)育的溶洞、溶縫和巖溶漏斗等，在隧道開挖過程中，經(jīng)超前地質預報，隧道洞身在約YDK211+900 右側探出一有壓含水溶洞。

2.2 計算模型

其中一種工況的數(shù)值計算模型如圖1 所示，水平方向取120m，豎向方向取60m，軸向方向取100m，溶洞位于隧道的右側，沿軸向50m 處。由于內含水且存在水源補給，使溶洞內有較高的水壓，溶洞周圍的孔隙水壓力為3MPa，在進行滲流分析時，假定滲流分析邊界的左右邊界、前后邊界及下邊界均為不透水邊界，隧道內表面為弱透水邊界。分析采用M-C 屈服準則。

根據(jù)地勘資料報告可知，巖土體的物理力學參數(shù)如表1 所示，溶洞大小編號如表2 所示。

表1 計算參數(shù)

表2 溶洞大小編號

2.3 隧道拱圈特征點圍巖應力隨溶洞半徑的變化特征

為研究不同半徑溶洞對隧道拱圈特征點圍巖壓力的影響，繪制拱圈的拱頂、左拱腰、右拱腰、左拱腳、右拱腳處的最大、中間、最小主應力及剪應力的影響變化曲線。隧道軸向50m 處拱圈特征點圍巖壓力變化圖如圖2 所示。

圖2 隧道軸向50m 處主應力變化曲線

由圖2 可知，溶洞大小變化對隧道左側圍巖主應力影響較小。右拱腳處產生應力集中，隨著溶洞大小的增大，右拱腳的應力也增大。拱頂應力及仰拱底應力在溶洞半徑較小時，其應力值都小于右拱腳的應力集中處，隨著溶洞的增大，拱頂應力及仰拱底主應力逐漸向右拱腳應力值逼近，拱頂及仰拱底最大剪應力在溶洞半徑為2.5m 時，三處的最大剪應力值相近；溶洞半徑在大于2.5m 后，拱頂及仰拱底最大剪應力超過右拱腳的應力值。左拱腳及右拱腰處的主應力值隨半徑的增大而增大，其增大的趨勢相同。

2.4 隧道拱圈特征點應力沿軸向變化曲線

為研究隧道拱圈不同特征點的應力變化情況，選擇離溶洞最近的右拱腰及沿軸向同一水平線的點，另還選擇出現(xiàn)應力集中的右拱腳作為分析對象，繪制出右拱腰及右拱腳的應力分布曲線，如圖3 所示。

圖3 隧道右拱腰處軸向主應力分布曲線

由圖3 可知，最大主應力沿軸向變化分三個階段：（1）不受影響階段，溶洞大小的變化對隧道拱圈特征點的最大主應力變化沒有影響或者影響很??；（2）應力減小階段，隧道沿軸向出現(xiàn)最大主應力減?。唬?）應力增大階段，隧道沿軸向出現(xiàn)應力變大的階段，拱圈特征點并在離溶洞最近的位置出現(xiàn)最大主應力的最大值。

2.5 不同半徑充水溶洞作用下對圍巖位移分析

由于隧道右側有一含水溶洞的存在，在隧道右側的特征點的位移明顯大于隧道左側特征點位移。為分析隧道拱圈特征點圍巖位移隨溶洞大小變化特征，將隧道沿軸向50m 處的隧道橫斷面上特征位置的位移隨溶洞大小的變化曲線繪制如圖4 所示。

由圖4 可以看出：隧道右側充水溶洞對隧道的右側拱腰的影響最大，在溶洞大小為1.5m 時，右側拱腰水平位移達到-3.2mm，且隨著溶洞的增大位移發(fā)生明顯的增大，在溶洞大小為4m 時，右側拱腰水平位移達到-29.7mm；右側充水溶洞的增大對隧道拱圈其他特征點的影響較小，對隧道的左拱腰水平位移幾乎沒有影響，左拱腰水平位移隨著溶洞的增大保持在同一水平。

2.6 隧道拱圈特征點水平位移沿軸向變化曲線

由圖4 可知，右側溶洞隧道的豎向位移影響較小，水平位移影響較大，為分析溶洞大小的變化對隧道軸向水平位移的影響及水平位移的影響范圍，選擇了右側拱腰及左側拱腰處的水平位移作為分析對象，對比分析右側溶洞對隧道左右側水平位移的影響趨勢并繪制出左右側拱腰隧道軸向水平位移曲線，如圖5 所示。

圖4 隧道軸向50m 處位移變化曲線

圖5 隧道軸向拱腰處水平位移

從圖5 可以看出：右側拱腰處的水平位移有向隧道內部移動的趨勢，溶洞隧道附近的拱腰有正向的影響，使溶洞附近的水平位移大小增大，隨著溶洞半徑的增大，溶洞附近的拱腰也增大；當溶洞半徑為1.5m 時，左側拱腰0 處的位移為-0.2mm，50m 處的位移-3mm，位移大小增大約15 倍；當溶洞半徑為4m 時，左側拱腰0 處的位移為-0.8mm，50m 處的位移-45mm，位移大小增大約56 倍。隨著溶洞的增大，溶洞對拱腰處的水平位移影響越大，影響范圍也隨著溶洞半徑的增大而增大。

3 結束語

（1）隧道一側有溶洞時，溶洞對隧道的應力影響沿隧道軸向大致成對稱分布，遠離溶洞位置影響較小，靠近溶洞處影響較大；（2）隨著溶洞的增大，溶洞對隧道的主應力影響大致成線形增長，溶洞越大，主應力大小越大，且隧道拱圈各特征點的主應力變化趨勢相近；（3）隧道一側有溶洞時，影響最大的是隧道離溶洞最近的位置且為水平位移的影響遠大于豎向位移，溶洞對拱圈其他位置的水平位移影響較小。