New technology of fiber leakage monitoring in long-distance water transfer channels in alpine and cold regions
-
摘要: 目前传统的监测技术在长距离输水渠道渗流监测上存在明显的不足,基于拉曼光时域反射原理,研发了一种新型渗流监测光缆,提出了可用于高寒区长距离输水渠道的分布式光纤拉曼渗流监测技术。针对空气、不同含水率的砂土、纯水等不同现场环境进行了室内试验模拟研究,对比分析光纤温度变化,结果表明在不同的环境中该新型光缆均可反映出环境特征,证实了光纤监测的可行性。在北疆干渠开展新型光纤渗漏监测技术现场原型试验,并在光纤埋设处对应埋设了传统的孔隙水压力计,对比分析并验证分布式光纤测渗漏技术应用于渠道渗漏监测的可行性及测量精度。本次现场测量光纤采样点间距为0.5 m,光纤测量的平均值与孔隙水压力计测值差小于0.5 m,差值在理论范围内,说明光纤渗漏监测技术能够适应现场测量条件,证实了其在高寒区长距离输水渠道应用的准确性、可行性和可靠性。Abstract: In the seepage monitoring of long-distance water transfer channels, the traditional monitoring technology has obvious shortcomings. Based on the principles of Raman optical time domain reflection, a new type of seepage monitoring optical cableh is developed, and a distributed optical fiber Raman seepage monitoring technology which can be used in the long-distance transfer channels in cold regions is proposed. The laboratory tests are carried out for different field environments such as air, sand with different water content, pure water, etc. The results show that the new optical fiber cable can reflect the environmental characteristics under different environments, which confirms the feasibility of the optical fiber monitoring. By using the new fiber-optic leakage monitoring technology, field prototype tests are carried out in the main canal of North Xinjiang, and the pore pressure cells are installed in the fiber-optic buried place. The feasibility and measurement accuracy of the distributed fiber-optic leakage monitoring technology for the cannel leakage monitoring has been analyzed and verified. The distance between the optical fiber sampling points in the field measurement is 0.5 m, and the difference between the average value of optical fiber measurement and the pore water pressure gauge measurement is less than 0.5 m, which is within the theoretical range. It is shown that the optical fiber leakage monitoring technology is suitable for the field measurement conditions, and the accuracy, feasibility and reliability of its application in the long-distance water transfer channels incold regions are validated.
-
Keywords:
- distributed fiber /
- channel in cold region /
- leakage monitoring /
- prototype test
-
0. 引言
软岩堆石料填筑碾压时,岩块易发生崩解,级配细化,且易受温度和干湿循环等环境因素影响,产生较大蠕变变形[1-3],这对高土石坝长期稳定是极其不利的。为了更准确地掌握软岩堆石坝运行期的稳定性,有必要对含软岩堆石料蠕变特性进行研究。
目前,许多学者对堆石料的蠕变特性进行了研究,得到大量有益的结论。梁军[4]探究了干密度对堆石料最终蠕变应变的影响,并建立了两者的经验关系式。孙振远等[5]通过试验发现提高干密度可减小堆石料的蠕变应变。上述研究主要考虑了干密度对堆石料蠕变特性的影响,但干密度不宜作为堆石坝填筑质量的控制指标。朱晟等[6]认为除了孔隙率,还应将相对密实度作为高土石坝填筑质量的控制指标。因此有必要研究不同相对密实度的含软岩堆石料的蠕变特性。
在本文中,以相对密实度为研究因素,采用大型压缩仪进行含软岩堆石料单轴蠕变试验,探讨不同相对密实度堆石料的蠕变规律。利用幂函数蠕变模型拟合试验结果,探究相对密实度对拟合参数的影响,进而提出考虑相对密实度的单轴蠕变模型,为高土石坝的设计和稳定性分析提供依据。
1. 试验仪器和标准
1.1 试验仪器
单轴蠕变试验采用成勘院组建的大型压缩仪进行,仪器由主机架、加载系统和计算机测控与数据采集系统等组成。加载系统采用交流伺服与液压加载方式自动加压,加载压力和变形分别采用精密力传感器和高精度变形传感器测量,计算机自动采集。仪器稳压效果好,测量精度高。仪器主要技术参数:试样尺寸Φ502 mm×252 mm,最大负荷1500 kN,最大变形50 mm。
1.2 试验标准
目前,堆石料蠕变试验尚无规范和标准可循,在本文的研究中,参考相关研究成果[7]制定蠕变起始时间判别标准和稳定标准。堆石料蠕变试验在加载1 h后,变形基本能平稳发展,因此选择加载1 h作为蠕变发生起始时间。蠕变稳定采用双标准控制,即24 h内变形量不大于0.05 mm且试验时间不少于252 h。
2. 试验土料和试验方案
2.1 试验土料
以某高心墙堆石坝工程的堆石料为研究对象,堆石料以板岩为主,夹少量条带状或透镜状变质砂岩,饱和抗压强度为25.7~140.0 MPa,软化系数为0.65~0.85。根据《GB 50021—2001岩土勘察规范》[8]分类,试验土料中含有部分软岩,遇水容易软化。
工程所用堆石料的最大粒径为800 mm,受限于试验仪器尺寸,需处理堆石料中粒径超过60 mm的颗粒。对超粒径土料采用混合法[9]处理,先将原级配的粒径缩小为1/2,再对超过60 mm的颗粒进行等量替代,试验用料的原级配与试验级配见表1。
表 1 堆石料原级配与试验级配Table 1. Grain-size distribution of rockfill materials土样级配 不同粒径(mm) 堆石料百分比/% 800~400 400~200 200~60 60~40 40~20 20~10 10~5 <5 原级配 17.5 44.0 24.5 5.5 9.0 6.5 5.5 5.0 试验级配 27.0 18.4 25.6 18.5 10.5 软岩料的压实密度受含水率影响,存在最优含水率和最大干密度[10]。制备不同含水率的含软岩堆石料,利用振动台法进行最大干密度试验。测得试验用料最大干密度为2.00 g/cm3,对应的含水率为4.6%。采用固定体积法测得最小干密度为1.57 g/cm3。
2.2 试验方案
依据文献[11]中关于土石坝填筑要求的规定,堆石料碾压时宜加水。因此,试验中的试样均为非饱和状态,采用最大干密度试验确定的含水率4.6%。一般高堆石坝堆石料的相对密实度为0.9左右,采用0.85,0.90和0.95三种相对密实度分别进行试验,每个相对密实度的试样都在0.4,0.8,1.6,3.2 MPa 4级轴向压力下进行单轴蠕变试验。
3. 结果及分析
3.1 相对密实度对堆石料轴向蠕变特性的影响
不同相对密实度单轴蠕变试验的轴向蠕变应变随时间变化的曲线见图1,图中最右侧的数值为轴向压力,压力相同的一组曲线中上、中、下分别为相对密实度0.85,0.90,0.95的曲线。在蠕变初期,蠕变应变曲线斜率较大,蠕变速率较快;随着时间增长,蠕变速率逐渐减慢,最终蠕变应变趋于一个稳定值。这是由于在蠕变初期,堆石料的蠕变应变主要由颗粒破碎和颗粒重排列共同产生。随着轴向压力作用时间延长,颗粒破碎速率逐渐减慢,直至不再发生破碎。说明在试验压力范围内,含软岩堆石料的蠕变为衰减蠕变。
随着相对密实度增大,试样的轴向蠕变应变逐渐减小,说明堆石料蠕变与其自身密实程度相关。密实程度低的堆石料孔隙率大,外力作用下,粗大颗粒棱角或软岩颗粒发生破碎、细化后,更容易滑移填充到孔隙中,产生较大蠕变应变。随着轴向压力增大,相同压力、不同相对密实度试样的蠕变应变曲线的差距逐渐增大。这是由于轴向压力较小时,颗粒表面接触应力较小,破碎颗粒数量少且粒径大,破碎颗粒填充受孔隙的影响较小,故不同密实程度的堆石料蠕变应变相差较小。随着轴向压力增大,颗粒表面接触应力增大,破碎颗粒数量增加且粒径变小,破碎颗粒滑移填充受孔隙的影响程度提高,故各试样间蠕变应变曲线差距变得明显。轴向蠕变应变曲线局部存在交叉现象,主要是由于堆石料颗粒互相咬合,若局部发生破碎或错动,变形量会有一定的激增,在蠕变应变曲线中表现较明显。
3.2 考虑轴向压力的蠕变模型
长江科学院提出的幂函数蠕变模型已成功应用于水布垭、双江口[12]的堆石料蠕变研究。说明该模型用于表达堆石料的蠕变特性具有较好的适用性,故本文的试验结果也采用该模型进行拟合。该蠕变模型为
, (1) 式中,
为轴向蠕变应变, 为最终蠕变应变, 为衰减指数,t为时间,单位为h。 利用Matlab分别对3.1节所述的12组试验获得的轴向蠕变应变曲线用式(1)拟合,得到各组试验的最终蠕变应变和衰减指数。由于拟合得到的未必是与实测蠕变应变过程线变化趋势吻合最佳的模型参数,因此需优化模型参数。利用Excel表格的试算功能,以拟合得到的两参数为基准,在±5%变化范围内分别间隔0.003%和0.001选取多个值,利用不同参数值组合试算,对比试算和试验结果,确定与实测蠕变应变过程线变化趋势吻合最佳的模型参数组合,结果见表2。
表 2 幂函数蠕变模型参数Table 2. Parameters of creep model of power function相对密实度Dr 轴向压力p/MPa 最终蠕变应变 /% 衰减指数 决定系数R2 0.85 0.4 0.295 0.090 0.9712 0.8 0.574 0.107 0.9638 1.6 0.938 0.121 0.9826 3.2 1.441 0.131 0.9729 0.90 0.4 0.268 0.089 0.9866 0.8 0.535 0.102 0.9844 1.6 0.837 0.121 0.9107 3.2 1.347 0.130 0.8501 0.95 0.4 0.237 0.087 0.9621 0.8 0.488 0.099 0.9870 1.6 0.731 0.120 0.9632 3.2 1.225 0.129 0.9829 利用表2的参数计算各相对密实度和轴向压力的轴向蠕变应变曲线,其结果与试验结果的比较见图2。由图2可知,曲线与离散点偏差小,变化趋势一致,说明幂函数蠕变模型能较好表现含软岩堆石料蠕变特性。
通过绘制最终蠕变应变随轴向压力变化的曲线可知,两者关系曲线形状近似双曲线:
, (2) 式中,p为轴向压力(MPa),a,b为拟合参数。故用双曲线来拟合两者关系,见图3,图例中T代表试验点,F代表拟合曲线。衰减指数随轴向压力的变化关系见图4,由图4可知,两者关系曲线形状近似幂函数:
, (3) 式中,c和d为拟合参数,拟合曲线也绘于图4中。
将式(2),(3)代入式(1),整理后可得考虑轴向压力p的单轴蠕变模型表达式:
。 (4) 不同相对密实度的最终蠕变应变和衰减指数随轴向压力变化曲线的拟合参数及决定系数R2如表3所示。由表3可知,R2均在0.95以上,拟合效果良好。
表 3 最终蠕变应变和衰减指数拟合曲线的参数及决定系数Table 3. Parameters and determination coefficients of fitting curve of final creep strain and attenuation index相对密实度Dr 最终蠕变应变 衰减指数 a b R2 c d R2 0.85 0.330 1.167 0.9994 0.109 0.174 0.9721 0.90 0.334 1.319 0.9973 0.107 0.184 0.9730 0.95 0.349 1.517 0.9934 0.105 0.193 0.9687 平均值 0.338 1.334 - 0.107 0.184 - 标准差 0.008 0.143 - 0.002 0.008 - 公式(4)中的各参数的标准差见表3最后一行。参数a,c,d的标准差均小于0.01,远小于参数b的标准差0.143,说明相对密实度对这三个参数的影响很小。为简化模型,将参数a,c,d取平均值,分别为0.338,0.107,0.184。由表2中不同相对密实度、相同轴向压力的衰减指数
值的对比,也可看出相对密实度对衰减指数的影响很小。 3.3 考虑轴向压力和相对密实度的蠕变模型
由表3可知参数b与相对密实度呈正相关,据此可进一步分析参数b与相对密实度的函数关系,建立考虑相对密实度的蠕变模型。将不同相对密实度的参数b绘于图5中,由图可知,其与相对密实度呈良好的线性关系,拟合的决定系数大于0.95。将图5中参数b与相对密实度的线性函数关系代入式(4),可得
, (5) 式中,a,b1,b2,c,d为模型参数,各参数取值见表4。
表 4 考虑轴向压力和相对密实度的蠕变模型参数Table 4. Parameters of creep model for rockfill considering axial pressure and relative density模型参数 a b1 b2 c d 取值 0.338 3.50 -1.816 0.107 0.184 利用考虑轴向压力和相对密实度的蠕变模型计算各轴向压力和相对密实度含软岩堆石料的轴向蠕变应变曲线,计算结果和试验结果的比较见图6。由图6可知,计算结果与试验结果偏差较小,且变化趋势一致。说明前文中求取模型参数的方法可行,模型参数可靠。
对于以相对密实度为填筑指标的堆石坝,可进行单轴蠕变试验获得堆石料轴向蠕变应变曲线。利用式(1)拟合求取各相对密实度堆石料最终蠕变应变和衰减指数,然后分别用式(2),(3)拟合两者随轴向压力变化的曲线,得到参数a,b,c,d。将参数a,c,d取平均值,将参数b表示为相对密实度的线性函数,得到参数b1和b2。为符合堆石料变形特性,参数b不能为负值,因此其与相对密实度的线性关系仅适用于密实情况下的堆石料。本文建立的蠕变模型在前文讨论的压力和相对密实度范围对含软岩堆石料适用,上述范围以外模型的适用性需进一步研究确定。单轴蠕变模型能较好表现堆石料的蠕变特性,预测堆石坝的长期变形特性。如梁军[4]利用单轴蠕变模型结合广义虎克定律,计算得到的关门山坝的沉降曲线与实测沉降过程线相似程度较高。
3.4 模型参数的物理意义
本文的模型是通过将式(1)中最终蠕变应变和衰减指数表示为轴向压力和相对密实度的函数关系建立的,因此模型参数的物理意义可通过最终蠕变应变和衰减指数表达式中各参数的物理意义来说明。由式(2)可知,1/a为轴向压力趋近于无穷大时最终蠕变应变。对式(2)求导,1/b为轴向压力为0时最终蠕变应变随轴向压力增长的速率。由图5可知,b1为参数b随相对密实度的变化速率,b2为堆石料相对密实度为0时的b值。由式(3)为幂函数可知,c为轴向压力为1 MPa时的衰减指数,d是双对数坐标下衰减指数与轴向压力关系曲线的斜率,表征衰减指数随轴向压力变化的速率。
4. 结论
本文对不同相对密实度的含软岩堆石料进行了单轴蠕变试验研究,得到以下结论:
(1)在试验压力范围内,含软岩堆石料的蠕变类型为衰减蠕变。随着相对密实度的增大,堆石料的轴向蠕变应变逐渐减小。轴向压力较低时,不同相对密实度的轴向蠕变应变差距较小;随着轴向压力提高,轴向蠕变应变的差距逐渐增大。
(2)利用幂函数蠕变模型拟合轴向蠕变应变随时间变化的曲线,得到最终蠕变应变和衰减指数。其分别与轴向压力满足双曲线和幂函数关系,由此建立了考虑轴向压力的含软岩堆石料的单轴蠕变模型。
(3)分析模型参数可知相对密实度对衰减指数的影响很小。依据与最终蠕变应变相关的模型参数b和相对密实度呈线性关系,提出了同时考虑轴向压力和相对密实度的单轴蠕变模型。利用分析确定的模型参数计算各轴向压力和相对密实度情况下堆石料的轴向蠕变应变,计算结果与试验结果较吻合。说明提出的单轴蠕变模型在本文讨论的压力和相对密实度范围对含软岩堆石料是适用的。
-
表 1 光纤测渗漏现场测值分析表
Table 1 Analysis of measured values of optical fiber leakage
日期 光纤在水位以上的长度/m 水面距离孔压的高度/m 孔压相应的测值/m 2019-07-22 温差法 56.5~66 3.25 3.175 85.5~93.5 3.61 116.5~122.5 4.03 144.5~152.5 3.62 2019-07-22 K值法 56~65.5 3.25 84.5~93 3.50 116.5~122.5 4.03 144.5~152 3.73 2019-09-19 温差法 50.5~70 1.02 0.685 82~101.5 1.04 109~130.5 0.57 140~157.5 1.50 2019-09-19 K值法 51~70 1.13 82.5~101 1.26 111~130 1.13 141.5~160.5 1.16 -
[1] 何宁, 丁勇, 吴玉龙, 等. 基于分布式光纤测温技术的堤坝渗漏监测[J]. 水利水运工程学报, 2015(1): 20-27. HE Ning, DING Yong, WU Yu-long, et al. Experimental study of distributed optical fiber temperature measurement technology for measuring lea[J]. Hydro-Science and Engineering, 2015(1): 20-27. (in Chinese)
[2] 翟明杰, 杨进新, 李启升, 等. 长距离正反向输水渠道设计若干问题探讨[J]. 人民黄河, 2019, 41(11): 97-101. doi: 10.3969/j.issn.1000-1379.2019.11.020 ZHAI Ming-jie, YANG Jin-xin, LI Qi-sheng. Discussion on engineering design of long-distance and bi-directional canal[J]. Yellow River, 2019, 41(11): 97-101. (in Chinese) doi: 10.3969/j.issn.1000-1379.2019.11.020
[3] 朱萍玉, 蒋桂林, 冷元宝. 采用分布式光纤传感技术的土坝模型渗漏监测分析[J]. 中国工程科学, 2011, 13(3): 82-85, 96. doi: 10.3969/j.issn.1009-1742.2011.03.013 ZHU Ping-yu, JIANG Gui-ling, LENG Yuan-bao. Analysis of seepage monitoring for an earth dam model with distributed optical fiber sensors[J]. Engineering Sciences, 2011, 13(3): 82-85, 96. (in Chinese) doi: 10.3969/j.issn.1009-1742.2011.03.013
[4] 肖衡林, 蔡德所. 基于温度测量的分布式光纤渗漏监测技术机理探讨[J]. 岩土力学, 2008, 29(增刊1): 550-554. XIAO Heng-ling, CAI De-suo. Discussion on mechanism of distributed fiber optical leakage monitoring technology based on temperate[J]. Rock and Soil Mechanics, 2008, 29(S1): 550-554. (in Chinese)
[5] 商峰, 范哲, 朱新民. 基于分布式光纤监测的压力钢管模拟泄漏室内试验[J]. 中国水利水电科学研究院学报, 2017, 15(3): 208-212. SHANG Feng, FAN Zhe, ZHU Xin-min. Distributed optical fiber monitoring in simulated leakage test on pressure steel pipe[J]. Journal of China Institute of Water Resources and Hydropower Research, 2017, 15(3): 208-212. (in Chinese)
[6] 李科, 施斌, 唐朝生, 等. 黏性土体干缩变形分布式光纤监测试验研究[J]. 岩土力学, 2010, 31(6): 1781-1785. LI Ke, SHI Bing, TANG Cao-sheng, et al. Feasibility research on soil deformation monitoring with distributed optical fiber sensing technique[J]. Rock and Soil Mechanics, 2010, 31(6): 1781-1785. (in Chinese)
[7] 张巍, 施斌, 索文斌, 等. 冻土瞬态温度场的分布式光纤监测方法及应用[J]. 岩土工程学报, 2007(5): 723-728. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200705015.htm ZHANG Wei, SHI Bing, SUO Weng-bin, et al. Monitoring and application of distributed optical fiber sensors in transient temperature field of frozen soil[J]. Chinese Journal of Geotechnical Engineering, 2007(5): 723-728. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200705015.htm
[8] 刘康康. 基于ROTDR的渠道渗漏监测技术研究[D]. 南京: 南京理工大学, 2018. LIU Kang-kang. Research on Monitoring Technology of Canal Leakage Based on ROTDR[D]. Nanjing: Nanjing University of Science & Technology, 2018. (in Chinese)
-
期刊类型引用(3)
1. 陈志敏,易明炀,杨志强,王壹敏. 埋深和水土含量对松散岩堆抗剪强度的影响. 铁道建筑. 2024(02): 112-117 . 百度学术
2. 张赓旺,宋嘉杰. 含水率及填充土含量对松散岩堆体抗剪强度的影响规律. 高速铁路技术. 2023(05): 34-39 . 百度学术
3. 祝艳波,刘振谦,李文杰,苗帅升,李红飞,兰恒星. 黄土-三趾马红土滑坡滑带土的长期强度影响因素研究. 水文地质工程地质. 2022(02): 148-156 . 百度学术
其他类型引用(5)