рефераты скачать

МЕНЮ


Математическое моделирование в сейсморазведке

Математическое моделирование в сейсморазведке

Лекция 1

ВВЕДЕНИЕ

МЕТОД МАТЕМАТИчЕСКОГО МОДЕЛИРОВАНИя ПОЛУчИЛ РАСПРОСТРАНЕНИЕ В

СЕЙСМОРАЗВЕДКЕ ПРИМЕРНО С СЕРЕДИНЫ 60-Х ГОДОВ. ИСПОЛЬЗОВАНИЕ СИНТЕТИчЕСКИХ

СЕЙСМОГРАММ (СС), яВЛяВШИХСя РЕЗУЛЬТАТОМ РЕШЕНИя ОДНОМЕРНОЙ ДИНАМИчЕСКОЙ

ЗАДАчИ, ИМЕЛО СЛЕДУЮЩИЙ ЦЕЛИ:

. анализ процесса формирования поля отраженных волн в тонкослоистой среде(

. оценку роли многократных волн в этом поле(

. определение сейсмических эффектов, обусловленных изменением литологии или

углеводородосодержания и др.

В целом это позволило получить важные для практики интерпретации

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

сейсмозаписи при изучении того или иного геологического объекта.

Переход к двумерному сейсмическому моделированию, т. е. к

использованию синтетических временных разрезов (СВР), означал не просто

увеличение количества синтезируемых трасс, а качественно новый уровень

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

применения математического моделирования непосредственно в процессе

интерпретации данных сейсмических наблюдений

К началу 80-х годов сложилась следующая классификация видов

сейсмомоделирования.

1. Структурное моделирование. Обычно такое моделирование осуществляется

путем прослеживания лучей, что позволяет воспроизвести истинный путь

сейсмических волн при пересечении границ напластований, установить точную

связь между временем и глубиной и понять причины своеобразного

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

структурного моделирования могут производиться оценка и учет влияния

вышележащих толщ на кинематику сейсмических волн в интересующем

(перспективном) диапазоне времен или глубин при решении стратиграфических

задач.

2. Стратиграфическое моделирование. Первоначально применялось с целью

получить оценку влияния залежей нефти и газа, выклинивающихся слоев, зон

литолого-фациального замещения и других неоднородностей на такие

характеристики сейсмической записи, как изменение амплитуд, искажение

вышезалегающих горизонтов, изменение полярности, понижение скорости,

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

менее важны, а упругие параметры горных пород, наоборот, являются очень

важными и решающими. В последнее время стратиграфическое моделирование

успешно применяется для обнаружения и подтверждения залежей

углеводородов, определения литологии пород, связанных с этими залежами,

границ распространения залежей и др.

3. Моделирование сейсмических скоростей. Вначале такое моделирование

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

отражающих и промежуточных границ и локальных неоднородностей на поле

сейсмических скоростей или, другими словами, для анализа отклонений

параметра VОГТ реальных скоростей в среде. Впоследствии были осуществлены

удачные опыты использования этого вида моделирования в качестве основы в

методах решения обратных кинематических задач для многослойных сред с

криволинейными границами раздела и с градиентами скоростей в слоях.

Кроме того, двумерное сейсмомоделирование стало эффективно

использоваться и на этапе обработки сейсмической информации для решения

таких задач, как:

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

верхней части разреза,

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

. синтез оптимальных графов обработки.

Общие принципы интерпретации данных сейсморазведки на основе

математического моделирования

1 Системный анализ проблемы интерпретации данных сейсмических наблюдений

В соответствии с методологическим принципом системного подхода

представим объект нашего изучения (процесс интерпретации данных

сейсмических наблюдений) в виде целостной системы взаимодействующих

элементов (верхняя часть рис. 1, а).

Будем называть интерпретацией данных сейсмических наблюдений процесс

построения сейсмогеологической модели, которая не противоречит имеющейся

априорной информации (наблюденному волновому полю, данным промысловой

геофизики, геологической информации) и опыту геофизика-интерпретатора. Из

этого определения следует несколько важных методологических выводов:

1) процесс интерпретации является целенаправленным и поэтому должен быть

управляемым;

2) в процессе интерпретации необходимо сопоставлять имеющуюся в данный

момент сейсмогеологическую модель с априорными данными (в первую очередь

с наблюденным волновым полем) на предмет анализа их противоречивости и

нахождения способов ее устранения;

3) ввиду невозможности непосредственного сопоставления таких разнородных

объектов, как сейсмогеологическая модель и наблюденное волновое поле, в

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

волновое поле по сейсмогеологической модели.

Таким образом, математическое моделирование становится неотъемлемой

частью технологии интерпретации.

Конкретизируя схему рис. 1, а, получаем схему интерпретации данных

сейсморазведки на основе математического моделирования, представленную на

рис. 1, б. Она включает операции шести уровней.

I уровень – получение исходной информации в результате геофизических

измерений и сбора априорных геологических данных.

II уровень – обработка и анализ указанной информации с различными

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

. годографов;

. горизонтальных спектров скоростей или графиков VОГТ;

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

и искажений и максимум объективной информации о строении среды.

Данные промысловой геофизики обрабатываются главным образом для

получения эффективной по сейсмическим критериям одномерной сейсмической

модели. Наконец, важнейшую роль, определяющую впоследствии все решения

геофизика-интерпретатора, играет предварительно выработанная гипотеза о

строении разреза, не противоречащая имеющимся геологическим представлениям.

III уровень состоит в создании исходной для итеративного процесса

интерпретации двумерной сейсмогеологической модели или модели нулевого

приближения. Эта операция в принципе неформальна и требует максимального

использования всей доступной информации I и II уровней. На этом же уровне

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

(моделирование сейсмического сигнала).

На IV уровне для получения модельных аналогов промежуточных и

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

прямые задачи сейсморазведки.

V уровень – операции сравнения промежуточных и окончательных

результатов обработки с их модельными аналогами, имеющие целью

количественную оценку сходства между ними.

VI уровень в рассматриваемой схеме представляют процессы принятий по

коррекции параметров в общем случае всех операций уровней II–V. В

частности, при наименее "глубокой" обратной связи корректируются параметры

сейсмомоделирования, т. е. сейсмогеологическая модель и модель импульса

падающей волны. Исходными данными для принятия таких решении являются

оценки сходства ("рассогласования"), получаемые на уровне V.

2 Теоретические вопросы автоматизированной интерпретации данных

сейсморазведки

Лекция 2

ТАБЛИЦА 1. ВЛИяНИЕ ПАРАМЕТРОВ ДВУМЕРНОГО СЕЙСМОМОДЕЛИРОВАНИя

НА ХАРАКТЕРИСТИКИ ОТРАЖЕНИЙ

|Кинематические и |Параметры |

|динамические | |

|характеристики | |

|отражений | |

|А. Определяемые по отдельным трассам синтетического временного разреза |

|1. Время отражения|Локальные мощности пластов вышележащей толщи |

| |Локальные скорости в пластах вышележащей толщи |

| |Геометрия отражающей и промежуточных границ |

|2. Амплитуда |Дифференциация скоростей и плотностей соседних слоев |

|отражения |Мощности слоев |

| |Количество слоев, участвующих в формировании отраженной|

| |волны |

| |Геометрия отражающей и промежуточных границ |

| |Частота исходного сигнала |

|3. Преобладающая |Частота исходного сигнала |

|частота отражения |Мощности слоев |

| |Количество слоев, участвующих в формировании отраженной|

| |волны |

| |Величины частотно-зависимого коэффициента поглощения |

|4. Полярность |Полярность исходного сигнала |

|отражения |Порядок чередования слоев |

| |Тип насыщающего флюида |

|5. Форма |Количество слоев, участвующих в формировании отраженной|

|отражения: |волны |

|а) длительность |Мощности слоев |

|волны, выраженная |Ширина спектра исходного сигнала |

|количеством фаз |Частота исходного сигнала |

|б) соотношение |Форма огибающей исходного сигнала |

|амплитуд |Количество слоев, участвующих в формировании отраженной|

|экстремумов (форма|волны |

|огибающей) |Дифференциация скоростей и плотностей соседних слоев |

| |Мощности слоев |

|Б. Определяемые по синтетическому временному разрезу |

|6. Поведение линий|Геометрия отражающей и промежуточных границ |

|t0 |Скорости и величины их градиентов в пластах вышележащей|

| |толщи |

| |Мощности пластов вышележащей толщи |

|7. Интерференция |Градиент изменения мощностей слоев, участвующих в |

|а) изменение |формировании отраженной волны |

|времени между |Градиент изменения скоростей слоев, участвующих в |

|соседними фазами |формировании отраженной волны |

|отражения | |

|б) изменения |Градиент изменения плотностей слоев, участвующих в |

|амплитуды |формировании отраженной волны |

|отдельных фаз |Криволинейность границ, участвующих в формировании |

|отражения |отраженной волны |

|(изменение формы | |

|огибающей) | |

|8. Когерентность |Градиент изменения мощностей слоев, участвующих в |

| |формировании отраженной волны |

| |Градиент изменения скоростей слоев, участвующих в |

| |формировании отраженной волны |

| |Градиент изменения плотностей слоев, участвующих в |

| |формировании отраженной волны |

| |Криволинейность границ, участвующих в формировании |

| |отраженной волны |

|9. Расположение и |Наличие и местоположение объектов дифракции (точки |

|интенсивность |выклинивания, примыкания; тектонические нарушения; |

|дифрагированных |резкие перегибы слоев, радиус кривизны которых меньше |

|волн |длины волны; участки резкого изменения пластовых |

| |параметров и т. п.) |

| |Дифференциация скоростей и плотностей в дифрагирующих |

| |телах и вмещающих породах |

Способы построения сейсмических моделей геологических сред

Предметом нашего рассмотрения являются волновые поля, образующиеся в

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

преимущественно продольные волны, наблюдения отраженных волн при достаточно

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

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

задавать в слоях модели следующие параметры: скорость продольных волн Vp,

плотность ( и коэффициент поглощения продольных волн (p. Поле продольных

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

параметрами, а распределение параметров поперечных волн не будет играть

существенной роли. Вследствие допущения о малых углах падения волны на

границы раздела анизотропия скоростей также не учитывается.

В большинстве случаев для построения двумерных моделей используется

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

геофизическая информация по разведочным скважинам и менее точная, но

существенно более плотная сейсмическая информация между скважинами. Первая

позволяет получить достоверные оценки физических свойств разреза в

отдельных точках, т. е. построить одномерные модели. С помощью второй

информации осуществляется переход к двумерным моделям.

1 Построение одномерных моделей

Исходная информация, т. е. значения детальных скоростей и плотностей,

для построения одномерных тонкослоистых моделей может быть получена

несколькими способами:

1. По данным акустического (АК), гамма-гамма (ГГК) или гравитационного

каротажей после соответствующей их обработки; обработка АК обычно

включает процедуры вычисления скоростей с учетом кавернометрии, коррекции

полученных скоростей по сейсмическому каротажу (СК), осреднения и др.;

ГГК дает сразу плотность, поэтому обработка его заключается только в

осреднении.

2. При отсутствии АК или ГГК, а также при низком их качестве акустические

свойства разреза прогнозируются с использованием других широко

распространенных промыслово-геофизических характеристик: кажущегося

сопротивления ((k), интенсивности первичного (ГК) и вторичного (НГК)

гамма-излучения и др.

3. Для приближенного задания акустических параметров тонких слоев иногда

используются нормальные или обобщенные зависимости скорости и плотности

от глубины для пород различной литологии.

Кроме того, информация о детальном распределении скоростей и

плотностей в разрезе может быть получена по данным изучения керна, однако

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

проводились в условиях, близких к пластовым.

Из перечисленных способов предпочтение следует отдать использованию

данных АК и ГГК.

Осреднение данных АК и ГГК

Большое количество данных АК, накопленное к настоящему времени,

подтверждает представления о тонкослоистой структуре реального скоростного

разреза. Практически все осадочные породы, за редким исключением (чистая

соль, лед), имеют тонкослоистую структуру с той или иной степенью

скоростной дифференциации.

Исходные непрерывные скоростные и плотностные разрезы,

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

одномерных моделей, по которым в дальнейшем предстоит построить двумерную

модель. Тем или иным способом производится их осреднение и построение

максимально упрощенной однородно-слоистой (или тонкослоистой) модели среды.

Такая модель представляется в виде серии тонких однородных пластов,

разделенных границами первого рода. При построении тонкослоистых моделей

предполагается, что акустическая неоднородность, обусловленная внутренней

изменчивостью пород пласта, незначительна по сравнению с межпластовой

акустической неоднородностью, связанной с изменением литологии или типа

насыщения.

Способ осреднения с порогом. Применение его позволяет получить

тонкослоистую модель в виде серии однородных слоев большей мощности по

сравнению с исходным разрезом. Все границы в модели представляются

границами первого рода. Сущность алгоритма осреднения в данном способе

заключается в том, что по заданным (V – величине значимой скоростной

дифференциации и ((min – минимальной временной мощности слоев из разреза

исключаются тонкие слои, время пробега в которых (ti < ((min, объединяются

слои с номерами i и i-1, если разница скоростей в них удовлетворяет условию

|Vi – Vi–1| ( (V

Значение скорости в объединенном слое вычисляется как среднее из Vi и

Vi-1. Пороговое значение скачка скорости (V может быть различным для разных

частей разреза.

Изменяя (V, можно менять число слоев в модели N, так как оно тем

меньше, чем больше (V. Это может быть использовано для автоматического

поиска моделей с числом слоев, находящихся в заданных пределах Nmin – Nmax.

2 Построение двумерных моделей

Рассмотрим методику построения двумерных сейсмогеологических моделей,

представляющих собой комбинацию толстослоистых толщ (покрывающей и

подстилающей) и собственно моделируемого интервала в виде совокупности

тонких слоев. Чтобы условия интерференции волн на верхней и нижней границах

моделируемого интервала не отличались от реальных, необходимо этот интервал

расширить вверх и вниз на величину не менее ( (длина волны). Пример

комбинированной модели представлен на рис. 8, д.

Такие модели используются, как правило, при решении стратиграфических

задач, в которых объектами исследования могут быть зоны выклинивания и

фациального замещения, залежи углеводородов и др. При этом моделируемый

интервал должен совпадать с объектом исследований. Желательно, чтобы в

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

данным глубоких скважин заданы одномерные модели. Когда на профиле или

вблизи него нет глубоких скважин, то в принципе возможно построение

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

Выбор комбинированного типа моделей для описания способов построения

самых разнообразных в целевом отношении двумерных моделей оправдан тем,

что:

. во-первых, такая модель получила наибольшее распространение в практике

моделирования и,

. во-вторых, излагаемые ниже способы пригодны как для построения

толстослоистых моделей (используемых при решении прямых и обратных

кинематических задач), так и для построения тонкослоистых моделей по

всему разрезу (используемых при решении прямых и обратных динамических

задач).

Однако на практике последние строятся очень редко из-за крайней

трудоемкости построения таких моделей в двумерном варианте. Поэтому тонкими

слоями задается ограниченный интервал, т. е. и в этом случае приходится

Страницы: 1, 2, 3, 4


Copyright © 2012 г.
При использовании материалов - ссылка на сайт обязательна.