КОМПЛЕКСНЫЙ АНАЛИЗ КИНЕМАТИКИ ДЛЯ НЕГИПЕРБОЛИЧЕСКОЙ КИНЕМАТИЧЕСКОЙ КОРРЕКЦИИ
COMPLEX ANALYSIS OF KINEMATICS FOR NON-HYPERBOLIC KINEMATIC CORRECTION

ОБЛАСТЬ ТЕХНИКИ
Настоящее изобретение в целом относится к структурному анализу сред, а именно к обработке сейсмических данных для анализа подземных структур.
УРОВЕНЬ ТЕХНИКИ
Известные способы построения сейсмических изображений используют систему сейсмических источников и приемников, для сбора данных о подповерхностных (т.е. подземных) структурах. В ответ на сейсмическое воздействие (например, взрыв или механический удар) на некотором источнике каждый из приемников создает сейсмотрассу, т.е. запись сейсмического сигнала на приемнике как функцию времени, за счет отражений воздействующей волны от подповерхностных слоев, располагающихся под системой. Трассы, полученные от множества приемников в ответ на воздействие источников возбуждения, расположенных в разных точках, подвергаются совместной обработке для построения изображения слоев.
В рамках этого процесса, чтобы улучшить отношение сигнал-шум, множество трасс, полученных от пар источник-приемник, окружающих общую глубинную точку (ОГТ) выравнивают по времени, а затем суммируют (подобного рода группу трасс принято называть выборкой, а процессы выравнивания и суммирования трасс известны как коррекция и суммирование нормальных кинематических поправок (НКП)). Выравнивание предназначено для учета разностей времен прохождения сейсмических волн между различными парами источник-приемник. Зависимость времени прохождения обычно представляют в виде гиперболической функции расстояния между источником и приемником, а временной сдвиг, применяемый для выравнивания трасс, чтобы скомпенсировать эту гиперболическую зависимость, вычисляют по несложной формуле, в которую входят расстояние между источником и приемником, а также скорость распространения волн в данной среде.
Известно, однако, что гиперболическая зависимость неточно отражает фактические характеристики подземных структур, и ее использование приводит к потере информации. Поэтому был предложен ряд альтернативных подходов для корректного выравнивания между выборками трасс. Например, в патенте US 5,103,429 раскрыт способ анализа таких структур с применением гомеоморфного отображения (в гомеоморфном образе каждый элемент анализируемого объекта отображается один к одному в соответствующий элемент его образа так, чтобы объект изучения и его образ были топологически эквивалентны). Такой способ, как утверждается, позволяет строить разнообразные виды суммарных разрезов и соответствующих образов без потерь разрешения.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
В описанных ниже вариантах реализации настоящего изобретения принципы гомеоморфного отображения применяются при обработке группы выборок трасс сейсмических сигналов, именуемой в дальнейшем супервыборкой. Кинематические поправки находятся адаптивно в процессе комплексного кинематического анализа, исходя из оценки параметров волновых фронтов, распространяющихся в подземной среде. Трассы суммируют с учетом этих кинематических поправок, для получения образа среды, который является более точным и подробным, чем те образы, которые могут быть получены известными способами. Описанные ниже способы не только повышают отношение сигнал/шум подповерхностных изображений, но и позволяют учесть такие факторы, как изменения скоростей, как для малоглубинных, так и глубокозалегающих подповерхностных областей, а также кривизну отражающих подповерхностных слоев, а кроме того, позволяют избежать эффектов неравномерного растяжения, свойственных известным способам кинематической поправки.
Хотя раскрытые ниже варианты реализации относятся именно к обработке сейсмических данных и построению образов подповерхностных структур, принципы данного изобретения можно аналогичным образом применять при построении изображений неоднородных сред иного рода.
Таким образом, в соответствии с одним из вариантов реализации настоящего изобретения, предлагается способ компьютерной обработки данных, включающий:
прием совокупности трасс, соответствующих сигналам, возникшим в результате отражения сейсмических волн от подповерхностных структур и принятым в различных точках в течение некоторого времени;
вычисление меры корреляции между трассами как функции набора параметров волнового фронта, которые задают соответствующие кинематические поправки, используемые впоследствии при выравнивании трасс;
составление матрицы, имеющей как минимум три измерения, причем одно из измерений соответствует времени пробега сейсмических волн, а как минимум два других измерения соответствующим образом относятся как минимум к двум параметрам волнового фронта, причем элементы этой матрицы содержат вычисленную меру корреляции;
нахождение с помощью матрицы значений параметров волнового фронта, дающих максимальное значение меры корреляции для, соответственно, каждого из множества времен пробега; и
составление сейсмического изображения подповерхностных структур путем выравнивания и интегрирования трасс с использованием кинематических поправок, заданных найденными величинами параметров волнового фронта.
В одном из раскрытых вариантов реализации получение совокупности трасс включает в себя получение супервыборки трасс, содержащей множество выборок трасс от различных соответствующих источников и приемников, находящихся в двухмерной или трехмерной геометрии системы наблюдений, причем нахождение этих значений включает обработку меры корреляции, вычисленной из множества выборок для того, чтобы найти значения параметров волнового фронта, дающих максимальное значение меры корреляции. Как правило, получение супервыборки включает отбор нескольких выборок, каждая из которых включает соответствующую серию трасс, для которой источники и приемники задают особую геометрию системы наблюдений, причем создание сейсмического изображения подповерхностных структур включает вычисление соответствующих сейсмических изображений в двух или трех измерениях, согласующихся с выборками.
Как правило, кинематические поправки находят с помощью аппроксимации, выбранной из набора аппроксимаций, состоящего из гиперболической, негиперболической и параболической аппроксимации.
В некоторых вариантах реализации как минимум два указанных параметра волнового фронта включают угловой параметр и, по крайней мере, один радиус кривизны волнового фронта, исходящего от подповерхностного отражателя. В одном из вариантов реализации угловой параметр включает в себя угол выхода, причем, по крайней мере, один указанный радиус кривизны включает в себя соответствующие радиусы общего отраженного участка и общего эволютного участка волновых фронтов.
Как правило, нахождение значений предусматривает составление множества многомерных матриц, имеющих как минимум три измерения и содержащих вычисленные значения меры корреляции для соответствующих различных значений радиуса общего эволютного участка, а также сложение этого множества многомерных матриц для составления итоговой многомерной матрицы, как минимум два из указанных измерений которой включают в себя первое измерение, относящееся к углу выхода, и второе измерение, относящееся к радиусу общего отраженного участка. В одном из раскрытых вариантов реализации сложение множества многомерных матриц включает внесение в элементы итоговой многомерной матрицы вычисленных значений меры корреляции, являющихся максимальными среди соответствующих элементов множества многомерных матриц. В дополнение к этому или вместо этого, нахождение значений может включать разделение итоговой многомерной матрицы на слои с получением множества двухмерных матриц и нахождение максимальных значений вычисленной меры корреляции в этих двухмерных матрицах.
В одном из раскрытых вариантов реализации как минимум два из указанных параметров волнового фронта включают угловой параметр и, по крайней мере, один параметр радиуса кривизны фронта волны и параметр скорости.
В некоторых вариантах реализации нахождение значений включает нахождение одного или нескольких максимальных значений вычисленной меры корреляции для каждого из времен пробега, а выравнивание трасс включает в себя выравнивание трасс согласно значениям параметров волнового фронта, соотносящихся с соответствующими максимальными значениями. Нахождение такого или нескольких соответствующих максимальных значений может включать сортировку максимальных значений по значению вычисленной меры корреляции для каждого из максимальных значений. В дополнение к этому или вместо этого, нахождение такого или нескольких соответствующих максимальных значений может включать группировку максимальных значений для каждого из времен пробега, для которых значения параметров волнового фронта различаются не более чем на заранее заданную пороговую величину, как относящихся к одиночной волне, а максимальных значений, для которых значения параметров волнового фронта различаются более чем на заранее заданную пороговую величину, как относящихся к разным волнам, при этом выравнивание трасс включает выравнивание трасс согласно значениям параметров волнового фронта для, по крайней мере, одной из этих разных волн.
В некоторых вариантах реализации нахождение значений включает отображение полученных данных, демонстрирующее вычисленную меру корреляции как функцию времен пробега и, по крайней мере, одного из параметров волнового фронта, а также получение входных данных от пользователя, который идентифицирует одно или несколько максимальных значений вычисленной меры корреляции. Как правило, получение таких входных данных от пользователя включает получение выбранного пользователем коридора, содержащего ряд максимальных значений для разных времен пробега, а также соответствующие значения, по крайней мере, одного из параметров волнового фронта, причем максимальные значения в коридоре соотносятся с волнами, отраженными, соответственно, от последующих подповерхностных слоев.
В другом варианте реализации способ включает определение некоего пространственного горизонта внутри некоторой области сейсмического изображения, и нахождение с помощью матрицы значений параметров волнового фронта, дающих максимальную меру корреляции для, соответственно, каждой из множества точек, расположенных вдоль пространственного горизонта.
Также, согласно одному из вариантов реализации настоящего изобретения, предлагается устройство для обработки данных, содержащее, в том числе:
интерфейс, присоединенный для приема совокупности трасс, соответствующих сигналам, возникшим в результате отражения сейсмических волн от подповерхностных структур и принятых в различных точках в течение некоторого времени; а также
процессор для обработки сигналов, сконфигурированный для вычисления меры корреляции между трассами как функции набора параметров волнового фронта, которые задают соответствующие кинематические поправки, вносимые при выравнивании трасс, для построения матрицы, имеющей как минимум три измерения, причем одно из этих измерений соответствует временам пробега сейсмических волн, и как минимум два измерения соответствуют должным образом как минимум двум из параметров волнового фронта, причем элементы матрицы содержат вычисленную меру корреляции, для нахождения, с помощью матрицы, значений параметров волнового фронта, дающих максимальное значение меры корреляции для, соответственно, каждого из множества времен пробега, для составления сейсмического изображения подповерхностных структур путем выравнивания и интегрирования трасс с использованием кинематических поправок, заданных найденными величинами параметров волнового фронта.
Согласно одному из вариантов реализации настоящего изобретения дополнительно предлагается компьютерный программный продукт, включающий машиночитаемый носитель, на котором хранятся команды программы, которые при считывании их компьютером заставляют компьютер получать совокупности трасс, соответствующих сигналам, возникшим в результате отражения сейсмических волн от подповерхностных структур и принятых на множестве точек в различных точках в течение некоторого времени, вычислить меру корреляции между трассами как функцию набора параметров волнового фронта, которые определяют соответствующие кинематические поправки, вносимые при выравнивании трасс, для составления матрицы, имеющей как минимум три измерения, причем одно из этих измерений соответствует временам пробега сейсмических волн, а как минимум два измерения соответствуют должным образом как минимум двум параметрам волнового фронта, причем элементы матрицы содержат вычисленную меру корреляции, для нахождения с помощью матрицы значений параметров волнового фронта, дающих максимальное значение меры корреляции для, соответственно, каждого из множества времен пробега, а также для составления сейсмического изображения подповерхностных структур путем выравнивания и интегрирования трасс с использованием кинематических поправок, заданных найденными величинами параметров волнового фронта.
Настоящее изобретение станет более понятно из нижеследующего подробного описания вариантов его реализации с прилагаемыми чертежами.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Фиг.1 представляет собой схематическую иллюстрацию системы построения сейсмических изображений, согласно одному из вариантов реализации настоящего изобретения.
Фиг.2 - 4 представляют собой схематические разрезы подповерхностной области, показывающие параметры, используемые для представления сейсмических волн, проходящих через эту область, согласно одному из вариантов реализации настоящего изобретения.
Фиг.5 представляет собой блок-схему алгоритма способа обработки сейсмических данных, согласно одному из вариантов реализации настоящего изобретения.
Фиг.6 представляет собой схематическую блок-диаграмму структур данных, используемых в рамках способа с фиг.5, согласно одному из вариантов реализации настоящего изобретения.
Фиг.7А представляет собой изображение суммарного сейсмического разреза (суммарный разрез), полученный одним из известных способов нормальной кинематической поправки.
Фиг.7В представляет собой суммарный разрез, полученный способом, основанным на оценке параметров волнового фронта, согласно одному из вариантов реализации настоящего изобретения.
Фиг.8 представляет собой блок-схему алгоритма способа комплексного кинематического анализа сейсмических трасс, согласно одному из вариантов реализации настоящего изобретения.
Фиг.9 представляет собой блок-схему алгоритма способа комплексного кинематического анализа сейсмических трасс, согласно одному из вариантов реализации настоящего изобретения.
Фиг.10 представляет собой схематическое графическое изображение структур данных, используемых в рамках способа с фиг.9, согласно одному из вариантов реализации настоящего изобретения.
Фиг.11 представляет собой схематическое графическое отображение структур данных, используемых в рамках другого способа комплексного кинематического анализа сейсмических трасс, согласно одному из вариантов реализации настоящего изобретения.
РАСКРЫТИЕ ВАРИАНТОВ РЕАЛИЗАЦИИ
Фиг.1 представляет собой схематическую иллюстрацию системы 20 для построения сейсмических изображений, согласно одному из вариантов реализации настоящего изобретения. Множество источников 22 и приемников 24 сейсмических сигналов размещены в ряде точек 26 на области 28 земной поверхности, лежащей над исследуемой под поверхностной областью. Как правило, источники 22 содержат заряды взрывчатки, которые подрывают в заданные моменты времени, тогда как приемники 24 содержат сейсмодатчики, генерирующие сейсмотрассы, на которых амплитуда сейсмических волн, дошедших до соответствующих точек 26, обозначена как функция времени. Хотя на фиг.1 для простоты показано, что источники и приемники расположены вдоль одной прямой линии (образуя, таким образом, двухмерную геометрию системы наблюдения), система 20, как правило, содержит ряд источников и приемников, распределенных по некоторой области земной поверхности и, таким образом, может использовать как двухмерную, так и трехмерную геометрию системы наблюдения.
Например, как показано на фиг.1, сейсмическое воздействие (например, взрыв или механический удар) на источнике 22 приводит к возникновению волн с лучевыми траекториями 30, распространяющимися вглубь Земли под поверхностью 28. Эти волны отражаются от подповерхностной структуры 32, например от стыка слоев различных пород. Отраженные волны по лучевым траекториям 34 возвращаются к поверхности 28, где их регистрируют приемники 24. Волны, отраженные от такой структуры, достигают различных датчиков в различные моменты времени в силу того, что под землей они должны преодолеть различные дистанции, а также в силу вариативности формы самой структуры 32 и наличия промежуточных слоев между отражающей структурой и поверхностью 28.
Процессор для обработки сигналов 36 принимает трассы от приемников 24 через предназначенный для этого интерфейс 38. Этот интерфейс может быть сконфигурирован для управления созданием воздействий, а также для приема в реальном времени сигналов с поверхности 28, или же он может получать запись трасс от соответствующей аппаратуры регистрации данных (не показана). Процессор 36 обычно содержит универсальный компьютер с программным обеспечением, сконфигурированным на исполнение описанных ниже функций. Такое программное обеспечение можно загружать в компьютер как электронным способом, например, по сети, так и с материального носителя данных, например оптического, магнитного или электронного запоминающего устройства. Как вариант, по крайней мере, часть функций процессора 36 может быть выполнена подходящим цифровым сигнальным процессором ЦСП (DSP), или другими видами специализированных или программируемых аппаратных устройств обработки данных. Обычно процессор 36 соединен с устройством вывода данных сейсмических данных и изображений, например, таких как дисплей 40. Для взаимодействия с пользователем к процессору 36 можно подключить одно или несколько устройств ввода, например, таких как клавиатура 42 и/или указательное устройство ввода 44.
Для получения сейсмического изображения области, лежащей под поверхностью 28, процессор 36 применяет к супервыборке сейсмограмм от приемников 24 принципы гомеоморфного построения изображений. Супервыборка может содержать совокупность выборок общей точки взрыва (то есть трасс, полученных результате единичного воздействия), выборок общей глубинной площадки (трасс, сгенерированных одним приемником в результате различных воздействий), и/или же сейсмограммы ОГТ (СМР), как описано выше. Принципы и правила, применяемые процессором при создании изображений подповерхностных структур, раскрыты ниже со ссылкой на фиг.2-4.
0
0
-
+
Фиг.2 представляет собой схематический разрез подповерхностной области под поверхностью 28, показывающий прохождение через эту область типичных сейсмических волн, согласно одному из вариантов реализации настоящего изобретения. На данной фигуре показаны некоторые из параметров, используемых при вычислении кинематических поправок в отношении структуры 32 для всех выборок супервыборки с центром в точке Х. Поправку здесь вычисляют для произвольных положений источника S и приемника R, расположенных на поверхности 28 на расстояниях, соответственно, X и X от точки Х. Хотя данный вариант реализации (как и все варианты в настоящей заявке) изображается для простоты в двухмерном пространстве X-Z, принципы данных реализации можно без затруднений применить к трехмерному пространству, в котором обычно производят построение сейсмических изображений.
0
0
Центральный луч 50 из точки X определяют как луч, падающий по нормали на структуру 32 в точке, обозначенной на фигуре как О. Луч 50 исходит от поверхности 28 и возвращается к ней под углом выхода β относительно нормали к поверхности. Сейсмический луч 52 из точки S отражается от структуры 32 в точке P и возвращается как отраженный луч 54 в точку R. В данном примере луч 52 пересекается с лучом 50 в точке F. Совокупность лучей 50, 52 и 54 можно, таким образом, смоделировать как сфокусированную волну, которая сходится от волнового фронта 56 к точке F, а затем отражается от структуры 32 и достигает поверхности 28 вдоль фронта 58 расходящейся волны. Волновые фронты 56 и 58 пересекаются с лучами 52 и 54 в точках S* и R* соответственно.
0
Лучи 52 и 54 представляют собой часть "супервыборки", содержащей множество выборок лучей, которые пересекают луч 50, причем их источники и приемники расположены с обеих сторон от начала центральной точки Х. Супервыборка такого рода может быть использована, как подробно раскрывается ниже, в технологии, применяемой для поиска значений параметров волнового фронта в окрестностях центрального луча 50. Эта технология определяет действительный и мнимый фокусы волн, отраженных от подповерхностных структур, и поэтому в дальнейшем именуется "многофокусной". В качестве альтернативы, для оценки значений этих параметров волнового фронта можно применять другие способы.
0
Если считать скорость распространения сейсмических волн во всей исследуемой области одинаковой, то времена пробега сейсмической волны на отрезке луча S*F будут равны времени пробега на отрезке луча FX, то есть , и аналогичным образом времена пробега на отрезке S*PR равны временам пробега на отрезке
0
0
X·ОX:
Время пробега волны от источника 5 до приемника R можно выразить как:
0
В этом уравнении первый член  равен времени пробега в обе стороны по центральному (с нулевым смещением) лучу XFО. Поправки времени  и  можно выразить через приповерхностную скорость в подповерхностной среде и параметры, относящиеся к центральному лучу, как раскрывается ниже.
0
CRE
CEE
Фиг.3 и фиг.4 представляют собой схематические разрезы подповерхностной области с фиг.2, показывающие параметры волнового фронта, используемые при вычислении локальных кинематических поправок в данной области, согласно одному из вариантов реализации настоящего изобретения. Для вычислений используют гомеоморфную модель, как раскрыто в вышеупомянутом патенте US 5,103,429. В следующих вариантах реализации в число этих параметров входят угол выхода β, показанный на фиг.2 (который равным образом можно определить и как угол входа); радиус R кривизны волнового фронта с общим участком отражения (ОУО (CRE)), показанный на фиг.3; и радиус R кривизны волнового фронта с общим эволютным участком (ОЭУ (СЕЕ) или нормальная волна), показанный на фиг.4.
0
CRE
CEE
Как показано на фиг.2 - фиг.4, параметры β, R и R относятся непосредственно к физическим характеристикам структуры 32. Как вариант, волновые фронты, используемые в рамках гомеоморфной модели, можно выразить через иные угловые и радиальные параметры, как будет очевидно для специалистов в данной области. Значения параметров волнового фронта не известны априорно, они скорее служат базисом для параметрической модели подповерхностных структур. Значения параметров задают подгонкой фактических данных трасс к этой модели в подробно описанном ниже процессе комплексного кинематического анализа.
CRE
0
0
0
Как показано на фиг.3, волна, отражающаяся в точке Р структуры 32, создает волновой фронт 60, достигающий поверхности 28. R есть радиус дуги 62 окружности (или сферы), которая аппроксимирует волновой фронт 60. Этот радиус может быть сопоставлен с обычной эффективной скоростью в подповерхностной области и по формуле, где V - средняя скорость сейсмических волн у поверхности 28, a t - время пробега сейсмической волны вдоль отрезка ХР.
CRE
CEE
Параметр R на фиг.4 есть радиус дуги 74 окружности (или сферы), который аппроксимирует ОЭУ (СЕЕ) волновой фронт 72, исходящий от отрезка 70 структуры 32 в сторону поверхности 28. ОЭУ (СЕЕ) волновой фронт 72 образован нормальными лучами, исходящими от различных точек отрезка 70 между точками a и b. Результирующий волновой фронт 72 между соответствующими точками A и B определяется, таким образом, формой структуры 32. R моделирует эту форму относительно виртуального центра 0, из которого исходит волна.
Теперь вернемся к уравнению (1). Поправки времени можно выразить следующим образом через геометрические параметры, показанные на фиг.2:
s
r
s
r
s
r
0
0
CRE
CEE
0
В этих выражениях R и R есть радиусы волновых фронтов 56 и 58 соответственно, а V и V - соответствующие приповерхностные сейсмические скорости в окрестностях источника и приемника. Если принять приповерхностную скорость постоянной, то есть V=V=V, то уравнения (2) и (3) можно объединить и преобразовать, как показано ниже, используя β, R и R, для получения кинематической поправки Δt для лучей вокруг центральной точки X. Приведенные ниже формулы были получены и опубликованы Берковичем в работе «The Multifocusing Method for Homeomorphic Imaging and Stacking of Seismic Data» (докторская диссертация, Университет Тель-Авива, 1995):
Параметр фокусировки σ зависит от положения точки F, в которой луч 52 или 54 пересекается с центральным лучом 50 (см. фиг.2). Этот параметр может быть выражен через радиусы как:
0
CRE
CEE
+
-
и меняется сравнительно медленно с изменениями β, R и R. От выбора параметра σ зависит применимое правило сортировки трасс. Например, для случаев X=0 и X=0 уравнение (7) дает σ=1 и σ=-1 соответственно. Эти случаи относят к схемам с общим источником и общим приемником. Промежуточные значения σ можно вычислить, решив уравнения (5) и (6).
±
±
a
а
Переменная Y есть несимметричный масштабный коэффициент, определяемый геометрией наблюдений. Он пропорционален функции распространения (которая пропорциональна X) и углу отражения волны Δθ: Y≈XΔθ, где  - постоянный коэффициент пропорциональности.
0
0
CRE
CEE
+
-
Для некоторой заданной выборки трасс вокруг центральной точки X значения X и X известны из координат источника и приемника для каждой трассы, тогда как величины β, R и R зависят от неизвестного местоположения и формы подповерхностных структур, от которых отражаются сейсмические волны.
0
CRE
CEE
0
CRE
CEE
0
0
CRE
CEE
+
-
Корректная комбинация величин β, R и R будет соответствовать данным по трассам в выборке, в том смысле, что если сдвиги по времени между трассами скорректированы согласно этим величинам, трассы будут иметь высокую степень когерентности (то есть высокое значение взаимной корреляции), поскольку отражения от структуры 32 будут выровнены во времени и, таким образом, хорошо скоррелированы. Процессор 36 производит поиск на интервалах комбинаций β, R и R с Х, X и X по всей супервыборке с применением процесса, именуемого здесь комплексным кинематическим анализом, для нахождения комбинаций, дающих наилучшее приближение для всех трасс. Результирующая карта β, R и R характеризует местоположение и форму структуры 32, а также других структур в исследуемой подповерхностной области. Этот процесс подробно раскрыт ниже.
0
Хотя приведенный выше вариант многофокусной модели и раскрытые ниже способы относятся непосредственно к ОГТ (СМР) выборкам вокруг центральной точки Х, эти способы можно сходным образом применять, с необходимыми изменениями, и к выборкам иных типов. Например, уравнение (4) и связанные с ним определения параметров можно переформулировать для выборок трасс с общей точкой взрыва или с общей глубинной площадкой.
Подразумевается, что такие альтернативные подходы не выходят за рамки настоящего изобретения.
Раскрытый выше многофокусный способ создает взаимно однозначное соответствие между параметрами волнового фронта и свойствами подповерхностных отражателей. Как вариант, можно использовать иные аппроксимирующие выражения и формулы для вычисления таких параметров для построения оптимального суммарного разреза. Некоторые из этих альтернативных способов можно рассматривать как особый случай раскрытого в настоящей заявке способа.
Например, Чира-Олива (Chira-Oliva) и другие предложили следующую формулу для вычисления кинематических поправок в работе «2-D ZO CR5 Stack by Considering an Acquisition Line with Smooth Topography» в журнале Revista Brasileira de Geofisica 23:1 (2005), страницы 1-18:
0
m
0
0
N
0
NIP
0
CEE
N
CRE
NIP
0
1
±
Здесь t - время пробега отраженной волны, t - время смещения нуля, - угол выхода нормального луча, x - координата средней точки, х - координата центральной точки, K - местная кривизна поверхности Земли в некоторой точке на линии сбора данных, h - половинное смещение, R - радиус кривизны в точке х гипотетического волнового фронта N (нормального), а R - радиус кривизны в точке x гипотетического волнового фронта из точки нормального падения NIP. Если переменные в уравнении (4) заменить их эквивалентными значениям в многофокусной модели (R=R, R=R, X=h, V=ν), и разложить уравнение (4) в ряд Тейлора, можно увидеть, что уравнение (8) есть попросту гиперболическая аппроксимация уравнения (4).
Другая формула для кинематических поправок предложена Хартвеком (Hertweck) и другими в публикации «Data Stacking Beyond CMP», The Leading Edge 26:7 (2007), страницы 818-827:
m
0
При замене переменных в уравнении (9) как , m=x-x, х=h, и , это уравнение также становится гиперболической аппроксимацией уравнения (4).
Теперь обратимся к фиг.5 и фиг.6, на которых схематически показан способ обработки сейсмических данных, согласно одному из вариантов реализации настоящего изобретения. Фиг.5 представляет собой блок-схему алгоритма, на которой показаны ключевые стадии способа, а фиг.6 представляет собой схематическую блок-диаграмму структур данных, используемых в рамках способа с фиг.5. Этот способ применяют при обработке двухмерных (2D) сейсмических данных, то есть данных, собранных посредством источников и приемников, расположенных вдоль одной линии. На практике, для построения трехмерной карты подповерхностных слоев, таким образом обрабатывают множество аналогичных двухмерных срезов, полученных на различных линиях. Как вариант, описанный ниже способ можно адаптировать для полностью трехмерного сейсмического исследования и обработки данных за счет добавления других параметров волнового фронта, в том числе поперечных углов и радиусов в трех направлениях. В этом случае описанные ниже трехмерные матрицы с фиг.6 будут заменены многомерными матрицами, имеющими более трех измерений.
В рамках способа с фиг.5 и фиг.6, на этапе 80 накопления данных, собирают множество сейсмических трасс. Трассы можно собирать, например, с использованием такого программного обеспечения системы, как показано на фиг.1, или с помощью любых других подходящих средств сбора сейсмических данных. Трассы образуют супервыборку 82, как описано выше.
CEE
CEE
0
0
CRE
0
0
0
CRE
CEE
0
1
2
j
J
1
2
l
L
Для анализа данных на этапе 84 определения ОЭУ (СЕЕ) задают интервалы возможных значений радиуса R. Формально, , где n - число подлежащих проверке радиусов. На этапе 86 вычисления степени когерентности, для каждого значения R процессор 36 вычисляет степени когерентности у трасс вокруг каждой центральной точки X для различных значений набора параметров β, R и t. Иначе говоря, для каждого отрезка времени t процессор берет разные комбинации значений β и R в пределах их соответствующих интервалов: β=[β,β,…β,…β] и , где J и K - соответствующие числа подлежащих проверке значений параметров, совместно с разными значениями R. Каждая комбинация  используется для вычисления соответствующей кинематической поправки Δτ для каждой трассы супервыборки относительно точки X, как определено в уравнении (4) с использованием формул (5), (6) и (7). Затем вычисляют степень s когерентности для всех трасс относительно центральной трассы для каждого значения Δτ=[Δτ,Δτ,…Δτ,…Δτ], где L - число трасс в супервыборке.
0
0
CRE
0
0
CRE
0
0
0
CRE
Таким образом, для каждой тройки значений  и каждого выбранного t, процессор 36 определяет одно соответствующее значение степени когерентности. Эти значения когерентности можно внести в кубы 88 данных, причем каждый куб отвечает своему значению . Каждый куб 88 заполнен элементами матрицы согласно их степени когерентности S(β, R, t) для соответствующего значения . Как показано на фиг.6, оси этих кубов - β, R и t, а число кубов равно N, как указано выше. Иначе говоря, всякий куб представляет собой трехмерную матрицу, в которой одно из измерений соответствует времени t пробега, а другие измерения отвечают параметрам β и R волнового фронта. Как вариант, степени когерентности можно поместить в одну четырехмерную матрицу или в набор кубов, измерения которых отвечают другим группам вышеописанных параметров волнового фронта.
0
CRE
0
Внутри каждого куба точки (β, R, t) с высокими значениями степени когерентности можно соотнести с действительными точками структуры 32 и других подземных слоев, от которых отразились сейсмические волны. Другие точки с высокой когерентностью, однако, могут соответствовать артефактам, обусловленным, например, многократными отражениями от различных слоев. Для корректного отображения подповерхностных слоев нужно выбрать и сгруппировать правильные точки с высокой степенью когерентности, а артефакты отбросить.
0
CRE
0
CEE
0
CRE
0
0
Чтобы облегчить выбор правильных точек в кубах 88, на этапе 90 консолидации по когерентности процессор 36 находит значения пары (β, R), дающие максимальные степени когерентности для каждого значения t по всем значениям R. Таким образом, процессор анализирует значения S(β, R, t) на интервале  для поиска таких пар , при которых S максимальна для каждого t, и в итоге создает единый результирующий куб 92, с вписанными величинами .
Процессор 36 рассекает куб на двухмерные матрицы 96, 98 и 100 на этапе 94 рассечения. Эти матрицы определены следующим образом:
- Матрицу 96, , получают извлечением максимальных значений когерентности по оси .
- Матрицу 100, , получают извлечением максимальных значений когерентности по оси .
0
0i
- Матрицы 98, , получают извлечением срезов степеней когерентности для множества различных значений t, обозначаемых как t для i=0…M.
На этапе 102 кинематического анализа матрицы 96, 98 и 100 используют в процессе комплексного кинематического анализа для отыскания наборов максимумов когерентности, отвечающих действительным отражениям от подповерхностных структур. Этот этап может исполняться процессором 36 либо автоматически, либо интерактивно с участием пользователя системы 20. Автоматические и интерактивные способы, которые могут быть применены на этапе 102, раскрываются ниже со ссылкой на фиг.8 и фиг.9. Используя найденные на этапе 102 максимумы, процессор 36 выравнивает и интегрирует трассы в супервыборке для того, чтобы получить законченное двухмерное или трехмерное изображение структур в исследуемой подповерхностной области. Это изображение, как правило, отражает поведение слоев под каждой из центральных точек на поверхности 28. Конкретные технологии для получения таких изображений на основе результатов анализа волнового фронта более подробно раскрыты в уже упоминавшейся докторской диссертации Берковича (Berkovitch).
Теперь обратимся к фиг.7А и фиг.7В, представляющим собой суммарные сейсмические разрезы, то есть сечения трехмерной подповерхностной среды, полученные на основе некоторого набора сейсмотрасс, собранного с применением системы наподобие показанной на фиг.1. На фиг.7А показан суммарный разрез, полученный известным способом ввода кинематических поправок и суммирования. На фиг.7В показан суммарный разрез, полученный путем вычисления гомеоморфных параметров и комплексного кинематического анализа, согласно одному из вариантов реализации настоящего изобретения. На фиг.7В структуры видны значительно подробнее, чем на фиг.7А. Процессор 36 достигает такого превосходного результата за счет обработки всей супервыборки, поскольку, как объяснялось выше, он использует более точный негиперболический ввод кинематических поправок, применяя параметры, связанные с направлением падения волнового фронта, скоростью в подповерхностной среде и кривизной отражателей.
0i
вх
m
m
m
Фиг.8 представляет собой блок-схему автоматического способа комплексного кинематического анализа, который можно применять на этапе 102, согласно одному из вариантов реализации настоящего изобретения. Этот способ осуществляют с использованием двухмерных матриц 98,  (или более многомерных матриц, приведенных к такому виду). На этапе 110 поиска максимума, для каждого значения t процессор 36 находит значения пар  и β, дающие максимальные значения степени когерентности. Затем, на этапе 112 группировки, процессор группирует смежные максимумы. Иначе говоря, если локальные максимумы функции когерентности достигаются при близких значениях  и β, процессор считает, что эти максимумы принадлежат к одной волне, то есть относятся к отражениям от одной и той же подземной поверхности. Чтобы определить, следует ли составлять пару из смежных максимумов, процессор вычисляет Δβ и , которые представляют собой разности между соответствующими значениями β и  на смежных максимумах. Если любая из этих разностей превышает некую предустановленную пороговую величину  или , то процессор делает заключение, что соответствующие максимумы принадлежат к разным волнам. Пороговые величины  и  определяют опытным путем.
0
CRE
CEE
0
0i
0i
0
После идентификации различных максимумов процессор 36 сортирует максимальные величины когерентности по убыванию на этапе 114 сортировки. Наибольшие степени когерентности относятся к волне с наибольшей корреляцией. Значение тройки (β, R, R) при этом максимуме определяет соответствующую кинематическую поправку Δτ для каждой трассы в супервыборке, согласно приведенному выше уравнению (4). На этапе 116 выравнивания процессор выравнивает трассы вокруг соответствующей центральной точки X согласно значениям этих поправок. Сумма выровненных трасс дает амплитуду центральной трассы для соответствующего отрезка времени t. Аналогично, процессор может выравнивать и суммировать трассы, используя значения поправок, выявленные по другим максимумам, для всех максимумов выше некоторого порога и/или до некоего заданного числа максимумов в порядке убывания. Процессор повторяет вышеописанную процедуру для всех значений t и всех центральных точек Х для построения законченного суммарного разреза.
1
2
3
0
0
0i
CRE
m
Теперь обратимся к фиг.9 и фиг.10, представляющим собой схематические графические изображения интерактивного способа комплексного кинематического анализа, который можно применять на этапе 102 согласно другому варианту реализации настоящего изобретения. Фиг.9 представляет собой блок-схему алгоритма способа, а фиг.10 - схематическое графическое изображение структур данных, используемых при применении этого способа. Для инициализации данного способа, на этапе вычисления матриц, процессор 36 вычисляет матрицы S, S и S (обозначенные на фиг.10 как, соответственно, матрицы 96, 100 и 98) для некоторого выбранного набора базовых центральных точек. Задействованные в процедуре базовые центральные точки могут быть либо выбраны пользователем, либо автоматически определены процессором и, как правило, совпадают с координатами центральных точек, используемыми в традиционном процессе. Результирующие матрицы для некоторой данной центральной точки Х показаны в виде полутонового изображения на фиг.10, причем степень когерентности в каждой точке данной матрицы представлена соответствующей яркостью изображения. Горизонтальные пунктирные линии на фигуре отвечают некоторому выбранному значению t, которое представляет собой значение t, используемое для матрицы 98. Вертикальные пунктирные линии на матрицах 96 и 100 отвечают местоположениям R и β для одного из максимумов (светлая зона).
1
1
st
0
1
Для удобства понимания, на этапе 122 пересчета матриц, можно выразить матрицу 96 через скорость, например эффективной скоростью. Для этого матрицу S когерентности пересчитывают как S(V,t) по, например, приведенной выше формуле . Новая матрица S, таким образом, описывает изменение скорости во времени.
1
st
CRE
0
На этапе 124 выбора скорости пользователь выбирает максимумы когерентности в матрице S. Обычно процессор 36 выводит на дисплей 40 графические данные, соответствующие изображению матрицы 96, как показано на фиг.10. Затем пользователь может посредством клавиатуры 42 и/или указательного устройства ввода 44 ввести данные, идентифицирующие "коридоры" максимумов на дисплее, например, таких как изогнутый столбец из диагональных светлых зон в матрице 96, видимый на фиг.10. Эти коридоры соответствуют последовательностям оптимальных величин V (или, эквивалентно, R) для соответствующих значений t. Пользователь без труда опознает и выделит единообразный строй последовательности максимумов в подобном коридоре. Единообразие последовательности максимумов является хорошим индикатором того, что максимумы отвечают действительным отражениям от последовательности приблизительно параллельных подземных слоев, а не артефактам. Как отмечено выше, пользователь опознает подобные коридоры на некоторых базовых точках, после чего процессор может способом интерполяции выявить соответствующие коридоры на прочих центральных точках.
CRE
0
CRE
Чтобы обеспечить выбор оптимальных R для каждого значения t, на этапе 128 выбора радиуса, процессор 36 может задействовать соответствующую матрицу 98 . На этом этапе процессору позволено выбрать наилучшее значение внутри коридора, заданного на этапе 124. В примере на фиг.10 оптимальное значение R соответствовало бы центру светлой зоны вдоль вертикальной оси матрицы 98.
0
0
0
Аналогично пользователь может опознать наилучшие значения β для соответствующих значений t и  на этапе 130 выбора угла. Для этого процессор 36 может отобразить матрицу 100 , а пользователь может опознать коридоры на дисплее, подобно тому, как это делалось на этапе 124. Коридор в матрице 100 отчетливо виден на фиг.10. Аналогично этапу 128, можно использовать матрицу 98 для подтверждения выбора наилучшего β.
0
CRE
0i
0
0i
0
Процессор 36 вычисляет соответствующие величины кинематических поправок исходя из выбранных значений β и R для каждого t, после чего, используя эти кинематические поправки, на этапе 134 выравнивает трассы вокруг соответствующей центральной точки Х. Проведение этого этапа аналогично описанному выше этапу 116, и подобным образом этот этап повторяется для всех значений t и всех центральных точек Х для построения законченного суммарного разреза.
CRE
CEE
0
CRE
Фиг.11 представляет собой схематическое графическое отображение структур данных, используемых в рамках другого способа комплексного кинематического анализа сейсмических трасс, согласно одному из вариантов реализации настоящего изобретения. В противоположность раскрытым выше способам анализа, которые относятся к кинематическому анализу параметров (β, R, R etc.) волнового фронта по оси времени (t), на фиг.11 показан анализ сейсмографических данных вдоль пространственного горизонта 140. Такой горизонт можно выбирать (автоматически или с участием пользователя системы 20) так, чтобы он соответствовал некой отражающей подповерхностной структуре. Для значений координаты X вдоль горизонта 140 процессор 36 строит куб 88, а затем извлекает матрицы 142 и 144 значений степени когерентности для того, чтобы найти максимумы R (или скорости) и, соответственно, β. Эти максимумы затем используют для выделения волн, отраженных от данного горизонта.
Хотя в раскрытых здесь фигурах и способах используются конкретные приемы комплексного кинематического анализа для опознания, выделения, сортировки и использования максимумов функции когерентности, аналогичным образом для этих целей к сейсмическим данным можно применить и другие автоматические и интерактивные приемы кинематического анализа, причем оставаясь в рамках настоящего изобретения.
Например, для поиска максимумов можно применить другие функции когерентности, такие как функции корреляции или критерий Стьюдента. Таким образом, в заявке и формуле настоящего изобретения ссылки на "меру корреляции" следует считать относящимися как к степени когерентности, так и ко всем другим сходным функциям корреляции и когерентности. Кроме того, при вычислении вышеописанных параметров можно использовать иные аппроксимации формул (2)-(7). Помимо этого вышеописанные технологии можно применять и к минимальным супервыборкам, содержащим только одну ОГТ (СМР)-выборку, или только одну выборку общей точки взрыва, или только одну выборку общей глубинной площадки.
Более того, хотя варианты реализации, раскрытые выше, относятся непосредственно к обработке сейсмических данных и построению изображений подповерхностных структур, принципы настоящего изобретения можно сходным образом применять при построении изображений негомогенных сред иного рода, например, в томографии или иных способах неразрушающего наблюдения. Определенные выше параметры, используемые при обработке сейсмических данных, можно также применять в сейсмотомографии, в рамках способа Острандера (AVO, амплитуда как функция удаления), для статической коррекции (коррекция влияния приповерхностных областей), и так далее.
Таким образом, подразумевается, что описанные выше варианты реализации приведены только для примера, в то время как настоящее изобретение не ограничивается непосредственно раскрытым в данном тексте и фигурах. Напротив, в рамки настоящего изобретения входят как комбинации и подкомбинации различных описанных выше признаков, так и вариации и модификации этих комбинаций, которые являлись бы очевидными для специалистов в данной области, но которые не раскрыты в известном уровне техники.