999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

相控陣ADCP編碼信號相移波束形成

2015-10-14 01:49:39李鵬飛彭東立
聲學技術(shù) 2015年6期
關(guān)鍵詞:信號系統(tǒng)

李鵬飛,彭東立

?

相控陣ADCP編碼信號相移波束形成

李鵬飛1,彭東立2

(1. 海軍裝備研究院艦艇部,北京100841;2. 中國科學院聲學研究所東海研究站,上海200032)

相對于活塞式聲學多普勒流速剖面儀,相控陣聲學多普勒流速剖面儀的換能器體積大大縮小,并且依托相控陣本身的物理特性,無需進行聲速補償。利用窄帶相控陣聲學多普勒流速剖面儀系統(tǒng)進行編碼信號的發(fā)射與接收,從而提高相控陣聲學多普勒流速剖面儀的設(shè)備性能。首先介紹了相控陣ADCP的波束形成方法,然后分別介紹了相控陣的相移波束形成和時延波束形成原理,對編碼信號相移相控波束形成情況下,相控陣發(fā)射和接收信號波束開角進行了對比,并對作用距離、測量精度、系統(tǒng)復(fù)雜度等系統(tǒng)性能進行了分析。分析結(jié)果表明:窄帶相控陣聲學多普勒流速剖面儀可以進行編碼信號的發(fā)射和接收,能夠提高系統(tǒng)的空間分辨率和測量精度,從而提高相控陣聲學多普勒流速剖面儀的性能。

聲學多普勒流速剖面儀;相控陣;波束形成;編碼信號

0 引言

聲學多普勒流速剖面儀(Acoustic Doppler Current Profile, ADCP)是一種利用聲學換能器發(fā)射的聲脈沖在隨水流運動的懸浮物質(zhì)中產(chǎn)生的多普勒頻移而進行流速測量的儀器。其突出特點是能夠測量不同水層的三維流速和流向,流速剖面范圍能夠達到上千米。聲學多普勒流速剖面儀系統(tǒng)的信號可分為窄帶和寬帶兩類。窄帶聲學多普勒流速剖面儀(Narrowband Acoustic Doppler Current Profile, NBADCP)發(fā)射單頻脈沖信號,優(yōu)點是流速剖面范圍較大;缺點是流速測量標準方差較大,時間和空間分辨率較低。寬帶聲學多普勒流速剖面儀(Broadband Acoustic Doppler Current Profile, BBADCP)發(fā)射由重復(fù)二進制偽隨機相位編碼調(diào)制的正弦脈沖信號,接收來自不同水層的散射回波信號,計算多普勒頻偏所在時間間隔內(nèi)的相位變化,進而獲得沿波束方向的流速。BBADCP采用脈沖編碼技術(shù)解決了窄帶ADCP空間分辨率和時間分辨率的矛盾,提高了分層測流的精度。

隨著要求流速測量范圍的增加,聲學多普勒流速剖面儀的工作頻率降低,換能器基陣的面積也急劇增大。如果仍然采用傳統(tǒng)的活塞式換能器基陣,則換能器基陣尺寸會大到無法忍受的地步。以38 kHz的ADCP為例,它的最大作用距離可以達到1000 m[1]。如果采用4°的波束開角,單個換能器的直徑會達到0.6 m。如果采用傳統(tǒng)的Janus配置的四波束換能器基陣,則整個換能器基陣的直徑會達到1.5 m的量級,遠遠超過了一般艦船所能承受的安裝尺寸。考慮到換能器基陣的適裝性,必須要減小換能器基陣的尺寸。

近年來,相控陣技術(shù)獲得了比較多的應(yīng)用。相控陣原理就是采用波束形成的方法,對多元陣陣元的發(fā)射或接收信號進行時延或相移補償,從而獲取預(yù)定方向波束的方法。相控陣ADCP采用的波束形成的方法可獲取Janus配置的四個收發(fā)合置的波束。相對于傳統(tǒng)的四波束活塞式換能器,在相同波束寬度和工作頻率的情況下,相控陣換能器能夠大大減小尺寸。由于相控陣換能器本身的物理特性,它還具有如下優(yōu)點:相控陣ADCP不需要進行聲速的補償和修正,大大簡化了傳統(tǒng)ADCP換能器陣利用聲透鏡等方法進行聲速補償?shù)姆爆嵅襟E;相控陣ADCP一般為圓形平面陣,很容易與載體共形,不需導(dǎo)流罩來抑制氣泡的產(chǎn)生;相控陣由多陣元組合而成,通過波束形成技術(shù)獲取Janus配置的四個波束,很容易保證波束指向、波束寬度、發(fā)射響應(yīng)、接收靈敏度等聲學特性參數(shù)的一致性。相控陣的這些優(yōu)點,使它在低頻、大量程測速方面獲得了廣泛的應(yīng)用。

常規(guī)相控陣聲學多普勒流速剖面儀采用窄帶信號體制,利用相移波束形成的方法形成發(fā)射和接收波束。由于信號相位的周期性,相移波束形成可在一定程度上大大簡化控制結(jié)構(gòu),所以窄帶相控陣波束形成易于實現(xiàn)。但是對寬帶系統(tǒng)來說,信號帶寬大大寬于窄帶系統(tǒng),又由于相移是頻率的函數(shù),寬帶系統(tǒng)的相移隨著頻率的改變而改變,所以寬帶系統(tǒng)一般情況下無法使用窄帶系統(tǒng)所用的相移波束形成方法。時延波束形成可做到與頻率無關(guān),所以寬帶系統(tǒng)一般采用時延波束形成方法。時延波束形成給寬帶系統(tǒng)帶來的問題就是波束控制結(jié)構(gòu)太復(fù)雜,每路子陣列均需進行時延控制,系統(tǒng)硬件復(fù)雜度大大增加。

本文在分析相移波束形成和時延波束形成基礎(chǔ)上,介紹了窄帶相控陣ADCP相移波束形成和寬帶相控陣ADCP時延波束形成的原理。基于相移波束形成,計算隨著頻率的偏移,波束形成的波束角與頻率變化的關(guān)系,從而獲取帶寬的偏移對波束寬度帶來的影響。然后從作用距離、分層厚度和測量精度以及系統(tǒng)復(fù)雜度三個方面分析利用相移波束形成發(fā)射編碼信號帶來的影響。由分析可知,窄帶相控陣聲學多普勒流速剖面儀可進行編碼信號的發(fā)射和接收,能夠提高系統(tǒng)的空間分辨率和測量精度,從而提高相控陣聲學多普勒流速剖面儀的性能。

1 相控陣的相移和時延波束形成

1.1 相控陣ADCP的波束形成

以TRDI公司的75 kHz平面相控陣ADCP為例,分析多元圓形平面陣的波束形成方法。相控陣發(fā)射信號中心頻率=75 kHz,取水中聲速=1500 m/s,則波長=20 mm,陣元間距=/2=10 mm,相控陣總直徑為480 mm,陣元數(shù)為988[2]。如圖1所示,黑色陣元產(chǎn)生左右兩個波束,白色陣元產(chǎn)生前后兩個波束,構(gòu)成Janus配置的ADCP的四個波束,如圖2所示[3]。由于白色陣元相當于黑色陣元旋轉(zhuǎn)90°而構(gòu)成,兩種顏色的陣元具有相似的性質(zhì),所以這里只針對黑色陣元產(chǎn)生的左右兩個波束進行分析。黑色陣元共有36列,組成一個36元子陣的線陣,每個子陣隨包含陣元數(shù)目的不同而具有不同的靈敏度。36列各子陣列輸出信號為:

式中:M為子陣列的接收靈敏度,隨子陣包含陣元數(shù)目的不同而不同;為信號角頻率;為相鄰陣元接收信號間的相位差,有

(2)

式中:為相鄰陣元間的信號時延;為波束角。則相控陣的總輸出可以表示為

1.2 相控陣的相移波束形成

一個元等間距線陣的歸一化自然指向性函數(shù)為[4]

在陣元之間插入相移使波束主極大方向控制到特定的方位的方法稱為相移波束形成。在窄帶應(yīng)用(一般在主動聲吶中)時,常用相移波束形成。在陣元之間插入相移后的元等間距線陣的歸一化自然指向性函數(shù)為

(5)

即只要相鄰陣元間均插入相移后,即波束主極大方向滿足

將式(6)代入式(2)可知,波束主極大方向為

(7)

由式(7)可知,相移波束形成下的波束主極大方向是頻率的函數(shù)。只有在窄帶情況下(相控陣信號系統(tǒng)帶寬很窄,僅收發(fā)單頻脈寬信號),才能假定波束主極大方向0與插入相移的關(guān)系不受頻率的影響。如圖3所示,當平面聲波從偏離基陣法線0= 30°方向入射,且陣元間距=/2時,相鄰兩列陣元的相位差為

如圖1所示的相控陣ADCP陣元排列模型中,每一列陣元作為一路子陣列。由于0、4、8、12…列陣元,1、5、9…列陣元,2、6、10…列陣元,3、7、11…列陣元同相位,各陣列間可以分別并聯(lián)起來,構(gòu)成平面陣的四路輸出線。改變四路輸出線之間的相移角,便可以構(gòu)成窄帶相移相控陣的發(fā)射和接收波束。這里只針對黑色陣元產(chǎn)生的左右兩個波束進行分析,接收波束1的相移波束形成如圖4所示,接收波束2的相移波束形成如圖5所示,發(fā)射波束的相移波束形成如圖6所示。相移波束形成使窄帶相控陣ADCP的相控波束形成大大簡化,只需要對四路波束進行相移控制,便能夠同時形成兩路發(fā)射波束和單獨的兩路接收波束。

圖3 線陣波束形成

Fig.3 Beamformer of linear array

由于信號相位的周期性,相移波束形成可在一定程度上大大簡化硬件控制結(jié)構(gòu),所以窄帶相控陣波束形成易于實現(xiàn),并且被廣泛應(yīng)用于聲學多普勒測速測流的產(chǎn)品中。

1.3 相控陣的時延波束形成

相鄰陣元間的插入相移與時延的關(guān)系為=2π,代入式(7),則波束主極大方向?qū)懗蓵r間延遲形式為

由式(9)可知,時延波束形成下的波束主極大方向與頻率無關(guān)。一般在寬帶情況下(相控陣信號系統(tǒng)帶寬很寬,甚至占到中心頻率的25%帶寬),可以使用時延波束形成的方法進行波束形成。

如圖1所示的相控陣ADCP陣元排列模型中,每一列陣元作為一路子陣列。如圖3所示,當聲波從0方向入射時:第1列陣元的輸出較第0列陣元的時延為=sin0/,第2列陣元的輸出較第1列陣元的時延為=sin0/,……,即相鄰陣列之間均需要插入相同的時延值。接收波束的時延波束形成如圖7所示,發(fā)射波束的時延波束形成如圖8所示。時延波束形成的優(yōu)點是波束主極大方向只與時延值有關(guān),通過改變時延值能夠靈活地改變波束主極大方向;缺點是增加了波束控制的復(fù)雜度,由于時延沒有周期性的特點,每個相鄰陣列間均需進行時延控制,波束控制的路數(shù)要遠大于相移波束形成方法的路數(shù)。

2 編碼信號的相移波束形成

2.1 編碼信號

寬帶聲學多普勒流速剖面儀采用重復(fù)二進制偽隨機相位編碼調(diào)制的正弦脈沖信號,編碼一般采用巴克碼、m序列等編碼方式[5]。編碼信號在時域表現(xiàn)為信號相位的[0,]反轉(zhuǎn),在頻域表現(xiàn)為信號能量帶寬的拓展,編碼信號的-3 dB寬度由每個編碼碼元填充正弦波周期數(shù)確定。如果碼元填充數(shù)為2,則編碼信號-3 dB帶寬為中心頻率的50%。圖9為正弦波由二進制偽隨機編碼相位調(diào)制的示意圖。編碼信號的帶寬隨著碼元正弦波填充數(shù)的增加而減小,當碼元填充數(shù)為10時,編碼信號-3 dB帶寬為中心頻率的10%。

2.2 編碼信號對相移波束形成的影響

相移波束形成的方法無法應(yīng)用在寬帶情況下,主要是由于隨著信號帶寬的增大,相對中心頻率產(chǎn)生偏移的頻段,形成波束的主極大方向也隨之發(fā)生偏移,使寬帶信號的波束寬度增大,能量分散,從而影響系統(tǒng)性能。為了估算編碼信號的波束寬度對系統(tǒng)性能的影響,計算頻率偏移與波束主極大方向的關(guān)系。

為了減小帶寬對波束形成的影響,假設(shè)編碼信號系統(tǒng)的帶寬中心頻率比為10%。對于75 kHz中心頻率的相控陣ADCP系統(tǒng)來說,則偏離中心頻率的上下限分別為±3.75 kHz。對相移波束形成系統(tǒng)來說,插入相移的元等間距線陣的歸一化自然指向性函數(shù)如式(5)所示[6]。利用相移波束形成,+30°方向接收波束形成如圖10所示。陣元間距與波長之比為0.5,陣元數(shù)目為36,波束角E為30°,波束開角=3.27°,-3 dB角度分別為28.38°和31.65°。

為了進行編碼信號發(fā)射,對頻率上下限±3.75 kHz的波束形成角度進行計算。當頻率升高5%時,波長減小,陣元間距不變,假設(shè)陣元間的相移仍然為90°,則波束形成如圖11所示,波束角E為28.445°,波束開角=3.07°,-3 dB角度分別為26.91°和29.98°。當頻率降低5%時,波長增大,則波束形成如圖12所示,波束角E為31.775°。波束開角=3.5°,-3 dB角度分別為30.02°和33.53°。因此,對帶寬10%的編碼信號來說,相控陣換能器的波束開角變?yōu)?3.53°-26.91°=6.62°,波束開角大約增大到窄帶信號波束開角的2倍。波束開角的增大引起換能器聚焦系數(shù)減小,發(fā)射聲源級減小,系統(tǒng)作用距離減小,但是由于發(fā)射信號變?yōu)閭坞S機編碼的編碼信號,相控陣系統(tǒng)的測量精度得到提高。

考慮到換能器的面積和發(fā)射信號的能量聚集程度,聲學多普勒流速剖面儀的波束寬度一般取為3°~4°。采用相移波束形成的方法進行限定帶寬的編碼信號收發(fā),僅僅增加信號收發(fā)的波束寬度,減小發(fā)射信號的能量聚集程度,優(yōu)點是能夠大大簡化波束控制復(fù)雜度,減少信號控制流程。

3 編碼信號相移波束形成的性能分析

3.1 作用距離

由第2節(jié)中計算可知,利用相移波束形成的方法進行編碼信號發(fā)射,使信號波束寬度增大,從而造成聲吶系統(tǒng)聲源級降低,減小了系統(tǒng)的作用距離。相控陣ADCP為主動聲吶的一種,聲吶方程可以表示為[7]:

其中:為體積散射強度,為聲源級,為傳播損失,為體積散射強度,為散射體積,為接收靈敏度。

相控陣ADCP發(fā)射10%帶寬的編碼信號,波束開角為發(fā)射窄帶信號的2倍。換能器的聚集系數(shù)減小6 dB,則聲源級和接收靈敏度均減小6 dB;照射體積增大一倍,則散射體積對應(yīng)分貝值增加6 dB;總回波強度減小6 dB。對75 kHz相控陣ADCP來說,如果均采用相移波束形成的方法,10%帶寬的編碼信號系統(tǒng)比窄帶信號系統(tǒng)的作用距離大概減小20%。

3.2 分層厚度和測量精度

窄帶ADCP系統(tǒng)中,由分層厚度決定發(fā)射信號脈寬,每層得到一個速度測量值。而在編碼ADCP系統(tǒng)中,編碼信號(如圖9所示)采用復(fù)相關(guān)算法進行計算,每一個重復(fù)的編碼信號都可以計算得到一個速度值,編碼的個數(shù)對應(yīng)計算得到的速度個數(shù),則每層測量得到的速度值個數(shù)可以由碼元數(shù)所確定,碼元寬度對應(yīng)編碼信號分層厚度,即編碼信號ADCP的分層厚度可遠小于窄帶ADCP的分層厚度。

多普勒測頻為隨機過程,窄帶ADCP系統(tǒng)和編碼信號ADCP系統(tǒng)分別采用單頻脈寬信號和二進制偽隨機編碼信號作為發(fā)射信號,由于信號體制的不同,兩種系統(tǒng)的測量精度也不相同。窄帶ADCP系統(tǒng)速度測量的標準方差為[5]

其中:為水中聲速;為波長;為分層厚度。

其中:為編碼和信號處理的非理想程度,一般取為1.5,M為碼元數(shù),為相關(guān)系數(shù)幅值,一般取0.5。

利用相控陣ADCP的相移波束形成,相控陣換能器的-3 dB帶寬大于中心頻率的10%,發(fā)射編碼信號,編碼信號碼元填充數(shù)為10,編碼信號的-3 dB帶寬為中心頻率的10%。當編碼長度為11位時,由式(11)和式(12)可知,編碼信號ADCP系統(tǒng)速度測量的標準方差是窄帶系統(tǒng)標準方差的40%,即精度提高了2.5倍。相對于發(fā)射單頻脈寬信號的窄帶系統(tǒng),發(fā)射編碼信號的ADCP系統(tǒng)能夠提高系統(tǒng)的測量精度性能。

3.3 系統(tǒng)復(fù)雜度

由第2.1節(jié)的分析可知,對于988個陣元的75 kHz相控陣ADCP,黑白陣元共構(gòu)成72路子陣列。采用相移波束形成的相控陣ADCP中,72路子陣列可以根據(jù)相移規(guī)律并成8路子陣,其中4路子陣形成左右兩個波束,另外4路子陣形成前后兩個波束,從而構(gòu)成Janus配置的四個波束。對8路相移通道進行控制,便可采用一套發(fā)射機和四路接收機實現(xiàn)相控陣ADCP的相移波束形成。

采用時延波束形成的相控陣ADCP中,由于時延信號特征,72路子陣列無法簡并。其中36路子陣列形成左右兩個波束,另外36路子陣列形成前后兩個波束,從而構(gòu)成Janus配置的四個波束。需要對72路子陣列分別進行時延控制,相對于相移波束形成系統(tǒng),增加了系統(tǒng)控制電路的復(fù)雜度,加大了系統(tǒng)的硬件開銷。

4 結(jié)論

在不改變硬件電路結(jié)構(gòu)的前提下,相控陣ADCP系統(tǒng)利用窄帶相移波束形成技術(shù)進行編碼信號的發(fā)射與接收,提高了系統(tǒng)的空間分辨率和測量精度(單呯測速標準方差降為窄帶系統(tǒng)的40%),從而提高了設(shè)備的性能;缺點是由于波束開角的增加,系統(tǒng)的作用距離減小。對比相控陣系統(tǒng)發(fā)射窄帶信號和編碼信號的優(yōu)缺點,在損失作用距離的情況下提高了空間分辨率和測量精度,滿足了流速測量主要指標的提高。下一步可以考慮進行系統(tǒng)試驗,驗證波束寬度的增大對系統(tǒng)的實際影響,為系統(tǒng)的實用化奠定基礎(chǔ)。

[1] Ocean Surveyor Datasheet. Teledyne RD Instruments. 2009

[2] 劉國勤, 彭東立. 相控陣多普勒計程儀的相控波束形成[J]. 聲學技術(shù), 2015, 34(4): 199-203.

LIU Guoqing, PENG Dongli. Beam former of Phased-Array Doppler Velocity Log[J]. Technical Acoustics, 2015, 34(4): 199-203.

[3] Yü Xiaolong, Steve E Bradley, Francis D Rowe. Two-dimensional array transducer and beamformer, US Patent No.5808967, Sep.15, 1998.

[4] 田坦. 聲吶技術(shù)(第二版)[M]. 哈爾濱: 哈爾濱工程大學出版社, 2009.

TIAN Tan. Technical of Sonar[M]. 2nd edition. Harbin: Harbin Engineering University Press, 2009.

[5] Blair H Brumley, Ramon G Cabrera, Kent L Deines, et al. Performance of a broad-band acoustic doppler current profile[J]. IEEE Journal of Oceanic Engineering, 1991, 16(4): 402-407.

[6] 田坦, 張殿倫. 相控陣多普勒測速技術(shù)研究[J]. 哈爾濱工程大學學報, 2002.

TIAN Tan, ZHANG Dianlun. Study of phased_array doppler velocity measurement technique[J]. Journal of Harbin Engineering University, 2002.

[7] Urick R J著. 洪申譯. 水聲原理[M]. 哈爾濱: 哈爾濱船舶工程學院出版社, 1990.

Urick R J. Principles of underwater sound[M]. 3rd edition, Harbin: Harbin Engineering University Press, 1990.

Coded signal phase shift beamforming of phased-array ADCP

LI Peng-fei1, PENG Dong-li2

(1. Naval equipment academy, Department of warship and submarine,Beijing 100841,China;2.Shanghai Acoustic Laboratory, Institute of Acoustics,Chinese Academy of Sciences,Shanghai 200032,China)

The volume of the transducer of phased-array ADCP(Acoustic Doppler Current Profiler) is much smaller than the volume of traditional piston ADCP, and compensation of sound velocity is not needed for physical characteristic of phased-array ADCP. The narrowband phased-array ADCP can transmit and receive coded signal to improve the equipment performance. This article starts with the beam forming method of phased-array ADCP, and then describes the theory of phase shift beamforming and time delay beamforming respectively. The beam angles of phased array ADCP are compared, and the working distance, measurement accuracy and system complexity are analyzed. The conclusion indicates that the narrowband Phased-array ADCP can transmit and receive coded signal to improve spatial resolution and measurement precision, and the performance of phased array ADCP is highly improved.

Acoustic Doppler Current Profiler(ADCP); phased-array; beam forming;coded signal

U666.7

A

1000-3630(2015)-06-0545-06

10.16300/j.cnki.1000-3630.2015.06.015

2015-07-20;

2015-10-19

李鵬飛(1978-), 男, 黑龍江哈爾濱人, 工程師, 研究方向為檢測技術(shù)。

彭東立, E-mail: bernou@163.com

猜你喜歡
信號系統(tǒng)
Smartflower POP 一體式光伏系統(tǒng)
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
WJ-700無人機系統(tǒng)
ZC系列無人機遙感系統(tǒng)
北京測繪(2020年12期)2020-12-29 01:33:58
完形填空二則
基于PowerPC+FPGA顯示系統(tǒng)
半沸制皂系統(tǒng)(下)
孩子停止長個的信號
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
基于LabVIEW的力加載信號采集與PID控制
主站蜘蛛池模板: …亚洲 欧洲 另类 春色| 亚洲天堂日韩在线| 日韩欧美中文字幕在线韩免费| 日韩久久精品无码aV| 秘书高跟黑色丝袜国产91在线| 中文字幕永久在线看| 欧美视频在线播放观看免费福利资源 | 99在线观看免费视频| 国产成人做受免费视频| 日韩欧美高清视频| 中文字幕日韩丝袜一区| 超清无码熟妇人妻AV在线绿巨人| 国产免费自拍视频| 色噜噜综合网| 婷婷伊人五月| 久久精品人妻中文系列| 久久久噜噜噜久久中文字幕色伊伊| 欧美日韩中文字幕在线| 久久综合九九亚洲一区| 人人爽人人爽人人片| 国产成人免费手机在线观看视频 | 日韩欧美中文字幕在线精品| 国产高潮流白浆视频| 亚洲精品麻豆| 国产在线一二三区| 精品人妻系列无码专区久久| 手机在线国产精品| 四虎国产永久在线观看| 99re在线免费视频| 亚洲欧美成人网| 一区二区无码在线视频| 成人亚洲天堂| 国产综合精品日本亚洲777| 97国产在线视频| 亚洲精品制服丝袜二区| 91免费精品国偷自产在线在线| 四虎在线观看视频高清无码 | 久久亚洲天堂| 亚洲视频二| 欧美 国产 人人视频| 亚洲AV无码一区二区三区牲色| 久久77777| 亚洲欧洲美色一区二区三区| 色哟哟国产精品| 亚洲床戏一区| 中文字幕首页系列人妻| 黄色三级网站免费| 国模极品一区二区三区| 四虎成人精品在永久免费| 首页亚洲国产丝袜长腿综合| 欧美翘臀一区二区三区| 久久这里只精品热免费99| 国产人人乐人人爱| 亚洲欧美国产视频| 秋霞午夜国产精品成人片| 国内老司机精品视频在线播出| 91系列在线观看| www精品久久| 久久国产高潮流白浆免费观看| av在线手机播放| 日韩欧美中文亚洲高清在线| 中文字幕在线永久在线视频2020| 亚洲最新地址| 亚洲综合香蕉| 国产精品99在线观看| 国产精品亚洲一区二区三区在线观看| 国产白浆视频| 国产麻豆精品久久一二三| 国产一区二区精品福利 | 美女国产在线| 无码专区国产精品一区| 美女一区二区在线观看| 国产精品开放后亚洲| 亚洲欧美日韩中文字幕在线一区| 亚洲无码A视频在线| 国产毛片基地| 久久一色本道亚洲| 99这里精品| 国产男女免费视频| 日本影院一区| 中文字幕 日韩 欧美| 99久久国产综合精品女同|