Главная » 2014 » Сентябрь » 27 » Скачать Континентальный рифтогенез и метаморфическая зональность как следствие термических процессов в литосфере. Полянский, Олег бесплатно
0:34 AM
Скачать Континентальный рифтогенез и метаморфическая зональность как следствие термических процессов в литосфере. Полянский, Олег бесплатно
Континентальный рифтогенез и метаморфическая зональность как следствие термических процессов в литосфере
Диссертация
Автор: Полянский, Олег Петрович
Название: Континентальный рифтогенез и метаморфическая зональность как следствие термических процессов в литосфере
Справка: Полянский, Олег Петрович. Континентальный рифтогенез и метаморфическая зональность как следствие термических процессов в литосфере : диссертация доктора геолого-минералогических наук : 25.00.03 Новосибирск, 2005 320 c. : 71 05-4/71
Объем: 320 стр.
Информация: Новосибирск, 2005
Содержание:
Глава ДИНАМИЧЕСКИЕ КВАЗИТРЕХМЕРНЫЕ МОДЕЛИ ДЕФОРМИРОВАНИЯ ЛИТОСФРЫ
11 Современные геодинамические модели и подходы к описанию деформирования литосферы
12 Управляющие уравнения и ограничения
13 Реология континентальной литосферы
14 Геодинамические причины рас1фытия Байкальского рифта
15 Модель развития осадочного бассейна типа пул-апарт
16 Модель колизионно-сдвиговых деформаций литосферы Алтае-Саянской складчатой области
Глава РАСТЯЖЕНИЕ КОНТИНЕНТАЛЬНОЙ ЛИТОСФЕРЫ И ФОРМИРОВАНИЕ РИФТОВЫХ БАССЕЙНОВ
21 Причины термомеханической дестабилизации континентальной литосферы
22 Реконструкция погружения путем процедуры декомпакции осадков и литостатической разгрузки
23 Модель растяжения
24 Пространственные и временные вариации погружения и седиментации в Кузнецком осадочном бассейне
25 Эволюция Енисей-Хатангского прогиба как пример рифтогенного осадочного бассейна
26 Реконструкция эволюции осадочных бассейнов восточного обрамления Сибирской платформы
Глава IIРОЦЕССЫ КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В ОСАДОЧНЫХ БАССЕЙНАХ
31 Теплоперенос при эволюции Байкальской рифтовой зоны
32 Конвекция флюида в коллекторах осадочного бассейна при тепловом воздействии даек и силлов
33 Модель фильтрации при засолении коллекторов
34 Моделирование теплового воздействия траппов на фазовое состояние и свойства углеводородного флюида
Глава НАРУШЕНИЕ ТЕРМОМЕХАНИЧЕСКОГО РАВНОВЕСИЯ В ЗЕМНОЙ КОРЕ ПРИ ДИАПИРИЗМЕ И ИНТРУЗИЯХ МАГМ
41 Диапиризм в коре и мантии при плотностной инверсии
42 Геодинамические модели гравитационно-неустойчивых систем
43 Математическая постановка задачи
44 Результаты численного моделирования диапиризма
45 Геологические приложения модели гравитационной неустойчивости
46 Моделирование метаморфической зональности с целью определения теплового источника в земной коре
461 Метаморфизм низких давлений/высоких температур в Тонгулакском комплексе (Горный Алтай)
462 Моделирование тепловой истории метаморфической зональности Коннемара (западная Ирландия)
Введение:
Объектом данного исследования являются рифтогенные осадочные бассейны и метаморфические зональные комплексы земной коры, представляющие структуры взаимодействия тектонических и термических процессов в континентальной литосфере.В работе рассматриваются следующие процессы: растяжение континентальной литосферы - рифтогенез и осадконакопление - диапиризм в коре и мантии - интрузии магмы. Выбор геологических структур определялся масштабом процесса (в порядке уменьшения): литосфера в целом, земная кора, осадочный чехол. Исследовались структуры, в основном относящиеся к континентальному Типу литосферы (Кузнецкий, Енисей-Хатангский осадочные бассейны). Рассмотрены примеры бассейнов типа пул-апарт (Тункинский, Чуйско-Курайский прогибы) и др. Основное внимание уделено регионам Центральной и Южной Сибири, входящим в состав Центрально-Азиатского геосегмента. Эта область является крупнейшим на Земле коллажем континентальных литосферных плит, и здесь проявляется большинство из вышеназванных геодинамических обстановок как современных, так и древних. Кроме того, рассматриваются также структуры, территориально не относящиеся к Сибирскому региону: рифтовые структуры Северного моря, Датский и Уэльский бассейны, Днепрово-Донецкий авлакоген.На рис. В.1 схематично показаны участки литосферы, которые в разное время испытали значительные деформации и тепловое воздействие, связанные с нарушением термомеханического равновесия в литосфере. Оконтуренными полями изображено положение ряда осадочных бассейнов Евразии, которые рассматривались в работе. Вторым источником геологической информации являются участки литосферы, испытавшие термические воздействия магматических интрузий и диапиров: траппы Сибирской платформы, метаморфические зональные комплексы Тонгулак (Алтай) и Коннемара (Ирландия), Баян-Хонгорский гранито-гнейсовый диапир (Монголия) (показаны ромбами).Нарушение термомеханического равновесия в континентальной литосфере может привести к формированию термических аномальных областей или к деструкции литосферы. С этой точки зрения и рифтогенез, и развитие метаморфической зональности являются результатом такого нарушения. Рифтогенез, как один из примеров дестабилизации литосферы, хорошо понятен для океанической литосферы, но менее изучен на континентах. Прогресс в понимании океанического этапа развития рифтов определяется более простым строением литосферы океанов по сравнению с континентами, более молодым возрастом современных океанических рифтов, отсутствием наложенных процессов. В отношении метаморфизма пород вблизи магматических интрузий слабо изучено влияние фазового перехода пар/жидкость на эффективность передачи тепла и флюидодинамический режим. Для изучения метаморфизма ранее применялись термические модели с постоянными теплофизическими коэффициентами, не всегда учитывался механизм конвективного теплопереноса норовым флюидом, вместо чего использовалось понятие «эффективная теплопроводность» пород. Как правило, модели ограничивались рассмотрением магматических тел простой формы, не учитывалась нерегулярность контактов с вмещающими породами, неоднородные свойства последних. В силу вышесказанного, актуальность темы представленной работы определяется необходимостью применения многомерных методов численного моделирования к изучению процессов нарушения термомеханического равновесия литосферы. Причем модели должны основываться на появившихся новых данных по глубокому бурению, учитывать возросший объем и новые методы интерпретации сейсмических данных, методы сиквенс-стратиграфии и т.д.Цель исследований - разработка нового интеграционного подхода к решению проблемы нарушения термомеханического равновесного состояния континентальной литосферы путем синтеза и совместной интерпретации математических многомерных сеточных моделей и геолого-геофизических данных (седиментационно-стратиграфический бассейновый анализ, оценки термодинамических параметров метаморфизма, наблюдений теплового потока, GPSизмерения деформаций литосферы и др.).Научные задачи исследований. Цель исследования достигается путем решения четырех взаимосвязанных задач: 1. Установить роль межплитных сил растяжения и мантийного диапиризма при формировании контине1ггальных рифтов; выяснить характер взаимодействия частей лрггосферы с хрупко-пластичной и вязкой реологией в зависимости от напряженного состояния в условиях растяжения, сжатия, сдвига.2. Определить влияние интенсивности растяжения континентальной литосферы на скорость осадконакопления и амплитуду погружения, на динамику и степень метаморфизма осадочных отложений.3. Количественно оценить масштаб и механизмы тепломассопереноса в осадочных бассейнах и вблизи интрузивных тел в зависимости от распределения теплофизических параметров, температуры теплового магматического источника, скорости конвективных течений порового флюида и его фазового состояния.4. Определить природу, форму и глубину магматического теплового источника по термодинамическим характеристикам метаморфических зональных комплексов низкого давления и высоких температур.Этапами решения сформулированных задач являлись: 1. Классификация и систематизация типов изученных осадочных бассейнов для понимания геодинамических причин их формирования.2. Усовершенствование и применение «back;stripping»-мeтoдa для реконструкции эволюции осадочных бассейнов, оценки максимальных температур при метаморфизме погружения, оценки длительности процессов рифтогенеза, а также для оценок толщины термической литосферы и величины растяжения/утонения коры и литосферной мантии в конкретных осадочных бассейнах.Определение количествишых характеристик взаимосвязи процессов: а) растяжения/утонения литосферы, б) осадконакопления и в) метаморфизма погружения.3. Разработка новых и адаптация существующих численных методов для решения геотектонических задач деформирования и напряженного состояния литосферы в различных геодинамических обстановках с учетом латеральнонеоднородной двуслойной литосферы, характеризующейся бимодальной нелинейной реологией; реализация математических моделей формирования континентальных рифтов на основе численных методов конечных элементов в квазитрехмерном приближении (модели в «глобальном» масштабе п* 100-1000 км); количественная оценка влияния контролирующих реологических параметров на характер и степень деформирования литосферы.4. Определение преобладающего типа тепломассопереноса (конвективный или кондуктивный) в континентальных рифтах и контактовых ореолах магматических интрузий; определение вклада внешних (тепловой поток) и внутренних источников тепла (базитовый вулканизм, интрузии в форме силлов и даек, конвективные течения флюидов при уплотнении осадков); учет фазового перехода «газ-жидкость» во флюиде в двумерных нестационарных моделях конвективного тепломассопереноса в пористой среде.Фактический материал и методы исследований.Теоретической основой решения проблемы термомеханического состояния литосферы являются принципы глубинной геодинамики - научного направления, сформировавшегося в течение последних десятилетий (Теркотг, Шуберт, 1985; Артюшков, 1979; 1993; Добрецов, Кирдяшкин, 1994; Хаин, 2001; Леонов и др., 2004; Лобковский и др., 2004).В основу диссертации положены фактические данные, собранные автором по комплексу геодинамических, метаморфических, геофизических, литолого-стратиграфических и структурно-геологических аспектов формирования осадочных бассейнов и метаморфических зональных комплексов в обрамлении Сибирской платформы, а также Алтая, Монголии, Байкальской и ЦентральноАтлантической рифтовых зон. Для сравнительного анализа, верификации и тестирования моделей использованы термодинамические данные по метаморфической зональности Коннемары (Ирландия) (Yardley, 1989, 1986), метаморфизму погружения в Уэльском бассейне (Robinson et all., 1987; 1999), литологостратиграфические данные по строению Днепрово-Донецкого авлакогена (Чекунов и др., 1990), Усть-Енисейского и Хатангского рифтов (Калинко, 1959; Енисей-Хатангская нефтегазоносная, 1974; Палеогеография, 1967), ЛеноВилюйской синеклизы и Верхоянского прогиба (Парфенов, Кузьмин, 2001). Для бассейнового анализа использованы опубликованные данные изучения керна и скважинные исследования (Конторович и др., 1986), измерения теплового потока в осадочных бассейнах Сибирской платформы, Алтае-Саянской складчатой области и Байкальской рифтовой зоны (Дучков и др., 1987; Лысак, Зорин, 1976; Голубев, 1979). Геологические задачи решены с использованием современных методов: палереконструкции осадконакопления и погружения фундамента бассейна, бекстриппинг-анализа, картирования зональных метаморфических комплексов, минералогической термобарометрии для оценок термодинамических параметров метаморфизма, углепетрографического метода определения степени диагенеза и катагенеза органического вещества по отражательной способности витринита и др.В основе исследований лежигг комплексный подход, сочетающий геолого-геофизические наблюдения и математическое моделирование с применением сеточных алгоритмов решения уравнений механики сплошной среды и уравнений нестационарного тепломассопереноса в пористой среде в двумерной и квазитрехмерной постановке. При математическом моделировании использованы конечноразностные и конечно-элементные методы: в задаче диапиризма применен переход к лагранжево-эйлеровой системе координат, который позволяет в отличие от методов сквозного счета определять границу раздела сред в явном виде (Полянский, 1989). Использованы разностные схемы для уравнений движения и теплопроводности, решаемые трех-точечной продольно-поперечной прогонкой и схема чередующихся направлений Писмена-Ракфорда, основанная на методе Яненко расщепления шага по времени. Для изз^ения деформирования литосферы применен подход П.Берда (Bird, 1989) к решению трехмерной задачи о статической деформации тонкой пластины в терминах «мгновенных» скоростей. Визуализация результатов моделирования осуществлена в картографической проекции распределения компонент скорости деформации, эффективных напряжений и др. Для расчета тепломассопереноса в бассейне с изменяющейся геометрией и накоплением осадков использована аппроксимация конечными треугольными элементами с линейными интерполяционными функциями (Kinzelbach, 1986). Применена программа, реализующая описанный алгоритм, на основе метода конечных элементов с автоматическим контролем временного шага, организацией и перестройкой конечно-элементной сетки и присваиванием узловых граничных условий (Bitzer, 1996; Polyansky, Poort, 1998). Для задачи о тепловом влиянии траппового магматизма численная реализация системы уравнений осуществлялась на основе конечно-разностного подхода (Faust, Mercer, 1979). Вычислительный алгоритм строится на основе неявной конечноразностной схемы с итерациями по методу Ньютона-Рафсона (НауЬа, Ingebritsen, 1997). Для решения линеаризованного матричного уравнения используется алгоритм последовательной верхней релаксации, встроенный в итерационный цикл.Автором разработаны собственные компьютерные программы (Полянский, 1987; 1988; Полянский, Волков 1990; Полянский 1998), модифицированы для задач бассейнового анализа программы BASTA (Friedinger, 1988; Фридингер и др., 1991) и BASIN (Bitzer, 1996; 1999; Polyansky, Poort, 2000; Poort, Polyansky 2002), a также использованы имеющиеся готовые некоммерческие пакеты Hydrotherm (Hayba, Ingebritsen 1997), PLATES (Bird, 1989). Последние были адаптированы для применения на современных персональных компьютерах.Как обязательный элемент верификации моделей проводился сравнительный анализ результатов моделирования и природных наблюдений параметров, полученных геологическими методами (палереконструкции осадконакопления, геотермобарометрии, углепетрографическим методом и др.).Защищаемые положения.1. Формирование континентальных рифтов вызвано напряжениями растяжения, возникающими в результате взаимодействия плит. Деформирование литосферы осуществляется в бимодальном реологическом режиме: в хрупкопластичном — преобладают смещения по разломам, в вязком - в виде течений континентальных литосферных масс с элементами вращательного и поступательного движений. При дальнейшем растяжении и утонении литосферы (в 5-6 раз) начинается плавление мантии с последующим вулканизмом и магматизмом в пределах рифтовых зон и переход к акгивной стадии рифтогенеза. Трапповый магматизм совпадает с этапом активизации растяжения, быстрого погружения сформированных бассейнов и может инициировать или завершать фазу рифтогенеза.3. Преобладающий механизм теплопереноса в литосфере определяется пороговой проницаемостью пород (порядка 5*10"'^ м )^ и переходит от конвективного в осадках бассейнов (в условиях метаморфизма погружения) к кондуктивному в консолидированной земной коре (в условиях метаморфизма средних и высоких ступеней). Прогрев пород проницаемого осадочного чехла траппами происходит путем конвекции порового флюида в двухфазном, либо в гомогенном состоянии в зависимости от глубины внедрения магмы. Скорость прогрева вмещающих пород при конвекщ1И флюида вблизи интрузий на один-два порядка выще, чем при кондуктивном механизме.Научная новизна. Личный вклад. В результате проведенных исследований получены следующие новые научные результаты.1. С использованием одномерной модели МакКензи (1978), описывающей процесс растяжения литосферы и образования осадочного бассейна разработан алгоритм реконструкции осадконакопления с учетом тектонического погружения и нагрузки осадочных отложений. Модель растяжения и утонения двуслойной литосферы учитывает уплотнение осадков и температурную релаксацию литосферы. Метод впервые применен для изучения эволюции ряда осадочных бассейнов Сибирской платформы и ее обрамления и позволил классифицировать их как структуры прогибания вследствие растяжения земной коры и мантии с погружением осадков до глубины 10 км и более.2. На основе совокупности данных палеореконструкций осадконакопления, геохронологических данных о возрасте магматизма, геофизических данных о глубинном строении сделан вывод о механизмах формирования ЕнисейХатангского и Кузнецкого бассейнов. В координатах «возраст осадков - глубина» форма эволюционных кривых погружения для этих бассейнов с большой долей вероятности соответствует рифтовому механизму формирования.5. В приближении тонкой пластины, трехмерная задача деформирования термально-неоднородной литосферы переменной мощности с сетью нерегулярных разломов сведена к квазитрехмерному моделированию. Разработана численная модель формирования сдвиговых осадочных бассейнов типа пул-апарт (pull-apart), контролируемая геометрией разломов, величиной скорости сдвига и реологией литосферы. Из сравнения модели с геологическими структурами следует, что впадины типа пул-апарт являются специфическим видом осадочных бассейнов, формирующихся в сдвиговой зоне под действием сверхлитостатических напряжений порядка 25-100 МПа.6. Адаптированы и модифицированы программы гидрогеологических расчетов Hydrotherm и Basin для моделирования теплового режима земной коры с формирующимся осадочным бассейном. Выполнены расчеты двухфазных течений водного порового флюида с учетом фазового перехода пар/жидкость. Сделан вывод о режимах газо-жидкостной конвекции в условиях вскипания и без вскипания: в случае гомогенного флюида (без вскипания) формируются регулярные конвективные ячейки, периодичные как в отношении скорости течеПИЙ, так и по температуре; в случае вскипания вблизи поверхности вместо ячеистой структуры течений адвекция тепла происходит путем восходящих потоков газовой или газожидкостной смеси в виде «гриба» или «плюма».7, На основе теории консолидации Терцаги, разработана полуэмпирическая модель уплотнения пористой среды. По результатам моделирования доказано, что увеличение теплоотдачи в краевых зонах Байкальского рифта происходит без участия мантийных источников тепла. Региональная циркуляция флюида в проницаемом осадочном чехле и по зонам трещиноватости фундамента может перераспределить тепло и создать локальные гидротермальные аномалии в бортах рифта с тешюпотоком до 300-400 мВт/м , что наблюдается в действительности.Теоретическое и практическое значение.Научные разработки автора могут использоваться при решении многих геолого-геофизических задач: геодинамические и кинематические построения, проверка гипотез тектоники плит, изучение температурных условий в осадочных бассейнах, реконструкция их термальной истории и прогноз нефтегазоносности, интерпретация результатов наблюдений над современными региональными деформациями и тепловыми потоками.Полученные результаты по термической истории осадочных бассейнов с наличием траппового магматизма могут применяться в связи с поисками новых, «нетрадиционных» геологических структур, содержащих месторождения углеводородов. Усовершенствованный метод палеореконструкции эволюции бассейнов может быть применен для оценки температур и давлений, соответствующим условиям «нефтяного окна», благоприятным созреванию углеводородов.В представленной работе средствами математического моделирования доказано, что процессы конвективного тепломассопереноса в нефтегазоносных осадочных бассейнах с траппами играют важную роль в миграции и преобразовании состава углеводородного флюида. Исследования показали, что магматическое тепло является причиной фазовых переходов жидких углеводородов в газовую фазу в объемах пород, сопоставимых с телами магматических интрузий (силлов и даек), что может быть использовано при прогнозной оценке углеводородного потенциала осадочных бассейнов с трапповым магматизмом. Для таких бассейнов важное прикладное значение имеют результаты моделирования фильтрации флюида для выявления детальных особенностей гидродинамического режима месторождений углеводородов. Используя разработанный аппарат, можно ответить на вопрос о сохранности углеводородов вблизи магматических тел: в открытой системе нефтяная залежь может быть разрушена, в гидравлически закрытой системе происходит сепарация, т.е. разделение легких и тяжелых фракций; нефтяная залежь при этом обогащается легкими углеводородами. Смоделированные механизмы доказывают возможность образования литологических ловушек углеводородов в соленосно-карбонатных толщах вследствие конвективного переноса и переотложения солей под воздействием тепла интрузий, что может применяться для построения моделей конкретных месторождений.Ряд разработок автора использовался в ходе совместных работ с Институтом геологии нефти и газа СО РАН. Публикации и апробация работы.Проведенные исследования вошли в состав приоритетного направления СО РАН №26 «Геодинамичекая и геохимическая эволюция литосферы и мантии Земли: тектоника, магматизм, флюидный режим и металлогения» по программе № 26.1 «Глубинная геодинамика на основе геолого-геофизических, петрологогеохимических данных и теплофизического моделирования» (акад. Добрецов Н.Л.). Работа проводилась согласно планам НИР Института минералогии и петрографии СО РАН в рамках научного проекта на 2004-2006 гг. № 26.1.6. «Фундаментальные физико-химические проблемы метаморфизма горных пород: экспериментальные и термодинамические исследования минеральных равновесий, кинетика реакций; моделирование метаморфизма как результата геодинамических событий в земной коре» (акад. Ревердатто В.В.). Исследования были поддержаны РФФИ: проекты № 96-05-66051, 99-05-64677, 00-05-65386, 01-05-65253 (руководитель), 02-05-64042, 02-05-65346, 03-05-64014, 04-05-64347 (руководитель), фондами Royal Society 1995-1997, INTAS №134 "Continental Rift Tectonics and Sedimentary Basins Evolution" (1996-1998), IGCP-420 "Continental growth in phanerozoic: evidence from central Asia", Международным научным фондом RPUOOO, RPU300, Президентским фондом поддержки ведущих научных школ № НШ-225.2003.5, интеграционными проектами СО РАН № 30,106.Подходы и результаты, полученные в работе, докладывались автором на 10 всесоюзных и российских и 9 международных конференциях. Среди конференций можно выделить следующие: «Физико-химические и геофизические проблемы образования и эволюции вещества верхней мантии Земли» (Звенигород, 1986); «Composition and processes of deep-seated zones of continental lithosphere» (Novosibu-sk, 1988); INTAS Workshop "Continental Rift Tectonics and Evolution of Sedimentary Basins", (Novosibirsk, 1996); EUG 9 (Strassbourg, 1997); EGS XXII General Assembly (Vienna, 1997); 8th Workshop of the ILP Task Force "Origin of sedimentary basins" (Palermo, 1997); 5* Int. Conf "Computer Aided Design of Advanced Materials and Technologies" (Tomsk, 1997); Int. Conf INTAS Project №134 "Active Tectonic Continental Basins" (Gent, 1998); «Физико-химические проблемы эндогенных геологических процессов» (Москва, 1999); Geofluids III (Barcelona, 2000); XIV Российское совещание по экспериментальной минералогии (Черноголовка, 2001); «Теория нафтидогенеза и органическая геохимия на рубеже веков» (Новосибирск, 2002); «Физика нефтяного пласта» (Новосибирск, 2002); «Геология, геохимия и геофизика на рубеже XX и XXI веков» (Иркутск, 2002); Geofluids IV (Utrecht, 2003); XXXVI Тектоническое совещание (Москва, 2003); «Напряженно-деформированное состояние и сейсмичность литосферы» (Иркутск, 2003); «Глубинные флюиды и геодинамика» (Москва, 2003); XXXVII Тектоническое совещание (Новосибрфск, 2004); 32nd International Geological Congress (Florence, 2004).Фактический материал и основные выводы диссертации изложены в 60 публикациях в отечественных и зарубежных периодических изданиях, а также в отчетах по проектам РФФИ и интеграционных проектов СО РАН. Cxpyiciypa и объем работы. Работа состоит из введения, четырех глав, заключения и списка литературы из 404 наименований. Полный объем диссертации - 320 страниц, 85 рисунков, 16 таблиц.Благодарности. Автор глубоко признателен своему учителю и наставнику академику В.В. Ревердатто, который научил творческому подходу к исследованию, пониманию физических явлений. Автор выражает благодарность академику Соболеву Н.В. создавшему благоприятные условия для творческой работы в Институте минералогии и петрографии ОШТМ СО РАН, а также заведующему лабораторией метаморфизма и метасоматоза д.г.-м.н. Шведенкову Г.Ю. Автор также глубоко признателен за полезные обсуждения результатов в ходе работы и конструктивную крипид' академикам Добрецову Н.Л. и Гольдину СВ., чл.-корр. Артюшкову Е.В., Берниковскому В.А., Склярову Е.В., д.г.-м.н.Лиханову И.И., Кирдяшкину А.Г., Владимирову А.Г., Суворову В.Д., Дучкову А.Д., Шарапову В.Н., Чикову Б.М., Голубеву В.А., Новикову И.С, Ан Ину, П. Берду, Дж. Поорту, К. Битцеру, к.г.-м.н. Соболеву СВ., Хлестову В.В., Бульбаку Т.А., д.ф.-м.н. Макарову П.В., к.ф.-м.н. Стефанову Ю.П., Перепечко Ю.В., Свердловой В.Г., и многим другим коллегам по работе.
