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

    基于PCA法的巖體震源破壞面產(chǎn)狀分析及工程應(yīng)用

    2021-10-20 09:55:02劉胤池李庶林周夢婧
    金屬礦山 2021年9期
    關(guān)鍵詞:產(chǎn)狀微震橢球

    劉胤池 李庶林 周夢婧 唐 超

    (廈門大學(xué)建筑與土木工程學(xué)院,福建廈門361005)

    近十多年以來,微震監(jiān)測技術(shù)被廣泛用于我國礦山生產(chǎn)、水電工程和隧道工程等安全監(jiān)測領(lǐng)域,在巖體災(zāi)害預(yù)警方面取得了一定的應(yīng)用效果[1-3]。研究工程巖體的震源機(jī)制,可以深刻認(rèn)識工程巖體的細(xì)觀破裂機(jī)理和宏觀破壞模式,對于災(zāi)害的監(jiān)測預(yù)警具有重要意義。震源機(jī)制涉及到巖體破裂源的物理力學(xué)機(jī)制,是微震監(jiān)測領(lǐng)域的一個理論核心問題,也是目前微震研究的一個熱點(diǎn)。破壞面產(chǎn)狀分析是震源機(jī)制解中很重要的一部分,常利用矩張量反演或應(yīng)力場反演來求得破壞面產(chǎn)狀[4-5]。

    針對求解構(gòu)造地震破壞面產(chǎn)狀的問題,MICHELINI等[6]首次提出了一種快速獲取破壞面產(chǎn)狀的方法,將主成分分析法運(yùn)用至地震學(xué)中,提出了“局部破裂橢球體”的概念,并根據(jù)震源發(fā)生的時(shí)間順序及空間分布推斷復(fù)雜地區(qū)微震活動的平均破壞面幾何形狀以及活動趨勢;POSADAS等[7]進(jìn)一步完善了該方法,補(bǔ)充了空間位置的主成分分析,用這兩種方法可以求得地震導(dǎo)致的主破壞面產(chǎn)狀以及破壞面的時(shí)空演化趨勢。上述方法在地震學(xué)中得到了較為廣泛的應(yīng)用,得到國內(nèi)外學(xué)者的進(jìn)一步深入分析、發(fā)展和運(yùn)用[8-11]。

    巖石工程中的巖體破裂與地震類似,在宏觀破壞面最終形成之前會有大量的微震事件產(chǎn)生,且這些微震事件具有空間和時(shí)間上的聚類特征,微震事件在時(shí)間以及空間上的分布可以提供關(guān)于破壞面空間產(chǎn)狀的相關(guān)信息。URBANCIC等[12]首次在微震監(jiān)測中利用主成分分析法快速求解破壞面產(chǎn)狀,其結(jié)果與原位測量結(jié)果相差較小,證明了該方法在礦山定位中的實(shí)用性;TRIFU等[13]利用多種方法對礦山的應(yīng)力狀態(tài)進(jìn)行分析,將主成分分析法得到的破壞面構(gòu)造與P、T軸方向的應(yīng)力反演結(jié)果對比,結(jié)果相近;COULSON[14]等運(yùn)用時(shí)間主成分技術(shù)分析煤柱屈服區(qū)的范圍和趨勢,通過煤柱中形成斷裂結(jié)構(gòu)的事件簇來表征破壞過程;GAMMANS[15]引入平面度與線性度兩個概念來描述局部破裂橢球體,并測試了PCA法求解破壞面產(chǎn)狀的穩(wěn)定性。

    在地震學(xué)中利用PCA法快速求解地震斷層面產(chǎn)狀信息及演化趨勢的應(yīng)用較多,然而在巖石工程中,僅利用PCA法獲取破壞面產(chǎn)狀信息,在獲取破壞面產(chǎn)狀演化趨勢方面鮮有研究。為此,本研究將PCA方法應(yīng)用至柿竹園礦一次典型的特大山體剪切滑移案例中,根據(jù)接收到的微震定位數(shù)據(jù)來獲取滑移破壞面的相關(guān)信息及演化趨勢,以進(jìn)一步探索PCA方法在巖體剪滑破壞模式中的應(yīng)用效果。

    1 主成分分析原理

    主成分分析(PCA)是一種量化數(shù)據(jù)分布的統(tǒng)計(jì)方法,利用正交變換將一系列可能線性相關(guān)的變量轉(zhuǎn)換為一組線性不相關(guān)的新變量(也稱為主成分),利用新變量展示數(shù)據(jù)的特征[16]。在工程微震定位事件分析中,對定位事件位置應(yīng)用PCA,可量化微震事件聚簇類的形狀(線、平面、球體)和方向的聚集程度。震源定位位置是在一個三維X-Y-Z坐標(biāo)系中進(jìn)行標(biāo)記的,可用一個最佳擬合橢球體來描述定位事件簇的形狀。定位事件在空間與時(shí)間上的分布會提供有關(guān)破壞面的幾何形狀信息,本研究從空間與時(shí)間兩個方面來對巖體破壞面進(jìn)行分析?;谡鹪纯臻g分布的主成分分析(PCA With a Spatial Window,SPCA)是根據(jù)微震信號的空間坐標(biāo)位置利用PCA獲取微震信號的空間分布特征,再進(jìn)一步獲得礦山的主要破壞面的產(chǎn)狀;基于震源時(shí)間分布的主成分分析(PCA With a Temporal Window,TPCA)將一個覆蓋有固定數(shù)量事件的時(shí)間窗口在按時(shí)間順序排列好的微震事件群上滑動,求得每個窗口代表的橢球體并將方向不同的橢球體進(jìn)行分離,從而能夠獲得不同時(shí)刻的破壞面產(chǎn)狀,利用該方法可以獲得破壞面的演化趨勢。

    1.1 基于震源空間分布的主成分分析

    對于空間一定區(qū)域內(nèi)的N個微震定位事件,利用主元分析方法對震源破裂面進(jìn)行分析,該事件群的協(xié)方差矩陣用S表示

    該矩陣在幾何學(xué)上可表示為由震源坐標(biāo)擬合的空間橢球,其特征向量方向給出了橢球軸的空間方向,特征值為橢球體半軸的長度,該橢球體稱為局部破裂橢球體[6,15]。

    橢球的橢圓度為e,公式為

    當(dāng)e≥2.5時(shí),橢球體可以看作是微震事件作用下所產(chǎn)生的破裂面;當(dāng)e<2.5時(shí),此時(shí)的橢球體不能視為平面,而是體積形態(tài),在這種情況下不能根據(jù)橢球體的特征來推斷破壞面的產(chǎn)狀。

    對于可以視為平面的橢球體,第1特征向量n1和第2特征向量n2構(gòu)成的平面即為潛在破壞面,第3特征向量n3為破壞面法向。根據(jù)破壞面法向中的3個元素n3(1)、n3(2)和n3(3),可以求得破壞面得出走向φ和傾角δ,公式為[17]

    1.2 基于震源時(shí)間分布的主成分分析

    為了推斷巖體潛在破壞面在時(shí)間序列上的演化進(jìn)程,本研究采用基于震源位置時(shí)間分布的主成分分析法對監(jiān)測到的微震定位數(shù)據(jù)進(jìn)行分析。

    對于N個震源定位事件,將其按時(shí)間順序排列,選定一個長為K的滑動窗口,K為一次性覆蓋微震數(shù)據(jù)的個數(shù),此時(shí)每K個事件可以看作一個橢球體,在N個事件上共有p=N-K+1個橢球體。為了減少窗口長度對結(jié)果的影響,需選擇最佳窗口長度,當(dāng)Q接近于常數(shù)時(shí),此時(shí)的K值為最佳窗口長度[6-7]:

    將窗口在震源序列上滑動,每個窗口為一個橢球體,橢球體之間主軸的相差角度可進(jìn)行如下計(jì)算:

    式中,α1h、β1h分別為第1個橢球體與第h個橢球體主軸相差的角度,(°);、分別為第h個橢球體的最大特征向量與最小特征向量。

    設(shè)定分離角度分別為δα與δβ,當(dāng)滿足下列條件時(shí),說明兩橢球體主軸方向幾乎相同,為分離出主軸不同的橢球體,窗口繼續(xù)滑動直至式(7)中一個不等式不成立為止,即找出不同主軸方向的橢球體。

    當(dāng)上式不滿足時(shí),取該時(shí)間段所有橢球體主軸方向向量的平均值作為該時(shí)刻橢球體主軸方向,分別用v1(K ,r)和v3(K ,s)表示:

    式中,r和s為滿足下列不等式的連續(xù)橢球體個數(shù)。在工程中的微震中,由于事件與事件間相距較近,因此δα與δβ需選擇比地震學(xué)中更為精確的角度來區(qū)分不同的橢球體,本研究兩者取值分別為5°和10°。

    本研究采用主成分分析法獲取破壞面產(chǎn)狀及破壞面演化進(jìn)程信息,計(jì)算流程如圖1所示。

    在兩種分析方法中,均按時(shí)間順序?qū)⒄鹪醋鴺?biāo)進(jìn)行排列,根據(jù)SPCA可以獲得局部破裂橢球體的主軸信息,據(jù)此可進(jìn)一步求得破壞面產(chǎn)狀;在TPCA中,利用式(5)確定合適的滑動窗口長度K,是為了降低因滑動窗口長度不同對破壞面產(chǎn)狀的影響[6],之后求解每個窗口的局部橢球體信息與SPCA類似,根據(jù)式(6)和式(7)可以分離出主軸方向不同的橢球體,從而獲得不同時(shí)刻的破壞面產(chǎn)狀。

    2 微震監(jiān)測案例分析

    2.1 工程概況

    為了監(jiān)測采空區(qū)、礦柱和上覆厚大巖層的穩(wěn)定性和安全性,柿竹園有色金屬有限責(zé)任公司所屬的多金屬采礦場于2010年引入加拿大ESG的36通道的全數(shù)字微震監(jiān)測系統(tǒng)。該系統(tǒng)全部采用單軸加速度傳感器,分別布置在514 m平面、558 m平面和630 m平面內(nèi)。2012年6月21日該礦進(jìn)行的一次大爆破落礦過程中,直接誘發(fā)了采場上覆巖層(山體)的特大規(guī)?;菩涂逅F(xiàn)象。微震監(jiān)測系統(tǒng)完整地監(jiān)測到了這次破壞的孕育、垮塌的全過程[18]。本次大爆破誘發(fā)的微震定位事件總共為85個,圖2中給出了山體垮塌區(qū)域的微震定位事件分布,主要集中于上部懸頂破壞區(qū)域的640~755 m高程范圍內(nèi)。這是一次非常典型的、由大爆破誘發(fā)的特大山體破壞案例。本研究根據(jù)記錄到的微震定位數(shù)據(jù),利用SPCA法和TPCA法,對礦山破壞面的相關(guān)信息進(jìn)行分析。

    2.2 案例分析

    根據(jù)1.1節(jié)和1.2節(jié)描述,利用SPCA法與TPCA法能夠快速計(jì)算出微震活動產(chǎn)生的破壞面產(chǎn)狀,并獲得破壞面的發(fā)展過程。將SPCA法與TPCA法應(yīng)用于2012年6月21日一次因大爆破產(chǎn)生的特大山體剪切滑移事件中,推斷可能產(chǎn)生的潛在破壞面相關(guān)信息。本研究選取上部懸頂巖體滑移定位事件進(jìn)行分析,去除定位誤差較大的事件后,最終選取52個定位事件(表1)進(jìn)行SPCA與TPCA分析。

    2.2.1 SPCA法求解破壞面產(chǎn)狀

    將微震監(jiān)測系統(tǒng)監(jiān)測到52個巖體破裂的震源定位事件的坐標(biāo)按時(shí)間順序進(jìn)行記錄,利用Matlab中的z-score函數(shù)對原始空間坐標(biāo)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理,結(jié)果見表2。

    通過對協(xié)方差矩陣進(jìn)行特征值-特征向量分解,并對求出的特征值進(jìn)行排序。由于前兩個主元向量是構(gòu)成所求潛在破裂平面的基向量,最后一個向量為該平面的法向量,因而使用該法向矢量和位于平面坐標(biāo)的均值擬合到平面上,計(jì)算出可能的破壞面產(chǎn)狀。

    案例中微震定位數(shù)據(jù)的協(xié)方差矩陣為

    求得的特征值從大到小排列分別為λ1=2.484 8,λ2=0.435 6,λ3=0.079 6,其對應(yīng)的特征向量分別為n1=(-0.610 8,-0.589 8,0.528 3),n2=(0.255 2,0.485 0,0.836 4),n3=(0.749 6,-0.645 7,0.145 7)。

    根據(jù)1.1節(jié)描述,因?yàn)棣?/λ3≥ 2.5,即可將該橢球體從幾何特征上視為一個平面,該平面即為所求的可能的震源破裂面由n1與n2兩軸確定,n3為破壞面法向。

    根據(jù)式(4)及傾向與走向之間的關(guān)系,可求得破壞面的傾向及傾角分別為319.26°與81.62°。在Wulff投影下,用SPCA法求出的主破壞面產(chǎn)狀由圖3所示;實(shí)際測量的破裂面產(chǎn)狀傾向?yàn)?27°和傾角為78°,如圖4和圖5所示。可見,用SPCA法可以較為快速且精準(zhǔn)地獲得破壞面產(chǎn)狀。

    2.2.2 TPCA法求解破壞面產(chǎn)狀演化趨勢

    首先通過Q值的變化對K取值,如圖6所示。當(dāng)K=36時(shí),Q值之后的變化較為緩慢,故選擇此時(shí)的K值作為滑動窗口長度。

    在求得每個窗口所代表的橢球體的主軸方向向量后,利用式(6)和式(7)將不同的橢球體分離出,根據(jù)式(8)求得不同橢球體的主軸方向,據(jù)此可以得到不同時(shí)刻的不同橢球體的破壞面產(chǎn)狀,從而獲得破壞面的演化趨勢。

    利用TPCA法獲得的破壞面變化信息見表3。

    由表3可知:此次滑移事故中,主破壞面傾向從NW316.14°向NW317.47°方向發(fā)展,破壞面傾角由87.96°向85.05°發(fā)展。

    3 結(jié)論

    本研究分析了基于主成分分析快速獲取主破壞面相關(guān)信息的方法,該方法分為SPCA和TPCA兩部分。根據(jù)震源坐標(biāo)的時(shí)間與空間分布,利用SPCA法可以獲得局部破裂橢球體的主軸信息,據(jù)此可進(jìn)一步求得破壞面產(chǎn)狀;利用TPCA法可以分離出主軸方向不同的橢球體,從而可獲得不同時(shí)刻破壞面的產(chǎn)狀,來判斷破壞面的演化趨勢。將該方法應(yīng)用于柿竹園礦山的一次山體滑移事件分析中,快速獲取了破壞面的產(chǎn)狀及演化趨勢,得到了如下結(jié)論:

    (1)由SPCA法確定了此次山體滑移產(chǎn)生的破壞面產(chǎn)狀,其傾向?yàn)镹W319.26°,傾角為81.62°,與實(shí)際測量結(jié)果(傾向?yàn)镹W327°,傾角為78°)較接近,證明了該方法的合理性,此外,SPCA法利用了微震信號的定位信息,因而計(jì)算過程較簡便。

    (2)利用TPCA法分析了破壞面的演化過程,破壞面產(chǎn)狀與SPCA法計(jì)算結(jié)果相近,其傾向變化方向由NW316.14°向NW317.47°方向發(fā)展,傾角由 87.96°向85.05°發(fā)展,可推斷該破壞傾向向北發(fā)展,且傾角逐漸減小,可見,TPCA法夠推斷出破壞面的演化趨勢,為災(zāi)害預(yù)防提供理論支撐。

    (3)當(dāng)定位數(shù)據(jù)充足時(shí),在求得最佳的窗口長度K后可以進(jìn)一步考慮最佳步距問題(本研究默認(rèn)為1),在保證結(jié)果準(zhǔn)確性的前提下提高分析效率。

    猜你喜歡
    產(chǎn)狀微震橢球
    基于產(chǎn)狀Fisher分布的巖石節(jié)理最小樣本容量確定方法
    獨(dú)立坐標(biāo)系橢球變換與坐標(biāo)換算
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    淺談KJ768煤礦微震監(jiān)測系統(tǒng)的應(yīng)用
    新疆鋼鐵(2021年1期)2021-10-14 08:45:48
    淺談砂巖儲層的巖石學(xué)特征
    激電聯(lián)合剖面在判斷矽卡巖型礦床礦體產(chǎn)狀中的應(yīng)用
    長平煤業(yè)5302 綜放工作面頂板巖層移動規(guī)律研究
    魅力中國(2020年46期)2020-02-02 11:50:08
    赤平投影法在邊坡穩(wěn)定性分析中的應(yīng)用
    橢球精加工軌跡及程序設(shè)計(jì)
    基于外定界橢球集員估計(jì)的純方位目標(biāo)跟蹤
    色视频www国产| 国产精品一区二区三区四区久久| 亚洲三级黄色毛片| 狠狠精品人妻久久久久久综合| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 久久精品人妻少妇| 人人妻人人看人人澡| 国产成人a区在线观看| 久久久色成人| 色综合站精品国产| 亚洲国产av新网站| 精品熟女少妇av免费看| 国产精品1区2区在线观看.| av国产久精品久网站免费入址| 精品一区二区免费观看| 色网站视频免费| 亚洲最大成人中文| 久久久a久久爽久久v久久| 亚洲国产精品成人综合色| 久久久久久久久久黄片| 亚洲欧美精品自产自拍| 国产亚洲av片在线观看秒播厂 | 成人一区二区视频在线观看| 欧美日韩国产mv在线观看视频 | 乱系列少妇在线播放| 国产免费一级a男人的天堂| 男女那种视频在线观看| 亚洲国产欧美人成| 国产精品久久久久久精品电影| 草草在线视频免费看| 亚洲av免费在线观看| 男女边吃奶边做爰视频| 久久久久久国产a免费观看| 男女边吃奶边做爰视频| 搡老乐熟女国产| 一级二级三级毛片免费看| 亚洲av.av天堂| 国产一级毛片七仙女欲春2| 久久久久久久久久久丰满| 亚洲av免费高清在线观看| 国产黄a三级三级三级人| 亚洲熟女精品中文字幕| 国产熟女欧美一区二区| 久久这里有精品视频免费| 国产一级毛片在线| 能在线免费看毛片的网站| 九色成人免费人妻av| 汤姆久久久久久久影院中文字幕 | 你懂的网址亚洲精品在线观看| 国产伦精品一区二区三区四那| 国产精品久久久久久久久免| 亚洲婷婷狠狠爱综合网| 青春草亚洲视频在线观看| 亚洲av男天堂| av在线观看视频网站免费| 三级男女做爰猛烈吃奶摸视频| 美女被艹到高潮喷水动态| 亚洲国产色片| 国产精品爽爽va在线观看网站| 久久久久精品久久久久真实原创| 国产成人精品久久久久久| 日韩成人伦理影院| 丝瓜视频免费看黄片| 久久久色成人| av在线天堂中文字幕| 国产精品一及| 91午夜精品亚洲一区二区三区| 免费无遮挡裸体视频| 亚洲在久久综合| 欧美激情国产日韩精品一区| 三级毛片av免费| 久久韩国三级中文字幕| 麻豆国产97在线/欧美| 只有这里有精品99| 人体艺术视频欧美日本| 欧美人与善性xxx| 夜夜爽夜夜爽视频| 亚洲精品国产av蜜桃| 熟女人妻精品中文字幕| 亚洲国产精品成人综合色| 中国国产av一级| 亚洲一级一片aⅴ在线观看| 日日啪夜夜撸| 成人亚洲精品一区在线观看 | 黄色一级大片看看| 亚洲熟妇中文字幕五十中出| 国产一级毛片七仙女欲春2| 亚洲无线观看免费| 精品人妻熟女av久视频| 国产黄色小视频在线观看| 最近最新中文字幕免费大全7| 亚洲成人精品中文字幕电影| 免费观看在线日韩| 亚洲最大成人中文| 欧美日本视频| 中文在线观看免费www的网站| 97精品久久久久久久久久精品| 能在线免费看毛片的网站| 一个人免费在线观看电影| 丰满乱子伦码专区| 欧美日韩精品成人综合77777| 亚洲av国产av综合av卡| 女人久久www免费人成看片| 国产成人一区二区在线| videos熟女内射| 女的被弄到高潮叫床怎么办| 亚洲欧美中文字幕日韩二区| 久久97久久精品| 久久精品综合一区二区三区| 九九久久精品国产亚洲av麻豆| 只有这里有精品99| 国内少妇人妻偷人精品xxx网站| av天堂中文字幕网| 黄色一级大片看看| 亚洲国产日韩欧美精品在线观看| 亚洲在线自拍视频| 成年av动漫网址| 中文字幕亚洲精品专区| 麻豆乱淫一区二区| 插阴视频在线观看视频| 听说在线观看完整版免费高清| 色吧在线观看| 精华霜和精华液先用哪个| 欧美精品一区二区大全| 久久久久久久久久久免费av| av在线亚洲专区| 亚洲无线观看免费| 可以在线观看毛片的网站| 久久久a久久爽久久v久久| 国产乱人偷精品视频| 久久99热6这里只有精品| 国产69精品久久久久777片| 国产美女午夜福利| av女优亚洲男人天堂| 国产精品嫩草影院av在线观看| 日本欧美国产在线视频| 精品久久久久久久末码| 国产黄片视频在线免费观看| 亚洲精品国产成人久久av| 久久人人爽人人爽人人片va| 97超视频在线观看视频| 亚洲18禁久久av| 色吧在线观看| 七月丁香在线播放| 亚洲国产av新网站| 3wmmmm亚洲av在线观看| 亚洲精品久久久久久婷婷小说| 日韩欧美 国产精品| 中文天堂在线官网| 国产爱豆传媒在线观看| 国产午夜福利久久久久久| 插逼视频在线观看| 国产精品久久久久久久电影| 高清午夜精品一区二区三区| 国产男女超爽视频在线观看| 国产男女超爽视频在线观看| 亚洲av二区三区四区| 女人久久www免费人成看片| 插逼视频在线观看| 成年女人看的毛片在线观看| 日韩av在线免费看完整版不卡| 26uuu在线亚洲综合色| av天堂中文字幕网| 日韩电影二区| 国产探花极品一区二区| 国产免费视频播放在线视频 | 韩国av在线不卡| 三级国产精品片| 亚洲成人久久爱视频| 乱系列少妇在线播放| 午夜爱爱视频在线播放| 国产亚洲精品久久久com| 午夜精品在线福利| 国产亚洲精品久久久com| 好男人在线观看高清免费视频| 国产成人aa在线观看| 午夜爱爱视频在线播放| av女优亚洲男人天堂| 三级毛片av免费| 丝瓜视频免费看黄片| 麻豆av噜噜一区二区三区| 久久精品久久久久久久性| 欧美成人精品欧美一级黄| 国产亚洲一区二区精品| 久久99热6这里只有精品| 国产v大片淫在线免费观看| 亚洲精品国产av蜜桃| 观看免费一级毛片| 久久精品熟女亚洲av麻豆精品 | 大陆偷拍与自拍| 日本-黄色视频高清免费观看| 亚洲av中文av极速乱| 精品人妻偷拍中文字幕| 日韩欧美 国产精品| 男女啪啪激烈高潮av片| 日韩人妻高清精品专区| 亚洲精品,欧美精品| av免费观看日本| 国产精品av视频在线免费观看| eeuss影院久久| 国产伦精品一区二区三区视频9| 热99在线观看视频| 国产男人的电影天堂91| 一边亲一边摸免费视频| 一级毛片电影观看| 夫妻午夜视频| 国产视频首页在线观看| 亚洲精品国产成人久久av| 深爱激情五月婷婷| 一区二区三区免费毛片| 男插女下体视频免费在线播放| 一本久久精品| 国产精品久久视频播放| 六月丁香七月| 国产亚洲av嫩草精品影院| 国产黄色视频一区二区在线观看| a级毛色黄片| 成人鲁丝片一二三区免费| 亚洲人与动物交配视频| 国产乱来视频区| 特大巨黑吊av在线直播| 九九爱精品视频在线观看| 丝袜喷水一区| 国产精品不卡视频一区二区| 日本熟妇午夜| 成年女人看的毛片在线观看| 亚洲精品aⅴ在线观看| 蜜臀久久99精品久久宅男| 美女主播在线视频| 免费播放大片免费观看视频在线观看| 日韩制服骚丝袜av| 欧美精品国产亚洲| 丰满乱子伦码专区| 少妇猛男粗大的猛烈进出视频 | 乱人视频在线观看| 啦啦啦中文免费视频观看日本| 国产亚洲av嫩草精品影院| 99热6这里只有精品| 免费看av在线观看网站| av线在线观看网站| 天堂网av新在线| 高清视频免费观看一区二区 | 青青草视频在线视频观看| 嫩草影院入口| 春色校园在线视频观看| 久久久久精品性色| 又爽又黄无遮挡网站| 在线观看人妻少妇| 午夜免费观看性视频| 国产又色又爽无遮挡免| 人妻制服诱惑在线中文字幕| 老司机影院毛片| 精品人妻偷拍中文字幕| 亚洲av成人精品一二三区| 日本猛色少妇xxxxx猛交久久| 日本av手机在线免费观看| 久久精品夜色国产| 精品久久久久久成人av| 又黄又爽又刺激的免费视频.| 亚洲欧美成人精品一区二区| av播播在线观看一区| 国产精品三级大全| 一个人观看的视频www高清免费观看| 午夜精品在线福利| 狠狠精品人妻久久久久久综合| 在线观看人妻少妇| 亚洲国产日韩欧美精品在线观看| 亚洲国产高清在线一区二区三| 人妻夜夜爽99麻豆av| 亚洲国产精品sss在线观看| 亚洲内射少妇av| 日韩成人av中文字幕在线观看| 久久久色成人| 久久久久久久大尺度免费视频| 中文字幕久久专区| 十八禁网站网址无遮挡 | 伦精品一区二区三区| 夜夜看夜夜爽夜夜摸| 国产午夜精品久久久久久一区二区三区| 久久这里只有精品中国| 丰满人妻一区二区三区视频av| 免费观看精品视频网站| 高清在线视频一区二区三区| 免费av观看视频| 成年女人看的毛片在线观看| 麻豆久久精品国产亚洲av| 97人妻精品一区二区三区麻豆| 丰满人妻一区二区三区视频av| 亚洲18禁久久av| 日本爱情动作片www.在线观看| 日韩欧美精品免费久久| 建设人人有责人人尽责人人享有的 | 99热网站在线观看| 欧美另类一区| 欧美xxxx黑人xx丫x性爽| 三级毛片av免费| 精品久久久久久久久久久久久| 国产午夜精品一二区理论片| 国产精品1区2区在线观看.| 久久久亚洲精品成人影院| 久久久亚洲精品成人影院| 久久久久久久久大av| 免费观看无遮挡的男女| 日韩 亚洲 欧美在线| 久久久久久久大尺度免费视频| 免费看不卡的av| 一级毛片黄色毛片免费观看视频| 欧美日韩精品成人综合77777| 高清av免费在线| 一二三四中文在线观看免费高清| 男人和女人高潮做爰伦理| 亚洲图色成人| 十八禁网站网址无遮挡 | 国产精品1区2区在线观看.| 亚洲精品自拍成人| 日韩大片免费观看网站| 免费观看性生交大片5| 国产精品美女特级片免费视频播放器| 高清毛片免费看| 不卡视频在线观看欧美| 亚洲精品视频女| xxx大片免费视频| 久久久久九九精品影院| 蜜桃亚洲精品一区二区三区| 秋霞伦理黄片| av.在线天堂| 国产成人91sexporn| 日本免费a在线| 久久久色成人| 国产色爽女视频免费观看| av福利片在线观看| 婷婷六月久久综合丁香| 嫩草影院新地址| 能在线免费观看的黄片| 天堂影院成人在线观看| 天天一区二区日本电影三级| 免费高清在线观看视频在线观看| 精品午夜福利在线看| 国产精品久久视频播放| 久久久久久久久久久丰满| 麻豆成人av视频| 亚洲精品乱码久久久久久按摩| 亚洲自偷自拍三级| 性色avwww在线观看| 别揉我奶头 嗯啊视频| 国产熟女欧美一区二区| 亚洲怡红院男人天堂| 日日干狠狠操夜夜爽| 成年免费大片在线观看| 在线播放无遮挡| 国产极品天堂在线| 青春草亚洲视频在线观看| 狂野欧美白嫩少妇大欣赏| 日本一本二区三区精品| 人体艺术视频欧美日本| 搡女人真爽免费视频火全软件| 久久精品久久精品一区二区三区| 午夜精品国产一区二区电影 | 亚洲国产精品成人久久小说| 亚洲在线观看片| 99视频精品全部免费 在线| 色综合亚洲欧美另类图片| 精品久久久精品久久久| 免费看日本二区| 国产伦在线观看视频一区| 国产黄频视频在线观看| 噜噜噜噜噜久久久久久91| 久久久久性生活片| 中国美白少妇内射xxxbb| 国内精品美女久久久久久| 日韩欧美三级三区| 免费观看的影片在线观看| 亚洲人成网站在线观看播放| 亚洲精品,欧美精品| 久久久国产一区二区| 成人性生交大片免费视频hd| 大又大粗又爽又黄少妇毛片口| 插逼视频在线观看| 免费av毛片视频| 亚洲va在线va天堂va国产| 在线观看人妻少妇| av.在线天堂| 久久久国产一区二区| 狂野欧美白嫩少妇大欣赏| 国产成人午夜福利电影在线观看| 日韩av在线免费看完整版不卡| 精品一区在线观看国产| 午夜激情欧美在线| 97超碰精品成人国产| av播播在线观看一区| 国产高清有码在线观看视频| 日本爱情动作片www.在线观看| 春色校园在线视频观看| 一级av片app| 尤物成人国产欧美一区二区三区| .国产精品久久| 久久这里只有精品中国| 青青草视频在线视频观看| 秋霞在线观看毛片| 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 一边亲一边摸免费视频| 久久人人爽人人片av| 免费av观看视频| 黄片无遮挡物在线观看| 嫩草影院新地址| 中文字幕亚洲精品专区| 人妻制服诱惑在线中文字幕| 亚洲国产精品sss在线观看| 嫩草影院新地址| 国产乱来视频区| 一个人看视频在线观看www免费| 亚洲欧美日韩无卡精品| 久久精品久久精品一区二区三区| 亚洲精品中文字幕在线视频 | 晚上一个人看的免费电影| 不卡视频在线观看欧美| 26uuu在线亚洲综合色| 久久韩国三级中文字幕| 亚洲精品一区蜜桃| 白带黄色成豆腐渣| 日韩精品有码人妻一区| 成人午夜高清在线视频| 午夜精品国产一区二区电影 | 国产av国产精品国产| 人妻系列 视频| 不卡视频在线观看欧美| 久久精品夜夜夜夜夜久久蜜豆| 色播亚洲综合网| 国产在视频线在精品| av.在线天堂| 成人特级av手机在线观看| 国产精品一二三区在线看| 人人妻人人看人人澡| 亚洲欧美精品专区久久| 日本免费a在线| 狠狠精品人妻久久久久久综合| 久久久久久久久久人人人人人人| 国产精品一及| 成年人午夜在线观看视频 | 欧美区成人在线视频| 久久久国产一区二区| 国产视频内射| 亚洲色图av天堂| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 亚洲欧美日韩东京热| 真实男女啪啪啪动态图| 国产免费又黄又爽又色| 国产三级在线视频| 亚洲婷婷狠狠爱综合网| 波多野结衣巨乳人妻| 亚洲成人中文字幕在线播放| 男插女下体视频免费在线播放| 亚洲国产日韩欧美精品在线观看| 亚洲精品乱久久久久久| 国产精品久久久久久久久免| 日韩av不卡免费在线播放| 婷婷色综合www| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 婷婷色av中文字幕| 亚洲av日韩在线播放| 日日摸夜夜添夜夜爱| 一级毛片aaaaaa免费看小| 只有这里有精品99| 97热精品久久久久久| 日韩欧美一区视频在线观看 | 精品熟女少妇av免费看| 国产国拍精品亚洲av在线观看| 色综合站精品国产| 国产91av在线免费观看| 日韩精品青青久久久久久| 色播亚洲综合网| 午夜精品一区二区三区免费看| 婷婷色麻豆天堂久久| 午夜福利在线观看吧| 尾随美女入室| 99久久中文字幕三级久久日本| 成年免费大片在线观看| 人妻夜夜爽99麻豆av| 69av精品久久久久久| 国产在线男女| 亚洲av成人精品一二三区| 一级av片app| 日韩欧美 国产精品| 亚洲欧美清纯卡通| 精品少妇黑人巨大在线播放| 床上黄色一级片| 久久精品国产亚洲av天美| 在线观看av片永久免费下载| 18禁在线播放成人免费| 免费少妇av软件| 日韩不卡一区二区三区视频在线| 亚洲精品第二区| 人人妻人人澡人人爽人人夜夜 | 日韩av在线大香蕉| 夜夜爽夜夜爽视频| 少妇被粗大猛烈的视频| 久久人人爽人人片av| 全区人妻精品视频| 欧美日韩视频高清一区二区三区二| 欧美一级a爱片免费观看看| 欧美日韩精品成人综合77777| 在线 av 中文字幕| 少妇熟女aⅴ在线视频| 床上黄色一级片| 亚洲av在线观看美女高潮| 91av网一区二区| 老司机影院成人| 色综合色国产| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 国产精品av视频在线免费观看| 乱系列少妇在线播放| 成人亚洲欧美一区二区av| 色哟哟·www| 中文字幕免费在线视频6| 亚洲欧美一区二区三区黑人 | 亚洲色图av天堂| 国产视频首页在线观看| 亚洲三级黄色毛片| 色综合亚洲欧美另类图片| 久久久久久久大尺度免费视频| 国产亚洲5aaaaa淫片| 欧美xxxx性猛交bbbb| 黄色一级大片看看| 青春草亚洲视频在线观看| 成人综合一区亚洲| 少妇人妻一区二区三区视频| 最近2019中文字幕mv第一页| 欧美日韩综合久久久久久| 啦啦啦中文免费视频观看日本| 波多野结衣巨乳人妻| 七月丁香在线播放| 亚洲久久久久久中文字幕| 国产午夜精品久久久久久一区二区三区| 国产精品嫩草影院av在线观看| 免费无遮挡裸体视频| 最近中文字幕高清免费大全6| 中文字幕av在线有码专区| 99久久中文字幕三级久久日本| 韩国高清视频一区二区三区| 国产免费一级a男人的天堂| videossex国产| 99热这里只有精品一区| 男女下面进入的视频免费午夜| 国产精品蜜桃在线观看| 男人舔女人下体高潮全视频| 99久久人妻综合| 日本wwww免费看| 99久久人妻综合| 99久久精品热视频| 亚洲在线自拍视频| 亚洲av成人av| 18+在线观看网站| xxx大片免费视频| 免费黄频网站在线观看国产| 男女啪啪激烈高潮av片| 国产视频内射| 欧美日韩在线观看h| 国产精品一二三区在线看| 亚洲av国产av综合av卡| 人妻系列 视频| 永久免费av网站大全| 亚洲国产色片| 午夜福利网站1000一区二区三区| 婷婷色综合www| 国产av码专区亚洲av| 午夜免费激情av| 寂寞人妻少妇视频99o| 在现免费观看毛片| ponron亚洲| 久久精品夜色国产| 亚洲精品久久久久久婷婷小说| av女优亚洲男人天堂| 国产有黄有色有爽视频| 一个人观看的视频www高清免费观看| 午夜免费男女啪啪视频观看| www.色视频.com| 国产v大片淫在线免费观看| 亚洲欧美精品专区久久| 大香蕉97超碰在线| 国产精品1区2区在线观看.| 超碰av人人做人人爽久久| 秋霞伦理黄片| 免费av毛片视频| 精品久久久久久成人av| 91aial.com中文字幕在线观看| 在线免费十八禁| 高清午夜精品一区二区三区| 精品久久国产蜜桃| 国产色爽女视频免费观看| 97超视频在线观看视频| 十八禁国产超污无遮挡网站| 国产 亚洲一区二区三区 | 99久久精品国产国产毛片| 伦精品一区二区三区| 久久久国产一区二区| 狠狠精品人妻久久久久久综合| 国产午夜精品久久久久久一区二区三区| 51国产日韩欧美| 久久精品人妻少妇| 国产成人aa在线观看| 寂寞人妻少妇视频99o| 18禁裸乳无遮挡免费网站照片| 黄色日韩在线| 日日干狠狠操夜夜爽| 熟妇人妻久久中文字幕3abv| 欧美一区二区亚洲| 国产成人一区二区在线| 日本猛色少妇xxxxx猛交久久| 国产淫语在线视频| xxx大片免费视频|