張 洪，黎海林，陳 震

（云南財(cái)經(jīng)大學(xué) 城市管理與資源環(huán)境學(xué)院，昆明650221）

在流域水環(huán)境變化原因的研究中，土地利用類型變化是一個(gè)非常重要的影響因素，如果能夠明確土地利用類型和水質(zhì)之間的關(guān)系，就可以通過(guò)適當(dāng)?shù)耐恋乩靡?guī)劃來(lái)確定合理的污染預(yù)防措施，從而解決水環(huán)境污染問(wèn)題，保護(hù)水質(zhì)。國(guó)外在這方面研究較多，Basta和Bower總結(jié)了14個(gè)土地利用變化對(duì)水量、水質(zhì)影響模型。比較詳細(xì)地描述了這些模型使用的土地類型、時(shí)空屬性、水文特征、水類型及測(cè)定的污染物。He將GIS和AGNPS模型結(jié)合，研究密歇根州一個(gè)農(nóng)業(yè)流域土壤侵蝕和N、P流失情況，并在模擬預(yù)測(cè)的結(jié)果基礎(chǔ)上提出了最佳的土地利用方法[1－3]。目前國(guó)內(nèi)也有這方面的研究，于興修以太湖上游浙江西苕溪流域?yàn)槔芯苛送恋乩茫脖蛔兓皩?duì)水質(zhì)的影響[4]，李文杰研究了梁子湖流域土地利用變化對(duì)流域水環(huán)境的影響[5]，張婷對(duì)白洋淀水質(zhì)時(shí)空變化及影響因子進(jìn)行了評(píng)價(jià)與分析[6]。但是云貴高原的滇池流域土地利用結(jié)構(gòu)及動(dòng)態(tài)變化對(duì)水質(zhì)的影響研究沒(méi)有，因此本文是LUCC及計(jì)量地理學(xué)理論在滇池流域的一次實(shí)證應(yīng)用，致力于研究滇池流域1974—2008年土地利用動(dòng)態(tài)變化及對(duì)其水環(huán)境的影響，旨在真實(shí)反映該區(qū)域的土地利用情況，同時(shí)得出這種變化對(duì)水環(huán)境的影響，為土地利用規(guī)劃、水污染治理及其它領(lǐng)域研究提供參考借鑒作用。

1 研究區(qū)域及研究方法

1.1 研究區(qū)域

滇池流域位于云貴高原中部，地處長(zhǎng)江、珠江和紅河三大水系分水嶺地帶，范圍包括昆明市盤(pán)龍、五華、西山、官渡四區(qū)及呈貢、晉寧、嵩明等縣市區(qū)部分地區(qū)。地理坐標(biāo)為東經(jīng)102°29′—103°01′，北緯24°29′—25°28′，流域面積為2 840.7km2。整個(gè)流域?yàn)槟媳遍L(zhǎng)、東西窄的湖盆地，地形可分為山地丘陵、淤積平原和滇池水域三個(gè)層次。滇池流域?qū)俦眮啛釒駶?rùn)季風(fēng)氣候，年均氣溫14.9℃，年降雨量931.8mm，全年無(wú)霜期227d，具有低緯高原的季風(fēng)氣候特征，冬無(wú)嚴(yán)寒、夏無(wú)酷暑，四季如春，干濕分明，垂直差異大的氣候特點(diǎn)。

1.2 研究方法及數(shù)據(jù)來(lái)源

1.2.1 研究方法 結(jié)合研究區(qū)域的土地利用方式特點(diǎn)，將土地利用按一級(jí)分類分為5大類，即：耕地、林地、建設(shè)用地、水域及其它地類。利用Edras遙感影像處理軟件和ArcGIS地理信息系統(tǒng)軟件對(duì)1974年、1988 年、1998 年、2008 年 4 個(gè) 時(shí) 期 Landsat－5（TM）遙感衛(wèi)片進(jìn)行校驗(yàn)配準(zhǔn)及增強(qiáng)處理，通過(guò)野外GPS實(shí)地考察，獲取遙感影像的解譯標(biāo)志，并對(duì)每個(gè)時(shí)期的影像進(jìn)行人機(jī)交互解譯，得出四個(gè)時(shí)期的土地利用現(xiàn)狀。利用ArcGIS地理信息系統(tǒng)軟件中的ArcMAP模塊進(jìn)行圖層疊加分析，根據(jù)馬爾柯夫模型建立土地利用空間動(dòng)態(tài)變化轉(zhuǎn)移矩陣。其中流域土地利用與水環(huán)境變化模擬采取單位根檢驗(yàn)與協(xié)整檢驗(yàn)來(lái)分析，并在此之前的基礎(chǔ)上，利用馬爾柯夫預(yù)測(cè)模型，通過(guò)Matlab軟件實(shí)現(xiàn)對(duì)后期的土地利用結(jié)果及對(duì)水環(huán)境的影響進(jìn)行預(yù)測(cè)。

1.2.2 數(shù)據(jù)來(lái)源 研究區(qū)數(shù)據(jù)來(lái)源于1974年、1988年、1998年及2008年Landsat－5（TM）遙感衛(wèi)片，采取TM＿Band＿4（R），TM＿Band＿3（G），TM＿Band＿2（B）波段合成影像，通過(guò)幾何校正及Edras遙感影像處理軟件和ArcGIS地理信息系統(tǒng)軟件人機(jī)交互解譯得來(lái)。

1.2.3 馬爾柯夫模型 馬爾柯夫方法是利用某一變量的現(xiàn)狀和動(dòng)向去預(yù)測(cè)未來(lái)的狀態(tài)及動(dòng)向的一種分析手段，數(shù)學(xué)表達(dá)為：

式中：P（n）——將來(lái)時(shí)刻n土地利用類型的狀態(tài)概率向量；P（0）——研究初始時(shí)刻土地利用類型的狀態(tài)概率向量，簡(jiǎn)稱初始狀態(tài)概率向量；Pij——土地類型i轉(zhuǎn)變?yōu)橥恋仡愋蚸的轉(zhuǎn)移概率向量[7]。

1.2.4 ADF單位根檢驗(yàn) ADF單位根檢驗(yàn)是檢驗(yàn)數(shù)據(jù)序列的平穩(wěn)性，如果存在同階平穩(wěn)的話，就可以對(duì)它們進(jìn)行協(xié)整檢驗(yàn)。由以下回歸方程構(gòu)成：

式中：β1——常數(shù)項(xiàng)；β2——t的系數(shù)；εt——純粹白噪音誤差項(xiàng)，而 ΔYt－1＝（Yt－1－Yt－2），ΔYt－2＝（Yt－2－Yt－3）等。假設(shè)δ＝0，Yt－1系數(shù)的估計(jì)t值服從τ統(tǒng)計(jì)量，即存在一個(gè)單位根，即時(shí)間序列是非平穩(wěn)的。假設(shè)不能接受δ＝0假設(shè)，可以使用通常的t檢驗(yàn)，即不存在單位根，該時(shí)間序列是平穩(wěn)的[8]。

1.2.5 協(xié)整檢驗(yàn) 采用Johansen檢驗(yàn)法進(jìn)行協(xié)整檢驗(yàn)，其基本思路是在多變量向量自回歸（VAR）系統(tǒng)回歸構(gòu)造兩個(gè)殘差的積矩陣，計(jì)算矩陣的有序本征值，根據(jù)本征值得出一系列的統(tǒng)計(jì)量判斷協(xié)整關(guān)系是否存在以及協(xié)整關(guān)系的個(gè)數(shù)[9]。

2 土地利用變化分析

2.1 土地利用現(xiàn)狀

根據(jù)前面所得數(shù)據(jù)，利用ArcGIS地理信息系統(tǒng)軟件中的面積統(tǒng)計(jì)功能，得到4個(gè)研究時(shí)期的5大地類的面積及其所占總面積的百分比（表1）。

從表1可以看出，4個(gè)研究時(shí)期的各類土地利用中林地所占比率最大，最大是1974年達(dá)到44.8%，這與流域的地質(zhì)地貌特點(diǎn)有很大關(guān)系。耕地面積和建設(shè)用地面積變化則較大，耕地由1974年的17.3%下降到2008年的9.3%，建設(shè)用地則從1974年的3.1%上升到2008年的14.0%。

表1 1974－2008年土地利用現(xiàn)狀結(jié)構(gòu)

2.2 土地利用類型變化

土地利用變化分析主要是對(duì)4個(gè)時(shí)期土地利用現(xiàn)狀空間信息在ArcGIS地理信息系統(tǒng)軟件下進(jìn)行空間疊加，利用馬爾柯夫轉(zhuǎn)移矩陣，通過(guò)統(tǒng)計(jì)分析、計(jì)算整理后得到各個(gè)時(shí)期的土地利用現(xiàn)狀圖及土地利用轉(zhuǎn)移矩陣（附圖9、表2）。

所得轉(zhuǎn)移矩陣表中，數(shù)據(jù)所在列代表1974年土地利用類型，行代表2008年土地利用類型。數(shù)字代表的是1974的土地利用類型轉(zhuǎn)變?yōu)?008年各類利用類型的面積。根據(jù)轉(zhuǎn)移矩陣不難看出，在1974—2008年期間，耕地轉(zhuǎn)化為建設(shè)用地的面積為212.86 km2，其中1988—2008年中，耕地轉(zhuǎn)化為建設(shè)用地的面積為為128.98km2；其它地類轉(zhuǎn)化為建設(shè)用地的面積為102.86km2，水域及林地呈遞減趨勢(shì)，但是減少的面積很小，整體變化不大。

3 基于馬爾柯夫模型的土地利用預(yù)測(cè)

根據(jù)馬爾柯夫預(yù)測(cè)模型，應(yīng)用Matlab軟件對(duì)各個(gè)時(shí)期的轉(zhuǎn)移概率進(jìn)行計(jì)算，并將預(yù)測(cè)相應(yīng)結(jié)果與實(shí)際值比較，結(jié)合實(shí)際情況，采用1988—2008這個(gè)時(shí)間段的里1988年作為初始狀態(tài)年，以這個(gè)時(shí)間段的轉(zhuǎn)移概率作為轉(zhuǎn)移概率向量，預(yù)測(cè)得到1998年、2008年的土地利用結(jié)構(gòu)與1998年、2008年實(shí)際的土地利用結(jié)構(gòu)進(jìn)行對(duì)比，模擬值和實(shí)際值十分接近，兩個(gè)時(shí)期的相對(duì)誤差平均值在0.05%內(nèi)。由此可見(jiàn)利用馬爾柯夫過(guò)程預(yù)測(cè)滇池流域土地利用結(jié)構(gòu)是可行的。得到的預(yù)測(cè)結(jié)果見(jiàn)表3。

表3 馬爾柯夫鏈預(yù)測(cè)的各土地利用類型面積 km2

馬爾柯夫模型預(yù)測(cè)結(jié)果表明，2015—2030年研究區(qū)內(nèi)建設(shè)用地面積進(jìn)一步增加，林地、水域和其他地類進(jìn)一步減少，但減少的幅度在逐漸變小，耕地面積也繼續(xù)減少。

4 流域土地利用演變與水環(huán)境變化的關(guān)系

土地利用是人類對(duì)自然環(huán)境最直接的一種作用方式，如果利用方式不當(dāng)將會(huì)導(dǎo)致大量的硝酸鹽、磷酸鹽等物質(zhì)未經(jīng)處理排入流域水體，進(jìn)一步造成對(duì)流域水環(huán)境的污染。在此以滇池流域土地利用變化的研究為基礎(chǔ)，結(jié)合查找相關(guān)年份的土地利用數(shù)據(jù)，利用1994—2008年以來(lái)可以獲得的水質(zhì)資料，分析滇池流域土地利用變化對(duì)水環(huán)境所帶來(lái)的影響。

4.1 流域土地利用與水環(huán)境變化模擬

水環(huán)境主要由地表水和地下水兩部分組成，本文只對(duì)滇池的地表水環(huán)境進(jìn)行相應(yīng)分析。滇池水環(huán)境的評(píng)價(jià)標(biāo)準(zhǔn)是參考中華人民共和國(guó)環(huán)境保護(hù)總局頒布的《中華人民共和國(guó)地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB3838－2002）。水質(zhì)監(jiān)測(cè)的常規(guī)項(xiàng)目包括：水溫、電導(dǎo)率、pH、溶解氧（DO）、總磷（TP）、總氮（TN）、氨氮（NH3－N）、高錳酸鹽指數(shù)、生化需氧量BOD，揮發(fā)酚、石油類、重金屬和細(xì)菌等指標(biāo)[6]。結(jié)合本次研究實(shí)際情況，對(duì)從監(jiān)測(cè)部門(mén)所獲取的1994—2008年的總氮（TN）、總磷（TP）、生化需氧量（BOD）及化學(xué)需氧量（COD）資料采用時(shí)間序列方法進(jìn)行分析，研究滇池流域土地利用變化對(duì)滇池水環(huán)境產(chǎn)生的影響。

4.1.1 單位根檢驗(yàn) 采用監(jiān)管部門(mén)調(diào)研得到的土地?cái)?shù)據(jù)，將滇池流域的土地類型分為城鎮(zhèn)村工礦建設(shè)用地（CH）、耕地（GE）、有林地（LIN）、牧草地（MU）、濕地（SHI）、園地（YUAN）6大類，分別采用各類土地類型占流域土地的比重，對(duì)1994—2008年滇池流域的總氮（TN）、總磷（TP）、生化需氧量（BOD）、化學(xué)需氧量（COD）以及各土地類型的時(shí)間序列采用ADF法進(jìn)行單位根檢驗(yàn)，計(jì)算結(jié)果如表4所示。

由表4可以看出總氮量與園地序列為0階單整，其它序列均為2階單整。

表4 各變量ADF檢驗(yàn)

4.1.2 協(xié)整檢驗(yàn) 通過(guò)Johansen檢驗(yàn)計(jì)算得到總磷（TP）、生化需氧量（BOD）、化學(xué)需氧量（COD）與城鎮(zhèn)村工礦建設(shè)用地（CH）、耕地（GE）、有林地（LIN）、牧草地（MU）、濕地（SHI）、園地（YUAN）之間存在協(xié)整關(guān)系。其關(guān)系為：

其中（1）（2）（3）式下括號(hào)內(nèi)數(shù)為t值，R2表示調(diào)整后的擬合優(yōu)度，DW為durbin watson檢驗(yàn)值?？梢缘玫皆?%顯著性水平下TP和BOD與土地利用類型之間存在相關(guān)關(guān)系，城鎮(zhèn)建設(shè)用地及獨(dú)立工礦用地對(duì)總氮量的彈性為15.31，即城鎮(zhèn)建設(shè)用地每增加666.7 m2，會(huì)對(duì)總氮量產(chǎn)生15.31mg／L的正向影響，同理，耕地對(duì)總氮量彈性為10.47，牧草地對(duì)總氮量彈性為－156.74。城鎮(zhèn)建設(shè)用地及獨(dú)立工礦用地對(duì)BOD的影響為543.79，耕地對(duì)BOD的影響為－272.49，有林地對(duì)BOD的影響為160.88。在20%顯著性水平下COD與土地利用類型之間存在相關(guān)關(guān)系，城鎮(zhèn)建設(shè)用地及獨(dú)立工礦用地、耕地、濕地對(duì)COD的彈性分別為173.88，－100.75和－4920.59。可以認(rèn)為，城鎮(zhèn)村工礦建設(shè)用地和耕地顯著提高了流域總磷量的排放，然而牧草地顯著減少了流域總磷量的影響；城鎮(zhèn)村工礦建設(shè)用地及有林地顯著提高了流域BOD的污染，而耕地有效地減少了流域BOD的污染；城鎮(zhèn)村工礦建設(shè)用地顯著提高了流域COD的排放，而耕地和濕地有效減少了COD的排放。

4.2 流域土地利用與水環(huán)境變化預(yù)測(cè)

采用馬爾柯夫鏈分別預(yù)測(cè)2015年、2020年、2025年、2030年4期土地利用狀況后，采用模型（1），（2），（3）分別對(duì)滇池流域TP，BOD，COD進(jìn)行預(yù)測(cè)，結(jié)果見(jiàn)表5。

表5 4個(gè)時(shí)期水質(zhì)污染狀況預(yù)測(cè) mg／L

由以上預(yù)測(cè)可以看出，隨著年份的增長(zhǎng)，水質(zhì)污染狀況指標(biāo)呈現(xiàn)持續(xù)惡化的趨勢(shì)，主要原因是由于隨著年份的增長(zhǎng)，流域城鎮(zhèn)村工礦建設(shè)用地所占比例不斷擴(kuò)大，而流域內(nèi)耕地、有林地、牧草地、濕地等地類不斷減少，導(dǎo)致了環(huán)境狀況的持續(xù)惡化。

5 結(jié)論

通過(guò)對(duì)1974—2008年4個(gè)時(shí)期的分析，揭示了滇池流域近30a土地利用／覆被變化情況以及這種變化對(duì)滇池水環(huán)境的影響。結(jié)果表明，滇池流域土地利用類型變化顯著，土地利用對(duì)滇池水質(zhì)有一定的影響。土地利用相互轉(zhuǎn)化過(guò)程中，耕地、建設(shè)用地和其它地類的變化幅度較大，水域和林地變化不明顯。其中，耕地轉(zhuǎn)化為建設(shè)用地的面積為212.86km2，其它地類轉(zhuǎn)化為建設(shè)用地的面積為102.86km2。各時(shí)間段中，1998—2008年研究期間的土地利用變化最大，同樣表現(xiàn)為耕地的大量減少和建設(shè)用地的急劇增加。

馬爾柯夫模型預(yù)測(cè)結(jié)果表明，2015—2030年研究區(qū)內(nèi)建設(shè)用地面積進(jìn)一步增加，林地、水域和其它地類進(jìn)一步減少，但減少的幅度在逐漸變小，耕地面積也繼續(xù)減少。

通過(guò)單位根和協(xié)整檢驗(yàn)可以看出土地利用類型的變化對(duì)滇池水環(huán)境變化具有明顯的影響作用，在5%顯著性水平下TP和BOD與土地利用類型之間存在相關(guān)關(guān)系，建設(shè)用地對(duì)總氮量的彈性為15.31，即城鎮(zhèn)建設(shè)用地每增加0.067hm2，會(huì)對(duì)總氮量產(chǎn)生15.31mg／L的正向影響，耕地對(duì)總氮量彈性為10.47，牧草地對(duì)總氮量彈性為－156.74。建設(shè)用地對(duì)BOD的影響為543.79，耕地對(duì)BOD影響為－272.49，有林地對(duì)BOD影響為160.88。在20%顯著性水平下COD與土地利用類型之間存在的相關(guān)關(guān)系，建設(shè)用地、耕地、濕地對(duì)COD的彈性分別為173.88，－100.75，－4 920.59?？梢哉J(rèn)為，城鎮(zhèn)村工礦建設(shè)用地和耕地顯著提高了流域總磷量的排放，然而牧草地顯著減少了流域總磷量的影響；城鎮(zhèn)村工礦建設(shè)用地及有林地顯著提高了流域BOD的污染，而耕地有效地減少了流域BOD的污染；城鎮(zhèn)村工礦建設(shè)用地顯著提高了流域COD的排放，而耕地和濕地有效減少了COD的排放。

[1]張殿發(fā)，王世杰.土地利用／土地覆被變化對(duì)長(zhǎng)江流域水環(huán)境的影響研究[J].地域研究與開(kāi)發(fā)，2003，2（1）：69－72.

[2]郭旭東，陳利頂.土地利用／土地覆被變化對(duì)區(qū)域生態(tài)環(huán)境的影響[J].環(huán)境科學(xué)進(jìn)展，1999，12（6）：66－75.

[3]趙米金，徐濤.土地利用／土地覆被變化環(huán)境效應(yīng)研究[J].水土保持研究，2005，12（1）：43－46.

[4]于興修，楊桂山.典型流域土地利用／覆被變化及對(duì)水質(zhì)的影響[J].長(zhǎng)江流域資源與環(huán)境，2003，12（3）：211－217.

[5]張婷，劉靜玲，王雪梅.白洋淀水質(zhì)時(shí)空變化及影響因子評(píng)價(jià)與分析[J].環(huán)境科學(xué)學(xué)報(bào)，2010，30（2）：261－267.

[6]李文杰.梁子湖流域土地利用變化對(duì)流域水環(huán)境的影響[D].武漢：華中師范大學(xué)，2009.

[7]劉瓊，歐名豪，彭曉英.基于馬爾柯夫過(guò)程的區(qū)域土地利用結(jié)構(gòu)預(yù)測(cè)研究[J].南京農(nóng)業(yè)大學(xué)學(xué)報(bào)，2005，28（3）：107－112.

[8]李子奈，葉阿忠.高等計(jì)量經(jīng)濟(jì)學(xué)[M].北京：高等教育出版社，2003.

[9]達(dá)莫達(dá)爾·N·古扎拉蒂.計(jì)量經(jīng)濟(jì)學(xué)基礎(chǔ)[M].4版.北京：中國(guó)人民大學(xué)出版社，2005.