齊立忠 ,張蘇, 張亞平, 孫小虎
(國網經濟技術研究院有限公司,北京市 102209)
在架空輸電線路工程設計階段,快速、準確地獲取設計導線與通道中植被、建筑、道路以及交跨線路等各類地物之間的空間距離,實現輸電線路通道隱患快速檢測對于設計方案尋優、工程造價計算以及通道清障工作具有重要意義。在傳統設計中,通道隱患排查僅考慮通道地物在導線風偏面內的平面距離是否滿足對應的安全距離要求,然而當導線弧垂變化較大時導線距離地物點的最近點往往不在風偏面上,導致最近距離計算變大,存在漏檢風險。隨著電力工程設計模式從傳統二維向三維轉變,為進一步實現快速、準確的架空輸電線路工程設計階段通道隱患檢測提供了數據基礎。
為了解決不同軟件之間的電網三維設計數據交換問題,國家電網有限公司制定了統一電力工程三維設計數據標準,即電網信息模型(grid information model,GIM),GIM技術成為了電力工程三維設計未來發展的主流方向。目前,已有大量學者開展了基于GIM模型的相關應用研究[1-8],但關于輸電線路通道隱患檢測方面的相關成果較少。在現有文獻中,由于無法直接求解三維點到GIM電力線曲線的最近距離,在部分研究中采用了近似方法進行求解。文獻[9]在地物點到電力線空間曲線最近距離求解問題上提出一種近似計算方法,將xOy平面上電力線空間曲線投影到地物點平面距離最小點視為電力線空間曲線到地物點的最近點,該方法與傳統二維設計隱患檢測原理一致,未發揮出三維設計的空間特性優勢,存在隱患點漏檢情況。與此同時,有學者從另一思路出發,將導線點集化后通過計算點到離散點集最近距離實現隱患檢測。文獻[10]在基于GIM的架空輸電線路建模中,基于懸鏈線方程對電力線模型進行離散點集化,實現了三維點到電力線的最近距離計算,然而懸鏈線方程無法直接求解指定線長下的電力線位置,該方法采用傳統等水平間隔的采樣方法,采樣點分布并不均勻,影響高山區、重大交跨等特殊檔隱患檢測精度。在優化通道隱患檢測效率方面,文獻[11]提出一種基于網格化的優化算法,對地物點進行規則網格化,認為當網格內最高點滿足安全距離要求時,該網格內不存在隱患點,跳過該網格內其他地物點的隱患判斷計算,以到達加快檢測效率的目的。但對于高大樹木、建筑而言,導線可能位于側方而非上方,此時該方法易導致隱患點漏檢。
為解決現有輸電線路GIM通道隱患檢測方法存在漏檢多檢、檢測效率低下的問題,提出一種基于等線長點集化和球體包含原理的通道隱患快速檢測方法。在本文算法中,首先,采用等線長間隔采樣方式將導線由空間曲線轉化為離散點集,使離散點集在電力線上均勻分布,解決傳統等水平間隔采樣方式離散點集分布不均勻、特殊檔精度較差的問題;其次,在地物點到電力線離散點最近距離計算過程中引入球體包含原理對地物點隱患可能進行嚴密判斷,對不存在隱患可能的地物點跳過鄰近點搜索計算,實現通道隱患快速檢測。應用表明,該方法在保證計算精度的前提下達到隱患點檢測效率優化的目的,將隱患檢測耗時優化為傳統逐點檢測耗時的6.9%,對GIM的進一步應用推廣具有重要作用。
GIM文件本質上是一種分層次的電力設施屬性文件,在輸電工程中分為5個層次級別,分別為全線級、分段級、耐張段級、設備組級以及設備級,具體信息如圖1所示。

圖1 輸電GIM文件信息結構Fig.1 Transmission GIM file information structure
圖1中設備組一級定義了各耐張段內桿塔、基礎、導地線以及對應控制工況信息,其中桿塔部分記錄各桿塔模型節點坐標與掛點坐標;氣象部分記錄各個設計工況的氣象條件;導線部分對應屬性信息,如表1所示。
表1中,各項電力線參數與溫度、覆冰、風速等氣象因素結合即可計算出自重比值、冰重比載以及風偏比載,其中自重比值、冰重比載之和為垂直比載,風偏荷載為水平比載。

表1 導線屬性信息表Table 1 Wire attribute information table
由于無法直接在三維空間內求解三維點到空間曲線的最小距離,在通道隱患檢測中往往采用近似方法進行求解,易發生隱患點多檢漏檢情況。相反,三維點到離散點集的最小距離求解相對簡單,因此本文基于導線點集化思想,將導線空間曲線轉化為具有一定間隔的離散點集,通過點到點集最小距離計算實現通道隱患檢測。
1.2.1任意氣象條件下的導線三維位置計算
在目前的電力線力學研究中,通常認為自然條件下電力線形態符合懸鏈線分布,其在掛點連線豎直平面投影分布如圖2所示。

圖2 掛點連線豎直平面內電力線形態分布示意圖Fig.2 Schematic diagram of vertical plane power line distribution form of hanging point connection
圖2中:A點與B點分別為電力線兩側掛點,lAB為掛點A與掛點B在水平面上的投影距離,h為掛點A與掛點B高程差值,點P為掛點AB連線任意一點,P′為點P在電力線上對應點,lAP為掛點A與點P在水平面上的投影距離,fvP為點P對應的弧垂值,β為高差角,由掛點A、B之間高差值h與水平距離lAB計算確定。
(1)
式中:fhP為點P對應的風偏值;γv、γh分別為指定氣象條件下導線的垂直荷載與水平荷載;σ0為指定氣象條件下導線最低點的應力值。γv、γh與σ0由導線設計參數結合指定的氣象參數通過狀態方程計算確定,具體理論與步驟參照文獻[11-17]。根據任意點P三維坐標,結合弧垂值fvP、風偏值fhP以及AB連線可推算出點P對應導線點的三維坐標。
1.2.2基于等線長間隔的電力線離散點集化
在現有研究中,通常通過調整掛點A與電力線上任意點P在水平面上的投影距離lAB,每隔一定水平間距dl求解一次lAP對應的弧垂與風偏值,進而得到導線空間曲線的離散點集,離散點集分布如圖3所示。

圖3 等水平間隔采用下電力線離散點集分布圖Fig.3 Distribution diagram of power line discrete point set with equal horizontal interval
水平間隔dl的選取直接影響通道隱患檢測的檢測效果,水平間距選取過大,電力線點集過于稀疏,易導致隱患檢測漏檢;水平間距選取過小,導線點集過于密集,影響最小距離計算效率,在現有研究中水平間距dl一般選取0.10 m。
從圖3中可以看出,采用等水平間距采樣方式,存在離散點集在電力線上分布并不均勻的問題。圖3中虛線框中線長間隔dS較實線框中線長間隔明顯較長,且隨著掛點A與掛點B檔距lAB、高程差值h以及線路設計最大弧垂的增大,虛線框中線長間隔dS將進一步增長,較大可能達到水平間距dl的1.5~2.0倍。此時,虛線框附近的地物點到電力線最近點的定位極限誤差控制在0.10 m左右,不符合精度控制要求,易導致隱患點漏檢。因此,本文采用等線長間隔的采樣方法,如圖4所示使得離散點集在電力線上呈均勻分布,采用0.10 m的間隔時,將地物點到電力線最近點的定位極限誤差控制在0.05 m以內,以保證隱患檢測精度。

圖4 等線長間隔采用下電力線離散點集分布圖Fig.4 Distribution diagram of discrete point set of lower power line with equal line length interval
由于電力線在較小的局部區域可視為沿空間直線分布,基于此本文利用待計算點之前的兩個采樣點構造空間直線方程,并將該空間直線視為待計算點與其前一采樣點之間的電力線分布,通過計算空間直線在掛點水平連線上的單位投影長度實現線長間隔到水平間隔之間轉換,具體步驟如下:
步驟1:由導線設計參數結合指定的氣象參數通過狀態方程計算得到γv、γh與σ0三個參數值;

步驟3:計算掛點A處的初始區域空間直線向量,取lAB=0.2dS代入式(1)計算臨時點Pt的風偏值與弧垂值,從而推算點Pt坐標為(Xt,Yt,Zt);
步驟4:令PS=Pt,PE=A,lAP=0,i=1;


(2)
步驟7:計算采樣點坐標,令lAP=lAP+dl代入式(1)計算采樣點的風偏值與弧垂值,得到采樣點Pi坐標為(Xi,Yi,Zi)并加入到電力線離散集中;
步驟8:令PS=PE,PE=Pi,i=i+1,重復步驟5—8,直至lAP>lAB,算法結束,得到電力線離散點集。
目前,點到點集最近距離的計算問題基本是通過鄰近搜索算法解決的,然而在輸電線路通道隱患檢測中,一個塔段通道的相關地物點數量級在千萬以上,在隱患檢測中需計算每個地物點到電力線點集的最近距離,單個塔段就需要進行上千萬次鄰近搜索運算,效率低下、無法滿足隱患檢測效率需求[18-20]。
文獻[10]在優化通道隱患檢測效率方面提出一種網格法的優化算法,對地物點進行規則網格化,認為當網格內最高點滿足安全距離要求時,該網格內不存在隱患點,跳過該網格內其他地物點的隱患判斷計算,以達到加快檢測效率的目的。但對于高大樹木、建筑而言,導線可能位于側方而非上方,此時該方法易導致隱患點漏檢。本文在該方法的基礎上引入球體包含原理,采用嚴密的數學原理來對網格內地物點是否存在隱患可能性進行判定,以避免高大樹木、建筑隱患點漏檢的問題,如圖5所示。

圖5 基于球體包含原理的隱患點可能性判定示意Fig.5 Schematic diagram of potential hazard point determination based on sphere inclusion principle
如圖5所示,格網G為線路通道劃分的平面格網,曲線L為電力線,現格網G內任意方格的最高點O1,其與電力線L的距離最近點為O2,最近距離為d12,可得以點O1為球心,d12為半徑構建的實線球體C1,可得C1區域中不含有電力線。
在O1、O2所在直線上取一點O3,O3與O2之間的距離d23設置為電力線安全距離閾值,此時以點O1為球心,以d12-d23為半徑構建得到虛線球體C2。由求體包含原理可知,C2區域內任意點到C1區域外任意點的距離均大于d23,因此C2區域內的地物點均滿足安全距離要求,無隱患可能。
為了便于網格內其他相對低的地物點進行隱患可能判斷,在球體C2內取投影與平面網格一致的最大內切長方體,如圖5中實線長方體所示。P1為網格內最大平面距離點,P2為長方體與球體C2的內切點,在進行其他相對低的地物點隱患可能性判斷時,僅需要與內切點P2的高程進行比較,當地物點高程大于P2時,該點位于內切長方體內,無隱患可能。以單個塔段為例,具體步驟如下:
步驟1:坐標變換,將輸電線路通道地物點由原有坐標(x,y,z)轉化為線路走廊坐標(u,v,z),線路走廊平面坐標系如圖6所示。點A為小號桿塔中心樁、B為大號桿塔中心樁,P為AB通道內任意一點,u為AP連線在AB上的投影距離,v為點P到與AB平行且在線路前方向右側50 m處直線的垂直距離,u、v的具體計算參照三角形原理。

圖6 線路走廊平面坐標系示意圖Fig.6 Schematic diagram of line corridor plane coordinate system
步驟2:將地物點按高程z值從高到低進行排序,并按地物點類別進行歸類,得到各類地物點的點集集合{PC1,PC2,…,PCn},n為地物點的類別格式。
步驟3:令j=1。
步驟4:點集PCj含有mj個點,對應種類的地物點到電力線的安全距離閾值為τj。
步驟5:令k=1,建立網格對應的高程矩陣H,其行數為100,列數為LAB+1,LAB為A、B之間的平面距離,由于地球上任意點高程值不超過10 000 m,因此矩陣H內各元素初始設置為10 000。
步驟6:隱患可能性判斷,點Pk為集合PCj中第k個點,根據Pk的坐標(uk,vk,zk),其所在方格的行號為[vk],列號為[uk],其對應的高程限差為zmin=H[vk]×[uk]。若zk>zmin,則該點無隱患可能,進行步驟10;反之則進入步驟7進行隱患檢測。
步驟7:隱患點判定,通過點到電力線離散點集的最近點搜索,計算出點Pk到電力線的最近距離為dk。若dk≤τj,則點Pk為隱患點,加入隱患點集,進行步驟10;反之進行步驟8來更新網格對應的高程限差。
步驟8:計算點Pk到其所在方格四邊的最大平面距離ηk,計算公式如下:

(3)

步驟10:令k=k+1,若k≤mj重復步驟6;反之進行步驟11;
步驟11:令j=j+1,若j≤n重復步驟3;反之則算法運行結束。
算法流程如圖7所示。

圖7 通道隱患快速檢測算法流程圖Fig.7 Flow chart of channel hidden danger fast detection algorithm
為了驗證上述算法的有效性,本文對算法中電力線離散點集化、通道隱患檢測兩個階段進行了對比實驗。
由于采樣點之間存在一定的采樣間隔,采用計算地物點到離散點集最近點距離的方式來解決點到電力線曲線最近距離求解問題,電力線上到地物點的最近點不可避免會偏移到其在離散點集中的最近采樣點,產生最近點定位偏差。
為了驗證等水平間隔采樣方式與等線長間隔采樣方式對最近點定位偏差的影響,本文采用檔距、掛點高差以及設計冰厚均不一致的兩檔220 kV GIM電力線模型數據進行對比分析,并按照0.1 m的采樣間隔分別按照等水平間隔、等線長間隔的采樣方式進行采樣得到電力線離散點集,GIM電力線詳細信息與采樣點數如表2所示。

表2 兩檔GIM電力線相關信息與采樣點數Table 2 Relevant information and sampling points of two GIM power lines
從表1中可以看出,電力線1跨越檔距較小、掛點高程相近且位于輕冰區,其設計線長與水平檔距大致一致,用以表征普通輸電線路的電力線分布;電力線2跨越檔距較大、掛點高程懸殊且位于重冰區,其設計線長與水平檔距存在較大差距,用以表征高山區、重大交跨等特殊檔的電力線分布。
本文采用密集采樣的方式來統計電力線各區域的最近點定位偏差,每隔0.001 m等水平間隔在電力線上取點視為地物點的最近點,分別計算其在等水平間隔、等線長間隔離散點集中的最近點定位偏差。為了便于展示分析,本文每5 m水平距離區間內統計一個最近點定位偏差最大值作為區域最近點定位限差,兩檔GIM電力線的區域最近點定位限差隨水平距離分布如圖8所示。
從表2、圖8中可以得出,等線長間距采樣方式在表征普通輸電線路的電力線分布的電力線1上采樣點數、區域最近點定位限差與等水平間距采樣基本保持一致,但在表征高山區、重大交跨等特殊檔電力線分布的電力線2上表現出明顯優勢。如圖8(b)所示,等水平間隔采樣方式在電力線2上區域最近點定

圖8 兩檔GIM電力線的最近點定位限差分布圖Fig.8 Distribution diagram of maximum value of nearest point positioning deviation of two GIM power lines
位限差達到0.08 m左右,然而等線長采樣方式仍能夠保持在0.05 m左右。
為了進一步分析最近點點位偏差的分布規律,在電力線2上截取了等水平間隔采樣方式下定位限差最大與最小的兩個區域最近點定位偏差數據,如圖9所示。

圖9 區域最近點定位偏差分布圖Fig.9 Distribution map of regional nearest point positioning deviation
從圖9中可以得出,等水平間隔采樣與等線長采樣兩種方式最近點定位偏差均符合0到區域最近點定位限差之間均勻分布,而區域最近點定位限差由區域采樣點之間的線長間隔決定。從圖9(a)與圖9(b)對比中可以看出,對于跨越檔距較大、掛點高程懸殊的特殊檔等水平間隔的采樣點在電力線上分布并不均勻,局部區域的線長間隔差異大,導致區域最近點定位限差達到設計的1.5~2.0倍,不利于精度控制。而等線長間距采樣方式無論在檔距較小、掛點高程相近的普通檔以及跨越檔距較大、掛點高程懸殊的特殊檔均能夠保證采樣點分布均與,將最近點定位限差控制在目標范圍內。
為驗證本文通道檢測算法的應用效果,選取某一220 kV架空輸電線路GIM模型中連續50檔數據進行通道隱患檢測,地物點數據通過激光雷達系統進行采集,各檔地物點統計如圖10所示。

圖10 50檔地物點點數統計Fig.10 Statistical chart of 50 ground feature points
基于GIM模型電力線采用等線長間隔采樣后得到電力線離散點集后,分別按照傳統逐點檢測、網格法過濾以及本文方法進行隱患檢測,分別對檢測速度以及檢測結果進行對比分析,如圖11所示。
從圖11中可以看出,網格法與本文方法在檢測效率上較傳統逐點檢測均有較大提升,各檔隱患檢測耗時縮短為原本耗時的1/20到1/8。根據統計,50檔通道檢測逐點檢測耗時1 169 s,而網格法耗時94 s,為逐點檢測耗時的8.0%;本文方法耗時80 s,為傳統逐點檢測耗時的6.9%。

圖11 3種隱患檢測方法檢測耗時對比Fig.11 Comparison of detection time of three hidden danger detection methods
表3為3種通道隱患檢測方法檢測出的隱患點數統計情況,表中標紅部分的網格法檢測結果與其余兩種方法檢測結果并不一致,存在漏點情況。而本文方法從表3中可以看出,各檔檢測結果仍與逐點檢測保持一致,在大大提升檢測效率的同時,保證檢測結果正確性。

表3 通道隱患檢測隱患點數統計Table 3 Statistics of hidden danger points of channel hidden danger detection
為進一步闡述本文方法與網格法檢測的差異性原因,截取了本文方法與網格法隱患檢測結果不同處的點云場景,如圖12所示。
圖12中白色點為電力線離散點、藍色點為地面點、綠色點為植被點、紅色點為通道隱患點。由俯視角、側視角對比可知,當導線從高大樹木側方通過時,由于側方地物高點滿足安全距離要求,網格法認為地物高點所在網格無隱患可能,跳過最近距離檢測導致漏檢。而本文提出的基于球體包含原理采用更為嚴密的隱患判斷依據,針對側方高大地物情況仍能夠保證正確的檢測結果。

圖12 本文方法與網格法隱患檢測結果對比Fig.12 Comparison of hidden danger detection results between this method and grid method
本文提出了一種基于等線長點集化和球體包含原理的通道隱患快速檢測方法,在保證隱患檢測精度的基礎上實現通道隱患檢測速度。
1)提出了基于等線長間隔的GIM電力線離散點集方法,將電力線局部區域視為空間直線分布實現線長間隔到水平間隔的改正計算,克服現有等水平間隔采樣離散點集在電力線上分布并不均勻,跨越檔距較大、掛點高程懸殊特殊檔最近點定位偏差增大的問題。
2)針對網格法檢測速度優化算法存在漏檢問題,在該方法的基礎上引入球體包含原理,采用嚴密的數學原理來對網格內地物點是否存在隱患可能性進行判定,以避免高大樹木、建筑隱患點漏檢的問題。
3)經實驗驗證,本文提出的等線長間隔采樣方法無論在檔距較小、掛點高程相近的普通檔以及跨越檔距較大、掛點高程懸殊的特殊檔均能夠保證采樣點分布均勻,將最近點定位限差控制在目標范圍內;基于球體包含原理的隱患快速檢測方法能夠在保證檢測精度的前提下將隱患檢測耗時優化為原本耗時的6.9%。