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

    一種水下聲成像算法研究及擴展

    2017-12-23 06:09:06作者廖鵬中國船舶重工集團公司第七一研究所
    電子制作 2017年18期
    關(guān)鍵詞:簡正波聲速聲場

    作者/廖鵬,中國船舶重工集團公司第七一〇研究所

    一種水下聲成像算法研究及擴展

    作者/廖鵬,中國船舶重工集團公司第七一〇研究所

    本文闡述了水下聲成像的基本原理與研究的重要意義,利用Kraken簡正波模型進行了水下聲場傳播模型的仿真,分析了水中聲波在距離上和深度上的衰減情況。推導(dǎo)了一種三陣元水聽器系統(tǒng)的成像方法并給出了仿真結(jié)果,最后將算法由二維成像擴展到三維成像。

    水下聲成像;Kraken簡正波模型;三維成像

    引言

    水下聲成像是聲成像的一個分支,在水聲探測中有著重要的作用,通過物體對聲波的后向散射作用進行水下定位,逐點畫圖成像。利用聲波進行成像,是在聲納測距與測向的基礎(chǔ)上發(fā)展而來的。聲納測距主要是利用接收回波與發(fā)射脈沖信號間的時間差,計算出目標(biāo)的距離;聲納測向主要是利用回波到達由多個換能器組成的水聽器系統(tǒng)的聲程差和相位差,計算出目標(biāo)的方位。

    在水下聲成像中,基本的聲成像技術(shù)有三種:聲透鏡技術(shù)、波束形成技術(shù)、聲全息技術(shù)。這三種聲成像方法都使用相同的操作:空間處理(從聲場中得到圖像)、換能(將聲能轉(zhuǎn)換成電能)、檢波(將高頻信號轉(zhuǎn)換成可觀測且接近直流的圖像信號)、顯示(以某種形式顯示為圖像)。這三種方法的不同之處在于對這些操作進行的先后順序的不同。

    本文使用的是一種與以上三種均不同的新的聲成像方法。三個換能器的空間位置形成各自距離相等正三角形,組成一個水聽器系統(tǒng)。利用回波信號到其中兩個等效陣元之間的相位差、時間差和振幅差值,測定各個目標(biāo)點的距離與方位。

    1.水下聲場仿真

    目前用于計算淺海混響的理論模型主要有射線模型、簡正波模型,本文使用簡正波模型進行水下聲場仿真。利用簡正波模型來計算水下混響的基本思想最早由Bucker和Morris于1968年提出,簡正波混響模型用數(shù)學(xué)公式描述為:

    其中,R(t)是t時刻的混響強度,I0是持續(xù)時間為τ0的聲源強度,Gm是指經(jīng)路徑m從聲源到散射元dA的傳播。Gn是指經(jīng)路徑n從散射元dA到接收器的傳播。Smn是指單位區(qū)域?qū)θ肷渎窂絤外向路徑n的散射,rmn是指從聲源或接收器到散射元dA的距離,dA基本的散射區(qū)域。求和是相對于連接聲源與接收器和散射元的所有路徑(射線或模態(tài))。積分是相對于t時刻對混響有貢獻的所有散射區(qū)域。

    水下聲場仿真主要由三部分組成,分別為目標(biāo)亮點回波信號產(chǎn)生、前置預(yù)處理和后置處理。

    (1)目標(biāo)亮點回波產(chǎn)生:主動發(fā)射一個載頻信號,經(jīng)信道傳輸,由反射點時延、相位時延形成反射后的回波信號,總目標(biāo)反射的回波信號實際上是由目標(biāo)亮點回波的相干疊加合成。

    (2)前置預(yù)處理:完成對信號的預(yù)處理,為后續(xù)波束形成和圖像生成提供可靠的數(shù)據(jù),功能主要包括對仿真信號的采用、高頻帶通濾波器、正交解調(diào)、低通濾波和低頻抽取等。

    (3)后置處理:完成前置預(yù)處理后仿真信號的復(fù)數(shù)字波束形成,對目標(biāo)方位、距離、強度等信息進行估計,再利用圖像色彩顯示控制技術(shù)最終生成瞬態(tài)回波圖像,即目標(biāo)二維亮點強度分布區(qū)。

    圖1 淺海聲速剖面圖

    圖2 深海聲速剖面圖

    海洋中的平均聲速近似等于1500m/s,聲速隨溫度、鹽度、壓力的增加而增加,其中以溫度的影響最顯著。聲速在海中的聲速剖面圖如圖1和圖2所示。圖1表示的是淺海(0~250米)的聲速剖面圖,圖2表示的是深海(0~4000米)的聲速剖面圖。在淺海,海洋表面受到陽光照射時水溫較高,聲速隨著深度增加而降低,呈現(xiàn)聲速負(fù)梯度;在深海,水溫比較低而且穩(wěn)定,聲速隨著深度的增加而增加,呈現(xiàn)海洋內(nèi)部的聲速正梯度。

    根據(jù)淺海的聲速梯度分布,利用Kraken簡正波模型對淺海水下聲傳播進行仿真。聲速剖面圖如圖1所示,假設(shè)海面為3級海況,水深為250m,聲源頻率為500Hz,聲源深度為60m,海水密度為1.03g/cm2,吸收系數(shù)等于0.15dB/kmHz。聲場仿真圖如下圖3所示。

    圖3 淺海水下聲場仿真圖

    圖3的左邊的圖表示聲波在深度和距離上的傳播損失,右圖表示的是在深度60m處,即聲源位置處開始,聲波隨距離的衰減趨勢。由圖3可以看出,越靠近聲源的位置聲強越大,但是水下聲傳播過程中的衰減并不是線性的。

    聲場中任意點的聲壓公式為:

    且p(r,z, ω) =S(ω)g(r,z,ξ),上式可化為:

    其中,S(ω)為聲源信號,g(r,z,ω)為聲場的傳播函數(shù)。

    2.聲成像算法研究

    本文假設(shè)的聲納系統(tǒng)共有三個收發(fā)合制的換能器,三個換能器的空間位置成正三角形,同時接收各個方向上的回波信號。DDS信號發(fā)生器產(chǎn)生22kHz~24kHz的線性調(diào)頻信號LFM,經(jīng)過三個發(fā)射換能器輪流向水中或者空氣中發(fā)射出去。三個換能器的空間位置是相互固定的,它們理論上可以根據(jù)接收到的回波信號確定空間任一點的位置。下面以平面內(nèi)兩個換能器確定平面內(nèi)任意一點為例,建立理想情況下的聲傳播模型(即認(rèn)為聲在介質(zhì)中的傳播路徑是一條直線)。圖4為平面內(nèi)信號的收發(fā)模型。

    圖4 平面內(nèi)兩換能器收發(fā)信號模型

    圖5 線性調(diào)頻信號

    換能器A、B發(fā)射的信號為線性調(diào)頻信號,表達式如下:

    圖5為此線性調(diào)頻信號在MATLAB下的波形圖。公式中信噪比SNR=10,載波頻率fc=5kHz,帶寬B=5kHz,線性調(diào)頻信號脈沖寬度TP=0.01s(只發(fā)射一段時間很短的脈沖),采樣頻率fs=30kHz。在已有發(fā)射信號LFM的基礎(chǔ)上產(chǎn)生的回波信號,可以表達為在原信號的基礎(chǔ)上疊加其延時衰減的分量。假設(shè)只有一個回聲的情況下,可簡化其模型為:

    其中α為反射系數(shù),N為延遲時間。為了簡化計算,假設(shè)平面一共有10個點,換能器A、B的間距為1個單位,水平面的長度和高度均為10個單位,換能器發(fā)射的信號傳播至每個點上都會產(chǎn)生不同時延的回波,每個換能器接收到的總的回波信號為各個點的回波信號的總和,表達式為:

    所以,LMF信號在每個點處的回波模型可以表達為:

    rn為點到換能器的直線距離,?t為每個回波信號的時延。換能器接收到的總的回波信號表達式為:

    對回波信號進行自相關(guān)運算后得到的函數(shù)存在一些極值點,根據(jù)相關(guān)函數(shù)的定義我們知道,這些極值點就原信號自相關(guān)、原信號與回波信號互相關(guān)、回波信號與原信號互相關(guān)的相關(guān)函數(shù)值,根據(jù)這些點的相對位置、采樣頻率和聲速,就可以確定反射點的位置了。

    在MATLAB仿真中,對發(fā)射信號加窗處理后,與接收到的信號進行卷積,得到的圖形的波峰位置處就是信號的相關(guān)函數(shù)值。

    圖6表示的是換能器A和B各自的發(fā)射信號與接收信號進行卷積運算的結(jié)果,同一個目標(biāo)點在兩個換能器中形成的波峰位置不同,根據(jù)兩個換能器中波峰點的相對位置與采樣頻率的大小,通過交叉運算就可以計算出目標(biāo)點的位置坐標(biāo)。圖7表示的是對10個目標(biāo)點方位估計的仿真結(jié)果。

    由仿真可以看到,兩個換能器就可以確定平面內(nèi)任意點的位置。拓展到三維空間中,則三個換能器理論上就可以確定空間任一點的位置。在具體的實驗中,首先利用三個換能器中的兩個畫出一系列二維圖像,在將這些圖像拼接成一個完整的三維圖像。

    圖6 發(fā)射信號與回波信號的卷積運算結(jié)果

    圖7 目標(biāo)點的方位估計

    3.聲成像算法擴展

    上文建立了一個理想的信號模型,利用兩個換能器畫出了平面內(nèi)10個目標(biāo)點的位置,但是在復(fù)雜的環(huán)境中,由于目標(biāo)點數(shù)非常的多,回波信號所含的信息量非常巨大,這種算法效果并不好。

    圖8 信號的發(fā)射與接收

    圖9 逐行掃描信號模型

    圖8表示的是信號的發(fā)射與接收波形示意圖,發(fā)射換能器發(fā)射的信號經(jīng)過一段時間后回來,換能器A與換能器B接收到的回波信號前沿的時延分別為T1、T2。由于兩個換能器的位置不同,接收到的信號也就不同,所以理論上可以利用兩個換能器接收到的信號差異確定目標(biāo)的方位。

    在這里用逐點逐行掃描的方法確定空間點的位置。如圖9所示,以平面內(nèi)兩個換能器A、B為例,目標(biāo)物體為L,利用這兩個換能器畫出L的平面圖。建立平面坐標(biāo)系,將換能器前方的區(qū)域用網(wǎng)格劃分為一系列的點,進行逐點掃描。

    首先建立網(wǎng)格中任一點的聲信號模型。換能器發(fā)射信號為x(t),是一段短脈沖,采用假設(shè)法,假設(shè)在S點出有目標(biāo)點,回波信號為y(t)。利用基于滑動窗口的快速傅立葉變換進行歸一化互相關(guān)匹配,采用滑動窗和快速傅立葉變換(FFT)與反變換(IFFT)對歸一化互相關(guān)系數(shù)進行計算,公式如下:

    其中u= 0,1,...M-1,M為回波信號總長度。用x(t-u)代表標(biāo)準(zhǔn)波形在t-u位置上的幅度值,x是標(biāo)準(zhǔn)波形序列x(t-u)的平均幅度值,yu(t)代表查找區(qū)域上以第u個位置為起點的與標(biāo)準(zhǔn)波形同樣長度區(qū)域內(nèi)的幅度值序列,yu則是此序列的平均幅度值。r(u)構(gòu)成了一個相關(guān)序列,序列上的點均位于[-1,+1]區(qū)間內(nèi),+1意味著兩個波形之間具有完全匹配的正向線性關(guān)系,-1意味著兩個波形之間具有完全匹配的反向線性關(guān)系,在匹配位置將出現(xiàn)峰值,峰值接近于1。

    根據(jù)兩個換能器接收到的回波信號與發(fā)射信號的相關(guān)性,大致可以確定假設(shè)點S是否存在。如果假設(shè)點存在,則匹配濾波的相關(guān)性強,峰值明顯,否則假設(shè)點就不存在,繼續(xù)查找下一點。圖中只有目標(biāo)L才能反射回波信號,所以只有在L處匹配濾波相關(guān)性最強。平面目標(biāo)方位估計仿真如圖10所示。

    圖10 平面目標(biāo)點的方位估計

    圖11 空間方位估計模型

    擴展到空間內(nèi),空間方位估計模型如圖11所示。和平面目標(biāo)方位估計一樣,以三個換能器中的兩個確定某一平面內(nèi)的目標(biāo)點的坐標(biāo),將產(chǎn)生的一系列平面內(nèi)的點結(jié)合起來,形成一個完整的三維圖像。

    4.結(jié)論

    本文闡述了水下聲成像的基本原理與研究的重要意義,推導(dǎo)了一種三基陣元水聽器系統(tǒng)的聲成像算法。利用Kraken簡正波模型,進行了水下聲場仿真,并在此基礎(chǔ)上建立了水下聲信號的信號模型。利用聲信號的理想模型,通過MATLAB軟件進行了仿真,成功實現(xiàn)了平面內(nèi)的目標(biāo)成像,并擴展到三維空間中成像。

    * [1]田坦,劉國枝,孫大軍.聲納技術(shù)[M].哈爾濱工程大學(xué)出版社,2000

    * [2]劉伯勝, 雷家煜主編.水聲學(xué)原理[M].哈爾濱工程大學(xué)出版社,1993.12

    * [3]唐利娜,孫大軍.水下三維成像及聲納圖像處理技術(shù)研究[D].哈爾濱工程大學(xué)碩士學(xué)位論文,2009

    * [4]Lee C.-C., Lee J.-H., Robust adaptive array beamforming under steering vector sensors. IEEE Trans.Antennas propaga te,1997

    * [5]Alex B. Gershman, Nicholas D. Sidiropoulos. Convex optimization-based beamforming.IEEE Signal Processing Magazine,2010

    猜你喜歡
    簡正波聲速聲場
    傾斜彈性海底條件下淺海聲場的簡正波相干耦合特性分析*
    基于BIM的鐵路車站聲場仿真分析研究
    探尋360°全聲場發(fā)聲門道
    一種高效的寬帶簡正波本征值計算方法
    一種快速求解寬頻簡正波的方法
    warping變換提取單模態(tài)反演海底衰減系數(shù)?
    聲速是如何測定的
    跨聲速風(fēng)洞全模顫振試驗技術(shù)
    機翼跨聲速抖振研究進展
    板結(jié)構(gòu)-聲場耦合分析的FE-LSPIM/FE法
    女警被强在线播放| 99久久久亚洲精品蜜臀av| 精品国产乱子伦一区二区三区| 老司机午夜十八禁免费视频| 美女cb高潮喷水在线观看 | 两人在一起打扑克的视频| 欧美色视频一区免费| 久久久久国内视频| 51午夜福利影视在线观看| 性色av乱码一区二区三区2| 一夜夜www| 91在线精品国自产拍蜜月 | 少妇裸体淫交视频免费看高清| 琪琪午夜伦伦电影理论片6080| 国产黄a三级三级三级人| 国产黄a三级三级三级人| 久久久久久国产a免费观看| 久久精品综合一区二区三区| 女生性感内裤真人,穿戴方法视频| 久久精品综合一区二区三区| 午夜福利免费观看在线| 亚洲中文av在线| 成人特级av手机在线观看| 国产一区二区激情短视频| 午夜福利在线观看免费完整高清在 | 午夜福利在线在线| 老熟妇仑乱视频hdxx| 国产一区二区三区视频了| 三级国产精品欧美在线观看 | 成人无遮挡网站| 又大又爽又粗| 18禁裸乳无遮挡免费网站照片| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久com| 国产亚洲精品久久久com| 一级作爱视频免费观看| 精品一区二区三区视频在线 | x7x7x7水蜜桃| 99久久精品国产亚洲精品| 午夜激情福利司机影院| 中文字幕精品亚洲无线码一区| 成人av在线播放网站| 国产成人一区二区三区免费视频网站| 久久午夜亚洲精品久久| 好男人在线观看高清免费视频| 成人性生交大片免费视频hd| 午夜激情欧美在线| 亚洲欧美日韩高清专用| 国产成人影院久久av| 国产亚洲av嫩草精品影院| 给我免费播放毛片高清在线观看| 国产欧美日韩一区二区精品| 可以在线观看的亚洲视频| 国产美女午夜福利| 国产97色在线日韩免费| 欧美成人免费av一区二区三区| 国产激情久久老熟女| 黄色丝袜av网址大全| 精品久久久久久久久久免费视频| 免费在线观看日本一区| 国产1区2区3区精品| 好看av亚洲va欧美ⅴa在| 亚洲最大成人中文| 亚洲精品一卡2卡三卡4卡5卡| 精品福利观看| 国语自产精品视频在线第100页| 黑人巨大精品欧美一区二区mp4| 亚洲中文字幕日韩| 亚洲 国产 在线| 午夜福利18| 国产不卡一卡二| 丰满的人妻完整版| 村上凉子中文字幕在线| 色噜噜av男人的天堂激情| 他把我摸到了高潮在线观看| 国产成年人精品一区二区| 国产黄a三级三级三级人| 99国产综合亚洲精品| 最近最新中文字幕大全免费视频| 久久精品91无色码中文字幕| 久久久久久久精品吃奶| 久久久久国内视频| 操出白浆在线播放| 久久久国产成人精品二区| 一卡2卡三卡四卡精品乱码亚洲| 级片在线观看| 一级毛片女人18水好多| 亚洲成av人片在线播放无| 亚洲欧美日韩东京热| 在线观看免费视频日本深夜| 免费在线观看影片大全网站| 久久久久亚洲av毛片大全| 五月伊人婷婷丁香| 在线永久观看黄色视频| 十八禁网站免费在线| 亚洲性夜色夜夜综合| 99久久无色码亚洲精品果冻| 禁无遮挡网站| 久久这里只有精品19| 成人亚洲精品av一区二区| 中文在线观看免费www的网站| 欧美最黄视频在线播放免费| 亚洲国产日韩欧美精品在线观看 | 国产精品1区2区在线观看.| 成年女人看的毛片在线观看| 日韩欧美国产一区二区入口| 美女午夜性视频免费| 88av欧美| 在线观看美女被高潮喷水网站 | 色综合婷婷激情| 熟妇人妻久久中文字幕3abv| 国产乱人伦免费视频| 麻豆久久精品国产亚洲av| 亚洲真实伦在线观看| 国产精品av视频在线免费观看| 久久久水蜜桃国产精品网| 我要搜黄色片| 男插女下体视频免费在线播放| 国模一区二区三区四区视频 | 一区二区三区高清视频在线| 丁香欧美五月| 美女高潮喷水抽搐中文字幕| 国产真人三级小视频在线观看| 精品国产三级普通话版| 欧美黄色片欧美黄色片| 亚洲性夜色夜夜综合| 欧美成人一区二区免费高清观看 | 久久国产精品影院| 伊人久久大香线蕉亚洲五| 国产一区在线观看成人免费| 99精品欧美一区二区三区四区| 淫秽高清视频在线观看| 国产成人精品无人区| 一本一本综合久久| 免费看日本二区| 高潮久久久久久久久久久不卡| 亚洲一区二区三区不卡视频| 不卡av一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 男女床上黄色一级片免费看| 亚洲成人中文字幕在线播放| 久久婷婷人人爽人人干人人爱| 欧美黑人欧美精品刺激| 深夜精品福利| 香蕉av资源在线| 他把我摸到了高潮在线观看| 草草在线视频免费看| 女生性感内裤真人,穿戴方法视频| 亚洲欧美日韩卡通动漫| 国产欧美日韩一区二区精品| 亚洲欧美精品综合久久99| 天堂网av新在线| 在线十欧美十亚洲十日本专区| 麻豆久久精品国产亚洲av| 国产一区在线观看成人免费| 可以在线观看毛片的网站| 免费无遮挡裸体视频| 18禁观看日本| 精品国产三级普通话版| a级毛片a级免费在线| 欧美成人一区二区免费高清观看 | 精品人妻1区二区| 久久久成人免费电影| 亚洲性夜色夜夜综合| 亚洲五月婷婷丁香| 美女高潮的动态| 午夜精品一区二区三区免费看| 色老头精品视频在线观看| 久久久久久九九精品二区国产| 男女下面进入的视频免费午夜| 国产一区二区激情短视频| 亚洲五月天丁香| av片东京热男人的天堂| 国产在线精品亚洲第一网站| 一个人看视频在线观看www免费 | 男人舔女人的私密视频| 一进一出好大好爽视频| 叶爱在线成人免费视频播放| 午夜精品在线福利| 亚洲第一欧美日韩一区二区三区| 精品久久久久久久久久免费视频| 又黄又爽又免费观看的视频| 大型黄色视频在线免费观看| 国产欧美日韩一区二区三| 观看美女的网站| 久久久色成人| 日韩欧美在线二视频| 中文资源天堂在线| 国内精品美女久久久久久| 国产又色又爽无遮挡免费看| 亚洲中文av在线| 在线观看午夜福利视频| 女人被狂操c到高潮| 丁香六月欧美| 一二三四在线观看免费中文在| 亚洲精品在线观看二区| 亚洲18禁久久av| 一级毛片高清免费大全| 午夜影院日韩av| 丰满的人妻完整版| 欧美成人性av电影在线观看| 午夜激情欧美在线| 男女那种视频在线观看| 国产精品国产高清国产av| 波多野结衣高清无吗| 丰满人妻一区二区三区视频av | 国产视频一区二区在线看| 巨乳人妻的诱惑在线观看| 欧美一区二区国产精品久久精品| 日本一二三区视频观看| 欧美中文综合在线视频| 一个人看视频在线观看www免费 | 又爽又黄无遮挡网站| 脱女人内裤的视频| 成人国产综合亚洲| 国产男靠女视频免费网站| 亚洲精品国产精品久久久不卡| 嫁个100分男人电影在线观看| 色视频www国产| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽| 午夜福利欧美成人| 综合色av麻豆| 精品熟女少妇八av免费久了| 中文字幕最新亚洲高清| 亚洲午夜理论影院| 国产爱豆传媒在线观看| 天堂动漫精品| 啪啪无遮挡十八禁网站| 又爽又黄无遮挡网站| 9191精品国产免费久久| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看| 51午夜福利影视在线观看| 亚洲成a人片在线一区二区| 99国产精品99久久久久| av在线蜜桃| 亚洲av美国av| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 99国产精品一区二区三区| 人人妻人人澡欧美一区二区| 视频区欧美日本亚洲| 可以在线观看的亚洲视频| 精品一区二区三区av网在线观看| 黄色日韩在线| 国产精品av视频在线免费观看| 最近最新中文字幕大全免费视频| 国产成人欧美在线观看| 91九色精品人成在线观看| 欧美+亚洲+日韩+国产| 女生性感内裤真人,穿戴方法视频| 最近在线观看免费完整版| 久久久久久久久免费视频了| 在线观看66精品国产| 国产真人三级小视频在线观看| 少妇裸体淫交视频免费看高清| 亚洲成a人片在线一区二区| 亚洲av成人一区二区三| 日韩欧美国产一区二区入口| 九九热线精品视视频播放| 无人区码免费观看不卡| 久久久精品大字幕| 亚洲中文字幕一区二区三区有码在线看 | 搡老妇女老女人老熟妇| 两人在一起打扑克的视频| 欧美日本视频| 亚洲精品一区av在线观看| 亚洲性夜色夜夜综合| 国产视频内射| 韩国av一区二区三区四区| 国产男靠女视频免费网站| 成在线人永久免费视频| 男插女下体视频免费在线播放| 国产av一区在线观看免费| 国产精品99久久99久久久不卡| 丁香欧美五月| 99国产精品一区二区蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 国产精品乱码一区二三区的特点| 日韩欧美精品v在线| 香蕉丝袜av| 日本黄色片子视频| 国产伦在线观看视频一区| 亚洲欧美日韩无卡精品| 18禁美女被吸乳视频| 99久久综合精品五月天人人| 天天一区二区日本电影三级| 99re在线观看精品视频| 亚洲aⅴ乱码一区二区在线播放| 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 色噜噜av男人的天堂激情| 嫩草影院入口| 中文在线观看免费www的网站| 久久精品91蜜桃| 一二三四社区在线视频社区8| 国产淫片久久久久久久久 | 亚洲真实伦在线观看| 国内毛片毛片毛片毛片毛片| 又爽又黄无遮挡网站| www日本在线高清视频| 精品国产三级普通话版| 一边摸一边抽搐一进一小说| 亚洲人成网站高清观看| 日韩欧美 国产精品| a级毛片a级免费在线| 亚洲一区高清亚洲精品| 成人国产综合亚洲| 国产精品久久久久久精品电影| 国产亚洲av嫩草精品影院| 精品国内亚洲2022精品成人| 成人国产一区最新在线观看| 国产精品自产拍在线观看55亚洲| av片东京热男人的天堂| 日本熟妇午夜| 午夜福利欧美成人| 亚洲五月天丁香| 九色国产91popny在线| 国产极品精品免费视频能看的| 亚洲国产欧美人成| 非洲黑人性xxxx精品又粗又长| 国产伦一二天堂av在线观看| 巨乳人妻的诱惑在线观看| 精品久久久久久久末码| 男女那种视频在线观看| 大型黄色视频在线免费观看| 久久中文字幕一级| 美女扒开内裤让男人捅视频| 曰老女人黄片| 日韩欧美三级三区| 欧美日韩一级在线毛片| 亚洲无线观看免费| 午夜a级毛片| 毛片女人毛片| 麻豆一二三区av精品| 亚洲真实伦在线观看| 国产成人精品久久二区二区91| 亚洲第一电影网av| 免费在线观看影片大全网站| 日韩精品中文字幕看吧| 久久这里只有精品中国| 嫁个100分男人电影在线观看| 美女被艹到高潮喷水动态| 亚洲乱码一区二区免费版| 亚洲国产精品成人综合色| 最近最新中文字幕大全电影3| 1000部很黄的大片| 成在线人永久免费视频| 精品久久久久久成人av| 精品国产超薄肉色丝袜足j| 亚洲av日韩精品久久久久久密| 国产97色在线日韩免费| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| 亚洲人成伊人成综合网2020| 国产精品一区二区三区四区免费观看 | 午夜精品久久久久久毛片777| 好看av亚洲va欧美ⅴa在| 国产高清激情床上av| 久久久国产成人免费| 特级一级黄色大片| 最新中文字幕久久久久 | 亚洲人成网站高清观看| 国产精品久久久久久亚洲av鲁大| 最近最新中文字幕大全免费视频| 欧美大码av| 国产精品综合久久久久久久免费| 俄罗斯特黄特色一大片| 亚洲第一电影网av| 精品电影一区二区在线| 51午夜福利影视在线观看| 两个人看的免费小视频| 日本五十路高清| 亚洲精华国产精华精| 一级毛片精品| 宅男免费午夜| 日本免费a在线| 欧美国产日韩亚洲一区| 国产精品综合久久久久久久免费| 免费观看人在逋| 最近最新中文字幕大全电影3| 亚洲av成人不卡在线观看播放网| 免费人成视频x8x8入口观看| 欧美最黄视频在线播放免费| 欧美日本亚洲视频在线播放| 18美女黄网站色大片免费观看| 中文字幕最新亚洲高清| 丰满的人妻完整版| 日韩欧美精品v在线| 99热6这里只有精品| 99精品欧美一区二区三区四区| 国产三级黄色录像| 一本精品99久久精品77| 成年女人永久免费观看视频| 丰满人妻熟妇乱又伦精品不卡| 99riav亚洲国产免费| 国产精品 国内视频| 女生性感内裤真人,穿戴方法视频| 国产一级毛片七仙女欲春2| 久久精品人妻少妇| 我要搜黄色片| 久久久国产成人精品二区| 白带黄色成豆腐渣| 国产精品香港三级国产av潘金莲| 18禁美女被吸乳视频| 久久热在线av| 97碰自拍视频| av在线天堂中文字幕| 国产精品 欧美亚洲| 12—13女人毛片做爰片一| a级毛片a级免费在线| 女人被狂操c到高潮| 99视频精品全部免费 在线 | 欧美日韩一级在线毛片| 欧美乱码精品一区二区三区| 丰满的人妻完整版| 一进一出抽搐动态| 国产精品1区2区在线观看.| 成人特级黄色片久久久久久久| 老汉色av国产亚洲站长工具| 欧美成狂野欧美在线观看| 欧美大码av| 欧美绝顶高潮抽搐喷水| 国产一级毛片七仙女欲春2| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av香蕉五月| 亚洲熟妇中文字幕五十中出| 国产黄色小视频在线观看| 黄色视频,在线免费观看| 日本一二三区视频观看| 在线看三级毛片| 午夜精品一区二区三区免费看| 欧美乱色亚洲激情| 亚洲欧美日韩东京热| 欧美三级亚洲精品| 国产一区二区三区视频了| 精品久久久久久,| a级毛片在线看网站| 国产视频内射| 日韩大尺度精品在线看网址| 成人高潮视频无遮挡免费网站| 男人舔女人的私密视频| 黄色女人牲交| 日韩欧美免费精品| 亚洲人成电影免费在线| 国产成人啪精品午夜网站| 日本一本二区三区精品| 免费在线观看视频国产中文字幕亚洲| 亚洲一区二区三区色噜噜| 日韩欧美精品v在线| 91麻豆精品激情在线观看国产| 在线永久观看黄色视频| 色综合亚洲欧美另类图片| 国产又色又爽无遮挡免费看| 中文字幕人成人乱码亚洲影| 亚洲,欧美精品.| 岛国在线免费视频观看| 香蕉国产在线看| 色综合欧美亚洲国产小说| 黑人欧美特级aaaaaa片| 欧美日韩国产亚洲二区| 国产乱人伦免费视频| 色综合亚洲欧美另类图片| 欧美xxxx黑人xx丫x性爽| 两个人看的免费小视频| 亚洲自拍偷在线| 麻豆国产av国片精品| 欧美乱码精品一区二区三区| 色综合站精品国产| 伦理电影免费视频| 男女午夜视频在线观看| 99久久国产精品久久久| 国产亚洲av高清不卡| 51午夜福利影视在线观看| 波多野结衣巨乳人妻| 亚洲精品色激情综合| 久久伊人香网站| 男女下面进入的视频免费午夜| 久久久国产成人精品二区| 久久精品人妻少妇| 精品熟女少妇八av免费久了| 亚洲av五月六月丁香网| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | 又大又爽又粗| 一本综合久久免费| 日韩欧美免费精品| 国产午夜精品论理片| 久久精品综合一区二区三区| 国产成人av激情在线播放| 日韩欧美一区二区三区在线观看| 国产成人aa在线观看| 一区二区三区激情视频| 精品免费久久久久久久清纯| 亚洲avbb在线观看| 日本黄色视频三级网站网址| 99久久无色码亚洲精品果冻| 国产免费男女视频| 久久精品aⅴ一区二区三区四区| 中国美女看黄片| 欧美日韩国产亚洲二区| 久久国产精品人妻蜜桃| 麻豆久久精品国产亚洲av| 精品国产乱码久久久久久男人| 久久伊人香网站| 美女大奶头视频| 一进一出抽搐动态| 国产精品乱码一区二三区的特点| 欧美午夜高清在线| 久久精品国产99精品国产亚洲性色| 一二三四在线观看免费中文在| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 两个人的视频大全免费| 综合色av麻豆| 法律面前人人平等表现在哪些方面| 99视频精品全部免费 在线 | 99久久综合精品五月天人人| 99久久久亚洲精品蜜臀av| 成在线人永久免费视频| 国产激情久久老熟女| 国产一区二区三区视频了| 免费看美女性在线毛片视频| 久久久久精品国产欧美久久久| 成人av一区二区三区在线看| 国产午夜精品久久久久久| 欧美丝袜亚洲另类 | 淫妇啪啪啪对白视频| 国产精品综合久久久久久久免费| 色综合站精品国产| 午夜成年电影在线免费观看| 国产高潮美女av| 亚洲性夜色夜夜综合| 亚洲av成人精品一区久久| 一a级毛片在线观看| 禁无遮挡网站| 亚洲一区二区三区不卡视频| 亚洲18禁久久av| 女人被狂操c到高潮| 在线视频色国产色| 热99re8久久精品国产| 久久久久久人人人人人| 欧美黑人巨大hd| 亚洲av中文字字幕乱码综合| av国产免费在线观看| a在线观看视频网站| 亚洲色图av天堂| a级毛片a级免费在线| 黑人欧美特级aaaaaa片| 色噜噜av男人的天堂激情| 久久久国产成人免费| 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 日韩大尺度精品在线看网址| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 又黄又粗又硬又大视频| 黄片小视频在线播放| 日本黄色片子视频| 成人特级黄色片久久久久久久| 国产精品精品国产色婷婷| 午夜久久久久精精品| 国产男靠女视频免费网站| a在线观看视频网站| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| 亚洲国产精品久久男人天堂| 亚洲真实伦在线观看| a级毛片在线看网站| 一本综合久久免费| 国产黄a三级三级三级人| 伦理电影免费视频| 色在线成人网| 成人性生交大片免费视频hd| 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区mp4| 香蕉国产在线看| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美在线二视频| 久久亚洲精品不卡| 成人18禁在线播放| 亚洲精品美女久久av网站| 看片在线看免费视频| 极品教师在线免费播放| 黄色成人免费大全| 国产亚洲av嫩草精品影院| 此物有八面人人有两片| www国产在线视频色| 精华霜和精华液先用哪个| 老司机福利观看| av欧美777| 91麻豆av在线| 草草在线视频免费看| 午夜久久久久精精品| 国产一区二区激情短视频| 看黄色毛片网站| 免费观看人在逋| 婷婷精品国产亚洲av| 1024手机看黄色片| 亚洲专区中文字幕在线| 午夜激情福利司机影院| 亚洲精品一区av在线观看| 成年免费大片在线观看| 亚洲乱码一区二区免费版| 精品久久蜜臀av无| 18禁黄网站禁片免费观看直播| 亚洲avbb在线观看| 国产免费av片在线观看野外av| 在线播放国产精品三级| 99精品欧美一区二区三区四区| 欧美最黄视频在线播放免费|