寧夏黃河水系干流水質監測斷面優化

侯凱+李亞斌+錢會+楊咪

摘要：黃河是寧夏地區最重要的地表水資源，所以整個黃河寧夏段水環境質量的優劣將影響寧夏社會經濟的發展。基于2010年-2014年寧夏環境監測站的統計資料，以污染最嚴重的指標COD、氨氮、總磷的監測資料為基礎，采取物元分析法、模糊聚類法及WASP模型對黃河干流寧夏段水質監測斷面的布設進行優化。結果表明：現有的6個國控監測斷面對黃河干流水質變化狀況的監測不完善，須增設2個監測斷面，分別位于下河沿與金沙灣國控斷面間及葉盛公路橋與銀古公路橋國控斷面間，地理位置分別處于石空鎮與臨河鎮，即105°40′31″E，37°34′04″N與106°18′36″E，38°16′11″N。

關鍵詞：水質監測；斷面優化；物元分析；模糊聚類；WASP模型

中圖分類號：X832 文獻標識碼：A 文章編號：1672-1683（2017）03-0100-08

Abstract：The Yellow River is the most important surface water resource in Ningxia，so the water environment quality will directly affect the social and economic development of Ningxia.Based on the statistical data from the environmental monitoring station in Ningxia from 2010 to 2014，using the monitoring data of the main pollutants such as chemical oxygen demand，ammonia nitrogen，and total phosphorus，we employed matter-element analysis，fuzzy clustering，and WASP model to optimize the layout of the Yellow River water quality monitoring cross-sections in Ningxia section.The results showed that the present six national monitoring cross-sections could not provide sufficient monitoring of the water quality variation of Yellow River，and two additional monitoring cross-sections should be added.One should be between Xiaheyan and Jinshawan national monitoring cross-sections，and the other should be between Yesheng Road Bridge and Yingu Road Bridge national monitoring cross-sections.They are respectively at Shikong county （105°40′31″E，37°34′04″N） and Linhe county （106°18′36″E，38°16′11″N）.

Key words： water quality monitoring；cross-section optimization；matter-element analysis；fuzzy clustering；WASP model

隨著現代社會經濟、工業的迅速發展，水環境保護已經成為社會高度關注的問題，其中江河作為重要的地表水資源，它的污染問題倍受關注。為了解我國河流水質的變化規律及污染物空間分布，需要對河流水質的監測斷面進行合理布置，使水質監測狀況更加合理有效，保證河流水環境的健康發展[1]。蔣艷君等[3]運用改進TOPSIS法對浦陽江干流水質監測斷面進行優化，將原來9個斷面優化為7個，證明此方法成立且計算結果準確有效。王輝等[4]運用改進的貼近度法對渾河干流水質監測斷面進行優化，其中引入重心距離算法對貼近度法進行改進，優化結果科學準確。蔣艷君等[5]也將改進的物元分析法應用到水質監測斷面的優化中，對浦陽江干流斷面進行優化。目前國內水系干流水質監測斷面優化的研究主要采用實測數據進行優化計算，由于個別監測斷面間距離過大，污染物在河流中遷移輸運規律復雜多變，不能了解其真實情況。故本文采用WASP模型對黃河寧夏段干流進行細分概化處理，對概化后的斷面進行優化研究，使結果更加準確合理，從而為寧夏水環境的監測、控制及水網建設提供更科學的保障，同時為寧夏的農業發展提供建設性的指導意義。

1 研究區概況

黃河寧夏段是從中衛縣南長灘開始，沿程至石嘴山市頭道坎，整個流程距離為397 km。由于工業廢水、生活污水和農田退水持續排入黃河，對沿河水環境造成極大污染。流域內河流TDS高、含沙量大，TDS一般為0.5～7.0 g/L，最高可達19 g/L，黃土丘陵區河流年均含沙量為100～380 kg/m3，年輸沙模數達到2 000～12 000 t/km2[2]。根據《2014年寧夏回族自治區環境狀況公報》，黃河干流寧夏段Ⅱ類水質斷面比例為66.7%，Ⅲ類水質斷面比例為33.3%，良好以上水質斷面達100%。

黃河寧夏段現有六個國控監測斷面，監控中衛、吳忠、銀川及石嘴山四市境內黃河水質狀況。其中中衛市境內黃河河道長度占總體35.8%，下河沿至金沙灣斷面間距較其他斷面間距離過大，監測數據不能反映其間水質真實情況；銀川市境內存在9條入黃排水溝（中干溝、永二干溝及第二排水溝等）及多個工業園區（寧東能源化工基地、望遠工業園區及銀川經濟開發區等），僅有一個監測斷面不能全面反映水質變化情況（見表1及圖1）。因此，對黃河寧夏段水質監測斷面的優化具有重要的實際意義。

2 研究方法與原理

2.1 模糊聚類法

聚類分析是采用數學方法判斷樣品之間的關系，客觀的對其進行分類，是數理統計中“物以類聚之”的多元分析方法[6-7]。由于水環境本身就是灰色系統，存在極大的模糊性，故將模糊數學的方法引入到聚類分析中，使分類更符合實際[8]。該方法在很多文獻中均有詳細介紹，這里不再贅述。

2.2 物元分析法

物元分析法是解決矛盾問題的方法，可將復雜問題抽象為形象化模型，并應用這些模型進行研究，提出相應的應用方法[9]。通過建立各點位物元模型和節域分析模型，計算各點位的線性關聯函數和綜合關聯函數，根據計算結果繪制點聚圖，根據圖中點位分布確定優化的點位。

2.2.1 確立各監測斷面物元模型

對全部監測斷面的污染指標進行對比，選出每項指標的最優值A、最差值B及期望值C，然后選擇各監測指標的最大值、最小值和平均值，用每一項指標的最優值A、最差值B及期望值C構造出2個物元矩陣：

3 WASP水質模型構建

基于黃河干流寧夏段地形特征、水環境特征及水動力學條件可知：黃河河道屬于寬淺型，水深遠小于河寬，且干流流量與污水流量之比遠大于10～20，可認為污染物短時間內在河道中混合均勻，其垂向和橫向濃度梯度變化可忽略不計，主要靠斷面縱向混合作用向前輸移，河流水力要素和水質狀況均可采用斷面平均值。綜合考慮以上水質模型適用性和可行性，選取WASP模型建立黃河干流水質模型。

3.1 模型概化

為保證研究的科學準確性，運用WASP水質模型將黃河干流寧夏段概化為兩部分，即衛寧段（WN）與青石段（QS）：衛寧段為下河沿至青銅峽段，劃分為27個斷面，各斷面間的距離約為4 km；青石段為青銅峽至石嘴山河段，劃分為44個斷面，各斷面間的距離約為4 km（其中將衛寧段最后一個斷面數據作為青石段初始條件輸入模型）。為增加模型精度，運用差值法對部分河段進行適當加密，分別計算出每段水體體積、相鄰分區間的特征距離及分區間剖面面積，作為模型輸入條件。圖2和圖3分別為衛寧段與青石段WASP模型概化后的劃分情況，其中由下河沿至金沙灣依次為WN1、WN2、…、WN27，青銅峽至石嘴山依次為QS1、QS2、…、QS44。

3.2 模型建立及參數率定

以2014年黃河干流水環境為基礎建立水質模型，模型選用EUTRO模塊，時間步長設為1，模擬指標為COD、氨氮、總磷，依次輸入河段初始濃度、邊界濃度、含氧量及水溫等相關參數。模型部分參數參考WASP用戶手冊進行設定，COD、氨氮、總磷自凈系數根據黃河上、中、下游水樣的室內降解實驗求得，分別為0.05/d、0.015/d和0.025/d，污染物縱向離散系數取50 m2/s，衛寧段與青石段河道綜合粗糙率分別為0.045和0.036。

3.3 相關性檢驗

通過MATLAB中corrcoef函數對實測數據與模擬數據進行相關性分析，表3為三種污染物的實測數據與模擬數據進行相關性分析后所得r值與p值，其中3種水質指標在6個監測斷面監測頻率均為每月一次，監測時間為上午9點至12點之內，采用水質自動監測技術。由表3可知：三種污染物實測數據與模擬數據的r值分別為0.978 9、0.993 1、0.986 9，均通過了95%的置信水平。圖4、圖5、圖6分別為COD、氨氮、總磷的實測數據與模擬數據對比，三幅圖中橫坐標序號（1至6）依次代表6個水質監測斷面，即下河沿、金沙灣、葉盛公路橋、銀古公路橋、平羅黃河大橋與麻黃溝監測斷面。由圖4-圖6可知：WASP水質模型的模擬結果與實際較符合，證明運用WASP水質模型對黃河干流的概化分段是科學的，能夠有效模擬河流水質的變化規律。

4 水質監測斷面優化過程

4.1 衛寧段斷面優化

4.1.1 衛寧段物元分析法

根據《地表水環境質量標準》（GB 3838-2002）中污染指標的標準值，運用指數超標法可得污染指標的權值Wj，經計算最后得歸一化權值的結果，表4為各污染指標的分級標準及歸一化權值，由表4可知：污染指標COD、氨氮、總磷的歸一化權值分別為0.100、0.425、0.475。將Wj代入式（7）、式（8），可得各斷面的綜合關聯函數，表5為代入公式后的各斷面綜合關聯函數。

在以KA為橫坐標、KB為縱坐標的平面上，繪制27個監測斷面綜合關聯函數的點聚圖，所有點均分布于第II、IV象限內。圖7（a）為第II象限內散點圖，由圖7（a）可知：在第II象限內KA>0，符合最佳斷面，斷面從左到右依次為WN1至WN10，其中WN2、WN10符合程度最高，選為最優斷面。圖7（b）為第IV象限內散點圖，由圖7（b）可知：在第IV象限內KB>0，符合最差斷面，斷面從左到右依次為WN11至WN27，其中WN11至WN15與WN24最為符合，是斷面中條件最優化的點，由于WN27有青銅峽斷面作為監測斷面，故可舍去。綜上所述，采用物元分析法將衛寧段優化為以下8個監測斷面：WN2、WN10、WN11、WN12、WN13、WN14、WN15、WN24。

4.1.2 衛寧段模糊聚類法

采用MATLAB進行模糊聚類分析，判斷各監測斷面的相關性。運用聚類函數分析法，對各監測斷面進行聚類分析，得到聚類樹圖（圖8），由圖8可知：衛寧段水質監測斷面總體上分為三部分。第一部分中的WN1至WN8八個斷面看作一類，WN9看作一類；第二部分中，WN10、WN11看作一類，WN12、作為一類，WN14至WN19看作一類；第三部分中，WN20至WN22、WN27作為一類，WN24至WN26劃為一類，WN23屬于單獨一類。4.1.3 綜合對比

表6為綜合物元分析法與模糊聚類法得到的優化結果，綜合兩種方法最終得到的優化后斷面為WN2、WN11、WN13、WN15、WN24。

4.2 青石段斷面優化

青石段水質監測斷面優化方法與衛寧段一致。表7為最后所得到的優化結果，綜合兩種方法最終得到的優化后斷面為QS5、QS7、QS10、QS25、QS30、QS44。4.3 優化結果

依據寧夏區域的地形地貌及沿河主要入黃排水溝、支流的分布情況，布置與優化監測斷面需考慮如下因素：（1）城鎮密集或工業園區流域內，需布設對照斷面、控制斷面和消減斷面；（2）污染嚴重河段可

根據主要入黃排污口分布及排污狀況，連續設置幾個控制斷面；（3）若河段內有較大支流匯入，須在交點上游與干流下游處設監測斷面；四、流入或流出行政區邊界處應設置監測斷面。

對于衛寧段，由于WN2距下河沿監測斷面過近，考慮去掉WN2，保留下河沿斷面；因WN11、WN13、WN15相鄰，且在它們之前有第一、第四排水溝及清水河匯入，考慮選擇其中一個斷面，鑒于監測的方便，最終選擇WN13，即石空鎮處設置監測斷面；WN24距金沙灣監測斷面過近，考慮去掉WN24，保留金沙灣斷面。對于青石段，由于QS5、QS7距葉盛公路橋斷面過近，考慮去掉它們，保留葉盛公路橋斷面；鑒于QS10處于臨河鎮上游，斷面之前聚集較多城鎮及羅家河排水溝的匯入，為了解此處水質受污染狀況，需在此設立監測斷面，即QS10（臨河鎮上游）處設置監測斷面；鑒于QS25、QS30距平羅黃河大橋較近，考慮去掉QS25、QS30，保留平羅黃河大橋斷面；因QS44距麻黃溝監測斷面較近，考慮去掉QS44，保留麻黃溝斷面。

綜合上述物元分析法、模糊聚類法、地理因素等，最終確定水質監測斷面為8個，較原來增設2個。表8為水質監測斷面最終優化結果，增設斷面分別位于石空鎮（靠近下河沿）與臨河鎮，其地理位置分別為105°40′31″E，37°34′04″N與106°18′36″E，38°16′11″N（見圖9）。對于增設的石空鎮監測斷面，可彌補因下河沿與金沙灣兩監測斷面之間距離過大，且兩斷面間有清水河及一系列入黃排水溝匯入，而造成對其間水質真實變化情況監測的空白；增設的臨河鎮監測斷面，可彌補因此段較多的城鎮及入黃排水溝匯入，而造成對其間水質真實變化情況掌握的空白。

5 結論

對黃河干流寧夏段國控水質監測斷面進行優化設置，可為寧夏水環境監測控制及水網建設提供建設性意見，且優化后監測斷面能更加科學合理的監測水質變化，為寧夏農業的發展提供保障。本文運用物元分析法、模糊聚類法及WASP水質模型對現有監測斷面進行優化布置，得以下結論。

（1）采用WASP水質模型對黃河寧夏段干流進行概化，將污染物的實測數據與模擬數據進行對比，得到其相關性檢驗的r值分別為0.978 9、0.993 1、0.986 9，均通過95%的置信水平，證明運用WASP水質模型對黃河干流寧夏段概化分段是合理有效的。

（2）為全面了解干流流域水質變化情況，監測斷面最終優化結果為：將原來6個國控水質監測斷面增加為8個，增設的水質監測斷面分別位于石空鎮與臨河鎮，其地理位置分別為105°40′31″E，37°34′04″N與106°18′36″E，38°16′11″N。

參考文獻（References）：

[1] 李瑞杰，王迪，張穎.河流水質監測控制斷面設置探討[J].內蒙古環境科學，2008，20（2）：41-42.（LI Rui-jie，WANG Di，ZHANG Ying.Approach on establishment of control section of river water quality monitoring[J].Inner Mongolian Environmental Sciences，2008，20（2）：41-42.（in Chinese））

[2] 寧夏回族自治區環境保護局.環境狀況公報[R].銀川：寧夏回族自治區環境保護局，2001-2003.（The environmental protection bureau of ningxia hui autonomous region.Environmental bulletin[R].Yinchuan：The environmental protection bureau of ningxia hui autonomous region，2001-2003.（in Chinese））

[3] 蔣艷君，謝悅波，黃旻.基于改進TOPSIS法的水質監測斷面優化研究[J].南水北調與水利科技，2016（5）：78-82.（JIANG Yan-jun，XIE Yue-bo，HUANG Min.Optimization of water quality monitoring sections based on improved TOPSIS method[J].South-to-North Water Transfers and Water Science & Technology，2016（5）：78-82.（in Chinese））

[4] 王輝，孫家君，孫麗娜，等.貼進度法的改進及在渾河干流水質監測斷面優化中應用[J].生態學雜志，2014（12）：3470-3474.（WANG Hui，SUN Jia-jun，SUN Li-na，et al.Modification of similarity method and its application in optimization design of river quality monitoring section in the mainstream of Hun River[J].Chinese Journal of Ecology，2014（12）：3470-3474.（in Chinese））

[5] 蔣艷君，謝悅波，黃旻.改進的物元分析法在水質監測斷面布設優化中的應用[J].水資源保護，2016，32（4）：136-141.（JIANG Yan-jun，XIE Yue-bo，HUANG Min.Application of improved matter element analysis method to optimization of water quality monitoring sections[J].Water Resources Protection，2016，32（4）：136-141.（in Chinese））

[6] 杜艷君，孫慶華，李湉湉.不同通勤方式下PM2.5個體暴露與固定站點監測的相關性研究[J].環境與健康雜志，2015（04）：319-323.（DU Yan-jun，SUN Qing-hua，LI Tian-tian.Relationship between personal exposure levels in different commuting patterns and corresponding air monitoring data of PM2.5[J].Journal of Environment and Health，2015（04）：319-323.（in Chinese））

[7] 安貝貝，蔣昌潭，劉蘭玉，等.水質監測斷面優化技術研究[J].環境科學與技術，2015（04）：107-111.（AN Bei-bei，JIANG Chang-tan，LIU Yu-lan，et al.Optimized approach to setting up water quality monitoring crosssections[J].Environmental Science & Technology，2015（04）：107-111.（in Chinese））

[8] 吳文強，陳求穩，李基明，等.江河水質監測斷面優化布設方法[J].環境科學學報，2010（08）：1537-1542.（WU Wen-qiang，CHEN Qiu-wen，LI Ji-ming，et al.Optimization of river water quality monitoring sections[J].Acta Scientiae Circumstantiae，2010（08）：1537-1542.（in Chinese））

[9] 袁正鋒.常見環境質量評價數學模型討論[J].生物技術世界，2012（02）：90-91.（YUAN Zheng-feng.Common environmental quality evaluation mathematical model for discussion[J].Biotech World，2012（02）：90-91.（in Chinese））

[10] 張荔，孫程，林金輝，等.WASP6水質模型在渭河流域水環境容量解析中的應用[J].水資源與水工程學報，2006，17（06）：12-14.（ZHANG Li，SUN Cheng，LIN Jin-hui，et al.Application of WASP6 model to water environmental capacity in Weihe River basin[J].Journal of Water Resources and Water Engineering，2006，17（06）：12-14.（in Chinese））

[11] 韋興荷.邕江南寧市區段水質時空變化規律分析[J].水資源保護，2005（06）：49-52.（WEI Xing-he.Analysis of the space-time variation rule of the water quality in Nanning section of the Yongjiang River[J].Water Resources Protection，2005（06）：49-52.（in Chinese））

[12] 王福偉.我國環境監測的現狀與發展[J].黑龍江科技信息，2016（24）.（WANG Fu-wei.Present situation and the development of our country′s environmental monitoring[J].Heilongjiang Science and Technology Information，2016（24）.（in Chinese））

[13] George Yu-zhu FU.Integration of USEPA WASP model in a GIS platform[J].Journal of Zhejiang University-Science A（Applied Physics & Engineering）.2010（12）.

[14] 左其亭，李可任.最嚴格水資源管理制度理論體系探討[J].南水北調與水利科技，2013，11（1）：34-38.（ZUO Qi-ting，LI Ke-ren.Discussion on theoretical system of the strictest water resources management system[J].South-to-North Water Transfers and Water Science & Technology，2013，11（1）：34-38.（in Chinese））

[15] 張志強，左其亭，馬軍霞.最嚴格水資源管理制度的和諧論解讀[J].南水北調與水利科技，2013，11（6）：133-137.（ZHANG Zhi-qiang，ZUO Qi-ting，MA Jun-xia.Harmony theory interpretation for the strictest water resources management syste[J].South-to-North Water tansfcrs and Water Scicnace &Tachnology，2013，11（6）：133-137.（in Chinese））

[16] 田靜宜，王新軍.基于熵權模糊物元模型的干旱區水資源承載力研究：以甘肅民勤縣為例[J].復旦學報：自然科學版，2013，52（1）：86-93.（TIAN Jing-yi，WANG Xin-jun.Assessment of the carrying capacity of water resources in arid areas based on an entropy fuzzy matter element model-A case study in Minqin county，Gansu province[J].Journal of Fudan University（Natural Science），2013，52（1）：86-93.（in Chinese））

[17] 王亞尼，周序協，張桂蓉，等.大茶湖浮游藻類調查與水質初步評價[J].華中農業大學學報，2013，32（03）：118-123.（WANG Ya-ni，ZHOU Xu-xie，ZHANG Gui-rong，et al.Investigation of phytoplankton and preliminary evaluationof water quality in Dacha Lake[J].Journal of Huazhong Agricultural University，2013，32（03）：118-123.（in Chinese））

[18] 富天乙，鄒志紅，王曉靜.基于多元統計和水質標識指數的遼陽太子河水質評價研究[J].環境科學學報，2014，34（02）：473-480.（FU Tian-yi，ZHOU Zhi-hong，WANG Xiao.Water quality assessment for Taizi River watershed in Liaoyang section based on multivariate statistical analysis and water quality identification index[J].Acta Scientiae Circumstantiae，2014，34（02）：473-480.（in Chinese））

[19] 崔芳.基于均值化的灰色關聯法在渭河寶雞段水質評價中的應用[J].中國農學通報，2016，32（25）：138-141.（CUI Fang.Water quality assessment of Weihe River Baoji section based on equalized grey correlation analysis[J].Chinese Agricultural Science Bulletin，2016，32（25）：138-141.（in Chinese））

[20] 王維，紀枚，蘇亞楠.水質評價研究進展及水質評價方法綜述[J].科技情報開發與經濟，2012，22（13）：129-131.（WANG Wei，JI Mei，SU Ya-nan.Reviews of the progress in the research of water quality evaluation and the methods for water quality evaluation[J].Ci-tech Information Development & Economy.2012，22（13）：129-131.（in Chinese））

[21] 章新，賀石磊，張雍照，等.水質評價的灰色關聯分析方法研究[J].水資源與水工程學報，2010，21（05）：117-119（ZHANG Xin，HE Shi-lei，ZHANG Yong-zhao，et al.Study on grey relational analysismethod for water quality assessment[J].Journal of Water Resouces & Water Engineering.2010，21（05）：117-119.（in Chinese））

[22] 弓曉峰，陳春麗，趙晉，等.鄱陽湖樂安河流域水質監測優化布點[J].湖泊科學，2006，18（05）：545-549.（GONG Xiao-feng，CHEN Chun-li，ZHAO Jin，et al.Optimized points selection of the water quality monitoring in the Le′an River basin of Lake Poyang[J].J.Lake Sci，2006，18（05）：545-549.（in Chinese））

[23] 翁建中，李繼影，梁柱，等.太湖藍藻水華時空分布與預警監測響應的分析[J].環境監控與預警，2010，02（03）：1-4.（WENG Jian-zhong，LI Jiying，LIANG Zhu，et al.Spatial and tanporal d istribution of cyanobacteria bloom in Taihu Lake and analysis of the response to early warning monitoring[J].Envirormental Monitoring and Forewarning，2010，02（03）：1-4.（in Chinese））

[24] 李茜，張建輝，羅海江.區域地表水監測斷面優化調整的技術方法[J].環境監控與預警，2012，04（01）：34-37.（LI Qian，ZHANG Jian-hui，LUO Hai-jiang.The study on the optimization and adjustment of regional surface water environment monitoring sections[J].Environmental Monitoring and Forewarning，2012，04（01）：34-37.（in Chinese））