Influences of height of leakage point on uncoordinated deformation of heightening canals in cold areas
-
摘要: 寒区渠道由于水位和温度的周期性变化,渠基土经历“湿干冻融”循环作用,这将加剧加高渠道的不协调变形。因此,以北疆输水渠道加高改造工程为例,结合有限元计算结果分析了不同渗漏点高度对渠坡不协调变形的影响,主要结论如下：①渠坡不协调变形最显著的时间点出现在“冻”的阶段,但是并非出现在温度最低点,而是出现在温度由最低温度开始向上回升时；②当渗漏点在渠坡中点附近时,加高层引起的渠坡不协调变形最显著,渠坡与衬砌之间的漏空变形最大,且随着“湿干冻融”循环次数的增大而增大,并有逐渐趋于稳定的趋势；③渗漏点高于渠坡中点时,虽然渠坡冻胀变形增大,但是不协调变形引起的渠坡与衬砌之间的漏空变形反倒降低。渗漏点低于渠坡中点时,渠坡冻胀变形和不协调变形都较低。可见,渠坡冻胀量大并不一定意味着加高渠坡不协调变形大,渗漏点位置较低时,冻胀量越小且加高渠道的不协调变形也最小,是一种最理想的工况。Abstract: Due to the periodic changes of water supply and temperature in canals in high cold regions, the canal soils will undergo the "wetting-drying and freeze-thaw" cycles, which will aggravate the uncoordinated deformation of the heightening canals. The influences of different heights of leakage point on the uncoordinated deformation of canal slopes is analyzed by taking the heightening canal in northern Xinjiang as an example. The main conclusions are as follows: (1) The most significant uncoordinated deformation of the canal slopes occurs at the "freezing" stage, but it does not occur at the lowest temperature point, it occures when the temperature starts to rise from the lowest temperature. (2) When the leakage point is near the midpoint of the canal slopes, the uncoordinated deformation of the canal slopes caused by the addition of high rise is the most significant. The distance between the canal slope and the linings is the maximum, which increases with the increasing number of "wetting-drying and freeze-thaw" cycles to a constant value. (3) When the seepage point is higher than the midpoint of the canal slopes, the frost-heave deformation of the canal slopes increases, but the distance between the canal slopes and the linings caused by the uncoordinated deformation decreases. When the seepage point is lower than the midpoint of the canal slopes, the frost-heave deformation and the uncoordinated deformation of the canal slopes are both lower. It can be seen that a larger frost heave in the canal slopes does not necessarily mean a larger uncoordinated deformation in the heightening canals. When the leakage point is lower, the frost heave is the minimum and the uncoordinated deformation caused by heightening layer is also the minimum, which is an ideal working condition.
-
Keywords:
- heightening canal /
- uncoordinated deformation /
- wetting-drying and freeze-thaw /
- leakage /
- cold area
-
0. 引言
长距离调水工程是实现国家水资源优化配置的重大战略举措,对西北高寒缺水地区而言更是名副其实的生命线工程。经济的发展导致用水量增加,寒区一些渠道的供水能力已渐显不足,在原有渠道上进行加高扩建成为提高渠道输水能力的重要手段之一。以北疆供水工程为例,采用贴坡加高的方式将渠道由原来的5.6 m加高到7.5 m,由于渠道是梯形断面,加高之后的渠道供水能力得到了显著提高。
但是,渠道加高后面临的一大工程问题是渠坡的不协调变形。寒区冬季气温低,可达-30℃；渠道输水主要集中在春夏季,因此,供水和气温都是周期性变化的。若渠道防渗体系完好,无渠水渗漏,则对渠基土而言只会经历“冻融”循环。事实上,渠道在实际运行过程中难免出现渗漏,根据渠道管理部门的监测数据,不同的断面渗漏点高度不同：最高渗漏点与渠水位齐平,最低渗漏点发生在渠底。可见,渗漏断面的渠基土同时也经历了“干湿”循环。以1 a为周期,渠基土经历了“干湿”循环和“冻融”循环的耦合作用,朱洵等[1]将其概括为“湿干冻融”耦合循环[1]。
渠道渗漏点高度决定了渠道浸润线的位置,一方面,渠水入渗本身会造成渠基土性质的劣化,使得渗漏点以上和以下的渠基土性质产生差异,加剧渠坡的不协调变形。另一方面,渠水入渗之后,浸润线以下土体逐渐饱和,使得冻胀作用更加显著,进一步加剧渠坡的不协调变形。渗漏点下部的土体经历“湿干冻融”耦合循环作用,劣化程度将显著高于上部土体[1],渠坡的不协调变形可能会持续增大,对坡面衬砌产生危害,甚至导致衬砌的架空、断裂及剥落等,严重影响渠道的正常运行[2-5]。
由此可见,渗漏点高度是影响加高渠道不协调变形的主要因素之一。但是,该问题在此前并未被重点关注,主要原因在于普通地区的输水渠道无需考虑“冻融”循环或者“冻融”循环的影响较小；且基本能够保证全年供水,“干湿”循环的影响也较小,渠基土的性质常年能保持较为稳定的状态,渠坡不协调变形不明显。
本文将以北疆寒区输水渠道加高改造工程为背景,结合水-热-力三场耦合计算理论,计算分析不同渗漏点高度下渠坡的变形特性,并重点分析不同渗漏点高度对渠坡不协调变形的影响,为工程设计、施工、运行维护以及相关研究提供科学依据。
1. 模型建立
1.1 边界条件
北疆供水渠道改扩建工程中,老渠道为挖方工程,主要土料为膨胀土,渠顶高度为5.6 m,边坡坡度1∶2,如图1所示；渠道加高改造是填方工程,采用就近取材的原则,主要填筑料为戈壁料,采用贴坡加高的方法,加高之后的渠顶高度为7.5 m,如图1所示。
由加高层引起的渠坡不协调变形主要集中在老渠坡中部以上,对于北疆输水渠道,该渠坡段是图1中的DE段。其中,B为坡面加高分界点,E为渠坡顶点,D点为加高后渠坡的中部。
将DE坡段均匀选取21个点为研究点,通过分析各点的法向变形来分析DE段渠坡的变形特性。其中,A为DB中点,C为BE中点。
北疆供水渠道并非全年供水,而是只集中在春夏季供水,冬季气温很低,因此,以一年为一个周期,渠道的水分和温度都是周期性变化的,呈现出“湿干冻融”循环的特点。以一个周期为例,有限元计算中的主要边界条件如图2所示。
根据对渠基土在“湿干冻融”循环下的室内单元试验所总结的规律,渠基土的劣化在经历5个循环周期后会逐渐趋于稳定。因此,本文计算了渠道在5个“湿干冻融”循环周期内的变形特性,每个周期内“湿干冻融”四个状态对应的时间节点如表1所示。
表 1 有限元计算边界条件变化的时间节点Table 1. Time points of change of boundary conditions in finite element analysis状态 第1年 第N年 边界条件 湿 6 360(N-1)+ 6 水位从0 m上升到最高水位 150 360(N-1)+ 150 持续供水到降水前 160 360(N-1)+ 160 水位下降到0 m 干 180 360(N-1)+ 180 温度下降到0 冻 240 360(N-1)+ 240 温度下降到-20℃ 330 360(N-1)+ 330 温度上升到0℃ 融 360 360(N-1)+ 360 融化阶段 有限元计算所采用的“水-热-力”三场耦合计算方法在笔者此前的研究中已有详细的介绍和验证,渠道正常运行时的水位为5.5 m,此处不再赘述。
1.2 计算方案
渠道计算有限元网格如图3所示。理论上讲,若渠道防渗体系完好,渠道不会发生渗漏,渠道土体将没有水分的持续补充。然而,根据管理部门提供的监测资料,较多断面出现了渗漏情况。最严重的渗流断面,渗漏点高度与渠水正常运行高度持平,都为5.5 m；有的渗漏则发生在渠底,即渗漏点高度为0 m。因此,本文设计了4种计算方案,渗漏点高度分别为5.5,4.0,2.0,0 m,如图3所示。
![]()
图 3 不同计算方案对应的渗漏点高度Figure 3. Heights of leakage point corresponding to different calculation schemes图3中,方案Ⅰ是典型工况,对应的是渗漏点与渠水位齐平,都为5.5 m。方案Ⅳ也是典型工况,对应的是渗漏点在渠底,高度为0 m。方案Ⅱ对应的渗漏点高度为4.0 m,高于D点（3.7 m）,方案Ⅲ对应的渗漏点高度为2.0 m,低于D点（3.7 m）。
综上可见,4种计算方案中,5.5,4.0 m这两种方案渗漏点高于D点,2.0 m和0 m这两种方案是渗漏点低于D点,计算方案设计较为合理。
2. 渗漏点高度影响分析
2.1 不协调变形指标
如何描述不协调变形的程度,笔者[6]提出用衬砌与渠基土之间产生的漏空变形来表征。由渠道加高层引起的渠坡不协调变形主要集中在渠道的上部,即图1中的DE段,DE段衬砌与土体之间在“湿干冻融”的作用下会产生漏空,其中最大漏空值即为所定义的最大漏空变形Δvmax。
以渗漏点高度为5.5 m计算结果为例,第280天时DE段坡面法向位移如图4（a）所示,由图2可见,第280天时,温度约为-15℃,渠道正处于“冻”的阶段。因此,渠道土体冻胀变形,产生垂直于坡面向上的法向位移,DE段衬砌被顶起。
![]()
图 4 渗漏点高度5.5 m时渠坡最大漏空变形Figure 4. Maximum values of Δvmax at height of leakage point of 5.5 m由于渗漏高度为5.5 m,老渠坡中点附近的D点处于渠道渗漏点以下,为饱和区域,如图4（c）所示；位于渠顶的E点是双向冻结,易产生冻胀变形的点。因此,DE两个端点在该冻胀时间段内分别成为老渠道和加高层冻胀位移最大的点,如图4（a）所示,将DE段衬砌顶起,使得渠道土体与衬砌之间产生了漏空区域,如图4（b）所示。最大漏空变形位于B点附近,即加高层与老渠道的分界面附近,为Δvmax=3.1 cm。
根据图4中定义的最大漏空变形Δvmax,以渗漏点高度为5.5 m工况的第一个周期为例,在360 d内DE段最大漏空变形Δvmax如图5所示。从0～180 d,最大漏空变形Δvmax较小,基本可以忽略,漏空主要来源于通水后老渠坡在水荷载作用下产生了微小变形,而加高层在水位以上,无水压力作用,因此DE段产生了微小漏空。180 d以后,温度变为负温,渠基土开始产生冻胀,随着冻胀时间的延续,冻胀量持续增大,最大漏空变形Δvmax也逐渐增大,在280 d左右达到最大,为3.1 cm。随后,随着温度的回升,冻胀量变小,最大漏空变形Δvmax逐渐小。
值得注意的是,280 d时温度为-15℃,并不是全年温度最低的时间点,温度最低点出现在第240天,为-20℃,如图2所示。但是,由图5可见,第240天时,最大漏空变形Δvmax并不是最大,略小于第280天。其原因在于,第280天介于温度最低点（240天,-20℃）至温度回升至0℃（330天）之间,与第330天相比,第280天的温度-15℃远低于第330天时的0℃；与第240天相比,第280天冻结时间更长,且温度较为接近。换言之,第280天可视为每个周期内冻结时间与冻结温度综合作用下冻胀效果最显著的时间节点,此时新老渠坡漏空变形Δvmax达到最大,新老渠坡不协调变形最为显著。因此,在以下研究中,将各个周期内的第280天作为考察最大不协调变形的时间节点。
2.2 渗漏点高度的影响
图6整理了4种计算方案在每个周期280 d时的最大漏空变形Δvmax。从图6中可以看出关于Δvmax的两个重要特征：
![]()
图 6 最大漏空变形Δvmax与渗漏点高度的关系Figure 6. Relationship between Δvmax and height of leakage point（1）各周期内的最大漏空变形Δvmax在渗漏点高度为0～2 m范围内时基本都相同,处在最低的水平；然后,Δvmax随着渗漏点高度的增大而增大,在渗漏点高度为4 m时,Δvmax达到最大值；渗漏点高度继续增大到5.5 m时,最大漏空变形Δvmax反倒降低。
（2）渗漏点高度为4 m时,最大漏空变形Δvmax最大,且随着“湿干冻融”循环次数的增大而增大,并有逐渐趋于稳定的趋势。
第（2）个特征比较容易理解,以渗漏点高度4.0 m为例,如图7所示,随着“湿干冻融”循环次数的增加,渠坡冻胀量会逐年积累,最大漏空变形Δvmax随着“湿干冻融”循环次数的增加而增大,由第一年的4.4 cm逐步升高到第5年的6.8 cm。但是,Δvmax不会一直增大,而是逐渐趋向于稳定,这是由于经过多次“湿干冻融”循环作用以后,渠基土的性质会逐渐趋于稳定[1],因此渠坡与衬砌之间的漏空也会逐渐趋于稳定。
图6中第（1）个规律则比较反常：直观印象中,渗漏点高度越高,渠坡饱和区域面积越大,则渠坡上的冻胀量越大,不协调变形应越显著,即最大漏空变形Δvmax应随着渗漏点高度的增大而增大。但是,图6中的计算结果则显示渗漏点高度为4 m时Δvmax最大,渗漏点高度为5.5 m时Δvmax比4 m时要低。
渗漏点高度越高,渠坡上的冻胀变形越大,这一规律是毋庸置疑的。图8绘制了DE渠坡段土体的平均冻胀变形,用vDE来表示。由图8可见,随着渗漏点高度增大,vDE是单调增大的,在渗漏点高度为5.5 m时,vDE达到最大值；且DE段的平均冻胀变形vDE是逐年增大的,增加幅度逐年放缓,有趋于稳定的趋势。
![]()
图 8 DE段渠坡平均冻胀法向位移与渗漏点高度的关系Figure 8. Relationship between normal displacement of canal slope of DE section and height of leakage points图8显示,渗漏点高度为5.5 m时,DE段渠坡的平均冻胀变形显著高于渗漏点4.4 m。但是,渠坡上的冻胀变形越大,并不一定造成DE段不协调变形越大,不协调变形主要根据渠坡上各点的相对法向位移来决定。以第1个周期280 d这一时间节点时渗漏点高度为5.5,4.0,2 m的计算结果为例进行具体分析,分别如图4,9,10所示。
![]()
图 9 渗漏点高度4.0 m时渠坡最大漏空变形Figure 9. Maximum values of Δvmax at height of leakage point of 4.0 m![]()
图 10 渗漏点高度2.0 m时渠坡最大漏空变形Figure 10. Maximum values of Δvmax at height of leakage point of 2.0 m由图4（a）,9（a）,10（a）可见,渗漏点高度为5.5,4.0,2.0 m时,位于渠顶处的E点法向冻胀变形基本相同,都约为14 cm左右,这是由于E点处于渠顶,高度为7.5 m,远在渗漏点以上,不受渗漏点高度的影响。此前已分析,280 d时的漏空都是由于衬砌被E、D两点顶起,因此,在E点冻胀变形都一样的情况下,最大漏空变形Δvmax主要取决于D点的冻胀变形和DE段渠基土的冻胀,3个方案对比分析如下：
（1）方案Ⅰ中,D点高度为3.7 m,渗漏点高度为5.5 m,D点处于渗漏点以下的饱和区域,冻胀变形也较大,约为10 cm,如图4所示；DB段也基本都处于浸润线以下,冻胀变形同样可观,DB渠坡与衬砌之间的漏空并不大。造成漏空区主要集中在加高层的BE段,最大漏空变形为3.1 cm。
（2）方案Ⅱ中,渗漏高度为4.0 m时,D点位于渗漏点以下,冻胀量较大；但是D点以上的渠坡AB段都处于渗漏点以上,如图9（b）所示,冻胀显著小于渠顶的E点和位于浸润线以下的D点,因此,AB段漏空区域较大、且最大漏空变形也较大,如图9（a）所示,最大漏空变形为4.4 cm。
（3）方案Ⅲ中,膜后水位为0 m,整个DE段都位于渗漏点以上,没有水分补给,在初始饱和度相同的情况下,除了靠近渠顶的加高层冻胀量稍微较大,渠坡上的冻胀量都较小,反而使得整个渠坡的变形没有出现明显的冻胀变形不均匀现象,仅结合面B点附近出现漏空区域,最大漏空变形为2.2 cm,如图10所示。
可见,渗漏点接近渠道中点（D点）时,如高度为4.0 m,加高分界面B点在渗漏点以上,B点附近的渠坡冻胀变形较小,而D点和E点两个端点的冻胀变形较大,反倒使得渠坡与衬砌之间的漏空最大。渗漏点高于渠道中点（D点）时,如高度为5.5 m,虽然DE段渠坡的冻胀变形最大,但是DE段渠坡上各点的冻胀变形都很大,使得渠坡冻胀较为均匀,渠坡与衬砌之间的漏空反倒降低。渗漏点低于渠道中点（D点）时,接近如渗漏点在0～2 m时,DE段都在渗漏点以上,没有水分补充,冻胀变形较为均匀,渠坡与衬砌之间的漏空最小。
综上所述,渠坡冻胀量大并不一定意味着加高渠坡不协调变形大,渗漏点位置较低甚至没有渗漏时,冻胀量最小且加高渠道引起的不协调变形也最小,是一种最理想的工况。
3. 结论
本文以北疆寒区加高渠道为例,结合水-热-力三场耦合计算理论,重点分析了不同渗漏点高度对渠坡不协调变形的影响,主要结论如下：
（1）每个“湿干冻融”周期内,渠坡不协调变形最显著是出现在“冻”的阶段,来自于渠坡的不均匀冻胀。但是时间点并非是在温度最低点,而是出现在温度由最低温度开始向上回升时。
（2）当渗漏点在渠坡中点附近时,加高层引起的渠坡不协调变形最显著,渠坡与衬砌之间的漏空变形最大,且随着“湿干冻融”循环次数的增大而增大,有逐渐趋于稳定的趋势。
（3）渗漏点高于渠坡中点时,虽然渠坡冻胀变形增大,但是不协调变形引起的渠坡与衬砌之间的漏空变形反倒降低。渗漏点高于渠坡中点时,渠坡冻胀变形和不协调变形都较低。
-
![]()
图 3 不同计算方案对应的渗漏点高度
Figure 3. Heights of leakage point corresponding to different calculation schemes
![]()
图 4 渗漏点高度5.5 m时渠坡最大漏空变形
Figure 4. Maximum values of Δvmax at height of leakage point of 5.5 m
![]()
图 6 最大漏空变形Δvmax与渗漏点高度的关系
Figure 6. Relationship between Δvmax and height of leakage point
![]()
图 8 DE段渠坡平均冻胀法向位移与渗漏点高度的关系
Figure 8. Relationship between normal displacement of canal slope of DE section and height of leakage points
![]()
图 9 渗漏点高度4.0 m时渠坡最大漏空变形
Figure 9. Maximum values of Δvmax at height of leakage point of 4.0 m
![]()
图 10 渗漏点高度2.0 m时渠坡最大漏空变形
Figure 10. Maximum values of Δvmax at height of leakage point of 2.0 m
表 1 有限元计算边界条件变化的时间节点
Table 1 Time points of change of boundary conditions in finite element analysis
状态 第1年 第N年 边界条件 湿 6 360(N-1)+ 6 水位从0 m上升到最高水位 150 360(N-1)+ 150 持续供水到降水前 160 360(N-1)+ 160 水位下降到0 m 干 180 360(N-1)+ 180 温度下降到0 冻 240 360(N-1)+ 240 温度下降到-20℃ 330 360(N-1)+ 330 温度上升到0℃ 融 360 360(N-1)+ 360 融化阶段 
下载: 导出CSV
-
[1] 朱洵, 蔡正银, 黄英豪, 等. 湿干冻融耦合循环作用下膨胀土力学特性及损伤演化规律研究[J]. 岩石力学与工程学报, 2019, 38(6): 1233-1241. ZHU Xun, CAI Zheng-yin, HUANG Ying-hao, et al. Research on mechanical properties and damage evolution law of expensive soils under cyclic action of coupling wetting-drying and freeze-thaw[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(6): 1233-1241. (in Chinese)
[2] 张晨, 蔡正银, 黄英豪, 等. 输水渠道冻胀离心模拟试验[J]. 岩土工程学报, 2016, 38(1): 109-117. ZHANG Chen, CAI Zheng-yin, HUANG Ying-hao, et al. Centrifuge modelling of frost-heave of canals[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(1): 109-116. (in Chinese)
[3] 王羿, 王正中, 刘铨鸿, 等. 寒区输水渠道衬砌与冻土相互作用的冻胀破坏试验研究[J]. 岩土工程学报, 2018, 40(10): 1799-1808. WANG Yi, WANG Zheng-zhong, LIU Quan-hong, et al. Experimental investigations on frost damage of canals caused by interaction between frozen soils and linings in cold regions[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(10): 1799-1808. (in Chinese)
[4] 王正中, 刘少军, 王羿, 等. 寒区弧底梯形衬砌渠道冻胀破坏的尺寸效应研究[J]. 水利学报, 2018, 49(7): 803-813. WANG Zheng-zhong, LIU Shao-jun, WANG Yi, et al. Size effect on frost heave damage for lining trapezoidal canal with arc-bottom in cold regions[J]. Journal of Hydraulic Engineering, 2018, 49(7): 803-813. (in Chinese)
[5] 李爽, 王正中, 高兰兰, 等. 考虑混凝土衬砌板与冻土接触非线性的渠道冻胀数值模拟[J]. 水利学报, 2014, 45(4): 497-503. LI Shuang, WANG Zheng-zhong, GAO Lan-lan, et al. Numerical simulation of canal frost heaving considering nonlinear contact between concrete lining board and soil[J]. Journal of Hydraulic Engineering, 2014, 45(4): 497-503. (in Chinese)
-
期刊类型引用(1)
1. 蔡正银，张晨，朱洵，黄英豪，王羿. 高寒区长距离供水工程能力提升与安全保障技术. 岩土工程学报. 2022(07): 1239-1254 . 
本站查看
其他类型引用(1)
计量
- 文章访问数: 164
- HTML全文浏览量: 23
- PDF下载量: 59
- 被引次数: 2
