close

Вход

Забыли?

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

?

Патент BY10844

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
BY (11) 10844
(13) C1
(19)
(46) 2008.06.30
(12)
(51) МПК (2006)
G 01N 15/02
G 01N 15/06
G 01N 21/00
G 01N 21/25
G 01N 21/47
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
СПОСОБ ОПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИИ ФОНОВОГО
АТМОСФЕРНОГО АЭРОЗОЛЯ
BY 10844 C1 2008.06.30
(21) Номер заявки: a 20070061
(22) 2007.01.24
(43) 2007.10.30
(71) Заявитель: Белорусский государственный университет (BY)
(72) Авторы: Кугейко Михаил Михайлович; Лысенко Сергей Александрович (BY)
(73) Патентообладатель: Белорусский государственный университет (BY)
(56) Sviridenkov M.A. etc. Comparison of
Black Carbon Content, Aerosol Optical
and Microphysical Characteristics in
Moscow and the Moscow Region. Fifteenth ARM Science Team Proceedings,
Daytona Beach, Florida, March 14-18,
2005.
RU 2148812 C1, 2000.
RU 2155954 C2, 2000.
SU 1437746 A1, 1988.
US 4525627, 1985.
(57)
Способ определения концентрации фонового атмосферного аэрозоля путем измерения
коэффициента рассеяния оптического излучения, установления функциональных соотношений между концентрацией и измеряемой оптической характеристикой, отличающийся
тем, что измеряют коэффициенты рассеяния D оптического излучения под углами, равными 5°, 20° и 165°, на длине волны, равной 0,532 мкм, а значения счетной N и объемной
Сυ концентраций транзитивной, аккумулятивной и грубодисперсной фракций фонового
атмосферного аэрозоля определяют из уравнений множественной регрессии:
N = A D,
где
[
Cυ = B D,
( )
( )
D = 1 log D 5o log D 20o
5,277 0,619 − 0,938
A = 2,555 0,824 − 0,937
1,373 1,305 − 0,861
 − 0,733 − 0,294 1,319
B =  − 1,894 1,378 − 0,445
− 0,285 1,255 − 0,480
( )] ,
log D 65o
0,934
0,873,
0,774
− 0,186
− 0,156,
0,473 
N = [log N 1 log N 2 log N 3 ] ,
T
Cυ = [log Cυ 1 log Cυ 2 log Cυ 3 ] ,
T
T
BY 10844 C1 2008.06.30
N1, N2, N3 – соответственно счетные концентрации транзитивной, аккумулятивной и
грубодисперсной фракций фонового атмосферного аэрозоля,
Сυ1, Сυ2, Сυ3 – соответственно объемные концентрации транзитивной, аккумулятивной
и грубодисперсной фракций фонового атмосферного аэрозоля.
Предлагаемое изобретение относится к области охраны окружающей среды, аналитического приборостроения.
Известны способы определения концентрации аэрозольных частиц [1-4] путем регистрации проходящего через исследуемый объем или рассеянного в нем оптического излучения. В [1] о концентрации аэрозольных частиц судят по глубине модуляции рассеянного
излучения, принятого на элементах приемника, соответствующих формируемым полосам
света. В [2] массовую концентрацию тонкодисперсных аэрозолей определяют путем зондирования их ультрафиолетовым излучением с длиной волны λ = 266 нм и регистрации интенсивностей сигналов люминесценции в области спектра 300-400 нм. О концентрации двух
фракций угольно-водного аэрозоля судят по измеряемому коэффициенту ослабления в видимой области спектра и степени линейной поляризации излучения, рассеянного аэрозолем
в направлении, составляющем угол 60 ± 20° направлением зондирующего излучения [3]. В
[4] масса аэрозольных частиц в газообразных средах определяется по разности интенсивностей излучения на длинах волн в диапазонах 3,8 ÷ 4,15 и 3,35 ÷ 3,5 мкм, проходящих через
исследуемую среду и эталонный образец. Однако ни один из данных способов не позволяет
определять счетные и объемные концентрации всех фракций атмосферного аэрозоля.
Наиболее близким к предлагаемому является способ определения концентрации аэрозоля путем измерения коэффициента рассеяния оптического излучения под углом 45° на
длине волны λ = 0,51 мкм и использования установленной экспериментально функциональной связи между ним и определяемой характеристикой [5]. Данный способ позволяет
определять объемную концентрацию атмосферного аэрозоля с радиусами r = 0,05-0,6 мкм.
К недостаткам метода следует отнести невозможность использования его при наличии
частиц грубодисперсной фракции (с радиусами r > 1 мкм); невозможность определения
фракционных концентраций атмосферного аэрозоля (транзитивной - r ∼ 0,01 мкм, аккумулятивной - r ~ 0,1 мкм, грубодисперсной - r > 1 мкм) по отдельности; может использовать
только в атмосферных ситуациях, соответствующих той, для которой установлена используемая функциональная связь.
Предлагаемое изобретение направлено на решение задачи одновременного определения как объемных, так и счетных фракционных концентраций фонового атмосферного аэрозоля в различных атмосферных ситуациях.
Для решения поставленной задачи в предлагаемом способе измеряют коэффициенты
рассеяния оптического излучения D(θ) под углами θ, равными 5°, 20° и 165° на длине
волны λ = 0,532 мкм, а значения счетной N и объемной Сυ концентраций транзитивной,
аккумулятивной и грубодисперсной фракций фонового атмосферного аэрозоля определяют из уравнений множественной регрессии:
r
r r
r
(1)
N = AD, C υ = BD,
где
r
r
N = [log N1 log N 2 log N 3 ]T ; C υ = [log C υ1 log C υ2 log C υ3 ]T ;
( )
( )
( )]
D = [1 log D 5o log D 20o log D 165o ;
5,277 0,619 − 0,938 0,934
 − 0,733 − 0,294 1,319 − 0,186


A = 2,555 0,824 − 0,937 0,873 , B =  − 1,894 1,378 − 0,445 − 0,156 ,
1,373 1,305 − 0,861 0,774
− 0,285 1,255 − 0,480 0,473 
2
T
BY 10844 C1 2008.06.30
N1, N2, N3 соответственно счетные концентрации транзитивной, аккумулятивной и грубодисперсной фракции фонового атмосферного аэрозоля,
Cυ1, Cυ2, Сυ3 - соответственно объемные концентрации транзитивной, аккумулятивной и
грубодисперсной фракции фонового атмосферного аэрозоля.
Отмеченная сущность изобретения поясняется фиг. 1, 2, 3, на которых изображены
соответственно: фиг. 1 - совокупность различных значений коэффициентов рассеяния под
углом 45° на длине волны λ = 532 нм и соответствующих им концентраций трех фракций;
фиг. 2 - угловые зависимости коэффициентов корреляции между logD(θ) и logNi, logCυi
(i = 1, 2, 3); фиг. 3 - гистограммы ошибок определения Сυi и Ni с использованием полученных уравнений множественной регрессии (1).
Действительно, при известных микрофизических параметрах (МП) аэрозоля (функция
распределения частиц на размерах f(x), показатель преломления m, концентрация) соотношения, связывающие Cυ и N аэрозоля с коэффициентом рассеяния под углом D(θ),
можно записать в следующем виде:
∞
C υ = D(θ) ⋅
2λ
3π ∞
3
∫ x f (x )dx
0
2
∫ x Q(m, x )f (x )dx
,
N=
0
8π 2
λ3
∞
D(θ)
,
2
∫ x Q(m, x )f (x )dx
0
где х = 2πr/λ - параметр Ми; λ - длина волны; Q(m,x) - фактор эффективности рассеяния
излучения частицей с радиусом r; θ - угол рассеяния.
Отметим, что в используемой нами в дальнейшем аэрозольной модели коэффициенты
µυ(λ) = D(θ)/Сυ и µN(λ) = D(θ)/N при каждой комбинации МП рассчитывались в приближении сферических частиц по известным формулам Ми. Такое упрощение обусловлено
несколькими причинами. Во-первых, объем экспериментальных данных по форме и внутренней структуре аэрозольных частиц явно недостаточен для построения соответствующих моделей. Во-вторых, методики и алгоритмы расчета оптических характеристик (ОХ)
несферических и неоднородных рассеивателей громоздки и имеют лишь ограниченную
область применимости. Однако для сфероидов и цилиндров с умеренным параметром
асферичности k (от 0,5 до 2,0) индикатриса рассеяния света, усредненная по различным
пространственным ориентациям, совпадает с индикатрисой шара того же объема. Отметим также, что для несферических частиц эквивалентные по площади и по объему радиусы не совпадают, однако для умеренно вытянутых (сплюснутых) сфероидов и цилиндров
(0,5 ≤ k ≤ 2) эта разница не превышает 5 %, что позволяет рассчитывать объемную концентрацию, не учитывая форму частиц.
Что касается использования аналитической зависимости f(r), то выбор последней малосуществен. Важны модальные радиусы частиц и крутизна спада концентрации частиц
при размерах частиц больше модального.
Для нахождения указанных корреляций с учетом взаимозависимых вариаций микрофизических параметров аэрозолей нами использовалась трехфракционная оптическая модель приземного слоя атмосферы Андреева-Ивлева [6]. Модальные размеры частиц
фракций, отличающихся по своему происхождению, составу и всему комплексу физикохимических свойств, оцениваются как rmi ≈ 0,03, 0,3 и 1 мкм при влажности воздуха
q = 50 % для транзитивной, аккумулятивной и грубодисперсной фракций соответственно
(в дальнейшем для обозначения данных фракций нами используются индексы 1, 2 и 3 соответственно). Для описания распределения по размерам частиц отдельных фракций аэрозолей используется гамма-распределение с отрицательными степенями свободы. Выбор
этой модели обусловлен тем, что ее входными параметрами являются метеорологическая
дальность видимости Sm (позволяет осуществить привязку используемой в расчете ОХ
концентрации аэрозолей к экспериментально измеренному аэрозольному ослаблению в
3
BY 10844 C1 2008.06.30
видимой области спектра) и влажность q (к которой привязываются вариации спектра
размеров и изменение показателя преломления mi(λ)), позволяющие в комплексе менять
МП атмосферы. При расчетах предполагается, что вода равномерно перемешана в частице
(микрокапиллярное внедрение и разбухание частицы), а зависимость размера частиц от q
описывается эмпирической формулой Кастена [7]. Таким образом, варьируя эти два параметра (что соответствует изменениям условий в атмосфере), можно отследить наличие
или отсутствие устойчивых корреляций между D(θ) и Ni, Cυi(i = 1, 2, 3).
Значения D и Cυ, полученные нами на основе модели [6], при различных комбинациях
Sm и q и рассчитанная нами по этим данным корреляционная кривая (Сυ = 0,87D0,98) хорошо совпадают с зависимостью, полученной в [5] (Cυ = D). Это свидетельствует о том, что
данная модель достаточно адекватно описывает взаимосвязь оптических и микрофизических характеристик аэрозоля и может быть использована для исследования в рамках поставленной задачи.
На фиг. 1 в логарифмическом масштабе изображена совокупность различных значений коэффициентов рассеяния под углом 45° на длине волны λ = 532 нм и соответствующих им концентраций трех фракций, отображающие разброс коэффициентов связи между
ними. Данная зависимость получалась следующим образом: задавалось случайным образом состояние среды, характеризуемое параметрами Sm и q из диапазонов их значений
5 ÷ 60 км, 10 ÷ 99 %, по теории Ми рассчитывался коэффициент D(45°, 0,532 мкм) и по
стандартным формулам математической статистики [8] проводилась линейная аппроксимация полученной совокупности точек. Расчет Ni и Cυi осуществлялся через весовые
функции фракций аэрозолей при различных значениях Sm и q, полученных в [6] путем сопоставления модельных и экспериментально измеренных спектров коэффициента ослабления аэрозолем. Помимо задания Sm и q по равномерному закону распределения на все
МП модели также накладывался 10 %-й случайный разброс (не зависящий от Sm и q), что
позволяет получить более сглаженные и устойчивые к ошибкам коэффициенты множественной регрессии и еще больше расширить диапазоны вариаций всех МП модели. Как
видно из фиг. 1 б, в, значение D(45°, 0,532 мкм) достаточно тесно коррелирует с Сυ2 и Сυ3,
что позволяет с приемлемой точностью определять данные концентрации по результатам
оптических измерений. Однако большой разброс точек на фиг. 1 а, г, д, е свидетельствует
о том, что D(45°, 0,532 мкм) практически не информативен относительно счетных концентраций всех трех аэрозольных фракций и объемной концентрации транзитивной фракции.
С целью нахождения наиболее перспективных, в плане повышения информативности
нефелометрических измерений, углов рассеяния оптического излучения θ нами были исследованы угловые зависимости коэффициентов парной корреляции ρNi(θ) и ρυi(θ) между
logD(θ) и logNi, logCυi(i = 1, 2, 3) соответственно, определяемых согласно [8] по формуле:
∑ (log Fij − log Fi )(log D j (θ) − log D(θ))
n
ρ Fi (θ) =
j =1
∑ (log Fi
n
j =1
j
− log Fi
) ∑ (log D (θ) − log D(θ))
2 n
j =1
;
2
j
где j = 1, 2,…, n = 1000 - реализация состояния модельной атмосферы, с соответствующими
ей коэффициентом направленного рассеяния Dj(θ) и фракционными концентрациями
Fij N ij или C υj i ; logD(θ ) ; log Fi - математические ожидания логарифмов данных величин.
Угловые зависимости данных коэффициентов представлены на фиг. 2. Как следует из
фиг. 2 а для определения счетных концентраций аэрозольных фракций предпочтительней
осуществлять измерения коэффициента рассеяния под углами из диапазона 160° ÷ 170°,
характеризуемыми высокими значениями соответствующих коэффициентов корреляции
для всех трех фракций. Обращают также на себя внимание высокие значения ρC1(θ) и
(
)
4
BY 10844 C1 2008.06.30
ρC2(θ) при малых углах θ ≈ 1° - для аккумулятивной и θ ≈ 5° - для грубодисперсной фракций аэрозоля, а также наличие максимума коэффициента корреляции для транзитивной
фракции ρC3(θ) в районе 20°.
Таким образом, принимая в качестве критерия оптимальности выбора углов рассеяния
максимальную информативность коэффициентов рассеяния под данными углами относительно искомых параметров (наибольшие значения соответствующих коэффициентов
корреляции) в сочетании с минимальным количеством измеряемых ОХ следует, что эффективнее проводить измерения D(θ) одновременно под малыми углами и углами θ = 20°,
160° ÷ 170°. Полученные нами с применением стандартных формул математической статистики [8] корреляционные соотношения между значениями D(θ) для углов θ из отмеченных диапазонов и проведенная оценка точности восстановления счетной и объемной
концентрации с их использованием показывают, что наименьшая погрешность соответствует углам θk = 5°, 20°, 165°. Для данных углов корреляционные соотношения между
D(θk) и счетными Ni, объемными Сυi концентрациями трех фракций имеют вид:
r
r r
r
(1)
N = A ⋅ D, C υ = B ⋅ D,
где
r
r
N = (log N1 log N 2 log N 3 )T ; C υ = (log C υ1 log C υ2 log C υ3 )T ;
r
T
D = (1 log D 5o log D 20o log D 165o ;
5,277 0,619 − 0,938 0,934
 − 0,733 − 0,294 1,319 − 0,186


A = 2,559 0,824 − 0,937 0,873 , B =  − 1,894 1,378 − 0,445 − 0,156 .
1,373 1,305 − 0,861 0,774
− 0,285 1,255 − 0,480 0,473 
Для оценки эффективности предлагаемой методики определения Сυi и Ni проведен
расчет данных характеристик в 1000 модельных атмосферных ситуациях, отличающихся
значениями Sm и q и при наложении на значения D(θk) случайной ошибки в пределах
10 %. Гистограммы ошибок определения Сυi и Ni с использованием уравнений (2) приведены на фиг. 3, где PCi и PNi - частоты повторяемости ошибок δСυi и δNi соответственно
(число состояний, для которых получены данные ошибки). На рисунках также указаны
среднеквадратичные ошибки определения отмеченных только что величин. Как видно,
среднеквадратичные ошибки измерений фракционных концентраций аэрозоля по описанной методике не превышают суммарные ошибки измерения коэффициентов рассеяния под
указанными углами, что удовлетворяет многим практическим потребностям.
Для реализации предложенного способа необходимо измерить ОХ коэффициенты рассеяния под углами θ, равными 5°, 20° и 165° на длине волны λ = 0,532 мкм, а значения
счетной N и объемной Сυ концентраций транзитивной, аккумулятивной и грубодисперсной фракций фонового атмосферного аэрозоля определить из уравнений множественной
регрессии (1). Для измерения коэффициентов рассеяния под указанными углами θ можно
использовать известные измерители коэффициентов рассеяния, индикатрисы рассеяния
[9-11]. Как следует из [10, 11], погрешность измерения коэффициентов рассеяния под углом θ, не превышающая 15 % в области малых углов, - 5° (из-за использования светофильтров при измерении световых потоков, изменяющихся на несколько порядков) и 5 %
в остальной части индикатрисы, что в сумме укладывается в погрешности, накладываемые
нами на измеряемые коэффициенты рассеяния под тремя углами θ при получении регрессионных соотношений (суммарная накладываемая погрешность составляла 30 %, можно
же измерять коэффициенты рассеяния с суммарной погрешностью 25 %).
Предлагаемый способ определения фракционных концентраций атмосферного аэрозоля может быть с успехом использован для оперативного и автоматизированного контроля
состояния атмосферы в сельской местности и в районах, где мал вклад антропогенных аэрозолей, вместо традиционных методик отбора проб с последующим анализом. Слабая
( )
( )
5
( ))
BY 10844 C1 2008.06.30
чувствительность способа к погрешностям измерений D(θk) делает его использование в
большинстве случаев гораздо более предпочтительным по сравнению с современными оптическими методами, требующими или проведения калибровочных измерений, или решения некорректных обратных задач, а значит, обладающих большим коэффициентом
усиления ошибок.
Таким образом, как следует из вышесказанного, предлагаемый способ позволяет определять концентрации одновременно всех трех фракций атмосферного аэрозоля в различных атмосферных ситуациях, к тому же не требует проведения дополнительных
(калибровочных) измерений по установлению коэффициентов связи между отдельными
характеристиками среды и измеряемыми коэффициентами рассеяния (используемые здесь
регрессионные уравнения получены для максимально возможного разброса микрофизических характеристик фонового атмосферного аэрозоля).
Источники информации:
1. RU 2148812 C1, 2000.
2. RU 2155954 С2, 2000.
3. SU 1437746 А1, 1988.
4. US 4525627, 1985.
5. Sviridenkov M.A. etc Comparison of Black Carbon Content, Aerosol Optical and Microphysical Characteristics in Moscow and Moscow Region. Fifteenth ARM Science Team Meeting
Proceedings, Daytona Bech, Florida, March 14-18, 2005.
6. Андреев С.Д., Ивлев Л.С. Моделирование оптических характеристик аэрозолей приземного слоя атмосферы // Оптика атмосферы и океана. - 1995. - Т. 8. - № 8. - С. 1236-1243.
7. Креков Г.М., Звенигородский С.Г. Оптическая модель средней атмосферы. - Новосибирск: Наука, 1990. - 278 с.
8. Колемаев В.А., Калинина В.Н. Теория вероятностей и математическая статистика. М.: ИНФРА-М, 1997. - 302 с.
9. Иванов А.П. Оптика рассеивающих сред. - Мн.: Наука и техника, 1969. - 592 с.
10. Hammer М., Yaroslavsky A.N., Schweitzer D. // Physics in Medicine and Biology. 2001. - V. 46. - N 65. - P. 134-142.
11. Александров А.Ф., Свистов Д.А., Тимофеев И.Б. // Биомедицинская радиоэлектроника. - 2000. - № 5. - С. 50-64.
Фиг. 1
6
BY 10844 C1 2008.06.30
Фиг. 2
Фиг. 3
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
7
Документ
Категория
Без категории
Просмотров
1
Размер файла
434 Кб
Теги
by10844, патент
1/--страниц
Пожаловаться на содержимое документа