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

    深海海底反射會(huì)聚區(qū)聲傳播特性?

    2019-01-25 09:53:28張鵬李整林吳立新張仁和秦繼興
    物理學(xué)報(bào) 2019年1期
    關(guān)鍵詞:聲線聲速聲源

    張鵬 李整林 吳立新 張仁和 秦繼興

    1)(中國科學(xué)院聲學(xué)研究所,聲場(chǎng)聲信息國家重點(diǎn)實(shí)驗(yàn)室,北京 100190)

    2)(中國科學(xué)院大學(xué)電子電氣與通信工程學(xué)院,北京 100049)

    (2018年9月26日收到;2018年11月9日收到修改稿)

    在深海聲道條件下,海水折射效應(yīng)會(huì)使得聲場(chǎng)出現(xiàn)會(huì)聚效應(yīng);在不完全聲道條件下,深海海底對(duì)聲場(chǎng)具有重要影響.利用在中國南海海域收集到的一次深海聲傳播實(shí)驗(yàn)數(shù)據(jù),研究了深海不完全聲道環(huán)境下的海底反射對(duì)聲傳播的影響.實(shí)驗(yàn)觀測(cè)到不同于深海會(huì)聚區(qū)的海底反射會(huì)聚現(xiàn)象,在直達(dá)聲區(qū)范圍內(nèi)的海底地形隆起可導(dǎo)致海底反射會(huì)聚區(qū)提前形成,并使得部分影區(qū)的聲強(qiáng)明顯提高.由于不平坦海底和海面的反射破壞了完全聲道環(huán)境下的會(huì)聚區(qū)結(jié)構(gòu),在60 km范圍內(nèi)存在兩個(gè)海底反射會(huì)聚區(qū),會(huì)聚區(qū)增益可達(dá)10 dB以上,同時(shí)在11 km附近的影區(qū)和51 km附近形成高聲強(qiáng)區(qū)域.當(dāng)接收深度與聲源深度相同時(shí),第二會(huì)聚區(qū)的增益高于第一會(huì)聚區(qū).在第一會(huì)聚區(qū)內(nèi),隨著接收深度的增加,聲線到達(dá)結(jié)構(gòu)趨于復(fù)雜,多途效應(yīng)更加明顯.使用拋物方程數(shù)值分析結(jié)合射線理論對(duì)深海海底反射會(huì)聚區(qū)現(xiàn)象產(chǎn)生的物理原因進(jìn)行了分析解釋.研究結(jié)果對(duì)于聲納在深海復(fù)雜環(huán)境下的性能分析具有重要的指導(dǎo)意義.

    1 引 言

    聲速剖面隨深度的變化而導(dǎo)致的聲場(chǎng)會(huì)聚區(qū)效應(yīng)是深海海區(qū)一個(gè)重要的水聲環(huán)境特性[1],利用聲傳播的會(huì)聚區(qū)效應(yīng),可以更好地實(shí)現(xiàn)遠(yuǎn)程水聲通信和探測(cè).低頻聲信號(hào)的會(huì)聚區(qū)傳播對(duì)實(shí)現(xiàn)遠(yuǎn)程水聲通信系統(tǒng)和潛艇低頻噪聲的探測(cè)性能分析和預(yù)報(bào)具有重要的意義.

    Hale[2]很久以前就在海上實(shí)驗(yàn)中觀測(cè)到很強(qiáng)的會(huì)聚區(qū)效應(yīng),并對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行了詳細(xì)地研究.Urick[3]使用射線理論分析了會(huì)聚區(qū)隨聲源深度變化的情況.布列霍夫斯基赫等[4]則從波動(dòng)理論出發(fā),給出了會(huì)聚區(qū)焦散線的聚焦因子公式.Williams和Horul[5]應(yīng)用簡(jiǎn)正波理論研究海洋聲場(chǎng)的會(huì)聚現(xiàn)象,指出會(huì)聚區(qū)是由大量同相簡(jiǎn)正波的疊加形成的.張仁和[6]在聲場(chǎng)簡(jiǎn)正波理論的基礎(chǔ)上,利用廣義相積分近似得出了會(huì)聚區(qū)增益的解析表達(dá)式,解決了經(jīng)典射線理論在反轉(zhuǎn)點(diǎn)會(huì)聚區(qū)發(fā)散的困難.這些理論都深刻地揭示了海洋聲傳播的規(guī)律,能夠較為準(zhǔn)確地預(yù)報(bào)會(huì)聚區(qū)的位置,但他們的討論都著重于聲波在典型深海SOFAR聲道中的傳播問題.聲波在SOFAR聲道中傳播時(shí),大部分聲線在未觸及海底之前就發(fā)生反轉(zhuǎn),因而海底的影響可以忽略.在深海不完全聲道中,當(dāng)海底處的海水聲速小于海面附近聲源處的海水聲速時(shí),其聲場(chǎng)主要由海底反射聲波所形成,這個(gè)時(shí)候海底的反射損失和海底剖面成為影響海洋聲傳播的重要因素,有必要仔細(xì)考慮海底邊界的聲學(xué)特性.

    張仁和等[7]應(yīng)用反轉(zhuǎn)點(diǎn)會(huì)聚區(qū)理論,引入最小掠射角的海底反射系數(shù)因子,計(jì)算討論了深海負(fù)梯度剖面下海底反射形成的會(huì)聚區(qū)的問題,并從實(shí)驗(yàn)數(shù)據(jù)中得到最小掠射角的海底反射系數(shù).但是該論文中并沒有直接描述海底的聲學(xué)特性,理論計(jì)算的會(huì)聚區(qū)位置和增益與實(shí)驗(yàn)結(jié)果也存在一定的誤差.王剛[8]應(yīng)用WKBZ簡(jiǎn)正波方法,對(duì)深海負(fù)梯度聲速剖面下的海底反射會(huì)聚區(qū)進(jìn)行了研究.范培勤等[9]利用海洋環(huán)境統(tǒng)計(jì)數(shù)據(jù)與反射會(huì)聚區(qū)特征參數(shù)計(jì)算模型,仿真實(shí)現(xiàn)了海底反射會(huì)聚區(qū)距離的快速預(yù)報(bào),但是缺乏相關(guān)實(shí)驗(yàn)數(shù)據(jù)的支持.針對(duì)海底地形變化以及沉積層的特性對(duì)聲傳播效應(yīng)的影響,許多水聲學(xué)者進(jìn)行了較為深入的研究[10-21].李文等[10]對(duì)深海海底山二維聲傳播規(guī)律進(jìn)行了研究,由于海底山的反射遮擋效應(yīng),深海會(huì)聚區(qū)結(jié)構(gòu)遭到破壞,傳播損失增大.秦繼興等[11]利用大陸坡海域的實(shí)驗(yàn)數(shù)據(jù),分析了下坡時(shí)聲源位于海面附近和淺海斜坡表面時(shí)斜坡對(duì)聲傳播的影響.胡治國等[12]、楊坤德等[13]研究了深海海底斜坡的復(fù)雜地形下的聲傳播規(guī)律,分析了不平坦海底引起的傳播損失異常的機(jī)理.

    本文利用一次南海深海聲傳播實(shí)驗(yàn)數(shù)據(jù),分析了存在復(fù)雜地形的不完全聲道條件下海底反射會(huì)聚區(qū)的聲傳播規(guī)律,通過射線模型分析脈沖到達(dá)結(jié)構(gòu)與本征聲線到達(dá)時(shí)間,解釋海底反射會(huì)聚現(xiàn)象.

    2 深海不完全聲道條件下聲傳播

    2018年4月,中科院聲學(xué)所聲場(chǎng)聲信息國家重點(diǎn)實(shí)驗(yàn)室在中國南海進(jìn)行了一次綜合性的海上實(shí)驗(yàn),其中主要內(nèi)容之一是研究深海復(fù)雜環(huán)境下的聲傳播特性.海上實(shí)驗(yàn)設(shè)備布放示意圖如圖1所示,采用的是實(shí)驗(yàn)船結(jié)合深海聲學(xué)接收潛標(biāo)的作業(yè)方式進(jìn)行,其中接收潛標(biāo)系統(tǒng)由20個(gè)自容式水聽器(USR)組成,USR以非等間距方式布放在85—3400 m深度范圍內(nèi),接收靈敏度為-170 dB,信號(hào)的采樣率為16 kHz.實(shí)驗(yàn)1號(hào)科考船以4節(jié)的航速拖曳發(fā)射換能器,并發(fā)射頻帶范圍為250—350 Hz的雙曲調(diào)頻(hyperbolic frequency modulated,HFM)信號(hào),拖曳換能器的深度在150 m左右.發(fā)射信號(hào)的時(shí)序如圖2所示,每次發(fā)射信號(hào)時(shí)間長(zhǎng)度為20 s,兩個(gè)信號(hào)間隔20 s,連續(xù)發(fā)射4次信號(hào)后再空50 s,一輪信號(hào)總共持續(xù)時(shí)間190 s.發(fā)射聲源級(jí)為192 dB.

    圖3給出了聲傳播實(shí)驗(yàn)路徑上距離接收潛標(biāo)180 km范圍內(nèi)的海深.由圖3可見,傳播路徑上的海底地形變化較為復(fù)雜,20 km內(nèi)有一個(gè)高度500 m的海底小山,海底地形起伏較大,20—80 km距離上海深變化較小,地形整體較為平坦,80 km處海深急劇減小到2800 m左右,80 km以后的海深穩(wěn)定在2700 m附近.實(shí)驗(yàn)期間使用拋棄式溫鹽深探頭(conductivity temperature depth,CTD)和拋棄式溫度探頭(expendable bathy thermograph,XBT)每隔10 km對(duì)海水聲速剖面進(jìn)行了一次測(cè)量,不同距離上的聲速剖面變化較小,其中接收位置處的聲速剖面如圖4所示.圖4中黑色圓點(diǎn)表示接收陣的20個(gè)水聽器布放深度.實(shí)測(cè)的聲速剖面最大深度為1200 m,由于聲道軸以下的聲速剖面隨時(shí)間變化較為穩(wěn)定,因此1200 m以下未測(cè)量深度的聲速剖面從數(shù)據(jù)庫查詢獲得,將其與實(shí)測(cè)聲速插值連接得到完整深度的聲速剖面.由圖4可見,實(shí)驗(yàn)測(cè)線海域?yàn)檩^典型的深海不完全聲道,聲道軸位于1100 m深度附近,海底最深處聲速為1515 m/s,小于海面處海水聲速1539 m/s.圖4中藍(lán)色虛線表示拖曳聲源所在深度150 m,藍(lán)色實(shí)線表示聲源深度處的聲速,與海底處的海水聲速接近,但是在10 km以后海底位置的聲速小于聲源位置處的聲速.因此,海底對(duì)于聲波的反射作用不可忽視.

    圖1 聲傳播實(shí)驗(yàn)設(shè)備布放示意圖Fig.1.Experimental con figuration.

    圖2 發(fā)射信號(hào)序列示意圖Fig.2.The acoustic emission signal during the experiment.

    圖3 聲傳播路徑上的海深變化Fig.3.The bathymetry along the propagation track.

    圖4 實(shí)驗(yàn)中測(cè)量的海水聲速剖面Fig.4.Sound speed pro file during the experiment.

    為了獲得海底聲學(xué)特性,實(shí)驗(yàn)中在接收潛標(biāo)位置附近進(jìn)行了海底底質(zhì)采樣,采集的海底樣品總長(zhǎng)度為190 cm.采樣分析結(jié)果如表1所列,底質(zhì)類型為粉砂質(zhì)黏土,海底表層處沉積層的聲速和密度分布較為均勻,平均聲速為1499 m/s,平均密度為1.36 g/cm3,底質(zhì)采樣測(cè)量結(jié)果為海底參數(shù)的選取提供了重要參考依據(jù).

    表1 潛標(biāo)位置附近的海底底質(zhì)采樣Table 1.Seabed sampling near the receiving array.

    將實(shí)驗(yàn)中接收到的聲信號(hào)記為x(t),對(duì)其進(jìn)行脈沖壓縮來提高信噪比,得到脈壓信號(hào)xc(t),再經(jīng)過離散傅里葉變換得到信號(hào)的頻譜Xi,并且在中心頻率f0=300 Hz的實(shí)驗(yàn)發(fā)射信號(hào)帶寬內(nèi)(100 Hz)求取平均,得到窄帶信號(hào)的平均能量為

    式中,Fs為采樣率,f1和f2為發(fā)射信號(hào)頻率的上下限,nf1和nf2分別為頻率上下限對(duì)應(yīng)的頻點(diǎn)數(shù).

    實(shí)驗(yàn)的傳播損失可以表示為

    式中,SL(f0)為發(fā)射換能器聲源級(jí),Mv為接收水聽器的靈敏度,EC為根據(jù)已知的發(fā)射聲源信號(hào)計(jì)算得到的脈沖壓縮時(shí)頻增益.

    最后,得到傳播路徑上實(shí)驗(yàn)的二維傳播損失如圖5(a)所示.從圖5(a)中可以看到,在60 km距離內(nèi)存在兩個(gè)明顯的會(huì)聚區(qū)結(jié)構(gòu),而60 km以外的會(huì)聚區(qū)結(jié)構(gòu)并不明顯,因此我們選取前60 km距離上的傳播損失(圖5(b)).在60 km距離內(nèi)可以觀察到兩個(gè)明顯的會(huì)聚區(qū)結(jié)構(gòu),它們分別位于20和40 km的距離附近.同時(shí),在11和51 km附近可以觀測(cè)到兩個(gè)較明顯的聲增強(qiáng)區(qū)域.此外需要說明:圖5中給出了整個(gè)深度的傳播損失結(jié)果,由于實(shí)驗(yàn)中在較大深度上水聽器布放間距較大,所以繪圖時(shí)兩個(gè)相鄰接收深度中間插值得到的傳播損失結(jié)果與實(shí)際聲場(chǎng)分布情況有偏差,但在有水聽器的接收深度上,實(shí)驗(yàn)傳播損失結(jié)果可以反映實(shí)際聲場(chǎng)分布情況.

    為了分析圖5所示的實(shí)驗(yàn)現(xiàn)象,使用了拋物方程近似聲場(chǎng)模型(range-dependent acoustic modelparabolic equation,RAM-PE[22,23])計(jì)算海深水平變化環(huán)境下的聲傳播損失,選擇與實(shí)驗(yàn)換能器發(fā)射信號(hào)相同的聲源中心頻率和帶寬,頻率間隔10 Hz,頻點(diǎn)數(shù)為11.海底參數(shù)選取兩層海底模型,其中沉積層厚度為20 m,根據(jù)采樣測(cè)量結(jié)果,沉積層聲速和密度分別選為1499 m/s和1.36 g/cm3,無限大基底的聲速取為1650 m/s,密度為1.8 g/cm3,沉積層和基底的衰減系數(shù)取[24]

    選取兩種海底地形進(jìn)行對(duì)比,第一種地形是實(shí)驗(yàn)中同步測(cè)量的海底地形(圖3),第二種地形是接收位置處海深的平坦海底.根據(jù)聲學(xué)互易原理,理論計(jì)算時(shí)將聲源置于0 km距離處不同水聽器所在深度,在拖曳聲源深度的不同距離處接收信號(hào).由于實(shí)驗(yàn)中在聲道軸以下布設(shè)的水聽器較為稀疏,所以圖5中大深度位置直達(dá)聲區(qū)與影區(qū)之間的聲傳播損失會(huì)因插值出現(xiàn)不平滑,為避免出現(xiàn)這種現(xiàn)象,模型計(jì)算時(shí)聲道軸以下的接收深度增加,最終得到兩種地形下的二維聲傳播損失(圖6).

    對(duì)比圖5(b)和圖6(a)可以發(fā)現(xiàn),實(shí)測(cè)地形下仿真結(jié)果與實(shí)驗(yàn)數(shù)據(jù)處理得到的傳播損失較為符合.對(duì)比圖6中實(shí)測(cè)地形(a)和平坦地形(b)下的仿真結(jié)果,實(shí)測(cè)地形下在60 km內(nèi)存在兩個(gè)明顯的會(huì)聚區(qū)結(jié)構(gòu),聲線經(jīng)過海底小山丘的反射作用在20 km附近形成海底反射會(huì)聚區(qū),有一部分直達(dá)聲線和第一會(huì)聚區(qū)海水折射到達(dá)聲線一起經(jīng)海底反射后在40 km附近形成第二個(gè)海底反射會(huì)聚區(qū).而平坦海底環(huán)境下,部分折射路徑的聲線不與海底相互作用,經(jīng)過海水折射后到達(dá)50 km附近形成會(huì)聚區(qū),在10—40 km的影區(qū)內(nèi),海底反射使得聲能量較高,但其隨空間的分布較圖6(a)均勻.

    圖5 實(shí)驗(yàn)傳播損失隨距離和深度變化 (a)0—180 km;(b)0—60 kmFig.5.Two dimensional experimental TLs:(a)0–180 km;(b)0–60 km.

    圖6 RAM-PE模型計(jì)算的聲傳播損失 (a)實(shí)測(cè)海深;(b)平坦海底Fig.6.Numerical TLs results from RAM-PE model:(a)Real bathymetry;(b) flat bottom.

    為了定量分析其中的實(shí)驗(yàn)現(xiàn)象,圖7給出了四個(gè)典型深度(85,145,1008和3021 m)上的傳播損失對(duì)比.圖7結(jié)果可以進(jìn)一步證實(shí)實(shí)驗(yàn)結(jié)果與實(shí)測(cè)海底環(huán)境下的模型計(jì)算結(jié)果較為一致.實(shí)測(cè)地形下在20和40 km附近出現(xiàn)了兩個(gè)會(huì)聚區(qū),隨著接收深度增加,到1008 m深度處出現(xiàn)會(huì)聚區(qū)分裂現(xiàn)象.在本文中,將會(huì)聚區(qū)與其附近影區(qū)內(nèi)傳播損失的差值稱作會(huì)聚區(qū)增益.與第一會(huì)聚區(qū)相比,相同接收深度上第二會(huì)聚區(qū)的距離跨度更大,兩個(gè)會(huì)聚區(qū)的增益都在10 dB左右,在145 m深度上的第二會(huì)聚區(qū)增益高于其他接收深度的會(huì)聚區(qū)增益,甚至高于第一會(huì)聚區(qū).此外,在145 m深度上,距離11和51 km位置附近可觀察到兩個(gè)聲強(qiáng)增大的區(qū)域,其增益小于20和40 km處的會(huì)聚區(qū)增益,但是85 m接收深度上僅在11 km處存在聲增強(qiáng)區(qū)域,并沒有51 km處的聲強(qiáng)增大區(qū)域.

    在平坦地形下,85 m深度上沒有出現(xiàn)會(huì)聚區(qū),但是在10—40 km內(nèi)存在隨空間均勻分布的聲能量增強(qiáng)區(qū)域,這個(gè)區(qū)域同樣在145 m深度上出現(xiàn),而且145 m深度上在51 km處可觀測(cè)到一個(gè)明顯的會(huì)聚區(qū),其增益達(dá)到25 dB以上,遠(yuǎn)大于實(shí)驗(yàn)環(huán)境下的會(huì)聚區(qū)增益.在1008 m深度上會(huì)聚區(qū)發(fā)生分裂,但是仍具有較高的增益.在靠近海底的3021 m處,平坦地形下的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果在35 km以內(nèi)較為符合,直達(dá)聲區(qū)延伸至20 km附近,在31 km附近存在一個(gè)明顯的聲增強(qiáng)區(qū)域,增益到達(dá)15 dB以上,所以其傳播損失整體小于淺深度位置的結(jié)果.在35 km以后,隨著距離的增加,實(shí)驗(yàn)傳播損失逐漸減小,而平坦海底下聲傳播損失逐漸增大.

    圖7 四個(gè)不同接收深度實(shí)驗(yàn)傳播損失與模型計(jì)算結(jié)果比較 (a)接收深度85 m;(b)接收深度145 m;(c)接收深度1008 m;(d)接收深度3021 mFig.7.Comparison of experimental TLs and numerical TLs at four different receivers:(a)85 m;(b)145 m;(c)1008 m;(d)3021 m.

    3 海底反射會(huì)聚區(qū)現(xiàn)象的理論解釋

    為了解釋圖5和圖7中實(shí)驗(yàn)觀測(cè)到的一些特殊深度和距離條件下的海底反射會(huì)聚現(xiàn)象,我們利用射線進(jìn)行分析.根據(jù)互易原理,圖8給出了兩種海底地形情況下聲源深度在145 m處的聲線圖,圖中0 km處黑色圓點(diǎn)表示聲源位置,黑色的虛線表示對(duì)應(yīng)的150 m的互易接收深度.其中:藍(lán)綠色聲線表示經(jīng)過多次海底海面反射的聲線,紅色聲線表示只與海底發(fā)生一次相互作用的聲線,綠色聲線表示只與海面發(fā)生作用的聲線,藍(lán)色聲線表示只經(jīng)過海水折射,不與海底、海面相互作用的聲線.在海底附近,海水聲速與聲源位置處的聲速接近,除了極少數(shù)小掠射角的聲線可以在到達(dá)海底之前發(fā)生反轉(zhuǎn),其他大部分聲線都不滿足發(fā)生反轉(zhuǎn)的條件,因此大部分聲線會(huì)與海底發(fā)生相互作用,這個(gè)結(jié)論與圖8(a)的仿真結(jié)果相符合.從圖8(a)可以看到,大多數(shù)聲線都經(jīng)過海底海面的多次反射,剩下一小部分小掠射角聲線經(jīng)過多次海底反射,不與海面相互作用.還有極少數(shù)的小掠射角聲線可以發(fā)生反轉(zhuǎn),不與海面和海底發(fā)生相互作用.圖8(b)中的一部分聲線不與海底海面相互作用,經(jīng)過水體折射后在51 km附近形成會(huì)聚區(qū),平坦海底情況下的海底附近聲速大于聲源位置處聲速,因此小掠射角的聲線在到達(dá)海底之前就可以發(fā)生反轉(zhuǎn).大掠射角的聲線經(jīng)海底海面多次反射或者只經(jīng)過一次海底反射,在第一影區(qū)范圍內(nèi)形成較高聲能量,但是并沒有形成會(huì)聚區(qū)結(jié)構(gòu).觀察圖8中聲線軌跡,發(fā)現(xiàn)大量的聲線經(jīng)過多次海底反射到達(dá)海底后會(huì)被完全吸收,這是由于海底沉積層聲速較小,對(duì)到達(dá)海底的聲吸收較強(qiáng),除了極少數(shù)小掠射角的聲線不與海底相互作用或者相互作用次數(shù)較少,大多數(shù)經(jīng)多次海底反射后的反射損失較大,這也可以解釋圖5(a)中60 km以后聲能量較弱,會(huì)聚區(qū)結(jié)構(gòu)不明顯的實(shí)驗(yàn)現(xiàn)象.

    圖8 145 m深度聲源的聲線圖 (a)實(shí)測(cè)地形;(b)平坦地形.黑色圓點(diǎn)表示聲源深度在145 m深度,黑色虛線表示150 m互易接收深度Fig.8.Rays traces for two different bathymetry:(a)The real bathymetry;(b)the flat bottom.Black dots mark the source depth at 145 m,and black dashed line indicates the receiver depth at 150 m.

    觀察圖8(a)中150 m深度上的聲線,可見實(shí)測(cè)地形下在20和40 km附近有大量聲線到達(dá),對(duì)應(yīng)圖7(b)中的兩個(gè)海底反射會(huì)聚區(qū).在11和51 km附近也有一部分聲線到達(dá),形成了聲能量相對(duì)較高的區(qū)域,其中51 km附近是由于少量經(jīng)海水折射后到達(dá)的聲線引起的聲能量增強(qiáng),這與傳播損失仿真的結(jié)果相符合.平坦海底(圖8(b))下聲線雖經(jīng)過多次海底反射,但沒有出現(xiàn)類似圖8(a)中由于海底小山的遮擋作用導(dǎo)致的聲線在某一區(qū)域比較集中到達(dá)的現(xiàn)象,所以在10—40 km范圍內(nèi)聲場(chǎng)分布較為均勻.與圖8(a)相比,有較多的小掠射角聲線會(huì)經(jīng)過海水折射后在51 km附近形成反轉(zhuǎn)點(diǎn)會(huì)聚區(qū).這些聲線仿真結(jié)果與圖7(b)中兩種地形情況下的傳播損失結(jié)果相匹配.

    圖9給出了85,1008,和3021 m這三個(gè)互易聲源深度上的聲線圖.圖9(a)和圖9(b)兩種地形條件下的85 m聲源深度上,由于其聲速大于海底附近聲速,不滿足聲線發(fā)生反轉(zhuǎn)的條件,聲線都會(huì)與海底相互作用.圖9(a)中聲線經(jīng)不平海底的一次反射后在11,20 km附近聚集,第一海底反射會(huì)聚區(qū)20 km處的聲線經(jīng)30 km處的海底反射在40 km附近聚集形成第二個(gè)海底反射會(huì)聚區(qū),在51 km沒有海水折射的聲線到達(dá),因此在圖7(a)中的85 m接收深度上的這個(gè)位置處并沒有觀察到聲增強(qiáng)區(qū)域.圖9(b)中聲線沒有出現(xiàn)圖9(a)中在部分距離上相對(duì)集中的到達(dá)形成的會(huì)聚區(qū)域,而是相對(duì)均勻地分布在10—40 km的距離上,形成了圖7(a)中的高聲強(qiáng)區(qū)域.圖9(c)在1008 m聲道軸附近的聲源深度上,20 km距離處,左側(cè)僅有經(jīng)過一次海底海面反射的聲線到達(dá),右側(cè)除了一次海底海面反射到達(dá)的聲線,還有一部分只經(jīng)過一次海底反射到達(dá)的聲線,因此右側(cè)的聲增益高于左側(cè);在40 km處也有明顯的分裂現(xiàn)象;51 km附近并沒有聲線到達(dá),因此在圖7(c)的1008 m接收深度上的51 km距離處沒有出現(xiàn)聲增強(qiáng)區(qū)域.圖9(d)的平坦地形下,40 km以內(nèi)的聲線到達(dá)相對(duì)均勻,在45和56 km處有部分只經(jīng)過海水折射的聲線到達(dá),形成圖7(c)中兩個(gè)高聲強(qiáng)的會(huì)聚區(qū)分裂結(jié)構(gòu).

    圖9(e)和圖9(f)在海底附近的聲源深度上,20 km以內(nèi)有大量的直達(dá)聲線和一次海底反射聲線,因此聲能量較高,在30 km處,大量的兩次海底反射聲線和部分直達(dá)聲線到達(dá),傳播損失較小.但是,由于海底反射損失的聲能量較大,因此增益小于20 km以內(nèi)的區(qū)域.40 km以外的距離上,實(shí)測(cè)地形下有多次海底反射聲線到達(dá),且到達(dá)數(shù)量隨著距離的增加而增多,平坦海底下40 km以后僅有較少數(shù)的多次海底反射聲線到達(dá),且隨著距離的增加,到達(dá)聲線與海底作用的次數(shù)逐漸增加,因此傳播損失小于實(shí)測(cè)地形且差距逐漸增大.

    圖9 三個(gè)不同聲源深度對(duì)應(yīng)的聲線圖 (a)實(shí)測(cè)海底聲線,聲源深度85 m;(b)平坦海底聲線,聲源深度85 m;(c)實(shí)測(cè)海底聲線,聲源深度1008 m;(d)平坦海底聲線,聲源深度1008 m;(e)實(shí)測(cè)海底聲線,聲源深度3021 m;(f)平坦海底聲線,聲源深度3021 mFig.9.Comparison of rays traces for different bathymetries at three source depths:(a)The real bathymetry and the source depth at 85 m;(b)the flat bottom and the source depth at 85 m;(c)the real bathymetry and the source depth at 1008 m;(d)the flat bottom and the source depth at 1008 m;(e)the real bathymetry and the source depth at 3021 m;(f)the flat bottom and the source depth at 3021 m.

    圖10 四個(gè)不同距離處的本征聲線和時(shí)間到達(dá)結(jié)構(gòu) (a)和(b)表示11 km聲源的本征聲線和到達(dá)結(jié)構(gòu);(c)和(d)表示20 km聲源的本征聲線和到達(dá)結(jié)構(gòu);(e)和(f)表示42 km聲源的本征聲線和到達(dá)結(jié)構(gòu);(g)和(h)表示51 km聲源的本征聲線和到達(dá)結(jié)構(gòu)Fig.10.Comparison of the eigenrays and arrivals at four different ranges,where(a)and(b)for the source range at 11 km,(c)and(d)for the source range at 20 km,(e)and(f)for the source range at 42 km,(g)and(h)for the source range at 51 km.

    為了進(jìn)一步研究分析距離在20,40 km處兩個(gè)海底反射會(huì)聚區(qū)以及11,51 km處兩個(gè)聲強(qiáng)較高區(qū)域?qū)嶒?yàn)的聲傳播損失變化特性,計(jì)算了聲源位于這四個(gè)距離附近處的S點(diǎn)的本征聲線及對(duì)應(yīng)的時(shí)間到達(dá)結(jié)構(gòu),其中接收點(diǎn)R深度為145 m.結(jié)果如圖10所示,其中:時(shí)間到達(dá)結(jié)構(gòu)中黑色線表示聲線掠射角小于10?,藍(lán)綠色線表示掠射角在10?—20?內(nèi)的時(shí)間到達(dá)結(jié)構(gòu),其他掠射角度的到達(dá)結(jié)構(gòu)用藍(lán)色表示.

    由圖10可見,在11 km處聲線只經(jīng)過一次海底反射到達(dá),能量損失較小,可以形成聲強(qiáng)較高的區(qū)域,但是該區(qū)域內(nèi)只有大掠射角的聲線到達(dá),多途到達(dá)產(chǎn)生的時(shí)間展寬相對(duì)20 km處較大;在20 km處,由于海底小山對(duì)第一海底反射會(huì)聚區(qū)的形成有較大影響,大量聲線與海底小山的斜坡面發(fā)生相互作用,經(jīng)過一次海底反射到達(dá)接收處.第二海底反射會(huì)聚區(qū)的聲線有兩次較為集中的到達(dá),一部分小掠射角的聲線經(jīng)過一次海底反射,另一部分大掠射角的聲線到達(dá)第一會(huì)聚區(qū)后再次經(jīng)過海底反射到達(dá)42 km處,但是第二次到達(dá)的聲壓幅度相對(duì)較低.對(duì)比圖10(d)和圖10(f)發(fā)現(xiàn),與第一海底反射會(huì)聚區(qū)相比,第二海底反射會(huì)聚區(qū)內(nèi)一次海底反射到達(dá)聲線的幅值更高,因?yàn)榈诙5追瓷鋾?huì)聚區(qū)到達(dá)聲線入射海底時(shí)的掠射角小于第一會(huì)聚區(qū),海底反射損失的聲能量相對(duì)更低,而且有部分來自兩次海底反射聲線的貢獻(xiàn),因此圖7(b)的實(shí)驗(yàn)傳播損失中觀察到第二會(huì)聚區(qū)相比第一會(huì)聚區(qū)有更高的增益,而且第二會(huì)聚區(qū)多途結(jié)構(gòu)更明顯.在51 km處,聲線經(jīng)過兩次海底反射到達(dá)接收位置,大量的聲能量在與海底、海面的相互作用中衰減掉,只有極少數(shù)小掠射角的聲線只經(jīng)過水體折射到達(dá)接收位置,因此51 km聲增強(qiáng)區(qū)域內(nèi)的增益小于前面的兩次海底反射會(huì)聚區(qū).

    對(duì)比圖9中平坦海底地形和實(shí)測(cè)海底地形條件下的聲線圖,海底小山的存在對(duì)海底反射會(huì)聚區(qū)的提前形成具有重要的貢獻(xiàn).為了分析海底小山丘對(duì)于會(huì)聚區(qū)形成的具體影響,圖11給出了在第一反射會(huì)聚區(qū)內(nèi)(20 km)四個(gè)不同接收深度處(85,406,1008,和3021 m)的本征聲線、時(shí)間到達(dá)結(jié)構(gòu)及實(shí)驗(yàn)測(cè)量的時(shí)域脈沖到達(dá)結(jié)構(gòu),可見四個(gè)深度上理論計(jì)算的時(shí)間到達(dá)結(jié)構(gòu)與實(shí)驗(yàn)結(jié)果較符合.隨著接收深度的增加,第一會(huì)聚區(qū)內(nèi)到達(dá)接收點(diǎn)的聲線與海底相互作用的范圍增大,聲線到達(dá)的路徑趨于復(fù)雜.在較淺的兩個(gè)接收深度上,僅有一次海底反射的聲線到達(dá)接收點(diǎn),1008 m接收深度上的到達(dá)聲線包含一次海底反射和兩次海底反射聲線,其中一次海底反射聲線分兩次集中到達(dá),多途結(jié)構(gòu)相對(duì)復(fù)雜.在接近海底的3021 m接收深度上,除了海底反射聲線,還有直達(dá)聲線到達(dá)接收點(diǎn).直達(dá)聲線的能量損失較小,其幅值遠(yuǎn)大于其他路徑到達(dá)聲線,也大于其他接收深度,因此20 km海底附近的傳播損失小于其他深度.

    圖11 第一海底反射會(huì)聚區(qū)不同接收深度本征聲線和時(shí)間到達(dá)結(jié)構(gòu)比較 (a)85 m接收深度;(b)406 m接收深度;(c)1008 m接收深度;(d)3021 m接收深度Fig.11.Comparison of eigenrays and arrivals at different receiver depths in the first bottom re flection convergence zone:(a)The receiver depth at 85 m;(b)the receiver depth at 406 m;(c)the receiver depth at 1008 m;(d)for the receiver depth of 3021 m.

    4 總 結(jié)

    本文利用2018年4月南海聲傳播實(shí)驗(yàn)數(shù)據(jù)以及同步海洋環(huán)境資料,對(duì)深海海底反射會(huì)聚區(qū)聲傳播特性進(jìn)行了研究.基于拋物方程數(shù)值分析結(jié)合射線理論對(duì)實(shí)驗(yàn)現(xiàn)象進(jìn)行分析,很好地解釋了水平變化的海底地形環(huán)境下海底反射會(huì)聚區(qū)形成的機(jī)理.結(jié)果表明:直達(dá)聲覆蓋距離范圍內(nèi)的海底地形變化對(duì)海底反射會(huì)聚區(qū)的形成及傳播特性具有很大的影響,在20和40 km距離附近形成兩個(gè)明顯的海底反射會(huì)聚區(qū)結(jié)構(gòu),其增益達(dá)10 dB.由于海底地形隨距離逐漸變淺,破壞了水平不變深海環(huán)境下原有的會(huì)聚區(qū)結(jié)構(gòu),折射到達(dá)的聲線路徑變小,使得對(duì)應(yīng)距離下的聲場(chǎng)能量小于海底反射會(huì)聚區(qū).當(dāng)聲源深度與接收深度都處于相對(duì)較淺的相同深度時(shí),隨著一次海底反射聲線海底反射能量損失的減小以及到達(dá)聲線數(shù)量的增加,第二會(huì)聚區(qū)的增益高于第一會(huì)聚區(qū).在第一反射會(huì)聚區(qū)內(nèi),隨著接收深度的增加,聲線到達(dá)的路徑逐漸增加,聲線到達(dá)結(jié)構(gòu)趨于復(fù)雜,多途效應(yīng)更加明顯.此外,由于海底小山的反射作用,大角度聲線經(jīng)海底反射導(dǎo)致部分深度上的海底反射會(huì)聚區(qū)可提前到11 km.研究結(jié)果為聲納在深海復(fù)雜地形下的應(yīng)用及性能環(huán)境效應(yīng)分析具有重要借鑒意義.

    感謝參與2018年春季南中國海聲傳播實(shí)驗(yàn)的全體工作人員,是他們的辛勤勞動(dòng)為本文提供了可靠的實(shí)驗(yàn)數(shù)據(jù).

    猜你喜歡
    聲線聲速聲源
    虛擬聲源定位的等效源近場(chǎng)聲全息算法
    水聲中非直達(dá)聲下的聲速修正方法①
    基于聲線法的特殊體育館模型中聲場(chǎng)均勻性分析
    基于GCC-nearest時(shí)延估計(jì)的室內(nèi)聲源定位
    電子制作(2019年23期)2019-02-23 13:21:12
    糾纏的曲線
    優(yōu)雅(2017年3期)2017-03-09 17:02:52
    聲速是如何測(cè)定的
    運(yùn)用內(nèi)積相關(guān)性結(jié)合迭代相減識(shí)別兩點(diǎn)聲源
    三維溫度梯度場(chǎng)中本征聲線軌跡的求取*
    跨聲速風(fēng)洞全模顫振試驗(yàn)技術(shù)
    機(jī)翼跨聲速抖振研究進(jìn)展
    午夜福利欧美成人| 后天国语完整版免费观看| 丁香欧美五月| 亚洲国产欧美一区二区综合| 国产av一区二区精品久久| 亚洲精品国产区一区二| 国产乱人伦免费视频| 国产精品日韩av在线免费观看 | 一二三四在线观看免费中文在| 亚洲av成人av| 亚洲欧美日韩高清在线视频| 看片在线看免费视频| 欧美一级a爱片免费观看看 | 亚洲少妇的诱惑av| 无遮挡黄片免费观看| 国产在线观看jvid| 十八禁网站免费在线| 国产精品久久久久久人妻精品电影| 欧美成人一区二区免费高清观看 | 精品久久久久久久人妻蜜臀av | 成人三级黄色视频| 久久精品影院6| 香蕉丝袜av| 在线观看午夜福利视频| 国内久久婷婷六月综合欲色啪| 91大片在线观看| 精品久久久久久久人妻蜜臀av | 精品国产美女av久久久久小说| 午夜免费观看网址| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美激情综合另类| 成人18禁在线播放| 又黄又粗又硬又大视频| 久久精品亚洲熟妇少妇任你| 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av在线| 久久亚洲精品不卡| 亚洲一区中文字幕在线| 国产成人一区二区三区免费视频网站| 9色porny在线观看| 精品国产一区二区久久| 日韩欧美免费精品| 女人高潮潮喷娇喘18禁视频| 亚洲av成人一区二区三| 精品一区二区三区四区五区乱码| 淫秽高清视频在线观看| 777久久人妻少妇嫩草av网站| 如日韩欧美国产精品一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲av嫩草精品影院| 亚洲精品一区av在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品91蜜桃| 亚洲精品在线美女| 88av欧美| 色精品久久人妻99蜜桃| 久久久久精品国产欧美久久久| 人妻久久中文字幕网| 成人三级做爰电影| 神马国产精品三级电影在线观看 | 99精品在免费线老司机午夜| 精品福利观看| 国产xxxxx性猛交| 日韩欧美一区二区三区在线观看| 满18在线观看网站| 非洲黑人性xxxx精品又粗又长| 国产成人啪精品午夜网站| 美国免费a级毛片| 久久久久久免费高清国产稀缺| av中文乱码字幕在线| 久久国产精品男人的天堂亚洲| 日本五十路高清| 国产av又大| 麻豆一二三区av精品| videosex国产| 少妇粗大呻吟视频| 级片在线观看| 午夜福利影视在线免费观看| 免费在线观看亚洲国产| 在线国产一区二区在线| 成在线人永久免费视频| 国产精品乱码一区二三区的特点 | 成人亚洲精品一区在线观看| 成年版毛片免费区| 国产精品 欧美亚洲| 国产精华一区二区三区| 成年女人毛片免费观看观看9| 亚洲情色 制服丝袜| 国产亚洲精品综合一区在线观看 | 黄色视频不卡| 欧美一级a爱片免费观看看 | 少妇被粗大的猛进出69影院| 男人舔女人下体高潮全视频| 久久久久亚洲av毛片大全| 亚洲专区中文字幕在线| 身体一侧抽搐| 国产一级毛片七仙女欲春2 | 免费在线观看视频国产中文字幕亚洲| 看黄色毛片网站| 9色porny在线观看| 久久人妻av系列| 成人18禁在线播放| 亚洲男人的天堂狠狠| 久久久国产成人精品二区| 黄片大片在线免费观看| 电影成人av| 久久人妻福利社区极品人妻图片| 欧美黄色片欧美黄色片| 成人国产一区最新在线观看| 精品少妇一区二区三区视频日本电影| 欧美日韩瑟瑟在线播放| 99国产精品99久久久久| 国产精品一区二区免费欧美| 首页视频小说图片口味搜索| 精品国产一区二区久久| 免费搜索国产男女视频| 一个人免费在线观看的高清视频| 精品一品国产午夜福利视频| 国产一区在线观看成人免费| 91在线观看av| 国内久久婷婷六月综合欲色啪| 精品一区二区三区四区五区乱码| 国产av精品麻豆| 精品乱码久久久久久99久播| 亚洲伊人色综图| 天天躁夜夜躁狠狠躁躁| 欧美国产精品va在线观看不卡| 国产激情久久老熟女| 国产熟女午夜一区二区三区| 午夜福利视频1000在线观看 | 黄片大片在线免费观看| 精品国产乱码久久久久久男人| 国产精品av久久久久免费| 欧美日韩乱码在线| 中文字幕另类日韩欧美亚洲嫩草| 国产三级在线视频| 黄片播放在线免费| 老鸭窝网址在线观看| 国产欧美日韩一区二区三| 欧美不卡视频在线免费观看 | 精品久久久久久,| 国产男靠女视频免费网站| 女人被躁到高潮嗷嗷叫费观| 侵犯人妻中文字幕一二三四区| 精品第一国产精品| 757午夜福利合集在线观看| 亚洲专区中文字幕在线| 久久久久久久午夜电影| 十分钟在线观看高清视频www| 日韩大尺度精品在线看网址 | 欧美在线一区亚洲| 久久这里只有精品19| 国产激情欧美一区二区| 欧美日韩瑟瑟在线播放| 9191精品国产免费久久| 别揉我奶头~嗯~啊~动态视频| 欧美乱妇无乱码| 精品欧美一区二区三区在线| 91精品国产国语对白视频| 一级毛片女人18水好多| √禁漫天堂资源中文www| av福利片在线| 日韩 欧美 亚洲 中文字幕| 亚洲美女黄片视频| 国产精品亚洲一级av第二区| 亚洲电影在线观看av| 国产欧美日韩一区二区三区在线| 精品国产一区二区三区四区第35| 免费人成视频x8x8入口观看| xxx96com| 啦啦啦免费观看视频1| 国产高清有码在线观看视频 | 久久精品国产综合久久久| 欧美国产精品va在线观看不卡| 在线观看免费视频日本深夜| 黄色成人免费大全| 日韩免费av在线播放| 99国产极品粉嫩在线观看| 少妇被粗大的猛进出69影院| 色播亚洲综合网| 国产精品久久电影中文字幕| 国产亚洲精品久久久久久毛片| 亚洲精品粉嫩美女一区| aaaaa片日本免费| 成人国产综合亚洲| 美女高潮喷水抽搐中文字幕| 男人操女人黄网站| 精品高清国产在线一区| 亚洲精品久久成人aⅴ小说| 色精品久久人妻99蜜桃| 久久青草综合色| 嫩草影视91久久| 一个人观看的视频www高清免费观看 | 在线永久观看黄色视频| 久久亚洲真实| 国产aⅴ精品一区二区三区波| 中文亚洲av片在线观看爽| 18禁国产床啪视频网站| 国产区一区二久久| 精品少妇一区二区三区视频日本电影| 国产一区二区在线av高清观看| 色av中文字幕| 黄色成人免费大全| 美女高潮到喷水免费观看| 美女扒开内裤让男人捅视频| 999久久久精品免费观看国产| 女生性感内裤真人,穿戴方法视频| 成人手机av| 欧美久久黑人一区二区| 国产亚洲精品第一综合不卡| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 婷婷六月久久综合丁香| 露出奶头的视频| 两人在一起打扑克的视频| 中文字幕高清在线视频| 啦啦啦免费观看视频1| 亚洲欧美一区二区三区黑人| 日韩欧美在线二视频| 美女国产高潮福利片在线看| 久久婷婷人人爽人人干人人爱 | 午夜日韩欧美国产| 亚洲av五月六月丁香网| 欧美老熟妇乱子伦牲交| 满18在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉精品热| 国产亚洲欧美在线一区二区| 国产精品久久电影中文字幕| 精品一区二区三区视频在线观看免费| 久久香蕉国产精品| 欧美一级毛片孕妇| 嫁个100分男人电影在线观看| 精品国产美女av久久久久小说| 一二三四社区在线视频社区8| 午夜亚洲福利在线播放| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看 | 国产精品av久久久久免费| 午夜福利免费观看在线| 久久中文字幕人妻熟女| 别揉我奶头~嗯~啊~动态视频| 亚洲av电影不卡..在线观看| 欧美大码av| 国产在线精品亚洲第一网站| 丁香六月欧美| 黑人巨大精品欧美一区二区mp4| 日韩精品中文字幕看吧| 此物有八面人人有两片| 中文字幕另类日韩欧美亚洲嫩草| 欧美色视频一区免费| 成人国语在线视频| 日本三级黄在线观看| 亚洲天堂国产精品一区在线| 亚洲va日本ⅴa欧美va伊人久久| 国产一级毛片七仙女欲春2 | 国产又爽黄色视频| 久久午夜亚洲精品久久| 午夜成年电影在线免费观看| 亚洲色图av天堂| 美女免费视频网站| 亚洲欧美激情在线| 国产精品免费视频内射| 亚洲av五月六月丁香网| 人人妻人人澡欧美一区二区 | 19禁男女啪啪无遮挡网站| 国产精品精品国产色婷婷| 免费看美女性在线毛片视频| 欧美乱妇无乱码| www国产在线视频色| 欧美久久黑人一区二区| 亚洲人成网站在线播放欧美日韩| or卡值多少钱| 国产精品秋霞免费鲁丝片| 亚洲五月色婷婷综合| 国产精品二区激情视频| 久久中文字幕人妻熟女| 好男人在线观看高清免费视频 | 人人妻人人澡欧美一区二区 | 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 国产精品久久久人人做人人爽| 在线观看日韩欧美| 亚洲色图av天堂| 天堂√8在线中文| 男人操女人黄网站| 中文亚洲av片在线观看爽| 亚洲熟妇熟女久久| 国产激情久久老熟女| 亚洲精品在线美女| 黄片大片在线免费观看| 欧美乱码精品一区二区三区| 国产成年人精品一区二区| 亚洲精品国产一区二区精华液| 亚洲成人久久性| 91av网站免费观看| 国产精品野战在线观看| 精品人妻在线不人妻| 非洲黑人性xxxx精品又粗又长| 国产精品1区2区在线观看.| 亚洲精品美女久久av网站| 日本vs欧美在线观看视频| 亚洲精品中文字幕在线视频| 日韩av在线大香蕉| 一区二区三区国产精品乱码| 国产一区二区激情短视频| 国产精品 欧美亚洲| 亚洲一区二区三区色噜噜| 亚洲久久久国产精品| 国产高清videossex| 黑人巨大精品欧美一区二区蜜桃| 亚洲av熟女| 久久午夜亚洲精品久久| 日本五十路高清| 亚洲最大成人中文| 国产精品1区2区在线观看.| 婷婷精品国产亚洲av在线| 岛国视频午夜一区免费看| 90打野战视频偷拍视频| 国产精华一区二区三区| 久久人妻av系列| 亚洲精品在线观看二区| 啪啪无遮挡十八禁网站| 亚洲免费av在线视频| 亚洲专区字幕在线| 美女国产高潮福利片在线看| 午夜免费观看网址| 日本 av在线| av超薄肉色丝袜交足视频| 亚洲成a人片在线一区二区| 韩国av一区二区三区四区| 99re在线观看精品视频| 制服丝袜大香蕉在线| 免费不卡黄色视频| 国产成+人综合+亚洲专区| 999久久久精品免费观看国产| 国产免费av片在线观看野外av| 一级,二级,三级黄色视频| 性欧美人与动物交配| 久久精品国产综合久久久| 成人三级黄色视频| 成在线人永久免费视频| 亚洲人成伊人成综合网2020| 咕卡用的链子| 精品人妻1区二区| 亚洲久久久国产精品| 国产熟女午夜一区二区三区| 看黄色毛片网站| 欧美一区二区精品小视频在线| av有码第一页| 两性夫妻黄色片| 两个人看的免费小视频| 999精品在线视频| 亚洲精品国产精品久久久不卡| av福利片在线| 国产1区2区3区精品| 日本免费一区二区三区高清不卡 | 亚洲激情在线av| 母亲3免费完整高清在线观看| 国产精品一区二区在线不卡| 久久精品91蜜桃| 亚洲精品久久成人aⅴ小说| a级毛片在线看网站| 人妻久久中文字幕网| 成人国产综合亚洲| 村上凉子中文字幕在线| 亚洲av第一区精品v没综合| √禁漫天堂资源中文www| 麻豆成人av在线观看| 女人高潮潮喷娇喘18禁视频| 国产成人精品在线电影| 午夜a级毛片| 波多野结衣一区麻豆| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三区在线| 日本 av在线| a级毛片在线看网站| 女警被强在线播放| 麻豆av在线久日| 午夜成年电影在线免费观看| 在线观看66精品国产| 亚洲,欧美精品.| 国产精品久久视频播放| 校园春色视频在线观看| 啦啦啦韩国在线观看视频| 在线观看www视频免费| 免费看a级黄色片| 成人免费观看视频高清| 人人妻人人澡人人看| 国内精品久久久久精免费| 亚洲国产日韩欧美精品在线观看 | 日韩免费av在线播放| 国产精品免费一区二区三区在线| 亚洲精华国产精华精| 国产片内射在线| 在线观看www视频免费| 嫩草影院精品99| 久久久久久亚洲精品国产蜜桃av| 一本综合久久免费| av天堂在线播放| av在线天堂中文字幕| 亚洲欧美激情综合另类| netflix在线观看网站| 国产成人欧美| 男女床上黄色一级片免费看| 午夜精品久久久久久毛片777| 欧美大码av| 99久久久亚洲精品蜜臀av| 露出奶头的视频| 日本撒尿小便嘘嘘汇集6| 久久久久国产精品人妻aⅴ院| 乱人伦中国视频| 久久久久久久精品吃奶| 日韩欧美国产在线观看| 纯流量卡能插随身wifi吗| 少妇裸体淫交视频免费看高清 | 国产亚洲精品av在线| 久久国产精品男人的天堂亚洲| 热99re8久久精品国产| 久久香蕉精品热| 国产高清有码在线观看视频 | 黄色女人牲交| 国产亚洲欧美98| 久久婷婷成人综合色麻豆| 欧美成人一区二区免费高清观看 | 亚洲国产毛片av蜜桃av| 国产av又大| 免费在线观看视频国产中文字幕亚洲| 波多野结衣高清无吗| 国产黄a三级三级三级人| 欧美日韩瑟瑟在线播放| 色在线成人网| 国产成人精品久久二区二区91| 亚洲第一欧美日韩一区二区三区| 黑人操中国人逼视频| 国产三级在线视频| 长腿黑丝高跟| 午夜精品在线福利| 日日爽夜夜爽网站| 巨乳人妻的诱惑在线观看| 色播亚洲综合网| 人人澡人人妻人| 免费看a级黄色片| 一本久久中文字幕| 日日夜夜操网爽| 国产又爽黄色视频| 亚洲国产精品成人综合色| 成人18禁在线播放| 国产高清视频在线播放一区| 亚洲狠狠婷婷综合久久图片| 首页视频小说图片口味搜索| 在线观看日韩欧美| 久久人妻av系列| 国产av在哪里看| 成人免费观看视频高清| 97超级碰碰碰精品色视频在线观看| 亚洲av五月六月丁香网| 国产精品久久久av美女十八| а√天堂www在线а√下载| www.熟女人妻精品国产| 又紧又爽又黄一区二区| 色综合站精品国产| 免费不卡黄色视频| 久久香蕉精品热| 午夜福利视频1000在线观看 | 校园春色视频在线观看| 亚洲全国av大片| 变态另类丝袜制服| www.www免费av| 国产熟女xx| 丰满的人妻完整版| 欧美国产精品va在线观看不卡| 91麻豆精品激情在线观看国产| 成人国语在线视频| 亚洲精品中文字幕在线视频| 满18在线观看网站| 日本在线视频免费播放| 亚洲三区欧美一区| 国产午夜精品久久久久久| www国产在线视频色| 欧美日韩黄片免| 黑丝袜美女国产一区| 国产日韩一区二区三区精品不卡| 九色亚洲精品在线播放| 亚洲av成人av| 性少妇av在线| 国产成+人综合+亚洲专区| 欧美黄色片欧美黄色片| 免费在线观看亚洲国产| 成年版毛片免费区| 一区二区三区高清视频在线| 老鸭窝网址在线观看| 亚洲欧洲精品一区二区精品久久久| 97人妻天天添夜夜摸| 亚洲全国av大片| 又大又爽又粗| 精品一区二区三区av网在线观看| 亚洲自拍偷在线| 亚洲成av人片免费观看| 亚洲精品av麻豆狂野| 欧美av亚洲av综合av国产av| 久久人妻熟女aⅴ| 免费在线观看日本一区| 99re在线观看精品视频| 国产97色在线日韩免费| 美女大奶头视频| 久久久久久亚洲精品国产蜜桃av| 成人国产一区最新在线观看| 久久精品成人免费网站| 久久午夜综合久久蜜桃| 女人被躁到高潮嗷嗷叫费观| 免费观看人在逋| 日韩欧美国产在线观看| 女人被躁到高潮嗷嗷叫费观| 在线观看午夜福利视频| 日日摸夜夜添夜夜添小说| 亚洲国产精品999在线| 久久久久国内视频| 国产精品 欧美亚洲| 999久久久国产精品视频| 日本a在线网址| 天堂影院成人在线观看| 啦啦啦免费观看视频1| 国产麻豆69| 在线观看66精品国产| 香蕉丝袜av| 人人澡人人妻人| 999久久久精品免费观看国产| 国产高清有码在线观看视频 | 成人欧美大片| 欧美黄色片欧美黄色片| 窝窝影院91人妻| 国产一区二区在线av高清观看| 日日摸夜夜添夜夜添小说| 日韩视频一区二区在线观看| 美国免费a级毛片| 久久久精品欧美日韩精品| 九色国产91popny在线| 在线天堂中文资源库| 亚洲av成人一区二区三| 亚洲精品中文字幕一二三四区| 国产乱人伦免费视频| 成人亚洲精品一区在线观看| 成人永久免费在线观看视频| 高清黄色对白视频在线免费看| 变态另类成人亚洲欧美熟女 | 一级毛片精品| 国产三级黄色录像| 男男h啪啪无遮挡| 丝袜美腿诱惑在线| 欧美 亚洲 国产 日韩一| 人妻久久中文字幕网| 国产一区二区在线av高清观看| 成人亚洲精品一区在线观看| 亚洲人成电影免费在线| 国产成+人综合+亚洲专区| 19禁男女啪啪无遮挡网站| 中文字幕久久专区| 丰满的人妻完整版| 变态另类丝袜制服| 免费久久久久久久精品成人欧美视频| 9191精品国产免费久久| 日本精品一区二区三区蜜桃| 亚洲国产精品久久男人天堂| 国产精品香港三级国产av潘金莲| 两人在一起打扑克的视频| 日本五十路高清| 操出白浆在线播放| 欧美一区二区精品小视频在线| 久久香蕉激情| 国产精品综合久久久久久久免费 | 露出奶头的视频| 精品一区二区三区av网在线观看| 免费高清视频大片| 久久热在线av| 日韩欧美三级三区| 最近最新中文字幕大全电影3 | 亚洲色图av天堂| 亚洲av第一区精品v没综合| 日日爽夜夜爽网站| 岛国视频午夜一区免费看| 黄色片一级片一级黄色片| 在线观看www视频免费| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲熟妇熟女久久| 欧美黄色淫秽网站| 最好的美女福利视频网| 亚洲人成伊人成综合网2020| 好男人电影高清在线观看| 亚洲五月婷婷丁香| 国产精品亚洲一级av第二区| 欧美成人一区二区免费高清观看 | 夜夜看夜夜爽夜夜摸| 久久久久九九精品影院| 亚洲成人免费电影在线观看| 久久性视频一级片| 亚洲 欧美 日韩 在线 免费| 久久精品91无色码中文字幕| 午夜激情av网站| 狠狠狠狠99中文字幕| av片东京热男人的天堂| 成年女人毛片免费观看观看9| 午夜视频精品福利| 欧美成人一区二区免费高清观看 | 精品欧美一区二区三区在线| 国产av在哪里看| 国产成人一区二区三区免费视频网站| ponron亚洲| 中文字幕人妻熟女乱码| 国产高清视频在线播放一区| 一进一出抽搐动态| 波多野结衣高清无吗| 国产精品秋霞免费鲁丝片| 国产精品久久久久久亚洲av鲁大| 黄色片一级片一级黄色片|