На правах рукописи
ГАЙНАНОВ ВАЛЕРИЙ ГАРИФЬЯНОВИЧ
РАЗРАБОТКА КОМПЬЮТЕРИЗОВАННОЙ ТЕХНОЛОГИИ ОДНОКАНАЛЬНЫХ И МНОГОКАНАЛЬНЫХ СЕЙСМОАКУСТИЧЕСКИХ ИССЛЕДОВАНИЙ НА АКВАТОРИЯХ
25.00.10 Геофизика, геофизические методы поисков
полезных ископаемых
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
доктора технических наук
Москва - 2009
Работа выполнена на кафедре сейсмометрии и геоакустики геологического факультета Московского государственного университета
имени М.В. Ломоносова
Официальные оппоненты:
| доктор технических наук Гогоненков Георгий Николаевич доктор технических наук, ст. научный сотрудник Кондратьев Игорь Константинович доктор технических наук Гуленко Владимир Иванович |
Ведущая организация: | ФГУ НПП Севморгео |
Защита состоится 18 марта 2009 г. в 14 час. 30 мин. на заседании диссертационного совета Д 501.001.64 при Московском государственном университете имени М.В.Ломоносова по адресу: 119991 ГСП-1, Москва, Ленинские горы, ГЗ МГУ, зона А, геологический факультет, аудитория 308
С диссертацией можно ознакомиться в библиотеке геологического факультета МГУ (ГЗ МГУ), зона А, 6-ой этаж
Автореферат разослан У Ф февраля 2009 года.
Ученый секретарь
диссертационного совета Никулин Б.А.
Общая характеристика работы
Актуальность темы. Сокращение запасов полезных ископаемых, в особенности нефти и газа, на суше ведет к расширению их разведки и добычи на акваториях морей и океанов. Для этого сооружаются крупные, в то же время аварийно опасные объекты: буровые платформы, подводные газопроводы, портовые сооружения. Эти проекты требуют серьезного инженерно-геологического обоснования, и не последнюю роль в этом играет метод сейсмоакустического профилирования, который позволяет изучать геологический разрез на требуемую для этих целей глубину с достаточно высокой детальностью. Метод широко используется также при региональных геологических исследованиях и геокартировании.
Повышение эффективности сейсмоакустического профилирования означает увеличение глубинности и разрешающей способности исследований, получение дополнительной информации по физико-механическим свойствам и литологии отложений, достоверного выделения зон повышенной опасности, таких как газонасыщенные отложения, грунты низкой несущей способности и т.д. Однако довольно долгое время первоначальная технология работ - одноканальная система наблюдений с получением временного разреза непосредственно в процессе профилирования - оставалась неизменной, что тормозило дальнейшее развитие метода.
Успехи электроники и цифровой техники обеспечили возможность цифровой регистрации и обработки данных сейсмоакустического профилирования. Это повысило качество и надежность получаемой информации, создало условия для реализации более сложных систем наблюдений.
Так как цифровая обработка и многоканальные наблюдения в сейсморазведке начали применяться гораздо раньше, чем в сейсмоакустическом профилировании, то казалось, что эти технологии могут быть перенесены и на сейсмоакустическое профилирование. Однако практика показала, что механический перенос способов цифровой обработки данных, разработанных для сейсморазведки, и переход к многоканальным наблюдениям, в сейсмоакустическом профилировании очень часто не дают желаемого эффекта. Причины неудач могли быть как в специфике технологии сейсмоакустического профилирования, так и в особенностях строения изучаемой этим методом верхней части разреза. Поэтому потребовались теоретические и экспериментальные исследования с целью установления границ применимости традиционных способов цифровой обработки к сейсмоакустическим данным, их усовершенствования и разработки новых способов. Необходим был также анализ и расчет систем многоканальных сейсмоакустических наблюдений и оценка их реальных возможностей.
Целью работы является теоретическое обоснование, разработка и совершенствование технических средств, технологии полевых исследований и способов обработки данных сейсмоакустического профилирования на акваториях, направленные на повышение разрешающей способности и глубинности исследований, надежности и качества получаемых данных.
Задачи исследований
- Разработка и совершенствование компьютеризованных цифровых регистрирующих комплексов для качественной записи данных сейсмоакустического профилирования в широком частотном и динамическом диапазонах, обеспечивающих реализацию как одноканальных, так и многоканальных систем наблюдений, а также одновременную работу с источиками разных типов.
- Теоретическое обоснование и практическая реализация новых методических приемов наблюдений при сейсмоакустических исследованиях на акваториях, повышающих информативность и надежность получаемых данных.
- Разработка алгоритмов и программ для обработки сейсмоакустических данных с целью повышения разрешающей способности и глубинности исследований, для повышения информативности исследований путем использования кинематических и динамических параметров сигнала.
- Внедрение разработанных технических средств, методических приемов и способов обработки в практику научно-исследовательских и производственных работ при решении инженерно-геологических, региональных и других задач.
Научная новизна
- На основании теоретических и экспериментальных исследований, а также обобщения опубликованной информации, сформулированы требования к компьютеризованным аппаратурным комплексам для одноканального и многоканального сейсмоакустического профилирования.
- Созданы технические средства и программное обеспечение для одноканального и многоканального сейсмоакустического профилирования на акваториях с одновременным использованием двух типов источников - спаркера и пьезоэлектрической антенны или спаркера и бумера.
- Впервые проведены теоретические и экспериментальные исследования многоканальных систем наблюдений для сейсмоакустического профилирования на мелководье с использованием различных типов источников, проанализированы возможности и ограничения таких систем для подавления основных типов помех и определения скорости волн в среде.
- Исследованы возможности ряда известных способов цифровой обработки применительно к сейсмоакустическим данным, разработаны новые оригинальные способы обработки, повышающие эффективность исследований при проведении работ в сложных сейсмогеологических и погодных условиях.
- Разработана технология оценки литологических и физико-механических характеристик отложений с использованием кинематических и динамических параметров сейсмоакустических записей, выявлены наиболее информативные для этих целей параметры.
- Доказана на практических примерах эффективность разработанной технологии проведения полевых работ и обработки данных сейсмоакустического профилирования для повышения разрешающей способности и глубинности исследований, для оценки физико-механических характеристик отложений, литологии, газонасыщенности.
Основные защищаемые положения
- Компьютеризованные регистрирующие системы для одноканального и многоканального сейсмоакустического профилирования, в том числе с одновременным использованием двух типов источников, существенно повышают эффективность исследований на акваториях, сокращают время и стоймость работ.
- Технология многоканальных наблюдений увеличивает глубинность сейсмоакустического профилирования на мелководных акваториях и повышает точность построения разрезов, позволяет определять сейсмические скорости и оценивать физико-механические характеристики отложений.
- Способ обработки данных, позволяющий избавиться от специфических помех, возникающих при сейсмоакустическом профилировании в условиях волнения моря, обеспечивает получение высококачественных временных разрезов даже при неблагоприятных погодных условиях.
- Усовершенствованный способ подавления многократных волн по Бакусу позволяет повысить степень подавления многократных волн и устойчивость работы процедуры в применении к сейсмоакустическим материалам, полученным на мелководье.
- Технология обработки с комплексным использованием кинематических и динамических параметров сейсмоакустических данных позволяет оценивать литологические и физико-механические характеристики, газонасыщенность отложений.
Практическая значимость
- Разработанные и изготовленные под руководством и при непосредственном участии автора компьютеризованные аппаратурные комплексы для одноканального и многоканального сейсмоакустического профилирования используются в ряде научно-исследовательских и производственных организаций в нашей стране (ИО РАН, ВСЕГЕИ, МГУ, ТОИ ДВО РАН, ГЕОХИ РАН) и за рубежом (Ханойский океанографический институт, Геофизическая экспедиция Министерства природных ресурсов - Вьетнам, Стамбульский университет - Турция, Фундаментпроект - Украина).
- Способы обработки данных реализованы в программных комплексах для обработки сейсмоакустических данных, которые используются в перечисленных выше организациях.
- Расчетные соотношения и результаты анализа систем наблюдений для многоканального сейсмоакустического профилирования позволяют заранее оценивать их эффективность для решения тех или иных задач, правильно выбирать параметры таких систем и способы обработки данных.
- Разработанные аппаратурные комплексы, технологии и системы обработки сейсмоакустических данных опробованы при решении задач региональной геологии, геокартирования, инженерно-геологических изысканий на акваториях с различными сейсмогеологическими условиями, и показали высокую эффективность.
- Автор постоянно внедряет свои теоретические и прикладные разработки в учебный процесс - он участник и руководитель учебных практик по морской сейсморазведке и сейсмоакустике на геологическом факультете МГУ, в том числе международных практик УTraining Through ResearchФ под эгидой ЮНЕСКО, им подготовлены учебные пособия и программы учебных курсов, читаемых на геологическом факультете МГУ, в университете Дубна, в Высшей Школе Инновационного Бизнеса МГУ.
ичный вклад автора. Работа является обобщением исследований, выполнявшихся начиная с 70-х годов прошлого века на кафедре сейсмометрии и геоакустики геологического факультета МГУ, а в отдельные годы - также при участии ИО РАН, ГИН РАН, ВСЕГЕИ, учебно-научного центра ЮНЕСКО при МГУ, географического факультета МГУ.
Все результаты получены автором лично, либо под его руководством и при непосредственном участии во всех этапах проектирования и проведения исследований, обработки и интерпретации результатов.
Разработка, изготовление и внедрение компьютеризованных комплексов для одноканального и многоканального сейсмоакустического профилирования осуществлялось совместно с А.С. Зверевым, разработка технологии многоканального сейсмоакустического профилирования и методов обработки данных - совместно с Н.А. Кузубом и М.Ю. Токаревым.
Полевые испытания разработанных аппаратурных комплексов и технологии работ проводились при поддержке и непосредственном участии Л.Р.Мерклина и О.В. Левченко (ИО РАН), М.А. Спиридонова и Ю.П. Кропачева (ВСЕГЕИ), М.К. Иванова (Центр ЮНЕСКО по морской геологии и геофизике при МГУ), Л.В. Поляка (Центр полярных исследования Университета штата Огайо, США), М.Л. Владова (кафедра сейсмометрии и геоакустики геологического факультета МГУ), руководства и сотрудников компании ДЕКО проект.
Апробация и публикации. Основные положения и результаты работы докладывались на ведомственных, всероссийских и международных конференциях: UNESCO-IOC-ESF 4-th Post Cruise Meeting УSedimentary Basins of the Mediterranean and Black SeasФ (Moscow,1996), International Earth Sciences Colloquium of the Aegean regions (IESCA-95, Izmir, 1996), 11-th Petroleum Congress and Exibition of Turkey (Ankara, 1996), Congress УGas and Fluids in Marine SedimentsФ (Amsterdam, 1997), Congress УCarbonate Mud Mounds and Cold Water ReefsФ (Gent, Belgium, 1998), Third Workshop on Land Ocean Interactions in the Russian Arctic (LOIRA) (Moscow, 2000), 3rd Balkan Geological Congress and Exibition (Sofia, Bulgaria, 2002), Ломоносовских чтениях. Секция геологии. (Москва 2003, 2005), Международной геофизической конференции Москва 2003, VII Международной научно-практической конференции Геомодель-2005 (Геленджик, 2005), Международной конференции Инженерная геофизика-2006 (Геленджик, 2006), Международной конференции Санкт-Петербург-2006, IV Всероссийском литологическом совещании ОСАДОЧНЫЕ ПРОЦЕССЫ: седиментогенез, литогенез, рудогенез (эволюция, типизация, диагностика, моделирование) ( Москва, 2006), Международной конференции Нефть и газ Арктического шельфа-2006 (Мурманск, 2006), Международной научно-практической конференции Инженерная и рудная геофизика-2007 (Геленджик, 2007), Конференции по инженерной геологии (Москва, 2007), Международной конференции Санкт-Петербург-2008, Международной научно-практической конференции Инженерная и рудная геофизика-2008 (Геленджик, 2008),
По результатам выполненных исследований автором опубликовано более 70 работ в виде статей и тезисов докладов, в том числе 25 статей в рецензируемых журналах. Прилагаемый в конце автореферата список основных трудов по теме диссертации содержит 50 наименований.
Структура и объем работы. Диссертация состоит из введения, 5 глав и заключения, а также списка литературы из __ наименований. Работа содержит __ страниц, в том числе __ рисунка.
Благодарности. Проведение исследований, сопряженных с аппаратурными разработками и полевыми экспедициями, не было возможно без поддержки и непосредственного участия многих людей.
Автор благодарен М.К. Иванову (Центр ЮНЕСКО по морской геологии и геофизике), Л.Р. Мерклину и О.В. Левченко (ИО РАН), М.А. Спиридонову и Ю.П. Кропачеву (ВСЕГЕИ), А.Б. Матвеенко и К.Г. Пухаеву (ДЕКО проект), Л. В. Поляку (Университет штата Огайо), оказавшим существенную поддержку в проведении исследований.
Автор искренне благодарен А.С. Звереву, без конструкторского таланта которого аппаратурные разработки не могли быть реализованы.
Автор благодарен В.В. Калинину, В.К. Хмелевскому, В.А. Богословскому, Е.А. Ефимовой, за ценные советы, внимание и поддержку при подготовке работы.
Автор глубо признателен своим коллегам из кафедры сейсмометрии и геоакустики геологического факультета МГУ: М.Л. Владову, А.В. Старовойтову, В.А. Стручкову, М.Ю. Токареву, Н.В. Шалаевой, а также Н.А. Девдариани и В.Н. Ефремову за многолетнее плодотворное сотрудничество.
Содержание работы
Введение содержит общую характеристику работы, включая решаемые задачи, научную новизну, защищаемые положения и практическую ценность.
Глава 1. Сейсмоакустическое профилирование на акваториях Ц история и современное состояние.
В этой главе дан анализ истории развития и современного состояния метода сейсмоакустического профилирования на акваториях, обоснованы цели и задачи диссертационной работы.
В первом параграфе вкратце изложена история развития метода, который появился в конце 50-х годов XX века. Специалистами Ламонтской геологической обсерватории (Колумбийский университет, США) и Массачусетсткого технологического института (США) для проведения морских сейсмических исследований было предложено использовать новую комбинацию технических средств - факсимильный регистратор эхолотного типа, источник многократного действия электроискрового (спаркер) или газовзрывного типа и буксируемую пъезосейсмографную косу [Beckmann et al.,1959, Ewing et al., 1961, McGuinness et al., 1962, Edgerton, 1963, Edgerton et al., 1964]. Это сочетание оказалось чрезвычайно удачным, и метод начал быстро развиваться в мире, в том числе и в России [Калинин и др., 1966 - 1983].
Теперь уже очевидно, что первоначально предложенная комбинация - электроискровой источник и буксируемая пъезосейсмографная коса - оказалась наиболее оптимальной для такого рода исследований. Ведь первоначально многие разработчики и исследователи сосредоточили свои усилия на совершенствовании традиционного эхолота для применения в геологических целях [Clay et al., 1962, 1964, Nawar, 1964, Дубров, 1967], где для возбуждения и приема упругих колебаний в воде используется гидроакустическая антенна - пъезоэлектрический или магнитострикционный излучатель-приемник, возбуждается многофазный узкополосный сигнал. Это направление также получило развитие - сейчас существуют акустические профилографы самого разного принципа действия [www.geoacoustics.co.uk, www.benthos.com, www.innomar.com], однако из-за высокой частоты колебаний (5-10 кГц) они пригодны в основном для исследования лишь рыхлых осадков, глубинность исследований редко превышает первые метры.
Изначально сейсмоакустическое профилирование развивалось как одноканальный метод как по причине заимствования регистрирующих систем от эхолотов, так и простоты и мобильности аппаратуры и технологии работ. Да в те годы из-за неразвитости еще технологии воспроизводимой записи и обработки данных и не было необходимости в многоканальных наблюдениях.
Несмотря на многочисленные попытки усоврешенствования источников упругих волн и приемных систем для сейсмоакустического профилирования, они конструктивно до настоящего времени практически не изменились. Решающие изменения произошли в регистрирующих системах - прогресс в области цифровой техники сделал возможным регистрацию сейсмоакустических данных практически в любом частотном и динамическом диапазонах без предварительного сжатия. Это внесло кардинальные изменения в конструкции регистрирующих систем в целом - теперь в них отсутствуют специализированные дорогостоящие аппаратные средства, практически нет механических узлов, наиболее подверженных поломкам. Основная часть регистратора теперь - персональный компьютер, куда вводятся оцифрованные данные и регистрируются на штатных носителях информации, а для контроля визуализируются на экране компьютера. Переход к полностью цифровой регистрации данных и увеличение производительности персональных компьютеров позволили по-новому подойти и к цифровой обработке сейсмоакустических данных - если раньше цифровой обработке подвергались лишь отдельные части данных, то сейчас возможна достаточно глубокая цифровая обработка всех данных.
Прогресс в развитии регистрирующих и обрабатывающих систем создал условия для реализации более сложных систем наблюдений при сейсмоакустических исследованиях на акваториях - теперь можно проводить непрерывное многоканальное сейсмоакустическое профилирование [Гайнанов и др., 2006-2008], а также одновременные наблюдения с использованием источников и приемных систем, работающих в разных частотных диапазонах [Безродных и др. 2002, Гайнанов и др. 2008,].
В заключение этого обзора хочу назвать имена исследователей - сотрудников лаборатории сейсмоакустики геологического факультета МГУ, внесших ощутимый вклад в развитие сейсмоакустического профилирования на акваториях, и которые также явились моими учителями: Ш.А. Азими, А.В. Калинин, В.В. Калинин, Б.Л. Пивоваров, И.Я. Ковальская, М.Х. Фаталиев.
Во втором параграфе дан краткий обзор современных технических средств для сейсмоакустического профилирования - источников возбуждения упругих волн и приемно-регистрирующих систем.
Довольно долгое время научно-исследовательские и производственные организации нашей страны для сейсмоакустического профилирования использовали аппаратурный комплекс Аквамарин, разработанный еще в советское время. Хотя регистрирующая часть комплекса давно уже морально устарела и не используется, высоковольтные источники СКАТ с накапливаемой энергией 10 кДж еще недавно можно было встретить в таких организациях, как ИО РАН и МАГЭ [Казанин и др., 2002].
Из современных систем известен многоканальный геофизический комплекс серии САК (АО "Моринжгеология", г. Рига), позволяющий осуществлять сейсмоакустическое профилирование одновременно по двум разночастотным каналам [Безродных и др., 2002]. В институте Океангеология и НПО Севморгео есть своя аппаратура для сейсмоакустического профилирования. В ООО ДЕКО проект разработаны компактный и мощный (1 кДж) электроискровой источник и 16-канальная регистрирующая система.
Производственные организации, проводящие инженерно-геофизические работы, часто по требованию заказчика должны использовать аппаратуру известных зарубежных компаний, которые для сейсмоакустического профилирования на акваториях предлагают:
- Компания GeoAcoustics Ltd (UK) [www.geoacoustics.co.uk] имеет GeoPulse Boomer/Sparker, в состав которого входят: высоковольтный источник GeoPulse 5420S с выходным напряжением 3750 В и энергией 105-455 Дж, бумер на катамаране или многоэлектродный спаркер, гидрофон, усилитель и регистратор GeoPulse Sonar Processing system.
- Applied Acoustic Engineering Ltd (UK) [www.appliedacoustics.com] предлагает широкую линейку высоковольтных источников от 50 Дж до 6000 Дж, для питания бумеров (AA200, AA300 Boomer Plates) и спаркеров разной мощности (Squid 500 Sparker, Squid 2000 Sparker, 6000J Delta Sparker). Выходное напряжение 2500-4000 В.
- Оригинальные электроискровые источники для очень высоко разрешающих исследований предлагает Geo-Resources Consulting B.V. (The Netherlands) [www.geo-resources.com]. Например: Geo-Spark 1000DT (1000 Дж, Deep Towed), как утверждают авторы, может обеспечивать разрешающую способность не хуже 10 см по вертикали. Компания изготовляет приемные косы до 48 каналов, одноканальные и многоканальные регистрирующие системы для высокоразрешающих наблюдений.
Западные компании изготовляют высоковольтные накопители электрической энергии теперь только с полупроводниковыми выключателями, применяют комбинированные системы безопасности, излучатели изготовляют из специальных сплавов. Это сильно повышает стоимость аппаратуры - по сравнению с зарубежными аналогами отечественные системы, имея сопоставимые функциональные возможности, стоят в несколько раз дешевле.
В конце 1990-х - начале 2000-х годов при активном участии автора в компании Geont Shelf разработаны компьютеризованные сейсмоакустические комплексы для различных видов сейсмоакустического профилирования на акваториях и скважинных сейсмоакустических работ [Гайнанов и Зверев, 2008], которые используются в ряде научно-исследовательских организаций у нас в стране (ИО РАН, ВСЕГЕИ, МГУ, ТОИ ДВО РАН, ГЕОХИ РАН) и за рубежом (Вьетнам, Турция, Украина). Принципы построения этих комплексов изложены в главе 2.
В третьем параграфе дан краткий обзор программных средств, применяемых для обработки данных сейсмоакустического профилирования, которые разделены на 2 группы:
1. Программные комплексы, предназначенные для обработки данных сейсморазведки МОВ-ОГТ типа ProMax, GeoFocus, RadExPro [RadExPro Plus 3.6., 2005]. Это - мощные системы, рассчитанные на обработку больших массивов данных многоканальной сейсморазведки, содержащие в себе практически все мыслимые процедуры обработки сейсмических данных. В принципе, они в состоянии обеспечить обработку сейсмоакустических данных по любым современным алгоритмам. Однако именно их достоинства превращаются в недостатки: они громоздки, очень часто только специалист высокой квалификации в состоянии разобраться в тонкостях реализации тех или иных процедур, наконец, по стоимости они недоступны для научно-исследовательских или инженерно-геологических организаций, которые в основном и используют сейсмоакустическое профилирование. Исключение - программа RadExPro, которая более приспособлена для таких работ, и по стоимости вполне доступна, хотя и у нее есть некоторые недостатки.
2. За рубежом распространены специальные программы для обработки данных сейсмоакустического профилирования и акустических профилографов, например, Coda Octopus [Octopus 361 Seismic Processing Toolkit], DrGeo [www.activesoft.net.au], SES "ISE 2.9" [www.innomar.com]. Эти программы компактны, просты в использовании и стоят недорого. Однако в них предусмотрена лишь достаточно простая обработка данных, как полосовая частотная фильтрация, регулировка амплитуд и т.п. В некоторых программах предусмотрена ручная пикировка горизонтов в целях интерпретации.
Опробование нами этих программ показало, что отнюдь не все идеи по обработке сейсмоакустических данных в их рамках могут быть реализованы. Поэтому автором создана своя система обработки данных сейсмоакустического профилирования [Гайнанов, 2008], которая обладает рядом достоинств. Принципы построения этой системы обработки рассматриваются в главе 3.
В четвертом параграфе рассматриваются достижения и проблемы внедрения цифровой технологии в сейсмоакустическое профилирование.
Хотя современная сейсморазведка немыслима без цифровой обработки данных, в сейсмоакустическом профилировании, несмотря на наличие цифровой записи данных, отдача от цифровой обработки не такая заметная, да и не всегда она применяется. Часто применение к сейсмоакустическим данным способов цифровой обработки, разработанных для сейсморазведки, не дает желаемого эффекта. Чтобы определить причины такого отставания и найти пути повышения эффективности цифровой обработки сейсмоакустических данных, нужно сделать:
- Оценить параметры регистрируемой при сейсмоакустическом профилировании информации, допустимые искажения при записи, и выработать требования к неискаженной регистрации данных.
- Проанализировать как принципиальные, так и реальные возможности известных способов обработки с учетом специфических особенностей сейсмоакустического профилирования, а при необходимости разработать новые способы обработки таких данных.
- Оценить возможности многоканальных и двухчастотных сейсмоакустических наблюдений для повышения разрешающей способности и глубинности исследований, для определения литологии и физико-механических свойств отложений.
Автор постоянно участвовал в решении указанных проблем, разрабатывая и совершенствуя цифровую регистрирующую аппаратуру для сейсмоакустического профилирования, создавая программное обеспечение для регистрирующих систем и обработки данных, проводя испытания и внедрение этих систем. Наиболее значимые достигнутые результаты представлены в последующих главах данной работы.
Глава 2. Разработка и совершенствование приемно-регистрирующих систем для сейсмоакустического профилирования на акваториях.
Эта глава посвящена научному обоснованию и разработке компьютеризованных систем для сейсмоакустического профилирования.
В первом параграфе сформулированы требования к приемно-регистрирующим системам для сейсмоакустического профилирования.
Частотный диапазон сигналов в сейсмоакустическом профилировании зависит от глубинности и разрешающей способности исследований, от типа применяемых источников. Наиболее низкочастотные источники - пневматические. Согласно нашим данным основная энергия возбуждаемых ими колебания обычно находится в диапазоне частот 20 - 250 Гц. Наиболее высокочастотные источники - бумеры, возбуждают колебания в диапазоне частот от 200 - 8000 Гц. Полоса возбуждаемых спаркером частот лежит внутри этого диапазона (100 - 1500 Гц).
Шумы буксируемой пьезокосы имеют разную природу. Но наиболее интенсивные шумы возникают из-за флуктуации гидростатического давления вокруг пьезоприемников, происходящего вследствие колебания глубины буксировки косы. Эти шумы имеют максимальную интенсивность в диапазоне низких частот - от долей герц до первых десятков герц. Способ борьбы с ними - оптимально подобранные фильтры высоких частот.
Таким образом, диапазон регистрируемых при сейсмоакустическом профилировании частот может быть ограничен полосой 20 - 8000 Гц.
Динамический диапазон принимаемых сигналов при сейсмоакустическом профилировании определяется рядом факторов: типом и энергией источника волн, длиной и количеством элементов в приемной косе, параметрами системы наблюдений и сейсмогеологическими условиями. Поэтому можно лишь примерно оценить максимальный и минимальный динамические диапазоны.
Наши оценки показывают, что полный динамический диапазон принимаемых сигналов может достигать до 5*106 раз или 130 дб [Полшков и др., 1976; Гайнанов, 1977]. Однако в подавляющем большинстве случаев динамический диапазон регистрируемых сигналов гораздо меньше. Например, при используемой обычно технологии сейсмоакустического профилирования наиболее интенсивной принимаемой волной оказывается не прямая, а отраженная от дна волна. Амплитуда ее при малых глубинах воды также может быть значительной, но при этом и уровень шумов высокий - шумы судна отражаются от дна и принимаются косой без ослабления. В результате динамический диапазон сигналов при сейсмоакустическом профилировании обычно не превышает 60 - 90 дб. Кроме того, он может быть сжат до 40 - 50 дб за счет применения линейной регулировки усиления.
Скорость потока информации при сейсмоакустическом профилировании большая из-за сравнительно высокой частоты принимаемых сигналов и частого повторения циклов возбуждения и приема. Например, при 16-канальной регистрации поток информации может быть выше 600 кбайт/с. Хотя для современных компьютеров эта скорость приема данных и вполне доступная, но с учетом того, что сигнал нужно принять, обработать, визуализировать на экране, записать на диск, распечатать на принтере, могут возникать трудности с использованием ресурсов компьютера. Это ставит ограничения на канальность регистрации и частоту циклов возбуждения. Например, в 16-канальной станции Нильма минимальный период регистрации равен 0,9 с.
В 60-е и 70-е годы прошлого века для контроля параметров аппаратуры и принимаемых сигналов требовались специальные осциллографы и генераторы. В современной аппаратуре предусматривают встроенные генераторы, которые подключаются к входу усилителя по командам с компьютера. Программа регистрации визуализирует сигнал в виде временного разреза на экране и на специальном окне в виде осциллограммы, тем самым компьютер объединяет в себе цифровой регистратор, факсимильный регистратор и осциллограф. Он же управляет запуском источников, принимает и регистрирует навигационную информацию.
Во втором параграфе представлены функциональные схемы цифровых регистрирующих трактов сейсмоакустической аппаратуры.
Согласно сформулированным требованиям сейсмоакустический регистрирующий канал должен содержать следующие функциональные узлы: предварительный усилитель с высокоомным дифференциальным входом, фильтр высоких частот (ФВЧ), программный регулятор усиления (ПРУ), основной усилитель, фильтр низких частот (ФНЧ) - антиаляйсинг-фильтр, согласующий каскад, аналого-цифровой преобразователь (АЦП).
В такой схеме, хотя суммарный динамический диапазон и достаточно большой (118 дб), мгновенный динамический диапазон определяется разрядностью АЦП, и в нашем случае (14 разрядов) составляет около 8000 раз (78 дб), что не всегда бывает достаточным. Еще до появления дельта-сигма АЦП с 24 разрядами мы разработали схему многокаскадного усилителя с одновременной оцифровкой сигнала с 3-х последовательно включенных каскадов усиления обычным 14-разрядным АЦП. Программа собирает единый сигнал с учетом коэффициентов усиления каскадов, тем самым обеспечивается динамический диапазон регистрации 8000*40*40=1280000 раз или 120 дб.
В третьем параграфе дано описание управляющей программы для компьютеризованных сейсмоакустических комплексов.
Программа обеспечивает настройку усилителей и фильтров, запуск источников, прием и регистрацию сейсмоакустических данных, навигационной информации, визуализацию данных в виде временных разрезов и как осциллограмм, вывод твердых копий на принтер и ряд вспомогательных функций.
Разработаны программы для одноканальных и многоканальных систем, а также для двухканального профилирования с двумя источниками, возбуждающими упругие импульсы разного частотного состава и с разным периодом возбуждения. Здесь приводится описание последнего варианта.
Запуск источников и регистрации осуществляется от ЭВМ через цифровой блок управления. Минимальный период запуска равен 0,25 с. Период запуска любого источника может увеличиваться кратно этому времени. Источники могут запускаться поочередно, или в любой последовательности.
Рис. 1. Окно программы регистрации с изображением временных разрезов от двух источников и всплывающее окно УОсциллографФ.
Рис. 2. Временные разрезы, полученные с бумером (а) и со спаркером (б).
Цифровое управление делает систему гибкой и удобной в работе. Теперь оператор с одного компьютера может управлять источниками, процессом регистрации и контролировать сигналы по каналам. На экране можно визуализировать временной разрез по первому или второму каналу по одному из источников, или два разреза одновременно один по бумеру, другой по спаркеру. На всплывающем окне можно наблюдать форму сигнала по обоим каналам в выбранном диапазоне времен и амплитуд (рис. 1). На экран выводится также вспомогательная информация: номера трасс, время записи, координаты по GPS, значения задержки и т.д.
Прилагаемая программа обработки позволяет читать данные по определенному каналу, а также отбирать их по признаку источника (рис. 2). Далее обработка данных производится раздельно.
В четвертом параграфе дано краткое описание разработанных под руководством автора компьютеризованных сейсмоакустических комплексов для одноканального и многоканального профилирования.
Исходя из своего многолетнего опыта сейсмоакустических исследований на акваториях и используя современную техническую базу, нам удалось сделать их компактными, надежными и доступными по стоимости, а самое главное, с широкими функциональными возможностями. Аппаратно комплексы могут быть изготовлены в разных вариантах:
Комплекс для региональных сейсмоакустических исследований с глубинностью до 500 - 1000 м имеет в составе, кроме усилительно-регистрирующей части, мощный электроискровой источник 7,5 кДж и косу длиной 60 м. 5-ти электродный электроискровой излучатель возбуждает колебания с относительно низкочастотным спектром (центральная частота около 80 Гц). Период возбуждения импульсов - 3 - 10 с.
Усилитель имеет дифференциальный вход, программную (линейную) регулировку усиления и установку максимального коэффициента усиления от компьютера в пределах 5 - 5000 раз, ФНЧ (500, 1000, 2000, 4000, 10000 Гц), ФВЧ (30, 60, 120, 240, 500, 1000 Гц). Время ПРУ - 1 с. Максимальная амплитуда выходного сигнала - 4 В.
Комплекс для высокоразрешающих исследований со спаркером - наиболее востребованный аппаратурный комплекс для сейсмоакустического профилирования на мелководных акваториях в инженерно-геологических целях. С этим комплексом может достигаться глубинность исследований до первых сотен метров при разрешающей способности около 2 - 5 м.
В составе комплекса - высоковольтный силовой блок 600 Дж, многоэлектродный электроискровой излучатель (30 - 100 электродов) и пьезосейсмографная коса длиной 10 м. Полоса регистрируемых частот 100 - 1500 Гц. Период возбуждения импульсов - 1 с.
Комплекс для высокоразрешающих исследований с бумером предназначен для инженерных исследований на мелководных акваториях, когда требуется разрешающая способность исследований 0,3 - 0,5 м. Глубинность при этом не превышает первых десятков метров.
В состав комплекса входят маломощный высоковольтный накопитель электрической энергии 300 Дж, излучатель типа бумер и небольшая косичка длиной 1 м. Полоса регистрируемых частот 100 - 8000 Гц.
Сейсмоакустический комплекс для двухчастотного профилирования включает в себя оба вида источников - спаркер и бумер - и двухканальную приемно-регистрирующую систему [Гайнанов и др., 2008]. Комплекс позволяет при однократном прохождении профиля получать высокое разрешение для верхней части разреза и достаточно большую глубинность исследований. В отличие от систем, использующих для этого два независимых комплекта аппаратуры, разработанный аппаратно-программный комплекс обеспечивает синхронизацию источников, исключающую их взаимное влияние, позволяет подбирать параметры и регистрировать данные на одном компьютере.
Комплекс для многоканального сейсмоакустического профилирования включает в себя многоканальную косу (16, 24 или 32 канала), многоканальный усилитель, 32 канальный АЦП и источник типа спаркер или бумер.
С учетом специфики сейсмоакустического профилирования (работа с маломерных судов, осуществление спуско-подъемных операций вручную) многоканальная коса изготавливается в максимально облегченном варианте - диаметр полиуретанового шланга всего 25 мм. Чтобы вся начинка могла размещаться в таком тонком шланге, применяются специальные малогабаритные пьезоприемники и блоки предварительных усилителей.
Программа сбора данных позволяет принимать и регистрировать сигналы с 32 входных каналов. Для контроля параметров сигналов в программе предусмотрены дополнительные окна - Одноканальный осциллограф и Многоканальный осциллограф. На основном окне изображается временной разрез по одному из выбранных каналов.
Для обработки многоканальных данных по способу ОГТ используется система УRadExProФ компании ДЕКО геофизика [RadExPro Plus 3.6, 2005].
Скважинный сейсмоакустический комплекс предназначен для проведения высокочастотных сейсмических исследований (ВСП и межскважинное просвечивание) в скважинах глубиной до 200 м, имеет в своем составе высоковольтный силовой блок 600 Дж со скважинным излучателем акустических колебаний и двухканальную приемно-регистрирующую систему со скважинной пьезоэлектрической косой.
Усилители имеют 3 ступени усиления, выбираемые в автоматическом режиме, как это описано выше. Полоса пропускания усилителей 20 - 5000 Гц.
Выводы.
- На основании теоретических и экспериментальных исследований установлены количественные закономерности для частотного и динамического диапазонов регистрируемых при сейсмоакустическом профилировании сигналов, разработаны функциональные и принципиальные схемы цифровых регистрирующих систем.
- Разработано программное обеспечение для одноканальных, многоканальных, а также двухчастотных систем сейсмоакустического профилирования на акваториях.
- Разработаны и изготовлены компьютеризованные сейсмоакустические комплексы для различных видов сейсмоакустического профилирования на акваториях и скважинных сейсмоакустических работ.
Глава 3. Исследование эффективности способов цифровой обработки сейсмоакустических данных и разработка новых способов.
В этой главе рассмотрены возможности и условия применения различных способов цифровой обработки данных сейсмоакустического профилирования, предложены и обоснованы новые способы.
Попытки повышения эффективности сейсмоакустических исследований на акваториях путем внедрения способов цифровой обработки данных, напрямую позаимствованных из сейсморазведки, часто не дают желаемого эффекта. Автором проведены теоретические и экспериментальные исследования с целью установления границ применимости тех или иных способов цифровой обработки данных, их усовершенствования и разработки новых способов.
В первом параграфе сформулированы цели и задачи цифровой обработки данных сейсмоакустического профилирования.
Цифровую обработку данных сейсмоакустического профилирования можно разделить на две части: 1) обработка данных с целью получения качественных временных разрезов; 2) определение кинематических и динамических характеристик волн с целью оценки физико-механических, литологических и других свойств отложений.
Первая часть обработки - это повышение отношения сигнал/помеха с целью увеличения глубинности и разрешающей способности исследований. Для этого применяются разные виды частотной фильтрации, в том числе предсказывающую деконволюцию и деконволюцию по форме сигнала. Для улучшения коррелируемости границ на временных разрезах, полученных в условиях волнения моря, применяют приемы обработки наподобие коррекции статических поправок. Для подавления кратных волн разработаны специальные приемы.
Вторая часть обработки - это определение таких сейсмических характеристик разреза, как пластовые скорости, коэффициенты отражения, коэффициенты поглощения. Так как вычисление их с достаточной точностью не всегда бывает возможно, то часто ограничиваются сравнительной оценкой таких параметров, как амплитуды отраженных волн, видимые периоды.
Обработка данных многоканального сейсмоакустического профилирования имеет свои особенности, связанные с сейсмогеологическими условиями изучаемого разреза, применяемыми системами наблюдений и характеристиками излучаемых сигналов. Технологии многоканального сейсмоакустического профилирования и обработки данных посвящена отдельная глава 4.
Во втором параграфе рассматриваются некоторые теоретические и практические аспекты применения цифровых фильтров при обработке данных сейсмоакустического профилирования.
На полевых примерах показано, например, что в одних случаях предсказывающая деконволюция дает неплохой эффект в повышении разрешающей способности и подавлении многократных волн, а других - может приводить к ложным результатам. Поэтому применять ее надо очень осторожно.
Деконволюция по форме сигнала в целях повышения разрешающей способности в приведенных примерах работает гораздо лучше предсказывающей деконволюции. Однако использовать ее можно только в тех случаях, когда форма сигнала практически не меняется от трассы к трассе, и ее можно вычислить по отраженной от дна волне.
В третьем параграфе описываются разработанные автором способы исключения влияния волнения моря на результаты сейсмоакустического профилирования.
При работах в открытом море сейсмоакустическое профилирование часто приходится проводить в условиях волнения моря, недопустимого с точки зрения методики высокоразрешающих исследований, что приводит к резкому ухудшению качества записей, а иногда делает их вообще неинтерпретируемыми (рис. 3, 4) [Гайнанов и др., 2007, 2008].
Существуют разные способы борьбы с такого рода помехами, но некоторые из них не применимы при сейсмоакустическом профилировании, другие - часто не дают желаемого эффекта.
Например, за рубежом в состав аппаратуры акустических профилографов включают специальные датчики колебательных движений - MRU (Motion Reference Unit), сигналы которых используются для введения поправок на вертикальные смещения антенны профилографа [Product Line SES-2000, www.innomar.com]. При сейсмоакустическом профилировании использование таких датчиков практически невозможно, так как источник и приемник буксируются за бортом судна каждый по отдельности, и совершают независимые колебания, в том числе и вращательные.
В обрабатывающей системе RadExPro есть несколько процедур, предназначенных для сглаживания границ на временных разрезах. Опробование их на полевых материалах показало, что с одной стороны их применение бывает недостаточно удобным, с другой стороны их помехоустойчивость невысокая, и в сложных случаях они перестают работать.
Автором разработаны свои оригинальные способы исключения влияния волнения моря на материалы сейсмоакустического профилирования, дающие положительные результаты иногда даже в кажущихся безнадежными ситуациях. Основные идеи их заключаются в следующем:
Анализ данных, полученных в самых разных условиях волнения моря, показал, что как в точке возбуждения, так и в точке приема, волнение оказывает 3 вида воздействия на сигнал: 1) время прихода отраженных волн меняется в связи с изменением высоты приемника и источника над отражающей границей; 2) амплитуда принимаемых волн меняется в связи с изменением глубины источника и приемника, а также кривизны поверхности воды - сильной отражающей границы; 3) по этим же причинам меняется форма и частотный состав принимаемого сигнала.
Полное исключение влияния волнения на сейсмоакустические данные означает корректировку всех вышеуказанных искажений сигнала. При этом мы исходим их предположения, что сигналы по соседним трассам должны быть коррелированны, т.е. резкие изменения времени вступления, амплитуды и формы сигнала вызваны исключительно влиянием волнения.
Рис. 3. Иллюстрация работы основного алгоритма исключения влияния волнения: a Ц фрагмент полевого временного разреза с изображением прослеженной линии дна (1) и высокочастотной составляющей этой линии (2); б - тот же фрагмент после исключения влияния волнения. Источник - бумер.
Рис. 4. Иллюстрация работы процедуры исключения влияния волнения с предварительным редактированием трасс: a Ц увеличенный фрагмент полевого временного разреза; б - тот же фрагмент после обработки. Источник - спаркер.
Основной алгоритм предусматривает прослеживание первых вступлений отраженной от дна волны по всему обрабатываемому интервалу профиля, высокочастотную фильтрацию прослеженной линии, и сдвиг трасс на оставшиеся после фильтрации значения времен (рис. 3). Однако волнение вызывает еще флуктуации амплитуды и формы сигнала и повышение уровня шума на записи, так что непосредственное прослеживание первых вступлений донной волны по исходной записи часто не удается. Поэтому нами предусмотрены специальные процедуры предварительного редактирования трасс и разработан алгоритм повышенной помехоустойчивости для прослеживания донной волны (рис. 4).
В четвертом параграфе рассматривается применимость различных способов подавления многократных волн при обработке сейсмоакустических данных и предлагается усовершенствованный способ подавления многократных волн по Бакусу, который позволяет повысить степень подавления многократных волн и устойчивость работы процедуры в применении к сейсмоакустическим материалам, полученным на мелководье.
В сейсморазведке методом отраженных волн на акваториях сильной помехой являются многократно отраженные в водном слое волны, и борьбе с ними уделяется много внимания. В многоканальной сейсморазведке сейчас для подавления кратных волн наряду с методом суммирования по способу ОГТ [Mayne, 1962; Мешбей, 1973] используют веерную фильтрацию [Кондратьев, 1972], преобразования Радона [Moore et al., 2002; Weglein, 1999], предсказание и вычитание на основе волновых уравнений (метод SRME - surface related multiple elimination) [Verschuur, 2006]. В одноканальном сейсмоакустическом профилировании эти способы, к сожалению, в принципе не применимы.
Разработано много разных способов одномерной фильтрации данных с целью подавления многократных волн [Backus, 1959; Kunetz & Fourmann, 1968; Middleton & Whittlesey, 1968; Калинин, 1976, 1978; Сильвиа и Робинсон, 1983; Verschuur, 2006]. Однако при применении их к данным сейсмоакустического профилирования часто достаточный уровень подавления кратных волн не достигается, или одновременно с кратными волнами с записи исчезают и однократные отражения. Причины разные - отношение сигнал/помеха на сейсмоакустических данных низкое, принятая модель сейсмограммы не соответствует реальным данным, и т.д.
На основе развития идей, предложенных М.Бакусом [Backus, 1959] автором разработан способ подавления многократных волн на записях сейсмоакустического профилирования, наиболее полно учитывающий влияние мешающих факторов, и обеспечивающий устойчивую настройку процедуры на реальные полевые данные. Алгоритм предусматривает постепенное приближение к конечному результату - наилучшему подавлению кратных волн.
Рис. 5. Подавление кратных волн: а фрагмент временного разреза до обработки (спаркер 600 Дж, Карское море); б тот же фрагмент после исключения влияния волнения и подавления кратных волн. 1- кратные и неполнократные волны.
Сначала наблюденные данные подгоняются под модель, принятую для описания образования кратных волн, затем производится первичное - грубое определение времени пробега донной волны и коэффициента отражения от дна k. После методом итераций вычисляются точные значения и k, осуществляется вычитание кратных волн.
Если требуемый уровень подавления кратных волн не достигается, то модель усложняется, например, коэффициент отражения от дна считается частотно зависимой. Тогда вычитаемая трасса кратных волн пропускается через фильтр с предварительно определенной характеристикой. Однако чем сложнее модель, тем труднее определить его параметры по относительно сильно зашумленным полевым данным. Поэтому на практике чаще лучшие результаты получаются по более простым способам.
Одним из преимуществ разработанного способа является то, что вне зависимости от того, лучше или хуже подавляются кратные волны, однократные отражения на временном разрезе остаются практически нетронутыми. Это помогает опытному интерпретатору отличить кратные волны от однократных, даже если они не полностью подавились в результате работы процедуры.
Пятый параграф посвящен технологии определения динамических параметров волн по данным сейсмоакустического профилирования.
Первое требование для успешного определения динамических параметров волн - вся предварительная обработка данных должна выполняться с сохранением Уистинных амплитудФ, и должны быть получены временные разрезы, где визуально уже можно заметить изменения динамических характеристик волн [Gainanov et al, 1998].
Использование динамических параметров сигнала в интерпретации предусматривает, во-первых, проведение преобразований Гильберта и получение разрезов мгновенных амплитуд, фаз и полярностей [Танер и Шериф, 1982], во-вторых, вычисление сейсмических атрибутов вдоль выделенных границ. В системе RadExPro для этого есть соответствующие процедуры. Однако процедуры для определения полярностей отражений нет. Дело в том, что принцип определения полярности по знаку амплитуды колебания на сейсмической трассе в момент локального максимума на трассе мгновенных амплитуд справедлив лишь для симметричного по форме (нуль-фазового) сигнала, отраженного от кровли достаточно мощного пласта. Отклонения формы сигнала от требуемой, интерференция отражений от тонких пластов могут привести к неправильному определению полярности отражений.
Автором составлена специальная процедура и разработана технология, обеспечивающая устойчивое определение полярности отражений в большинстве встречающихся на практике случаев. На множестве примеров показано, что полярность отражений часто оказывается наиболее надежным параметром в расшифровке природы отражающих границ, например, лярких пятен на сейсмических разрезах (рис. 6) [Gainanov, 1998; Гайнанов, 2008].
Рис. 6. Фрагмент сейсмоакустического временного разреза с ляркими пятнами (а) и разрез полярностей (б).
При интерпретации необходимо совместное использование динамических параметров сейсмоакустической записи (амплитуда, полярность, частотный состав) с кинематическими параметрами (скорость волн, конфигурация границ). Во многих случаях это позволяет найти наиболее достоверные варианты интерпретации материалов, даже когда отсутствуют скважинные данные о вещественном составе отложений [Гайнанов, 2008].
В шестом параграфе дается описание разработанного программного комплекса для обработки данных сейсмоакустического профилирования.
Многолетний научно-технический и производственный опыт сейсмоакустических работ на акваториях, а также знакомство с отечественными и зарубежными программными средствами для обработки сейсморазведочных и сейсмоакустических данных, позволили автору создать свою систему обработки данных сейсмоакустического профилирования [Гайнанов, 2008], которая обладает определенными преимуществами, а именно: 1) она компактна и проста в использовании; 2) работает быстро; 3) реализует, наряду с общеизвестными процедурами обработки сейсмических данных, процедуры, разработанные автором специально для сейсмоакустического профилирования (Глава 3, з 3, з 4).
Интерфейс программы подобен интерфейсу широко распространенных Windows программ и является интуитивно понятным. После открытия файла данных и загрузки заданного количества трасс, на экране появляется изображение временного разреза. Простым щелчком мыши на выбранной точке разреза можно посмотреть осциллографическое изображение данной трассы, производить редактирование трасс, посмотреть частотный спектр трасс.
Специальное меню открывает список имеющихся процедур обработки данных. После выбора процедуры и задания параметров программа, переходит в интерактивный режим подбора параметров, если такой режим предусмотрен, или производит обработку всех трасс и выводит результаты на экран. Обработанные данные могут быть сохранены для дальнейшей обработки, или же обработанный фрагмент временного разреза может быть сохранен в виде изображения для последующей обработки в графических системах.
На рис. 7 приводится пример улучшения качества материала последовательным применением различных обрабатывающих процедур.
Рис. 7. Фрагмент сейсмоакустического временного разреза: а) до обработки; б) после полосовой частотной фильтрации (700 - 2000 Гц); в) после фильтрации, мьютинга, редактирования трасс, исключения влияния волнения.
Выводы.
- Исследованы известные, предложены и практически реализованы новые оригинальные способы обработки данных, повышающие эффективность сейсмоакустического профилирования на акваториях.
- Разработана технология оценки литологических и физико-механических характеристик, газонасыщенности отложений с использованием динамических параметров сейсмоакустических записей.
- Разработан программный комплекс для обработки сейсмоакустических данных, существенно улучшающий качество материалов за счет применения как известных, так и оригинальных способов обработки.
Глава 4. Расчет систем наблюдений и исследование эффективности многоканального сейсмоакустического профилирования на модельных и полевых примерах.
В этой главе анализированы возможности и ограничения многоканального сейсмоакустического профилирования для подавления основных типов помех и определения скорости волн в среде.
В первом параграфе описана технология многоканального сейсмоакустического профилирования на акваториях, разработанная при участии автора в компании ДЕКО Проект [Гайнанов и др., 2006].
Для проведения работ по методу МОГТ используется 16-канальный приемно-регистрирующий комплекс и источник типа бумер с центральной частотой излучения 1Ц2 кГц или спаркер с центральной частотой излучения 200Ц400 Гц. Шаг между приемниками - 2 м, общая длина приемной системы - 30 м. Буксировка излучателя и приемной косы осуществляется на глубине 0,4 м от поверхности воды в случае бумера и на глубине 0,8 м в случае спаркера. Кратность перекрытия составляет в среднем 16 раз при бинировании через 2 м.
Обработка данных проводится в системе RadExPro [RadExPro Plus 3.6, 2005] и включает такие процедуры, как частотная фильтрация, регулировка амплитуд, скоростной анализ, суммирование по ОГТ. По мере необходимости осуществляется также ввод и коррекция статических поправок, деконволюция и двухмерная фильтрация данных [Гайнанов и др., 2007].
Многоканальные наблюдения практически во всех случаях позволяют получить более качественные временные разрезы. Это достигается не только за счет суммирования многоканальных записей, но также благодаря возможности предварительной выбраковки шумных каналов или варьированию количества суммируемых каналов в зависимости от глубины отражающих границ. Во многих случаях удается достаточно точно определять скорости сейсмических волн в осадках и строить глубинные разрезы даже в отсутствии скважин на исследуемых площадках.
Однако результаты не всегда оправдывают ожидания. Например, степень подавления кратных волн суммированием часто оказывается недостаточной для уверенного прослеживания глубоких границ, иногда глубинность исследований не превышает глубинность одноканальных наблюдений. Для того чтобы разобраться в причинах этих успехов и неудач, мы решили провести анализ результатов полевых работ, сопоставляя их с теоретическими расчетами, основы которых даны в работах [Гогоненков и др., 1975; Гольдин, 1974; Гольцман, 1964; Козлов и др., 1973; Мешбей, 1973], и проверяя собственными модельными примерами [Гайнанов, 2007; Гайнанов и Токарев, 2008].
Во втором параграфе рассматриваются свойства многоканальной косы как настраиваемой группы для подавления помех.
Частотная характеристика группы с равномерным распределением приемников для плоских волн выражается формулой [Гольцман, 1964]
,
где - частота; - кажущаяся скорость; - запаздывание волны на последнем приемнике, - длина группы, n - число приемников в группе.
В область пропускания группы попадают те волны, для которыхне превышает четверти периода. Волны, для которых больше периода, группой подавляются. Так как волны, отраженные от границ на разных глубинах имеют разные запаздывания, то фиксированная группа, оптимальная для приема одних волн, не подходит для приема других.
В работе на полевых примерах показано, как при регистрации сигнала с каждого приемника отдельно можно управлять параметрами группы, подбирая для каждого интервала глубин требуемое число суммируемых каналов.
По-настоящему преимущество многоканальной регистрации проявляется при суммировании по ОГТ с вводом кинематических поправок. Лучшее подавление помех достигается веерной фильтрацией, но для исключения аляйсинг-эффекта перед веерной фильтрацией необходимо спрямить оси синфазности помехи специальными поправками [Гайнанов и Токарев, 2008].
В третьем параграфе анализируется эффективность многоканальных сейсмоакустических систем в борьбе с кратными волнами.
Получены рассчетные соотношения для определения оптимальных параметров систем наблюдений по способу ОГТ в зависимости от глубины исследуемых границ, от отношения эффективных скоростей однократных и многократных волн, от частотного состава волн (рис. 8).
На модельных и полевых примерах показана возможность подавления кратных волн до 2,5 - 4 раз. Как ожидалось, степень подавления кратных волн падает при глубинах границ, существенно превышающих длину расстановки. Но она резко падает и на малых глубинах из-за растяжения сигнала на дальних каналах после введения кинематических поправок - здесь линия допустимого растяжения сигнала (2) ограничивает длину расстановки гораздо раньше, чем она достигнет длины (1), необходимой для подавления кратных волн.
Рис. 8. Графики зависимости минимальной длины расстановки от времени вступления отраженных волн для разных соотношений скоростей многократных (V1) и однократных (V2) волн при f=1 кГц (бумер) и f=250 Гц (спаркер) (1); линии относительного растяжения сигнала на 20% (2).
Длина расстановки, необходимая для подавления кратных волн, определяется исходя из условия, чтобы разница времен на последних суммируемых каналах была не меньше преобладающего периода волн
. (1)
иния фиксированного растяжения сигнала S, согласно [Kleyn, 1983] определяется формулой
(2)
Для проверки полученных выводов в условиях, максимально приближенных к реальным, были рассчитаны синтетические многоканальные сейсмограммы для моделей сред, близких по параметрам к сейсмогеологическим разрезам в районах работ. Эти сейсмограммы затем обрабатывались в системе RadExPro по тому же графу, что и полевые сейсмограммы, и оценивалось подавление помех. Результаты опубликованы в работах [Гайнанов, 2007; Гайнанов и Токарев, 2008].
Для большего подавления кратных волн рекомендуется применять совместно с суммированием по ОГТ одноканальные способы подавления кратных волн (рис. 9).
Рис. 9. Пример комплексного подавления кратных волн: а) временной разрез ОГТ; б) суммирование по ОГТ после потрассного подавления кратных волн в водном слое. 1 - однократная волна; 2 - многократные волны.
В четвертом параграфе оценивается эффективность скоростного анализа по многоканальным сейсмоакустическим данным.
По теоретической оценке чувствительность этой процедуры, так же, как и степень подавления кратных волн, определяется частотным составом сигнала и длиной расстановки. Но есть некоторые ньюансы:
Расчеты на модельных сейсмограммах показали низкую разрешающую способность скоростного анализа для малых глубин, что объясняется сильным растяжением сигнала с удалением и ограничением числа суммируемых каналов. В то же время оказалось, что для глубин, сравнимых с длиной расстановки, многоканальная система имеет относительно высокую избирательность по эффективным скоростям при не очень больших способностях к подавлению кратных волн. Это можно объяснить тем, что для разделения волн по скоростям достаточно, чтобы разделились соответствующие максимумы на графиках, а для подавления волн необходимо, чтобы их спектры полностью попали в полосу подавления.
На полевых данных максимумы на спектрах скоростей получаются более расплывчатыми из-за разброса времен вступлений, но в целом скорости определяются достаточно уверенно.
Выводы.
- Многоканальная регистрация данных позволяет эффективно использовать интерференционные свойства многоэлементной сейсмоакустической косы.
- Теоретические и модельные расчеты показывают, что при суммировании по ОГТ возможно ослабление многократных волн до 2,5Ц4 раз при достаточно большой длине косы и высокочастотном спектре сигнала.
- На глубине менее 10Ц20 м для бумера и менее 40Ц70 м для спаркера невозможно эффективное подавление многократных волн суммированием по ОГТ вследствие сильного растяжения сигнала на дальних каналах после ввода кинематических поправок, возможно лишь некоторое их ослабление.
- Скоростной анализ по сейсмограммам ОГТ возможен в более широком интервале глубин. Однако на малых глубинах точность скоростного анализа падает так же по причине сильного растяжения сигнала на дальних каналах.
- Использованная в полевых наблюдениях коса (N=16, =2 м) оптимальна для исследования разреза до глубин 30Ц50 м.
Глава 5. Технологии сейсмоакустических исследований на примерах решения региональных и инженерно-геологических задач.
В этой главе показана эффективность разработанных технологий сейсмоакустических исследований на практических примерах.
Представленные в данной работе аппаратура, методика исследований и способы обработки данных разрабытывались и совершенствовались в процессе самых разных геолого-геофизических исследований. Результаты исследований, доказывающие эффективность этой технологии для повышения разрешающей способности и глубинности сейсмоакустического профилирования на акваториях, для оценки физико-механических характеристик отложений, литологии, газонасыщенности, докладывались на научных конференциях, опубикованы в научной литературе, список которых прилагается. Здесь приводятся некоторые наиболее характерные примеры.
Первый параграф посвящен региональным исследованиям.
Автор участвовал в исследованиях конусов выноса больших рек [Гайнанов и др., 1995], грязевых вулканов и диапиров Черного и Средиземного морей [Лимонов и др., 1992; Gaynanov & Cifci, 1996; Cifci et al, 1997; Gaynanov et al, 1998; Гайнанов, 2008], карбонатных горок в Атлантическом океане [Gaynanov, 1998; Гайнанов, 2008].
Использованная технология сейсмоакустического профилирования:
Вариант 1. Возбуждение колебаний осуществлялось группой из 6 электроискровых источников с энергией 5 кДж (центральная частота спектра 80 Гц), прием - на одноканальную косу длиной 50 м. Усилитель - с линейной регулировкой усиления, с полосой пропускания 20 - 2000 Гц.
Источник буксировался на глубине 5 м на расстоянии 50 м от кормы судна, коса - на такой же глубине на расстоянии 100 м. Профилирование проводилось на скорости судна 6 - 8 узлов. Возбуждение колебаний производилось с интервалом 10 с. Регистрация данных осуществлялась цифровой компьютеризованной системой с 12-разрядным АЦП.
Вариант 2. Пневмопушка объемом 3 л (центральная частота спектра 50 Гц), 6-канальная коса общей длиной 100 м, 6-канальный сейсмический усилитель с полосой пропускания 15 - 250 Гц.
Обработка данных включала частотную фильтрацию, регулировку амплитуд, приведение записей к единой задержке, предсказывающую деконволюция, а также вычисление мгновенных амплитуд и определение полярности отражений (рис. 10).
Рис. 10. Сейсмоакустический временной разрез через массив карбонатных горок на плато Роколл (а), разрез полярности отражений (б). Отражения от границы в основании горок (показаны желтой стрелкой) имеют отрицательную полярность, что свидетельствует о том, что карбонатные постройки растут на неконсолидированных осадочных отложениях.
Во втором параграфе приводятся примеры высокоразрешающих региональных и инженерно-геологических исследований в Карском море и в центральной части Каспийского моря.
Технология. Электроискровой источник 600 Дж, многоэлектродный излучатель (центральная частота спектра 400 Гц). Приемная коса 10 м (10 пьезоприемников). Буксировка источника и приемника на глубине 1 м на скорости 6 узлов. Интервал возбуждения 1,5 с. Регистрация данных в полосе частот 30 - 1000 Гц.
Обработка данных проводилась примерно по тому же графу, что и для региональных исследований. Но в связи с тем, что при высокоразрешающем сейсмоакустическом профилировании качество данных существенно ухудшается при волнении моря, то здесь большой эффект дала процедура исключения влияния волнения моря (Глава 3, з 3), а в некоторых случаях - процедура подавления многократных волн в водном слое (Глава 3, з 4).
Рис. 11. Ледниковые отложения на эродированной поверхности коренных пород перекрыты сверху слоистыми постледниковыми осадками.
Рис. 12. Осадочный бассейн, заполненный осадками с субпараллельной слоистостью, которые облекают неровную поверхность коренных пород и ледниковых отложений.
Данная технология оказалась особенно эффективной при исследовании следов глобальных оледенений в Баренцевом и Карском морях, обеспечивая при благоприятных сейсмогеологических условиях глубинность исследований до сотен метров при разрешающей способности по вертикали 2 - 5 м. Это позволяло просветить весь разрез ледниковых отложений, иногда даже сквозь мощную толщу постледниковых осадков, захватывая также и доледниковые коренные породы (рис. 11, 12).
Региональные и инженерно-геологические исследования под постановку полупогружных плавучих буровых установок (ППБУ) в глубоководной части Каспийского моря были проведены ИО РАН на НИС Рифт в 2004 году при участии сотрудников МГУ. Результаты сейсмоакустических исследований представлялись на всероссийских и международных конференциях [Гайнанов и др., 2005 - 2006; Левченко и др., 2006; Поляков и др., 2006].
Рис. 13. Временной разрез по профилю Махачкала-Ялама (а) и его увеличенный фрагмент (б).
Инженерные изыскания проводились на 2-х площадках - Ялама-Самур и Центральная. Для понимания общей геологической ситуации на прилегающей территории были отработаны несколько региональных профилей, временной разрез по одному из которых представлен на рис.13.
В третьем параграфе приводятся примеры многоканальных и комбинированных сейсмоакустических исследований в инженерных целях.
При выполнении инженерно-геологических изысканий под строительство крупных сооружений перед сейсмоакустическими исследованиями ставятся повышенные требования - необходимо изучать разрез на глубину до 50Ц60 м с разрешающей способностью 0,5 м и лучше, нужно строить не только временные, а и глубинные разрезы, требуется оценить физико-механические свойства отложений.
Здесь, безусловно, преимущества имеют многоканальные наблюдения, позволяющие определять скоростные характеристики отложений, а также двухчастотное профилирование, обеспечивающее высокую разрешающую способность. Параллельно проводится гидролокация бокового обзора (ГЛБО), которая позволяет получить сведения о характере донных отложений не только вдоль профиля, но и по площади.
Технология многоканальных наблюдений описана в главе 4 з 3, а двухчастотного профилирования - в главе 2 з 3.
Инженерно-геологические изыскания в заливе Св. Петра проводились под строительство портовых сооружений. Требовалось расчленить разрез на верхнюю рыхлую толщу и коренные отложения, выделить разломы, зоны залегания рыхлых УполужидкихФ илов и газоносных отложений.
Работы проводились при сильном волнении моря, и получение качественных временных разрезов оказалось возможным именно благодаря многоканальным наблюдениям и цифровой обработке данных. Многоканальность наблюдений позволила определить скорости сейсмических волн в осадках и построить более точные глубинные разрезы. Благодаря относительно высокому отношению сигнал/помеха на суммарных разрезах ОГТ удалось выделить множество амплитудных аномалии типа ляркое пятно и определить полярности отражений. Комплексная интерпретация динамических и кинематических параметров сейсмоакустической записи позволила вполне достоверно идентифицировать разные по литологии и физико-механическим характеристикам отложения при отсутствии скважинных данных (рис. 14).
Рис. 14. Пример кинематической и динамической обработки данных по одному из профилей: a) график средней амплитуды отраженной от дна волны; б) суммарный временной разрез ОГТ (указаны значения пластовых скоростей по результатам скоростного анализа); в) разрез полярности отражений.
При инженерных изысканиях под строительство трубопроводных переходов и мостов через реки использовалось многоканальное сейсмоакустическое профилирование на двух частотных диапазонах и гидролокация бокового обзора [Гайнанов и др., 2008].
На сейсмоакустическом разрезе через русло большой реки (рис. 15) четко выделяется подошва аллювиальных отложений, накопленных у левого берега. На профиле, пройденном вдоль реки, а также на мозаике ГЛБО, наблюдаются песчаные валы, накопленные в русле.
Рис.15. Сейсмоакустический разрез через русло реки (а) и вдоль русла (б). Источник - бумер, суммирование по ОГТ 10 каналов. Песчаные валы хорошо видны также на мозаике ГЛБО (в).
По результатам скоростного анализа известняки, залегающие в глубоководной части русла выделяются достаточно большими значениями пластовых скоростей - 2000 - 2200 м/с, в то время как аллювиальные отложения характеризуются значениями около 1600-1700 м/с. Тем не менее, эти значения очень низкие для известняков, что свидетельствует о их разрушенности и низкой несущей способности.
В четвертом параграфе приведены примеры сейсмоакустических исследований в скважинах с использованием разработанной автором аппаратуры и методики работ.
Электроискровой источник оказался удобным и достаточно эффективным источником возбуждения продольных волн в водонаполненных скважинах. Поэтому техника, разработанная для морских исследований, после некоторой модификации может применяться и для исследований в скважинах [Гайнанов и др., 1984; Владов и др., 1988; Владов, 2003].
Наиболее эффективной оказалась технология межскважинного сейсмоакустического просвечивания для изучения физико-механических характеристик под основанием строящихся и уже сооруженных больших энергетических объектов. Первоначально предусматривалось только томографическое восстановление скоростного разреза между скважинами по временам пробега прямых волн. Однако наличие на некоторых сейсмограммах записей отраженных волн подсказало идею использовать и эти волны для получения дополнительной информации о разрезе [Гайнанов, 2003; Гайнанов и др., 2005, 2006].
Выводы.
- Эффективность разработанных технических средств, методических приемов и программ обработки доказана на многочисленных примерах решения региональных и инженерно-геологических задач.
- Комбинация одноканальных и многоканальных систем наблюдений, а также источников типа бумер или спаркер позволяет исследовать разрез на достаточно большую глубину с высокой разрешающей способностью.
- Цифровая обработка позволяет построить качественные временные и глубинные разрезы, оценить физико-механические и литологические свойства отложений.
- Оригинальная процедура исключения влияния волнения моря во многих случаях позволяет получить качественные временные разрезы в таких погодных условиях, при которых раньше получение качественных данных было невозможно.
Заключение.
Основные результаты диссертационной работы следующие:
- На основании теоретических и экспериментальных исследований, а также обобщения опубликованной информации, сформулированы требования и разработаны компьютеризованные регистрирующие системы для одноканального и многоканального сейсмоакустического профилирования, в том числе с одновременным использованием двух типов источников, существенно повышающие эффективность исследований на акваториях, сокращающие время и стоймость работ.
- Впервые проведены теоретические и экспериментальные исследования многоканальных систем наблюдений для сейсмоакустического профилирования на мелководье с использованием различных типов источников, анализированы возможности и ограничения таких систем для подавления основных типов помех и определения скорости волн в среде.
- Исследованы возможности ряда известных способов цифровой обработки применительно к сейсмоакустическим данным, разработаны новые оригинальные способы обработки, повышающие эффективность исследований при проведении работ в сложных сейсмогеологических и погодных условиях.
- Разработана технология оценки литологических и физико-механических характеристик отложений с использованием кинематических и динамических параметров сейсмоакустических записей, выявлены наиболее информативные для этих целей параметры.
- Показана на большом количестве практических примеров эффективность разработанной технологии проведения полевых работ и обработки данных сейсмоакустического профилирования для повышения разрешающей способности и глубинности исследований, для оценки физико-механических характеристик отложений, литологии, газонасыщенности.
Направления повышения эффективности сейсмоакустических исследований на акваториях в будущем видятся автору в внедрении новых методических приемов, таких как 3-D наблюдения, что в свою очередь, невозможно без разработки компактных и надежных систем контроля глубин и координат источников и приемников. Еще не исчерпаны возможности совершенствования способов обработки, например, перспективной может оказаться совместная динамическая обработка данных двухчастотного профилирования для определения поглощающих свойств разреза.
Список основных публикаций по теме диссертации
Статьи в реферируемых журналах, рекомендуемых ВАК:
- Опыт оптического анализа данных непрерывного сейсмического профилирования на море // Вестн. Моск. ун-та. Сер. 4. Геология. 1976. № 1, С. 116Ц120.
- Опыт определения эффективных скоростей при непрерывном сейсмическом профилировании с двумя приемными системами // Океанология. 1978. т. 18, вып. 3 (соавторы: Ельников И.Н., Стручков В.А.)
- Сейсмические свойства горных пород в районе Крымского полигона по данным скважинных исследований // Вестн. Моск. ун-та. Сер. 4. Геология. 1983. № 3 (соавтор: Рыбачук Г.О.)
- Опыт применения динамических параметров сигналов в НСП для целей геологического картирования на акваториях // Вестн. Моск. ун-та. Сер. 4. Геология. 1989. № 4.
- Исследование возможностей повышения разрешенности сейсмических данных при изучении глубоководных конусов выноса // Вестн. Моск. ун-та. Сер. 4. Геология. 1995. № 4, С. 85Ц87 (соавторы: Коротков И.П., Тихоцкий С.В., Трофимов С.В.)
- Сейсмоакустические исследования следов покровных оледенений в Карском море // Вестн. Моск. ун-та. Сер. 4. Геология. 2005. № 1, C. 38-44. (соавторы: Поляк Л.В., Гатауллин В.Н., Зверев А.С.)
- Использование отраженных волн при межскважинном сейсмическом просвечивании // Вестн. Моск. ун-та. Сер. 4. Геология. 2005. № 5, С. 83Ц85. (соавтор: Рыковская Н.В.)
- Новые данные о сейсмостратиграфии и процессах седиментогенеза на западном склоне Среднего Каспия // Докл. РАН. 2006. т. 411, № 5. С. 663Ц666. (соавторы: Левченко О.В., Мерклин Л.Р., Поляков А.С., Росляков А.Г.)
- Опыт сейсмоакустического профилирования с многократным перекрытием: возможности и ограничения // Разведка и охрана недр.2006, № 12. с. 21-24 (соавторы: Кузуб Н.А., Токарев М.Ю., Клещин С.М.)
- Комплексное использование технологий скважинных сейсмических наблюдений при инженерно-геологических исследованиях // Там же. с. 30-32. (соавтор: Скворцов А.Г.)
- Многоканальное сейсмоакустическое профилирование на разных частотных диапазонах - реальные возможности // Разведка и охрана недр. 2008. № 1, С. 35Ц38. (соавторы: Токарев М.Ю., Зверев А.С., Росляков А.Д.)
- Возможности и ограничения многоканального сейсмоакустического профилирования в инженерных целях: теория и практика. // Вестн. Моск. ун-та. Сер. 4. Геология. 2008. № 5 (соавтор: Токарев М.Ю.)
- Об использовании динамических параметров записи при сейсмоакустическом профилировании. // Вестн. Моск. ун-та. Сер. 4. Геология. 2008. № 6.
- Сейсмоакустический комплекс для двухчастотного профилирования на акваториях // Океанология. 2008. № (соавтор: Зверев А.С.)
- Программный комплекс для обработки данных сейсмоакустического профилирования // Океанология. 2008..№
- Примеры использования динамических параметров сейсмоакустической записи при интерпретации // Разведка и охрана недр. 2008. № 12.
- Сейсмоакустические методы при инженерно-геологических изысканиях на реках // Там же. (соавторы: Старовойтов А.В., Баскакова Г.В.)
- Mud volcanoes and dome-like structures at the Eastern Mediterranean Ridge // Mar. Geophys. Res., 1997, 19, p. 421 - 438. (co-authors: Cifci G., Limonov A., Dimitrov L.)
- Seismic evidence for gas accumulation related to the area of mud volcanism in the deep Black Sea // Geo-Marine Letters, 1998, 18, p. 139-145. (co-authors: Bouriak S., Ivanov M.)
- Kara Sea expedition yields insight into LGM ice sheet extent // Eos. 2002. 83, No. 46, p. 525, 529. (co-authors: Polyak L., Gataullin V. et al.)
- The eastern extent of the Barents-Kara Ice Sheet during the Last Glacial Maximum based on seismic-reflection data from the eastern Kara Sea // Polar Research (APEX Special Volume), 2008. (co-authors: Polyak L., Niessen F., Gataullin V.)
Другие публикации:
- К вопросу о путях развития способов регистрации и обработки в методе сейсмоакустического профилирования. Экспресс-информация ВИЭМС, сер. "Морск. геол.и геофиз.", 1976. (соавторы: Полшков М.К., Калинин А.В., Калинин В.В.)
- Регистрирующая система для сейсмоакуст. проф. на базе микро-ЭВМ "Электроника 60". Тез. докл. всесоюзн. школы "Тех. средства и методы исслед. Мирового океана". М.,1987.
- Экспериментальные исследования зависимости динамических параметров сейсмического сигнала от литологического состава донных осадков. В сб. "Геол. конт. террасы окраинных и внутр. морей". Изд. МГУ, 1989. (соавторы: Старовойтов А.В., Девдариани Н.А., Рыбалко А.Е.)
- Сейсмоакустические исследования следов глобальных оледенений в Баренцевом и Карском морях. (Тезисы) - Ломонос. чтения. Секц. геол. М. 2003. (соавтор: Зверев А.С.)
- Высокоразрешающее межскважинное сейсмическое просвечивание с использованием отраженных волн. (Расшир. тезисы) - Междун. геофиз. конферен. Москва 2003. М. 2003.
- Высокоразрешающие сейсмоакустические исследования в центральной части Каспийского моря. (Расшир. тезисы) - Ломонос. чтения. Секц. геологии. М. 2005. (соавторы: Левченко О.В., Мерклин Л.Р. и др.)
- Комплексные сейсмоакустические исследования площадок под установку буровых платформ в Каспийском море. Тез. VII-ой Международной научно-практической конференции Геомодель-2005. (соавторы: Левченко О.В., Мерклин Л.Р., Соколов С.Ю.)
- Опыт сейсмоакустического профилирования с многократным перекрытием - возможности и ограничения. Тез. Междун. конферен. Инженерная геофизика-2006. Геленджик, 2006. (соавторы: Кузуб Н.А., Лупырь Р.Р.)
- Комплексное использование различных технологий скважинных сейсмических наблюдений при инженерно-геологических исследованиях. Там же. (соавтор: скворцов а.г.)
- Многокан. сейсмоакуст. профилир. на акваториях - ожидания и результаты. Расшир. тез. Междун. конферен. Санкт-Петербург-2006. (соавторы: Кузуб Н.А., Лупырь Р.Р.)
- Сейсмоакустические исследования ледниковых отложений в Карском море. Расшир. тез. Междун. конферен. Нефть и газ Арктического шельфа-2006. Мурманск, 2006. (соавторы: Поляк Л. В., Гатауллин В. Н., Зверев А. С.)
- Многоканальное сейсмоакустическое профилирование в Байдарацкой губе. Там же. (соавторы: Кузуб Н.А., Клещин С.М. и др.)
- Некоторые результаты сейсмоакустических исследований в центральной части Каспийского моря. Там же. (соавторы: Левченко О.В., Мерклин Л.Р.и др.)
- Сейсморазведка (Руководство к практическим занятиям по курсу сейсморазведка). Издательство МГУ. 2006. 148 с.
- Многоканальное сейсмоакустическое профилирование в инженерных целях - больше возможностей, больше информации. Тез. Междун. научно-практ. конферен. Инженерная и рудная геофизика-2007. Геленджик, 2007. (соавторы: Кузуб Н.А., Токарев М.Ю.)
- Комбинированное сейсмоакустическое профилирование спаркер+бумер - реальные возможности. Там же. (соавторы: Зверев А.С., Токарев М.Ю., Росляков А.Д.)
- Расчет многоканальных систем сейсмоакустического профилирования на моделях. Сб. трудов каф. общ. и прикл. геофизики. Университет Дубна. М., РАЕН, 2007. с. 55 - 68.
- Особенности обработки данных многоканального сейсмоакустического профилирования на акваториях. Там же. с. 42 - 54 (соавторы: Кузуб Н.А., Токарев М.Ю.)
- О природе лярких пятен на временных разрезах сейсмоакустического профилирования. ГЕОРазрез, Электронный научный журнал университета Дубна, 2008.
- Технология сейсмоакустических исследований на акваториях - от двухчастотного профилирования к многоканальным системам. Расшир. тез. Междун. конферен. Санкт-Петербург-2008. (соавторы: Токарев М.Ю., Зверев А.С., Росляков А.Д.)
- Характеристики лярких пятен на сейсмоакустических разрезах. Тез. междун. научно-практ. конф. Инженерная и рудная геофизика-2008. Геленджик, 2008.
- Seismic characteristics of the Black Sea and the Eastern Mediterranean mud volcanoes. Proc. for the Intern. Earth Sc. Colloq. of the Aegean regions. Izmir. 1996. (co-author: cifci g.)
- Investigations of the mud volcanoes at the Eastern Extension of the Mediterranean Ridge. Proc. of the 11-th Petrol. Congress and Exibition of Turkey, Ankara, Turkey, 1996, p. 49 - 57. (co-authors: Cifci G., Dimitrov L., Limonov A.)
- Features of the Black Sea mud volcanoes and mud diapirs according to seismic data. (Abstr.) - Congress УGas and Fluids in Marine SedimentsФ. Amsterdam, 1997.
- Seismic Investigations of Carbonate Mounds in the North Atlantic to the West of Ireland. (Abstr.) - Congress УCarbonate Mud Mounds and Cold Water ReefsФ. Gent, Belgium. 1998.
- General characteristics of the Black Sea mud volcanoes and gas hydrates УshieldФ in the eastern Mediterr. sea // Annales Geophysicae. Suppl. I to V. 16. 1998. P. 296. (co-author: Cifci G.)
- High-resolution seismic survey in the Barents and Kara Seas: implication for the last-glaciation history. (Abstr.) - Third Workshop on Land Ocean Interactions in the Russian Arctic (LOIRA). Moscow. 2000. (co-authors: Levchenko O., Merklin L. et al.)
- Seismic mapping on the Thrace shelf (SW Black Sea); implications for recent sea level fluctuations (Abstr.) - 3rd Balkan Geological Congress and Exibition. Sofia, Bulgaria. 2002. (co-authors: Alpar B., Yaltirak C.)
- Late Pleistocene glaciation history of the southwestern Kara Sea. XVI INQUA Congress, Program with Abstracts, 2003, p. 71-72. (co-authors: Gataullin V., Polyak L. et al.)