本文首先對(duì)各觀測(cè)站（地區(qū)）的污染狀況進(jìn)行了評(píng)價(jià)，綜合考慮了DO，的確定方法，基本思想是將支流并入到干流處理，支流與干流在各個(gè)觀測(cè)點(diǎn)匯進(jìn)行了分析，最終得到整個(gè)長(zhǎng)江流域的各項(xiàng)水質(zhì)評(píng)價(jià)指標(biāo)的綜合評(píng)估結(jié)果：近兩年多的水質(zhì)類型總體上表現(xiàn)為Ⅱ類CODMnNH3-N的排放量和自然降解因在預(yù)測(cè)分析未來(lái)十年長(zhǎng)江的水質(zhì)污染發(fā)展趨勢(shì)時(shí)，先用逐步回歸法研究時(shí)廢水排放總量是唯一的對(duì)各類河長(zhǎng)百分比影響顯著的變量后，運(yùn)用多項(xiàng)式擬合得到了未來(lái)十年每年的廢水排放總量和長(zhǎng)江總流量。最后，采用動(dòng)態(tài)規(guī)劃和排隊(duì)論思想（見(jiàn)圖十五利用平衡轉(zhuǎn)移速率圖流域水文年在200~2014年各類水河長(zhǎng)百分比，其中205年Ⅰ類占0.7937%，Ⅱ類水占18.266%，Ⅲ類水占21.6882%，Ⅳ類水占在2005~2014年各類水河長(zhǎng)百分比矩陣。運(yùn)用矩陣迭代運(yùn)算2005~2014年每年的污水處理可行性方案（具體數(shù)據(jù)見(jiàn)表七、八、九：時(shí)間滯后性、逐步回歸法、動(dòng)態(tài)規(guī)劃、排隊(duì)論、平衡轉(zhuǎn)移速率基本情況

問(wèn)題提了相關(guān)部門(mén)和專家們的高度重視。斷提前，長(zhǎng)江斷流日益近；水質(zhì)嚴(yán)重，危及沿江許多城市的飲用水，相關(guān)信息2：國(guó)標(biāo)(GB3838-2002)4個(gè)主要項(xiàng)信息3：過(guò)去10長(zhǎng)江流域在不同時(shí)段的水質(zhì)情況需要解決的問(wèn)101020%以內(nèi)，且沒(méi)有劣Ⅴ類水,那么每年需要處理多少污水？問(wèn)題分、需要將已有的歷史數(shù)據(jù)進(jìn)行合理歸類、修正，然后對(duì)長(zhǎng)江的水質(zhì)作出分析、問(wèn)題一的分析2727個(gè)月獲得的水質(zhì)的各項(xiàng)指標(biāo)，包括溶解氧，高錳酸鹽濃度，氨氮濃度，PH值。712個(gè)月的水文資料，包括水流量，流速以及各觀測(cè)這里采用國(guó)標(biāo)(GB3838-2002)和指標(biāo)評(píng)價(jià)原則作為對(duì)長(zhǎng)江水質(zhì)情況的2的分站處排入江河的污染物的數(shù)量可以假設(shè)就是在該觀測(cè)站處檢測(cè)算得的污染物的問(wèn)題三的分10枯水季節(jié)全流域，干流，支流的水質(zhì)狀況（各類水的比例豐水季節(jié)全流域，干流，支流的水質(zhì)狀況（各類水的比例水文年全流域，干流，支流的水質(zhì)狀況（各類水的比例從附表3中我們發(fā)現(xiàn)，與各類水河長(zhǎng)百分比可能相關(guān)的因素有時(shí)段、評(píng)問(wèn)題四的分(1)長(zhǎng)江在過(guò)去十年的水質(zhì)狀況問(wèn)題五的分3、長(zhǎng)江水質(zhì)綜合評(píng)價(jià)和各地區(qū)水質(zhì)污染情況分1基本假設(shè)t3符號(hào)設(shè)定i:數(shù)據(jù)采樣時(shí)所處的月份（200361j:（3k:各項(xiàng)水質(zhì)指標(biāo)的編號(hào)（PH1;DO2;CODMn:3;NH3-N4XijkjkikkWjj模型建立（3.4各地區(qū)水質(zhì)總體指標(biāo)的確對(duì)兩年多來(lái)每個(gè)觀測(cè)點(diǎn)的各個(gè)水質(zhì)指標(biāo)（PH,DO,CODMn,NH3-N的測(cè)Xjk1l

X

（其中的為當(dāng)?shù)厮|(zhì)的類別。例如：當(dāng)某測(cè)試點(diǎn)兩年來(lái)DO的均值是6.8（符合Ⅱ，CODMn，NH3-各水質(zhì)指標(biāo)的歸類函數(shù)HjkGXjk，其中GXjk)H11為Ⅱ類。HjminHk長(zhǎng)江水質(zhì)綜合評(píng)價(jià)函數(shù)根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)其中兩個(gè)指標(biāo)CODMn和NH3-N的評(píng)價(jià)標(biāo)準(zhǔn)WHWT*H權(quán)重的確定3可以得到長(zhǎng)江干流上每個(gè)觀測(cè)點(diǎn)的水流量，結(jié)合各觀測(cè)點(diǎn)兩年來(lái)CODMnNH3-N濃度的數(shù)據(jù)，則污染物的濃度與水流量的積可以看作是單位 供水量V=V2-供水量V=V3- 這個(gè)假設(shè)可以避免我們對(duì)支流水質(zhì)情況的討論因?yàn)閮捎^測(cè)點(diǎn)間的各支流模型求解為行文的方便起見(jiàn)這里僅給出了在運(yùn)用計(jì)算平均值方法下采用水質(zhì)評(píng)（：、各地區(qū)水質(zhì)污染狀況7個(gè)城市均達(dá)到了二級(jí)標(biāo)準(zhǔn)，V類型，特別差，查找程序運(yùn)行時(shí)得出的數(shù)據(jù)可以4.6332mg/L，與國(guó)家標(biāo)準(zhǔn)對(duì)比表發(fā)現(xiàn)該值確實(shí)是嚴(yán)重。、對(duì)長(zhǎng)江水質(zhì)的定量綜合評(píng)價(jià)1（DO=7.7340mg/L,（CODMn=2.4661氨氮（NH3-N）0.2144mg/L模型的進(jìn)一步討論分析研究長(zhǎng)江干流主要污染源在哪些地補(bǔ)充假設(shè)如下有關(guān)符號(hào)的約定iijjVijijvijijmijijPjjMijijrT：?jiǎn)挝粫r(shí)間（月Sj,j1：在第i個(gè)月時(shí)第j觀測(cè)點(diǎn)到第j+1個(gè)觀測(cè)點(diǎn)的距vij,j1：在第i個(gè)月時(shí)水流從第j觀測(cè)點(diǎn)到第j+1個(gè)觀測(cè)點(diǎn)的速tij,(j1)ijj+1Cijij模型的建ij（即Mij） 月的水流量V（m3s）與該觀測(cè)點(diǎn)該月的污染物濃度C（mgL=g/m3 個(gè)月的時(shí)間T（30×24×3600s）

Vij*T

MijMij

下面來(lái)確定tijj1。我們假定水流速vij相鄰觀測(cè)點(diǎn)間的時(shí)間可由定積分S

ti

(s)dt確定。

(sLagrange 0 出水流經(jīng)任意兩相鄰觀測(cè)點(diǎn)間的時(shí)間tiijj+1vijj1 vijj1(vijvij1)/

tij,(j1)=Sj,j1/[(vijvij1)/

綜合以上四個(gè)式子可以得到在第ij

ij附表3有長(zhǎng)江各干流觀測(cè)點(diǎn)12個(gè)月各污染物的濃度，計(jì)算各觀測(cè)各月份排對(duì)時(shí)間滯后性的研究時(shí)間間隔為1個(gè)月的話某月在上游觀測(cè)點(diǎn)中得到的數(shù)據(jù)和相鄰下游觀測(cè)點(diǎn)得到的數(shù)據(jù)是不能夠運(yùn)用上述模型所獲得的（5）式進(jìn)行運(yùn)用的，因?yàn)榈趇個(gè)過(guò)上游觀測(cè)點(diǎn)的水流在第i+1個(gè)月才到達(dá)相鄰下游觀測(cè)點(diǎn)，而相鄰下游觀測(cè)點(diǎn)第i個(gè)月的水流是第i-1個(gè)過(guò)上游觀測(cè)點(diǎn)的水流。因此（5）式應(yīng)當(dāng)作進(jìn)一步改

1rtij1,

i1,j這僅僅是一個(gè)簡(jiǎn)單的討論，因?yàn)樗髁鹘?jīng)兩相鄰地區(qū)的時(shí)間間隔不是固定安慶-4.16736.27675.01545.69941.69014.3096模型的求解我們求得12個(gè)月各觀測(cè)點(diǎn)水速以及各觀測(cè)點(diǎn)污染物濃度的期望，以此代表（1（2（3（4）測(cè)點(diǎn)排放污染物的數(shù)量，運(yùn)用程序求解得到了所要求的主要污染物的（2）高錳酸鹽污染源的主要分布地區(qū)氨氮污染源的主要分布地區(qū)說(shuō)明：在用附件2求氨氮污染源的主要分布地區(qū)時(shí)，矩陣C應(yīng)換成各觀測(cè)點(diǎn)對(duì)長(zhǎng)江未來(lái)水質(zhì)污染發(fā)展趨勢(shì)的預(yù)測(cè)分補(bǔ)充假設(shè)符號(hào)設(shè)定1x1

,x2

x3 2xx

xx

,x63、x7：該年長(zhǎng)江總流量4、x8：該年廢水排放總量模型的建逐步回歸法研究各個(gè)自變量對(duì)各類水河長(zhǎng)百分比的影我們認(rèn)為所列各項(xiàng)因素 x8不一定每一項(xiàng)都對(duì)Ⅰ類河長(zhǎng)的百分比有顯百分比與各項(xiàng)因間的模型，這個(gè)模型應(yīng)盡量簡(jiǎn)單和有效，并且由此能給長(zhǎng)江（見(jiàn)文獻(xiàn)[2]323-324）利用統(tǒng)計(jì)工具箱中的逐步回歸命令stepwise，以某類河長(zhǎng)百分比因變量，各項(xiàng)因素 x8為自變量，得到如圖窗口圖七、8StepwisePlotNextstep“Move

圖八、僅含x8模型的Stepwise窗（5.3.3）預(yù)測(cè)未來(lái)十年每年的長(zhǎng)江總流觀察圖九可以發(fā)現(xiàn)，1998得到剔除掉異常點(diǎn)后的擬合曲線（圖十預(yù)測(cè)未來(lái)十年每年的廢水排放總表三、未來(lái)十年每年的廢水排放總量（單位：億噸趨勢(shì)，其中，200420051.9運(yùn)用動(dòng)態(tài)規(guī)劃和排隊(duì)論思想預(yù)測(cè)長(zhǎng)江未來(lái)水質(zhì)污染趨勢(shì)一、選取過(guò)去十年水文年全流域各類河長(zhǎng)百分比數(shù)據(jù)構(gòu)造初始矩陣A，并A 3.13.41.64.1A Ai1995i5.36.810 再選取過(guò)去十年第i類水每相鄰兩年的變化率作為未來(lái)十年第i+1類水的每相鄰頁(yè)水水水水水從矩陣A的第一列元素可以看出在過(guò)去十年中Ⅰ類水從所占的25.8%1.2%，顯然如果對(duì)廢水排放總量不加控制的話，未來(lái)十年中Ⅰ類水從所水水水水水20041995百分比十分相似（19525.8%20426.9%102004年水體循環(huán)系統(tǒng)的狀況：水水水水圖十六、2004水水水2014年的水體循環(huán)系統(tǒng)狀況B：A中第六列元素，因?yàn)榱英躅愃荒軌蛟傧騽e的水類型轉(zhuǎn)化，得A。 A' A'第二行的數(shù)據(jù)除以第一行的數(shù)據(jù)，得到各類水在未來(lái)一年所占總BA'BBA'（B4得到）得到轉(zhuǎn)移矩陣BB=用初始條件矩陣T.*B（TB相應(yīng)元素相乘，就可以得到未來(lái)十2005200620072008200920102011201220132014表四、200510年的Ⅱ~200521II表五、200510I22005~2014I1，我們利用歸一化思想進(jìn)行標(biāo)準(zhǔn)化處理。5%2005200620072008200920102011201220132014表六、2005105.4模型的進(jìn)一步討論6、給定控制條件下的對(duì)每年需要處理污水?dāng)?shù)量的分模型求1：2005200620072008200920102011201220132014以避免在求轉(zhuǎn)移矩陣時(shí)出現(xiàn)除數(shù)為0的情形。2：（200520%2005~2014行性方案。20052005199620052006201520062006920052014從上面表格數(shù)據(jù)我們可以直接得到2005（數(shù)據(jù)的意義為劣Ⅴ類水占總河長(zhǎng)的百分比00表七、2005之后，用處理后的2005年的數(shù)據(jù)作為轉(zhuǎn)移矩陣中第一行元素去計(jì)算得到20060表八、20062007~201400000000000表九、2007~2014模型的評(píng)價(jià)和改進(jìn)方向水水水V7、對(duì)解決長(zhǎng)江水質(zhì)污染的切實(shí)可行性建議和意養(yǎng)育4億人，南水北調(diào)后便是8億人。長(zhǎng)江沿岸GDP占總量54%。長(zhǎng)江是中國(guó)資源帶、、產(chǎn)業(yè)帶，長(zhǎng)江就是我國(guó)經(jīng)濟(jì)脈[4]。但是長(zhǎng)江現(xiàn)在卻著前所未有的六大：森林覆蓋率嚴(yán)重下降，泥沙含量增加，生態(tài)環(huán)境急??；枯水期不斷提前，長(zhǎng)江斷流日益近；水質(zhì)嚴(yán)重惡化，危及沿江許多城市的飲用水，各種疾病肆虐沿江城鄉(xiāng)；物種受到，珍稀根據(jù)預(yù)測(cè)：長(zhǎng)江已陷入深度，若不及時(shí)拯救，10年之內(nèi)，長(zhǎng)江水系生態(tài)將瀕臨。解決長(zhǎng)江水污染問(wèn)題已經(jīng)刻不容緩了。決長(zhǎng)江的水質(zhì)污染的最的問(wèn)題就是要控制污染源。2的老路。針對(duì)此情況在當(dāng)?shù)厝罕姰?dāng)保教育是非常迫切的。只有從觀念上堵住污染的，才可以徹底控制污染。相應(yīng)的措施，將污染事故產(chǎn)生的降到最低。在長(zhǎng)江沿岸建設(shè)場(chǎng)，防治洪水季節(jié)被沖入江中污染江水。對(duì)港口和船舶實(shí)行污染責(zé)任制，向江中排放污染物鼓勵(lì)農(nóng)民使用生物和有機(jī)肥這樣可以減少和氨氮類化肥對(duì)將每個(gè)地區(qū)的環(huán)保改善狀況列入當(dāng)?shù)氐恼?jī)掛鉤使他們加大對(duì)[1]胡知能、運(yùn)籌學(xué)（第二版，：科學(xué)，2004年月姜啟源謝金星葉俊數(shù)學(xué)模型（第三版高等教育，2004年8雷功炎數(shù)學(xué)模型講義，，2004年1周末報(bào)網(wǎng)“長(zhǎng)江還有十年”htt 2005年9月18日附件：1：%各地區(qū)水質(zhì)污染情況和長(zhǎng)江水質(zhì)綜合評(píng)價(jià)的實(shí)現(xiàn)%其中矩陣Aii個(gè)月份各地區(qū)各水質(zhì)指標(biāo)的觀測(cè)值，這里沒(méi)有列出。fori=1:17if(A(i,1)>=6&&A(i,1)<=9)elseB(i,1)=0;if(A(i,2)>=7.5)elseif(A(i,2)>=6)elseif(A(i,2)>=5)elseif(A(i,2)>=3)elseif(A(i,2)>=2)elseif(A(i,2)>=0)B(i,2)=0;

elseif(A(i,3)<=4)elseif(A(i,3)<=6)elseif(A(i,3)<=10)elseif(A(i,3)<=15)

if(A(i,4)<=0.15)elseif(A(i,4)<=0.5)elseif(A(i,4)<=1)elseif(A(i,4)<=1.5)elseif(A(i,4)<=2)

fori=1:17forif(B(i,j)<C(i,1)) 攀枝花 朱沱 宜昌 岳陽(yáng) 九江 安慶 樂(lè)山 宜賓 瀘州沱江 丹江口水庫(kù) 湘江 洞庭湖 漢江gtext('劣 贛江 鄱陽(yáng)湖 2：%M2表示近一年多單位時(shí)間內(nèi)流經(jīng)第j個(gè)觀測(cè)點(diǎn)的高錳酸鹽污染物的平均質(zhì)量，單位為g/s%V2表示近一年多單位時(shí)間內(nèi)流經(jīng)第j個(gè)觀測(cè)點(diǎn)的高錳酸鹽污染物的平均水流速單位為m/s%V3表示相鄰兩個(gè)觀測(cè)點(diǎn)間的平均水流%M22表示近一年多在單位時(shí)間內(nèi)由各觀測(cè)點(diǎn)注入的高錳酸鹽污染物的平均質(zhì)%C表示一年多各個(gè)月份各觀測(cè)點(diǎn)高錳酸鹽的濃%U%v表示一年多各個(gè)月份各觀測(cè)點(diǎn)高錳酸鹽的水流C=[2.33.52.23.32.92.83322422222 951712612623642 951712612