详情

原文传递一种海洋荧光探测半解析蒙特卡罗模拟方法

专利名称：	一种海洋荧光探测半解析蒙特卡罗模拟方法
摘要：	本发明提供了一种海洋荧光探测半解析蒙特卡罗模拟方法，该方法基于荧光物质垂向密度分布，利用蒙特卡罗方法生成荧光物质光源，模拟荧光光子在海水随机介质中的传输，模型考虑海水和水中悬浮物质的散射和吸收特性；结合半解析方法，快速计算每次散射过程中水上探测器感应的能量，统计每个光子的光程并形成时域探测信号。本发明可用于分析不同水体的水色信息，也有助于开展海洋光学散射特性的研究，节省了大量成本，从而为遥感探测和抗干扰算法设计、参数选择、性能评估提供参考。
专利类型：	发明专利
国家地区组织代码：	上海;31
申请人：	上海无线电设备研究所
发明人：	朱凌轩;王晓冰;王彪;霍熠炜;张元
专利状态：	有效
申请日期：	2018-11-22T00:00:00+0800
发布日期：	2019-06-07T00:00:00+0800
申请号：	CN201811398218.1
公开号：	CN109856092A
代理机构：	上海元好知识产权代理有限公司
代理人：	徐雯琼;张妍
分类号：	G01N21/64(2006.01);G;G01;G01N;G01N21
申请人地址：	200090 上海市杨浦区黎平路203号
主权项：	1.一种海洋荧光探测半解析蒙特卡罗模拟方法，用于实现荧光物质垂向分布不均匀的水体受光激发的荧光能量传输特性仿真，其特征在于，所述方法，包含以下过程：步骤1：将海水划分为深度等间距的光子多层传输媒介，归一化海水垂向深度的荧光物质分布谱，并将连续分布谱转化为等间距深度采样的离散概率分布；步骤2：统计实例海域中荧光物质的平均浓度，依据荧光物质浓度与激发光波长的经验公式关系计算受光激发的荧光峰波长；基于该荧光峰波长，利用Mie理论球形粒子近似，计算海水各层媒介的固有光学特性，包括散射系数、吸收系数、不对称因子及折射率；步骤3：基于荧光物质离散概率分布，利用蒙特卡罗方法随机生成荧光光源；考虑水下悬浮粒子及水体的散射、吸收效应，不同折射率的媒介层直接的反射、透射效应，模拟荧光光子的传输过程；结合解析法，计算每次散射后，荧光光子被水上探测器感应的能量；步骤4：基于光子输运过程，进行统计运算，转化为其他所需的传输信息；输出探测器收集的荧光光子能量，生成时域探测信号。 2.如权利要求1所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：步骤1中，荧光物质分布谱为等深度间距采样的离散概率分布谱；其中每单位间距水深的海水视为一层光子传输媒介，在层内的固有光学特性包含荧光物质浓度C、吸收系数σa、散射系数σb不随着深度变化。 3.如权利要求2所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：步骤2中，海水被处理为多层光子传输媒介，其固有光学特性依据荧光峰波长计算。 4.如权利要求3所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：步骤2：计算荧光物质受太阳光激发的荧光峰，根据Gitelson的经验回归公式得到荧光峰 λpeak＝a+b·Cchl (1) 其中Cchl为该海域荧光物质平均浓度，a和b为经验参数，依据荧光物质不同而变化。 5.如权利要求4所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：荧光物质来源于海洋中的叶绿素时，a在661.7-666.4nm之间变化，b在0.001-0.03之间变化。 6.如权利要求3所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：步骤3中，蒙特卡罗方法模拟每次光子散射后，水上探测器直接感应的能量的半解析计算公式为其中Ej为第j次散射时光子剩余能量，初始能量被设定为1，f(α)为Henyey-Greenstein相函数，为光子坐标指向探测器坐标的归一化向量，μti为第i层媒介的消光系数，Li为第i层的传输路径长度，ΔS为探测器感应面积，dj为探测器到光子的垂直距离，Ti,i-1(θ)为基于Fresnel方程的从i到第i-1层的光子透过率。 7.如权利要求3所述的海洋荧光探测半解析蒙特卡罗模拟方法，其特征在于：步骤3进一步包括以下步骤：步骤3.1：荧光光源位置的抽样；通过蒙特卡罗直接抽样法，生成[0,1]区间均匀分布随机数ξ1，判断随机数落在的概率区间其中Pi为第i层海水传输媒介的荧光物质分布概率密度，ZI为第I层海水的上边界深度；假设光源在层内服从均匀分布，直接抽样得到光子的出射位置(x,y,z)如下其中ξi是[0,1]区间均匀分布的随机数；步骤3.2：荧光光子传输方向的抽样；光子的传输方向由方位角和偏向角θ表示，在各向同性介质中通过蒙特卡罗采样得到转化为三维向量形式步骤3.3：解析法计算光子直接被探测的能量；探测器位于水平面上方M处(xm,ym,zm)，忽然多次折射效应，荧光光子不经过散射直接被探测的能量为其中Ej为第j次散射时光子剩余能量，初始能量被设定为1，f(α)为Henyey-Greenstein相函数，为光子坐标指向探测器坐标的归一化向量，μti为第i层的消光系数，Li为第i层的传输路径长度，ΔS为探测器感应面积，dj为探测器到光子的垂直距离，Ti,i-1(θ)为基于Fresnel方程的从i到第i-1层的光子透过率；统计完直接探测的能量后进入步骤3.4；步骤3.4：生成光子传输步长；首先判断光子是否有剩余步长Sleft，如果有则将剩余步长赋值给光子传输步长，否则根据光子在介质中的自由程概率密度函数计算，可以得到进一步，判断在现有的光子传输方向下，以该传输步长移动是否会碰撞媒介层边界；如果光子不会碰到边界进入步骤步骤3.5，否则进入步骤3.7；步骤3.5：如果光子不会碰到边界，则以该步长和方向更新光子的坐标步骤3.6：光子吸收和散射；光子经过传输步长到一个新位置后，光子的能量会被吸收，剩余光子能量为同时，光子的传输方向受到散射点影响也会改变；其中方位角在[0,2π]中均匀分布，偏转角服从Henyey-Greenstein相函数概率分布，重复步骤3.2更新光子的传输方向，并进入步骤3.10；步骤3.7：如果光子会碰到边界，则根据Fresnel方程判断光子在边界会发生反射或是透射过程：其中，θi和θt分别为光子的入射角和透射角；透射角可以由Snell公式计算；取[0,1]区间均匀分布的随机数ξ，如果ξ≤R(θ)，光子发生反射，否则发生透射；然后，根据Snell公式更新光子运动方向和光子所在媒介层；如果发生透射进入步骤3.8，发生反射则进入步骤3.9；步骤3.8：判断光子所在位置是否位于多层媒介的最外层，如果是则进入步骤3.11，否则进入步骤3.9 步骤3.9：将光子依据传输方向移动步长s'到该层媒介的边界，并计算剩余步长 sleft＝s-s' (12) 步骤3.10：判断光子是否死亡；如果光子死亡则进入步骤3.11，否则进一步判断光子剩余能量是否大于阈值；如果光子能量大于阈值则进入步骤3.3，否则进行Russian roulette轮盘判定光子是否继续输运过程，如果继续光子输运则进入步骤3.3，否则进入步骤3.11；步骤3.11：判断现在处理的光子是否为程序设定的最后一个光子，如果还有光子剩余，则重复步骤3.1-1.11，否则程序结束。
所属类别：	发明专利