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

海底聲散射特性研究進展

2017-07-07 13:19:36劉保華闞光明裴彥良楊志國于凱本于盛齊
海洋學報 2017年7期
關鍵詞:界面測量模型

劉保華,闞光明,裴彥良,楊志國,于凱本,于盛齊

(1.國家深?;毓芾碇行?,山東 青島 266061;2.青島海洋科學與技術國家實驗室 海洋地質過程與環境功能實驗室,山東 青島 266237;3.國家海洋局第一海洋研究所 海洋沉積與環境地質國家海洋局重點實驗室,山東 青島 266237)

?

海底聲散射特性研究進展

劉保華1,2,闞光明2,3*,裴彥良2,3,楊志國1,于凱本1,2,于盛齊1

(1.國家深?;毓芾碇行?,山東 青島 266061;2.青島海洋科學與技術國家實驗室 海洋地質過程與環境功能實驗室,山東 青島 266237;3.國家海洋局第一海洋研究所 海洋沉積與環境地質國家海洋局重點實驗室,山東 青島 266237)

海底是水下聲場的重要邊界,其聲散射特性對水下聲場空間結構及分布規律具有至關重要的影響。對目前國際上海底聲散射特性研究方面的進展進行了系統的分析和總結,從海底聲散射測量技術、海底聲散射特性及機理、海底聲散射預測模型3個方面進行了論述,并提出了未來研究的方向、研究重點與難點。該工作對于充分了解和認識海底聲散射研究的目前現狀和未來發展趨勢具有很好的借鑒和指導意義。

海底聲散射特性;散射模型;海底混響;水下聲場

1 引言

水下聲場空間結構及分布規律研究在水下目標探測、水下通訊和導航、水文測量和海底地形地貌勘測等軍事、海洋科學及生產領域具有重要的應用價值。海底是水下聲場的一個重要邊界,包括海底聲散射特性在內的海底聲學特性對水下聲場的空間結構和分布規律具有至關重要的影響。海底聲散射是海洋混響的重要來源,而海洋混響則是主動聲吶探測的主要干擾源。因此,研究海底聲散射有利于包括海底混響在內的海洋混響的準確預報,對于提高水下聲吶性能和目標探測的精度具有重要意義。另外,海底聲散射研究還與淺海聲場預報、匹配場定位、潛艇著底隱蔽地點選擇等應用具有密切的關系。

海底聲散射特性與海底粗糙度、沉積物非均勻性以及聲波掠射角密切相關,對于不同的測量頻率,其相關性變化規律也不相同。而且,不同的海底環境條件以及不同的聲波頻率,其聲散射機理也不相同。因此,深入認識海底聲散射特性、影響機理以及預測模型,是進行水下聲場空間結構研究、水下聲場準確預報和水下目標精確探測的重要環節。鑒于海底聲散射特性的重要性,早在20世紀50年代,國外科學家就開始了海底聲散射特性的測量和研究[1]。之后,不同研究者分別從測量技術、特性及機理分析、預測模型等方面開展了大量的研究工作。與國外相比,目前國內無論是在海底聲散射測量技術還是在散射機理和模型等方面開展的研究還非常少,存在很大的差距。

本文對目前國際上海底聲散射特性測量技術及機理和預測模型方面的研究成果進行了總結,闡述了國際上在海底聲散射測量技術研發、海底聲散射特性測量及機理分析、海底聲散射預測模型研究等方面的發展現狀,并提出了目前研究工作中尚未解決的問題及未來研究的方向、研究重點與難點,以期能夠對我國將來海底聲散射特性的研究工作提供借鑒和指導。

2 海底聲散射特性測量技術研究進展

美國海底聲散射測量開始于20世紀50年代。Urick[1]將一收發合置的圓柱活塞換能器固定在木桿上,借助于駁船將換能器放置在靠近海底處,進行了最早期的海底反向散射測量。測量的頻率范圍為10~60 kHz,通過水平和垂直方向旋轉木桿,來獲得不同方位角和掠射角的測量數據。Urick和Saling[2]又采用炸藥聲源對水深為4 400 m的海底進行了反向散射測量,聲源的頻帶約為500~8 000 kHz,炸藥的爆炸深度和水聽器沉放深度均為15 m。Wong和Chestermax[3]采用一個磁致伸縮的換能器作為聲源和接收水聽器,在香港的近岸海域進行了海底的反向散射測量,測量的頻率為48 kHz,發射聲波的脈沖寬度為0.4 ms和2.8 ms兩種。總之,在20世紀50-60年代,海底聲散射測量剛剛起步,實驗所使用的均是非常簡單的裝置,還未開展專業的海底聲散射測量技術研究。

20世紀70年代末,美國華盛頓大學應用物理實驗室的Barry等[4—5]研制出一種拖曳式海底反向聲散射測量裝置(圖1)。通過一個類似于球狀萬向節(ball-in-socket)的調節機構,將發射和接收合置的聲學換能器安裝在拖曳平臺上,調節機構可以方便實現換能器掠射角的調整,平臺上還安裝由測深和姿態傳感器。Jackson等[6]采用該設備對粉砂、砂質和礫石海底度進行了反向聲散射測量,通過安裝3個不同頻帶的收發合置平面換能器獲得了20~85 kHz的頻帶覆蓋范圍,換能器垂直方向發射和接收波束角為20°~40°,水平波束角為10°~20°,聲源級為205 dB re.μPa@1 m,實驗中平臺的拖曳速度為2~5 kn。Stanic等[7]研制出一種適用于淺水的坐底式海底聲散射系統(圖2),整套系統由兩個坐底式的塔式支架組成,兩個支架安裝由16個水聽器組成的T型接收陣,其中一個還安裝有高頻參量陣聲源,聲源的差頻頻率為20~180 kHz,差頻聲源級為187 dB re.μPa@20 kHz和214 dB re.μPa@180 kHz,該系統可進行高頻的海底反向聲散射和前向聲散射測量。除上述兩套海底聲散射測量設備外,Boehme等[8]將發射和接收換能器安裝在一個高4 m的三角架上,將該裝置放置在海底進行了頻率范圍為30~95 kHz的海底反向聲聲散射測量。從上述的分析可以看出,在20世紀70-80年代,海底聲散射測量技術得到快速發展,研發出了多臺套的專門用于海底聲散射測量的專業設備,測量精度得到很大提高。本階段的測量頻率主要集中在20 kHz以上的高頻。

圖1 Barry等研發的拖曳式海底聲散射測量裝置(改自文獻[5-6])Fig.1 Towed sonar system for bottom scattering measurement(modified from reference[5-6])

圖2 Stanic等研發的坐底式海底聲散射測量系統(改自文獻[7])Fig.2 Bottom-mounted bottom scattering measurement system(modified from reference[7])

20世紀90年代中,Greaves和Stephen[9]采用由10個低頻彎張換能器組成的垂直線陣聲源和由128個水聽器組成的水平接收陣在大西洋中脊進行了海底聲散射測量。垂直線陣聲源的彎張換能器的排放間距為2.29 m,通過各換能器的時延相控發射,聲源可以形成俯角為9°的相控波束,在實驗中,聲源進行LMF掃頻發射,掃頻寬度為200~255 Hz,信號長度為5 s。接收陣列的基元間距為2.5 m,通過波束形成技術,產生了126個具有不同方向的接收波束,但0°~30°和150°~180°兩個波束方向的數據因波束太寬且受接收陣姿態變化和船舶噪音的影響大而無法使用[9]。20世紀90年代末和21世紀初,美國海軍研究辦公室聯合華盛頓大學、Scripps海洋研究所、意大利NATOSACLANT海底科學研究中心等科研機構開展了兩個綜合的海底聲學實驗,分別為:SAX99(Sediment Acoustic Experiment-1999)和SAX04(Sediment Acoustic Experiment-2004)[10—12]。在SAX99實驗中采用了BAMS(Benthic Acoustic Measurement System)、STMS(Sediment Transmission Measurement System)、XBAMS(Accelerated Benthic Acoustic Measurement System)3種系統進行海底聲散射測量[10-11]。BAM系統的工作頻率為40 kHz和300 kHz,XBAMS的工作頻率為300 kHz,二者的換能器均可以按一定的波束角在水平方向上步進旋轉,以獲得不同方位角的海底聲散射強度(圖3)。STMS為一個可用于海底沉積物聲衰減測量、海底聲透射和海底聲散射測量的綜合測量系統,其海底聲散射測量的工作頻率為20~150 kHz,換能器不能夠自動旋轉,需要潛水員通過移動設備來進行不同區域的海底散射測量。SAX04項目使用的設備類似于STMS的固定式換能器,但較SAX99具有如下兩點改進:(1)測量頻帶拓寬為20~500 kHz;(2)將聲學測量設備安裝在一個長約28 m的鋪設在海底的導軌上,聲學測量設備在導軌上按照指定步長自動移動,以減少潛水員人工移動設備對測量產生的影響[11]??v觀20世紀90年代的海底聲散射測量技術發展,主要有如下兩個特點:(1)新技術不斷被應用到海底聲散射測量,如:低頻彎張換能器技術,時延相控發射技術,多基元接收波束形成技術,步進自動控制技術等;(2)同步開展海底粗糙度、沉積物非均勻性等環境參數測量以及相關技術研發,以便能夠建立精細的海底聲散射預測模型。

進入21世紀,海底聲散射測量及相關技術研發在很多國家得到廣泛重視。2000年,意大利SACLANT海底研究中心的Holland等[12]采用如圖4所示垂直接收陣和組合換能器聲源對400~4 000 Hz頻帶范圍內的海底反向聲散射進行了測量。實驗采用的ITC4001換能器聲源由3個換能器組成,通過不同組合間隔產生頻率分別為1 200 Hz、1 800 Hz、3 600 Hz的指向性聲波,彎張換能器聲源由兩個間距為1.25 m的低頻彎張換能器組合產生頻率為600 Hz的指向性聲波。接收采用間距為18 cm的32基元垂直接收陣。加拿大國防研究發展中心的Hines等[13]研發了一種用于淺海中頻小掠射角海底聲散射測量的設備,對頻率為4 kHz和8 kHz、掠射角為3°~15°的砂質海底反向聲散射特性進行了測量。系統主要由參量陣聲源、超指向性線列陣、聲強立體接收陣、安裝平臺等部分組成(圖5)。參量陣聲源由9個基元組成,4 kHz和8 kHz的發射響應分別為185 dB和192 dB(測量距離為5 m),水平和垂直波束寬度為4°~7°。超指向性接收陣由6個小型全向性水聽器組成,基元間距為16 cm。安裝參量陣和超指向性線列陣的支架可以360°旋轉,水平和垂直轉角的測量精度為±1°。2010年,韓國漢陽大學的La和Choi采用單個的全向性聲源和全向性水聽器在韓國南部的近岸淺水海域開展了頻率為8 kHz的海底聲散射測量[14]。Manik等[15]使用定量回聲測深儀進行了海底反向散射測量。近幾年的海底聲散射測量的顯著特點是:研究人員將研究重點轉向了10 kHz以下的中低頻海底聲散射的測量和研究,主要原因在于中低頻聲吶在水聲通訊、水下探測等方面的廣泛應用。

圖3 SAX99實驗海底聲散射測量觀測系統示意圖(改自文獻[11])Fig.3 Geometry of bottom scattering in SAX99(modified from reference[11])

圖4 中低頻(400~4 000 Hz)海底聲散射測量裝置(改自文獻[12])Fig.4 The device used for bottom scattering measurement at low-mid frequency(400-4 000 Hz) (modified from ref-erence[12])

圖5 淺海小掠射角海底反向散射測量實驗(改自文獻[13])Fig.5 Geometry of bottom backscattering at low grazing angle in shallow water(modified from reference[13])

在國內,金國亮等[16]在放置于海底的支架上安裝10個換能器,輪流以其中一個作為發射,其余作為接收,進行了頻率為10 kHz的海底聲散射測量[16]。近幾年,宋磊[17]和薛婷[18]研究了采用具有指向性平面換能器和T型接收陣進行海底散射系數的實驗測量方法,并在淺海進行了初步的實驗測量。曹正良等[19]采用T矩陣方法對平面海底界面上球體目標的聲散射建模進行了研究??傮w來說,國內目前還沒有用于海底聲散射測量的專業設備,研究主要集中于理論分析和數值仿真。

3 海底聲散射特性及機理研究進展

20世紀50-60年代,海底聲散射特性研究主要是分析海底反向聲散射與掠射角、聲波頻率、發射脈沖長度、海底底質類型等參數的關系[1—3,20]。主要結論如下:(1)海底聲散射總體上隨著掠射角的增大而增加,但對于不同的海底類型和掠射角范圍,二者的函數關系則不同;(2)絕大部分研究人員認為海底聲散射強度不存在明顯的頻率依賴性,或僅存在很弱的頻率依賴性;(3)海底聲散射強度與發射脈沖長度不存在明顯的相關性;(4)雖然散射強度與海底沉積物顆粒粒徑不存在明顯變化規律,但對于不同類型的海底,聲散射強度還是存在一些普遍的趨勢,即:砂質和巖石等硬質海底聲散射強度一般大于黏土和粉砂等軟質海底的聲散射強度。在聲散射機理方面,研究人員得出一些初步的認識,普遍認為海底聲散射主要是由海底的粗糙度或微起伏所引起的。但對于海底沉積物顆粒對聲散射的貢獻,不同研究者給出了不同的結論。Urick[1]認為海底聲散射主要由海底粗糙性(即不規則性)引起,而不是沉積物顆粒對聲波的散射。Wong和Chestermax[3]則指出,對于48 kHz的聲波,砂質或更大粒徑的顆粒是小掠射角時海底聲散射的主要散射體,在較大掠射角和近垂直入射時,海底粗糙散射是主要機制。McKinney和Anderson[20]指出沉積物的顆粒屬性也是引起海底聲散射的一個重要方面,但顆粒和粗糙性并不是相互獨立的兩個方面,沉積物顆粒堆積在一起形成與聲波波長尺寸相當的散射體,這是海底散射與顆粒具有一定相關性的原因;另一方面,這些顆粒堆積體也正是形成海底微起伏(即粗糙性)或沉積層結構的重要因素。除此之外,Urick和Saling[2]基于中低頻(500~8 000 Hz)聲散射數據指出海底沉積物中沉積層反射可能是引起海底聲散射增加的一個因素。

20世紀70-80年代,海底聲散射測量主要集中在淺海(水深小于50 m)和高頻(20~180 kHz)反向聲散射測量,海底底質包括淤泥、粉砂、細砂、含貝殼層、礫石和巖石等多種類型。此時期主要的研究進展如下:(1)海底聲散射強度隨掠射角的增大而增強,二者關系可以采用Lambert法則公式較好地擬合,即:BS=10lgμ+10lg(sin2θ),θ為掠射角,10lgμ為垂直入射時的聲散射[21—24];(2)聲散射存在微弱的頻率依賴性,但不同研究者依據不同的底質類型和測量頻率分別給出了不同的變化關系,總體來說,聲散射與聲波頻率的關系很難用一個簡單的函數來表述,其與海底底質類型、海底粗糙度、測量頻率等多個參數有關[6,21—24];(3)開展了聲散射強度與方位角關系的研究,Stanic等[22—23]在Panama城附近海域的砂質海底和Jacksonville附近海域的含有粗貝殼層海底的測量數據均未顯示出明顯的方位角依賴性,而Boehme等[21]通過在水平方向上以不同的方位角掃描具有良好分界的細砂和粗砂海底,發現了在分界線處海底散射強度的明顯變化;(4)在開展海底聲散射測量的同時,采用側掃聲吶、水下攝像、立體照相、高分辨率測深、淺地層剖面、沉積物巖心分析等技術對海底粗糙性和非均勻性等底質環境進行了測量,為深入研究散射機理和模型奠定了基礎[21];(5)雖然海底粗糙性是引起海底聲散射的主要機制,但多個海區的數據表明海底散射強度與均方根高度沒有明顯的相關關系,這說明對于海底聲散射來說海底均方根粗糙度不是海底粗糙性的有效表述參數,研究者開始采用海底粗糙度譜來表征海底粗糙度[23]。在海底聲散射機理研究方面,在目前研究的頻率范圍內(20~180 kHz),顆粒散射不是海底聲散射的主要機制,這一點得到普遍共識。Jackson等[6]綜合多種數據指出,對于顆粒粒徑大小相似的海底散射強度相差10~15 dB。研究人員普遍認為,海底粗糙散射和體積散射是海底散射的主要機制,但對于兩種散射機制對海底聲散射的貢獻程度及適用條件,不同研究研究人員給出了不同的結論:Bunchuk和Zhitkovskii[24]認為在淺水區主要是體積非均勻性而不是界面粗糙性來主導海底聲散射;而Jackson等[25]則認為,對于淤泥和粉砂等軟質海底,在除了很小和很大掠射角之外的中等掠射角范圍內,體積散射占主導作用,而對于粗砂等硬質海底,在很寬的掠射角范圍內粗糙散射均占主導地位。Jackson等開發了包括粗糙散射和體積散射在內的復合粗糙度散射模型,但該模型并未給出體積散射強度的計算公式,而是采用一個自由變量來代替,需通過數據擬合來確定[25]。如何能夠更好地揭示粗糙散射和體積散射兩種機制,則需要更完善的模型以及更多的聲學和底質參數測量數據。

20世紀90年代,海底聲散射特性測量和研究具有如下進展和特點:(1)除傳統的高頻(20~300 kHz)測量外,研究人員對1 000 Hz以下的低頻海底聲散射特性進行了測量和散射機理研究:Tang等[26]采用頻率為250~650 Hz的深拖地震技術測量了大西洋中脊附近沉積物盆地的海底聲散射特性;Greaves和Stephen[9]采用200~255 Hz的聲波對大西洋中脊地殼特征進行了測量,試圖建立聲散射強度與地殼傾角的相關關系以便能夠采用聲散射數據來識別具有不同傾角的地殼類型,但未能給出滿意的結果。(2)測量和研究了海底淺表層沉積物中存在的氣泡所引起的聲散射特性及模型:Tang等[27]對比了砂質海底與含氣泡層的粉砂海底的聲散射特性,指出砂質海底的聲散射主要是由海底界面散射引起,而粉砂海底聲散射主要由淺表層沉積物中的氣泡對聲波的散射引起;Chu等[28]進一步建立了一個雙層的模型來模擬由海底氣泡引起的雙基地散射強度,模型假設氣泡為扁圓球形,且不同氣泡散射不存在相干性,均為一次散射。(3)開展了高頻的收發分置的雙基地海底聲散射的測量和研究:Stanic等[29]采用如圖2所示裝置測量了頻率為20~180 kHz的海底雙基地散射特性;Williams和Jackson[30]分別測量了頻率為40 kHz的粗砂海底和含氣的粉砂海底前向散射強度,并將測量數據與基于反向散射模型推導出的前向散射模型的預測結果進行了對比。(4)SAX99和SAX04綜合聲學實驗開始實施,實驗的最大特點是同步開展了精細的海底環境特性測量和精確的高頻海底聲散射同步測量,對海底粗糙度和沉積物非均勻性特征進行了詳細測量和表征,為高頻海底聲散射特征及聲散射模型綜合研究提供了大量的基礎數據[10—11,31—32]。在聲散射機理研究方面,除海底粗糙界面引起的聲散射之外,研究者開始關注海底之下沉積物特性引起的聲散射,并給出了多種海底之下聲散射的機制[25,27—28,33—35]。Jackson等建立的復合粗糙度散射模型雖然指出沉積物體積散射的存在,并沒有分析引起沉積物體積散射的具體機制[25];Hines認為海底沉積物孔隙度的波動是引起體積散射的重要機制[33];Lyons等[34]將海底之下聲散射歸結為兩種機制:一是與Hine提出的孔隙度波動相類似的沉積物的非均勻性,另一種是海底沉積物層引起的粗糙界面聲散射;Tang等[27]和Chu等[28]認為,對于淺表層沉積物含氣泡的粉砂質海底,氣泡對聲波的散射是海底散射的主要機制;Jackson和Ivakin[35]認為除了密度和縱波速度(聲速)非均勻性之外,剪切波速度非均勻性也是引起體積散射的主要因素。總之,到目前為止,研究人員普遍認識到海底散射主要由海底界面散射和來自海底沉積物的聲散射組成。海底的粗糙性是引起界面散射的主要機制,可采用海底粗糙度譜來表征海底的粗糙性。海底沉積物的聲散射來源于多種機制,包括由海底沉積物非均勻性(密度、孔隙度、聲速、剪切波速度等非均勻性)引起的體積散射、海底氣泡引起的聲散射、海底沉積物層(或基底)粗糙性引起的聲散射、非連續沉積層引起的體積散射等。因此,在進行海底聲散射特性研究時,應針對不同的聲波頻率和特定的海底環境條件,認真分析引起聲散射的機制,以便更好地對聲散射特性進行建模。

21世紀以來,中低頻的海底聲散射特性測量和研究獲得廣泛關注。Holland等[12]采用圖4所示的方法分別在泥質和巖漿巖海底測量了頻率為400~4 000 Hz的海底反向散射強度,掠射角為10°~40°。測量結果表明,在該頻段內巖漿巖海底的散射強度均大于泥質海底。機理分析認為,泥質海底聲散射主要來自于海底之下沉積層的體積散射,而海底界面散射非常?。粠r漿巖海底聲散射在低于臨界掠射角時為海底界面散射,大于臨界掠射角時主要為體積散射。Soukup和Gragg[36]采用由全向性換能器組成的線性聲源陣和9基元垂直接收陣對石灰巖海底進行了頻率為2~3.5 kHz的海底反向聲散射測量,并分析了此頻率段內海底反向聲散射強度與掠射角的關系。Hines等[13]采用如圖5所示裝置在兩個砂質海底站位測量了頻率為4 kHz和8 kHz的海底反向散射,掠射角為3°~15°。總體來說,目前中低頻海底聲散射特性測量逐步開始,但其散射機理及預測模型還未開展系統性的研究。

4 海底聲散射模型研究進展

4.1 Lambert模型

Lambert模型(或稱為Lambert定律)能夠很好地描述粗糙面上反向散射強度以及在相對光滑界面上(組成界面的小散射體的幾何尺度遠小于入射聲波波長)較小掠射角情況下的反向散射強度與掠射角的關系[37]。但是,Lambert定律是半經驗公式,缺乏具體的物理含義,而且,Lambert定律僅僅描述了散射強度與掠射角的關系,與海底粗糙度參數無關。通過該定律只能獲得入射聲線、散射聲線及散射面法線位于同一平面時的散射強度。Ellis[38],Caruthers和Novarini[39]對Lambert定律進行改進,給出了一種與Lambert定律類似的三維散射強度計算模型。

4.2 基爾霍夫(Kirchhoff)近似模型

基爾霍夫近似模型適合粗糙界面起伏變化比較緩慢的情況(Rc?λ/πsin3θ,Rc為界面的曲率半徑,λ和θ為入射波局部掠射角和波長),而對海底起伏高度和界面斜率沒有要求。Echart最早將基爾霍夫近似用于海面的聲散射研究,Jackson等[25]將基爾霍夫近似用于鏡面反射方向附近的頻率為10~100 kHz的海底聲散射的計算,并將近垂直的方向散射測量數據與基爾霍夫模型預測進行了對比。Moustier[40]將基爾霍夫近似應用于多波束數據的近垂直入射時聲散射與海底底質屬性關系的研究。Dacol[41]曾將基爾霍夫近似用于具有隨機粗糙特性的彈性海底的聲散射研究。Williams和Jacson[30]將基爾霍夫近似用于收發分置的雙基地散射模型的研究。對于高頻聲散射,Jackson和Richardson[42]等給出了高頻極限條件下的基爾霍夫近似的散射截面表達式?;鶢柣舴蚪浦形纯紤]影區效應和多次散射,目前基爾霍夫近似主要用于傾角變化緩慢的粗糙海底近垂直方向附近的海底聲散射場的計算。

4.3 小粗糙度微擾近似模型

微擾近似理論適用于粗糙界面起伏很小(小于入射波波波長)且界面斜率足夠小的情況,因此稱為小粗糙度微擾近似模型。微擾理論認為當粗糙界面起伏比較小時,它對聲場的影響比較小,可以將聲場看成由光滑界面產生的場和一個擾動散射場疊加組成。小粗糙度微擾近似是在海底聲散射場計算中得到廣泛應用的模型,根據導出散射截面表達式時所基于的波動方程類型的不同,又可以分為小粗糙度微擾近似流體模型、小粗糙度微擾近似彈性模型以及小粗糙度微擾近似多孔彈性模型。Kuo[43]給出了無損耗流體沉積物情況下的小粗糙度微擾流體近似模型表達式,Jackson和Briggs[44]將表達式推廣到有限衰減的情況,并將模型預測與砂和粉砂等不同類型的海底聲散射測量數據進行了對比,但該模型假設沉積物各向同性,且不存在任何形式的分層或梯度。后來,Moe和Jackson[45]將微擾流體模型應用于具有聲速梯度的海底,給出了此種情況下的海底聲散射強度和散射截面的一般表達式。Kuperman等研究人員將微擾流體模型擴展至分層海底的情況,并將剪切效應納入到海底聲散射的研究,建立了各向同性小粗糙度微擾近似彈性模型,用于彈性海底聲散射的研究[46—48]。實際海底是由流體和固體骨架組成的兩相多孔彈性介質,Williams等[49]將Biot模型與小粗糙度微擾近似相結合來研究海底聲散射問題,建立了小粗糙度微擾多孔彈性模型。Jackson和Richardson[42]對該模型運算進行了改進,采用矩陣方法來表示散射截面,從而降低代數運算的復雜性。

4.4 小斜率近似模型

小斜率近似模型最早為研究海面聲散射而建立的,后來被廣泛應用于海底聲散射的研究,其表達式是關于界面斜率的一系列展開式[50—52]。普遍認為,小斜率近似比小粗糙度微擾近似和基爾霍夫近似更為精確,且一種近似方法幾乎可以涵蓋幾乎所有的掠射角范圍(基爾霍夫近似適合鏡像反射方向附近的散射,小粗糙度微擾近似適合鏡像反射附近之外的散射)。Broschat和Thorsos[53]詳細推導了適合Dirichlet邊界條件的粗糙界面小斜率近似模型的散射截面表達式,并詳細研究了小斜率近似用于海底聲散射計算時的精度,小斜率近似的精度與界面均方根傾角、界面相關長度和入射角有關。目前,小斜率近似模型是應用最廣泛的海底聲散射模型。Gragg等[54]采用小斜率近似模型研究包含均勻散射體的隨機粗糙海底的聲散射,推導出此種條件下的雙基地散射公式。Soukup等[55],Jackson[56]將小斜率近似應用于包含剪切效應的彈性海底聲散射的研究,Jackson將其推廣到層狀海底聲散射研究。

4.5 復合粗糙度近似模型

復合粗糙度近似模型的核心思想是將海底粗糙度劃分為大尺度和小尺度兩部分,充分利用基爾霍夫近似和小粗糙度微擾近似在適用條件上的互補性,對于大尺度粗糙度采用基爾霍夫近似計算海底散射,對于小尺度粗糙度采用小粗糙度微擾近似。早期主要用于海面聲散射的研究,后來被Jackson應用于海底聲散射研究,其將垂直入射附近的大尺度粗糙度采用基爾霍夫近似計算海底散射,小粗糙度采用微擾近似計算海底散射[23]。復合粗糙度近似模型的問題在于用于區分大小尺度粗糙的截止波數的選擇很難十分清楚地確定,一般通過多次的試算和比較來確定。Thorsos[57]指出,對于小掠射角的聲散射計算,復合粗糙度近似模型的精度偏低。研究人員在早期的海底聲散射模型研究中采用復合粗糙度模型,目前復合粗糙度模型逐漸被小斜率近似模型代替。

4.6 其他模型

針對具有分層的海底中低頻聲散射,Jackson等[58]建立了GABIM模型。GABIM模型將海底視為多層的流體,綜合運用基爾霍夫近似、一階微擾近似和經驗公式來計算海底界面粗糙散射,采用一階微擾近似和經驗公式來計算體積散射,該模型對100~10 000 Hz的中頻聲波海底散射強度具有較好的預測效果。國內研究人員在海底聲散射和海底混響計算方法和模型方面也開展了卓有成效的研究,張仁和和金國亮[59],李風華等[60]采用簡正波理論對淺海環境下的混響特性進行了研究。彭朝暉等[61]結合Ivankin等提出的射線管積分法和Hines的復射線法,推導出了一種計算隨機非均勻海底和粗糙界面引起的平面內海底散射模型,此模型中考慮了折射波和衰減波的貢獻。

5 發展趨勢及研究難點

5.1 海底中低頻聲散射測量及模型研究

從上述分析來看,國內外對于海底聲散射特性研究主要集中在大于10 kHz的高頻段,而對頻率小于10 kHz的中低頻聲散射研究相對較少,目前所開展的相關測量也只是局限在中頻段的幾個離散頻率點,缺少系統性。中低頻聲散射研究可以彌補目前海底聲散射研究在頻段上的局限性。而且,中低頻聲波在海洋聲場測量與預報、海底埋藏物聲學探測、水聲通訊以及海洋地球物理探測等領域具有重要的應用價值。海底中低頻聲散射的系統性測量以及相關模型研究已成為海底聲散射研究發展的趨勢和研究熱點。相對于高頻聲散射而言,海底中低頻聲散射無論是聲散射特性測量,還是散射機理分析和預測模型研究,均更為復雜。測量的復雜性在于常規中頻換能器聲源的指向性差,易受海面聲散射的干擾,影響測量精度;散射機理和模型的復雜性在于其可能同時受海底粗糙界面散射和沉積物內非均勻體散射的共同作用,從而導致其預測模型更為復雜。研發聲學參量陣等中低頻高指向性聲源和多基元接收陣列是解決中低頻聲散射精確測量的途徑之一。在機理分析方面,與中低頻聲散射相匹配的海底粗糙度和海底非均勻等參數的測量和表征是一個很大的挑戰,也是今后需要重點解決的難點之一。而且,海底中低頻聲散射受海底粗糙面散射、海底體積散射、海底沉積層散射等多種機制共同作用,徹底弄清各種機理及其對聲散射的影響和貢獻,也是今后海底聲散射機理研究的難點和需要重點關注的研究課題。在中低頻模型研究方面,需要重點解決經過高頻聲散射數據驗證了的高頻聲散射模型在中低頻段的適用性問題,通過數據和模型對比,研究建立適合中低頻海底聲散射特性預測的地聲模型。

5.2 小掠射角和各向異性粗糙海底等特殊條件下的聲散射特性研究

小掠射角(小于10°)聲散射是引起遠距離混響的重要因素,在遠程聲探測方面具有重要的應用價值,是在今后的研究中需要重點關注的研究課題。淺海小掠射角聲散射測量的難點在于海面散射和多途聲傳播的影響。在散射機理研究方面,當聲波以小于臨界掠射角的小掠射角入射到海底時,有時在近海底沉積物中產生漸逝耗散波,這增加了小掠射角海底聲散射機理的復雜性,也是小掠射角聲散射研究的難點之一。在模型研究方面,目前常用的散射模型在小掠射角條件下的預測能力還未進行足夠的測試,其適用性還有待研究。另外,目前海底聲散射測量和研究往往將海底看成是統計意義上的各向同性粗糙海底,主要關注海底聲散射隨掠射角的變化。而當海底存在明顯方向性的地貌特征時(比如方向性的海底沙波或方向性明顯的海底沖刷溝槽等),海底粗糙度表現出明顯的各向異性特征。由各向異性粗糙海底引起的聲散射強度,不僅與掠射角有關,還與方位角有關。各向異性粗糙海底聲散射測量和模型研究主要難點在于海底三維粗糙特性的測量和表征以及三維聲散射場的精確測量,這將會成為未來發展的方向。

5.3 多層海底聲散射模型研究

目前的海底散射模型一般將海底看作為未分層的半空間介質,而實際的海底沉積物中往往存在分層結構。在與海底相互作用過程中,聲波往往透射進入海底界面之下(特別是對于中低頻聲波),被沉積物中的非均勻體和粗糙沉積層界面再次散射。采用更接近實際海底條件的多層聲散射模型進行海底聲散射特性預測將成為未來發展的趨勢。雖然Jackson等在多層聲散射模型方面開展了一些研究,但建立的GABIM模型將多層海底簡化成流體模型,且模型的預測結果未進行實際數據的驗證[58]。在多層海底散射模型方面今后還需要開展大量的針對實際海底情況的多層散射模型普適化和模型驗證工作。多層海底聲散射模型研究的難點一方面在于如何甄別海底界面粗糙散射、沉積物中體積散射、海底沉積層界面粗糙散射、沉積物聲傳播衰減等多種機制對散射總場的影響,另一方面是如何測量和表征海底以下多個沉積層界面粗糙度特性和界面之間沉積物的體積非均勻性。

5.4 雙基地和前向聲散射特性研究

目前國內外研究主要集中在收發合置(單基地)海底散射強度的測量和建模。相對于單基地聲吶,雙(多)基地聲吶具有探測范圍大、隱蔽性好、抗干擾能力強等優點。隨著雙(多)基地聲納的廣泛應用,研究不同入射和散射掠射角以及不同入射和散射方位角時的海底雙基地前向散射成為一種十分迫切的需要。研究人員在高頻海底前向聲散射測量方面做了一些嘗試試驗,但總體研究水平較低,特別是目前還沒有開展中低頻海底前向聲散射測量。對于中低頻前向聲散射測量,研發滿足測量要求的指向性聲源和指向性接收換能器或接收陣列是需要解決的關鍵技術。對于淺海中低頻前向海底聲散射測量,如何避免海面聲散射影響和聲波的多途傳播,是未來研究的難點。在散射模型研究方面,目前的海底聲散射模型是否還適用于海底前向聲散射特性的計算和預測,還需要開展大量的研究工作。

[1] Urick R J. The Backscattering of sound from a harbor bottom[J]. The Journal of the Acoustical Society of America, 1954, 26(2): 231-235.

[2] Urick R J, Saling D S. Backscattering of explosive sound from the deep-sea bed[J]. The Journal of the Acoustical Society of America, 1962, 34(6): 1721-1724.

[3] Wong H K, Chestermax W D. Bottom backscattering near grazing incidence in shallow water[J]. The Journal of the Acoustical Society of America, 1968, 44(6): 1713-1718.

[4] Barry W, Jackson D, Schultz J. A flexible towed sonar for ocean acoustic measurements[C]//Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, Tulsa, OK, USA, 1978:152-154.

[5] Barry W A, Jackson D R. Split-beam towed sonar for ocean acoustic measurement[C]//OCEANS’80 Conference Proceedings, Seattle, WA, USA, 1980:267-271.

[6] Jackson D R, Baird A M, Crisp J J, et al. High-frequency bottom backscatter measurement in shallow water[J]. The Journal of the Acoustical Society of America, 1986, 80(4): 1188-1199.

[7] Stanic S, Eckstein B E, Williams R L, et al. A high-frequency shallow water acoustic measurement system[J]. IEEE Journal of Oceanic Engineering, 1988, 13(3): 155-161.

[8] Boehme H, Chotiros N P, Rolleigh L D, et al. Acoustic backscattering at low grazing angles from the ocean bottom. Part I. Bottom backscattering strength[J]. The Journal of the Acoustical Society of America, 1985, 77(3): 962-974.

[9] Greaves R J, Stephen R A. Seafloor acoustic backscattering from different geological provinces in the Atlantic Natural Laboratory[J]. The Journal of the Acoustical Society of America, 1997, 111(1): 193-208.

[10] Thorsos E I, William K L, Chotiros N P, et al. An overview of SAX99: acoustic measurements[J]. IEEE Journal of Oceanic Engineering, 2001, 26(1): 4-25.

[11] Williams K L, Jackson D R, Tang D J, et al. Acoustic backscattering from a sand and a sand/mud environment: experiments and data/model comparisons[J]. IEEE Journal of Oceanic Engineering, 2009, 34(4): 388-398.

[12] Holland C W, Hollett R, Troiano L. Measurement technique for bottom scattering in shallow water[J]. The Journal of the Acoustical Society of America, 2000, 108(3): 997-1011.

[13] Hines P C, Osler J C, MacDougald D J. Acoustic backscatter measurements form littoral seabeds at shallow grazing angles at 4 and 8 kHz[J]. The Journal of the Acoustical Society of America, 2005, 117(6): 3504-3516.

[14] La H, Choi J W. 8-kHz bottom backscattering measurements at low grazing angles in shallow water[J]. The Journal of the Acoustical Society of America, 2010, 127(4): 160-165.

[15] Manik H M, Furusawa M, Amakasu K. Measurement of sea bottom surface backscattering strength by quantitative echo sounder[J]. Fisheries Science, 2006, 72(3): 503-512.

[16] 金國亮, 吳承義, 張國華, 等. 淺海二維海底散射系數的測量[J]. 聲學學報, 1987, 12(3):227-231.

Jin Guoliang, Wu Chengyi, Zhang Guohua, et al. The measurement of two-dimensional bottom backscattering coefficients at shallow water[J]. Acta Acustica, 1987, 12(3):227-231.

[17] 宋磊. 海底散射系數測量方法研究[D]. 哈爾濱: 哈爾濱工程大學, 2007.

Song Lei. Study on measurement method of ocean bottom acoustic scattering coefficient[D]. Harbin: Harbin Engineering University, 2007.

[18] 薛婷. 基于T型乘積陣的海底散射系數測量方法研究[D]. 哈爾濱: 哈爾濱工程大學, 2008.

Xue Ting. Study on measurement method of ocean bottom acoustic scattering coefficient based on T-type multiplicative arrays[D]. Harbin: Harbin Engineering University, 2008.

[19] 曹正良,杜栓平,周士弘,等. 平面海底界面上球體目標的聲散射建模研究[J]. 地球物理學報, 2010, 53(2): 401-410.

Cao Zhengliang, Du Shuanping, Zhou Shihong, et al. Modeling of acoustic scattering by sphere on a planar seabed[J]. Chinese Journal of Geophysics, 2010, 53(2): 401-410.

[20] McKinney C M, Anderson C D. Measurements of backscattering of sound from the ocean bottom[J]. The Journal of the Acoustical Society of America, 1964, 36(1): 158-163.

[21] Boehme H, Chotiros N P, Rolleigh L D, et al. Acoustic backscattering at low grazing angles from the ocean bottom. Part I. Bottom backscattering strength[J]. The Journal of the Acoustical Society of America, 1985, 77(3): 962-974.

[22] Stanic S, Briggs K B, Fleischer P, et al. Shallow-water high-frequency bottom scattering off Panama City, Florida[J]. The Journal of the Acoustical Society of America, 1988, 83(6): 2134-2144.

[23] Stanic S, Briggs K B, Fleischer P, et al. High-frequency acoustic backscattering from a coarse shell ocean bottom[J]. The Journal of the Acoustical Society of America, 1989, 85(1): 125-136.

[24] Bunchuk A V, Zhitkovskii Y Y. Sound scattering by the ocean bottom in shallow-water region(review)[J]. Soviet Physics Acoustics, 1980, 26:363-370.

[25] Jackson D R, Winebrenner D P, Ishimaru A. Application of the composite roughness model to high-frequency bottom backscattering[J]. The Journal of the Acoustical Society of America, 1986, 79(5): 1410-1422.

[26] Tang D J, Frisk G V, Sellers C J, et al. Low-frequency acoustic backscattering by volumetric inhomogeneities in deep-ocean sediments[J]. The Journal of the Acoustical Society of America, 1995, 98(1): 508-518.

[27] Tang D J, Jin G L, Jackson D R, et al. Analyses of high-frequency bottom and subbottom backscattering for two distinct shallow environments[J]. The Journal of the Acoustical Society of America, 1994, 96(5): 2930-2936.

[28] Chu D Z, Williams K L, Tang D J, et al. High-frequency bistatic scattering by sub-bottom gas bubbles[J]. The Journal of the Acoustical Society of America, 1997, 102(2): 806-814.

[29] Stanic S, Kennedy E, Ray R I. Variability of shallow-water bistatic bottom backscattering[J]. The Journal of the Acoustical Society of America, 1991, 90(1): 547-553.

[30] Williams K L, Jackson D R. Bistatic bottom scattering: Model, experiments, and model/data comparison[J]. The Journal of the Acoustical Society of America, 1998, 103(1): 169-181.

[31] Briggs K B, Tang D J, Williams K L. Characterization of interface roughness of rippled sand off fort Walton Beach, Florida[J]. IEEE Journal of Oceanic Engineering, 2002, 27(3): 505-514.

[32] Wang C C, Hefner B T, Tang D J. Evaluation of laser scanning and stereo photography roughness measurement systems using a realistic model seabed surface[J]. IEEE Journal of Oceanic Engineering, 2009, 34(4): 466-475.

[33] Hines P C. Theoretical model of acoustic backscatter from a smooth seabed[J]. The Journal of the Acoustical Society of America, 1990, 88(1): 324-334.

[34] Lyons A P, Anderson A L, Dwan F S. Acoustic scattering from the seafloor∶modeling and data comparison[J]. The Journal of the Acoustical Society of America, 1994, 95(5): 2441-2451.

[35] Jackson D R, Ivakin A N. Scattering from elastic sea beds: First-order theory[J]. The Journal of the Acoustical Society of America, 1998, 103(1): 336-345.

[36] Soukup R J, Gragg R F. Backscatter from a limestone seafloor at 2-3.5 kHz: measurements and modeling[J]. The Journal of the Acoustical Society of America, 2003, 113(5): 2501-2514.

[37] 胡劍忠. 海底三維散射強度的測量和模型計算[D]. 哈爾濱: 哈爾濱工程大學, 2009.

Hu Jianzhong. Measurement and model calculation of sea bottom three-dimensional scattering strength[D]. Harbin: Harbin Engineering University, 2009.

[38] Ellis D D. Bistatic reverberation calculations using a three-dimensional scattering function[J]. The Journal of the Acoustical Society of America, 1991, 89(5): 2207-2214.

[39] Caruthers J W, Novarini J C. Modeling bistatic bottom scattering strength including a forward scatter lobe[J]. IEEE Journal of Oceanic Engineering, 1993, 18(2): 100-107.

[40] Moustier C D. Beyond bathymetry: mapping acoustic backscattering from the deep seafloor with sea beam[J]. The Journal of the Acoustical Society of America, 1986, 79(2): 316-331.

[41] Dacol D K. The kirchhoff approximation for acoustic scattering from a rough fluid-elastic solid interface[J]. The Journal of the Acoustical Society of America, 1990, 88(2): 978-983.

[42] Jackson D R, Richardson M D. High-frequency seafloor acoustics[M]. New York: Springer Science+Business Media, LLC, 2007.

[43] Kuo E Y. Wave scattering and transmission at irregular surfaces[J]. The Journal of the Acoustical Society of America, 1964, 36(11): 2135-2142.

[44] Jackson D R, Briggs K B. High-frequency bottom backscattering: roughness versus sediment volume scattering[J]. The Journal of the Acoustical Society of America, 1992, 92(2): 962-977.

[45] Moe J E, Jackson D R. First-order perturbation solution for rough surface scattering cross section including the effects of gradients[J]. The Journal of the Acoustical Society of America, 1994, 96(3): 1748-1754.

[46] Kuperman W A, Schmidt H. Rough surface elastic wave scattering in a horizontally stratified ocean[J]. The Journal of the Acoustical Society of America, 1986, 79(6): 1767-1977.

[47] Kuo E Y. Acoustic wave scattering from two solid boundary at the ocean bottom: reflection loss[J]. IEEE Journal of Oceanic Engineering, 1992, 17(1): 159-170.

[48] Essen H H. Scattering from a rough sedimental seafloor containing shear and layering[J]. The Journal of the Acoustical Society of America, 1994, 95(3): 1299-1310.

[49] Williams K L, Grochocinski J M, Jackson D R. Interface scattering by poroelastic seafloors: first-order theory[J]. The Journal of the Acoustical Society of America, 2001, 110(6): 2956-2963.

[50] Berman D H. Simulation of rough interface scattering[J]. The Journal of the Acoustical Society of America, 1991, 89(2): 623-636.

[51] Yang T Q, Broschat S L. Acoustic scattering from a fluid-elastic-solid interface using the small slope approximation[J]. The Journal of the Acoustical Society of America, 1994, 96(3): 1796-1804.

[52] Gragg R F, Wurmser D, Gauss R C. Small-slope scattering from rough elastic ocean floors: general theory and computational algorithm[J]. The Journal of the Acoustical Society of America, 2001, 110(6): 2878-2901.

[53] Broschat S L, Thorsos E I. An investigation of the small slope approximation for scattering from rough surface. Part Ⅱ. Numerical studies[J]. The Journal of the Acoustical Society of America, 1997, 101(5): 2615-2625.

[54] Gragg R F, Wurmser D, Gauss R C. Small-slope scattering from rough elastic ocean floors:General theory and computational algorithm[J]. The Journal of the Acoustical Society of America, 2001, 110(6): 2878-2901.

[55] Soukup R J, Canepa G, Simpson H J, et al. Small-slope simulation of acoustic backscatter from a physical model of an elastic ocean bottom[J]. The Journal of the Acoustical Society of America, 2007, 122(5): 2551-2559.

[56] Jackson D. The small-slope approximation for layered seabeds[J]. The Journal of the Acoustical Society of America, 2013, 133(5): 3251.

[57] Thorsos E I. Acoustic scattering from a “Pierson-Moskowitz” sea surface[J]. The Journal of the Acoustical Society of America, 1990, 88(1): 335-349.

[58] Jackson D R, Odom R I, Boyd M L, et al. A geoacoustic bottom interaction model(GABIM)[J]. IEEE Journal of Oceanic Engineering, 2010, 35(3): 603-617.

[59] 張仁和, 金國亮. 淺海平均混響的簡正波理論[J]. 聲學學報, 1984, 9(1): 12-20.

Zhang Renhe, Guo Jinliang. Normal-mode theory of the average reverberation intensity in shallow water[J]. Acta Acustica, 1984,9(1):12-20.

[60] 李風華,金國亮,張仁和. 淺海相干混響理論與混響強度的振蕩現象[J]. 中國科學:A輯, 2000, 30(6):560-566.

Li Fenghua, Jin Guoliang, Zhang Renhe. Correlation reverberation theory in shallow water and oscillation phenomena of reverberation intensity [J]. Science in China (Series A), 2000,30(6):560-566.

[61] 彭朝暉,周紀潯,張仁和. 非均勻海底和粗糙度界面引起的平面內海底散射[J]. 中國科學:G輯, 2004, 34(4):378-391.

Peng Zhaohui, Zhou Jixun, Zhang Renhe. Bottom scattering in a plane from heterogeneous ocean bottom and rough interface[J]. Science in China (Series G), 2004, 34(4): 378-391.

A review on the progress of bottom acoustic scattering research

Liu Baohua1,2, Kan Guangming2,3, Pei Yanliang2,Yang Zhiguo1,Yu Kaiben1,2, Yu Shengqi1

(1.NationalDeepSeaCenter,StateOceanicAdministration,Qingdao266061,Chnia; 2.LaboratoryforMarineGeology,QingdaoNationalLaboratoryforMarineScienceandTechnology,Qingdao266237,China; 3.KeyLaboratoryofMarineSedimentologyandEnvironmentalGeology,TheFirstInstituteofOceanography,StateOceanicAdministration,Qingdao266237,China)

Seafloor is an important boundary for underwater sound field, and its acoustic scattering property has significant influence to spatial structure and distribution of underwater sound field. In this paper, we systematically analyzed and summarized the present progress on the measuring techniques, mechanism and modeling of bottom acoustic scattering property in the world. In addition, some research directions and topics which should be paid more attention to in the future were put forwarded as well as some research difficulties. This work will have great contribution to further studying and understanding the present status and the future trend of the research on bottom acoustic scattering property.

bottom acoustic scattering property; scattering model; bottom reverberation; underwater sound field

2016-08-04;

2016-12-12。

國家自然科學基金(41330965, 41527809);青島海洋國家實驗室開放基金(QNLM20160RP0209);海洋公益性科研專項資金項目(201405032);泰山學者工程專項經費項目(TSPD20161007)。

劉保華(1960-),男,山東省嘉祥縣人,研究員,博士生導師,主要從事海洋地球物理研究。E-mail:bhliu@fio.org.cn

*通信作者:闞光明,男,副研究員,主要從事海底聲學研究。E-mail:kgming135@fio.org.cn

10.3969/j.issn.0253-4193.2017.07.001

P733.23;P736.21

A

0253-4193(2017)07-0001-11

劉保華, 闞光明, 裴彥良, 等. 海底聲散射特性研究進展[J]. 海洋學報, 2017, 39(7): 1-11,

Liu Baohua, Kan Guangming, Pei Yanliang, et al. A review on the progress of bottom acoustic scattering research[J]. Haiyang Xuebao, 2017, 39(7): 1-11, doi:10.3969/j.issn.0253-4193.2017.07.001

猜你喜歡
界面測量模型
一半模型
重要模型『一線三等角』
國企黨委前置研究的“四個界面”
當代陜西(2020年13期)2020-08-24 08:22:02
重尾非線性自回歸模型自加權M-估計的漸近分布
把握四個“三” 測量變簡單
滑動摩擦力的測量和計算
滑動摩擦力的測量與計算
基于FANUC PICTURE的虛擬軸坐標顯示界面開發方法研究
人機交互界面發展趨勢研究
3D打印中的模型分割與打包
主站蜘蛛池模板: 91精品最新国内在线播放| 亚洲区欧美区| 亚洲精品无码日韩国产不卡| 四虎国产永久在线观看| 久久人妻xunleige无码| 色天堂无毒不卡| 99精品在线视频观看| 国产欧美日韩专区发布| 波多野结衣视频一区二区| 青青操国产视频| 伊人网址在线| 国产精品成人一区二区不卡| 成年A级毛片| 日本国产在线| 亚洲中文字幕在线观看| 色男人的天堂久久综合| 欧美日韩综合网| 日韩色图在线观看| 久久久噜噜噜久久中文字幕色伊伊| 久无码久无码av无码| av免费在线观看美女叉开腿| 视频二区中文无码| 国产欧美日韩视频怡春院| 91亚洲免费| www亚洲天堂| 国产香蕉97碰碰视频VA碰碰看| 无码av免费不卡在线观看| 超碰91免费人妻| 国产无码网站在线观看| 精品1区2区3区| 片在线无码观看| 国产精品免费电影| 99re热精品视频国产免费| 欧美自慰一级看片免费| 青青青视频蜜桃一区二区| 亚洲国产精品久久久久秋霞影院| 国产精品自在拍首页视频8| 超清无码熟妇人妻AV在线绿巨人| 任我操在线视频| 全部免费毛片免费播放| 国产一级二级三级毛片| 91午夜福利在线观看| 欧美精品在线看| 久久夜色精品国产嚕嚕亚洲av| 成人免费一级片| 国产黄色免费看| 97无码免费人妻超级碰碰碰| 欧美综合中文字幕久久| 内射人妻无套中出无码| 97超碰精品成人国产| 国产福利不卡视频| 精品国产成人av免费| 日本高清成本人视频一区| 国产视频资源在线观看| 成人国产三级在线播放| 亚洲精品第五页| 日本91视频| a级毛片网| 亚洲综合网在线观看| 亚洲综合极品香蕉久久网| 啊嗯不日本网站| 国产精品部在线观看| 欧美特黄一免在线观看| 亚洲最大福利网站| 国产又黄又硬又粗| 久久毛片免费基地| 亚洲精品中文字幕无乱码| 日韩高清中文字幕| 青青青伊人色综合久久| 亚洲性影院| 国产精品成人免费视频99| 日韩精品高清自在线| 国产制服丝袜91在线| 免费在线一区| 内射人妻无套中出无码| 四虎影视8848永久精品| 亚洲专区一区二区在线观看| 国产精品 欧美激情 在线播放| 日韩a级毛片| 老司机久久99久久精品播放| 99热免费在线| 亚洲香蕉在线|