CN107345893A - 一种粒子散射相函数测量装置及测量方法 - Google Patents
一种粒子散射相函数测量装置及测量方法 Download PDFInfo
- Publication number
- CN107345893A CN107345893A CN201710607284.4A CN201710607284A CN107345893A CN 107345893 A CN107345893 A CN 107345893A CN 201710607284 A CN201710607284 A CN 201710607284A CN 107345893 A CN107345893 A CN 107345893A
- Authority
- CN
- China
- Prior art keywords
- mrow
- phase function
- particle
- light intensity
- measured
- 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.)
- Granted
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000002245 particle Substances 0.000 claims abstract description 162
- 238000012937 correction Methods 0.000 claims abstract description 24
- 238000002474 experimental method Methods 0.000 claims abstract description 17
- 239000000523 sample Substances 0.000 claims description 53
- 230000003068 static effect Effects 0.000 claims description 19
- 239000000725 suspension Substances 0.000 claims description 12
- OVBPIULPVIDEAO-LBPRGKRZSA-N folic acid Chemical compound C=1N=C2NC(N)=NC(=O)C2=NC=1CNC1=CC=C(C(=O)N[C@@H](CCC(O)=O)C(O)=O)C=C1 OVBPIULPVIDEAO-LBPRGKRZSA-N 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 2
- 239000011521 glass Substances 0.000 abstract description 8
- 239000004793 Polystyrene Substances 0.000 description 11
- 229920002223 polystyrene Polymers 0.000 description 11
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 10
- 239000000377 silicon dioxide Substances 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000004587 chromatography analysis Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 239000012153 distilled water Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 239000004005 microsphere Substances 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 239000002612 dispersion medium Substances 0.000 description 2
- 230000004313 glare Effects 0.000 description 2
- 239000008187 granular material Substances 0.000 description 2
- 239000002609 medium Substances 0.000 description 2
- 235000012239 silicon dioxide Nutrition 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 239000005356 container glass Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000007430 reference method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
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
-
- 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/51—Scattering, i.e. diffuse reflection within a body or fluid inside a container, e.g. in an ampoule
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (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)
Abstract
一种粒子散射相函数测量装置及测量方法，本发明涉及粒子散射相函数测量装置及测量方法。本发明的目的是为了解决现有方法将探测器放置在透明容器外部时玻璃容器表面反射及内部界面多重反射的影响引起的测量误差大以及现有测量方法不能在悬浮颗粒样品量少的条件下准确获得其散射相函数的问题。过程为：一、实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；二、得到不同散射角下标准粒子系散射相函数的修正系数；三、实验测量待测粒子系的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系的散射光强分布，得到待测粒子系的散射相函数。本发明用于粒子散射相函数测量。
Description
技术领域
本发明涉及粒子散射相函数的测量装置及测量方法。
背景技术
粒子的散射相函数与吸收系数、散射系数等构成颗粒弥散介质的基本辐射特性参数，在大气遥感、海洋光学、燃烧过程模拟、光生物反应器、太阳能利用等涉及颗粒弥散介质辐射传输的领域均有重要应用。粒子散射相函数的测量有间接测量和直接测量两种方法。间接测量方法一般基于辐射传递正问题的求解结合反演算法获得散射相函数。这种方法需要事先假定散射相函数的参数化模型，一般不能得到颗粒的准确散射相函数。对于直接测量方法，一般通过周向布置的探测器来进行测量，通过探测器检测不同散射角度的散射光信号。传统的直接测量方法为浸入式测量，即将探测器浸入在液体内部，因为可以避免容器的界面反射干扰从而较准确的获得颗粒的散射相函数。但是该模式需要大量样品，不能在悬浮颗粒样品量少的条件下准确获得其散射相函数，同时探测器有受到样品污染的风险。虽然将探测器放置在透明容器外部时可不需要大的样品池，探测器也不会被污染，测试方便，但是会受到玻璃容器表面反射及内部界面多重反射的影响而引起很大测量误差。
发明内容
本发明的目的是为了解决现有方法将探测器放置在透明容器外部时玻璃容器表面反射及内部界面多重反射的影响引起的测量误差大以及现有测量方法不能在悬浮颗粒样品量少的条件下准确获得其散射相函数的问题，而提出一种粒子散射相函数测量装置及测量方法。
一种粒子散射相函数测量装置包括激光器(1)、斩光器(2)、分束器(3)、静态探测器(探测器1)(4)、动态探测器(探测器2)(5)、圆形比色皿(6)、平凸透镜(7)、前置放大器(8)和锁相放大器(9)；
激光器(1)发射激光，激光通过斩光器(2)到达分束器(3)，激光经分束器(3)后分为2束激光，2束光成90度角，其中一束激光法向照射装有粒子混悬液的圆形比色皿(6)，动态探测器(探测器2)(5)测量圆形比色皿(6)中粒子混悬液的粒子系光强信号；
动态探测器(探测器2)(5)与圆形比色皿(6)之间设置2个平凸透镜(7)；
动态探测器(探测器2)(5)与步进电机信号连接，步进电机控制动态探测器(5)(探测器2)移动，对圆形比色皿(6)中粒子混悬液的粒子系光强信号进行多角度测量；
另一束激光法向照射静态探测器(4)(探测器1)，静态探测器(4)测量激光的光强信号；
动态探测器(探测器2)(5)与静态探测器(4)分别连接前置放大器(8)和锁相放大器(9)，将测量得到的粒子系光强信号和激光的光强信号传递给前置放大器(8)和锁相放大器(9)。
一种粒子散射相函数测量方法具体过程为：
步骤一、实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；
步骤二、对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值，求解标准粒子系在不同散射角下的散射光强分布与标准粒子系散射相函数理论值的比值，得到不同散射角下标准粒子系散射相函数的修正系数；
步骤三、实验测量盛装在透明圆形比色皿中的待测粒子系在不同散射角下的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系在不同散射角下的散射光强分布，得到待测粒子系的散射相函数。
发明效果
由于本发明充分考虑了玻璃对测量产生的影响，采用分步计算的方式，即首先实验测量盛装在透明圆形比色皿中标准粒子的散射光强分布，然后用Lorenz-Mie模拟计算得到已知粒径分布的标准粒子系散射相函数理论值，标准粒子散射光强分布与理论散射相函数值对比得到不同散射角下标准粒子系散射相函数的修正系数，最后实验测量待测粒子的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子的散射光强分布，得到待测粒子的散射相函数，为一种非浸入式测量方法，装置简单，同时测量时需要的样品量少，降低了测量误差。
如图4是聚苯乙烯粒子系散射相函数实验测量值与Lorenz-Mie理论值随散射角的分布。聚苯乙烯混悬液(供应商：天津市倍思乐色谱技术开发中心)的基液为蒸馏水，单分散聚苯乙烯标准微球直径为0.5μm。聚苯乙烯在532nm波长下的折射指数是1.59。如图所示，修正并归一化后聚苯乙烯粒子系散射相函数实验测量值与Lorenz-Mie理论值两曲线吻合较好，进而说明该方法能够消除器皿的影响，提高粒子系散射相函数测量的精度。
附图说明
图1是具体实施方式一提出的粒子系散射相函数测量实验原理图；
图2是实施例提出的二氧化硅粒子系散射相函数图；
图3是实施例提出的散射相函数的修正系数图；
图4是实施例提出的聚苯乙烯粒子系散射相函数图；
图5是具体实施方式三提出的一种粒子散射相函数测量方法流程图。
具体实施方式
具体实施方式一：结合图1说明本实施方式，本实施方式的一种粒子散射相函数测量装置包括激光器(1)、斩光器(2)、分束器(3)、静态探测器(探测器1)(4)、动态探测器(探测器2)(5)、圆形比色皿(6)、平凸透镜(7)、前置放大器(8)和锁相放大器(9)；
激光器(1)发射激光(参考激光)，激光通过斩光器(2)到达分束器(3)，激光经分束器(3)后分为2束激光，2束光成90度角，其中一束激光法向照射装有粒子混悬液的圆形比色皿(6)，动态探测器(探测器2)(5)测量圆形比色皿(6)中粒子混悬液的粒子系光强信号；
动态探测器(探测器2)(5)与圆形比色皿(6)之间设置2个平凸透镜(7)；
动态探测器(探测器2)(5)与步进电机信号连接，步进电机控制动态探测器(5)(探测器2)移动，对圆形比色皿(6)中粒子混悬液的粒子系光强信号进行多角度测量，就可以获得一定波长在多角度下的散射光强分布的实验结果；
另一束激光法向照射静态探测器(4)(探测器1)，静态探测器(4)测量激光的光强信号；以静态探测器(探测器1)的测量信号作为参考光强信号，
动态探测器(探测器2)(5)与静态探测器(4)分别连接前置放大器(8)和锁相放大器(9)，将测量得到的粒子系光强信号和激光的光强信号传递给前置放大器(8)和锁相放大器(9)。
实验装置中探测器1和探测器2的测量信号接入前置放大器和锁相放大器，这可以显著提高对微弱散射信号的探测灵敏度。
具体实施方式二：本实施方式与具体实施方式一不同的是：所述2个平凸透镜(7)的凸面相对设置，两个平凸透镜的平面分别面向圆形比色皿(6)和动态探测器(探测器2)(5)设置。
该方式只能接收垂直探头方向的光线，从而可消除杂散光的影响。
先根据需要选取一定波长的激光器(或其它形式的单色光源)，激光器发射激光，激光通过斩光器到达分束器，激光经分束器后分为2束激光，其中一激光法向照射装有粒子混悬液的圆形比色皿池，然后通过步进电机调整动态探测器(探测器2)的角度进行多角度测量，就可以获得一定波长在多角度下的散射光强分布的实验结果；
另一束激光照射静态探测器(探测器1)，以静态探测器(探测器1)的测量信号作为参考光强信号，
探测器2与圆形比色皿池之间放置两个平凸透镜，且两个平凸透镜的凸相对设置，两个平凸透镜的平面分别面向圆形比色皿池和探测器2设置，(动态探测器(探测器2)的前面放置两个平凸透镜)，该方式只能接收垂直探头方向的光线，从而可消除杂散光的影响。实验装置中探测器1和探测器2的测量信号接入前置放大器和锁相放大器，这可以显著提高对微弱散射信号的探测灵敏度。
动态探测器的周向移动通过步进电机驱动，提供水平入射方向从0到180°的动态测量范围；
其它步骤及参数与具体实施方式一相同。
具体实施方式三：结合图5说明本实施方式，本实施方式的一种粒子散射相函数测量方法具体是按照以下步骤实施的：
步骤一、实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；
步骤二、对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值，求解标准粒子系在不同散射角下的散射光强分布与标准粒子系散射相函数理论值的比值，得到不同散射角下标准粒子系散射相函数的修正系数；
步骤三、实验测量盛装在透明圆形比色皿中的待测粒子系在不同散射角下的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系在不同散射角下的散射光强分布，得到待测粒子系的散射相函数；即完成了一种测量粒子系散射相函数的参比法。
本实施方式效果：
本实施方式充分考虑了透明玻璃器皿对于粒子系散射相函数测量的影响，弥补了现存方法中测量粒子系散射相函数的缺陷，提高了测量精度如图4。图4是修正后结果与Lorenz-Mie计算得到的散射相函数理论值进行对比，而本发明方法得到的数值结果与理论值吻合较好，消除了玻璃器皿对于散射相函数测量的影响。
由于本实施方式充分考虑了玻璃对测量产生的影响，采用分步计算的方式，即首先实验测量盛装在透明圆形比色皿中标准粒子的散射光强分布，然后用Lorenz-Mie模拟计算得到已知粒径分布的标准粒子系散射相函数理论值，标准粒子散射光强分布与理论散射相函数值对比得到不同散射角下标准粒子系散射相函数的修正系数，最后实验测量待测粒子的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子的散射光强分布，对修正后的待测粒子散射光强分布进行归一化处理，得到待测粒子的散射相函数。本发明是要解决现有测量方法测量时忽略玻璃器皿影响从而产生偏差的问题，提出一种简洁方便且又能获得高精度粒子系散射相函数的测量方法。
具体实施方式四：本实施方式与具体实施方式三不同的是：所述步骤一中实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；具体过程为：
实验测量粒子的散射光强分布时，由于散射相函数的数量级跨度较大(104～10-3)，所以当探测弱信号时需要用锁相放大器放大其强度，且需要在不同放大档位进行多次测量来提高其精度。其它步骤及参数与具体实施方式一相同。
如图1所示，探测器2与圆形比色皿池之间放置两个平凸透镜，且两个平凸透镜的凸相对设置，两个平凸透镜的平面分别面向圆形比色皿池和探测器2设置，计算并确定焦距位置后将它们放入透镜套管中，从而保证探测器2只能接收垂直方向O点的光线。
在单次散射条件下，待测粒子在O点的光通量表示为：
式中：q0是入射光通量；tcuv是比色皿容器玻璃的透射比；LOA是A点到O点的直线距离；A点为激光通过圆形比色皿池O点与比色皿池的交点，O点为圆形比色皿池和中心位置，βeff是介质等效衰减系数；在散射角为Θ时，探测器2探测到的待测粒子光通量为：
式中：Φ(Θ)是待测粒子散射角为Θ的散射相函数；σ是待测粒子介质的散射系数；LOB是B点到O点的直线距离；B点为比色皿池O点与探测器2连线在比色皿池上的交点；ΔΩ是探测器2相对于O点的立体角；G(Θ)是圆形比色皿皿池与介质相互影响的修正系数；公式(2)转化为：
以标准粒子作为参考时，标准粒子的散射相函数为：
式中：qD,Ref(Θ)为在散射角为Θ时，待测粒子探测器2探测到的待测粒子的光通量，σRef为待测粒子介质的散射系数，Ref为待测粒子；
把公式3和4相除可得：
式中：X为一个定值系数；从公式(5)可得，待测粒子和参考粒子的散射相函数比值大小只与不同散射角的探测器接收信号数值和X大小相关；由此可知，使用单分散标准粒子作为参考系，能消除玻璃器皿对粒子系散射相函数测量产生的影响。
测量系统中，散射强度分布的实验测量关系为：标准粒子系在不同散射角下的散射光强分布SEXP(Θ)为：
式中：qD，1(Θ)为静态探测器(探测器1)(4)测量标准粒子系得到的光强信号；qD，2(Θ)为动态探测器(探测器2)(5)探测器2(散射信号探测器)测量标准粒子系得到的光强信号；a为分束系数，当光路确定以后，该系数应为常数(当光路确定后，激光光源强度确定时，激光经分束器后，一条激光法向量直接入射到探测器1和另一条激光法向量入射到探测器2，探测器1与探测器2的比值为分束系数)；ΔΩ为动态探测器(探测器2)(5)相对于O点的立体角；O点为透明圆形比色皿的中心位置。
其它步骤及参数与具体实施方式三相同。
具体实施方式五：本实施方式与具体实施方式三或四不同的是：所述步骤二中对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值，求解标准粒子系在不同散射角下的散射光强分布与标准粒子系散射相函数理论值的比值，得到不同散射角下标准粒子系散射相函数的修正系数；具体过程为：
在单散射条件下，对标准粒子的粒径分布进行测定(显微镜拍照片，用ImageJ软件)，通过米散射理论计算标准粒子系散射相函数理论值ΦMie(Θ)，求解标准粒子系在不同散射角下的散射光强分布SEXP(Θ)与标准粒子系散射相函数理论值ΦMie(Θ)的比值，得到不同散射角下标准粒子系散射相函数的修正系数U(Θ)，公式为：
其它步骤及参数与具体实施方式三或四相同。
具体实施方式六：本实施方式与具体实施方式三至五之一不同的是：所述步骤三中实验测量盛装在透明圆形比色皿中的待测粒子系在不同散射角下的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系在不同散射角下的散射光强分布，得到待测粒子系的散射相函数；具体过程为：
待测粒子系在不同散射角下的散射光强分布SEXP(Θ)′为：
式中，qD,2(Θ)′为动态探测器(探测器2)(5)(散射信号探测器)测量待测粒子系得到的光强信号；qD,1(Θ)′为静态探测器(探测器1)(4)测量待测粒子系得到的光强信号；
利用得到的标准粒子系散射相函数的修正系数U(Θ)去修正待测粒子的散射光强分布SEXP(Θ)′，得到待测粒子系散射相函数，表示为：
对得到的待测粒子的散射相函数Φ(Θ)进行归一化处理，得到待测粒子的散射相函数。
步骤三中计算散射相函数修正系数时，参考粒子的粒径分布对其至关重要，因此需要精确给出。
其它步骤及参数与具体实施方式三至五之一相同。
采用以下实施例验证本发明的有益效果：
实施例一：
本实施例是一种粒子散射相函数测量方法具体是按照以下步骤实施的：
二氧化硅标准微球(供应商：天津市倍思乐色谱技术开发中心)被选为验证样本，其基液为蒸馏水。二氧化硅标准微球直径为0.5μm，且待测样本为单分散微球。盛放待测样本圆形比色皿池的玻璃壁厚为1.5mm。特别注意的是实验测量的时候光斑要控制的尽量小(光斑直径小于或者等于1mm)，使光斑入射圆形比色皿时近似于平面入射。激光器的波长为532nm，二氧化硅在此波长下的折射指数是1.46。
图2是二氧化硅粒子实验测量散射光强分布与Lorenz-Mie散射相函数理论值随散射角的分布。在散射角从0°变化至160°时散射光强从103数量级降低至10-2数量级，由于探测器本身的遮挡，在临近散射角180°的附近5°数据无法有效探测到。图3是散射相函数的修正系数。在每个散射角下的修正系数是通过实验测量的散射光强值除以Lorenz-Mie散射相函数理论值得到，但由于遮挡作用，临近散射角180°的附近5°数据没有给出。
图4是聚苯乙烯粒子系散射相函数实验测量值与Lorenz-Mie理论值随散射角的分布。聚苯乙烯混悬液(供应商：天津市倍思乐色谱技术开发中心)的基液为蒸馏水，单分散聚苯乙烯标准微球直径为0.5μm。聚苯乙烯在532nm波长下的折射指数是1.59。如图所示，修正并归一化后聚苯乙烯粒子系散射相函数实验测量值与Lorenz-Mie理论值两曲线吻合较好，进而说明该方法能够消除器皿的影响，提高粒子系散射相函数测量的精度。
以上实例的结果表明，本发明中的方法能够消除玻璃容器的影响，经修正后其结果更加准确有效。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (6)
1.一种粒子散射相函数测量装置，其特征在于：所述装置包括激光器(1)、斩光器(2)、分束器(3)、静态探测器(4)、动态探测器(5)、圆形比色皿(6)、平凸透镜(7)、前置放大器(8)和锁相放大器(9)；
激光器(1)发射激光，激光通过斩光器(2)到达分束器(3)，激光经分束器(3)后分为2束激光，2束光成90度角，其中一束激光法向照射装有粒子混悬液的圆形比色皿(6)，动态探测器(5)测量圆形比色皿(6)中粒子混悬液的粒子系光强信号；
动态探测器(5)与圆形比色皿(6)之间设置2个平凸透镜(7)；
动态探测器(5)与步进电机信号连接，步进电机控制动态探测器(5)移动，对圆形比色皿(6)中粒子混悬液的粒子系光强信号进行多角度测量；
另一束激光法向照射静态探测器(4)，静态探测器(4)测量激光的光强信号；
动态探测器(5)与静态探测器(4)分别连接前置放大器(8)和锁相放大器(9)，将测量得到的粒子系光强信号和激光的光强信号传递给前置放大器(8)和锁相放大器(9)。
2.根据权利要求1所述一种粒子散射相函数测量装置，其特征在于：所述2个平凸透镜(7)的凸面相对设置，两个平凸透镜的平面分别面向圆形比色皿(6)和动态探测器(5)设置。
3.一种粒子散射相函数测量方法，是基于一种粒子散射相函数测量装置实现的，其特征在于：一种粒子散射相函数测量方法具体过程为：
步骤一、实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；
步骤二、对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值，求解标准粒子系在不同散射角下的散射光强分布与标准粒子系散射相函数理论值的比值，得到不同散射角下标准粒子系散射相函数的修正系数；
步骤三、实验测量盛装在透明圆形比色皿中的待测粒子系在不同散射角下的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系在不同散射角下的散射光强分布，得到待测粒子系的散射相函数。
4.根据权利要求3所述一种粒子散射相函数测量方法，其特征在于：所述步骤一中实验测量盛装在透明圆形比色皿中的标准粒子系在不同散射角下的散射光强分布；具体过程为：
标准粒子系在不同散射角下的散射光强分布SEXP(Θ)为：
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mi>X</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>aq</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>&Omega;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中：qD，1(Θ)为静态探测器(4)测量标准粒子系得到的光强信号；qD，2(Θ)为动态探测器(5)测量标准粒子系得到的光强信号；a为分束系数；ΔΩ为动态探测器(5)相对于O点的立体角；O点为透明圆形比色皿的中心位置。
5.根据权利要求4所述一种粒子散射相函数测量方法，其特征在于：所述步骤二中对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值，求解标准粒子系在不同散射角下的散射光强分布与标准粒子系散射相函数理论值的比值，得到不同散射角下标准粒子系散射相函数的修正系数；具体过程为：
在单散射条件下，对标准粒子的粒径分布进行测定，通过米散射理论计算标准粒子系散射相函数理论值ΦMie(Θ)，求解标准粒子系在不同散射角下的散射光强分布SEXP(Θ)与标准粒子系散射相函数理论值ΦMie(Θ)的比值，得到不同散射角下标准粒子系散射相函数的修正系数U(Θ)，公式为：
<mrow>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mi>X</mi>
<mi>P</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>M</mi>
<mi>i</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
6.根据权利要求5所述一种粒子散射相函数测量方法，其特征在于：所述步骤三中实验测量盛装在透明圆形比色皿中的待测粒子系在不同散射角下的散射光强分布，利用得到的标准粒子系散射相函数的修正系数去修正待测粒子系在不同散射角下的散射光强分布，得到待测粒子系的散射相函数；具体过程为：
待测粒子系在不同散射角下的散射光强分布SEXP(Θ)′为：
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mi>X</mi>
<mi>P</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<msub>
<mi>aq</mi>
<mrow>
<mi>D</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>&prime;</mo>
</msup>
<mi>&Delta;</mi>
<mi>&Omega;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，qD,2(Θ)′为动态探测器(5)测量待测粒子系得到的光强信号；qD,1(Θ)′为静态探测器(4)测量待测粒子系得到的光强信号；
利用得到的标准粒子系散射相函数的修正系数U(Θ)去修正待测粒子的散射光强分布SEXP(Θ)′，得到待测粒子系散射相函数，表示为：
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>E</mi>
<mi>X</mi>
<mi>P</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>&Theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710607284.4A CN107345893B (zh) | 2017-07-24 | 2017-07-24 | 一种粒子散射相函数测量装置及测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710607284.4A CN107345893B (zh) | 2017-07-24 | 2017-07-24 | 一种粒子散射相函数测量装置及测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107345893A true CN107345893A (zh) | 2017-11-14 |
| CN107345893B CN107345893B (zh) | 2020-07-07 |
Family
ID=60257004
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710607284.4A Active CN107345893B (zh) | 2017-07-24 | 2017-07-24 | 一种粒子散射相函数测量装置及测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107345893B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108287126A (zh) * | 2018-03-23 | 2018-07-17 | 中国计量科学研究院 | 纳米颗粒粒径测量系统 |
| CN109073549A (zh) * | 2018-06-28 | 2018-12-21 | 深圳市汇顶科技股份有限公司 | 散射角测量装置和散射角测量方法 |
| CN110132874A (zh) * | 2019-05-27 | 2019-08-16 | 哈尔滨工业大学 | 基于多角度测量的弥散介质光学参数场探测装置及方法 |
| CN110618109A (zh) * | 2019-10-31 | 2019-12-27 | 中国科学院长春光学精密机械与物理研究所 | 一种用于液态油中杂质颗粒和游离水进行标定和测量的装置 |
| CN111239079A (zh) * | 2020-03-09 | 2020-06-05 | 上海交通大学 | 一种定光学深度的时变浑浊场模拟装置 |
| CN111650100A (zh) * | 2020-06-16 | 2020-09-11 | 辽东学院 | 一种基于Mie氏散射理论的颗粒粒度测量设备 |
| CN113155683A (zh) * | 2021-05-17 | 2021-07-23 | 北京市人工影响天气办公室 | 一种云粒子散射函数的探测系统、方法和装置 |
| CN114609093A (zh) * | 2022-03-24 | 2022-06-10 | 合肥工业大学 | 基于鱼眼镜头的成像浊度计的散射相函数的计算方法 |
| CN115032176A (zh) * | 2022-05-13 | 2022-09-09 | 中国科学院南海海洋研究所 | 一种广角体散射函数及衰减系数测量方法 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08101113A (ja) * | 1994-09-30 | 1996-04-16 | Shimadzu Corp | 雨滴粒度分布測定装置 |
| EP1464966A1 (en) * | 2003-04-03 | 2004-10-06 | C.S.E.M. Centre Suisse D'electronique Et De Microtechnique Sa | Apparatus and method for determining the concentration and motility of light scattering particles |
| CN101140222A (zh) * | 2007-10-19 | 2008-03-12 | 天津炜辐医疗科技有限公司 | 测定包括混浊介质材料全部光学参数的光谱仪系统及方法 |
| CN101535803A (zh) * | 2006-09-25 | 2009-09-16 | 兴和株式会社 | 凝胶化测定装置以及试样盒 |
| EP2169433A1 (en) * | 2008-09-26 | 2010-03-31 | Siemens Aktiengesellschaft | A method and a device for analyzing a forest enviroment by using scattered beams of electromagnetic radiation |
| CN102692367A (zh) * | 2012-06-14 | 2012-09-26 | 南京中医药大学 | 纳米粒子辨识系统装置及其识别方法 |
| CN202486052U (zh) * | 2011-12-15 | 2012-10-10 | 华中科技大学 | 一种三维体散射函数的测量系统 |
| CN103454203A (zh) * | 2013-09-09 | 2013-12-18 | 中国科学院合肥物质科学研究院 | 一种大气颗粒物粒径和化学成分的实时在线测量系统及测量方法 |
| CN104777080A (zh) * | 2015-04-10 | 2015-07-15 | 中国石油大学(北京) | 输气管道内颗粒与液滴的在线检测装置及方法 |
| CN106018346A (zh) * | 2016-05-09 | 2016-10-12 | 山东省科学院海洋仪器仪表研究所 | 气溶胶散射相函数观测装置及观测方法 |
| CN106018286A (zh) * | 2016-05-18 | 2016-10-12 | 哈尔滨工业大学 | 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法 |
| CN106383072A (zh) * | 2016-10-17 | 2017-02-08 | 哈尔滨工业大学 | 基于多角度光散射‑透射法的球形颗粒光学常数与粒径分布同时测量方法 |
| CN106644852A (zh) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学 | 基于超短脉冲激光辐照同时获取球形颗粒光学常数与粒径分布的测量方法 |
-
2017
- 2017-07-24 CN CN201710607284.4A patent/CN107345893B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08101113A (ja) * | 1994-09-30 | 1996-04-16 | Shimadzu Corp | 雨滴粒度分布測定装置 |
| EP1464966A1 (en) * | 2003-04-03 | 2004-10-06 | C.S.E.M. Centre Suisse D'electronique Et De Microtechnique Sa | Apparatus and method for determining the concentration and motility of light scattering particles |
| CN101535803A (zh) * | 2006-09-25 | 2009-09-16 | 兴和株式会社 | 凝胶化测定装置以及试样盒 |
| CN101140222A (zh) * | 2007-10-19 | 2008-03-12 | 天津炜辐医疗科技有限公司 | 测定包括混浊介质材料全部光学参数的光谱仪系统及方法 |
| EP2169433A1 (en) * | 2008-09-26 | 2010-03-31 | Siemens Aktiengesellschaft | A method and a device for analyzing a forest enviroment by using scattered beams of electromagnetic radiation |
| CN202486052U (zh) * | 2011-12-15 | 2012-10-10 | 华中科技大学 | 一种三维体散射函数的测量系统 |
| CN102692367A (zh) * | 2012-06-14 | 2012-09-26 | 南京中医药大学 | 纳米粒子辨识系统装置及其识别方法 |
| CN103454203A (zh) * | 2013-09-09 | 2013-12-18 | 中国科学院合肥物质科学研究院 | 一种大气颗粒物粒径和化学成分的实时在线测量系统及测量方法 |
| CN104777080A (zh) * | 2015-04-10 | 2015-07-15 | 中国石油大学(北京) | 输气管道内颗粒与液滴的在线检测装置及方法 |
| CN106018346A (zh) * | 2016-05-09 | 2016-10-12 | 山东省科学院海洋仪器仪表研究所 | 气溶胶散射相函数观测装置及观测方法 |
| CN106018286A (zh) * | 2016-05-18 | 2016-10-12 | 哈尔滨工业大学 | 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法 |
| CN106383072A (zh) * | 2016-10-17 | 2017-02-08 | 哈尔滨工业大学 | 基于多角度光散射‑透射法的球形颗粒光学常数与粒径分布同时测量方法 |
| CN106644852A (zh) * | 2016-10-17 | 2017-05-10 | 哈尔滨工业大学 | 基于超短脉冲激光辐照同时获取球形颗粒光学常数与粒径分布的测量方法 |
Non-Patent Citations (5)
| Title |
|---|
| (美)布朗（BRAUN,R.D.）著,北京大学化学系 等译: "《最新仪器分析技术全书》", 30 September 1990, 化学工业出版社 * |
| M. P. MENGOC ET AL.: "A STEP-FUNCTION APPROXIMATION FOR THE EXPERIMENTAL DETERMINATION OF THE EFFECTIVE SCATTERING PHASE FUNCTION OF PARTICLES", 《J.QUANT.SPECTROSC.RADIAT.TRANSFER》 * |
| 本书编委: "《辞海 数学•物理•化学分册》", 31 December 1987, 上海辞书出版社 * |
| 王保华: "《生物医学测量与仪器》", 28 February 2009, 复旦大学出版社 * |
| 齐利敏: "生物组织散射相函数实验的研究", 《河北师范大学学报(自然科学版)》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108287126A (zh) * | 2018-03-23 | 2018-07-17 | 中国计量科学研究院 | 纳米颗粒粒径测量系统 |
| CN109073549A (zh) * | 2018-06-28 | 2018-12-21 | 深圳市汇顶科技股份有限公司 | 散射角测量装置和散射角测量方法 |
| WO2020000319A1 (zh) * | 2018-06-28 | 2020-01-02 | 深圳市汇顶科技股份有限公司 | 散射角测量装置和散射角测量方法 |
| CN110132874B (zh) * | 2019-05-27 | 2021-11-19 | 哈尔滨工业大学 | 基于多角度测量的弥散介质光学参数场探测装置及方法 |
| CN110132874A (zh) * | 2019-05-27 | 2019-08-16 | 哈尔滨工业大学 | 基于多角度测量的弥散介质光学参数场探测装置及方法 |
| CN110618109A (zh) * | 2019-10-31 | 2019-12-27 | 中国科学院长春光学精密机械与物理研究所 | 一种用于液态油中杂质颗粒和游离水进行标定和测量的装置 |
| CN111239079A (zh) * | 2020-03-09 | 2020-06-05 | 上海交通大学 | 一种定光学深度的时变浑浊场模拟装置 |
| CN111650100A (zh) * | 2020-06-16 | 2020-09-11 | 辽东学院 | 一种基于Mie氏散射理论的颗粒粒度测量设备 |
| CN113155683A (zh) * | 2021-05-17 | 2021-07-23 | 北京市人工影响天气办公室 | 一种云粒子散射函数的探测系统、方法和装置 |
| CN113155683B (zh) * | 2021-05-17 | 2024-03-15 | 北京市人工影响天气中心 | 一种云粒子散射函数的探测系统、方法和装置 |
| CN114609093A (zh) * | 2022-03-24 | 2022-06-10 | 合肥工业大学 | 基于鱼眼镜头的成像浊度计的散射相函数的计算方法 |
| CN114609093B (zh) * | 2022-03-24 | 2024-04-19 | 合肥工业大学 | 基于鱼眼镜头的成像浊度计的散射相函数的计算方法 |
| CN115032176A (zh) * | 2022-05-13 | 2022-09-09 | 中国科学院南海海洋研究所 | 一种广角体散射函数及衰减系数测量方法 |
| CN115032176B (zh) * | 2022-05-13 | 2024-06-11 | 中国科学院南海海洋研究所 | 一种广角体散射函数及衰减系数测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107345893B (zh) | 2020-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107345893A (zh) | 一种粒子散射相函数测量装置及测量方法 | |
| CN104020083B (zh) | 一种确定水中悬浮颗粒物散射特性的方法 | |
| CN104266945B (zh) | 动态光散射颗粒检测装置的检测方法 | |
| CN103852078B (zh) | 空间光学姿态敏感器杂散光保护角的测量方法及装置 | |
| CN104792672A (zh) | Pm2.5质量浓度的高度分布探测装置及其探测方法 | |
| Li et al. | Monitoring optical properties of aerosols with cavity ring‐down spectroscopy | |
| CN103983571B (zh) | 探测器像素响应非均匀误差校正装置及其校正的方法 | |
| CN204556500U (zh) | 一种基于激光照射的透明液体折射率测量装置 | |
| CN107727007A (zh) | 测量两轴之间对中偏差量的方法 | |
| Crosland et al. | Instantaneous in-flame measurement of soot volume fraction, primary particle diameter, and aggregate radius of gyration via auto-compensating laser-induced incandescence and two-angle elastic light scattering | |
| Zhang et al. | Sizing of ambient particles from a single‐particle soot photometer measurement to retrieve mixing state of black carbon at a regional site of the North China plain | |
| CN105928850B (zh) | 一种光散射法颗粒粒径分布反演的索套回归方法 | |
| Sader | Turbidity measurement: A simple, effective indicator of water quality change | |
| CN103499393B (zh) | 光谱的测量方法 | |
| CN101819146B (zh) | 非接触式汽轮机蒸汽湿度和水滴直径的测量方法及其装置 | |
| CN109883930A (zh) | 一种基于偏振光散射的气溶胶颗粒折射率测量方法 | |
| CN103234942B (zh) | 利用天空散射光测量大气水平能见度的方法及装置 | |
| CN103759675A (zh) | 一种用于光学元件非球面微结构的同步检测方法 | |
| CN107894395A (zh) | 一种基于腔衰荡光谱技术的大气气溶胶消光系数测量方法 | |
| CN103389272B (zh) | 基于脉冲激光的半透明介质衰减系数和散射反照率的快速测量方法 | |
| CN105928905A (zh) | 颗粒物后向散射系数偏振敏感性测量装置 | |
| CN103344612A (zh) | 消除盛装容器对水表面brdf测量影响的测量容器和测量方法 | |
| Xie et al. | Investigating directional reflection characteristics of anisotropic machined surfaces using a self-designed scatterometer | |
| CN106546562A (zh) | 一种基于浊度检测的微生物传感器及浊度检测方法 | |
| CN111272241A (zh) | 一种基于散射光测量管道流量的系统及方法 |
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 |