• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Super Resolution Perception for Improving Data Completeness in Smart Grid State Estimation

    2020-10-20 08:12:52GoqiLingGuolongLiuJunhuZhoYnliLiuJinjinGuGungzhongSunZhoyngDong
    Engineering 2020年7期

    Goqi Ling,Guolong Liu,Junhu Zho,c,*,Ynli Liu*,Jinjin Gu,Gungzhong Sun,Zhoyng Dong

    a School of Science and Engineering, Chinese University of Hong Kong, Shenzhen, Shenzhen 518172, China

    b School of Computer Science and Technology, University of Science and Technology of China, Hefei 230027, China

    c Shenzhen Institute of Artificial Intelligence and Robotics for Society, Shenzhen 518172, China

    d School of Electrical and Information Engineering, Tianjin University, Tianjin 300072, China

    e School of Electrical Engineering and Telecommunications, University of New South Wales, Sydney, NSW 2052, Australia

    Keywords:State estimation Low-frequency data High-frequency data Super resolution perception Data completeness

    A B S T R A C T The smart grid is an evolving critical infrastructure, which combines renewable energy and the most advanced information and communication technologies to provide more economic and secure power supply services.To cope with the intermittency of ever-increasing renewable energy and ensure the security of the smart grid,state estimation,which serves as a basic tool for understanding the true states of a smart grid, should be performed with high frequency. More complete system state data are needed to support high-frequency state estimation.The data completeness problem for smart grid state estimation is therefore studied in this paper. The problem of improving data completeness by recovering highfrequency data from low-frequency data is formulated as a super resolution perception (SRP) problem in this paper. A novel machine-learning-based SRP approach is thereafter proposed. The proposed method, namely the Super Resolution Perception Net for State Estimation (SRPNSE), consists of three steps: feature extraction, information completion, and data reconstruction. Case studies have demonstrated the effectiveness and value of the proposed SRPNSE approach in recovering high-frequency data from low-frequency data for the state estimation.

    1. Introduction

    The smart grid is a critical infrastructure that continuously provides a secure and economic electricity supply to modern society.State estimation in the smart grid plays a vital role in system monitoring and control, which helps the system operator to perceive the system’s operational states and make accurate control decisions [1]. At present, the smart grid development poses new requirements for state estimation. On the one hand, with the fast-increasing penetration of renewable energy sources and new power appliances—such as wind power, solar power, and electric vehicles—much greater uncertainties are being introduced into the smart grid[2].In order to mitigate the adverse impacts of intermittent renewable energy, the system operator needs to perceive the system’s operational states more frequently and shorten the dispatch interval, which requires the support of high-frequency state estimation. On the other hand, with the rapid improvement of computational capability,big data technologies are being widely applied to discover hidden knowledge by analyzing system measurement data [3]. State estimation, which directly monitors voltage,current, and real and reactive power values, is a basic tool for the perception of system working states.Therefore,more frequent state estimation results are helpful for discovering more hidden knowledge, which is beneficial for improved system security and efficiency.

    Performing high-frequency state estimation poses new practical challenges and difficulties.First,the state estimation monitors system states based on thousands of meters deployed on nodes,generators, transmission lines, and so forth. In the existing supervisory control and data acquisition (SCADA) system, the majority of meters in a smart grid are traditional sensors that were deployed years ago. Traditional meters collect measurements with a relatively low frequency, such as every few minutes [4,5]. Second,some high-frequency sampling meters, such as phasor measurement units(PMUs),can gather measurement data at a much higher frequency than traditional meters [6–8]. However, a single measurement is insufficient for state estimation; at any specific time point, the input of the state estimation must be a vector of measurements.PMUs are extremely expensive;it is financially impractical to replace all traditional meters with PMUs [7]. Third, even if all traditional meters can be replaced by PMUs, the capability of system state perception is very restricted by the capacity of communication channels. The high-frequency sampling data collected by PMUs cannot be transferred completely to the control center,but are usually stored at remote locations temporarily [5]. Therefore, PMUs’ high-frequency data cannot be truly utilized. In summary, due to technical restrictions, performing high-frequency state estimation is still a difficult task.

    Furthermore, both the data collected by traditional meters and those collected by PMUs may be lost or manipulated due to communication faults or cyber-attacks [9], such as false data injection attacks(FDIAs)[10,11],cyber topology attacks[12–14],and cyber–physical attacks [15,16]. Real-world cyber-attack events, such as the Ukraine blackout in December 2015 [17] and the Venezuelan power outage in March 2019 [18], indicate that cyber-attacks can result in serious consequences. In the literature, many studies have been conducted to detect abnormal data caused by cyber attackers. For example, Ashok et al. [19] proposed the detection of anomalies by describing a statistical characterization of the variation between SCADA-based state estimates and predicted system states based on load forecast information, generation schedules,and the synchro-phasor data. Esmalifalak et al. [20] utilized a distributed support vector machine (SVM)-based method to distinguish between attacked data and normal measurement data.However, so far there is very limited research on effective approaches to recover these manipulated data. Missing and tampered data also introduce significant challenges to accurate state estimation.

    To tackle the challenges discussed above, novel methods are needed to support the high-frequency perception of system operational states based on the existing metering infrastructure. In this paper, we consider this problem as a data completeness improvement problem.We consider the original measurements received in the control center as incomplete data (i.e., low-frequency data),and the data that the system operator is expected to use to make high-frequency decisions as complete data (i.e., high-frequency data). Therefore, this problem is equivalent to the question of how to recover high-frequency data from low-frequency data in order to achieve data completeness improvement. Approaches may vary according to different research fields and different data quality attributes; it is of utmost importance to explore appropriate approaches to improve data completeness for smart grid state estimation.

    In this paper,we achieve this goal by applying super resolution(SR) technology. SR is a technology that can recover highresolution data from low-resolution data from both temporal and spatial perspectives[21,22].Currently,the most effective methods to obtain high-resolution data from low-resolution data are mainly based on interpolation, reconstruction, and machine learning,respectively[23].The machine-learning-based SR,which attempts to obtain a priori mapping between low-resolution and highresolution image blocks by given samples, has become a hot topic in recent years due to its good performance [24]. For example,Dong et al.[25]proposed a novel deep learning approach—namely,a super-resolution convolutional neural network (SRCNN)—to learn the end-to-end mapping between low-resolution and highresolution images; this approach shows superior performance in restoring high-resolution images. Wang et al. [26] proposed an enhanced super-resolution generative adversarial network (ESRGAN) to recover images, which achieves consistently better perceptual quality than other SR methods.

    The motivations for this paper are obvious.High-frequency perception of the system’s operation status is of great importance for the development of the smart grid. However, traditional meters’low sampling frequency,PMUs’high investment cost, the capacity limitation on communication channels,and the abnormal statuses caused by communication mistakes or cyber-attacks present obstacles and challenges in practical situations.The purpose of this paper,therefore,is to develop an super resolution perception(SRP)approach to improve the data completeness for smart grid state estimation.This paper makes the following two key contributions:

    (1)We are among the first to study the data completeness problem for smart grid state estimation in this paper.

    (2)We are among the first to propose an effective SRP approach to recover high-frequency data from low-frequency data for state estimation. This paper has proved the effectiveness and value of the proposed Super Resolution Perception Net for State Estimation(SRPNSE) approach.

    This paper is organized as follows. Section 2 provides some background information on smart grid state estimation and its data completeness problem. Section 3 describes the SRP problem and presents the network structure and solving framework of a novel deep learning approach:SRPNSE.Section 4 demonstrates the effectiveness of the proposed approach by simulations.Finally,Section 5 provides conclusions and future work discussions.

    2. The data completeness problem for state estimation

    In this section, we give a brief introduction to smart grid state estimation and its data completeness problem.

    2.1. Smart grid state estimation

    In the smart grid, the key for the perception of system operational states is to obtain system measurements—that is,the vector of the steady-state voltage (magnitude and angle) at each bus of the network.Once the voltage information is grasped,all other system state variables can be readily calculated using power flow equations [4,5]. However, not all nodes’ voltage magnitudes and phase angles can be easily telemetered. Other information, such as the real and reactive power flows of some transmission lines,or some real and reactive power injections, need to be monitored so as to satisfy different control purposes (e.g., providing alerts for emergency situations). In addition, not all telemetered data are reliable, due to measurement errors caused by disturbances or cyber-attacks.State estimation is a tool to estimate system state variables from all available system measurements.Therefore,state estimation in modern power systems plays a vital role in the online monitoring, analysis, and control of smart grids.

    Usually, state estimation is a module embedded in the energy management system(EMS)of smart grids.In addition to necessary communication networks, the overall state-estimation-related modules contain three main components: sensors, a state estimator, and a bad data detector. Sensors measure system states, such as bus voltage magnitude, real and reactive power injections, real and reactive power flow, and so forth, at a certain sampling frequency [1]. The state estimator utilizes all collected data to estimate system state variables in order to obtain a snapshot of the power system in the steady state.The bad data detector thereafter detects and eliminates obvious errors in measurements.

    The state estimation process can be considered as a generalized power flow calculation. As shown in the following model, z represents the vector of measurement data with size[m,1];x represents the vector of system state variables with size[n,1],and m >n;h(x)represents the functional relationship between measurement values and system state variables; and e represents the vector of noises with size [m, 1].

    The goal is to calculate the ^x vector from the z vector, making the estimated^x as close as possible to the actual x based on certain estimate criteria. A widely used method for estimating ^x is the maximum likelihood (ML) estimation:

    where R is the diagonal matrix of the measurement error variances.

    2.2. Data completeness for state estimation

    Data quality is a critical problem for solving any industrial problems. Without high-quality data, the performance of both analytical and data-driven models will be seriously compromised. When based on poor-quality data,model outcomes will encounter unpredictable deviations that can cause substantial economic losses and security risks. Most of the existing data quality studies focus on database systems [27]. Data completeness is an important attribute for assessing the quality of a dataset [27]. Completeness is usually defined as whether there are any gaps in the data from what was actually collected and what was expected to be collected—that is, whether there are missing data, damaged data, or manipulated data [28].

    In this paper, the low-frequency measurements actually received by the control center are considered as incomplete. The control center perceives a big picture of the system state by processing the measurement vector z. Let the discrete sequence zk(·)represent the value of the kth row of vector z; that is, zk(·) is the measurement value of the kth meter. The corresponding mathematical expression can be represented by a time series, as shown in Eq. (4).

    where fk(t) represents the function of the measurement of the kth meter at time t.

    Due to the low sampling frequency of traditional meters, the capacity limitation on communication channels, and communication mistakes or even cyber-attacks,some part of zk(t)will be missing. Similarly, let us define another time series~zk(t ), in which its time labels t differ from zk(t), as shown in Eq. (5).

    For example, the time series zk(t) and~zk(t ) are shown in Figs.1(a) and (b), respectively. The aggregated series is then shown in Fig. 1(c).

    Data completeness for state estimation in this paper mainly refers to the completeness of the measurement vector z in the temporal dimension.The data completeness improvement is to generate more data based on available information so that it is as close as possible to the actual measurements—that is, to recover the missing time series~zk(t ).

    3. SRP model and the solving method

    In this section, the SRP problem is first proposed. Second, the problem of data completeness improvement for state estimation is formulated. Third, due to the powerful feature extraction capability of deep neural networks,a deep learning approach—namely,SRPNSE—is proposed.Finally,the optimization algorithms for estimating the model parameters of SRPNSE are introduced.

    3.1. SRP modeling

    Regardless of the size of the system, the input of state estimation is a vector of the measurements collected by many meters at a specific time point. Given vector z, the data completeness improvement problem aims to recover each meter’s missing data.We can solve this problem one meter after another. Recovering the missing data of each specific meter is equivalent to a SRP problem. The SRP problem can be expressed as follows:

    where Lβrepresents the low-frequency data actually collected by meters;H represents the original high-frequency data;e represents the vector of noises, where the noises here are caused by the meters; and ↓βrepresents the down-sampling function, in which β is the down-sampling factor.For example,suppose H is a measurement vector sampled at 60 Hz; if β is equal to 10, then Lβwill be a measurement vector based on H with a sampling rate of 6 Hz.

    In this paper, the SRP problem is formulated as a maximum a posteriori (MAP) estimation problem. Based on the MAP estimation, the goal is to estimate an Hβto maximize the posterior probability, shown as follows:

    Fig. 1. The relationship of discrete sequences for meter k. Time series (a) zk(t) and(b) ~zk (t ); (c) aggregated series Zk(t).

    where p(Lβ, ↓βH) is the likelihood function and p(↓βH) is the prior probability of ↓βH. In this situation, the prior probability mainly relates to the original sampling error caused by meters. Since the state estimation model assumes that the measurement noises follow the normal distribution,in this paper,we assume that the prior probability follows Gaussian distribution as well.

    3.2. Problem formulation

    The problem of improving data completeness for state estimation is equivalent to the problem of generating high-frequency data from low-frequency data. The high-frequency data is considered as complete data because it can recover the information lost in the incomplete data. The SRP problem formulation is shown in Fig.2.Given a set of original data D, two down-sampling data sets Lα and Hβare generated with down-sampling factors α and β. The objective of SRP is to take the lower frequency data Lα as the input and generate a set of estimated higher frequency data H^βthat is as close to the real down-sampling data Hβas possible. ?2-norm can be used to measure the difference between Hβand H^β.

    3.3. SRPNSE framework

    The network structure of the proposed SRPNSE method is shown in Fig. 3. The SRPNSE network directly uses the lowfrequency data as the input, and then performs information enhancement to output the estimated high-frequency data.SRPNSE implements data quality improvement through the following three steps: feature extraction, information completion, and data reconstruction.

    Fig. 2. Diagram of the SRP problem.

    In the feature extraction stage,three one-dimensional(1D)convolutional layers [29] are used to extract features from lowfrequency historical data. After obtaining the abstract features,the second part—the information completion stage—will supply higher resolution of features,based on the knowledge learned from the relationship between the low-frequency data and the highfrequency data. The information completion stage consists of several SRPNSE blocks that are implemented by a residual structure.The residual structure consists of a big global residual connection and a number of local residual blocks. The global residual connection forces the network to learn the missing information rather than form the signal itself,and the local residual blocks provide the possibility to train deeper networks [29]. For better performance, this research used a total of 22 local residual blocks in the information completion stage. Fig. 4. shows the structure of the residual block used in this paper,where g represents the output of the previous layer; the rectified linear unit(ReLU)function[29]is used as the activation function;and identity represents the identity mapping. Then, the higher resolution features that contain more details of the system patterns are used to reconstruct the targeted high-frequency data in the third part—the data reconstruction stage—which is implemented by three 1D convolutional layers. In this part, the feature vectors outputted by the information completion are integrated into β sub-sequences with length l. The β sub-sequences are then rearranged into the reconstructed high-frequency sequence with length β × l.

    When training the proposed SRPNSE network, the mean squared error(MSE)is chosen as the loss function,which is shown as follows:

    where H^iβand Hiβare the ith data of H^βand Hβ,respectively.N is the size of the vector.

    In this paper, we consider both the mean absolute percentage error (MAPE) [30] and the signal-to-noise ratio (SNR) [31] as evaluation metrics. The MAPE represents the degree of average absolute error compared with the actual value. A higher MAPE value means a larger difference between the actual value and the test one. The MAPE can be calculated as follows:

    In the signal processing field,the SNR represents the ratio of the average power of the signal to the average power of the noise. A higher SNR value indicates a smaller noise that the test value contains. The SNR can be calculated as follows:

    Fig. 3. Network structure of the proposed SRPNSE framework.

    Fig. 4. Network structure of the SRPNSE block i, where i = [1, 2, ..., K].

    3.4. Optimization method

    In this paper, a novel neural network, SRPNSE, is proposed to solve the SRP problem. As a deep neural network, it consists of multiple residual blocks and has strong expressive power. That is, given a specific nonlinear function, the neural network can asymptotically approximate this function by appropriately adjusting its parameters. However, due to the existence of activation functions and multiple hidden layers in the network, the function underlying SRPNSE is highly nonlinear and nonconvex.The nonlinearity and nonconvexity of SRPNSE makes its parameter optimization problem extremely difficult. In this paper, we investigate efficient optimization algorithms for estimating the parameters of SRPNSE.

    In the existing literature, gradient-based algorithms (as shown in Algorithm 1) remain the mainstream methods for the parameter estimation of deep neural networks. Typical examples are the batch gradient descent (BGD), stochastic gradient descent (SGD),and mini-batch gradient descent (Mini-BGD) algorithms [32]. In these algorithms, parameter updating is performed according to the following formulas:

    where W represents the weights; b is the bias; αlis the learning rate;and dW and db are the partial derivatives of cost function with resp

    ect to variables W and b, respectively.

    The main difference between these gradient-based algorithms is that when updating parameters in one iteration, BGD trains the network based on all batches of training data and SGD stochastically selects only one batch for training, while Mini-BGD selects only a portion of the batches.The biggest drawback of BGD is that the convergence speed is slow, especially when the number of batches is large, because it solves the gradient by calculating all batches. On the other hand, since SGD’s and Mini-BGD’s updating of the gradient direction is dependent on one or only a few data batches,their convergence trajectories are very unstable,resulting in continuous oscillation and local optima. Considering the drawbacks of traditional methods, based on Eqs. (13) and (14), many other methods such as Momentum [32], the adaptive gradient algorithm (Adagrad) [32], the root mean square prop (RMSProp)[32,33],and the adaptive moment estimation(ADAM)[32,34]have been proposed for improving the training process.In this paper,we will investigate the effectiveness of two new algorithms—RMSProp and ADAM—for estimating the parameters of SRPNSE.

    3.4.1. Root mean square prop

    One main disadvantage of the dominant gradient descent methods is that the learning rate αlis a fixed value. Choosing a proper learning rate can be difficult. If it is too small, the convergence speed will be very slow,whereas if it is too large,the loss function will oscillate or even deviate significantly from the minimum value.RMSProp(as shown in Algorithm 2)is a variant of the dominant gradient descent methods that overcomes this shortcoming.Compared with Eqs. (13) and (14), RMSProp achieves an excellent adaptation of the learning rate by adding a moving average of the squared gradient over adjacent mini-batches [34]. For each iteration, as shown in Eqs. (17) and (18), the given learning rate αlis dynamically adjusted by the root mean square. The root mean square is actually the root of the exponential moving averages of squared past gradients. According to Eqs. (15) and (16), RMSProp limits the reliance of the update to only the past few gradients

    [34]. The root mean square in RMSProp aims to balance the oscillation amplitude of different dimensions. When the parameter space is relatively flat, the partial derivation is small; then, the exponential moving average is small and the learning rate speeds up as a result. When the parameter space is relatively steep, the partial derivation is large; then, the exponential moving average is large and the learning rate slows down as a result.

    ?

    Algorithm 2. The pseudo-code for the RMSProp algorithm.Require: f(W, b) Objective function with parameters W, b Require:αl Learning rate Require: β2 Exponentially decaying average of squared gradients Require: W0, b0, S0(dW), S0(db) Initial parameters For each epoch t, while {not converged} do:1. Calculate the gradient of the objective function with respect to the current parameter:dW, db = ?f(W, b)2. Calculate biased the second moment estimate from historical gradients:SdW := β2SdW + (1-β2)(dW)2, Sdb := β2Sdb + (1-β2)(db)2 3. Calculate the descending gradient at the current moment:ηW =-αl ■ +ε, ηb = -αl dW■ +ε 4. Update based on the descending gradient:W :=W +ηW, b:=b +ηb End Return w, b■■■■■■db SdW■■■■■Sdb

    where(dW)2and(db)2are the square of the gradient;β2represents the exponential decay rate,which is usually set as 0.9 or 0.999;SdWand Sdbrepresent the exponential moving averages of squared past gradients; and ε is a very small number, say 10-8, to prevent the denominator from being 0.

    3.4.2. Adaptive moment estimation

    The other disadvantage of the dominant gradient descent methods is that the current gradient is the only factor to determine the descent direction. Once the current gradient is pointing in the opposite direction of the previous gradient, the loss function will oscillate or even deviate from the minimum value. Momentum

    [32] is a variant of the dominant gradient descent methods that overcomes this shortcoming. Compared with Eqs. (13) and (14),Momentum achieves the stability for faster learning by adding the accumulation of the exponential moving average of past gradients and then moving in that direction [32].

    The ADAM algorithm (as shown in Algorithm 3) combines the ideas of both the Momentum and RMSProp algorithms. As shown in Eqs. (19) and (20), the exponentially decaying average of the gradients is calculated, which is the Momentum and is called the‘‘first-order moment estimation” in ADAM; as shown in Eqs. (21)and(22),the exponentially decaying average of the squared gradients is calculated,which is the RMSProp and is called the‘‘secondorder moment estimation” in ADAM. In addition, as shown in Eqs.(23) and (24), ADAM computes a bias-corrected first-order moment estimate and second-order moment estimate to offset the deviation caused by the initialized zero vectors. As shown in Eqs. (25) and (26), ADAM not only updates the descent direction by an exponentially decaying average of gradients,but also divides the learning rate by an exponentially decaying average of squared gradients. As a result, faster convergence and reduced oscillation are gained [34].

    Algorithm 3. The pseudo-code for the ADAM algorithm.Require: f(W, b) Objective function with parameters W, b Require:αl Learning rate Require: β1, β2 Exponentially decaying average of squared gradients Require: W0, b0, S0(dW), S0(db),υ0(dW),υ0(db) Initial parameters For each epoch t, while {not converged} do:1. Calculate the gradient of the objective function with respect to the current parameter:dW, db = ?f(W, b)2. Calculate biased first and second moment estimate from historical gradients:υdW :=β1υdW + 1-β1( )dW, υdb :=β1υdb +(1-β1)db SdW:=β2SdW +1-β2( ) db( )2 3. Update the bias-corrected of first and second order moment estimate:( ) dW( )2,Sdb:=β2Sdb +1-β2 υc dW :=υdW 1-βt 1,υc db := υdb 1-βt 1,Sc dW :=SdW 1-βt 2,Sc db :=Sdb 1-βt 2 4. Update based on the descending gradient:W :=W -αl υc d■ +ε,b:=b-αl ■ +ε End Return W, b■■■■■■W υc d b Scd W■■■■Scd b

    4. Case studies

    In this section,we conduct case studies based on a 9-bus system

    [35]. As shown in Fig. 5, three meters are deployed on nodes 5, 7,and 9 to record the real power loads; three meters are deployed on nodes 1, 2, and 3 to measure the real and reactive power outputs of the generators; five meters are deployed on the‘‘from-end” of branches 1–4, 5–6, 6–7, 8–2, and 9–4 (e.g., the‘‘from-end”of branch 1–4 is bus 1)to measure the real and reactive power flows; and four meters are deployed on the ‘‘to-end” of branches 4–5, 3–6, 7–8 and 8–9 (e.g., the ‘‘to-end” of branch 4–5 is bus 5) to measure the real and reactive power flows.

    In order to better simulate the state estimation scenarios, the input values of state estimation in this paper are assumed to be some of the results of optimal power flow(OPF)calculations based on measured loads. As shown in Fig. 5, the first three meters, providing three real powers of loads, are considered as the input values of OPF;the remaining 12 meters,providing overall 24 real and reactive powers,are the input values of the state estimation for the system; and the values of the 12 meters come from the results of OPF.

    It is assumed that on each load node, 1 MW electricity is set to supply approximately 200 households.Each household contains 11 types of appliances, such as air conditioners, heaters, washing machines, microwaves, and so on; each appliance’s waveform comes from the plug load appliance identification dataset (PLAID)

    Fig. 5. Topology structure of the 9-bus system. G: generator.

    [36]. The PLAID samples 11 different types of appliances at 3 × 104Hz, which is down-sampled to 100 Hz in this paper.In this paper, we simulate 900, 1000, and 1250 households, with 100 times magnification for the load on nodes 5, 7, and 9,respectively. The super resolution perception state estimation dataset (SRPSED) was designed for testing the proposed SRPNSE; it has been released and can be found at https://www.zhaojunhua.org/SRP/SRPSE/dataset/. In generating this dataset, the user behavior of a normal office worker was adopted for each household; the user behavior has also been released with the dataset. The SRPSED contains a total of 60 d of high-frequency data at a frequency of 100 Hz, in which the first 45 d are used for training and validation,and the last 15 d are used for testing.Since the data completeness improvement problem for every single meter is independent in this paper, the proposed SRPNSE approach can be applied to larger systems in a similar manner.

    4.1. Performance of SRPNSE in state estimation

    Interpolation[37]is a popular way to fill vacancies and replace wrong data in many areas. In this paper, linear interpolation and cubic interpolation are applied as comparisons to the proposed SRPNSE approach. State estimation calculates the 9-bus system’s state—that is,the voltage magnitude and angle—based on the given measurements. In this paper, we execute the state estimation based on the real down-sampled data,SRPNSE data,linear interpolation data, and cubic interpolation data, respectively. Four case studies with a total of 16 scenarios were conducted. For each scenario, we calculated the MAPE and SNR of the voltage magnitude and angle.

    4.1.1. Performance evaluated with MAPE

    The MAPE values for both voltage magnitude and angle for all scenarios under SRPNSE, linear interpolation, and cubic interpolation are shown in Appendix A Tables S1 and S2, respectively (also see Figs.6 and 7).It should be noticed that each subfigure contains four scenarios. For example, Fig. 6(a) represents scenarios of recovering data from 1/60, 1/300, 1/600, and 1/900 Hz with η =5, respectively.

    4.1.2. Performance evaluated with SNR

    The SNR values for both voltage magnitude and angle for all scenarios under SRPNSE, linear interpolation, and cubic interpolation are shown in Appendix A Tables S3 and S4, respectively(also see Figs. 8 and 9).

    Fig. 6. MAPE comparisons among SRPNSE, linear interpolation, and cubic interpolation about the voltage magnitude. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    Fig. 7. MAPE comparisons among SRPNSE, linear interpolation, and cubic interpolation about the voltage angle. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    Table 1 The down-sampling factors and SR factors used in this paper.

    4.2. Performances of the SRPNSE on load nodes

    4.2.1. Performance evaluated with MAPE

    The MAPE values in load nodes 5,7,and 9 for all scenarios under SRPNSE, linear interpolation,and cubic interpolation are shown in Appendix A Tables S5, S6, and S7, respectively. Here, we take the scenarios on node 5 as a representative (Fig. 10).

    4.2.2. Performance evaluated with SNR

    The SNR values in load nodes 5, 7,and 9 for all scenarios under SRPNSE, linear interpolation,and cubic interpolation are shown in Appendix A Tables S8, S9, and S10, respectively. Here, we take the scenarios on node 5 as a representative (Fig. 11).

    Fig. 8. SNR comparisons among SRPNSE, linear interpolation, and cubic interpolation about the voltage magnitude. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    Fig. 9. SNR comparisons among SRPNSE, linear interpolation, and cubic interpolation about the voltage angle. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    Fig. 10. MAPE comparisons among SRPNSE, linear interpolation, and cubic interpolation on load node 5. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    Fig. 11. SNR comparisons among SRPNSE, linear interpolation, and cubic interpolation on load node 5. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    4.3. Visualized comparison of state estimation

    The 9-bus system has a total of three generators, nine nodes,and nine branches. Four case studies with the state estimation for a total of 16 scenarios were conducted; here, we randomly selected case 4, node 3, and branch 5–6 to show visualized comparisons of the state estimation results in more detail, as a representative. Specifically, we randomly chose a time period to show the fluctuation of voltage magnitude on node 3, voltage angle on node 3, power flow on branch 5–6, and generator output on node 3 with η =5,η =10,η =50,and η =100,respectively.In each figure,the true data,SRPNSE data,linear interpolation data,and cubic interpolation data are visualized for comparison.

    4.3.1. Case 4 with η =5: Voltage magnitude on node 3

    In this scenario α =90 000 and η =5;that is,the low-frequency data received in the control center is 1/900 Hz, and the goal is to restore 1/180 Hz data from the 1/900 Hz data.Based on the recovered data, state estimation is executed. Here, the voltage magnitude on node 3 is drawn. A comparison of the voltage magnitude after state estimation between the real down-sampled data and the estimated ones by SRPNSE, linear interpolation, and cubic interpolation is shown in Fig. 12.

    4.3.2. Case 4 with η =10: Voltage angle on node 3

    In this scenario, α =90 000 and η =10; that is, the lowfrequency data received in the control center is 1/900 Hz, and the goal is to restore 1/90 Hz data from the 1/900 Hz data. Based on the recovered data, state estimation is executed. Here, the voltage angle on node 3 is drawn. A comparison of the voltage angle after state estimation between the real down-sampled data and the estimated ones by SRPNSE, linear interpolation, and cubic interpolation is shown in Fig. 13.

    Fig.12. Voltage magnitude(Mag)on node 3 after state estimation,recovering data from 1/900 to 1/180 Hz.(a)True data;(b)SRPNSE;(c)linear interpolation;(d)cubic interpolation. p.u.: per unit.

    4.3.3. Case 4 with η =50: Power flow on branch 5–6

    In this scenario, α =90 000 and η =50; that is, the lowfrequency data received in the control center is 1/900 Hz, and the goal is to restore 1/18 Hz data from the 1/900 Hz data. Based on the recovered data, state estimation is executed. Here, the power flow on branch 5–6 is drawn. A comparison of the power flow on branch 5–6 after state estimation between the real down-sampled data and the estimated ones by SRPNSE, linear interpolation, and cubic interpolation is shown in Fig. 14.

    4.3.4. Case 4 with η =100: Generator output on node 3

    In this scenario, α =90 000 and η =100; that is, the lowfrequency data received in the control center is 1/900 Hz, and the goal is to restore 1/9 Hz data from the 1/900 Hz data. Based on the recovered data,state estimation is executed.Here,the generator output on node 3 is drawn. A comparison of the generator output on node 3 after state estimation between the real downsampled data and the estimated ones by SRPNSE, linear interpolation, and cubic interpolation is shown in Fig. 15.

    4.4. Visualized comparison of load nodes

    Four case studies with a total of 16 scenarios were conducted.Here, we randomly selected case 1 on load node 5 and case 4 on load node 9 and discussed them in greater detail. In each figure,the true data, SRPNSE data, linear interpolation data, and cubic interpolation data are visualized for comparison.

    Fig.13. Voltage angle(Ang) on node 3 after state estimation,recovering data from 1/900 to 1/90 Hz. (a) True data; (b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 14. Power flow on branch 5–6 after state estimation, recovering data from 1/900 to 1/18 Hz. (a) True data; (b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    4.4.1. Load node 5: Data completeness improvement from 1/60 Hz

    In this case, α =6000 and η =5, η =10, η =50, and η =100;that is, the low-frequency data received in the control center is 1/60 Hz, and the goal is to restore 1/12, 1/6, 5/6, and 5/3 Hz data from the 1/60 Hz data, respectively. Here, we randomly chose a time period. Comparisons between the real 1/12, 1/6, 5/6, and 5/3 Hz down-sampled data and the estimated ones by SRPNSE, linear interpolation,and cubic interpolation are shown in Figs.16–19.

    4.4.2. Load node 9: Data completeness improvement from 1/900 Hz

    In this case,α =90 000 and η =5,η =10,η =50,and η =100;that is, the low-frequency data received in the control center is 1/900 Hz, and the goal is to restore 1/180, 1/90, 1/18, and 1/9 Hz data from the 1/900 Hz data, respectively. Here, we randomly chose a time period. Comparisons between the real 1/180, 1/90,1/18, and 1/9 Hz down-sampled data and the estimated ones by SRPNSE, linear interpolation, and cubic interpolation are shown in Figs. 20–23.

    4.5. Comparison of the SGD, RMSProp, and ADAM algorithms for solving the SRPNSE framework

    Tables S11 and S12 in Appendix A provide comparisons of the MAPE and SNR values using the ADAM and RMSProp algorithms compared with the SGD algorithm for solving the proposed SRPNSE framework. Here, the MAPE values on load node 5 for three algorithms are shown as a representative (Fig. 24).

    Fig. 15. Generator output on node 3 after state estimation, recovering data from 1/900 to 1/9 Hz. (a) True data; (b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 16. Load measurements, recovering data from 1/60 to 1/12 Hz. (a) True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 17. Load measurements, recovering data from 1/60 to 1/6 Hz. (a) True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 18. Load measurements, recovering data from 1/60 to 5/6 Hz. (a) True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 19. Load measurements, recovering data from 1/60 to 5/3 Hz. (a) True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig.20. Load measurements,recovering data from 1/900 to 1/180 Hz.(a)True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig.21. Load measurements, recovering data from 1/900 to 1/90 Hz.(a)True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig.22. Load measurements,recovering data from 1/900 to 1/18 Hz.(a)True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    Fig. 23. Load measurements, recovering data from 1/900 to 1/9 Hz. (a) True data;(b) SRPNSE; (c) linear interpolation; (d) cubic interpolation.

    We also compared the loss function in iterations using the RMSProp and ADAM algorithms compared with the SGD algorithm.Here, the MSE value for the scenario of case 4 with η =100 is selected as a representative(Fig.25(a)).Fig.25(b)provides magnified views of the first 80 iterations and the iterations from 300 to 380, respectively.

    4.6. Result analysis

    Tables S1–S4, Figs. 6–9, and Figs. 12–15 provide comparisons focused on the state estimation results. Tables S5–S10, Figs. 10 and 11, and Figs. 16–23 provide comparisons focused on the measurements of load nodes. It should be noted that the latter comparisons are SR results of the meters’measurements,while the former comparisons illustrate the state estimation results, which are based on the SR results of the meters’ measurements.

    Fig. 24. MAPE comparisons using RMSProp and ADAM algorithm compared with SGD. (a)η = 5; (b)η = 10; (c)η = 50; (d)η = 100.

    First, it is clear from Tables S1–S4 and Tables S5–S10 that the proposed SRPNSE significantly outperforms the linear and cubic interpolation methods. This indicates that the data supplemented by the SRPNSE approach can achieve a more accurate estimate of the actual situation, and thus helps to achieve a more accurate state estimation result.More importantly,the differences between using the SRPNSE, linear interpolation, and cubic interpolation methods on the state estimation are obvious. As shown in Tables S1–S4 and Figs. 6–9, the value differences are as high as one or two orders of magnitude. This indicates that the linear and cubic interpolations are weak in recovering lost information from relatively low-frequency data, while the proposed SRPNSE approach performs well.

    Second,it is obvious from Tables S5–S7 that,no matter what the SR factor is, the MAPE values of the SRPNSE and interpolation methods keep increasing when the sampling frequency drops;and, more importantly, the lower the frequency is, the higher the MAPE difference between the SRPNSE and interpolation methods is.The reason for the small MAPE differences between the SRPNSE and interpolation methods, as in case 1, is that relatively highfrequency data already contains enough information, which helps to improve the accuracy of the interpolations.

    Third, by comparing the scenarios of different SR factors for a specific case, such as case 3 with η =5, η =10, η =50, and η =100, it is clear from Tables S5–S7 and Tables S8–S10 that a smaller SR factor will usually lead to better performance for both the SRPNSE and the interpolation methods.

    Fourth,from Fig.24 and Tables S11 and S12,it is clear that most MAPE values based on the ADAM algorithm are lower than those based on the RMSProp and SGD algorithms. As shown in Fig. 25(a), the SGD algorithm achieves a slow convergence speed and encounters continuous oscillation; as shown in Fig. 25(b),compared with RMSProp, the ADAM algorithm is slightly slow in convergence speed, but its stability is outstanding. The result is consistent with Section 3.4: The ADAM algorithm not only uses a dynamically adjusted learning rate to speed up convergence, but also uses an accumulated momentum to stay stable. Therefore,the ADAM algorithm performs better than the RMSProp and SGD algorithms in solving the proposed SRPNSE framework.

    5. Conclusions and future works

    In this article,we proposed a novel machine-learning-based SRP approach to improve data completeness for smart grid state estimation.The case studies demonstrated the effectiveness and value of the proposed approach.

    Concerning the applicability of the SRPSNE approach in a larger system, please note that the SRPNSE is an algorithm that recovers high-frequency data for a single meter.In other words,when solving the SRP problem, the SRPNSE approach is applied to recover one meter after another, without using any information from neighboring meters. Therefore, when the SRPNSE approach is applied to a larger system, it is still possible to solve each meter one by one.Although we used the 9-bus system for the case study,the load data generated by this test system is big.The training data size of the 9-bus system used in the current case study is almost ten gigabytes. Adding a larger testing system into the case study would require substantial computational resources, which would need further investment in the hardware (including GPUs and larger memory).We would therefore like to leave this as our future work.

    Fig.25. Loss function comparisons using the RMSProp and ADAM algorithms compared with the SGD algorithm for case 4 with η=100.MSE value in interations(a)[1,500];(b) [1, 80], and (c) [300, 380].

    Furthermore, we will also perform trials on relatively higher frequency data in the future, such as recovering data from 100,10, or 1 Hz. In fact, the SRPNSE approach can not only be applied in state estimation, but also in many other important modules in smart grids.The SRPNSE approach can help to improve data quality and thus overcome the obstacles and challenges caused by deployed meters, communication channels, and abnormal data intrusion. By applying the SRPNSE approach, the efficiency and security of existing industrial systems may be improved based on poor-quality data in practical situations without further investment and upgrading.

    Acknowledgements

    This work was partially supported by the Training Program of the Major Research Plan of the National Natural Science Foundation of China (91746118), partially supported by the Shenzhen Municipal Science and Technology Innovation Committee Basic Research project (JCYJ20170410172224515), partially supported by funding from Shenzhen Institute of Artificial Intelligence and Robotics for Society, and partially supported by Youth Innovation Promotion Association of Chinese Academy of Sciences.

    Compliance with ethics guidelines

    Gaoqi Liang, Guolong Liu, Junhua Zhao, Yanli Liu, Jinjin Gu,Guangzhong Sun, and Zhaoyang Dong declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2020.06.006.

    91aial.com中文字幕在线观看| 午夜福利视频精品| av国产久精品久网站免费入址| 免费观看性生交大片5| 国精品久久久久久国模美| 久久国产亚洲av麻豆专区| 大片免费播放器 马上看| 91老司机精品| 国精品久久久久久国模美| 国产精品久久久久成人av| 成年人午夜在线观看视频| 丰满饥渴人妻一区二区三| 一边摸一边做爽爽视频免费| 一二三四在线观看免费中文在| 秋霞在线观看毛片| 老司机靠b影院| 男女午夜视频在线观看| 一个人免费看片子| 久久婷婷青草| 国产成人一区二区在线| 久久精品国产亚洲av涩爱| 国产精品一区二区在线观看99| 欧美激情高清一区二区三区 | a级毛片黄视频| 超碰97精品在线观看| 狠狠精品人妻久久久久久综合| 国产淫语在线视频| 男女之事视频高清在线观看 | 亚洲,一卡二卡三卡| 最近中文字幕高清免费大全6| 午夜日韩欧美国产| 丁香六月天网| 精品人妻熟女毛片av久久网站| 国产男女超爽视频在线观看| 色吧在线观看| 男人爽女人下面视频在线观看| 国产一区有黄有色的免费视频| 女人久久www免费人成看片| 日本色播在线视频| 亚洲成人手机| 一个人免费看片子| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 精品少妇一区二区三区视频日本电影 | 日韩熟女老妇一区二区性免费视频| 男人操女人黄网站| 精品少妇内射三级| 少妇精品久久久久久久| 中文天堂在线官网| 午夜久久久在线观看| 一边摸一边抽搐一进一出视频| 亚洲第一av免费看| 婷婷色麻豆天堂久久| 免费av中文字幕在线| 水蜜桃什么品种好| 国产毛片在线视频| 成人国产麻豆网| 久久 成人 亚洲| 国产精品蜜桃在线观看| 欧美最新免费一区二区三区| 色综合欧美亚洲国产小说| 国产精品熟女久久久久浪| 国产成人系列免费观看| 色精品久久人妻99蜜桃| 激情五月婷婷亚洲| 飞空精品影院首页| 日韩熟女老妇一区二区性免费视频| 又大又黄又爽视频免费| 精品国产国语对白av| 亚洲国产精品国产精品| 欧美激情高清一区二区三区 | 欧美黄色片欧美黄色片| 波多野结衣一区麻豆| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品古装| 高清视频免费观看一区二区| 最黄视频免费看| 高清不卡的av网站| 丝袜喷水一区| 欧美在线一区亚洲| 国产精品女同一区二区软件| 黑人巨大精品欧美一区二区蜜桃| 精品国产乱码久久久久久男人| 亚洲免费av在线视频| 无限看片的www在线观看| 少妇精品久久久久久久| 欧美日韩成人在线一区二区| 91精品国产国语对白视频| 亚洲熟女毛片儿| 日本91视频免费播放| 一区二区日韩欧美中文字幕| 亚洲,欧美精品.| 波野结衣二区三区在线| 最近最新中文字幕大全免费视频 | 成人手机av| 少妇的丰满在线观看| 久久综合国产亚洲精品| 在现免费观看毛片| 如何舔出高潮| 一区福利在线观看| 亚洲成色77777| 只有这里有精品99| 日本vs欧美在线观看视频| 午夜91福利影院| 国产一区有黄有色的免费视频| 交换朋友夫妻互换小说| 亚洲 欧美一区二区三区| 嫩草影视91久久| 国产精品av久久久久免费| 爱豆传媒免费全集在线观看| 亚洲欧美精品综合一区二区三区| 成人免费观看视频高清| 一级片'在线观看视频| 国产 精品1| 久久亚洲国产成人精品v| 国产精品一国产av| 99久久综合免费| 国产亚洲精品第一综合不卡| 男的添女的下面高潮视频| 国产片内射在线| 蜜桃国产av成人99| 精品免费久久久久久久清纯 | 天堂俺去俺来也www色官网| 在线天堂中文资源库| 国产国语露脸激情在线看| 王馨瑶露胸无遮挡在线观看| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 91精品国产国语对白视频| 丁香六月天网| 丝袜人妻中文字幕| 亚洲欧美一区二区三区久久| 另类亚洲欧美激情| 精品一区二区三区四区五区乱码 | 丰满饥渴人妻一区二区三| 岛国毛片在线播放| 日韩av免费高清视频| 黄网站色视频无遮挡免费观看| 精品酒店卫生间| 最近最新中文字幕免费大全7| 国产精品熟女久久久久浪| 热99久久久久精品小说推荐| 精品一区二区三区四区五区乱码 | 老汉色∧v一级毛片| 亚洲成色77777| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆乱淫一区二区| 免费观看人在逋| svipshipincom国产片| 黄色 视频免费看| 波野结衣二区三区在线| 久久午夜综合久久蜜桃| 亚洲av成人不卡在线观看播放网 | 人体艺术视频欧美日本| 青春草亚洲视频在线观看| 中文字幕亚洲精品专区| 在线天堂最新版资源| 男女国产视频网站| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠久久av| 9色porny在线观看| 国产xxxxx性猛交| 欧美日韩成人在线一区二区| 亚洲四区av| 午夜免费男女啪啪视频观看| 菩萨蛮人人尽说江南好唐韦庄| 老司机亚洲免费影院| 国产不卡av网站在线观看| 99热国产这里只有精品6| 老司机深夜福利视频在线观看 | 美女中出高潮动态图| av在线观看视频网站免费| 亚洲,一卡二卡三卡| 精品人妻一区二区三区麻豆| 亚洲精品aⅴ在线观看| 成人影院久久| 国产一级毛片在线| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 国产精品久久久av美女十八| 亚洲一卡2卡3卡4卡5卡精品中文| 日本欧美国产在线视频| 深夜精品福利| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 男女之事视频高清在线观看 | 一级毛片我不卡| 少妇被粗大的猛进出69影院| 日日撸夜夜添| 国产精品蜜桃在线观看| 一级片免费观看大全| 久久 成人 亚洲| avwww免费| 国产熟女午夜一区二区三区| 欧美精品亚洲一区二区| 波野结衣二区三区在线| 天美传媒精品一区二区| 国产色婷婷99| 极品人妻少妇av视频| 少妇 在线观看| 亚洲中文av在线| 国产男女内射视频| 国产欧美日韩一区二区三区在线| 色精品久久人妻99蜜桃| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 99久久精品国产亚洲精品| 最近2019中文字幕mv第一页| 久久久久久久久久久免费av| 久久综合国产亚洲精品| 国产日韩欧美在线精品| 一区二区av电影网| 一二三四在线观看免费中文在| 久久精品人人爽人人爽视色| 精品一区在线观看国产| 国产日韩欧美亚洲二区| 一区在线观看完整版| 色婷婷久久久亚洲欧美| 岛国毛片在线播放| 欧美另类一区| www.自偷自拍.com| 成人手机av| 热99国产精品久久久久久7| 久久av网站| 免费不卡黄色视频| 午夜福利影视在线免费观看| 亚洲精品国产区一区二| 我要看黄色一级片免费的| 国产一区有黄有色的免费视频| 在线天堂中文资源库| 国产成人免费观看mmmm| 欧美日韩av久久| 成人午夜精彩视频在线观看| 国产精品成人在线| 啦啦啦在线观看免费高清www| 欧美日韩av久久| 欧美激情 高清一区二区三区| 黄片无遮挡物在线观看| 日韩成人av中文字幕在线观看| 黄色 视频免费看| 90打野战视频偷拍视频| 日韩制服骚丝袜av| 亚洲欧美一区二区三区久久| 亚洲成人免费av在线播放| 美女脱内裤让男人舔精品视频| 国产精品秋霞免费鲁丝片| 毛片一级片免费看久久久久| 狠狠精品人妻久久久久久综合| 久久女婷五月综合色啪小说| 女人精品久久久久毛片| 免费在线观看视频国产中文字幕亚洲 | 亚洲av电影在线进入| 免费在线观看黄色视频的| 国产1区2区3区精品| videos熟女内射| 美女脱内裤让男人舔精品视频| 国产精品亚洲av一区麻豆 | 亚洲欧美一区二区三区黑人| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 男女国产视频网站| 日韩中文字幕欧美一区二区 | 在线精品无人区一区二区三| 最新的欧美精品一区二区| 亚洲av电影在线观看一区二区三区| 热re99久久国产66热| 侵犯人妻中文字幕一二三四区| 男女午夜视频在线观看| 欧美激情 高清一区二区三区| 哪个播放器可以免费观看大片| 校园人妻丝袜中文字幕| 久久人人97超碰香蕉20202| 女人久久www免费人成看片| 欧美日韩一级在线毛片| 一二三四中文在线观看免费高清| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲欧美一区二区av| 中文字幕色久视频| 国产 精品1| 19禁男女啪啪无遮挡网站| 亚洲美女视频黄频| 9热在线视频观看99| 亚洲成av片中文字幕在线观看| 午夜av观看不卡| 综合色丁香网| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看| 日韩一卡2卡3卡4卡2021年| av卡一久久| 老熟女久久久| 亚洲国产精品一区二区三区在线| www.精华液| 亚洲欧美精品综合一区二区三区| 亚洲欧美成人精品一区二区| 18在线观看网站| 欧美久久黑人一区二区| 国产欧美日韩一区二区三区在线| 丝袜美腿诱惑在线| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 欧美老熟妇乱子伦牲交| 波多野结衣一区麻豆| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品古装| 国产一区二区三区综合在线观看| 一区二区av电影网| 9热在线视频观看99| 男女床上黄色一级片免费看| 免费在线观看黄色视频的| 女人爽到高潮嗷嗷叫在线视频| 日韩电影二区| 中文字幕人妻熟女乱码| 99久久99久久久精品蜜桃| 美女午夜性视频免费| 丝袜人妻中文字幕| av天堂久久9| 国产深夜福利视频在线观看| 人妻一区二区av| 青草久久国产| 色视频在线一区二区三区| 黄网站色视频无遮挡免费观看| 久久久精品免费免费高清| 国产福利在线免费观看视频| 可以免费在线观看a视频的电影网站 | 亚洲国产精品成人久久小说| 欧美黑人精品巨大| 街头女战士在线观看网站| 久久人人爽av亚洲精品天堂| 亚洲三区欧美一区| 免费看不卡的av| 国产极品天堂在线| 亚洲情色 制服丝袜| 美女视频免费永久观看网站| 久久鲁丝午夜福利片| 日本一区二区免费在线视频| 又黄又粗又硬又大视频| 欧美黄色片欧美黄色片| 在线观看免费午夜福利视频| 亚洲自偷自拍图片 自拍| 日韩欧美一区视频在线观看| 成人手机av| 人妻人人澡人人爽人人| 久久久欧美国产精品| 国产国语露脸激情在线看| 女性被躁到高潮视频| 黄片小视频在线播放| 看非洲黑人一级黄片| 在现免费观看毛片| 欧美黑人欧美精品刺激| 午夜老司机福利片| 老司机亚洲免费影院| 中文字幕人妻丝袜一区二区 | 久久久久久人妻| 欧美在线一区亚洲| 亚洲精品av麻豆狂野| 久久国产精品男人的天堂亚洲| 香蕉国产在线看| 国语对白做爰xxxⅹ性视频网站| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 少妇人妻精品综合一区二区| 大香蕉久久网| svipshipincom国产片| 丝瓜视频免费看黄片| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 少妇被粗大的猛进出69影院| 看非洲黑人一级黄片| a 毛片基地| 人人妻人人澡人人爽人人夜夜| 多毛熟女@视频| 自线自在国产av| 麻豆av在线久日| 99久久综合免费| 久久亚洲国产成人精品v| 久久久久视频综合| 国产精品久久久久久精品古装| 久久久久久久久久久久大奶| 天天躁日日躁夜夜躁夜夜| 久久天堂一区二区三区四区| 国产极品粉嫩免费观看在线| 一级a爱视频在线免费观看| 老汉色av国产亚洲站长工具| 午夜激情久久久久久久| 亚洲精品在线美女| 男女免费视频国产| 久久人人爽av亚洲精品天堂| 欧美精品高潮呻吟av久久| 只有这里有精品99| 国产在视频线精品| 欧美人与善性xxx| 777米奇影视久久| 超碰97精品在线观看| 国产成人精品久久久久久| 天美传媒精品一区二区| 欧美97在线视频| 又大又黄又爽视频免费| 丝袜脚勾引网站| 性高湖久久久久久久久免费观看| 一边亲一边摸免费视频| 中文字幕亚洲精品专区| 中文天堂在线官网| 大片免费播放器 马上看| 国产午夜精品一二区理论片| 亚洲国产av新网站| 777久久人妻少妇嫩草av网站| 午夜福利在线免费观看网站| 国产日韩一区二区三区精品不卡| 日日啪夜夜爽| 大香蕉久久成人网| 在线 av 中文字幕| 在线天堂中文资源库| 免费av中文字幕在线| 国产野战对白在线观看| 国产黄频视频在线观看| 无限看片的www在线观看| 久久精品国产亚洲av涩爱| 亚洲欧洲日产国产| 亚洲少妇的诱惑av| 制服诱惑二区| 曰老女人黄片| 黑丝袜美女国产一区| 久久99热这里只频精品6学生| 久久精品熟女亚洲av麻豆精品| 操出白浆在线播放| 久久人人爽人人片av| 免费女性裸体啪啪无遮挡网站| 国精品久久久久久国模美| 精品一区二区三卡| 国产精品秋霞免费鲁丝片| 国产精品一二三区在线看| 亚洲欧美一区二区三区黑人| 久久韩国三级中文字幕| 热99久久久久精品小说推荐| 十分钟在线观看高清视频www| 久久久精品区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲,一卡二卡三卡| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲精品久久午夜乱码| 日韩欧美一区视频在线观看| 国产麻豆69| 黄色视频在线播放观看不卡| 九色亚洲精品在线播放| 成年人免费黄色播放视频| 久久精品久久久久久噜噜老黄| 日日撸夜夜添| 欧美亚洲日本最大视频资源| 伊人亚洲综合成人网| 王馨瑶露胸无遮挡在线观看| 亚洲成av片中文字幕在线观看| 亚洲欧洲日产国产| 涩涩av久久男人的天堂| 99香蕉大伊视频| 下体分泌物呈黄色| 最近中文字幕2019免费版| 国产亚洲欧美精品永久| 中文字幕制服av| 婷婷色麻豆天堂久久| 在线 av 中文字幕| 最近中文字幕2019免费版| 黄色怎么调成土黄色| 亚洲人成网站在线观看播放| 久久这里只有精品19| 天堂俺去俺来也www色官网| av网站免费在线观看视频| 亚洲国产最新在线播放| 岛国毛片在线播放| 黑人欧美特级aaaaaa片| 人妻 亚洲 视频| 成人三级做爰电影| 女人高潮潮喷娇喘18禁视频| av福利片在线| 精品少妇一区二区三区视频日本电影 | 久久久精品免费免费高清| 免费不卡黄色视频| 丝袜在线中文字幕| 国产黄频视频在线观看| 亚洲精品第二区| 午夜福利网站1000一区二区三区| 最近中文字幕高清免费大全6| 久久亚洲国产成人精品v| 多毛熟女@视频| 欧美97在线视频| 亚洲成国产人片在线观看| 日韩制服丝袜自拍偷拍| 精品国产一区二区久久| 国产欧美日韩一区二区三区在线| 欧美日韩一级在线毛片| 韩国av在线不卡| 男女下面插进去视频免费观看| 咕卡用的链子| 久久久国产精品麻豆| 色视频在线一区二区三区| 国产免费又黄又爽又色| 最近最新中文字幕免费大全7| 嫩草影视91久久| 精品国产乱码久久久久久小说| 狂野欧美激情性bbbbbb| 久久99一区二区三区| 超碰成人久久| 精品酒店卫生间| av国产精品久久久久影院| a级片在线免费高清观看视频| 亚洲精品国产一区二区精华液| 搡老岳熟女国产| 亚洲av成人精品一二三区| 亚洲国产欧美日韩在线播放| 街头女战士在线观看网站| 侵犯人妻中文字幕一二三四区| 老鸭窝网址在线观看| 9191精品国产免费久久| 国产精品.久久久| 国精品久久久久久国模美| 十八禁网站网址无遮挡| 校园人妻丝袜中文字幕| 在线精品无人区一区二区三| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品 | 下体分泌物呈黄色| 满18在线观看网站| 99久国产av精品国产电影| 欧美日韩一区二区视频在线观看视频在线| 久久精品国产综合久久久| 亚洲av日韩精品久久久久久密 | e午夜精品久久久久久久| 99久久人妻综合| 美女主播在线视频| 国产欧美日韩一区二区三区在线| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三区在线| 中文字幕人妻熟女乱码| 如何舔出高潮| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 婷婷色综合大香蕉| 飞空精品影院首页| 精品国产乱码久久久久久小说| 18禁观看日本| 又大又黄又爽视频免费| 丝瓜视频免费看黄片| 亚洲国产av影院在线观看| 人人妻,人人澡人人爽秒播 | av.在线天堂| 2021少妇久久久久久久久久久| 亚洲欧美一区二区三区黑人| av不卡在线播放| 欧美xxⅹ黑人| 99久国产av精品国产电影| 大话2 男鬼变身卡| 久久久久精品性色| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 一级,二级,三级黄色视频| 一级黄片播放器| 超色免费av| 亚洲精品中文字幕在线视频| av免费观看日本| 嫩草影视91久久| 亚洲伊人久久精品综合| 伊人久久国产一区二区| 亚洲成国产人片在线观看| 丝袜喷水一区| 成人国语在线视频| 久久国产精品男人的天堂亚洲| av网站在线播放免费| 午夜免费观看性视频| 熟女少妇亚洲综合色aaa.| 成年动漫av网址| 成人影院久久| 国产在线视频一区二区| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 久久97久久精品| 国产男女超爽视频在线观看| 夜夜骑夜夜射夜夜干| 天天操日日干夜夜撸| 亚洲欧美日韩另类电影网站| 午夜日本视频在线| 亚洲欧美激情在线| 超色免费av| 女人被躁到高潮嗷嗷叫费观| 日本一区二区免费在线视频| 97在线人人人人妻| 精品久久蜜臀av无| 国产精品香港三级国产av潘金莲 | av一本久久久久| 久久av网站| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 午夜日本视频在线| 国产精品久久久久久久久免| 日韩不卡一区二区三区视频在线| 男女午夜视频在线观看| 成人三级做爰电影| 欧美老熟妇乱子伦牲交| 日日爽夜夜爽网站| 精品国产露脸久久av麻豆| 男人爽女人下面视频在线观看| 国产精品 国内视频| 亚洲av综合色区一区| 日本一区二区免费在线视频| av在线app专区| 午夜久久久在线观看| 久久天堂一区二区三区四区| 看免费av毛片| 精品亚洲乱码少妇综合久久| 日本wwww免费看| 男女无遮挡免费网站观看|