«Санкт-Петербургский государственный электротехнический университет
«ЛЭТИ» им. В.И.Ульянова (Ленина)»
(СПбГЭТУ «ЛЭТИ»)
Направление
11.03.04 – Электроника и наноэлектроника
Профиль
Квантовая и оптическая электроника
Факультет
ФЭЛ
Кафедра
Фотоники
К защите допустить
Зав. кафедрой
д.т.н. Тарасов С.А.
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА
БАКАЛАВРА
Тема: ИССЛЕДОВАНИЕ ПРИНЦИПОВ ПОСТРОЕНИЯ
ПОРТАТИВНЫХ ВЫСОКОТОЧНЫХ ЦИФРОВЫХ КОМПЛЕКСОВГАЗОДЕТЕКТОРОВ ДЛЯ ПРОМЫШЛЕННЫХ ПРИМЕНЕНИЙ
Студент
Карась Б.Ю.
подпись
Руководитель
к.т.н.
Редька Д.Н.
(Уч. степень, уч. звание)
Консультант
подпись
к.т.н.
Степанова О.С.
(Уч. степень, уч. звание)
Консультант
подпись
к.т.н., доц.
Буканин В.А.
(Уч. степень, уч. звание)
подпись
Санкт-Петербург
2020
1
ЗАДАНИЕ
НА ВЫПУСКНУЮ КВАЛИФИКАЦИОННУЮ РАБОТУ
Утверждаю
Зав. кафедрой Фотоники
____________ Тарасов С.А.
«20» апреля 2020 г.
Студент
Тема
Карась Б.Ю.
работы:
Исследование
Группа 6204
принципов
построения
портативных
высокоточных цифровых комплексов-газотекторов для промышленных
применений
Место выполнения ВКР: СПбГЭТУ «ЛЭТИ», каф. Фотоники.
Исходные данные (технические требования): обосновать использование
метода комбинационного рассеяния для построения комплексовгазотекторов, привести энергетический расчет системы, провести
моделирование монохроматора для комплекса.
Содержание ВКР: актуальность работы, обоснование метода КР,
энергетический расчет системы, моделирование монохроматора в
программном комплексе COMSOL Multiphysics.
Перечень отчетных материалов: пояснительная записка, иллюстративный
материал.
Дополнительные разделы: Безопасность жизнедеятельности.
Дата выдачи задания
Дата представления ВКР к защите
«20» апреля 2020 г.
«11» июня 2020 г.
Студент
Карась Б.Ю.
подпись
Руководитель
к.т.н.
Редька Д.Н.
подпись
2
КАЛЕНДАРНЫЙ ПЛАН ВЫПОЛНЕНИЯ
ВЫПУСКНОЙ КВАЛИФИКАЦИОННОЙ РАБОТЫ
Утверждаю
Зав. кафедрой Фотоники
____________ Тарасов С.А.
«20» апреля 2020 г.
Студент
Карась Б.Ю.
Группа 6204
Тема работы: Исследование принципов построения портативных
высокоточных цифровых комплексов-газотекторов для промышленных
применений
№
п/п
1
2
3
4
5
6
Наименование работ
Обзор литературы по теме работы
Подбор модулей для расчета
Энергетический расчет
Моделирование монохроматора
Оформление пояснительной записки
Оформление иллюстративного материала
Студент
Срок
выполнения
20.04 – 01.05
01.05 – 05.05
05.05 – 15.05
15.05 – 25.05
25.05 – 30.05
30.05 – 03.06
Карась Б.Ю.
подпись
Руководитель
к.т.н.
(Уч. степень, уч. звание)
Редька Д.Н.
подпись
3
РЕФЕРАТ
Пояснительная записка 87 стр., 45 рис., 24 табл., 52 ист.
РАМАНОВСКАЯ
СПЕКТРОСКОПИЯ,
ГАЗОВЫЙ
АНАЛИЗ,
МОНОХРОМАТОР ЧЕРНИ-ТЕРНЕРА
Объектом
исследования
является
портативный
рамановский
спектрометр для определения концентраций УВГ газов.
Цель работы: рассмотреть перспективы использования приборов на
основе метода комбинационного рассеяния при разведке нефти и газа,
предложить конфигурацию портативного рамановского газоанализатора,
произвести энергетический расчет и моделирование монохроматора.
В ходе работы были рассмотрены основные физические эффекты, на
которых строят приборы для газового анализа, проведены их недостатки в
сравнении с методом комбинационного рассеяния света. Определены
необходимые углеводородные газы для анализа при разведке нефти и газа.
Предложена схема построения портативного рамановского газового
анализатора. Произведен подбор лазерных модулей и фотоприемного
устройства и произведен энергетический расчет системы, для определения
минимальных
концентраций
метана.
Произведено
моделирование
монохроматоров по схеме Черни-Тернера в программном комплексе COMSOL Multiphysics.
4
ABSTRACT
The paper examined the prospects of using gas analysis to search for oil and
gas. This method is not presented as a replacement for traditional search methods,
but is considered as an addition to reduce the cost of oil and gas exploration. An
analysis is made of the existing physical principles on which gas analyzers are
built, and the main advantages of using the Raman scattering method are highlighted. A configuration of a portable Raman gas analyzer is proposed. The energy calculation according to the lidar equation was carried out, and the modeling of a
monochromator in the program COMSOL Multiphysics.
5
СОДЕРЖАНИЕ
Определения, обозначения и сокращения .......................................................... 8
Введение ............................................................................................................... 9
1. Актуальность работы ..................................................................................... 10
1.1 Перспективы использования газового анализ при разведки нефти и газа 10
1.2. Обзор и анализ существующих методов и приборов для газового анализа
............................................................................................................................. 17
1.2.1. Полупроводниковые сенсоры .................................................................. 18
1.2.2.Электрохимические сенсоры .................................................................... 19
1.2.3.Термокаталитеские .................................................................................... 20
1.2.4.Фотоионизационный сенсор ..................................................................... 22
1.2.5Хромотография ........................................................................................... 23
1.2.6. Физические основы метода ДП. .............................................................. 26
1.2.6.Физические основы метода комбинационного рассеяния. ..................... 28
Выводы по главе................................................................................................. 33
2. Энергетический расчет .................................................................................. 35
2.1. Концепция разрабатываемой системы ....................................................... 35
2.2 Решение лидарного уравнения .................................................................... 39
Выводы по главе................................................................................................. 54
3. Моделирование............................................................................................... 56
3.1. Оптический канал ........................................................................................ 56
3.2 Телескоп ........................................................................................................ 58
3.3 Монохроматор .............................................................................................. 60
Выводы по главе................................................................................................. 76
4. Безопасность жизнедеятельности .................................................................. 79
4.1 Общие положения по использованию лазеров в рамках данной работы .. 79
4.2 Анализ опасностей при работе с лазерами ................................................. 79
6
4.3 Предельно допустимые уровни лазерного излучения при воздействии на
глаза и кожу ........................................................................................................ 81
4.4 Меры и средства защиты от лазерного излучения ..................................... 82
Заключение ......................................................................................................... 84
Литература ............................................................Error! Bookmark not defined.
7
ОПРЕДЕЛЕНИЯ, ОБОЗНАЧЕНИЯ И СОКРАЩЕНИЯ
В настоящей пояснительной записке применяют следующие термины с
соответствующими определениями:
ДП – дифференциальное поглощение,
КР – комбинационное рассеяние,
ПДК – предельно допустимая концентрация,
ПДУ – предельно допустимый уровень,
УВГ – углеводородные газы.
8
ВВЕДЕНИЕ
Разведка
и
доразведка
нефтегазовых
месторождений
является
актуальной задачей отрасли, методы используемые сегодня крайне дороги
либо слабо информативны, один из возможных способов решение данной
проблемы это использование высокоточных комплексов для газового анализа
с целью определения УВГ аномалий.
Также растет количество и общая протяженность нефтегазопроводов, в
связи с чем, востребованы инновационные высокоточные методы и
технологий для осуществления экологического мониторинга в целях
обнаружения
утечек
добывающих
и
из
нефтегазопроводов
перерабатывающих
и
нефтезагрязнений
предприятий.
В
виду
от
низкой
концентрации в воздухе, для обнаружения в воздухе паров углеводородов
нефти и нефтепродуктов необходимо использование высокочувствительной
аппаратуры.
Существующие на сегодняшний день технические средства не могут в
полном объеме обеспечить качественное решение поставленных задач.
Поэтому требуется создание совершенно новых приборов основанных на
методе рамановской спектроскопии. Такой прибор способен проводить
детальный анализ исследуемых объектов и компонентов окружающей среды
по их спектральным характеристикам в условиях существенно ниже
предельного уровня пространственного разрешения существующих систем.
9
1. АКТУАЛЬНОСТЬ РАБОТЫ
1.1 Перспективы использования газового анализ при разведке
нефти и газа
Ежегодный
рост
количества
и
общей
протяженности
нефтегазопроводов, а также необходимость разведки новых и доразветки
существующих месторождений нефти и газа, повышают в востребованность
инновационных высокоточных методов и технологий для осуществления
газового анализа, экологического мониторинга для обнаружения утечек из
нефтегазопроводов
и
нефтезагрязнений
от
добывающих
и
перерабатывающих предприятий. В виду низкой концентрации в воздухе, для
обнаружения в воздухе паров углеводородов нефти и нефтепродуктов
необходимо использование высокочувствительной аппаратуры.
Доминирующее место среди методов разведки нефти и газа занимает
сейсморазведка, суть метода заключается в возбуждении упругих волн
искусственным способом (взрыв), анализ распространения волн по земной
коре
и
выявления
аномалий.
Порядка
70%
от
общего
объема
разведывательных работ осуществляется данным методом, но стоимость
разведки оценивается в несколько тысяч долларов за км2. Также используют
метод гравиаразведки, основан на сравнении
плотностей пород, области
содержащие нефть менее плотные, но данные получаемые таким методом
сложно интерпретировать. Электрический поиск основан на пониженной
электропроводности у нефтяных пластов, а магниторазведка основана на
неоднородности магнитной проницаемости слоев, содержащих нефтяные
карманы, оба метода нуждаются в высокочувствительной аппаратуре.
Главная проблема в поисках нефти и газа заключается в малой
площади перспективных областей по отношению к общему объему
исследуемого пространства. Регионы перспективные на поиск нефти и газа
зачастую
имеют
сложный
ландшафт,
10
покрытый
лесами
и
густой
растительностью, которые, в свою очередь,
затрудняют исследования
классическими методами. Поэтому в дополнение к традиционным методам
разведки следует вводить дополнительный метод, за счет которого при
низких финансовых затратах было бы возможно выявлять потенциально
перспективные участки. Так для этого можно использовать газовый анализ,
метод основан на эффекте миграции из недр углеводородов и формирования
облака над поверхностью земли, детектируя такое облако можно заключить о
нефтегазоносности подстилающей поверхности.
Процесс миграции углеводородов из залежей нефти и газа делят на два
процесса: миграцию посредством фильтрации и диффузию. Оба процесса
определяются геологическими условиями, первый обусловлен перепадом
давления вверх по залежи, а неоднородность концентрации вызывает
диффузию.
Диффузия
является
постоянным
процессом
миграции
углеводородов из залежи, но ее процесс снижается с ростом молекулярной
массы углеводородов.
Теоретически и экспериментально доказано: почти все горные породы
проницаемы для углеводородных газов; если залежь УВ-сырья перекрыта
плотными породами, то миграция осуществляется диффузно. Поток
вертикально
мигрирующих
стабилизируется
и
УВ
достигает
обеспечивает
появление
дневной
зон
с
поверхности,
повышенной
газоконцентрацией на участках проекции нефтегазовых скоплений, скрытых
на глубине [15].
11
Рисунок 1.1 – Физико-химическая модель залежи (по А. В. Петухову)
а - формы и виды геохимических аномалий: 1- апикальная аномалия
сорбированных
углеводородных газов;
2-
фильтрационные
аномалии
углеводородных газов и минеральных образований; 3- аномалии жидких и
газообразных УВ, фильтрующихся по зонам макро- и микротрещиноватости;
4- кольцевые аномалии газов биохимического происхождения
б
–
обобщенный
минералого-геохимический
разрез:
1-3
–
продуктивные отложения (1- залеж; 2- ВНК; 3- покрышка); 4-6 – комплекс
надпродуктивных отложений (4- нижняя геохимическая зона (комплекс
подбарерных отложений); 5- верхняя геохимическая зона (комплекс
надбарерных отложений); 6- первый регионально выдержанный водоупор
(литологический барьер); 7 – зона фильтрации УВ; 8 – барьерная зона.
Тяжелые газообразные углеводороды С2-С4 характерны только для
залежей. Они практически не образуются бактериями, и их генерация
12
органическим
веществом
в
зоне
геохимического
зондирования
незначительна. Поэтому при анализе полей концентрации этан, пропан и
бутан являются ведущими показателями. Углеводороды C5-C8 и более
тяжелые типичны не только для нефти, но их миграционные способности
ограничены в связи с большим молекулярным весом. Традиционно
считается, что они могут быть встречены только в районах с интенсивными
процессами вертикальной миграции – фильтрация по зонам тектонических
нарушений [16].
Рисунок 1.2 – Модель формирования углеводородных аномалий над залежей:
1 – горизонт: а - сингенетичной геохимической обстановки, б эпигенетичной геохимической обстановки; 2- шлейф миграции УВ; 3- УВ аномалия; 4 - вектор диффузии; 5 - зональность приповерхностного поля
концентрации УВ: а - кольцевого типа, б - линейного типа; 6 - залежь УВ; 7разрывные нарушения; 8 - кумулятивные кривые: а - ∑CH4, б - ∑ТУ, в –
усредненная [24]
Согласно работе [18] анализ закономерностей УВГ аномалий позволил
сделать вывод: чем сложнее строение органического вещества, тем больше
13
вероятность его прямого биосинтеза (углеводороды, жиры, белки), и тем
меньше вероятность его эмиграции из глубоких недр на поверхность.
Анализ эффективности применения геохимических методов при
поисках месторождений нефти и газа за рубежом показывает, что
отрицательная оценка перспектив нефтегазоностности локальных объектов
по результатам геохимических исследований подтверждается на 100%, а
положительная на 60-80 % [17]
Рисунок 1.3 – Результаты положительного прогноза на Алимовском
поднятии Мензелино-Актанышской зоны: а – карта аномалий естественного
электрического поля; б – карта аномалий магнитного поля; в – карта
аномалий геохимического поля [25]
Сам по себе состав природного газа индивидуален для каждого
месторождения, компонентные составы газа некоторых месторождений
приведены в таблицах 1.1, 1. 2, 1.3, 1.4.
Таблица 1.1 – Состав природного газа некоторых газовых месторождений [14]
Месторождение
Состав газа, % мас.
CH4
C2H6
C3H8
C4H10
C5H12
CO2
N2
Ленинградское
88,2
4,5
2,3
0,7
0,3
0,6
3,4
Уренгойское
96,00
0,09
0,01
-
0,01
0,49
3,4
14
Продолжение таблицы 1.1
Степновское
92,2
2,6
1,0
0,5
0,3
1,2
1,9
Медвежье
99,20
0,08
0,01
0,07
0,02
0,06
0,57
Ямбургское
95,20
0,04
0,01
-
0,01
0,30
4,00
Вуктылское
75,7
9,1
3,1
0,7
7,5
0,2
3,7
Губкинское
98,6
0,12
0,01
0,01
0,01
0,1
1,2
Заполярное
98,5
0,20
0,05
0,01
-
0,5
0,70
Марковское
81,5
7,9
2,9
0,9
0,4
0,2
5,4
Таблица 1.2 – Состав газа некоторых газоконденсаторных месторождений
Месторождение
Состав газа, % мас.
CH4
C2H6
C3H8
C4H10 C5H12
CO2
N2
H2S
Астраханское
54,15
5,54
1,68
0,93
1,57
21,55
1,98
12,60
Оренбургское
81,70
4,50
1,80
1,00
3,55
2,35
1,1
4,00
Карачаганское
75,31
5,45
2,62
1,37
5,98
4,79
0,79
3,69
Вуктыльское
75,00
9,00
3,90
1,80
5,20
0,10
5,00
-
Шатылыкское
95,70
1,70
0,23
0,04
0,02
1,24
1,4
-
Состав газов нефтяных месторождений и нефтяных попутных газов
заметно отличаются по составу от месторождений природного газа:
Таблица 1.3 – Состав газа некоторых нефтяных месторождений
Месторождение
Состав газа, % мас.
CH4
C2H6
C3H8
C4H10 C5H12
CO2
N2
H2S
Ромашкинское
39,00 20,00 18,50
6,20
4,70
0,10
11,50
-
Туймазинское
41,00 21,00 17,40
6,80
4,60
0,10
7,10
2,00
Шкаповское
47,00 14,10 27,20
9,50
5,20
-
-
-
Жирновское
82,00
3,50
1,00
5,00
-
1,50
Мухановское
30,10 20,20 23,60 10,60
4,80
1,50
6,80
2,40
6,00
3,00
15
Таблица 1.4 – Состав ПНГ различных месторождений Западной Сибири [13]
Состав газа, % мас.
Месторожден
ие
CO2
N2
CH4
C2H6
C3H8 C4H10 C5H12
Самотлорское
0,59
1,48
60,64
4,13
13,05
4,04
2,52
2,65
Варьеганское
0,69
1,51
59,33
8,31
13,51
4,05
2,2
1,8
Аганское
0,5
1,53
46,94
6,89
17,37
4,47
3,36
3,88
Советское
1,02
1,53
51,89
5,29
15,57
5,02
2,99
3,26
Губкинское
0,2
1,70
97,70
0,27
0,01
0,01
0,001
0,11
Уренгойское
0,17
1,10
89,60
4,36
2,01
0,90
1,70
0,16
Лугинецкое
0,81
6,84
85,66
4,03
1,89
0,61
0,14
0,02
Талинское
1,48
0,54
48,71 10,05 16,27
7,43
2,77
12,75
Прочие
Проводя анализ информации указанных выше исследований, а также
компонентного состава месторождений нефти и газа можно определить
список индикаторных углеводородов необходимых для проведения газового
анализа с целью поиска месторождений нефти и газа:
Таблица 1.5 – Список газов индикаторов нефтегазоносности
Вещество
Формула
Метан
CH4
Этан
C2H6
Пропан
C3H8
Бутан
C4H10
Изо-бутан
iC4H10
Пентан
С5Н12
Изо-пентан
iC5H12
Главной проблемой при таком методе разведки является низкая
концентрация паров нефти, что требует применения инновационной
16
аппаратуры высокой точности. На сегодняшний день разработано большое
количество
различных
систем
для
осуществления
эколого-
метеорологического мониторинга, но, как показывает практика, большинство
прошедших
лабораторные
испытания
методов
оказываются
малоэффективными при эксплуатации в рабочих условиях. Помимо геолого
разведывательных работ одной из важнейших проблем нефтяной отрасли
является экологический мониторинг в целях обнаружения утечек из
нефтегазопроводов
и
нефтезагрязнений
от
добывающих
и
перерабатывающих предприятий. Для комплексного решения данных задач
перспективным является использование лазерных систем мониторинга, их
достоинствами
являются:
дистанционность,
высокая
высокая
скорость,
помехоустойчивость,
возможность
точность,
одновременного
определения широкого набора химических элементов и соединений.
1.2. Обзор и анализ существующих методов и приборов для
газового анализа
Для осуществления газового анализа существует множество методов,
на сегодняшний день существует тенденция к миниатюризации систем для
газового анализа.
Существует множество физических принципов, на базе которых
созданы приборы для газового анализа, однако большинство из них обладает
существенным недостатком: отсутствием селективности, в виду чего
приходится комбинировать в одном приборе множество сенсорных систем.
Главным методом для анализа химического состава газа, попутных
газов,
газового
конденсата
и газов
нефтепереработки
является
хроматография. Данный считается самым достоверным, а также позволяет
определять широкий набор веществ: тяжелые углеводороды от метана до
пентанов,
кислород, азот, водород, гелий,
другие.
17
окислы углерода, и многие
В классе личных портативных устройств используются персональные
газовые анализаторы, которые рассчитывают суммарную концентрацию
горючих веществ в атмосфере, и служат для предупреждения угрозы жизни и
анализа помещений на предмет утечек.
Для разведки месторождений и мониторинга объектов нефтегазовой
отрасли
широкое
распространения
получили
оптические
системы,
основанные на методе дифференциального поглощения, которые отличаются
высокой чувствительностью и точностью измерений.
Рассмотрим подробнее основные типы сенсоров и физические
принципы, на которых построены современные приборы для газового
анализа.
1.2.1. Полупроводниковые сенсоры
Полупроводниковые сенсоры работают за счет реакции хемосорбциии
молекул детектируемого компонента на нагретом чувствительном элементе
(рисунок 1.4). В современных сенсорах на основе полевых транзисторов
используется принцип действия МДП-транзистора с индуцированным
каналом [5]. Конструкция датчика представлена на рисунке 1.5: нагретый
элементом 8 полупроводниковый датчик 1 состоит из чувствительной
мембраны 6, нанесенной на диэлектрик, который отделяет ее от контактов 7,
мембрана управляет током стока, протекающим между электродами истока и
стока [5]. Проволоки 2 соединяют чувствительный элемент с электродами 4,
сенсор защищен корпусом 5, в котором имеется сетка, через которую газ
попадет внутрь.
Сигнал,
снимаемый
с
полупроводникового
сенсора,
имеет
логарифмическую зависимость от концентрации детектируемого газа, что
пагубно сказывается на точности измерений и уменьшает пределы
измерений.
Сильное
воздействие
18
на
чувствительность
оказывают
«отравляющие вещества» такие как щелочные или кислотные соединения,
силиконы, тетраэтилсвинец, серные соединения, галогенные соединения [6].
Рисунок 1.4 – Конструкция полупроводникового датчика
1.2.2. Электрохимические сенсоры
Принцип
работы
электрохимического
сенсора
основан
на
окислительно-восстановительной реакции протекающей на электроде, в
результате чего изменяется аналоговый сигнал в цепи (рисунок 1.5). Газ
попадает внутрь датчика через капиллярное отверстие 1, проходит через
газопроницаемую мембрану 2 к рабочему электроду 3, на котором
происходит реакция окисления либо восстановления. На вспомогательном
электроде 5 происходит зеркальная рабочему электроду реакция, подводя к
электродам
напряжение,
в
цепи
появляется
электрический
ток,
пропорциональный концентрации детектируемого газа. Для поддержания
постоянного напряжения на рабочем электроде в схему добавляют электрод
сравнения 4, все три электрода помещены в электролит 6.
Зависимость
тока
от
концентрации
практически
линейна,
электрохимические сенсоры высоко чувствительны и селективны, но не
19
пригодны для определения состава многокомпонентных газовых смесях,
поэтому часто используются в мультисенсорных приборах.
Рисунок 1.5 – Структура электрохимического газового датчика
1.2.3.Термокаталитеские
Достаточно большое количество газоанализаторов, представленных на
рынке, являются термокаталитическими. Сигнал снимается со специального
сенсора, в настоящее время серийно выпускают лишь объемные (3D), т. е.
сенсоры пеллисторного типа [4].
Пеллисторный сенсор схематически представлен на рисунке 1.6.
Платиновая спираль 1
покрывается
сверху пористым керамическим
материалом 2, на который наносят катализатор, спираль включается в общую
схему 3.
На катализаторе адсорбируется кислород, и при прохождении горючего
газа
через
поры
керамики
происходит
сопровождающаяся выделением тепла:
𝐶𝐻4 + 2𝑂2 = 2𝐻2 𝑂 + 𝐶𝑂2 + 𝑄
20
реакция
окисления,
Рисунок 1.6 – Структура пеллисторного сенсора
Спираль, разогретая электрическим током, дополнительно нагревается
в результате реакции окисления, изменяется ее сопротивление, и сигнал в
цепи соответственно.
Величина
изменения
сигнала
пропорциональна
концентрации детектируемого газа, для метана диапазон измеряемых
объемных концентраций обычно составляет 0–4,4 % [4]. Для компенсации
условий среды используют второй пеллисторный сенсор, подключенный по
мостовой схеме, с той лишь разницей, что на него не наносят катализатор.
Термокаталитические сенсоры привлекательны своей низкой ценой и
простотой в сравнении с оптическими датчиками, имеют линейную
градуировочную
характеристику,
способны
детектировать
большое
количество углеводородов и водород. Но имеют ряд существенных
недостатков: катализатор имеет свойство устаревать, из-за чего со временем
падает чувствительность сенсора, также катализатор может быть отравлен
соединениями свинца, фосфора и серы, вследствие чего некорректно
работать. Для реакции окисления необходим кислород, его содержание в
воздухе должно быть не менее 12%, в противном случае пеллисторный
сенсор попросту не будет работать. Данные газоанализаторы как правило
21
используют для регулярного мониторинга объектов, где возможны утечки
горючих веществ, с целью предупреждения об опасности.
1.2.4.Фотоионизационный сенсор
Принцип работы приборов основан на измерении электрического тока,
обусловленного
ионизацией
детектируемого
вещества
посредством
воздействия источника вакуумного ультрафиолетового излучения, общая
схема представлена на рисунке 1.7.
Рисунок 1.7 – Структура фотоионизационного сенсора
Ультрафиолетовое излучение, излучаемое лампой 1, воздействует на
газовую смесь, попадающую через фильтр-мембрану 3 в межэлектродное
пространство 2. Молекулы энергия ионизации которых, меньше энергии
излучения ионизируются, и под действием электрического поля между
электродами
2
обуславливают
ток,
пропорциональный
концентрации
молекул детектируемого вещества.
Основные компоненты воздуха такие как: O2, N2, Ar, CO2 и др. имеют
значительно большую энергию ионизации чем большинство углеводородов.
Аналоговый сигнал, снимаемый с электродов, усиливается и преобразуется в
значение массовой концентрации компонента.
22
Таблица 1.6 – Энергия ионизации некоторых веществ
Вещество
Энергия ионизации, эВ
Метан
12,98
Бутан
10,63
Гексан
10,18
Пропан
11,07
Гептан
10,07
Изобутелен
9,43
Пентан
10,33
При нахождении в воздухе паров нескольких детектируемых веществ,
прибор способен определять лишь суммарную концентрацию, т.е. является
не селективным и для анализа состава газа может быть использован в
комплексе с другими датчиками.
1.2.5 Хромотография
Газовая хроматография – основной метод анализа компонентного
состава природного газа, попутных газов, газового конденсата и газов
нефтепереработки. С его помощью можно разделять и определять
количественно все компоненты природного газа – постоянные газы (водород,
гелий,
кислород,
азот,
окись
и
двуокись
углерода),
газообразные
углеводороды от метана до пентанов, тяжелые углеводороды (гексаны и
выше), соединения серы (сероводород, меркаптаны, сульфиды). Основные
проблемы качественного и количественного анализов природного газа –
широкие
пределы
концентраций
компонентов,
различные
сочетания
отношений компонентов газа, невозможность определения всех компонентов
за один анализ (постоянные и углеводородные газы, сернистые соединения,
23
вода, привнесенные в результате обработки газа примеси, например, метанол,
гликоли, амины и др.) [8].
Хроматография – это метод разделения и определения веществ,
основанный на распределении компонентов между двумя фазами –
подвижной и неподвижной [3]. В основе лежит разность сорбируемости
веществ смеси, в процессе ее движения по сорбенту. Адсорбция может быть
обусловлена неспецифическими (ориентационными,
дисперсионными)
и
специфическими
индукционными и
взаимодействиями
(комплексообразованием, либо образованием водородной связи) и зависит от
природы адсорбента и сорбата [3]. Подвижная фаза - газ или жидкость,
неподвижная фаза – жидкость или твёрдое тело, при разделении компоненты
смеси движутся с разными скоростями.
Существует несколько видов сорбционных процессов: абсорбция,
адсорбция и хемосорбцию (при получении химических соединений на
адсорбенте). Все три процесса могут быть реализованы одновременно, но в
основном используют один из трех, в зависимости от основного процесса,
получается
и
название:
адсорбция
-
метод
разделения
называют
адсорбционной хроматографией; абсорбция - то метод разделения называют
распределительной хроматографией.
Рассмотрим схему разделения двухкомпонентной смеси (рисунок 1.8),
1: подача исследуемой смеси A+B в подвижную фазу С, движущуюся по
хроматографической колонне, там смесь контактирует с неподвижной фазой
D, в результате процессов адсорбции скорость компонентов смеси начинает
розниться, получаем зоны A+C и B+C, чем более компоненты продвигаются
по колонне, тем больше зона чистой подвижной фазы меду ними. Так
получается
последовательность
компонентов,
детектором на стадии 3.
24
которая
регистрируется
Рисунок 1.8 – Схема хроматографического анализа двухкомпонентной смеси
Измерение площади пика проводят с использованием интегрирующих
измерительных приборов, площадь пика пропорциональна массе вещества,
градуировочная характеристика зависит от того, какой параметр принимают
за выходную характеристику: высоту пика h или его площадь S:
ℎ = 𝑘ℎ 𝑞 или 𝑠 = 𝑘𝑠 𝑞 ,
(1.1)
где 𝑘ℎ и 𝑘𝑠 чувствительность хроматографа.
Каждый компонент проходит колонну за разные промежутки времени,
поэтому временной параметр является качественной характеристикой в
распознавании. Время от момента ввода дозы в хроматограф до момента
получения максимума хроматографического пика называется временем
удерживания τR. В колонке оно складывается из двух составляющих: времени
τ0 нахождения вещества в подвижной фазе и времени 𝜏𝑅′ , в течение которого
вещество находится в сорбенте [9]:
𝜏𝑅 = 𝜏0 + 𝜏𝑅′ .
(1.2)
Время удерживания зависит от скорости подвижной фазы, поэтому на
практике используют понятие удерживаемого объёма:
𝑉𝑅 = (𝜏𝑅 − 𝜏0 )𝑄 ,
где Q – объёмный расход подвижной фазы.
25
(1.3)
1.2.6. Физические основы метода ДП
В основе методе измерений лежит закон Ламберта-Бера, согласно
поглощение электромагнитного излучения веществом или газом можно
записать в виде:
𝐼 (𝜆) = 𝐼0(𝜆)𝑒𝑥𝑝(−𝑘𝜆 𝑅 ) .
(1.4)
Здесь I0(λ) – интенсивность зондирующего источника, 𝑘𝜆 = 𝜎(𝜆)𝑁 –
коэффициент поглощения: σ(λ) - сечение поглощения на длине волны λ
молекул концентрации N. Сечение поглощения является индивидуальной
характеристикой
исследуемого
вещества.
Тогда
можно
рассчитать
концентрацию молекул исследуемого вещества:
𝐼 (𝜆)
𝑁 = ln ( 0
𝐼(𝜆)
)
1
𝜎(𝜆)𝑅
.
(1.5)
Метод ДПР заключается в сравнении двух обратно рассеянных
лазерных импульсов, частоты генерации которых лежат близко друг к другу.
На
практике
для
достижения
высокой
чувствительности
и
пространственного разрешения сравнивают обратно рассеянные сигналы.
Длины волн излучения подбираются таким образом, чтобы одна из них
соответствовала максимуму линии поглощения исследуемой молекулы,
вторая на крыле линии поглощения исследуемой молекулы.
Рисунок 1.9 – Схема лидарного метода дифференциального поглощения
26
Лидарное уравнение для однократного упругого обратного рассеяния
можно записать в виде:
𝑃 (𝜆, 𝑅 ) = 𝑃0(𝜆)𝐾1𝐴0 𝑇 2(𝜆, 𝑅) 𝜌⁄𝑅 2,
(1.6)
где P(λ, r) – мощность сигнала обратного рассеяния на фотоприёмнике на
длине волны лазера λ, приходящая с расстояния R; P0(λ) – мощность лазера;
К1 – постоянная лидара; А0 – площадь зеркала приёмного телескопа; Т(λ) –
пропускание атмосферы на длине волны λ; ρ – суммарный коэффициент
упругого рассеяния Ми и молекулярного рассеяния Релея.
Вся информация о концентрации молекул детектируемого газа
содержится в сомножителе Т(λ), который в общем случае может быть
представлен в виде
𝑅
𝑇(𝜆, 𝑅 ) = 𝑒𝑥𝑝 [− ∫0 𝑘 (𝜆, 𝑅 )𝑑𝑅 ] ,
(1.7)
где коэффициент ослабления k(λ,r) определяется соотношением
𝑘 (𝜆, 𝑅 ) = 𝑘𝐴 (𝜆, 𝑅 ) + 𝑁(𝑅 ) ∗ 𝜎0(𝜆),
(1.8)
в котором первое слагаемое является коэффициентом ослабления атмосферы
на длине волны лазерного излучения за вычетом молекул детектируемого
газа, а второе – произведение концентрации молекул детектируемого газа на
сечение поглощения этих молекул.
Определение концентрации детектируемых молекул происходит при
делении уравнений для двух зондирующих длин волн. Которые подбираются,
так чтобы коэффициенты обратного рассеяния и ослабления атмосферы не
зависели от длины волны. Тогда можно будет пренебречь отдельными
слагаемыми
в
выражении
для
определения
концентрации
молекул
детектируемого газа. При этом выражение для концентрации упрощается к
виду:
1
𝑃0 (𝑅+∆𝑅)𝑃1 (𝑅)
2
𝑃1 (𝑅+∆𝑅)𝑃0 (𝑅) ∆𝜎∆𝑅
𝑁 = − 𝑙𝑛 (
)
1
,
(1.9)
здесь Δσ = σ(λ0) - σ(λ1) – сечение дифференциального рассеяния; σ(λ0) и σ(λ1) –
сечения
поглощения
на
соответствующих
27
длинах
волн
λ0
и
λ1,
соответственно; P0(R), P0(R + ΔR) и P1(R), P1(R + ΔR) – мощности
отражённых сигналов для соответствующих длинах волн на расстояниях R и
(R + ΔR); ∆𝑅 =
𝑐∆𝑡
2
– пространственное разрешение лидарной системы; с –
скорость света; Δt – длительность лазерного импульса.
Для
оценки
чувствительности
лидарного
метода
ДПР
удобно
пользоваться формулой:
𝑁мин (𝑅 ) =
1
2
∆𝜎∆𝑅 𝑆𝑁𝑅(𝑅)
,
(1.10)
где SNR(R) – отношение сигнал/шум на расстоянии R от лидара.
На чувствительность метода дифференциального поглощения в
значительной мере могут влиять: ширина линии лазера, доплеровское
уширение, температурные вариации, поглощение «мешающими» газами,
сдвиг линии поглощения, нестабильность длины волны лазера.
1.2.6.Физические основы метода комбинационного рассеяния.
При регистрации рассеянного на частицах излучения наблюдаются
серии боковых частот
сдвинутые
на
величины равные
частотам
колебательно – вращательных переходов облучаемых молекул. При этом не
происходит ни поглощения, ни генерации фотона, оболочка молекулы
поляризуется под действием излучения.
Боковые частоты с меньшей частотой относительно зондирующего
излучения называют стоксовые (рисунок 1.10 а), а с большей антистоксовые
(рисунок 1.10 б), стоксовая частота i-ой компоненты определяются по
формуле:
𝑣𝑖 = 𝑣0 −𝑣кв .
28
(1.11)
Рисунок 1.10 – Схемы энергетических переходов
Метод комбинационного рассеяния позволяет определять отношение
концентрации любой составляющей к концентрации некоторой опорной
компоненты при условии, что известны соответствующие сечения.
На рисунке 1.11 представлена схема, отображающая принцип работы
рамановского спектрометра. Зондирующий пучок, генерируемый лазером,
направляется на объект исследования, обратно рассеянное излучение
улавливается
телескопической
системой.
Спектр
рассеянного
назад
излучения состоит из компонентов рэлеевского и Ми-рассеяния с частотой v0,
совпадающей с
частотой лазерного
излучения,
а
также
из
линий
комбинационного рассеяния с частотами v1, v2, ..., vn и различных
флуоресцентных компонентов с частотами vI, v II, ..., v N, расположенных
полосами или непрерывно с учетом разнообразия газов, составляющих
атмосферу. Спектральные компоненты выделяются и регистрируются
одновременно с помощью анализатора спектра вместе с оптическими
фильтрами и блоком чувствительных фотодетекторов.
.
29
Рисунок 1.11 – Принципиальная схема рамановского спектрометра
В схеме, использующей комбинационное рассеяние, не требуется
селективного возбуждения отдельных частот, ибо там рассеянный в
направлении назад сигнал содержит спектры комбинационного рассеяния
всех активных к этому процессу составляющих воздуха с интенсивностями
линий, пропорциональными их концентрациям [1]
𝑃 (𝜆, 𝑅 ) = 𝑃𝐿 𝐾𝑙 ∆𝑅𝐴0𝑇(𝜆𝐿 , 𝑅)𝑇(𝜆, 𝑅) (
𝑑𝜎
) 𝑁𝑎 ⁄𝑅 2 ,
𝑑𝛺
где PL -мощность лазерного излучения,
λL – длина волны зондирующего излучения,
K l – постоянная лидара,
∆𝑅 =
𝑐𝜏𝐿
2
– пространственное разрешение,
c – скорость света,
с
𝜏𝐿 = (𝑡имп + 𝑡вз + 𝑡эс )–время одного измерения,
2
tимп – длительность лазерного импульса,
30
(1.12)
tвз – время оптического взаимодействия данной компоненты с
излучением,
т.к.
комбинационное
рассеяние
происходит
практически
мгновенно, им можно пренебречь,
tэс – постоянная времени регистрирующего тракта,
А0 – площадь приемной аппаратуры,
𝑇(𝜆𝐿 , 𝑅 ), 𝑇(𝜆, 𝑅) – пропускание атмосферы на соответствующих
длинах волн
𝑅
𝑇(𝜆𝐿 , 𝑅 ) = 𝑒𝑥𝑝 [− ∫ 𝑘 (𝜆𝐿 , 𝑅 )𝑑𝑅 ],
(1.13)
0
𝑅
𝑇(𝜆, 𝑅 ) = 𝑒𝑥𝑝 [− ∫ 𝑘 (𝜆, 𝑅 )𝑑𝑅 ],
(1.14)
0
где 𝑘 (𝜆, 𝑟) = 𝑘𝑚𝑜𝑙 + 𝑘aer + 𝑘abs, объёмный коэффициент ослабления за счёт
рассеяния
на
аэрозолях,
молекулах
и
поглощения
атомарными
и
молекулярными составляющими атмосферы.
(
𝑑𝜎
𝑑𝛺
)-дифференциальное сечение КРС исследуемой молекулой на
длине волны лазерного излучения
Na – концентрация молекул
По значениям спектральной яркости солнечного излучения для длин
волн полос КРС исследуемыми молекулами Sb(λ), можно оценить степень
фоновых шумов:
𝑃𝑏 (𝜆, 𝑅 ) = 𝑆𝑏 (𝜆)𝑇(𝜆, 𝑅 )𝐾𝑏 𝜉𝑏 (𝜆)𝐴0𝛺(𝑅 )∆𝜆,
(1.15)
где 𝛺(𝑅 ) = 𝐴0⁄𝑅 2 - телесный угол поля зрения приемного телескопа, ∆λспектральная ширина приемного тракта.
Благодаря достаточному сдвигу частоты исходной лазерной линии
относительно частоты линии КР газа, выделение сигналов в лидаре КР легко
проводится традиционными фильтрами. Существенным недостатком этого
метода является малое сечение КР, т.е. низкая чувствительность и малое
отношение сигнал/шум при малых концентрациях газа [10].
31
Метод рамановской спектроскопии
имеет ряд преимуществ над
методом дифференциального поглощения:
1. В методе ДП для достижения высокого числа детектируемых газов
приходится
использовать
сложные
системы
перестраиваемых
лазеров с узкой линией генерации, что повышает стоимость и время
измерений. В свою очередь метод КР позволяет использовать
многоканальную систему регистрации интенсивности спектральных
линий, поэтому возможно определение широкого пула компонент
при единичном измерении с высокой скоростью.
2. На чувствительность оптических систем, применяемых в методе
ДПР, влияет не только рассеянное солнечное излучение днём, но и
тепловой фон Земли, особенно в летних условиях.
3. Многоканальная система регистрации позволяет создавать приборы
нечувствительные
к
флуктуациям
некоторых
параметров
окружающей среды (температура и давление).
Помимо объективных плюсов метод КР имеет ряд трудностей
связанных со спецификой анализа УВГ:
1. Метан, концентрация которого очень высока, имеет КР спектре две
интенсивные
полосы,
деформационным
соответствующие
колебаниям
C-H
с
валентным
богатой
и
вращательной
структурой.
2. УВГ молекулы построены из одного и того же набора атомов с
одинаковыми связями, поэтому наиболее интенсивные колебания
C-H связей близки друг к другу
3. Для
количественного
анализа
легких
углеводородов
(этана,
пропана, бутана) можно использовать менее интенсивные линии
валентных колебаний C-C связей, но для более тяжелых молекул
они начинают сливаться
32
Выводы по главе
Проведенный анализ показывает перспективность метода КР при
проведении нефтеразведовательных работ. Также можно составить список
газов-индикаторов необходимых для проведения газового анализа в
вышеуказанных целях:
Таблица 1.7 – Список необходимых газов-индикаторов
Вещество
Формула
ПДК, мг/м3
ПДК, ppm
Метан
CH4
300
450
Этан
C2H6
50
40
Пропан
C3H8
300
164
Бутан
C4H10
200
83
Изо-бутан
iC4H10
300
-
Пентан
С5Н12
100
33
Бутилен
С4Н8
10
-
Метилмеркаптан
CH4SH
Этилмеркаптан
C2H5SH
Азот
N2
-
-
Кислород
О2
-
-
Геохимический способ разведки является более дешевым и мобильным
методом, но он имеет ряд недостатков: поскольку аномалии могут возникать
из-за тектонических нарушений, высока сложность определения источника
аномалии относительно нее самой, также трактовка аномалии может быть
ошибочной
т.к.
довольно
трудно
определить
пороговые
значения,
характеризующие нефтегазонесущие области. Ввиду этого метод газового
анализа интересен в качестве выявления перспективных участков, за счет
чего можно существенно сократить расходы на разведку.
33
2. ЭНЕРГЕТИЧЕСКИЙ РАСЧЕТ
2.1. Концепция разрабатываемой системы
В случае рамановской спектроскопии информативными параметрами
полезного сигнала являются: интенсивность, которая несет информацию о
числе частиц, и волновой сдвиг, который индивидуален для каждого
вещества и позволяет определять детектируемый газ. Из этого вытекают
основные требуемые параметры системы:
1. Высокое соотношение сигнал/шум
2. Высокое разрешение монохроматора
Общая схема лазерных измерений показана на рисунке 2.1. Источник
излучения направляет лазерный луч в среду распространения, где он
ослабевает за счет поглощения и рассеяния, попадет на исследуемый объект
и отражается в обратном направлении, попадая на приемник излучения, где
сигнал преобразуется в аналоговый, усиливается и обрабатывается.
Рисунок 2.1 – Обобщенная схема лазерных иземерений [28]
35
Прием оптического сигнала всегда происходит на фоне оптических
помех: обратного рассеяния, фоновых, флуктуационных, в условиях, когда
полезный сигнал значительно ослаблен и искажен из-за влияния среды
распространения.
Для
эффективной
борьбы
с
фоновым
излучением
можно
воспользоваться двумя источниками излучения в двух разных спектральных
диапазонах, тогда суммарное соотношение сигнал/шум можно оценить
выражением:
𝑆𝑁𝑅 = √𝑃12 + 𝑃22 ,
(2.1)
где P1 и P2 – величины фона соответственно. Использование двух источников
и, как следствие, получение двух спектрограмм, за счет цифровой обработки
и сравнения позволит существенно увеличить точность системы.
Правильный выбор длины волны источников излучения также
существенно влияет на итоговую предельную способность обнаружения
газоанализатора. На рисунке 2.2 представлен энергетический спектр
солнечного излучения до его входа в атмосферу и на уровне моря при
отсутствии облачности. На этом же рисунке дан энергетический спектр
излучения абсолютно черного тела при температуре поверхности Солнца
(~5900 К). Атмосфера оказывает сильное влияние на общий спектр
солнечного
излучения,
появляются
резкие
провалы,
обусловленные
селективным поглощением излучения атмосферными газами. Это явление
создает солнечно слепые области спектра, которые весьма перспективны для
лазерного зондирования, без солнечной засветки можно получить очень
высокое соотношение полезный сигнал/шум.
36
Рисунок 2.2 - Энергетический спектр Солнца и полосы поглощения
атмосферы Земли [27]
Интенсивность стоксовой компоненты в рамановском спектре можно
описать следующей формулой:
𝐼𝑐 = 𝐾𝐼0𝜎𝑁 [1 − exp(−
𝐸𝑖
)],
𝑘𝑇
(2.2)
где К-константа; I0-интенсивность зондирующего излучения; σ – сечение КР;
N- число рассеивающих частиц в 1 см 3; Ei –энергия колебательно
вращательных уровней; k-постоянная Больцмана; T- газовая температура.
Перспективным является ультрафиолетовый диапазон 220 – 290 нм,
помех в данном диапазоне почти нет, интенсивность рамановских сдвигов
получается максимальной, зондирование можно осуществлять, используя 4ую гармонику Nd:YAG лазера, на длине волны 266 нм.
В качестве второго диапазона прекрасно подойдет длина волны 785 нм,
во-первых, самые интенсивные сдвиги углеводородов имеют значения
волнового числа колебательного перехода более 2800 см-1, тогда спектры
37
будут располагаться в диапазоне от 1 мкм до 1.1 мкм, где не наблюдается ни
одной сильной полосы поглощения атмосферными газами. Во- вторых, при
зондировании на источниками 266 нм и 785 нм, общий рабочий
спектральный диапазон будет от 270 нм до 1100 нм, данная область является
чувствительной для большинства кремниевых матриц, поэтому возможна
реализация при использовании одного приемника для двух сигналов
(рисунок 2.3). Лазер 266 нм необходим для получения сигнала в режиме
реального времени, а лазер 785 должен работать в режиме накопления
импульса.
Рисунок 2.3 – Функциональная схема рамановского газоанализатора
Также в системе необходимо реализовать канал для определения
концентрации
в
атмосфере
реперного
газа.
Метод
рамановской
спектроскопии позволяет определять концентрации любых других газовых
компонент по концентрации опорного газа. Таким газом может являться азот,
так его концентрация в атмосфере является постоянной.
38
Отношение мощности
концентрации
сигнала от определяемого компонента к его
пропорционально
отношению
мощности
сигнала
от
атмосферного азота к его концентрации соответственно. Тогда, зная
мощность сигнала от азота Pt и его концентрацию Nt, а также мощность
сигнала от исследуемого вещества Ps, можно определить его концентрацию
согласно формуле:
𝑁𝑠 =
𝑃𝑠 𝑁𝑡
.
𝑃𝑡
(2.3)
2.2 Решение лидарного уравнения
В случае, если поле зрения приемной системы лидара соответствует
поперечному сечению зондирующего луча на объекте наблюдения или
перекрывает это сечение, то используемая в работе оценка мощности
принимаемого лидарного эхо-сигнала согласно формуле для атмосферы
может быть представлена в виде:
𝑃(𝜆, 𝑅 ) = 𝑃𝐿 𝐾𝑙 ∆𝑅𝐴0 𝑇(𝜆𝐿 , 𝑅)𝑇(𝜆, 𝑅) (
𝑑𝜎
) 𝑁𝑎 ⁄𝑅 2,
𝑑𝛺
(2.4)
где PL -мощность лазерного излучения, λL – длина волны зондирующего
излучения, Kl – постоянная лидара, ∆𝑅– пространственное разрешение, А0 –
площадь приемной аппаратуры, 𝑇(𝜆𝐿 , 𝑅 ), 𝑇(𝜆, 𝑅) – пропускание атмосферы
на соответствующих длинах волн, (
𝑑𝜎
𝑑𝛺
) - дифференциальное сечение КРС
исследуемой молекулой на длине волны лазерного излучения, Na –
концентрация молекул.
Как указывалось выше, для газоанализатора выбраны две зондирующие
длины волны 266 нм и 785 нм, энергетический расчет проводится для
каждого лазера в качестве регистрируемых газов берется реперный газ – азот
и основной газ для определения – метан. Для проведения расчета в качестве
источников и приемника излучения брались существующие модули: Nd:YAG
лазер на 4-ой гармонике от компании Quantum Composers из линейки
39
SpectroLite Series (таблица 1), 785 нм лазер Cobolt 06-01 Series и CCD
линейка Hamamatsu S10420.
Таблица 2.1 – параметры лазера на длине волны 266 нм
Параметр
Значение
Тип
Nd:YAG
Режим работы
Импульсный
Длительность импульса
4 нс
Частота следования импульсов
1 Гц
Расходимость пучка
< 2 мрад
Энергия в импульсе 4 гармоники, 266 нм
10 мкДж
Таблица 2.2 – параметры лазера на длине волны 785 нм
Параметр
Значение
Тип
DPSS
Расходимость пучка
< 2 мрад
Выходная энергия излучения
250 мВ
Таблица 2.3 – параметры ПЗС линейки
Параметр
Значение
Тип
CCD
Количество пикселей
2068*64
Размер пикселя
14 мкм
Эффективность в диапазоне 270 нм – 290 нм
~0.5
Эффективность в диапазоне 1 мкм – 1.1 мкм
~0.15
Метод КР позволяет одновременно регистрировать концентрации
нескольких
веществ,
каждая
молекула
обладает
индивидуальными
особенностями, соответственно и параметры в лидарном уравнении (2.4) для
40
каждого газа будут свои. Дифференциальное сечение колебательного
комбинационного рассеяния для случая линейно поляризованного лазерного
излучения и наблюдения рассеянного назад излучения обеих поляризаций
суммарно, можно определить по формуле [32]:
𝑑𝜎
( ) 𝑗~
𝑑𝛺
16𝜋 4 𝑏𝑗2 𝑔𝑗
𝜆4 (1 − 𝑒𝑥𝑝 (
ℎ𝑐
))
𝜆𝑘𝑇
(𝛼̅𝑗2 +
7 2
𝛾̅ ),
45 𝑗
где bj - амплитуда нулевых колебаний j-ой моды, gj
вырождения,
𝛼̅𝑗
и
𝛾̅𝑗
след
и
анизотропия
(2.5)
- степень ее
тензора
производной
поляризуемости молекулы по нормальной координате qj соответственно, T колебательная температура молекул; k, h - постоянные Больцмана и Планка
соответственно, c - скорость света. При рассмотрении формулы (N) видно,
что она зависит только λ, поэтому опираясь на экспериментальные
результаты исследований сечения поглощения на длине волны 337.1 нм [33],
можно рассчитать сечения поглощения для длин волн 266 нм и 785 нм.
Стоксовая частота i-ой компоненты определяются по формуле:
𝑣𝑖 = 𝑣0−𝑣кв .
(2.6)
Таблица 2.4 – Параметры молекул при зондировании на длине волны 266 нм
Частота
Молекула
собственных
колебаний
(
λкр, нм
∗ 10
молекул, см-1
Азот N2
Метан CH4
𝑑𝜎
)
𝑑𝛺 337.1 нм
−30
см−2
ср
(
𝑑𝜎
)
𝑑𝛺 266 нм
∗ 10
−30
см−2
ср
2330,7
283,58
3.5
9,02
v1 - 2914
288,35
32,2
82,96
v3 - 3017
289,20
14
36,06
41
Таблица 2.5 – Параметры молекул при зондировании на длине волны 785 нм
Частота
Молекула
собственных
колебаний
(
λкр, нм
∗ 10
молекул, см-1
Азот N2
Метан CH4
𝑑𝜎
)
𝑑𝛺 337.1 нм
−30
см−2
ср
(
𝑑𝜎
)
𝑑𝛺 785 нм
∗ 10
−31
см−2
ср
2330,7
960,79
3.5
1,19
v1 - 2914
1017,83
32,2
10,93
v3 - 3017
1028,61
14
4,76
При зондировании атмосферы на результаты оказывают сильное
влияние ее параметры, их можно разделить на несколько групп:
1. Ослабление
за
счет
селективного
поглощения
газовыми
компонентами и рассеяние на молекулах и аэразольных частицах
2. Искажения из-за турбулентных искажений в воздухе, при этом
создается случайная пространственно-временная структура поля
показателя преломления атмосферы, ввиду чего меняется фаза и
амплитуда лазерного излучения.
3. Нелинейные эффекты, которые вносят случайные и системные
ошибки, зависящие от длины волны источника, параметров
оптического тракта, длины трассы и др.
Полная энергия молекулы имеет набор возможных дискретных
значений, которые соответствуют различным электронно-колебательнымвращательным состояниям.
∆𝐸𝑖𝑗 = 𝐸𝑖 − 𝐸𝑗 = ℎ𝜈𝑖𝑗 ,
(2.7)
где Ei и Ej - начальное и конечное значения энергии соответственно, h постоянная
Планка.
Данная
формула
является
решением
уравнения
Шредингера и определяет разрешенные значения E. Этим значениям
соответствуют частоты νij, определяющие положения спектральных линий
42
молекулы, таким образом формируется спектр поглощения молекулы. В
Земной атмосфере реальные спектральные линии поглощения имеют
определенный контур, полуширина которого зависит от параметров
атмосферы (давления и температуры).
Рисунок 2.4 – Спектр поглощения атмосферы Земли: а — на уровне земной
поверхности, б — на высоте 11 км [28]
На
рисунке
видно,
что
в
ультрафиолетовой
области
спектра
поглощение лазерного излучения земной атмосферой определяется в
основном кислородом и озоном. Кислород имеет интенсивные полосы
поглощения в ультрафиолетовых областях спектра: 0,17-0,19 мкм и 0,24-0,26
мкм. В рамках данной работы интересен диапазон 0,26-0,3 мкм, где
преимущественный вклад в поглощение вносит озон.
Главная полоса поглощения излучения озоном, называемая полосой
Хартли, простирается от 0,2 до 0,3 мкм. В этой электронной полосе
поглощения на общий континуум накладываются многочисленные слабые
43
дополнительные полосы, отстоящие друг от друга по длине волны примерно
на 0,001 мкм [27].
Таблица 2.6 - Сечения поглощения озона при 298 К [30]
λ, нм
σ, 10-20
см2
λ, нм
σ, 10-20
λ, нм
см2
σ, 10-20
см2
λ, нм
σ, 10-20
см2
260,0
1080,0
270,0
798,0
280.0
400.1
290.0
136.5
260,5
261,0
261,5
262,0
262,5
263,0
263,5
264,0
264,5
265,0
265,5
266,0
266,5
267,0
267,5
268,0
268,5
269,0
269,5
1114,0
1094,0
1057,0
1057,0
1066,0
1022,0
995,9
1006,0
1007,0
965,7
965,0
948,5
922,1
884,1
893,8
875,4
836,3
810,4
815,1
270,5
271,0
271,5
272,0
272,5
273,0
273,5
274,0
274,5
275,0
275,5
276,0
276,5
277,0
277,5
278,0
278,5
279,0
279,5
763,6
741,5
726,6
714,7
689,9
669,8
646,0
614,0
596,8
591,3
560,6
545,0
522,1
509,6
490,6
466,8
451,6
432,6
415,8
280.5
281.0
281.5
282.0
282.5
283.0
283.5
284.0
284.5
285.0
285.5
286.0
286.5
287.0
287.5
288.0
288.5
289.0
289.5
372.7
367.3
340.4
325.0
312.2
302.5
296.3
271.2
259.8
246.5
232.2
223.8
210.3
203.4
191.2
175.0
167.3
158.5
142.9
290.5
291.0
291.5
292.0
292.5
293.0
293.5
294.0
294.5
295.0
295.5
296.0
296.5
297.0
297.5
298.0
298.5
300.0
300.5
129.9
122.5
111.8
104.3
99.69
93.75
86.37
81.48
78.24
72.70
65.98
62.02
58.65
54.69
50.02
47.14
44.49
36.16
33.96
Коэффициент поглощения озона (км-1) может быть рассчитан по
формуле [35]:
𝑘 (𝜆) = 𝐴𝑣 (𝜆)𝐷,
(2.8)
где 𝐴(𝜆) — экспоненциальный коэффициент поглощения озона (см -1), D –
эквивалентная концентрация озона, D=3,56*10-3 см/км [35].
Таблица 2.7 - эквивалентный коэффициент поглощения озона [35]
λ (мкм)
𝐴𝑣 , см−1
0,27
2,10*102
0,28
1,06*102
44
0,30
1,01*101
0,32
8,98*10-1
0,34
6,40*10-2
0,36
1,80*10-3
Общее ослабление за счет поглощения озоном будет описываться
законом Бугера [29]:
𝐼 (𝜆) = 𝐼0(𝜆)𝑒𝑥𝑝(−𝑘 (𝜆)𝐿),
(2.9)
где k(λ) —объемный коэффициент поглощения излучения озоном; I0(λ) и I(λ)
—
соответственно
интенсивность
излучения
зондирования
и
регистрируемого, а L длинна трассы.
Помимо молекулярного поглощения излучение ослабевает за счет
молекулярного
(Релеевского)
и
аэрозольного
рассеяния.
Показатель
молекулярного рассеяния σR имеет зависимость от длины волны [42]:
8𝜆2 (𝑛2 − 1)2 6 + 3𝛿
𝜎𝑅 (𝜆) =
,
3𝑁𝜆4
6 − 7𝛿
(2.10)
где λ — длина волны излучения; n — показатель преломления воздуха; N —
число молекул в единице объема; δ = 0,035 — фактор деполяризации
рассеянного излучения.
Для расчетов часто применяют эмпирическую формулу[30]:
𝜎𝑅 = 0.0119км
где
величина
0.0119км−1
−1
0,55 4
∗(
) ,
𝜆
соответствует
показателю
(2.11)
молекулярного
рассеяния в приземном слое атмосферы на длине волны 0,55 мкм, поэтому
для расчета по данной формуле длина волны берется в микрометрах.
Вклад молекулярного рассеяния на трассе L можно оценить законом
Бугера:
𝐼 (𝜆) = 𝐼0(𝜆)𝑒𝑥𝑝(−𝜎𝑅 (𝜆)𝐿).
45
(2.12)
Рисунок 2.5 - Зависимость показателя ослабления излучения вследствие
молекулярного рассеяния и поглощения [49]
Аэрозольное рассеяние вносит серьезный факт в ослабление лазерного
излучения, но данных измерения показателя аэрозольного рассеяния в
ультрафиолетовой
области
спектра
мало
[29].
Значения
показателя
аэрозольного рассеяния земной атмосферой в спектральной области 0,2. . .
0,5 мкм лежат (для различных атмосферных условий) в диапазоне единицы
— десятые доли километра в минус первой степени [28]. Для расчетов можно
использовать эмпирическое выражение [28]:
0,55 𝑚𝑎
𝜎𝑎 = 𝜎𝑎0,55 мк ∗ (
𝜆
)
,
(2.13)
где ma — показатель Ангстрема: для земной атмосферы над сушей ma = 0,92,
среднее значение для открытого океана ma = 0,55, для прибрежных районов
ma = 0,7.
𝜎𝑎0,55 мк =
3,91
𝑆м
− 0,0119,
где Sm – метеорологическая дальность в км.
46
(2.14)
Таблица 2.8 – Параметры ослабления в атмосфере и спектральная постоянная
на исследуемых длинах волн при зондировании 266 нм
𝑘𝑚𝑜𝑙 (𝜆),
𝜎𝑅 (𝜆),
𝜎𝑎 (𝜆),
𝑘 (𝜆),
км−1
км−1
км−1
км−1
266
0,833
0,218
0,270
1,321
283,58
0,304
0,168
0,255
0,727
288,35
0,171
0,158
0,251
0,580
289,2
0,154
0,156
0,250
0,560
𝜆, нм
В диапазоне 1 мкм не наблюдается сильных полос поглощения,
поэтому
для
расчетов
можно
ограничится
только
коэффициентами
аэрозольного и молекулярного рассеяния.
Таблица 2.9 – Параметры ослабления в атмосфере и спектральная постоянная
на исследуемых длинах волн при зондировании 785 нм
𝜆, нм
𝜎𝑅 (𝜆), км−1
𝜎𝑎 (𝜆), км−1
𝑘 (𝜆), км−1
785
0,003
0,132
0,135
960,79
0,001
0,11
0,111
1017,83
0,001
0,104
0,105
1028,61
0,001
0,103
0,104
Таблица 2.10 – Концентрации метана для расчета
Молекула
Метан: CH4
ПДК,
Молекулярная
мг/м
масса, г/моль
3
300
16
Nпдк, см-3
N2ppm, см-3
1,13*1016
5*1013
Прием лазерного сигнала всегда производится на фоне шумов, которые
ограничивают предельные возможности детектировании. Оценить фоновые
шумы можно по формуле [31]:
47
𝑃ф = 𝐿𝜆 𝐾2𝜋𝑟 2𝛺∆𝜆,
(2.15)
где K2 - коэффициент пропускания оптической системы приемного канала
определяется выражением = K1ξλ где K1 - коэффициент пропускания
оптической системы, ξλ – спектральная чувствительность ФПУ на длине
волны λ, Ω=πα2 - телесный угол поля зрения приемной оптической системы,
∆λ – полоса пропускания. Параметр 𝐿𝜆 – спектральная яркость, в случае
горизонтального зондирования может быть рассчитана по формуле [32]:
𝐿𝜆 = 0,25𝜆𝑠 𝑆𝜆 (𝜆)𝜒(𝛾)
cos 𝜃0
𝜏0 − 𝜏
𝜏0 − 𝜏
(exp (−
) − exp (−
)),
cos 𝜃 − cos 𝜃0
cos 𝜃
cos 𝜃0
(2.16)
(2.17)
𝑐𝑜𝑠𝛾 = − cos 𝜃 cos 𝜃0 + sin 𝜃 sin 𝜃0 cos(𝜑−𝜑0 ),
где 𝑆𝜆 (𝜆) – спектральная солнечная постоянная на длине волны λ, 𝜆𝑠 –
отношение
показателя
рассеяния
аэрозольной
составляющей
земной
атмосферы к показателю ослабления (в видимой области спектра 𝜆𝑠 ≅ 1), 𝜏0оптическая толща всей атмосферы, 𝜏 – оптическая толща атмосферы меду
подстилающей поверхностью и приемником, 𝜃 и 𝜑 – зенитный угол и азимут
направления наблюдения, 𝜃0 и 𝜑0 – зенитный угол и азимут направления
Солнца, 𝜒(𝛾) – индикатриса рассеяния атмосферы, параметр обозначающий
плотность вероятности рассеяния на угол 𝛾.
Индикатриса релеевского рассеяния атмосферы [32]:
3
𝜒𝑅 (𝛾) = (1 + 𝑐𝑜𝑠 2𝛾).
4
(2.18)
Индикатриса аэрозольного рассеяния атмосферы [32]:
𝜒𝑎 (𝛾) =
1−𝑔2
3
(1+𝑔 2 −2𝑔𝑐𝑜𝑠𝛾)2
,
(2.19)
где g- параметр характеризующий вытянутость индикатрисы 0 ≤ 𝑔 ≤ 1
Оптическая толща атмосферы безразмерный параметр, который
складывается из составляющих толщь: τR - молекулярного и τa - аэрозольного
рассеяния, а также τmol - молекулярного поглощения.
Как указывалось выше в случае УФ области поглощение идет
преимущественно за счет озона, поэтому толщу молекулярного поглощения
48
можно рассчитать с использованием не экспоненциального, а десятичного
объемного коэффициента поглощения [28]:
𝜏𝑚𝑜𝑙 = 𝑘(𝜆)𝑥,
(2.20)
где 𝑘(𝜆) – десятичный объемный коэффициент поглощения излучения
озоном, а х – толщина однородного слоя озона по направлению
распространения излучения, х=0,3 см – средний слой озона от нижней до
верхней границы атмосферы.
Рисунок 2.6 - Спектральная зависимость десятичного коэффициента
поглощения излучения озоном в ультрафиолетовой области спектра [28]
Толщи
молекулярного
и
аэрозольного
рассеяния
могут
быть
рассчитаны по эмпирическим формулам [28]:
0.55 4
𝜏𝑅 = 0.0952 (
𝜏𝑎 = 0.2 (
𝜆
) ,
0.55 𝑚𝐴
𝜆
)
,
(2.21)
(2.22)
где λ берется в микрометрах, а ma — показатель Ангстрема: для земной
атмосферы над сушей ma = 0,92 [28].
49
Усредненная по всей земной атмосфере индикатриса рассеяния
определяется формулой [28]:
𝜒(𝛾) =
𝜒𝑅 (𝛾)𝜏𝑅 +𝜒𝑎 (𝛾)𝜏𝑎
𝜏𝑅 +𝜏𝑎
.
(2.23)
Таблица 2.11 - Оптические толщи атмосферы при зондировании 266 нм
𝑘 (𝜆), см−1
𝜆, нм
𝜏𝑚𝑜𝑙
𝜏𝑎
𝜏𝑅
𝜏0
30,4
283,58
6,08
0,367877
1,347055
7,794932
17,1
288,35
3,42
0,362274
1,260108
5,042382
15,4
289,2
3,08
0,361294
1,245359
4,686653
Таблица 2.12 - Оптические толщи атмосферы при зондировании 785 нм
𝜆, нм
𝜏𝑎
𝜏𝑅
𝜏0
960,79
0,119714
0,010223
0,129937
1017,83
0,113528
0,008117
0,121645
1028,61
0,112433
0,007782
0,120215
Таблица 2.13 – Расчет фоновых шумов
𝜆, нм
S(λ),
м2
Вт
∗ мкм
𝑃ф ,∗ 10−15 Вт
𝜆, нм
S(λ),
м2
Вт
∗ мкм
𝑃ф ,∗ 10−12Вт
283,58
0,03
2,64
960,79
0,0803
4,13
288,35
0,038
32,7
1017,83
0,0752
3,64
289,2
0,04
57,6
1028,61
0,0752
3,62
Расчет производился программном комплексе MathCad, для метана
были выбраны две наиболее интенсивные полосы 3017 см-1 и 2914 см-1.
50
Рисунок 2.7 – Зависимость величины полезного сигнала от дальности
зондирования для азота при λ=266 нм и концентрации молекул N=1016 cм-3
Рисунок 2.8 – Зависимость величины полезного сигнала от дальности
зондирования для азота при λ=785 нм, времени накопления 30 мс и
концентрации молекул N=1016 cм-3
51
Рисунок 2.9 – Зависимость величины полезного сигнала от дальности
зондирования для метана при λ=266 нм и концентрации по ПДК: график P(r)
соответствует волновому числу 2914 см-1, P2(r) волновому числу 3017 см-1
Рисунок 2.10 – Зависимость величины полезного сигнала от дальности
зондирования для метана при λ=785 нм, времени накопления 30 мс и
концентрации по ПДК: график P(r) соответствует волновому числу 2914 см-1,
P2(r) волновому числу 3017 см-1
52
Рисунок 2.13 – Зависимость величины полезного сигнала от дальности
зондирования для метана при λ=266 нм и концентрации 2ppm: график P(r)
соответствует волновому числу 2914 см-1, P2(r) волновому числу 3017 см-1
Рисунок 2.14 – Зависимость величины полезного сигнала от дальности
зондирования для метана при λ=785 нм, времени накопления 500 мс и
концентрации 2 ppm: график P(r) соответствует волновому числу 2914 см-1,
P2(r) волновому числу 3017 см-1
53
Выводы по главе
Решение лидарного уравнения задача весьма непростая, итоговые
показатели по дальности детектирования на порядок могут отличаться от
расчетных. Определение уровня шумов является достаточно сложной
задачей, исходные параметры зависят от большого числа факторов, и если
рассмотреть иные погодные условия, то уровни фонового излучения будут
иными, для УФ диапазона это не так значительно, поглощение озоном вносит
настолько сильный вклад, что изменение других параметров атмосферы
несущественно влияют на дальность зондирования.
В ИК диапазоне картина иная, самые интенсивные полосы УВГ лежат в
диапазоне 1015 нм – 1040 нм (приложение 1), где не наблюдается полос
поглощения атмосферными газами. Поэтому ИК канал сильно зависим от
погодных явлений, и дальность зондирования будет меняться от случая к
случаю. Устранить данный изъян можно изменяя время зондирования, так
как ИК канал работает в режиме накопления можно увеличивать, как
дальность зондирования, так и предел детектирования за счет увеличения
времени накопления. Но тут важно, что помимо полезного сигнала
накапливаются и шумы, поэтому при длительном накоплении ПЗС детектор
может перейти в насыщение и потерять разрешающую способность.
По дальности детектирования метана имеется достаточный запас,
поэтому использование для дальнейшего рассмотрения параметры лазеров и
приемной аппаратуры оправданы.
Полученные
при энергетических расчетах графики показывают
высокие параметры детектирования метана, который является основным
газом для регистрации. Также видно, что при зондировании метана хорошо
можно наблюдать не одну его линию, но это является и некоторой
проблемой, влияющей на работоспособность системы. Как отмечалось в
выводах первой главы, у углеводородов самые интенсивные линии
расположены близко друг к другу, при расчетах метана видно, что две
54
близкие интенсивные линии различимы на одинаковом расстоянии. Если
учесть, что у всех УВГ самые интенсивные линии лежат в диапазоне 280 нм 290 нм, нужно предъявлять высокие требования к оптическому тракту, дабы
линии не сливались на детекторе.
При необходимости столь высокой точности использование второго
источника полностью себя оправдывает, к тому же в расстояние между
сдвигами в ИК диапазоне заметно выше чем в УФ, а значит требования к
регистрирующей системе для ИК значительно ниже.
55
3. МОДЕЛИРОВАНИЕ
3.1. Оптический канал
Функциональные
возможности
любого
оптико-электронного
спектрального прибора определяются характеристиками его оптической
схемы. Важно оценить общие потери потока излучения на всех ее элементах,
данные потери имеют различную природу и сложно оценимы, такой
характеристикой оптического тракта является коэффициент пропускания,
численно равный отношению потока монохроматического излучения на его
выходе к потоку на входе:
𝜏=
𝛷𝑠′
𝛷𝑠
,
(3.1)
где 𝛷𝑠′ - поток излучения, измеренный в плоскости входной щели, Φ𝑠 - поток
излучения, измеренный на поверхности изображения.
Одной из основных характеристик оптической системы является
угловая
дисперсия
диспергирующего
устройства,
характеризующая
способность решетки отклонять излучение различных длин волн. Она
определяется как производная угла дифракции по длине волны 𝑑𝜑′/𝑑𝜆.
Линейная дисперсия 𝑑𝑙/𝑑𝜆 характеризуется как отношение расстояния между
двумя близкими спектральными компонентами в изображении 𝑑𝑙 к
соответствующему малому интервалу длин волн 𝑑𝜆 [43].
В задачах спектрометрии важным является различимость спектральных
линий
на
детекторе,
монохроматических
наименьшая
спектральных
разность
линий
длин
равной
волн
двух
интенсивности,
наблюдаемых раздельно, называют спеткральным пределом разрешения.
Также существует линейный предел разрешения, это наименьшее расстояние
между двумя различаемыми монохроматическими линиями в плоскости
изображения спектра. Спектральной разрешающей способностью называют
отношение:
56
𝑅=
𝜆
∆𝜆
(3.2)
.
Плоская отражательная дифракционная решетка является одним из
часто применяемых диспергирующих элементов. Она представляет собой
пластинку, с нанесенным на нее рядом параллельных равноотстоящих
штрихов. Расcтояние между штрихами называют постоянной решетки, как
правило решетку характеризуют числом штрихов на единицу длины:
1
𝑁= .
(3.3)
𝑙
При падении на решетку пучка лучей, параллельного главному
сечению,
на
каждой
полоске
происходит
дифракция.
Пучки,
дифрагированные на отдельных штрихах интерферируют между собой на
плоскости регистрации. Дифракция на решетке описывается известным
уравнением, называемым уравнением решетки:
𝑠𝑖𝑛𝜑 ′ − 𝑠𝑖𝑛𝜑 = 𝑘𝜆𝑁,
(3.4)
где 𝜑 - угол падения излучения на решетку, 𝜑′ - угол дифракции для длины
волны 𝜆, 𝑁 - частота штрихов решетки, 𝑘 - порядок дифракции.
Существует несколько видов дифракционных решеток, но особый
интерес представляют голограммные дифракционные решетки с коррекцией
аберраций, т.к. они отличаются высоким качеством изображения при
минимальных
габаритах,
имеют
достаточно
высокие
относительное
отверстие и спектральное разрешение [44]. При этом сохраняются такие
преимущества, как конструктивная простота, надежность функционирования
и технологичность в изготовлении.
Голографические решетки представляют собой зарегистрированную на
светочувствительном материале картину интерференции двух световых
пучков. Параметры голографических решеток можно изменять в широком
диапазоне с помощью схемы записи и формы поверхности, на которой
регистрируется решетка.
Угловая дисперсия:
57
𝑑φ
dλ
=
k
,
b·cos φ
(3.5)
где b – период решетки, 𝜑 - угол дифракции для длины волны 𝜆, 𝑘 - порядок
дифракции.
Предельная разрешающая способность:
𝜆
∆𝜆
=
𝑁𝑏
𝜆
(3.6)
.
3.2 Телескоп
Оптический канал начинается с телескопа, существуют рефракторные
(линзовые) и рефлекторные (зеркальные) телескопы. По сравнению с
рефракторами рефлекторные телескопы полностью лишены хроматической
и сферической аберрации. В приемной
системе лидаров наиболее часто
применяют рефлекторные схемы Ньютона (рисунок 3.1 а) и Кассегрена
(рисунок 3.1 б) [45].
Рисунок 3.1 – Схемы телескопов а – Ньютона и б – Кассегрена
При одинаковых размерах главных зеркал схема Кассегрена имеет
значительный выигрыш по длине, но астигматизм в схеме Кассегрена больше
58
астигматизма
Ньютоновского
параболоида.
Также
схема
Ньютона
значительно более проста и дешева в изготовлении, что делает ее более
приемлемой для массового производства. Выигрыш по длине стоит
рассматривать применительно ко всей системе, учитывая стыковку с
полихроматором, тогда при использовании схемы Кассегрена увеличиваются
продольные размеры системы, а поперечные меняются не значительно [43].
Для
дальнейшего
рассмотрения
выбрана
схема
Ньютона,
для
моделирования ее геометрии можно воспользоваться готовым модулем в
программном комплексе Comsol. Главной характеристикой является входная
апертура телескопа, поскольку задачей стоит рассмотрение портативного
прибора апертура бралась равной 5 см.
В модели лучи света распространяются от источников, расположенных
на бесконечности, в телескоп. Лучи отражаются от параболического зеркала,
на плоское зеркало, и попадают в фокальную плоскость, формируя
изображения объектов-источников (рисунок 3.2).
Таблица 3.1 - Параметры телескопа
Параметр
Значение
Диаметр главного зеркала
50 мм
Фокусное расстояние главного зеркала
200
Толщина главного зеркала
7 мм
Относительное отверстие
1:4
Диаметр парабалического зеркала
30 мм
Толщина парабалического зеркала
3 мм
Положение плоскости изображения
(относительно оптической оси)
59
40 мм
Рисунок 3.2 – Лучевая диаграмма телескопа
3.3 Монохроматор
Как уже было отмечено выше отражательная дифракционная решетка
на сегодняшний день является самым распространенным диспергирующим
элементом в спектральных приборах различного назначения. Решетка
является пусть и главным, но не единственным элементом системы
регистрации оптического сигнала. Монохроматор может быть выполнен по
различным схемам, в зависимости от используемого типа дифракционной
решетки, самыми распространенными схемами являются: схема ПашенаРунге, основанная на использовании вогнутой решетки, на ней изготовлено
большинство приборов с фокусным расстоянием от 500 мм и больше. При
разработке малогабаритных приборов (фокусное расстояние около 100 мм),
как правило, используют схему Черни-Тернера, основанную на плоской
решетке (в некоторых случаях схему Черни-Тернера строят и на вогнутых
решетках).
поверхность
Схему Черни
регистрации,
- Тернера
отличает протяженная
обеспечивающая
одинаковое
разрешение во всем рабочем спектральном диапазоне [46].
60
плоская
спектральное
Рисунок 3.3 - Ассиметричная скрещенная схема Черни-Тернера, s расстояние между входной щелью и коллиматорным зеркалом, θ0 – угол
падения главного луча, θi - угол дифракции, l1 - расстояние от решетки до
фокусирующего зеркала, l2 - а расстояние от фокусирующего зеркала до ПЗС
детектора [47]
Использование ассиметричной скрещенной схемы Черни-Тернера на
основе плоской дифракционной решетки помогает устранить рассеянный
свет
и
добиться
минимальных
размеров
монохроматора.
Типичная
оптическая схема асимметричного скрещенного спектрометра Черни –
Тернера показано на рисунке 3.3. Оптические элементы включают в себя
входную щель, сферические коллиматорное и фокусирующее зеркало,
плоскую дифракционную
решетку и ПЗС матрицу для регистрации
оптических сигналов.
В этой оптической схеме световой пучок попадет в пространство
монохроматора через входную щель, а затем его собирает коллимирующее
зеркало. На коллимирующем зеркале формируется плоский волновой фронт
и направляется на решетку. Отражаясь от решетке, под несколькими углами,
61
зависящими от параметров решетки и длины волны, лучи фокусируются на
ПЗС детекторе фокусирующим зеркалом.
Рассмотрим геометрическую модель схемы Черни-Тернера (рисунок
3.4), модель создана таким образом, что центр дифракционной решетки
находится в точке 0, световой луч выходит из входной щели и попадает на
коллимирующее зеркало с центром в точке Qc, наклоненное под углом θc.
После коллимирующего зеркала луч падает на решетку, повернутую на угол
θg, под углом α. На решетке происходит дифракция, угол β – угол дифракции,
луч падает на фокусирующее зеркало, попадая в точку Pi, центр
фокусирующего зеркала находится в точке Qi под углом θi. Затем луч
приходит в точку Pd на ПЗС детекторе, с центром в точке Qd под углом θd.
Рисунок 3.4 – Геометрическая модель скрещённой схемы Черни-Тернера
62
Дифракция света на решетке определяется уравнением решетки:
𝑑(𝑠𝑖𝑛𝛼 + 𝑠𝑖𝑛𝛽 ) = 𝑘𝜆,
(3.7)
где k порядок дифракции, должен быть целым числом; λ - длина волны, α угол падения света на дифракционную решетку, а β - угол дифракции.
Итоговые параметры монохроматора зависят от его геометрии:
размеров входного отверстия, размеров и углов поворота зеркал, а также от
расстояния между всеми элементами. Также влияние оказывают параметры
дифракционной решетки: ее постоянная и размеры заштрихованной области,
а также параметры ПЗС детектора: размер пикселей, их количество и общий
размер рабочей поверхности.
Используемая
в
работе
модель
является
готовым
модулем
в
программном комплексе Comsol, его основной задачей является построение
взаимосвязи между геометрическими параметрами системы. Поэтому при
моделировании мы предполагаем, что падающий луч представляет собой
плоский волновой фронт, а не сферический, за счет чего не возникает
аберрации, все размеры намного больше, чем длина волны; это позволяет нам
заменить волны пучком и рассмотреть распространение света как траектории
луча. Ввиду принятых выше допущений итоговые параметры монохроматора
определяются только геометрией схемы.
Исходная концепция прибора заключалась в том, чтобы проводить
спектральный анализ на одном приемнике, но при таком условии даже в
теории получить высокое разрешение анализируя полосу 840 нм (от 260 нм
до 1100 нм) невозможно. Поэтому было принято решение достаточно сильно
развести оптические элементы между собой, таким образом, чтобы можно
было использовать две матрицы: с одной вести анализ в УФ диапазоне, а
другой ИК, для реализации данного метода готовый модуль нужно было
изменить введением в него второго детектора.
Рамановские спектрометры требовательны к светосиле, поэтому
относительное оптическое отверстие бралось 1:4, главная задача получить
63
высокое спектральное разрешение, для этого обратимся к таблице
рамановских сдвигов (приложение 1), откуда можно сделать вывод: для УФ
диапазона (270 нм – 300 нм) необходимое разрешение менее 0,08 нм, а для
диапазона ИК (850 нм – 1050 нм ) менее 1 нм. Параметры детекторов брались
как и при энергетических расчетах для ПЗС линейки Hamamatsu S10420.
Итоговые результаты представлены на рисунке 3.5.
Общие размеры схемы получились 120 мм на 120 мм, что приемлемо
для портативного прибора, но его спектральные характеристики далеки от
требуемых, в чем можно убедится, посмотрев на графики: рисунок 3.6 и 6.7.
Таблица 3.1 – Параметры дифракционной решетки и детектора
Параметр
Значение
Относительное отверстие
1:4
Постоянная решетки
700 мм-1
Размер решетки
30 мм
Количество пикселей детектора
2064
Размер пикселя
14 нм
Рисунок 3.5 – Спектральное разрешение первой матрицы
64
Рисунок 3.6 - Монохроматор Черни-Тернера с двумя приемниками
излучения, спектральный диапазон 259 нм – 1110 нм
Рисунок 3.7 – Спектральное разрешение второй матрицы
65
В ИК диапазоне спектральные характеристики получились пусть и не
достаточные, но близкие к требуемым, возможно немного подстроив
параметры, можно добиться нужных значений. А в УФ диапазоне параметры
слишком далеки от необходимых, любые попытки улучшить их за счет
изменений углов элементов или их положений успехом не увенчались. К
тому же при такой конфигурации схема в общем становится слишком
сложной, пытаясь получить высокие параметры на одном из приемников,
сбиваются параметры на другом, кажется что получить по такой схеме
одинаково высокое разрешение в обоих диапазонах просто не возможно.
Единственным
возможным
вариантом
получения
высокого
спектрального разрешения в УФ диапазоне остается увеличение числа
штрихов на решетке. Используя схему, представленную на рисунке 3.5,
максимальное количество штрихов в см ограничивается приблизительно
1000 мм-1, уже при значении 700 мм-1 размеры фокусирующего зеркала
достаточно большие, а при значениях свыше 700 мм -1 углы дифракции на
столько большие, что сфокусировать их просто невозможно.
Отказываться от концепции одного монохроматора решительно не
хотелось, но при текущей постоянной решетки это невозможно. Поэтому
необходимо использовать другую схему монохроматора, реализовать схему с
постоянной решетки 1500 мм-1 можно ипользуя схему Пашена-Рунге.
Схематически данная схема представлена на рисунке 3.8, главное ее
правило заключается в том, что центр вогнутой решетки, щель и все ее
монохроматические изображения находятся на поверхности кругового
цилиндра радиусом в два раза меньше радиуса решетки, в двумерном
представлении цилиндр является кругом, называемым кругом Роуланда [43].
66
Рисунок 3.8 – Схема Пашена-Рунге, R - радиус круга Роуланда, r –
расстояние между входной щелью и дифракционной решеткой, s –
расстояние от решетки до изображения, α – угол падения, β – угол дифракции
Для регистрации изображения при использовании схемы ПашенаРунге необходимо использовать дополнительные коллимирующие линзы, в
целом для того, чтобы провести оценку спектрального разрешения делать
столь сложную модель не обязательно, попробуем для начала создать модель
на круге Роулонда и оценить спектральное разрешение в фокусе. Готовых
модулей для данной задачи не существует, но несложно преобразовать
модуль со схемой Черни-Тернера в схему Пашена-Рунге, убрав лишние
зеркала
и
изменив
геометрию
дифракционной
решетки.
моделирования представлены на рисунках 3.9, 3.10 и 3.11.
Таблица 3.2 – Параметры дифракционной решетки и детектора
Параметр
Значение
Относительное отверстие
1:4
Постоянная решетки
1500 мм-1
67
Результат
Продолжение таблицы 3.2
Размер решетки
30 мм
Радиус решетки
200 мм
Количество пикселей детектора
2064
Размер пикселя
14 нм
Рисунок 3.9 – Монохроматор Пашена-Рунге спектральный диапазон 258 нм 1120 нм
68
Рисунок 3.10 – Геометрия схемы Пашена- Рунге
Рисунок 3.11 – Спектрально разрешение в диапазоне 270 нм - 300 нм
69
При моделировании второй ПЗС детектор не вводился намерено, так
как задачей стояло оценить спектральное разрешение именно в УФ
диапазоне, чтобы понять имеет ли смысл более углубленно рассмотреть
схему Пашена-Рунге. Предельная постоянная решетки, которую удалось
получить это 1500 мм-1, при дальнейшем увеличении углы дифракции
начинают перекрывать входную щель. При радиусе 100 мм и спектральном
диапазоне 30 нм приближение о том, что все длины волн фокусируются на
прямом
детекторе,
не
вносит существенной погрешности.
Итоговое
разрешение получилось заметно выше чем, при использовании схемы ЧерниТернера, но все же далеким от требуемых значений. Требуется дальнейшее
увеличения числа штрихов на решетке, а при этом наблюдать столь широкий
спектральный диапазон физически невозможно. Поэтому приходится
заключить, что для поставленной задачи
одного монохроматора не
достаточно.
Тогда придется делить входное излучение на два монохроматора и
проводить спектральный анализ каждого диапазона отдельно. Проблемой
при таком анализе является деление излучения, разделить излучение можно
при помощи оптического клина, но при этом происходит потеря сигнала в
два
раза.
Более
подходящими
выглядят
схемы
с
оптоволоконным
переключателем или схема при использовании вращающегося клина. Даная
инженерная задача будет рассмотрена в дальнейшем, при развитии работы, в
рамках данной работы такой задачи не ставится.
С учетом того, что необходимо максимально сохранить портативность
системы,
нужно
использовать
два
монохроматора
Черни-Тренера.
Наибольшее спектральное разрешение нужно получить в УФ диапазоне,
добиться этого попытаемся с использованием дифракционной решетки с
постоянной 3600 мм-1. Моделирование проведем в том же модуле, что и
ранее, но не внося в него изменений.
70
Таблица 3.3 – Параметры монохроматора Черни-Тернера для УФ диапазона
Параметр
Значение
Относительное отверстие
1:4
Постоянная решетки
3600 мм-1
Ширина заштрихованной области
15 мм
Угол поворота дифракционной решетки
19,5°
Размеры коллимирующего зеркала
15 мм на 3 мм
Радиус коллимирующего зеркала
125 мм
Угол поворота коллимирующего зеркала
11°
Размеры фокусирующего зеркала
20 мм * 3 мм
Радиус фокусирующего зеркала
95 мм
Угол поворота фокусирующего зеркала
78°
Количество пикселей детектора
2064
Размер пикселя
14 нм
Угол поворота детектора
1°
71
Рисунок 3.12 - Геометрия монохроматора Черни-Тернера для УФ диапазона
Рисунок 3.13 – Монохроматор Черни-Тернера для УФ диапазона
72
Рисунок 3.14 - Освещенность пикселей на ПЗС детекторе
Рисунок 3.15 – спектральное разрешение для УФ диапазона
Таблица 3.3 – Параметры монохроматора Черни-Тернера для ИК диапазона
Параметр
Значение
Относительное отверстие
1:4
Постоянная решетки
800 мм-1
Ширина заштрихованной области
15 мм
Угол поворота дифракционной решетки
14,5°
Размеры коллимирующего зеркала
15 мм на 3 мм
Радиус коллимирующего зеркала
100 мм
Угол поворота коллимирующего зеркала
11°
Размеры фокусирующего зеркала
35 мм * 3 мм
Радиус фокусирующего зеркала
110 мм
Угол поворота фокусирующего зеркала
78.7°
Количество пикселей детектора
2064
Размер пикселя
14 нм
73
Угол поворота детектора
-2°
Рисунок 3.16 - Геометрия монохроматора Черни-Тернера для ИК диапазона
74
Рисунок 3.17 – Монохроматор Черни-Тернера для ИК диапазона
75
Рисунок 3.18 - Освещенность пикселей на ПЗС детекторе
Рисунок 3.19 – спектральное разрешение монохроматора для ИК диапазона
Выводы по главе
Использование двух монохроматоров дает нужные плоды: в УФ
диапазоне
спектральное
разрешение
получилось
более
чем
удовлетворительным, особенно стоит обратить внимание на диапазон 280 нм
- 290 нм где спектральное разрешение порядка 0,028 нм, стоит напомнить что
именно в данном диапазоне лежат самые интенсивные рамановские линии
исследуемых углеводородов. Разрешение в ИК диапазоне, также превосходит
ожидания, основные сдвиги лежат в диапазоне 1015 нм – 1040 нм,
спектрального разрешения 0,6 нм достаточно для того, чтобы эффективно
регистрировать сигналы. Также стоит отметить, что оба монохроматора
получились достаточно компактными 65 мм на 55 мм.
Но существует и обратная сторона медали, самое важное это то, что
придется разводить сигнал на два монохроматора, соответственно вдвое
76
терять полезный сигнал. Также использование двух монохроматоров
усложняет оптическую систему в целом, придется добавлять дополнительные
элементы в схему, из-за чего в схему вносятся дополнительные аберрации, и
как итог уменьшается пропускная способность. Рамановский сигнал и так
крайне слаб на фоне общего излучения, поэтому его потеря в оптическом
тракте может существенно сказать я на конечных параметрах спектрометра.
Схема Пашена-Рунге примечательна тем, что имеет малое количество
оптических элементов, а значит меньше рассеянного света в схеме, но
получить малогабаритный прибор с высоким разрешением используя эту
схему сложно.
Рисунок 3.20 - S – образная схема монохроматора, 1 – входная щель, 2
– вогнутая дифракционная решетка, 3 – дифрагирующие лучи, 4 –
фокусирующее зеркало, 5 – фотодетектор [44]
Решить данную проблему можно за счет введения дополнительных
оптических компонентов. За счет использования двух отражений и большего
числа свободных коррекционных параметров в подобной схеме удается
77
достичь высокой светосилы и относительно высокого качества изображения
при минимальных габаритах [48]. Схема, представленная на рисунке,
получила название S-образной, подробно рассматривается в работах [44] и
[48], из-за меньшего числа оптических элементов чем в схеме Черни-Тернера
в схеме меньше рассеянного света, а это существенно для спектроскопии КР,
поэтому в дальнейшем стоит рассмотреть данный тип монохроматора в
качестве альтернативы схемы Черни-Тернера.
78
4. БЕЗОПАСНОСТЬ ЖИЗНЕДЕЯТЕЛЬНОСТИ
4.1 Общие положения по использованию лазеров в рамках данной
работы
Дипломная работа является лишь первым этапом в цепочке разработки
описываемого прибора, и, несмотря на то, что в рамках данной работы
экспериментов с
лазерными установками не
производилось,
вопрос
рассмотрения безопасности работы с лазерами остается актуальным для
дальнейшей работы над проектом.
В главе 2 проводился анализ перспективных для использования
лазерных модулей, предполагается использовать два лазерных модуля:
первый, работающий в ультрафиолетовой области спектра, Nd:YAG лазер на
4-ой гармонике, с длинной волны 266 нм и энергией в импульсе 10 мкДж,
второй, работающий в инфракрасной области на длине волны 785 нм и
мощностью 250 мВт.
4.2 Анализ опасностей при работе с лазерами
Лазеры это огромное количество различных приборов с широким
диапазоном характеристик, в вопросах обеспечения безопасности при работе
с лазерами особый интерес представляет спектральный диапазон, так как
чувствительность органов человека к лазерному излучению различных длин
волн разная. По рабочему спектральному диапазону различают приборы
1. Вакуумного ультрафиолета 1-185 нм
2. Ближнего ультрафиолета 185-400 нм
3. Видимой области 400-700 нм
4. Ближней инфракрасной области 0,7-2,5 мкм
5. Средней инфракрасной области 2,5-50 мкм
6. Дальней инфракрасной области 50-1000 мкм
79
Помимо спектрального диапазона на характер и величину нанесенного
вреда оказывают влияние следующие факторы: направление лазерного луча
по
отношению
к
человеку,
длительность
импульса
излучения,
пространственное распределение энергии в луче. Также на степень
поражения лазерным излучением влияют индивидуальные особенности
чувствительности органов каждого человека.
На
организм
человека
оказывает воздействие лазерное излучение любой длины, но в основном
принято выделять воздействие на кожу и органы зрения. Так же разделяют
прямое воздействие лазерного излучения и отраженное от зеркальных
поверхностей.
При попадании лазерного излучения в глаз энергия лазерного пучка
поглощается пигментным эпителием, за счет чего течение очень короткого
времени
повышается
происходит
вызывая
внутриглазная
температура,
термокоагуляция
и,
прилегающих
как
результат,
тканей
или
хориоретинальный ожог. Повреждения центральной ямки особенно опасны,
так как нарушения в ее функционировании могут привести к критическим
нарушениям центрального зрения.
Поражения кожи индивидуальны, так как зависят от пигментации и
состояние кровообращения. Светлая кожа более проницаема для лазерного
излучения, а пигментация увеличивает поглощение излучения на коже.
Основными нормативными документами, регулирующими правила
эксплуатации, лазеров являются: ГОСТ 12.1.040-83, ГОСТ Р 50723-94, ГОСТ
12.0.003-2015 и СанПиН 5804-91.
Согласно ГОСТ Р 50723-94 выделяют
следующие классы лазерных изделий по степени опасности:
1. К первому классу относят лазерные изделия полностью безопасные
при предполагаемых условиях эксплуатации.
2. Второй класс это лазерные изделия, генерирующие видимое
излучение в диапазоне длин волн от 400 до 700 нм. Защита глаз
80
обеспечивается
естественными
реакциями,
включая
рефлекс
мигания.
3. Третий класс делят на подклассы А и Б, лазеры подкатегории А
работают в диапазоне от 400 до 700 нм защита обеспечивается
естественными реакциями, но наблюдение таких лазеров с помощью
оптики опасно. Подкласс Б опасен при непосредственном видении
излучения, но рассеянное излучение безопасно
4. Четвертый класс это лазеры, рассеянное излучение которых
представляет опасность для органов человека.
4.3 Предельно допустимые уровни лазерного излучения при
воздействии на глаза и кожу
Предельно
допустимые
уровни
(ПДУ)
лазерного
излучения
устанавливаются для двух условий облучения - однократного и хронического
для трех диапазонов длин волн:
1. 180 нм < λ < 380 нм
2. 380 нм < λ < 1400 нм
3. 1400 нм < λ < 105 нм
Нормируемыми
энергетическая
параметрами
экспозиция
H
и
лазерного
облученность
излучения
E,
являются
усредненные
по
ограничивающей апертуре.
В импульсном режиме планируется использовать ультрафиолетовый
лазер, рассчитаем предельно допустимые параметры при однократном
воздействии на глаза и кожу коллимированного или рассеянного излучения в
диапазоне 180 нм < λ < 380 нм. Ограничивающая апертура в данном случае
принимается равной 1,1*10-3 м, длительность импульса берем согласно
параметрам лазера (таблица 2.1).
𝐻пду = 25 Дж ∗ м−2
81
𝐸пду =
25
25
=
= 6,25 ∗ 109 Вт ∗ м
−9
𝑡
4 ∗ 10
𝑊пду = 𝐻пду ∗ 10−6 = 25 ∗ 10−6 Дж
𝑃пду = 𝐸пду ∗ 10−6 = 6,25 ∗ 103 Вт
Теперь произведем расчеты для ИК лазера, он лежит во втором
диапазоне, его планируется использовать в непрерывном режиме (время
накопления 30 мс). При расчете хронического воздействия на глаза
соответствующие величины, рассчитанные для одиночного воздействия,
уменьшаются на 10. При воздействии на кожу:
5
𝐻пду
2 ∗ 104 5√𝑡 2 ∗ 104 √30 ∗ 10−3
=
=
= 9,92 ∗ 103 Дж ∗ м−2
10
10
𝑊пду = 𝐻пду ∗ 10−6 = 9,92 ∗ 10−3 Дж
При воздействии на глаз:
3
3
3 ∗ 10−4 √𝑡 2 3,0 ∗ 10−4 √(30 ∗ 10−3)2
𝑊пду =
=
= 2,9 ∗ 10−6 Дж
10
10
Из проведенных выше расчетов следует, что рассматриваемые в работе
лазерные источники не превышают предельно допустимых параметров по
мощности излучения. Тем не менее, не стоит подвергать себя лишней
опасности при работе с лазерными установками и соблюдать меры
предосторожности.
4.4 Меры и средства защиты от лазерного излучения
Для предотвращения пагубного воздействия лазерного излучения на
органы человека существует целый ряд мер инженерно-технического,
планировочного, организационного, санитарно-гигиенического характера.
Повышенные меры безопасности нужно применять к лазерам 2-го
класса, рассматриваемые в работе модули принадлежат ему, необходимо
экранировать зону рабочей зоны лазеров.
82
К самым серьезным ограничениям нужно прибегать используя лазеры
четвертого класса, для обеспечения полной безопасности их размещают в
отдельных изолированных помещениях и обеспечивают дистанционным
управлением их работой.
К индивидуальным средствам защиты, обеспечивающим безопасные
условия труда при работе с лазерами, относятся специальные огни, щитки,
маски, снижающие облучение глаз до ПДУ.
83
ЗАКЛЮЧЕНИЕ
В
работе
проводилось
исследование
перспектив
использования
приборов на основе метода рамановской спектроскопии при разведке нефти
и газа. Проведен анализ иных физических принципов, на которых строят
приборы для газового анализа основным преимуществом метода КР является
высокая селективность при одноразовом измерении и чувствительность,
необходимая для регистрации УВГ аномалий.
Метод газового анализа не заменяет существующие методы разведки
нефти и газа, а дополняет их, с целью снижения затрат на производство.
Поэтому важна портативность прибора для проведения исследований на
местности.
В работе предложена концепция прибора, использующая два лазерных
источника на длинах волн 263 нм и 785 нм, данная конфигурация позволить
заметно повысить точность прибора, но, в свою очередь, повышает
сложность системы. При выбранных параметрах лазеров и приемника
получилось добиться высоких характеристик детектирования при расчете
лидарного уравнения.
К
сожалению,
добиться
высокого
разрешения
системы
при
использовании одного монохроматора не получилось, было сделано две
попытки моделирования: по схеме Черни-Тернера и Пашена-Рунге, обе
успехом не увенчались. Поэтому пришлось отказаться о концепции одного
монохроматора
для
двух
спектральных
диапазонов,
и
провести
моделирование отдельно.
Разработка спектрометра достаточно сложная задача, в данной работе
освещены лишь первые аспекты будущего прибора, разработка которого
станет предметом будущих научных работ.
84
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Аблязов Э. К., Шеманин В. Г. Решение лидарного уравнения для
контроля углеводородов в атмосфере //Научно-технические ведомости
Санкт-Петербургского
государственного
политехнического
университета. Физико-математические науки. – 2009. – №. 2 (77).
2. Бажанов Ю. В. и др. Количественный анализ газовых сред методом
спектроскопии комбинационного рассеяния света //Аналитика и
контроль. 1998.№ 3/4. – 1998. – С. 65-74.
3. Шаповалова Е. Н., Пирогов А. В. Хроматографические методы анализа
//Методическое пособие для специального курса. М. – 2007.
4. Карелин А. П., Миронов С. М. Развитие сенсорных технологий и
техники
мониторинга
взрывоопасности
углеводородо-воздушных
смесей.
5. Зятьков И. И., Максимов А. И., Мошников В. А. Сенсоры на основе
полевых транзисторов. – СПб. : ЛЭТИ, 2002.
6. ГОСТ IEC 60079-29-2—2013
7. Артемова А. И., Прокофьев И. В., Суханов А. В. Цифровой газовый
датчик
для
портативных
и
беспроводных
газоанализаторов
//Инженерный вестник Дона. – 2018. – №. 2 (49).
8. Яшин Я.И., Яшин Е.Я., Яшин А.Я. Газовая хроматография. М., 2009.
528 с.
9. Дивин А. Г., Пономарев С. В., Мозгова Г. В. Методы и средства
измерений, испытаний и контроля //Учебное пособие/АГ Дивин, СВ
Пономарев. – 2011. – №. 5.
10. Васильев Б. И., Маннун У. М. ИК лидары дифференциального
поглощения для экологического мониторинга окружающей среды
//Квантовая электроника. – 2006. – Т. 36. – №. 9. – С. 801-820.
85
11. Коваленко А. А., Елисеев А. А. Спектроскопия комбинационного
рассеяния //Методическая разработка Московского государственного
университета имени МВ Ломоносова факультета наук о материалах. –
2011. – С. 37.
12. Петькина Е. Д., Марченко Л. С. Развитие методов определения
компонентного состава природного газа //Научный альманах. – 2016. –
№. 12-2. – С. 225-227.
13. Ященко И. Г. Попутный нефтяной газ Западной Сибири.
14. Ситтиг М. Процессы окисления углеводородного сырья: Пер. с англ. –
Химия, 1970.
15. Ванисов А. М., Клопов А. Л. Особенности газогеохимических съемок в
Западной Сибири (в комплексе дистанционных и геофизических
методов
локального
нефтепрогнозирования)
//Вестник
недропользователя. – 2005. – №. 16.
16. Справочник по геохимии нефти и газа / Ред. С.Г. Неручев. Спб.: Недра,
1986. 200с.
17. Жильцова А. А., Исаев В. И., Коржов Ю. В. Миграция углеводородных
флюидов и геохимический метод индикации залежей //Региональные
проблемы. – 2010. – Т. 13. – №. 1.
18. Зуев В. А. К вопросу о прямых геохимических поисках нефти и газа
//Известия Томского политехнического университета. Инжиниринг
георесурсов. – 2002. – Т. 305. – №. 8.
19. Астахов А. Анализ физико-химических свойств природного газа
//Аналитика. – 2013. – Т. 8. – №. 1. – С. 40-45.
20. Чупин В. В., Жильцов И. Н. Сравнительный обзор современных
средств измерений компонентного состава природного газа //Газовая
промышленность. – 2011. – №. 4. – С. 13-16.
21. Гордадзе Г. Н., Кошелев В. Н., Гируц М. В. Углеводороды нефти и их
анализ методом газовой хроматографии. – 2010.
86
22. Бажанов Ю. В. и др. Количественный анализ газовых сред методом
спектроскопии комбинационного рассеяния света //Аналитика и
контроль. 1998.№ 3/4. – 1998. – С. 65-74.
23. Gordon R. J. et al. 8.25 Hydrocarbon Analyzers //Process Measurement and
Analysis. – 2003. – С. 1358.
24. Боровиков
В.
Н.
Системный
базис
парагенезиса
геополей
продуктивного и надпродуктивного комплексов //Геология нефти и
газа. – 2005. – №. 6.
25. Войтович С. Е. и др. Развитие геофизических и геохимических методов
исследования для геологоразведочных работ поиска и разведки
месторождений нефти и газа //Георесурсы. – 2013. – №. 1 (51).
26. Прохоров А. М. "Справочник по лазерам" М.: Сов. радио. - – 1978. том 1.
27. Ортенберг Ф. С., Трифонов Ю. М. Озон: взгляд из космоса //М.:
Знание. – 1990. – Т. 64.
28. Козинцев В. И. и др. Расчет яркости фона и ослабления лазерного
излучения в ультрафиолетовой области спектра. – 2011.
29. Николаев И. В. и др. Измерение концентрации озона в атмосфере по
поглощению излучения УФ-светодиода //Краткие сообщения по
физике Физического института им. ПН Лебедева Российской Академии
Наук. – 2013. – №. 2.
30. Molina L. T., Molina M. J. Absolute absorption cross sections of ozone in
the 185‐to 350‐nm wavelength range //Journal of Geophysical Research:
Atmospheres. – 1986. – Т. 91. – №. D13. – С. 14501-14508.
31. Козинцев В. И. и др. Основы импульсной лазерной локации. – 2010.
32. Привалов В. Е., Фотиади А. Э., Шеманин В. Г. Лазеры и экологический
мониторинг атмосферы //Учебное пособие. Санкт-Петербург: Лань. –
2013. – Т. 288.
33. Межерис Р. М. Лазерное дистанционное зондирование. – Мир, 1987.
87
34. Malicet J. et al. Ozone UV spectroscopy. II. Absorption cross-sections and
temperature dependence //Journal of atmospheric chemistry. – 1995. – Т. 21.
– №. 3. – С. 263-273.
35. Elterman L. Rayleigh and extinction coefficients to 50 km for the region
0.27 μ to 0.55 μ //Applied Optics. – 1964. – Т. 3. – №. 10. – С. 1139-1147.
36. Городничев В. А. и др. Оценка дальности обнаружения систем
лазерной локации УФ спектрального диапазона //Машиностроение и
компьютерные технологии. – 2014. – №. 11.
37. Тимофеев Ю. М., Васильев А. В. Теоретические основы атмосферной
оптики. – СПб : Наука, 2003. – С. 474.
38. RCA Corporation. RCA electro-optics handbook. – Rca, 1974.
39. Privalov V. E., Shemanin V. G. Lidars for control and measurement
//International Workshop on New Approaches to High-Tech Materials:
Nondestructive Testing and Computer Simulations in Materials Science and
Engineering. – International Society for Optics and Photonics, 1998. – Т.
3345. – С. 6-10.
40. Свердлов Л. М., Ковнер М. А., Крайнов Е. П. Колебательные спектры
многоатомных молекул. – Наука. Гл. ред. физ.-мат. лит., 1970.
41. Bucholtz A. Rayleigh-scattering calculations for the terrestrial atmosphere
//Applied Optics. – 1995. – Т. 34. – №. 15. – С. 2765-2773.
42. Козинцев В. И. и др. Оптико-электронные системы экологического
мониторинга природной среды. – федеральное государственное
бюджетное
образовательное
учреждение
высшего
образования
Московский государственный технический университет имени НЭ
Баумана (национальный исследовательский университет), 2002.
43. Пейсахсон И. В. Оптика спектральных приборов. – Машиностроение,
1975.
44. Павлычева Н. К., Ахметгалеева Р. Р., Сибгатулина Д. Ш. Оптические
схемы
малогабаритных
спектроанализаторов
88
для
мониторинга
гидротехнических сооружений //Вестник Казанского государственного
технического университета им. АН Туполева. – 2019. – Т. 75. – №. 3. –
С. 121-127.
45. Алимов
С.
В.
и
др.
Авиационный
рамановский
лидар
с
ультраспектральным разрешением //Оптический журнал. – 2009. – Т.
76. – №. 4. – С. 41-51.
46. Зарубин И. А., Лабусов В. А., Бабин И. А. Характеристики
малогабаритных спектрометров с дифракционными решетками разных
типов //Заводская лаборатория. Диагностика материалов. – 2019. – Т.
85. – №. 1 (II). – С. 117-121.
47. ang M. et al. General study of asymmetrical crossed Czerny–Turner spectrometer //Applied optics. – 2015. – Т. 54. – №. 33. – С. 9966-9975.
48. Akhmetgaleeva R. R., Muslimov É. R., Pavlycheva N. K. Designing compact S-shaped optical layouts of spectrographs //Journal of Optical Technology. – 2016. – Т. 83. – №. 8. – С. 475-480.
49. Белов М. Л., Городничев В. А., Пашенина О. Е. Сравнительный анализ
мощности входных сигналов лазерных систем локации и видения
ультрафиолетового диапазона //Машиностроение и компьютерные
технологии. – 2013. – №. 08.
50. Comsol
ресурс].
[Электронный
//Model
URL:
https://www.comsol.com/model (дата обращения 15.04.2020)
51. Hamamatsu
ресурс].
[Электронный
//Products
https://www.hamamatsu.com/jp/en/product/index.html
(дата
URL:
обращения
1.05.2020)
52. Findlight
ресурс].
[Электронный
https://www.findlight.net/lasers/solid-state-lasers
3.05.2020)
89
//Lasers
(дата
URL:
обращения
ПРИЛОЖЕНИЕ 1
Волновые числа комбинационных сдвигов исследуемых молекул [40]
Молекула
Колебательный
переход
Волновое
число Длина
перехода νкол, см–1
волны Длина
сдвига λ266
волны
сдвига λ785
N2
0–1
2330,7
283,58
960,79
O2
0–1
1556
277,48
894,22
ν1
2914
288,35
1017,83
ν2
1533,6
277,31
892,44
ν3
3017
289,21
1028,61
ν1
2953,8
288,68
1021,97
ν2
994.8
273,23
851,49
ν12
1190
274,70
865,89
ν3
2875
288,02
1013,80
ν5
1460
276,75
886,61
ν7
1145
274,35
862,53
ν12
2976
288,87
1024,29
ν30
1480
276,90
888,19
ν15
1260
275,22
871,17
ν2
2892
288,17
1015,55
ν3, ν4, ν27
2866
287,95
1012,88
ν9
1346
275,88
877,74
ν10, ν33
1300
275,52
874,21
ν47
1387
276,18
880,91
ν3
2576
285,56
983,98
ν1
3000
289,07
1026,81
ν11
1057
273,69
856,03
ν1
3050
289,48
1032,11
ν9
1097
273,99
858,97
ν10
978
273,10
850,28
CH4
C2H6
C3H8
C4H10
C5H12
C6H14
CH3SH
C2H5SH
90
ПРИЛОЖЕНИЕ 2
Рисунок 1 - Среднее высотное распределение содержание озона
на различныхширотах:1 — 9°; 2 — 30°; 3 — 43°; 4 — 59° [29]
91
ПРИЛОЖЕНИЕ 3
Таблица 1 – Состав сухого воздуха в приземном слое атмосферы [27]
Газовый компонент
Объемное содержание, %
Молекулярная масса
N2
78.084± 0.004
28,016
O2
20.948±0.002
32,000
Ar
0.934±0.001
39,944
CO2
0.033±0.003
44,010
Ne
(1.821±0.004)*10-3
20,183
He
(5.239±0.05)*10-4
4,0030
Kr
(1.14±0.1)*10-4
83,800
Xe
(8.7±0.1)*10-6
131,30
H2
5*10-5
2,0160
CH4
2*10-4
16,000
N 2O
5*10-5
44,013
Таблица 2 - Значения показателя молекулярного рассеяния и полной
оптической толщи земной атмосферы для спектрального диапазона 0,2 ...0,3
мкм [29]
𝜆, мкм
𝜎𝑅 (𝜆), км-1
𝜏𝑅 (𝜆)
-1
0,20
9,202-10
7,826
-1
0,21
7,225-10
6,145
-1
0,22
5,781-10
4,917
-1
0,23
4,691-10
3,990
-1
0,24
3,859-10
3,282
-1
0,25
3,208-10
2,728
-1
0,26
2,690-10
2,228
-1
0,27
2,277-10
1,937
-1
0,28
1,940-10
1,650
-1
0,29
1,665-10
1,416
-1
0,30
1,437-10
1,222
92
ПРИЛОЖЕНИЕ 4
Рисунок 1 - Значения коэффициентов ослабления, рассчитанные для
излучения, распространяющегося горизонтально на уровне моря в
стандартной прозрачной атмосфере: 1 – суммарный коэффициент
ослабления, 2 – рассеяние аэрозолями, 3 – релеевское рассеяние, 4 –
поглощение озоном[26]
93
Отзывы:
Авторизуйтесь, чтобы оставить отзывПроныра, озорник, Любитель книг, Ловкач, игрок, Жизнь между строк. И потому Открыт ему Незримый путь В любую суть. Танец злобного гения На страницах произведения Это игра без сомнения Обречённый ждёт поражения. Подсыпать в душу яд Всегда он рад Всего за час Прочтёт он вас. Он волен взять И поменять Строку и с ней Смысл темы всей.Танец злобного гения На страницах произведения Это игра без сомнения Обречённый ждёт поражения. Открыт роман Читатель пьян Разлив вино - Шагнул в окно. Танец злобного гения На страницах произведения Это игра без сомнения Обречённый ждёт поражения. Танец злобного гения На страницах произведения Это игра без сомнения Обречённый ждёт поражения.
Танец Злобного Гения КиШ
А теперь я скину тексты своих любимых песен для этого:)
и хорошего настроения
удачи
успехов в конкурсе
Наверное было затрачено много времени и труда на работу
Продолжай свое исследование
Админам респект
И продвижения статьи в топы?
Как на счет взаимных комментариев под работами?)
Красиво написанная работа
Так держать
Молодец
Интересная работа!