Numerical Simulation of Negative Pressure Transport of Extra Large Particles of Suction Earthquake Rescue Machine
-
摘要: 针对某抽吸式地震救援机气力输送系统工作时存在的易堵塞和能耗高等问题,选用不同输送管径开展特大颗粒物料气固两相流数值模拟,研究风机以额定功率工作时特大颗粒在不同管道中的运动特性和管道内压损失。利用计算流体力学和离散单元法进行耦合求解,在传统耦合接口的基础上更改体积分数算法实现颗粒粒径大于流体域网格单元的气力输送数值仿真,利用软件Fluent进行连续相气体的流动分析和曳力计算,利用软件EDEM进行离散相颗粒的受力分析和运动轨迹计算。研究结果表明,管径为280 mm的管道输送性能最佳,管中颗粒运动速度最大为7.1 m/s,0~25 s内的运输量达66.4 kg,管道内压损失在5.5 s时达到稳定,为3.04 kPa,与试验结果基本吻合,可为抽吸式地震救援机气力输送系统设计和优化提供参考。Abstract: Aiming at the problems of easy blockage and high energy consumption in the pneumatic conveying system of a pneumatic suction seismic rescue machine, the gas-solid two-phase flow of large particles in different pipes with different diameters was numerically simulated to study the movement characteristics and pressure loss of large particles in different pipes when the fan works at rated power. Computational fluid mechanics and discrete element method were used to solve the coupling, and the volume fraction algorithm was changed on the basis of the traditional coupling interface to realize the numerical simulation of pneumatic conveying with the particle size larger than the fluid domain grid element. The flow analysis and drag calculation of continuous phase gas were carried out by using the software Fluent. The force analysis and motion trajectory calculation of discrete phase particles are carried out by using Edem software. The results show that the pipe with a diameter of 280 mm has the best transport performance, the particle movement velocity in the pipe is up to 7.1 m/s, the transport volume within 0~25 s is up to 66.4 kg, and the pressure loss in the pipe reaches a stable level of 3.04 kPa at 5.5 s, which is basically consistent with the test results. It can provide reference for the design and optimization of pneumatic conveying system of suction earthquake rescue aircraft.
-
引言
跨断层形变监测是一种观测物理意义明确、测量精度高、便于组织实施以及数据产品稳定可靠的地震前兆监测方法之一(薄万举等, 1994, 1998;李杰等,2010;张晶等,2011;张希等,2014;马伶俐等,2015;李腊月等,2016;孙启凯等,2017)。通过多期次、短距离的跨断层水准和基线测量,能够准确记录断层活动的位移量,监测活动断裂的运动特征,更好地反映断裂带及其所处区域的形变和运移特征(牛安福等,2013)。针对川滇地块的地壳形变、运动学特征及其形成机制开展的研究工作产出了大量成果,但是,目前关于这些问题也存在不同认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015;白玉柱等,2016)。本文利用川滇地块南部附近的跨断层观测资料,对区域断裂带的现今形变变化进行分析,研究川滇地块南部的整体运动学特征,综合提取断层活动的异常信息,为区域地震危险性分析提供一定的参考。
1. 区域地质概况
川滇地块位于青藏高原的东缘,南北地震带的南段,大地构造上处于扬子板块、松潘-甘孜地块、兰坪-思茅地块以及冈底斯-腾冲地块等构造单元的拼合汇聚地带,在地质历史上经历了多次不同时期不同规模的构造运动与岩浆活动,构造属性复杂。该区域地震活动强烈而频繁,地震与构造运动之间存在复杂的因果关系,是研究地震活动与断层运动关系的热点地区(图 1)。
自新生代以来,印度板块与亚欧板块的相互碰撞拼贴,造成地壳强烈缩短增厚,青藏高原得以快速隆升(Tapponnier等,2001;许志琴等,2011)。伴随着两大板块的俯冲拼合,高原也发生着同构造挤压伸展,物质和能量向北西、北东、东、南东等多个方向转移释放。这种持续的俯冲碰撞使得两大板块拼合区域的喜马拉雅造山带、雅鲁藏布江缝合带以及高原内部地块与周缘区域构造活动强烈,并形成一系列大断裂带,如北西走向的鲜水河断裂、红河断裂及近南北走向的小江断裂、安宁河-则木河断裂带等。很多学者对此区域开展了研究,并取得了大量成果(Wang等,1998;Burchfiel等,2003;陈长云等,2016)。其中,处于川滇地块南部的活动断裂,是本文关注的重点。
2. 跨断层监测概况
川滇地块南部断裂所在的云南地区的跨断层监测场地从20世纪80年代开始观测,至今30余年,主要为流动观测场地,共有11处,分别是楚雄、峨山、建水、剑川、丽江、石屏、通海、下关、永胜、宜良和羊街,主要分布于金沙江-红河断裂带、小江断裂带、南华-楚雄-建水断裂带及其次级断裂(剑川断裂、洱海断裂、通海-峨山断裂、建水-石屏断裂、丽江断裂)上(表 1),多为综合观测场地。每个场地包括2条短水准和2条短基线测线,分别各有1条直交和斜交测线,通过这种组合方式来共同监测场地所跨断层的三维运动变化。从场地观测历史来看,该区域的跨断层观测基本是每月1期,个别年份每2个月观测1期,数据连续且完整。
表 1 研究区跨断层监测场地概况Table 1. Basic information of the cross-fault monitoring sites in the area场地 经度/°E 纬度/°N 所跨断层 断层产状 断层性质 楚雄 101.52 25.03 楚雄断裂 NE∠58° 右旋挤压 峨山 102.50 24.12 通海-峨山断裂 S∠60° 右旋挤压 建水 102.81 23.63 石屏-建水断裂 S∠70° 左旋挤压 剑川 99.93 26.63 剑川断裂 W∠75° 左旋拉张 丽江 100.25 26.88 丽江断裂 SW∠70° 左旋挤压 石屏 102.45 23.77 石屏-建水断裂 NE∠70° 右旋拉张 通海 102.72 24.00 通海-峨山断裂 NE∠68° 右旋挤压 下关 100.30 25.62 红河断裂 NE∠60° 右旋拉张 永胜 100.72 26.70 程海断裂 W∠35° 左旋挤压 宜良 103.03 25.02 马街-南羊街断裂 SE∠68° 左旋挤压 羊街 103.10 25.55 小江断裂西支 SE∠73° 左旋挤压 3. 数据计算及分析
将收集到的研究区30余年的跨断层观测数据进行预整理,对数据格式做统一标准化处理。首先,消除跨断层观测资料中明确因干扰引起的突跳变化,对观测资料进行连续性检验,剔除重复的数据。然后,对于日常观测中因各种原因造成缺测的数据采用线性插值法修补,得到周期统一的、可靠的形变数据。对数据的处理分析采用以下方法:
云南区域内的跨断层场地同时布设有两条同桩的基线和水准测线,可采用下面3个公式来计算断层活动的三维分量:
$$\Delta S=\frac{\Delta {{L}_{1}}\text{sin}{{\alpha }_{2}}-\Delta {{L}_{2}}\sin {{\alpha }_{1}}}{\sin ({{\alpha }_{2}}-{{\alpha }_{1}})}$$ (1) $$\Delta R=\frac{\Delta {{L}_{1}}-\Delta S\cos {{\alpha }_{1}}}{\text{sin}{{\alpha }_{1}}}$$ (2) $$\Delta H=\frac{1}{2}(\Delta {{H}_{1}}+\Delta {{H}_{2}})$$ (3) 其中,$\Delta H$为水准垂直变化量(以上盘相对于下盘的下降为正),$\Delta L$为基线伸缩量(以伸长为正),$\Delta S$为断层活动水平走滑量(两盘做顺扭运动为正),$\Delta R$为断层活动水平张压量,${{\alpha }_{1}}$和${{\alpha }_{2}}$分别代表由断层线方向逆时针转动至与两条测线方向重合或者平行时的角度。当$\Delta S$为正值时,表示右旋活动,$\Delta R$为正值时,表示张性活动,若均为负值则表示断层活动性质相反(薄万举等,1998)。
利用MATLAB对计算得出的数据绘制曲线图,单幅曲线图上方均标注有其代表的场地名及运动方式(图 2)。通过曲线图可以发现跨断层场地所观测的不同断裂表现出了跟地质背景相同或者不同的运动学特征,其中丽江、石屏和下关场地所跨断层的现今活动特征跟地质背景一致;楚雄、剑川和永胜场地所跨断层在走滑和拉张特征上都表现出跟背景特征不一致的现象;峨山、建水和通海场地所观测的断层在走滑性质上和地质背景相符,而在张压性质上相反。根据现今运动特征,可将场地分为3种类型:
(1)运动特征与地质背景一致。丽江场地所跨的丽江次级断裂,在地质背景上是一处表现为左旋走滑特征的活动断裂,兼具挤压特点,现今计算的三维活动量,依然表现出相同特征,其中断层的水平张压量在2012年后出现大幅异常变化。石屏场地所跨的石屏-建水断裂,在地质背景上是一处右旋拉张特征的活动断裂,现今右旋运动速率在2013年后虽然有所放缓,但是仍然与地质背景一致。下关场地所处的是一条表现为右旋走滑兼有拉张特征的深断裂,2003年水平走滑量曾出现右旋剪切转折,现今曲线重回向上趋势,继续保持右旋走滑特征。羊街和宜良场地所处的小江断裂带西支断裂,现今数据仍显示出以继承性的左旋走滑兼具挤压特征为主。
(2)运动特征与地质背景不一致。楚雄场地所跨断层为具有右旋挤压背景特征的楚雄断裂,在2008年后现今运动特征以强烈左旋拉张为主,但自2014年左旋和拉张的运动速率都出现大幅放缓的异常现象。剑川场地所处断裂为红河断裂带北段的剑川断裂,是一处左旋走滑兼具拉张特征的活动断裂,现今三维运动特征表现为右旋走滑兼具挤压特征,其中的垂直位移量自2015年至今出现了大幅抬升的异常变化;水平走滑量自2011年后也出现大幅抬升,右旋运动加速。永胜场地所监测的断裂是具有左旋挤压特征的程海深断裂,从计算所得的三维活动量曲线可以看出,该断裂现在运动学特征为右旋拉张,其垂直位移量曲线从2014年开始大幅异常抬升,水平走滑量自1994年开始稳步抬升。
(3)与地质背景相比,走滑特征一致,张压特征相反。峨山场地和通海场地共同监测曲江断裂的不同部位,该断裂在地质背景上具有右旋挤压特征,两者现今走滑特征与背景一致,但是张压量曲线表现出相反的拉张运动的异常特点。其中,峨山场地走滑分量在1998年后大幅抬升,之后趋于稳定的右旋运动状态;通海场地的垂直位移量在1992年后出现转折,断层表现为张性活动,走滑分量在2003年大幅下降,断层剪切活动速率减小,之后趋于稳定,依然显示右旋运动。建水场地所跨的石屏-建水断裂的建水段,断层性质为左旋挤压,从计算得出的水平走滑量参数来看,现今的走滑运动特征与地质背景一致,垂直位移量在1999年后发生趋势改变,断层活动由压性转变为相反的拉张特征。
4. 讨论与结论
在印度洋板块和亚欧板块的互相作用下,巨大的青藏高原造就了地球的“第三极”,这种作用力至今仍在持续。板块的碰撞拼合产生了强烈的构造活动,许多大规模断裂带应运而生(Xu等,2016),同时,也孕育了一系列强震。对此一些学者做过许多研究,并取得大量有价值的科研成果(Tapponnier等,2001;许志琴等,2011;陈长云等,2016)。尽管目前对青藏高原的陆内侧向运动形式、时间和逃逸量存在着争议,但是对位于其东缘的川滇地块发生的晚新生代南东向的逃逸和围绕东喜马拉雅构造结顺时针的旋转有着一致的认识(乔学军等,2004;武艳强等,2011;洪敏等,2014;吴中海等,2015)。
通过对位于川滇地块南部主要断裂上的跨断层场地三维形变量的分析,可以看出部分断裂表现出继承性特征,如现今下关和石屏场地的三维活动量表现出跟以往相同的持续平稳的活动速率。而另一部分断裂表现出了异常于背景活动特征的现象,如楚雄、永胜和剑川场地在走滑分量上都表现出了反向活动特点。研究区具有以水平运动为主、垂直运动为辅的背景特征,现今看来,依然如此。通过计算分析得到的宜良、羊街、楚雄以及丽江场地的形变结果来看,羊街、宜良场地所控制的小江断裂带分支为川滇地块东部边界断裂,现呈现左旋走滑兼具压性的特征;楚雄场地监控的楚雄-南华断裂位于川滇地块的南边部,表现出左旋拉张的运动特征;丽江场地所处的丽江-小金河断裂作为川滇地块南北部的分界断裂,现今以左旋挤压运动为主。可以看出,川滇地块南部区域整体上表现出以水平顺时针旋转为主的运动方式。根据GPS观测资料得到的该区域现今水平运动速度场也显示了相同的运动特点(乔学军等,2004;武艳强等,2011;洪敏等,2014;王伶俐等,2015)。造成这种运动方式的深层次动力学机制,应该源于印度洋板块俯冲至亚欧板块下,碰撞拼合后巨大高耸的青藏高原崛起,在超级重力作用下,物质开始向周缘,尤其是向东缘流动,在遇到横亘的北东向龙门山造山带时受阻,流动速率下降,被消化、吸收,再次转移(张家声等,2003;崔效锋等,2006)。其中,相当一部分能量和物质向南东运移,规模相对较小的川滇活动地块随之发生南东向旋转,同时起着调节周围各大板块和地块之间运动平衡的“交通环岛”的作用。另外,印度洋板块至今持续强烈地北向运动,跟扬子板块形成了大尺度上的右旋运动,也推动着川滇地块南东向旋转运移。
根据观测场地所跨断层的三维运动特征及计算所得曲线可以看出,在左旋的小江断裂带、右旋的红河断裂带和右旋的南华-楚雄-建水断裂带共同汇聚的峨山、通海和建水场地区域,原先的地壳形变表现出压性特征,现今由于断裂性质的异常变化,该区域场地集体表现出拉张的运动特点。滇西北区域的永胜、丽江和剑川场地所跨断层也表现出了一定的异常运动特征,剑川断裂和程海断裂呈水平反向运动,除此之外,包括丽江场地所跨的丽江断裂在内,3条断裂的垂向运动特征也存在大幅异常变化。云南地区中强震发生前的一定时段或震时,区域内跨断层场地观测到的水准和基线数据多数会出现较明显的变化,且部分场地形变量会发生速率的急变或转折,如文山5.3级、孟连西7.3级、姚安6.5级、丽江7.0级地震震前和震时,部分场地就表现出了加速继承性运动和转折性运动的特征(张兴华等,1996;付虹等,1997;李忠华等,1998;施顺英等,2005;李瑞莎等,2009;刘强等,2010;王永安等,2011)。同时,云南地区2015年至今存在长时间缺少4级以上中强地震的背景,不符合区域以往的地震活动特征,现今应力积累较高。自2016年云龙5.0级地震后至今,区域中强震活动已平静300天左右;滇西北地区自2001年永胜6.0级地震至今,已14.7年未发生6级以上地震;滇中南的石屏、建水等区域从2001年江川5.1级地震起算,5级地震平静期已长达15.5年。基于此,该区域在中短期尺度可能存在一定中强震的危险。
综上所述,得到以下几点结论:
(1)川滇地块南部跨断层场地的形变特征表现出继承和异常于地质背景特征的现象。① 丽江、石屏、下关、宜良和羊街场地所跨断层的现今活动特征与地质背景一致。丽江断裂和小江断裂表现出持续地左旋挤压运动特征;石屏-建水断裂石屏段表现出右旋拉张特征;红河断裂表现出跟地质背景相同的右旋拉张特征。② 楚雄、剑川和永胜场地所跨断层在走滑和张压特征上都出现与背景特征相背的现象。楚雄断裂、剑川断裂和程海断裂均表现出相反的运动特征。③ 峨山、建水和通海场地所监测的断层现今活动状态与地质背景相比,在走滑性质上相符,张压性质上相反。通海-峨山断裂现今运动特征,表现出与地质背景相同的右旋走滑和与地质背景相反的拉张特征。石屏-建水断裂建水段现今依然呈左旋走滑运动,但是张压性质发生反转,表现出张性运动的活动特点。
(2)川滇地块南部整体上表现出以水平顺时针旋转为主的运动方式。
(3)根据跨断层观测资料所得出的断层三维运动特征及曲线可知,滇西北及滇中的峨山、通海、石屏和建水等区域地壳形变特征存在异常,应关注该区域中短期尺度内可能存在的中强震危险。
致谢: 对编辑和审稿老师的辛苦审阅和宝贵建议,表示衷心感谢。 -
表 1 数值仿真过程参数
Table 1. Corresponding parameters of soil layers
种类 参数 数值 气相性质 密度/(kg·m−3)
粘度/(kg·m−1·s−1)
质量入口/(kg·s−1)
压力出口/Pa1.169
1.55×10−5
4.6
0颗粒性质 密度/(kg·m−3)
泊松比
剪切模量/Pa
粒径/mm
生成速率/(kg·s−1)1 800
0.3
1.0×107
101.6
3.3壁面性质 密度/(kg·m−3)
泊松比
剪切模量/Pa1 210
0.4
7.0×107颗粒-颗粒 碰撞恢复系数
静摩擦系数
滚动摩擦系数0.4
0.5
0.04颗粒-管壁 碰撞恢复系数
静摩擦系数
滚动摩擦系数0.5
0.2
0.03表 2 不同直径管道内气体流速与运输量
Table 2. Gas flow velocity and transport volume in different diameters
管径/mm 气体质量流量/(kg·s−1) 空管时管内气体平均流速/(m·s−1) 0~25 s内运输颗粒总质量/kg 230 4.6 94.6 73.9 250 4.6 80.1 72.2 280 4.6 63.9 66.4 300 4.6 55.7 0 表 3 不同颗粒生成速率下的运输量
Table 3. The amount of transport at different generation rates
颗粒生成速率/(kg·s−1) 0~25 s内运输颗粒总质量/kg 3.3 0 3.0 41.5 2.5 34.0 2.0 18.3 表 4 试验测量与仿真结果对比
Table 4. Experimental results and simulation results
管径/mm 试验测量管道
内压损失/Pa仿真得到管道
内压损失/Pa仿真结果
相对误差/%试验测量运输
速率/(t·h−1)试验测量0~25 s内
运输量/kg仿真得到0~25 s内
运输量/kg仿真结果
相对误差/%230 8253 7543 −8.6 9.7 67.6 73.9 9.3 250 6028 5558 −7.8 9.3 64.2 72.2 12.5 280 3352 3004 −10.4 8.3 57.3 66.4 13.7 -
崔益华, 2021. 水平管道负压气力输送CFD-DEM数值模拟. 港口装卸, (1): 63—67Cui Y. H. , 2021. CFD-DEM numerical simulation of horizontal pipeline negative pressure pneumatic conveying. Port Operation, (1): 63—67. (in Chinese) 杜俊, 2015. 基于CFD-DEM方法的稀相气力输送数值模拟研究. 武汉: 武汉大学.Du J., 2015. Simulation of dilute pneumatic conveying by CFD-DEM. Wuhan: Wuhan University. (in Chinese) 封凯, 聂伟, 陈凤官等, 2021. 气力输送系统中弯管的易磨损位置及其机理分析. 化学反应工程与工艺, 37(2): 106—112Feng K. , Nie W. , Chen F. G. , et al. , 2021. Positioning the erosion wear of bend pipes in the pneumatic conveying system. Chemical Reaction Engineering and Technology, 37(2): 106—112. (in Chinese) 胡国明, 2010. 颗粒系统的离散元素法分析仿真. 武汉: 武汉理工大学出版社.Hu G. M., 2010. Analysis and simulation of granular system by Discrate element method using EDEM. Wuhan: Wuhan University of Technology Press. (in Chinese) 胡昱, 黄小华, 陶启友等, 2019. 基于CFD-EDM的自动投饵饲料颗粒气力输送数值模拟. 南方水产科学, 15(3): 113—119Hu Y. , Huang X. H. , Tao Q. Y. , et al. , 2019. Simulation of automatic pneumatic conveying process of automatic feeding system based on computational fluid dynamics-discrete element method. South China Fisheries Science, 15(3): 113—119. (in Chinese) 柳波, 吴中鼎, 尹高冲等, 2018. 供料器内颗粒碰撞条件下的两相流压力特性仿真分析. 现代制造工程, (3): 93—98Liu B. , Wu Z. D. , Yin G. C. , 2018. Simulation analysis of the pressure characteristics of two-phase flow in feeding device based on particle collision. Modern Manufacturing Engineering, (3): 93—98. (in Chinese) 刘德天, 傅旭东, 王光谦, 2017. CFD-DEM耦合计算中的体积分数算法. 清华大学学报(自然科学版), 57(7): 720—727Liu D. T. , Fu X. D. , Wang G. Q. , 2017. Volume fraction allocation using characteristic points for coupled CFD-DEM calculations. Journal of Tsinghua University (Science and Technology), 57(7): 720—727. (in Chinese) 庞林祥, 贾志勇, 2018. 我国地震救援装备现状与发展趋势. 水利水电快报, 39(6): 27—30. 王国强, 郝万军, 王继新, 2010. 离散单元法及其在EDEM上的实践. 西安: 西北工业大学出版社.Wang G. Q. , Hao W. J. , Wang J. X. , 2010. Discrete element method and its practice on EDEM. Xi’an: Northwestern Polytechnical University Press. (in Chinese) 周甲伟, 上官林建, 许兰贵等, 2022. 粗重颗粒轻介共流气力输送机理及特性. 机械工程学报, 58(14): 308—319Zhou J. W. , Shangguan L. J. , Xu. L. G. , 2022. Mechanism and characteristics of light medium mixed flow pneumatic conveying for coarse particle. Journal of Mechanical Engineering, 58(14): 308—319. (in Chinese) Bharadwaj R. , Ketterhagen W. R. , Hancock B. C. , 2010. Discrete element simulation study of a Freeman powder rheometer. Chemical Engineering Science, 65(21): 5747—5756. doi: 10.1016/j.ces.2010.04.002 Laín S. , Sommerfeld M. , 2019. Numerical prediction of particle erosion of pipe bends. Advanced Powder Technology, 30(2): 366—383. doi: 10.1016/j.apt.2018.11.014 Yang D. L. , Wang Y. X. , Li J. P. , et al. , 2020. Simulation study on interaction coefficient of DEM in non–spherical large size (5–30 mm) coal particles. Particuology, 53: 142—153. doi: 10.1016/j.partic.2020.03.001 -