长江水质的评价和预测

长江水质的评价和预测摘要长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。考察团对长江沿线城市做了实地考察，收集了长江干流和支流17个地区近两年水质指标的检测数据，揭示了长江的污染情况。本文首先根据沿江17个地区近两年多的水质检测数据，利用主成分分析法，对长江水质情况做了定量的综合评价，并且给出了各地区水质综合排名。然后按照一维水质模型对长江污染情况进行处理，通过高锰酸盐指数（CODMn）和氨氮（NH3-N）在各地区分布以及随时间的变化情况，得出干流各地区污染物贡献量，以此来判断污染源。其次，根据所给10年的主要统计数据，利用改进后的灰色预测模型，在不做污染治理情况下，预测出未来10年长江水质情况，包括各年的污水量和各类水质所占百分比，若未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,可以计算出每年需要处理污水量。最后从国家政策、企业单位和个人三个角度，对解决长江水质污染问题提出了一些建议和意见。关键字:主成分分析法、一维水质模型、灰色预测一问题重述长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”（附件），并发出了“拿什么拯救癌变长江”的呼唤（附件2）。附件3给出了长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上个观测站近一年多的基本数据（站点距离、水流量和水流速）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.10.5之间，比如可以考虑取0.2(单位：1/天)。附件4是“19952004年长江流域水质报告”给出的主要统计数据。下面的附表是国标(GB3838-2002) 给出的地表水环境质量标准中4个主要项目标准限值，其中、类为可饮用水。请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。（4）根据你的预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年需要处理多少污水？ （5）你对解决长江水质污染问题有什么切实可行的建议和意见。附表: 地表水环境质量标准（GB38382002）中4个主要项目标准限值 单位：mg/L序号 分 类标准值 项 目类类类类类劣类1溶解氧(DO)7.5（或饱和率90%）6532 02高锰酸盐指数(CODMn) 24610153氨氮（NH3-N） 0.150.51.01.52.04PH值（无量纲）6-9二模型假设1. 一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。2. 长江水质污染属于一维水质模型。3. 监测站之间水流速取平均值。三符号说明W 单位时间排污量（单位：g/s） 干流观测站当前污染物浓度(单位：mg/L) 观测站水流量（单位： ）上游观测站当前污染物浓度(单位：mg/L)上游观测站水流量（单位： ）K 降解系数（本文取0.2/天）l观测站之间距离（单位：km）v两个观测站之间的平均水流速（单位：m/s）四问题（1）的求解4.1问题分析附件3给出了长江沿线17个观测站（地区）从2003.6月至2005.9月主要水质指标的检测数据，包括监测项目（pH*DOCODMnNH3-N）水质类别和主要污染指标，通过这些数据可知各地区水质的变化情况以及沿江水质分布情况。考虑到监测指标的多种因素，采用主成分分析法1，对多种因素进行综合处理，得出各监测断面水质综合评价结果，这样可以对长江近两年多的水质情况做出定量的综合评价，并分析出各地区水质的污染状况。4.2问题处理4.2.1指标选择各类水质的PH值范围为6-9，并且各地区PH变化不大，因此PH不作为计算指标，而重点选择溶解氧(DO)高锰酸盐指数(CODMn)和氨氮（NH3-N）。4.2.2综合评价函数 各成分参数经因子分析处理后的结果； 各成分系数若存在两个主成分则 各成分的贡献率4.2.3评价结果根据综合评价函数，计算各监测断面和长江水质污染综合得分，计算结果见表4，给出水质污染程度的定量化描述，得分越大表明污染越严重。表1各监测断面水质综合评价结果序号点位名称断面情况水质污染综合得分得分排名1四川攀枝花干流0.99978882重庆朱沱干流（川-渝省界）0.99984973湖北宜昌南津关干流（三峡水库出口）1.00017944湖南岳阳城陵矶干流1.00001655江西九江河西水厂干流（鄂-赣省界）-0.034248099166安徽安庆皖河口干流0.999454107江苏南京林山干流（皖-苏省界）0.999437118四川乐山岷江大桥岷江（与大渡河汇合前）1.00037129四川宜宾凉姜沟岷江（入长江前）0.9894399331310四川泸州沱江二桥沱江（入长江前）1.000007611湖北丹江口胡家岭丹江口水库（库体）0.006808741512湖南长沙新港湘江（洞庭湖入口）0.995599091213湖南岳阳岳阳楼洞庭湖出口1.000405114湖北武汉宗关汉江（入长江前）1.000368315江西南昌滁槎赣江（鄱阳湖入口）0.9891759691416江西九江蛤蟆石鄱阳湖出口0.999752917江苏扬州三江营夹江（南水北调取水口）-0.06618083917通过对照地表水环境质量标准以及上表的结果，可知长江中下游如湖南岳阳岳阳楼，四川乐山岷江大桥，湖北武汉宗关，湖北宜昌南津关，湖南岳阳城陵矶等地区水质较差，水质在类之间。从水质监测数据来看，主成分分析结果较真实的反映了实际情况。五问题（2）的求解5.1问题分析通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，反映江河自然净化能力的指标称为降解系数。可以认为长江是一维水质模型2，已知干流上个观测站近一年多的基本数据（站点距离、水流量和水流速）以及高锰酸盐指数和氨氮的检测数据。通过这些数据可以知道各地区在单位时间的排污量，称之为排污贡献量，排污贡献量越大，代表该地区是污染源。5.2模型建立检测点污染物浓度: （5.1）由式（5.1）可推导出各地区排污贡献量: （5.2）5.2 模型求解通过模型公式（5.2）可以计算出各观测站主要污染物高锰酸盐和氨氮的排污贡献量。结果见下表。表2各观测地区污染物贡献量（单位：g/s）月份污染指标四川攀枝花重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京2004.04CODMn8487483004620084480814908260044700NH3-N553.52344.52772986931-1097-58942004.05CODMn159962882059400738009238010880062100NH3-N260.43537594065606556850031052004.06CODMn10025426007714079100767005136066200NH3-N160.42556690279104425545719862004.07CODMn111845412072640101220702005423057780NH3-N186.42952363286764320701816052004.08CODMn216922120074400103600738306346059670NH3-N37401802456085477704534417552004.09CODMn38308209440181900209820138320118719.1186300NH3-N565.2761615515161401820010387.8864802004.10CODMn2608259207067078050471206510072960NH3-N260.83726439380284712527015362004.11CODMn4200163402014031200321203740035280NH3-N902205.925444320321220409802004.12CODMn1521.685151776037450303602397029800NH3-N76.083340.5118433171584197414902005.01CODMn854.45226868341769261603198036000NH3-N49.842170.81005.43439.82289246031682005.02CODMn550.86485.4902025536319303699034730NH3-N91.81981.65811.83431.41339246646812005.03CODMn685.390061087828864257405778038700NH3-N161.982559.6984.22956.8157349224302005.04CODMn706.273001296020996302004646033150NH3-N64.22299.57023040.8211430304420通过上表数据可以画出各检测地区污染物柱状图.图表1 各地区2004.04-2005.04内CODMn变化情况图表2 各地区2004.04-2005.04内NH3-N变化情况由上表CODMn和NH3-N变化情况可知，湖北宜昌、湖南岳阳和江西九江污染物贡献量较大，可以判定这些地区是污染源，即长江主要污染物高锰酸盐指数和氨氮的污染源集中在中游地区。六、问题（3）（4）的求解6.1传统的模型传统的模型的建模过程如下所示：原始数据记为，为增强数据的规律性，对数据进行依次累加，得到的累加生成数列记为，那么和满足： -(1)称为原始数列的1-AGO数列。取的加权均值，则，为确定参数，一般情况取，记，如果将生成数列的时刻看成连续的变量，又将生成数列看成关于时间的函数，即，那么只要生成数列对的变化率有影响，就可以建立下面的白化微分方程模型： -(2)其中是发展灰度，是内生控制灰度，由于，取为灰导数，为背景值，则将方程(2)相应的灰微分方程为：即矩阵形式为：其中，用最小二乘法求得参数的估计值为于是方程(2)有特解 则 - (3)将时间参数代入（3）式，即可得出原始数列的模拟值及预测值。6.2 模型的改进二次拟合传统的模型往往计算误差较大，无法满足要求，通过验算也发现其不适用于长江排污量及水变化质等数据的预测。而通过二次模拟改进后的预测模型，则大大提高了计算精度。以下是二次模拟灰色预测模型的建模过程：将一次模拟所得的模型记为- (4)根据第一次模拟所得的a值及原始1-AGO数列对M和N进行估计。由于写成矩阵形式即为其中- (5)由最小二乘法可得- (6)求出M和N后即可得到更精确的二次拟合参数模型- (7)则- (8)6.3模型的建立利用灰色理论，建立改进的二次灰色预测模型。把数据中不明确的随机量看成灰色数，利用数据分析处理的方法去寻找数据里面隐含的规律。6.4模型的求解6.4.1数据的预处理表3 1995-2004年长江流量及排污量统计表总流量（亿立方米）废水排放总量（亿吨）940528519948.3443.946031729930.51674.678399439847.05912.508952249698.351157.65605759484.851410.34479769206.971670.80718378865.141939.28235988459.772216.0168397991.292501.264682107460.12795.287818表4 19952004年长江水质所占百分比统计表水文期类类类类类劣类199525.842.624.73.930199615.320.249.89.71.93.1199712.224.943.613.32.63.4199811.524.152.88.31.71.619995.239.835.29.56.24.120005.632.835.616.64.45.320015.933.134.7145.56.820024.44428.3103.21020034.741.531.36.45.810.320041.226.939.914.85.911.3表3及表4中的各项数据均波动较大，须经过预处理才能使用。本文采用弱化因子对原数据序列进行弱化处理。原始数据记为：弱化序列记为：即其中，- (9)将表3及表4中的各项数据代入（9）式进行预处理，处理结果见表5。年度2005200620072008200920102011201220132014总流量9894.11 9970.67 10027.88 10150.26 9654.13 9682.36 9621.95 9865.00 9692.50 9405.00 排污量219.75 224.83 230.56 237.36 245.42 253.10 257.88 270.33 277.50 285.00 类9.187.336.345.54.54.364.053.432.951.2类32.9931.9233.3934.636.3535.6636.3837.4734.226.90类37.5939.0237.6836.8334.1733.9633.5533.1735.639.9类10.6511.411.6111.3711.8812.3611.310.410.614.8类4.024.1340414.675.174.965.14.975.855.9劣类5.596.216.67.067.978.749.6010.5310.811.36.4.2二次拟合模型的求解将表5中经过弱化处理后的数据经过两次模拟，即可求得预测公式(8)，并可据此得出19952004年各项数据的模拟值，并可预测未来黄河总流量、排污量、水质断面比例的各项数据。长江总流量、排污量及各类水质断面比例关于时间的函数关系模型见表6，对2005-2014年十年间各项数据的预测成果见表7、表8。对灰色预测模型进行检验，以类水质所占百分比为例，均方差比值(C)为0.2087，小误差概率(p)为1.09，可认为达到一般要求，由此结果可看出灰色预测模型的预测值较为贴近真实值，误差较小。表6 二次拟合得到的函数关系模型项目（M，N）函数关系模型总流量(27288,7389)27288exp(0.1587t)+(7389)排污量(593.9032，-154.8429)593.9032exp(0.1742t)+(-154.8429)各类水质断面所占比例(-55.3307,64.2834）-55.3307exp(-0.1453t)+(64.2834)类（-6624.1,6653.5）-6624.1exp(-0.0054t)+(6653.5)类（-6951.7，6992.2）-6951.7exp(-0.0052t)+(6992.2)类（931.672,-920.334）931.672exp(0.0118t)+(-920.334)类（105.5166,-101.507）105.5166exp(0.0396t)+(-101.507)劣类（77.6697,-72.3957）77.6697exp(0.0782t)+(-72.3957) 表7 二次拟合模型得到的各项指标预测值年度2005200620072008200920102011201220132014排污量267.08134.52160.12190.59226.86270.03321.41382.58455.38542.03类15.246.475.64.844.183.623.132.712.342.02类38.1735.4835.2935.134.9134.7234.5434.3534.1733.98类36.6635.8735.6835.535.3135.1334.9534.7734.5934.41类7.611.1911.3211.4611.5911.7311.8712.0112.1512.30类2.374.434.614.84.995.25.415.625.856.09劣类0.296.837.397.998.649.3410.110.9210.8112.77 对表7中各类水质断面所占比例进行归一化处理（即：某一年某一类水质的预测值/当年各类水质预测值之和*100%），即得出给类水质的百分比，见表8。表8 二次拟合模型得到的各类水质比例预测值（归一化处理后）类15.196.455.64.854.23.633.132.72.321.99类38.0535.3835.3335.2135.0434.8134.5434.2233.8633.46类36.5435.7735.7235.6135.4435.2234.9534.6334.2833.88类7.5711.1611.3411.4911.6411.7611.8711.9712.0412.11类2.364.424.624.825.015.215.415.605.805.99劣类0.296.817.48.018.679.3610.1010.8811.7012.57 图3 2005-2014年期间长江排污量 图4 20052014年期间长江各类水质所占比例发展趋势图从表7及图3可以看出，未来十年内长江排污量有先减少后逐年增加的趋势。从表8及图4可以看出，未来十年内长江优质水质(类、类及类水质)所占比例呈逐年下降趋势，劣质水质(类、类及劣类水质)所占比例呈逐年上升趋势。这与排污量逐年递增的规律保持一致。排污量增加速率较快，虽然类、类及类水质断面所占比例在未来十年内仍然保持在较高的水平，但是优质水质所占比例也在逐年递减，并且劣质水质所占比例在逐年增加，到2014年类、类及劣类水质总比例将增加到40%。说明若不采取有力的治理措施，未来十年长江水质将面临巨大威胁，未来长江污染治理的任务将十分艰巨。综合上述分析结果，未来十年长江年排污量呈增加趋势，优质水质断面所占比例逐年降低，而劣质水质断面所占比例逐年上升，长江水质正逐渐变坏。这说明近年若不加紧治理，长江水质令人担忧。6.5由上述预测可知未来10年长江年废水排放量如下表：表9 未来10年长江年废水排放量（单位：亿吨）2005200620072008200920102011201220132014267.08134.52160.12190.59226.86270.03321.41382.58455.38542.03要求未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水,那么每年需要处理的污水量为：（类和类占的百分比-20%+劣类占的百分比）*年废水排放量有预测的未来10年类类和劣类水所占比例及废水排放量，带入上述公式可得每年需处理的污水量如下表：表10 未来10年需处理的污水量（单位：亿吨）20052006200720082009201020112012201320140.789.1611.8415.2719.6725.2932.4741.6353.368.15七解决长江污染问题的一些建议从上述统计数据可知，若不采取相关措施治理长江污染，任凭“一江浊水向东流”，10年后的长江就会变得面目全非，生态恶化，直接影响沿江地区的经济发展以及生命健康。为治理长江污染，本文有如下建议：（1） 从国家的角度来说，应该出台相关政策，严厉打击长江沿岸违规建厂违规排污行为，同时加大环保的投入，对有关工厂处理污水设施给予补助。（2） 从工厂及社会的角度来说，工厂应该按规定处理污水，不达标不排放，增加科技投入，高效利用资源，减少污染；对于社会媒体及相关单位，应该加大宣传力度，对长江上游水土流失，中游污染严重等现象进行曝光和呼吁改善。（3） 从个人角度，尤其沿江居民。更应该保护长江，呼吁对长江污染的治理，监督举报污染长江的行为。同时长江上游地区居民应该加大植树造林力度，减少水土流失。八参考文献1 吴孟达,成礼智.数学建模教程.高等教育出版社.2011-8-19.2 余江,许刚,康建雄.河流一维水质预测模型在污染贡献值计算中的应用.科学环境与技术.第29卷.32-33.2006-8.3 石志广,唐玲艳.MATLAB工程计算.清华大学出版社.2008-7九附录MATLAB程序说明：1、求灰度值function c7fun73X0=;%弱化后的数据AU=c7fun72(X0); a=AU(1); u=AU(2); m2=length(X0);for k=0:1:m2-1 xx1(k+1)=(X0(1)-u/a)*exp(-a*k)+u/a;ends=0;xx0(1)=X0(1);for jj=2:1:m2; xx0(jj)=xx1(jj)-xx1(jj-1);enddisp(a); %灰度值AU(1)function au=c7fun72(X0)m=length(X0);s1=0;for jj=1:1:m; X1(jj)=s1+X0(jj); s1=X1(jj); endfor ii=1:1:m-1; B(ii)=-(X1(ii)+X1(ii+1)/2; endB=B(:),ones(m-1,1);y=X0(2:m);au=inv(B*B)*B*y;2、对M、N求解function result=calc()n=10;g=zeros(n,2);a1= ; %所求数组a=-0.0901;for i=1:n g(i,1)=exp(-a*(i-1); g(i,2)=1;endresult=inv(g*g)*g*a1;end长江流域主要城市水质检测报告发布日期;2003-06序号点位名称断面情况主要监测项目(单位:mg/L)水质类别主要污染指标pH*DOCODMnNH3-N本月上月1四川攀枝花干流7.66.80.20.1II2重庆朱沱干流（川-渝省界）7.638.412.80.34II3湖北宜昌南津关干流（三峡水库出口）7.077.815.80.55III4湖南岳阳城陵矶干流7.586.472.90.34II5江西九江河西水厂干流（鄂-赣省界）7.346.191.70.13II6安徽安庆皖河口干流7.526.543.20.22II7江苏南京林山干流（皖-苏省界）7.786.93.10.11II8四川乐山岷江大桥岷江（与大渡河汇合前）7.664.25.80.53IV溶解氧9四川宜宾凉姜沟岷江（入长江前）8.017.632.40.25II10四川泸州沱江二桥沱江（入长江前）7.634.023.61.06IV溶解氧、氨氮11湖北丹江口胡家岭丹江口水库（库体）8.6310.21.80.1I12湖南长沙新港湘江（洞庭湖入口）7.426.454.30.99III13湖南岳阳岳阳楼洞庭湖出口7.736.261.40.21II14湖北武汉宗关汉江（入长江前）86.432.40.17II15江西南昌滁槎赣江（鄱阳湖入口）6.645.181.10.92III16江西九江蛤蟆石鄱阳湖出口7.286.872.70.15II17江苏扬州三江营夹江（南水北调取水口）7.296.91.60.15II发布日期;2003-07序号点位名称断面情况主要监测项目(单位:mg/L)水质类别主要污染指标pH*DOCODMnNH3-N本月上月18四川攀枝花龙洞干流8.38.15.60.1IIIII19重庆朱沱干流（川-渝省界）7.478.541.60.25IIII20湖北宜昌南津关干流（三峡水库出口）7.737.895.20.22IIIIII21湖南岳阳城陵矶干流7.787.573.90.31IIII22江西九江河西水厂干流（鄂-赣省界）7.066.252.30.18IIII23安徽安庆皖河口干流7.636.7330.28IIII24江苏南京林山干流（皖-苏省界）7.416.372.50.1IIII25四川乐山岷江大桥岷江（与大渡河汇合前）7.614.634.50.92IVIV溶解氧26四川宜宾凉姜沟岷江（入长江前）8.098.0880.37IVII高锰酸盐指数27四川泸州沱江二桥沱江（入长江前）7.94.76.40.37IVIV高锰酸盐指数、溶解氧28湖北丹江口胡家岭丹江口水库（库体）8.179.811.90.09II28湖南长沙新港湘江（洞庭湖入口）7.646.932.60.4IIIII30湖南岳阳岳阳楼洞庭湖出口7.866.493.30.26IIII31湖北武汉宗关汉江（入长江前）85.383.70.14IIIII32江西南昌滁槎赣江（鄱阳湖入口）6.714.8511.01IVIII溶解氧33江西九江蛤蟆石鄱阳湖出口7.686.082.60.19IIII34江苏扬州三江营夹江（南水北调取水口）7.314.491.60.32IVII溶解氧发布日期;2003-08序号点位名称断面情况主要监测项目(单位:mg/L)水质类别主要污染指标pH*DOCODMnNH3-N本月上月35四川攀枝花龙洞干流8.248.10.50.15IIII36重庆朱沱干流（川-渝省界）7.738.551.80.2IIII37湖北宜昌南津关干流（三峡水库出口）7.86.652.80.31IIIII38湖南岳阳城陵矶干流7.897.922.60.32IIII39江西九江河西水厂干流（鄂-赣省界）8.576.8830.08IIII40安徽安庆皖河口干流7.617.0720.32IIII41江苏南京林山干流（皖-苏省界）7.596.4720.12IIII42四川乐山岷江大桥岷江（与大渡河汇合前）7.575.494.71.93VIV氨氮43四川宜宾凉姜沟岷江（入长江前）7.778.6740.34IIIV44四川泸州沱江二桥沱江（入长江前）7.197.162.40.39IIIV45湖北丹江口胡家岭丹江口水库（库体）7.729.032.50.07III46湖南长沙新港湘江（洞庭湖入口）6.294.342.90.92IVII溶解氧47湖南岳阳岳阳楼洞庭湖出口7.258.323.70.37IIII48湖北武汉宗关汉江（入长江前）7.915.873.60.25IIIIII49江西南昌滁槎赣江（鄱阳湖入口）6.785.351.72.18劣VIV氨氮50江西九江蛤蟆石鄱阳湖出口6.76.673.50.16IIII51江苏扬州三江营夹江（南水北调取水口）7.096.023.80.19IIIV发布日期;2003-09序号点位名称断面情况主要监测项目(单位:mg/L)水质类别主要污染指标pH*DOCODMnNH3-N本月上月52四川攀枝花龙洞干流8.348.81.10.11II53重庆朱沱干流（川-渝省界）7.638.581.50.24IIII54湖北宜昌南津关干流（三峡水库出口）7.5510.63.60.36IIII55湖南岳阳城陵矶干流7.938.363.10.3IIIII56江西九江河西水厂干流（鄂-赣省界）7.67.432.30.14IIII57安徽安庆皖河口干流7.676.950.39IIIII58江苏南京林山干流（皖-苏省界）7.836.611.50.04IIII59四川乐山岷江大桥岷江（与大渡河汇合前）8.026.323.61.41IVIV氨氮60四川宜宾凉姜沟岷江（入长江前）7.969.122.60.27IIII61四川泸州沱江二桥沱江（入长江前）7.447.733.40.53IIIIII62湖北丹江口胡家岭丹江口水库（库体）7.527.112.20.12IIII63湖南长沙新港湘江（洞庭湖入口）7.66.72.50.75IIIIII64湖南岳阳岳阳楼洞庭湖出口6.978.83.50.31IIII65湖北武汉宗关汉江（入长江前）7.956.322.50.35IIIII66江西南昌滁槎赣江（鄱阳湖入口）6.784.040.773.06劣V劣V氨氮67江西九江蛤蟆石鄱阳湖出口6.977.142.80.26IIII68江苏扬州三江营夹江（南水北调取水口）8.035.761.10.41IIIII发布日期;2003-10序号点位名称断面情况主要监测项目(单位:mg/L)水质类别主要污染指标pH*DOCODMnNH3-N本月上月69四川攀枝花龙洞干流8.378.63.10.17I70重庆朱沱干流（川-渝省界）7.729.022.10.28II71湖北宜昌南津关干流（三峡水库出口）7.2211.92.40.26II72湖南岳阳城陵矶干流7.918.0160.27II73江西九江河西水厂干流（鄂-赣省界）7.417.663.30.13II74安徽安庆皖河口干流7.757.652.90.28III75江苏南京林山干流（皖-苏省界）7.736.891.90.16II76四川乐山岷江大桥岷江（与大渡河汇合前）7.975.825.70.76IV77四川宜宾凉姜沟岷江（入长江前）7.9910.12.60.32II78四川泸