Development and initial application of test devices for leakage of earth-rockfill dams
-
摘要: 为准确模拟渗漏病险发生条件下土石堤坝渗流特性，研制一种土石堤坝渗漏病险试验装置。该装置能够精准控制试验槽坡度、输水流量、库水温度等试验变量，可实现边壁及底板平整度的自动校核、流控系统与温控系统的自动化控制、多种供吸水组合方式的水循环、库水的分级与线性升温和降温、坝体温度场和体积变形以及孔隙水压力的动态监测等功能。选用黏土心墙堆石坝为研究对象，初步开展了土石堤坝渗漏病险试验研究。研究表明：渗漏病险区将改变坝体的流热双场，其周围的温度变幅大于孔压变幅；随着渗漏病险区逐步下移，所在位置处的温度值显著增大，浸润线的波动情况与坝体渗流速率也随之增大，影响坝体渗流稳定性的不利情况逐渐增强。该装置能够准确模拟各工况真实情况以及多场耦合效应对土石堤坝渗漏病险发生发展过程的影响，能够实时监测渗漏病险发展过程中各物理场演变情况，为深入研究土石堤坝渗流特性提供一种有效的技术工具。Abstract: To accurately simulate the seepage characteristics of earth-rockfill dams during leakage, a test device is developed. The device can precisely control the test variables such as platform slope, water supply discharge, water temperature, etc. The device can realize the automatic check of the flatness of the side wall and the bottom plate, the automatic control of the flow control system and the temperature control system, the water circulation of various water supply and suction combinations, and piecewise or linear temperature rise and drop control of the upstream reservoir water, dynamic monitoring of dam body temperature field, volume deformation and pore water pressure. A clay core wall rockfill dam is selected as the research object, and the leakage tests on earth-rockfill dams are conducted. The results show that the leakage will change the seepage field and temperature field of the dam body, causing the surrounding temperature value to change significantly. With the gradual downward movement of the leakage area, the temperature value at the location increases significantly, the fluctuation of the saturation line and the seepage rate of the dam also increase, and the adverse conditions affecting the seepage stability of the dam body gradually increase. The test device can accurately simulate the real situation of various working conditions and the influences of multiphysics coupling effects on the leakage occurrence and development process of earth-rockfill dams, and it can monitor the evolution of various physical fields in the leakage development process in real time, providing an effective technical tool for the in-depth study on the seepage characteristics of earth-rockfill dams.
-
Keywords:
- earth-rockfill dam /
- leakage /
- test device /
- test method /
- temperature field /
- seepage field
-
0. 引言
近年来地下综合管廊逐渐被应用在城市地下空间布局中,改变了以往地下管线的布置方式。过去人们总是认为地下结构具有良好的抗震性能,但是自1995年日本阪神地震以后,这一观点开始被打破。在后来的几次大地震中,也出现过地下结构破坏的报道,比如1999年台湾集集地震、2008年汶川地震和2011年日本东北大地震等[1]。鉴于此,20多年来不断有国内外的学者开始在地下结构抗震方面开展一些研究工作[2-8]。
有学者试图从理论解析的角度来研究地下结构的抗震性能,比如Gil等[3]和Tateishi[5],其中Tateishi在研究过程中,假设结构和近场的加速度响应与自由场情况下在相同位置是一样的,但这有待考究。振动台试验目前在地震动研究中被广泛应用[9-11],有学者单从数值角度来进行研究[1, 12],也有学者将数值技术与试验结合起来,探讨地下结构的地震响应特征[13],但也有少数学者开展离心机地震模拟试验研究[14]。有学者从振动台试验方法、试验结果以及数值模拟方面,对地下管廊做了大量研究[15-17]。另外,很多学者对振动台试验装备与技术,特别是对叠层式土箱,提出了很多宝贵意见[18-20]。总结之前的研究,主要存在以下问题：①大多只研究单仓结构形式的地下管廊；②弯矩响应分析不足；③部分研究人员只进行了数值分析,缺乏试验和实测数据支持。
本研究采用重庆大学小型振动台模拟系统,开展了双仓矩形截面地下管廊的抗震性能研究,主要研究地震作用下管廊结构和周围土体产生的地震响应规律。根据刘晶波等[21]的研究,选取结构的最不利时刻对管廊的弯矩特征进行了研究,分析过程简单明了,结果可靠。同时,利用数值软件ABAQUS开展数值模拟研究,对试验结果进行对比检验,表明试验结果可靠。
1. 试验设计及传感器布置
图1展示的是小型振动台试验系统,台面尺寸为1.2 m×1.2 m,配备叠层式钢制剪切土箱,以模拟地震作用下岩层的剪切效应,并消除边界处复杂的边界效应[18, 20],其最大载重量为10 kN,可实现水平与垂直两个方向的振动,10 kN载重下振动加速度可达到1.2g,5 kN时为2.0g。
模拟双仓管廊的截面尺寸为6.7 m×4.5 m,需根据相似模拟试验研究[15],对试验模型的几何尺寸、物理力学性质以及时间参数的相似比进行设计,如表1所示。根据1∶15的几何相似比,使用特定水灰比配置的微粒混凝土和镀锌铁丝制成的双仓管廊模型尺寸为80 cm×44.67 cm×30 cm（长×宽×高）,实测该微粒混凝土的弹性模量为6062 MPa。应实际工程要求,采用福建厦门的标准砂近似模拟砂土场地（或者说Ⅱ、Ⅲ类场地）,试验中通过分层夯实处理后,模拟砂土场地的密度近似为1740 kg/m3。但考虑到试验条件,未能同时满足所有条件,只能根据实际情况抓主要因素,忽略次要因素。最终完整的试验模型如图2所示。
表 1 相似比设计Table 1. Design of similarity ratios变量 计算公式 相似比 几何Sl 选定 1∶15 应变Sε 选定 1∶1 应力Sσ 选定 1∶3 弹性模量SE SE= Sσ/Sε 1∶3 加速度Sa 选定 5∶1 密度Sρ Sρ= SE/(Sl·Sa) 1∶1 时间St St= √Sl/Sa 1∶8.67 配筋 按照配筋率一致 根据谢礼立等[22]和翟长海等[23]的研究成果,1952年美国加利福利亚地震的Taft地震波比较适合Ⅱ、Ⅲ类场地的抗震试验研究,其时程曲线和傅里叶谱频谱如图3所示,持时为54.00 s,峰值加速度为0.016g。为充分考虑不同PGA的情况,将原始地震波的加速度峰值进行调整,生成场地峰值加速度PGA分别为0.2g,0.4g,0.8g和1.2g的地震波,研究不同PGA的影响规律。试验中振动波从试验系统底部向上输入,其振动方向垂直于管廊结构的延长方向。
![]()
图 3 原始Taft地震波时程曲线及其傅里叶谱Figure 3. Time-histories and Fourier spectra of original acceleration of Taft earthquake spectra考虑到结构的对称性,以管廊截面纵轴线为对称轴,土压力计和加速度计各布置一边,如图4所示。为研究管廊结构的弯矩响应规律,在管廊中部截面外壁均匀布置了一系列应变片,如图5所示。
2. 试验结果与讨论
2.1 最大土压力响应
试验中土压力计测出的值是土压力增量,正值表示动土压力相对于初始土压力增大。如图6所示,不同输入PGA（0.2g,0.4g,0.8g和1.2g）时,最大土压力响应沿深度总是呈现倒立的“W”形分布,这种分布特征可能是由于振动过程中管廊结构与土体之间相互作用,引起了水平向的土拱效应所致。随着输入PGA的增大,各个测点的最大土压力响应也不断增大,只是在1.2g时,Taft地震波下管廊结构底部土压力响应略小于0.8g下的响应,说明随着输入PGA的增大,PGA的影响相较于管廊结构逐渐减小。
![]()
图 6 管廊侧壁邻近土体沿深度最大动土压力响应Figure 6. Maximum dynamic earth pressure responses along depth in soils adjacent to side wall图7反映了在Taft地震波下,不同输入PGA时（0.2g,0.4g,0.8g和1.2g）,管廊侧壁沿深度从上到下3个土压力测点（P3、P4和P5）测出的土压力响应时程曲线。可以看出,地震过程中,管廊底部的土压力响应一般大于其他两点的响应,中部最小,特别是在最大土压力响应发生的时刻,这进一步验证了之前关于土拱效应的推断。且在地震过后,结构周围土体的应力场相比较于地震发生前已经发生了变化。
![]()
图 7 管廊侧壁周边土压力响应时程曲线Figure 7. Time-histories of dynamic earth pressure for surrounding soils adjacent to side wall由表2可知,相比较于初始土压力场,地震过后是减小还是增大,没有特别明显的规律,但在管廊中下部,往往会在震后变小。这说明在地震过后,管廊结构与周围土体可能出现局部脱离、局部挤密的现象,数据显示在结构中下部易出现局部脱离,而在中上部会易出现局部挤密,然而这引起的是与静态作用下截然相反的土压力分布模式。震后土压力相比于最大动土压力的作用小得多,对于顶部P3测点,其(Pr/Pdmax)×100%的值为1.38%～9.15%；对于中部P4测点,其(Pr/Pdmax)×100%的值为3.16%～18.87%；底部P5测点的值为1.21%～4.14%。
表 2 振动过程中最大动土压力与震后土压力的比较Table 2. Comparison between maximum dynamic earth pressure and earth pressure at end of earthquake测点 最大动土压力响应Pdmax/Pa 震后土压力Pr/Pa (Pr/Pdmax)/% P3 P4 P5 P3 P4 P5 P3 P4 P5 Taft-0.2g 1942 906 5700 54 -29 -69 2.77 3.16 1.21 Taft-0.4g 1942 906 5700 83 104 -236 4.25 11.49 4.14 Taft-0.8g 5395 1736 26735 -494 209 -322 9.15 12.07 1.21 Taft-1.2g 6411 3472 26197 89 655 826 1.38 18.87 3.15 2.2 加速度响应
由于场地性质、地层深度、结构与岩土体的材料性质等差异,同一地震作用下不同地点的地震响应也会有较为明显的不同[24]。图8反映了在Taft地震波下,不同输入PGA时（0.2g,0.4g,0.8g和1.2g）,管廊侧壁邻近土体沿深度的加速度响应情况。从整体来看,不同输入PGA下最大加速度响应均沿深度逐渐减小,各深度处随输入PGA的增大而增大。结合管廊结构在土体中的位置分析,最大加速响应随着深度先减小,然后在管廊侧壁中下部转而逐渐增大,这可能是由土与结构相互作用所引起的。
![]()
图 8 管廊侧壁邻近沿深度土体最大加速度响应Figure 8. Maximum acceleration responses of soils adjacent to side wall along depth同时,根据最大加速度沿深度的发展特征可知,管廊侧壁土体的加速度响应沿深度可以分3个区段,一区段是管廊结构顶部以上土体,二区段是邻近管廊侧壁土体,三区段是管廊结构底部以下土体,而且随着输入PGA的增大,该区段的划分越发明显,也即结构对于周边土体的地震影响影响愈发突出。
如图9所示,当输入PGA在0.2g～1.2g变化时,加速度放大系数在0.5～1.5变化,沿深度的分布规律与图10中的最大加速度响应类似,且加速度放大系数随着输入PGA的增大而减小。所以,当输入峰值加速度加大,所处位置较深时,加速度放大系数多数小于1.0。
![]()
图 9 加速度放大系数在不同地震波型下沿深度的分布Figure 9. Distribution of acceleration amplification of the surrounding soils adjacent to side wall along depth图10清楚反映了在Taft地震波下,加速度放大系数在不同输入PGA和深度所构成的域中的分布特征—在地表处以及输入PGA较大时,加速度放大系数可能会有最大值。例如PGA为0.2g而深度为7.5 cm时,加速度放大系数为1.44；PGA为1.2g而深度为52.5 cm时,为0.56。此发现与Yan等[9]、Chen等[6]的研究结论一致。而地震来临时,地震波的PGA已经确定,唯一能够控制地下结构加速度响应的只有结构埋深。因此在城市地下综合管廊建设时,应当尽可能的深埋,以保证地下设施在地震作用下不被破坏。
为了研究在地震过程中管廊结构与周围土体的加速度响应差异,在管廊侧壁内壁对应于外壁土体中的加速度测点（A8、A15和A11）,相应布置了3个加速度测点（A10、A13和A5）。不同输入PGA下（0.2g,0.4g,0.8g）,上述6个测点的加速度响应时程曲线和傅里叶谱如图11所示。
![]()
图 11 土体和管廊结构加速度响应时程曲线与傅里叶谱对比Figure 11. Comparison of time histories of acceleration and Fourier spectra between structures and surrounding soils从时程曲线来看,在各个输入PGA的Taft地震波下,管廊侧壁与周围土体的加速度时程曲线总是基本保持一致,史晓军等[25]通过振动台试验也得出了类似的结论。另外,由图11可以看出,随着输入峰值加速度PGA的增加而加速度响应不断增大。从上述现象说明,在地震过程中,管廊结构与周围土体还是保持基本一致的运动模式,二者之间没有发生明显的脱离现象。
由傅里叶谱可知,在同一PGA地震波下,管廊结构及其周围土体的振幅在相同频段处也基本保持一致—同一地震波下最大振幅随着深度稍减小,同一位置随输入PGA的增大而增大。只是在15～30 Hz时,土体中的振幅稍大于结构的振幅,尤其是在管廊结构顶部,此现象更为明显。这是由于土体和管廊结构的自然频率不同,同时也说明了地震波的频率成分也会对结构和土体的地震响应造成差异。
2.3 弯矩响应
根据应变片测点数据,在同一模拟地震作用下,提取应变片测点测出的最大应变响应,以此响应对应时刻为本次地震的最不利时刻,计算出本次模拟地震下管廊结构的弯矩响应,最终得出了双仓管廊结构在地震作用下的弯矩响应情况,如图12所示。由结构周围的弯矩包络图可知,在地震波作用下,管廊结构中拐角处弯矩较大,而各侧边中点的弯矩就小得多,几乎为0。另外,随着输入PGA的增大,弯矩包络图的范围也随之增大,转角处测点在不同PGA下的弯矩数据如表3所示。
![]()
图 12 地震作用下管廊结构弯矩响应Figure 12. Responses of bending moment of utility tunnel structures under earthquake表 3 不同PGA的Taft地震波下管廊角点处弯矩值Table 3. Values of bending moment at corners under Taft earthquake with various PGAs(N·m) 测点编号 0.2g 0.4g 0.8g 1.2g S13 -4.74 -10.80 -19.72 -23.02 S15 4.48 9.46 15.36 17.58 S16 5.02 11.38 20.11 22.89 S18 -5.52 -12.15 -20.72 -23.59 S19 -5.71 -13.14 -21.19 -23.83 S21 5.90 12.65 21.47 24.50 S22 7.72 15.98 25.71 28.92 S24 -8.26 -17.35 -27.71 -30.95 注： 负号表示结构的内侧以及内隔墙的右侧受拉。对弯矩数据进一步分析,表4给出了管廊结构在Taft地震波下,随着PGA的增大弯矩响应增量,由表4可知,4个拐角处的增量随输入PGA的增大而逐渐减小。4个拐角在输入峰值加速度分别从0.2g增大到0.4g,0.4g增大到0.8g以及0.8g增大到1.2g时,其4个拐角处的弯矩增幅依次为118.46%,67.99%以及13.69%。
表 4 Taft地震波下管廊4个角点弯矩随PGA增大的增量Table 4. Incremental values of bending moment at corners with increasing PGAs under Taft earthquake(%) PGA S13 S15 S16 S18 S19 S21 S22 S24 均值 0.2g ~0.4g 127.77 111.17 126.90 120.18 129.96 114.55 107.12 110.04 118.46 0.4g ~0.8g 82.65 62.32 76.70 70.60 61.33 69.73 60.87 59.69 67.99 0.8g ~1.2g 16.75 14.43 13.79 13.83 12.45 14.08 12.48 11.69 13.69 由此说明当震级越来越大,PGA越来越高时,结构中的弯矩响应不会无限增大。这是因为当地震引起的地面运动越剧烈时,结构随着周围土体整体运动的趋势越发明显,此时管廊结构的破坏更主要的可能是在延长方向上,由于地层差异引起管廊结构在延长方向上的纵向弯矩过大,从而导致地下管廊在某一位置发生断裂。
3. 数值模拟验证
3.1 数值模型的建立
考虑到管廊属于无限长结构,可以应用平面应变模型进行分析,结合本次试验目的,本次数值模型为2D模型。数值模拟中,模拟剪切波引起地层产生剪切变形是至关重要的一步。根据Chen等[26]的研究,控制模型左右边界在任意深度处的位移自由度相一致,即可实现上述要求。同时,考虑到模拟的方便性,直接模拟原型管廊,综合管廊截面尺寸为6.7 m×4.5 m,其埋深为2.5 m,壁厚为0.3 m,模拟土体范围为48 m×28 m。地震波选取Taft原始波,但将其PGA调整为0.2g,0.4g,0.8g和1.2g,地震波从模型底部垂直向上输入。整体数值模拟网格划分模型及地震输入方式如图13所示。
![]()
图 13 双仓地下综合管廊数值模拟网格划分模型Figure 13. Meshing model for numerical simulation of double-box utility tunnel数值模拟中应用的材料参数如表5所示,包括密度
ρ 、弹性模量E、泊松比ν 、黏聚力c、内摩擦角φ和阻尼比D。其中土体弹性模量E由于和剪切模量G存在理论关系,而动剪切模量G与动剪应变γ 有关,所以30 MPa是土体在振动前的初始弹性模量,振动模拟过程中随动剪切模量的变化而变。本次模拟中应用的动剪切模量G和阻尼比D与动剪应变γ 的关系如图14所示。对于土体,应用莫尔-库仑理想弹塑性的本构模型,对于管廊结构应用混凝土塑性损伤模型。表 5 数值模拟材料参数Table 5. Parameters of material in numerical simulation材料 ρ/(kg·m-3) E/MPa ν c/kPa φ/(°) D 土体 1740 30(初始) 0.3 2000 35 — 管廊 2400 31500 0.2 — — 0.05 ![]()
图 14 动剪切模量和阻尼比与动剪应变的关系曲线Figure 14. Relationship among dynamic shearing modulus, damping ratio and dynamic shearing strain3.2 数值模拟结果分析
提取不同工况下管廊侧壁邻近土体（如图13中左侧虚线所示）的加速度响应数据,发现加速度响应沿深度与试验结果有非常类似的规律,如图15所示。首先从整体上来看,水平加速度峰值沿深度逐渐减小,各深度的水平加速度峰值随输入PGA的增大而增大。结合管廊结构在土体中的位置来看（如图15中虚线框所示）,也出现了从管廊顶部先减小、在管廊中下部又转而增大的现象。
![]()
图 15 管廊侧壁邻近土体水平加速度响应Figure 15. Maximum acceleration responses of soils adjacent to side wall of utility tunnel图16反映了结构侧壁邻近土体（如图13中右侧黄色虚线所示）的最大水平向土压力响应。从图15中可以看出,最大土压力响应在管廊结构附近呈现了非常明显的倒立的“W”形分布,这充分说明了在地震动下,管廊结构侧壁与邻近土体相互作用形成了土拱效应,这也非常符合之前由试验结果做出的猜想,同时也进一步说明前述“三区段”分析的正确性。随着输入PGA的增大,最大土压力响应也随之增大,但增加幅度在管廊结构周围逐渐减小。此外,由于本次数值研究的土体深度较试验相对更深,所以反映出了另一重要特征,即在远离结构的深部土体,其最大土压力并不是逐渐减小或逐渐增大的,而是在从管廊结构先开始增大,而后在距管廊底部约三倍管廊高度的深度位置逐渐减小。
![]()
图 16 管廊侧壁邻近土体水平土压力响应Figure 16. Responses of maximum horizontal earth pressure of soils adjacent to side wall of utility tunnel由图15,16给出的加速度响应和土压力响应特征来看,数值模拟的结果与试验的结果在规律上的有很高的吻合度,由此证明了试验结果的可信性。
4. 结论
基于相似理论,利用振动台试验系统模拟研究了在不同地震波下,管廊结构及其周围土体的物理力学以及运动响应特征,同时结合数值模拟对试验结果进行了对比,由此主要得到如下4点结论。
（1）管廊侧壁沿深度最大土压力响应呈现倒立的“W”形分布,这可能是振动引起土拱效应所致,且经过地震作用后,结构周围土体的土压力场也发生了改变。
（2）管廊侧壁最大加速度响应沿深度先随深度减小,在管廊中下部转而增大。加速度放大系数随着输入PGA增大而减小,PGA为0.2g～1.2g时,其值在0.5～1.5变化。在地震过程中,结构及其周围土体基本保持一致的运动模式。
（3）在地震荷载作用下,在结构拐角处会产生较大的弯矩响应,随着输入PGA的增大而增大。
（4）随着输入PGA的增大,管廊结构对于周围土体的土压力以及加速度响应的影响越发明显,但管廊结构自身弯矩则不会无限增大。
-
![]()
图 7 光纤光栅温度计与光纤测温结果对比图
Figure 7. Comparison of results of fiber grating thermometer and fiber temperature measurement
![]()
图 8 各工况坝体温度分布云图
Figure 8. Cloud map of temperature distribution of dam under various working conditions
![]()
图 9 不同工况下温度监测结果
Figure 9. Monitoring results of temperature under different working conditions
![]()
图 10 各工况心墙中心线监测结果
Figure 10. Monitoring results of center line of core wall under various working conditions
![]()
图 11 不同工况下孔隙水压力监测结果
Figure 11. Monitoring results of pore water pressure under different working conditions
![]()
图 12 不同工况下坝体单宽渗漏速率时序曲线
Figure 12. Variation of single-width leakage rate of dam body with time under different working conditions
![]()
图 13 各工况浸润线监测结果
Figure 13. Monitoring results of infiltration line under various working conditions
表 1 各分区土料参数
Table 1 Parameters of soils in each area of dam body
坝体分区 渗透系数Κ/(m·s-1) 密度
ρ/(g·cm-3)含水率w/% 堆石区 4.71×10-2 1.69 1.99 过滤区 2.83×10-5 1.79 6.87 病险区 2.83×10-5 1.79 6.87 心墙区 5.16×10-7 2.03 8.81 
下载: 导出CSV
表 2 各工况设置详情
Table 2 Setting details of each working condition
工况 入渗水温/
℃入渗水头/
cm距坝基高度/cm 备注 1 28 55 — 无病险 2 28 55 40 高位 3 28 55 30 中位 4 28 55 15 低位
下载: 导出CSV
-
[1] SALMASI F, NOROUZI R, ABRAHAM J, et al. Effect of inclined clay core on embankment dam seepage and stability through LEM and FEM[J]. Geotechnical and Geological Engineering, 2020, 38(6): 6571-6586. doi: 10.1007/s10706-020-01455-7
[2] 苏怀智, 周仁练. 土石堤坝渗漏病险探测模式和方法研究进展[J]. 水利水电科技进展, 2022, 42(1): 1-10, 39. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD202201001.htm SU Huaizhi, ZHOU Renlian. Research progress and prospect of earth-rockfill dam leakage detection modes and method[J]. Advances in Science and Technology of Water Resources, 2022, 42(1): 1-10, 39. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD202201001.htm
[3] SCHMITZ V, ERPICUM S, EL KADI ABDERREZZAK K, et al. Overtopping-induced failure of non-cohesive homogeneous fluvial dikes: effect of dike geometry on breach discharge and widening[J]. Water Resources Research, 2021, 57(7): 029660.
[4] 钟登华, 时梦楠, 崔博, 等. 大坝智能建设研究进展[J]. 水利学报, 2019, 50(1): 38-52, 61. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201901006.htm ZHONG Denghua, SHI Mengnan, CUI Bo, et al. Research progress of the intelligent construction of dams[J]. Journal of Hydraulic Engineering, 2019, 50(1): 38-52, 61. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201901006.htm
[5] 彭铭, 毕竞超, 朱艳, 等. 存在高渗透区的黏土心墙土石坝渗流稳定性分析[J]. 水利学报, 2020, 51(11): 1347-1359. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202011005.htm PENG Ming, BI Jingchao, ZHU Yan, et al. Seepage stability analysis of earth and rockfill dam with clay core by considering high permeability zones[J]. Journal of Hydraulic Engineering, 2020, 51(11): 1347-1359. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202011005.htm
[6] DONG H Z, CHEN J S, LI X Y. Delineation of leakage pathways in an earth and rockfill dam using multi-tracer tests[J]. Engineering Geology, 2016, 212: 136-145. doi: 10.1016/j.enggeo.2016.08.003
[7] RATIAT A, KHETTAL T, MEDDI M. The piezometric and isotopic analysis of leaks in earth dams: the case of the fountain of Gazelle Dam, Biskra, Algeria[J]. Environmental Earth Sciences, 2020, 79(6): 138. doi: 10.1007/s12665-020-8886-8
[8] LUO Y L, NIE M, XIAO M. Flume-scale experiments on suffusion at bottom of cutoff wall in sandy gravel alluvium[J]. Canadian Geotechnical Journal, 2017, 54(12): 1716-1727. doi: 10.1139/cgj-2016-0248
[9] ROCHIM A, MAROT D, SIBILLE L, et al. Effects of hydraulic loading history on suffusion susceptibility of cohesionless soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2017, 143(7): 04017025. doi: 10.1061/(ASCE)GT.1943-5606.0001673
[10] NGUYEN C D, BENAHMED N, ANDÒ E, et al. Experimental investigation of microstructural changes in soils eroded by suffusion using X-ray tomography[J]. Acta Geotechnica, 2019, 14(3): 749-765. doi: 10.1007/s11440-019-00787-w
[11] 孔纲强, 孙学谨, 刘汉龙, 等. 孔隙液体对透明土渗透特性影响对比试验[J]. 水利学报, 2017, 48(11): 1303-1310. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201711005.htm KONG Gangqiang, SUN Xuejin, LIU Hanlong, et al. Contrast experiments on permeability of transparent soil influenced by pore fluids[J]. Journal of Hydraulic Engineering, 2017, 48(11): 1303-1310. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB201711005.htm
[12] 梁越, 陈鹏飞, 林加定, 等. 基于透明土技术的多孔介质孔隙流动特性研究[J]. 岩土工程学报, 2019, 41(7): 1361-1366. doi: 10.11779/CJGE201907022 LIANG Yue, CHEN Pengfei, LIN Jiading, et al. Pore flow characteristics of porous media based on transparent soil technology[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(7): 1361-1366. (in Chinese) doi: 10.11779/CJGE201907022
[13] 谷敬云, 罗玉龙, 张兴杰, 等. 基于平面激光诱导荧光的潜蚀可视化试验装置及其初步应用[J]. 岩石力学与工程学报, 2021, 40(6): 1287-1296. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202106019.htm GU Jingyun, LUO Yulong, ZHANG Xingjie, et al. A suffusion visualization apparatus based on planar laser induced fluorescence and the preliminary application[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(6): 1287-1296. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202106019.htm
[14] HUNTER R P, BOWMAN E T. Visualisation of seepage-induced suffusion and suffosion within internally erodible granular media[J]. Géotechnique, 2018, 68(10): 918-930.
[15] NAN S H, REN J, NI F, et al. Heat tracing of embankment dam leakage: laboratory experiments and 2D numerical modelling[J]. Journal of Hydrology, 2022, 608: 127663.
[16] 王新建, 潘纪顺. 堤坝多集中渗漏通道位置温度探测研究[J]. 岩土工程学报, 2010, 32(11): 1800-1805. http://www.cgejournal.com/cn/article/id/8471 WANG Xinjian, PAN Jishun. Location detection of concentrated-leakage passages in dam by groundwater temperature[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(11): 1800-1805. (in Chinese) http://www.cgejournal.com/cn/article/id/8471
[17] 谷艳昌, 王士军, 庞琼, 等. 土坝温度场反馈渗流场可行性研究[J]. 岩土工程学报, 2014, 36(9): 1721-1726. doi: 10.11779/CJGE201409020 GU Yanchang, WANG Shijun, PANG Qiong, et al. Feasibility of temperature field feeding back seepage field for earth dams[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(9): 1721-1726. (in Chinese) doi: 10.11779/CJGE201409020
[18] CUONG B Q, ZHOU Y H, ZHAO C J. Seepage evaluation in embankment dam based on short-term temperature observation and heat injection[J]. Electronic Journal of Geotechnical Engineering, 2016, 21(26): 10493-10506.
[19] GHAFOORI Y, MAČEK M, VIDMAR A, et al. Analysis of seepage in a laboratory scaled model using passive optical fiber distributed temperature sensor[J]. Water, 2020, 12(2): 367.
[20] 王家琛, 朱鸿鹄, 王静, 等. 基于主动加热光纤法的毛细阻滞入渗模型试验研究[J]. 岩土工程学报, 2021, 43(1): 147-155. doi: 10.11779/CJGE202101017 WANG Jiachen, ZHU Honghu, WANG Jing, et al. Laboratory model tests on capillary barrier infiltration using actively heated fiber optic method[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(1): 147-155. (in Chinese) doi: 10.11779/CJGE202101017
[21] SMITH M, KONRAD J M. Analysis of the annual thermal response of an earth dam for the assessment of the hydraulic conductivity of its compacted till core[J]. Canadian Geotechnical Journal, 2008, 45(2): 185-195.
[22] GHAFOORI Y, VIDMAR A, ŘÍHA J, et al. A review of measurement calibration and interpretation for seepage monitoring by optical fiber distributed temperature sensors[J]. Sensors, 2020, 20(19): 5696.
[23] 何斌, 徐剑飞, 何宁, 等. 分布式光纤传感技术在高面板堆石坝内部变形监测中的应用[J]. 岩土工程学报, 2023, 45(3): 627-633. doi: 10.11779/CJGE20220089 HE Bin, XU Jianfei, HE Ning, et al. Application of distributed optical fiber sensing technology in internal deformation monitoring of high face rockfill dam[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(3): 627-633. (in Chinese) doi: 10.11779/CJGE20220089
[24] SU H Z, TIAN S G, CUI S S, et al. Distributed optical fiber-based theoretical and empirical methods monitoring hydraulic engineering subjected to seepage velocity[J]. Optical Fiber Technology, 2016, 31: 111-125.
[25] SU H Z, CUI S S, WEN Z P, et al. Experimental study on distributed optical fiber heated-based seepage behavior identification in hydraulic engineering[J]. Heat and Mass Transfer, 2019, 55(2): 421-432.
[26] CHEN J A, CHENG F, XIONG F, et al. An experimental study: Fiber Bragg grating–hydrothermal cycling integration system for seepage monitoring of rockfill dams[J]. Structural Health Monitoring, 2017, 16(1): 50-61.
-
期刊类型引用(17)
1. 义扬，肖映雄，余科. 任意多边形骨料混凝土细观模型的建立与数值模拟. 山东大学学报(工学版). 2025(01): 97-107 . 
百度学术
2. 朱艳贵. 不等向固结时混合土剪切行为的离散元分析. 工业建筑. 2024(03): 167-173 . 百度学术
3. 李冠鹏，陈占扬，李洪博，武海荣，黄鹏，刘子墨，杜文韬. 基于三维细观的混凝土靶板侵彻仿真研究. 河南城建学院学报. 2024(06): 48-55 . 百度学术
4. 王泽华，李昺，邢磊，龚文平. 基于PFC~(3D)的滑坡与建筑物相互作用过程研究. 安全与环境工程. 2023(01): 107-118+191 . 百度学术
5. 杨忠平，刘浩宇，李进，李绪勇，刘新荣. 土石混合料–基岩接触面剪切力学特性及剪切带变形特征研究. 岩石力学与工程学报. 2023(02): 292-306 . 百度学术
6. 朱艳贵. 考虑粗颗粒真实形状的混合土不排水剪切特性离散元研究. 公路交通科技. 2023(02): 28-35 . 百度学术
7. 茹晓军. 考虑细观特征的土石混合体边坡破坏机制分析. 铁道建筑技术. 2023(04): 1-4+10 . 百度学术
8. 贾聿颉，李冬冬. 基于筛分试验的土石混合体三维数值建模方法研究. 华北水利水电大学学报(自然科学版). 2023(02): 97-103 . 百度学术
9. 张宏虎，白伟，孙明祥，邓涛. 闽东地区含砾花岗岩残积土的细观剪切特性研究. 水利与建筑工程学报. 2022(02): 48-53+191 . 百度学术
10. 吴尚杰. 粉土填料细观剪切特性的颗粒离散元数值模拟研究. 水利与建筑工程学报. 2022(03): 129-133+175 . 百度学术
11. 李泽闯，刘志斌，程培峰，张昊，蔡启源. 基于颗粒流方法的含粗粒滑带土宏细观力学特性. 辽宁工程技术大学学报(自然科学版). 2022(02): 121-129 . 百度学术
12. 胡峰，李志清，刘琪，胡瑞林. 土石混合体的剪应力波动和跌落行为机制. 水文地质工程地质. 2021(03): 90-101 . 百度学术
13. 杨忠平，李进，蒋源文，胡元鑫，赵亚龙. 含石率对土石混合体–基岩界面剪切力学特性的影响. 岩土工程学报. 2021(08): 1443-1452 . 本站查看
14. 金扬国. 交通荷载下碎石土回填路基沉降分析. 安徽建筑. 2021(11): 143-144 . 百度学术
15. 杨忠平，蒋源文，李诗琪，李进，胡元鑫. 土石混合体—基岩界面剪切力学特性试验研究. 岩土工程学报. 2020(10): 1947-1954 . 本站查看
16. 俞隽，孙洪浩，郑霄阳. 块石含量对土石混合体剪切力学特性的影响. 南通大学学报(自然科学版). 2020(03): 83-89 . 百度学术
17. 董辉，朱宪明，陈立，罗正东，蒋秀姿. 基于阻抗谱的土石体电阻相似性孔隙特征研究. 仪器仪表学报. 2020(11): 119-128 . 百度学术
其他类型引用(26)
-
其他相关附件
计量
- 文章访问数: 361
- HTML全文浏览量: 46
- PDF下载量: 107
- 被引次数: 43
