Способ определения коэффициента ослабления рассеивающей среды

 

Использование: изобретение относится к метеорологии и гидрологии, охраны окружающей среды и может быть использовано для определения оптических характеристик рассеивающих сред. Сущность изобретения: повышение точности определения коэффициентов ослабления в средах с резкими границами раздела путем коррекции сигнала обратного рассеяния на величину перепада лидарного отношения. Сигналы обратного рассеяния дополнительно накапливают на двух равных участках до границы раздела и двух равных участках после границы раздела, по измеренным величинам сигналов до и после границы раздела определяют коэффициент коррекции сигналов обратного рассеяния от участков среды после границы раздела и о величине коэффициента ослабления/ судят по значениям сигналов обратного рассеяния в точке определения Д накопленным значениям сигналов от калибровочной точки до точки определенияД опорным значениям/ и зна-. чениям коэффициента коррекции сигналов от участков после границы раздела. 4 ил. ел с

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК

ГОСУДАРСТВЕННОЕ ПАТЕНТНОЕ

ВЕДОМСТВО СССР (ГОСПАТЕНТ СССР) ОПИСАНИЕ ИЗОБРЕТЕНИЯ

К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ (21) 4844462/25 (22) 29,05,90 (46) 15.04.93, Бюл. N. 14 (71) Белорусский государственный университет. им. В.И.Ленина (72) М.M.Êóãåéêî, И.А,Малевич и Д,З.Шиперко (56) Авторское свидетельство СССР

N. 390401, кл, G 01 W 1/00, 1973.

Креков Г.М., Кавкянов С.И„Крекова

М.М. Интерпретация сигналов оптического зондирования атмосферы. Новосибирск, Наука, 1987, с, 108.

Зуев В.Е., Креков Г,М., Крекова IVl,M., Наац И.Э, Теоретические аспекты проблемы лазерного зондирования облаков. В кн„

Вопросы лазерного зондирования атмосферы. Новосибирск, Наука, 1976: с. 3-33, (54) СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ОСЛАБЛЕНИЯ РАССЕИВАЮЩЕЙ

СРЕДЫ (57) Использование: изобретение относится к метеорологии и гидрологии, охрайы окруИзобретение относится к области метеорологии и гидрологии, охраны окружающей среды и может быть. использовано для определения оптических характеристик рассеивающих сред, Наиболее близким к предлагаемому является способ, заключающийся в посылке импульсов света малой длительности, пре) образовании рассеянного в обратном направлении света в электрические сигналы, накопление этих сигналов в течение. заданного времени, в зависимости QT местоположения калибровочной точки и точки определения коэффициента ослабления,,, Я2 „1809408 А1 (я)я G 01 W 1/00, G 01 N 21/47 жающей среды и может быть использовано для определения оптических характеристик рассеивающих сред. Сущность изобретения: повышение точности определения коэффициентов ослабления в средах с резкими границами раздела путем коррекции сигнала обратного рассеяния на величину перепада лидэрного отношения.

Сигналы обратного рассеяния дополнительно накапливают на двух равных участках до границы раздела и двух равных участках после границы раздела, по измеренныМ величинам сигналов до и после границы раздела определяют коэффициент коррекции сигналов обратного рассеяния от участков среды после границы раздела и о величине коэффициента ослабления Р судят по значениям сигналов обратного рассеяния в точке определения Р, накопленным значениям сигналов от калибровочной точки до точки определенияР, опорным значениямP и значениям коэффициента коррекции сигналов от участков после границы раздела. 4 ил. при увеличении усиления принятых сигналов пропорционально квадрату текущего времени, отсчитываемого с момента посылки зондирующего импульса и определении независимым способом калибровочного (опорного значения коэффициента ослабле- . ния, Недостатком данного способа является то, что реализация данного способа требует постоянства лидарного отнэшения по трассе зондирования, Более того, при зондировании по наклонным трассам, а также на большом удалении от лидара практически невозможно получение значения коэффици1809408 ента ослабления в удаленной от лидара калибровочной (опорной) точке.

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

Цель изобретения — повышение точности определения коэффициента ослабления в средах с резкими границами раздела путем учета изменения лидарного отношения на границах раздела.

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

Производится также регистрация сигнала обратного рассеяния из исследуемой точки среды..

Для величины сигнала обратного рассеяния в исследуемой точке Zj можно записать следующее выражение:

Р(21) = А21 2Т (0;Zo) р (2 )Т (Zo, Zg), (1) где Р(2 ) — регистрируемый сигнал обратного рассеяния в исследуемой точке Z ; А — аппаратурная константа лидарэ, Pz(2 ) =

= g (2)P(Zj) — коэффициент обратного рассеяния в исследуемой точке Z; g(2 )ô {2 )— соответственно лидарное отношение и коэффициент ослабления в точке 2 ; Т (0;Zo), 2

Т.(Zo,21) — прозрачность участков трассы

2 зондирования (0;Zo), (Zo, 2 ), где Zo, Zj — точки по трассе зондирования, Из (1) следует, что для решения уравнения необходимо ис5 пользовать какие-либо допущения о среде или априорную информацию ввиду многопараметрической зависимости. Решение уравнения (1) относительно/3{21) для экспериментально определяемой функции

$(Z)) = P(Z))Z Ац (Z))T (0;2o), 2 -1 -2 как показано имеет вид:

Z Z

P(Z ) = ф(2)(- 2 1 1/г(2)б2), Zp (2) где ф(2к) = Р(2к)2 Ag (ZK)T (0;Zo) — экспериментально определяемая функция в точке

2к, которая является калибровочной (опорной), В этой же точке 2к другим независимым способом определяют и значение коэффициента ослабления P (ZK). Удобнее всего точку ZK располагать вблизи лидара, так как на удалении от него (особенно для наклонных или вертикальных трасс) измерять P(2„) затруднительно или даже невозможно.

Следует отметить, что решение (2) для практической реализации требует предположения о неизменности лидарного отношения вдоль трассы зондирования, В случае зондирования слоисто-неоднородных сред, в которых наряду с изменением концентрации рассеивающих частиц по трассе зондирования происходит изменение их, микрофизических характеристик, а значит и лидарного отношения g(Z), использование алгоритмов, основанных на неизменнрсти

g(Z) становится некорректным. Очевидно и следствие из этого — рост ошибки определения Р (7) Данного недостатка можно избежать, если произвести корректировку

45 экспериментально измеряемой функции, на перепадлидарного отношения после границы раздела. Действительно, запишем функцию ф (Z) для точки Zl, лежащей до границы раздела, и точки 7I+1, лежащей после грани50 цы раздела; ф {21) = S(ZI)CI, ф (Zl+1) = S(ZI+1)CI+1 = Я(21+1)СЩЦ+1, (3) С = А Т {0;Zo)g(zl)

Cl+1= А Т (0;Zo)g(zl+I) 5Д где S(7I) = P(ZI)ZI; S(ZI+1) = P(Z;+1)ZI+1 — сигналы обратного рассеяния, восстановленные на квадраты расстояний ZI и 21+1.

1809408

Как видно из (3), при переходе от точки

i к точке i + 1 (переходе через границу раздела) происходит изменение функции tP (Z) за счет изменения лидарного отношения.

Величина этого изменения равна с!),(+1, так как Сн-1 = Ci.с (,(+1, где ц(,(+1 = gl/gl+1. А это означает, что при расчете /3 (Z) в точке Zl+1 необходимо производить корректировку в константе Сь используемой в точке Zi va величину с!ц+1, Физический смысл необходимости корректировки состоит в использовании одной константы для всей многослойной трассы зондирования, что означает сведение алгоритма к модели с постоянством лидарного отношения по всей трассе. Другими словами, для более корректного использования. методик, требующих предположения g(Z) = const, в случае многослойных сред необходима подстройка (корректировка) сигнала обратного рассеяния под используемое предположение, С учетом (3) выражение (2) можно переписать в виде;

p(Zj) = S(Zj)(р @- - 2(l)(i+ !(((:!ц+1), (4)

S (2к1

Zi где IK = j S(Z)d2 — сигнал, накопленный на

Ек участке от калибровочной точки 2» до границы раздела, Точка ZI(выбирается из-заудобства вблизи лидара и в ней определяется опорное значение коэффициента ослабления ((3(2к), Iij — сигнал, накопленный на участке от границы раздела Zi до исследуемой точки Zj.

Для определения же величины коэффициента коррекции ц!.и-1 необходимо измерение накопленных сигналов на двух участках произвольной ширины ЛЕ(до границы раздела Zi и двух участках шириной ЛZ после границы раздела, Запишем выражение для накопленных сигналов в виде функционала I, который. имеет вид (3):

l(Z, Z + хх 2) = — Aggl, 2 +EZ)T (О, 2) х

z+Az х (1 - ехр(-2 f  (Z)dz}).

При выполнении на участках Л2ь hZ условий;

Ххххх;Г хх - х. х2;х =

g(ZZI, Zl х- Z) =.g(ZI+ . Zl+ 2 02)= gl+I

T(Zl - 2 Л Ъ, 2> - Л Zi) = T(Zl - hZ,Z() =. Т(hZi)

T(Zl, Z + Л Z) = T(Zl + BZ, Zl + 2 Л Z) =

= T(h,Z) (5) для измеряемых накопленных сигналов до границы раздела !ь2 = I(Z — 2 Л Zi; Zi - Л21), Ii-1-- l(Zi - AZ(, Zi) и после границы раздела

Il = !(Ъ; Zl + hZ), II+1 = l(Zi + Л Z; Zl + 2 Л Z) можно записать систему уравнений

Il-2 = Ац)Т (О, Zi - 2 Л Zi)(1. - Т (Л2!)) ! (-1 = Ац(Т (О; Zl-2 AZ()T (ЛЪ)(1-Т (AZi)) ! (= Ацн-1Т (О; Zl - 2 Л Zi)T (Л Zi)(1- Т (hZ))

Ir+l = А0х-1Т (О, Zl -2 AZl)T (ЛZl)T х х (Ы) (1 - Т (Л Z)), (б) 10

Il — 1 Ii — Il +1

tII,1+1 = gl/gl+1 = (7)

ll — 2 Ii — 1

Параметр корректировки о(,н-1 зависит только от микрофизических характеристик рассеивающей среды и не зависит ни от концентрации рассеивающих частиц, ни от аппаратурных констант, Даже если одновременно со скачком лидарного отношения на границе раздела происходит и скачкообразное изменение коэффициента ослабления (концентрации), то и в этом случае величина (1(,1+1 показывает, какую часть в изменении сигнала обратного рассеяния составляет величина изменения лидарного отношения g(Z), т.е. изменение качественного состава, т,к, выражение (7) не содержит

Т (h,Zi) и T (Л2) и при выводе (7) никаких, поедположений о поведении перепада

Т (ЛЪ) и Т (ЛЕ) не делалось, Таким образом, как следует иэ (4), (7), измеряя накопленные сигналы от калибровочной точки до границы раздела. от границы раздела до точки определения коэффициента ослабления и на двух равных участках, прилегающих к границе раздела с двух сторон определяется значение коэффициента ослабления в исследуемой точке с использованием корректировки сигналов обратного рассеяния на перепад лидарного . отношения, Устранение при этом используемого в известных способах предположения g = const приводит к исключению погрешностей, обусловленн (х этим, Формула изобретения

Способ определения коэффициента ослабления рассеивающей среды путем облуВ где черта над g означает усреднения на участках накопления (Zi - 2 Я; Zi - Л 2(), (Zi - -Л Zl;

Zi), прилегающих к границе раздела слева, и (Ъ; Ъ+ ЛЕ), (Zi+ Л2, Zi+ 2 Л2) — прилегающих к границе раздела справа, Решение системы (6) относительно

20 ц(/ц(+! имеет вид:

1809408

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

P (ZI) = $(2Я вЂ” 2(!и + !Цыц+1), Z1 где IKt = f S(Z)dz — сигнал, накопленный от

2к калибровочной точки Z< до границы раздела Zi;

S{Z) = P(Z)Z — сигнал обратного расстояния от точки Z;

2)

ItI = J S(Z)dZ — сигнал, накопленный от

Zi границы раздела Z; до исследуемой точки ZI;

S(Z) = P(Z)Z;

Я(2к) = P(Z„)Z„

P(ZI) P(Z„) — сигналы обратного рассеяния соответственно из точек ZI и 7»,.

II — 1 (I| — II +1)

РЦ+1 =

zi — Azi zI !

А f P(Z)Z dZ; lt-1 = (P{Z)Z 07—

21 — 2 Azl zi — Л2( сигналы обратного рассеяния, накопленные на двух участках длиной AZi до границы раздела;

z) +Az zi +2 Az

1i = ) P(Z)Z dZ; li+> = f P(Z)Z dZ—

zl z) +Л сигналы обратного рассеянйя, накопленные на двух участках длиной AZ после границы раздела;

Z, Z< — расстояния до рассеивающего объема и калибровочной точки соответственно;

P (ZK) — коэффициент ослабления в калибровочной точке, 1809408

Фиг. 2

1809408

1809408

Ф р ч с Ю

Составитель М. Кугейко

Техред М,Моргентал Корректор Л. Филь

Редактор

Производственно-издательский комбинат "Патент", г. Ужгород, ул.Гагарина, 101

Заказ 1285 Тираж Подписное

ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР

113035, Москва, Ж-35, Раушская наб., 4/5