鄧 榮，金曉宇，宋春霞

(海軍蚌埠士官學(xué)校，蚌埠 233012)

?

風(fēng)成海面激光反射特性分析

鄧 榮，金曉宇，宋春霞

(海軍蚌埠士官學(xué)校，蚌埠 233012)

為了計算低風(fēng)速下海面的激光反射特性，采用北海聯(lián)合海浪計劃(JONSWAP)非穩(wěn)態(tài)海譜模型，結(jié)合線性濾波方法，將入射激光光斑所覆蓋的海面分割為多個微小平面，再采用幾何光學(xué)方法，對入射到每個微小平面上的光線的反射光進(jìn)行理論計算，通過仿真，取得了二維風(fēng)成海面對激光反射的仿真結(jié)果。結(jié)果表明反射角在鏡面反射方向和入射反方向的±30°范圍內(nèi)，反射方位角在風(fēng)向或者反風(fēng)向2個方向，海面反射光能量較為集中。

大氣與海洋光學(xué)；反射特性；幾何光學(xué)；風(fēng)成海面

0 引 言

目前國內(nèi)外眾多學(xué)者在海面模擬以及海面的電磁散射方面已經(jīng)做了深入研究，發(fā)表了大量的文章[1-2]，而關(guān)于海面光反射的研究卻不多。同時他們的研究主要集中在激光雷達(dá)的應(yīng)用方面，即發(fā)射與接收在相同方向時的海面的電磁散射特性，而關(guān)于發(fā)射與接收系統(tǒng)分置時的海面光反射特性較少有人研究。本文采用線性濾波方法隨機(jī)生成二維風(fēng)成海面，利用幾何光學(xué)方法計算出海面激光的反射特性，研究了不同風(fēng)速下的計算結(jié)果。該研究對于風(fēng)成海面背景下的激光目標(biāo)探測具有一定的參考價值。

1 二維海面仿真

海洋表面是一種復(fù)雜的隨機(jī)粗糙表面，其表面的波浪起伏是在大尺度的近似周期性的波浪上疊加小尺度的波紋、泡沫和浪花。海面仿真通常是通過海譜實現(xiàn)的，目前存在多種海譜，分別從不同的角度對海面進(jìn)行描述，可以根據(jù)需要選擇合適的海譜來對海面進(jìn)行仿真[3]。本文選擇北海聯(lián)合海浪計劃(JONSWAP)非穩(wěn)態(tài)海譜模擬二維海面，其源于1968～1969年的聯(lián)合北海波浪計劃，是英、荷、美、德等國對北海進(jìn)行長期系統(tǒng)觀測后得出的，其具體表達(dá)式參看文獻(xiàn)[4]。

由海譜出發(fā)，應(yīng)用線性濾波方法可以產(chǎn)生隨機(jī)粗糙海面[5]。其基本思路是對海譜進(jìn)行濾波得到隨機(jī)海面高度起伏函數(shù)的傅里葉變換對，再對其進(jìn)行傅里葉逆變換得到隨機(jī)海面的高度起伏函數(shù)。因此，面積為Lx×Ly的二維隨機(jī)海面上每一點的高度起伏函數(shù)f(x,y)可由下式得到：

(1)

(2)

式中：M和N為等間隔離散點數(shù)；km和kn為離散波數(shù)，其表達(dá)式為km=2πm/Lx,kn=2πn/Ly；S(km,kn)為二維離散海譜；ωm,n為km與kn對應(yīng)的海譜時域頻率；N(0,1)表示均值為0、方差為1的正態(tài)分布的隨機(jī)數(shù)。

2 海面激光反射特性的仿真計算

采用上述方法模擬出二維隨機(jī)粗糙海面之后，將其劃分為許多微小面元，每一面元的局部坐標(biāo)系可以用2個相鄰邊的方向矢量x′、y′和其法向矢量z′來表示，即用x′、y′、z′組成局部坐標(biāo)系，其中每一面元的法向矢量由兩相鄰邊方向矢量叉乘求得。為了方便計算，需要將每個面元的局部坐標(biāo)系變換為本地基準(zhǔn)坐標(biāo)系，變換公式為：

(3)

式中：x、y、z為本地基準(zhǔn)坐標(biāo)系3個坐標(biāo)軸方向的單位矢量；α和β為x′、y′與水平面的夾角。

求解式(3)，則可以得到該面元的法向矢量為：

n=z′=xsinαcosβ+ycosαsinβ+zcosαcosβ

(4)

假設(shè)k為入射光方向的矢量，r為反射光方向的矢量，入射角和入射方位角分別為θi和φi，反射角和反射方位角分別為θr和φr，則入射方向單位矢量在基準(zhǔn)坐標(biāo)系中表示為：

(5)

根據(jù)反射定律，反射方向單位矢量可以由面元法向矢量與入射方向矢量求得，其表達(dá)式為：

r=k-2(k·n)n=xsinθrcosφr+

ysinθrsinφr+zcosθr

(6)

將入射激光光束看作由大量光線組成，依據(jù)公式(4)～(6)計算出每條光線的反射方向，將結(jié)果綜合起來就可以得到總體的隨機(jī)海面激光反射特性。

3 仿真分析

為簡化計算，本文不考慮入射激光光束的發(fā)散角，即將其看做平行且均勻的光束,構(gòu)造出面積為2 m×2 m的二維隨機(jī)粗糙海面，并將激光光斑處的海面分成邊長為1 mm的微小平面，再利用上節(jié)所述方法就可以計算隨機(jī)海面的激光反射特性。

影響海面激光反射特性的因素有很多，本文分析的是風(fēng)成海面的激光反射特性，主要以風(fēng)速對海面激光反射特性的影響為例進(jìn)行分析。需要說明的是，當(dāng)風(fēng)速較大時，海面的遮擋效應(yīng)會非常強(qiáng)烈，在不計算二次反射的情況下，計算結(jié)果誤差較大。因此本文只分析中低風(fēng)速下的海面激光反射特性。

當(dāng)取風(fēng)速U分別為5 m/s和10 m/s，風(fēng)向φ為0°，入射角θi為30°，入射方位角φi為0°，激光光斑直徑為100 mm時，分別仿真計算3次，計算結(jié)果如圖1、圖2所示。

如圖1所示，從反射角來看，在-60°～-10°和10°～60°的范圍內(nèi)反射光能量較為集中。從反射方位角來看，在0°和180° 2個方向出現(xiàn)了明顯的反射尖峰，反射光能量集中在這2個方向附近。

如圖2所示，隨著風(fēng)速的增大，海面高度起伏變化更加劇烈，因此從反射角來看，反射光的能量分布更加分散，沒有明顯的反射尖峰，但是在-60°～-10°和10°～60°的范圍內(nèi)反射光的能量仍然較為集中。從反射方位角來看，由于海面的遮擋效應(yīng)增強(qiáng)，在0°和180° 2個方向的反射尖峰更加明顯，反射光能量在這2個方向附近更加集中。

綜合圖1和圖2可以看出，反射角在鏡面反射方向和入射反方向的±30°范圍內(nèi)，反射方位角在風(fēng)向或者反風(fēng)向2個方向，海面反射光能量較為集中。同時隨著風(fēng)速的增大，反射光的反射角分布更加分散，而反射方向角卻更加集中。

圖1 U=5 m/s時的仿真結(jié)果

圖2 U=10 m/s時的仿真結(jié)果

4 結(jié)束語

本文采用幾何光學(xué)方法仿真計算了中低風(fēng)速下隨機(jī)風(fēng)成海面的激光反射特性，計算結(jié)果表明反射角在鏡面反射方向和入射反方向的±30°范圍內(nèi)，反射方位角在風(fēng)向或者反風(fēng)向2個方向，海面反射光能量較為集中。本文的計算結(jié)果對中低風(fēng)速下海面背景的激光目標(biāo)探測具有一定參考價值。

[1] 陳勇，侯德亭，齊國雷，刑召偉.分形海面電磁散射特性的基爾霍夫近似[J].信息與電子工程,2009,7(5):409

-412.

[2]ELFOUHAILYTANOSM，JOHNSONJOELT.Anewmodelforroughsurfacescattering[J].IEEETransactionsonGeoscienceandRemoteSensing,2007,45(7):2300-2308.

[3] 郭立新,王蕊,吳振森.隨機(jī)粗糙面散射的基本理論與方法[M].北京:科學(xué)出版社,2009:11-16.

[4]HASSELMANNDE.DirectionalwavespectraobservedduringJONSWAP1973[J].JournalofPhsicalOceanography,1980,10(8):1264-1280.

[5] 楊惠珍，康鳳舉，褚彥軍，聶衛(wèi)東.基于海浪譜的隨機(jī)海浪仿真及驗證[J].系統(tǒng)仿真學(xué)報,2005,17(10):2324-2326.

Laser Reflection Characteristics Analysis of Windy Sea Surface

DENG Rong,JIN Xiao-yu,SONG Chun-xia

(Bengbu Navy Petty Officer Academy,Bengbu 233012,China)

In the interest of calculating the laser reflection characteristics of sea surface under low wind speed,this paper adopts the unsteady sea spectrum model of joint north sea wave project (JONSWAP)to divide the sea surface covered by the incident laser beam into many small planes combining with linear filtering method,then uses geometrical optics method to calculate the reflected light of each incident light to each small plane theoretically,through the simulation,obtains the simulation results of the laser reflected light in 2-D windy sea surface.The results show that：the reflection angle is in ± 30° of the vertical direction or incident direction;the reflective azimuth is in the wind direction or against the wind direction;the reflected light energy is relatively concentrated.

atmospheric and ocean optics;reflection characteristics;geometrical optics;windy sea surface

2016-05-03

TN24

A

CN32-1413(2016)04-0025-04

10.16426/j.cnki.jcdzdk.2016.04.006