Исследование распространения волн в нецилиндрических скважинах в анизотропной породе методом спектральных элементов

УДК 534.213.4

Г. С. Щелик

Московский физико-технический институт (государственный университет) Московский научно-исследовательский центр Шлюмберже

Исследование распространения волн в нецилиндрических скважинах в анизотропной породе методом спектральных элементов

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

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

1. Введение

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

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

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

вдоль направления минимальных напряжений. В работах [6,7] было проведено исследование влияния эллиптической формы скважины на дисперсию основных мод с применением метода граничных элементов и конечно-разностного моделирования в двухмерном случае. Другой подход, основанный на теории возмущений, был впервые предложен в работе [8] для оценки изменений в осесимметричных модах. Метод позволяет производить быстрые расчеты дисперсионных кривых нормальных мод, когда параметры возмущенной модели не сильно отличаются от значений известного начального решения. При этом было показано, что границы применимости теории возмущений могут быть значительно расширены при более удачном выборе этого начального приближения [9].

Анизотропия горных пород может являться следствием ряда факторов. Некоторые виды известняков и сланцев (например, Austin chalk и Bakken shale) обладают внутренней анизотропией и могут считаться однородными на масштабах, характерных для каротажа. Также анизотропия может быть вызвана горизонтально или наклонно расположенными слоями с толщинами, существенно меньшими, чем характерная длина волны источника. Кроме того, наличие трещин в окрестности ствола скважины может нарушать симметрию при распространении волн. За последние десятилетия было проведено множество исследований влияния анизотропии на данные сейсмических и акустических измерений. Уже в первых работах по моделированию распространения волн в скважинах в анизотропной породе было показано существование трех основных головных волн, соответствующих продольной (qP) и двум поперечным (SH и qSV) волнам [10]. Однако большая часть распространяющейся энергии содержится в собственных нормальных модах скважины. Достаточно много внимания было уделено исследованию дисперсионных характеристик волн в анизотропной породе, возникающих при проведении монопольного и дипольного каротажа [11,12]. К слабоанизотропным породам также может быть применен метод теории возмущений [8,13]. В последнее время для получения достоверной картины физических полей все чаще используются конечно-разностные методы [14,15], которые имеют широкие возможности по моделированию любого вида анизотропии, хотя и требуют больших вычислительных затрат.

В данной работе проводится моделирование волнового поля методом спектральных элементов (SEM), который ранее успешно применялся для расчета задач геофизики [16] и моделирования акустического каротажа [17]. Численный метод основан на разновидности метода конечных элементов высокого порядка. Реализация [16] предназначена для решения уравнений упругого анизотропного тела с возможностью включения заполненных жидкостью областей сложной геометрии. Метод позволяет с высокой точностью моделировать волны как в однородных областях, так и поверхностные волны на границах раздела сред, аккуратное воспроизведение которых играет ключевую роль при решении. В работе проведена необходимая адаптация алгоритма к рассматриваемым задачам, разработан блок построения сеток и задания свойств среды, написаны необходимые алгоритмы обработки данных. В тестовых расчетах трассы с приемников получены для случая дипольного источника в скважине в изотропных и трансверсально-изотропных породах. Рассматриваются скважины эллиптической и более сложной геометрии. Ось симметрии анизотропной породы имеет произвольные углы ориентации по отношению к стволу и направлению деформации скважины. Для получения дисперсионных кривых нормальных мод из сигнатур приемников используется модификация метода Пропп [18].

Структура статьи следующая. В параграфе 2 приведены сведения из теории распространения упругих волн в анизотропных средах и скважинах, необходимые для описания и обсуждения основных результатов работы. В разделах 3 и 4 проводится оценка точности расчетного метода и алгоритмов обработки на основе сравнения модельных дисперсионных кривых с доступными аналитическими решениями для цилиндрических скважин и результатами, полученными ранее методом граничных элементов. Ключевые результаты по исследованию дисперсии изгибных мод в нецилиндрических скважинах в анизотропной породе обсуждаются в разделе 5.

2. Теория

2.1. Волны в однородной анизотропной среде

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

ри Т

V- Т + {, С : Vи,

где и = (и\,и2,из) — вектор смещений, р — плотность, £ — вектор внешних сил, Т — тензор напряжений, имеющий линейную связь с малыми деформациями, определяемую тензором

СС ства породы и в случае общей анизотропии содержит 21 независимую компоненту.

Рассмотрим распространение плоских волн в однородной анизотропной породе. Для этого подставим в систему (1) решение для плоской волны вида и^ = А^ ехр(шщх1/У), где щ — компоненты вектора направления распространения, А^ — компоненты единичного вектора поляризации, Х[ — координаты. Тогда уравнения относительно фазовой скорости V можно записать в виде

рУ2 = Ау АкСутк1 птп1, ], к,т,1 = 1, 2,3.

Как обычно, по всем повторяющимся индексам производится суммирование. Полученное выражение известно как уравнение Кельвина-Кристофеля [19]. Оно представляет собой задачу на собственные значения и имеет три независимых решения, соответствующие трем типам плоских волн с различной поляризацией. Компоненты вектора поляризации А^ могут быть найдены из (2) для каждого собственного значения V2. По характеру колебаний эти решения принято соотносить с продольной (Р-волна) и двумя поперечными волнами (б'-волны). При этом важно отметить, что в случае общей анизотропии волны не могут иметь чисто продольную или поперечную поляризации.

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

С¿7 =

где С66 = (Си - С\2)/2.

В такой породе поляризация одной из поперечных волн, которую принято обозначать БН, будет всегда перпендикулярна направлению распространения волны. Две другие волны для большинства направлений будут иметь квазипродольную и квазипоперечную поляризации и обозначаются дР и дБУ соответственно. Фазовые скорости этих волн в силу симметрии могут быть вычислены для любого направления аналитически [19].

Для оценки степени анизотропии породы часто удобно пользоваться параметрами, предложенными Томсеном [20]:

(Си С\2 С\з 0 0 0

С12 с 11 С\з 0 0 0

С\з С\з Сзз 0 0 0

0 0 0 С44 0 0

0 0 0 0 С44 0

V 0 0 0 0 0 Сбб)

С\\ — С33

7

Сб6 — С44

6 =

(С\з + С44)2 — (Сзз — С44)2

2Сзз ' 2С'44 2С'зз

от направления распространения. В слабоанизотропных породах (при е,^,5 < 0.2) собственные значения и собственные векторы уравнения (2) незначительно отличаются от соответствующих значений, получаемых в изотропном случае. Такие породы могут быть аппроксимированы эффективными изотропными моделями [13].

2.2. Волны в скважине

Поиск решения (1) для системы, включающей скважину, является более сложной задачей. Как правило, в задачах, связанных с каротажем, используют математическую модель бесконечной цилиндрической скважины в изотропной породе. В предположении идеально упругой твердой породы, заполненной слабосжимаемой жидкостью, такая постановка допускает аналитическое решение. Аналитическое решение для анизотропного случая найдено только для трансверсально-изотропной породы со скважиной, ориентированной в направлении оси симметрии (УТ1) [11,12].

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

щ = и^гу^+^е-^, з = 1,2,3.

При подстановке (3) в граничные условия получим дисперсионное уравнение, которое определяет свойства распространяющихся в скважине волн. Дисперсионное уравнение имеет бесконечное количество решений, которые принято классифицировать двумя целыми числами тип, определяющими степень азимутальной симметрии и порядок функций Бесселя, входящих в решение для функций и™(г) [9].

Нормальные моды с т = 0 являются полностью симметричными, поле смещений и давления в жидкости не зависит от азимута. Первые два решения имеют названия трубной волны, или волны Стоунли (п = 0), и псевдорэлеевской волны (п = 1). Волна Стоунли локализована на поверхности раздела жидкой и твердой фаз и вызывает относительное движение флюида в пористой среде. Поэтому по скорости её распространения можно определять проницаемость окружающей породы. Следует отметить, что трубная волна существует на любых частотах, в том числе и очень низких. В данной работе влияние формы скважины и анизотропии породы на симметричные моды подробно не рассматривается.

Моды ст = 1 принято называть изгибными волнам,и. Они имеют ярко выраженную дисперсию, а скорость их распространения не превышают скорости поперечных волн в среде. В симметричной скважине в изотропной среде каждая из волн с определенной дисперсионной кривой является суперпозицией двух решений с одинаковыми скоростями [8]. При нарушении симметрии задачи, как будет видно в дальнейшем, происходит разделение этих волн и дисперсионных кривых на нечетные и четные ветви. При некоторой частоте (порядка 1-2 кГц для первой изгибной моды) дисперсионные кривые асимптотически выходят на значение скорости поперечной волны в породе, а энергия волны резко уменьшается. Это свойство используется в каротаже для определения скорости поперечной волны в породе, особенно в случаях, когда это невозможно сделать посредством регистрации головных волн. Следуя принятым в геофизике обозначениям, далее также будем называть изгибные моды дипольными, по названию типа акустического источника, которым они возбуждаются.

Рассмотрение всех мод со всеми типами симметрии и порядковыми номерами не требуется, поскольку на практике источники акустических колебаний возбуждают лишь моды с т = 0,1, 2 (так называемые монопольный, дипольный и квадрупольный источники). Также частотный диапазон большинства источников находится в пределах от 1 до 15 кГц и, как правило, не затрагивает далее второй моды по радиусу (п ^ 2).

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

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

3. Нецилиндрические скважины в изотропной породе

В качестве первого этапа проверки точности расчетного метода и алгоритмов обработки были рассмотрены скважины с нецилиндрической геометрией в изотропной породе. Методом спектральных элементов (SEM) было проведено моделирование распространения волн от дипольного источника в двух поляризациях в скважинах эллиптической (рис. 1) и несимметричной формы (замочная скважина) (рис. 2). При построении модели используются 20-узловые элементы второго порядка кривизны поверхности, функции внутри элементов интерполируются с полиномами 8 порядка.

2а

Рис. 1. Геометрия скважины эллиптической фор- Рис. 2. Геометрия скважины несимметричной мы с разбиением на элементы формы с разбиением на элементы

Т а б л и ц а 1

Свойства изотропных материалов, используемых в расчетах

Название материала Плотность Vp Vs кг/м3 м/с м/с

Жидкость 1000 1500 -

Быстрая порода 2160 4848 2601

Медленная порода 2000 2545 1018

Дисперсионные кривые для изгибных мод в эллиптических скважинах с отношениями большой и малой полуосей ( = а/Ъ = 1.5 и ( = 2 были сопоставлены с результатами, полученными методом граничных элементов в работе [6]. В сравнении использовались данные для мод с продольной поляризацией. В качестве наполняющей скважину жидкости взята чистая вода, значения свойств быстрой и медленной пород, используемых при вычислениях, приведены в табл. 1. Результаты сравнения дисперсионных кривых, полученных двумя методами, приведены на рис. 3 и рис. 4. Кривые имеют хорошее совпадение почти во всем диапазоне частот. Относительное расхождение результатов в основном диапазоне частот

не превышает 0.5% в медленной породе и 1% в быстрой. Незначительное повышение этого показателя на низких частотах связано с ростом погрешности вычислений этими методами в этой области.

700

03

о £

да ^

в

03

£

Четные

х - 10 х 10 см ое- 15 х 10 см   - 20 х 10 см

Нечетные

4 6

Частота, кГц

1120

х - 10 х 10 см ое- 15 х 10 см   - 20 х 10 см

Четные

2 3 4

Частота, кГц

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

Рис. 4. Сравнение результатов моделирования (маркеры) с дисперсионными кривыми, полученными методом граничных элементов (сплошные линии), в эллиптических скважинах различного размера в медленной породе

03 щ

о £

да ^

в

03

I

700

650

600

550

500

450

400

350

• - Несимметричная  - Эллиптическая.

20 х 20 см •

10 х 10 см

CJ

03 щ

о £

да ^

в

03

I

46 Частота, кГц

350

4 6

Частота, кГц

Рис. 5. Дисперсионные кривые дипольных мод Рис. 6. Дисперсионные кривые дипольных мод в несимметричной скважине в быстрой породе, в несимметричной скважине в быстрой породе. Продольная поляризация Поперечная поляризация

1

2

8

5

2

8

2

8

Эллиптическая форма является достаточно простой и распространенной нецилиндрической геометрией для скважин. Однако в реальных условиях набор возможных геометрий гораздо богаче. Для оценки влияния небольших отклонений от идеальной формы поперечного сечения на дисперсию дипольных мод был проведен ряд расчетов для сильно искаженной скважины. В качестве такой формы было взято сечение, ограниченное кривой вида г (в) = го + ri [(1 — {)/(1 — £ cos &)]р, где значения параметров г0 = 10 см, г1 = 20 см, £ = 0.65 и р = 8. Форма напоминает замочную скважину (рис. 2). Нарушения поперечного

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

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

¡1 1000

980

2 3 4 Частота, кГц

»Т

03 щ

о £

да ^

в

03

• - Несимметричная  - Эллиптическая

-20 х 20 см

¡1 1000

980

2 3 4 Частота, кГц

Рис. 7. Дисперсионные кривые дипольных мод в несимметричной скважине в медленной породе. Продольная поляризация

Рис. 8. Дисперсионные кривые дипольных мод в несимметричной скважине в медленной породе. Поперечная поляризация

1

6

0

1

5

5

4. Цилиндрические скважины в анизотропной породе

В анизотропной породе скорость волны зависит от направления её распространения и поляризации, поэтому собственные моды скважины могут испытывать расщепление на несколько волн, как и в случае нецилиндрических скважин. В расчетах использовалась трансверсально-изотропная модель материала, часто применяемая для описания осадочных горных пород. Ранее было показано [15], что вид дисперсионных кривых дипольных мод в таких породах определяется взаимной ориентацией оси симметрии материала и оси скважины, характеризуемой углом в (рис. 9). Если эти направления не совпадают {в > 0°), то каждая из дисперсионных кривых для дипольных мод (т = 1) разделится на две ветви, которые в низкочастотном пределе выходят на значения скоростей поперечных волн в среде в различных направлениях (дБУ- и БН- волны). Физически это эквивалентно существованию двух схожих по характеру колебаний волн с ортогональными направлениями поляризации, ориентированными вдоль главных направлений анизотропной породы. Если

направление ствола скважины совпадает с осью симметрии трансверсально-изотропной породы (в = 0°), то скорость поперечных волн в плоскости сечения в силу симметрии не зависит от направления колебаний, и структура дипольных мод имеет вид, аналогичный изотропному случаю. Помимо в положение оси симметрии определяется азимутальным углом поворота р (рис. 9). Угол характеризует ориентацию оси анизотропии относительно направления излучения источника колебаний или геометрии скважины.

Рис. 9. Схема ориентации осей симметрии трансверсально-изотропной породы

Для тестовых расчетов были взяты две трансверсально-изотропные породы: известняковая Austin Chalk и сланцеватая Cotton Valev Shale. Эти материалы ранее использовались для анализа в ряде работ (например, [15,20]). Упругие параметры и плотности приведены в табл. 2. В качестве жидкости в стволе были взяты параметры воды с плотностью 1000 кг/м3 и скоростью звуковых волн 1500 м/с. По соотношению скорости звука в жидкости и характерной скорости поперечных волн в Austin Chalk и Cotton Valev Shale, эти породы классифицируется как медленная и быстрая соответственно.

Т а б л и ц а 2

Свойства анизотропных материалов, используемых в расчетах

Название материала Плотность Vp V0 e 7 S кг/м3 м/с м/с

Austin Chalk 2200 2522 1044 0.2857 0.1458 0.2241

Cotton Valev Shale 2640 4721 2890 0.1350 0.1800 0.2050

В этом исследовании будем рассматривать частный случай наклонной трансверсаль-ной анизотропии с в = 60°. Ранее было показано, что в слабоанизотропных породах дисперсионные кривые могут быть аппроксимированы эффективными изотропными моделями [13]. Значения упругих постоянных (параметры Ламе) эквивалентной изотропной модели для кривой, соответствующей горизонтальной поляризации (БН), определяются выражениями ^вн = рУ'н и ^ви = р{Удр — 2^82н). Аналогичные выражения для вертикально-поляризованной волны {двУ) имеют вид р.дву = Р^'у и ^ву = р(УдР — ). На рис. 10 и рис. 11 приведены дисперсионные кривые, полученные из данных расчета для анизотропной породы в цилиндрической скважине. Для сравнения на графики нанесены кривые для эффективных изотропных моделей (сплошные линии), построенные на основании аналитических оценок фазовых скоростей. Результаты подтверждают наличие в обоих случаях двух типов волн медленной и быстрой. Как можно заметить, дисперсионные кривые для медленных волн лучше приближаются изотропной моделью, но даже для быстрых волн расхождение в скоростях не превышает 2%, что обычно является достаточным для корректной работы алгоритмов обработки для цилиндрических скважинах. Но, как будет показано в следующем параграфе, при искажении формы поперечного сечения вид дисперсионных кривых будет существенно отличаться от приближения изотропными моделями.

и

К

о S

да

<3

03

2 4 6 8

Частота, кГц

Рис. 10. Дисперсионные кривые для быстрой и медленной изгибной волны в быстрой трансверсально-изотропной породе. Сплошные линии обозначают приближение эффективными изотропными моделями для каждой из волн

1020 1000

§г 980

А

940

0

а

1 900

¡1 880

О - Медленная   - Быстрая

У " 0<Х><>Х пП /  о ✓ -  / У пП0-nnD о ^^  ^^ VqSV .

о / о / о / о / -  / / "  Vsh '

, , , , ,

860

1 2 3 4 5 6 7 Частота, кГц

Рис. 11. Дисперсионные кривые для быстрой и медленной изгибной волны в медленной трансверсально-изотропной породе. Сплошные линии обозначают приближение эффективными изотропными моделями для каждой из волн

5. Нецилиндрические скважины в анизотропной породе

Деформации формы скважины в анизотропной породе, аналогично изотропному случаю, приводят к изменению вида дисперсионных кривых каждой из существующих ди-польных мод. На характер этих изменений помимо геометрических параметров скважины влияет также взаимная ориентация направлений деформации ствола и анизотропии породы. На рис. 13 приведены дисперсионные кривые для эллиптической скважины в медленной трансверсально-анизотропной породе (Austin Chalk) в случае, когда ось симметрии анизотропной породы повернута относительно направлений деформации скважины {<р = 45°). Сравнение с результатами расчетов для цилиндрических скважин в анизотропной породе показывает, что дисперсионные кривые лежат между решениями для скважин с радиусами соответствующими длинам большой и малой полуосей эллипса в анизотропной породе. В низкочастотной асимптотике полученные дипольные моды будут иметь скорость поперечных волн в породе, а на высоких частотах приближаться к решению для цилиндрической скважины с наиболее близким радиусом кривизны. При этом в зависимости от угла поворота р эффекты могут как суммироваться, так и вычитаться, тем самым влияя на удаление медленной и быстрой мод друг от друга в заданных пределах.

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

Для определения параметров анизотропии окружающей скважину породы используется метод кросс-дипольного каротажа. Он основывается на данных измерений приходящего сигнала от двух дипольных источников, ориентированных в перпендикулярных направлениях. Определение положения оси по азимуту р может быть выполнено на основе данных таких измерений с помощью алгоритма Alford rotation [21]. При моделировании в качестве входных данных для алгоритма используются результаты расчетов распространения волн для двух ортогональных начальных ориентаций дипольного источника, имитирующих реальные измерения в скважине. Если скважина в анизотропной породе не является цилиндрической, можно говорить о существовании двух независимых факторов, влияющих на

дипольные волны. С точки зрения алгоритмов обработки это эквивалентно влиянию некоторой дополнительной анизотропии в породе, связанной с искаженной формой скважины, что может приводить к ошибкам в определении направления истинной анизотропии. На рис. 1 1а I приведены результаты определения главных направлений для цилиндрических и эллиптических скважин в анизотропной породе. Как можно видеть, отклонение от истинного направления в 45° увеличивается с ростом параметра эллиптичности (. При этом в быстрой породе, как уже было сказано выше, даже при относительно небольших деформациях фактор геометрии будет являться определяющим для направлений распространения главных волн (рис. 14а-г).

700

со

о £

да §

со

4 6

Частота, кГц

1050

I 1000

со со о £

да §

со <и

950

900

850

20 х 20 см

хх£++

Х+++++++

• хх ++;.••• \хХх**х

+ • X •

• х • ,

10 10 см

х - Быстрая + - Медленная

1 2 3 4 5 6 Частота, кГц

Рис. 12. Сравнение дисперсионных кривых дипольных мод в быстрой анизотропной породе в эллиптической скважине 20 х 10 см (черные маркеры) с результатами для цилиндрических скважин радиусом 10 и 20 см

Рис. 13. Сравнение дисперсионных кривых дипольных мод в медленной анизотропной породе в эллиптической скважине 20 х 10 см (черные маркеры) с результатами для цилиндрических скважин радиусом 10 и 20 см

7

2

8

Рис. 14. Результаты определения главных осей анизотропии в эллиптических скважинах 15 х 10 20 х 10

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

6. Заключение

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

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

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

Работа выполнена на базе Московского научно-исследовательского центра Шлюмберже. Автор благодарен своему научному руководителю И. Л. Софронову за конструктивную критику и помощь в подготовке статьи, а также признателен Т. Жарникову и М. Шарара за консультации и помощь в проведении численных расчетов.

Литература

1. Cheng N., Cheng С. Н., Toksdz М. N. Borehole wave propagation in three dimensions // Journal of the Acoustical Society of America. — 1995. — V. 97, N. 6. — P. 3483-3493.

2. Wang Т., Tang X. Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach // Geophysics. — 2003. — V. 68, N. 5. — P. 1749-1755.

3. Bolshakov A. O., Zheng Y., Schmitt D. P. Study of Flexural-mode Dispersion With Anisotropic Formations and Non-circular Boreholes Using the Finite Element Method // 71st EAGE Conference & Exhibition. - 2009. - P. 8-11.

4. Mallan R. K., Torres-Verdín С., Ma J. Simulation of borehole sonic waveforms in dipping, anisotropic, and invaded formations // Geophysics. — 2011. — V. 76, N. 4.

5. Simsek E., Sinha В. K., Zeroug S., et aZ.Influence of breakouts on borehole sonic dispersions // 2007 SEG Annual Meeting. - 2007. - P. 313-317.

6. Randall C. J. Modes of noncircular fluid-filled boreholes in elastic formations // Journal of the Acoustical Society of America. - 1991. - V. 89, N. 3. - P. 1002-1016.

7. Randall C. J. Multipole acoustic waveforms in nonaxisvmmetric boreholes and formations // Journal of the Acoustical Society of America. — 1991. — V. 90, N. 3. — P. 1620-1631.

8. Ellefsen K. J., Cheng С. H., Toksdz M. N. Applications of perturbation theory to acoustic logging // Journal of Geophysical Research. — 1991. — V. 96. — P. 537-549.

9. Simsek E., Sinha В. K. Analysis of noncircular fluid-filled boreholes in elastic formations using a perturbation model // Journal of the Acoustical Society of America. — 2008. — V. 124, N. 1. - P. 213-7.

10. Kurkjian A. L., Chang S. K. Acoustic multipole sources in fluid-filled boreholes // Geophysics. - 1986. - V. 51, N. 1. - P. 148-163.

11. White J. E., Tongtaow C. Cylindrical waves in transversely isotropic media // Journal of the Acoustical Society of America. — 1981. — V. 70, N. 4. — P. 1147-1155.

12. Chan A. K., Tsang L. Propagation of acoustic waves in a fluid-filled borehole surrounded by a concentrically layered transversely isotropic formation // Journal of the Acoustical Society of America. - 1983. - V. 74, N. 5. - P. 1605-1616.

13. Sinha B. K., Norris A. N., Chang S. K. Borehole flexural modes in anisotropic formations // Geophysics. - 1994. - V. 59, N. 7. - P. 1037-1052.

14. Leslie H. D., Randall C. J. Multipole sources in boreholes penetrating anisotropic formations: Numerical and experimental results // Journal of the Acoustical Society of America. - 1992. - V. 91. — P. 12-27.

15. Sinha B. K., Simsek E., Liu Q. H. Elastic-wave propagation in deviated wells in anisotropic formations // Geophysics. - 2006. - V. 71, N. 6. - P. D191-D202.

16. Komatitsch D., Tromp J. Introduction to the spectral element method for three dimensional seismic wave propagation // Geophysical Journal International. — 1999. — V. 139. — P. 806-822.

17. Charara M., Vershinin A., Deger E., et al. 3D spectral element method simulation of sonic logging in anisotropic viscoelastic media // SEG Technical Program Expanded Abstracts 2011. - 2011. - P. 432-437.

18. Ekstrom M. P. Dispersion estimation from borehole acoustic arrays using a modified matrix pencil algorithm // Proceedings of ASILOMAR-29. - 1995. - P. 449-453.

19. Musgrave M. J. P. Crystal Acoustics: Introduction to the Study of Elastic Waves and Vibrations in Crystals. — San Francisco: Holden-Dav, 1970. — 288 p.

20. Thomsen L. Weak elastic anisotropv // Geophysics. — 1986. — V. 51, N. 10. — P. 19541966.

21. Alford R. M. Shear data in the presence of azimuthal anisotropv: Dillev Texas // 1986 SEG Annual Meeting. — Society of Exploration Geophvsicists, 1986. — P. 476-479.

Поступим в редакцию 29.10.2013.