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

    ANISOTROPIC PERMEABILITY EVOLUTION MODEL OF ROCK IN THE PROCESS OF DEFORMATION AND FAILURE*

    2012-05-11 06:54:42WANGHuanling

    WANG Huan-ling

    Key Laboratory of Coastal Disaster and Defence, Ministry of Education, Hohai University, Nanjing 210098, China, E-mail: whl_hm@163.com

    CHU Wei-jiang

    Hydrochina Huadong Engineering Corporation, Hangzhou 310014, China

    HE Miao

    Jiangsu Transportation Research Institute CO., LTD, Nanjing 211112, China

    ANISOTROPIC PERMEABILITY EVOLUTION MODEL OF ROCK IN THE PROCESS OF DEFORMATION AND FAILURE*

    WANG Huan-ling

    Key Laboratory of Coastal Disaster and Defence, Ministry of Education, Hohai University, Nanjing 210098, China, E-mail: whl_hm@163.com

    CHU Wei-jiang

    Hydrochina Huadong Engineering Corporation, Hangzhou 310014, China

    HE Miao

    Jiangsu Transportation Research Institute CO., LTD, Nanjing 211112, China

    (Received April 29, 2011, Revised September 30, 2011)

    The rock permeability is an important parameter in the studies of seepage and stress coupling. The micro-cracks and pores can initiate and grow on a small scale and coalesce to form large-scale fractures and faults under compressive stresses, which would change the hydraulic conductivity of the rock, and therefore, the rock permeability. The rock permeability is, therefore, closely related with the micro-cracking growing, coalescence, and macro new fracture formation. This article proposes a conceptual model of rock permeability evolution and a micro kinematics mechanism of micro-cracking on the basis of the basic theory of micromechanics. The applicability of the established model is verified through numerical simulations of in situ tests and laboratory tests. The simulation results show that the model can accurately forecast the peak permeability evolution of brittle rock, and can well describe the macro-experimental phenomenon before the peak permeability evolution of brittle rock on a macro-scale.

    rock hydraulics, micro-mechanics, permeability, model validation

    Introduction

    The stresses will be re-distributed in a rock mass, under external forces or other engineering situations, which in turn will significantly change the permeability within the rock mass. This process is now a hot issue in the study of the coupling between the seepage field and the stress field in rock masses. With the development of new methods and new equipment in rock mechanics tests, especially, the application of the seepage coupling experimental method in the stressstrain process, it is widely recognized that the permeability of rock is closely related with the micro-crack growth and coalescence, and the macro-fracture formation[1]. A large number of experimental curves related with stress-strain relations and permeability evolutions[2]show that the rock stress state and the deformation degree in different directions in the stressstrain process are not the same, therefore, the evolution of micro-crack opening, density, and connectivity are different and the permeability is not the same in different directions, that is, the permeability shows a significant anisotropy[3]. So the rock deformation not only changes the permeability value, but also brings about a significant anisotropy, which should be considered in rock engineering.

    Extensive researches were carried out on the impact of the load (the stress state) on the rock permeability, which may be summarized as follows: (1) The empirical formulas obtained through extensive experiments. Early studies focused on establishing the relationship between seepage and stress by using direct experiments or indirect analyses. Among them, the most important result is on the effect of the normalstress on the permeability of a single fracture. The representative formulas were obtained by Barton (1995) and Oda et al.[4], and the Barton’s formula is widely used at present. The studies of the impact of the shear stress on the permeability of a single fracture lagged behind and with fewer significant results. But it was confirmed by many studies that the shear stress has a significant effect on the rock permeability. The most influential results were obtained by Olsson and Barton[5], Lee and Cho[6], Makurat[7]. With the development of technology in the rock mechanics tests, the seepage characteristics of the natural fracture rock under the three-dimensional stress state were studied by Tang et al.[8], Tan et al.[9], Liu et al.[10]. Useful conclusions were obtained, which may serve as a guide for practical projects. (2) The theoretical formulas based on the conceptual model. Liu et al.[10]derived the permeability formula for fractured rock mass under complex conditions based on the relationships between the stress-seepage and the strain-seepage. In these studies, a common problem is a strong hypothesis towards the anisotropy of permeability. For instance, Liu et al.[10]assumed that the principal directions of seepage and stress were the same, which, strictly speaking, is unsubstantiated, and thus the practicability of the established model is limited. (3) Statistical laws based on a large number of numerical experiments. Many numerical experiments were carried out by using Discrete Fracture Network (DFN) and Discrete Element Model (DEM). The relationship between permeability and stress state, and loading path was obtained by statistical analyses of calculation results. This method is the most accurate in its consideration of anisotropy as compared with other methods, but so far they are limited to two-dimensional cases.

    Fig.1 Conceptual model of permeability evolution caused by microscopic micro-cracks changes in the process of deformation of brittle rock

    If the permeability is not adequately considered in the coupling analysis of the seepage field and the stress field, the practicability of the established coupled model will be limited and the true coupling rules of rock can not be revealed. This study focuses on the anisotropy of the permeability evolution in the process of rock deformation from a micromechanics view, to develop a reasonable mathematical model of fluidsolid coupling. The applicability of the proposed anisotropy permeability model is validated by simulating the in situ experiments and the laboratory experiments.

    1. Anisotropy permeability model

    1.1 Conceptual model

    Extensive rock micro-experiments and theoretical analyses[11,12]show that the micro-evolution of rock micro-cracks in the process of deformation is the main influencing factor for the increase of permeability within the rock system, and micro-cracks are the main channels of water flow. Based on a great number of test results of the stress-strain curves and the permeability curves under the triaxial compression[13,14], a conceptual model of the permeability evolution in different stress and strain processes is established, and the permeability evolution involves the following four phases (Fig.1):

    (1) Elastic compression stage. The compression force tends to close the initial micro-cracks and the deformation of micro-cracks is reversible at this stage, that is, the deformation of micro-cracks disappears after unloading, and the permeability decreases at this stage.

    (2) Stable stage of compression. The initial opening of micro-cracks is the minimum at this stage.The micro-cracks do not grow or slip and the permeability at this stage changes very little and can be approximately considered as a constant.

    (3) Rapid increase stage for permeability. The micro-cracks slip and expand with the increase of deformation, the fracture opening increases, and the fracture network connectivity is enhanced at this stage. Meanwhile, the new fissures within the rock system begin to crop up. These factors increase the water conductivity of rock by 2-4 orders of magnitude.

    (4) Post-peak phase. The material becomes softened at this stage, and the rock permeability is changed closely according to the confining pressure during the test. When the confining pressure is high, the formed macro-cracks have a tendency of compaction, thus the permeability decreases, when the confining pressure is low, the permeability further increases with the increase of rock deformation, but the increase rate is very small.

    In order to study the macroscopic phenomenon of anisotropic permeability changes in the stress and strain processes using the established model, one must establish the micro-mechanical properties of microcracks based on the following rules and assumptions:

    (1) The change of the stress (deformation) state will change the opening, density, and connectivity of micro-cracks, and thus will change the rock permeability, that is, the change in rock permeability is related with the opening, density, and connectivity of microcracks.

    (2) The water flow within the micro-cracks obeys the generalized flow law.

    (3) The critical starting gradient is not considered.

    (4) The rock is regarded as a binary system composed of the rock matrix and micro-cracks. The rock permeability depends on the following two factors: the single fracture permeability, and the connectivity degree of the micro-cracks.

    (5) The rock mass has a Representative Element Volume (REV).

    (6) The shape of the rock fracture is assumed to follow the Baecher disk model.

    (7) The interaction between micro-cracks is neglected, and the mechanical effect of all micro-cracks is in line with the superposition principle.

    1.2 Permeability of a single micro-crack

    The generalized flow law can be expressed as follows when the hydraulic opening of a single fracture is more than 50 μm

    where Q is the flow in unit time at the cross-section of micro-cracks, HΔ is the hydraulic gradient,wρ is the density of water, g is the gravity acceleration, μ is the viscosity of water, b is the mechanical opening of micro-cracks, ξ is the index of mechanical opening of micro-cracks, and f is the roughness of the surface of micro-cracks.

    It can be seen from Eq.(1) that the generalized flow law would be turned into a cubic law when ξ is equal to three. Equation (1) can be expressed in the form of Darcy’s law as

    Equation (2) can be expressed in the tensor form as

    whereijδ is the unit vector, andin andjn are normal unit vectors of fracture surfaces.

    Therefore, the permeability coefficient K of the single micro-cracks is expressed as

    The relationship between the permeability coefficient and the permeability is as follows

    Thus, the permeability k of a single microcrack can be written in the following form

    By comparing Eq.(4) with Eq.(6), it can be seen that the permeability k is only related to the transport characteristics of fracture networks, while the permeability coefficient K is related to fluid properties in addition to the fracture network.

    The macro-permeability of a single group microcrack can be obtained by the homogenization method:according to assumption (6), and for describe the geometric properties of a group of micro-cracks, the following three variables are required, that is,,, and aτ, whereis the normal vector of the ith group microcrack (τ=1,n), and Ντand aτdenote the numbers and radii of micro-cracks in unit volume, respectively. Here, the fracture density can be expressed as dτ=Nτaτ3. The volume integral is defined as follows for the flow velocity within all micro-cracks in REV, in order to obtain the average velocity within a group of micro-cracks

    where viis the actual velocity vector of water within the rock mass,is the actual velocity vector of water flow along the fracture within the rock mass, and Ωcis the volume occupied by fractures.

    As for the group of micro-cracks with the unit normal vector ofn, the average openings of microcracks in the same direction in the REV are assumed to have the same value, and thus Ωccan be written as

    Equations (3) and (8) are substituted into Eq.(7), which is then integrated on unit spheres, to obtain the following equation

    Therefore, the macro-permeability and permeability coefficients of a group of micro-cracks can be expressed as:

    1.3 Permeability of groups of micro-cracks

    Generally, the fractures can be classified according to whether they are connected or disconnected as far as the permeable fractures containing water are concerned. The connected fractures would form a water flow channel, and disconnected fractures are not involved in the water flow and the cycle alternation and need not be considered in the seepage calculation. Therefore, it is necessary to introduce a parameter to describe the overall connectivity of micro fracture networks in the permeability tensor.

    We assume that λ is the fracture connectivity factor with values in the range [0,1]. =0λ indicates that various micro-cracks are isolated and they cannot form a water flow network, while =1λ indicates that the micro-cracks are fully developed, and each micro-crack is involved in the water flow. Obviously, the fracture connectivity is related to the fracture density: the greater the density, the greater the connectivity will be. Accordingly, the connectivity factor can be expressed as

    where0λ is the initial connectivity factor of microcrack groups,0d is the initial density of micro-crack groups, and d is the average density of micro-crack groups under the current load.

    Fig.2 The test curves of permeability of Lac du Bonnet granite (σc=5 MPa )

    Fig.3 The test curves of permeability of Lac du Bonnet granite (σc=10 MPa )

    Table 1 Parameters of numerical simulation of Lac du Bonnet granite in seepage test

    In order to determine the superimposed effect of micro-cracks in all directions, the density of microcracks in all directions is integrated over the unit sphere2=,=1 Sn n, and then the average density {} of micro-cracks in the macroscopic scale is obtained as

    where m denotes the number of integration points, and wτis the micro-crack integral weight in the nτdirection, which can be obtained by a look-up table using the integral method[15].

    According to assumption (4), combining Eq.(10) with Eq.(12) and following the superposition principle, the permeability coefficient and the permeability of multiple groups of micro-cracks can be obtained, respectively, as:

    It can be seen from Eq.(15) that the permeability of rock is not a constant during the process of loading, but rather a variable related to the micro-crack state within the rock mass, the density, the connectivity, and other factors.

    Fig.4 Numerical simulation of permeability behaviors of Lac du Bonnet granite (σc=5 MPa )

    2. Validation for the proposed model of permeability evolution

    To test whether the established evolution equation matches the real macro test characteristics of brittle rock, according to in situ and experimental data, the evolution law of the permeability tensor with the changes of stress and strain under load is numerically simulated.

    Fig.5 Numerical simulation of permeability behaviors of Lac du Bonnet granite (σc=10 MPa )

    Fig.6 The test curve of permeability-strain of Los ancient basalt

    2.1 Validation through in situ test

    In order to evaluate the simulation effectiveness of the proposed model under more complex conditions, the classic and universally accepted in situ test results are used to validate the permeability evolution model. For example, the Lac du Bounet granite, stored in the Nuclear Waste Underground Research Laboratory (URL) in Canada, was extensively studied in literature[13]. The permeability evolution curves of two groups of Lac du Bonnet granite samples, taken froma depth of 240 m underground in the Canadian URL, are shown in Figs.2 and 3. In figures, the (σ3?σ1) at the vertical axis indicates the deviatoric stress.

    Table 2 Parameters of numerical simulation of Los ancient basalt in seepage test

    Because the two groups of rock samples are taken from deep underground, most of the initial micro-cracks within the rock are in the closed state, as verified by images of electron microscopy. Thus the initial permeability of the rock is very low, within a range from 10 m2-18 m2to 10 m2-19 m2. It can be seen from Figs.2 and 3 that the permeability decreases by 1-2 magnitude in the elastic compressive stage and increases by 2-3 magnitude at the stage of slip and expansion of the micro-cracks. Because the initial permeability and the changes of permeability of the two groups of rock samples are quite different, two sets of parameters are used to simulate the permeability curves, as shown in Table 1.

    In Table 1, Esand vsare the elastic modulus and Poisson’s ratio of the rock matrix, respectively, η is a material constant used to describe the material’s expansibility for resisting the tensile stress, kcis used to describe the material’s surface energy when the micro-cracks are in the states of compression and shear, c is the rock cohesion, ρ is the critical friction of the micro-crack surface when the rock is in the states of compression and shear, f is used to describe the roughness of the fracture surface, Kniis the initial compressibility when the micro-cracks are in the state of elastic compression, b0and bcare the initial opening and the maximum closure opening of the micro-cracks, respectively, and λ0is the initial connectivity rate of micro-crack groups. The simulation results are shown in Figs.4 and 5.

    From Fig.4 and Fig.5, it is shown that the curves of the numerical simulation agree well with the test data with a clear indications of the four test stages of elastic compression, stable stage of permeability, dramatic increase of permeability, and slowing down of increase of permeability.

    2.2 Validation through laboratory test

    According to the permeability test under loading for the Los ancient basalt samples in Liangshan, Sichuan Province, the results are shown in Fig.6 and the test was carried out on MTS815.02 servo testing machine.

    Because the two rock samples are selected from the same basalt, their physical properties are similar and the mechanical laws they follow are relatively uniform. The parameters used in the numerical simulation are shown in Table 2 and the simulation results are shown in Fig.7 and Fig.8.

    Fig.7 Numerical simulations of permeability behaviors of Los ancient basalt (σc=5 MPa )

    Fig.8 Numerical simulations of permeability behaviors of Los ancient basalt (σc=10 MPa )

    From Fig.7 and Fig.8, it is shown that the proposed model can well describe the permeability evolution behaviors before the peak stage of Los ancient basalt, and also the relation between permeability evolution and confining pressure, but the post-peak stage of the permeability evolution behaviors can not be well described by this model, because the changes of permeability at the post-peak stage are very complex, which is not adequately reflected in the proposed model. In fact, the micro to macro-break phenomena of rock under loading are difficult to be described by the continuum mechanics alone, other angles should be considered.

    3. Conclusion

    The seepage field and the stress field coupling in the rock mass is an important issue in rock mechanics. The established mathematical models of seepagestress coupling mainly involve the relationships between the permeability and the stress or strain, without a due consideration of the anisotropy and nonhomogeneity of the permeability, which has restricted their applicability. In this study, based on the micromechanics theory and the evolution law of microcrack opening, density, connectivity, and the loading process, an anisotropic permeability evolution model is established and is simultaneously homogenized to the macroscopic scale by mathematical methods, and then the established evolution equation is validated by experimental data. Numerical simulation results show that the proposed model can well describe the permeability evolution law of different brittle rocks under different confining pressures, and especially, the before-peak macro-test phenomenon of permeability evolution behaviors on the macro-scale.

    [1] SCHULZE O., POPP T. and KERN H. Development of damage and permeability in deforming rock salt[J]. Engineering Geology, 2001, 61(2-3): 163-180.

    [2] WANG Xiao-dong, ZHOU Ying-fang and LUO Wanjing. A study on transient fluid flow of horizontal wells in dual-permeability media[J]. Journal of Hydrodynamics, 2010, 22(1): 44-50.

    [3] WANG J. A., PARK H. D. Fluid permeability of sedimentary rocks in a complete stress-strain process[J]. Engineering Geology, 2002, 63(3-4): 291-300.

    [4] ODA M. T., TAKEMURA A. and AOKI T. Damage growth and permeability change in triaxial compression tests of Inada granite[J]. Mechanics of Materials, 2002, 34(6): 313-331.

    [5] OLSSON R., BARTON N. An improved model for hydro mechanical coupling during shearing of rock joint[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(3): 317-329.

    [6] LEE H. S., CHO T. F. Hydraulics characteristics of rough fractures in linear flow under normal and shear load[J]. Rock Mechanics and Rock Engineering, 2002, 35(4): 299-318.

    [7] SEILER K. P., WOHNLICH S. New approaches characterizing groundwater flow[M]. Rotterdam, The Netherlands: A. A. Balkema, 2001, 845-846.

    [8] TANG C. A., THAM L. G. and LEE P. K. K. et al. Coupled analysis of flow, stress and damage (FSD) in rock failure[J]. International Journal of Rock Mechanics and Mining Science, 2002, 39(4): 477-489.

    [9] TAN Ye-fei, ZHOU Zhi-fang and HUANG Yong. Solute transport in natural fractures based on digital image technology[J]. Journal of Hydrodynamics, 2009, 21(2): 219-227.

    [10] LIU J., ELSWORTH D. and BRADY B. H. et al. Strain-dependent fluid flow defined through rock mass classification schemes[J]. Rock Mechanics and Rock Engineering, 2000, 33(3): 75-92.

    [11] TAKEMURA T., GOLSHANI A. and ODA M. et al. Preferred orientation of open microcrack in granite and their relation with anisotropic elasticity[J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(3): 443-454.

    [12] ZHOU J. J., SHAO J. F., XU W. Y. Coupled modeling of damage growth and permeability variation in brittle rocks[J]. Mechanics Research Communications, 2003, 33(4): 450-459.

    [13] SOULEY M., HOMAND F. and PEPA S. et al. Damage-induce permeability change in granite: A case study at the URL in Canada[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(2): 297-310.

    [14] WANG H. L., XU W. Y. and YANG S. Q. Experimental investigation on permeability evolution law during course of deformation and failure of rock specimen[J]. Rock and Soil Mechanics, 2006, 27(10): 101-106(in Chinese).

    [15] PENSEE V., KONDO D. and DORMIEUX L. Micromechanical analysis anisotropic damage in brittle materials[J]. Journal of Engineering Mechanics, 2002, 128(8): 889-897.

    10.1016/S1001-6058(11)60215-1

    * Project supported by the Natural National Science Foundation of China (Grant Nos. 51009052, 11172090), the Three Gorge Research Center for Geo-Hazards, the Ministry of Education (Grant No. TGRC201026), the Fundamental Research Funds for the Central Universities (Grant No. 2010B02514), and the Key Laboratory of Coastal Disasters and Defence of Ministry of Education, Hohai University (Grant No. 2010016).

    Biography: WANG Huan-ling (1976-), Female,

    Ph. D., Associate Professor

    2012,24(1):25-31

    免费看十八禁软件| 国产精品一区二区免费欧美| 757午夜福利合集在线观看| 12—13女人毛片做爰片一| 一a级毛片在线观看| 亚洲熟妇熟女久久| 久久久国产成人免费| 午夜福利成人在线免费观看| 熟妇人妻久久中文字幕3abv| 亚洲国产精品999在线| 在线天堂中文资源库| 久久久久九九精品影院| 伦理电影免费视频| 熟女电影av网| 久久久久久久久久黄片| 制服诱惑二区| 久久国产精品人妻蜜桃| 白带黄色成豆腐渣| 欧美av亚洲av综合av国产av| 免费女性裸体啪啪无遮挡网站| 热99re8久久精品国产| 国产精品电影一区二区三区| 男人舔女人的私密视频| 淫秽高清视频在线观看| 亚洲成a人片在线一区二区| 窝窝影院91人妻| 日本黄色视频三级网站网址| 99热只有精品国产| a在线观看视频网站| 欧美中文日本在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 久久中文字幕人妻熟女| 免费在线观看亚洲国产| 又黄又粗又硬又大视频| 人人妻人人澡欧美一区二区| 亚洲精品中文字幕一二三四区| 人人妻人人澡人人看| 午夜免费鲁丝| 日韩欧美 国产精品| 美女高潮喷水抽搐中文字幕| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 国产在线精品亚洲第一网站| 亚洲熟女毛片儿| 欧美黑人精品巨大| 满18在线观看网站| 欧美另类亚洲清纯唯美| 一边摸一边做爽爽视频免费| 国产一卡二卡三卡精品| 好男人在线观看高清免费视频 | 亚洲成av片中文字幕在线观看| 啪啪无遮挡十八禁网站| 日本三级黄在线观看| 婷婷亚洲欧美| 欧美午夜高清在线| 一边摸一边抽搐一进一小说| 成人18禁高潮啪啪吃奶动态图| 熟女电影av网| 变态另类丝袜制服| 国内精品久久久久久久电影| 不卡av一区二区三区| 亚洲美女黄片视频| av在线天堂中文字幕| 国内精品久久久久久久电影| 757午夜福利合集在线观看| 美女午夜性视频免费| 亚洲精品中文字幕一二三四区| 99久久综合精品五月天人人| 99国产精品一区二区三区| 亚洲自拍偷在线| 少妇裸体淫交视频免费看高清 | 伦理电影免费视频| 动漫黄色视频在线观看| 日韩三级视频一区二区三区| 亚洲欧美激情综合另类| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 中亚洲国语对白在线视频| svipshipincom国产片| 亚洲欧美精品综合一区二区三区| 亚洲男人天堂网一区| 国产亚洲av嫩草精品影院| 91老司机精品| 好男人电影高清在线观看| 啦啦啦 在线观看视频| 99精品久久久久人妻精品| 日日爽夜夜爽网站| 日韩成人在线观看一区二区三区| 亚洲av熟女| 色尼玛亚洲综合影院| 亚洲五月色婷婷综合| 在线免费观看的www视频| 91国产中文字幕| 日本一本二区三区精品| 少妇粗大呻吟视频| 免费观看人在逋| 中出人妻视频一区二区| 精品少妇一区二区三区视频日本电影| 亚洲熟妇熟女久久| 亚洲成国产人片在线观看| 天天一区二区日本电影三级| 国产又色又爽无遮挡免费看| 大香蕉久久成人网| 亚洲av中文字字幕乱码综合 | 成人三级做爰电影| 亚洲电影在线观看av| 国产精品一区二区三区四区久久 | 变态另类成人亚洲欧美熟女| 亚洲熟妇熟女久久| 狂野欧美激情性xxxx| 日韩精品中文字幕看吧| 国产精品久久久人人做人人爽| www.www免费av| 久久久久久久久免费视频了| 午夜视频精品福利| 亚洲精品中文字幕一二三四区| 久久伊人香网站| 国产亚洲精品久久久久久毛片| 亚洲av美国av| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 男女做爰动态图高潮gif福利片| 午夜免费激情av| 高潮久久久久久久久久久不卡| 一区二区三区高清视频在线| 狠狠狠狠99中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 国产黄片美女视频| 啦啦啦 在线观看视频| 日韩一卡2卡3卡4卡2021年| 无限看片的www在线观看| 久久狼人影院| 麻豆国产av国片精品| 中文字幕高清在线视频| 午夜a级毛片| 亚洲精品久久国产高清桃花| 精品高清国产在线一区| 国产真人三级小视频在线观看| 精品久久久久久久末码| 看片在线看免费视频| 亚洲美女黄片视频| 两性夫妻黄色片| 中文字幕高清在线视频| 国产又色又爽无遮挡免费看| 99精品久久久久人妻精品| 国产亚洲精品一区二区www| 在线观看一区二区三区| 一本一本综合久久| 搞女人的毛片| 哪里可以看免费的av片| 亚洲avbb在线观看| 精品少妇一区二区三区视频日本电影| 国内毛片毛片毛片毛片毛片| 亚洲天堂国产精品一区在线| 国产亚洲精品第一综合不卡| 又大又爽又粗| 女人爽到高潮嗷嗷叫在线视频| 村上凉子中文字幕在线| 一本久久中文字幕| 国产精品亚洲美女久久久| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 12—13女人毛片做爰片一| 嫩草影院精品99| 狂野欧美激情性xxxx| 国产精品亚洲一级av第二区| 免费高清在线观看日韩| 中文在线观看免费www的网站 | 长腿黑丝高跟| 成年女人毛片免费观看观看9| 激情在线观看视频在线高清| 在线天堂中文资源库| aaaaa片日本免费| 天天一区二区日本电影三级| 亚洲中文av在线| 可以免费在线观看a视频的电影网站| 一级黄色大片毛片| or卡值多少钱| 日韩欧美国产在线观看| 91成人精品电影| 精品久久久久久久末码| 欧美一级a爱片免费观看看 | 69av精品久久久久久| √禁漫天堂资源中文www| 色哟哟哟哟哟哟| av在线天堂中文字幕| 国内久久婷婷六月综合欲色啪| 12—13女人毛片做爰片一| а√天堂www在线а√下载| 91大片在线观看| 国产av在哪里看| 99热这里只有精品一区 | 黄片播放在线免费| 夜夜看夜夜爽夜夜摸| 999久久久国产精品视频| 久久精品影院6| 男人的好看免费观看在线视频 | 夜夜躁狠狠躁天天躁| 国产视频一区二区在线看| 欧美激情 高清一区二区三区| 99riav亚洲国产免费| 九色国产91popny在线| 欧美日韩黄片免| 1024香蕉在线观看| 狠狠狠狠99中文字幕| 国内精品久久久久久久电影| 国内久久婷婷六月综合欲色啪| tocl精华| 看免费av毛片| 一区二区三区国产精品乱码| 少妇粗大呻吟视频| 成人免费观看视频高清| 99热6这里只有精品| 欧美乱色亚洲激情| 午夜影院日韩av| 亚洲熟妇中文字幕五十中出| 欧美激情 高清一区二区三区| 欧美色视频一区免费| 99精品在免费线老司机午夜| 18禁观看日本| 成人免费观看视频高清| 亚洲五月婷婷丁香| 狂野欧美激情性xxxx| 国产99久久九九免费精品| 成人国产综合亚洲| 久久国产精品影院| 亚洲成人精品中文字幕电影| 淫秽高清视频在线观看| 俺也久久电影网| 欧美中文综合在线视频| 国产精品电影一区二区三区| 国产99白浆流出| 久久国产精品男人的天堂亚洲| 久久中文字幕人妻熟女| 国产国语露脸激情在线看| 一级a爱视频在线免费观看| 色综合站精品国产| 精品久久久久久久久久免费视频| 午夜福利一区二区在线看| 久久这里只有精品19| 久久香蕉精品热| 黄频高清免费视频| 欧美日本视频| 国产91精品成人一区二区三区| 99久久精品国产亚洲精品| 夜夜看夜夜爽夜夜摸| 欧美成人性av电影在线观看| av电影中文网址| 亚洲成av片中文字幕在线观看| 色综合欧美亚洲国产小说| 久久久久久九九精品二区国产 | 日本一区二区免费在线视频| 日韩中文字幕欧美一区二区| 国产精品久久久久久人妻精品电影| 看黄色毛片网站| 亚洲激情在线av| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 精品欧美国产一区二区三| 免费无遮挡裸体视频| 亚洲狠狠婷婷综合久久图片| 免费一级毛片在线播放高清视频| 中文字幕人妻丝袜一区二区| 国产aⅴ精品一区二区三区波| 亚洲国产精品999在线| 免费在线观看黄色视频的| 亚洲专区国产一区二区| 满18在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 国产午夜精品久久久久久| 国产成人啪精品午夜网站| 人人妻人人澡欧美一区二区| 亚洲人成伊人成综合网2020| ponron亚洲| 18美女黄网站色大片免费观看| av电影中文网址| 国产一区二区激情短视频| 国产av一区在线观看免费| 波多野结衣av一区二区av| 一本一本综合久久| 动漫黄色视频在线观看| av超薄肉色丝袜交足视频| 女人爽到高潮嗷嗷叫在线视频| 19禁男女啪啪无遮挡网站| 正在播放国产对白刺激| 日本 av在线| 欧美中文日本在线观看视频| 久久久国产欧美日韩av| 国产97色在线日韩免费| 成人国产一区最新在线观看| 日本一区二区免费在线视频| 精品电影一区二区在线| 真人一进一出gif抽搐免费| 婷婷精品国产亚洲av| 中文字幕另类日韩欧美亚洲嫩草| netflix在线观看网站| 变态另类丝袜制服| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 日本免费一区二区三区高清不卡| 一本一本综合久久| 身体一侧抽搐| 成年人黄色毛片网站| 婷婷精品国产亚洲av| 两个人视频免费观看高清| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 中文字幕最新亚洲高清| 国产精品久久久人人做人人爽| 成人三级做爰电影| 一a级毛片在线观看| 在线看三级毛片| 好男人电影高清在线观看| 后天国语完整版免费观看| 午夜两性在线视频| 禁无遮挡网站| 国产精品亚洲一级av第二区| 日韩高清综合在线| 国产精品免费视频内射| 丝袜人妻中文字幕| 国产精品99久久99久久久不卡| 两个人视频免费观看高清| 国产精品亚洲美女久久久| 天天添夜夜摸| 99久久99久久久精品蜜桃| 啦啦啦 在线观看视频| 国产精品一区二区免费欧美| 国产高清激情床上av| 国产欧美日韩精品亚洲av| 亚洲欧洲精品一区二区精品久久久| 美女免费视频网站| 99久久综合精品五月天人人| 天天躁狠狠躁夜夜躁狠狠躁| 99久久综合精品五月天人人| 欧美在线黄色| 人妻久久中文字幕网| 最好的美女福利视频网| 亚洲成国产人片在线观看| 人人妻人人澡欧美一区二区| 在线av久久热| 亚洲一区二区三区色噜噜| 少妇裸体淫交视频免费看高清 | 久久香蕉国产精品| 日本 av在线| 亚洲国产欧美一区二区综合| 波多野结衣巨乳人妻| avwww免费| 久久这里只有精品19| 日韩中文字幕欧美一区二区| 国产成人精品无人区| 久久 成人 亚洲| 久久精品国产亚洲av香蕉五月| 无人区码免费观看不卡| 国产黄片美女视频| 国产精品自产拍在线观看55亚洲| 国产又黄又爽又无遮挡在线| 91成人精品电影| 最近在线观看免费完整版| 亚洲午夜精品一区,二区,三区| 中文字幕人妻丝袜一区二区| 欧美一级毛片孕妇| 国内少妇人妻偷人精品xxx网站 | 日韩欧美 国产精品| 最近最新中文字幕大全免费视频| 一级毛片女人18水好多| 亚洲美女黄片视频| 麻豆国产av国片精品| 国产视频一区二区在线看| 免费观看人在逋| 级片在线观看| 成人国产一区最新在线观看| 99国产综合亚洲精品| 午夜久久久久精精品| 女人高潮潮喷娇喘18禁视频| 好男人在线观看高清免费视频 | 免费高清视频大片| 国产97色在线日韩免费| 亚洲三区欧美一区| 国产亚洲av嫩草精品影院| 国产午夜福利久久久久久| 亚洲色图 男人天堂 中文字幕| 黄色 视频免费看| 男人操女人黄网站| 999久久久国产精品视频| 国产成人av教育| 自线自在国产av| 美女 人体艺术 gogo| 午夜免费成人在线视频| www国产在线视频色| 99re在线观看精品视频| 久久伊人香网站| 1024香蕉在线观看| 欧美最黄视频在线播放免费| 久久久久久大精品| 国产伦在线观看视频一区| 91av网站免费观看| 国产乱人伦免费视频| 最近最新免费中文字幕在线| 亚洲熟女毛片儿| 精品欧美国产一区二区三| 欧美一区二区精品小视频在线| 在线观看免费午夜福利视频| 久久人人精品亚洲av| 99久久综合精品五月天人人| 他把我摸到了高潮在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品九九99| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 亚洲av五月六月丁香网| 日韩高清综合在线| 一区二区三区高清视频在线| 亚洲免费av在线视频| 国产成人系列免费观看| 久久久久精品国产欧美久久久| 亚洲成人国产一区在线观看| 久久青草综合色| 琪琪午夜伦伦电影理论片6080| 欧美乱码精品一区二区三区| 亚洲av第一区精品v没综合| www日本在线高清视频| 欧美黑人巨大hd| 久久精品91蜜桃| 99在线视频只有这里精品首页| 97碰自拍视频| 特大巨黑吊av在线直播 | 免费一级毛片在线播放高清视频| 香蕉丝袜av| 白带黄色成豆腐渣| 欧美成狂野欧美在线观看| 天天躁夜夜躁狠狠躁躁| 在线av久久热| 国产精品久久久久久精品电影 | 黄片小视频在线播放| www.熟女人妻精品国产| 亚洲国产欧洲综合997久久, | 日日干狠狠操夜夜爽| 亚洲午夜精品一区,二区,三区| 又大又爽又粗| 亚洲 国产 在线| 99久久精品国产亚洲精品| 国产真人三级小视频在线观看| 久久久久精品国产欧美久久久| 一二三四社区在线视频社区8| 成人一区二区视频在线观看| 国产精品av久久久久免费| 欧美国产精品va在线观看不卡| 男女视频在线观看网站免费 | 黄片小视频在线播放| 日日爽夜夜爽网站| 制服丝袜大香蕉在线| 欧美精品亚洲一区二区| 欧美色视频一区免费| 99国产极品粉嫩在线观看| 成年女人毛片免费观看观看9| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女| 亚洲色图av天堂| 久久午夜亚洲精品久久| 久久久久亚洲av毛片大全| 人妻久久中文字幕网| 免费在线观看完整版高清| 看片在线看免费视频| 51午夜福利影视在线观看| 女同久久另类99精品国产91| 中文字幕人妻熟女乱码| 国产精品久久久av美女十八| 亚洲男人的天堂狠狠| 亚洲,欧美精品.| 精品高清国产在线一区| 人成视频在线观看免费观看| 黄色 视频免费看| 亚洲免费av在线视频| 久久精品国产亚洲av香蕉五月| 啦啦啦观看免费观看视频高清| 国产精品精品国产色婷婷| 777久久人妻少妇嫩草av网站| 亚洲色图 男人天堂 中文字幕| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 中文字幕人妻丝袜一区二区| 美国免费a级毛片| 麻豆av在线久日| 12—13女人毛片做爰片一| 女人被狂操c到高潮| 久久99热这里只有精品18| 亚洲精品av麻豆狂野| 精品国产乱码久久久久久男人| 日本成人三级电影网站| 国产麻豆成人av免费视频| 午夜福利欧美成人| 亚洲在线自拍视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成人免费电影在线观看| 在线看三级毛片| 久久99热这里只有精品18| 中文字幕av电影在线播放| 18禁美女被吸乳视频| 欧美成人午夜精品| 午夜老司机福利片| 国产高清激情床上av| 成人欧美大片| 免费人成视频x8x8入口观看| 国内少妇人妻偷人精品xxx网站 | 久久久国产成人精品二区| 91麻豆精品激情在线观看国产| 国产成人精品久久二区二区免费| 一级片免费观看大全| 日韩免费av在线播放| 久久精品国产99精品国产亚洲性色| 国产精品,欧美在线| 欧美色视频一区免费| 国产亚洲精品第一综合不卡| 久久国产精品人妻蜜桃| 国内毛片毛片毛片毛片毛片| 老司机在亚洲福利影院| 日韩av在线大香蕉| 国语自产精品视频在线第100页| 免费看十八禁软件| 欧美大码av| 身体一侧抽搐| 黄色视频,在线免费观看| 美女国产高潮福利片在线看| a级毛片在线看网站| 黄网站色视频无遮挡免费观看| 亚洲aⅴ乱码一区二区在线播放 | 啦啦啦观看免费观看视频高清| 欧美性猛交黑人性爽| 无限看片的www在线观看| 免费人成视频x8x8入口观看| 老熟妇乱子伦视频在线观看| 国产成人影院久久av| 国产亚洲精品久久久久久毛片| 中亚洲国语对白在线视频| 精品一区二区三区视频在线观看免费| 99久久无色码亚洲精品果冻| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片 | 久久中文字幕人妻熟女| 又黄又粗又硬又大视频| 亚洲七黄色美女视频| 免费无遮挡裸体视频| 日韩高清综合在线| 身体一侧抽搐| 无人区码免费观看不卡| 巨乳人妻的诱惑在线观看| 伦理电影免费视频| 亚洲av成人av| 国产伦在线观看视频一区| 18美女黄网站色大片免费观看| 精品国产超薄肉色丝袜足j| 高清毛片免费观看视频网站| 精品国产超薄肉色丝袜足j| 操出白浆在线播放| 国产野战对白在线观看| 精品日产1卡2卡| 18禁国产床啪视频网站| 我的亚洲天堂| 他把我摸到了高潮在线观看| 在线播放国产精品三级| 女生性感内裤真人,穿戴方法视频| 美国免费a级毛片| 久久精品成人免费网站| 欧美 亚洲 国产 日韩一| 欧美zozozo另类| 婷婷精品国产亚洲av| 好男人在线观看高清免费视频 | 国产日本99.免费观看| 一区二区日韩欧美中文字幕| 男女那种视频在线观看| 久久久水蜜桃国产精品网| 国产亚洲精品av在线| 成人午夜高清在线视频 | 欧美日韩福利视频一区二区| 人人妻,人人澡人人爽秒播| 国产成人av教育| av在线播放免费不卡| 免费观看精品视频网站| 国产伦人伦偷精品视频| 欧美日韩瑟瑟在线播放| 青草久久国产| 久久久国产欧美日韩av| 91大片在线观看| 国产精品爽爽va在线观看网站 | www.999成人在线观看| 成人永久免费在线观看视频| 欧美成人一区二区免费高清观看 | 欧美最黄视频在线播放免费| 国产伦人伦偷精品视频| 日韩精品中文字幕看吧| 男人操女人黄网站| 日本五十路高清| 久久伊人香网站| videosex国产| 一级a爱视频在线免费观看| 1024手机看黄色片| 波多野结衣高清无吗| 欧洲精品卡2卡3卡4卡5卡区| 午夜影院日韩av| 久久久久精品国产欧美久久久| 成人亚洲精品一区在线观看| 婷婷精品国产亚洲av| 美女 人体艺术 gogo| 国产蜜桃级精品一区二区三区| 看免费av毛片| 18禁美女被吸乳视频| 亚洲国产高清在线一区二区三 | 变态另类成人亚洲欧美熟女| 中文字幕精品亚洲无线码一区 |