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

Тип работы:
Реферат
Предмет:
Геофизика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 550. 348. 432
В. В. Сказка, С. В. Сердюков, А. С. Сердюков
МОДЕЛИРОВАНИЕ МИКРОСЕЙСМИЧЕСКОГО ШУМА ПОРОДНЫХ МАССИВОВ БЛОЧНОГО СТРОЕНИЯ
Предложена модель микросейсмического шума в горных породах блочного строения. Приведены результаты численных исследований формы сигналов и спектрального состава микросейсмического шума продуктивных пластов блочного строения, зависимость радиуса корреляции среднечастотных микросейсмов от положения приемников по отношению к пласту, генерирующему шум. Даны рекомендации по алгоритмам выделения полезного сигнала и локации эндогенных сейсмических событий на основе решения обратных кинематических задач.
The model of microseismic noise in block structuring rocks is offered. Results of numerical researches of the signal form and spectral structure of microseismic noise in productive block structuring layers, dependence of the correlation radius of mid-frequency noise on the receivers position in relation to a layer generating noise are resulted. Recommendations about algorithms of the useful signal allocation and a location of internal seismic events on the basis of the inverse kinematic problems decision are made.
Ключевые слова: массив горных пород, блочная среда, микросейсмиче-ский шум.
Key words: rock massif, block media, microseismic noise.
Введение
Для породных массивов характерны дислокации различных масштабов. Их система задает блочное строение геологической среды [1], особенности ее деформирования и выделения упругой энергии в виде эндогенного микросейсмического излучения [2- 3].
Сейсмические методы исследования продуктивных пластов предусматривают выделение полезных сигналов на фоне помех, связанных с релаксацией блочной среды к равновесному состоянию, нарушенному процессом разработки [4]. В ряде случаев, например для карбонатных нефтегазовых коллекторов, выраженная система трещин ведет к определенным закономерностям в пространственном распределении источников шума по частотам и энергии [5]. Для решения модельных задач в таких средах необходимо учитывать влияние блочной структуры на характеристики шума. Простые подходы добавления к сигналу случайных помех не применимы.
В статье рассмотрена модель микросейсмического шума, учитывающая блочное строение массива горных пород.
131
© Сказка В. В., Сердюков С. В., Сердюков А. С., 2013
Вестник Балтийского федерального университета им. И. Канта. 2013. Вып. 4. С. 131 — 139.
132
1. Модель микросейсмического события
Пусть микросейсмическое событие генерируется разрывом смещений вдоль поверхности площадью А, нормаль к которой задана единичным вектором V. Скачок смещений т параллелен плоскости разрыва. Площадка находится в точке с координатами % = (^1,^2,^3). Прием сейсмического сигнала производят в точке х = (х1,х2,х3). Введем следующие обозначения: г = |х -?| - расстояние от источника до приемника- у = (у1,у2,у3), у = (х{ - ^) /г- Ур, У5 — скорости распространения соответственно продольных и поперечных волн- и — коэффициент Ламе- р — плотность среды.
Компоненты волнового поля и = (м1,и2,и3) в однородной среде от микросейсмического события будем вычислять по формуле [6]
и =
(30у у у V — 6v у — 66 у V
/ п р ц ц П'- р пр ц ц
4лрг4
(12 у у у V — 2у у — 26 у V Л
І пі р ц ц п ц пр ц ц
| С г -У
иАт/?р (і -т'-)Ст+
4жрУр г2
(12улуу — 3кгг — 36пРгуа Л
4жрУ^ г2
иАй,(() —
4ярУр г р Т7
иАтР V — -) -У
(1)
Ур
(2уууv-vу-6 у V Л
/ п& gt- р ц ц п& gt- р пр ц а
4жрУ^г
иАтр (^ - у),
где t — время- 8^ - символ Кронекера.
Сигнала микросейсмического события зададим импульсом Риккера:
ДО = (1 — 2р?)ехр (-р1?), (2)
где р — затухание- р / п = о0 — доминантная частота упругих колебаний. Поглощения сигнала средой моделируем функцией:
1 & quot-г '-е-Нате°*
Р (г, ?) = ^= [ /(аУИагвшс1а, & lt-2я -і
(3)
где г — расстояние от источника до приемника- а — коэффициент поглощения- /(о) = 1_ I& quot- /(?)е-ойЬ — преобразование Фурье исходно-
у11ж о
го сигнала (1) — со — круговая частота колебаний, Гц.
На рисунке 1 показаны графики функций (1) и (3), а также их амплитудные спектры при: г = 1000 м, р = 40^ Гц,
а = 0,15/ (1400) с/м.
0,006
0,004
о 0,002
0,000
20
40
«I---------1---------Т---------1---------Т---------1---------Т---------1---------Т---------1 со, Гц
80
100
Рис. 1. Форма импульсов (а) и их амплитудные спектры (б):
— Т (г, І) — г = 1000 м
Интеграл (3) выразим в явном виде
F (r, t) =
(-1 / 8 -1 / 4 p2 r2a2 +1 / 4 p212 -1 / 2 lp2 rat)
M-t2 +2 irat+ r 2a2
erf (pra + ipt) +
+(-1 /4p2r2a2 +1 /2ip2rat +1 /4p2t2 -1 /8)e (-2 p^t2 + 2 p2жr2a2 +ж+ 4 їp2жrat
p (-t2 -2 irat+r a1
'-'-erf (ra-ipt) +
'-l-t2 +2 irat +r 2a2
(4)
+1 / 8
+1 / 8
(-.2 2−2 22 -.2 p2 f-t2 — 2 irat+r 2a2
ж-2p2жt2 + 2p^ra2 -4ip^rat)e (-
-1 / 2
pra
у[ж
1 Г 2
где erf (z) =. I exp (-z)dz
комплексный интеграл вероятности.
о
Будем полагать, что функция Ф (0 в (1) имеет следующий вид:
^) = Iи -?(г, ^), (5)
где Т (г, Ї) определена в (3) — векторы 1и и V взаимно ортогональны.
2. Модель шума и результаты численных расчетов
Пусть все микросейсмические события приурочены к продуктивному пласту мощностью 20 м, имеющему блочное строение и расположенному горизонтально в однородной среде со следующими значениями сейсмических волн: Ур = 2400 м/с, У5 = 2400 /л/3 м/с.
а
б
ж
133
134
Будем полагать, что объемы дислокаций разного масштаба равны между собой [3], их толщины пропорциональны линейному размеру [7]. Считаем скорость подвода энергии к дислокациям разного масштаба одинаковой и постоянной во времени.
Тогда энергия и период следования сигналов отдельного микро-сейсмического источника будут пропорциональны площади дислокации, а пространственная плотность событий — обратно пропорциональна линейному размеру блоков в третьей степени, что схоже с известным сейсмологическим законом повторяемости землетрясений [8].
В соответствии с экспериментальными данными [5] зададим доминантные частоты блоков в виде дискретного набора:
а0к = 20 х Ур/К10- К = 1,2,…, 12,
где О = 0,022 [7], 10 — минимальный размер блока (дислокации), отвечающего за верхнюю частоту генерируемых колебаний. В численных экспериментах был задан диапазон сейсморазведочных частот от 10 Гц (К = 12) до 120 Гц (К = 1), представляющий практический интерес.
Будем полагать, что источники шума равномерно распределены в объеме продуктивного пласта. Для каждой группы источников с одинаковым значением доминантной частоты зададим свою пространственно-временную плотность событий 1/60 х (®0 /(00г), (м3Хс), учитывающую указанную выше зависимость от линейного размера блоков. Энергию и время следования событий, генерируемых дислокациями различного размера, устанавливаем пропорционально значению параметра (®01 /®0к), отражающего отношение их площадей.
Направление пары ортогональных векторов п и /и, фигурирующих в (1), (5), выбираем случайным образом. Временные последовательности микросейсмических сигналов генерируем независимо для каждой группы источников с равномерным распределением по времени в период наблюдений длительностью 60 с.
Рассматривались два случая: 1) приемники расположены внутри продуктивного пласта блочного строения, 2) регистрирующая система находится в однородной среде, вмещающей пласт (рис. 2).
Ориентация трех компонентных приемников согласована с системой координат.
В первом эксперименте приемники расположены внутри продуктивного пласта на расстоянии в один метр между собой.
Микросейсмические шумы, регистрируемые приемниками, рассчитывались по формуле (1) по каждой группе источников отдельно. Дискретизация по времени составляла 5×10−4 с. Полученные по группам результаты были просуммированы. Пример итоговой синтетической сейсмограммы Х-компоненты микросейсмического шума с разной степенью детализации временной шкалы показан на рисунке 3.
Отметим, что пики большой амплитуды на сейсмограмме (рис. 3) связаны с сигналами от близлежащих источников.
5
4
3
2
1
5 м
15 м
-> X дни йшт 12 3 4 Приемники в пласте 5
& gt- & lt- г Пласт
135
Рис. 2. Схема размещения приемников в продуктивном пласте
її ііпі 1
. 1 '- 1 '- 1 '- 1» '- 1 р|'-«1 Т & quot-Тії і '- і '- ч- '- 1 —
І С
28 29 30 31 32 33 34 35 36 37 38
І С
С
Рис. 3. Синтетическая сейсмограмма Х-компоненты микросейсмического шума внутри продуктивного пласта (эксперимент 1)
Для определения взаимосвязи сигналов, принятых разными приемниками, были вычислены коэффициенты корреляции между соответствующими компонентами сигналов в разных точках приема. Полученные значения представлены в таблице 1.
У
X
Таблица 1
Корреляция микросейсмов в продуктивном пласте
Х-компонента
136
№ приемника 1 2 3 4 5
1 1 0,5 025 281 0,3 302 378 0,2 320 549 0,1 723 080
2 0,5 025 281 1 0,6 789 085 0,4 667 565 0,3 342 566
3 0,3 302 378 0,6 789 085 1 0,7 373 485 0,5 239 655
4 0,2 320 549 0,4 667 565 0,7 373 485 1 0,7 986 688
5 0,1 723 080 0,3 342 566 0,5 239 655 0,7 986 688 1
У-компонента
№ приемника 1 2 3 4 5
1 1 0,4 243 855 0,2 492 039 0,1 677 603 0,1 041 544
2 0,4 243 855 1 0,4 993 286 0,3 127 862 0,1 848 669
3 0,2 492 039 0,4 993 286 1 0,6 221 576 0,3 526 777
4 0,1 677 603 0,3 127 862 0,6 221 576 1 0,6 698 635
5 0,1 041 544 0,1 848 669 0,3 526 777 0,6 698 635 1
Х-компонента
№ приемника 1 2 3 4 5
1 1 0,4 953 246 0,2 873 339 0,1 910 331 0,1 251 473
2 0,4 953 246 1 0,5 179 037 0,3 260 742 0,2 162 353
3 0,2 873 339 0,5 179 037 1 0,5 996 217 0,3 523 680
4 0,1 910 331 0,3 260 742 0,5 996 217 1 0,6 488 788
5 0,1 251 473 0,2 162 353 0,3 523 680 0,6 488 788 1
Видно, что радиус корреляции микросейсмов в продуктивном пласте меньше 1 м.
Во втором эксперименте профиль наблюдений расположен перпендикулярно пласту (рис. 2). Расстояние между приемниками составляет 5 м, минимальное расстояние от приемников до кровли пласта — 15 м. Из-за большего расстояния между приемниками и источниками микросейсмических событий время наблюдения увеличено до 600 с.
Пример синтетической сейсмограммы Х-компоненты микросейс-мического шума, зарегистрированного приемником № 2, показан на рисунке 4 с разной степенью детализации временной шкалы.
Полученные графики имеют более сглаженный характер по сравнению с сейсмограммой на рисунке 3, не содержат выраженных пиков сигнала, что особенно заметно на коротких временных интервалах.
На рисунке 5 представлены графики спектральной плотности микросейсмического шума в продуктивном пласте и вне его пределов. Внутри пласта спектр обогащен высокочастотными составляющими, менее сглажен за счет большего влияния близкорасположенных источников.
228
и С
137
Рис. 4. Синтетическая сейсмограмма Х-компоненты микросейсмического шума вне продуктивного пласта (эксперимент 2, приемник № 2)
5
4
3
2
1
ю, Гц
а
о, Гц
б
Рис. 5. Спектральная плотность Х-компоненты микросейсмического шума в продуктивном пласте (а) и вне продуктивного пласта (б) (приемник № 2)
0
Значения коэффициентов корреляции между соответствующими компонентами сигналов в разных точках приема регистрирующей системы вне пласта (рис. 2) приведены в таблице 2.
138
Отметим высокие значения коэффициентов корреляции 2-компоненты микросейсмического сигнала между соседними приемниками при выносе сейсмоприемников вне зоны генерации шумов, что может быть использовано для локации таких зон в геологической среде.
Предложенная модель микросейсмического шума хорошо описывает форму сигналов и спектральный состав микросейсмического шума продуктивных пластов блочного строения, зависимость радиуса корреляции среднечастотных микросейсмов от положения приемников по отношению к пласту, генерирующему шум.
Таблица 2
Корреляция микросейсмов вне продуктивного пласта
Х-компонента
№ приемника 1 2 3 4 5
1 1 0,6 443 400 0,1 476 195 -0,1 706 539 -0,2 785 052
2 0,6 443 400 1 0,6 494 048 0,1 440 907 -0,1 885 829
3 0,1 476 195 0,6 494 048 1 0,6 543 913 0,1 423 415
4 -0,1 706 539 0,1 440 907 0,6 543 913 1 0,6 586 707
5 -0,2 785 052 -0,1 885 829 0,1 423 415 0,6 586 707 1
У-компонента
№ приемника 1 2 3 4 5
1 1 0,6 477 117 0,1 534 663 -0,1 678 226 -0,2 803 487
2 0,6 477 117 1 0,6 518 418 0,1 475 161 -0,1 875 984
3 0,1 534 663 0,6 518 418 1 0,6 556 661 0,1 433 333
4 -0,1 678 226 0,1 475 161 0,6 556 661 1 0,6 587 450
5 -0,2 803 487 -0,1 875 984 0,1 433 333 0,6 587 450 1
Х-компонента
№ приемника 1 2 3 4 5
1 1 0,8 524 379 0,5 550 823 0,2 615 874 0,366 206
2 0,8 524 379 1 0,8 476 838 0,5 426 787 0,2 464 935
3 0,5 550 823 0,8 476 838 1 0,8 424 351 0,5 288 800
4 0,1 910 331 0,5 426 787 0,8 424 351 1 0,8 366 847
5 0,366 206 0,2 464 935 0,5 288 800 0,8 366 847 1
При размещении приемников в продуктивном пласте для выделения полезного сигнала следует применять алгоритмы подавления слабо коррелируемого шума.
Для задач локации эндогенных сейсмических событий на основе решения обратных кинематических задач предпочтительно устанавливать сейсмические приемники вне продуктивного пласта.
Работа выполнена при поддержке Министерства образования и науки РФ и Российского фонда фундаментальных исследований (грант № 12−05−31 493).
Список литературы
1. Садовский М. А. Естественная кусковатость горной породы //ДАН СССР. 1979. Т. 247, № 4. С. 829−831.
2. Курленя М. В., Миренков В. Е., Сердюков С. В. Взгляд на природу напряженно-деформированного состояния недр Земли и техногенные динамические явления // Горный информационно-аналитический бюллетень. 2008. № 8. С. 5−19.
3. Родионов В. Н., Сизов И. А., Цветков В. М. Основы геомеханики. М., 1986.
4. Сердюков С. В., Курленя М. В. Механизм стимуляции добычи нефти сейсмическими полями малой интенсивности // Акустический журнал. 2007. Т. 53, № 5. С. 703 — 714.
5. Курленя М. В., Сердюков С. В. Низкочастотные резонансы сейсмической люминесценции горных пород в вибросейсмическом поле малой энергии // ФТПРПИ. 1999. № 1. С. 3−7.
6. Аки К., Ричардсон П. Количественная сейсмология. Теория и методы. Т. 1. М., 1983.
7. Николаевский В. Н. Механизм вибровоздействия на нефтеотдачу месторождений и доминантные частоты // ДАН СССР. 1989. Т. 307, № 3. С. 570−575.
8. Ризниченко Ю. В. Проблемы сейсмологии. М., 1985.
Об авторах
Валерий Всеволодович Сказка — д-р физ. -мат. наук, вед. науч. сотр., Институт математики им. С. Л. Соболева СО РАН, Новосибирск.
E-mail: vskazka@ngs. ru
Сергей Владимирович Сердюков — д-р техн. наук, гл. науч. сотр., Балтийский федеральный университет им. И. Канта, Калининград.
E-mail: ss3032@yandex. ru
Александр Сергеевич Сердюков — канд. физ. -мат. наук, мл. науч. сотр., Институт горного дела СО РАН, Новосибирск.
E-mail: aleksanderserdyukov@yandex. ru
About authors
Dr Valeriy Skazka — prof., senior research fellow, Sobolev Institute of Mathematics (Siberian Branch of the Russian Academy of Sciences), Novosibirsk.
E-mail: vskazka@ngs. ru
Dr Sergey Serdyukov — prof., principal research fellow, I. Kant Baltic Federal University, Kaliningrad.
E-mail: ss3032@yandex. ru
139
Alexandr Serdyukov — PhD, junior research fellow, Institute of Mining (Siberian Branch of the Russian Academy of Sciences), Novosibirsk.
E-mail: aleksanderserdyukov@yandex. ru

ПоказатьСвернуть
Заполнить форму текущей работой