非点源污染产出估算方法研究李丹，郝振纯，薛联青河海大学水资源环境学院，江苏南京（210098）摘要：根据非点源难以在其源头处监测的特点，对非点源污染的产出估算方法进行了初步探讨，提出在已有监测资料的基础上，利用成熟的数学模拟技术反演非点源污染负荷，从而分析得到非点源在流域的分布概况，其结果既可作为非点源污染模型应用的基础资料，同时也能指导非点源污染控制。
论文以实例论证了该方法的可行性。
关键词：非点源污染，产出估算，反演1. 引言随着城市污水处理率的提高，点源污染控制技术已经比较成熟。
然而在水体污染中占比重较大的非点源污染却仍是一大难题。
至今国内外尚无成熟和标准化的控制和监测技术1。
非点源污染发生的随机性、排放的间歇性以及污染负荷的时空变化大等特点决定了其监测、控制和处理都相当困难而复杂2。
国外对非点源污染的研究真正起步于60年代。
首先在美、英、日等一些发达国家率先开展，随后在世界各地逐渐受到重视。
40多年来，西方国家在这方面的研究工作从概念、研究方法到新技术应用，管理手段逐步发展，已经相对完善。
相比而言，我国对非点源污染的研究起步就晚得多。
北京城市径流污染及20世纪80年代初全国湖泊、水库富营养化调查和河流水质规划的研究标志着我国开始了真正意义上的非点源污染研究。
此后，我国的非点源污染研究开始变得活跃。
刘枫3于1988年对流域非点源污染量化识别进行了初步研究，并成功应用于于桥水库流域；李怀恩等4-5以非点源污染物的迁移转化机理为基础，从宏观角度和较大尺度上直接研究野外实际流域的非点源污染发生过程与特点，建立了次暴雨非点源污染负荷数学模型，并在相关流域做了验证；朱萱等6（1985）通过研究农田暴雨径流污染特征及污染物输出规律提出了采用统计技术的区域径流－污染负荷模型。
到目前为止，虽然在非点源污染定量化研究方面已经做了大量工作并有了一定进展，但是真正能在实践中推广应用的成果并不多见。
大部分是引进并改造国外已有的非点源污染模型7-10，然后这些模型对资料的要求太高，且需要率定大量的参数，根据我国目前的现状，非点源污染试验与监测条件都达不到国外的先进水平，几乎没有系统的长系列非点源污染监测资料，从而限制了这些模型的应用并造成研究上的困难及计算结果的不精确、不可靠。
鉴于此，本文尝试利用已有的水质监测资料和成熟的水质数学模型来反演非点源污染负荷，通过分析确定非点源污染在流域的分布概况，以此作为非点源污染模型运行的基础资料。
这不仅对模型的实际应用具有指导作用，同时也有利于区域面源污染的控制和管理，对少资料或无资料区域非点源污染模拟和预测具有重要意义。
2. 思路与方法2.1思路传统上按污染的发生类型将水环境污染源分为点源（point source）和非点源（non-point source）。
点源通常具有固定的排放位置，集中排放，如工业废水、污水处理厂的出水以及生活废水等。
非点源则是向环境中排放不连续且不能由一般的污水处理方法获得水质改善的污染源。
也可简单定义为晴天积累、雨天排放的不固定排放源11。
以往的定量研究都是直接利用非点源污染模型预测进入水体的非点源污染负荷过程，然后再选用合适的水质模型来评价其影响。
这种方法需要具有连续多年的流量和浓度同步监测资料，然而在我国，由于受试验和监测条件的限制，完整的野外小区实验资料很少，仅在一些流域的出口断面处测有水量水质同步资料，因而严重制约了模型的应用，并影响了模型的模拟预测精度。
针对这种情况，我们设想在已有监测资料的基础上，利用点源反演面源。
首先选取成熟的水质数学模型模拟点源污染在流域的分布，然后通过剔除点源污染负荷来反求非点源污染负荷，分析确定非点源污染在流域的分布概况，从而估算非点源污染对流域的环境污染贡献率。
如果将流域各控制段的总污染负荷量（W T ）表示为：T p n W W W =+p Q =p pe e Q =) （1）式中：和 分别为点源和非点源污染负荷（t/a ）。
p W np W 其中有W C ，W C ，式中为流域控制段的平均流量（m T e e e Q 3/s ）；和分别污染物平均浓度以及点源污染物平均浓度（mg/L ）。
e C pe C 则由（1）式可得到控制区域非点源污染负荷量为：(np e pe e W C C Q =− （2）式中各项符号意义同前。
2.2方法2.2.1 点源污染负荷的确定降雨径流的冲刷是产生非点源污染的原动力，降雨径流又是非点源污染物的载体。
如果没有地表径流的产生，非点源污染物就很难进入受纳水体。
因此可以认为非点源污染主要是由汛期地表径流引起的，而枯水季节的污染主要来自点源，此时非点源污染所占比重很小，可以忽略不计。
污染物排入水体后，在自净作用下其浓度会降低。
为计算简便，可用一维对流扩散方程表示： C ux∂=−∂KC （3） 式中：C —污染物平均浓度（mg/L ）；u —平均流速（m/s ）；K —污染物综合自净系数（1/d ）；x —距离（km ）。
上式的解为：()0exp x C t C Ku ⎛⎞=−⎜⎝⎠⎟，其中为初始背景浓度。
0C 利用上述一维模型模拟得到水体中点源污染物浓度变化过程，结合流量过程即可算出各控制段点源污染负荷：（4）,,p i p i i W C Q =式中：—第段点源污染负荷（t/a ）；—第i 段点源污染物平均浓度（mg/L ）；—第i 段平均流量（m ,p i W i ,p i C i Q 3/s ）。
2.2.2 非点源污染负荷的估算汛期的水体污染是点源与非点源污染综合作用的结果，根据此时已有的水质监测资料，由计算得到流域各控制段的总污染负荷，通过扣除前面计算得到的点源污染负荷，由（2）式即可求出各相应控制段的非点源污染负荷量。
T e W C Q =e 3. 实例应用3.1流域概况子牙河流域东临渤海，西靠云中山，包括河北南部和山西东部地区。
子牙河系主要由滹沱河（685km ）、釜阳河（410km ）、釜阳新河（132km ）、子牙河（177km ）和子牙新河（143km ）等组成，面积46868km 2。
滏阳河是子牙河两大支流之一，发源于河北省邯郸市峰峰矿区和村附近，流经磁县、邯郸、永年、曲周、鸡泽后出境入邢台，在磁县境内建有大Ⅱ型东武仕水库。
1956-1989年平均降水量为558.8mm ，其中汛期6-9月的降水量占74％～76％，降水年内分配很不均匀。
该流域易在7、8月发生大到暴雨，带来突发性水体污染事件，这与非点源污染密切相关。
点源污染则主要来自于沿河生活污水和工业废水的排放。
为了研究上的方便，同时也受资料限制，本文做如下假设：首先，只研究流域干流的污染物混合输移过程，支流均作为点源输入，采用一维稳态水质模型模拟点源的沿程变化。
其次，选取具有代表性且监测资料系列相对完整的氨氮作为模拟因子。
根据已有的水质监测资料，在研究河段上分别选取四个控制断面张庄桥、联纺桥、苏里桥和莲花口。
其中支流作为点源排入。
考虑到点源排污主要集中在邯郸县内的张庄桥－苏里段，而苏里－莲花口段则无点源排污，可以利用计算结果进行点源与非点源污染影响分析3.2设计时段流量及点源污染负荷确定据1956-1989年共34年流量资料，经水文统计计算得到50％、75％和95％保证率下的年径流量分别19370万m 3、9285万m 3、2057万m 3。
利用（3）式对枯水年各河段的点源污染变化过程进行模拟，率定得到氨氮综合自净系数K 为2.02。
然后用验证过的模型预测2003年丰水年各控制断面点源污染物浓度变化过程，并由（4）式计算得到各河段点源污染负荷，结果见表1。
表1 控制断面丰水期点源污染物浓度及污染负荷量Table 1 Point source concentration and pollution load of each section排污断面 点源污染物浓度(mg/L) 流量（m 3/s ）污染负荷（t/a ） 1 5.05 8.731605.31 2 4.32 10.081435.92 3 1.5710.54 521.853.3非点源污染负荷计算根据2003年各断面水质监测资料，计算出各河段总污染负荷，结合以上点源污染负荷值由（2）式计算得到非点源污染负荷量，具体结果见表2，各段污染负荷比较示意如图1所示。
表2 各控制段非点源污染负荷量计算结果（单位：t/a）Table 2 Non-point source pollution load of each section河段点源负荷总负荷非点源负荷非点源比例（％）2867.35 1262.04 441 1605.314705.46 3269.54 692 1435.923 521.85 3722.76 3200.91 867732.49 68∑3563.08 11295.57图1 控制河段污染负荷比较示意图Fig.1 Pollution load of each section and the total comparison4. 结果分析由表2可知，2003年研究河段总污染负荷量为11295.57t/a，其中非点源污染负荷为7732.49t/a，占总污染负荷的68％，各控制段非点源污染平均水平已超过50％。
根据2003年实测资料，该研究区段水体越往下游水质越差，可以发现这与我们只考虑点源排污时的情况不太符合。
经过分析得出，尽管沿河企业和生活污水排量不大，但河流水体水质仍然得不到改善，其主要原因就是汛期有大量的非点源污染排放进入受纳水体，造成水体水质的进一步恶化。
经实地调查，邯郸县内水土流失类型以面蚀为主，其中轻度流失面积为116.69 km2，中度流失面积为38.34 km2。
耕地面积532740亩，亩均氮肥和农药使用量分别为20.90kg和0.5kg，可能产生的面源污染负荷约为8978t/a，约占总污染负荷的74％。
与本文的计算结果比较接近。
因此，要从根本上改善釜阳河水质，除了对点源污染进行控制削减外，还要对该流域可能产生的非点源污染采取一定的防治措施。
5. 结语由于非点源污染监测困难，资料难以获取，限制了非点源污染模型的应用，因此本文提出了根据有限监测资料反演流域非点源污染负荷量的简便而有效的方法。
该方法避免了直接利用非点源污染模型要率定大量参数的困难。
对资料的要求不高，仅利用已有的点源排污和水质监测资料分析得到非点源污染在流域的分布概况，对非点源污染的控制具有十分重要的指导意义。
同时可以将估算得到的非点源污染负荷量作为运用非点源污染模型的基础资料，有利于少资料或无资料地区的非点源污染的模型及预测。
本方法在典型流域的应用还有待进一步的研究。
参考文献[1] 张维理，等. 中国农业面源污染形势估计及控制对策. 中国农业科学，2004，37（7）：1008-1017[2] 王晓燕，等. 非点源污染及其管理. 北京：海洋出版社，2003[3] 刘枫. 流域非点源污染的量化识别及其在于桥水库流域的应用. 地理学报1988，43（4）[4] 李怀恩，沈晋. 非点源污染数学模型. 西安：西北工业大学出版社，1996[5] 李怀恩. 估算非点源污染负荷的平均浓度法及其应用. 环境科学学报，2000，20（4）：397-400[6] 朱萱，等. 农田径流非点源污染特征及负荷定量化方法探讨. 环境科学，1985，6（5）：6-11[7] Knisel W G. A field-scale model for chemicals runoff and erosion from agricultural management system. Proc 13th Conf Modeling and Simulation, 1982, 4[8] Young R A. et al. AGNPS: A non-point source pollution modeling for evaluation agricultural watersheds. J of Soil and Water Conservation, 1989[9] Lane L J. et al. The water erosion prediction Project: Model overview. In Proc Natl Water Conf ASCE, 1989[10] Bwasley D B. Applying distributed parameter modeling techniques to watershed hydrology and nonpoint source pollution. Proc 13th Conf Modeling and Simulation, 1982, 4[11] 金相灿. 湖泊富营养化控制和管理技术. 北京：化学工业出版社，2001.Evaluation method of non-point source pollutioncontributionLi Dan, Hao Zhenchun, Xue LianqingHohai University, College of Water Resources & Environment, Nanjing Jiangsu （210098）AbstractAs it’s difficult to monitor non-point sources directly, some research has been done on the evaluation method of non-point source pollution contribution. It is proposed to estimate the non-point source pollution load by inversion method on the basis of existing data and perfect modeling technology, and to know the distribution of it consequently, which is either the basic data for the non-point source pollution modeling or the guideline for its control. The feasibility of this method is illustrated by an example in this paper.Keywords: non-point source pollution，contribution estimate， inversion method。