CN111537425B - 生物气溶胶强消光特性颗粒形态的快速计算方法 - Google Patents
生物气溶胶强消光特性颗粒形态的快速计算方法 Download PDFInfo
- Publication number
- CN111537425B CN111537425B CN202010654585.4A CN202010654585A CN111537425B CN 111537425 B CN111537425 B CN 111537425B CN 202010654585 A CN202010654585 A CN 202010654585A CN 111537425 B CN111537425 B CN 111537425B
- Authority
- CN
- China
- Prior art keywords
- bioaerosol
- photon
- calculating
- target transmittance
- particle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims abstract description 84
- 230000008033 biological extinction Effects 0.000 title claims abstract description 49
- 239000000443 aerosol Substances 0.000 title abstract description 16
- 238000004364 calculation method Methods 0.000 title description 8
- 238000002834 transmittance Methods 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000015271 coagulation Effects 0.000 claims description 9
- 238000005345 coagulation Methods 0.000 claims description 9
- 244000005700 microbiome Species 0.000 claims description 9
- 230000009191 jumping Effects 0.000 claims description 6
- 239000002028 Biomass Substances 0.000 abstract description 2
- 241000228212 Aspergillus Species 0.000 description 16
- 230000000694 effects Effects 0.000 description 3
- 230000005670 electromagnetic radiation Effects 0.000 description 2
- 230000000813 microbial effect Effects 0.000 description 2
- 241001052560 Thallis Species 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
- G01N21/53—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
- G01N21/534—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke by measuring transmission alone, i.e. determining opacity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/49—Scattering, i.e. diffuse reflection within a body or fluid
- G01N21/53—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke
- G01N21/538—Scattering, i.e. diffuse reflection within a body or fluid within a flowing fluid, e.g. smoke for determining atmospheric attenuation and visibility
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N2015/0023—Investigating dispersion of liquids
- G01N2015/0026—Investigating dispersion of liquids in gas, e.g. fog
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1022—Measurement of deformation of individual particles by non-optical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/103—Particle shape
Landscapes
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出一种生物气溶胶强消光特性颗粒形态的快速计算方法,包括:获取生物气溶胶的目标透过率T和粒子群消光系数α;当所述生物气溶胶的目标透过率的最大值不小于设定阈值时,利用朗伯‑比尔定律推算生物气溶胶的实用化形态优选参数;小于设定阈值时,利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数;获取的生物气溶胶的实用化形态优选参数所对应的形态值,得到目标透过率条件下的生物气溶胶强消光的特性形态。本发明采用衰减物的最低质量作为形态实用化优选参数,可以对生物种质形态不同所带来的消光性能差异进行量化评判,指导生物气溶胶电磁衰减物的释放。
Description
技术领域
本发明属于光学技术领域,尤其涉及一种生物气溶胶强消光特性颗粒形态的快速计算方法。
背景技术
在大气中,生物气溶胶是其一个重要组成部分。生物气溶胶主要由微生物菌体、微生物孢子以及花粉等有机生物颗粒组成,由自然界或人类生产、生活产生,对人类的活动产生重要影响。生物气溶胶在自然界的大量存在以及其对光的衰减作用引起了越来越多人的关注和研究。实验观察和研究表明,大部分生物气溶胶粒子并不是单一的,有球形、有杆状、有丝状还有凝聚体等。研究强散射条件下生物气溶胶粒子形态确定的快速计算方法有助于加深对生物气溶胶强消光特性的理解,同时有助于对生物气溶胶的开发和利用。
当前对生物气溶胶消光特性的计算通常是正向的,即在确定形态下,采用适用的理论去计算生物气溶胶的消光性能。但是对于一定消光性能要求下,生物气溶胶最佳形态的确定方法和计算方法尚未发现相关文献资料。
因此,现有技术需要改进。
发明内容
为解决上述技术问题,本发明提出了一种生物气溶胶强消光特性颗粒形态的快速计算方法。
基于本发明实施例的一个方面,公开一种生物气溶胶强消光特性颗粒形态的快速计算方法,包括:
获取生物气溶胶的目标透过率T和粒子群消光系数α;
当所述生物气溶胶的目标透过率的最大值不小于设定阈值时,利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin,所述实用化形态优选参数mmin中的mmin为生物气溶胶中干扰物的最低质量,用于对生物种质形态不同所带来的消光性能差异进行量化评判;
当所述生物气溶胶的目标透过率的最大值小于设定阈值时,利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin;
根据所述实用化形态优选参数mmin获取生物气溶胶所对应的形态值,从而得到目标透过率条件下的生物气溶胶强消光的特性形态;
所述利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin包括:
步骤S1,初始化光子,并计算光子的行进步长,更新光子的位置;
步骤S2,判断光子是否在微生物凝聚粒子群中;
步骤S3,如果光子不在微生物凝聚粒子群中,则更新目标透过率数据,并判断该光子是否为最后一个光子;
步骤S4,如果步骤S3中判断是最后一个光子,则算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果步骤S3中判断不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
步骤S5,如果步骤S2中判断光子在微生物凝聚粒子群中,则更新该位置的光子的权重;步骤S6,判断光子的权重是否小于设定权重阈值;
如果步骤S6中光子的权重小于设定权重阈值,则执行步骤S3,更新目标透过率数据,判断该光子是否为最后一个光子;
如果是最后一个光子,则执行步骤S4,算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;步骤S7,如果步骤中种光子权重大于设定权重阈值,则更新光子的行进方向,跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
通过仿真算法计算粒子群的浓度C,计算mmin=C·ρv,式中,ρ为生物粒子的密度,v为生物粒子的体积。
基于本发明的生物气溶胶强消光特性颗粒形态的快速计算方法的另一个实施例中,还包括:优化计算生物气溶胶的形态;
所述优化计算生物气溶胶的形态包括:
获取生物气溶胶电磁衰减应用背景下的种质及其对应的形态;
计算生物气溶胶各种质的形态在其透光率小于设定阈值时的种质衰减物的最低质量;
获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态。
基于本发明的生物气溶胶强消光特性颗粒形态的快速计算方法的另一个实施例中,所述生物气溶胶的目标透过率的设定阈值为36.78%。
基于本发明的生物气溶胶强消光特性颗粒形态的快速计算方法的另一个实施例中,所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件包括:
生物气溶胶粒子群的光学厚度f(λ)小于1或目标透过率T(λ)大于20%;
当所述生物气溶胶粒子群的光学厚度f(λ)小于1时,所述目标透过率T(λ)大于36.78%;所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件为目标透过率T(λ)≥36.78%。
基于本发明的生物气溶胶强消光特性颗粒形态的快速计算方法的另一个实施例中,所述利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin包括:
朗伯-比尔定律的算法公式为:
T(λ)=exp[-α(λ)cL];
式中,T(λ)为目标透过率,α为粒子群消光系数,单位是m2/g,c为粒子群的质量浓度,单位为g/m3;L为光程,单位为m;
在利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin时,T(λ)、α、L均为已知参数,根据公式反推c,c为粒子群的质量浓度,单位为g/m3,由于实际计算的空间体积V也是确定的,因此,可以得到mmin=cV。
基于本发明的生物气溶胶强消光特性颗粒形态的快速计算方法的另一个实施例中,所述获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态包括:
已知某生物气溶胶种质W具有n种形态,为了达到T≤Tmax,Tmax为最大设定目标透过率,所需抛洒的1、2、3、……、n种形态的W种质衰减物的最低质量分别记为m1,min,m2,min,m3,min,……,mn,min;
mmin=min(mi,min)(i=1,2,3,……,n);
式中,i形态即为种质W在该电磁衰减应用背景下的最优形态,mmin为实用化形态优选参数。
与现有技术相比,本发明具有如下优点:
采用本发明的生物气溶胶强消光特性颗粒形态的快速计算方法,采用衰减物的最低质量mmin作为形态实用化优选参数,可以对生物种质形态不同所带来的消光性能差异进行量化评判,并且mmin可以作为工程应用中很有用的一个参数,用来指导生物气溶胶电磁衰减物的释放,本发明有广阔的应用范围,能够快速确定电磁衰减背景下生物气溶胶衰减物的最优形态,指导生物衰减物制备过程中结构参数的设定。
附图说明
图1为本发明提出的生物气溶胶强消光特性颗粒形态的快速计算方法的一个实施例的流程图;
图2为本发明提出的生物气溶胶强消光特性颗粒形态的快速计算方法的一个实施例的利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数的流程图;
图3为本发明提出的生物气溶胶强消光特性颗粒形态的快速计算方法的一个实施例的生物气溶胶的粒子群消光系数。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例只是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图和实施例对本发明提供的一种生物气溶胶强消光特性颗粒形态的快速计算方法进行更详细地说明。
图1为本发明提出的生物气溶胶强消光特性颗粒形态的快速计算方法的一个实施例的流程图,如图1所示,所述生物气溶胶强消光特性颗粒形态的快速计算方法包括:
10,获取生物气溶胶的目标透过率T和粒子群消光系数α;
20,当所述生物气溶胶的目标透过率的最大值不小于设定阈值时,利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin;
30,当所述生物气溶胶的目标透过率的最大值小于设定阈值时,利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin;
40,根据所述实用化形态优选参数mmin获取生物气溶胶所对应的形态值,从而得到目标透过率条件下的生物气溶胶强消光的特性形态;
50,优化计算生物气溶胶的形态;
所述优化计算生物气溶胶的形态包括:
获取生物气溶胶电磁衰减应用背景下的种质及其对应的形态;
计算生物气溶胶各种质的形态在其透光率小于设定阈值时的种质衰减物的最低质量;
获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态。
所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件包括:
生物气溶胶粒子群的光学厚度f(λ)小于1或目标透过率T(λ)大于20%;
当所述生物气溶胶粒子群的光学厚度f(λ)小于1时,所述目标透过率T(λ)大于36.78%;
所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件为目标透过率T(λ)≥36.78%。
所述利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin包括:
朗伯-比尔定律的算法公式为:
朗伯-比尔定律的算法公式为:
T(λ)=exp[-α(λ)cL];
式中,T(λ)为目标透过率,α为粒子群消光系数,单位是m2/g,c为粒子群的质量浓度,单位为g/m3;L为光程,单位为m;
在利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin时,T(λ)、α、L均为已知参数,根据公式反推c,c为粒子群的质量浓度,单位为g/m3,由于实际计算的空间体积V也是确定的,因此,可以得到mmin=cV。
如图2所示,当所述生物气溶胶的目标透过率的最大值小于设定阈值Tmax<36.78%时,利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin包括:
步骤S1,初始化光子,并计算光子的行进步长,更新光子的位置;
步骤S2,判断光子是否在微生物凝聚粒子群中;
步骤S3,如果光子不在微生物凝聚粒子群中,则更新目标透过率数据,并判断该光子是否为最后一个光子;
步骤S4,如果步骤S3中判断是最后一个光子,则算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果步骤S3中判断不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
步骤S5,如果步骤S2中判断光子在微生物凝聚粒子群中,则更新该位置的光子的权重;
步骤S6,判断光子的权重是否小于设定权重阈值;
如果步骤S6中光子的权重小于设定权重阈值,则执行步骤S3,更新目标透过率数据,判断该光子是否为最后一个光子;
如果是最后一个光子,则执行步骤S4,算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;步骤S7,如果步骤中种光子权重大于设定权重阈值,则更新光子的行进方向,跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
通过仿真算法计算粒子群的浓度C,计算mmin=C·ρv,式中,ρ为生物粒子的密度,v为生物粒子的体积。
所述获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态包括:
已知某生物气溶胶种质W具有n种形态,为了达到T≤Tmax,Tmax为最大设定目标透过率,所需抛洒的1、2、3、……、n种形态的W种质衰减物的最低质量分别记为m1,min,m2,min,m3,min,……,mn,min;
mmin=min(mi,min)(i=1,2,3,……,n);
式中,i形态即为种质W在该电磁衰减应用背景下的最优形态,mmin为实用化形态优选参数。
具体的,本发明中强消光一般指气溶胶透过率低于10%,但不同应用场景指标会有所不同。不论“强消光”的指标是多少,专利里的计算方法都是通用的。
实施例1:
针对某生物气溶胶电磁衰减应用背景,达到“衰减有效”的最低标准是电磁波在体积为V的生物气溶胶粒子群中衰减,当入射电磁波行进光程L后,目标透过率不能超过Tmax。假设某生物气溶胶种质W具有以下形态:a球形;b球形凝聚粒子;c菌丝。那么,为了达到T≤Tmax,所需抛洒的a、b、c三种形态的W种质衰减物的最低质量分别记为ma,min,mb,min,mc,min;mmin=min(mi,min)(i=a,b,c)
上式成立,则i形态即为种质W在该电磁衰减应用背景下的最优形态,mmin称为实用化形态优选参数。
生物气溶胶粒子均匀弥散在一定空间内形成的粒子群,会对入射电磁波造成衰减。当生物气溶胶的密度比较低,或形成的粒子群较薄的情况下,只考虑光源的一次散射时,这种衰减符合朗伯-比尔定律:
T(λ)=exp[-α(λ)cL];
式中,T(λ)为目标透过率,α为粒子群消光系数,单位是m2/g,c为粒子群的质量浓度,单位为g/m3;L为光程,单位为m;
对于朗伯-比尔定律来说,粒子群的光学厚度f(λ)是十分重要的一个参数:
f(λ)=α(λ)c;
T(λ)=exp[-f(λ)];
图3为本发明提出的生物气溶胶强消光特性颗粒形态的快速计算方法的一个实施例的生物气溶胶的粒子群消光系数,如图2所示,以某曲霉种质为例,设现在有某曲霉种质的三种结构:a孢子、b凝聚粒子、c菌丝;在0.24-14μm波段某曲霉三种形态的粒子群消光系数见附图3。对本实施例有影响的波段有:紫外波段0.24-0.4μm、可见光波段0.39-0.76μm、中红外波段3-5μm、远红外波段8-14μm。根据附图3中数据,求某曲霉三种气溶胶形态在紫外、可见光、中红外和远红外波段的平均粒子群消光系数αλ1~λ2,如下表所示;
假设现在需要对紫外、可见光、中红外和远红外波段进行衰减,要求某曲霉种质的生物气溶胶在V=30m3的空间内,在L=4m的光程范围,对入射电磁波衰减后透过率不超过40%,则各形态生物气溶胶衰减物最少用量如下表;
若要求某曲霉种质的生物气溶胶对入射电磁波衰减后,透过率不超过10%,则各形态生物气溶胶衰减物最少用量如下表;
因此,在30m3的空间内,4m的光程范围,电磁辐射穿过某曲霉气溶胶目标透过率小于40%,紫外和可见光波段某曲霉的最优衰减形态是单个孢子,最少需要释放某曲霉孢子的质量分别是6.8909g和6.2860g;中红外和远红外波段某曲霉的最优衰减形态是菌丝体,最少需要释放某曲霉菌丝体的质量分别是4.1457g和12.0978g。若要求电磁辐射穿过某曲霉气溶胶目标透过率小于10%,紫外和可见光波段某曲霉的最优衰减形态是单个孢子,最少需要释放某曲霉孢子的质量分别是8.6583g和7.8983g;中红外和远红外波段某曲霉的最优衰减形态是菌丝体,最少需要释放某曲霉菌丝体的质量分别是5.2090g和15.2007g。
对于本领域技术人员而言,显然本发明实施例不限于上述示范性实施例的细节,而且在不背离本发明实施例的精神或基本特征的情况下,能够以其他的具体形式实现本发明实施例。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明实施例的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本发明实施例内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。此外,显然“包括”一词不排除其他单元或步骤,单数不排除复数。系统、装置或终端权利要求中陈述的多个单元、模块或装置也可以由同一个单元、模块或装置通过软件或者硬件来实现。第一,第二等词语用来表示名称,而并不表示任何特定的顺序。
最后应说明的是,以上实施方式仅用以说明本发明实施例的技术方案而非限制,尽管参照以上较佳实施方式对本发明实施例进行了详细说明,本领域的普通技术人员应当理解,可以对本发明实施例的技术方案进行修改或等同替换都不应脱离本发明实施例的技术方案的精神和范围。
Claims (6)
1.一种生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,包括:
获取生物气溶胶的目标透过率T和粒子群消光系数α;
当所述生物气溶胶的目标透过率的最大值不小于设定阈值时,利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin,所述实用化形态优选参数mmin中的mmin为生物气溶胶中干扰物的最低质量,用于对生物种质形态不同所带来的消光性能差异进行量化评判;
当所述生物气溶胶的目标透过率的最大值小于设定阈值时,利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin;
根据所述实用化形态优选参数mmin获取生物气溶胶所对应的形态值,从而得到目标透过率条件下的生物气溶胶强消光的特性形态;
所述利用Monte Carlo概率仿真推算生物气溶胶的实用化形态优选参数mmin包括:
步骤S1,初始化光子,并计算光子的行进步长,更新光子的位置;
步骤S2,判断光子是否在微生物凝聚粒子群中;
步骤S3,如果光子不在微生物凝聚粒子群中,则更新目标透过率数据,并判断该光子是否为最后一个光子;
步骤S4,如果步骤S3中判断是最后一个光子,则算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果步骤S3中判断不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
步骤S5,如果步骤S2中判断光子在微生物凝聚粒子群中,则更新该位置的光子的权重;
步骤S6,判断光子的权重是否小于设定权重阈值;
如果步骤S6中光子的权重小于设定权重阈值,则执行步骤S3,更新目标透过率数据,判断该光子是否为最后一个光子;
如果是最后一个光子,则执行步骤S4,算法结束,获取入射光通过生物气溶胶目标透过率,得到生物粒子的密度ρ;
如果不是最后一个光子,则跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
步骤S7,如果步骤中种光子权重大于设定权重阈值,则更新光子的行进方向,跳转到步骤S1,重新对光子进行初始化,并计算光子的行进步长;
通过仿真算法计算粒子群的浓度C,计算mmin=Cρv,式中,ρ为生物粒子的密度,v为生物粒子的体积。
2.根据权利要求1所述的生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,还包括:优化计算生物气溶胶的形态;
所述优化计算生物气溶胶的形态包括:
获取生物气溶胶电磁衰减应用背景下的种质及其对应的形态;
计算生物气溶胶各种质的形态在其透光率小于设定阈值时的种质衰减物的最低质量;
获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态。
3.根据权利要求1所述的生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,所述生物气溶胶的目标透过率的设定阈值为36.78%。
4.根据权利要求3所述的生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件包括:
生物气溶胶粒子群的光学厚度f(λ)小于1或目标透过率T(λ)大于20%;
当所述生物气溶胶粒子群的光学厚度f(λ)小于1时,所述目标透过率T(λ)大于36.78%;
所述朗伯-比尔定律推算生物气溶胶的实用化形态优选参数的条件为目标透过率T(λ)≥36.78%。
5.根据权利要求3所述的生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,所述利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin包括:
朗伯-比尔定律的算法公式为:
T(λ)=exp[-α(λ)cL];
式中,T(λ)为目标透过率,α为粒子群消光系数,单位是m2/g,c为粒子群的质量浓度,单位为g/m3;L为光程,单位为m;
在利用朗伯-比尔定律推算生物气溶胶的实用化形态优选参数mmin时,T(λ)、α、L均为已知参数,根据公式反推c,c为粒子群的质量浓度,单位为g/m3,由于实际计算的空间体积V也是确定的,因此,可以得到mmin=cV。
6.根据权利要求2所述的生物气溶胶强消光特性颗粒形态的快速计算方法,其特征在于,所述获取种质衰减物的质量最低的生物气溶胶的种质,则该生物气溶胶的种质结构即为该电磁衰减应用背景下的最优形态包括:
已知某生物气溶胶种质W具有n种形态,为了达到T≤Tmax,Tmax为最大设定目标透过率,所需抛洒的1、2、3、……、n种形态的W种质衰减物的最低质量分别记为m1,min,m2,min,m3,min,……,mn,min;
mmin=min(mi,min)(i=1,2,3,……,n);
式中,i形态即为种质W在该电磁衰减应用背景下的最优形态,mmin为实用化形态优选参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010654585.4A CN111537425B (zh) | 2020-07-09 | 2020-07-09 | 生物气溶胶强消光特性颗粒形态的快速计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010654585.4A CN111537425B (zh) | 2020-07-09 | 2020-07-09 | 生物气溶胶强消光特性颗粒形态的快速计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111537425A CN111537425A (zh) | 2020-08-14 |
CN111537425B true CN111537425B (zh) | 2020-11-10 |
Family
ID=71974738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010654585.4A Active CN111537425B (zh) | 2020-07-09 | 2020-07-09 | 生物气溶胶强消光特性颗粒形态的快速计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111537425B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111811998B (zh) * | 2020-09-01 | 2020-12-01 | 中国人民解放军国防科技大学 | 一种目标波段下强吸收性生物颗粒组分的确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6453060B1 (en) * | 1999-06-29 | 2002-09-17 | Tri Path Imaging, Inc. | Method and apparatus for deriving separate images from multiple chromogens in a branched image analysis system |
CN101788719A (zh) * | 2009-12-31 | 2010-07-28 | 华中科技大学 | 一种获取衰光率连续且小于偏振片消光比的方法及衰减器 |
CN105137229A (zh) * | 2015-08-19 | 2015-12-09 | 中国人民解放军电子工程学院 | 一种生物颗粒跨波段综合电磁衰减能力的快速判定方法 |
CN106323920A (zh) * | 2015-07-10 | 2017-01-11 | 中国科学院遥感与数字地球研究所 | 基于蒙特卡洛算法的气溶胶多次散射模拟方法及系统 |
WO2018039497A1 (en) * | 2016-08-26 | 2018-03-01 | Spheryx, Inc. | Holographic characterization using hu moments |
CN109614851A (zh) * | 2018-10-26 | 2019-04-12 | 上海市环境监测中心(上海长三角区域空气质量预测预报中心) | 大气污染颗粒物形态参数的图像检测方法 |
-
2020
- 2020-07-09 CN CN202010654585.4A patent/CN111537425B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6453060B1 (en) * | 1999-06-29 | 2002-09-17 | Tri Path Imaging, Inc. | Method and apparatus for deriving separate images from multiple chromogens in a branched image analysis system |
CN101788719A (zh) * | 2009-12-31 | 2010-07-28 | 华中科技大学 | 一种获取衰光率连续且小于偏振片消光比的方法及衰减器 |
CN106323920A (zh) * | 2015-07-10 | 2017-01-11 | 中国科学院遥感与数字地球研究所 | 基于蒙特卡洛算法的气溶胶多次散射模拟方法及系统 |
CN105137229A (zh) * | 2015-08-19 | 2015-12-09 | 中国人民解放军电子工程学院 | 一种生物颗粒跨波段综合电磁衰减能力的快速判定方法 |
WO2018039497A1 (en) * | 2016-08-26 | 2018-03-01 | Spheryx, Inc. | Holographic characterization using hu moments |
CN109614851A (zh) * | 2018-10-26 | 2019-04-12 | 上海市环境监测中心(上海长三角区域空气质量预测预报中心) | 大气污染颗粒物形态参数的图像检测方法 |
Non-Patent Citations (2)
Title |
---|
The impact of aerosols on polarized sky radiance: model development, validation, and applications;C. Emde等;《Atmospheric Chemistry and Physics》;20100618;第10卷;第383-396页 * |
抗红外烟幕中固体消光材料的研究进展;霸书红等;《含能材料》;20181231;第26卷(第4期);第364-372页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111537425A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111537425B (zh) | 生物气溶胶强消光特性颗粒形态的快速计算方法 | |
CN115169969B (zh) | 生态环境状况综合评价方法、装置、电子设备及存储介质 | |
CN107256341B (zh) | 一种生态环境健康质量评估方法 | |
Hebert et al. | DISEASE OUTBREAKS IN PLANT-VECTOR-VIRUS MODELS WITH VECTOR AGGREGATION AND DISPERSAL. | |
Carneseca et al. | Association between particulate matter air pollution and monthly inhalation and nebulization procedures in Ribeirão Preto, São Paulo State, Brazil | |
CN113222854A (zh) | 一种小波变换的ct图像去噪方法 | |
CN112161815A (zh) | 一种车辆路噪主观评价值预测方法 | |
Rodrigues et al. | Morphological characterization of active germoplasm bank fig tree accessions | |
CN115187688A (zh) | 基于大气光偏振正交盲分离的雾图重构方法及电子设备 | |
CN111999268B (zh) | 一种大气消光系数湿度订正方法 | |
Li et al. | Sources and processes of organic aerosol in non-refractory PM1 and PM2. 5 during foggy and haze episodes in an urban environment of the Yangtze River Delta, China | |
Bowler | Accounting for the effect of observation errors on verification of MOGREPS | |
Kastak et al. | Onset, growth, and recovery of in-air temporary threshold shift in a California sea lion (Zalophus californianus) | |
Naoe et al. | State of mixture of atmospheric submicrometer black carbon particles and its effect on particulate light absorption | |
CN115865708B (zh) | 一种基于sir-d模型的复杂社交网络信息处置方法 | |
Pouliot et al. | Integration of the Model of Emissions of Gases and Aerosols from Nature (MEGAN) into the CMAQ Modeling System | |
Majumder et al. | The effect of noise in an HIV infection model with cytotoxic T-lymphocyte impairment | |
CN116663433B (zh) | 一种基于差分吸收光谱仪的差分优化方法 | |
Fei et al. | Channel simulation analysis of seawater laser communication | |
Wang et al. | Effects of forest gaps on soil microbial diversity in a Cunninghamia lanceolata Stand after a Severe Ice Storm | |
CN117766049A (zh) | 基于归一化敏感度通道选择的大气co2比值反演方法 | |
CN110187499A (zh) | 一种基于神经网络的片上集成光功率衰减器的设计方法 | |
Maroso et al. | Fitness consequences and ancestry loss in the Apennine brown bear after a simulated genetic rescue intervention | |
CN113989748A (zh) | 基于频域卷积神经网络的火焰探测方法 | |
Cavazzoni et al. | Modeling the effect of diffuse light on canopy photosynthesis in controlled environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |