基于海洋安全的大氣波導觀測方法研究

汪文杰, 許佳立, 黃賢青, 賈東寧, 2

基于海洋安全的大氣波導觀測方法研究

汪文杰1, 許佳立1, 黃賢青1, 賈東寧1, 2

(1. 青島海洋科學與技術試點國家實驗室, 山東 青島 266237; 2. 中國海洋大學, 山東 青島 266010)

海上船舶一旦遇到大氣波導, 可使船載電子信息系統之間的短波通信距離增大數倍, 同時通信距離范圍內的部分區域會出現電磁盲區。本文圍繞海洋安全保障需求, 開展了大氣波導協同觀測技術研究, 闡述了海上大氣波導對超視距傳輸與探測的重要性。通過對大氣波導的形成機制和影響機理進行研究, 系統化概述了蒸發波導、表面波導和抬升波導的觀測原理及常規觀測方法。通過分析研究傳統大氣波導的觀測方法, 發現其在觀測手段、觀測精度、觀測范圍和觀測成本方面存在局限性。基于此, 本文提出一種多物理場協同探測的大氣波導觀測方法, 通過集合海上浮標網、船聯網、星聯網和岸基網等各類觀測儀器和設備, 可實現海上大氣參數的實時獲取, 構建海上大氣波導預報系統, 形成海上大氣波導大范圍、高精度、實時性預報, 為海上超視距雷達探測提供關鍵技術支撐。

大氣波導; 海洋安全; 電磁傳播; 超視距探測; 協同觀測

大氣波導是一種雷達波超視距傳輸的現象, 是由一定的對流層氣象條件產生的具有垂直分布和水平擴展性的一種特殊大氣層結構。在海洋環境中, 由于受海面以上水汽垂直分布以及大氣邊界層穩定度的影響, 使得大氣波導在海洋中出現的頻率和持久度要遠高于其在陸地環境中出現的情況[1-2], 例如海上蒸發波導, 其出現的概率高達85%以上。隨著世界各國對海洋經濟的發展越來越重視, 海洋安全逐漸成為海洋經濟持續穩定發展的重要支撐。因此, 利用海上大氣波導來實現目標的超視距探測與超視距接收將有利于為海洋安全保障。

本文主要利用大氣波導的雙面性, 從超視距雷達探測和利用盲區隱身突防兩個角度解析了海上蒸發波導、表面波導和懸空波導對海洋安全的重要性, 并詳細闡述了這幾種大氣波導的探測原理與傳統觀測方法, 通過對比分析發現現有觀測方法存在一定的局限性, 無論是觀測精度還是觀測范圍都無法滿足超視距傳輸需求。為此, 本文提出了一種基于多物理場協同觀測的海上大氣波導觀測方法, 該方法集合了海上眾多觀測儀器設備, 不僅提高了觀測精度, 還實現了大范圍觀測, 解決了傳統觀測方法存在的局限性。通過實時獲取目標海域的大氣波導參數, 如波導高度、強度、厚度和范圍等, 可指導海上船舶進行超視距雷達精準探測。

1 國內外研究現狀

關于大氣波導的研究最早可追溯至20世紀30年代,當時雷達技術已得到廣泛應用, 由大氣波導引發的奇怪現象也逐漸受到世界各國的重視。20世紀70年代, 隨著軍事競爭的進一步需求, 世界上又掀起了研究大氣波導的高潮。由于大氣波導具有超視距發現、易形成電磁盲區、跳躍盲區[3]并增加雷達測高、測距、測角、測速的誤差以及增強雷達雜波等特點, 如圖1所示, 因此, 各國開始競相開展海上大氣波導觀測技術相關研究, 以便掌握大氣波導的強度和位置, 進而為海上目標探測、發現和識別提供技術支持, 把握海上安全主動權。

目前, 國外對大氣波導的探測方法主要包括衛星探測、低空系留氣球探測、氣象梯度塔探測、微波折射儀探測、微波輻射計探測、多普勒擴頻雷達探測等方式。其中, 美國于1976年開發出集成大氣折射率預報系統。美國在20世紀80年代投入的典型雷達診斷系統包括EMPE、AREPS、RPO以及RFPPAS等[4], 此類系統依托大氣波導的有利傳播條件來提升艦載雷達系統、偵測系統、岸基雷達系統以及通信系統等裝備的性能。美國對大氣波導的研究一直致力于其在軍事上的應用, 部分學者通過X波段信標接收機陣列系統推導蒸發波導的高度, 并與大氣垂直表面層模型一起使用, 以獲得蒸發波導高度的精確值[5]。1995年, 美國Micro Lab1號[6]低軌衛星上天, 首次實現了利用地球掩星對大氣波導的觀測。前蘇聯也在各大洋中做了大量實驗性工作, 結果表明在特定情況下, 由海上蒸發波導引起的電磁特殊效應可使雷達探測范圍延伸至1 000 km以外, 澳大利亞[7]通信愛好者也多次收聽到1 000 km外的超短波信號。除此之外, 日本、英國、法國、荷蘭和印度等國家都相繼開展了大氣波導觀測技術研究, 并取得一定成果。

圖1 大氣折射示意圖

中國在大氣波導方面的研究工作起步較晚, 最初的研究工作其主要方式是通過釋放探空氣球或發射探空火箭進行觀測, 從而獲取大氣溫度、濕度、氣壓等相關參數隨高度的變化情況, 進而探明對應觀測點出現波導傳播的條件。2000年前后, 中國學者開展基于模型的研究, 焦林等[8]提出了一種基于拋物線方程法的波導盲區位置確定方法, 并開發了盲區預測系統。康士峰等[9]提出了一種基于目標函數的雷達天線高度優化方法, 針對特定區域不同目標函數給出了雷達天線高度的優化結果。左雷等[3]對蒸發波導和表面波導條件下的雷達盲區形成機理及特性進行了分析, 并形成了大氣波導對艦載雷達盲區影響的初步成果。近年來也有相關科研機構開始應用地基雷達、機載雷達等觀測設備開展雷達波反演與仿真計算研究[10-11], 相關研究結果表明國外大氣波導的計算模型并不適應中國海域[12], 應該根據中國海域的季節特點來選擇最佳的波導計算模型并進行參數優化, 開發滿足中國海域特殊環境特點的波導模型。

雖然各國對大氣波導的探測方法呈現多元化, 但由于不同的探測方式所獲取的信息參數不同, 采用的計算模型也不同, 故均存在一定局限性。其中, 衛星探測大氣波導是通過地面接收站接收的衛星信號來探測大氣波導存在與否; 低空系留氣球和氣象梯度塔探測方式采用探測大氣濕度、溫度、氣壓等參數的方式間接獲得大氣折射率; 而微波折射儀和微波輻射計等設備則是直接探測大氣折射率來獲取大氣波導相關參數[13]。隨著衛星觀測技術和相關探測儀器的進一步發展, 國際上出現了一些新的大氣波導探測方法:

1) 利用GPS氣象學獲得海上氣象參數廓線, 開展大氣波導探測研究, 即利用低軌衛星的星載GPS接收機和岸基雷達接收站的探測結果反演氣溫、氣壓、濕度等參數, 用于大氣波導探測;

2) 利用雷達接收的海雜波反演大氣波導剖面;

3) 利用激光雷達測得大氣環境參數, 用于海上大氣波導探測;

4) 利用數值天氣預報和數據融合與同化技術, 結合海洋遙感衛星觀測數據進行大氣波導預測預報;

5) 利用岸基已有的電站, 如電視、電臺和移動通信基站等接收的信號對大氣波導進行反演, 從而獲得大氣波導相關參數[14]。

2 大氣波導傳統觀測方法研究

海洋大氣環境中一般會出現三種形式的大氣波導: 蒸發波導、表面波導和抬升波導。由于海洋表面大氣濕度在垂直方向變化急劇, 故蒸發波導出現的概率最高, 而表面波導和抬升波導出現的概率相對較低。因此, 掌握這三種大氣波導的形成機理和探測原理對大氣波導的進一步研究至關重要[15]。下面將詳細分析蒸發波導、表面波導和抬升波導的探測原理, 通過分析各種觀測方法的優缺點進一步完善大氣波導觀測理論體系, 同時提出一種多物理場協同觀測的探測方法, 該方法集合各類海洋觀測儀器、設備以及衛星數據, 可實現大氣波導的高精度、全覆蓋、全天候、實時性觀測與預報, 從整體上把握中國制海權。

2.1 蒸發波導

由于海面上蒸發的水汽會在貼近海平面一定高度的氣體層達到飽和, 導致其濕度在很小的高度范圍內產生很大的垂直梯度, 從而形成波導傳播的條件, 由此產生的波導稱為蒸發波導。這是一類比較特殊的表面波導, 不受時間和空間的限制, 可以在任何時間、所有海域產生蒸發波導, 是出現概率最高的一種大氣波導。蒸發波導高度通常在離海面40 m以下[16], 這一高度正好覆蓋了大多數艦載微波雷達的天線高度, 可以有效利用蒸發波導進行遠距離目標探測[17-18]。

目前, 蒸發波導的探測方法主要包括直接測量法、遙感反演法和數值模擬法。其中, 直接測量法目前主要以激光探測為主, 利用幾何光學理論對大氣中電磁波進行射線追蹤, 從而進行大氣波導的推導; 遙感反演法主要通過獲取雷達海雜波和GPS信號, 并分析其中所包含的大氣參數要素, 進而推導出蒸發波導折射情況; 數值模擬法則主要通過波導理論模型來計算蒸發波導的存在, 采用的波導模型主要有: Paulus-Jeske模型[19]、MGB模型[20]、偽折射率模型[21]、RSHMU模型和拋物線模型等[22-23]。

激光探測主要是利用激光的電磁波特性, 其在大氣中傳播時會受到大氣折射率的影響發生折射, 激光折射率與激光波長、大氣溫度、大氣濕度和大氣壓強有關, 其關系如下:

=1+(,,,). (1)

由于激光在大氣中的折射率和消光系數的變化率與上述大氣參數密切相關, 因此可以通過激光的折射率和消光系數變化率來推導蒸發波導的特性, 基本理論框架如圖2所示。通過建立激光在大氣中傳輸特性和蒸發波導之間的邏輯關系, 便可利用激光的回波信號對蒸發波導進行預測。同時, 由于激光傳輸具有單向性的特點, 使得激光探測蒸發波導更具隱蔽性, 在海上復雜電磁環境下, 可以在雷達保持靜默條件下實現對蒸發波導的探測。采用雷達進行海上目標探測易被干擾, 而激光的單向傳輸特性使其不易被干擾和截獲。

圖2 激光探測大氣波導理論框架

遙感反演法通過分析接收到的電磁波信號所承載的大氣要素信息來求解大氣折射狀況。其中, 雷達海雜波反演信號主要通過艦載雷達和岸基雷達等設備獲得, 因為此類設備安裝高度大部分位于蒸發波導高度范圍之內, 能夠準確捕獲海雜波信息, 通過將接收到信號和電磁波傳播模式正向模擬結構作比較, 即可獲得大氣折射情況, 從而推導出蒸發波導的存在。其中最常采用的算法和技術包括最小二乘法、模擬退火法、最大后驗概率算法、降維技術、遺傳算法和射線追蹤技術[24]。

GPS反演則是利用岸基固定GPS探測系統進行蒸發波導推導, 通過分析GPS探測系統接收的電磁信號, 求解出GPS信號因蒸發波導存在而產生的時延, 將時延作為參數輸入到蒸發波導計算模型中, 通過分析時延與氣象要素剖面關系進行蒸發波導反演。因岸基GPS系統分布廣泛、位置固定, 因而可以實現蒸發波導全天候、高可靠、大范圍觀測。

數值模擬探測主要是通過一種基于相似理論的經驗模型, 首先通過海上的智能浮標、漂流浮標、船載設備、低空系留氣球以及梯度塔等相關探測儀器獲取海上大氣濕度、溫度和壓強等氣象參數, 并以此作為輸入導入到相關波導模型中推導出大氣折射率, 最后通過大氣折射率預測大氣波導及其他特征參數[25]。具體可分別根據(2)、(3)式計算出大氣折射率并診斷波導是否存在, 當滿足(3)式時即可判別存在大氣波導, 最后根據(4)式來判斷大氣波導強度, 公式如下:

其中,為考慮到地球曲率的影響加以訂正的大氣折射指數;為海表大氣溫度, 單位為K;為水汽壓, 單位為hPa;為大氣壓, 單位為hPa;為海拔高度, 單位為m;0為地球平均半徑, 取值6.371×106m。Δ表示的是波導強度, 即波導陷獲層頂部與波導底部大氣修正折射指數之差。

由于蒸發波導的波導高度一般在40 m以下, 通過地基雷達接收站觀測蒸發波導受波導高度影響較大, 相關學者開始通過改進的離散混合傅里葉變換方法求解用于蒸發波導的拋物線模型, 并基于此模型開發了大規模的波導損耗模型, 最終獲取了接收天線高度和波導高度對波導傳播損耗的影響關系[26]。

2.2 低空波導(表面波導和抬升波導)

不同于蒸發波導貼近海面的特殊性, 表面波導和抬升波導的性質及觀測方法基本相同, 只在形成機理上略有區別: 表面波導是由大陸上暖而干燥的氣團進入到相對冷而潮濕的海面形成的, 其波導厚度一般在300 m以下; 抬升波導主要是由副熱帶高壓氣體的運動形成的, 在副熱帶和熱帶海域的高層大氣中存在大面積的強烈下沉運動, 當干而熱的氣團下降到一定高度碰到冷而濕的海氣時, 會直接覆蓋在海氣邊界層上方形成逆溫, 從而形成抬升波導, 抬升波導的高度一般在3 000 m以下。因表面波導向上爬升一定高度后即成為抬升波導, 因此, 本文不再區分表面波導和抬升波導的觀測區別, 一并進行研究, 并統稱為低空波導。

由于低空波導的下邊界層發生變化, 不再是海平面, 而是具有一定高度的大氣層, 電磁波被陷獲在一定高度的大氣層中傳播, 地基接收站無法接收到電磁信號, 因此低空波導無法像蒸發波導那樣利用地基雷達接收海雜波來反演大氣波導。目前, 對于低空波導主要利用GNSS系統進行觀測, GNSS系統主要包括美國的GPS系統、俄羅斯的GLONASS系統和中國的北斗系統等, 可通過GNSS機接收到的衛星仰角、方位角等信息來判定大氣波導的存在。原理為通過地面接收站接收到的衛星的仰角信息來判斷此時的衛星是否在視距范圍之內, 如果在視距范圍之外便接收到衛星的信息則說明發生了波導傳播, 再根據它們之間的射線幾何關系便可判斷大氣波導的長度, 原理如圖3所示。根據圖3中的幾何關系可以看出: 當<時, 為可視范圍; 當>時, 若接收站接收到了衛星信號, 則說明在此空間范圍內存在大氣波導[27]。

根據圖3的幾何關系可以推導出:

圖3 大氣波導長度計算示意圖

其中,為地球球心,為地球半徑,為天線高度,為GNSS衛星,為天線,為GNSS衛星到地球中心距離,為GNSS衛星仰角,為截止角度。

2.3 大氣波導傳統觀測方法存在不足

1) 觀測精度不高

在大氣波導探測方面普遍存在的現象是探測精度不高, 對于蒸發波導主要通過低空系留氣球和梯度塔等設備進行觀測, 只能對定點高度進行探測, 無法保證在垂直方向上的探測精度。目前, 世界各國正在努力提高大氣波導的探測精度, 其中不萊梅大學的von Engeln等[28]研究人員利用歐洲的中尺度大氣預報模式數據進行了一次全球大氣波導研究, 并推導出大氣波導在不同區域出現的概率、波導高、波導厚度等參數, 由于該大氣模式預報和探空數據的精度只有200～300 m, 只能發現一些厚度較大的波導層, 因此無法進行精細化預報。在衛星數據反演大氣波導技術上當受到不同高度的波導影響時也會出現精度偏低的情況, 因為單波段雷達雜波只對特定高度范圍內的波導反演具有較高精度, 當遇到復合波導層時則須利用多波段雷達雜波反演技術推導大氣波導[29]。

2) 觀測范圍有限

目前對于大氣波導的觀測儀器主要包括探空氣球、梯度塔以及遙感衛星等, 而探空氣球和梯度塔作為常規觀測工具雖然精度較衛星稍高, 但其觀測的范圍非常小, 屬于點對點觀測。且梯度塔和低空系留氣球多建設在海岸線, 只能對陸地和近海海域的大氣環境進行觀測, 而對于遠海海域還是一片空白, 無法對遠海的電磁環境進行預測, 使得遠海環境無法得到保障。雖然采用遙感衛星觀測大氣波導可實現大范圍觀測, 但其垂直分辨相對較低, 因此考慮到精度的制約, 該技術一直停留在定性分析上, 并未得到實際應用。而對于其他地基雷達觀測站一般存在分布間距大、觀測參數單一等缺點, 亦無法滿足海上大氣波導大范圍觀測需求。

3) 觀測成本高

原始的大氣波導觀測方法都是通過探空氣球進行觀測的, 多采用投放式, 大多無法回收, 導致收益比較低。而若用無線電探空儀、微波折射儀和地基雷達觀測站等設備進行觀測, 則其測量和制作成本又很高, 不適合大面積布放。采用衛星信息參數反演大氣波導可實現大面積觀測, 但其運行維護成本高, 存在技術不成熟和擴展性差等缺點, 且一般只供軍方使用, 民用衛星普及率還很低, 亦無法滿足觀測需求[30]。

綜上, 大氣波導傳統觀測方法雖然技術已進入實踐應用, 但仍存在觀測精度低、觀測范圍小和觀測成本高等缺陷, 急需一種簡單、高效、實時和低成本的觀測方法實現對大氣波導的觀測。

3 大氣波導多物理場協同觀測方法研究

針對傳統大氣波導觀測方法存在的不足, 同時為滿足海上作業需求, 提出一種以浮標網、船聯網、星聯網和岸基網為觀測設備的多物理場協同觀測方法, 海上作業指揮流程圖如圖4所示。該立體觀測網絡充分利用海上智能漂流浮標、船載智能終端設備、GNSS系統、地基雷達接收站和電臺等觀測設備, 構建了從點到線、從線到面、再從面到體的“空-天-地-海”立體觀測體系, 形成了海氣界面以上3 000 m內的大氣空間“溫-濕-壓”等氣象參數網絡。同時, 本研究利用船載計算終端和海洋試點國家實驗室的超級計算機, 充分發揮邊緣計算和云計算優勢, 實現了大氣波導觀測方式的設備協同、數據協同和計算協同。因協同觀測方法充分利用了海洋試點國家實驗室海洋星簇、海氣表面、海洋模擬器等大科學計劃, 獲得了涵蓋西太平洋、南印度洋以及中國海等海域的大量海洋氣象數據, 海洋模擬器可做到大氣預報精度27 km×27 km, 洋流預報精度9 km×9 km。其中為助力我國“海上絲綢之路”的發展, 針對“兩洋一海”區域進行了精細化預報, 可做到大氣預報精度9 km×9 km, 洋流預報精度3 km×3 km, 充分解決了傳統單點式觀測方法精度不高、觀測范圍有限和觀測成本高等缺點。

圖4 超視距探測流程圖

大氣波導傳統觀測方法和協同觀測方法的對比分析如表1所示。從表1可以看出, 由于協同觀測方法依托“空-天-地-海”立體觀測網絡, 開展大氣波導仿真模擬, 因此其在預報精度和范圍方面都優于傳統觀測方法。同時, 協同觀測方法不增設大氣波導專用觀測設備, 而是利用現有觀測網絡采集數據, 只在數據傳輸方面有所損耗, 因此觀測成本相對較低。大氣波導協同觀測方法的初始場數據均通過終端設備進行數據清洗之后再實時回傳至服務器, 數據的時效性較高, 且依托全球海洋模式系統可對未來15 d的波導條件進行預報, 為海上航行安全提供了強有力的支撐。

表1 大氣波導觀測方法對比分析

截止目前, 海洋試點國家實驗室海洋大數據平臺擁有40 PB的存儲容量, 且實驗室組建的科考船共享航次、浮標網、潛標網等觀測網絡體系持續不斷的將數據實時回傳至大數據中心, 通過緊密結合現有存量數據和實時增量數據, 完成了海洋氣象數據智能獲取與傳輸、重構與融合、存儲與管理關鍵技術研究, 實現大數據技術對氣象數據的服務管理。鑒于大氣波導協同觀測方法獲取的數據分布范圍廣、價值密度低等特點, 通過人工智能技術對獲取的數據進行智能過濾、智能篩選和智能壓縮, 提高了數據傳輸效率, 降低了數據傳輸成本。最后, 利用神經網絡技術, 以協同觀測網絡獲取的“兩洋一海”數據為輸入, 通過預先設定好的網絡神經單元和預期誤差, 對大氣波導模型進行強化訓練, 不斷優化大氣波導預測預報系統, 解決不同波導模型在中國海域不適用問題[12], 為適合中國海域的特殊波導模型研究提供了數據支撐和模型支撐。

最終, 利用大數據、人工智能和神經網絡技術對大氣參數進行了實時有效的分析和整合, 為海上大氣波導觀測提供了高時空分辨率數據支撐[31], 構建了大氣波導協同觀測網絡架構, 為海上超視距雷達探測提供指導, 系統架構如圖5所示。

圖5 大氣波導協同觀測系統框架

其中, 對于波導高度的求解, 本文提出的協同觀測方法與傳統波導模型的預測方法略有不同。傳統波導高度預測大多通過構建大氣折射率廓線來推導波導高度, 對氣象要素瞬時值依賴較高, 而氣象瞬時值往往難以準確快速獲得, 導致波導模型預測與實際情況偏差較大。而大氣波導協同方法則是根據海洋大氣近地層相似理論, 利用海氣界面宏觀觀測物理量來計算蒸發波導及其特征量, 并引入風速、溫度以及濕度普適函數來進行波導高度計算。最終通過求解溫度、大氣壓和水汽壓的廓線并結合大氣修正折射率公式即可確定波導高度。

溫度垂直廓線如下式:

大氣壓廓線由靜力學公式和理想氣體定律聯合積分如下:

水汽壓垂直廓線如下:

該協同觀測方法深度融合了立體觀測網絡的物理場數據, 而浮標網、船聯網、星聯網和岸基網作為多物理場協同觀測網絡的重要組成部分, 主要負責獲取前端空間立體網絡的感知數據, 并利用自組網技術對獲取的信息進行數據融合, 用于海上大氣波導預測預報, 多物理場協同感知網絡如圖6所示。

圖6 多物理場協同感知網絡

1) 浮標網

以海上智能漂流浮標、氣象浮標、固定觀測浮標為載體, 利用人工智能、大數據、自組網等新一代信息技術, 基于多傳感器同步觀測技術、智能邊緣計算技術、多源數據智能感知與預處理技術、智能數據壓縮編碼技術、自適應通信組網技術以及海氣耦合與數值計算等關鍵技術研究成果, 構建多物理場、自組網、低功耗、低成本、長時序智能浮標觀測系統, 為蒸發波導檢測提供設備支撐[32]。

2) 船聯網

開展基于IEEE1451標準協議的網絡化智能化標準接口的設計、船基數據傳輸鏈路的構建方式研究, 并設計完成海洋觀測智能終端系統, 該系統集合了全球船舶體量龐大的優勢, 形成全海域覆蓋的船基觀測網絡, 該智能終端設備不僅自身具備觀測能力, 還可通過自組網接收其他海洋裝備的觀測數據, 實現對海上大氣溫度、濕度、氣壓等關鍵大氣波導參數的全面采集, 為大氣波導預測提供堅實的數據支撐。

3) 星聯網

利用已有GNSS系統構建衛星觀測網絡, 因GNSS接收機可獲取多星信息, 故可將美國的GPS、俄羅斯的GLONASS系統和中國的北斗(BD)系統接收的信息進行有效融合, 使得信息參數更加準確。同時, 海上各類觀測設備均內置GPS定位模塊, 可在獲取溫、濕、風、壓等參數的同時記錄坐標位置, 最終結合GNSS接收機數據為低空波導提供大范圍、高精度及實時性氣象參數。

4) 岸基網

為減少觀測成本, 還可充分利用廣泛分布的岸基電站、電臺等信號反演大氣波導特征參數。由于沿海岸線分布的各類電站和電臺在收發信號時會受到海上大氣波導的影響, 使得FM或者TM電臺傳播方式由直線傳播變為弧線傳播或階躍傳播, 這將導致電磁信號的收發間隔和收發距離發生改變, 通過分析電磁信號強弱及電臺實際分布情況, 即可得知該信號是由對流層散射導致還是大氣波導影響導致。通過岸基觀測網絡對捕獲到的異常電磁信號進行反演, 推導電磁信號中含有的大氣波導特種參量, 即可用于大氣波導預報系統的分析研究。

最后, 本文基于Jeske模型、P-J模型、A模型、偽折射率模型及MGB模型等多種大氣波導數值模型及海雜波與電臺信號反演計算, 依托前端構建的立體感知網絡實時獲取的大氣波導氣象特征參數, 開展了大氣波導數據實時感知、參數智能篩選、模型深度優化、波導精準預報等關鍵技術研究, 并利用海洋試點國家實驗室海洋模式預報系統對未來15 d內海洋氣象條件進行精細化模擬, 模擬系統包括海洋大氣、降水、海溫和海流四個方面, 對海上復雜氣象條件進行預判預估, 如圖7所示。依托海洋模式預報系統, 結合多物理場協同觀測體系, 構建了大氣波導預測預報系統, 實現未來海上電磁偵察環境早感知、早預報和早部署, 如圖8所示。

圖7 海洋模式精細預報系統

圖8 大氣波導預報系統

4 結論

本文以大氣波導在海上航行安全的應用為背景, 分析了波導效應對海上超視距探測與利用電磁盲區的重要性, 以及如何利用大氣波導為海上空對空探測、預警、通信、截收等方面提供決策依據。通過對傳統大氣波導觀測方法的深入研究, 本文對蒸發波導、表面波導和抬升波導的觀測原理進行了對比分析, 研究了傳統觀測方法在觀測精度、觀測范圍和觀測成本上存在的局限性, 提出了一種新的大氣波導多物理場協同觀測方法, 充分利用了海洋試點國家實驗室立體觀測網絡對海上垂向3 000 m內全偵察環境實施高精度多物理量協同觀測網絡, 依托大氣波導高通量精細化預報系統, 實現了對海上大氣波導進行高精度、大范圍、實時性觀測, 強化了海上船舶的超視距探測能力, 為海上安全提供保障。

[1] ZHOU T, SUN T Y, HU H L. Analysis and prediction of 100 km-scale atmospheric duct interference in TD-LTE networks[J]. Journal of Communications and Information Networks, 2017, 2(1): 66-80.

[2] ZHAO X F, WANG D X, HUANG S X, et al. Statistical estimations of atmospheric duct over the South China Sea and the tropical eastern Indian Ocean[J]. Chinese Science Bulletin, 2013, 58(23): 2794-2797.

[3] 左雷, 涂擁軍, 姚燦, 等. 海上大氣波導環境下艦載超視距雷達盲區研究[J]. 火力與指揮控制, 2011, 36(10): 165-168.

ZUO Lei, TU Yongjun, YAO Can, et al. Preliminary investigation on the blind of shipborne OTH radar ba-sed on sea atmospheric duct[J]. Fire Control and Com-mand Control, 2011, 36(10): 165-168.

[4] 張玉生, 郭相明, 趙強, 等. 大氣波導的研究現狀與思考[J]. 電波科學學報, 2020, 35(6): 813-831.

ZHANG Yusheng, GUO Xiangming, ZHAO Qiang, et al. Research status and thinking of atmospheric duct[J]. Chinese Journal of Radio Science, 2020, 35(6): 813-831.

[5] POZDERAC J, JOHNSON J, YARDIM C. X-band beacon-receiver array evaporation duct height estimation[J]. IEEE Transactions on Antennas and Propagation, 2018, 66(5): 2545-2556.

[6] 林偉, 察豪, 史建偉, 等. 淺談海上大氣波導形成及其測量預報技術[J]. 微計算機信息, 2006, 22(1/2): 156-158.

LIN Wei, CHA Hao, SHI Jianwei, et al. The oceanic atmospheric duct form and its measure and forecast technology[J]. Microcomputer Information, 2006, 22(1/2): 156-158.

[7] 王華, 劉欣, 黃小毛. 表面波導對超短波通信對抗的影響[J]. 艦船科學技術, 2006, 28(6): 96-99.

WANG Hua, LIU Xin, HUANG Xiaomao. The effect of duct condition on ultrashortwave communication coun-ter-measure[J]. Ship Science and Technology, 2006, 28(6): 96-99.

[8] 焦林, 張永剛. 大氣波導條件下雷達電磁盲區的預報研究[J]. 西安電子科技大學學報, 2007, 34(6): 989- 994.

JIAO Lin, ZHANG Yonggang. Prediction of the electromagnetic shadow zone under the atmospheric duct[J]. Journal of Xidian Universiy, 2007, 34(6): 989-994.

[9] 康士峰, 曹仲晴, 王紅光, 等. 基于目標函數的微波超視距雷達天線高度優化方法[J]. 電子與信息學報, 2014, 36(7): 1767-1770.

KANG Shifeng, CAO Zhongqing, WANG Hongguang, et al. Antenna height optimization based on object function for microwave over-the-horizon radar[J]. Journal of Elec--tronics & Information Technology, 2014, 36(7): 1767-1770.

[10] GUO X M, ZHAO D L, ZHANG L J, et al,. C band transho-rizon signal characterisations in evaporation duct propagation environment over Bohai Sea of China[J]. Antennas Propagation IET Microwaves, 2019, 13(4): 407-413.

[11] 張鵬, 張守寶, 張利軍, 等. 區域海面蒸發波導預報與監測試驗對比分析[J]. 電波科學學報, 2017, 32(2): 215-220.

ZHANG Peng, ZHANG Shoubao, ZHANG Lijun, et al. Comparison between predicted and experimental result of regional evaporation duct over sea[J]. Chinese Journal of Radio Science, 2017, 32(2): 215-220.

[12] LIANG X, TIAN B, LV X S, YAN J W. Study on the best calculation scheme of equivalent evaporation duct feature[J]. The Journal of Engineering, 2019, 19: 5661- 5664.

[13] WANG Y J, GUO L X, LI Q L. A narrow-angle parabolic equation model in atmospheric ducts[C]// 2016 11th International Symposium on Antennas, Propagation and EM Theory (ISAPE), 2016: 404-407.

[14] WANG H, SU S P, TANG H C. Atmospheric duct detection using wind profiler radar and RASS[J]. Journal of Atmospheric and Oceanic Technology, 2019, 36(4): 557-565.

[15] 王海斌, 張利軍, 王紅光. 南海海區低空大氣波導氣候學分析[J]. 電波科學學報, 2019, 34(5): 633-642.

WANG Haibin, ZHANG Lijun, WANG Hongguang. The climatological analysis of the lower atmospheric ducts in South China Sea[J]. Chinese Journal of Radio Science, 2019, 34(5): 633-642.

[16] ESR Jinadasa, 田斌, 郭鵬. 蒸發波導模型概述[J]. 船電技術, 2018, 38(10): 13-15.

ESR Jinadasa, TIAN Bin, GUO Peng. Overview of evaporation duct models[J]. Marine Electric and Electronic Engineering, 2018, 38(10): 13-15.

[17] SHI Y, ZHANG Q, WANG S W, et al. Impact of typhoon on evaporation duct in the northwest Pacific Ocean[J]. IEEE Access, 2019, 7: 109111-109119.

[18] 熊威, 車夢虎. 海洋大氣波導環境下反艦導彈作戰效能評估[J]. 計算機測量與控制, 2019, 27(11): 252-254.

XIONG Wei, CHE Menghu. Combat effectiveness evaluation of anti-ship missile in ocean atmospheric duct environment[J]. Computer Measurement & Control, 2019, 27(11): 252-254.

[19] 宋偉, 田斌, 周沫, 等. 蒸發波導預測模型研究[J]. 華中科技大學學報(自然科學版), 2013, 41(5): 52-56.

SONG Wei, TIAN Bin, ZHOU Mo, et al. Research on evaporation duct predicting model[J]. Journal of Hua-zhong University of Science and Technology (Natural Science Edition), 2013, 41(5): 52-56.

[20] 田斌, 孔大偉, 周沫, 等. 蒸發波導迭代MGB模型適用性研究[J]. 電波科學學報, 2013, 28(3): 590-594.

TIAN Bin, KONG Dawei, ZHOU Mo, et al. Study on the applicability of the iterative MGB evaporation duct model[J]. Chinese Journal of Radio Science, 2013, 28(3): 590-594.

[21] 田斌, 察豪, 李杰, 等. 蒸發波導偽折射率模型特性[J]. 火力與指揮控制, 2010, 35(12): 73-75.

TIAN Bin, CHA Hao, LI Jie, et al. Study on the characteristic of pseudo-refractivity evaporation duct mo-del[J]. Fire Control and Command Control, 2010, 35(12): 73-75.

[22] DINC E, AKAN O B. Channel model for the surface ducts: large-scale path-loss, delay spread, and AOA[J]. IEEE Transactions on Antennas & Propagation, 2015, 63(6): 2728-2738.

[23] ZHU X Y, LI J C, JIANG Z H. An evaporation duct height prediction method based on deep learning[J]. IEEE Geoscience and Remote Sensing Letters, 2018, 15(9): 1307- 1311.

[24] 楊超, 陳競, 王一旨, 等. 雷達海雜波反演大氣波導的改進回溯搜索算法[J]. 系統工程與電子技術, 2018, 40(8): 1743-1749.

YANG Chao, CHEN Jing, WANG Yizhi, et al. Inversion of atmospheric duct from the radar sea clutter with the improved backtracking[J]. Systems Engineering and Electronics, 2018, 40(8): 1743-1749.

[25] YAN X D, YANG K D, MA Y L. Calculation method for evaporation duct profiles based on artificial neural network[J]. IEEE Antennas and Wireless Propagation Letters, 2018, 17(12): 2274-2278.

[26] YANG C, GUO L X. Inferring the atmospheric duct from radar sea clutter using the improved artificial bee colony algorithm[J]. International Journal of Microwave & Wireless Technologies, 2018, 10(4): 437-445.

[27] SHI Y, ZHANG Q, WANG S W, et al. A comprehensive study on maximum wavelength of electromagnetic pro-pagation in different evaporation ducts[J]. IEEE Access, 2019, 7: 82308-82319.

[28] VON ENGELN A, TEIXEIRA J. A ducting climatology derived from the European Centre for Medium-Range Weather Forecasts global analysis fields[J]. Journal of Geophysical Research，2004, 109：D18104，doi：10.1029/2003JD004380.

[29] 郝曉靜, 李清亮, 郭立新, 等. 基于MM5V3的大氣波導預報可靠性和準確度分析[J]. 電波科學學報, 2017, 32(6): 725-734.

HAO Xiaojing, LI Qingliang, GUO Lixin, et al. Reliability and accuracy to atmospheric duct prediction based on MM5V3[J]. Chinese Journal of Radio Science, 2017, 32(6): 725-734.

[30] DINC E, ALAGOZ F, AKAN Q B. Path-loss and correlation analysis for space and polarization diversity in surface ducts[J]. IEEE Transactions on Antennas and Propagation, 2016, 64(10): 4498-4503.

[31] ZHAO X F and HUANG S X. Radar sea clutter power modeling under the atmospheric duct propaga-tion- conditions[J]. Acta Physica Sinica, 2013, 62(9): 099204.

[32] 閆西蕩, 楊坤德. 蒸發波導環境下風浪對電磁波傳播影響的數值模擬研究[J]. 海洋與湖沼, 2020, 51(1): 13-19.

YAN Xidang, YANG Kunde. Numerical simulation of the influence of wind waves on electromagnetic wave propagation in evaporation duct environment[J]. Ocea-nologia et Limnologia Sinica, 2020, 51(1): 13-19.

Observation methods of atmospheric duct based on marine safety

WANG Wen-jie1, XU Jia-li1, HUANG Xian-qing1, JIA Dong-ning1, 2

(1. Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China; 2. Ocean University of China, Qingdao 266010, China)

The short-wave communication distance among shipborne electronic information systems will increase several times once the sea battlefield is exposed to the atmospheric duct. Meanwhile, the electromagnetic blind zone will appear within the communication distance range. Based on the strategic requirements of ocean security, this study carried out the reseacrch on the cooperative observation technology of atmospheric duct, and expounded the importance of atmospheric duct for over-the-horizon transmission and detection. We systematically summarize the observation principle of evaporation duct, surface duct and lifting duct and their conventional observation methods by studying the formation and influence mechanism of the atmospheric duct. In addition, we studid the traditional observation methods of atmospheric duct in depth and find that these methods have some limitations in observation means, accuracy, range and cost. Therefore, an atmospheric duct observation method based on multi-physical field cooperative detection is proposed in this paper, which can obtain real-time atmospheric parameters at sea and build an marine atmospheric duct prediction system by integrating various observation instruments and equipment, such as buoy network, ship network, satellite network and shore-based network. The large-scale, high-precision and real-time atmospheric duct forecasting system formed by our method is greatly significant to master air superiority in over-the-horizon radar detection and stealth penetration operations in the sea.

atmospheric duct; marine safety; electromagnetic transmission; transhorizon detection; collaborative observation

Dec. 30, 2020

O451

A

1000-3096(2022)1-0112-11

10.11759/hykx20201230002

2020-12-30;

2021-04-07

山東省重點研發計劃項目(2019JMRH0109)

[Key Research and Development Plan of Shandong Province, No. 2019JMRH0109]

汪文杰(1988—), 碩士, 工程師, 主要從事海洋科學、海工裝備研發以及海洋大數據等領域研究, E-mail: wangwenjiebkd@ 163.com; 賈東寧(1978—),通信作者, 博士, 高級工程師, 主要從事計算機信息科學、高性能計算與仿真以及國產超算典型應用等領域研究, E-mail: wjwang@qnlm.ac

(本文編輯: 趙衛紅)