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

    基于產(chǎn)狀Fisher分布的巖石節(jié)理最小樣本容量確定方法

    2024-02-02 02:02:14楊梓楓夏華宗杜玉芳王曉明
    河北地質(zhì)大學學報 2024年1期
    關鍵詞:離散性樣本容量產(chǎn)狀

    楊梓楓,夏華宗,杜玉芳,王曉明,3

    1.河北地質(zhì)大學 城市地質(zhì)與工程學院,河北 石家莊 050031;2.河北地礦建設工程集團有限責任公司,河北 石家莊 050000;3.河北省地下人工環(huán)境智慧開發(fā)與管控技術創(chuàng)新中心,河北 石家莊 050031

    0 引言

    節(jié)理是巖體中普遍發(fā)育的一種結構面,這些結構面將完整巖體切割成大小不等、形狀各異的巖石塊體,而且節(jié)理通常強度低、易變形,同時能夠為地下水提供存儲空間和運移通道,對巖體的完整性、力學特性、滲流特性及穩(wěn)定性起著控制作用。節(jié)理產(chǎn)狀是影響巖體結構及其工程性質(zhì)的關鍵參數(shù)之一,其平均值及離散性對邊坡、地下洞室等工程巖體的穩(wěn)定性及可靠性具有顯著影響[1,2]。因此,準確計算節(jié)理的平均產(chǎn)狀及其離散性對掌握巖體的工程特性進而開展工程巖體的評價、設計及施工具有重要意義。

    節(jié)理平均產(chǎn)狀的準確性受樣本容量(即采樣數(shù)量)和離散性的影響,產(chǎn)狀離散性越大,則節(jié)理的采樣數(shù)量應越大,反之亦然。為保證統(tǒng)計結果的準確性,節(jié)理樣本應不少于某一數(shù)量,該數(shù)量即為最小樣本容量。在實際工程中,受采樣條件的限制,如節(jié)理露頭較少或測量困難而造成采樣數(shù)量較少,平均產(chǎn)狀的統(tǒng)計準確性勢必會受到影響。因此,確定節(jié)理的最小樣本容量對于節(jié)理統(tǒng)計分析具有現(xiàn)實意義。

    巖體中的節(jié)理往往具有一定的方向性,通常會沿2~5個優(yōu)勢方向發(fā)育,產(chǎn)狀相近的節(jié)理應劃分為同一組。節(jié)理的產(chǎn)狀一般用傾向(或走向)和傾角表示,同一組節(jié)理的產(chǎn)狀通常服從一定的分布形式。常見的節(jié)理分布類型包括Fisher分布、Bingham分布和雙正態(tài)分布等[3-6]。其中,F(xiàn)isher分布由于僅包含1個參數(shù)且為可積分函數(shù)便于實現(xiàn)Monte Carlo模擬,得到了廣泛應用[7-12]。此外,一些研究表明,實測節(jié)理產(chǎn)狀能夠較好地擬合Fisher分布[13,14]。由于節(jié)理發(fā)育的隨機性和離散性,在工程實踐中通常將節(jié)理進行分組統(tǒng)計來確定每組節(jié)理的平均產(chǎn)狀[15-17]。

    基于節(jié)理產(chǎn)狀服從Fisher分布,本文探討產(chǎn)狀離散性對節(jié)理采樣數(shù)量的影響。通過Monte Carlo隨機模擬技術產(chǎn)生具有不同離散性的產(chǎn)狀隨機數(shù),在保證準確性的條件下確定節(jié)理的最小樣本容量,并建立最小樣本容量與產(chǎn)狀離散性的經(jīng)驗關系,最后通過算例驗證其適用性。

    1 Fisher分布及產(chǎn)狀離散性

    該分布于1953年由Fisher提出,是一種將節(jié)理產(chǎn)狀視為球面上關于平均產(chǎn)狀對稱的分布形式[3]。Fisher分布的概率密度函數(shù)為:

    (1a)

    (1b)

    式中,θ′和φ′分別是當把坐標系z軸旋轉至節(jié)理的平均矢量方向(θm,φm)后在新坐標系中的方向角;參數(shù)κ反映了節(jié)理產(chǎn)狀分布的離散性或集中程度,κ越大表明節(jié)理的產(chǎn)狀圍繞平均產(chǎn)狀越集中,即離散性小。如果κ=0,則Fisher分布轉化為球面上的均勻分布。均質(zhì)巖體中的構造節(jié)理產(chǎn)狀通常較為集中,具有較大的κ值;而非均質(zhì)巖體中的節(jié)理產(chǎn)狀往往具有較大的離散性,其κ值通常較小。κ可按下式計算[3]:

    (2)

    式中,N為節(jié)理的數(shù)量;rN為節(jié)理的和向量,其方向即為節(jié)理的平均產(chǎn)狀,其計算過程將在2.2節(jié)中詳述。Watson的研究認為當κ>3時,公式(2)的計算結果具有足夠精度[18]。眾多研究表明[10-12],同一組節(jié)理的κ一般大于3,因此公式(2)具有較好的適用性。

    2 研究方法

    2.1 Fisher分布產(chǎn)狀的隨機生成

    節(jié)理的方向是以其單位法向量表示的,一般取指向上方的法向量。在直角坐標系中(圖1),假設OP為某節(jié)理的單位法向量,點N為P點在xOy平面上的垂直投影,θ為OP與z軸正方向的夾角,φ為x軸正方向沿逆時針方向與ON的夾角。

    Fisher分布是當把z軸旋轉至節(jié)理的平均矢量方向(θm,φm)之后的形式,此時它的兩個變量θ′和φ′是獨立的,其概率密度函數(shù)見公式(1)。其中φ′服從(0,2π)之間的均勻分布,可直接利用公式(3)產(chǎn)生,u為(0,1)范圍內(nèi)的均勻分布隨機數(shù)。

    φ′=2πu

    (3)

    θ′的密度函數(shù)是一個可積分函數(shù),其累計分布函數(shù)為:

    (4)

    這樣可利用逆變換法產(chǎn)生Fisher分布的隨機數(shù),即:

    (5)

    得到:

    (6)

    如果稱旋轉前(z軸鉛直向上)的坐標系為舊坐標系,它的變量為(φ,θ);稱旋轉后(z軸與節(jié)理的平均矢量一致)的坐標系為新坐標系,它的變量為(φ′,θ′),則兩個坐標系的變化關系為[7,9]:

    式中:(l,m,n)和(l′,m′,n′)分別為對應(φ,θ)和(φ′,θ′)的方向余弦。l,m,n分別為:

    (8)

    地質(zhì)學中,節(jié)理的方向是用傾向和傾角(α,β)表示的,而非(φ,θ)。如果采用圖3中的坐標系,(α,β)和(φ,θ)的關系可以用下式計算:

    (9)

    綜合上述分析,產(chǎn)生一組N個Fisher分布的隨機節(jié)理產(chǎn)狀的過程如圖2所示。

    圖2 Fisher分布隨機節(jié)理產(chǎn)狀的產(chǎn)生流程Fig.2 Generation flowchart of random joint orientation obeying Fisher distribution

    2.2 節(jié)理平均產(chǎn)狀計算

    節(jié)理平均產(chǎn)狀是指節(jié)理組中每一條節(jié)理面的單位法向量之和所在的方向。假設一組節(jié)理包含N條節(jié)理,每一條節(jié)理的傾向和傾角分別為φi和θi,通過公式(8)可以得到每條節(jié)理的單位法向量。由于節(jié)理產(chǎn)狀的離散性,節(jié)理產(chǎn)狀的分布范圍有時較大,當采用上半球赤平極射投影時可能出現(xiàn)“超越上半球投影的產(chǎn)狀”(OEBEUHP)[9,13]。圖3為200個隨機產(chǎn)生的服從Fisher分布的節(jié)理產(chǎn)狀,其平均產(chǎn)狀為45°∠60°,κ=10,其中11個節(jié)理極點為OEBEUHP。

    圖3 平均產(chǎn)狀為45°∠60°的節(jié)理產(chǎn)狀極點圖(κ=10)Fig.3 Pole diagram of joint orientations with mean orientation of 45°∠60°(κ=10)

    由于存在OEBEUHP,同一組節(jié)理向上的法向量之間的夾角可能大于90°,所以在計算節(jié)理平均產(chǎn)狀時,不能將節(jié)理的法向量ni簡單相加,應考慮采用ni還是-ni。鄭俊等提出了一種識別和調(diào)整OEBEUHP的策略,通過計算ni與合向量rN的夾角來識別OEBEUHP,并用-ni重新計算rN[9,13]。實際上,初始rN的計算已經(jīng)存在由OEBEUHP引起的問題,因此這種策略需要適當調(diào)整。本文在此基礎上,提出了一種計算rN的修正方法,具體過程見圖4。根據(jù)得到的rN,可利用公式(8)和(9)將其轉換為節(jié)理的平均產(chǎn)狀。節(jié)理產(chǎn)狀的隨機生成和平均產(chǎn)狀的計算均通過Python程序實現(xiàn)。

    圖4 節(jié)理和向量rN的計算過程Fig.4 Calculation process of joint and vector rN

    3 最小節(jié)理采樣數(shù)量的確定

    由于節(jié)理產(chǎn)狀的離散性,節(jié)理平均產(chǎn)狀的計算結果受樣本大小的影響。產(chǎn)狀離散性越大,即κ值越小,節(jié)理的樣本容量應越大才能保證平均產(chǎn)狀計算的準確性。本文以平均產(chǎn)狀為45°∠60°的節(jié)理為例,通過隨機生成不同數(shù)量、不同κ值的Fisher分布產(chǎn)狀來研究樣本容量對平均產(chǎn)狀計算準確性的影響,κ值分別取5,10,15,20,25,30,35,40,45,50,55,60,70,80,90和100。對于每一個κ值,共進行1 000次抽樣,樣本容量為ns。圖5繪制了一次隨機實現(xiàn)時不同κ值的節(jié)理極點分布圖,樣本容量ns為250。

    若每次抽樣的平均產(chǎn)狀與45°∠60°的夾角小于5°,則認為滿足準確性要求。將滿足以上要求的抽樣次數(shù)記為np,那么準確率p可以定義為:

    (10)

    當κ≤40時,每次抽樣的樣本容量ns為5的倍數(shù);當κ>40時,ns為自然數(shù)。

    圖6反映了不同κ值時樣本大小對平均產(chǎn)狀準確率的影響。樣本容量對平均產(chǎn)狀的統(tǒng)計準確性影響顯著,樣本容量越大,平均產(chǎn)狀的統(tǒng)計準確率越高;產(chǎn)狀離散性(即κ值)對統(tǒng)計準確率也有明顯影響,離散性越大即κ值越小,達到一定準確率所需的樣本數(shù)越多。一般情況下,準確率p為0.95可滿足工程要求,因此將p=0.95對應的樣本容量定義為最小樣本容量。當樣本大小超過最小樣本容量時,可以保證產(chǎn)狀的統(tǒng)計準確率不低于0.95。圖7為具有不同κ值的節(jié)理產(chǎn)狀所應滿足的最小樣本容量ns,min,可以看出,ns,min隨著κ值的增大迅速減小,當κ值達到40以后,ns,min的變化很小且趨于水平。為方便在工程實踐中確定節(jié)理的最小樣本容量,采用反比例函數(shù)(公式(11))對ns,min和κ值的關系進行擬合,得到了圖7中的擬合曲線。公式(11)的確定系數(shù)R2為0.994,表明該式能準確反映ns,min與κ值的反比例關系,可用于確定節(jié)理的最小樣本容量。

    圖6 不同κ值時樣本容量對統(tǒng)計準確率的影響Fig.6 Effect of sample size on statistical accuracy with different κ values

    圖7 不同κ值對應的最小樣本容量Fig.7 Minimum sample size corresponding to different κ values

    (11)

    實際上,在對節(jié)理產(chǎn)狀進行采樣之前無法預知節(jié)理組的κ值。但隨著采樣過程的進行,會得到越來越多的產(chǎn)狀數(shù)據(jù),可以利用這些數(shù)據(jù)對κ值進行估算,從而判斷最終的采樣數(shù)量ns是否滿足公式(11)的要求。如果ns不小于ns,min,則采樣數(shù)量滿足要求,平均產(chǎn)狀的準確率能夠得到保障。本文的研究基于產(chǎn)狀服從Fisher分布,實際巖體中大多數(shù)構造節(jié)理均成組發(fā)育且服從Fisher分布,因此公式(11)對構造節(jié)理具有較好的適用性。對于孤立的隨機節(jié)理,其發(fā)育沒有規(guī)律性,不適合采用公式(11)進行計算。

    4 驗證算例

    為驗證公式(11)的適用性,對三峽工程地下廠房的上游邊墻和烏東德水電站平洞PD5-2節(jié)理產(chǎn)狀進行了抽樣統(tǒng)計分析。三峽水電站是世界上最大的水利樞紐之一,位于中國湖北省宜昌市夷陵區(qū)。地下電站位于長江右岸白巖尖山體內(nèi),主廠房洞室為圓拱直墻型,全長 311.3 m,軸線走向 223.5°,最大跨度 32.6 m,最大高度 87.3 m。廠房圍巖為堅硬的花崗巖,巖體完整~較完整,整體穩(wěn)定性較好。地下電站廠房洞室的開挖揭露了708條節(jié)理,分組后的節(jié)理極點分布情況見圖8,每組節(jié)理產(chǎn)狀的統(tǒng)計參數(shù)見表1。烏東德水電站平洞PD5-2位于金沙江烏東德壩址左岸850 m高程,為順江支洞。平洞洞形為城門洞形,斷面尺寸約2 m×2 m。本平洞共揭露節(jié)理354條,分組后的節(jié)理極點分布情況見圖9,每組節(jié)理產(chǎn)狀的統(tǒng)計參數(shù)見表2。

    表1 三峽地下廠房上游邊墻節(jié)理產(chǎn)狀統(tǒng)計結果Table 1 Statistical results of the orientation of joints in the upstream side wall of the Three Gorges underground powerhouse

    表2 烏東德水電站平洞PD5-2節(jié)理產(chǎn)狀統(tǒng)計結果Table 2 Statistical results of the orientation of joints in the adit PD5-2 of the Wudongde hydropower station

    圖8 三峽地下廠房上游邊墻節(jié)理極點分布圖Fig.8 Distribution diagram of joint poles in the upstream side wall of the Three Gorges underground powerhouse

    圖9 烏東德水電站平洞PD5-2節(jié)理極點分布圖Fig.9 Distribution diagram of joint poles in the adit PD5-2 of the Wudongde hydropower station

    根據(jù)公式(11),計算出三峽地下廠房上游邊墻每組節(jié)理的最小樣本容量應達到97、97和105。接下來,對每組節(jié)理產(chǎn)狀進行模擬抽樣,抽樣次數(shù)均為1 000次,樣本大小為最小樣本容量。最后利用公式(10)計算抽樣結果的準確率。經(jīng)程序計算,每組節(jié)理平均產(chǎn)狀的準確率分別為0.988、0.992和0.994,均滿足準確率要求。該研究區(qū)每組節(jié)理的采樣數(shù)量均大于相應的最小樣本容量,滿足采樣要求。因此,三峽工程地下廠房的上游邊墻的節(jié)理產(chǎn)狀具有統(tǒng)計學意義。

    根據(jù)公式(11),計算得到烏東德水電站平洞PD5-2每組節(jié)理的最小樣本容量應達到87、100和104。假設每組節(jié)理采樣數(shù)分別為80、90和95,未達到相應組別的最小樣本容量。利用公式(10),計算出每組抽樣結果的準確率分別為0.955、0.948和0.946。結果表明,除第一組采樣能達到準確率以外,其余兩組未能滿足準確率要求。抽樣過程具有一定的隨機性和偶然性,在抽樣數(shù)量不足時,仍有較小可能達到準確率,但是為保證抽樣結果具有較高的可信度,節(jié)理抽樣應達到最小樣本容量的要求。本節(jié)通過正反兩個算例驗證了公式(11)的可靠性,表明節(jié)理最小樣本容量ns,min與κ值的經(jīng)驗關系適用于確定節(jié)理產(chǎn)狀的最小樣本容量。

    5 結論

    本文基于節(jié)理產(chǎn)狀服從Fisher分布,實現(xiàn)了節(jié)理產(chǎn)狀的隨機生成。為解決“超越上半球投影的產(chǎn)狀”(OEBEUHP)引起平均產(chǎn)狀計算不準確的問題,提出了一種計算平均產(chǎn)狀的修正方法。在此基礎上,隨機生成不同數(shù)量、不同離散性(κ值)的Fisher分布產(chǎn)狀來研究節(jié)理樣本容量對平均產(chǎn)狀準確性的影響,并建立了最小樣本容量與κ值的經(jīng)驗關系。研究結果表明:節(jié)理樣本容量對平均產(chǎn)狀的統(tǒng)計準確性影響顯著,樣本容量越大,平均產(chǎn)狀的統(tǒng)計準確率越高;產(chǎn)狀離散性對統(tǒng)計準確率也有明顯影響,離散性越大即κ值越小,達到一定準確率所需的樣本數(shù)越多。為保證平均產(chǎn)狀的統(tǒng)計準確率不低于0.95,節(jié)理的最小樣本容量隨著Fisher分布常數(shù)κ值的增大迅速減小;二者之間具有明顯的反比例關系,反比例系數(shù)為1 007.23,經(jīng)驗公式的確定系數(shù)為0.994。對三峽地下廠房上游邊墻和烏東德水電站平洞PD5-2的節(jié)理統(tǒng)計分析從正反兩方面驗證了本文方法的適用性,節(jié)理采樣具有一定的隨機性,在滿足最小樣本容量的條件下,產(chǎn)狀統(tǒng)計分析結果具有較高的可靠性。本文提出的方法適用于確定節(jié)理產(chǎn)狀的最小采樣數(shù)量,判斷節(jié)理的采樣數(shù)量是否具有統(tǒng)計意義,能有效指導節(jié)理的采樣工作。

    猜你喜歡
    離散性樣本容量產(chǎn)狀
    淺談砂巖儲層的巖石學特征
    激電聯(lián)合剖面在判斷矽卡巖型礦床礦體產(chǎn)狀中的應用
    “三點解析法”估算地質(zhì)體產(chǎn)狀及應用
    采用無核密度儀檢測壓實度的樣本容量確定方法
    有關Cliford矩陣群的幾個離散性判別準則
    赤平投影法在邊坡穩(wěn)定性分析中的應用
    熱水表的應用研究及與其他熱計量方式的比較分析
    均值與方差例析
    一類具冪指積系數(shù)微分算子譜的離散性
    廣義高斯分布參數(shù)估值與樣本容量關系
    測繪通報(2013年2期)2013-12-11 07:27:44
    国产v大片淫在线免费观看| 亚洲av电影不卡..在线观看| 午夜日韩欧美国产| 亚洲 国产 在线| 国产精品一区二区三区四区久久| 18禁裸乳无遮挡免费网站照片| 成人鲁丝片一二三区免费| 淫妇啪啪啪对白视频| 日韩精品中文字幕看吧| 搡女人真爽免费视频火全软件 | 91av网一区二区| 最近视频中文字幕2019在线8| 亚洲va在线va天堂va国产| 亚洲精品456在线播放app | 亚洲av熟女| 熟女人妻精品中文字幕| 国产精品人妻久久久久久| 在线观看66精品国产| 少妇丰满av| 亚洲七黄色美女视频| 久久99热6这里只有精品| 可以在线观看毛片的网站| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 国产淫片久久久久久久久| 国产在线精品亚洲第一网站| 99久久九九国产精品国产免费| 久久欧美精品欧美久久欧美| 91久久精品国产一区二区三区| 1024手机看黄色片| 午夜视频国产福利| 精品人妻1区二区| 国产一区二区激情短视频| 毛片女人毛片| 联通29元200g的流量卡| 又粗又爽又猛毛片免费看| 亚洲精品粉嫩美女一区| 一进一出抽搐动态| 日韩在线高清观看一区二区三区 | 亚洲va在线va天堂va国产| 亚洲精品乱码久久久v下载方式| 国产精品乱码一区二三区的特点| 成人国产麻豆网| 免费大片18禁| 男人狂女人下面高潮的视频| 成年版毛片免费区| 一进一出好大好爽视频| 日本一二三区视频观看| 国产av麻豆久久久久久久| 精品一区二区三区人妻视频| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验 | 欧美日韩黄片免| 国产91精品成人一区二区三区| 国产精品久久久久久av不卡| 精品99又大又爽又粗少妇毛片 | 日本欧美国产在线视频| a级一级毛片免费在线观看| 1000部很黄的大片| 国产精品一区二区三区四区久久| 狂野欧美白嫩少妇大欣赏| 国产精品野战在线观看| 欧美日本亚洲视频在线播放| 校园人妻丝袜中文字幕| 欧美中文日本在线观看视频| h日本视频在线播放| 久久久成人免费电影| 啪啪无遮挡十八禁网站| 欧美另类亚洲清纯唯美| 欧美日韩综合久久久久久 | 91精品国产九色| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av香蕉五月| 91久久精品国产一区二区成人| 精品久久久久久成人av| 国产视频一区二区在线看| 日韩欧美一区二区三区在线观看| 日本欧美国产在线视频| 免费av观看视频| 神马国产精品三级电影在线观看| 日本爱情动作片www.在线观看 | 可以在线观看的亚洲视频| 精品久久久久久成人av| 日韩一区二区视频免费看| 人人妻人人看人人澡| 在线观看66精品国产| 午夜福利在线在线| 草草在线视频免费看| 亚洲内射少妇av| 制服丝袜大香蕉在线| 久久九九热精品免费| 成人av在线播放网站| 欧美丝袜亚洲另类 | 成人永久免费在线观看视频| 久久久久免费精品人妻一区二区| 亚洲图色成人| 美女被艹到高潮喷水动态| 国产精品98久久久久久宅男小说| ponron亚洲| 午夜a级毛片| 久久久久久大精品| 亚洲精品影视一区二区三区av| 非洲黑人性xxxx精品又粗又长| 久久国产精品人妻蜜桃| 欧美激情久久久久久爽电影| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久国产a免费观看| 久久九九热精品免费| 免费人成在线观看视频色| 日本免费一区二区三区高清不卡| 国产一区二区亚洲精品在线观看| 中出人妻视频一区二区| 午夜福利在线观看免费完整高清在 | 全区人妻精品视频| 悠悠久久av| 国产精品三级大全| 国产视频内射| 国产视频一区二区在线看| 国产精品人妻久久久影院| 麻豆国产av国片精品| 人人妻人人看人人澡| 国产午夜精品久久久久久一区二区三区 | 国产高清视频在线观看网站| 麻豆精品久久久久久蜜桃| 欧美绝顶高潮抽搐喷水| 在线观看一区二区三区| 亚洲欧美日韩高清专用| 伦精品一区二区三区| 欧美日韩精品成人综合77777| 夜夜爽天天搞| 精品久久久久久久久久免费视频| 国产精品福利在线免费观看| 最后的刺客免费高清国语| 99热精品在线国产| 日韩中字成人| 男人和女人高潮做爰伦理| 免费一级毛片在线播放高清视频| 亚洲精品一区av在线观看| 最近视频中文字幕2019在线8| 久久草成人影院| 国产精品一区二区三区四区免费观看 | 免费不卡的大黄色大毛片视频在线观看 | 真人一进一出gif抽搐免费| 亚洲中文日韩欧美视频| 欧美bdsm另类| 日本一二三区视频观看| 伦精品一区二区三区| 九九久久精品国产亚洲av麻豆| 国产亚洲精品久久久com| 可以在线观看毛片的网站| 午夜影院日韩av| 毛片一级片免费看久久久久 | 成人国产麻豆网| 禁无遮挡网站| 欧美人与善性xxx| 最新中文字幕久久久久| 简卡轻食公司| 欧美日韩瑟瑟在线播放| 88av欧美| 中出人妻视频一区二区| 99久久久亚洲精品蜜臀av| 我要看日韩黄色一级片| 中文在线观看免费www的网站| 身体一侧抽搐| 欧美性感艳星| 高清在线国产一区| 亚洲真实伦在线观看| 亚洲aⅴ乱码一区二区在线播放| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| av视频在线观看入口| 最近最新中文字幕大全电影3| 亚洲专区中文字幕在线| 色精品久久人妻99蜜桃| 中国美白少妇内射xxxbb| 九九久久精品国产亚洲av麻豆| 久久久国产成人免费| 色5月婷婷丁香| 久久久久久伊人网av| 国产精品一区www在线观看 | 精品人妻熟女av久视频| 亚洲av二区三区四区| 国产亚洲91精品色在线| 老司机福利观看| 动漫黄色视频在线观看| 国内毛片毛片毛片毛片毛片| 成人综合一区亚洲| 久久精品久久久久久噜噜老黄 | 欧美一区二区精品小视频在线| 美女xxoo啪啪120秒动态图| 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 色精品久久人妻99蜜桃| 欧美潮喷喷水| 99久久精品热视频| 男人狂女人下面高潮的视频| 欧美一区二区精品小视频在线| 日日干狠狠操夜夜爽| 久久久精品大字幕| 色噜噜av男人的天堂激情| 1024手机看黄色片| 亚洲专区国产一区二区| 国产精品精品国产色婷婷| 直男gayav资源| 老熟妇仑乱视频hdxx| 久久久久久久精品吃奶| 精品久久久久久久人妻蜜臀av| avwww免费| 好男人在线观看高清免费视频| 欧美+日韩+精品| 一级a爱片免费观看的视频| 欧美zozozo另类| 99热这里只有精品一区| 精品福利观看| 97超视频在线观看视频| 久9热在线精品视频| 毛片一级片免费看久久久久 | 蜜桃亚洲精品一区二区三区| 内射极品少妇av片p| 一区福利在线观看| 一夜夜www| 99久国产av精品| 伊人久久精品亚洲午夜| 少妇丰满av| 91久久精品国产一区二区成人| 女人被狂操c到高潮| 国模一区二区三区四区视频| 少妇裸体淫交视频免费看高清| 一进一出抽搐动态| 国产伦精品一区二区三区四那| 91久久精品电影网| 国产欧美日韩精品亚洲av| 美女cb高潮喷水在线观看| 中文字幕免费在线视频6| 日韩国内少妇激情av| 午夜福利在线在线| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 久久亚洲真实| 成人亚洲精品av一区二区| 网址你懂的国产日韩在线| 一级a爱片免费观看的视频| 欧美激情在线99| 欧美高清性xxxxhd video| 99国产精品一区二区蜜桃av| 欧美日韩精品成人综合77777| 国产精品女同一区二区软件 | 国产色婷婷99| 午夜福利高清视频| 五月玫瑰六月丁香| 日本免费a在线| 51国产日韩欧美| 日韩欧美国产在线观看| 麻豆成人av在线观看| 一进一出好大好爽视频| 99热6这里只有精品| 国产亚洲精品综合一区在线观看| 免费电影在线观看免费观看| 欧美色欧美亚洲另类二区| 999久久久精品免费观看国产| 亚洲七黄色美女视频| 看片在线看免费视频| 亚洲美女视频黄频| 女人被狂操c到高潮| 日韩欧美 国产精品| 色精品久久人妻99蜜桃| 国产精品一区二区三区四区免费观看 | 国产蜜桃级精品一区二区三区| 美女 人体艺术 gogo| 国产老妇女一区| 午夜亚洲福利在线播放| 欧美高清成人免费视频www| 国产精品免费一区二区三区在线| 校园春色视频在线观看| 国产精品无大码| 熟女人妻精品中文字幕| 国产免费一级a男人的天堂| 99久久无色码亚洲精品果冻| 看免费成人av毛片| 少妇裸体淫交视频免费看高清| 国产av麻豆久久久久久久| 亚洲综合色惰| 日本在线视频免费播放| 午夜a级毛片| 成年免费大片在线观看| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 成人国产综合亚洲| 狠狠狠狠99中文字幕| a级一级毛片免费在线观看| 国产精品一及| 国产在视频线在精品| 亚洲自拍偷在线| 日韩 亚洲 欧美在线| 欧美三级亚洲精品| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 精品免费久久久久久久清纯| 91在线观看av| 色综合站精品国产| 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 别揉我奶头 嗯啊视频| 国产精品无大码| 亚洲人与动物交配视频| 亚洲第一区二区三区不卡| 久久亚洲真实| 国产乱人伦免费视频| 在线观看美女被高潮喷水网站| 国产精品电影一区二区三区| 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 五月玫瑰六月丁香| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 九九爱精品视频在线观看| 91av网一区二区| 日本撒尿小便嘘嘘汇集6| 午夜福利高清视频| 亚洲人与动物交配视频| 久久精品人妻少妇| 日韩 亚洲 欧美在线| 欧美日韩黄片免| 久久久久国内视频| 成年版毛片免费区| 夜夜爽天天搞| 少妇高潮的动态图| 一个人观看的视频www高清免费观看| 国产精品99久久久久久久久| 深爱激情五月婷婷| 国产欧美日韩精品一区二区| 99在线人妻在线中文字幕| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 精品久久久久久久久久久久久| 热99在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| 校园春色视频在线观看| 热99在线观看视频| 久久精品国产鲁丝片午夜精品 | 国产在线精品亚洲第一网站| 国产精品人妻久久久影院| 国产又黄又爽又无遮挡在线| 国产精品国产三级国产av玫瑰| 日本三级黄在线观看| 18禁黄网站禁片免费观看直播| 性欧美人与动物交配| 国产三级在线视频| 免费看av在线观看网站| 日韩中字成人| 亚洲久久久久久中文字幕| 日本色播在线视频| 亚洲久久久久久中文字幕| 精华霜和精华液先用哪个| 久久精品国产99精品国产亚洲性色| av在线观看视频网站免费| 非洲黑人性xxxx精品又粗又长| 亚洲专区中文字幕在线| 一级黄片播放器| 嫩草影视91久久| 热99re8久久精品国产| 日韩人妻高清精品专区| 免费在线观看成人毛片| 午夜影院日韩av| 日本黄色片子视频| 亚洲成人中文字幕在线播放| 欧美另类亚洲清纯唯美| 国产免费av片在线观看野外av| 一卡2卡三卡四卡精品乱码亚洲| 搡老岳熟女国产| 亚洲精品456在线播放app | 九色国产91popny在线| 国产私拍福利视频在线观看| 国产成人a区在线观看| 国产高清视频在线观看网站| 国产欧美日韩精品一区二区| 嫩草影院精品99| 色综合婷婷激情| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| a级毛片a级免费在线| 午夜视频国产福利| 我要搜黄色片| 伦理电影大哥的女人| 欧美zozozo另类| 欧美不卡视频在线免费观看| 一边摸一边抽搐一进一小说| 久久精品91蜜桃| 禁无遮挡网站| 午夜福利高清视频| 麻豆av噜噜一区二区三区| av中文乱码字幕在线| 桃色一区二区三区在线观看| 男人舔奶头视频| 麻豆成人午夜福利视频| 免费av不卡在线播放| 欧美日韩中文字幕国产精品一区二区三区| 99久国产av精品| 狂野欧美白嫩少妇大欣赏| 麻豆久久精品国产亚洲av| 乱系列少妇在线播放| 97碰自拍视频| 日本黄色片子视频| 一区福利在线观看| 国产又黄又爽又无遮挡在线| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添av毛片 | 亚洲精品一卡2卡三卡4卡5卡| 99久久无色码亚洲精品果冻| 中文亚洲av片在线观看爽| 国产精品嫩草影院av在线观看 | 春色校园在线视频观看| 人妻夜夜爽99麻豆av| 一进一出抽搐动态| 国产精品无大码| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女视频黄频| 国产av一区在线观看免费| 国产精品一及| 国产欧美日韩精品一区二区| 超碰av人人做人人爽久久| 欧美精品啪啪一区二区三区| 中国美女看黄片| 国产精品久久久久久久电影| 久99久视频精品免费| 日本五十路高清| 最新在线观看一区二区三区| 日本 av在线| 乱码一卡2卡4卡精品| 国产男人的电影天堂91| 国产在线男女| 亚洲av第一区精品v没综合| 亚洲色图av天堂| 一区二区三区免费毛片| 日韩欧美 国产精品| 国产中年淑女户外野战色| 最好的美女福利视频网| 变态另类成人亚洲欧美熟女| 性色avwww在线观看| 国产一级毛片七仙女欲春2| 精品人妻熟女av久视频| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 亚州av有码| 99久久精品一区二区三区| av在线天堂中文字幕| 人妻制服诱惑在线中文字幕| 欧美xxxx黑人xx丫x性爽| 亚洲中文字幕日韩| 欧美xxxx性猛交bbbb| 精品久久久久久,| 免费黄网站久久成人精品| 国产精品,欧美在线| 啦啦啦观看免费观看视频高清| 久久久久久大精品| 99久久九九国产精品国产免费| www.www免费av| 悠悠久久av| 99热这里只有是精品在线观看| 3wmmmm亚洲av在线观看| 观看免费一级毛片| 亚洲精品国产成人久久av| 内射极品少妇av片p| 精品久久久久久久末码| ponron亚洲| 干丝袜人妻中文字幕| 草草在线视频免费看| 日本一本二区三区精品| 中文字幕熟女人妻在线| 婷婷亚洲欧美| 国产成人影院久久av| 一区二区三区高清视频在线| 免费看av在线观看网站| 99久久中文字幕三级久久日本| 99久久精品一区二区三区| 国产伦人伦偷精品视频| 在线观看免费视频日本深夜| 51国产日韩欧美| 3wmmmm亚洲av在线观看| 黄色欧美视频在线观看| 美女xxoo啪啪120秒动态图| 白带黄色成豆腐渣| 1024手机看黄色片| 国产 一区精品| 天堂影院成人在线观看| 国产精品爽爽va在线观看网站| 人妻制服诱惑在线中文字幕| 亚洲最大成人手机在线| 无遮挡黄片免费观看| 直男gayav资源| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 国内精品宾馆在线| 国产一区二区三区av在线 | 成人鲁丝片一二三区免费| 久久99热这里只有精品18| 国产欧美日韩精品一区二区| 日日干狠狠操夜夜爽| 亚洲国产高清在线一区二区三| 亚洲欧美激情综合另类| 久久精品夜夜夜夜夜久久蜜豆| 又黄又爽又刺激的免费视频.| 琪琪午夜伦伦电影理论片6080| 日本熟妇午夜| 亚洲欧美日韩高清专用| 亚洲欧美日韩无卡精品| 国产 一区 欧美 日韩| 久久久久久久午夜电影| 国内精品一区二区在线观看| 真人一进一出gif抽搐免费| 久久人人精品亚洲av| 亚洲第一电影网av| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品合色在线| 99在线视频只有这里精品首页| 无遮挡黄片免费观看| 国产高清三级在线| 国产精品女同一区二区软件 | 久久精品国产亚洲网站| 成人国产麻豆网| 99久久精品一区二区三区| 两个人视频免费观看高清| av.在线天堂| 麻豆精品久久久久久蜜桃| 亚洲av.av天堂| 免费av观看视频| 国模一区二区三区四区视频| 日日啪夜夜撸| 一本精品99久久精品77| 黄色视频,在线免费观看| 熟女电影av网| 亚洲欧美激情综合另类| 色吧在线观看| 欧美日本视频| 久久天躁狠狠躁夜夜2o2o| 日韩一本色道免费dvd| 欧美另类亚洲清纯唯美| 国国产精品蜜臀av免费| 成年版毛片免费区| 一进一出抽搐动态| 久久99热这里只有精品18| 成人精品一区二区免费| 亚洲精品久久国产高清桃花| 午夜福利高清视频| 日韩欧美在线二视频| a级毛片a级免费在线| 亚洲专区国产一区二区| 久久99热6这里只有精品| 美女黄网站色视频| 91狼人影院| www.www免费av| 久久人人精品亚洲av| 国产 一区 欧美 日韩| 少妇的逼好多水| 在线国产一区二区在线| 男女之事视频高清在线观看| 永久网站在线| 久久久国产成人精品二区| 成人国产麻豆网| 国产 一区 欧美 日韩| 91久久精品国产一区二区三区| 在线观看免费视频日本深夜| 麻豆av噜噜一区二区三区| 国产精品一区二区三区四区久久| 日韩,欧美,国产一区二区三区 | 亚洲成人精品中文字幕电影| 亚洲人成网站高清观看| 狠狠狠狠99中文字幕| 国产爱豆传媒在线观看| 欧美高清性xxxxhd video| 九九热线精品视视频播放| 欧美最新免费一区二区三区| 99热这里只有是精品在线观看| 深夜精品福利| 不卡视频在线观看欧美| 人妻少妇偷人精品九色| 最新中文字幕久久久久| 免费av观看视频| 美女cb高潮喷水在线观看| 亚洲精品成人久久久久久| 午夜精品一区二区三区免费看| 校园人妻丝袜中文字幕| 久久热精品热| 日韩欧美国产在线观看| 日本爱情动作片www.在线观看 | 综合色av麻豆| 少妇高潮的动态图| aaaaa片日本免费| 深夜精品福利| 免费无遮挡裸体视频| 很黄的视频免费| 天堂动漫精品| 国产乱人伦免费视频| 久久久久性生活片| 99久久中文字幕三级久久日本| 国产探花极品一区二区| 别揉我奶头~嗯~啊~动态视频| 久久这里只有精品中国| 成人永久免费在线观看视频| 国产淫片久久久久久久久| 国产精品久久久久久精品电影| 国产伦一二天堂av在线观看| 成年人黄色毛片网站| 色噜噜av男人的天堂激情| 亚洲国产精品成人综合色| 又粗又爽又猛毛片免费看| 久久久久久久亚洲中文字幕| 一个人观看的视频www高清免费观看| 免费看a级黄色片| 日日夜夜操网爽| 中亚洲国语对白在线视频|