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

    Spectral characteristic of laser-induced plasma in soil

    2020-01-10 07:40:50YangYU余洋NanjingZHAO趙南京ZhigaoLAN蘭智高DeshuoMENG孟德碩andMingjunMA馬明俊
    Plasma Science and Technology 2020年1期
    關(guān)鍵詞:南京

    Yang YU (余洋), Nanjing ZHAO (趙南京), Zhigao LAN (蘭智高),Deshuo MENG (孟德碩) and Mingjun MA (馬明俊)

    1 College of Physics and Electronic Information, Huang Gang Normal University, Huanggang 438000,People's Republic of China

    2 Key Laboratory of Environmental Optics and Technology,Anhui Institute of Optics and Fine Mechanics,Chinese Academy of Sciences, Hefei 230031, People's Republic of China

    Abstract The spectral characteristic of laser-induced plasma in soil was studied in this work,laser-induced breakdown spectroscopy was used to analyze the spectral characteristic of plasma under the condition of different time delays and irradiances.Moreover,the time evolution characteristics of plasma temperature and electron density were discussed.Within the time delay range of 0-5 μs,the spectral intensity of the characteristic lines of Si I: 288.158 nm, Ti I: 336.126 nm, Al I:394.400 nm and Fe I: 438.354 nm of the four main elements in two kinds of national standard soil decayed exponentially with time. The average lifetime of the spectral lines was nearly 1.56 μs.Under the condition of different time delays,the spectral intensity of Pb I:405.78 nm in soil increased linearly with laser energy. However, the slope between the spectral intensity and laser energy decreased exponentially with the increase in time delay, from 4.91 to 0.99 during 0-5 μs. The plasma temperature was calculated by the Boltzmann plot method and the electron density was obtained by inversion of the full width at half maximum of the spectrum. The plasma temperature decreased from 8900 K to 7800 K and the electron density decreased from 1.5 × 1017 cm?3 to 7.8 × 1016 cm?3 in the range of 0-5 μs.

    Keywords: LIBS, soil, delay time, laser energy, plasma temperature, electron density(Some figures may appear in colour only in the online journal)

    1. Introduction

    With the acceleration of industrialization,discharges of heavy metals from electroplating, mining, chemical and mechanical manufacturing industries are becoming increasingly prominent. This has posed a serious threat to the ecological environment and human health.Rapid and accurate detection of heavy metals is key to the prevention and control of heavy metal discharge. At present, heavy metals in soil are commonly detected using laboratory methods such as inductively coupled plasma mass spectrometry[1]and atomic absorption spectrometry [2]; these methods have high sensitivity and accuracy for heavy metal detection. However, they also have some disadvantages, for example, complex sample pretreatment process,high time consumption,and expensive and high requirements for operators. Thus, real-time and fast analysis of heavy metals cannot be realized using these methods.

    Laser-induced breakdown spectroscopy(LIBS)[3,4]has several advantages such as fast analysis, no sample pretreatment, easy operation and analysis of different specimens. In recent years, researchers have carried out related research work on environment [5-9], steel [10-13], food [14-16],ocean [17-19], energy [20, 21], space exploration [22, 23],and biology [24, 25]. However, LIBS is also faced with difficulties of low quantitative accuracy and high detection limit.

    Spectral characteristic analysis is the basis of quantitative analysis in LIBS, which is closely related to the experimental parameters. Haider et al measured the radiative lifetime of excited neutral nitrogen atoms by time-resolved LIBS; the values were largely consistent with the calculated value of lifetime given by Hibbert [26]. Sirven et al used LIBS to measure standard metallic samples; the measurement was found to linearly depend on the pulse energy over the range of investigation, and this behavior was related to the increase of the number of vaporized atoms when the pulse energy was increased [27]. Dong et al used LIBS to analyze spectral characterization of nitrogen in farmland soil;the effects of laser energy on spectral characteristic and the time evolution characteristic of LIBS were studied to determine the basis of parameter optimization for the measurement[28].Pardini et al recorded Balmer lines of hydrogen under the condition of different time delays and evaluated the electron density of plasma through the measurement of the Stark broadening; the experimental results were compared with theoretical predictions [29]. However, investigation of the time evolution characteristic of laser-induced plasma in different types of soil has been seldom presented,and the variation law of slope between spectral intensity and laser energy has not yet been reported.

    In this study,the time evolution characteristic of LIBS in soil was investigated and the lifetime of characteristic lines were also determined. The influence of laser energy on the spectral characteristic was studied under the condition of different time delays,and the time evolution of slope between spectral intensity and laser energy was further discussed.The evolutions of plasma temperature and electron density of plasma were also analyzed in this study.

    2. Experimental system and samples

    2.1. Experimental system

    An Nd:YAG pulsed laser was used as the excitation source with a 1064 nm central wavelength and a 5 ns pulse width(Brilliant type, Quantel Company, France); its energy was adjustable. The pulse laser was reflected by the 1064 nm total reflector; a spot with high peak irradiance was produced by the convergence of a flat convex lens(focal length 50 mm), which acted on the surface of the sample instantaneously and produced a plasma. The spectral signal emitted by the plasma was coupled to the fiber by the convex lens (focal length 30 mm). The spectral signal was separated by a medium-step grating spectrometer (Mechell 5000, Andor Company, UK); the spectral detection range was 200-850 nm, and the spectral resolution was λ /Δλ = 5000.The signal was eventually detected by ICCD (intensified charge coupled device, i Star, Andor Company,UK) with a pixel size of ICCD is 1024 × 1024. The samples were placed on a rotating two-dimensional mobile platform; the time delay of spectral signal acquisition was controlled by a four-channel digital delay pulse generator(DG535, SRS Company, America), and the time series between the laser pulse and two-dimensional platform rotation was also controlled by DG535. The experimental system is shown in figure 1.

    2.2. Samples

    The standard solution was prepared with Pb(NO3)2and deionized water; samples containing Pb were prepared by mixing the standard solution with national standard soil(GBW-07405,GBW-07408). The samples were air-dried, grinded and sifted(100 mesh); each sample weighed 2.8 g as measured by an electronic balance and was pressed into a round cake with a diameter of 30 mm and a thickness of 2.6 mm.

    3. Results and discussion

    3.1. Spectral characteristic of plasma under the condition of different time delays

    Plasma was produced by the interaction between a pulsed laser with high peak irradiance and the sample; the plasma gradually decayed and disappeared. Four main elements, Si,Ti, Al, and Fe, in the national standard soil were selected as analytical objects. In the 0-5 μs of plasma evolution, the samples were measured five times under the condition of each time delay (with interval of 0.5 μs). The spectra of each measurement were accumulated by 10 laser pulses, and the final result was the average value of the five measurements.The time evolution characteristic of the spectral intensity of the four characteristic lines, Si I: 288.158 nm, Ti I:336.126 nm, Al I: 394.400 nm, and Fe I: 438.354 nm, are shown in figures 2(a)-(d).

    The results in figures 2(a)-(d) indicated that the spectral intensity of four characteristic lines, Si I: 288.158 nm, Ti I:336.126 nm,Al I:394.400 nm,and Fe I:438.354 nm,differed considerably in the initial stage.With gradual evolution of the characteristic lines, the difference between the spectral intensity of the characteristic lines also diminished gradually.The evolution law of spectral intensity was basically the same, and the spectral intensity of the characteristic lines decreased exponentially with the time delay. The spectral intensity I can be expressed as:

    Table 1.Lifetime of four characteristic lines in two national standard soils.

    where I0and A are constants, andτ is the lifetime of the characteristic line. The lifetime of spectral intensity of the four characteristic lines was determined by exponential decay function fitting; the corresponding results are shown in table 1.

    The above results showed that the lifetime of the characteristic spectral lines of the four main elements in the soil samples was different and ranged over 1.0-2.0 μs. For two different types of national standard soils,the average lifetime of the four characteristic lines was approximately 1.56 μs.This finding has important guiding significance for the determination of time delay in simultaneous measurement of multiple elements.

    3.2. Spectral characteristic of plasma under the condition of different irradiances

    In the measurement and analysis of LIBS, laser energy has a significant influence on the spectral characteristic of plasma radiation. The process of plasma excitation, emission and evolution is influenced by the laser energy.With the increase of laser energy, the number of ablated samples will increase;therefore, the spectral intensity of plasma radiation will increase. However, the evolution of spectral intensity of characteristic lines with laser energy under the condition of different time delays has not been thoroughly studied.Within the time delay of 0-5 μs, the samples were measured five times under the condition of each time delay(with intervals of 0.5 μs); the spectra of each measurement were accumulated by 10 laser pulses, and the final result was the average value of the five measurements. The dependence of the spectral intensity of the characteristic line Pb I:405.78 nm on the laser energy under the condition of different time delays is shown in figures 3(a)-(f).

    The results in figures 3(a)-(f)indicated that the spectral intensity of Pb I: 405.78 nm increased linearly with the increase of laser energy in the range of 50-280 mJ.Figure 4 shows the relationship between the slope of spectral intensity versus laser energy and time delay in the range of 0-5 μs.

    The results in figure 4 showed that the slope of spectral intensity of the characteristic line Pb I: 405.78 nm versus laser energy decreased exponentially with the increase of the time delay.The slope of spectral intensity of the characteristic line versus laser energy differed significantly for different time delays. The slope decreased from 4.91 to 0.99 during 0-5 μs. This illustrates that the change trend of spectral intensity with laser energy became smaller with the evolution of plasma, i.e., the influence of laser energy on the spectral characteristic of plasma radiation became weaker.The reason was that the spectral intensity of plasma radiation decreased significantly with the increase of time delay. The slope of spectral intensity versus laser energy was large during the initial period of plasma evolution. As the time delay increased, the variation of spectral intensity was relatively flat, i.e., the variation of spectral intensity with laser energy became smaller.

    3.3.Time evolution of plasma temperature and electron density

    The spectral characteristic of plasma radiation is closely related to the plasma temperature and electron density.Therefore, it is necessary to measure and study the evolution of plasma temperature and electron density for analyzing and understanding the spectral characteristic of plasma radiation.

    3.3.1. Plasma temperature. Under the assumption of local thermal equilibrium,plasma temperature can be calculated by the Boltzmann plot method as:

    In the above formula, I is the intensity of spectral line, h is the Planck constant,N is the density of particle number,g is the degeneracy factor of upper level, A is the transition probability,λ is the wavelength of spectral line, Z is the partition function,E is the energy of upper level,and k is the Boltzmann constant.

    Considering the wide distribution and large number of characteristic lines of Fe element in the samples, 13 characteristic lines of Fe in the wavelength range of 360-450 nm were selected for the analysis [30], the energy levels and related parameters were quoted from the website of NIST[3].Under the condition of laser energy 174.5 mJ,time delay td= 1.0 μs,and acquisition gate width tg= 0.5 μs,and average temperatures are measured with a gate-width of 0.5 μs.The Boltzmann diagram of Fe is shown in figure 5(a),and the time evolution characteristic of plasma temperature is shown in figures 5(b).

    3.3.2.Electron density. The broadening mechanism involves natural broadening,Doppler broadening and Stark broadening in LIBS plasma spectroscopy. However, the most important contribution to linewidth is Stark broadening (the instrument width is 0.01 nm,the Doppler width is 0.005 nm).The electric field causing Stark effect is mainly contributed by the collision between electrons. Therefore, Stark broadening is related to the electron density. Ignoring the effect of ion broadening, the relationship between the full width at half maximum of the Lorentz fitting curve, Δλ1/2,and electron density in plasma, Ne,is given by the following formula:

    In the above formula,ω is the electron collision parameter [31]. The characteristic line Ca I: 422.67 nm was fitted by the Lorentz function; the Lorentz fitting results and time evolution of electron density are shown in figures 6(a)-(b).

    The results shown in figures 5 and 6 indicated that the plasma temperature decreased linearly with time;it decreased from 8900 K to 7800 K in the range of 0-5 μs. Thus, the plasma temperature decreased by 1100 K(nearly 12.4%).The electron density decreased exponentially with time; it decreased from 1.5 × 1017cm?3to 7.8 × 1016cm?3within the range of 0 ~ 5 μs. Hence, the electron density decreased by approximately 50%. Plasma was produced by the interaction between pulsed laser and matter. The plasma expanded and cooled with time, and the degree of collision between particles in plasma became weaker. The number of particles per unit volume reduced, so plasma temperature decreased and the electron density became lesser. This was also the essential reason why the spectral intensity of plasma radiation decreased with time.

    4. Conclusions

    In this study, the time evolution of spectral intensity of two kinds of national soil was analyzed by LIBS; the average lifetime of spectral lines of four main elements was nearly 1.56 μs. This observation is helpful for the determination of time delay in the simultaneous measurement of multiple elements.The spectral intensity increased with the increase of laser energy over the range of 50-280 mJ (corresponding to peak irradiance of 3.18 × 1012-1.78 × 1013W cm?2; spot radiusω was approximately 10 μm order of magnitude by simulation[32],pulse widthΔt was 5 ns,and focal length of lens was 50 mm, the peak irradiance was calculated by the formula E / (Δ t · πω2)); the peak irradiance of laser(TW cm-2) was far beyond the threshold of optical breakdown of soil(GW cm-2).The slope of spectral intensity of Pb I:405.78 nm versus laser energy decreased from 4.91 to 0.99 during the time delay of 0-5 μs. This illustrated that the influence of irradiance on spectral characteristic became weaker as the time delay increased.Within 0-5 μs,the plasma temperature and electron density decreased by nearly 12.4%and 50%, respectively.

    Acknowledgments

    This work is supported by the Special Project of Knowledge Innovation of Hubei Province (No. 2018CFC862), National Natural Science Foundation of China (No. 61705238), the National High Technology Research and Development Program of China (No. 2014AA06A513).

    猜你喜歡
    南京
    南京比鄰
    “南京不會(huì)忘記”
    南京大闖關(guān)
    江蘇南京卷
    南京·九間堂
    金色年華(2017年8期)2017-06-21 09:35:27
    南京·鴻信云深處
    金色年華(2017年7期)2017-06-21 09:27:54
    南京院子
    電影(2017年1期)2017-06-15 16:28:04
    又是磷復(fù)會(huì) 又在大南京
    南京:誠(chéng)實(shí)書店開張
    南京、南京
    99精国产麻豆久久婷婷| 美女主播在线视频| 黑人猛操日本美女一级片| 欧美日韩国产mv在线观看视频 | 久久99热这里只有精品18| 综合色丁香网| a级毛色黄片| 日本欧美视频一区| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 青青草视频在线视频观看| 久久精品夜色国产| 国产高清三级在线| 夫妻性生交免费视频一级片| 男女下面进入的视频免费午夜| 亚洲精品国产色婷婷电影| 亚洲国产色片| 国产伦精品一区二区三区四那| 亚洲欧美精品自产自拍| 久久久成人免费电影| 久久99热这里只频精品6学生| 在线免费观看不下载黄p国产| 久久精品国产a三级三级三级| 成人高潮视频无遮挡免费网站| 一区二区三区精品91| 自拍欧美九色日韩亚洲蝌蚪91 | 尤物成人国产欧美一区二区三区| 熟女电影av网| 九九久久精品国产亚洲av麻豆| 亚洲,一卡二卡三卡| 在线免费观看不下载黄p国产| 青春草视频在线免费观看| 国产精品一区二区在线观看99| 3wmmmm亚洲av在线观看| 校园人妻丝袜中文字幕| 韩国av在线不卡| 日韩视频在线欧美| 日本免费在线观看一区| 成人国产av品久久久| 我的女老师完整版在线观看| 我要看黄色一级片免费的| 国产日韩欧美亚洲二区| 大陆偷拍与自拍| 国产伦精品一区二区三区视频9| 亚洲欧美日韩卡通动漫| 纯流量卡能插随身wifi吗| 免费看光身美女| 精品人妻视频免费看| 亚洲av电影在线观看一区二区三区| 国产美女午夜福利| 人妻 亚洲 视频| 永久免费av网站大全| xxx大片免费视频| 一本一本综合久久| 在线观看一区二区三区激情| 寂寞人妻少妇视频99o| 日韩大片免费观看网站| 美女cb高潮喷水在线观看| 新久久久久国产一级毛片| 黑丝袜美女国产一区| 国产精品秋霞免费鲁丝片| 亚洲高清免费不卡视频| 国产乱人偷精品视频| 国产男女超爽视频在线观看| 日韩中文字幕视频在线看片 | 亚洲欧美精品专区久久| 18禁在线播放成人免费| 超碰97精品在线观看| 麻豆成人av视频| 国产真实伦视频高清在线观看| 青青草视频在线视频观看| 超碰97精品在线观看| 一个人看视频在线观看www免费| 国产精品欧美亚洲77777| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 国产黄频视频在线观看| 观看美女的网站| 极品教师在线视频| 国产精品.久久久| 少妇人妻 视频| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 亚洲一区二区三区欧美精品| 少妇裸体淫交视频免费看高清| 欧美3d第一页| 女人久久www免费人成看片| 国产精品秋霞免费鲁丝片| 精品久久久久久久末码| 老师上课跳d突然被开到最大视频| 一级毛片黄色毛片免费观看视频| 亚洲精品国产色婷婷电影| 国产亚洲一区二区精品| 人妻制服诱惑在线中文字幕| 亚洲一级一片aⅴ在线观看| 一本一本综合久久| 亚洲av中文字字幕乱码综合| 日韩中字成人| 一区二区三区免费毛片| 久久久a久久爽久久v久久| 高清黄色对白视频在线免费看 | 日韩一区二区三区影片| 欧美人与善性xxx| 3wmmmm亚洲av在线观看| 亚洲国产最新在线播放| 日本欧美视频一区| 亚洲av国产av综合av卡| 视频中文字幕在线观看| 国产91av在线免费观看| 欧美区成人在线视频| 各种免费的搞黄视频| 麻豆成人午夜福利视频| 国产精品人妻久久久久久| 亚洲美女黄色视频免费看| 啦啦啦中文免费视频观看日本| 熟妇人妻不卡中文字幕| 夜夜看夜夜爽夜夜摸| 久久久久视频综合| 大片免费播放器 马上看| 亚洲人成网站在线播| 中文资源天堂在线| 日本免费在线观看一区| 日本色播在线视频| 国产亚洲欧美精品永久| 国产av国产精品国产| 97精品久久久久久久久久精品| 一个人免费看片子| 国产一区二区三区综合在线观看 | 边亲边吃奶的免费视频| 日本黄大片高清| 日韩av在线免费看完整版不卡| 亚洲精品日韩在线中文字幕| 3wmmmm亚洲av在线观看| 亚洲婷婷狠狠爱综合网| 99热这里只有是精品50| 日韩人妻高清精品专区| 人人妻人人添人人爽欧美一区卜 | av播播在线观看一区| 国产高清有码在线观看视频| 亚洲第一av免费看| 亚洲av中文av极速乱| 人人妻人人添人人爽欧美一区卜 | 国产精品无大码| 国产免费视频播放在线视频| 中文在线观看免费www的网站| 成人毛片60女人毛片免费| 国产亚洲欧美精品永久| 老司机影院成人| 一级爰片在线观看| 黄色日韩在线| 免费观看a级毛片全部| 99九九线精品视频在线观看视频| 国产精品熟女久久久久浪| 又爽又黄a免费视频| 国产探花极品一区二区| 久久精品人妻少妇| 国产男女内射视频| 最后的刺客免费高清国语| 国产精品伦人一区二区| 蜜桃在线观看..| 交换朋友夫妻互换小说| 亚洲av男天堂| 久久 成人 亚洲| 美女视频免费永久观看网站| 人妻一区二区av| 22中文网久久字幕| 大码成人一级视频| 免费观看的影片在线观看| av专区在线播放| 国产一区二区三区av在线| 欧美日韩国产mv在线观看视频 | 久久人人爽av亚洲精品天堂 | 国产精品欧美亚洲77777| 男人狂女人下面高潮的视频| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级 | 2018国产大陆天天弄谢| 色视频www国产| 一个人看的www免费观看视频| 日韩中文字幕视频在线看片 | 美女国产视频在线观看| 少妇人妻精品综合一区二区| 国产av码专区亚洲av| 国产精品欧美亚洲77777| 亚洲av不卡在线观看| av.在线天堂| 丰满迷人的少妇在线观看| 直男gayav资源| 日本黄大片高清| 一级毛片久久久久久久久女| 热re99久久精品国产66热6| 成年女人在线观看亚洲视频| 男的添女的下面高潮视频| 十分钟在线观看高清视频www | 亚洲内射少妇av| 国产在视频线精品| av黄色大香蕉| 免费观看av网站的网址| 日日撸夜夜添| 久久国内精品自在自线图片| 免费少妇av软件| 成人综合一区亚洲| 美女脱内裤让男人舔精品视频| 欧美日韩精品成人综合77777| av一本久久久久| 在现免费观看毛片| 涩涩av久久男人的天堂| h日本视频在线播放| 亚洲真实伦在线观看| 精品久久国产蜜桃| 国产精品熟女久久久久浪| av黄色大香蕉| 亚洲欧美精品自产自拍| 一本色道久久久久久精品综合| 久久97久久精品| 干丝袜人妻中文字幕| 另类亚洲欧美激情| 亚洲精品aⅴ在线观看| 亚洲精品乱久久久久久| 欧美丝袜亚洲另类| av.在线天堂| 国产综合精华液| tube8黄色片| 日本与韩国留学比较| 亚洲综合精品二区| 99re6热这里在线精品视频| 国产精品人妻久久久影院| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 久久人人爽av亚洲精品天堂 | 青春草亚洲视频在线观看| 国产乱人偷精品视频| 国产极品天堂在线| 最后的刺客免费高清国语| 尾随美女入室| 九色成人免费人妻av| 国产欧美日韩一区二区三区在线 | 夫妻性生交免费视频一级片| 欧美一区二区亚洲| tube8黄色片| 日本色播在线视频| 人妻少妇偷人精品九色| 亚洲色图av天堂| 91久久精品国产一区二区三区| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 亚洲av免费高清在线观看| 国产成人freesex在线| 自拍欧美九色日韩亚洲蝌蚪91 | 精品酒店卫生间| 有码 亚洲区| 亚洲精品一区蜜桃| 汤姆久久久久久久影院中文字幕| 在线免费十八禁| 伦理电影大哥的女人| 夫妻午夜视频| 国产亚洲av片在线观看秒播厂| 亚洲美女黄色视频免费看| 色婷婷久久久亚洲欧美| 一个人看的www免费观看视频| 成人影院久久| 一个人免费看片子| 少妇丰满av| 免费在线观看成人毛片| 一级毛片黄色毛片免费观看视频| 啦啦啦视频在线资源免费观看| 99久久精品一区二区三区| 日韩伦理黄色片| 国产美女午夜福利| 熟女电影av网| 大片免费播放器 马上看| 国产精品不卡视频一区二区| 精品久久国产蜜桃| 老女人水多毛片| 3wmmmm亚洲av在线观看| 少妇熟女欧美另类| 亚洲av日韩在线播放| 美女福利国产在线 | 啦啦啦在线观看免费高清www| 极品教师在线视频| 免费高清在线观看视频在线观看| 国产美女午夜福利| 久久人妻熟女aⅴ| 欧美zozozo另类| 在线播放无遮挡| 在线观看三级黄色| av在线app专区| av免费观看日本| 51国产日韩欧美| 国产有黄有色有爽视频| 欧美成人a在线观看| 亚洲欧美成人精品一区二区| 国产精品欧美亚洲77777| 老女人水多毛片| 久久精品久久久久久噜噜老黄| 99热全是精品| 国产乱人视频| 久久久久精品性色| 久久人妻熟女aⅴ| 天天躁日日操中文字幕| 久久久久久久久久成人| 男女国产视频网站| 成人一区二区视频在线观看| 九九久久精品国产亚洲av麻豆| 老司机影院成人| 免费看不卡的av| 国产 一区精品| 欧美zozozo另类| 乱码一卡2卡4卡精品| 国产女主播在线喷水免费视频网站| a级毛片免费高清观看在线播放| 久久人人爽人人片av| 伊人久久精品亚洲午夜| 国产成人91sexporn| 久久久色成人| 伊人久久精品亚洲午夜| 成人国产av品久久久| 男人爽女人下面视频在线观看| 欧美97在线视频| 欧美亚洲 丝袜 人妻 在线| 亚洲,一卡二卡三卡| 观看免费一级毛片| 午夜福利在线在线| 国产免费一区二区三区四区乱码| 国产成人午夜福利电影在线观看| 色婷婷av一区二区三区视频| 国产精品国产三级国产专区5o| 天天躁夜夜躁狠狠久久av| 国产老妇伦熟女老妇高清| 你懂的网址亚洲精品在线观看| 欧美成人a在线观看| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 岛国毛片在线播放| 亚洲色图综合在线观看| 人妻制服诱惑在线中文字幕| 日日摸夜夜添夜夜爱| 国产在线视频一区二区| 国产高清有码在线观看视频| 五月开心婷婷网| 欧美zozozo另类| 久久久久国产网址| av在线播放精品| 国产一区二区三区综合在线观看 | 日本免费在线观看一区| 午夜福利影视在线免费观看| 永久网站在线| 国产片特级美女逼逼视频| 国产成人精品久久久久久| 18禁裸乳无遮挡免费网站照片| 一个人看的www免费观看视频| 91精品国产国语对白视频| 国产在线免费精品| 国产成人免费无遮挡视频| 午夜免费观看性视频| 中文精品一卡2卡3卡4更新| 国产亚洲av片在线观看秒播厂| 亚洲综合色惰| 少妇精品久久久久久久| 亚洲精品乱久久久久久| 成人美女网站在线观看视频| 免费久久久久久久精品成人欧美视频 | 91精品一卡2卡3卡4卡| 亚洲精品国产色婷婷电影| 久久久欧美国产精品| av一本久久久久| 美女中出高潮动态图| 大码成人一级视频| 免费人成在线观看视频色| 2018国产大陆天天弄谢| 亚洲av中文字字幕乱码综合| 亚洲国产欧美人成| videos熟女内射| 一区在线观看完整版| 99re6热这里在线精品视频| 欧美xxxx黑人xx丫x性爽| 大又大粗又爽又黄少妇毛片口| 色婷婷av一区二区三区视频| 精品99又大又爽又粗少妇毛片| 成人毛片a级毛片在线播放| 美女cb高潮喷水在线观看| 熟女av电影| 美女主播在线视频| 久久久色成人| 成人高潮视频无遮挡免费网站| 久久人人爽人人爽人人片va| 国产亚洲av片在线观看秒播厂| 欧美日韩视频精品一区| 国产淫片久久久久久久久| 亚洲久久久国产精品| 一边亲一边摸免费视频| 蜜桃久久精品国产亚洲av| 婷婷色麻豆天堂久久| 最近手机中文字幕大全| 在线精品无人区一区二区三 | 免费观看无遮挡的男女| 人人妻人人看人人澡| 国产男女内射视频| 午夜精品国产一区二区电影| 午夜福利高清视频| 日日摸夜夜添夜夜爱| 乱码一卡2卡4卡精品| 婷婷色麻豆天堂久久| 丰满少妇做爰视频| 高清视频免费观看一区二区| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 中文字幕av成人在线电影| 午夜老司机福利剧场| 精品少妇黑人巨大在线播放| 国国产精品蜜臀av免费| 九九爱精品视频在线观看| 汤姆久久久久久久影院中文字幕| 国产精品嫩草影院av在线观看| h日本视频在线播放| 18禁在线播放成人免费| 成人综合一区亚洲| 99热网站在线观看| 欧美 日韩 精品 国产| 午夜日本视频在线| 五月天丁香电影| 插逼视频在线观看| 日韩欧美精品免费久久| 三级国产精品欧美在线观看| 久久久久国产精品人妻一区二区| av播播在线观看一区| 18禁动态无遮挡网站| 精品一区二区免费观看| 一级爰片在线观看| 国产精品爽爽va在线观看网站| 欧美+日韩+精品| 欧美zozozo另类| 国产av精品麻豆| 女性被躁到高潮视频| 男女边摸边吃奶| 成人二区视频| 久久久久久伊人网av| 中文字幕久久专区| 亚洲av电影在线观看一区二区三区| 日本黄色日本黄色录像| 成人无遮挡网站| 日本黄大片高清| 国产中年淑女户外野战色| 免费观看a级毛片全部| 久久精品国产a三级三级三级| 亚州av有码| 国产午夜精品一二区理论片| 黄色怎么调成土黄色| 国产男女超爽视频在线观看| 在线 av 中文字幕| 国产午夜精品久久久久久一区二区三区| 欧美成人精品欧美一级黄| 亚洲综合色惰| 午夜视频国产福利| 天美传媒精品一区二区| 熟女电影av网| 一级黄片播放器| 黄片wwwwww| av在线app专区| 国产精品免费大片| 国产成人freesex在线| 最后的刺客免费高清国语| 边亲边吃奶的免费视频| 国产精品伦人一区二区| 日韩大片免费观看网站| 亚洲精华国产精华液的使用体验| 伦理电影大哥的女人| 国产精品爽爽va在线观看网站| 日韩,欧美,国产一区二区三区| 日韩欧美精品免费久久| 男女边摸边吃奶| 18禁动态无遮挡网站| 丰满人妻一区二区三区视频av| 午夜免费观看性视频| 黄色怎么调成土黄色| 99精国产麻豆久久婷婷| 久久精品夜色国产| 最新中文字幕久久久久| 联通29元200g的流量卡| 偷拍熟女少妇极品色| 波野结衣二区三区在线| 国产人妻一区二区三区在| av播播在线观看一区| 久久精品国产a三级三级三级| 一级片'在线观看视频| 亚洲性久久影院| 欧美日韩国产mv在线观看视频 | 亚洲精品456在线播放app| 日日撸夜夜添| 亚洲精华国产精华液的使用体验| 久久 成人 亚洲| 身体一侧抽搐| 亚洲精品456在线播放app| 中国三级夫妇交换| 亚洲欧美成人精品一区二区| 亚洲婷婷狠狠爱综合网| 男女边摸边吃奶| 国产无遮挡羞羞视频在线观看| 国产精品秋霞免费鲁丝片| 亚洲丝袜综合中文字幕| 欧美zozozo另类| av在线蜜桃| 妹子高潮喷水视频| 色吧在线观看| 亚洲精品一二三| videossex国产| 亚洲av欧美aⅴ国产| 久久精品久久久久久久性| 99热网站在线观看| 欧美老熟妇乱子伦牲交| 精品久久久精品久久久| 免费高清在线观看视频在线观看| av黄色大香蕉| 青春草亚洲视频在线观看| 久久女婷五月综合色啪小说| 国产成人aa在线观看| 最新中文字幕久久久久| 男人添女人高潮全过程视频| 久久97久久精品| 日本黄大片高清| 男人狂女人下面高潮的视频| 99热网站在线观看| 国产午夜精品一二区理论片| 美女内射精品一级片tv| 亚洲精品自拍成人| 蜜桃久久精品国产亚洲av| 亚洲国产欧美人成| 人人妻人人澡人人爽人人夜夜| 黄色视频在线播放观看不卡| 欧美bdsm另类| 免费人成在线观看视频色| 久久6这里有精品| 两个人的视频大全免费| 国产v大片淫在线免费观看| 亚洲av福利一区| 免费看光身美女| 成人亚洲欧美一区二区av| 久久久久国产精品人妻一区二区| 久久99热这里只有精品18| 99热这里只有是精品50| 久久精品久久精品一区二区三区| 亚洲精品aⅴ在线观看| 久久久久久久国产电影| 国产精品免费大片| 三级经典国产精品| 99热6这里只有精品| 99久久精品热视频| 另类亚洲欧美激情| 性色avwww在线观看| 国产极品天堂在线| 极品少妇高潮喷水抽搐| 亚洲av电影在线观看一区二区三区| 国产av码专区亚洲av| 欧美3d第一页| 久久人人爽人人爽人人片va| 国产精品.久久久| 黄色视频在线播放观看不卡| 亚洲怡红院男人天堂| 亚洲精华国产精华液的使用体验| 日韩免费高清中文字幕av| 国产一区二区三区av在线| 成人综合一区亚洲| 久久韩国三级中文字幕| 汤姆久久久久久久影院中文字幕| 免费黄频网站在线观看国产| 嫩草影院新地址| 欧美精品人与动牲交sv欧美| 久久精品国产鲁丝片午夜精品| 中国国产av一级| 欧美xxⅹ黑人| 国产91av在线免费观看| 日韩强制内射视频| 久久人妻熟女aⅴ| 中文字幕制服av| 丰满乱子伦码专区| 全区人妻精品视频| 欧美日韩亚洲高清精品| 男女国产视频网站| 欧美日韩一区二区视频在线观看视频在线| 免费看日本二区| 男女边吃奶边做爰视频| 最近中文字幕2019免费版| 国产精品久久久久久精品古装| 一边亲一边摸免费视频| 老司机影院毛片| 亚洲av综合色区一区| 国产成人精品婷婷| 波野结衣二区三区在线| 热re99久久精品国产66热6| 啦啦啦在线观看免费高清www| 国产白丝娇喘喷水9色精品| 97精品久久久久久久久久精品| 最近最新中文字幕免费大全7| 伊人久久精品亚洲午夜| 狂野欧美激情性bbbbbb| 18禁裸乳无遮挡动漫免费视频| 久久久久国产网址| 少妇熟女欧美另类| 在线观看三级黄色| 波野结衣二区三区在线| 免费看av在线观看网站| 亚洲欧洲日产国产| 欧美最新免费一区二区三区| 人妻系列 视频| 大片电影免费在线观看免费| 久久国产精品大桥未久av | 久久久久久久久久成人| 国产 精品1| 一本色道久久久久久精品综合| 国产成人精品福利久久| 麻豆精品久久久久久蜜桃| 国产成人精品婷婷| 大陆偷拍与自拍| 久久99热6这里只有精品| 日韩av免费高清视频|