• <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
    嫩草影视91久久| av有码第一页| 亚洲无线在线观看| 又大又爽又粗| 最新在线观看一区二区三区| 亚洲av中文字字幕乱码综合| 亚洲精品在线美女| 十八禁网站免费在线| 久久久久国内视频| 亚洲天堂国产精品一区在线| 在线十欧美十亚洲十日本专区| 一进一出抽搐gif免费好疼| 国产主播在线观看一区二区| 搡老妇女老女人老熟妇| 精品午夜福利视频在线观看一区| 中文字幕av在线有码专区| 久久婷婷成人综合色麻豆| 午夜福利欧美成人| 亚洲黑人精品在线| 久久伊人香网站| 在线十欧美十亚洲十日本专区| 亚洲乱码一区二区免费版| 五月伊人婷婷丁香| 国产成年人精品一区二区| 香蕉av资源在线| 真人做人爱边吃奶动态| 看片在线看免费视频| 一级毛片女人18水好多| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品综合一区在线观看 | 欧美最黄视频在线播放免费| 99精品在免费线老司机午夜| 亚洲精品一卡2卡三卡4卡5卡| 国产成人系列免费观看| 亚洲欧美一区二区三区黑人| 成人亚洲精品av一区二区| 亚洲一区二区三区不卡视频| 亚洲乱码一区二区免费版| 国产一区二区三区在线臀色熟女| 久久久久性生活片| 少妇熟女aⅴ在线视频| 久久国产乱子伦精品免费另类| or卡值多少钱| 99热6这里只有精品| 国产1区2区3区精品| 悠悠久久av| 中文字幕人妻丝袜一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人精品一区久久| 国产精品久久久久久久电影 | 美女大奶头视频| 日本免费a在线| 欧美色欧美亚洲另类二区| 国产亚洲精品综合一区在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 午夜两性在线视频| 国产成年人精品一区二区| 毛片女人毛片| 欧美高清成人免费视频www| 好男人电影高清在线观看| 久久香蕉国产精品| 久久久久久亚洲精品国产蜜桃av| 国产精品乱码一区二三区的特点| 五月伊人婷婷丁香| 1024视频免费在线观看| 欧美一级a爱片免费观看看 | 成熟少妇高潮喷水视频| 久热爱精品视频在线9| 亚洲美女视频黄频| 成年人黄色毛片网站| 黄色丝袜av网址大全| 久热爱精品视频在线9| 久久精品国产清高在天天线| 欧美色视频一区免费| 看免费av毛片| 亚洲真实伦在线观看| 黄色视频,在线免费观看| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 国产又色又爽无遮挡免费看| 在线观看午夜福利视频| 国产激情欧美一区二区| 精品久久久久久,| 免费观看人在逋| 老司机午夜福利在线观看视频| 不卡av一区二区三区| 18禁黄网站禁片免费观看直播| 成人午夜高清在线视频| 亚洲国产中文字幕在线视频| 精品国产超薄肉色丝袜足j| 国产高清视频在线观看网站| 亚洲人与动物交配视频| 国产精品乱码一区二三区的特点| 久久天堂一区二区三区四区| 欧美日韩福利视频一区二区| 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影| 俺也久久电影网| 国内精品久久久久精免费| 我的老师免费观看完整版| 哪里可以看免费的av片| 首页视频小说图片口味搜索| 国产精品98久久久久久宅男小说| 深夜精品福利| 久久久国产成人免费| 成人精品一区二区免费| 午夜福利欧美成人| 一本一本综合久久| 一区二区三区国产精品乱码| 美女免费视频网站| 中国美女看黄片| 成人18禁高潮啪啪吃奶动态图| 老司机靠b影院| 国产亚洲欧美在线一区二区| 久久久久久亚洲精品国产蜜桃av| 又紧又爽又黄一区二区| 99精品欧美一区二区三区四区| 国产熟女午夜一区二区三区| 香蕉久久夜色| 一进一出好大好爽视频| x7x7x7水蜜桃| 手机成人av网站| 精品国产美女av久久久久小说| 欧美乱妇无乱码| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 全区人妻精品视频| 国产精品日韩av在线免费观看| 国产区一区二久久| 一区二区三区激情视频| 免费看美女性在线毛片视频| 午夜久久久久精精品| 免费在线观看影片大全网站| 欧美日韩亚洲综合一区二区三区_| 成人三级做爰电影| 亚洲人成网站在线播放欧美日韩| 90打野战视频偷拍视频| 欧洲精品卡2卡3卡4卡5卡区| 国产高清视频在线播放一区| 最近最新中文字幕大全免费视频| 最近最新中文字幕大全电影3| 欧美一区二区精品小视频在线| 欧美性猛交黑人性爽| 国产精品亚洲av一区麻豆| 婷婷亚洲欧美| 久久中文看片网| 欧美日韩一级在线毛片| 曰老女人黄片| av有码第一页| 黄色丝袜av网址大全| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩乱码在线| 国内久久婷婷六月综合欲色啪| 精品国产亚洲在线| 久久人妻福利社区极品人妻图片| 伦理电影免费视频| av中文乱码字幕在线| 久久精品国产亚洲av高清一级| 亚洲第一电影网av| 久久婷婷成人综合色麻豆| 国产精品精品国产色婷婷| 人妻久久中文字幕网| 亚洲人成网站在线播放欧美日韩| 国产精品精品国产色婷婷| 日韩欧美在线二视频| 国产一区二区三区视频了| 亚洲av电影在线进入| 精品欧美一区二区三区在线| 一区二区三区高清视频在线| 国产成人精品久久二区二区91| 中文在线观看免费www的网站 | 成人三级做爰电影| 国内少妇人妻偷人精品xxx网站 | 欧美日韩亚洲国产一区二区在线观看| 成人特级黄色片久久久久久久| 欧美日韩中文字幕国产精品一区二区三区| 99热这里只有精品一区 | 亚洲精品在线美女| 五月玫瑰六月丁香| 亚洲国产欧美一区二区综合| xxxwww97欧美| 国产亚洲精品综合一区在线观看 | 久久性视频一级片| 99热这里只有精品一区 | 两个人看的免费小视频| 欧美三级亚洲精品| 国产亚洲精品久久久久5区| 天堂影院成人在线观看| 欧美日本亚洲视频在线播放| 91国产中文字幕| 国产成人精品久久二区二区免费| 亚洲全国av大片| 亚洲人成77777在线视频| 狂野欧美白嫩少妇大欣赏| 亚洲 欧美一区二区三区| av福利片在线观看| www日本黄色视频网| 51午夜福利影视在线观看| 最新在线观看一区二区三区| 午夜久久久久精精品| 在线a可以看的网站| 国产精品av视频在线免费观看| 国产成人精品无人区| 丰满人妻熟妇乱又伦精品不卡| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久视频播放| 无限看片的www在线观看| 亚洲精品色激情综合| 97超级碰碰碰精品色视频在线观看| av免费在线观看网站| 国产一区二区在线av高清观看| 欧美成人免费av一区二区三区| 久久婷婷成人综合色麻豆| 亚洲欧美一区二区三区黑人| 人人妻人人看人人澡| 国产精品久久久人人做人人爽| 国产成人一区二区三区免费视频网站| 国产视频一区二区在线看| 国产精品爽爽va在线观看网站| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 亚洲五月婷婷丁香| 少妇熟女aⅴ在线视频| 亚洲熟妇熟女久久| 最新在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲av成人av| 在线观看免费午夜福利视频| 国产成+人综合+亚洲专区| 日韩欧美三级三区| 午夜福利成人在线免费观看| 色老头精品视频在线观看| 久久精品aⅴ一区二区三区四区| 久久亚洲真实| 国产精品永久免费网站| 曰老女人黄片| 久久这里只有精品中国| 亚洲一码二码三码区别大吗| svipshipincom国产片| 欧美激情久久久久久爽电影| 亚洲美女视频黄频| 青草久久国产| 一二三四在线观看免费中文在| 欧美性猛交黑人性爽| 天堂影院成人在线观看| 国产主播在线观看一区二区| 中出人妻视频一区二区| 亚洲精品一区av在线观看| 亚洲五月婷婷丁香| 在线观看66精品国产| 嫩草影院精品99| 久久人妻av系列| 亚洲真实伦在线观看| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久久久久| 欧美乱码精品一区二区三区| 精华霜和精华液先用哪个| 91av网站免费观看| 免费搜索国产男女视频| 美女扒开内裤让男人捅视频| 白带黄色成豆腐渣| 久久婷婷人人爽人人干人人爱| 后天国语完整版免费观看| 亚洲午夜精品一区,二区,三区| 国产一区二区激情短视频| 正在播放国产对白刺激| av在线播放免费不卡| 狂野欧美白嫩少妇大欣赏| 亚洲成人免费电影在线观看| 中国美女看黄片| 国产精品久久久久久人妻精品电影| 在线国产一区二区在线| 国产一区二区在线av高清观看| www.熟女人妻精品国产| 免费看a级黄色片| tocl精华| 国产精品1区2区在线观看.| 国产精品精品国产色婷婷| 国产探花在线观看一区二区| 两个人免费观看高清视频| 首页视频小说图片口味搜索| 国产av在哪里看| 国产熟女午夜一区二区三区| 欧美日韩黄片免| 中亚洲国语对白在线视频| 一区二区三区国产精品乱码| 国产一区二区激情短视频| 久久人妻av系列| 久久精品夜夜夜夜夜久久蜜豆 | 成人特级黄色片久久久久久久| 99国产精品一区二区蜜桃av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲中文av在线| 国产激情久久老熟女| 国产一区二区三区视频了| 亚洲熟妇中文字幕五十中出| 99热6这里只有精品| 国产亚洲精品av在线| 18禁黄网站禁片免费观看直播| 免费看十八禁软件| 日日夜夜操网爽| 亚洲成人中文字幕在线播放| 成年免费大片在线观看| 97碰自拍视频| 在线观看免费午夜福利视频| 国产精品综合久久久久久久免费| 别揉我奶头~嗯~啊~动态视频| 中文亚洲av片在线观看爽| 欧美zozozo另类| 久久人妻av系列| av在线播放免费不卡| 亚洲免费av在线视频| 又粗又爽又猛毛片免费看| 好男人在线观看高清免费视频| 国产视频一区二区在线看| 亚洲黑人精品在线| 白带黄色成豆腐渣| 一级片免费观看大全| 成人高潮视频无遮挡免费网站| 嫩草影院精品99| 国产精品野战在线观看| 波多野结衣高清作品| 久久午夜综合久久蜜桃| 国产高清激情床上av| 国产单亲对白刺激| 此物有八面人人有两片| 别揉我奶头~嗯~啊~动态视频| 在线观看www视频免费| 日韩大码丰满熟妇| 一区二区三区激情视频| 香蕉国产在线看| 亚洲中文字幕日韩| 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 久久中文看片网| 国产区一区二久久| 国产精品亚洲一级av第二区| 国产欧美日韩精品亚洲av| 欧美3d第一页| 国产亚洲精品久久久久久毛片| 天天添夜夜摸| 日韩欧美国产在线观看| 亚洲全国av大片| 日韩高清综合在线| 91九色精品人成在线观看| 亚洲国产欧美人成| 在线国产一区二区在线| 国产精品一区二区三区四区久久| 国产亚洲欧美98| 久久久久久免费高清国产稀缺| 18禁美女被吸乳视频| 国产av一区在线观看免费| 久久婷婷成人综合色麻豆| av视频在线观看入口| 精品国内亚洲2022精品成人| 一级a爱片免费观看的视频| 欧美人与性动交α欧美精品济南到| 成人国产综合亚洲| 国产亚洲欧美在线一区二区| 蜜桃久久精品国产亚洲av| 久久久久久人人人人人| 一个人观看的视频www高清免费观看 | 神马国产精品三级电影在线观看 | 一级黄色大片毛片| 久久久精品欧美日韩精品| 男女那种视频在线观看| 午夜福利欧美成人| 日本在线视频免费播放| 免费电影在线观看免费观看| 精品国产乱子伦一区二区三区| 日韩欧美三级三区| 成人一区二区视频在线观看| 国产成年人精品一区二区| 黄色a级毛片大全视频| 天堂动漫精品| 欧洲精品卡2卡3卡4卡5卡区| 色在线成人网| 久久久精品国产亚洲av高清涩受| 在线观看免费视频日本深夜| 亚洲av日韩精品久久久久久密| 法律面前人人平等表现在哪些方面| 国产精品综合久久久久久久免费| 免费人成视频x8x8入口观看| 看黄色毛片网站| 日韩欧美一区二区三区在线观看| 99热6这里只有精品| 国产一级毛片七仙女欲春2| 69av精品久久久久久| 婷婷亚洲欧美| 十八禁人妻一区二区| 日本黄色视频三级网站网址| 在线观看日韩欧美| 精品人妻1区二区| 又爽又黄无遮挡网站| 亚洲国产精品sss在线观看| 可以免费在线观看a视频的电影网站| 亚洲成人中文字幕在线播放| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 国产探花在线观看一区二区| 精品久久久久久久毛片微露脸| 欧美日韩国产亚洲二区| 亚洲国产精品成人综合色| 777久久人妻少妇嫩草av网站| 一边摸一边抽搐一进一小说| 亚洲 欧美 日韩 在线 免费| 日韩有码中文字幕| 精品人妻1区二区| 免费看日本二区| 一级a爱片免费观看的视频| 美女扒开内裤让男人捅视频| 久久中文看片网| 脱女人内裤的视频| 视频区欧美日本亚洲| 久久久久久久久中文| 午夜福利免费观看在线| 国产伦一二天堂av在线观看| 狠狠狠狠99中文字幕| 亚洲一区中文字幕在线| 亚洲五月婷婷丁香| 91字幕亚洲| 亚洲成a人片在线一区二区| 欧美av亚洲av综合av国产av| 国产精品永久免费网站| 国产欧美日韩精品亚洲av| 精品高清国产在线一区| 亚洲欧美一区二区三区黑人| 国产精品一区二区精品视频观看| 一本一本综合久久| 精品人妻1区二区| 亚洲第一欧美日韩一区二区三区| 两个人免费观看高清视频| 两个人看的免费小视频| 午夜福利免费观看在线| 亚洲国产精品合色在线| 丝袜人妻中文字幕| 亚洲一区二区三区色噜噜| 日本一区二区免费在线视频| 午夜免费成人在线视频| 91麻豆av在线| 免费av毛片视频| 啦啦啦观看免费观看视频高清| 亚洲国产欧美人成| 日本熟妇午夜| 男人的好看免费观看在线视频 | 久久草成人影院| 国内揄拍国产精品人妻在线| 亚洲av成人精品一区久久| av有码第一页| 亚洲欧美日韩东京热| 1024视频免费在线观看| 国产精品综合久久久久久久免费| 国产成人欧美在线观看| 中文字幕人成人乱码亚洲影| 欧美大码av| 国产一区二区激情短视频| 成人永久免费在线观看视频| 久久精品91无色码中文字幕| 悠悠久久av| 夜夜夜夜夜久久久久| 色噜噜av男人的天堂激情| 精品第一国产精品| 淫秽高清视频在线观看| 老熟妇乱子伦视频在线观看| 国产亚洲av嫩草精品影院| 久久精品人妻少妇| 色av中文字幕| 亚洲熟妇中文字幕五十中出| 亚洲一区二区三区不卡视频| 国模一区二区三区四区视频 | 麻豆国产97在线/欧美 | 听说在线观看完整版免费高清| 色老头精品视频在线观看| 在线观看免费日韩欧美大片| 国产精品影院久久| 亚洲国产精品久久男人天堂| 黄色视频不卡| 亚洲va日本ⅴa欧美va伊人久久| 在线观看午夜福利视频| 国产三级中文精品| 国产精品久久电影中文字幕| 国产av不卡久久| 欧美乱色亚洲激情| 这个男人来自地球电影免费观看| 男人舔女人的私密视频| 亚洲精品国产一区二区精华液| 久久精品国产清高在天天线| 日韩欧美国产在线观看| 亚洲激情在线av| 中文字幕最新亚洲高清| 国产av一区二区精品久久| 欧美日韩亚洲国产一区二区在线观看| 成人三级黄色视频| 久久久久久亚洲精品国产蜜桃av| 人人妻人人澡欧美一区二区| av片东京热男人的天堂| 国产av在哪里看| 国产亚洲欧美98| 国产av在哪里看| 色尼玛亚洲综合影院| 19禁男女啪啪无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久久九九精品影院| 国产男靠女视频免费网站| 高清在线国产一区| 90打野战视频偷拍视频| 国产精品免费视频内射| 亚洲五月婷婷丁香| 亚洲一区二区三区色噜噜| 99在线视频只有这里精品首页| 不卡一级毛片| 国产69精品久久久久777片 | tocl精华| 美女免费视频网站| 中文字幕久久专区| а√天堂www在线а√下载| 精品一区二区三区视频在线观看免费| 欧美丝袜亚洲另类 | 午夜精品久久久久久毛片777| 黄色成人免费大全| 又大又爽又粗| 国产一区二区三区视频了| 男插女下体视频免费在线播放| 老熟妇乱子伦视频在线观看| 啪啪无遮挡十八禁网站| 色老头精品视频在线观看| 久久久久久久精品吃奶| 亚洲欧美日韩东京热| 久久久久久久精品吃奶| 成年版毛片免费区| 国产精品 国内视频| 黑人巨大精品欧美一区二区mp4| 亚洲国产欧美一区二区综合| 国产视频一区二区在线看| aaaaa片日本免费| 1024香蕉在线观看| 国产伦在线观看视频一区| 在线观看免费日韩欧美大片| 欧美av亚洲av综合av国产av| 午夜精品一区二区三区免费看| 18禁美女被吸乳视频| 人妻夜夜爽99麻豆av| 在线看三级毛片| 国产成人系列免费观看| 蜜桃久久精品国产亚洲av| 亚洲自偷自拍图片 自拍| 大型黄色视频在线免费观看| 欧美日韩一级在线毛片| 欧美一区二区精品小视频在线| 18禁黄网站禁片午夜丰满| 国产精品久久久久久亚洲av鲁大| 天堂影院成人在线观看| 久久人妻福利社区极品人妻图片| 一区二区三区激情视频| 中出人妻视频一区二区| 女生性感内裤真人,穿戴方法视频| 精品福利观看| xxxwww97欧美| 亚洲18禁久久av| 不卡av一区二区三区| 国产真实乱freesex| 久久精品aⅴ一区二区三区四区| 色综合婷婷激情| 嫁个100分男人电影在线观看| 欧美av亚洲av综合av国产av| 成人一区二区视频在线观看| 日本 欧美在线| 人人妻人人澡欧美一区二区| 淫妇啪啪啪对白视频| 99国产精品一区二区蜜桃av| 此物有八面人人有两片| 欧美大码av| 精品免费久久久久久久清纯| 国产爱豆传媒在线观看 | 国产主播在线观看一区二区| netflix在线观看网站| 久久久久九九精品影院| 18禁黄网站禁片免费观看直播| 丁香六月欧美| 日韩大码丰满熟妇| 亚洲av成人不卡在线观看播放网| 久久久久性生活片| 日韩精品免费视频一区二区三区| 日本三级黄在线观看| 国产私拍福利视频在线观看| 日日爽夜夜爽网站| 欧美日韩乱码在线| 一级黄色大片毛片| 麻豆成人午夜福利视频| 日韩精品青青久久久久久| 精品欧美一区二区三区在线| 97碰自拍视频| 天堂√8在线中文| 国产精品野战在线观看| or卡值多少钱| 老司机深夜福利视频在线观看| 一个人免费在线观看电影 | 免费一级毛片在线播放高清视频| 久久久久久久午夜电影| 中文字幕高清在线视频| 免费看美女性在线毛片视频| 国产精品一区二区三区四区免费观看 | 深夜精品福利| 黄色毛片三级朝国网站| 国产精品日韩av在线免费观看| 国产熟女xx| 69av精品久久久久久| 熟女电影av网| 久久久国产成人免费| 又粗又爽又猛毛片免费看| 国产精品亚洲av一区麻豆| 中文字幕人成人乱码亚洲影| www.熟女人妻精品国产|