結(jié)構(gòu)健康監(jiān)測系統(tǒng)光纖Bragg 光柵應變計的異常診斷

孫 曉，王清梅，李振偉，楚敬敬

（1. 青島科技大學自動化與電子工程學院，山東青島266061；2. 中科院國家天文臺，北京 100012）

1 引 言

為確保結(jié)構(gòu)的安全可靠，結(jié)構(gòu)健康監(jiān)測技術(shù)［1］廣泛應用于土木工程結(jié)構(gòu)、鋼結(jié)構(gòu)等大型工程中，通過布設(shè)各類傳感器，監(jiān)測并評估結(jié)構(gòu)的健康狀態(tài)。光纖布拉格光柵（Fiber Bragg Grat?ing，F(xiàn)BG）應變計是一種光纖傳感器［2］，它利用光纖的光敏性在纖芯上建立折射率調(diào)制，光柵的周期和折射率隨安裝點的溫度和應變變化，引起反射光波長的變化，通過測量反射光波長的變化，來測量目標的應變信息。FBG 因其光信號穩(wěn)定，單根光纖可串聯(lián)多個光柵測點［3］，不受電磁干擾，適合分布傳感等特點［4］，大量應用在工程結(jié)構(gòu)的應變監(jiān)測中。

500 米口徑球形射電望遠鏡（Five-hundredmeter Aperture Spherical Telescope，F(xiàn)AST）為具有主動反射面的球冠狀大型射電望遠鏡［5］。主動反射面系統(tǒng)結(jié)構(gòu)特殊且空間跨度大，為保證主動反射面可靠工作，建設(shè)了一套主動反射面健康監(jiān)測系統(tǒng)［6］，用于反射面的故障診斷與預警。其中，用于索力監(jiān)測及圈梁鋼結(jié)構(gòu)監(jiān)測的FBG 應變計共計416 個。這些傳感器長期工作于復雜的室外環(huán)境中，易發(fā)生安裝失效、蠕變、漂移、老化和線路故障等問題，使其監(jiān)測數(shù)據(jù)錯誤或失效。傳感器故障會直接影響健康監(jiān)測系統(tǒng)的性能，必須進行及時的檢查與維護，識別并排除異常，保證數(shù)據(jù)的可靠性。如果采用人工巡檢排查的方法，不僅工作量大，而且許多傳感器的安裝位置不具備人工觀察的條件，無法及時有效地發(fā)現(xiàn)與維護故障測點，從而形成安全隱患。

國內(nèi)外學者針對FBG 傳感器的異常診斷展開了研究。Cazzulani［7］等基于測量估計殘差的分析，識別FBG 傳感器鏈中不同的故障類型，排除故障傳感器數(shù)據(jù)。張曉麗［8］等借助光波導數(shù)值分析軟件對FBG 傳感器性能退化的信號特征進行了分析，為FBG 傳感器的自診斷提供依據(jù)。姜婷睿［9］等利用小波包和主元分析結(jié)合提取船用FBG 應變計故障特征，實現(xiàn)了故障的基本診斷。傳感器異常自診斷可有效提高傳感器的工作效率與可靠性，不同傳感器的異常診斷方法存在共通性［10］，廣泛應用于機電［11］、驅(qū)動控制［12］、航空航天發(fā)動機［13］以及數(shù)據(jù)采集等系統(tǒng)中。

傳感器的異常診斷方法主要有基于模型、基于知識和基于信號處理方法等［14］。基于模型的方法需獲得系統(tǒng)的精確解析模型，且系統(tǒng)的非線性及不確定性會對診斷結(jié)果產(chǎn)生很大影響；基于知識的方法需有足夠的故障原因與特征的先驗知識積累轉(zhuǎn)化為推理規(guī)則；基于信號處理的方法只需對傳感器的信號特性進行分析，不需建立精確的模型或足夠的先驗知識，利用信號特性結(jié)合主成分分析、人工神經(jīng)網(wǎng)絡(luò)等算法就可完成對傳感器的故障診斷［15］。鑒于此，本文基于信號處理方法，利用結(jié)構(gòu)健康監(jiān)測中FBG 應變計的信號特征進行異常診斷。

2 原 理

2.1 光纖光柵應變計原理

FBG 應變計的原始輸出信號為光柵反射中心波長，當光柵外界應力或溫度變化時，反射中心波長隨之變化，輸出波長為：

式中：KT為光柵溫度靈敏系數(shù)，Kε為光柵應變靈敏系數(shù)，ΔT為光柵溫度變化值，Δε為光柵處應變變化值，λB0為中心波長的初始值，KεT為光柵溫度應變交叉敏感系數(shù)。可見FBG 應變計的原始輸出信號易受溫度影響，在實際應用中常采取相應措施，補償溫度變化對應變的影響［16］。

2.2 異常診斷原理

FAST 的主動反射面為索網(wǎng)結(jié)構(gòu)，主索網(wǎng)與支承鋼結(jié)構(gòu)的圈梁應力變化存在傳遞性。在主動反射面不變形期間，結(jié)構(gòu)應力發(fā)生變化的主要原因為環(huán)境荷載，F(xiàn)BG 應變計所承受的主要周期荷載為溫度荷載。所有應變計測點的工作環(huán)境及監(jiān)測對象承受的荷載變化規(guī)律大致相同。分析發(fā)現(xiàn)，正常FBG 應變計測點輸出信號，其波形幅值等存在相似性，統(tǒng)計不同傳感器的多種特征值，其數(shù)據(jù)分布呈聚合狀態(tài)。傳感器故障會導致信號特征發(fā)生變化，如表1 所示。所以，基于傳感器測點歷史數(shù)據(jù)，提取每個測點的關(guān)鍵特征值，利用其數(shù)據(jù)分布規(guī)律尋找特征值偏差較大的離群測點，可對所有FBG 傳感器進行自動化異常診斷。

表1 傳感器故障與特征變化Tab.1 Failure types and characteristic change of sensors

3 數(shù)據(jù)預處理

FAST 主動反射面健康監(jiān)測系統(tǒng)傳感器的數(shù)據(jù)特點如下：（1）時間上存在斷層，F(xiàn)AST 對電磁環(huán)境極為敏感，健康監(jiān)測系統(tǒng)在工作過程中并非持續(xù)監(jiān)測；（2）部分測點數(shù)據(jù)缺失，由于軟件設(shè)置或者硬件設(shè)備的檢修，部分傳感器數(shù)據(jù)通道被關(guān)閉，需加以篩選；（3）數(shù)據(jù)量大，采樣周期為1 s，單測點每月數(shù)據(jù)可達250 多萬條，附帶傳感器其他信息，單測點每月數(shù)據(jù)可達200 MB 以上。在進行故障特征提取時，需對數(shù)據(jù)庫內(nèi)的歷史數(shù)據(jù)進行預處理。

3.1 數(shù)據(jù)抽取

根據(jù)傳感器的編號及標識符，循環(huán)讀取健康監(jiān)測系統(tǒng)中單測點規(guī)定時間周期內(nèi)的歷史數(shù)據(jù)。由于傳感器的主要荷載為溫度荷載，其變化周期以天為單位。為豐富信號特征，抽取7 天左右的歷史數(shù)據(jù)，去除與診斷無關(guān)的信息，僅保留診斷所需的原始波長信號，作為原始樣本λ(n)，如圖1所示。

圖1 測點輸出原始信號樣本示例Fig.1 Original samples of measuring point example

3.2 壓縮濾波

根據(jù)FAST 的運行特點，結(jié)構(gòu)主要荷載的變化周期較長，為壓縮樣本數(shù)據(jù)量，便于轉(zhuǎn)換與計算，對原始樣本進行采樣壓縮，壓縮采樣周期選為1 000 s，滿足采樣定理的要求。再對數(shù)據(jù)低通濾波，濾除高次諧波的干擾，壓縮采樣公式為：

式中：N為原始樣本長度，壓縮樣本長度為N/1 000。采用三階低通數(shù)字濾波器濾波，濾波公式為：

式中：bi，aj為濾波器系數(shù)，x為輸入信號，y為輸出信號。壓縮并濾波后的數(shù)據(jù)作為壓縮樣本λ*(m)，示例如圖 2 所示。

圖2 壓縮采樣后信號示例Fig.2 Samples after compression and filtering of measur?ing point example

4 特征提取

正常FBG 應變計近似的信號統(tǒng)計指標主要有信號長度、標準差、時域能量值、頻域主成分周期、幅值、方差、波形、峰值和脈沖等。這些指標的實質(zhì)評價作用存在重復，經(jīng)分析選擇前4 項作為特征進行診斷。

4.1 原始樣本長度

FAST 健康監(jiān)測系統(tǒng)軟件采集數(shù)據(jù)時，測量數(shù)值波動超過一定閾值則會舍棄，所以部分異常測點的數(shù)據(jù)表現(xiàn)為大量缺失。軟件錯誤，偶發(fā)故障等原因，也會造成數(shù)據(jù)缺失。按時間段抽取數(shù)據(jù)，每個測點所得的樣本長度不一致，所以測點最終讀取的樣本長度是評判傳感器狀態(tài)的重要特征之一。樣本長度特征記為：

4.2 原始樣本標準差

測點信號標準差能反映一個數(shù)據(jù)集的離散程度，可作為評判測點信號穩(wěn)定程度的依據(jù)。為保留更多有效信息，使用測點原始樣本，計算公式為：

4.3 壓縮樣本時域能量值

時域能量值反映測點輸出信號的波動強度，對輸出幅值無明顯變化或者增益過大有較強的指示作用。為避免異常值干擾，使用壓縮樣本，先對數(shù)據(jù)去均值化，再計算該測點能量值，計算公式為：

式中：M=N/1 000 取正整數(shù)，為壓縮濾波后樣本數(shù)據(jù)的長度為壓縮樣本平均值。

4.4 信號頻域主成分周期

在FAST 健康監(jiān)測系統(tǒng)中，造成結(jié)構(gòu)應力變化較大的主要荷載是環(huán)境溫度。正常FBG 應變計輸出的原始波長信號的主要荷載成分受溫度影響，波動周期在一天左右。計算信號頻譜，獲取頻域信號的最大峰值，由它對應頻率換算得信號頻域主成分周期，可作為信號周期特征是否正常的評價依據(jù)。

使用壓縮樣本，先對數(shù)據(jù)去均值化，信號頻譜函數(shù)的計算公式為：

圖3 壓縮信號頻譜示例Fig.3 Frequency characteristics example of compressed signal

提取峰值對應的頻率fmax，換算頻域主成分對應的周期特征值，為使數(shù)據(jù)更直觀，換算單位為小時：

循環(huán)提取所有測點特征值，若有K個測點正常提取到數(shù)據(jù)并存在特征值，可得特征值列表為：

5 異常診斷

5.1 診斷方法

分析式（9）可知，正常傳感器的特征值在數(shù)據(jù)上較為接近，呈聚合狀態(tài)，而異?；蚬收蟼鞲衅鲃t呈離群狀態(tài)。為識別測點對應特征的離群程度，評判傳感器的異常狀態(tài)，需獲取特征值中心點。

5.1.1 獲取聚合中心點

為避免離群過大的數(shù)值對特征中心的影響，中心點不使用數(shù)據(jù)均值，而是采用循環(huán)迭代的方式搜尋附近聚合點數(shù)最多的測點特征值作為一維數(shù)據(jù)的聚合中心點。以樣本長度特征為例，定義數(shù)據(jù)收斂系數(shù)為：

式中：a為收斂系數(shù)倍數(shù)，可通過調(diào)節(jié)a來調(diào)節(jié)判斷靠近聚合點的距離，取值在0.5~1.5 為佳，此處a取0.5；K為特征值列表長度。

兩點之間的歐式距離為：

式中i，j為兩不同測點序號。

若D(i，j)

圖4 循環(huán)迭代流程Fig.4 Flow chart of loop iteration

數(shù)組num 中最大值，即附近聚合點數(shù)最多的數(shù)據(jù)點，對應腳標imax所對應特征值L（imax）作為中心點數(shù)值，記作L（center）。若有多個數(shù)據(jù)點聚合數(shù)目均為相等最大值，則取這多個數(shù)據(jù)點對應特征值的均值，作為特征值中心點。同理，可獲取其他三個特征值的中心點，得［L（center），σ（center），E（center），T（center）］。

5.1.2 傳感器綜合異常指數(shù)

離群程度評價的目的是使傳感器每類特征具備同等權(quán)重，避免因為特征值本身數(shù)量級差別造成的影響。計算特征離群程度時，引入Z-score思想對偏離距離進行標準化，每個特征向量標準差記作：

定義傳感器綜合異常指數(shù)為：

式中：i=1，2，…，K，為測點編號。

當傳感器發(fā)生卡死時，輸出恒定不變，提取主成分周期結(jié)果為空值，進一步處理時需將空值替換為 0，且將d（i）+3 作為故障標記。

式中：k1，k2為故障系數(shù)與異常系數(shù)，1

5.2 仿真驗證

取FAST 健康監(jiān)測系統(tǒng)中FBG 傳感器2019年7 月份1~9 號測點作為正常傳感器樣本，并模擬幾種傳感器的故障數(shù)據(jù)作為故障樣本。

5.2.1 傳感器數(shù)據(jù)大量缺失

對3 號測點數(shù)據(jù)進行采樣壓縮，壓縮4 倍，作為模擬10 號測點，得到：

5.2.2 傳感器恒增益

測量信號呈倍數(shù)增長，對3 號測點數(shù)據(jù)的輸出幅值進行放大，增益4 倍，作為模擬11 號測點，公式為：

5.2.3 傳感器卡死

傳感器卡死表現(xiàn)為輸出為某一固定值疊加微弱的噪聲信號。實際情況中，由于傳感器本身精度原因，以及采集系統(tǒng)在數(shù)據(jù)記錄時預處理，微弱噪聲往往已被消除，且本方法會進行去均值處理。采用與3 號測點同長度的數(shù)據(jù)，每個數(shù)據(jù)點替換為0，作為模擬12 號測點。

5.2.4 輸出信號隨機

使用生成的隨機數(shù)據(jù)與3 號測點幅值的乘積作為模擬數(shù)據(jù)，作為模擬13 號測點。

故障測點4 個，正常測點9 個，為增加正常測點比例，將測點1~9 號數(shù)據(jù)重復一次，作為14~22 號測點。使用診斷方法對仿真模擬數(shù)據(jù)進行診斷，診斷結(jié)果如圖5 所示，特征值與異常指數(shù)見表2。測試證明，當故障測點占總測點數(shù)量的20%以內(nèi)時，可有效識別故障傳感器。異常測點過多，會導致特征值分布過于離散，診斷效果不佳，此時需調(diào)整故障和異常系數(shù)，進一步嚴格指數(shù)偏離程度。仿真結(jié)果顯示，當故障測點比例每增加 10%，故障系數(shù)k1大約減少（k1-1）/2，如故障測點占比30%時k1調(diào)整為2，40%時k1調(diào)整為1.5，異常系數(shù)仍選擇中值；當故障測點比例占50%時，故障系數(shù)低至1.1，此時傳感器故障過半，正常傳感器與異常傳感器界定不明顯，方法無法診斷。

圖5 仿真實驗結(jié)果Fig.5 Simulation experiment results

表2 仿真?zhèn)鞲衅鞯奶卣骷爱惓Ｖ笖?shù)Tab.2 Eigen values and abnormal indexes of simulated sensors

5.3 實際應用

以FAST 健康監(jiān)測系統(tǒng)2019 年7 月份的FBG 應變計數(shù)據(jù)為樣本，根據(jù)FAST 的觀測需求，健康監(jiān)測系統(tǒng)的實際工作時間約為7 天，數(shù)據(jù)長度滿足方法診斷需求。遍歷416 個傳感器數(shù)據(jù)庫，共成功提取317 個測點對應特征值，異常診斷結(jié)果如圖6 所示，特征值與異常指數(shù)見表3。

圖6 傳感器異常診斷結(jié)果Fig.6 Result of sensor abnormal diagnosis

表3 故障傳感器的特征及異常指數(shù)Tab.3 Eigen values and abnormal indexes of fault sensors

診斷共識別出4 個故障測點與14 個異常測點。人工對故障與異常測點進行核查，缺失的99個測點未獲取數(shù)據(jù)的原因為光纖光柵解調(diào)儀設(shè)備離線檢修，無監(jiān)測數(shù)據(jù)。故障測點中，38 與265號測點因大量傳感器數(shù)據(jù)返回錯誤被丟棄，傳感器數(shù)據(jù)長度遠低于正常水平。163 與305 號測點為數(shù)據(jù)存在大量跳變，引起多項特征離群，其中305 號測點更為嚴重。14 個異常測點的主要問題為監(jiān)測信號中存在較強的高頻信號干擾，明顯異于其他正常傳感器的平滑信號。

由此可見，診斷結(jié)果對傳感器的工作狀態(tài)有較強的指示作用，可根據(jù)診斷結(jié)果對故障傳感器通道進行關(guān)閉并維修，防止錯誤數(shù)據(jù)影響健康監(jiān)測系統(tǒng)的監(jiān)控功能。而對異常傳感器，需進行進一步排查，例如高頻干擾，多由測點附近的機械振動或光纖通道干擾造成，應及時排除，防止傳感器狀態(tài)進一步惡化。

6 結(jié) 論

本文基于健康監(jiān)測系統(tǒng)中光纖光柵應變計原始輸出波長信號具有統(tǒng)一規(guī)律的原理，使用信號處理的方法提取歷史數(shù)據(jù)中心反映傳感器狀態(tài)的特征值，建立特征列表。通過循環(huán)迭代尋找一維特征值的聚合中心，并根據(jù)特征值的離群程度定義綜合異常指數(shù)，實現(xiàn)了傳感器的故障與異常診斷。仿真結(jié)果表明，當故障測點占總測點數(shù)量的20%以內(nèi)時，可有效識別故障；當異常測點過多致使特征值過于離散時，可進一步調(diào)節(jié)異常與故障系數(shù)。利用FAST 健康監(jiān)測系統(tǒng)2019 年9 月份的數(shù)據(jù)對416 個FBG 傳感器進行了診斷，共識別出4 個故障測點，14 個異常測點，人工核查無誤。

該方法是一種自動化的傳感器異常診斷方法，不需要大量先驗知識的訓練，根據(jù)傳感器的歷史數(shù)據(jù)對傳感器狀態(tài)進行評估，對多種故障有較強的指示作用。它能夠及時、自動化地識別故障傳感器測點，節(jié)省大量人力，提高系統(tǒng)可靠性。不過，傳感器的具體故障與異常原因仍需結(jié)合人工核查，在故障類別的自動判別上存在不足，需進一步研究。