Compressive Sensing-based Optimal Reactive Power ... - IEEE Xplore

0 downloads 0 Views 1MB Size Report
Sep 1, 2017 - Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE ...... pp. 1-4, 2014. [18] O. Alizadeh Mousavi and R. Cherkaoui, "Maximum ...
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301

Compressive Sensing-based Optimal Reactive Power Control of a Multi-area Power System Irfan Khan, Student Member, IEEE, Yinliang Xu, Member, IEEE, Soummya Kar, Member, IEEE and Hongbin Sun, Senior Member IEEE  Abstract—To deal with the increasing load demand and environmental effects of conventional power devices, power system has become an enlarged complex network with the integration of distributed generators. Real-time control of power system now needs a massive amount of information to be transmitted between each local device and the central controller. The transmission of a huge quantity of information data poses great challenge to the communication network. To deal with this issue in reactive power control, this paper proposes a novel realtime compressive sensing based optimal reactive power control of a multi-area interconnected power system. The objective is to minimize the power loss, voltage deviation and reactive power generation cost simultaneously. According to the proposed scheme, the measured data in each control area is compressed before being transmitted through the communication network, and then recovered accurately by the discrete central controller. Orthogonal matching pursuit algorithm is adopted to recover the compressed data, owing to its fast convergence speed. Simulation results demonstrate the effectiveness of the proposed compressive sensing based approach by significantly reducing the data size of the transmitted data. Index Terms—Optimal reactive power control, compressive sensing, orthogonal matching pursuit algorithm, power loss, voltage deviation.

F

I. INTRODUCTION

or the reliable and stable operation of the power system, information of real-time parameters, state, and control variables of the whole power system shall be collected and processed. With the recent large-scale integration of Distributed Generation (DG) and introduction of smart meters in the power system, the volume of data that needs to be processed, has grown drastically [1]. According to [2], in 2014, 94 million Manuscript received June. 19, 2017; accepted Sept 1, 2017. This work was supported by the National Natural Science Foundation of China (Grant No. 51537006, 51507193), (Corresponding Author: Yinliang Xu). I. Khan is with Electrical and Computer Engineering Department, Carnegie Mellon University, Pittsburgh, PA, 15213 USA, and School of Electronics and Information Technology, Sun Yet-sen University, Guangzhou, 510275, China. (e-mail: [email protected]) Y. Xu and H. Sun are with Tsinghua-Berkeley Shenzhen Institute (TBSI), Shenzhen, Guangdong, 518055 China, and the Department of Electrical Engineering, State Key Laboratory of Power Systems, Tsinghua University, Beijing, 100084 China (e-mail: [email protected], [email protected]). S. Kar is with Electrical and Computer Engineering Department, Carnegie Mellon University, Pittsburgh, PA, 15213 USA (email: [email protected]).

smart meters were shipped worldwide and it is predicted that by 2022, the number of smart meters will reach 1.1 billion. When more and more devices are being connected to the same network segments, bandwidth becomes more limited with more network participants [2]. The synchrophasor data is reported on an N frames/s basis, where N can be selected between 1 and 120. Choosing a higher frame rate would further accelerate a demand for more bandwidth [3]. Power utilities often find it difficult to arrange large investments to build/upgrade communication networks to meet the increased bandwidth demand [4]. Therefore, in this context, one must focus to achieve an efficient use of the existing communication channel bandwidth. It is exhibited in [5] that the energy consumption by the communication channel increases significantly as the size of message to be transmitted increases in various applications of power grid. Furthermore, if the data package exceeds the maximum allowable data size, it needs to be divided into several smaller data packets, which may lead to considerable time delay. Therefore, it is advisable to decrease the size of realtime data for obvious reasons. Reactive power control plays a vital role on power loss reduction and voltage control of the power system [6]. Under high penetration of new DG units, more and more reactive power control devices are being integrated into power grids. The large-scale deployment of these devices requires upgrading the existing reactive power control solutions [7]. Many conventional centralized algorithms have been employed to solve the optimal reactive power control (ORPC) problem offline [8, 9]. However, only few attempts have been made to control reactive power in real-time applications, which needs to upgrade the existing communication networks to meet with the increasing bandwidth demand. To address the issue of fast two-way communication, decentralized ORPC techniques, which rely on local controllers only without explicit communication, have been proposed in [10, 11]. An approximate measure for power loss is considered by controlling the reactive power generation from each DG in [10]. The authors in [11] demonstrated that total voltage deviation is higher using decentralized ORPC than that of using centralized ORPC. It is because of the reason that using decentralized control, each local DG controller tries to satisfy its local load demand while load buses without having their own reactive power generation sources, may happen to have large voltage deviation. Thus, due to the lack of broader available information, decentralized control may not be effective to coordinate all available resources in the network [12].

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 Recently, distributed ORPC techniques [7, 13, 14], which take advantage of sparse communication network, have been very popular among researchers. In [13], voltage regulation and power loss were minimized by adjusting the reactive power generation using a feedback strategy based distributed control algorithm. Minimization of the approximate power loss was proposed by controlling the bus voltages in [14], whereas in [7], authors proposed distributed multiple agent system based ORPC to optimize power loss and voltage deviation. It is shown that the data loss and communication delays may slow down the overall converging speed. Moreover, distributed algorithms require synchronized communication between the neighboring entities i.e. failing of the communication or delaying of information from even a single bus may slow down the distributed algorithm. When comparing distributed and centralized schemes, it is found that centralized protocols show better performance in data loss recovery than the distributed protocols [15]. Also, the optimal solution of the centralized optimal control is found to be better than that of distributed optimal control [16, 17]. The reason may be that in distributed control, local information is shared among neighboring entities only, whereas in the centralized control, all the information of the system is known which ensures the better optimal control of reactive power. Also, most of the existing reactive power control solutions are centralized [7] and it is economical to improve the existing communication techniques than to completely reshape the overall communication paradigm. It may be expensive to deploy a controller on each DG bus. Reactive power control service is mainly contemplated as a local control service because the reactive power may not be transmitted efficiently through long distances in transmission networks. Thus, it is usually suggested to control the voltage by using control devices dispersed throughout the power system. Therefore, system operators normally provide voltage control services from resources within their own controlled area [18, 19]. Although voltage control is primarily a local control problem, the recent two decades widespread blackouts have demonstrated that the voltage instability and collapse could be considered as an important factor in major power outages worldwide. It may involve several areas of an interconnected system and increase the scale of blackouts and even affect the intact areas [20-21]. The multi-area power system without inter-area voltage coordination may be operated in a non-optimum state which means less security and stability margin. For instance, [22] reveals that in multi-area power systems, the optimization solution of an area for reducing the active power losses in its own region might lead to increase the losses globally in the interconnected system. Additional literature about the coordinated reactive power control in multi-area power system can be found in [23-29]. For centralized schemes, ORPC problem for a large power system requires a huge amount of data to be transmitted to the central controller. If the bandwidth of the communication network is not sufficiently large to support transmission of such massive amount of data, it may deteriorate the performance of the controller. To address the problem of limited bandwidth, recently, compressive sensing (CS) which takes advantage of

sparse signal structures is considered as a promising joint data compression and reconstruction method to transmit information efficiently [30]. As the data obtained from observing a physical system like power system, inherently contains spatial and temporal correlation [31], it can be compressed before transmitting through the unreliable communication channel and recovered at the receiving end. A CS-based data loss recovery algorithm to improve the communication efficiency by reducing the communication burden and thus, saving the power demand is presented in [32], which is further improved in [33] with low space complexity, low floating-point calculations, and low time complexity robust CS framework. Recently, a data reduction scheme by removing redundant phasor information before transmission is proposed in [3] in order to reduce the amount of data. To the authors’ best knowledge, few attempts have been made to deploy advanced information processing techniques in the ORPC problem. In this paper, a centralized real-time ORPC is proposed for multi-area interconnected power system. Realtime data of bus voltages and line currents in each area is first sampled and compressed before being transmitted through the communication network. At the receiving end, it is first decoded and recovered accurately by the central controller. The size of the transmitted data can be significantly reduced with the use of CS. The proposed approach is reliable, robust and its effectiveness is validated through simulations. Major contributions of this paper are as follows: 1) Real-time optimal reactive power control of multi-area power system is solved to minimize power loss, voltage deviation and reactive power cost, simultaneously. 2) To meet with the challenges of high bandwidth requirement for real-time transmission of massive data of power system, compressive sensing technique is utilized to first compress the data package before transmitting it through the communication network and then recover it at the end of communication channel. 3) To deal with large control areas with huge amount of data, an Orthogonal Matching Pursuit algorithm is used to improve the reconstruction speed of the transmitted data. 4) Scalability of the proposed approach is investigated. The rest of the paper is organized as follows. Section II formulates the CS based ORPC problem. Section III explains the working of compressive sensing. Section IV presents the proposed algorithm. Section V discusses the simulation results of the proposed control approach, and Section VI concludes the paper. II. PROBLEM FORMULATION Optimal reactive power control of DGs plays a key role in the power system control and operation, which can lead to improved voltage profile, minimal active power loss and controllable reactive power generation cost. Therefore, the objective function to be optimized is formulized as: f1 (QGi ,Vi , i )  W1 Ploss  W2 DV  W3CQ min (1) QGi

where W1, W2 and W3 are the weight coefficients. Ploss, DV and CQ are the power loss, voltage deviation and cost of reactive

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 power generation, respectively. QGi is the total reactive power generation from each bus i. Vi and δi are the bus voltage magnitude and its angle on bus i. Expressions for the three sub-functions are given as: n

Ploss 

n

 V V Y cos( i

j ij

ij

  j   i ) [34]

(2a)

i 1 j 1

n

DV 

 (V  V ) i

* 2 i

(2b)

i 1

n

CQ 

a Q i

i 1

2 Gi

 bi QGi  ci [35]

(2c)

where Yij ij is the element of admittance matrix. Vi* is the reference voltage on bus i, whereas ai, bi and ci are the reactive power cost coefficients. To minimize (1), admittance matrix, reference voltage and reactive power generation cost coefficients are usually fixed whereas, voltage phasors, and real & reactive power flows are the state measurements of the control area of power system as given in (3). x i = Vi  i PGi QGi  . (3) For minimization of (1) by a centralized approach, the state information is required to be transmitted to the central controller through communication channel. However, for real time ORPC, the communication channel bandwidth is usually limited and cannot execute the data transfer beyond its limit. To deal with this issue, in this paper, a compressive sensing based ORPC algorithm is proposed where state measurements data is first compressed before sending on to the communication channel and is recovered at the end of communication channel by the central controller. The theory of encoding and decoding of measurement signals using CS is explained in Fig. 1. At each time step k, the state measurements x(k) are taken from n sensors of the power system control area and compressed as (4) y ( k )  Φ ( k ) x( k ) (4) where x(k) is the state information data of the control area, y(k) is the data after compression and Φ(k)∈Rm*n , is the measurement matrix at time step k where m is the number of samples after encoding while n is the number of samples before encoding. The value of n/m defines the compression ratio (ρ) of the CS. DATA DECODING

RECONSTRUCTION ALGORITHM

DATA ENCODING

= Compressed Samples mx1

Measurement Matrix mxn

Recovered Data nx1

Sending Data nx1

Fig. 1 Encoding and decoding of measurement signal using CS

At the end of communication channel the original n samples of x(k) can be recovered from m samples of y(k) measurements by solving ℓ1-norm minimization problem as given in (5) [36]

xˆ ( k )  arg min x( k ) x(k )

1

. (5) y ( k )  Φ ( k ) x( k ) s.t. where xˆ (k ) is the recovered state information data. A number of methods can be used to solve the optimization problem in equation (5) [37, 38]. For real-time data recovery, matching pursuit methods [38] are preferred to prioritize the computational speed. Matching pursuit approximates the solution to Eqn. (5) by solving with Orthogonal Matching Pursuit Algorithm (OMP) as given in Eqn. (6) 2

xˆ (k )  arg minx( k ) y(k )  Φ(k )x(k ) 2 s .t . x 0   (6) where  is the bound on the ℓ0-norm of the sparse vector x. Now recovered xˆ (k ) will be used by the central controller to generate ORPC control input for the control area of power system. x(k)

Compression

y(k)

y(k)

Recovery Algorithm

Communication Channel with Limited Bandwidth

Control Area of power system

u(k )

xˆ(k)

Central Controller

u(k)

Fig.2 CS based centralized control scheme

The working of CS based centralized control strategy is explained in Fig. 2, which shows the compression of x(k) before transmitting it through communication channel. At the end of communication channel it can be recovered using OMP recovery algorithm as proposed in (6). III. DETAILED WORKING OF COMPRESSIVE SENSING In this paper, matrices are represented by bold and italic letters, vectors are represented by bold letters whereas scalar quantity is neither bold nor italic. The theory of CS states that a sparse or compressible signal can be recovered with high probability from a few measurements, which is far smaller than the size of the original signal. Power system state signals, having spatial-temporal correlation, thus, can be regarded as compressible. For a control area with n sensors, a data window X(k) which contains w≥1 consecutive readings of n sensors at the time step k is given as  x1 ( k  w  1)  x1 ( k )  X ( k )       (7)  x n ( k  w  1)  x n ( k )  Denote the ith row of X(k) by xr,i (k) =[xi(k-w+1),…, xi(k)]T, which contains readings of sensor i at time steps {k-w+1,…, k}. Similarly, denote the jth column of X(k) by xc,j(k-w+j) =[x1(kw+j) ,…, xn(k-w+j)]T, which contains the readings of all sensors i at the time step k-w+j. If there exists a basis ΨS ∈ Rn*n and a basis ΨT ∈ Rw*w for spatial and temporal domains respectively, X(k) can be written in a compressible expression using spatial basis as (8) X (k )  Ψ S θ S (k ) (8) θ S (k )  [θ S (k - w +1),...., θ S (k )] nw (9) where θS(k) ∈Rn contains the spatial transform coefficients of n sensors at the time step k and θS(k) is its matrix for the whole

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 data window of w. Similarly, X(k) can also be written in a compressible expression using its temporal basis as X T (k ) = Ψ T θT (k ) (10)

θ T (k )  [θ T,1 (k )....θ T,n (k )] wn (11) where θT,i (k) ∈Rw contains the temporal transform coefficients of sensor i for all w steps. Kronecker sparsifying bases can succinctly combine the individual sparsifying bases of each signal dimension into a single transformation matrix [36, 39]. Hence, we can merge the transformations in (8) and (10), and represent X(k) as x(k )  vec( X (k ))  vec(Ψ S θS (k ))

 vec(Ψ S θS (k )Ψ T-T Ψ TT )  vec(Ψ S Ζ (k )Ψ TT )  (Ψ S  Ψ TT )vec( Ζ (k ))

(12)

 Ψz (k )

is the Kronecker sparsifying where Ψ = (ΨT ⊗ ΨS) ∈ R basis and z(k) contains the joint transformation domain coefficients. x(k)=[xTc,1(k-w+1), …, xTc,n(k)]T is vector– reshaped data window of X(k), Z(k) can be interpreted as the matrix representation of spatial domain θS(k) in temporal basis ΨT, i.e. θTS(k)= ΨT ZT(k). It should be noted that any universal data independent basis can be considered as ΨS and ΨT. Some of the transformations found suitable for revealing the underlying sparsity of many natural signals include the discrete Fourier transform (DFT) basis, discrete cosine transform (DCT) basis and the discrete wavelet transform (DWT), etc. x(k) can be represented by ΨT z(k). in (4) and (6) for encoding and decoding of the system data respectively. For recovery of the state measurements by (6), initially zˆ (k ) will be recovered and then xˆ (k ) can be reconstructed using (13) xˆ (k )   zˆ (k ) (13) nw×nw

Start

End

Choose Compression Ratio (

Sent to the central controller

Input measurement data X(k) and vectorize it as x(k) Select

and

Reconstruct (k) as per (13)

Recover the data (k) as per (6)

2. Signal to error ratio (SNR): SNR  20(log x  xˆ ) x

(15)

ˆ is the recovered signal. The where x is the original signal, x compression ratio (ρ) decides the size of the measurement matrix. It decides the number of rows (m) in the measurement matrix as given in (4). Higher the compression ratio, lower the amount of data to be transmitted on the communication channel. IV. PROPOSED ORPC ALGORITHM Optimization of the multi-objective function of (1) can be performed using optimal reactive power control of the generators. However, the inequality constraints of reactive power generation limits should be included to ensure that the reactive power generator limits are not violated. Here, reactive power optimization can be solved using Interior Point Method (IPM) algorithm [40, 41], as given in (16) [42]. min f (x) x

n

PGi  PLi  Vi

(16)

x £ xu x ³ xl

s.t.

V Y

j ij

cos( i   j  ij )  0

(16a)

j 1

n

QGi  QLi  Vi

V Y sin(   j ij

i

j

 ij )  0 (16b)

j 1

where xui and xli are the maximum and minimum rated 0 is the capacities of state variable for bus i. x=[x1, …, xn], vector of barrier parameters which keeps on decreasing until the optimal solution is achieved. (16a & b) are the power flow constraints which need to be satisfied during the control operation and can be written as c(x) in a compact form. These constraints can be included in the formulation using lagrangian approach. Lagrangian of (16) can be written as (17)

L(x,  w l , w u )  f (x)  c(x) T   w Tl (x  x l )  w Tu (x u  x) (17) (16) consists of nonlinear equality constraints. It is convenient to utilize a Primal-dual based Interior Point method (PDIPM) to solve nonlinear optimization problems with equality constrains. Hence, the primal-dual equation of (16) is given as (18) [43] g ( x)  A ( x)   w l  w u  0 c ( x)  0 ( X  X l )w l  e l  0

(18) (18a) (18b)

( X u  X )w u  e u  0 (18c) Transmit y(k) through Compress x(k) to where g(x) is the dual of f1(x).  is the vector of equality limited bandwidth calculate y(k) as per (4) communication channel constraint dual variables. w l and w u are the vectors of nonFig. 3 Proposed CS based data transmission algorithm negative dual variables of inequality constraints. Xu and Xl are the matrix of maximum and minimum limits of state variables Flowchart of the proposed CS based data transmission for a specified window length as given in (7). e and e are two l u algorithm is shown in Fig. 3, which shows the steps required to indicator variables to restrict states within feasible regions. transmit the data using CS. For a fixed u at each iteration and counted by k, the solution The performance of CS is mainly evaluated by two criterions: of (16) can be achieved by Newton method where search 1. Compression ratio (ρ): directions of states can be calculated as (19) [44] size of data being transmitted without CS  (14) size of data being transmitted with CS

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301  H (k ) A(k )   d x ( k )  f (x(k ))  A(k ) (k )    A T (k )  (19) 0  d  (k )   c(k )   where d x (k ) and d  (k ) are vectors of search direction of x and λ . Here, H(k ) is a matrix given by (20)

H (k )  V (k )  Wl (k )( X (k )  X l )) 1 + Wu (k )( X u  X (k )) 1 (20) where V (k ) is the Hessian of Langrangian of (17). X(k) is the matrix of state information as given in (7). Wl and Wu are the matrices of non-negative dual variables of inequality constraints for a window length of w. In case of using CS to transmit network data to the central controller, rate of compression and recovery is usually very high. Thus, it is reasonable to approximate the gradient and Hessian of function of f1 numerically as L L( k )  L( k  1) (k )  (21) x x( k )  x( k  1) L L (k )  ( k  1)  2L (22)  x x (k )  x 2 x( k )  x( k  1) Similarly, search direction of the dual variables can be calculated as given in (23a) & (23b) dwl (k)  (k)(X(k)  Xl )1el Wl (k)  wl (k)(X(k)  Xl )1dx(k)) (23a) dwu (k)  (k)(Xu  X(k))1eu wu (k) Wu (k)(Xu  X(k))1dx(k)) (23b) Finally, update of next iteration of PDIPM can be obtained by (24(a-d)) x(k +1)= x(k ) +  (k )d x (k ) (24a)

λ (k +1)= λ (k ) +  (k )d λ (k )

(24b)

w l (k +1)= w l (k ) +  (k )d w l (k )

(24c)

w u (k +1)= w u (k ) +  (k )d w u (k )

(24d)

where  (k ) 0,1 is the step size. In classical PDIPM,  is decreased at each iteration by using a fixed decrement ratio which may slow down the convergence speed of the algorithm. In this paper, authors propose an improved version of PDIPM which decreases the barrier parameter at a super linear rate by using (25)   (25)  (k  1)  max  tol , min{   (k ),  (k )   ( k ) }    where  tol is the given tolerance,  5,15 ,   0,1 , and

reactive power control application. Using improved PDIPM, the central controller needs to collect the real-time data of voltage phasors and power flows which in this case is being achieved by using CS based communication network. After obtaining the real-time data, the central controller generates the optimal control input for each reactive power source in the control area. The optimal control input is the optimal reactive power generation and it is decided based on the minimization of (16). Start

Initialize iteration counter k

Calculate optimal solution as per (1)

Yes

as per (26) No

k=k+1

Update

as per (25)

End

Find the search direction as per (19) and (23)

Determine states as per (24)

Fig. 4 Proposed Improved PDIPM algorithm

Finally, a comprehensive flowchart of the proposed CS based ORPC algorithm is shown in Fig. 5. At the beginning of the algorithm, PDIPM and CS preset parameters are stored in the central controller. Measured data from control area of a power system is compressed as per Fig. 3, before transmitting it through the communication channel. For compression of measured data, two matrices Kronecker sparsifying basis (Ψ) and measurement matrix Φ(k) are required. Ψ is achieved by taking Kronecker tensor product of Discrete Cosine Transform basis. Measurement matrix Φ(k) depends on the selected compression ratio for each simulation. If compression ratio rises, number of the rows in the measurement matrix decreases. It is obtained dynamically for each iteration by randomly selecting a sparse matrix with only one non-zero value in each row. At the end of communication channel, zˆ (k ) is recovered using OMP algorithm applied on (27).

(27) y(k )  Φ(k ) zˆ (k ) It is an under-determined system of equations and is solved using the method of least squares to get an estimate of the solution. However, the information of the sparsity for zˆ ( k )   ( k ) 1,2 . To check the convergence of the proposed does help to derive a more deterministic solution set by improved PDIPM, error can be calculated by finding the incorporating λ with the optimization problem as given in (28) maximum violation in the primal-dual equations as given in (26) [45]. 2 min Φ ( k ) zˆ ( k )  y ( k ) 2  g(x)  A(x)  wl  wu   zˆ ( k ) , c(x)  ,   (28) sd E  max   (26) s.t . zˆ ( k ) 0    (x(k)  xl )wl  el scl , (xu  x)wu  eu scu      where sd, scl and scu are the scaling parameters. Detailed working of the proposed improved PDIPM algorithm is provided in Fig. 4. It is proved to be an efficient technique for nonlinear constrained optimization and thus, suitable for our optimal

where  is the number of non-zeros in the least squares problem

and optimal

 is calculated by solving (29) T   2 Φ (k )  y (k ) 

(29)

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 START

Initialization QGi,max, Ԑtol , β,κμ, ʈμ , ρ

Measuring real-time power system data Vi, δi, QGi, PGj

Compression of data package

connecting 3 control areas with 6- bus radial system in each area as shown in Fig. 7. The impedances for all the tie-lines are set to 0.02+j0.07 p.u. In the first control area, bus 1 is a slack bus of the power system and it is connected to the main grid. Three generators with reactive power ranges -3.75 to 3.75MVAR, 3.75 to 3.75 MVAR and -5.0 to 5.0 MVAR are attached at buses 4, 5 and 6, respectively. Reactive power generation from these DGs is utilized to minimize the power loss and voltage deviation of the system, taking the reactive power generation cost into account simultaneously. Reference voltages for 6 buses are assigned as [1.04 1.025 1.025 1.0 1.01 1.015] in an ascending order. Weight coefficients for Power loss, voltage deviation and reactive power cost are set as W1=0.2, W2=100 and W3=1.0, respectively. To analyze the real-time dynamic behavior of CS based ORPC approach, a sequence of load changing events are launched as given in Table. I.

Communication Link 1

Recovery and Reshaping the data at shown in Fig. 3.

7

6

Calculation of the control input, QGi (k+1) as shown in Fig. 4.

III Fig. 6. A power system of 18 buses with 6-bus in each control area

Transmit QGi (k+1) back to power system

END Fig. 5. Comprehensive flowchart of the proposed CS based ORPC approach

Calculation of λ is very important in the recovery algorithm of CS and performance of CS widely depends on the selection of λ. SNR ratio of the recovered signal decreases if we choose value of λ away from optimal value. After zˆ (k ) is achieved using above recovery algorithm, it is reshaped using (13) as explained in chapter III. Then at the central controller, control input QGi(k+1) is generated using improved PDIPM as explained in Fig. 4. Finally, the generated control input is sent back to the control area of power system to improve the objective function (1).

Fig. 7. 6-bus radial distribution network used as one control area TABLE.I EVENT SEQUENCES OF LOAD CHANGES ON 18-BUS POWER SYSTEM Events Time stamp Bus No Load Type Load Change Event1 5,11,17 Reactive Load 2.0×Initial 12s Event2 6,12,18 Reactive Load 0.5×Event1 24s Event3 4,10,16 Real Load 1.1×Event2 36s Event4 5,11,17 Real Load 0.95×Event3 48s

In this section, 3 different case studies of 18 bus, 129 bus and 486 bus power system are presented to test effectiveness of CS for various sizes of control area. Each power system is obtained by combing multiple control areas through tie lines to test the proposed algorithm for multi-area power systems. A. 18-bus System The proposed algorithm is applied to a multi-area power system of 18-bus network as shown in Fig. 6. It is built by

Q(p.u)

V. SIMULATION STUDIES

Fig. 8. Reactive power generation updates of 9 generators for 18-bus system

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 the network. A significant decrease in the reactive power cost and power loss is observed at 48s when real power load is reduced. According to the power system theory, the less real power demand is, the lower real power loss will be, which in turn, decreases the reactive power demand in the network.

Fig. 9.Voltage improvement of the non-PV buses for 18-bus system

Fig. 11. Updates of the individual sub-function 2.2 2

F(p.u)

1.8

Fig. 10. Error in the voltage updates with CS

To begin with, the first load change is initiated on bus 5, 11 and 17 by doubling the reactive power load at 12s. The increased load is compensated by the reactive power generation source at bus 5, 11 and 17 as shown in Fig. 8. Similarly, decrease in reactive power load and change in real power load is counteracted by their respective generators. Fig. 9 shows that none of the above changes in load has introduced significant long-term voltage deviation on any of the bus voltage and all the voltages return to their original optimal values, immediately after the load change. Buses regained their initial optimal values because controllable reactive power sources adjusted their reactive power generation to maintain the optimal bus voltages. Convergence of the proposed algorithm in case of varying loading conditions shows that CS works even when the events of load change occurs in the multi-area power system. CS has been used to transmit the network data to the central controller. Difference in the voltage updates of two cases: with and without CS is shown in Fig. 10, in which only a small error (of the order of 10-3) is observed in voltage updates during the load variation period. However, after the algorithm reaches the steady-state load condition, error in the voltage vanishes. It is inferred that the proposed CS based approach does not affect the optimal control setting of the reactive power; at the same time, decreasing the data size. Commonly, the power system operates in the steady state condition with only few dynamics in the whole day. It can be concluded that CS is an excellent means of reducing the data size of power system real-time data without affecting the control settings. Figs. 11 & 12 show the updates of each sub-function and the overall objective function, respectively. The reactive power generation increases as reactive power load rises on bus 5, 11 and 17 at 12s. In comparison, the reactive power generation decreases when reactive power load declines at 24s. A rise in real power load increases reactive power cost and power loss of

1.6 1.4 1.2 1 0

10

20

30 40 time step unit: s

50

60

Fig. 12. Minimization of the overall objective function

To analyze the effect of compression ratio (ρ) on the recovered signal, OPRC algorithm was initially implemented without CS. After that, CS based ORPC algorithm was simulated with varying compression ratios. Results from both algorithms are compared in Fig. 13 which shows the error in voltage of bus 18 of area 3 by changing the compression ratio. It is clear that higher the compression ratio, more is the error in the voltage. For large compression ratios, proposed algorithm takes more time to converge. Thus, large compression ratio makes the algorithm slow down but converges to same optimal value. It is noteworthy that in small power system, the achievable compression ratio is small.

Fig. 13. Error in V18 with the use of CS

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 B. Effect of random disturbance 1.04 1.02 1 0.98

V9

V10

V11

V12

0.96 0.94

0

10

20

30 40 time step unit: s

50

60

Fig. 14 Voltage updates of area 2 with white noise of SNR=50 on V10 2.5

states have spatial correlation in them and thus are compressible. To analyze the behavior of proposed ORPC algorithm in presence of random disturbances, reactive power generation updates are shown in Fig. 15. It shows that all generators in area 2, particularly generator on bus 10 act by changing their generation to deal with any abrupt random disturbances. To study the effect of random disturbances of higher magnitudes on the recovered data, SNR ratio was decreased from 60 to 30, as shown in Fig. 16. Fig. 16 exhibits that when the random disturbance increases in the system states, error in the recovered data rises. Similar response is observed in the reactive power generation updates as shown in Fig. 17, which shows reactive power generation output of bus 10 for different SNR. C. 129-Bus power system

2 1

1.5

44

43

1 Q10

Q11

Q12

0.5 0

0

10

20

30 40 time step unit: s

50

60

Fig. 15 Reactive power updates of area 2 with white noise of SNR=50 in V10 0.02 0.01 0 SNR=60 SNR=50 SNR=40 SNR=30

-0.01 -0.02 -0.03 0

10

20 time step unit: s

30

40

Fig. 16 Error in V10 for different SNR ratios

Fig. 17 Reactive power generation updates of bus 10 for different SNR ratios

Effect of random disturbance is studied by introducing random white noise in the voltage of bus 10 in area 2 from 25s to 33s as shown in Fig.14 where white noise with SNR=50 is added. It is observed that adding random disturbance in the voltage of one bus introduces disturbance in the neighboring buses as well. It validates our assumption that power system

III Fig. 18. A power system of 129 buses with 43-bus in each control area

CS is used to deal with the signals from the 129-bus power system as shown in Fig. 18 formed by connecting three 43-bus control areas, with the line data and bus data as given in [46]. In each area, 14 controllable reactive power generators attached on [3, 5, 7, 12, 14, 16, 18, 20, 24, 26, 27, 38, 41, 42] buses are used for ORPC. Voltages of non-PV buses and line currents are taken as the data for transmission through communication channel using CS. Real-time load changes are introduced on several buses, as shown in Table. II. Reactive power generation updates, as shown in Fig. 19, converge to their optimal values before a dynamic load change occurs at 33s, when reactive power generation increases, as exhibited in Figs. 19 & 22. This rise in reactive power load causes the voltage dip for a very short interval and it immediately recovers its optimal voltage profile by applying the proposed approach, as shown in Fig. 20. Similar behavior is observed when real power load demand is reduced by 0.6 times of the initial load. Reactive power generation cost decreases as can be observed in Fig. 22. Notice that, the fall in real power demand at 66s reduces the current flow through the network, which in turn, reduces power loss, reactive power generation requirement and hence reduces the overall objective function, as shown in Fig. 23. The result in Fig. 21 exhibits the effectiveness of the proposed CS based approach as the difference for the state variables for the two cases: with and without using CS, converges to zero within a few seconds of any dynamic change. Results in Figs. 19-22 of reactive power generation, voltage magnitudes and objective function, all are obtained using CS based control scheme, which validates the effectiveness of using compressive sensing for ORPC.

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 TABLE.II EVENT SEQUENCES OF LOAD CHANGES ON 43-BUS CONTROL AREA Events Load Load Time Bus No. Type Change Event1 [12,14,16,20,24,38,41,42,55, React. 1.6*Initial 33s 57,59,63,67,81,83,85,98,100, Load 102,106,110,124,126,128] Event2 [12,14,16,20,24,38,42,55,57, Real 0.5*Event 66s 59,63,67,81,85,98,100,102, Load 1 106,110,124,126,128]

Fig. 23. Minimization of the overall objective function

Fig. 19. Reactive power updates of 14 generators

Fig. 20.Voltage improvement of the 33 non-PV buses in area 2

D. 648- bus power system In order to test the proposed algorithm’s performance with a large scale control area, a test case on 648-bus power system, built by connecting 4 areas with 162-buses each as shown in Fig 24, is conducted here. Parameters of 162-bus system can be found in [47]. The resistance and reactance for all the tie-lines are set to 0.028 p.u and 0.096 p.u., respectively. The overall system has 1144 branches and 200 reactive power generation controllers, with 50 reactive power controllable sources in each control area. 153 buses are non-PV buses, and their state information data needs to be transmitted to the central controller for controlling the 50 reactive power sources. Reactive power loads on buses [15, 21, 22, 23, 27, 36, 67, 68, 84, 89, 94, 98, 137, 142, 148] are increased to 1.6 times of their initial values at 100s to test the system’s dynamic behavior under the proposed approach. 1

163

8 91 151

170 253 313

495 487 637

333 325 475

577

415

Fig. 24.Large power system of 648 buses with 162-bus in each control area

Fig. 21. Error in the voltage updates of the 33 non-PV buses

Update of reactive power generation using CS for data transmission is shown in Fig. 25, which shows increase in its generation at 100 s when the reactive power load is increased. Fig. 26 shows the convergence of the error to zero, immediately after the load change. Subjected to an abrupt increase in the reactive load, bus voltages deviate from their reference voltages and thus voltage deviation is increased. Similar increase is observed in power loss when load increases at 100s, as shown in Fig. 27.

Fig. 22. Updates of the individual sub-functions for 129-bus power system

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 corresponding SNR value of the recovered signal is measured as shown in Fig. 28. It shows that each power system has an optimal value of λ for which SNR value is the highest. For rest of the λs, quality of the recovered signal decreases. It also exhibits that λ changes with the size of power system. For large power systems, amount of data transfer is large and thus optimal λ will be large for high quality signal recovery.

1 0.8 0.6 0.4 0.2 0

0

50

100 time step unit: s

150

200

 

Fig. 25. Reactive power updates of 50 reactive power controllers 10-3

2 0 -2 -4 -6

F. Effect of the compression ratio on the recovered data Each of the three case studies is tested for various compression ratios and it is observed that by increasing the compression ratio, the quality of recovered signal is affected. Fig. 29 shows the curves between SNR and the compression ratio for the three case studies, which exhibits that SNR drops as the compression ratio increases. It is also shown that as the network size increases, the value of achievable compression ratio grows i.e. the values of achievable compression ratio for 6-bus, 43-bus and 162-bus control areas are 2.4, 5.445 and 9.3636, respectively. Thus, the proposed CS based approach is promising for the large system applications where data size can be reduced to 1/9 of the original size.

-8 -10

0

50

100 time step unit: s

150

200

SNR

Fig. 26. Error in the voltage updates

10

F(p.u)

8 6 4 Dv

2 0

0

50

P loss

100 time step unit: s

Fig. 29. Curves showing the relation between ρ and SNR for the three cases

CQ

150

VI. CONCLUSIONS

200

Fig. 27.Updates of three sub-objective functions

 

SNR

E. Analysis of choosing optimal λ and its effect on the recovered data

Fig. 28 Effect

of λ on the quality of the recovered signal

λ is the number of non-zeros used in the least squares problem of data recovery and it has an optimal value for each power system. If λ used in OMP problem changes from the optimal λ, the recovered signal may be of poor quality. To visualize its effect on the signal quality, various values of λ are used and

Issue of limited-bandwidth communication network in ORPC has been addressed by using CS for data transmission from dispersed control areas to the central controller. OMP based algorithm, which can significantly accelerate the recovery speed, is utilized to recover the data at the receiving end of communication channel. It is shown that data size can be reduced 1/9 times with the help of CS without deteriorating the performance of the proposed approach using CS. Three case studies are presented, ranging from small to large power systems, to validate the effectiveness of the proposed CS based approach for ORPC problem. Various load changes are introduced to different buses in the network to test the performance of the proposed approach on the real time power system. It is observed that difference in the state variables (with and without using CS) approaches to zero, immediately after the occurrence of load changes. The value of achievable compression ratio is higher for the large control areas than that of the small control areas, which indicates the promising applications of the proposed approach in the large systems. Error in the recovered data depends on the compression ratio of the measurement matrix. It also varies with number of non-zero values taken in the recovery algorithm. Recovery algorithm has an optimal number of non-zeros for each power system. Any

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 deviation from this value will cause rise in the error of the recovered data. VII. REFERENCES [1] [2] [3] [4] [5] [6]

[7] [8]

[9] [10] [11] [12] [13] [14]

[15]

[16] [17]

[18]

[19] [20] [21]

W. Zheng, W. Wu, B. Zhang, et al., "A fully distributed reactive power optimization and control method for active distribution networks," IEEE Trans. Smart Grid, vol. 7, no. 2, pp. 1021-1033, March 2016. Verizon. State of the Market: The Internet of Things 2015. http://www.verizonenterprise.com/state-of-the-market-internet-ofthings/, 2015. [Online; accessed 2016-12-22]. H. Gharavi and B. Hu, "Scalable synchrophasors communication network design and implementation for real-time distributed generation grid," IEEE Trans. Smart Grid, vol. 6, no. 5, pp. 2539-2550, Sept. 2015. S. Das and T. Singh Sidhu, “Application of compressive sampling in synchrophasor data communication in WAMS,” IEEE Trans. Ind. Inform., vol. 10, no. 1, pp. 450– 460, 2014. M. Qiu, W. Gao, M. Chen, J. W. Niu and L. Zhang, "Energy efficient security algorithm for power grid wide area monitoring system," IEEE Trans. Smart Grid, vol. 2, no. 4, pp. 715-723, Dec. 2011. S. Chen, W. Hu, C. Su, X. Zhang and Z. Chen, "Optimal reactive power and voltage control in distribution networks with distributed generators by fuzzy adaptive hybrid particle swarm optimization method," IET Gener., Transm. & Distrib., vol. 9, no. 11, pp. 1096-1103, 2015 W. Zhang, W. Liu, X. Wang, L. Liu and F. Ferrese, "Distributed multiple agent system based online optimal reactive power control for smart grids," IEEE Trans. Smart Grid, vol. 5, no. 5, pp. 2421-2431, Sept. 2014. A. Colet-Subirachs, A. Ruiz-Alvarez, O. Gomis-Bellmunt, et al., "Centralized and distributed active and reactive power control of a utility connected microgrid using IEC61850," IEEE Syst. Journal, vol. 6, no. 1, pp. 58-67, March 2012. K. Turitsyn, P. Šulc, S. Backhaus, M. Chertkov, "Use of reactive power flow for voltage stability control in radial circuit with photovoltaic generation", Proc. IEEE Power Eng. Soc. Gen. Meeting, pp. 1-6, Jul. 2010. M. E. Elkhatib, R. E. Shatshat and M. M. A. Salama, "Decentralized Reactive Power Control for Advanced Distribution Automation Systems," IEEE Trans. Smart Grid, vol. 3, no. 3, pp. 1482-1490, Sept. 2012. H. K. Nguyen, H. Mohsenian-Rad, A. Khodaei and Z. Han, "Decentralized reactive power compensation using nash bargaining solution," IEEE Trans. Smart Grid, 2016, to be published. Y. Xu, W. Zhang and W. Liu, "Distributed dynamic programming-based approach for economic dispatch in smart grids," IEEE Trans. Ind. Inform., vol. 11, no. 1, pp. 166-175, Feb. 2015. S. Bolognani, R. Carli, G. Cavraro and S. Zampieri, "Distributed reactive power feedback control for voltage regulation and loss minimization," IEEE Trans. Automatic Control, vol. 60, no. 4, pp. 966-981, April 2015. A. Maknouninejad and Z. Qu, "Realizing Unified Microgrid Voltage Profile and Loss Minimization: A Cooperative Distributed Optimization and Control Approach," IEEE Trans. Smart Grid, vol. 5, no. 4, pp. 16211630, July 2014. M. S. Lacher, J. Nonnenmacher and E. W. Biersack, "Performance comparison of centralized versus distributed error recovery for reliable multicast," IEEE/ACM Trans. Networking, vol. 8, no. 2, pp. 224-238, Apr 2000. Qiang Wan, W. Zhang, Y. Xu and I. Khan, "Distributed control for energy management in a microgrid," 2016 IEEE/PES Transm. & Dist. Conf. & Expo. (T&D), pp. 1-5, 2016. N. Jumpasri, K. Pinsuntia, K. Woranetsuttikul, T. Nilsakorn and W. Khan-ngern, "Comparison of distributed and Centralized control for partial shading in PV parallel based on Particle Swarm Optimization Algorithm," 2014 Int. Electr. Engineering Congress (IEECON), Chonburi, pp. 1-4, 2014. O. Alizadeh Mousavi and R. Cherkaoui, "Maximum voltage stability margin problem with complementarity constraints for multi-area power systems," IEEE Trans. Power Systems, vol. 29, no. 6, pp. 2993-3002, Nov. 2014. . O. Alizadeh Mousavi, R. Cherkaoui, Literature survey on fundamental issues of voltage and reactive power control, 2011, [online] Available: http://mars.ethz.ch/en/researh-and-publications/publications.html. Yannick Phulpin., “Coordination of reactive power scheduling in a multiarea power system operated by independent utilities”, Electric power. Georgia Institute of Technology, 2009. English. Mousavi, Omid Alizadeh, and Rachid Cherkaoui. "On the inter-area optimal voltage and reactive power control." Int. J. Electr. Power Energy Syst., no. 52, pp. 1-13, 2013.

[22] M. Alizadeh. “On Voltage and Frequency Control in Multi-Area Power Systems Security,” Diss. ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE, 2014. [23] A. Zhang H. Li F. Liu H. Yang "A coordinated voltage/reactive power control method for multi-TSO power systems" Int. J. Electr. Power Energy Syst. vol. 43 pp. 20-28 2012. [24] Y. Phulpin M. Begovic D. Ernst "Coordination of voltage control in a power system operated by multiple transmission utilities" Proc. IREP Symp. 2010. [25] Y. Phulpin, “Coordination of reactive power scheduling in a multi-area power system operated by independent utilities,” Georgia Inst. Technol. 2009. [26] Y. Phulpin M. Begovic M. Petit D. Ernst "A fair method for centralized optimization of multi TSO power system" Int. J. Electr. Power Energy Syst. vol. 31 pp. 482-488 2009. [27] Y. Phulpin M. Begovic M. Petit D. Ernst "On the fairness of centralized decision making strategies in multi TSO power systems" Proc. Power Syst. Comput. Conf., 2008. [28] Y. Phulpin M. Begovic M. Petit "External Network modeling for MVAr scheduling in multi area power systems", Proc. PowerTech 2007. [29] Y. Phulpin M. Begovic M. Petit J. B. Heyberger D. Ernst "Evaluation of network equivalents for voltage optimization in multi-area power systems" IEEE Trans. Power Syst. vol. 24 no. 2 pp. 729-743 May 2009 [30] E. J. Candes and M. B. Wakin, "An Introduction To Compressive Sampling," IEEE Signal Processing Magazine, vol. 25, no. 2, pp. 21-30, March 2008. [31] J. Haupt, W. U. Bajwa, M. Rabbat and R. Nowak, "Compressed Sensing for Networked Data," IEEE Signal Processing Magazine, vol. 25, no. 2, pp. 92-101, March 2008. [32] Z. Zou, Y. Bao, H. Li, B. F. Spencer and J. Ou, "Embedding Compressive Sensing-Based Data Loss Recovery Algorithm Into Wireless Smart Sensors for Structural Health Monitoring," IEEE Sensors Journal, vol. 15, no. 2, pp. 797-808, Feb. 2015. [33] Y. Yu, F. Han, Y. Bao and J. Ou, "A Study on Data Loss Compensation of WiFi-Based Wireless Sensor Networks for Structural Health Monitoring," IEEE Sensors Journal, vol. 16, no. 10, pp. 3811-3818, May15, 2016. [34] A. Ozdemir, "Loss-constrained economic power dispatch," IEE Proc. Gener., Transm. and Distrib., vol. 143, no. 5, pp. 387-392, Sep. 1996 [35] H. Nezamabadi and M. Setayesh Nazar, "Arbitrage strategy of virtual power plants in energy, spinning reserve and reactive power markets," in IET Gener., Transm. & Distrib. vol. 10, no. 3, pp. 750-763, 2 18 2016. [36] D. Jiang, L. Nie, Z. Lv and H. Song, "Spatio-Temporal Kronecker Compressive Sensing for Traffic Matrix Recovery," IEEE Access, vol. 4, no. , pp. 3046-3053, 2016. [37] D. Needell and J. A. Tropp, “CoSaMP: Iterative signal recovery from incomplete and inaccurate samples. Applied and computational harmonic,” Appl. Comput. Harmonic Anal., vol. 26, no. 3, pp. 301–321, 2009. [38] J. Tropp and A. C. Gilbert, “Signal recovery from random measurements via orthogonal matching pursuit,” IEEE Trans. Inf. Theory, vol. 53, no. 12, pp. 4655–4666, 2007. [39] M. Leinonen, M. Codreanu and M. Juntti, "Sequential Compressed Sensing With Progressive Signal Reconstruction in Wireless Sensor Networks," IEEE Trans. Wireless Commun., vol. 14, no. 3, pp. 1622-1635, March 2015. [40] Jizhong Zhu, "Reactive Power Optimization,” in Optimization of Power System Operation, 1, Wiley-IEEE Press, pp.409-454, 2009. [41] F. Meng; B. Chowdhury; M. Chamana, "Three-phase optimal power flow for market-based control and optimization of distributed generations, "IEEE Trans. on Smart Grid, to be published. [42] C. Florin, G. Mevludin, E. Damien and W. Louis, “Interior-point based algorithms for the solution of optimal power flow problems,” Electric Power Syst. Research, vol. 77, no. 5–6, pp. 508-517, April 2007. [43] Yu, Z., and D. Pu. “A New Nonmonotone Line Search Technique for Unconstrained Optimization.” Journal of Computational and Applied Mathematics 219 (1): 134–144, 2008. [44] Liwei Wang, Xinggao Liu & Zeyin Zhang, “An efficient interior-point algorithm with new non-monotone line search filter method for nonlinear constrained programming”, Engineering Optimization, 49:2, 290-310, 2017. [45] Pati, Yagyensh Chandra, Ramin Rezaiifar, and P. S. Krishnaprasad. "Orthogonal matching pursuit: Recursive function approximation with applications to wavelet decomposition." Signals, Systems and Computers, 1993. 1993 Conference Record of The Twenty-Seventh Asilomar Conference on. IEEE, 1993.

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.

This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/ACCESS.2017.2752178, IEEE Access

IEEE ACCESS-2017-04301 [46] C. Ahn and H. Peng, "Decentralized voltage control to minimize distribution power loss of microgrids," IEEE Trans. Smart Grid, vol.4, no.3, pp.1297-1304, Sept. 2013 [47] Power system test case archive, 1993, [Online], Available: http://www.ee.washington.edu/research/pstca/.

Irfan Khan (S’ 16) received a B.Sc. in Electrical Engineering from University of Engineering & Technology, Lahore, Pakistan, in June 2009 and a M.Sc. in Electrical Power Engineering from the University of Greenwich, London, UK, in 2011. He is currently working towards his Ph.D. in Electrical and Computer Engineering from Carnegie Mellon University, Pittsburgh, PA, USA. His research interests include control and optimization of smart energy networks. Yinliang Xu (M’ 13) received the B.S. and M.S. degrees in Control Science and Engineering from Harbin Institute of Technology, China, in 2007 and 2009, respectively, and the Ph.D. degree in Electrical and Computer Engineering from New Mexico State University, Las Cruces, NM, USA, in 2013. He is now an Assistant Professor with TsinghuaBerkeley Shenzhen Institute (TBSI), Lab of Smart Grid and Renewable Energy, Shenzhen, Guangdong, China. He is also an adjunct faculty with Department of Electrical and Computer Engineering at Carnegie Mellon University, Pittsburgh, PA, USA. Prior to joining TBSI, he was with School of Electronics and Information Technology, Sun Yat-Sen University during Nov. 2013 to Aug. 2017. From Nov. 2013 to Aug. 2014, he was a Visiting Scholar with the Department of Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, PA, USA His research interests include distributed control and optimization of power systems, renewable energy integration and microgrid modeling and control.

Soummya Kar (M’ 10) received B. Tech. in electronics and electrical communication engineering from the Indian Institute of Technology, Kharagpur, India, in May 2005 and Ph.D. in electrical and computer engineering from Carnegie Mellon University, Pittsburgh, PA, in 2010. From June 2010 to May 2011, he was with the Electrical Engineering Department, Princeton University, Princeton, NJ, USA, as a Postdoctoral Research Associate. He is currently an Associate Professor of Electrical and Computer Engineering at Carnegie Mellon University, Pittsburgh, PA, USA. His research interests include decision-making in large-scale networked systems, stochastic systems, multiagent systems and data science, with applications to cyber-physical systems and smart energy systems. He has published extensively in these topics with more than 140 articles in journals and conference proceedings and holds multiple patents. Recent recognition of his work includes the 2016 O. Hugo Schuck Best Paper Award from the American Automatic Control Council, the 2016 Dean's Early Career Fellowship from CIT, Carnegie Mellon, and the 2011 A.G. Milnes Award for best PhD thesis in Electrical and Computer Engineering, Carnegie Mellon University.

Hongbin Sun (SM’ 12) received his double B.S. degrees from Tsinghua University in 1992, the Ph.D from Dept. of E.E., Tsinghua University in 1996. From 2007.9 to 2008.9, he was a visiting professor with School of EECS at the Washington State University in Pullman. He is now ChangJiang Scholar Chair professor in Dept. of E.E. and the director of energy management and control research center, Tsinghua University. He is currently IEEE senior member and IET Fellow. He also serves as the editor of the IEEE TSG, associate editor of IET RPG, and member of the Editorial Board of four international journals and several Chinese journals. His technical areas include electric power system operation and control with specific interests on the Energy Management System, System-wide Automatic Voltage Control, and Energy System Integration. He has published more than 400 peer-reviewed papers, within which over 60 are IEEE and IET journal papers, and 4 books. He has been authorized 5 US Patents of Invention and more than 100 Chinese Patents of Invention.

2169-3536 (c) 2017 IEEE. Translations and content mining are permitted for academic research only. Personal use is also permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.