# Patent

## Title
Detecting power quality events in power distribution networks

## Abstract

Power quality (PQ) events in a waveform in a power distribution network are detected by sampling the waveform. For each sample, a normalized sample is determined by subtracting an undistorted waveform from the sample. An accumulated weighted log-likelihood ratio (LLR) is determined for all normalized samples, as well as a minimum of the accumulated weighted LLR. If the difference between the accumulated weighted LLR an the minimum LLR is greater than a predetermined threshold, then the PQ event is signaled. The network can be a smart grid.

## Background
FIELD OF THE INVENTION
This invention relates generally to monitoring power in power distribution networks, and more particularly to detecting power quality events.
BACKGROUND OF THE INVENTION
In the United States, power outages and power quality problems annually cost hundreds of billions each year. Being motivated by this concern, one of the defining characteristics of the emerging smart grid is to support more stable and higher-quality power supply by using information technology.
To assess power quality (PQ), it is a common practice to monitor the quality of voltage and current waveforms by analyzing samples acquired by sensors installed in the power distribution networks.
In contrast with the sinusoidal power waveform generated by electric utilities, power waveforms over transmission lines and at consumer equipment are often distorted. Generally, distortions can be classified a PQ variations and PQ events. While PQ variations are characterized by small and gradual deviations from the sinusoidal voltage and current waveforms, PQ events incur large waveform deviations. PQ events are more detrimental to the power distribution network because the events can potentially inflict severe damages to the electrical equipment and injuries the consumers. Consequently, the occurrence of the PQ events has to be detected accurately and timely to allow appropriate actions.
In practice, PQ event monitoring includes detection and classification. During detection, a PQ event is declared when the waveform distortion exceeds a pre-defined threshold. Then, the distorted waveform can be classified to identify the cause of the PQ event before further analysis is performed.
Three primary PQ event detection methods are known. The first method tracks the root mean squared (RMS) value of the voltage waveform over a moving time window. The likelihood of the occurrence of the PQ event is evaluated based on the RMS change across windows. The RMS-based method is effective in detecting amplitude-related distortions. The second method detects the distortion in the frequency domain by transforming the time waveform into a frequency waveform using either wavelets or a short-time Fourier transform (STFT). The third method decomposes the waveform into a sum of damped sinusoids using super-resolution spectral analysis techniques, such as signal estimation via rotational invariance techniques, or multiple signal classification. The distorted waveform is detected by comparing the decomposed frequency-domain components of a monitored waveform with those of the normal waveform.
The moving time window segments the waveform into blocks before any transformation or decomposition is applied. Therefore, the time resolution of all three conventional methods is restricted by the size of the moving window.
If the window size is sufficiently large to meet the detection rate and the false alarm rate requirements, the delay in detecting the PQ window is increased, and real-time detection is not possible.
The invention corrects this problem.
SUMMARY OF THE INVENTION
The embodiments of the invention provide a method and apparatus for detecting power quality (PQ) events in a waveform of a power distribution network. The method is based on change-point detection. Generally, change-point detection discovers time points in time series-samples at which properties of time-series samples change, such as changes in the waveform due to the PQ events.
The method uses a sequential cumulative sum (CUSUM) procedure for detecting the PQ events in real-time. The CUSUM procedure evaluates weighted likelihood ratios based on instantaneous and long-term samples of the power waveform. The invented method achieves a significant performance gain over conventional PQ event detection schemes.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a block diagram of a method for detecting a power-quality (PQ) event in a waveform in a power distribution network according to embodiments of the invention;
FIG. 2 is graph of a transient voltage event detected by the method according to embodiments of the invention;
FIG. 3 is a graph of a temporal-frequency using SIFT with respect to the transient event detected according to embodiments of the invention;
FIG. 4 is a conventional spectral estimates for a transient PQ event;
FIG. 5 is a graph of a conventional sample-by-sample RMS of the transient PQ event;
FIG. 6 is a graph of a logarithm of the weighted likelihood ratio of the transient PA event detected by the method according to embodiments of the invention; and
FIG. 7 is a graph of a voltage sag event detected by the method according to embodiments of the invention.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
The embodiments of the invention provide a method detecting a power quality (PQ) event in a waveform in a power distribution networks. In one embodiment the network is a smart grid. The method periodically samples the waveform.
If the PQ event occurs at time t=te, a goal of the method is to detect the PQ event with a minimum delay and a highest detection accuracy. The method can also to detect an end of the PQ event. The focus of this description is on detecting the start of the PQ event. It would be obvious to those of arinary skill in the art to extend the method to detecting the end of the event.
Waveform Before PQ Event
Samples of the continuous-time waveform before the PQ event are represented by

y(t)=sθ0(t)+n(t), 0≦t<te,  (1)

where n(t) is additional white Gaussian noise (AWGN) with zero-mean and variance σn2, denoted by N(0,σn2), and

sθ0(t)=a0·sin(2πf0t+φ0),  (2)

is an undistorted power waveform with θ0[a0, f0, φ0]T, where a0=1 is a signal amplitude gain, and f0 and φ0 are a fundamental frequency and an initial phase of the power waveform, respectively.

Waveform After PQ Event
Samples of the power waveform after the PQ event are represented by

y(t)=sθ1(t)+n(t), t≧te,  (3)

where θ1[a1, f1, φ1, φT]T and

sθ1(t)=a1·sin(2πf1t+φ1)+ξ(t),  (4)

with ξφ(t) being the additive distortion parameterized by φ. Eq. (3) represents a generalized power waveform that takes a typical PQ events into consideration.

For example, voltage sags can be modeled as sudden decrease in the waveform amplitude gain with a1<a0 while setting ξφ(t)=0. In contrast, a transient voltage event can be described by non-zero ξφ(t) with f0=f1 and φ0=φ1.
A before probability density functions (PDF) pθ0(y) and an after PDF pθ1(y) represent the distribution of the sample y before and after the PQ event, respectively. The before PDF pθ0(y) can be estimated due to the fact that {a0, f0, φ0} are deterministic and σn2 can be accurately measured.
In contrast, the after PDF pθ1(y) depends on a specific type of PQ event. It is generally difficult to fully characterize the after PDF pθ1(y) before the occurrence of the PQ event. This is a problem for a conventional statistical hypothesis test methods that use Neyman-Pearson hypothesis testing.
As a result, most conventional PQ event detection methods are designed to only use the instantaneous changes in the waveform, such as changes in amplitude, frequency or phase without utilizing associated long-term statistics.
For example, the conventional RMS method concentrates on amplitude changes by sampling and determining the RMS of the waveform. Let yk be the k-th current sample of the waveform. The conventional RMS method tracks the RMS of the samples over a moving window of size N, where N usually covers one cycle of the power-system frequency. The m-th RMS is




  
    
      Y
      rms
    
    ⁡
    
      (
      m
      )
    
  
  =
  
    
      
        
          1
          N
        
        ⁢
        
          
            ∑
            
              k
              =
              
                m
                -
                N
                +
                1
              
            
            m
          
          ⁢
          
            
          
          ⁢
          
            y
            k
            2
          
        
      
    
    .
  




Conventionally, the PQ event is detected when the current RMS value change is larger than a predetermined threshold. In addition to the time-resolution problem associated with the size of the moving window, conventional methods are sub-optimal due to the fact that they do not use statistical distributions of the samples before and after the PQ event.
The PQ event dection method according to embodiments of the invention uses a cumulative sum (CUSUM) procedure. The PDF pθ0(y) before the PQ event is known while the post-event PDF pθ1(y) needs to be determined. To circumvent the uncertainty of the after PDF, the weighted CUSUM procedure is used to replace the conventional log-likelihood ratio (LLR) test. In the following, a0=1, f0 is either 50 or 60 Hz, and φ0 and σn2 are known.
PQ Event Detection Method
FIG. 1 shows the steps of the method for detecting and signaling a PQ event 109 in the waveform 102 of the power distribution network 100 according to embodiments of the invention. The steps of the method can be performed in a processor 105 including memory and input/output interfaces as known in the art. It is understood that the method can also be implemented in hardware, e.g., an application-specific integrated circuit (ASIC). The steps are performed for each current sample y 101 sensed by a sensor. The sensor can be designed to sample, voltage, amplitude, phase frequency or combinations thereof.
Pre-Event PDF
Because {a0, f0, φ0 } are known, the undistorted waveform sθ0 before the PQ event is also known. Each sample 100 y(t) in Eq. (1) is transformed 110 into a normalized sample z(t) by subtracting the undistorted waveform from the current sample to obtain the normalized sample

z(t)=y(t)−sθ0(t)=n(t), 0≦t<te.  (6)

Thus, the before PDF of the normalized sample z is

pθ0(z)=N(0,σn2),  (7)

Post-Event PDF
The after PDF is derived using a change-point detection procedure as known in the art, see generally U.S. Pat. No. 7,016,797 “Change point detection apparatus, method and program therefor,” issued to Takeuchi et al. on Mar. 21, 2006, and incorporated herein by reference.
The change point detection procedure according to embodiments of the invention uses a weighted CUSUM procedure. Therefore, to obtain the normalized sample z(t), Eq. (3) is converted to subtract the undistorted wavefrom from the current sample according to

z(t)=y(t)−sθ0(t)=x(t)+w(t),  (8)

where

x(t)=a1·sin(2πf1t+φ1),  (9)

w(t)=ξφ(t)−sθ0(t)+n(t).  (10)

Because the after θ1 is unknown, rather than evaluating the LLR ratio





  
    
      p
      
        θ
        1
      
    
    ⁡
    
      (
      
        z
        i
      
      )
    
  
  
    
      p
      
        θ
        0
      
    
    ⁡
    
      (
      
        z
        i
      
      )
    
  




directly, the weighted LLR is determined with the weighted CUSUM method as





  
    
      
        
          
            s
            i
          
          =
          
            ln
            ⁡
            
              [
              
                
                  ∫
                  
                    θ
                    1
                  
                  
                    
                  
                
                ⁢
                
                  
                    
                      
                        p
                        
                          θ
                          1
                        
                      
                      ⁡
                      
                        (
                        
                          z
                          i
                        
                        )
                      
                    
                    
                      
                        p
                        
                          θ
                          0
                        
                      
                      ⁡
                      
                        (
                        
                          z
                          i
                        
                        )
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    
                      F
                      ⁡
                      
                        (
                        
                          θ
                          1
                        
                        )
                      
                    
                  
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        11
        )
      
    
  




The PDF of w of is approximated as N(0,σn2) using the well known central limit theorem, where σw2=σξ2+σn2+a02. Furthermore, recall that x(t) is approximately uniformly distributed over [−|a1|,+|a1|]. Thus, it is straightforward to show that the after PDF of the normalized samples is




  
    
      
        
          
            
              p
              
                θ
                1
              
            
            ⁡
            
              (
              z
              )
            
          
          =
          
            
              1
              
                4
                ⁢
                
                  
                  
                    a
                    1
                  
                  
                
              
            
            ⁡
            
              [
              
                
                  erf
                  (
                  
                    
                      z
                      +
                      
                        
                        
                          a
                          1
                        
                        
                      
                    
                    
                      
                        2
                      
                      ·
                      
                        σ
                        w
                      
                    
                  
                  )
                
                -
                
                  erf
                  (
                  
                    
                      z
                      -
                      
                        
                        
                          a
                          1
                        
                        
                      
                    
                    
                      
                        2
                      
                      ·
                      
                        σ
                        w
                      
                    
                  
                  )
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        12
        )
      
    
  




where erf is an error function.

With the assumption that x(t) and w(t) are statistically independent, F(θ1) can be expressed as

F(θ1)=F(a1)·F(σw).  (13)

As a result, the weighted LLR in Eq. (11) becomes




  
    
      
        
          s
          i
        
        =
        
          
            ln
            [
            
              
                ∫
                
                  a
                  1
                
              
              ⁢
              
                
                  ∫
                  
                    σ
                    w
                  
                
                ⁢
                
                  
                    
                      
                        p
                        
                          θ
                          1
                        
                      
                      ⁡
                      
                        (
                        
                          z
                          i
                        
                        )
                      
                    
                    
                      
                        p
                        
                          θ
                          0
                        
                      
                      ⁡
                      
                        (
                        
                          z
                          i
                        
                        )
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    
                      F
                      ⁡
                      
                        (
                        
                          σ
                          w
                        
                        )
                      
                    
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    d
                    
                      F
                      ⁡
                      
                        (
                        
                          a
                          1
                        
                        )
                      
                    
                  
                
              
            
            ]
          
          .
        
      
    
    
      
        (
        14
        )
      
    
  




The distributions F(.) are generally in the form of a Gaussian distribution.

From the before and after PDFs in (7) and (12) respectively, the weighted LLR of the sample at time I can be evaluated according to Eq. (14).
From the LLR, the method 120 determines the accumulated weighted LLR for all samples as




  
    
      
        
          
            S
            k
          
          =
          
            
              ∑
              
                i
                =
                1
              
              k
            
            ⁢
            
              
            
            ⁢
            
              s
              i
            
          
        
        ,
      
    
    
      
        (
        15
        )
      
    
  




that is, the accumulation is from the first data sample i=1 to the current sample k. The accumulated weighted LLR is based on the ratios of the before and after PDFs 125.

Next, the method determines 130 a minimum of the accumulated weighted LLR Sj for j=1, 2, . . . , k as




  
    
      
        
          
            m
            k
          
          =
          
            
              min
              
                1
                ≤
                j
                ≤
                k
              
            
            ⁢
            
              S
              j
            
          
        
        ,
      
    
    
      
        (
        16
        )
      
    
  




where the function min returns a minimum value.

Next, the method determines 140 a difference between the current accumulated weighted LLR and the minimum accumulated weighted LLR

gk=Sk−mk.  (17)

Then, the method determines 150 if the difference gz is greater than a predetermined threshold h. If the difference gk is greater than the threshold, i.e., gk>h, then the method signals a PQ event 109. Otherwise, the next sample is processed.
FIG. 2 shows a transient voltage event 200 starting at ime te=0.06 s.
FIGS. 3 and 4 respectively show temporal-frequency plots using STFT and spectral estimates over multiple windows, respectively. As shown in FIGS. 3 and 4, it is difficult to detect the PQ event directly from these plots.
FIGS. 5 and 6 show the sample-by-sample RMS generated by the conventional RMS scheme, and the logarithm of the weighted likelihood ratio by the CUSUM method according to the embodiments of the invention. The CUSUM method clearly has much stronger indication on the PQ event occurrence at t=0.06 s.
Voltage Sags
FIG. 7 shows a voltage sag detected by the method according to embodiments of the invention. The sag is due to a temporary ground fault.
EFFECT OF THE INVENTION
The invention provides a method for detecting a power quality event using a change-point detection theoretic framework. The method analysis a difference of statistical distributions of power waveforms before and after the PQ event occurrence. Because the proposed scheme performs sample-by-sample evaluation, it can achieve the detection task with the smallest time resolution. Furthermore, capitalizing on the change-point detection theory, the proposed scheme provides accurate detection performance, irrespective of the availability of the post-event statistical distribution. Although the change-point detection theory has been successfully applied to a wide variety of applications such as spectrum sensing in cognitive radio networks, as best s is known, the method according to the embodiments of the invention is the first use the change-point detection theory to the PQ monitoring problem.
Although the invention has been described by way of examples of preferred embodiments, it is to be understood that various other adaptations and modifications may be made within the spirit and scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true spirit and scope of the invention.

## Claims

1. A method for detecting a power quality (PQ) event in a waveform in a power distribution network by sampling the waveform, comprising for each sample the steps of:
subtracting an undistorted waveform from the sample to obtain a normalized sample;
determining an accumulated weighted log-likelihood ratio (LLR) for all normalized samples;
determining a minimum of the accumulated weighted LLR;
determining a difference between the accumulated weighted LLR and the minimum LLR;
signaling the PQ event if the difference is greater than a predetermined threshold, where the steps are performed in a processor.

2. The method of claim 1, wherein the network is a smart grid.

3. The method of claim 1, wherein the accumulated weighted LLR uses change-point detection.

4. The method of claim 1, wherein the accumulated weighted LLR uses a sequential cumulative sum procedure.

5. The method of claim 1, wherein the samples of the waveform before the PQ event are represented by


y(t)=sθ0(t)+n(t), 0≦t<te,


where n(t) is additional white Gaussian noise (AWGN) with zero-mean and a variance σn2, denoted by N(0,σn2), and


sθ0(t)=a0·sin(2πf0t+φ0),


is an undistorted power waveform with θ0[a0, f0, φ0]T, where a0=1 is a signal amplitude gain, and f0 and φ0 are a fundamental frequency and an initial phase of the waveform, respectively, and wherein the samples of the waveform after the PQ event are represented by


y(t)=sθ1(t)+n(t), t≧te,


where θ1[a1, f1, φ1, φT]T and


sθ1(t)=a1·sin(2πf1t+φ1)+ξφ(t),


with ξφ(t) being an additive distortion parameterized by φ.

6. The method of claim 1, wherein the PQ event is a voltage sag.

7. The method of claim 1, further comprising:
detecting an end of the PQ event.

8. The method of claim 5, wherein a before probability distribution function (PDF) of the samples before the PQ event is pθ0 (y), and wherein the accumulated weighted LLR is based on the before PDF and the after PDF.

9. The method of claim 8, wherein the normalized sample


z(t)=y(t)−sθ0(t)=n(t), 0≦t<te.

10. An apparatus for detecting a power quality (PQ) event in a waveform in a power distribution network, comprising:
a sensor configured to sample the waveform;
means for subtracting an undistorted waveform from each sample to obtain a normalized sample;
means for determining an accumulated weighted log-likelihood ratio (LLR) for all normalized samples;
means for determining a minimum of the accumulated weighted LLR;
means for determining a difference between the accumulated weighted LLR and the minimum LLR; and
means for signaling the PQ event if the difference is greater than a predetermined threshold.