Введение к работе
Актуальность проблемы. Разработка способов изучения околоскважинного пространства, включая- исследования анизотропии геологического разреза, является одним из важнейших направлений развития современной сквакинной сейсморазведки, применяемой при поисках месторождений нефти и газа. Места проведения таких геофизических исследований располагаются в различных нефтегазоносных регионах от шельфа морей и океанов до платформенных и рифтовых зон на континенте
Разработка методики изучения околоскважинного пространства проводилась в рамках совершенствования метода скважинной сейсморазведки в- процессе- бурения* В основе предлагаемой методики лежит возможность использования энергии породоразрушавдего инструмента в качестве источника непрерывных колебаний. Обработка регистрируемых на поверхности сейсмических сигналов ведется по аналогии с методом вибросейсморазведки на основе процедуры корреляции следящего сигнала посылки с полевой записью. Данная методика позволяет изучить пройденный бурением разрез и оперативно прогнозировать физические свойства геологической среды ниже забоя скважины, допуская практически любые пространственные схемы расположения сейсмоприемников на земной поверхности.
Повышению информативности и детальности геофизических исследований околоскважинного пространства может способствовать изученііе азимутальных изменений кинематических и динамических характеристик волнового поля, регистрируемого методом околоскважинного сейсмического просвечивания как по стандартной методике, так и по методике сейсморазведки в процессе бурения. Данные об анизотропных свойствах пород геологического разреза,
рассматриваемые в рамках гипотезы об упорядоченной трещиноватости среды, представляют большую ценность при изучении вещественного, состава, текстуры и нефтегазоносности горных поррд-. Эти данные позволяют прогнозировать возможные пути миграции углеводородов, а также учитывать закономерности распространения проницаемых зон, что особенно актуально при разработке нефтегазовых месторождений.
Целью работы является совершенствование сейсмических методов изучения околоскважинного пространства в том числе ноеого метода сейсморазведки в процессе оурения скважины (СПБ), как для решения задач структурной сейсморазведки, так и для оценки эффектов сейсмической анизотропии.
В соответствии с поставленной целью автотэом решается ряд конкретных задач, основными из которых являются:
1.Разработка методики продольной и разноазимутальной регистрации по методу СПБ на основе физических предпосылок и. специфического характера волнового поля , возбуждаемого долотом в скважине.
2. Разработка алгоритмического и программного аппарата для
обработки и интерпретации кинематических и динамических
характеристик волнового поля СПБ.
3. Проведение численного моделирования процесса
распространения реверберационнр- импульсного сигнала в геологи
ческой среде для проверки эффективности алгоритмов и программ.
4. Разработка интерпретационной базы для изучения
анизотропных свойств горных пород.
-
Опробование методики СПБ в реальных условиях.
-
Обработка и интерпретация данных скважинных сейсмических наблюдений .полученных как по традиционной (ВСП.СОГ) на шельфе Баренцева моря, так и по новой методике СПБ в пределах
Днепрово- Донецкой впадины.
Научная новизна.Разработана оригинальная методика обработки данных сейсморазведки в процессе бурения. Получена ноеэя аппроксимационная формула для анизотропии скорости продольных волн, позволяющая находить объемные плотности насыщенных и ненасыщенных флюидом трещин. Выявлен новый подход к интерпретации данных анизотропии кинематических и динамических характеристик упругих волн. Получены новые данные об анизотропных свойствах готаых пород в Днепрово-Донецкой впадине и на шельфе Баренцева моря.
Практическая ценность. Создана методика сейсморазведки, которая позволяет непосредственно в процессе бурения скЕажины на основе изучения характеристик проходящих или отраженных волн определять анизотропные свойства разреза, с последующим выявлением по этим данным направлений повышенной трещиноватости осадочного чехла. Полученная информация может бить использована при установлению! зон и направлений благоприятных для вертикальной и горизонтальной миграции углеводородов в новых нефтегазоносных районах и при выборе рациональной методики разработки, месторождений нефти и газа.
Реализация работы в производстве. По разработанной методике СПБ в 1991 г. были проведены полевые работы на скважине Скворцовская №1 на северном борту Днепрово-Донецкой впадины. Данные скважинных наблюдений были обработаны по оригинальным программам и были проинтерпретированы с учетом учетом анизотропии пород геологического разреза, связанной с разломной тектоникой данного района.
Методика интерпретации данных изучений анизотропии' пород околосквамшного пространства опробована на материалах треста "Севморнефтегазгеофизразведка" в Баренцевом море. Результаты
исследований вощли в научные отчеты ГАНГ им. И.М.Губкина по договору с трестом "Севморнефтегазгеофизразведка" ПО "Союзморгео" (1989-1991).
Апробация работы. Основные научные и практические результаты диссертационной работы обсуждались на Всесоюзной конференции "Роль молодежи в решении конкретных научно-технических проблем нефтегазового комплекса страны" (Красный Курган,1989), на Всесоюзном совещании "Современные метода геологической интерпретации геофизических данных при решении задач поисков и разведки залежей нефти и газа"(Краснодар, 1989), на научно- технической конференции "Построение физико-геологической модели и системный подход при истолковании результатов геофизических исследований"' ( Пермь, 1990). Теоретические предпосылки и предлагаемая методика скважинных геофизических исследований докладывалась на рабочих совещаниях специалистов в Восточно-Украинской Геофизической Разведочной Экспедиции ( ВУГРЭ Полтава, .1990, 1991 ), Киевской геофизической .экспедиции (Киев, 1990), а. также руководству и инженерной службе ИГО "Укргазгеофизика" (Киев* 1990 ).
Публикации. По- теме диссертации опубликовано 4 печатных работы.
Объем работы. Диссертация, состоит из введения, трех глав и заключения; содержит 110 страниц машинописного текста, 40 рисунков. Список литературы содержит 127 наименований.
Работа выполнена в Государственной академии нефти и газа им. И.М. Губкина на кафедре полевой геофизики в 1989-1991 годах под научным руководством профессора, д.г.-м.н. А.К.Урупова и доцента;:к.т.н. , A.M. Жукова , которым автор выражает глубокую благодарность.
Автор благодарит руководство и сотрудников ВУГРЭ и ПО
"Союзморгео" за проявленный интерес к результатам работы и за помощь при проведении опытных полевых ра.бот по новой методике. СОДЕРЖАНИЕ РАБОТЫ
В первой главе выполнен обзор сейсмических способов изучения околоскважинного пространства'. В- многочисленнных публикациях по скважинной сейсморазведке, которые, включают в себя такие традиционные методы как сейсмический каротаж (Пузырев Н.Н. 1957), вертикальное сейсмическое профилирование (Гальперин Е.И. 1971), межскважинное сейсмическое прозвучивание (Карус Е.В., Кузнецов О.Л., Файзулин И.О. 1986), метод обращенных годографов ( Теплицкий В.А. 1973 ), глубинное сейсмическое торпедирование (Силаев В.А. 1983 ) и другие, указывается , что с помощью этих методов можно решать достаточно широкий круг геологических задач при поисках и разведке месторождений нефти и газа. Отметим главные особенности проведения таких работ. Скважина , как правило, должна быть закончена' строительством и в качестве генераторов сейсмических колебаний на поверхности используются обычные взрывные или невзрывные источники упругих, волн. Проведение сейсморазведки в скважинах выполняется с применением либо специальных скважинных зондов, либо особых скважинных источников. Интерпретация данных обычно основывается на подходах взятых из традиционной сейсморазведки.
Вопросы изучения геологического разреза сейсмическими методами в процессе бурения сквакин , несмотря на казалось бы свою многолетнюю историю (РукаЕицын В.Н. 1980), тем не менее, разработаны далеко не полностью. Это относится как к методике проведения полевых наблюдений на скважине, так и к методике обработки данных. Анализ особенностей развития этого направления приводит к необходимости, разработки новых, более
эффективных способов исследования геологического строения околоскважинного пространства в процессе бурения, включая, изучение анизотропных свойств горных пород.
Вторая глава предваряется кратким обзором сейсмических способов исследования околоскважинного пространства в процессе бурения (СПБ) В основе методов СПБ лежит принцип использования вибрации, создаваемой долотом в процессе бурения, которое является источником упругих волн. Констатируется, что волновое поле СПБ характеризуется наличием интенсивных специфических волн-помех, которые не позволяют проводить выделение сейсмического сигнала стандартными средствами обработки. Из рассмотрения физических предпосылок методики СПБ и анализа волн- помех можно сделать . следующие выводы о характере возникающих колебаний и возможностях выделения полезных сигналов :
1) Частотные диапазоны колебаний, возбуждаемых наземным
оборудованием и породоразрушающим инструментом, практически
совпадают и лежат в пределах от.5 до 80 Гц.
2) Несмотря на различие волн- помех по кажущимся скоростям,
выделение полезных волн от долота путем -волновой фильтрации в
случае СПБ не эффективно , так как наземное оборудование
буровой работает как непрерывный поверхностный источник.
Исходя из вышеназванных особенностей колебаний, в качестге ' основного :способа обработки предлагается применить процедуру корреляции данных СПБ с записью следящего сигнала посылки (свип-сигнала) в надежде на то, что эта процедура позволит выделить из приходящих, непрерывно сейсмических волн сигналы, вызванные ударами долота, и подавить на коррелограммах-присутствие случайного шума.
Для целей исследования сейсмической анизотропии требуются
значительные затраты на отработку всех точек возбуждения, находящихся в разных азимутальных направлениях, и на разных расстояниях от скважины. Поэтому в данной работе отдается предпочтение обращенным методам скважишшх исследований, поскольку замена положения источника приемником имеет то преимущество, что приемники проще, и. дешевле размещать и работать с ними на земной поверхности. Пункты приема должны' быть размещены на земной поверхности в зоне благоприятной для регистрации колебаний, вызванных долотом при бурении. При наличии акустически жестких границ, могут быть,-прослежены и отраженные волны.
Выполняя регистрацию на нескольких глубинах в скважине, можно получить систему временных задержек падавдих и отраженных волн, характеризующих строение геологического разреза.
Поскольку бурильная колонна .являясь механическим каналом передачи упругих колебаний, обладает достаточно широкой полосой пропускания и практически без искажения передает вибрацию на поверхность, следящий сигнал посылки .(свип-сигнал) целесообразно регистрировать в верхней части буровой колонны на вертлюге.
Для оценки анизотропии кинематических и динамических характеристик упругих волн в процессе бурения автором разработана методика полевых наблюдений и регистрации, а с учетом поставленное задачи, предложена практическая схема расположения сейсмических профилей В рамках договора о научно- техническом сотрудничестве между ГАНГ, ВУГРЭ и "Укргазгеофизикой" под руководством аиторз были проведены опытные полевые работы на скважине Пкворцовская J61, расположенной на северном борту Лдепрово-Донецкой впадашы в 50 км западнее Харькова.
Перед обработкой стояла задача выделения прямых волн от
долота и отраженных волн от геологических границ чиже забоя скважины при наблюдениях на радиальной продольной расстановке. Критерием такого обнаружения может служить согласованность значений времен прихода и кажущихся скоростей с априорными сведениями о геологическом разрезе в районе проведения работ. При условии успешного решения данной задачи, было запланировано провести серию ^регистрации на кольцевой расстановке для обнаружения эффектов анизотропии в. волновом поле, связанных, предположительно, с трещиноватостью горных пород. :
Реально была отработаны1 продольная расстановка длиной 1380 м и полукольцо с радиусом 300 м. -Сейсмические наблюдения, в процессе роторного бурения'на глубине 2173 м, выполнялись с использованием стандартной сеисморазведочной аппаратуры "Прогресс-3".
Для обработки данных СПБ была разработана специальная методика «состоящая из двух основных этапов. .Первый этап предварительной обработки включал следующие процедуры :
1. Корреляция виброграмм" со следящим сигналом посылки
(свип-сигналом)..
2. Подавление шумов корреляции с помощью специально
разработанных программ.
3. Накопление коррелограмм с целью улучшения соотношения
сигнал-помеха.
4. Ввод статпоправки для учета задержки на время однократного
пробега по трубе от долота до датчика свип-сигнала наверху
буровой колонны.
Предобработка полевых данных выполнялась на ПЭВМ IBM PC/AT. Применению программ предшествовало их опробование путем численного моделирования процесса распространения реверберационно- импульсного сигнала в среде.
На втором- заключительном этапе обработки к сейсмограмме продольного профиля была .применена процедура веерной фильтрации подавляющая волны-помехи с низкими кажущимися скоростями. Затем применялась многоканальная обратная фильтрация, выоранная на осноеє трсг-лрования различных программ: селективной предсказывающей, нулъфазовой и многоканальной обратной фильтрации Расчет оператора деконволюции производился по прямой волне от долота.
В результате обработки получены коррелограммы, на которых удалось выделить оси синфазности с высокими кажущимися скоростями ( более 8000 м/с ), которые были идентифицированы как прямая и отраженная волн"
По результатам обработки были рассчитаны' вертикальные спектры скоростей, на которых четко выделяются .максимумы, соответсвующие прямой и отраженной-волнам. На спектрах видны кратные волны , вызванные реверберацией, обусловленной двойным пробегом волны по буровой колонне. Также выделяются кратные волны, вызванные реверберацией в слое ЗМС.
По вертикальным спектрам определены -эффективные скорости прямых и отраженных волн и расстояние до отражающей границы. Сделана попытка оценить -интервальную скорость между точкой забоя и границей. Эффективная скорость прямой волны сотавила V =2920м/с ,а эффективная.скорость отраженной волны VQ =3400 м/с. Рассчитанная глубина отражающей границы была . равна потр= 2610 м, что соответствует отраженному горизонту В-12 временного разреза ОГТ , проходящего через скважину. Величина интервальной скорости составила V „_= 5400 м/с, что в общих чертах не противоречит геологическим данным.
В третьей главе рассмотрены как теоретические, так и экспериментальные работы по . изучению сейсморазведкой
анизотропных свойств горных пород. Несмотря на то, что подобного рода работы ведутся давно, однако разработка способов интерпретации данных по анизотрогаш отстает от современных требований, а факты практического их использования пока малочисленны.
Известные способы интерпретации основаны на результатах математического моделирования. В. основе моделирования азимутальных изменений скоростей лежит процедура выявления фундаментальных связей, которая заключается в нахождении для заданного физического тела (модели) упругих констант матрицы Гука, как функции скорости волны в изотропной матрице, и некоторых параметров модели, ^например. Объемной плотности включений и др.). Выполняя подстановку найденных компонент упругости в формулу азимутального изменения скорости (Bacus 1964), переходим к модели анизотропии скоростей.
Интерпретация полевых измерений заключается в том, что, имея данные о скоростях упругих волн под разными азимутами и, задавшись той или иной моделью анизотропии, можно найти соответствие между модельными данными и полевыми .т.е. подобрать параметры соответствующей модели анизотропии скоростей. В итоге можно пытаться решить следующие задачи:
-
Определить направления простирания носителей анизотропии ( например,трещин Это не. вызывает особых затруднений, при условии, что; само явление .реально существует.
-
Определить параметры модели анизотропного тела (например, раскрцтость и флюидонасыщенность трещин -), рассчитать величину коэффициента анизотропии.
Первая задача решается сравнительно легко. Как известно из результатов из теоретических и экспериментальных исследований, выделение направлений анизотропии для модели с одной системой
трещиноватости основывается на следующих критериях: совпадение направлений носителей анизотропии (трещин) с максимумами, наблюдаемыми на индикатрисах скорости" , амплитуды и энергии; наличие 2-х взаимно-перпендикулярных плоскостей симметрии,, одна из которых проходит через Названное направление. При 2-х системах трещин, в целом, эти критерии, сохраняются, однако претерпевают некоторые изменения.
Для решения второй задачи использовалэсь поперечно-изотропная модель Гарбина-Кнопова ( Garbln.Knopoff 1974 ) и усовершенствованная С.Кремпиным (Crampln.S., Мс. Gonlgle R., Bamford D. 1980 ). По Гарбину-Кнопову -, модель представляет изотропную матрицу породы с низкой концентрацией тонких, имеющих форму монеты трещин, с преобладающей их ориентировкой в одном направлении."Считается , что тонкие, трещины имеют диаметр значительно меньший, чем длина данной сейсмической волны, в то время как общий объем породы с трещинами может быть соизмерим с длиной еолны. По С.Кремшшу в модели, действуют два фактора упорядоченной неоднородности: сухие трещины' и трещины, насыщенные флюидом.
Формулы азимутального изменения нормальной скорости продольных волн в случае сухих трещин и насыщенных флюидом имеют следующий вид: для сухих трещин:
Vp= Ко <1-(71/21)є-(8/3) Є COS28 +(1-/21.)ЄС034в) (1)
и для трещин, насыщенных .флюидом*
Vp = Vpo (1-(8/21)єг(8/21) є C0S46) (2)
где VQ - скорости Р и S волн в матрице анизотропной породы без трещин, V*o= 3V*o , е= N a / V - объемная плотность N сферических трещин с радиусом а в объеме породы V. Угол 8-измеряется от направления простирания трещины,, которое является осью лучевой -симметрии.
Совместное влияние сухих и насыщенных трещин на скорость упругих волн определяется формулой среднего времени, когда действуют два различных механизма замедления скорости пробега (Ї) волны.
1 1-р р
у у" у*
где объемные доли материала с насыщенными и "сухими" трещинами будут составлять р и 1-р соответственно.
Формула Кремпина была нами модифицирована путем уменьшения числа параметров в модели, что позволило упростить алгортим без потери точности аппроксимации.
г 8Є2
V =V /(1+-(1+ sln229+28 cds29+28) +— з1п*2Є) (4)
г ро 21 21
где ев - .объемная плотность, сухих трещин , es - объемная
плотность насыщенных трещин.
Создана программа, позволяющая находить аппроксимирующие
кривые и определять параметры модели в терминах: объемной
плотности сухих и насыщенных флюидом трещиь.
Правомерность .применения рассмотренной модели и модифицированной формулы к , реальным данным подтверждается
результатами физического моделирования-.
По результатам ультразвукового прозвучивания дефектоскопом УК-10 ПМС стеклянного диска с параллельными трещинами, выполненного Э.Б.Трушниковым (1991) для б моделей с различным числом трещин (1,3,7,15,31,63) нами' по формуле (1) были найдены параметры объемной- плотности сухих трещин.. Полученные в результате эксперимента скорости продольных волн И теоретические аппроксимирующие- кривые имеют хорошее СХОДСТЕО, которое увеличивается при увеличении числа трещин Установлено, что эффекты изменения скорости из-за анизотропии превышают .предельные погрешности измерений в 4-6 раз. Автором показана зависимость параметра объемной плотности трещин от их количества, которая хорошо аппроксимируется линейной регрессией, что также соответствует исходным условиям модели.
Для того, чтобы показать практическую значимость предложенного подхода к интерпретации были обработаны данные, полученные методом околоскважинного просвечивания на скв. Мурманская № 26, расположенной на одноименной структуре шельфа Баренцевого моря і
Морские сейсморазведочные работы на- скважине Ж6 проводились по разработанной для этого метода схеме и заключались в отработке судном-взрывпунктом серии кольцевых профилей с центром на скважине и радиусами 1 и 2 км. Прием упругих колебаний осуществлялся четырехточечным скЕажинным прибором с уровня горизонта на глубине 2460 м. В ходе работ контроль за положением судна-взривпункта обеспечивался гидрофоном, установленным у устья скважины, по временам прихода прямой водной волны от источника. На волновой картине, полученной на сейсмограммах кольцооал. профилей. были зарегистрированы как прлмая волна , так, и серия однократных и
многократных отраженных воля. По временным полям проходящей Р-_ волны для двух кольцевых профилей радиусов 1 и 2 км были построены азимутальные индикатрисы лучевых и нормальных скоростей, а такхе азимутальные индикатрисы коэффициентов, анизотропии скорости в вертикальной плоскости. Индикатрисы, рассчитывались из предположения об эллиптичности анизотропии с использованием материалов ВСП ближнего пункта взрыва для определения значения скорости, по вертикали. Необходимо отметить, что расчет всех кинематических характеристик выполнялся с учетом приходов водной волны, позволяющим более точно определить значения удаления судва-взрывпункта от скважины, но сравнению с использованием аналогичных данных навигационной системы.
Спектр сигнала' прямой проходящей Р-волны находится в диапазоне от 4 до 62 Гц. Расчет динамических параметров проводился, по спрямленной прямой волне в окне анализа 500 мс. Были построены индикатрисы максимума спектра мощности прямой волны, что соответствует частоте 24 Гц. Индикатрисы значений спектров, нормированных на собственную энергию волнового процесса, характеризуются более высокой разрешенностью по .сравнению с ненормированными, что находит подтверждение при их сопоставлении с индикатрисами кинематических параметров (скорости, коэффициентов анизотропии). На индикатрисах кинематических и динамических характеристик можно отметить два направления максимальных значений параметров: ССЗ-ИВ И ССВ-ШЗ.
Обработка морских сейсмических данных проходила в системе "SKGHPECC-BGII", разработанной сотрудниками кафедры полевой геофй'^ррси ( М.Б.Рашюпорт, В.И.Рыжков и др. -1990). Для обработки данных скоростей на интерпретационном этапе автором была разработана программа для ПЭВМ, которая позволяла в
интерактивном режиме проводить некоторые специальные^ процедуры обработки. Так для исключения влияния неоднородностей применялось: а) медианное сглаживание в скользящем окне; б) осреднение по взаимопротивоположным азимутам, которое признано правомерным, поскольку все теоретические вида анизотропии имеют лучевую симметрию.
Производился пересчет значений лучевых скоростей в нормальные , так как. теоретический зависимости (1)-(4) были получены для нормальных скоростей.
Для определения параметров, модели выполнялась аппроксимация данных теоретическими кривыми, рассчитываемыми по' формуле (4) , учитнвающей совместное влияние сухих и насыщенными флюидом трещин. В результате установлено, что наилучшая аппроксимация достигается при 7%'для сухих и 2;2 .% для насыщенных флюидом трещин для круговых наблюдений радиусом 2 км и соответственно 5,3% и 3% для круга радиусом 1 км.
На азимутальных индикатрисах скоростей.(V(a) ) и амплитуд (А(а)) отмечается два ортогональных направления, согласующиеся, с аналогичные направлениями,' выделенными на индикатрисах, энергии: направление ССЗ-ПВ (азимуты 300,315) и направление ГОЗ-ССВ (азимуты 30.45).
Непосредственно по результатам теоретической аппроксимации полевых индикатрис моделью (4) можно, сделать предположение, что анизотропия, пермотриасовой толщи Мурманской структури обусловлена, в основном, системой трещин ССЗ-ПВ направления, а штянутость "индикатрис в направлении. ШМЗСВ связана с частичным насыщением .этих,трещин флюидом;
. Полученные результаты в района шельфа Баренцева моря, в общих чертах, согласуются с направлениями вертикальной трещиноватости пород Кольского 'п-ова. крт'орый непосредственно
примыкает к району исследований. Имеется также согласованность полученных результатов с даннными исследований анизотропии на проходящих и преломленных волнах в Баренцевом и Карском морях ( А.К.Урупов, А.В.Степанов, В.И.Богоявленский, С.В.Добрынин 1989,1590,1991)
При исдледованиях анизотропии околоскважиннпго пространства в процесс.е бурения на скважине Скворцовская JM в Днепрово-Донецкой впадине, были получены эллипсовидные азимутальные индикатрисы значений мощности амплитудного спектра для гармоник в полосе полезного сигнала (частоты 18,38 Гц) и близкие к окружности для гармоник в полосе помехи (10 Гц). Индикатриса полезного сигнала образует фигуру, вытянутую в направлении С-Ю, а перепад значений мощности по отношению. З-Ф составляет 2-4 раза, что свидетельствует об анизотропных, свойствах горных пород в районе северного борта ДПВ.
Отмечается -хорошее соответствие выделенного направления анизотропии зоне древнейщих .разломов меридионального' простирания. Это же направление прослеживается на карте геофизических аномалий, выделенных по материалам электро- и гравиразведки.
Следует отметить, что при прозвучиваниг образца аргиллита (С^), отобранного с глубины 2950 м на этой же1 скважине, были получены" также .эллипсовидные индикатрисы кинематических и динамических характеристик, которые свидетельствуют о выраженных анизотропных свойствах пород, но на другом уровне -микроуровне. Коэффициент азимутальной.анизотропии скорости при этом составил 3%, а отношение максимальной амплитуды к минимальной оказалось больше двух.