周 宇 ，馬宏忠 ，李 凱 ，許洪華 ，蔣本洲

（1.河海大學(xué) 能源與電氣學(xué)院，江蘇 南京 211100；2.國網(wǎng)江蘇省電力公司 南京供電公司，江蘇 南京 210019；3.南京立業(yè)電力變壓器有限公司 電氣試驗站，江蘇 南京 211800）

0 引言

變壓器是電力傳輸?shù)臉屑~設(shè)備。運輸、安裝、短路電流沖擊等會導(dǎo)致變壓器繞組壓緊力減?。?-2］，使得繞組松動［3］。繞組壓緊力減小會減弱繞組的抗短路沖擊能力，嚴(yán)重威脅著變壓器的正常運行［4-6］。

變壓器油箱振動是對繞組振動和鐵芯振動等振動源的響應(yīng)。研究表明：鐵芯和繞組分別受到磁致伸縮力［7］和電動力的作用，產(chǎn)生的振動分別正比于電壓的平方［8］和電流的平方［2］。 振動的傳播介質(zhì)為變壓器油、固體構(gòu)件等。振動傳遞至油箱表面的過程中，其幅值衰減［9］。繞組松動、匝間短路等會影響變壓器振動信號，因此振動信號可以用于變壓器故障診斷。目前，國內(nèi)外可用于評估變壓器繞組壓緊狀態(tài)的有效方法并不多［10-14］。文獻(xiàn)［3］指出變壓器振動幅值變化與固有頻率有關(guān)。文獻(xiàn)［13］深入研究了空載狀況下油箱表面振動加速度幅值和繞組壓緊力之間的關(guān)系，提出利用空載時油箱表面測點振動加速度信號隨壓緊力減小的變化檢測繞組壓緊狀態(tài)。研究表明［1］，在短路情況下，預(yù)緊力減小時，器身上的振動有一個先減小，然后隨預(yù)緊力進(jìn)一步減小而增大的過程。文獻(xiàn)［15］中的實驗是在短路情況下進(jìn)行的，實驗電壓較低而電流較大，所以振動源主要是繞組振動，在其后的分析中考察了振動信號與電流之間的線性關(guān)系以及繞組松動后振動的變化。

在對變壓器振動信號進(jìn)行時頻域分析時會遇到頻譜泄漏、噪聲影響等問題。阻尼、剛度等參數(shù)難以準(zhǔn)確測定，所以建立動力學(xué)模型也存在較大難度。相空間重構(gòu)是根據(jù)有限的數(shù)據(jù)來重構(gòu)吸引子以研究系統(tǒng)動力學(xué)的方法。繞組壓緊狀態(tài)影響著變壓器繞組振動傳遞能力、變壓器固有頻率以及漏磁場分布等，進(jìn)而影響振動。本文通過相空間重構(gòu)方法分析振動信號，目的是識別松動狀態(tài)下振動信號的相空間軌線特征，該特征可以體現(xiàn)繞組松動對振動信號的影響。本文設(shè)計短路實驗的原因是：短路實驗的振動源主要是繞組振動，而空載實驗的振動源主要是鐵芯振動，相對于鐵芯振動，繞組振動是繞組壓緊狀態(tài)更加直接的反應(yīng)；變壓器運行電壓在額定電壓附近，而電流隨負(fù)荷變化，因此運行中的變壓器振動信號的影響因素主要是電流和繞組狀態(tài)。

1 相空間重構(gòu)

1.1 概念

相空間重構(gòu)的基本思想是：系統(tǒng)中任一分量的演化都是由與之相互作用的其他分量所決定的，相關(guān)分量的信息隱含在任一分量的發(fā)展過程中，即用系統(tǒng)的一個觀察量可以重構(gòu)出原動力系統(tǒng)模型。文獻(xiàn)［16］證明了可以找到一個合適的嵌入維數(shù)重構(gòu)相空間，可以把相空間內(nèi)有規(guī)律的軌跡（吸引子）恢復(fù)出來。

設(shè)有一時間序列｛xi｝（i=1，2，…，N0），對其重構(gòu)的相空間Rm的元素為：

其中，m為重構(gòu)相空間的維數(shù)；τ為延遲時間間隔數(shù)，為延遲時間與采樣時間間隔之比；p=N0-（m-1）τ，為時間序列嵌入相空間的向量個數(shù)；N0為時間序列的數(shù)據(jù)點數(shù)。

1.2 重構(gòu)參數(shù)的確定方法

由重構(gòu)相空間理論［17-20］可知，選取合適的時間延遲τ和嵌入維數(shù)m是進(jìn)行相空間重構(gòu)的關(guān)鍵。

考慮到相圖的可視性，選取m=2或3。被測的振動信號時間序列是有限長的信號并含有環(huán)境、儀器等產(chǎn)生的噪聲，由Takens定理知，τ不能隨意選取。一般根據(jù)經(jīng)驗來選擇τ，其基本思想是使xj和xj＋τ具有某種程度的獨立但又不完全無關(guān)，以便它們在重構(gòu)相空間中作為獨立的坐標(biāo)處理。求取τ的方法有自相關(guān)函數(shù)法、平均位移法、復(fù)自相關(guān)法和互信息法［19-20］等。本文采用平均位移法。

在L2范數(shù)（指先求向量各元素的平方和然后求平方根）下定義平均位移S（m，τ）為：

平均位移S（m，τ）表示重構(gòu)狀態(tài)空間的軌線從狀態(tài)空間主對角線打開的程度。當(dāng)τ較小時，重構(gòu)后的吸引子會被壓縮在主對角線一帶，而隨著τ的增大，吸引子的形狀會逐漸展開。最后求得的τ應(yīng)該保證重構(gòu)的相空間軌線從相空間的主對角線向外擴展，而又不發(fā)生折疊現(xiàn)象。文獻(xiàn)［17］建議取S（m，τ）-τ曲線的斜率第一次減少到小于初始斜率40%時的τ為最佳延遲時間間隔。本文應(yīng)用該方法時，求取最佳延遲時間間隔后進(jìn)行了小范圍的調(diào)整，以達(dá)到最佳的重構(gòu)效果。

2 實驗系統(tǒng)與實驗方法

對南京立業(yè)電力變壓器有限公司的一臺110 kV變壓器進(jìn)行不同壓緊狀態(tài)下的短路實驗，測取油箱表面的振動加速度信號。

2.1 實驗系統(tǒng)

實驗系統(tǒng)包括發(fā)電機、變壓器、振動加速度傳感器、采集儀和液壓系統(tǒng)等。

a.發(fā)電機用于向變壓器高壓側(cè)繞組提供50 Hz交流電。

b.變壓器的型號為SFZ10-31500／110，聯(lián)接組標(biāo)號為YNd11，低壓側(cè)額定電壓為10.5 kV，額定電流IN為1732 A。實驗時，令變壓器的低壓側(cè)繞組三相短路。

c.振動加速度傳感器的型號為JF2020，其輸出電壓與振動加速度的對應(yīng)關(guān)系為20 mV／g（g為重力加速度）。振動加速度傳感器靠磁座吸附在變壓器油箱表面。

d.采集儀選用Nicolet采集儀，采樣頻率設(shè)置為10 kHz，每次采樣時間為10 s。由采集儀向振動加速度傳感器施加2 mA激勵。

e.液壓系統(tǒng)可以通過油壓泵精確控制繞組壓緊力。

2.2 繞組壓緊力設(shè)置方法

根據(jù)實驗方案人為地對繞組壓緊力大小進(jìn)行設(shè)置。對變壓器進(jìn)行抽油、吊罩操作，然后通過液壓系統(tǒng)控制繞組壓緊力，實現(xiàn)繞組不同壓緊狀態(tài)的設(shè)置。設(shè)置過程為如下。

a.使用扳手?jǐn)Q松壓緊螺母，然后將液壓千斤頂放置在繞組和結(jié)構(gòu)件之間；液壓表顯示液壓千斤頂所受壓力，28 MPa對應(yīng)繞組未松動狀態(tài)的額定壓緊力，0 MPa對應(yīng)松動狀態(tài)的繞組壓緊力。

b.控制液壓，待壓力穩(wěn)定后，擰緊壓緊螺母，使得液壓千斤頂可以抽出；然后再將變壓器罩恢復(fù)、注油。

針對設(shè)計的不同壓緊力方案，每次均按上述過程進(jìn)行抽油、吊罩、調(diào)整壓緊力等操作。本實驗共設(shè)置3種壓緊狀態(tài)：壓緊（三相繞組壓緊力均為28 MPa）、松動（三相繞組壓緊力均為0 MPa）和B相半松狀態(tài)（B相壓緊力為14 MPa而其余相繞組壓緊力為28 MPa）。在尋找繞組松動特征時，選取振動信號差別較大的壓緊狀態(tài)和松動狀態(tài)下的振動信號進(jìn)行分析，找出特征并制定繞組松動識別方法。然后利用B相半松狀態(tài)的振動信號進(jìn)行繞組松動識別方法的驗證。

2.3 振動信號測量方法

圖1為振動信號測點布置圖，振動信號測點布置在圖中的1、2、3號位置，位于變壓器油箱頂面，測點正下方分別對應(yīng)著C、B、A三相繞組。繞組和鐵芯通過結(jié)構(gòu)件與油箱底部相連，與油箱頂部緊密接觸。同一實驗條件（同一種壓緊狀態(tài)和同一種電流狀態(tài)）下的實驗數(shù)據(jù)測3次。電流狀態(tài)有6種，分別為20%額定電流、40%額定電流、60%額定電流、80%額定電流、100%額定電流和110%額定電流。

3 振動信號相空間重構(gòu)

首先對處于油箱頂部中間位置的2號測點的振動信號進(jìn)行相空間重構(gòu)：求取S（m，τ）-τ 曲線，并依據(jù)平均位移法求得最佳延遲時間，繪制不同電流狀態(tài)和不同壓緊狀態(tài)下的相圖。然后運用同樣的方法分析1號測點和3號測點的振動信號。

3.1 振動信號相空間重構(gòu)參數(shù)的確定

本文進(jìn)行實驗數(shù)據(jù)處理的目的是把相空間內(nèi)有規(guī)律的軌跡恢復(fù)出來，通過相圖識別繞組松動缺陷。重構(gòu)參數(shù)m和τ的選取尤為關(guān)鍵，本文選取2號測點未松動時額定電流下的第1次測量振動信號進(jìn)行計算，確定重構(gòu)參數(shù)。本文取m=2，考察二維相圖（相平面）。由式（2）計算振動信號的平均位移，繪制S（m，τ）-τ 折線圖，如圖2 所示。

圖2 2號測點振動信號 S（m，τ）-τ 折線圖Fig.2 Curve ofS（m，τ） vs.τ for vibration signal of measuring point No.2

根據(jù)平均位移法，計算初始斜率，當(dāng)曲線斜率第一次下降到初始斜率的40%以下時，取對應(yīng)的τ值為最佳延遲時間間隔。由圖2可以看出在τ=20時曲線斜率約為0，所以最佳延遲時間間隔在20以內(nèi)。為方便計算，運用MATLAB曲線擬合工具箱中的Fourier對延遲時間間隔在20以內(nèi)的S（m，τ）-τ散點圖進(jìn)行擬合。擬合結(jié)果如圖3所示。擬合曲線的解析式為：

其中，ω=0.09377。

圖3S（m，τ）-τ 散點擬合Fig.3 Fitted curve ofS（m，τ） vs.τ

從x=1開始計算導(dǎo)數(shù)f′（x），由計算結(jié)果可知，使得 f′（x）／f′（1）＜40% 成立的第一個整數(shù)點為 x=10（f′（10）／f′（1）=38.12%）。 附近點的斜率與初始斜率的比值為：f′（9）／f′（1）=41.17%，f′（11）／f′（1）=34.1%。 MATLAB 顯示結(jié)果為：R-square=0.9816，該值越接近1，表明擬合越理想；均方根誤差（RMSE）為0.001318，較小。根據(jù)擬合曲線，求得最佳延遲時間間隔τ為10，因為采樣頻率是10 kHz，所以對應(yīng)的延遲時間為 0.001 s。

3.2 振動信號相空間重構(gòu)

將2號測點在不同電流狀態(tài)和不同壓緊狀態(tài)下的振動信號按m=2和τ=10進(jìn)行相空間重構(gòu)，發(fā)現(xiàn)同一實驗條件下3次測量振動信號得到的相平面幾乎相同，本文取第1次測量結(jié)果進(jìn)行分析，如圖4所示（由于篇幅所限，第2次和第3次測量信號對應(yīng)的相平面本文不再示出）。對于本文中的相空間重構(gòu)，其主對角線是 y=x［22］。 圖4（a）、（b）中，波形由上至下分別為 110%、100%、80%、60%、40%、20% 額定電流下的相平面，橫軸為振動向量v（i），縱軸為振動向量v（i+τ）。為了防止圖形隨橫縱坐標(biāo)的比例而變化，將圖4的橫軸和縱軸的尺度比例設(shè)置成1∶1，下文的所有相圖均采用該設(shè)置。

圖4 2號點振動信號的延遲10相平面Fig.4 Phase planes of vibration signals of measuring point No.2 whenτ=10

從圖4中可以得出以下結(jié)論。

a.松動狀態(tài)和壓緊狀態(tài)下振動信號的相平面不同，這是因為繞組松動后影響了變壓器繞組的振動傳遞能力、變壓器固有頻率以及漏磁場分布等，導(dǎo)致松動狀態(tài)和壓緊狀態(tài)下的振動信號不同，其重構(gòu)出的相平面也不相同。

b.觀察松動狀態(tài)下，電流為80%、100%、110%的額定電流時的相平面軌線，可以發(fā)現(xiàn)：軌線近似為空心橢圓形狀，即每個相點均在空心橢圓上并且每個相點與下一個相點的連線也在空心橢圓上；振動信號重構(gòu)得到的相平面軌線從空間主對角線打開，證明選取的維數(shù)和最佳延遲時間參數(shù)效果較好；壓緊狀態(tài)下的振動信號相平面軌線交叉混疊，未沿著主對角線打開。根據(jù)動力學(xué)理論［21］，系統(tǒng)振動所得相平面軌線的閉合性可以用于判斷系統(tǒng)振動的周期性，軌線是否有畸變可以用于判定系統(tǒng)處于線性振動或非線性振動狀態(tài)。相平面軌線是閉合的橢圓，表明變壓器短路實驗中的穩(wěn)態(tài)振動是周期振動，其振動基頻為100 Hz［2］；橢圓形狀有畸變，表明變壓器振動系統(tǒng)是非線性系統(tǒng)。非線性關(guān)系存在于磁感應(yīng)強度與磁場強度之間以及電流與磁通之間。

c.觀察20%、40%、60%的額定電流狀態(tài)下的相平面可見，松動狀態(tài)和壓緊狀態(tài)下的軌線均未沿著主對角線打開，2種狀態(tài)下的相平面區(qū)分度不大。這是因為較小電流下，振動信號不強，隨機性較大，信噪比小，容易受噪聲影響，其振動信號與大電流下的振動信號差別較大，嵌入維數(shù)和最佳時間延遲由額定電流狀態(tài)下的振動信號得到，運用該參數(shù)對較小電流狀態(tài)下（60%及以下的額定電流）的振動信號進(jìn)行相空間重構(gòu)，并不能達(dá)到滿意的效果。

因為所處位置均在頂面，所以1、3號測點的振動信號與2號測點振動信號存在一定的相似性。將由2號測點的振動信號得到的重構(gòu)參數(shù)運用到1、3號測點振動信號的相空間重構(gòu)中，發(fā)現(xiàn)重構(gòu)效果并不理想，原因是2號測點振動信號與1號、3號測點的振動信號也存在一定的差異性，存在差異的原因是分接開關(guān)、儲油罐等結(jié)構(gòu)對3個測點的振動信號有不對稱的影響。所以在τ=10附近尋找合適的時間延遲。首先進(jìn)行初選：以10為中心，每隔5個時間延遲間隔確定1個延遲進(jìn)行相空間重構(gòu)；當(dāng)呈現(xiàn)出有規(guī)律的軌線時，確定當(dāng)前延遲為初選延遲。然后進(jìn)行細(xì)選：以初選延遲為中心，每隔1個時間延遲間隔確定1個延遲進(jìn)行相空間重構(gòu)。如果采用某個延遲后，松動狀態(tài)下1、3號測點的振動信號在盡可能多的電流狀態(tài)下的相平面軌線可以沿主對角線打開而又盡量不發(fā)生折疊現(xiàn)象，則確定當(dāng)前延遲為最佳時間延遲。最終得到1號測點的最佳時間延遲為τ=6，3號測點的最佳時間延遲為τ=16。同一實驗條件下的3次測量振動信號得到的相平面幾乎相同，取第1次測量結(jié)果進(jìn)行分析。1、3號測點的振動信號相平面分別如圖5、6所示，各波形與電流、壓緊狀態(tài)的對應(yīng)關(guān)系與圖4相同。限于篇幅，1號測點的第2次和第3次測量信號對應(yīng)的相平面，以及3號測點的第2次和第3次測量信號對應(yīng)的相平面本文不再示出。

圖5 1號點振動信號的延遲6相平面Fig.5 Phase planes of vibration signals of measuring point No.1 whenτ=6

按照對圖4的分析方法，分析圖5、6，分析結(jié)果如表1所示。

表1表明：對于在100%、110%額定電流狀態(tài)下運行的變壓器，可以通過測繪其頂面3個測點振動信號相平面，根據(jù)識別的松動狀態(tài)相平面特征（軌線近似為空心橢圓），對繞組壓緊狀態(tài)進(jìn)行判斷；在80%、100%、110%額定電流狀態(tài)下運行的變壓器，可以通過測繪其頂面2號和3號測點振動信號相平面，判斷繞組壓緊狀態(tài)。

4 應(yīng)用振動信號相空間重構(gòu)方法識別變壓器繞組松動

4.1 方法內(nèi)容

依據(jù)前文對差別較大的壓緊狀態(tài)和松動狀態(tài)下振動信號的分析可知，能夠根據(jù)相圖特征識別變壓器繞組松動，步驟如下。

a.尋找某種類型（容量相同且結(jié)構(gòu)類似）變壓器的繞組松動對應(yīng)的相圖特征。對該類型的某臺變壓器進(jìn)行短路實驗，測取油箱上表面的振動信號。為保證采集的信號中變壓器振動信號是穩(wěn)定的主要成分，電流應(yīng)該設(shè)置為80%及以上的額定電流。

圖6 3號點振動信號的延遲16相平面Fig.6 Phase planes of vibration signals of measuring point No.3 whenτ=16

表1 變壓器油箱頂面3個測點相平面分析Table1 Analysis of phase planes for 3 measuring points

b.根據(jù)平均位移法求取最佳延遲時間，在小范圍內(nèi)調(diào)整最佳延遲時間，使得繪制出的松動狀態(tài)下的振動信號2維或3維相圖具有明顯的特征。

c.根據(jù)發(fā)現(xiàn)的相圖特征，監(jiān)測同種類型的變壓器是否存在繞組松動缺陷。若圖形滿足相圖特征，則判斷變壓器存在繞組松動缺陷。

步驟a、b僅針對某一類型的變壓器相空間特征的尋找過程，對不同類型的變壓器應(yīng)有不同的相圖數(shù)據(jù)庫。步驟c為變壓器的監(jiān)測過程。

對于與本文實驗對象同類型的變壓器，根據(jù)表1，在80%、100%、110%額定電流下運行的該類型變壓器繞組松動的識別方法為：對于在100%、110%額定電流下運行的變壓器，若其頂面1、2、3號測點的振動信號相平面近似為空心橢圓，則判定其存在繞組松動缺陷；對于在80%、100%、110%額定電流下運行的變壓器，若其頂面2、3號測點的振動信號相平面近似為空心橢圓，則判定其存在繞組松動缺陷；若變壓器頂面振動信號相平面未沿主對角線打開，未呈現(xiàn)出空心橢圓形狀，則判定其不存在繞組松動缺陷。

4.2 方法驗證

利用B相半松狀態(tài)下的振動信號對繞組松動的識別方法進(jìn)行驗證。選用表1中的重構(gòu)參數(shù)對B相半松狀態(tài)下1、2、3號測點的振動信號進(jìn)行相空間重構(gòu)。在同一實驗條件下測取3次振動信號，取第1次測量結(jié)果進(jìn)行分析，見圖7，圖中從左至右分別對應(yīng)110%、100%、80%額定電流。B相半松狀態(tài)下1、2、3號測點第2次和第3次測量信號的相平面與第1次測量信號相平面幾乎相同，限于篇幅，不再顯示。

圖7 B相半松狀態(tài)下1—3號測點振動信號相平面Fig.7 Phase planes of vibration signals of 3 measuring points，when phase-B winding is half loosened

由圖7可以得出以下結(jié)論。

a.1、2、3號測點的軌線在一定程度上有空心橢圓的形狀，與圖4（a）、圖5（a）和圖6（a）中的交叉混疊的繞組壓緊狀態(tài)軌線之間存在差異，可以判斷測試變壓器存在繞組松動缺陷。

b.2號測點相平面軌線沿主對角線打開的程度大于1、3號測點，原因是2號測點更靠近發(fā)生松動的B相繞組，其振動信號與壓緊狀態(tài)下振動信號的差別更大。根據(jù)這一特點可以進(jìn)一步實現(xiàn)繞組松動相的定位。

c.對比圖4（b）、圖5（b）和圖6（b）的繞組松動狀態(tài)軌線和B相半松狀態(tài)軌線（圖7）可以發(fā)現(xiàn)，松動（三相繞組全部松動）狀態(tài)軌線比半松動狀態(tài)軌線的形狀特征更加明顯，畸變較小，這為進(jìn)一步區(qū)分繞組松動程度提供了思路。

5 結(jié)論

繞組壓緊狀態(tài)影響著變壓器繞組的振動傳遞能力、固有頻率以及漏磁場分布，進(jìn)而影響振動。繞組松動前后的振動信號重構(gòu)相圖的差異體現(xiàn)了繞組松動對振動信號的影響。

恰當(dāng)?shù)剡x取重構(gòu)參數(shù)m和τ，可以把振動信號相空間內(nèi)有規(guī)律的軌線恢復(fù)出來。結(jié)合試湊方法，應(yīng)用平均位移法求取最佳時間延遲間隔。松動狀態(tài)下，較大電流（80%、100%、110%額定電流）下的軌線為一閉合的空心畸變橢圓，這表明變壓器振動是非線性周期振動。

本文中的實驗變壓器并未裝載風(fēng)扇，對于帶風(fēng)冷的變壓器，其振動信號含有風(fēng)扇振動成分，與本文測取的振動信號存在差異，需要重新選定重構(gòu)參數(shù)。本文實驗對象為鐘罩式結(jié)構(gòu)的大型變壓器，而小型變壓器（如10 kV變壓器）為吊芯式結(jié)構(gòu)，振動信號抗干擾能力較弱，其振動信號測點布置以及重構(gòu)效果有待進(jìn)一步研究。

參考文獻(xiàn)：

［1］謝坡岸.振動分析法在電力變壓器繞組狀態(tài)監(jiān)測中的應(yīng)用研究［D］.上海：上海交通大學(xué)，2008.XIE Poan.Study on application of vibration analysis to the condition monitoring of power transformers windings［D］.Shanghai：Shanghai Jiao Tong University，2008.

［2］馬宏忠，耿志慧，陳楷，等.基于振動的電力變壓器繞組變形故障診斷新方法［J］.電力系統(tǒng)自動化，2013，37（8）：1-7.MA Hongzhong，GENG Zhihui，CHEN Kai，et al.A new fault diagnosismethod forpowertransformerwinding deformation based on vibration［J］.Automation of Electric Power Systems，2013，37（8）：1-7.

［3］陸瑾.應(yīng)用機械振動法對變壓器繞組狀態(tài)進(jìn)行在線測試的研究［D］.上海：上海交通大學(xué)，2008.LU Jin.The studyon the on-linediagnosisoftransformer winding status with the method of mechanical vibration［D］.Shanghai：Shanghai Jiao Tong University，2008.

［4］姚森敬，歐陽旭東，林春耀.電力變壓器繞組變形診斷分析［J］.電力系統(tǒng)自動化，2005，29（18）：95-98.YAO Senjing，OUYANG Xudong，LIN Chunyao.Analysison winding deformation diagnosing of power transformers［J］.Automation of Electric Power Systems，2005，29（18）：95-98.

［5］馬宏彬，何金良，陳青恒.500 kV單相電力變壓器的振動與噪聲波形分析［J］.高電壓技術(shù)，2008，34（8）：1599-1604.MA Hongbin，HE Jinliang，CHEN Qingheng.Vibration and sound waveform analysis of 500 kV single phase power transformer［J］.High Voltage Engineering，2008，34（8）：1599-1604.

［6］劉軍，張安紅.電力變壓器繞組短路動穩(wěn)定能力的仿真和評估［J］.變壓器，2012，49（6）：14-25.LIU Jun，ZHANG Anhong.Simulation and evaluation of short circuit dynamic stability of windings in power transformer［J］.Transformer，2012，49（6）：14-25.

［7］祝麗花，楊慶新，閆榮格，等.考慮磁致伸縮效應(yīng)電力變壓器振動噪聲的研究［J］.電工技術(shù)學(xué)報，2013，28（4）：1-6.ZHU Lihua，YANG Qingxin，YAN Rongge，etal.Research on vibration and noise ofpowertransformer cores including magnetostriction effects［J］.Transaction of China Electrotechnical Society，2013，28（4）：1-6.

［8］馬宏忠，趙宏飛，陳楷，等.基于振動的變壓器鐵芯松動判定方法［J］.電力系統(tǒng)自動化，2013，37（14）：101-106.MA Hongzhong，ZHAO Hongfei，CHEN Kai，etal.Determine method for transformer core looseness based on vibration ［J］.Automation of Electric Power Systems，2013，37（14）：101-106.

［9］沈培鋒，趙宏飛，李凱，等.基于振動的變壓器鐵芯松動定位［J］.中國電力，2015，47（9）：66-70.SHEN Peifeng，ZHAO Hongfei，LI Kai，etal.Positioning of transformer core loosening based on the vibration signals ［J］.Electric Power，2015，47（9）：66-70.

［10］SANZ-BOBI M A，GARCIA-CERRADA A，PALACIOS R，et al.Experiences learned from the on-line internal monitoring of the behavior of a transformer［C］∥Proceedings of the IEEE International Electric Machines and Drives Conference Record.Milwaukee，WI，USA：IEEE，1997：TC3-11.1-TC3-11.3.

［11］GARCA B，BURGOS J C，áGEL A.Winding deformations detection in power transformers by tank vibrations monitoring［J］.Electric Power Systems Research，2005，74（1）：129-138.

［12］BERLER Z，GOLUBEV A，RUSOV V，et al.Vibro-acoustic of transformer clamping pressure monitoring［C］∥Conference Record of the 2000 IEEE International Symposium on Electrical Insulation.Anaheim，CA，USA：IEEE，2000：263-266.

［13］謝坡岸，金之儉，饒柱石，等.振動法檢測空載變壓器繞組的壓緊狀態(tài)［J］.高電壓技術(shù)，2007，33（3）：188-189.XIE Poan，JIN Zhijian，RAO Zhushi，et al.Windings clamped state detected by vibration method on unloadeded transformer［J］.High Voltage Engineering，2007，33（3）：188-189.

［14］程錦，李延沐，汲勝昌，等.振動法在線監(jiān)測變壓器繞組及鐵心狀況［J］.高電壓技術(shù)，2005，31（4）：43-48.CHENG Jin，LIYanmu，JIShengchang，etal.Application of vibration method on monitoring the winding and core condition of transformer［J］.High Voltage Engineering，2005，31（4）：43-48.

［15］趙宏飛，馬宏忠，陳楷，等.基于振動信號的變壓器繞組松動實驗研究［J］.中國電力，2015，47（1）：13-16.ZHAO Hongfei，MA Hongzhong，CHEN Kai，et al.Experimental study oftransformerwinding loosenessbased on vibration signal［J］.Electric Power，2015，47（1）：13-16.

［16］TAKENS F.On the numerical determination of the dimension of an attractor［J］.Lecture Notes in Mathematics，1985，1125：99-106.

［17］ROSENSTEIN M T，COLLINS JJ，DE LUCA C J，etal.Reconstruction expansion as a geometry-based framework for choosing proper delay times［J］.Nonlinear Phenomena，1994，73（1-2）：82-98.

［18］顧潔.電力系統(tǒng)短期負(fù)荷預(yù)測方法研究［J］.電力自動化設(shè)備，2002，22（2）：20-22.GU Jie.Study on short-term load forecast method for power system［J］.Electric Power Automation Equipment，2002，22（2）：20-22.

［19］安學(xué)利，蔣東翔.風(fēng)力發(fā)電機組運行狀態(tài)的混沌特性識別及其趨勢預(yù)測［J］.電力自動化設(shè)備，2010，30（3）：15-19，24.AN Xueli，JIANG Dongxiang.Chaotic characteristics identification and trend prediction of running state for wind turbine ［J］.Electric Power Automation Equipment，2010，30（3）：15-19，24.

［20］李慧杰，劉亞男，衛(wèi)志農(nóng)，等.基于相關(guān)向量機的短期風(fēng)速預(yù)測模型［J］.電力自動化設(shè)備，2013，33（10）：28-32.LIHuijie，LIU Yanan，WEIZhinong，etal.Short-term wind speed forecasting model based on relevance vector machine［J］.Electric Power Automation Equipment，2013，33（10）：28-32.

［21］司馬文霞，陳莉珺，楊慶，等.基于過電壓時間序列的鐵磁諧振在線建模與反饋控制［J］.高電壓技術(shù)，2014，40（7）：1948-1956.SIMA Wenxia，CHEN Lijun，YANG Qing，et al.Online modeling and feedback control of ferroresonance based on overvoltage time series［J］.High Voltage Engineering，2014，40（7）：1948-1956.

［22］林嘉宇，黃芝平，王躍科，等.語音信號相空間重構(gòu)中時間延遲選擇的改進(jìn)的平均位移法［J］.國防科技大學(xué)學(xué)報，1999，21（3）：59-62.LIN Jiayu，HUANG Zhiping，WANG Yueke，et al.Modification of average displacement method for selection of time-delay in phase space reconstruction of speech signals［J］.Journal of National University of Defense Technology，1999，21（3）：59-62.