close

Вход

Забыли?

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

?

Патент BY11892

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2009.04.30
(12)
(51) МПК (2006)
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
G 01J 3/44
СПОСОБ ОПРЕДЕЛЕНИЯ ПОКАЗАТЕЛЯ ОСЛАБЛЕНИЯ
РАССЕИВАЮЩЕЙ СРЕДЫ С СИЛЬНО ВЫТЯНУТОЙ
ИНДИКАТРИСОЙ РАССЕЯНИЯ
(21) Номер заявки: a 20070328
(22) 2007.03.29
(43) 2008.10.30
(71) Заявитель: Государственное научное
учреждение "Институт физики имени Б.И.Степанова Национальной
академии наук Беларуси " (BY)
(72) Авторы: Науменко Елена Константиновна (BY)
BY 11892 C1 2009.04.30
BY (11) 11892
(13) C1
(19)
(73) Патентообладатель: Государственное
научное учреждение "Институт физики имени Б.И.Степанова Национальной академии наук Беларуси " (BY)
(56) ИВАНОВ А.П. Оптика рассеивающих
сред. - Мн.: Наука и техника, 1969. С. 80, 462, 469.
BY 3670 С1, 2000.
SU 1300302 А1, 1987.
SU 1430839 А1, 1988.
JP 2004163299 A, 2004.
US 4986656 A, 1991.
НАУМЕНКО Е.К. Журнал прикладной
спектроскопии. - 1996. - Т. 63. - № 2. С. 202-209.
(57)
Способ определения показателя ослабления рассеивающей среды с сильно вытянутой
индикатрисой рассеяния, заключающийся в том, что берут два образца исследуемой рассеивающей среды; для первого образца задают величину a1 = (C V )1 L1 ,
где (C V )1 - объемная концентрация частиц в первом образце, равная доле объема, занимаемого частицами в единице объема рассеивающей среды,
L1 - геометрическая толщина слоя первого образца;
для второго образца задают величину a 2 = (C V ) 2 L 2 ,
где (C V ) 2 – объемная концентрация частиц во втором образце,
L 2 - геометрическая толщина слоя второго образца;
Фиг. 1
BY 11892 C1 2009.04.30
при освещении образцов коллимированным пучком света измеряют коэффициент пропускания для первого образца T1 и коэффициент пропускания для второго образца T2 ; определяют оптическую толщину τ1 слоя толщиной L1 из уравнения:
 a
 a
T2 = exp − 2 τ1  + 2 [T1 − exp(− τ1 )] ;
 a1  a1
затем определяют показатель ослабления ε′ рассеивающей среды, нормированный на
единичную объемную концентрацию частиц в среде, из выражения:
τ
τ1
ε′ = 1 =
;
a1 (C V )1 L1
причем объемную концентрацию частиц в образцах и геометрическую толщину слоя образцов выбирают с возможностью выполнения условий однократного рассеяния.
Изобретение относится к экспериментальной физике, в частности к определению показателя ослабления рассеивающих сред с сильно вытянутой индикатрисой рассеяния.
Показатель ослабления является одной из важнейших оптических характеристик, который широко используется для оптической диагностики структуры и состава природных
и технических рассеивающих сред (атмосферный аэрозоль, гидрозоли внутренних и открытых водоемов, суспензии биологических клеток, порошкообразных материалов, пигментов, используемых в печатной и лакокрасочной промышленности и мн.др.). Точность
решения обратных задач оптики рассеивающих сред в значительной степени определяется
(зависит) точностью измерения оптических характеристик, по величине которых восстанавливается дисперсный состав рассеивающей среды и (или) оптические свойства
вещества рассеивающих частиц. Наличие даже незначительных ошибок приводит к неустойчивости систем уравнений, разбалтыванию решений и неоднозначным результатам
восстановления искомых параметров. Поэтому вопрос о точности измерений характеристик рассеянного излучения, в том числе и показателя ослабления, имеет важное значение
для оптической диагностики дисперсных сред.
Известен нефелометрический способ измерения показателя ослабления для непоглощающих и слабопоглощающих рассеивающих сред (при этих условиях показатель ослабления практически равен показателю рассеяния), основанный на измерении интенсивности рассеянного излучения в некотором заданном направлении и определении
показателя ослабления с использованием корреляционного соотношения связи между
интенсивностью рассеянного излучения, рассеянного в заданном направлении, и показателем ослабления [1]. Основной недостаток способа состоит в том, что необходима априорная информация о наличии корреляционного соотношения. Обычно такую информацию
можно получить, выполнив предварительно большой объем измерений [1] или вычислений [2], моделирующих рассеивающие свойства, для интересующего типа рассеивающей
среды с целью получения устойчивых, статистически обоснованных соотношений.
Наиболее близким к заявляемому является способ определения показателя ослабления, основанный на измерении коэффициента пропускания коллимированного пучка света
слоем рассеивающей среды и определении показателя ослабления с использованием закона Бугера [3] (прототип).
Недостаток данного способа состоит в наличии принципиально неустраняемых ошибок, обусловленных тем, что при регистрации пропускания приемником с конечной угловой апертурой α помимо прямо проходящего ослабленного излучения в приемник
попадает поток излучения, рассеянного в пределах апертурного угла приемника. Следовательно, из измерений пропускания получаем заниженные значения показателя ослабления.
2
BY 11892 C1 2009.04.30
Вопросу анализа возможных ошибок измерения показателя ослабления способом прототипа уделено внимание во многих экспериментальных работах применительно к разным
рассеивающим средам и условиям эксперимента. В [3-5], например, было показано, что
ошибка измерений показателя ослабления среды, состоящей из крупных частиц (больше
10 мкм) стандартными фотометрами может достигать 100 %. В общем случае рассеянный
поток, попадающий в приемник при измерениях пропускания излучения слоем частиц,
зависит от большого числа параметров, которые априори не известны. По этой причине
оценить ошибки измерения показателя ослабления базисным методом не представляется
возможным.
Для рассеивающих сред с сильно вытянутой индикатрисой рассеяния вопрос о влиянии рассеянного излучения на точность измерения показателя ослабления (в частности,
для суспензий биологических клеток) становится критическим. Результаты вычислений
согласно модели, предложенной в [6] для описания рассеивающих свойств эритроцитов,
показали, что интенсивность излучения, рассеянного образцами разбавленной крови в направлениях, близких к направлению распространения коллимированного пучка света, на
4-6 порядков больше интенсивности излучения, рассеянного в обратных направлениях
(90-180°). В такой ситуации даже при малых апертурных углах в приемник попадает значительная доля рассеянного излучения, что приводит к большим неконтролируемым
ошибкам измерения показателя ослабления.
Задачей изобретения является повышение точности определения показателя ослабления рассеивающей среды с сильно вытянутой индикатрисой рассеяния.
Точность измерения показателя ослабления зависит от многих геометрических и оптических параметров: апертурного угла приемника, геометрии эксперимента и оптических
характеристик рассеивающей среды, большинство из которых априори не известны. Исключить вклад рассеянного излучения принципиально невозможно, поэтому стремятся
минимизировать его влияние изменением геометрии эксперимента и характеристик оптической системы, в первую очередь, уменьшая угловую апертуру приемника и увеличивая
базу экспериментальной установки.
Поставленная задача решается следующим образом. В способе определения показателя ослабления рассеивающей среды с сильно вытянутой индикатрисой рассеяния берут
два образца исследуемой рассеивающей среды;
для первого образца задают величину a1 = (CV)1L1,
где (CV)1 - объемная концентрация частиц в первом образце, равная доле объема, занимаемого частицами в единице объема рассеивающей среды,
L1 - геометрическая толщина слоя первого образца;
для второго образца задают величину a2 = (CV)2L2,
где (CV)2 - объемная концентрация частиц во втором образце,
L2 - геометрическая толщина второго образца;
при освещении образцов коллимированным пучком света измеряют коэффициент
пропускания для первого образца Т1 и коэффициент пропускания для второго образца Т2;
определяют оптическую толщину τ1 рассеивающего слоя толщиной L1 из уравнения
 a
 a
T2 = exp − 2 τ1  + 2 [T1 − exp(− τ1 )] ,
 a1  a1
затем определяют показатель ослабления ε' рассеивающей среды, нормированный на единичную объемную концентрацию частиц в среде, из выражения
ε′ = τ1 / a1 = τ1 /(C V )1 L1 .
причем объемную концентрацию рассеивающих частиц в образцах и толщину образцов
выбирают с возможностью выполнения условия однократного рассеяния.
Приготавливают два образца суспензии исследуемых частиц, один из которых имеет
заданную величину a1 = (CV)1L1, где (CV)1 - объемная концентрация частиц, равная доле
3
BY 11892 C1 2009.04.30
объема, занимаемого частицами в единице объема суспензии, L1 - геометрическая толщина слоя, второй образец имеет величину a2 = (CV)2L2, где (CV)2 - объемная концентрация
частиц, L2 - геометрическая толщина слоя, измеряют коэффициенты пропускания T1 и Т2
соответствующих образцов при освещении коллимированным пучком света, определяют
оптическую толщину τ1 = a1ε' рассеивающего слоя толщиной L1 из решения уравнения
 a
 a
T2 = exp − 2 τ1  + 2 [T1 − exp(− τ1 )] , затем определяют показатель ослабления рассеи a1  a1
вающих сред, нормированный на единичную объемную концентрацию, по формуле
ε′ = τ1 / a1 = τ1 /(C V )1 L1 . При выборе значений СV и L необходимо выполнение условий
однократного рассеяния.
Предлагаемый способ измерения показателя ослабления позволяет повысить точность
определения показателя ослабления при измерениях коэффициента пропускания приемником с априори неизвестными значениями: угловой апертуры, оптических и микроструктурных параметров рассеивающей среды путем учета влияния вклада рассеянного
излучения.
Сущность предлагаемого изобретения поясняется чертежами, где:
на фиг. 1 показана принципиальная схема эксперимента по измерению показателя ослабления,
на фиг. 2 - угловое распределение интенсивности рассеянного излучения,
на фиг. 3 - доля рассеянного потока, попадающего в приемник, в зависимости от угловой апертуры приемника,
на фиг. 4 - сравнение спектральных зависимостей показателя ослабления, полученных
известным (прототип) и заявляемым способом,
на фиг. 5 - относительные отклонения значений показателя ослабления, полученных
известным (прототип) и заявляемым способом.
Способ измерения показателя ослабления осуществляют следующим образом.
Приготавливают для измерений два образца:
один образец суспензии исследуемых частиц в жидкой среде с заданной величиной
a1 = (CV)1L1, где (CV)1 - объемная концентрация частиц, равная доле объема, занимаемого
частицами в единице объема суспензии, L1 - геометрическая толщина слоя,
второй образец суспензии исследуемых частиц в жидкой среде с заданной величиной
a2 = (CV)2L2, где (CV)2 - объемная концентрация частиц, L2 - геометрическая толщина слоя,
измеряют коэффициент пропускания Т1 для первого образца, соответствующего значению a1,
затем измеряют значение коэффициента пропускания Т2 для второго образца, соответствующего значению а2,
определяют оптическую толщину τ1 = a1ε' рассеивающего слоя толщиной L1 из решения уравнения
 a
 a
T2 = exp − 2 τ1  + 2 [T1 − exp(− τ1 )] ,
 a1  a1
затем определяют показатель ослабления, нормированный на единичную объемную концентрацию, по формуле
ε′ = τ1 / a1 = τ1 /(C V )1 L1 .
Объемную концентрацию и толщину слоя выбирают таким образом, чтобы коэффициенты
пропускания T1 и T2 попадали в интервал значений 0,08-0,1 ≤ Т ≤ 0,8-0,85 и максимально
отличались по величине.
При измерении пропускания излучения слоем рассеивающей среды в приемник попадает поток излучения, представляющий собой сумму двух потоков:
F = FB + FSC,
(1)
4
BY 11892 C1 2009.04.30
FB - поток прямо проходящего коллимированного пучка света, ослабленного в соответствии с законом Бугера
(2)
FB = F0 exp(-τ) ,
FSC - поток излучения, рассеянного освещенным объемом исследуемого образца:
(3)
FSC = F0CVV∆f = F0CVLS∆f ,
F0 - поток освещающего коллимированного пучка света,
τ - оптическая толщина рассеивающего слоя
(4)
τ = CVLε' ,
L - геометрическая толщина слоя,
СV - объемная концентрация частиц, равная доле единичного объема, занимаемого
частицами в единице объема суспензии,
ε' - показатель ослабления, нормированный на единичную объемную концентрацию
частиц исследуемой рассеивающей среды,
V = SL - освещенный объем рассеивающего слоя,
S - сечение коллимированного освещающего пучка света,
∆f - величина пропорциональна интегралу от индикатрисы рассеяния по телесному углу, равному угловой апертуре приемника.
Индикатриса рассеяния - функция, описывающая угловое распределение рассеянного
излучения с длиной волны λ для системы частиц с функцией распределения частиц по
размерам f(r) и оптическими постоянными вещества частиц n(λ) и k(λ):
r2
I(β) = A ∫ i(n , k , λ, f (r ))dr = AI' (β)
(5)
r1
где n(λ) и k(λ) - значения действительной и мнимой частей комплексного показателя преломления вещества частиц по отношению к показателю преломления окружающей среды,
А - константа, включающая геометрические и оптические параметры измерительного
прибора.
На фиг. 2 показаны зависимости интенсивности рассеянного излучения от угла рассеяния I' (β) для длины волны λ = 550 нм, вычисленные для модели образца крови, разбавленного собственной плазмой. Параметры модели следующие: показатель преломления
плазмы крови npl = 1,36, действительная часть комплексного показателя преломления
эритроцитов ner = 1,41, мнимая часть комплексного показателя преломления эритроцитов
ker = 0,0083. Свойства эритроцитов аппроксимировались сферами эквивалентного объема.
Значения модального (или наиболее вероятного) радиуса задавались равными R01 = 2,285
мкм (кривая 1) и R02 = 2,965 мкм (кривая 2). Система эритроцитов описывалась гаммараспределением частиц по размерам с параметром полуширины µ = 1.
Интенсивность излучения, рассеянного суспензией эритроцитов в интервале углов
рассеяния до 30°, на 3-4 порядка больше, чем для углов 30-40° и на 5-6 порядков, чем в
направлениях 90-180°.
Поток рассеянного излучения, попадающего в апертуру приемника при измерении коэффициента пропускания слоем суспензии частиц, находится путем интегрирования интенсивности рассеянного излучения
β*
∆f = ∫ I(β) sin β dβ ,
(6)
0
β* - угол рассеяния, численно равный половине плоского угла, соответствующего
угловой апертуре приемника. Из формул (5)-(6) следует, что при заданной длине волны
освещающего излучения λ величина рассеянного потока, попадающего в апертуру приемника, является сложной функцией размеров, комплексного показателя преломления веще5
BY 11892 C1 2009.04.30
ства исследуемых частиц и показателя преломления окружающей их среды. В экспериментальных исследованиях большинство указанных параметров не известны и, следовательно, определить величину потока ∆f практически невозможно.
Для широкого круга рассеивающих сред с малым поглощением излучения можно оценить возможные ошибки определения показателя ослабления величиной ∆F, равной отношению рассеянного потока, попадающего в апертуру приемника, к потоку, рассеянному
во всех направлениях:
β*
π
0
0
δF = ∫ I(β) sin β dβ / ∫ I(β) sin β dβ .
(7)
На фиг. 3 показаны зависимости доли рассеянного потока δF (%), попадающего в приемник, от угловой апертуры приемника для модели образца разбавленной крови с параметрами, приведенными выше. Из данных фиг. 3 следует, что из измерений коэффициента
пропускания (прототип) приемником с угловой апертурой α = 4-5° доля рассеянного потока, не дающая вклада в ослабление прямо проходящего излучения и, следовательно, занижающая значения показателя ослабления, составляет 20-25 %. При угловой апертуре
α = 8-10° значения ∆F достигают 40-50 %. Относительная ошибка определения показателя
ослабления зависит от соотношения между FB и FSC (l)-(3) и увеличивается с увеличением
оптической толщины рассеивающего слоя (4).
При выполнении условий однократного рассеяния относительное угловое распределение рассеянного излучения не изменяется при изменении толщины слоя. Согласно [3], для
сильно вытянутых индикатрис рассеяния при небольших значениях оптической толщины
слоя τ ≤ 3 - 4 угловое распределение рассеянного излучения в пределах углов рассеяния
β ≤ 10° также слабо зависит от геометрической толщины слоя. Следовательно, при условиях, близких к условиям однократного рассеяния, интегральные потоки рассеянного излучения ∆f для двух образцов будут практически одинаковыми.
Тогда коэффициенты пропускания для двух образцов с различными значениями а
можно записать в виде
(8)
T1 = F1/F0 = exp(-τ1) + a1∆F ,
С учетом соотношения
T2 = F2/F0 = exp(-τ2) + a2∆F ,
(9)
∆F = 4πAS(∆f) .
(10)
a2
τ1
(11)
a1
из (8), (9) получаем уравнение относительно неизвестной величины τ1
 a
 a
T2 = exp − 2 τ1  + 2 [T1 − exp(− τ1 )] .
(12)
 a1  a1
Решая уравнение (12) методом наименьших квадратов, находим значение оптической
толщины τ1 и определяем значение показателя ослабления ε', нормированного на единичную объемную концентрацию, по формуле
ε′ = τ1 / a1 = τ1 /(C V )1 L1 .
(13)
Пример экспериментальной реализации способа.
Определение показателя ослабления образцами разбавленной крови. Образцы готовились из крови здорового человека следующим образом: из донорской крови путем центрифугирования была получена плазма, затем к 1000 мкл плазмы было добавлено 20 мкл
цельной крови. Объемная концентрация эритроцитов в цельной крови была равна 0,42
(т.е. показатель гематокрита цельной крови составлял 42 %). Объемная концентрация
τ2 =
6
BY 11892 C1 2009.04.30
эритроцитов СVer в разбавленных образцах равна отношению объема эритроцитов, содержащихся в 20 мкл цельной крови к полному объему приготовленного образца
СV = 20 × 0,42 / 1020 = 0,00824. Для измерений использовались две кюветы разной толщины L1 = 1 мм и L2 = 2 мм. Измерения коэффициента пропускания выполнены на стандартном спектрофотометре Cary 500 [7]. При измерениях использовалась кювета сравнения,
наполненная плазмой крови без эритроцитов. Результаты измерений коэффициента пропускания в интервале длин волн 860-700 нм приведены в таблице: T1 для кюветы с толщиной L1 = 1 мм - в графе 2 и T2 для кюветы с толщиной L2 = 2 мм - в графе 3.
Соответствующие значения показателя ослабления, нормированные на единичную объемную концентрацию эритроцитов, определяемые известным способом (прототип), приведены в графе 4 для L2 = 1 мм и в графе 5 для L2 = 2 мм. В графе 6 приведены значения τ1,
определяемые заявляемым способом из измерений T1 и T2 и соответствующие значения
показателя ослабления - в графе 7.
1
λ, нм
860
850
840
820
810
800
790
780
770
760
750
740
730
720
710
700
2
Т1
0,3314
0,3245
0,3150
0,3022
0,2863
0,2830
0,2821
0,2790
0,2765
0,2740
0,2709
0,2681
0,2651
0,2624
0,2599
0,2570
3
Т2
0,1474
0,1447
0,1407
0,1344
0,1264
0,1245
0,1231
0,1213
0,1195
0,1178
0,1159
0,1142
0,1125
0,1108
0,1093
0,1076
4
ε (Т1), см-1
1340
1366
1402
1452
1518
1532
1536
1549
1560
1571
1585
1598
1611
1624
1635
1649
5
ε(Т2), см-1
1162
1173
1190
1218
1255
1264
1271
1280
1289
1298
1308
1317
1326
1335
1343
1353
6
τ1, (12)
1,191
1,217
1,256
1,302
1,364
1,376
1,377
1,389
1,397
1,406
1,417
1,428
1,440
1,450
1,460
1,470
7
ε (13)
1445
1477
1524
1580
1655
1670
1671
1686
1695
1706
1720
1733
1748
1760
1772
1784
Спектральные зависимости показателей ослабления, определяемых по способу прототипа (кривая 1 из измерений пропускания при толщине слоя 1 мм и кривая 2 из измерений
при толщине слоя 2 мм) и по заявляемому способу (кривая 3), показаны на фиг. 4. Как
отмечалось выше, по способу прототипа получаем заниженные значения показателя
ослабления. При этом чем меньше толщина слоя, тем меньше помехи от рассеянного излучения, попадающего в апертуру приемника, тем ближе значения показателя ослабления
к значениям, определяемым по заявляемому способу.
На фиг. 5 приведены спектральные зависимости отклонения значений показателя ослабления, определяемых по способу прототипа, по отношению к значениям, определяемым
заявляемым способом: кривая 1 - δε1 = [ε - ε(T1)]*100/ε, кривая 2 - δε2 = [ε - ε(T2)]*100/ε. Из
данных фиг. 5 видно, что при толщине слоя 1 мм значения, определяемые по способу прототипа, занижены на 8-9 %, а при толщине слоя 2 мм - на 24-25 %.
Применение заявляемого способа увеличивает точность определения показателя ослабления рассеивающих сред и обеспечивает возможность использования стандартных
спектрофотометров с неизвестной априори угловой апертурой для измерений спектральных зависимостей показателя ослабления.
7
BY 11892 C1 2009.04.30
Источники информации:
1. Горчаков Г.И., Розенберг Г.В. Изв. АН СССР: Сер. физика атмосферы и океана, 3. № 6. - 1967. - С. 611-617.
2. Науменко Е.К., Пришивалко А.П. , О нефелометрическом методе измерения прозрачности рассеивающих сред ЖПС. - 1972. - Т. 17. - № 4. - С. 687-694.
3. Иванов А.П. Оптика рассеивающих сред. - Мн.: Наука и техника, 1969. - 592 с.
4. Аникин П.П. О вкладе рассеянного излучения в спектральное пропускание полупрозрачных облаков: Сб. Радиационные свойства перистых облаков / Под ред. Е.М. Фейгельсон. - М.: Наука, 1989. - С. 156-162.
5. Зеге Э.П. , Кацев И.Л., Полонский И.Н., Прихач А.С., Влияние рассеянного излучения на точность измерения оптической толщины тонких облаков солнечным фотометром,
Изв.АН СССР, сер.физ.атм.и океана. - 1994. - Т. 30. - № 3. - С. 328-336.
6. Науменко Е.К. Влияние агрегации эритроцитов на рассеивающие свойства крови,
Журнал прикладной спектроскопии. - Т. 70. - № 3. - С. 375-380.
7. Naumenko Е.K., Korolevich A.N. Spectral transmission by the specimens of bloodglucose mixture. Optical Technologies in Biophysics and Medicine VII, edited by V.V. Tuchin,
Proceed ings of SPIE Vol. 6163., (SPIE, Bellingham, WA, 2006) Article CID 61630Q.
Фиг. 2
Фиг. 3
8
BY 11892 C1 2009.04.30
Фиг. 4
Фиг. 5
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
9
Документ
Категория
Без категории
Просмотров
0
Размер файла
342 Кб
Теги
by11892, патент
1/--страниц
Пожаловаться на содержимое документа