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

    The characteristics of negative corona discharge and radio interference at different altitudes based on coaxial wire-cylinder gap

    2024-03-19 02:36:58WanglingHE何旺齡HongyuWEl魏宏宇YinluZHANG張胤祿YongcongLlU劉永聰YunpengLlU劉云鵬BaoquanWAN萬保權(quán)YangZHOU周陽YanzhaoWANG王延召andZheyuanGAN干喆淵
    Plasma Science and Technology 2024年2期

    Wangling HE (何旺齡),Hongyu WEl (魏宏宇),Yinlu ZHANG (張胤祿),Yongcong LlU (劉永聰),Yunpeng LlU (劉云鵬),Baoquan WAN (萬保權(quán)),Yang ZHOU (周陽),Yanzhao WANG (王延召) and Zheyuan GAN (干喆淵)

    1 Hebei Key Laboratory of Green and Efficient New Electrical Materials and Equipment,North China Electric Power University,Baoding 071000,People’s Republic of China

    2 State Grid Suzhou Power Supply Company,Suzhou 215000,People’s Republic of China

    3 State Key Laboratory of Power Grid Environmental Protection (China Electric Power Research Institute),Wuhan 430000,People’s Republic of China

    Abstract The corona discharge from transmission lines in high-altitude areas is more severe than at lower altitudes.The radio interference caused thereby is a key factor to be considered when designing transmission lines.To study the influence of altitude on negative corona characteristics,an experimental platform comprising a movable small corona cage was established: experiments were conducted at four altitudes in the range of 1120-4320 m,and data on the corona current pulse and radio interference level of 0.8-mm diameter fine copper wire under different negative voltages were collected.The experimental results show that the average amplitude,repetition frequency and average current of the corona current pulse increase with increasing altitude.The dispersion of pulse amplitude increases with increase in altitude,while the randomness of the pulse interval decreases continuously.Taking the average current as an intermediate variable,the relationship between radio interference level and altitude is obtained.The result of this research has some significance for understanding the corona discharge characteristics of ultra-highvoltage lines.

    Keywords: altitude,negative corona,corona current pulse,radio interference level

    1.lntroduction

    China’s terrain is high in the west and low in the east,and the energy in the west is inversely distributed with the load center in the east.Areas with an altitude of over 2000 m in China account for 37% of the land area,and the average elevation of the Qinghai-Tibet Plateau is around 4000 m.It is an inevitable choice for China to build an ultra-high-voltage direct current (UHV DC) transmission network to meet its energy needs,and an UHV DC system has the advantages of long distance,large capacity and low loss [1,2].At the same time,corona discharge on the surface of transmission lines will cause some electromagnetic environmental problems,such as audible noise,radio interference (RI) and corona loss.RI can affect the normal operation of radio facilities around the transmission line,so its impact needs to be taken into account when erecting transmission lines.

    With increase in altitude,the density of air decreases,the electron mean free path increases,effective collisional ionization is more likely to occur,the corona discharge phenomenon is more severe and the RI of transmission lines in high-altitude areas is consequently more serious [3].A negative pole line is important component for bipolar operation of a DC transmission system,so it is necessary to discuss the negative discharge and its level of RI.

    Researchers have conducted a lot of work on negative corona discharges.In 1938,Trichel first studied the regular repetitive current pulses that occur in a negative corona discharge and proposed the reasons for the generation of these repetitive current pulses [4].In 1974,Lama and Gallo conducted a comprehensive analysis of the characteristics of negative polarity corona pulses for the first time [5],while in 1981,Henson measured and theoretically derived a mathematical model for the needle plate gap in liquid helium [6].In 1985,Morrow proposed a one-dimensional simulation model for negative corona current pulses [7].In 2002,Akishevet alobtained a series of Trichel pulses based on a two-dimensional model of a negative corona discharge in air[8].In 2004,Yamada measured the negative corona current using a ‘point to grid’ electrode system and obtained an empirical formula for the corona volt-ampere characteristics[9].In 2007,Hui studied the influence of air pressure on the characteristics of a negative corona pulse based on a needleplate electrode in an artificial climate chamber.Hui believed that the decrease in air pressure and the increase in electron density were the main reasons for the increase in current pulse amplitude [10].In 2014,Yinet alsimulated Trichel pulses in a negative corona discharge based on a timedomain finite-volume approach [11].In 2014,Denget alobtained the effect of humidity on negative corona Trichel pulses [12],while in 2021,Liuet alconducted experiments using small corona cages at different altitudes to study the characteristics of AC negative half-cycle corona current pulses at different altitudes [13].

    As for the RI on transmission lines,Adams first proposed the concept of an RI generation function in 1956 [14].Moreau and Gary studied the effect of weather conditions and conductor surface condition on RI on transmission lines [15],while Heet almeasured the RI excitation function values of different bundle conductors at different altitudes,proposing a correction term for the RI level which fitted well for regions with different altitudes below 4300 m[16].Based on a large corona cage,Peiet almeasured the RI current of various conductor bundles under conditions of heavy rain.The proposed excitation function is suitable for the prediction of RI level of UHV alternating current (AC)transmission lines in plains [17].

    Most of these studies used numerical simulation methods or model experiments based on artificial climate chambers,and few experiments have been conducted on discharge characteristics under actual environmental scenarios.Therefore,it is necessary to conduct corona discharge experiments at different altitudes,further study the influence of altitude on negative corona discharge and obtain negative corona discharge current characteristics and their connection to RI.

    A movable small corona cage and a high-frequency current acquisition system were used to test negative corona characteristics at four different altitudes,and the RI levels were measured by an EMI receiver.The influence of altitude on the negative corona characteristics was analyzed,and the relationship between the average corona current and the RI level constructed.The results can provide reference values for understanding the negative discharge and its relationship with RI in high-altitude areas.

    2.Experimental principle

    2.1.Experimental platform

    A mobile small corona cage and high-frequency current measurement system were used to conduct corona discharge experiments at different altitudes (figure 1).The corona cage was fabricated in metal,with the cage wall being a singlelayer network structure with a circular cross section of 0.4 m diameter and length of 1 m.The support frame and wire support rod of the corona cage were made of epoxy resin materials to ensure reliable insulation and prevent leakage current from affecting the measurement.Fine copper wire with a diameter of 0.8 mm was selected as the experimental wire.To obtain stable corona discharge,the surface of the wire was polished with abrasive paper.

    Due to the end effect,the electric field on the surface of the conductor at the end of the corona cage will be distorted.The electric field intensity at the distortion is large,and the corona discharge is severe,which affect the accuracy of the whole corona current measurement system.Therefore,a grading ring was inserted between the conductor support rod and the small corona cage to eliminate the end effect.

    Figure 1.Layout of experimental setup.

    2.2.Measurement system and experimental method

    To ensure the safety of the experimental personnel,the highand low-potential parts of the measurement system must be reliably isolated.Therefore,a high-frequency current acquisition system was used in the present research to ensure electrical isolation between high and low potential through the optical fiber transmission measurement signal.The experimental circuit diagram is shown in figure 2.The inductive trap can block the high-frequency current generated by the corona discharge of the wire into the power supply end,and can also prevent interference on the power supply side; the coupling capacitor can couple the high-frequency component of the conductor corona current to the low-voltage side;the impedance at the end of the coupling loop prevents signal reflection between the measurement signal and the receiver,and its value was set to 50 Ω.

    To ensure the effectiveness of corona current pulse measurement,the current acquisition system is required to distinguish the pulse with a minimum rising edge of less than 10 ns.The acquisition card is a U2702 A-USB modular oscilloscope (Agilent,USA) with a 12 V battery to supply power separately.The acquisition card has a bandwidth of 200 MHz,a single-channel sampling frequency of 5 × 108samples s-1,a storage capacity of 16 million samples and a minimum sampling interval of 2 ns.At the same time,to eliminate the influence of resistance,lead inductance and stray capacitance on the ground of the current sensor,the current sensor used in this article adopted coaxial parallel connection of multiple non-inductive resistors and the total resistance value of the current sensor was 400 Ω.

    Corona discharge is caused by the electric field strength on the surface of the wire being greater than the corona field strength,so the electric field strength on the surface of the wire was adjusted by changing the voltage applied to the wire.The corona current pulse and RI level data were measured at different altitudes and different applied voltages.

    After the onset of conduction,data were measured for every 5 kV step,and the following negative DC voltages were applied: -15 kV,-20 kV,-25 kV and -30 kV.After the corona discharge state stabilized at each applied voltage,10 current pulse measurements were repeated and the waveforms recorded.

    2.3.Experimental site and environment

    Figure 2.Experimental circuit diagram.

    To study the performance of the corona at different altitudes,experiments were conducted at four different altitudes in Tianshui,Gansu Province (altitude 1120 m),Baiyin,Gansu Province (altitude 1680 m),Lhasa,Tibet Autonomous Region (altitude 3660 m) and Yangbajing,Tibet Autonomous Region (altitude 4320 m).

    All experiments were carried out at an indoor site with sufficient space.The relative humidity of the measurement environment was maintained at 40%-60%,and the ambient temperature was maintained between 20 °C and 30 °C.Studies have shown that temperature and humidity changes in this small range have little effect on the characteristics of negative corona discharge [10,18],therefore,in this work,the influences of temperature and humidity on the experimental results are ignored and the effect of altitude is mainly considered.

    3.Measurement and analysis of corona current pulses at different altitudes

    3.1.Typical corona current pulse sequence

    The typical corona current pulse sequences at different altitudes when a voltage of -20 kV was applied to the copper wire is shown in figure 3.

    Taking the pulse waveform at an altitude of 1120 m as an example,the characteristics of the corona current pulse were investigated (figure 4).The corona current pulse was fitted by the least squares method,as shown by the imaginary line.Here,Ipis the pulse peak andTis the duration of a single pulse.

    3.2.Characteristics of the corona current pulse

    3.2.1.Average value of the pulse amplitude and pulse repetition frequencyThe corona current pulse data measured by the experiment were analyzed,and the average value of the corona current pulse amplitude and the current pulse repetition frequency were obtained and discussed (figures 5 and 6).The negative corona pulse amplitude and the average current were taken as absolute values.

    Figure 3.Corona current pulse sequences at different altitudes.

    Figure 4.Pulse waveform parameter.

    It can be seen from figures 5 and 6 that the average value of the pulse amplitude and repetition frequency increase with increasing altitude.For a voltage of -20 kV applied to the conductor,when the altitude increases from 1120 to 4320 m the average value of the pulse amplitude increases from 0.95 to 1.5 mA and the pulse repetition frequency increases from 18.2 to 40.2 kHz.The change of pulse repetition frequency is greater than the change in average value of the pulse amplitude.By analyzing the development of the corona discharge pulse,it can be seen that with increase in altitude the density of electrons and positive ions near the electrode increases when the pulse occurs,resulting in more severe spatial distortion of the electric field,a higher electric field intensity on the surface of the conductor and a greater corona pulse amplitude.At the same time,the increase in altitude makes the positive and negative ions migrate faster in the dead time after the end of the pulse,so the electric field intensity on the surface of the wire is better able to return to the state at the beginning of the corona discharge and the corona pulse repetition frequency increases.

    3.2.2.Average value of corona currentTraditional research on corona current characteristics is mostly based on the average current measured in experiments,and the volt-ampere characteristics of the corona discharge have been studied [19].In the present work,the measured current pulse data were processed to allow calculation of the average value of the corona current,as shown in equation (1)

    whereIaveis the average value of corona current (inμA),Tis the total duration of the measured current pulse fragment (in s) andiis the measured instantaneous value of corona current (inμA).

    The change in the average value of the corona current with altitude is plotted in figure 7.

    Figure 5.Average value of the pulse amplitude at different altitudes.

    Figure 6.Pulse repetition frequency at different altitudes.

    Figure 7 shows that the average value of corona current increases with increasing altitude.When the altitude is below 2000 m,the average current increases slowly when the applied voltage is low.At altitudes above 2000 m,the average current increases rapidly,and the greater the altitude,the greater the increase.When the voltage applied to the conductor is -20 kV and the altitude increases from 1120 to 1680 m,the average current increases from 1.9 to 2.49μA.When the altitude increases from 3660 to 4320 m,the average current increases from 5.71 to 8.58μA.With increasing altitude,the amplitude of the corona current pulse increases,the repetition frequency increases and,ultimately,the average value of the current increases.

    3.3.Statistical analysis of corona current pulses

    To estimate the random distribution of negative DC corona current pulse eigenvalues,it is necessary to evaluate the probability density distribution of pulse amplitude and pulse interval.

    3.3.1.Probability density distribution of pulse amplitudeThe probability distribution of the amplitude of the negative corona current pulse at different altitudes when the voltage applied to the conductor is -20 kV is shown in figure 8.At low altitude (1120 m and 1680 m),the pulse amplitude is mainly concentrated at about 1 mA,which conforms to the characteristics of a normal distribution.As the altitude increases,the pulse amplitude gradually shifts to the right.When the altitude is higher than 3660 m,a larger part of the pulse amplitude appears,and the second peak appears at about 2 mA; this is characterized by a bimodal distribution,indicating that the discharge pulse may have a new form with a higher amplitude.

    The standard deviation of corona current pulse amplitude at different altitudes for an applied voltage of -20 kV was calculated and is given in table 1.

    For a constant voltage applied to the wire,as the altitude increases the randomness of the corona current pulse amplitude distribution increases,and the standard deviation of the pulse amplitude increases.This is because as the altitude increases,pulses with higher amplitudes will appear while pulses with lower amplitudes still exist,resulting in an increase in the dispersion of pulse amplitudes.

    3.3.2.Probability density distribution of pulse intervalThe probability distribution of the negative pulse interval at different altitudes for an applied voltage of -20 kV is shown in figure 9.

    Figure 9 shows that the probability distribution of the corona current pulse interval shows a trend of high on the left and low on the right,with exponential attenuation consistent with the probability density function curve of the exponential distribution.At an altitude of 1120 m,the maximum time interval between corona current pulses is close to 300μs,and their distribution is random.With increasing altitude,the maximum time interval between pulse gradually decreases,the overall distribution shifts to the left and the randomness of the distribution decreases.Above an altitude of 1680 m,the maximum value of the pulse time interval is less than 200μs and the probability density of small time intervals continues to increase.

    The standard deviation of the pulse interval is shown intable 2: with increasing altitude,the dispersion of the pulse time interval distribution decreases and the standard deviation decreases.This is because when the voltage applied to the wire is constant,as the altitude increases the number of corona current pulses per unit time increases,the corona discharge tends to be more stable and regularity is enhanced,so the time interval distribution is more concentrated.

    Figure 8.Histogram of the probability distribution of pulse amplitude at different altitudes: (a) 1120 m,(b) 1680 m,(c) 3660 m,(d)4320 m.

    3.4.Variation of corona current with altitude

    To quantify the variation of corona current characteristics with altitude,the influence of altitude on the negative DC corona characteristics of conductors is further analyzed.In most research,the empirical formula describing the corona current characteristics at different altitudes is derived using a data-fitting method.The average value of corona current contains information such as the voltage applied to the wire,the amplitude of the corona current pulse and the repetition frequency.Therefore,use of the average value of the current as a characteristic can more comprehensively and accurately describe the influence of altitude on the characteristics of the corona current pulse.

    Under the same applied voltage,change in the electrode radius and gap distance will cause the electric field intensity to change,and the corona current characteristics will also change.Therefore,no accurate fitting formula can be obtained by taking the applied voltage as the independent variable,and the electric field strength on the surface of the conductor is directly related to the corona discharge characteristics.In addition,because the present research is based on scale-model experiments,the voltage of the actual line is quite different but the surface field strengths of the two are similar.Taking the electric field strength as the independent variable,the equivalence of the scale model to the actual line can be realized and simulation can be done.

    Based on the principle of electrostatic field,according to equation (2) the applied voltage of the wire is converted into the field strength on the surface of the wire.Taking the surface field strength of the wire as the independent variable,the change in average corona current with electric field strength at different altitudes is described (see figure 10)

    whereUis the conductor voltage (in kV),Eis the electric field intensity on the conductor surface (in kV cm-1),ris the radius of the wire (in cm) andRis the outer electrode radius(in cm).

    It can be seen from figure 10 that at different altitudesthe average corona current increases with increase in the surface field strength on the conductor,and the average corona current increases with applied voltage as a power function.Through the pre-treatment and preliminary fitting of the experimental data,it can be found that the relationship between the average corona current and the surface field strength of the conductor takes the formK(E-E0)n,withnbeing about 1.6,which is consistent with previous results [20].

    Table 2.Standard deviation of the pulse interval.

    Figure 9.Histogram of pulse interval time distribution at different altitudes: (a) 1120 m,(b) 1680 m,(c) 3660 m,(d) 4320 m.

    Figure 10.Relationship between the average corona current and the field strength on the wire surface.

    The relationship between the average corona current and the electric field intensity of the conductor surface is fitted by equation (3)

    whereIaveis the average value of corona current (inμA),Krepresents a linear coefficient of proportionality,Eis the electric field intensity on the surface of the conductor (in kV cm-1) andE0is the corona onset field strength (in kV cm-1).When the surface electric field intensity of the wire is close toE0,the average corona current tends to 0 and the corona discharge tends to stop.

    The results for fitting the altitude andKare shown in equation (4) and figure 11

    wherehis the altitude (in m).

    The relationship between the average corona current and the altitude is obtained by equation (5)

    whereIaveis the average value of corona current (inμA) andEis the surface field strength of the conductor (in kV cm-1).

    4.Relationship between corona pulse and Rl level

    Figure 11.The coefficient K in relation to altitude.

    The variation of the average corona current with altitude is given above: this covers the information on the applied voltage of the wire,the amplitude of the corona current pulse,the repetition frequency and so on.The average value of the corona current can be used to quantify the relationship between the corona current and RI level.Therefore,the average value of the corona current can be used as an intermediate variable to assess the influence of altitude on the negative corona effect.

    Regression analysis was used to establish the correlation between the average corona current and the RI level [21].Based on the calculated results for average corona current,the relationship between the negative corona current and RI level at different altitudes was analyzed with an RI level of 0.5 MHz taken as the characteristic quantity.The RI level after eliminating background interference is summarized in table 3.

    The level of RI and corona current pulse measurement data at different altitudes were analyzed.As shown in figure 12,at different altitudes the RI level increases with increase in the corona current and generally shows a logarithmic relationship.Altitude itself has no obvious relationship with RI level.

    Based on the non-linear least squares method,the relationship between the RI level and the average corona current is fitted to quantify the relationship between them (figure 12)and determine the characteristics of the power function curve(table 4,figure 13).

    In table 4,aandbare the optimal parameter values of the fitting curve,R2is the coefficient of determination,with a value range of [0,1].The closerR2is to 1,the more accurate the fitted formula.RMSE is the fitting standard deviation of the regression system.The closer the RMSE is to 0,the better the model fits the data.

    Table 3.0.5 MHz RI level at different altitudes (unit dB).

    Figure 12.The relationship between RI and the average value of corona current.

    Table 4.Fitting results for RI level and average corona current.

    Figure 13.Fitted curve of RI versus average corona current.

    In summary,the functional relationship between corona current pulse and RI level can be obtained as shown in equation (6).That is,at the 0.5 MHz RI level,RI0.5and the average corona currentIare roughly in a power function relationship

    whereRI0.5is the 0.5 MHz RI level (in dB) andIis the average corona current (inμA).

    By substituting equation (5) into equation (6),the relationship between RI level and altitude can be obtained (equation (7))

    whereEis the field strength of the conductor surface (in kV cm-1),Kis a function of altitudehand the relationship between them is given by equation (4).

    5.Conclusions

    The negative corona discharge of a conductor was tested at four different altitudes using a movable small corona cage,and the influence of altitude on negative corona current pulse and RI was discussed.The following conclusions were obtained:

    (1) The average amplitude,repetition frequency and average current of the corona current pulse increase with increasing altitude.When the altitude increases from 1120 to 4320 m,the average amplitude of the corona current pulses increases from 0.95 to 1.5 mA,the repetition frequency increases from 18.2 to 40.2 kHz and the average current increases from 1.9 to 8.58μA.

    (2) The amplitude of the pulse follows a lognormal distribution.As the altitude increases,the dispersion of the amplitude distribution increases.At an altitude of 3660 m,some high-amplitude pulses appear,resulting in two peaks appearing in the probability distribution.The pulse time interval follows an exponential distribution.As the altitude increases,the standard deviation decreases and the probability density distribution becomes more concentrated.

    (3) The average corona currentIaveand the electric field strengthEof the conductor surface were selected as the characteristic quantities for fitting the relationship between the two at different altitudes.The results show thatIaveis proportional to (E-E0)1.6.The constant of proportionalityKlinking the two is a quadratic function of altitudeh.

    (4) Through regression analysis,the relationship between RI and the average value of corona current shows power function growth (0 <α< 1).Combined with the relationship between the average value of corona current and the change in wire surface field strength and altitude,the relationship of RI level to altitude is obtained.

    Acknowledgments

    This work was supported by the Science and Technology Project of State Grid Corporation of China (No.5200-202155587A-0-5-GC).

    成人国产av品久久久| 欧美日韩精品网址| 精品久久久精品久久久| 亚洲国产精品成人久久小说| 老汉色∧v一级毛片| 大片免费播放器 马上看| www.精华液| 亚洲国产av影院在线观看| 免费不卡黄色视频| 天堂8中文在线网| 黄片播放在线免费| 国产亚洲精品一区二区www | 老司机午夜十八禁免费视频| 欧美日韩亚洲国产一区二区在线观看 | 视频区欧美日本亚洲| 肉色欧美久久久久久久蜜桃| 男人操女人黄网站| 高清av免费在线| 最新的欧美精品一区二区| 老司机福利观看| 啦啦啦 在线观看视频| 欧美精品高潮呻吟av久久| 欧美激情久久久久久爽电影 | 美女高潮喷水抽搐中文字幕| 咕卡用的链子| 少妇人妻久久综合中文| 久久精品亚洲av国产电影网| 国产精品自产拍在线观看55亚洲 | 成年动漫av网址| 91麻豆精品激情在线观看国产 | 丁香六月欧美| 国产成+人综合+亚洲专区| 欧美精品亚洲一区二区| 国产高清视频在线播放一区 | 成年动漫av网址| 国产在线视频一区二区| 永久免费av网站大全| 在线观看免费日韩欧美大片| 777米奇影视久久| 日韩免费高清中文字幕av| av线在线观看网站| 1024香蕉在线观看| 国内毛片毛片毛片毛片毛片| 国产精品自产拍在线观看55亚洲 | 午夜福利视频精品| 美女国产高潮福利片在线看| 韩国精品一区二区三区| 亚洲成人国产一区在线观看| 99久久人妻综合| 一二三四社区在线视频社区8| 九色亚洲精品在线播放| 国产成+人综合+亚洲专区| 免费av中文字幕在线| 制服诱惑二区| 日韩视频一区二区在线观看| 国产av精品麻豆| 亚洲欧美激情在线| 9色porny在线观看| 国产亚洲欧美精品永久| 日本vs欧美在线观看视频| 大型av网站在线播放| 男女免费视频国产| 69精品国产乱码久久久| 一本色道久久久久久精品综合| 欧美日韩视频精品一区| 在线天堂中文资源库| 久久人人97超碰香蕉20202| 天天躁日日躁夜夜躁夜夜| 我要看黄色一级片免费的| 欧美乱码精品一区二区三区| 丁香六月欧美| 久久国产亚洲av麻豆专区| 亚洲欧洲日产国产| 99久久99久久久精品蜜桃| 久久久久国产精品人妻一区二区| 乱人伦中国视频| 在线永久观看黄色视频| 一二三四社区在线视频社区8| 久9热在线精品视频| 人人妻人人添人人爽欧美一区卜| 国产精品熟女久久久久浪| 亚洲精品中文字幕在线视频| 自线自在国产av| 在线观看免费视频网站a站| 久久久久精品人妻al黑| 免费高清在线观看视频在线观看| 一级毛片电影观看| 一本久久精品| 自线自在国产av| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区在线观看99| 高清黄色对白视频在线免费看| 国产成人一区二区三区免费视频网站| 欧美日韩国产mv在线观看视频| 亚洲少妇的诱惑av| 99精国产麻豆久久婷婷| 人妻久久中文字幕网| 水蜜桃什么品种好| 在线永久观看黄色视频| 亚洲人成电影免费在线| 无限看片的www在线观看| 亚洲国产毛片av蜜桃av| 男女边摸边吃奶| 久久ye,这里只有精品| 亚洲熟女毛片儿| 久久av网站| 欧美日韩视频精品一区| 热re99久久国产66热| videos熟女内射| 久久女婷五月综合色啪小说| 丝袜人妻中文字幕| 日本wwww免费看| 日本vs欧美在线观看视频| 精品一区在线观看国产| 欧美久久黑人一区二区| 真人做人爱边吃奶动态| 成人影院久久| 一级毛片电影观看| 亚洲情色 制服丝袜| 色播在线永久视频| 水蜜桃什么品种好| 成人国产一区最新在线观看| 久热这里只有精品99| 岛国毛片在线播放| 亚洲欧美成人综合另类久久久| 久久精品国产综合久久久| 91成人精品电影| 亚洲国产精品一区二区三区在线| 两个人看的免费小视频| 一个人免费看片子| 又大又爽又粗| 9色porny在线观看| 一区二区三区四区激情视频| 国产精品 国内视频| 伊人亚洲综合成人网| 99久久国产精品久久久| bbb黄色大片| 国产精品免费大片| av国产精品久久久久影院| 丝袜脚勾引网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品美女久久久久99蜜臀| 人人妻人人澡人人看| 人妻人人澡人人爽人人| 国产一卡二卡三卡精品| 国产麻豆69| 超碰成人久久| 欧美午夜高清在线| 老司机午夜十八禁免费视频| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区综合在线观看| 亚洲国产av新网站| 国产精品1区2区在线观看. | 青春草亚洲视频在线观看| 亚洲精品中文字幕在线视频| 天天添夜夜摸| 欧美另类亚洲清纯唯美| 成年av动漫网址| 精品免费久久久久久久清纯 | 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品影院| 国产成人精品在线电影| 男女国产视频网站| 老熟妇仑乱视频hdxx| 大陆偷拍与自拍| 人妻久久中文字幕网| 激情视频va一区二区三区| 亚洲欧美日韩另类电影网站| 精品久久久久久久毛片微露脸 | 亚洲av成人不卡在线观看播放网 | 一边摸一边做爽爽视频免费| 最近中文字幕2019免费版| 后天国语完整版免费观看| 欧美国产精品一级二级三级| 69av精品久久久久久 | 欧美亚洲日本最大视频资源| 丁香六月天网| 国产精品久久久av美女十八| 色94色欧美一区二区| 桃红色精品国产亚洲av| 欧美中文综合在线视频| 国产精品久久久久久精品电影小说| 另类亚洲欧美激情| 99久久精品国产亚洲精品| 老司机影院成人| kizo精华| 日韩大片免费观看网站| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品电影小说| 精品亚洲乱码少妇综合久久| 亚洲精品美女久久久久99蜜臀| 99国产综合亚洲精品| 中文字幕人妻熟女乱码| 看免费av毛片| 亚洲少妇的诱惑av| 男人操女人黄网站| 欧美性长视频在线观看| 日韩,欧美,国产一区二区三区| 在线av久久热| 12—13女人毛片做爰片一| 美女脱内裤让男人舔精品视频| 国产成人啪精品午夜网站| 飞空精品影院首页| 一级片免费观看大全| 一级毛片精品| 久久天堂一区二区三区四区| 精品国产一区二区久久| 一二三四在线观看免费中文在| 大香蕉久久网| 国产精品欧美亚洲77777| 12—13女人毛片做爰片一| 男女免费视频国产| 99久久综合免费| 亚洲 国产 在线| 精品一区二区三卡| 在线观看舔阴道视频| 99国产精品99久久久久| 欧美成人午夜精品| 一区二区日韩欧美中文字幕| 亚洲av日韩精品久久久久久密| 国产成人精品无人区| 中文字幕人妻熟女乱码| 久久ye,这里只有精品| 国产精品久久久久久精品电影小说| 啦啦啦免费观看视频1| 丝袜美腿诱惑在线| 久久精品久久久久久噜噜老黄| 精品国产乱子伦一区二区三区 | av欧美777| 99国产综合亚洲精品| 最黄视频免费看| 黄片播放在线免费| 99热全是精品| 搡老岳熟女国产| 亚洲性夜色夜夜综合| tocl精华| 国产1区2区3区精品| 自线自在国产av| 国产av一区二区精品久久| 精品人妻熟女毛片av久久网站| 国产av精品麻豆| 香蕉国产在线看| 国产片内射在线| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品古装| 国产精品久久久人人做人人爽| 美女大奶头黄色视频| 首页视频小说图片口味搜索| 久久国产亚洲av麻豆专区| 国产日韩欧美视频二区| 国产日韩欧美视频二区| 女人高潮潮喷娇喘18禁视频| 欧美老熟妇乱子伦牲交| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 999久久久精品免费观看国产| 国产精品香港三级国产av潘金莲| videosex国产| av福利片在线| 窝窝影院91人妻| 侵犯人妻中文字幕一二三四区| 日韩视频在线欧美| 99精国产麻豆久久婷婷| 日韩中文字幕视频在线看片| 超碰成人久久| 亚洲精品久久成人aⅴ小说| 欧美xxⅹ黑人| 一区在线观看完整版| www.熟女人妻精品国产| 国产精品久久久久久精品电影小说| 久久久精品区二区三区| tube8黄色片| 日韩三级视频一区二区三区| 又大又爽又粗| 久久影院123| 看免费av毛片| 国产精品1区2区在线观看. | 午夜福利免费观看在线| 丝袜人妻中文字幕| av一本久久久久| 不卡av一区二区三区| 国产精品麻豆人妻色哟哟久久| 香蕉丝袜av| 午夜激情久久久久久久| 国产日韩欧美视频二区| 欧美激情极品国产一区二区三区| 精品人妻熟女毛片av久久网站| 久久人妻福利社区极品人妻图片| 各种免费的搞黄视频| 国产av国产精品国产| 自线自在国产av| 91成年电影在线观看| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www | 制服诱惑二区| 看免费av毛片| 男女午夜视频在线观看| 人妻 亚洲 视频| 丝袜人妻中文字幕| 久久热在线av| 国产精品一区二区精品视频观看| 久久久精品免费免费高清| 中国美女看黄片| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 亚洲精品国产精品久久久不卡| 午夜老司机福利片| 夜夜骑夜夜射夜夜干| tube8黄色片| 亚洲人成电影免费在线| 国产高清视频在线播放一区 | 免费在线观看影片大全网站| 亚洲色图 男人天堂 中文字幕| 97精品久久久久久久久久精品| 中文字幕另类日韩欧美亚洲嫩草| 日韩免费高清中文字幕av| 日韩熟女老妇一区二区性免费视频| 我的亚洲天堂| 一本一本久久a久久精品综合妖精| 亚洲精品乱久久久久久| 久热爱精品视频在线9| 午夜影院在线不卡| 亚洲九九香蕉| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产av成人精品| 亚洲欧美色中文字幕在线| 丝袜人妻中文字幕| 午夜两性在线视频| a级毛片在线看网站| 色老头精品视频在线观看| 女性生殖器流出的白浆| av一本久久久久| 国产亚洲av片在线观看秒播厂| 日韩制服骚丝袜av| av天堂在线播放| 欧美日韩福利视频一区二区| 午夜视频精品福利| 免费人妻精品一区二区三区视频| 午夜两性在线视频| 中文欧美无线码| 黄网站色视频无遮挡免费观看| 老鸭窝网址在线观看| 老司机影院成人| 成年美女黄网站色视频大全免费| 我要看黄色一级片免费的| 丝袜美腿诱惑在线| av不卡在线播放| 99国产综合亚洲精品| 午夜91福利影院| 一个人免费在线观看的高清视频 | 另类亚洲欧美激情| 高清在线国产一区| 国产不卡av网站在线观看| 国产成人啪精品午夜网站| 男女无遮挡免费网站观看| 久久久久视频综合| 久久精品熟女亚洲av麻豆精品| 婷婷丁香在线五月| 精品一区二区三卡| 免费不卡黄色视频| 手机成人av网站| 成年av动漫网址| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 午夜福利视频精品| 亚洲欧美色中文字幕在线| 一级片免费观看大全| 黄色怎么调成土黄色| www.自偷自拍.com| 国产有黄有色有爽视频| 黑人猛操日本美女一级片| 自线自在国产av| 男女边摸边吃奶| 亚洲成人手机| 男女无遮挡免费网站观看| 波多野结衣一区麻豆| 国产亚洲精品一区二区www | 日本精品一区二区三区蜜桃| 久久久久久久久免费视频了| 日本精品一区二区三区蜜桃| 久久 成人 亚洲| 久久久国产精品麻豆| 亚洲av片天天在线观看| 免费女性裸体啪啪无遮挡网站| 99久久99久久久精品蜜桃| 精品人妻熟女毛片av久久网站| 国产欧美日韩综合在线一区二区| 最近最新免费中文字幕在线| 国产成人啪精品午夜网站| 成人av一区二区三区在线看 | 国产97色在线日韩免费| 欧美日韩黄片免| 日韩熟女老妇一区二区性免费视频| 欧美国产精品一级二级三级| 人人妻,人人澡人人爽秒播| 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 国产精品久久久久久精品古装| 国产亚洲精品久久久久5区| 两个人看的免费小视频| 在线亚洲精品国产二区图片欧美| 成人手机av| 国产精品国产三级国产专区5o| 黄频高清免费视频| 中国国产av一级| 欧美97在线视频| 丝袜人妻中文字幕| av免费在线观看网站| 日韩视频一区二区在线观看| 高潮久久久久久久久久久不卡| 国产精品麻豆人妻色哟哟久久| 国产一区二区三区av在线| 99热全是精品| 亚洲欧美一区二区三区黑人| av国产精品久久久久影院| 午夜福利乱码中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 12—13女人毛片做爰片一| 久久久久国内视频| 五月天丁香电影| 在线天堂中文资源库| 一级片'在线观看视频| 丝袜在线中文字幕| 欧美精品一区二区大全| 久久久久久久精品精品| 搡老岳熟女国产| 国产精品成人在线| 亚洲欧美精品综合一区二区三区| 国产精品九九99| 国产成人啪精品午夜网站| 天天躁夜夜躁狠狠躁躁| 亚洲一码二码三码区别大吗| 午夜91福利影院| 久久久久精品人妻al黑| 制服人妻中文乱码| 免费少妇av软件| 少妇的丰满在线观看| 丝袜喷水一区| 国产欧美日韩综合在线一区二区| 建设人人有责人人尽责人人享有的| 日本五十路高清| 精品人妻熟女毛片av久久网站| 精品一区在线观看国产| 欧美在线黄色| 十分钟在线观看高清视频www| 亚洲 国产 在线| 国产一区二区三区av在线| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美一区二区三区黑人| 一本大道久久a久久精品| 老司机影院毛片| 欧美乱码精品一区二区三区| 亚洲中文av在线| 亚洲国产精品一区三区| 久久av网站| 亚洲国产精品一区三区| 国产不卡av网站在线观看| 老司机福利观看| 国产又爽黄色视频| 午夜老司机福利片| 亚洲 国产 在线| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区mp4| 99国产综合亚洲精品| 我要看黄色一级片免费的| 日韩电影二区| 十八禁人妻一区二区| 我要看黄色一级片免费的| 中文字幕人妻丝袜一区二区| 久久久久国产一级毛片高清牌| 亚洲精品国产av蜜桃| 久久免费观看电影| 一本大道久久a久久精品| 捣出白浆h1v1| 国产高清国产精品国产三级| av在线老鸭窝| 久久性视频一级片| 亚洲专区字幕在线| 天堂中文最新版在线下载| 美女福利国产在线| 欧美午夜高清在线| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 欧美精品一区二区大全| 午夜福利影视在线免费观看| 日韩视频一区二区在线观看| 美女福利国产在线| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产一区二区三区四区第35| 中国国产av一级| 91成年电影在线观看| 男人爽女人下面视频在线观看| 国产色视频综合| 国产日韩欧美亚洲二区| 成人影院久久| 黄色视频,在线免费观看| 九色亚洲精品在线播放| 国产欧美日韩综合在线一区二区| 精品国产一区二区三区四区第35| 日本av手机在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文字幕日韩| 菩萨蛮人人尽说江南好唐韦庄| 午夜免费观看性视频| 一级毛片精品| 精品福利永久在线观看| 三上悠亚av全集在线观看| av网站在线播放免费| 99精品欧美一区二区三区四区| 王馨瑶露胸无遮挡在线观看| 免费看十八禁软件| av天堂在线播放| 精品一区在线观看国产| 免费高清在线观看日韩| 侵犯人妻中文字幕一二三四区| 视频区图区小说| 日韩欧美一区视频在线观看| 午夜福利视频精品| 美女大奶头黄色视频| 亚洲精品国产av蜜桃| 美女扒开内裤让男人捅视频| 大香蕉久久网| 九色亚洲精品在线播放| 久久这里只有精品19| 在线观看免费视频网站a站| 啦啦啦在线免费观看视频4| 最近最新中文字幕大全免费视频| 12—13女人毛片做爰片一| 亚洲综合色网址| 亚洲国产日韩一区二区| 色老头精品视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲国产精品一区三区| 国产有黄有色有爽视频| 激情视频va一区二区三区| av片东京热男人的天堂| 欧美午夜高清在线| 咕卡用的链子| 两性夫妻黄色片| 国产日韩欧美在线精品| 人人妻人人添人人爽欧美一区卜| 亚洲av日韩在线播放| 午夜视频精品福利| 国产亚洲午夜精品一区二区久久| 91成年电影在线观看| 又紧又爽又黄一区二区| 国产成人免费观看mmmm| 国产伦理片在线播放av一区| 国产成人a∨麻豆精品| 亚洲精品自拍成人| 亚洲精品一区蜜桃| 国产精品免费大片| av有码第一页| 亚洲精品第二区| 久久久久久亚洲精品国产蜜桃av| 深夜精品福利| 欧美激情久久久久久爽电影 | 午夜91福利影院| 在线观看免费日韩欧美大片| 无遮挡黄片免费观看| 9色porny在线观看| 在线看a的网站| 国产成人欧美| 狂野欧美激情性xxxx| 欧美午夜高清在线| 又黄又粗又硬又大视频| 后天国语完整版免费观看| 人人妻,人人澡人人爽秒播| 亚洲av电影在线观看一区二区三区| 肉色欧美久久久久久久蜜桃| 少妇粗大呻吟视频| 国产激情久久老熟女| 成人免费观看视频高清| 亚洲专区国产一区二区| 日韩大片免费观看网站| 国产一区二区三区综合在线观看| 欧美精品亚洲一区二区| 99热网站在线观看| 91精品三级在线观看| 黄色视频在线播放观看不卡| 人人妻人人澡人人看| 宅男免费午夜| 欧美人与性动交α欧美软件| av福利片在线| 91麻豆精品激情在线观看国产 | 国产精品偷伦视频观看了| 天天躁狠狠躁夜夜躁狠狠躁| 日日夜夜操网爽| 中国美女看黄片| 国产一区有黄有色的免费视频| 日日爽夜夜爽网站| 18在线观看网站| 丝袜喷水一区| 一级毛片精品| 18在线观看网站| 啦啦啦免费观看视频1| 中文字幕高清在线视频| 亚洲精品久久久久久婷婷小说| 午夜影院在线不卡| 国产成人a∨麻豆精品| 国产精品偷伦视频观看了| 日韩电影二区| 搡老熟女国产l中国老女人| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看 | 国产成人一区二区三区免费视频网站| 黄片小视频在线播放| 美女国产高潮福利片在线看| 99久久精品国产亚洲精品| 一二三四社区在线视频社区8| 亚洲精品一卡2卡三卡4卡5卡 | 黑人巨大精品欧美一区二区蜜桃| 搡老乐熟女国产|