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

    Evaluation of gases’ molecular abundance ratio by fiber-optic laser-induced breakdown spectroscopy with a metal-assisted method

    2022-01-10 14:51:50YuZHU朱瑜PingweiZHOU周平偉andShengfuLI李生福
    Plasma Science and Technology 2021年12期

    Yu ZHU (朱瑜), Pingwei ZHOU (周平偉) and Shengfu LI (李生福)

    Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, People’s Republic of China

    Abstract A metal-assisted method is proposed for the evaluation of gases’ molecular abundance ratio in fiber-optic laser-induced breakdown spectroscopy(FO-LIBS).This method can reduce the laser ablation energy and make gas composition identification possible.The principle comes from the collision between the detected gases and the plasma produced by the laser ablation of the metal substrate.The interparticle collision in the plasma plume leads to gas molecules dissociating and sparking, which can be used to determine the gas composition.The quantitative relationship between spectral line intensity and molecular abundance ratio was developed over a large molecular abundance ratio range.The influence of laser ablation energy and substrate material on gas quantitative calibration measurement is also analyzed.The proposed metal-assisted method makes the measurement of gases’molecular abundance ratios possible with an FO-LIBS system.

    Keywords: gases’ molecular abundance ratio, collision of plasma, fiber optical laser induced breakdown spectroscopy, metal assisted method

    1.Introduction

    Laser-induced breakdown spectroscopy (LIBS) refers to focusing a pulsed laser onto the surface of target to create a high-temperature plasma.The emission spectrum of plasma can provide information about the elemental composition of the target[1–3].LIBS has been used for element analysis in a wide variety of fields, such as material analysis [4, 5],environmental monitoring[6],biological identification[7,8],forensics [9], and security, due to the advantages of online,in situ, and simultaneous multi-element analysis [10].It is a potential method for atomic spectrum analysis for species diagnosis.Several books and review articles have provided a complete description of its theoretical foundations and recent applications in industrial and analytical situations [11–13].

    The measurement of gases’ molecular abundance ratios has garnered great attention in combustion and energy-related research.The optical methods are of great importance due to their effectiveness in many situations where intrusive measurements are inconvenient.The real-time,in situnature of LIBS techniques makes them suitable to many applications requiring online data acquisition [14–16].The LIBS technique has been used to investigate gases’molecular abundance ratio in gas mixtures such as the in-air detection of Cl[17],F[18],and pollutant concentrations [19], as well as the deuterium/hydrogen ratio in gas reaction products[20].Of particular note,it has also been used for equivalence ratio measurement in a sparkignited engines and methane–air premixed flames,which arises from the desire to quantify real-time, localized combustion performance [21–24].Although LIBS has been considered a promising method for elemental analysis in gaseous detection,the conventional LIBS system (without using optical fibers to deliver the laser pulses) is limited when the environment is too hazardous or access space is restricted [25].Fiber-optic LIBS(FO-LIBS)systems,in which the laser beam is transferred to the target by a transmission fiber, provide a solution to the above problem.Researchers have developed many fiber-coupled laser-induced breakdown spectroscopy instruments for analysis of target elements[26–28].However,compared with traditional LIBS technology,the laser energy delivered to the target surface is limited by the aperture and the damage threshold of the fiber,which is lower than the breakdown threshold of the gas[29,30].In FO-LIBS, it is impossible to detect the gas by directly focusing the laser pulse delivered by the optical fiber on the gas,while increasing the incident laser energy would lead to inevitable optical fiber damage [31].Using optical fiber LIBS to realize gas composition identification is very important for harsh or narrow environment gas identification; however, to our knowledge,there has been no report on the measurement of gas composition using FO-LIBS.

    It is worth mentioning that surface-enhanced LIBS,based on a metal substrate,which was proposed by M A Aguirre to improve the detection limit of LIBS with the assistance of a metal substrate, exhibits superior spectral enhancement and repeatability [32, 33].The strong collision between the plasma and the detected particles in the plasma plume results in enhanced luminescence of the detected particles [34, 35].This gives us inspiration that the gas composition detection could be realized below the gas breakdown threshold in FOLIBS with the assistance of a metal substrate.The excitation radiation induced by the collision between metal particles in plasma and surrounding gas molecules could provide information about the elemental composition of the surrounding gaseous sample.

    In the current work, we propose a metal-assisted FOLIBS method to determine gas composition, and we take a methane mixture as an example to analyze the molecular abundance ratio over a wide range of mixture fractions.Atomic emissions from laser-induced plasma were observed and the ratios of the elemental line intensity were used to infer the molecular abundance ratio in mixture fractions.The advantage of this method is that gas detection could be realized using FO-LIBS as the laser energy transmitted by the fiber is lower than the breakdown threshold of the gas, while the fiber is not damaged.Compared to the traditional gas detection LIBS method,the metal substrate-assisted FO-LIBS is effective when the environment is too hazardous or the access space is restricted.Meanwhile,metal substrate-assisted FO-LIBS expands the field of applications of FO-LIBS from solid detection to gas detection [26, 31].

    2.Experiment

    For gas detection, a sample chamber was designed as shown in figure 1(a).Before filling, the sample chamber was evacuated to 0.01 Pa by a two-stage rotary vane vacuum pump and the gas pressure was measured by a capacitance pressure gauge.We first inflated a single gas component N2in the sample chamber at a pressureP1, and refilled the chamber with the second gas component CH4to the last immobile gas pressureP2=90 kPa; the final gas pressure was slightly below atmospheric pressure in order to prevent mingling of atmosphere.The gas flow of all experiments was set at 10 l min?1, and the purity of N2and CH4used in the experiment was more than 99.99%.Follow-up spectral testing will not be carried out until the two gases are mixed uniformly, and the waiting time is about 1 min.

    The schematic diagram of the metal-assisted FO-LIBS system for measurement of gases’ molecular abundance ratios is shown in figure 1(a)and the principle schematic diagram of the metal-assisted FO-LIBS method was shown in figure 1(b).A thin aluminum plate was fixed located near the chamber window as the metal substrate.The metal substrate was ablated using a Q-switch Nd:YAG laser (pulse width: 9 ns; repetition rate: 10 Hz; wavelength: 1064 nm).The laser beams were coupled into a 3 m long optical fiber using a single lens with a focal length of 100 mm.A fused silica fiber is used.The core diameter of the optical fiber is 1 mm.At the distal end of the fiber, the output beam from the fiber is focused on the metal substrate surface by a probe head.The probe consists of two lenses with a focal length of 100 mm and 50 mm respectively,which is used to focus the laser beam on the surface of metal substrate through the 10 mm thick chamber window.The probe head is mounted outside the window.The distance between the probe and the metal substrate is about 50 mm.

    As shown in table 1,the delivery efficiency of the fiber as the coupled laser energy ranges from 20 mJ to 44 mJ is about 86%.The maximum output energy at the distal end of the fiber probe is approximate 40 mJ,which was measured at the output end of the 50 mm focal length lens.To protect the fiber from damage,the largest energy used in this work was 30 mJ;no fiber damage was observed after more than 5 min of laser transmission.The emission from the plasma was collected in the optical fiber by the probe head, reflected by a dichroic mirror to the spectrograph slit, and then analyzed by the intensified charge-coupled device (ICCD) detector.The spectrometer equipped with a home built 1024×1024-pixel intensified charge-coupled device (ICCD) detector analyzed the plasma spectroscopy.The grating used here was 150 grooves/mm and the spectrometer has a focus length of 750 mm.The detection system has a spectral range of approximately 100 nm with a resolution of 0.35 nm.The experiment timing includes a laser Q-switch,and ICCD delay time control was performed through a TTL signal generator.The ICCD detector started signal integration with a delay time 400 ns relative to the ablation of aluminum sample, and the gate time was set to 1000 ns.A spectrum was recorded for each laser shot and 10 spectrum data were recorded to analyze the relative standard deviation (RSD).

    3.Results and discussion

    3.1.Metal-assisted FO-LIBS used for gas detection

    We first compare the gas detection capability of the direct focusing and metal substrate-assisted methods.The fractional molecular equilibrium ratios of nitrogen and methane are 12.1%and 87.9%respectively,while the concentrations error is about 1%.The laser energy for both was 30 mJ at the front of the probe head.The spectrum of the laser direct focusing method is shown as the solid line in figure 2(a),which just is a straight line and does not exhibit any spectral peak signal.The gases at the focus are not excited by the laser or transformed into plasma because the laser energy is lower than the gas breakdown threshold,while increasing the laser energy would lead to the fiber damage.The dashed solid line in figure 2(a)presents the spectrum of the gas sample with metal substrateassisted FO-LIBS.The spectrum ranging from 780 nm to 860 nm contains the line spectral of N 818.8 nm,N 821.6 nm,N 824.2 nm, C 833.5 nm, O 846 nm, N 856.7 nm, N 859.4 nm, and N 862.9 nm.The substrate used in the experiment is a thin aluminum bulk and the surface of the aluminum bulk was simply polished by sandpaper following cleaning by alcohol and deionized water.Compared with the direct focusing method, the metal-assisted method can obtain a clear gas composition spectrum at the same laser energy.It is obvious that metal-assisted method makes gas detection possible with the FO-LIBS system.

    The ablation trace of the metal substrate is measured as the size of the laser focus spot,which is 400 μm.The laser energy under the focusing lens is 30 mJ, corresponded to an average power 2.6 GW cm?2.While the ablated threshold of aluminum is about 1–5 J cm?2for a spot size of approximately 160 μm,the corresponding power is 0.1–0.5 GW cm?2[36] and the effective breakdown threshold of gas reaches approximately 100–1000 GW cm?2[37].The laser pulse energy used in this experiment is larger than the ablated threshold of aluminum but lower than the gas breakdown threshold.We were able to observe the atom line spectrum, shown as solid arrows in figure 2(a),when undertaking metal-assisted FO-LIBS with an aluminum substrate.However, we could not observe an atom line spectrum, even the continuous spectral, during direct focusing.The difference between direct focusing and metal substrate-assisted methods may be explained as follows.The laser energy is not enough to induce the gas breakdown into the plasma state in direct gas ablation,whereas the plasma particles produced by the ablation of the aluminum bulk substrate collided with the surrounding detected gas molecules and induced the gas breakdown.The collision between the detected gases particles and fast particles, include ions and electrons in the aluminum plasma, could lead to gas molecules dissociating to excited atoms and entering an ion state, as the schematic diagram in figure 1(b)shows.As the plasma cools,excited atoms and ions relax back into their ground states,emitting light at the characteristic atomic frequencies.Therefore, gas component detection could be achieved as the laser energy is lower than the breakdown threshold in the metal-assisted FO-LIBS system.

    Comparisons between the different mixing ratio gases were useful for molecular abundance ratio detection.The spectral comparison is shown in figure 2(b).The molecular abundance ratios of CH4are 0, 12.1%, and 22.7% respectively, and the total pressure was the same for all situations.When the CH4molecular abundance ratios rise from 0% to 22.7%, a significant increase of the C 833 nm spectral line intensity was noticed as the internal magnificent plate showed, while the N emission decreased as the CH4molecular abundance ratios rise.In all mixed sample tests,the laser energy and other parameters are the same, so the plasma temperature and electron concentration should be equal.The number of total particles that collide with plasma particles(include atom, ion and electron) is constant.The number of excited C atoms increases as the density of the CH4molecule increases, enhancing the collision probability between CH4molecule and aluminum particles when the CH4molecular abundance ratios rise.Finally, the C atom spectral line intensity increases and the N atom spectral line intensity decreases as the CH4molecular abundance ratios rises.The quantitative molecular abundance ratio measurement could be achieved based on the above process.

    3.2.Quantitative relationships between spectral line intensity and molecular abundance ratio

    The potential contamination presented on the metal substrate surface and the gas filling progress is inevitable.To eliminate the interference of carbon and nitrogen atoms sticking to the metal substrate surface and mixing in the gas filling progress,we firstly addressed the contamination quantity of carbon and nitrogen atoms before the analysis of quantitative relationships between spectral line intensity and molecular abundance ratio.The spectrum of the metal substrate in pure N2and CH4surroundings was shown in figure 3.The solid line was the metal substrate spectrum of pure CH4,where the spectral line intensity of N was used for background elimination of nitrogen atoms.The dashed line was the metal substrate spectrum of pure N2, where the spectral line intensity of C was used for background elimination of carbon atoms.To obtain the atom line spectral intensity of the detected gases,the nearby spectrum peak C 833 nm and N 818 nm at different CH4molecular abundance ratio were fitted by the single peak and multi-peak Lorentz functions; then the spectral peak intensities of C 833 nm and N 818 nm were extracted from the fitting parameter.Finally, the spectral peak intensity that represented the surrounding gases was obtained from subtracting the above fitting spectral peak intensity from the background C and N spectral intensity.

    To study the quantitative relationship between spectral line intensity and molecular abundance ratio, the mean values and standard deviations of the spectral peak intensity ratios were calculated and given in the calibration curves as shown in figure 4(a),where CH4molecular abundance ratios range from 0%to 22.7%.A function y =b+a·b·exp (-s·x) ,wherexdenotes the relatively molecular abundance ratio of CH4/N2,was fitted to the data points in the curve of figure 4(a).This fitting leads to the C/N calibration curve versus the molecular abundance ratio of CH4/N2, approximately corresponding toIC833/N818=0.005 +0.3·x·exp (1.25x).Meanwhile, for the low CH4abundance ratio range, the molecular abundance ratio of CH4and N2could be fitted to the linear functiony=0.00446+0.368xas the CH4abundance ratios range from 1%to 6%,and the fitting result was shown in figure 4(b).Just like most trace element detection, there exist lower concentration limits,which are known as the linear dynamic range[38].The calibration curves at high concentrations deviated from linearity,and the most likely reason is the self-absorption phenomenon,which comes from the fact that the emissions of the central hotter regions are absorbed by the colder atoms surrounding the plasma [39].

    In order to clearly understand the influence of the gas concentration on quantitative concentration measurements,the C 833 nm and N 818 nm spectral line intensities were shown in figures 4(c) and (d) as the molecular abundance ratios increasing at three laser pulse energy levels.It is clear that the evolution trend is different between the C 833 nm and N 818 nm spectral line intensities as the fractional molecular abundance ratio increases.The C 833 nm spectral line intensity showed a saturation phenomenon, which becomes more and more obvious with the decrease of the incident laser pulse energy,which is evidently self-absorption.While the N 818 nm spectral line intensity showed an exponential-like function increase as the N2molecular abundance ratio increased, this seems to be an avalanche mechanism.The collision between the plasma particles produced by aluminum ablation and the N2molecule induced the latter to the excited N atom and ion state.Excited atoms and ions get enough energy that they could even re-collide and re-ionize other ground state N2molecules.Finally,more CH4molecules lead to the increment ofIC833stagnating due to the self-absorption phenomenon, while more N2molecules lead to an increasingIN818avalanche due to secondary particle collision between N atoms and ions.The calibration curves deviate from a linear function at larger abundance ratios,which is the joint result of the C self-absorption phenomenon and the N avalanche mechanism.

    3.3.The incident laser pulse energy

    The incident laser pulse energy is a key factor in quantitative concentration measurements in the experimental calibration of LIBS.Figure 5 illustrates three calibration curves of CH4and N2mixture for incremental laser pulse energy from 20 mJ to 30 mJ with an aluminum plate substrate.The spectral line intensity ratio between C 833 nm and N 818 nm is extracted for analysis of the molecular abundance ratios in the CH4and N2gas mixture.Three calibration curves show a similar trend:the peak valueIC833/IN818positive correlation with the fractional molecular abundance ratio of CH4/N2.In a similar manner,we fitted the C/N calibration curve with the functionIC833/N818=b+a·x·exp (s·x),wherebdenotes the background of plasma luminescence,arepresents the linear relation between intensityIC833/N818and the fractional molecular abundance ratio of CH4/N2, andsrepresents the degree of the deviation from linearity.The fitted parametersis 2.25,1.93,and 1.51 as the laser energy is 20 mJ,25 mJ,and 30 mJ, respectively.The fitting curve presents a lager curvature as the incident laser pulse decreases, which indicates that the degree of the deviation from linearity is weakening along with the increase of the incident laser pulse energy.If the linear relationship between the spectral line intensity and molecular immunity ratio is satisfied,the molecular immunity ratio of gas can be inferred directly from spectral line intensity without considering the influence of self-absorption or other effects.To acquire suitable linear quantitative concentration measurements of gases in metal-assisted FO-LIBS, the incident laser pulse energy should be increased to produce more plasma particles at higher plasma temperatures.

    Figure 1.(a) Schematic diagram of metal-assisted FO-LIBS for measurement of gases’ molecular abundance ratios.The sample chamber connects with the CH4 and N2 mixing gas tank.The laser pulse output from the Nd:YAG lasers was transferred to the target by a fused silica fiber and focused by a probe head.(b) Principle schematic diagram of a metal-assisted FO-LIBS system for measurement of gases’ molecular abundance ratios.The plasma particle (PP) produced by the laser ablation of aluminum plate (AP)collision with gases’ molecule (GM) leads to the gas molecule excitation.

    Figure 2.(a)The compared spectra of gases with the same composition between directly focusing on the gas and the metal-assisted method with FO-LIBS.(b)The compared spectra between the different gases’molecular abundance ratios.The inset shows a detailed amplification of line spectral C 833 nm.

    Figure 3.The background spectral line intensity assessment using the method of spectroscopy measurement of subtraction in a pure CH4 and N2 surrounding.The backgrounds of the C, N, O atom spectral lines were marked corresponding.

    Figure 4.The quantitative relationship between spectral line intensity and the molecular abundance ratio as the molecular abundance ratio of CH4/N2 (a) ranges from 0 to 30% and (b) ranges from 0 to 6%.(c) The intensity of the C 833 nm line spectral as the fractional atomic abundance of CH4 ranges from 0 to 0.2 for laser ablation at three energy levels.(d)The intensity of the N 818 line spectral as the fractional atomic abundance of N2 ranges from 0.7 to 1.0 during laser ablation at three energy levels.

    Figure 5.The influence of incident laser energy on the quantitative calibration of molecular abundance ratio.

    Figure 6.The quantitative concentration measurements of the CH4/N2 mixture over large range of molecular abundance ratios.

    3.4.Quantitative concentration measurements at wide abundance ratio range

    Figure 6 gives the quantitative concentration measurement results when the molecular abundance ratios of CH4/N2range from 0 to 14 and the experimental data were fitted with a linear function.It is clear that the values for all experiment points ofIC833/N818are larger than the fitted data as the molecular abundance ratio of CH4/N2is greater than 5.Additionally,the intensity ofIC833/N818decreases with the increase of the incident laser energy at the same molecular absorption ratio,which is consistent with the range of the molecular abundance ratio from 0% to 22.7%.Several factors combined lead the quantitative curve to deviate from linear function, such as the selfabsorption and collision avalanche of gas molecules.Even though the experiment points of the intensity ofIC833/N818deviate from a linear relation at a large CH4, the linear calibration analysis is relatively effective for the molecular abundance ratio of CH4/N2over wide range as the coefficient of the calibration curveR2of all fitted curves is larger than 0.99.Regardless, the proposed metal-assisted FO-LIBS is an efficient method for the measurement of gases’ molecular abundance ratios over a wide abundance ratio range.

    Figure 7.Comparison of molecular abundance ratio measurement using aluminum and copper as metal substrates.

    3.5.The influence of different metal material substrates

    To further study the effects of different metal substrates on the analysis of the component abundance ratios in a mixing gas with metal-assisted FO-LIBS, quantitative analyses of molecular abundance ratio of CH4with different metal substrates were carried out.The quantitative curve is shown in figure 7.Obviously, when Al is used as the metal substrate, the quantitative curve between the spectral line intensity and molecular immunity ratio is more linear.Comparing the two cases where Cu and Al are used as metal substrates, the total pressure of mixing gas in the two cases is the same as that of the molecular abundance ratio,and the laser energy is enough to induce metal substrate ablation; thus, the primary reason that causes the difference between the quantitative curve in the two cases should be the difference in matrix elements of the substrate.The boiling points of Cu and Al are 3200 K and 2792 K.The ablation threshold of the target was positively correlated with the boiling point of the target and the ablation amount of the target was negatively correlated with the boiling point[34,40].As a result,the lower the boiling point of the metal substrate,the higher the temperature and electron concentration of the plasma produced by laser ablation under the same conditions,which may have caused the more-intense particle collisions in the plasma.Eventually, the low boiling point of metal subtracts results in a more linear quantitative curve of the quantitative analyses of the gases.It is favorable that the boiling point of the metal substrate is lower in metalassisted FO-LIBS for gas detection.

    4.Conclusions

    To perform gas detection with the FO-LIBS system, a metalassisted laser-induced breakdown spectroscopy method was proposed for the measurement of gases’molecular abundance ratios as the laser energy is delivered by optical fiber.The quantitative relationship between the spectral line intensity and gases’molecular abundance ratio was studied over a large range.The influences of ablation laser energy and the substrate material on the quantitative calibration were also analyzed, which demonstrate that producing more plasma particles at higher plasma temperatures and electron densities is beneficial for undertaking the linear quantitative calibration of gas measurement.The proposed metal-assisted method makes the measurement of gases’molecular abundance ratios possible with the FO-LIBS system, which is a practical strategy for gas detection when the environment is too hazardous or the access space is restricted.

    Acknowledgments

    This work was supported by the National Key R&D Program of China (No.2017YFC1200400), the Development Fund of Institute of Fluid Physics, China Academy of Engineering Physics (No.SFZ20150302).

    99热全是精品| 嫩草影院新地址| 亚洲国产欧美在线一区| 99久久精品热视频| av国产免费在线观看| 欧美又色又爽又黄视频| 欧美精品一区二区大全| 一级av片app| 欧美高清成人免费视频www| 在线国产一区二区在线| 女同久久另类99精品国产91| eeuss影院久久| 九色成人免费人妻av| 国产高清有码在线观看视频| 麻豆av噜噜一区二区三区| 国产色爽女视频免费观看| 中国美白少妇内射xxxbb| 国产亚洲精品久久久久久毛片| 波多野结衣高清无吗| 国产精品久久久久久精品电影| 精品午夜福利在线看| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 丝袜喷水一区| 一级黄色大片毛片| 成人国产麻豆网| 国产av不卡久久| 久久韩国三级中文字幕| 嘟嘟电影网在线观看| av在线老鸭窝| 波野结衣二区三区在线| 亚洲欧美日韩无卡精品| 乱人视频在线观看| 午夜激情欧美在线| 美女被艹到高潮喷水动态| 三级经典国产精品| 日本黄大片高清| 99久久久亚洲精品蜜臀av| 2022亚洲国产成人精品| 成人二区视频| 日韩高清综合在线| 国产精品美女特级片免费视频播放器| 久久人妻av系列| 国产成人一区二区在线| 自拍偷自拍亚洲精品老妇| 国产一区亚洲一区在线观看| 好男人视频免费观看在线| 少妇熟女欧美另类| 天堂网av新在线| 老熟妇乱子伦视频在线观看| 亚洲成人中文字幕在线播放| 校园人妻丝袜中文字幕| 亚洲av不卡在线观看| 欧美日本视频| 久久人人精品亚洲av| 国产精品三级大全| 最近2019中文字幕mv第一页| 在线国产一区二区在线| 男女啪啪激烈高潮av片| 色哟哟·www| 性色avwww在线观看| 亚洲精品乱码久久久v下载方式| 国产伦一二天堂av在线观看| 欧美日韩在线观看h| 看非洲黑人一级黄片| 欧美激情久久久久久爽电影| 亚洲精品日韩av片在线观看| 国产av麻豆久久久久久久| 久久草成人影院| 麻豆av噜噜一区二区三区| 欧美色视频一区免费| 美女黄网站色视频| 18禁黄网站禁片免费观看直播| 麻豆久久精品国产亚洲av| 变态另类成人亚洲欧美熟女| 看免费成人av毛片| 国产女主播在线喷水免费视频网站 | 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 插逼视频在线观看| 在线观看一区二区三区| 精品免费久久久久久久清纯| 欧美另类亚洲清纯唯美| 人人妻人人澡人人爽人人夜夜 | 中文资源天堂在线| 深夜a级毛片| 狂野欧美激情性xxxx在线观看| 99在线视频只有这里精品首页| 中文字幕人妻熟人妻熟丝袜美| 欧美不卡视频在线免费观看| 搡女人真爽免费视频火全软件| 两个人视频免费观看高清| 51国产日韩欧美| 男女那种视频在线观看| 边亲边吃奶的免费视频| 国国产精品蜜臀av免费| 亚洲在线观看片| 18+在线观看网站| 亚洲中文字幕日韩| 国产成人91sexporn| 久久午夜福利片| 99热这里只有精品一区| 中文资源天堂在线| 大香蕉久久网| 99久久精品热视频| 国产精品一区www在线观看| 99热这里只有是精品50| ponron亚洲| 精品午夜福利在线看| 久久九九热精品免费| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 婷婷色综合大香蕉| 精品久久久久久久末码| 真实男女啪啪啪动态图| 99久久久亚洲精品蜜臀av| 国产成人福利小说| 亚洲成a人片在线一区二区| 久久久久久久久中文| 国产亚洲5aaaaa淫片| 青青草视频在线视频观看| 麻豆av噜噜一区二区三区| 亚洲性久久影院| 只有这里有精品99| 亚洲欧美日韩卡通动漫| 日韩欧美国产在线观看| 日韩欧美 国产精品| 久久人人爽人人爽人人片va| 最近的中文字幕免费完整| 成人欧美大片| 久久99蜜桃精品久久| 久久综合国产亚洲精品| 欧美高清性xxxxhd video| 看十八女毛片水多多多| 免费电影在线观看免费观看| 国产片特级美女逼逼视频| 国产欧美日韩精品一区二区| 欧美日本亚洲视频在线播放| 中国国产av一级| 日本撒尿小便嘘嘘汇集6| 国产成人午夜福利电影在线观看| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 欧美色欧美亚洲另类二区| 国产成人一区二区在线| 久久久久久久久久黄片| 国产私拍福利视频在线观看| 天堂影院成人在线观看| 中文字幕精品亚洲无线码一区| 女人被狂操c到高潮| 亚洲精品456在线播放app| 亚洲精品456在线播放app| 观看免费一级毛片| 久久久久久久久中文| 日韩国内少妇激情av| 欧美日韩一区二区视频在线观看视频在线 | 亚洲va在线va天堂va国产| 国产精品一区www在线观看| 夜夜夜夜夜久久久久| 免费看日本二区| 久久精品国产自在天天线| 久久99热这里只有精品18| 一个人免费在线观看电影| 精品人妻熟女av久视频| 夜夜夜夜夜久久久久| 午夜视频国产福利| 色哟哟哟哟哟哟| 91午夜精品亚洲一区二区三区| 26uuu在线亚洲综合色| 亚洲国产欧洲综合997久久,| 一级毛片久久久久久久久女| 日本欧美国产在线视频| 青春草国产在线视频 | 在线免费十八禁| 欧美色欧美亚洲另类二区| 亚洲第一电影网av| 在线国产一区二区在线| 亚洲国产精品合色在线| 两个人的视频大全免费| 亚洲自拍偷在线| 变态另类丝袜制服| 一级毛片久久久久久久久女| 给我免费播放毛片高清在线观看| 性欧美人与动物交配| 可以在线观看毛片的网站| 美女脱内裤让男人舔精品视频 | 欧美成人免费av一区二区三区| 国产高清不卡午夜福利| 综合色av麻豆| 国产久久久一区二区三区| 不卡一级毛片| 麻豆国产97在线/欧美| 99热这里只有是精品在线观看| 亚洲最大成人中文| 亚洲内射少妇av| 成年免费大片在线观看| 一级黄色大片毛片| 超碰av人人做人人爽久久| 九九爱精品视频在线观看| 精品一区二区三区人妻视频| 九九热线精品视视频播放| 欧美高清性xxxxhd video| 1000部很黄的大片| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| 午夜福利高清视频| 美女内射精品一级片tv| 国产精品人妻久久久久久| 黄片wwwwww| 男插女下体视频免费在线播放| 国产成人午夜福利电影在线观看| 3wmmmm亚洲av在线观看| 精品人妻熟女av久视频| 99视频精品全部免费 在线| 色综合站精品国产| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 97在线视频观看| 国产精品一二三区在线看| 网址你懂的国产日韩在线| 99热只有精品国产| 日本欧美国产在线视频| 国产毛片a区久久久久| 亚洲自偷自拍三级| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 午夜免费激情av| 亚洲av中文字字幕乱码综合| 日韩制服骚丝袜av| 乱系列少妇在线播放| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 国产成人freesex在线| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 亚洲经典国产精华液单| 欧美性猛交╳xxx乱大交人| 别揉我奶头 嗯啊视频| 日韩欧美 国产精品| 精品免费久久久久久久清纯| 欧美变态另类bdsm刘玥| 国产私拍福利视频在线观看| 特大巨黑吊av在线直播| 22中文网久久字幕| 午夜老司机福利剧场| 18+在线观看网站| 国产蜜桃级精品一区二区三区| 三级毛片av免费| 网址你懂的国产日韩在线| 99在线人妻在线中文字幕| 国产一级毛片在线| 国内少妇人妻偷人精品xxx网站| 69人妻影院| 精品久久久噜噜| 亚洲精品日韩av片在线观看| 久99久视频精品免费| 免费无遮挡裸体视频| 色哟哟·www| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 国产亚洲精品久久久com| 国产午夜精品论理片| 波多野结衣高清作品| 久久精品国产99精品国产亚洲性色| 国内精品美女久久久久久| 两个人视频免费观看高清| 免费观看精品视频网站| 久久久久久久久大av| 亚洲成a人片在线一区二区| 日本黄色视频三级网站网址| 亚洲七黄色美女视频| 午夜精品一区二区三区免费看| 91av网一区二区| av女优亚洲男人天堂| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 亚洲av不卡在线观看| 午夜福利高清视频| 99久国产av精品国产电影| 内射极品少妇av片p| 中文精品一卡2卡3卡4更新| 老司机福利观看| ponron亚洲| 一区二区三区免费毛片| 伦精品一区二区三区| 国产av不卡久久| 蜜桃久久精品国产亚洲av| 日本五十路高清| 18禁裸乳无遮挡免费网站照片| 亚洲,欧美,日韩| 色哟哟·www| 精品人妻一区二区三区麻豆| 久久精品夜夜夜夜夜久久蜜豆| 99热全是精品| 久久久久久九九精品二区国产| 久久韩国三级中文字幕| 久久精品国产亚洲av涩爱 | 精品少妇黑人巨大在线播放 | 国产成人aa在线观看| 国产精品一区二区在线观看99 | 久久久午夜欧美精品| videossex国产| 特级一级黄色大片| 九九热线精品视视频播放| 精品久久久久久久久久免费视频| 99久久中文字幕三级久久日本| 又爽又黄无遮挡网站| 亚洲成人久久性| 嫩草影院入口| 亚洲最大成人中文| 青青草视频在线视频观看| 日韩一区二区视频免费看| 国产高清激情床上av| 99热这里只有是精品在线观看| 不卡一级毛片| 禁无遮挡网站| 伦理电影大哥的女人| 我要看日韩黄色一级片| 中文字幕av在线有码专区| 色吧在线观看| 高清在线视频一区二区三区 | 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 日本黄大片高清| 亚洲精品影视一区二区三区av| 亚洲无线观看免费| 久久精品综合一区二区三区| 亚洲av成人精品一区久久| 国产精品av视频在线免费观看| 精品无人区乱码1区二区| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 岛国在线免费视频观看| 国产精品乱码一区二三区的特点| 18禁在线播放成人免费| 亚洲av中文字字幕乱码综合| 国产成人福利小说| 欧美最黄视频在线播放免费| 亚洲一级一片aⅴ在线观看| 人妻久久中文字幕网| 国产精品无大码| 三级毛片av免费| 麻豆久久精品国产亚洲av| 日韩高清综合在线| 欧美日韩综合久久久久久| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 97在线视频观看| 中文字幕久久专区| 亚州av有码| 久久精品国产鲁丝片午夜精品| 少妇被粗大猛烈的视频| 一级av片app| 久久99精品国语久久久| 秋霞在线观看毛片| 日日啪夜夜撸| a级毛色黄片| 22中文网久久字幕| 成人性生交大片免费视频hd| 日本黄大片高清| 亚洲四区av| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 色哟哟哟哟哟哟| 麻豆国产av国片精品| 国产高潮美女av| 麻豆国产av国片精品| 精品一区二区三区人妻视频| 一进一出抽搐动态| 国产亚洲5aaaaa淫片| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 蜜桃亚洲精品一区二区三区| 久久精品久久久久久噜噜老黄 | 国语自产精品视频在线第100页| 日韩av在线大香蕉| 欧美区成人在线视频| 午夜福利在线观看吧| 婷婷精品国产亚洲av| 美女内射精品一级片tv| 欧美三级亚洲精品| 乱系列少妇在线播放| 婷婷亚洲欧美| 国产精品久久视频播放| 日本成人三级电影网站| 美女cb高潮喷水在线观看| 你懂的网址亚洲精品在线观看 | 级片在线观看| 国产男人的电影天堂91| 99久久精品一区二区三区| 精品久久国产蜜桃| 国产精品不卡视频一区二区| 搡老妇女老女人老熟妇| 伦理电影大哥的女人| 久久午夜福利片| 九九在线视频观看精品| 男人和女人高潮做爰伦理| 国产极品精品免费视频能看的| 久久午夜亚洲精品久久| 三级国产精品欧美在线观看| 在线免费十八禁| 22中文网久久字幕| 桃色一区二区三区在线观看| 精品人妻视频免费看| 欧美3d第一页| 在线观看一区二区三区| 亚洲人成网站在线播| av在线观看视频网站免费| 欧美色视频一区免费| 老师上课跳d突然被开到最大视频| 高清日韩中文字幕在线| 国产色爽女视频免费观看| 成人三级黄色视频| 成人欧美大片| 日本一本二区三区精品| 欧美不卡视频在线免费观看| 久久亚洲精品不卡| 精品无人区乱码1区二区| 97超碰精品成人国产| 国产三级中文精品| 国产黄色小视频在线观看| 成年av动漫网址| 边亲边吃奶的免费视频| 综合色av麻豆| 夜夜夜夜夜久久久久| 成年版毛片免费区| 成年女人永久免费观看视频| 国产日本99.免费观看| 看黄色毛片网站| 国产三级在线视频| 综合色丁香网| 成人综合一区亚洲| 亚洲精品自拍成人| 亚洲无线在线观看| 久久人人爽人人片av| 身体一侧抽搐| 久久久久九九精品影院| 欧美激情国产日韩精品一区| 最近中文字幕高清免费大全6| 成人美女网站在线观看视频| 亚洲自偷自拍三级| 日本在线视频免费播放| 日韩欧美一区二区三区在线观看| 国产一级毛片在线| 亚洲精品亚洲一区二区| 国产精品久久视频播放| 91精品国产九色| 亚洲性久久影院| 成年女人永久免费观看视频| 天堂av国产一区二区熟女人妻| 日产精品乱码卡一卡2卡三| 精品无人区乱码1区二区| 欧美一级a爱片免费观看看| 一区福利在线观看| 亚洲图色成人| 久久久午夜欧美精品| 国产人妻一区二区三区在| 插逼视频在线观看| 国产淫片久久久久久久久| 欧美最新免费一区二区三区| 性色avwww在线观看| 搞女人的毛片| 久久6这里有精品| 麻豆一二三区av精品| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 亚洲va在线va天堂va国产| 看十八女毛片水多多多| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 男人和女人高潮做爰伦理| 99视频精品全部免费 在线| 能在线免费观看的黄片| 天天躁夜夜躁狠狠久久av| 亚洲人与动物交配视频| 男人舔女人下体高潮全视频| 精品久久久久久成人av| 国产精品人妻久久久影院| 亚洲人成网站在线播| 国产精品一区二区三区四区久久| 国产色爽女视频免费观看| 狂野欧美激情性xxxx在线观看| 在线国产一区二区在线| 日日干狠狠操夜夜爽| 国产精品.久久久| av在线天堂中文字幕| 波野结衣二区三区在线| 欧美成人一区二区免费高清观看| 日本免费一区二区三区高清不卡| 成人av在线播放网站| 亚洲精品日韩av片在线观看| 中出人妻视频一区二区| av又黄又爽大尺度在线免费看 | 51国产日韩欧美| 波多野结衣高清无吗| 国产精品久久电影中文字幕| 少妇丰满av| 亚洲精品456在线播放app| 国产精品,欧美在线| 精品人妻一区二区三区麻豆| 一级毛片电影观看 | 免费人成在线观看视频色| 一边亲一边摸免费视频| 内射极品少妇av片p| 一夜夜www| 国产熟女欧美一区二区| 毛片女人毛片| av天堂在线播放| 亚洲成a人片在线一区二区| 久久久成人免费电影| 久久精品国产鲁丝片午夜精品| 欧美日韩精品成人综合77777| 日本一本二区三区精品| 免费不卡的大黄色大毛片视频在线观看 | 国产熟女欧美一区二区| 老司机影院成人| 日本色播在线视频| 国产v大片淫在线免费观看| av又黄又爽大尺度在线免费看 | 国产成人精品一,二区 | 在线播放国产精品三级| 国产精品,欧美在线| 亚洲无线观看免费| 亚洲七黄色美女视频| 不卡一级毛片| 男插女下体视频免费在线播放| 亚洲av中文av极速乱| 日韩一区二区三区影片| 国产三级在线视频| 成人二区视频| 久久人人爽人人爽人人片va| 小蜜桃在线观看免费完整版高清| 18禁在线无遮挡免费观看视频| 中文亚洲av片在线观看爽| 麻豆乱淫一区二区| 亚洲精品日韩av片在线观看| 免费观看精品视频网站| 哪里可以看免费的av片| 九草在线视频观看| 男女视频在线观看网站免费| 大型黄色视频在线免费观看| 熟女人妻精品中文字幕| 黄片wwwwww| 亚洲av电影不卡..在线观看| 欧美精品国产亚洲| 小蜜桃在线观看免费完整版高清| 亚洲高清免费不卡视频| av.在线天堂| 成年免费大片在线观看| 欧美日韩精品成人综合77777| 啦啦啦韩国在线观看视频| 久久久久久大精品| 搞女人的毛片| 美女内射精品一级片tv| 免费av观看视频| 国模一区二区三区四区视频| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 亚洲中文字幕一区二区三区有码在线看| 国产精品一区www在线观看| 亚洲美女视频黄频| 国产色爽女视频免费观看| 一夜夜www| 岛国毛片在线播放| 三级毛片av免费| 久久国产乱子免费精品| 国产色婷婷99| 日韩强制内射视频| 精品熟女少妇av免费看| 99久久精品国产国产毛片| www日本黄色视频网| 欧美高清性xxxxhd video| 狂野欧美激情性xxxx在线观看| 国产精品不卡视频一区二区| 99久久精品国产国产毛片| 麻豆久久精品国产亚洲av| 特大巨黑吊av在线直播| 99热网站在线观看| 欧美+日韩+精品| 91在线精品国自产拍蜜月| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| www.色视频.com| 国产色爽女视频免费观看| 观看美女的网站| 国产免费男女视频| 亚洲成人久久爱视频| 国产午夜福利久久久久久| 久久精品国产99精品国产亚洲性色| 在线观看一区二区三区| 亚洲成av人片在线播放无| 一级二级三级毛片免费看| 天堂中文最新版在线下载 | 丝袜喷水一区| 一个人免费在线观看电影| 欧美日韩综合久久久久久| 嘟嘟电影网在线观看| 国产成人91sexporn| 精品少妇黑人巨大在线播放 | 国产成人a∨麻豆精品| 久久鲁丝午夜福利片| 夜夜爽天天搞| 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| 美女大奶头视频| 天堂av国产一区二区熟女人妻| 国产精品精品国产色婷婷| 1000部很黄的大片| 国产中年淑女户外野战色| 97超视频在线观看视频|