基于慣性測(cè)量的輸電導(dǎo)線覆冰運(yùn)動(dòng)特征提取

楊金顯， 晁麗君， 李雙磊

(1. 河南理工大學(xué) 電氣工程與自動(dòng)化學(xué)院，河南 焦作 454003；2.國(guó)網(wǎng)山東臨朐供電公司，山東 臨朐 262600)

由于不同氣候條件影響，輸電導(dǎo)線會(huì)出現(xiàn)不同的運(yùn)動(dòng)狀態(tài)，從振幅和頻率看，主要為高頻微幅的微風(fēng)振動(dòng)，中頻中幅的次檔距振蕩和低頻大振幅的覆冰運(yùn)動(dòng)。其中，覆冰運(yùn)動(dòng)是對(duì)輸電線路危害最大的一種運(yùn)動(dòng)。覆冰運(yùn)動(dòng)公知的幅值范圍約為導(dǎo)線直徑的5倍～300倍、頻率范圍約為0.1～3 Hz?？偟膩?lái)說(shuō)輸電導(dǎo)線覆冰運(yùn)動(dòng)包括平動(dòng)(水平和垂直運(yùn)動(dòng)統(tǒng)稱為平動(dòng))和扭動(dòng)等，覆冰初期導(dǎo)線在平衡位置附近做微幅振動(dòng)，隨著覆冰程度增大和風(fēng)激勵(lì)的影響，水平振動(dòng)向垂直振動(dòng)轉(zhuǎn)變，并逐漸趨于穩(wěn)定。當(dāng)幅值和扭動(dòng)達(dá)到一定程度時(shí)，形成以垂直方向運(yùn)動(dòng)為主的低頻大振幅舞動(dòng)，舞動(dòng)使導(dǎo)線弧垂度增大、軸向應(yīng)力增加，疲勞極限降低，極易造成導(dǎo)線斷股；同時(shí)會(huì)使絕緣子串、橫擔(dān)及桿塔荷載變大，造成絕緣子擺動(dòng)、橫擔(dān)變形、桿塔傾斜等嚴(yán)重災(zāi)害，對(duì)電網(wǎng)安全運(yùn)行構(gòu)成極大威脅，對(duì)國(guó)民經(jīng)濟(jì)造成直接損失[1]。因此導(dǎo)線運(yùn)動(dòng)阻尼器的設(shè)計(jì)和安裝就變得尤為重要，但對(duì)于某個(gè)具體地區(qū)來(lái)說(shuō)，導(dǎo)線運(yùn)動(dòng)存在主要的幅值和頻率，要想針對(duì)性設(shè)計(jì)出適合某地區(qū)的運(yùn)動(dòng)阻尼器就需要識(shí)別出某地區(qū)導(dǎo)線的具體運(yùn)動(dòng)特征信息，進(jìn)而根據(jù)運(yùn)動(dòng)特征信息設(shè)計(jì)阻尼器參數(shù)。為有效預(yù)警和阻尼某地區(qū)覆冰運(yùn)動(dòng)，采集輸電導(dǎo)線三軸加速度(平動(dòng))和三軸角速度(扭動(dòng))運(yùn)動(dòng)信息，經(jīng)過(guò)預(yù)處理后計(jì)算多維時(shí)間序列的相關(guān)系數(shù)矩陣，對(duì)相關(guān)系數(shù)矩陣進(jìn)行奇異值分解提取運(yùn)動(dòng)特征，然后通過(guò)相關(guān)分析法估計(jì)導(dǎo)線運(yùn)動(dòng)的幅值和頻率，從而為阻尼器的設(shè)計(jì)提供理論參數(shù)。

近年來(lái)，由于MEMS傳感器體積小、功耗低、成本低、抗振性能好和可靠性高等優(yōu)點(diǎn)，已有學(xué)者采用MEMS傳感器對(duì)輸電線路運(yùn)動(dòng)進(jìn)行監(jiān)測(cè)研究[2]。Wang等[3-4]采用MEMS加速度計(jì)多點(diǎn)監(jiān)測(cè)導(dǎo)線振動(dòng)加速度，根據(jù)異常加速度信息檢測(cè)舞動(dòng)并進(jìn)行等級(jí)劃分；邵穎彪等[5]提出采用九軸傳感器測(cè)量導(dǎo)線舞動(dòng)的扭轉(zhuǎn)和平動(dòng)，用來(lái)克服只用加速度計(jì)監(jiān)測(cè)導(dǎo)線舞動(dòng)只能得到位移信息的缺陷；Huang等[6-7]考慮導(dǎo)線自激振蕩和導(dǎo)線扭轉(zhuǎn)引起的導(dǎo)線舞動(dòng)，設(shè)計(jì)了一種舞動(dòng)慣性測(cè)量單元，避免導(dǎo)線舞動(dòng)時(shí)計(jì)算的相對(duì)位移與導(dǎo)線扭轉(zhuǎn)引起的實(shí)際位移的偏差；汪滔等[8]提出一種改進(jìn)的MEMS陀螺儀和加速度計(jì)互補(bǔ)濾波算法，設(shè)計(jì)了舞動(dòng)軌跡的解算方法，對(duì)陀螺計(jì)算出的角度進(jìn)行誤差修正。在上述裝置和方法中：加速度計(jì)會(huì)隨導(dǎo)線發(fā)生扭轉(zhuǎn)，對(duì)通過(guò)積分計(jì)算得到的輸電線運(yùn)動(dòng)軌跡等參數(shù)會(huì)造成計(jì)算不準(zhǔn)，盡管慣性傳感器能夠克服單獨(dú)使用加速度傳感器帶來(lái)的那些問(wèn)題，但實(shí)現(xiàn)這種功能需要導(dǎo)航級(jí)的慣性測(cè)量系統(tǒng)，體積大且成本昂貴，況且姿態(tài)解算需要精確的初始值，而且產(chǎn)生的累計(jì)誤差；另外采用單一特征對(duì)導(dǎo)線運(yùn)動(dòng)分析，會(huì)割裂數(shù)據(jù)整體的內(nèi)在聯(lián)系，無(wú)法準(zhǔn)確識(shí)別輸電導(dǎo)線的整體運(yùn)動(dòng)特征。

在以上研究的基礎(chǔ)上，提出一種MIMS多維時(shí)間序列覆冰運(yùn)動(dòng)特征提取方法：通過(guò)MEMS慣性測(cè)量系統(tǒng)測(cè)量導(dǎo)線運(yùn)動(dòng)加速度和角速度等信息，不通過(guò)姿態(tài)解算，而是直接利用原始信息進(jìn)行處理，對(duì)多維時(shí)間序列的相關(guān)矩陣進(jìn)行奇異值分解，構(gòu)建加權(quán)特征矩陣作為輸電導(dǎo)線整體運(yùn)動(dòng)特征，然后采用近鄰分類法進(jìn)行覆冰舞動(dòng)識(shí)別。同樣，該算法也適用于其它導(dǎo)線運(yùn)動(dòng)形式檢測(cè)，當(dāng)待識(shí)別的加權(quán)特征矩陣與輸電導(dǎo)線其他運(yùn)動(dòng)形式時(shí)的加權(quán)特征矩陣相似度較高時(shí)，可判定導(dǎo)線以該形式運(yùn)動(dòng)。本文所提方法降低了數(shù)據(jù)的維數(shù)，提高了處理速度和運(yùn)動(dòng)特征識(shí)別精度。

1 MIMS數(shù)據(jù)采集與處理

1.1 信號(hào)采集與預(yù)處理

由于導(dǎo)線附近具有較強(qiáng)的電磁干擾，MEMS慣性測(cè)量系統(tǒng)不能直接安裝于導(dǎo)線上，為此采用一定長(zhǎng)度的防磁中空彈性桿，一端固定于絕緣子串懸掛點(diǎn)，另一端安裝MEMS慣性測(cè)量節(jié)點(diǎn)，來(lái)測(cè)量導(dǎo)線運(yùn)動(dòng)的三軸加速度和三軸角速度等信息。彈性桿內(nèi)放置電源線、地線和數(shù)據(jù)線，并且數(shù)據(jù)線封裝于屏蔽層內(nèi)，采用電流互感器整流變壓后直接給測(cè)量節(jié)點(diǎn)供電，為避免信號(hào)受到干擾，把采集的數(shù)據(jù)通過(guò)電力光纜傳輸給監(jiān)控中心進(jìn)行數(shù)據(jù)處理[9]。

目前對(duì)于輸電導(dǎo)線運(yùn)動(dòng)的測(cè)量主要是對(duì)單一特征的識(shí)別，但是由于MIMS各軸加速度和角速度之間存在著內(nèi)在耦合關(guān)系，單一軸向的加速度或角速度無(wú)法表征整個(gè)輸電導(dǎo)線的運(yùn)動(dòng)特征，而多維時(shí)間序列既可以描述某一時(shí)刻導(dǎo)線三個(gè)軸向的轉(zhuǎn)動(dòng)運(yùn)動(dòng)信息，又包含了不同時(shí)刻導(dǎo)線平動(dòng)以及這些運(yùn)動(dòng)的相互關(guān)系。為了準(zhǔn)確識(shí)別覆冰運(yùn)動(dòng)，考慮外部因素影響及導(dǎo)線運(yùn)動(dòng)的整體關(guān)系，將三軸加速度和三軸角速度信號(hào)由6維時(shí)間序列表示：

X={ax,ay,az,ωx,ωy,ωz}

(1)

式中:ax,ay,az為三軸加速度序列，ωx,ωy,ωz為三軸角速度序列，將上式記為X={x1,x2,…,xi,…,x6}，xi={ai1,ai2,…,aiN}T，N為序列長(zhǎng)度。

MEMS慣性傳感器在工作過(guò)程中有周期漂移項(xiàng)的存在，會(huì)對(duì)運(yùn)動(dòng)分析結(jié)果造成影響；另外隨機(jī)干擾信號(hào)也是系統(tǒng)主要誤差源之一，傳感器在工作中極易產(chǎn)生低頻周期性信號(hào)[10]。由于這些周期趨勢(shì)項(xiàng)的存在，通過(guò)二次積分所得的運(yùn)動(dòng)位移軌跡可能有較大程度的偏離實(shí)際，甚至產(chǎn)生畸變或失真，因此需要消除趨勢(shì)項(xiàng)和隨機(jī)噪聲，需采取合適的方法減少信號(hào)中的干擾成分。常用的去噪方法有小波去噪和Kalman去噪，雖然這兩種方法都能在一定程度上平滑MEMS信號(hào)，但采用小波去噪方法閾值選取比較困難，且依賴于信號(hào)噪聲統(tǒng)計(jì)特性；基于Kalman濾波的去噪方法，只適用于噪聲特性易知的場(chǎng)合，然而由于外界環(huán)境影響，MEMS信號(hào)噪聲較為復(fù)雜，且具有時(shí)變性，動(dòng)態(tài)系統(tǒng)模型的狀態(tài)噪聲和觀測(cè)噪聲統(tǒng)計(jì)特性往往具有不確定性。因此上述依賴噪聲統(tǒng)計(jì)特性的方法，在輸電導(dǎo)線運(yùn)動(dòng)信號(hào)處理中很難獲得理想的去噪效果。由于五點(diǎn)三次平滑法可減小時(shí)域信號(hào)中的隨機(jī)噪聲，有效處理周期性的干擾信號(hào)或不規(guī)則的隨機(jī)干擾信號(hào)所引起的尖峰、毛刺現(xiàn)象，適用于輸電導(dǎo)線運(yùn)動(dòng)的MEMS信號(hào)去噪場(chǎng)合，且計(jì)算簡(jiǎn)單方便，因此本文采用五點(diǎn)三次平滑法對(duì)數(shù)據(jù)進(jìn)行濾波處理，其計(jì)算公式為：

(2)

1.2 信號(hào)處理效果評(píng)價(jià)

為了描述輸電導(dǎo)線在不同時(shí)刻的運(yùn)動(dòng)加速度和角速度的變化關(guān)系及其變化趨勢(shì)，建立一多維時(shí)間序列X的矩陣模型，以便于后續(xù)運(yùn)動(dòng)特征提取。為消除偶然測(cè)量誤差，采用移動(dòng)滑窗的方法對(duì)數(shù)據(jù)進(jìn)行分段，將長(zhǎng)度為L(zhǎng)的時(shí)間窗在多維時(shí)間序列上滑動(dòng)并取平均。對(duì)處理過(guò)的多維時(shí)間序列建立時(shí)變測(cè)量矩陣Xn×6，其中，n=N/L。

(3)

信號(hào)預(yù)處理后，通過(guò)信噪比來(lái)驗(yàn)證信號(hào)處理效果。信噪比是測(cè)量信號(hào)中的噪聲量的傳統(tǒng)方法，常被用作去噪效果評(píng)價(jià)指標(biāo)。信噪比定義如下，單位為dB：

SNR=10ln(Ps/Pn)

(4)

NRR=1/(SNR+1)+ε

(5)

式中：NRR為噪聲剩余率，ε為誤差因子。當(dāng)NRR∈(0,0.1)時(shí)，說(shuō)明剩余噪聲量較少，去噪效果較好；當(dāng)NRR∈(0.1,0.3)時(shí)，說(shuō)明剩余噪聲量較多，去噪效果一般；當(dāng)NRR∈(0.3,1)時(shí)說(shuō)明剩余大量噪聲，去噪效果不好。

由于加速度和角速度量綱和運(yùn)動(dòng)特征方面具有差異，根據(jù)max-min法將導(dǎo)線運(yùn)動(dòng)三軸加速度和三軸角速度的值歸一化得到新的各維的序列值：

(6)

式中：xi max和xi min分別為各維加速度或角速度測(cè)量序列的最大值和最小值。

2 運(yùn)動(dòng)特征提取與識(shí)別

2.1 信號(hào)特征提取

盡管輸電覆冰線路運(yùn)動(dòng)形式復(fù)雜，幅度和頻率都在變化，但其基本運(yùn)動(dòng)模式為平動(dòng)和扭動(dòng)。導(dǎo)線平動(dòng)時(shí)，X、Y和Z軸加速度、角速度信號(hào)與各自軸向相似程度較高，自相關(guān)系數(shù)較大；當(dāng)導(dǎo)線的扭轉(zhuǎn)頻率與其橫向運(yùn)動(dòng)頻率相接近時(shí)，覆冰導(dǎo)線產(chǎn)生偏心慣性，誘使導(dǎo)線發(fā)生扭動(dòng)，X軸加速度和角速度向Y軸或Z軸轉(zhuǎn)化，三軸發(fā)生交叉耦合運(yùn)動(dòng)，X、Y和Z軸信號(hào)與其他軸向相似程度較高，互相關(guān)系數(shù)較大。

導(dǎo)線運(yùn)動(dòng)多維時(shí)間序列觀測(cè)值數(shù)目不同時(shí)，無(wú)法采用傳統(tǒng)的歐氏距離進(jìn)行相似度分析，為了體現(xiàn)MIMS各軸加速度和角速度之間的相互關(guān)系，同時(shí)解決MEMS測(cè)量數(shù)目不同的問(wèn)題，計(jì)算運(yùn)動(dòng)矩陣X′的相關(guān)系數(shù)矩陣RX′[11]，對(duì)信號(hào)整體內(nèi)在關(guān)系進(jìn)行分析。

(7)

采用相關(guān)系數(shù)矩陣對(duì)輸電導(dǎo)線運(yùn)動(dòng)進(jìn)行分析，既可以表征導(dǎo)線運(yùn)動(dòng)的加速度和角速度相互關(guān)系，又將運(yùn)動(dòng)特征矩陣固定為6×6維，便于矩陣之間相似性度量。由于噪聲干擾的存在，會(huì)對(duì)輸電導(dǎo)線運(yùn)動(dòng)的識(shí)別結(jié)果造成影響，因此對(duì)導(dǎo)線運(yùn)動(dòng)的相關(guān)系數(shù)矩陣RX′進(jìn)行奇異值分解提取導(dǎo)線的運(yùn)動(dòng)特征[12]。

RX′=UΛUT

(8)

2.2 加權(quán)矩陣計(jì)算

奇異值可以反映出MIMS采集的信號(hào)中的有用信號(hào)和噪聲的能量分布特點(diǎn)，在奇異值分解中，大的奇異值對(duì)應(yīng)的特征向量包含了更多導(dǎo)線運(yùn)動(dòng)信息，具有較高的鑒別性，特征矩陣中能夠表征MIMS信號(hào)特征的有用信號(hào)主要由前k個(gè)較大的奇異值反映出來(lái)，而后6-k個(gè)較小的奇異值對(duì)輸電導(dǎo)線運(yùn)動(dòng)影響不大。較大的奇異值對(duì)應(yīng)的奇異向量在矩陣特征識(shí)別中所占權(quán)重較大，而奇異值所占比例大小可以反映出特征向量所占權(quán)重。通過(guò)對(duì)特征矩陣進(jìn)行加權(quán)，進(jìn)行輸電導(dǎo)線的運(yùn)動(dòng)特征識(shí)別[13]。采用歸一化奇異值的方法計(jì)算加權(quán)矩陣：

(9)

(10)

利用對(duì)多維時(shí)間序列的相關(guān)系數(shù)矩陣奇異值分解后的特征值構(gòu)造加權(quán)特征矩陣，通過(guò)計(jì)算兩個(gè)多維時(shí)間序列加權(quán)特征矩陣之間的相似度，識(shí)別導(dǎo)線運(yùn)動(dòng)形式，大大降低了計(jì)算量，并提高了運(yùn)動(dòng)形式識(shí)別精度。

2.3 覆冰運(yùn)動(dòng)識(shí)別

(11)

對(duì)于輸電導(dǎo)線的覆冰運(yùn)動(dòng)特征識(shí)別，首先對(duì)MIMS采集的原始多維時(shí)間序列進(jìn)行預(yù)處理并建立時(shí)變測(cè)量矩陣；然后計(jì)算相關(guān)系數(shù)矩陣并進(jìn)行奇異值分解，根據(jù)特征向量和特征值計(jì)算加權(quán)特征矩陣，對(duì)導(dǎo)線運(yùn)動(dòng)特征進(jìn)行提取；最后通過(guò)度量當(dāng)前序列的加權(quán)特征矩陣和基準(zhǔn)特征矩陣之間的相似性，采用近鄰分類法識(shí)別覆冰運(yùn)動(dòng)，具體算法流程如圖1所示。

圖1 覆冰運(yùn)動(dòng)識(shí)別算法流程Fig.1 Motion identification algorithm procedure of iced conductor

2.4 覆冰運(yùn)動(dòng)幅值和頻率估計(jì)

由于輸電網(wǎng)絡(luò)覆冰運(yùn)動(dòng)的復(fù)雜性，擬合覆冰運(yùn)動(dòng)軌跡比較困難且容易出現(xiàn)解算誤差。輸電導(dǎo)線覆冰運(yùn)動(dòng)時(shí)，平動(dòng)與扭動(dòng)伴隨而生，通過(guò)對(duì)加速度和角速度信號(hào)進(jìn)行相關(guān)分析，進(jìn)而估計(jì)覆冰導(dǎo)線幅值和頻率變化程度。首先對(duì)MEMS信號(hào)平滑分段并歸一化處理，將加速度和角速度特征量尺度變換到相同的范圍，以便于統(tǒng)一計(jì)算。將輸電線路靜止或微風(fēng)振動(dòng)情況下的時(shí)變測(cè)量矩陣作為初始矩陣，初始矩陣第i時(shí)間段的信號(hào)記為xi，通過(guò)計(jì)算覆冰測(cè)量矩陣與初始矩陣在同一時(shí)間段的相關(guān)系數(shù)，度量MEMS信號(hào)的變異程度，進(jìn)而根據(jù)相似程度預(yù)測(cè)覆冰導(dǎo)線幅值變化程度(相關(guān)系數(shù)越小，數(shù)據(jù)相似程度越低，幅值變化程度越大)。覆冰測(cè)量矩陣與初始矩陣第i時(shí)間段的相關(guān)系數(shù)Ra為：

(12)

覆冰導(dǎo)線運(yùn)動(dòng)信號(hào)的變化頻率即導(dǎo)線運(yùn)動(dòng)頻率。通過(guò)計(jì)算覆冰測(cè)量矩陣當(dāng)前時(shí)間段與上一時(shí)間段的相關(guān)系數(shù)，度量信號(hào)的變異速率，進(jìn)而根據(jù)變異速率預(yù)測(cè)覆冰導(dǎo)線頻率變化程度(相關(guān)系數(shù)越小，數(shù)據(jù)相似程度越低，頻率變化程度越快)。相鄰時(shí)間段信號(hào)的相關(guān)系數(shù)Rc為：

(13)

分別記d為幅值變異程度，f為頻率變異速率，由于數(shù)據(jù)存在一定誤差，考慮誤差因素，則[15]：

d=(1-Ra)+ε1

(14)

f=(1-Rc)+ε2

(15)

式中：ε1、ε2為誤差因子。

輸電導(dǎo)線微風(fēng)振動(dòng)時(shí)為高頻微幅運(yùn)動(dòng)，次檔距振蕩時(shí)為中頻中幅運(yùn)動(dòng)，覆冰舞動(dòng)時(shí)為低頻大振幅運(yùn)動(dòng)。為進(jìn)一步驗(yàn)證覆冰識(shí)別算法的有效性，按照上述算法估計(jì)導(dǎo)線在微風(fēng)振動(dòng)和次檔距振蕩運(yùn)動(dòng)狀態(tài)下幅值和頻率變異程度，并和覆冰運(yùn)動(dòng)進(jìn)行比較。

3 實(shí) 驗(yàn)

為驗(yàn)證上述輸電線路覆冰運(yùn)動(dòng)識(shí)別及幅頻估計(jì)算法，按一定比例設(shè)計(jì)了輸電線路模型，將MIMS按1.1節(jié)方案安裝于模型，搭建如圖2所示輸電線路運(yùn)動(dòng)識(shí)別系統(tǒng)，采集輸電導(dǎo)線在微風(fēng)振動(dòng)和次檔距振蕩時(shí)的數(shù)據(jù)作為基準(zhǔn)數(shù)據(jù)并進(jìn)行特征提取，通過(guò)監(jiān)控主機(jī)或智能終端監(jiān)測(cè)導(dǎo)線運(yùn)動(dòng)并保存數(shù)據(jù)。

圖2 輸電網(wǎng)絡(luò)運(yùn)動(dòng)測(cè)試系統(tǒng)Fig.2 Motion test system of power transmission network

由于白蠟密度與冰相似，將導(dǎo)線拉直偏置于橢圓模具中，將白蠟融化注入模具覆在導(dǎo)線上，在白蠟上涂上一層膠，在硬度上也和冰保持一致，模擬偏心覆冰導(dǎo)線，采集MIMS信號(hào)作為覆冰基準(zhǔn)數(shù)據(jù)并進(jìn)行特征提取，模擬覆冰導(dǎo)線,如圖3所示。

采集任意一組輸電導(dǎo)線覆冰運(yùn)動(dòng)時(shí)的MEMS信號(hào)通過(guò)五點(diǎn)三次平滑法預(yù)處理，如圖4所示。由圖4可知，預(yù)處理后可有效減少隨機(jī)噪聲、尖峰和毛刺現(xiàn)象。

取L=30，將預(yù)處理數(shù)據(jù)滑動(dòng)平均分段。為定量分析數(shù)據(jù)處理效果，根據(jù)式(4)和式(5)分別計(jì)算三軸加速度和三軸角速度信號(hào)的信噪比和噪聲剩余率，噪聲剩余率如表1所示。由表1可知，三軸加速度和三軸角速度的噪聲剩余率都在(0,0.1)之間，說(shuō)明去噪效果良好。

圖3 模擬覆冰導(dǎo)線Fig.3 Analog iced conductor

表1 噪聲剩余率Tab.1 Noise residual ratio

表2 運(yùn)動(dòng)特征矩陣與三種基準(zhǔn)矩陣距離Tab.2 Distance between motion characteristic matrix and three reference matrices

采用相關(guān)分析法，對(duì)同一時(shí)間段覆冰測(cè)量矩陣與初始測(cè)量矩陣MEMS信號(hào)進(jìn)行相關(guān)分析，估計(jì)覆冰導(dǎo)線幅值變異程度；對(duì)相鄰時(shí)間段覆冰測(cè)量矩陣MEMS信號(hào)進(jìn)行相關(guān)分析，估計(jì)頻率變異速率。同時(shí)，估計(jì)導(dǎo)線微風(fēng)振動(dòng)和次檔距振蕩時(shí)的幅值和頻率，并與覆冰運(yùn)動(dòng)進(jìn)行比較，三種形式運(yùn)動(dòng)特征的幅值和頻率變化分別如圖6和圖7所示。

由圖6和圖7可知，微風(fēng)振動(dòng)時(shí)導(dǎo)線近似做高頻微幅運(yùn)動(dòng)，次檔距振蕩時(shí)導(dǎo)線近似做中頻中幅運(yùn)動(dòng)，覆冰運(yùn)動(dòng)時(shí)導(dǎo)線近似做低頻大振幅運(yùn)動(dòng)。不同運(yùn)動(dòng)狀態(tài)下導(dǎo)線幅值頻率曲線走勢(shì)進(jìn)一步說(shuō)明了算法的有效性。根據(jù)所估計(jì)的幅值頻率變化，可以獲取本地區(qū)覆冰輸電導(dǎo)線運(yùn)動(dòng)特性，為后續(xù)振動(dòng)阻尼器設(shè)計(jì)提供參數(shù)。

圖4 原始數(shù)據(jù)及預(yù)處理結(jié)果Fig.4 Raw data and preprocessing results

圖5 歸一化后三軸加速度和三軸角速度Fig.5 Normalized triaxial acceleration and angular velocity

圖6 三種運(yùn)動(dòng)形式幅值變異程度對(duì)比圖

Fig.6 Amplitude variation degree comparison diagram of three motion modes

4 總 結(jié)

針對(duì)現(xiàn)有輸電導(dǎo)線運(yùn)動(dòng)檢測(cè)方法的缺點(diǎn)，從新的角度對(duì)導(dǎo)線覆冰運(yùn)動(dòng)進(jìn)行檢測(cè)：采用MEMS慣性傳感器，既兼顧了導(dǎo)線自激振蕩和扭轉(zhuǎn)引起的運(yùn)動(dòng)對(duì)覆冰測(cè)量的干擾，也解決了單獨(dú)使用加速度傳感器產(chǎn)生扭轉(zhuǎn)誤差的問(wèn)題；采用多維時(shí)間序列解決了一維時(shí)間序列或合加速度序列進(jìn)行識(shí)別割裂了輸電導(dǎo)線整體運(yùn)動(dòng)信息；采用對(duì)相關(guān)系數(shù)矩陣進(jìn)行加權(quán)奇異值分解的方法，降低了維數(shù)，提高了覆冰運(yùn)動(dòng)識(shí)別的準(zhǔn)確性。實(shí)驗(yàn)結(jié)果表明，使用本文所提算法能夠識(shí)別導(dǎo)線運(yùn)動(dòng)情況并估計(jì)幅頻變化程度，且實(shí)施方案簡(jiǎn)單易行，無(wú)累計(jì)誤差，可為長(zhǎng)期測(cè)量某地區(qū)輸電線路運(yùn)動(dòng)狀態(tài)提供參考，進(jìn)一步為輸電線路阻尼器的設(shè)計(jì)提供參數(shù)。

圖7 三種運(yùn)動(dòng)形式頻率變異速率對(duì)比圖

Fig.7 Frequency variation degree comparison diagram of three motion modes