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

    炭質(zhì)泥頁巖剪切破壞聲發(fā)射特性及其分形特征

    2024-05-23 14:07:50代鑫胡斌李京曹建軍崔喜興
    科學(xué)技術(shù)與工程 2024年12期
    關(guān)鍵詞:振鈴炭質(zhì)維數(shù)

    代鑫, 胡斌*, 李京, 曹建軍, 崔喜興

    (1.武漢科技大學(xué)資源與環(huán)境工程學(xué)院, 武漢 430081; 2.冶金礦產(chǎn)資源高效利用與造塊湖北省重點(diǎn)實(shí)驗(yàn)室, 武漢 430081)

    含炭質(zhì)泥頁巖軟弱夾層廣泛分布于中國(guó)中南、西南地區(qū)的礦山邊坡中,易受自然環(huán)境變化或人為擾動(dòng)影響而失穩(wěn),對(duì)邊坡工程和隧道圍巖等工程的穩(wěn)定性存在嚴(yán)重威脅,但巖石的失穩(wěn)破壞過程可通過聲發(fā)射技術(shù)進(jìn)行監(jiān)測(cè)和預(yù)報(bào),聲發(fā)射參數(shù)中蘊(yùn)含著巖石破壞過程信息,同時(shí)聲發(fā)射分形理論可以通過分形維數(shù)的變化來反映巖石破壞過程中微裂紋的演化規(guī)律和破壞前兆特征[1]。因此,可通過炭質(zhì)泥頁巖破壞的聲發(fā)射參數(shù)及其分形特征來探究含炭質(zhì)泥頁巖軟弱夾層的破壞機(jī)理和失穩(wěn)前兆特征。

    近些年來,許多學(xué)者將聲發(fā)射技術(shù)和分形理論應(yīng)用到巖石破裂機(jī)制的研究中,王創(chuàng)業(yè)等[2]探究了巖石和類巖石材料在單軸壓縮破壞過程中的聲發(fā)射特性及其差異。胡英等[3]運(yùn)用時(shí)頻分析方法,提出一種新的特征參數(shù),基于此得到巖石破裂前兆信息。高安森等[4]通過紅砂巖單軸壓縮聲發(fā)射試驗(yàn),探究了紅砂巖在臨界加速失穩(wěn)階段的破壞特征和失穩(wěn)前兆。劉文德等[5]通過對(duì)灰?guī)r開展單軸壓縮試驗(yàn),證明了灰?guī)r破裂所產(chǎn)生的聲發(fā)射事件數(shù)的時(shí)間序列存在分形特征。龔囪等[6]對(duì)紅砂巖進(jìn)行了蠕變聲發(fā)射實(shí)驗(yàn),運(yùn)用分形理論,得出聲發(fā)射震源值演化規(guī)律。孫博等[7]進(jìn)行不同層理傾角的板巖單軸壓縮試驗(yàn),分析了聲發(fā)射計(jì)數(shù)的多重分形特征,得出板巖的破壞前兆和預(yù)警時(shí)間。龔囪等[8]研究了鋼筋混凝土受剪梁聲發(fā)射參數(shù)在不同剪跨比和粗骨粒徑的條件下的發(fā)展規(guī)律。王偉等[9]通過對(duì)砂巖進(jìn)行常規(guī)三軸和三軸循環(huán)加載聲發(fā)射試驗(yàn),分析了砂巖聲發(fā)射特性和聲發(fā)射計(jì)數(shù)分形特征,得到了巖石裂痕發(fā)育的演化規(guī)律。楊東杰等[10]設(shè)計(jì)了巖石預(yù)損試驗(yàn),闡述了預(yù)損后花崗巖的振鈴計(jì)數(shù)分形演化。以上研究表明巖石、類巖石等材料的聲發(fā)射參數(shù)都具有分形特征且蘊(yùn)含著失穩(wěn)破壞過程信息,但大多數(shù)學(xué)者只采用一種聲發(fā)射參數(shù)進(jìn)行分形特征的研究,所得出的失穩(wěn)前兆信息單一且難以得出不同參數(shù)中所蘊(yùn)含失穩(wěn)前兆信息之間的差別,同時(shí)軟弱夾層(炭質(zhì)泥頁巖)作為礦山邊坡中極易造成滑坡事故的危險(xiǎn)源,卻少有學(xué)者對(duì)其進(jìn)行基于聲發(fā)射的失穩(wěn)前兆研究。

    考慮到上述現(xiàn)有研究的不足,對(duì)少有學(xué)者關(guān)注的軟弱夾層(炭質(zhì)泥頁巖)進(jìn)行剪切破壞聲發(fā)射試驗(yàn),考慮了不同聲發(fā)射參數(shù)之間的聯(lián)系與差別,選取多個(gè)參數(shù)進(jìn)行分形特征研究,綜合確定炭質(zhì)泥頁巖失穩(wěn)前兆信息,以避免研究結(jié)果的依據(jù)單一,同時(shí),對(duì)比分析,探索更適合作為炭質(zhì)泥頁巖破壞前兆信息的參數(shù),旨在為巖石失穩(wěn)破壞全程裂隙演化和失穩(wěn)預(yù)測(cè)提供理論依據(jù)。

    1 聲發(fā)射試驗(yàn)

    1.1 巖樣制備

    試驗(yàn)所用炭質(zhì)泥頁巖取自四川省峨眉山地區(qū)二疊系堅(jiān)硬石灰?guī)r層中的軟弱夾層,夾層厚約為0.1~0.8 m,產(chǎn)狀與上下巖層一致。將所取巖樣切割為尺寸為150 mm(長(zhǎng))×75 mm(寬)×75 mm(高)的長(zhǎng)方體試件(圖1),加工過程中盡量保證試樣不受擾動(dòng)。

    圖1 炭質(zhì)泥頁巖巖樣Fig.1 Carbonaceous shale

    1.2 試驗(yàn)設(shè)備及方法

    巖石剪切試驗(yàn)在團(tuán)隊(duì)自主研發(fā)的微機(jī)伺服控制巖石剪切-滲流耦合試驗(yàn)機(jī)上進(jìn)行,設(shè)備如圖2所示。主要適用于巖石剪切試驗(yàn)、剪切滲流和剪切蠕變?cè)囼?yàn)。該設(shè)備采用高剛度整體框架主機(jī)、全數(shù)字測(cè)控器,具有多個(gè)測(cè)控通道可對(duì)任意一個(gè)通道進(jìn)行閉環(huán)控制。

    圖2 WAJS-100型微機(jī)伺服控制巖石剪切-滲流耦合試驗(yàn)機(jī)Fig.2 WAJS-100 microcomputer servo control shear-seepage coupling testing machine for rock

    此次試驗(yàn)采用美國(guó)物理聲學(xué)公司生產(chǎn)的PCI-2型聲發(fā)射測(cè)試分析系統(tǒng)(圖3),選擇的傳感器為RS-2A型,前置放大器可以將信號(hào)放大100倍。試驗(yàn)時(shí)將4個(gè)探頭兩兩一組,貼于巖樣兩個(gè)相對(duì)的側(cè)面,并使4個(gè)探頭不在同一平面,采樣速率設(shè)為3 MHz,門檻電壓值為100 mV。

    圖3 聲發(fā)射測(cè)試分析系統(tǒng)、傳感器、前置放大器Fig.3 Acoustic emission system of analysis, sensors and preamplifiers

    對(duì)炭質(zhì)泥頁巖進(jìn)行剪切破壞試驗(yàn)。法向加載采用力控制,切向加載選擇位移加載,加載速率設(shè)置0.2 mm/min。先施加法向力,穩(wěn)定30 s后,再施加剪切力。試驗(yàn)過程中采集時(shí)間、剪切位移、聲發(fā)射振鈴計(jì)數(shù)、能量、振幅等試驗(yàn)數(shù)據(jù)。

    2 聲發(fā)射試驗(yàn)結(jié)果分析

    一般地,當(dāng)巖石在剪切荷載下發(fā)生破壞,引起巖石內(nèi)部部分能量以應(yīng)力波的形式擴(kuò)散,即形成聲發(fā)射信號(hào)[11]。選擇振鈴計(jì)數(shù)、能量、振幅、上升時(shí)間與幅度的比值RA、振鈴計(jì)數(shù)與持續(xù)時(shí)間的比值A(chǔ)F對(duì)試驗(yàn)結(jié)果進(jìn)行分析。圖4~圖6為剪切力、時(shí)間和不同聲發(fā)射參數(shù)之間的關(guān)系曲線。

    圖4 剪切力-時(shí)間-振鈴計(jì)數(shù)曲線Fig.4 Stress-time-counts curve

    圖5 剪切力-時(shí)間-能量曲線Fig.5 Stress-time-energy curve

    圖6 剪切力-時(shí)間-振幅曲線Fig.6 Stress-time-amplitude curve

    從圖4~圖6可以看出,在實(shí)驗(yàn)剛開始時(shí)巖石內(nèi)就發(fā)生了聲發(fā)射活動(dòng),這是因?yàn)榇藭r(shí)巖石處于壓密階段(0~150 s),隨著原有裂隙被壓實(shí),巖石顆粒間發(fā)生錯(cuò)動(dòng),同時(shí)也伴有少量微裂紋的生成,聲發(fā)射特性表現(xiàn)為小事件較多,能量值和振幅值波動(dòng)較大,新產(chǎn)生的微破裂強(qiáng)度大小不一。隨試驗(yàn)的進(jìn)行,聲發(fā)射信號(hào)值都逐漸下降,聲發(fā)射能量和振幅逐漸稀疏,這是因?yàn)閹r石進(jìn)入彈性變形階段(150 ~400 s),原有裂隙已被壓實(shí),巖石發(fā)生彈性變形,體積減小,聲發(fā)射特性表現(xiàn)為能量逐漸降低,振幅相對(duì)減小。在巖石進(jìn)入塑性變形階段(400~500 s),新裂紋不斷產(chǎn)生、匯集、貫通,聲發(fā)射特性表現(xiàn)為聲發(fā)射振鈴計(jì)數(shù)和能量值在總體趨勢(shì)上逐漸上升,在臨近破壞時(shí)達(dá)到峰值。

    圖7為炭質(zhì)泥頁巖臨近破壞時(shí)裂紋的聲發(fā)射三維定位演化圖,可以看出,紅色球體及其在x、y、z面上投影數(shù)量上越來越多,形狀也越來越大,說明在此時(shí)巖石內(nèi)部已產(chǎn)生大量裂紋并向四周擴(kuò)展,印證了此時(shí)聲發(fā)射振鈴計(jì)數(shù)和能量值極具驟增的變化。值得注意的是,此階段巖石出現(xiàn)約30 s的聲發(fā)射平靜期(450~480 s),平靜期內(nèi)振鈴計(jì)數(shù)、能量值均有明顯的大幅減弱,平靜期結(jié)束后聲發(fā)射活動(dòng)又瞬間激增,緊接著巖石發(fā)生破壞,平靜期這一特殊現(xiàn)象可作為巖石破壞前兆信息。

    圖7 炭質(zhì)泥頁巖臨近破壞時(shí)裂紋的聲發(fā)射三維定位演化圖Fig.7 Three-dimensional localization evolution of cracks in carbonaceous shale near destruction

    聲發(fā)射RA、AF作為巖石破裂機(jī)制的依據(jù),常用來區(qū)分巖石破壞類型,判定方法如圖8所示,高RA低AF為聲發(fā)射剪切波特征,對(duì)應(yīng)巖石剪破壞,高AF值低RA則對(duì)應(yīng)張拉破壞[12]。RA為上升時(shí)間與幅度的比值,AF為振鈴計(jì)數(shù)與持續(xù)時(shí)間的比值,巖石RA-AF分布如圖9所示,其中高AF低RA的分布明顯少于高RA低AF,故而炭質(zhì)泥頁巖樣的破壞類型是為以剪切破壞為主的張剪復(fù)合型破壞。

    圖8 破壞類型判定方法Fig.8 Criteria for fracture classification

    圖9 RA-AF分布圖Fig.9 RA-AF distribution

    3 聲發(fā)射信號(hào)分形特征研究

    3.1 分形理論和關(guān)聯(lián)維數(shù)計(jì)算模型

    一般地,常用的聲發(fā)射表征參數(shù)都可以構(gòu)成單獨(dú)的聲發(fā)射時(shí)間序列[13],分形維數(shù)是分形理論的基本量,它可以定量描述分形[14],其中關(guān)聯(lián)維數(shù)是最常用的分形維數(shù)之一。G-P算法可將聲發(fā)射基本參數(shù)時(shí)間序列作為研究對(duì)象,每個(gè)聲發(fā)射基本參數(shù)時(shí)間序列對(duì)應(yīng)容量為n的序列集Xn={x1,x2,…,xn},然后在該時(shí)間序列的基礎(chǔ)上構(gòu)造m維的相空間(m

    (1)

    對(duì)于給定尺度r,可表示為

    r=kr0

    (2)

    式(2)中:r0可表示為

    (3)

    對(duì)每一個(gè)給定的尺度,都可得到一個(gè)C(r)。在雙對(duì)數(shù)坐標(biāo)系中可得到一系列點(diǎn):[lnCr(k)],lnr(k)],對(duì)這些點(diǎn)進(jìn)行數(shù)據(jù)擬合,擬合直線的斜率就是聲發(fā)射參數(shù)序列的關(guān)聯(lián)維數(shù)D,即

    (4)

    3.2 聲發(fā)射分形特征的判定和相空間維數(shù)的確定

    選擇振鈴計(jì)數(shù)、能量、振幅3個(gè)參數(shù)進(jìn)行關(guān)聯(lián)維數(shù)的計(jì)算和分析。首先編寫該算法的程序,然后將聲發(fā)射參數(shù)、時(shí)間序列、應(yīng)變等導(dǎo)入程序,實(shí)現(xiàn)對(duì)聲發(fā)射參數(shù)的歸一化處理、相空間重構(gòu)和關(guān)聯(lián)維數(shù)的計(jì)算。

    G-P算法重構(gòu)了巖石損傷破壞過程動(dòng)力系統(tǒng)吸引子的多維空問,但對(duì)于系統(tǒng)空間的吸引子的性質(zhì),相空間維數(shù)m的取值至關(guān)重要[16]。在對(duì)m取值時(shí),不可將m盲目縮小,這樣會(huì)導(dǎo)致重構(gòu)的相空間處于低維狀態(tài),使得兩個(gè)原本不接近的相點(diǎn)產(chǎn)生相鄰的假象。也不可將m過度放大,這樣會(huì)使得原本接近的相點(diǎn)造成相離的假象,導(dǎo)致相空間吸引子分散。所以,應(yīng)選取合適的m進(jìn)行計(jì)算。

    m取值一般采用幾何不變量法[17],在MATLAB環(huán)境下將m依次取2、4、6、8、10、12、14、16、18、20,在600個(gè)觀測(cè)尺度r下分別分析3個(gè)聲發(fā)射參數(shù)的r-C(r)雙對(duì)數(shù)關(guān)系。

    如圖10所示,對(duì)于振鈴計(jì)數(shù)、能量、振幅三個(gè)聲發(fā)射參數(shù)在同一m下,當(dāng)r很小時(shí),每?jī)蓚€(gè)散點(diǎn)的距離都比較大,即每?jī)蓚€(gè)點(diǎn)的C(r)相差很大,相點(diǎn)距離都比較遠(yuǎn),所以對(duì)于r的小幅增大,C(r)并沒有呈現(xiàn)連續(xù)的增大,而是跳躍式增加,即不存在分形現(xiàn)象。當(dāng)r取值很大時(shí),C(r)趨近于1,lnC(r)趨近于0,散點(diǎn)構(gòu)成曲線雖然連續(xù)但是已成為一條直線,隨r的增大,關(guān)聯(lián)維數(shù)一直為0,不存在分形現(xiàn)象。而當(dāng)lnr取值為2~4時(shí),振鈴計(jì)數(shù)、能量、振幅3個(gè)聲發(fā)射參數(shù)在10個(gè)m下的觀測(cè)尺度r和關(guān)聯(lián)函數(shù)C(r)的對(duì)數(shù)曲線會(huì)部分接近一次線性關(guān)系,即此時(shí)聲發(fā)射參數(shù)是具有分形特征的,而r-C(r)雙對(duì)數(shù)散點(diǎn)曲線線性回歸的斜率就是關(guān)聯(lián)維數(shù)。

    圖10 不同聲發(fā)射參數(shù)在不同m值下r-C(r)雙對(duì)數(shù)散點(diǎn)圖Fig.10 r-C(r) double logarithmic scatter plot for different m values and parameters

    如圖11所示,3個(gè)參數(shù)的關(guān)聯(lián)維數(shù)D隨m的增大先呈線性增加,在m增大到某一點(diǎn)后,曲線斜率逐漸下降且D增長(zhǎng)速度逐漸穩(wěn)定。一般取D收斂的第一點(diǎn)作為合適的相空間維數(shù)m,即在計(jì)算3種聲發(fā)射參數(shù)的關(guān)聯(lián)維數(shù)時(shí)m分別取10、12、14。

    圖11 不同聲發(fā)射參數(shù)在不同m值下的關(guān)聯(lián)維數(shù)DFig.11 The correlation dimension D of different acoustic emission parameters at different values of m

    3.3 聲發(fā)射分形演化研究

    通過聲發(fā)射參數(shù)計(jì)算出來的關(guān)聯(lián)維數(shù)可以描述巖石聲發(fā)射分形演化行為[18]?;?.2節(jié)對(duì)分形特征的判定和確定的m,在MATLAB環(huán)境下,每隔2%的應(yīng)變比對(duì)聲發(fā)射參數(shù)進(jìn)行一次關(guān)聯(lián)維數(shù)D的計(jì)算,直至應(yīng)變比達(dá)到100%。

    圖12為在不同應(yīng)變比下的分形維數(shù)和對(duì)應(yīng)計(jì)算位置剪應(yīng)力的曲線,可以看出,在加載初期,振鈴計(jì)數(shù)、能量、振幅的關(guān)聯(lián)維數(shù)不斷升高,進(jìn)入升維狀態(tài),此時(shí)巖石內(nèi)部的原生裂隙被壓實(shí)閉合,而后又隨壓力的增加,這種被壓密的有序性逐步提升,故而關(guān)聯(lián)維數(shù)逐漸上升。隨試驗(yàn)進(jìn)行,能量的關(guān)聯(lián)維數(shù)率先發(fā)生變化,由上升趨勢(shì)轉(zhuǎn)為上下波動(dòng),表明能量的關(guān)聯(lián)維數(shù)進(jìn)入了動(dòng)態(tài)波動(dòng)階段,而后振鈴計(jì)數(shù)與振幅的關(guān)聯(lián)維數(shù)才相繼進(jìn)入動(dòng)態(tài)波動(dòng)階段,在此階段巖石內(nèi)部不斷出現(xiàn)微破裂的發(fā)育和閉合,但微破裂的發(fā)育和閉合隨時(shí)間不均勻出現(xiàn),處于一種“混沌”狀態(tài),所以關(guān)聯(lián)維數(shù)呈現(xiàn)出無規(guī)則的上下波動(dòng)。隨壓力進(jìn)一步升高,能量的關(guān)聯(lián)維數(shù)最先進(jìn)入降維階段且降維幅度較大,持續(xù)時(shí)間較長(zhǎng),而振鈴計(jì)數(shù)和振幅的關(guān)聯(lián)維數(shù)分別在塑性變形階段中期和后期才進(jìn)入降維階段且降維幅度較小,持續(xù)時(shí)間較短。關(guān)聯(lián)維數(shù)的降維現(xiàn)象意味著此時(shí)聲發(fā)射信號(hào)變得活躍,巖石內(nèi)部開始有大量的裂紋產(chǎn)生、擴(kuò)展和貫通,且降維現(xiàn)象出現(xiàn)在巖石破壞前,所以關(guān)聯(lián)維數(shù)的降維現(xiàn)象可作為巖石破壞的前兆特征。

    圖12 不同參數(shù)的關(guān)聯(lián)維數(shù)-應(yīng)變-剪應(yīng)力曲線Fig.12 Correlation dimension-strain-shear stress curve of different parameters

    關(guān)聯(lián)維數(shù)在炭質(zhì)泥頁巖剪切破壞全過程中呈現(xiàn)出“升維-動(dòng)態(tài)波動(dòng)-降維”的演化規(guī)律。能量的關(guān)聯(lián)維數(shù)比其他參數(shù)的關(guān)聯(lián)維數(shù)更先進(jìn)入動(dòng)態(tài)波動(dòng)和降維階段,說明關(guān)聯(lián)維數(shù)的演化對(duì)聲發(fā)射能量的變化更敏感,且能量的關(guān)聯(lián)維數(shù)的降維程度比其他參數(shù)更顯著,持續(xù)時(shí)間更長(zhǎng),所以能量關(guān)聯(lián)維數(shù)的降維現(xiàn)象更適合用來預(yù)警炭質(zhì)泥頁巖失穩(wěn)破壞的前兆信息。

    4 結(jié)論

    通過對(duì)炭質(zhì)泥頁巖剪切試驗(yàn),對(duì)采集到的聲發(fā)射數(shù)據(jù)及其分形特征進(jìn)行研究,得出以下結(jié)論。

    (1)聲發(fā)射活動(dòng)在炭質(zhì)泥頁巖的不同變形階段活躍程度不同,在塑性變形階段最強(qiáng),其次是壓密階段,最后是彈性變形階段,且聲發(fā)射活動(dòng)在巖石破壞前存在平靜期。裂紋在巖樣中心產(chǎn)生向四周不斷擴(kuò)張,破壞模式為以剪破壞為主的拉剪復(fù)合型破壞。

    (2)炭質(zhì)泥頁巖的聲發(fā)射振鈴計(jì)數(shù)、能量、振幅的時(shí)間序列都具有分形特征,振鈴計(jì)數(shù)、能量、振幅的關(guān)聯(lián)維數(shù)在破壞全程呈現(xiàn)出“升維-動(dòng)態(tài)波動(dòng)-降維”的演化規(guī)律,

    (3)聲發(fā)射平靜期可作為炭質(zhì)泥頁巖破壞前兆信息。相比于振幅和振鈴計(jì)數(shù),關(guān)聯(lián)維數(shù)的演化對(duì)聲發(fā)射能量的變化更敏感,且能量的降維現(xiàn)象更顯著,持續(xù)時(shí)間更長(zhǎng),更適合作為炭質(zhì)泥頁巖失穩(wěn)破壞的前兆特征。

    (4)若將工程監(jiān)測(cè)手段與平靜期、聲發(fā)射能量的分形維數(shù)演化特征相結(jié)合,對(duì)邊坡巖體失穩(wěn)預(yù)警具有一定應(yīng)用價(jià)值。

    猜你喜歡
    振鈴炭質(zhì)維數(shù)
    β-變換中一致丟番圖逼近問題的維數(shù)理論
    基于敏感性分析的炭質(zhì)板巖引水隧洞支護(hù)結(jié)構(gòu)優(yōu)化研究
    高速公路炭質(zhì)頁巖高邊坡穩(wěn)定性分析及施工方案
    振鈴現(xiàn)象對(duì)心電圖分析和診斷的影響
    基于連續(xù)-非連續(xù)單元方法的炭質(zhì)板巖隧道圍巖穩(wěn)定分析
    一類齊次Moran集的上盒維數(shù)
    家庭網(wǎng)關(guān)振鈴業(yè)務(wù)配置策略研究
    干燥和飽水炭質(zhì)板巖流變力學(xué)特性與模型研究
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    涉及相變問題Julia集的Hausdorff維數(shù)
    久久人人爽人人爽人人片va| 免费av观看视频| 色综合婷婷激情| 欧美激情久久久久久爽电影| 波多野结衣高清无吗| 热99在线观看视频| 日韩中字成人| 亚洲国产精品sss在线观看| 久久精品人妻少妇| 我的老师免费观看完整版| 成年女人看的毛片在线观看| 成年版毛片免费区| h日本视频在线播放| 窝窝影院91人妻| 午夜久久久久精精品| 尤物成人国产欧美一区二区三区| 国产精品久久久久久久久免| 亚洲精品456在线播放app | 99久久九九国产精品国产免费| 免费无遮挡裸体视频| 亚洲电影在线观看av| 免费人成在线观看视频色| 尾随美女入室| 日韩一本色道免费dvd| 久久精品夜夜夜夜夜久久蜜豆| 极品教师在线免费播放| 中文字幕精品亚洲无线码一区| 人妻久久中文字幕网| 国产精品电影一区二区三区| 麻豆久久精品国产亚洲av| 免费av毛片视频| 观看美女的网站| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 大又大粗又爽又黄少妇毛片口| 伦理电影大哥的女人| 99久久久亚洲精品蜜臀av| 香蕉av资源在线| 日本免费a在线| 人妻少妇偷人精品九色| 亚洲图色成人| 日韩国内少妇激情av| 国产精品亚洲一级av第二区| 麻豆成人av在线观看| 亚洲avbb在线观看| av女优亚洲男人天堂| 99精品久久久久人妻精品| 久久久色成人| 国产亚洲av嫩草精品影院| 免费搜索国产男女视频| 最近中文字幕高清免费大全6 | 国产亚洲91精品色在线| 九色成人免费人妻av| 狂野欧美激情性xxxx在线观看| 久久6这里有精品| 午夜免费激情av| 国产成人影院久久av| 97热精品久久久久久| 午夜亚洲福利在线播放| 校园人妻丝袜中文字幕| 亚洲av一区综合| 国产一区二区亚洲精品在线观看| 91久久精品国产一区二区三区| 99热网站在线观看| 国产成年人精品一区二区| 黄色视频,在线免费观看| 嫩草影院新地址| www.www免费av| 桃红色精品国产亚洲av| 国产成人影院久久av| 国产亚洲av嫩草精品影院| 亚洲不卡免费看| 欧美+日韩+精品| 欧美区成人在线视频| 日韩精品有码人妻一区| 欧美成人性av电影在线观看| 女的被弄到高潮叫床怎么办 | 91麻豆精品激情在线观看国产| 一级毛片久久久久久久久女| 亚洲avbb在线观看| 亚洲色图av天堂| 成人三级黄色视频| 天天躁日日操中文字幕| 精品久久久久久成人av| 国产精品98久久久久久宅男小说| 精品国产三级普通话版| 男女视频在线观看网站免费| 精品久久久久久成人av| 99久久九九国产精品国产免费| 国产男靠女视频免费网站| 99热只有精品国产| 男女做爰动态图高潮gif福利片| 尤物成人国产欧美一区二区三区| 亚洲av不卡在线观看| 日本在线视频免费播放| 国产精品一区二区免费欧美| 亚洲在线观看片| 看片在线看免费视频| 欧美绝顶高潮抽搐喷水| av天堂中文字幕网| 精品一区二区三区人妻视频| 午夜免费成人在线视频| 国产高清激情床上av| 床上黄色一级片| 久久国内精品自在自线图片| 亚洲欧美日韩高清在线视频| 午夜福利在线在线| 中文字幕熟女人妻在线| 美女 人体艺术 gogo| 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 麻豆国产97在线/欧美| 亚洲国产色片| 在线国产一区二区在线| 午夜a级毛片| 美女xxoo啪啪120秒动态图| 精品无人区乱码1区二区| 欧洲精品卡2卡3卡4卡5卡区| 99精品在免费线老司机午夜| 日韩国内少妇激情av| 国产亚洲精品综合一区在线观看| 欧美精品啪啪一区二区三区| 亚洲电影在线观看av| 嫩草影视91久久| 一级a爱片免费观看的视频| 日本黄大片高清| 麻豆久久精品国产亚洲av| 日韩精品青青久久久久久| 少妇人妻一区二区三区视频| 久久99热这里只有精品18| 99久久中文字幕三级久久日本| 欧美一区二区亚洲| 午夜老司机福利剧场| 亚洲成人久久爱视频| 别揉我奶头 嗯啊视频| 麻豆一二三区av精品| 搞女人的毛片| 动漫黄色视频在线观看| 91麻豆精品激情在线观看国产| 成人午夜高清在线视频| 国产精品一区二区免费欧美| 神马国产精品三级电影在线观看| 久久精品影院6| 非洲黑人性xxxx精品又粗又长| 亚洲精华国产精华液的使用体验 | 成人一区二区视频在线观看| 精品久久久久久久久亚洲 | 成人国产一区最新在线观看| 91久久精品国产一区二区成人| 日韩欧美国产在线观看| 国产三级在线视频| 久久九九热精品免费| 精品久久久久久久人妻蜜臀av| 美女免费视频网站| 简卡轻食公司| 亚洲av五月六月丁香网| 亚洲乱码一区二区免费版| 久久婷婷人人爽人人干人人爱| 中文字幕熟女人妻在线| 亚洲图色成人| 欧美最新免费一区二区三区| 少妇的逼好多水| 色视频www国产| 看十八女毛片水多多多| or卡值多少钱| 色尼玛亚洲综合影院| 婷婷亚洲欧美| 午夜福利18| 色哟哟·www| 性色avwww在线观看| 亚洲成a人片在线一区二区| 中文字幕久久专区| 麻豆成人av在线观看| 国产不卡一卡二| 天美传媒精品一区二区| 免费在线观看成人毛片| 联通29元200g的流量卡| 亚洲中文字幕一区二区三区有码在线看| 亚洲 国产 在线| 国产乱人伦免费视频| 久久精品人妻少妇| 长腿黑丝高跟| 亚洲七黄色美女视频| 国产精品一区www在线观看 | 精品人妻偷拍中文字幕| 日韩欧美三级三区| 少妇猛男粗大的猛烈进出视频 | 精品不卡国产一区二区三区| 少妇人妻一区二区三区视频| 中文字幕熟女人妻在线| 国产美女午夜福利| 天堂av国产一区二区熟女人妻| 日本-黄色视频高清免费观看| 成年女人永久免费观看视频| 天堂影院成人在线观看| 久久国产乱子免费精品| 成人av在线播放网站| 亚洲最大成人中文| 我的女老师完整版在线观看| 韩国av在线不卡| 亚洲经典国产精华液单| 国产一区二区激情短视频| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产真实乱freesex| 又爽又黄a免费视频| 国产毛片a区久久久久| 97碰自拍视频| 一a级毛片在线观看| 亚洲无线在线观看| 国产精品精品国产色婷婷| 久久午夜亚洲精品久久| 色视频www国产| 午夜精品一区二区三区免费看| 黄色配什么色好看| 97碰自拍视频| 97超视频在线观看视频| 欧美zozozo另类| 无人区码免费观看不卡| 久久国产精品人妻蜜桃| 亚洲色图av天堂| av中文乱码字幕在线| 国产精品久久久久久精品电影| 一个人免费在线观看电影| 午夜爱爱视频在线播放| 日韩国内少妇激情av| 色吧在线观看| 一级黄片播放器| 日韩精品中文字幕看吧| 亚州av有码| 亚洲精品456在线播放app | 嫁个100分男人电影在线观看| 99久久精品一区二区三区| 女同久久另类99精品国产91| 男女边吃奶边做爰视频| 在线播放国产精品三级| 一进一出抽搐gif免费好疼| 亚洲av熟女| av在线老鸭窝| 人人妻人人澡欧美一区二区| 日本欧美国产在线视频| 免费av观看视频| 国产av在哪里看| 男插女下体视频免费在线播放| 中文字幕av在线有码专区| 亚洲成人中文字幕在线播放| 中文字幕熟女人妻在线| 精华霜和精华液先用哪个| 国产成人福利小说| 欧美激情在线99| 日韩 亚洲 欧美在线| av在线老鸭窝| 国产精品久久久久久久电影| 国产高清激情床上av| 亚洲av第一区精品v没综合| 韩国av在线不卡| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在 | 日本-黄色视频高清免费观看| 国产蜜桃级精品一区二区三区| 少妇裸体淫交视频免费看高清| 热99re8久久精品国产| 日韩欧美在线乱码| 精华霜和精华液先用哪个| 午夜福利视频1000在线观看| 天天一区二区日本电影三级| 舔av片在线| 亚洲第一电影网av| 亚洲精品一区av在线观看| 国产精品嫩草影院av在线观看 | 神马国产精品三级电影在线观看| 一边摸一边抽搐一进一小说| 国产精品99久久久久久久久| 麻豆国产97在线/欧美| 免费看日本二区| 日韩一本色道免费dvd| 嫩草影院入口| 国产亚洲精品综合一区在线观看| 久久精品综合一区二区三区| 国内久久婷婷六月综合欲色啪| 男人舔奶头视频| 久久国内精品自在自线图片| 伦理电影大哥的女人| 亚洲熟妇熟女久久| 国内精品久久久久精免费| 人人妻人人澡欧美一区二区| 国产日本99.免费观看| 亚洲专区国产一区二区| 一级毛片久久久久久久久女| 波多野结衣巨乳人妻| 国产伦精品一区二区三区视频9| 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| а√天堂www在线а√下载| 国产成人a区在线观看| 久久这里只有精品中国| 久久久国产成人精品二区| 1000部很黄的大片| 欧美xxxx黑人xx丫x性爽| 99久久精品一区二区三区| 亚州av有码| 极品教师在线视频| 人妻夜夜爽99麻豆av| www.www免费av| 韩国av一区二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 俺也久久电影网| 人妻制服诱惑在线中文字幕| 亚洲av日韩精品久久久久久密| 精品午夜福利视频在线观看一区| 欧美+亚洲+日韩+国产| 嫩草影院入口| 两性午夜刺激爽爽歪歪视频在线观看| 色哟哟·www| 淫妇啪啪啪对白视频| 黄片wwwwww| 亚洲人成伊人成综合网2020| 日韩中文字幕欧美一区二区| 国产精品久久久久久av不卡| 欧美绝顶高潮抽搐喷水| 日韩强制内射视频| 国产探花在线观看一区二区| 国产男人的电影天堂91| 乱码一卡2卡4卡精品| 少妇被粗大猛烈的视频| 成人二区视频| 亚洲精品在线观看二区| x7x7x7水蜜桃| 免费在线观看影片大全网站| 欧美色视频一区免费| 国产高清不卡午夜福利| 亚洲av.av天堂| 日本在线视频免费播放| eeuss影院久久| 一级av片app| 12—13女人毛片做爰片一| 日本欧美国产在线视频| 在线观看免费视频日本深夜| 久久亚洲真实| 成人高潮视频无遮挡免费网站| 少妇裸体淫交视频免费看高清| 国产av麻豆久久久久久久| 久久久久久伊人网av| 亚洲国产精品久久男人天堂| 3wmmmm亚洲av在线观看| 亚洲国产精品久久男人天堂| 黄片wwwwww| 老司机午夜福利在线观看视频| 大型黄色视频在线免费观看| 精品一区二区三区视频在线| 三级国产精品欧美在线观看| 国产成人影院久久av| 搞女人的毛片| 69av精品久久久久久| 日韩欧美在线乱码| 亚洲av日韩精品久久久久久密| 午夜精品一区二区三区免费看| 亚洲一级一片aⅴ在线观看| 亚洲男人的天堂狠狠| 亚洲乱码一区二区免费版| 国产视频内射| 午夜福利在线在线| 男人狂女人下面高潮的视频| 国产av在哪里看| 亚洲欧美日韩高清专用| 91久久精品电影网| 给我免费播放毛片高清在线观看| 免费观看在线日韩| 日韩高清综合在线| 99久久九九国产精品国产免费| 两个人视频免费观看高清| 亚洲午夜理论影院| 听说在线观看完整版免费高清| 国产色婷婷99| 一进一出好大好爽视频| 高清在线国产一区| 日本 av在线| 男人狂女人下面高潮的视频| 精品人妻视频免费看| 91精品国产九色| ponron亚洲| 精品一区二区三区视频在线| 色哟哟哟哟哟哟| 国产淫片久久久久久久久| 国产视频一区二区在线看| 午夜老司机福利剧场| 波野结衣二区三区在线| 免费搜索国产男女视频| 在线a可以看的网站| 91在线精品国自产拍蜜月| 亚洲欧美日韩高清在线视频| 五月玫瑰六月丁香| 淫妇啪啪啪对白视频| 九色成人免费人妻av| 国产成人a区在线观看| АⅤ资源中文在线天堂| 九色国产91popny在线| 久久久久久久久大av| 日本三级黄在线观看| 婷婷色综合大香蕉| 成熟少妇高潮喷水视频| 日本 av在线| 国产精品久久电影中文字幕| 少妇丰满av| 免费看日本二区| 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看| 久久精品国产亚洲av香蕉五月| 91精品国产九色| 国产主播在线观看一区二区| 麻豆成人午夜福利视频| 黄色一级大片看看| 久久久久久久久久黄片| 亚洲人成网站在线播放欧美日韩| 97人妻精品一区二区三区麻豆| 欧美成人a在线观看| 在线看三级毛片| 欧美最新免费一区二区三区| 欧美黑人欧美精品刺激| 午夜精品久久久久久毛片777| 韩国av一区二区三区四区| 九九久久精品国产亚洲av麻豆| 99久久精品国产国产毛片| 日韩欧美国产在线观看| 国产爱豆传媒在线观看| 日韩欧美精品免费久久| 少妇的逼水好多| 18禁黄网站禁片免费观看直播| 一级毛片久久久久久久久女| 中文字幕久久专区| 少妇被粗大猛烈的视频| 成人亚洲精品av一区二区| 久久这里只有精品中国| 久久久久精品国产欧美久久久| 美女免费视频网站| 级片在线观看| 亚洲av免费在线观看| 欧美潮喷喷水| 国国产精品蜜臀av免费| 国产免费男女视频| 日本a在线网址| 在线观看66精品国产| 少妇猛男粗大的猛烈进出视频 | 亚洲精品国产成人久久av| 中国美白少妇内射xxxbb| 我要看日韩黄色一级片| 国产毛片a区久久久久| 亚洲欧美日韩卡通动漫| 中文亚洲av片在线观看爽| 淫秽高清视频在线观看| 国产午夜精品论理片| 我要搜黄色片| 亚洲国产精品合色在线| 成人亚洲精品av一区二区| 国产男靠女视频免费网站| 成人鲁丝片一二三区免费| 欧美区成人在线视频| 男插女下体视频免费在线播放| 久久热精品热| 免费在线观看成人毛片| 欧美bdsm另类| 九九爱精品视频在线观看| 亚洲精品色激情综合| 天堂网av新在线| 我要看日韩黄色一级片| 国产精品亚洲美女久久久| 亚洲五月天丁香| 黄色配什么色好看| 嫩草影院入口| 22中文网久久字幕| 少妇人妻精品综合一区二区 | 少妇人妻精品综合一区二区 | 亚洲在线观看片| 精品午夜福利在线看| 国产av一区在线观看免费| 丰满的人妻完整版| 亚洲av成人精品一区久久| a级毛片a级免费在线| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 成人特级黄色片久久久久久久| 中文字幕精品亚洲无线码一区| av中文乱码字幕在线| 日韩欧美在线二视频| 日日撸夜夜添| 赤兔流量卡办理| 偷拍熟女少妇极品色| 黄色丝袜av网址大全| 床上黄色一级片| 色综合婷婷激情| 丝袜美腿在线中文| 久久香蕉精品热| 一级黄片播放器| 国产精品无大码| 午夜激情福利司机影院| 在线观看av片永久免费下载| 国产亚洲欧美98| avwww免费| 国产av一区在线观看免费| 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 性插视频无遮挡在线免费观看| 久久香蕉精品热| 亚洲人成伊人成综合网2020| 91av网一区二区| 国产精品一区二区免费欧美| 国产一区二区亚洲精品在线观看| 老熟妇乱子伦视频在线观看| 桃红色精品国产亚洲av| 成人特级黄色片久久久久久久| 非洲黑人性xxxx精品又粗又长| 人人妻人人澡欧美一区二区| 成人三级黄色视频| 蜜桃亚洲精品一区二区三区| 亚洲美女视频黄频| 少妇的逼水好多| 观看美女的网站| 亚洲第一区二区三区不卡| 能在线免费观看的黄片| 男女下面进入的视频免费午夜| 中文字幕人妻熟人妻熟丝袜美| 露出奶头的视频| 国产免费男女视频| 给我免费播放毛片高清在线观看| 日韩一本色道免费dvd| 好男人在线观看高清免费视频| 精品人妻1区二区| 亚洲欧美精品综合久久99| 五月玫瑰六月丁香| 日韩在线高清观看一区二区三区 | 亚洲人成网站高清观看| 国产一区二区亚洲精品在线观看| 网址你懂的国产日韩在线| 嫁个100分男人电影在线观看| 亚洲精品一区av在线观看| 蜜桃久久精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 欧美人与善性xxx| 国产精华一区二区三区| 日韩人妻高清精品专区| 亚洲精品乱码久久久v下载方式| 国产精品国产三级国产av玫瑰| 国产亚洲精品av在线| 国产精品久久久久久av不卡| 伦理电影大哥的女人| 国产真实乱freesex| 午夜影院日韩av| 黄色丝袜av网址大全| 变态另类丝袜制服| 国产乱人伦免费视频| 中国美女看黄片| 日韩精品青青久久久久久| 特大巨黑吊av在线直播| АⅤ资源中文在线天堂| 伦精品一区二区三区| 精品欧美国产一区二区三| av中文乱码字幕在线| av在线老鸭窝| 狂野欧美激情性xxxx在线观看| 国产高清不卡午夜福利| 能在线免费观看的黄片| 国产免费av片在线观看野外av| 免费人成在线观看视频色| 男女边吃奶边做爰视频| 伊人久久精品亚洲午夜| 久久午夜福利片| 午夜福利高清视频| 国产69精品久久久久777片| 天堂动漫精品| 韩国av一区二区三区四区| 嫁个100分男人电影在线观看| 婷婷精品国产亚洲av在线| 国产高清三级在线| 国产伦人伦偷精品视频| ponron亚洲| 99久久成人亚洲精品观看| 超碰av人人做人人爽久久| 麻豆久久精品国产亚洲av| 亚洲久久久久久中文字幕| 狂野欧美激情性xxxx在线观看| 夜夜看夜夜爽夜夜摸| 亚洲欧美精品综合久久99| 国产一区二区亚洲精品在线观看| 日本成人三级电影网站| 国产视频内射| 国产欧美日韩精品亚洲av| 亚洲经典国产精华液单| 午夜精品在线福利| 亚洲精品在线观看二区| 亚洲avbb在线观看| 丰满人妻一区二区三区视频av| 大又大粗又爽又黄少妇毛片口| 亚洲内射少妇av| 国产精品人妻久久久久久| 欧美成人a在线观看| 成年女人看的毛片在线观看| 久99久视频精品免费| 亚洲自拍偷在线| 欧美国产日韩亚洲一区| 国产亚洲91精品色在线| 91狼人影院| 久久久久国内视频| 免费在线观看影片大全网站| 两人在一起打扑克的视频| 日本色播在线视频| 国产精品永久免费网站| 久久久久久伊人网av| 国产成人a区在线观看| 久久久久久大精品| 色综合婷婷激情| 久久久精品欧美日韩精品| 成人国产综合亚洲| 欧美成人a在线观看| 人妻制服诱惑在线中文字幕|