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

    An investigation on improving the homogeneity of plasma generated by linear microwave plasma source with a length of 1550 mm

    2021-02-27 09:16:54JichengZHOU周繼承WeiXU徐偉andTechaoCHEN陳特超
    Plasma Science and Technology 2021年2期
    關(guān)鍵詞:徐偉

    Jicheng ZHOU (周繼承), Wei XU (徐偉) and Techao CHEN (陳特超)

    1 School of Energy Science and Engineering, Central South University, Changsha 410083, People’s Republic of China

    2 Hunan Red Solar Photoelectricity Science and Technology Co.,Ltd,Changsha 410083,People’s Republic of China

    Abstract To develop a larger in-line plasma enhanced chemical vapor deposition (PECVD) device, the length of the linear microwave plasma source needs to be increased to 1550 mm.This paper proposes a solution to the problem of plasma inhomogeneity caused by increasing device length.Based on the COMSOL Multiphysics, a multi-physics field coupling model for in-line PECVD device is developed and validated.The effects of microwave power, chamber pressure, and magnetic flux density on the plasma distribution are investigated, respectively, and their corresponding optimized values are obtained.This paper also presents a new strategy to optimize the wafer position to achieve the balance between deposition rate and film quality.Numerical results have indicated that increasing microwave power and magnetic flux density or decreasing chamber pressure all play positive roles in improving plasma homogeneity,and among them,the microwave power is the most decisive influencing factor.It is found that the plasma homogeneity is optimal under the condition of microwave power at 2000 W,chamber pressure at 15 Pa, and magnetic field strength at 45 mT.The relative deviation is within ?3.7% to 3.9%,which fully satisfies the process requirements of the equipment.The best position for the wafer is 88 mm from the copper antenna.The results are very valuable for improving the quality of the in-line PECVD device.

    Keywords: in-line PECVD, linear microwave plasma source, multi-physics field, COMSOL simulation, homogeneity

    1.Introduction

    As a mature technology, plasma enhanced chemical vapor deposition (PECVD) is widely used for the fabrication of various thin-film materials for photovoltaics.At present,there are three main types of PECVD equipment used in solar cell production lines: the planar RF-PECVD, the tubular RFPECVD, and the in-line MW-PECVD [1–3].Owing to its outstanding features such as high deposition rate, large area uniformity, high safety, and excellent controllability [4], the in-line PECVD device is developed for the deposition of various thin films, like antireflection layers of silicon nitride,back passivation layers of AlOxand amorphous silicon layers for heterojunction cells [5–8].

    The linear microwave plasma source is the kernel of the in-line PECVD device.Petaschet al[9]studied on the coaxial line-coupled cavities,and developed a plasma generator called Duo-Plasmaline.Liehret al[10]optimized the Duo-Plasmaline by spirally winding a copper tube around the inner conductor surface to facilitate the transfer of microwave energy and increase the plasma density.Kaiseret al[11] redesigned the outer cavity of the linear plasma source to improve the uniformity and electron density of the linear microwave plasma source.Katoet al[12]suggested that an external magnetic field could effectively increase plasma density.This was because the magnetic field could increase the confinement and the efficient collision frequency of the charged particles under Lorenz force.Based on this, Schlemmet al[13] developed the linear microwave plasma source used in the in-line PECVD device for silicon nitride deposition of solar cells, and realized its application in industry.The uniform plasma with large area can be generated easily from multiple linear microwave plasma sources that are aligned in parallel.To keep pace with the development of the PV (Photovoltaic) industry, the linear plasma sources need to be evolved towards better performance(>5% uniformity over large areas) and higher productivity in order to reduce production costs and match new cell concepts and technologies.

    It is known that higher productivity can be achieved by improving deposition rate or enlarging equipment size.There have been a few comprehensive works focusing on improving deposition rate, including optimization of the structure and process parameters [14–16].In fact, enlarging the size of the device is supposed to be a more efficient approach.In recent years,the length of linear microwave plasma source in in-line PECVD device has increased from 1100 mm to 1340 mm(with a 4×6 graphite frame),which has been widely applied in the industry.It is aimed to develop the next generation linear microwave plasma source with a length of 1550 mm(with a 4×7 graphite frame).That means the production capacity would be increased by 17%.However, the most prominent problem is the decrease in plasma uniformity,causing inhomogeneity of the target film, which cannot meet the PECVD process requirement.

    As indicated by Hübneret aland Ganesanet al[17,18],the uniformity of plasma generated from linear microwave source is mainly influenced by the device structure, microwave power,chamber pressure,and the magnetic flux density.However,previous studies have mainly focused on laboratory equipment, and both the size and process conditions are dramatically different from linear microwave plasma sources in the in-line PECVD device.This will result in different or even opposite effects on plasma uniformity influenced by certain parameters, such as chamber pressure.Besides, the ultimate difference from smaller devices is that the linear microwave plasma source with a length of 1550 mm is difficult to optimize by adjusting a single parameter, but requiring adjustment of multiple parameters simultaneously.Therefore, this paper aims to study the influence of microwave power,chamber pressure,and magnetic flux density on plasma uniformity and to obtain an optimal parameter recipe by adjusting the parameters mentioned above.

    In the industry, obtaining the optimal parameter recipe mainly relies on extensive trial and experience work, which would cost much time and money.Therefore, numerous numerical studies have been adopted to explore the performance of linear microwave plasma sources.Rahimiet al[19]presented a two-dimensional(2D)model of a linear coaxial plasma source and investigated the propagation of electromagnetic wave.Chenet al[20] developed a selfconsistent 2D models based on the COMSOL Multiphysics.Their work unveiled the characteristics and available propagation mechanism of linear coaxial microwave plasma.Besides, Obrusníket al[21] proposed a mixed 2D/3D fluid model and demonstrated the influence of the microwave shielding on the plasma and particles.In this paper, a multiphysics field coupling model for in-line PECVD device is developed based on the simulation software, COMSOL Multiphysics.It is robust in multi-field coupling with specialized plasma module.More information about COMSOL Multiphysics can be found in [22].

    This paper is generally formed as below.The numerical model and simulation conditions are presented in section 2,with verifying the validity of this model.The section 3 explores the plasma distribution by increasing the length of the linear microwave plasma source to 1550 mm.Then the effects of microwave power,chamber pressure,and magnetic flux density on plasma homogeneity are investigated,obtaining an optimal parameter recipe.The summary is given in section 4.

    2.Numerical simulation and verification

    2.1.Model simplification

    Figure 1(a) presents a schematic diagram of the in-line PECVD.The microwave at 2.45 GHz, generated by two microwave magnetrons, passes the microwave mode converter which is designed to change the microwave from TE10mode to TEM mode into the chamber.The working gas (e.g.argon)enters the chamber through the gas shower,then being excited into plasma around the quartz tube.The magnetic field,generated by the external magnet,increases the collision frequency of charged particles, contributing to improve the plasma density and uniformity [23].The working gas, under the action of high-energy electrons in the plasma, forms a large number of highly reactive surface reaction precursors,which eventually form a film on the wafer surface by chemisorption.The wafers are placed on the graphite frame.The heating plate is used to regulate the surface temperature of the wafer to ensure the quality of the film [24], and the vacuum pump is used to maintain chamber pressure.

    The model needs to be simplified while keeping it feasible to investigate the effect of relevant parameters on plasma homogeneity.First of all,it is assumed that microwave is fed directly from the ports on both sides of the quartz tube in TEM mode[20],and the microwave power used in the model refers to average microwave power.After that, due to the subtle influence of the gas velocity field on the plasma discharge characteristics, the flow field and the region under the heating plate are not calculated in this model [21].Finally,argon is employed instead of the complex process gas.Ultimately a 2D simplified model is developed, as presented in figure 1(b).

    Figure 1.(a) Schematic configuration of in-line PECVD.(b) Schematic diagram of the simplified model.

    2.2.Mathematical model formulas

    The 2D self-consistent argon plasma model is developed based on COMSOL Multiphysics.The whole model containing microwave, plasma, magnetic, and heat field are solved self-consistently in a frequency-transient study.The electromagnetic waves are solved in the frequency domain,and all plasma variables are solved in the time domain.

    The electric field satisfies the following wave equation[25]:

    whereErepresents the amplitude of electric field and ω represents the angular frequency of the microwave.k0means the wavenumber in free space,andσare permeability, permittivity, and conductivity of the material, respectively.μ0andε0represent permeability,permittivity in vacuum.μrandεrare relative permeability and relative permittivity.Equation (1) is solved in air, quartz, and plasma.Both air and quartz are non-conductive (σ= 0) and non-magnetic (μr= 1).The relative permittivities areandrespectively, and the electromagnetic properties of the plasma are calculated from the plasma and magnetic field module.

    For magnetized high-frequency plasma, the plasma conductivity σ needs to be expressed in full tensor form.The inverse of the plasma conductivity is denoted as follows:

    whereνenand j represent the frequency of electron–neutral particle collisions and the imaginary symbol, respectively;Bsx,Bsy,Bszare three external magnetic field components in thex,y,andzdirections, respectively.The relative permittivity and permeability of plasma are given in equations (4) and (5) [26]

    In the plasma modules, the drift-diffusion equations describing electron density and electron energy density are given as follows:

    Table 1.Reactions scheme for the argon plasma.

    and

    whereRedenotes the electron rate expression.Esrepresents the electric field,which is given by equation(1),umeans the neutral fluid velocity vector anduis set to zero in this model.Qdescribes heating of the electrons.μeandDeare the electron mobility and diffusivity, respectively.μεandεDare the electron energy mobility and diffusivity, respectively.μeshould represent in full tensor form:

    Plasma potential calculation by Poisson’s equation:

    Here,V,ρa(bǔ)re potential and space charge density, respectively.The space charge density is calculated from the number density of electrons and other charged particles:

    In the relation above, the indexkruns over all species that enter the reactions.

    Figure 2.Comparisons between the numerical results calculated from the model and the experimental results in the literature [9].

    The mass fraction of heavy species is derived by solving the mixed-average diffusion formula:represent the mass fraction of heavy species, diffusive flux vector, rate expression for speciesk, respectively.ρrepresents the density of the mixture.

    The gas temperature distributionTgis obtained from the heat equation in the form:

    whereλis the thermal conductivity, andQhis the heat sources from neutrals.

    All the reactions that are considered in the model are presented in table 1.To simplify the model, assuming that only four species exist in the argon plasma, i.e.electrons (e),ground state atoms (Ar), metastable particles (Ars), and ions(Ar+).The rate constant is computed from cross-section data by the following equation:

    Figure 3.Under the standard parametric recipe, linear microwave plasma source length increased to 1550 mm.(a) Plasma density,(b) electron temperature, (c) gas temperature distribution.

    whereσkis the collision cross-section, which can be obtained from database [31].frepresents the EEDF with Maxwellian distribution.

    2.3.Model validation

    In this paper,the model accuracy is validated by comparing with the extensive and detailed experimental results from literature[9].In order to verify the accuracy of the model,the simulation conditions need to be kept consistent with literature[9].It means that the linear microwave plasma source length is set to 600 mm,the microwave power is 210 W and the chamber pressure is 50 Pa.Figure 2 shows the comparison between the simulation results and the experimental results.By comparing the scale and distribution of electron density along the axial direction at distances of 10,20,and 30 mm from the quartz tube,it is found that the simulation results are in good agreement with the experimental results.Moreover, the simulation results show that the radial electron density attenuates with increasing radial distance,which is also in accordance with the experimental results in literature [17].The results above demonstrate that the model is valuable for improving the chamber design and process parameter of linear microwave plasma source.

    3.Results and discussion

    Figure 4.Plasma density distribution and relative deviation at the monitoring line (r = 40 mm).

    After validating the accuracy of the model, we simulate the plasma distribution by increasing the length of the linear microwave plasma source to 1550 mm with a standard parametric recipe.For example, microwave power is set at 1800 W, chamber pressure at 25 Pa, and magnetic flux density at 30 mT(For existing 1340 mm length linear microwave plasma sources).The results of plasma density, electron temperature, and gas temperature distribution are shown in figure 3.From figure 3(a),it can be seen that the plasma density decreases axially after the length of the linear microwave plasma source being increased up to 1550 mm under the standard parametric recipe.Particularly in the intermediate regions of the device, the plasma is significantly insufficient.That is because the microwave is emitted from the magnetron at both ends of the device, then being fed from both ends of the quartz tube after passing the waveguide.Subsequently,the microwave interacts with the gas in the chamber, causing the ionization of gas.As microwave propagates to the intermediate region of the device, the microwave energy is consumed and gradually decays under the interaction between microwave and plasma.Thus, the probability of ionization of the neutral particles is reduced by direct absorption of the microwave energy.On the other hand,due to the lower energy of the electrons,the probability of ionization of the neutral particles by collision with the electrons is greatly reduced.Comparing figures 3(a) and (b), it is known that a high plasma density corresponds to a high electron temperature,which indicates that the increase in electron temperature is mainly due to inelastic collisions between electrons and heavy particles.In addition,it can be seen that both electron density and electron temperature decay along the radial direction, but the attenuation amplitude of electron temperature is much larger than that of electron density.Its deeper physical mechanism can be explained as follows: as the gas around the quartz tube discharges by absorbing microwave energy, the microwave energy decays exponentially owing to being absorbed and consumed along the normal direction.Therefore, when the gas molecules are close to the quartz tube,they will absorb more microwave energy, so that the ionization degree, the electron energy, as well as the probability of secondary ionization through electron collision are high, thus resulting in higher electron density and temperature.On the contrary,when the gas molecules are away from the quartz tube,they will maintain the discharge mainly by ionization through collision with the electrons.Therefore, the ionization degree of electrons with less energy will be relatively low, resulting in lower electron density and temperature.Figure 3(c) shows the temperature distribution of gas inside the chamber.Such distribution in the chamber is far more uniform than the electron temperature, with a maximum value of about 650 K, which is consistent with the actual PECVD process temperature range.

    To quantify the variation of plasma homogeneity, the monitoring line is selected asr= 40 mm, introducing the relative deviation:

    nedenotes the plasma density at the point of measurement,andis the average plasma density on the monitoring line.Figure 4 presents the corresponding results, and it can be seen that the relative deviation of the plasma density ranges from?90%to 13.8%,which is unacceptable for the PECVD process.

    3.1.Limitations on microwave power

    The relationship between microwave power and plasma homogeneity is shown in figure 5.As shown in figures 5(a)and(b), along with the increase in microwave power, the plasma keeps extending from the two sides of the device to the intermediate region, and the density and uniformity of the plasma also keep increasing.It means that increasing microwave power is beneficial for plasma homogeneity, which is consistent with our experience.According to the work of Chenet al[32], the generation and propagation of the plasma in linear microwave plasma source mainly depend on the interaction of the penetration part,transmission part and absorption part of the microwave electromagnetic field.Therefore,higher microwave power generates higher electron density, making the reflectivity of the electromagnetic field of plasma higher.That is to say,the plasma behaves more like a metal,allowing greater fractions of the microwave energy to propagate along the quartz tube surface.Consequently, increasing the microwave power can significantly improve the plasma uniformity.Figure 5(c) presents the average density of the plasma and relative deviation at different microwave powers.It is evident that with the increase in microwave power,the average plasma density is almost proportional to the microwave power,and the range of relative deviation of the plasma is also significantly narrowed.Plasma homogeneity is optimal when the microwave power increases to 2200 W, with a relative deviation of?10.2% to 6.1%.However, due to the limitations of the microwave generator, the microwave power can only be increased to 2000 W at most.In this case,the relative deviation of plasma density is ?23.1% to 11.9%, which is unable to satisfy the process requirements of the in-line PECVD device.

    3.2.Influence of chamber pressure

    Figures 6(a)and(b)illustrate the plasma distribution both in the chamber and at the monitoring line under different chamber pressures, respectively.The simulation results for a chamber pressure of 25 Pa have been shown in figure 3(a), so they are omitted in figure 6(a)to avoid repetition.It can be found that as the chamber pressure reduces from 25 to 15 Pa, the plasma becomes more uniform in the axial direction while becoming more concentrated in the radial direction.When the chamber pressure is increased to 30 Pa, the plasma homogeneity will decrease.Thus,reducing the chamber pressure helps to generate a more uniform plasma,but results in a lower plasma density.It can be explained by the fact that due to the decrease in pressure,the collision frequency between electron–neutral particles decreases.On one hand,the loss of energy caused by collisions also decreases, which is beneficial for electron diffusion, leading to excellent uniformity.On the other hand, the number of emerging electrons caused by collisions decreases,resulting in a lower plasma density.Due to the limitations of the vacuum pump, the chamber pressure could not be reduced less than 15 Pa.Corresponding relative deviation of the plasma density is calculated to be within ?11.3% to 8.2%.

    Figure 5.Plasma density distribution (a) in the chamber and (b) at the monitoring line, (c) average plasma density and relative deviation under different microwave powers.

    Figure 6.Plasma density distribution (a) in the chamber and (b) at the monitoring line under different pressures.

    3.3.Effects of magnetic flux density

    Figure 7.Plasma density distribution (a) in the chamber and (b) at the monitoring line at different magnetic field strengths.

    Table 2.Average plasma density and relative deviation for different conditions.

    In this section,we turn our attention to the effect of magnetic flux density on the uniformity of the plasma generated by the linear microwave plasma source.The magnetic flux density varies from 15 to 60 mT.Figures 7(a) and (b) show the corresponding changes in plasma distribution for different magnetic flux densities inside the chamber and at the monitoring line.As shown in figure 7(a), when the magnetic flux density is 15 mT, the plasma density is concentrated at the two sides of the chamber.When the magnetic flux density is adjusted to 45 mT, the plasma distribution in the chamber becomes more uniform,and the plasma density in the middle region increases significantly at the same time.But when the magnetic flux density continues to be increased to 60 mT,the plasma distribution will show a tendency to extend to the two sides again.As shown in figure 7(b),the plasma density raises along with the increasing of magnetic flux density, because higher magnetic flux density can reduce the complexation of charged particles at the walls.Increasing the magnetic flux density appropriately contributes to improving plasma uniformity, but excessive magnetic flux density leads to poor uniformity instead.It may be due to the fact that electromagnetic waves gradually attenuate as they are propagated in the axial direction.But an appropriately higher magnetic flux density helps to suppress the decay of electromagnetic waves,which means that the electromagnetic waves can travel further in the axial direction.As a result, increasing the magnetic field appropriately can improve the plasma uniformity to some extent, and 45 mT is the best choice.In this case, the relative deviation of plasma is within ?14.2% to 16.3%.

    By the analysis above, after optimizing the microwave power, chamber pressure, and magnetic flux density, respectively, the properties of plasma are given in table 2.According to the data in table 2, it is shown that higher microwave power, lower chamber pressure, and relatively higher magnetic field all have positive effects on improving plasma uniformity.Among them,the microwave power has a prominent enhancement on plasma uniformity.It is because the microwave power directly affects the transport of microwaves in the plasma, which determines the characteristics of the plasma.It is also found that only optimizing one parameter cannot make plasma homogeneity satisfy the process requirements of the in-line PECVD device.

    3.4.Optimal parameter recipe and wafer position

    In this section, the variation of plasma homogeneity for different parameter combinations is discussed to explore the optimal parameter recipe.It should be noted that the adjustments of parameter are based on previously standard parametric recipe (i.e.microwave power is set at 1800 W,chamber pressure at 25 Pa, and magnetic flux density at 30 mT),the results are shown in figures 8(a)–(d).First of all,as described in section 3.1,the microwave power is increased up to 2000 W, and the corresponding plasma density distribution is plotted in figure 8(a).It is shown that the plasma density in the intermediate region is significantly lower than the average,with relative deviations ranging from ?23.1%to 11.9%, which indicates that increasing the microwave power to 2000 W is inadequate to sustain the energy attenuation during microwave transmission over a length of up to 1550 mm.However,in figure 8(b),after decreasing the chamber pressure to 15 Pa, the plasma homogeneity improved dramatically, especially in the intermediate region where the plasma density is higher than the average.The relative deviation of plasma density ranges from ?4.2% to 6.1%.

    Figure 8.Plasma density distribution and relative deviation when parameters are adjusted to (a) microwave power 2000 W, (b) microwave power 2000 W and chamber pressure 15 Pa,(c)microwave power 2000 W and chamber pressure 15 Pa and magnetic flux density 45 mT,(d)chamber pressure 15 Pa and magnetic flux density 45 mT.

    On the basis of figures 8(b), (c) presents the plasma density distribution when the magnetic flux density is adjusted to 45 mT.It can be seen that the plasma density is distributed very evenly, the plasma density at two sides of the device is slightly higher than that in the intermediate region.The relative deviation is in the range of ?3.7% to 3.9%, which fully satisfies the requirements of process uniformity of in-line PECVD devices (±4%).By the way,the average plasma density is the highest in this case.When the microwave power is fixed, if the chamber pressure is adjusted to 15 Pa and the magnetic flux density to 45 mT,the corresponding distribution of plasma density is shown in figure 8(d).The plasma density relative deviation is within ?7% to 8%.Consequently, the optimal parameter recipe is to change the microwave power to 2000 W, the chamber pressure to 15 Pa,and the magnetic flux density to 45 mT simultaneously.

    The position of the wafer in the in-line PECVD device is crucial for the film quality.If the wafers are too close to the linear microwave plasma source, massive charged ion bombardment will damage the film quality.But if they are placed too far away, the lower density of excited state atoms would lead to a lower film deposition rate.Therefore, the best process position for silicon wafers should be where the value of the density of excited state atoms divided by electron density is maximized.Figure 9 shows the distribution of electron density and excited state atoms density atz= 800 mm under the optimal parameter recipe.It can be seen from figure 9 that the valuenwArs/neis maximized atr= 88 mm.It means that when the wafer is 88 mm away from the copper antenna, a proper balance between deposition rate and film quality can be achieved.

    Figure 9.Distribution of electron density and excited state atom density under the optimal parameter recipe (at z = 800 mm).

    4.Conclusion

    Enlarging the length of the linear microwave plasma source in the in-line PECVD from 1340 to 1550 mm could increase production capacity by 17%.Based on the COMSOL Multiphysics, a numerical model of coupling microwave fields,magnetic fields, plasma fields as well as thermal fields has been developed and validated.Optimization of the plasma uniformity is achieved by adjusting the microwave power,chamber pressure, and magnetic flux density.The main findings and conclusions are summarized as follows:

    (1) Under the standard parameter recipe,increasing the length of the linear microwave plasma source to 1550 mm results in a dramatic decrease in plasma homogeneity with a relative deviation within ?90% to 13.8%, which is unacceptable for the PECVD process.Either increasing the microwave power or decreasing the chamber pressure contributes to the homogeneity of the plasma generated by the linear microwave plasma source.An appropriate increase in the magnetic flux density facilitates homogeneous plasma generation, but excessive magnetic flux density will have negative influence on the plasma homogeneity.Among them, the microwave power is the most decisive factor.

    (2) For the linear microwave plasma source with a length of 1550 mm, plasma homogeneity is optimal when simultaneously adjusting the microwave power to 2000 W, the chamber pressure to 15 Pa,and the magnetic flux density to 45 mT,and the relative deviation of the plasma density is within ?3.7%to 3.9%,which fully satisfies the process uniformity requirements of the in-line PECVD device.

    (3) Under the optimal parameter recipe, the optimum position of the wafer is 88 mm from the copper antenna.

    Acknowledgments

    This work was funded by Hunan Provincial Science and Technology Department ‘Key Field Research and Development Project’(No.2019WK2011).We thank B X Zhao,Y Z Lu,H B Zhou and H Yi for discussions and experimental and technical assistance.

    猜你喜歡
    徐偉
    徐偉:野菊
    Most probable transition paths in eutrophicated lake ecosystem under Gaussian white noise and periodic force
    高階時頻變換理論與應(yīng)用
    Study of the tungsten sputtering source suppression by wall conditionings in the EAST tokamak
    武漢工程大學(xué)土木工程與建筑學(xué)院 徐偉 水彩作品選登
    跨越三十年的師生情
    Theory and method of dual-energy x-ray grating phase-contrast imaging?
    Quantification of CP4-EPSPS in genetically modifiedNicotianatabacum leaves by LC-MS/MS with 18O-labeling
    NOVEL HIGH-SPEED FPGA-BASED FFT PROCESSOR
    “2006大學(xué)生年度人物”徐偉,一個貧困大學(xué)生的快樂經(jīng)
    做人與處世(2007年7期)2007-06-30 03:04:36
    97人妻天天添夜夜摸| 国产深夜福利视频在线观看| 好男人视频免费观看在线| 亚洲欧美日韩另类电影网站| 蜜桃国产av成人99| 国产精品国产三级国产专区5o| 亚洲欧美清纯卡通| 欧美成人午夜精品| 午夜老司机福利剧场| 免费观看无遮挡的男女| 色哟哟·www| 日日摸夜夜添夜夜爱| 日本av手机在线免费观看| 久久久国产一区二区| 久久综合国产亚洲精品| 久久久久精品性色| 97人妻天天添夜夜摸| 国产乱来视频区| 美女视频免费永久观看网站| 妹子高潮喷水视频| 啦啦啦啦在线视频资源| 国产精品一区www在线观看| 欧美xxxx性猛交bbbb| 欧美精品亚洲一区二区| 少妇 在线观看| 亚洲五月色婷婷综合| 久久精品久久久久久久性| 女的被弄到高潮叫床怎么办| 国产深夜福利视频在线观看| tube8黄色片| a级毛色黄片| 国产黄色视频一区二区在线观看| 九色亚洲精品在线播放| 国产乱来视频区| 极品少妇高潮喷水抽搐| 一个人免费看片子| 久久婷婷青草| 亚洲,欧美,日韩| 肉色欧美久久久久久久蜜桃| 人成视频在线观看免费观看| 国产成人精品无人区| 国产精品无大码| 欧美bdsm另类| av又黄又爽大尺度在线免费看| 国产男女内射视频| 一级a做视频免费观看| 亚洲国产欧美在线一区| xxxhd国产人妻xxx| 日韩欧美精品免费久久| 欧美 亚洲 国产 日韩一| 国产乱来视频区| 在线亚洲精品国产二区图片欧美| 久久久久久久亚洲中文字幕| 欧美激情 高清一区二区三区| 久久女婷五月综合色啪小说| 欧美激情极品国产一区二区三区 | 国产色爽女视频免费观看| 免费少妇av软件| 国产av一区二区精品久久| 欧美激情 高清一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 人人澡人人妻人| 日本vs欧美在线观看视频| 五月玫瑰六月丁香| 激情五月婷婷亚洲| 中文精品一卡2卡3卡4更新| 人人妻人人添人人爽欧美一区卜| 中文字幕av电影在线播放| 国产免费视频播放在线视频| 午夜福利影视在线免费观看| 日韩免费高清中文字幕av| 日韩免费高清中文字幕av| 建设人人有责人人尽责人人享有的| 成年av动漫网址| 欧美97在线视频| 高清视频免费观看一区二区| 色94色欧美一区二区| 91午夜精品亚洲一区二区三区| 久久女婷五月综合色啪小说| 精品酒店卫生间| 伊人亚洲综合成人网| 国产 一区精品| 免费久久久久久久精品成人欧美视频 | 久久婷婷青草| 久久精品国产a三级三级三级| 亚洲国产精品国产精品| 内地一区二区视频在线| 亚洲图色成人| 日日啪夜夜爽| 少妇精品久久久久久久| 男女无遮挡免费网站观看| 少妇被粗大的猛进出69影院 | 亚洲激情五月婷婷啪啪| 高清黄色对白视频在线免费看| av福利片在线| 美女国产视频在线观看| 国产一区二区在线观看av| 精品久久国产蜜桃| 制服人妻中文乱码| 99九九在线精品视频| 少妇 在线观看| 老熟女久久久| 久久国产精品男人的天堂亚洲 | 国产视频首页在线观看| 日韩欧美精品免费久久| 另类精品久久| 伊人久久国产一区二区| 免费黄频网站在线观看国产| 免费在线观看完整版高清| 国产在线视频一区二区| 久久久久久久久久久免费av| 最近中文字幕2019免费版| 国产福利在线免费观看视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品色激情综合| 久久99热这里只频精品6学生| 婷婷色综合www| 亚洲精品第二区| 晚上一个人看的免费电影| 精品一品国产午夜福利视频| 看免费成人av毛片| 国产有黄有色有爽视频| 交换朋友夫妻互换小说| 日韩中文字幕视频在线看片| 最后的刺客免费高清国语| 亚洲四区av| 在线亚洲精品国产二区图片欧美| 亚洲性久久影院| av有码第一页| 亚洲第一区二区三区不卡| 亚洲欧美一区二区三区黑人 | 一级毛片 在线播放| 日韩欧美精品免费久久| videos熟女内射| 久久国产精品大桥未久av| 精品一区二区免费观看| 性色avwww在线观看| 日本91视频免费播放| 在线天堂最新版资源| 极品人妻少妇av视频| 最后的刺客免费高清国语| 精品亚洲成a人片在线观看| 99精国产麻豆久久婷婷| 亚洲av国产av综合av卡| 日本wwww免费看| 国产精品人妻久久久影院| 大片免费播放器 马上看| 国产一级毛片在线| 国产乱人偷精品视频| 久久久久国产网址| 91国产中文字幕| 亚洲av男天堂| 丝袜美足系列| 建设人人有责人人尽责人人享有的| 精品亚洲成国产av| 欧美激情国产日韩精品一区| 人体艺术视频欧美日本| 欧美精品人与动牲交sv欧美| 中文字幕人妻熟女乱码| 丰满乱子伦码专区| av不卡在线播放| 999精品在线视频| 久久久a久久爽久久v久久| 宅男免费午夜| 免费高清在线观看视频在线观看| 少妇的丰满在线观看| 一级毛片我不卡| 国产精品久久久久久精品古装| 夜夜爽夜夜爽视频| 岛国毛片在线播放| 成人手机av| 9色porny在线观看| 亚洲欧美一区二区三区国产| 久久久久久久亚洲中文字幕| 在线精品无人区一区二区三| 亚洲欧美精品自产自拍| 一区二区av电影网| √禁漫天堂资源中文www| 日本爱情动作片www.在线观看| 激情视频va一区二区三区| 成年动漫av网址| 中文字幕另类日韩欧美亚洲嫩草| 精品久久蜜臀av无| 又黄又粗又硬又大视频| 久久久久久人妻| 精品亚洲成国产av| 欧美日本中文国产一区发布| freevideosex欧美| 精品一区在线观看国产| 国产精品99久久99久久久不卡 | 精品国产国语对白av| 国产成人欧美| 国产又色又爽无遮挡免| 久久久久久久久久久久大奶| 免费av中文字幕在线| 午夜激情av网站| 内地一区二区视频在线| 人妻 亚洲 视频| 最近最新中文字幕大全免费视频 | 亚洲天堂av无毛| 国产男女内射视频| 色婷婷久久久亚洲欧美| 国产一级毛片在线| 亚洲精品美女久久av网站| 国产精品一区二区在线观看99| 亚洲熟女精品中文字幕| 多毛熟女@视频| 亚洲成人手机| 国产日韩欧美视频二区| 国产免费又黄又爽又色| 国产探花极品一区二区| 91午夜精品亚洲一区二区三区| 人人妻人人澡人人看| 欧美国产精品一级二级三级| 日韩成人av中文字幕在线观看| 毛片一级片免费看久久久久| 久久久精品区二区三区| 久久狼人影院| 18+在线观看网站| 高清黄色对白视频在线免费看| av播播在线观看一区| 精品第一国产精品| 黄色视频在线播放观看不卡| 久久精品国产鲁丝片午夜精品| 男女无遮挡免费网站观看| 成年人免费黄色播放视频| 国产精品 国内视频| www.色视频.com| 高清不卡的av网站| 色婷婷久久久亚洲欧美| 国产精品久久久久久久久免| videos熟女内射| 精品酒店卫生间| 22中文网久久字幕| 免费在线观看黄色视频的| 色婷婷av一区二区三区视频| 免费观看性生交大片5| 一二三四在线观看免费中文在 | 高清毛片免费看| 在线天堂中文资源库| 日产精品乱码卡一卡2卡三| 欧美人与性动交α欧美精品济南到 | 中文字幕精品免费在线观看视频 | 美女内射精品一级片tv| 人妻系列 视频| 婷婷成人精品国产| 自线自在国产av| 精品人妻在线不人妻| av.在线天堂| 精品国产乱码久久久久久小说| 国产av国产精品国产| 亚洲成人一二三区av| 亚洲,欧美精品.| 亚洲av电影在线进入| 国产一区二区三区av在线| 亚洲欧洲日产国产| 最后的刺客免费高清国语| 欧美3d第一页| 在线观看免费视频网站a站| 欧美人与善性xxx| 国产黄色免费在线视频| 国产高清国产精品国产三级| 色婷婷久久久亚洲欧美| 国产一区有黄有色的免费视频| 精品国产露脸久久av麻豆| 少妇人妻 视频| 久久午夜福利片| 久久久久久久国产电影| 国产探花极品一区二区| 精品国产乱码久久久久久小说| 在线观看美女被高潮喷水网站| 建设人人有责人人尽责人人享有的| 欧美激情 高清一区二区三区| 18在线观看网站| 狂野欧美激情性xxxx在线观看| 少妇熟女欧美另类| 国产成人午夜福利电影在线观看| 国产黄色免费在线视频| 欧美国产精品一级二级三级| 国产精品国产三级国产专区5o| 久久亚洲国产成人精品v| 日韩av在线免费看完整版不卡| 男女高潮啪啪啪动态图| 久久久a久久爽久久v久久| 欧美日韩视频精品一区| 国产亚洲一区二区精品| 亚洲精品视频女| 国产精品久久久久成人av| 99re6热这里在线精品视频| 午夜激情久久久久久久| 丝袜喷水一区| 日本wwww免费看| 丝袜在线中文字幕| 精品视频人人做人人爽| 在线精品无人区一区二区三| 女人精品久久久久毛片| 久久精品国产鲁丝片午夜精品| 中文欧美无线码| 最近2019中文字幕mv第一页| 永久网站在线| 亚洲四区av| 国产白丝娇喘喷水9色精品| 久久影院123| 日韩视频在线欧美| 一本色道久久久久久精品综合| 婷婷色综合www| 九九爱精品视频在线观看| 一本久久精品| 精品国产一区二区三区四区第35| 母亲3免费完整高清在线观看 | 国产成人一区二区在线| 国产精品无大码| 2021少妇久久久久久久久久久| 在线天堂中文资源库| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码 | 99视频精品全部免费 在线| 久久人人爽av亚洲精品天堂| 久久毛片免费看一区二区三区| 少妇的丰满在线观看| 高清毛片免费看| 如日韩欧美国产精品一区二区三区| 九色亚洲精品在线播放| 亚洲精品视频女| 久久久久网色| 精品国产一区二区三区久久久樱花| 91精品国产国语对白视频| 高清视频免费观看一区二区| 免费不卡的大黄色大毛片视频在线观看| 美女国产高潮福利片在线看| 全区人妻精品视频| 一级毛片我不卡| 久久精品人人爽人人爽视色| 久久国产亚洲av麻豆专区| 久久久久久久亚洲中文字幕| 只有这里有精品99| 国产日韩欧美视频二区| 国产免费一级a男人的天堂| 在线观看国产h片| 久久国产精品男人的天堂亚洲 | 日韩电影二区| av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 国产精品蜜桃在线观看| 国产极品天堂在线| 日韩一区二区三区影片| 婷婷色麻豆天堂久久| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 人妻少妇偷人精品九色| 亚洲激情五月婷婷啪啪| 侵犯人妻中文字幕一二三四区| 亚洲av日韩在线播放| 亚洲精品国产av成人精品| 免费av不卡在线播放| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 人人妻人人添人人爽欧美一区卜| 欧美变态另类bdsm刘玥| 免费久久久久久久精品成人欧美视频 | 久久亚洲国产成人精品v| 欧美国产精品一级二级三级| 欧美变态另类bdsm刘玥| 成人漫画全彩无遮挡| 国产午夜精品一二区理论片| 两个人看的免费小视频| 亚洲欧美成人精品一区二区| 久久热在线av| 免费看光身美女| 欧美性感艳星| 99国产精品免费福利视频| 成人18禁高潮啪啪吃奶动态图| 国产午夜精品一二区理论片| 两个人看的免费小视频| 在线精品无人区一区二区三| 久久青草综合色| 欧美+日韩+精品| 两个人看的免费小视频| 看免费成人av毛片| 少妇猛男粗大的猛烈进出视频| 捣出白浆h1v1| 久久 成人 亚洲| 大香蕉久久网| 99久久综合免费| 夫妻性生交免费视频一级片| 黄色一级大片看看| 免费黄网站久久成人精品| 久久久久久久久久人人人人人人| 校园人妻丝袜中文字幕| 三上悠亚av全集在线观看| 你懂的网址亚洲精品在线观看| 久久午夜综合久久蜜桃| 久久久久久人妻| 亚洲一码二码三码区别大吗| 亚洲色图综合在线观看| 亚洲国产看品久久| 在线观看免费日韩欧美大片| 色婷婷久久久亚洲欧美| 女人久久www免费人成看片| 午夜久久久在线观看| 99久久综合免费| 中文字幕精品免费在线观看视频 | 波多野结衣一区麻豆| 一区二区三区乱码不卡18| 九草在线视频观看| 97在线视频观看| 免费黄网站久久成人精品| 欧美最新免费一区二区三区| 啦啦啦在线观看免费高清www| 成年av动漫网址| 午夜日本视频在线| 18禁在线无遮挡免费观看视频| 人妻 亚洲 视频| 99热网站在线观看| 丝袜在线中文字幕| 男女免费视频国产| 午夜福利影视在线免费观看| 最近最新中文字幕大全免费视频 | 亚洲精品,欧美精品| 日韩一区二区三区影片| 最黄视频免费看| 国产成人欧美| 免费在线观看黄色视频的| 内地一区二区视频在线| 99久久中文字幕三级久久日本| 国产一区有黄有色的免费视频| 高清av免费在线| 午夜视频国产福利| 久久久久人妻精品一区果冻| 欧美激情极品国产一区二区三区 | 欧美精品人与动牲交sv欧美| 97精品久久久久久久久久精品| 欧美日韩成人在线一区二区| 女的被弄到高潮叫床怎么办| 少妇人妻 视频| 久久久a久久爽久久v久久| 大码成人一级视频| 亚洲欧美成人精品一区二区| 啦啦啦中文免费视频观看日本| 九草在线视频观看| 一本大道久久a久久精品| 亚洲欧美精品自产自拍| 下体分泌物呈黄色| 丝袜人妻中文字幕| 婷婷色av中文字幕| 久久精品久久久久久噜噜老黄| 丁香六月天网| 日韩av免费高清视频| 最近中文字幕2019免费版| 国产免费视频播放在线视频| 免费看光身美女| 成人18禁高潮啪啪吃奶动态图| 国产深夜福利视频在线观看| 在线观看人妻少妇| 欧美人与性动交α欧美精品济南到 | 亚洲欧美精品自产自拍| 日韩av免费高清视频| 国产欧美另类精品又又久久亚洲欧美| 最黄视频免费看| 男女边摸边吃奶| 美女主播在线视频| 亚洲国产欧美在线一区| 欧美国产精品va在线观看不卡| 在线观看一区二区三区激情| av.在线天堂| 亚洲综合色网址| 国产精品久久久久久精品电影小说| 成年女人在线观看亚洲视频| 亚洲av欧美aⅴ国产| 国产精品人妻久久久影院| 51国产日韩欧美| 人妻一区二区av| 一区二区av电影网| 国产精品一区www在线观看| 蜜桃国产av成人99| 欧美激情国产日韩精品一区| 国产激情久久老熟女| 熟女人妻精品中文字幕| 午夜视频国产福利| 男人操女人黄网站| 国产女主播在线喷水免费视频网站| av免费观看日本| 欧美人与性动交α欧美软件 | 亚洲av电影在线观看一区二区三区| 国产欧美日韩一区二区三区在线| videosex国产| 亚洲内射少妇av| 丰满饥渴人妻一区二区三| 90打野战视频偷拍视频| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 美女脱内裤让男人舔精品视频| 亚洲国产欧美在线一区| 自线自在国产av| 成人毛片60女人毛片免费| 我要看黄色一级片免费的| 亚洲欧美一区二区三区黑人 | 国产乱人偷精品视频| 亚洲av.av天堂| av黄色大香蕉| 国产一区亚洲一区在线观看| 男女无遮挡免费网站观看| 亚洲欧洲日产国产| 国产伦理片在线播放av一区| 亚洲 欧美一区二区三区| 国产成人精品福利久久| 午夜91福利影院| 寂寞人妻少妇视频99o| 黄色视频在线播放观看不卡| 日本av免费视频播放| 纯流量卡能插随身wifi吗| 精品久久久久久电影网| 天天影视国产精品| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 久久精品久久精品一区二区三区| 亚洲av电影在线进入| 老司机亚洲免费影院| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 男人操女人黄网站| 91成人精品电影| 国产精品欧美亚洲77777| 中国美白少妇内射xxxbb| 国产1区2区3区精品| 国产av码专区亚洲av| 美女国产高潮福利片在线看| 免费高清在线观看视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲av男天堂| 秋霞伦理黄片| 免费大片18禁| 美女脱内裤让男人舔精品视频| a 毛片基地| 日本与韩国留学比较| 中文字幕免费在线视频6| 女人精品久久久久毛片| 天天影视国产精品| 美女xxoo啪啪120秒动态图| 久久鲁丝午夜福利片| 性色avwww在线观看| 交换朋友夫妻互换小说| 国产69精品久久久久777片| 日韩精品有码人妻一区| 一级爰片在线观看| 亚洲精品国产色婷婷电影| 国产在线视频一区二区| 中文字幕精品免费在线观看视频 | 国产精品蜜桃在线观看| 少妇的丰满在线观看| 国产成人一区二区在线| 午夜老司机福利剧场| 日韩av不卡免费在线播放| 晚上一个人看的免费电影| 亚洲精品乱久久久久久| 熟女av电影| 国产成人a∨麻豆精品| 国产女主播在线喷水免费视频网站| 久热久热在线精品观看| 女人精品久久久久毛片| 尾随美女入室| 最后的刺客免费高清国语| 国产精品国产三级专区第一集| 成年av动漫网址| 久久久国产一区二区| 成年人免费黄色播放视频| 亚洲精品国产av成人精品| 18在线观看网站| 黑人巨大精品欧美一区二区蜜桃 | 看免费成人av毛片| 9热在线视频观看99| 欧美日韩视频精品一区| 国产免费又黄又爽又色| 一区二区三区四区激情视频| 中文字幕人妻丝袜制服| 亚洲精品成人av观看孕妇| 99热国产这里只有精品6| 精品一区二区三卡| 黄色怎么调成土黄色| 亚洲国产精品一区二区三区在线| 热re99久久精品国产66热6| 亚洲一级一片aⅴ在线观看| 国产又色又爽无遮挡免| 欧美日韩国产mv在线观看视频| 免费高清在线观看视频在线观看| 在线观看人妻少妇| 欧美日韩视频高清一区二区三区二| 国产有黄有色有爽视频| 韩国高清视频一区二区三区| 少妇的丰满在线观看| 亚洲精品美女久久av网站| 久久久久精品性色| 少妇被粗大猛烈的视频| 80岁老熟妇乱子伦牲交| 嫩草影院入口| 免费高清在线观看视频在线观看| 国产成人精品一,二区| 日本av手机在线免费观看| 国产免费福利视频在线观看| 亚洲精品456在线播放app| 一边亲一边摸免费视频| 国产一区二区激情短视频 | 18+在线观看网站| 午夜福利影视在线免费观看| 一区二区三区精品91| 老女人水多毛片| 综合色丁香网| 国产女主播在线喷水免费视频网站| 成年av动漫网址| av天堂久久9| 欧美xxⅹ黑人| 亚洲国产欧美日韩在线播放| 大片电影免费在线观看免费|