• <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ù)
    熟妇人妻不卡中文字幕| 亚洲精品第二区| 国产精品无大码| 九草在线视频观看| 舔av片在线| 国产精品秋霞免费鲁丝片| 日产精品乱码卡一卡2卡三| 边亲边吃奶的免费视频| 国产男女超爽视频在线观看| 国产欧美日韩精品一区二区| 国产无遮挡羞羞视频在线观看| 国产av国产精品国产| 97超视频在线观看视频| 2022亚洲国产成人精品| .国产精品久久| 亚洲在久久综合| 国产精品久久久久成人av| 国产男女内射视频| 2022亚洲国产成人精品| 只有这里有精品99| www.av在线官网国产| 国产精品久久久久久精品古装| 国产中年淑女户外野战色| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 亚洲av福利一区| 好男人视频免费观看在线| 美女视频免费永久观看网站| 一个人免费看片子| 亚洲中文av在线| 大陆偷拍与自拍| 高清黄色对白视频在线免费看 | 国产精品秋霞免费鲁丝片| 日本爱情动作片www.在线观看| 女性被躁到高潮视频| videos熟女内射| 人人妻人人澡人人爽人人夜夜| 亚洲美女黄色视频免费看| 97在线视频观看| 亚洲国产色片| av视频免费观看在线观看| 欧美xxⅹ黑人| 国产精品人妻久久久影院| 性色avwww在线观看| 亚洲av在线观看美女高潮| 伦理电影大哥的女人| 亚洲精品日韩av片在线观看| 嫩草影院新地址| 久久精品久久久久久久性| 少妇的逼水好多| av天堂中文字幕网| 99久久精品热视频| 啦啦啦视频在线资源免费观看| 欧美 日韩 精品 国产| 欧美激情国产日韩精品一区| 国产精品偷伦视频观看了| 国产视频首页在线观看| 少妇熟女欧美另类| 欧美高清成人免费视频www| 男女啪啪激烈高潮av片| 久久人人爽av亚洲精品天堂 | 精品国产乱码久久久久久小说| 久久精品久久精品一区二区三区| 亚洲精品乱码久久久v下载方式| 不卡视频在线观看欧美| 99热国产这里只有精品6| 亚洲欧美日韩卡通动漫| 老熟女久久久| 黑丝袜美女国产一区| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 女的被弄到高潮叫床怎么办| 黄色日韩在线| 亚州av有码| 99久久精品国产国产毛片| 国产色爽女视频免费观看| 国产av精品麻豆| 色婷婷av一区二区三区视频| 在线 av 中文字幕| 小蜜桃在线观看免费完整版高清| 男人舔奶头视频| 多毛熟女@视频| 最近的中文字幕免费完整| 国产av一区二区精品久久 | 国产在视频线精品| 大又大粗又爽又黄少妇毛片口| 亚洲欧美一区二区三区黑人 | 热99国产精品久久久久久7| av在线老鸭窝| 亚洲av日韩在线播放| 国产亚洲91精品色在线| 99热6这里只有精品| 赤兔流量卡办理| 亚洲精品亚洲一区二区| 人人妻人人添人人爽欧美一区卜 | 欧美精品一区二区大全| 中国国产av一级| 日韩一区二区三区影片| 黄色配什么色好看| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 日韩视频在线欧美| 久久99精品国语久久久| 91久久精品国产一区二区成人| 欧美极品一区二区三区四区| 欧美三级亚洲精品| 一级a做视频免费观看| 国产黄色视频一区二区在线观看| 亚洲最大成人中文| 一级毛片 在线播放| 夜夜骑夜夜射夜夜干| 成人美女网站在线观看视频| 国产黄色免费在线视频| 校园人妻丝袜中文字幕| 国产91av在线免费观看| 亚洲第一区二区三区不卡| 2018国产大陆天天弄谢| 美女高潮的动态| 在线精品无人区一区二区三 | 夜夜爽夜夜爽视频| 国产精品国产三级专区第一集| 欧美性感艳星| 直男gayav资源| 国产精品久久久久久久电影| 精品一区在线观看国产| 免费看日本二区| 色综合色国产| 国产成人免费观看mmmm| 丝袜脚勾引网站| kizo精华| 国产成人91sexporn| 日本一二三区视频观看| 亚洲国产精品999| 婷婷色综合大香蕉| 六月丁香七月| 男女无遮挡免费网站观看| 午夜日本视频在线| 久久青草综合色| 久久久久久久久久久丰满| 肉色欧美久久久久久久蜜桃| 偷拍熟女少妇极品色| 特大巨黑吊av在线直播| 大话2 男鬼变身卡| 久久青草综合色| 伊人久久国产一区二区| 午夜福利在线在线| 亚洲国产最新在线播放| 欧美高清成人免费视频www| 亚洲成人中文字幕在线播放| 男女国产视频网站| 毛片女人毛片| av黄色大香蕉| 亚洲四区av| 国产精品女同一区二区软件| 国产精品嫩草影院av在线观看| h视频一区二区三区| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 亚洲美女搞黄在线观看| 三级国产精品欧美在线观看| 91久久精品国产一区二区三区| 国产视频内射| 免费av不卡在线播放| 日韩三级伦理在线观看| 麻豆精品久久久久久蜜桃| 久久久久久人妻| av在线app专区| 99久久精品一区二区三区| 三级国产精品欧美在线观看| 亚洲自偷自拍三级| 久久6这里有精品| 性色av一级| 久久鲁丝午夜福利片| 亚洲怡红院男人天堂| 成年人午夜在线观看视频| 欧美老熟妇乱子伦牲交| 成人国产av品久久久| 少妇人妻久久综合中文| 亚洲国产色片| 黄色日韩在线| 免费久久久久久久精品成人欧美视频 | 在线观看国产h片| 男人狂女人下面高潮的视频| 在线看a的网站| 久久影院123| 亚洲精品国产av蜜桃| 久久人妻熟女aⅴ| 亚洲内射少妇av| av国产久精品久网站免费入址| 久久精品久久久久久久性| 精品视频人人做人人爽| 交换朋友夫妻互换小说| 高清av免费在线| 一区二区三区乱码不卡18| 极品教师在线视频| 日韩成人伦理影院| 两个人的视频大全免费| 亚洲内射少妇av| 国产亚洲精品久久久com| 内地一区二区视频在线| 免费av不卡在线播放| 亚洲天堂av无毛| 国产亚洲最大av| 精品一区二区三区视频在线| 国产午夜精品一二区理论片| 亚洲精品久久久久久婷婷小说| 日本猛色少妇xxxxx猛交久久| 毛片女人毛片| 秋霞伦理黄片| 2022亚洲国产成人精品| 精品久久久精品久久久| 色5月婷婷丁香| 最近最新中文字幕免费大全7| 青春草国产在线视频| 插逼视频在线观看| 国产成人免费无遮挡视频| 成人一区二区视频在线观看| 日本午夜av视频| 女性被躁到高潮视频| 看免费成人av毛片| 成人18禁高潮啪啪吃奶动态图 | 99久国产av精品国产电影| 人人妻人人爽人人添夜夜欢视频 | 一级毛片久久久久久久久女| 国产欧美日韩精品一区二区| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 欧美xxⅹ黑人| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 十分钟在线观看高清视频www | 极品少妇高潮喷水抽搐| 日本色播在线视频| 秋霞在线观看毛片| 亚洲av欧美aⅴ国产| 亚洲成人手机| 欧美激情国产日韩精品一区| 久久人人爽人人爽人人片va| 国产乱人偷精品视频| 自拍偷自拍亚洲精品老妇| 国产精品久久久久成人av| 亚洲欧美日韩东京热| 自拍欧美九色日韩亚洲蝌蚪91 | 黄色日韩在线| 一级毛片久久久久久久久女| 另类亚洲欧美激情| 国产男女内射视频| 国产成人精品一,二区| 3wmmmm亚洲av在线观看| 黄色视频在线播放观看不卡| 身体一侧抽搐| 一本久久精品| 亚洲欧美日韩东京热| 亚洲国产欧美人成| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 免费看av在线观看网站| 免费高清在线观看视频在线观看| 2021少妇久久久久久久久久久| 日本黄色日本黄色录像| 国产高清国产精品国产三级 | 亚洲国产精品一区三区| 午夜免费男女啪啪视频观看| 黑人高潮一二区| 免费黄网站久久成人精品| freevideosex欧美| 97精品久久久久久久久久精品| 日本免费在线观看一区| 午夜福利视频精品| 秋霞伦理黄片| 国产精品无大码| 亚洲,一卡二卡三卡| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 国产成人精品婷婷| 在线 av 中文字幕| 精品少妇黑人巨大在线播放| 中国三级夫妇交换| 免费在线观看成人毛片| 国产黄片美女视频| 国产一区二区三区av在线| 国产亚洲91精品色在线| 99热网站在线观看| 亚洲国产精品成人久久小说| 国产精品一二三区在线看| 最黄视频免费看| 国产精品无大码| 噜噜噜噜噜久久久久久91| 丰满少妇做爰视频| 啦啦啦啦在线视频资源| 国产精品蜜桃在线观看| 中文乱码字字幕精品一区二区三区| 18禁动态无遮挡网站| 男人添女人高潮全过程视频| 欧美人与善性xxx| 精品熟女少妇av免费看| 黄片wwwwww| 午夜老司机福利剧场| 国产欧美日韩一区二区三区在线 | 亚洲av.av天堂| 亚洲精品国产av成人精品| 日韩视频在线欧美| 美女脱内裤让男人舔精品视频| 97超视频在线观看视频| 亚洲欧美中文字幕日韩二区| 久热这里只有精品99| 大话2 男鬼变身卡| 亚洲av电影在线观看一区二区三区| 国产精品国产三级国产专区5o| 国产在线视频一区二区| 亚洲美女黄色视频免费看| 大片免费播放器 马上看| 在线观看人妻少妇| 91久久精品电影网| 久久人人爽av亚洲精品天堂 | 26uuu在线亚洲综合色| 99久久综合免费| 99久久精品一区二区三区| 中文欧美无线码| 日韩免费高清中文字幕av| 国产精品国产av在线观看| 汤姆久久久久久久影院中文字幕| 国产高清有码在线观看视频| 又爽又黄a免费视频| 日韩av不卡免费在线播放| 日本爱情动作片www.在线观看| 久久国产精品男人的天堂亚洲 | 日韩中文字幕视频在线看片 | 亚洲精品乱码久久久久久按摩| 久久久久人妻精品一区果冻| 国产v大片淫在线免费观看| 亚洲经典国产精华液单| 天堂8中文在线网| 国产 精品1| 亚洲国产精品999| 亚洲国产av新网站| 亚洲国产最新在线播放| 亚洲欧美日韩东京热| 18禁裸乳无遮挡动漫免费视频| 中文乱码字字幕精品一区二区三区| 久久毛片免费看一区二区三区| 又大又黄又爽视频免费| 啦啦啦啦在线视频资源| 夜夜看夜夜爽夜夜摸| 夫妻午夜视频| 乱系列少妇在线播放| 国产一区有黄有色的免费视频| 欧美日韩在线观看h| 国产高清有码在线观看视频| 水蜜桃什么品种好| 亚洲国产高清在线一区二区三| 亚洲av成人精品一二三区| 国产精品99久久99久久久不卡 | 啦啦啦啦在线视频资源| 丰满少妇做爰视频| 在线播放无遮挡| 女性被躁到高潮视频| 亚洲人与动物交配视频| 精品人妻一区二区三区麻豆| 韩国高清视频一区二区三区| 久热久热在线精品观看| 九九爱精品视频在线观看| 免费大片黄手机在线观看| 深夜a级毛片| 中国国产av一级| 18禁裸乳无遮挡免费网站照片| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看 | 特大巨黑吊av在线直播| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 99久久精品国产国产毛片| 一个人免费看片子| 久热这里只有精品99| 国产一区有黄有色的免费视频| 纵有疾风起免费观看全集完整版| 最近中文字幕高清免费大全6| 三级国产精品片| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| 亚洲精品第二区| 一个人看视频在线观看www免费| 久久久午夜欧美精品| a级毛片免费高清观看在线播放| 一级黄片播放器| 天堂中文最新版在线下载| 成人18禁高潮啪啪吃奶动态图 | 大香蕉久久网| 男女国产视频网站| av线在线观看网站| 在线天堂最新版资源| 欧美少妇被猛烈插入视频| 一级av片app| 欧美最新免费一区二区三区| 18禁裸乳无遮挡动漫免费视频| 国产精品99久久99久久久不卡 | 交换朋友夫妻互换小说| 国产精品精品国产色婷婷| 如何舔出高潮| 91久久精品电影网| 欧美成人一区二区免费高清观看| 国产成人a区在线观看| 亚洲精品日韩在线中文字幕| 国内精品宾馆在线| 国产成人freesex在线| 国国产精品蜜臀av免费| 高清在线视频一区二区三区| 2022亚洲国产成人精品| 亚洲欧美日韩卡通动漫| 另类亚洲欧美激情| 亚洲av男天堂| 午夜精品国产一区二区电影| 国产精品一区二区性色av| 久久 成人 亚洲| 少妇的逼水好多| 精品国产乱码久久久久久小说| 亚洲综合精品二区| 精品少妇黑人巨大在线播放| 国产精品免费大片| 国产免费又黄又爽又色| 日本wwww免费看| 干丝袜人妻中文字幕| 日韩一本色道免费dvd| av不卡在线播放| h日本视频在线播放| 久久久久久久精品精品| 欧美高清成人免费视频www| 高清av免费在线| 亚洲美女搞黄在线观看| 日产精品乱码卡一卡2卡三| 日韩免费高清中文字幕av| 精品少妇黑人巨大在线播放| 成人18禁高潮啪啪吃奶动态图 | 亚洲熟女精品中文字幕| 下体分泌物呈黄色| 免费av不卡在线播放| 青青草视频在线视频观看| 秋霞伦理黄片| av免费在线看不卡| 人妻系列 视频| 欧美xxⅹ黑人| 欧美zozozo另类| 日韩在线高清观看一区二区三区| 一本色道久久久久久精品综合| 亚洲国产精品成人久久小说| 欧美97在线视频| 国产熟女欧美一区二区| 99热网站在线观看| 国产精品一及| 下体分泌物呈黄色| av线在线观看网站| 国产在线免费精品| 波野结衣二区三区在线| 深夜a级毛片| 偷拍熟女少妇极品色| 干丝袜人妻中文字幕| 香蕉精品网在线| 网址你懂的国产日韩在线| 91aial.com中文字幕在线观看| 偷拍熟女少妇极品色| 亚洲欧美成人综合另类久久久| 黑人高潮一二区| 大香蕉97超碰在线| 熟女av电影| 中国三级夫妇交换| 亚洲不卡免费看| 纵有疾风起免费观看全集完整版| 18禁在线播放成人免费| 又爽又黄a免费视频| 日本黄色片子视频| 伦理电影大哥的女人| 亚洲国产精品一区三区| 久热这里只有精品99| 少妇人妻久久综合中文| 久热久热在线精品观看| 欧美精品一区二区大全| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩东京热| 久久影院123| 成年女人在线观看亚洲视频| 亚洲精品乱久久久久久| 直男gayav资源| 最黄视频免费看| 人妻系列 视频| 亚洲一区二区三区欧美精品| 精品99又大又爽又粗少妇毛片| 亚洲一级一片aⅴ在线观看| 亚洲人成网站在线观看播放| 亚洲,一卡二卡三卡| 日本免费在线观看一区| 人妻系列 视频| 亚洲精华国产精华液的使用体验| 男女免费视频国产| 我要看黄色一级片免费的| 伦理电影大哥的女人| 男男h啪啪无遮挡| 女人久久www免费人成看片| 91在线精品国自产拍蜜月| 免费大片黄手机在线观看| 在线观看免费高清a一片| 麻豆成人av视频| 97在线视频观看| 免费看不卡的av| 99九九线精品视频在线观看视频| 伦精品一区二区三区| 我要看日韩黄色一级片| h视频一区二区三区| 亚洲欧洲日产国产| 秋霞伦理黄片| 国产精品免费大片| av视频免费观看在线观看| 又粗又硬又长又爽又黄的视频| 国产高清三级在线| 又大又黄又爽视频免费| 欧美成人精品欧美一级黄| 亚洲欧美日韩无卡精品| 久久99热6这里只有精品| 汤姆久久久久久久影院中文字幕| 国产成人aa在线观看| 男的添女的下面高潮视频| 精品人妻熟女av久视频| 婷婷色麻豆天堂久久| 蜜桃久久精品国产亚洲av| 少妇的逼好多水| www.av在线官网国产| 老司机影院成人| 久久国产乱子免费精品| 一级爰片在线观看| 观看美女的网站| 肉色欧美久久久久久久蜜桃| 成年免费大片在线观看| 国产精品国产三级国产av玫瑰| 日韩在线高清观看一区二区三区| 最后的刺客免费高清国语| 汤姆久久久久久久影院中文字幕| 极品少妇高潮喷水抽搐| 我要看日韩黄色一级片| 高清午夜精品一区二区三区| 天堂俺去俺来也www色官网| 亚洲av.av天堂| 亚洲无线观看免费| 婷婷色综合www| av播播在线观看一区| 新久久久久国产一级毛片| 国产精品精品国产色婷婷| 一级毛片黄色毛片免费观看视频| 亚洲成人一二三区av| 777米奇影视久久| 国产在线视频一区二区| 五月天丁香电影| 精品国产乱码久久久久久小说| 国产欧美日韩精品一区二区| 一二三四中文在线观看免费高清| 男女啪啪激烈高潮av片| 久久久久久久久久久丰满| 狂野欧美激情性xxxx在线观看| 黑丝袜美女国产一区| 一级黄片播放器| 成人国产av品久久久| 午夜福利视频精品| 精品人妻偷拍中文字幕| 99久久人妻综合| 噜噜噜噜噜久久久久久91| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久精品古装| 国产欧美另类精品又又久久亚洲欧美| 亚洲最大成人中文| 亚洲欧美一区二区三区国产| 亚洲av男天堂| 一个人看视频在线观看www免费| 国产老妇伦熟女老妇高清| 欧美bdsm另类| 亚洲精华国产精华液的使用体验| 免费观看的影片在线观看| 亚洲天堂av无毛| 一级毛片aaaaaa免费看小| 亚洲婷婷狠狠爱综合网| 国产黄频视频在线观看| 成人美女网站在线观看视频| 国产av国产精品国产| 久久热精品热| 乱码一卡2卡4卡精品| 一级毛片久久久久久久久女| 高清黄色对白视频在线免费看 | 成人免费观看视频高清| 99热全是精品| 一本色道久久久久久精品综合| 久久久久国产网址| 久久精品夜色国产| 久久久久视频综合| 九草在线视频观看| 免费在线观看成人毛片| 中国三级夫妇交换| 岛国毛片在线播放| 亚洲精品成人av观看孕妇| 一级毛片 在线播放| 国产在线视频一区二区| 狂野欧美白嫩少妇大欣赏| 中国三级夫妇交换| 色综合色国产| 亚洲欧美日韩另类电影网站 | 亚洲人成网站在线播| 久久ye,这里只有精品| 寂寞人妻少妇视频99o| 丝瓜视频免费看黄片| 天堂8中文在线网| 日韩欧美 国产精品| 国产深夜福利视频在线观看| 亚洲av电影在线观看一区二区三区| 国产免费一级a男人的天堂| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 制服丝袜香蕉在线| 亚洲精品乱码久久久v下载方式| 一级毛片久久久久久久久女|