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

    基于地震波頻譜分析與線性疊加的電子雷管延時(shí)優(yōu)選

    2021-12-18 09:09:10謝先啟韓傳偉鐘冬望王洪剛黃小武黃炳林
    金屬礦山 2021年11期
    關(guān)鍵詞:單孔延時(shí)監(jiān)測(cè)點(diǎn)

    何 理 謝先啟 韓傳偉 鐘冬望 王洪剛 黃小武 黃炳林

    (1.武漢科技大學(xué)冶金工業(yè)過(guò)程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430065;2.爆破工程湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430056;3.武漢爆破有限公司,湖北 武漢 430056)

    鉆爆法施工是現(xiàn)今最為經(jīng)濟(jì)、高效的破巖手段之一,被廣泛應(yīng)用于交通隧道建設(shè)、礦產(chǎn)資源開采與水利水電工程建設(shè)中。研究和應(yīng)用實(shí)踐均表明,使用電子雷管精確延時(shí)可以降低爆破振動(dòng),改善爆破效果,在城鎮(zhèn)爆破開挖及隧道爆破掘進(jìn)中優(yōu)勢(shì)明顯[1-3]。鐘冬望等[4-5]研究指出,電子雷管合理延期時(shí)間往往是一個(gè)或多個(gè)時(shí)間區(qū)間,并且隨著爆心距的變化合理延期時(shí)間并非定值。傅洪賢等[6]結(jié)合隧道穿孔爆破實(shí)踐給出了不同隧道圍巖級(jí)別情況下掏槽孔、周邊孔及輔助孔的建議延期時(shí)間。田振農(nóng)等[7]采用時(shí)—頻分析揭示了隧道爆破振動(dòng)的一般特征,提出了電子雷管錯(cuò)相減震機(jī)理。李清等[8]結(jié)合隧道爆破實(shí)踐,研究指出延期時(shí)間應(yīng)不小于新自由面形成所需要的時(shí)間。凌同華等[9]基于小波變換從實(shí)測(cè)微差爆破振動(dòng)信號(hào)中分離出分段振波,通過(guò)不同延期時(shí)間疊加各分段振波達(dá)到優(yōu)選孔間延期時(shí)間的目的。趙明生等[10]、楊年華等[11]基于實(shí)測(cè)單段波形線性疊加分析,優(yōu)選孔間延期時(shí)間,實(shí)現(xiàn)了地震波的干擾相消。然而迄今為止,對(duì)孔間延期時(shí)間的設(shè)置尚未形成統(tǒng)一的認(rèn)識(shí),各種延期時(shí)間的優(yōu)選方法普適性均有所不足,施工現(xiàn)場(chǎng)更多依賴經(jīng)驗(yàn)和多次試驗(yàn)來(lái)確定電子雷管延期時(shí)間[12-13]。并且,往往爆區(qū)待保護(hù)對(duì)象有多個(gè),但是電子雷管延期時(shí)間設(shè)置需要考慮爆心距因素,現(xiàn)有研究成果多是針對(duì)特定建(構(gòu))筑物提出,因此,亟需解決不同爆心距處被保護(hù)對(duì)象延期時(shí)間的優(yōu)選問(wèn)題。

    結(jié)合國(guó)內(nèi)某紅砂巖石方開挖爆破工程,本研究設(shè)計(jì)開展電子雷管爆破振動(dòng)監(jiān)測(cè)試驗(yàn),通過(guò)頻譜分析討論地震波主振頻率演變規(guī)律;基于地震波線性疊加原理,運(yùn)用MATLAB編制計(jì)算程序,研究疊加波列數(shù)與孔間延期時(shí)間對(duì)合成波形峰值振動(dòng)速度的影響機(jī)制,建立電子雷管延期時(shí)間的確定方法。研究成果可為城市復(fù)雜環(huán)境下巖石開挖爆破電子雷管延期時(shí)間設(shè)置提供參考。

    1 爆破振動(dòng)測(cè)試試驗(yàn)與分析

    1.1 場(chǎng)地概況

    爆破振動(dòng)測(cè)試試驗(yàn)在國(guó)內(nèi)某紅砂巖石方開挖爆破工程一平整場(chǎng)地進(jìn)行,由于該工程部分工作面緊鄰村莊、學(xué)校等建(構(gòu))筑物,因此需要重點(diǎn)關(guān)注爆破有害效應(yīng),尤其是實(shí)現(xiàn)爆破振動(dòng)的精細(xì)控制極為必要。試驗(yàn)前,在工程現(xiàn)場(chǎng)鉆孔取芯選用密實(shí)均勻和完整性良好的紅砂巖制備試件,按照巖石力學(xué)的試驗(yàn)測(cè)試要求,試件尺寸為50 mm×50 mm(直徑×高度),并且對(duì)試樣端面進(jìn)行打磨,確保端面法線偏差小于0.25°,試樣端面不平整度小于0.02 mm。室內(nèi)試驗(yàn)測(cè)得紅砂巖試件的密度為2.55 g/cm3,縱波波速為2 121 m/s,單軸抗壓強(qiáng)度為60 MPa。

    1.2 爆破參數(shù)與測(cè)振方案

    根據(jù)實(shí)際生產(chǎn)爆破情況,單孔爆破試驗(yàn)時(shí)垂直下向鉆孔,選用連續(xù)裝藥結(jié)構(gòu),通過(guò)單發(fā)電子雷管孔底反向起爆。炮孔直徑為90 mm,炮孔深度為2.5 m,單孔裝藥量4 kg,堵塞長(zhǎng)度為1.5 m。炮孔裝藥結(jié)構(gòu)如圖1所示。

    為獲取可真實(shí)反映紅砂巖地質(zhì)條件屬性的單孔爆破振動(dòng)信號(hào),選取較為平齊、完整的工程區(qū)域范圍,在爆心距分別為10、15、20、25 m處共布置了4個(gè)振動(dòng)監(jiān)測(cè)點(diǎn),每個(gè)振動(dòng)監(jiān)測(cè)點(diǎn)放置2臺(tái)振動(dòng)測(cè)試儀器,且4個(gè)振動(dòng)監(jiān)測(cè)點(diǎn)均布置在同一條測(cè)線上,目的是盡可能減小傳播路徑差異(傳播介質(zhì)屬性)對(duì)振動(dòng)測(cè)試波形的影響。試驗(yàn)完成后分別相互對(duì)比分析4個(gè)振動(dòng)監(jiān)測(cè)點(diǎn)處2臺(tái)振動(dòng)測(cè)試儀的監(jiān)測(cè)指標(biāo)(峰值振動(dòng)速度、振動(dòng)持續(xù)時(shí)間、振動(dòng)主頻),如果任一指標(biāo)間的相對(duì)誤差超過(guò)5%,則重復(fù)進(jìn)行爆破振動(dòng)監(jiān)測(cè)試驗(yàn),直到監(jiān)測(cè)指標(biāo)趨于穩(wěn)定。

    1.3 測(cè)試結(jié)果與分析

    試驗(yàn)測(cè)得的單孔爆破振動(dòng)波形如圖2所示。

    由于爆破試驗(yàn)在平整場(chǎng)地開展,無(wú)需考慮高程差因素對(duì)爆破振動(dòng)傳播的影響,故選取如下形式的PPV預(yù)測(cè)模型:

    式中,PPV為質(zhì)點(diǎn)峰值振動(dòng)速度,cm/s;R為爆心距,m;Q為最大單響藥量,kg;k、α是場(chǎng)地條件相關(guān)系數(shù),與爆破方式、裝藥結(jié)構(gòu)、爆破點(diǎn)至監(jiān)測(cè)點(diǎn)間的地形、地質(zhì)條件密切相關(guān)[14];SD為比例距離[15],m/kg1/3。

    同時(shí),利用MATLAB編制計(jì)算程序,對(duì)振動(dòng)信號(hào)進(jìn)行頻譜分析,得到不同爆心距處振動(dòng)信號(hào)的主振頻率,如圖3所示。

    爆破地震波峰值振動(dòng)速度PPV與振動(dòng)主頻F隨比例距離的變化關(guān)系如圖4所示。

    根據(jù)圖4,可擬合得到紅砂巖場(chǎng)地爆破振動(dòng)速度的預(yù)測(cè)模型為

    數(shù)據(jù)擬合相關(guān)系數(shù)達(dá)到0.951 4,表明該模型回歸擬合效果良好,式(2)中的場(chǎng)地條件相關(guān)系數(shù)可真實(shí)表征施工現(xiàn)場(chǎng)的地質(zhì)條件屬性。

    并且,可得到爆破振動(dòng)主頻隨比例距離增加近似呈冪函數(shù)形式衰減,振動(dòng)主頻F衰減公式為

    2 單孔爆破振動(dòng)曲線擬合及其線性疊加原理

    單孔爆破獲得的不同爆心距處爆破地震波形中包含了爆源到振動(dòng)監(jiān)測(cè)點(diǎn)傳播介質(zhì)的地質(zhì)屬性。假設(shè)單孔爆破地震波形是由無(wú)數(shù)個(gè)正弦諧波組成,諧波間相互干擾疊加,正、正相位間波形幅值疊加相長(zhǎng),致使地震波幅值增加,正、負(fù)相位間波形幅值干擾相消,致使地震波幅值降低[16-19]。群孔爆破地震波可視為由特定數(shù)目的單孔地震波形線性疊加組成[20],疊加方法如下式:

    式中,V(t)為線性疊加原理預(yù)測(cè)得到的群孔爆破地震波形;t為時(shí)間;n為爆破地震波列數(shù);vi(t)為單孔爆破地震波形;i為地震波形段號(hào),i=1,2,3,…,n;Δt為段間延期時(shí)間。

    利用MATLAB編程計(jì)算程序,基于Fourier級(jí)數(shù)擬合得到單孔地震波形的擬合函數(shù)?(t)[21],函數(shù)形式為

    式中,?(t) 為單孔地震波形擬合函數(shù);a0、aj、bj為Fourier擬合系數(shù);ω為基頻;l為Fourier擬合級(jí)數(shù)。

    式(5)中,a0、aj、bj以及基頻ω的計(jì)算公式為

    式中,T為波形截?cái)鄷r(shí)間;M為總采樣點(diǎn)數(shù);ym為第m個(gè)采樣值。

    將擬合函數(shù)?()t擴(kuò)展至?xí)r間全域,得到最終單孔地震波形函數(shù)v()t為

    3 波列數(shù)與延期時(shí)間對(duì)合成波形峰值振動(dòng)速度的影響機(jī)制

    3.1 波列數(shù)對(duì)合成波形峰值振動(dòng)速度的影響規(guī)律

    生產(chǎn)爆破實(shí)踐中,現(xiàn)場(chǎng)實(shí)測(cè)波形往往是多個(gè)炮孔爆炸后各波形的疊加結(jié)果。假設(shè)每段產(chǎn)生的單孔波形相同,當(dāng)段數(shù)不同時(shí),相同延期時(shí)間下合成(疊加)波形的峰值振動(dòng)速度也可能存在差異[17]。因此,利用地震波線性疊加法分析延期時(shí)間和合成波形峰值振動(dòng)速度間的關(guān)系時(shí),首先消除疊加波列數(shù)對(duì)合成波形峰值振動(dòng)速度的影響極為必要。本研究基于地震波線性疊加原理,利用MATLAB編制計(jì)算程序,疊加波列數(shù)取為n∈[ ]2,40,且為整數(shù);各波列間延期時(shí)間取為Δt∈[ ]1,80ms,且為整數(shù)。

    以圖2中10 m監(jiān)測(cè)點(diǎn)處單孔地震波形為例,當(dāng)疊加波列數(shù)為25時(shí),合成波形峰值振動(dòng)速度隨延期時(shí)間的變化關(guān)系如圖5所示。

    根據(jù)圖5可知,單孔爆破地震波形在一定波列數(shù)疊加情況下,不同延期時(shí)間疊加產(chǎn)生的合成波形具有不同的峰值振動(dòng)速度。延期時(shí)間為1 ms或2 ms時(shí),合成波形峰值振動(dòng)速度小于單孔振波峰值振動(dòng)速度;當(dāng)延期時(shí)間超過(guò)特定值后,合成波形的峰值振動(dòng)速度逐漸趨近并最終等于單孔爆破峰值振動(dòng)速度。合成波形的峰值振動(dòng)速度最小值與最大值分別出現(xiàn)在延期時(shí)間Δt=1 ms和Δt=19 ms時(shí),此時(shí)的合成波形如圖6所示。

    通過(guò)前述分析計(jì)算方法,依次計(jì)算并統(tǒng)計(jì)出任意延期時(shí)間下合成波形的峰值振動(dòng)速度最大值及最小值,得到合成波形峰值振動(dòng)速度隨波列數(shù)的變化關(guān)系,如圖7所示。

    依據(jù)圖7可以得到,延時(shí)1~80 ms范圍內(nèi),各監(jiān)測(cè)點(diǎn)合成波形峰值振動(dòng)速度最大值隨波列數(shù)增加,呈現(xiàn)出先增大后減小,并最終趨于穩(wěn)定值的變化規(guī)律。對(duì)于合成波形峰值振動(dòng)速度的最小值,除了10 m爆心距處出現(xiàn)最小值隨波列數(shù)上下波動(dòng),而后逐漸趨近單孔振波峰值振動(dòng)速度外,其它監(jiān)測(cè)點(diǎn)處最小值始終維持在單孔振波的峰值振動(dòng)速度水平。進(jìn)一步分析發(fā)現(xiàn),最小值小于單孔振波峰值的情況僅僅出現(xiàn)在Δt=1 ms和Δt=2 ms處。然而在工程實(shí)踐中,為取得良好的巖石破碎效果,通常設(shè)置段間延時(shí)不小于7 ms[16],故本研究在分析波列數(shù)對(duì)合成波形峰值振動(dòng)速度影響時(shí),可忽略其對(duì)合成波形峰值振動(dòng)速度最小值的影響。

    根據(jù)圖7統(tǒng)計(jì)出的合成波形峰值振動(dòng)速度最大值趨于穩(wěn)定的最小臨界波列數(shù),見表1。

    ?

    由圖7和表1可知:合成波形峰值振動(dòng)速度最大值趨于穩(wěn)定的臨界波列數(shù)隨監(jiān)測(cè)點(diǎn)爆心距的增大逐漸增加,試驗(yàn)條件下當(dāng)波列數(shù)n≥23時(shí),在監(jiān)測(cè)范圍內(nèi)可忽略其對(duì)合成波形峰值振動(dòng)速度的影響,此時(shí)可僅考慮延期時(shí)間對(duì)合成波形峰值振動(dòng)速度的影響。臨界波列數(shù)隨比例距離的變化關(guān)系如圖8所示。

    由圖8可擬合得到合成波形峰值振動(dòng)速度趨于穩(wěn)定的臨界波列數(shù)N的趨勢(shì)線方程為

    式中,N≥2且取整數(shù)。

    利用式(8)可計(jì)算不同監(jiān)測(cè)點(diǎn)處的臨界波列數(shù),研究結(jié)果可為紅砂巖地質(zhì)情況下,采用地震波線性疊加法預(yù)測(cè)爆破振動(dòng)速度峰值及優(yōu)選合理延期時(shí)間提供理論參考。

    3.2 延期時(shí)間對(duì)合成波形峰值振動(dòng)速度的影響規(guī)律

    以15 m監(jiān)測(cè)點(diǎn)處的單孔地震波形為例,以不同延期時(shí)間疊加23個(gè)波列,可得到合成波形峰值振動(dòng)速度隨延期時(shí)間的變化關(guān)系如圖9所示。

    由圖9可知:針對(duì)露天臺(tái)階爆破工程實(shí)踐,忽略段間延時(shí)小于7 ms的情況[16],當(dāng)孔間延時(shí)為12 ms時(shí),合成波形峰值振動(dòng)速度首次達(dá)到最小值3.739 cm/s;當(dāng)孔間延時(shí)為9 ms和18 ms時(shí),合成波形峰值振動(dòng)速度達(dá)到較小值分別為4.027 cm/s、3.96 cm/s;當(dāng)孔間延期時(shí)間為15 ms時(shí),合成波形峰值振動(dòng)速度達(dá)到最大值4.767 cm/s。

    當(dāng)孔間延期時(shí)間分別為12 ms和15 ms時(shí),合成波形如圖10所示。

    由圖10可知:當(dāng)孔間延時(shí)為12 ms時(shí),合成波形峰值振動(dòng)速度出現(xiàn)在25.5 ms處,即第1和第2個(gè)波形產(chǎn)生了較弱的干擾疊加,峰值振動(dòng)速度基本與單孔振波峰值振動(dòng)速度相當(dāng),而后續(xù)波形產(chǎn)生了干擾相消。當(dāng)孔間延期時(shí)間為15 ms時(shí),合成波形峰值振動(dòng)速度出現(xiàn)在28.5 ms處,表明第1和第2個(gè)波形產(chǎn)生了較強(qiáng)程度的干擾疊加,此時(shí)峰值振動(dòng)速度達(dá)到單孔振波峰值的1.28倍,后期波形干擾疊加程度略微降低。

    利用上述方法進(jìn)行分析計(jì)算,依次統(tǒng)計(jì)出各監(jiān)測(cè)點(diǎn)處合成波形峰值振動(dòng)速度出現(xiàn)時(shí)刻,結(jié)果見表2。

    ?

    由表2可知:各監(jiān)測(cè)點(diǎn)處合成波形峰值振動(dòng)速度最小值和較小值所對(duì)應(yīng)的延期時(shí)間值多為相應(yīng)測(cè)點(diǎn)單孔地震波形主振周期的0.4~0.6倍(0.5倍左右),說(shuō)明通過(guò)設(shè)置合理的延期時(shí)間,利用電子雷管錯(cuò)相降振是可行的。然而,合成波形峰值振動(dòng)速度最小值通常僅為單孔波形的峰值振動(dòng)速度。當(dāng)合成波形峰值振動(dòng)速度出現(xiàn)最大值時(shí),對(duì)應(yīng)的延期時(shí)間值約為單孔振波主振周期的0.9~1.1倍(1倍左右),此時(shí)分段振波間存在振速疊加增強(qiáng)現(xiàn)象。

    另外,當(dāng)孔間延時(shí)達(dá)到主振周期約2.6倍以上時(shí),合成波形峰值振動(dòng)速度逐漸趨于穩(wěn)定值(單孔振波峰值振動(dòng)速度),各分段振波主震相之間的疊加干擾效應(yīng)微乎其微,不會(huì)對(duì)合成波形峰值振動(dòng)速度造成影響,工程實(shí)踐中可借助這一原則對(duì)地震波主震相持續(xù)時(shí)間進(jìn)行估算。進(jìn)一步可得到確保合成波形峰值振動(dòng)速度趨于穩(wěn)定的主振周期倍數(shù)β隨比例距離的變化關(guān)系如圖11所示。

    4 電子雷管延期時(shí)間優(yōu)選方法

    根據(jù)單孔爆破地震波形主振頻率隨比例距離的衰減公式,可以推導(dǎo)出孔間延期時(shí)間Δt的計(jì)算公式:

    可以明確的是,當(dāng)孔間延期時(shí)間為0.5倍的主振周期時(shí),并不能完全確保合成波形峰值振動(dòng)速度最小,只能說(shuō)延期時(shí)間為主振周期的0.4~0.6倍時(shí),合成波形峰值振動(dòng)速度可能會(huì)出現(xiàn)最小值或較小值。然而,當(dāng)延期時(shí)間為單孔振波主振周期的1倍左右時(shí),合成波形必然會(huì)由于干擾疊加產(chǎn)生振速增強(qiáng)現(xiàn)象[17]。為了增強(qiáng)分段振波干擾降振的可靠性,使合成波形峰值振動(dòng)速度盡可能降低,本研究提出了一種不同監(jiān)測(cè)點(diǎn)處孔間最優(yōu)延時(shí)確定方法,計(jì)算流程如圖12所示。

    計(jì)算步驟為:首先通過(guò)單孔爆破試驗(yàn)獲取可表征現(xiàn)場(chǎng)地質(zhì)條件與震動(dòng)屬性的單孔爆破地震波形,基于Fourier級(jí)數(shù)擬合得到單孔爆破地震波形的擬合函數(shù)v(t),通過(guò)地震波線性疊加原理計(jì)算不同延期時(shí)間下的合成波形,從而初步選取孔間延期時(shí)間值Δt;然后,根據(jù)判斷準(zhǔn)則1(孔間延期時(shí)間Δt若滿足確定該值為孔間最優(yōu)延期時(shí)間Δt優(yōu);否則,根據(jù)判斷準(zhǔn)則2(孔間延期時(shí)間值Δt若滿足確定該值為孔間最優(yōu)延期時(shí)間Δt優(yōu)。如果判斷準(zhǔn)則1和判斷準(zhǔn)則2均不滿足,則返回前述步驟,直至得到孔間最優(yōu)延期時(shí)間值Δt優(yōu)為止。

    5 結(jié)論

    (1)試驗(yàn)研究得到合成波形峰值振動(dòng)速度趨于穩(wěn)定值的最小臨界波列數(shù)N為23,并建立了臨界波列數(shù)N的趨勢(shì)線方程,其隨比例距離近似呈對(duì)數(shù)形式增加,可為紅砂巖地質(zhì)情況下爆破振動(dòng)速度峰值預(yù)測(cè)及合理延期時(shí)間優(yōu)選提供理論依據(jù)。

    (2)當(dāng)孔間延時(shí)為主振周期的0.4~0.6倍時(shí),合成波形峰值振動(dòng)速度可能會(huì)出現(xiàn)最小值或較小值;當(dāng)孔間延時(shí)為單孔振波主振周期的1倍左右時(shí),合成波形必然會(huì)由于干擾疊加產(chǎn)生振速增強(qiáng)現(xiàn)象。當(dāng)孔間延時(shí)達(dá)到單孔振波主振周期的2.6倍以上時(shí),各分段振波主震相間疊加干擾效應(yīng)減弱,合成波形峰值振動(dòng)速度接近單孔爆破峰值振動(dòng)速度,實(shí)踐中可據(jù)此對(duì)地震波主震相持續(xù)時(shí)間進(jìn)行估算。

    (3)考慮爆心距因素對(duì)孔(段)間延時(shí)優(yōu)選的影響,提出了一種不同監(jiān)測(cè)點(diǎn)處孔(段)間最優(yōu)延期時(shí)間確定方法,擺脫了傳統(tǒng)方法僅能估算特定監(jiān)測(cè)點(diǎn)處最優(yōu)延時(shí)的限制,工程適用性更強(qiáng)。

    猜你喜歡
    單孔延時(shí)監(jiān)測(cè)點(diǎn)
    天津南港LNG接收站沉降監(jiān)測(cè)點(diǎn)位布設(shè)
    煤氣與熱力(2022年4期)2022-05-23 12:44:56
    撫河流域綜合治理監(jiān)測(cè)布局優(yōu)化
    基于級(jí)聯(lián)步進(jìn)延時(shí)的順序等效采樣方法及實(shí)現(xiàn)
    全站儀極坐標(biāo)法監(jiān)測(cè)點(diǎn)穩(wěn)定性分析方法研究
    我省舉辦家畜血吸蟲病監(jiān)測(cè)點(diǎn)培訓(xùn)班
    單孔腹腔鏡手術(shù)后臍窩創(chuàng)面的簡(jiǎn)單縫合術(shù)
    單孔腹腔鏡手術(shù)切除左位膽囊1例報(bào)告
    經(jīng)臍單孔腹腔鏡在普外手術(shù)中的應(yīng)用分析
    Two-dimensional Eulerian-Lagrangian Modeling of Shocks on an Electronic Package Embedded in a Projectile with Ultra-high Acceleration
    單孔后腹腔鏡下輸尿管切開取石術(shù)的初步應(yīng)用(附23例報(bào)告)
    最近最新免费中文字幕在线| 黄色女人牲交| 亚洲国产精品sss在线观看| 国产午夜福利久久久久久| 亚洲av五月六月丁香网| 欧美激情国产日韩精品一区| 国产久久久一区二区三区| 成人鲁丝片一二三区免费| 能在线免费观看的黄片| 久久久国产成人免费| 又爽又黄a免费视频| 日韩大尺度精品在线看网址| 人妻丰满熟妇av一区二区三区| 欧美国产日韩亚洲一区| 成年女人毛片免费观看观看9| xxxwww97欧美| 一进一出抽搐动态| 亚洲欧美激情综合另类| 国产亚洲精品综合一区在线观看| 人人妻人人澡欧美一区二区| 露出奶头的视频| 欧美激情在线99| 亚洲aⅴ乱码一区二区在线播放| 国产一区二区亚洲精品在线观看| a级毛片免费高清观看在线播放| 深夜精品福利| 免费不卡的大黄色大毛片视频在线观看 | 国产三级在线视频| 91久久精品电影网| 欧美最新免费一区二区三区| 天天一区二区日本电影三级| 免费观看的影片在线观看| 欧美成人一区二区免费高清观看| 精品久久久久久久久久久久久| 国产色爽女视频免费观看| 99久久成人亚洲精品观看| 免费看av在线观看网站| 久久精品国产99精品国产亚洲性色| 91麻豆av在线| 欧美高清性xxxxhd video| 国产高清三级在线| 老熟妇乱子伦视频在线观看| 午夜福利欧美成人| 69人妻影院| 国产成人aa在线观看| 能在线免费观看的黄片| 国产精品久久视频播放| 国产黄色小视频在线观看| av在线天堂中文字幕| 午夜爱爱视频在线播放| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在 | 精品久久久久久久末码| 国产真实乱freesex| 一区二区三区四区激情视频 | 国产精品久久电影中文字幕| 三级男女做爰猛烈吃奶摸视频| 日本 av在线| 69av精品久久久久久| 国产精品久久久久久久久免| 免费无遮挡裸体视频| 嫩草影视91久久| 久久久午夜欧美精品| 深夜a级毛片| 日本色播在线视频| 国产高清激情床上av| 大型黄色视频在线免费观看| 91在线精品国自产拍蜜月| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 日韩强制内射视频| 淫秽高清视频在线观看| 久久久久久大精品| 丝袜美腿在线中文| 日韩欧美在线二视频| 国产白丝娇喘喷水9色精品| 午夜a级毛片| 亚洲最大成人av| 欧美在线一区亚洲| 亚洲成人久久爱视频| 精品日产1卡2卡| 在线播放国产精品三级| 波多野结衣巨乳人妻| 日韩,欧美,国产一区二区三区 | 午夜久久久久精精品| 观看免费一级毛片| 啦啦啦韩国在线观看视频| 在线观看舔阴道视频| 亚洲一区二区三区色噜噜| 亚洲精品成人久久久久久| 久久亚洲精品不卡| 久久精品国产亚洲av香蕉五月| 亚洲精品国产成人久久av| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 国产三级中文精品| 69av精品久久久久久| 久久久久久久久久久丰满 | 精品久久国产蜜桃| 一边摸一边抽搐一进一小说| 欧美色视频一区免费| 春色校园在线视频观看| 女人被狂操c到高潮| 在线播放国产精品三级| 观看免费一级毛片| 不卡一级毛片| 亚洲国产欧美人成| 欧美xxxx性猛交bbbb| 久久久久久久久中文| 在线观看av片永久免费下载| 久久久久久国产a免费观看| 亚洲七黄色美女视频| 人妻少妇偷人精品九色| 精品午夜福利视频在线观看一区| 一本精品99久久精品77| 能在线免费观看的黄片| 国产精品99久久久久久久久| 日本 av在线| 一区二区三区免费毛片| 成人一区二区视频在线观看| 成人三级黄色视频| 嫁个100分男人电影在线观看| av黄色大香蕉| 色5月婷婷丁香| 变态另类成人亚洲欧美熟女| 久久精品人妻少妇| 亚洲美女黄片视频| 99久久精品热视频| 日本 欧美在线| 国产aⅴ精品一区二区三区波| 欧美三级亚洲精品| 欧美日韩中文字幕国产精品一区二区三区| 久久久国产成人精品二区| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 在线观看一区二区三区| 国产私拍福利视频在线观看| 天天躁日日操中文字幕| 九色国产91popny在线| 两个人的视频大全免费| 深爱激情五月婷婷| 亚洲avbb在线观看| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 三级毛片av免费| 国产又黄又爽又无遮挡在线| 老师上课跳d突然被开到最大视频| 中文字幕高清在线视频| 一个人看的www免费观看视频| 男女做爰动态图高潮gif福利片| 最后的刺客免费高清国语| 欧美一级a爱片免费观看看| 欧美+亚洲+日韩+国产| 日本爱情动作片www.在线观看 | 精品国内亚洲2022精品成人| 日本黄大片高清| 日韩中字成人| 深夜精品福利| 免费人成视频x8x8入口观看| 18禁黄网站禁片午夜丰满| 99热精品在线国产| 少妇被粗大猛烈的视频| 久99久视频精品免费| 国产伦在线观看视频一区| 男女做爰动态图高潮gif福利片| 男人和女人高潮做爰伦理| 中文字幕高清在线视频| 色在线成人网| 亚洲av美国av| 岛国在线免费视频观看| 午夜精品在线福利| 欧美中文日本在线观看视频| 亚洲精华国产精华液的使用体验 | 身体一侧抽搐| 久久久久久久精品吃奶| 可以在线观看毛片的网站| 美女免费视频网站| 别揉我奶头~嗯~啊~动态视频| 99久久无色码亚洲精品果冻| 在线免费十八禁| 精品久久久久久久久久免费视频| 日本黄色视频三级网站网址| 日韩欧美精品v在线| 免费看美女性在线毛片视频| 给我免费播放毛片高清在线观看| 老司机福利观看| 中文字幕高清在线视频| 国产高清视频在线播放一区| 国产成年人精品一区二区| 欧美色视频一区免费| 搞女人的毛片| 国产黄a三级三级三级人| 婷婷丁香在线五月| 亚洲国产高清在线一区二区三| 露出奶头的视频| 亚洲图色成人| 国产午夜精品久久久久久一区二区三区 | 男女边吃奶边做爰视频| 深爱激情五月婷婷| 又爽又黄无遮挡网站| 日韩一区二区视频免费看| 欧美绝顶高潮抽搐喷水| 午夜福利在线观看免费完整高清在 | 免费av不卡在线播放| 床上黄色一级片| 久久草成人影院| 身体一侧抽搐| 嫩草影院入口| 男人狂女人下面高潮的视频| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 国产成人a区在线观看| 国产三级在线视频| 国产精品一区www在线观看 | 51国产日韩欧美| 婷婷精品国产亚洲av| 国产精品日韩av在线免费观看| 亚洲男人的天堂狠狠| 国产一区二区三区视频了| 91在线精品国自产拍蜜月| 啦啦啦观看免费观看视频高清| 国产精品无大码| 亚洲av免费高清在线观看| 国产免费男女视频| 国国产精品蜜臀av免费| 亚洲久久久久久中文字幕| 国产精品综合久久久久久久免费| 亚洲精品久久国产高清桃花| 久久久久国产精品人妻aⅴ院| 国内精品一区二区在线观看| 麻豆国产97在线/欧美| 内地一区二区视频在线| 亚洲专区中文字幕在线| 成年人黄色毛片网站| 又爽又黄a免费视频| 日韩精品青青久久久久久| 亚洲精华国产精华精| 亚洲成av人片在线播放无| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 老熟妇仑乱视频hdxx| 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久午夜电影| 国产成人a区在线观看| eeuss影院久久| 亚洲国产色片| 亚洲中文字幕一区二区三区有码在线看| 天美传媒精品一区二区| 免费av观看视频| 88av欧美| 老师上课跳d突然被开到最大视频| 色在线成人网| 久久久久性生活片| 一区二区三区高清视频在线| 毛片女人毛片| 九九在线视频观看精品| 91久久精品国产一区二区成人| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 国产色婷婷99| 精品乱码久久久久久99久播| 国产av一区在线观看免费| 国产精品无大码| 午夜爱爱视频在线播放| 成年女人毛片免费观看观看9| 日本黄大片高清| 毛片一级片免费看久久久久 | 国产成年人精品一区二区| 欧美极品一区二区三区四区| 午夜福利在线在线| 国产成人av教育| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 黄色欧美视频在线观看| 色综合亚洲欧美另类图片| 亚洲成a人片在线一区二区| 精品一区二区三区人妻视频| 国产精品久久久久久av不卡| 日韩精品中文字幕看吧| 亚洲七黄色美女视频| 18+在线观看网站| 亚洲av美国av| 国产 一区 欧美 日韩| 一a级毛片在线观看| 无人区码免费观看不卡| 1000部很黄的大片| 在现免费观看毛片| 亚洲中文字幕一区二区三区有码在线看| 夜夜爽天天搞| 成年版毛片免费区| 亚洲av熟女| 在线观看舔阴道视频| 床上黄色一级片| 最近最新中文字幕大全电影3| av国产免费在线观看| 精品无人区乱码1区二区| av.在线天堂| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 观看美女的网站| 中文字幕av在线有码专区| 国产欧美日韩精品亚洲av| av黄色大香蕉| 亚洲av.av天堂| 亚洲久久久久久中文字幕| 国产一区二区三区av在线 | 久久午夜亚洲精品久久| 精品一区二区三区视频在线观看免费| 午夜老司机福利剧场| 日韩,欧美,国产一区二区三区 | 国产亚洲精品久久久久久毛片| 丝袜美腿在线中文| 白带黄色成豆腐渣| 男女视频在线观看网站免费| 欧美日韩瑟瑟在线播放| 女人十人毛片免费观看3o分钟| 国产成人一区二区在线| 深爱激情五月婷婷| 欧美人与善性xxx| or卡值多少钱| 日日摸夜夜添夜夜添av毛片 | 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| 在线播放无遮挡| 国产一区二区亚洲精品在线观看| 免费一级毛片在线播放高清视频| 波多野结衣高清作品| 尤物成人国产欧美一区二区三区| 他把我摸到了高潮在线观看| 欧洲精品卡2卡3卡4卡5卡区| 日本 av在线| 色av中文字幕| 最新中文字幕久久久久| 亚洲精品乱码久久久v下载方式| 少妇被粗大猛烈的视频| 蜜桃久久精品国产亚洲av| 日韩中字成人| 国产精品久久电影中文字幕| 九九久久精品国产亚洲av麻豆| 色av中文字幕| 亚洲美女视频黄频| 日本在线视频免费播放| 国产精品一及| 国产色爽女视频免费观看| 亚洲av.av天堂| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 特大巨黑吊av在线直播| 久久精品久久久久久噜噜老黄 | 久久人人爽人人爽人人片va| 成人欧美大片| 国国产精品蜜臀av免费| 日韩在线高清观看一区二区三区 | 国内精品久久久久精免费| 亚州av有码| 真实男女啪啪啪动态图| 我的女老师完整版在线观看| 国产高清三级在线| АⅤ资源中文在线天堂| 十八禁国产超污无遮挡网站| 国产精品人妻久久久影院| 中文在线观看免费www的网站| 久久这里只有精品中国| 国产大屁股一区二区在线视频| 国产精品av视频在线免费观看| 欧美成人免费av一区二区三区| 精华霜和精华液先用哪个| 免费看日本二区| 少妇人妻一区二区三区视频| 成人鲁丝片一二三区免费| 亚洲自拍偷在线| 18禁黄网站禁片午夜丰满| 久久热精品热| 日本撒尿小便嘘嘘汇集6| 国内精品美女久久久久久| 日本爱情动作片www.在线观看 | 不卡一级毛片| 久久久久久久精品吃奶| 色综合色国产| 精品久久久久久久末码| 国产精品久久久久久久电影| 禁无遮挡网站| 国产极品精品免费视频能看的| 国模一区二区三区四区视频| 国产午夜福利久久久久久| 国产精品嫩草影院av在线观看 | 精品久久久久久久末码| 成人特级黄色片久久久久久久| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 99久久精品国产国产毛片| 国产极品精品免费视频能看的| 99在线视频只有这里精品首页| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 国产高清视频在线播放一区| 国产成人av教育| 神马国产精品三级电影在线观看| 桃色一区二区三区在线观看| 少妇被粗大猛烈的视频| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 级片在线观看| 日本精品一区二区三区蜜桃| 成人精品一区二区免费| 国内毛片毛片毛片毛片毛片| 成年版毛片免费区| 国产精华一区二区三区| 永久网站在线| 国产精品美女特级片免费视频播放器| 午夜影院日韩av| 极品教师在线视频| 一进一出好大好爽视频| 99热精品在线国产| 亚洲最大成人av| av在线亚洲专区| 69人妻影院| 国产精品1区2区在线观看.| 99riav亚洲国产免费| 身体一侧抽搐| 国产成人aa在线观看| 欧美+亚洲+日韩+国产| bbb黄色大片| 国产精品综合久久久久久久免费| 两人在一起打扑克的视频| 亚洲国产精品sss在线观看| 亚洲欧美精品综合久久99| 三级国产精品欧美在线观看| 97超级碰碰碰精品色视频在线观看| 又粗又爽又猛毛片免费看| 97超级碰碰碰精品色视频在线观看| 亚洲自拍偷在线| 欧美另类亚洲清纯唯美| 一进一出好大好爽视频| 窝窝影院91人妻| 91久久精品国产一区二区三区| 国产亚洲精品av在线| 最近中文字幕高清免费大全6 | 日韩欧美免费精品| 国产精品98久久久久久宅男小说| 99riav亚洲国产免费| 身体一侧抽搐| 亚洲成a人片在线一区二区| 欧美日韩精品成人综合77777| 色尼玛亚洲综合影院| 深夜精品福利| 亚洲午夜理论影院| 欧美在线一区亚洲| 夜夜看夜夜爽夜夜摸| 人妻少妇偷人精品九色| 国产真实伦视频高清在线观看 | 亚洲午夜理论影院| 久久久久久久精品吃奶| av.在线天堂| 看黄色毛片网站| 久久久久久国产a免费观看| 久久草成人影院| 又黄又爽又刺激的免费视频.| 99精品在免费线老司机午夜| 亚州av有码| 免费观看精品视频网站| 五月玫瑰六月丁香| 久久久精品大字幕| 久久人人爽人人爽人人片va| 又黄又爽又刺激的免费视频.| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 97人妻精品一区二区三区麻豆| aaaaa片日本免费| 欧美又色又爽又黄视频| xxxwww97欧美| 少妇人妻一区二区三区视频| 欧美极品一区二区三区四区| 最新中文字幕久久久久| 日本在线视频免费播放| 久久精品久久久久久噜噜老黄 | 日日摸夜夜添夜夜添av毛片 | 国产精品电影一区二区三区| 高清日韩中文字幕在线| 精品久久久久久久久久免费视频| 国产亚洲91精品色在线| 国产成人a区在线观看| 久久精品影院6| 美女被艹到高潮喷水动态| 女生性感内裤真人,穿戴方法视频| 联通29元200g的流量卡| 99在线人妻在线中文字幕| 色av中文字幕| 午夜福利欧美成人| 性欧美人与动物交配| 国产 一区精品| 日本黄色片子视频| 精华霜和精华液先用哪个| 一级a爱片免费观看的视频| 久久99热6这里只有精品| 久久这里只有精品中国| 日韩一区二区视频免费看| 欧美成人性av电影在线观看| 国产午夜福利久久久久久| 两人在一起打扑克的视频| 久久久成人免费电影| 99热精品在线国产| 一区二区三区四区激情视频 | 国产精品乱码一区二三区的特点| 在线观看午夜福利视频| 亚洲av中文av极速乱 | 亚洲熟妇熟女久久| 亚洲国产精品合色在线| 国内精品宾馆在线| 色在线成人网| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 亚洲人成网站在线播放欧美日韩| 久久久久性生活片| 久久久久久久久大av| 色综合婷婷激情| 丝袜美腿在线中文| 久久久久久九九精品二区国产| 老司机午夜福利在线观看视频| 美女高潮喷水抽搐中文字幕| 国产午夜精品论理片| 高清毛片免费观看视频网站| 2021天堂中文幕一二区在线观| 久久久久精品国产欧美久久久| 九九热线精品视视频播放| 国产精品98久久久久久宅男小说| 三级毛片av免费| 久久这里只有精品中国| 小蜜桃在线观看免费完整版高清| 最近最新免费中文字幕在线| 久久欧美精品欧美久久欧美| 又爽又黄a免费视频| 日韩强制内射视频| 亚洲人成网站在线播放欧美日韩| 国产亚洲欧美98| 精品乱码久久久久久99久播| 成人永久免费在线观看视频| 亚洲av免费高清在线观看| 国产 一区精品| 亚洲欧美日韩高清在线视频| 一夜夜www| 我的老师免费观看完整版| 一级av片app| 日本精品一区二区三区蜜桃| 国产精品乱码一区二三区的特点| av女优亚洲男人天堂| 在线观看一区二区三区| 真人做人爱边吃奶动态| 夜夜爽天天搞| 婷婷精品国产亚洲av在线| 久久久久精品国产欧美久久久| 最近视频中文字幕2019在线8| 亚洲精品一卡2卡三卡4卡5卡| 人妻少妇偷人精品九色| 国产伦精品一区二区三区视频9| 国产伦精品一区二区三区四那| 国产精品久久久久久av不卡| 久久久午夜欧美精品| 国产综合懂色| 搡老妇女老女人老熟妇| 身体一侧抽搐| 免费黄网站久久成人精品| 欧美一区二区国产精品久久精品| 又紧又爽又黄一区二区| 亚洲一区高清亚洲精品| 免费高清视频大片| 国产毛片a区久久久久| 熟女人妻精品中文字幕| 丰满的人妻完整版| 波多野结衣巨乳人妻| 久久午夜福利片| av天堂在线播放| 亚洲人成网站高清观看| 欧美激情在线99| 成人国产一区最新在线观看| 老师上课跳d突然被开到最大视频| 免费观看人在逋| 99久久久亚洲精品蜜臀av| 亚洲精品乱码久久久v下载方式| 嫩草影院精品99| 精品一区二区三区视频在线观看免费| 女的被弄到高潮叫床怎么办 | 色综合色国产| 韩国av在线不卡| 99热网站在线观看| 欧美日韩国产亚洲二区| 日韩亚洲欧美综合| 国产精品av视频在线免费观看| 91在线观看av| 成年人黄色毛片网站| 最新中文字幕久久久久| netflix在线观看网站| 久久午夜亚洲精品久久| 3wmmmm亚洲av在线观看| 亚洲av中文字字幕乱码综合| 一本久久中文字幕| 婷婷色综合大香蕉| 在线国产一区二区在线| 亚洲av免费高清在线观看| 国产主播在线观看一区二区| 久久久久免费精品人妻一区二区| 麻豆国产av国片精品| 黄色女人牲交| 久久久久久久久久成人| 乱系列少妇在线播放| 色吧在线观看| 日韩,欧美,国产一区二区三区 | 九九久久精品国产亚洲av麻豆| 国产伦在线观看视频一区| 成人性生交大片免费视频hd| 久久精品国产亚洲av香蕉五月| 特级一级黄色大片| 欧美xxxx性猛交bbbb| 深夜精品福利|