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

    煤巖體振動破壞試驗及微震信號特征

    2015-02-21 01:50:35李成武孫曉元高天寶解北京徐曉萌
    煤炭學(xué)報 2015年8期

    李成武,孫曉元,2,高天寶,解北京,徐曉萌,4

    (1.中國礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院,北京 100083;2.太原科技大學(xué)環(huán)境與安全學(xué)院,山西太原 030024;3.天津城建大學(xué)能源與安全工程學(xué)院,天津 300384;4.科廷大學(xué)石油工程系,西澳大利亞州珀斯 6151)

    ?

    煤巖體振動破壞試驗及微震信號特征

    李成武1,孫曉元1,2,高天寶3,解北京1,徐曉萌1,4

    (1.中國礦業(yè)大學(xué)(北京)資源與安全工程學(xué)院,北京 100083;2.太原科技大學(xué)環(huán)境與安全學(xué)院,山西太原 030024;3.天津城建大學(xué)能源與安全工程學(xué)院,天津 300384;4.科廷大學(xué)石油工程系,西澳大利亞州珀斯 6151)

    摘 要:通過所搭建的大型振動試驗裝置和微震監(jiān)測系統(tǒng),探究不同激勵加速度和頻率條件下整體煤巖試件的振動破壞特性。利用希爾伯特-黃變換(HHT)對煤巖體受迫振動過程中的監(jiān)測信號進(jìn)行了模態(tài)分離和信號重構(gòu),并結(jié)合振動力學(xué)知識著重闡釋了共振放大效應(yīng)產(chǎn)生的機(jī)理。研究結(jié)果表明:微震傳感器監(jiān)測信號具有明顯的載波特征;當(dāng)激勵加速度峰值達(dá)到一定水平(0.5g)時,煤巖試件出現(xiàn)有效微震信號,首次破裂發(fā)生在煤層中,信號隨傳播距離的增大而衰減;隨著裂隙的產(chǎn)生和擴(kuò)展,煤巖試件自振頻率降低,當(dāng)激勵頻率等于煤體自振頻率時,試件會產(chǎn)生共振放大效應(yīng),微震信號明顯增強;激勵加速度和激勵頻率對微震信號均有影響,激勵加速度峰值越大、激勵頻率越接近自振頻率,則微震信號越密集、幅值越大,且出現(xiàn)了優(yōu)勢頻段向低頻轉(zhuǎn)移的趨勢;理論分析解釋了上述現(xiàn)象產(chǎn)生的原因。

    關(guān)鍵詞:振動破壞;激振響應(yīng);共振放大效應(yīng);微震監(jiān)測;HHT

    責(zé)任編輯:畢永華

    李成武,孫曉元,高天寶,等.煤巖體振動破壞試驗及微震信號特征[J].煤炭學(xué)報,2015,40(8):1834-1844.doi:10.13225/ j.cnki.jccs.2014.1153

    煤巖破壞的發(fā)生機(jī)制包括準(zhǔn)靜載作用下的流變破壞和動載作用下的瞬態(tài)破壞[1],不同試驗條件下煤巖體所表現(xiàn)出的破壞特性和規(guī)律有所差異。國內(nèi)外學(xué)者深入分析了單軸壓縮、三軸壓縮和拉伸剪切等室內(nèi)準(zhǔn)靜載條件下煤巖體破壞的物理力學(xué)性質(zhì)[2]。當(dāng)前,以分離式霍普金森壓桿(SHPB)為主要裝置的沖擊動載破壞試驗也成為研究的熱點[3-6],SHPB試驗可以較好地闡釋煤巖體在掘進(jìn)爆破、石門揭煤等承受強烈沖擊載荷時的動力學(xué)特性[7]。然而,除直接沖擊破壞外,工作面來壓、頂板斷裂、機(jī)械振動及相鄰區(qū)域的干擾等動力擾動對煤巖失穩(wěn)破壞的影響同樣不容小視。這些廣泛存在的擾動一方面使得煤巖體及周邊應(yīng)力急劇變化,形成局部應(yīng)力集中,另一方面擾動產(chǎn)生的振動波在煤巖體中傳播時可促進(jìn)裂紋的產(chǎn)生與擴(kuò)展[8],使煤巖體內(nèi)摩擦角減小,抗剪強度降低[9],進(jìn)而造成煤巖體失穩(wěn)破壞并最終導(dǎo)致煤巖動力災(zāi)害的發(fā)生。

    國內(nèi)外學(xué)者從理論分析和數(shù)值模擬角度進(jìn)行了大量動力擾動誘發(fā)煤巖動力災(zāi)害的探討[10-11]。就試驗研究而言,20世紀(jì)90年代潘一山等[12]曾利用改裝后的萬能試驗機(jī)進(jìn)行局部煤體振動試驗,并結(jié)合分形理論揭示了初應(yīng)力水平和振動次數(shù)對煤體裂隙分布和力學(xué)性能的影響。然而試驗中并未涉及振動頻率和振動加速度等因素對整體煤巖試件破壞的影響。研究表明,煤巖破壞主要與其內(nèi)部裂紋的產(chǎn)生和擴(kuò)展有關(guān)[13],而裂紋演變時往往伴隨微震信號的產(chǎn)生[14]。因此,筆者引入用于地震作用下建筑、巖體與邊坡穩(wěn)定性分析的大型振動系統(tǒng),以微震監(jiān)測為手段,分析整體煤巖試件在不同激勵頻率和加速度條件下的動力學(xué)特性,著重探討從微震信號特征角度反映煤巖體破壞程度的可行性,并給出理論解釋。研究成果將有助于完善煤巖動力災(zāi)害的發(fā)生機(jī)理,提高微震預(yù)測預(yù)警技術(shù)的準(zhǔn)確率。

    1 試驗系統(tǒng)與參數(shù)設(shè)置

    1.1 受迫振動系統(tǒng)與試件制備

    系統(tǒng)在外界控制的持續(xù)激勵作用下產(chǎn)生的振動稱為受迫振動。為達(dá)到整體煤巖試件受迫振動時激勵參數(shù)(頻率和加速度)均可精確調(diào)制的目的,采用臺面尺寸為3.0 m×3.0 m的大型振動臺來進(jìn)行試驗。該振動臺最大有效載重質(zhì)量為10 t,激勵輸入頻率范圍為0.2~50 Hz,滿載時最大激勵加速度為1.2g。振動過程中產(chǎn)生的最大偏心力矩100 kN·m,最大傾覆力矩300 kN·m,振動方式為一維水平往復(fù)振動。

    圖1 一維水平往復(fù)振動臺Fig.1 Single-dimensional horizontal and reciprocating vibrating table

    由于受迫振動系統(tǒng)中模型箱的容積有限,且其會跟隨振動臺運動,故無法準(zhǔn)確反映涉及半無限自由場的煤巖體實際條件。因此,采用能較好消除邊界影響的疊層剪切模型箱[15]。該模型箱由10層獨立的方型鋼制框架疊合而成,每層框架為槽型鋼焊接,斷面尺寸為60 mm×30 mm×3 mm(高×寬×厚),所圍成的框架內(nèi)部尺寸為400 mm×400 mm×60 mm。框架間隙6 mm,其中設(shè)置4個滾軸軸承,以減小層間摩擦和各層框架的撓度[16]。在模型箱自下而上的第2,5和9層的框架上每邊分別開3個45 mm×45 mm的方孔,以便布設(shè)微震傳感器和觀察煤巖試件表面的裂紋演變情況。

    圖2 制作完成后的模型試件Fig.2 Model specimen after completion

    為降低真實煤巖體非均質(zhì)、各向異性及原生結(jié)構(gòu)面對微震信號采集的影響,煤巖試件選用中砂和煤粉作為骨料,并輔以水泥作為膠結(jié)材料[17]。通過相似比論證,選擇幾何相似常數(shù)為50,底板厚度設(shè)為70 mm,煤層厚度為80 mm,頂板厚度為500 mm。選擇容重相似常數(shù)為1.4,由此確定的強度相似比為70。為區(qū)分煤層和巖層的力學(xué)性質(zhì),煤體配比時適當(dāng)降低水泥的比例,并摻入部分石灰材料。按照上述參數(shù),經(jīng)多次試配后所獲得的理想材料配合比見表1。在層間均勻鋪設(shè)滑石粉,以模擬振動時產(chǎn)生的層面相對運動。

    表1 材料配合比Table 1 Mix proportion material

    1.2 微震傳感器與ZDKT-1型微震監(jiān)測系統(tǒng)

    微震傳感器采用COLIBRYS公司生產(chǎn)的SF1500MEMS電容性加速度計。該傳感器具有動態(tài)范圍高、變形率低、耐沖擊、穩(wěn)定性好的優(yōu)點,尤其適用于動力擾動作用下的波形監(jiān)測,其頻響范圍為0~1 500 Hz,對弱信號可提高5 kHz,動態(tài)范圍為120 dB,正常工作溫度為-40~+125℃,可受沖擊限制為1 500 g。通過預(yù)先進(jìn)行的微振動信號測試實驗,發(fā)現(xiàn)傳感器靈敏度和準(zhǔn)確性很高,能較好地滿足試驗要求。

    微震監(jiān)測系統(tǒng)為中國礦業(yè)大學(xué)(北京)自主研發(fā)的ZDKT-1型煤巖動力災(zāi)害實驗?zāi)M系統(tǒng)。該系統(tǒng)集成了微震信號采集芯片,共12路通道,最大采樣頻率可達(dá)51 kHz,試驗設(shè)置的采集頻率為17 067 Hz。利用配套的MKJC煤巖動力監(jiān)測預(yù)警系統(tǒng)采集軟件,可實現(xiàn)采集數(shù)據(jù)的多通道實時監(jiān)測,監(jiān)測系統(tǒng)及數(shù)據(jù)采集軟件如圖3所示。

    圖3 微震監(jiān)測系統(tǒng)與數(shù)據(jù)采集軟件Fig.3 Micro-seismic monitoring system and data acquisition software

    圖4(a)為煤巖試件振動破壞試驗系統(tǒng),試驗中共布設(shè)了7個微震傳感器,利用改性丙烯酸酯將其貼于層疊箱框架方孔處的煤巖體表面上。微震傳感器Z1和Z4貼于煤層表面,Z3和Z5,Z7和Z6分別布置于不同高度的上覆巖層之上,其中Z1,Z3和Z7沿振動方向布置,Z4,Z5和Z6則垂直于振動方向布置。為分析煤巖試件振動過程中豎向方向上的微震信號變化特征,在其頂部布置微震傳感器Z2。

    圖4 試驗系統(tǒng)及微震傳感器布置Fig.4 Sketch diagram of experimental system and arrangement plan of micro-seismic sensors

    1.3 試驗參數(shù)設(shè)置

    為研究不同激勵頻率和加速度對煤巖試件破壞的影響,需要確定合理的試驗參數(shù)。劉保縣等指出,擾動頻率與煤巖體自振頻率之比小于為阻尼比)為振動誘發(fā)煤巖動力災(zāi)害的必要條件[18],故在確定試驗參數(shù)之前,需首先測定煤巖體試件的自振頻率。

    試驗選用錘擊法測定煤巖體試件的初始自振頻率,方法為使用脈沖力錘敲擊試件受測點,如果力錘敲擊時間夠短,則試件實際上受到一個沖擊脈沖的作用,并在該脈沖力的作用下以其自身固有的低階自振頻率做自由振動[19]。通過固定在試件頂部的速度傳感器接收振動信號,經(jīng)電荷電壓濾波放大器加工后,對傳感器數(shù)據(jù)進(jìn)行處理可得出試件的自振頻率。煤巖試件固有頻率測試系統(tǒng)如圖5所示。

    圖5 試件自振頻率測試系統(tǒng)Fig.5 Natural frequency test system

    試驗中自振頻率測試系統(tǒng)的采樣頻率設(shè)置為1 300 Hz,對速度傳感器采集的信號進(jìn)行快速傅里葉變換(FFT)。限于篇幅,以2次測試實驗為例,其振動信號及頻譜曲線如圖6所示。

    圖6 錘擊法測定的自由振動曲線與頻譜Fig.6 Free vibration curve and spectrogram by hammering method

    從圖6(b),(d)的振動測試頻譜圖中可看出,兩次捶擊法測得的自振頻率分別為47.1和50.2 Hz。不同煤巖介質(zhì)的阻尼比相差較大,文獻(xiàn)[20]通過試驗得出其介于2.64%~5.34%之間。由此計算可知振動臺的頻率應(yīng)以上述測試參數(shù)為基準(zhǔn)。

    振動臺啟動后,其激勵設(shè)置為交替輸入白噪聲掃描和多種正弦波信號。其中白噪聲掃描是指涵蓋所有頻率、各頻率均具有相同能量且功率譜密度在整個頻域內(nèi)均勻分布的隨機(jī)噪聲,目的在于實時測定煤巖試件的自振頻率,白噪聲掃描的加速度設(shè)置為0.02g。正弦波信號按0.02g,0.1g,0.5g,0.75g, 1.5g和2.4g六組逐級提高的峰值加速度等級依次輸入,每組加速度又分別對應(yīng)著不同的激勵頻率。整個試驗共設(shè)置了83個工況點,具體參數(shù)設(shè)置見表2。

    表2 試驗參數(shù)設(shè)置Table 2 Setting of experimental parameters

    續(xù)表

    2 試驗結(jié)果與分析

    2.1 未破裂時監(jiān)測信號特征

    當(dāng)激振加速度較小時,微震傳感器并未監(jiān)測到有效的煤巖破裂信號。以工況點27為例,此時激勵頻率為34 Hz,加速度峰值0.1g,振動周期數(shù)120個。圖7(a)是該工況點時部分傳感器所監(jiān)測到的原始信號。

    圖7 工況點27時所監(jiān)測到的部分傳感器信號及處理結(jié)果Fig.7 Part of the sensor signals and results process of working condition 27

    為記錄振動全過程,微震監(jiān)測系統(tǒng)先于振動臺啟動,工況點27中提前時間為2.8 s。由圖7(a)可知,傳感器的監(jiān)測數(shù)據(jù)為環(huán)境噪聲和試件振動等多種信號疊加的結(jié)果,這一特性在沿振動方向布置的Z1和Z7監(jiān)測數(shù)據(jù)中體現(xiàn)的尤為顯著。對于Z4,盡管其布置方向與振動方向垂直,理論上不應(yīng)包含振動響應(yīng),但在試驗中完全垂直是無法做到的,故Z4監(jiān)測數(shù)據(jù)中也不可避免的包含了振動信號。此外,對比布置煤巖試件底部Z1和上部Z7的傳感器監(jiān)測幅值,可知試件上部對激勵的振動響應(yīng)更為明顯。

    由于微震監(jiān)測數(shù)據(jù)為典型的非線性信號,故可采用希爾伯特-黃變換(HHT)對其進(jìn)行處理。圖7(b)為對Z7傳感器進(jìn)行整體經(jīng)驗?zāi)B(tài)分解(EEMD)的分析結(jié)果,信號被分解為16個固有模態(tài)函數(shù)(IMF)分量和一個殘余分量R,各IMF分量按頻率從高到低依次排列[21],其中IMF5和IMF6為典型的振動響應(yīng)信號,其能量占比分別達(dá)到了94.3%和5.5%。將上述振動響應(yīng)信號剔除后的重構(gòu)信號如圖7(c)所示,可見信號非常密集且幅值較小,沒有體現(xiàn)出典型微震信號所具備的特征[22],由此判斷該信號為環(huán)境與EEMD計算過程中所附帶的白噪聲,說明該工況點下煤巖試件未出現(xiàn)裂紋演變。

    2.2 煤巖試件破裂時的微震信號特征

    在工況29時,振動激勵加速度峰值達(dá)到0.5g,頻率為2 Hz,振動周期30個?,F(xiàn)對該工況下的監(jiān)測信號進(jìn)行詳細(xì)分析。

    圖8(a)為工況點29時所監(jiān)測到的部分傳感器信號,與圖7(a)相比,傳感器Z1和Z4監(jiān)測信號的幅值大為增加,且第1次出現(xiàn)了較為典型的微震信號特征,這說明煤巖試件在此工況點時首次產(chǎn)生裂紋演變。對Z1監(jiān)測信號進(jìn)行EEMD分解,可知該信號也同樣包含了具有較高頻率的煤巖破裂微震信號和較低頻率的試件振動信號,具有典型的載波特征。在圖8(b)中,IMF9為典型的正弦波信號,其能量占比達(dá)到了95.9% ,結(jié)合各IMF的快速傅里葉變換(FFT)和能量占比的計算結(jié)果,確定圖8(b)中具有典型的微震信號特征且能量占比相對較大的IMF1~IMF4為有效信號。利用上述分量進(jìn)行微震信號重構(gòu)。類似地,對布置于煤層表面的Z4傳感器與上覆巖層表面Z3, Z7傳感器的監(jiān)測數(shù)據(jù)也進(jìn)行重構(gòu),其有效微震信號如圖8(c)所示。由于煤巖試件的體積相對較小,各傳感器接收到微震信號的時刻未能發(fā)現(xiàn)明顯差異。通過對比各傳感器的監(jiān)測幅值,可以看出Z1,Z4監(jiān)測到的信號最大,位于巖層中部的Z3次之,貼于巖層上部的傳感器Z7最小,由此可以判定試件的首次破裂發(fā)生于煤層之中,且微震信號隨傳播距離的增加而逐漸衰減。

    圖8(d)為利用HHT計算得出的Z1傳感器重構(gòu)信號的EEMD分解、各IMF分量頻譜、Hilbert能量譜和邊際譜分析結(jié)果。通過EEMD分解可知IMF1~IMF3為微震信號的主體部分,其能量占比分別為69.5%,21.9%和7.56%,三者之和達(dá)到了有效微震信號總能量的98.96%。從各IMF分量的FFT計算結(jié)果可以看出,IMF1的頻帶范圍較寬,從500 Hz到4 500 Hz均有分布,且主要集中在1 000~1 500 Hz部分,峰值位于1 360 Hz;IMF2和IMF3的頻率范圍分別為200~1 500 Hz和200~1 000 Hz。由于上述3個IMF分量為微震信號主體部分,故可認(rèn)為其整體優(yōu)勢頻段為200~4 500 Hz。Hilbert能量譜簡潔明了地反映微震信號時間、頻率和能量幅值之間的關(guān)系,圖中絕大部分微震信號在4 000 Hz以下,且主要集中于500 Hz左右。Hilbert邊際譜的分析結(jié)果也驗證了這一點,圖8(d)中的有效微震信號集中在200~4 500 Hz區(qū)間內(nèi),主頻分別為500,1 180和2 800 Hz。

    類似地,對其他傳感器的有效微震信號進(jìn)行HHT分析,發(fā)現(xiàn)處理結(jié)果與Z1基本相同,故可認(rèn)為在工況點29(加速度0.5g、頻率2 Hz)時,煤巖體首次發(fā)生了破裂,裂隙產(chǎn)生于煤層位置,各傳感器的監(jiān)測信號均具有載波特性,重構(gòu)信號的優(yōu)勢頻段分布與Z1相近。

    2.3 煤巖試件共振時的微震信號特征

    在工況點29之后,振動臺輸入的激勵加速度峰值均設(shè)置在0.5g及以上,且每個工況點都能監(jiān)測到有效微震信號,說明此階段煤巖試件內(nèi)發(fā)生著裂隙演化,微震信號特征也與工況點29基本類似。但在工況點37(振動激勵所輸入的加速度峰值達(dá)到0.5g,振動頻率20 Hz,振動周期數(shù)120個)時,傳感器監(jiān)測信號出現(xiàn)明顯不同。

    對比圖8工況點29和圖9工況點37重構(gòu)的有效微震信號,可以看出后者極為密集。對Z1有效微震信號進(jìn)行HHT分析,發(fā)現(xiàn)能量占比最大的并非IMF1,而是相對頻率較低的IMF3(49.4%)。通過能量譜和邊際譜也可以看出,其優(yōu)勢頻段位于200~1 400 Hz,峰值為300和1 250 Hz,顯然此頻率要比工況29時低得多。此外,由Hilbert能量譜圖可知,其優(yōu)勢IMF分量的能量大于工況點29。

    在工況點37時各傳感器的信號處理結(jié)果與Z1基本相同,試驗過程中也發(fā)現(xiàn)煤巖試件的響應(yīng)異常強烈,這說明此工況點條件下煤巖試件出現(xiàn)了共振放大效應(yīng)。在此階段裂隙演化明顯加劇,對應(yīng)著微震信號的密集出現(xiàn)且振幅開始增加,微震信號的優(yōu)勢頻段向低頻段移動。根據(jù)陸菜平、竇林名等在文獻(xiàn)[23]中所總結(jié)的現(xiàn)場研究成果,這是煤巖結(jié)構(gòu)破壞和動力災(zāi)害發(fā)生的重要前兆。

    圖8 工況點29時所監(jiān)測到的部分傳感器信號及處理結(jié)果Fig.8 Part of the sensor signals and results process of working condition 29

    圖9 工況點37部分有效微震信號及處理結(jié)果Fig.9 Part of the effective micro-seismic signals and results in process of working condition 37

    需要說明的是,共振放大效應(yīng)在試驗中發(fā)生了多次。事實上,隨著裂紋的不斷產(chǎn)生與擴(kuò)展,煤巖試件的自振頻率呈逐步下降趨勢。這也從側(cè)面闡釋了煤礦開采在經(jīng)歷多次擾動之后,某些較低頻率的擾動源也會導(dǎo)致煤巖結(jié)構(gòu)破壞的原因。

    在工況點83(激勵加速度峰值2.4g,頻率5 Hz) 時,試件產(chǎn)生了猛烈搖晃,并由于裂隙的快速開展和閉合而向外界噴出大量的煤巖灰塵和煙霧,分析認(rèn)為此時煤巖試件已經(jīng)遭到嚴(yán)重破壞。停止試驗并拆除模型箱,發(fā)現(xiàn)上覆巖層基本完好,未出現(xiàn)較大的裂隙,而煤層產(chǎn)生了大尺度的宏觀裂縫,煤體內(nèi)部十分松散,并已完全破壞(圖10)。這說明振動破壞易發(fā)生于煤巖試件的薄弱位置,試驗中各傳感器監(jiān)測到的微震信號主要是煤層破裂產(chǎn)生的。

    圖10 煤體破裂Fig.10 Coal burst figure

    圖11 不同激勵頻率時的有效微震信號(加速度0.5g)Fig.11 Effective micro-seismic signals of different vibration frequencies (Acceleration is 0.5g)

    2.4 激勵加速度和激勵頻率對微震信號影響

    圖11和圖12分別描述了相同激勵加速度不同頻率(工況點29,35,37,41和47)和相同頻率不同加速度(工況點43,53和63)時傳感器Z1所監(jiān)測到的有效微震信號。圖11中各工況點的激勵加速度均為0.5g,激勵頻率分別為2,15,20,27和34 Hz。通過上節(jié)的分析已知煤巖試件在激勵頻率為20 Hz時發(fā)生了共振放大現(xiàn)象,并可以看出在此之前監(jiān)測到的有效微震信號較少,共振后試件所產(chǎn)生和擴(kuò)展的裂紋明顯增多,表明此時煤巖體已經(jīng)進(jìn)入非線性動力學(xué)階段[24]。結(jié)合各信號的HHT分析結(jié)果可知,激勵頻率對煤巖體的裂隙演化和結(jié)構(gòu)破壞具有十分顯著的影響,其越接近煤體的自振頻率,煤巖體產(chǎn)生較大幅值的微震信號越多,優(yōu)勢頻段向低頻轉(zhuǎn)移的趨勢越明顯。反之激勵頻率越遠(yuǎn)離自振頻率,越有利于保持煤巖結(jié)構(gòu)的穩(wěn)定性。

    圖12顯示的是激勵頻率均為30 Hz,加速度峰值為0.5g,0.75g和1.5g時由Z1傳感器監(jiān)測數(shù)據(jù)重構(gòu)的有效微震信號。從圖12可以看出,隨激勵加速度的增大,煤巖試件所產(chǎn)生的微震信號明顯增多,信號幅值也逐漸增大,這反映出煤巖試件的裂隙演化程度不斷加劇。此外,對各有效信號進(jìn)行HHT分析,同樣發(fā)現(xiàn)了優(yōu)勢頻段向低頻轉(zhuǎn)移的現(xiàn)象。

    圖12 不同激勵加速度時的有效微震信號(頻率30 Hz)Fig.12 Effective micro-seismic signals of different vibrationaccelerations (frequency is 30 Hz)

    3 一維振動響應(yīng)理論

    圖13 一維振動條件下的變形與受力示意Fig.13 Deformation and force under the condition of singledimensional vibration

    煤巖振動破壞過程中的微震等信號表現(xiàn)出復(fù)雜的特征和現(xiàn)象[16],且受激勵頻率和激勵加速度雙重作用的影響,下面著重對上述動力響應(yīng)和微震信號特征進(jìn)行理論上的分析。由于試驗采用了一維水平往復(fù)振動臺,故整體煤巖試件的運動方式可視為一維振動條件下對激勵源的動力響應(yīng),簡化的變形及受力情況如圖13所示。

    設(shè)某時刻振動臺經(jīng)過位移xg達(dá)到振幅位置,由于層疊箱底部的框架與振動臺通過螺栓固定在一起,故煤巖試件底部的O點也隨之停止運動。此時其上部O′點處的質(zhì)點m仍在慣性作用下繼續(xù)移動。設(shè)質(zhì)點m經(jīng)歷時間t的位移為xg+x,則在這一過程中質(zhì)點m受到了慣性力fI、彈性恢復(fù)力fr及阻尼力fc的作用,即存在公式:

    其中,c為質(zhì)點m周圍煤巖介質(zhì)的阻尼系數(shù);k為其恢復(fù)平衡位置的抗側(cè)移剛度,根據(jù)力的平衡關(guān)系可將式(1)改寫為

    顯然,式(2)為二階非齊次線性微分方程。由振 動力學(xué)理論,令系統(tǒng)的固有頻率阻尼比將其代入式(2),得到

    針對方程(3),本文重點分析其在激勵為正弦波信號時的特解,設(shè)正弦信號的表達(dá)形式為

    其中,A為振動臺的振幅;ωg為輸入激勵的圓頻率。將式(4)代入式(3)中得到

    利用對應(yīng)系數(shù)法計算式(5)的特解為

    式(6)可簡化為圓頻率ωg為周期的運動,即

    其中,B為質(zhì)點m在振動臺激勵作用下受迫振動的振幅;φ為相位角,其表達(dá)式為

    則質(zhì)點振幅與振動臺振幅的比值即振幅放大系數(shù)為

    式(9)反映出β的計算結(jié)果主要取決于頻率比ωg/ω和阻尼比ξ。特別地,當(dāng)振動臺所輸入激勵的圓頻率ωg與煤體試件的固有角頻率ω相同時,振幅放大系數(shù)達(dá)到最大值。

    顯然,β一定時,振動臺幅值A(chǔ)越大,煤巖體的響應(yīng)幅值B也隨之增大。此外,從圖14可以看出,振動臺輸入激勵的圓頻率與煤體試件的固有角頻率相接近時,振幅放大系數(shù)迅速增大。當(dāng)ξ介于2.64%~5.34%之間時[20],其振幅放大系數(shù)為10~20倍,從而表現(xiàn)為較小的激勵即能使煤巖試件產(chǎn)生較大的動力響應(yīng)。即在振動臺加速度不變的情況下,調(diào)整激勵圓頻率ωg的大小,就可以達(dá)到煤巖試件共振的效果。此時其動力響應(yīng)將顯著增強,裂隙演變加劇,易于造成結(jié)構(gòu)的失穩(wěn)與破壞。

    圖14 一維簡諧受迫振動振幅放大系數(shù)曲線Fig.14 Magnification coefficient curve of single-dimensional harmonic forced vibration

    4 結(jié) 論

    (1)在一維水平往復(fù)振動試驗中,不同方向布置的微震傳感器監(jiān)測數(shù)據(jù)中均包含煤巖試件本身振動所產(chǎn)生的信號,監(jiān)測結(jié)果具有典型的載波特征,可借助HHT來進(jìn)行模態(tài)分解和信號重構(gòu)。

    (2)當(dāng)激勵加速度峰值達(dá)到一定水平(0.5g)時,煤巖試件產(chǎn)生有效微震信號,首次破裂發(fā)生在煤層中,信號隨傳播距離的增大而衰減。

    (3)當(dāng)激勵頻率等于煤體的自振頻率時,煤巖試件發(fā)生共振放大效應(yīng),裂隙演變加劇,微震信號增強,此后煤巖體處于非線性階段,且自振頻率也隨裂隙的演變而逐漸降低。

    (4)激勵加速度和激勵頻率對微震信號具有明顯的影響。激勵加速度越大、激勵頻率越接近煤體的自振頻率,微震信號越明顯、幅值越大,且出現(xiàn)優(yōu)勢頻段向低頻轉(zhuǎn)移的趨勢,這是煤巖破壞的重要前兆之一。

    (5)利用一維振動響應(yīng)理論闡釋了動力響應(yīng)現(xiàn)象和微震信號特征的產(chǎn)生原因,試驗現(xiàn)象與理論分析的結(jié)果基本一致。

    參考文獻(xiàn):

    [1]王 振,胡千庭,金洪偉.煤巖失穩(wěn)的發(fā)生機(jī)制研究[A].自主創(chuàng)新與持續(xù)增長第十一屆中國科協(xié)年會論文集(1)[C].中國科學(xué)技術(shù)協(xié)會,重慶市人民政府,2009:754-759.

    [2]謝和平,陳忠輝.巖石力學(xué)[M].北京:科學(xué)出版社,2004.

    [3]Hopkinson B.A method of measuring the pressure in the deformation of high explosives by the impact of bullets [J].Philos.Trans.R.Soc.A,1914,213:437-452.

    [4]黎立云,徐志強,謝和平,等.不同沖擊速度下巖石破壞能量規(guī)律的實驗研究[J].煤炭學(xué)報,2011,36(12):2007-2011.

    Li Liyun,Xu Zhiqiang,Xie Heping,et al.Failure experimental study on energy laws of rock under differential dynamic impact velocities [J].Journal of China Coal Society,2011,36(12):2007-2011.

    [5]金解放,李夕兵,殷志強,等.軸壓和循環(huán)沖擊次數(shù)對砂巖動態(tài)力學(xué)特性的影響[J].煤炭學(xué)報,2012,37(6):923-930.

    Jin Jiefang,Li Xibing,Yin Zhiqiang,et al.Effects of axial pressure and number of cyclic impacts on dynamic mechanical characteristics of sandstone[J].Journal of China Coal Society,2012,37(6):923-930.

    [6]李夕兵.巖石動力學(xué)基礎(chǔ)與應(yīng)用[M].北京:科學(xué)出版社,2014.[7]解北京.煤沖擊破壞動力學(xué)特性及磁場變化特征實驗研究[D].北京:中國礦業(yè)大學(xué)(北京),2013.

    [8]姜耀東,趙毅鑫,宋彥琦,等.放炮震動誘發(fā)煤礦巷道動力失穩(wěn)機(jī)理分析[J].巖石力學(xué)與工程學(xué)報,2005,24(17):3131-3136.

    Jiang Yaodong, Zhao Yixin, Song Yanqi, et al.Analysis of blasting tremor impact on roadway stability in coal mining[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(17):3131-3136.

    [9]呂秀江.煤巷掘進(jìn)影響區(qū)動態(tài)應(yīng)力響應(yīng)及對動力災(zāi)害影響研究[D].北京:中國礦業(yè)大學(xué)(北京),2014.

    [10]高天寶.煤體受迫振動響應(yīng)及破壞特征試驗研究[D].北京:中國礦業(yè)大學(xué)(北京),2013.

    [11]Xue Sheng,Yuan Liang,Wang Junfeng,et al.A coupled DEM and LBM model for simulation of outbursts of coal and gas[J].In-ternational Journal of Coal Science & Technology,2015,2(1): 22-29.

    [12]潘一山,杜廣林,張永利,等.煤體振動后力學(xué)性質(zhì)變化規(guī)律的試驗研究[J].巖土工程學(xué)報,1998,20(5):41-43.

    Pan Yishan,Du Guangling,Zhang Yongli,et al.An experimental study on the mechanical properties of coal mass after vibrating[J].Journal of Geotechnical Engineering,1998,20(5):41-43.

    [13]萬志軍,周楚良.巖石聲發(fā)射源機(jī)理的激振模型研究[J].礦山壓力與頂板管理,1998(4):68-70,80.

    Wan Zhijun,Zhou Chuliang.Study of shock generating model of uncovering mechanism of rock acoustic emission[J].Ground Pressure and Strata Control,1998(4):68-70,80.

    [14]姜福興,楊淑華.微地震監(jiān)測揭示的采場圍巖空間破裂形態(tài)[J].煤炭學(xué)報,2003,28(4):357-360.

    Jiang Fuxing,Yang Shuhua.Spatial fracturing progresses of surrounding rock masses in long wall face monitored by micro-seismic monitoring techniques [J].Journal of China Coal Society, 2003,28(4):357-360.

    [15]孫海峰,景立平,王寧偉,等.振動臺多功能疊層剪切箱研制[J].巖石力學(xué)與工程學(xué)報,2011,30(12):2498-2506.

    Sun Haifeng, Jing Linping, Wang Ningwei, et al.Development of multifunctional laminar shear container for shaking table test [J].Chinese Journal of Rock Mechanics and Engineering,2011, 30(12):2498-2506.

    [16]Li Chengwu,Hu Po,Gao Tianbao,et al.An experiment monitoring signals of coal bed simulation under forced vibration conditions [J].Shock and Vibration,doi:10.1155/2015/693612.

    [17]李曉紅,盧義玉,康 勇,等.巖石力學(xué)實驗?zāi)M技術(shù)[M].北京:科學(xué)出版社,2007.

    [18]劉??h,鮮學(xué)福,劉新榮,等.爆破激發(fā)煤瓦斯突出的研究[J].中國礦業(yè),2000(2):94-96.

    Liu Baoxian,Xian Xuefu,Liu Xinrong,et al.Coal gas outburst induced by blasting[J].China Mining Magazine,2000(2):94-96.

    [19]劉習(xí)軍,賈啟芬,張文德,等.工程振動與測試技術(shù)[M].天津:天津大學(xué)出版社,2004.

    [20]劉建鋒,徐 進(jìn),李青松,等.循環(huán)荷載下巖石阻尼參數(shù)測試的試驗研究[J].巖石力學(xué)與工程學(xué)報,2010,29(5):1036-1041.

    Liu Jianfeng,Xu Jin,Li Qingsong,et al.Study on damping parameters of rock under cyclic loading[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(5):1036-1041.

    [21]李成武,解北京,楊 威,等.基于HHT法的煤沖擊破壞SHPB測試信號去噪[J].煤炭學(xué)報,2012,37(11):1796-1801.

    Li Chengwu, Xie Beijing, Yang Wei, et al.Coal impact damage SHPB testing signal denoising based on HHT method[J].Journal of China Coal Society,2012,37(11):1796-1801.

    [22]李夕兵,凌同華,張義平.爆破震動信號分析理論與技術(shù)[M].北京:科學(xué)出版社,2009.

    [23]陸菜平,竇林名,吳興榮,等.煤巖沖擊前兆微震頻譜演變規(guī)律的試驗與實證研究[J].巖石力學(xué)與工程學(xué)報,2008,27(3): 519-525.

    Lu Caiping,Dou Linming,Wu Xingrong,et al.Experimental and empirical research on frequency-spectrum evolvement rule of rockburst precursory micro-seismic signals of coal-rock [J].Chinese Journal of Rock Mechanics and Engineering,2008,27(3):519-525.

    [24]陳 颙,黃庭芳,劉恩儒.巖石物理學(xué)[M].合肥:中國科學(xué)技術(shù)大學(xué)出版社,2009.

    Li Chengwu,Sun Xiaoyuan,Gao Tianbao,et al.Coal and rock vibration failure and the characteristics of micro-seismic signals[J].Journal of China Coal Society,2015,40(8):1834-1844.doi:10.13225/ j.cnki.jccs.2014.1153

    Coal and rock vibration failure and the characteristics of micro-seismic signals

    LI Cheng-wu1,SUN Xiao-yuan1,2,GAO Tian-bao3,XIE Bei-jing1,XU Xiao-meng1,4

    (1.School of Resource and Safety Engineering,China University of Mining and Technology (Beijing),Beijing 100083,China;2.College of Environment and Safety,Taiyuan University of Science and Technology,Taiyuan 030024,China;3.School of Energy & Safety Engineering,Tianjin Chengjian University,Tianjin 300384,China;4.Department of Petroleum Engineering,Curtin University,Perth 6151,Australia)

    Abstract:By means of a large vibration apparatus and a micro-seismic observation system,the vibration and failure behavior of coal and rock specimen was analyzed under the conditions of different vibration accelerations and vibration frequencies.Using the Hilbert-Huang Transform (HHT) method,the mode decomposition and signal reconstruction on the monitor signals were made in the process of coal and rock vibration.Combined with the vibration mechanics knowledge,the paper explained the mechanism of resonant amplification effect.The result shows that the monitor signals have obvious carrier characteristics.When the peak acceleration reaches 0.5g,effective micro-seismic signals appear, fractures occur in coal seam for the first time and signals begin to decay with the increase of propagation distance.The natural frequencies are decreased along with the generation and expansion of cracks.When the vibration frequency is equal to the natural frequency,resonant amplification effect happens to specimen and micro-seismic signals grow noticeably.Vibration acceleration and vibration frequency both play roles in the micro-seismic signals.The greater the peak acceleration is,the closer vibration frequency gets to the natural frequency,which makes micro-seismic signals more intensive with the larger amplitude,and has a tendency to shift to low frequency.Theoretical analysis may help

    reveal the causes of these phenomenon.

    Key words:vibration failure;stimulation response;resonance amplification effect;micro-seismic observation;HHT

    作者簡介:李成武(1969—),男,黑龍江延壽人,研究員,博士生導(dǎo)師。E-mail:lcw@ cumtb.edu.cn

    基金項目:國家自然科學(xué)基金資助項目(51274206);國家自然科學(xué)青年基金資助項目(51404277);中央高?;究蒲袠I(yè)務(wù)費專項資金資助項目(2014QZ05)

    收稿日期:2014-09-18

    中圖分類號:TU45

    文獻(xiàn)標(biāo)志碼:A

    文章編號:0253-9993(2015)08-1834-11

    国产精品野战在线观看| 亚洲激情在线av| 九九在线视频观看精品| 国产黄片美女视频| 真实男女啪啪啪动态图| 亚洲成人久久爱视频| 亚洲av第一区精品v没综合| 久久久成人免费电影| 三级国产精品欧美在线观看| 在线国产一区二区在线| 老司机福利观看| 国模一区二区三区四区视频| 亚洲熟妇熟女久久| 搡老岳熟女国产| 午夜福利在线观看免费完整高清在 | 亚洲av中文字字幕乱码综合| 热99re8久久精品国产| 亚洲精品影视一区二区三区av| 我的女老师完整版在线观看| 一区二区三区高清视频在线| 网址你懂的国产日韩在线| 可以在线观看毛片的网站| 成年免费大片在线观看| 久久天躁狠狠躁夜夜2o2o| 黄色视频,在线免费观看| 久久国产乱子免费精品| 久久久久国产精品人妻aⅴ院| 久久久精品大字幕| 亚洲第一欧美日韩一区二区三区| 男插女下体视频免费在线播放| 久久午夜福利片| a在线观看视频网站| 欧美性感艳星| 亚洲国产欧美人成| www日本黄色视频网| 久久久久久九九精品二区国产| 成人永久免费在线观看视频| 国产综合懂色| 国产人妻一区二区三区在| 日本a在线网址| 国产高清视频在线观看网站| 中文字幕人成人乱码亚洲影| 91久久精品电影网| 国产免费男女视频| 高潮久久久久久久久久久不卡| 日韩欧美在线二视频| 51国产日韩欧美| 天天躁日日操中文字幕| 一二三四社区在线视频社区8| 久久人人精品亚洲av| 欧美激情久久久久久爽电影| 免费在线观看日本一区| 蜜桃久久精品国产亚洲av| 一级黄片播放器| 欧美成狂野欧美在线观看| 国产美女午夜福利| 中亚洲国语对白在线视频| 国产亚洲欧美在线一区二区| 久久精品国产自在天天线| 欧美一区二区精品小视频在线| 欧美激情久久久久久爽电影| 国产在线男女| 又爽又黄a免费视频| 两人在一起打扑克的视频| 禁无遮挡网站| 精品免费久久久久久久清纯| 亚洲av电影在线进入| 欧美精品国产亚洲| 少妇丰满av| 国产野战对白在线观看| 嫩草影视91久久| 人人妻人人看人人澡| 欧美一级a爱片免费观看看| 亚洲欧美日韩高清在线视频| 亚洲欧美精品综合久久99| 51国产日韩欧美| 日本免费a在线| 国产91精品成人一区二区三区| 亚洲,欧美,日韩| 欧美国产日韩亚洲一区| 夜夜夜夜夜久久久久| 国产黄a三级三级三级人| 人人妻人人澡欧美一区二区| 搡老熟女国产l中国老女人| 欧美一级a爱片免费观看看| 直男gayav资源| 欧美成人免费av一区二区三区| 日韩 亚洲 欧美在线| 精品免费久久久久久久清纯| 亚洲av一区综合| 精品午夜福利在线看| 少妇被粗大猛烈的视频| 日本五十路高清| 国产精品一区二区性色av| 亚洲专区国产一区二区| 天天一区二区日本电影三级| 精品国内亚洲2022精品成人| 人妻久久中文字幕网| 亚洲成av人片免费观看| 国产成人福利小说| 深夜精品福利| 亚洲无线在线观看| 久久中文看片网| 国产免费av片在线观看野外av| 午夜福利18| 午夜免费成人在线视频| 成人毛片a级毛片在线播放| 久久香蕉精品热| 亚洲av五月六月丁香网| 久久久久免费精品人妻一区二区| 啦啦啦韩国在线观看视频| 国产在视频线在精品| 色av中文字幕| 蜜桃亚洲精品一区二区三区| 国产亚洲欧美98| 国产久久久一区二区三区| 久久人妻av系列| 9191精品国产免费久久| 一进一出抽搐gif免费好疼| 色噜噜av男人的天堂激情| 国产亚洲精品av在线| 欧美日韩综合久久久久久 | 国产视频内射| 亚洲国产精品久久男人天堂| 成人美女网站在线观看视频| 日本三级黄在线观看| 黄色日韩在线| 日韩国内少妇激情av| 国产亚洲精品av在线| 亚洲中文日韩欧美视频| 国产精品野战在线观看| 天天一区二区日本电影三级| 国产免费男女视频| 国产免费男女视频| 看十八女毛片水多多多| 特级一级黄色大片| 久久国产精品影院| 亚洲专区国产一区二区| 在线观看免费视频日本深夜| 婷婷精品国产亚洲av| 国产白丝娇喘喷水9色精品| 日本熟妇午夜| 亚洲无线观看免费| 亚洲最大成人中文| 国产成+人综合+亚洲专区| 国产精品国产高清国产av| 天天一区二区日本电影三级| 欧美日韩瑟瑟在线播放| 一个人看的www免费观看视频| 伊人久久精品亚洲午夜| 精品人妻偷拍中文字幕| 在线观看66精品国产| 午夜两性在线视频| 91久久精品电影网| 少妇被粗大猛烈的视频| 亚洲欧洲国产日韩| 亚洲国产精品专区欧美| 国产精品久久久久久精品电影| 精品久久久久久电影网| 国产成人免费观看mmmm| 99久久中文字幕三级久久日本| 国产精品三级大全| 天美传媒精品一区二区| 亚洲国产精品成人综合色| 国产成人a∨麻豆精品| 综合色丁香网| 一区二区三区四区激情视频| 能在线免费看毛片的网站| 色视频www国产| 香蕉精品网在线| 亚洲人成网站在线播| 高清在线视频一区二区三区| 久久久久久久亚洲中文字幕| 成年人午夜在线观看视频| 国模一区二区三区四区视频| 国产乱来视频区| 久久久久久久亚洲中文字幕| 精品久久久噜噜| 国产又色又爽无遮挡免| 十八禁网站网址无遮挡 | 亚洲欧美日韩无卡精品| 久久99蜜桃精品久久| 国产午夜精品久久久久久一区二区三区| 午夜福利在线在线| 只有这里有精品99| 少妇人妻一区二区三区视频| 99热国产这里只有精品6| 日日撸夜夜添| 国产精品人妻久久久影院| 久久久久久久午夜电影| 久久久久久久精品精品| 大片免费播放器 马上看| 好男人视频免费观看在线| 日韩欧美精品免费久久| 在现免费观看毛片| 99热这里只有是精品在线观看| 国产亚洲一区二区精品| av线在线观看网站| 久久久a久久爽久久v久久| 成人欧美大片| 尾随美女入室| 99久久九九国产精品国产免费| 人体艺术视频欧美日本| 精品熟女少妇av免费看| 亚洲激情五月婷婷啪啪| 深夜a级毛片| 最新中文字幕久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 又爽又黄无遮挡网站| 国产一级毛片在线| 老女人水多毛片| 亚洲精华国产精华液的使用体验| 大香蕉久久网| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频 | 成人二区视频| 在线亚洲精品国产二区图片欧美 | 九草在线视频观看| 久久久久久久亚洲中文字幕| 国产精品伦人一区二区| 97热精品久久久久久| 成人特级av手机在线观看| 久久久久久伊人网av| 中文欧美无线码| 亚洲欧美一区二区三区黑人 | 亚洲自拍偷在线| 精品国产一区二区三区久久久樱花 | 亚洲国产精品专区欧美| 国产男人的电影天堂91| 综合色丁香网| 国产精品熟女久久久久浪| 国产精品熟女久久久久浪| 美女高潮的动态| 又爽又黄a免费视频| 日韩成人伦理影院| 人体艺术视频欧美日本| 男人舔奶头视频| 99久久精品热视频| 成人特级av手机在线观看| videossex国产| 一级黄片播放器| 中文字幕久久专区| 国产欧美亚洲国产| 男人舔奶头视频| 亚洲怡红院男人天堂| 日韩强制内射视频| 大香蕉97超碰在线| 女的被弄到高潮叫床怎么办| 国产淫语在线视频| 女人久久www免费人成看片| 精品人妻偷拍中文字幕| 少妇人妻 视频| 综合色av麻豆| 男插女下体视频免费在线播放| 国产精品久久久久久精品古装| 一级爰片在线观看| 日韩成人伦理影院| 免费观看在线日韩| 91aial.com中文字幕在线观看| 午夜激情久久久久久久| 午夜福利视频精品| 少妇人妻精品综合一区二区| 97在线人人人人妻| 日日摸夜夜添夜夜添av毛片| 精品国产一区二区三区久久久樱花 | 大香蕉97超碰在线| 亚洲欧洲日产国产| 色吧在线观看| 亚洲怡红院男人天堂| 美女主播在线视频| 精品国产乱码久久久久久小说| 欧美日韩视频精品一区| 久久热精品热| 午夜福利网站1000一区二区三区| 高清午夜精品一区二区三区| .国产精品久久| 天天一区二区日本电影三级| 免费av观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲av不卡在线观看| 久久精品夜色国产| 日本欧美国产在线视频| 性插视频无遮挡在线免费观看| 国产精品精品国产色婷婷| 日本猛色少妇xxxxx猛交久久| 亚洲av.av天堂| 看非洲黑人一级黄片| 亚洲欧美精品专区久久| 国产爱豆传媒在线观看| 日本-黄色视频高清免费观看| 日韩av在线免费看完整版不卡| 一级毛片 在线播放| 欧美日韩视频高清一区二区三区二| 免费观看的影片在线观看| 国产69精品久久久久777片| 欧美高清性xxxxhd video| 色5月婷婷丁香| 永久免费av网站大全| 亚洲欧美成人综合另类久久久| 国产精品精品国产色婷婷| 在线 av 中文字幕| 久热久热在线精品观看| 男人添女人高潮全过程视频| 亚洲av二区三区四区| 欧美日韩视频精品一区| 成人特级av手机在线观看| 欧美xxxx黑人xx丫x性爽| 超碰97精品在线观看| 亚洲av中文字字幕乱码综合| 欧美日韩一区二区视频在线观看视频在线 | 人妻一区二区av| 国产高清不卡午夜福利| 亚洲熟女精品中文字幕| 久久人人爽人人爽人人片va| 丝袜美腿在线中文| 美女高潮的动态| 日韩欧美一区视频在线观看 | 成人漫画全彩无遮挡| 国产淫片久久久久久久久| 久久精品国产鲁丝片午夜精品| 99re6热这里在线精品视频| 久久精品人妻少妇| 七月丁香在线播放| 在线观看免费高清a一片| 啦啦啦啦在线视频资源| 少妇的逼水好多| 你懂的网址亚洲精品在线观看| 久久精品国产自在天天线| 一级毛片久久久久久久久女| 青春草国产在线视频| 联通29元200g的流量卡| 大陆偷拍与自拍| 亚洲av成人精品一二三区| 黄色怎么调成土黄色| 性色av一级| 91狼人影院| 亚洲精品,欧美精品| 噜噜噜噜噜久久久久久91| 好男人在线观看高清免费视频| 极品少妇高潮喷水抽搐| 三级男女做爰猛烈吃奶摸视频| 一区二区三区精品91| 激情五月婷婷亚洲| 亚洲欧美中文字幕日韩二区| 老女人水多毛片| 人体艺术视频欧美日本| 美女国产视频在线观看| 色视频在线一区二区三区| 久久久久性生活片| 在线播放无遮挡| 亚洲av国产av综合av卡| 亚洲欧洲日产国产| 好男人在线观看高清免费视频| 国产成人精品一,二区| 人人妻人人澡人人爽人人夜夜| 2021天堂中文幕一二区在线观| 一级a做视频免费观看| av女优亚洲男人天堂| 国产视频内射| 亚洲精品国产av蜜桃| 国产日韩欧美亚洲二区| 国产成人午夜福利电影在线观看| 内射极品少妇av片p| 成人亚洲精品一区在线观看 | 777米奇影视久久| 欧美国产精品一级二级三级 | 黄色怎么调成土黄色| 久久6这里有精品| 99久久精品国产国产毛片| 少妇熟女欧美另类| 久久久久久久国产电影| 99热这里只有是精品50| 欧美+日韩+精品| 三级国产精品欧美在线观看| 国产男人的电影天堂91| 51国产日韩欧美| 在线精品无人区一区二区三 | 婷婷色麻豆天堂久久| 偷拍熟女少妇极品色| 亚洲精品中文字幕在线视频 | 亚洲精品亚洲一区二区| 狂野欧美激情性xxxx在线观看| 少妇熟女欧美另类| 国产成人午夜福利电影在线观看| 蜜臀久久99精品久久宅男| 狂野欧美激情性xxxx在线观看| 三级国产精品片| 久久久久久久久久成人| 婷婷色综合大香蕉| 午夜福利视频精品| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久久电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品嫩草影院av在线观看| 在线观看一区二区三区激情| 永久网站在线| 夫妻午夜视频| 性色av一级| 黄色一级大片看看| 久久久久久伊人网av| 国产精品久久久久久精品电影小说 | 成人漫画全彩无遮挡| 干丝袜人妻中文字幕| h日本视频在线播放| 日日撸夜夜添| 精品一区二区三区视频在线| 美女被艹到高潮喷水动态| 亚洲av国产av综合av卡| 18禁在线播放成人免费| 18禁动态无遮挡网站| 亚洲成色77777| 国产一区亚洲一区在线观看| 国内精品宾馆在线| 国产精品人妻久久久久久| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 日日啪夜夜爽| 你懂的网址亚洲精品在线观看| 国产国拍精品亚洲av在线观看| 乱系列少妇在线播放| 亚洲国产精品999| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 国产色婷婷99| 欧美日韩精品成人综合77777| 久久精品久久精品一区二区三区| 亚洲av在线观看美女高潮| 欧美日韩综合久久久久久| 一级毛片 在线播放| 少妇 在线观看| 免费大片黄手机在线观看| 99久久精品热视频| 亚洲成人久久爱视频| 伦精品一区二区三区| 视频区图区小说| 在线 av 中文字幕| 人人妻人人爽人人添夜夜欢视频 | 99九九线精品视频在线观看视频| 特大巨黑吊av在线直播| 人妻夜夜爽99麻豆av| 天美传媒精品一区二区| 一级毛片aaaaaa免费看小| 在线观看一区二区三区| 久久人人爽人人片av| 网址你懂的国产日韩在线| 少妇的逼好多水| 久久精品国产鲁丝片午夜精品| 免费高清在线观看视频在线观看| 男插女下体视频免费在线播放| 亚洲最大成人中文| 精品久久久久久久久av| 91精品国产九色| 精品熟女少妇av免费看| av在线天堂中文字幕| av线在线观看网站| 久久久久久九九精品二区国产| 亚洲国产精品成人综合色| 少妇的逼水好多| 99久久中文字幕三级久久日本| 三级国产精品片| 久久女婷五月综合色啪小说 | 精品一区二区免费观看| 亚洲精品久久午夜乱码| 尤物成人国产欧美一区二区三区| 成人毛片a级毛片在线播放| 久久韩国三级中文字幕| 可以在线观看毛片的网站| 亚洲精品日韩av片在线观看| 久久久久网色| 日韩国内少妇激情av| 又粗又硬又长又爽又黄的视频| 嘟嘟电影网在线观看| 日韩一本色道免费dvd| 日韩欧美精品免费久久| 亚洲欧美日韩东京热| 一级二级三级毛片免费看| 久久久久久久亚洲中文字幕| 亚洲av成人精品一二三区| 99热全是精品| 亚洲电影在线观看av| 美女xxoo啪啪120秒动态图| 99热这里只有是精品50| 欧美日韩精品成人综合77777| 最近的中文字幕免费完整| 日本-黄色视频高清免费观看| 一区二区三区四区激情视频| 天堂网av新在线| 精品久久久久久电影网| 美女高潮的动态| 在线看a的网站| 亚洲,一卡二卡三卡| 99久久精品一区二区三区| www.色视频.com| 性色av一级| 亚洲av欧美aⅴ国产| 王馨瑶露胸无遮挡在线观看| 欧美日本视频| 国产成年人精品一区二区| 日本熟妇午夜| 国产男女内射视频| 涩涩av久久男人的天堂| 亚洲精品成人久久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲电影在线观看av| 最近中文字幕高清免费大全6| 黄色欧美视频在线观看| 久久久成人免费电影| 国产精品一二三区在线看| 精品熟女少妇av免费看| 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 欧美精品国产亚洲| 日日啪夜夜爽| 欧美日韩精品成人综合77777| 狂野欧美激情性bbbbbb| 欧美三级亚洲精品| 六月丁香七月| 日韩,欧美,国产一区二区三区| 搡老乐熟女国产| 天堂网av新在线| 国产精品秋霞免费鲁丝片| 女人被狂操c到高潮| 欧美精品一区二区大全| 最近手机中文字幕大全| 最近中文字幕2019免费版| 午夜福利视频精品| 国产一区二区三区综合在线观看 | 在线观看一区二区三区激情| 免费av毛片视频| 国产欧美日韩一区二区三区在线 | 日本色播在线视频| 亚洲人与动物交配视频| 免费观看av网站的网址| 久久午夜福利片| 国产精品一区二区三区四区免费观看| av在线天堂中文字幕| 国产精品国产av在线观看| 一边亲一边摸免费视频| 天堂网av新在线| 99久久中文字幕三级久久日本| 99久久人妻综合| 国产精品av视频在线免费观看| 色视频www国产| 日产精品乱码卡一卡2卡三| 亚洲成人久久爱视频| 大话2 男鬼变身卡| 六月丁香七月| 白带黄色成豆腐渣| 91在线精品国自产拍蜜月| 少妇的逼水好多| 国产亚洲午夜精品一区二区久久 | 久久精品久久久久久噜噜老黄| 欧美3d第一页| 国产精品偷伦视频观看了| 少妇人妻精品综合一区二区| 亚洲精品视频女| 国产美女午夜福利| 人妻夜夜爽99麻豆av| 在线免费十八禁| 寂寞人妻少妇视频99o| 国产亚洲精品久久久com| 大陆偷拍与自拍| 一二三四中文在线观看免费高清| av在线亚洲专区| 日韩成人av中文字幕在线观看| 97超碰精品成人国产| 久久精品综合一区二区三区| 看黄色毛片网站| 日韩成人av中文字幕在线观看| 免费av观看视频| 国产黄片视频在线免费观看| 看免费成人av毛片| 亚洲精品久久午夜乱码| 亚洲无线观看免费| 91在线精品国自产拍蜜月| 特级一级黄色大片| 国产精品女同一区二区软件| 中文资源天堂在线| 免费大片18禁| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 麻豆成人av视频| 久久久久久久精品精品| 精品酒店卫生间| 欧美zozozo另类| 少妇高潮的动态图| 亚洲国产精品成人综合色| 一个人看的www免费观看视频| 色播亚洲综合网| 国产精品一区二区在线观看99| 久久午夜福利片| 我的老师免费观看完整版| 国产成人aa在线观看| 国产精品国产三级专区第一集| 亚洲精品456在线播放app| 97人妻精品一区二区三区麻豆| 亚洲精品中文字幕在线视频 | 国产一区二区亚洲精品在线观看| 精品视频人人做人人爽| 又粗又硬又长又爽又黄的视频| 国产探花在线观看一区二区| 麻豆国产97在线/欧美| 小蜜桃在线观看免费完整版高清| 久久ye,这里只有精品| 国产免费又黄又爽又色| 国产国拍精品亚洲av在线观看| 日韩一区二区视频免费看| 激情 狠狠 欧美| 免费av观看视频| 免费大片18禁| 女人久久www免费人成看片| 大片免费播放器 马上看| 欧美区成人在线视频| 亚洲综合精品二区| 国产精品久久久久久久电影|