RU2413930C1 - Способ определения оптических характеристик однородного рассеивающего вещества - Google Patents
Способ определения оптических характеристик однородного рассеивающего вещества Download PDFInfo
- Publication number
- RU2413930C1 RU2413930C1 RU2009147946/28A RU2009147946A RU2413930C1 RU 2413930 C1 RU2413930 C1 RU 2413930C1 RU 2009147946/28 A RU2009147946/28 A RU 2009147946/28A RU 2009147946 A RU2009147946 A RU 2009147946A RU 2413930 C1 RU2413930 C1 RU 2413930C1
- Authority
- RU
- Russia
- Prior art keywords
- peak
- sample
- optical characteristics
- coefficient
- maximum
- Prior art date
Links
Images
Abstract
Изобретение относится к области измерения оптических характеристик твердых, жидких и газообразных рассеивающих веществ и может найти применение в промышленности и медицине, в процедурах контроля качества рассеивающих веществ путем измерения их оптических характеристик, а именно путем измерения фактора анизотропии и коэффициентов рассеяния и поглощения вещества. Способ заключается в том, что образец исследуемого вещества облучают ультракороткими лазерными импульсами, регистрируют прошедшее через образец исследуемого вещества излучение, образующее бимодальное временное распределение, после чего определяют следующие оптические характеристики: коэффициент экстинкции; коэффициент поглощения; редуцированный коэффициент рассеяния; коэффициент рассеяния; фактор анизотропии, что позволяет повысить точность и информативность измерения и определения оптических характеристик рассеивающих сред. 2 з.п. ф-лы, 4 ил.
Description
Изобретение относится к области измерения оптических характеристик твердых, жидких и газообразных рассеивающих веществ. Изобретение может найти применение в промышленности и медицине, в процедурах контроля качества рассеивающих веществ путем измерения их оптических характеристик, а именно путем измерения фактора анизотропии и коэффициентов рассеяния и поглощения вещества. К области применения изобретения относятся полиграфическая, нефтяная, нефтеперерабатывающая, газовая и пищевая промышленность, включая пивоварение и молочную промышленность; промышленное и бытовое водоснабжение, а также любые области, где необходимо измерение оптических характеристик вещества.
Известен способ определения оптических характеристик однородного рассеивающего вещества по измерению ослабления лазерного излучения на различных (как минимум, двух) длинах волн и на различных (как минимум, двух) расстояниях от источника до детектора излучения [1]. Указанный способ позволяет определять такие характеристики исследуемого вещества, как коэффициент поглощения µа и редуцированный коэффициент рассеяния µs'. В то же время взаимодействие оптического излучения с рассеивающим веществом характеризуется, как минимум, тремя физическими величинами: коэффициентом поглощения µa, коэффициентом рассеяния µs и фактором анизотропии (средним косинусом угла рассеяния) g. Эти величины связаны друг с другом следующим соотношением:
Таким образом, известный способ [1] недостаточно информативен, поскольку не позволяет определять значения фактора анизотропии и коэффициента рассеяния. Кроме того, предполагается, что величина редуцированного коэффициента рассеяния и длина свободного пробега фотона не зависят от длины волны излучения. Это предположение снижает точность получаемых результатов и ограничивает область применения способа.
Известен способ определения оптических характеристик [2]. Этот способ состоит в том, что исследуемое вещество облучают оптическим излучением и по измерению ослабления непрерывного излучения для двух различных длин оптического пути определяют коэффициент ослабления. При этом для описания взаимодействия излучения с исследуемым веществом используют модифицированный закон Бугера-Ламберта-Бэра:
где I(d) - интенсивность зарегистрированного сигнала, d - длина оптического пути; I0 - интенсивность излучения источника; К - коэффициент, учитывающий степень загрязнения окошка источника излучения, чувствительность детектора излучения и другие факторы; µ - коэффициент ослабления излучения (коэффициент экстинкции).
Для измерения оптических свойств вещества проводят измерение ослабления для двух оптических длин (d и d+x), после чего вычисляют отношение интенсивностей R:
По полученному значению R вычисляют коэффициент ослабления излучения.
Основным недостатком такого способа является то, что для описания взаимодействия излучения с веществом и определения его оптических характеристик используется закон Бугера-Ламберта-Бэра, не учитывающий процесс рассеяния излучения в среде. Частично закон Бугера-Ламберта-Бэра можно использовать в рассеивающей среде, если регистрировать только баллистические фотоны, а коэффициент поглощения заменить коэффициентом экстинкции. Точность результатов, получаемых с помощью модифицированного закона Бугера-Ламберта-Бэра, падает по мере уменьшения вклада баллистических фотонов в зарегистрированный сигнал.
Известен способ определения оптических характеристик однородного рассеивающего вещества [3] (прототип), основанный на измерении временных распределений прошедших через образец ультракоротких импульсов оптического излучения и сравнении измеренных временных распределений с эталонными, рассчитанными для известных значений µa и µs' путем решения диффузионного уравнения методом Монте-Карло. Однако данный способ не предусматривает измерения фактора анизотропии исследуемой среды, что снижает информативность измерений. Кроме того, он зависит от полноты библиотеки эталонных временных распределений, что снижает его точность.
Задача изобретения - повышение точности и информативности измерения и определения оптических характеристик рассеивающих сред.
Это достигается тем, что образец исследуемого вещества облучают ультракороткими лазерными импульсами, регистрируют прошедшее через образец исследуемого вещества излучение, что измеряют значение 
максимума баллистического пика Ib(t), измеряют положение 
на временной оси, измеряют значение максимума 
рассеянного пика Is(t), после чего определяют следующие оптические характеристики: коэффициент экстинкции µ определяют из выражения 
, где 
- максимум интенсивности излучения источника, l - толщина образца исследуемого вещества; коэффициент поглощения µa и редуцированный коэффициент рассеяния µs ' определяют методом наименьших квадратов так, чтобы форма нормированного рассеянного пика 
измеренного временного распределения наилучшим образом соответствовала форме нормированного теоретического временного распределения 
, где Is(t) - рассеянный пик измеренного временного распределения, - максимум Is(t), 
- максимум It(t) 
- коэффициент диффузии, v - скорость распространения света в образце исследуемого вещества, t - время, отсчитываемое от положения максимума баллистического пика; коэффициент рассеяния µs определяют из выражения 
, фактор анизотропии g определяют из выражения 
.
Для того чтобы получить бимодальное временное распределение, подбирают соответствующую толщину образца, а именно, если при имеющейся толщине зарегистрирован только баллистический пик, толщину образца увеличивают, а если зарегистрирован только рассеянный пик, толщину образца уменьшают.
Также можно провести дополнительные измерения полуширины рассеянного пика d и расстояния между максимумами баллистического и рассеянного пиков τ, после чего коэффициент поглощения µa и редуцированный коэффициент рассеяния µs ' определяют методом наименьших квадратов так, чтобы рассчитанные значения полуширины рассеянного пика dt и расстояния между максимумами баллистического и рассеянного пиков τt наилучшим образом соответствовали измеренным.
Измеренное бимодальное временное распределение можно предварительно обработать для исключения влияния аппаратной функции, характеризующей временное разрешение регистрирующей аппаратуры, путем решения уравнения свертки. При этом в качестве аппаратной функции можно взять баллистический пик временного распределения. Тогда 
, где 
- рассеянный пик временного распределения для дальнейшей обработки, а 
- операция свертки двух функций.
Графические изображения
На фиг.1 приведено схематическое изображение процесса прохождения тонкого лазерного пучка через однородное рассеивающее вещество. 1 -баллистические фотоны, 2 - приосевые фотоны, 3, 4 - рассеянные фотоны.
На фиг.2 изображены временные распределения ультракороткого лазерного импульса, прошедшего через исследуемый образец: унимодальное распределение с регистрацией только баллистических фотонов 1 (фиг.1,а), унимодальное распределение с регистрацией только рассеянных фотонов 2 (фиг.1,б), бимодальное распределение с регистрацией баллистических 1 и рассеянных фотонов 2 (фиг.1,в). Временное распределение характеризуется максимумом интенсивности баллистического пика 
и положением этого максимума на временной оси, принимаемым за точку отсчета.
Фиг.3 иллюстрирует процесс измерения временных распределений при регистрации только баллистических фотонов - форма исходного импульса источника (а) совпадает с формой импульса, зарегистрированного детектором (б), при этом амплитуда зарегистрированного импульса уменьшается в соответствии с законом Бугера-Ламберта-Бэра.
На фиг.4 приведено идеализированное бимодальное временное распределение, 1 - положение максимума баллистического пика, 2 - рассеянный пик, τ - смещение максимума рассеянного пика относительно положения максимума баллистического пика, d - полуширина рассеянного пика.
При прохождении оптического излучения через рассеивающую среду (СРС) наблюдается следующая физическая картина (фиг.1). Часть фотонов, называемая баллистическими, проходит через среду без взаимодействия с ней. Часть фотонов поглощается и выбывает из рассмотрения. Часть фотонов, после одного или нескольких актов рассеяния, меняет направление своего движения, при этом фотоны, движущиеся вблизи оси излучения источника, называются приосевыми, а фотоны, существенно изменившие направление движения, называют внеосевыми. При этом детектор, расположенный на оси источника, может регистрировать различные типы сигнала: только баллистический пик (фиг.2,а), только рассеянный пик (фиг.2,б) или сигнал, представляющий собой суперпозицию баллистических и приосевых фотонов, так называемое бимодальное временное распределение (фиг.2,в).
Баллистическая часть бимодального временного распределения по форме повторяет импульс источника (фиг.3,а,б) и может быть описана с помощью модифицированного закона Бугера-Ламберта-Бэра. Поскольку форма баллистической компоненты совпадает с формой импульса источника, коэффициент экстинкции исследуемой среды µ=µs+µa может быть определен из следующего выражения:
Описание рассеянной части временного распределения является существенно более сложной задачей. Основным способом описания прохождения оптического излучения через биологическую среду в настоящее время является нестационарное уравнение переноса излучения (УПИ) [4-6].
УПИ в общем виде не имеет аналитического решения. В связи с этим необходимы дополнительные допущения о виде самого УПИ, виде решения, свойствах СРС и т.п., позволяющие упростить исходное выражение и получить приближенное решение. В качестве одного из вариантов такого решения может быть использована уточненная диффузионная модель переноса излучения (УДМ) [7]. В отличие от широко используемой классической диффузионной модели [4-6] УДМ учитывает исходную мононаправленность лазерного излучения, что повышает ее точность.
В УДМ для временного распределения интенсивности рассеянного пика:
Из выражения (5) и результатов измерений можно определить коэффициент поглощения µа и редуцированный коэффициент рассеяния µs'. Для этого необходимо методом наименьших квадратов минимизировать отклонение расчетных кривых временного распределения от измеренных. Поскольку форма временного распределения достаточно хорошо описывается такими параметрами, как смещение максимума рассеянного пика относительно положения максимума баллистического пика τ и полуширина рассеянного пика d (фиг.4), процедура поиска коэффициента поглощения и редуцированного коэффициента рассеяния может заключаться в минимизации отклонения расчетных значений dt и τt от измеренных путем подбора соответствующих значений µа и µs'.
Коэффициент рассеяния исследуемой среды определяют из выражения:
Фактор анизотропии g определяют по полученным значениям коэффициента рассеяния и редуцированного коэффициента рассеяния из выражения
Для повышения точности измерений и более качественного разделения баллистической и рассеянной составляющих экспериментально полученного временного распределения можно использовать процедуру коррекции влияния аппаратной функции. Суть такой коррекции заключается в следующем. Полученное экспериментально временное распределение I(t) можно представить как результат свертки идеального временного распределения I(0)(t) с некоторым ядром h(t), которое принято называть аппаратной функцией:
Для уменьшения влияния аппаратной функции необходимо, во-первых, найти аппаратную функцию, а во-вторых, решить свертку (8) относительно I(t). Вторая задача решается известным способом с помощью преобразования Фурье и теоремы о свертке. При этом решение 
свертки (8) можно представить в виде:
где ℑ1{•} - преобразование Фурье.
Для бимодального временного распределения в качестве аппаратной функции h(t) можно выбрать часть временного распределения, соответствующую баллистическим фотонам Ib[t).
Предлагаемый способ определения оптических характеристик однородного рассеивающего вещества основан на использовании существующих методов измерения оптического излучения, прошедшего через рассеивающие среды, и математической обработки данных. Таким образом, реализация предложенного способа определения оптических характеристик рассеивающих сред не требует разработки новых методов обработки данных или новых технических устройств и может быть осуществлена в настоящее время. Повышение информативности определения оптических характеристик рассеивающих сред позволит существенно повысить эффективность методов контроля качества жидких и газообразных рассеивающих сред.
Источники информации
1. Патент США 5,529,065, МПК A61B 500.
2. Патент США 4,981,362, МПК G01N 2100; G01N 1506; G01J 132.
3. Патент США 5,867,807, МПК G06F 1900.
4. Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах. - М., Мир, 1981. - T.1.
5. Кейз К., Цвайфель П. Линейная теория переноса. - М., Мир, 1972.
6. Терещенко С.А. Методы вычислительной томографии. - М., Физматлит, 2004.
7. Данилов А.А., Терещенко С.А., Подгаецкий В.М. Уточненная диффузионная модель для описания взаимодействия лазерного излучения с биологической тканью. // Оптика и спектроскопия, 2007. - Т.102, №5. - С.849-854.
Claims (3)
1. Способ определения оптических характеристик однородного рассеивающего вещества, включающий облучение образца исследуемого вещества ультракороткими лазерными импульсами, регистрацию прошедшего через образец исследуемого вещества излучения и определение оптических характеристик однородного рассеивающего вещества, отличающийся тем, что измеряют значение максимума баллистического пика Ib(t), измеряют положение на временной оси, измеряют значение максимума рассеянного пика Is(t), после чего определяют следующие оптические характеристики: коэффициент экстинкции µ определяют из выражения
где - максимум интенсивности излучения источника;
l- толщина образца исследуемого вещества,
коэффициент поглощения µa и редуцированный коэффициент рассеяния µs' определяют методом наименьших квадратов так, чтобы форма нормированного рассеянного пика измеренного бимодального временного распределения наилучшим образом соответствовала форме нормированного теоретического временного распределения
где Is(t) - рассеянный пик бимодального измеренного временного распределения;
- максимум Is(t);
- максимум It(t);
где - коэффициент диффузии;
v - скорость распространения света в образце исследуемого вещества;
t - время, отсчитываемое от положения максимума баллистического пика;
коэффициент рассеяния µs определяют из выражения µs=µ-µа;
фактор анизотропии g определяют из выражения
где
l- толщина образца исследуемого вещества,
коэффициент поглощения µa и редуцированный коэффициент рассеяния µs' определяют методом наименьших квадратов так, чтобы форма нормированного рассеянного пика
где Is(t) - рассеянный пик бимодального измеренного временного распределения;
где
v - скорость распространения света в образце исследуемого вещества;
t - время, отсчитываемое от положения максимума баллистического пика;
коэффициент рассеяния µs определяют из выражения µs=µ-µа;
фактор анизотропии g определяют из выражения
2. Способ определения оптических характеристик однородного рассеивающего вещества по п.1, отличающийся тем, что для получения бимодального временного распределения подбирают соответствующую толщину образца, а именно, если при имеющейся толщине зарегистрирован только баллистический пик, толщину образца увеличивают, а если зарегистрирован только рассеянный пик, толщину образца уменьшают.
3. Способ определения оптических характеристик однородного рассеивающего вещества по п.1, отличающийся тем, что измеренное бимодальное временное распределение предварительно обрабатывают для исключения влияния аппаратной функции путем решения уравнения свертки 
где 
- рассеянный пик временного распределения для дальнейшей обработки.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2009147946/28A RU2413930C1 (ru) | 2009-12-24 | 2009-12-24 | Способ определения оптических характеристик однородного рассеивающего вещества |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2009147946/28A RU2413930C1 (ru) | 2009-12-24 | 2009-12-24 | Способ определения оптических характеристик однородного рассеивающего вещества |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2413930C1 true RU2413930C1 (ru) | 2011-03-10 |
Family
ID=46311202
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2009147946/28A RU2413930C1 (ru) | 2009-12-24 | 2009-12-24 | Способ определения оптических характеристик однородного рассеивающего вещества |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2413930C1 (ru) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2533538C1 (ru) * | 2013-08-19 | 2014-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тверской государственный университет" | Способ раздельного определения вероятностей поглощения и рассеяния фотонов на единицу пути в твердых оптических материалах |
| RU2700722C1 (ru) * | 2018-11-06 | 2019-09-19 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МО РФ | Способ исследований температурных зависимостей оптических характеристик полупроводниковых материалов |
-
2009
- 2009-12-24 RU RU2009147946/28A patent/RU2413930C1/ru not_active IP Right Cessation
Non-Patent Citations (1)
| Title |
|---|
| ДАНИЛОВ А.А. и др. Экспериментальное определение коэффициентов рассеяния и поглощения излучения в однородном слое сильнорассеивающей биологической среды. Медицинская техника, №4, 2006, с.17-20. * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2533538C1 (ru) * | 2013-08-19 | 2014-11-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Тверской государственный университет" | Способ раздельного определения вероятностей поглощения и рассеяния фотонов на единицу пути в твердых оптических материалах |
| RU2700722C1 (ru) * | 2018-11-06 | 2019-09-19 | Федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия Ракетных войск стратегического назначения имени Петра Великого" МО РФ | Способ исследований температурных зависимостей оптических характеристик полупроводниковых материалов |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3171160B1 (en) | Apparatus and method for detecting microbes or bacteria | |
| US20160309068A1 (en) | Spatial frequency domain imaging using custom patterns | |
| CA2988263C (en) | Nanoparticle analyzer | |
| JP6775660B2 (ja) | 多位置拡散スペクトルデータの処理、モデリング、予測方法および処理装置 | |
| Duadi et al. | Dependence of light scattering profile in tissue on blood vessel diameter and distribution: a computer simulation study | |
| CN105277497A (zh) | 连续光谱下检测农产品的光学特性参数的装置与方法 | |
| JPWO2006085403A1 (ja) | 実時間テラヘルツ・トモグラフィー装置および分光イメージング装置 | |
| JPWO2020105682A1 (ja) | 血中脂質濃度計測装置及びその方法 | |
| ElMasry et al. | Prediction of meat spectral patterns based on optical properties and concentrations of the major constituents | |
| Varkentin et al. | Simple model to simulate oct-depth signal in weakly and strongly scattering homogeneous media | |
| Singh et al. | Discriminating turbid media by scatterer size and scattering coefficient using backscattered linearly and circularly polarized light | |
| CN105334204B (zh) | 一种基于傅里叶变换拉曼光谱分析方法 | |
| EP3797687A1 (en) | Blood vessel detection device and method therefor | |
| RU2413930C1 (ru) | Способ определения оптических характеристик однородного рассеивающего вещества | |
| JP7077175B2 (ja) | 自動分析装置、自動分析方法、および、プログラム | |
| WO2018151150A1 (ja) | 散乱体濃度計測装置及びその方法 | |
| Borozdova et al. | Laser Doppler anemometer signal processing for blood flow velocity measurements | |
| Singh et al. | Spatial helicity response metric to quantify particle size and turbidity of heterogeneous media through circular polarization imaging | |
| JP3902999B2 (ja) | 光学的散乱特性推定装置およびその作動方法 | |
| JP2013217909A (ja) | 屈折率算出方法及び装置、屈折率算出用物質、及びトモグラフィ装置 | |
| Akbarzadeh et al. | Experimental validation of a spectroscopic Monte Carlo light transport simulation technique and Raman scattering depth sensing analysis in biological tissue | |
| WO2019123559A1 (ja) | 脂質計測装置及びその方法 | |
| Sato et al. | Estimating the absorption coefficient of the bottom layer in four-layered turbid mediums based on the time-domain depth sensitivity of near-infrared light reflectance | |
| Di Rocco et al. | CW laser transilluminance in turbid media with cylindrical inclusions | |
| Wang et al. | A sequential method for estimating the optical properties of two-layer agro-products from spatially-resolved diffuse reflectance: Simulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20161225 |