• <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.

    欧美日韩亚洲国产一区二区在线观看| 欧美激情高清一区二区三区| 99riav亚洲国产免费| 在线看a的网站| 这个男人来自地球电影免费观看| 成人黄色视频免费在线看| 精品久久蜜臀av无| 淫妇啪啪啪对白视频| 久久久久久久久中文| 咕卡用的链子| 久久人妻av系列| 一区二区三区国产精品乱码| 成年版毛片免费区| 黑人巨大精品欧美一区二区mp4| 99国产精品一区二区蜜桃av| 亚洲五月色婷婷综合| 最新美女视频免费是黄的| 国产亚洲欧美98| 国产成人啪精品午夜网站| 中文字幕最新亚洲高清| 可以免费在线观看a视频的电影网站| 人人妻人人澡人人看| 中文欧美无线码| 99riav亚洲国产免费| 69精品国产乱码久久久| 国产国语露脸激情在线看| 精品国产一区二区久久| av片东京热男人的天堂| 精品欧美一区二区三区在线| svipshipincom国产片| 人人澡人人妻人| 中文字幕高清在线视频| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片 | 午夜影院日韩av| 亚洲av电影在线进入| 99热国产这里只有精品6| xxxhd国产人妻xxx| 成人国产一区最新在线观看| 国产一区二区激情短视频| 亚洲色图 男人天堂 中文字幕| 成人影院久久| 在线观看www视频免费| 亚洲精品一二三| 18美女黄网站色大片免费观看| 精品久久久久久成人av| 国产亚洲精品久久久久久毛片| 两个人免费观看高清视频| 69av精品久久久久久| 亚洲一区高清亚洲精品| 咕卡用的链子| 亚洲五月色婷婷综合| 99久久久亚洲精品蜜臀av| 女人高潮潮喷娇喘18禁视频| 一a级毛片在线观看| 黄色女人牲交| 亚洲熟妇中文字幕五十中出 | 怎么达到女性高潮| 好男人电影高清在线观看| 国产午夜精品久久久久久| 欧美日韩福利视频一区二区| av电影中文网址| 国产男靠女视频免费网站| 国产免费现黄频在线看| 老司机午夜福利在线观看视频| 又黄又粗又硬又大视频| 国产成人系列免费观看| 最近最新中文字幕大全免费视频| 午夜成年电影在线免费观看| 国产无遮挡羞羞视频在线观看| 免费在线观看完整版高清| 大香蕉久久成人网| 看片在线看免费视频| www.自偷自拍.com| 亚洲七黄色美女视频| 色综合欧美亚洲国产小说| 美女国产高潮福利片在线看| 88av欧美| 男人操女人黄网站| 脱女人内裤的视频| 男女午夜视频在线观看| 亚洲精华国产精华精| 男女高潮啪啪啪动态图| 免费在线观看影片大全网站| 极品教师在线免费播放| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女 | 国产日韩一区二区三区精品不卡| 老熟妇仑乱视频hdxx| 亚洲av日韩精品久久久久久密| av福利片在线| 色婷婷久久久亚洲欧美| 男人舔女人下体高潮全视频| 18禁裸乳无遮挡免费网站照片 | 制服诱惑二区| 亚洲国产欧美网| 精品国产亚洲在线| 久久久精品欧美日韩精品| av片东京热男人的天堂| 亚洲精品av麻豆狂野| av国产精品久久久久影院| 国产精品一区二区免费欧美| 97人妻天天添夜夜摸| 日韩中文字幕欧美一区二区| avwww免费| 亚洲精品国产精品久久久不卡| 亚洲九九香蕉| 大香蕉久久成人网| 中文字幕色久视频| 亚洲av片天天在线观看| 国产高清视频在线播放一区| 两个人免费观看高清视频| 成在线人永久免费视频| 成人手机av| 99久久精品国产亚洲精品| 午夜免费鲁丝| 国产aⅴ精品一区二区三区波| 高潮久久久久久久久久久不卡| 人人澡人人妻人| 久久精品国产综合久久久| 在线观看www视频免费| 免费少妇av软件| 成人免费观看视频高清| 亚洲人成电影观看| 国产亚洲欧美精品永久| 99久久99久久久精品蜜桃| 老司机靠b影院| 国产成人系列免费观看| 欧美黄色片欧美黄色片| 搡老熟女国产l中国老女人| 亚洲午夜精品一区,二区,三区| 色哟哟哟哟哟哟| 亚洲国产看品久久| 在线天堂中文资源库| 久久天躁狠狠躁夜夜2o2o| e午夜精品久久久久久久| 夜夜躁狠狠躁天天躁| 一边摸一边抽搐一进一出视频| 波多野结衣高清无吗| 日本黄色日本黄色录像| 亚洲精华国产精华精| 99久久综合精品五月天人人| 色播在线永久视频| 老汉色∧v一级毛片| 精品国产美女av久久久久小说| 欧美日韩乱码在线| 国产激情久久老熟女| 天堂俺去俺来也www色官网| 国产黄a三级三级三级人| 欧美黄色淫秽网站| 亚洲一区二区三区色噜噜 | 久久久精品国产亚洲av高清涩受| 欧美日韩瑟瑟在线播放| 国产一区二区在线av高清观看| 国产亚洲欧美98| 亚洲精品一卡2卡三卡4卡5卡| 99国产精品免费福利视频| 日本vs欧美在线观看视频| 色婷婷久久久亚洲欧美| 五月开心婷婷网| 在线av久久热| 色哟哟哟哟哟哟| 俄罗斯特黄特色一大片| 黄片播放在线免费| 午夜免费成人在线视频| 久久精品亚洲av国产电影网| 国产成人精品久久二区二区免费| 久热这里只有精品99| 成人黄色视频免费在线看| 美国免费a级毛片| 婷婷六月久久综合丁香| 极品人妻少妇av视频| 亚洲国产精品合色在线| 国产免费现黄频在线看| 国产日韩一区二区三区精品不卡| 亚洲国产精品sss在线观看 | 777久久人妻少妇嫩草av网站| 欧美午夜高清在线| 中国美女看黄片| 国产精品一区二区免费欧美| 精品人妻1区二区| 免费人成视频x8x8入口观看| 搡老乐熟女国产| 看黄色毛片网站| 在线视频色国产色| 欧美激情高清一区二区三区| 成熟少妇高潮喷水视频| 亚洲专区中文字幕在线| 一夜夜www| 欧美大码av| 欧美+亚洲+日韩+国产| 大型av网站在线播放| 亚洲情色 制服丝袜| 国产精品影院久久| 国产精品偷伦视频观看了| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 色婷婷av一区二区三区视频| 精品国产亚洲在线| 国产精品久久久久久人妻精品电影| 亚洲情色 制服丝袜| 午夜影院日韩av| 超色免费av| 国产熟女xx| 久久中文看片网| 国产亚洲欧美98| 夜夜看夜夜爽夜夜摸 | 亚洲在线自拍视频| 免费少妇av软件| 一区在线观看完整版| 国产精品乱码一区二三区的特点 | 久久香蕉激情| 国产精品亚洲一级av第二区| 黄频高清免费视频| 日日爽夜夜爽网站| 男女下面进入的视频免费午夜 | 校园春色视频在线观看| www.自偷自拍.com| 国产av又大| 中文字幕精品免费在线观看视频| 亚洲成av片中文字幕在线观看| 性欧美人与动物交配| 欧美日韩精品网址| 欧美激情高清一区二区三区| 久久 成人 亚洲| 久久亚洲精品不卡| 成人国语在线视频| 男人舔女人的私密视频| 欧美乱妇无乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜免费观看网址| 色综合欧美亚洲国产小说| ponron亚洲| 夜夜爽天天搞| 亚洲自偷自拍图片 自拍| 男女床上黄色一级片免费看| 日本三级黄在线观看| www.自偷自拍.com| 黄片大片在线免费观看| 亚洲专区字幕在线| 淫秽高清视频在线观看| av天堂久久9| 亚洲成人精品中文字幕电影 | 久久国产精品影院| 久久久久久久午夜电影 | 免费高清视频大片| 18禁黄网站禁片午夜丰满| 久久国产亚洲av麻豆专区| 无限看片的www在线观看| 校园春色视频在线观看| 亚洲一区中文字幕在线| 久久久水蜜桃国产精品网| 日本wwww免费看| 怎么达到女性高潮| 91成人精品电影| 国产有黄有色有爽视频| www国产在线视频色| 91在线观看av| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久男人| 中出人妻视频一区二区| 亚洲色图av天堂| 日韩欧美一区二区三区在线观看| av超薄肉色丝袜交足视频| 欧美大码av| 精品一区二区三区视频在线观看免费 | 黄网站色视频无遮挡免费观看| 99在线人妻在线中文字幕| 亚洲人成伊人成综合网2020| 欧美 亚洲 国产 日韩一| 欧美黑人精品巨大| 久久午夜亚洲精品久久| 一级毛片高清免费大全| 婷婷六月久久综合丁香| 成人亚洲精品av一区二区 | 欧美另类亚洲清纯唯美| 亚洲成人免费av在线播放| 别揉我奶头~嗯~啊~动态视频| 欧美精品亚洲一区二区| 满18在线观看网站| 一本大道久久a久久精品| 五月开心婷婷网| 国产精品一区二区免费欧美| 亚洲激情在线av| 老司机福利观看| 精品人妻1区二区| 真人一进一出gif抽搐免费| 手机成人av网站| 国产av精品麻豆| 丁香六月欧美| 亚洲熟妇熟女久久| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 久久久久久亚洲精品国产蜜桃av| 乱人伦中国视频| tocl精华| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av高清一级| 色播在线永久视频| 90打野战视频偷拍视频| 真人做人爱边吃奶动态| 久久香蕉国产精品| 亚洲色图综合在线观看| 国产欧美日韩精品亚洲av| 国产片内射在线| 欧美在线黄色| 一级片'在线观看视频| 国产亚洲精品第一综合不卡| 国产亚洲欧美98| 欧美+亚洲+日韩+国产| tocl精华| 12—13女人毛片做爰片一| 桃红色精品国产亚洲av| 日韩中文字幕欧美一区二区| 国产精品偷伦视频观看了| av免费在线观看网站| 中文字幕人妻丝袜制服| 亚洲自偷自拍图片 自拍| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 老司机在亚洲福利影院| 亚洲精品一区av在线观看| 久久精品人人爽人人爽视色| 视频区图区小说| 中文字幕精品免费在线观看视频| 久久人人精品亚洲av| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| 亚洲国产毛片av蜜桃av| 丝袜人妻中文字幕| 99国产精品99久久久久| 国产又色又爽无遮挡免费看| 国产精品成人在线| 韩国av一区二区三区四区| 精品一品国产午夜福利视频| 欧美色视频一区免费| 亚洲专区字幕在线| 欧美老熟妇乱子伦牲交| 极品人妻少妇av视频| 黄色成人免费大全| 午夜福利在线免费观看网站| 夫妻午夜视频| 韩国精品一区二区三区| 精品无人区乱码1区二区| 老熟妇仑乱视频hdxx| 后天国语完整版免费观看| 久久香蕉精品热| 这个男人来自地球电影免费观看| 在线观看66精品国产| 日韩人妻精品一区2区三区| 一边摸一边抽搐一进一出视频| 男女之事视频高清在线观看| 岛国在线观看网站| 大型av网站在线播放| 欧美日韩瑟瑟在线播放| 母亲3免费完整高清在线观看| 久久精品亚洲熟妇少妇任你| 久久99一区二区三区| 欧美性长视频在线观看| 免费人成视频x8x8入口观看| 一级a爱视频在线免费观看| cao死你这个sao货| 亚洲专区国产一区二区| 日韩大码丰满熟妇| 动漫黄色视频在线观看| 国产熟女午夜一区二区三区| 999久久久国产精品视频| 国产欧美日韩一区二区精品| 999久久久国产精品视频| 欧美日韩一级在线毛片| 自线自在国产av| 国产欧美日韩综合在线一区二区| 欧美+亚洲+日韩+国产| 亚洲人成电影免费在线| 女人被躁到高潮嗷嗷叫费观| 一级作爱视频免费观看| 韩国精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 人人妻,人人澡人人爽秒播| 精品久久久久久电影网| 亚洲国产毛片av蜜桃av| 亚洲男人的天堂狠狠| 欧美黑人精品巨大| 最新美女视频免费是黄的| 黄片小视频在线播放| 久久人妻av系列| 99精品欧美一区二区三区四区| 黄色丝袜av网址大全| 亚洲视频免费观看视频| 我的亚洲天堂| videosex国产| 一本大道久久a久久精品| 一夜夜www| 美女国产高潮福利片在线看| www.www免费av| 久久精品91蜜桃| 久久九九热精品免费| 精品一区二区三区视频在线观看免费 | 婷婷丁香在线五月| 日韩欧美一区二区三区在线观看| 亚洲情色 制服丝袜| 一二三四社区在线视频社区8| 国产精品电影一区二区三区| 久久久精品欧美日韩精品| 久久精品亚洲精品国产色婷小说| 亚洲黑人精品在线| 亚洲成人免费电影在线观看| 99精品久久久久人妻精品| 91字幕亚洲| 免费不卡黄色视频| 午夜福利欧美成人| 91大片在线观看| 日韩欧美免费精品| 欧美日韩精品网址| 欧美精品一区二区免费开放| 欧美亚洲日本最大视频资源| 久久精品成人免费网站| 亚洲自偷自拍图片 自拍| 久久中文看片网| 首页视频小说图片口味搜索| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站 | 婷婷六月久久综合丁香| 日本vs欧美在线观看视频| 一级黄色大片毛片| 久久久水蜜桃国产精品网| a级毛片在线看网站| 午夜91福利影院| 19禁男女啪啪无遮挡网站| 国产高清国产精品国产三级| 久久国产乱子伦精品免费另类| 久9热在线精品视频| 水蜜桃什么品种好| 18禁美女被吸乳视频| netflix在线观看网站| 国产精品一区二区三区四区久久 | 久久人人精品亚洲av| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 久久亚洲真实| 精品久久久精品久久久| 亚洲av第一区精品v没综合| 国产精品1区2区在线观看.| 成人手机av| 免费在线观看影片大全网站| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 脱女人内裤的视频| 精品国产乱码久久久久久男人| 9热在线视频观看99| a在线观看视频网站| 欧美日韩亚洲高清精品| 亚洲成国产人片在线观看| 日韩欧美一区视频在线观看| 高潮久久久久久久久久久不卡| 国产精品爽爽va在线观看网站 | 国产精品久久电影中文字幕| 岛国视频午夜一区免费看| 国产成人免费无遮挡视频| 叶爱在线成人免费视频播放| 成年版毛片免费区| 一区在线观看完整版| 国产av在哪里看| 成熟少妇高潮喷水视频| 99精品欧美一区二区三区四区| 男男h啪啪无遮挡| 欧美激情高清一区二区三区| 大型av网站在线播放| 色老头精品视频在线观看| 亚洲国产欧美一区二区综合| 麻豆久久精品国产亚洲av | 午夜福利一区二区在线看| 黄色a级毛片大全视频| 露出奶头的视频| 国产精品国产高清国产av| 黄色丝袜av网址大全| 婷婷六月久久综合丁香| 亚洲精品av麻豆狂野| 免费女性裸体啪啪无遮挡网站| 一级黄色大片毛片| 国产国语露脸激情在线看| 亚洲国产精品sss在线观看 | 超碰成人久久| 亚洲精品粉嫩美女一区| 久久伊人香网站| 精品人妻在线不人妻| 欧美日韩亚洲国产一区二区在线观看| 天天影视国产精品| 91麻豆av在线| 国产主播在线观看一区二区| 首页视频小说图片口味搜索| 婷婷精品国产亚洲av在线| 国产欧美日韩综合在线一区二区| 国产精品久久久av美女十八| 人人澡人人妻人| 俄罗斯特黄特色一大片| 国产精品亚洲一级av第二区| 中文字幕人妻熟女乱码| 免费少妇av软件| 9色porny在线观看| 叶爱在线成人免费视频播放| 黄色成人免费大全| 国产黄色免费在线视频| 亚洲成人免费电影在线观看| bbb黄色大片| 欧美在线黄色| svipshipincom国产片| 久久精品国产亚洲av香蕉五月| 久久午夜综合久久蜜桃| 国产成人精品在线电影| 他把我摸到了高潮在线观看| 久久这里只有精品19| 一a级毛片在线观看| 一级毛片女人18水好多| 18禁美女被吸乳视频| 动漫黄色视频在线观看| 性欧美人与动物交配| 国产精品免费一区二区三区在线| netflix在线观看网站| 午夜福利在线观看吧| 自线自在国产av| 午夜成年电影在线免费观看| 久久久久久久午夜电影 | 亚洲精品久久成人aⅴ小说| 免费观看精品视频网站| 亚洲一码二码三码区别大吗| 女人被躁到高潮嗷嗷叫费观| 国产亚洲欧美精品永久| 国产一区二区激情短视频| 亚洲一区中文字幕在线| 国产欧美日韩综合在线一区二区| 欧美黑人精品巨大| 亚洲中文日韩欧美视频| 国产精品 国内视频| 精品电影一区二区在线| 久久伊人香网站| 成人亚洲精品一区在线观看| 久久国产精品影院| 99久久精品国产亚洲精品| 亚洲精品一二三| 91在线观看av| 在线观看免费日韩欧美大片| 成人三级黄色视频| 一级毛片女人18水好多| 日韩欧美免费精品| 亚洲中文字幕日韩| 国产单亲对白刺激| 久久精品影院6| 女性被躁到高潮视频| 亚洲一卡2卡3卡4卡5卡精品中文| av欧美777| 中文字幕高清在线视频| 99riav亚洲国产免费| 亚洲精品国产精品久久久不卡| 亚洲色图av天堂| 久久久精品国产亚洲av高清涩受| 精品国产美女av久久久久小说| 韩国av一区二区三区四区| 国产精华一区二区三区| 欧美av亚洲av综合av国产av| 多毛熟女@视频| 中文字幕人妻丝袜一区二区| 91精品国产国语对白视频| 久久久国产精品麻豆| 精品高清国产在线一区| 欧美激情高清一区二区三区| 又紧又爽又黄一区二区| tocl精华| 国产男靠女视频免费网站| 在线av久久热| 啦啦啦免费观看视频1| 男女高潮啪啪啪动态图| 亚洲成av片中文字幕在线观看| 久久久久久久久中文| 十八禁人妻一区二区| 在线视频色国产色| 国产激情欧美一区二区| 免费看十八禁软件| 亚洲欧美日韩无卡精品| 精品国产乱子伦一区二区三区| 免费人成视频x8x8入口观看| 亚洲精品在线观看二区| 亚洲在线自拍视频| 交换朋友夫妻互换小说| 精品人妻1区二区| 我的亚洲天堂| 啪啪无遮挡十八禁网站| 女人高潮潮喷娇喘18禁视频| 夜夜躁狠狠躁天天躁| 交换朋友夫妻互换小说| av中文乱码字幕在线| 黄片播放在线免费| 亚洲黑人精品在线| av中文乱码字幕在线| 亚洲欧美日韩高清在线视频| 亚洲国产欧美一区二区综合| 国产三级在线视频| 脱女人内裤的视频| 午夜福利在线观看吧| 亚洲国产毛片av蜜桃av| 真人做人爱边吃奶动态| 黄片大片在线免费观看| 男人舔女人下体高潮全视频| 亚洲午夜理论影院| 欧美日韩视频精品一区| 男女做爰动态图高潮gif福利片 | 少妇 在线观看| 中出人妻视频一区二区| 久久精品人人爽人人爽视色| 日本黄色视频三级网站网址| 美女扒开内裤让男人捅视频| 黄色视频,在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品国产色婷婷电影| 亚洲七黄色美女视频|