漫衰减系数常用于计算水体浑浊度和划分水体等级,是水中生物进行光合作用等生态过程中的重要参数。本文假设水体固有光学特性(IOPs)不随垂直剖面变化,水体表面为镜面。通过Monte Carlo方法模拟追踪光子在水中的运动过程,统计大量光子在不同深度上的能量分布,计算出下行漫衰减系数Kd值。结果表明:在450 nm,550 nm和650 nm波长处模拟的Kd值基本上小于实测Kd值,对应后向散射概率为3%的散射相函数模拟的相对误差最小,适合用于构建太湖水体的Kd值反演模型。进一步分析入射天顶角、b/a值对Kd值变化的影响,结果发现:在b/a较大的情况下,散射相函数成为影响Kd值变化的主要因素。当b/a值大于9时,入射天顶角从10°变化到80°时,Kd值的变化率小于10%,可忽略入射天顶角的变化对Kd值的影响。利用单一散射相函数模拟不同波长处的Kd值,发现悬浮颗粒物浓度较高且以无机颗粒物为主体时,450—550 nm内模拟值大于实测值,550—650 nm内则相反,说明随着波长的增加,水体中颗粒物的散射呈现出前向散射减弱、后向散射增强的变化趋势。在反演不同波长处的Kd值时需考虑到散射相函数的光谱变化所造成的误差。