甘文霞GAN Wen-xia；潘俊杰PAN Jun-jie；肖安娜XIAO An-na；楊如琴YANG Ru-qin

（①武漢工程大學(xué)土木工程與建筑學(xué)院，武漢 430205；②湖北交投智能檢測(cè)股份有限公司，武漢 430050；③武漢市自然資源和規(guī)劃信息中心，武漢 430014）

0 引言

武漢市作為百湖之市﹐其湖泊水環(huán)境直接影響著周邊居民生活、城市生態(tài)與經(jīng)濟(jì)和諧發(fā)展。近年來(lái)，隨著武漢市水環(huán)境綜合治理規(guī)劃實(shí)施，已基本實(shí)現(xiàn)“污水進(jìn)管網(wǎng)、雨水入湖塘”目標(biāo)，當(dāng)前的城市內(nèi)部湖泊水體污染主要來(lái)源為雨水徑流[1]，對(duì)降雨事件前后湖泊水質(zhì)的影響因素及其特點(diǎn)進(jìn)行分析具有重要的理論意義和實(shí)用價(jià)值。

學(xué)者們對(duì)水環(huán)境的影響因素研究已經(jīng)做了很多工作，其中針對(duì)土地利用因素對(duì)雨水徑流問(wèn)題的影響研究較多[2，3]，研究結(jié)果表明不同土地利用的差異會(huì)導(dǎo)致水質(zhì)濃度有規(guī)律地差異分布。伴隨著道路建設(shè)蓬勃發(fā)展，交通活動(dòng)日趨頻繁，由路面徑流沖刷污染物進(jìn)入受納水體而產(chǎn)生污染也引起了人們重視[4，5]。但已有研究主要從水質(zhì)的整體狀態(tài)展開(kāi)分析，將研究對(duì)象視為整體對(duì)其水質(zhì)或污染源進(jìn)行分析[6]，但當(dāng)水體面積較大時(shí)，其水質(zhì)分布及其影響特點(diǎn)具有空間異質(zhì)的特點(diǎn)，整體性的分析難以揭示影響因素作用的空間分布異質(zhì)的特點(diǎn)，缺少對(duì)典型降雨事件帶來(lái)的水質(zhì)影響特征的精細(xì)分析。地理統(tǒng)計(jì)學(xué)中的地理加權(quán)回歸（GWR，Geographic weighted regression）模型可探究不同自變量對(duì)因變量的影響強(qiáng)度和作用方向，并反應(yīng)影響因素的空間異質(zhì)特點(diǎn)[7]，在不同行業(yè)中有著廣泛應(yīng)用，如生態(tài)指數(shù)、空氣質(zhì)量與碳排放指數(shù)等環(huán)境空間異質(zhì)性分析。趙晶婭等[8]運(yùn)用美國(guó)洛杉磯市的人口、路網(wǎng)、交通等數(shù)據(jù)，利用地理加權(quán)回歸模型分析了多個(gè)影響因素對(duì)二氧化氮排放量的影響；陳曉輝[9]等利用地理加權(quán)回歸，分析福州市不同建成區(qū)域環(huán)境對(duì)生態(tài)指數(shù)的影響，定量分析了生態(tài)指數(shù)的空間異質(zhì)性；Wang[10]等使用GWR 模型，分析路網(wǎng)交通特征對(duì)城市空氣質(zhì)量的空間異質(zhì)性影響；袁長(zhǎng)偉[11]等基于我國(guó)各省運(yùn)輸業(yè)規(guī)模、技術(shù)以及結(jié)構(gòu)三大影響因子，構(gòu)建各省的碳排放量的地理加權(quán)回歸模型，進(jìn)而實(shí)現(xiàn)各地區(qū)碳減排壓力指數(shù)的劃分。

綜合背景及現(xiàn)狀，本文以武漢市東湖為例，采用地理加權(quán)回歸方法，對(duì)降雨事件前后水體懸浮物濃度變化與不同影響因子之間的影響特征進(jìn)行定量分析。武漢市東湖作為典型的大型城中湖，與城市和人群的地理交互關(guān)系密切，通過(guò)在東湖這一研究區(qū)域開(kāi)展降雨事件分析研究，有助于揭示不同影響因素對(duì)水質(zhì)影響的空間異質(zhì)性，為區(qū)域水環(huán)境治理提供依據(jù)。

1 研究區(qū)域與數(shù)據(jù)

東湖位于湖北省武漢市城區(qū)東部，具體可細(xì)分為水果湖、郭鄭湖、湯菱湖、廟湖、團(tuán)湖、后湖、喻家湖等子湖[12]，全湖面積在正常高水位19.78 米時(shí)為31.75 平方公里，是我國(guó)最大的城中湖泊（圖1），是武漢市的重要水源地，和知名生態(tài)旅游風(fēng)景區(qū)。西部水域與武昌區(qū)連接緊密，東、南、北三面為洪山區(qū)與青山區(qū)所圍繞。

1.1 遙感影像與實(shí)測(cè)水質(zhì)數(shù)據(jù)

本研究選擇衛(wèi)星多光譜遙感影像Sentinel-2 MSI 數(shù)據(jù)作為遙感水質(zhì)反演影像數(shù)據(jù)源，分別選取降雨前后2021 年1 月19 日、2021 年1 月29 日兩個(gè)時(shí)間的影像數(shù)據(jù)進(jìn)行了懸浮物濃度反演，東湖懸浮物反演模型使用盧東爍[13]構(gòu)建的東湖水質(zhì)參數(shù)波段組合反演模型。

1.2 影響因子數(shù)據(jù)

從研究區(qū)域的社會(huì)、經(jīng)濟(jì)、基礎(chǔ)設(shè)施等實(shí)際狀況出發(fā)，結(jié)合已有研究中對(duì)雨水徑流及其影響因素的探究，以及降雨通過(guò)雨水沖刷進(jìn)而影響湖泊水質(zhì)的過(guò)程分析，構(gòu)建了城市湖泊水質(zhì)初選影響因子集合如表1 所示，其數(shù)據(jù)來(lái)源包括網(wǎng)絡(luò)公開(kāi)數(shù)據(jù)及解譯數(shù)據(jù)[14]，研究中還將采用因子篩選方法進(jìn)一步挑選主要影響因子。

表1 城市湖泊水質(zhì)初選影響因子

2 研究方法

本研究中首先統(tǒng)計(jì)了東湖懸浮物濃度的空間自相關(guān)性，之后通過(guò)相關(guān)性分析和共線性檢驗(yàn)，并結(jié)合事件特點(diǎn)，進(jìn)行了影響因子的篩選，最后對(duì)懸浮物濃度與2.2 中不同因子進(jìn)行了地理加權(quán)回歸建模，基于建模結(jié)果分析不同因素對(duì)降雨前后懸浮物濃度變化的影響程度和空間特征，具體方法分別如下。

2.1 空間自相關(guān)分析

空間自相關(guān)性可以反映空間變量在各個(gè)區(qū)域上的相互依賴程度[15]。本文采用在地理統(tǒng)計(jì)研究中經(jīng)典的莫蘭指數(shù)探究水質(zhì)的全局空間自相關(guān)性，作為空間異質(zhì)性判定依據(jù)[16]，其計(jì)算公式為[17]：

式中：I 為自相關(guān)度量指標(biāo)；N 為研究區(qū)域內(nèi)的懸浮物濃度數(shù)據(jù)點(diǎn)位數(shù)量；zˉ為懸浮物濃度的平均值；zi、zj分別為懸浮物濃度數(shù)據(jù)的空間位置，wij為zi、zj的距離。

2.2 地理加權(quán)回歸模型

地理加權(quán)回歸（Geographically Weighted Regression）模型，將數(shù)據(jù)的空間位置嵌入到回歸參數(shù)中，建模得到不同地理位置各自變量與因變量間不同的空間關(guān)系[18]，模型如公式（2）所示：

在本文研究中，式中（ui，vi）為第i 個(gè)懸浮物濃度反演數(shù)據(jù)的坐標(biāo)（經(jīng)緯度投影坐標(biāo)，單位為m），β0（ui，vi）為第i個(gè)懸浮物濃度反演數(shù)據(jù)的常數(shù)項(xiàng)，βk（ui，vi）為第i 個(gè)懸浮物濃度反演數(shù)據(jù)內(nèi)第k 個(gè)自變量的估計(jì)系數(shù)，p 表示顯著影響水質(zhì)分布的自變量個(gè)數(shù)，εi為第i 個(gè)懸浮物濃度反演數(shù)據(jù)的隨機(jī)誤差。

2.3 相關(guān)性分析及共線性檢驗(yàn)

地理加權(quán)回歸模型求解時(shí)，如果因子之間存在局部多重共線性，地理加權(quán)回歸模型的結(jié)果將變得不穩(wěn)定，因此需對(duì)備選影響因子間的相關(guān)性分析與多重共線性檢驗(yàn)，以保證解釋變量的獨(dú)立性，提高模型準(zhǔn)確度。本文采用容忍度（Tolerance）和方差膨脹系數(shù)（variance inflation factor，VIF）指標(biāo)來(lái)共同度量影響因子間的多重共線性[19]，方差膨脹系數(shù)為容忍度的倒數(shù)，其計(jì)算公式為：

式中Ri表示該影響因子與其他影響因子之間的相關(guān)系數(shù)，i 為不同因子的編號(hào)。

3 結(jié)果和討論

3.1 典型降雨事件前后路域水質(zhì)路域水質(zhì)空間自相關(guān)分析

圖2 分別展示了降雨前后懸浮物濃度變化量空間自相關(guān)分析的局部莫蘭指數(shù)散點(diǎn)圖以及LISA 聚類地圖。降雨前后懸浮物濃度差異的莫蘭指數(shù)為0.482，表明其在空間上呈顯著的正空間自相關(guān)（P<0.01），即其在空間上趨于集聚。從LISA 聚類圖可發(fā)現(xiàn)其空間聚類特征明顯，東湖北部和東部呈現(xiàn)出高-高聚類，西部則呈現(xiàn)出低-低聚類，具有顯著的空間異質(zhì)性。

圖2 降雨前后懸浮物變化空間自相關(guān)分析圖

3.2 地理加權(quán)回歸模型因子相關(guān)性分析及共線性檢驗(yàn)

表1 所列初選影響因子之間的相關(guān)系數(shù)如表2 所示，從表2 可以發(fā)現(xiàn)，在選定的八個(gè)因子中，X4（道路面積）、X5（道路寬度）、X6（道路長(zhǎng)度）以及X8（車輛速度）具有顯著相關(guān)性，即這幾種初選影響因子間存在緊密聯(lián)系，特別是道路寬度與道路速度，兩者的相關(guān)系數(shù)為0.992，說(shuō)明車道寬度與車速和道路通行能力之間具有高相關(guān)性。

表2 初選影響因子相關(guān)系數(shù)

表3 為初選影響因子共線性診斷，八個(gè)初選因子中，建筑面積（BR）和道路面積（RR）初選因子的VIF 小于10，分別為8.946 與5.005，說(shuō)明這兩種因子與其他初選因子不存在多重共線性。森林面積（FR）、綠地面積（GR）、道路長(zhǎng)度（RL）和交叉口數(shù)（RC）這四種初選因子的VIF 大于10 小于100，這表明這四種因子與其他初選因子存在較強(qiáng)的多重共線性，需結(jié)合相關(guān)性分析結(jié)果判斷其是否保留。

表3 初選影響因子共線性診斷

綜合以上初選影響因子相關(guān)系數(shù)及多重共線性檢驗(yàn)結(jié)果，本研就選擇保留對(duì)于VIF 較大但相關(guān)性小于0.7 的影響因子，最終確定選取了建筑面積（BR）、森林面積（FR）、綠地面積（GR）、道路面積（RR）以及交叉口數(shù)（RC）這五個(gè)因子作為地理加權(quán)回歸的影響因子。

3.3 地理加權(quán)回歸模型結(jié)果分析

地理加權(quán)回歸R2以及影響因子系數(shù)分布圖如圖3 所示，總模型擬合精度整體上呈現(xiàn)由外向內(nèi)擴(kuò)散性遞減的圈層關(guān)系，即湖岸擬合效果好，湖心擬合效果差。主要原因是湖岸水質(zhì)受影響因子影響程度高，湖心水體由于水的流動(dòng)性，其水質(zhì)分布較為均勻，導(dǎo)致地理加權(quán)回歸模型無(wú)顯著關(guān)系。

圖3 地理加權(quán)回歸R2 分布圖

地理加權(quán)回歸影響因子系數(shù)分布圖如圖4 所示，總體來(lái)看，各影響因子的回歸系數(shù)在地理空間上存在明顯的非平穩(wěn)性，各影響因子對(duì)于水質(zhì)的影響能力和影響效應(yīng)有所不同，進(jìn)一步證明了影響因素的空間異質(zhì)性。

圖4 地理加權(quán)回歸自變量系數(shù)分布圖

①道路面積、交叉口數(shù)與懸浮物濃度在水體邊緣區(qū)域內(nèi)總體呈正相關(guān)，即道路面積越大、道路交叉口數(shù)量越多，該區(qū)域的懸浮物濃度越高。主要原因是道路面積越大，車流量越大，道路殘留污染物越多。

②建筑物面積與懸浮物濃度在水體邊緣區(qū)域內(nèi)總體呈正相關(guān)，即建筑物占地面積越大，該區(qū)域的懸浮物濃度越高。主要原因是建筑物高密度地區(qū)勢(shì)必會(huì)產(chǎn)生較多的屋面雨水徑流排放，進(jìn)而導(dǎo)致水體固態(tài)懸浮物污染加劇。

③森林、綠地面積與懸浮物濃度在水體邊緣區(qū)域內(nèi)均呈負(fù)相關(guān)。主要原因是森林與綠地對(duì)雨水具有截留、吸收以及凈化作用，減少了雨水徑流排放對(duì)水體帶來(lái)的污染。呈正相關(guān)部分區(qū)域周圍為建筑物、道路密集區(qū)，路面硬化面積大，樹(shù)木、綠地對(duì)雨水徑流的減少作用不明顯。

4 結(jié)論

本文以降雨事件前后，城市湖泊懸浮物變化為研究對(duì)象，采用地理加權(quán)回歸方法對(duì)其影響因子進(jìn)行分析。結(jié)果表明：

①考慮空間異質(zhì)性的GWR 模型對(duì)降雨前后懸浮物濃度變化與影響因子間具有較好的擬合度。

②各影響因子系數(shù)在地理空間上存在明顯的非平穩(wěn)性，進(jìn)一步證明了影響因素的空間異質(zhì)性。

③建筑物面積與懸浮物濃度在水體邊緣區(qū)域內(nèi)總體呈正相關(guān)，即建筑物占地面積越大，該區(qū)域的懸浮物濃度越高；森林、綠地面積與懸浮物濃度在水體邊緣區(qū)域內(nèi)均呈負(fù)相關(guān)；道路交叉口數(shù)、道路面積與懸浮物濃度在水體邊緣區(qū)域內(nèi)呈正相關(guān)。

在本研究的實(shí)驗(yàn)過(guò)程中仍存在一些不足之處，對(duì)影響因子的選取可加以擴(kuò)展，如從自然、社會(huì)、建成環(huán)境、交通系統(tǒng)等方面出發(fā)，構(gòu)建城市湖泊水質(zhì)影響指標(biāo)體系。盡可能全面衡量其影響因素，使得所建模型精度更高，擬合效果更理想，研究更加科學(xué)、準(zhǔn)確。