本發明專利技術涉及一種短時閃變實時計算方法及裝置，屬于電力系統測量技術領域。本發明專利技術只需要將短時閃變計算的所需的累計概率轉換成與自適應方法適配的概率值，將該概率值作為分位數概率，再利用得到的分位數概率、瞬時閃變度值的所屬范圍以及瞬時閃變度值計算各個分位數概率對應的分位數值；最后基于得到的分位數概率以及對應的分位數值就可以進行短時閃變實時計算。本發明專利技術不需要對檢測時間內的數據進行累積概率分布進行統計，也不需分段線性插值，只需要利用瞬時閃變度值的所屬范圍以及瞬時閃變度值計算各個分位數概率對應的分位數值即可，計算方式簡單，也不需要存儲大量數據；同時能夠根據瞬時閃變值的動態變化自適應調整估計值，保證估計精度。
【技術實現步驟摘要】
本專利技術涉及一種短時閃變實時計算方法及裝置，屬于電力系統測量。
技術介紹
1、電壓波動造成燈光照度不穩定(燈光閃爍)的人眼視感反應稱為閃變。換言之，閃變反映了電壓波動引起的燈光閃爍對人的視感產生的影響。短時閃變值(short?termseverity,pst)是衡量短時間(若干分鐘)內閃變強弱的一個統計量值。通常，pst＝1為閃變引起視感刺激性的限值。閃變作為電能質量監測中的一項重要技術指標，可反應電壓的實時波動情況，對于閃變的研究越來越多，但基本上都是需要緩存一段時間的瞬時視感度數據，然后進行統計計算，譬如短時閃變值需要累計10分鐘計算一次，按400hz采樣率計算，需要緩存240000個數據，不僅對存儲空間要求比較大，而且計算量大影響計算效率。
2、為此，有人提出了能夠減小存儲量的短時閃變計算方法，申請公布號為cn102928703a的中國專利申請文件公開一種測量短時閃變的方法及系統，該方法通過對瞬時閃變值進行分級，記錄各個級別下的瞬時閃變值的個數，可以在數據存儲空間只存儲各個級別下的瞬時閃變值的個數，能夠減少計算過程中需要存儲的數據量。該方案將計算短時閃變所需存儲的數據個數由240000降至分級數(按iec標準至少分為64級)。但是該方法需預先設置分級門限，對預設時間內的所有瞬時閃變值逐個按照預設門限進行分級，需存儲分級的每一級個數，即實時計算瞬時閃變值分布直方圖，而現實環境中的閃變值往往是無法預知的，一旦預設門限值不符合瞬時閃變值分布規律，計算得到的分布直方圖也將失真，無法得到準確的短時閃變值。
/>技術實現思路
1、本專利技術的目的是提供一種短時閃變實時計算方法及裝置，以解決目前短時閃變實時計算中需預先設置分級門限導致的計算精度差的問題。
2、本專利技術為解決上述技術問題而提供一種短時閃變實時計算方法，該方法包括：
3、1)獲取檢測時間內的瞬時閃變度值；
4、2)將短時閃變計算所需的累計概率轉換成與自適應方法適配的概率值，稱為分位數概率，各分位數概率等于1減去對應的累計概率；
5、3)利用得到的分位數概率、瞬時閃變度值的所屬范圍以及瞬時閃變度值計算各個分位數概率對應的分位數值；
6、4)將各分位數概率與對應的分位數值相乘后求和，以得到檢測時間內的短時閃變值。
7、進一步地，所述步驟3)中各個分位數概率對應的分位數值采用的計算公式為：
8、
9、其中，pm(n)為第m個分位數概率對應的分位數值，p'm為第m個分位數概率，r(n)為第n個采樣值對應的瞬時閃變度值的所屬范圍，s(n)為第n個采樣值對應的瞬時閃變度值。
10、進一步地，瞬時閃變度值的所屬范圍是實時更新得到的，采用的計算公式為：
11、r(n)＝rmax(n)-rmin(n)
12、
13、其中rmax(n)、rmin(n)分別是瞬時閃變度值的所屬范圍的上、下界，r(n)為第n個采樣值對應的瞬時閃變度值的所屬范圍大小，s(n)為第n個采樣值對應的瞬時閃變度值，α和β均為大于0小于1的系數。
14、進一步地，所述步驟4)中檢測時間內的短時閃變值采用的計算公式為：
15、
16、其中p0′.999、p0′.99、p0′.97、p0′.9和p0′.5分別為檢測時間內累計概率為0.1％、1％、3％、10％和50％對應的分位數概率，p0.999、p0.99、p0.97、p0.9和p0.5分別為檢測時間內分位數概率p0′.999、p0′.99、p0′.97、p0′.9和p0′.5對應的分位數值。
17、進一步地，檢測時間內的瞬時閃變度值是由檢測時間內的電壓采樣信號經過平方化處理、帶通濾波、視感度加權濾波和平滑加權濾波得到。
18、進一步地，所述的帶通濾波采用的帶通濾波器由截止頻率為0.05hz的高通濾波器和截止頻率為35hz的低通濾波器組成。
19、進一步地，視感度加權濾波采用的視感度加權濾波器的傳遞函數為：
20、
21、其中，w1＝2π*9.15494，w2＝2π*2.27979，w3＝2π*1.22535，w4＝2π*21.9，k＝1.74802，λ1＝2π*4.05981。
22、進一步地，所述的視感度加權濾波器采用matlab的bilinear函數轉換為數字濾波器實現。
23、進一步地，所述的平滑加權濾波采用的傳遞函數為：
24、
25、其中τ選用300ms。
26、本專利技術還提供了一種短時閃變實時計算裝置，包括有處理器，所述處理器用于執行計算機程序指令以實現本專利技術的短時閃變實時計算方法。
27、本專利技術的有益效果是:作為改進型專利技術創造，本專利技術只需要將短時閃變計算的所需的累計概率轉換成與自適應方法適配的概率值，將該概率值作為分位數概率，再利用得到的分位數概率、瞬時閃變度值的所屬范圍以及瞬時閃變度值計算各個分位數概率對應的分位數值；最后基于得到的分位數概率以及對應的分位數值就可以進行短時閃變實時計算。本專利技術不需要對檢測時間內的數據進行累積概率分布進行統計，也不需分段線性插值，只需要利用瞬時閃變度值的所屬范圍以及瞬時閃變度值計算各個分位數概率對應的分位數值即可，計算方式簡單，也不需要存儲大量數據。同時能夠根據瞬時閃變值的動態變化自適應調整估計值，保證估計精度。
本文檔來自技高網...
【技術保護點】
1.一種短時閃變實時計算方法，其特征在于，該方法包括：
2.根據權利要求1所述的短時閃變實時計算方法，其特征在于，所述步驟3)中各個分位數概率對應的分位數值采用的計算公式為：
3.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，瞬時閃變度值的所屬范圍是實時更新得到的，采用的計算公式為：
4.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，所述步驟4)中檢測時間內的短時閃變值采用的計算公式為：
5.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，檢測時間內的瞬時閃變度值是由檢測時間內的電壓采樣信號經過平方化處理、帶通濾波、視感度加權濾波和平滑加權濾波得到。
6.根據權利要求5所述的短時閃變實時計算方法，其特征在于，所述的帶通濾波采用的帶通濾波器由截止頻率為0.05hz的高通濾波器和截止頻率為35hz的低通濾波器組成。
7.根據權利要求5所述的短時閃變實時計算方法，其特征在于，視感度加權濾波采用的視感度加權濾波器的傳遞函數為：
8.根據權利要求7所述的短時閃變實時計算方法，其特征在于，所述的視感度加權濾波器采用matlab的bilinear函數轉換為數字濾波器實現。
9.根據權利要求5所述的短時閃變實時計算方法，其特征在于，所述的平滑加權濾波采用的傳遞函數為：
10.一種短時閃變實時計算裝置，包括有處理器，其特征在于，所述處理器用于執行計算機程序指令以實現如權利要求1-9任一項所述的短時閃變實時計算方法。
...
【技術特征摘要】
1.一種短時閃變實時計算方法，其特征在于，該方法包括：
2.根據權利要求1所述的短時閃變實時計算方法，其特征在于，所述步驟3)中各個分位數概率對應的分位數值采用的計算公式為：
3.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，瞬時閃變度值的所屬范圍是實時更新得到的，采用的計算公式為：
4.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，所述步驟4)中檢測時間內的短時閃變值采用的計算公式為：
5.根據權利要求1或2所述的短時閃變實時計算方法，其特征在于，檢測時間內的瞬時閃變度值是由檢測時間內的電壓采樣信號經過平方化處理、帶通濾波、視感度加權濾波和平滑加權濾波得到。
6.根據權利要求5所述的短時...
【專利技術屬性】
技術研發人員：陳飛飛，吳義文，黃梅瑩，孫胤杰，樓紅偉，
申請(專利權)人：深圳智微電子科技股份有限公司，
類型：發明
國別省市：
還沒有人留言評論。發表了對其他瀏覽者有用的留言會獲得科技券。