TWI425226B - Method and system for fault detection, identification and location in high-voltage power transmission networks - Google Patents
Method and system for fault detection, identification and location in high-voltage power transmission networks Download PDFInfo
- Publication number
- TWI425226B TWI425226B TW100123915A TW100123915A TWI425226B TW I425226 B TWI425226 B TW I425226B TW 100123915 A TW100123915 A TW 100123915A TW 100123915 A TW100123915 A TW 100123915A TW I425226 B TWI425226 B TW I425226B
- Authority
- TW
- Taiwan
- Prior art keywords
- fault
- voltage
- neural network
- signal
- distribution network
- Prior art date
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
- Locating Faults (AREA)
Description
本發明有關於高壓輸配電網路,特別是有關於高壓輸配電網路之故障偵測、辨識與定位之方法及系統。The invention relates to a high voltage transmission and distribution network, in particular to a method and a system for fault detection, identification and positioning of a high voltage transmission and distribution network.
傳輸線的電力品質在現今是被視為使發電廠能夠達到持續供電給終端使用者的主要因素。由於故障防護系統能避免經濟損失,故對於高壓輸配電網路而言,故障防護系統是基本且重要的。故障防護系統會處理電壓或電流信號,並據此判斷是否有故障產生。若有故障的發生,則故障防護系統會判斷故障種類與故障位置,並決定所需採取之行動,以據此移除輸電線系統中的故障。用於繼電器的傳統傳輸線系統保護方法是在時域與頻域監視電壓和電流信號的失真(distortion)。The power quality of transmission lines is now seen as a major factor in enabling power plants to achieve continuous power supply to end users. Since the fail-safe system can avoid economic losses, the fail-safe system is fundamental and important for high-voltage transmission and distribution networks. The fail-safe system processes the voltage or current signal and determines if a fault has occurred. In the event of a fault, the fail-safe system determines the type of fault and the location of the fault and determines the action required to remove the fault in the power line system. A conventional transmission line system protection method for relays monitors the distortion of voltage and current signals in the time domain and the frequency domain.
防護系統所採用的趨近方式通常可以分為三個類型,其分別為基於量化模式法(quantitative model-based approaches)、基於品質模式法(qualitative model-based approaches)與資料為中心之方法(data-driven approaches)。基於量化模式法與基於品質模式法在模擬研究時具有較好的效果。然而,在故障之電壓與電流信號有雜訊的情況下,上述這兩個方法的效果比較不佳。因此,量測雜訊和系統雜訊是這兩個方法表現好壞的關鍵因素。The approach approach adopted by the protection system can usually be divided into three types, which are based on quantitative model-based approaches, qualitative model-based approaches and data-centric approaches. -driven approaches). Based on the quantitative mode method and the quality model based method, it has a good effect in simulation research. However, in the case where the voltage and current signals of the fault have noise, the above two methods are not effective. Therefore, measuring noise and system noise is a key factor in the performance of these two methods.
在實際的應用上,資料為中心之方法在實施上較具彈性,且可針對不同系統,進行故障偵測與故障分類,因此資料為中心之方法可以使故障防護系統更具能抵抗雜訊之強健性。資料為中心之方法可以例如採用人工智慧技術,因此資料為中心之方法適用於具有未知程度之量測雜訊與系統雜訊的不同系統。In practical applications, the data-centric approach is more flexible in implementation and can be used for fault detection and fault classification for different systems. Therefore, the data-centric approach can make the fault-protection system more resistant to noise. Robust. The data-centric approach can be, for example, artificial intelligence, so the data-centric approach is applicable to different systems with unknown levels of noise and system noise.
在電力傳輸線上所產生的故障可能導致電力供應的中斷,而造成電力連續服務的兩個主要指數(index)的惡化。此兩個主要指數是系統平均中斷頻率指數與系統平均中斷時間指數。為了降低這兩個主要指數,多種方法已經被提出用來偵測故障、分類故障或定位故障位置。當故障發生時,故障防護系統須分類故障類型與辨識問題相位。A failure generated on a power transmission line may cause an interruption in power supply, resulting in deterioration of two main indices of power continuous service. The two main indices are the system average interruption frequency index and the system average interruption time index. In order to reduce these two main indices, various methods have been proposed to detect faults, classify faults or locate fault locations. When a fault occurs, the fail-safe system must classify the fault type and identify the problem phase.
故障類型的特徵可以透過轉換器或濾波器的處理而被得知,濾波器或轉換器例如為傅立業(Fourier)轉換器、卡爾曼濾器(Kalman filter)、碎形(fractal)轉換器或小波(wavelet)轉換器等。然而,故障種類的數量眾多,故傳統故障防護系統可能無法同時保護各種故障類型。簡單地說,傳統故障防護系統在故障發生後的短時間內,可能無法完成故障事件的偵測、故障種類的判定與故障位置的定位。The characteristics of the fault type can be known by the processing of a converter or a filter such as a Fourier converter, a Kalman filter, a fractal converter or a wavelet. (wavelet) converters, etc. However, the number of types of faults is large, so traditional fail-safe systems may not be able to protect various fault types at the same time. Simply put, the traditional fault protection system may not be able to complete the detection of fault events, the determination of the fault type and the location of the fault location within a short time after the fault occurs.
本發明實施例提供一種高壓輸配電網路之故障偵測、辨識與定位之方法,此方法包括以下步驟。首先,將故障之一個週期內的電壓信號與電流信號進行多階層小波轉換(multi-level wavelet transform),並據此產生複數個高頻信號。然後,對複數個高頻信號進行主成分分析(Principle Component Analysis,PCA),並據此產生至少一個信號特徵值(signal characteristic value)接著,使用支援向量機(Support Vector Machine,SVM)分類信號特徵值,以得到故障類型。再來,使用對應故障類型之可適性結構類神經網路(Adaptive Structure Neural Networks,ASNN)偵測故障位置。Embodiments of the present invention provide a method for fault detection, identification, and location of a high voltage transmission and distribution network, and the method includes the following steps. First, a multi-level wavelet transform is performed on a voltage signal and a current signal in one cycle of the fault, and a plurality of high-frequency signals are generated accordingly. Then, Principal Component Analysis (PCA) is performed on a plurality of high frequency signals, and at least one signal characteristic value is generated accordingly, and then the signal characteristics are classified using a Support Vector Machine (SVM). Value to get the type of fault. Then, the fault location is detected using an Adaptive Structure Neural Networks (ASNN) corresponding to the fault type.
本發明實施例還提供一種高壓輸配電網路之故障偵測、辨識與定位之系統,此高壓輸配電網路之故障偵測、辨識與定位之系統包括至少一個量測單元、故障診斷裝置與波形記錄器。量測單元透過匯流排連接電性連接高壓輸配電網路,用以獲得三個相位的電壓與電流信號。故障診斷裝置接收量測單元所獲得的三個相位的電壓與電流信號,且故障診斷裝置進行故障偵測、分類與定位。波形記錄器用以記錄故障發生時的三個相位的電壓與電流信號。The embodiment of the invention further provides a system for fault detection, identification and location of a high voltage transmission and distribution network, wherein the system for detecting, identifying and locating the fault of the high voltage transmission and distribution network comprises at least one measuring unit and a fault diagnosis device. Waveform recorder. The measuring unit is electrically connected to the high-voltage transmission and distribution network through the bus bar connection to obtain voltage and current signals of three phases. The fault diagnosis device receives the voltage and current signals of the three phases obtained by the measurement unit, and the fault diagnosis device performs fault detection, classification and positioning. The waveform recorder is used to record the voltage and current signals of the three phases when the fault occurs.
綜上所述,本發明實施例所提供的高壓輸配電網路之故障偵測、辨識與定位之方法及其系統可偵測故障、判定故障類別與估計故障位置。另外,所述高壓輸配電網路之故障偵測、辨識與定位之方法及其系統自故障發生後,其完成故障偵測、分類與定位的時間遠短於臨界故障清除時間(critical fault clearing time)。In summary, the method, system and system for detecting, identifying, and locating a high voltage transmission and distribution network provided by the embodiments of the present invention can detect a fault, determine a fault category, and estimate a fault location. In addition, the method for detecting, identifying and locating the fault of the high-voltage transmission and distribution network and the system after the fault occurs, the time for completing the fault detection, classification and positioning is much shorter than the critical fault clearing time (critical fault clearing time). ).
為使能更進一步瞭解本發明之特徵及技術內容,請參閱以下有關本發明之詳細說明與附圖,但是此等說明與所附圖式僅係用來說明本發明,而非對本發明的權利範圍作任何的限制。The detailed description of the present invention and the accompanying drawings are to be understood by the claims The scope is subject to any restrictions.
請參照圖1,圖1為本發明實施例之高壓輸配電網路之故障偵測、辨識與定位之方法的流程圖。高壓輸配電網路之故障偵測、辨識與定位之方法可以使用硬體電路來實現,或者使用計算機裝置配合軟體的方式來實現。Please refer to FIG. 1. FIG. 1 is a flow chart of a method for detecting, identifying, and locating a fault in a high voltage transmission and distribution network according to an embodiment of the present invention. The method of fault detection, identification and location of the high-voltage transmission and distribution network can be implemented by using a hardware circuit or by using a computer device and a software.
首先,在步驟S1中,使用負序轉換(Negative Sequence Component,NSC)偵測高壓輸配電網路是否有故障。若無故障發生,則重複進行步驟S1。若有故障發生,則進行步驟S2。然後,在步驟S2中,將故障之一個週期內的電壓信號與電流信號進行多階層小波轉換,並據此產生複數個高頻信號。First, in step S1, a Negative Sequence Component (NSC) is used to detect whether the high voltage transmission and distribution network is faulty. If no fault occurs, step S1 is repeated. If a fault occurs, step S2 is performed. Then, in step S2, the voltage signal and the current signal in one cycle of the fault are multi-level wavelet converted, and a plurality of high-frequency signals are generated accordingly.
接著,在步驟S3中,對高頻信號進行主成分分析,並據此產生至少一個信號特徵值。再來,在步驟S4中,使用支援向量機分類信號特徵值,以得到故障類型。然後,在步驟S5中,使用可適性結構類神經網路偵測故障位置。Next, in step S3, principal component analysis is performed on the high frequency signal, and at least one signal feature value is generated accordingly. Further, in step S4, the signal vector values are classified using the support vector machine to obtain the fault type. Then, in step S5, the fault location is detected using an adaptive structural neural network.
在高壓輸配電網路中,電壓與電流信號的失真(distortion)主要是歸咎於故障。故障可能產生暫態的現象,且會造成電壓信號與電流信號之三相不平衡。無論在平衡或非平衡的三相系統,負序(negative sequence)電壓信號與負序電流信號是敏感於故障與誤接所產生的信號。In high-voltage transmission and distribution networks, the distortion of voltage and current signals is mainly due to faults. A fault may cause a transient phenomenon and cause a three-phase imbalance between the voltage signal and the current signal. Negative sequence voltage signals and negative sequence current signals are sensitive to faults and misconnections, whether in balanced or unbalanced three-phase systems.
若三相電壓與三相電流表示成V abc 與I abc ,其中V abc =(V a ,V b ,V c )且I abc =(I a ,I b ,I c ),且a ~c 分別代表不同相位。三個對稱的電壓與電流信號的分量可以整理為V 012 =(V 0 ,V 1 ,V 2 )與I 012 =(I 0 ,I 1 ,I 2 ),其中0、1與2分別代表零、正與負序分量。電壓信號的負序分量(V 2, a ,V 2, b ,V 2, c )與電流信號的負序分量(I 2 , a ,I 2. b ,I 2, c )被用來作為故障指標。If the three-phase voltage and the three-phase current are expressed as V abc and I abc , where V abc = ( V a , V b , V c ) and I abc = ( I a , I b , I c ), and a to c respectively Represents different phases. The components of the three symmetrical voltage and current signals can be organized into V 012 =( V 0 , V 1 , V 2 ) and I 012 =( I 0 , I 1 , I 2 ), where 0, 1 and 2 represent zero respectively. Positive and negative sequence components. The negative sequence components of the voltage signal ( V 2, a , V 2, b , V 2, c ) and the negative sequence components of the current signal ( I 2 , a , I 2. b , I 2, c ) are used as faults index.
然而,在非平衡的三相輸配電系統中,因為不同的系統阻抗與相異的負載分佈使得負序電壓與負序電流可能會具有相當大的值。為了辨識非零的電壓負序分量V 2 與電流負序分量I 2 ,將電壓負序分量V 2 與電流負序分量I 2 對時間t 做偏微分,如下面的方程式,However, in an unbalanced three-phase power transmission and distribution system, the negative sequence voltage and the negative sequence current may have considerable values because of different system impedances and different load distributions. In order to identify the non-zero voltage negative sequence component V 2 and the current negative sequence component I 2 , the voltage negative sequence component V 2 and the current negative sequence component I 2 are differentiated from the time t , as in the following equation,
接著,為了降低錯誤警報的機率,將偏微分的電壓負序分量 V 2 / t 與電流負序分量 I 2 / t 作迴旋積分,以獲得綜合故障指示(joint fault indicator)。綜合故障指示對於故障事件的偵測較為可靠,且綜合故障指示可以如式(1)表示,Then, in order to reduce the probability of false alarms, the differential voltage component of the partial differential will be V 2 / t and current negative sequence component I 2 / t is used for the cyclotron integral to obtain a joint fault indicator. The comprehensive fault indication is more reliable for detecting fault events, and the comprehensive fault indication can be expressed as equation (1).
,其中D (t )為在時間t 量測的接點故障指示,R (τ)可以為三角面積函數R (t )=(1-100t )‧(H (t )-H (t -0.01)),其中H (t )為階梯函數。Where D ( t ) is the contact fault indication measured at time t , and R (τ) can be a triangular area function R ( t )=(1-100 t )‧( H ( t )- H ( t -0.01 )), where H ( t ) is a step function.
請參照圖2,圖2為本發明實施例之接點故障指示之參數R (t )的波形圖。R (t )用以作為綜合故障指示D (t )作迴旋積分的函數。例如,可以定義當綜合故障指示D (t )大於門限值(例如為1)時,則代表偵測到故障發生。在本實施例中,以圖2所示的R (t )的條件下,若在0.01秒的連續時間內,偏微分的電壓負序分量 V 2 / t 與電流負序分量 I 2 / t 的迴旋積分之和大於1,則可以視為偵測到故障發生。Please refer to FIG. 2. FIG. 2 is a waveform diagram of a parameter R ( t ) of a contact failure indication according to an embodiment of the present invention. R ( t ) is used as a function of the integrated fault indication D ( t ) as a cyclotron integral. For example, it can be defined that when the integrated fault indication D ( t ) is greater than a threshold (for example, 1), it indicates that a fault has been detected. In the present embodiment, under the condition of R ( t ) shown in FIG. 2, if the voltage is negatively sequenced in a continuous time of 0.01 second. V 2 / t and current negative sequence component I 2 / If the sum of the whirling integrals of t is greater than 1, it can be considered that a fault has been detected.
另外,綜合故障指示D (t )的參數R (t )可以被設計來調整偵測故障的敏感度。值得一提的是,此實施例的高壓輸配電網路之故障偵測、辨識與定位之方法使用綜合故障指示D (t )來判斷故障發生,故能避免故障之電壓或電流信號的頻率偏差與振幅變化所產生之誤判,其中頻率偏差與振幅變化導因於高壓輸配電網路中之操作參數的變化。In addition, the parameter R ( t ) of the integrated fault indication D ( t ) can be designed to adjust the sensitivity of detecting faults. It is worth mentioning that the method for detecting, identifying and locating the fault of the high-voltage transmission and distribution network of this embodiment uses the comprehensive fault indication D ( t ) to judge the occurrence of the fault, so that the frequency deviation of the voltage or current signal of the fault can be avoided. Misjudgment with amplitude variations, where frequency deviation and amplitude variation are caused by changes in operating parameters in the high voltage transmission and distribution network.
請參照圖3A至圖3E,圖3A與圖3B分別為本發明實施例之電壓信號與電壓負序分量的波形圖。圖3C與圖3D分別為本發明實施例之電流信號與電流負序分量的波形圖。圖3E為本發明實施例之接點故障指示的波形圖。於此實施例中,故障為在時間t 約為1秒時發生,且其為在距離安裝同步相量量測裝置(phasor measurement unit,PMU)的匯流排0.95p.u.(per unit)的位置所觸發三相接地的故障。Please refer to FIG. 3A to FIG. 3E . FIG. 3A and FIG. 3B are waveform diagrams of a voltage signal and a voltage negative sequence component respectively according to an embodiment of the present invention. 3C and 3D are waveform diagrams of a current signal and a current negative sequence component, respectively, according to an embodiment of the present invention. FIG. 3E is a waveform diagram of a contact fault indication according to an embodiment of the present invention. In this embodiment, the fault occurs when the time t is about 1 second, and is triggered by a position of 0.95 pu (per unit) in the busbar from which the phasor measurement unit (PMU) is installed. Three-phase ground fault.
在圖3A與圖3C中,相位a 的電壓與電流信號以虛線表示,相位b 的電壓與電流信號以實線表示,相位c 的電壓與電流信號以點線表示。由圖3A與圖3C可知,故障約在時間t 為1.00034秒時被成功偵測,電壓負序分量與電流負序分量也同時增加(請對應參照圖3B與圖3D)。由圖3B與圖3D可知,電壓負序分量與電流負序分量在故障發生時並不穩定,使得電壓負序分量與電流負序分量無法直接作為故障的指示。綜合故障指示D (t )則可作為故障分類與故障位置估計的有效的觸發。In FIGS. 3A and 3C, the voltage and current signals of phase a are indicated by broken lines, the voltage and current signals of phase b are indicated by solid lines, and the voltage and current signals of phase c are indicated by dotted lines. As can be seen from FIG. 3A and FIG. 3C, the fault is successfully detected when the time t is 1.00034 seconds, and the voltage negative sequence component is detected. Negative sequence component Also increase at the same time (please refer to FIG. 3B and FIG. 3D correspondingly). As can be seen from FIG. 3B and FIG. 3D, the voltage negative sequence component Negative sequence component It is unstable when the fault occurs, making the voltage negative sequence component Negative sequence component Cannot be directly used as an indication of a fault. The integrated fault indication D ( t ) can be used as an effective trigger for fault classification and fault location estimation.
於步驟S2中所產生的高頻信號包括第一階與第二 階小波細節信號,且甚至可以包括更高階的小波細節信號。小波轉換適用於非穩態與非週期信號在時域與頻域的分析,並據此汲取信號特徵。The high frequency signal generated in step S2 includes first order and second The order wavelet detail signal, and may even include higher order wavelet detail signals. Wavelet transform is applied to the analysis of the unsteady and aperiodic signals in the time domain and the frequency domain, and the signal characteristics are extracted accordingly.
在本實施例中,所述小波轉換可以是多階層連續小波轉換(multi-level continuous wavelet transform,ML-CWT),且多階層連續小波轉換的轉換公式可以表示為式(2),
須要注意的是,多階層連續小波轉換僅是實現小波轉換的其中一個實施例,然而,本發明並不因此限定。實際應用時,多階層連續小波轉換可以改用多階層離散小波轉換(Multi-level DWT,ML-DWT)來實施。據此,式(2)的多階層連續小波轉換可以離散化為式(3),
多階層離散小波轉換用以遞迴地分解三相電壓信號V abc 與三相電流信號I abc ,以獲得對應的小波近似係數與小波詳細係數,此過程又可稱為多重解析小波分析(multi-resolution analysis,MRA)。小波近似係數與小波詳細係數分別代表三相電壓信號V abc 與三相電流信號I abc 的低頻信號與高頻信號。The multi-level discrete wavelet transform is used to recursively decompose the three-phase voltage signal V abc and the three-phase current signal I abc to obtain a corresponding wavelet approximation coefficient and wavelet detailed coefficient. This process can also be called multi-analytic wavelet analysis (multi- Resolution analysis, MRA). The wavelet approximation coefficient and the wavelet detailed coefficient represent the low frequency signal and the high frequency signal of the three-phase voltage signal V abc and the three-phase current signal I abc , respectively.
請參照圖4,圖4為本發明實施例之多重解析小波分析的馬雷(Mallat)式離散小波分解過程的示意圖。式(3)的多階層離散小波轉換可以用馬雷演算法分解。在本實施例中,第一階與第二階小波詳細係數可以表示為(),且在步驟S3中,第一階與第二階小波細節信號會被拿來進行主成分分析。Please refer to FIG. 4. FIG. 4 is a schematic diagram of a Mallat-type discrete wavelet decomposition process of multiple analytical wavelet analysis according to an embodiment of the present invention. The multi-level discrete wavelet transform of equation (3) can be decomposed by the Marley algorithm. In this embodiment, the first-order and second-order wavelet detailed coefficients can be expressed as ( And, in step S3, the first order and second order wavelet detail signals are taken for principal component analysis.
母小波Ψ (.)的選擇與故障分類的準確性有關,在本實施例中,多貝西(Daubechies)小波的Db4小波是被選為多階層離散小波轉換(ML-DWT)的母小波。多貝西(Daubechies)小波的Db4小波已被證實是用於輸配電系統的故障分析的最有效的母小波,如參考文獻,「A.H. Osman and O.P. Malik,“Protection of parallel transmission lines using wavelet transform,”IEEE Trans. Power Deliv .,Vol. 19,no. 1,pp. 49-55,2004.」。舉例來說,若以取樣頻率3840Hz(電力系統的頻率為60Hz,每週期取樣64次)來量測電壓與電流,第一階小波詳細係數的與第二階小波詳細係數的頻率分別為960-1920Hz與480-960Hz。除此之外,透過上述兩個在高頻帶的第一階與第二階小波係數,將可以取得故障之暫態現象的資訊。And fault classification accuracy selection mother wavelet [Psi] (.) Is related, in the present embodiment, Daubechies (of Daubechies) Db4 wavelet wavelet multi-stratum has been selected as discrete wavelet transform (ML-DWT) the mother wavelet. The Db4 wavelet of the Daubechies wavelet has been proven to be the most effective mother wavelet for fault analysis of transmission and distribution systems, as in the reference, "AH Osman and OP Malik, "Protection of parallel transmission lines using wavelet transform, IEEE Trans. Power Deliv ., Vol. 19, no. 1, pp. 49-55, 2004.". For example, if the sampling frequency is 3840 Hz (the frequency of the power system is 60 Hz, 64 times per cycle is sampled) to measure the voltage and current, the frequency of the first-order wavelet detailed coefficient and the second-order wavelet detailed coefficient are respectively 960- 1920Hz and 480-960Hz. In addition, through the above two first-order and second-order wavelet coefficients in the high frequency band, information on the transient phenomenon of the fault can be obtained.
由於從實際高壓輸配電網路所量測到的信號充滿多餘的資訊,且這些多餘的資訊也存在於轉換後的小波係數中。在本實施例中,高頻信號的主成分為第一階小波細節信號與第二階小波細節信號。在故障的電壓和電流信號經過多階層離散小波轉換(ML-DWT)後,電壓和電流信號的第一階小波細節信號與第二階小波細節信號可以表示成下述的向量,Since the signal measured from the actual high-voltage transmission and distribution network is full of redundant information, and this redundant information is also present in the converted wavelet coefficients. In this embodiment, the main components of the high frequency signal are a first order wavelet detail signal and a second order wavelet detail signal. After the faulty voltage and current signals pass the multi-level discrete wavelet transform (ML-DWT), the first-order wavelet detail signal and the second-order wavelet detail signal of the voltage and current signals can be expressed as the following vectors.
其中C a 、C b 與C c 分別代表傳輸線的相位a ~c 的高頻特徵。在本實施例中,每個週期的取樣次數為64次,所以C a 、C b 與C c 分別包含64個向量(k =1~64)。C a 、C b 與C c 的平均為式(7)與式(8),Where C a , C b and C c represent the high frequency characteristics of the phases a to c of the transmission line, respectively. In this embodiment, the number of samples per cycle is 64, so C a , C b and C c respectively contain 64 vectors ( k =1 to 64). The average of C a , C b and C c is formula (7) and formula (8),
,其中x {a ,b ,c },y {d 1,d 2}。, where x { a , b , c }, y { d 1, d 2}.
接著,將C a 、C b 與C c 做置中平減(mean-centered)得到式(9),Next, C a , C b and C c are subjected to mean-centered to obtain equation (9).
其中x {a ,b ,c },代表傳輸線的相位x 的電壓與電流信號之置中平減的小波詳細係數。為了找到在的每一個欄(column)最大的可能投影向量e ,得找到一組正交向量e x (k )使得下式的值為最大,Where x { a , b , c }, A wavelet detail coefficient that represents the sum of the voltage of the phase x of the transmission line and the current signal. In order to find Each column of the largest possible projection vector e , find a set of orthogonal vectors e x ( k ) such that the value of the following formula is the largest,
而且正規化(orthonormality)限制是(l )e x (k )=δ x (l,k ),其中x {a ,b ,c },δ x (l,k )是大小為l ×k 的單位矩陣。向量e x (k )與純量λ x (k )是共變異數(covariance)Cov x =的特徵向量與特徵值。And the orthonormality limit is ( l ) e x ( k )= δ x ( l,k ), where x { a , b , c }, δ x ( l, k ) is an identity matrix of size l × k . The vector e x ( k ) and the scalar λ x ( k ) are covariances Cov x = The eigenvectors and eigenvalues.
為了獲得向量e x (k )與純量λ x (k ),假設d (k )與μ (k )是共變異數的特徵向量與特徵值,亦即 d (k )=μ k d k 。將 d (k )=μ k d k 兩邊乘上而成為,因此共變異數Cov x 的首M -1個特徵向量e x (k )與特徵值λ x (k )可以由 d (k )與μ k 得到,其中x {a ,b ,c }。In order to obtain the vector e x ( k ) and the scalar λ x ( k ), it is assumed that d ( k ) and μ ( k ) are the eigenvectors and eigenvalues of the co-variation, ie d ( k )= μ k d k . will d ( k )= μ k d k multiplied by both sides And become Therefore, the first M -1 eigenvector e x ( k ) of the covariance number Cov x and the eigenvalue λ x ( k ) can be d ( k ) and μ k are obtained, where x { a , b , c }.
需要注意的是,為了讓 d (k )與e x (k )相同, d (k )得被歸一化。對應於共異變數Cov x 的非零特徵值的特徵向量產生正規化基底,此基底形成的子空間使得大多數的傳輸線的電壓與電流的小波細節信號僅具有少量誤差。特徵向量e x (k )可以依據特徵向量的特徵值λ x (k )以遞減的方式排列。具有最大特徵值的特徵向量可用以反應小波細節信號中的最大的變化量,因為經過排列後的特徵值是以指數的方式遞減,使得大約百分之九十的變化量是在前面百分之五到百分之十的維度內。Need to pay attention to, in order to let d ( k ) is the same as e x ( k ), d ( k ) has to be normalized. Mutation corresponding to a total number of non-zero eigenvalues of Cov x normalized feature vector generation substrate, a substrate formed of this subspace such that the voltage and current of the detail wavelet most signal transmission line having only a small error. The feature vector e x ( k ) may be arranged in a decreasing manner according to the feature value λ x ( k ) of the feature vector. The feature vector with the largest eigenvalue can be used to reflect the largest amount of variation in the wavelet detail signal, since the aligned eigenvalues are exponentially decremented such that approximately ninety percent of the change is in the front percent Within five to ten percent of the dimensions.
透過計算Ωx =[υ x . 1 ,υ x . 2 ,...,υ x . M ' ],其中υ x.k ' =(k ')(k '),k '=1,2,...,M ',小波細節信號可以被投影在M '(<<64)維度的向量空間中。υ x.k ' 是在新的向量空間的小波細節信號的第k '個座標,且υ x.k ' 也是傳輸線在相位x (x {a ,b ,c })的電壓和電流信號的主要成分。簡言之,在步驟S3中,依據傳輸線的相位a 、b 、c 的高頻特徵C a 、C b 與C c ,而得到信號特徵值υ x.k ' ,其中x {a ,b ,c }。By calculating Ω x =[υ x . 1 ,υ x . 2 ,...,υ x . M ' ], where υ xk ' = ( k ') ( k '), k '=1, 2,..., M ', the wavelet detail signal can be projected in the vector space of the M '(<<64) dimension. υ xk ' is the k'th coordinate of the wavelet detail signal in the new vector space, and υ xk ' is also the transmission line at phase x ( x { a , b , c }) The main components of the voltage and current signals. In short, in step S3, the signal characteristic value υ xk ' is obtained according to the high frequency features C a , C b and C c of the phases a , b , c of the transmission line, where x { a , b , c }.
信號特徵值υ x.k ' 的首兩個維度使電壓與電流信號在傳輸線的每個相位的小波細節信號(=[υ a ,1 ,υ a ,2 ]、=[υ b . 1 ,υ b . 2 ]與=[υ c .1 ,υc .2 ])投影在一個特徵平面。在本實施例中,步驟S4所述的支援向量機用以分類此特徵空間的首兩個維度的主成分值,並據此分類出不同類型的故障。然而需要說明的是,本發明並不因此限定,支援向量機也可以分類此特徵空間的首兩個以上的維度的主成分值,並據此分類出不同類型的故障。支援向量機是一種通用學習演算法,具有單輸出且可以藉由統計理論的理論結果處理線性與非線性分離模型。因此,支援向量機(SVM)可以依據分類訓練樣本而得到的不同類型的故障,從而用以分類新產生的故障。The first two dimensions of the signal characteristic value υ xk ' make the wavelet and detail signals of the voltage and current signals at each phase of the transmission line ( =[υ a ,1 ,υ a ,2 ], =[υ b . 1 ,υ b . 2 ] =[υ c .1 , υ c .2 ]) Projected in a feature plane. In this embodiment, the support vector machine described in step S4 is used to classify the principal component values of the first two dimensions of the feature space, and classify different types of faults accordingly. It should be noted, however, that the present invention is not limited thereto. The support vector machine may also classify the principal component values of the first two or more dimensions of the feature space, and classify different types of faults accordingly. The support vector machine is a general learning algorithm with a single output and can handle linear and nonlinear separation models by theoretical results of statistical theory. Therefore, the support vector machine (SVM) can classify newly generated faults based on different types of faults obtained by classifying the training samples.
請參照圖5,圖5為本發明實施例之故障類型的示意圖。在圖5中,電阻R f 代表故障電阻,電阻R g 代表接地電阻,電阻HIF 代表高阻抗故障電阻。在本實施例中的故障類型可以分為以下六類:(a)單相接地的故障,“相位a 接地”、“相位b 接地”與“相位c 接地”;(b)雙相位接地的故障,“相位及b 接地”、“相位b 及c 接地”與“相位a 及c 接地”;(c)三相接地的故障,“相位a 、b 、c 接地”;(d)雙相位間的故障,“相位a 及相位b ”、“相位b 及相位c ”與“相位a 及相位c ”;(e)三相故障,“相位a 、b 、c 間的故障”;以及(f)高阻抗故障(high-impedance fault,HIF)。Please refer to FIG. 5. FIG. 5 is a schematic diagram of a fault type according to an embodiment of the present invention. In FIG. 5, the resistance R f represents a fault resistance, the resistance R g represents a ground resistance, and the resistance HIF represents a high impedance fault resistance. The types of faults in this embodiment can be classified into the following six categories: (a) single-phase ground fault, "phase a ground", "phase b ground" and "phase c ground"; (b) double phase ground fault , "phase and b ground", "phase b and c ground" and "phase a and c ground"; (c) three-phase ground fault, "phase a , b , c ground"; (d) between two phases Faults, "phase a and phase b ", "phase b and phase c " and "phase a and phase c "; (e) three-phase fault, "fault between phases a , b , c "; and (f) High-impedance fault (HIF).
為了正確的分類以上故障類型,一組六個支援向量機(SVM#1~SVM#6)以及六個支援向量機(SVM#1~SVM#6)所對應的六個可適性結構類神經網路(Adaptive Structure Neural Networks,ASNN)被整合用來分類故障與定位。六個支援向量機(SVM#1~SVM#6)分別對應於前述的六個故障類型(a)~(f)。支援向量機依據訓練樣本產生特徵空間與決定邊界,且依據故障所對應的信號特徵值在特徵空間中的位置,配合決定邊界,以分類信號特徵值。藉此,可以使用支援向量機對故障進行分類。另外,支援向量機依據訓練樣本而產生特徵空間與決定邊界的方式如下述。In order to correctly classify the above fault types, a set of six support vector machines (SVM#1 to SVM#6) and six support vector machines (SVM#1 to SVM#6) correspond to six adaptive structural neural networks. Adaptive Structure Neural Networks (ASNN) are integrated to classify faults and locations. The six support vector machines (SVM#1 to SVM#6) correspond to the aforementioned six failure types (a) to (f), respectively. The support vector machine generates the feature space and the decision boundary according to the training sample, and according to the position of the signal feature value corresponding to the fault in the feature space, the decision boundary is matched to classify the signal feature value. Thereby, the fault can be classified using the support vector machine. In addition, the manner in which the support vector machine generates the feature space and determines the boundary based on the training sample is as follows.
給與一組已知知識庫(prior knowledge vector)=[υ x ,y x ],其中x {a ,b ,c },且y x 是關聯於υ x 的標籤(label),也就是已知故障類別。支援向量機(SVM)的分辨函數可表示成,其中i =1,2,...,N ,N 表示超平面(hyperplane)的數目,是權重向量,b i 是偏移項。每個超平面將特徵空間分成兩個部份,而回歸函數f i ()所具有的正負符號用以表示特徵點是在超平面的其中一側。Give a set of known knowledge bases =[υ x , y x ], where x { a , b , c }, and y x is the label associated with υ x , which is the known fault category. The resolution function of the support vector machine (SVM) can be expressed as , where i =1, 2,..., N , N represents the number of hyperplanes, Is the weight vector, and b i is the offset term. Each hyperplane divides the feature space into two parts, and the regression function f i ( The positive and negative signs are used to indicate that the feature point is on one side of the hyperplane.
非線性的分類機(classifier)可以用於分類故障信號的非線性的暫態現象,但同時得使用非線性的轉移函數φ以從輸入空間至特徵空間來繪製特徵點。分辨函數可以是,其中。在此情況下,分類機的決定邊界可以寫成。這個轉換給決定邊界增加了彈性,使得非線性決定邊界可以被達成。A nonlinear classifier can be used to classify the nonlinear transient phenomena of the fault signal, but at the same time a nonlinear transfer function φ is used to draw the feature points from the input space to the feature space. The resolution function can be ,among them . In this case, the decision boundary of the sorter can be written as . This transformation adds flexibility to the decision boundary so that nonlinearly determined boundaries can be achieved.
再者,假設權重向量可以用訓練樣本()的線性組合來表示成。然後,決定邊界可以表示成。在特徵空間(在本實施例中,特徵空間為特徵平面)中,決定邊界可以寫成。以變數α x 為項的表示法是被稱為決定邊界的雙表徵(dual representation)。Furthermore, assume that the weight vector can use training samples ( Linear combination to represent . Then, the decision boundary can be expressed as . In the feature space (in this embodiment, the feature space is a feature plane), the decision boundary can be written as . A notation in which the variable α x is a term is a dual representation called a decision boundary.
與b i 的值可被用來求解懲罰風險函數(regularized risk function)(式(10)), The value of b i can be used to solve the regularized risk function (Eq. (10)),
其中懲罰風險函數被限制於“與ξ x 0”的條件中。ξ x 代表惰變數(slack variable),惰變數允許特徵點在超平面的界限(0ξ x 1,已知為界限誤差),或者特徵點被誤分類(ξ x 1),而C為懲罰參數。Where the penalty risk function is limited to " With ξ x In the condition of 0". ξ x represents the slack variable, and the idling allows the feature point to be at the boundary of the hyperplane (0 ξ x 1, known as the limit error), or the feature points are misclassified (ξ x 1), and C is the penalty parameter.
藉由拉格朗日乘數(Lagrange multipliers)的方法可以從式(10)得到以變數α x 為項的表示法的式(11),By the method of Lagrange multipliers, the equation (11) with the representation of the variable α x can be obtained from the equation (10).
其中式(11)被限制於“”的條件下。Where equation (11) is limited to " "Under conditions.
請參照圖6A與圖6B,圖6A與圖6B分別為本發明實施例之相位接地故障與高阻抗故障的特徵平面圖。圖6A與6B中的虛線是分別對應於相位a 、b 、c 的決定邊界。依據訓練樣本的信號特徵值,每一個訓練樣本在圖6A與圖6B上產生一個點。每個點旁邊所標示的位置是透過可適性類神經網路所估計的傳輸線上的故障位置。0 p.u.與1 p.u.代表從匯流排至故障發生處的距離。例如,0 p.u.為0公里,1 p.u.為1000公里。Please refer to FIG. 6A and FIG. 6B. FIG. 6A and FIG. 6B are respectively a plan view of a phase ground fault and a high impedance fault according to an embodiment of the present invention. The broken lines in Figs. 6A and 6B are decision boundaries corresponding to the phases a , b , and c , respectively. Each training sample produces a point on Figures 6A and 6B based on the signal characteristic values of the training samples. The location indicated next to each point is the location of the fault on the transmission line estimated by the adaptive neural network. 0 pu and 1 pu represent the distance from the busbar to the point where the fault occurred. For example, 0 pu is 0 km and 1 pu is 1000 km.
請參照圖7,圖7為本發明實施例使用可適性結構類神經網路偵測故障事件的位置的子步驟的流程圖。圖1中的步驟S5更可以包括下述的子步驟。首先,在子步驟S51中,進行可適性結構類神經網路神經元權重的建構演算法,以獲得可適性結構類神經網路。然後,在子步驟S52中,輸入故障之電壓與電流信號的主成分至可適性結構類神經網路。接著,在子步驟S53中,使用可適性結構類神經網路對故障之電壓與電流信號的主成分進行分析,以獲得故障位置。Please refer to FIG. 7. FIG. 7 is a flow chart showing the sub-steps of detecting the location of a fault event using an adaptive structural neural network according to an embodiment of the present invention. Step S5 in Fig. 1 may further include the sub-steps described below. First, in sub-step S51, a construction algorithm for adaptive structural neural network neuron weights is performed to obtain an adaptive structural neural network. Then, in sub-step S52, the principal components of the fault voltage and current signals are input to the adaptive structural neural network. Next, in sub-step S53, the principal components of the fault voltage and current signals are analyzed using an adaptive structural neural network to obtain the fault location.
請同時參照圖與圖7與圖8,圖8為本發明實施例之可適性結構類神經網路的架構圖。於步驟S51中所獲得之可適性結構類神經網路(ASNN)如圖8所示,此可適性結構類神經網路包括輸入層、隱藏層與輸出層。Please refer to the same figure and FIG. 7 and FIG. 8. FIG. 8 is a structural diagram of an adaptive structural neural network according to an embodiment of the present invention. The adaptive structural neural network (ASNN) obtained in step S51 is as shown in FIG. 8. The adaptive structural neural network includes an input layer, a hidden layer and an output layer.
輸入層之資料輸入節點X 1 、X 2 …X p 會在故障發生時,輸入小波轉換後的主成分。在本實施例中,主成分為第一階小波細節信號與第二階小波細節信號。隱藏層中的神經元節點N 1 、N 2 …N p 用以連結資料輸入節點或自身外之神經元節點N 1 、N 2 …N n ,以依照輸入節點-神經元節點權重W X (1、2… p ),(1、2… n ) 、神經元節點-神經元節點權重W (1、2… n ).(1、2… n ) 與神經元節點-輸出節點權重W Y (1、2… n ),(1、2… q ) 估測故障位置。上述的輸入節點-神經元節點權重W X (1、2… p ),(1、2… n ) 與神經元節點-輸出節點權重W Y (1、2… n ),(1、2… q ) 是經由可適性結構類神經網路神經元權重建構演算法所求得,利用此演算法可以使可適性結構類神經網路之架構能符合環境需求。The data input nodes X 1 , X 2 ... X p of the input layer input the wavelet transformed principal component when a fault occurs. In this embodiment, the principal component is a first order wavelet detail signal and a second order wavelet detail signal. The neuron nodes N 1 , N 2 ... N p in the hidden layer are used to link the data input nodes or the neuron nodes N 1 , N 2 ... N n outside the node to the input node-neuron node weight W X (1 , 2... p ), (1, 2... n ) , neuron node-neuron node weight W (1, 2... n ). (1, 2... n ) and neuron node-output node weight W Y (1 , 2... n ), (1, 2... q ) Estimate the fault location. The above input node-neuron node weights W X (1, 2... p ), (1, 2... n ) and neuron node-output node weights W Y (1, 2... n ), (1, 2... q It is obtained through the adaptive structure-based neural network neuron weight reconstruction algorithm, and the algorithm can make the architecture of the adaptive structure-like neural network meet the environmental requirements.
請參照圖9,圖9為本發明實施例之可適性結構類神經網路神經元權重的建構演算法的流程圖。獲得可適性結構類神經網路的步驟(步驟S51)包括以下子步驟。首先,在步驟S511中,隨機產生隱藏層之神經元個數、權重以及複數個神經元連線。然後,在步驟S512中,計算輸出之誤差率。再來,在步驟S513中,判斷輸出之誤差率是否可被接受。若輸出之誤差率可以被接收,則完成可適性結構類神經網路神經元權重的建構,並進行後續的步驟S52。若輸出之誤差率不可以被接受,則進行步驟S514。Please refer to FIG. 9. FIG. 9 is a flowchart of a constructing algorithm for the weight of the neural network of the adaptive structure-like neural network according to an embodiment of the present invention. The step of obtaining an adaptive structural neural network (step S51) includes the following sub-steps. First, in step S511, the number of neurons of the hidden layer, the weight, and a plurality of neuron connections are randomly generated. Then, in step S512, the error rate of the output is calculated. Further, in step S513, it is judged whether or not the error rate of the output is acceptable. If the error rate of the output can be received, the construction of the adaptive structure-like neural network neuron weights is completed, and the subsequent step S52 is performed. If the error rate of the output is not acceptable, step S514 is performed.
在步驟S514中,調整複數個神經元節點-輸出節點權重。然後,在步驟S515中,調整複數個輸入節點-神經元節點權重。然後,在步驟S516中,調整複數個神經元節點-神經元節點權重。接著,在步驟S517中,建立至少一個新神經元或刪除至少一個神經元。再來,在步驟S518中,刪除複數個神經元間之連線中的至少一個或建立至少一個神經元連線。然後,重新進行步驟S512,即計算輸出之誤差率。如此重複步驟S512~S518,直至在步驟S513中判斷輸出之誤差率可以被接受,則完成可適性結構類神經網路神經元權重的建構,並進行後續的步驟S52。In step S514, a plurality of neuron node-output node weights are adjusted. Then, in step S515, a plurality of input node-neuron node weights are adjusted. Then, in step S516, a plurality of neuron node-neuron node weights are adjusted. Next, in step S517, at least one new neuron is established or at least one neuron is deleted. Then, in step S518, at least one of the connections between the plurality of neurons is deleted or at least one neuron connection is established. Then, step S512 is performed again, that is, the error rate of the output is calculated. Steps S512 to S518 are repeated in this manner until it is judged in step S513 that the error rate of the output can be accepted, and the construction of the adaptive structure-like neural network neuron weight is completed, and the subsequent step S52 is performed.
須要注意的是,可適性結構類神經網路神經元權重的建構即是對可適性結構類神經網路的訓練,而不需要重新設定可適性結構類神經網路。因為電力公司可能無法負擔關閉電力網格的偵測系統太久的風險及成本,故此種採用可適性結構類神經網路來偵測故障位置的方式更可以滿足目前電力公司的需求。It should be noted that the construction of adaptive structural neural network neuron weights is the training of adaptive structural neural networks without the need to reset the adaptive structural neural network. Because power companies may not be able to afford the risk and cost of shutting down the power grid detection system for too long, such an adaptive structural neural network to detect fault locations can meet the needs of current power companies.
復參照圖1,透過上述的步驟S11~S15,故障類型與故障位置可以被獲得,藉此可以進行防護決策。防護決策可以包括採取行動將故障排除,以減低故障所造成的損害。防護決策也可以包括記錄分析結果以作進一步的分析。總而言之,電力公司可以自訂對應各故障類型的防護決策。Referring back to FIG. 1, through the above steps S11 to S15, the fault type and the fault location can be obtained, whereby the protection decision can be made. Protection decisions can include taking action to troubleshoot to reduce the damage caused by the failure. Protection decisions may also include recording analysis results for further analysis. All in all, the power company can customize the protection decisions for each type of failure.
請參照圖10,圖10是本發明實施例之高壓輸配電網路之故障偵測、辨識與定位之系統的架構圖。高壓輸配電網路具有三個相位的輸配線與用以電性連接負載的匯流排。為了方便說明,本實施例之高壓輸配電網路包括三個相位(S a 、S b 、S c )的輸配線與匯流排Bus#1~5,且每個匯流排可以電性連接至負載。如圖10所示,匯流排Bus#1~4分別電性連接至負載。Please refer to FIG. 10. FIG. 10 is a structural diagram of a system for detecting, identifying, and locating faults in a high voltage transmission and distribution network according to an embodiment of the present invention. The high-voltage power transmission and distribution network has three phases of transmission lines and bus bars for electrically connecting the loads. For convenience of description, the high-voltage transmission and distribution network of the embodiment includes three phases ( S a , S b , S c ) of the transmission line and the bus bars Bus #1 to 5, and each bus bar can be electrically connected to the load. . As shown in FIG. 10, the bus bars Bus#1~4 are electrically connected to the load respectively.
高壓輸配電網路之故障偵測、辨識與定位之系統1包括量測單元1a、1c、1e、故障診斷裝置10與波形記錄 器11。量測單元1a、1c、1e分別電性連接至匯流排Bus#1、Bus#3、Bus#5。故障診斷裝置10透過廣域網路(Wide Area Network,WAN)接收量測單元1a、1c、1e所量測到的電壓與電流信號。波形記錄器11透過廣域網路(WAN)連接至故障診斷裝置10。The system 1 for detecting, identifying and locating faults of a high-voltage transmission and distribution network includes measurement units 1a, 1c, 1e, fault diagnosis device 10 and waveform recording Device 11. The measuring units 1a, 1c, and 1e are electrically connected to the bus bars Bus#1, Bus#3, and Bus#5, respectively. The fault diagnosis apparatus 10 receives the voltage and current signals measured by the measurement units 1a, 1c, and 1e through a Wide Area Network (WAN). The waveform recorder 11 is connected to the fault diagnostic device 10 through a wide area network (WAN).
量測單元1a、1c、1e用以獲得三個相位(S a 、S b 、S c )的電壓與電流信號。故障診斷裝置10接收量測單元1a、1c、1e所獲得的三個相位的電壓與電流信號,且故障診斷裝置10偵測、分類與定位故障。波形記錄器11用以記錄故障發生時的三個相位的電壓與電流信號。The measuring units 1a, 1c, 1e are used to obtain voltage and current signals of three phases ( S a , S b , S c ). The fault diagnosis apparatus 10 receives the voltage and current signals of the three phases obtained by the measurement units 1a, 1c, and 1e, and the fault diagnosis apparatus 10 detects, classifies, and locates the fault. The waveform recorder 11 is used to record voltage and current signals of three phases when a fault occurs.
請參照圖11,圖11是本發明實施例之故障診斷裝置的方塊圖。故障診斷裝置10包括信號量測模組12、故障診斷模組13與圖像顯示模組14。信號量測模組12包括類比/數位轉換器(AD converter)15~17。故障診斷模組13包括故障偵測模組(未繪示)、故障分類模組(未繪示)與故障定位模組(未繪示)。在本實施例中,故障診斷模組13由可程式化邏輯閘陣列(Field-Programmable Gate Array,FPGA)實施,使得故障偵測模組、故障分類模組與故障定位模組的功能整合於可程式化邏輯閘陣列(FPGA),但本發明並不因此限定。故障偵測、故障分類與故障定位的功能也可以用其他硬體電路或者軟體程式來實施。Please refer to FIG. 11. FIG. 11 is a block diagram of a fault diagnosis apparatus according to an embodiment of the present invention. The fault diagnosis device 10 includes a signal measurement module 12, a fault diagnosis module 13, and an image display module 14. The signal measurement module 12 includes analog to digital converters (AD converters) 15-17. The fault diagnosis module 13 includes a fault detection module (not shown), a fault classification module (not shown), and a fault location module (not shown). In this embodiment, the fault diagnosis module 13 is implemented by a Field-Programmable Gate Array (FPGA), so that the functions of the fault detection module, the fault classification module, and the fault location module are integrated. A programmatic logic gate array (FPGA), but the invention is not so limited. The functions of fault detection, fault classification and fault location can also be implemented by other hardware circuits or software programs.
信號量測模組12的輸入端接收高壓輸配電網的三個相位的電壓與電流信號。信號量測模組12的輸出端電性連接故障診斷模組13的輸入端。故障診斷模組13的輸出端電性連接至圖像顯示模組14。The input of the signal measurement module 12 receives the voltage and current signals of the three phases of the high voltage transmission and distribution network. The output end of the signal measurement module 12 is electrically connected to the input end of the fault diagnosis module 13. The output end of the fault diagnosis module 13 is electrically connected to the image display module 14.
信號量測模組12之類比/數位轉換器(AD converter)15~17用以對電壓與電流信號進行類比/數位轉換。在本實施例中類比/數位轉換器(AD converter)15~17為14位元的類比/數位轉換器,但本發明並不因此限定。另外,為了使類比/數位轉換器(AD converter)15~17所輸出的14位元的數位信號(電壓信號或電流信號)便於處理,此14位元的數位信號可以轉換成浮點數位信號。所以信號量測模組12更可包括另一數位信號轉換器,用以將數位信號轉換成浮點數位信號。例如:將此14位元的數位信號轉換為16位元的半精度格式。The analog/digital converters 15-17 of the signal measurement module 12 are used for analog/digital conversion of voltage and current signals. In the present embodiment, the analog-to-digital converters (AD converters) 15 to 17 are 14-bit analog/digital converters, but the present invention is not limited thereto. In addition, in order to facilitate the processing of the 14-bit digital signal (voltage signal or current signal) output by the analog/digital converters (15 to 17), the 14-bit digital signal can be converted into a floating-point digital signal. . Therefore, the signal measurement module 12 can further include another digital signal converter for converting the digital signal into a floating point digital signal. For example: Convert this 14-bit digital signal to a 16-bit half-precision format.
故障診斷模組13接收數位化的電壓信號與電流信號,以進行負序轉換以判斷故障是否發生,接著進行多階層小波轉換、主成分分析與使用支援向量機對故障進行分類,最後再使用可適性結構類神經網路對故障進行定位。The fault diagnosis module 13 receives the digitized voltage signal and the current signal for negative sequence conversion to determine whether the fault occurs, and then performs multi-level wavelet transform, principal component analysis, and uses the support vector machine to classify the fault, and finally uses the fault. The adaptive structural neural network locates the fault.
負序轉換可以透過加法器與乘法器的組合以計算出如前一實施例的電壓負序分量 V 2 / t 與電流負序分量 I 2 / t ,進而產生綜合故障指示D (t )。當接點故障指示D (t )大於門限值時,則代表故障被偵測到,藉此故障診斷模組13進行後續的小波轉換。有關於綜合故障指示D (t )的敘述可以參照前一實施例的說明。The negative sequence conversion can be performed by a combination of an adder and a multiplier to calculate a voltage negative sequence component as in the previous embodiment. V 2 / t and current negative sequence component I 2 / t , which in turn produces a comprehensive fault indication D ( t ). When the contact failure indication D ( t ) is greater than the threshold value, the fault is detected, whereby the fault diagnosis module 13 performs subsequent wavelet conversion. For the description of the integrated fault indication D ( t ), reference may be made to the description of the previous embodiment.
故障診斷模組13利用多階層離散小波轉換(ML-DWT)進行小波轉換。多階層離散小波轉換(ML-DWT)可以用馬雷(Mallat)演算法計算。基於馬雷(Mallat)演算法第j 階層的小波係數可以用(j -1)階層的小波係數表示,如下列方程式:The fault diagnosis module 13 performs wavelet transform using multi-level discrete wavelet transform (ML-DWT). Multi-level discrete wavelet transform (ML-DWT) can be calculated using the Mallat algorithm. The wavelet coefficients based on the jth level of the Mallat algorithm can be expressed by the wavelet coefficients of the ( j -1) hierarchy, such as the following equation:
其中H (z )與G (z )為小波低通濾波係數與小波高通濾波係數,為電壓在j 階層的小波近似信號與小波細節信號。為電流在j 階層的小波近似信號與小波細節信號。為了提高運算的速度,可以使用串列輸入並列輸出(serial-in parallel-out,SIPO)的架構實施上述的運算。Where H ( z ) and G ( z ) are wavelet low-pass filter coefficients and wavelet high-pass filter coefficients. Wavelet approximation signals and wavelet detail signals for voltages in the j- level. Wavelet approximation signal and wavelet detail signal for current in the j- level. In order to increase the speed of the operation, the above operation can be implemented using a serial-in parallel-out (SIPO) architecture.
完成多階層小波分析後,故障診斷模組13進行主成分分析,用以減少多階層離散小波轉換(ML-DWT)所產生的小波係數的維度。主成分分析的方法可以是前一實施例所敘述的方法。主成分分析的計算最密集的部份是計算特徵值與排序。After completing the multi-level wavelet analysis, the fault diagnosis module 13 performs principal component analysis to reduce the dimension of the wavelet coefficients generated by the multi-level discrete wavelet transform (ML-DWT). The method of principal component analysis may be the method described in the previous embodiment. The most computationally intensive part of principal component analysis is the calculation of eigenvalues and ordering.
以故障診斷模組13由可程式化邏輯閘陣列(FPGA)實施為例,實施主成分分析可以使用多階段的管線架構,而且使用大量數目的多階段的管線可以減少計算的等待時間。例如,使用具有48982個邏輯元件的Stratix-III FPGA,其中Stratix-III FPGA之最大的運算時脈可以是76MHz,且此運算時脈受益於大量的管線化電路設計(pipe stage)。Taking the fault diagnosis module 13 as an example of a programmable logic gate array (FPGA) implementation, the implementation of principal component analysis can use a multi-stage pipeline architecture, and the use of a large number of multi-stage pipelines can reduce the computational latency. For example, a Stratix-III FPGA with 48982 logic elements is used, where the maximum operational clock of a Stratix-III FPGA can be 76 MHz, and this operational clock benefits from a large number of pipeline stages.
再來,故障診斷模組13處理支援向量機。支援向量機包括兩個階段,學習階段與故障分類階段。在學習階段,如同前一實施例所提及的,在特徵空間中,決定邊界可以寫成。變數α x 與b i 的值可以事先計算獲得,且儲存於故障診斷模組12的儲存單元。在故障分類階段,代表主成分值的首兩個座標的小波細節信號()描述故障信號的特徵,藉此故障可依據決定邊界而被分類。Then, the fault diagnosis module 13 processes the support vector machine. The support vector machine consists of two phases, a learning phase and a fault classification phase. In the learning phase, as mentioned in the previous embodiment, in the feature space, the decision boundary can be written as . The values of the variables α x and b i can be calculated in advance and stored in the storage unit of the fault diagnosis module 12. In the fault classification phase, the wavelet detail signal representing the first two coordinates of the principal component value ( Describe the characteristics of the fault signal, whereby the fault can be classified according to the decision boundary.
然後,故障診斷模組13處理可適性結構類神經網路。可適性結構類神經網路的處理步驟可以參照前一實施例的說明,在此不再贅述。須要注意的是,在實施時由於硬體資源的限制,神經元的數目與每一個神經元的參數(包括權重、神經元連線數等)也是有限的。神經元個數、權重以及複數個神經元連線可以儲存於可程式化邏輯閘陣列(FPGA)的查找表(lookup table)。在計算過程中所產生的新的神經元個數、權重以及複數個神經元連線可以儲存於可程式化邏輯閘陣列(FPGA)的暫存器。The fault diagnostic module 13 then processes the adaptive structural neural network. For the processing steps of the adaptive structural neural network, reference may be made to the description of the previous embodiment, and details are not described herein again. It should be noted that the number of neurons and the parameters of each neuron (including weights, number of connected neurons, etc.) are also limited due to limitations of hardware resources during implementation. The number of neurons, the weights, and the plurality of neuron connections can be stored in a lookup table of a programmable logic gate array (FPGA). The number of new neurons, weights, and complex neuron connections generated during the calculation can be stored in a register of a programmable logic gate array (FPGA).
圖像顯示模組14用以顯示故障偵測、辨識與定位的結果。圖像顯示模組14也可以顯示故障的波形。圖像顯示模組14可以是顯示器,例如是液晶螢幕。另外,圖像顯示模組14也可以是觸控螢幕,使得處理故障事件的操作員或工程師可以透過觸控的方式點選所要觀看的故障資訊。The image display module 14 is configured to display the results of fault detection, identification, and positioning. The image display module 14 can also display the waveform of the fault. The image display module 14 can be a display, such as a liquid crystal screen. In addition, the image display module 14 can also be a touch screen, so that an operator or an engineer who handles a fault event can select the fault information to be viewed through touch.
根據本發明實施例,上述的高壓輸配電網路之故障偵測、辨識與定位之方法及其系統可偵測故障、判定故障類別與估計故障位置。另外,所述高壓輸配電網路之故障偵測、辨識與定位之方法及其系統自故障發生後,其完成故障偵測、辨識與定位的時間遠短於臨界故障清除時間。According to an embodiment of the invention, the above method and system for detecting, identifying and locating a high voltage transmission and distribution network can detect a fault, determine a fault category and estimate a fault location. In addition, the method for detecting, identifying and locating the fault of the high-voltage transmission and distribution network and the system after the fault occurs, the time for detecting, identifying and locating the fault is much shorter than the critical fault clearing time.
以上所述僅為本發明之實施例,其並非用以侷限本發明之專利範圍。The above description is only an embodiment of the present invention, and is not intended to limit the scope of the invention.
S1~S5、S51~S53、S511~S518...步驟流程S1~S5, S51~S53, S511~S518. . . Step flow
1‧‧‧高壓輸配電網路之故障偵測、辨識與定位之系統1‧‧‧System for fault detection, identification and location of high voltage transmission and distribution networks
10‧‧‧故障診斷裝置10‧‧‧Diagnostic device
11‧‧‧波形記錄器11‧‧‧ Waveform recorder
1a、1c、1e‧‧‧量測單元1a, 1c, 1e‧‧‧ measurement unit
12‧‧‧信號量測模組12‧‧‧Signal Measurement Module
13‧‧‧故障診斷模組13‧‧‧Troubleshooting Module
14‧‧‧圖像顯示模組14‧‧‧Image display module
15~17‧‧‧類比/數位轉換器15~17‧‧‧ Analog/Digital Converter
圖1為本發明實施例之高壓輸配電網路之故障偵測、辨識與定位之方法的流程圖。1 is a flow chart of a method for detecting, identifying, and locating a fault in a high voltage transmission and distribution network according to an embodiment of the present invention.
圖2為本發明實施例之接點故障指示之參數R (t )的波形圖。2 is a waveform diagram of a parameter R ( t ) of a contact failure indication according to an embodiment of the present invention.
圖3A與圖3B分別為本發明實施例之電壓信號與電壓負序分量的波形圖。3A and 3B are waveform diagrams of a voltage signal and a voltage negative sequence component, respectively, according to an embodiment of the present invention.
圖3C與圖3D分別為本發明實施例之電流信號與電流負序分量的波形圖。3C and 3D are waveform diagrams of a current signal and a current negative sequence component, respectively, according to an embodiment of the present invention.
圖3E為本發明實施例之接點故障指示的波形圖。FIG. 3E is a waveform diagram of a contact fault indication according to an embodiment of the present invention.
圖4為本發明實施例之多重解析小波分析的馬雷(Mallat)式離散小波分解過程的示意圖。4 is a schematic diagram of a Mallat-type discrete wavelet decomposition process of multiple analytical wavelet analysis according to an embodiment of the present invention.
圖5為本發明實施例之故障類型的示意圖。FIG. 5 is a schematic diagram of a fault type according to an embodiment of the present invention.
圖6A與圖6B分別為本發明實施例之相位接地故障與高阻抗故障的特徵平面圖。6A and 6B are respectively a plan view of a phase ground fault and a high impedance fault according to an embodiment of the present invention.
圖7為本發明實施例使用可適性結構類神經網路偵測故障事件的位置的示意圖。FIG. 7 is a schematic diagram of a location of detecting a fault event using an adaptive structural neural network according to an embodiment of the present invention.
圖8為本發明實施例之可適結構性類神經網路的架構圖。FIG. 8 is a structural diagram of an adaptive structural neural network according to an embodiment of the present invention.
圖9為本發明實施例之可適性結構類神經網路神經元權重的建構演算法的流程圖。FIG. 9 is a flowchart of a construction algorithm for the weight of a neural network element of an adaptive structure based on an embodiment of the present invention.
圖10為本發明實施例之高壓輸配電網路之故障偵測、辨識與定位之系統的架構圖。FIG. 10 is a structural diagram of a system for detecting, identifying, and locating a fault of a high voltage transmission and distribution network according to an embodiment of the present invention.
圖11為本發明實施例之故障診斷模組的方塊圖。FIG. 11 is a block diagram of a fault diagnosis module according to an embodiment of the present invention.
S1~S5...步驟流程S1 ~ S5. . . Step flow
Claims (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100123915A TWI425226B (en) | 2011-07-06 | 2011-07-06 | Method and system for fault detection, identification and location in high-voltage power transmission networks |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW100123915A TWI425226B (en) | 2011-07-06 | 2011-07-06 | Method and system for fault detection, identification and location in high-voltage power transmission networks |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201303319A TW201303319A (en) | 2013-01-16 |
TWI425226B true TWI425226B (en) | 2014-02-01 |
Family
ID=48138021
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW100123915A TWI425226B (en) | 2011-07-06 | 2011-07-06 | Method and system for fault detection, identification and location in high-voltage power transmission networks |
Country Status (1)
Country | Link |
---|---|
TW (1) | TWI425226B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104699073A (en) * | 2015-02-05 | 2015-06-10 | 哈尔滨工程大学 | Autonomous type underwater robot fault identification method based on wavelet fractal |
TWI554765B (en) * | 2014-03-19 | 2016-10-21 | 國立高雄師範大學 | Method for power line outage identification |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454549B (en) * | 2013-09-18 | 2016-03-09 | 福州大学 | A kind of looped network distribution system short trouble early detection, identification and routing method |
CN104655980B (en) * | 2014-11-05 | 2019-04-09 | 云南电网公司楚雄供电局 | A fast identification method for short-circuit faults of adjacent lines with non-local voltage levels based on db4 wavelet decomposition |
CN105021925B (en) * | 2015-07-15 | 2018-06-29 | 深圳市双合电气股份有限公司 | Power monitoring system and method based on cement industry |
CN105137354B (en) * | 2015-08-27 | 2018-04-06 | 江苏科技大学 | One kind is based on neutral net electrical fault detection method |
CN105425150B (en) * | 2015-11-09 | 2018-08-31 | 江苏科技大学 | A kind of Method of Motor Fault Diagnosis based on RBF and PCA-SVDD |
CN105974265B (en) * | 2016-04-29 | 2018-11-27 | 北京四方继保自动化股份有限公司 | A kind of electric network fault cause diagnosis method based on svm classifier technology |
CN106950446A (en) * | 2017-03-22 | 2017-07-14 | 国网上海市电力公司 | The abnormal quick detection of power network based on PCA and localization method |
EP3921658B1 (en) * | 2019-02-06 | 2023-04-05 | Lisa Dräxlmaier GmbH | Method and testing device |
CN109613402B (en) * | 2019-02-14 | 2020-11-27 | 福州大学 | Detection method of high-resistance grounding fault in distribution network based on wavelet transform and neural network |
CN110988590B (en) * | 2019-11-25 | 2022-06-14 | 云南电网有限责任公司临沧供电局 | PCA-SVM model-based distribution network line selection method and system |
CN111537830A (en) * | 2020-03-26 | 2020-08-14 | 南方电网科学研究院有限责任公司 | Power distribution network fault diagnosis method based on cloud edge architecture and wavelet neural network |
CN111458602B (en) * | 2020-05-27 | 2022-05-10 | 南京信息工程大学 | Single-phase earth fault detection method based on symbiotic multi-functional calculation |
CN111650472B (en) * | 2020-06-05 | 2023-04-18 | 众诚开源电气科技(成都)有限公司 | Method for positioning voltage sag source |
CN112462193B (en) * | 2020-11-05 | 2023-10-24 | 重庆邮电大学 | Automatic reclosing judgment method for power distribution network based on real-time fault filtering data |
CN112763853B (en) * | 2020-12-29 | 2022-03-15 | 福州大学 | System for detecting and positioning short-circuit fault of alternating-current micro-grid in grid-connected mode and working method thereof |
CN113466624A (en) * | 2021-06-30 | 2021-10-01 | 云南电网有限责任公司电力科学研究院 | Method and system for detecting fault area of multi-terminal hybrid direct-current transmission line |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW437138B (en) * | 1996-12-31 | 2001-05-28 | Abb Transmit Oy | System for locating faults and estimating fault resistance in distribution networks with tapped loads |
TW200914846A (en) * | 2007-09-20 | 2009-04-01 | Univ Far East | Detection method of electronic signal using wavelet transformation |
US20110074437A1 (en) * | 2009-09-30 | 2011-03-31 | Korea Electric Power Corporation | System and method for detecting a location of fault in a cable |
-
2011
- 2011-07-06 TW TW100123915A patent/TWI425226B/en active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW437138B (en) * | 1996-12-31 | 2001-05-28 | Abb Transmit Oy | System for locating faults and estimating fault resistance in distribution networks with tapped loads |
TW200914846A (en) * | 2007-09-20 | 2009-04-01 | Univ Far East | Detection method of electronic signal using wavelet transformation |
US20110074437A1 (en) * | 2009-09-30 | 2011-03-31 | Korea Electric Power Corporation | System and method for detecting a location of fault in a cable |
Non-Patent Citations (1)
Title |
---|
"基于小波变换和支持向量机的电力电子故障诊断", 胡双俊, 2008.04 "應用人工智慧於電力系統諧波源與位置偵測", 杜耿邦, 2003.06。 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI554765B (en) * | 2014-03-19 | 2016-10-21 | 國立高雄師範大學 | Method for power line outage identification |
CN104699073A (en) * | 2015-02-05 | 2015-06-10 | 哈尔滨工程大学 | Autonomous type underwater robot fault identification method based on wavelet fractal |
CN104699073B (en) * | 2015-02-05 | 2017-06-20 | 哈尔滨工程大学 | Autonomous type underwater robot fault identification method based on small wave fractal |
Also Published As
Publication number | Publication date |
---|---|
TW201303319A (en) | 2013-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI425226B (en) | Method and system for fault detection, identification and location in high-voltage power transmission networks | |
Ferreira et al. | A survey on intelligent system application to fault diagnosis in electric power system transmission lines | |
Hajian et al. | A new hybrid pattern recognition scheme for automatic discrimination of power quality disturbances | |
Rizeakos et al. | Deep learning-based application for fault location identification and type classification in active distribution grids | |
Gilanifar et al. | Multi-task logistic low-ranked dirty model for fault detection in power distribution system | |
CN109975673A (en) | A kind of photovoltaic micro DC side fault electric arc recognition methods | |
CN108535572B (en) | Metering system secondary circuit monitoring method and device based on fundamental wave zero sequence characteristics | |
CN103018660A (en) | Multi-fault intelligent diagnosing method for artificial circuit utilizing quantum Hopfield neural network | |
CN117272143A (en) | Transmission line fault identification method and device based on Gram's angle field and residual network | |
CN114397569A (en) | A fault arc detection method for circuit breakers based on VMD parameter optimization and sample entropy | |
Cui et al. | High-impedance fault detection method based on sparse data divergence discrimination in distribution networks | |
CN109782158B (en) | An Analog Circuit Diagnosis Method Based on Multilevel Classification | |
CN118885944B (en) | Multi-mode fault diagnosis and early warning method | |
Lan et al. | Unsupervised power system event detection and classification using unlabeled PMU data | |
de Alencar et al. | A fault recognition method for transmission systems based on independent component analysis and convolutional neural networks | |
Dang et al. | seq2graph: Discovering dynamic non-linear dependencies from multivariate time series | |
CN118410415A (en) | Power System Fault Diagnosis Method Based on MP-Convformer Parallel Network | |
CN104655991A (en) | Power system fault matching method based on mutant point dejection combinational algorithm | |
CN118193986A (en) | A method for early fault prediction in distribution network based on periodic time series | |
Fayyad et al. | Transmission line protection using high-speed decision tree and artificial neural network: A hardware co-simulation approach | |
CN117630797A (en) | Electric meter health status detection method, system and storage medium based on working current | |
Sinha et al. | Power system fault detection using image processing and pattern recognition | |
Soni et al. | Deep learning technique for recurrence plot-based classification of power quality disturbances | |
CN114690070A (en) | Fault identification method based on Hilbert-Huang transform and support vector machine | |
Yu et al. | Knowledge reasoning of transmission line component detection using CRITIC and TOPSIS approaches |