基于Apriori算法的云南省雷電活動與地表覆蓋類型關(guān)聯(lián)性分析

馬儀，潘浩，周仿榮，何錦強(qiáng)，李銳海，廖永力

(1. 電力遙感技術(shù)聯(lián)合實驗室(云南電網(wǎng)有限責(zé)任公司電力科學(xué)研究院)，昆明 650217；2. 南方電網(wǎng)科學(xué)研究院，廣州 510663)

0 引言

隨著我國電網(wǎng)朝著高電壓等級、大輸送容量的方向發(fā)展，架空線路的桿塔高度和走廊寬度逐漸增加，架空線路遭受雷擊的風(fēng)險日益增大。統(tǒng)計表明，雷害是架空線路跳閘的最主要誘因之一，南方電網(wǎng)區(qū)域110 kV及以上電壓等級線路雷擊跳閘占總跳閘次數(shù)的比例高達(dá)60.5%[1 - 2]。雷雨天氣多發(fā)于每年的4月至10月，雷電活動的強(qiáng)度受氣候、地形地貌等多種因素影響[3 - 6]。

目前，國內(nèi)外學(xué)者對于雷電活動與地表覆蓋類型等地形地貌因素的相關(guān)性研究大多基于數(shù)理統(tǒng)計法和回歸分析方法，Orville等對美國近10年的雷電活動數(shù)據(jù)進(jìn)行了統(tǒng)計，得到了地閃數(shù)據(jù)隨年、月和日分布的柱狀圖[7]。康鵬等采用數(shù)理統(tǒng)計方法來探究青藏鐵路輸電線路近50年的雷電活動特點，認(rèn)為不同地形或緯度地區(qū)的地閃活動有一定差異[8]。趙生昊等采用ArcGIS空間分析和經(jīng)典數(shù)理統(tǒng)計方法來探究重慶市2007—2016年雷電活動與海拔、坡向和坡度之間的相關(guān)性，認(rèn)為雷電流幅值與海拔、坡度之間沒有明顯關(guān)聯(lián)性，而隨坡向變化有較大的差異[9]。余建華等采用關(guān)聯(lián)性挖掘方法來探究江西省雷電流幅值均值、地閃密度極大值與地表覆蓋類型之間的相關(guān)性[10]。Ezcurra等選用5 km×5 km的網(wǎng)格為統(tǒng)計單元對西班牙巴斯克地區(qū)的地閃密度和雷電日進(jìn)行了分析，得到該區(qū)域的雷電隨著地形地貌的分布特征規(guī)律[11]。以上研究表明一個地區(qū)雷電活動一定程度上依賴于當(dāng)?shù)氐牡乩憝h(huán)境，且與多個因素息息相關(guān)，但使用數(shù)理統(tǒng)計與回歸方法難以得到雷擊與地表覆蓋類型等諸多因素之間的定量關(guān)系。V.Bourscheidt等[12]比較了巴西南部地閃密度以及海拔、坡度的參數(shù)之間的關(guān)系，認(rèn)為該區(qū)域閃電密度與海拔高度無關(guān)，但與坡度之間存在明顯的線性關(guān)系。我國學(xué)者司馬文霞、李永福等[13]在2011—2012年通過改進(jìn)的網(wǎng)格法對地閃密度和雷電流幅值進(jìn)行了統(tǒng)計分析，優(yōu)化了傳統(tǒng)網(wǎng)格法存在的空間分辨率不夠和部分雷電參數(shù)物理意義不明確的問題。文獻(xiàn)[14]以浙江省為例對雷擊的地閃密度與氣象、地形、地貌等要素進(jìn)行Pearson相關(guān)系數(shù)分析，結(jié)果表明地閃密度與月均氣溫、降水量、濕度、風(fēng)速間具有較高的相關(guān)性，但并未得出定量的關(guān)聯(lián)結(jié)果。

近年來，數(shù)據(jù)挖掘技術(shù)廣泛應(yīng)用[15 - 17]，包括神經(jīng)網(wǎng)絡(luò)、依賴性分析、遺傳算法、關(guān)聯(lián)規(guī)則等。部分學(xué)者嘗試基于關(guān)聯(lián)規(guī)則算法探究雷電活動參數(shù)與地形地貌的相關(guān)性[18 - 22]。吳巍巍等用聚類算法分析得到雷電地閃活動的聚集點，再結(jié)合Alpha Shapes輪廓識別算法分析雷電運動走廊，發(fā)現(xiàn)雷電地閃容易聚集于平地、平原或河流等地勢較低的區(qū)域，并有沿著河流走向運動的趨勢[23]。潘健利用聚類算法來探究江蘇省雷電流幅值分布特點與地表覆蓋類型之間的相關(guān)性，認(rèn)為林地與大雷電流幅值具有較大相關(guān)性[24]。喬飛等利用神經(jīng)網(wǎng)絡(luò)技術(shù)模擬輸電走廊沿線地閃密度分布情況，因未充分考慮地理地形、土地類型、季節(jié)等因素預(yù)測模型效果不佳[20]。Apriori關(guān)聯(lián)性算法可以從海量的數(shù)據(jù)中，挖掘出具有研究價值、有代表性的數(shù)據(jù)，被挖掘出的數(shù)據(jù)之間存在著一定的關(guān)聯(lián)關(guān)系。其特性可以滿足雷電活動數(shù)據(jù)與地表覆蓋數(shù)據(jù)關(guān)聯(lián)性的分析。

目前關(guān)于地表覆蓋類型與雷電活動關(guān)聯(lián)關(guān)系的研究尚未成熟，存在較多爭議。云南省地表覆蓋類型極具多樣化，本文基于云南省2010—2019年雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)和地表覆蓋類型數(shù)據(jù)，采用統(tǒng)計方法和Apriori算法對地表覆蓋類型與正閃比、地閃密度、雷電流幅值、雷電流幅值累積概率分布等雷電活動參數(shù)的相關(guān)性進(jìn)行分析，結(jié)果對輸電線路建設(shè)、桿塔選址、防雷改造等工作具有指導(dǎo)意義。

1 地表覆蓋類型與雷電活動的統(tǒng)計特性

1.1 數(shù)據(jù)來源

云南省已建成覆蓋全省的雷電定位系統(tǒng)(lightning location system，LLS)，并與相鄰省份實現(xiàn)數(shù)據(jù)互聯(lián)，雷電探測效率超過90%[21]?？紤]到太陽黑子活動周期一般為10 a，為減小太陽活動引起的氣候變化對雷電活動的影響，本文選取云南地區(qū)2010—2019年的雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)來進(jìn)行雷電活動參數(shù)的統(tǒng)計與分析，共17 254 429條雷擊數(shù)據(jù)，包括雷擊發(fā)生的時間，緯度、經(jīng)度、雷電流幅值以及回?fù)舻葏?shù)。

本文選用中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心發(fā)布的2020年版地表覆蓋類型數(shù)據(jù)。該版數(shù)據(jù)是基于Landsat-8遙感影像，將地球表面陸地分為6大類和25小類，具體如表1所示。

表1 地表覆蓋類型Tab.1 Types of land cover

云南省的地表覆蓋類型分布情況如圖1所示，可以看到云南省包含大部分地表覆蓋類型，僅不含灘涂、沙地、戈壁和鹽堿地。忽略大類“未利用土地”中的小類“其他”，進(jìn)一步分類處理，得到云南省地表覆蓋面積的柱形圖，如圖2所示。云南省地表覆蓋類型占地面積最大的是林地，其中有林地、灌木林面積較大，均大于80 000 km2。云南省占地面積最小的地表覆蓋類型是裸土地和沼澤地，兩者占地面積均只有50 km2左右。

圖1 云南省地表覆蓋類型Fig.1 Types of land cover in Yunnan Province

圖2 云南省土地覆蓋類型面積統(tǒng)計圖Fig.2 Statistical areas of land cover types in Yunnan Province

1.2 統(tǒng)計結(jié)果

1.2.1 地閃強(qiáng)度

地閃指云內(nèi)荷電中心與大地和地物之間的放電過程。地閃強(qiáng)度指地閃發(fā)生時監(jiān)測到的雷電流幅值。經(jīng)統(tǒng)計，正地閃次數(shù)遠(yuǎn)少于負(fù)地閃次數(shù)，平均只占總地閃次數(shù)的5%左右。按正、負(fù)地閃極性分別計算各地表覆蓋類型區(qū)域地閃的平均雷電流幅值。

負(fù)地閃強(qiáng)度的統(tǒng)計結(jié)果如圖3左邊所示，可以看出，云南省負(fù)地閃的雷電流幅值分布的方差較小。其中，中覆蓋度草地和裸巖石質(zhì)地的雷電流幅值絕對值最大，均超過30 kA；沼澤地和裸土地的雷電流幅值絕對值最小，均接近20 kA。

正地閃強(qiáng)度統(tǒng)計結(jié)果如圖3右邊所示，可以看出，不同地表覆蓋類型下正地閃強(qiáng)度的差異較大，且正地閃的雷電流幅值比負(fù)地閃的雷電流幅值平均值更大。其中永久性冰川雪地和低覆蓋度草地的正雷電流幅值均高于60 kA；雷電流幅值最低的灘地，其正雷電流幅值值也達(dá)到了30 kA。

圖3 正、負(fù)地閃強(qiáng)度Fig.3 Positive and negative ground flash intensity

對雷電流幅值累積概率分布進(jìn)行正態(tài)分布擬合，大致上服從式(1)，使用Matlab軟件采用最小二乘法，對圖4曲線進(jìn)行擬合來求出參數(shù)a和參數(shù)b?？紤]到曲線擬合需要大量數(shù)據(jù)來支撐，一級地表覆蓋類型數(shù)據(jù)量遠(yuǎn)大于二級地表覆蓋類型，為保證擬合結(jié)果的精確性，本文采用采用6種一級地表覆蓋類型進(jìn)行曲線擬合，擬合結(jié)果表2所示。

表2 地表覆蓋類型對應(yīng)雷電流幅值a和概率b分布參數(shù)Tab.2 Probability distribution parameter of lightning current amplitude of land cover types

圖4 雷電雷電流幅值/kA流幅值雷擊概率曲線擬合Fig.4 Curve fittings of lightning strike probability of lightning current amplitude

(1)

式中：I為雷電流幅值，kA；P(>I)為雷電流幅值超過I的概率；a為待定參數(shù)，物理意義為中值電流，即p(>a)=0.5；b為待定參數(shù)，物理意義為曲線陡度，b值越大則曲線越陡、衰減越快。

由表2可見，草地雷電流幅值分布較集中，方差較大，雷電流幅值累積概率分布曲線減小得最快；未利用土地雷電流幅值的累計概率曲線減小得最慢，雷電流幅值方差較小，分布更分散。城鄉(xiāng)、工礦、居民用地的中值電流最小，為20.318 7 kA；草地的中值電流最大，為23.235 8 kA。

由于草地出現(xiàn)高雷電流幅值的可能性是最高的，故發(fā)生雷擊跳閘事故的可能性是最高的。因此在防雷設(shè)計中，建議著重關(guān)注穿越草地的輸電線路桿塔防雷措施，以此來減少線路的雷擊跳閘率。

1.2.2 地閃密度

地閃密度指某區(qū)域每平方公里每年地面落雷次數(shù)，單位為次/(km2·a)。對于每種地表覆蓋類型區(qū)域，分別統(tǒng)計落在該區(qū)域的雷擊次數(shù)，進(jìn)而計算地閃密度，統(tǒng)計結(jié)果如圖5所示。城鎮(zhèn)的地閃密度是最大的，超過了7次/(km2·a)，其次較大的有裸土地、其他林地、灌木林等。永久性冰川雪地和裸巖石質(zhì)地的地閃密度最低，均小于2次/(km2·a)。

圖5 地閃密度分布情況Fig.5 Distribution of lightning density

1.2.3 正閃比

正閃比指正極性雷數(shù)量占所有雷擊數(shù)量的比例。正閃比例分布情況的統(tǒng)計結(jié)果如圖6所示，正閃比最高的地表覆蓋類型是裸土地，達(dá)到6.60%，其次為其他林地，達(dá)到6.26%。正閃比較低的地表覆蓋類型是沼澤地、永久性冰川雪地、裸巖石質(zhì)地、城鎮(zhèn)用地、中覆蓋草地等。由于裸土地、其他林地的正閃比較大，因此在裸土地、其他林地附近的輸電線路設(shè)計和安裝，建議考慮正閃電對此地區(qū)線路防雷的影響。

圖6 正閃比分布情況Fig.6 Distribution of positive flash ratio

2 地表覆蓋類型與雷電活動關(guān)聯(lián)性分析

2.1 算法簡介

本文利用Apriori算法對雷電數(shù)據(jù)和地表覆蓋類型數(shù)據(jù)進(jìn)行深度挖掘，分析雷電參數(shù)與地表覆蓋類型間的關(guān)聯(lián)性，并與統(tǒng)計結(jié)果相互對比驗證，進(jìn)一步說明地表覆蓋類型對雷電活動的影響。

Apriori算法是一種經(jīng)典的基于關(guān)聯(lián)規(guī)則的挖掘算法，其主要思想是：通過尋找頻繁項集，從海量的數(shù)據(jù)中，挖掘出具有研究價值、有代表性的數(shù)據(jù)，被挖掘出的數(shù)據(jù)之間存在著一定的關(guān)聯(lián)關(guān)系。算法中涉及的具體概念介紹如下。

1)項與項集

項是數(shù)據(jù)集合中的基本元素，項集是項的集合。項集一般用L表示，滿足最小支持度的項集稱為頻繁項集。

2)事務(wù)

事務(wù)是指在給定的數(shù)據(jù)集合中有相同的數(shù)據(jù)。

3)支持度與置信度

項集A、B同時發(fā)生的概率P(A∪B)稱為關(guān)聯(lián)規(guī)則的支持度(也稱相對支持度)dSupport(A∪B)， 如式(2)所示。

dSupport(A∪B)=P(A∪B)

(2)

項集A發(fā)生，則項集B發(fā)生的概率P(A|B)為關(guān)聯(lián)規(guī)則的置信度dConfidence(A?B)， 如式(3)所示。

dConfidence(A?B)=P(A|B)

(3)

4)最小支持度和最小置信度

最小支持度min_sup與最小置信度min_conf是用戶或者專家定義的衡量支持度的一個閾值，表示項目集在統(tǒng)計意義上的最低重要性與關(guān)聯(lián)規(guī)則的最低可靠性。

dSupport(A→B)≥dmin_sup

(4)

dConfidence(A→B)≥dmin_conf

(5)

5)強(qiáng)關(guān)聯(lián)規(guī)則

對于同時滿足用戶或?qū)＜翌A(yù)設(shè)的最小支持度與最小置信度的關(guān)聯(lián)規(guī)則，稱之為強(qiáng)關(guān)聯(lián)規(guī)則。我們一般感興趣的就是強(qiáng)關(guān)聯(lián)規(guī)則。

Apriori算法的挖掘過程主要包含4個步驟[22 - 23]。

1)掃描所有的事務(wù)，獲得出現(xiàn)過的每個項，對于每個項都單獨進(jìn)行數(shù)量統(tǒng)計。根據(jù)人為設(shè)定的最小支持度選擇是否剔除，從而生成頻繁1項集L1。

2)對L1進(jìn)行連接步驟，產(chǎn)生候補(bǔ)2項集的集合C2，同理，掃描數(shù)據(jù)庫中的所有事物，對C2每個項集進(jìn)行單獨計數(shù)，再根據(jù)最小支持度從C2中刪除不滿足要求的項集，從而獲得頻繁2項集L2。

3)同樣的，再產(chǎn)生頻繁3項集L3。

4)以此類推，對Lk-1經(jīng)連接步驟產(chǎn)生的集合執(zhí)行剪枝，產(chǎn)生候補(bǔ)k項集Ck，然后掃描所有事務(wù)，對每個項集進(jìn)行統(tǒng)計數(shù)量，再根據(jù)最小支持度執(zhí)行剔除策略，獲得頻繁k項集Lk。這樣，算法迭代完成后就能獲得所有的頻繁項集。算法的流程圖如圖7所示。

圖7 Apriori 算法流程圖Fig.7 Flow chart of Apriori algorithm

利用一個包含10條事務(wù)數(shù)的簡單數(shù)據(jù)庫為例進(jìn)行分析計算，說明算法的原理。算例數(shù)據(jù)庫如表3所示，包括編號1—10的10條事務(wù)數(shù)據(jù)，a、b、c、d、e、f這6種不同屬性，每條事務(wù)包含2種到5種不同的屬性，設(shè)定最小支持度為0.4，最小置信度為0.7。

表3 算例數(shù)據(jù)庫Tab.3 Example database

首先找出1項候選集如表4所示，剔除不滿足最小支持度的候選集，得到頻繁1項集如表5所示。

表4 1項候選集Tab.4 1-candidate set

表5 頻繁1項集Tab.5 Frequent 1-candidate set

由頻繁1項集組合得到2項候選集如表6所示,剔除不滿足最小支持度的候選集，得到頻繁2項集如表7所示。

表6 2項候選集Tab.6 Two-candidates set

表7 頻繁2項集Tab.7 Frequent two-candidates set

以此類推，可以找出所有頻繁k項集，進(jìn)而計算頻繁項集的置信度，篩選得到我們感興趣的關(guān)聯(lián)規(guī)則。

2.2 基于Apriori算法的數(shù)據(jù)庫建立

Apriori算法的運行是事務(wù)型數(shù)據(jù)庫為基礎(chǔ)的，首先利用網(wǎng)格法將云南省劃分為大小相等的矩形網(wǎng)格，并對網(wǎng)格進(jìn)行編號，統(tǒng)計每個網(wǎng)格的正閃比、雷電流幅值和地閃密度。本文選取的網(wǎng)格均為1 km×1 km大小，既有較高的精度，又可以保證與地表覆蓋類型有較好的對應(yīng)關(guān)系。具體做法如下。

1)對柵格數(shù)據(jù)進(jìn)行“柵格轉(zhuǎn)點”操作，使用了ArcGIS自帶的柵格轉(zhuǎn)點工具來完成，把柵格數(shù)據(jù)轉(zhuǎn)換成了一個個矢量點。每個矢量點的位置是對應(yīng)柵格的中點。

2)每個柵格的編號取矢量點數(shù)據(jù)的屬性值，再利用ArcGIS自帶的點轉(zhuǎn)柵格工具，再把矢量點轉(zhuǎn)變?yōu)闁鸥駭?shù)據(jù)。

3)再次用ArcGIS自帶的系統(tǒng)工具“extract multiple values to points”，對每條雷電數(shù)據(jù)進(jìn)行取屬性值操作，即得到了初始數(shù)據(jù)庫，示例如圖8所示。

圖8 初始數(shù)據(jù)庫示例Fig.8 Initial database example

該數(shù)據(jù)庫建立時注意新增一個名為“柵格編號”的屬性，表明該雷擊位置所處柵格的柵格編號。

輸入Access(Microsoft Office Access)軟件是由微軟公司發(fā)布的關(guān)系數(shù)據(jù)庫管理系統(tǒng)，將數(shù)據(jù)存儲于基于Microsoft Jet Database 數(shù)據(jù)庫里，還可以直接導(dǎo)入或者連接數(shù)據(jù))中進(jìn)行分類統(tǒng)計，求得每個柵格的正閃比、地閃密度、雷電流幅值3個參數(shù)。再拿該表與原數(shù)據(jù)庫進(jìn)行連接查詢，得到每起雷擊發(fā)生位置所處柵格的正閃比、地閃密度和雷電流幅值。至此，初始的Apriori算法的數(shù)據(jù)庫就建立完成了。每一條雷擊數(shù)據(jù)就代表一個事務(wù)。

Apriori算法要求輸入的數(shù)據(jù)為離散化數(shù)據(jù)，離散化的數(shù)據(jù)有助于得到有效的布爾型關(guān)聯(lián)規(guī)則。ksdensity函數(shù)可用于計算一維或二維核密度或分布估計，本文利用ksdensity函數(shù)對正閃比、地閃密度、地閃強(qiáng)度數(shù)據(jù)的分布情況進(jìn)行分析，三者均近似呈正態(tài)分布，采用分級法把正閃比、地閃密度和地閃強(qiáng)度分別概化為低、中、高3個等級，具體分級標(biāo)準(zhǔn)如表8所示。

表8 正閃比、地閃密度、地閃強(qiáng)度數(shù)據(jù)概化Tab.8 Generalization of positive flash ratio, ground flash density and ground flash intensity

將20種小類地表覆蓋類型編號為1～20，將正閃比_低、正閃比_中和正閃比_高分別編號為21、22和23，同理也對地閃密度3個等級編號為24、25和26；對雷電流幅值3個等級編號為27、28和29。完成數(shù)據(jù)預(yù)處理后得到能被算法直接調(diào)用的事務(wù)數(shù)據(jù)庫，部分?jǐn)?shù)據(jù)如表9所示。

表9 Apriori事務(wù)數(shù)據(jù)庫部分?jǐn)?shù)據(jù)示例Tab.9 Examples of transaction database in Apriori

表9中每行數(shù)據(jù)都是一個事務(wù)，每條數(shù)據(jù)包含4個屬性：地表覆蓋類型、正閃比、地閃密度和地閃強(qiáng)度。

2.3 地表覆蓋類型與雷電活動關(guān)聯(lián)規(guī)則挖掘

利用關(guān)聯(lián)分析模型開展關(guān)聯(lián)規(guī)則挖掘，由于云南省地表覆蓋類型中城鎮(zhèn)用地、沼澤地等的占地面積較小，故為了防止在算法迭代過程中被剔除，本文設(shè)定最小支持度為5%，最小置信度為45%，篩選得到的結(jié)果如表10所示，如2.1節(jié)所述，規(guī)則前件與規(guī)則后件同時發(fā)生的概率稱為關(guān)聯(lián)規(guī)則的支持度；規(guī)則前件發(fā)生，則規(guī)則后件發(fā)生的概率為關(guān)聯(lián)規(guī)則的置信度。

表10 關(guān)聯(lián)規(guī)則結(jié)果Tab.10 Association rule results

強(qiáng)關(guān)聯(lián)規(guī)則為：1)城鎮(zhèn)用地→地閃密度_高；2)其他林地→正閃比_高；3)灌木林→地閃密度_高；4)中覆蓋度草地→正閃比_低。也就是說，對規(guī)則進(jìn)一步進(jìn)行解讀：城鎮(zhèn)用地的地閃密度比其他地表覆蓋類型的地閃密度高的多，應(yīng)該注意防范城市雷擊災(zāi)害；其他林地區(qū)域普年正地閃發(fā)生概率較大，在其他林地的輸電桿塔應(yīng)注意防范正極性雷擊；灌木林的地閃密度也普遍較高；中覆蓋草地的正閃比顯著低于其他幾種地表覆蓋類型；地表覆蓋類型和地閃強(qiáng)度的關(guān)聯(lián)性并不大，這是由于各個地表覆蓋類型的雷電流幅值均值分布較為集中。

從表10可以發(fā)現(xiàn)，林地地閃密度容易較高。這主要是因為尖端放電原理，在強(qiáng)電場作用下，物體曲率大的地方附近，等電位面密，電場強(qiáng)度劇增，致使這里空氣被電離而產(chǎn)生氣體放電現(xiàn)象，稱為電暈放電。而尖端放電為電暈放電的一種，專指尖端附近空氣電離而產(chǎn)生氣體放電的現(xiàn)象。當(dāng)雷雨云過境時，云的中下部是強(qiáng)大負(fù)電荷中心，云的下墊面是正電荷中心，于是在云與地面間形成強(qiáng)電場。在樹木、山頂草、林木、巖石等尖端附近，等電位面就會很密集，這里電場強(qiáng)度極大，空氣發(fā)生電離，因而形成從地表向大氣的尖端放電。

城鎮(zhèn)用地的地閃密度容易較高，其主要受3點因素影響：1)熱島效應(yīng)。在城市熱島效應(yīng)下，城市中心空氣不斷升溫，城中區(qū)的氣壓降低，導(dǎo)致周圍郊區(qū)冷空氣聚集到城市來補(bǔ)足壓強(qiáng)差。城中區(qū)的熱空氣上升到一定高度后會下沉到周圍的郊區(qū)，來補(bǔ)充郊區(qū)的氣流，從而形成了一種氣流循環(huán)，即熱島循環(huán)。熱島循環(huán)使大氣結(jié)構(gòu)極其不穩(wěn)定，增大了對流產(chǎn)生的概率和強(qiáng)度，容易形成對流云和降水[25]；2)城市污染。人類活動使得城市污染加劇，城市空氣中的塵埃顆粒是遠(yuǎn)多于郊區(qū)的。這些顆粒進(jìn)入大氣層，有利于水蒸氣的凝結(jié)，并且使云層更加容易攜帶電荷而產(chǎn)生雷云；3)同林地分析類似，城市有許多高聳構(gòu)筑物，雷電具有先導(dǎo)特性，城區(qū)中的高層建筑、通信塔等高聳構(gòu)筑物相較周圍有很強(qiáng)的引雷作用。

3 結(jié)論

基于云南省地表覆蓋類型數(shù)據(jù)以及多年累積雷電定位系統(tǒng)監(jiān)測數(shù)據(jù)，結(jié)合統(tǒng)計方法和Apriori算法對地表覆蓋類型與正閃比、地閃密度、地閃強(qiáng)度的相關(guān)性進(jìn)行分析。主要結(jié)論如下。

1)中覆蓋度草地和裸巖石質(zhì)地的負(fù)地閃強(qiáng)度較大，雷電流幅值均值均超過30 kA；冰川雪地和低覆蓋度草地的正地閃強(qiáng)度很大，均值均高于60 kA。應(yīng)注意這些區(qū)域線路的防雷工作。

2)由于熱島效應(yīng)、人類活動以及城區(qū)高聳構(gòu)筑物的影響，城鎮(zhèn)用地的正、負(fù)地閃密度均比其他地區(qū)要高。

3)林地的正閃比例較高，因此林地區(qū)域輸電線路的雷電防護(hù)措施需要考慮正極性雷電的影響。