太湖梅梁湾下行漫衰减系数的Monte Carlo 模拟及分析 - 遥感学报.PDF

太湖梅梁湾下行漫衰减系数的Monte Carlo 模拟及分析 1 2 姜勃 ，赵巧华 1.南京信息工程大学地理与遥感学院，南京 210044 2.南京信息工程大学水文气象学院，南京 210044 摘要：假设水体光学性质均一，固有光学特性（IOPs）不随深度变化，水体表面为镜面，通过Monte Carlo 方法模拟光子在水中的运动过程，统计大量光子在不同深度上的能量分布，计算出下行漫衰减系数 Kd 值。结果 表明：在450,550 和650nm 波长处模拟的Kd 值基本上小于实测值，对应后向散射概率为3%的散射相函数 （SPF1） 模拟的相对误差最小，适合用于构建太湖水体的Kd 值遥感估算模型。进一步分析入射天顶角、散射与吸收系数 的比值（b/a ）对Kd 值变化的影响，结果发现：在b/a 较大的情况下，散射相函数成为影响Kd 值变化的主要因素。 0 0 当b/a 值大于9 时，入射天顶角从10 变化到80 时，Kd 值的变化率小于10%，可忽略入射天顶角的变化对Kd 值 的影响。利用单一散射相函数模拟不同波长处的Kd 值，发现悬浮颗粒物浓度较高且以无机颗粒物为主体时，450～ 550nm 内模拟值大于实测值，550～650nm 内则相反，说明随着波长的增加，水体中颗粒物的散射呈现出前向散射 减弱、后向散射增强的变化趋势。在估算不同波长处的Kd 值时需考虑到散射相函数的光谱变化所造成的误差。 关键词：Monte Carlo, 数值模拟,漫衰减系数，入射天顶角，散射相函数 1.引 言 根据辐射传输理论，在光学性质均一的水体中，光的辐射强度随深度增加呈近似指数型衰减， 下行漫衰减系数（K ）表征了整个下行光场的分布情况，是一重要的表观光学量（AOP ）。无论 d 是计算水体的浑浊度、真光层深度，划分水体等级（Jerlov ，1976）；还是研究例如浮游植物光合 作用等一系列的生物、化学过程（Marra et al.，1995；Platt et al.，1988），Kd 值都有着不可或缺 的作用。海洋清洁水体的Kd 值的反演主要以经验算法为主，这种算法的劣势在于不能指出影响 Kd 值变化的因素，且在二类水体上反演精度较差。为了克服经验算法的局限性，许多学者做了 大量关于漫衰减系数（K ）与固有光学量（IOPs ）之间关系的定量研究，并构建了一系列半分析 d 估算模型 （Gordon，1989；Kirk ，1981，1984，1991；Lee et al.，2005 ；Liu et al.，2002 ）。由于 内陆湖泊水体的光学特性存在区域差异，半分析模型仍存在参数上的不确定性，移植性较差（Pan and Richard,2010；Wang et al.，2009 ）。国内主要集中研究太湖水体漫衰减系数在不同季节、不 同湖区以及不同深度上的漫衰减系数的差异 （黄昌春等，2009；张运林等，2003）。在计算 Kd 值时主要根据不同深度上的Ed 值进行推算，缺乏漫衰减系数数值模拟方面的研究。 Monte Carlo 方法采用一种较为灵活的方式给出辐射传输方程的数值解，它利用计算机模拟 大量个光子的运动轨迹，再将其转化为对应的辐射能量。在国外已经有大量研究利用Monte Carlo 方法来研究水下光场分布 （Kirk，1981,1984,1994,2003；Gordon,1975）,但研究对象主要是海 洋与近岸水体，有关内陆富营养化湖泊水体方面的研究较少。 太湖属于典型的内陆富营养浅水湖泊，由于风浪作用导致水体悬浮泥沙含量较高。悬浮泥沙 等无机悬浮颗粒物对光的衰减主要体现为散射，所以太湖水体的散射尤其是后向散射特性与一类 水体相比有较大差异 （乐成峰等，2009）。过去在使用辐射传输模型模拟一类或近海水体的漫衰 减系数时主要考虑到水体各光学组分的吸收和入射天顶角对Kd 值的影响，关于散射相函