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

    Research on Evaluation of Multi-Timescale Flexibility and Energy Storage Deployment for the High-Penetration Renewable Energy of Power Systems

    2023-01-25 02:52:32HongliangWangJiahuaHuDanhuangDongCenfengWangFeixiaTangYizhengWangandChangsenFeng

    Hongliang Wang,Jiahua Hu,Danhuang Dong,Cenfeng Wang,Feixia Tang,Yizheng Wang and Changsen Feng,★

    1Economic Research Institute of State Grid,Hangzhou,310016,China

    2College of Information Engineering,Zhejiang University of Technology,Hangzhou,310023,China

    ABSTRACT With the rapid and wide deployment of renewable energy,the operations of the power system are facing greater challenges when dispatching flexible resources to keep power balance.The output power of renewable energy is uncertain,and thus flexible regulation for the power balance is highly demanded.Considering the multi-timescale output characteristics of renewable energy,a flexibility evaluation method based on multi-scale morphological decomposition and a multi-timescale energy storage deployment model based on bi-level decision-making are proposed in this paper.Through the multi-timescale decomposition algorithm on the basis of mathematical morphology,the multi-timescale components are separated to determine the flexibility requirements on different timescales.Based on the obtained flexibility requirements,a multi-timescale energy resources deployment model based on bi-level optimization is established considering the economic performance and the flexibility of system operation.This optimization model can allocate corresponding flexibility resources according to the economy,flexibility and reliability requirements of the power system,and achieve the trade-offbetween them.Finally,case studies demonstrate the effectiveness of our model and method.

    KEYWORDS Multi-timescale;morphological decomposition;flexibility evaluation;energy storage deployment

    1 Introduction

    Vigorously developing renewable energy and establishing a clean,low-carbon,safe and efficient energy system have been the theme of global energy development since the 21st century[1,2].Highpenetration renewable energy power access to grid is one of the most important and clear development trends of China’s power grid in the next decade.Renewable energy represented by wind power,photovoltaic power and hydropower will maintain a rapid development momentum and have an increasingly significant impact on the power grid.Compared with traditional generator sets,the output of wind power and photovoltaic is affected by wind speed,sunshine,weather,season and geographical location,with randomness,fluctuation and intermittence.The uncertainty of wind and photovoltaic power output will impact the power grid in many aspects,such as peak shaving and frequency modulation,dispatching operation,voltage stability and relay protection[3,4].Therefore,in order to improve the proportion of new energy and in the meantime meet the system’s requirements of economy,flexibility and reliability,a reasonable flexibility evaluation method and energy storage allocation scheme are needed to explore the flexibility demand according to the output characteristics of new energy,deploy corresponding flexibility resources and improve power flexibility.

    At present,experts and scholars have conducted extensive research on power system flexibility evaluation.In the early research,the system flexibility resource is often scored according to the capacity and output characteristics of each flexibility resource to measure the system flexibility.Although the scoring method has the advantages of simple calculation process and intuitive evaluation results,it cannot reflect the relationship between flexibility resources and flexibility requirements within the system.From the perspective of probability method,Lannoye et al.[5] established the cumulative probability distribution curve of system upward and downward flexible adjustment capacity representing the flexibility resources of the system according to the time series output sequence of each power supply,and then compares the cumulative distribution curve with the time series output curve of system net load,The probability of insufficient ramping resource probability(RRP) and insufficient ramping resource expectation (IRRE) are obtained,and the evaluation of system flexibility is completed from the perspective of probability method.Nosair et al.[6]simulated the sequential operation of each unit in the system,obtains the output of each unit including wind turbine unit and conventional unit,analyses the load shedding and wind power limitation of the system,and quantitatively characterizes the system flexibility by using the probability of insufficient system flexibility and the expectation of insufficient system flexibility.Through production simulation,this method can fully reflect the flexible regulation ability of each unit in the wind-solar access system,but its high computational complexity makes it difficult to apply the method in large-scale power grid system.Li et al.[7] considered the directional characteristics of flexibility,establishes upward flexibility index and downward flexibility index,respectively,and proposes a flexibility evaluation method considering the impact of wind power and load forecasting errors on system flexibility based on Monte Carlo simulation and economic dispatching model.

    To sum up,the existing research on flexibility evaluation mainly reflects the flexibility of power system on a specific scale,and it is difficult to quantitatively evaluate the flexibility on multiple different time scales at the same time.Moreover,the consideration of multi-timescale in the above research is based on the first-order difference of net load curve,which cannot independently reflect the differences of flexibility requirements of different time scales.Accordingly,in terms of flexibility resources,it is difficult to effectively distinguish the regulation characteristics of system controllable units on different time scales.The resulting flexibility evaluation method lacks the differentiation of regulation rate,which is not conducive to the programming of flexible power supply.Therefore,how to comprehensively consider the corresponding relationship between flexibility requirements and flexibility resources on multiple time scales,and establish an evaluation method that can reflect the multi-timescale flexibility of the system is a problem worthy of research.Reference [8] proposed a flexibility evaluation method based on multi-scale morphological decomposition,which obtains multiscale components through multi-scale morphological filtering to reflect the multi-time scale flexibility requirements of the system.The decomposition principle of multi-scale morphological filtering can be referred to reference[9].

    As an important flexible regulation power supply,energy storage has been widely studied in the related applications of new energy access power system,and the energy storage technologies represented by pumped storage,battery energy storage,flywheel energy storage and super capacitor energy storage have also been more mature.Based on the method of empirical formula,reference[10]defined the peak load regulating capacity ratio as the ratio of energy storage adjustable capacity to maximum installed capacity,generalizes the concept to the power grid,and then optimizes the allocated capacity of the pumped storage unit with the objective of optimizing the comprehensive peak regulating ability of the power grid.From the perspective of power grid managers,reference[11]studied the energy storage configuration under large-scale photovoltaic access,and establishes a bi-level planning model with the objective of optimizing the improvement effect of system voltage fluctuation and minimizing the energy storage allocation capacity,considering the power grid operation constraints and energy storage operation constraints.From the perspective of power supply side,reference [12]optimized the allocation of energy storage battery capacity with the goal of smoothing the fluctuation of photovoltaic output and the longest utilization time under the unit cost of energy storage battery.Reference [13] considered the influence of the charging and discharging times and charging and discharging depth of the energy storage battery on its service life,establishes the life attenuation model of the energy storage battery,considers the optimization of the battery life in the energy storage operation strategy,and carries out the energy storage configuration of wind power.Reference [14]analysed and quantifies the impact of energy storage on the peak shaving and valley filling capacity,power quality level and feeder power active regulation capacity of active distribution network,and constructs the optimal allocation model of energy storage aiming at the optimal impact of energy storage on these three aspects.Reference[15]respectively aimed at the optimal operation stability and economy of the system,and solves the multi-objective problem based on the bi-level decision-making model.By comparing the cycle life and allocation cost of super capacitor and battery in different scenarios,reference[16]proved that the complementary allocation of two energy storage can effectively improve the operating environment of each energy storage and improve the cycle life of energy storage.Considering the complementarity of output of pumped storage and new battery energy storage on long-time scale and short-time scale,and based on the constraint of output fluctuation rate of wind power grid connected system,reference [17] established a hybrid energy storage allocation model composed of pumped storage and battery energy storage with the goal of minimizing energy storage investment cost and operation cost.Reference [8] proposed a multi-scale energy storage allocation model based on bi-level programming,and established a hybrid energy storage allocation model composed of pumped storage,liquid flow battery and super capacitor with the goal of minimizing energy storage investment cost and operation cost.

    To sum up,at present,there are relatively mature methods for building a single energy storage deployment model.The research on participation of multiple energy storage in deployment is mainly focused on battery and super capacitor energy storage deployment in the context of small-scale wind power and photovoltaic access to micro-grid.Large-scale wind power and photovoltaic access to large power grid is an irresistible development trend.Energy storage system’s economy and operation characteristics of large power grid are different from those of the micro-grid.Therefore,the research on hybrid energy storage deployment of the power grid with large-scale wind and photovoltaic access to is a positive exploration for the future power grid.

    Under the above background,this article proposes a multi-scale evaluation strategy to evaluate the output characteristics of wind power and photovoltaic,establishes a multi-scale flexibility evaluation method based on morphologically decomposed flexibility requirements.Then according to the multitimescale characteristics of flexibility demand of power system with renewable energy,a multi-scale energy storage system framework is constructed and a multi-scale energy storage deployment model based on bi-level decision-making is established considering the economy of energy storage allocation and the flexibility of system operation.

    2 Net Load Forecasting Based on Bayesian Neural Network

    2.1 Principle of Bayesian Deep Neural Network

    Compared with other deep networks,Bayesian neural network has a special probabilistic layer structure,which can achieve a more stable prediction model according to a small amount of data.At the same time,the weight and bias of neurons in the probabilistic layer obey a certain probability distribution,which gives it the ability to describe uncertain variables.Bayesian neural network can be regarded as a conditional distribution model.X,Y and W are used to represent the neural network input,predicted value and weight.Then the principle of Bayesian reasoning can be described by Eq.(1):

    where(X,Y)is training set data,p(W)is the priori probability of W;p(Y|X,W)is the probability distribution of output Y given parameters W and X.Because the dataset is known,namelyp(Y|X)is constant,the principle of Bayesian neural network can be described by Eq.(2),that is,the posterior probabilityp(W|X,Y)is calculated by the likelihoodp(Y|X,W)of data under given parameters:

    For general neural networks,it is necessary to design the objective function or loss function,and obtain the optimal neural network parameters through optimization methods such as random gradient descent algorithm or adaptive motion estimation(Adam).Bayesian neural network carries out weighted average on multiple groups of data,that is,integrating the whole weight space.In addition,the training of general neural network needs a lot of data,so it is difficult to obtain the mean variance ofp(W|X,Y).Bayesian neural network uses a relatively simple distribution,such as normal distributionq(μw,σw),and use the method of variational inference to approximatep(W|X,Y),whereμwandσware the parameters of the neural network that need to be updated and adjusted.Therefore,the training objective of Bayesian neural network is to minimize the relative entropy of two probability distributions.Relative entropy can measure the distance between two random distributions.When two random distributions are the same,their relative entropy is zero.The specific training algorithm of Bayesian neural network is as follows:

    Figure 1:Structure of the probability layer of the Bayesian neural network

    As shown in the Fig.1,ωis the weight of neurons in the probabilistic layer,which follows the normal distribution with the expected value ofμand the variance ofρ.This is a priori distribution,andσis the a priori variance of y.

    2.2 Net Load Forecasting Based on Bayesian Deep Neural Network

    The net load distribution probability is determined by the PV output distribution and the wind turbine output distribution.Therefore,firstly,the light intensity with the most significant impact on photovoltaic output and the wind speed with the most important impact on wind turbine output are preprocessed for neural network input.Secondly,due to the time coupling characteristics of photovoltaic,wind turbine and load,the output of photovoltaic and wind turbine at the current time should also be input in the neural network.Finally,the net load is the difference between the supply power and the demand power,and the load power at time t should be added to the input.Based on this,the probability distribution of net load can be obtained in Bayesian neural network.The predicted results of net load probability distribution on a typical day are as follows:

    Figure 2:Prediction results of net load

    It can be seen from Fig.2 that the net load forecasting model based on Bayesian depth network has high accuracy.The larger prediction interval at 7:00 is due to the obvious sudden change of photovoltaic output at this time.After 17 o’clock,the photovoltaic output is basically 0,and the impact of net load prediction results comes from the wind turbine.

    3 Multi-Timescale Flexibility Evaluation

    Considering the impact of multi-timescale flexibility demand caused by multi-timescale output of wind power and photovoltaic,a power system flexibility evaluation method based on multi-scale morphology is proposed to effectively distinguish multi-timescale flexibility requirement.

    3.1 Flexibility Requirement Analysis Based on Multi-Scale Morphology

    The power fluctuation of wind power,photovoltaic power and other new energy power generation has multi-timescale characteristics.New energy output and load can both be regarded as the power disturbance source of the power grid.The net load curve formed by the superposition of new energy output and grid load shows the fluctuation characteristics on multiple time scales different from the original load curve,which puts forward new requirements for flexible regulating ability of the gridcontrolled power supply.Starting from demand analysis,this article uses mathematical morphology algorithm to decompose the fluctuation of different frequency bands of the net load curve,and extract upward and downward flexibility requirements under different fluctuation rates,so as to realize the classification,evaluation and deployment of available flexibility resources.

    Expansion operation and corrosion operation are the most basic operations of mathematical morphology.Let the input signalf (n)be a discrete function defined on the definition domainF= {0,1,···,N-1},the structural elementg(n)be a discrete function defined on the definition domainG= {0,1,···,M-1},and N ≥m,then the corrosion and expansion operations off (n)ong(n)are respectively defined as:

    where ?represents corrosion operation and ⊕represents expansion operation.

    Generally,corrosion and expansion are irreversible operations.The results of corrosion before expansion and expansion before corrosion are usually different.Therefore,new morphological transformations are produced,namely opening operation and closing operation.The opening and closing operations are expressed as

    where°represents opening operation,·represents closing operation.

    Different cascading of opening and closing operations can form a variety of filters,of which OCCO filter is more commonly used.In this article,a cluster weighted multi-scale morphology filter(WMMF)is constructed.The filtering results of the target signal using multi-scale structural elements are weighted and superimposed,and the waveform decomposition under custom divided fluctuation period is realized.The mathematical expression of WMMF filter is

    whereis the variance of difference before and after filtering under scalesi,ωsiis the weight of filtering results on each scale with respective structural elements,which is determined by the variance value of the filter noise on each scale.

    The steps of mathematical morphological decomposition algorithm based on WMMF filter mainly include:

    (1) Construct different structural elements of each time scale;

    (2) Calculate the OCCO filtering results on each time scale;

    (3) The filtering results on each time scale are weighted and superimposed to form WMMF filtering results.

    Combined with the fluctuation characteristics of wind power and photovoltaic and the regulation rate distribution of conventional power supply,the net load fluctuation can be divided into three frequency bands:fast,medium and slow,and their flexibility requirements are evaluated respectively.Therefore,two series WMMF filters are designed to decompose the input signal.The input net load time series is first separated into high-frequency fluctuation series and residual component series with medium and low frequency fluctuations by the first stage WMMF filter.After the latter passes through the second stage WMMF filter,the low frequency and intermediate frequency fluctuation sequences are obtained,respectively.The multi-scale morphological decomposition effect is shown in Fig.3.

    Figure 3:Multi-scale morphological decomposition schematics

    System flexibility evaluation includes two directions:upward adjustment and downward adjustment.After decomposing the fluctuation components of different time scales,the fluctuations of different scales correspond to different regulation rate requirements.Through waveform identification,the fluctuation curve is divided into upward and downward climbing subsets.The amplitude of the climbing section in the climbing subset represents the flexibility demand of the climbing section,and the duration represents the fluctuation period of the climbing section.The length of the fluctuation period under different fluctuation components should correspond to its specified time scale.Taking the intermediate frequency component obtained by the above filtering result as an example,the flexibility demand analysis is carried out,and the effect diagram of the flexibility demand is shown in Fig.4.

    Figure 4:Multi-scale flexibility requirements schematics

    3.2 Multi-Timescale Flexible Resource Analysis

    The size of system flexibility resources is relative with fluctuation period,operation state of each unit and output regulation characteristics of the unit itself.The same flexible resource has different output regulation characteristics under different time scales,and the output regulation characteristics of different flexible resources are also different under the same time scale.The regulation capacity models of different types of controllable units in the system under different fluctuation time scales are shown in Table 1.

    Table 1: Time scale of flexible resource participation in regulation

    Table 1 (continued)

    3.3 Flexibility Analysis Steps

    (1) The daily net load curve is forecasted according to the output curves of load,wind power and photovoltaic,and a large number of sample data are simplified into a few representative typical modes.

    (2) Based on the multi-scale morphological algorithm,the daily net load curve is decomposed in multiple time scales and divided into net load component curves under multiple fluctuation frequency bands.The upward/downward flexibility requirements of each time scale are obtained by dividing the upward climbing set and downward climbing set through waveform identification.

    (3) Analyze the upward/downward flexible resources under different fluctuation time scales according to the operation status,output regulation characteristics and fluctuation period of each controllable unit.

    (4) Through the matching analysis of flexible resources and requirements on the same time scale,guidance is provided for multi-timescale energy storage deployment afterwards.

    4 Multi-Timescale Energy Storage Deployment

    In order to meet the flexibility requirements on multiple time scales,according to the technical characteristics and application scenarios of various energy storage units,energy storage capacity is programmed for pumped storage unit,all vanadium redox flow battery on hourly,15-minute scales,respectively.A multi-scale energy storage optimal deployment method based on bi-level decisionmaking is proposed.

    4.1 Energy Storage Capacity Deployment Model

    In the power and capacity deployment of energy storage system,planning and operation are nested.In addition to the flexibility requirements of the system operation,it is also necessary to consider the economy of deployment results,including energy storage investment cost and system operation cost,so as to achieve balance in the conflict between the two.Therefore,it is necessary to establish a multi-scale energy storage capacity deployment model with double-layer nested optimization.

    The outer layer of the model is responsible for making decisions on the power and capacity of the energy storage deployment,considering the impact of the energy storage allocation scheme on the investment cost of the system,the operation cost and lack of flexibility penalty caused in the system operation,in order to optimize the power and capacity deployment of the energy storage system with the goal of optimal economy.The inner layer of the model is the flexible operation layer on multiple time scales.Based on the deployment scheme obtained by the outer layer optimization,the flexible operation of different types of units on each scale of the system is considered,and the results are output to the outer layer to realize the nested optimization of the inner and outer layers.Fig.5 is the structure diagram of energy storage deployment model.

    Figure 5:Structure diagram of energy storage deployment model

    4.1.1 Outer Layer Optimization Model

    The outer model is the planning layer of energy storage deployment.According to the system operation output from the inner optimization model,considering the energy storage investment cost,system unit operation cost and lack of flexibility penalty,the planning problem of energy storage capacity and power allocation is solved with the goal of minimizing the total flexibility investment cost of power system.

    Objective function:

    whereCis the total cost of flexibility investment,which consists of investment cost,operation cost and penalty for insufficient flexibility.Cinvis the cost of storage and gas unit investment,Copris operation cost (including maintenance costCm) andClackis the penalty for insufficient flexibility,which are transmitted from the inner layer.

    whereηp,ηe,ηmare unit capacity cost,unit power cost and unit maintenance cost of energy storage units,respectively refer to deployed rated power,rated capacity and life expectancy of energy storage under the time scale i.ηp.gasis unit power cost of gas unit,Ps.gasis deployed rated power of gas unit.

    Constraints:

    The rated power and rated capacity allocated for each time scale shall be within the upper and lower limits of energy storage output power and output capacity.

    4.1.2 Inner Layer Optimization Model

    Based on the energy storage allocation scheme on each time scale calculated by the outer model,the inner model optimizes the flexible operation of various types of units in the system on each scale.

    In the inner multi-scale optimization model,the units on the short time scale need to adjust their output on the base of long time scale operation scheme,and the unit operation scheme on hourly timescale will affect the upward and downward adjustment space of gas units on 15-minute time scale.Firstly,according to the multi-scale morphological filter,the output curve of new energy is decomposed into two different time scales:1-hour level,15-minute level.Then,the scenarios on the hourly timescale are optimized,the operation of flexible units on 15-minute timescale is adjusted on the basis of the units’operation scheme on hourly timescale.Thus,the optimal operation of multi-scale flexible power supply is realized,and the operation results of each scale are transmitted to the outer layer.

    (1)Hourly scale

    On the hourly timescale,the system operation with the participation of conventional thermal power unit and pumped storage unit is optimized in order to meet the balance of system supply and demand,improve the economy of system operation.

    Constraints of coal-fired and gas unit:

    The output of gas unit is constrained by the upper and lower limits of unit output and climbing capacity.(t)is the output of coal-fired unit at time t,is upper and lower limits of coalfired unit output,is upward and downward climbing capacity of coal-fired unit.(t)is the output of gas unit at time t.

    Constraints of energy storage unit:

    The output of the energy storage unit is constrained by its stored power and the upper and lower limits of charge and discharge.are the maximum power of discharge and charging of the energy storage unit.When(t)is greater than zero,the energy storage device discharges.When less than zero,the energy storage device is charged.In order to ensure that the energy storage has the same charge and discharge capacity in adjacent operation cycles,the conservation of charge and discharge in energy storage days is constrained.Set

    Constraints of insufficient flexibility:

    When the upward flexible regulation capacity of the system power supply is insufficient to track the downward fluctuation of new energy output,cut off the power load of the system to maintain the balance between supply and demand of the system,which results in corresponding punishment for insufficient upward flexibility.When the downward flexible adjustment ability of the system power supply is insufficient to meet the requirements of rapid load reduction,the supply and demand balance of the system is maintained by cutting off the wind and photovoltaic power,causing the corresponding punishment for insufficient downward flexibility.Load(t)is current load,Pw(t),Pp(t)are current output of wind and photovoltaic power.

    Constraint of electric power balance:

    The sum of the load output,wind power and photovoltaic power is equal to the sum of the conventional unit output,pumped storage unit output,load shedding power,wind power and photovoltaic power.

    (2)15-minute scale

    Based on the power operation scheme of hourly timescale and according to the net load fluctuation on 15-minute time scale,the regulated output of 15-minute flexibility unit including gas unit and all vanadium redox flow battery storage device is optimized.

    Constraints of gas unit:

    The output of gas unit is constrained by the upper and lower limits of unit output and climbing capacity.(t)is output of gas unit at time t,t*represents the hourly timescale optimizing time corresponding to the adjustment time t in the 15 min timescale,are the upward and downward climbing capacity of gas unit.△(t)is the regulation amount of gas unit power on 15-minute scale relative to hourly scale output.

    Constraints of energy storage unit:

    The output of the energy storage unit is constrained by its stored power,the upper and lower limits of charge and discharge.are the maximum discharge and charging power of energy storage unit.When(t)is greater than zero,the energy storage device discharges.And when it is less than zero,the energy storage device is charged.In order to ensure that the energy storage has the same charge and discharge capacity in adjacent operation cycles,the conservation of energy storage charge and discharge in a day is constrained.

    Constraints of insufficient flexibility:

    Constraint of electric power balance:

    The sum ofNetload(i)(t),wind and photovoltaic limitation is equal to the sum of gas unit power adjustment,energy storage unit output and load shedding power.

    5 Case Study and Analysis

    5.1 Flexibility Demand Analysis of a Typical Day

    The net load curve of a typical day is forecasted based on Bayesian deep neural network,and the mathematical morphology algorithm is used to decompose the net load fluctuation of different frequency bands,as shown in Fig.6.Extract the upward and downward flexibility requirements under different fluctuation rates.Fig.7 shows the resulting low-frequency requirements.

    Figure 6:Multi-scale morphological decomposition of typical daily net load curve

    On the basis of flexibility demand analysis and according to the technical characteristics of various energy storage units,the energy storage units are divided into two types of flexible regulation units acting on the system’s flexibility requirements of different time scales,namely,1-hour timescale flexible regulation units represented by pumped storage units and 15-minute time scale flexible regulation units represented by all vanadium redox flow batteries,The flexible supply-demand balance of the system on the scale of 1 h and 15 min is achieved through the optimal deployment of multi-scale energy storage.

    5.2 Energy Storage Deployment Analysis

    Although both thermal power unit and pumped storage unit can realize flexible adjustment of various time scales within a day,there are great differences between them.First,thermal power units include gas-fired unit and coal-fired unit.Coal fired thermal power unit accounts for a large proportion in China’s major power grids.Therefore,as the main power supply,it undertakes the peak shaving task of the power system.The gas unit has high flexibility and usually participates in the flexibility adjustment of medium and short time scales within hours.The coal-fired unit after flexibility transformation can participate in the flexibility adjustment of hourly time scale.It takes about 2 h for a coal-fired thermal power unit to start in hot state.The pumped storage unit starts for about tens of seconds to several minutes at a time.The hot start-up of combined cycle units in gas-fired thermal power unit takes 1 h,and the hot start-up of pure gas-fired unit takes about 30 min.It takes about 2 h for coal-fired thermal power units from grid connection to full load.After the pumped storage unit is started,it only takes a few minutes to run from zero load to full load.It takes about half an hour for gas-fired thermal power units to go from grid connection to full load.Pumped storage unit has low loss and strong adaptability,and its participation in peak shaving has more advantages than coal-fired thermal power unit.Although the proportion of pumped storage units and gas-fired thermal power units is small at present,their wider application is the development direction in the future.

    Figure 7:Intermediate and low frequency flexibility requirements of typical day

    As a high-quality clean energy,natural gas is obviously superior to coal-fired power generation in terms of carbon emission.Compared with coal-fired units,natural gas power generation has high start-up success rate,good power generation quality,high electric energy conversion rate,fast electric regulation reaction speed and low carbon emission.When the power generation cost is reduced,it may replace coal-fired units as the main peak shaving resources in our province.However,restricted by the high price of natural gas,the price of natural gas power generation is significantly higher than that of coal and other power generation types.Under the same generation hours and peak shaving demand,the regulation cost per kilowatt and average annual investment cost of natural gas power generation are the highest relative to pumped storage and electrochemical energy storage.Considering natural gas units in the hybrid energy storage configuration scheme of power system is more conducive to achieve the low-carbon goal and further improve the peak shaving capacity of the system.

    Pumped storage and electrochemical energy storage are high-quality peak shaving resources through night storage and day generation.Pumped storage technology is mature,flexible and reliable,with daily regulation capacity and full power response in minutes.It has the characteristics of large construction capacity,strong regulation capacity,long service life and stable operation.As the most commonly used peak shaving power supply at present,pumped storage unit can flexibly adjust the water volume of reservoir according to the load change,realize the function of peak shaving and valley filling,and provide corresponding flexible resources according to its operation state at each time scale in a day.The research and application of electrochemical energy storage are also relatively mature.Its full power response is in the millisecond level,with fast frequency modulation speed,high regulation accuracy,high charge and discharge efficiency,short construction cycle and flexible configuration.It can stabilize load fluctuation and improve power quality,and can be used for flexible regulation in the medium frequency band.

    In the operation of power system,the application scenarios of 15-minute scale and hourly scale are more practical.Therefore,this paper analyzes the optimal deployment of energy storage based on the intermediate frequency and low frequency components of net load curve.Combined with the unit scale,unit output characteristics,energy storage technical characteristics and whether the practical application is mature,the coordination of three flexible resources of thermal power unit,pumped storage unit and battery energy storage are mainly considered in case study.The power supply structure of conventional thermal power,wind power and photovoltaic in the case study is shown in Table 2.Related parameters in case study is shown in Table 3.

    Table 2: Structure of power unit

    Table 3: Example related parameters

    Table 3 (continued)

    Table 4 shows the optimal deployment results of hybrid energy storage.

    Table 4: Optimal deployment results of hybrid energy storage

    The daily average output fluctuation of wind power has the characteristic of inverse peak regulation,and the downhill section of photovoltaic noon peak has the characteristic of inverse peak regulation.Due to the joint access of wind power and photovoltaic,the two characteristics are superimposed,which greatly increases the peak regulation demand of the system.Therefore,the deployed capacity of pumped storage units on the hourly scale is large.For fluctuations on the 15-minute scale,flexible adjustment is accomplished by configuring all vanadium redox flow battery.In terms of investment cost,large scale of pumped storage units leads to high investment cost.In order to reduce the allocation demand of pumped storage units,the impact of wind and photovoltaic power on system peak shaving can be improved by optimizing the access ratio of wind and photovoltaic power,and other flexible resources such as controllable load can also be introduced for collaborative regulation.The optimal deployment results of single energy storage are shown in Table 5.

    Table 5: Optimal deployment results of single energy storage

    In the single storage scheme of pumped storage unit,the configured power in each scale is close to that in the hybrid energy storage scheme,which can meet the flexibility requirements of the system in each scale and reduce the occurrence of wind,photovoltaics and load rejection.However,the system investment cost increases compared with the hybrid energy storage scheme.Due to frequent regulation of pumped storage unit on a short timescale,it is operation life is significantly reduced,which increases the average daily investment cost of the pumped storage unit,so that the total investment cost is greater than the result of hybrid energy storage.In the single storage scheme of vanadium flow battery,due to the large requirements for energy storage capacity and the high cost of unit capacity and power,the investment of energy storage unit is significantly increased.Moreover,in fact the flow battery unit cannot meet flexibility requirements with weak regulation ability on the hourly scale.Additionally,it’s service life is short.Therefore,the comparison of examples reflects that according to the output characteristics of energy storage,the use of hybrid energy storage can not only ensure the flexibility of the system,but also reduce the energy storage deployment cost and improve the economy of the system.

    Combined with the peak and valley distribution of typical load curve,the controllable load is adjusted in different periods.During 1:00–9:00,the unit adjustment cost is 300 yuan/MWh; During the peak hours of 11:00–13:00 and 17:00–22:00,the unit adjustment cost is 300 yuan/MWh.Because different controllable loads have different response times and durations,the response loads on hourly and 15-minute time scales account for 6%and 3%of the maximum total load respectively.The energy storage allocation results with the participation of controllable load are shown in Table 6.

    Table 6: Energy storage deployment results with controllable load participation

    Added the controllable load,energy storage capacity of each scale is reduced,indicating that controllable load participates in the peak shaving task and reduces the demand for energy storage deployment.During the downward climbing period of load in the early morning,the output of thermal unit is limited by its climbing constraints,and part of the downward climbing load is borne by the controllable load.In the evening peak period,the system alleviates the upward climbing pressure of thermal unit and pumped storage unit by reducing part of the controllable load.Therefore,the flexible regulation of the controllable load in multiple periods reduces the peak shaving pressure of the system.The addition of controllable load provides a relatively economic choice for system peak shaving.

    6 Conclusion

    Based on the output characteristics and complementary characteristics of wind power and photovoltaic,this article forecasts net load through a Bayesian deep neural network,and proposes a power system flexibility evaluation method considering multi-scale decomposition of flexibility demand and a bi-level energy storage deployment optimization model.The multi-scale output characteristics,flexibility evaluation and energy storage deployment of new energy are considered,and the following valuable conclusions are drawn.

    Due to the different fluctuation characteristics of wind and photovoltaics output on each time scale,high-penetration renewable energy has multi-timescale effects on the flexibility of the power system.The multi-scale morphological algorithm can accurately obtain separated multiscale components in flexibility scale division,which is conducive to the determination of flexibility requirement in a specific fluctuation period.Considering the impact of multi-timescale flexibility demand caused by multi-timescale output of wind and photovoltaic power,a power system flexibility demand evaluation method based on multi-scale morphology is proposed to effectively distinguish multi-timescale flexibility requirement.

    According to the multi-scale characteristics of flexibility demand in power system with renewable energy,considering the economy of energy storage deployment and the flexibility of system operation,a multi-scale energy storage deployment model based on bi-level decision-making is established.This model can balance flexibility requirements and flexibility resources on each time scale according to the situation of wind and photovoltaic power output in the system,and improve the flexibility of power supply.

    Funding Statement:This research was supported by the National Natural Science Foundation of China(Grant No.52107129).

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    久久精品亚洲av国产电影网| 国产极品粉嫩免费观看在线| 天天操日日干夜夜撸| 最新的欧美精品一区二区| 亚洲少妇的诱惑av| 中文字幕精品免费在线观看视频| 纵有疾风起免费观看全集完整版| 亚洲男人天堂网一区| 乱人伦中国视频| 午夜日韩欧美国产| 日韩大片免费观看网站| 大话2 男鬼变身卡| 熟女av电影| 国产在线一区二区三区精| 丝袜脚勾引网站| 久久99热这里只频精品6学生| 亚洲国产av新网站| 欧美日韩一区二区视频在线观看视频在线| 欧美在线黄色| 国产成人欧美| 亚洲综合精品二区| 成人18禁高潮啪啪吃奶动态图| 国产成人91sexporn| 视频在线观看一区二区三区| 国产一区二区激情短视频 | 亚洲一区二区三区欧美精品| 国产精品一区二区在线观看99| 一区福利在线观看| 女的被弄到高潮叫床怎么办| www日本在线高清视频| 午夜免费鲁丝| 可以免费在线观看a视频的电影网站 | 久久99精品国语久久久| 国产av精品麻豆| 免费在线观看完整版高清| 国产精品麻豆人妻色哟哟久久| 国产无遮挡羞羞视频在线观看| 亚洲自偷自拍图片 自拍| 久久久久久久久免费视频了| 18禁裸乳无遮挡动漫免费视频| 制服诱惑二区| 亚洲精品aⅴ在线观看| 老司机影院毛片| 国产不卡av网站在线观看| 久久精品国产a三级三级三级| 在线免费观看不下载黄p国产| 亚洲国产最新在线播放| 亚洲国产精品一区三区| 一本色道久久久久久精品综合| 午夜av观看不卡| 国产精品久久久久久人妻精品电影 | 在线观看一区二区三区激情| 飞空精品影院首页| 狂野欧美激情性bbbbbb| 大片免费播放器 马上看| 黑丝袜美女国产一区| 在线精品无人区一区二区三| 久久狼人影院| 亚洲专区中文字幕在线 | 亚洲第一区二区三区不卡| 十八禁网站网址无遮挡| 在线观看免费视频网站a站| 男女午夜视频在线观看| 美女视频免费永久观看网站| 一二三四在线观看免费中文在| 99九九在线精品视频| xxxhd国产人妻xxx| 黑人猛操日本美女一级片| 国产97色在线日韩免费| 欧美日韩一级在线毛片| 亚洲一区二区三区欧美精品| 久久精品国产亚洲av高清一级| 国产在线视频一区二区| 久久毛片免费看一区二区三区| 搡老岳熟女国产| 自线自在国产av| 亚洲av综合色区一区| 日韩一区二区三区影片| 国产精品久久久久成人av| 毛片一级片免费看久久久久| 在线观看国产h片| 国产亚洲午夜精品一区二区久久| 精品一区二区免费观看| 免费观看a级毛片全部| 老鸭窝网址在线观看| 天堂中文最新版在线下载| 国产精品一区二区在线观看99| 久久精品国产亚洲av涩爱| 两个人免费观看高清视频| 国产日韩欧美在线精品| 午夜福利视频在线观看免费| 熟妇人妻不卡中文字幕| 飞空精品影院首页| 美女视频免费永久观看网站| 久久精品久久久久久久性| 国产精品久久久久久精品电影小说| 亚洲自偷自拍图片 自拍| 女人爽到高潮嗷嗷叫在线视频| 午夜日韩欧美国产| 91精品伊人久久大香线蕉| 久久99热这里只频精品6学生| 国产日韩一区二区三区精品不卡| www.熟女人妻精品国产| 美女扒开内裤让男人捅视频| 中国三级夫妇交换| 菩萨蛮人人尽说江南好唐韦庄| 日本午夜av视频| 亚洲美女黄色视频免费看| 18在线观看网站| 最新的欧美精品一区二区| 久久天躁狠狠躁夜夜2o2o | 母亲3免费完整高清在线观看| 永久免费av网站大全| 国产免费福利视频在线观看| 男女之事视频高清在线观看 | 国产爽快片一区二区三区| 丝袜在线中文字幕| 国产成人欧美| 欧美成人精品欧美一级黄| 亚洲精品成人av观看孕妇| 亚洲精品第二区| 免费黄色在线免费观看| 秋霞在线观看毛片| 少妇人妻精品综合一区二区| 97精品久久久久久久久久精品| 欧美日韩精品网址| 美女国产高潮福利片在线看| 女性被躁到高潮视频| 不卡av一区二区三区| 日韩一区二区三区影片| 成人午夜精彩视频在线观看| 国产在线视频一区二区| 久久久久久久国产电影| 日韩一区二区三区影片| 一区二区三区激情视频| 一区二区三区乱码不卡18| 少妇猛男粗大的猛烈进出视频| 中文字幕av电影在线播放| 国产精品一区二区在线不卡| 一区二区av电影网| 欧美日韩一级在线毛片| 香蕉国产在线看| 亚洲色图 男人天堂 中文字幕| 成人免费观看视频高清| 两个人免费观看高清视频| 亚洲视频免费观看视频| 两个人免费观看高清视频| 国产精品国产av在线观看| 男女午夜视频在线观看| 国产一区有黄有色的免费视频| 亚洲第一av免费看| 亚洲少妇的诱惑av| 男女午夜视频在线观看| 久久久国产欧美日韩av| 在线观看一区二区三区激情| 日韩视频在线欧美| 欧美97在线视频| 涩涩av久久男人的天堂| 啦啦啦啦在线视频资源| 国产又爽黄色视频| 欧美 日韩 精品 国产| 久久久久久久大尺度免费视频| 晚上一个人看的免费电影| 亚洲精品国产区一区二| 悠悠久久av| 国产成人欧美| 国产欧美日韩一区二区三区在线| 午夜精品国产一区二区电影| 18禁观看日本| 国产野战对白在线观看| 久热这里只有精品99| 国产又色又爽无遮挡免| 午夜福利,免费看| 美女高潮到喷水免费观看| 性高湖久久久久久久久免费观看| 黄频高清免费视频| 国产精品 国内视频| 制服诱惑二区| 免费不卡黄色视频| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡人人看| 操美女的视频在线观看| 日韩中文字幕视频在线看片| 亚洲国产中文字幕在线视频| 欧美精品av麻豆av| 啦啦啦在线免费观看视频4| 一级毛片我不卡| 久久久久网色| 国产精品一二三区在线看| 18禁观看日本| 午夜影院在线不卡| av视频免费观看在线观看| 欧美变态另类bdsm刘玥| 99国产精品免费福利视频| 国产高清不卡午夜福利| 纵有疾风起免费观看全集完整版| 日韩大码丰满熟妇| 免费观看a级毛片全部| 一区二区三区精品91| 永久免费av网站大全| 欧美日韩亚洲综合一区二区三区_| 国产在线视频一区二区| 最近中文字幕2019免费版| 亚洲第一区二区三区不卡| 又大又爽又粗| 国产极品粉嫩免费观看在线| 国产精品人妻久久久影院| 大片电影免费在线观看免费| 欧美久久黑人一区二区| 热re99久久国产66热| 最近的中文字幕免费完整| 悠悠久久av| 国产精品久久久av美女十八| 80岁老熟妇乱子伦牲交| 精品酒店卫生间| 天堂中文最新版在线下载| 欧美日韩成人在线一区二区| 国产成人免费无遮挡视频| 一边摸一边做爽爽视频免费| 久久久久久久大尺度免费视频| 免费在线观看完整版高清| 亚洲欧美精品自产自拍| 免费高清在线观看视频在线观看| 国产一区二区三区av在线| 人人妻人人澡人人看| 亚洲精品乱久久久久久| 国产成人啪精品午夜网站| av片东京热男人的天堂| 欧美另类一区| av不卡在线播放| 免费黄网站久久成人精品| 国产一卡二卡三卡精品 | 国产精品成人在线| 久久久国产一区二区| 精品少妇内射三级| 国产亚洲一区二区精品| 国产成人系列免费观看| 十八禁网站网址无遮挡| 日韩人妻精品一区2区三区| av片东京热男人的天堂| 久久久久久久久免费视频了| 精品少妇久久久久久888优播| 亚洲国产最新在线播放| 99久久99久久久精品蜜桃| 99九九在线精品视频| 侵犯人妻中文字幕一二三四区| 韩国精品一区二区三区| 国产淫语在线视频| 十八禁网站网址无遮挡| 在线看a的网站| 美女主播在线视频| 国产成人91sexporn| 亚洲av电影在线观看一区二区三区| 国产成人免费无遮挡视频| 日本wwww免费看| 成人国产麻豆网| 国产亚洲午夜精品一区二区久久| 美女福利国产在线| 欧美成人午夜精品| 波多野结衣av一区二区av| 国产精品二区激情视频| 精品国产一区二区三区久久久樱花| 亚洲第一区二区三区不卡| 成人亚洲欧美一区二区av| 亚洲成人av在线免费| 制服丝袜香蕉在线| 亚洲精品,欧美精品| 久久久国产一区二区| 精品少妇内射三级| 日韩中文字幕欧美一区二区 | 久久天堂一区二区三区四区| 国产精品国产三级专区第一集| 日韩一卡2卡3卡4卡2021年| 激情五月婷婷亚洲| 99久久99久久久精品蜜桃| 免费看av在线观看网站| 少妇人妻精品综合一区二区| 精品少妇内射三级| 婷婷色av中文字幕| 美女中出高潮动态图| 成人亚洲欧美一区二区av| 一区在线观看完整版| 精品一区在线观看国产| 人妻人人澡人人爽人人| 亚洲成人一二三区av| 午夜福利一区二区在线看| 国产av精品麻豆| 日韩伦理黄色片| 91老司机精品| 高清欧美精品videossex| 亚洲专区中文字幕在线 | 日韩中文字幕视频在线看片| 卡戴珊不雅视频在线播放| 久久久久久免费高清国产稀缺| 国产av一区二区精品久久| 欧美在线一区亚洲| 久久这里只有精品19| 日本爱情动作片www.在线观看| a 毛片基地| 丝袜喷水一区| 国产免费现黄频在线看| 亚洲精品日韩在线中文字幕| 久久久久国产一级毛片高清牌| 99久久综合免费| 女人被躁到高潮嗷嗷叫费观| 人妻一区二区av| 电影成人av| 日日撸夜夜添| 99久久综合免费| 欧美日韩亚洲高清精品| 久久天堂一区二区三区四区| 交换朋友夫妻互换小说| 精品久久蜜臀av无| 国产免费一区二区三区四区乱码| 久久久精品区二区三区| 久久久精品免费免费高清| 免费观看av网站的网址| 日本wwww免费看| 美女福利国产在线| 女人被躁到高潮嗷嗷叫费观| 亚洲,欧美,日韩| 日韩人妻精品一区2区三区| 一级片'在线观看视频| 人妻人人澡人人爽人人| 亚洲av在线观看美女高潮| 亚洲中文av在线| 亚洲美女黄色视频免费看| 成人18禁高潮啪啪吃奶动态图| 免费黄色在线免费观看| 精品人妻在线不人妻| 亚洲av欧美aⅴ国产| 啦啦啦在线观看免费高清www| 2021少妇久久久久久久久久久| 丁香六月天网| 人体艺术视频欧美日本| 99久久精品国产亚洲精品| 高清不卡的av网站| 超碰97精品在线观看| 亚洲欧美成人精品一区二区| 麻豆av在线久日| av天堂久久9| 麻豆乱淫一区二区| 国产精品国产av在线观看| 视频在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久 | xxxhd国产人妻xxx| 婷婷成人精品国产| 不卡视频在线观看欧美| 一区在线观看完整版| 久久精品久久精品一区二区三区| 一区二区三区激情视频| 久久毛片免费看一区二区三区| 亚洲精品自拍成人| 夫妻午夜视频| 亚洲精品第二区| 男女无遮挡免费网站观看| 久久精品人人爽人人爽视色| 精品一区二区免费观看| 精品少妇久久久久久888优播| 久久久久久免费高清国产稀缺| 80岁老熟妇乱子伦牲交| 欧美精品一区二区大全| 免费av中文字幕在线| 亚洲欧美精品综合一区二区三区| 久久久精品免费免费高清| 亚洲av国产av综合av卡| av天堂久久9| 超碰97精品在线观看| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 少妇猛男粗大的猛烈进出视频| 久久精品国产a三级三级三级| 欧美日韩综合久久久久久| 色吧在线观看| 深夜精品福利| 欧美精品一区二区免费开放| 高清不卡的av网站| 人人妻,人人澡人人爽秒播 | 国产精品久久久av美女十八| 最近最新中文字幕大全免费视频 | 一边亲一边摸免费视频| 老司机深夜福利视频在线观看 | 国产亚洲最大av| 免费看av在线观看网站| www.av在线官网国产| 丰满乱子伦码专区| 久久人人爽人人片av| 免费在线观看完整版高清| 99精品久久久久人妻精品| 色播在线永久视频| 国产精品免费视频内射| 欧美黄色片欧美黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| 高清不卡的av网站| 极品少妇高潮喷水抽搐| 国产毛片在线视频| 制服丝袜香蕉在线| 国产高清国产精品国产三级| 国产乱来视频区| 黄色视频不卡| 免费人妻精品一区二区三区视频| 久久青草综合色| 十八禁高潮呻吟视频| 美女脱内裤让男人舔精品视频| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 久久天堂一区二区三区四区| 亚洲av欧美aⅴ国产| 欧美久久黑人一区二区| 美女国产高潮福利片在线看| 宅男免费午夜| 99热国产这里只有精品6| 久久精品久久精品一区二区三区| 亚洲伊人色综图| 纵有疾风起免费观看全集完整版| 搡老乐熟女国产| 国产日韩欧美视频二区| 另类亚洲欧美激情| 99香蕉大伊视频| 咕卡用的链子| 成人午夜精彩视频在线观看| 哪个播放器可以免费观看大片| 国产成人a∨麻豆精品| 搡老岳熟女国产| 久久天躁狠狠躁夜夜2o2o | 久久久久精品国产欧美久久久 | 色网站视频免费| 午夜激情久久久久久久| 久久久久久免费高清国产稀缺| www.av在线官网国产| 大码成人一级视频| 侵犯人妻中文字幕一二三四区| 欧美激情极品国产一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久精品电影小说| 黄片小视频在线播放| 极品人妻少妇av视频| 在线 av 中文字幕| 欧美人与善性xxx| 极品人妻少妇av视频| 最近的中文字幕免费完整| 99国产精品免费福利视频| e午夜精品久久久久久久| 啦啦啦 在线观看视频| bbb黄色大片| 欧美日韩av久久| 久久毛片免费看一区二区三区| 免费在线观看黄色视频的| 老司机亚洲免费影院| 欧美黄色片欧美黄色片| av片东京热男人的天堂| 国产日韩欧美视频二区| 色网站视频免费| 久久午夜综合久久蜜桃| 伊人久久国产一区二区| 日韩av在线免费看完整版不卡| 久久 成人 亚洲| 久久精品aⅴ一区二区三区四区| 欧美日韩视频精品一区| 一本久久精品| 亚洲成av片中文字幕在线观看| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 老司机影院毛片| 欧美另类一区| 丰满饥渴人妻一区二区三| 丝袜美腿诱惑在线| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 免费人妻精品一区二区三区视频| 又黄又粗又硬又大视频| 欧美黄色片欧美黄色片| 国产男人的电影天堂91| 亚洲精品国产av蜜桃| 欧美xxⅹ黑人| 国产成人免费无遮挡视频| 亚洲人成网站在线观看播放| 成人国产麻豆网| 国产亚洲av高清不卡| 天天添夜夜摸| 观看美女的网站| 久热这里只有精品99| 国产成人免费观看mmmm| 国产一级毛片在线| 男女之事视频高清在线观看 | 国产极品粉嫩免费观看在线| 一二三四在线观看免费中文在| 国产熟女欧美一区二区| 欧美日韩视频高清一区二区三区二| 男人添女人高潮全过程视频| 男女床上黄色一级片免费看| 日韩一卡2卡3卡4卡2021年| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 最近中文字幕高清免费大全6| 国产淫语在线视频| 少妇人妻久久综合中文| 婷婷成人精品国产| 日韩人妻精品一区2区三区| 国语对白做爰xxxⅹ性视频网站| 精品人妻熟女毛片av久久网站| 午夜久久久在线观看| 国产在线一区二区三区精| 亚洲精品久久久久久婷婷小说| av不卡在线播放| 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 80岁老熟妇乱子伦牲交| 成人国语在线视频| 韩国精品一区二区三区| 如日韩欧美国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 中文字幕人妻丝袜制服| 毛片一级片免费看久久久久| 欧美人与善性xxx| 男女边吃奶边做爰视频| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区久久| 老汉色∧v一级毛片| 午夜日本视频在线| 亚洲四区av| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| 亚洲欧美成人精品一区二区| 国产激情久久老熟女| 人人妻,人人澡人人爽秒播 | 亚洲男人天堂网一区| 男女国产视频网站| 成年动漫av网址| 日韩制服骚丝袜av| 亚洲成人一二三区av| 天天躁夜夜躁狠狠久久av| 热99国产精品久久久久久7| 美女脱内裤让男人舔精品视频| 国产免费视频播放在线视频| 久久精品国产亚洲av涩爱| 精品午夜福利在线看| 在线观看人妻少妇| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| 中文字幕人妻熟女乱码| 亚洲三区欧美一区| 汤姆久久久久久久影院中文字幕| 97人妻天天添夜夜摸| 丰满饥渴人妻一区二区三| av免费观看日本| 高清欧美精品videossex| 国产欧美亚洲国产| 国产亚洲最大av| 久久久精品区二区三区| 亚洲一区二区三区欧美精品| 最近中文字幕高清免费大全6| 亚洲国产精品999| 51午夜福利影视在线观看| 午夜免费鲁丝| 乱人伦中国视频| av.在线天堂| 精品卡一卡二卡四卡免费| avwww免费| 一区二区av电影网| 热re99久久国产66热| 免费黄网站久久成人精品| 叶爱在线成人免费视频播放| 色综合欧美亚洲国产小说| 欧美精品高潮呻吟av久久| 国产精品偷伦视频观看了| 国产精品国产三级国产专区5o| 久久精品人人爽人人爽视色| 毛片一级片免费看久久久久| av有码第一页| 精品一区在线观看国产| 国产精品久久久人人做人人爽| 国产成人91sexporn| 九九爱精品视频在线观看| 久久天躁狠狠躁夜夜2o2o | 色婷婷av一区二区三区视频| 99热全是精品| 亚洲精品视频女| 美女视频免费永久观看网站| 男的添女的下面高潮视频| 欧美精品亚洲一区二区| av有码第一页| 国产精品久久久久成人av| 国产免费福利视频在线观看| 建设人人有责人人尽责人人享有的| bbb黄色大片| 国产深夜福利视频在线观看| av网站在线播放免费| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 1024香蕉在线观看| 亚洲综合色网址| 精品国产超薄肉色丝袜足j| 欧美xxⅹ黑人| 久久久久网色| 色网站视频免费| 久久久久网色| 人成视频在线观看免费观看| 日韩欧美精品免费久久| 亚洲欧美精品自产自拍| 亚洲精品av麻豆狂野| 国产熟女欧美一区二区| 中文乱码字字幕精品一区二区三区| 你懂的网址亚洲精品在线观看| 99久国产av精品国产电影| 国产免费现黄频在线看| 久久精品国产综合久久久| 欧美日本中文国产一区发布| av在线播放精品| 另类亚洲欧美激情| 9色porny在线观看| 国产老妇伦熟女老妇高清| 在线天堂最新版资源| 男人爽女人下面视频在线观看| 欧美激情高清一区二区三区 |