何云玲，熊巧利，羅賢 ，李同艷，余嵐

1. 云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院，云南 昆明 650091；2. 云南大學(xué)亞洲國(guó)際河流中心，云南 昆明 650091；3. 云南省國(guó)際河流與跨境生態(tài)安全重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650091

水質(zhì)是水體的物理（色度、懸浮物等）、化學(xué)（有機(jī)和無(wú)機(jī)物含量）和生物（細(xì)菌、微生物、浮游生物等）方面的特征及其組成狀況。生活和工農(nóng)業(yè)生產(chǎn)廢水中大量氮、磷等營(yíng)養(yǎng)物質(zhì)進(jìn)入淡水水體后導(dǎo)致藻類（主要由藍(lán)藻、綠藻組成）大量繁殖的“水華”現(xiàn)象，是湖泊富營(yíng)養(yǎng)化的典型特征。水華一直是評(píng)價(jià)湖泊水體水質(zhì)狀況的重要指標(biāo)之一。自上個(gè)世紀(jì) 90年代以來(lái)，中國(guó)太湖、巢湖、滇池等湖泊藍(lán)藻水華暴發(fā)相繼出現(xiàn)，對(duì)飲用水源安全造成較大影響（李旭文等，2011）。在采取一系列治理措施后，“三湖”水質(zhì)有所改善，富營(yíng)養(yǎng)化程度較歷史峰值均有一定降低，然而以藍(lán)藻為優(yōu)勢(shì)類群的水華現(xiàn)象卻依然時(shí)有發(fā)生（Zhang et al.，2015；Zhang et al.，2016）。如何監(jiān)測(cè)和控制藻類水華，是湖泊富營(yíng)養(yǎng)化研究中的重點(diǎn)和難點(diǎn)（楊桂山等，2010；蔣大林，2015）。

水質(zhì)遙感是指利用遙感技術(shù)對(duì)水質(zhì)指標(biāo)進(jìn)行檢測(cè)，并根據(jù)遙感波段的信息與水質(zhì)指標(biāo)的光譜特征，建立反演模型對(duì)水質(zhì)指標(biāo)的濃度進(jìn)行預(yù)測(cè)（冼翠玲等，2017）。相比于傳統(tǒng)的水質(zhì)監(jiān)測(cè)手段，遙感監(jiān)測(cè)技術(shù)具有掃描面積大、資料更新周期短、成本相對(duì)較低和便于進(jìn)行長(zhǎng)期實(shí)時(shí)監(jiān)測(cè)等優(yōu)點(diǎn)。自20世紀(jì) 70年代初期開始，國(guó)外就廣泛地開展水質(zhì)遙感反演研究（Song et al.，2012）。國(guó)內(nèi)的相關(guān)研究在最近幾年發(fā)展迅速，聞建光等（2006）基于Hyperion星載高光譜遙感影像，運(yùn)用波段比值法、一階微分處理技術(shù)構(gòu)建了太湖葉綠素a濃度的反演模型。Shi et al.（2015）發(fā)現(xiàn)MODIS Aqua影像與總懸浮物在645 nm波長(zhǎng)處具有極強(qiáng)相關(guān)性，據(jù)此發(fā)展了適應(yīng)性極強(qiáng)的總懸浮物經(jīng)驗(yàn)反演模型，并在太湖取得了良好效果。歸一化植被指數(shù)（Normalized Difference Vegetation Index，NDVI），是遙感影像中近紅外波段的反射值與紅光波段的反射值之差與兩者之和的比值，常被用于確定目標(biāo)區(qū)綠色植物覆蓋程度的指標(biāo)值（陳旭等，2015）；徐京萍等（2008）發(fā)現(xiàn)利用MODIS數(shù)據(jù)的NDVI能夠確定湖泊藍(lán)藻的時(shí)空分布；周立國(guó)等（2008）基于MODIS衛(wèi)星影像，采用比值模型方法確定藍(lán)藻分布的相對(duì)濃度信息；藍(lán)藻暴發(fā)時(shí)引發(fā)水體溫度、色度和透明度等一系列物理性質(zhì)發(fā)生變化，進(jìn)而導(dǎo)致水體反射波譜特性也發(fā)生變化，利用這一特性，就可以基于NDVI值對(duì)藍(lán)藻水華分布范圍和面積進(jìn)行識(shí)別（魯韋坤等，2009）。

滇池位于云南省會(huì)昆明市的南部，是該市重要的水源地，支持著昆明市工業(yè)、農(nóng)業(yè)、旅游業(yè)等發(fā)展，歷史上素有“高原明珠”之稱。但隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展，大量的工農(nóng)業(yè)廢水和生活污水排入滇池，加之入湖支流水系較多，出流水系較少，湖泊換水周期較長(zhǎng)，致使滇池水質(zhì)惡化，水體富營(yíng)養(yǎng)化問題突出（王杰等，2018），成為全國(guó)湖泊富營(yíng)養(yǎng)化治理的重點(diǎn)對(duì)象。已有學(xué)者對(duì)滇池水質(zhì)的時(shí)空分布特征做過研究，主要集中在滇池水質(zhì)傳統(tǒng)方法采樣監(jiān)測(cè)（梁中耀等，2014）、葉綠素 a濃度分析（韋力元等，2016）、NDVI值域的劃分（謝國(guó)清等，2010）、水質(zhì)等級(jí)預(yù)測(cè)模型（徐玉妃等，2018）、以及氣象條件與水質(zhì)之間的關(guān)系方面（張虎才等，2017）。本研究基于2000－2017年的NDVI不同值域?qū)Φ岢厮{(lán)藻水華信息進(jìn)行提取，并從逐年、逐月到逐旬詳細(xì)分析滇池水華的時(shí)空分布規(guī)律以及變化特征，以期為滇池的治理及水資源的開發(fā)利用提供數(shù)據(jù)支持，也為建立實(shí)時(shí)監(jiān)測(cè)和預(yù)警系統(tǒng)奠定基礎(chǔ)。

1 滇池水域概況

滇池，亦稱昆明湖，是云南省最大的淡水湖，位于昆明市主城區(qū)的西南方向，湖面海拔1886 m，平均水深5 m，最深8 m；湖水在西南?？谛钩?，稱螳螂川，為長(zhǎng)江上游干流金沙江支流普渡河上源，屬長(zhǎng)江流域金沙江水系，為地震斷層陷落型的湖泊。南北長(zhǎng)39 km，東西最寬為13 km，流域面積約為2920 km2，容水量為15.7億立方米。滇池流域地處低緯度、高海拔地區(qū)，屬于熱帶高原季風(fēng)氣候；年均溫 14.7 ℃，日溫差較大，降雨年內(nèi)分配特點(diǎn)明顯，分為干濕分明的雨季和旱季。流域內(nèi)水系眾多，有 30余條河流呈向心狀注入滇池，其中徑流面積大于100 km2的河流有盤龍江、寶象河、洛龍河、撈魚河、大河、柴河、東大河。滇池流域的地理位置及基本概況如圖1所示。

2 數(shù)據(jù)來(lái)源與處理

2.1 數(shù)據(jù)來(lái)源

圖1 滇池流域地理位置及基本概況圖Fig. 1 Geographical location and basic profile of Dianchi Basin

利用中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：//www.resdc.cn/）2015年中國(guó)土地利用遙感監(jiān)測(cè)數(shù)據(jù)提取滇池水體邊界，分辨率為30 m×30 m。DEM 數(shù)據(jù)來(lái)自于美國(guó)國(guó)防部國(guó)家測(cè)繪局公布的SRTM-DEM 數(shù)據(jù)（http：//www.cgiar-csi.org），分辨率為90 m×90 m。MODIS數(shù)據(jù)來(lái)自于美國(guó)國(guó)家航空航天局 NASA（National Aeronautics and Space Administration）（https：//ladsweb.modaps.eosdis.nasa.gov/)的Terra遙感衛(wèi)星2000年2月18日－2017年12月31日的MOD13Q1數(shù)據(jù)產(chǎn)品，選用16 d合成250 m×250 m分辨率的NDVI數(shù)據(jù)產(chǎn)品；該數(shù)據(jù)產(chǎn)品已經(jīng)完成影像數(shù)據(jù)的輻射定標(biāo)、大氣校正、幾何校正、氣溶膠、臭氧吸收訂正及去云等預(yù)處理。

2.2 數(shù)據(jù)處理

將MODIS影像數(shù)據(jù)導(dǎo)入MRT軟件中進(jìn)行影像拼接并進(jìn)行UTM投影轉(zhuǎn)換。借助ENVI 5.1、ArcGIS 10.2軟件進(jìn)行空間裁剪、掩膜、重分類得到 2000－2017年的時(shí)間序列遙感數(shù)據(jù)集。并采用最大值合成法（maximum value composite，MVC）建立NDVI時(shí)間序列柵格數(shù)據(jù)集，分別求出月、季、年 NDVI值，公式如下：式中，i為月序號(hào)，取值范圍1－12月；NDVIi為第 i月 NDVI最大值；NDVIk為月上下旬 NDVI值，k取值為1－24；j為季序號(hào)，取值范圍為1－4，NDVIj為第j季NDVI最大值；y為年序號(hào)，取值范圍為1－18，NDVIy為第y年NDVI最大值，yNDVIj為第y年第j季NDVI最大值。

采用最小二乘法計(jì)算數(shù)據(jù)集中所有像元的NDVI與時(shí)間的回歸斜率Slope，若Slope＞0，表示隨時(shí)間變化NDVI呈增加趨勢(shì)，并且數(shù)值越大增幅越大；若Slope＜0則代表NDVI呈下降趨勢(shì)。計(jì)算公式為：

式中，i為年序號(hào)；n為年跨度，即18年；NDVIi為第i年NDVI值。

采用謝國(guó)清等（2010）在滇池藍(lán)藻水華光譜特征研究中關(guān)于NDVI閾值劃分水華的等級(jí)方法，即NDVI值大于-0.1的區(qū)域，為藍(lán)藻水華區(qū)域，其中-0.1≤NDVI≤0.2為輕度水華，表征像元內(nèi)水華覆蓋度為0－30%；0.2＜NDVI≤0.4為中度水華，表征像元內(nèi)藍(lán)藻水華覆蓋度為 31%－80%；NDVI＞0.4為重度水華，表征像元內(nèi)藍(lán)藻水華覆蓋度為 81%－100%。應(yīng)用GetData Graph Digitizer 2.26軟件對(duì)以往實(shí)測(cè)研究中滇池富營(yíng)養(yǎng)化指數(shù)變化特征進(jìn)行數(shù)值提取，對(duì)比分析及驗(yàn)證基于遙感識(shí)別基礎(chǔ)上NDVI表征的滇池水華特征。

3 研究結(jié)果與分析

3.1 年際變化

從2000－2017年滇池NDVI的逐年變化趨勢(shì)圖中可以看出（圖2），通過實(shí)測(cè)研究得到的2008－2014年滇池富營(yíng)養(yǎng)化指數(shù)與同時(shí)段基于 NDVI表征的水華特征變化趨勢(shì)相比，出現(xiàn)較高值和較低值的波動(dòng)變化趨勢(shì)相似。其次，除了秋季有緩慢增加的趨勢(shì)外，其他季節(jié)和全年的NDVI整體上均呈減少趨勢(shì)，但均未通過P＜0.05的顯著性檢驗(yàn)。據(jù)統(tǒng)計(jì)，冬季的NDVI減少趨勢(shì)最明顯（-0.0028/a），夏季次之（-0.0022/a），春季變化趨勢(shì)不明顯。另外，夏季NDVI與年NDVI的變化趨勢(shì)較為一致，18年間NDVI經(jīng)歷幾次增加和減少的波動(dòng)?？傮w而言，從2000年開始滇池年NDVI呈減少趨勢(shì)，藍(lán)藻水華覆蓋度呈降低趨勢(shì)。

圖2 2000－2017年滇池NDVI值與富營(yíng)養(yǎng)化指數(shù)變化圖Fig. 2 Annual NDVI and eutrophication index in Dichi from 2000 to 2017

通過數(shù)據(jù)分析得到 2000－2017年滇池 NDVI的年平均值空間等級(jí)分布及變化斜率（圖 3），多年平均的輕度水華面積占滇池總面積的54.99%，中度和重度水華面積分別占總面積的 27.42%和17.59%，從南到北，滇池水華呈現(xiàn)從輕度到中度、重度的變化格局（圖3a）。滇池大部分區(qū)域NDVI變化斜率主要分布在 0－-0.1之間，NDVI整體上呈現(xiàn)不同程度的減小趨勢(shì)（圖3b），呈減小趨勢(shì)的區(qū)域占研究區(qū)總面積55.64%；但有區(qū)域差異，NDVI減小，表征藍(lán)藻水華覆蓋度減少的地區(qū)主要分布在滇池北岸的斷橋，滇池中部的觀音山西、白漁口、?？谖鞯鹊?；NDVI增大，表征藍(lán)藻水華覆蓋度增加的地區(qū)主要分布在滇池北岸的草海、羅家營(yíng)，滇池東南部的觀音山東、滇池南等地。

圖3 2000－2017年滇池年平均NDVI空間分布（a）及變化斜率（b）圖Fig. 3 Spatial distribution of average NDVI and variation slope in Dianchi from 2000 to 2017

3.2 四季變化

按照上述 NDVI閾值劃分藍(lán)藻水華等級(jí)的方法，統(tǒng)計(jì)各季節(jié)滇池水華的空間分布特征，仍然表現(xiàn)出“北重南輕”的分布態(tài)勢(shì)（圖 4）；春、夏、秋三季滇池以輕度水華為主，分別占滇池總面積的86.49%、72.11%和 83.23%；冬季則以無(wú)水華區(qū)域所占總面積比例最大（69.34%），冬季的重度水華區(qū)域僅占1.57%；而中度和重度水華則主要出現(xiàn)在夏季，所占比例最高，分別占滇池總面積的13.89%和14%。總體而言，滇池的重度水華區(qū)域主要分布在滇池的北部（斷橋、草海）和滇池的邊界；中度水華區(qū)域主要分布在滇池的中部羅家營(yíng)和觀音山西；輕度水華區(qū)域主要集中在滇池南部的觀音山中、海口西和滇池南。

圖4 不同季節(jié)滇池水華特征的空間變化Fig. 4 Spatial distribution of water bloom in four seasons in Dianchi

3.3 月變化

通過實(shí)測(cè)得到的滇池逐月富營(yíng)養(yǎng)化指數(shù)與基于NDVI表征的逐月水華特征變化趨勢(shì)相比較（圖5），可以看出二者出現(xiàn)較大值與較小值的波動(dòng)變化相吻合，均是在7月達(dá)到最大值，12月達(dá)到最小值；5－8月滇池富營(yíng)養(yǎng)化指數(shù)為高值區(qū)；同時(shí)，NDVI也表現(xiàn)為高值區(qū)，代表藍(lán)藻水華覆蓋度大，富營(yíng)養(yǎng)化程度嚴(yán)重。

圖5 滇池逐月NDVI值與富營(yíng)養(yǎng)化指數(shù)變化圖Fig. 5 Monthly NDVI value and eutrophication index in Dianchi

再?gòu)闹鹪碌牡岢豊DVI值空間變化特征可以看出（圖6），每年的1－3月除了靠近滇池北岸附近區(qū)域以外的大部分地區(qū)NDVI值很小，從4月開始，NDVI開始逐漸增大，在不同地區(qū)出現(xiàn)水華；7－8月份NDVI值較大，滇池水華較嚴(yán)重；9月份開始滇池NDVI值又逐漸減小，水華面積縮減。

統(tǒng)計(jì)不同程度水華出現(xiàn)的面積比例在各月的變化（圖7），滇池11－4月以無(wú)水華區(qū)域?yàn)橹鳎娣e比例占絕大優(yōu)勢(shì)（尤其11－3月占80%以上）；5－10月則以輕度水華為主（尤其5－9月所占面積比例為80%左右）（圖7a）。中度水華主要出現(xiàn)在5－10月之間，所占總面積比例均小于12%；重度水華則出現(xiàn)在 7－9月之間，所占總面積比例較?。ň?%以下），尤其重度水華出現(xiàn)的時(shí)間從5月開始到8月達(dá)到峰值（月平均值為7.61%），10月開始逐漸消失（圖7b）。

3.4 旬變化

從上述滇池水華特征的四季和月變化分析中可知，藍(lán)藻水華的暴發(fā)主要集中在夏季，進(jìn)一步對(duì)夏季各旬的NDVI值空間變化進(jìn)行研究（圖8）。藍(lán)藻水華的暴發(fā)時(shí)間7月上旬和8月下旬，其重度水華面積分別高達(dá)11.4%和14.0%。水華較嚴(yán)重的區(qū)域仍集中在滇池的草海、斷橋和滇池的邊界地區(qū)。

4 討論與結(jié)論

4.1 討論

滇池在維持流域生態(tài)平衡、滿足生產(chǎn)生活用水、減輕洪澇災(zāi)害和提供豐富水產(chǎn)品等方面發(fā)揮著不可替代的作用。20世紀(jì)60年代其水質(zhì)為Ⅱ類，水產(chǎn)資源豐富，70年代為Ⅲ類，80年代為Ⅴ類，90年代為劣Ⅴ類，“九五”開始，滇池被列為全國(guó)重點(diǎn)治理湖泊。伴隨遙感監(jiān)測(cè)技術(shù)的進(jìn)步和迅速發(fā)展，對(duì)滇池的相關(guān)研究也逐漸開展，謝國(guó)清等（2010）利用假彩色合成法以及 NDVI閾值進(jìn)行了滇池藍(lán)藻水華劃分，但未通過時(shí)間序列NDVI值對(duì)水華時(shí)空變化特征進(jìn)行研究；王杰等（2018）通過傳統(tǒng)采樣方法對(duì)滇池 2008－2014年富營(yíng)養(yǎng)化程度進(jìn)行分析，但是數(shù)據(jù)樣本較少；蔣大林（2015）通過浮游藻指數(shù)（FAI）對(duì)滇池水華時(shí)空分布特征及其驅(qū)動(dòng)因子進(jìn)行研究，主要集中在FAI閾值提取和水華面積的簡(jiǎn)單統(tǒng)計(jì)方面。本研究基于遙感數(shù)據(jù)，借鑒上述學(xué)者研究成果，對(duì)NDVI值域進(jìn)行了分類，對(duì)18年來(lái)滇池水華的時(shí)空分布特征進(jìn)行了更深入、詳細(xì)的分析研究，研究結(jié)果和侯秀麗等（2018）、余麗燕等（2016）、余佑金等（2017）的結(jié)論相一致，表明基于NDVI值域判斷滇池水華時(shí)空變化特征具有一定可行性。

滇池地處昆明城市下游，是昆明盆地最低凹地帶，接納生活污水、工業(yè)廢水和含有農(nóng)藥化肥的農(nóng)業(yè)污水，造成研究區(qū)藍(lán)藻水華暴發(fā)格局呈現(xiàn)“北重南輕”的分布態(tài)勢(shì)。此外，昆明地區(qū) 7－8月氣溫較高，風(fēng)速較小，氣候條件加速了藍(lán)藻水華的暴發(fā)，導(dǎo)致了滇池水華在月變化上存在顯著差異。2000年開始，昆明市開始大力推進(jìn)滇池污染治理的各項(xiàng)工作，通過鋪設(shè)截污管道、建污水處理廠、牛欄江-滇池補(bǔ)水工程通水，成效逐步顯現(xiàn)，水質(zhì)逐漸提升。這和本研究中2000－2017年滇池年平均NDVI總體呈現(xiàn)減少趨勢(shì)，藍(lán)藻水華覆蓋度有所降低的結(jié)論是一致的。滇池北部的草海、斷橋、暉彎和觀音山地區(qū)過去是滇池流域富營(yíng)養(yǎng)化程度較重、水質(zhì)較差的地區(qū)，但是經(jīng)過整治，這些地區(qū)的水質(zhì)有變好趨勢(shì)，反而是過去富營(yíng)養(yǎng)化程度較輕、水質(zhì)較好的滇池南部羅家營(yíng)南端部分地區(qū)水華有加重趨勢(shì)，這在未來(lái)滇池水質(zhì)綜合防治中值得關(guān)注。

圖6 滇池水華的月變化空間分布圖Fig.6 Spatial distribution of monthly water bloom in Dianchi

雖然基于MODIS遙感數(shù)據(jù)研究水華變化有其局限性，水體識(shí)別主要依賴于光譜信息導(dǎo)致NDVI無(wú)法區(qū)分藍(lán)藻區(qū)域和高渾濁水體，不能有效識(shí)別低密度藍(lán)藻水體，研究結(jié)果易擴(kuò)大或縮小藍(lán)藻的分布范圍（郭鵬等，2012）。然而，不可否認(rèn)的是，MODIS數(shù)據(jù)以其波段窄、信噪比高、波譜敏感性好以及可免費(fèi)獲取等優(yōu)點(diǎn)可為實(shí)時(shí)、準(zhǔn)確、經(jīng)濟(jì)地獲取藍(lán)藻水華信息提供便捷。在今后的研究中需探索更精確混合像元估算藍(lán)藻水華面積的方法；需充分利用不同遙感數(shù)據(jù)的優(yōu)點(diǎn)，提高數(shù)據(jù)時(shí)空分辨率。不同學(xué)者對(duì)NDVI值域劃分有不同的劃分依據(jù)，如謝國(guó)清等（2010）認(rèn)為NVDI≥-0.1為藍(lán)藻水華區(qū)域，孫佳麗等（2010）認(rèn)為NDVI＞0為藍(lán)藻水華區(qū)域。本研究采用的是文獻(xiàn)引用量達(dá)高頻次的分級(jí)標(biāo)準(zhǔn)，在今后的研究中可進(jìn)一步分析通過 NDVI值域變化判斷藻類水華的時(shí)空分布特征；以及進(jìn)一步分析水華暴發(fā)與當(dāng)?shù)貧庀髼l件、經(jīng)濟(jì)社會(huì)相關(guān)因子的關(guān)系（王菁晗等，2018；魯韋坤等，2017；歐陽(yáng)志宏等，2015；盛海燕等，2015）。目前基于MODIS數(shù)據(jù)建立相應(yīng)的藍(lán)藻預(yù)警預(yù)測(cè)模型已有相關(guān)報(bào)道（魏清宇等，2008），隨著衛(wèi)星遙感技術(shù)的發(fā)展，傳感器性能的進(jìn)一步提高，直接通過遙感手段開展湖泊水質(zhì)機(jī)理層面的分析以及預(yù)測(cè)預(yù)警，將是今后重要的研究方向。

圖7 滇池輕度和無(wú)水華區(qū)域（a）中度和重度水華區(qū)域（b）所占面積比例圖Fig. 7 Monthly variation of water bloom area ratio in Dianchi

圖8 滇池水華的旬變化空間分布圖Fig. 8 Spatial distribution of water bloom in every half month in Dianchi

4.2 結(jié)論

本研究基于滇池 2000－2017年時(shí)間序列的遙感影像數(shù)據(jù)，利用 GIS空間分析技術(shù)，對(duì)其近 18年來(lái) NDVI的時(shí)間變化和空間分布特征進(jìn)行了分析，探討了藍(lán)藻水華不同等級(jí)的分布面積和空間變化特征，得出以下結(jié)論：

（1）2000－2017年滇池NDVI總體呈現(xiàn)減少趨勢(shì)（趨勢(shì)系數(shù)未達(dá)到P＜0.05的顯著性檢驗(yàn)），呈減小趨勢(shì)的區(qū)域占研究區(qū)總面積55.64%，藍(lán)藻水華覆蓋度有所降低。

（2）7－8月滇池NDVI值最大，水華最嚴(yán)重，9月開始水華面積縮減，11－4月以無(wú)水華區(qū)域?yàn)橹鳎娣e比例 80%左右），5－10月以輕度水華為主（80%左右），重度水華出現(xiàn)在7－9月（8%以下）。藍(lán)藻水華的暴發(fā)時(shí)間為7月上旬和8月下旬，其重度水華面積分別高達(dá)11.4%和14.0%。

（3）滇池藍(lán)藻水華的暴發(fā)格局總體呈現(xiàn)“北重南輕”的分布態(tài)勢(shì)，水華較嚴(yán)重的區(qū)域集中在北部草海、斷橋、暉彎和滇池邊界。

（4）滇池NDVI減小表征藍(lán)藻水華覆蓋度減少，富營(yíng)養(yǎng)化程度減輕的地區(qū)主要分布在滇池北岸的斷橋，中部的觀音山西、白漁口、?？谖鞯鹊兀籒DVI增大表征藍(lán)藻水華覆蓋度增加，主要分布在滇池北岸的草海、羅家營(yíng)，滇池東南部的觀音山東、觀音山中、滇池南等地。