趙飛飛，張顯云

（1.貴州大學(xué) 礦業(yè)學(xué)院，貴州 貴陽 550025；2.貴州省測繪資料檔案館，貴州 貴陽 550004；3.貴州省自然資源衛(wèi)星應(yīng)用中心，貴州 貴陽 550004）

0 引 言

隨著“碳達(dá)峰、碳中和”戰(zhàn)略的提出，可再生能源必然會在我國經(jīng)濟(jì)發(fā)展中扮演著越來越重要的角色。風(fēng)能作為一種清潔、無公害的可再生能源，早已被各國所青睞。但由于長期運(yùn)作和外力影響，風(fēng)力發(fā)電機(jī)塔筒極可能發(fā)生彎曲、傾斜，甚至倒塌。垂直度檢測作為風(fēng)力發(fā)電機(jī)日常運(yùn)維的一項(xiàng)重要指標(biāo)，決定著機(jī)器的正常運(yùn)行。如果能精確、快速地檢測出風(fēng)電塔垂直度，便能高效監(jiān)控其安全性，大大減少損失[1-2]。

三維激光掃描技術(shù)作為一種無接觸式的測量技術(shù)，不僅具有高精度、高效率等特點(diǎn)，還能自動獲取海量點(diǎn)云，生成高精度三維模型[3]。這為風(fēng)機(jī)塔筒垂直度檢測提供了很好的技術(shù)支撐。目前，三維激光技術(shù)廣泛應(yīng)用于構(gòu)筑物監(jiān)測等領(lǐng)域。史陽軍、?；⒘諿1-2]等人均采用三維激光技術(shù)檢測風(fēng)電塔垂直度，并將結(jié)果與全站儀測量結(jié)果做對比，驗(yàn)證了該方法能滿足工作需求。聞亞[4]提出采用三維激光技術(shù)將古亭進(jìn)行分層，分別提取古亭的底部、中部和頂部的中心點(diǎn)，采用最小二乘法擬合中心點(diǎn)，分析古亭的垂直度和偏移量。王巍[5]利用三維激光掃描技術(shù)分別對電力鐵塔垂直度進(jìn)行了分析。以上學(xué)者的研究均證明了三維激光掃描技術(shù)在測量中的可靠性。本文針對風(fēng)力發(fā)電機(jī)塔筒，利用三維激光掃描技術(shù)對其進(jìn)行垂直度檢測研究。

1 測區(qū)與儀器

本文研究對象是位于貴州大學(xué)西校區(qū)的風(fēng)力發(fā)電機(jī)。該風(fēng)力發(fā)電機(jī)西面靠山，周圍樹木茂密，東面為校內(nèi)交通要道，來往車輛行人偏多，測量干擾源較多，不適合激光掃描儀的架設(shè)。其南北方向較為空曠，無大型遮擋物，干擾較少，適合架設(shè)儀器?，F(xiàn)場具體情況如圖1 所示。

圖1 風(fēng)力發(fā)電機(jī)現(xiàn)場情況

2 點(diǎn)云數(shù)據(jù)采集

2.1 點(diǎn)云數(shù)據(jù)的獲取

本次實(shí)驗(yàn)采用的儀器是Riegl VZ-1000 地面三維激光掃面儀。該激光掃描儀具有掃描速度快、拼接時間短、產(chǎn)品質(zhì)量好、配套軟件多等特點(diǎn)。其掃描精度為5 mm，掃描速度為每秒30 萬點(diǎn)，視角范圍0°～360°。進(jìn)行數(shù)據(jù)掃描前，需實(shí)地查看，進(jìn)行儀器校驗(yàn)、選取測站點(diǎn)位、確定掃描方式和分辨率精度、選取控制點(diǎn)等操作。本次實(shí)驗(yàn)數(shù)據(jù)處理將采用點(diǎn)云拼接方式，所以需根據(jù)風(fēng)電發(fā)電機(jī)周圍環(huán)境，滿足兩站之間具有高重復(fù)率數(shù)據(jù)，最終設(shè)置4 個測站。每站掃描時間約為6 ～9 min，確保掃描數(shù)據(jù)的完整性。

2.2 點(diǎn)云數(shù)據(jù)處理

本次實(shí)驗(yàn)數(shù)據(jù)采用RiScan PRO 軟件進(jìn)行點(diǎn)云數(shù)據(jù)的處理。處理步驟包括點(diǎn)云拼接、點(diǎn)云去噪及點(diǎn)云切片等。

2.2.1 點(diǎn)云拼接

將多站掃描數(shù)據(jù)導(dǎo)入RiScan PRO 軟件，建立新的拼接站點(diǎn)，選擇重疊部分明顯的特征點(diǎn)，如房屋角點(diǎn)、窗子的邊線交點(diǎn)等，進(jìn)行點(diǎn)云匹配。每個測站至少選擇3 至4 個同名點(diǎn)進(jìn)行拼接。在拼接過程中，通過尋找多站之間風(fēng)機(jī)周圍各種特征較為明顯、清楚的同名點(diǎn)，多次拼接，最終采用拼接效果最佳的結(jié)果進(jìn)行后續(xù)操作。

2.2.2 點(diǎn)云去噪

獲取點(diǎn)云數(shù)據(jù)時，由于受到各種因素，如行人、車輛、建筑物及植物等的影響，采集的點(diǎn)云數(shù)據(jù)會存在大量噪聲。對于較為明顯的噪聲，直接采用手動刪除的方式進(jìn)行處理。對于靠近風(fēng)機(jī)塔筒的散亂點(diǎn)云，可采用軟件自帶的函數(shù)進(jìn)行去除。圖2 為去噪后的風(fēng)機(jī)塔筒圖。

圖2 去噪后的風(fēng)機(jī)塔筒

2.2.3 點(diǎn)云切片

平面截面法是將選取的基準(zhǔn)面作為截面，以相同高度間隔截取橫斷面，將橫斷面的中心點(diǎn)作為計算物體垂直度的特征點(diǎn)。為獲取風(fēng)力發(fā)電機(jī)塔筒垂直度的情況，本次實(shí)驗(yàn)每間隔2 m 做一次切片，依次獲取各個截面的點(diǎn)云數(shù)據(jù)，并導(dǎo)出各截面上各個點(diǎn)坐標(biāo)。圖3 為風(fēng)機(jī)塔筒切片結(jié)果。

圖3 風(fēng)機(jī)塔筒切片結(jié)果

通過1 mm 平面切片每間隔2 m 截取一個平面，實(shí)驗(yàn)共截取了16 個平面數(shù)據(jù)，獲得的點(diǎn)云數(shù)據(jù)個數(shù)范圍為28 ～464。如表1 所示，該風(fēng)力發(fā)電機(jī)塔筒外形為椎體，所以由底部到頂部點(diǎn)云個數(shù)呈遞減趨勢。

表1 各截面點(diǎn)云個數(shù)統(tǒng)計

3 垂直度計算原理

3.1 切片圓心

詹慶明[6]指出，最小二乘法比較適用于圓形物體的特征提取。設(shè)圓周圍上有個點(diǎn)，任意兩個點(diǎn)坐標(biāo)為(xi,yi)，(xj,yj)，0

通過變換可得到二元一次方程：

可列出n-1 個方程，由此組成一個方程組，該方程的誤差方程為：

式中：A為序數(shù)矩陣，l為真實(shí)值與測量值之差。按照最小二乘原理，各個觀測點(diǎn)到擬合曲線距離的平方和最小，即值最小。將式（3）代入其中，并對x求導(dǎo)使其等于0，得到，即ATA=0，得到：

3.2 垂直度

截取的部分立柱分段圖如圖4 所示。點(diǎn)A為頂部截面的中心點(diǎn)即擬合圓心點(diǎn)，點(diǎn)B為底面截面的中心點(diǎn)即底面擬合圓心點(diǎn)，點(diǎn)C為點(diǎn)A在底面的垂直投影。當(dāng)點(diǎn)B與點(diǎn)C不重合，說明這個截面段是傾斜的，此時就會出現(xiàn)水平偏移量。點(diǎn)B與點(diǎn)C的橫坐標(biāo)之差稱為?X，點(diǎn)B與點(diǎn)C的縱坐標(biāo)之差稱為?Y。計算偏移量為

圖4 垂直度計算示意圖

對應(yīng)的傾斜度為

式中：h為頂部截面與底部截面的高差值。

4 數(shù)據(jù)分析

4.1 數(shù)據(jù)分析

在垂直度計算中，傾斜度和偏移量是衡量構(gòu)筑物安全性的重要指標(biāo)。按照上述原理，利用Matlab軟件編程擬合各截面圓心，并以2 m 為底，計算其與其他高度截面圓心坐標(biāo)的偏移量和傾斜度。各截面擬合圓心及相關(guān)參數(shù)如表2 所示。

表2 各截面圓心坐標(biāo)和偏移量

為更直觀地觀察風(fēng)力發(fā)電機(jī)塔架垂直度情況，以截面高度為橫軸，以偏移量和傾斜度為縱軸，分別繪制柱狀圖，如圖5 所示。

圖5 不同截面高度的偏移量和傾斜度

從表2 和圖5（a）可看出，隨著截面高度上升，風(fēng)力發(fā)電機(jī)塔筒偏移量整體呈上升趨勢。當(dāng)截面高度到達(dá)34 m 時，偏移量到達(dá)最大值6.29 mm。但截面高度為20 m 和22 m 時，偏移量反而減小，未隨截面高度的增加而增加。從圖5（b）可看出，塔筒傾斜度與偏移量變化基本一致，整體趨勢為，隨著截面高度上升，傾斜度逐漸增大，其值在2.18′處上下浮動。同樣地，當(dāng)截面高度為34 m 時，傾斜度到達(dá)最大值2.83′。當(dāng)截面高度為8 m 和22 m 時，塔筒傾斜度結(jié)果與其他相差較大。對此，可能有以下幾點(diǎn)原因：

（1）風(fēng)力發(fā)電機(jī)塔筒本身有變形，在數(shù)據(jù)處理時未修復(fù)；

（2）掃描過程中，隨著高度的增加，掃描精度變?nèi)酰沂墉h(huán)境限制，掃描結(jié)果會受影響；

（3）拼接過程中，采用特征點(diǎn)手動拼接，拼接成果可能存在偏差；

（4）去噪過程中，在某些高度上會有其他噪聲點(diǎn)未被刪除，以至于影響最終結(jié)果；

（5）切片過程中，由于切平面是有一定厚度的，所謂的“切片”，實(shí)際上是“切圓柱”，一定程度上也會影響圓心坐標(biāo)擬合。

4.2 結(jié)果討論

針對以上原因，在之后的研究中，可從以下幾方面進(jìn)行完善。

（1）設(shè)計數(shù)據(jù)采集方案時，可考慮利用標(biāo)靶球作為標(biāo)志點(diǎn)，提高拼接成果精度。本次實(shí)驗(yàn)采用RiScan PRO 軟件進(jìn)行拼接，軟件中的算法在某些方面還不夠精密，可根據(jù)情況設(shè)計拼接算法。

（2）可分別針對風(fēng)機(jī)塔筒形變、去噪分別設(shè)計算法，修復(fù)風(fēng)機(jī)塔筒，去除噪聲。

（3）嘗試采用穩(wěn)健估計等多種算法進(jìn)行圓心坐標(biāo)的擬合。

5 結(jié) 語

本文針對風(fēng)力發(fā)電機(jī)塔筒，基于三維激光掃描技術(shù)對風(fēng)力發(fā)電機(jī)塔筒進(jìn)行建模，采用最小二乘法擬合各截面圓心坐標(biāo)，計算其偏移量、傾斜度，實(shí)現(xiàn)對垂直度的檢測。本文驗(yàn)證了三維激光掃描技術(shù)方法在垂直度檢測工作中具有可行性，滿足工程需要，對三維激光技術(shù)在其他領(lǐng)域的應(yīng)用有一定參考意義。