廣東省入海河流的水質(zhì)評價實證分析案例目錄TOC\o"1-3"\h\u30722廣東省入海河流的水質(zhì)評價實證分析案例 111236第一章課題來源 129949第二章廣東省入海河流基本概括 1412.1廣東省基本概況 1115672.1.1地理位置 1298742.1.2氣候特征 2173672.1.3地貌地形 3289392.1.4區(qū)域社會經(jīng)濟概括 4139442.2廣東省入海河流基本概況 5266942.2.1水文地質(zhì)條件 5207452.2.2水資源概括 5202012.2.3水功能區(qū) 521694第三章廣東省入海河流水質(zhì)評價 733473.1水質(zhì)單因子評價分析 7281713.1.1水質(zhì)評價方法—單因子評價法 7140623.1.2水質(zhì)監(jiān)測斷面及指標(biāo)的選取 8153123.1.3監(jiān)測斷面單因子評價指數(shù)計算 9260393.1.4單因子評價指數(shù)計算結(jié)果分析 10309263.2水質(zhì)綜合評價分析 1276353.2.1水質(zhì)評價方法——灰色關(guān)聯(lián)評價法 1286603.2.2灰色關(guān)聯(lián)評價法計算過程 12212023.2.3灰色關(guān)聯(lián)評價結(jié)果分析 21課題來源任務(wù)書廣東省入海河流基本概括2.1廣東省基本概況2.1.1地理位置廣東省地處我國大陸最南端，自東向西依次與福建省、江西省、湖南省、廣西壯族自治區(qū)接壤；毗鄰香港、澳門特別行政區(qū)；西南端隔瓊州海峽與海南省相望。陸地最東端至饒平縣大埕鎮(zhèn)，最西端至廉江市高橋鎮(zhèn)，廣東省最東端至最西段最大跨度800千米。最北端為樂昌市白石鎮(zhèn)，最南端至徐聞縣角尾鎮(zhèn)，南北跨度約為600千米。根據(jù)2018年土地變更調(diào)查統(tǒng)計數(shù)據(jù)，全省土地總面積為1797.25萬公頃，約占全國陸地面積的1.87%。海域面積42萬平方千米，是陸域面積的2.3倍。大陸海岸線長4114千米，居全國首位。有海島1963個，總面積1513.17平方千米，在全國沿海?。▍^(qū)、市）中位列第二，僅次于浙江省。圖2-1廣東省政區(qū)圖2.1.2氣候特征廣東省從北向南共有三種氣候類型，分別為中亞熱帶氣候、南亞熱帶氣候和熱帶氣候，全省光、熱、水資源均較為豐沛，該地區(qū)降水時間主要集中在4—9月份并且雨熱同期。全省年平均氣溫21.8℃。年平均氣溫分布呈南高北低特征，年平均氣溫最高23.8℃，最低18.9℃。月平均氣溫最低出現(xiàn)在1月為13.3℃，最高出現(xiàn)在7月為28.5℃。年平均降水量為1789mm，降水量最少為1314.1mm，最多達2254.1mm。并且年降水年量規(guī)律為分布不均，呈多中心分布。恩平年平均降水量超過2500mm、海豐接近2500mm和龍門為2100mm，恩平、海豐、龍門為3個多雨中心。海豐地區(qū)暴雨頻率最高，年平均暴雨周期達13.5天。月平均降水量最少出現(xiàn)在12月份，最多出現(xiàn)在6月份。最大日降水量640.6mm，出現(xiàn)在清遠。廣東省自北向南年平均日照時數(shù)逐漸增加，平均日照時數(shù)最少為不足1500小時，最多高于2300小時。由于廣東省地理位置特殊，各種氣象災(zāi)害頻繁發(fā)生，主要災(zāi)害有雷擊、高溫、洪澇災(zāi)害、熱帶氣旋、強對流天氣、干旱及低溫陰雨、寒露風(fēng)、寒潮和冰（霜）凍等低溫災(zāi)害，呈現(xiàn)災(zāi)害種類多，災(zāi)害時間長，災(zāi)害頻率高，災(zāi)害受損程度嚴重的狀態(tài)[35]。2.1.3地貌地形廣東省由于受到地殼運動、褶皺和斷裂構(gòu)造以及外力作用等諸多因素的綜合影響，導(dǎo)致其地貌類型復(fù)雜多樣，山地面積占全省土地總面積的33.7%、丘陵占全省土地總面積的24.9%、臺地占全省土地總面積的14.2%、平原占全省土地總面積的21.7%，河流和湖泊等只占全省土地總面積的5.5%。廣東省地勢總體呈現(xiàn)南低北高的規(guī)律，南部多為平原和臺地，而北部多為山地和丘陵。全省山脈走向大多是北東―南西，山脈之間分布有有大小谷地和盆地。珠江三角洲是廣東省面積最大的平原，潮汕是廣東省面積第二的平原，此外還有高要、清遠等諸多的沖積平原。臺地以雷州半島―電白―陽江一帶和海豐―潮陽一帶分布較多。在廣東省地貌類型中花崗巖較為普遍是構(gòu)成各類地貌的基巖巖石，而砂巖和變質(zhì)巖也較多，另外部分地區(qū)還有紅色巖系地貌構(gòu)成了奇特景觀，如丹霞山和金雞嶺等。并且廣東省擁有優(yōu)質(zhì)海灣數(shù)量眾多，在雷州半島西南岸的珊瑚礁也成為了非常重要的地貌旅游資源。圖2-2廣東省地勢圖2.1.4區(qū)域社會經(jīng)濟概括經(jīng)國家統(tǒng)計局統(tǒng)一核算，2020年廣東實現(xiàn)地區(qū)生產(chǎn)總值（初步核算數(shù)）110760.94億元，同比增長2.3%。其中，第一產(chǎn)業(yè)增加值4769.99億元，同比增長3.8%，對地區(qū)生產(chǎn)總值增長的貢獻率為6.4%；第二產(chǎn)業(yè)增加值43450.17億元，同比增長1.8%，對地區(qū)生產(chǎn)總值增長的貢獻率為33.7%；第三產(chǎn)業(yè)增加值62540.78億元，同比增長2.5%，對地區(qū)生產(chǎn)總值增長的貢獻率為59.9%。第一、第二、第三產(chǎn)業(yè)結(jié)構(gòu)比重為4.3：39.2：56.5，第三產(chǎn)業(yè)所占比重同比提高0.7個百分點。新經(jīng)濟增加值27862.23億元，增長3.0%，占地區(qū)生產(chǎn)總值的25.2%[36]。圖2-3廣東省地區(qū)生產(chǎn)總值2.2廣東省入海河流基本概況2.2.1水文地質(zhì)條件廣東省地處珠江流域下游，境內(nèi)河流眾多，水系發(fā)達，集雨面積大于50平方千米的河流有1211條，河流總長3.67萬千米，主要江河有東江、西江、北江、韓江、鑒江等，海岸線長4114千米，為全國最長。2018年，全省年平均降水量為1843.1mm，較2017年增多6.0%，較常年增多4.1%，屬平水年。全省降水呈現(xiàn)時空分布不均的特點。在空間分布上，沿海和珠江三角洲部分地區(qū)年降水量大于內(nèi)陸地區(qū)。在時間分布上，汛期（4—9月）降水量1542.2mm；前汛期（4—6月）降水量647.6mm；后汛期（7—9月）降水量894.6mm。從鄰省流入省內(nèi)的總?cè)刖乘?022.9億立方米，從省內(nèi)流出鄰近省份的水量24.5億立方米，全省入海水量3746.9億立方米。2.2.2水資源概括廣東省境內(nèi)河流眾多，以珠江流域、韓江流域、粵東諸河和粵西諸河為主，集水面積占全省面積的99.8%，其余屬于長江流域的鄱陽湖和洞庭湖水系。全省流域面積在100平方千米以上的各級干支流共有614條。獨流入海河流52條，較大的有韓江、榕江、漠陽江、鑒江、九洲江等。全省多年平均水資源總量1830億立方米，其中地表水資源量1820億立方米，地下水資源量450億立方米。除省內(nèi)降水外，還有來自珠江、韓江等上游從鄰省入境水量2361億立方米。廣東省水資源呈現(xiàn)時空分布不均的特點，夏秋水量多，冬春較為干旱。由于沿海臺地和低丘陵區(qū)較難蓄水，所以缺水現(xiàn)象尤其突出。部分河流中下游河段由于工業(yè)廢水和生活污水的排放造成污染，存在水質(zhì)性缺水問題。2018年，全省地表水資源量1885.2億立方米，地下水資源量460.6億立方米，水資源總量1895.1億立方米。全省38座大型水庫和339座中型水庫年末蓄水總量168.0億立方米，較年初減少32.8億立方米。全省用水總量為421億立方米；萬元國內(nèi)生產(chǎn)總值用水量43.3立方米，萬元工業(yè)增加值用水量26.4立方米。2.2.3水功能區(qū)2018年廣東省全省劃定水功能區(qū)3520個，廣東省水功能區(qū)劃定采用了兩級體系。一級功能區(qū)劃分的依據(jù)是從宏觀上解決水資源開發(fā)利用與保護的問題，主要目的為協(xié)調(diào)地區(qū)間用水關(guān)系；二級區(qū)劃分的依據(jù)主要是為了協(xié)調(diào)各市和市內(nèi)用水部門之間的關(guān)系。一級功能區(qū)分為四類，包括保護區(qū)、保留區(qū)、開發(fā)利用區(qū)、緩沖區(qū)；二級功能區(qū)分重點在一級區(qū)所劃的開發(fā)利用區(qū)內(nèi)進行，分七類，包括飲用水源區(qū)、工業(yè)用水區(qū)、農(nóng)業(yè)用水區(qū)、漁業(yè)用水區(qū)、景觀娛樂用水區(qū)、過渡區(qū)、排污控制區(qū)[37]。

廣東省入海河流水質(zhì)評價3.1水質(zhì)單因子評價分析3.1.1水質(zhì)評價方法—單因子評價法單因子評價法是最常見的水質(zhì)評價方法，該方法將單項指標(biāo)作為判斷該水體水質(zhì)類別的判斷依據(jù)，通過水質(zhì)實測值與標(biāo)準值的比值來判定。單因子評價方法的優(yōu)點在于其計算簡便、能以較快的速度判斷目標(biāo)水質(zhì)的水質(zhì)類別,應(yīng)用起來簡單方便。但僅使用單因子評價方法評價水質(zhì)過于局限,影響水體質(zhì)量的因素眾多，往往是各種水質(zhì)指標(biāo)綜合作用的結(jié)果，使用單因子評價法來進行水質(zhì)評價未考慮到水質(zhì)指標(biāo)的綜合影響不夠全面缺乏可信性和科學(xué)性[38]。單因子評價方法具體如下所示:單因子評價指數(shù)越大，表明評價因子對目標(biāo)評價水體的污染越大，反之則污染越小。單因子評價指數(shù)公式為：（3.1）（3.1）式中：特別的，對于溶解氧（DO）指標(biāo)，其單因子評價計算公式為：（3.2）（3.3）（3.2）（3.3）式中：3.1.2水質(zhì)監(jiān)測斷面及指標(biāo)的選?。?）廣東省入海河流監(jiān)測斷面的選擇據(jù)數(shù)據(jù)顯示，廣東省入海河流共有52條，但由于數(shù)據(jù)缺乏并且有些河流較小，故本文針對廣東省設(shè)置的27個國考監(jiān)測斷面進行研究，選擇全部27個監(jiān)測斷面進行研究，各監(jiān)測斷面名稱及地理位置如表3-1所示。表3-1廣東省入海河流監(jiān)測斷面概況序號監(jiān)測斷面名稱河流名稱所在城市1鳳江橋黃岡河潮州市2沙田泗盛東江南支流東莞市3蓮花山珠江廣州河段廣州市4洪奇瀝珠江廣州河段廣州市5蕉門珠江廣州河段廣州市6虎爪斷橋淡澳河惠州市7蒼山渡口潭江江門市8寨頭河出海口寨頭河茂名市9蓮陽橋閘韓江東溪汕頭市10外砂韓江外砂河汕頭市11升平韓江梅溪河汕頭市12地都榕江汕頭市13海門灣橋閘練江汕頭市14東溪水閘黃江河汕尾市15海豐西閘黃江河汕尾市16陸豐半灣水閘螺河汕尾市17烏坎烏坎河汕尾市18深圳河口深圳河深圳市19埠場漠陽江陽江市20尖山漠陽江陽江市21那格漠陽江陽江市22黃坡鑒江湛江市23大山江袂花江湛江市序號監(jiān)測斷面名稱河流名稱所在城市24營仔九洲江湛江市25中山港碼頭橫門水道中山市26雞啼門大橋雞啼門水道珠海市27珠海大橋磨刀門水道珠海市（2）水質(zhì)評價指標(biāo)的選取綜合考慮廣東省內(nèi)水質(zhì)污染情況以及以往各類水質(zhì)評價研究，最終確定選取5項水質(zhì)評價指標(biāo)如下表3-2所示。表3-2本研究選取的水質(zhì)評價指標(biāo)序號指標(biāo)1溶解氧（DO）2高錳酸鹽指數(shù)3生化需氧量（BOD）4氨氮（NH3-N）5總磷（TP）3.1.3監(jiān)測斷面單因子評價指數(shù)計算在廣東省入海河流27個監(jiān)測斷面上選取溶解氧（DO）、高錳酸鹽指數(shù)、生化需氧量（BOD）、氨氮（NH3-N）、總磷（TP）５個指標(biāo)作為主要水質(zhì)污染因子指標(biāo)，依據(jù)《地表水環(huán)境質(zhì)量標(biāo)準》（GB3838-2002）中Ⅲ類水水質(zhì)標(biāo)準，對廣東省入海河流進行水質(zhì)評價，并判定監(jiān)測斷面的水質(zhì)類別和主要污染指標(biāo)。表3-3地表水環(huán)境質(zhì)量標(biāo)準基本項目標(biāo)準限值單位：mg/L序號項目Ⅰ類Ⅱ類Ⅲ類Ⅳ類Ⅴ類1溶解氧（DO）飽和率90%（或7.5）65322高錳酸鹽指數(shù)24610153生化需氧量（BOD）3346104氨氮（NH3-N）0.150.511.525總磷（TP）0.020.10.20.30.4單因子評價指數(shù)計算結(jié)果見附錄53.1.4單因子評價指數(shù)計算結(jié)果分析對廣東省入海河流2019-2020年27個監(jiān)測斷面的五項水質(zhì)指標(biāo)初步分析，溶解氧（DO）、高錳酸鹽指數(shù)、生化需氧量（BOD）、氨氮（NH3-N）、總磷（TP）均有未達Ⅲ類水體標(biāo)準情況，2019年全省27個入海河流水質(zhì)類型Ⅰ-Ⅲ為66.7%（18個），14.8%（4個）為Ⅳ類，7.4%（2個）為Ⅴ類，11.1%（3個）為劣Ⅴ類。2020年全省27個入海河流斷面水質(zhì)類型Ⅰ-Ⅲ為74.1%，Ⅳ類占22.2%，消除劣Ⅴ類水質(zhì)類型。圖3-12019、2020廣東省入海河流水質(zhì)類別比例其中虎爪斷橋、寨頭河出?？?、海門灣橋閘三個監(jiān)測斷面為劣Ⅴ類水質(zhì)類型，水質(zhì)情況較差。虎爪斷橋監(jiān)測斷面主要污染指標(biāo)為氨氮、總磷、生化需氧量、溶解氧；寨頭河出?？诒O(jiān)測斷面主要污染指標(biāo)為氨氮、總磷、高錳酸鹽指數(shù)、溶解氧、生化需氧量；海門灣橋閘監(jiān)測斷面主要污染指標(biāo)為氨氮、總磷、高錳酸鹽指數(shù)、溶解氧、生化需氧量。虎爪斷橋監(jiān)測斷面2019年4月及7月氨氮指標(biāo)處于較大值，2019年2月總磷出現(xiàn)最大值，生化需氧量和溶解氧均處于較高水平，到2020年虎爪斷橋監(jiān)測斷面水質(zhì)有所改善，已基本消除劣Ⅴ類?；⒆鄻虮O(jiān)測斷面位于廣東省惠州市淡澳河，根據(jù)相關(guān)資料顯示淡澳河河流域污水處理設(shè)施、配套管網(wǎng)設(shè)施較為落后，重污染企業(yè)數(shù)量較多，養(yǎng)殖場數(shù)量較多，故該斷面常年處于劣Ⅴ類水體，由于惠州市開展淡澳河流域水污染防治攻堅戰(zhàn)，加大環(huán)境違法案件查處力度、清理整治非法養(yǎng)殖場、擴建污水處理廠等一系列措施，使得虎爪斷橋監(jiān)測斷面水質(zhì)改善明顯。寨頭河出海口監(jiān)測斷面氨氮指標(biāo)最大值、總磷指標(biāo)最大值、生化需氧量最大值和高錳酸鹽指數(shù)最大值均出現(xiàn)在2019年12月；溶解氧最大值出現(xiàn)在2019年4月。寨頭河出?？诒O(jiān)測斷面位于廣東省茂名市寨頭河，根據(jù)相關(guān)資料顯示寨頭河周邊以農(nóng)業(yè)耕地、村莊為主，并有數(shù)個大型樓盤和城郊住宅區(qū)。河道主要污染源為農(nóng)村的生活污水和養(yǎng)殖廢水；部分河段受到磚廠以及混凝土攪拌廠等生活污水的污染；河道內(nèi)垃圾多、淤積嚴重、駁岸松散、自凈能力差，故常年為劣Ⅴ類水體。由于我國開展黑臭水體整治工作，寨頭河被列入重點整治問題，經(jīng)過河底清淤、設(shè)置生態(tài)圍隔、河床底泥原位消減、曝氣增氧等一系列措施，到2020年寨頭河出海口監(jiān)測斷面水質(zhì)達到地表水Ⅴ類。海門灣橋閘監(jiān)測斷面溶解氧指標(biāo)最大值出現(xiàn)在2019年5月；高錳酸鹽指數(shù)最大值出現(xiàn)在2019年11月；生化需氧量指標(biāo)最大值出現(xiàn)在2019年4月；氨氮指標(biāo)最大值出現(xiàn)在2019年3月；總磷指標(biāo)最大值出現(xiàn)在2019年2月。海門灣橋閘監(jiān)測斷面位于汕頭市練江，根據(jù)相關(guān)資料顯示練江被廣東省環(huán)境保護廳定性為“全省污染最嚴重的河流”，練江天然水源不足、自凈能力較弱，紡織制衣是當(dāng)?shù)氐闹еa(chǎn)業(yè)，印染企業(yè)排放的工業(yè)廢水成為練江主要污染源之一，并且練江流域生活污水處理能力不足、配套管網(wǎng)不完善、污水處理率低，導(dǎo)致練江常年處于劣Ⅴ類水體。由于汕頭市練江流域綜合整治項目開展，不斷推進印染紡織環(huán)保綜合處理中心、污水處理設(shè)施及配套管網(wǎng)和流域環(huán)境整治等重點項目，到2020年練江流域水質(zhì)得到改善，海門灣橋閘監(jiān)測斷面水質(zhì)達到地表水Ⅴ類。3.2水質(zhì)綜合評價分析3.2.1水質(zhì)評價方法——灰色關(guān)聯(lián)評價法由于水環(huán)境質(zhì)量數(shù)據(jù)監(jiān)測的時間以及空間都存在一定的局限性，信監(jiān)測信息不夠準確，因此可將水環(huán)境系統(tǒng)視為一個灰色系統(tǒng)，就是部分信息已知，部分信息未知的系統(tǒng)[39]，由于單因子評價只是根據(jù)某一水質(zhì)指標(biāo)進行水質(zhì)評級，缺乏綜合評價，而灰色關(guān)聯(lián)評價法是根據(jù)待評價樣本的水質(zhì)與水環(huán)境質(zhì)量評價標(biāo)準在整體上的關(guān)聯(lián)度從而判定水質(zhì)類型，因此能減少一些不確定因素，較好地為區(qū)域水環(huán)境指標(biāo)進行綜合評價分級，據(jù)此對水環(huán)境進行綜合評價，具體步驟如下：（1）確定分析序列。分析數(shù)列分為參考數(shù)列和比較數(shù)列，本文選擇水體實測值作為參考數(shù)列X0k（k=1,2,3，……n），水體質(zhì)量的分級標(biāo)準為比較數(shù)列Xik（k=1,2,3，……n）。（2）數(shù)據(jù)預(yù)處理。將樣本數(shù)據(jù)進行無量綱處理，在水質(zhì)評價中使用的無量綱化方法主要有、均值化變換、初值化變換，本文采用標(biāo)準化變換。（3）求關(guān)聯(lián)系數(shù)。計算公式如下：（3.4（3.4）計算關(guān)聯(lián)度。計算公式如下：（3.5）（3.5）根據(jù)關(guān)聯(lián)度的大小進行排序。Ri的值越大，表示參考序列與比較序列的關(guān)聯(lián)度越高，關(guān)聯(lián)度最大則可判定該水體水質(zhì)類別[40]。3.2.2灰色關(guān)聯(lián)評價法計算過程1、確定分析序列本文根據(jù)廣東省生態(tài)環(huán)境廳發(fā)布的2019年和2020年廣東省入海河流27個斷面監(jiān)測數(shù)據(jù)年平均值為參考數(shù)列X0k（k=1,2,3，……n），根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準》（GB3838-2002），確定5級水質(zhì)分級標(biāo)準為標(biāo)較數(shù)列Xik（k=1,2,3，……n）。表3-42019年廣東省入海河流監(jiān)測數(shù)據(jù)年均值單位：mg/L序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷1鳳江橋7.042.982.970.360.122沙田泗盛3.733.283.040.760.123蓮花山4.283.182.090.440.144洪奇瀝7.051.800.720.120.085蕉門6.892.540.730.120.106虎爪斷橋5.125.183.685.820.397蒼山渡口6.322.501.070.080.078寨頭河出?？?.256.965.853.381.199蓮陽橋閘8.292.021.240.070.0510外砂8.361.631.230.040.0511升平6.432.331.640.520.1012地都5.003.382.880.320.0813海門灣橋閘7.445.736.833.040.2214東溪水閘7.645.683.120.480.0915海豐西閘8.164.293.130.790.1316陸豐半灣水閘8.002.111.530.150.0517烏坎6.383.341.230.230.07序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷18深圳河口3.565.122.421.570.2319埠場5.992.631.930.640.2020尖山6.393.032.640.830.1821那格5.672.302.030.650.1722黃坡6.754.412.150.170.1223大山江6.354.382.730.410.1624營仔5.425.293.110.260.1125中山港碼頭7.201.771.180.310.0826雞啼門大橋7.461.930.900.050.0627珠海大橋7.301.890.620.050.07表3-52020年廣東省入海河流監(jiān)測數(shù)據(jù)年均值單位：mg/L序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷1鳳江橋7.053.052.390.320.102沙田泗盛4.793.671.680.350.093蓮花山5.663.181.130.110.094洪奇瀝7.631.750.630.100.065蕉門7.372.450.820.100.066虎爪斷橋5.734.712.411.170.197蒼山渡口6.751.880.850.080.108寨頭河出?？?.955.532.991.380.249蓮陽橋閘8.031.741.120.050.0210外砂7.631.731.100.050.0311升平6.662.791.680.610.1212地都4.913.691.980.140.0713海門灣橋閘6.867.274.641.390.1514東溪水閘8.607.233.410.280.0715海豐西閘8.604.572.630.310.0716陸豐半灣水閘8.131.801.580.060.0317烏坎5.973.641.450.300.0618深圳河口4.955.313.020.810.1819埠場7.113.042.400.510.17序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷20尖山6.243.432.730.740.1421那格5.902.742.500.530.1622黃坡7.564.982.630.090.0723大山江6.155.163.420.440.1624營仔7.257.133.560.220.1025中山港碼頭6.911.480.980.300.0726雞啼門大橋7.281.870.940.090.0727珠海大橋7.441.840.570.060.06表3-6地表水環(huán)境質(zhì)量標(biāo)準（GB3838-2002）單位：mg/L水質(zhì)級別溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷Ⅰ7.5230.150.02Ⅱ6430.50.1Ⅲ56410.2Ⅳ31061.50.3Ⅴ2151020.42、數(shù)據(jù)預(yù)處理由于各個指標(biāo)的量級不同，必須在進行關(guān)聯(lián)度計算之前對各指標(biāo)的原始數(shù)據(jù)進行預(yù)處理，本文將采用標(biāo)準化變化，分以下兩種情況：當(dāng)評價指標(biāo)為正向指標(biāo)時（3.6）（3.6）當(dāng)評價指標(biāo)為負向指標(biāo)時（3.7）（3.7）將原始數(shù)據(jù)帶入以上公式，處理后得到如下數(shù)據(jù)，見表3-7：表3-72019年廣東省入海河流監(jiān)測數(shù)據(jù)無量綱值序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷1鳳江橋0.720.750.620.940.942沙田泗盛0.030.690.610.880.943蓮花山0.150.710.760.930.924洪奇瀝0.730.970.980.990.985蕉門0.690.830.980.990.966虎爪斷橋0.330.330.510.000.707蒼山渡口0.580.840.930.990.988寨頭河出海口0.140.000.160.420.009蓮陽橋閘0.990.930.900.991.0010外砂1.001.000.901.001.0011升平0.600.870.840.920.9612地都0.300.670.640.950.9813海門灣橋閘0.810.230.000.480.8514東溪水閘0.850.240.600.920.9715海豐西閘0.960.500.600.870.9316陸豐半灣水閘0.920.910.850.981.0017烏坎0.590.680.900.970.9818深圳河口0.000.350.710.740.8419埠場0.510.810.790.900.8720尖山0.590.740.670.860.8921那格0.440.870.770.900.8922黃坡0.670.480.750.980.9423大山江0.580.480.660.940.9024營仔0.390.310.600.960.9525中山港碼頭0.760.970.910.950.9826雞啼門大橋0.810.940.951.000.9927珠海大橋0.780.951.001.000.98表3-82020年廣東省入海河流監(jiān)測數(shù)據(jù)無量綱值序號監(jiān)測斷面名稱溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷1鳳江橋0.590.730.550.800.642沙田泗盛0.000.620.730.770.683蓮花山0.230.710.860.950.704洪奇瀝0.750.950.990.960.825蕉門0.680.830.940.960.836虎爪斷橋0.250.440.550.170.237蒼山渡口0.510.930.930.980.658寨頭河出?？?.300.300.400.010.009蓮陽橋閘0.850.950.871.000.9810外砂0.750.960.871.000.9711升平0.490.770.730.580.5512地都0.030.620.650.930.7613海門灣橋閘0.540.000.000.000.4014東溪水閘1.000.010.300.830.7915海豐西閘1.000.470.500.810.7616陸豐半灣水閘0.880.940.751.000.9817烏坎0.310.630.780.810.8118深圳河口0.040.340.400.430.2919埠場0.610.730.550.660.3220尖山0.380.660.470.480.4521那格0.290.780.530.640.3722黃坡0.730.390.500.970.7623大山江0.360.360.300.710.3624營仔0.650.020.270.870.6525中山港碼頭0.561.000.900.810.7826雞啼門大橋0.650.930.910.970.7727珠海大橋0.700.941.001.000.82表3-9地表水環(huán)境質(zhì)量標(biāo)準無量綱值水質(zhì)級別溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷Ⅰ1.001.001.001.001.00Ⅱ0.730.851.000.810.79Ⅲ0.550.690.860.540.53Ⅳ0.180.380.570.270.26Ⅴ0.000.000.000.000.003、求關(guān)聯(lián)系數(shù)根據(jù)表3-7、表3-8和表3-9的數(shù)據(jù)，監(jiān)測數(shù)據(jù)組成一個參考數(shù)列，水環(huán)境質(zhì)量的評價標(biāo)準中的5個水質(zhì)分級組成一個比較數(shù)列，然后將參考數(shù)列逐個與比較數(shù)列進行組合，求出絕對差Δi（k）見表4，依次求出2019年和2020年廣東省入海河流27個監(jiān)測斷面監(jiān)測數(shù)據(jù)的minminΔi（k）和maxmaxΔi（k）。根據(jù)公式3-4求出關(guān)聯(lián)系數(shù)ξi（k），根據(jù)公式3-5算出關(guān)聯(lián)度。由于各個年份對應(yīng)的影響因子數(shù)據(jù)不同，首先計算2019年鳳江橋斷面與比較數(shù)列的絕對值差，根據(jù)公式，以鳳江橋斷面：溶解氧=0.72，高錳酸鹽指數(shù)=0.75，生化需氧量=0.62，氨氮=0.94，總磷=0.94為特征序列，其余序列是以水環(huán)境質(zhì)量標(biāo)準等級(Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ)中各個因子溶解氧、高錳酸鹽指數(shù)、生化需氧量、氨氮、總磷5個指標(biāo)不同水質(zhì)標(biāo)準值進行灰色關(guān)聯(lián)分析，計算過程如下:絕對值差序列取兩級最小差和最大差結(jié)合上述計算結(jié)果帶入公式3-4，取分辨系數(shù)ζ=0.5得：把計算出的差序列依次帶入上述公式中得出2019年鳳江橋斷面待評價因子與水質(zhì)標(biāo)準等級的關(guān)聯(lián)函數(shù)，見表3-10：表3-102019年鳳江橋斷面待評價因子與比較序列的關(guān)聯(lián)系數(shù)關(guān)聯(lián)系數(shù)溶解氧高錳酸鹽指數(shù)生化需氧量氨氮總磷ξ10.640.670.560.910.91ξ21.000.840.560.800.77ξ30.750.910.680.550.55ξ40.480.570.920.420.42ξ50.400.390.440.340.344、求關(guān)聯(lián)度根據(jù)以上求出的關(guān)聯(lián)系數(shù)帶入公式得2019年鳳江橋斷面待評價樣本對各級別的關(guān)聯(lián)度：R1=0.74，R2=0.8，R3=0.69，R4=0.56，R5=0.38，根據(jù)關(guān)聯(lián)度最大原則，R2最大，2019年鳳江橋斷面水環(huán)境質(zhì)量等級為Ⅱ類。根據(jù)以上步驟分別求出2019年和2020年廣東省入海河流27個斷面的關(guān)聯(lián)系數(shù)及關(guān)聯(lián)度，如下表3-11：表3-112019年廣東省入海河流監(jiān)測斷面待評價樣本對各級別的關(guān)聯(lián)度序號監(jiān)測斷面名稱R1R2R3R4R51鳳江橋0.740.800.690.560.382沙田泗盛0.640.670.650.630.503蓮花山0.700.720.710.640.464洪奇瀝0.900.840.640.460.355蕉門0.860.890.700.490.376虎爪斷橋0.460.560.620.750.627蒼山渡口0.970.730.530.420.348寨頭河出?？?.370.430.540.690.829蓮陽橋閘0.940.760.640.460.3510外砂0.970.740.620.450.3411升平0.800.850.770.520.3812地都0.720.700.700.650.4513海門灣橋閘0.540.630.610.570.5914東溪水閘0.740.700.600.630.45序號監(jiān)測斷面名稱R1R2R3R4R515海豐西閘0.760.730.640.630.4116陸豐半灣水閘0.890.770.660.470.3517烏坎0.800.790.790.530.3918深圳河口0.560.660.640.690.5619埠場0.740.850.800.560.4020尖山0.720.850.800.600.4121那格0.750.820.740.560.4022黃坡0.750.760.710.610.4223大山江0.710.760.740.650.4424營仔0.680.660.640