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

На правах рукописи

Гошко Елена Юрьевна

ВЫЧИСЛИТЕЛЬНАЯ ТЕХНОЛОГИЯ ИЗУЧЕНИЯ ГЕТЕРОГЕННЫХ СРЕД ЗЕМНОЙ КОРЫ ПО ДИНАМИЧЕСКИМ ХАРАКТЕРИСТИКАМ ЛОКАЛЬНЫХ ВОЛНОВЫХ ПАКЕТОВ (по данным цх>фильных глубинных сейсмических наблодений МОВ-ОП)

Специальность 05.13.18 - Математическое моделирование,

численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени кандидата технических наук

Новосибирск- 2006

Работа выполнена в Новосибирском государственном архитектурно- строительном

университете НГАСУ (Сибстрвн)

Научный руководитель: Доктор технических наук, профессор

С.М.Зеркаль

Официальные оппоненты: Доктор технических наук

Глинский ЕМ.

Кандидат физико-математических наук ЧевердаВ.А

Ведущая организация: Югорский научно-исследовательский

институт информационных технологий

Защита диссертации состоится 19 декабря 2006 года в 10 часов на заседании диссертационного совета Д 003.061.02 при Институте вычислительной математики и математической геофизики СО РАН по адресу: г. Новосибирск, пр. акад. Лаврентьева, 6.

С диссертацией можно ознакомиться в читальном зале библиотеки Института вычислительной математики и математической геофизики СО РАН

Автореферат разослан 17 ноября 2006 г.

Ученый секретарь диссертационного совета

д. ф.-м. и. Сорокин СБ.

Общая характеристика работы Актуальность работы

Бурное развитие средств вычислительной техники и связанных с вей информационных технологий обеспечивает рост эффективности современной производстве ни ой базы геофизических исследований. Сейсмические методы исследования имеют преимущество перед другими геофизическими методами по детальности и достоверности, получаемой при их использовании, информации о структуре слоев земной коры. Открытие большинства нефтяных и газовых месторождений на суше и в морских акваториях стало возможным благодаря высокой эффективности сейсморазведки. Роль вычислительной техники и информационных технологий в сейсмических исследованиях трудно переоценить учитывая, что уже во второй половине прошлого столетия объемы обрабатываемых сейсмических данных оценивались астрономическими числами [Алексеев A.C., Цибульчик Г.М Актуальные проблемы вычислительной математики и математического моделирования. - Новосибирск: Наука, 1985г.].

Результатами сейсморазведочиых работ и обобщающими геологическими исследованиями последних лег доказана принципиальная возможность открытия промышленных углеводородных скоплений на больших глубинах (более 4-5 км), в том числе значительных по запасам. Дальнейшая задача прогноза и поисков пгубокозалегающих залежей УВ зависит от степени изученности условий и факторов формирования и размещения таких объектов. В этой связи проблема выявления особенностей строения земной коры, тектонических обстановок в ее глубоких частях, влияющих на формирование залежей УВ на больших глубинах, приобретает особую актуальность. Полученные в последнее время многочисленные глубинные сейсмические данные показали, что строение консолидированной земной коры, обусловленное неоднородностями различной геологической природы, отличается от горизонтально слоистого. На глубинных разрезах ОГТ не установлено отражающих границ, расположенных ниже 4-5 сек и имеющих региональную протяженность. Наблюдаются лишь хаотически расположенные отражающие площадки с длиной, соизмеримой с первой зоной Френеля. Они образуют вертикально и латерально неоднородные зоны насыщения коры отражателями. В этом случае важно выбрать модель, соответствующую выявленному строению геологической среды, и Применять адекватные алгоритмы обработка сейсмических данных с целью получения информативных характеристик земной коры.

Общим подходом к обработке глубинных сейсмических данных в настоящее время является то, что в основу закладывается гетерогенная сейсмическая модель земной коры [Караев H.A., Рабинович ГЛ. Рудная сейсморазведка. Москва: ЗАО «ГеоинформмаркяДЮО!,

Процесс обработки сейсмических данных при этом ориентируется на статистические методы вычисления интегральных характеристик разреза: динамических, энергетических, спектральных. Эта характеристики используются дат последующего многофакгорното разделения полей сейсмических атри&ггов не области однородности. Как было отмечено [Пузырен Н.Н. Методы в объекты сейсмических исследований. Н-ск Иэд-во СО РАН, 1997], «...в сейсмологии отчетливо обозначилась тенденция возможно полного использования динамических характеристик волн (амплитуд, спектров, импульсов и ¿р). Это направление, безусловно, должно получить еще большее развитие в будущем». Изучению динамических характеристик сейсмических данных способствует в настоящее время наличие хорошо калиброванных сейсмограмм, предоставляющих данные об отражениях б истинных амплитудах. Именно это обстоятельство делает возможным волновой подход, при котором в качестве исходных данных используются: форма отраженного волнового пакета, время и знак вступления, максимальная амплитуда, спектральные и энергетические характеристики отраженной волны.

Одним из ключевых способов получения информации о свойствах объектов неоднородной среды является выделение и исследование формы отраженных волновых пакетов (сигналов) на разрезе ОГТ. Чувствительным параметром, характеризующим свойства неоднородных горных пород, принято считать затухание амплитуды сигнала. Это явление связано с повышенной 1р«ц1£Щ>ватостък>, разломами и разрушениями в среде, а также заполнением перового пространства флюидами различного происхождения, что делает затухание амплитуды сигнала, поглощение энергии на высоких частотах актуальным явлением с точки зрения прогноза месторождений. 5 настоящей работе предложен алгоритм и вычислительная технология расчета куба спектральных амплитуд сейсмических сигналов в пространстве временной я латеральной координат сейсмического профиля и частоты (t, х, f). Особенностью предложенного технологического решения является вычисление локального волнового пакета произвольной формы на разрезе ОГТ в бегущем окне с размером, согласованным с масштабом выявляемых неоднородности. Визуализация куба спектральных амплитуд в виде среза на определенной частоте / - Cortst или в интересующем интервале частот Л £ / £ f, позволяет исследовать эффекты затухания амплитуд сейсмических сигналов. Обработка данных выполняется с использованием разработанного специализированного программного комплекса StrcamSDS [3, 4, б]. В настоящее время изучение глубинного строения Восточной территории России дня развития ее минерально-сырьевой базы стало приоритетной задачей Министерства Природных Ресурсов Российской Федерации. С 1994 г. и по настоящее время в России создается государственная сеть опорных геофизических профилей для обеспечения широкого круга недропользователей информацией

о строения недр. Актер диссертационной работа принимает участие в глубинных исследованиях, проводящихся в настоящее время на опорных профилях: 2-ДВ, 3-СБ.

В диссертации представлены результаты применения разработанного автором алгоритма н вычислительной технологии БЬ-еатЙОЗ .шея диагностики зон аномального затухания амплитуд спектров сейсмических сигналов в условиях гетерогенных ерел земной коры. Полученные результат! хорошо согласуются с прежними представлениями о геологии исследуемых территорий, а также добавляют новые данные о строении земной коры и верхней мантии, интерес к которым обусловлен современными работами на нефть, газ и некоторые виды твердых полезных ископаемых [5.9].

Цель исследований

Изучение строения земной коры по глубинным сейсмическим разрезам ОГТ на региональных и опорных профилях требует разработки специальных приемов обработки и анализа волновых нолей, формирующихся в сложных гетерогенных срезах.

Целью работы является сознание алгоритма вычисления локальных пакетов отраженных волн по данным профильных глубинных сейсмических наблюдений МОВ-ОГТ и разработка компьютерной технологии анализа их динамических характеристик, в частности, одной из наиболее важных, - затухания амплитуды волн, являющейся прогнозным признаком при локализации месторождений полезных ископаемых.

Задачи исследований

1. Разработка и обоснование алгоритма вычисления локальной фермы пакета отраженных, дифрагированных и рассеянных волн по данным профильных глубинных сейсмических наблюдений МОВ-ОГТ.

2. Численное исследование алгоритма определения формы локального волнового пакета на основе математического моделирования волновых образов сложных геологических сред с помощью синтетического разреза ОГТ.

3. Разрабопса технологии изучения динамических характеристик локальных волновых пакетов на глубинных разрезах ОГТ в условиях сложных гетерогенных сред земной коры.

4. Создание специализированного вычислительного программного комплекса ЯЬеатБОЗ и практическое использование его на материалах глубинных шорных профилей: 2-ДВ (г. Магадан - мыс Эммытаген), 3-СБ (Алтай - Северная Земля).

Фактический материал и методы исследований

Разработанное автором алгоритмическое и компьютерное обеспечение для исследования динамических параметров отраженных волн было опробовано: — на дальневосточном опорном профиле 2-ДВ (г. Магадан - мыс Эммытаген, 2000

погонных километров, длительность регистрации наблюдений 25 сек); — на опорном профиле Восточной Сибири 3-СБ (рассечка к профилю Алтай-Северная Земля, 300 пог. км, длительность регистрации наблюдений 24 сек);

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

H.ИБерзон, Н.Н.Пузыревым, Аки и Ричардсом, Н.АКараевым, О.Л.Кузнецовым, методы математического анализа и статистик».

Основные положений, выносимые на защиту

I. Алгоритм вычисления локалыюК формы пакета отраженных, дифрагированных и рассеянных вола, формирующейся на разномасштабных неоднородностях в условиях гетерогенных сред земной коры, по данным профильных глубинных сейсмических наблюдений МОВ-ОГТ.

2. Технология исследования динамических (спектральных, энергетических) характеристик локальных волновых пакетов в условиях сложных гетерогенных сред земной коры.

3. Специализированный вычислительный программный комплекс SireamSDS, позволяющий:

— создавать и визуализировать матрицу локальных волновых пакетов,

— формировать куб амплитудных спектров локальных волновых пакетов в пространстве трех переменных t, х, f,

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

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

4. Результаты применения специализированного программного комплекса на фрагменте дальневосточного опорного профиля 2-ДВ, в а профиле 3-СБ в Восточной Сибири.

Научная новизна работы

Автором разработан алгсритм вычисления локальной фермы волнового пакета но данным профильных глубинных сейсмических наблюдений МОВ-ОГТ. Это позволило впервые получить цифровые данные об изменении формы волнового пакета отраженных,

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

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

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

Практическое использование разработанного специализированного программного комплекса StreamSDS на материалах глубинных опорных профилей: 2-ДВ (г. Магадан - мыс ЭмыытагенХ 3-СБ (Алтай - Северная Земля) позволило впервые составить предварительные геолого-геофизические и тектонические модели строения земной коры этих территорий.

Внедрение результатов работы

Специализированный программный комплекс StreamSDS внедрен в ФГУП «СНИИГТИМС», что подтверждено Актом о внедрении (Приложение 1).

Проводящиеся в настоящее время автором работы по обработке и интерпретации материалов СГ-ОГТ опорного геофизического профиля 2-ДВ (г. Магадан - мыс Эммытаген, 2000 погонных километров) с использованием защищаемых алгоритмических и технологических разработок имеют положительные отзывы ведущей организации ВСЕПЕЩ-.Санкг-Петербург. Рекомендация к дальнейшему применению вычислительного программного комплекса StreamSDS зафиксирована в Решениях Комиссии Министерства природных ресурсов РФ по надзору за выполнением работ (от 22.09,06 г.) по Договору о научном сопровождении региональных геодого-геофизпческих работ на территории Магаданской области и Чукотского автономного округа.

Созданная автором вычислительная технология, реализованная в комплексе StreamSDS применяется в настоящее время для выполнения договоров с Федеральным агентством по недропользованию РФ, Территориальным агентством по недропользованию по Красноярскому краю по созданию опорного профиля 3-СБ «Алтай — Северная Земля».

Апробация работы

Основные результаты по теме диссертационной работы доложены на:

1. Всероссийском се мин ере «Методы, технические средства, методика обработки н интерпретации геолого-геофизических исследований ори создании государственной сргц опорных геофизических профилей», г. Новосибирск, 10-13 сентября 1999 г.,

2. Геофизической научно-практической конференции: «Проблемы повышения эффективности применения геофизических исследований при поисках, разведке, разработке и эксплуатации месторождений нефти и газа в Западной Сибири»! Г, Тюмень, 16-18 сентября 2003г.

3. 62-Й научно-технической конференции, посвященной 75-летию НГАСУ (СИЕСТРИН). - Новосибирск НГАСУ, 2005

4. VII Международной конференции «Новые идеи в науках о земле» в Москве, 2005

5. 2-ом Международном симпозиуме «Активный геофизический мониторинг литосферы земли, 12-16 сентября, Академгородок, Новосибирск, 200$ г.

6. Международной конференции «Вычислительные к информационные технологии в науке, технике и образовании» Казахстан, Павлодар, 20-22 сентября 2006 г.

Публикации

По теме диссертации опубликовано 11 печатных научных работ. Из них одна монография и три статьи в рекомендованных ВАК научных журналах.

Личный вклад автора диссертации

В работах 3, 4 научному руководителю Зеркалю С М. принадлежит постановка задач и руководство работой. В работе б Маркову В,М. принадлежит программная реализация разработанных Гошко ЕЮ. алгоритмов. В работе 9 Мшурскому A.B. принадлежит геологическая и тектоническая интерпретация результатов применения специализированного программного комплекса StteamSDS на материалах опорного профиля 2-ДВ. В остальных перечисленных работах соавторам Рудницкой ДИ., Сагайдачной О.М., Сальникову принадлежит предоставление цифровых сейсмических данных по глубинным разрезам OTT опорных профилей 2-ДВ и 3-СБ и геологический анализ результатов. Работы 1, 2, 10, 11 посвящены созданию я применению системы РЕАПАК; руководящая роль принадлежит Рудницкой Д.И.. алгоритмическое обеспечение отдельных блоков системы, тестирование программ и их использование — Гошко Е.Ю.

Структура и объем работы

Диссертация состоит из введения, 4-х глав, заключения и 2-х приложений, содержит 133 страницы машинописного тенета и 44 рисунка. Библиография включает 79 наименований.

Введение

Во введении излагается современное состояние исследуемой проблемы, обосновывается актуальность. Научная новизна й важность полученных в диссертации результатов, определяется их место среди близких научнцх исследований. Глава 1. Физические и математические основы сейсмических исследований строения земной коры

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

Сейсмические модели по существу являйтся математическими. Изложим основные положения теории, следуя [Аки. 3]. Физические предпосылки для создания модели содержат:

— уравнение движения под действием объемной силы / = f{X,t): Э1 3

р—j-Uj = /(где и «u{x,í) - функция смещения, г, - компоненты тензора ct dxj ' '

Напряжений в точке (X,í);

— современное обобщение закона Гука, заключающееся в утверждении, что каждая компонента тензора напряжений есть линейная комбинация всех компонент тензора деформаций: = с^ ■ ен, где ев - компоненты тензора деформаций в точке (Л\г), величины ст - упругие постоянные, независимые от деформаций, Наиболее продуктивной моделью в настоящее время считается модель Земли как вергнкаявдо-неоднородйого изшроиного полупространства или шара. В этом случае модель неоднородна в вертикальном направлении и упругие постоянные си выраясаются через константы Ламе ém =А8дЗя -t-fÁS^Sft при этом упругие постоянные и плотность среды являются Изменяющимися функциями глубины: Л = ft = я(*)> р =f />(*) ■

— первый закон термодинамики, утверждающий, что тело обладает внутренней энергией Е, которая может изменяться при его деформировании: dE =dR±if-de^, где (¡R-приток тепла извне. То есть, скорость увеличения внутренней энергии = скорость разогрева + скорость выполнения механической работы. Если приток тепла пренебрежимо мал, то dE^T-dS+Vg -ífey, где Т—абсолютная температура, S - энтропия в единице объема. Таким образом, энтропНй и компоненты деформаций являются параметрами, определяющими внутреннюю энергию среды.

Земля - неидеальва упругая среда а том смысле, что малые движения частиц, вызываемые распространяющейся волной, ведут к необратимым потерям энергии волн и>за наличия множества различных механизмов диссипации. Вызванное диссипацией затухшие распространяющейся волны удобно описывать с помощью параметра <Хдобротности):

1 М -

—г—г =--, или если амплитуда волны А=Е3 (при линейном соотношении между

у (да) 1жЕ

напряжением и деформацией) —т-г = . Получаем экспоненциальное затуханиие

<Л®) хА

амплитуд сигнала во времени и в пространстве Л(х)-*Лц

Глава 2, Обратная задача в сейсмике отраженных волн

В данной главе дается краткое описание известной методики суммирования по ОГТ и ее модификации в случае сверхглубинных данных, содержащих наблюдения о приходе волн до времени 50 сек. Пригодится точное решение обратной задачи сейсмики в случае нормального падения волн на плоско-параллельную однородную изотропную среду, а также формулы для вычисления коэффициентов отражения и преломления продольной волны на границе двух сред. Описан граф современной динамической обработки для получения итогового суммарного разреза ОГТ «в истинных амплитудах».

В $ 2.3 представлено решение обратной задачи сейсмики на отраженных волнах в случае их нормального падеиия на однородную изотропную сроду.

Уравнение движения для одномерного случая: р(х)■ = ~'где плотность,

В(х)- упругая постоянная, равная дня продольных волн сжатия А(х)+2р(х), а Л(ог)-посгоянная сжатия Ламе, постоянная сдвига Ламе.

Если среда состоит из отдельных однородных слоев, то для ¡-го слоя уравнение

движения примет вид " ® 1 где В,,/?, - константы слоя. Тогда решение

запишется в виде

Удобно выбрать мощность слоев &х, таким образом, чтобы время пробега в каждом слое At, оставалось постоянным: А/, = const. Такая дискретная среда удобна в том отношении, что смещение н(г)можно выразить в виде последовательности импульсов, существующих только в дискретные моменты времени, разделенные промежутком Ai:

ы(г)=^ «„<?(/ -2 лД/). Преобразование Фурье будет нмегь вид

я

U(m) ~ |«(/)ехр(ш)Л = Xм» ехр(2/лй>Д/). Если положить j=ехр(2/«>Д(), то получим

выражение £/(z)=, называющееся z-преобразованием. ■

В § 2.5 рассматривается преобразование данных в современном процессе динамической обработки сейсмического разреза ОГТ. Предполагая, что в случае близвертикальных отражений все тины волн имеют одинаковую форму, записано математическое выражение сейсмической трассы ОГТ:

где av = а^, -«у - независимые от времени коэффициенты нормировки, (i) •

функции, выражающие затухание энергии полезных отраженных волн и волн-помех соответственно, K^iO.K^ (l) - коэффициенты отражения для сднократпо-отраженных и многократных, дифрагированных, рассеянных, обменных воли всех типов соответственно.

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

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

>.,<i) = S0(0*K,(0. w x = (i+J)/2.

Глава 3. Изучение гетерогенных сред земной коры по динамическим характеристикам локальных волновых пакетов

В данной главе доказывается необходимость совершенствования математических моделей, применяемых для гетерогенных сред. Модель выступает как средство адекватного качественного описания изучаемых объектов и как мера количества информации, извлекаемой из них. Приведены примеры отображения иа сверхглубинном сейсмическом разрезе МОВ-ОГТ гетерогенной земной коры. Вводится понятие локального пакета отраженных, дифрагированных и рассеянных волн. Предложен алгоритм вычисления локальной формы волнового пакета на разрезе ОГТ. Основные идеи разработанного алгоритма определения локальной формы пакета отраженных, дифрагированных и рассеянных волн:

— Алгоритм основав на синфазном суммировании отдельных волновых пакетов f(t), выявленных в пределах участка сейсмического разреза (окна размером , х0), соответствующего области D неоднородности с характерным размером а, которая служит объектом геологического интереса. Алгоритм исключает зависимость полученной итоговой формы волнового пакета 5.(0 от взаимного геометрического расположения отдельных волновых пакетов /(<) в пределах фрагментов трасс /(/) сейсмического разреза, так как вое волновые пакеты/(/) суммируются синфазно.

— В качестве точек синфазного суммирования волновых пакетов /(() во фрагменте сейсмической трассы /(г) выбираются положении максимумов огибающей модуля сейсмической трассы [/(/) {r,:¡ /(!,) |-»тах}. Локальные максимумы

огибающей соответствуют максимальным амплитудам отдельных волновых пакетов. Их разделение обеспечивается наличием окрести ости у каждого локального максимума. Временная длительность волнового пакета /(í) задается окрестностью, шгорая окружает точку суммирования в равных или неравных пропорциях: (/, — с,;/, где е,,ег-радиусы окрестности.

— Выделенный волновой пакет /(/) нормируется на максимальную амплитуду

шах/((,)

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

1 . —

волнового пакета на локальном участке сейсмического разреза S„(l) = —(0 • гае п -

п м

количество волновых пакетов в пределах окна вычислений.

— Далее в цикле осуществляется смещение окна вычисления локальной формы волнового пакета с шагом Дх по координате х и Д( по времени t вдоль сейсмического разреза. Размер шагов смещения окна задается с учетом обеспеченна необходимой гладкости изменения фермы

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

Алгоритм позволяет вычислить локальную форму отраженной волны в пределах некоторой области D разреза ОГТ и по ее спектральным и энергетическим характеристикам судить об отражающих свойствах гетерогенной среды в пределах области неоднородности на разрезе ОГТ.

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

Известен способ определения формы элементарного сигнала, реализованный в системе РеапакРД [Рудницкая Д.И. Методические рекомендации по обработке сейсмических записей. Новосибирск; СНИИГГиМС, 1983; Рудницкая Д.И., Гошко В.Ю. и др., 1], Способ позволяет оценить форму сигнала по временному разрезу ОГТ с помощью суммирования волновых пакетов, выделенных на трассах. Основанием для выделения момента вступления сейсмического сигнала является тот факт [Берзон И.С. Динамические характеристики сейсмических волн. М.: Иэд-во АН СССР, 1962], что сейсмический сигнал характеризуется плавным нарастанием дервого отклонения, то есть начальная скорость смещения частиц среды равна нулю, а ускорение или производные от смещения еще более высокого порядка претерпевают скачок.

Способ, реализованный в системе РеапакРД обладает радом ограничений. Выявление момента вступления сигнала по критерию максимального градиента амплитуд трассы оправдало лить для минимально-фазовых сигналов. Для сигналов других типов, например, близких к симметричным (нульфазовым) или автнеимметричным, такой способ приводит к ошибкам в оценке фазовой характеристики сигнала и, следовательно, к ошибкам в определении момента и амплитуды вступления волны [Кондратьев ПК. Опенка точности д разрешающей способности используемых на практике способов решения обратной динамической задачи сейсморазведки,- Москва; Обзор ВИЭМС, 1990].

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

Далее в параграфах этой главы приводятся результаты моделирования сейсмического разреза ОГТ и модальное тестирование разработанного алгоритма вычисления формы локальных волновых пакетов.

Глава 4. Вычислительный программный комплекс «StreamSDS» и результаты его практического использования

Б данной главе дается описание вычислительной технологии «StreamSDS» диагностики зон затухания энергии отраженных ваян но глубинным разрезам ОГТ на опорных профилях, приводится «Руководство пользователя» пакетом программ. Результаты проведенных вычислений формы локальных волновых пакетов в режиме «скользящего окна» на материалах оперного профиля 2-ДВ убедительно доказывают существование значительных ее изменений в пределах сверхглубинного разреза ОГТ.

Реализация технологии StreamSDS на материалах профиля 2-ДВ для детального исследования аномального затухания амплитуд локальных волновых пакетов заключалась в анализе. спектрального куба для частот от f » 18 Гц до f = 40 Гц с шагом по частоте 2 Гц. Построенные частотные срезы спектрального куба были визуализированы в виде цветовых полей. Анализ этих изображений позволил определить интервал частот <18 - 23 Гц) и контуры зон на сейсмическом разрезе, где уверенно проявляется обнаруженный эффект затухания амплитуд. Приведена предварительная геолого-тектоиическая модель фрагмента геограверса 2-ДВ {1100 - 1460 км), построенная на основании анализа поля амплитуд в интервале средних и высоких частот Мигурским A.B.

Технология StreamSDS реализована также на материалах глубинного сейсмического разреза ОГТ рассечка к опорному профилю 3-СБ (Алтай-Северная Земля), пересекающего крупные структурно-тектонические элементы Сибирской платформы. Использование поля амплитуд средних и высоких частот спектра локальных волновых пакетов помогает более уверенно выделить и проследить отражающие и проводящие горизонты внутри кристаллического фундамента, в рифейском и палеозойском структурно-тектонических комплексах осадочного чехла; изучить строение земной коры, заяартировать поверхность Мохоровичича, увязать их с процессами нефгегаэообразова-ни* и рудообразования в осадочном чехле.

Заключение

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

Численные расчеты, проведенные с помощью описанного алгоритма и разработанного на его основе специализированного вычислительного комплекса StreamSDS на фрагменте опорного сейсмического профиля ОГТ, показывают адекватные (с точки зрения

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

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

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

Приложение

В приложении 1 представлен Акт о внедрении вычислительного программного комплекса StreamSDS в ФГУП «СНИИГТиМС», подтверждающий его апробирование на реальных сейсмических данных и использования при выполнении Договоров Института с Министерством Природных Ресурсов РФ. В приложении 2 представлено Свидетельство об официальной регистрации программы для ЭВМ «РеапакРД», одним из правообладателей которого является автор диссертационной работы [1].

Основные результаты диссертационной работы

1. Разработан алгоритм вычисления локальной формы отраженного волнового пакета по данным профильных глубинных сейсмических наблюдений МОВ-ОГТ. Ои позволяет изучать изменяющуюся форму сейсмического отражения в сложных гетерогенных средах.

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

волновых пакетов на синтетическом разрезе ОГТ позволило восстановить заданную форму сейсмических импульсов и их частотные характеристики. Исследована устойчивость алгоритма, показано, что достаточно около ста суммирований в пределах окна вычислений ДНЯ получения стабильной формы локального отражения.

3. Разработана технология изучения динамических (спектральных, энергетических) характеристик локальных волновых пакетов в условиях сложных гетерогенных сред земной коры

4. Создан вычислительный программный комплекс ЗиеатБОЗ, позволяющий:

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

— формировать куб амплитудных спектров локальных волновых пакетов в пространстве трех переменных: временной н латеральной координат сейсмического профиля и частоты (I,

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

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

5. Разработанное автором алгоритмическое н компьютерное обеспечение для исследования динамических параметров отраженных волн было внедрено во ФГУП «СНИШТиМС» в соответствии с Актом о внедрении (см. Приложение 1). Вычислительный комплекс ЗйеатЗОЗ был применен для обработки и интерпретации материалов по оперным сейсмическим профилям Восточной Сибири и Дальнего Востока в рамках выполнения договоров ФГУП «СНИИГТиМС» с Министерством Природных ресурсов, Агентством по недропользованию РФ:

— опорному профилю 2-ДВ (г. Магадан - мыс Эммытаген), 2000 погонных километров, длительность регистрации наблюдений 25 сек;

— рассечки к профилю 3-СБ «Алтай-Северная Земля», 300 погонных километров, длительность регистрации наблюдений 24 сек.

Благодарности

Автор выражает глубокую благодарность и признательность своему научному руководителю доктору технических наук, профессору С.М. 3 ер калю. За сотрудничество в области программных разработок автор выражает глубокую благодарность В.М. Маркову, За постоянную поддержку в работе автор признателен к.г.-м.н. А.С.Сальникову, аа ценные советы, консультации, критические замечания благодарит к.т.н. О.М. Сагайдачную. Автор благодарен к.г.-и.н. ДИ. Рудницкой за переданный геофизический опыт во время многолетней совместной работы. Автор выражает благодарность за труд геологического осмысления получетных в диссертации результатов д.г.-м.н. AB. Мигурскому, дг.-м.н. С. С. Долгушину.

Автор признателен сотрудникам отдела сейсморазведки и отдела тектоники ФГУП СНИИГГвМС за совместный труд и помощь в реализации разработок.

Основное содержание диссертации опубликовано в следующих работах:

1. Авторское свидетельство Ks 990128 «свидетельство об официальной регистрации программы для ЭВМ», - Рудницкая ДИ., Горячев Д.Н., Корнилов М.В., Гошко ЕЮ, в др. -Регистрация 12.03.99 г.

2. Васильев В,П., Рудницкая ДИ., Горячев ДН., Гошко Е.Ю., Корнилов М.В., Лурье АХ., Сагайдачная О.М., Стражникова И.Н., Технология и результаты построения цифровой себсмогеологической модели методами системы РЕАПАК по геотраверсу УралсеЙс. -Глубинное строение и геодинамика Южного Урала (проект УралсеЙс). Монография. Тверь: ИздательствоГЕРС,2001.-е. 35-42

3. Гошко Е.Ю., Зеркаль СМ. «Вычислительная технология диагностики зон затухания анергии отраженных волн по данным глубинных сейсмических разрезов ОГТ». - Труды международной конференции «Вычислительные и информационные технологии в науке, технике и образовании». • Павлодар, Казахстан, 20 - 22 сентября 2006 г, I том, С. 359-371.

4. Гошко ЕЮ., Зеркаль C.tvt Алгоритмическое обеспечение исследования гетерогенных сред земной коры по динамическим характеристикам локальных волновых пакетов. Вестник НГУ. Том 4, Выпуск 1, Серия: Информационные технологии, 2006 г., с. 33-43.

5. Гошко Е.Ю., Рудницкая ДЛ, Сагайдачная О.М, Сальников A.C. Сейемоакустический образ земной коры по геотраверсу 2-ДВ в систше «РЕАПАК» - VII Меж1(увародная конференция «Новые идеи в науках о земле»: Материалы докладов. Т.2. - М.: КДУ, 2005. -335 с. КДУ—книжный дом университет, с. 288.

6. Е.Ю. Гошко, BJví. Марков, ДИ. Рудяинкая, О.М, Сагайдачная, A.C. Сальников,

Компьютерная технология расчета нолей сейсмических атрибутов в системе РЕАПАК доя создания сейсмоакустических образов глубинных разрезов земной коры. Вестник НГУ, Том 2, Выпуск I, Серия: Информационные технологии, 2005 г., с. 11-19.

7. ЕЮ. Гошко, Д.И Рудницкая, О.М. Сагайдачная, А.С. Сальников. Компьютерная технология обработки глубинных сейсмических данных с использованием методики РеапакРД - Тезисы докладов 62-й научно-технической конференции, посвященной 75-летию НГАСУ (СИБСТРИН). -Новосибирск: НГАСУ. 2005., (с. 119-124).

8. ЕЮ. Гошко, ДИ. Рудницкая, О.М. Сагайдачная, А.С. Сальников. Построение сейсноакустических образов глубинных разрезов земной коры в системе РеапакРД (на примере опорного профиля 2-ДВ). - 2-й Международный симпозиум «Активный геофизический мониторинг литосферы земли»: Материалы международной конференции 1216 сентября, Академгородок, Новосибирск, 2006 г., с 242-247.

9. Мигурский А В., Гошко Е.Ю., Рудницкая Д.И., Сагайдачная ОМ. Сейсмогеологическая модель земной коры по геотраверсу Уралсейс. (Результаты интерпретации в системе РЕАПАК). - Глубинное строение и геодинамика Южного Урала (проект Уралсейс). Монография. Тверь: Издательство ГЕРС, 2001. - с. 209-214

10. Рудницкая Д.И., Биезайс Я,Я., Гошко ЕЮ. «Применение системы РЕАПАК в решении задач региональной сейсморазведки». - «Методы, технические средства, методика обработки в интерпретации геолого-геофизических исследований при создании государственной сети опорных геофизических профилей», Доклады Всероссийского семинара, г. Новосибирск, 10-13 сентября 1999 г, отв. Ред. Сурков B.C., Васильев В.П,

11. Рудницкая Д. И., Корнилов М.В., Гошко ЕЮ. «Опыт использования системы РЕАПАК при изучении месторождений нефти и газа в Западной Сибири». — Научно-технический журнал ЕАГО, Геофизика, 1996 г.,№ 3, с. 19-24.

Новосибирский государственный архитектурно-строительный университет

(Сибстрин)

630008, г.Новосибирск, ул. Ленинградская, 113 Отпечатано мастерской оперативной полиграфии НГАСУ (Сибстрин)

Тираж И Заказ ш

Введение.

Глава 1. Физические и математические основы сейсмических исследований строения земной коры.

§1.1. Физические предпосылки существующей модели сейсмических данных и физическая сущность используемых параметров.

§ 1.2. Функция упругой энергии в случае плоских волн.

§ 1.3. Распространение волн в поглощающей среде.

§ 1.4. Классификация проблем рассеяния.

Глава 2. Обратная задача в сейсмике отраженных волн.

§ 2.1. Суммирование в методе общей глубинной точки (ОГТ).

§ 2.2. Сверхглубинная модификация метода ОГТ.

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

§ 2.4. Коэффициенты отражения и преломления (прохождения) продольной волны на границе двух сред.

§ 2.5. Преобразование данных в современном процессе динамической обработки сейсмического разреза ОГТ.

Глава 3. Алгоритмическое обеспечение исследования гетерогенных сред земной коры по динамическим характеристикам локальных волновых пакетов.

§ 3.1. Гетерогенность земной коры и ее отображение на сверхглубинном сейсмическом разрезе МОВ-ОГТ.

§ 3.2. Алгоритм вычисления локальной формы волнового пакета на разрезе

§ 3.3. Процедуры моделирования сейсмического разреза ОГТ и модельное тестирование алгоритма вычисления формы локальных волновых пакетов.

Глава 4. Вычислительный программный комплекс «StreamSDS» и результаты его практического использования.

§ 4.1. Вычислительная технология «StreamSDS» диагностики зон затухания энергии отраженных волн по разрезам ОГТ на опорных профилях.

§ 4.2. Специализированный программный комплекс «StreamSDS».

Руководство пользователя.

§ 4.3. Реализация технологии StreamSDS на материалах глубинных сейсмических разрезов ОГТ опорных профилей 2-ДВ и 3-СБ.

§ 4.4. Сейсмоакустический разрез как материал для исследования физических и геометрических свойств неоднородных объектов.

§ 4.5. Библиотека физических и геометрических атрибутов сейсмоакустического разреза.

§ 4.6. Расчет полей физических и геометрических атрибутов сейсмоакустических разрезов по опорным профилям 2-ДВ и 3-СБ.

Бурное развитие средств вычислительной техники и связанных с ней информационных технологий обеспечивает рост эффективности современной производственной базы геофизических исследований. Сейсмические методы исследования имеют преимущество перед другими геофизическими методами по детальности и достоверности, получаемой при их использовании, информации о структуре слоев земной коры. Открытие большинства нефтяных и газовых месторождений на суше и в морских акваториях стало возможным благодаря высокой эффективности сейсморазведки. Роль вычислительной техники и информационных технологий в сейсмических исследованиях трудно переоценить учитывая, что уже во второй половине прошлого столетия объемы обрабатываемых сейсмических данных оценивались астрономическими числами [Алексеев А.С., Цибульчик Г.М. , 4.]. Полевые сейс-моразведочные работы относятся к категории дорогостоящих, полученные в их результате данные являются уникальными, т.к. полностью воспроизвести сейсморазведочный эксперимент не представляется возможным и, кроме того, объем производимых при обработке сейсморазведочной информации вычислительных работ весьма велик. В этой связи важно иметь объективные оценки качества и строгие формулировки физико-математического и геологического содержания полученных результатов. При этом следует иметь в виду, что факты обнаружения месторождений могут быть результатом случайного успеха, за счет большой статистики и высокой территориальной плотности изыскательских работ. Таким образом, важна оценка достоверности геологической интерпретации полученных результатов и факта открытия месторождений как стргого--следствия адекватности наших геологических знаний об исследуемом регионе и верности нашей физико-математической модели сейсмических процессов. Необходимость отмеченных оценок определяется сегодняшними тенденциями рассматривать информационно-вычислительные технологии, разрабатываемые и эксплуатируемые в различных областях человеческой деятельности, как наиболее естественное средство накопления знаний и дальнейшего уточнения математических моделей изучаемых объектов и явлений.

Физико-математической основой решения сейсмических задач и формулировки их математических моделей является идея о взаимнооднохначном соответствии между внутренним механическим строеним среды (Земли) и режимом колебаний ее поверхности при заданных источниках зондирующего сигнала. Эта идея реализуется в законе, записываемом в форме системы дифференциальных уравнений динамической теории упругости - уравнений Ламе. Модели образуются постановкой задач для указанных уравнений и производных от них в различных ситуациях достаточно типичных, чтобы составить определенный класс.

Результатами сейсморазведочных работ и обобщающими геологическими исследованиями последних лет доказана принципиальная возможность открытия промышленных углеводородных (УВ) скоплений на больших глубинах (более 4-5 км), в том числе значительных по запасам [34]. Дальнейшая задача прогноза и поисков глубокозалегающих залежей УВ зависит от степени изученности условий и факторов формирования и размещения таких объектов. В этой связи проблема выявления особенностей строения земной коры, тектонических обстановок в ее глубоких частях, влияющих на формирование залежей углеводородов на больших глубинах, приобретает особую актуальность.

В настоящее время происходит активное вовлечение в сейсморазведоч-ные работы глубоких горизонтов земной коры со сложным геологическим строением. Соответственно этому усложняются алгоритмы обработки и растет ее объем. В среднем методы обработки требуют выполнения около 105 операций; имеется тенденция к росту годового объема информации в десятки раз за каждые пять лет [Алексеев А.С., Цибульчик Г.М., 4]. В последнее время главными операциями в сейсморазведке стали: цифровая фильтрация, спектральный, спектрально-временной, энергетический анализ, автоматическое выделение и прослеживание формы сигналов.

Ведущее место как по информативности, так и по уровню промышленного внедрения в комплексе геофизических методов сохраняет сейсморазведка. Необходимо подчеркнуть, что в рамках существующих теоретических представлений, вследствие сложности математического аппарата, не все задачи могут быть удовлетворительно решены. Разработка более совершенного теоретического аппарата, которая уже успешно ведется, безусловно, позволит обосновать новые модификации сейсмических исследований, способные решать более тонкие задачи фундаментального и прикладного характера. К таким модификациям следует, в частности, отнести интенсивно развивающиеся направления в сейсмике - сейсмическую томографию (точнее ее кинематическую составляющую, как наиболее разработанную [20, 21]) и развиваемый в данной диссертации подход, основанный на исследовании динамических характеристик волновых пакетов. Динамический анализ в сейсмике стал возможным благодаря повсеместному использованию мощных интерференционных систем наблюдений, обеспечивающих эффективное накапливание полезной информации, широкому внедрению цифровой регистрирующей аппаратуры и современных ЭВМ с развитым программным и методическим обеспечением обработки и интерпретации сейсмических данных.

До недавнего времени сейсмические исследования, направленные на изучение глобальных границ раздела в земной коре, проводились в основном методом ГСЗ с использованием разреженных систем наблюдений. Однако, начиная с середины 80-х годов, в мировой практике для изучения глубинных границ в земной коре стали использоваться наблюдения отраженных волн в области, близкой к источнику возбуждения, с привлечением технологии МОВ-ОГТ - близвертикальные отражения. Сейсмический разрез, полученный этим методом, основан на многократном (порядка 100) суммировании записей от различных сейсмоприемников, позволяющим значительно усилить соотношение сигнал/помеха. В настоящее время современная регистрирующая аппаратура и методика проведения региональных сейсморазведоч-ных работ и работ по созданию опорных сейсмических профилей ОГТ позволяет регистрировать сейсмические волны, распространяющиеся на глубину, значительно большую, чем протяженность системы наблюдений на поверхности земли. Такой сейсмический разрез назван глубинным разрезом ОГТ. В настоящее время глубинные разрезы ОГТ содержат данные об отражениях в земной коре с временем прихода на ее поверхность, равным 50 и более сек. При этом «плотность» регистрации отражений, то есть кратность суммирования сейсмических трасс настолько высока, что мы можем рассчитывать на информативность, а не случайность отражений в трассе ОГТ даже при большом времени их регистрации.

В настоящее время регистрируются сейсмические наблюдения отраженных волн в земной коре с временем прихода до 50 сек и более. Современные глубинные разрезы ОГТ содержат информацию о строении всей консолидированной коры и верхней части мантии. Исследования последних лет позволили создать методики высокоразрешающей сейсморазведки МОВ-ОГТ, обеспечивающей существенное расширение полосы эффективно используемых частот на динамических разрезах ОГТ, повысить детальность изучения разреза по вертикали, установить тонкие особенности динамики волнового поля. С помощью этих методик осуществляется построение динамического разреза ОГТ, наилучшим образом соответствующего модели среды, возбуждаемой плоскими волнами при нормальном падении, и трансформация его в глубинный динамический разрез с восстановлением истинных соотношений интенсивности отражающих границ.

Полученные в последнее время многочисленные глубинные сейсмические данные показали, что строение консолидированной земной коры, обусловленное неоднородностями различной геологической природы, отличается от горизонтально слоистого. На глубинных разрезах ОГТ не установлено отражающих границ, расположенных ниже 4 - 5 сек и имеющих региональную протяженность. Наблюдаются лишь хаотически расположенные отражающие площадки с длиной, соизмеримой с первой зоной Френеля. Они образуют вертикально и латерально неоднородные зоны насыщения коры отражателями. В этом случае важно выбрать модель, соответствующую выявленному строению геологической среды, и применять адекватные алгоритмы обработки сейсмических данных с целью получения информативных характеристик земной коры.

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

Общим подходом к обработке глубинных сейсмических данных [Караев, Петров, Кузнецов О. Л., 24 - 28, 32, 48 - 52 ] в настоящее время является то, что в основу понимания строения глубинных зон закладывается гетерогенная модель земной коры. Гетерогенные среды, в отличие от слоистых, характеризуются размерами неоднородностей, соизмеримыми или меньшими, чем длина волны. При этом возникают дифрагированные и расеянные волны. В физике упругих колебаний эти волны называют соответственно рассеянным отражением резонансного типа и рассеянным отражением рэлеевского типа, в отличие от зеркального отражения, возникающего в слоистой среде

Кузнецов О.JI, 33]. Идея сейсмической гетерогенности отвечает неоднородному строению коры. Процесс обработки сейсмических данных при этом ориентируется на статистические методы вычисления интегральных характеристик разреза: динамических, энергетических, спектральных. Эти характеристики используются для последующего многофакторного разделения полей сейсмических атрибутов на области однородности.

В сейсмике параметры, характеризующие распространение сейсмических волн, принято разделять на динамические и кинематические. Кинематические параметры связаны с изучением фронтов, лучей, измерением времен распространения волн. К динамическим относятся параметры, характеризующие форму отраженных волн, во всем многообразии аспектов ее рассмотрения. Изучению динамических характеристик сейсмических данных способствует в настоящее время наличие хорошо калиброванных сейсмограмм, предоставляющих данные об отражениях в истинных амплитудах. Именно это обстоятельство делает возможным волновой подход, при котором в качестве исходных данных используются: форма отраженного волнового пакета, время и знак вступления, максимальная амплитуда, спектральные и энергетические характеристики отраженной волны. В связи с этим появилась возможность исследования качественно новых характеристик волновых полей - ускорений, напряжений, деформаций, наклонов. Одним из ключевых способов получения информации о свойствах объектов неоднородной среды является выделение и исследование формы отраженных волновых пакетов (сигналов) на разрезе ОГТ.

Опыт исследований последних лет позволил установить совокупность диагностических признаков - атрибутов сейсмических записей, которые являются количественными характеристиками наблюдаемых волновых полей. К настоящему времени [Chorpa S., 70] выделено более 50 сейсмических атрибутов, которые используются при интерпретации получаемых материалов для более полного описания геологических структур, свойств пород и флюидов. В эволюции атрибутов выделено несколько периодов, которые неразрывно связаны с развитием и компьютеризацией сейсморазведки. Период 1980 - 1990 гг. - применение цифровых регистрирующих и обрабатывающих систем привело к использованию динамических характеристик регистрируемых волн и соответствующих атрибутов. Были разработаны методики разделения разреза на однородные зоны и их количественной оценки, основанные на вычислении как интегральных, осредненных параметров поля и акустических свойств разреза (например, корреляционная методика прямых поисков, способы определения пластовых и интервальных скоростей и оценки вариаций статистических характеристик волновых полей и др.), так и дифференциальных, базирующихся на анализе узкой области волнового поля (например, «яркое пятно», итеративное сейсмогеологическое моделирование локальной зоны нефтегазонасыщения, псевдоакустический каротаж и др. [Авербух А.Г., 1; Михальцев А.В., Мушин И.А., Погожев В.М., 43].

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

Настоящий период 2000 - 2005 гг. характеризуется значительным увеличением объемов трехмерной сейсморазведки и сверхглубинной сейсморазведки на опорных и региональных профилях, сопровождающимся ростом числа атрибутов, связанных с углубленной характеристикой наблюдаемых волновых полей.

В сложных сейсмогеологических условиях объекты земной коры сильно маскируются влиянием различных глубинных и поверхностных факторов. Поэтому при обработке необходимо применять приемы подавления ложных аномалий, а при интерпретации - решать задачу определения надежности и распознавания природы аномалий. В последние годы созданы методики обработки и анализа волновых полей, позволяющие исключить или ослабить влияние некоторых факторов на оценки кинематических и динамических характеристик отраженных волн (итеративные алгоритмы коррекции динамики поля за счет влияния неоднородностей верхней части разреза, способы исключения влияния кривизны границ на оценки амплитуд волн и др. [Сергеев B.JL, Михальцев А.В., 19, Кондратьев O.K., 31, 32, Гогоненков Г.Н., 58]). Точная оценка кинематических и динамических параметров сейсмического разреза ОГТ является основой для создания достоверной модели строения гетерогенной геологической среды. При этом кинематические параметры в большей степени связаны со структурным взаиморасположением различных областей гетерогенности, а динамические - с изменением вещественного состава, физических свойств и напряженного состояния вещества внутри этих областей. Динамическая обработка разреза ОГТ с восстановлением истинного соотношения амплитуд позволяет перейти к разделению сейсмических разрезов на основе кинематико-динамических параметров.

Среди существующих в настоящее время методик обработки и интерпретации глубинного ОГТ как основные можно выделить следующие.

Комплекс программ OKS (О.К.Кондратьев, [31, 32]) позволяет осуществить интерактивную обработку сложных сейсмических материалов с соблюдением следующих принципов: объектная ориентированность обработки с возможностью получения и проверки разноплановых вариантов разреза и количественной характеристики их достоверности; обеспечение корректности выполняемых процедур при минимальных искажениях полезной части волнового поля и исключении ложной информации путем использования в программах блоков анализа ситуации и принятия логических решений.

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

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

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

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

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

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

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

В этом смысле представляет интерес использование при интерпретации сейсмических материалов компьютерной технологии КОСКАД 3D [48 - 52], включающей реализацию широкого спектра оригинальных методов обработки геолого-геофизической информации.

Технология КОСКАД 3D находит все большее применение при обработке сейсмической параметризованной информации алгоритмами статистического и спектрально-корреляционного анализа, обнаружения слабых аномалий, распознавания образов и классификации. Применение алгоритмов программного комплекса в обработке данных сейсморазведки стало возможным благодаря появлению высокопроизводительной вычислительной техники.

Адаптации современных статистических алгоритмов технологии КОСКАД 3D в обработке сейсмической информации послужило то, что большинство из них ориентировано на анализ разнородной геолого-геофизической информации, т. е. информации, не зависящей от природы анализируемых полей. Такие сейсмические параметры, как мгновенная амплитуда, частота, коп герентность и фаза, хотя и обладают специфическими особенностями, но вполне могут быть предметом анализа методами вероятностно-статистического подхода, реализованными в технологии КОСКАД 3D.

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

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

Наиболее распространенным предположением относительно свойств наблюдаемых геополей при построении линейных оптимальных фильтров является предположение о стационарности спектрально-корреляционных характеристик полезного сигнала и помехи. Реально наблюдаемые геополя, в том числе данные 2D- и ЗВ-сейсморазведки, довольно редко удовлетворяют этому жесткому требованию. Таким образом, применение линейных оптимальных фильтров не всегда приводит к качественным конечным результатам без предварительной обработки исходной информации.

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

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

Настоящее исследование посвящено созданию специальных приемов обработки глубинных сейсмических данных по региональным и опорным профилям, позволяющих изучать динамические свойства волновых полей, вычисляемые в рамках в рамках характерного размера локальных неодно-родностей гетерогенной среды. Динамические характеристики неоднородных объектов являются более чувствительными параметрами в гетерогенных средах, чем кинематические. Кинематические параметры определяются расстоянием и скоростью распространения волн в среде и используются, в основном, для определения скоростной характеристики толщ пород и пространственного положения изучаемых объектов. Значения динамических параметров упругих волн зависят от неоднородности элементарного объема геосреды, где формируется упругая или акустическая волна. Важным является также и характер неоднородности, к числу которых можно отнести неравномерное распределение слоистости, трещиноватости, пористости, насыщенности и т.п. Именно для определения этих свойств гетерогенной среды используются динамические параметры.

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

Например, красноярскими геофизиками развивается способ сейсмической разведки объектов, рассеивающих упругие волны, включающий возбуждение сейсмического сигнала, регистрацию сейсмического поля и обработку данных, направленную на выделение отраженных и рассеянных волн с помощью моделирования волнового поля от ожидаемых объектов рассеяния [Кочнев В.А., Поляков B.C., 33]. В данном способе с помощью моделирования подбираются такие параметры обработки сейсмических данных, которые позволяют наилучшим образом выделить рассеивающие объекты. Известно, что отраженные, дифрагированные и рассеянные волны, формирующиеся на границах неоднородных объектов в земной коре, дают общее интерференционное поле отраженных волн. В указанном выше способе на основе гипотезы о расположении ожидаемых объектов рассеяния проводят направленное выделение как отраженных, так и рассеянных волн путем моделирования волнового поля. Контроль результата и уточнение модели производится по обобщенным зависимостям амплитуд спектров сейсмических трасс от времени и частоты. Аномалии спектров сейсмических трасс используют как настройку для выделения объектов, рассеивающих упругие волны.

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

Известен также способ обработки сейсмических данных, позволяющий добиться улучшения количественной оценки и визуализации отражения от тонкого пласта и других боковых нарушений сплошности горной породы, основанный на выявлении специальных эффектов в амплитудных спектрах участков сейсмических трасс, заключающих в себе указанные геологические объекты [Патрика Г.А., Гридлей Д.М, 48]. Авторами этого способа замечено, что отражение от тонкого пласта имеет характерную форму выражения в частотной области. Его амплитудно-частотный спектр содержит периодическую последовательность отметок (всплесков), отстоящих друг от друга на расстояние, обратно пропорциональное «временной толщине» тонкого пласта. Это характерное выражение используется для отслеживания отражений тонкого пласта, а также для оценки его толщины и протяженности.

Способ включает в себя получение отображения набора сейсмических трасс, распределённых по заданному объёму земли, выбор интересующей части указанного объема, преобразование выбранных фрагментов сейсмических трасс с использованием дискретного ортонормального преобразования (например, Фурье, Уолша), которое даёт множество коэффициентов преобразования (вычисление амплитудного и фазового спектров фрагментов трасс). Полученные коэффициенты ортонормального преобразования организуют в куб настройки, который используют для идентификации и локализации углеводородов. Известно, что амплитудный спектр всей сейсмической трассы близок к спектру импульса источника возбуждения упругих волн, вследствие случайного характера последовательности отражений, возникающей в геологической среде. Однако спектры, вычисленные для более коротких фрагментов сейсмических трасс, уже не являются стационарными в указанном выше смысле и в них отражены свойства геологической среды, заключенной в пределах изучаемого фрагмента. Тонкий пласт действует как фильтр в частотной области и вносит собственный вклад в спектр импульса источника возбуждения: возникают узкие полосы ослабленных частот, имеющие характерный внешний вид, а расстояние между ними равно обратной «временной толщине» тонкого пласта. Таким образом, ослабленные частоты в амплитудном спектре используются в изобретении-прототипе для идентификации отражения от тонкого пласта и для измерения его толщины. Достижением известного способа является переход к спектральному анализу целевого участка сейсмической трассы, а не всей трассы в целом. Это позволяет добиться большей детальности при изучении свойств конкретных геологических объектов. Авторам рассматриваемого способа удалось выявить и использовать для идентификации спектральный образ тонкого пласта и близких к нему ситуаций выклинивания, бокового нарушения спошности горной породы. Однако, геологические объекты чрезвычайно многообразны и столь же многообразны их проявления в амплитудных и фазовых спектрах. Способ - прототип направлен главным образом на улучшение количественной оценки и на настройку визуализации для обнаружения эффектов от тонкого сейсмического пласта и других малых геологических объектов: линз, областей выклинивания пластов, боковых нарушений сплошности пород, являющиеся актуальными при поиске и локализации неантиклинальных ловушек углеводородов в осадочных толщах, что ограничивает область его применения.

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

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

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

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

Разработанный автором специализированный программный комплекс StreamSDS [Гошко Е.Ю., 13 - 15 ] позволяет исследовать один из наиболее важных динамических эффектов с точки зрения прогноза месторождений полезных ископаемых в условиях сложных гетерогенных типов сред, мало различающихся по скоростям, - поглощение энергии отраженных волн, проявляющееся в нелинейности коэффициента затухании амплитуд отраженных волн в области средних и высоких частот. Эта характеристика гетерогенной среды связана с важными параметрами прогноза полезных ископаемых на больших глубинах: повышенной трещиноватостью, разломами и разрушениями, а также заполнением порового пространства флюидами различного происхождения.

Глубинные разрезы ОГТ по опорным и региональным профилям содержат наиболее полную информацию о строении и состоянии пород в глубоко залегающих зонах земной коры и мантии. В настоящее время изучение глубинного строения Восточной территории России для развития ее минерально-сырьевой базы стало приоритетной задачей Министерства Природных Ресурсов Российской Федерации. С 1994 г. и по настоящее время в России создается государственная сеть опорных геофизических профилей для обеспечения широкого круга недропользователей информацией о строении недр. В настоящее время автор данной работы принимает участие в глубинных исследованиях, проводящихся на опорных профилях: 1-СБ, 3-СБ, 2-ДВ. В диссертации представлены результаты развития автором алгоритмов и вычислительных технологий обработки и интерпретации глубинных сейсмических разрезов в условиях гетерогенных сред земной коры.

Применение предложенной вычислительной технологии на материалах опорных сейсмических профилей ОГТ дало результаты, хорошо согласующиеся с прежними представлениями о строении исследуемых территорий, но и добавило новые, интересные с геологической точки зрения, данные о строении земной коры [Гошко Е.Ю.,16 - 18, Мигурский А.В., Гошко Е.Ю., 38 -40].

Основные результаты диссертационной работы

1. Разработан алгоритм вычисления локальной формы отраженного волнового пакета по данным профильных глубинных сейсмических наблюдений МОВ-ОГТ. Он позволяет изучать изменяющуюся форму сейсмического отражения в сложных гетерогенных средах.

2. Проведено численное исследование алгоритма выделения формы сейсмического импульса на основе компьютерного моделирования волновых образов геологических сред с помощью синтетического разреза ОГТ. Применение алгоритма выделения локальных волновых пакетов на синтетическом разрезе ОГТ позволило восстановить заданную форму сейсмических импульсов и их частотные характеристики. Исследована устойчивость алгоритма, показано, что достаточно около ста суммирований в пределах окна вычислений для получения стабильной формы локального отражения.

3. Разработана методика изучения динамических (спектральных, энергетических) характеристик локальных волновых пакетов в условиях сложных гетерогенных сред земной коры.

4. Создан вычислительный программный комплекс StreamSDS, позволяющий: создавать и визуализировать матрицу локальных волновых пакетов, вычисленную на глубинном разрезе ОГТ с помощью вышеназванного алгоритма, позволяющую анализировать динамические характеристики отражений от локальных неоднородностей среды; формировать куб амплитудных спектров локальных волновых пакетов в пространстве трех переменных: временной и латеральной координат сейсмического профиля и частоты (t, х, f); получать и визуализировать двумерные срезы спектрального куба по любой из вышеперечисленных координат: t, х, f для анализа и разделения глубинных сейсмических данных на области, однородные по спектральным и энергетическим характеристикам локальных волновых пакетов, отраженных от неоднородностей геологической среды; диагностировать зоны затухания энергии отраженных волн на основе исследования срезов спектрального куба по частоте, позволяющего обнаружить эффекты, связанные с поглощением высоких частот, рассеянием и другими потерями энергии при прохождении сейсмического импульса в неоднородной толще пород.

5. Разработанное автором алгоритмическое и компьютерное обеспечение для исследования динамических параметров отраженных волн было внедрено во ФГУП «СНИИГГиМС» в соответствии с Актом о внедрении (см. Приложение 1). Вычислительный комплекс StreamSDS был применен для обработки и интерпретации материалов по опорным сейсмическим профилям Восточной Сибири и Дальнего Востока в рамках выполнения договоров ФГУП «СНИИГГиМС» с Министерством Природных ресурсов, Агентством по недропользованию РФ: опорному профилю 2-ДВ (г. Магадан - мыс Эммытаген), 2000 погонных километров, длительность регистрации наблюдений 25 сек; рассечки к профилю 3-СБ «Алтай-Северная Земля», 300 погонных километров, длительность регистрации наблюдений 24 сек; рассечки к опорному профилю 1-СБ «Батолит», 50 пог. км, длительность регистрации наблюдений 18 сек.

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

Заключение

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

Численные расчеты, проведенные с помощью описанного алгоритма и разработанного на его основе специализированного вычислительного комплекса StreamSDS на фрагменте опорного сейсмического профиля ОГТ, показывают адекватные (с точки зрения геологической интерпретации) результаты. Предлагаемый алгоритм, его программная реализация, методика анализа и интерпретации амплитудных полей высокочастотной части спектра отраженного сейсмического сигнала являются новым инструментом исследования нелинейных эффектов поглощения высоких частот, связанных со свойствами среды. Особенно это актуально в глубинных зонах, не доступных бурению.

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

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

Исследование сейсмоакустического разреза, получаемого в системе РеапакРД [Рудницкая Д.И., 2], с помощью сейсмоакустических атрибутов значительно увеличивает обоснованность и достоверность геолого-геофизическиз моделей изучаемых сред земной коры. Программы блока ATRIBUTES [Гошко Е.Ю., 15] были успешно применены для вычисления физических и геометрических характеристик сейсмоакустических разрезов по геотраверсам 2-ДВ и 3-СБ.

1. Авербух А.Г. Изучение состава и свойств пород при сейсморазведке. - М.: Недра, 1982, - 232 с.

2. Авторское свидетельство № 990128 «свидетельство об официальной регистрации программы для ЭВМ». Рудницкая Д.И., Горячев Д.Н., Корнилов М.В., Гошко Е.Ю. и др. - Регистрация 12.03.99 г.

3. Аки К., Ричарде П. Количественная сейсмология. Т. 1, 2. М.: Мир, 1983.

4. Алексеев А.С., Цибульчик Г.М. Математические модели сейсморазведки. В сб. «Актуальные проблемы вычислительной математики и математического моделирования» Отв. ред. Алексеев А.С. Новосибирск: Наука, 1985 г., с. 91-108.

5. Алексеев А.С. О проблемах развития информационно-вычислительных технологий в геологоразведке. Доклады научно-практической конференции «Инновационные технологии для нефтегазового комплекса», г. Новосибирск, 2005 г.

6. Берзон И.С., Епинатьева A.M., Парийская Г'.Н., Стародубровская

7. С.П. Динамические характеристики сейсмических волн в реальных средах.122

8. Изд. АН СССР, Москва, 1962.

9. Биезайс Я.Я., Рудницкая Д.И., Романов Н.Я, 2000 Строение земной коры в районе Мирнинского кимберлитового поля по данным СГ-ОГТ: Геофизика,№1, 40-42.

10. Васильев В.П., Рудницкая Д.И., Сагайдачная О.М., Горячев Д.Н. Результаты интерпретации данных СГ-ОГТ методами системы РЕАПАК. Геофизика, № 3, 1999. 10 13.

11. Ведерников Г.В., Ермолаева Г.М., Петров А.В., Поздняков В.А., Рудницкая Д.И., Сергеев В.Л. Экспертная обработка материалов СГ-ОГТ. -Геофизика, № 5, 2003 г.

12. Гошко Е.Ю., Зеркаль С.М. Алгоритмическое обеспечение исследования гетерогенных сред земной коры по динамическим характеристикам локальных волновых пакетов. Вестник НГУ, Том 4, Выпуск 1, Серия: Информационные технологии, 2006 г.

13. Гошко Е.Ю., Рудницкая Д.И., Сагайдачная О.М., Сальников А.С.

14. Гутаров А.П., Жданович В.В., Михальцев А.В., Потапов О.А. Возможности и ограничения динамической обработки сейсморазведочных данных. М., 1985, 45 е., - Разведочная геофизика: Обзор ВНИИ экон. Минер. Сырья и геологоразвед. Работ. (ВИЭМС)

15. Зеркаль С.М. Определение непрозрачных зон в Земле методом компьютерной томографии в кинематической постановке. ДАН СССР. 1991. Т. 317, №2, с. 330-333.

16. Зеркаль С.М. Численное решение обратной трехмерной кинематической задачи сейсмики в линеаризованной постановке. Геология и геофизика. 1988. № 11, с. 126- 133.

17. Жерняк Г. Ф., Сергеев В. JI. и др. Система высокоразрешающей обработки, анализа и моделирования сейсмических данных: Геофизика, спецвыпуск Сибнефтегеофизика, 2001, с. 53 -55.

18. Караев Н. А., Анисимов А. А., Кашкевич В. И., Травинская Т. К. Сейсмическая гетерогенность земной коры и ее отображение в поле рассеянных волн: Геофизика, 2,1998, с. 29-39.

19. Караев Н. А., Скотт Смитсон, 1995, Многокомпонентные сейсмические исследования в районе Кольской сверхглубокой скважины: Геофизика, 1, 1995.32 -40.

20. Караев Н.А., Лебедкин П.А. Сейсмическая гетерогенность земной коры и модель волнового поля. Исследование литосферы в работах Петербургских геофизиков. СПб.: ВИРГ - «Рудгеофизика», ВНИИОкеанология, 2003 г.

21. Козлов Е.А., Гогоненков Г.Н., Лернер Б.Я. и др. Цифровая обработка сейсмических данных. - М.: Недра, 1973, - 312 с.

22. Кондратьев И.К., Киселев Ю.А., Ковригина Е.И., Крылов Д.Н.

23. Оценка точности и разрешающей способности используемых на практике способов решения обратной динамической задачи сейсморазведки. Региональная и морская геофизика: Обзор ВИЭМС. - М., 1990.

24. Кондратьев О. К. Идеология и средства обработки сейсмических материалов: Геофизика, 5, 1998, с. 3 -11.

25. Кондратьев О. К., Говорова Л. А., Гуляева М. М., Скоростная параметризация земной коры по отраженным волнам: Глубинное строение и геодинамика Южного Урала. Проект "Уралсейс" Спецгеофизика: Тверь, 2001, с. 62-70

26. Кочнев В.А., Поляков B.C. и др, патент РФ № 2 248 014, G 01 V 1/00, Красноярск, 2004 г.

27. Кузнецов O.JL, Чиркин И.А., Курьянов Ю.А. и др. Сейсмоакусти-ка пористых и трещиноватых геологических сред. Экспериментальные исследования. В трех томах. - М.: Государственный научный центр РФ -ВНИИГеосистем, 2004.

28. Кунин Н.Я., Иогансон Л.И. Геофизическая характеристика и строение земной коры Западной Сибири. М.: ИФЗ АН СССР. 1984. 218 с.

29. Кунин Н.Я., Шейх-Заде Э.Р. Исследования литосферы докрити-ческими отраженными волнами. М.: Наука, 1993.

30. Максимов С.П., М.И. Лоджевская, Р.Г. Самвелов, Б.А. Соловьев, Б.Д. Гончарешсо. Геологические условия нефтегазоносности на больших глубинах. Международный геологический конгресс, XXVIII сессия. Докл.сов. геологов. -М., 1988.

31. Мигурский А.В., Гошко Е.Ю. «Сейсмогеологическая модель земной коры по геотраверсу Уралсейс». Монография: «Результаты исследований по международной программе «Уралсейс». Тверь, 2002 г.

32. Михалъцев А.В., Вилкова Э.С., Давыдова JI.H. и др. Влияние не-однородностей верхней части разреза на оценки амплитуд отраженных волн. Прикладная геофизика, вып. 105. М.: Недра, 1982, с. 19-33.

33. Михальцев А.В, Киричек М.А., Петухов А.В. Прямое прогнозирование месторождений нефти и газа с применением комплекса геофизических и геохимических методов. Доклады геологического конгресса в Москве, Секция «Геофизика», Т. 8, М.: Наука, 1984, с. 134 153.

34. Михальцев А.В., Мушин И.А., Погожев В.М. Алгоритмы и графыцифровой обработки амплитуд отражений в структурно-формационной сейсморазведке. Прикладная геофизика, 1980, вып. 97, с. 24-44.

35. Мячкин В. И., 1978, Процессы подготовки землетрясений: Наука.

36. Никитин А. А., 1979, Статистические методы выделения геофизических аномалий: ML, Недра.

37. Павленкова Н.И. Развитие представлений о сейсмических моделях земной коры. Геофизшса, № 4, 1996. ,11 18.

38. Патрика Г.А., Гриддей Д.М. (США, БП Корпорейшн Норт Америка инк.), патент РФ № 2 187 828, G 01 V 1/30,1998 г

39. Петров А. В., 1996, Адаптивная фильтрация геополей: Геоинформатика,

40. Петров А. В., 1996, Методы адаптивной фильтрации и многомерного дисперсионного анализа в компьютерной технологии КОСКАД 3D: Автореферат диссертации на соискание ученой степени

41. Петров А. В., Никитин А. А. Классификация комплексных геополей на однородные области. Геология и разведка. Изв. вузов, № 3,1999

42. Поздняков В, А., Ледяев А. И., 2001, Программная реализация сейсмоголографического фокусирующего преобразованияв РгоМАХ: Новые технологии в геофизике: Уфа, МинэнергоРФ 247 248.

43. Поздняков В. А., Сафонов Д. В., Чеверда В. А., 2000, Оптимизация параметров фокусирующих преобразований с использованием численного моделирования: Геология и геофизика, 41, 6, 930 938.

44. Применение цифровой сейсморазведки для прямых поисков нефтегазовых залежей/ А.Г.Авербух, В.А.Гельфанд, Г.Н.Гогоненков и др. Обзор. М.:ВИЭМС, 1979. - 69с.

45. Пузырев Н.Н. Методы и объекты сейсмических исследований (введение в общую сейсмологию). Новосибирск: Изд-во СО РАН, НИЦ ОИГГМ, 1997 г.

46. Результаты интерпретации данных СГ-ОГТ методами системы РЕАПАК / В. П. Васильев, Д. И. Рудницкая, О. М. Сагайдачная, Д. Н. Горячев // Геофизика. М.: ГЕРС, 1999. - № 3.

47. Рудницкая Д.И., Биезайс Я.Я., Гошко Е.Ю. «Применение системы130

48. Рудницкая Д.И., Корнилов М.В., Гошко Е.Ю. «Опыт использования системы РЕАПАК при изучении месторождений нефти и газа в Западной Сибири». Научно-технический журнал ЕАГО, Геофизика, 1996 г., № 3, с. 19-24.

49. Скрипий А.А., Структура глубинного оползания Южного Урала по данным «Уралсейс-95»: Тезисы совещания «Современные проблемы шарьяжно-надвиговой тектоники»: Уфа. 1997

50. Соборнов К.О., 1996, Структурные сегменты уральского пояса надвигов и его нефтегазоносность: Геология нефти и газа, 3.

51. Солодилов JI. Н., Попова О. Г., Коновалов Ю. Ф., Кухмазов С. У., Минина Н. А., Кадурина JI. С., 1999, Возможность мониторинга состоянияземной коры с использованием далеких землетрясений: Разведка и охрана недр, 11, 51 53.

52. Страхов В. Н., 1989, К новой парадигме сейсмологии: Природа, 12, 4-9. Кондратьев О. К., 2000, Новая парадигма прогноза землетрясений: Наука в России, 4, 45 51.

53. Суворов В. Д. Глубинные сейсмические исследования в Якутской кимберлитовой провинции. Новосибирск: Наука, 1993.

54. Янг Р. П., Хилл Дж. Дж. Статистический анализ спектров сейсмических записей с целью определения свойств горных пород. В сб.: Анализ и выделение сейсмических сигнлов. -М.: Мир, 1986.

55. Benin R., Oncken О., Knapp J. Н., Perez-Estaun A., Hismatulin T.,Yunusov N., Lipilin A., 1996, Orogenic evolution of the UralMountains: results from an integrated seismic experiment: Science, 274/5285, 220 221.

56. Berzin R., Oncken 0., Knapp J.H., Perez-Estaun A.,Hismatulin Т., Yunusov N.,Lipilin A., 1996, Orogenic evolution of the Ural Mountains: results from an integrated seismic experiment: Science, 274/5285, pp. 220—221.

57. Carbonell R., Perez-Estaun A., Gallar J., Dia J., Kashubin S., Mechie J., Stadtlander R., Schulze A., Knapp J.H., Morozov A., 1996, Crustal root beneath the Urals: wide-angle seismic evidence: Science, 274/5285, pp. 222—224.

58. Carbonell R., Perez-Estaun A., Gallar J., Dia J., Kashubin S., Mechie J., Stadtlander R., Schulze A., Knapp J. H., Morozov A., 1996, Crustal root beneath the Urals: wide-angle seismic evidence: Science, 274/5285, 222 224.

59. Chorpa S., Marfurt K.J. Seismic attributes A historical perspective. Geophysics, Vol. 70, № 5, 2005.

60. Hurich C. A., 1996, Statistical description of seismic reflection wave-fields: a step towards quantitative interpretation of deep seismic reflection profiles: Geophys.J.Int, 125, pp. 719 728.

61. Knopoff L.Q. Review of Geophysics, 2,625 660,1964

62. Solomon S.C. Seismic-wave attenuation and partial melting in the upper mantle of North America. Journal of Geophysical Research. 1972, 78, 6044 -6059.

63. Spenser T.W. The method of generalized reflection and transmission coefficient. Geophysics, 25, 625 641, 1960.