Способ определения трещинной пористости пород

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

 

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

Известен способ оценки трещинной пористости в шлифах под микроскопом (Гмид Л.П. Методическое руководство по литолого-петрографическому и петрохимическому изучению осадочных пород-коллекторов. СПб: ВНИГРИ, 2009, с.160), в котором измеряется площадь шлифа, длина следов трещин и их ширина. По данным замеров производится подсчет параметров трещиноватости по формулам, полученным экспериментальным путем на опытах в открытых щелях, имитирующих трещины разной раскрытости и ориентировки. Недостатком указанного способа является его низкая точность и достоверность полученных данных, обусловленная тем, что исследуется только малый объем горной породы в шлифе.

Наиболее близким к предложенному способу (прототипом) является способ определения трещинной пористости пород (патент РФ №2012021, G01V 1/40, опубл. 30.04.1994), заключающийся в проведении в изучаемом разрезе волнового акустического и гамма-гамма каротажа. По данным каротажа определяют коэффициент сжимаемости пород для двух значений плотности заполняющей скважину промывочной жидкости. При этом плотность увеличивают на 15-20% в зависимости от глубины скважины и допустимой величины давления гидроразрыва пород. Коэффициент пористости пород определяют с учетом коэффициентов сжимаемости пород, определенных по двум замерам, коэффициента сжимаемости матрицы (блока) и изменения плотности бурового раствора перед повторным исследованием. Недостатком известного способа является невысокая точность, обусловленная отсутствием достоверных данных о коэффициенте сжимаемости матрицы и методов его определения в реальных условиях залегания пород, а также отсутствием надежных данных о зависимости коэффициента сжимаемости пласта от изменений плотности бурового раствора в скважине.

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

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

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

Кп тр=[100-1,6Кп общ - 100(Vp изм/Vp ск)]/20,4,

где Кп общ - экспериментально определенная общая пористость образца;

Vp изм - измеренная скорость распространения продольной волны в образце;

Vp ск - скорость распространения продольной волны в минеральном скелете исследуемой породы,

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

В горной породе поры и трещины образуют общую пористость

К п о б щ = К п п о р + К п т р , ( 1 )

где Кп пор - поровая пористость горной породы, %;

Кп тр - трещинная пористость горной породы, %.

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

Отношение измеренной скорости распространения продольных волн в образце исследуемой породы к скорости распространения продольных волн в минеральном скелете исследуемой породы (при Кп общ, равной нулю), выраженное в процентах, называется добротностью Q (Мори В. Механика горных пород применительно к проблемам разведки и добычи нефти. М.: Мир, 1994, с.176-184) и характеризует воздействие пор и трещин на породу

Q = 100 ( V р и з м / V р с к ) , ( 2 )

где Vp изм - измеренная скорость распространения продольной волны в образце исследуемой породы, км/с;

Vp ск - скорость распространения продольной волны в минеральном скелете исследуемой породы, км/с.

Как вытекает из выражения (2), при добротности, равной 100%, порода не имеет ни трещин, ни пор. Уменьшение значения добротности отражает наличие в горной породе трещин и пор.

Известна зависимость добротности Q от поровой пористости Кп пор

Q = 100 1,6 К п п о р , ( 3 )

где Q - добротность горной породы, %;

а также - зависимость добротности от трещинной пористости горной породы:

Q = 100 22 К п т р , ( 4 )

откуда вытекает зависимость добротности от общей пористости:

Q = 100 1,6 К п п о р 22 К п т р , ( 5 )

решая известное уравнение (5) относительно Кп тр, получаем формулу:

К п т р = ( 1 0 0 - 1 , 6 К п п о р - Q ) / 22. ( 6 )

Подставляя в формулу (6) выражение для Q по формуле (2) и выражение для Кп пор по формуле (1), получаем конечную формулу для вычисления трещинной пористости:

К п т р = [ ( 100 1,6 К п о б щ 100 ( V р и з м / V р с к ) ] / 20,4. ( 7 )

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

На фиг.2 - зависимость скорости распространения продольной волны от общей пористости для образцов исследуемой породы.

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

Способ осуществляют следующим образом.

- Формируют набор образцов исследуемой породы.

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

- Определяют общую пористость каждого из образцов в условиях, моделирующих пластовые. Общую пористость определяют посредством измерений объема жидкости, вытесненной из порового пространства образца при увеличении эффективного давления от 0,1 МПа до давления в пласте (обычно более 15 МПа), и с учетом объема образца по формуле

К п о б щ = К п о б щ а т м - ( Δ V п о р / V о б р ) 100, ( 8 )

где Кп общ пл - общая пористость образца в условиях, моделирующих пластовые, %;

Кп общ атм - общая пористость образца в атмосферных условиях, %;

ΔVпор - объем жидкости, вытесненный из порового пространства образца при переходе от атмосферных условий к условиям, моделирующим пластовые, см3;

Vобр - объем образца, см3.

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

Общую пористость (Кп общ пл) и скорость распространения продольной волны (Vp пл) в условиях, моделирующих пластовые, определяют экспериментально с помощью установки ПУМА-650 или «Экопласт - Гео», или другой подобной установки, позволяющей моделировать пластовые условия и определять общую пористость и скорость распространения упругой продольной волны в исследуемой породе. На установке моделирования пластовых условий изменяют напряженное состояние образцов исследуемой породы путем создания всестороннего давления Рве, равного литостатическому давлению, и порового давления Рпор, равного давлению флюида (газ, жидкость) в пласте. При этом создают эффективное давление Рэф, равное их разности. При достаточно больших значениях эффективного давления (40,0 МПа и более) скорость распространения продольной волны в образцах достигает максимума (фиг.1).

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

V р и з м = А К п о б щ + V р с к , ( 9 )

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

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

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

Пример осуществления способа.

- Сформировали набор из 22 образцов песчаника.

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

- С помощью установки ПУМА-650 определили общую пористость и скорость распространения продольной волны для каждого из образцов в термобарических условиях, моделирующих пластовые (всестороннее давление Рвс=50 МПа, поровое давление Рпор=13 МПа, температура Т=22°С), при этом точность определения: пористости - ±0,01%, скорости распространения продольной волны - ±0,002 км/с.

- Определили скорость распространения продольной волны в минеральном скелете исследуемого песчаника (Vp ск), используя линейную зависимость (9) (при А=0,11) и полученные экспериментально значения скорости распространения продольной волны в образцах и общей пористости всех образцов в условиях, моделирующих пластовые (фиг.3). Полученная Vp ск=5,8824 км/с.

- Определили величину трещинной пористости (Кп тр) для каждого из образцов исследуемого песчаника, используя формулу (7).

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

Полученные экспериментальные и расчетные данные приведены в табл.1.

Точность определения трещинной пористости исследуемых образцов определяется погрешностью, допущенной при определении общей пористости в условиях, моделирующих пластовые, и составляет±0,01%.

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

Таблица 1
Кп общ атм, % Кп общ пл, % Vp пл, км/с Кп тр пл, % Кп пор пл, %
1 8.9 7.95 4.65 0.41 7.55
2 7.5 6.69 5.18 0.06 6.63
3 6.9 6.45 5.15 0.11 6.34
4 10.9 10.20 4.53 0.33 9.86
5 11.1 10.42 4.90 0.01 10.41
6 11.1 10.31 4.74 0.14 10.17
7 12.9 12.20 4.46 0.23 11.96
8 10.3 9.75 4.90 0.06 9.69
9 10.5 9.66 4.82 0.13 9.53
10 12.8 12.15 4.53 0.17 11.98
11 12.6 11.98 4.60 0.13 11.85
12 11.9 11.22 4.67 0.13 11.09
13 10.5 9.88 4.81 0.12 9.76
14 8.4 7.77 5.05 0.08 7.68
15 6.9 6.32 5.23 0.05 6.27
16 5.5 4.86 5.33 0.08 4.77
17 11.2 10.74 4.73 0.11 10.62
18 10.1 9.61 4.89 0.07 9.53
19 10.0 9.39 4.89 0.09 9.30
20 10.1 9.52 4.73 0.21 9.31
21 10.7 10.17 4.82 0.09 10.08
22 4.7 4.09 5.51 -0.01 4.11

Способ определения трещинной пористости пород, включающий измерение скорости распространения продольной волны в исследуемой породе и выполнение расчета с использованием полученных данных, отличающийся тем, что предварительно формируют набор образцов исследуемой породы, экспериментально определяют общую пористость каждого из упомянутых образцов в атмосферных условиях, также экспериментально определяют скорость распространения продольной волны и общую пористость в образцах исследуемой породы в условиях, моделирующих пластовые условия, после чего определяют скорость распространения продольной волны в минеральном скелете исследуемой породы с использованием полученной экспериментально зависимости скорости распространения продольной волны в образцах исследуемой породы от их общей пористости, определенных в условиях, моделирующих пластовые условия, далее рассчитывают величину трещинной пористости (Кп тр) для каждого из образцов исследуемой породы по формуле
Кп тр=[100-1,6Кп общ - 100(Vp изм/Vp ск)]/20,4,
где Кп общ - экспериментально определенная общая пористость образца;
Vp изм - измеренная скорость распространения упругой продольной волны в образце;
Vp ск - скорость распространения продольной волны в минеральном скелете исследуемой породы,
после чего определяют поровую пористость, как разницу между общей пористостью и трещинной пористостью.



 

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

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

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

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

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

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

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

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

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

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

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

Изобретение относится к области геофизики и может быть использовано при проведении сейсморазведочных работ. Согласно заявленному способу проводится попарное непрерывное сопоставление множества трасс сейсмического разреза или куба. Для каждой пары сейсмических трасс разреза или каждой пары трасс всех ин-лайнов и всех кросс-лайнов вычисляется двумерная функция взаимного различия фрагментов сейсмических записей трасс. Значения этой функции вычисляются для всех возможных сочетаний времени (глубины) данной пары трасс, которые могут принадлежать одному сейсмостратиграфическому уровню. Вычисляют функции оптимального соответствия времени/глубины от i-той к i+k-той сейсмической трассе t i + k = φ + k ( t i ) и t i = φ − k ( t i + k ) , а затем вычисляют функции перехода для временной или глубинной области всех соседних сейсмических трасс. Данные функции перехода позволяют создать единую непрерывную двумерную для разреза и трехмерную для куба сейсмостратиграфическую модель. Технический результат - повышение точности и достоверности данных картирования горизонтов и восстановления параметров геологической среды. 6 ил.

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

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

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

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

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