Анализ информационных свойств взаимных фазовых спектров отраженных сейсмических волн Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

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

При завершении использования объектов определенного типа работа специализированного агента мониторинга завершается.

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

Таким образом, агент мониторинга выступает сервисом, который заказывает агент координатор. На схеме также представлены блоки [3]:

• контейнер - агентная среда на конечном рабочем месте в автоматизированной информационной системе;

• платформа - группа контейнеров, объединенная одним блоком управления, реализующим агентную парадигму.

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

СПИСОК ЛИТЕРАТУРЫ

1. Ткаченко Н.И., Спирин Н.А. Применение сервис-ориентированной архитектуры при интеграции систем управления технологическими процессами // Известия Томского политехнического университета. - 2010. - Т. 317. - № 5. - С. 61-67.

2. Новиков Д.А. Математические модели формирования и функционирования команд. - М.: Изд-во физико-математической литературы, 2008. - 184 с.

3. Bellifemine F., Caire G. Developing Multi-Agent Systems with JADE. - Chichester: John Wiley & Sons Ltd., 2007. - 286 p.

Поступила 21.09.2012 г.

УДК 550.8.053

АНАЛИЗ ИНФОРМАЦИОННЫХ СВОЙСТВ ВЗАИМНЫХ ФАЗОВЫХ СПЕКТРОВ ОТРАЖЕННЫХ

СЕЙСМИЧЕСКИХ ВОЛН

В.П. Иванченков, А.И. Кочегуров, М.А. Черкасова*

Томский политехнический университет *ЗАО «Гринатом», г. Северск E-mail: kai@cc.tpu.edu.ru

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

Ключевые слова:

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

Key words:

Linear system, model of layered absorbing column, mutual phase spectrum of signals, reflection and wave refraction factors.

При решении задач прогноза геологического разреза (ПГР), в том числе прогноза нефтегазонос-ности осадочных толщ, по данным сейсмических наблюдений наиболее широко используются в качестве диагностических признаков динамические характеристики отраженных волн, непосредственно связанные с их амплитудой и энергией [1, 2]. Информация о свойствах фазочастотных характеристик (ФЧХ) сейсмических волн до последнего

времени практически не использовалась. Между тем в фазу сейсмических сигналов, а точнее в сложный закон изменения их фазовых спектров, заложена информация, позволяющая в условиях априорной неопределенности надежно обнаруживать и разрешать сигналы на фоне интенсивных помех, производить оценку их кинематических параметров [3, 4]. Как показано в [5, 6], текущие фазовые спектры отраженных сейсмических волн могут

1DS

Рис. 1. Модель плоско-параллельной слоистой поглощающей толщи

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

В рамках предложенного ранее системного подхода [6] предварительно остановимся нарассмо-трении простой модели слоистой поглощающей толщи с горизонтальными границами раздела выделяемых слоев (рис 1, а). В лучевом приближении при нормальном падении волны 50(0 процесс формирования отраженных волн от кровли 51(0 и подошвы 52(0 линейно-неупругого поглощающего слоя II можно представить в виде эквивалентной системы, показанной на рис. 1, б.

При этом в отсутствии помех взаимный спектр (ВСП) для сигналов 51(0 и 52(0 можно определить:

й,(/) = [ 5 ’,(/) • £,(/)] = |е,2( / )| ' >,

где ф12(/) - взаимный фазовый спектр сигналов 51(0 и 52(о. 51(0=^51(0}, 52(0=^52(0} - спектры волн, отраженных от кровли и подошвы 11-го слоя, ^ -оператор прямого преобразования Фурье. В соответствии с рис. 1, б, спектры 51(/) и 52(/) запишутся:

5,(/) = 50(/)ки(/) = |5,(/)| е^</

5 2 (/) = 5о(/)ТП(/)Н (/)к 23 (/)г21(Л =

= |^2 (/)| ^>, где 50(0=^50(0}, к(/), г(/) - коэффициенты отражения и преломления для границ раздела соответствующих слоев, Н(/) - частотная характеристика системы, определяющая свойства поглощающего слоя II.

Тогда взаимный спектр 012(/) можно представить:

е^/) = |5о(/)Г /) *23 (/) Н /);

|а2(/)| = |5о(/)Г • М/)|• М/)|х

х| Н (/)| •• \га(/)| • |гм( / )|;

Фп(„ () =

= ф*23(/)-Ф (/) + Фн (/) + ф„(/) +1, (/)- (1)

где фк(/}, фг(/} - аргументы комплексных коэффициентов отражения и преломления, фд{/) - фазочастотная характеристика частной системы, определяющей распространение волны в поглощающем слое II в двух направлениях.

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

Важно отметить, что ВФС волн, отраженных от кровли и подошвы исследуемой толщи, не зависит от фазового спектра падающей волны 50(0, что является одним из положительных факторов, который позволяет повысить надежность дальнейшего прогноза свойств изучаемого геологического разреза.

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

Н (/) = г12(/) Н (/)к2з(/)г21(/) + к12(/) =

,.М: (Р)

|Н £ (/ )|,

Рис. 2. Два изотропных полупространства, разделенных плоской границей

Согласно [6, 7], для двух изотропных полупространств, разделенных плоской границей, аргументы коэффициентов отражения и преломления запишутся как (рис. 2):

4пуи(/)(52(/)-5Х(/))

фк = arctg

4п2(1 -г12(/)) + ф2( /) -52( /) Гі22.( /))

; (2)

фг = аг^

2п(8х( /) -82( /))

4п +5! (/ШГ) + у12(/)(4п +522(/))

, (3)

где ъ(/)=рУ(/), (г'=1,2); Ы/) - - отноше-

ние акустических жесткостей контактирующих сред, а V/ и 8/ соответственно определяют скорость распространения продольных волн и декремент поглощения в среде на частоте / При этом

а(/ )У (/)

5(/) = -

/

где а(/) - коэффициент поглощения.

По результатам различных исследований часто принимается, что коэффициент поглощения а(/) связан с частотой линейной зависимостью [1]:

а( /) = в/,

где в - параметр поглощения.

Исходя из (1), взаимный фазовый спектр отраженных волн зависит от ФЧХ системы, определяющей свойства поглощающего слоя II. Принимая для поглощающего слоя минимально-фазовую модель среды и учитывая ранее сделанное допущение о линейной зависимости коэффициента поглощения от частоты /, выражения для частотной характеристики Н(/) можно представить как [6]:

Н (/) = |Н (/)| е-УН (/), |н (/)| = е-2а( /) А,

1

Ун (/) = 4п/

V ( /о)

(4)

наблюдаются даже при очень малом изменении (порядка 0,05) декремента поглощения. Это свидетельствует о высокой чувствительности (рк/ к небольшим изменениям поглощающих свойств отложений. Фазовые набеги за счет преломления на границах раздела (рг/ в поглощающих средах существенно меньше, чем за счет отражения. С ростом частоты при любых сочетаниях у12(/0) и 512(/0) аргумент коэффициента преломления уменьшается и имеет отрицательную величину.

где V (/0) - фазовая скорость распространения волны на некоторой базовой частоте /0, к - мощность слоя.

Из соотношений (1-4) видно, что взаимный фазовый спектр волн, отраженных от кровли и подошвы анализируемой толщи, непосредственно зависит от акустических и поглощающих свойств контактирующих толщ. На рис. 3 в качестве примера приведены графики аргументов коэффициентов отражения и преломления волн, определенных при различных значениях отношений декрементов поглощения 8П/ и акустических жесткостей у12/ контактирующих сред. Исследования показали, что аргумент коэффициента отражения во всех случаях, когда декременты поглощения контактирующих сред не равны между собой, зависит от частоты. Влияние поглощения на аргумент коэффициента отражения тем значимее, чем меньше дифференциация акустических жесткостей контактирующих сред. Так, при у12(/0)=0,9 и 512/0)=0,05-1

/

изменение фазы по частоте — - 2 составляет

/о

1,5 рад., а в случае 712/0)=0,62 при тех же самых условиях значение уменьшается до 1 рад.

При этом следует отметить, что в этом случае достаточно большие фазовые сдвиги («0,25 рад.)

У12=0,62 г Уп=0,9

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

По аналогии с изложенным выше может быть построена модель поглощающей толщи, состоящей из произвольного числа слоев и определена частотная характеристика Н/ ее эквивалентной системы [6]. На основе определенной Н/ может быть рассчитана синтетическая сейсмотрасса для любой точки поверхности наблюдения:

^-1{5о( /) Н (/)},

где Т- - оператор обратного преобразования Фурье.

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

Предварительно исследования свойств ВФС проводились для 4-х моделей слоистых поглощающих сред. Общая структура принятых моделей приведена на рис. 4.

Аргиллиты

Продуктивная толща

Аргиллиты

Рис. 4. Структура моделей слоистых поглощающих сред

Каждая модель включает три толщи, где вторая толща - песчаный коллектор (продуктивная тол-

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

У 12 = °ДГ23 = 0,8

уХ1 = 0,96; у = 0,825

0 0,5 1 1,5

!:£„ = 0,167; &3 = 5

Ф.2(/)

К1— 3

2

^ 3 2

х 4 1 0

ЧУ о

0 0,5 1 1,5

3:812 = 0,038; = 11,667

Т.8и = 0,064; £23 = 13,33; 4: £12 = 0,026; = 31,3

а Модель I

/12 = 0,897; /23 = 0,929 /12 = 0,922; /13 = 0,966

1:512 = 0,316; 81Ъ = 2,633; 4: 8Х1 = 0,162; 81Ъ = 5,133 2:812 = 0,24; 81Ъ = 3,467; 5: 811 = 0,14; 81Ъ = 5,967 3 : 812 = 0,194; 523 = 4,3; 6: 812 = 0,117;523 = 7,133

б Модель II

уп = 0,765;^23 =1,4

ф,:(/)

уи =0,771; у23 =1,5 ф ■=(/>

о

0,5

-1

-0,5

Л.

3

V \ 4

\6_

0 0,5 1 1,5

= 0,266; &3 = 3,133

Ис

1: 81г = 0,391; 8гъ = 2,133; 4: 81 2:81г = 0,338; 8г3 = 2,467; 5: 81г = 0,24; 8гг = 3,467 3: 81г = 0,298; 523 = 2,8; 6:512 = 0,21;523 = 3,967

в Модель III

у 12 = 0,5;г23 =1,5

ф,=(/)

3

2

1

0

А Л

1 / г~

у 1П =0,513;^,, =1,58 ф,:(/)

0 0,5 1 1,5

1:812 = 0,962; 8,

1,5 ,3,4,5 ,6

/

и

0,5 1 1,5

= 0,446; 523 = 1,867 = 2,2 = 2,433

Ио

При исследованиях скорости распространения продольных волн, плотности пород и параметры поглощения во всех слоях выбирались исходя из обобщенных данных, характерных для верхнеюрской продуктивной толщи Крапивинской площади Томского Приобья. Мощность (толщина) второго слоя задавалась 15 м. На рис. 5 показаны графики взаимных фазовых спектров волн, отраженных от кровли и подошвы П-го слоя, рассчитанных для различных отношений декрементов поглощения 8п(/), 823(/ и акустических жесткостей 7п(/), УпО) контактирующих сред.

Анализ графиков показывает, что насыщение породы-коллектора газом или нефтью непосредственно приводит к значимому изменению взаимных фазовых спектров (модель I и II на рис. 5). При этом наибольшие отклонения взаимного фазового спектра при определенных сочетаниях отношений акустических жесткостей и декрементов поглощения контактирующих сред могут достигать 1,8 рад. В случае насыщения песчаного коллектора водой (модель III на рис. 5) в диапазоне изменения отношения частот// от 0,5 до 2 ВФС волн слабо зависят от частоты и их значения сосредоточены в области от 0...-1 рад. При отсутствии коллектора (модель IV на рис. 5) взаимные фазовые спектры отраженных волн практически не зависят от частоты и их значения при различных сочетаниях петро-физических параметров контактирующих слоев близки к 0. Из приведенных зависимостей также видно, что графики ВФС в случае насыщения песчаного коллектора газом или нефтью во всем диапазоне анализируемых частот обладают большей кривизной и имеют больший интервал изменения значений, чем в случае отсутствия коллектора или его насыщения водой. Это позволяет при решении задач прогноза нефтегазоносности осадочных толщ ввести в качестве диагностических признаков параметры ВФС, учитывающие отклонение их значений от средних. В частности, для этих целей может быть использован центральный момент 2-го порядка ВФС:

Ё Фі2( /)-Ф2( / ))2

п -1

(5)

где ф12 (/) =—ЁФ12 (/) - среднее значение взаим-

п І=1

ного фазового спектра; п - число гармоник, учитываемых при определении взаимного фазового спектра отраженных волн.

Наряду с взаимным фазовым спектром волн информацию о дисперсионности сред и поглощающих свойствах исследуемых толщ также несут взаимные фазовые и групповые задержки:

(/)=П • /)=-1/■

Степень их отклонения может быть соответственно определена:

Ио

0,867; 4:й,2 2: 81г = 0,694; 523 = 1,2; 5: 82г = 0,379; 83Ъ

3:81г = 0,543; 83 = 1,533; 6:82г = 0,342; 83 г Модель IV

Рис. 5. Взаимные фазовые спектры волн, отраженных от кровли и подошвы II слоя: а) газонасыщенного песчаника, б) нефтенасыщенного песчаника, в) водонасыщенного песчаника, г) карбонизированного песчаника

2

I=1

=

п

Аргиллиты (Баженовская пачка)

Аргиллиты (Баженовская пачка)

песчаник песчаник

Карбонизированный | Газонасыщенный | Карбонизированный Водонасыщенный | Газонасыщенный | Водонасыщенный

Аргиллиты Аргиллиты

а Модель I б Модель III

Рис. 6. Структура моделей слоистых поглощающих сред при изменении их свойств по латерали

п -1

м г, ) - и2

п -1

(6)

(7)

С целью исследования информативности параметров (5-7) дополнительно были рассмотрены 4 модели слоистых сред, состоящие из трех толщ. Покрывающая и подстилающие толщи являются, как и ранее, слоями аргиллитов, а второй слой задавался как продуктивная толща, свойства которой изменялись по латерали (рис 6, а, б).

Рассматривались следующие изменения по латерали структуры 11-го слоя:

1) Карбонизированный песчаник-газонасыщен-ный песчаник-карбонизированный песчаник (рис. 6, а);

2) Карбонизированный песчаник-нефтенасы-щенный песчаник-карбонизированный песчаник;

3) Водонасыщенный песчаник-газонасыщенный песчаник-водонасыщенный песчаник (рис. 6, б);

4) Водонасыщенный песчаник-нефтенасыщен-ный песчаник-водонасыщенный песчаник.

Петрофизические параметры слоев (скорости распространения волн, плотности пород, параметры поглощения) были выбраны из тех же условий, что и для рассмотренных ранее моделей сред. Мощность И-го слоя была задана равной 50 м. Во всех случаях при исследовании варьировались петрофи-зические параметры И-го слоя, при неизменных параметрах покрывающей и подстилающей толщ. Для каждой из принятых моделей по разработанной программе осуществлялся расчет и построение наблюдаемых волновых полей (рис. 7), на которых автоматически прослеживались волны, отраженные от кровли и подошвы П-го слоя.

Для них определялся их взаимный фазовый спектр (рис. 8) и оценивались прогнозные параметры (5-7).

На рис. 9 в качестве примера приведены графики изменений по латерали центральных моментов 2-го порядка фазовой и групповой задержек для моделей сред, показанных на рис. 6.

Из рис. 9 видно, что на участках насыщения песчаного коллектора газом наблюдаются аномальные изменения указанных прогнозных параметров. При этом, как следовало ожидать, при структуре И-го слоя типа «вода-газ-вода» (рис. 6, б) относительное изменение параметров на различных участках профиля оказалось меньшим, чем в слу-

2

і=1

* =

Рис. 8. Взаимные фазовые спектрыы волн, отраженных от кровли и подошвыы 11-го слоя: а) карбонизированный; б) газонасыы-щенный; в) водонасыщенный песчаник

2

ст і

0,2

0,15

ОД

0,05

0

Модель I

500 1x10 1.5x10

Модель I

а

2

0,25

0,2

0,15

б

Модель III

500 1x10 1,5x10

Модель III

Рис. 9. Центральный момент 2-го порядка: а) взаимной фазовой задержки; б) групповой задержки

чае, когда структура 11-го слоя по латерали была задана в виде, приведенном на рис. 6, а. Аналогичные закономерности в изменении указанных прогнозных параметров также наблюдались при насыщении породы коллектора нефтью (модели II и IV). Результаты проведенных исследований показали, что наиболее информативными параметрами являются центральные моменты взаимной фазовой и групповой задержек. Центральный момент взаимного фазового спектра при определенных сочетаниях петрофизических параметров пород продуктивного слоя имеет незначительные аномалии при задании структуры модели 11-го слоя по латерали в виде «вода-нефть-вода».

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

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

Рис. 10. а) фрагмент разреза ОГТ, профиль ПР 1; б) Фрагмент объектно-ориентированного фазового разреза ПР 1

Рис. 11. Изменение прогнозных параметров вдоль профиля наблюдения ПР 1: а) центрального момента 2-го порядка взаимной фазовой задержки; б) центрального момента 2-го порядка взаимной групповой задержки

а приведены фрагменты временных разрезов ОГТ, полученных после стандартной обработки, для двух сейсмических профилей (условное обозначение ПР 1 и ПР 2), на пересечении которых расположена скважина Р-1 (условное обозначение) с дебитом нефти 130 м3/сут. Следует отметить, что на временных разрезах ОГТ не всегда удается проследить оси синфазности отраженных волн, отождествляемых по стратиграфической привязке с продуктивным пластом Ю13-4 (рис. 10, а), что не позволяет провести вычисление оценок взаимных фазовых спектров целевых волн, привязанных к подошве Bg и пласту Ю13-4. В этой связи на основе ранее предположенных методов фазочастотного прослеживания (ФЧП) [4] была прове-

дена детальная обработка временных разрезов ОГТ и построены объектно-ориентированные разрезы, на которых на всем интервале выделяются целевые волны (рис. 10, б, 12, б). На основе данных разрезов была проведена оценка временного положения целевых волн и определение их взаимных фазовых спектров и указанных выше прогнозных параметров.

На рис. 11, 13 приведены графики принятых прогнозных параметров, определённых для указанных фрагментов разрезов ПР 1 и ПР 2. Из них видно, что на пикетах наблюдения обоих профилей, где расположена скважина Р-1, прогнозные параметры имеют ярко выраженные аномалии, свидетельствующие о наличии нефтяной залежи.

Р1

ПК

Вё

ю,

Рис. 12. а) фрагмент разреза ОГТ, профиль ПР 2; б) фрагмент объектно-ориентированного фазового разреза ПР 2

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

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

СПИСОК ЛИТЕРАТУРЫ

1. Гурвич И.И., Боганик Г.Н. Сейсморазведка. - Тверь: Изд-во АИС, 2006. - 774 с.

2. Птецов С.Н. Анализ волновых полей для прогнозирования геологического разреза. - М.: Недра, 1989. - 135 с.

3. Иванченков В.П., Кочегуров А.И. Определение временного положения сейсмических сигналов по оценкам их фазочастотных характеристик // Геология и геофизика. - 1988. - № 9. -С. 77-83.

4. Иванченков В.П., Кочегуров А.И., Орлов О.В. Исследование разрешающей способности методов фазочастотного прослеживания сейсмических сигналов // Известия Томского политехнического университета. - 2012. - Т. 320. - № 5. - С. 80-85.

5. Иванченков В.П., Козлов А.А. Классификация типов геологического разреза методами кластерного анализа поданным сейсмических наблюдений // Известия Томского политехнического университета. - 2011. - Т. 318. - № 1. - С. 54-59.

6. Иванченков В.П., Кочегуров А.И., Орлов О.В. Фазочастотные характеристики сейсмических волн и основные предпосылки их применения для решения задач прогноза геологического разреза // Известия Томского политехнического университета. - 2008. - Т. 313. - № 5. - С. 132-138.

7. Авербух А.Г, Трапезникова Н.А. Отражения и преломления плоских волн при нормальном падении на границы // Физика земли. - 1972. - №9. - С. 74-83.

Поступила 10.04.2013 г.