Sedimentation characteristics of gravel in discharge pipe of slurry shield tunneling and its impact on slurry
-
摘要: 为研究泥水盾构掘进风化岩层时排浆管内砾石的输送特性，基于CFD-DEM方法建立排浆管砾石输送模型，分析砾石输送过程中的速度分布、受力特征、运动形式，以及砾石床层对泥浆流动的影响规律，研究砾石尺寸、砾石形状、泥浆流速和泥浆相对质量密度对砾石沉积特性的影响。研究表明，砾石颗粒以沙丘状床层间歇式推移。砾石速度沿高度分层，表层砾石速度快，运动以碰撞和滚动为主；近壁层砾石速度慢，运动以滑移为主。近壁层和表层砾石间存在过渡层。砾石床层使管道内泥浆流速和压力局部增大。砾石床层速度随泥浆流速增大呈现线性增大。砾石床层高度随砾石尺寸增大而增大。泥浆相对质量密度越大，砾石受到曳力越大。砾石形状系数越小，输送难度越大。建立砾石沉积风险评价指标，通过传速系数α和增压系数β对砾石沉积风险评价，并提出针对性调控方法，为预防和处理砾石沉积堵塞提供参考。Abstract: The slurry shield tunneling in weathered rock layers is prone to gravel accumulation within the slurry pipe. In order to study the transportation characteristics of gravel in the slurry pipe, a model of gravel transportation in the slurry pipe of the shield tunneling was established based on the computational fluid dynamics-discrete element (CFD-DEM) coupling method. The rheological characteristics of the slurry were considered through the rheological tests, and the accuracy of the model is verified through the gravel sedimentation tests. The velocity distribution, force characteristics and motion forms of gravel during transportation as well as the influences of gravel bed deposition on the slurry were analyzed. The effects of gravel size, gravel shape, slurry flow rate and slurry density on the sedimentation characteristics of gravel are explored. This study indicates that gravel particles move intermittently in a sand dune-like bed. The gravel velocity is layered along the height, and the surface gravel velocity is fast, with collision and rolling as the main movements. The velocity of gravel in the near wall layer is slow, and its movement is mainly dominated by sliding. There is a transition layer between the near wall layer and the surface gravel layer. The gravel bed layer locally increases the flow velocity and pressure of the slurry in the pipeline. The velocity of the gravel bed increases linearly with the increase of the slurry flow rate. The height of the gravel bed increases with the size of the gravel. The higher the density of the slurry, the greater the drag force on the gravel. The smaller the gravel shape factor, the greater the difficulty of transportation. Finally, the risk assessment indices are established for gravel sedimentation, and the targeted control methods are proposed to provide reference for preventing and treating sedimentation blockage of gravel.
-
0. 引言
中国属多震国家,很多大型水利工程兴建在强震地区,在地震作用下,它们可能发生裂隙、下沉、坍塌甚至破坏,不仅将对工程本身造成损失,而且会带来巨大的次生灾害,后果十分严重,因此,有关地震方面的研究历来受到重视。
分析表明,目前研究地震破坏的方法大致可分为3类：一是数值模拟,其计算结果受计算参数和本构数学模型的影响很大；二是现场灾后调查,通过对地震现场破坏现象的调查,研究分析地震的破坏机理及抗震措施,由于实际地震的时空不确定性和复杂性,这类研究比较表观,不能重复,且无法取得地震全过程实测资料与数值分析结果进行相互验证；三是振动台试验方法,常规地面上的振动台模型试验由于无法模拟岩土材料的重力作用,因此多用于刚性材料的结构试验,而采用离心机振动台,则可以在原型应力条件下,在模型底部产生可控制的地震波,从而可以通过各种监测手段直接获得地震引起的岩土结构物的动力变形和稳定特性。因此土工离心机振动台被国内外岩土工程界认为是最有效的地震模拟试验手段,通过与数值模拟分析相互补充和验证,可以提高土工抗震研究水平,解决相关的岩土工程抗震问题[1-6]。
鉴于离心机振动台在抗震研究方面的突出优势,南京水利科学研究院在2002年自行研制了中国第一台NS-1型电液式土工离心机振动台[7],并利用该振动台研究了某混凝土面板堆石坝体的地震加速度反应特性、面板应变反应特性、坝体的变形特性等[4]。试验结果表明,该离心机振动台能满足岩土工程地震问题研究的要求,提供了一个非常有效的研究手段,可在工程实际中大力推广应用。
为了更好地满足科研和工程需要,南京水利科学研究院在此基础上研制了性能更可靠、技术指标更先进的NS-2型离心机振动台系统。本文主要介绍了该设备的设计关键技术和性能指标,给出其主要功能特点和技术参数。通过开展某沥青混凝土心墙砂砾石坝动态离心模型试验,初步验证了离心机振动台系统性能指标。
1. 离心机振动台系统概况
1.1 系统总体结构
NS-2型离心机振动台是由南京水利科学研究院和中国工程物理研究院联合研制,采用电液伺服液压驱动,实现水平单向振动（垂直于离心机转臂）,如图1所示。主要由机械系统、液压系统、伺服控制系统、动态数据采集系统及离心模型试验辅助设备等组成,其中液压系统是系统的动力机构,机械系统提供负载运动支承以及试验体安放平台,伺服控制系统完成地震波形的模拟加载控制,系统组成如图2所示。
1.2 主要技术参数
根据离心机振动台模型试验研究的国内现状和所配置的土工离心机运行情况,同时考虑到土木和水利工程的实际情况,设计采用如下技术参数：
a）振动方向：一维水平振动（垂臂向）；
b）驱动方式：液压（蓄能器）驱动；
c)最大振动加速度：amax=20g；最大速度：Vmax=0.5 m/s；最大振动幅值：Amax=5 mm；
d）最大振动时间：3 s
e）频率范围：20～200 Hz；
f）最大有效负载（含模型箱、模型、以及其它辅助装备）：Mmax=500 kg；
g）振动台控制：利用以太网通过滑环实现通讯；
h）激振波形：地震波、随机波；
i）最大离心加速度：80g。
2. 主要部件的结构
2.1 机械系统
机械结构包括：滚柱导轨、振动台面、剪切模型箱、以及激振系统安装底板。
安装底板作为激振系统的安装基础,向下连接离心机吊篮、向上固定系统结构；滚柱导轨作为运动支承安装在底板上；振动台面通过滑块安装在滚柱导轨上,模型及模型箱作为振动负载安装在振动台面上。
2.2 液压系统
液压系统包括：伺服液压缸一套、蓄能器组一套、21 MPa油源、以及离心机上的供油系统等。
液压缸安装在吊篮底板中轴线上,背负一台MOOG电液伺服阀为液压缸供油；在吊篮内布置一组小容量蓄能器作为伺服阀的进、回油及供油稳压配套设备,在离心机大臂上布置一组大容量蓄能器作为激振系统主油路的供油设备,重新配置21 MPa油源。如图3所示。
2.3 伺服控制系统
伺服控制系统包括：伺服驱动系统、上位控制系统、数据测量采集系统。伺服驱动系统布置在离心机上仪器仓内,包括伺服驱动器、运动控制器以及配套的供电设备,实现伺服液压缸的位移跟随控制,具有预设凸轮曲线跟踪、自动调零、自动标定等功能；上位控制系统放置在实验控制间,安装有地震波形再现软件一套,实现地震波加速度振动波形的信号处理、地震波响应测量及控制,并具有试验数据和试验曲线的实时测量、屏显、存储、数据管理等功能,如图4；试验数据测量采集系统由传感器及信号采集设备组成,为伺服控制、安全保护、试验数据分析提供支持。
3. 性能评价
3.1 试验介绍
在南京水利科学研究院NHRI400 g·t大型土工离心机和NS-2型离心机振动台上开展某沥青混凝土心墙砂砾石坝的动态离心模型试验,限于篇幅,这里只介绍坝体模型试验的部分结果。振动台和离心机总体布置见图5。
从目前来看即使用世界最大的离心机要等比尺模拟全部的坝体仍然是不可能的,因此,本项研究采用不等比尺的离心模型试验方法,即模型几何比尺ηl小于加速度比尺的倒数1/ηg,振动台模型箱净尺寸为700 mm×350 mm×650 mm（长×宽×高）。模型设计时取大坝标准剖面,按平面问题进行试验。模拟范围：竖向从建基面2176 m高程到坝顶2304.5 m高程,即整个坝高128.50 m,顺河向取坝轴线上游160 m、下游155 m,共315 m宽。标准剖面的上游坝坡为1∶2.75～1∶2.25,下游坝坡为1∶1.8～1∶1.6；而模型模拟范围内上游坝坡为1∶2.75～1∶2.5,下游坝坡为1∶1.8～1∶1.7,不考虑马道和上坝公路,试验将上、下游坝坡简化成同一坡度,上游坝坡取1∶2.75,下游坝坡取1∶2.6。模型几何比尺ηl =1∶450,离心加速度比尺ηg =40。模型布置见图6。
本次试验输入波形为实际地震波,地震峰值加速度有2个：设计地震波516.5g（50 a超越概率2%场地波）,对应无抗震措施模型M1和有抗震措施模型M3；校核地震波643.3g（100年超越概率2%场地波）,对应无抗震措施模型M2和有抗震措施模型M4。模型试验所需筑坝材料取自现场,根据以往粒径效应研究成果,模型土料的限制粒径应小于土作用构件最小边长的1/15～1/30,原型砂砾石料的最大粒径为300 mm,离心振动台模型砂砾石料的限制粒径取为20 mm,存在一个对超粒径颗粒进行处理的问题。依据《水电水利工程粗粒土试验规程》DL/T5356—2006,把原级配缩制成试验级配最常用的有相似级配法和等量替代法。根据现场检测级配曲线,用相似级配法与等量替代法相结合的方法确定砂砾石料的试验级配,如图7所示,其不均匀系数和曲率系数与现场料相近。沥青砼心墙采用有机玻璃板模拟,心墙厚度按抗弯刚度相似条件确定；模型护坡采用1～2 mm厚的砼,由于厚度较薄,不会影响坝体的变形或破坏；有抗震措施是指在筑坝材料区域铺设土工格栅,土工格栅采用纱窗模拟,抗拉强度要求满足相似条件,为了模拟实际工程中筑坝材料嵌入格栅的情形,模型格栅在铺设前,需要在上下面各刷一层胶,使其能够粘合少量的堆石料,从而可以模拟筑坝材料嵌入格栅的效果。
为了量测坝体地震加速度反应,在坝轴线上沿高程2202.4,2224.8,2246.9,2269.8,2290.1 m共布置5个加速度计,同时在坝顶布置1个激光位移传感器,测试坝体地震震陷。
3.2 结果分析
以下给出的所有试验结果均已根据模型相似律换算至原型。
（1）无抗震措施模型试验
a）坝体加速度反应
图8和图9分别给出了模型M1和M2对应的建基面原型场地波和试验实际输入波时程线,前者目标波峰值分别为504.514g和-510.21g,输入波峰值分别是501.937g和-523.127g；后者目标波峰值分别为645.023g和-644.977g,输入波峰值分别是646.388g和-649.306g,从中可以看出,波形的频响特性较为吻合,试验输入波形较好地再现了目标场地波。
![]()
图 8 模型M1目标地震波和实际输入波时程曲线Figure 8. Time histories of target seismic waves and actual input waves for model M1![]()
图 9 模型M2目标地震波和实际输入波时程曲线Figure 9. Time histories of target seismic waves and actual input waves for model M2图10给出了无抗震措施情况下坝体在2种场地波条件下的最大加速度和放大系数随坝高的变化。
![]()
图 10 坝体坝轴线最大加速度和放大系数沿坝高分布Figure 10. Distribution of maximum acceleration and amplification coefficient of acceleration along dam height由此可知：①在基岩输入地震加速度的作用下,坝体地震加速度反应随着高程的增加而相应增大,呈现出明显的放大效应；②坝体加速度反应随坝高的变化可以按约2/3坝高为界,大致分成两个线性变化段,上部的加速度放大效应强于下部；③随着基岩输入地震加速度的增加,坝体加速度放大系数总体上略有减小；④无抗震措施情况下坝体坝顶地震加速度放大系数约为2.7～3.0。
b）坝体破坏模式
图11和图12分别表示无抗震措施情况坝体在设计地震波作用下（M1）和校核地震波作用下（M2）的破坏情况。
从上述试验照片可以看出：①无抗震措施情况坝体在设计地震波条件下,下游未发现明显破坏,上游出现了一定程度的坝体破坏；在蓄水位处可见明显的护坡开裂；②无抗震措施情况坝体在校核地震波条件下,下游未发现明显破坏,上游出现了较大程度的坝体破坏；在蓄水位处和坝顶处均有明显的护坡开裂。
（2）有抗震措施模型试验
图13,14分别给出了2种场地波条件下有抗震措施和无抗震措施坝体的最大加速度和放大系数随坝高的变化。
![]()
图 13 设计地震波作用下坝轴线加速度及放大系数分布Figure 13. Distribution of acceleration and amplification coefficient of acceleration under design seismic waves![]()
图 14 校核地震波作用下坝轴线加速度及放大系数分布Figure 14. Distribution of acceleration and amplification coefficient of acceleration under check seismic waves由上述图中可以看出：①在基岩输入地震加速度的作用下,有抗震措施坝体地震加速度反应随着高程的增加而相应增大,呈现出明显的放大效应；②有抗震措施坝体加速度反应随坝高的变化可以按约2/3坝高为界,大致分成两个线性变化段,上部的加速度放大效应强于下部；③在抗震措施加固范围以下高程,有抗震措施和无抗震措施坝体的加速度放大效应大体一致；在抗震措施加固范围内,有抗震措施坝体的加速度放大效应略强于无抗震措施坝体；④有抗震措施坝体坝顶地震加速度放大系数约为3.0～3.4。
从试验结果分析可知,离心机振动台较好地再现了在地震作用下坝体内部加速度动态响应规律及坝体失稳破坏模式,为工程设计部门提供了参考依据,进一步验证了NS-2型离心机振动台设计性能指标完全满足工程科研需要。
4. 结论
（1）介绍南京水利科学研究院40 g·t离心机振动台设计关键技术和性能指标,主要包括机械系统、液压系统、伺服控制系统,可为土工离心机振动台建设和规划提供参考。
（2）NHRI-40 g·t离心机振动台的技术指标科学合理,设备的研制思路正确,波形控制系统可实现对目标波形的跟踪、调节和再现,目标波形和再现波形拟合度高；基于Labview编制的数据采集系统包括16个应变采集通道和24个电压采集通道,对试验数据进行实时测量、显示和存储。
（3）通过该离心机振动台测得国内某大坝的地震加速度响应分布表明,振动台的主要技术指标和性能达到了预期目的,完全能满足工程试验研究的需求。
-
![]()
图 1 安徽某隧道工程泥水盾构排渣情况
Figure 1. Slag discharge situation of slurry shield tunneling in a tunnel project in Anhui Province
![]()
图 10 砾石床层对管道内泥浆的影响示意图
Figure 10. Schematic diagram of impact of gravel bed on slurry inside pipeline
![]()
图 11 第2模拟组泥浆速度和压力横断面分布图
Figure 11. Cross section distribution of slurry velocity and pressure in second simulation group
![]()
图 14 不同尺寸砾石床层的速度分布图
Figure 14. Velocity distribution of gravel beds with different sizes
![]()
图 15 砾石尺寸和床层流态的关系图
Figure 15. Diagram of relationship between gravel size and bed flow state
![]()
图 16 砾石尺寸对管道内泥浆增压和增速的影响
Figure 16. Effects of gravel size on increase of slurry pressure and flow velocity in pipelines
![]()
图 17 不同形状砾石床层的速度分布图
Figure 17. Velocity distribution of gravel beds with different shapes
![]()
图 19 砾石形状对管道内泥浆增压和增速的影响
Figure 19. Effects of gravel shape on increase of slurry pressure and flow velocity in pipelines
![]()
图 20 不同相对质量密度泥浆输送砾石床层的速度分布图
Figure 20. Velocity distribution of gravel bed transported by slurry with different densities
![]()
图 21 泥浆相对质量密度和床层流态的关系图
Figure 21. Relationship between slurry density and bed flow state
![]()
图 22 泥浆相对质量密度对管道内泥浆增压和增速的影响
Figure 22. Effects of slurry specific gravity on increase of slurry pressure and flow velocity in pipelines
![]()
图 23 不同流速泥浆输送砾石床层的速度分布图
Figure 23. Velocity distribution of gravel bed transported by slurry with different flow rates
![]()
图 24 泥浆流速和床层流态的关系图
Figure 24. Relationship between slurry flow rate and bed flow state
![]()
图 25 泥浆流速对管道内泥浆增压和增速的影响
Figure 25. Effects of slurry flow rate on increase of slurry pressure and flow velocity in pipelines
![]()
图 26 1~16模拟组沉积风险分区图
Figure 26. Zoning map of sedimentary risk for simulation groups 1~16
表 1 岩石模拟参数表
Table 1 Simulation parameters of rocks
密度/
(g·cm-3)弹性模量/
GPa泊松比 抗压强度/
MPa2.70 17.5 0.16 42.5~116.2 
下载: 导出CSV
表 2 仿真接触参数表
Table 2 Contact parameters for simulation
相互作用 恢复系数 静摩擦系数 滚动摩擦系数 颗粒-颗粒 0.4 0.32 0.012 颗粒-管道 0.5 0.28 0.01
下载: 导出CSV
表 3 仿真模拟方案表
Table 3 Simulation plans
模拟
组号砾石
形状砾石尺寸/
cm砾石
数量泥浆流速/
m/s泥浆相对质量密度 1 纺锤形 2.5 2500 2.0 1.15 2 纺锤形 3.0 2500 2.0 1.15 3 纺锤形 3.5 2500 2.0 1.15 4 纺锤形 4.0 2500 2.0 1.15 5 圆形 — 2500 2.0 1.15 6 纺锤形 — 2500 2.0 1.15 7 三角形 — 2500 2.0 1.15 8 条形 — 2500 2.0 1.15 9 三角形 3.0 2500 1.5 1.15 10 三角形 3.0 2500 2.0 1.15 11 三角形 3.0 2500 2.5 1.15 12 三角形 3.0 2500 3.0 1.15 13 纺锤形 2.5 2500 2.0 1.05 14 纺锤形 2.5 2500 2.0 1.15 15 纺锤形 2.5 2500 2.0 1.25 16 纺锤形 2.5 2500 2.0 1.35
下载: 导出CSV
-
[1] 金大龙, 袁大军, 郑浩田, 等. 高水压条件下泥水盾构开挖面稳定离心模型试验研究[J]. 岩土工程学报, 2019, 41(9): 1653-1660. JIN Dalong, YUAN Dajun, ZHENG Haotian, et al. Centrifugal model tests on face stability of slurry shield tunnels under high water pressures[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(9): 1653-1660. (in Chinese)
[2] 徐涛, 史庆锋, 章定文, 等. 泥水盾构开挖面泥膜渗透特性与压力传递机制[J]. 岩土工程学报, 2023, 45(9): 1878-1887. XU Tao, SHI Qingfeng, ZHANG Dingwen, et al. Permeability characteristics of filter cake and pressure transfer on face during slurry shield tunnelling[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(9): 1878-1887. (in Chinese)
[3] 李承辉, 贺少辉, 刘夏冰. 粗粒径砂卵石地层中泥水平衡盾构下穿黄河掘进参数控制研究[J]. 土木工程学报, 2017, 50(增刊2): 147-152. LI Chenghui, HE Shaohui, LIU Xiabing. Study on main parameters control of tunneling through the Yellow River by a slurry balance shield in sandy gravel stratum with some large-size grains[J]. China Civil Engineering Journal, 2017, 50(S2): 147-152. (in Chinese)
[4] 孔玉清. 泥水盾构环流系统及排泥管携碴能力分析与应用[J]. 现代隧道技术, 2018, 55(3): 205-213. KONG Yuqing. Analysis of the circulation system of a slurry shield and the muck carrying ability of a dredging pipe[J]. Modern Tunnelling Technology, 2018, 55(3): 205-213. (in Chinese)
[5] 黄子木, 李兴高, 杨益. 泥水盾构泥浆输送渣土冲蚀磨损力学性能的研究[J]. 土木工程学报, 2020, 53(增刊1): 25-30. HUANG Zimu, LI Xinggao, YANG Yi. Mechanical properties of pipeline erosion duringslurry transportation of slurry shield[J]. China Civil Engineering Journal, 2020, 53(S1): 25-30. (in Chinese)
[6] MESSA G V, MALAVASI S. Improvements in the numerical prediction of fully-suspended slurry flow in horizontal pipes[J]. Powder Technology, 2015, 270: 358-367.
[7] WASP B E J, KENNY J P, GANDHI R L. Solid-Liquid Flow Slurry Pipeline Transportation[M]. Houston: Gulf Pub Co, 1979.
[8] JAYASUNDARA C T, YANG R Y, GUO B Y, et al. Effect of slurry properties on particle motion in IsaMills[J]. Minerals Engineering, 2009, 22(11): 886-892.
[9] 夏毅敏, 王洋, 吴遁, 等. 泥水盾构环流系统管道输送特性[J]. 中南大学学报(自然科学版), 2017, 48(11): 2889-2896. doi: 10.11817/j.issn.1672-7207.2017.11.037 XIA Yimin, WANG Yang, WU Dun, et al. Transport characteristics of shield slurry system[J]. Journal of Central South University (Science and Technology), 2017, 48(11): 2889-2896. (in Chinese) doi: 10.11817/j.issn.1672-7207.2017.11.037
[10] 吴遁. 卵石地层中泥水盾构管道出渣影响因素研究[J]. 人民黄河, 2020, 42(5): 157-161. WU Dun. Study on the influencing factors of the ballast carrying capacity of the\r slurry shield pipe under the sandy pebble formation[J]. Yellow River, 2020, 42(5): 157-161. (in Chinese)
[11] MIN F L, SONG H B, ZHANG N. Experimental study on fluid properties of slurry and its influence on slurry infiltration in sand stratum[J]. Applied Clay Science, 2018, 161: 64-69. doi: 10.1016/j.clay.2018.03.028
[12] 朱红钧, 林元华, 谢龙汉. FLUENT流体分析及仿真实用教程[M]. 北京: 人民邮电出版社, 2010. ZHU Hongjun, LIN Yuanhua, XIE Longhan. Practical Course of Fluent Fluid Analysis and Simulation[M]. Beijing: Posts & Telecom Press, 2010. (in Chinese)
[13] YANG D, XIA Y M, WU D, et al. Numerical investigation of pipeline transport characteristics of slurry shield under gravel stratum[J]. Tunnelling and Underground Space Technology, 2018, 71: 223-230. doi: 10.1016/j.tust.2017.08.011
[14] ADNAN M, SUN J, AHMAD N, et al. Comparative CFD modeling of a bubbling bed using a Eulerian–Eulerian two-fluid model (TFM) and a Eulerian-Lagrangian dense discrete phase model (DDPM)[J]. Powder Technology, 2021, 383: 418-442.
[15] BADR H M, HABIB M A, BEN-MANSOUR R, et al. Numerical investigation of erosion threshold velocity in a pipe with sudden contraction[J]. Computers and Fluids, 2005, 34(6): 721-742.
[16] 张新卓. 异形粗颗粒管道输送流态及动力学特性研究[D]. 武汉: 武汉理工大学, 2019. ZHANG Xinzhuo. Study on Flow Pattern and Dynamic Characteristics of Special-Shaped Coarse Particles Transported in Pipeline[D]. Wuhan: Wuhan University of Technology, 2019. (in Chinese)
[17] 王洋. 盾构环流系统石碴运动特性及管道使用寿命预测研究[D]. 长沙: 中南大学, 2023. WANG Yang. Study on Movement Characteristics of Ballast in Shield Circulating System and Prediction of Pipeline Service Life[D]. Changsha: Central South University, 2023. (in Chinese)
[18] UZI A, LEVY A. Flow characteristics of coarse particles in horizontal hydraulic conveying[J]. Powder Technology, 2018, 326: 302-321.
-
期刊类型引用(1)
1. 王永志，贾仕林，韩俊伟，王体强，张雪东，汤兆光，周燕国，张宇亭. 迭代算法与柔性结构对离心机振动台控制性能影响. 岩土工程学报. 2024(S1): 22-26 . 
本站查看
其他类型引用(4)
-
其他相关附件
计量
- 文章访问数: 223
- HTML全文浏览量: 26
- PDF下载量: 21
- 被引次数: 5
