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

    堤防土體滲透參數(shù)的概率分布研究

    2009-12-27 23:16:12李少龍朱國勝定培中張文三
    長江科學(xué)院院報 2009年4期

    李少龍,朱國勝,定培中,張文三

    堤防土體滲透參數(shù)的概率分布研究

    李少龍,朱國勝,定培中,張文三

    (長江科學(xué)院水利部巖土力學(xué)與工程重點實驗室,武漢 430010)

    確定滲透參數(shù)的概率分布是堤防滲透穩(wěn)定可靠性分析的基礎(chǔ)。闡述了土體滲透系數(shù)和臨界水力比降的概率分布研究概況。通過分布參數(shù)的極大似然估計和A-D法分布擬合檢驗,分析了堤防粉質(zhì)粘土和粉質(zhì)壤土的滲透系數(shù)概率分布;基于太沙基公式,將土粒相對密度和孔隙比作為隨機變量,推導(dǎo)了臨界水力比降的概率密度函數(shù)。分析表明,滲透系數(shù)的概率分布符合對數(shù)正態(tài)分布,臨界比降的概率分布為近似正態(tài)分布。

    滲透系數(shù);臨界水力比降;概率分布;堤防

    我國堤防多是經(jīng)過一定歷史時期不斷加高培厚而形成的,且其地質(zhì)情況復(fù)雜,堤線較長,受不確定性因素影響眾多,往往存在滲透穩(wěn)定問題。堤防土體多具有不均勻性,與堤防滲流密切相關(guān)的滲透系數(shù)與臨界水力比降等參數(shù)具有不同程度的變異性。然而,用傳統(tǒng)的方法進行堤防滲流安全評估時,通常假定土體分區(qū)均勻,將滲透系數(shù)與臨界水力比降在分區(qū)內(nèi)作為常數(shù)。近些年發(fā)展起來的可靠性分析方法較好地考慮了堤防土體參數(shù)的隨機性和變異性。

    在滲透穩(wěn)定可靠性分析中,確定滲透參數(shù)的概率分布是必不可少的工作??煽慷确治龇椒?,如目前應(yīng)用較多的一次二階矩方法,需要根據(jù)基本隨機變量的概率分布進行當(dāng)量正態(tài)化來計算可靠指標(biāo)和失效概率。隨機參數(shù)的概率分布對可靠度計算的最終計算結(jié)果具有重要影響。參數(shù)概率分布規(guī)律是利用可靠度理論進行堤防滲透穩(wěn)定分析的基礎(chǔ),然而目前對堤防土體滲透參數(shù)特別是抗?jié)B強度概率分布的研究較少。本文分析土體滲透系數(shù)與抗?jié)B強度的概率分布特征及分析方法,為堤防滲透穩(wěn)定的可靠性分析提供基礎(chǔ)。

    1 土體滲透參數(shù)的隨機特征

    江河堤防的土質(zhì)組成極其復(fù)雜,堤基土多為第四系沖積層,常含有粘土、壤土、細(xì)砂、砂礫石等。各種土的特性指標(biāo)相差較大,無粘性土和粘性土的滲透變形特性有顯著不同。同時,在宏觀上土體性質(zhì)被視為均勻的土層,其物理力學(xué)性質(zhì)也會隨著位置的不同而有所差異。土的顆?;蚩紫对诳臻g上作隨機排列,顆粒之間的作用力也是隨機分布的,土體存在非均勻性和隨機性,同時土體所處的物理狀態(tài)及環(huán)境條件可能不同,因而土性參數(shù)具有較大的隨機變異性。鑒于堤防土體的復(fù)雜性,研究其滲透參數(shù)概型分布的推斷方法及概率分布規(guī)律,可加深認(rèn)識巖土介質(zhì)的隨機特性,對堤防滲透穩(wěn)定可靠性分析具有重要意義。

    隨機水文地質(zhì)領(lǐng)域的研究表明,滲透系數(shù)在空間上的分布變化是極其復(fù)雜的,含水層的滲透系數(shù)在空間上的變化可達(dá)幾個數(shù)量級[1]。R.A.Freeze[2]統(tǒng)計大量滲透系數(shù)空間變異性的結(jié)果,認(rèn)為滲透系數(shù)可用對數(shù)正態(tài)分布來描述。這一結(jié)論被后來的研究者在隨機地下水運動的研究中廣泛應(yīng)用。E.A.Sudicky[3]在加拿大Borden含水層進行野外試驗,經(jīng)過分析得出滲透系數(shù)服從對數(shù)正態(tài)分布的結(jié)論。眾多試驗研究和隨機模擬都認(rèn)可滲透系數(shù)服從對數(shù)正態(tài)分布,即對數(shù)滲透系數(shù)服從正態(tài)分布,這在數(shù)學(xué)處理上也帶來較大方便,其統(tǒng)計特征可用均值和協(xié)方差函數(shù)完全描述。此外,根據(jù)中心極限定理,受大量相互獨立的隨機因素影響,每一個別因素在總影響中所起作用不太大的隨機變量,仍服從或近似服從正態(tài)分布。因此,正態(tài)分布或?qū)?shù)正態(tài)分布在隨機研究中得到廣泛應(yīng)用。

    抗?jié)B強度反映了土體抵抗?jié)B透破壞的能力,可用臨界水力比降來表征。臨界水力比降是滲透穩(wěn)定性分析中最基本的參數(shù),受多種不確定性因素影響,土體的臨界比降也是一個不確定性量。文獻(xiàn)[4]指出土層的抗?jié)B強度各處不同,同時孔隙的大小不一,分布也不均勻,這就使得各處的滲流坡降也大小不一,有的地方滲流坡降相對增大,或者抗?jié)B強度比較弱,就有可能發(fā)生管涌。在防洪搶險的實踐中發(fā)現(xiàn),距堤腳一定范圍內(nèi)的不同位置都可能出現(xiàn)規(guī)模不一的險情。這就從一個方面說明堤基土層分布不均勻,實際發(fā)生的滲透比降和臨界比降具有變異性。將臨界比降作為隨機變量,進行堤防工程可靠性分析時需要知道它的概率分布。然而關(guān)于土體滲透變形參數(shù)的隨機特征研究不多,在建立堤防滲流穩(wěn)定可靠性分析模型時往往對其概率分布進行某種假設(shè)。王卓甫等[5]建立了防洪堤滲透變形的風(fēng)險計算模型,假定臨界比降服從三角分布。周小文等[6]在堤防自動化監(jiān)測預(yù)警系統(tǒng)中,建立了滲流安全概率評價模型,假設(shè)臨界比降服從正態(tài)分布。張士辰和李雷[7]通過粗砂滲透變形20組室內(nèi)實驗,對比分析了流土型粗砂滲透系數(shù)與抗?jié)B強度概型分布,認(rèn)為對數(shù)正態(tài)分布的擬合效果最好。

    2 基于試驗數(shù)據(jù)的滲透系數(shù)概率分布推斷

    在土工概率計算時,應(yīng)該先要對土性參數(shù)試驗數(shù)據(jù)進行必要的統(tǒng)計分析,充分利用有限的資料得到盡可能多的規(guī)律性認(rèn)識。隨機土性參數(shù)的均值、方差、分布類型等統(tǒng)計資料往往需要較長時間的積累。如果有足夠的試驗數(shù)據(jù),那么可以根據(jù)觀測樣本,選用適宜的理論概率分布,用統(tǒng)計檢驗的方法選擇最佳的分布類型。

    2.1 概率分布類型

    對于土體滲透參數(shù),用某個適當(dāng)?shù)姆植紒砻枋鏊?,根?jù)對參數(shù)概率分布問題的分析和經(jīng)驗,可以大體確定出幾種可能的分布概型。巖土工程中比較常用的分布類型包括:正態(tài)分布、對數(shù)正態(tài)分布和極值Ⅰ型分布等,它們的概率密度函數(shù)分別如下:

    式中:μ和σ分別為正態(tài)分布的均值和標(biāo)準(zhǔn)差;a和b分別為極值Ⅰ型分布的尺度參數(shù)和位置參數(shù)。

    2.2 參數(shù)的極大似然估計

    推斷總體的分布類型時,通常它的某些分布參數(shù)也是未知的,這時要對未知的分布參數(shù)作出推斷,即參數(shù)估計,一般可用矩法和極大似然法確定。極大似然估計法的計算較復(fù)雜,但是能充分利用分布函數(shù)對參數(shù)所提供的信息,因而具有優(yōu)良性質(zhì)。根據(jù)已知或假定的分布類型,以樣本的概率作為似然函數(shù),再以樣本似然函數(shù)最大為準(zhǔn)則,確定相應(yīng)的分布參數(shù)。對于連續(xù)型隨機變量,其密度函數(shù)為f(x;θ),θ∈Θ為待估參數(shù),樣本的觀測值是x1,x2,…,xn,其似然函數(shù)為L(θ)= ∏f(xi;θ),似然函數(shù)的最大值點即為參數(shù)的極大似然估計值。根據(jù)以上原理和方法可以得到:

    正態(tài)分布參數(shù)的極大似然估計為

    對數(shù)正態(tài)分布參數(shù)的極大似然估計為

    極值Ⅰ型分布參數(shù)的極大似然估計為

    正態(tài)分布和對數(shù)正態(tài)分布的參數(shù)極大似然估計值可用樣本數(shù)據(jù)根據(jù)式(4)和(5)直接計算;由式(6)可見,極值分布參數(shù)的極大似然估計值滿足非線性方程,不能直接計算其解析解,本文采用牛頓-拉夫森迭代法求解。

    2.3 分布擬合檢驗的A-D法

    在實際問題中,不知道總體具體服從什么類型的分布,這就需要根據(jù)樣本來檢驗關(guān)于分布的假設(shè),看是否在概率統(tǒng)計意義上合理,即分布擬合檢驗。常用的檢驗方法有:χ2檢驗法、K-S檢驗法和A-D檢驗法等。χ2檢驗法是用劃分區(qū)間的方法來考慮樣本頻率與理論概率之間的偏差,依賴于區(qū)間的劃分,而且要求樣本容量足夠大。K-S法和A-D法是在每一點上考慮經(jīng)驗分布與理論分布之間的偏差,因而沒有依賴于區(qū)間劃分的缺點。可靠度計算中隨機變量的尾部分布較為重要,A-D法能較好反映所檢驗的經(jīng)驗分布和母體分布在尾部的擬合情況,檢驗精度較高,對大子樣檢驗和小子樣檢驗都適用[8]。本文采用A-D檢驗法,其步驟如下:

    (1)將樣本觀測值x1,x2,…,xn按由小到大次序進行重排得到x(1)≤x(2)≤…x(n);

    (2)計算統(tǒng)計量

    (3)在顯著性水平α下,如果A2小于臨界值A(chǔ)2α,則接受原假設(shè);否則拒絕原假設(shè)。

    2.4 滲透系數(shù)概率分布

    對于堤防工程,由于堤線較長,現(xiàn)場和室內(nèi)試驗不會在一個較小范圍內(nèi)取樣,基本上沿堤線布置取樣,而少量的試驗數(shù)據(jù)又不滿足統(tǒng)計分析的要求,因而可考慮在一個較大的范圍內(nèi)按土性分類收集測試數(shù)據(jù),同一類土其基本特性應(yīng)當(dāng)是相同或相近似的,試驗條件也應(yīng)基本相同。這樣就易于滿足統(tǒng)計分析的要求,也充分利用了已有資料。此外,不同土類滲透系數(shù)的取值范圍不同,其統(tǒng)計特征也有一定差異,可以檢驗它們是否來自同一母體。

    應(yīng)用上述原理方法,對長江中游某段堤防的粉質(zhì)粘土和粉質(zhì)壤土的滲透系數(shù)概率分布進行擬合檢驗,這兩類土的試驗數(shù)據(jù)較多,分別有73組和54組數(shù)據(jù)。表1給出了它們的分布參數(shù)極大似然估計和A-D法檢驗結(jié)果。分析可見,滲透系數(shù)在對數(shù)正態(tài)分布的假設(shè)條件下計算統(tǒng)計量小于臨界值,而在正態(tài)分布和極值Ⅰ型分布的假定下計算統(tǒng)計量大于臨界值,說明滲透系數(shù)的概率分布符合對數(shù)正態(tài)分布。采用t檢驗法和F檢驗法分別對兩類土的滲透系數(shù)均值和方差是否來自同一母體進行檢驗,結(jié)果表明在1%顯著水平下接受兩類土的方差來自同一母體,但是拒絕接受兩類土的均值來自同一母體,說明它們的統(tǒng)計特征存在一定差異。

    3 隨機變量函數(shù)的概率密度推導(dǎo)

    土的臨界水力比降一般通過試驗測定,所獲得的參數(shù)相對準(zhǔn)確可靠,但是需要花費大量的人力物力和財力,并非所有的工程項目都具備開展大量試驗的條件,有時根據(jù)類似工程或經(jīng)驗來選取參數(shù),存在一定的主觀性。實際工程中經(jīng)常遇到樣本容量不足甚至沒有試驗數(shù)據(jù)的情況。半經(jīng)驗半理論的計算公式為土體臨界比降的確定提供了一種方法,同時也為其概率分布的推斷提供了一條途徑。

    對于滲流自下而上的情況,由單位體積土在水中的浮重和作用于該土的滲透力相平衡的原理,可得到流土的臨界水力比降,即太沙基公式

    式中:J為臨界水力比降;G為土粒相對質(zhì)量密度;e為土的孔隙比。堤防工程地質(zhì)勘察規(guī)程(SL188-2005)[9]中推薦該公式計算無粘性土的流土臨界水力比降。除了上述太沙基公式之外,國內(nèi)外研究機構(gòu)和學(xué)者還提出了其它計算臨界水力比降的數(shù)學(xué)公式[10,11],可供采用。計算公式反映了土體臨界比降的主要影響因素,能夠滿足工程需要,這些因素又各有隨機性,我們可以通過計算公式根據(jù)影響因子的分布特征來分析臨界水力比降的概率分布。下面基于太沙基公式推導(dǎo)臨界水力比降的概率密度函數(shù)。

    太沙基公式中有2個自變量,即土粒相對密度和孔隙比,通常認(rèn)為土粒相對密度的變化范圍不大而將其作為常量。有關(guān)研究表明,土體孔隙比服從正態(tài)分布。不失一般性,本文考慮2種情況,即僅將孔隙比作為隨機變量以及將兩者都作為隨機變量。

    根據(jù)式(8),由分布函數(shù)的定義有

    于是得到J的概率密度表達(dá)式為

    由孔隙比的正態(tài)分布概率密度可得臨界水力比降的概率密度函數(shù)為

    將土粒相對密度和孔隙比都作為隨機變量,可以看到,土粒相對密度是土顆粒與同體積純水的質(zhì)量比值,而孔隙比是孔隙體積與固體顆粒體積之比,從概念上來理解兩者相關(guān)性不大,此處認(rèn)為兩隨機變量相互獨立,聯(lián)合概率密度fG,e(g,e)=fG(g)fe(e)。

    表1 滲透系數(shù)概型分布擬合A-D法檢驗表Table 1 A-D test of probability distribution of permeability

    構(gòu)造隨機變量G和e的兩個函數(shù)U和V,即U=(G-1)/(1+e),V=G,其反函數(shù)為 G=V,e=(V-1)/U-1,易知該變換的雅克比行列式為J=(V-1)/U2,由二維隨機變量函數(shù)的分布公式可得臨界比降的概率密度為

    上式關(guān)于臨界比降的概率密度表達(dá)式含有正態(tài)分布的積分,由于形式復(fù)雜較難給出顯式的解析表達(dá),在實際運用中也會有一些困難。為此可考慮將臨界比降公式按照泰勒級數(shù)展開,即

    式中:a0=(μG-1)/(1+μe);a1=1/(1+μe);a2=-(μG-1)/(1+μe)2;μG和μe分別為土粒相對密度和孔隙比的均值。相互獨立的正態(tài)隨機變量之和仍服從正態(tài)分布,于是可得臨界水力比降的分布為

    式中,σ2G和σ2e分別為土粒相對密度和孔隙比的方差。

    4 結(jié)語

    滲透參數(shù)的概率分布規(guī)律是利用可靠度理論進行堤防滲透穩(wěn)定分析的基礎(chǔ)。本文對土體滲透參數(shù)的概率分布研究進行了小結(jié),針對巖土工程中常用的正態(tài)分布、對數(shù)正態(tài)分布和極值Ⅰ型分布,通過分布參數(shù)的極大似然估計和A-D法分布擬合檢驗,分析了堤防土體滲透系數(shù)的概率分布。極大似然法能充分利用分布函數(shù)對參數(shù)所提供的信息;A-D法能較好地反映所檢驗的經(jīng)驗分布和母體分布在尾部的擬合情況,適合可靠度分析中的分布檢驗。分析表明粉質(zhì)粘土和粉質(zhì)壤土的滲透系數(shù)服從對數(shù)正態(tài)分布。基于太沙基公式,將孔隙比作為隨機變量直接得到臨界比降的概率密度,將孔隙比和土粒相對密度都作為隨機變量,通過泰勒級數(shù)推導(dǎo)出臨界水力比降的近似正態(tài)分布概率密度函數(shù)。研究結(jié)果可為堤防工程滲透穩(wěn)定的可靠性分析提供基礎(chǔ)。

    [1] 楊金忠,蔡樹英,黃冠華,等.多孔介質(zhì)中水分及溶質(zhì)運移的隨機理論[M].北京:科學(xué)出版社,2000.

    [2] FREEZE R A.A Stochastic Conceptual Analysis of One Dimensional Groundwater Flow in Nonuniform Homogeneous Media[J].Water Resour.Res.,1975,11(5):725-741.

    [3] SUDICKY E A.A Natural Gradient Experiment on Solute Transport in a Sand Aquifer:Spatial Variability of Hydraulic Conductivity and Its Role in the Dispersion Process[J].Water Resources Research,1986,22(13):2069-2082.

    [4] 曹敦侶,曹 罡,鄒火元,等.水工建筑物滲流管涌的Monte-Carlo模擬[J].人民長江,1997,28(6):11-13.

    [5] 王卓甫,章志強,楊高升.防洪堤結(jié)構(gòu)風(fēng)險計算模型探討[J].水利學(xué)報,1998,29(7):64-67.

    [6] 周小文,包偉力,吳昌瑜.基于堤防監(jiān)測數(shù)據(jù)的安全評價模型研究[J].大壩觀測與土工測試,2001,27(5):17-20.

    [7] 張士辰,李 雷.粗砂滲透系數(shù)與抗?jié)B強度概型分布[J].水利水運工程學(xué)報,2004,(3):58-61.

    [8] 吳世偉.結(jié)構(gòu)安全度與可靠度分析論文集[M].南京:河海大學(xué)出版社,1988.

    [9] SL188-2005,堤防工程地質(zhì)勘察規(guī)程[S].

    [10]劉 杰.土的滲透穩(wěn)定與滲流控制[M].北京:水利電力出版社,1992.

    [11]毛昶熙.滲流計算分析與控制[M].北京:水利電力出版社,1990.

    Probability Distribution of Levee Foundation Seepage Parameters

    LI Shao-long,ZHU Guo-sheng,DING Pei-zhong,Zhang Wen-san
    (Key Laboratory of Geotechnical Mechanics and Engineering of the Ministry of Water Resources,Yangtze River Scientific Research Institute,Wuhan 430010,China)

    The determination of probability distribution of seepage parametersisa basisof reliability analysis of leveeseepage stability.The research development of soil permeability and critical hydraulic gradient is presented.The probability distribution of levee soil permeability is analyzed by maximum likelihood estimation of distribution parameters and A-D test method for distribution fitting.On the basis of Terzaghi formula of flowing soil,the probability density function of critical hydraulic gradient is derived with soil relative density and void ratio as random variables.Results show that the distributions of the permeability of silty clay and silty loam are both subject to lognormal distribution,the soil critical hydraulic gradient obeys approximately normal distribution.

    permeability;critical hydraulic gradient;probability distribution;levee

    TB114

    A

    1001-5485(2009)04-0036-04

    2008-06-24;

    2008-10-17

    長江科學(xué)院中央級公益科研院所基金資助(YWF0733/YT03);長江科學(xué)院博士科研啟動基金資助(YJJ0790/TY30);國家科技支撐計劃項目基金資助(2006BAC05B05-01)

    李少龍(1979-),男,湖北嘉魚人,工程師,主要從事水工滲流及地下水環(huán)境研究,(電話)027-82927243(電子信箱)lisl3000@163.com。

    (編輯:劉運飛)

    欧美又色又爽又黄视频| 蜜桃久久精品国产亚洲av| 老鸭窝网址在线观看| 欧美性猛交╳xxx乱大交人| 久久久久久久午夜电影| 国产午夜精品久久久久久| 免费电影在线观看免费观看| 日韩成人在线观看一区二区三区| 亚洲自拍偷在线| 人人妻,人人澡人人爽秒播| 国内少妇人妻偷人精品xxx网站 | 日韩欧美精品v在线| 免费看a级黄色片| 嫩草影视91久久| 神马国产精品三级电影在线观看| 免费在线观看影片大全网站| 国产一区二区三区视频了| 我的老师免费观看完整版| 亚洲精品在线美女| 九九在线视频观看精品| 亚洲无线观看免费| 小蜜桃在线观看免费完整版高清| 一区二区三区国产精品乱码| 国内精品久久久久精免费| 亚洲欧洲精品一区二区精品久久久| 一卡2卡三卡四卡精品乱码亚洲| 免费在线观看亚洲国产| 国内精品久久久久久久电影| 在线免费观看不下载黄p国产 | 亚洲七黄色美女视频| cao死你这个sao货| 欧美在线一区亚洲| 国产精品久久久久久人妻精品电影| 色精品久久人妻99蜜桃| 精品熟女少妇八av免费久了| 无人区码免费观看不卡| 麻豆成人午夜福利视频| 久久天躁狠狠躁夜夜2o2o| 男人的好看免费观看在线视频| 久久久色成人| 欧美日韩福利视频一区二区| 99精品久久久久人妻精品| av片东京热男人的天堂| 老司机午夜福利在线观看视频| 成人av在线播放网站| 日韩高清综合在线| 久久九九热精品免费| 午夜福利视频1000在线观看| 99精品久久久久人妻精品| 色视频www国产| 美女高潮的动态| 国产久久久一区二区三区| 啦啦啦观看免费观看视频高清| 国产伦精品一区二区三区四那| 99国产精品99久久久久| 色综合婷婷激情| h日本视频在线播放| 国产一区二区三区视频了| 黄色 视频免费看| 2021天堂中文幕一二区在线观| 久久精品aⅴ一区二区三区四区| 国产高清视频在线观看网站| 九九热线精品视视频播放| 女生性感内裤真人,穿戴方法视频| 香蕉丝袜av| 亚洲精华国产精华精| tocl精华| 久久久久性生活片| 精品日产1卡2卡| 天堂影院成人在线观看| 日本一本二区三区精品| 精品久久久久久久人妻蜜臀av| 欧美日韩国产亚洲二区| 曰老女人黄片| 亚洲成人免费电影在线观看| 久久久久亚洲av毛片大全| 亚洲成人中文字幕在线播放| 国产欧美日韩一区二区三| 深夜精品福利| 成人特级av手机在线观看| 岛国在线免费视频观看| 18禁裸乳无遮挡免费网站照片| 午夜精品久久久久久毛片777| 国产成人一区二区三区免费视频网站| 久久久国产成人精品二区| 精品99又大又爽又粗少妇毛片 | 日日摸夜夜添夜夜添小说| 亚洲自偷自拍图片 自拍| 中文亚洲av片在线观看爽| 日韩精品青青久久久久久| av国产免费在线观看| 老汉色av国产亚洲站长工具| 精品国产乱码久久久久久男人| 啪啪无遮挡十八禁网站| 99精品久久久久人妻精品| 久久久久国产一级毛片高清牌| 男女那种视频在线观看| 国产成人精品久久二区二区免费| 天天一区二区日本电影三级| 久久精品91蜜桃| 波多野结衣高清无吗| 亚洲美女视频黄频| 琪琪午夜伦伦电影理论片6080| 两个人视频免费观看高清| 18禁裸乳无遮挡免费网站照片| bbb黄色大片| 亚洲国产精品久久男人天堂| 亚洲av中文字字幕乱码综合| 婷婷亚洲欧美| 1000部很黄的大片| 99re在线观看精品视频| 岛国视频午夜一区免费看| 国产精品av久久久久免费| 色精品久久人妻99蜜桃| 日韩有码中文字幕| 人妻久久中文字幕网| 两个人看的免费小视频| 韩国av一区二区三区四区| 亚洲欧美精品综合一区二区三区| 老汉色∧v一级毛片| 一a级毛片在线观看| АⅤ资源中文在线天堂| 性色av乱码一区二区三区2| 久久久久九九精品影院| 亚洲成av人片免费观看| 国产v大片淫在线免费观看| 亚洲在线观看片| 国产午夜精品论理片| 欧美不卡视频在线免费观看| 亚洲精品在线观看二区| 亚洲美女视频黄频| 琪琪午夜伦伦电影理论片6080| 欧美成人性av电影在线观看| 午夜精品在线福利| 又爽又黄无遮挡网站| 亚洲在线观看片| 天天一区二区日本电影三级| 国产高清视频在线播放一区| 日韩 欧美 亚洲 中文字幕| 九色成人免费人妻av| 亚洲av第一区精品v没综合| 在线观看免费午夜福利视频| 1024手机看黄色片| 很黄的视频免费| 亚洲精品456在线播放app | 中文在线观看免费www的网站| aaaaa片日本免费| 亚洲国产中文字幕在线视频| 成人av在线播放网站| 最新美女视频免费是黄的| 亚洲av片天天在线观看| 国产精品98久久久久久宅男小说| 日本免费一区二区三区高清不卡| 久久精品夜夜夜夜夜久久蜜豆| 哪里可以看免费的av片| 99国产精品一区二区蜜桃av| 国产亚洲av嫩草精品影院| 午夜激情福利司机影院| 熟女少妇亚洲综合色aaa.| 亚洲中文日韩欧美视频| 日本黄色视频三级网站网址| 日本与韩国留学比较| 婷婷六月久久综合丁香| 成年女人永久免费观看视频| 精品免费久久久久久久清纯| 真人一进一出gif抽搐免费| 成年女人毛片免费观看观看9| 欧美中文综合在线视频| 一级毛片高清免费大全| 免费av不卡在线播放| 一级a爱片免费观看的视频| 成年女人永久免费观看视频| 一边摸一边抽搐一进一小说| 亚洲av电影在线进入| 亚洲在线观看片| 身体一侧抽搐| 亚洲精华国产精华液的使用体验| 菩萨蛮人人尽说江南好唐韦庄 | 日韩高清综合在线| 精品无人区乱码1区二区| 国产高清三级在线| 亚洲欧美日韩东京热| 日产精品乱码卡一卡2卡三| 变态另类丝袜制服| 亚洲婷婷狠狠爱综合网| 成年免费大片在线观看| 18禁动态无遮挡网站| 国产午夜精品一二区理论片| 免费观看在线日韩| 久久亚洲精品不卡| 最近视频中文字幕2019在线8| 一个人观看的视频www高清免费观看| 最近2019中文字幕mv第一页| 国产亚洲av片在线观看秒播厂 | 免费看日本二区| 国产综合懂色| 精品少妇黑人巨大在线播放 | 一本一本综合久久| 晚上一个人看的免费电影| 一级毛片aaaaaa免费看小| 免费无遮挡裸体视频| 边亲边吃奶的免费视频| 国产成人91sexporn| 一夜夜www| 男人和女人高潮做爰伦理| 国产精品久久电影中文字幕| 欧美日韩国产亚洲二区| 纵有疾风起免费观看全集完整版 | 波野结衣二区三区在线| 亚洲欧美日韩无卡精品| 精品久久久久久久久av| av国产久精品久网站免费入址| 亚洲欧美精品综合久久99| 亚洲国产日韩欧美精品在线观看| 国产高潮美女av| 色吧在线观看| 国产免费视频播放在线视频 | 国产真实乱freesex| 国产又黄又爽又无遮挡在线| 亚洲无线观看免费| 男女啪啪激烈高潮av片| 日韩三级伦理在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品人妻熟女av久视频| 久久人妻av系列| 老司机影院毛片| 少妇人妻一区二区三区视频| 国产亚洲最大av| 国产精品熟女久久久久浪| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 男女啪啪激烈高潮av片| 国产高清视频在线观看网站| 一本一本综合久久| 日韩欧美国产在线观看| 亚洲欧美日韩东京热| 久久热精品热| 国产亚洲91精品色在线| 99久久精品国产国产毛片| 尤物成人国产欧美一区二区三区| 国产成人a区在线观看| .国产精品久久| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 99热这里只有是精品在线观看| 可以在线观看毛片的网站| 麻豆久久精品国产亚洲av| 九色成人免费人妻av| 天天躁日日操中文字幕| 极品教师在线视频| 亚洲,欧美,日韩| 村上凉子中文字幕在线| 国产精品日韩av在线免费观看| 插阴视频在线观看视频| 蜜桃久久精品国产亚洲av| 两性午夜刺激爽爽歪歪视频在线观看| 久久欧美精品欧美久久欧美| 色综合站精品国产| 好男人在线观看高清免费视频| 在线播放国产精品三级| 亚洲国产精品合色在线| 国产一区二区三区av在线| 少妇裸体淫交视频免费看高清| 久久人人爽人人片av| 赤兔流量卡办理| 亚洲精品自拍成人| 中文精品一卡2卡3卡4更新| 能在线免费看毛片的网站| 国产在视频线在精品| 亚洲中文字幕日韩| 七月丁香在线播放| 国产在视频线精品| 真实男女啪啪啪动态图| 国内精品一区二区在线观看| 欧美日本亚洲视频在线播放| 91久久精品国产一区二区三区| 国产黄色视频一区二区在线观看 | 美女国产视频在线观看| 亚洲高清免费不卡视频| 99久久成人亚洲精品观看| 一个人看的www免费观看视频| 久久这里只有精品中国| 国产免费视频播放在线视频 | 免费看av在线观看网站| 中文字幕人妻熟人妻熟丝袜美| 内射极品少妇av片p| 久久精品久久久久久噜噜老黄 | 国产精品人妻久久久久久| 日韩人妻高清精品专区| 中文字幕熟女人妻在线| 男人和女人高潮做爰伦理| 成人美女网站在线观看视频| 国产一级毛片七仙女欲春2| 国产欧美另类精品又又久久亚洲欧美| 性色avwww在线观看| 久久精品国产亚洲av天美| 免费看光身美女| 色噜噜av男人的天堂激情| 毛片一级片免费看久久久久| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| 在线观看一区二区三区| 亚洲国产精品成人综合色| 欧美丝袜亚洲另类| 国产伦精品一区二区三区四那| 日韩人妻高清精品专区| 亚洲av电影在线观看一区二区三区 | videossex国产| 直男gayav资源| 国内精品美女久久久久久| 长腿黑丝高跟| 久久精品国产亚洲网站| 秋霞伦理黄片| 寂寞人妻少妇视频99o| 97在线视频观看| 国产69精品久久久久777片| 国产在线男女| 欧美日韩在线观看h| av在线播放精品| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 亚洲四区av| 99久久精品一区二区三区| 人人妻人人澡欧美一区二区| 深爱激情五月婷婷| 免费黄网站久久成人精品| 女人久久www免费人成看片 | 又爽又黄无遮挡网站| 国产伦在线观看视频一区| www日本黄色视频网| 卡戴珊不雅视频在线播放| 中文在线观看免费www的网站| 99视频精品全部免费 在线| 色哟哟·www| 99热6这里只有精品| 日韩国内少妇激情av| 久久亚洲精品不卡| 一本一本综合久久| 97超碰精品成人国产| 国产一区二区三区av在线| 精品久久久久久久末码| 国产成人精品一,二区| 精品一区二区三区人妻视频| 男女国产视频网站| 极品教师在线视频| 熟女电影av网| 床上黄色一级片| 国产精品一二三区在线看| 欧美激情国产日韩精品一区| 亚洲欧美成人综合另类久久久 | 全区人妻精品视频| 亚洲最大成人av| 亚洲久久久久久中文字幕| 国产极品精品免费视频能看的| 久久久精品94久久精品| 国产69精品久久久久777片| 在线观看66精品国产| 色哟哟·www| 中文字幕制服av| 国产午夜福利久久久久久| 国产精品三级大全| 美女高潮的动态| 中文亚洲av片在线观看爽| 亚洲精品国产成人久久av| 超碰97精品在线观看| 欧美区成人在线视频| 久久久久久久午夜电影| 欧美极品一区二区三区四区| 亚洲欧美精品自产自拍| 久久久精品欧美日韩精品| 内地一区二区视频在线| 三级国产精品欧美在线观看| 久久综合国产亚洲精品| av线在线观看网站| 又粗又爽又猛毛片免费看| 色播亚洲综合网| 国产欧美日韩精品一区二区| 热99在线观看视频| 久久99热6这里只有精品| 大香蕉久久网| 久久久成人免费电影| 欧美xxxx性猛交bbbb| 99久久九九国产精品国产免费| 久久人人爽人人爽人人片va| a级毛色黄片| 人人妻人人澡人人爽人人夜夜 | 91午夜精品亚洲一区二区三区| 日日干狠狠操夜夜爽| 色综合色国产| 久久久久久伊人网av| 少妇熟女欧美另类| 亚洲成人av在线免费| 免费不卡的大黄色大毛片视频在线观看 | 99久久精品热视频| 国产亚洲av嫩草精品影院| 搞女人的毛片| 亚洲av免费在线观看| 成人二区视频| 国产精品一区二区性色av| 精品免费久久久久久久清纯| 亚洲在线自拍视频| 国产精品无大码| 春色校园在线视频观看| 成人午夜精彩视频在线观看| 久久精品久久久久久噜噜老黄 | 欧美激情久久久久久爽电影| 国产男人的电影天堂91| 干丝袜人妻中文字幕| 寂寞人妻少妇视频99o| 欧美性感艳星| 九九爱精品视频在线观看| 天天躁夜夜躁狠狠久久av| 国产免费视频播放在线视频 | 日韩欧美在线乱码| 男女下面进入的视频免费午夜| 欧美色视频一区免费| 亚洲欧美日韩东京热| 国产男人的电影天堂91| 亚洲久久久久久中文字幕| 午夜福利高清视频| 午夜亚洲福利在线播放| 精品久久久久久久久久久久久| 国产精品国产三级专区第一集| 亚洲国产最新在线播放| kizo精华| 在线免费观看的www视频| 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 日韩强制内射视频| 乱码一卡2卡4卡精品| 国产精华一区二区三区| 亚洲欧美中文字幕日韩二区| 久久亚洲国产成人精品v| 97超视频在线观看视频| 女人被狂操c到高潮| 国产成人精品久久久久久| 亚洲欧美日韩卡通动漫| 精品久久久久久久人妻蜜臀av| 日韩精品有码人妻一区| 国产真实乱freesex| 欧美激情久久久久久爽电影| 国产精品三级大全| av黄色大香蕉| 久久99蜜桃精品久久| 99在线视频只有这里精品首页| 国产黄片美女视频| 我要搜黄色片| 日韩视频在线欧美| 精品久久久久久电影网 | 欧美三级亚洲精品| 一夜夜www| 午夜免费男女啪啪视频观看| 国产欧美另类精品又又久久亚洲欧美| 欧美成人精品欧美一级黄| 22中文网久久字幕| 亚洲va在线va天堂va国产| 日韩国内少妇激情av| 国产亚洲午夜精品一区二区久久 | 久久婷婷人人爽人人干人人爱| 麻豆av噜噜一区二区三区| 亚洲国产高清在线一区二区三| 丝袜美腿在线中文| 亚洲欧美成人精品一区二区| 亚洲精品乱久久久久久| 嘟嘟电影网在线观看| 亚洲精品亚洲一区二区| 18禁在线无遮挡免费观看视频| 国产伦一二天堂av在线观看| 亚洲av成人av| 亚洲图色成人| 亚洲国产精品sss在线观看| 简卡轻食公司| 日本黄色视频三级网站网址| 亚洲五月天丁香| 内射极品少妇av片p| 亚洲av成人av| 久久人妻av系列| 婷婷六月久久综合丁香| 免费黄色在线免费观看| 欧美区成人在线视频| 久久99热这里只频精品6学生 | 中文亚洲av片在线观看爽| 免费观看性生交大片5| 成人午夜精彩视频在线观看| 日韩大片免费观看网站 | 成人午夜精彩视频在线观看| 久久久久久伊人网av| 如何舔出高潮| 国产三级在线视频| 成人av在线播放网站| av.在线天堂| 欧美成人a在线观看| 少妇的逼水好多| 高清av免费在线| 国产日韩欧美在线精品| 久久久久久伊人网av| 中文字幕制服av| 精品人妻视频免费看| 99热全是精品| 午夜日本视频在线| 最近2019中文字幕mv第一页| 人妻夜夜爽99麻豆av| 国产亚洲精品久久久com| 免费av不卡在线播放| 亚洲av一区综合| 亚洲欧美成人精品一区二区| 久久久久久久午夜电影| 午夜视频国产福利| 男女国产视频网站| 成人综合一区亚洲| 日日撸夜夜添| 亚洲av电影不卡..在线观看| 国产乱人视频| 简卡轻食公司| 永久免费av网站大全| 亚洲精品国产av成人精品| 少妇熟女aⅴ在线视频| 久久久久久久午夜电影| 啦啦啦韩国在线观看视频| 国产高清视频在线观看网站| 国产又黄又爽又无遮挡在线| 美女国产视频在线观看| 国产精品久久久久久久电影| 国产精品女同一区二区软件| 亚洲国产日韩欧美精品在线观看| 九色成人免费人妻av| 久久精品熟女亚洲av麻豆精品 | 国产老妇女一区| 亚洲aⅴ乱码一区二区在线播放| 国产午夜精品一二区理论片| 99久久中文字幕三级久久日本| 99热这里只有精品一区| av免费观看日本| 久久久国产成人精品二区| 自拍偷自拍亚洲精品老妇| 免费黄色在线免费观看| av.在线天堂| 禁无遮挡网站| 亚洲中文字幕日韩| 久久午夜福利片| 一级黄色大片毛片| 国产精品久久久久久av不卡| 久久鲁丝午夜福利片| 69人妻影院| 岛国毛片在线播放| 精品免费久久久久久久清纯| 精品少妇黑人巨大在线播放 | 中文字幕av在线有码专区| 国产亚洲午夜精品一区二区久久 | 丰满乱子伦码专区| 中文字幕久久专区| 国产国拍精品亚洲av在线观看| 亚洲在线自拍视频| 亚州av有码| 久久久久久久亚洲中文字幕| 国产精品,欧美在线| 亚洲国产精品国产精品| 国产精品一及| 国产大屁股一区二区在线视频| 国产精品国产三级国产av玫瑰| 真实男女啪啪啪动态图| 欧美精品国产亚洲| 日本一本二区三区精品| 黑人高潮一二区| 午夜精品一区二区三区免费看| 国产男人的电影天堂91| 亚洲丝袜综合中文字幕| 又粗又硬又长又爽又黄的视频| 国产精品电影一区二区三区| 精品久久久久久久久av| 精品国内亚洲2022精品成人| 欧美精品国产亚洲| 夜夜爽夜夜爽视频| 精品久久久久久久久av| 亚洲三级黄色毛片| 成人特级av手机在线观看| 美女被艹到高潮喷水动态| 欧美一区二区国产精品久久精品| 99久久精品一区二区三区| 中文字幕熟女人妻在线| 国产精品99久久久久久久久| av视频在线观看入口| 国产高清不卡午夜福利| 日本一二三区视频观看| 身体一侧抽搐| 午夜激情欧美在线| 亚洲精品日韩av片在线观看| 亚洲不卡免费看| 亚洲精品色激情综合| 国产精品野战在线观看| 欧美性猛交╳xxx乱大交人| 日韩,欧美,国产一区二区三区 | 久久久久免费精品人妻一区二区| 午夜亚洲福利在线播放| 少妇高潮的动态图| av国产久精品久网站免费入址| 最近中文字幕2019免费版| 99久久精品国产国产毛片| 欧美精品一区二区大全| a级毛色黄片| 婷婷色综合大香蕉| 色尼玛亚洲综合影院| 亚洲av免费在线观看| 99久久精品一区二区三区| 乱人视频在线观看| 国产伦一二天堂av在线观看| 免费搜索国产男女视频| 国产毛片a区久久久久| 一级毛片我不卡| 1000部很黄的大片| 国产视频内射| 午夜a级毛片| 日本一本二区三区精品| 有码 亚洲区|