RU2758416C1 - Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные - Google Patents
Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные Download PDFInfo
- Publication number
- RU2758416C1 RU2758416C1 RU2020131555A RU2020131555A RU2758416C1 RU 2758416 C1 RU2758416 C1 RU 2758416C1 RU 2020131555 A RU2020131555 A RU 2020131555A RU 2020131555 A RU2020131555 A RU 2020131555A RU 2758416 C1 RU2758416 C1 RU 2758416C1
- Authority
- RU
- Russia
- Prior art keywords
- azimuth
- cavernous
- fractured
- seismic
- waves
- Prior art date
Links
- 230000004069 differentiation Effects 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 title claims description 43
- 238000010276 construction Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000004807 localization Effects 0.000 claims description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 3
- 239000003208 petroleum Substances 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 7
- 230000035508 accumulation Effects 0.000 description 6
- 238000009825 accumulation Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 230000005284 excitation Effects 0.000 description 4
- 230000002547 anomalous effect Effects 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- 150000002430 hydrocarbons Chemical class 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 150000001491 aromatic compounds Chemical class 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 125000003118 aryl group Chemical group 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 150000004649 carbonic acid derivatives Chemical class 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Изобретение относится к области сейсмической разведки и может быть использовано при поиске нефтяных и газовых месторождений со сложно построенными карбонатными трещиновато-кавернозными коллекторами. По стандартным сейсмическим изображениям среды, полученным с использованием отраженных волн, проводят анализ поведения отражающих горизонтов на временных разрезах. Определяют углы наклона границ целевых горизонтов, по которым рассчитывают оптимальные параметры асимметричного суммирования сейсмических данных - угол и азимут трассирования Гауссовых пучков, которые определяют весовые множители и апертуру, обеспечивающие максимально надежное построение поля энергии рассеянных волн. С полученными параметрами проводят асимметричное суммирование сейсмограмм после сигнальной обработки и куба средних скоростей. Полученные селективные изображения в области углов наклона и азимутов трассирования суммируют и из полученного полноазимутального поля энергии рассеянных волн проводят расчет атрибутов рассеяния, основанных на анализе селективных изображений в зависимости от азимута и типа дифракции/рассеяния, и на основании имеющейся на площади фактической геолого-геофизической информации осуществляют реконструкцию тонкой структуры геологических объектов и их дифференциацию на зоны кавернозности и трещиноватости. Технический результат - повышение точности и детальности реконструкции структуры карбонатных коллекторов при выявлении и локализации зон трещиноватости и кавернозности, а также возможность дифференцирования этих объектов, различая скопления трещин и каверн. 5 ил.
Description
Изобретение относится к области сейсмической разведки и может быть использовано при поиске нефтяных и газовых месторождений со сложно построенными карбонатными трещиновато-кавернозными коллекторами.
В настоящее время при интерпретации данных сейсморазведочных работ используются в основном только регулярные отраженные волны, обеспечивающих построение как скоростной модели, так и геометрии границ раздела и реконструкции свойств, связанных с ними слоев горных пород. Однако фундаментальные физические свойства сейсмических волновых полей накладывают весьма жесткие ограничения на разрешающую способность методов, опирающихся на использование отраженных волн. В то же время, уже само по себе присутствие рассеянных волн несет информацию о наличии в среде микроструктуры, характерные размеры которой существенно меньше доминирующей длины волны. Следовательно, использование таких волн открывает возможность значительного повышения информативности и разрешающей способности сейсмических методов изучения внутреннего строения среды.
Известен способ исследования вертикальных зон трещиноватости, состоящий в возбуждении и последующей регистрации сейсмических колебаний в горных породах исследуемой зоны с целью выделения из них дуплексных волн и использование их амплитуды для вынесения суждения о присутствии вертикальной границы, разделяющей среды с различными свойствами с использованием дополнительных сейсмических измерений с помощью установки, содержащей, по меньшей мере, две пары «источник -приемник», размещенные на одинаковом расстоянии с противоположных сторон выявленной вертикальной границы, измерение расстояния от каждого источника и каждого приемника до выявленной вертикальной границы, выделение дуплексных волн, отраженных от точек, расположенных на выявленной вертикальной границе, определение их скорости и амплитуды, выявление зоны, ограниченной двумя вертикальными границами, определение конфигурации указанной зоны с последующим ее картированием, отнесение выделенной зоны к зоне трещиноватости при расстоянии между двумя вертикальными границами, ее ограничивающими, меньшем длины дуплексной волны, оценку степени трещиноватости в выявленной зоне путем сравнения амплитуд дуплексных волн, зарегистрированных при дополнительных сейсмических измерениях, с их эталонными значениями. Регистрацию дуплексных волн при основных и дополнительных сейсмических измерениях проводят послойно и раздельно для слоев, характеризующихся различными интенсивностью и преобладающими направлениями трещиноватости. Шаг послойного измерения дуплексных волн не менее их длины. RU 2415448 C1, опубл. 27.03.2011.
Недостаток способа состоит в чрезвычайно узком спектре применения, ограниченном весьма частными случаями протяженной вертикальной трещиноватости. Кроме того, при этом требуется регулярность таких трещин в вертикальном направлении, то есть отсутствие шероховатости берегов трещин. Дело в том, что только в этом случае возникает дуплексная волна, то есть волна, распространяющаяся вдоль вертикальной или субвертикальной трещины. К сожалению, вертикальной трещиноватости, как правило, не присуща выдержанность по вертикали. Длина регулярных трещин составляет, как правило, первые метры, другими словами редко превышает десятые доли доминирующей длины волны. Таким образом, область применения метода дуплексных волн ограничена геологическими средами с протяженными вертикальными системами трещин и неприменим к изучению мелкомасштабных неоднородностей карбонатных резервуаров.
Известен способ прямого прогноза залежей углеводородов, который основан на выделении аномалий микросейсмической активности и концентраций химических соединений ароматического ряда углеводородов по одним и тем же профилям. Проводится сравнение аномальных распределений между собой и по участкам совпадения аномалий микросейсмической эмиссии и концентраций ароматических соединений со структурными особенностями на сейсмических временных разрезах, производится прогноз пространственного положения залежи. Для прогноза проводится дообработка данных сейсморазведки путем формирования поля случайных волн не только до вступления регулярных волн, но и в последующих частях сейсмограмм, полученных по технологии общей глубинной точки, где время регистрации превышает появление отраженных волн от глубинных частей разреза, путем создания совокупности трасс равных удалений от пунктов возбуждения и сводных временных полей, объединенных единым пунктом приема, выделяют участки аномальных значений энергии сейсмической эмиссии, нормированной на энергию общего поля случайных колебаний. Согласно изобретению, волновое поле случайных волн формируют на основе сейсмограмм, полученных при профильных и площадных наблюдениях отраженных волн, путем создания серии временных полей, составленных из совокупности трасс равных удалений от пункта взрыва по участкам, свободным от записи регулярных волн, и сводных временных полей, составленных из сейсмограмм многократных перекрытий, объединенных единым пунктом приема. Выделяют участки аномальных значений энергии низкочастотной составляющей, определяют ее долю в общем поле микросейсм, рассчитывают амплитудно-частотные спектры и их градиенты в низкочастотной области. Далее строят карты распределения аномальных концентраций ароматических соединений по площади. Наносят на эти карты распределение аномалий микросейсмической активности и по результатам совпадения аномалий прогнозируют наличие залежи и дают экспертную оценку. RU 2454687 C1, опубл. 27.06.2012.
Недостатком этого решения является необходимость проведения геохимических наблюдений.
Известен способ комплексной обработки геофизических данных, который включает последовательное накопление измерительной информации от измерителей параметров геофизических полей, обработку измеренных данных, а также анализ и интерпретацию данных, причем измеренные данные обрабатывают последовательно в несколько этапов. Накопление информации осуществляют в базе данных (БД) априорной геологической информации, в БД акустического каротажа (АК) опорных скважин и в БД сейсмограмм общей глубинной точки (ОГТ) 2D/3D в окрестности опорных скважин. Данные обрабатывают последовательно в семь основных этапов: на первом этапе обрабатывают данные АК и метода ОГТ и формируют априорные скорости ОГТ, на втором этапе формируют горизонты отражающих границ, поля скоростей ОГТ и временные поля сейсмограмм метода ОГТ, на третьем этапе формируют детальное поле скоростей ОГТ с увеличенной латеральной разрешенностью, на четвертом этапе формируют среднеслоистую модель скоростей упругих волн, на пятом этапе формируют тонкослоистую модель скоростей упругих волн, на шестом этапе формируют тонкослоистую модель упругих параметров, на седьмом этапе формируют тонкослоистые модели фильтрационно-емкостных свойств и тип флюидонасыщения в межскважинном и околоскважинном пространстве. По данным обработки измерений на первом - седьмом этапах проводят анализ и комплексную интерпретацию совокупности данных с вынесением суждения о наличии объектов углеводородов малой мощности (менее 15-20 м), целесообразности их разработки, мониторинга и оптимизации размещения эксплуатационных скважин на исследованной площади. Отличием способа также является то, что БД априорной геологической информации и БД акустического каротажа АК опорных скважин формируют в виде баз данных, которые содержат данные геофизического исследования скважин (ГИС), причем БД априорной геологической информации содержит данные литографических колонок опорных скважин, стратиграфические разбивки и данные лабораторного анализа керна, а БД АК опорных скважин содержит данные измерений акустического каротажа АК, данные гамма-гамма плотностного каротажа и данные кавернометрии, причем проводят корректировку данных ГИС с возможностью статистического формирования поправок с использованием петрофизических зависимостей. RU 2490677 C2, опубл. 20.08.2013.
К недостаткам способа следует отнести привлечение стандартных методов повышения информативности и разрешающей способности сейсмических данных на основе изучения корреляционных зависимостей, которые не гарантируют корректную локализацию объектов.
Наиболее близким к изобретению относится способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения, основанный на использовании численного моделирования трехмерных волновых полей с целью калибровки результатов выделения рассеянной составляющей. Само выделение производится путем асимметричного суммирования данных многократного перекрытия в пределах апертур источников и приемников. Принцип асимметричного суммирования данных является общей чертой прототипа и предлагаемого способа. RU 2563323 C1, опубл. 20.09.2015.
Главными недостатками прототипа является недостаточная разрешающая способность и отсутствие четкого алгоритма для определения апертур источников и приемников, обеспечивающих максимально возможную разрешающую способность способа.
Технической задачей предлагаемого изобретения является разработка способа реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные, позволяющего прогнозирование коллекторских свойств геологического разреза трещиновато-кавернозных карбонатных месторождений по данным сейсморазведки.
Технический результат, достигаемый изобретением, заключается в повышении точности и детальности реконструкции структуры карбонатных коллекторов при выявлении и локализации зон трещиноватости и кавернозности, а также в возможности дифференцирования этих объектов, различая скопления трещин и каверн.
Технический результат достигается способом реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные, согласно которому, по стандартным сейсмическим изображениям среды, полученным с использованием отраженных волн, проводят анализ поведения отражающих горизонтов на временных разрезах, определяют углы наклона границ целевых горизонтов, по которым рассчитывают оптимальные параметры асимметричного суммирования сейсмических данных - угол и азимут трассирования Гауссовых пучков, которые определяют весовые множители и апертуру, обеспечивающие максимально надежное построение поля энергии рассеянных волн, с полученными параметрами проводят асимметричное суммирование сейсмограмм после сигнальной обработки и куба средних скоростей, получают селективные изображения в области углов наклона и азимутов трассирования, которые суммируют и из полученного полно азимутального поля энергии рассеянных волн проводят расчет атрибутов рассеяния, основанных на анализе селективных изображений в зависимости от азимута и типа дифракции/рассеяния, затем на основании имеющейся на площади фактической геолого-геофизической информации, осуществляют реконструкцию тонкой структуры геологических объектов и их дифференциацию на зоны кавернозности и трещиноватости.
Указанные признаки изобретения существенны.
Изобретение основано на специальной обработке данных трехмерной сейсморазведки для построения изображений в рассеянных волнах, на основе которых выполняется реконструкция тонкой структуры геологических объектов месторождения и дифференциация на трещиноватые и кавернозные объекты.
Возможность реконструкции тонкой структуры геологических объектов месторождения за счет глубокого использования рассеянных волн, дифференциация объектов на зоны кавернозности и трещиноватости, существенно отличает заявленный способ от прототипа.
Предлагаемый способ свободен от недостатков, которые присутствуют у описанных выше способов: он не ограничен геологическими средами с протяженными вертикальными системами трещин, для него отсутствует необходимость привлечения геохимических методов, не использует корреляционные зависимости, он гарантирует корректную локализацию и дифференциацию трещиноватых и кавернозных объектов, при этом обеспечивается более высокая разрешающая способность.
Основой изобретения является применение способа обработки данных трехмерной сейсморазведки в целях построения изображений в рассеянных/дифрагированных волнах. В ее основе лежит асимметричное суммирование данных многократного перекрытия с использованием специальным образом вычисляемых весовых множителей. Именно использование таких весов и позволяет обеспечить высокую разрешающую способность построенных изображений и сохранение на них истинной амплитуды отражения/рассеяния. Такие свойства определяют за счет разложения полного волнового поля по Гауссовым пучкам - наиболее остронаправленным специально построенным сейсмическим волнам. Из этого можно получать изображения заранее выбранных геологических объектов, от протяженных регулярных границ до зон скопления микронеоднородностей, таких как, зоны повышенной трещиноватости и кавернозности горных пород.
Таким образом, ключевыми шагами реализации данного изобретения являются:
1. Выделение поля энергии рассеянных волн за счет разложения полного волнового поля по Гауссовым пучкам, получение сейсмических атрибутов, основанных на анализе изображений в зависимости от азимута и типа дифракции/рассеяния.
2. Реконструкция тонкой структуры геологических объектов и их дифференциация на трещиноватые и кавернозные.
Способ иллюстрирован следующими фигурами:
На фиг. 1 представлены принципы формирования апертур асимметричного суммирования предлагаемого способа (а) и прототипа (б).
На фиг. 2 представлены изображения разрезов поля рассеянных волн, полученных предлагаемым способом (а), способом прототипа (б) и графики их спектров (в).
На фиг. 3 представлены изображения каркаса цифровой модели-двойника реального месторождения (а) и результирующий синтетический трехмерный куб (б).
На фиг. 4 приведено сравнение горизонтального сечения модели резервуара месторождения (а) и его реконструкция тонкой структуры геологических объектов в поле рассеянных волн (б).
На фиг. 5 показана дифференциация объектов резервуара модели месторождения на зоны трещиноватости и кавернозности, при помощи атрибутов дифракции/рассеяния. Способ реализуют последовательным выполнением следующих приемов:
1. Привлекают стандартные сейсмические изображения среды, полученные одним из методов обработки сейсмических данных, ориентированным на использование отраженных волн (например, реализации миграции до или после суммирования). Проводят анализ поведения отражающих горизонтов на временных разрезах, рассчитывают утлы наклона границ. Полученные оценки используют для определения оптимальных параметров асимметричного суммирования сейсмических данных: угла и азимута трассирования Гауссовых пучков, которые определяют весовые множители и апертуру, обеспечивающие максимально надежное построение поля энергии рассеянных волн.
2. С полученными параметрами проводят асимметричное суммирование сейсмограмм после сигнальной обработки и куба средних скоростей, с использованием программного обеспечения. В результате получают селективные изображения в области углов наклона и азимутов.
3. Селективные изображения суммируют по всем углам и азимутам трассирования, получают полноазимутальное поле рассеянных волн.
4. Из полноазимутального поля энергии рассеянных волн с использованием программного обеспечения выполняют расчет сейсмических атрибутов, основанных на анализе изображений в зависимости от азимута и типа дифракции/рассеяния.
5. На заключительном этапе по полноазимутальному полю рассеянных волн проводят пространственную локализацию областей с повышенным уровнем энергии рассеянных волн и ее оценку на основании имеющейся на площади фактической геолого-геофизической информации (материалы геофизических исследований скважин (ГИС) и данных лабораторного анализа керна), проводят реконструкцию тонкой структуры геологических объектов. Атрибуты дифракции (рассеяния) используют для дифференциации трещиноватых и кавернозных геологических объектов по типу, на основании фактических данных по скважинам.
С точки зрения теории и физики процесса предлагаемый способ заключается в следующем.
Сейсмические волны несут информацию как о строении регулярных протяженных границ раздела (отраженные волны), так и о распределении в среде скоплений мелкомасштабных неоднородностей (дифрагированные-рассеянные волны), что особенно важно для изучения тонкой структуры кавернозно-трещиноватых коллекторов. Построение трехмерных дифракционных изображений неоднородностей в упругих средах обусловлено тем, что рассеяние и отражение волн учитывается корректно именно в трехмерном пространстве, то есть там, где происходят реальные физические процессы. Предложенный способ во многом опирается на ранее выполненные исследования по изучению дифракции и рассеяния сейсмических волн в упругих средах. На которых, в том числе основан способ прототипа - фокусирующие преобразования. Основное отличие заключается в использовании для фокусировки сейсмической энергии специальных асимптотических решений динамической теории упругости - Гауссовых пучков.
Результат достигается за счет инновационного способа построения изображений путем выполнения асимметричного суммирования данных сейсморазведки с использованием специальным образом вычисляемых весовых множителей на основе использования возможности представления полного волнового поля в виде суперпозиции Гауссовых пучков. Благодаря выполняемому асимметричному суммированию со специальным образом построенными весами, удается подавить регулярную часть волнового поля, задаваемую падающими-отраженными-преломленными Гауссовыми пучками и, тем самым, подчеркнуть его рассеянную компоненту. При этом, применяемое подавление регулярной составляющей является наиболее полным как раз из-за того, что Гауссовы пучки определяют наиболее остро направленные волновые поля, чем достигается повышение разрешающей способности. В расчете апертур заключается основное различие между прототипом и предлагаемом способом. Апертура в предложенном способе формируется за счет трассирования Гауссовых пучков к конкретным пунктам возбуждения и приема. Апертура прототипа формируется за счет суммирования пунктов возбуждения и приема на фиксированных базах, которые остаются постоянными для любых точек, в которых строится изображение. Данные различия принципов формирования апертур для предложенного способа, и способа прототипа проиллюстрированы на фиг. 1. Этим, прежде всего, достигается высокая разрешенность способа, что на практике приводит к получению более детальных изображений, с широким спектром, наглядно продемонстрированных на фиг. 2.
Одним из преимуществ изображений в дифрагированных (рассеянных) волнах является возможность подчеркивать на разрезе такие важные с геологической точки зрения объекты, как разломы (дифракция), зоны повышенной трещиноватости (рассеяние) и другие структурные нарушения. Пространственная структура сингулярных объектов находит свое отражение в пространственном распределении порождаемых ими рассеянных/дифрагированных волн. Приведем классическую классификацию сингулярных объектов и соответствующие им типы рассеянных-дифрагированных волн:
• Локальное возмущение параметров среды, все размеры которого меньше длины волны, так называемый точечный дифрагирующий/рассеивающий объект (каверна) порождает дифрагированную/рассеянную волну, распространяющуюся равномерно по всем направлениям и на изображениях в рассеянных волнах, представляется в виде локального разрастания энергии;
• Локальные возмущения, протяженные в пространстве, такие как скопления трещин, геологические разломы. Порождают рассеянные/дифрагированные волны, распространяющиеся в двух направлениях: наиболее интенсивная энергия по нормали к направлению протяженности объекта, а вторичный максимум вдоль сингулярности.
В результате выделения рассеянной компоненты из полного волнового поля, появляются возможности построения локальных объектов акустических неоднородностей в разных комбинациях азимутальности:
• построение изображения полноазимутального поля рассеянных волн (сумма изображений, полученных по всему набору азимутов) для оценки интенсивности рассеяния;
• построение изображений азимутальных сумм рассеянного поля в заданных азимутальных секторах;
• получение сейсмических атрибутов:
- атрибут «структурной дифракции» есть амплитуда полного поля рассеянных волн, умноженная на амплитуду «структурного» рассеяния. Этот атрибут указывает на присутствие дифрагирующих объектов с определенной структурой, прежде всего это ориентированная трещиноватость и разломы.
- атрибут - «азимут структурной дифракции» являющийся азимутальным направлением, перпендикулярным азимутальному направлению, в котором амплитуда рассеяния является максимальной. Этот атрибут указывает на азимутальное направление объектов «структурной дифракции», прежде всего направление трещин.
- атрибут «точечной дифракции» - амплитуда полного поля рассеянных волн, умноженная на амплитуду «точечного» рассеяния. Атрибут указывает на присутствие точечных дифрагирующих объектов без азимутальной ориентации, это могут быть вариации плотности трещин, каверны, изломы разломов.
Заявляемый способ апробирован на различных участках с карбонатным типом коллектора. Для практической реализации способа и в качестве примера была создана модель - сейсмогеологический цифровой двойник одного из характерных месторождений Восточной Сибири характеризующегося сложным геологическим строением. Наибольший интерес на месторождении представляют залежи, приуроченные к сложным по своему геологическому строению рифейскому карбонатному коллектору, с которыми связаны основные объемы запасов нефти и газа. Пустотность коллектора представлена каверново-трещинным типом. При построении сеточной цифровой модели месторождения, каркас был описан множеством границ раздела слоев, полученных по данным сейсморазведки и скважинных наблюдений. Слои были заполнены упругой средой на основе данных ГИС в скважинах. В модели преобладает субгоризонтальное залегание слоев при наличии структурных нарушений - разломов. По данным исследований в скважинах, были учтены неоднородности строения разломов и приразломных зон, осложненными участками высоких деформаций. Помимо того в рефейский резервуар модели включает элементы тонкой структуры, такие трещиноватость, длинной от 5 до 300 м, образующих реалистичные коридоры трещиноватости, и кавернозность, с пустотностью до 20% и толщинами в 1-10 м, подтвержденные результатами ГИС и лабораторного анализа керна, в пределах месторождения. Все описанные особенности модели наглядно проиллюстрированы на фиг. 3а.
Полномасштабное численное моделирование системы многократного перекрытия проводилось для фиксированной системы трехкомпонентных приемников, расположенных на всей апертуре 8 на 10 км с равномерным шагом 25 на 25 м. Источники типа центра расширения (взрыв) располагались на линиях с шагом 50 м, сами линии отстояли друг от друга на 300 м. В качестве сигнала в источнике использовали импульс Рикера с доминирующей частотой 40 Гц. В результате, была сформирована реальная полевая система, обеспечивающая полноазимутальную запись. По сейсмогеологическому двойнику месторождения были рассчитаны полные синтетические сейсмические данные, достоверно отображающие как регулярные, так и рассеянные волновые поля, вызванные наличием мелкомасштабных неоднородностей, полностью идентичные реальным данным. В результате расчетов получен синтетический куб, который представлен на фиг. 3б.
Для практической реализации способа было создано программное обеспечение для ЭВМ «Построение трехмерных изображений рассеянных волн во временной области на основе Гауссовых пучков» (свидетельство о государственной регистрации №2019613485). Оно использует следующие входные данные: сейсмограммы и скоростная модель среды в сейсмическом формате SEG-Y. В результате расчета программы строится трехмерное изображение геологической среды в рассеянных волнах в сейсмическом формате SEG-Y.
Этап 1. По суммарным временным разрезам, полученным после миграции по синтетическим данным, проведен анализ поведения отражающих горизонтов на временных разрезах для выбора параметров асимметричного суммирования - углы и азимуты для трассирования Гауссовых пучков в направлении источников и приемников, обеспечивающих отсутствие регулярных границ раздела в области частичного восстановления. Поскольку залегание границ - субгоризонтальное, был выбран угол подхода, достаточный для разрушения плоских волн, в интервале от 20 до 30°, с шагом 5°, для статистического накопления.
Этап 2. Используя полученные параметры, с помощью программного обеспечения «Построение трехмерных изображений рассеянных волн во временной области на основе Гауссовых пучков», используя в качестве входных данных сейсмограммы и скоростную модель синтетических данных, были вычислены селективные изображения рассеянных волн в области углов и азимутов трассирования.
Этап 3. Селективные изображения были просуммированы по всем углам и азимутам трассирования, в результате чего получено общее полноазимутальное поле рассеянных волн.
Этап 4. Из общего полноазимутального поля рассеянных волн с помощью программного обеспечения «Построение трехмерных изображений рассеянных волн во временной области на основе Гауссовых пучков» рассчитаны атрибуты дифракции/рассеяния: атрибут «структурной дифракции», «азимут структурной дифракции» и «точечной дифракции».
Этап 5. На заключительном этапе, по полноазимутальному полю энергии рассеянных волн проведена пространственная локализация областей с повышенным уровнем энергии и на основании имеющейся на месторождении априорной скважинной информации, была проведена реконструкция структуры целевого рифейского резервуара. Сопоставив горизонтальные срезы цифровой модели резервуара и ее изображение в рассеянных волнах, полученного в результате работы способа, отмечается схожесть между моделью и изображением - существует почти однозначное соответствие: местоположение всех дифрагирующих объектов, включая разломы, трещиноватость, кавернозность и их относительная амплитуда реконструированы корректно, что наглядно продемонстрировано на фиг. 4.
Затем на основе атрибутов дифракции/рассеяния была произведена дифференциация геологических объектов целевого резервуара по типу дифракции. Зоны развития направленной трещиноватости на площади, подтвержденные микросканерами ГИС в скважинах №1 и №2, уверенно локализуются на изображении рассеянных волн атрибута «структурной дифракции», а азимут простирания трещиноватости, соответствует азимуту атрибута «азимут структурной дифракции». Имеющаяся на площади зона интенсивной кавернозности, фиксируется на изображении атрибута «точечной дифракции», что соответствует данным ГИС и исследования керна в скважине №3, где выявлено интенсивное выщелачивание карбонатов. Вышеописанные результаты проиллюстрированы на фиг. 5.
Как видно из приведенного примера, заявленный способ позволяет четко локализовать мелкомасштабные неоднородности карбонатного рифейского коллектора и дифференцировать их по типу дифракции, а именно, на зоны трещиноватости и кавернозности.
Claims (1)
- Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные, характеризующийся тем, что по стандартным сейсмическим изображениям среды, полученным с использованием отраженных волн, проводят анализ поведения отражающих горизонтов на временных разрезах, определяют углы наклона границ целевых горизонтов, по которым рассчитывают оптимальные параметры асимметричного суммирования сейсмических данных - угол и азимут трассирования Гауссовых пучков, которые определяют весовые множители и апертуру, обеспечивающие максимально надежное построение поля энергии рассеянных волн, с полученными параметрами проводят асимметричное суммирование сейсмограмм после сигнальной обработки и куба средних скоростей, получают селективные изображения в области углов наклона и азимутов трассирования, которые суммируют по всем полученным углам и азимутам трассирования Гауссовых пучков, из полученного общего полноазимутального изображения поля энергии рассеянных волн рассчитывают сейсмические атрибуты дифракции/рассеяния, основанные на анализе селективных изображений в зависимости от азимута и типа дифракции/рассеяния, затем по полноазимутальному полю энергии рассеянных волн проводят пространственную локализацию областей с повышенным уровнем энергии рассеянных волн и на основании ее оценки и имеющейся на площади фактической геолого-геофизической информации осуществляют реконструкцию тонкой структуры геологических объектов, и, используя атрибуты дифракции/рассеяния и фактические данные по скважинам, осуществляют дифференциацию трещиноватых и кавернозных геологических объектов по типу.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2020131555A RU2758416C1 (ru) | 2020-09-25 | 2020-09-25 | Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2020131555A RU2758416C1 (ru) | 2020-09-25 | 2020-09-25 | Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2758416C1 true RU2758416C1 (ru) | 2021-10-28 |
Family
ID=78466625
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2020131555A RU2758416C1 (ru) | 2020-09-25 | 2020-09-25 | Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2758416C1 (ru) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2789759C1 (ru) * | 2022-06-20 | 2023-02-09 | Общество с ограниченной ответственностью "Сахалинская Энергия (ООО "Сахалинская энергия") | Способ определения границ трещиноватой зоны |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2490677C2 (ru) * | 2011-11-28 | 2013-08-20 | Александр Алексеевич Архипов | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления |
| RU2563323C1 (ru) * | 2014-04-02 | 2015-09-20 | Открытое акционерное общество "Нефтяная компания "Роснефть" | Способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения |
| CN103984012B (zh) * | 2014-04-16 | 2016-06-29 | 刘豫宝 | 基于叠前高斯束深度偏移的绕射波场分离方法 |
| CN109100783A (zh) * | 2017-06-20 | 2018-12-28 | 中国石油化工股份有限公司 | 一种方位反射角度域高斯束层析反演方法及系统 |
-
2020
- 2020-09-25 RU RU2020131555A patent/RU2758416C1/ru active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2490677C2 (ru) * | 2011-11-28 | 2013-08-20 | Александр Алексеевич Архипов | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления |
| RU2563323C1 (ru) * | 2014-04-02 | 2015-09-20 | Открытое акционерное общество "Нефтяная компания "Роснефть" | Способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения |
| CN103984012B (zh) * | 2014-04-16 | 2016-06-29 | 刘豫宝 | 基于叠前高斯束深度偏移的绕射波场分离方法 |
| CN109100783A (zh) * | 2017-06-20 | 2018-12-28 | 中国石油化工股份有限公司 | 一种方位反射角度域高斯束层析反演方法及系统 |
Non-Patent Citations (2)
| Title |
|---|
| Кутовенко М.П., Протасов М.И., Чеверда В.А. " Использование Гауссовых пучков для построения сейсмических изображений в истинных амплитудах по многокомпонентным данным", Технологии сейсморазведки, 2010, номер 4, с. 3 - 13. * |
| ПЕТРОВ Д.А., МЕЛЬНИК А.А., ШИЛИКОВ В.В., ТУЗОВСКИЙ А.А. и др. "Выявление трещиновато-кавернозных коллекторов на основе интерпретации сейсмических рассеянных волн методом гауссовых пучков", Нефтяное хозяйство, 2019, номер 1, с.6-10. * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2789759C1 (ru) * | 2022-06-20 | 2023-02-09 | Общество с ограниченной ответственностью "Сахалинская Энергия (ООО "Сахалинская энергия") | Способ определения границ трещиноватой зоны |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wuestefeld et al. | A strategy for automated analysis of passive microseismic data to image seismic anisotropy and fracture characteristics | |
| RU2563323C1 (ru) | Способ реконструкции тонкой структуры геологического объекта и прогноза его флюидонасыщения | |
| EP3140676B1 (en) | System and method for analyzing geologic features using seismic data | |
| Lee et al. | Delineation of gas hydrate reservoirs in the Ulleung Basin using unsupervised multi-attribute clustering without well log data | |
| Jones et al. | Improved microseismic event location by inclusion of a priori dip particle motion: a case study from Ekofisk | |
| RU2451951C2 (ru) | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам | |
| WO2016032353A1 (ru) | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам | |
| Takougang et al. | Extraction and characterization of faults and fractures from 3D VSP data in a carbonate reservoir: A workflow | |
| Barthwal et al. | Passive seismic tomography using recorded microseismicity: Application to mining-induced seismicity | |
| Sanda et al. | The integrated approach to seismic attributes of lithological characterization of reservoirs: case of the F3 Block, North Sea-Dutch Sector | |
| Protasov et al. | 3D diffraction imaging attributes and their application for fault and fracture localization and characterization | |
| EP2917857B1 (en) | Fracture characterization from refraction travel time data | |
| RU2732035C1 (ru) | Способ определения трещинной пористости пород | |
| Jones et al. | Characterization of fractures and faults: a multi‐component passive microseismic study from the Ekofisk reservoir | |
| RU2758416C1 (ru) | Способ реконструкции тонкой структуры геологических объектов и их дифференциации на трещиноватые и кавернозные | |
| El-Emam et al. | Interbed multiple prediction and attenuation: Case history from Kuwait | |
| CN113075656A (zh) | 一种探地雷达和三维地震综合预测煤矿工作面断裂方法 | |
| Engelsfeld et al. | Investigation of underground cavities in a two‐layer model using the refraction seismic method | |
| CN104345337B (zh) | 一种用于地震反演的时控储层参数建模方法 | |
| RU2300126C1 (ru) | Способ геофизической разведки для выявления малоамплитудных тектонических нарушений нефтегазопродуктивных горных пород в трехмерном межскважинном пространстве | |
| Nesterkina et al. | Assessment of the seismic impact of industrial explosions in the central part of the East European Platform | |
| Xiang et al. | Improving the resolution of impedance inversion in karst systems by incorporating diffraction information: A case study of Tarim Basin, China | |
| Telegin | Possibilities of seismic exploration for crystalline basement study | |
| Carcione et al. | Seismic modeling study of the Earth's deep crust | |
| Cheverda et al. | Three-Dimensional Model of Oil and Gas Reservoirs Based on Gaussian Beam Processing of Scattered Seismic Waves |