Finite element method for analyzing environmental vibration without apparent sources and its application in Beijing High-Energy Photon Source
-
摘要: 无源环境微振动是指工程场地中没有明显振源时存在的微小振动,其可能导致高精密科学装置无法正常工作,如何准确分析无源条件下结构的环境微振动是此类工程亟需攻克的难题。系统介绍了一种无源环境条件下微振动的有限元分析方法。首先采用实测与理论方法得到地表任意点的振动时程,再利用传递矩阵法计算土层瑞利波弥散曲线及各频率振动沿深度衰减因子,再结合地表振动时程可计算出任意深度处的振动时程。将上述方法计算得到的每个节点的振动时程施加到有限元模型边界的相应节点上,使有限元模型的边界振动场与实际环境微振动场趋于一致。最终将该方法应用于国家“十三五”重大科技基础设施北京高能同步辐射光源工程(HEPS)地基基础的环境微振动分析,结果表明该方法能很好地反映无源振动的特性,为类似高精密科学装置的无源环境微振动分析提供了重要参考。Abstract: The environmental vibration without apparent sources refers to the vibration with very small amplitude in the field where there are no apparent excitation sources. It may significantly affect the normal operation of high-precision scientific apparatuses, and therefore appropriate methods are needed to analyze such environmental vibration precisely. A dynamic finite element method is introduced to analyze the environmental vibration without apparent sources in this study. Firstly, the time histories of vibration at the soil surface are obtained through the field measurements. Then, the transfer-matrix method is utilized to calculate the dispersion curve of the Rayleigh waves and the vibration attenuation factors with depth for each frequency. With the measured time histories on the ground and the attenuation factors, the time histories at different depths can be calculated. The calculated time history at each node is applied to the corresponding node on the finite element model boundary, which ensures the similarity of vibration conditions in the model and the reality. Finally, the proposed method is adopted to analyze the environmental vibration of the raft foundation of High-Energy Photon Source (HEPS). The results indicate that the proposed method can successfully capture the characteristics of vibration without apparent sources and provide useful reference for similar projects.
-
0. 引言
随着经济的不断发展和城市人口的增长,中国城市内河的污染问题日益凸显。据统计中国有80%以上的城市内河受到了不同程度污染[1]。河湖底泥的污染物主要分为有机质、氮磷营养盐、重金属三大类。
底泥污染物含量通常受底泥物理-化学性质的影响较大,大部分的有机质与黏土矿物结合在一起,并随细颗粒含量增加而增加,且有机质含量与底泥的比表面积线性相关[2];矿物表面的吸附过程对于有机质的保存起着重要作用,不同黏土矿物对有机质的吸附机理不同[3];底泥液限、塑限和塑性指数与黏粒中有机碳含量、黏粒含量与蒙脱石含量显著相关[4];不同矿物和不同有机质对底泥的物理性状影响显著不同[5]。
底泥中的有机质、营养盐和各种重金属影响底泥的物理性质和工程性质,进而影响底泥的处理技术与效果。因此,分析河湖底泥污染物与底泥物理–化学性质的相关性,对污染底泥的处置以及资源化利用具有重要的工程意义。
已有的研究成果表明了河湖底泥污染物与底泥物理–矿物成分密切相关,但是底泥物理–矿物成分–污染性状关联性的实例研究较少。本研究针对福州市晋安区河道的5处代表性污染底泥,进行了物性指标、矿物成分与污染物含量试验,同时搜集已有的国内外不同底泥污染物含量数据,分析了底泥物理–矿物成分与底泥污染物含量的关联性,并且探讨了底泥中不同污染物含量的相关关系,为河湖底泥的污染治理和处理技术选择提供科学依据。
1. 材料与方法
选取福州晋安区水系5个代表性点位,分别用A,B,C,D,E表示,如图 1所示。该水系有两条干流,分别是凤坂河和浦东河,浦东河有3条支流,分别是福兴河、新厝河、淌洋河,浦东河干流的最下游处为一个公园内的人工湖。B位于凤坂河干流的中段,A位于浦东河干流的上游点,C位于新厝河与浦东河干流的汇集处,D位于淌洋河与浦东河干流的汇集处,E位于浦东河下游的人工湖处。
采集A,B,C,D,E共5处0~10 cm深度的表层底泥,测定有机质(OM)、总氮(TN)、总磷(TP)和重金属(Cu,Zn,Ni,Pb),测定方法见表 1。底泥的颗粒组成、黏土矿物组成和界限含水率见表 2,采用筛分法和密度计法对底泥进行颗粒分析,分别采用Casagrande法和搓条法测定液限wL和塑限wP,采用X射线衍射法测定底泥黏粒中主要矿物成分,包括伊利石(I)、高岭石(K)、绿泥石(C)、蒙脱石(S)的含量。
表 1 本研究底泥污染物测定方法Table 1. Method for determination of pollutants in sediments序号 测试项目 测试方法 试验标准 1 OM 烧失量法 ASTM D2974 2 TN 凯氏法 HJ717—2014 3 TP 钼锑抗分光光度法 HJ 632—2011 4 重金属 ICP-MS法 US EPA 3050B 表 2 本研究底泥颗粒组成、黏土矿物组成及界限含水率Table 2. Particle composition, clay mineral composition and atterberg limits of sediments(%) 底泥 颗粒组成 矿物组成 wL wP Clay Silt Sand I K C S A 42.4 47.5 10.2 19 48 24 9 79.2 35.3 B 12.3 78.8 8.9 22 59 19 0 44.0 31.5 C 12.3 82.9 4.8 33 47 20 0 38.5 23.3 D 26.0 64.0 10.0 27 53 20 0 83.6 35.0 E 28.3 62.4 9.3 24 46 30 0 111.9 44.0 表 3显示了来源于文献的具有不同颗粒级配、界限含水率、黏土矿物组成,以及不同污染物含量的河道底泥数据,结合本文的试验数据,分析底泥物理–矿物成分–污染性状的关联性。
表 3 不同文献收集的底泥数据Table 3. Database of sediment pollutants compiled from literatures序号 颗粒级配 界限
含水率黏土矿物 污染物 主要污染来源 参考文献 OM TN TP Cu Zn Ni Pb 1 √ — — √ √ √ — — — — — 魏岚等[6] 2 √ — — √ √ √ — — — — — Xia等[7] 3 — — — √ √ √ — — — — 生活污水 孙广垠等[8] 4 √ — — — √ √ — — — — 废水、肥料 余成等[9] 5 √ — — — — — √ √ √ √ 废水 El-Sayed等[10] 6 √ — — — — — √ √ √ √ 养殖场 Wang等[11] 7 — — √ √ — — — — — — — Khim[12] 8 — — √ √ — — — — — — — Andrade等[13] 9 — — — √ — — √ √ √ √ 生活污水 Nguyen等 [14] 10 — — — √ — — √ √ √ √ 生活污水 牛红义等[15] 11 — — — — √ √ √ √ √ √ 废水 严玉林[16] 12 — √ — √ — — — — — — — 徐日庆等[17] 13 — √ — √ — — — — — — — Stanchi等 [18] 14 — √ — — — — √ — — — — Phanija等 [19] 15 — √ — — — — — √ — — — 储亚等[20] 16 — √ — — — — — — — √ — Ayodele等 [21] 17 — √ — — — — √ √ — √ — 吕伟豪[22] 2. 底泥试验结果与分析
2.1 污染物与底泥颗粒级配的关系
(1)有机质与底泥颗粒级配的关系
底泥有机质与细颗粒含量的关系绘制于图 2中,可以发现底泥的细颗粒与有机质之间具有较强的相关性,有机质含量随细颗粒含量的增加而增加。底泥有机质含量随细颗粒含量的关系曲线的斜率不同,斜率越大表明底泥中的细颗粒对有机质的吸附作用越强。细颗粒具有较大的比表面积,有利于对有机质的吸附和聚集。底泥中有机质不仅与颗粒级配有关,还有底泥附近的污染源和环境有关。本研究河道底泥位于城市居民区,周围有大量排污管道将居民生活废水排入河道中,使得底泥中含有较高的有机质,文献[6,7]的样品分别取自水库底泥和海湾底泥中,周围没有人为污染源,由于水库的流动性小于海湾,使得水库底泥的有机质含量>海湾底泥的有机质含量。
(2)总氮、总磷与底泥颗粒级配的关系
底泥细颗粒含量与总氮、总磷含量的关系绘制于图 3中,可以发现同一河道底泥的总氮、总磷含量随底泥细颗粒含量的增加而增加,这与有机氮、有机磷易于吸附在细颗粒上有关。底泥周边环境,黏土矿物成分的不同造成了总氮、总磷含量与细颗粒含量的关系曲线的斜率不同。
(3)重金属含量与底泥颗粒级配的关系
底泥细颗粒含量与重金属含量的关系绘制于图 4中。由图 4可见,底泥中重金属的积累受底泥颗粒级配的影响,底泥的Cu,Zn,Ni含量随底泥细颗粒含量的增加而增加,由于粒度影响底泥的比表面积、孔隙体积以及活性组分,使得底泥细颗粒具有强吸附能力,有利于重金属元素的汇集。同时,底泥粒径越细,所含有机质也越多,对重金属的吸附络合作用也越强。本研究的底泥重金属含量较高,与沿河汽车修理厂等工厂废水的长期污染有关,且本研究底泥中的有机质含量较高,使得重金属元素大量累积。
2.2 有机质与黏土矿物的关系
为了研究黏土矿物组成与底泥有机质含量的关系,选取文献[12,13]黏粒含量在40%±1%范围内的底泥,其有机质含量与各黏土矿物组成的关系绘制于图 5中,可以发现,对于相同黏粒含量的底泥,不同黏土矿物对有机质含量有不同的影响,蒙脱石与高岭土对有机质的吸附和储存能力较强,且与有机质含量呈现明显的正相关,伊利石与有机质含量呈弱负相关。虽然伊利石的比表面积大于高岭石,但是本研究对比发现高岭石含量高的底泥中有机质含量较伊利石多,其原因可能是黏土矿物对有机质存在选择性的吸附,不同的黏土矿物保存着不同的有机组分,高岭石易于吸附有机质中的—CH2基团,而在底泥中含量较多有机质是胡敏酸,—CH2是其主要官能团,易与高岭石吸附结合。这一现象有待今后积累更多的试验数据,开展进一步的探讨。
2.3 底泥各污染物之间的关系
(1)底泥有机质与总氮、总磷的关系
底泥有机质含量与总氮、总磷含量的关系绘制于图 6中。由图 6(a)可以发现底泥总氮含量随有机质含量的增加而增加,由图 6(b)可以看出底泥中总磷含量随有机质含量增加的规律不明显,有机质与总磷含量的相关性较有机质与总氮含量的相关性弱。底泥中的氮素有95%以上存在于有机物质中,因此总氮含量与有机质含量呈显著正相关。
(2)底泥有机质与重金属含量的关系
底泥有机质与重金属含量的关系绘制于图 7中。可以发现不同底泥中的重金属含量差异巨大,某些重金属元素的含量甚至相差100倍以上,本研究底泥的重金属含量明显远大于文献[14,15]底泥,这与重金属污染源有关,本研究底泥河道沿线有不锈钢加工厂、汽修厂等众多污染源,造成底泥中重金属污染严重,文献[14,15]底泥的主要污染源为生活污水,因此文献[14,15]底泥的重金属污染程度较本研究底泥轻。底泥重金属含量随有机质含量的增加而增加,不同文献底泥关系曲线的斜率不同,重金属污染源对曲线斜率的大小影响很大。
(3)底泥总氮、总磷含量与重金属含量的关系
底泥总氮、总磷含量与重金属含量的关系绘制于图 8中。可以发现底泥的重金属含量随底泥总氮、总磷含量的增加而增加。有机质与总氮总磷的同源性,以及有机质对重金属的吸附和络合作用,使得底泥重金属含量与底泥总氮总磷含量同样具有正相关的关系。
2.4 污染物与底泥界限含水率的关系
(1)有机质与底泥界限含水率的关系
底泥有机质含量与液限、塑限和塑性指数的关系绘制于图 9中,可以发现底泥中液限、塑限及塑性指数随底泥有机质含量的增加而增加,底泥的有机质含量与液塑限及塑性指数之间具有较强的相关性,液塑限、塑性指数与有机质的关系式列于图中。有机物对液限和塑限的影响是通过改变土颗粒结合水膜的厚度来实现的,有机物具有较高的比表面积和较强的持水能力,可吸附在黏土矿物表面,形成较厚的结合水膜,从而提高底泥的液塑限。
(2)重金属含量与底泥界限含水率的关系
底泥重金属含量与液限、塑限及塑性指数的关系绘制于图 10中。可以发现文献[19~22]底泥的液塑限随着重金属含量的增加而减小。重金属离子对底泥液塑限的影响主要是引起了黏土矿物的聚集和双电层厚度的改变。本研究底泥的液塑限及塑性指数则随着重金属含量的增加而增加,这是因为本研究底泥中有机质含量较高,而文献[19~22]底泥中几乎不含有机质,重金属含量会随着有机质含量的增加而增加,且有机质对液塑限的增加作用大于重金属对液塑限的减小作用。
3. 底泥污染物及物理性质相关性分析
对底泥中的污染物与物理性质数据进行相关性分析,为了减少底泥所处环境因素对相关性分析的影响,对底泥污染物和物理性质按区域进行分析,后取其相关系数平均值。底泥污染物与物理性质的相关系数,如表 4所示。有机质、总氮及Cu,Zn,Ni重金属含量与底泥中黏粒含量和粉粒含量的相关性较高,相关系数均大于0.5以上;污染物与黏土矿物相关性不强,相关系数均小于0.5,这是由于黏土矿物在整个底泥颗粒中所占比重较小,影响力有限;有机质含量与各污染物含量之间的相关系数均较高,表明底泥污染性状与有机质含量密切相关;有机质对底泥液限、塑限、塑性指数的相关系数分别为0.915,0.916,0.797,这表明有机质对底泥的物理性质有着重要的影响,其他污染物与底泥物理性质的相关系数均较低,对底泥物理性质的影响较小。
表 4 底泥污染物及物理性质的相关系数Table 4. Correlation coefficients of sediment pollutants and physical properties污染物与底泥颗粒级配 污染物与黏土矿物 污染物与污染物 污染物与底泥物理性质 Clay OM 0.231 Illite OM -0.099 OM TN 0.809 OM wL 0.915 Silt OM 0.524 Kaolinite OM -0.185 OM TP 0.456 OM wP 0.916 Sand OM -0.485 Chlorite OM -0.194 TN TP 0.623 OM IP 0.797 Clay+Silt OM 0.717 Smectite OM 0.020 OM Cu 0.636 TN wL 0.254 Clay TN 0.686 Illite TN -0.249 OM Zn 0.794 TN wP 0.242 Clay TP 0.439 Illite TP -0.356 OM Ni 0.490 TN IP 0.161 Silt TN -0.026 Kaolinite TN -0.216 OM Pb 0.777 TP wL -0.009 Silt TP -0.097 Kaolinite TP -0.269 TN Cu 0.452 TP wP -0.005 Sand TN -0.763 Chlorite TN 0.485 TN Zn 0.603 TP IP -0.033 Sand TP 0.321 Chlorite TP 0.264 TN Ni 0.511 Cu wL -0.342 Clay+Silt TN 0.763 Smectite TN 0.197 TN Pb 0.433 Cu wP 0.300 Clay+Silt TP 0.321 Smectite TP 0.054 TP Cu 0.335 Cu IP -0.350 Clay Cu 0.355 Illite Cu -0.216 TP Zn 0.577 Zn wL -0.331 Clay Zn 0.363 Illite Zn -0.382 TP Ni 0.203 Zn wP -0.317 Clay Ni 0.335 Illite Ni -0.270 TP Pb 0.501 Zn IP 0.275 Clay Pb 0.335 Illite Pb -0.286 Cu Zn 0.655 Pb wL 0.067 Silt Cu 0.572 Kaolinite Cu -0.229 Cu Ni 0.610 Pb wP -0.365 Silt Zn 0.558 Kaolinite Zn -0.329 Cu Pb 0.551 Pb IP 0.112 Silt Ni 0.639 Kaolinite Ni -0.245 Zn Ni 0.539 Silt Pb 0.006 Kaolinite Pb 0.047 Zn Pb 0.729 Sand Cu -0.555 Chlorite Cu 0.258 Ni Pb 0.406 Sand Zn -0.593 Chlorite Zn 0.482 Sand Ni -0.641 Chlorite Ni 0.445 Sand Pb -0.292 Chlorite Pb 0.362 Clay+Silt Cu 0.554 Smectite Cu -0.137 Clay+Silt Zn 0.591 Smectite Zn -0.214 Clay+Silt Ni 0.635 Smectite Ni -0.159 Clay+Silt Pb 0.295 Smectite Pb 0.103 4. 结论
基于福州晋安东区五处河道底泥系列试验研究结果,结合搜集的独立试验数据,进行了底泥物理-矿物成分–污染性状关联性分析,得出4点结论。
(1)在细颗粒含量较高的底泥中,有利于污染物的吸附积累,随着细颗粒含量的增加,污染物含量近似呈线性增加的趋势。而在砂粒含量较高的底泥中,则缺少这种吸附能力,底泥中污染物含量低。
(2)底泥中不同黏土矿物对有机质含量有不同的影响,蒙脱石与高岭土与有机质含量呈现明显的正相关。
(3)底泥中有机质与氮磷营养盐一般具有同源性,同时有机质对重金属具有络合作用,底泥中有机质含量与总氮总磷含量,有机质含量与重金属含量,总氮总磷含量与重金属含量,均具有良好的线性关系。
(4)底泥中的液限、塑限和塑性指数随着底泥有机质含量的增加而增加,相比于有机质,重金属对底泥界限含水率的影响较小。
-
表 1 高精密科学装置对基础振动的控制标准[2, 8, 10]
Table 1 Control criteria of vibration of foundation of high-precision scientific equipments [2, 8, 10]
装置名称 1~100 Hz竖向振动位移均方根/nm 上海光源 < 150(安静时段),< 300(嘈杂时段) 软X射线自由电子激光 < 500 硬X射线自由电子激光 < 150 北京高能同步辐射光源 < 25 注:安静时段系指0:00—4:00;嘈杂时段系指4:00—24:00。 表 2 北京光源工程地基土层分布及参数
Table 2 Soil profile and parameters at HEPS site
土层 厚度/m 密度
/(kg·m-3)压缩波速/(m·s-1) 剪切波速/(m·s-1) 泊松比 弹性
模量
/MPa①素填土 1 1910 636 173 0.460 166.9 ②1卵石 11 2200 952 344 0.425 742.0 ②2卵石 6 2200 1451 492 0.435 1528.4 ④全风化基岩 8 2250 1962 777 0.407 3822.5 ⑤强风化基岩 4 2350 2253 1011 0.374 6600.7 ⑥中风化基岩 21 2590 3616 1547 0.388 17206.8 ⑦微风化基岩 49 2680 4302 1936 0.373 27583.3 -
[1] 夏禾. 交通环境振动工程[M]. 北京: 科学出版社, 2010. XIA He. Traffic Induced Environment Vibrations and Controls[M]. Beijing: Science Press, 2010. (in Chinese)
[2] 余宽原. 考虑土体弹性各向异性的环境微振动分析研究[D]. 上海: 同济大学, 2021. YU Kuan-yuan. Analysis of Environmental Vibrations in Transverse Isotropic Elastic Ground[D]. Shanghai: Tongji University, 2021. (in Chinese)
[3] 电子工业防微振工程技术规范: GB50176)—2015[S]. 北京: 中国计划出版社, 2015. Technical Code for Anti-Microvibration Engineering of Electronics Industry: GB50176—2015[S]. Beijing: China Planning Press, 2015. (in Chinese)
[4] AMIRIKAS R, EHRLICHMANN H, BIALOWONS W, et al. Ground motion and comparison of various sites[R]. Hamburg: Deutsches Elektronen-Synchrotron, 2005.
[5] 黄茂松. 上海光源工程微振动数值模拟计算研究报告[R]. 上海: 同济大学, 2005. HUANG Mao-song. Research Report on Numerical Simulations of Micro-vibration in Shanghai Synchrotron Radiation Facility (SSRF) [R]. Shanghai: Tongji University, 2005. (in Chinese)
[6] 黄茂松, 任青, 周仁义, 等. 层状地基中瑞利波随深度的衰减特性[J]. 岩土力学, 2009, 30(1): 113-117, 122. doi: 10.3969/j.issn.1000-7598.2009.01.018 HUANG Mao-song, REN Qing, ZHOU Ren-yi, et al. Attenuation characters of Rayleigh wave in layered soils[J]. Rock and Soil Mechanics, 2009, 30(1): 113-117, 122. (in Chinese) doi: 10.3969/j.issn.1000-7598.2009.01.018
[7] 周仁义. 分层地基中环境微振动的衰减特性[D]. 上海: 同济大学, 2006. ZHOU Ren-yi. Attenuation of Environmental Vibration with Depth[D]. Shanghai: Tongji University, 2006. (in Chinese)
[8] 冯忠俊. 上海自由电子激光工程地基振动测试研究及数值模拟[D]. 上海: 中国科学院上海应用物理研究所, 2009. FENG Zhong-jun. Ground Vibration Measurement and Numerical Simulation on Shanghai Free Electron Laser Project[D]. Shanghai: Shanghai Institute of Applied Physics, 2009. (in Chinese)
[9] 岳建勇, 姚激, 周春, 等. 软土地区特殊精密装置基础微变形控制设计与实践[J]. 岩土工程学报, 2010, 32(增刊2): 267–270. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2010S2066.htm YUE Jian-yong, YAO Ji, ZHOU Chun, et al. Design and practice of differential settlement control for soft foundation by using a special sensitive facility[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(S2): 267–270. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2010S2066.htm
[10] 岳建勇. 软土地基精密装置基础微振动控制技术与工程应用[J]. 建筑科学, 2020, 36(增刊1): 57–67. https://www.cnki.com.cn/Article/CJFDTOTAL-JZKX2020S1011.htm YUE Jian-yong. Technology of the micro-vibration control for the foundation of a special sensitive facility on soft ground with engineering applications[J]. Building Science, 2020, 36(S1): 57–67. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JZKX2020S1011.htm
[11] RICHART F E, HALL J R, WOODS R D. Vibrations of Soils and Foundations[M]. New Jersey: Prentice-Hall Inc., 1970.
[12] YU K Y, GU X Q, HUANG M S, et al. Experimental, numerical and analytical studies on the attenuation of maglev train-induced vibrations with depth in layered soils[J]. Soil Dynamics and Earthquake Engineering, 2021, 143: 106628. doi: 10.1016/j.soildyn.2021.106628
[13] THOMSON W T. Transmission of elastic waves through a stratified solid medium[J]. Journal of Applied Physics, 1950, 21(2): 89–93. doi: 10.1063/1.1699629
[14] HASKELL N A. The dispersion of surface waves on multilayered media[J]. Bulletin of the Seismological Society of America, 1953, 43(1): 17–34. doi: 10.1785/BSSA0430010017
[15] 北京市勘察设计研究院有限公司. 高能同步辐射光源项目岩土工程勘察报告[R]. 北京: 北京市勘察设计研究院有限公司, 2018. BGI Engineering Consultants Ltd. The Geotechnical Investigation Report of HEPS Project[R]. Beijing: BGI Engineering Consultants Ltd, 2018. (in Chinese)
[16] 吴世明. 土介质中的波[M]. 北京: 科学出版社, 1997. WU Shi-ming. Wave Propagation in Soils[M]. Beijing: Science Press, 1997. (in Chinese)
[17] ZERWER A, CASCANTE G, HUTCHINSON J. Parameter estimation in finite element simulations of Rayleigh waves[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2002, 128(3): 250–261. doi: 10.1061/(ASCE)1090-0241(2002)128:3(250)
[18] 杨永斌. 高速列车所引致之土壤振动分析[R]. 台湾: 台湾大学, 1996. YANG Yong-bin. Analysis on Soil Vibration Induced by High-Speed Railways[R]. Taiwan: Taiwan University, 1996. (in Chinese)
-
期刊类型引用(1)
1. 张闯,王苒,尚晓雨,杨阳,陈龙,董硕,刘泽,杨科. 水域地表基质调查方法探索——以衡水湖为例. 国土资源导刊. 2023(02): 109-113 . 百度学术
其他类型引用(2)
-
其他相关附件
-
PDF格式
审稿意见与作者答复 点击下载(159KB)
-