本發(fā)明涉及微脈沖波形分析，具體為一種變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法及系統(tǒng)。

背景技術(shù)：

1、通常采用阻尼振蕩來等效開關(guān)操作引起的變電站二次設備端口所受的電磁干擾，并在電磁兼容測試時用其對二次設備進行電磁兼容抗擾度考核。在iec61000-4-18阻尼振蕩波抗擾度測試標準中，明確定義了針對gis開關(guān)引起的二次設備阻尼振蕩波波形參數(shù)，并給出了二次設備抗擾度要求和測試方法。標準僅列出共模電壓的試驗等級，假定差模電壓是共模電壓的一半，即為等級1:0.25kv；等級2:0.5kv；等級3:1kv；等級4:2kv。上述定義中，重復率＝微脈沖個數(shù)/宏脈沖持續(xù)時間：按照pk10衰減到小于pk1值的50％的衰減程度，假設再pk20衰減到0，即認為10個振蕩周期為振蕩波的持續(xù)時間。

2、現(xiàn)有技術(shù)主要統(tǒng)計了峰值電壓和主頻特性，對特高壓等級下的端口騷擾水平和頻率范圍有了大致的了解，在實際應用中，還需對騷擾波形的具體特征進一步分析和研究。因此，亟需一種對騷擾微脈沖波形的自動分析方法。

技術(shù)實現(xiàn)思路

1、鑒于上述現(xiàn)有存在的問題，提出了本發(fā)明。

2、因此，本發(fā)明提供了一種變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法及系統(tǒng)的目的是能夠?qū)Σ杉冸娬舅矐B(tài)電磁騷擾波形數(shù)據(jù)進行自動分析，生成相應的幅值閾值，結(jié)合生成的幅值閾值，采用全局搜索的方式生成微脈沖序列，并對微脈沖序列自動分析，得到相應的波形參數(shù)。

3、為實現(xiàn)上述技術(shù)目的，本發(fā)明提供如下技術(shù)方案：

4、第一方面，本發(fā)明提供了一種變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法，包括：采集變電站瞬態(tài)電磁波形數(shù)據(jù)并提取波形數(shù)據(jù)特征，對所述波形數(shù)據(jù)特征進行分析得到端口類型和背景噪聲幅值，并獲取第一幅值閾值；搜索序列全局最大值，將所述序列全局最大值作為最大微脈沖峰值點，根據(jù)事先確定的距離閾值提取得到所述最大微脈沖峰值點的位置及附近對稱的距離閾值個點的波形數(shù)據(jù)序列，作為其中一個微脈沖提取到二維數(shù)組中；在原序列中將提取出的所述波形數(shù)據(jù)序列置零，并重復提取及置零的步驟直到找出所有幅值大于所述第一幅值閾值的微脈沖；對所述二維數(shù)組中的微脈沖序列進行分析，得到變電站瞬態(tài)電磁波形數(shù)據(jù)對應的微脈沖波形參數(shù)。

5、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：所述端口類型和背景噪聲幅值的獲取包括：

6、基于卷積神經(jīng)網(wǎng)絡構(gòu)建波形識別模型，采集歷史變電站瞬態(tài)電磁波形數(shù)據(jù)，針對每一個變電站瞬態(tài)電磁波形數(shù)據(jù)標注端口類型和背景噪聲幅值，生成樣本數(shù)據(jù)，采用所述樣本數(shù)據(jù)對所述波形識別模型進行訓練；

7、將實時采集的變電站瞬態(tài)電磁波形數(shù)據(jù)導入訓練后的波形識別模型，提取波形數(shù)據(jù)特征，并對所述波形數(shù)據(jù)特征進行分析得到端口類型和背景噪聲幅值；

8、根據(jù)所述端口類型和背景噪聲幅值設置第一幅值閾值。

9、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：大于所述第一幅值閾值的微脈沖的獲取包括：

10、微脈沖的持續(xù)時間為us級，整個宏脈沖的時間為ms級，相鄰兩脈沖間隔時間遠大于微脈沖本身持續(xù)時間，基于此特點搜索序列全局最大值，并將所述序列全局最大值作為最大微脈沖峰值點；

11、根據(jù)事先確定的距離閾值提取得到所述最大微脈沖峰值點的位置及附近對稱的距離閾值個點的波形數(shù)據(jù)序列，將所述波形數(shù)據(jù)序列另存為新的序列，即第一峰值微脈沖，并將所述第一峰值微脈沖作為其中一個微脈沖提取到二維數(shù)組中；

12、將原序列中提取第一峰值微脈沖位置處的數(shù)值全部置零后再次進行全局搜索找出最大值，按上述步驟提取波形數(shù)據(jù)序列另存為新的序列，即第二峰值微脈沖，并將所述第二峰值微脈沖作為其中一個微脈沖提取到二維數(shù)組中；

13、重復置零和提取的步驟直到找出所有幅值大于所述第一幅值閾值的微脈沖。

14、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：將所有幅值大于所述第一幅值閾值的微脈沖全部提取到二維數(shù)組中，所述二維數(shù)組的每一行都是一維微脈沖波形數(shù)據(jù)。

15、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：所述微脈沖波形參數(shù)的獲取包括：

16、將所述二維數(shù)組中的每一行數(shù)據(jù)利用小波變換進行降噪處理后得到新的二維數(shù)組pulse1，計算所述新的二維數(shù)組pulse1中每一行數(shù)據(jù)的變化率；

17、將所述每一行數(shù)據(jù)的變化率與事先確定的變化率閾值進行比較，將所述每一行數(shù)據(jù)的變化率第一個大于所述變化率閾值的絕對時間記為pstart，所述每一行數(shù)據(jù)的變化率最后一個大于所述變化率閾值的絕對時間記為pend，并將pstart和pend的差值作為微脈沖持續(xù)時間；

18、將所述二維數(shù)組中最大峰值微脈沖對應的位置的時間記為tmax，再從所述最大峰值微脈沖對應的位置向零時刻進行搜索，直至搜索到首個與所述最大峰值微脈沖異號的點，將所述首個與所述最大峰值微脈沖異號的點對應的時間記為tmin，再將tmax與tmin之差乘以0.8得到微脈沖上升時間。

19、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：所述微脈沖波形參數(shù)的獲取還包括：

20、將所述二維數(shù)組中所有負值置零得到新的二維數(shù)組pulsep，按照時間搜索滿足條件的點作為波形中過零點zeros，將波形中過零點的位置記為zp，按照所述波形中過零點的位置在原波形數(shù)據(jù)中搜索兩兩過零點間的最大值或者最小值，將所述最大值或者最小值中絕對值大于第二幅值閾值的值作為脈沖峰值pk，并生成脈沖峰值序列；

21、將所述脈沖峰值序列按照時間進行排序形成正負交錯的脈沖峰值序列，將峰值pk1和pk3之間的時間差作為脈沖振蕩周期，則微脈沖衰減即為pk5/pk1，pk10/pk1。

22、作為本發(fā)明所述的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法的一種優(yōu)選方案，其中：所述按照時間搜索滿足條件的點作為波形中過零點zeros，其中所述條件具體包括：

23、搜索到的點本身等于零，且同時前一個位置的點等于零，后一個位置的點大于零；或者搜索到的點本身等于零，且同時前一個位置的點大于零，后一個位置的點等于零。

24、第二方面，本發(fā)明提供了一種變電站瞬態(tài)電磁騷擾微脈沖波形自動分析系統(tǒng)，包括：

25、波形采集模塊，用于采集變電站瞬態(tài)電磁波形數(shù)據(jù)；

26、閾值設置模塊，用于提取變電站瞬態(tài)電磁波形數(shù)據(jù)的特征，對所述波形數(shù)據(jù)特征進行分析得到端口類型和背景噪聲幅值，并獲取第一幅值閾值；

27、全局搜索模塊，用于搜索序列全局最大值，將所述序列全局最大值作為最大微脈沖峰值點，根據(jù)事先確定的距離閾值提取得到所述最大微脈沖峰值點的位置及附近對稱的距離閾值個點的波形數(shù)據(jù)序列，作為其中一個微脈沖提取到二維數(shù)組中；

28、置零搜索單元，用于在原序列中將提取出的所述波形數(shù)據(jù)序列置零，并重復提取及置零的步驟直到找出所有幅值大于所述第一幅值閾值的微脈沖；

29、參數(shù)分析模塊，用于對所述二維數(shù)組中的微脈沖序列進行分析，得到變電站瞬態(tài)電磁波形數(shù)據(jù)對應的微脈沖波形參數(shù)。

30、第三方面，本發(fā)明提供了一種計算設備，包括：

31、存儲器和處理器；

32、所述存儲器用于存儲計算機可執(zhí)行指令，所述處理器用于執(zhí)行所述計算機可執(zhí)行指令，該計算機可執(zhí)行指令被處理器執(zhí)行時實現(xiàn)所述方法的步驟。

33、第四方面，本發(fā)明提供了一種計算機可讀存儲介質(zhì)，其存儲有計算機可執(zhí)行指令，該計算機可執(zhí)行指令被處理器執(zhí)行時實現(xiàn)所述方法的步驟。

34、與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果為：本發(fā)明的變電站瞬態(tài)電磁騷擾微脈沖波形自動分析方法和系統(tǒng)，能夠?qū)Σ杉冸娬舅矐B(tài)電磁騷擾波形數(shù)據(jù)進行自動分析，生成相應的幅值閾值，結(jié)合生成的幅值閾值，采用全局搜索的方式生成微脈沖序列，并對微脈沖序列自動分析，得到相應的波形參數(shù)。