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

Изобретение относится к области геофизики и может быть использовано для определения упругих свойств горных пород по сейсмическим данным. Заявлен способ определения упругих свойств горных пород на основе пластовой адаптивной инверсии сейсмических данных, характеризующийся применением пластовых моделей среды, в которых минимальные временные мощности τmin пластов соответствуют реальной разрешающей способности сейсморазведки и геологии осадконакопления и вычисляются согласно формуле:

τmin(мс)= 1 4 1000 Δ f , где Δf - рабочая полоса частот. В качестве исходной модели используют откорректированную геоакустическую модель импедансов соответствующей скважины, либо трассу импедансов полученного ранее сейсмоакустического разреза по секущему профилю. При построении рабочей модели преобразуют сейсмический временной разрез в детальную пластовую модель акустических импедансов (сейсмоакустический разрезом), по которой свидетельствуют о диапазонах изменения акустических параметров слоев. Технический результат - повышение точности оценки упругих свойств горных пород. 1 ил.

 

Изобретение относится к способам расчета упругих свойств горных пород: продольного импеданса (Ip), поперечного импеданса (Is) и плотности (ρ) по сейсмическим данным.

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

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

Прогнозирование свойств коллекторов по данным сейсморазведки выполняется способами динамической интерпретации (ДИ), основным элементом которой является сейсмическая инверсия, осуществляющая оптимизационный подбор модели среды. На сегодняшний день практически во всех имеющихся отечественных и зарубежных системах используются способы инверсии, основанные на тонкослоистых моделях среды. В этих моделях временная мощность пластов постоянна и равна шагу дискретизации, что соответствует толщинам слоев порядка 3-4 м. Поскольку вследствие регуляризации получаемые такими способами кривые импедансов и плотностей являются сглаженными, уместно назвать эти способы термином «непрерывная инверсия». Теоретические исследования и практические примеры опробования показали, что непрерывная инверсия дает надежные оценки акустических свойств достаточно мощных (15-20 м) однородных коллекторов, тогда как тонкие пласты-коллекторы мощностью до 10 м либо не отображаются на результатах инверсии, либо их акустические свойства оцениваются с недопустимо большими погрешностями (до 30-40%).

Для обработки 2D сейсмических данных нами разработана система динамической интерпретации OTRI (в русской аббревиатуре ОТДИ), в состав которой входит оригинальная процедура акустической инверсии, основанная на применении пластовой модели среды («пластовая инверсия»). Эта технология позволяет оценивать акустические свойства пластов с высокой точностью (единицы процентов). Как показали наши исследования в случае различия акустических свойств коллекторов и вмещающих пород, применение пластовой инверсии позволяет оценивать эффективные мощности коллекторов, начиная с Нэф=2-3 м. Данная точность достаточна для построения детальных геологических моделей, используемых не только на этапе разведки, но и на этапе эксплуатации месторождений, например при дизайне гидроразрыва пласта.

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

Выявление коллекторов и изучение распределения их емкостных свойств по площади, что необходимо для разведки и рациональной эксплуатации месторождений, осуществляется в результате применения динамической интерпретации (ДИ) сейсморазведки в комплексе с геофизическим исследованием скважин (ГИС). Основным элементом динамической интерпретации сейсморазведки является сейсмическая инверсия, осуществляющая трансформацию сейсмических волновых полей в детальные модели среды, описывающие распределение в пространстве акустических параметров: скоростей Vp и Vs продольных и поперечных волн, а также плотностей ρ. Совокупная интерпретация полученных по данным сейсморазведки прогнозных акустических параметров и имеющихся данных ГИС позволяет прогнозировать литологию пород, емкостные свойства коллекторов и характер их насыщения (газ или жидкость).

В практике интерпретации данных сейсморазведки используются разнообразные способы одномерной (акустической) и двумерной (упругой) сейсмических инверсий. Одномерная сейсмическая инверсия обрабатывает трассы суммарных сейсмических разрезов с целью получения акустических импедансов Ip=ρ*Vp и используется для оценки емкостных свойств коллекторов. Двумерная инверсия обрабатывает исходные сейсмограммы (до суммирования) продольных волн, получает все три акустических параметра и используется для прогнозирования литологического состава пород и характера насыщения коллекторов.

Все программы сейсмической инверсии, содержащиеся в широко применяемых зарубежных интерпретационных системах, основаны на тонкослоистых моделях среды, в которых минимальные мощности пластов составляют 3-4 м. Казалось бы, что именно такие программы являются наиболее эффективными для прогнозирования свойств маломощных коллекторов, имеющих подобные толщины. Однако вследствие неминуемого ограничения частотного диапазона сейсморазведки, ее реальная разрешающая способность (обычно 10-15 м) не позволяет надежно оценивать параметры отдельных пластов, слагающих используемые тонкослоистые модели [Veeken P., Da Silva. M., 2004]. Теоретические исследования и многочисленные практические эксперименты показали, что эти программы инверсии позволяют с высокой точностью (единицы процентов) определять параметры однородной группы тонких слоев, объединяющихся в более мощные пласты (от 15-20 м и более) [Buland A., Omre H., 2003, Bosch et al., 2010]. То есть такие способы инверсии дают хорошие результаты при прогнозировании свойств достаточно мощных коллекторов, но не пригодны для исследования коллекторов, мощность которых не превышает 10-15 м. На полученных практических результатах видно, что ошибки оценки акустических свойств отдельных тонких пластов с резко отличными акустическими свойствами, мощность которых не превышает 5-10 м, могут достигать 30% и более.

Выход из казалось бы неразрешимой ситуации с исследованием маломощных коллекторов состоит в использовании при выполнении сейсмической инверсии пластовых моделей среды, в которых минимальная мощность пластов соответствует реальной разрешающей способности сейсморазведки в конкретных сейсмогеологических условиях. При этом имеется возможность оценить с высокой точностью акустические параметры пластов, мощность которых составляет 15-20 м. Это оказывается достаточным для того, чтобы выявить наличие в таких пластах маломощных коллекторов (от 2-3 м), отличающихся по акустическим свойствам от вмещающих пород, а также количественно оценить их эффективные толщины Нэф, а в некоторых ситуациях - их линейные емкости Е=Нэфпп - коэффициент пористости).

Нами разработана оригинальная отечественная система одномерной инверсии ОТДИ, основанная на использовании пластовой модели среды [Кондратьев И.К. и др., 2005, Рыжков В.И. и др., 2008]. Эта система при практическом применении с использованием 2D сейсморазведки на десятках месторождений, расположенных в России и за рубежом (Казахстан, Узбекистан, Китай, Алжир), показала высокую эффективность прогнозирования маломощных коллекторов, в том числе и в сложных сейсмогеологических условиях Восточной Сибири, а также в подсолевых отложениях Прикаспия [Кондратьев И.К. и др., 2010].

Однако на сегодня существует только одномерная пластовая инверсия, применяемая для обработки данных 2D сейсморазведки. Именно из-за этого она не находит широкого применения, и нефтяные компании вынуждены использовать импортные программные продукты. Необходимо развить подход, реализующий пластовую инверсию, на вариант двумерной инверсии и сделать пластовую инверсию пригодной для обработки данных 3D сейсморазведки. Необходимо также опробовать новые варианты пластовой инверсии в основных нефтегазоносных провинциях России и разработать методические рекомендации по ее использованию для эффективного прогнозирования маломощных коллекторов.

Известен «СПОСОБ ПОИСКА, РАЗВЕДКИ, ИССЛЕДОВАНИЯ, ОЦЕНКИ И ПРОГНОЗИРОВАНИЯ РАЗРАБОТКИ ЗАЛЕЖИ И МЕСТОРОЖДЕНИЯ ПОЛЕЗНЫХ ИСКОПАЕМЫХ» (RU 2001120889, опубл. 2004, Миколаевский Э.Ю.).

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

Предлагаемый нами метод предназначен для изучения пластов малой мощности и основывается на модели физики пласта (rock physics).

Известна «АВТОМАТИЗИРОВАННАЯ СИСТЕМА ПОИСКА, РАЗВЕДКИ, ИССЛЕДОВАНИЯ, ОЦЕНКИ И ПРОЕКТИРОВАНИЯ РАЗРАБОТКИ ЗАЛЕЖИ И МЕСТОРОЖДЕНИЙ ПОЛЕЗНЫХ ИСКОПАЕМЫХ» (RU 22831U, Миколаевский Э.Ю.).

Данная система предназначена для реализации способа поиска и разведки, описанного в заявке RU 2001120889. Автоматическая обработка большого количества разнородных параметров требует мощного логического блока, который является центральным звеном данной системы. Устойчивые статистические связи разнородных параметров характерны только для массивных однородных пластов. Сложноустроенные коллекторы малой мощности не находят отклик в большинстве геофизических полей, регистрируемых на поверхности, и имеют сложные корреляционные связи геофизического параметр-петрофизического свойства.

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

Известен патент CN 2577296 «UNDERGROUND MULTI-FREQUENCY ACOUSTIC COMBINED DETECTOR», HOU CHUNHUI [CN]; CHU ZEHAN [CN]; LI JIANHAO [CN].

В данном решении для измерения акустических свойств пластов предлагается многоканальный многочастотный зонд. Данный прибор позволяет изучать свойства пласта-коллектора с высоким разрешением (десятки сантиметров), но на расстоянии не более 1 м от ствола скважины.

Наш метод позволяет изучать акустические свойства пласта с меньшим разрешением (первые единицы метров), но на существенном удалении от скважины (десятки километров), что необходимо при неоднородном строении пласта.

Известен патент EP 0889331 «METHODE POUR MODELISER EN 3D L'IMPEDANCE D'UN MILIEU HETEROGENE», Grizon, Laurent, Leger Michel, Richard Vincent и др.

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

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

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

Известен «СПОСОБ ДЛЯ ОБНАРУЖЕНИЯ ПОТЕНЦИАЛЬНОГО ВОДНОГО ПОТОКА С МАЛОЙ ГЛУБИНОЙ ЗАЛЕГАНИЯ» (RU 2319983, за авт. Надир Дуттаб, Субхашис Маллик).

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

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

Известна международная заявка WO 2010/092084 «FINITE ELEMENT MODELLING OF BOREHOLE SEISMIC RESPONSES IN LAYERED ANISOTROPIC FORMATIONS AND ITS USE FOR ELASTIC INVERSION», JORGENSEN Ole.

Данное решение также является ближайшим аналогом и посвящено решению прямой задачи (моделирование волнового поля при известных свойствах среды) в слоистых анизотропных средах.

Наша работа посвящена способу решения обратной задачи (прогноз неизвестных свойств среды на основе зарегистрированного волнового поля).

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

Для достижения поставленной цели необходимо решить следующие задачи:

a) разработка алгоритмов и методики пластовой акустической инверсии данных 3D сейсморазведки;

b) разработка алгоритмов и методики пластовой синхронной инверсии данных 3D сейсморазведки;

c) разработка методических рекомендаций оптимального комплексирования акустической и синхронной пластовой инверсии.

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

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

τmin(мс)= 1 4 1000 Δ f , где Δf - рабочая полоса частот, в которой, после всех примененных процедур обработки, уровень сигналов превышает уровень; причем в качестве исходной модели используют откорректированную геоакустическую модель импедансов соответствующей скважины, либо трассу импедансов полученного ранее сейсмоакустического разреза по секущему профилю; а при построении рабочей модели преобразуют сейсмический временной разрез в детальную пластовую модель акустических импедансов (сейсмоакустический разрезом), по которой свидетельствуют о диапазонах изменения акустических параметров слоев.

Благодаря тому что способ основан на применении пластовых моделей среды, в которых минимальные временные мощности τmin пластов соответствуют реальной разрешающей способности сейсморазведки и геологии осадконакопления, эмпирическая формула, связывающая τmin с качеством исходных сейсмических материалов (отношение сигнал/помеха, реальная полоса частот сейсмических сигналов), применение пластовой сейсмической инверсии обеспечивают максимальную точность оценки упругих свойств горных пород. Другие способы сейсмической инверсии основаны на применении эквидистантных моделей среды с постоянной мощностью пластов, равной интервалу дискретности Δt сейсмической записи. Реализуемое в пластовой инверсии согласование детальности получаемых сейсмоакустических разрезов с качеством исходных материалов и границ геологических пластов позволяет эффективно использовать метод для прогнозирования коллекторов в сложных сейсмогеологических условиях.

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

Осуществление изобретения

Перечень исходных данных, требуемых для реализации способа:

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

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

t 0 k ( x ) - линия t0(х) для k-го реперного пласта, полученная в результате корреляции отражений от кровли соответствующего пласта на временном разрезе;

I p k min , I p k max - граничные значения продольных импедансов для k-го реперного пласта;

I s k min , I s k max - граничные значения поперечных импедансов для k-го реперного пласта;

ΔIpmax,ΔIsmax - максимальное значение изменения импедансов на одну итерацию.

Для задания τmin используют эмпирическую формулу, характеризующую оптимальное значение этой величины по критерию точности аппроксимации «истинной» детальной модели соответствующей трассой полученного сейсмоакустического разреза:

τmin(мс)= 1 4 1000 Δ f

В этой формуле Δf - рабочая полоса частот, в которой, после всех примененных процедур обработки, уровень сигналов превышает уровень.

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

Главным вопросом в задании линий t 0 k ( x ) является выбор реперных горизонтов, осуществляемых в соответствии со следующими условиями. Во-первых, реперные пласты должны быть распространены на всей площади и характеризоваться выдержанными упругими свойствами (обычно плащевидные глины или маркирующие карбонаты). Во-вторых, желательно, чтобы эти пласты располагались на границах (или вблизи от них) контрастирующих по акустическим свойствам мощных толщ (сотни метров). В случае выполнения первого условия, диапазоны изменения импеданса ΔIp= I p p max I p p min в реперных пластах составляют обычно 800 - 1500 м/с·г/см3(ед). Диапазоны изменения Ip в межреперных пространствах ( I p k , k + 1 min , I p k , k + 1 max ) определяются в результате анализа кривых Ip(H) по всем имеющимся скважинам. Величина ΔIpmax - максимально допускаемая корректировка Ip в пластах модели при одном шаге интерактивной оптимизации, задается обычно в диапазоне (0.35-0.5) от минимальной величины ΔIp (диапазонов изменения Ip в реперных пластах). Аналогично ограничения задаются для поперечных импедансов Is.

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

Постановка задачи пластовой акустической инверсии

Введем следующие обозначения:

s ¯ =[s1, s2, …, sN]Т - сейсмическая трасса, представленная вектором - столбцом отсчетов s(t) с шагом дискретизации Δt, в интервале tнач, tкон;

p ¯ =[Ip1, Ip2, …, Ipm, t1, t2, …, tm]Т= [ I p ¯ ,   t ¯ ] - произвольный вектор параметров модели, имеющий ту же структуру, что и q ¯ ;

m ¯ ( p ¯ ) =[m1, m2, …, mN]Т - модельная трасса для вектора p ¯ ;

y ¯ = s ¯ m ¯ ( q ¯ ) - разностная трасса между s ¯ и модельной трассой для начального приближения q ¯ ;

x ¯ = p ¯ q ¯ - корректирующий вектор, приводящий q ¯ к вектору p ¯ = q ¯ + x ¯ .

Целью акустической инверсии является преобразование сейсмического временного разреза s(t, х) в детальную пластовую модель акустических импедансов Iр(t, х), называемую сейсмоакустическим разрезом. Основой нашего алгоритма является обобщенная линейная инверсия, которая модифицирована нами для введения априорной информации о диапазонах изменения акустических параметров слоев. Ниже приводится описание математической основы модифицированного алгоритма обобщенной линейной инверсии (МОЛИ).

Базовые уравнения обобщенной линейной инверсии

Реальную трассу s ¯ представим в виде:

s ¯ = m ¯ ( p ¯ ) + n ¯ ( p ¯ ) , (1)

где n ¯ ( p ¯ ) - помеха, под которой понимается все то, что отличает s ¯ от m ¯ ( p ¯ ) при выбранной модели p ¯ . Если полагать, что p ¯ близка к q ¯ , можно считать, что m ¯ ( p ¯ ) m ¯ ( q ¯ ) линейно зависит от ( p ¯ q ¯ ) . Для того чтобы выразить это математически, используют разложение m ¯ ( p ¯ ) m ¯ ( q ¯ ) в ряд Тэйлора с отбрасыванием всех членов, кроме первого, что дает:

m ¯ ( p ¯ ) m ¯ ( q ¯ ) + i m ¯ ( q ¯ ) q i ( p i q i ) . (2)

Уравнение (2) можно представить в матричной форме:

m ¯ ( p ¯ ) = m ¯ ( q ¯ ) + A x ¯ , (3)

или:

[ m 1 ( p ¯ ) m 2 ( p ¯ ) . . m N ( p ¯ ) ] = [ m 1 ( q ¯ ) m 2 ( q ¯ ) . . m N ( q ¯ ) ] + [ m 1 q 1 m 1 q 2 . . m 1 q M m 2 q 1 m 2 q 2 . . m 2 q M . . . . . . . . m N q 1 m N q 2 . . m N q M ] [ p 1 q 1 p 2 q 2 . . p M q M ] (4)

Далее под p ¯ будем понимать оптимальное решение, которое находится в результате применения нашего алгоритма. Одно из условий оптимальности:

m ¯ ( p ¯ ) s ¯ 2 = min (5)

( u ¯ означает норму пространства L2, т.е. u ¯ 2= i u i 2 - сумма квадратов компонент u ¯ ).

Второе условие оптимальности - определенные ограничения на отклонение компонент p ¯ , т.е. параметров pi, от соответствующих компонент qi исходного приближения q ¯ . В результате в алгоритме МОЛИ решается задача минимизации следующего функционала:

F ( p ¯ ) = m ¯ ( p ¯ ) s ¯ 2 + i α i ( p i q i ) 2 = Ф 1 + Ф 2 (6)

Подставляя m ¯ ( p ¯ ) из (3) и учитывая выражение для у ¯ в приведенных вначале раздела обозначениях, получим тот же функционал в другой форме:

Ф ( х ¯ ) = А ( х ¯ ) у ¯ 2 + i α i х i 2 = Ф 1 ' + Ф 2 ' (7)

Минимизируя Ф( x ¯ ), мы находим корректирующий вектор x ¯ , дающий оптимальное решение p ¯ = q ¯ + x ¯ . Для нахождения минимума функционала Ф( х ¯ ) составляем систему уравнений:

Ф ( х ¯ ) х 1 = 0 , Ф ( х ¯ ) х 2 = 0 , …, Ф ( х ¯ ) х M = 0 ,

которая в матричной форме имеет вид:

[ A Т A + E α ¯ ] х ¯ = А Т у ¯ , (8)

где Е α ¯ - диагональная матрица, в которой отличны от нуля лишь диагональные элементы ( α 1 , α 2 , …, α M ).

Расчет демпфирующих множителей α i

Основной особенностью модифицированного алгоритма инверсии является способ задания демпфирующих множителей α i , задачей которого является не позволить параметрам р i оптимальных моделей выйти далеко за границы установленных диапазонов. Как следует из (6), (7), величины α i регулируют степень свободы изменения q i , т.к. чем больше α i , тем сильнее вклад корректировки q i , т.е. величины x i 2 = ( p i x i ) 2 , в функционалы, которые необходимо минимизировать. Величины α i выбираются в МОЛИ, исходя из условий:

А. Оба слагаемых в функционалах (6), (7) должны быть примерно одинаковы по величине.

В. Величины α i должны быть тем больше, чем ближе q i к заданным ограничениям q г р .

Предварительно заметим, что в (6) и (7) Ф 1 ' 1 и Ф 2 ' 2, т.к. (7) получено из (6) путем простой замены буквенных обозначений.

Для того чтобы определить α i , постараемся вначале определить величину Ф 1 = m ¯ ( p ¯ ) s ¯ 2 . Она заранее не известна, т.к. p ¯ - это модель, которая получится после выполнения соответствующей итерации в процессе оптимизации. Однако можно предположить, что энергия ошибки Ф1 после оптимизации будут представлять некоторую часть Cn(Cn<1) энергии начальной ошибки m ¯ ( q ¯ ) s ¯ 2 , где n - номер итерации. Предполагая, что амплитуды ошибки (точнее, помехи n ¯ в формуле (1)) в результате оптимизации в среднем уменьшатся в 2 раза на первом шаге итерации, задаем Сn=0,25. Поскольку на последующих шагах итерации, по мере приближения к истинной модели, уменьшение ошибки будет замедляться, задаем следующую последовательность величин Сn:

С1=0,25; С2=0,5; С3=0,7. (9)

Тогда величину Ф1 можно оценить формулой:

Ф 1 = С n m ¯ ( q ¯ ) s ¯ 2 = C n y ¯ 2 . (10)

В предыдущем разделе говорилось о том, что оптимизация выполняется на каждом шаге итерации последовательно: вначале по вектору t ¯ , затем по вектору I p ¯ . Соответственно, по-разному оцениваются величины α i . Описание способов получения этих оценок начнем с процесса оптимизации акустических параметров I p ¯ (далее для упрощения обозначений - величина q i ).

Вначале рассмотрим частный случай, когда значения х i корректирующего вектора x ¯ одинаковы по модулю, а именно:

| x i | = x с р = 0.5 Δ q max . (11)

В данном случае Δqmax=ΔIpmax - максимальное изменение q i на 1 итерацию

Тогда, согласно условию А, требуя равенство Ф1 и Ф2 в формуле (7) и учитывая (10), получим одинаковые значения α i = α с р в виде:

α с р = С n y ¯ 2 M x с р 2 , (12)

где М - число слоев в модели.

Теперь перейдем к общему случаю, задавая α i в соответствии с условием В, ориентируясь на пример конкретного положения оптимизируемого параметра q i на оси q (см. Фиг.1).

На Фиг.1:

q i - начальное приближение параметра Iр в i-ом слое,

q г р min , max - I p min , I p max для любого, в т.ч. нереперного слоя,

Δ q i = min { ( q i q г р min ) , ( q г р max q i ) } , (13)

δ - минимальное допускаемое расстояние между q i и любым из q г р , обеспечивается путем редактирования q ¯ (см. ниже), определяется формулой

δ= 1 6 Δ q max , (14)

где Δ q max - максимальное изменение Ip на 1 итерацию.

Теперь, для выполнения условия В, зададим α i со структурой, аналогичной (12), но с учетом реального положения q i :

α i = β C n y ¯ 2 M Δ q i 2 . (15)

В выражении (15) β находится путем тестирования при обработке конкретных реальных материалов по отдельным профилям, исходя из условия обеспечения устойчивости результата при внесении вариаций в исходную модель q ¯ (β=1 ÷ 4).

Посмотрим на отдельных примерах, как меняется α i в зависимости от положения q i . При q i = q с р , Δ q i = 0.5 ( q г р max q г р min ) - мы имеем минимальное α i . В ситуации, когда Δ q i = x с р = 0.5 Δ q max (полагаем, что это наиболее вероятная ситуация), α i = α с р (для выбранного фиксированного β). При минимальном Δ q i =δ (см. (14)), Δ q i = 1 3 х с р , следовательно, α 1 будет примерно на порядок больше, чем α с р .

Рассмотрим теперь способ определения демпфирующих множителей α i при расчете корректирующих значений Δti вектора t ¯ , входящего в начальное приближение q ¯ = [ I p ¯ , t ¯ ] . Учитывая, что Δti не должен быть велик изначально (Δti ≅ Δt, Δt - интервал дискретности), полагаем, что все α i одинаковы и, по аналогии с (12), выражаются формулой:

α i = α с р = C n y ¯ 2 M θ 2 . (16)

В этом выражении θ представляет собой ожидаемый средний сдвиг границ модели при оптимизации, зависит от детальности модели, определяемой величиной τmin. Опытным путем для величины θ найдено следующее выражение

θ = [ τ min 1 3 ] + 1 , (17)

где [] означает целую часть.

Рассчитанные в результате решения уравнения (8) величины Δti (компоненты вектора x ¯ являются рациональными числами, не кратными Δt. Величины Δti ограничиваются согласно условию

| Δ t i | 0.4 τ min , (18)

после чего они добавляются в значения ti исходной модели q ¯ для получения оптимизированной модели р ¯ . Таким же образом при оптимизации Iр, рассчитанные в результате решения (8) значения ΔIpi, прибавляются к исходным значениям Ipi вектора q ¯ .

Редактирование вектора q ¯ , поступающего на оптимизацию

Учитывая итерационность оптимизированного процесса, полученные после некоторого шага оптимизации векторы р ¯ подаются вновь на вход процесса оптимизации, т.е. становятся вновь начальным приближением q ¯ . Все начальные приближения, в том числе и исходные модели, подвергаются редактированию перед процессом оптимизации.

Редактирование вектора t ¯ начального приближения q ¯

Компоненты вектора t ¯ - величины ti редактируются либо путем сдвига по оси t, либо в результате объединения пластов, при этом:

а) для реперных пластов с заданными временами tk (k=1, 2, …, K) допускается объединение только с подстилающими пластами;

б) при объединении пластов скорость нового пласта рассчитывается как средневзвешенная скоростей компонент с весами, равными временной мощности τi объединяемых пластов.

Редактирование вектора t ¯ проводится следующим образом:

- проверяется выполнение условия τi>Δt, если нет, то два (или более) соседних пластов объединяются до выполнения этого условия;

- границы всех пластов выставляются на ближайшие точки шкалы дискретности;

- пласты с временной мощностью τimin объединяются с одним из двух соседних пластов, у которого Ip ближе к Ipi (в отношении реперных пластов смотри а)).

Редактирование вектора I p ¯ начального приближения q ¯

Выполняется проверка: если I p i > I p i max δ , то задаем

I p i = I p i max δ ,

а где I p i max - граничное значение Ip для i-го пласта, δ задано формулой (14). Аналогично, если I p i < I p i min + δ , то задаем

I p i = I p i min + δ .

Редактирование вектора I p ¯ в окончательной модели p ¯

Окончательную модель р ¯ получаем после выполнения последнего шага итерации. Редактируются значения Ipi согласно условиям:

если I p i > I p i max , то задаем

I p i = I p i max ;

если I p i < I p i min , то задаем

I p i = I p i min .

Алгоритм пластовой синхронной инверсии

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

s ¯ = [ s 1 ¯ , s 2 ¯ ,..., s L ¯ ] T , (19)

где компоненты s ¯ представляют собой трассы s(t, ξ1), s(t, ξ2), …, s(t, ξL). Соответственно, развернутое выражение для начального приближения q ¯ имеет вид (М-слойная модель)

q ¯ =[Vp1, Vp2, …,VpM, Vs1, Vs2, …, VsM, ρ1, ρ2, …, ρM]=

=[q1, q2, …, qM, qM+1, qM+2, …, q2M, q2M+1, q2M+2, …, q3M]. (20)

Вопрос о структуре вектора q ¯ требует специального объяснения. Синхронная пластовая инверсия выполняется как второй этап после проведения акустической инверсии. При этом временные мощности пластов Δti в точке х профиля отсчитываются с сейсмоакустических разрезов как фиксированные параметры, не требующие оптимизации. Это объясняет отсутствие вектора t ¯ в структуре q ¯ . Теперь надо объяснить, откуда берутся значения Vp, Vs, ρ пластов в начальном приближении. Для получения приближенных оценок этих параметров исследуемая пачка пластов разбивается на несколько интервалов с близкими литофизическими свойствами (карбонатные, терригенные - преимущественно глинистые или преимущественно песчанистые). Для каждого из этих интервалов, используя данные ГИС по региону, где находится площадь исследования, определяют средние величины плотностей ρ, которые приписываются всем пластам интервала. С полученного ранее сейсмоакустического разреза отсчитывают значения Ip пластов и вычисляют Vp=Ip/ρ. Также по анализу имеющейся геолого-геофизической информации для каждого из выделенных интервалов задают величины γ=Vs/Vp, что позволяет рассчитать по значениям Vp значения Vs для каждого пласта модели.

По аналогии с акустической инверсией, далее нам требуется рассчитать модельную сейсмограмму m ¯ ( q ¯ ) ,

m ¯ ( q ¯ ) = [ m 1 ¯ ( q ¯ ) , m 2 ¯ ( q ¯ ) ,..., m L ¯ ( q ¯ ) , ] T , (21)

трассы которой m(t, ξ1), m(t, ξ2), …, m(t, ξL), рассчитываются для тех же удалений ξ, что и в соответствующих трассах реальной сейсмограммы s ¯ .

Для расчета m ¯ ( q ¯ ) потребуется форма сейсмического сигнала h(t, ξ). Для ее нахождения вначале определяет форму h(t, 0) на нулевом удалении взрыв-прибор. Для этого рассчитывают трассу коэффициентов отражений r(t, 0) по модели Ip(t), выбираемой с сейсмоакустического разреза Ip(t, x) в нужной точке х профиля. Сигнал h(t, 0) определяется по известным формулам как оптимальный винеровский фильтр, преобразующий r(t, 0) в трассу s(t, ξ1), где ξ1 - минимальное удаление в сейсмограмме s(t, ξ). Затем, по известным из обработки кинематическим поправкам для кровли Δtкр(ξ) и подошвы Δtпод(ξ) исследуемой пачки слоев, растягиваем сигналы h(t, 0) согласно формуле:

h ( t , ξ ) = h [ Δ t п о д ( ξ ) Δ t к р ( ξ ) t ,0 ] . (22)

Учитывая сложность и времяемкость вычислительного процесса синхронной инверсии, мощность исследуемой пачки слоев ΔНобщ ограничивают сравнительно небольшим интервалом, включающим основные продуктивные пласты. При этом ΔНобщ обычно на порядок меньше глубины Н до кровли пачки. В этой ситуации можно пренебречь преломлением лучей внутри пачки и считать, что угол падения θ одинаков для всех плоскопараллельных границ пластов, слагающих исследуемую пачку. При этом реальное двойное время пробега волн в i-ом пласте выражается формулой:

Δ t i ( θ ) = Δ t i ( 0 ) cos θ , (23)

где Δti(0) - временная мощность пласта в модели Ip(t), выбираемой из сейсмоакустического разреза. Как сказано в разделе 1.1, в исходную сейсмограмму s ¯ при обработке вводились статические и кинематические поправки, приводящие годографы отражений от каждой i-ой границе к константе: ti(ξ)=const=ti. Это условие должно выполняться и для модельной сейсмограммы m ¯ ( q ¯ ) . При этом Δti(θ) заменяются на Δti(0), тогда используемую для расчета m ¯ ( q ¯ ) импульсную трассу на удалении взрыв-прибор ξ можно определить выражением:

r ( t , θ ) = i r i ( θ ) δ ( t t i ) , (24)

где ri(θ) - коэффициент отражения от кровли i-го пласта при падении луча под углом θ,

ti - время прихода отражения от кровли i-го пласта в модели Ip(t),

δ(t) - единичный импульс.

При обычных для практики небольших углах наклона β отражающих границ (до 10°) угол θ определяется по формуле

θ ( ξ ) = a r c t g ξ t 1 V с р , (25)

где Vср - средняя скорость до кровли пачки, определяется в процессе предшествующих структурных построений.

Коэффициент отражения Ri(θ) описывается формулой

R i ( θ ) = 1 2 ( Δ V p V ¯ p + Δ ρ ρ ¯ ) + ( Δ V p 2 V ¯ p 4 V s 2 V p 2 Δ V s V ¯ s 2 V s 2 V p 2 Δ ρ ρ ¯ ) sin 2 θ + 1 2 Δ V p V ¯ p sin 2 θ tan 2 θ (26)

где Δ V = V ( t i ) V ( t i 1 ) , V ¯ = 1 2 ( V ( t i ) + V ( t i 1 ) )

После этих промежуточных вычислений, трассы m(t, ξ) - компоненты вектора m ¯ ( q ¯ ) (21), рассчитываются по формуле:

m(t, ξ)=h(t, ξ) * r[t, θ(ξ)]. (27)

Введем обозначения:

y ¯ = s ¯ m ¯ ( q ¯ ) - разность между исходной сейсмограммой и модельной сейсмограммой для параметров исходной модели q ¯ ,

x ¯ = p ¯ q ¯ - корректирующий вектор, приводящий q ¯ к оптимальному решению p ¯ = q ¯ + x ¯ .

Применяя в случае синхронной инверсии тот же метод обобщенной линейной инверсии, как и для акустической инверсии, корректирующий вектор x ¯ находится в результате решения системы уравнений (8), имеющий вид:

[ A Т A + E α ¯ ] х ¯ = А Т у ¯ (28)

В данном случае матрица А представляется в виде:

А = [ m 1 ¯ ( q ¯ ) q 1 m 1 ¯ ( q ¯ ) q 2 . . . m 1 ¯ ( q ¯ ) q 3 M m 2 ¯ ( q ¯ ) q 1 m 2 ¯ ( q ¯ ) q 2 . . . m 2 ¯ ( q ¯ ) q 3 M . . . . . . . . . . . . m L ¯ ( q ¯ ) q 1 m L ¯ ( q ¯ ) q 2 . . . m L ¯ ( q ¯ ) q 3 M ] (29)

В этом выражении m ¯ ( q ¯ ) q i - частная производная -ой трассы сейсмограммы m ¯ ( q ¯ ) по i-ому параметру вектора q ¯ (20). По аналогии со случаем акустической инверсии, она представляется в виде:

m ¯ ( q ¯ ) q i = Δ r i ' [ θ ( ξ ) ] h ( t t i ' , ξ ) + Δ r i ' + 1 [ θ ( ξ ) ] h ( t t i ' + 1 , ξ ) Δ q i , (30)

где i' - номер слоя (i'=1, 2, …, M), соответствующего параметру qi - см. (20),

Δri', Δri'+1 - изменения коэффициентов отражения от кровли и подошвы слоя i' при замене qi на qi+Δqi,

ti', ti'+1 - положение на оси t кровли и подошвы слоя i' в модели Ip(t), полученной в результате акустической инверсии.

Для составления уравнения (28), так же, как и в случае акустической инверсии, нужно задать демпфирующие множители α i . Они определяются точно таким же методом, как и для акустической инверсии. По аналогии с (15), в данном случае выражение для α i можно представить в виде:

α i = β C n y ¯ 2 3 M Δ q i , (31)

где Δ q i = min ( | q i q i min | , | q i q i max | ) , а граничные значения для qi заданы в исходной информации.

Способ определения упругих свойств горных пород на основе пластовой адаптивной инверсии сейсмических данных, характеризующийся применением пластовых моделей среды, в которых минимальные временные мощности τmin пластов соответствуют реальной разрешающей способности сейсморазведки и геологии осадконакопления и вычисляются согласно формуле:
,
где
Δf - рабочая полоса частот;
причем в качестве исходной модели используют откорректированную геоакустическую модель импедансов соответствующей скважины, либо трассу импедансов полученного ранее сейсмоакустического разреза по секущему профилю; а при построении рабочей модели преобразуют сейсмический временной разрез в детальную пластовую модель акустических импедансов (сейсмоакустический разрез), по которой свидетельствуют о диапазонах изменения акустических параметров слоев.



 

Похожие патенты:

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

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

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

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

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

Изобретение относится к области измерительной техники, в частности к средствам мониторинга технического состояния различных сооружений, и может быть использовано для текущей оценки и прогноза безопасной эксплуатации зданий и/или сооружений при возможных неблагоприятных воздействиях на объект.

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

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

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

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

Изобретение относится к области геофизики и может быть использовано при проведении акустического каротажа при бурении подземных формаций. Способ проведения измерений акустического каротажа включает группирование полученных форм акустических сигналов в одну из множества групп. При этом каждая такая репрезентативная группа соответствует некоторым измеренным параметрам состояния буровой скважины (например, диапазон измеренных значений отклонения и/или диапазон измеренных азимутальных углов). Формы акустических сигналов, сохраненные, по меньшей мере, в одной из групп, накладываются одна на другую для получения усредненной формы сигнала. Впоследствии такая усредненная форма сигнала может подвергаться обработке, например, с использованием алгоритма определения меры когерентности для получения, по меньшей мере, одного значения замедления акустической волны. Технический результат - повышение точности каротажных данных. 4 н. и 13 з.п. ф-лы, 4 ил.

Изобретение относится к области геофизики и может быть использовано для получения сейсмических разрезов изображений геологической среды. Способ включает последовательные действия, при которых получают и подготавливают данные методов общей глубинной точки, сейсмического каротажа, вертикального сейсмического профилирования, акустического каротажа, плотностного гамма-гамма каротажа и проверяют качество этих данных, а также получают эталонные значения интервальных скоростей. Получают исходный годограф и рассчитывают синтетическую сейсмограмму. Затем проводят контроль качества и вводят постоянную временную поправку для посадки на верхний опорный горизонт литолого-стратиграфического комплекса. Затем вновь рассчитывают синтетическую сейсмограмму и вновь проводят контроль качества. Вслед за этим рассчитывают и вводят поправку для посадки на нижний опорный горизонт литолого-стратиграфического комплекса. После этого вновь рассчитывают синтетическую сейсмограмму и осуществляют контроль качества. Переносят точки полученного годографа на ближайшие акустически слабые границы. Повторно рассчитывают синтетическую сейсмограмму с последующим контролем качества и получают априорный годограф. Технический результат - повышение достоверности и точности соответствия горизонтов временного разреза и геологических отметок скважины. 10 з.п. ф-лы, 2 ил.

Изобретение относится к области геофизики и может быть использовано для обработки данных сейсморазведки. Заявлен способ преобразования сейсмических данных для получения модели объемного модуля упругости или плотности геологической среды. Градиент целевой функции вычисляют (103), используя сейсмические данные (101) и модель (102) вмещающей геологической среды. Коэффициент освещения источником и коэффициент освещения приемника вычисляют (104) в модели вмещающей среды. Объем сейсмического разрешения вычисляют (105), используя скорости из модели вмещающей среды. Градиент преобразуют (106) в разностные параметры модели геологической среды, используя коэффициент освещения источником, коэффициент освещения приемника, объем сейсмического разрешения и модель вмещающей геологической среды. Указанные величины являются масштабными коэффициентами, используемыми для компенсации сейсмических данных после миграции путем обратной миграции во временной области, которые затем можно связать с моделью объемного модуля упругости геологической среды. В случае итерационной инверсии разностные параметры (106) модели геологической среды используют (107) в качестве предобусловленных градиентов. Технический результат - повышение точности получаемых данных. 20 з.п. ф-лы, 7 ил.

Изобретение относится к области геофизики и может быть использовано при обработке сейсмических данных при поиске месторождений углеводородов. Заявленный способ идентификации геологических особенностей из геофизических или атрибутивных данных предполагает использование выполняемого в окне анализа главных компонент или анализа независимых компонент, либо диффузионного картирования. Едва уловимые особенности становятся идентифицируемыми в частичных или остаточных массивах данных. Остаточные массивы данных создаются путем исключения данных, не захваченных самыми заметными главными компонентами. Частичные массивы данных создаются путем проецирования данных на выбранные главные компоненты. Геологические особенности также можно идентифицировать из анализа образов или массивов аномалий, сформированных с помощью матрицы подобия данных с переменным масштабом . Технический результат - повышение точности данных прогноза нефтегазоносности исследуемого участка. 4 н. и 27 з.п. ф-лы, 11 ил.

Изобретение относится к области геофизики и может быть использовано при обработке данных сейсмических исследований. Заявлен способ перестроения моделей (110) Q геологической среды на основании сейсмических данных (10) путем осуществления лучевой Q томографии сдвига центроидных частот. Амплитудный спектр волнового сигнала сейсмического источника аппроксимируют (40) частотно-взвешенной экспоненциальной функцией частоты, имеющей два подбираемых параметра для приведения в соответствие данным о сдвиге частот. В результате чего обеспечивают лучшее соответствие различным асимметричным амплитудным спектрам источника. Боксовые ограничения могут использоваться при выполнении процедуры оптимизации, а многоиндексный способ активных множеств, используемый при томографии скорости, является предпочтительным способом для реализации (100) боксовых ограничений. Технический результат - повышение точности данных сейсмических исследований. 2 н. и 18 з.п. ф-лы, 12 ил.

Изобретение относится к области сейсморазведки и может быть использовано для поиска и разведки углеводородов (УВ). Согласно способу оценки низкочастотной резонансной эмиссии (НРЭ) для поиска УВ прогнозирование УВ осуществляется в процессе анализа геодинамического шума непосредственно по временному разрезу метода общей глубинной точки (МОГТ) в широком диапазоне частот (5-130 Гц). Технический результат - расширение функциональных возможностей при повышении информативности, точности и достоверности данных по поиску и разведке УВ непосредственно по временному разрезу МОГТ. 4 ил.

Изобретение относится к сейсмической разведке и может использоваться при разведке нефтяных и газовых залежей. Согласно заявленному решению выбирают и устанавливают фиксированную приемную базу, располагают источники возбуждения сейсмических колебаний и приемники на этой базе симметрично относительно ее центра, принятого за начало координат. Поочередно возбуждают сейсмические волны и регистрируют их на всей заданной фиксированной приемной базе с получением сейсмограммы для каждого пункта возбуждения (ПВ). Определяют время годографа дифрагированных волн (ДВ). Затем суммируют на каждой сейсмограмме сейсмические волны по годографам ДВ с получением суммотрассы для каждого ПВ(L) и определяют временные статические сдвиги (Δtc) ДВ для каждого ПВ(L) относительно времен регистрации на центральном ПВ(L=0). Суммируют на основании временных статических сдвигов суммотрассы для каждого ПВ(L) по общим точкам дифракции и получают итоговую суммотрассу ДВ. Затем перемещают фиксированную приемную базу на следующий пикет профиля, совмещая его с началом координат этой базы, и весь процесс повторяют. Технический результат - повышение достоверности и качества отображения изучаемой геологической среды. 6 ил.

Изобретение относится к области геофизики и может быть использовано в процессе проведения микросейсмического мониторинга. Настоящее изобретение предусматривает волоконно-оптическую систему сейсмического мониторинга, включающую в себя источник света, который возбуждает оптическое волокно, расположенное в стволе скважины. По меньшей мере одним датчиком света анализируется свет рэлеевского обратного рассеяния для получения акустических сигналов для каждой из множества точек на всем протяжении ствола скважины. Множество точек измерения связаны с различными подмассивами, включающими в себя по меньшей мере один подмассив перед кривой вдоль траектории ствола скважины и по меньшей мере один подмассив после кривой. Процессор определяет угол падения для каждого из различных подмассивов и находит посредством триангуляции местоположение микросейсмического события на основании определенных углов падения. Технический результат - повышение точности данных исследований. 2 н. и 14 з.п. ф-лы, 13 ил.

Изобретение относится к области геофизики и может быть использовано при проведении сейсморазведочных работ. Измеренные сейсмические данные принимаются от сейсмического датчика. Также принимаются данные вращения. Данные вращения представляют вращение относительно, по меньшей мере, одной определенной оси. Причем данные вращения относительно определенной оси обеспечивают получение опорного шумового сигнала для поверхностных шумовых волн, распространяющихся по поверхности земли, на которой установлен сейсмический датчик. Далее осуществляют комбинирование, используя адаптивную фильтрацию, данных вращения с измеренными сейсмическими данными с целью понижения, по меньшей мере, части компоненты шума, содержащей поверхностные шумовые волны в измеренных сейсмических данных. Адаптивная фильтрация содержит адаптивное вычитание опорного шумового сигнала, полученного с использованием данных вращения, из сейсмических данных. Технический результат - повышение точности получаемых данных. 3 н. и 20 з.п. ф-лы, 6 ил.

Изобретение относится к области сейсмической разведки и может быть использовано при поиске нефтяных и газовых месторождений со сложно построенными кавернозно-трещиновато-пористыми коллекторами. Заявлен способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения на основе выделения и специализированной обработки рассеянной компоненты сейсмического волнового поля с привлечением метода декомпозиции эмпирических мод и полномасштабного численного моделирования в целях формирования высокоточных и достоверных поисковых признаков флюидонасыщения карбонатного коллектора. Технический результат - повышение точности и достоверности получаемых данных. 7 ил.
Наверх