способ определения пространственного распределения физических параметров газовой среды

Формула изобретения

1. Способ определения пространственного распределения физических параметров газовой среды путем зондирования электромагнитным излучением с фиксированной длиной волны исследуемой области среды, измерения интенсивности сигналов взаимодействия излучения со средой на второй длине волны, направление приема сигналов взаимодействия выбирают под определенным углом к направлению зондирования, проведения приема сигналов взаимодействия из заданного количества точек исследуемой области среды, отличающийся тем, что осуществляют зондирование среды в прямом и обратном направлениях, измеряют интенсивности сигналов взаимодействия излучения со средой на третьей длине волны для обоих направлений зондирования, осуществляют зондирование среды на первой длине волны в поперечном направлении в плоскости, нормальной к плоскости, образованной зондирующей и приемными оптическими осями при сохранении прежних направлений приема и измерении интенсивностей сигналов на третьей длине волны, пространственное распределение физических параметров среды определяют по формуле



где _n - физический параметр в n-й точке сечения исследуемой среды;

L - коэффициент пропорциональности;

P^п_n^р(₂) - интенсивность измеренного сигнала из n-й точки на второй длине волны ₂ при зондировании в прямом направлении;

P^п_n^р(₃) и P^о_n^бр(₃) - интенсивности измеренного сигнала из n-й точки на третьей длине волны ₃ в прямом и обратном направлениях соответственно;

P_n(₃) - интенсивность измеренного сигнала в n-й точке сечения исследуемой среды на второй длине волны ₃ при зондировании в поперечном направлении в плоскости, нормальной к плоскости, образованной зондирующей и приемными оптическими осями.

2. Способ по п.1, отличающийся тем, что принимают сигналы взаимодействия излучения со средой на длине волны _3, равной длине волны зондирующего излучения _1.

3. Способ по п.1, отличающийся тем, что принимают сигналы взаимодействия излучения со средой на длинах волн ₂ и _3, равных длине волны зондирующего излучения ₁.

Описание изобретения к патенту

Изобретение относится к измерительной технике для диагностики атмосферы.

Известен способ определения физических параметров газовой среды путем облучения ее пучком излучения и последующего выделения и регистрации сигнала, обусловленного взаимодействием зондирующего пучка и исследуемой среды (а.с. СССР N 1288561, МКИ кл. G 01 N 21/64, 1984). Недостатком способа является низкая точность определения в условиях оптической неоднородности исследуемой среды.

Наиболее близким к предлагаемому является способ определения пространственного распределения физических параметров газовой среды, заключающийся в облучении его пучком излучения, выделении сигнала, обусловленного взаимодействием зондирующего пучка и исследуемой среды, и регистрации пространственного распределения этого сигнала (а.с. СССР N 1665285, МКИ кл. G 01 N 21/64, 1991). Оптическое излучение просвечивает объект исследования, результатом взаимодействия излучения и исследуемой среды могут быть лазерно-индуцированная флуоресценция, комбинационное рассеяние, также являющиеся излучением. Выделение сигнала, обусловленного этим взаимодействием для разных областей воздушной среды, позволяет зарегистрировать его пространственное распределение.

Недостатком этого способа является низкая точность определения в условиях неоднородности объекта измерений, поскольку неравномерность оптической плотности воздушной среды при зондировании разных точек выбранного сечения непосредственно влияет на результат измерения.

Решаемая техническая задача - повышение точности определения пространственного распределения физических параметров газовой среды.

Решаемая техническая задача в способе определения пространственного распределения физических параметров газовой среды, заключающемся в зондировании электромагнитным излучением с фиксированной длиной волны исследуемой области среды, измерении интенсивности сигналов взаимодействия излучения со средой на второй длине волны, направление приема сигналов взаимодействия выбирают под углом к оси зондирования, проведении приема сигнала взаимодействия из заданного количества точек исследуемой области среды, достигается тем, что зондирование среды осуществляют в прямом и обратном направлениях, измеряют интенсивности сигналов взаимодействия излучения со средой на третьей длине волны для обоих направлений зондирования, осуществляют зондирование среды на первой длине волны в поперечном направлении в плоскости, нормальной к плоскости, образованной зондирующей и приемными оптическими осями при сохранении прежних направлений приема, и измерении интенсивности сигналов взаимодействия излучения со средой на третьей длине волны, пространственное распределение физических параметров среды определяют по формуле:



где _n - физический параметр в n-й точке сечения исследуемой среды; L - коэффициент пропорциональности; P^п_n^р(₂) - интенсивность сигнала из n-й точки среды на второй длине волны ₂ при зондировании в прямом направлении; P^п_n^р(₃) и P^о_n^бр(₃) - интенсивности сигнала из n-й точки среды на третьей длине волны ₃ в прямом и обратном направлениях соответственно, интенсивность сигнала из n-й точки сечения исследуемой среды на третьей длине волны ₃ при зондировании в плоскости, ортогональной к плоскости основного измерения, образованной оптической осью основного зондирования и приемной оптической осью.

Возможен прием сигналов взаимодействия излучения со средой на длине волны ₃, равной длине волны зондирующего излучения ₁.

Возможен прием сигналов взаимодействия излучения со средой на длинах волн ₂ и ₃, равных длине волны зондирующего излучения ₁.

На чертеже изображена блок-схема устройства, реализующего способ.

Устройство содержит: блок 1 управления, блоки передающей оптики 2 и 3 с диаграммами направленности, расположенными вдоль одной оптической оси АВ; формирователь 4 плоского "веерообразного" пучка с диаграммой направленности, охватывающей все исследуемые точки среды в поперечном направлении в плоскости, ортогональной плоскости, образованной осью зондирования АВ и приемными оптическими осями; линии (каналы) передачи зондирующего излучения 5, 6 и 7; оптический коммутатор 8; излучатель 9 с источником накачки; исследуемую среду 10, приемный объектив 11; блок оптических фильтров 12, включающий селективный элемент 13 на длину волны излучения ₂, обусловленного полезным эффектом взаимодействия зондирующего пучка и объекта исследования, селективный элемент 14 на длину волны ₃ другого взаимодействия и блок замены селективных элементов 15; координатно-чувствительное фотоприемное устройство 16 с блоком предварительной обработки; вычислитель 17.

Излучатель 9 электрически соединен с блоком 1 управления и оптически связан через оптический коммутатор 8 с линиями передачи светового излучения 5, 6 и 7 и блоками передающей оптики 2, 3 и формирователем веерообразного пучка 4. Диаграммы направленности блоков передающей оптики 2 и 3 расположены вдоль одной оптической оси, на которой измеряется пространственное распределение физического параметра, например концентрации молекул объекта 10. Координатно-чувствительное фотоприемное устройство 16 служит для приема сигналов взаимодействия излучения с исследуемой средой.

Координатно-чувствительное фотоприемное устройство 16, селективные элементы 13 и 14 оптически связаны с приемным объективом 11. Блок 1 управления соединен с управляющими входами излучателя 9, снабженного источником накачки, оптического коммутатора 8, управляющим входом координатно-чувствительного фотоприемного устройства 16, блоком 15 замены селективных элементов, управляющим входом вычислителя 17. Выходы координатно- чувствительного фотоприемного устройства 16 связаны с информационными входами вычислителя 17.

Способ осуществляют с помощью устройства следующим образом.

Блок 1 управления запускает лазер 9 с источником накачки. Импульс зондирующего излучения на длине волны ₁, пройдя через оптический коммутатор 8, линию передачи светового излучения 5 и блок передающей оптики 2, направляется из точки А в подлежащую исследованию область 10 газовой среды. При взаимодействии излучения со средой могут иметь место явления лазерно-индуцированной флуоресценции, комбинационного рассеяния, поглощения и другие. Результат взаимодействия несет в себе информацию о таких параметрах исследуемой среды, как концентрация частиц, их скорость, давление и т.д., и также является излучением. Часть этого излучения воздушной среды, являющегося результатом взаимодействия, например из точки M₁ собирается объективом 11, выделяется селективным элементом 13 (например, интерференционным фильтром) на длине волны ₂ и регистрируется координатно-чувствительным фотоприемным устройством 16. После этого по сигналу с блока управления 1 система 15 поворота блока селективных элементов устанавливает перед блоком 16 селективный элемент 14 на длину волны ₃. Прием излучения из той же точки M₁ ведется на длине волны ₃.

Далее блок управления 1 при помощи оптического коммутатора 8 переключает излучение лазера 9 на линию передачи светового излучения 6. Зондирующее излучение через блок передающей оптики 3 направляется в исследуемую область 10 из точки В. Координатно-чувствительное фотоприемное устройство 16 регистрирует сигнал на длине волны ₃ из той же точки M₁, но направление облучения среды противоположно рассмотренному ранее. Далее блок управления 1 переключает излучение при помощи оптического коммутатора 8 на линию передачи излучения 7, при помощи передающей оптики формируется луч просвечивающий исследуемые точки пространства в поперечном направлении в плоскости, нормальной к плоскости, образованной зондирующей и приемными оптическими осями. Прием излучения ведется из той же точки M₁ на длине волны ₃. Сигналы с выхода координатно-чувствительного фотоприемного устройства 16 подаются на вычислитель 17 для определения искомой характеристики в точке M₁ объекта исследования.

После этого по сигналу с блока 1 координатно-чувствительное фотоприемное устройство 16 перестраивается на прием излучения из новой точки исследуемой области воздушной среды, например из точки М₂, затем из точки М₃ и т.д. Число точек в данном случае определяется количеством фотоэлементов в координатно-чувствительном фотоприемном устройстве (например, число исследуемых точек равно 10). При выполнении условия малости расстояния между точками M₁ и M₂ и между M₂ и М₃ по сравнению с расстояниями M₁C, М₂С и М₃С можно считать угол постоянным для всех точек М₁, М₂ и М₃ (см. чертеж).

В результате снимается пространственное распределение физического параметра (например, концентрации молекул какого-либо газа) в выбранном сечении среды.

Оптический коммутатор 8 может быть реализован на основе механического дефлектора дискретного отклонения с системой поворотных зеркал либо электрооптического дефлектора. "Веерообразный" пучок излучения может быть получен путем сканирования лазерного луча в выбранной плоскости зондирования при помощи дефлектора непрерывного отклонения.

Покажем, что предлагаемый способ позволяет достичь поставленной решаемой технической задачи - повысить точность измерения пространственного распределения физических параметров исследуемой среды.

При определении пространственного распределения физических параметров исследуемой среды (например, концентрации газа) по результатам измерения интенсивности сигналов взаимодействия зондирующего излучения с исследуемой средой, как это делается в прототипе, значительной оказывается погрешность, вносимая различным ослаблением сигналов в направлении измерения пространственного распределения от точек M_i до точки А и до точки С. Пренебрежение этой погрешностью может сильно исказить измеряемое распределение.

Рассмотрим возможность повышения точности измерений на примере комбинационного рассеяния и флуоресценции.

Особенностью метода комбинационного рассеяния является то, что длина волны лазерного излучения не связана жесткими требованиями резонанса частоты излучения лазера с частотами переходов атомов и молекул. Это позволяет использовать один источник накачки для одновременной регистрации большого числа молекулярных примесей, выделяя рассеянные стоксовые (иногда - антистоковые) компоненты на различных частотах ^c_i = _I-^м_i(^а_i^с= ₁+^м_i - для антистоксовых), где ₁ - частота лазерного излучения; ^м_i - частота, соответствующая колебательному кванту молекулы i-й примеси.

Как концентрации N, так и (в значительно меньшей степени) сечения комбинационного рассеяния ^кр различных компонент могут сильно различаться. Если для двух молекулярных компонент справедливо N₂^к₂^р N₃^к₃^р, то принимаемый на длине волны ₂ сигнал комбинационного рассеяния, соответствующий частоте ^c₂ = ₁-^м₂, может быть во много раз меньше, чем сигнал комбинационного рассеяния, принятый на длине волны ₃, соответствующей частоте ^c₃ = ₁-^м₃. В силу этого относительная погрешность измерений на ₃ оказывается значительно ниже, чем на ₂.

Поэтому если для измерений пространственного распределения физического параметра исследуемой среды по результатам измерений интенсивности сигнала комбинационного рассеяния на длине волны ₂ суметь ввести корректирующие коэффициенты, устраняющие влияние различного ослабления сигналов от точек М_i в направлении измерения пространственного распределения, которые получены по результатам более точных измерений интенсивности комбинационного рассеяния на длине волны ₃, то можно повысить точность измерений пространственного распределения первой компоненты. При этом длины волн ₂ и ₃ должны отличаться незначительно, чтобы не проявилась возможная спектральная изменчивость оптической плотности исследуемой среды. В то же время спектральное разрешение измерительной аппаратуры должно обеспечивать возможность уверенного раздельного измерения сигналов на этих длинах волн.

Приводить количественные соотношения для комбинационного рассеяния не будем - рассмотрим их ниже на примере флуоресценции. Однако следует иметь в виду, что они полностью применимы для комбинационного рассеяния, в случае которого ₁,₂,₃ - различны, тогда как для флуоресценции удобно принять ₃= ₁.

Метод лазерной флуоресценции основан на регистрации спонтанного излучения молекул и атомов примеси, подвергнутых действию излучения, длина волны которого находится в резонансе с одним из разрешенных переходов с переводом молекул в возбужденное состояние. При этом релаксация в основное состояние происходит излучательно.

Пусть методом лазерно-индуцированной флуоресценции измеряется относительный профиль распределения концентрации исследуемых молекул в точках M₁, M₂, М₃ (см. чертеж), расположенных на одной прямой; длина волны зондирующего излучения ₁, длина волны флуоресценции ₂.

Прозрачности среды от приемника излучения (от точки С) до точек M₁, М₂, М₃ равны соответственно T(CM₁), Т(СМ₂), Т(СМ₃).

При приеме излучения из точки M₁, три сигнала флуоресценции P^п₁^р(₂), рассеянный сигнал P^п₁^р(₁) при облучении среды в прямом направлении, а также рассеянный сигнал P^о₁^бр(₁) при облучении среды в обратном направлении (см. Приложение, формулы (7п) и (8п)) равны соответственно:

P^п₁^р(₂) = kT(CM₁)T(M₁A)₁ (1)

P^п₁^р(₁) = kT(CM₁)T(M₁A)⁰₁i() (2)



Здесь и далее понимается: _n= ^ф_n (для комбинационного рассеяния - ^к_n^р1 - характеризует основное взаимодействие; ⁰_n = ^p_n (для комбинационного рассеяния - ^к_n^р2) - вспомогательное.

Аналогично при приеме сигналов из точек M₁ и М₂

P^п₂^р(₂) = kT(CM₂)T(M₁M₂)T(M₁A)₂ (4)

P^п₂^р(₁) = kT(CM₂)T(M₁M₂)T(M₁A)⁰₂i() (5)







P^о₃^бр(₁) = kT(CM₃)T(M₃B)⁰₃i(-) (9)

Определение относительного распределения концентрации N по измерениям только сигналов P^п₁^р(₂),P^п₂^р(₂),P^п₃^р(₂) приведет к значительным погрешностям, которые определяются различным ослаблением сигналов от точек M₁, М₂, М₃ до точки А и точки С.

₁:₂:₃= P^п₁^р(₂):P^п₂^р(₂):P^п₃^р(₂) (10)

Где = N^ф, N - концентрация молекул; ^ф - сечение флуоресценции. Чем больше отличаются обычно не учитываемые оптические плотности между точками M₁, М₂, М₃ и точкой A (T(M₁A), Т(М₂А), Т(М₃А)) и точкой С (T(CM₁), Т(СМ₂), Т(СМ₃)), тем ниже точность измерения пространственного распределения по результатам измерения величин P^п_n^р(₂), n = 1, 2, 3,...

Формула (10) характеризует пространственное распределение физических параметров исследуемой среды, например концентрации какой-либо молекулярной компоненты, получаемое прототипом и другими традиционными методами.

Из (1) и (2), (4) и (5), (7) и (8) следует:

_n= P^п_n^р(₂)⁰_ni()/P^п_п^р(₁), n = 1,2,3,... (11)

Из (2) и (3) вытекает, что произведение принятых сигналов P^п_n^р(₁) и P^о_n^бр(₁) равно:



Аналогичные выражения можно получить из (5) и (6), (8) и (9).

В общем виде:



Здесь: Т(АВ) - прозрачность среды на всем участке зондирования от точки А до точки В.

Из (12) получаем:



Из (13) с учетом (11)



Для устранения влияния ослаблений среды T(CM₁), Т(СМ₂), Т(СМ₃) нами производится зондирование среды в плоскости, ортогональной плоскости основного измерения, образованной зондирующей (АВ) и приемными (СМ_i) оптическими осями.

Будем считать, что прозрачности среды (см. чертеж) от точки D излучателя до точек M₁, M₂, М₃ (T(DM₁), T(DM₂), T(DM₃)) отличаются незначительно и равны Т.

При приеме излучения из точек M₁, М₂, М₃ сигнал флуоресценции на = ₁:



Откуда



Из (14) и (16) получаем:

(17)

Делая допущение, что ⁰₁ = ⁰₂ = ⁰₃ вдоль исследуемого направления, получаем соотношение между параметрами _n в разных точках



Формула (18) описывает пространственное распределение физического параметра. Параметр _n в произвольной n-й точке сечения:



L - коэффициент пропорциональности. Выражение (19) сразу дает профиль концентрации, поскольку N = /.

Прием и регистрация, кроме флуоресценции, рассеянных сигналов на длине волны ₁ зондирующего излучения позволяют значительно повысить точность измерения. Хотя _{1 2}, но это отличие незначительно, так что оптические характеристики (в частности, прозрачность) близки с высокой точностью. В то же время сечение флуоресценции ^ф мало по сравнению с сечением упругого рассеяния на частицах. При изучении воздушной среды различные частицы, находящиеся в исследуемом сечении, способны весьма существенно ослабить (за счет рассеяния и поглощения) возбужденный лазером сигнал флуоресценции. Поэтому в силу ^ф измерение на = ₁ рассеянных сигналов от облучения среды с двух сторон, а так же зондирование среды в плоскости, ортогональной плоскости основного измерения, образованной зондирующей (АВ) и приемными (СМ_i) оптическими осями, позволяет получить более высокую точность измерения соотношения прозрачностей участков между исследуемыми точками сечения, так как значительно возрастает по сравнению с флуоресценцией отношение сигнал/шум.

Ход приведенных выше рассуждений и окончательные формулы справедливы и для определения распределения физических параметров среды методом комбинационного рассеяния, которое регистрируется на ₂ и ₃. Если, как отмечалось ранее, неравенство N₁^к₁^р N₂^к₂^р выполняется довольно сильно, то точность измерения физического параметра (например, концентрации молекул N) в условиях неоднородной среды может быть существенно повышена. Для количественного описания пространственного распределения, измеренного методом комбинационного рассеяния, в формуле (19) под P^п_n^р(₂) следует понимать интенсивность принятого сигнала на ₂, а под интенсивности сигналов на ₃, полученных при облучении среды в прямом, обратном направлениях и в поперечном направлении в плоскости, ортогональной плоскости основного измерения.

Поэтому в данном случае



Оценим выигрыш в точности измерений по предлагаемому способу на примере использования эффекта комбинационного рассеяния. Пусть концентрация измеряемых молекул в точках M₁, M₂, М₃, а значит и параметры ₁,₂,₃ относятся как 1: 1: 1. Предположим, что излучение ослабляется средой, и прозрачность среды Т между точками M₁ и М₂, а также между точками М₂ и М₃ одинакова и равна, например, 0,8. T(M₁M₂) = T(M₂М₃) = 0,8. Ослабление излучения на пути М₁А и М₃В (соответственно Т(М₁А) и Т(М₃В)) учитывать не будем, так как при измерении относительного профиля пространственного распределения T(M₁A) и Т(М₃В) сокращаются. Прозрачности среды от приемника излучения (точка С) до точек M₁, М₂, М₃ равны соответственно T(CM₁) = 0,9; Т(СМ₂) = 0,8; Т(СМ₃) = 0,7.

Без учета аппаратурных и геометрических констант, влияние которых легко учесть, сигнал комбинационного рассеяния на = ₂ из точки М₁ при облучении среды в прямом направлении из точки А

P^п₁^р(₂) = 1T(CM₁) = 0,9

Аналогично





Относительный профиль распределения концентрации молекул, полученный по традиционной методике, повторяет распределение принятых сигналов Р^пр_n



Как видно, это распределение сильно искажено по сравнению с исходным 1: 1:1.

Согласно предлагаемому способу измеряются и сигналы комбинационного рассеяния на длине волны ₃ из тех же точек M₁, М₂, М₃ при облучении среды в прямом направлении (от точки А), в обратном (от точки В), а также зондирование среды в плоскости, ортогональной плоскости основного измерения, образованной зондирующей (АВ) и приемными (CM_i) оптическими осями:

P^п₁^р(₃) = 0,9;





P^п₂^р(₃) = 1T(M₁M₂)T(CM₂) = 0,64;











Согласно предложенному алгоритму



Поэтому



Таким образом, исходное распределение восстановлено правильно. Предлагаемый способ определения пространственного распределения физического параметра устраняет погрешность традиционных методов, связанную с ослаблением излучения в исследуемой среде.

Изобретение можно проиллюстрировать следующими примерами, в которых использован принцип комбинационного рассеяния: зондирующий лазер - аргоновый: ₁ - 0,5145 мкм, молекулы с близкими комбинационными сдвигами - СО (комбинационный сдвиг - 2145 см^-1, _кр= 0,5713 мкм) и N₂ (комбинационный сдвиг - 2331 см ^-1, _кр= 0,5762 мкм); другая возможная пара молекул: SO₂ (комбинационный сдвиг - 1151 см^-1, _кр= 0,5450 мкм) и N₂O (комбинационный сдвиг - 1290 см^-1, _кр- 0,5486 мкм).

Случай, когда длины волн ₁,₂ и ₃ различны, отвечает использованию метода комбинационного рассеяния для измерения пространственного распределения физических параметров.

Случай, когда ₃= ₁ соответствует использованию лазерно-индуцированной флуоресценции:,



Случай, когда ₁= ₂= ₃ реализуется при измерении упругого рассеяния на длине волны зондирующего лазера ₁. В этом случае и _n= L(P^о_n^бр(₁)P^п_n^р(₁))^1/2/P_n(₁).$