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

    A Storage and Transmission Joint Planning Method for Centralized Wind Power Transmission

    2021-12-14 10:29:40XiuyuYangQiGuoJianzhongGuiRenyongChaiandXueyuanLiu
    Computers Materials&Continua 2021年7期

    Xiuyu Yang,Qi Guo,Jianzhong Gui,Renyong Chai and Xueyuan Liu

    1Key Laboratory of Modern Power System Simulation and Control&Renewable Energy Technology(Northeast Electric Power University),Jilin,132012,China

    2Electrical Engineering-PHD program at the University of Idaho,Idaho,83843,America

    3State Grid Jiangxi Power Supply Company,Xinyu,338025,China

    Abstract: Centralized delivery has become the main operation mode under the scaled development of wind power.Transmission channels are usually the guarantee of out-delivered wind power for large-scale wind base.The configuration of transmission capacity,which has the features of low utilization and poor economy, is hardly matching correctly due to the volatility and low energy density of wind.The usage of energy storage can mitigate wind power fluctuations and reduce the requirement of out-delivery transmission capacity,but facing the issue of energy storage cost recovery.Therefore, it is necessary to optimize the allocation of energy storage while considering the problem of wind power transmission.This paper studies the joint optimization of large-scale wind power transmission capacity and energy storage, reveals the mechanism of energy storage in order to reduce the power fluctuation of wind power base and slow down the demand of transmission.Then,analyze the multi-functional cost-sharing mode of energy storage, improve the efficiency of energy storage cost recovery.Constructs the coordination optimization configuration model to deal with the problem of large-scale wind power transmission capacity and energy storage,and realizes the transmission capacity optimization coordination and optimization with energy storage.The proposed method is verified by a wind base located in Northeast China.

    Keywords: Wind farm group; energy storage configuration; transmission capacity; coordination and optimization

    1 Introduction

    Wind power is developed vigorously in the last 10 years, which is a major measure of many countries to promote the transition to low-carbon energy and respond to climate change.The overall capacity of all wind turbines installed worldwide by the end of 2019 reached 65,100 MW,which is 5 times more than the capacity in 2009 according to statistics showed in [1].Taking China as an example, 10 tens of millions of kilowatt wind power bases are successively built in many provinces, such as Gansu, Xinjiang, Hebei, Jilin, Jiangsu, and Shandong.The installed wind power capacity is predicted to exceed 1.1 billion kW in 2030 and reach 2.4 billion kW in 2050 [2].Wind power is developed on large scale based on wind energy distribution in China,and the centralized delivery of wind power is the main operating mode to transmit power to the substations outside the wind base.

    The wind power base usually is far from the load center with small local power demand in the weak grid.Supporting power transmission projects are constructed to ensure the effective transmission of wind power [2,3].Compared to the conventional power supply, wind power includes the obvious features of volatility, intermittent, and low energy density.Two flexible options have been used in the power grid to configure the transmission capacity.One solution is that transmission capacity is allocated according to the installed capacity of the wind power base,which causes over-allocation of transmission capacity, resulting in increased costs of transmission assets.Another one is to build limited transmission capacity, which introduces the waste risk of wind power due to transmission congestion, even with related low transmission investment costs.The remaining issues mentioned above are going to be solved by the Energy Storage System(Energy Storage System, ESS).The reason is that ESS has the ability to transfer power and stored energy in time.The working mode of charging at high wind power generation times and discharging during low wind power generation time has the advantages of reducing peak wind power, slowing down the demand for transmission resources, increasing the utilization rate of transmission channels, and improving the economics of transmission projects.Therefore, the collaborative planning of transmission channel capacity and energy storage is of great significance for large-scale wind power under the background of the rapid development of renewable energy and energy storage technology.

    Nowadays, several research results have been achieved in the planning of power out-delivery for large-scale wind power.The authors in [4] proposed a coordinated two-level planning method for a large-scale wind farm integrated system and related regional transmission network.The contradiction between the expansion speed of transmission capacity and the speed of developing installed wind capacity is solved under the background of large-scale centralized development of intermittent wind energy.Due to a lack of ability to handle uncertainty in the traditional deterministic planning model, the authors in [5] consider the variability of wind speed and the uncertainty of demand and propose an algorithm to solve the transmission expansion planning problem in the generation period of large-scale wind power.Reference [6] modeled many uncertain factors and proposed a new expected value model of transmission network expansion planning based on DC power flow.

    The optimized configuration of energy storage is an effective way to deal with the fluctuation of renewable energy output and insufficient system flexibility [7], which has been a hot topic for research.Energy storage plays a critical role in the power system, such as wind power fluctuation suppression [8], frequency response [9,10], spinning reserve [11], peak shaving [12,13] as well as mitigating congestion [14].The energy storage system has changed the energy distribution of the original power system.Reference [15–17] provides the method of energy management and distribution.In addition, extensive and fruitful research has been conducted on energy storage and transmission grid planning for renewable energy integration.An approximate optimization method is used in [18] to achieve an optimal configuration of distributed energy storage site selection and constant capacity in the transmission grid.The long-term and short-term uncertainties of high-permeability renewable energy are solved by a joint planning method proposed in [19]for energy storage and transmission grids.The authors in [20] build a new robust minimum–max–minimum cost optimization model with the consideration of topology, aiming for joint transmission and energy storage expansion planning where wind farms are connected to the power system considering topology optimization.Reference [21–26] provides algorithmic ideas for model optimization.Considering many costs, such as the cost of transmission grid expansion planning, energy storage construction and operation, unit combination, and wind abandonment,authors achieve the goal of comprehensive planning for wind farm installation and location,energy storage, and transmission grid expansion in [27].However, the interaction between planning variables is still unclear from the above references.From the perspective of wind curtailment,the various functions of energy storage and transmission grids are verified in [28,29], and it is shown that energy storage is both a supplement to and a substitute for transmission lines,and the combined effect of the two is better.In the electricity market, the economy of largecapacity energy storage technology is evaluated in [30].The fact that energy storage can mitigate transmission congestion and bring benefits, but the benefits are reduced or eliminated after the expansion of the transmission grid, is explained.Considering the energy storage life cycle benefits and return on investment helps to establish a framework for joint planning of energy storage site selection and capacity determination and transmission grid in [31].Moreover, a dynamic joint planning method is proposed for energy storage and transmission grid under consideration of the life and performance of energy storage in [32].

    All the above mentioned, there is multi-angle consideration of the problems between energy storage and transmission grids, but rarely including transmission capacity matching of a wind power base.Therefore, this paper studies the joint optimization problem of large-scale wind power transmission capacity and energy storage.The main contributions are as follows:

    (1) Coordinated planning of energy storage and capacity of outgoing transmission channels of wind power bases can better deal with wind power with high amplitude and small probability, thus improving the economy of transmission engineering;

    (2) The 365-day scenario is divided into high wind, stroke, and low wind conditions.The peakvalley arbitrage of energy storage in low wind conditions can reduce the unit cost of the energy storage system and overcome the problem of the high cost to a certain extent.

    The rest of the paper is organized as follows.The coordination and optimization principle of transmission capacity and energy storage is introduced in Section 2.The model of coordination optimization configuration of transmission capacity and energy storage is proposed in Section 3.Section 4 presents the coordination and optimization simulations with compared results using the proposed method and other methods.Conclusions and a discussion of future work are provided in Section 5.

    2 Coordination and Optimization Principle of Transmission Capacity and Energy Storage

    Large-scale wind power and conventional power sources both rely on transmission channel networking, but their output characteristics are completely different from conventional generation units.There is a contradiction between the wind power transmission rate and the economics of the transmission channel.

    2.1 Economic Analysis of Transmission Capacity of Wind Power Base

    The output power with volatility is less than the installed capacity in the wind power base,as shown in Fig.1 (left).Due to the temporal and spatial differences between the wind farms in the wind power base, the output power volatility will become weaker as the scale of the wind farm group increasing, and the maximum per-unit value decreases with the increasing capacity of the wind power base.The phenomenon is called the “convergence effect” [33].The long-term fluctuation characteristics of large-scale wind power are characterized by using the power curve of wind power, as shown in Fig.1.When planning the external transmission channel of the wind farm group, if the transmission capacity Plineequals to the installed capacity of the wind farm group PN, the transmission capacity between PNand the maximum wind power Pmaxwill not be used, which is called invalid capacity, resulting in the waste of investment, and is deteriorating as the scale of the wind farm group increasing.If Plineequals to Pmax, the wind power in the neighborhood will be below Pmaxwith large magnitude and short duration (only a few hours to tens of hours throughout the year).This transmission capacity with lower utilization is called“inefficient capacity.” If Plineis lowered manually, it is easy to cause transmission congestion and a large amount of wind abandonment.

    Figure 1:Economics analysis of wind power transmission capacity

    2.2 Economic Analysis of Energy Storage Configuration in Outward Transmission Project

    The energy storage system at the exit of the wind farm group is shown in Fig.2.Using energy storage system to capture the “l(fā)arge value and small probability” wind power between Plineand Pmax, it can not only reduce the transmission capacity, save the investment in transmission projects,but also alleviate the blockage and abandon wind, and better solve the dilemma of insufficient transmission capacity allocation of large-scale wind power transmission channels.

    Figure 2:The sketch diagram of centralized grid-connected wind power

    However, the reason why energy storage is not applied on large scale is that it has a high cost and single energy storage function for planning with difficulty in cost recovering.The economics of energy storage systems is analyzed as following when applied in large-scale centralized delivery.

    Taking the wind farm group shown in Fig.1 as an example, the maximum output value is 0.7 p.u., and the corresponding time series of output features are shown in Fig.3.As the transmission capacity Plinedecreasing, the decreasing trend of the blocked abandonment power value ΔP (is Pmax?Pline), the abandonment duration Tdur, and the abandonment power ΔE can be noticed directly.Three types of typical operating conditions of wind power can be divided based on the relationship between wind power and transmission capacity from the output feature curve.(a) Strong wind conditions:the minimum value Pminof wind power for the whole day is bigger than Pline.There is no discharge space for the energy storage system, and wind power exceeding Plinewill be abandoned.(b) Moderate wind conditions:the maximum value of wind power Pmaxis bigger than Pline, and the minimum value of wind power Pminis smaller than Pline.The energy storage system captures the abandoned wind.(c) Weak wind conditions:the maximum value of wind power Pmaxis smaller than Plinethroughout one day, and there is no blocking and abandoning wind, as shown in Fig.3.The occurrence times of the three-type working conditions are T1, T2, T3.

    Figure 3:Three typical working conditions of wind power

    When the transmission capacity Plinedecreases from 0.65 to 0.40 by the step of 0.05, ΔP, ΔE,Tdur, T1, T2, and T3are shown in Tab.1.

    From Tab.1, when the transmission capacity is 0.4, the duration time of wind abandonment is 485h, accounting for only 5.54% of the annual time, and the wind abandonment conditions(strong wind conditions, moderate conditions) is lasting a total of 89 days, accounting for 24.38%of the whole year.If the energy storage system only captures blocking and abandoning wind in the power transmission project, its working time is only T2.The utilization rate of energy storage assets is extremely low, and the cost is difficult to recover.When coordinating the planning with the transmission capacity, the high cost of energy storage per unit time results in a small amount of energy storage configuration, or even no energy storage.

    Table 1:ΔP, ΔE, Tdur, T1, T2, and T3 under different Pline

    2.3 Energy Storage Cost Allocation Model

    To improve the economic feasibility of the energy storage system, reduce energy storage cost in unit time, a cost-sharing model based on the time ratio of the multi-function tasks of the energy storage system is proposed, in which the energy storage system can participate in the auxiliary peak shaving service (or peak-valley arbitrage) of the power system in the weak wind condition, except for the blocking and abandonment of the wind in the moderate condition,improving the recovery efficiency of energy storage cost.The formula of cost allocation model is as follows:

    where,Itot(Pess,Eess) is the total investment cost of the energy storage system, andPessis the charge/discharge power of the energy storage system;Eessis the capacity of the energy storage system, andKpis the purchase cost of the unit power converter;Keis the purchase cost per unit capacity of the energy storage system;λ1is the proportion of the total cost of the energy storage system under strong wind conditions;λ2is the proportion of the cost of the energy storage system to capture and abandon the wind under moderate conditions;λ3is the cost ratio of energy storage auxiliary peak shaving (peak-valley arbitrage) under weak wind conditions;RN?0(Pess,Eess) is the net income of the energy storage system under strong wind conditions.Because energy storage cannot be dispatched under strong wind conditions, the net income is negative and its cost;RN?L(Pess,Eess)is to capture the net income of obstructed abandonment,RL(Pess,Eess) is to capture the income of abandonment;RN?P(Pess,Eess) is the net income of peak and valley arbitrage,andRP(Pess,Eess) is the peak and valley arbitrage income.

    2.4 Transmission Capacity and Energy Storage Coordination Principle of Optimal Allocation

    The essence of coordination and optimization of the transmission capacity for wind power base and energy storage system is to use energy storage to capture “l(fā)arge value and small probability” peak wind power, reduce the dependence of wind power on transmission capacity, thereby delaying transmission line investment.Namely, in the case of maximizing the capture of abandoned wind, it’s a game between the net income of transmission and the net income of energy storage, whose goal is to maximize the comprehensive income.The variables of transmission capacity, energy storage charge/discharge power, and energy storage capacity should be optimized.

    The coordination and optimization procedure of energy storage and transmission capacity is shown in Fig.4.

    Figure 4:Coordinated optimal planning of outgoing

    Firstly, given the transmission capacity Pline, use the continuous output curve of the wind farm group to calculate the transmission project investment and transmission income, and then obtain the net transmission income;

    Secondly, to prevent the effect of the energy storage system from being affected when switching between adjacent operating conditions, the energy storage system meets the full charge and discharge constraints within the day, that is, the state of charge (SOC) should meet:SOC(0)=SOC(24).The energy storage is configured to maximize the blocking and abandonment of wind in moderate conditions.Then, according to the energy storage cost allocation model, calculate the energy storage to capture the net income of wind power.The configuration capacity Eessof the energy storage system is the maximum value of the curtailed wind power captured in all moderate conditions.Among them, the energy storage for a single moderate condition is the smaller one between the charging space and the discharge space during the day.And the energy storage system charge/discharge power Pessis the maximum value of charge/discharge power in all moderate conditions, The charge/discharge power of a single moderate condition is the difference between the wind power and the transmission capacity corresponding to the smaller one between the charging space and the discharge space.The equations are as follows:

    where, Pess?iand Eess?iare the energy storage power and energy storage capacity on the i-th day in the annual moderate condition (capture of abandoned wind power); Wc?i, Wdc?iare the charging space and discharge space on the i-th day in the moderate condition; Pmax?i, Pmin?iare the maximum and minimum wind power on thei-th day in the moderate condition throughout the year.

    Thirdly, use the above-mentioned energy storage configuration to carry out peak-to-valley arbitrage under weak wind conditions to increase energy storage cost recovery efficiency.Among them, the energy storage is charged during the low load period and discharged during the peak load period.There will be no restriction when charging, and the remaining transmission capacity of the transmission line during the peak load period during discharge.Therefore, the peak-to-valley arbitrage power is the smaller of the remaining transmission space and the energy storage capacity, that is:

    where Wp?iis the remaining transmission power space during the peak load period on thei-th day in weak wind conditions.

    Finally, the comprehensive income is calculated using variables such as power transmission income, income from capturing wind abandonment, peak-to-valley arbitrage income, and wind abandonment loss, iteratively modify the Plineand finally obtain a joint optimization plan for the transmission channel capacity and energy storage.

    3 Coordination and Optimal Allocation Model of Transmission Capacity and Energy Storage

    To facilitate the analysis of the relationship between transmission power of transmission line,energy storage power and wind power, the transmission power of transmission line, energy storage charge/discharge power, and energy storage capacity are set to continuous values, and the energy storage system capacity-to-power ratio is set as an integer.The actual project is rounded according to the optimal solution to obtain a feasible project plan.

    The factors such as transmission revenue, transmission investment, energy storage system income, energy storage investment and wind abandonment loss compensation are all considered,and a coordinated optimization model of energy storage and transmission capacity is built to maximize the comprehensive benefits of transmission projects and energy storage systems.

    where, R(Pline) is the power grid transmission income; I(Pline) is the investment in the transmission project; L(Pess, Eess) is the value of the abandonment loss that cannot be captured by energy storage, and the wind power enterprise is fully compensated according to the wind power grid price.

    Using wind power continuous power curve calculation:R(Pline), and using wind power time series power curve to calculate RL(Pess, Eess), Rp(Pess, Eess) and L(Pess, Eess) and other variables.These variables are calculated as the following:

    (1) R(Pline) is the income generated by the transmission of wind power through the transmission channel, and the calculation formula is

    where Kris the revenue of the unit of wind power delivered by the power transmission company;Tsis the payback period of the transmission investment; Gwindis the annual out-delivery of wind power.The continuous power curve of the used year is calculated as follows:

    where Tlineis the continuous output time when the output power of the wind power base is higher than Pline; Tendis the total continuous output time of the wind power base.

    (2) Transmission project investment I(Pline) can be calculated as:

    where Kcis the transmission project cost per unit capacity and unit length; l is the transmission distance.

    (3) Energy storage system captures the benefits of blocking wind curtailment RL(Pess, Eess).Wind power exceeding the transmission capacity will be discarded due to transmission congestion,and this part of the wind power will be filled by thermal power.The energy storage system captures the curtailment of wind while also reducing corresponding carbon emissions.Therefore,the energy storage system will produce the benefits of increasing wind power generation and reducing carbon emissions.We can find the captured benefit using the following equation.

    where Ksis the on-grid electricity price for wind power; Kcis the cost of carbon emissions per unit of electricity;Teis the static recovery period of energy storage investment; Gcapis the abandoned wind power captured by energy storage.

    where Eess?iis the captured and abandoned wind power on the i-th day in moderate condition.

    (4) Loss value of wind curtailment that energy storage cannot capture L(Pess, Eess):

    where Glostis the abandonment of wind power caused by transmission congestion, and the equation is as follows:

    where Pw·maxis the maximum output power of the wind power base.

    (5) The cost of the energy storage system is as Eq.(1).

    (6) Benefits of energy storage system peak-valley arbitrage.

    The energy storage system uses the peak-to-valley price difference to conduct peak and valley arbitrage in weak wind conditions, which can be calculated as follows:

    where Kpis the peak-to-valley price difference,Wiis the peak-to-valley charge and discharge capacity of the i-th day in weak wind conditions, calculated as Eq.(7).

    4 Case Study

    4.1 Calculation Conditions

    Taking wind power of the wind power base located in the northeast of China as the study case, The wind power base is composed of 15 wind farms, with a total installed capacity of 2013 MW.The continuous output curve of the wind power base is shown in Fig.1.The maximum wind power is 1403 MW, concentrated interconnection and energy storage configuration diagram,as shown in Fig.5.

    Figure 5:Wind farm cluster centralized grid

    The conditions for calculation are as following:

    (1) Charges for transmission of unit wind power by transmission utility.

    Kr=0.009 dollars/(kW·h),

    (2) 220 kV integrated transmission project unit cost

    Kc=1.513 million dollars/(100 MW·km), Transmission distancel=200 km,

    (3) Power grid transmission congestion caused by wind abandon the wind loss compensation unit price (according to the wind power generation unit price calculation Ks= 0.0908 dollars/(kW·h).Unit power carbon charge Kc=0.0348 dollars (kW·h),

    (4) Unit electricity peak valley price difference KP= 0.0953 dollars/(kW · h), local peak electricity price 0.1468 dollars, the valley of electricity price 0.0514 dollars,

    (5) Transmission investment recovery period TS=20a,

    (6) Power inverter unit cost 0.7567 million dollars/MW,

    (7) Lithium-ion batteries are divided into unit capacity costs 2.27 million dollars/(MW·h),Energy storage efficiency of 80%,

    (8) Investment recovery period of energy storage system Te=10a.

    4.2 Coordination Optimization Results

    In the case of wind power base size, the wind sources are approximately the same every year, and the relationship between wind power and transmission capacity is relatively stable.The proposed coordination optimization model can meet the requirements of planning for transmission capacity and energy storage.Transmission capacity and energy storage configuration optimization results are shown in Tab.2.

    Table 2:Results of coordinated allocation of transmission capacity and energy storage

    The changing trend of each variable is shown in Fig.6 during the optimization process.The comprehensive income increases with transmission capacity increasing, and it reaches the maximum value before transmission capacity equals 1115 MW, then it decreases with the increasing transmission capacity.Transmission income also increases with the increase in transmission capacity.When the capacity increases to the maximum wind power, the income remains constant.Transmission investment is increasing linearly.With the growth of transmission capacity,the demand for energy storage is reduced, and the working space is also compressed; when the maximum wind power of 1403 MW is reached, all wind power will be transmitted by the transmission line, without the need to configure energy storage.

    Figure 6:Sketch map of optimal transmission capacity for wind power output

    The optimal solution for transmission capacity and energy storage coordination is thatPlineequals 1115 MW, in which strong wind condition is 0 day, moderate wind condition is on 27th day, and weak wind condition is on the 338th day, as shown in Fig.7

    Figure 7:Wind power conditions at Pline=1115 MW

    With the development of energy storage technology, the cost drops further.The influence of energy storage continuous cost reduction on transmission capacity is analyzed as follows.Setting the cost of energy storage unit to three levels:(1) High-0.3027 million dollars/MW·h,(2) Medium-0.227 million dollars/MW · h, and (3) Low-0.1513 million dollars/MW · h.Three optimization results are calculated by using the proposed coordination planning model, shown in Tab.3 and Fig.8.From the comparison among the three schemes, energy storage power/capacity configuration increases as energy storage cost and transmission capacity gradually decrease.

    Table 3:Optimization results of different energy storage costs

    4.3 Methods of Comparison

    Under the same calculation conditions, the results calculated by the proposed coordination optimization planning model are compared with the two solutions calculated by the two methods in [33].From Tab.4, it shows that the full capacity planning method proposed in [33] is based on the installed capacity to plan the capacity of the transmission channel.Although no wind abandonment occurs, it has caused the over-allocation of transmission capacity and the worst economy.This method is optimized between the loss of wind curtailment and transmission investment, at the expense of a little wind curtailment in exchange for great savings in transmission investment.The method proposed in this paper adds energy storage to the planning of the power transmission project, which not only reduces the transmission capacityPline, but also alleviates the occurrence of wind curtailment.At the same time, it provides a certain auxiliary peak shaving capability for the system, and the economy is optimal.

    Figure 8:Optimal transmission capacity for different energy storage costs

    Table 4:Technical and economic comparison of different methods

    5 Conclusion

    The economic impact of energy storage is analyzed in the condition where energy storage participates in the optimization of transmission capacity for large-scale wind power base.The coordination configuration optimization model of large-scale wind power transmission capacity and energy storage is proposed, in which a multi-functional cost allocation mode of energy storage is applied based on the impact analysis.The proposed method is verified by the comparison of optimal solutions with other methods.

    The conclusions are drawn as follows:

    (1) A multi-functional cost-sharing model for energy storage is constructed depending on the different operating tasks of energy storage, which improves the efficiency of energy storage cost recovery and overcomes the problem of high energy storage cost and difficulty in cost recovery.

    (2) Compared with the method in literature [33], the coordinated optimization method for the capacity and energy storage of the outgoing transmission channel of large-scale wind power base proposed in this paper reduces the transmission capacity by 35 MW, reduces the wind abandonment loss by 575 million yuan, and increases the overall economic benefit by 11.229 billion yuan.The proposed method not only reduces the transmission capacity but also alleviates the occurrence of wind abandonment.At the same time, it provides a certain auxiliary peak regulation ability for the power system and obtains better transmission engineering benefits.

    (3) The coordination and optimization model of the capacity and energy storage of largescale wind power bases is proposed, which not only reduces the transmission capacity, but also alleviates the occurrence of wind curtailment.At the same time, it provides a certain amount of auxiliary peak regulation capability for the power system, and gains better income for the power transmission project.

    (4) With the further reduction of energy storage costs in the future, the advantages of the proposed method will become more prominent, providing an effective way for the expansion planning of large-scale wind power bases.

    Acknowledgement:I’d like to express my sincere gratitude to Brother Xu Minghong.I was particularly moved and encouraged by his detailed comments and suggestions when revising my thesis.Without his supervision and guidance, I couldn’t have finished the task.

    Funding Statement:This study was supported by the National Key Research and Development Program (2016YFB0900100)

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

    欧美高清成人免费视频www| 亚洲国产欧美人成| 黄片wwwwww| 国产真实伦视频高清在线观看| 赤兔流量卡办理| 国产欧美日韩精品一区二区| 午夜福利视频1000在线观看| 好男人视频免费观看在线| 六月丁香七月| 麻豆成人午夜福利视频| 亚洲国产精品久久男人天堂| 久久久久久久国产电影| 国产 一区 欧美 日韩| 国产亚洲av片在线观看秒播厂 | www.av在线官网国产| 国产白丝娇喘喷水9色精品| 美女脱内裤让男人舔精品视频| 日本免费在线观看一区| 亚洲电影在线观看av| 日韩欧美 国产精品| 亚洲乱码一区二区免费版| 国内精品美女久久久久久| 日本黄大片高清| 精品人妻偷拍中文字幕| 婷婷色麻豆天堂久久 | 国产真实乱freesex| 亚洲成人久久爱视频| 国产精品美女特级片免费视频播放器| 韩国av在线不卡| 色视频www国产| 一边摸一边抽搐一进一小说| 综合色丁香网| 国产一区二区在线观看日韩| 色吧在线观看| 亚洲欧美日韩卡通动漫| 欧美成人精品欧美一级黄| 国产午夜精品一二区理论片| 久久久亚洲精品成人影院| 国产精品人妻久久久影院| 欧美成人一区二区免费高清观看| 色尼玛亚洲综合影院| 亚洲成色77777| 在线观看一区二区三区| av国产久精品久网站免费入址| 91精品国产九色| 亚洲,欧美,日韩| 亚洲精品日韩在线中文字幕| av在线天堂中文字幕| 狂野欧美白嫩少妇大欣赏| 久久久久久久亚洲中文字幕| 国产午夜福利久久久久久| 亚洲综合色惰| 欧美激情在线99| 日本黄色片子视频| 少妇高潮的动态图| 久久久久久久久久黄片| 亚洲欧美中文字幕日韩二区| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 超碰av人人做人人爽久久| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 男女啪啪激烈高潮av片| 青青草视频在线视频观看| 三级国产精品片| 亚洲欧美成人精品一区二区| 免费观看a级毛片全部| 日本免费a在线| 久久久国产成人精品二区| 麻豆成人av视频| 中文欧美无线码| 一级二级三级毛片免费看| 国产中年淑女户外野战色| 中文字幕av在线有码专区| 亚洲精品久久久久久婷婷小说 | 精品人妻一区二区三区麻豆| 国产免费一级a男人的天堂| 国产亚洲91精品色在线| 秋霞在线观看毛片| 久久久久免费精品人妻一区二区| eeuss影院久久| 欧美激情在线99| 高清毛片免费看| 亚洲av二区三区四区| 亚洲av成人精品一二三区| av卡一久久| 国产精品永久免费网站| 老司机影院毛片| 国产精品久久久久久精品电影| 97热精品久久久久久| 在线播放国产精品三级| 国产精品一区二区性色av| 国产av不卡久久| 在线免费十八禁| 日韩一区二区视频免费看| 精品人妻偷拍中文字幕| 免费观看性生交大片5| 成人亚洲精品av一区二区| 伦理电影大哥的女人| 热99在线观看视频| 国内精品宾馆在线| 国产亚洲精品av在线| 亚洲精品日韩av片在线观看| 久久国内精品自在自线图片| 亚洲欧美成人精品一区二区| 九九久久精品国产亚洲av麻豆| 久久午夜福利片| 乱码一卡2卡4卡精品| 国产午夜福利久久久久久| 成年女人永久免费观看视频| 久久精品国产99精品国产亚洲性色| 亚洲在久久综合| 色噜噜av男人的天堂激情| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄 | 大香蕉久久网| 国产国拍精品亚洲av在线观看| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 日本午夜av视频| 卡戴珊不雅视频在线播放| 久久精品国产鲁丝片午夜精品| 成年女人看的毛片在线观看| 床上黄色一级片| 午夜激情欧美在线| 成年版毛片免费区| 高清av免费在线| 51国产日韩欧美| 日日摸夜夜添夜夜添av毛片| 亚洲怡红院男人天堂| 亚洲欧美清纯卡通| 中国美白少妇内射xxxbb| 日本色播在线视频| 亚洲国产欧洲综合997久久,| 少妇猛男粗大的猛烈进出视频 | 午夜福利成人在线免费观看| 免费观看在线日韩| 偷拍熟女少妇极品色| 嘟嘟电影网在线观看| 看十八女毛片水多多多| 国产一区二区亚洲精品在线观看| 亚洲av电影不卡..在线观看| 一级黄色大片毛片| 激情 狠狠 欧美| 91精品国产九色| 小说图片视频综合网站| 国产精品女同一区二区软件| 午夜精品国产一区二区电影 | 免费观看精品视频网站| 国产一区二区三区av在线| 免费搜索国产男女视频| 免费一级毛片在线播放高清视频| 国产不卡一卡二| 插阴视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 日韩制服骚丝袜av| 国产中年淑女户外野战色| 大香蕉久久网| 91久久精品国产一区二区三区| 好男人在线观看高清免费视频| 中文乱码字字幕精品一区二区三区 | 精品久久久久久电影网 | 国产成人精品一,二区| 波多野结衣高清无吗| 天美传媒精品一区二区| 亚洲在线自拍视频| 中文欧美无线码| 国模一区二区三区四区视频| 国产精品不卡视频一区二区| 国产一区亚洲一区在线观看| 极品教师在线视频| 国产一级毛片七仙女欲春2| 男女国产视频网站| 又黄又爽又刺激的免费视频.| 午夜福利高清视频| 久久午夜福利片| 国产激情偷乱视频一区二区| www日本黄色视频网| 午夜精品在线福利| 大香蕉久久网| 我要看日韩黄色一级片| 天堂网av新在线| 国产av在哪里看| 性插视频无遮挡在线免费观看| 成人综合一区亚洲| 男女下面进入的视频免费午夜| 亚洲四区av| 在线播放无遮挡| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 在线a可以看的网站| 综合色av麻豆| 午夜久久久久精精品| 日韩av在线免费看完整版不卡| 天美传媒精品一区二区| 欧美97在线视频| 美女大奶头视频| 男人的好看免费观看在线视频| 成人一区二区视频在线观看| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类| 在线天堂最新版资源| 欧美性猛交╳xxx乱大交人| av线在线观看网站| 2021天堂中文幕一二区在线观| 嫩草影院入口| 国产伦在线观看视频一区| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的 | 五月伊人婷婷丁香| 亚洲欧美中文字幕日韩二区| 99视频精品全部免费 在线| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 欧美精品一区二区大全| 97热精品久久久久久| 看片在线看免费视频| 久久久亚洲精品成人影院| 亚洲成人av在线免费| 中文字幕亚洲精品专区| 亚洲国产日韩欧美精品在线观看| 成人亚洲欧美一区二区av| 七月丁香在线播放| 毛片女人毛片| 午夜精品一区二区三区免费看| 成人毛片60女人毛片免费| 久久久久久大精品| 亚洲人与动物交配视频| 国产伦一二天堂av在线观看| 蜜臀久久99精品久久宅男| 99九九线精品视频在线观看视频| 91精品国产九色| 九九在线视频观看精品| 色哟哟·www| 日韩一区二区三区影片| 亚洲成av人片在线播放无| 中文欧美无线码| 岛国在线免费视频观看| 18禁动态无遮挡网站| 色5月婷婷丁香| 又粗又爽又猛毛片免费看| 久久久久久大精品| 超碰97精品在线观看| h日本视频在线播放| 亚洲美女视频黄频| 亚洲国产最新在线播放| 国产单亲对白刺激| 岛国毛片在线播放| 少妇被粗大猛烈的视频| av卡一久久| 久99久视频精品免费| 亚洲av中文字字幕乱码综合| 久久久久久九九精品二区国产| 日韩强制内射视频| 欧美3d第一页| 久久久久国产网址| 色综合亚洲欧美另类图片| 禁无遮挡网站| 精品久久久噜噜| 午夜福利在线观看免费完整高清在| 久久久久久久久久黄片| 日韩 亚洲 欧美在线| 特级一级黄色大片| 精品久久久久久久久亚洲| 少妇猛男粗大的猛烈进出视频 | 午夜a级毛片| 亚洲国产最新在线播放| 亚洲在线观看片| 国产午夜福利久久久久久| 欧美一区二区国产精品久久精品| 免费看av在线观看网站| 免费看光身美女| 亚洲怡红院男人天堂| 午夜免费激情av| 又黄又爽又刺激的免费视频.| 2022亚洲国产成人精品| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 欧美三级亚洲精品| 亚洲av成人av| 国产片特级美女逼逼视频| 久久这里有精品视频免费| ponron亚洲| 国内揄拍国产精品人妻在线| 欧美zozozo另类| 亚洲av电影不卡..在线观看| 亚洲av日韩在线播放| 激情 狠狠 欧美| 成年av动漫网址| 夜夜看夜夜爽夜夜摸| 一边亲一边摸免费视频| 亚洲欧美中文字幕日韩二区| 久久久久久久久中文| 精品人妻偷拍中文字幕| 狠狠狠狠99中文字幕| 成人二区视频| 人妻夜夜爽99麻豆av| 床上黄色一级片| 欧美极品一区二区三区四区| 精品不卡国产一区二区三区| 国国产精品蜜臀av免费| 成人亚洲精品av一区二区| 国产精品不卡视频一区二区| 久久99精品国语久久久| 免费观看性生交大片5| 九九爱精品视频在线观看| 欧美三级亚洲精品| 欧美精品国产亚洲| 国产精品久久久久久精品电影小说 | 欧美一级a爱片免费观看看| 国产av在哪里看| 国产精品久久久久久av不卡| 麻豆精品久久久久久蜜桃| 亚洲在线观看片| 国产亚洲5aaaaa淫片| 国产精品爽爽va在线观看网站| 在现免费观看毛片| 精品99又大又爽又粗少妇毛片| 欧美不卡视频在线免费观看| 成人国产麻豆网| 午夜免费激情av| 国产成人精品一,二区| 三级经典国产精品| 中文字幕av成人在线电影| 乱人视频在线观看| 免费看av在线观看网站| 人人妻人人看人人澡| 2022亚洲国产成人精品| 亚洲一级一片aⅴ在线观看| 卡戴珊不雅视频在线播放| 国产视频内射| 真实男女啪啪啪动态图| 国产精品人妻久久久久久| 国产精品综合久久久久久久免费| 99热网站在线观看| 插阴视频在线观看视频| 国产乱来视频区| 国产精品人妻久久久久久| 少妇人妻一区二区三区视频| 免费观看人在逋| 国产亚洲av片在线观看秒播厂 | 亚洲精品456在线播放app| 亚洲精华国产精华液的使用体验| 免费观看a级毛片全部| 51国产日韩欧美| 日本免费一区二区三区高清不卡| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 成人av在线播放网站| 免费观看在线日韩| 男人和女人高潮做爰伦理| 午夜福利在线观看免费完整高清在| 国产在视频线精品| 永久网站在线| 老司机影院成人| 寂寞人妻少妇视频99o| 一本久久精品| 精品一区二区三区人妻视频| 能在线免费观看的黄片| 男女啪啪激烈高潮av片| 韩国av在线不卡| 日本免费在线观看一区| 国产黄片视频在线免费观看| 一级av片app| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲伊人久久精品综合 | 久久精品久久精品一区二区三区| 蜜桃久久精品国产亚洲av| 边亲边吃奶的免费视频| 久久6这里有精品| 精品久久久久久久人妻蜜臀av| 草草在线视频免费看| 韩国av在线不卡| 亚洲中文字幕日韩| 观看美女的网站| 日本一本二区三区精品| 国产一区二区三区av在线| 国产片特级美女逼逼视频| 少妇高潮的动态图| 国产亚洲最大av| 99久久人妻综合| 久久欧美精品欧美久久欧美| 亚洲欧美清纯卡通| 熟妇人妻久久中文字幕3abv| 亚洲av免费在线观看| 干丝袜人妻中文字幕| 91aial.com中文字幕在线观看| 国产视频内射| 久久精品国产自在天天线| 精品酒店卫生间| 一二三四中文在线观看免费高清| 小说图片视频综合网站| 99久国产av精品国产电影| 三级国产精品欧美在线观看| 久99久视频精品免费| 日本免费一区二区三区高清不卡| 乱系列少妇在线播放| 亚洲精品亚洲一区二区| 婷婷色av中文字幕| 国产不卡一卡二| 99久久人妻综合| 日韩大片免费观看网站 | av播播在线观看一区| 最近2019中文字幕mv第一页| 亚洲成人av在线免费| 18禁在线播放成人免费| 日韩欧美国产在线观看| 国产av在哪里看| 国产视频首页在线观看| 久久精品久久久久久久性| 亚洲三级黄色毛片| 亚洲激情五月婷婷啪啪| 97超视频在线观看视频| 国产精品美女特级片免费视频播放器| 免费看光身美女| 久久久久久久久久黄片| 18禁动态无遮挡网站| av天堂中文字幕网| 精品久久久久久久久av| av女优亚洲男人天堂| 免费av观看视频| 国内揄拍国产精品人妻在线| 在线免费十八禁| 特级一级黄色大片| 男女下面进入的视频免费午夜| 日本一二三区视频观看| 色综合色国产| 69人妻影院| 久久精品国产亚洲av天美| 99视频精品全部免费 在线| 亚洲熟妇中文字幕五十中出| 中文乱码字字幕精品一区二区三区 | 日日摸夜夜添夜夜爱| 婷婷色av中文字幕| 亚洲,欧美,日韩| 久久精品影院6| 国产综合懂色| av免费在线看不卡| 国产综合懂色| av免费在线看不卡| 国产在线一区二区三区精 | 亚洲欧美清纯卡通| 天天一区二区日本电影三级| 亚洲av免费在线观看| 高清在线视频一区二区三区 | 久久久国产成人免费| 特级一级黄色大片| 久久婷婷人人爽人人干人人爱| 国内精品宾馆在线| 亚洲人成网站在线观看播放| 2021天堂中文幕一二区在线观| 中文精品一卡2卡3卡4更新| 久久国内精品自在自线图片| 波多野结衣巨乳人妻| 免费看a级黄色片| 一二三四中文在线观看免费高清| 国产精品伦人一区二区| 亚洲av二区三区四区| 久久久精品欧美日韩精品| a级一级毛片免费在线观看| 国产成人精品婷婷| 国产乱来视频区| 岛国毛片在线播放| 97在线视频观看| 国产探花在线观看一区二区| 99热网站在线观看| 国产精品一区www在线观看| 国产在线男女| 午夜a级毛片| 国产麻豆成人av免费视频| 一个人看视频在线观看www免费| 日韩视频在线欧美| 精品不卡国产一区二区三区| 噜噜噜噜噜久久久久久91| 亚洲图色成人| 只有这里有精品99| 国产亚洲av嫩草精品影院| 不卡视频在线观看欧美| 五月玫瑰六月丁香| 午夜精品国产一区二区电影 | 天堂√8在线中文| 久久久国产成人精品二区| 国产亚洲5aaaaa淫片| 人妻少妇偷人精品九色| 国产白丝娇喘喷水9色精品| 青春草亚洲视频在线观看| 成人一区二区视频在线观看| 美女黄网站色视频| videos熟女内射| 国产不卡一卡二| 成人av在线播放网站| 欧美一区二区精品小视频在线| 99久国产av精品国产电影| 在线天堂最新版资源| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| 免费观看人在逋| 一个人看视频在线观看www免费| 一二三四中文在线观看免费高清| 国内精品一区二区在线观看| 嫩草影院新地址| 欧美人与善性xxx| 国产精品精品国产色婷婷| 人人妻人人看人人澡| 噜噜噜噜噜久久久久久91| 人人妻人人看人人澡| 国产精品爽爽va在线观看网站| 看黄色毛片网站| 天堂√8在线中文| 欧美潮喷喷水| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 国产亚洲最大av| 亚洲色图av天堂| 成人av在线播放网站| 婷婷色综合大香蕉| 国产淫语在线视频| 免费观看性生交大片5| 日本黄大片高清| 蜜桃亚洲精品一区二区三区| 亚州av有码| 亚洲国产日韩欧美精品在线观看| 国语对白做爰xxxⅹ性视频网站| 汤姆久久久久久久影院中文字幕 | 国产免费又黄又爽又色| 69人妻影院| av在线蜜桃| 国产单亲对白刺激| 亚洲国产欧洲综合997久久,| 不卡视频在线观看欧美| 国产老妇伦熟女老妇高清| 中文欧美无线码| 亚洲人成网站高清观看| 最近中文字幕高清免费大全6| av.在线天堂| 成人二区视频| 国产亚洲91精品色在线| 国产精品国产高清国产av| 午夜爱爱视频在线播放| 色噜噜av男人的天堂激情| 大又大粗又爽又黄少妇毛片口| 两个人视频免费观看高清| 欧美精品一区二区大全| 丝袜喷水一区| 我要看日韩黄色一级片| 欧美日韩一区二区视频在线观看视频在线 | 不卡视频在线观看欧美| 久久久国产成人免费| 26uuu在线亚洲综合色| 精品国内亚洲2022精品成人| 亚洲av熟女| av在线蜜桃| 极品教师在线视频| 国产av不卡久久| 看免费成人av毛片| 亚洲欧美一区二区三区国产| 国产三级在线视频| 亚洲久久久久久中文字幕| 久久久精品94久久精品| 99视频精品全部免费 在线| 国产69精品久久久久777片| 免费电影在线观看免费观看| 亚洲丝袜综合中文字幕| 成人特级av手机在线观看| 特级一级黄色大片| 欧美日韩综合久久久久久| 国产精品人妻久久久影院| 欧美变态另类bdsm刘玥| 国产黄色小视频在线观看| 九草在线视频观看| 午夜福利在线观看吧| 国产免费福利视频在线观看| 热99re8久久精品国产| 成人毛片a级毛片在线播放| 一级黄片播放器| 人人妻人人澡欧美一区二区| 国产精品一二三区在线看| 哪个播放器可以免费观看大片| 欧美成人午夜免费资源| 我的女老师完整版在线观看| 久久久a久久爽久久v久久| 国产精品三级大全| 麻豆成人av视频| 性插视频无遮挡在线免费观看| 免费观看在线日韩| 三级毛片av免费| 美女大奶头视频| 级片在线观看| 精品久久国产蜜桃| 两性午夜刺激爽爽歪歪视频在线观看| 美女内射精品一级片tv| 天堂中文最新版在线下载 | 国产免费男女视频| 丰满少妇做爰视频| 日本黄大片高清| 99热这里只有是精品在线观看| 国内揄拍国产精品人妻在线| 国产精品麻豆人妻色哟哟久久 | 嫩草影院入口| 亚洲成人精品中文字幕电影| 欧美潮喷喷水| 18禁裸乳无遮挡免费网站照片| 亚洲精品成人久久久久久| 久久久欧美国产精品| 一本一本综合久久| 亚洲怡红院男人天堂| 搞女人的毛片| 最近的中文字幕免费完整| 精品久久久久久电影网 | 99久久精品国产国产毛片| 午夜a级毛片| 伦精品一区二区三区| 国产美女午夜福利|