水位驟降條件下粘土邊坡滲流特征分析

范會芳

（河南省商丘市虞城縣水利局，河南 虞城 476300）

1 邊坡概況

臨水庫動水環境條件下的粘土邊坡坡度為75°，邊坡高13 m，坡底長15 m，水庫蓄水水位在該邊坡位置處的常見水位深度為19 m，水庫排水后水位在該邊坡位置處的常見水位深度為13 m，屬典型的臨水粘土邊坡工程。構成邊坡坡體的土體主要物理力學參數和水力學參數如表1所示。

表1 粘土的物理力學和水力學參數表

2 數值模型建立

針對岸坡在水位驟降條件下的滲流特征分析，由于其主要考慮邊坡形態的不同產生的影響，在岸坡延伸方向上的岸坡形態近似一致，基底高度相同，水位驟降幅度也基本一致，因此該類岸坡工程問題可簡化為平面應變問題，針對二維條件下位驟降條件下的滲流特征分析，Geostudio軟件中SEEP模塊無疑可作為首選分析軟件，其可在建立分析模型時采用設置水位線的方式指定蓄水水位，用于對靜態水位條件下的邊坡穩定性求解和分析，文中的水位驟降條件下滲流特征分析，通過設定水壓線的方式設置蓄水水位及排水后的低水位，劃定分析模型區域，經過材料參數設置等步驟，按照表1 中給出的土體水力和物理力學參數對邊坡土體材料進行賦值，最后創建的用于岸坡滲流特征研究的分析模型如圖1所示。

圖1 分析模型圖

邊界條件設置：首先設定蓄水水位為19 m，然后通過構建函數的方式使水位經過10 d時間驟降為13 m，除上部邊界外，其它外部邊界均設置為固定約束，設置模型網格邊長尺寸為1 m，兼顧計算效率和求解效果。

3 計算結果分析

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，獲得了水位驟降前后的邊坡水壓圖以及壓力水頭與含水量與土體水壓關系曲線圖、水力傳導度與土體水壓關系曲線圖。

3.1 初始水位水壓特征

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中獲得的初始水位的邊坡水壓圖，如圖2所示。

圖2 初始水位水壓云圖

由圖2 所示初始水位條件下粘土臨水邊坡的土體內的水壓圖可知，水壓自上而下逐漸增大，與土體距離X 軸的距離成反比，最大水壓位于邊坡體的最底部，約為190 kPa，在與初始水位齊平位置處的邊坡土體中的水壓為0 kPa，邊坡頂部土體的水壓出現負值，表現為負孔隙水壓力，最大的負孔隙水壓力約為24 kPa，在坡體底部的水壓約為120 kPa這與靜水壓力大小一致，說明土體內的水壓力大小與理論結果一致，該邊坡在初始水位時能夠保持穩定。

3.2 初始水位壓力水頭特征

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中獲得的初始水位的壓力水頭圖。結果顯示初始水位條件下壓力水頭自上而下逐漸增大，與土體距離X軸的長度成反比，最大壓力水頭位于邊坡體的最底部，約為19 m，在與初始水位齊平位置處的邊坡土體中的壓力水頭為0 m，邊坡頂部土體的壓力水頭出現負值，表現為負壓力水頭，最大的負壓力水頭約為1.50 m，在坡體底部的壓力水頭約為12 m 這與蓄水深度一致，說明土體內的壓力水頭大小與理論結果一致，同樣說明該邊坡在初始水位時能夠保持穩定。

3.3 水位驟降6 m的水壓特征

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中獲得的水位驟降6 m 后的邊坡水壓圖，結果顯示水位驟降6 m 后的粘土臨水邊坡的土體內的水壓自上而下逐漸增大，與土體距離X 軸的距離成反比，最大水壓位于邊坡體的最底部，約為127 kPa，在與水位齊平位置處的邊坡土體中的水壓為0 kPa，邊坡頂部土體的水壓出現負值，表現為負孔隙水壓力，最大的負孔隙水壓力約為48 kPa，在坡體底部的水壓約為57 kPa這與靜水壓力大小一致，說明土體內的水壓力大小與理論結果一致，該邊坡在水位驟降6 m后能夠保持穩定。

3.4 水位驟降6 m的壓力水頭特征

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中獲得的水位驟降6 m 后的壓力水頭圖，結果顯示水位驟降6 m 后的粘土臨水邊坡土體內壓力水頭自上而下逐漸增大，與土體距離X 軸的長度成反比，最大壓力水頭位于邊坡體的最底部，約為13 m，在與水位齊平位置處的邊坡土體中的壓力水頭為0 m，邊坡頂部土體的壓力水頭出現負值，表現為負壓力水頭，最大的負壓力水頭約為1.50 m，在坡體底部的壓力水頭約為6 m這與蓄水深度一致，說明土體內的壓力水頭大小與理論結果一致，同樣說明該邊坡在水位驟降6 m后能夠保持穩定。

3.5 體積含水量與水壓關系

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中獲得的含水量與土體水壓關系曲線圖，如圖3所示。

圖3 體積含水量與水壓關系曲線圖

由圖3 所示的粘土臨水邊坡的土體含水量與土體水壓關系曲線圖可知，隨土體水壓的不斷增大，土體的體積含水量也隨之增大，初始時增大速率小，在越過-300 kPa時，增大速率迅速提高，整體呈凹形曲線上升的增大變化趨勢，最大體積含水量為0.45，此時的土體水壓為0 kPa，也即土體飽和，最小體積含水量為0.13，此時的土體水壓為-1e3 kPa，也即達到最大負水壓值，土體中的自由水已完全流失說明土體在對水環境具有良好的滲流作用，土體未發生破壞。

3.6 水力傳導度與水壓關系

為了研究該岸坡在水位驟降條件下的滲流特征，通過求解計算和繪制相應的初始水位和驟降水位后的水壓圖以及壓力水頭圖的方式，其中水力傳導度與土體水壓關系曲線圖，結果顯示，隨土體水壓的不斷增大，土體的水力傳導度也隨之增大，初始時增大速率小，同樣在越過-300 kPa時，增大速率迅速提高，整體呈凹形曲線上升的增大變化趨勢，最大水力傳導度為1e-6 m/s，此時的土體水壓為0 kPa，也即土體飽和，最小水力傳導度為2.6e-11 m/s，此時的土體水壓為-1e3 kPa，也即達到最大負水壓值，水力傳導度降低至最小值，土體中的自由水已完全流失，說明土體在對水環境具有良好的傳導水的作用，土體未發生破壞。

4 結論

文章依托山西省長治市長治北地區一水庫的實體岸坡工程，開展了對該實體岸坡工程在庫水位發生驟降條件下的滲流特性研究，得到結論：①水位驟降6 m后最大水壓為127 kPa，在水位線處為0 kPa，水位線以上的水壓表現為負值，初始水位時有同樣的規律，說明該邊坡在水位變化前后能夠保持穩定。②初始水位壓力水頭值最大為19 m，在水位線處為0 m，水位線以上的壓力水頭表現為負值，水位驟降后有同樣的規律，說明該邊坡在水位變化前后能夠保持穩定。③水位變化前后的水力傳導度和體積含水率隨土體水壓的增大而增大，且具有相似的變化規律，水在土體中具有良好的滲流作用，土體未發生破壞。