• <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
    春色校园在线视频观看| 欧美+日韩+精品| 国产女主播在线喷水免费视频网站 | 日韩av不卡免费在线播放| 色综合亚洲欧美另类图片| 俄罗斯特黄特色一大片| 狂野欧美白嫩少妇大欣赏| 久久久a久久爽久久v久久| 老熟妇乱子伦视频在线观看| 村上凉子中文字幕在线| 日韩欧美免费精品| 91精品国产九色| 日韩在线高清观看一区二区三区| 国产老妇女一区| 男女下面进入的视频免费午夜| 激情 狠狠 欧美| 久久久国产成人免费| 国产又黄又爽又无遮挡在线| 欧美三级亚洲精品| 黑人高潮一二区| 69av精品久久久久久| 又爽又黄无遮挡网站| 麻豆一二三区av精品| 永久网站在线| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 亚洲图色成人| 白带黄色成豆腐渣| 国产精品久久久久久精品电影| 久久6这里有精品| 成人鲁丝片一二三区免费| 国产在线男女| 波多野结衣高清无吗| 亚洲内射少妇av| 99久国产av精品国产电影| 精品久久久久久久久久久久久| 白带黄色成豆腐渣| 国产私拍福利视频在线观看| 亚洲精品影视一区二区三区av| 亚洲av不卡在线观看| 熟妇人妻久久中文字幕3abv| 国产精品无大码| 中国国产av一级| 最近的中文字幕免费完整| 99精品在免费线老司机午夜| 波多野结衣高清无吗| 欧美性猛交黑人性爽| 午夜爱爱视频在线播放| a级一级毛片免费在线观看| 两个人的视频大全免费| avwww免费| 亚洲av熟女| 中出人妻视频一区二区| 熟妇人妻久久中文字幕3abv| 男女做爰动态图高潮gif福利片| 亚洲七黄色美女视频| 男人狂女人下面高潮的视频| h日本视频在线播放| 少妇丰满av| 一个人免费在线观看电影| 99热这里只有是精品在线观看| 18+在线观看网站| 91午夜精品亚洲一区二区三区| 国产蜜桃级精品一区二区三区| 深夜a级毛片| 又爽又黄a免费视频| 午夜精品国产一区二区电影 | 极品教师在线视频| 直男gayav资源| 国产久久久一区二区三区| 亚洲国产精品久久男人天堂| 99热这里只有是精品在线观看| 成熟少妇高潮喷水视频| 国产久久久一区二区三区| 丝袜喷水一区| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品合色在线| 国产精品久久久久久久电影| 欧美一区二区精品小视频在线| 黄片wwwwww| 亚洲av美国av| 18禁在线播放成人免费| 国产精品久久视频播放| 日韩 亚洲 欧美在线| 老师上课跳d突然被开到最大视频| 成人一区二区视频在线观看| 99久久久亚洲精品蜜臀av| 亚洲内射少妇av| 在线免费观看不下载黄p国产| 成人精品一区二区免费| 亚洲真实伦在线观看| av在线天堂中文字幕| 老司机午夜福利在线观看视频| 天天一区二区日本电影三级| 国产高清激情床上av| 国产精品久久视频播放| 成人一区二区视频在线观看| 亚洲成人中文字幕在线播放| 成人无遮挡网站| 天堂影院成人在线观看| 午夜福利高清视频| 免费大片18禁| 99久久成人亚洲精品观看| 九九在线视频观看精品| 深夜精品福利| 99热6这里只有精品| 国产亚洲精品av在线| 午夜影院日韩av| 少妇裸体淫交视频免费看高清| 91精品国产九色| 黄色日韩在线| 国产一区二区激情短视频| 在线免费观看不下载黄p国产| 欧美日韩国产亚洲二区| 99久久中文字幕三级久久日本| 美女cb高潮喷水在线观看| 日本一本二区三区精品| 不卡一级毛片| 天美传媒精品一区二区| 国产aⅴ精品一区二区三区波| 亚洲专区国产一区二区| 国产黄a三级三级三级人| 成人性生交大片免费视频hd| 精品少妇黑人巨大在线播放 | 欧美三级亚洲精品| 久久精品久久久久久噜噜老黄 | 国产精品久久久久久亚洲av鲁大| 婷婷精品国产亚洲av在线| 久久久久久九九精品二区国产| 国产av在哪里看| 少妇丰满av| 亚洲性夜色夜夜综合| 亚洲精品国产av成人精品 | 久久久午夜欧美精品| 欧美日韩精品成人综合77777| 亚洲国产欧洲综合997久久,| 搡老妇女老女人老熟妇| 在线免费观看不下载黄p国产| av黄色大香蕉| 国内精品一区二区在线观看| 国产精品久久久久久亚洲av鲁大| 男女之事视频高清在线观看| 又黄又爽又刺激的免费视频.| 1000部很黄的大片| 69av精品久久久久久| 男女那种视频在线观看| 97人妻精品一区二区三区麻豆| 国产精品久久久久久av不卡| 又黄又爽又刺激的免费视频.| 日韩精品青青久久久久久| 久久综合国产亚洲精品| 国产综合懂色| 精品人妻视频免费看| 嫩草影院入口| 成人精品一区二区免费| 3wmmmm亚洲av在线观看| 亚洲人与动物交配视频| 黄色配什么色好看| 黑人高潮一二区| 99热全是精品| 国产精品亚洲一级av第二区| 我要搜黄色片| 久久久精品大字幕| 婷婷色综合大香蕉| 在线观看午夜福利视频| 日韩欧美在线乱码| 丰满人妻一区二区三区视频av| 简卡轻食公司| .国产精品久久| 久久精品影院6| 国产精品一二三区在线看| 网址你懂的国产日韩在线| 日韩三级伦理在线观看| 岛国在线免费视频观看| 欧美成人一区二区免费高清观看| 欧美最黄视频在线播放免费| 男女下面进入的视频免费午夜| 日韩亚洲欧美综合| 免费看光身美女| 免费不卡的大黄色大毛片视频在线观看 | 天天一区二区日本电影三级| 国产亚洲91精品色在线| 午夜福利成人在线免费观看| 精品人妻偷拍中文字幕| 日日摸夜夜添夜夜爱| av在线老鸭窝| www日本黄色视频网| 国产精品久久久久久av不卡| a级一级毛片免费在线观看| 少妇被粗大猛烈的视频| 久久久国产成人免费| 国产免费男女视频| 国产精品免费一区二区三区在线| 免费在线观看成人毛片| av福利片在线观看| 久久久精品欧美日韩精品| 免费黄网站久久成人精品| 精品一区二区三区人妻视频| 丝袜美腿在线中文| 久久久欧美国产精品| 精品一区二区免费观看| 99热全是精品| 成年女人看的毛片在线观看| 亚洲av熟女| 免费人成在线观看视频色| www日本黄色视频网| 久久久国产成人免费| 国内精品一区二区在线观看| 我要看日韩黄色一级片| 国产在线精品亚洲第一网站| 欧美bdsm另类| 欧美色视频一区免费| 国产精品野战在线观看| 亚洲av美国av| 久久久久久久久久久丰满| 日本免费一区二区三区高清不卡| 少妇猛男粗大的猛烈进出视频 | 亚洲真实伦在线观看| 插阴视频在线观看视频| 久久久久久久久中文| 国产成人91sexporn| 美女 人体艺术 gogo| 亚洲经典国产精华液单| 免费观看人在逋| 99热这里只有是精品50| 成人国产麻豆网| 在线免费观看不下载黄p国产| 亚洲成人av在线免费| 露出奶头的视频| 午夜久久久久精精品| av天堂在线播放| 美女cb高潮喷水在线观看| 精品少妇黑人巨大在线播放 | 国产私拍福利视频在线观看| 免费黄网站久久成人精品| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| av免费在线看不卡| 亚洲国产精品久久男人天堂| 国产 一区精品| 成人性生交大片免费视频hd| 女同久久另类99精品国产91| 综合色丁香网| 亚洲人成网站在线观看播放| 国产一区亚洲一区在线观看| 欧美国产日韩亚洲一区| av在线观看视频网站免费| 给我免费播放毛片高清在线观看| 看十八女毛片水多多多| 舔av片在线| 真实男女啪啪啪动态图| 亚洲美女视频黄频| 亚洲精品456在线播放app| 我的老师免费观看完整版| 成人午夜高清在线视频| 尾随美女入室| 国产精品三级大全| 久久久久久久亚洲中文字幕| 亚洲成人中文字幕在线播放| 成人鲁丝片一二三区免费| 国产高清视频在线播放一区| 久久久久九九精品影院| 日本熟妇午夜| 午夜亚洲福利在线播放| 日韩在线高清观看一区二区三区| 久久久久久国产a免费观看| 国内久久婷婷六月综合欲色啪| 日韩精品青青久久久久久| 国产亚洲av嫩草精品影院| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站| 日韩大尺度精品在线看网址| av免费在线看不卡| 99久久成人亚洲精品观看| 亚洲精品日韩av片在线观看| 国产日本99.免费观看| 日本-黄色视频高清免费观看| 亚洲av第一区精品v没综合| 精品国内亚洲2022精品成人| 免费不卡的大黄色大毛片视频在线观看 | 大型黄色视频在线免费观看| 在现免费观看毛片| 欧美区成人在线视频| 校园人妻丝袜中文字幕| 亚洲欧美成人精品一区二区| 哪里可以看免费的av片| av在线老鸭窝| 日韩欧美三级三区| 国产乱人偷精品视频| 精品人妻视频免费看| 热99re8久久精品国产| 三级男女做爰猛烈吃奶摸视频| 免费大片18禁| 国产伦一二天堂av在线观看| 噜噜噜噜噜久久久久久91| 亚洲在线自拍视频| 久久久久久久久久黄片| 十八禁国产超污无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 免费看a级黄色片| 国产精品嫩草影院av在线观看| 变态另类成人亚洲欧美熟女| 丝袜美腿在线中文| 国产精品99久久久久久久久| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| av免费在线看不卡| 老女人水多毛片| 婷婷六月久久综合丁香| 中文字幕免费在线视频6| 国产真实伦视频高清在线观看| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 一本久久中文字幕| 亚洲电影在线观看av| 欧美成人a在线观看| 亚洲欧美日韩东京热| 欧美一区二区亚洲| 亚洲国产日韩欧美精品在线观看| 黄色视频,在线免费观看| 欧美一级a爱片免费观看看| 国产精品一及| 国产大屁股一区二区在线视频| 97人妻精品一区二区三区麻豆| a级毛色黄片| 国内少妇人妻偷人精品xxx网站| 亚洲av免费高清在线观看| 欧美区成人在线视频| 美女内射精品一级片tv| 欧美日韩综合久久久久久| 亚洲真实伦在线观看| 日韩 亚洲 欧美在线| 悠悠久久av| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 18禁在线播放成人免费| 欧美人与善性xxx| a级毛片a级免费在线| 天堂av国产一区二区熟女人妻| 深夜精品福利| 免费人成在线观看视频色| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 欧美精品国产亚洲| 国产精品野战在线观看| 狂野欧美激情性xxxx在线观看| 99精品在免费线老司机午夜| 亚洲欧美成人综合另类久久久 | 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 免费av不卡在线播放| 综合色av麻豆| 在线播放国产精品三级| 国产午夜精品久久久久久一区二区三区 | 夜夜爽天天搞| 国产精品人妻久久久久久| 高清日韩中文字幕在线| 色视频www国产| 精品熟女少妇av免费看| 天美传媒精品一区二区| 欧美一级a爱片免费观看看| 男人舔女人下体高潮全视频| 久久午夜亚洲精品久久| 联通29元200g的流量卡| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区成人| 青春草视频在线免费观看| 在线观看av片永久免费下载| 91在线精品国自产拍蜜月| 别揉我奶头 嗯啊视频| 小说图片视频综合网站| 国产黄a三级三级三级人| 最近的中文字幕免费完整| 国产精品国产高清国产av| 哪里可以看免费的av片| 亚洲第一电影网av| 国产精品一二三区在线看| 成人欧美大片| 久久人人精品亚洲av| 久久午夜亚洲精品久久| 国产探花极品一区二区| 国产亚洲91精品色在线| 99热这里只有精品一区| 男女啪啪激烈高潮av片| 亚洲aⅴ乱码一区二区在线播放| 国产精品永久免费网站| 在线观看av片永久免费下载| 国产成人影院久久av| 国产伦一二天堂av在线观看| 国产成年人精品一区二区| 美女高潮的动态| 两个人的视频大全免费| 69人妻影院| av在线老鸭窝| 国产av在哪里看| 久久午夜亚洲精品久久| 亚洲人成网站在线观看播放| 亚洲天堂国产精品一区在线| 99视频精品全部免费 在线| 欧美性猛交╳xxx乱大交人| 99久久成人亚洲精品观看| 一级a爱片免费观看的视频| 麻豆国产av国片精品| 少妇的逼水好多| 婷婷六月久久综合丁香| 亚洲国产精品成人久久小说 | 99热这里只有是精品50| 国产高清视频在线播放一区| 欧美日韩乱码在线| 久久精品国产鲁丝片午夜精品| 成人特级黄色片久久久久久久| av在线天堂中文字幕| 久久草成人影院| 精品一区二区三区视频在线观看免费| 久久精品综合一区二区三区| 老司机福利观看| 99久久久亚洲精品蜜臀av| 久久午夜福利片| 国产亚洲精品综合一区在线观看| 久久精品国产清高在天天线| 色视频www国产| 国产午夜福利久久久久久| 欧美在线一区亚洲| 麻豆一二三区av精品| 免费电影在线观看免费观看| 成年女人毛片免费观看观看9| 欧美一区二区精品小视频在线| 高清毛片免费观看视频网站| 久久久久久久久久成人| 国产精品三级大全| 午夜精品一区二区三区免费看| 欧美精品国产亚洲| 日日摸夜夜添夜夜添av毛片| av视频在线观看入口| 国产黄片美女视频| 欧美成人一区二区免费高清观看| 国产高潮美女av| 午夜激情欧美在线| 少妇被粗大猛烈的视频| 久久久久久久亚洲中文字幕| 噜噜噜噜噜久久久久久91| 久久国产乱子免费精品| 欧美高清性xxxxhd video| 久久久久久久久大av| 国产探花在线观看一区二区| 欧美性感艳星| 波野结衣二区三区在线| 在现免费观看毛片| 日本五十路高清| 老司机午夜福利在线观看视频| 欧美bdsm另类| 亚洲在线观看片| a级毛片a级免费在线| 国产精品乱码一区二三区的特点| 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 久久久久性生活片| 国产不卡一卡二| 亚洲天堂国产精品一区在线| 色综合站精品国产| 久久久久国产精品人妻aⅴ院| 精品久久久久久久末码| ponron亚洲| 99热这里只有是精品在线观看| 在线观看一区二区三区| 国产在视频线在精品| 精品久久国产蜜桃| 少妇熟女aⅴ在线视频| 久久天躁狠狠躁夜夜2o2o| 我的老师免费观看完整版| 99热这里只有是精品在线观看| 亚洲精品成人久久久久久| 国产aⅴ精品一区二区三区波| 精品久久国产蜜桃| 国产精品乱码一区二三区的特点| 国产精华一区二区三区| 国产三级中文精品| 亚洲欧美精品综合久久99| 国产精品免费一区二区三区在线| 在线国产一区二区在线| 女人十人毛片免费观看3o分钟| 国产亚洲精品久久久久久毛片| av女优亚洲男人天堂| 日韩欧美精品v在线| 国产久久久一区二区三区| 深夜精品福利| 国产精品亚洲美女久久久| 国模一区二区三区四区视频| 亚洲激情五月婷婷啪啪| 国产精品久久电影中文字幕| 少妇的逼好多水| 久久精品国产亚洲av涩爱 | 国内久久婷婷六月综合欲色啪| 少妇的逼好多水| 欧美一级a爱片免费观看看| 俄罗斯特黄特色一大片| 九九在线视频观看精品| 国产日本99.免费观看| 久久精品国产亚洲av香蕉五月| 亚洲丝袜综合中文字幕| 欧美zozozo另类| 三级国产精品欧美在线观看| 男人狂女人下面高潮的视频| 亚洲精品色激情综合| 色5月婷婷丁香| 日韩欧美精品免费久久| 精品人妻偷拍中文字幕| 亚洲久久久久久中文字幕| 少妇的逼好多水| 少妇熟女aⅴ在线视频| 久久久色成人| 亚洲精品成人久久久久久| 国产大屁股一区二区在线视频| 日韩强制内射视频| 国产黄片美女视频| 久久精品夜色国产| 久久久成人免费电影| 亚洲中文日韩欧美视频| 亚洲国产精品成人综合色| 亚洲一区二区三区色噜噜| 成人二区视频| 丰满的人妻完整版| 午夜视频国产福利| 少妇猛男粗大的猛烈进出视频 | 乱人视频在线观看| 国产精品亚洲一级av第二区| 欧美色欧美亚洲另类二区| 婷婷六月久久综合丁香| 午夜日韩欧美国产| 日韩高清综合在线| 欧美又色又爽又黄视频| 乱码一卡2卡4卡精品| 成人漫画全彩无遮挡| 在线天堂最新版资源| 国产69精品久久久久777片| 国产精品国产三级国产av玫瑰| 亚洲性久久影院| 国产精华一区二区三区| 久久国产乱子免费精品| 精品国内亚洲2022精品成人| 卡戴珊不雅视频在线播放| 国产高潮美女av| 精品久久久久久成人av| 麻豆国产97在线/欧美| 国产白丝娇喘喷水9色精品| 国产麻豆成人av免费视频| avwww免费| 天堂av国产一区二区熟女人妻| 天天躁日日操中文字幕| 男女啪啪激烈高潮av片| av黄色大香蕉| 日韩欧美免费精品| 三级国产精品欧美在线观看| 精品99又大又爽又粗少妇毛片| 变态另类丝袜制服| 婷婷亚洲欧美| 在线免费观看不下载黄p国产| 99热网站在线观看| 亚洲,欧美,日韩| 99久久精品国产国产毛片| 97人妻精品一区二区三区麻豆| 大又大粗又爽又黄少妇毛片口| 亚洲七黄色美女视频| ponron亚洲| 99久久成人亚洲精品观看| 欧美性猛交╳xxx乱大交人| 精品熟女少妇av免费看| 老司机福利观看| 久久久久久久午夜电影| 国产精品亚洲一级av第二区| 国产精品伦人一区二区| 精品久久国产蜜桃| 在线观看66精品国产| 美女高潮的动态| 长腿黑丝高跟| 少妇高潮的动态图| 少妇熟女欧美另类| 国产精品,欧美在线| 亚洲,欧美,日韩| 国产人妻一区二区三区在| 国内揄拍国产精品人妻在线| 91在线观看av| 好男人在线观看高清免费视频| 日韩精品有码人妻一区| 午夜视频国产福利| 亚洲欧美清纯卡通| 欧美bdsm另类| 日韩在线高清观看一区二区三区| 国内精品宾馆在线| 日韩中字成人| 欧美高清成人免费视频www| 国产亚洲av嫩草精品影院| 久久精品国产自在天天线| 午夜日韩欧美国产| 国内精品宾馆在线| 日本 av在线| 91在线观看av| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看| 中文资源天堂在线| 精品一区二区三区视频在线| 国产精品久久久久久精品电影| 精品少妇黑人巨大在线播放 | 在现免费观看毛片| 久久这里只有精品中国| 国产精品亚洲美女久久久| 人妻少妇偷人精品九色| 色吧在线观看| 少妇熟女aⅴ在线视频| 秋霞在线观看毛片| 国产亚洲91精品色在线| 在线a可以看的网站|