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

    極化格雷互補碼探地雷達探測地下管線分布*

    2022-10-22 03:36:08張寧李靜霞徐航趙明瑞劉麗王冰潔
    電子器件 2022年4期
    關鍵詞:旁瓣極化格雷

    張寧李靜霞*徐航趙明瑞劉麗王冰潔

    (1.太原理工大學新型傳感器與智能控制教育部和山西省重點實驗室,山西 太原 030024;2.太原理工大學物理與光電工程學院,山西 太原 030024)

    城市地下管線的探測與繪制對于促進城建發(fā)展和保障民生至關重要[1-2]。對于有限的城市地下空間,管線大都以平行且走向不同的方式鋪設。在此種鋪設方式下,準確探測地下管線的位置與分布一直是管線探測的重點與難點。探地雷達作為探測管線的一種重要手段,其是利用電磁波在地下傳輸遇到目標時發(fā)生反射的原理進行探測的,具有快速、安全、簡單、低成本和無損檢測等特點[3]。

    現(xiàn)有探地雷達主要采用調(diào)頻連續(xù)波信號[4]、脈沖信號[5]、步進頻率信號[6]和隨機信號[7]作為探測信號。相較于前三種信號,基于隨機信號的探地雷達受益于良好的隨機性,可實現(xiàn)強抗電磁干擾測量[8]和無模糊探測[9]。隨機信號具體包括噪聲信號、混沌信號和偽隨機信號。其中,偽隨機序列因產(chǎn)生方法簡單且易集成,具有良好的應用前景。然而,基于單一的偽隨機序列會產(chǎn)生較高的旁瓣,導致強目標的旁瓣會將弱目標的主瓣覆蓋,從而增大漏警概率。故為減少旁瓣,應當尋找一種低旁瓣甚至零旁瓣的發(fā)射信號。格雷互補碼突破了傳統(tǒng)的單一序列,基于兩列信號自相關函數(shù)之和實現(xiàn)旁瓣抵消,從而實現(xiàn)弱目標增強。Alejos A V 等人[10-11]利用格雷互補碼作為探測信號,結(jié)合單邊帶傳輸對現(xiàn)有噪聲雷達系統(tǒng)進行改進,實現(xiàn)更好的動態(tài)范圍檢測和距離估計。作者所在課題組[12]采用格雷互補碼作為探地雷達的發(fā)射信號,實現(xiàn)了地下管線的高信噪比探測,進一步證明了格雷互補碼有助于提高地下目標的探測深度與精度。

    上述常規(guī)探地雷達通常采用單一極化方式對地下目標進行探測,接收到的回波數(shù)據(jù)中僅攜帶部分信息,難以實現(xiàn)目標的準確識別。已有研究發(fā)現(xiàn),電磁波的極化具有方向性,可用來確定目標的方向、傾角和形狀等屬性。此外,電磁波的不同極化方式對目標信號的濾波、增強和抗干擾等方面具有優(yōu)勢并且可以獲得穩(wěn)健的檢測能力,提高物體的分辨能力。因此,探地雷達的極化屬性被廣泛關注并逐步從單一極化測量轉(zhuǎn)到多極化乃至混合極化測量。Georgios P T 等人[13]通過不同的天線極化方式獲得多極化數(shù)據(jù),證明了多極化測量可以準確識別多方向性目標體的位置信息。吉林大學曾昭發(fā)課題組[14]提出多極化多輸入多輸出(Multiinput-Multioutput,MIMO)步進頻率探地雷達系統(tǒng),克服了目標截面和天線輻射的影響,有效提高了目標檢測精度。吉林大學馮晅課題組[15]將極化步進頻探地雷達系統(tǒng)與極化熵相結(jié)合,對地下不同特征目標進行識別。Shao Z L 等[16]提出了時域激發(fā)極化法,達到對地下管線的重建與定位的目的。Boniger U 等[17]將極化屬性分析法與主成分分析法相結(jié)合,實現(xiàn)了對地下管線的三維重建。以上方法均證明了電磁波的極化屬性在精確探測地下目標體上的潛力。盡管如此,現(xiàn)有方法對于走向分布復雜管線探測中存在的目標成像模糊以及虛假目標較多的問題仍未得到良好的解決且無實驗驗證。

    本文提出以格雷互補碼作為探測信號并結(jié)合極化組件獲取極化數(shù)據(jù)。通過將一組格雷互補碼的相關結(jié)果相加,實現(xiàn)旁瓣疊加為零,峰值信號加倍。同時結(jié)合地下管線的極化信息,實現(xiàn)地下管線的多極化測量,從而達到對地下管線的定位與分布成像。

    1 極化格雷互補碼探地雷達探測原理

    1.1 格雷互補碼的產(chǎn)生及特性分析

    假設一組序列包含A=a(n),B=b(n),碼長為N,非周期自相關函數(shù)分別為:

    當除零位置點外,其自相關函數(shù)的和為零,即滿足式(3)

    則稱序列A和B是一對格雷互補序列。其中,單列碼的旁瓣值約占峰值的10%,而通過將一對格雷互補碼的相關函數(shù)進行相加,可以實現(xiàn)完全消除旁瓣,且峰值增大一倍,是一種理想的雷達信號。

    進一步地,采用公式(4)的方式對碼長為N的A碼和B碼進行K次迭代,則可以產(chǎn)生N×2K長度的格雷互補碼序列。

    根據(jù)仿真和實驗環(huán)境,本文選擇碼長和頻率分別為1 024 和1 GHz 的一對格雷互補碼序列作為探測信號。圖1(a)、(b)分別為A碼和B碼的時序,其脈沖幅值為0.1 V,最窄脈沖寬度為1 ns。圖1(c)、(d)分別代表A碼和B碼的功率譜,帶寬為1 GHz。格雷互補碼的自相關曲線如圖1(e)所示,由自相關曲線可以看出,疊加后的自相關曲線峰值擴大一倍,旁瓣完全抵消。

    圖1 格雷互補碼特性

    1.2 極化格雷互補碼探地雷達測量原理

    基于格雷互補碼的探地雷達測量原理如下:雷達發(fā)射信號被功分成兩路,分別作為參考信號和探測信號。參考信號直接送回至接收機,探測信號則經(jīng)過雷達系統(tǒng)發(fā)射,再經(jīng)目標反射后返回接收機中。因為信號經(jīng)過目標體反射會產(chǎn)生時間延遲,因此通過將回波信號與參考信號做互相關可以計算出時延,進而求出管線的深度信息。

    格雷互補序列A=a(n)和B=b(n)的互相關表達式分別如式(5)和式(6)所示:

    上式中,?為相關算子,τ表示回波信號相對于參考信號的時延,k為相關系數(shù),VAcorr(t)、VAref(t)、VAech(t)、VBcorr(t)、VBref(t)和VBech(t)分別代表A碼和B碼的相關信號、參考信號以及回波信號。采用相關法求得A碼、B碼的互相關信息后,將A碼與B碼的互相關函數(shù)相加,可使二者時間延遲τ所對應的峰值加倍并消除旁瓣。格雷互補序列A和B互相關函數(shù)相加公式為:

    進一步地,沿測線(即雷達掃描軌跡)進行測量,可以獲得二維成像結(jié)果。探測時,采用天線組合的方式,獲取不同極化下的數(shù)據(jù)。不同天線的極化方式如圖2 所示,其中,T和R分別代表發(fā)射天線和接收天線。當發(fā)射方向為X,接收方向為X 時,代表天線的XX 極化方式,如圖2(a)所示。同理可依次獲得天線YY、XY、YX 極化方式,如圖2(b)~(d)所示。

    圖2 不同天線極化方式

    2 極化格雷互補碼探地雷達仿真分析

    2.1 一維軌線仿真分析

    仿真基于時域有限差分法(Finite-Difference-Time-Domain method,F(xiàn)DTD)進行,對埋設于2 m×1 m×2 m(長×寬×高)干沙中的金屬管線進行探測。金屬管線直徑、埋深、介電常數(shù)和磁導率分別為10 cm、0.5 m、300 m·N/C 和200 H/m~400 H/m,干沙的介電常數(shù)和磁導率分別為3 m·N/C 和1 H/m。信號源采用1 GHz 的格雷互補碼信號,天線緊貼地表掃描,等步長2 cm 沿X 軸正方向軌線進行測量。

    圖3 顯示了埋設單根金屬管線模型。本仿真中,將天線設置在管線正上方并旋轉(zhuǎn)360°,每5°測量一次。圖4(a)、(b)顯示了單根金屬管線的全角度響應圖。

    圖3 單根管線仿真模型

    圖4 單根管線全角度響應圖

    從全角度響應圖的分布可以看出,采用共極化模式時,隨著管線與天線之間角度的變化,散射能量變化較為均勻,因此在各個角度上均能夠探測到目標響應。而當采用交叉極化模式時,散射能量隨角度變化很大,當角度為45°或其倍數(shù)時散射能量達到最大。因此在探測不同走向和方向管線時可能會產(chǎn)生目標丟失,弱目標探測不明顯的現(xiàn)象。由此說明,當探測地下管線等線性目標體時,共極化的探測方式優(yōu)于交叉極化探測方式。

    進一步地,研究不同走向分布的雙根金屬管線的極化響應,模型如圖5 所示。將其中一根管線設置為與X 軸夾角為90°,另一根管線與X 軸為不同交角,設為α,分別采用XX 極化、XY 極化、YX 極化和YY 極化四種方式對其進行探測,結(jié)果如圖6所示。

    其中,圖6(a1)~(a4)、(b1)~(b4)、(c1)~(c4)、(d1)~(d4)分別代表另一根金屬管線與X 軸的夾角α為0°、30°、45°、60°的B-Scan 響應圖。從B-Scan 響應圖可以得出以下結(jié)論:對于材料相同、走向分布相同的管線,XX 極化和YY 極化中獲得的信號響應強于XY 極化和YX 極化的信號響應且YY 極化的信號響應強于XX 極化。而對于同一種極化探測方式,不同走向分布管線的雷達信號強度不同。例如采用XX 極化方式進行探測時,管線在0°、30°、45°、60°延伸時管線的信號響應均強于90°,其中0°延伸時信號響應最強。信號響應強度隨管線與X 軸的交角變化而呈現(xiàn)規(guī)律性變化。這說明當探測具有方向性的線性目標時,交叉極化會出現(xiàn)數(shù)據(jù)丟失的情況,考慮到散射效應,細長的圓柱線性目標如地下管線使用相互平行的共極化方式測量會產(chǎn)生更好的效果。因此,在進行管線分布探測時,應當優(yōu)先考慮共極化方式探測。以下仿真和實驗中,均采用XX 極化和YY 極化進行探測。

    圖6 金屬管線與X 軸夾角分別為0°、30°、45°、60°的B-Scan 響應圖

    2.2 二維軌線仿真分析

    基于以上一維軌線仿真的成像分析,將天線掃描范圍擴展為二維軌線。仿真模型如圖7 所示,分別對L 型金屬和非金屬管線進行探測。金屬管線與周圍干沙的基本參數(shù)與一維軌線仿真中設置相同。非金屬管線的直徑、埋深、介電常數(shù)和磁導率分別為10 cm、1 m、3.3 m·N/C 和1 H/m。天線緊貼地表掃描,天線之間的距離為2 cm,并以2 cm 的移動距離沿X軸正方向軌線運動,共獲得20 道剖面數(shù)據(jù),從每一道剖面數(shù)據(jù)的B-Scan 中可以確定管線的位置信息。將所獲得剖面數(shù)據(jù)進行三維數(shù)據(jù)插值即可得到三維多極化數(shù)據(jù),再將三維多極化數(shù)據(jù)沿B-Scan 頂點做水平切片,即可得到地下管線的分布成像。

    圖7 不同走向管線仿真模型2

    成像結(jié)果如圖8 所示,圖8(a1)、(a2)和(b1)、(b2)分別代表L 形金屬管線和L 形非金屬管線的XX、YY極化方式,圖8(a3)、(a4)和(b3)、(b4)分別代表L 形金屬管線和L 形非金屬管線的極化合成圖與三維顯示圖。從圖中可以看到,不同極化方式對不同方向管線的散射能量不同,如圖8(a1)、(a2)、(b1)和(b2)所示。在管線相交處散射強度大于任何一個地方的散射強度。非金屬管線與金屬管線類似,但由于非金屬管線的介電常數(shù)遠小于金屬管線且與干沙的介電常數(shù)極為相近,所以整體散射能量較低。以上仿真結(jié)果表明:以格雷互補碼作為發(fā)射信號并結(jié)合天線的不同極化方式可以對管線的分布清晰成像。

    圖8 金屬管線與非金屬管線不同極化顯示圖

    3 極化格雷互補碼探地雷達實驗分析

    仿真結(jié)果顯示了極化格雷互補碼探地雷達對地下管線走向分布精確成像的能力。在此基礎上,本文采用極化格雷互補碼探地雷達對地下管線進行了實驗研究。圖9 所示為實驗裝置圖。在發(fā)射端,信號源產(chǎn)生的格雷互補碼信號先經(jīng)過放大器1 放大并與信號發(fā)生器產(chǎn)生的3.4 GHz 正弦信號混頻后,再經(jīng)放大器2 放大,最后向地下發(fā)射信號對管線進行探測。在接收端,接收信號先經(jīng)過放大器3 放大,再經(jīng)過功分器2 分成兩路,隨后分別經(jīng)過放大器4 和5 放大,下變頻進行IQ 分解,最后經(jīng)過I 路和Q 路的低噪放大器放大后,示波器接收回波數(shù)據(jù)。其中,IQ 分解是為了獲得更大的動態(tài)范圍、提高測量精度和抗擊共模干擾。

    圖9 實驗裝置圖

    實驗中,T 型非金屬管線埋設于2.0 m×1.2 m×0.8 m(長×寬×高)的干沙箱中,實驗場景如圖10 所示。其中,干沙的介電常數(shù)為4 m~6 m N/C,T 型非金屬管線的埋深為38 cm,直徑為20 cm,介電常數(shù)為9 m N/C。天線采用XX極化與YY極化兩種極化方式對地下目標進行探測,發(fā)射和接收天線的間距為2 cm。

    圖10 實驗場景圖

    實驗結(jié)果如圖11(a1)~(a3)所示,分別代表XX極化、YY極化以及極化合成圖。由圖可知,目標響應與管線實際位置一致。當采用XX極化時,橫向位置所獲取的極化信息多于縱向位置;當采取YY極化時,縱向位置所獲得的極化信息多于橫向位置。當將多極化信息結(jié)合時,可以反映出管線的分布特性,如圖11(a3)所示,目標響應與周圍介質(zhì)區(qū)分明顯,實驗結(jié)果與仿真結(jié)果一致。

    為了分析極化格雷互補碼探地雷達的探測性能,將其與極化步進頻率雷達進行了比較。實驗中,將矢量網(wǎng)絡分析儀用作步進頻率雷達,信號頻率范圍為1 GHz~5 GHz,其他實驗參數(shù)與極化格雷互補碼探地雷達的參數(shù)相同。圖12(a1)~(a3)所示為基于極化步進頻率雷達的探測結(jié)果。將其與圖11(a1)~(a3)進行對比可得,極化格雷互補碼探地雷達與步進頻探地雷達相比,后者的成像結(jié)果中虛假目標較多,目標響應較弱,因此,管線不能與周圍背景介質(zhì)進行較好的區(qū)分,進而無法準確對管線的走向分布進行辨別。而極化格雷互補碼探地雷達可以實現(xiàn)對管線的走向分布進行明顯的識別。此外,極化格雷互補碼探地雷達的探測結(jié)果中目標位置與周圍環(huán)境的色差更大,這意味著其具備探測更深目標的能力。

    圖11 極化互補格雷碼探地雷達極化響應圖

    圖12 極化步進頻率探地雷達極化響應圖

    進一步地,采用四種無參考指標對極化格雷互補碼與極化步進頻率雷達系統(tǒng)成像結(jié)果進行量化分析。其中圖像的信息熵[18]常被用來評價一個系統(tǒng)中的有效信息量。其定義如公式(8)所示:

    式中:H(X)代表圖像矩陣信息熵值。p(xi)表示圖像矩陣中灰度值為i的像素數(shù)與總像素數(shù)的比值,n為圖像的灰度值。應用到探地雷達成像上,熵值的大小可以用來反映圖像中平均信息含量的多少,即目標體極化信息表現(xiàn)的多少。信息熵越大代表圖像包含信息越多、圖像質(zhì)量越好。

    平均梯度[19]是對圖像局部色差的累計,平均梯度越大,代表圖像細節(jié)反差越強烈,圖像越清晰。其計算公式如式(9)所示:

    公式(9)中M、N代表圖像f(x,y)的行數(shù)與列數(shù)。

    空間頻率[20]作為一種圖像矩陣評價指標,用于測量圖像的整體活動性,其值越高代表圖像所攜帶的信息量越多,細節(jié)表現(xiàn)能力越強。其定義為:

    式(10)、式(11)、式(12)中,Z(n,m)代表圖像在(n,m)處的灰度值。RZ和CZ代表圖像的行頻率和列頻率。

    邊緣強度[21]是凸顯圖像局部變化的幅值,實質(zhì)上是邊緣點之間的梯度幅值,邊緣強度越大,則凸顯邊緣越清晰,失真越小。其計算公式如式(13)所示:

    式中:E(x,y)代表像素點邊緣強度大小,f(x+m,y+n)表示相鄰像素點的灰度值,M×N為所選窗口大小,mean(x-m:x+m,y-n:y+n)為所選區(qū)域的均值。以上量化指標計算結(jié)果如表1 所示:

    表1 極化格雷互補碼探地雷達與極化步進頻探地雷達參考指標對比圖

    由表1 可知,極化格雷互補碼探地雷達的信息熵、平均梯度、空間頻率和邊緣強度四個量化指標均優(yōu)于極化步進頻探地雷達。該結(jié)果表明前者所攜帶的目標信息多于后者,成像更加清晰。為了更加直觀的對兩種體制雷達進行對比分析,圖13(a)、(b)分別給出了極化格雷互補碼探地雷達與極化步進頻探地雷達的三維成像顯示圖,從圖中可以看出,相較于極化步進頻探地雷達,采用極化格雷互補碼探地雷達可以對管線的走向分布進行精準辨別與定位,且能夠有效提高探測深度。

    圖13 三維成像顯示圖

    4 結(jié)論

    提出以格雷互補碼與極化測量相結(jié)合的方法探測地下管線,實現(xiàn)對管線走向分布的精準探測。將格雷互補碼作為探測信號,其具有良好的相關特性并且具有抗干擾性強和分辨率高的特點。而收發(fā)天線采用不同的極化方式可以根據(jù)目標屬性和方向增強目標響應。本文分析了極化格雷互補碼探地雷達的基本理論,并基于FDTD 對管線進行了一維和二維軌線掃描仿真。一維軌線仿真結(jié)果表明,格雷互補碼作為探測信號結(jié)合極化探測可以對地下不同角度管線進行識別,不同極化方式可以表征管線的不同狀態(tài)。二維軌線仿真結(jié)果表明,采用多極化方式可以對地下管線走向進行識別。除此之外,與傳統(tǒng)的步進頻信號探地雷達進行了對比,分別從圖像信息熵、平均梯度、空間頻率和邊緣強度四個指標進行量化分析。實驗結(jié)果表明:在同等條件下,極化格雷互補碼探地雷達的探測效果更好,所獲得的目標信息量更多,成像更清晰,同時目標極化信息損失較少,細節(jié)信息表現(xiàn)明顯。本文為城市地下管線的分布探測提供了一種新方法。

    猜你喜歡
    旁瓣極化格雷
    基于圓柱陣通信系統(tǒng)的廣義旁瓣對消算法
    認知能力、技術進步與就業(yè)極化
    一種基于線性規(guī)劃的頻率編碼旁瓣抑制方法
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    我們生活在格雷河畔
    基于加權積分旁瓣最小化的隨機多相碼設計
    基于四項最低旁瓣Nuttall窗的插值FFT諧波分析
    氯吡格雷治療不穩(wěn)定型心絞痛臨床觀察
    基于PWM控制的新型極化電源設計與實現(xiàn)
    電源技術(2015年1期)2015-08-22 11:16:18
    《道林·格雷的畫像》中的心理解讀
    城市地理(2015年24期)2015-08-15 00:52:57
    久久亚洲精品不卡| 国产精华一区二区三区| 女同久久另类99精品国产91| 国产成人影院久久av| 久久久久久九九精品二区国产| 麻豆av噜噜一区二区三区| 观看美女的网站| 国产精品一区二区三区四区免费观看 | 在线国产一区二区在线| 久久久久久久久大av| 88av欧美| 又紧又爽又黄一区二区| 禁无遮挡网站| 成人国产一区最新在线观看| 韩国av一区二区三区四区| 天天躁日日操中文字幕| 亚洲内射少妇av| 亚洲一区二区三区不卡视频| 夜夜夜夜夜久久久久| 国产91精品成人一区二区三区| 亚洲精品日韩av片在线观看| 91在线观看av| 男人和女人高潮做爰伦理| 在线观看av片永久免费下载| 国产精品99久久久久久久久| 成人欧美大片| 桃色一区二区三区在线观看| 99久久无色码亚洲精品果冻| 在线a可以看的网站| 国产亚洲av嫩草精品影院| 亚洲午夜理论影院| 男插女下体视频免费在线播放| 高潮久久久久久久久久久不卡| 制服丝袜大香蕉在线| 九九久久精品国产亚洲av麻豆| 国产精品1区2区在线观看.| 国产亚洲精品综合一区在线观看| av女优亚洲男人天堂| 国产一区二区在线观看日韩| 国产男靠女视频免费网站| .国产精品久久| 亚洲男人的天堂狠狠| 少妇裸体淫交视频免费看高清| 99久久九九国产精品国产免费| 成年人黄色毛片网站| 亚洲国产精品成人综合色| 中文字幕熟女人妻在线| 色综合婷婷激情| 看黄色毛片网站| 亚洲人与动物交配视频| 波多野结衣高清作品| 日韩欧美 国产精品| 久久精品国产99精品国产亚洲性色| 丰满人妻一区二区三区视频av| 757午夜福利合集在线观看| 日韩欧美在线二视频| 亚洲人成伊人成综合网2020| 一区福利在线观看| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 久久久久久久久久成人| 欧美黑人巨大hd| 国产白丝娇喘喷水9色精品| 乱人视频在线观看| 一个人看视频在线观看www免费| 亚洲无线在线观看| 全区人妻精品视频| 又爽又黄a免费视频| 国产亚洲欧美在线一区二区| 麻豆久久精品国产亚洲av| 国产精品永久免费网站| 99热只有精品国产| 99久国产av精品| 日本三级黄在线观看| 极品教师在线免费播放| 中文资源天堂在线| 亚洲成人久久爱视频| 色在线成人网| 日日夜夜操网爽| 久久久久久国产a免费观看| 亚洲美女黄片视频| 日本成人三级电影网站| 国产高清三级在线| 亚洲最大成人手机在线| 观看美女的网站| а√天堂www在线а√下载| a在线观看视频网站| 午夜老司机福利剧场| 久久精品综合一区二区三区| 一级毛片久久久久久久久女| 网址你懂的国产日韩在线| 老女人水多毛片| 偷拍熟女少妇极品色| 亚洲成人久久爱视频| 国产精品伦人一区二区| 国内久久婷婷六月综合欲色啪| 午夜久久久久精精品| 一二三四社区在线视频社区8| 欧美bdsm另类| 国产大屁股一区二区在线视频| 夜夜看夜夜爽夜夜摸| 特大巨黑吊av在线直播| 91午夜精品亚洲一区二区三区 | 久99久视频精品免费| 精品欧美国产一区二区三| 国内揄拍国产精品人妻在线| 国产精品国产高清国产av| 能在线免费观看的黄片| 午夜福利欧美成人| 看片在线看免费视频| 亚洲aⅴ乱码一区二区在线播放| 美女高潮喷水抽搐中文字幕| av女优亚洲男人天堂| 午夜福利在线观看免费完整高清在 | 麻豆国产97在线/欧美| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 欧美激情国产日韩精品一区| 色视频www国产| 搡女人真爽免费视频火全软件 | 色综合站精品国产| 不卡一级毛片| 亚洲国产日韩欧美精品在线观看| 亚洲av不卡在线观看| 国模一区二区三区四区视频| 草草在线视频免费看| 脱女人内裤的视频| 简卡轻食公司| 黄片小视频在线播放| 国产精品国产高清国产av| 麻豆av噜噜一区二区三区| 嫁个100分男人电影在线观看| 人人妻人人澡欧美一区二区| 丁香六月欧美| 成人无遮挡网站| 亚洲国产精品999在线| 日日摸夜夜添夜夜添av毛片 | 国产黄a三级三级三级人| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 亚洲黑人精品在线| 午夜福利18| 日本熟妇午夜| 欧美极品一区二区三区四区| 久久亚洲精品不卡| 伊人久久精品亚洲午夜| 中文字幕精品亚洲无线码一区| 女人十人毛片免费观看3o分钟| 真人做人爱边吃奶动态| 国产成人aa在线观看| 国产黄a三级三级三级人| 男女下面进入的视频免费午夜| 国产精品不卡视频一区二区 | 国产成+人综合+亚洲专区| 国产成人aa在线观看| 草草在线视频免费看| 久久精品国产亚洲av涩爱 | 性色av乱码一区二区三区2| 久久99热6这里只有精品| 又粗又爽又猛毛片免费看| 听说在线观看完整版免费高清| 欧美xxxx黑人xx丫x性爽| 黄色日韩在线| 久久久久久九九精品二区国产| 免费大片18禁| 国产一级毛片七仙女欲春2| 国产日本99.免费观看| 久久国产乱子免费精品| 伦理电影大哥的女人| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美激情综合另类| 最近最新免费中文字幕在线| 黄片小视频在线播放| 免费无遮挡裸体视频| 国产亚洲欧美98| 啪啪无遮挡十八禁网站| 日韩中字成人| 在线观看66精品国产| x7x7x7水蜜桃| 91久久精品电影网| 日韩欧美一区二区三区在线观看| 伦理电影大哥的女人| 日韩高清综合在线| 少妇丰满av| 国产午夜精品论理片| av视频在线观看入口| 少妇高潮的动态图| 免费人成在线观看视频色| 国产av麻豆久久久久久久| 国产精品乱码一区二三区的特点| 91九色精品人成在线观看| 人人妻,人人澡人人爽秒播| 亚洲男人的天堂狠狠| 成熟少妇高潮喷水视频| 天堂av国产一区二区熟女人妻| 国产日本99.免费观看| 首页视频小说图片口味搜索| 精品99又大又爽又粗少妇毛片 | 亚洲一区二区三区色噜噜| 老司机深夜福利视频在线观看| 午夜视频国产福利| 国产一区二区激情短视频| 久久久久亚洲av毛片大全| 欧美黄色淫秽网站| 久久婷婷人人爽人人干人人爱| 欧美国产日韩亚洲一区| 欧美精品国产亚洲| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美日韩乱码在线| 2021天堂中文幕一二区在线观| a在线观看视频网站| 麻豆成人av在线观看| 噜噜噜噜噜久久久久久91| 99久久99久久久精品蜜桃| 免费在线观看亚洲国产| 少妇人妻精品综合一区二区 | 久久精品久久久久久噜噜老黄 | 欧美黄色片欧美黄色片| 国产精品亚洲美女久久久| 久久久久免费精品人妻一区二区| 此物有八面人人有两片| 亚洲第一电影网av| eeuss影院久久| 夜夜爽天天搞| 99国产精品一区二区三区| 99久久99久久久精品蜜桃| 成熟少妇高潮喷水视频| 久久性视频一级片| 国模一区二区三区四区视频| 深爱激情五月婷婷| 两个人视频免费观看高清| 色吧在线观看| 两个人的视频大全免费| 少妇的逼水好多| 1024手机看黄色片| 国产av在哪里看| 午夜日韩欧美国产| 色av中文字幕| 日本 欧美在线| 国产精品影院久久| 国产毛片a区久久久久| 色av中文字幕| 美女 人体艺术 gogo| 少妇被粗大猛烈的视频| 成人av一区二区三区在线看| 在线十欧美十亚洲十日本专区| 国产中年淑女户外野战色| 免费电影在线观看免费观看| 在线观看66精品国产| 久久草成人影院| 在线观看午夜福利视频| 有码 亚洲区| 好男人电影高清在线观看| 国产在视频线在精品| 久久精品夜夜夜夜夜久久蜜豆| 欧美色欧美亚洲另类二区| 日日干狠狠操夜夜爽| 亚洲av不卡在线观看| 国产精品野战在线观看| 欧美又色又爽又黄视频| 成人午夜高清在线视频| 成人美女网站在线观看视频| 2021天堂中文幕一二区在线观| 日本黄色视频三级网站网址| 欧美+日韩+精品| 好看av亚洲va欧美ⅴa在| 1024手机看黄色片| 亚洲国产色片| 床上黄色一级片| 亚洲精品一区av在线观看| 动漫黄色视频在线观看| 熟女电影av网| 黄色配什么色好看| 国产精品不卡视频一区二区 | 亚洲av成人av| 国内精品一区二区在线观看| 高清毛片免费观看视频网站| 伦理电影大哥的女人| 国产精品久久久久久精品电影| 天堂网av新在线| 久久亚洲真实| 99久久久亚洲精品蜜臀av| 人人妻人人看人人澡| 欧美高清性xxxxhd video| 在线免费观看不下载黄p国产 | 国产乱人视频| 身体一侧抽搐| 国产欧美日韩精品一区二区| 亚洲五月天丁香| 我的女老师完整版在线观看| 91九色精品人成在线观看| 国产高清三级在线| xxxwww97欧美| 桃色一区二区三区在线观看| 亚洲五月婷婷丁香| 久久亚洲真实| 精品久久国产蜜桃| 亚洲熟妇熟女久久| 毛片女人毛片| 最近中文字幕高清免费大全6 | 淫秽高清视频在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲avbb在线观看| 亚洲熟妇中文字幕五十中出| av在线老鸭窝| 中文字幕熟女人妻在线| 中文字幕av成人在线电影| 国产精品伦人一区二区| 久久久久亚洲av毛片大全| 免费观看的影片在线观看| 波野结衣二区三区在线| 色吧在线观看| 最近最新中文字幕大全电影3| 性色av乱码一区二区三区2| 亚洲欧美日韩无卡精品| 国产老妇女一区| 搡老妇女老女人老熟妇| 久久香蕉精品热| 看十八女毛片水多多多| 国产 一区 欧美 日韩| 久99久视频精品免费| 国产视频一区二区在线看| 成人无遮挡网站| 在线播放无遮挡| 亚洲综合色惰| 色尼玛亚洲综合影院| 最新中文字幕久久久久| 色av中文字幕| 国产爱豆传媒在线观看| 精品人妻一区二区三区麻豆 | 国产人妻一区二区三区在| 国产精品久久久久久亚洲av鲁大| 别揉我奶头 嗯啊视频| 亚洲精品粉嫩美女一区| 国产成人av教育| 99久久99久久久精品蜜桃| 亚洲成av人片免费观看| 美女高潮喷水抽搐中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成av人片免费观看| 国内揄拍国产精品人妻在线| 欧美丝袜亚洲另类 | 国产毛片a区久久久久| 精品午夜福利视频在线观看一区| 色播亚洲综合网| 亚洲国产高清在线一区二区三| 乱码一卡2卡4卡精品| 在线观看午夜福利视频| 久久久久久久午夜电影| 国产伦精品一区二区三区视频9| 宅男免费午夜| 变态另类丝袜制服| 男女做爰动态图高潮gif福利片| 黄色配什么色好看| 国产精品嫩草影院av在线观看 | 99在线视频只有这里精品首页| 青草久久国产| 亚洲精品色激情综合| 亚洲国产欧美人成| 亚洲av美国av| 精品不卡国产一区二区三区| 在线播放无遮挡| 亚洲av美国av| 久久久精品大字幕| 欧美日韩亚洲国产一区二区在线观看| 一级a爱片免费观看的视频| 在线播放无遮挡| 国产精品久久久久久久久免 | 亚洲av二区三区四区| 直男gayav资源| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久免费视频| 色播亚洲综合网| 亚洲成人久久性| 男插女下体视频免费在线播放| 在线观看66精品国产| av女优亚洲男人天堂| 人妻夜夜爽99麻豆av| 一个人看视频在线观看www免费| 亚洲人成网站在线播| 欧美乱妇无乱码| 亚洲自偷自拍三级| 日本成人三级电影网站| 色哟哟·www| 欧美日本亚洲视频在线播放| 欧美日韩综合久久久久久 | 在线观看av片永久免费下载| 久久久精品大字幕| 美女高潮的动态| 久久久精品欧美日韩精品| 欧美日韩中文字幕国产精品一区二区三区| 91字幕亚洲| netflix在线观看网站| 黄色丝袜av网址大全| 亚洲成av人片免费观看| 亚洲内射少妇av| 又粗又爽又猛毛片免费看| 精品久久久久久,| 国产午夜精品久久久久久一区二区三区 | 欧美xxxx黑人xx丫x性爽| 国产激情偷乱视频一区二区| 国产精品亚洲一级av第二区| 日本免费一区二区三区高清不卡| 一级a爱片免费观看的视频| 岛国在线免费视频观看| 精品日产1卡2卡| 最新在线观看一区二区三区| 亚洲人成网站高清观看| 日韩精品青青久久久久久| avwww免费| а√天堂www在线а√下载| 免费在线观看亚洲国产| 国产大屁股一区二区在线视频| 一区二区三区激情视频| 日本与韩国留学比较| 国产精品日韩av在线免费观看| 美女 人体艺术 gogo| 男人狂女人下面高潮的视频| 精品人妻偷拍中文字幕| 精品久久久久久久末码| 看黄色毛片网站| 在线观看一区二区三区| 一区二区三区免费毛片| 老熟妇乱子伦视频在线观看| 最近中文字幕高清免费大全6 | 欧美日韩综合久久久久久 | 99热精品在线国产| 国产精品一及| 欧美丝袜亚洲另类 | 日本黄色片子视频| 日本与韩国留学比较| 国产日本99.免费观看| 永久网站在线| 1024手机看黄色片| 999久久久精品免费观看国产| 亚洲狠狠婷婷综合久久图片| a级一级毛片免费在线观看| 长腿黑丝高跟| 真人做人爱边吃奶动态| 久久热精品热| 嫩草影院精品99| 直男gayav资源| 精品久久久久久久久亚洲 | 99热这里只有是精品在线观看 | 亚洲最大成人手机在线| 99视频精品全部免费 在线| 超碰av人人做人人爽久久| 欧美精品啪啪一区二区三区| av视频在线观看入口| 久久久久久久久中文| 日本a在线网址| 国产午夜精品论理片| 色综合站精品国产| 成人欧美大片| 女生性感内裤真人,穿戴方法视频| 99久久精品一区二区三区| 一a级毛片在线观看| 欧美丝袜亚洲另类 | 亚洲中文字幕日韩| 免费大片18禁| 日本熟妇午夜| 亚洲最大成人av| 啦啦啦观看免费观看视频高清| 精华霜和精华液先用哪个| 99热这里只有精品一区| 午夜福利18| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 3wmmmm亚洲av在线观看| 国产 一区 欧美 日韩| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 精品久久久久久,| 欧美高清性xxxxhd video| 老熟妇仑乱视频hdxx| 精品一区二区三区视频在线| 成人国产一区最新在线观看| 成年女人毛片免费观看观看9| 午夜福利18| 99久久成人亚洲精品观看| 国产免费av片在线观看野外av| 啦啦啦韩国在线观看视频| 五月玫瑰六月丁香| 久久久久久久久久成人| 两个人视频免费观看高清| 亚洲国产高清在线一区二区三| 久久亚洲真实| 在线观看免费视频日本深夜| 国产精品久久视频播放| 日韩免费av在线播放| 亚洲五月婷婷丁香| 精品一区二区三区人妻视频| 日本 欧美在线| 精品国产亚洲在线| 男人和女人高潮做爰伦理| 91午夜精品亚洲一区二区三区 | 91在线精品国自产拍蜜月| 欧美精品啪啪一区二区三区| 在现免费观看毛片| 一进一出抽搐gif免费好疼| 欧美黄色淫秽网站| 久久国产乱子免费精品| 欧洲精品卡2卡3卡4卡5卡区| 国产私拍福利视频在线观看| 91在线观看av| 日本a在线网址| 全区人妻精品视频| 成年人黄色毛片网站| 国产三级中文精品| 免费在线观看成人毛片| 成人三级黄色视频| 国产久久久一区二区三区| 久久这里只有精品中国| 日日夜夜操网爽| 他把我摸到了高潮在线观看| 精品久久久久久久人妻蜜臀av| 九九热线精品视视频播放| 黄色一级大片看看| 欧美黑人巨大hd| 国产精品1区2区在线观看.| 男人和女人高潮做爰伦理| 欧美日韩黄片免| 真人一进一出gif抽搐免费| 欧美区成人在线视频| 一级作爱视频免费观看| 中文字幕av成人在线电影| 欧美日韩乱码在线| 免费看光身美女| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片 | 免费观看的影片在线观看| 白带黄色成豆腐渣| 欧美成人a在线观看| 久久精品人妻少妇| 在线免费观看不下载黄p国产 | 一级毛片久久久久久久久女| 国产美女午夜福利| 国产精品亚洲一级av第二区| 少妇裸体淫交视频免费看高清| 天堂av国产一区二区熟女人妻| 丰满人妻一区二区三区视频av| 老司机午夜福利在线观看视频| 成人国产一区最新在线观看| 日本 欧美在线| 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| av专区在线播放| 两个人视频免费观看高清| 波野结衣二区三区在线| 一边摸一边抽搐一进一小说| 午夜福利在线在线| 色综合欧美亚洲国产小说| 亚洲,欧美精品.| 真人做人爱边吃奶动态| 国产黄色小视频在线观看| 午夜视频国产福利| 嫩草影院新地址| netflix在线观看网站| 露出奶头的视频| 免费大片18禁| 桃色一区二区三区在线观看| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 又紧又爽又黄一区二区| 99riav亚洲国产免费| 亚洲av二区三区四区| 性色av乱码一区二区三区2| 特级一级黄色大片| 国产精品三级大全| 真实男女啪啪啪动态图| 97人妻精品一区二区三区麻豆| 久久99热6这里只有精品| 成人国产一区最新在线观看| 精品人妻熟女av久视频| 午夜两性在线视频| 亚洲欧美日韩卡通动漫| 在线十欧美十亚洲十日本专区| 91久久精品国产一区二区成人| 老司机深夜福利视频在线观看| 欧美性猛交╳xxx乱大交人| 精品久久久久久久人妻蜜臀av| 我要看日韩黄色一级片| 人人妻人人看人人澡| 国产一区二区激情短视频| 91九色精品人成在线观看| 乱人视频在线观看| 18+在线观看网站| 亚洲,欧美精品.| 亚洲avbb在线观看| 一级a爱片免费观看的视频| 国产单亲对白刺激| 美女免费视频网站| 可以在线观看的亚洲视频| 小说图片视频综合网站| 亚洲av一区综合| 亚洲国产日韩欧美精品在线观看| 久久久久久国产a免费观看| 99热精品在线国产| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 性欧美人与动物交配| 日本与韩国留学比较| 亚洲欧美激情综合另类| 亚洲精品在线观看二区| 禁无遮挡网站| 国产乱人视频| bbb黄色大片| 中文字幕av成人在线电影| 精品一区二区三区视频在线| 亚洲中文字幕一区二区三区有码在线看| 成年人黄色毛片网站| 亚洲人成网站在线播放欧美日韩|