Monte-Carlo方法在大气辐射中的应用.pptVIP

Monte-Carlo 方法  在大气辐射传输中的应用 张礼春 2007202139 主要内容： 一.Monte Carlo(蒙特卡罗)方法介绍 二.一个三维Monte-Carlo地气耦合辐射传输模式（3DMC） 三. Win X-ray软件的简单使用 一.Monte Carlo(蒙特卡罗)方法 1.Monte Carlo方法是计算数学中的一种计算方法，以概率与统计中的理论与方法为基础，以是否适于在计算机上使用为重要标志。Monte Carlo是摩纳哥的一个著名城市，以赌博闻名于世。 Monte Carlo方法借用这一城市的名称，属于象征性的。 2. Monte Carlo方法用于计算辐射传输时，即通过随机变量的统计试验求其近似解。 主要有两个特点： （1） Monte Carlo方法不受平面分层大气的限制，它可以解水平非均匀的问题，也可以解球面大气中的问题。 （2）Monte Carlo方法是模拟辐射在大气中传输的过程。它考察一个个光子在大气中传输时发生的散射、吸收，从统计角度来推算最终应有的辐射强度分布。 二.一个三维Monte-Carlo地气耦合辐射传输模式（3DMC） 1. 3DMC对太阳光子随机运动过程的模拟 模式中，太阳光子从大气上界进入后，可直接到达地表；也可与大气中的分子或粒子发生碰撞，然后能量完全消失或向上逸出大气层或到达地表；到达地表的光子可被地表完全吸收，也可被反射；经地表反射后的光子可直接逸出大气层；也可与大气中的分子或粒子发生一次或多次碰撞，然后能量完全消失或逸出大气层或重新到达地表。 判断光子在地气系统中传输时是发生散射、地表反射或逸出大气，要看光子前进一个随机位移后的位置。光子前进一个随机位移后，可能逸出大气层；如果仍在大气层内，则发生散射，此时是否存在吸收要看单次散射反照率是否为1；如果光子到达地表，当所在位置的反射率不为0时发生反射，否则被地表完全吸收。 2. 3DMC设计了均匀分布太阳入射能量的模拟方案 在水平非均一反射地表情况下，大气上界的视反射率表现为水平位置的函数，故模拟太阳光子入射时不能与均一地表（平面平行大气状态假设）情况下一样采用等效的定点入射方式，而必须模拟光子的随机水平位置入射方式。 由于实际太阳能量是水平均匀分布的（通常情况下可以这样假定），模拟实验证明，如果想当然地用两个随机数简单地分别定出入射光子的x坐标和y 坐标，就不能保证能量分布的均匀性，会得到错误的模拟结果。 3DMC根据概率分布原理设计了有效的随机入射方案。 假设已知光子入射范围的大小，这是一个大气上界以接收点为圆心的水平圆形区域，圆半径为R0，由于光子落在该圆内的概率为1，则落在圆内部任意半径R的同心圆中的概率P1与R和R0之间有如下关系： 其中，P1 为0~1之间的随机数。 光子的方位角为： 其中，P2为0~1之间的随机数。转化为直角坐标系后，光子释放位置（z取模拟大气上界高度）为 该方案保证了水平入射区域内辐射能量的均匀分布，符合实际物理机制。 3. 3DMC中便捷的光子三维空间定位方案 采用直角坐标系绝对位置定位方式, 方便模拟各种真实大气状况, 比如垂直分层大气和水平非均一大气。模式中地表目标物、太阳和探测器的位置关系见图1 , 图中地表目标物中心位置为A (0 , 0 , 0) ; θS为太阳天顶角, 太阳方位角取0 , θV为探测器观测天顶角, 为探测器与太阳的相对方位角。 图1 　3DMC 中太阳、探测器和地表目标物的几何方位示意图 3DMC 对光子的定位方法是在模式的初始化部分, 即读入大气消光系数的垂直分布廓线后,将其转换为大气光学厚度的垂直分布廓线, 同时采用光学厚度来度量光子运动步长。设Δτ为斜程光学厚度步长,Δτz 为垂直方向光学厚度步长, P为0～1 之间的随机数, θ为光子发生位移之前在前进方向的天顶角。则: Δτ= - ln P ; Δτz =Δτ cos θ。这样可以便捷地随时定出光子在地气系统中的空间位置。 其中, 光子散射方向模拟方法如下: 光子散射后的方位角是随机的, 而散射角与相函数有关。设相函数为P (θ) , P为随机数, 则散射角θ由下式确定: 4. 3