Определение направления трещины ГРП с дневной поверхности с использованием тензора сейсмического момента

Биряльцев Евгений Васильевич; Рыжов В.А.; Феофилов Сергей Анатольевич; Шарапов И.Р.; Камилов Марсель Робертович; Рыжов Дмитрий Александрович; Мокшин Евгений Владимирович

УДК 622.245

ОПРЕДЕЛЕНИЕ НАПРАВЛЕНИЯ ТРЕЩИНЫ ГРП С ДНЕВНОЙ ПОВЕРХНОСТИ С ИСПОЛЬЗОВАНИЕМ ТЕНЗОРА СЕЙСМИЧЕСКОГО МОМЕНТА

Е.В. Биряльцев1, В.А. Рыжов1*, С.А. Феофилов1, И.Р. Шарапов1, М.Р. Камилов2,

Д.А. Рыжов1, Е.В. Мокшин2

1ЗАО «Градиент», Казань, Россия 2ООО «Градиент технолоджи», Казань, Россия

Мониторинг микросейсмических событий с дневной поверхности применяется при разработке трудноизвле-каемых запасов, особенно в процессе гидроразрыва пласта (ГРП). Приводится сравнение методов микросейсмического мониторинга ГРП с дневной поверхности: дифракционного суммирования, метода моделирования в обратном времени, спектральный метод. В (Aki & Richards 1980) показано, что достоверность локализации сейсмических событий на фоне коррелированных шумов существенно увеличивается при применении метода максимального правдоподобия. Метод максимума правдоподобия позволяет исключить коррелированную компоненту шума, а также позволяет определить не только координаты, но и тензор сейсмического момента при локализации с дневной поверхности.

Оценка тензора сейсмического момента сейсмического события позволяет определять тип события: «Центр взрыва» (EXP), «Трещина разрыва» (TC), «Сдвиг» (DC) и «Компенсированный линейный диполь» (CLVD), а также ориентацию трещины каждого события, даже когда нет явной асимметрии пространственного распределения облака событий.

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

Ключевые слова: микросейсмический мониторинг, тензор сейсмического момента, трещиноватость, сейсмическое событие, метод максимального правдоподобия

DOI: https://doi.Org/10.18599/grs.19.3.13

Для цитирования: Биряльцев Е.В., Рыжов В.А., Феофилов С.А., Шарапов И.Р., Камилов М.Р., Рыжов Д.А., Мокшин Е.В. Определение направления трещины ГРП с дневной поверхности с использованием тензора сейсмического момента. Георесурсы. 2017. Т. 19. № 3. Ч. 1. С. 229-233. DOI: https://doi.Org/10.18599/grs.19.3.13

Введение

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

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

*Ответственный автор: Василий Александрович Рыжов E-mail: v. [email protected]

Технология локации

Направление распространения трещины обычно определяют по ориентации облака событий, сопровождающих процесс образования трещины ГРП. Для локализации микросейсмических событий применяются различные техники наблюдения и обработки микросейсмической информации. Наиболее хорошо известна техника дифракционного суммирования, которая применяется для локализации микросейсмических событий более 50 лет (Кгеу, 1952; Hagedoom, 1954). Главный принцип метода дифракционного суммирования заключатся в расчете временных задержек, соответствующих времени хода сигнала от анализируемых точек геосреды до точек приема. После применения вычисленных задержек амплитуды записей суммируются с учетом рассчитанных задержек.

Известен также метод расчетов в обратном времени (Gajewski & Tessmer, 2005), данный метод базируется на численном моделировании процесса распространения упругих волн в интересуемой части среды. Сигналы, принимаемые во время мониторинга, инверсируются по шкале времени и используются в качестве источников возбуждения в месте расположения соответствующих датчиков. За время, равное времени хода от источника до наблюдаемой группы приемников, импульсы локализуются в месте возникновения.

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ ■НН^Н

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

В (Aki & Richards, 1980) показано, что достоверность локализации сейсмических событий на фоне коррелированных шумов существенно увеличивается при применении метода максимального правдоподобия. Авторы используют его для локации одного очага землетрясения на фоне другого. Однако, наличие коррелированных шумов не характерно для сейсмологии, учитывая значительный разнос сейсмологических станций. При локации микросейсмических событий, напротив, коррелированные шумы от работающего оборудования (флот ГРП, нефтегазовая инфраструктура) составляют основную часть шумов. Как показано в (Birialtsev et al., 2017), метод максимума правдоподобия позволяет исключить коррелированную компоненту шума, а также позволяет определить не только координаты, но и тензор сейсмического момента при локализации с дневной поверхности.

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

Для применения метода максимального правдоподобия необходимо знать форму полезного сигнала. В общем случае полезным сигналом является полноволновой отклик среды на импульсное воздействие. Рассчитать форму полезного сигнала в геологической среде можно с помощью применения полноволнового 3D численного моделирования (Биряльцев и др., 2008). Для расчета тензора сейсмического момента необходимо

моделирование 6 типов импульсных воздействий различного типа (Birialtsev et а1., 2017).

Полноволновое 3D численное моделирование и локация событий методом максимального правдоподобия требуют значительных вычислительных затрат, поэтому для получения результатов в приемлемое время применяются суперкомпьютерные вычисления (Галимов и др. 2010).

Технология полноволновой локации, характеризуется рядом особенностей:

1. Регистрация при мониторинге ГРП выполняется с дневной поверхности независимыми комплектами широкополосных высокочувствительных сейсмометров, устанавливаемых в самых тихих местах территории (Рыжов и др., 2015);

2. Благодаря применению полноволнового 3D численного моделирования, при локации используется полная информация о сигнале в местах установки датчиков по трем компонентам (полноволновой отклик, включая волны сжатия, сдвига, обменные, переотраженные) от единичных импульсных воздействий;

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

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

В то же время, технологический скачок в развитии

Рис. 1. Пример результата мониторинга ГРП с определением тензора сейсмического момента событий

23

Ш GEDRESURSY

суперкомпьютерных вычислительных систем (Галимов и др., 2010; Демидов и др., 2010; Биряльцев и др., 2015; Беляева и др., 2017), позволил применять ресурсоемкие методы локации (Birialtsev et al., 2017), использующие максимально полную информацию о сейсмическом событии.

Полноволновая локация применяется для решения задач в нефтегазовой промышленности с 2011 года, вышло несколько статей с результатами её применения (Biryaltsev et al.,, 2016; Рыжов и др., 2015; Хисамов и др., 2015; Шабалин и др. 2013).

Результаты локации

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

Следовательно, выявленные в процессе мониторинга ГРП зоны повышенной микросейсмической активности, могут быть ассоциированы с процессами, протекающими в пласте под влиянием ГРП и неразрывно включающими в себя следующее:

1. Образование/смыкание трещин;

2. Раскрытие трещин в процессе заполнения их проппантом;

3. Деформаций породы в областях с предкритическим состоянием за счет продвижения фронта давления по естественным каналам фильтрации;

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

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

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

Результатом локации является набор событий с координатами в пространстве и тензором сейсмического момента, на основе которого оценивается степень принадлежности события к базовым типам событий. Выделяется несколько базовых типов событий: «Центр взрыва» (EXP), «Трещина разрыва» (TC), «Сдвиг» (DC) и «Компенсированный линейный диполь» (CLVD). Ориентация трещины разрыва оценивается только для событий с высоким весом «Трещина разрыва» (TC). Для событий типа «Центр взрыва» (EXP) говорить об азимуте нецелесообразно, так как все его три собственных вектора равнозначны и параметр азимута в данном случае определяется неустойчиво. Для событий смешанного типа, например, 45% TC и 40% EXP, оценка азимута трещины

справедлива, однако с меньшей достоверностью, чем для более выраженной трещины разрыва 80% TC. Для событий типа DC и CLVD определение ориентации трещины не выполняется, напротив, такие события исключаются из обработки, как события, не связанные с изменением объема (раскрытием/смыканием трещины ГРП).

На рисунке 1 представлены результаты локации событий с поверхности технологией полноволновой локации для вертикальных глубин порядка 2.7 км. В результате имитационного моделирования получены оценки точности локации для данных условий: погрешность локации не более 35м, погрешность определения азимута трещины не более 15 градусов. На результирующей карте события в области порта имеют достаточно большой разброс, не позволяющий уверенно выделить направление трещины только по их местоположению, однако благодаря розе азимутов каждого события, через тензор сейсмического момента стало возможным оценить направление трещины.

Заключение

Таким образом, технология локализации микросейсм от ГРП с использованием метода максимального правдоподобия позволяет определять направление трещин ГРП даже в условиях неуверенной пространственной локализации микросейсмических событий. Неуверенная локализация событий может быть вызвана трудно устранимыми причинами: высоким уровнем шума техногенной активности на поверхности, низкочастотным рабочим диапазоном (из-за затухания высокочастотной компоненты сигнала вследствие большой глубинности), а также самим сложным характером распространения трещины, как fracture fabric (Cipolla, 2011), образующую сеть параллельных трещин. Определение тензора сейсмического момента и направления трещины с его использованием более устойчиво к перечисленным факторам, что оправдывает применение вычислительно более сложного метода максимума правдоподобия в условиях сложной геологической обстановки.

Литература

Беляева А.А., Биряльцев Е.В., Галимов М.Р., Демидов Д.Е., Елизаров А.М., Жибрик О.Н. Кластерная архитектура программно-технических средств организации высокопроизводительных систем для нефтегазовой промышленности. Программные системы: теория и приложения. 2017. Т. 8. № 1. С. 151-171. http://psta.psiras.ru/read/psta2017_7_151-171.pdf

Биряльцев Е.В., Богданов П.Б., Галимов М.Р., Демидов Д.Е., Елизаров А.М. Программно-техническая платформа высокопроизводительных вычислений для нефтегазовой промышленности. Программные системы: теория и приложения. 2015. Т. 7. № 1. С. 15-27. http://psta. psiras.ru/read/psta2016_1_15-27.pdf

Биряльцев Е.В., Бережной Д.В., Биряльцева Т.Е., Храмченков М.Г. Результаты численного моделирования распространения природных микросейсм в зонах залегания нефтегазовых залежей. Конференция EAGE, ГЕОМОДЕЛЬ-2008. Геленджик. 2008.

Галимов М.Р., Биряльцев Е.В. Некоторые технологические аспекты применения высокопроизводительных вычислений на графических процессорах в прикладных программных системах. Вычислительные методы и программирование. 2010. Т. 11. С. 77-93

Исламов Д.Э. Совершенствование методов проектирования операций по гидроразрыву пластов для повышения продуктивности скважин. Дисс. канд. тех. наук. Тюмень. 2015. 124 с.

Рыжов В.А., Шарапов И.Р., Биряльцев Е.В., Феофилов С.А., Рыжов Д.А., Камилов М.Р., Степанов А.И. Опыт применения метода полноволновой локации при микросейсмическом мониторинге МГРП горизонтальной скважины в Западной Сибири. EAGE—Horizontal Wells. Problems and Prospects. 2015.

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ HHHB :

ШШШШягЗЬ'

Е.В. Биряльцев, В.А. Рыжов, С.А. Феофилов и др.

E.V. Birialtsev, V.A. Ryzhov, S.A. Feofilov et al.

Сметанников О.Ю., Кашников Ю.А., Ашихмин С.Г., Шустов Д.В. Численная модель развития трещины при повторном гидроразрыве плата. Вычислительная механика сплошных сред. 2015. Т. 8. № 2. С. 208-218 Хисамов Р.С., Ахметшина А.С., Таипова В.А., Салихов М.М., Шарапов И.Р. Пассивный наземный микросейсмический мониторинг гидроразрыва пласта в нагнетательных скважинах ОАО «Татнефть»: результаты и их достоверность. Нефтяное хозяйство. 2015. № 7. С. 34-38 Шабалин Н.Я., Биряльцев Е.В., Рыжов В.А., Мокшин Е.В., Феофилов С.А., Шарапов И.Р., Рыжов Д.А. Мониторинг многостадийного ГРП с дневной поверхности. Теоретические подходы и практические результаты. Экспозиция НЕФТЬ ГАЗ. 2013. № 6. С. 40-43.

Aki, K. and Richards, P.G. Quantitative seismology. Freeman and Co. 1980. Birialtsev E.V., Demidov D.E., Mokshin E.V. Determination of moment tensor and location of microseismic events under conditions of highly correlated noise based on the maximum likelihood method. Geophysical Prospecting. 2017. doi: 10.1111/1365-2478.12485

Biryaltsev E.V., Shabalin N.Y., Ryzhov V.A., Sharapov I.R. Application of Full Wave Location technology with determining seismic moment tensor of events for Hydraulic Fracture Monitoring and Natural Fractures. Sixth EAGE Workshop on Passive Seismic: From Data to Decisions. 2016.

Cipolla C., Weng X., Mack M., Ganguly U., Gu H., Kresse O., Cohen C. Integrating microseismic mapping and complex fracture modeling to characterize fracture complexity. Society of Petroleum Engineers. 2011. doi:10.2118/140185-MS.

Demidov D., Ahnert K., Rupp K. and Gottschling P. Programming CUDA and OpenCL: A case study using modern C++ libraries. SIAM Journal on Scientific Computing. 2013. 35(5). Pp. 453-472.

Hagedoorn J.G. A process of seismic reflection interpretation. Geophysical Prospecting. 1954. 2. Pp. 85-127.

Gajewski, D. and Tessmer, E. Reverse modelling for seismic event characterization. Geophys. J. Int. 2005. 163(1). Pp. 276-284.

Kushnir A., Varypaev A., Dricker I., Rozhkov M. and Rozhkov N. Passive surface microseismic monitoring as a statistical problem: location of weak microseismic signals in the presence of strongly correlated noise. Geophys. J. Int. 2014. 198 (2). Pp. 1186-1198.

Krey T. The significance of diffraction in the investigation of faults. Geophysics. 1952. 17. Pp. 843-858.

Shapiro S.A. Fluid-Induced Seismicity. Cambridge University Press. 2015. 289 p.

Сведения об авторах

Евгений Васильевич Биряльцев - канд. тех. наук, зам. генерального директора по науке и новым технологиям, ЗАО «Градиент»