• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

高寒区长距离输水渠道光纤渗漏监测新技术研究

李登华, 张桂荣, 丁勇, 孔洋, 何宁

李登华, 张桂荣, 丁勇, 孔洋, 何宁. 高寒区长距离输水渠道光纤渗漏监测新技术研究[J]. 岩土工程学报, 2020, 42(S2): 83-87. DOI: 10.11779/CJGE2020S2015
引用本文: 李登华, 张桂荣, 丁勇, 孔洋, 何宁. 高寒区长距离输水渠道光纤渗漏监测新技术研究[J]. 岩土工程学报, 2020, 42(S2): 83-87. DOI: 10.11779/CJGE2020S2015
LI Deng-hua, ZHANG Gui-rong, DING Yong, KONG Yang, HE Ning. New technology of fiber leakage monitoring in long-distance water transfer channels in alpine and cold regions[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(S2): 83-87. DOI: 10.11779/CJGE2020S2015
Citation: LI Deng-hua, ZHANG Gui-rong, DING Yong, KONG Yang, HE Ning. New technology of fiber leakage monitoring in long-distance water transfer channels in alpine and cold regions[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(S2): 83-87. DOI: 10.11779/CJGE2020S2015

高寒区长距离输水渠道光纤渗漏监测新技术研究  English Version

基金项目: 

国家重点研发计划项目 2017YFC0405100

国家自然科学基金面上基金项目 51979174

详细信息
    作者简介:

    李登华(1980— ),男,教授级高级工程师,博士,主要从事岩土工程安全监测和特殊土改良领域的科研工作。E-mail:dhli@nhri.cn

  • 中图分类号: TU431

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.
  • 软岩堆石料填筑碾压时,岩块易发生崩解,级配细化,且易受温度和干湿循环等环境因素影响,产生较大蠕变变形[1-3],这对高土石坝长期稳定是极其不利的。为了更准确地掌握软岩堆石坝运行期的稳定性,有必要对含软岩堆石料蠕变特性进行研究。

    目前,许多学者对堆石料的蠕变特性进行了研究,得到大量有益的结论。梁军[4]探究了干密度对堆石料最终蠕变应变的影响,并建立了两者的经验关系式。孙振远等[5]通过试验发现提高干密度可减小堆石料的蠕变应变。上述研究主要考虑了干密度对堆石料蠕变特性的影响,但干密度不宜作为堆石坝填筑质量的控制指标。朱晟等[6]认为除了孔隙率,还应将相对密实度作为高土石坝填筑质量的控制指标。因此有必要研究不同相对密实度的含软岩堆石料的蠕变特性。

    在本文中,以相对密实度为研究因素,采用大型压缩仪进行含软岩堆石料单轴蠕变试验,探讨不同相对密实度堆石料的蠕变规律。利用幂函数蠕变模型拟合试验结果,探究相对密实度对拟合参数的影响,进而提出考虑相对密实度的单轴蠕变模型,为高土石坝的设计和稳定性分析提供依据。

    单轴蠕变试验采用成勘院组建的大型压缩仪进行,仪器由主机架、加载系统和计算机测控与数据采集系统等组成。加载系统采用交流伺服与液压加载方式自动加压,加载压力和变形分别采用精密力传感器和高精度变形传感器测量,计算机自动采集。仪器稳压效果好,测量精度高。仪器主要技术参数:试样尺寸Φ502 mm×252 mm,最大负荷1500 kN,最大变形50 mm。

    目前,堆石料蠕变试验尚无规范和标准可循,在本文的研究中,参考相关研究成果[7]制定蠕变起始时间判别标准和稳定标准。堆石料蠕变试验在加载1 h后,变形基本能平稳发展,因此选择加载1 h作为蠕变发生起始时间。蠕变稳定采用双标准控制,即24 h内变形量不大于0.05 mm且试验时间不少于252 h。

    以某高心墙堆石坝工程的堆石料为研究对象,堆石料以板岩为主,夹少量条带状或透镜状变质砂岩,饱和抗压强度为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~400400~200200~6060~4040~2020~1010~5<5
    原级配17.544.024.55.59.06.55.55.0
    试验级配   27.018.425.618.510.5
    下载: 导出CSV 
    | 显示表格

    软岩料的压实密度受含水率影响,存在最优含水率和最大干密度[10]。制备不同含水率的含软岩堆石料,利用振动台法进行最大干密度试验。测得试验用料最大干密度为2.00 g/cm3,对应的含水率为4.6%。采用固定体积法测得最小干密度为1.57 g/cm3

    依据文献[11]中关于土石坝填筑要求的规定,堆石料碾压时宜加水。因此,试验中的试样均为非饱和状态,采用最大干密度试验确定的含水率4.6%。一般高堆石坝堆石料的相对密实度为0.9左右,采用0.85,0.90和0.95三种相对密实度分别进行试验,每个相对密实度的试样都在0.4,0.8,1.6,3.2 MPa 4级轴向压力下进行单轴蠕变试验。

    不同相对密实度单轴蠕变试验的轴向蠕变应变随时间变化的曲线见图1,图中最右侧的数值为轴向压力,压力相同的一组曲线中上、中、下分别为相对密实度0.85,0.90,0.95的曲线。在蠕变初期,蠕变应变曲线斜率较大,蠕变速率较快;随着时间增长,蠕变速率逐渐减慢,最终蠕变应变趋于一个稳定值。这是由于在蠕变初期,堆石料的蠕变应变主要由颗粒破碎和颗粒重排列共同产生。随着轴向压力作用时间延长,颗粒破碎速率逐渐减慢,直至不再发生破碎。说明在试验压力范围内,含软岩堆石料的蠕变为衰减蠕变。

    图  1  不同相对密实度的轴向蠕变应变随时间变化
    Figure  1.  Change in axial creep strain with time under different relative densities

    随着相对密实度增大,试样的轴向蠕变应变逐渐减小,说明堆石料蠕变与其自身密实程度相关。密实程度低的堆石料孔隙率大,外力作用下,粗大颗粒棱角或软岩颗粒发生破碎、细化后,更容易滑移填充到孔隙中,产生较大蠕变应变。随着轴向压力增大,相同压力、不同相对密实度试样的蠕变应变曲线的差距逐渐增大。这是由于轴向压力较小时,颗粒表面接触应力较小,破碎颗粒数量少且粒径大,破碎颗粒填充受孔隙的影响较小,故不同密实程度的堆石料蠕变应变相差较小。随着轴向压力增大,颗粒表面接触应力增大,破碎颗粒数量增加且粒径变小,破碎颗粒滑移填充受孔隙的影响程度提高,故各试样间蠕变应变曲线差距变得明显。轴向蠕变应变曲线局部存在交叉现象,主要是由于堆石料颗粒互相咬合,若局部发生破碎或错动,变形量会有一定的激增,在蠕变应变曲线中表现较明显。

    长江科学院提出的幂函数蠕变模型已成功应用于水布垭、双江口[12]的堆石料蠕变研究。说明该模型用于表达堆石料的蠕变特性具有较好的适用性,故本文的试验结果也采用该模型进行拟合。该蠕变模型为

    εc=εf(1tλ), (1)

    式中,εc为轴向蠕变应变,εf为最终蠕变应变,λ为衰减指数,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最终蠕变应变εf/%衰减指数λ决定系数R2
    0.850.40.2950.0900.9712
    0.80.5740.1070.9638
    1.60.9380.1210.9826
    3.21.4410.1310.9729
    0.900.40.2680.0890.9866
    0.80.5350.1020.9844
    1.60.8370.1210.9107
    3.21.3470.1300.8501
    0.950.40.2370.0870.9621
    0.80.4880.0990.9870
    1.60.7310.1200.9632
    3.21.2250.1290.9829
    下载: 导出CSV 
    | 显示表格

    利用表2的参数计算各相对密实度和轴向压力的轴向蠕变应变曲线,其结果与试验结果的比较见图2。由图2可知,曲线与离散点偏差小,变化趋势一致,说明幂函数蠕变模型能较好表现含软岩堆石料蠕变特性。

    图  2  幂函数蠕变模型计算的轴向蠕变应变和试验结果的对比
    Figure  2.  Comparison between axial creep strain calculated by creep model of power function and test results

    通过绘制最终蠕变应变随轴向压力变化的曲线可知,两者关系曲线形状近似双曲线:

    εf=pap+b, (2)

    式中,p为轴向压力(MPa),a,b为拟合参数。故用双曲线来拟合两者关系,见图3,图例中T代表试验点,F代表拟合曲线。衰减指数随轴向压力的变化关系见图4,由图4可知,两者关系曲线形状近似幂函数:

    λ=cpd, (3)
    图  3  最终蠕变应变-轴向压力关系
    Figure  3.  Relationship between the final creep strain and axial pressure
    图  4  衰减指数-轴向压力关系
    Figure  4.  Relationship between attenuation index and axial pressure

    式中,cd为拟合参数,拟合曲线也绘于图4中。

    将式(2),(3)代入式(1),整理后可得考虑轴向压力p的单轴蠕变模型表达式:

    εc=pap+b(1tcpd) (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最终蠕变应变εf衰减指数λ
    abR2cdR2
    0.850.3301.1670.99940.1090.1740.9721
    0.900.3341.3190.99730.1070.1840.9730
    0.950.3491.5170.99340.1050.1930.9687
    平均值0.3381.3340.1070.184
    标准差0.0080.1430.0020.008
    下载: 导出CSV 
    | 显示表格

    公式(4)中的各参数的标准差见表3最后一行。参数a,c,d的标准差均小于0.01,远小于参数b的标准差0.143,说明相对密实度对这三个参数的影响很小。为简化模型,将参数a,c,d取平均值,分别为0.338,0.107,0.184。由表2中不同相对密实度、相同轴向压力的衰减指数λ值的对比,也可看出相对密实度对衰减指数的影响很小。

    表3可知参数b与相对密实度呈正相关,据此可进一步分析参数b与相对密实度的函数关系,建立考虑相对密实度的蠕变模型。将不同相对密实度的参数b绘于图5中,由图可知,其与相对密实度呈良好的线性关系,拟合的决定系数大于0.95。将图5中参数b与相对密实度的线性函数关系代入式(4),可得

    εc=pap+(b1Dr+b2)(1tcpd), (5)
    图  5  参数b与相对密实度的关系
    Figure  5.  Relationship between parameter b and relative density

    式中,a,b1,b2,c,d为模型参数,各参数取值见表4

    表  4  考虑轴向压力和相对密实度的蠕变模型参数
    Table  4.  Parameters of creep model for rockfill considering axial pressure and relative density
    模型参数ab1b2cd
    取值0.3383.50-1.8160.1070.184
    下载: 导出CSV 
    | 显示表格

    利用考虑轴向压力和相对密实度的蠕变模型计算各轴向压力和相对密实度含软岩堆石料的轴向蠕变应变曲线,计算结果和试验结果的比较见图6。由图6可知,计算结果与试验结果偏差较小,且变化趋势一致。说明前文中求取模型参数的方法可行,模型参数可靠。

    图  6  最终蠕变模型计算的轴向蠕变应变与试验结果的对比
    Figure  6.  Comparison between axial creep strain calculated by final creep model and test results

    对于以相对密实度为填筑指标的堆石坝,可进行单轴蠕变试验获得堆石料轴向蠕变应变曲线。利用式(1)拟合求取各相对密实度堆石料最终蠕变应变和衰减指数,然后分别用式(2),(3)拟合两者随轴向压力变化的曲线,得到参数a,b,c,d。将参数a,c,d取平均值,将参数b表示为相对密实度的线性函数,得到参数b1b2。为符合堆石料变形特性,参数b不能为负值,因此其与相对密实度的线性关系仅适用于密实情况下的堆石料。本文建立的蠕变模型在前文讨论的压力和相对密实度范围对含软岩堆石料适用,上述范围以外模型的适用性需进一步研究确定。单轴蠕变模型能较好表现堆石料的蠕变特性,预测堆石坝的长期变形特性。如梁军[4]利用单轴蠕变模型结合广义虎克定律,计算得到的关门山坝的沉降曲线与实测沉降过程线相似程度较高。

    本文的模型是通过将式(1)中最终蠕变应变和衰减指数表示为轴向压力和相对密实度的函数关系建立的,因此模型参数的物理意义可通过最终蠕变应变和衰减指数表达式中各参数的物理意义来说明。由式(2)可知,1/a为轴向压力趋近于无穷大时最终蠕变应变。对式(2)求导,1/b为轴向压力为0时最终蠕变应变随轴向压力增长的速率。由图5可知,b1为参数b随相对密实度的变化速率,b2为堆石料相对密实度为0时的b值。由式(3)为幂函数可知,c为轴向压力为1 MPa时的衰减指数,d是双对数坐标下衰减指数与轴向压力关系曲线的斜率,表征衰减指数随轴向压力变化的速率。

    本文对不同相对密实度的含软岩堆石料进行了单轴蠕变试验研究,得到以下结论:

    (1)在试验压力范围内,含软岩堆石料的蠕变类型为衰减蠕变。随着相对密实度的增大,堆石料的轴向蠕变应变逐渐减小。轴向压力较低时,不同相对密实度的轴向蠕变应变差距较小;随着轴向压力提高,轴向蠕变应变的差距逐渐增大。

    (2)利用幂函数蠕变模型拟合轴向蠕变应变随时间变化的曲线,得到最终蠕变应变和衰减指数。其分别与轴向压力满足双曲线和幂函数关系,由此建立了考虑轴向压力的含软岩堆石料的单轴蠕变模型。

    (3)分析模型参数可知相对密实度对衰减指数的影响很小。依据与最终蠕变应变相关的模型参数b和相对密实度呈线性关系,提出了同时考虑轴向压力和相对密实度的单轴蠕变模型。利用分析确定的模型参数计算各轴向压力和相对密实度情况下堆石料的轴向蠕变应变,计算结果与试验结果较吻合。说明提出的单轴蠕变模型在本文讨论的压力和相对密实度范围对含软岩堆石料是适用的。

  • 图  1   光缆结构示意图

    Figure  1.   Structural diagram of optical cable

    图  2   室内试验装置构成图

    Figure  2.   Composition of laboratory test unit

    图  3   1号光纤在不同环境下温升对比图

    Figure  3.   Comparison of temperature rises of fiber No. 1 under different environments

    图  4   2号光纤在不同环境下温升对比图

    Figure  4.   Comparison of temperature rises of fiber No. 2 under different environments

    图  5   非运行期时新疆某输水渠道图片

    Figure  5.   Picture of a water transfer in Xinjiang during non-operational period

    图  6   试验段光缆布置剖面图

    Figure  6.   Cross-section of cable layout during field tests

    图  7   渗漏光缆和孔隙水压力计布设位置示意图

    Figure  7.   Schematic diagram of placement of pore water pressure gauge and cable

    图  8   渗漏光缆和孔隙水压力计埋设现场

    Figure  8.   Field photos in test area

    图  9   点加热定位技术测量数据

    Figure  9.   Measured data by using point heating location technology

    图  10   典型时间点K值曲线图

    Figure  10.   Values of K at typical time points

    表  1   光纤测渗漏现场测值分析表

    Table  1   Analysis of measured values of optical fiber leakage

    日期光纤在水位以上的长度/m水面距离孔压的高度/m孔压相应的测值/m
    2019-07-22 温差法56.5~663.253.175
    85.5~93.53.61
    116.5~122.54.03
    144.5~152.53.62
    2019-07-22 K值法56~65.53.25
    84.5~933.50
    116.5~122.54.03
    144.5~1523.73
    2019-09-19 温差法50.5~701.020.685
    82~101.51.04
    109~130.50.57
    140~157.51.50
    2019-09-19 K值法51~701.13
    82.5~1011.26
    111~1301.13
    141.5~160.51.16
    下载: 导出CSV
  • [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)

图(10)  /  表(1)
计量
  • 文章访问数:  168
  • HTML全文浏览量:  25
  • PDF下载量:  78
  • 被引次数: 8
出版历程
  • 收稿日期:  2020-08-06
  • 网络出版日期:  2022-12-07
  • 刊出日期:  2020-10-31

目录

/

返回文章
返回