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

    土剪破壞次聲監(jiān)測試驗研究1)

    2020-06-10 06:22:20陳吉龍張少杰周儒夫韋方強(qiáng)
    力學(xué)學(xué)報 2020年3期
    關(guān)鍵詞:土樣滑坡剪切

    陳 喬 張 闊 陳吉龍 張少杰 呂 君? 周儒夫 韋方強(qiáng),3)

    ?(礦山地質(zhì)災(zāi)害成災(zāi)機(jī)理與防控重點實驗室,西安 710054)

    ?(中國科學(xué)院重慶綠色智能技術(shù)研究院,重慶 400714)

    ??(四川省川建勘察設(shè)計院,成都 610041)

    ??(中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,成都 610041)

    ???(中國科學(xué)院聲學(xué)研究所,北京 100190)

    ???(荊江水文水資源勘測局,湖北荊州 434099)

    引言

    滑坡是一種全球性的地質(zhì)災(zāi)害,它對世界上很多國家都造成了巨大的危害,促使這些國家的研究人員不斷探索滑坡災(zāi)害監(jiān)測和預(yù)警的技術(shù)手段.追溯滑坡監(jiān)測的研究歷史,學(xué)界比較公認(rèn)的時間起點為日本學(xué)者齋藤迪孝[1]在1963 年提出滑坡預(yù)報經(jīng)驗公式時,之后的幾十年各種高新技術(shù)在滑坡監(jiān)測領(lǐng)域中得到應(yīng)用.在滑坡監(jiān)測技術(shù)的發(fā)展過程中,有不少專家學(xué)者已經(jīng)從不同角度以不同觀點進(jìn)行了系統(tǒng)的總結(jié)[2-3],并對其發(fā)展趨勢做出了展望.而就目前我國的監(jiān)測技術(shù),預(yù)警預(yù)報的成功率仍然較低,如2016 年共發(fā)生地質(zhì)災(zāi)害9710 起,而全國共成功預(yù)報地質(zhì)災(zāi)害只有676 起[4],因此開展滑坡監(jiān)測新技術(shù)的探索具有較高的參考價值.

    次聲波是一種在自然界非常常見的聲波,其頻率在0~20 Hz 之間,低頻率和長波長的特點使得次聲波具有傳播距離遠(yuǎn)、衰減小、抗干擾能力強(qiáng)及穿透能力強(qiáng)等優(yōu)點,在傳播過程中大氣等介質(zhì)對它的吸收系數(shù)很小,適用于遠(yuǎn)距離觀測,被廣泛的用于火山噴發(fā)[5-8]、地震[9-11]、泥石流 [12-15] 等自然災(zāi)害監(jiān)測中.前期研究表明巖土體在變形破裂過程中因能量釋放必然會產(chǎn)生彈性波和在空氣中傳播的聲波,滑坡過程會產(chǎn)生明顯的次聲波信號[16],因此次聲監(jiān)測可以作為判斷滑坡的一種技術(shù)手段.

    近幾年,國內(nèi)外學(xué)者圍繞基于次聲波的滑坡監(jiān)測技術(shù)展開了大量研究.在次聲滑坡監(jiān)測預(yù)警系統(tǒng)研究方面.Li 等[17]基于譜熵和支持向量機(jī)對次聲信號進(jìn)行分類,可用于滑坡的次聲實時監(jiān)測.任際周等[18-19]從軟硬件設(shè)計的角度研究了基于次聲波技術(shù)的滑坡監(jiān)測預(yù)警系統(tǒng),通過對近年來我國數(shù)十起重大滑坡災(zāi)害臨滑預(yù)警中收集到的數(shù)據(jù)進(jìn)行分析,建立了滑坡預(yù)警模型.由于該模型是利用信號分析和統(tǒng)計的方法建立起來的,缺乏對滑坡次聲波產(chǎn)生機(jī)理的考慮,因此預(yù)警模型的準(zhǔn)確度依賴于大量現(xiàn)場數(shù)據(jù)的統(tǒng)計,而各個區(qū)域滑坡因素的多樣性決定了該模型的現(xiàn)場應(yīng)用具有一定的局限性.

    在滑坡破壞過程中次聲波監(jiān)測機(jī)理研究方面.Chai 等[20]通過單軸壓縮試驗發(fā)現(xiàn)花崗巖破裂過程中也產(chǎn)生次聲信號.朱星等[21-22]以巖體穩(wěn)定性監(jiān)測為目的,對6 種典型的巖石試樣單軸全過程加載試驗過程中的次聲波信號進(jìn)行了實時監(jiān)測,并對數(shù)字次聲波信號進(jìn)行了處理與分析.徐洪等[23]對巖石變形破壞次聲異常的能量特征進(jìn)行分析,研究結(jié)果為巖石破壞前兆預(yù)警提供了重要依據(jù).賈炳等[24-30]以次聲波技術(shù)應(yīng)用于監(jiān)測、預(yù)警礦井煤巖動力災(zāi)害為目的,相繼開展了大量的研究,通過進(jìn)行煤巖試樣的單軸加載試驗,對采集的次聲信號的相對能量和頻率進(jìn)行分析,發(fā)現(xiàn)相對能量和頻率具有明顯的階段性.除了進(jìn)行實驗室物理模型試驗次聲波的研究外,朱星等[31]還對大光包滑坡側(cè)緣拉裂壁巖崩進(jìn)行次聲監(jiān)測,并分析了巖崩次聲的基本特征.上述研究結(jié)果為巖質(zhì)滑坡次聲波監(jiān)測提供了理論參考,但對土質(zhì)滑坡次聲波監(jiān)測的研究鮮有報道.基于此,考慮到滑坡的主要力學(xué)破壞方式為剪切破壞,開展土剪破壞次聲監(jiān)測試驗研究,分析了剪切過程中的微觀聲學(xué)響應(yīng)[32-34],探索土質(zhì)滑坡中次聲產(chǎn)生機(jī)理,為實現(xiàn)土質(zhì)滑坡次聲波監(jiān)測提供理論參考.

    1 土體直剪實驗次聲?力?位移聯(lián)合監(jiān)測系統(tǒng)

    在實驗中,土體直剪實驗力?聲聯(lián)合監(jiān)測系統(tǒng)主要由變尺寸直剪盒、力學(xué)數(shù)據(jù)采集模塊、聲學(xué)數(shù)據(jù)采集模塊以及位移采集模塊四部分組成,如圖1 所示.

    1.1 變尺寸直剪盒設(shè)計

    傳統(tǒng)的應(yīng)變控制式直剪儀在進(jìn)行土體剪切實驗時,采用的土樣規(guī)格一般為面積30 cm2高2 cm 的圓柱,可以實時獲取剪切應(yīng)力與剪切位移數(shù)據(jù),但不能采集聲波信號,而在此尺寸下的土樣在剪切過程中所產(chǎn)生的聲波能量又是非常微小,因此本研究需要設(shè)計一個大尺寸且尺寸可變的剪切儀器,在本論文中稱之為變尺寸直剪盒.

    如圖2 所示,整個變尺寸剪切盒包括上部剪切盒、上部擋土板、固定上部擋土板的螺紋桿、下部擋土板、導(dǎo)向輪(滾珠花鍵)、下部剪切盒、固定上下剪切盒的螺紋桿、固定上部擋土板的螺紋桿,其中的上部剪切盒、上部擋土板、下部剪切盒、下部擋土板是4 個主要結(jié)構(gòu)體.在設(shè)計變尺寸直剪盒的各個部件和整體的協(xié)調(diào)組合時,根據(jù)研究條件重點考慮了基本的剪切運動、土樣尺寸可變、防銹蝕和防變形的性能.圖3 為加工廠制作后的實際成品.傳統(tǒng)直剪儀的固有結(jié)構(gòu)不利于安裝導(dǎo)波管放大次聲信號,故本次研究所設(shè)計的變尺寸直剪儀更適合土體直剪實驗過程中的次聲?力?位移聯(lián)合監(jiān)測.

    圖2 變尺寸剪切盒整體設(shè)計圖Fig.2 Overall design of the variable size shear box

    圖3 變尺寸剪切儀成品圖Fig.3 Finished product of variable size shear

    1.2 力學(xué)數(shù)據(jù)采集模塊

    在實驗中需要對土樣剪切過程中的推力進(jìn)行實時采集,因此需要一套合適的力學(xué)數(shù)據(jù)采集系統(tǒng).本研究中使用的力學(xué)采集儀器主要為平膜盒式力學(xué)傳感器,除此之外,還有配套的數(shù)字轉(zhuǎn)換器以及數(shù)據(jù)處理軟件等,圖4 為整個系統(tǒng)的示意圖.

    圖4 力學(xué)數(shù)據(jù)采集系統(tǒng)示意圖Fig.4 Schematic diagram of mechanical data acquisition system

    1.3 次聲數(shù)據(jù)采集模塊

    次聲數(shù)據(jù)采集模塊包括剛性導(dǎo)波管、次聲傳感器、次聲數(shù)據(jù)存儲器以及相應(yīng)的數(shù)據(jù)處理軟件,其中最核心的部件當(dāng)屬次聲傳感器.本次實驗所使用的次聲傳感器來自于中國科學(xué)院聲學(xué)研究所研制的IDS2016 型次聲傳感器,該傳感器具有體積小、靈敏度高、對機(jī)械振動不敏感、頻響好的優(yōu)點.其中3 dB 平穩(wěn)帶寬保持在0.5~200 Hz 范圍內(nèi),靈敏度為50 mV/Pa,在1~100 Hz 頻率范圍內(nèi)靈敏度曲線平直,可以完全覆蓋剪切過程中產(chǎn)生次聲波的整個頻段.

    1.4 位移采集模塊

    針對傳統(tǒng)的基于傳感器的接觸式測量方法來測量剪切位移,存在儀器安裝困難的現(xiàn)狀,故結(jié)合本研究的實際情況,實驗采用自行研發(fā)的非接觸式光學(xué)測量方法結(jié)合數(shù)字圖像識別技術(shù)來進(jìn)行位移數(shù)據(jù)的采集.

    測量的原理是通過計算物體位移變化前后呈現(xiàn)的兩幅圖像之間的相關(guān)性來計算位移值.根據(jù)本實驗的實際情況,利用此測量方法時需要在上部剪切盒表面貼一張長方形白紙(白色與鋼板的深色在二值化后會形成鮮明反差)作為定位標(biāo)記,然后使用數(shù)碼攝像機(jī)錄制整個剪切盒移動過程,得到一串?dāng)?shù)字圖像序列,再利用計算機(jī)程序?qū)γ恳桓眻D像進(jìn)行預(yù)處理,包括圖形矯正、二值化處理、選取特征像素行(特征像素行指,圖像序列中的所有圖像中,作為定位標(biāo)記的像素數(shù)量都相等的那一行像素,之所以直接取行像素作為計算量,是因為在本次研究的實驗中,上部剪切盒只做水平直線運動)等步驟(流程如圖5所示).

    圖像預(yù)處理的每一步需使用相同的變換參數(shù).顯然,預(yù)處理后的圖像的特征像素行是一個一維二進(jìn)制行向量(稱為特征行向量).如果將第一張圖像的特征像素行右移(這里的移位為循環(huán)移位,即移位后末位變首位),每次右移后得到新的特征行向量序列,將這一特征行向量序列依次與其他圖像的特征行向量進(jìn)行相關(guān)性計算(見圖6),每次計算后都會得到一組相關(guān)性系數(shù),在一組相關(guān)性系數(shù)中的每一個值都對應(yīng)一個移位的像素數(shù),而移位的像素數(shù)又對應(yīng)一個位移距離,因此,每一組相關(guān)性系數(shù)序列中的最大值對應(yīng)的位移距離就是這幅圖像中被測物體在當(dāng)前時刻的位移值.

    由于位移數(shù)據(jù)的采集最初來自數(shù)碼攝像機(jī),所以本測量方法的精度主要依賴于攝像機(jī)的分辨率.本次研究采用的4K 攝像機(jī)橫向分辨率為3840 個像素,在白紙所在的焦距處整個屏幕所囊括的總距離在所有實驗中平均約為1000 mm,由此可以算得測量精度為每像素0.26 mm,在之后的研究中可以采用更高分辨率的攝像機(jī)來提高測量精度.整個流程可以通過Photoshop 圖像處理軟件和Matlab 編程實現(xiàn).

    圖5 圖像序列第1 幀和第n 幀的處理流程圖Fig.5 Processing flowchart of the first and nth frames of the image sequence

    圖6 相關(guān)性系數(shù)計算示意圖Fig.6 Schematic diagram of correlation coefficient calculation

    2 實驗材料

    本次實驗選用的土樣均來自三峽庫區(qū)[35-36]重慶市奉節(jié)縣,共采集了5 種土類,包括棕壤(黏土)、黃壤(黏土)、黃棕壤(黏土)、粗骨土(含礫石黏土)、紫色土(含礫石黏土),部分土樣如圖7 所示.其中棕壤取自云霧鄉(xiāng)羅圈村,黃壤取自安坪鄉(xiāng)新浦滑坡后緣,黃棕壤取自長安鄉(xiāng)保安村,粗骨土取自草堂鎮(zhèn)廟灣,紫色土取自康樂鎮(zhèn)平皋村,取樣點位置見圖8,重塑土樣是原狀的黃壤土經(jīng)烘干、碾碎而成,各類土的相關(guān)物理力學(xué)性質(zhì)見表1.

    圖7 土樣類別圖Fig.7 Types of soil samples

    3 實驗方案

    3.1 實驗流程

    首先,將野外取回的所有土樣分批放入烘箱中,以105?C 的溫度持續(xù)24 h 鼓風(fēng)干燥,并將烘干后的土樣研磨并過2 mm 篩備用; 接下來,按具體實驗?zāi)康?將變尺寸剪切儀的上部擋土板和下部擋土板固定在規(guī)定的位置,采用分層壓實的方法將處理好的土樣堆填于剪切盒中;然后,在土樣頂部插入剛性導(dǎo)波管,在導(dǎo)波管頂部安裝次聲傳感器探頭,將力學(xué)傳感器放置在千斤頂與上部剪切盒之間,將記錄位移用的攝像頭擺放在合適位置;最后,記錄實驗開始時間,并打開數(shù)據(jù)采集軟件,在記錄30 s 背景噪音的基礎(chǔ)上,開始按壓千斤頂為上部剪切盒加載直至土樣剪切破壞.

    圖8 土樣取樣點位置圖Fig.8 Location of soil sampling points

    表1 各類土的相關(guān)物理力學(xué)性質(zhì)Table 1 Related physical and mechanical properties of various soils

    3.2 分析方法研究

    3.2.1 背景噪聲分析與處理

    在室內(nèi)試驗過程中,會出現(xiàn)很多外在條件干擾的情況.例如開、關(guān)門和抽風(fēng)機(jī)運轉(zhuǎn)時會導(dǎo)致空氣不穩(wěn)定流通,實驗用水所發(fā)出的超低頻聲等.針對上述可控的低頻聲源,本試驗時間均選定在凌晨,關(guān)閉與試驗無直接關(guān)系的電器設(shè)備,以減少外界帶來的干擾.在預(yù)處理的基礎(chǔ)上,考慮到本次研究中關(guān)注的重點是聲波信號的次聲部分,所以需要對聲波信號進(jìn)行了低通濾波處理.在文中利用巴特沃斯濾波器對信號進(jìn)行處理,設(shè)置通帶截止頻率為20 Hz,阻帶截止頻率為25 Hz,通帶內(nèi)最大衰減為0.5 dB,阻帶衰減為40 dB.同時,利用背景噪聲強(qiáng)度較低的現(xiàn)狀,設(shè)置聲壓強(qiáng)度閾值為0.000 5 Pa,最后,對閾值處理后小波信號進(jìn)行小波逆變換.

    3.2.2 次聲事件拾取

    為了提取次聲事件,將整個時域按次聲事件分段,本文提出了一種按次聲事件由計算機(jī)程序自動劃分時域的方法.首先根據(jù)每個采樣點的信噪比初步判斷該采樣點是否位于一個次聲事件的時段內(nèi),再歸并比較密集的有效采樣點形成一個次聲事件,具體步驟如下.

    第一步,初步判斷一個采樣點是否位于一個次聲事件的時段內(nèi).

    根據(jù)離散聲波信號的能量公式

    其中,E表示信號能量;s表示聲波通過的面積,可以認(rèn)為是常量; ρ 表示聲波傳播介質(zhì)的密度,這里即指空氣的密度;v表示音速;I表示采樣點總數(shù);x(i) 表示第i個采樣點的聲壓值,可以得到單個采樣點所代表的能量

    根據(jù)功率的定義(單位時間內(nèi)的能量),本文定義離散聲波信號中第i個采樣點的功率為

    本文對“均方功率”的定義由式(1)右側(cè)去掉s,ρ 和v再開方得到,因為式中的s,ρ 和v為常量,所以去掉這3 個量不影響信號點間的相對取值且避免了對它們的測量,再開方的目的是為了使最后得到的值與聲壓值量綱相同.

    由式(3)定義一個判據(jù)式

    其中,Pj表示第j個采樣點的功率,Pmax表示前30 s背景噪聲的最大功率,L表示用于計算信噪比的時窗長度.

    式(4)左側(cè)的意義為以第i個采樣點為中心的L時長范圍內(nèi)的信號能量,右側(cè)的意義為在長度同樣為L的時段內(nèi)噪聲的最大能量,整個式子表示將局部的信號能量與相同長度噪聲能量進(jìn)行比較.針對新信號中的每一個點代入式(5) 中判斷,結(jié)果為真時,把該點標(biāo)記為1,否則標(biāo)記為0,這樣就形成一個由0 和1 組成的序列,稱為第一輪標(biāo)記序列,其長度等于原始信號.于是整個時域按信噪比被初步劃分出來.此時的劃分是不夠的,因為第一輪標(biāo)記序列中的1 值點過于零碎,需要將間隔很短的信號段再一次歸并.

    第二步,對第一輪標(biāo)記序列進(jìn)行歸并整合.這一步會用到一個所有值均為1 的行向量r作為參考向量,把這個參考向量作為時窗在第一輪標(biāo)記序列的時域中平移,并與和時窗長度相同的局部按下式做計算

    其中,ci表示判別該采樣點(第i個點)是否需要重新標(biāo)記為1 或0 的值,J表示移動時窗的采樣點序列總數(shù),sj表示第一輪標(biāo)記序列中的第j個點,??表示向上取整,??表示向下取整.

    由式(6)計算得到結(jié)果ci后,再判別是否需要重置采樣點的聲壓值,當(dāng)ci> 0 時,說明第一輪標(biāo)記序列中與時窗長度相同的局部向量與作為參考樣本的向量相互匹配,該點可以標(biāo)記為1,否則當(dāng)ci< 0 時應(yīng)該標(biāo)記為0.最后再次得到一個基于時域的由若干0 和1 分別連續(xù)分布的新序列,稱為次聲事件標(biāo)記序列,其中連續(xù)的1 指明該時段內(nèi)存在能量明顯大于背景噪音的次聲信號,而連續(xù)的0 指明該時段內(nèi)不存在有效的次聲信號.

    式(5)中的時窗長度L和式(6)中時窗長度J是兩個對最終劃分結(jié)果影響較大的計算因子.其中L應(yīng)設(shè)置一個較小值,若設(shè)置值過大,一些能量較小的次聲事件將不會被識別出.而時窗長度J影響最終次聲事件劃分出來的數(shù)量,若J取值過小,次聲事件會劃分的過密,無疑會增加后期分析的工作量,而J若取值過大,對次聲事件的劃分又太過于簡單,一些間隔很遠(yuǎn)的有效信號也被歸并為一個次聲事件.經(jīng)過反復(fù)實驗,本文認(rèn)為L=5,J=100 是合適的.最后根據(jù)式(5)和式(6),使用MATLAB 編制了對信號按次聲事件自動劃分時域的程序.

    4 試驗結(jié)果分析

    4.1 應(yīng)變軟化類黏性土的漸進(jìn)性破壞過程

    本次共完成了不同尺寸、不同土性和不同飽和度共13 組土剪切實驗(見表2).圖9 展示了所有實驗的“推力?時間”曲線和聲波信號時域圖土樣的剪切力與位移原始關(guān)系.從推力變化的整體趨勢來看,每組推力都隨時間先陡然升高,超過峰值后逐漸降低.通過對推力進(jìn)行規(guī)格化處理,統(tǒng)一在0~7000 N 范圍內(nèi),可以進(jìn)一步反映推力在時域上的變化趨勢.以具有代表性的sy3 實驗為例(如圖10),圖中實線為實測推力值,推力曲線呈現(xiàn)鋸齒狀是因為本次實驗加載方式采用的是千斤頂人工加載,所以導(dǎo)致加載不連續(xù),反映出了該類土的應(yīng)力松弛特點.

    人工加載對剪切速度的影響可以認(rèn)為是微乎其微的,在本實驗中推力與剪切帶受到的剪力幾乎相等,由牛頓第二定律

    其中,f是土體受到的剪力,單位kN;F是千斤頂?shù)耐屏?單位kN;m是土體上半部分質(zhì)量,實測sy3 的土體總質(zhì)量49.24 kg;a為位移加速度,單位為m/s2.經(jīng)計算ma的值在?0.046 6 到0.049 9 之間,遠(yuǎn)遠(yuǎn)小于推力大小,因此可大致認(rèn)為在千斤頂加載的過程中上部剪切盒保持勻速運動.

    表2 實驗組別一覽Table 2 List of experimental groups

    圖9 “推力?時間”曲線和聲波信號時域圖Fig.9 “Thrust-time”curve and time domain diagram of acoustic signal

    圖10 sy3 推力?位移曲線Fig.10 Thrust-displacement curve of sy3

    同時,推力的趨勢線隨位移變化呈現(xiàn)升高后降低的規(guī)律,當(dāng)sy3 發(fā)生剪切破壞時,首先在某一點處超過抗剪強(qiáng)度而發(fā)生破裂,然后從此破裂點開始向整個剪切面蔓延擴(kuò)展,隨著剪切面的蔓延擴(kuò)展,整個土體的抗剪強(qiáng)度也在降低,應(yīng)變局部化現(xiàn)象增強(qiáng),直至整個剪切面貫通,整體表現(xiàn)出明顯的應(yīng)變軟化特性.

    4.2 剪切過程的次聲響應(yīng)

    根據(jù)漸進(jìn)破壞理論和推力趨勢線的形態(tài),本文將其劃分為了彈性變形階段、彈塑性變形階段以及峰值后的應(yīng)變軟化階段3 個階段.為了探究與聲學(xué)信號的關(guān)系,將推力?位移曲線在時域上重新描述,繪制出的“推力?時間”曲線圖,以sy3 為例(見圖11),推力趨勢線在經(jīng)過前30 s 的背景噪音錄制階段后隨時間逐漸增大,加載后先是呈直線狀態(tài)增大,處于彈性變形階段;在74 s 后曲線斜率開始減小,曲線本身呈向下彎曲狀,處于彈塑性變形階段; 在118 s 處推力趨勢線達(dá)到峰值6.1 kN 后開始減少,隨后隨著時間推移逐漸減小,對應(yīng)著應(yīng)變軟化階段.

    sy3 在次聲事件拾取的基礎(chǔ)上獲取了26 個次聲事件(見圖11),將每個事件對應(yīng)的次聲信號均方功率峰值用虛線連接起來形成均方功率的包絡(luò)線,隨著時間推移包絡(luò)線開始逐漸上升,在105 s 處均方功率達(dá)到0.047 9 Pa 后開始下降,之后偶有起伏但都沒有超過105 s 處的峰值,并且很明顯的,該包絡(luò)線的時域形態(tài)與推力趨勢線較為相似.

    圖11 sy3 均方功率圖Fig.11 sy3 mean square power diagram

    從圖11 來看,剪切試驗過程中的時域信號在105 s 處能量達(dá)到了峰值,同時應(yīng)力曲線也達(dá)到峰值,實際上從現(xiàn)場結(jié)果來看試樣在達(dá)到應(yīng)力峰值時已經(jīng)破壞,但此時應(yīng)力并沒有大幅度下降,表明土樣還有一定殘余強(qiáng)度.破壞過程中的次聲信號強(qiáng)度比破裂后剪切面滑移所產(chǎn)生的次聲信號強(qiáng)度大,說明土樣內(nèi)部裂隙發(fā)展所產(chǎn)生的次聲信號與土樣顆粒摩擦所產(chǎn)生的次聲信號有一定的區(qū)別.

    在此基礎(chǔ)上,將所有組實驗的均方峰值與推力峰值之間的時間差進(jìn)行計算,從結(jié)果來看(見表3 和圖12),次聲信號的均方功率峰值都在推力峰值之前,提前的時間分布在2.30~47.11 s,平均值為22.95 s.該結(jié)論可為土質(zhì)滑坡次聲預(yù)警模型研究提供理論參考,具有重大意義和工程應(yīng)用價值.

    表3 各組實驗次聲事件次數(shù)與時間差計算結(jié)果Table 1 Calculation results of the number and time difference of infrasound events in each group of experiments

    圖12 所有實驗均方功率圖Fig.12 Mean square power plots for all experiments

    在“推力?位移” 曲線的3 個階段中,彈性變形階段推力趨勢線呈近直線型,斜率幾乎不變,彈塑性變形階段趨勢線的增量逐漸降低,線形向下凹,應(yīng)變軟化階段趨勢線經(jīng)過峰值后逐漸降低,但較彈性變形階段下降趨勢更為平緩.根據(jù)應(yīng)變軟化類黏性土在剪切破壞中滿足漸進(jìn)性破壞理論,剪切面破壞面的形成呈漸進(jìn)性,另外參考一些對大型直剪實驗剪切帶形態(tài)的相關(guān)研究[37],本文認(rèn)為彈性變形階段對應(yīng)于圖13 中的剪切前到初始剪切之間的狀態(tài),彈塑性變形階段和應(yīng)變軟化階段對應(yīng)初始剪切到剪切完成之間的狀態(tài).

    圖13 剪切過程示意圖(陰影區(qū)域表示壓密區(qū),1#、2#和3#三根粗直線表示鋁絲,用于表征土樣平行于剪切方向的截面的形變)Fig.13 Schematic diagram of the cutting process

    在剪切過程初期,土體整體變形以上部試樣中壓密區(qū)土體的彈性變形為主,剪切帶中的剪切變形為輔,因此推力趨勢線呈近直線狀,在此階段,土顆粒主要發(fā)生彈性變形,顆粒間接觸面相互擠壓產(chǎn)生的聲波能量較小,但隨著剪切過程的發(fā)展,土體整體變形進(jìn)入到以剪切帶塑性變形為主的階段,此階段中主要以黏土膠體團(tuán)粒本身的變形為主,這種膠體團(tuán)粒的變形比彈性變形所產(chǎn)生的聲波能量要大,從以土顆粒彈性變形為主到黏土團(tuán)粒本身變形為主,后者產(chǎn)生的聲波信號成分也逐漸超過前者,所以在宏觀上表現(xiàn)為彈塑性階段次聲信號幅度值比彈性階段更大(見圖14),均方功率峰值一般也發(fā)生在此階段.在推力峰值時,土體剪切面開始漸進(jìn)性破壞,并進(jìn)入到應(yīng)變軟化階段,在此階段剪切面從某一點發(fā)生破裂并逐步蔓延到整個剪切面,此時剪切帶的運動形式在微觀上發(fā)展成以黏土團(tuán)粒拉張破裂后土顆粒間的摩擦為主,黏土團(tuán)粒的破裂為輔,此時的信號能量也逐漸減弱,偶爾局部團(tuán)粒破裂活動較為集中時,在時域信號上還會產(chǎn)生相應(yīng)的脈沖,但包含此脈沖的次聲事件的幅值一般不會再超過彈性階段的幅值.

    圖14 不同階段聲波幅值變化情況Fig.14 Changes in sound wave amplitude at different stages

    5 結(jié)論

    (1)針對常規(guī)土體應(yīng)變式直剪儀用于聲發(fā)射實驗時聲學(xué)信號微弱的缺陷,研發(fā)了專門用于探測次聲信號的變尺寸直剪儀,設(shè)計了非接觸式位移數(shù)據(jù)的自動采集模塊,搭建了一套大尺寸土體直剪實驗力?聲?位移聯(lián)合采集系統(tǒng),為土質(zhì)滑坡次聲監(jiān)測技術(shù)室內(nèi)實驗研究奠定了基礎(chǔ).

    (2) 將聲學(xué)信號結(jié)合力學(xué)和位移數(shù)據(jù)進(jìn)行分析,發(fā)現(xiàn)次聲信號在時域上的脈沖形態(tài)與直剪實驗中的“推力?時間” 曲線具有明顯的相關(guān)性,并在應(yīng)變軟化類黏性土直剪過程中剪切帶的漸進(jìn)性破壞理論的基礎(chǔ)上,分析了剪切過程中彈性階段、彈塑性階段和應(yīng)變軟化階段的微觀聲學(xué)響應(yīng),認(rèn)為黏性土剪切破壞時發(fā)出的次聲信號主要來源于彈性階段中黏土顆粒的相互擠壓、彈塑性階段和應(yīng)變軟化階段中由黏土顆粒構(gòu)成的膠體團(tuán)粒本身的變形以及黏土膠體團(tuán)粒拉張破裂后包含其中的土顆粒間的相互摩擦三種形式.

    (3)在土剪不同階段,次聲信號的幅值大小呈現(xiàn)出應(yīng)變軟化階段< 彈性階段< 彈塑性階段的規(guī)律,同時,次聲信號的均方功率峰值都在推力峰值之前,平均提前為22.95 s,這說明將次聲技術(shù)作為土質(zhì)滑坡預(yù)警手段的具有理論可行性.

    本研究對于土質(zhì)邊坡的動態(tài)監(jiān)測和穩(wěn)定性預(yù)測都具有重要的指導(dǎo)意義,可為進(jìn)一步基于次聲的土質(zhì)滑坡預(yù)警模型研究提供一定的理論參考.

    猜你喜歡
    土樣滑坡剪切
    灌區(qū)渠道基土工程水敏性試驗研究
    檸檬酸對改良紫色土中老化銅的淋洗研究
    滑坡推力隱式解與顯式解對比分析——以河北某膨脹土滑坡為例
    寬厚板剪切線控制系統(tǒng)改進(jìn)
    山東冶金(2018年5期)2018-11-22 05:12:46
    膨脹土干濕交替作用下殘余強(qiáng)度試驗方案分析
    治淮(2018年6期)2018-01-30 11:42:44
    淺談公路滑坡治理
    北方交通(2016年12期)2017-01-15 13:52:59
    混凝土短梁斜向開裂后的有效剪切剛度與變形
    基于Fluent的滑坡入水過程數(shù)值模擬
    “監(jiān)管滑坡”比“渣土山”滑坡更可怕
    山東青年(2016年3期)2016-02-28 14:25:50
    土-混凝土接觸面剪切破壞模式分析
    精品人妻视频免费看| 欧美日韩精品成人综合77777| av国产免费在线观看| 精品久久久久久久久av| 欧美zozozo另类| 亚洲av免费高清在线观看| 极品教师在线视频| 伊人久久精品亚洲午夜| 免费观看a级毛片全部| 中文精品一卡2卡3卡4更新| 校园人妻丝袜中文字幕| 欧美3d第一页| 精品久久久久久久人妻蜜臀av| 夜夜看夜夜爽夜夜摸| 国产不卡一卡二| a级毛片a级免费在线| 午夜福利在线观看吧| 国产v大片淫在线免费观看| 少妇的逼水好多| 精品一区二区三区人妻视频| 国产黄a三级三级三级人| 午夜老司机福利剧场| 中文字幕制服av| av在线亚洲专区| 精品人妻偷拍中文字幕| 久久久久久久午夜电影| 中国美白少妇内射xxxbb| 国产老妇女一区| 色噜噜av男人的天堂激情| 偷拍熟女少妇极品色| 熟女人妻精品中文字幕| 成人午夜精彩视频在线观看| 在线观看av片永久免费下载| 国产成人精品婷婷| 在线免费十八禁| 蜜桃亚洲精品一区二区三区| 亚洲,欧美,日韩| 午夜福利高清视频| 国产视频首页在线观看| 久久精品久久久久久久性| 国产激情偷乱视频一区二区| 国产成人精品久久久久久| 国产精品久久久久久精品电影| 成人午夜高清在线视频| 日本三级黄在线观看| 日韩,欧美,国产一区二区三区 | 波多野结衣高清无吗| 国产精品一及| 在现免费观看毛片| 久久久a久久爽久久v久久| 亚洲av二区三区四区| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 少妇熟女欧美另类| 欧美变态另类bdsm刘玥| 99国产极品粉嫩在线观看| 97在线视频观看| 久久久国产成人免费| 国产成人精品一,二区 | 久久99热6这里只有精品| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| 一本久久精品| 亚洲在线自拍视频| 国产女主播在线喷水免费视频网站 | 一夜夜www| 国产亚洲欧美98| 国产精品不卡视频一区二区| 精华霜和精华液先用哪个| 观看美女的网站| 性欧美人与动物交配| 91狼人影院| 国产三级在线视频| 女人十人毛片免费观看3o分钟| 特大巨黑吊av在线直播| 人体艺术视频欧美日本| 长腿黑丝高跟| 99在线视频只有这里精品首页| 悠悠久久av| 免费av不卡在线播放| 在线免费观看不下载黄p国产| 国产成人福利小说| 精品久久久久久久久亚洲| 亚洲熟妇中文字幕五十中出| 岛国在线免费视频观看| 亚洲精品久久国产高清桃花| 国模一区二区三区四区视频| 日韩强制内射视频| 国产av不卡久久| 日韩成人av中文字幕在线观看| 美女cb高潮喷水在线观看| 99热全是精品| 欧美成人精品欧美一级黄| 日本一本二区三区精品| 国产欧美日韩精品一区二区| 精品久久久久久久久久免费视频| 成人欧美大片| 久久久久国产网址| 欧美精品一区二区大全| 久久精品国产清高在天天线| 国产一区亚洲一区在线观看| 免费人成视频x8x8入口观看| 女人被狂操c到高潮| 国产色婷婷99| 两性午夜刺激爽爽歪歪视频在线观看| 在线观看av片永久免费下载| 成人毛片60女人毛片免费| 精品不卡国产一区二区三区| 欧美一区二区精品小视频在线| 亚洲自拍偷在线| 色吧在线观看| 青青草视频在线视频观看| 男人的好看免费观看在线视频| 亚洲成a人片在线一区二区| 麻豆乱淫一区二区| 亚洲国产精品久久男人天堂| 2021天堂中文幕一二区在线观| 两个人的视频大全免费| 精品人妻视频免费看| 久久久精品94久久精品| 给我免费播放毛片高清在线观看| 国产v大片淫在线免费观看| 免费人成在线观看视频色| 一区二区三区高清视频在线| 嫩草影院精品99| 久久99热6这里只有精品| 天天一区二区日本电影三级| 97超视频在线观看视频| 美女 人体艺术 gogo| 日韩在线高清观看一区二区三区| 色综合亚洲欧美另类图片| 一本久久精品| 欧美日本视频| 国产一区二区三区av在线 | 天堂网av新在线| 国产爱豆传媒在线观看| 亚洲最大成人av| 夜夜夜夜夜久久久久| 亚洲在久久综合| 99久久人妻综合| 九草在线视频观看| 在线免费观看的www视频| 日本欧美国产在线视频| 日韩一区二区视频免费看| 国产老妇女一区| 国产一区二区在线av高清观看| 1024手机看黄色片| 国产探花极品一区二区| 久久热精品热| 久久这里只有精品中国| 日日干狠狠操夜夜爽| 又爽又黄无遮挡网站| 午夜精品一区二区三区免费看| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久久人妻蜜臀av| 少妇的逼水好多| 神马国产精品三级电影在线观看| 精品久久久久久久久av| 久久午夜亚洲精品久久| 男女那种视频在线观看| 一级黄色大片毛片| 午夜a级毛片| 中国国产av一级| 国产真实乱freesex| 晚上一个人看的免费电影| 亚洲av成人精品一区久久| 久久精品91蜜桃| 我的老师免费观看完整版| 亚洲欧美成人综合另类久久久 | 三级国产精品欧美在线观看| 国产乱人视频| 2021天堂中文幕一二区在线观| 亚洲欧洲国产日韩| 中文字幕熟女人妻在线| 国产亚洲精品久久久com| 人妻少妇偷人精品九色| 老司机福利观看| 嘟嘟电影网在线观看| 麻豆一二三区av精品| 国产精品电影一区二区三区| 看黄色毛片网站| 亚洲最大成人中文| 1000部很黄的大片| 九九久久精品国产亚洲av麻豆| 久久久久久大精品| 日韩欧美精品免费久久| 成熟少妇高潮喷水视频| 赤兔流量卡办理| 国产一级毛片七仙女欲春2| 日本色播在线视频| 99久久精品热视频| 免费看美女性在线毛片视频| 美女cb高潮喷水在线观看| 校园春色视频在线观看| 国产探花极品一区二区| 一区福利在线观看| 亚洲av一区综合| 久久精品国产鲁丝片午夜精品| 国产午夜福利久久久久久| 免费大片18禁| 亚洲欧美日韩高清在线视频| 九九久久精品国产亚洲av麻豆| 亚洲人与动物交配视频| 久久精品91蜜桃| 亚洲av一区综合| 男女那种视频在线观看| avwww免费| 亚洲国产色片| 亚洲经典国产精华液单| 有码 亚洲区| 亚洲欧美清纯卡通| 蜜桃久久精品国产亚洲av| 99久久九九国产精品国产免费| 三级国产精品欧美在线观看| 欧美日韩一区二区视频在线观看视频在线 | 麻豆乱淫一区二区| 欧美高清性xxxxhd video| 伦精品一区二区三区| 亚洲精品乱码久久久久久按摩| 国产精品一及| 天堂√8在线中文| 99在线视频只有这里精品首页| 久久久久久久久久久丰满| 熟妇人妻久久中文字幕3abv| 国产精品综合久久久久久久免费| 亚洲一级一片aⅴ在线观看| 18禁在线无遮挡免费观看视频| 久久草成人影院| 少妇的逼水好多| 一级黄色大片毛片| 男女做爰动态图高潮gif福利片| 一区二区三区四区激情视频 | 亚洲精品久久久久久婷婷小说 | 久久人人爽人人爽人人片va| 免费黄网站久久成人精品| 身体一侧抽搐| 日日摸夜夜添夜夜爱| 黄色欧美视频在线观看| 日韩精品有码人妻一区| 日本黄色片子视频| 亚洲精品自拍成人| 久久午夜亚洲精品久久| 欧美区成人在线视频| 国产视频首页在线观看| 久久人妻av系列| 自拍偷自拍亚洲精品老妇| 国内精品宾馆在线| 欧美日韩综合久久久久久| 我要看日韩黄色一级片| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 国产一区亚洲一区在线观看| 亚洲精品国产成人久久av| 中文在线观看免费www的网站| 国产精品一区二区在线观看99 | 91久久精品国产一区二区三区| 天天一区二区日本电影三级| 日韩,欧美,国产一区二区三区 | 国产黄片美女视频| 波多野结衣高清无吗| 国产高清视频在线观看网站| 插阴视频在线观看视频| 精品日产1卡2卡| 免费人成在线观看视频色| 两个人的视频大全免费| 久久久久久久午夜电影| 青青草视频在线视频观看| 狠狠狠狠99中文字幕| 老司机福利观看| 亚洲av中文字字幕乱码综合| 日本-黄色视频高清免费观看| 丰满人妻一区二区三区视频av| 欧美成人a在线观看| 联通29元200g的流量卡| 国产亚洲5aaaaa淫片| 久久久久久久久久黄片| 日韩亚洲欧美综合| 99久国产av精品| 久久鲁丝午夜福利片| 日韩制服骚丝袜av| 久久国内精品自在自线图片| 欧美高清性xxxxhd video| 免费av不卡在线播放| 精品一区二区免费观看| 51国产日韩欧美| 亚洲色图av天堂| 国产av在哪里看| 长腿黑丝高跟| 美女高潮的动态| 少妇高潮的动态图| 99久久中文字幕三级久久日本| 免费人成视频x8x8入口观看| 亚洲图色成人| 日韩国内少妇激情av| 国产一区二区在线av高清观看| 99riav亚洲国产免费| 美女高潮的动态| 亚洲人成网站在线观看播放| 日韩精品有码人妻一区| av女优亚洲男人天堂| 搞女人的毛片| 免费一级毛片在线播放高清视频| 午夜福利在线在线| 国产高清有码在线观看视频| 精品欧美国产一区二区三| 国产乱人偷精品视频| 国产精品不卡视频一区二区| 中文亚洲av片在线观看爽| 亚洲18禁久久av| 熟女人妻精品中文字幕| 精品久久久久久久久久免费视频| 亚洲最大成人中文| 国产精品,欧美在线| 一区二区三区免费毛片| 观看免费一级毛片| 天堂中文最新版在线下载 | 男女下面进入的视频免费午夜| 中国美女看黄片| 亚洲av成人精品一区久久| 亚洲美女视频黄频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品久久久久久婷婷小说 | 如何舔出高潮| 亚洲自拍偷在线| 免费在线观看成人毛片| 男人舔奶头视频| 男人和女人高潮做爰伦理| 岛国毛片在线播放| 亚洲国产精品成人综合色| 国产伦理片在线播放av一区 | 男女那种视频在线观看| 99国产极品粉嫩在线观看| 欧美性感艳星| 一本久久精品| 午夜精品在线福利| 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 又爽又黄a免费视频| 99久久精品一区二区三区| 亚洲经典国产精华液单| 中文精品一卡2卡3卡4更新| 亚洲国产色片| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 欧美zozozo另类| 久久久国产成人精品二区| 白带黄色成豆腐渣| 我的女老师完整版在线观看| 男插女下体视频免费在线播放| 国产精品电影一区二区三区| 亚洲欧美日韩无卡精品| 一个人免费在线观看电影| 一本精品99久久精品77| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 99视频精品全部免费 在线| 国产高清视频在线观看网站| 少妇熟女欧美另类| 人人妻人人看人人澡| 99热这里只有精品一区| 亚洲经典国产精华液单| 97在线视频观看| 十八禁国产超污无遮挡网站| 日韩精品有码人妻一区| 国产69精品久久久久777片| 麻豆久久精品国产亚洲av| 毛片一级片免费看久久久久| av在线老鸭窝| 久久99蜜桃精品久久| 最近中文字幕高清免费大全6| 直男gayav资源| 一区二区三区高清视频在线| 丰满乱子伦码专区| 如何舔出高潮| 边亲边吃奶的免费视频| 免费看美女性在线毛片视频| 人妻久久中文字幕网| 欧美日本视频| 99热网站在线观看| 在线免费观看不下载黄p国产| 亚洲精华国产精华液的使用体验 | 国产精品美女特级片免费视频播放器| 亚洲国产精品国产精品| 91aial.com中文字幕在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 精品国产三级普通话版| 亚洲精品久久国产高清桃花| 免费人成在线观看视频色| 国产一区二区三区av在线 | 熟女电影av网| 亚洲最大成人av| 国产日韩欧美在线精品| 只有这里有精品99| 午夜免费激情av| 波野结衣二区三区在线| 最近手机中文字幕大全| 一进一出抽搐gif免费好疼| 国产精品美女特级片免费视频播放器| 亚洲真实伦在线观看| 干丝袜人妻中文字幕| 日本与韩国留学比较| 非洲黑人性xxxx精品又粗又长| 日韩人妻高清精品专区| 少妇熟女欧美另类| av国产免费在线观看| 欧美3d第一页| 搡女人真爽免费视频火全软件| 99热这里只有精品一区| 我的女老师完整版在线观看| 欧美最黄视频在线播放免费| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 国产黄片视频在线免费观看| 亚洲色图av天堂| 能在线免费看毛片的网站| 欧美三级亚洲精品| 51国产日韩欧美| 日本与韩国留学比较| 国产69精品久久久久777片| 精华霜和精华液先用哪个| 国产精品久久电影中文字幕| 99久久精品热视频| 成人综合一区亚洲| 丰满乱子伦码专区| 不卡视频在线观看欧美| 99久久无色码亚洲精品果冻| 精品熟女少妇av免费看| 亚洲成人久久爱视频| 亚洲真实伦在线观看| 国产精品久久久久久久久免| 美女国产视频在线观看| 成人美女网站在线观看视频| 久久久久久国产a免费观看| 美女被艹到高潮喷水动态| 中文字幕熟女人妻在线| 国产成人freesex在线| 国产成人a∨麻豆精品| 亚洲中文字幕日韩| 狂野欧美白嫩少妇大欣赏| 欧美日本亚洲视频在线播放| 亚洲美女搞黄在线观看| 丰满的人妻完整版| 日韩大尺度精品在线看网址| 国产精品精品国产色婷婷| 人妻制服诱惑在线中文字幕| 一本久久中文字幕| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区| 身体一侧抽搐| 国产片特级美女逼逼视频| 中文字幕av成人在线电影| 啦啦啦啦在线视频资源| 久久这里只有精品中国| 啦啦啦观看免费观看视频高清| 国产一区二区激情短视频| 尾随美女入室| 丝袜美腿在线中文| 成人特级黄色片久久久久久久| 亚洲欧美日韩高清专用| 久久久午夜欧美精品| 婷婷色av中文字幕| 日日撸夜夜添| 美女被艹到高潮喷水动态| 久久综合国产亚洲精品| 在线观看66精品国产| 亚洲四区av| 神马国产精品三级电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| 日韩一区二区视频免费看| 在线观看美女被高潮喷水网站| 一个人看的www免费观看视频| 在现免费观看毛片| 久久人人精品亚洲av| 国产久久久一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | www日本黄色视频网| 亚洲国产精品国产精品| 成年女人永久免费观看视频| 国产精品99久久久久久久久| 99riav亚洲国产免费| 久久热精品热| 校园春色视频在线观看| 又粗又硬又长又爽又黄的视频 | 精品久久国产蜜桃| 国内精品一区二区在线观看| 午夜精品在线福利| 91aial.com中文字幕在线观看| 国内久久婷婷六月综合欲色啪| 一边亲一边摸免费视频| 国产精品99久久久久久久久| 国产成人影院久久av| 久久精品人妻少妇| 99热这里只有精品一区| 亚洲精品粉嫩美女一区| 亚洲在久久综合| 久久久久九九精品影院| av在线观看视频网站免费| 欧美另类亚洲清纯唯美| 九色成人免费人妻av| 日韩国内少妇激情av| 日本一二三区视频观看| 国产精品福利在线免费观看| 乱码一卡2卡4卡精品| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 亚洲精品国产av成人精品| 久久欧美精品欧美久久欧美| 亚洲图色成人| 久久久久久大精品| av福利片在线观看| 中文字幕精品亚洲无线码一区| 日韩欧美精品v在线| 久久人妻av系列| 中国国产av一级| 成人亚洲欧美一区二区av| 狂野欧美白嫩少妇大欣赏| 国内揄拍国产精品人妻在线| 日韩一区二区视频免费看| 欧美激情在线99| 两个人的视频大全免费| 特大巨黑吊av在线直播| АⅤ资源中文在线天堂| 亚洲成av人片在线播放无| 国产成人a∨麻豆精品| 欧美一区二区精品小视频在线| 特大巨黑吊av在线直播| 变态另类丝袜制服| 老熟妇乱子伦视频在线观看| 国产成人aa在线观看| 久久精品夜夜夜夜夜久久蜜豆| 日韩视频在线欧美| 久久6这里有精品| 日韩欧美 国产精品| 又爽又黄无遮挡网站| 深夜精品福利| 国产午夜精品论理片| 亚洲欧美清纯卡通| 日韩欧美精品v在线| 亚洲在线自拍视频| 中文字幕熟女人妻在线| 插逼视频在线观看| 啦啦啦啦在线视频资源| 亚洲国产欧美在线一区| 日韩,欧美,国产一区二区三区 | 精品久久久久久久人妻蜜臀av| 国产成人a∨麻豆精品| 亚洲精品亚洲一区二区| 日韩欧美国产在线观看| 久久综合国产亚洲精品| 日本黄大片高清| 欧美3d第一页| 长腿黑丝高跟| 成年版毛片免费区| 99热6这里只有精品| 亚洲av男天堂| 国产伦在线观看视频一区| 国内精品宾馆在线| 国产一区二区在线观看日韩| 午夜精品一区二区三区免费看| 两个人视频免费观看高清| 久久久久久九九精品二区国产| 亚洲不卡免费看| 国产亚洲91精品色在线| 成年版毛片免费区| 秋霞在线观看毛片| 男女做爰动态图高潮gif福利片| 日韩一区二区视频免费看| 18禁在线无遮挡免费观看视频| 99国产极品粉嫩在线观看| 男女那种视频在线观看| 一区二区三区四区激情视频 | 精品日产1卡2卡| 亚洲欧洲日产国产| 精品少妇黑人巨大在线播放 | 看片在线看免费视频| 欧美成人一区二区免费高清观看| 又粗又爽又猛毛片免费看| 国产久久久一区二区三区| 最近视频中文字幕2019在线8| 免费观看在线日韩| 成人鲁丝片一二三区免费| 久久精品影院6| 一个人看的www免费观看视频| 日本撒尿小便嘘嘘汇集6| 国产一区亚洲一区在线观看| 国产毛片a区久久久久| 男人和女人高潮做爰伦理| 国内揄拍国产精品人妻在线| 欧美高清成人免费视频www| 99久久精品国产国产毛片| 少妇被粗大猛烈的视频| 欧美最黄视频在线播放免费| 亚洲婷婷狠狠爱综合网| 此物有八面人人有两片| 国产在视频线在精品| 嫩草影院新地址| 国产午夜福利久久久久久| 中文字幕免费在线视频6| 久久久精品欧美日韩精品| 国内精品美女久久久久久| 不卡视频在线观看欧美| 午夜爱爱视频在线播放| 亚洲av男天堂| 久久鲁丝午夜福利片| 精品无人区乱码1区二区| av又黄又爽大尺度在线免费看 | 亚洲中文字幕一区二区三区有码在线看| 欧美极品一区二区三区四区| 激情 狠狠 欧美| 麻豆久久精品国产亚洲av| 久久久久久久久久久丰满|