close

Вход

Забыли?

вход по аккаунту

?

Патент BY14094

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2011.02.28
(12)
(51) МПК (2009)
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
G 01N 15/02
G 01N 15/06
G 01N 21/00
G 01N 21/25
G 01N 21/47
СПОСОБ ОПРЕДЕЛЕНИЯ МАССОВОЙ КОНЦЕНТРАЦИИ
РЕСПИРАБЕЛЬНОЙ ФРАКЦИИ АТМОСФЕРНОГО АЭРОЗОЛЯ
(21) Номер заявки: a 20090452
(22) 2009.03.27
(43) 2009.08.30
(71) Заявитель: Белорусский государственный университет (BY)
(72) Авторы: Кугейко Михаил Михайлович; Лысенко Сергей Александрович (BY)
BY 14094 C1 2011.02.28
BY (11) 14094
(13) C1
(19)
(73) Патентообладатель: Белорусский государственный университет (BY)
(56) BY 10845 C1, 2008.
BY 10844 C1, 2008.
RU 2148812 C1, 2000.
SU 1420474 A1, 1988.
JP 2005049228 A, 2005.
US 4525627, 1985.
(57)
Способ определения массовой концентрации респирабельной фракции атмосферного
аэрозоля с границами по верхнему размеру частиц x, равными 2,5 и 10 мкм, заключающийся в том, что измеряют сигналы обратного рассеяния из точек r зондируемой трассы
на длинах волн λ1 = 0,35 мкм, λ2 = 0,532 мкм, λ3 = 1,06 мкм, по значениям сигналов обратного рассеяния с использованием регрессионных соотношений между коэффициентами
ослабления εa(r, λi) и обратного рассеяния βπа(λi):
εa(λ1) = 1,3820·βπа(λ1) + 75,9625·βπа(λ2) + 43,7665·βπа(λ3),
εa(λ2) = 1,3738·βπа(λ1) + 52,8682·βπа(λ2) + 25,3304·βπа(λ3),
εa(λ3) = 0,6273·βπа(λ1) + 4,5588·βπа(λ2) + 28,0760·βπа(λ3),
определяют значения коэффициентов обратного рассеяния βπа(r, λi) для указанных длин
волн и определяют массовую концентрацию респирабельных фракций PM2,5(r) и PM10(r)
атмосферного аэрозоля из выражения:
PM x (r ) = ρ H x (r ) ,
где ρ - средняя плотность аэрозольного вещества,
Hx(r) - объемная концентрация аэрозольных частиц размером до x в точках r,
lnH2,5(r) = 1,4545 + 0,0353·lnβπа(r, λ1) + 0,4011·lnβπа(r, λ2) + 0,3937·lnβπа(r, λ3),
lnH10(r) = 2,1362 + 0,0379·lnβπа(r, λ1) + 0,2423·lnβπа(r, λ2) + 0,6111·lnβπа(r, λ3).
Предлагаемое изобретение относится к области охраны окружающей среды, аналитического приборостроения.
Известны способы определения концентрации аэрозольных частиц [1-4] путем регистрации проходящего через исследуемый объем или рассеянного в нем оптического излучения. В [1] о концентрации аэрозольных частиц судят по глубине модуляции рассеянного
BY 14094 C1 2011.02.28
излучения, принятого на элементах приемника, соответствующих формируемым полосам
света.
В [2, 3] о параметрах частиц аэрозоля судят по ослаблению пучка излучения на трех
длинах волн, а в [4] концентрация аэрозоля определяется по яркости принятого приемником рассеянного излучения. Известен также способ определения концентрации фонового
атмосферного аэрозоля путем измерения коэффициентов рассеяния под тремя углами [5].
Однако ни один из данных способов не позволяет определять массовые концентрации
респирабельных фракций атмосферного аэрозоля.
Наиболее близким к предлагаемому изобретению является способ определения концентрации фонового атмосферного аэрозоля путем измерения оптических характеристик коэффициентов ослабления, коэффициентов обратного рассеяния на длинах волн
λ1 = 0,35 мкм; λ2 = 0,532 мкм; λ3 = 1,06 мкм, определение его счетной и объемной концентрации транзитивной, аккумулятивной и грубодисперсной фракции из уравнений множественной регрессии [6].
Для практических потребностей санитарных служб наиболее важным является контроль массовой или объемной концентрации атмосферных аэрозолей и, в первую очередь,
их респирабельной фракции, оказывающей отрицательное воздействие на здоровье человека. Согласно стандарту Американского агентства по защите окружающей среды (EPA) и
европейскому стандарту ISO 7708 [7, 8], определены две границы отделения респирабельной фракции от общей массы аэрозоля - PM2,5 и PM10, соответствующие ограничению по
верхнему размеру частиц 2,5 и 10 мкм. Оценка загрязнения воздуха в соответствии с требованиями PM10 определяет массу частиц, проникающих в дыхательную систему человека
за гортань. PM2,5 позволяет оценить массу частиц, попадающих в альвеолярные отделы
дыхательной системы. Требуемые для санитарных служб диапазоны массовой концентрации PM2,5 и PM10 с использованием прототипа не определяются. Более того, используемые
в прототипе [6] регрессионные соотношения не соответствуют отмеченным только что
диапазонам, а также городскому аэрозолю, для которого в наибольшей степени требуется
оперативный дистанционный контроль.
Предлагаемое изобретение направлено на решение задачи оперативного дистанционного контроля респирабельных фракций как фонового, так и городского аэрозоля.
Для решения поставленной задачи в предлагаемом способе измеряют сигналы обратного рассеяния P(r, λi) в точках r зондируемой трассы на длинах волн λ1 = 0,35 мкм;
λ2 = 0,532 мкм; λ3 = 1,06 мкм, по значениям измеренных сигналов обратного рассеяния
P(r, λi) с использованием регрессионных соотношений между коэффициентами ослабления
εa(r, λi) и коэффициентами обратного рассеяния βπa(r, λi), соответствующих как фоновому,
так и городскому аэрозолю, определяют значения коэффициентов обратного рассеяния
βπa(r, λi), а массовую концентрацию респирабельных фракций PM2,5 и PM10 определяют из
выражения PM x (r ) = ρ H x (r ) ,
где x = 2,5 и 10,
ρ - средняя плотность аэрозольного вещества,
lnH2,5(r) = 1,4545 + 0,0353⋅lnβπa(r, λ1) + 0,4011⋅lnβπa(r, λ2) + 0,3937⋅lnβπa(r, λ3),
lnH10(r) = 2,1362 - 0,0379⋅lnβπa(r, λ1) + 0,2423⋅lnβπa(r, λ2) + 0,6111⋅lnβπa (r, λ3).
Hx(r) - объемная концентрация аэрозольных частиц для диапазона размеров
x < 2,5 мкм и x < 10 мкм в точках r.
Отмеченная сущность изобретения поясняется фиг. 1, 2, 3, на которых изображены
соответственно: фиг. 1 - оценки погрешностей определения εa(λi) с использованием уравнений регрессии (2); (а-в): гистограммы погрешностей восстановления εа(λi) на длинах
волн 355 (а), 532 (б) и 1064 нм (в) для 1000 реализаций параметров статистической микрофизической модели аэрозоля; (г-е): результаты восстановления εа(λi), соответствующие
данным по оптическим характеристикам аэрозоля, полученным на станциях AERONET;
2
BY 14094 C1 2011.02.28
фиг. 2 - результаты восстановления профилей βπa(r, λi) на длинах волн 355 (а), 532 (б) и
1064 нм (в) с использованием локальной калибровки в конце трассы зондирования (1 - заданные профили βπa(r, λi); 2 и 3 - соответственно профили βπa(r, λi), восстановленные в
предположении Sa(r, λi) = const и при использовании уравнений регрессии (2)); фиг. 3 - результаты восстановления концентраций H2,5 и H10 с использованием уравнений (5) при наложении на значения βπa(λi), рассчитанные для 1000 реализаций варьируемых параметров
модели, случайных погрешностей в пределах 10 %.
Для численного моделирования вариаций аэрозольных оптических характеристик и
получения регрессионных соотношений между ними использовалась статистическая микрофизическая модель аэрозоля, содержащая все необходимые для расчетов оптических
характеристик, микрофизические параметры аэрозоля, диапазоны и модель их вариаций.
Модель основана на общепринятом представлении, что частицы атмосферных аэрозолей разного происхождения различаются по своим свойствам и сравнительно слабо взаимодействуют между собой [11]. На основе такого предположения функция распределения
по размерам частиц f(a) реальных атмосферных аэрозолей представляется в виде суммы
функций распределения, описывающих отдельные фракции. Каждая фракция характеризуется концентрацией частиц в единице объема Nj, комплексным показателем преломления вещества частиц (КПП) mj = nj + iχj и функцией распределения частиц по размерам
fj(a). Многочисленные экспериментальные исследования показывают, что наиболее близко к реальной картине описание распределения аэрозольных частиц с помощью суммы
трех функций, описывающих распределение мелко-, средне- и грубодисперсной фракции
[9-12]. Для описания функции распределения частиц фракций по размерам нами использовалось логонормальное распределение [12]:
  ln a − ln a  2 
dN j
Nj
1
0j  
=
f j (a ) =
exp − 
,
(1)
 
 2  ln σ j
d ln a
2π ln σ j
 

где Nj - счетная концентрация частиц, a0j - модальный радиус, σj - стандартное отклонение, j = 1 соответствует аэрозолям естественного конденсационного происхождения, 2 фотохимическим аэрозолям, в основном антропогенного происхождения, 3 - аэрозолям
дисперсного происхождения. Первые две фракции образуются в основном в результате
фотохимических реакций и последующих конденсационных процессов. Третья фракция
образуется главным образом в результате выветривания и разрушения почв.
Таким образом, трехфракционная модель аэрозоля определяется пятнадцатью микрофизическими параметрами (по 5 для каждой фракции: вещественная и мнимая части
КПП - nj(λ) и χj(λ) соответственно, параметры a0j и σj функции распределения частиц по
размерам (1), концентрация частиц фракции Nj). Так как вторая фракция образуется в результате эволюции первой фракции, то наблюдается определенное сходство химического
состава и КПП вещества частиц первых двух фракций [9-12]. Полагая равенство КПП частиц мелко- и среднедисперсной фракций аэрозоля (тонкодисперсная фракция), число параметров модели можно сократить до 13. Диапазоны изменений указанных параметров,
выбранные на основании анализа многочисленных литературных данных [9-16], охватывают значения параметров аэрозолей как в регионах с малой антропогенной нагрузкой
(фоновый аэрозоль), так и в крупных индустриальных районах. При фиксированных значениях параметров модели можно по известным формулам теории Ми вычислять оптические характеристики каждой фракции и находить общие характеристики.
При каждом конкретном акте использования модели по известным правилам моделирования равномерно распределенных случайных величин для каждого набора из 13 параметров рассчитываются требуемые оптические аэрозольные характеристики. Однако при
получении выборки микрофизических параметров путем их случайных вариаций могут
получаться не встречающиеся в реальности расчетные значения как оптических, так и не-
3
BY 14094 C1 2011.02.28
которых интегральных микрофизических характеристик, например объемных концентраций, эффективных размеров r32, лидарного отношения, коэффициента ослабления, его
спектрального хода и т.д. Чтобы "отфильтровать" выборку, т.е. исключить не встречающиеся в реальности значения отмеченных выше параметров, нами накладывались ограничения на значения:
1) коэффициента ослабления: 0,005 ≤ ε(0,55 мкм) ≤ 1 км-1;
2) отношения коэффициентов обратного рассеяния: 0,7 ≤ βπ(0,55)/βπ(1,02) ≤ 2,5;
ln ε(λ1 ) − ln ε(λ 2 )
3) параметра Ангстрема ν =
(λ1 = 355 нм, λ2 = 532 нм): 0,6 ≤ ν ≤ 2,2;
ln λ 2 − ln λ1
4) отношения объемных концентраций частиц грубодисперсной и тонкодисперсной
фракций: 0,1 ≤ C V3
2
∑ CVi ≤ 10 ;
i =1
5) эффективного размера частиц аэрозоля rt =
a max
∫
a min
a
3
3
rmax
i =1
rmin
∑ fi (a )d ln a
∫
3
a 2 ∑ f i (a )d ln a , а
i =1
также его тонкодисперсной rf и грубодисперсной rc фракций: 0,1 мкм ≤ rt ≤ 1 мкм;
0,1 мкм ≤ rf ≤ 0,33 мкм; 1 мкм ≤ rc ≤ 5 мкм;
Данные ограничения выбраны на основании анализа многочисленных данных сети
глобального мониторинга на основе сканирующих солнечных фотометров AERONET [17,
18].
На основе сгенерированного ансамбля из 1000 реализаций микрофизических параметров аэрозоля и расчетов его коэффициентов ослабления εa(λi) (км-1) и обратного рассеяния
βπa(λi) (км-1ср-1) на длинах волн λ1 = 355 нм; λ2 = 532 нм и λ3 = 1064 нм, с использованием
методов математической статистики, получены уравнения множественной регрессии между ними:
ε0(λ1) = 1,3820⋅βπa(λ1) + 75,9625⋅βπa(λ2) + 43,7665⋅βπa(λ3),
(2)
ε0(λ2) = –1,3738⋅βπa(λ1) + 52,8682⋅βπa(λ2) + 25,3304⋅βπa(λ3),
εа(λ3) = 0,6273⋅βπa(λ1) + 4,5588⋅βπa(λ2) + 28,0760⋅βπa(λ3).
Для каждой реализации по уравнениям регрессии (2), с учетом добавления в расчетные значения βπa(λi) случайных погрешностей измерения (в пределах 10 %), рассчитывались значения εa(λi), которые сравнивались с их точными значениями. Соответствующие
гистограммы погрешностей определения εa(λi) из уравнений (2) представлены на фиг. 1, а-в.
На рисунках также указаны значения средних по выборке микрофизических параметров
аэрозоля погрешностей восстановления εa(λi).
Следует отметить, что, в силу широких диапазонов вариаций аэрозольных микрофизических характеристик, регрессионные уравнения (2) справедливы как для городских, так
и фоновых аэрозолей. Очевидно, что при меньших диапазонах (соответствующих региональным особенностям аэрозоля) коэффициенты ослабления можно определять с большей
точностью.
Интерпретация данных многоволнового аэрозольного лидара в задачах исследования
пространственно-временной структуры атмосферного аэрозоля основывается на обращении системы уравнений лазерного зондирования, связывающих мощность принимаемых
сигналов обратного рассеяния P(r, λi) с параметрами среды. В приближении однократного
рассеяния данные уравнения имеют вид [19]:

 r
−2
P(r, λ i ) = G (λ i ) ⋅ r ⋅ [βπa (r, λi ) + βπm (r, λ i )] ⋅ exp− 2 ∫ [εa (r, λ i ) + ε m (r, λ i )]dx ,
(3)

 r0
где r - текущая координата по трассе зондирования; λi - длина волны; G(λi) - константа калибровки; βπa(r, λi) и βπm(r, λi) - соответственно профили коэффициентов обратного аэро-
4
BY 14094 C1 2011.02.28
зольного и молекулярного рассеяния; εa(r, λi) и εm(r, λi) - профили коэффициентов аэрозольного и молекулярного ослабления; расстояние r0 соответствует концу "теневой зоны"
лидара.
Для определения константы калибровки лидара необходима априорная информация
об оптических характеристиках среды. Исходя из общего свойства устойчивости решения
(3), опорное значение восстанавливаемой оптической характеристики более предпочтительно задавать в конце трассы [19]. При известном опорном значении βπa(rk, λi) восстановление профилей βπa(r, λi) удобно осуществлять с использованием итерационного
метода Фернальда [20], требующего задания профиля лидарного отношения Sa(r, λi) =
εa(r, λi)/βπa(r, λi) и использующего предположение о его постоянстве по трассе зондирования.
Профили β(π1a) (r, λ i ) , полученные при первоначальном допущении S(a1) (r, λ i ) = const ,
рассматриваются как первое приближение. Далее в соответствии с уравнениями (2) вычисляются профили лидарных отношений Sa( 2) (r, λi ) следующей итерации и проводится
уточнение профилей β(π2a) (r, λi ) с использованием получаемых на основе метода Фернальда соотношений:
β(πna) (r, λ i ) = −β πm (r, λ i ) +
P ( r, λ i ) Y ( n ) ( r , λ i )
r
,
P(rk , λ i )Y ( n ) (rk , λ i ) Sa( n ) (r, λ i )
(n )
− 2Sa (r, λ i ) ∫ P( x , λ i )Y ( n ) ( x , λ i )dx
(n)
(n)
β πa (rk , λ i ) + β πm (rk , λ i ) Sa (rk , λ i )
rk
[
]
(4)
 r

где Y ( n ) (r, λ i ) = S(an ) (r, λ i ) ⋅ r 2 ⋅ exp− 2 ∫ (S(an ) ( x , λ i ) − Sm )β πm ( x , λ i )dx .
 r0

Критерием завершения итерационной процедуры является выполнение условия:
rmax
∫
β(πna) ( x , λ i ) − β(πna−1) ( x , λ i ) dx
r0
rmax
∫
β(πna) ( x , λ i )dx ≤ δ,
r0
где δ - заданная малая величина.
Определение характеристик H2,5(r) и Н10(r) в точках r исследуемой трассы основано на
множественной корреляции между ними и значениями коэффициента обратного рассеяния на указанных длинах волн, установленной нами на основе рассмотренной выше микрофизической модели аэрозоля с использованием аппарата математической статистики
для выборки из 1000 реализаций варьируемых параметров модели:
lnH2,5(r) = 1,4545 + 0,0353⋅lnβπa(r, λ1) + 0,4011⋅lnβπa(r, λ2) + 0,3937⋅lnβπa(r, λ3),
(5)
lnH10(r) = 2,1362 - 0,0379⋅lnβπa(r, λ1) + 0,2423⋅lnβπa(r, λ2) + 0,6111⋅lnβπa(r, λ3).
В (5) используются логарифмы величин, т.к. в этом случае мы имеем дело с числами
одного порядка, а абсолютная ошибка логарифма дает оценку относительной ошибки самой
величины. Использование уравнений (5) позволяет осуществлять оперативный контроль
концентрации респирабельной фракции аэрозоля по данным трехчастотного зондирования
без решения некорректных обратных задач, использования дополнительной априорной информации и сопутствующих измерений. Широкие диапазоны вариаций микрофизических
параметров модели позволяют говорить о справедливости уравнений (5) для аэрозолей как в
сельской местности, так и в крупных индустриальных районах.
На фиг. 2 представлены результаты восстановления характеристик H2,5 и H10 для
1000 реализаций параметров модели аэрозоля, при наложении на расчетные значения
βπa(λi) случайных погрешностей в пределах 10 %. Средние погрешности их восстановления, равные для Н2,5 и H10 соответственно 35 и 33 %, ввиду широких диапазонов возможных значений искомых характеристик (в два порядка), являются вполне приемлемыми для
удовлетворения потребностей ряда задач экологического мониторинга окружающей среды (например, по превышению ПДК).
5
BY 14094 C1 2011.02.28
Из выражений PM x (r ) = ρ H x (r ) , где x = 2,5 и 10, ρ - средняя плотность аэрозольного
вещества,
lnH2,5(r) = 1,4545 + 0,0353⋅lnβπa(r, λ1) + 0,4011⋅lnβπa(r, λ2) + 0,3937⋅lnβπa(r, λ3),
lnH10(r) = 2,1362 - 0,0379⋅lnβπa(r, λ1) + 0,2423⋅lnβπa(r, λ2) + 0,6111⋅lnβπa(r, λ3),
Hx - объемная концентрация аэрозольных частиц для диапазона размеров x < 2,5 мкм и
x < 10 мкм,
легко определить массовые концентрации респирабельных фракций атмосферного аэрозоля.
Измерения сигналов обратного рассеяния на трех отмеченных длинах волн и использование уравнений регрессии, связывающих значения коэффициентов аэрозольного ослабления и обратного рассеяния на длинах волн 355, 532 и 1064 нм позволяет значительно
повысить точность восстановления профилей коэффициентов обратного рассеяния. Использование же множественной корреляции между объемными концентрациями H2,5 и H10
и спектральными значениями коэффициента аэрозольного обратного рассеяния позволяет
проводить достаточно точные оценки экологически значимых характеристик аэрозоля по
данным трехчастотного зондирования атмосферы без решения некорректных обратных
задач и использования дополнительной априорной информации, что также повышает точность определения H2,5 и H10, поскольку данные задачи чувствительны к априорной информации и погрешностям определения используемых оптических характеристик.
Предлагаемый способ может эффективно использоваться для организации оперативного
экологического мониторинга как городского, так и фонового аэрозолей.
Источники информации:
1. RU 2148812 C1, 2000.
2. SU 1420474 A1, 1988.
3. JP 2005049228 А, 2005.
4. US 4525627, 1985.
5. BY 10844 C1, 2008.
6. U.S.EPA Guidance for using continuous monitors in PM2.5 monitoring networks.
OAQPS EPA-454/R-98-012, May 1998.
7. European Standard EN 12341. Air quality. Determination of the PM10 fraction of suspended particulate matter. Reference method and fied test procedure to demonstrate reference
equivalence of measurment methods, 1998.
8. Патент РБ 10845 C1, 2008.
9. Кондратьев К.Я., Ивлев Л.С., Крапивин В.Ф. Свойства, процессы образования и последствия воздействий атмосферного аэрозоля: от нано- до глобальных масштабов. СПб.: BBM, 2007. - С. 807.
10. Ивлев Л.С., Андреев С.Д. Оптические свойства атмосферных аэрозолей. - Ленинград: ЛГУ, 1986. - С. 359.
11. Андреев С.Д., Ивлев Л.С. Моделирование оптических характеристик аэрозолей
приземного слоя атмосферы в области 0,3-15 мкм. Ч. I. Принципы построения модели.
Ч. II. Выбор параметров модели. Ч. III. Результаты моделирования // Опт. атмосф. и океана. - 1995. - T. 8. - № 5. - С. 788-795. - № 8. - С. 1227-1235; 1236-1243.
12. Зуев В.Е., Креков Г.М. Оптические модели атмосферы. - Л.: Гидрометеоиздат,
1986. - С. 256.
13. Васильев А.В., Ивлев Л.С. Оптическая статистическая модель атмосферы для района Ладожского озера // Опт. атмосф. и океана. - 2000. - T. 13. - № 2. - С. 198-203.
14. Ивлев Л.С., Васильев А.В., Белан Б.Д., Панченко М.В., Терпугова С.А. Оптикомикрофизические модели городских аэрозолей. В сб. "Третья международная конферен-
6
BY 14094 C1 2011.02.28
ция. Естественные и антропогенные аэрозоли, сборник трудов. Санкт-Петербург, 24.0927.09.2001"/ Отв.ред. Л.С. Ивлев, НИИ Химии СПбГУ, 2003. - С. 161-170.
15. Chunlei L., Smith M.H. Urban and rural aerosol particle optical properties // Atmos. Environ. - 1995. - V. 29. No. 1-2. - P. 3293-3301.
16. Иванов В.П. Прикладная оптика атмосферы в тепловидении. - Казань: Новое Знание, 2000. - С. 357.
17. Cattrall C., Reagan J., Thome K., Dubovic O. Variability of aerosol spectral lidar and
backscatter and extinction ratio of key aerosol types derived from selected Aerosol Robotic Network locations // J. Geophys. Res. - 2005. - Vol. 110, D10S11. doi: 10.1029/2004 JD005124.
18. http://aeronet.gsfc.nasa.gov.
19. Креков Г.М., Кавкянов С.И., Крекова М.М. Интерпретация сигналов оптического
зондирования атмосферы. - Новосибирск: Наука, 1987. - С. 184.
20. Fernald F.G. Analysis of atmospheric lidar observation: some comments // Appl. Opt.
1984. - V. 23. - No. 5. - P. 652-653.
Фиг. 1
Фиг. 2
7
BY 14094 C1 2011.02.28
Фиг. 3
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
8
Документ
Категория
Без категории
Просмотров
0
Размер файла
226 Кб
Теги
by14094, патент
1/--страниц
Пожаловаться на содержимое документа