基于生态系统净化与人为调控耦合作用的水环境承载力评估

2020-08-18严长安李金城

生态学报 2020年14期

关键词：滇池人为承载力

高伟,陈岩,严长安,李金城

1 云南大学生态学与环境学院,云南省高原山地生态与退化环境修复重点实验室,昆明 650091

2 生态环境部环境规划院,北京 100012

3 昆明市环境科学研究院,昆明 650032

4 北京大学环境科学与工程学院,北京 100871

随着经济社会的发展,水资源短缺、水污染、水生态破坏等水环境问题日益突出,已成为制约人类社会可持续发展的瓶颈。我国最严格水管理制度的提出从水资源开发、用水效率和水功能区纳污三个角度设定管理目标,推动经济社会与水环境承载能力协调发展[1-2]。水环境承载力(Water Environment Carrying Capacity,WECC)作为衡量区域经济社会与水环境协调发展程度,判定区域可持续发展状态的重要指标,对区域水资源管理、水环境优化、水生态恢复具有指导意义。

目前,水环境承载力评价尚未形成统一完整的评价体系,常用的评价方法有多目标规划模型、模糊综合评价模型、结构方程模型、压力-状态-响应模型等,由于研究出发的角度不同,其构建的评价指标体系也存在诸多差异。赵卫等人以水环境系统所能承载的最大经济规模为目标,采用多目标优化模型评价了辽河流域水环境承载力,运用情景分析法探讨水环境承载力优化途径,其在排污结构优化的过程中未将生态系统净化功能与人为调控污染削减功能区分[3]。蓝希等人以水资源环境、水污染控制和社会经济承载3个子系统构建城市水环境承载力指标评价体系,采用熵权法和均方差决策法确定指标权重,评价武汉市城市水环境承载力,其所构建的评价指标体系的23个指标中包含了城镇污水处理等人为调控指标,但缺乏湿地等生态系统污染净化相关指标[4]。崔丹等人以水环境承载力和水环境压力准则层确定水环境承载力综合评价指标,采用结构方程模型确定各指标权重,评估了湟水流域小峡桥断面以上的水环境承载力,其评估指标中包含水质净化能力相关指标,但缺乏人为调控污染削减相关指标[5]。徐志清等人从社会经济、水资源和水生态3个子系统中选取18个指标构建水环境承载力评价指标体系并采用熵权法确定指标权重,利用模糊综合评价模型对南京市水环境承载力进行评价,其指标体系中涉及了污水处理率、绿化覆盖率等人为调控指标,缺乏生态净化相关指标[6]。众多专家学者在水环境承载力计算过程中未能突显生态系净化作用与人为调控作用在流域环境承载力中的重要性。截止2015年,我国共建设污水处理厂5343座,污水处理量达到507.78亿m3/a,污水处理、再生水利用、湿地工程等环保设施不断完善,已成为改善流域水环境质量,提高水环境承载力最为重要的手段[7-8]。第二次全国湿地调查结果显示,全国湿地总面积 5360.26万hm2,湿地率为5.58%,人工湿地生态系统对污水有十分重要的净化作用,已成为改善我国地表水环境措施中重要的一环[9- 11]。如何在水环境承载力评估方法体系中凸显生态系统和人为调控双重作用显得十分关键。

针对当前水环境承载力评估方法体系中未突显生态净化作用与人为调控作用重要性的问题,本论文提出基于生态系统净化—人为调控耦合作用的流域水环境承载力评估模型,并以滇池流域为案例做简要分析,以期为流域水环境承载力评价提供新的思路与方法。

1 材料与方法

1.1 研究区域概况

滇池流域(24°29′—25°28′N,102°29′—103°01′E)地处长江、红河、珠江三大水系分水岭地带,地势北高南低,平均海拔1900 m,面积2920 km2。流域所处区域为亚热带湿润季风气候区,多年平均降水量978.8 mm(1951—2015年,昆明气象站)。流域全部位于昆明市境内,2015年常住人口为407.34万人,其中城镇369.68万人,城镇化率高达90.76%,流域国内生产总值达到3003.13亿元,占昆明市的75.65%。

2015年滇池流域土地利用中面积最大的是林地,为992.09 km2,占总面积的34.43%,其次是建筑用地,耕地,草地和水域,分别占总面积的19.72%,18.66%,15.52%和11.67%。滇池多年平均水资源量5.30×108m3,2015年人均水资源量仅244 m3,远低于国际公认的人均500 m3的“极度缺水”标准[4],因此,滇池流域水质受到了自然条件和经济社会的双重压力。

1.2 流域水环境承载力模型

从流域生态系统净化和人为调控两个角度出发,综合物质流分析和生态经济评估模,提出流域水环境承载力评估模型,公式如下：

WECC=ESV+HAV

(1)

式中,WECC为流域环境承载力(t/a)；ESV(Watershed Ecosystem Pollution Purification Capacity)为流域生态系统污染净化能力(t/a)；HAV(Watershed Human Pollution Reduction Capacity)为流域人为调控污染削减能力(t/a)。

本研究采用水环境承载力综合指数作为评价流域是否超载的定量指标,公式如下：

(2)

式中,WECSI(Comprehensive Index of Water Environment Carrying Capacity)为水环境承载力综合指数；WPD(Watershed Pollutant Discharge)为污染物排放量(t/a)；WECC为水环境承载力(t/a)；下标i、j分别为污染源类型和污染物类型。

1.3 流域生态系统污染净化能力评估方法

污染物在生态环境中降解是一个十分复杂的过程,本研究采用物料平衡法评估流域生态系统污染净化能力。物料平衡法又称物料衡算法,根据质量守恒定律产生,原理是在生产过程中,物质的输入和输出保持不变,其广泛应用于河流上下游监测断面间污染物变化,准确得到区域内污染物的排放和降解情况[12]。物料平衡法在国外以成功应用于Calumet Lake[13]、Black River[14]等流域,国内成功应用于三峡水库氮迁移转化[15]、黄海近岸多环芳烃源解析[16]等。物料平衡法原理如图1所示,流域污染物径流输出量与降解量之和等于流域污染物排放量与径流输入量之和[17-18]。

图1 基于物料平衡法生态系统净化能力原理示意

根据物料平衡原理,流域多年污染物蓄变量为零,即污染物输入与输出相等。因此,生态系统污染净化能力计算公式如下：

(3)

式中,ESV(Degradation of Pollutants in Ecosystems)为生态系统污染物降解量(t/a)；WPD为污染物排放量(t/a)；UPI(Pollutant Input of Upstream Runoff)为上游径流污染物输入量(t/a)；WTI(Input Pollutant from Outer Watershed Water Transfer)为外流域调水污染物输入量(t/a)；HAV为人为调控污染物削减量(t/a)；OPA(Output Pollutant of Runoff)为污染物径流输出量(t/a)。

1.3.1污染物排放量计算方法

污染物的排放量可分为点源和非点源,本研究中点源包括：城镇生活、工业、第三产业、规模化畜禽养殖,非点源包括：农村生活、散养畜禽、农业种植。其中,城镇生活、第三产业、农村生活污染物排放量计算参照《生活源产排污系数及使用说明》[19],散养畜禽污染物排放量计算参照《中国有机肥料养分志》[20],农业种植污染物排放采用标准农田修正法,参照《全国水环境容量核定技术指南》[21],规模化畜禽养殖点源污染物排放量可从污染源普查数据中获得。

1.3.2污染物输出量计算方法

本研究采用 The Load Estimator Model (LOADEST模型)估算河流出境断面的污染物通量。该模型基于实测污染物浓度、流量数据和监测日期三个变量,建立污染物浓度与其他两个变量的回归关系,以此为基础,计算出境断面污染物通量[22- 24]。当水质监测数据信息不完整,出现缺失数据时,采用修正最大似然估算法消除这种误差；当模型残差符合正态分布时,使用最大似然法和修正最大似然法；当模型残差不符合正态分布时,采用最小绝对偏差法。

1.3.3流域人为调控污染削减能力计算方法

本研究从减排角度,提出流域人为调控能力等于流域污染物削减能力,即流域内污染控制工程所发挥的污染物削减功能,计算公式如下：

HAVj=HAVPj+HAVNPj

(4)

式中,HAV为流域人类调控污染削减能力(t/a)；HAVP为点源污染物削减能力(t/a)；HAVNP为非点源污染物削减能力(t/a)；下标j为污染物类型。

点源污染削减能力计算包括污水处理厂类削减量和非污水类排放源的削减量,计算公式如下：

(5)

式中,WT为污水处理厂等点源污染削减工程的污水处理量(万m3/a)；IC和OC分别为污水处理设施的进口和出口污染物浓度(mg/L)；PT为规模化畜禽养殖厂的设施污染物削减量(t/a)。

畜禽养殖业的污染排放方式以固体(粪便)和液体(粪尿)为主,削减能力的计算方式也有所差别。本研究基于《畜禽养殖业污染治理工程技术规范(HJ497—2009)》和《畜禽养殖业污染物排放标准(GB 18596—2001)》中提供的处理方式和参数提出规模化畜禽养殖的污染削减计算公式,公式如下：

(6)

式中,LVS为畜禽养殖规模(头或只或羽)；LSR和LLR分别为粪便和尿液中的污染物排放系数(kg 只-1a-1)；LVSS为畜禽存栏规模；LPW为污水排放系数(L 只-1d-1)；LWC为畜禽养殖场排放污水中的污染物浓度(mg/L)；r为畜禽类型。

非点源主要包括散养畜禽、农村生活、农业种植等。非点源污染削减计算相对困难,本研究主要计算可量化评估的工程措施。非点源工程措施分为污水处理设施类、湿地工程类和其他类3个部分,计算公式如下：

HAVNPj=NPWWj+NPWTj+NPOTj

(7)

式中,NPWW为污水处理类工程的污染物削减量(t/a)；NPWT为湿地类工程的污染物削减量(t/a)；NPOT为其他工程的污染物削减量(t/a)。

污水处理设施类工程是指农村生活污水处理设备,其计算方法与城镇污水处理厂类似,计算公式如下：

(8)

式中,RWT为农村生活污水处理厂等非点源污染削减工程的污水处理量(万m3/a)；RIC和ROC分别为污水处理设施的进口和出口年均污染物浓度(mg/L)。

湿地类工程是指河口或湖岸的湿地,其削减量的计算方法与污水处理设施计算方法相同。

其他污染源是指散养畜禽。散养畜禽的污染处理尚无相关规范和标准,农村地区也难以建设集中式污染处理设备,本研究主要考虑散养畜禽粪便的综合利用,如沼气和堆肥。

(9)

式中,LVS为畜禽养殖规模(头或只或羽)；LSR和LLR分别为粪便和尿液中的污染物排放系数(kg 只-1a-1)；MT为粪便处理量(t/a)；LSE为畜禽粪污系数(kg 只-1a-1)；RT为畜禽养殖周期(d)；k为研究区的畜禽粪便处理设施编号。

2 案例研究

2.1 滇池流域生态系统污染净化能力评估

2.1.1污染物排放量

根据上文中污染物排放量计算方法,对滇池流域的主要污染源排放量进行计算,具体结果如表1所示。滇池流域的化学需氧量(Chemical Oxygen Demand, COD)、氨氮(Ammoniac Nitrogen, NH3-N)、总磷(Total Phosphorus, TP)排放量分别达到120740、14228、2342 t/a。污染物排放以点源排放为主,COD、NH3-N、TP排放量分别占总排放量的95%、94%、82%,而点源排放中以城镇生活排放为主,分别占总排放量的80.46%、90.18%、77.07%。表1中未对滇池流域规模化畜禽养殖污染物排放进行统计,原因是昆明市颁布《在“一湖两江”流域禁止畜禽养殖的规定》[25],自2008年6月30日起,禁止任何单位和个人在滇池流域2920km2范围内进行大规模畜禽养殖。

表1 滇池流域主要污染源污染物排放量/(t/a)

2.1.2流域出境污染物量

滇池流域出口有2个,海口河与西园隧洞,基于2000—2015年逐日流量和逐月水质监测数据,采用LOADEST模型,估算出流域出境逐日输出负荷,并统计多年平均污染物输出负荷量,最后将2个出口的负荷进行加和处理,作为流域出境污染物通量值。经计算,滇池流域COD、NH3-N、TP出境污染物多年平均值分别为17693.43、899.42、401.63 t/a。其中海口河COD逐日输出负荷多年平均值为9137.06 kg/d,变化范围为501—68503 kg/d,NH3-N逐日输出负荷的多年平均值为574.03 kg/d,变化范围为14—1451 kg/d,TP逐日输出负荷的多年平均值为141.91 kg/d,变化范围为10—1753 kg/d。具体结果如图2所示,海口河逐日污染物负荷量年际和年内波动较大,主要与径流量有关,其中2010—2013年滇池流域出现较严重的干旱,导致该时期的污染物输出量较低。

西园站COD逐日输出负荷多年平均值为8556.49 kg/d,变化范围为0—47049 kg/d,NH3-N逐日输出负荷的多年平均值为325.39 kg/d,变化范围为0—749 kg/d,TP逐日输出负荷的多年平均值为259.72 kg/d,变化范围为0—862 kg/d。具体结果如图2所示,西园隧洞逐日污染物负荷变化规律与海口河相似,年际和年内变化较大,与海口河不同的是西园隧洞出流量与负荷在干旱期间未出现显著下,主要原因是该区域污水来自于城镇生活。

图2 海口站、西园站逐日流量及负荷

2.1.3外流域调水污染物输入量

滇池流域于2007年启动掌鸠河供水工程,年均调水量1.75亿m3,2012年启动清水海供水工程,年均调水量0.68亿m3。2013年启动牛栏江调水工程,2015年向滇池流域补水5.27亿m3,滇池2011—2015年调水量如图3中所示,2014年后牛栏江调水工程已超过总调水量的50%。滇池流域2011—2015年调水水质如图3中所示,调水水质年际变化大,COD、NH3-N含量整体呈下降趋势。外流域调水的一方面缓解了流域内水资源短缺和水污染问题,另一方面也带来污染负荷。经计算,滇池流域2015年COD、NH3-N、TP污染物输入负荷分别为1039.76、524.29、167.33 t/a。

图3 滇池流域2011—2015年调水量及水质

根据1.3节公式(3),计算得到滇池流域2015年COD、NH3-N、TP三种污染物生态系统净化能力分别为17935.21、6187.58、700.37 t/a。

2.2 滇池流域人为调控污染削减能力评估

2015年底,滇池流域已建成23座污水处理厂,污水处理能力达到199万m3/d。未纳入主城污水处理厂纳污范围的集镇建设污水处理站10个,处理能力为1.23万m3/d,开发区、工业园区建成5个园区污水处理厂,处理能力为9.5万m3/d。全流域已建成885个村庄生活污水收集处理设施,484座再生水利用实施。根据流域人为调控能力计算方法,计算得到2015年滇池流域的COD、NH3-N和TP的削减量分别为86151、7665、1407 t/a。各子流域污染物削减能力分布如图4所示,整体呈现北高南低,其中COD削减能力大于NH3-N大于TP,因昆明市主城区分布在滇池北岸,污水处理设施完善,污染物削减量高。南部地区为昆明市农田耕种区,污水收集处理率低,污染物削减量少。

图4 滇池流域子流域污染物削减量

2.3 滇池流域水环境承载力评估

根据水环境承载力计算模型,滇池流域环境承载力由滇池流域生态系统污染净化能力和滇池流域人为调控污染削减能力。滇池流域2015年COD、NH3-N、TP环境承载力分别为104086.21、13852.58、2107.37t/a,其中生态系统污染净化能力分别为17935.21、6187.58、700.37t/a,人为调控削减能力分别为86151、7665、1407 t/a,分别占总承载力的82.77%、55.33%、66.77%,人为调控污染削减能力已经超过生态系统污染削减能力,成为流域环境管理重要的手段。滇池流域各子流域环境承载力如图5所示,3种污染物承载能力分布上均呈北高南低,与人为调控污染削减能力分布类似。北部地区以森林、草地为主,产流能力高,南部耕地面积多,产流能力低,相应的生态系统污染净化能力北部比南部高。

图5 滇池流域子流域环境承载力

滇池流域COD、NH3-N和TP水环境承载指数分别为1.16、1.09和1.11,均大于1,处于超载状态。根据流域水环境承载力综合指数为所有指标中的最大值,滇池流域2015年水环境承载力综合指数为1.16,超载16%。

滇池流域15个主要入湖河流子流域中有6个子流域水环境承载力综合指数大于1,占流域陆域面积的54%,具体分布如图6所示,其中最大值是草海子流域,为1.98,最小值是海河子流域,为0.22,水环境承载力综合指数空间分布上呈现北部高,中部低的特征,这与区域经济发展模式和污染源分布有关。从污染物排放格局中发现,排放强度以城区为中心呈放射状分布。各子流域水环境承载力制约因子分布如图6所示,整体上COD制约因子的流域面积多于TP,其中北部子流域以COD制约因子为主,南部子流域以TP制约因子为主,主要原因是为北部为昆明市主城区,城镇生活污染压力大,南部为昆明市耕地集中区,农业面源污染严重,且南部为富磷区,因此TP为制约因子。

图6 滇池流域子流域水环境承载力综合指数与制约因子

本研究构建的流域水环境承载力模型计算得到的水环境承载力会高于目前常用模型,如石建屏等人从水资源承载力和水质量承载力两个角度采用层次分析法和指标体系法评估了滇池流域2003—2010年水环境承载力,其计算的水环境承载力低于本研究[26]。邓伟明等人以水资源和水环境承载力作为广义水环境承载力,采用水环境承载率评估滇池流域2008年水环境承载状态,其计算的水环境承载力低于本研究[27]。造成差异的主要原因是本研究水环境承载力中包含人为调控污染削减能力,而目前常用的水环境承载力评估模型在计算污染物排放量时将削减量直接扣除,即压力中包含了污染削减能力,导致水环境承载力计算结果低于本研究。2015年滇池流域开展牛栏江调水工程以改善滇池水污染现状,调水水质COD、NH3-N、TP含量年平均值明显低于滇池出流水质,调入清洁水资源增大了滇池流域水环境承载力。根据调水水量水质与流域出流水量水质,计算得到调水出境携带COD、NH3-N、TP分别为16654、375、235t,相当于水环境承载力的16.00%、2.71%、11.13%,调水工程对滇池流域污染状况的改善具有一定的推动作用。