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

    Impact of consolidation pressure on contaminant migration in clay liner

    2013-06-22 13:25:26ZhihongZHANGZhaogangXUXiuliDUHongyanLI
    Water Science and Engineering 2013年3期

    Zhi-hong ZHANG*, Zhao-gang XU, Xiu-li DU, Hong-yan LI

    The Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, P. R. China

    Impact of consolidation pressure on contaminant migration in clay liner

    Zhi-hong ZHANG*, Zhao-gang XU, Xiu-li DU, Hong-yan LI

    The Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, P. R. China

    Consolidation deformation occurs in clay liners under the self-weight of wastes at a simple garbage dump or dredged sediment dump, which leads to a decrease in the porosity. However, the migration of contaminants in clay liners is influenced by the porosity. Thus, the impact of consolidation deformation of clay liners on the migration of contaminants cannot be ignored. Based on Biot’s consolidation theory, the contaminant migration theory, and consideration of the three kinds of migration mechanisms of convection, diffusion, and adsorption, a one-dimensional migration model of contaminants in deforming porous media was established, and the finite difference method was adopted to obtain the numerical solutions for an established initial-boundary value problem. The impact of consolidation pressure on the migration law of a contaminant was studied. The results show that, regardless of adsorption modes, different consolidation pressures have similar impacts on the migration law of the contaminant. Namely, over a certain migration time, the greater the consolidation pressure is, the smaller the migration depth of the contaminant. The results also show that, while the migration time increases, the impact of a certain increment of consolidation pressure on the variation of contaminant concentration with the depth increases gradually and, while the migration depth increases, the impact of a certain increment of consolidation pressure on the variation of the contaminant concentration with time increases gradually.

    clay liner; consolidation pressure; consolidation deformation; contaminant migration

    1 Introduction

    Natural clay layers are widely used as the liners at simple garbage dumps and dredged sediment dumps, without consideration of impermeable treatment or environmental protection. Compared with sanitary landfills, contaminants in leachate or in sediment penetrate through the clay liner more easily, and soil and groundwater may be polluted. Thus, this type of dump poses a high risk of secondary pollution to the surrounding environment. In addition, due to the large dump area, large dumped amount, and high absolute content of waste, as time passes on, if the clay liner is penetrated and leakage occurs, the influence on the surroundingenvironment lasts a long time and spreads over large areas. Therefore, it is of significance to study whether the contaminant will penetrate the liner and cause secondary pollution to the surrounding soil and groundwater over a certain time, which is related to the migration of the contaminants in the liner. When studying the migration of contaminants in a liner, soil is generally assumed to be a rigid porous media, without consideration of the impact of deformation of the liner on the contaminant migration (El-Zein 2008; Li and Cleall 2010), which cannot fully reflect the actual contaminant migration in the liner. On one hand, the contaminant migration in the liner is affected by multiple complex mechanisms, among which the convection and dispersion effects are directly related with soil porosity. On the other hand, due to the facts that the service life of these dumps is generally long and the dumped amount is huge, the clay liner will undergo time-dependent consolidation deformation under the self-weight of garbage or dredged sediment, leading to a decrease in the porosity, which is bound to affect the migration process of contaminants. Hence, when studying the migration of contaminants in the clay liner, the impact of consolidation deformation of the clay liner on the migration process of the contaminants cannot be ignored.

    Studies on contaminant migration in deforming soil have been performed recently by some scholars. Smith (2000) developed a one-dimensional model of contaminant transport in a deforming porous medium based on the conservation of mass including both fluid and soil, and found that the deformation of soil had a great impact on the contaminant transport. Peters and Smith (2002) studied one-dimensional contaminant transport through a deforming porous medium based on the research of Smith (2000). Alshawabkeh et al. (2005) and Alshawabkeh and Rahbar (2006) studied one-dimensional solute transport in a capped sediment under consolidation based on Terzaghi’s consolidation theory, and found that consolidation accelerated the breakthrough of the solute to the overlying cap. Fox (2007a, 2007b) developed a model, coupled with large strain consolidation and solute transport, to study solute transport in soil with large deformation, and found that consolidation had a long-term impact on solute transport. Lee et al. (2009) and Lee and Fox (2009) studied consolidation-induced solute transport considering nonlinear non-equilibrium sorption, and found that consolidation-induced advection has an important contribution to solute transport. Lewis et al. (2009a, 2009b) developed a one-dimensional large deformation model of coupled mechanical consolidation and solute transport to study solute transport in deforming clay liners, and found that consolidation is capable of accelerating solute transport. Zhang et al. (2008) developed a one-dimensional model for contaminant transport combining Biot’s consolidation theory with contaminant transport theory, and an analytical solution of the model was obtained based on simplified conditions. It was found that consolidation had an important impact on contaminant transport. Huang et al. (2010) developed a contaminant transport model considering the influence of liner consolidation on transport coefficients to study the contaminant transport in a deforming liner, and found that the impact of consolidation on contaminant transport should beconsidered when consolidation pressure is high.

    Based on the investigation of Zhang et al. (2008), a one-dimensional contaminant transport model coupling Biot’s consolidation theory with the contaminant migration theory was improved. Without any simplification in parameters of the model, a numerical solution for the model was achieved, and the impact of consolidation pressure on the migration of a contaminant was studied for different adsorption modes.

    2 Model establishment

    2.1 Basic assumptions

    The basic assumptions of the migration model are as follows:

    (1) The soil is in an isothermal condition.

    (2) The soil is an isotropic, saturated, and linear elastic continuum with a small amount of deformation.

    (3) The soil particles and pore water are incompressible.

    (4) The consolidation deformation of soil caused by the self-weight of clay liners is not considered.

    (5) Soil deformation and contaminant migration occur only in the vertical direction.

    2.2 Biot’s one-dimensional consolidation model

    The one-dimensional model of Biot’s consolidation theory, taking no account of body force, can be written as

    whereGis the shear modulus of soil,νis Poisson’s ratio of soil,uis the excess pore water pressure,wis the displacement of soil in thezdirection, andzis the spatial coordinate taking downwards as positive.

    There are two unknown parametersuandw, in the one-dimensional consolidation equation. Therefore, one more equation is needed to solve this problem. Due to the fact that pore water is incompressible, the continuity equation of water flow can be used during the soil consolidation process. Namely, the changing rates of water and the whole soil volume are equal in a representative elementary volume (REV). Thereby, the one-dimensional flow continuity equation can be written as

    whereρwis the density of water;gis the gravitational acceleration;εzis the vertical strain of soil;kis the permeability coefficient of soil, which changes with the consolidation deformation; andtis time.

    Applying the generalized Hooke’s law in Eqs. (1) and (2), one obtains

    whereσzis the total stress of soil in thezdirection.

    An empirical equation (Southeast University et al. 2005) related to the porosity ratio was adopted to describe the change law ofk, which can be defined as

    wherek0is the initial permeability coefficient,eis the void ratio of soil,e0is the initial void ratio of soil, andckis the permeability index and is defined as 0.5e0.

    As the pore water and soil particles are assumed to be incompressible, the changed soil volume is equal to the changed pore volume during the consolidation process. Th us, the computational equations of porosity and the void ratio can be expressed, respectively, as follows:

    wherenis the porosity of soil,n0is the initial porosity, andpis the consolidation pressure on the clay liner.

    2.3 One-dimensional migration model of contaminant in deforming porous media

    While contaminants migrate in soil, a portion of the contaminants are adsorbed on the soil skeleton, and this portion of the contaminants moves together with the soil particles during the soil consolidation compression process; other contaminants move alon g with the pore water. Based on the mass conservation law of contaminants in the liquid phase and in the solid phase, and only taking account of convection, diffusion, and adsorption mechanisms, a one-dimensional migration model of contaminants in deforming porous media was established in the space coordinates, as follows:

    wherecis the concentration of a contaminant in pore water;ρsis the density of soil particles;Sis the adsorption concentration on the soil skeleton, which reflects the adsorption properties of soil particles for the contaminant and can be described by the adsorption equation;vis the Darcy seepage speed;vsis the velocity of the soil skeleton; andDis the hydrodynamic dispersion coefficient, composed of the effective molecular diffusion coefficient and the mechanical diffusion coefficient, and can be obtained from Eqs. (8), (9), and (10):

    whereDeis the effective molecular diffusion coefficient;DLis the mechanical dispersion coefficient;D0is the free diffusion coefficient in an open water body;vtis the actual seepage speed of pore water;αLis the vertical dispersity, whose value is generally in the range of 0.1 cm to 5 m in indoor tests (Wang 2008); andτis the tortuosity factor of the porous media, related to the pore properties of soil, and can be defined by porosity asτ=nm, wheremis an empirical parameter (Liu et al. 2004).

    Compared with the conventional contaminant migration equation, the novelty of Eq. (7) is the introduction of the term describing the movement of contaminants adsorbed by soil particles. Biot’s consolidation equation and the contaminant migration equation are linked by this term and parameters affected by consolidation deformation, and, thus, the consolidation equation and the migration equation are coupled. Because of consolidation deformation, the parameters such as the porosity, the Darcy seepage speed, the velocity of the soil skeleton, and the hydrodynamic dispersion coefficient, are variable.

    2.4 Definite conditions

    To obtain the numerical solution of Eq. (7), the following assumptions were established for the model:

    (2) The soil layer discharges water from the upper surface.

    (3) The concentration of contamination sourcec0is a constant and does not change with time.

    (4) The background contaminant concentration in the soil layer is not taken into consideration.

    The definite conditions of the model can be obtained as follows from the above assumptions:

    For the consolidation equation:

    For the migration equation:whereHis the thickness of the soil layer. Combined with the definite conditions, Eq. (11) through (16), a numerical solution for Eqs. (1), (2), and (7) can be worked out with the finite difference method. The migration law of the contaminant can be obtained with consideration of the small deformation of soil.

    3 Results and discussion

    3.1 Parameter values

    Ammonia nitrogen was selected as the simulated contaminant. The concentration of the contamination source wasc0=500 mg/L and the free diffusion coefficient of ammonia nitrogen in the open water body was=1.76× 10?9m2/s. The studied soil layer is a natural clay layer whose thickness isH= 2 m, initial void ratio ise0=0.8, initial permeability coefficient iscm/s, density of soil particles is, Poisson’s ratio isν=0.3, shear modulus isG=2.6× 103kPa, empirical parameter ism= 2.5, and vertical dispersity isαL= 0.001 m.

    The program Matlab was used to obtain the solution of the migration model, and the effects of the consolidation pressure on the temporal and spatial distribution of the ammonia nitrogen concentration in the clay liner were studied for different adsorption modes.

    3.2 Results of numerical calculation and discussion

    Three adsorption modes were considered to describe the adsorption property of soil particles: the linear adsorption mode, Freundlich adsorption mode, and Langmuir adsorption mode. The equations of different adsorption modes and the values of adsorption parameters are shown in Table 1.

    Table 1 Adsorption modes and main parameters

    3.2.1 Impact of consolidation pressure on variation of contaminant concentration with depth

    As shown in Figs. 1, 2, and 3, regardless of the kinds of adsorption modes, for the specific migration years, the greater the consolidation pressure is, the smaller the migration depth of the contaminant in the soil. It is also shown in the figures that the frontal surface of the contaminant concentration shifts toward the left when the consolidation pressure increases. The contaminant concentration values on the frontal surfaces at a certain depth decrease with the increase of consolidation pressure. At the same time, it can also be observed that along with the increase of the migration time, the contaminant concentration frontal surfacescorresponding to different consolidation pressures become sparse, which shows that the impact of a certain increment of the consolidation pressure on contaminant migration increases gradually with the migration time. In other words, the increment of the corresponding migration depth under a certain increment of consolidation pressure is increasingly large with the migration time. In addition, regardless of the kinds of adsorption isotherm equations applied, the migration depth of the contaminant will increase gradually with the increase of migration time under the effect of a certain consolidation pressure.

    The phenomenon that the migration depth and contaminant concentration at a certain depth decrease as the consolidation pressure increases is due to different compressed volumes of soil under different consolidation pressures, which result in different impact of consolidation pressures on the migration of the contaminant. During the compression process of soil under consolidation pressure, pore water is discharged and porosity decreases, which narrows the migration channel of the contaminant, reduces the permeability, weakens the convection and diffusion effects, and retards the migration of the contaminant. In addition, the greater the consolidation pressure is, the greater the amount of soil compression deformation and reduction of porosity, and the stronger the retardation of the contaminant migration, so the migration depth and contaminant concentration at a certain depth are smaller when the consolidation pressure is greater.

    Fig. 1 Variation of contaminant concentration with depth at different consolidation pressures for linear adsorption

    Fig. 2 Variation of contaminant concentration with depth at different consolidation pressures for Freundlich adsorption

    Fig. 3 Variation of contaminant concentration with depth at different consolidation pressures for Langmuir adsorption

    3.2.2 Impact of consolidation pressure on variation of contaminant concentration with time

    As shown in Figs. 4, 5, and 6, regardless of the kinds of adsorption modes, the impacts of consolidation pressure on the variation of contaminant concentration with time are similar. That is, for any depth profile in the clay liner, the contaminant concentration increases as time goes on under different consolidation pressures, until it finally reaches the value that is equal to the contamination source. The closer to the contamination source the depth profile is, the faster the contaminant concentration accumulates, and the more quickly the contaminant concentration will reach a value that is equal to that of the contamination source. In other words, the contaminant concentration at a depth profile that is far away from the contaminant source accumulates more slowly. However, as long as the migration time is long enough, with a constant concentration of the contamination source and a continuous supply, the concentration value will also eventually reach a value equal to that of the contamination source. For any depth profile studied, at a fixed time, the concentration on the frontal surfaces corresponding to different consolidation pressures decreases as the consolidation pressure increases. At the same time, it can be observed that the contaminant concentration frontal surfaces corresponding to different consolidation pressures gradually become sparse when the migration depth increases, which means that with the increase of the depth, the impact of a certain increment of consolidation pressure on the variation of the contaminant concentrationwith time increases gradually. In addition, no matter what the adsorption mode is, for a certain consolidation pressure, as the depth studied increases, the frontal surface of the contaminant concentration gradually becomes gentle, and the farther the depth profile is away from the contamination source, the more obvious this trend is.

    Fig. 4 Variation of contaminant concentration with time under different consolidation pressures for linear adsorption

    Fig. 5 Variation of contaminant concentration with time under different consolidation pressures for Freundlich adsorption

    Fig. 6 Variation of contaminant concentration with time under different consolidation pressures for Langmuir adsorption

    The reason that the concentration values on the frontal surfaces corresponding to different consolidation pressures at a fixed time decrease as the consolidation pressure increases is that the greater the consolidation pressure is, the greater the amount of deformation of soil and the porosity reduction, and the slower the contaminant migrates, resulting in less contaminant accumulation at a certain depth profile, and thus a lower concentration in the same migration time when the consolidation pressure is higher.

    3.2.3 Contaminant migration depth at different consolidation pressures

    The main purposes of studying the contaminant migration in a clay liner are to judge whether the contaminant would penetrate through the clay liner in a certain migration period, to provide a scientific basis for risk assessment of landfills or dredged sediment dumps, and to carry out a final evaluation to investigate the risk of secondary pollution to the surrounding environment and the groundwater.

    According toQuality Standard for Ground Water(GB/T14848-93) of China, potable water must meet the standard requirements of Class III, in which the concentration of ammonia nitrogen should be no more than 0.2 mg/L. If this value is set as a standard for different adsorption modes, the migration depth of ammonia nitrogen at different times and under different consolidation pressures from Figs 1, 2, and 3 can be determined.

    For both the linear adsorption and Langmuir adsorption, within 20 years, the contaminant concentration at the bottom of the investigated soil layer did not exceed the standard value. The relation curves of the migration depths versus consolidation pressures at different migration times are shown in Fig. 7. However, for the Freundlich adsorption, the contaminant concentration at the bottom of the investigated soil layer at the migration time of 20 years already exceeded the standard value, and the relation curve of the contaminant concentration at the bottom of the soil layer under different consolidation pressures can be seen in Fig. 8.

    As shown in Fig. 7, regardless of the modes of adsorption, the migration depths of contaminants at different migration times decreased when the consolidation pressure increased, and there was a satisfactory linear relation between the migration depth and the consolidation pressure. The relational expressions and correlation coefficients are shown in Table 2. For the Freundlich adsorption, the contaminant concentration at the bottom of the investigated soil layer at the migration time of 20 years already exceeded the standard value, so the relation curve of the contaminant migration depth and consolidation pressure at the migration time of 20 years is not given for this mode.

    As shown in Fig. 8, for the Freundlich adsorption, the contaminant concentration at the bottom of the investigated soil layer decreased while the consolidation pressure increased atthe migration time of 20 years; they were in an exponential relation, where the relational expression can be written asc=35.04e?0.001pand the correlation coefficient wasR=0.998. The contaminant concentration at the bottom of the investigated soil layer under any consolidation pressure at the migration time of 20 years can be obtained through this fitting formula.

    Fig. 7 Relation curves of contaminant migration depth versus consolidation pressure at different times for different adsorptions

    Fig. 8 Relation curve of contaminant concentration at bottom of soil layer versus consolidation pressure at migration time of 20 years for Freundlich adsorption

    Table 2 Fitted formulas of migration depth and consolidation pressure

    4 Conclusions

    Because of the movement of the contaminant adsorbed to soil particles and parameters affected by consolidation deformation, Biot’s consolidation theory and contaminant migration theory were coupled to establish a migration model of contaminants in deforming porous media. The numerical solution for the established model was obtained with the finite difference method. The impact of consolidation pressure on the migration law of the contaminant for different adsorption modes was studied. From this study, the following conclusions can be obtained:

    (1) The migration depths of the contaminant under different consolidation pressures were compared and analyzed. The results show that regardless of the adsorption modes, the greater the consolidation pressure is, the lesser the migration depth and the weaker the penetration ability of the contaminant. The reason is that the greater the consolidation pressure is, the more compacted the soil layer and the weaker the convection diffusion effect, and thus, the lesser the migration depth of the contaminant at the same migration time.

    (2) The consolidation pressure had a direct impact on the accumulation velocity of the contaminant at any depth profile. The greater the consolidation pressure is, the slower the concentration accumulates, which results in a low concentration of contaminants at the same migration time.

    (3) The variation curve shapes of the contaminant concentration with depth under different consolidation pressures are similar, and so are the variation curves of contaminants with time, which shows that there is no intrinsic impact of the variation of consolidation pressure on the migration law of the contaminant, but there is an impact on the migration speed of the contaminant. This impact will increase gradually as migration time goes on.

    (4) Generally, a basically linear relationship is shown between the migration depth of the contaminant and consolidation pressure for all the three adsorption modes.

    Alshawabkeh, A. N., Rahbar, N., and Sheahan, T. 2005. A model for contaminant mass flux in capped sediment under consolidation.Journal of Contaminant Hydrology, 78(3), 147-165. [doi:10.1061/j. jconhyd.2005.05.003]

    Alshawabkeh, A. N., and Rahbar, N. 2006. Parametric study of one-dimensional solute transport in deformable porous media.Journal of Geotechnical and Geoenvironmental Engineering, 132(8), 1001-1010. [doi:10.1061/(ASCE)1090-0241(2006)132:8(1001)]

    El-Zein, A. 2008. A general approach to the modelling of contaminant transport through composite landfill liners with intact or leaking geomembranes.International Journal for Numerical and Analytical Methods in Geomechanics, 32(3), 265-287. [doi:10.1002/nag.626]

    Fox, P. J. 2007a. Coupled large strain consolidation and solute transport, I: Model development.Journal of Geotechnical and Geoenvironmental Engineering, 133(1), 3-15. [doi:10.1061/(ASCE)1090-0241(2007) 133:1(3)]

    Fox, P. J. 2007b. Coupled large strain consolidation and solute transport, II: Model verification and simulation results.Journal of Geotechnical and Geoenvironmental Engineering, 133(1), 16-29. [doi:10.1061/ (ASCE)1090-0241(2007)133:1(16)]

    Huang, L., Zhao, C. G., Zhang, X. D., and Liu, Y. 2010. 3D contaminant migration model with consolidation dependent transport coefficients.Chinese Journal of Geotechnical Engineering, 32(3), 420-427. (in Chinese)

    Lee, J., and Fox, P. J. 2009. Investigation of consolidation-induced solute transport, II: Experimental and numerical results.Journal of Geotechnical and Geoenvironmental Engineering, 135(9), 1239-1253. [doi: 10.1061/(ASCE)GT.1943-5606.0000048]

    Lee, J., Fox, P. J., and Lenhart, J. J. 2009. Investigation of consolidation-induced solute transport, I: Effect of consolidation on transport parameters.Journal of Geotechnical and Geoenvironmental Engineering, 135(9), 1228-1238. [doi:10.1061/(ASCE)GT.1943-5606.0000047]

    Lewis, T. W., Pivonka, P., Fityus, S. G., and Smith, D. W. 2009a. Parametric sensitivity analysis of coupled mechanical consolidation and contaminant transport through clay barriers.Computers and Geotechnics, 36(1-2), 31-40. [doi:10.1061/j.compgeo.2008.04.003]

    Lewis, T. W., Pivonka, P., and Smith, D. W. 2009b. Theoretical investigation of the effects of consolidation on contaminant transport through clay barriers.International Journal for Numerical and Analytical Methods in Geomechanics, 33(1), 95-116. [doi:10.1002/nag.708]

    Li, Y. C., and Cleall, P. J. 2010. Analytical solutions for contaminant diffusion in double-layered porous media.Journal of Geotechnical and Geoenvironmental Engineering, 136(11), 1542-1554. [doi:10.1002/nag.903]

    Liu, J. G., Wang, H. T., and Nie, Y. F. 2004. Fractal model for predicting effective diffusion coefficient of solute in porous media.Advances in Water Science, 15(4), 458-462. (in Chinese)

    Peters, G. P., and Smith, D. W. 2002. Solute transport through a deforming porous medium.International Journal for Numerical and Analytical Methods in Geomechanics, 26(7), 683-717. [doi:10.1002/nag.219]

    Smith, D. W. 2000. One-dimensional contaminant transport through a deforming porous medium: Theory and a solution for a quasi-steady-state problem.International Journal for Numerical and Analytical Methods in Geomechanics, 24(8), 693-722. [doi:10.1002/1096-9853(200007)24:8<693::AID-NAG91>3.0.CO;2-E]

    Southeast University, Zhejiang University, Hunan University, and Suzhou University of Science and Technology. 2005.Soil Mechanics. Beijing: China Architecture and Building Press. (in Chinese)

    Wang, H. T. 2008.Dynamics of Fluid Flow and Contaminant Transport in Porous Media. Beijing: Higher Education Press. (in Chinese)

    Zhang, Z. H., Li T., Zhao, C. G., Yu, Y., and Tang, B. R. 2008. Model of contamination transport considering consolidation deformation of soils.Rock and Soil Mechanics, 29(6), 1435-1439. (in Chinese)

    (Edited by Yun-li YU)

    This work was supported by the National Natural Science Foundation of China (Grant No. 51008008).

    *Corresponding author (e-mail:zhangzh2002@bjut.edu.cn)

    May 21, 2012; accepted Oct. 12, 2012

    一a级毛片在线观看| 国内精品宾馆在线| 免费看av在线观看网站| 日韩欧美在线乱码| 国产午夜福利久久久久久| 久久国产乱子免费精品| 乱人视频在线观看| 99精品久久久久人妻精品| 老女人水多毛片| 亚洲av二区三区四区| 亚洲成人精品中文字幕电影| 偷拍熟女少妇极品色| 国产真实乱freesex| 国产在视频线在精品| 国产成人一区二区在线| 精品人妻1区二区| 麻豆国产97在线/欧美| 午夜福利在线观看吧| 亚洲av电影不卡..在线观看| 日本熟妇午夜| 久久国产乱子免费精品| 一夜夜www| 韩国av在线不卡| 国产av麻豆久久久久久久| 变态另类成人亚洲欧美熟女| 搡女人真爽免费视频火全软件 | 一个人观看的视频www高清免费观看| 国产69精品久久久久777片| 国产在线男女| 成人毛片a级毛片在线播放| 免费不卡的大黄色大毛片视频在线观看 | 97超视频在线观看视频| 成人无遮挡网站| 内地一区二区视频在线| 欧美激情久久久久久爽电影| 亚洲成人中文字幕在线播放| 一夜夜www| 男人狂女人下面高潮的视频| 特级一级黄色大片| 国产真实伦视频高清在线观看 | 麻豆成人av在线观看| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩卡通动漫| 色av中文字幕| 成人国产综合亚洲| 精品国产三级普通话版| 久久亚洲真实| 欧美高清性xxxxhd video| 少妇高潮的动态图| 男女下面进入的视频免费午夜| 日本爱情动作片www.在线观看 | 国内精品一区二区在线观看| eeuss影院久久| 天美传媒精品一区二区| 中出人妻视频一区二区| 色5月婷婷丁香| 久久人人精品亚洲av| 国产精品精品国产色婷婷| 中文字幕精品亚洲无线码一区| 精品久久久久久久人妻蜜臀av| 12—13女人毛片做爰片一| 欧美精品国产亚洲| 简卡轻食公司| 村上凉子中文字幕在线| 国产主播在线观看一区二区| 亚洲美女黄片视频| 午夜视频国产福利| 亚洲午夜理论影院| 老司机深夜福利视频在线观看| 亚洲专区国产一区二区| 国产精品综合久久久久久久免费| 国产乱人视频| 1000部很黄的大片| 亚洲真实伦在线观看| 国产主播在线观看一区二区| 国产一区二区激情短视频| 18禁在线播放成人免费| 亚洲av第一区精品v没综合| 嫩草影院新地址| 亚洲电影在线观看av| 国产精品综合久久久久久久免费| 国产精品一区二区性色av| 精华霜和精华液先用哪个| 久久亚洲精品不卡| 琪琪午夜伦伦电影理论片6080| 性插视频无遮挡在线免费观看| 亚洲精品456在线播放app | 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 黄色女人牲交| 国产精品爽爽va在线观看网站| 午夜老司机福利剧场| 麻豆国产av国片精品| 啦啦啦观看免费观看视频高清| 免费人成在线观看视频色| 成人美女网站在线观看视频| 亚洲无线观看免费| 亚洲第一区二区三区不卡| 国产爱豆传媒在线观看| 亚洲成人久久性| 亚洲av.av天堂| 国产一区二区三区在线臀色熟女| a在线观看视频网站| 国产真实乱freesex| 精品无人区乱码1区二区| 色尼玛亚洲综合影院| 窝窝影院91人妻| 黄色欧美视频在线观看| 国产亚洲精品综合一区在线观看| 亚洲国产精品合色在线| 999久久久精品免费观看国产| 亚洲精品亚洲一区二区| 日本五十路高清| 少妇被粗大猛烈的视频| 美女高潮喷水抽搐中文字幕| 极品教师在线免费播放| 一区二区三区免费毛片| 国产成人av教育| 精品久久久久久久久久久久久| 成人精品一区二区免费| 成人国产综合亚洲| 男女啪啪激烈高潮av片| 搡女人真爽免费视频火全软件 | 嫩草影视91久久| 免费av毛片视频| 女生性感内裤真人,穿戴方法视频| 一级毛片久久久久久久久女| 国产成人aa在线观看| 久久国产精品人妻蜜桃| 亚洲欧美日韩东京热| 悠悠久久av| 久久精品国产自在天天线| 成年女人永久免费观看视频| 精品人妻熟女av久视频| 在线免费十八禁| 女的被弄到高潮叫床怎么办 | 91麻豆精品激情在线观看国产| 校园人妻丝袜中文字幕| 久久久国产成人免费| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久久久| 日韩欧美国产一区二区入口| 国产一区二区在线观看日韩| 五月玫瑰六月丁香| 狂野欧美激情性xxxx在线观看| 非洲黑人性xxxx精品又粗又长| 观看免费一级毛片| 国产伦一二天堂av在线观看| 亚洲中文字幕日韩| 深夜a级毛片| 哪里可以看免费的av片| 国产精品久久久久久av不卡| 免费看av在线观看网站| 十八禁国产超污无遮挡网站| 国产伦在线观看视频一区| 亚洲一区二区三区色噜噜| 欧美bdsm另类| 国产成年人精品一区二区| 俄罗斯特黄特色一大片| 久久久国产成人免费| 日本在线视频免费播放| 中文资源天堂在线| 成年女人看的毛片在线观看| 亚洲精品一区av在线观看| 欧美精品国产亚洲| 午夜免费激情av| 亚洲av电影不卡..在线观看| 国产精品不卡视频一区二区| 动漫黄色视频在线观看| 国产精品一及| 日本欧美国产在线视频| av在线亚洲专区| 草草在线视频免费看| 国产淫片久久久久久久久| 免费在线观看影片大全网站| 中文字幕熟女人妻在线| 国产久久久一区二区三区| 91久久精品电影网| 搡老妇女老女人老熟妇| 国产大屁股一区二区在线视频| 日韩人妻高清精品专区| 亚洲黑人精品在线| 日日夜夜操网爽| 五月伊人婷婷丁香| 国产真实乱freesex| 不卡视频在线观看欧美| 国产精品人妻久久久影院| 日本撒尿小便嘘嘘汇集6| 网址你懂的国产日韩在线| 国产男人的电影天堂91| 如何舔出高潮| 欧美最新免费一区二区三区| 久久久久久大精品| 在线播放国产精品三级| 日日啪夜夜撸| 一边摸一边抽搐一进一小说| 日日摸夜夜添夜夜添av毛片 | 国产精品电影一区二区三区| 男人舔女人下体高潮全视频| 国产淫片久久久久久久久| 亚洲欧美日韩东京热| 国产精品一区二区三区四区久久| 欧美中文日本在线观看视频| 国内精品美女久久久久久| 制服丝袜大香蕉在线| 国产又黄又爽又无遮挡在线| 成人三级黄色视频| 99在线视频只有这里精品首页| 亚洲欧美日韩东京热| 国产亚洲91精品色在线| 中文字幕免费在线视频6| 悠悠久久av| av视频在线观看入口| 成人综合一区亚洲| 久久精品国产自在天天线| 国产人妻一区二区三区在| 亚洲欧美日韩东京热| АⅤ资源中文在线天堂| 别揉我奶头 嗯啊视频| 国内少妇人妻偷人精品xxx网站| 女人十人毛片免费观看3o分钟| 成人一区二区视频在线观看| 久久这里只有精品中国| 桃色一区二区三区在线观看| 999久久久精品免费观看国产| 日本爱情动作片www.在线观看 | 午夜福利高清视频| videossex国产| 国产精品人妻久久久久久| 久久精品国产清高在天天线| 成人二区视频| 亚洲专区中文字幕在线| 最后的刺客免费高清国语| 在线国产一区二区在线| 99热只有精品国产| 一卡2卡三卡四卡精品乱码亚洲| 日韩人妻高清精品专区| 亚洲欧美清纯卡通| 欧美激情久久久久久爽电影| 嫩草影院精品99| 午夜久久久久精精品| 国产av在哪里看| 国产亚洲av嫩草精品影院| 日本熟妇午夜| 噜噜噜噜噜久久久久久91| 午夜久久久久精精品| 精品久久久久久,| 简卡轻食公司| 麻豆成人av在线观看| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩高清在线视频| 国产 一区精品| 一个人观看的视频www高清免费观看| 级片在线观看| 成人国产一区最新在线观看| 制服丝袜大香蕉在线| 亚洲精品乱码久久久v下载方式| 91久久精品电影网| 国产成人影院久久av| 欧美激情国产日韩精品一区| 有码 亚洲区| 尤物成人国产欧美一区二区三区| 久久国产乱子免费精品| 一级黄片播放器| 全区人妻精品视频| 亚洲aⅴ乱码一区二区在线播放| 免费一级毛片在线播放高清视频| 99久久中文字幕三级久久日本| 极品教师在线视频| 免费看a级黄色片| 精品欧美国产一区二区三| 亚洲av五月六月丁香网| 午夜亚洲福利在线播放| 波多野结衣巨乳人妻| 97超视频在线观看视频| 久久精品国产鲁丝片午夜精品 | 精品人妻一区二区三区麻豆 | 亚洲av美国av| 别揉我奶头 嗯啊视频| 精品久久久久久,| 亚洲成人久久爱视频| 亚洲中文字幕日韩| 一区二区三区免费毛片| 伦理电影大哥的女人| or卡值多少钱| 99热精品在线国产| 丰满人妻一区二区三区视频av| 亚洲av电影不卡..在线观看| 精品人妻视频免费看| 久99久视频精品免费| 国产视频一区二区在线看| 精品免费久久久久久久清纯| 亚洲无线观看免费| 午夜精品在线福利| 天天一区二区日本电影三级| 日韩强制内射视频| 波多野结衣巨乳人妻| 国产高清视频在线观看网站| 国产高清三级在线| 日日撸夜夜添| 精品久久久久久成人av| 18禁黄网站禁片免费观看直播| 尤物成人国产欧美一区二区三区| 婷婷色综合大香蕉| 啦啦啦啦在线视频资源| 久久欧美精品欧美久久欧美| 别揉我奶头~嗯~啊~动态视频| 又紧又爽又黄一区二区| 中亚洲国语对白在线视频| 久久精品国产亚洲av涩爱 | 天堂动漫精品| 亚洲人成网站高清观看| 亚洲一级一片aⅴ在线观看| 亚洲av电影不卡..在线观看| 精华霜和精华液先用哪个| 97超级碰碰碰精品色视频在线观看| 欧美一区二区精品小视频在线| 禁无遮挡网站| 亚洲国产高清在线一区二区三| 夜夜爽天天搞| 看十八女毛片水多多多| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 少妇裸体淫交视频免费看高清| 看黄色毛片网站| 少妇的逼好多水| 精品一区二区三区av网在线观看| 国产成年人精品一区二区| 精品久久久久久久久久久久久| 一a级毛片在线观看| 女人被狂操c到高潮| 神马国产精品三级电影在线观看| bbb黄色大片| 1000部很黄的大片| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 美女高潮喷水抽搐中文字幕| 国产成人a区在线观看| 一进一出好大好爽视频| 一本精品99久久精品77| 免费av毛片视频| 欧美高清性xxxxhd video| 美女免费视频网站| 老司机福利观看| 日本与韩国留学比较| 国产欧美日韩一区二区精品| 三级国产精品欧美在线观看| 麻豆国产av国片精品| 亚洲中文字幕一区二区三区有码在线看| videossex国产| 国产精品98久久久久久宅男小说| 亚洲aⅴ乱码一区二区在线播放| 两性午夜刺激爽爽歪歪视频在线观看| av在线天堂中文字幕| 精品久久久久久久人妻蜜臀av| 少妇的逼好多水| 国产探花在线观看一区二区| АⅤ资源中文在线天堂| av专区在线播放| АⅤ资源中文在线天堂| 国产色婷婷99| 女生性感内裤真人,穿戴方法视频| 欧美最黄视频在线播放免费| 久久精品综合一区二区三区| 亚洲精华国产精华液的使用体验 | 国产精品人妻久久久影院| 男女那种视频在线观看| 免费高清视频大片| 天堂av国产一区二区熟女人妻| 国产av不卡久久| 午夜精品一区二区三区免费看| 日本精品一区二区三区蜜桃| 99精品在免费线老司机午夜| 久久久久久大精品| 国产精品亚洲美女久久久| 偷拍熟女少妇极品色| 欧美人与善性xxx| 一个人免费在线观看电影| 色5月婷婷丁香| 我要搜黄色片| 婷婷亚洲欧美| 中亚洲国语对白在线视频| 免费av不卡在线播放| www.www免费av| 在线天堂最新版资源| 男人的好看免费观看在线视频| 国产视频一区二区在线看| 国产成人影院久久av| 欧美三级亚洲精品| 精品人妻熟女av久视频| 亚洲精品一区av在线观看| 欧美3d第一页| 久久国产乱子免费精品| 日本三级黄在线观看| 久久婷婷人人爽人人干人人爱| 又爽又黄a免费视频| 日韩欧美三级三区| 国产探花在线观看一区二区| 亚洲狠狠婷婷综合久久图片| 日韩国内少妇激情av| 久久6这里有精品| 一个人看的www免费观看视频| 欧美成人一区二区免费高清观看| 不卡一级毛片| 一进一出抽搐动态| 日韩欧美精品v在线| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 中文在线观看免费www的网站| 高清日韩中文字幕在线| 联通29元200g的流量卡| 淫秽高清视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 国产黄片美女视频| 日本一二三区视频观看| 国产女主播在线喷水免费视频网站 | 少妇人妻精品综合一区二区 | 欧美绝顶高潮抽搐喷水| 欧美成人a在线观看| 熟女电影av网| 国产精品一区www在线观看 | 国产av一区在线观看免费| 日韩欧美三级三区| 波多野结衣巨乳人妻| 欧美一区二区亚洲| 免费无遮挡裸体视频| 99riav亚洲国产免费| 麻豆久久精品国产亚洲av| 狂野欧美白嫩少妇大欣赏| 国产三级在线视频| 亚洲avbb在线观看| 精品乱码久久久久久99久播| 波多野结衣巨乳人妻| 我要搜黄色片| 在线a可以看的网站| 国产精品嫩草影院av在线观看 | 一本一本综合久久| 午夜免费成人在线视频| 精品久久久久久成人av| 国产精品一区二区三区四区久久| 亚洲成av人片在线播放无| 日本免费a在线| 观看免费一级毛片| 色哟哟·www| 欧美日本视频| 久久久久久久精品吃奶| 午夜福利在线观看免费完整高清在 | 亚洲精品粉嫩美女一区| 国产亚洲精品综合一区在线观看| 91在线精品国自产拍蜜月| 国产高清不卡午夜福利| 久99久视频精品免费| 村上凉子中文字幕在线| 直男gayav资源| 精品午夜福利在线看| 日韩欧美在线乱码| 99国产精品一区二区蜜桃av| 日韩中文字幕欧美一区二区| 三级毛片av免费| 亚洲欧美日韩无卡精品| 日韩欧美国产在线观看| 久久久久久国产a免费观看| 欧美激情国产日韩精品一区| 午夜老司机福利剧场| 欧美一区二区精品小视频在线| 日本黄色片子视频| 99热网站在线观看| 小蜜桃在线观看免费完整版高清| 91在线精品国自产拍蜜月| 久久香蕉精品热| 18禁黄网站禁片午夜丰满| 日韩欧美精品v在线| 国产av一区在线观看免费| 九色成人免费人妻av| 欧美日本视频| 一边摸一边抽搐一进一小说| 中文字幕av成人在线电影| 女人被狂操c到高潮| 亚洲综合色惰| 真人做人爱边吃奶动态| 亚洲欧美日韩卡通动漫| 精品一区二区三区人妻视频| 美女高潮的动态| 他把我摸到了高潮在线观看| 国产在线男女| 毛片一级片免费看久久久久 | 亚洲精品久久国产高清桃花| ponron亚洲| 美女高潮的动态| 免费av毛片视频| 精品久久久久久久人妻蜜臀av| а√天堂www在线а√下载| 夜夜看夜夜爽夜夜摸| 热99re8久久精品国产| 美女免费视频网站| 精品久久久噜噜| 免费看美女性在线毛片视频| 国产亚洲精品久久久久久毛片| 欧美xxxx性猛交bbbb| 欧美bdsm另类| 欧美潮喷喷水| 女人十人毛片免费观看3o分钟| 永久网站在线| 又爽又黄无遮挡网站| 国产精品免费一区二区三区在线| 中文字幕av在线有码专区| 99国产精品一区二区蜜桃av| 国产高清视频在线播放一区| 成熟少妇高潮喷水视频| 亚洲熟妇熟女久久| 色5月婷婷丁香| 国产成人影院久久av| 欧美又色又爽又黄视频| 在线播放国产精品三级| 国产一区二区三区在线臀色熟女| 麻豆成人av在线观看| 99热6这里只有精品| 国产精品av视频在线免费观看| 国产精品野战在线观看| 免费搜索国产男女视频| 男女啪啪激烈高潮av片| 直男gayav资源| 欧美性猛交黑人性爽| 国产黄片美女视频| 午夜激情福利司机影院| 精品一区二区免费观看| 日韩欧美在线乱码| 午夜老司机福利剧场| 婷婷精品国产亚洲av在线| 黄色丝袜av网址大全| 久久久精品欧美日韩精品| 精品不卡国产一区二区三区| 久久99热6这里只有精品| 99久久中文字幕三级久久日本| 韩国av一区二区三区四区| 久久久久久伊人网av| 亚洲中文日韩欧美视频| 亚洲精品影视一区二区三区av| 国内少妇人妻偷人精品xxx网站| 亚洲人成网站在线播| 国产黄a三级三级三级人| 亚洲国产精品合色在线| 久99久视频精品免费| 国产主播在线观看一区二区| 中文字幕高清在线视频| 国模一区二区三区四区视频| 自拍偷自拍亚洲精品老妇| 久久人妻av系列| 观看美女的网站| 我的老师免费观看完整版| 成熟少妇高潮喷水视频| 国国产精品蜜臀av免费| 国产乱人伦免费视频| 热99re8久久精品国产| 色综合亚洲欧美另类图片| 91在线观看av| 国产亚洲精品久久久com| 搡老岳熟女国产| 窝窝影院91人妻| 嫩草影院入口| 免费看a级黄色片| 亚洲av成人精品一区久久| 一夜夜www| 日本黄大片高清| 精品99又大又爽又粗少妇毛片 | 网址你懂的国产日韩在线| 国产中年淑女户外野战色| 亚洲综合色惰| 熟妇人妻久久中文字幕3abv| 欧美色视频一区免费| 国产不卡一卡二| 亚洲中文日韩欧美视频| 午夜免费激情av| 97人妻精品一区二区三区麻豆| 日韩大尺度精品在线看网址| 欧美黑人巨大hd| 免费看a级黄色片| 午夜精品一区二区三区免费看| 久久久久久久精品吃奶| 男人舔奶头视频| 国产色爽女视频免费观看| 免费大片18禁| 永久网站在线| 亚洲久久久久久中文字幕| 久久久久久大精品| 动漫黄色视频在线观看| 欧美极品一区二区三区四区| 日韩强制内射视频| 久久精品91蜜桃| 91麻豆av在线| 淫秽高清视频在线观看| 日本-黄色视频高清免费观看| 久久草成人影院| 色综合站精品国产| 夜夜看夜夜爽夜夜摸| 亚洲av不卡在线观看| 狂野欧美激情性xxxx在线观看| 动漫黄色视频在线观看| 精品久久久久久成人av| 哪里可以看免费的av片| 一级黄色大片毛片| 男人狂女人下面高潮的视频| 久久久久久国产a免费观看| 国产亚洲精品久久久久久毛片| 极品教师在线视频| 综合色av麻豆| 亚洲中文字幕一区二区三区有码在线看| 中文字幕免费在线视频6| 我要看日韩黄色一级片| 三级国产精品欧美在线观看| 欧美国产日韩亚洲一区| 天美传媒精品一区二区| 欧美+亚洲+日韩+国产| 麻豆一二三区av精品|