粗糙海面激光透射特性研究

李 恪,姚崇斌,谢宝蓉，周 勃

(上海航天电子技术研究所,上海 201109)

粗糙海面激光透射特性研究

李 恪,姚崇斌,谢宝蓉，周 勃

(上海航天电子技术研究所,上海 201109)

为改善激光在粗糙海面上透射特性仿真结果的准确性以提高卫星海洋水色遥感应的效率，对不同海况下激光在粗糙海面的透射特性进行了研究。根据激光海面透射原理，用基尔霍夫近似(KA)原理求解激光海面透射场，给出了计算模型。基于Gaussian随机模型对三、四级海况的粗糙海面进行了三维模拟仿真，用解析法计算了激光在海面上的透射特性。结果发现：波长0.532 μm的绿激光对海水有高透射率，入射角小于40°时透射率大于94%；海况对激光的透射率影响较小；与下行透过率相比，上行透射率受入射角的影响较大，入射角大于40°时激光的上行透射率迅速衰减，入射角大于50°时上行透射率几乎为0。

海洋水色遥感应； 粗糙海面； 激光； 海色； 透射特性； 基尔霍夫近似； 双站透射系数； 入射角

0 引言

海水颜色可作为各种能吸收和散射太阳光的海水组分(如浮游植物、有色可溶有机物、矿物粒子、海中气泡)含量的一个指示[1]。目前，卫星海色遥感中未考虑海中气泡这一重要的光学散射组分。气泡在上层海洋中无处不在，由于风浪破碎而形成的微气泡大量存在于水下0～10 m深度处，且海中气泡对蓝绿激光具强后向散射特性，因此针对不同海况，得出准确的激光海面透射系数，对提高卫星海色遥感的准确率有重要的意义[2-3]。之前的相关研究中，为便于计算，常将三维粗糙海面简化为二维粗糙海面，并用典型的随机函数模型近似海面的波动。因其建立的二维粗糙界面与现实中的粗糙界面差别较大，故仿真计算结果的准确性较差[4-5]。针对这种不足，本文采用了三维Gaussian随机模型，结合典型海况的特征参数，对不同海况的波动海面建立了三维仿真模型，并基于KA原理本文对不同海况下粗糙海面激光透射特性进行了研究，以提高激光在粗糙海面上透射特性仿真结果的准确性。

1 激光海面透射原理

图1 粗糙海面上激光的几何关系Fig.1 Geometry of laser diffusion from rough sea surface

假设入射激光的入射场Ei可表示为

(1)

(2)

(3)

式中：ω为激光的角频率；r，r′分别为接收点和海面光斑某点处的位置矢量；E(r′)，H(r′)分别为粗糙海面上r′∈S处的入射电场和磁场；G0(r,r′)，G1(r,r′)分别为空气和海水介质中的并矢格林函数，且

(4)

(5)

此处：k0，k1分别为激光在空气和海水中的波数；I0，I1分别为粗糙海面投影在x-y平面上各点的散射矢量和透射矢量的积分[6]。

2 基尔霍夫近似

粗糙海面的激光透射场Et(r)为随机场，可将其分解成平均场Etm(r)和起伏部分εt(r)，即

Et(r)=Etm(r)+εt(r).

(6)

根据KA原理，可将透射场表示为

(7)

3 数值计算与分析

因粗糙海面具随机分布特性，根据Gaussian随机模型对不同海况的随机海面进行模拟。Gaussian随机模型可表示为

(8)

式中：C″(0)为海面光斑两点间距离ρ=0时相关函数的二阶导数；h为粗糙海面的均方高度[10]。

基于Gaussian随机模型，分别对三级和四级海况进行了模拟，结果如图2所示。图2中：ξ为海浪高度。

图2 不同海况的粗糙海面模型Fig.2 Rough surface model in different sea condition

海面下方的局部法向矢量

(9)

将粗糙海面模型中的参数代入进行仿真计算，得激光束在不同入射角度的海面透射率计算模型为

(10)

式中：θi为入射角；η0，η1分别为空气和海水介质中的阻抗；a，b分别为透射光和入射光的偏振状态；Wab为激光束在不同介质中传输时的偏振系数[11]。不同海况时的激光海面透射率计算结果如图3所示。

图3 不同海况时的激光海面透射率Fig.3 Laser transmission ratio into rough sea surface indifferent sea conditions

对仿真计算结果进行分析，可发现：

a)波长为0.532 μm的绿激光对海水有高的透射率。当θi小于40°时，透射率可达94%以上。

b)海况对激光的透射率影响较小。这主要是因为相对海面波动来说，激光波长极短，激光光斑在任何海况下的任意时刻都可视作一个平面。

c)与下行透射率相比，上行透射率受θi的影响很大。这主要是因为激光束的下行通道是从光疏介质透射至光密介质，不会发生全反射，而上行通道是从光密介质透射至光疏介质，当θi大于临界角时会发生光的全反射。因此，当θi>40°时，激光的上行透射率开始迅速衰减；当θi>50°时，上行透射率几乎为0。

4 结束语

本文用Gaussian随机模型模拟不同海况的三维粗糙海面，根据KA原理对激光在不同海况下的海面透射率进行了仿真计算。研究发现：用Gaussian随机模型模拟不同海况的三维粗糙海面，可使仿真结果更符合实际情况，提高了仿真结果的准确性。采用KA原理与三维Gaussian粗糙海面模型结合的仿真计算方法，可在提高准确率的同时，极大地提高计算效率。激光能量在海面的透射率几乎不受海况影响，但受入射角的影响很大，当入射角大于某临界角(40°)时，激光的透射率将迅速减小。对激光海面透射率的研究对提高卫星海洋水色遥感效率，有一定的参考价值。后续研究中，将结合海洋实际情况，设定边界条件，用合理的近似方法简化数值计算的步骤，以进一步提高仿真计算的效率和准确性，为实际工程应用奠定基础。

[1] BELLRENFELD M J, RANDERSON J T, MECLAIN C R, et al. Biospherie primary production during an ENSO transition[J]. Seience, 2001, 291: 2594-2597.

[2] MONAHAN E C, LU M. Acoustically relevant bubble assemblages and their dependence on meteorological parameters[J]. IEEE Journal of Oceanic Engineering, 1990, 15: 340-349.

[3] ZEDEL L, FARMER D M. Organized structures in subsurface bubble clouds: Langmuir circulation in the open ocean[J]. GeoPhys Res, 1991, 100: 13309-13320.

[4] 任新成， 郭立新. 基于改进二维分形海面模型的分层海面电磁散射分析[J]. 上海航天, 2009, 26(4): 1-6.

[5] 闫沛文, 童创明, 邓发升. 不同风速下海洋粗糙面散射系数的计算[J]. 上海航天, 2006, 23(4): 23-26.

[6] 刘春香, 程传福, 任晓荣, 等. 随机表面散射光场的格林函数法与基尔霍夫近似的比较[J]. 物理学报, 2004, 53(2): 427-435.

[7] 李恪, 王江安, 潘琛. 激光海面漫反射解析特性研究[J]. 激光与红外, 2011. 41(5): 525-529.

[8] 田炜, 任新成. 基于基尔霍夫驻留相位近似法的深粗糙度高斯型粗糙面光散射[J]. 上海航天, 2011, 28(2): 1-6.

[9] 田炜, 姜志立. 高斯型粗糙面光透射的微扰法研究[J]. 科学技术与工程, 2012, 12(23): 5754-5757.

[10] 任新成, 郭立新. 基尔霍夫近似下高斯粗糙面透射波散射系数的特征研究[J]. 上海航天, 2008, 25(2): 1-5.

[11] TSANG L, KONG J A. Scattering of electromagnetic waves-numerical simulation[M]. New York: John Wiley, 2001.

Study on Laser Transmission Characteristics into Rough Sea Surfaces

LI Ke, YAO Chong-bin, XIE Bao-rong, ZHOU Bo

(Shanghai Aerospace Electronic Technology Institute, Shanghai 201109, China)

To improve the simulation accuracy of the laser transmission characteristic into rough sea surface for increasing the efficiency of remote sensing of ocean color, the laser transmission characteristics into rough sea surfaces in different sea conditions were studied in this paper. According to the principle of laser transmission into the sea surfaces, the laser transmission characteristics into rough sea surfaces in different sea conditions were computed by Kirchhoff approximation. The computation mode was derived. The rough sea surfaces in 3 and 4 level sea conditions were modulated based on the Gaussian random function. The laser transmission characteristics in the two sea conditions were calculated by the analytical method. It was found that the green laser with the wavelength 0.532 μm had high transmission into the seawater which was bigger than 94% when the incident angle was smaller than 40°. The influence of the sea condition on the transmission was small. The influence of the incident angle on the up reflectance was bigger than down reflectance. The up transmission would be attenuated quickly when the incident angle was bigger than 40° and near to zero when the incident angle was bigger than 50°.

Watercolor remote sensing laser; Rough sea surfaces; Laser; Ocean color; Transmission characteristics; Kirchhoff approximation; Bistatic transmission coefficient; Incident angle

1006-1630(2016)06-0078-04

2016-02-18；

2016-09-17

国家自然科学青年科学基金资助(61605110)

李 恪(1980—)，男，博士，主要研究方向为光电遥感技术。

TN249

A

10.19328/j.cnki.1006-1630.2016.06.011