环境科学学报  2020, Vol. 40 Issue (8): 2785-2794
北运河上游合流制管网溢流污染特性研究    [PDF全文]
海永龙1,2,3, 郁达伟1,2, 刘志红4, 钱晓辉4, 柳蒙蒙1,2,3, 狄斐1,2, 王拓1,2,3, 岳文慧1,2,3, 胡大洲1,2,3, 梁玉帅1,2, 魏源送1,2,3    
1. 中国科学院生态环境研究中心, 环境模拟与污染控制国家重点联合实验室, 北京 100085;
2. 中国科学院生态环境研究中心, 水污染控制实验室, 北京 100085;
3. 中国科学院大学, 北京 100049;
4. 昌平区水务局, 北京 102200
摘要:以北运河沙河水库周边合流制管网为研究对象,选取浊度为主要指示指标,通过监测典型溢流排口,考察了2019年4次合流制管网溢流污染物的变化过程,研究了降雨事件之间的干旱天数和降雨强度对溢流污染的影响,并分析浊度和典型污染物之间的相关性.结果表明,汛期降雨强度较大、历时较长,是溢流事件发生的主要时期.当单场次累积降雨量达到15 mm和单场次平均降雨强度达到1.4 mm·h-1时,溢流开始发生.其中,4月24日(第一次溢流)和汛期7月22日(干旱天数为23 d,降雨强度最强)的初期溢流污染最为严重,这两次溢流污水中TN、TP、TCOD都与浊度显著相关(p < 0.01);5月26日(非汛期典型降雨)的溢流污染物与浊度相关性不显著(p>0.05),污染负荷较低.这说明北运河上游沙河水库周边合流制管网溢流污染主要受降雨强度和干旱天数的影响.因此,以年总溢流污染负荷为控制目标时,应优先控制年度初次溢流和主汛期初次溢流的初期溢流污染.
关键词北运河上游    合流制溢流污染    降雨强度    浊度    
Characteristics of the combined sewer overflows pollution in the upper North Canal of Beijing
HAI Yonglong1,2,3, YU Dawei1,2, LIU Zhihong4, QIAN Xiaohui4, LIU Mengmeng1,2,3, DI Fei1,2, WANG Tuo1,2,3, YUE Wenhui1,2,3, HU Dazhou1,2,3, LIANG Yushuai1,2, WEI Yuansong1,2,3    
1. State Key Joint Laboratory of Environment Simulation and Pollution Control, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085;
2. Laboratory of Water Pollution Control Technology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085;
3. University of Chinese Academy of Sciences, Beijing 100049;
4. Changping Water Authority, Beijing 102200
Received 16 January 2020; received in revised from 6 March 2020; accepted 10 March 2020
Abstract: A study was carried out to investigate the combined sewer overflows (CSOs) pollution at the upper North Canal, and turbidity was selected as the main indicator for CSOs pollution. Through investigating changes of CSOs pollutants during 4 rainfall events in 2019, the effect of dry days and rainfall intensity on CSOs pollution was analyzed, as well as the correlation between turbidity and typical pollutants in this study. Results suggested that the CSOs mainly occurred in the rain season due to the higher rainfall intensity and longer rainfall duration. The CSOs occurred at the accumulated precipitation over 15 mm and the average intensity up to 1.4 mm·h-1. The heavier polluted initial CSOs occurred at the first overflow on Apr. 24 and the overflow on July 22 due to the intensive rainfall after 23 days rain free period. During these two CSOs events, concentrations of TN, TP, and TCOD significantly correlated with the turbidity (p < 0.01). However, there is no significant correlation (p>0.05) observed between concentrations of TN, TP, and TCOD and turbidity due to low pollutants loading of the CSOs occurred on May 26. In a word, the CSOs pollution was mainly affected by rainfall intensity and rain free period in the upper North Canal, Beijing. Therefore, it is of great importance to control the initial overflow pollutants of both the first CSOs in the year and the first CSOs in the main flood season, with aiming to control annual CSOs pollution loading.
Keywords: upper North Canal    combined sewer overflows    rainfall intensity    turbidity    
1 引言(Introduction)

合流制管网(Combined Sewer System, CSS)主要收集城镇生活污水、雨水乃至工业废水(邢丽贞等, 2011).旱季时管网流量较小, 污水中的颗粒物沉积在管网中并吸附富集有机物等污染物(Tait et al., 2003);雨季时管网流量较大, 但沉积物会影响管网的输送能力, 形成合流制管网溢流(Combined Sewer Overflows, CSOs), 导致管网中部分污水携带管道污泥直接溢流排放到附近的受纳水体, 从而产生溢流污染.溢流污染物包含颗粒态污染物, 特别是有机物、氮磷富营养化物质及潜在的病原微生物乃至有毒有害物质(杨雪等, 2008; Holeton et al., 2011; Passerat et al., 2011; Madoux-Humery et al., 2013), 这些污染物直接进入受纳水体, 由此造成水环境污染.Franco等(2018)监测了普布亚诺湖的167次CSOs, 发现河流中39%的总磷来源于CSOs;Casadio等(2010)考察了CSOs对河流城市段的影响, 发现即使在低降雨强度下, CSOs也会造成河流水质的明显变化.由此可见, 合流制管网溢流是雨季城市水体污染的关键潜在风险.为了合理管控CSOs对水体的污染风险, 美国环境保护署(EPA)在《清洁水法》等相关法律中制订了专门的CSOs章节(Salzman et al., 2001; Postel et al., 2005), 明确了溢流污染特征分析规程和管控目标, 我国目前也形成了技术导则等指导政策, 但还未出台明确的针对性法律.

合流制管网溢流事件的发生时间受降雨影响大, 较难把握, 污染监测和特征识别是其中的难点.现有研究表明, 浊度与颗粒态污染物有较强的相关性(Lacour et al., 2009).考虑到浊度易于监测, 通过监测浊度对溢流污染进行定量评估已成为一种广受认可的方法(Viviano et al., 2014; Bersinger et al., 2015; Zhang et al., 2018).北京市合流制管网约有766 km(潘国庆, 2008), 新近发布的《北京市进一步加快推进城乡水环境治理工作三年行动方案(2019年7月—2022年6月)》明确要求加强对合流制溢流污染的管控.因此, 本研究以北运河上游沙河水库北岸沿岸的城郊合流制排水管网为研究对象, 并选取浊度为主要指示指标, 通过监测典型溢流排口, 考察2019年4次合流制管网溢流污染物的变化过程, 探究降雨事件之间的干旱天数和降雨强度对溢流污染的影响, 并分析浊度与典型污染物之间的相关性, 以期为合流制管网溢流污染控制提供科学依据.

2 材料和方法(Materials and methods)) 2.1 研究区域概况

沙河水库位于北运河上游的源头区(图 1a), 水源主要来源于东沙河、南沙河和北沙河的汇水, 总流域面积为1125 km2 (张伟等, 2012), 其中, 东沙河、南沙河和北沙河的流域面积分别为265、263、597 km2 (王辰等, 2002; 刘博等, 2011; 朱光旭等, 2013).沙河水库为沙河闸控制的河道式水库, 水库总面积为1.8 km2, 总库容量为2045×104 m3(孙文等, 2019).沙河水库的汇水区是北运河源头的山水林田湖生态涵养区, 其排水分区也是区域人类活动集中区, 对北运河水质改善和水生态修复至关重要.

图 1 研究区域及监测点位(a.研究区域, b.流量、降雨监测点) Fig. 1 Investigated region and monitoring point (a.Study area, b.Flow and rainfall monitoring points)

根据课题组于2017—2018年的现场调查结果, 库区北沙河和南沙河的沿岸分别有9个、11个合流制溢流管网溢流口, 管网管径为300~1200 mm, 管网最小坡度为3‰.根据2017—2018年的现场监测结果, 本研究选择典型溢流口P(116°20′14.79″N, 40°00′26.82″E)作为在线连续监测位点.北沙河南岸总共有3个溢流口, 其中, 监测位点P所在排水分区的合流制溢流管网主要服务于巩华城所在区域的居民区, 其服务面积约为0.77 km2.巩华城所在区域属于快速发展的规划高端综合服务区, 该排水分区的服务(汇水区)面积为1.8 km2, 不透水地表面比例为43%~55%, 干管为DN500~DN1200, 设计截留倍数为2倍.监测点采用槽堰式截流井, 实测持续截留能力为0.162 m3·s-1(图 1b).

2.2 研究方法 2.2.1 采样和监测方法

本研究分别监测了2019年4月24日、5月26日、7月22日、7月28日4场降雨和溢流过程.降雨时, 采用自动采样器(型号:ZY-WLG-II)从管网和溢流口采集出流污水, 并保存于500 mL的聚乙烯瓶中, 在24 h内测定完毕.降雨期间的样品采集频率主要分为4个阶段:①降雨前1 h, 每5 min采集一个样品;②降雨中期1~2 h, 每10 min采集一个样品;③降雨后期2~5 h, 每30 min采集一个样品;④降雨末期5 h以上, 每2 h采集一个样品, 采样时长根据降雨历时和溢流过程确定.7月22日和7月28日的降雨信息通过现场安装的雨量计(型号:QS-3 000, 邯郸开发区清易电子科技有限公司)监测, 雨量统计频率为5 min·次-1, 而由于现场仪器故障, 未能监测到4月24日和5月26日的降雨, 故本研究采用中国气象局提供的昌平站点(ID 54499)降雨量数据.

图 1b所示, 流量数据监测采用远传式明渠流量计系统(北京泽源伟业环保科技有限公司安装), 其中, 溢流量监测采用薄壁三角堰, 管网流量监测采用矩形薄壁堰, 通过堰口模型尺寸, 建立水位-流量数学模型(水利部, 1991), 流量数据采集频率为10 min·次-1.

2.2.2 污染物分析测试方法

样品的采集、保存和监测方法均参考《水和废水监测分析方法指南》(中华人民共和国国家质量监督检验检疫总局等, 2002), 采用国标法(环境保护部, 2009)测定样品的TN、NH3-N、TP、COD和浊度, 其中, TN测定采用过硫酸钾-紫外分光光度法, NH3-N测定采用纳氏试剂分光光度法, TP测定采用钼酸铵分光光度法, 所用仪器均为紫外可见光光度计(TU-1901型, 北京普析通用仪器有限责任公司), COD测定采用分光光度方法(DR2800, HACH Inc., USA), 浊度采用便携式浊度计(Turb 550, WTW GmbH, Germany)测定.

2.3 数据处理及分析方法

本研究中降雨场次的划分标准:从某时刻降雨开始算起在2 h内降雨量≥1 mm, 即为一场降雨, 场次平均降雨强度为一场降雨量与降雨历时之比(Bersinger et al., 2018), 最大降雨强度为采集频率下的最大降雨量, 瞬时降雨强度为采集频率下的降雨量.利用SPSS 25.0软件对污染物指标与浊度进行Pearson相关性分析.

3 结果与讨论(Results and discussion) 3.1 降雨历时变化特征

2019年4次降雨的溢流监测结果表明, 当单场次累积降雨量超过15.0 mm和单场次平均降雨强度>1.4 mm·h-1后才开始发生溢流.因此, 当单场次累积降雨量 < 15.0 mm和单场次平均降雨强度 < 1.4 mm·h-1时, 合流制管网能够截留区域降水, 通常不会发生溢流.按照《降水等级标准(GB/T 28592—2012)》划分(国家气象中心, 2012), 4月24日、5月26日、7月22日、7月28日的降雨分别为中雨、中雨、暴雨、暴雨.这4次降雨特征如表 1图 2所示, 4月24日的降雨强度小、历时长, 5月26日的降雨强度小、历时短, 7月22日的降雨强度大、历时短, 7月28日的降雨强度小、历时长.

表 1 2019年4场降雨特征 Table 1 Characteristics of the four rainfalls in 2019

图 2 4场降雨的历时变化 Fig. 2 Changes of the four rainfall events along with time
3.2 管道和溢流的流量历时变化

图 23所示, 2019年这4场降雨的降雨-流量-溢流连续在线监测结果表明, 当4月24日、5月26日、7月22日、7月28日的累积降雨量分别达到约16.0、23.6、15.0和20.8 mm时, 溢流开始发生, 从开始降雨到发生溢流历时分别约为10.5、5.0、0.83和15.0 h, 管道流量的峰值分别滞后于瞬时降雨强度的峰值30、260、60和30 min.在溢流发生前的4场降雨平均降雨强度分别为1.5、4.7、16.6、1.4 mm·h-1, 因此, 溢流前的降雨强度会影响从降雨开始到溢流发生所需时间, 且溢流前的降雨强度越大, 所需时间越小.这种趋势与降雨产汇流过程的趋势一致, 降雨强度越大, 地面径流形成越快, 导致合流制管网溢流也发生的更快.根据图 3的监测结果, 开始溢流时管网流量瞬时值为0.162~1.531 m3·s-1, 瞬时值波动较大, 可能与管网压力波动和下游管网负荷等情况有关;溢流后较稳定的管网流量临界值为0.162~0.280 m3·s-1.

图 3 次累积降雨量与管道、溢流的流量之间的变化关系 Fig. 3 Relationship between pipeline flow, overflow flow and accumulated rainfall per event

值得注意的是, 合流制管网溢流的初期溢流污染不同于降雨产汇的初雨径流污染.本研究的监测结果表明:从开始降雨到管道流量增加(t1)和开始溢流(t2), 分别经过了t1=27~400 min和t2=110~900 min, 这说明研究区合流制管网通常可以截留大部分初雨径流(前15~30 min产流), 但由于管网截留能力、污水处理厂处理能力等限制, 导致初期溢流污染仍不可避免.如图 3所示, 管道流量迅速达到峰值导致溢流发生, 形成了初期溢流污染.因此, 可通过采取海绵城市的屋面雨水断流、植草沟缓滞等措施, 合理减缓管道流量的增速, 使初期溢流污染有望得到控制.

3.3 管道污水和溢流的水质变化过程

图 4~7展示了这4场降雨过程中合流制管网溢流端口的污染物浓度输出过程, 在溢流流量升高的过程中, 溢流污水的污染物浓度逐渐升高, 并在溢流中后期达到稳定.图 4是2019年4月24日第一次发生溢流事件时溢流口断面的流量和污染物输出浓度情况, 虽然两次降雨事件之间的干旱天数只有5 d, 但4月19日降雨量仅为6.5 mm, 而降雨历时却为13 h, 此次降雨强度小, 对管道污泥和地面的冲刷较弱, 没有发生溢流.而4月24日的降雨却导致了溢流事件, 由于降雨初期雨水对管道污泥和地面的冲刷, 在溢流流量达到峰值之前, 污染物浓度随着溢流流量的升高而升高;溢流流量达到溢流峰值以后, 由于稀释作用, 溢流中污染物浓度迅速降低, 780 min以后污染物浓度趋于稳定.由图 5可知, 5月26日的降雨发生了两次溢流, 第一次溢流过程中, 虽然溢流流量逐渐升高, 但溢流中TCOD和SCOD的浓度变化很小, 浊度逐渐减小, TN、NH3-N和TP的浓度逐渐升高, 由于此次降雨前的干旱天数为5 d, 因此, 污染物浓度变化幅度较小;同时, 由于雨水稀释作用, 导致第二次溢流过程中各污染物浓度均在溢流流量达到峰值之前降低, 并在流量达到峰值后随着流量的降低而升高.如图 6所示, 由于此次降雨(7月22日)前的干旱天数为23 d, 故溢流污水中各污染物浓度在溢流初期相对较高, 110 min以后污染物浓度变化相对较小.如图 7所示, 由于溢流时间较短, 采样个数比较少, 污染物浓度未有明显的变化趋势.根据监测点P的溢流量在线监测结果计算, 4次典型降雨溢流量占全年溢流量的72%, 根据溢流量-浓度时间分布关系合理外推的溢流污染负荷占比约为75%~83%.

图 4 4月24日溢流口断面的流量和污染物输出浓度示意图 Fig. 4 Changes of cross-section flow and pollutants output at the overflow port on April 24

图 5 5月26日溢流口断面流量和污染物输出浓度示意图 Fig. 5 Changes of cross-section flow and pollutants output at the overflow port on May 26

图 6 7月22日溢流口断面流量和污染物输出浓度示意图 Fig. 6 Changes of cross-section flow and pollutant output on the overflow port on July 22

图 7 7月28日溢流口断面流量和污染物输出示意图 Fig. 7 Changes of cross-section flow and pollutant output on the overflow port
3.4 溢流污水污染物特征及其相关性分析

为分析浊度与污染物指标之间的相关性, 本研究对这4次降雨进行了主成分分析(PCA), 结果如图 8所示.由图可知, TP、TN、TCOD与浊度的相关性较好, 而溶解性污染物SCOD、NH3-N与浊度之间的相关性较差, 这表明浊度与颗粒态污染物有较强的关系.

图 8 4次降雨中污染物与浊度的主成分分析 Fig. 8 Principal component analysis of pollutants and turbidity in four rainfalls

本研究对溢流过程中污染物形态与浊度之间进行了Person相关性分析(表 2), 并对4月24日和7月22日两次溢流事件进行了线性回归分析(图 9~图 10).分析结果表明, 除SCOD之外, 两次溢流的其他污染物与浊度的变化趋势较为相似, 例如, 两次溢流中TN、NH3-N、TP、TCOD与浊度相关性显著(p < 0.01), 因此, 对于溢流污染物的控制, 控制浊度就相当于控制TN、NH3-N、TP、TCOD.但5月26日的降雨事件中, TN、NH3-N、TP、TCOD、SCOD与浊度之间的相关性不显著(p>0.05)(表 2), 深入分析认为, 此次降雨前干旱天数仅为7 d, 干旱天数较短, 在管道和地面上沉积的颗粒较少, 降雨后对地面和管道冲刷产生的污染物贡献较小, 而且此次降雨事件平均降雨强度仅为4.7 mm·h-1, 降雨强度较小, 雨水对地面和管道污泥的冲刷强度较小, 故污染物主要来源于管道污水, 这与现有研究得出的结果相似(周美成等, 2015).

表 2 2019年4次溢流事件污染物浓度与浊度之间的相关性 Table 2 Correlation between pollutants concentration and turbidity at 4 overflow events in 2019

图 9 4月24日溢流污水污染物浓度与浊度之间的拟合结果(a.TCOD、SCOD, b.TP、TN、NH3-N) Fig. 9 Correlation between turbidity and TCOD, SCOD(a) and TP, TN, NH3-N(b) of overflow sewage on April 24

图 10 7月22日溢流污水污染物浓度与浊度之间的拟合结果(a.TCOD、SCOD, b.TP、TN、NH3-N) Fig. 10 Correlation between turbidity and TCOD, SCOD(a) and TP, TN, NH3-N(b) of overflow sewage on July 22
4 结论(Conclusions)

1) 2019年4次典型溢流污染特征的调查表明, 当单场次累积降雨量达到15 mm和单场次平均降雨强度达到1.4 mm·h-1时, 溢流开始发生;溢流发生受溢流前降雨强度影响, 溢流前降雨强度越大, 溢流发生越快.

2) 初期溢流污染浓度主要受干旱天数和降雨强度的影响.溢流初期的污染负荷最高;干旱天数越长, 初期溢流污染物浓度越高, 降雨强度会影响冲刷.因此, 控制年度初次溢流、汛期初次溢流对于年溢流污染总体控制最为关键.

3) 相关性分析结果表明, 污染物与浊度的相关性主要受干旱天数和降雨强度影响, 干旱天数长和降雨强度大的溢流, 污染物与浊度的显著相关性(p < 0.01), 污染物主要以颗粒态存在;反之, 则污染物与浊度间的相关性不显著(p>0.05), 污染物主要以溶解态存在.

参考文献
Bersinger T, Bareille G, Pigot T, et al. 2018. Online monitoring and conditional regression tree test:Useful tools for a better understanding of combined sewer network behavior[J]. Science of the Total Environment, 625: 336-343. DOI:10.1016/j.scitotenv.2017.12.239
Bersinger T, Le Hécho I, Bareille G, et al. 2015. Assessment of erosion and sedimentation dynamic in a combined sewer network using online turbidity monitoring[J]. Water Science and Technology, 72(8): 1375-1382. DOI:10.2166/wst.2015.350
Casadio A, Maglionico M, Bolognesi A, et al. 2010. Toxicity and pollutant impact analysis in an urban river due to combined sewer overflows loads[J]. Water Science and Technology, 61(1): 207-215. DOI:10.2166/wst.2010.809
Casadio A, Maglionico M, Bolognesi A, et al. 2018. Urbanization and climate change impacts on surface water quality:Enhancing the resilience by reducing impervious surfaces[J]. Water Research, 144: 491-502.
国家环境保护总局. 2002. 水和废水监测分析方法(第4版)[M]. 北京: 中国环境科学出版社, 211-213.
Holeton C, Chambers P A, Grace L. 2011. Wastewater release and its impacts on Canadian waters[J]. Canadian Journal of Fisheries and Aquatic Sciences, 68(10): 1836-1859. DOI:10.1139/f2011-096
环境保护部.2009.中华人民共和国国家环境保护标准[S].北京: 中国环境科学出版社
刘博, 徐宗学. 2011. 基于SWAT模型的北京沙河水库流域非点源污染模拟[J]. 农业工程学报, 27(5): 52-61. DOI:10.3969/j.issn.1002-6819.2011.05.009
Lacour C, Joannis C, Gromaire M C, et al. 2009. Potential of turbidity monitoring for real time control of pollutant discharge in sewers during rainfall events[J]. Water Science and Technology, 59(8): 1471-1478. DOI:10.2166/wst.2009.169
Madoux-Humery A S, Dorner S, Sauvé S, et al. 2013. Temporal variability of combined sewer overflow contaminants:evaluation of wastewater micropollutants as tracers of fecal contamination[J]. Water Research, 47(13): 4370-4382. DOI:10.1016/j.watres.2013.04.030
潘国庆.2008.不同排水体制的污染负荷及控制措施研究[D].北京: 北京建筑工程学院
Passerat J, Ouattara N K, Mouchel J M, et al. 2011. Impact of an intense combined sewer overflow event on the microbiological water quality of the Seine River[J]. Water research, 45(2): 893-903. DOI:10.1016/j.watres.2010.09.024
Postel S L, Thompson Jr B H. 2005. Watershed protection:Capturing the benefits of nature's water supply services[J]. Natural Resources Forum, 29(2): 98-108. DOI:10.1111/j.1477-8947.2005.00119.x
Salzman J, Thompson Jr B H, Daily G C. 2001. Protecting ecosystem services:Science, economics and law[J]. Stanford Environmental Law Journal, 20: 309-332.
水利部.1991.堰槽测流规范: SL 24-1991[S].北京: 中国水利水电出版社
孙文, 王理明, 刘吉宝, 等. 2019. 北运河沙河水库沉积物营养盐分布特征及其溯源分析[J]. 环境科学学报, 39(5): 1581-1589.
Tait S J, Ashley R M, Verhoeven R, et al. 2003. Sewer sediment transport studies using an environmentally controlled annular flume[J]. Water Science and Technology, 47(4): 51-60. DOI:10.2166/wst.2003.0219
Viviano G, Salerno F, Manfredi E C, et al. 2014. Surrogate measures for providing high frequency estimates of total phosphorus concentrations in urban watersheds[J]. Water Research, 64: 265-277.
王辰, 刘全儒. 2002. 北京北沙河湿地维管植物区系的研究[J]. 北京师范大学学报(自然科学版), 38(6): 815-819. DOI:10.3321/j.issn:0476-0301.2002.06.022
邢丽贞, 王恩革, 常青霞. 2011. 合流制污水溢流治理技术综述[J]. 山东建筑大学学报, 26(2): 166-169. DOI:10.3969/j.issn.1673-7644.2011.02.017
杨雪, 车伍, 李俊奇, 等. 2008. 国内外对合流制管道溢流污染的控制与管理[J]. 中国给水排水, 24(16): 7-11. DOI:10.3321/j.issn:1000-4602.2008.16.002
张伟, 张洪, 单保庆. 2012. 北运河源头区沙河水库沉积物重金属污染特征研究[J]. 环境科学, 33(12): 4284-4290.
Zhang P, Cai Y, Wang J. 2018. A simulation-based real-time control system for reducing urban runoff pollution through a stormwater storage tank[J]. Journal of cleaner production, 183: 641-652. DOI:10.1016/j.jclepro.2018.02.130
中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会国家气象中心.2012.降水量等级: GB/T 28592-2012[S].北京: 中国标准出版社
周美成, 陈俊, 董良飞, 等. 2015. 苏南老城区合流制管网的溢流污染特征[J]. 环境工程学报, 9(7): 3159-3164.
朱光旭, 郭庆军, 陈同斌, 等. 2013. 北京市南沙河沉积物重金属污染特征及风险评价[J]. 生态学杂志, 32(8): 2148-2153.