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

    海底大地基準(zhǔn)建設(shè)技術(shù)及其研究進(jìn)展

    2022-11-03 08:49:00劉焱雄李夢昊何秀鳳陳冠旭張林虎唐秋華
    海洋科學(xué)進(jìn)展 2022年4期
    關(guān)鍵詞:測量

    劉焱雄,李夢昊,劉 楊,何秀鳳,陳冠旭,張林虎,唐秋華

    (1.自然資源部 第一海洋研究所,山東 青島 266061;2.自然資源部 海洋測繪重點實驗室,山東 青島 266061;3.河海大學(xué) 地球科學(xué)與工程學(xué)院,江蘇 南京 211100)

    我國是海洋大國,海域面積十分遼闊。認(rèn)識海洋、經(jīng)略海洋,探索海洋奧秘、發(fā)展海洋經(jīng)濟(jì)、深耕藍(lán)色國土,需要建設(shè)海底觀測網(wǎng)等基礎(chǔ)設(shè)施;海底觀測網(wǎng)已是繼地面/洋面、空間之后,觀測地球系統(tǒng)的第三個平臺[1]。海底大地基準(zhǔn)網(wǎng)是海底觀測網(wǎng)的重要組成,是國家大地基準(zhǔn)網(wǎng)由陸域向海域的自然延伸[2-3],也是構(gòu)建陸??仗煲惑w化空間基準(zhǔn)的國家基礎(chǔ)設(shè)施。海洋空間基準(zhǔn)是一切海洋活動的前提和基礎(chǔ),大力發(fā)展海洋大地測量技術(shù),加快布設(shè)海底大地基準(zhǔn)網(wǎng),滿足國防安全保障和經(jīng)濟(jì)社會發(fā)展需求,對推動海洋導(dǎo)航定位技術(shù)進(jìn)步、推進(jìn)我國海洋強國戰(zhàn)略實施具有重要支撐作用。

    海底大地基準(zhǔn)網(wǎng)建設(shè)是當(dāng)今世界大國必爭的高技術(shù)戰(zhàn)略領(lǐng)域。美國、加拿大、俄羅斯等海洋強國早已開啟海底大地基準(zhǔn)網(wǎng)的研究[2,4-6],日本也建立了海底大地基準(zhǔn)網(wǎng)[7-8]。目前,我國僅僅在南海3 000 m 水深的海域開展了海底大地基準(zhǔn)網(wǎng)試驗,但尚未大規(guī)模布設(shè)海底大地基準(zhǔn)網(wǎng)[2]。研究海洋大地測量技術(shù),完善自主海底大地基準(zhǔn)基礎(chǔ)設(shè)施,對支撐我國2035 年前建成國家綜合PNT(Positioning,Navigation and Timing)系統(tǒng)意義重大。

    海底大地基準(zhǔn)網(wǎng)由若干海底大地基準(zhǔn)站組成,需要首先確定這些基準(zhǔn)站的準(zhǔn)確位置。聯(lián)合全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)和水下聲學(xué)定位系統(tǒng)(GNSS-聲學(xué)或GNSS-A),可實現(xiàn)海底大地基準(zhǔn)站的位置標(biāo)定,可以將全球統(tǒng)一的時空基準(zhǔn)傳遞到海底[3]。穩(wěn)定的海底大地基準(zhǔn)站是海底大地基準(zhǔn)網(wǎng)的基礎(chǔ)和構(gòu)成,海底大地基準(zhǔn)站建設(shè)需要歷經(jīng)站址勘選、設(shè)備布放回收、站址位置標(biāo)定等流程。鑒于此,本文聚焦海底大地基準(zhǔn)建設(shè)技術(shù),首先梳理水下聲學(xué)導(dǎo)航定位系統(tǒng)及技術(shù)背景,然后總結(jié)海底大地基準(zhǔn)建設(shè)的關(guān)鍵技術(shù)要點,最后探討GNSS-A 定位數(shù)據(jù)處理方法,并對未來海底大地基準(zhǔn)網(wǎng)建設(shè)技術(shù)及應(yīng)用進(jìn)行了展望。

    1 水聲導(dǎo)航定位系統(tǒng)及技術(shù)

    1.1 水下聲學(xué)導(dǎo)航定位系統(tǒng)

    水下聲學(xué)定位系統(tǒng)依據(jù)聲學(xué)單元的距離分為超短基線系統(tǒng)(<1 m)、短基線系統(tǒng)(1~50 m)、長基線(100~6 000 m)系統(tǒng)等類型[9]。國外聲學(xué)定位系統(tǒng)的研發(fā)較為成熟,產(chǎn)品呈現(xiàn)多樣化和系統(tǒng)化。挪威Kongsberg Simrad 公司研發(fā)了HiPAP 系列和μPAP 系列水下聲學(xué)定位系統(tǒng),作用距離可達(dá)10 000 m,測距精度達(dá)到0.02 m[10];法國iXBlue 公司的GAPS 系列產(chǎn)品集成了慣導(dǎo)設(shè)備,定位精度可達(dá)到0.2%斜距,Posidonia 超短基線定位系統(tǒng)最大作用距離超過10 000 m[11];英國Sonardyne 公司研發(fā)的水下聲學(xué)定位系統(tǒng),如Fusion 系列、Ranger 系列、Scout 系列以及Marksman 系列等,囊括了長基線、短基線、超短基線定位系統(tǒng)以及組合定位系統(tǒng)[12]。

    美國、俄羅斯等國家已開啟新型水下導(dǎo)航系統(tǒng)研發(fā)[13]。2015 年,美國國防高級研究計劃局(DARPA)提出構(gòu)建“深海導(dǎo)航定位系統(tǒng)”,開始研究在海床上安裝聲學(xué)信號源,組成類似全球定位系統(tǒng)(GPS)的水下GPS;2016 年,美國在菲律賓海開展海洋聲學(xué)深水計劃,驗證水下GPS;俄羅斯也開展了水下導(dǎo)航定位系統(tǒng)研制,并通過低頻、被動接收的水聲定位方式進(jìn)行水下導(dǎo)航誤差校準(zhǔn)。

    我國水聲定位技術(shù)研究和系統(tǒng)研發(fā)起步較晚,目前進(jìn)入快速發(fā)展期。哈爾濱工程大學(xué)、中國科學(xué)院聲學(xué)研究所、西北工業(yè)大學(xué)和自然資源部第一海洋研究所等單位對水下聲學(xué)導(dǎo)航定位技術(shù)進(jìn)行了深入研究[6,9,14],經(jīng)過近20 a 的努力,國內(nèi)具備了全海深高精度聲學(xué)導(dǎo)航定位能力。在國家“十一五”“國家高技術(shù)研究發(fā)展計劃”的支持下,中國測繪科學(xué)研究院和中船重工715 研究所研制了“水下GPS 高精度導(dǎo)航定位系統(tǒng)”。哈爾濱工程大學(xué)牽頭、自然資源部第一海洋研究所合作研制了“長程超短基線定位系統(tǒng)”,工作水深超過3 700 m,作用距離達(dá)到8.6 km,定位精度為0.2%~0.3%斜距,達(dá)到國際先進(jìn)水平,并獲得國家技術(shù)發(fā)明二等獎。

    1.2 GNSS-A 定位技術(shù)研究背景及進(jìn)展

    全球?qū)Ш叫l(wèi)星系統(tǒng)和聲學(xué)組合的定位方法(GNSS-A)由美國斯克利普斯海洋研究所(Scripps Institution of Oceanography,SIO)率先提出[4,15],是一種將海面平臺上的動態(tài)GNSS 定位和水下聲學(xué)測距(海面平臺與海底聲學(xué)信標(biāo)之間)相組合的技術(shù)。通過GNSS-A 定位,可以實現(xiàn)國際橢球參考框架(International Terrestrial Reference Frame,ITRF)下海底信標(biāo)的位置測定(圖1)。美國Scripps 海洋研究所采用GPS-A(GNSS-A 的早期形式)定位系統(tǒng)同時對3 個海底信標(biāo)進(jìn)行聲學(xué)測量,24~48 h 連續(xù)測量的點位精度達(dá)到厘米級,80 h 連續(xù)測量的點位誤差則小于1 cm[16-17]。日本東京大學(xué)于1987 年首次嘗試在Sagami 灣進(jìn)行海底定位實驗[18]。在20世紀(jì)90 年代中期,日本海上保安廳海洋水文部開始研發(fā)GPS-A 定位系統(tǒng),2000 年在日本南部海槽2 000 m水深的熊野盆地部署了海底大地基準(zhǔn)點[19]。目前,日本已采用GNSS-A 定位系統(tǒng),獲取了許多重要的海底大地測量結(jié)果,包括探測地震之間的關(guān)系、同震以及震后與地震周期有關(guān)的俯沖區(qū)變形,和山脊變形邊界附近的板塊運動等[5,20]。近年來,在國家重點研發(fā)計劃項目“海洋大地測量基準(zhǔn)與海洋導(dǎo)航新技術(shù)”的支持下,楊元喜院士帶領(lǐng)國內(nèi)相關(guān)科研團(tuán)隊,開展了海洋大地基準(zhǔn)的技術(shù)攻關(guān),采用GNSS-A 進(jìn)行了海底大地基準(zhǔn)站位置標(biāo)校試驗,并布設(shè)了3 000 m 水深的海底大地基準(zhǔn)點,實現(xiàn)了分米級精度的海底定位和米級精度的水下聲學(xué)導(dǎo)航[2],取得了系列研究成果。

    圖1 GNSS-A 定位原理Fig.1 Schematic of GNSS-A positioning

    目前,我國GNSS-A 定位技術(shù)方面主要集中于定位模型和算法研究,包括誤差處理、隨機模型和函數(shù)模型優(yōu)化等。為了提高定位效率,Yang 等[21]、Zhao 等[22]研究了利用聲速剖面的先驗信息反演聲速進(jìn)行定位解算的方法;孫文舟等[23]利用經(jīng)驗正交函數(shù)(Empirical Orthogonal Function,EOF)反演了聲速剖面信息,張林虎等[24]研究了基于分層EOF 函數(shù)的區(qū)域聲速場模型構(gòu)建方法。為了提高定位精度,Chen 等[25]針對桿臂矢量垂向偏心誤差與海底基準(zhǔn)點高程間的耦合問題,提出了基于樣本搜索法的偏心誤差解算方 法;Chen 等[26]、Chen 等[27]、趙建虎等[28]、曾安敏等[29]分析了圓形航跡的定位優(yōu)勢。針對聲速結(jié)構(gòu)時空變化,Yang 等[30]結(jié)合Fujita 等[31]方法,從定位殘差中擬合聲速長期時域變化項、聲速周期性變化項;劉楊等[32-33]通過構(gòu)建聲速時變引起的水聲垂向總時延,聯(lián)合估計了聲速時域變化和海底大地基準(zhǔn)點位置,估計改正的聲速與實測聲速的偏差小于0.13 m/s,提高了海底大地基準(zhǔn)點坐標(biāo)的準(zhǔn)確度;Wang 等[34]采用兩步估計法,分別估計了海底大地基準(zhǔn)點的位置、系統(tǒng)誤差,以及聲速長周期項誤差。關(guān)于隨機模型優(yōu)化,Zhao 等[35-36]、王薪普等[37]提出了與聲線入射角相關(guān)的分段指數(shù)隨機模型,馬越原等[38]指出了海洋環(huán)境的復(fù)雜性影響隨機模型的效果??紤]函數(shù)模型對定位的影響,Li 等[39]、閆鳳池等[40]研究了基于單向、雙向聲學(xué)傳播時間的定位觀測方程;鄺英才等[41-42]研究了聯(lián)合解算船載換能器與海底大地基準(zhǔn)點位置的方法;Xu 等[43]基于GNSS 差分定位思路,提出了水下定位差分觀測方程;Chen 等[26]、孫文舟等[44]研究了附加深度約束的差分定位算法;Xue 等[43]進(jìn)一步證明了水下定位差分解和非差解的等價性。趙建虎等[45-46]提出了基于絕對標(biāo)校+相對測量的海底網(wǎng)平差思路,實現(xiàn)海底大地基準(zhǔn)網(wǎng)的整體解算。面對復(fù)雜海洋環(huán)境,為了增加海底大地基準(zhǔn)測量的可用性、可靠性、精確性,需要進(jìn)一步完善定位算法,建立海洋場景的自適應(yīng)彈性隨機模型、彈性函數(shù)模型[30]。

    2 海底大地基準(zhǔn)建設(shè)關(guān)鍵技術(shù)

    海底大地基準(zhǔn)站包含水下潛標(biāo)和海床基兩類。水下潛標(biāo)基站一般由浮球、觀測設(shè)備、纜繩、聲學(xué)釋放器、配重等組成[47-48]。隨著潛標(biāo)離海底高度的增加,可以有效擴大基準(zhǔn)站的服務(wù)范圍。由于水下潛標(biāo)隨著海水運動而晃動,潛標(biāo)通常需要搭載壓力計、姿態(tài)傳感器等設(shè)備以實時修正其位置,因此這種類型的海底基準(zhǔn)站技術(shù)難度大、位置精度低。海床基是主要的海底基準(zhǔn)站布放方式,具有穩(wěn)定性好、集成度高等優(yōu)勢[49]。海床基基站通常包括坐底平臺(包括配重)、聲學(xué)釋放器、水下電源、回收浮體,以及搭載的各種傳感器等。常見的淺海海床基基站外觀主要以圓形和多邊形為主,圓形尺寸一般以直徑1~6 m 為主,如加拿大“海王星”海底觀測網(wǎng)(NEPTUNE)使用的海床基直徑約6 m;意大利PROTECOSUB 公司的海床基以圓形結(jié)構(gòu)為主,美國伍茲霍爾海洋研究所研制的海床基則以多邊形結(jié)構(gòu)為主。深海海床基基站通常配備支撐框架,框架中間部分集成有觀測儀器、電池倉等設(shè)備,框架支撐腿與腳部配重盤則作為坐底的配重。

    穩(wěn)定的海底大地基準(zhǔn)站是海底大地基準(zhǔn)網(wǎng)的基礎(chǔ),“放得穩(wěn)、測得準(zhǔn)、待得久”是海底大地基準(zhǔn)站的建設(shè)目標(biāo)。“待得久”主要涉及電源能量供應(yīng)問題,這里暫不討論;本文按照“放得穩(wěn)”和“測得準(zhǔn)”的要求,從海底基準(zhǔn)的站址勘選、布放回收、位置標(biāo)校等關(guān)鍵技術(shù),討論海底大地基準(zhǔn)建設(shè)的研究進(jìn)展。

    2.1 站址勘選

    海底大地基準(zhǔn)站的站址尤為關(guān)鍵。由于海底環(huán)境復(fù)雜,極易造成基準(zhǔn)站失穩(wěn),導(dǎo)致海底大地基準(zhǔn)網(wǎng)失準(zhǔn),進(jìn)而影響水下導(dǎo)航定位精度。為此,海底大地基準(zhǔn)站的站址選擇一般遵循大范圍(面)→局部區(qū)域 →布放點(點)的總原則,需要開展桌面研究和現(xiàn)場勘測等工作,涉及到海底地形、地貌、底質(zhì)、水文、聲學(xué)等觀測技術(shù)。由于水深地形、地貌底質(zhì)、水文環(huán)境調(diào)查屬于通用技術(shù),這里不再贅述,只介紹相關(guān)工作和技術(shù)要求。

    2.1.1 桌面研究

    收集歷史調(diào)查資料,包括水深地形地貌、工程地質(zhì)、水文環(huán)境、聲速環(huán)境等資料,分析海底地形地貌特征,考慮海底淺表層底質(zhì)分布及其工程地質(zhì)特征,了解海底底流和聲學(xué)傳播狀況,獲取海底平整度、淺表地層活動性與穩(wěn)定性,以及水動力特征,同時依據(jù)海底大地基準(zhǔn)站的外部形態(tài)、重量、工作方式與性能指標(biāo)[2,50]。初步甄選海床自身穩(wěn)定性較好的局部區(qū)域作為海底大地基準(zhǔn)站布放海區(qū),保證站址選擇的合理性與科學(xué)性。桌面研究主要考慮如下因素。

    1)地形坡度

    收集水深地形數(shù)據(jù),生成海區(qū)坡度DEM(Digital Elevation Model)和坡度等值線。以3°和5°為閾值進(jìn)行劃分,按坡度將海區(qū)劃分為小于3°區(qū)域、3°~5°區(qū)域、大于5°區(qū)域。盡量選取坡度小于3°海區(qū)。

    2)動力地貌類型

    選取海區(qū)海底地貌應(yīng)具有地形平坦、起伏小等特征,如深海平原地貌,其底部流場相對穩(wěn)定,且具有地面平坦或無起伏等特征。

    3)表層沉積物承載力

    選取海區(qū)表層沉積物類型應(yīng)以粉砂和黏土為主;以粉砂(包含粉質(zhì)黏土和黏土質(zhì)粉砂等類型)的承載力滿足地層穩(wěn)定和承載力條件。

    4)地質(zhì)災(zāi)害

    不管是具有活動能力的破壞性地質(zhì)災(zāi)害,還是不具有活動能力的地質(zhì)災(zāi)害,都是對海底穩(wěn)定性造成影響的潛在因素。因此,海區(qū)的選擇要避開各種海底地質(zhì)災(zāi)害的分布區(qū)域,如斷層、滑坡、沖刷槽、海底峽谷、沙波沙丘、麻坑、淺層氣、易液化砂層、軟弱地層等。

    5)試驗成本

    除了上述從自然環(huán)境條件進(jìn)行布放位置的考慮,同時從交通便利性角度出發(fā),對距離和通航條件進(jìn)行分析。考慮試驗成本,確定試驗海區(qū)。

    2.1.2 現(xiàn)場勘測及定址

    現(xiàn)場勘測采用多波束測深系統(tǒng)(Multi-beam Echo Sounder System)、側(cè)掃聲吶儀(Sidescan sonar)、淺地層剖面儀(Sub-bottom Profiler)等高分辨聲學(xué)探測系統(tǒng),利用高頻和低頻聲吶,對目標(biāo)海域進(jìn)行全覆蓋水深、海底地形地貌測量、工程地質(zhì)調(diào)查和水文環(huán)境調(diào)查,查明目標(biāo)海域的詳細(xì)環(huán)境信息[50]。其中,多波束測深系統(tǒng)從海面換能器發(fā)射聲波束條帶,通過水聽器接收海底反射回波,以獲取高精度水深信息、海底地形數(shù)據(jù),如德國L-3 ELAC Nautik 公司的SeaBeam 3012 全海深型多波束,發(fā)射頻率為12 kHz,測量水深為50~11 000 m[51-52];多波束測深系統(tǒng)還能獲取海底的反向散射強度數(shù)據(jù),用于海底底質(zhì)分類[53]。側(cè)掃聲吶的發(fā)射頻率范圍為50~500 kHz,可用于海底地形地貌測量[54]。淺地層剖面儀的發(fā)射頻率較低,具有較強的穿透力,能夠有效地穿透海底數(shù)十米的地層,通過反射波的走時、振幅、頻率等信息,連續(xù)探測水下淺部地層結(jié)構(gòu)、構(gòu)造和底質(zhì)等信息[55]。

    勘測期間,配合多波束測深、側(cè)掃聲吶掃測等需要,采用聲學(xué)多普勒流速剖面儀(Acoustic Doppler Current Profiler,ADCP)、重力取樣儀等進(jìn)行海洋流場測量、底質(zhì)取樣調(diào)查也是重要的勘測內(nèi)容。海底大地基準(zhǔn)站的背景聲速場是站址勘選需要考慮的重要因素,需要開展聲速剖面現(xiàn)場測量。

    對現(xiàn)場勘測數(shù)據(jù)進(jìn)行精細(xì)處理和系統(tǒng)分析,獲得調(diào)查區(qū)域的水深地形地貌、沉積物類型、工程地質(zhì)特征、水文環(huán)境和聲速環(huán)境等成果。根據(jù)這些成果,選取海底起伏度、海底坡度、沉積物類型分布、構(gòu)造穩(wěn)定性、工程地質(zhì)力學(xué)性質(zhì)、地層沖淤特征與地層穩(wěn)定性等作為站址勘選主要評價指標(biāo)。在遵循選址總原則的基礎(chǔ)上,基準(zhǔn)站需布放安全、運行穩(wěn)定,滿足作業(yè)要求,由此確定海底大地基準(zhǔn)站站址勘選的5 級標(biāo)準(zhǔn),其中1 級最優(yōu)、3 級適宜、5 級最差(表1)。

    表1 海底大地基準(zhǔn)站站址勘選指標(biāo)集與評價標(biāo)準(zhǔn)Table 1 Index set and evaluation criteria of station site survey

    2.2 布放與回收

    海底大地基準(zhǔn)站的布放包括直接投放方法、吊裝布放法、ROV(Remote Operated Vehicle)布放法等。直接投放方法針對淺海及簡單水下環(huán)境的場景,通常在貼近水面上方將脫鉤打開,利用基準(zhǔn)站自重下沉到海底[56]。該方法操作簡單、投放速度快,但受海水動力影響,基準(zhǔn)站的實際坐底地點與預(yù)期位置容易出現(xiàn)較大偏差。吊裝布放法是在調(diào)查船甲板由絞車進(jìn)行吊裝布放[48,56](圖2a);布放時絞車以0.5~1.0 m/s 的速度下放基準(zhǔn)站,為避免觸碰海底,基準(zhǔn)站距海底約10~20 m 時,由聲學(xué)釋放器甲板單元向水下發(fā)送釋放指令,命令聲學(xué)釋放器脫鉤,基準(zhǔn)站自由落體坐底。吊裝布放法容易受到涌浪及海流的共同作用,吊放狀態(tài)下基準(zhǔn)站的擺動幅度可達(dá)十多米,無法實現(xiàn)精準(zhǔn)定點布放。ROV 布放法可將海底大地基準(zhǔn)站精準(zhǔn)布放于預(yù)定位置[57],該方法首先通過ROV 將海底大地基準(zhǔn)站運送到近海底位置并懸停,然后操作ROV 移動到目標(biāo)位置后釋放海底大地基準(zhǔn)站。加拿大科學(xué)潛水設(shè)施的海洋科學(xué)遙控操作平臺(Remotely Operated Platform for Ocean Science,ROPOS),采用ROV 布放法,先后實現(xiàn)了加拿大和美國的觀測網(wǎng)設(shè)備的精確定位布放。ROV 布放法操作復(fù)雜,且需要多種設(shè)備同時作業(yè),成本較高。

    圖2 海底大地基準(zhǔn)站布放示例Fig.2 Example of seafloor geodetic station deployment

    海底大地基準(zhǔn)站回收主要有2 種方式[48-49]。一種方式是水面甲板單元向海床基聲學(xué)釋放器發(fā)送釋放指令,海床基釋放器打開,丟掉配重,浮體攜帶儀器上浮海面并打撈回收。另一種方式是利用浮體攜帶的回收繩索將海床基打撈回收,這種方式只適合于淺海,通常,儀器設(shè)備回收后可拋棄海床基。

    2.3 位置標(biāo)校

    位置標(biāo)校即獲得海底大地基準(zhǔn)站的絕對位置,需要通過GNSS-A 定位系統(tǒng)將陸地的絕對基準(zhǔn)傳遞到海底大地基準(zhǔn)站。采用海床基方式布放的海底大地基準(zhǔn)站,一般需在海床基上搭載一個聲學(xué)應(yīng)答器,應(yīng)答器的相位中心即為海底大地基準(zhǔn)站位置(圖2b)。另外,海底大地基準(zhǔn)站也需配備電池以供長期、穩(wěn)定地運行[58],也可搭載其他傳感器,如自容式溫鹽深測量儀、壓力傳感器等設(shè)備[59]。海底大地基準(zhǔn)站在布放時,一般按照組網(wǎng)的方式,通常由3~6 個基準(zhǔn)站組成三角形或方形的基準(zhǔn)網(wǎng)[60-62];位置標(biāo)校時通常采用圓形標(biāo)校線路,圓形經(jīng)驗半徑一般采用1.414 倍水深。

    2.3.1 GNSS-A 定位系統(tǒng)

    GNSS-A 定位系統(tǒng)主要包括GNSS 天線/接收機、水下聲學(xué)傳感器、姿態(tài)傳感器等[60,63]。位置標(biāo)校首先通過秒脈沖信號或GPS 時間完成GNSS-A 定位系統(tǒng)的時間同步,然后使用GNSS 測量海面平臺上GNSS 天線的位置,利用GNSS 天線和聲學(xué)換能器之間的相對位置,以及海面平臺的姿態(tài)(航向、俯仰和橫滾角),確定海面平臺上聲學(xué)傳感器的位置,接著基于聲學(xué)傳感器發(fā)射聲學(xué)信號并接收來自海底聲學(xué)應(yīng)答器的反射信號,獲得海面聲學(xué)傳感器和海面應(yīng)答器之間的往返傳播時間(Time of Flight,TOF),最后,通過最小化TOF觀測值和計算值之間的偏差,估計海底應(yīng)答器的位置。

    GNSS-A 定位系統(tǒng)使用的海面觀測平臺種類較多,主要包括調(diào)查船、無人艇、浮標(biāo)等觀測平臺[64-67]。

    1)調(diào)查船觀測平臺

    調(diào)查船觀測平臺根據(jù)觀測方式主要包括船舷懸掛方式、船底固定方式兩種[60]。船舷懸掛方式將設(shè)備安裝在船舷的一根桿子上,GNSS 天線和姿態(tài)傳感器安裝在桿的頂部,聲學(xué)換能器安裝在桿的底部。每次采集數(shù)據(jù)時,連接桿需要單獨安裝;連接聲學(xué)換能器的一端放入海水中,為避免干擾、保證觀測質(zhì)量,需盡量遠(yuǎn)離螺旋槳和發(fā)動機,遠(yuǎn)離船底1 m 以上。為了避免船只螺旋槳引起的噪音和水中負(fù)荷過大而導(dǎo)致的桿子變形,要求在船只漂流或低速運行時進(jìn)行聲學(xué)測量[31]。船舷懸掛方式需要花費2~4 d,才能夠獲得足夠多的聲學(xué)觀測數(shù)據(jù)。船底固定方式將設(shè)備固定安裝在船上,如聲學(xué)傳感器安裝在船底月池內(nèi),GNSS 天線安裝在主桅桿頂部。船底固定方式能夠在船只沿預(yù)定航跡航行的情況下進(jìn)行聲學(xué)觀測,并且觀測時間可縮短至16~24 h。為了進(jìn)一步提高觀測頻次、減少觀測時間,開發(fā)了多聲測距換能器,可在一個觀測序列中進(jìn)行多次信號發(fā)射和接收[60]。這種新的換能器連續(xù)發(fā)射聲學(xué)信號,一次性接收所有的回波信號,可以在3~4 h 內(nèi)完成聲學(xué)觀測。

    調(diào)查船觀測平臺能達(dá)到厘米級海底定位精度,可以獲取高分辨率、長周期的海底基準(zhǔn)站位置變化,但需投入大量資金和人力,且作業(yè)期間一旦船體安裝的換能器發(fā)生故障,無法及時進(jìn)入船塢修復(fù),使得長期連續(xù)觀測變得困難[68]。由于觀測頻次的不足,調(diào)查船觀測平臺難以檢測到海底基準(zhǔn)站短期位置變化,需要提高GNSS-A 觀測的時間分辨率,理想解決方案是使用海面無人平臺進(jìn)行連續(xù)觀測,如浮標(biāo)、海面無人艇等[58,65-66,69]。

    2)無人艇觀測平臺

    Kido 等[65]測試了使用海面無人艇(Autonomous Surface Vehicle,ASV)搭載的GNSS-A 定位系統(tǒng),其配備導(dǎo)航和聲學(xué)測距系統(tǒng)的ASV,可以獲得與調(diào)查船觀測平臺同等質(zhì)量和精度的結(jié)果,有助于降低調(diào)查成本、增加GNSS-A 觀測頻次[70-71]。與調(diào)查船觀測平臺相比,ASV 動力系統(tǒng)由大型電池提供電力,并由船上的柴油發(fā)電機充電,推進(jìn)器噪音可以忽略不計;且ASV 在定點位置可保持3 m 以內(nèi),沿著預(yù)定航跡航行的路徑差異小于2 m。因此,ASV 可作為GNSS-A 定位的候選平臺。Sakic 等[72]將GNSS-A 定位系統(tǒng)安裝在美國L3 Harris 公司設(shè)計的小型雙體船上,系統(tǒng)主要包括GNSS 天線/接收機,超短基線模塊,并集成慣性系統(tǒng)來校正船體位置;該GNSS-A 定位系統(tǒng)對淺水應(yīng)答器的定位實現(xiàn)了5 cm 的可重復(fù)性。Linuma 等[68]證明了波浪滑翔器(Wave Glider,WG)同樣具備搭載GNSS-A 定位系統(tǒng)的能力。波浪滑翔器可以依靠太陽能電池板產(chǎn)生足夠的電力,實現(xiàn)自動導(dǎo)航、衛(wèi)星通信和GNSS-A 觀測;該GNSS-A 定位系統(tǒng)的設(shè)計目標(biāo)是進(jìn)行實時數(shù)據(jù)處理,目前,仍需完善聲學(xué)傳播時間的實時檢測問題[73]。無人艇觀測平臺受海面風(fēng)浪的影響較大,平臺姿態(tài)精密測量、及由此帶來的聲信號傳播與定位的影響也是需要考慮的問題。

    3)浮標(biāo)觀測平臺

    錨系浮標(biāo)作為海面平臺,能利用GNSS-A 定位系統(tǒng)對海底地殼形變進(jìn)行連續(xù)的實時測量[66]。Imano 等[67]驗證了使用浮標(biāo)進(jìn)行3 000 m 深度海底定位的準(zhǔn)確性,但是由于浮標(biāo)是錨系的,在黑潮洋流影響下,浮標(biāo)在半徑4 000 m 的范圍內(nèi)漂移。Tadokoro 等[58]基于一個直徑為8 m 的大型錨系浮標(biāo)開發(fā)了GNSS-A 定位系統(tǒng),浮標(biāo)的漂移范圍在150~200 m 內(nèi)。該系統(tǒng)的GNSS 天線和衛(wèi)星通信天線安裝在浮標(biāo)的頂部,聲學(xué)傳感器用一根不銹鋼柱子連接并安裝在水下約1.7 m 處。使用浮標(biāo)觀測平臺進(jìn)行了106 d 測試,GNSS-A 定位系統(tǒng)能夠獲得高質(zhì)量的聲學(xué)數(shù)據(jù)。日本東北大學(xué)在2010 年采用小型浮標(biāo)進(jìn)行了海上試驗,由于電力供應(yīng)問題,試驗僅持續(xù)了2 d 時間[65]。在此基礎(chǔ)上,日本東北大學(xué)改用大型錨系浮標(biāo),并在熊野灘進(jìn)行了2 次實時的GNSS-A 觀測試驗,成功獲取了近10 個月的觀測數(shù)據(jù)。目前,浮標(biāo)觀測平臺仍需解決電氣系統(tǒng)問題,以滿足系統(tǒng)運行的電力需求。

    2.3.2 數(shù)據(jù)采集方法

    GNSS-A 定位系統(tǒng)的測量值為聲學(xué)信號從換能器到應(yīng)答器的往返傳播時間。聲學(xué)信號傳播時間與海水中的聲波傳播速度密切相關(guān),需要利用聲速剖面儀(Sound Velocity Profiler,SVP)、溫鹽深測量儀(Conductivity-Temperature-Depth profiler,CTD)和拋棄式溫鹽深測量儀(Expendable Conductivity-Temperature-Depth profiler,XCTD)等設(shè)備采集現(xiàn)場聲速剖面,每隔幾小時進(jìn)行一次。系統(tǒng)同時收集動態(tài)GNSS 數(shù)據(jù),船只的姿態(tài)由姿態(tài)傳感器測量,GNSS 天線和換能器之間的相對位置在船塢通過地面測量確定,用于確定換能器相對于GNSS 天線的坐標(biāo)。

    1)測量策略

    GNSS-A 定位系統(tǒng)的測量策略按海面觀測平臺與海底大地基準(zhǔn)站的位置關(guān)系,可分為靜態(tài)測量和動態(tài)測量。靜態(tài)測量將海面觀測平臺維持在海底大地基準(zhǔn)網(wǎng)中心上方連續(xù)進(jìn)行聲學(xué)測量,其主要應(yīng)用代表為美國Scripps 海洋研究所和日本東北大學(xué)[17,74]。靜態(tài)測量需要事先確定海底大地基準(zhǔn)站之間的相對位置和深度,只能獲得海底大地基準(zhǔn)網(wǎng)中心的水平位移[4]。在此基礎(chǔ)上,日本海上保安廳、東京大學(xué)以及名古屋大學(xué)開發(fā)了動態(tài)測量策略,通過海面平臺按預(yù)定航跡圍繞海底大地基準(zhǔn)站動態(tài)觀測聲學(xué)數(shù)據(jù),動態(tài)測量策略還可以進(jìn)行垂直方向定位[19]。與GNSS 衛(wèi)星分布對定位影響類似,空間分布良好的聲學(xué)數(shù)據(jù)可以減少估計基準(zhǔn)站的位置偏差。通過航行觀測收集幾何對稱的數(shù)據(jù),GNSS-A 定位精度和觀測效率都得到了提高[60,63]。為了提高定位的穩(wěn)定性,通常假設(shè)在整個觀測期海底大地基準(zhǔn)網(wǎng)的幾何形狀是恒定的[75-76]。Honson 等采用了兼容的測量策略,同時使用靜態(tài)和動態(tài)測量模式采集GNSS-A 數(shù)據(jù)來確定海底大地基準(zhǔn)網(wǎng)中心位移[74,77-78]。對于海底大地基準(zhǔn)網(wǎng)的標(biāo)校,也可采用靜態(tài)和動態(tài)測量策略??紤]標(biāo)校效率和精度,可以采用對單個海底大地基準(zhǔn)站進(jìn)行絕對校準(zhǔn),對海底大地基準(zhǔn)網(wǎng)(站)進(jìn)行相對測量,即絕對校準(zhǔn)+相對測量策略[45-46,79]。前者實現(xiàn)絕對基準(zhǔn)從海面?zhèn)鬟f到海底大地基準(zhǔn)站(時間長且受全水深聲速誤差影響),海底大地基準(zhǔn)網(wǎng)點之間相對測距(海底等溫層,聲速的影響非常小,相對測量精度高),類似于GNSS 控制網(wǎng)測量方法,采用網(wǎng)平差技術(shù)實現(xiàn)整體解算。

    2)聲學(xué)觀測

    GNSS-A 定位的核心是聲學(xué)觀測,需要精確測量聲學(xué)信號在換能器和應(yīng)答器之間的往返傳播時間。海面換能器向海底應(yīng)答器發(fā)送聲學(xué)信號,海底應(yīng)答器作為信號轉(zhuǎn)發(fā)射器,接收并返回從海面平臺發(fā)出的聲學(xué)信號。東京大學(xué)開發(fā)的GNSS-A 定位系統(tǒng)中,換能器發(fā)射2 個雙相調(diào)制的聲學(xué)信號,分別用于信號識別和測距信號,載波頻率為10 kHz;系統(tǒng)中使用的調(diào)制序列碼是8 階、9 階M 序列碼(M-sequence)。名古屋大學(xué)開發(fā)的GNSS-A 定位系統(tǒng)則采用了5 階M 序列碼[58,60]。如果應(yīng)答器識別信號,應(yīng)答器就會記錄下后續(xù)的測距信號,并在設(shè)定的間隔時間后,應(yīng)答器將重新編碼的測距信號與新的識別信號一起發(fā)送海面。在聲學(xué)測距期間,利用M 序列碼的特性,可以通過檢測合成信號和接收的聲學(xué)信號之間的互相關(guān)函數(shù)(Cross-Correlation Function,CCF)峰值,確定精確的信號傳播時間[19,73]。

    GNSS-A 測量的聲學(xué)信號傳播時間需要轉(zhuǎn)換為換能器與應(yīng)答器之間的空間距離??臻g距離計算值由聲線跟蹤確定,需要海水中的聲速剖面,利用SVP 測量數(shù)據(jù),或者利用CTD 測量的溫鹽深數(shù)據(jù)采用經(jīng)驗公式轉(zhuǎn)換為聲速[80]。采用線性或方形方式布設(shè)聲速剖面測量站位,聲速剖面測量采用定點SVP 為主、CTD 為輔的方式,適當(dāng)補充XCTD 測量。調(diào)查區(qū)內(nèi)影響聲速的水文條件(溫度、鹽度等)變化較大時,需增加聲速剖面的測量次數(shù),在XCTD 作業(yè)時,船速必須降至適當(dāng)航速。

    2.3.3 標(biāo)校線路設(shè)計

    空間分布良好的聲學(xué)數(shù)據(jù)可以提高GNSS-A 定位精度。靜態(tài)測量策略中,海面觀測平臺保持在海底大地基準(zhǔn)網(wǎng)中心上方,觀測平臺與基準(zhǔn)站形成一個圓錐形對稱結(jié)構(gòu);動態(tài)測量策略中,海面觀測平臺通常圍繞海底大地基準(zhǔn)站沿圓形航跡航行,觀測平臺與基準(zhǔn)站形成倒圓錐形對稱結(jié)構(gòu),圓形航跡被認(rèn)為是定位精度最高的標(biāo)校線路[28,81]??紤]到船舶姿態(tài)和海水湍流對聲學(xué)測量的影響,海面平臺最大航行速度為6~7 n mile/h[82]。因此,往往需要相對較長的觀測時間,而且很難通過提高船速來減少觀測時間。

    為了提高定位精度和效率,趙建虎等[27]利用幾何精因子(Geometric Dilution Of Positioning,GDOP)表示定位精度,通過尋求GDOP 的最小值,給出圓形航跡最佳半徑。如果修正了聲速誤差,航行半徑為海底大地基準(zhǔn)站深度的1.414 倍時,定位精度最高;若考慮聲速影響,航行半徑為深度的1.045 倍時,定位精度最高。另外,海面觀測平臺與海底大地基準(zhǔn)站形成的幾何結(jié)構(gòu)決定了Fisher 信息量,航行半徑為深度的1.414 倍時,F(xiàn)isher 準(zhǔn)則矩陣表征的定位精度也達(dá)到最佳[26]。

    由于GNSS-A 定位中平面方向的幾何對稱性,平面定位精度有了很大的提高[81],但海面觀測平臺與海底大地基準(zhǔn)站的高度差幾乎是相同的,導(dǎo)致GNSS-A 定位在垂直方向上存在固有的幾何缺陷,需要增加垂直方向的約束,如海底深度傳感器測量的深度[26,79]。另外,增加一個過頂?shù)氖纸徊婧桔E也可以提高垂直方向定位精度[26]。考慮聲速空間變化,同時利用多條聲學(xué)路徑可以實現(xiàn)高精度的聲學(xué)觀測[71]。因此,設(shè)計空間對稱的圓形和十字交叉標(biāo)校線路是合理和科學(xué)的(圖3)。

    圖3 標(biāo)校線路設(shè)計Fig.3 Design of survey line

    3 GNSS-A 數(shù)據(jù)處理技術(shù)

    3.1 誤差來源分析

    GNSS-A 定位誤差源包括2 類:一類是海面換能器有關(guān)的誤差,如GNSS 定位誤差、GNSS 天線與海面換能器相對位置偏差、姿態(tài)偏差等;另一類是聲學(xué)信號傳播有關(guān)的誤差。聲學(xué)換能器位置通過GNSS 定位獲取,如Hexagon 公司的VeriPos、NavCom 公司的Starfire、Fugro 集團(tuán)的OmniSTAR 等實時定位服務(wù),以及后處理精密單點定位(Precise Point Positioning,PPP),前者可提供分米量級的實時導(dǎo)航定位服務(wù),后者則能夠?qū)崿F(xiàn)厘米級的動態(tài)定位精度[83-84];慣性測量單元的姿態(tài)測量精度可達(dá)0.01°[72];GNSS 天線與海面換能器相對位置可通過全站儀精確測定,也可作為待估參數(shù)參與定位解算[25]。聲學(xué)信號傳播有關(guān)的誤差有2 種:一種是傳播時間的量測誤差,一種是聲速誤差。

    GNSS-A 定位系統(tǒng)可以精確測量直達(dá)波聲學(xué)信號傳播時間,時間量測誤差約為3~10 μs,海水中聲學(xué)信號若以1 500 m/s 的速度傳播,聲學(xué)測距的分辨率優(yōu)于1.5 cm[43,59,80]。直達(dá)波聲學(xué)信號經(jīng)常受到海面反射波的影響,在互相關(guān)波形時序中,波峰與直達(dá)信號到達(dá)時刻存在偏差[64]。反射波表現(xiàn)為多重峰值,若將反射波錯誤識別為直達(dá)波,會導(dǎo)致聲學(xué)信號傳播時間的量測誤差[58,85]。Honsho 等[86]考慮了聲學(xué)信號入射角與互相關(guān)波形的關(guān)系,提出了一種相位互相關(guān)法(Phase-Only Correlation,POC),以減少傳播時間測量的不確定性;Tadokoro 等[58]引入了能量比概念(Energy Ratio,ER)檢測識別直接波。由于Tadokoro 等設(shè)計了浮標(biāo)平臺搭載GNSS-A 定位系統(tǒng),直接波和海面反射波之間存在1.3~1.8 ms 的時間延遲,即使聲學(xué)信號被海面反射污染,也可以有效地分辨出反射波。此外,動態(tài)測量策略中海面觀測平臺與海底大地基準(zhǔn)站的相對運動會引起聲學(xué)信號的多普勒頻移,通常在聲學(xué)信號的探測與識別過程中需移除多普勒效應(yīng)影響[63,87]。

    聲速誤差是GNSS-A 定位的主要誤差,來源于沿聲學(xué)信號傳播路徑的海洋聲速結(jié)構(gòu)(Sound Speed Structure,SSS)的時空變化。海洋聲速主要與海水溫度、鹽度和靜壓力有關(guān)[88]。海表面混合層的存在改變了層內(nèi)的溫度和鹽度分布,而溫度和鹽度主要影響上層海洋的聲速值,隨著深度增加,海水混合作用減弱,溫度和鹽度分布趨于穩(wěn)定,壓力隨著深度的增加而增加,聲速也呈單向增加趨勢。因此,從海表面到海底的SSS 主要呈垂向分層[89],SSS 的時空變化主要體現(xiàn)在海洋上層。此外,SSS 具有混合的時空異質(zhì)性,聲速的時間、空間變化尺度均有所不同[61,90]。相比聲學(xué)信號傳播時間,海洋聲速變化具有更大的空間和時間尺度,難以實現(xiàn)聲速的連續(xù)測量以表征聲速時空變化,而且聲速測量準(zhǔn)確度通常是相對的,由于聲速測量設(shè)備的校準(zhǔn)誤差,可能每秒會出現(xiàn)幾米的偏差,一個恒定的聲速偏差也可作為未知參數(shù)參與定位解算[89]?;趯崪y聲速數(shù)據(jù),使用經(jīng)驗正交函數(shù)反演的聲速存在時空分辨率不足等問題[23]。聲學(xué)信號傳播時間包含了聲學(xué)信號路徑的聲速變化信息[89],因此,在GNSS-A 定位中,可以同時估計聲速的變化和海底大地基準(zhǔn)站的位置[31,77,89]。

    3.2 參數(shù)估計方法

    結(jié)合GNSS 和聲學(xué)的觀測數(shù)據(jù)可以估計海底大地基準(zhǔn)站的位置。理想情況下,如果海底大地基準(zhǔn)站位置的估計值與真實值一致,那么計算出的兩個單程傳播時間之和應(yīng)該等于觀測的往返傳播時間;當(dāng)海底大地基準(zhǔn)站位置的估計值與真實值有偏差時,計算的傳播時間與觀測時間也存在偏差。海底大地基準(zhǔn)站位置最優(yōu)估計的過程,是使傳播時間測量值和計算之間的時間偏差(Observation Minus Computed,OMC)平方和最小的過程[31,87,91],即

    3.2.1 解算策略

    GNSS-A 定位用于獲取海底大地基準(zhǔn)站的位置,其參數(shù)估計策略分為整體解和歷元解兩類。整體解算方法通過長期觀測的GNSS-A 測量數(shù)據(jù)來完成定位解算[31,76-77],24 h 的觀測數(shù)據(jù)可獲得±1.5 cm 的定位精度,3~4 d 的觀測數(shù)據(jù)可獲得亞厘米級定位精度[5]。歷元解算則利用1 組聲學(xué)觀測數(shù)據(jù)完成定位解算[61],由于海洋聲速結(jié)構(gòu)的時空變化,歷元解算存在幾十厘米的定位誤差[91-92]。通常利用靜態(tài)測量策略的觀測數(shù)據(jù)進(jìn)行定位解算,即在海底大地基準(zhǔn)網(wǎng)的中心上方采集GNSS-A 數(shù)據(jù),其定位精度隨距離海底大地基準(zhǔn)網(wǎng)中心越遠(yuǎn)會逐漸下降[64,67,71]。基于海洋聲速結(jié)構(gòu)具有垂向分層特征假設(shè),利用對稱分布的長時間聲學(xué)觀測數(shù)據(jù),可以削弱聲速時域變化對海底大地基準(zhǔn)站位置估計的影響,提高海底大地基準(zhǔn)網(wǎng)中心的水平方向位置精度,但是靜態(tài)測量策略無法獲取精確的垂直方向位置[31,87]。因此,靜態(tài)測量策略需要額外的觀測,如壓力傳感器觀測的深度信息,以增加垂直方向的約束。將聲學(xué)信號到達(dá)方向和聲學(xué)信號傳播時間觀測數(shù)據(jù)同時納入最優(yōu)化反演,可以提高定位解算精度[72]。動態(tài)測量策略按預(yù)定航跡動態(tài)采集GNSS-A 數(shù)據(jù),可以獲得海底大地基準(zhǔn)站的水平和垂直方向位置。部分學(xué)者通過GNSS-A 動態(tài)測量策略的觀測數(shù)據(jù)估計單個海底大地基準(zhǔn)站的位置[32-33,77,87],進(jìn)而確定海底大地基準(zhǔn)網(wǎng)幾何形狀。在假設(shè)海底大地基準(zhǔn)網(wǎng)幾何形狀是恒定的情況下,約束基準(zhǔn)網(wǎng)的幾何形狀,然后根據(jù)靜態(tài)測量策略的觀測數(shù)據(jù),可精確估計基準(zhǔn)網(wǎng)的中心位置[4,75,77,91]。對于海底大地基準(zhǔn)網(wǎng)的絕對校準(zhǔn)+相對測量策略,可以采用網(wǎng)平差技術(shù)實現(xiàn)整體解算。受海底大地基準(zhǔn)網(wǎng)的網(wǎng)形結(jié)構(gòu)影響,海底大地基準(zhǔn)站的垂線解精度不高,甚至出現(xiàn)解算不穩(wěn)定問題,因此同樣需要引入外部約束,如壓力傳感器測量的深度/深度差[45-46]。

    3.2.2 聲速誤差估計

    參數(shù)估計策略受到海洋聲速結(jié)構(gòu)時域變化的影響。在聲速結(jié)構(gòu)具有垂向分層特征的基礎(chǔ)上,F(xiàn)ujita 等[31]、Ikuta 等[87]提出了同時估計海底大地基準(zhǔn)站的位置和聲速時域變化的方法,前者利用定位殘差提取多項式擬合的聲速時域變化,迭代修正位置參數(shù),后者則基于B 樣條函數(shù)參數(shù)化表示平均聲速時域變化,聯(lián)合位置參數(shù)迭代計算;進(jìn)一步地,GNSS 定位中“天頂對流層延遲”(Zenith Total Delay,ZTD)[93]的概念被引入,Kido 等[89]、Honsho 等[77]提出了用“垂向總延遲”(Nadir Total Delay,NTD)表示平均聲速時域變化,聯(lián)合解算NTD 和位置參數(shù)的方法。

    實際海洋聲速結(jié)構(gòu)同時存在時間和空間變化。Yasuda 等[94]假設(shè)1 000 m 以淺水層的聲速結(jié)構(gòu)受黑潮的影響向一個方向傾斜,構(gòu)建了包含聲速時域及其水平梯度變化的模型;Yokota 等[95]從直接估計的聲速時域變化中提取了淺水層的聲速水平梯度變化,但與海底大地基準(zhǔn)站位置有關(guān)的深水層梯度變化無法用聲速時變模型表示,結(jié)合Fujita 等[31]的方法從定位殘差中可以提取這部分聲速變化;Honsho 等[78]在NTD 的基礎(chǔ)上,考慮了更為普遍的方向性聲速梯度(Direction NTD Gradient),主要與海底大地基準(zhǔn)站位置有關(guān)。Tomita 等[61]、Watanabe 等[96]分別測試了靜態(tài)、動態(tài)測量策略的海底大地基準(zhǔn)站位置、聲速時域變化和聲速水平梯度參數(shù)聯(lián)合估計方法。

    海洋聲速結(jié)構(gòu)還具有混合的時空異質(zhì)性,包括短周期(短波長)的異質(zhì)性以及長周期(長波長)的異質(zhì)性。通過空間對稱分布的長時間GNSS-A 數(shù)據(jù),整體解算策略可以成功消除聲速結(jié)構(gòu)的長周期異質(zhì)性[78]。而受到聲速結(jié)構(gòu)的短周期異質(zhì)性影響,歷元解算策略的定位結(jié)果較差[74]。劉楊等[32-33]指出定位的聲速殘余誤差主要來自海流、內(nèi)波等引起的聲速水平梯度變化,通過聯(lián)合估計聲速垂向結(jié)構(gòu)和聲速水平梯度的時域變化,降低了聲速殘余誤差,提高了海底大地基準(zhǔn)站坐標(biāo)精度。目前仍缺乏深入的海流、內(nèi)波引起聲速變化物理機制探究,需要進(jìn)一步了解其時空特征(如時間、空間范圍)[71,90]。除了更密集頻繁的船載GNSS-A 定位,結(jié)合多個觀測平臺的協(xié)同觀測方式,可以更加有效地建立短周期異質(zhì)性模型[60-61,90]。

    3.3 精度評定方法

    GNSS-A 定位結(jié)果精度評價通過比較不同解算結(jié)果的可重復(fù)性,如考慮不同策略觀測之間的重復(fù)性[63,97]。Sato 等[63]對比了不同觀測策略的定位結(jié)果,動態(tài)觀測策略的重復(fù)性約為2 cm,比靜態(tài)觀測策略平均降低了30%。GNSS-A 定位精度也可以通過對斜距殘差的統(tǒng)計分析間接描述[98]。而Chen 等[59]設(shè)計了1 套海底聲學(xué)應(yīng)答器系統(tǒng)(Seafloor Acoustic Transponder System,SATS)來直接評估GNSS-A 定位的準(zhǔn)確性。SATS 有3 個應(yīng)答器和1 個姿態(tài)傳感器,可以提供應(yīng)答器基線的真實長度和真實姿態(tài)信息,以確保對GNSS-A 定位評估的可靠性和有效性。使用GNSS-A 定位數(shù)據(jù)估計SATS 上3 個應(yīng)答器的位置,并計算基線長度和姿態(tài)角,通過比較真實基線和姿態(tài)與計算結(jié)果,可以直接評估GNSS-A 海底定位的準(zhǔn)確性。Chen 等[59]還進(jìn)行了敏感性分析,以研究GNSS-A 定位結(jié)果對聲速變化的穩(wěn)健性。在聲速偏差很大的情況下會嚴(yán)重惡化GNSS-A 定位的質(zhì)量,敏感性分析同樣可以證實GNSS 和聲學(xué)測量是可靠的。

    4 應(yīng)用及展望

    4.1 海底大地基準(zhǔn)建設(shè)技術(shù)主要應(yīng)用

    海底大地基準(zhǔn)建設(shè)技術(shù)除滿足水下潛航器導(dǎo)航定位需求外,目前主要應(yīng)用于海底精密工程測量、海底形變監(jiān)測。近幾十年來,陸地大地測量手段,如合成孔徑雷達(dá)(InSAR)和全球定位系統(tǒng)(GNSS),提供了近地空間、海洋表面高時空分辨率的定位和監(jiān)測能力[5,99-100]。然而,InSAR 或GNSS 無法適用于水下、海底場景。采用GNSS 和聲學(xué)測距的海底定位技術(shù)可以應(yīng)用于海底精密工程測量,包括大型海洋平臺安裝和穩(wěn)定性監(jiān)測、海底大型沉管隧道對接等。監(jiān)測海底地殼形變對于掌握全球地質(zhì)構(gòu)造運動至關(guān)重要,自20 世紀(jì)80 年代提出聲學(xué)方法監(jiān)測海底形變的概念以來[4],基于GNSS-A 定位技術(shù)已實現(xiàn)了在全球參考框架內(nèi)對海底點的精確測量[5],海底水平構(gòu)造運動監(jiān)測的分辨率達(dá)到了厘米級[17],并成功應(yīng)用于揭示大洋構(gòu)造板塊的運動和變形、俯沖帶和其他板塊的地震過程,以及海底火山和擴張中心的變形[101-102]。

    基于海底大地基準(zhǔn)測量技術(shù)的海底形變監(jiān)測應(yīng)用在日本取得了較好的發(fā)展[8]。日本水文和海洋局(Hydrographic and Oceanographic Department,Japan Coast Guard)已經(jīng)在日本太平洋一側(cè)的海底部署了15 個以上的海底大地基準(zhǔn)站,用以檢測和監(jiān)測由大洋板塊俯沖引起的海底地殼形變[31]。2011 年日本東北9.0 級(Mw)地震后,日本開發(fā)了一種新型的海底大地基準(zhǔn)站,可以在5 000 m 以深的海底工作,其聲學(xué)通信范圍大于15 km[65]。2012 年沿日本海溝建造了20 個新的GNSS-A 海底大地基準(zhǔn)站,大部分站點位于海溝軸線兩側(cè),深度大于5 000 m。

    4.2 海底大地基準(zhǔn)網(wǎng)建設(shè)展望

    海底大地基準(zhǔn)網(wǎng)是新型海洋觀測平臺。雖然以GNSS-A 為代表的海底大地測量技術(shù)已取得了顯著的進(jìn)步,海底大地基準(zhǔn)網(wǎng)的建設(shè)工作仍然任重道遠(yuǎn),面對新一代國家綜合PNT 系統(tǒng)的建設(shè)需求和水下導(dǎo)航定位實時應(yīng)用需求,需要進(jìn)一步開展海底大地測量技術(shù)的海洋聲速誤差處理、時間同步及標(biāo)校方法等方面的技術(shù)研究,也需要在分級組網(wǎng)和立體觀測方面開展應(yīng)用研究。

    4.2.1 技術(shù)層面

    1)聲速誤差處理

    GNSS-A 定位精度受聲速時空變化影響,需要精細(xì)化處理聲速。通過海洋環(huán)境信息構(gòu)建區(qū)域/全球海洋聲速場模型,采用同步觀測方式反演聲速誤差、生成并播發(fā)聲速誤差改正產(chǎn)品,提高水下導(dǎo)航定位精度。

    2)時間同步

    海底大地基準(zhǔn)網(wǎng)的空間基準(zhǔn)可統(tǒng)一至國際橢球參考框架,但是在時間基準(zhǔn)統(tǒng)一方面差距巨大,需要重點解決海底大地基站間的時間同步以及水下用戶與基準(zhǔn)網(wǎng)之間的時間同步,海底大地基準(zhǔn)網(wǎng)的時間基準(zhǔn)統(tǒng)一方法將是建設(shè)海洋綜合PNT 的難點和關(guān)鍵。

    3)標(biāo)校方式

    調(diào)查船沿預(yù)定航跡航行標(biāo)校海底大地基準(zhǔn)站位置觀測時間長、觀測頻次低,且耗費大量人力物力,需要進(jìn)一步優(yōu)化GNSS-A 定位方法。一方面可發(fā)展多種觀測手段,如聯(lián)合船舶、浮標(biāo)、無人艇、水下滑翔器等平臺,優(yōu)化觀測航跡及觀測采樣頻率,實現(xiàn)快速實時標(biāo)校;另一方面可借助立體組網(wǎng)優(yōu)勢,綜合多源信息,提高GNSS-A 位置標(biāo)校的效率、精度,增加海底大地基準(zhǔn)測量技術(shù)的可用性、可靠性、精確性。

    4.2.2 應(yīng)用層面

    1)分級布設(shè)

    考慮到我國海底大地基準(zhǔn)網(wǎng)建設(shè)的需要和巨大的經(jīng)費需求,需進(jìn)一步探索海底大地基準(zhǔn)網(wǎng)的分級布設(shè)方法。由于海底大地基準(zhǔn)網(wǎng)使用率低且不易長期保存,在布設(shè)海底大地基準(zhǔn)站時,應(yīng)采用“分級布設(shè)”原則,先在大范圍內(nèi)布設(shè)首級網(wǎng),然后在建成區(qū)域或急需提供水下導(dǎo)航定位服務(wù)區(qū)域布設(shè)次級網(wǎng)。這樣既能達(dá)到合理布網(wǎng)、滿足急需急用的目的,又能方便大地基準(zhǔn)網(wǎng)的擴展和加密??紤]不同用戶的應(yīng)用需求,在布設(shè)海底大地基準(zhǔn)站時,還應(yīng)采用“低頻引導(dǎo),高頻定位”原則,通過在海底大地基站搭載低頻信號聲吶設(shè)備提供遠(yuǎn)距離位置引導(dǎo)服務(wù),布設(shè)高頻信號聲吶設(shè)備提供近距離高精度導(dǎo)航定位服務(wù)。

    2)立體組網(wǎng)

    面向大范圍、高精度、多領(lǐng)域的海洋觀測需求,需進(jìn)一步發(fā)展海底大地基準(zhǔn)組網(wǎng)和服務(wù)方式,也需要綜合海底觀測網(wǎng)和海底大地基準(zhǔn)網(wǎng)的統(tǒng)一建設(shè)問題。通過構(gòu)建海面、水下、海底一體化立體觀測網(wǎng)絡(luò),進(jìn)一步優(yōu)化網(wǎng)型設(shè)計、布放回收、位置標(biāo)校等方法,擴大服務(wù)范圍、提高定位精度;同時,支持多傳感器協(xié)同作業(yè),實現(xiàn)大范圍、多領(lǐng)域的海洋觀測。

    猜你喜歡
    測量
    測量重量,測量長度……
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計算
    滑動摩擦力的測量與計算
    測量的樂趣
    二十四節(jié)氣簡易測量
    日出日落的觀察與測量
    滑動摩擦力的測量與計算
    測量
    測量水的多少……
    观看免费一级毛片| 精品久久久久久久久久久久久| 欧美丝袜亚洲另类 | 一个人免费在线观看电影| 精品一区二区三区视频在线| 淫秽高清视频在线观看| 老司机深夜福利视频在线观看| 乱码一卡2卡4卡精品| 精品99又大又爽又粗少妇毛片 | 欧美日韩福利视频一区二区| 成年人黄色毛片网站| 日韩有码中文字幕| 老司机午夜十八禁免费视频| 非洲黑人性xxxx精品又粗又长| 很黄的视频免费| 国产精品久久久久久久久免 | 欧美激情久久久久久爽电影| 精品不卡国产一区二区三区| 夜夜看夜夜爽夜夜摸| 久久欧美精品欧美久久欧美| 亚洲成人中文字幕在线播放| 日韩欧美三级三区| 最近视频中文字幕2019在线8| 69人妻影院| 男人的好看免费观看在线视频| 国内揄拍国产精品人妻在线| 免费人成在线观看视频色| 99久久精品一区二区三区| 午夜视频国产福利| 欧美日韩国产亚洲二区| 国产伦在线观看视频一区| 国产探花极品一区二区| 亚洲av二区三区四区| 国产精品亚洲美女久久久| 成人美女网站在线观看视频| eeuss影院久久| 国产伦精品一区二区三区视频9| 久久久成人免费电影| 97超视频在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 久久久精品欧美日韩精品| 乱人视频在线观看| 欧美3d第一页| 搞女人的毛片| 久久欧美精品欧美久久欧美| 天美传媒精品一区二区| 亚洲男人的天堂狠狠| 成人高潮视频无遮挡免费网站| 亚洲乱码一区二区免费版| 一个人免费在线观看的高清视频| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 男女做爰动态图高潮gif福利片| 啪啪无遮挡十八禁网站| 尤物成人国产欧美一区二区三区| 九九在线视频观看精品| 亚洲av免费在线观看| 久久九九热精品免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | www.色视频.com| 99久久无色码亚洲精品果冻| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| 哪里可以看免费的av片| 午夜精品久久久久久毛片777| 天堂网av新在线| 日本黄大片高清| 一二三四社区在线视频社区8| 日韩亚洲欧美综合| 精品久久久久久久人妻蜜臀av| 国产乱人视频| 乱人视频在线观看| 久久伊人香网站| 久久国产精品影院| 国产国拍精品亚洲av在线观看| 真人一进一出gif抽搐免费| 天堂网av新在线| 两个人的视频大全免费| 欧美性猛交黑人性爽| 日韩有码中文字幕| 国产精品99久久久久久久久| 丰满的人妻完整版| 在现免费观看毛片| 成人精品一区二区免费| 国产三级在线视频| 亚洲五月天丁香| 高清在线国产一区| 久久精品91蜜桃| 亚洲精品色激情综合| 真人一进一出gif抽搐免费| 窝窝影院91人妻| 亚洲精品亚洲一区二区| 男女之事视频高清在线观看| 精品久久久久久久久久久久久| 夜夜看夜夜爽夜夜摸| 国产欧美日韩精品亚洲av| 久9热在线精品视频| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| av在线天堂中文字幕| 国产精品伦人一区二区| 久久香蕉精品热| 97超视频在线观看视频| 久久99热6这里只有精品| 欧美激情国产日韩精品一区| 免费av毛片视频| 日本一二三区视频观看| 亚洲成av人片在线播放无| 国产三级在线视频| 国产成人影院久久av| 最近最新中文字幕大全电影3| 动漫黄色视频在线观看| 国产黄a三级三级三级人| x7x7x7水蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 国产91精品成人一区二区三区| 亚洲欧美日韩高清专用| 真人做人爱边吃奶动态| 午夜福利欧美成人| 高清毛片免费观看视频网站| 亚洲午夜理论影院| 男人舔女人下体高潮全视频| 免费搜索国产男女视频| 天天躁日日操中文字幕| 99视频精品全部免费 在线| 久久精品91蜜桃| 91午夜精品亚洲一区二区三区 | av天堂在线播放| 久久久久久久久大av| 免费大片18禁| 99在线视频只有这里精品首页| 成人亚洲精品av一区二区| 午夜免费男女啪啪视频观看 | 亚洲人与动物交配视频| 亚洲欧美日韩高清在线视频| 国产精品98久久久久久宅男小说| 亚洲国产欧洲综合997久久,| 国产成人欧美在线观看| 亚洲中文字幕一区二区三区有码在线看| 有码 亚洲区| 国产69精品久久久久777片| 亚洲va日本ⅴa欧美va伊人久久| 久久婷婷人人爽人人干人人爱| 久久午夜亚洲精品久久| 日韩国内少妇激情av| 国产伦精品一区二区三区视频9| 午夜精品在线福利| 免费看光身美女| 我的女老师完整版在线观看| 免费在线观看成人毛片| 变态另类成人亚洲欧美熟女| 亚洲最大成人av| 老司机福利观看| 亚洲电影在线观看av| 国产亚洲欧美在线一区二区| 久久草成人影院| 小说图片视频综合网站| 国产精品自产拍在线观看55亚洲| 免费黄网站久久成人精品 | 日本撒尿小便嘘嘘汇集6| 欧美日韩瑟瑟在线播放| 亚洲熟妇中文字幕五十中出| 成人亚洲精品av一区二区| 成人永久免费在线观看视频| 欧美日韩国产亚洲二区| 国产在视频线在精品| 亚洲欧美日韩东京热| 国产午夜精品久久久久久一区二区三区 | 日本黄大片高清| 国产免费av片在线观看野外av| 精品免费久久久久久久清纯| 99热6这里只有精品| 国产av麻豆久久久久久久| 久久精品91蜜桃| 免费电影在线观看免费观看| 欧美激情在线99| 蜜桃亚洲精品一区二区三区| 亚洲午夜理论影院| 久久精品人妻少妇| 精品久久国产蜜桃| 麻豆国产av国片精品| av在线老鸭窝| 村上凉子中文字幕在线| 性色avwww在线观看| 赤兔流量卡办理| 一a级毛片在线观看| 久久久色成人| 老熟妇乱子伦视频在线观看| 国产一区二区三区在线臀色熟女| 国产在线精品亚洲第一网站| 最好的美女福利视频网| 麻豆成人午夜福利视频| 又黄又爽又免费观看的视频| 国产黄a三级三级三级人| 特大巨黑吊av在线直播| 动漫黄色视频在线观看| 天美传媒精品一区二区| 免费一级毛片在线播放高清视频| 国产免费av片在线观看野外av| 国产成人aa在线观看| 伦理电影大哥的女人| 十八禁网站免费在线| 又爽又黄无遮挡网站| 亚洲第一区二区三区不卡| 超碰av人人做人人爽久久| 精品人妻一区二区三区麻豆 | netflix在线观看网站| 国产免费av片在线观看野外av| 亚洲人与动物交配视频| 永久网站在线| 亚洲五月天丁香| 色综合婷婷激情| 婷婷六月久久综合丁香| 一级黄片播放器| 欧美不卡视频在线免费观看| 亚洲av电影在线进入| 直男gayav资源| 嫩草影视91久久| 欧美xxxx黑人xx丫x性爽| 日韩欧美在线乱码| 国产高清视频在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 精品一区二区三区人妻视频| 亚洲精品影视一区二区三区av| 人人妻人人看人人澡| 久久99热这里只有精品18| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜| 亚洲av第一区精品v没综合| 天堂网av新在线| 国产一区二区在线观看日韩| 97超级碰碰碰精品色视频在线观看| 久久国产精品影院| 久久久久久久精品吃奶| 欧美乱色亚洲激情| 成人特级av手机在线观看| 久久伊人香网站| 国产男靠女视频免费网站| 久久久久九九精品影院| 婷婷精品国产亚洲av| 国产一区二区三区在线臀色熟女| 99精品在免费线老司机午夜| 亚洲精品色激情综合| 如何舔出高潮| 亚洲精品一区av在线观看| 免费观看的影片在线观看| 天天一区二区日本电影三级| av在线天堂中文字幕| 精品一区二区三区av网在线观看| 他把我摸到了高潮在线观看| 99热这里只有精品一区| 精华霜和精华液先用哪个| 18禁黄网站禁片午夜丰满| 欧美最新免费一区二区三区 | 日日夜夜操网爽| www.色视频.com| 美女高潮的动态| 国产午夜福利久久久久久| 欧美成人性av电影在线观看| 欧美+日韩+精品| 少妇的逼好多水| 亚洲不卡免费看| 色哟哟哟哟哟哟| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 成人美女网站在线观看视频| 欧美区成人在线视频| 美女高潮喷水抽搐中文字幕| 午夜精品久久久久久毛片777| 日韩欧美在线乱码| 精华霜和精华液先用哪个| 最后的刺客免费高清国语| 亚洲在线自拍视频| 亚洲专区中文字幕在线| av在线天堂中文字幕| 成人亚洲精品av一区二区| 在线a可以看的网站| 欧美黄色片欧美黄色片| 亚洲成人久久爱视频| 国产精品99久久久久久久久| 黄色女人牲交| 亚洲人成电影免费在线| 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 91麻豆av在线| 亚洲自拍偷在线| 日韩 亚洲 欧美在线| 久久精品91蜜桃| 啪啪无遮挡十八禁网站| 少妇高潮的动态图| 国产精品女同一区二区软件 | 国产高潮美女av| 久久午夜福利片| 搞女人的毛片| 亚洲av一区综合| 国产精品野战在线观看| 深夜精品福利| 日本一本二区三区精品| 欧美黄色片欧美黄色片| 男女那种视频在线观看| 亚洲av.av天堂| 欧美极品一区二区三区四区| 91av网一区二区| 午夜激情欧美在线| 日韩人妻高清精品专区| 美女cb高潮喷水在线观看| 亚洲av二区三区四区| а√天堂www在线а√下载| 免费av观看视频| 国产国拍精品亚洲av在线观看| 最近最新中文字幕大全电影3| 精品国产三级普通话版| 精品久久久久久成人av| 久久亚洲精品不卡| 日本在线视频免费播放| 久久婷婷人人爽人人干人人爱| 老司机福利观看| 亚洲自拍偷在线| 很黄的视频免费| 999久久久精品免费观看国产| 欧美在线一区亚洲| 亚洲av五月六月丁香网| 91久久精品国产一区二区成人| 韩国av一区二区三区四区| 国产一区二区亚洲精品在线观看| 精品国内亚洲2022精品成人| 床上黄色一级片| 国产色婷婷99| 免费高清视频大片| 免费观看精品视频网站| 国产成人福利小说| 免费人成在线观看视频色| 深夜精品福利| 精品国内亚洲2022精品成人| 色哟哟·www| 怎么达到女性高潮| АⅤ资源中文在线天堂| 成人一区二区视频在线观看| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 琪琪午夜伦伦电影理论片6080| 成人一区二区视频在线观看| 久久久久久久午夜电影| 久久久久久久久中文| 精品久久久久久久久av| 国产欧美日韩精品一区二区| 欧美性猛交黑人性爽| 看片在线看免费视频| 长腿黑丝高跟| 欧美潮喷喷水| 欧美日韩中文字幕国产精品一区二区三区| 看黄色毛片网站| 亚洲精品粉嫩美女一区| 看片在线看免费视频| 精品午夜福利在线看| 亚洲av不卡在线观看| 午夜福利成人在线免费观看| 国产aⅴ精品一区二区三区波| 国产精品亚洲一级av第二区| 偷拍熟女少妇极品色| 日本三级黄在线观看| 在线十欧美十亚洲十日本专区| 哪里可以看免费的av片| 亚洲成a人片在线一区二区| 欧美日韩黄片免| 免费搜索国产男女视频| 久久久久免费精品人妻一区二区| 亚洲av成人av| 国产精品久久电影中文字幕| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 亚洲成av人片在线播放无| 亚洲自偷自拍三级| 成人亚洲精品av一区二区| 国产免费一级a男人的天堂| 国产精品三级大全| 少妇人妻精品综合一区二区 | 内地一区二区视频在线| 欧美一区二区精品小视频在线| 成人午夜高清在线视频| 99久久精品一区二区三区| 国产伦精品一区二区三区视频9| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 欧美一区二区亚洲| 男女那种视频在线观看| 亚洲精品成人久久久久久| 久久九九热精品免费| 少妇熟女aⅴ在线视频| 国产黄a三级三级三级人| 国产av一区在线观看免费| 人人妻人人澡欧美一区二区| 精品熟女少妇八av免费久了| av天堂在线播放| 欧美黄色淫秽网站| 国产一区二区在线观看日韩| 99久久成人亚洲精品观看| 色吧在线观看| 高潮久久久久久久久久久不卡| www.999成人在线观看| 午夜精品一区二区三区免费看| 久久精品影院6| 99久久精品一区二区三区| a级毛片a级免费在线| 91在线观看av| 女人十人毛片免费观看3o分钟| 欧美不卡视频在线免费观看| a在线观看视频网站| 日本黄色视频三级网站网址| 成人av在线播放网站| 日韩欧美国产一区二区入口| 亚洲精品456在线播放app | 69av精品久久久久久| 免费无遮挡裸体视频| 欧美xxxx黑人xx丫x性爽| 日本黄色视频三级网站网址| 国产美女午夜福利| 日韩欧美精品v在线| 91麻豆av在线| 国产91精品成人一区二区三区| 国产成人a区在线观看| 波多野结衣高清作品| 美女被艹到高潮喷水动态| or卡值多少钱| 黄色一级大片看看| 桃红色精品国产亚洲av| 日韩人妻高清精品专区| 国产精品伦人一区二区| 免费av毛片视频| 成人美女网站在线观看视频| 90打野战视频偷拍视频| 亚洲va日本ⅴa欧美va伊人久久| 99久久无色码亚洲精品果冻| 国产老妇女一区| 精品欧美国产一区二区三| 欧美一区二区精品小视频在线| 很黄的视频免费| 丁香欧美五月| 国产aⅴ精品一区二区三区波| 免费av毛片视频| 欧美另类亚洲清纯唯美| 精品人妻1区二区| 亚洲国产精品sss在线观看| 免费大片18禁| 乱码一卡2卡4卡精品| 成人毛片a级毛片在线播放| 最后的刺客免费高清国语| 欧美激情久久久久久爽电影| 极品教师在线免费播放| 热99re8久久精品国产| 国产爱豆传媒在线观看| 一级av片app| 久久婷婷人人爽人人干人人爱| 国产在线精品亚洲第一网站| 国产一区二区三区在线臀色熟女| 99久久九九国产精品国产免费| 桃红色精品国产亚洲av| 此物有八面人人有两片| 国产精品久久久久久久久免 | 国产亚洲精品久久久com| 99久久久亚洲精品蜜臀av| 国产中年淑女户外野战色| www日本黄色视频网| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲欧美精品综合久久99| 国产欧美日韩一区二区三| 亚洲欧美日韩高清在线视频| 一级毛片久久久久久久久女| 婷婷亚洲欧美| 午夜福利视频1000在线观看| 国内少妇人妻偷人精品xxx网站| bbb黄色大片| 亚洲欧美日韩无卡精品| 亚洲不卡免费看| 免费电影在线观看免费观看| 成人一区二区视频在线观看| 国产精品自产拍在线观看55亚洲| 国产精品乱码一区二三区的特点| 久久午夜福利片| 国产极品精品免费视频能看的| 麻豆成人午夜福利视频| 色在线成人网| 级片在线观看| 床上黄色一级片| 国产午夜精品论理片| 啪啪无遮挡十八禁网站| 亚洲第一区二区三区不卡| 国产淫片久久久久久久久 | 一夜夜www| 中文资源天堂在线| 露出奶头的视频| 99热6这里只有精品| 村上凉子中文字幕在线| 亚洲不卡免费看| 深夜精品福利| 免费av观看视频| 69人妻影院| 小说图片视频综合网站| 国产真实乱freesex| 精品人妻1区二区| 精品午夜福利视频在线观看一区| 久久久久久久久大av| 国产探花在线观看一区二区| 99国产综合亚洲精品| 伦理电影大哥的女人| 桃红色精品国产亚洲av| 国内少妇人妻偷人精品xxx网站| 两人在一起打扑克的视频| 婷婷六月久久综合丁香| 丁香六月欧美| 制服丝袜大香蕉在线| 精品福利观看| 午夜精品在线福利| 精品人妻1区二区| 欧美极品一区二区三区四区| 日韩欧美三级三区| 91在线观看av| 一本久久中文字幕| 午夜福利18| 久久国产精品人妻蜜桃| 日韩欧美免费精品| 在线天堂最新版资源| 99久久精品国产亚洲精品| 丰满的人妻完整版| 国产亚洲精品av在线| 在线观看午夜福利视频| 女同久久另类99精品国产91| 最后的刺客免费高清国语| 一级黄片播放器| 动漫黄色视频在线观看| 桃红色精品国产亚洲av| 国产极品精品免费视频能看的| 别揉我奶头~嗯~啊~动态视频| 又粗又爽又猛毛片免费看| 能在线免费观看的黄片| 国内少妇人妻偷人精品xxx网站| 亚洲国产色片| 日韩欧美 国产精品| 成人av在线播放网站| 嫩草影视91久久| 两人在一起打扑克的视频| 日韩成人在线观看一区二区三区| 极品教师在线视频| 一a级毛片在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久大av| 99热这里只有是精品50| 亚洲熟妇熟女久久| 日韩欧美在线乱码| 亚洲av一区综合| 少妇高潮的动态图| 欧美成人一区二区免费高清观看| 日韩中文字幕欧美一区二区| 99在线视频只有这里精品首页| 好男人电影高清在线观看| 日本撒尿小便嘘嘘汇集6| 久久国产乱子免费精品| 免费黄网站久久成人精品 | 久久久久性生活片| av国产免费在线观看| 精品久久久久久久久久久久久| 精品乱码久久久久久99久播| 九色国产91popny在线| 国产欧美日韩精品亚洲av| 亚洲在线自拍视频| av福利片在线观看| 桃红色精品国产亚洲av| 中文字幕av成人在线电影| 精品午夜福利视频在线观看一区| 国产精品,欧美在线| 99久久精品一区二区三区| 制服丝袜大香蕉在线| 99久久九九国产精品国产免费| 午夜精品在线福利| 欧美+亚洲+日韩+国产| 精品午夜福利视频在线观看一区| 亚洲av.av天堂| 真人做人爱边吃奶动态| 偷拍熟女少妇极品色| 精品福利观看| 亚洲成av人片免费观看| 又粗又爽又猛毛片免费看| 国产精品伦人一区二区| 亚洲激情在线av| 午夜免费男女啪啪视频观看 | 国产高清视频在线播放一区| 麻豆av噜噜一区二区三区| 国产主播在线观看一区二区| 国产 一区 欧美 日韩| 亚洲狠狠婷婷综合久久图片| 看十八女毛片水多多多| 欧美色视频一区免费| av天堂在线播放| 欧美丝袜亚洲另类 | 日韩欧美三级三区| 国产av麻豆久久久久久久| 美女xxoo啪啪120秒动态图 | 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 成人性生交大片免费视频hd| 午夜福利在线在线| 精品久久久久久成人av| 亚洲自偷自拍三级| 中文资源天堂在线| 老鸭窝网址在线观看| 1000部很黄的大片| 成人鲁丝片一二三区免费| 他把我摸到了高潮在线观看| 美女被艹到高潮喷水动态| 亚洲av一区综合| 欧美区成人在线视频| 18禁黄网站禁片午夜丰满| 亚洲精品成人久久久久久| 又黄又爽又刺激的免费视频.| 露出奶头的视频|