一种内混合气溶胶光散射特性的计算方法技术

本发明专利技术涉及一种内混合气溶胶光散射特性的计算方法，其特征在于：包括如下步骤：1)内混合气溶胶单次散射特性计算：根据内混合气溶胶模型，计算各类型气溶胶粒子单次散射特性；2)光子状态初始化：确定光子首次发生散射的位置P1点的坐标；3)光子位置和光能量传输跟踪；4)光子接收判断；5)计算内混合气溶胶的各次散射特性。本发明专利技术的内混合气溶胶光散射特性的计算方法，可根据不同成分气溶胶组成的不同内混合模型，应用蒙特卡洛法对光在内混合气溶胶中的传输过程进行模拟仿真，并进一步考虑了内混合气溶胶粒子对光的散射和吸收作用，即可得到不同成分混合比例下组成的内混合气溶胶的各次光散射特性。

A Calculating Method for Light Scattering Characteristics of Internal Mixed Aerosols

The present invention relates to a method for calculating the light scattering characteristics of internal mixed aerosols, which includes the following steps: 1) calculating the single scattering characteristics of internal mixed aerosols: calculating the single scattering characteristics of various aerosol particles according to the internal mixed aerosol model; 2) initializing the photon state: determining the first scattering of photons. Position P1 coordinates; 3) photon position and light energy transmission tracking; 4) photon reception judgment; 5) calculation of scattering characteristics of mixed aerosols. The method for calculating the light scattering characteristics of internal mixed aerosols in the present invention can simulate the transmission process of light in internal mixed aerosols by Monte Carlo method according to different internal mixing models of aerosols with different components, and further consider the scattering and absorption of light by internal mixed aerosol particles. The light scattering characteristics of the internal mixed aerosols with different mixing ratios were obtained.

【技术实现步骤摘要】
一种内混合气溶胶光散射特性的计算方法
本专利技术涉及大气光学和激光传输领域的气溶胶光散射特性计算方法，特别是一种内混合气溶胶光散射特性的计算方法。
技术介绍
随着中国雾霾天气的越来越多，尤其是在我国北方供暖期，雾霾天气较为严重。雾霾主要是由矿物燃料和化石燃料的燃烧引起的，由此而产生的主要气溶胶成分是硫酸盐和黑碳。大量的研究指出，硫酸盐气溶胶是吸湿性粒子，其往往与大气中的水结合形成均匀混合的内混合粒子；而黑碳气溶胶是非吸湿性的粒子，但是随着时间慢慢老化，往往会作为核心部分与硫酸盐等水溶性气溶胶形成分层球形的内混合粒子。此时黑碳与硫酸盐气溶胶粒子的散射性质就会发生改变，气溶胶粒子的散射特性是影响光传输的关键因子。当光在大气中传输时，会受到上述内混合气溶胶粒子散射特性的影响，尤其是当大气中气溶胶粒子密度较大时，内混合气溶胶粒子的多次散射特性不可忽略，多次散射特性也是影响光传输的关键因子。近年来，国内外相关学者从不同的角度采用不同的方法对气溶胶多次散射特性进行了研究。对于多次散射特性，人们已经提出很多研究方法，如逐次迭代方法，以及建立于辐射传输方程的离散坐标法、球谐函数法、倍加累加法和蒙特卡洛法等。蒙特卡洛法作为处理光传输问题的常用方法，被广泛应用于光在多种介质的传输问题中，只要光子选取的足够多，便能够准确地模拟光子的多次散射问题。目前应用蒙特卡洛法研究内混合气溶胶多次散射特性时，为了方便起见，常假定其为单一均匀的粒子，不区分内混合的模型，就会产生较大的误差。因此，本专利技术结合实际大气中气溶胶内混合的模型，考虑单个内混合气溶胶粒子的散射特性，应用蒙特卡洛法对光在内混合气溶胶中的传输过程进行模拟仿真，即可得到内混合气溶胶的光散射特性，该方法将方便所有研究内混合气溶胶散射特性的学者进行应用，并对研究大气传输特性具有十分重要的现实意义。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种在考虑单个内混合气溶胶粒子散射特性基础上，计算内混合气溶胶光散射特性的方法，该方法将方便研究光在内混合气溶胶中的传输特性。本专利技术解决其技术问题是通过以下技术方案实现的：1)内混合气溶胶单次散射特性计算：根据内混合气溶胶模型，计算各类型气溶胶粒子单次散射特性：内混合气溶胶模型包括均匀混合模型及分层球混合模型，对于均匀混合模型，先计算粒子混合后等效折射率，再应用均匀球米散射理论计算其散射特性；对于分层球混合模型应用分层球米散射理论计算该内混合气溶胶粒子的非对称因子和单次散射反照率；2)光子状态初始化：光子从发射机发出，初始偏转角θ0在光束发散角2θ1内均匀发射，cosθ0在[cosθ1,1]之间均匀分布，则初始偏转角θ0可以抽样表示为：θ0＝arccos[1-r(1-cosθ1)](1)式中r为[0,1]区间内均匀分布的随机数；初始方位角在[0,2π]间均匀分布，可以抽样表示为：光子的初始传输方向D0由θ0和决定，可表示为光子的随机运动步长lm可根据比尔定律抽样表示为：式中σ为大气的消光系数,因此光子首次发生散射的位置P1点的坐标可表示为：P1(x1,y1,z1)＝O(0,0,0)+lmD0(5)3)光子位置和光能量传输跟踪：选取HG散射相函数确定每次碰撞之后的散射角，HG散射相函数的表达式为：式中g为光子所碰到的内混合气溶胶粒子的非对称因子，对散射相函数进行抽样即可得到散射角θs的表达式，当g≠0时：当g＝0时：θs＝arccos(2r-1)(8)同样方位角在[0,2π]间均匀分布，可抽样表示为：假设光子第m次碰撞时的坐标位置已知，为Pm(xm,ym,zm)，光子从Pm点散射到Pm+1的方向余弦为：则光子第m次碰撞后的迁移方向为：光子在某点与气溶胶粒子碰撞，在其第m次碰撞后的光子位置坐标为：Pm+1(xm+1,ym+1,zm+1)＝Pm(xm,ym,zm)+lDm′(12)光与气溶胶粒子碰撞之后，有一部分能量被气溶胶吸收，则第m次碰撞后的光能量的改变为：Em+1＝wEm(13)式中w为内混合气溶胶粒子的单次散射反照率，假设每一个初始光子对光能量的贡献相同都为a，即E0＝a；4)光子接收判断：满足以下两个判断条件，则认为光子被接收：光子恰好到达接收探测圆面，即光子从Pm点散射到Pm+1点，Pm+1点恰好在接收机探测圆面上；或光子的迁移轨迹穿越探测圆面，即光子的迁移轨迹PmPm+1与接收机探测截面有交点；当光子被接收，将停止对该光子的模拟跟踪，进行下一个光子的模拟，直到模拟完所有发射的N个光子；5)计算各次散射特性：对接收到的光能量进行统计，得到光在内混合气溶胶中传输时的各次散射光能量占总接收光能量的比值：直接透射光能量占总接收光能量的比值：一次散射光能量占总接收光能量的比值：二次散射光能量占总接收光能量的比值：三次散射光能量占总接收光能量的比值：式中I0为没有与气溶胶粒子发生碰撞的光子对光能量的贡献，即所有直接透射光子的贡献E0之和；I1为与气溶胶粒子碰撞一次的光子对光能量的贡献，即经过一次散射光子的贡献E1之和；I2为与气溶胶粒子碰撞两次的光子对光能量的贡献，即经过两次散射光子的贡献E2之和；I3为与气溶胶粒子碰撞三次的光子对光能量的贡献，即经过三次散射光子的贡献E3之和；以此类推，可得到各次散射光能量占总接收光能量的比值。本专利技术的优点和有益效果为：本专利技术的内混合气溶胶光散射特性的计算方法，可以根据不同类型气溶胶组成的不同内混合模型，计算不同混合比例下组成的内混合气溶胶的各次散射特性，而且应用蒙特卡洛法对光子的传输过程进行模拟，能够很好的揭示光子的多次散射现象，并进一步考虑了内混合气溶胶粒子对光的散射和吸收作用，可得到不同成分混合比例下组成的内混合气溶胶的各次光散射特性，该计算方法将方便所有研究气溶胶散射特性的学者进行应用，并且对研究实际大气传输特性具有十分重要的现实意义。附图说明图1(a)是内混合气溶胶模型：硫酸盐与水组成的均匀混合模型；图1(b)硫酸盐与黑碳组成的分层球内混合模型；图2是本专利技术的计算方法的流程图；图3是模拟计算的黑碳与硫酸盐组成的分层球内混合气溶胶粒子的各次散射特性的仿真结果图。具体实施方式下面通过具体实施例对本专利技术作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本专利技术的保护范围。一种内混合气溶胶光散射特性的计算方法，其特征在于：包括如下步骤：1)内混合气溶胶单次散射特性计算：根据内混合气溶胶的混合模型，分别应用不同方法计算该气溶胶粒子散射特性；所述内混合气溶胶模型包括均匀混合模型及分层球混合模型，对于均匀混合模型，先计算粒子混合后等效折射率，再应用均匀球米散射理论计算其散射特性；对于分层球混合模型应用分层球米散射理论计算该内混合气溶胶粒子的非对称因子和单次散射反照率。如图1所示，为内混合气溶胶的混合模型，其中图1(a)是硫酸盐和水组成的均匀混合模型，图1(b)是硫酸盐与黑碳组成的分层球混合模型。对于均匀混合模型，通常先计算粒子混合后的等效折射率，然后应用均匀球米散射理论计算其散射特性；对于分层球混合模型，根据内外混合半径的比值，应用分层球米散射理论计算其散射特性。本专利技术以硫酸盐与黑碳组成的分层球内混合模型为例，其中内核黑碳气溶胶粒子的半径假设为a，内混合气溶胶粒子的总半径为b，根据两种本文档来自技高网...

【技术保护点】
1.一种内混合气溶胶光散射特性的计算方法，其特征在于：包括如下步骤：1)内混合气溶胶单次散射特性计算：根据内混合气溶胶模型，计算各类型气溶胶粒子单次散射特性：内混合气溶胶模型包括均匀混合模型及分层球混合模型，对于均匀混合模型，先计算粒子混合后等效折射率，再应用均匀球米散射理论计算其散射特性；对于分层球混合模型应用分层球米散射理论计算该内混合气溶胶粒子的非对称因子和单次散射反照率；2)光子状态初始化：光子从发射机发出，初始偏转角θ0在光束发散角2θ1内均匀发射，cosθ0在[cosθ1,1]之间均匀分布，则初始偏转角θ0可以抽样表示为：θ0＝arccos[1‑r(1‑cosθ1)]  (1)式中r为[0,1]区间内均匀分布的随机数；初始方位角

【技术特征摘要】
1.一种内混合气溶胶光散射特性的计算方法，其特征在于：包括如下步骤：1)内混合气溶胶单次散射特性计算：根据内混合气溶胶模型，计算各类型气溶胶粒子单次散射特性：内混合气溶胶模型包括均匀混合模型及分层球混合模型，对于均匀混合模型，先计算粒子混合后等效折射率，再应用均匀球米散射理论计算其散射特性；对于分层球混合模型应用分层球米散射理论计算该内混合气溶胶粒子的非对称因子和单次散射反照率；2)光子状态初始化：光子从发射机发出，初始偏转角θ0在光束发散角2θ1内均匀发射，cosθ0在[cosθ1,1]之间均匀分布，则初始偏转角θ0可以抽样表示为：θ0＝arccos[1-r(1-cosθ1)](1)式中r为[0,1]区间内均匀分布的随机数；初始方位角在[0,2π]间均匀分布，可以抽样表示为：光子的初始传输方向D0由θ0和决定，可表示为光子的随机运动步长lm可根据比尔定律抽样表示为：式中σ为大气的消光系数,因此光子首次发生散射的位置P1点的坐标可表示为：P1(x1,y1,z1)＝O(0,0,0)+lmD0(5)3)光子位置和光能量传输跟踪：选取HG散射相函数确定每次碰撞之后的散射角，HG散射相函数的表达式为：式中g为光子所碰到的内混合气溶胶粒子的非对称因子，对散射相函数进行抽样即可得到散射角θs的表达式，当g≠0时：当g＝0时：θs＝arccos(2r-1)(8)同样方位角在[0,2π]间均匀分布，可抽样表示为：假设光子第m次碰撞时的坐标位置已知，为Pm(xm,ym,zm)，光子从Pm点散射到Pm+1的方向余弦为：则光子第m次碰撞...

【专利技术属性】
技术研发人员：马愈昭，刘文荣，高慧靓，熊兴隆，
申请(专利权)人：中国民航大学，
类型：发明
国别省市：天津,12