СИСТЕМА И СПОСОБ ВЕРТИКАЛЬНОГО СЕЙСМИЧЕСКОГО ПРОФИЛИРОВАНИЯ С ПРЕДСТАВЛЕНИЕМ РАЗВЕДОЧНЫХ ДАННЫХ В ВИДЕ КОМБИНАЦИИ ПАРАМЕТРИЗОВАННЫХ КОМПРЕССИОННОГО, СДВИГОВОГО И ДИСПЕРСИВНОГО ВОЛНОВЫХ ПОЛЕЙ
SYSTEM AND METHOD OF VERTICAL SEISMIC PROFILING WITH PRESENTATION OF SURVEY DATA IN FORM OF COMBINATION OF PARAMETERISED COMPRESSION, SHEAR AND DISPERSIVE WAVE FIELDS

Уровень техники
Разведка посредством вертикального сейсмического профилирования (ВСП) эффективно используется для измерения свойств геологических формаций, окружающих буровую скважину. В одной из технологий разведки посредством ВСП используется набор сейсмических датчиков, устанавливаемых в приблизительно вертикальной буровой скважине. Сейсмический источник создает сейсмические волны, исходящие из разных точек возбуждения на поверхности. Реакция датчиков на каждое из возбуждений регистрируется и анализируется для извлечения требуемой информации о свойствах формации.
Скорость компрессионных волн представляет собой одно из свойств формаций, для измерения которого широко применяется такой способ. Следует отметить, что компрессионные волны также могут называться волнами сжатия, продольными волнами, волнами давления, первичными волнами или P-волнами. Хотя для обозначения данного свойства широко используется термин «скорость», измеряемая величина обычно представляет собой скалярную величину, т.е. соответствует модулю скорости. Такой модуль скорости или «скорость» могут быть адекватно выражены в единицах медленности, т.е. величины, обратной скорости (другими словами, произведение модуля скорости и медленности равно единице).
Либо сам сейсмический источник, либо результат взаимодействия энергии сейсмической волны с разломами, границами между формациями и границами между твердыми и текучими средами (например, на поверхности земли или в буровой скважине) также могут порождать и волны других типов. В число таких волн входят сдвиговые волны и дисперсивные волны. Сдвиговые волны часто называются поперечными волнами, вторичными волнами или S-волнами. Дисперсивные волны представляют собой волны, частота которых не пропорциональна их волновому числу, т.е. волны с разными длинами волн распространяются с разной скоростью. В число примеров дисперсивных волн, некоторые из которых могут совпадать, входят каналовые волны, граничные волны, волны Лэмба, волны Лява, Q-волны, рэлеевские волны, волны Шолте, поверхностные волны, волны Стоунли и трубные волны.
В общем случае компрессионные волны имеют более высокие скорости и большие амплитуды, что облегчает их обнаружение и измерение, в частности потому, что волны этого типа первыми доходят до набора датчиков. Однако модуль сдвига представляет собой ключевое свойство формации, которое может быть определено только по измерениям скорости сдвиговой волны. К сожалению, приходящие дисперсивные волны и задержанные компрессионные волны (например, задержанные в результате отражения и преломления) могут маскировать достигающие набора датчиков сдвиговые волны, что затрудняет такие измерения и снижает их достоверность.
Краткое описание чертежей
На прилагаемых чертежах и в подробном описании представлены конкретные варианты осуществления систем разведки посредством вертикального сейсмического профилирования (ВСП), позволяющие разделить разведочные данные на компрессионные, сдвиговые и дисперсивные волновые поля. На чертежах:
Фиг. 1 иллюстрирует пример конфигурации разведки посредством ВСП.
На фиг. 2 представлен пример системы разведки посредством ВСП.
На фиг. 3 представлен пример сейсмических трасс.
На фиг. 4 приведена диаграмм примера способа разведки посредством ВСП.
На фиг. 5 представлен пример многокомпонентных разведочных данных ВСП.
На фиг. 6-8 представлены примеры соответственно компрессионного, сдвигового и дисперсивного волновых полей.
На фиг. 9а-9c представлены извлеченные параметры дисперсивного волнового поля.
На фиг. 10 представлен пример данных ВСП, восстановленных по волновым полям.
На фиг. 11 представлен пример компонент остаточного шума.
Однако следует иметь в виду, что конкретные варианты осуществления изобретения, представленные в чертежах и подробном описании, не накладывают каких-либо ограничений. Напротив, они образуют основу, позволяющую специалисту в данной области предусмотреть альтернативные варианты, аналоги и модификации, не выходящие за пределы объема охраны одного или нескольких вариантов осуществления изобретения, определенных в прилагаемой формуле изобретения.
Осуществление изобретения
Системы и способы по изобретению станут более ясны из проиллюстрированного примера их применения. На фиг. 1 представлен пример конфигурации разведки посредством вертикального сейсмического профилирования (ВСП), в котором операторы располагают набор сейсмических датчиков 102 в вертикальной буровой скважине 104 на расстоянии друг от друга. Датчики могут быть прикреплены к стенке скважины или зацементированы в нее для обеспечения возможности многокомпонентных измерений. Датчики 102 осуществляют беспроводную связь или связь через кабель с модулем 106 сбора данных, который принимает, обрабатывает и сохраняет данные сейсмических сигналов, собранные датчиками. Оператор включает сейсмический источник 110 (например, вибратор на шасси грузового автомобиля) в нескольких местоположениях («точках возбуждения») на поверхности 108 земли для формирования сейсмических волн, распространяющихся сквозь землю 112. Такие волны отражаются от разрывов акустического импеданса и достигают датчиков 102. В число примеров таких разрывов входят разломы, границы между пластами формаций и границы текучей среды. В представлении строения геологической среды, получаемом на основе данных сейсмических сигналов, такие разрывы могут иметь вид ярких пятен.
Кроме того, фиг. 1 иллюстрирует пример строения геологической среды. На схеме толща земли имеет три сравнительно плоских пласта подземных формаций и один наклонный пласт различного состава и, следовательно, с разными скоростями звука. Поры формаций могут быть заполнены газом, водой или нефтью, что также влияет на скорость распространения звука через формации.
Фиг. 2 иллюстрирует пример системы регистрации данных разведки ВСП, в которой датчики 102 соединены с шиной 302 для передачи цифровых сигналов в контур 306 регистрации данных. Информация о местоположении датчиков и другие параметры, используемые для интерпретации зарегистрированных данных, могут быть определены другими датчиками 304 и переданы в контур 306 регистрации данных для сохранения. Как вариант, такая дополнительная информация может включать точные местоположения датчиков и источников возбуждения, характеристики возбуждаемых источником волн, параметры оцифровывания, обнаруженные в системе разломы и т.д.
Каждый сейсмический датчик 102 может включать многоосевые акселерометры и/или геофоны и, в некоторых условиях, гидрофоны, каждый из которых может быть выполнен с возможностью осуществления выборок с высоким разрешением (например, с размером от 16 до 32 битов) с программируемой частотой выборки (например, от 400 Гц до 1 кГц). Контур 306 регистрации сохраняет потоки данных от датчиков 102 на долговременном носителе информации, например в массиве оптических или магнитных дисков. Данные сохраняются в виде сейсмических трасс (возможно, в сжатом виде), причем каждая из трасс представляет собой сигнал, который обнаружил и произвел из него выборку данный датчик в ответ на данное возбуждение. Местоположение источника возбуждения и датчика для каждой из трасс также сохраняются с указанием их принадлежности к данной трассе.
Система 308 обработки данных общего назначения принимает собранные данные разведки ВСП от контура 306 регистрации данных. В некоторых случаях система 308 обработки данных общего назначения физически соединена с контуром регистрации данных и обеспечивает возможность конфигурирования контура регистрации и осуществления предварительной обработки на месте. Однако чаще система обработки данных общего назначения расположена в центральном вычислительном узле с достаточными вычислительными ресурсами для интенсивной обработки данных. Передача данных разведки в центральный узел может быть осуществлена на физических носителях или через компьютерную сеть. Система 308 обработки содержит пользовательский интерфейс, содержащий графический дисплей и клавиатуру или другие средства пользовательского ввода, что обеспечивает возможность просмотра и анализа пользователями изображений и другой информации, полученной из данных разведки ВСП.
На фиг. 3 представлен пример сейсмических сигналов, которые могут быть зарегистрированы системой по фиг. 2. Эти сигналы отражают некоторую меру энергии сейсмических волн (например, смещение, скорость, ускорение, давление) как функцию от времени. Однако следует отметить, что число трасс обычно составляет несколько десятков, если не сотен, в связи с чем отдельное, изолированное отображение набора трасс, связанных с каким-либо определенным возбуждением, аналогично фиг. 3 на практике, как правило, невозможно. Вместо этого сигналы обычно представляют в «каскадном» формате, представленном на фиг. 5, в котором каждый из сигналов несколько сдвинут относительно сигналов, соответствующих соседним датчикам, но все они представлены в виде кривых, которые могут перекрывать друг друга. Перекрывающиеся кривые образуют рисунок, отражающий тенденции, имеющиеся в данных; так, например, наклонные линии указывают на приход сейсмических волн к набору датчиков.
На фиг. 4 представлена диаграмма способа разведки посредством ВСП, который может быть осуществлен с использованием системы по фиг. 2. На шаге 402 система получает данные разведки ВСП, как указано выше. На шаге 404 система формирует параметризованную модель волновых полей, имеющую параметры для по меньшей мере компрессионных волн, сдвиговых волн и дисперсивных волн. Структура такой модели подробно описана ниже. На шаге 406 система обеспечивает соответствие модели данным («подгоняет модель к данным»), используя способ нелинейной оптимизации для определения параметров, обеспечивающих наилучшее соответствие модели данным. По меньшей мере в некоторых из вариантов осуществления изобретения система использует для обеспечения наилучшего соответствия алгоритм Левенберга-Маркварта, но могут быть использованы и другие известные способы оптимизации, в том числе метод Гаусса-Ньютона, градиентный спуск, имитация отжига или оптимизация методом роя частиц. Предполагается, что в число параметров, определяемых для каждого из волновых полей, входят медленность и угол падения на набор датчиков.
Во многих случаях значения медленности волновых полей могут обеспечивать достаточную информацию для восстановления записей о требуемых свойствах формаций (например, модуля сдвига как функции глубины). В других случаях волновые поля параметризованной модели используются для дальнейшей обработки, так как уровень шумов в ней существенно ниже, чем в собранных данных. В связи с этим диаграмма на фиг. 4 включает шаг 408, на котором система строит на основе волновых полей параметризованной модели изображение геологической среды. Основные принципы построения сейсмических изображений хорошо известны и изложены в различных руководствах, в том числе в книге Jon F. Claerbout, Imaging the Earth′s Interior, Blackwell Scientific Publications, Oxford, 1985. Полученные таким образом изображения, профили или другие представления полученных свойств формаций отображаются системой на конечном шаге способа ВСП.
P
S
disp
P
S
disp
Углы падения волн P-волнового поля, S-волнового поля и дисперсивного волнового поля представлены в модели переменными θ, θ и θ соответственно. Эти углы падения вызывают распределение сейсмической энергии по вертикальной и радиальной компонентам сигнала согласно поляризационным векторам d, d и d.
P
S
disp
Тогда, если волновые поля имеют в частотной области волновую форму w(ω), w (ω) и w(ω), то двухкомпонентное смещение на (эталонном) датчике 0 может быть записано в виде:
Измерения на соседних датчиках могут быть выражены операторами временного сдвига в частотной области соответственно для P-волнового поля, S-волнового поля и дисперсивного волнового поля:
Z
P
S
phase
group
0
где Δ - расстояние между соседними датчиками, q и q - величины медленности (обратной скорости) P-волнового поля и S-волнового поля, q и q - величины фазовой и групповой медленности дисперсивного волнового поля, а ω - центральная частота волны дисперсивного волнового поля. Для набора из четырех датчиков уравнения модели имеют вид:
Разумеется, могут быть добавлены дополнительные датчики (и уравнения для датчиков). В обобщенном виде эти уравнения можно выразить следующим образом:
u(ω)=G(ω)w(ω),
P
P
S
S
disp
disp
0
phase
group.
где u(ω) - измеренные данные датчиков, w(ω) - вектор волновых полей, а G(ω) - параметризованная модель. Требуется определить восемь параметров: θ, q, θ, q, θ, q, ω, q и q Исходя из оцененного набора параметров можно найти соответствующее волновое поле с использованием метода наименьших квадратов:
причем расхождение между наблюдаемыми и модельными данными равно:
Задача алгоритма нелинейной оптимизации состоит в нахождении значений параметров, минимизирующих это расхождение. Может быть применен подход скользящего окна, в рамках которого одновременно анализируются сигналы от, например, 9 соседних датчиков. Помимо того что данный подход обеспечивает уменьшение объема вычислений, он также обеспечивает возможность изменения значений параметров при изменении местоположения для учета потенциальных изменений волновых полей в зависимости от глубины.
Кроме того, следует отметить, что вектор волнового поля может быть расширен для учета нескольких волновых полей каждого типа. Например, в уравнениях могут быть учтены восходящее P-волновое поле, нисходящее P-волновое поле, восходящее S-волновое поле, нисходящее S-волновое поле и нисходящее дисперсивное волновое поле. В соответствии с опытом и предположениями пользователя может быть выбрано большее или меньшее число волновых полей.
Фиг. 5-11 иллюстрируют пример применения систем и способов по изобретению. На фиг. 5 представлены измеренные при помощи геофонов сигналы разведки ВСП в терминах вертикального и радиального смещения. Наклонные линии, соответствующие нисходящим и восходящим волновым фронтам, ясно видны в обеих компонентах, хотя перекрывание волновых фронтов и образование ими интерференционных картин затрудняет их интерпретацию.
На фиг. 6-8 представлены волновые поля, извлеченные из данных по фиг. 5 с использованием вышеописанного способа. На фиг. 6 представлены нисходящее и восходящее P-волновые поля. На фиг. 7 представлены нисходящее и восходящее S-волновые поля. На фиг. 8 представлено дисперсивное волновое поле. Можно отметить, что интерференция между волновыми полями значительно уменьшена.
На фиг. 9a-9c представлены извлеченные значения параметров дисперсивного волнового поля по фиг. 8. На фиг. 9а представлена фазовая скорость как функция глубины. На фиг. 9b представлена групповая скорость как функция глубины. На фиг. 9c представлена центральная частота волны как функция глубины. Можно наблюдать постепенное уменьшение групповой скорости и центральной частоты по мере увеличения глубины. Фазовая скорость демонстрирует значительные колебания, но, по-видимому, не имеет систематической зависимости от глубины.
На фиг. 10 представлены восстановленные вертикальная и радиальная компоненты сигнала, полученные путем суммирования волновых полей по фиг. 6-8. Как и ожидалось, имеется значительное сходство с исходными данными с фиг. 5. На фиг. 11 представлены вертикальная и радикальная компоненты шумов, полученные путем вычитания восстановленных сигналов по фиг. 10 из исходных данных по фиг. 5. Видны слабые остаточные следы наиболее мощных компонент волновых полей (нисходящей P-волны и дисперсивной волны), которые можно объяснить наличием не учтенных в модели нелинейностей.
Для специалиста в данной области, ознакомившегося с вышеизложенным описанием, должны быть очевидны различные изменения и модификации, которые могут быть внесены в решение по изобретению. Подразумевается, что нижеприведенная формула изобретения охватывает все такие изменения и модификации.