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

    Dynamic propagation velocity of a positive streamer in a 3 m air gap under lightning impulse voltage

    2024-04-24 06:17:02ZhiweiLI李志偉TingLEI雷挺YuSU蘇宇XiuyuanYAO姚修遠(yuǎn)BingxueYANG楊冰雪DelongLIU劉德龍FangchengLV律方成andYujianDING丁玉劍
    Plasma Science and Technology 2024年4期
    關(guān)鍵詞:方成冰雪

    Zhiwei LI (李志偉),Ting LEI (雷挺),Yu SU (蘇宇),Xiuyuan YAO (姚修遠(yuǎn)),Bingxue YANG (楊冰雪),Delong LIU (劉德龍),Fangcheng LV (律方成) and Yujian DING (丁玉劍),*

    1 North China Electric Power University,Beijing 102206,People’s Republic of China

    2 China Electric Power Research Institute,Beijing 100192,People’s Republic of China

    Abstract Streamers represent an important stage in the initiation of gap discharge.In this work,we used an eight-frame intensified charge-coupled device camera to capture the streamer development process when a lightning impulse voltage of 95%–100% U50% was applied in a 3 m rod–plate gap and the streamer velocity was analyzed.Analysis of the observations shows that streamer velocity can be defined by three stages: rapid velocity decline (stage 1),rapid velocity rise (stage 2)and slow velocity decline (stage 3).The effects of electrode shape,applied voltage and gap breakdown or withstanding on streamer velocity were analyzed.The electrode with a larger radius of curvature will result in a higher initial velocity,and a higher voltage amplitude will cause the streamer to propagate faster at stage 3.Gap withstanding or breakdown has no obvious effect on streamer velocity.In addition,the experimental results are compared with previous results and the statistical characteristics of the primary streamer discharge are discussed.

    Keywords: lightning impulse,3 m air gap,development process,streamer velocity

    1.Introduction

    Discharge is a key and basic issue in the fields of high-voltage engineering and plasma research.The generation and development of streamers is an important stage in the gap discharge process [1–5].In-depth research into the microscopic mechanism of streamer discharge is helpful for improving the theory of gas discharge and reducing the thickness of the insulation during the design of power transmission and transformation equipment [6,7].Velocity is an important parameter in the development of streamers,so scholars have analyzed it by different methods.

    In 1972,Les Renardières’ group reported some interesting findings concerning the developmental properties of streamers [8].They estimated the velocity of a streamer by analyzing the results from a streak camera and found that the streamer starts at a high velocity which then decreases because the local field decreases further down the gap.However,due to the limitations of the experimental conditions,only the trend of streamer velocity development was obtained and the estimated results are not accurate.

    Grangeet al[9] measured the streamer velocity in a 1 cm rod–plate gap using a photomultiplier tube.When the applied voltage was 8 kV,the velocity decreased from the initial 0.5 mμs–1to 0.2 mμs–1with the development of the streamer.Wanget al[10] used a Blumlein pulse generator and a streak camera to study the development of positive and negative streamer discharges.When a positive voltage of 43–90 kV was applied,the streamer velocity was 0.8–1.2 mμs–1;when a negative voltage of 96 kV was applied,the streamer velocity was about 0.6 mμs–1.The accuracy of the above research results has been improved,but there is a lack of research under conditions of long gaps and multiple influencing factors.

    In recent years,more and more researchers have adopted charge-coupled devices (CCDs) and intensified chargecoupled devices (ICCDs) to carry out related studies and obtain more accurate parameters.Wanget al[11] observed the development of streamer discharge in a wire-cylinder electrode using an ICCD camera.It was found that the maximum propagation velocity of the streamer discharge was 1.8–3.3 mμs–1at 72–91 kV and the field strength on the wire surface significantly affected the propagation velocity.Brielset al[12] used a CCD camera to observe the streamer in a point–plate air gap of length 4 cm and found that no streamers were generated for an applied voltage of less than 5 kV.The streamer velocity increased from 0.1 mμs–1to 4 mμs–1for an applied voltage of 5–96 kV.Moreover,the streamer velocity was greatly affected by voltage,voltage rise time and polarity.

    Zeng and Chen [13] used 3a four-frame ICCD camera to study the dynamic velocity of a streamer in a 57 cm air gap under lightning impulse voltage.When the voltage amplitude was 250 kV,the streamer velocity decreased from 8.3 ±2.7 mμs–1to 0.4 ± 0.1 mμs–1under a hemispherical electrode.For different electrode shapes,the streamer velocity wasv(spherical)>v(hemispherical)>v(conical).Moreover,the streamer velocity increased with increasing rate of the voltage rise.References [14–19] have also demonstrated important progress in the simulation of the effects of applied voltage,background electric field,voltage polarity,electrode shape and experimental environment on streamer velocity.However,due to the restriction of experimental conditions,most studies on streamer discharge under lightning impulse were carried out in short air gaps with a length of several centimeters to tens of centimeters.There have been fewer experiments on streamer development in meter-scale long air gaps.Also,the scope of observations and studies cannot cover the whole process of streamer development.

    In this work,an eight-frame ICCD camera was used to observe the dynamic development characteristics of a streamer in a 3 m rod–plate gap under lightning impulse.The influences of electrode shape,applied voltage and gap withstanding or breakdown on the streamer velocity are analyzed and compared with previous results.The relationships between the inception voltage and the peak current of the primary streamer discharge and the peak value of the impulse voltage are discussed.

    2.Experimental setup and analytical method

    2.1. Experimental setup

    The experiment site is in Beijing (55 m above sea level),and a schematic diagram of the experimental layout is shown in figure 1.A positive lightning impulse voltage with a wavefront time of 1.3μs is generated by a 20-stage 6000 kV/900 kJ Marx generator.The rod–plate gap is 3 m in length.The hollow rod electrode is 1.5 m long with a measuring device placed inside,and the outer metal casing is insulated from the end electrode tip by an insulating layer to prevent the displacement current generated by the shielded casing from flowing through the measuring device via the electrode tip [20].The plate electrode is a 2 m × 2 m piece of galvanized iron sheeting.Two types of electrode tips were applied in the experiments,as shown in figure 1(b).One is a conical electrode tip with a bottom cross-section radius of 10 mm and a tip radius of 0.25 mm and the other is a hemispherical electrode tip with a radius of 10 mm.Under such experimental conditions,we haveU50%=1670 kV.

    The image recording device of the experimental observation system is a SIMX8 ICCD camera.The SIMX8 has eight image sensor channels and can take eight photos in a row with a maximum image resolution of 1280 × 960 pixels.The interval between two adjacent frames can be set from 0 ns to 20 ms and the exposure time can be set from 3 ns to 10 ms.The above two parameters are independently variable between two adjacent frames and the minimum setting amplitude is 1 ns,so we can observe the dynamic development of the streamer more accurately.The voltage was measured by a capacitive voltage divider,and the voltage signal was captured by an oscilloscope.The current recording device,embedded in the rod electrode,mainly comprises a 2 Ω non-inductive resistor and a voltage recorder [20].The sample rate of the current measuring device was 100 MHz.When a lightning impulse current flows through the noninductive resistor,the voltage recorder collects the voltage signals at both ends of the non-inductive resistor and transmits them to the controller in the control room through the optical fiber.

    When the oscilloscope receives the voltage signal,it immediately sends a trigger signal to the synchronous trigger DG645,which receives the trigger signal and synchronously triggers the image and current acquisition device.Meanwhile,the time delay of the cable,photoelectric conversion and other equipment that transmits the trigger signals should be taken into account.After synchronous calibration,the recording time errors of the image,voltage and current signals were controlled within 20 ns.As shown in figure 2(a),the synchronous observation system starts at time 0 and a current pulse appears at the instantI1,which is caused by the movement of charge carriers in the space at the beginning of the streamer discharge.This instant is identified as the initial time of the streamer.After a preset starting time delay,the ICCD starts shooting at the instantI2.Various color bands in the picture correspond to the shooting instants and exposure time of the ICCD.

    2.2. Methods for analyzing streamer velocity

    When the gap discharges,the atomic energy level transition in the streamer space produces photons.By utilizing an ICCD for shooting and setting different starting shooting instants,exposure times and shooting intervals for the experiments,and then integrating the results based on the development time,more comprehensive spatial motion characteristics of streamer development can be obtained.Figure 2(b)shows the projection of the streamer captured by an ICCD on a two-dimensional plane.The velocityvis calculated as the axial development distance Δyof the streamer between two adjacent frames divided by the time intervaltibetween two adjacent frames

    It is worth noting that the exposure time of the ICCD in this experiment is mostly set at tens of nanoseconds,whereas the entire streamer development process takes about 3μs.So,it is necessary to adjust the ICCD’s settings in order to allow the calculated value to cover all stages of streamer development.The advantages of such a setting are that the trend of computing results will be more obvious and have a wider coverage of the timeline.However,the disadvantage is that the computational results are more scattered.At the same time,the wider gap in this experiment compared with previous ones will increase the dispersion of the streamer’s initial time,which also leads to a greater dispersion of the computed velocity results.The experimental results are therefore mostly described in scatter charts.

    3.Experiments and results

    3.1. Dynamic propagation velocity of streamers

    Figure 3 shows the typical development of a streamer at different spatial ranges of observation and different time scales when 95%U50%is applied to the conical electrode.The streamer size is small during the initial stage of propagation.In order to improve the shooting accuracy at this stage,the focal length of the ICCD was adjusted to set the spatial range of observation to 0.4 m.In the later stage the streamer was about to penetrate the gap,and the experiment was repeated after the spatial range of observation was set to 3 m.

    The relationship between streamer velocity and time is shown in figure 4.The figure is the result of integration of multiple sets of data under the same experimental conditions and different shooting modes.In our research,the streamer velocity is defined as having three stages.In stage 1 (initiation–0.88μs),the streamer velocity develops at a high initial value and then decreases sharply.The axial propagation velocity reaches the lowest value whent=0.88μs.The inset in the upper right of figure 4 shows the result in this stage zoomed in along the timeline.As the streamer develops very quickly during the first stage and is difficult to capture,it is not possible to obtain a stable initial velocity.In stage 2 (0.88–1.30μs) of streamer development,the streamer velocity increases sharply.Att=1.30μs,the velocity peaks at about 1.6 mμs–1.In stage 3 (1.30–3μs),the streamer velocity decreases slowly until it penetrates the gap.It takes the streamer about 3μs to develop onto the plate electrode.

    Throughout the streamer development process,the velocity shows a trend of ‘declining rapidly,rising rapidly and then declining slowly’.Specifically,streamer development is fastest at the beginning and its velocity reaches its lowest value at the boundary of stage 1 and stage 2 (t=0.88μs).The streamer has the longest propagation time during stage 3,about 1.7μs.

    Next,the above results are analyzed.Figure 5 shows the different stages of development of a streamer.As the streamer discharge region increases during stage 1,the electric field distribution gradually becomes uniform and the electric field in the streamer development area decreases,resulting in a rapid decrease in the streamer velocity.In stage 2,the rising voltage generates a high electric field around the anode and the secondary streamer discharge starts at a high initial velocity and continues to decrease until it converges with the primary streamer [21].During this stage,initiation of the secondary streamer enhances the electric field in the primary streamer head.Meanwhile,the increased electric field accelerates the propagation of the primary streamer,resulting in an increase in the measured value of the streamer velocity.In stage 3,the secondary streamer completely covers the primary streamer.Thereafter,there is no difference between the primary streamer and secondary streamer.The electric field distribution gradually becomes uniform and the streamer velocity gradually decreases.

    3.2. Streamer velocity for different electrode shapes

    As shown in figure 6,an interesting phenomenon was found when using different electrodes in the experiments.Different electrodes will lead to different streamer morphologies.In the experiment,a lightning impulse voltage of 100%U50%was applied.The exposure time was 20 ns at an observation space range of 0.3 m and 500 ns at an observation space range of 3 m.

    It can be seen from the figure that the entire streamer discharge area under the two electrodes comprises a sector or a pyramid shape.In the early stage of development,the streamer under the hemispherical electrode with a bigger curvature radius has obvious branches,and its light trajectory has a needle shape.The streamer under the conical electrode with a smaller curvature radius has no obvious branches,and its light trajectory is in a dispersed state.At a late stage of development,the streamer morphology under both electrodes is dispersed and shows no significant difference within a 3 m observation space range.We will focus on those parts where there are differences in the initial stages of the streamer under the two electrodes.

    Figure 7 shows the refined measurement results for the initial stage of a streamer under both electrodes.The figure integrates the refined measurement results of the initial stage velocity of the streamer under the same experimental conditions and different shooting modes.From the figure,it can be found that the streamer velocity under the hemispherical electrode starts rapidly and decreases rapidly;the streamer velocity under the conical electrode starts slowly and decreases slowly.Att=0.44μs,the streamer velocity is the same under both hemispherical and conical electrodes.

    3.3. Streamer velocity in the case of gap withstanding and breakdown

    Figure 8 shows the variation of streamer length and velocity with time for gap withstanding and breakdown under the same conditions att=0.4–1.8μs.The figure shows the results of a single shot,where the linear plot is the length of streamer development and the scatter plot is the streamer velocity.It can be observed that the values and trends for streamer development length and velocity do not differ much in the two cases.

    Att=1.196μs,the difference in development length is the largest at about 5.2 cm;att=1.118μs,the velocity difference is the largest at about 0.705 mμs–1.Therefore,it can be inferred that there is no significant effect on streamer velocity for gap discharge withstanding or breakdown.

    3.4. Streamer velocity under different voltages

    Figure 9 shows the variation of streamer velocity when different voltages are applied.The figure contains the integrated results of multiple sets of data under the same experimental conditions and different shooting modes.The difference between the velocity of the streamer when 95%–100%U50%was applied is not significant in stages 1 and 2.In stage 3,the higher the applied voltage,the faster the streamer develops.Att=2.25μs,the velocity under 100%U50%is about 0.25 mμs–1higher than that under 95%U50%.

    4.Discussion

    4.1. Comparison with previous findings

    In this section,our experimental results are compared with the typical results of previous work.Morrow and Lowke [22]obtained a positive streamer development law in a 50 mm point–plane gap at an applied voltage of 20 kV by means of simulation.Based on their simulation results,the streamer velocity was divided into three stages: rapid streamer propagation (6–40 ns),slow streamer propagation (40–130 ns) and streamer termination (130–210 ns).In our experiment,the streamer velocity is defined as comprising three stages based on measured and calculated results: rapid velocity decline(stage 1,initiation–0.88μs),rapid velocity rise (stage 2,0.88–1.30μs) and slow velocity decline (stage 3,1.30–3.00μs).Compared with the simulation results of Morrow and Lowke,our experiment uses a high-precision ICCD to observe the process of streamer velocity recovery when a positive lightning impulse voltage is applied.

    Zeng and Chen [13] studied streamer discharge development in a 0.57 m rod–plate gap at an experimental voltage amplitude of 250 kV.For hemispherical and conical electrodes,the streamer velocity decreased rapidly from 3.2 mμs–1to 0.80 mμs–1in 50 ns and from 1.7 mμs–1to 0.5 mμs–1in 60 ns,respectively.In our experiment,the process of development of a 3 m rod–plate gap streamer was studied.When a voltage of 1670 kV was applied to the hemispherical and conical electrodes,the streamer velocity decreased rapidly from 4.8 mμs–1to 0.80 mμs–1in 0.25μs and from 3.5 mμs–1to 0.5 mμs–1in 0.4μs,respectively.The results of this study and previous studies show that,compared with a conical electrode,a hemispherical electrode has a faster initial velocity and a greater falling rate.

    Brielset al[12] at TU/e University studied streamer development at a 0.04 m point–plane gap.The applied voltage was increased from 5 kV to 96 kV,and the streamer velocity increased from 0.1 mμs–1to 4 mμs–1.In Creyghtonet al’s experiment [23],the streamer velocity increased from 2 mμs–1to 3.5 mμs–1when a pulse voltage of 20–25 kV was applied.The present experimental measurements found that the higher the applied voltage at stage 3 the faster the streamer velocity when applying 95%–100%U50%.Hence,a higher voltage amplitude leads to a greater velocity.References [24–26] also obtained the same results through experiments and simulations.

    The voltage rise rate is another important factor affecting streamer velocity.Komuroet al[27] simulated that when the applied voltage of a 13 mm rod–plate gap was 24 kV and the voltage rise rates were 0.11 kV ns–1and 0.52 kV ns–1,the average streamer velocities were 0.49 mμs–1and 0.69 mμs–1,respectively.Winandset al[28] found that when the voltage rise rate increased in the range of 1.5–2 kV ns–1,the streamer velocity increased from 0.5 mμs–1to 2.5 mμs–1.In our experiment,the voltage rise rates were 1.22 kV ns–1and 1.28 kV ns–1,and the average streamer development velocities were 0.9 mμs–1and 1 mμs–1,respectively.Yagiet al[29]and Yoshinagaet al[30] also did a similar study and concluded that the propagation velocity of the streamer heads increased with a higher rise rate of the applied voltage.

    It should be noted that although the initial streamer velocity in this experiment is large,this value is underestimated due to the limitations of the experimental equipment and methodology.In future research,we aim to extend the experiment to extreme environments such as high altitude and low temperature.At the same time,we will explore the universal model of streamer propagation based on the existing data and known influencing factors to obtain a more accurate picture of dynamic streamer development.

    4.2. Statistical characteristics of the primary streamer discharge

    In this work,the statistical characteristics of the primary streamer discharge were analyzed for lightning impulse voltages of 90%U50%,95%U50%and 100%U50%(U50%=1670 kV;the experiment applying 90%U50%was a supplementary experiment).Figures 10–12 show the integrated results for multiple sets of data under the same experimental conditions.

    Figure 10 shows the relationship betweenui,the inception voltage of the primary streamer discharge,andthe peak value of the impulse voltage.In this figure,the data points represent the results of each experiment and the line segments represent the linear fitting values forui.The statistical results were aggregated into three groups.Regardless of the voltage range,uishows a dispersion

    Figure 1.(a) Simultaneous observation systems.(b) Conical electrode tip and hemispherical electrode tip.

    Figure 2.(a) Voltage and current waveforms and ICCD operating states in a 3 m rod–plane gap at 100% U50% (synchronous observation system starting time 0 μs,ICCD starting time 0.396 μs,exposure time 0.16 μs).(b) The second,fourth,sixth and eighth images captured by the ICCD in figure 2(a).

    Figure 3.Typical development of a streamer under different observation spaces and time scales when 95% U50% is applied to a cone electrode: (a) observation space 0.4 m,exposure time 30 ns;(b) observation space 1 m,exposure time 50 ns;(c) observation space 3 m,exposure time 160 ns;(d) observation space 3 m,exposure time 160 ns.The shooting interval is 160 ns.

    Figure 4.Dynamic development of streamer velocity.The time zero is the start time of the synchronous observation system.

    Figure 5.The different stages of development of a streamer.

    Figure 6.Comparison of streamer development morphology for different electrode shapes.The applied voltage is 100% U50%: (a)hemispherical electrode and (b) conical electrode.

    Figure 8.Variation of streamer length and velocity with time in the case of gap withstanding and breakdown (cone electrode,95%U50%).

    Figure 9.Variation of streamer velocity when 100% U50% and 95%U50% are applied to a cone electrode.

    Figure 10.Relationship between the inception voltage of the primary streamer discharge and the peak value of the impulse voltage (correlation coefficient R2=0.036).

    It can be seen from the statistical results that the discharge inception voltageuiof the primary streamer and the deviation of inception voltage δuiremain approximately constant;the average inception voltageis almost independent of the peak value of the impulse voltageWith increase in the statistical samples,shows minor statistical changes.

    Figure 11 shows the relationship between the peak currentipof the primary streamer and the peak value of impulse voltage.Although there is still dispersion,it is obvious thatipincreases with increase in

    Figure 11.Relationship between the peak current of the primary streamer discharge and the peak value of the impulse voltage(correlation coefficient R2=0.34).

    Equation (3) is the fitting result of a linear function.Due to the influence of data dispersion,it is difficult to obtain definite values fork1andk2.Statistical studies have found that the value range ofk1is [0.004881,0.006942] and that ofk2is [-5.392,-2.017].The fitted curve in figure 12 uses the optimal values of both,withk1=0.005912,k2=-3.705.

    Figure 12 shows the relationship between the peak current valueipand the discharge inception voltageuiof the primary streamer.It can be seen that there is no obvious correspondence betweenipandui.

    Figure 12.Relationship between the peak current value and the inception voltage of the primary streamer discharge (correlation coefficient R2=0.014).

    5.Conclusions

    In this paper,an experimental study of positive streamer discharge in a 3 m rod–plate gap under lightning impulse voltage is described.Based on an ICCD synchronous observation system,the dynamical development of a streamer in the range of 0.3–3.0 m is obtained.The following conclusions are obtained from analyzing the experimental results:

    (1) When 95%–100%U50%is applied to the two types of electrodes,the streamer velocity shows three stages: rapid velocity decline (stage 1),rapid velocity rise (stage 2) and slow velocity decline (stage 3).

    (2) The influence of the electrode on streamer development is studied.When a hemispherical electrode is used,the streamer form gradually evolves from a needle shape at the initial stage to a state of dispersion.When a conical electrode is used,the streamer shows a state of dispersion throughout the development process.The initial streamer velocity of the hemispherical electrode is fast and the falling rate is high,while the initial streamer velocity of the conical electrode is slow and the falling rate is low.

    (3) Gap withstanding or breakdown has no obvious effect on streamer development.In addition,when a 100%U50%positive lightning impulse voltage is applied,the streamer velocity is faster than that when 95%U50%is applied in stage 3,but the difference is not obvious in stages 1 and 2.

    (4) The statistical characteristics of the primary streamer discharge are discussed.There is no obvious correspondence between the inception voltage and the peak current of the primary streamer discharge or the peak impulse voltage.The peak current of the primary streamer discharge is positively correlated with the peak impulse voltage.

    Acknowledgments

    This work was supported by the Beijing Science Fund for Distinguished Young Scholars (No.JQ22009) and National Natural Science Foundation of China (No.51977198).

    猜你喜歡
    方成冰雪
    逐夢的“冰雪一代”
    走向世界(2022年3期)2022-04-19 12:39:14
    擁抱冰雪向未來
    走向世界(2022年3期)2022-04-19 12:39:02
    點燃“冰雪”
    走向世界(2022年3期)2022-04-19 12:39:00
    漫畫大師方成
    ——煩心事也能『幽它一默』
    厚德載物 方成大道
    冰雪“十八盤”
    冰雪奇旅
    王萬奇:吃苦耐勞方成大事
    金色年華(2016年7期)2016-02-28 01:39:29
    冰雪幽默來襲
    方成漫畫選
    婷婷色av中文字幕| 亚洲国产精品成人综合色| 精品酒店卫生间| 国产在线一区二区三区精| 欧美成人午夜免费资源| 99热这里只有精品一区| 在线亚洲精品国产二区图片欧美 | 亚洲av中文av极速乱| 一个人看视频在线观看www免费| 欧美潮喷喷水| 日本三级黄在线观看| 热re99久久精品国产66热6| 亚洲自偷自拍三级| 亚洲av成人精品一二三区| 成人黄色视频免费在线看| 欧美日韩在线观看h| 哪个播放器可以免费观看大片| 亚洲欧美日韩东京热| 久久久久国产精品人妻一区二区| 日本爱情动作片www.在线观看| 天天躁日日操中文字幕| 亚洲高清免费不卡视频| 欧美xxxx黑人xx丫x性爽| 97人妻精品一区二区三区麻豆| 在线播放无遮挡| 国产精品伦人一区二区| 国产中年淑女户外野战色| 国产男女内射视频| 久久久亚洲精品成人影院| 少妇高潮的动态图| 又黄又爽又刺激的免费视频.| 国产精品一区www在线观看| 欧美少妇被猛烈插入视频| 久久精品久久久久久噜噜老黄| 高清视频免费观看一区二区| 日本一二三区视频观看| 黑人高潮一二区| 国产老妇女一区| 日韩免费高清中文字幕av| 黄片无遮挡物在线观看| 亚洲av中文av极速乱| 国产极品天堂在线| 久久鲁丝午夜福利片| 午夜亚洲福利在线播放| 欧美日韩精品成人综合77777| 亚洲精品自拍成人| 自拍偷自拍亚洲精品老妇| 3wmmmm亚洲av在线观看| 中文资源天堂在线| 青春草国产在线视频| 女人十人毛片免费观看3o分钟| 久久综合国产亚洲精品| videossex国产| 久久6这里有精品| 边亲边吃奶的免费视频| av专区在线播放| 蜜桃久久精品国产亚洲av| 亚洲国产精品国产精品| 亚洲av成人精品一区久久| 欧美老熟妇乱子伦牲交| 一个人看的www免费观看视频| 人妻 亚洲 视频| av国产精品久久久久影院| 女人十人毛片免费观看3o分钟| 精品久久久久久久末码| 18禁在线无遮挡免费观看视频| 亚洲精品,欧美精品| 五月开心婷婷网| 亚洲丝袜综合中文字幕| a级一级毛片免费在线观看| 国产毛片在线视频| 欧美日韩综合久久久久久| 午夜福利高清视频| 黄片wwwwww| 一边亲一边摸免费视频| 日日摸夜夜添夜夜爱| 97热精品久久久久久| 国产亚洲av嫩草精品影院| 国产精品精品国产色婷婷| 亚洲成人中文字幕在线播放| 边亲边吃奶的免费视频| 亚洲经典国产精华液单| 欧美高清成人免费视频www| 久久久a久久爽久久v久久| 欧美丝袜亚洲另类| 禁无遮挡网站| 能在线免费看毛片的网站| 九色成人免费人妻av| 精品亚洲乱码少妇综合久久| 午夜精品国产一区二区电影 | 2021天堂中文幕一二区在线观| 久久久精品免费免费高清| 免费看不卡的av| 白带黄色成豆腐渣| 国产免费福利视频在线观看| 2021天堂中文幕一二区在线观| 五月伊人婷婷丁香| 亚洲欧美精品自产自拍| 在线观看av片永久免费下载| 97人妻精品一区二区三区麻豆| 国产精品国产三级专区第一集| 夜夜看夜夜爽夜夜摸| 国产乱人视频| av福利片在线观看| 一级毛片aaaaaa免费看小| 午夜视频国产福利| 亚洲无线观看免费| 美女脱内裤让男人舔精品视频| 免费高清在线观看视频在线观看| 亚洲美女视频黄频| av网站免费在线观看视频| 亚洲av中文字字幕乱码综合| 97在线人人人人妻| kizo精华| 在线天堂最新版资源| 国产午夜精品久久久久久一区二区三区| 日本熟妇午夜| 精品一区二区三卡| 欧美极品一区二区三区四区| 91精品国产九色| 亚洲国产精品国产精品| av女优亚洲男人天堂| 女人十人毛片免费观看3o分钟| 欧美性感艳星| 久久久久久久午夜电影| 久久人人爽人人爽人人片va| 国产成年人精品一区二区| 少妇人妻久久综合中文| 少妇丰满av| av免费观看日本| 1000部很黄的大片| 夜夜看夜夜爽夜夜摸| 人体艺术视频欧美日本| 亚洲国产精品999| 欧美97在线视频| 日韩,欧美,国产一区二区三区| 人体艺术视频欧美日本| 91狼人影院| 可以在线观看毛片的网站| 午夜激情福利司机影院| 秋霞在线观看毛片| 高清毛片免费看| 欧美成人午夜免费资源| 热99国产精品久久久久久7| 69人妻影院| 91在线精品国自产拍蜜月| 欧美成人一区二区免费高清观看| 久久国内精品自在自线图片| 人妻系列 视频| 日韩视频在线欧美| 欧美+日韩+精品| 亚洲国产成人一精品久久久| 狂野欧美白嫩少妇大欣赏| 亚洲美女搞黄在线观看| 99久久中文字幕三级久久日本| 五月伊人婷婷丁香| 男女下面进入的视频免费午夜| 九九久久精品国产亚洲av麻豆| 亚洲最大成人av| 在线 av 中文字幕| 高清在线视频一区二区三区| 国产伦在线观看视频一区| 在线a可以看的网站| 在线看a的网站| 亚洲三级黄色毛片| tube8黄色片| 又黄又爽又刺激的免费视频.| 又黄又爽又刺激的免费视频.| 成人美女网站在线观看视频| 亚洲三级黄色毛片| 日韩一区二区视频免费看| av福利片在线观看| 国产成人freesex在线| 国产精品久久久久久久电影| 欧美区成人在线视频| 熟女电影av网| 大话2 男鬼变身卡| 国产爱豆传媒在线观看| 国产乱人视频| 蜜臀久久99精品久久宅男| 91aial.com中文字幕在线观看| 美女主播在线视频| 久久精品久久久久久噜噜老黄| 久久久成人免费电影| 亚洲精品乱久久久久久| 国产精品一二三区在线看| 人妻系列 视频| 成人亚洲精品av一区二区| 国产日韩欧美在线精品| 内地一区二区视频在线| 内地一区二区视频在线| 男女那种视频在线观看| 色综合色国产| 亚洲精品久久久久久婷婷小说| 国产成人免费观看mmmm| 丝袜美腿在线中文| 亚洲精品456在线播放app| 尤物成人国产欧美一区二区三区| 国产精品蜜桃在线观看| 久久久a久久爽久久v久久| 男女无遮挡免费网站观看| 免费av观看视频| 欧美xxⅹ黑人| 草草在线视频免费看| 免费观看性生交大片5| 建设人人有责人人尽责人人享有的 | 又爽又黄a免费视频| 高清视频免费观看一区二区| 久久精品国产亚洲av涩爱| 永久免费av网站大全| 爱豆传媒免费全集在线观看| 校园人妻丝袜中文字幕| 日韩免费高清中文字幕av| 狂野欧美白嫩少妇大欣赏| 中文字幕人妻熟人妻熟丝袜美| 国产毛片a区久久久久| 免费高清在线观看视频在线观看| 有码 亚洲区| 国产av国产精品国产| 综合色丁香网| 在线观看免费高清a一片| 在线免费十八禁| 真实男女啪啪啪动态图| 亚洲成色77777| 日本色播在线视频| 精品一区二区免费观看| 777米奇影视久久| 嫩草影院精品99| 最近中文字幕高清免费大全6| av女优亚洲男人天堂| 日本一本二区三区精品| 成年女人在线观看亚洲视频 | 3wmmmm亚洲av在线观看| 午夜免费鲁丝| 国产白丝娇喘喷水9色精品| 22中文网久久字幕| 亚洲国产欧美人成| 大香蕉97超碰在线| 中文字幕人妻熟人妻熟丝袜美| 女人久久www免费人成看片| 精品人妻视频免费看| 99久久精品热视频| 亚洲国产欧美人成| 美女内射精品一级片tv| 人妻系列 视频| 欧美潮喷喷水| 美女脱内裤让男人舔精品视频| 国产精品国产三级专区第一集| 少妇丰满av| 日韩大片免费观看网站| 综合色av麻豆| 人妻系列 视频| 日日摸夜夜添夜夜添av毛片| 亚洲激情五月婷婷啪啪| 激情五月婷婷亚洲| 麻豆精品久久久久久蜜桃| 人人妻人人爽人人添夜夜欢视频 | 婷婷色麻豆天堂久久| 国产精品秋霞免费鲁丝片| 91精品伊人久久大香线蕉| 久久久久久久国产电影| 男女那种视频在线观看| kizo精华| 国产精品一区二区性色av| 少妇丰满av| 听说在线观看完整版免费高清| 国产淫片久久久久久久久| 国产在线男女| 亚洲精品亚洲一区二区| 美女xxoo啪啪120秒动态图| 成人二区视频| 看十八女毛片水多多多| 免费播放大片免费观看视频在线观看| 久久久久久国产a免费观看| 一区二区三区四区激情视频| 久久久久精品性色| 午夜精品国产一区二区电影 | 成人免费观看视频高清| 国产精品福利在线免费观看| 国产黄片视频在线免费观看| 国产69精品久久久久777片| 免费黄频网站在线观看国产| av在线亚洲专区| 寂寞人妻少妇视频99o| 久久精品久久久久久噜噜老黄| 欧美国产精品一级二级三级 | 人妻制服诱惑在线中文字幕| 91久久精品电影网| 亚洲国产av新网站| 久久综合国产亚洲精品| 亚洲国产精品成人久久小说| 国产精品麻豆人妻色哟哟久久| 午夜免费观看性视频| 大片免费播放器 马上看| 亚洲精品国产av蜜桃| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍三级| 人妻少妇偷人精品九色| 看黄色毛片网站| 免费观看在线日韩| 亚洲欧美精品专区久久| 国产精品国产三级国产av玫瑰| 日本av手机在线免费观看| 91精品国产九色| 亚洲经典国产精华液单| av又黄又爽大尺度在线免费看| 久久精品夜色国产| 91狼人影院| 久久99热这里只有精品18| 18禁在线无遮挡免费观看视频| 国产色婷婷99| 精品国产一区二区三区久久久樱花 | 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频 | 国产国拍精品亚洲av在线观看| 国产人妻一区二区三区在| 国产伦精品一区二区三区四那| 免费观看av网站的网址| 最近最新中文字幕免费大全7| 丰满少妇做爰视频| 成人一区二区视频在线观看| 2022亚洲国产成人精品| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 我的老师免费观看完整版| 亚洲国产欧美人成| 三级国产精品片| 久久女婷五月综合色啪小说 | 国产免费视频播放在线视频| 老女人水多毛片| 久久久久久久久久人人人人人人| 亚洲真实伦在线观看| 男女国产视频网站| 在线免费观看不下载黄p国产| 夜夜爽夜夜爽视频| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 少妇熟女欧美另类| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 麻豆精品久久久久久蜜桃| 成年版毛片免费区| 国国产精品蜜臀av免费| 国产亚洲一区二区精品| 日韩一区二区视频免费看| 22中文网久久字幕| 18禁在线播放成人免费| 免费看av在线观看网站| 欧美xxⅹ黑人| 久久精品久久精品一区二区三区| 美女被艹到高潮喷水动态| 少妇人妻 视频| 人妻系列 视频| 最近中文字幕2019免费版| 一级二级三级毛片免费看| 涩涩av久久男人的天堂| 寂寞人妻少妇视频99o| 亚洲自偷自拍三级| 精品人妻熟女av久视频| 国产淫语在线视频| 亚洲精品,欧美精品| 男插女下体视频免费在线播放| 成人综合一区亚洲| 日韩欧美精品v在线| 精品久久久精品久久久| 少妇的逼水好多| 日本色播在线视频| 少妇猛男粗大的猛烈进出视频 | 国产男女超爽视频在线观看| 日韩不卡一区二区三区视频在线| 午夜福利视频精品| 久久精品国产a三级三级三级| 夜夜看夜夜爽夜夜摸| 国产精品精品国产色婷婷| 精品99又大又爽又粗少妇毛片| 亚洲欧美一区二区三区黑人 | 久久久国产一区二区| 欧美xxxx黑人xx丫x性爽| 日日啪夜夜爽| 寂寞人妻少妇视频99o| 黄色怎么调成土黄色| 观看免费一级毛片| 免费观看在线日韩| 国产欧美另类精品又又久久亚洲欧美| 午夜福利视频精品| 亚洲国产精品国产精品| 国产午夜福利久久久久久| 国产高清国产精品国产三级 | 国产精品嫩草影院av在线观看| 日本一本二区三区精品| 插逼视频在线观看| 国产黄a三级三级三级人| 久久久久性生活片| 久久久久久久久久久免费av| 麻豆国产97在线/欧美| 99热全是精品| 免费看光身美女| 国产亚洲5aaaaa淫片| 欧美日韩视频高清一区二区三区二| 女的被弄到高潮叫床怎么办| 插逼视频在线观看| 一级毛片久久久久久久久女| 亚洲人成网站在线播| 国产午夜福利久久久久久| 2018国产大陆天天弄谢| 特大巨黑吊av在线直播| 国产午夜福利久久久久久| 亚洲av免费在线观看| 亚洲四区av| 成人欧美大片| 日本与韩国留学比较| 91久久精品国产一区二区成人| 国产欧美日韩精品一区二区| 九九久久精品国产亚洲av麻豆| 97精品久久久久久久久久精品| 波野结衣二区三区在线| 91久久精品国产一区二区三区| 日日啪夜夜爽| 国产精品福利在线免费观看| 国产成人a∨麻豆精品| 久热久热在线精品观看| h日本视频在线播放| 国产国拍精品亚洲av在线观看| 99热6这里只有精品| 精品久久久久久久久亚洲| 国产在视频线精品| 精品少妇黑人巨大在线播放| 大片电影免费在线观看免费| 狠狠精品人妻久久久久久综合| 久久久久久九九精品二区国产| 春色校园在线视频观看| 国产伦在线观看视频一区| 欧美精品国产亚洲| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在| 国国产精品蜜臀av免费| 青春草国产在线视频| 伦理电影大哥的女人| 国产精品人妻久久久影院| 亚洲精品成人av观看孕妇| 九草在线视频观看| 国产精品一及| 国产成人aa在线观看| 不卡视频在线观看欧美| 色吧在线观看| 国产中年淑女户外野战色| 欧美激情久久久久久爽电影| 成人国产av品久久久| 欧美国产精品一级二级三级 | 中文字幕av成人在线电影| 毛片一级片免费看久久久久| 成人综合一区亚洲| 国产成人91sexporn| 青春草国产在线视频| 黄色日韩在线| www.色视频.com| 亚洲av不卡在线观看| 亚洲真实伦在线观看| 51国产日韩欧美| 国产午夜福利久久久久久| 久久97久久精品| 可以在线观看毛片的网站| 精品国产三级普通话版| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站| 亚洲综合精品二区| 一个人看的www免费观看视频| 观看美女的网站| 人妻 亚洲 视频| 在线观看三级黄色| 又爽又黄a免费视频| 精品久久久久久久久亚洲| 国产精品一二三区在线看| 22中文网久久字幕| 一级毛片我不卡| 亚洲,欧美,日韩| 各种免费的搞黄视频| 性色av一级| 全区人妻精品视频| 深爱激情五月婷婷| 真实男女啪啪啪动态图| 亚洲欧美精品专区久久| 男插女下体视频免费在线播放| 色5月婷婷丁香| 18禁裸乳无遮挡动漫免费视频 | 欧美日韩一区二区视频在线观看视频在线 | 精品国产乱码久久久久久小说| 一本一本综合久久| 天天一区二区日本电影三级| 又粗又硬又长又爽又黄的视频| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站| 91在线精品国自产拍蜜月| 成年女人在线观看亚洲视频 | 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| 亚洲国产欧美人成| 国产乱人偷精品视频| 日日啪夜夜撸| 99久久精品国产国产毛片| 精品少妇黑人巨大在线播放| 啦啦啦在线观看免费高清www| 久久久久久久大尺度免费视频| 成人毛片a级毛片在线播放| videossex国产| 99re6热这里在线精品视频| 777米奇影视久久| 内地一区二区视频在线| 少妇人妻一区二区三区视频| 亚洲精品中文字幕在线视频 | 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 欧美zozozo另类| 久久久久久久久久成人| 有码 亚洲区| 中文在线观看免费www的网站| 深夜a级毛片| 欧美xxxx性猛交bbbb| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频| 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 亚洲最大成人中文| 日本av手机在线免费观看| 亚洲,欧美,日韩| 一级毛片久久久久久久久女| 日韩在线高清观看一区二区三区| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 视频区图区小说| 国产爽快片一区二区三区| 日韩精品有码人妻一区| 别揉我奶头 嗯啊视频| 色吧在线观看| 亚洲精品一区蜜桃| 亚洲欧洲国产日韩| 国产在线一区二区三区精| 九九久久精品国产亚洲av麻豆| 亚洲不卡免费看| www.色视频.com| 精品人妻一区二区三区麻豆| 国产成人91sexporn| 涩涩av久久男人的天堂| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 日日撸夜夜添| 男女那种视频在线观看| 久久久国产一区二区| 日韩成人伦理影院| 国产av国产精品国产| 97在线人人人人妻| 国产高潮美女av| 欧美丝袜亚洲另类| 久久精品国产自在天天线| 免费黄频网站在线观看国产| 欧美少妇被猛烈插入视频| 色5月婷婷丁香| 日日摸夜夜添夜夜爱| 在线天堂最新版资源| 下体分泌物呈黄色| 只有这里有精品99| 久久久久久久国产电影| 一级毛片 在线播放| 一级毛片aaaaaa免费看小| 舔av片在线| 九草在线视频观看| 日本一二三区视频观看| 亚洲一区二区三区欧美精品 | 插逼视频在线观看| 亚洲成人久久爱视频| 国产美女午夜福利| 亚洲精品久久久久久婷婷小说| 美女视频免费永久观看网站| 99热网站在线观看| 日韩一区二区视频免费看| 久久人人爽av亚洲精品天堂 | 肉色欧美久久久久久久蜜桃 | 亚洲精品乱码久久久久久按摩| 性色av一级| 男女下面进入的视频免费午夜| 丝袜喷水一区| 一级a做视频免费观看| 国产伦理片在线播放av一区| 日日摸夜夜添夜夜爱| 国产精品久久久久久久电影| 国产美女午夜福利| 黄片wwwwww| 涩涩av久久男人的天堂| 色哟哟·www| 国产成人freesex在线| 精品久久久精品久久久| 午夜福利高清视频| 啦啦啦啦在线视频资源| 日韩制服骚丝袜av| 久热这里只有精品99| 校园人妻丝袜中文字幕| av天堂中文字幕网| 春色校园在线视频观看| 久久精品人妻少妇| 在线亚洲精品国产二区图片欧美 | 男人爽女人下面视频在线观看| 免费看a级黄色片| 免费黄网站久久成人精品| 一本色道久久久久久精品综合| 精品久久久久久久人妻蜜臀av| 国国产精品蜜臀av免费| www.av在线官网国产| 国产亚洲精品久久久com| 国产黄色免费在线视频| 成年av动漫网址| 日本猛色少妇xxxxx猛交久久| 午夜视频国产福利| 综合色丁香网| av天堂中文字幕网| 丝袜脚勾引网站|