辐射传递理论partokb

第四章 海洋辐射传递理论

第一节 引言

海洋辐射传递，顾名思义，即为辐射在海水介质中受到散射与吸收所导致的辐射场变化。海洋光学辐射传递理论即是定量地研究辐射能通过海洋水体，受到多次散射和光谱吸收后，辐射场的空间分布及光谱分布的变化。

海洋辐射传递理论是海中能见度、对比度传输、水中图象传输、激光水中传输、海洋激光雷达、海面向上光谱辐射、海洋光学遥感、海水光学参数测量等应用研究的理论基础。它与近代光学技术、激光、光学遥感探测海洋的应用研究密切相关。因此海洋辐射传递理论是海洋光学基本理论和理论核心。辐射传递又是天体物理和大气光学的重要理论工具，因为电磁波（包括核辐射）与物质相互作用的研究是近代物理的重要组成部分，故辐射传递是近代物理的重要工具，因此海洋辐射传递的研究同时也具有更普遍的理论意义。 海洋辐射传递的基础问题大致可划分为：

1）经典问题也称为辐射传递正问题，即已知海中空间各点的固有光学性质和边界面的辐射场,求海中的辐射场分布。

2)第二类问题(又称“逆问题”)，即已知海中辐射场分布,求海水固有光学性质的参数。它是遥测海表层光学参数的理论基础，也是光学遥感测定海中叶绿素、悬移质和有机溶解物的基础。

3)窄光束问题，主要是求解高方向性激光束在海中的传输。它是海洋激光雷达、激光水下--空中通讯应用的理论基础。

4)海洋--大气系统辐射传递问题，即在建立海洋--大气系统辐射传递模型基础上, 根据大气顶所接收到的辐射推算海表面辐射。

5）水下图象传输问题，研究水下目标通过水体后图象的模糊和变化，或归结为海中点扩展函数和光学传递函数理论问题。

按照大气光学、海洋光学中的辐射传递模型，辐射传递方程可写为

⎰+-=πωθθβθ4''),()',(d r L cL dr

dL （4-1） 这里，L 为辐亮度，c 为体积衰减系数，β为海洋水体的体积散射函数，图4.1为海洋中辐射传递物理模型的示意图。显然，方程（1）是一种微—积分方程，因为β函数的复杂性，方程难以解析求解。目前国际上对辐射传递问题的研究，主要有三种方法：

1、近代解析求解

2、分离坐标法（主要是球谐函数方法）

3、蒙特卡罗方法（Monte Carlo Method ）

这三种方法在国际上一直延用至今，比较有效的方法是蒙特卡罗方法，也是当前受人关注的方法。第一种方法一般作为理论条件下的研究，后两种方法都必须进行数值计算，计算量较大。这三种方法均未突破辐射传递积一微分方程所固有的解方程的困难，当前国际上海洋光学、大气物理、天体物理、中子迁移物理所进行的辐射传递研究几乎都徘徊于如何用数值模拟方法对方程直

接求解。如果能建立一套更筒洁、更精确的物理模型，这将是海洋辐射传递理论的重要发展。

将方程（4-1）的各项用球谐函数展开，进行数值求解是一种好的方法。近年来，因为计算机的发展，用Monte Carlo 方法模拟海洋中的辐射传递过程，取得了较好的结果。几十年来，有关海洋光学和大气光学的辐射传递问题，发表了大量的论文。但是，辐射传递实质上是一种信息传递过程，这一重要的概念一直未受到注意。如果运用线性系统理论、Fourier 光学方法讨论海洋中辐射传递问题，对海洋光学中许多较复杂的问题，如传递函数、海洋辐射传递逆问题、垂直衰减K 函数等，大大简化，并可得到很好的结果。

第二节 海洋两流辐射传递理论

一、海洋两流辐射传递微分方程

海洋水体一般认为是一种水平平面分层介质。两流辐射传递理论模型简单地将通过水平分层的辐射通量分为向上辐照度E u (z )和向下辐照度E d (z )两个方向的光子流，海洋光学中将此称为两流辐射传递模型（见图4.2）。 我们来讨论向下辐照度E d (z )随深度的变化， z 取向下为正。通过水层△z ， 因为海水的吸收所造成的辐照度E d (z )的衰减可表示为

⎰-=π

πθθθπ2sin cos 2)(d dr aL z dE da (4-2)

因为dz dr =θcos (r 矢量的方向取光子流方向)，所以

⎰-=-=-=ππθθπ2

0)()(sin 2)(z aDE z aE d aL dz z dE d d da (4-3) 同理，海水回向散射所造成的辐照度E d (z)的衰减可表示为

)()(z bDE dz

z dE d db -=(4-4) 可见，向下辐照度通过水层△z 受到了散射、吸收两个过程的作用而衰减。另外，因为向上辐照度E u (z)受水体的回向散射作用，成为向下辐照度E d (z)的增量，即bE u (z)d z ，因此，向下辐照度和向上辐照度通过△z 水层的变化率可分别表示为

())()()(z bE z E b aD dz

z dE u d d ++-= (4-5) ())()()(z bE z E b aD dz

z dE d u u ++-=- (4-6) 以上两式称为两流辐射传递微分方程。

当水深足够深时，其解为：

Kz d d e E z E -=)0()( (4-7)

Kz u u e E z E -=)0()((4-8) 其中：)2(b aD aD K +=，为辐照度衰减系数，是表征海中辐照度随深度增加而衰减的比例因子。D 为分布函数，它表征辐射场分布的漫射特性。辐射场分布的斜射光越强，则D 越大。由此可见，水中向上、向下辐照度随深度 z 而呈指数衰减。同理，水中标量辐照度随深度z 也为指数衰减其通解为

kz kz d e g m e g m z E -+--++=)( (4-9)

kz kz u e g m e g m z E ---+++=)((4-10)

式中m 为由边界条件确定的常数。其中

)2(b aD aD k +=

k

aD g ±=±1 二.海底无限深时的水中辐照度

因向上、向下辐照度是有限值，当水深z →∞时，必然要求式(4-9)和(4-10)中的m +=0，故当海底为无限深时，

)ex p()(kz g m z E d -=+-(4-11)

)ex p()(kz g m z E u -=--(4-12)

其辐照比为

aD

k aD k g g z E z E R d u +-===∞+-)()()( (4-13) 即)(∞R 与z 无关。

水中向上、向下辐照度随深度 z 而呈指数衰减

)

exp()0()()

exp()0()(kz E z E kz E z E d d u u -=-=(4-14)

水中标量辐照度随深度z 也为指数衰减 [])ex p()0()()()()()(0000kz E z E z E D z E z E z E d u d u -=+=+= (4-15) 水中辐照比及透射率的几个极限值为：