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

    HVSR法數(shù)據(jù)精細化處理中若干干擾和校正問題的解決方法研究

    2021-01-26 03:41:56周世昌王斌戰(zhàn)徐元璋苑益軍王勝侯
    工程地球物理學(xué)報 2021年1期
    關(guān)鍵詞:時窗干擾信號零點

    周世昌,王斌戰(zhàn),邱 波,徐元璋,劉 磊,苑益軍,王勝侯

    (1.湖北省地質(zhì)局 地球物理勘探大隊物探所,湖北 武漢 430056;2.中國地質(zhì)大學(xué) 地球物理與信息技術(shù)學(xué)院,北京 100083;3.中國地質(zhì)大學(xué) 教育部構(gòu)造與石油資源重點實驗室,湖北 武漢 430074)

    1 引 言

    天然地震領(lǐng)域的成熟地球物理方法,通常是解決大深度,廣區(qū)域的地質(zhì)問題,如果將其應(yīng)用到相對淺部的礦產(chǎn)勘探和工程勘察領(lǐng)域,可以解決許多能源、礦產(chǎn)、工程等生產(chǎn)問題。微動技術(shù)是利用天然產(chǎn)生的瑞雷面波或人工激發(fā)的瞬態(tài)瑞雷面波中不同頻率成分的穿透深度和傳播速度的不同,用其速度和深度(半波長)來對地下界面進行識別的地球物理勘探技術(shù)[1-6]。

    基于微動的H/V法是一種更加便捷的微動方法[7-12],又稱為三分諧振或HVSR方法,它最早由日本地震學(xué)家中村(Nakamura)1989年提出,它是一種估算地表層振動共振頻率和放大的技術(shù),其計算出微動信號的水平分量和垂直分量之比,典型的H/V譜比曲線具有一個明顯的峰值頻率fm。

    鑒于該方法劃分土石分界面的應(yīng)用效果良好,且布點靈活機動,故該方法在城市地質(zhì)勘察中應(yīng)用廣泛。然而城市中震動干擾因素眾多,如車輛震動干擾、高壓線的電磁干擾、儀器本身的零點漂移等,均讓采集的一手數(shù)據(jù)無法直接進行譜比計算,現(xiàn)今也少有文章全面且詳細地討論譜比計算時常見干擾因素的去除和畸變因素的校正,以及如何獲得優(yōu)質(zhì)的譜比曲線。

    本文將以武漢市江夏區(qū)法泗鎮(zhèn)巖溶塌陷區(qū)采集的三分量地震數(shù)據(jù)為例,詳細分析其中的干擾因素和畸變因素,并對如何進行去噪、校正、優(yōu)質(zhì)譜比曲線計算做細致介紹。

    2 瞬態(tài)干擾

    對采集的原始數(shù)據(jù)進行處理之前必須對原始信號進行信噪分析,從中剝離掉瞬態(tài)干擾信號。在城市道路旁采集時瞬態(tài)干擾往往源于行駛車輛和環(huán)境聲音,且振幅較大,剝離瞬態(tài)干擾后余留下的有效信號才能進行頻譜計算。本文標記瞬態(tài)干擾信號采用STA/LTA方法,傳統(tǒng)的STA/LTA算法在微震監(jiān)測和強震檢測中用于識別突變信號(劉晗等,2014;楊黎薇等,2017)[16,20],本文用該方法識別出突變信號(突變信號在城區(qū)施工時被認為是瞬態(tài)干擾信號)STA/LTA方法的理論原理如下:

    式(1)中,i為STA和LTA的計算時刻;NS的值決定了節(jié)選的信號樣點數(shù)量;Z(j)為j時刻垂直分量的振幅值;N(j)為j時刻北向分量的振幅值;E(j)為j時刻東向分量的振幅值,λ為“觸發(fā)閾值”。

    STA值計算的樣點數(shù)NS對應(yīng)捕捉干擾信號的時間窗,因此時窗越短,就對短周期的干擾信號捕捉越有效,LTA值是用于衡量時間窗內(nèi)的平均噪聲。STA/LTA就可以根據(jù)周圍環(huán)境噪聲程度自適應(yīng)地調(diào)整其對于某一類型干擾信號的敏感度。STA時間窗越短,越敏感,LTA時間窗越長,越敏感。對于局域干擾事件的捕捉,特別是從城市道路旁采集的數(shù)據(jù)中捕捉瞬態(tài)干擾,STA時間窗的典型值在1~10 s之間,LTA初始值一般在25 s左右,也可以適當(dāng)增大。

    “觸發(fā)閾值”λ越小,對于干擾信號越敏感,但也容易帶來誤拾??;相反地,閾值越大,被識別的干擾信號就越少。通過參數(shù)試驗發(fā)現(xiàn),本次采集的數(shù)據(jù)STA時窗采用5 s,LTA時窗采用25 s,“觸發(fā)閾值”λ選1.3,可有效識別出瞬態(tài)干擾事件(野外采集時發(fā)現(xiàn)這些干擾多由測線旁邊高架路上的行駛車輛產(chǎn)生)。如圖1中的“紅色凸起”即為按照STA/LTA法標記的干擾事件,在預(yù)處理環(huán)節(jié)可以按照這些標識剔除掉干擾信號,保留余下的有效信號。

    圖1 采用STA/LTA方法標記出原始數(shù)據(jù)中的瞬態(tài)干擾(紅色波峰)

    3 零點漂移

    由于儀器制造工藝等客觀原因影響,現(xiàn)今的三分量檢波器大多容易受到儀器低頻噪聲、環(huán)境背景信號、人為處理誤差及初始加速度等的影響,由地震加速度記錄積分得到的地震波位移時間曲線普遍會出現(xiàn)基線漂移現(xiàn)象(即零點漂移),造成振蕩起跳點并非位于0值。如果不進行零點漂移校正,會造成低頻信號的頻譜值失真,從而造成低頻信號的譜比值失真。

    如圖2所示的工區(qū)某道三分量信號,圖中可以看出Z分量、N分量、E分量均存在負的零點漂移(圖中綠色實線為零點基線,Z、N、E三個分量的振動曲線的中心基線均在0值往下)。用圖2中信號計算頻譜,結(jié)果如圖3所示,頻譜圖中0 Hz的振幅值非常突出。原始數(shù)據(jù)中的零漂幅值可以視為一個低頻分量,故造成頻譜中0 Hz存在顯著凸起“針狀尖峰”,而有效信號幅值相對零漂幅值(即“針狀尖峰”)較小,所以在頻譜圖中展示不夠明顯,故圖3中2~10 Hz左右的有效信號頻譜就顯得尤為“矮小”。

    圖2 原始三分量數(shù)據(jù)存在顯著的零點漂移(綠線為0值線)

    本文為了解決數(shù)據(jù)中存在的零點漂移,采用“分段標記重構(gòu)法”,即認為在局部時間段里數(shù)據(jù)的零點漂移值為固定常數(shù)。事實上儀器的零點漂移是緩慢進行的,如果儀器零點漂移劇烈,且幅值變化較大,說明儀器本身存在故障,相應(yīng)的采集數(shù)據(jù)也不可靠,必須舍棄掉。“分段標記重構(gòu)法”,就是抓住零點漂移的緩慢性,通過對局部信號進行傅里葉變換獲得該信號直流分量的值,并將該值用于這段信號的零點漂移校正。在實際處理步驟如下:

    1)設(shè)定好零點漂移值計算時間窗長度(例如:時窗長度1 min)。

    2)采用傅里葉變換計算時間窗內(nèi)信號的直流分量大小。

    3)按照固定步長滑動時間窗,計算下一個時間窗內(nèi)的直流分量大小(例如:滑動步長為1 s)。

    4)重復(fù)步驟“1)”,直到時間窗的尾部滑動到信號結(jié)尾為止。

    通過如上步驟可以按照固定時間間隔對整個地震信號進行零點漂移值的估算,再對估算的零點漂移值按照信號采樣率進行插值,即可獲得整段信號的每個采樣點的零點漂移校正值,用原始信號減去零點漂移校正值即完成了信號的零點漂移校正。圖4展示了采用“分段標記重構(gòu)法”對圖2中數(shù)據(jù)進行零點漂移值估算的結(jié)果,圖中紅色實線表示估算零點漂移線,從中可以看出紅色實線和振動曲線的零點漂移起伏規(guī)律十分吻合。

    圖5展示了圖4中數(shù)據(jù)按照估算零點漂移線進行零點漂移校正的結(jié)果,校正之后數(shù)據(jù)的振動基線已經(jīng)校正到0值位置。圖6展示了圖5中數(shù)據(jù)的頻譜結(jié)果,通過圖3和圖6對比可以看出零點漂移校正之后0 Hz位置的凸起“針狀尖峰”被去除,同時有效信號的頻譜被突顯出來。

    圖4 采用“分段標記重構(gòu)法”對原始數(shù)據(jù)進行零點漂移值估算(紅線為零漂線)

    圖5 進行零點漂移校正之后的數(shù)據(jù)

    圖6 進行零點漂移校正之后的頻譜

    4 諧波干擾

    諧波干擾是原始地震數(shù)據(jù)中一種常見的干擾波,它具有振幅固定、頻率單一的特點[21-25]。其主要有三個產(chǎn)生機制:①當(dāng)檢波器附近存在高壓輸電線時,檢波器會受到工業(yè)交流電產(chǎn)生的50 Hz電磁場干擾,以及150 Hz和250 Hz等頻率的伴隨電磁場干擾。②當(dāng)檢波器附近存在頻率固定的震源時(例如:停滯態(tài)汽車、勻速轉(zhuǎn)動的馬達、工作的發(fā)電機),檢波器會受到單一頻率的機械震動干擾。③當(dāng)檢波器內(nèi)部漏電時,檢波器也會記錄下諧波電流干擾。諧波干擾一般貫穿于整個地震記錄,其振幅往往明顯大于有效反射信號的振幅,從而會降低地震數(shù)據(jù)的信噪比,在計算譜比值時往往在諧波干擾頻率處產(chǎn)生畸變。因此,必須對地震數(shù)據(jù)中的諧波干擾波予以去除。

    圖7為數(shù)據(jù)中存在諧波干擾時的頻譜圖,從圖中可以看到在7.5 Hz左右存在一個顯著的“凸起頻率”,如果不對諧波干擾進行消除會使譜比曲線的相應(yīng)位置發(fā)生畸變。

    圖7 受諧波干擾影響的頻譜

    本文介紹一種“自適應(yīng)諧波干擾降噪技術(shù)”來去除該干擾。在時間域,原始地震信號可以看作有效信號和諧波干擾信號的疊加,如公式(2)所示:

    X=S+Y

    (2)

    (3)

    YT=(y1,y2,…,yi,…yk)

    公式(3)中XT為地震信號采樣序列,ST為有效信號采樣序列,YT為諧波干擾采樣序列,i為采樣點編號,k為總的采樣點個數(shù),T表示向量的轉(zhuǎn)置?,F(xiàn)需要在已知原始采樣序列XT的情況下,估算出干擾序列YT。諧波信號用正弦函數(shù)表示如下:

    yi=Lsin2πf(i+τ)Δt

    (4)

    式(4)中,L表示振幅;f表示頻率;Δt表示采樣時間間隔;τ表示初始相位參數(shù)。

    (5)

    由原始采樣序列估算出諧波干擾序列,必須建立合理的目標函數(shù):

    (6)

    其中,Q表示原始采樣序列減去諧波干擾序列后的振幅能量值。隨著構(gòu)造的諧波干擾序列的變化,Q也會變化,而使得Q取極小值的重構(gòu)序列就是實際的諧波干擾序列。將目標函數(shù)式(6)改寫為向量形式如下:

    要獲得Q的極小值,可以對Q關(guān)于各變量求導(dǎo)。假設(shè)頻率f和初始相位參數(shù)τ已知,對Q關(guān)于振幅L求導(dǎo),并使導(dǎo)數(shù)為0,可以求出頻率和初始相位一定時的最佳振幅,如公式(9)所示:

    (9)

    式(9)就是正弦函數(shù)法重構(gòu)諧波干擾時的振幅計算公式,如果確定了諧波干擾的頻率和初始相位,則可由該公式計算出相應(yīng)的振幅值。

    將式(9)代入式(7),目標函數(shù)可轉(zhuǎn)化為:

    式(10)和式(11)中,X是已知量,C是未知量。因此,使Q獲得極小值只需要R獲得極大值即可,R表達式的上方(XTC)2為原始地震序列和諧波干擾序列互相關(guān)的平方,R表達式的下方CTC為諧波干擾序列的自相關(guān)。所以,目標函數(shù)式(10)可以重新定義為原始地震序列和諧波干擾序列的歸一化互相關(guān):

    (12)

    使得Corr為極大值的f與τ就是重構(gòu)諧波干擾的最佳頻率和初始相位參數(shù)。再由式(9)計算出振幅值,根據(jù)式(4)可以重構(gòu)出諧波干擾。從原始地震信號中減去重構(gòu)的干擾信號進行去噪。

    圖8是對圖7中數(shù)據(jù)采用了“自適應(yīng)諧波干擾降噪技術(shù)”處理的結(jié)果,可以看出,7.5 Hz左右的“凸起頻率”被消除,而其他有效頻率成分的信號未有變化,即該方法在去除諧波干擾的同時,不存在傷害其它有效信號的“副作用”,在諧波干擾消除之后相應(yīng)頻率位置的譜比值畸變也被消除。

    圖8 采用“自適應(yīng)諧波干擾降噪技術(shù)”處理后的頻譜

    5 譜比曲線的優(yōu)化

    原始數(shù)據(jù)通過有針對性的預(yù)處理后便可開始頻譜和譜比計算,但是頻譜計算時窗該如何選擇?獲得的頻譜和譜比值存在‘鋸齒狀毛刺’該如何消除?這些問題解決不好,依然不能獲得品質(zhì)良好的譜比曲線。

    本文采用了滑動時窗計算多個譜比曲線進行疊加來改善譜比曲線品質(zhì),在計算時采用了歐洲SESAME提出的參數(shù)確定準則(Bard and SESAME-Team,2005):

    準則1 峰值頻率fm和時窗長度L應(yīng)該滿足:fm>10/L。

    準則2 時窗數(shù)量N、時窗長度L、峰值頻率fm應(yīng)滿足:L×N×fm>200。

    準則3 如果fm>0.5 Hz,當(dāng)0.5fm

    當(dāng)然在數(shù)據(jù)處理的初始情況下是不知道峰值頻率fm的,無法從一開始就獲得較為合理的時窗長度L,實際處理時L按照一定的間隔逐步擴大,計算每一個時窗的fm值,伴隨著時窗的逐步擴大,fm的值也趨于穩(wěn)定,這時再選用穩(wěn)定的fm值按照上述準則確定優(yōu)選時窗和步長。當(dāng)然也有更為簡單直接的時窗選擇方法(一般的三分量地震記錄的總時間長度會比時窗長度L大得多),直接選用較長的時窗也可估算fm的值。

    圖9為工區(qū)某道三分量地震信號分別按照10 s、30 s、60 s時窗計算出的N/Z譜和E/Z譜,通過該實驗可以確定在時窗長度大于10 s的情況下該道譜比曲線的峰值頻率穩(wěn)定在6 Hz左右,按照參數(shù)確定“準則1”可以確定本信號頻譜計算時窗長度大于5/3 s即可,最終選定為10 s作為該道的處理時窗。

    圖9 三分量地震信號按照10 s、30 s、60 s時窗計算出譜比曲線

    判斷“準則3”中某個頻率樣點fi(0.5fm

    第一步:用該頻率樣點fi全部時窗的譜比值計算標準差。

    第二步:按照“準則3”判斷標準差是否合適,如果合適進行疊加求平均。

    第三步:如果“第二步”中判斷標準差不合適,則舍棄掉離平均值最大的譜比值,再進入“第一步”,通過多次迭代即可獲得滿足要求的疊加結(jié)果。

    采用滑動時窗進行疊加計算得到的譜比曲線其連續(xù)性和平滑性會顯著提高,后續(xù)還可以選擇性地采用帕曾窗進一步進行平滑處理(視疊加效果而定)。帕曾窗平滑處理是按照一定的頻率帶寬求取帶寬內(nèi)振幅平均值作為帶寬中值頻率的振幅。

    圖10(a)展示了工區(qū)內(nèi)某道三分量數(shù)據(jù)采用單一時窗(節(jié)選了10 s長度)計算的譜比曲線,該譜比曲線存在顯著的“鋸齒狀”特征,圖10(b)是采用該時窗(10 s長度)按照1 s的滑動步長進行疊加計算,同時采用帕曾窗按照1 Hz帶寬進行曲線平滑的結(jié)果,可見“鋸齒狀”特征得到顯著改善,譜比曲線得到有效平滑。

    圖10 經(jīng)過平滑處理前后的譜比值曲線

    6 土石分界面和峰值頻率的對應(yīng)關(guān)系

    第四系土層和基巖分界面之間存在對應(yīng)關(guān)系,土石分界面深度h和H/V譜比曲線峰值頻率fm之間存在冪函數(shù)關(guān)系為(Seht and Wohlenberg,1999):

    (13)

    (14)

    lgh=lga+blgfm

    (15)

    如果令:lgh=y,lgfm=x,則式(13)可表達為y=lga+bx,收集武漢市法泗地區(qū)前期開展的H/V方法獲得的峰值頻率和工程鉆孔揭露的土石分界面深度關(guān)系如表1:

    表1 武漢市法泗地區(qū)鉆孔揭露土石分界面深度和H/V峰值頻率對照關(guān)系

    對表中l(wèi)gh和lgfm進行擬合如圖11所示,并對其進行直線擬合可解得b=-1.858 7(擬合直線的斜率),lga=2.641 6(擬合直線的截距),a=438.127

    圖11 峰值頻率和深度的擬合關(guān)系

    7 原始數(shù)據(jù)精細處理前后效果對比

    武漢法泗鎮(zhèn)長虹村、八壇村在2014年9月5日曾發(fā)生過大規(guī)模巖溶地陷,形成了9個大小不一的錐形深坑,后用黏土和碎石對深坑進行了回填?,F(xiàn)依托巖溶調(diào)查相關(guān)項目對塌陷回填區(qū)(塌陷坑位置已知)開展了HVSR法和其它物探方法試驗。下文將展示同一塌陷部位HVSR法和高密度電法的應(yīng)用效果。

    圖12為原始數(shù)據(jù)未進行精細化處理時獲得的譜比剖面(選用全時段信號直接進行譜比計算),圖13為原始數(shù)據(jù)進行精細化處理后獲得的譜比剖面。兩者對比發(fā)現(xiàn),精細化處理前的剖面低頻部分存在許多虛假異常,且譜比值的分布極不均勻,峰值頻率的連續(xù)性和均一性也不好(如果不是在已知塌陷坑開展工作,僅憑該剖面判斷,容易把小號40~50 m處也劃定為塌陷坑),而精細化處理后的剖面土石分界面(峰值頻率連線)的連續(xù)性和均一性明顯提高,巖溶塌陷坑的位置和邊界都更加清晰。

    圖12 塌陷坑回填位置的譜比法剖面(未精細化處理)

    圖13 塌陷坑回填位置的譜比法剖面(精細化處理后)

    圖14為塌陷部位高密度電法探測的結(jié)果,該方法對黏土層、砂層、基巖面分界線反映清晰,剖面中靠近大號部位的中間砂層中存在明顯的一處低阻異常,該異常為原塌陷坑所在部位,與地表塌陷位置一致。

    圖14 塌陷坑回填位置的高密度電法剖面

    通過圖12和圖13的橫向?qū)Ρ日f明:對HVSR方法原始數(shù)據(jù)有針對性地進行數(shù)據(jù)去噪、校正、譜比優(yōu)化計算是十分必要的。通過圖13和圖14兩種技術(shù)方法成果圖的對比,以及它們和鉆孔資料的對應(yīng)關(guān)系可以看出,塌陷區(qū)土石分界面在HVSR法剖面上表現(xiàn)為低比值,在高密度電法剖面上表現(xiàn)為低阻。同時HVSR剖面未發(fā)現(xiàn)兩層顯著峰值,這說明HVSR法對黏土層、砂層的區(qū)分不如高密度電法清晰(此條結(jié)論僅限于本工區(qū)),但是對基巖面的響應(yīng)良好(圖13中卓越頻率聯(lián)線起伏即反映了基巖面的起伏),這也證明了該方法確實在劃分土石分界面時具有優(yōu)勢。

    8 結(jié) 論

    1)野外采集得到的原始三分量地震數(shù)據(jù)中往往存在許多干擾和失真因素,如果采用原始數(shù)據(jù)直接進行頻譜和譜比計算無法得到優(yōu)質(zhì)可靠的譜比剖面,需要有針對性地對其進行預(yù)處理,為譜比計算提供品質(zhì)良好的振幅數(shù)據(jù)。

    2)針對原始三分量數(shù)據(jù)中普遍存在的零點漂移現(xiàn)象和諧波干擾現(xiàn)象,本文采用“分段標記重構(gòu)法”和“自適應(yīng)諧波干擾降噪技術(shù)”可以很好地解決相應(yīng)問題,為后續(xù)的頻譜計算和譜比計算提供品質(zhì)良好的振幅數(shù)據(jù)。

    3)在頻譜計算和譜比計算時為提高相應(yīng)曲線的可靠性,還需從初始信號中剝離出瞬態(tài)干擾信號,而采用STA/LTA方法可以有效標記出瞬態(tài)干擾信號,按照標記可以剝離掉這部分信號。

    4)為提高頻譜曲線和譜比曲線的平滑性可以采用滑動時窗進行疊加(相應(yīng)計算參數(shù)可遵循歐洲SESAME項目提出的準則),同時可以選用帕曾窗對單個頻譜曲線和譜比曲線進行平滑處理。

    5)在數(shù)據(jù)處理前需要對數(shù)據(jù)品質(zhì)進行綜合分析(通常從時間域和頻率域兩個層面進行分析),數(shù)據(jù)預(yù)處理時先進行諧波干擾去除(如果頻率域分析發(fā)現(xiàn)存在該干擾),接著進行零點漂移校正,然后進行瞬態(tài)干擾標記與剝離,最后才能進行實質(zhì)性的頻譜和譜比優(yōu)化計算。

    猜你喜歡
    時窗干擾信號零點
    GRAPES-GFS模式2 m溫度預(yù)報的最優(yōu)時窗滑動訂正方法
    氣象科技(2021年5期)2021-11-02 01:36:28
    一種基于改進時窗法的爆炸沖擊波檢測方法
    正弦采樣信號中單一脈沖干擾信號的快速剔除實踐方法
    2019年高考全國卷Ⅱ文科數(shù)學(xué)第21題的五種解法
    基于粒子群算法的光纖通信干擾信號定位方法
    一類Hamiltonian系統(tǒng)的Abelian積分的零點
    一道高考函數(shù)零點題的四變式
    淺析監(jiān)控干擾信號的優(yōu)化處置措施
    不同介入時窗和療程對高壓氧治療急性脊髓損傷的影響研究
    UWB搜救生物雷達分段時窗探測系統(tǒng)的實現(xiàn)
    午夜福利在线观看吧| 脱女人内裤的视频| 一边摸一边做爽爽视频免费| 国产真人三级小视频在线观看| 国产精品.久久久| 在线观看免费日韩欧美大片| 欧美一级毛片孕妇| 欧美 亚洲 国产 日韩一| 在线观看www视频免费| 日韩欧美一区视频在线观看| 日本a在线网址| 免费不卡黄色视频| 99久久99久久久精品蜜桃| 高清欧美精品videossex| 在线播放国产精品三级| 99精品欧美一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 国产av又大| 久久久精品区二区三区| 亚洲欧美日韩另类电影网站| 变态另类成人亚洲欧美熟女 | 亚洲精品粉嫩美女一区| 国产午夜精品久久久久久| 人人妻人人澡人人爽人人夜夜| 欧美激情极品国产一区二区三区| 老司机影院毛片| 国产成人精品久久二区二区91| 午夜福利在线观看吧| 69精品国产乱码久久久| 最新在线观看一区二区三区| 亚洲精品国产一区二区精华液| 丰满迷人的少妇在线观看| 一级毛片女人18水好多| 久久国产亚洲av麻豆专区| 看片在线看免费视频| www日本在线高清视频| 又紧又爽又黄一区二区| 国产在线观看jvid| av网站免费在线观看视频| 中文字幕色久视频| 亚洲自偷自拍图片 自拍| 国产亚洲精品第一综合不卡| 亚洲九九香蕉| 亚洲免费av在线视频| 精品少妇久久久久久888优播| 美女福利国产在线| 国产成人免费无遮挡视频| 黄色毛片三级朝国网站| 久久精品亚洲熟妇少妇任你| 午夜两性在线视频| 嫩草影视91久久| 久久精品熟女亚洲av麻豆精品| 咕卡用的链子| 中文字幕av电影在线播放| 亚洲欧美一区二区三区久久| 老司机亚洲免费影院| 麻豆成人av在线观看| 大型黄色视频在线免费观看| 日韩成人在线观看一区二区三区| 亚洲熟妇中文字幕五十中出 | 大码成人一级视频| 精品福利观看| 另类亚洲欧美激情| 丝袜人妻中文字幕| 极品人妻少妇av视频| 亚洲第一青青草原| 色婷婷久久久亚洲欧美| 亚洲va日本ⅴa欧美va伊人久久| 丝袜美足系列| av不卡在线播放| 午夜影院日韩av| 日韩欧美免费精品| √禁漫天堂资源中文www| 午夜激情av网站| 精品人妻在线不人妻| 日韩人妻精品一区2区三区| 欧美中文综合在线视频| 激情视频va一区二区三区| 亚洲第一青青草原| 无人区码免费观看不卡| 老司机深夜福利视频在线观看| 日本欧美视频一区| 亚洲人成电影观看| 日日夜夜操网爽| 一进一出抽搐gif免费好疼 | 国产精品久久久久久人妻精品电影| 亚洲成国产人片在线观看| av不卡在线播放| 国产淫语在线视频| 在线观看午夜福利视频| 精品第一国产精品| 国产日韩一区二区三区精品不卡| 黄色a级毛片大全视频| 亚洲欧美一区二区三区久久| 国产99白浆流出| 99国产极品粉嫩在线观看| 久久国产精品大桥未久av| 久久久久久久精品吃奶| 亚洲国产中文字幕在线视频| 一二三四在线观看免费中文在| av线在线观看网站| 免费在线观看完整版高清| 亚洲av欧美aⅴ国产| 校园春色视频在线观看| 丰满饥渴人妻一区二区三| 美女高潮喷水抽搐中文字幕| 麻豆国产av国片精品| 国产精品电影一区二区三区 | 久久久国产成人精品二区 | 国产成+人综合+亚洲专区| 变态另类成人亚洲欧美熟女 | 亚洲精品乱久久久久久| 亚洲精品久久午夜乱码| a级毛片黄视频| 亚洲 国产 在线| 少妇被粗大的猛进出69影院| 免费高清在线观看日韩| 黄片大片在线免费观看| 首页视频小说图片口味搜索| 欧美最黄视频在线播放免费 | 中文字幕av电影在线播放| 久久国产精品影院| 精品亚洲成国产av| 亚洲av成人av| 免费看十八禁软件| 精品午夜福利视频在线观看一区| 久久青草综合色| 亚洲欧美精品综合一区二区三区| 欧美精品亚洲一区二区| 一边摸一边抽搐一进一出视频| 日韩欧美国产一区二区入口| 国产精品免费一区二区三区在线 | 国产精品九九99| 国产成人精品无人区| 日韩一卡2卡3卡4卡2021年| 精品第一国产精品| 成年版毛片免费区| 久久久精品区二区三区| 久久久精品国产亚洲av高清涩受| 黑人猛操日本美女一级片| 欧美日韩瑟瑟在线播放| 新久久久久国产一级毛片| 国产午夜精品久久久久久| 嫁个100分男人电影在线观看| 在线观看日韩欧美| 亚洲第一欧美日韩一区二区三区| svipshipincom国产片| 人人妻人人添人人爽欧美一区卜| 欧美乱色亚洲激情| 侵犯人妻中文字幕一二三四区| 曰老女人黄片| 国产成人一区二区三区免费视频网站| 成人三级做爰电影| 色综合婷婷激情| 国产99白浆流出| 精品人妻在线不人妻| 欧美日韩亚洲综合一区二区三区_| 最近最新免费中文字幕在线| 亚洲欧美激情综合另类| 久久精品亚洲熟妇少妇任你| 深夜精品福利| 一边摸一边抽搐一进一小说 | 人人妻人人澡人人爽人人夜夜| 99国产精品一区二区蜜桃av | 国产乱人伦免费视频| 中文字幕制服av| 精品久久久久久,| 黄片小视频在线播放| 日韩视频一区二区在线观看| 69精品国产乱码久久久| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 亚洲国产毛片av蜜桃av| 大陆偷拍与自拍| 欧美色视频一区免费| 国产高清国产精品国产三级| 热99国产精品久久久久久7| 99久久人妻综合| 啦啦啦 在线观看视频| 91国产中文字幕| 男女高潮啪啪啪动态图| 色婷婷久久久亚洲欧美| 久久国产乱子伦精品免费另类| 啦啦啦 在线观看视频| 最新在线观看一区二区三区| 精品一区二区三区视频在线观看免费 | 国产亚洲精品久久久久久毛片 | 手机成人av网站| 国产精品二区激情视频| 老司机福利观看| 亚洲精品中文字幕在线视频| 久久中文字幕一级| 亚洲色图综合在线观看| 久久国产乱子伦精品免费另类| 免费看十八禁软件| 91九色精品人成在线观看| 另类亚洲欧美激情| 少妇的丰满在线观看| 欧美成人午夜精品| 俄罗斯特黄特色一大片| 欧美最黄视频在线播放免费 | 亚洲男人天堂网一区| 精品一区二区三区四区五区乱码| 少妇 在线观看| 成人亚洲精品一区在线观看| av一本久久久久| 亚洲专区中文字幕在线| 这个男人来自地球电影免费观看| 国产成+人综合+亚洲专区| 亚洲,欧美精品.| 99久久精品国产亚洲精品| 国产精品99久久99久久久不卡| 亚洲熟妇熟女久久| 久久ye,这里只有精品| 久久久久久久国产电影| 国产精品久久久av美女十八| 精品第一国产精品| 久热爱精品视频在线9| 亚洲一区二区三区不卡视频| 精品视频人人做人人爽| 日韩大码丰满熟妇| 欧美最黄视频在线播放免费 | 亚洲欧美激情综合另类| 中文字幕av电影在线播放| 国产精品综合久久久久久久免费 | 国产成人啪精品午夜网站| 老司机深夜福利视频在线观看| 精品亚洲成a人片在线观看| 国产一区二区三区在线臀色熟女 | 久久久久国产一级毛片高清牌| 天天躁日日躁夜夜躁夜夜| 新久久久久国产一级毛片| 中文字幕人妻丝袜制服| 老司机福利观看| 建设人人有责人人尽责人人享有的| 大香蕉久久网| 欧美日韩亚洲国产一区二区在线观看 | 一本大道久久a久久精品| 日本精品一区二区三区蜜桃| 法律面前人人平等表现在哪些方面| 男人舔女人的私密视频| 中国美女看黄片| 色综合婷婷激情| 亚洲精品一二三| 久久精品国产99精品国产亚洲性色 | 亚洲熟女精品中文字幕| 黄色成人免费大全| 国产成人精品无人区| 免费一级毛片在线播放高清视频 | 亚洲美女黄片视频| av福利片在线| 999久久久精品免费观看国产| 性色av乱码一区二区三区2| 免费观看精品视频网站| 亚洲成人国产一区在线观看| 高清av免费在线| 日韩免费av在线播放| 午夜福利在线观看吧| 成人永久免费在线观看视频| 成年人免费黄色播放视频| 国产精品 欧美亚洲| av线在线观看网站| 99热国产这里只有精品6| 九色亚洲精品在线播放| 一夜夜www| 啦啦啦视频在线资源免费观看| 一区二区日韩欧美中文字幕| 国产成人欧美在线观看 | 亚洲第一欧美日韩一区二区三区| 国产在线一区二区三区精| 亚洲五月婷婷丁香| 女同久久另类99精品国产91| 国产成人精品无人区| 亚洲午夜理论影院| 亚洲av日韩精品久久久久久密| 欧美精品av麻豆av| 80岁老熟妇乱子伦牲交| 性少妇av在线| 大型黄色视频在线免费观看| 国产欧美日韩一区二区精品| 日本黄色视频三级网站网址 | 丝袜在线中文字幕| 十八禁人妻一区二区| 国产精品av久久久久免费| 9色porny在线观看| 1024香蕉在线观看| 黄色丝袜av网址大全| 色在线成人网| 日韩有码中文字幕| 一进一出抽搐gif免费好疼 | av国产精品久久久久影院| 一区二区日韩欧美中文字幕| av网站在线播放免费| 午夜免费成人在线视频| 国产麻豆69| 在线观看www视频免费| 老鸭窝网址在线观看| 美国免费a级毛片| 久久久久精品人妻al黑| 丰满人妻熟妇乱又伦精品不卡| 午夜福利视频在线观看免费| 亚洲aⅴ乱码一区二区在线播放 | 一二三四社区在线视频社区8| 国产成人欧美在线观看 | 久久香蕉国产精品| 狂野欧美激情性xxxx| a级毛片黄视频| 一本一本久久a久久精品综合妖精| 视频区图区小说| 欧美性长视频在线观看| 国产精品国产av在线观看| 如日韩欧美国产精品一区二区三区| 一区二区三区国产精品乱码| 无遮挡黄片免费观看| 欧美精品人与动牲交sv欧美| 国产成人av激情在线播放| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 午夜免费鲁丝| 十八禁高潮呻吟视频| 麻豆av在线久日| 亚洲国产精品合色在线| 性色av乱码一区二区三区2| 色播在线永久视频| 岛国毛片在线播放| 国精品久久久久久国模美| 欧美老熟妇乱子伦牲交| 91大片在线观看| 80岁老熟妇乱子伦牲交| 黄色丝袜av网址大全| 精品一品国产午夜福利视频| 免费在线观看亚洲国产| 在线观看66精品国产| 18禁美女被吸乳视频| 亚洲av片天天在线观看| 国产av一区二区精品久久| 在线观看免费午夜福利视频| 久久久久久免费高清国产稀缺| 视频在线观看一区二区三区| 国产亚洲精品久久久久久毛片 | 久久婷婷成人综合色麻豆| 视频在线观看一区二区三区| 欧美成人午夜精品| 国产精品成人在线| 国产亚洲一区二区精品| 后天国语完整版免费观看| 日韩中文字幕欧美一区二区| 两性夫妻黄色片| 丰满饥渴人妻一区二区三| 日韩大码丰满熟妇| 亚洲精品在线美女| 久久精品亚洲av国产电影网| 国产精品久久久久久精品古装| 丝袜人妻中文字幕| 精品一区二区三区视频在线观看免费 | 99国产综合亚洲精品| 久久青草综合色| 一二三四社区在线视频社区8| 黄片播放在线免费| 黄片小视频在线播放| 久久午夜综合久久蜜桃| 欧美黑人欧美精品刺激| 99精品在免费线老司机午夜| 国产精品.久久久| 免费一级毛片在线播放高清视频 | 91麻豆精品激情在线观看国产 | 一边摸一边做爽爽视频免费| 青草久久国产| 午夜视频精品福利| 国产精品美女特级片免费视频播放器 | 日韩一卡2卡3卡4卡2021年| 咕卡用的链子| 欧美在线黄色| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 亚洲精品在线观看二区| 又大又爽又粗| 一本综合久久免费| 一进一出好大好爽视频| 人人妻人人澡人人看| 久热爱精品视频在线9| 国产精品自产拍在线观看55亚洲 | 精品一品国产午夜福利视频| 女警被强在线播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲人成77777在线视频| 国产高清videossex| 精品人妻1区二区| 高清毛片免费观看视频网站 | 50天的宝宝边吃奶边哭怎么回事| 国产xxxxx性猛交| 欧美日韩福利视频一区二区| 水蜜桃什么品种好| 国产有黄有色有爽视频| 精品高清国产在线一区| 日韩欧美免费精品| 手机成人av网站| 亚洲五月天丁香| 黄色成人免费大全| 欧美成狂野欧美在线观看| 美女视频免费永久观看网站| 午夜老司机福利片| 无遮挡黄片免费观看| 在线av久久热| 日韩精品免费视频一区二区三区| 国内久久婷婷六月综合欲色啪| 欧美黑人精品巨大| 9色porny在线观看| 女人久久www免费人成看片| 中文字幕高清在线视频| 久久人妻熟女aⅴ| 国产精品国产高清国产av | 中文字幕制服av| 美女 人体艺术 gogo| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色 | 国产精品 欧美亚洲| 一进一出好大好爽视频| 久久青草综合色| 久久久国产精品麻豆| 精品久久久精品久久久| 亚洲av欧美aⅴ国产| 亚洲精品成人av观看孕妇| 两性夫妻黄色片| 中文字幕制服av| 很黄的视频免费| 国产野战对白在线观看| 人妻久久中文字幕网| 国产免费av片在线观看野外av| 久久国产精品人妻蜜桃| 五月开心婷婷网| 亚洲欧美一区二区三区久久| 欧美在线一区亚洲| 丰满饥渴人妻一区二区三| 久久香蕉激情| 91成年电影在线观看| 最新美女视频免费是黄的| 国产不卡av网站在线观看| 人成视频在线观看免费观看| 麻豆成人av在线观看| 如日韩欧美国产精品一区二区三区| 国产欧美日韩一区二区三区在线| 免费不卡黄色视频| 黑人猛操日本美女一级片| 欧美精品高潮呻吟av久久| 可以免费在线观看a视频的电影网站| 国产一区在线观看成人免费| av天堂久久9| 19禁男女啪啪无遮挡网站| 亚洲成人手机| cao死你这个sao货| 黄网站色视频无遮挡免费观看| 一级作爱视频免费观看| 日韩成人在线观看一区二区三区| 高清在线国产一区| 亚洲精品粉嫩美女一区| 欧美 亚洲 国产 日韩一| 久久久久久久久久久久大奶| 婷婷成人精品国产| 亚洲免费av在线视频| 一区福利在线观看| 亚洲精品久久午夜乱码| 欧美日韩亚洲综合一区二区三区_| 美女国产高潮福利片在线看| 欧美精品一区二区免费开放| av天堂在线播放| 国产亚洲欧美精品永久| 亚洲,欧美精品.| 在线观看免费高清a一片| 免费女性裸体啪啪无遮挡网站| 一二三四在线观看免费中文在| 国产色视频综合| 日韩欧美国产一区二区入口| 国产亚洲欧美精品永久| 99国产精品免费福利视频| 男男h啪啪无遮挡| 久久精品国产亚洲av高清一级| 国产精品亚洲一级av第二区| 99久久国产精品久久久| 亚洲第一av免费看| 国产精品一区二区免费欧美| 大陆偷拍与自拍| 9热在线视频观看99| 欧美成狂野欧美在线观看| 乱人伦中国视频| 国产男女内射视频| 亚洲男人天堂网一区| 亚洲国产欧美网| 免费人成视频x8x8入口观看| 妹子高潮喷水视频| 国产成人免费观看mmmm| 欧美午夜高清在线| 亚洲成人手机| 国产精品 国内视频| 欧美 日韩 精品 国产| 国产精品亚洲av一区麻豆| 国产精品国产av在线观看| 啦啦啦视频在线资源免费观看| 黄频高清免费视频| 黄网站色视频无遮挡免费观看| 精品国产一区二区久久| 欧美黑人精品巨大| 国产成人免费无遮挡视频| 男女免费视频国产| 婷婷成人精品国产| 精品久久蜜臀av无| 亚洲久久久国产精品| 黑人操中国人逼视频| 久久亚洲真实| 男人舔女人的私密视频| 午夜久久久在线观看| 99re6热这里在线精品视频| av电影中文网址| 日韩大码丰满熟妇| 女人高潮潮喷娇喘18禁视频| 欧美精品av麻豆av| 老司机影院毛片| 九色亚洲精品在线播放| 女警被强在线播放| 国产高清激情床上av| 看免费av毛片| 久久精品国产亚洲av高清一级| 精品一区二区三区视频在线观看免费 | 欧美精品人与动牲交sv欧美| 国产精品综合久久久久久久免费 | 久久99一区二区三区| 丰满饥渴人妻一区二区三| e午夜精品久久久久久久| 精品人妻在线不人妻| 在线观看免费高清a一片| 人成视频在线观看免费观看| 日韩免费av在线播放| 亚洲avbb在线观看| 在线免费观看的www视频| 欧美日韩国产mv在线观看视频| 亚洲精品自拍成人| 国产深夜福利视频在线观看| 色尼玛亚洲综合影院| 久久精品人人爽人人爽视色| 一级a爱视频在线免费观看| 午夜老司机福利片| 国产欧美日韩一区二区三| 王馨瑶露胸无遮挡在线观看| 一级,二级,三级黄色视频| ponron亚洲| 婷婷丁香在线五月| 老汉色av国产亚洲站长工具| 亚洲一区中文字幕在线| 女警被强在线播放| 日韩欧美一区二区三区在线观看 | 中文字幕人妻丝袜一区二区| 欧美性长视频在线观看| 久久久国产成人免费| 午夜福利欧美成人| 亚洲精品av麻豆狂野| 人成视频在线观看免费观看| 好男人电影高清在线观看| 男女高潮啪啪啪动态图| 大陆偷拍与自拍| 国产亚洲一区二区精品| 久久草成人影院| 捣出白浆h1v1| a级毛片黄视频| 国产国语露脸激情在线看| 美女午夜性视频免费| 最近最新中文字幕大全电影3 | 成人18禁在线播放| 黄色视频,在线免费观看| 精品国产亚洲在线| 久久久久国产一级毛片高清牌| 欧美日本中文国产一区发布| www日本在线高清视频| 久99久视频精品免费| 久久精品国产亚洲av香蕉五月 | 中亚洲国语对白在线视频| 国产在线一区二区三区精| 国产aⅴ精品一区二区三区波| 91av网站免费观看| 12—13女人毛片做爰片一| √禁漫天堂资源中文www| 视频区欧美日本亚洲| 国产视频一区二区在线看| 久久久国产成人精品二区 | 亚洲人成电影免费在线| 黄色女人牲交| 天天躁日日躁夜夜躁夜夜| 国产精品成人在线| 久久狼人影院| 精品卡一卡二卡四卡免费| 久久精品亚洲av国产电影网| 日本欧美视频一区| 久久国产亚洲av麻豆专区| 国产黄色免费在线视频| 欧美精品高潮呻吟av久久| 丁香六月欧美| 18禁裸乳无遮挡免费网站照片 | 国产国语露脸激情在线看| 午夜免费观看网址| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看黄色视频的| 午夜免费观看网址| 在线十欧美十亚洲十日本专区| 亚洲人成电影免费在线| 国产不卡av网站在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品一区二区三区四区五区乱码| 精品乱码久久久久久99久播| 岛国在线观看网站| 免费在线观看亚洲国产| www.999成人在线观看| 激情在线观看视频在线高清 | 久久久精品免费免费高清| 午夜福利影视在线免费观看| 别揉我奶头~嗯~啊~动态视频| 一进一出好大好爽视频|