长江水质污染-数学建模

1、2005高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的 ，如果引用别人的成果或其他公开的 资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规 则的行为，我们将受到严肃处理。我们参赛的题目是： 长江水质的评价和预测我们的参赛报名号为（如果赛区设置报名号的话）：所属学校

2、（请填写完整的全名）：江苏大学参赛队员（打印并签名）：1. 滕成龙2. 周振华3. 黄旭升指导教师或指导教师组负责人（打印并签名）：日期：2005年_9_月19日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：长江水质的评价和预测摘要本文分别对长江水质建立了评价预测模型和污水治理模型，提出了合理化建议。针对问题一，应用标准指数评价法对单项水质进行参数评价，在此基础上建立综合 评价指数，给出了多项水质的综合评价检验，并利用统计数据合理地分析了各地区的水 质污染

3、状况。评价结果显示长江综合水质等级为三类，并通过建立退出标准系数，求得 了主要污染因素。该模型既对水质进行了综合评价，又指出主要污染因素，全面合理， 重点突出。针对问题二，利用质量守恒定律，对长江干流观测点水域进行物料衡算，并建立假 想水段模型。根据自身污染物排放量找出主要污染物高锰酸盐指数和氨氮指数的污染源 所在地区，结论是主要污染源为湖南、湖北、江西（排名有先后）三省地区，同时，江 苏省的高锰酸盐污染也较为显著。针对问题三，通过建立废水排放量随年度变化的曲线拟合函数和各类水百分比与废 水单位浓度的对数拟合模型，进而对长江未来水质污染的发展趋势做出预测分析，可知 非饮用水所占比率逐年增加，如

4、果不加以控制，到2014年非饮用水所占比率将超过57.9%，其中IV类水和V类水之和将达到30.8%,长江流域面临生存危机。针对问题四，根据预测模型，通过一维搜索算法，求得长江干流在未来 10年内， 每年都满足“ IV类和V类水的比例控制在20%以内，且没有劣V类水”的要求下的年 污水处理量如下表：表1-1各年份的年处理废水量（单位：亿吨）年份2005200620072008200920102011201220132014处理量122145169196224253284317352388针对问题五，应用扩展费用-效益分析法论述了长江污染的利益驱动因素，以及将 造成的社会和生态不良后果，建立改进

5、的治污费用最小化优化模型。为解决长江水质污 染问题，提出充分利用现有设备重点治污、循环利用和动态监测等合理化建议。在问题中运用相似度判断法、反推法对模型进行了检验，证明以上模型和计算结果 的合理性、准确性和可行性。1一 问题重述长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重。相关政府部门 和专家们高度重视， 由全国政协与中国发展研究院联合组成 “保护长江万里行” 考察团， 从长江上游宜宾到下游上海，对沿线 21 个重点城市做了实地考察，揭示了一幅长江严 重污染的真实画面。为了对长江水质进行科学合理的评价预测，题目附件给出了长江沿 线 17 个观测站（地区）近两年多主要水质指标的检测

6、数据，以及干流上 7 个观测站近 一年多的基本数据（站点距离、 水流量和水流速）和“ 19952004 年长江流域水质报告” 的主要统计数据。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和 上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环 境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净 化能力的指标称为降解系数。 事实上，长江干流的自然净化能力可以认为是近似均匀的， 根据检测可知 ，主要污染物高锰酸盐指数和氨氮的降解系数通常介于 0.10.5 之间，可 以考虑取 0.2 （单位： 1/天）。根据国标 （GB3838-200

7、2） 给出的地表水环境质量标准 中4个主要项目标准限值，可将地表水分为六类，其中I、U、川类为可饮用水，IV、V类和劣V类水为非饮用水。我们要对下列问题进行研究分析：（ 1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染 状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要 在哪些地区。（3）假如不采取更有效的治理措施，依照过去 10 年的主要统计数据，对长江未来 水质污染的发展趋势做出预测分析，研究未来 10 年的情况。（4）根据上述的预测分析，计算出如果未来10年内每年都要求长江干流的IV类和 V类水的比例控制在20%以内，且没有劣V类水,每年

8、需要处理的污水量。（5）对解决长江水质污染问题提出切实可行的建议和意见。二 问题假设1 长江水质平稳变化，检测数据真实可靠； 2认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水； 3认为污水集中在观测点集中排放，不考虑沿途排污问题；4影响地表水质量的主要项目个数取 4,分别为溶解氧（DO）、高锰酸盐指数（CODMn）、 氨氮（NH3-N和PH值（无量纲）；5 溶解氧（DO）为7.5mg/L时，认为其饱和率为 90%; 6主要污染物高锰酸盐指数和氨氮的降解系数取0.2（单位： 1/天） ；7污染物降解后不会造成二次污染；三 符号说明SDO , j 溶解氧在第 j 个观测点的标准

9、指数；DO f 饱和溶解氧的浓度，单位： mg/L ；DOs 溶解氧的评价标准，单位： mg/L；DO j 溶解氧在第 j 个观测点的浓度，单位： mg/L；SpH,j pH 在第 j 个观测点的标准指数；pH sd 评价标准中规定的 pH 值下限；pHsu 评价标准中规定的 pH 值上限；Si, j单项水质参数i在第j点的标准指数；G,ji污染物在第j个观测点的浓度，单位：mg/L;Cs,j i 污染物的水质评价标准，单位： mg/L；Sjj点的综合评价指数；m影响地表水质量的主要项目个数，m=4；Ai , j 近一年多 i 污染物在第 j 个观测点的自身排放总质量，单位： mgQj 近一年

10、多第 j 个观测点的水流量，单位： m3/ s；D j 近一年多第 j 个观测点的与源点(攀枝花)的距离，单位： kmU j 近一年多第 j 个观测点的水流速，单位： m/s；Q(n)第n年长江总流量，单位：108m3 ；W(n)第n年沿江经济体排放污水总质量，单位：亿吨；G(n)第n年国内生产总值(GDP)，单位：亿元；X(n)第n年长江沿岸GDP占当年全国GDP的百分数；D(n)第n年每亿吨污水所带来经济收益，单位：亿元；B!,i(n)第n年每亿吨污水处理到第i级合格水所需要的费用，单位：亿元;B2(n)第n年地下合格水的价格，单位：元/吨;T（n）第n年国家税收收入占全年 GDP的比率;

11、Wi第i个评价因子对水污染的权重;四模型建立与求解4. 1问题一的模型建立和求解要对长江水质作出定量的综合评价，先要选取合适的评价因子，建立合理的评价标准。根据本题的特点，选择国标（GB3838-2002）给出的地表水环境质量标准中 4个 主要项目溶解氧（DO）、高锰酸盐指数（CODMn）、氨氮（NH3-N和PH值（无量纲）作为 水质评价因子，构建各个评价因子的标准指数和综合评价指数，分别进行总体评价和地 区评价。4. 1.1评价模型的建立因为4个水质评价因子物理化学性质各有区别，水质等级评价的标准限值也有很大 差别，其检测数据的绝对值之间不具有可比性，所以需要在检测值与标准限值之间相对 联系

12、，建立统一的，具有相互可比性的标准指数。由文献1可以得到：高锰酸盐指数（CODMn和氨氮（NH3-N的标准指数为：Si,jCi,jCs,j5#溶解氧的标准指数为:SDO,jDO f -DO jDOf-DOs(DOj - DOs)SDO,j= 19-DOjDOs(DOj : DOs)#DOf =7.5 "0.9 ： 8.33mg/LpH的标准指数为:SpH ,j7.0 - pHj7.0 - pHsdSpH,jpH j -7.0pHsu 一 7.0(pH j "0)(pHj 7.0)参照地表水环境质量标准（GB3838-2002）中的标准限值，对近两年来各检测点 检测值取算术平

13、均作为长江对各评价因子的总体检测值，分别求取各因子的标准指数。 标准指数越大，说明检测值与假定水级的标准限值相差越大，水质越坏。若某个单项水 质因子的标准指数大于 1则表明该水质参数超过了假定水质等级标准限值，不满足该 类等级。由各指数的计算公式可以看出，单项水质参数评价能清楚地指示出单个污染因 子对水质变化的影响程度，判断主要污染因子的污染时段和水域范围，有助于直接了解 水质状况与评价标准之间的关系。各评价因子对水质的影响是不均衡的， 这种差别可以在某个评价因子引发水质重要 污染事故发生的频率上得到体现，根据频率的不同赋予各评价因子以相应的水质影响权 数，对于水质的综合评价具有重要意义。单项

14、水质评价标准指数的增加和评价因子权重 的增加都将对水质综合评价带来不利影响，综上分析，给出水质综合评价指数：mmSj =z WjSjj送 Wj =1i 壬i =1近两年多检测结果显示长江发生重度污染 73次，因溶解氧（DO）造成17次，高锰酸盐指 数（CODMn造成12次，氨氮（NH3-N造成43次，pH造成1次，所以可得：Wnh =0.58，Wco=0.16，Wdo =0.23，WpH =0.03用S,j确定水的等级，根据Sj的值与1的接近程度确定水质的检测结果与某一等级水的 接近程度，Sj和Si,j共同作为长江水质综合评价的量化指数，对长江水质做出定量的综 合评价。单项水质参数评价，能清楚

15、地指示出单个污染因子对水质变化的影响程度，判 断主要污染因子的污染时段和水域范围，有助于直接了解水质状况与评价标准之间的关 系。二者互相协调，能对水质做出全面合理的评价。4. 1. 2评价模型的求解求解原则：先假定水质符某一类水相应标准，进行单项水质参数评价，当各评价因 子标准指数均小于1时，认为假设成立，否则向下退一个等级（到 V类结束），直到所 有指数均小于1为止。如果V类水标准都不符合，就认为是劣 V类水。最后求取水质综 合评价指数Sj，结合单项水质评价指数综合评价水质等级。由此编制程序，求解可得下表推证过程:表4-1假设PH值标准指数溶解氧标准指数高锰酸盐标准指数氨氮标准指数I0.37

16、481.8029 >11.4897 >14.1534 >1类Sj =3.0732,综合指数不符合I类水标准II 类0.37480.79690.74481.2460 >1Sj =1.0364,综合指数也不符合II类水标准III 类0.37480.52970.49660.6230Sj =0.5739,各标准指数小于1,故综合指数符合III类水标准，为判断主要污染因素，引入“退出标准指数”概念：退出标准指数：推证过程中，最大标准指数刚达到小于1的目标时，全部标准指数的值。由此值可以清楚判断主要污染因素，标准指数最接近1的即为最主要污染因素。由上表可知，长江近两年的综合水质情况

17、符合III类地表水标准，且最主要的污染 指数为氨氮指数（因为其为最接近1的退出标准指数）。当假定水质为II类时，求得相应的Sj值为1.0364，接近1，故可认为长江综合水质 虽然为III类，尚偏向II类水。但是如果任由污染进行下去将会使长江的全流域综合水 质完全进入III类，甚至预测等级变为不可饮用水，这是十分危险的。因为长江流域极 广，全流域综合值如此，局部地区的污染状况将更惊人。长江治理，刻不容缓。下面分析各地区水质的污染状况，仍然采用标准指数判别法，求解判断过程参照整 个长江流域的综合评价过程。通过MATLABS程并求解可得：表4-2各地区水质污染等级及退出标准指数表（平均值）退出标准指

18、数PH值标准 指数溶解氧标 准指数高锰酸盐 标准指数氨氮标准 指数年水质等级四川攀枝花龙洞H0.3000P 0.65710.05000.2000II重庆朱沱0.31500.03290.70000.6800II湖北宜昌南津关0.03500.15700.96670.5500III湖南岳阳城陵矶10.2900P 0.79860.72500.6800II江西九江河西水厂0.17000.91860.42500.2600II安徽安庆皖河口0.2600P 0.76860.80000.4400II江苏南京林山0.3900P 0.61430.77500.2200II四川乐山岷江大桥0.33000.77500.5

19、8000.3533IV四川宜宾凉姜沟厂0.5050P 0.30140.60000.5000II四川泸州沱江二桥0.31500.80880.36000.7067IV湖北丹江口胡家岭0.81500.80000.45000.2000II湖南长沙新港0.2100P 0.56500.71670.9900III湖南岳阳岳阳楼0.36500.88860.35000.4200II湖北武汉宗关0.50000.81570.60000.3400II江西南昌滁槎:0.3600P 0.94600.18330.9200III江西九江蛤蟆石0.14000.62710.67500.3000II江苏扬州二江营0.14500.6

20、1430.40000.3000II各地区年综合水质等级可以很好反映总体情况，但是另一方面也会掩盖某些特定月份的水质恶化情况，比如攀枝花的年综合水质等级为II，但是波动明显（见下图）。这一般与排污企业的季节性工作有关，单纯求其全年平均值会忽略某些地点的季节性排 污。譬如造纸厂，其原材料的供给是局限在某几个月的，因而可能会造成某些地点只有 在一年中的某几个月份才会造成较大污染，如果年年如此，那么那个地点也应该算作污 染源，全年平均就会忽略掉这一因素。为避免这一情况，需对全部数据进行审阅。经审 阅结果发现，并不存在季节性排污的地点，因次本文的处理方法是合理的。图4-1攀枝花月水质等级变化条形图4.

21、2问题二的模型建立与求解一个地区的水质污染主要来自本地区的排污和上游的污水，在判断污染源的时候，必 须排除上游未被降解污水的影响，所以判断以高锰酸盐指数和氨氮为主要污染物的主要 污染源的依据是当前观测点的自身排污总量 A,j的值。4. 2. 1模型建立选择上游第二个观测点为A,j计算的起始基点，依次建立假设水段，水段中污水持续时间为上一观测点到计算观测点经过的时间 t,2® Dj）55各水段总净化系数为：（1_0.2）td （td为t将单位变换到天的四舍五入值）。通过计算观测点前的污染物总量经过了前段距离的净化。由于支流的引入，流量越 往下游是呈递增关系的，考虑地域广阔的因素，可以认

22、为其符合均匀分布，因而流入的 计算流量可由上游和下游的流量取均值得到，由本观测点数据可得到观测污染总量。二 者之差即为本地区自身污染量（即污染源量）。如下图所示：对长江干流各检测点近一年多的检测数据分类处理，进行水段物料衡算。由分析可知第j个观测点总检测污染量为103G,jtQj,j，净化后的上游污染量为lO'GgQavty -0.2）rd，自身净排污量=总检测污染量-净化后的上游污染量，所以可得 物料衡算的假想水段模型为：忙=103Ci,jtQi,j -103Ci,j.Qavt （1-0.2）rd对每月求得该观测点的自身污染总量之后，再求其全年的平均值。通过MATLABS程并求解可得

23、：表4-3污染源净排污总量（单位：吨）重庆朱沱湖北宜昌南 津关湖南岳阳城 陵矶江西九江河西 水厂安徽安庆皖 河口江苏南京林山COD Mn7998192401986014724199211141NH3-N879.71688.31823.51264.5155.928.2由此表可以看出湖南岳阳城陵矶、江西九江河西水厂、湖北宜昌南津关（排名有先后） 是主要的污染源观测点。由此可以得知，这两种污染物的污染源主要是在湖南、湖北、江西三省地区。但是江苏省的CODM污染也较为显著。对“非饮用水百分比矩阵”与“ VI,V类、劣V类和矩阵”进行相似性检验，即求 得每对元素的方差，方差值越小，相似度越高。由此可推证

24、问题三模型是一个精确模型， 误差很小。4. 3问题三的模型建立与求解为通过近十年的废水排放情况来预测未来十年的水质情况，应考虑到各量间的关 系。通过观察附表4所给各项数据可发现，除了少数奇点的影响外，各类数据都呈规律 性变化，故可采用曲线拟合模型。虽然奇点的影响会明显降低拟合精度，但是采用Matlab 内置的曲线拟合工具箱cftool可以较好地剔除无关项，或进行散点平滑处理，可以有效 的减小误差。曲线拟合有很多种类型，对于一般情况常采用多项式拟合，通过对比其各项误差系 数判断优劣。而一些有上限值的特殊量则需要采用参数拟合，譬如百分比就可采用对数 拟合模型。通常判断拟合结果是否可信的原则是：保证

25、均方差尽量小的情况下使得各项回归系 数尽量接近1,误差平方和系数尽量接近0。为了能够建立好的预测模型，需做如下几点说明：1 .由于预测是以全年平均值为基准的， 故在对未来十年状况预测时采用水文年的数据。2. 定义单位流量的污水排放总量（即污水单位浓度）M p （单位：吨）Mp皿P Q（n）3. 由于1998年是百年难遇的特大汛期，导致其水流量急剧增长，因此，该年的数据是 一个奇异点。明显影响到拟合过程，故在拟合分析中予以剔除。4. 除1998年以外，可认为年总流量近似符合均匀分布（可用 Matlab概率统计工具箱 函数检验），故采用除1998年流量以外的流量算术平均值作为预测年总流量。为了清楚

26、地反映模型，故先建立模型的建立步骤：1. 利用最近十年的废水排放总量向量，曲线拟合预测未来十年的废水排放总量。2. 求得平均的年总流量，进而求出未来十年的废水浓度。3. 利用过去十年的废水浓度向量，曲线拟合各类水的百分比。4. 由已经预测的未来十年浓度代入第三步的拟合式，以预测未来十年的各类水百分比模型的建立:步骤一：由于曲线多项式拟合包含线性和不同次数的非线性拟合，故需要对不同次数的曲线拟合进行比较，通过误差分析选择更优者。线性拟合结果：f (x) =12.86x 149对于拟合结果进行误差分析可知: 回归系数：0.9407拟合度回归系数：0.9332 均方差：10.37拟合图像如下:图4-

27、3二次非线性拟合结果：f(x) =0.8352x2 3.676x 167.4对于拟合结果进行误差分析可知：回归系数：0.966拟合度回归系数：0.9563均方差:8.39通过对两种拟合结果的误差分析可知，二次非线性拟合明显优于线性拟合。拟合图像如下:图4-4步骤二：由步骤一可得：表4-4未来十年的废水排放总量（单位：亿吨）年份2005200620072008200920102011201220132014排放量308.9331.78356.34382.56410.46440.03471.26504.17538.75575对近十年年流量剔除1998年后的矩阵求算术平均, 可得到流量预测值（即均值

28、）为 9534.9由此可求得:表4-5未来十年废水浓度值（单位：吨/立方米）年份2005200620072008200920102011201220132014浓度0.0323960.0347960.0373720.0401220.0430480.0461490.0494250.0528770.0565030.060305由于各类水的比率分布具有较大的随机性，为了获得稳定性数据，将水类指标划分为可 饮用水和非饮用水两类分析。为下文参照考虑和进行数据比较，再将非饮用水分为“IV,V ”类水和“劣V”类水两类进行拟合分析。步骤三:预测非饮用水的百分比，用废水单位浓度向量与非饮用水百分比向量进行拟合

29、分析:如果进行线性拟合，非饮用水的线性拟合结果为：由此可得：表4-6线性拟合数据(不可取)预测值2005200620072008200920102011201220132014非饮用水百分比35.639.443.547.952.657.562.768.274.080.1可饮用水百分比64.460.656.552.147.442.537.331.826.019.9由此表可以看出非饮用水的增长速率非常快，如果继续预测甚至会出现负值，显然不符 合事实。因为本质上非饮用水百分比不可能无限制增长，所以线性拟合模型是不可取的。考虑到百分比在增长到一定值之后会趋于一个极限值，因此对数据进行对数关系拟合。对数

30、拟合模型为：f(x) =aln(x) b非饮用水的对数拟合结果为：f (x) =382n(x) 165.2对拟合结果进行误差分析可得：回归系数：0.7505均方差:3.888由此可知，拟合结果是可接受的。拟合曲线如下图所示：图4-5同理，V,l类水的对数拟合结果为：f(x) =16.271 n(x) 76.52劣V类水的对数拟合结果为：f(x) =21.93ln(x) 88.66步骤四：由此拟合结果可得，问题三的模型结果为：表4-7未来十年污染趋势表(水文年，全流域)预测值2005200620072008200920102011201220132014不可饮用水百分比34.136.939.54

31、2.345.047.750.352.955.457.9可饮用水百分比65.863.160.457.755.052.349.747.144.642.1IV,V类水百分比20.721.923.024.225.326.527.628.729.830.8劣V类水百分比13.415.016.518.119.721.222.724.225.627.1很显然，对数拟合的结果是比较符合事实状况的4. 4问题四的模型建立与求解首先采用与问题三相同的分析方法可得干流的水文年污染趋势预测。(同理，需要剔除1998年项)由此可得：IV,V类水百分比对数拟合结果：f(x) =31.761 n(x) T35.3劣V类水百

32、分比对数拟合结果：f(x) =16.461 n(x) 64.66进而由程序可以得到:表4-8未处理时的预测百分比值预测值(干流)2005200620072008200920102011201220132014IV,V类水百分比26.428.630.933.235.437.639.841.944.046.1劣V类水百分比8.29.410.611.712.914.015.216.317.418.4未处理污水前显然不符合标准，因而需要通过搜索算法确定各年的处理污水量由此编制一维搜索程序，可得 问题四的模型结果为:表4-9各年份的年处理废水量(单位：亿吨)年份200520062007200820092

33、0102011201220132014废水量308.9331.78356.34382.56410.46440.03471.26504.17538.75575处理量122145169196224253284317352388(百分比的值为负数即表示实际百分比为零) 根据模型四的求解结果，可求得：表4-10治理后各类水预测百分比:处理后2005200620072008200920102011201220132014IV,V类水百分比10.410.410.410.510.510.410.510.510.410.4劣V类水百分比0000000000由上表可见经过模型给定水量的治理之后，水质已经完全符合

34、“W类和V类水的比例控 制在20%以内，且没有劣V类水”标准，故认为本模型是准确、可行的。4. 5问题五的模型建立与求解长江水污染与经济效应、社会效应和生态效用联系密切，为了全面科学评价长江水 污染的原因及影响，建立长江水污染扩展费用-效益分析(扩展BCA)的评价模型，利 用治污费用最小化原则建立长江水污染治理的优化模型。4. 5. 1可货币化费用-效益分析设沿江经济体第n年排放污水W(n)，当年每亿吨污水所带来经济收益 D(n)，沿江经济体第n年GDP占当年全国GDP (设为G(n)百比分X(n)，当年税收收入占GDP的比率为T(n),用当年沿江经济体的GDP的税收收入表征对政府的排污收益可

35、得：D( n)W( n) = G( n)X( n)T( n)每亿吨污水处理到第i级合格水所需要的费用Bh，n)，污水直排必将抽取地下水补充使用，地下合格水的价格 B2(n)元/吨，则排放W(n)污水对政府引起的费用：B(n) = Bi,i(n) B?(n)W(n)，总收益函数：R(n) =D(n)W(n)T(n) B1,i (n)+B2(n)W(n) =G(n)X(n) -B1,i(n) + B2 (n)W(n)每亿吨污水的单位收益：股=D(n)T(n) -BJn) B?(n)W (n)上式中，Bi,i(n)和B2(n)在污水排放逐年增加、地下水逐年枯竭的情况下将迅速增加，而D G(n)X(n

36、)T(n)，在G(n)相对稳定增长和X(n)相对稳定的情况下将随着 W(n)的 W( n)迅速增加而降低，从而导致单位收益下降，在不加控制的情况下甚至可能出现负增长。参考文献3，代入2004年的相关数据：X(n)=54%，B1,i(n)=0.3元 /吨，T(n)=19%，B2(n) -1.00元/吨，G(n)=1.35615 1013元，W(n) =2.85 1010吨，解得：-R(n) =48.82 _1.3 =47.52元 / 吨；W( n)R(n) =1.35 1012 元；政府收入中的部分是地方政府收入，可见这是地方政府的巨大利益驱动所在。4. 5. 2不可货币化的费用-效益分析不可货

37、币化的长江污染效益：(1) 因国民经济增长带来的投资信心；(2) .对其它地区的经济输出带来的经济繁荣；不可货币化的长江污染费用：(1) .水源枯竭，危及产业发展；(2) .江水污染无法饮用，地下水污染加剧，危及城乡饮用水；(3) .水质污染对人体伤害造成疾病发生率增长，危及人类健康；(4) .湿地面积日益缩减，水的自洁功能日益丧失；(5) .物种受到威胁，珍稀水生物日益灭绝；4. 5. 3对长江污染的综合评价以污染长江为代价换来的短期经济繁荣是一种短视行为，将对长江造成很难恢复的伤害，危及长江的长远综合利用。4. 5. 4对长江治污费用的最小化优化模型参考文献1的优化模型并加以改进得到: m

38、in F - ' ' Fji X jij i二二 C j,l,i > X jl - bj lX jl <1F治理费用的总和，单位：元；Fjl第j个污染源某丨个治理方案的费用，单位：元；Xjl第j个污染源某I个治理方案的权重；Cjl第j个污染源某l个治理方案后在某i个控制点上的浓度，单位：mg/Lbi第i个控制点上的环境目标值，单位： mg/L用上述方法进行治理方案的选取可以反映整个水域的经济合理性，有很好的经济效益、社会效益和环境效益，同时赋予污染源以决策权重，有利于排污单位在平等的市场 经济条件下，开展竞争，提高了排污单位自身的生产效率和挖潜改造的积极性。4. 5

39、. 5建议和意见(1) 建立长江水质的动态监测网，跟踪监测主要污染物的变化，重点监测重污染区，重点治理重度污染区。(2) 改善环保立法，加大执法力度，保证污水处理设备开工率，保证废水排放的合格排 放，从根源上减轻污染。沿江共有 351座污水处理厂，合理的设备开工率，可以从源头 上减轻污染。(3) 开发先进污水处理工艺，借鉴国外先进经验，大力推行水循环经济，将污水处理后供应工农业生产使用，如灌溉农田，工业再利用等。(4) 对长江的污染状况进行定期的科学全面的评估分析，实施长期的水污染防治战略。(5) 在江边设置监管委员会。实行公从参与，其会员最好是由民间的志愿者，或者是大 学里与环保有关系的学生，同时把自来水、矿泉水公司和食品选制企业都组织起来(他 们对水质最敏感)，全面高效率的监督有水污染的企业。(6) 做好绿化及其管理工作。尤其是对上游的绿化及管理，规定离江边一定距离为禁止 放牧区，禁止砍伐区。(7)