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

    Characteristics of laser-induced breakdown spectroscopy of liquid slag

    2024-03-19 02:37:02ChangyanDONG董長言HongxiaYU于洪霞LanxiangSUN孫蘭香YangLI李洋XiuyeLIU劉修業(yè)PingZHOU周平andShaowenHUANG黃少文
    Plasma Science and Technology 2024年2期
    關(guān)鍵詞:周平李洋

    Changyan DONG (董長言),Hongxia YU (于洪霞),Lanxiang SUN (孫蘭香),Yang LI (李洋),Xiuye LIU (劉修業(yè)),Ping ZHOU (周平) and Shaowen HUANG (黃少文)

    1 State Key Laboratory of Robotics,Shenyang Institute of Automation,Chinese Academy of Sciences,Shenyang 110016,People’s Republic of China

    2 Key Laboratory of Networked Control Systems,Chinese Academy of Sciences,Shenyang 110016,People’s Republic of China

    3 Institutes for Robotics and Intelligent Manufacturing,Chinese Academy of Sciences,Shenyang 110169,People’s Republic of China

    4 Shenyang University of Technology,Shenyang 110870,People’s Republic of China

    5 Shenyang Ligong University,Shenyang 110142,People’s Republic of China

    6 Shandong Iron and Steel Company Ltd.,Jinan 271105,People’s Republic of China

    Abstract Rapid online analysis of liquid slag is essential for optimizing the quality and energy efficiency of steel production.To investigate the key factors that affect the online measurement of refined slag using laser-induced breakdown spectroscopy (LIBS),this study examined the effects of slag composition and temperature on the intensity and stability of the LIBS spectra.The experimental temperature was controlled at three levels: 1350 °C,1400 °C,and 1450 °C.The results showed that slag composition and temperature significantly affected the intensity and stability of the LIBS spectra.Increasing the Fe content and temperature in the slag reduces its viscosity,resulting in an enhanced intensity and stability of the LIBS spectra.Additionally,42 refined slag samples were quantitatively analyzed for Fe,Si,Ca,Mg,Al,and Mn at 1350 °C,1400 °C,and 1450 °C.The normalized full spectrum combined with partial least squares (PLS) quantification modeling was used,using the Ca II 317.91 nm spectral line as an internal standard.The results show that using the internal standard normalization method can significantly reduce the influence of spectral fluctuations.Meanwhile,a temperature of 1450 °C has been found to yield superior results compared to both 1350 °C and 1400 °C,and it is advantageous to conduct a quantitative analysis of the slag when it is in a “water-like” state with low viscosity.

    Keywords: laser-induced breakdown spectroscopy (LIBS),slag,temperature,composition,viscosity,internal standard normalization,partial least squares (PLS)

    1.Introduction

    With continuous advancements in science and technology,improving the quality of steel has become a significant objective of steel production [1].Slag,a key by-product of the steelmaking process,is a valuable input for numerical simulations of the steelmaking process and directly influences the quality of the steel liquid and the final steel product [2-4].Real-time analysis of slag composition is essential for determining the end of smelting,ensuring the quality of smelting products,and reducing energy consumption [5].The alkalinity,viscosity,melting point,and other characteristics of slag are determined by its composition [6,7].The main components of slag are SiO2,CaO,Al2O3,MgO,FeO,MnO,P2O5,etc[8-10].Different component concentrations have different effects on the slag.Real-time analysis of slag components during industrial smelting processes can improve the control of subsequent process flows [11,12].

    Traditional quantitative analysis methods,such as chemical analysis [13],direct reading spectroscopy [14],and X-ray fluorescence spectroscopy (XRF) [15],are known for their high accuracy in component detection.However,their complex pre-processing procedures and long analysis time render them unsuitable for real-time online analysis.In contrast,laser-induced breakdown spectroscopy (LIBS)offers the benefits of not requiring sample preparation,fast analysis times,and online detection [16,17].Thus,LIBS is a promising technique for online analysis and detection of high-temperature liquid substances.Nollet alstudied the effects of laser radiation illuminance,propagation direction,and wavelength on the stability of LIBS detection [18].Hussainet alinvestigated the use of LIBS for rapid material analysis,analyzed elements such as Ca,Mg,and Fe in iron ore slag,and compared the results to those obtained using inductively coupled plasma (ICP) spectroscopy.The results were consistent,and the detection time was significantly reduced [19].Palagaset alstudied the parameters that impact slag analysis,including laser pulse energy,accumulated laser pulses,laser beam focusing conditions,and detector pass conditions [20].Sturmet alemployed a LIBS system to perform online analysis of solid-liquid mixed slag in a steel-making plant with a surface temperature range of(600-1400) °C.The analyzed components included FeO,SiO2,CaO,MgO,Al2O3,and MnO.The analysis took only two minutes,and the root mean square errors (RMSE) of Fe,Si,Ca,Mg,Al,and Mn in the liquid slag were 2.1%,1.4%,1.68%,0.66%,0.35%,and 0.39%,respectively [21].Maet alused slag as the experimental material and studied the effects of temperature on the plasma signal intensity,stability,and acquisition delay time in the temperature range of (50-1200) °C.The spectral intensity of the sample heated at 1000 °C is significantly higher than that of the sample heated at 200 °C [22].

    However,the temperature of slag in steel-making plants can reach 1400 °C or higher under actual field conditions,and the effect of temperature fluctuations at approximately 1400 °C on the LIBS signals has not yet been reported.In response to the need for online analysis of refined slag,experiments were conducted at three different temperatures(1350 °C,1400 °C,and 1450 °C) on slags with different component concentrations to investigate the effects of the component and temperature on spectral intensity and stability.Based on these findings,the sources of measurement error in the liquid slag were analyzed,laying the foundation for the subsequent detection of the components in the liquid slag.

    2.Experiment

    2.1.Measurement setup and parameters

    The LIBS liquid slag experimental device employed the SIALIBS analyzers developed by the Shenyang Institute of Automation Chinese Academy of Sciences.The principal diagram and photographs are shown in figure 1.The plasma excitation light source was a Nd: YAG (Neodymium-doped Yttrium Aluminium Garnet) pulse laser with an output wavelength of 1064 nm for each laser.After optimizing the experimental parameters,the pulse energy used was 100 mJ,with a pulse frequency of 5 Hz,and a spectral detection delay time of 1μs.The laser was passed through a focusing lens with a focal length of 547 mm and collected on the surface of the sample to form plasma.The wavelength range of the spectrometer was 243-361 nm with a resolution of approximately 0.1 nm.The heating device of this experiment was box type electric furnace,the heating element was silicon molybdenum rod and the temperature was adjusted from room temperature to 1500 °C.

    2.2.Slag samples and data collection

    Figure 1.Schematic diagram (a) and photograph (b) of LIBS setup for liquid slag measurement.

    Table 1.Summary of dataset.

    The slag samples used in the experiment were provided by a steel plant,and included 42 refined slags.Table 1 lists the mass fraction ranges of each component,and the label data was obtained through X-ray fluorescence analysis.During the sample preparation process,the block slag samples were crushed and ground into a fine powder.Subsequently,the sample powder was sieved using a screen with a pore size of 0.28 mm.The sieved powder was then transferred into a crucible and subjected to heating in a drying oven at 105 °C for approximately 40 min.The dried sample was then placed in a furnace,its position was adjusted,and the furnace lid was closed to begin the heating.To investigate the effects of different temperatures on the viscosity,spectral intensity,and spectral stability of the liquid slag,slag samples were heated in stages.The experimental temperature was conducted at 1350 °C,1400 °C,and 1450 °C.To ensure that the spectral quality is not affected by the deep ablation pits formed due to prolonged laser irradiation at the same location,five different positions on the sample surface were selected during the acquisition of spectral data.Each position was subjected to 400 consecutive laser shots,resulting in a total of 2000 spectra for subsequent data analysis.

    3.Result and discussion

    3.1.Effect of temperature on the viscosity of liquid slag

    In addition to its composition,the viscosity and temperature changes of slags are the main factors that cause inconsistencies in its spectrum.To detect the viscosity status of slag samples at different temperatures,this experiment used a quartz rod with calibrated scales,which was in contact with the surface of the steel slag sample through a fine rope.The experiment was repeated several times,and the depth and speed of the quartz rod sinking were recorded during each trial.If the viscosity of the sample is low,the quartz rod will sink deeper and its sinking speed will be faster over the same period of time.Conversely,if the viscosity of the sample is high,the sinking depth of the quartz rod will be shallower and the sinking speed will be slower over the same period of time.The viscosity of the slag is divided into three states in this experiment.The first state is a high-viscosity “asphaltlike” state,in which the quartz rod basically does not sink when in contact with the sample.The second state is a medium-viscosity “paste-like” state,in which the quartz rod sinks slowly when in contact with the sample.The third state is a low-viscosity “water-like” state,in which the quartz rod sinks quickly when in contact with the sample.Additionally,to investigate the impact of viscosity on the quantitative analysis,this experiment recorded the sample states at three temperatures (1350 °C,1400 °C,and 1450 °C) for 42 slag samples [23,24].Figure 2 shows that at 1350 °C,88.10% of the samples were “asphalt-like”,14% were “paste-like”,and 4.76% were “water-like”,resulting in a generally high viscosity for the samples at this temperature.At 1400 °C,38.10% of the samples were “asphalt-like”,42.86% were“paste-like”,and 19.04% were “water-like”,resulting in a large variation in viscosity among the samples,which may affect the precision of the quantitative analysis.At 1450 °C,83.33% of the samples were “water-like”,14.29% were“paste-like”,and 2.38% were “asphalt-like”,resulting in a generally low viscosity and better fluidity for the slag samples at this temperature.To study the relationship between temperature and slag viscosity,four slag samples with varying component contents were selected from the 42 samples for experimental analysis,as shown in table 2.

    Figure 2.Percentage of slag state at 1350 °C,1400 °C and 1450 °C.

    The experiment revealed that sample A was in a highviscosity “asphalt-like” state at 1350 °C and 1400 °C,and the viscosity decreased to “water-like” when the temperature increased to 1450 °C.Sample B was in a high-viscosity“asphalt-like” state at 1350 °C,and the viscosity decreased to “paste-like” state when the temperature increased to 1400 °C and further decreased to “water-like” when the temperature increased to 1450 °C.Sample C was in a medium-viscosity “paste-like” state at 1350 °C,and had relatively low viscosity,approaching “water-like”,at 1400 °C and 1450 °C.Sample D was in a high-viscosity “asphalt-like”state at 1350 °C,and was close to “water-like” at 1400 °C and 1450 °C.

    According to table 2,the Fe contents in slag samples A and B are relatively low,whereas those in slag samples C and D are relatively high.The free O2-ions in FeO can break down the silicate network structure SiO4-,lowering the viscosity of the slag.In contrast,SiO2can form more complex structures with the silicate network structure SiO4-,increasing the viscosity of the slag [25,26],Furthermore,as the Fe content in the slag increases,the viscosity of the slag decreases more rapidly with increasing temperature.Theseexperimental results indicate that temperature and composition are the primary factors influencing the viscosity of the slag,with the Fe content having a particularly significant impact.

    Table 2.Chemical composition concentration of slag samples.

    3.2.Effect of temperature on the spectrum

    To reduce large fluctuations in the slag spectrum,certain spectra with extreme intensity values were removed during preprocessing.This was done by using the peak intensity at Ca II 317.91 nm as a reference and removing spectra with intensities greater than 16000 or less than 100.Figure 3 shows the combination plots of the peak intensity at Ca II 317.91 nm for the four slag samples under different temperature conditions.Thex-axis represents the temperature and they-axis represents the spectral intensity.Each combination plot was composed of a box plot,histogram,and standard normal curve,representing the combination plots under three different temperature conditions from left to right: 1350 °C,1400 °C and 1450 °C.The yellow dashed line inside the box plot represents the median intensity,and the yellow pentagon represents the mean intensity.

    Figure 3 and table 3 show that the spectral intensities of slag samples C and D were significantly higher than those of samples A and B,and there were fewer outlier values above the box plots.The spectral stability of samples C and D was also better than that of samples A and B,particularly evident during the low-temperature phase.The effects of the temperature on the four slag samples were also different.For samples A and B,the intensity increased significantly when the temperature reached 1450 °C,and the number of abnormal points decreased compared to those at 1350 °C and 1400 °C,leading to improved stability.The effect of temperature on slag sample C was insignificant,whereas for slag sample D,the spectral line intensity increased significantly at 1400 °C,and the RSD value decreased significantly.Above 1400 °C,the change was relatively smooth.

    Figure 3.Ca II peak intensity at 317.91 nm for the four slag samples of (a) A #,(b) B #,(c) C #,(d) D #.

    Table 3.RSD of Ca II peak intensity at 317.91 nm for different slag samples (%).

    By analyzing the relationship between temperature and viscosity of liquid slag in section 3.1,in conjunction with data presented in figure 3,we can observe a strong correlation between the spectral intensity of liquid slag and its viscosity.When the slag is in a “paste-like” or “watery-like”state with a low viscosity,the spectral intensity is stronger and more stable.In contrast,when the slag is in an “asphaltlike” state with a higher viscosity,the spectral intensity is weaker,and there are many outlier points in the box plot.Furthermore,during the experiment,it was observed that when the slag was in an “asphalt-like” state,the spectral line was initially strong but gradually weakened each time the sample moved.This was likely owing to the high viscosity of the slag,which caused the laser to quickly generate unhealable “pits” on the surface of the slag,resulting in large spectral fluctuations and a decrease in intensity.In contrast,when the slag is in a “paste-like” or “water-like” state with low viscosity,the “pits” heal quickly,resulting in more stable spectra.Because spectral intensity of “asphalt-like”slag is only strong initially after it is moved,these points are considered outliers when creating the box plot.

    Internal standard methods are often used to reduce the impact of conditional fluctuations during quantitative analysis.To investigate whether the effect of the temperature can be compensated for using this method,we used the spectral line intensity of Ca II at 317.91 nm as an internal standard and investigated the trend of the spectral line intensity ratio.Figure 4 shows the combination plot of the ratio of the Fe and Ca spectral line intensities for the different slag samples.Thex-axis is the same as that in figure 3,and they-axis represents the ratio of the Fe II 259.90 nm and Ca II 317.91 nm spectral peak intensities.The corresponding RSD values are presented in table 4.

    It can be seen from table 4 that when using the internal standard normalization method,the RSD of the four slag samples was significantly reduced.As shown in figure 4,the data of slags samples A and B are more concentrated at the mean,and more obedient to a normal distribution at 1350 °C and 1400 °C.Compared to the results without pretreatment,the RSD was reduced by more than 50%.At a temperature of 1450 °C,the RSD was also decreasing from above 40% to 10%-20%.At the same time,the RSD of slag sample C was decreased from above 30% to within 10% at 1350 °C,1400 °C and 1450 °C.For slag sample D,a change in the RSD can also be observed in the asphalt state at 1350 °C,and it is more likely to follow a normal distribution.Furthermore,the RSD was reduced to less than 4% at 1400 °C and above.

    Figure 4.Peak intensity ratio of Fe II 259.90 nm to Ca II 317.9 nm in different slag samples.

    Table 4.Variation of the RSD of Fe II 259.90 nm peak intensity normalized by internal standard of Ca II 317.91 nm in different slag samples before and after normalization.

    The above results show that internal normalization can effectively reduce the impact of spectral fluctuations,particularly when viscosity is low.This method can also compensate for temperature fluctuations.However,when the viscosity was high,although the RSD value decreased significantly,it remained high and could not mitigate the effects of the temperature fluctuations.The key factor influencing viscosity is the Fe content.When the Fe content was relatively high,the flow ability of the liquid slag improved,and internal normalization effectively reduced the effect of temperature on the analysis results.

    3.3.Quantitative analysis of slag composition

    To investigate the effects of viscosity and temperature on the quantitative analysis of liquid slag,42 slag samples were used at temperatures of 1350 °C,1400 °C and 1450 °C,and 2000 spectra were collected for each sample.The upper and lower limits of the box plot of the Ca II 317.91 nm peak intensity were used as boundary lines,and the spectra between the upper and lower limits were selected for mean normalization.The spectral data were obtained for subsequent analyses by averaging the spectra of each sample.To minimize the effect of experimental condition fluctuations on data stability,we employed internal normalization to normalize the spectra using the Ca II 317.91 nm peak intensity as the internal standard.Of the 42 slag samples,34 were selected as the training set,and the remaining eight were used as the test set.To ensure that the training and test sets were representative,we divided them into different concentration intervals,and randomly selected one or more samples from each interval to ensure that the labels of the test set samples were evenly distributed.Because the partial least squares (PLS) algorithm is effective for handling issues of small sample size,high-dimensional variables,and multicollinearity among variables [27-29],this experiment used Ca internal normalization in combination with the PLS method to quantitatively analyze the TFe,SiO2,CaO,MgO,Al2O3,and MnO components in molten slag at three temperature levels of 1350 °C,1400 °C and 1450 °C.

    To determine the optimal number of latent variables in the model,we preprocessed the 34×2048 spectral data obtained through Ca internal normalization and used it as the input for the PLS model,with the corresponding mass fractions of the sample elements as the output.The number of latent variables was evaluated by calculating the minimum root mean square error (RMSE) of leave-one-out cross-validation for each element.

    Figure 5.Relationship between the RMSE of leave-one-out cross-validation and the number of latent variables for each element of the PLS model (1450 °C).

    As shown in figure 5,when the number of latent variables for Fe,Si,Ca,Mg,Al,and Mn in the slag at 1450 °C was 4,6,6,6,4,and 5,respectively,the corresponding RMSE was minimized.To investigate the comparability of the quantitative analysis results at 1350 °C,1400 °C and 1450 °C,we uniformly set the number of latent variables for each element at each temperature,while ensuring that the other experimental parameters and the test set are consistent.Finally,we selected the optimal number of latent variables for each element to establish the PLS calibration model.The quantitative analysis results of slag at three temperature stages are shown in figure 6,the blue pattern represents the training set and the red pattern represents the test set,each element corresponds to different geometric shapes.The evaluation metrics shown in the graph are the total coefficient of determination (R2) for the six major elements of slag in the test and training sets,the mean root mean square error of cross-validation (RMSECV),the mean root mean square error of the training set (RMSETR) and the mean root mean square error of the test set (RMSEP) for overall elements.

    According to figure 6,the average RMSEP values for elements in slag at temperatures of 1350 °C,1400 °C and 1450 °C are 2.38,2.54,and 1.54,respectively.It is evident that the overall predictive performance is optimal at a temperature of 1450 °C.Combining this with table 5 and figure 7,it can be seen that for individual element predictions,theR2values are also highest at 1450 °C.This is because the slag exhibits better fluidity at 1450 °C,with over 80% of the slag samples showing lower viscosity “waterlike”.As discussed in section 3.2,experimental observations indicate that slag spectra are stronger and more stable in a low viscosity state,especially when the slag viscosity resembles that of “water-like”,resulting in higher stability.This effect is particularly evident in the measurement of MgO content.At 1350 °C and 1400 °C,the correspondingR2values are only 0.49 and 0.55,whereas at 1450 °C,theR2value increases to a higher value of 0.89.This is mainly attributed to the fact that MgO has the highest melting point among the oxides in the slag,and at 1350 °C and 1400 °C,MgO may not have completely melted,leading to significant spectral fluctuations and affecting the accuracy of quantitative analysis.

    Based on the above,it can be concluded that both temperature and viscosity of the slag significantly impact the accuracy of quantitative analysis.Compared to temperatures of 1350 °C and 1400 °C,a temperature of 1450 °C offers greater advantages for precise quantitative analysis of the slag.Furthermore,it is highly advantageous to conduct a quantitative analysis of the slag when it is in a “water-like”state with low viscosity.

    Figure 6.Mass fraction analysis results of slag composition at 1350 °C (a),1400 °C (b),and 1450 °C (c).

    Table 5.Quantitative analysis information of each element at different temperature stages.

    Figure 7.RMESP of six elements at 1350 °C,1400 °C and 1450 °C.

    4.Conclusions

    In this study,the LIBS technique was used to investigate the effects of slag composition and temperature on spectral intensity and stability.The results indicate that both factors have a significant impact on the spectral properties of slag.The viscosity of the slag exhibits a strong correlation with its Fe content and temperature,particularly when the Fe content and temperature are elevated.This results in a reduction in viscosity,which in turn leads to an improved spectral intensity and stability.These findings indicate that slag with low Fe content has a weaker spectrum and poorer spectral stability,whereas slag with relatively high Fe content has a stronger spectrum and better spectral stability,making it more suitable for quantitative analysis.

    In addition,42 slag samples were quantitatively analyzed for TFe,SiO2,CaO,MgO,Al2O3,and MnO at temperatures of 1350 °C,1400 °C and 1450 °C.The mean RMSEP values for the six elements under the three temperature conditions are 2.38,2.54,and 1.54,respectively.The results show that the quantitative analysis of each element in the slag yielded more favorable results at a temperature of 1450 °C when compared to temperatures of 1350 °C and 1400 °C.Furthermore,it is advantageous to conduct a quantitative analysis of the slag when it is in a “water-like” state with low viscosity.

    Acknowledgments

    This work was financially supported by the National Key R&D Program Projects of China (No.2021YFB3202402)and National Natural Science Foundation of China (No.62173321).

    猜你喜歡
    周平李洋
    Time-dependent variational approach to solve multi-dimensional time-dependent Schr¨odinger equation
    拼湊
    揚子江(2023年5期)2023-09-11 12:34:39
    偏心父親生生把孝順兒子逼成了殺人犯
    婦女生活(2022年12期)2023-01-20 03:12:29
    李洋、龔有月、賴姝、程曉春作品
    大觀(2022年12期)2022-02-23 00:44:30
    李洋論
    中華詩詞(2021年1期)2021-12-31 07:51:46
    New teacher—student Relationship in Junior English Class Teaching
    魅力中國(2018年51期)2018-04-08 09:09:34
    希式男友
    特殊的美食比賽
    三月三(2015年8期)2015-08-12 16:49:35
    “奇怪”的算術(shù)題
    遺產(chǎn)作誘餌“替夫征婚”,單身女子落入圈套
    百姓生活(2012年2期)2012-04-29 00:44:03
    精品国产国语对白av| 免费在线观看日本一区| 一区二区av电影网| 免费人妻精品一区二区三区视频| 麻豆av在线久日| 两性夫妻黄色片| 女警被强在线播放| 99热国产这里只有精品6| 国产精品秋霞免费鲁丝片| av网站在线播放免费| 精品第一国产精品| 狂野欧美激情性bbbbbb| 精品乱码久久久久久99久播| 国产一区二区三区av在线| 欧美精品av麻豆av| 亚洲精品成人av观看孕妇| avwww免费| 大型av网站在线播放| 在线av久久热| av电影中文网址| 韩国精品一区二区三区| 一个人免费在线观看的高清视频 | 夜夜骑夜夜射夜夜干| 欧美+亚洲+日韩+国产| 国产欧美亚洲国产| 99久久国产精品久久久| 午夜福利视频精品| 国产野战对白在线观看| 日本a在线网址| 最新在线观看一区二区三区| 美国免费a级毛片| 国产在线一区二区三区精| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕高清在线视频| 亚洲专区中文字幕在线| 精品人妻熟女毛片av久久网站| 超色免费av| 欧美另类一区| 日韩一区二区三区影片| 国产精品二区激情视频| 69精品国产乱码久久久| 啦啦啦在线免费观看视频4| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 亚洲五月婷婷丁香| 日日爽夜夜爽网站| 99久久人妻综合| 久久久久视频综合| 国产精品熟女久久久久浪| 午夜福利视频在线观看免费| 国产成人精品无人区| 91麻豆av在线| 97人妻天天添夜夜摸| 日韩电影二区| 18在线观看网站| 久久久精品94久久精品| 不卡一级毛片| 午夜福利免费观看在线| 最新在线观看一区二区三区| 曰老女人黄片| 久久亚洲国产成人精品v| 男女边摸边吃奶| svipshipincom国产片| 一边摸一边抽搐一进一出视频| 免费少妇av软件| 成年美女黄网站色视频大全免费| www.精华液| 9热在线视频观看99| 欧美在线黄色| 汤姆久久久久久久影院中文字幕| 免费在线观看日本一区| 新久久久久国产一级毛片| 免费高清在线观看日韩| 桃花免费在线播放| 欧美日韩av久久| 亚洲国产看品久久| 一区二区三区精品91| 久久人人爽av亚洲精品天堂| 自线自在国产av| 亚洲人成77777在线视频| 国产成人一区二区三区免费视频网站| 欧美少妇被猛烈插入视频| 亚洲精品国产色婷婷电影| 熟女少妇亚洲综合色aaa.| 日韩有码中文字幕| 国产黄色免费在线视频| 一区在线观看完整版| 99久久精品国产亚洲精品| 午夜久久久在线观看| 精品第一国产精品| 国产精品麻豆人妻色哟哟久久| 欧美精品亚洲一区二区| 欧美日韩av久久| 久久人人爽人人片av| 久久久久国产精品人妻一区二区| 国产精品国产三级国产专区5o| 91麻豆精品激情在线观看国产 | 亚洲精品一卡2卡三卡4卡5卡 | 在线av久久热| 国产精品免费视频内射| 国产精品九九99| 精品国产一区二区久久| 老司机影院成人| 午夜福利乱码中文字幕| 后天国语完整版免费观看| 国产99久久九九免费精品| 午夜精品国产一区二区电影| 黄频高清免费视频| 9色porny在线观看| 久久精品国产a三级三级三级| 国产精品自产拍在线观看55亚洲 | 中文字幕人妻熟女乱码| 又黄又粗又硬又大视频| 宅男免费午夜| 日韩制服骚丝袜av| av免费在线观看网站| 国产av一区二区精品久久| av超薄肉色丝袜交足视频| 亚洲精品国产色婷婷电影| 热re99久久精品国产66热6| 国产一区二区激情短视频 | 久热这里只有精品99| 亚洲精品久久久久久婷婷小说| 国产一区二区 视频在线| 亚洲国产欧美日韩在线播放| 天堂中文最新版在线下载| 高清av免费在线| 九色亚洲精品在线播放| 19禁男女啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 亚洲欧美精品综合一区二区三区| 亚洲国产日韩一区二区| 亚洲第一欧美日韩一区二区三区 | 午夜视频精品福利| 亚洲欧美精品自产自拍| 久久久久久久久免费视频了| 国产野战对白在线观看| 国产精品偷伦视频观看了| 丁香六月欧美| 精品视频人人做人人爽| 亚洲天堂av无毛| 精品少妇内射三级| 午夜福利乱码中文字幕| 黄色片一级片一级黄色片| 欧美 亚洲 国产 日韩一| 一区福利在线观看| 欧美精品一区二区免费开放| 欧美人与性动交α欧美软件| 欧美黑人欧美精品刺激| 国产97色在线日韩免费| 在线观看www视频免费| 亚洲中文av在线| 一本—道久久a久久精品蜜桃钙片| 黄色视频在线播放观看不卡| 成人免费观看视频高清| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 久久人妻熟女aⅴ| 伦理电影免费视频| 亚洲伊人色综图| 久久99一区二区三区| av天堂久久9| 三上悠亚av全集在线观看| 国产精品av久久久久免费| 2018国产大陆天天弄谢| 欧美日韩亚洲高清精品| 国产成人精品久久二区二区91| 久久ye,这里只有精品| 18禁观看日本| 一级片免费观看大全| 天堂8中文在线网| 亚洲伊人色综图| 伊人亚洲综合成人网| 色视频在线一区二区三区| 精品人妻一区二区三区麻豆| 十八禁人妻一区二区| 少妇猛男粗大的猛烈进出视频| 天堂8中文在线网| 自线自在国产av| cao死你这个sao货| 精品高清国产在线一区| 久久毛片免费看一区二区三区| 日韩有码中文字幕| 9热在线视频观看99| 亚洲国产日韩一区二区| 男女国产视频网站| 丝袜喷水一区| av福利片在线| 在线 av 中文字幕| 人人妻人人澡人人爽人人夜夜| 欧美中文综合在线视频| 丁香六月天网| 69精品国产乱码久久久| 久久国产精品影院| 国产一区二区三区在线臀色熟女 | 午夜福利,免费看| 老司机午夜十八禁免费视频| 国产精品免费视频内射| 国产伦人伦偷精品视频| 狂野欧美激情性xxxx| 777米奇影视久久| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产欧美在线一区| 大片免费播放器 马上看| 成年av动漫网址| 成在线人永久免费视频| 国产亚洲欧美精品永久| 97人妻天天添夜夜摸| av欧美777| 人妻一区二区av| 国产精品一区二区在线观看99| bbb黄色大片| 午夜福利视频精品| 日韩欧美国产一区二区入口| 女人久久www免费人成看片| 国产成人一区二区三区免费视频网站| 我的亚洲天堂| 天天躁日日躁夜夜躁夜夜| 如日韩欧美国产精品一区二区三区| 欧美精品一区二区大全| 又大又爽又粗| 婷婷丁香在线五月| 岛国在线观看网站| av欧美777| 母亲3免费完整高清在线观看| www.熟女人妻精品国产| 久久人人爽人人片av| 免费在线观看完整版高清| 久久性视频一级片| 青青草视频在线视频观看| 另类精品久久| www.999成人在线观看| 如日韩欧美国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 飞空精品影院首页| 99re6热这里在线精品视频| 免费看十八禁软件| 国产男人的电影天堂91| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲国产毛片av蜜桃av| 久久久水蜜桃国产精品网| 免费观看a级毛片全部| 91精品伊人久久大香线蕉| 午夜影院在线不卡| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区在线不卡| 黄色 视频免费看| 黑丝袜美女国产一区| av网站在线播放免费| 久热爱精品视频在线9| 捣出白浆h1v1| 后天国语完整版免费观看| 精品国产国语对白av| 搡老岳熟女国产| a在线观看视频网站| 丰满少妇做爰视频| 精品一区二区三卡| 午夜福利在线观看吧| 国产高清国产精品国产三级| 黑人巨大精品欧美一区二区蜜桃| 欧美xxⅹ黑人| 国产精品一区二区免费欧美 | 欧美日韩福利视频一区二区| 91老司机精品| 午夜视频精品福利| 俄罗斯特黄特色一大片| 99久久精品国产亚洲精品| 无限看片的www在线观看| 一区二区日韩欧美中文字幕| 久久性视频一级片| 国产1区2区3区精品| 丁香六月天网| 天天操日日干夜夜撸| 亚洲精品国产色婷婷电影| 不卡一级毛片| 90打野战视频偷拍视频| 爱豆传媒免费全集在线观看| 九色亚洲精品在线播放| 亚洲久久久国产精品| 日本wwww免费看| 亚洲专区字幕在线| 国产欧美日韩一区二区三 | 黑人操中国人逼视频| 少妇被粗大的猛进出69影院| 两性夫妻黄色片| 老司机影院成人| 天天添夜夜摸| 亚洲av成人一区二区三| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 国产精品国产av在线观看| 亚洲成国产人片在线观看| 自线自在国产av| 如日韩欧美国产精品一区二区三区| 国产精品99久久99久久久不卡| 大陆偷拍与自拍| 看免费av毛片| 十八禁网站网址无遮挡| 免费av中文字幕在线| 国产精品麻豆人妻色哟哟久久| 国产精品秋霞免费鲁丝片| 欧美日韩黄片免| 欧美精品人与动牲交sv欧美| 精品国产国语对白av| 国产又爽黄色视频| 国产成人免费观看mmmm| 狠狠精品人妻久久久久久综合| h视频一区二区三区| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 电影成人av| 亚洲专区字幕在线| 肉色欧美久久久久久久蜜桃| 女人爽到高潮嗷嗷叫在线视频| 欧美日本中文国产一区发布| 精品一区二区三区av网在线观看 | 国产伦理片在线播放av一区| 成人18禁高潮啪啪吃奶动态图| 一区福利在线观看| 法律面前人人平等表现在哪些方面 | 久久av网站| 欧美另类一区| 亚洲第一av免费看| 国产成人欧美在线观看 | 老汉色∧v一级毛片| 亚洲一区二区三区欧美精品| 免费在线观看视频国产中文字幕亚洲 | 久久热在线av| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三区在线| 丁香六月天网| 国内毛片毛片毛片毛片毛片| 成人国语在线视频| 色老头精品视频在线观看| 建设人人有责人人尽责人人享有的| 精品乱码久久久久久99久播| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲av一区麻豆| 久久久久久久久久久久大奶| 日韩大片免费观看网站| 久久久久网色| 欧美日本中文国产一区发布| 欧美黄色淫秽网站| 国产成人a∨麻豆精品| 丝袜喷水一区| 交换朋友夫妻互换小说| av在线播放精品| 亚洲国产毛片av蜜桃av| 交换朋友夫妻互换小说| 制服人妻中文乱码| 亚洲精品国产区一区二| 精品免费久久久久久久清纯 | 久久久欧美国产精品| 日本vs欧美在线观看视频| 麻豆乱淫一区二区| 亚洲性夜色夜夜综合| 91大片在线观看| 下体分泌物呈黄色| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看www视频免费| 日韩欧美免费精品| 久久综合国产亚洲精品| 国产主播在线观看一区二区| 在线观看免费日韩欧美大片| 老司机亚洲免费影院| 国产三级黄色录像| 国产精品亚洲av一区麻豆| 亚洲精品国产av成人精品| 欧美97在线视频| 久久99一区二区三区| 一级毛片女人18水好多| 欧美精品啪啪一区二区三区 | 国产一区二区三区综合在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线 av 中文字幕| 在线精品无人区一区二区三| 国产精品偷伦视频观看了| 中国国产av一级| 十八禁高潮呻吟视频| 国产精品二区激情视频| 天天躁狠狠躁夜夜躁狠狠躁| 日韩人妻精品一区2区三区| 真人做人爱边吃奶动态| 人人妻,人人澡人人爽秒播| 叶爱在线成人免费视频播放| 婷婷丁香在线五月| 999久久久精品免费观看国产| 狠狠婷婷综合久久久久久88av| 我的亚洲天堂| h视频一区二区三区| 日韩制服骚丝袜av| 啦啦啦啦在线视频资源| 最近最新中文字幕大全免费视频| 丁香六月欧美| 黄色视频,在线免费观看| 99久久99久久久精品蜜桃| 美女脱内裤让男人舔精品视频| 国产在视频线精品| 男女免费视频国产| 亚洲五月色婷婷综合| 最新的欧美精品一区二区| 国产精品一二三区在线看| 日本av免费视频播放| 亚洲国产欧美日韩在线播放| 日韩一卡2卡3卡4卡2021年| 精品亚洲成国产av| 老司机福利观看| 一区福利在线观看| 久久热在线av| 国产精品久久久人人做人人爽| 国产区一区二久久| 建设人人有责人人尽责人人享有的| 大型av网站在线播放| 国产在视频线精品| 国产精品.久久久| 精品国产乱子伦一区二区三区 | 午夜老司机福利片| 丰满饥渴人妻一区二区三| 日韩人妻精品一区2区三区| 岛国在线观看网站| 亚洲人成77777在线视频| 欧美少妇被猛烈插入视频| 91字幕亚洲| 交换朋友夫妻互换小说| 热re99久久国产66热| 在线永久观看黄色视频| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 超色免费av| 丁香六月欧美| 一区二区日韩欧美中文字幕| 国产成人一区二区三区免费视频网站| 少妇 在线观看| 国产免费福利视频在线观看| 午夜91福利影院| 五月开心婷婷网| 天堂8中文在线网| 精品人妻在线不人妻| av又黄又爽大尺度在线免费看| 亚洲伊人色综图| 女警被强在线播放| 每晚都被弄得嗷嗷叫到高潮| 极品人妻少妇av视频| 免费少妇av软件| 亚洲精品国产色婷婷电影| 精品少妇黑人巨大在线播放| 亚洲成人手机| cao死你这个sao货| 不卡一级毛片| 女人精品久久久久毛片| 考比视频在线观看| 中文字幕av电影在线播放| 大陆偷拍与自拍| 手机成人av网站| 国产黄频视频在线观看| 丝瓜视频免费看黄片| 国产成人欧美在线观看 | 麻豆国产av国片精品| 亚洲av成人不卡在线观看播放网 | 精品乱码久久久久久99久播| 热re99久久国产66热| 亚洲欧美精品自产自拍| 建设人人有责人人尽责人人享有的| 在线观看免费午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久ye,这里只有精品| 亚洲精品粉嫩美女一区| 中文字幕精品免费在线观看视频| 又大又爽又粗| 婷婷成人精品国产| 日韩大片免费观看网站| 国产不卡av网站在线观看| 欧美国产精品va在线观看不卡| 亚洲欧美一区二区三区黑人| 国产亚洲一区二区精品| 日本vs欧美在线观看视频| 日韩精品免费视频一区二区三区| 亚洲欧美日韩高清在线视频 | 女人精品久久久久毛片| 日韩欧美一区二区三区在线观看 | 三级毛片av免费| 国产精品久久久av美女十八| 一级片'在线观看视频| 午夜福利,免费看| 99九九在线精品视频| 一级片免费观看大全| 飞空精品影院首页| 久久久久国产一级毛片高清牌| 18禁国产床啪视频网站| 国产免费一区二区三区四区乱码| a级毛片在线看网站| 蜜桃国产av成人99| 操美女的视频在线观看| 国产av一区二区精品久久| 亚洲精品久久午夜乱码| 日本a在线网址| 韩国高清视频一区二区三区| 两人在一起打扑克的视频| 午夜福利一区二区在线看| 免费在线观看日本一区| 狠狠狠狠99中文字幕| 午夜两性在线视频| 操美女的视频在线观看| 91麻豆精品激情在线观看国产 | 老鸭窝网址在线观看| 少妇粗大呻吟视频| av网站在线播放免费| 99re6热这里在线精品视频| 美国免费a级毛片| 欧美精品人与动牲交sv欧美| 侵犯人妻中文字幕一二三四区| 桃红色精品国产亚洲av| 91国产中文字幕| 成年人午夜在线观看视频| 天天操日日干夜夜撸| 欧美日韩精品网址| 色精品久久人妻99蜜桃| 亚洲精品第二区| 欧美成狂野欧美在线观看| 婷婷色av中文字幕| 嫁个100分男人电影在线观看| 色播在线永久视频| 中文字幕人妻丝袜一区二区| 人人澡人人妻人| 国产欧美日韩精品亚洲av| 国产免费视频播放在线视频| 最近中文字幕2019免费版| av电影中文网址| 老鸭窝网址在线观看| 亚洲三区欧美一区| 亚洲人成电影免费在线| 精品亚洲成国产av| 777米奇影视久久| 欧美黑人精品巨大| 精品第一国产精品| 久久久水蜜桃国产精品网| 国产成人免费无遮挡视频| 99精国产麻豆久久婷婷| 咕卡用的链子| 最近中文字幕2019免费版| 国产99久久九九免费精品| 女人久久www免费人成看片| 国产免费一区二区三区四区乱码| 大陆偷拍与自拍| 老司机深夜福利视频在线观看 | 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 一级毛片精品| 国产成人精品在线电影| 一区二区三区乱码不卡18| 国产伦人伦偷精品视频| 国产精品一二三区在线看| 成人国语在线视频| 国产欧美日韩精品亚洲av| 久久久久网色| 亚洲欧美成人综合另类久久久| 日韩欧美一区二区三区在线观看 | 在线av久久热| www.999成人在线观看| 成年人午夜在线观看视频| 日本一区二区免费在线视频| 久久毛片免费看一区二区三区| 丰满少妇做爰视频| 欧美成狂野欧美在线观看| 国产日韩欧美在线精品| 日日爽夜夜爽网站| 电影成人av| av天堂在线播放| 久久人人爽av亚洲精品天堂| 在线看a的网站| 成年av动漫网址| 久久ye,这里只有精品| 香蕉丝袜av| 亚洲精品国产区一区二| 国产精品免费大片| 天天添夜夜摸| 婷婷色av中文字幕| www.精华液| 黄色视频,在线免费观看| 妹子高潮喷水视频| 熟女少妇亚洲综合色aaa.| 国产精品一区二区在线观看99| 亚洲成人免费电影在线观看| 成年女人毛片免费观看观看9 | 老熟女久久久| 国产精品一区二区免费欧美 | 国产一级毛片在线| 亚洲欧洲日产国产| 俄罗斯特黄特色一大片| 黄色 视频免费看| 一级毛片电影观看| 中文字幕另类日韩欧美亚洲嫩草| 女人被躁到高潮嗷嗷叫费观| 男女国产视频网站| 国产伦人伦偷精品视频| 18禁观看日本| 久久国产精品人妻蜜桃| 男女免费视频国产| 黑人巨大精品欧美一区二区蜜桃| 国产精品一区二区免费欧美 | 成人av一区二区三区在线看 | 成人亚洲精品一区在线观看| 建设人人有责人人尽责人人享有的| 亚洲综合色网址| videosex国产| 妹子高潮喷水视频| √禁漫天堂资源中文www| 亚洲九九香蕉| 狠狠精品人妻久久久久久综合| av天堂久久9| 亚洲国产精品一区二区三区在线| 丝瓜视频免费看黄片| 国产成人av激情在线播放| 欧美精品人与动牲交sv欧美|