Влияние электрической анизотропии горных пород на электромагнитное поле в скважине
Решение прямых задач электродинамики о поле наклонного переменного магнитного диполя и линии переменного тока в присутствии модели анизотропного пласта конечной мощности. Определение электрических параметров горных пород при измерениях в скважинах.
| Рубрика | Геология, гидрология и геодезия |
| Вид | автореферат |
| Язык | русский |
| Дата добавления | 02.03.2018 |
| Размер файла | 471,9 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
На правах рукописи
УДК 550.8.013:550.832.7
25.00.10 - геофизика, геофизические методы поисков полезных ископаемых
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора физико-математических наук
ВЛИЯНИЕ ЭЛЕКТРИЧЕСКОЙ АНИЗОТРОПИИ ГОРНЫХ ПОРОД НА ЭЛЕКТРОМАГНИТНОЕ ПОЛЕ В СКВАЖИНЕ
Каринский А.Д.
Москва 2008
Работа выполнена в Российском государственном геологоразведочном университете имени Серго Орджоникидзе (РГГРУ).
Официальные оппоненты:
академик РАН, доктор технических наук, профессор Эпов Михаил Иванович.
член кор. РАЕН, доктор физико-математических наук, профессор Юдин Михаил Николаевич.
доктор технических наук, профессор Безрук Игорь Андреевич.
Ведущая организация
Геологический факультет Московского Государственного Университета имени М. В. Ломоносова.
Защита состоится “18” декабря 2008 г. в 1500 часов на заседании Диссертационного Совета Д.212.121.07 при Российском государственном геологоразведочном университете имени Серго Орджоникидзе.
Адрес: 117485, ул. Миклухо-Маклая, д. 23, РГГРУ.
С диссертацией можно ознакомиться в библиотеке Российского государственного геологоразведочного университета.
Автореферат разослан «____» _______________ 2008 г.
И. о. ученого секретаря Диссертационного совета Д.212.121.07 при РГГРУ
доктор технических наук, профессор Бондаренко В. М.
1. Общая характеристика работы
Состояние проблемы. Актуальность темы.
В последние годы проявляется все больший интерес к изучению анизотропии физическ их параметров горных пород и ее влияния на результаты геофизических исследований. Анизотропия, то - есть различие в значениях того или иного физического параметра вещества по разным направлениям, свойственна в первую очередь многим осадочным горным породам. Это связано с особенностями формирования осадочных пород, обычно сопровождающегося многократными изменениями режима осадконакопления с соответствующими изменениями состава и дисперсности материала осадков. Этот процесс часто приводит к формированию слоистой или тонкослоистой структуры осадочных толщ и анизотропии их физических параметров.
Анизотропия проявляется, в частности, в электрических параметрах горных пород, характеризующих их способность проводить электрический ток или поляризоваться в электрическом поле. Наиболее часто применяемой для характеристики электрических параметров анизотропных пород моделью среды является среда с осевой анизотропией (одноосно - анизотропная или трансверсально- изотропная среда). Каждый электрический параметр такой среды характеризуют два значения: по оси анизотропии n (ориентированной по нормали к напластованию пород) и по любому, ортогональному этой оси направлению t. В частности, удельное электрическое такой среды характеризуют два значения: поперечное удельное электрическое сопротивление n по направлению n и продольное удельное электрическое сопротивление t по направлению t. Параметры электрической анизотропии содержат важную геолого-геофизическую информацию о строении и свойствах пород слагающих анизотропную толщу, так как продольное и поперечное удельные электрические сопротивления по разному зависят от характера насыщения пород, их фильтрационных свойств и других факторов.
По сравнению с наземными геофизическими измерениями условия при геофизических исследованиях скважин (ГИС) являются более благоприятными для изучения электрических параметров анизотропных горных пород, так как в этом случае возбудители и измерители поля находятся непосредственно в исследуемой толще. Однако, применяемые ныне методы ГИС, не способны обеспечить получение достаточно полной информации об электрических параметрах анизотропных пород, в первую очередь - о значении n. В обычных условиях результаты измерений в таких широко применяемых методах, как индукционный каротаж (ИК), каротаж кажущихся сопротивлений (КС) и некоторых других, зависят, в основном, от влияния одного параметра анизотропных пород - t . В то же время, параметр n является более информативным для определения, например, такой важнейшей характеристики анизотропной толщи, как тип насыщения.
Следует заметить, что определение параметров анизотропии горных пород с помощью наблюдений в скважинах является сложной задачей. Для ее решения необходимо определить оптимальные типы возбудителей поля и измеряемые компоненты поля, наиболее тесно связанные с каждым из параметров анизотропных пород, установить характер этой связи. Необходимо также исследовать влияние реальных условий измерений в скважинах (присутствие скважины, зоны проникновения, границ анизотропных и изотропных пластов, в некоторых случаях - буровой колонны или обсадной трубы) на характер связи измеряемых величин с параметрами анизотропной среды.
Несмотря на научную и прикладную важность определения параметров электрически анизотропных горных пород с помощью измерений в скважинах этой проблеме до последнего времени не уделялось должного внимания. В настоящей работе делается попытка частично восполнить этот пробел путем теоретического анализа проблемы.
Отметим также, что выполненные ранее исследования по теории электромагнитных методов ГИС в анизотропных средах касались, как правило, связей характеристик электрического или электромагнитного поля с удельными электрическими сопротивлениями анизотропных сред. Но на характеристики высокочастотного электромагнитного поля в проводящей среде оказывает влияние не только , но и диэлектрическая проницаемость е. Теоретическое изучение связи характеристик высокочастотного электромагнитного поля (и определяемых по этим характеристикам значений кажущейся диэлектрической проницаемости) с диэлектрическими проницаемостями моделей анизотропных сред также нашло отражение в этой работе.
В работе приведены также результаты исследований по проблеме, связанной с правомерностью замены (при моделировании и интерпретации результатов измерений) макроанизотропных горных пород микроанизотропными средами с теми же, что у пород, значениями n и t.
Объектом исследований является электромагнитное поле различных его сторонних возбудителей в моделях анизотропной среды и связь характеристик поля с параметрами такой среды.
Цель работы - развитие теоретических основ электромагнитного каротажа анизотропных сред; теоретический анализ различных методов возбуждения и измерения электромагнитного поля в скважине, пройденной в анизотропных горных породах и выявление способов определения параметров электрически анизотропных пород по этим данным.
Основные задачи исследований. Для достижения цели исследований потребовалось решение следующих задач.
1. Получение аналитических выражений для компонент электромагнитного поля различных его сторонних возбудителей для модели однородной анизотропной среды. Вывод асимптотических выражений для компонент электромагнитного поля в ряде практически значимых случаев и их анализ.
2. Аналитическое решение прямых задач электродинамики о поле различных его сторонних возбудителей для кусочно-однородных моделей среды с коаксиальными цилиндрическими границами. Однородные части таких моделей могут отвечать анизотропным горным породам и зоне проникновения в проницаемых породах, изотропным скважине, буровой колонне, обсадной трубе. Составление программ для численных расчетов на основе полученных решений.
3. Аналитическое решение прямых задач электродинамики о поле наклонного переменного магнитного диполя и линии AB переменного тока в присутствии модели анизотропного пласта конечной мощности. Составление программ для численных расчетов на основе полученных решений.
4. Проведение численных расчетов на основе аналитических решений прямых задач для модели однородной анизотропной среды и 1D- моделей анизотропной среды. Анализ полученных результатов расчетов.
5. Разработка алгоритмов численных расчетов электромагнитного поля различных его сторонних возбудителей в осесимметричных кусочно-однородных моделях анизотропной среды с коаксиально- цилиндрическими и плоско- параллельными границами. Примененные алгоритмы расчетов основаны на методе конечных разностей. Составление программ для расчетов на ЭВМ, оценка погрешности численных расчетов.
6. Математическое моделирование для соответствующих условиям геофизических исследований скважин 2D- моделей анизотропной и изотропной среды при различных способах возбуждения электромагнитного поля. Сравнительный анализ результатов моделирования для микроанизотропных и макроанизотропных моделей среды.
7. Обобщение полученных результатов моделирования с целью обоснования методических рекомендаций по определению электрических параметров анизотропных горных пород при измерениях в скважинах и применении различных возбудителей электромагнитного поля и приемников (датчиков), позволяющих измерить различные характеристики поля.
Метод исследования. На разных этапах работы применялись различные методы исследований. Сначала были получены аналитические решения прямых задач электродинамики для модели однородной анизотропной среды и 1D- моделей анизотропной среды для различных сторонних возбудителей поля, проанализированы полученные решения и результаты численных расчетов. Затем с целью численного решения прямых задач (для приближающихся к реальным условиям измерений в скважине) осесимметричных 2D- моделей анизотропной среды разработаны алгоритмы и программы для ЭВМ, основанные на методах конечных разностей. Проведены математическое моделирование для 2D- моделей среды при различных возбудителях поля и анализ результатов моделирования.
Научная новизна работы. Личный вклад. Основная часть исследований, результаты которых приведены в настоящей работе, является оригинальной и выполнена лично автором данной работы. Среди новых результатов, полученных при исследованиях, можно отметить следующие.
1. Впервые получены аналитические решения прямых задач электродинамики для некоторых сторонних возбудителей поля в модели однородной анизотропной среды. Для такой модели среды получены выражения для поля тороидальной антенны бесконечно - малых и конечных размеров, бесконечно-длинного кабеля (БДК), для электрической составляющей поля переменного магнитного диполя и плотности источников кулоновой составляющей этого поля.
2. Получены аналитические решения прямых задач электродинамики для моделей анизотропной среды с цилиндрическими границами. Эти модели соответствуют условиям измерений в скважинах при неограниченной мощности пластов. Решения получены для случаев, когда сторонними возбудителями поля являются элементарная тороидальная антенна, линия AB переменного тока, кабель с точечным электродом, БДК, тороидальная антенна конечных размеров.
3. На основе метода конечных разностей разработаны алгоритмы численного решения осесимметричных 2D- прямых задач электродинамики для соответствующих условиям ГИС моделей анизотропной среды с коаксиальными цилиндрическими и плоско-параллельными границами для различных сторонних возбудителей поля.
4. Анализ полученных выражений и результатов математического моделирования позволил выявить ряд неизвестных ранее особенностей влияния параметров анизотропной среды на характеристики электромагнитного поля при различных способах его возбуждения, дать физическое истолкование полученным данным моделирования и обоснование новым методикам ГИС, применение которых может быть перспективным при изучении анизотропных пород. В частности, в рамках решения проблемы определения значений n и коэффициента анизотропии =(n/t)1/2 установлено следующее:
а) плотность электрических зарядов (и их поле), индуцируемых полем переменного магнитного диполя, в ближней зоне зависит лишь от одного электрического параметра среды - коэффициента анизотропии ;
б) в однородной анизотропной среде поле тороидальной антенны испытывает сильное влияние n и , но при расположении стороннего возбудителя в изотропной скважине поля элементарной тороидальной антенны и переменного электрического диполя испытывают одинаковое влияние окружающей скважину среды, в частности ее анизотропии;
в) при применении в условиях каротажа в процессе бурения тороидальных антенн и при соответствующем выборе измеряемых характеристик и длин зондов результаты измерений существенно зависят от величины n;
г) в присутствии границ анизотропных пластов индукционная составляющая ImEz низкочастотного электрического поля кабеля с токовым электродом (а при применении БДК - обе составляющие- ReEz и ImEz) вблизи этих границ испытывают значительное влияние n, а кривые для этих компонент поля против пластов имеют простую форму.
Основные защищаемые положения.
1. В однородной анизотропной среде (без учета влияния скважины) для большинства применяемых зондов электрического и электромагнитного каротажа при их ориентации по оси анизотропии результаты измерений зависят от продольного удельного электрического сопротивления t. Для получения информации о других параметрах анизотропной среды (n, л) предложены специальные зонды, отличающиеся способами возбуждения и измерения поля.
2. Наличие скважины, как правило, не ведет к существенному изменению характера влияния параметров окружающей ее анизотропной среды на характеристики поля различных возбудителей по сравнению с однородной средой. Однако при возбуждении поля тороидальной антенной условия измерений в скважине кардинальным образом изменяют это влияние.
3. При применении зондов с тороидальными антеннами в условиях каротажа в процессе бурения и при определенном выборе длин зондов и измеряемых характеристик поперечное удельное электрическое сопротивление n оказывает существенное влияние на получаемые результаты.
4. В приближающихся к реальным условиям при геофизических исследованиях скважин 2D- моделях анизотропной среды с цилиндрическими и плоскими границами при определенных способах возбуждения поля влияние n может быть существенно выше, чем в однородной анизотропной среде или в модели анизотропного пласта неограниченной мощности в присутствии скважины.
Практическая значимость работы. Анализ полученных решений и данных моделирования позволил обосновать новые методики ГИС, применение которых позволяет получить практически важные данные о параметрах анизотропных пород. На основе полученных решений прямых задач для различных сторонних возбудителей поля в 1D- моделях анизотропной среды и разработанных на основе методов конечных разностей алгоритмов численных расчетов для 2D- моделей анизотропной среды создано программное обеспечение, позволяющее проводить расчеты характеристик поля различных сторонних возбудителей. Эти программные средства могут послужить при разработке конкретных методик скважинных измерений с целью изучения параметров электрически анизотропных пород и при интерпретации результатов таких измерений.
Апробация работы. Основные результаты по мере их получения обсуждались на ряде международных конференций и симпозиумов. Среди них- III Латиноамериканская геофизическая конференция и IX Симпозиум Мексиканского общества геофизиков (III Conferencia Latinoamericana de Geofisica y IX Simposium de Geofisica), Бийа-Эрмоса, Табаско, Мексика, 2000г.; 42-ой ежегодный симпозиум профессиональных исследователей- каротажников (Annual Logging Symposium, Society of Professional Well Log Analysts - SPWLA), Хьюстон, США, 2001г.; II международная конференция Геофизического Союза Мексики (UNION GEOFISICA MEXICANA), Пуэрта Байарта, Мексика, 2001г.; 43-ой ежегодный симпозиум профессиональных исследователей- каротажников (SPWLA), Ойсо, Япония, 2002г.; VI, VII, и VIII международные конференции «Новые идеи в науках о Земле», Москва, Россия, 2003, 2005 и 2007 гг.; конференция РГГРУ «Современные геофизические и геоинформационные системы», Москва, Россия, 2008г.
Публикации. По теме диссертации опубликовано 27 работ, в том числе 10 статей в рецензируемых научных журналах, 6 докладов.
Объем и структура работы.
Диссертационная работа состоит из введения, шести глав и заключения, содержит 228 страниц текста, 53 рисунка. Библиография содержит 120 наименований.
Работа выполнена на кафедре электрических, гравитационных и магнитных методов поисков и разведки месторождений полезных ископаемых Российского государственного геологоразведочного университета имени Серго Орджоникидзе (г. Москва). В работе также нашли отражение результаты исследований автора в период его работы в Мексиканском институте нефти (IMP, г. Мехико) в 1999-2002гг.
Автор выражает искреннюю благодарность своему учителю и старшему коллеге, профессору Д. С. Даеву, который на всех этапах исследований способствовал написанию этой работы и детально ознакомился с рукописью диссертации, сделав ряд ценных замечаний.
С чувством глубокой благодарности автор обращается к памяти профессора Л. М. Альпина, одного из основателей отечественной школы каротажников и электроразведчиков, общение и совместная работа с которым оказала на автора огромное влияние.
Автор благодарен также сотрудникам Мексиканского Института Нефти (IMP) А. А. Мусатову, В. А. Шевнину и P. Anguiano за плодотворные обсуждения некоторых результатов исследований.
2. Содержание работы
Введение.
Дан краткий обзор состояния проблемы. Охарактеризованы актуальность проблемы и цель исследований, основные задачи, научная и практическая значимость работы, представительность ее апробации. Представлены основные защищаемые положения.
В области изучения влияния удельных электрических сопротивлений анизотропных горных пород на стационарное электрическое поле или квазистационарное электромагнитное поле выполнены теоретические исследования рядом отечественных и иностранных ученых. Среди таких исследований - работы Л.М. Альпина, М.Н. Бердичевского, В.Р. Бурсиана, Л.Л. Ваньяна, В.В. Вержбицкого, Ю.А. Дашевского, А.Е. Кулинковича, А.С. Семёнова, А.И. Сидорчука, Е.В Чаадаева, С.М. Шейнмана, М.И. Эпова, B.I. Anderson, T.D. Barber, A.B. Cheryauka, S. Gianzero, T. Hagiwara, W.D. Kennedy, K.S. Kunz, M.G. Lьling, J.H. Moran, L.А. Tabarovsky, M.S. Zhdanov и другие. Обоснованию методов решения прямых задач электродинамики в анизотропных средах при соответствующих нормированиях (калибровках) электродинамических потенциалов были посвящены работы А.Н. Тихонова, Д.Н. Четаева, Л.Л. Ваньяна.
Для условий измерений в скважинах теоретические исследования, с целью изучения влияния электрической анизотропии горных пород на результаты таких измерений, как правило, проводились применительно к теории двух широко применяемых на практике методов геофизических исследований скважин (ГИС): каротажа кажущихся сопротивлений (КС) и индукционного каротажа (ИК). В качестве простейших сторонних возбудителей поля в теории этих методов рассматривают соответственно источник стационарного электрического поля в виде точечного электрода и такой возбудитель гармонически - меняющегося квазистационарного электромагнитного поля, как переменный магнитный диполь.
Исследованию влияния анизотропии горных пород на постоянное электрическое поле расположенного в скважине точечного электрода были посвящены работы Л.М. Альпина, В.В. Вержбицкого, Ю.А. Дашевского, А.Е. Кулинковича, А.И. Сидорчука, Л.А. Табаровского, Е.В. Чаадаева и др. Работы Л.М. Альпина, в частности, содержат идею о том, как можно получить данные о «вертикальном» (поперечном) удельном электрическом сопротивлении n анизотропных пород путем перемещения измерительного диполя MN в поле неподвижного удаленного от анизотропных пластов токового электрода А. В работах А.Е. Кулинковича, В.Н. Румянцева, А.И. Сидорчука, Е.В. Чаадаева было высказано утверждение, что наличие анизотропии в пластах может проявиться в усложнении формы диаграмм градиент- зондов каротажа КС.
В работах А.Д. Гайдаша, В.А. Пантюхина, К.Л.Санто, Е.В.Чаадаева, М.И. Эпова, S., Gianzero, S. Graciet, T. Hagiwara, K.S.Kunz, J.H. Moran, L.C. Shen и других ученых предметом исследований являлось влияние анизотропных сред на компоненты низкочастотного магнитного поля переменного магнитного диполя. В последнее время в США различными исследователями выполнен большой объем работ с целью обоснования применения многокомпонентной аппаратуры индукционного каротажа (Multi-Component Induction Logs) при изучении анизотропных горных пород, и разработки методики интерпретации результатов измерений с этой аппаратурой. Зондовое устройство этой аппаратуры состоит из нескольких различно ориентированных генераторных и измерительных переменных магнитных диполей, что позволяет измерять компоненты низкочастотного магнитного поля, возбуждаемого магнитными диполями с взаимно- ортогональными моментами. Эти работы в США проводятся в компаниях Шелл, Бейкер, Шлюмберже, в Университетах штатов Юта, Техас и др. В проведении этих исследований принимает участие большое число ученых, в частности, A. Bespalov, B. Corley, A. Gribenko, S. Fang, O. Fanini, S. Forgang, Guo-Zhong Gao, S.W. Itskovich, B. Kriegshauser, G. Merchant, J. Morrison, E. Quint, M. Rabinovich, L.A. Tabarovsky, C. Torres, L. Yu, Z. Zhang, M. Zhdanov и др. В этих исследованиях также принимают участие и российские ученые (М. И. Эпов).
Но эти работы, как и практически все выполнявшиеся до последних лет исследования с целью определения электрической анизотропии пройденных скважиной горных пород, ориентированы на изучение влияния удельных электрических сопротивлений анизотропной среды на поле его сторонних возбудителей, традиционно применяемых в электрических и электромагнитных методах геофизических исследований скважин. Такими возбудителями поля являются токовый электрод и электрический диполь, как источники стационарного электрического поля, а также переменные магнитный и электрический диполи, как источники квазистационарного электромагнитного поля.
Оригинальный подход к исследованию анизотропии горных пород был предложен S. Gianzero. В его работе, опубликованной в 1999 г., высказано утверждение, что зонд из соосных генераторного электрического диполя и измерительной тороидальной антенны должен быть эффективным при определении «вертикального» (поперечного) удельного электрического сопротивления n. Этот вывод был сделан на основе анализа выражений для азимутальной компоненты магнитного поля переменного электрического диполя, лежащего на оси анизотропии однородной среды.
Результаты наших исследований показали, что при геофизических измерениях в скважинах с целью изучения электрических параметров анизотропных горных пород может быть перспективным применение некоторых методик измерений, отличающихся от тех, которыми пользуются в настоящее время в методах ГИС. Это, в частности, связано с использованием таких возбудителей электромагнитного поля, как тороидальная антенна, протяженная линия переменного тока или кабель с переменным током; измерение электрических компонент поля, возбуждаемого переменным магнитным диполем и др.
Глава 1. Уравнения переменного электромагнитного поля в анизотропной среде. Электродинамические потенциалы.
Выражения, приведенные в этой главе, послужили основой при аналитических и численных решениях прямых задач электродинамики в следующих главах настоящей работы.
Раздел 1.1. Представлены уравнения электромагнитного поля и уравнения связи в немагнитной анизотропной среде. Даны основные определения. Детальнее рассмотрены уравнения гармонически меняющегося поля (и некоторые другие выражения) в одноосно - анизотропной среде. В частности, соответствующие направлениям n и t волновые числа kn,t в немагнитной среде определяет равенство:
,
где еn,t - диэлектрические проницаемости по направлениям n и t, i - мнимая единица, 0=1/(0с2)(10?9/36), Ф/м - электрическая постоянная, c3108, м/с - электродинамическая постоянная, равная скорости света в вакууме, 0=410?7, Гн/м - магнитная постоянная, =2f - круговая частота, f - частота колебаний. При записи выражений для гармонически меняющихся величин и полей временной фактор был выбран в виде: exp(?iщt), где t - время. В случае гармонически меняющегося поля коэффициент анизотропии =kt/kn. Для стационарного электрического поля или квазистационарного электромагнитного поля (ще0еn,t<<1/сn,t) коэффициент анизотропии , где уn=1/сn и уt=1/сt - удельные электропроводности по направлениям n и t.
Отмечено, что в анизотропной среде (в отличие от изотропной среды) объемные электрические заряды могут существовать и в том случае, когда такая среда - однородна. Если ось Z направлена по оси анизотропии n, то плотность объемных зарядов
где Ez - компонента электрического поля по оси анизотропии n. То есть плотность зарядов 0, если 1 и компонента поля по оси n меняется по направлению этой оси.
Раздел 1.2. Для гармонически меняющегося поля приведен вывод дифференциальных уравнений, которым в однородных областях анизотропной среды удовлетворяет векторный электродинамический потенциал электрического типа. Представлены выражения, связывающие магнитное () и электрическое () поля с потенциалом .
Раздел 1.3. Приведен вывод дифференциальных уравнений, которым (при калибровке Тихонова) удовлетворяет векторный электродинамический потенциал магнитного типа в однородной анизотропной среде. Даны выражения, связывающие составляющие и гармонически меняющегося электромагнитного поля с этим потенциалом.
Глава 2. Поле различных сторонних токов (генераторных антенн) в однородной анизотропной среде.
На основе приведенных в главе 1 уравнений для векторных электродинамических потенциалов электрического и магнитного типа и выражений, связывающих эти потенциалы с составляющими электромагнитного поля, получены выражения для компонент гармонически меняющегося поля в однородной одноосно-анизотропной среде. Они определяют компоненты поля различных его сторонних возбудителей: переменного электрического диполя, переменного магнитного диполя, тороидальной антенны бесконечно- малых разметов, линии AB переменного тока, бесконечно- длинного кабеля (БДК), линейной тороидальной антенны конечных размеров. Анализ полученных выражений и результатов численных расчетов позволил установить некоторые закономерности, проявляющиеся во влиянии параметров однородной анизотропной среды на разные компоненты электромагнитного поля различных его возбудителей.
Раздел 2.1. Приведены выражения для электрической составляющей поля произвольно ориентированного относительно оси анизотропии n (Z) переменного электрического диполя. Проанализированы выражения для компонент поля, определяющих результаты измерений при применении дипольно-осевых и дипольно-экваториальных зондов (установок) при их различной ориентации относительно оси n. В частности, как известно, при применении дипольно-осевого зонда (из генераторного и измерительного электрических диполей), ориентированного по оси анизотропии, компонента Ez электрического поля, которой пропорционально измеряемое напряжение в измерительном диполе MN, не зависит от волнового числа kn. При выполнении условий квазистационарности поля и ближней зоны в соответствии с парадоксом анизотропии (Бурсиан В.Р и др.) синфазная току I=I0?exp(?iщt) в генераторном электрическом диполе компонента ReEz пропорциональна продольному удельному электрическому сопротивлению t анизотропной среды. Запаздывающая по фазе на р/2 относительно I компонента ImEz при этих условиях не зависит от электрических параметров среды.
В разделе 2.2 получены и проанализированы выражения для компонент магнитного и электрического поля переменного магнитного диполя, а также - для плотности д объемных зарядов, индуцируемых электрической составляющей поля переменного магнитного диполя. Рассмотрены результаты численных расчетов для некоторых характеристик поля.
Если центр, магнитного диполя расположен в начале 0 декартовых координат с осью Z по оси анизотропии n, а момент диполя, лежит в плоскости y=0, то =Mx+Mz, где , - единичные векторы (орты). Из равенства (1) и выражений, полученных для компонент электрического поля переменного магнитного диполя, следует, что в такой системе координат плотность объемных электрических зарядов где . В ближней зоне () имеем: . То есть при возбуждении поля в однородной среде переменным магнитным диполем в ближней зоне плотность д объемных электрических зарядов (и их поле) зависит лишь от одного электрического параметра среды - коэффициента анизотропии л.
Известно, что в однородной изотропной среде в точке И на оси l магнитного диполя электрическое поле =0, а магнитное поле направлено вдоль этой оси. В анизотропной среде из-за индуцированных в ней зарядов векторы и направлены не так, в изотропной среде. На рис. 1 показано распределение плотности д зарядов в плоскости S, ортогональной оси l. Абсолютная величина |д| на этом рисунке тем больше, чем темнее оттенок серого цвета. Здесь же показаны проекции на плоскость S векторов электрического поля, создаваемого этими зарядами. При л?1 и остром угле б между осью l магнитного диполя и осью анизотропии n ортогональная оси l компонента ImEy электрического поля отлична от нуля в точке И на оси l диполя. В ближней зоне при заданных моменте и расстоянии L между магнитным диполем и точкой наблюдения И эта компонента зависит лишь от коэффициента анизотропии л и от угла б.
Этим можно воспользоваться для определения коэффициента анизотропии л. На рис. 2 приведены примеры палеток для определения л при известной величине ImEy. Компонента ImEy на рис. 2 нормирована на абсолютную величину Hп напряженности первичного магнитного поля магнитного диполя в точке И (Hп=M/2рL3). Величина кажущегося удельного электрического сопротивления ск на рис. 2,б определена по осевой компоненте ImHl магнитного поля.
На основе полученных выражений для поля магнитного диполя выполнены численные расчеты различных характеристик компоненты Hl магнитного поля в зависимости от угла между осью l магнитного диполя и осью анизотропии n, длины L зонда, частоты f и др. Расчеты проведены для частот, применяемых в индукционном каротаже (ИК), высокочастотном каротаже проводимости (ВКП), волновом диэлектрическом каротаже (ВДК). Анализ результатов расчетов показал, в частности, что в случае квазистационарного поля при близких к прямому углах , или достаточно больших острых углах влияние n и t на характеристики осевой компоненты Hl магнитного поля на оси l зависит от параметра и от измеряемой характеристики (фаза, разность фаз, отношение амплитуд и др.). Это означает, что при исследовании горизонтальных и наклонных скважин зондами ВКП с соосными генераторными и измерительными магнитными диполями и при соответствующем выборе характеристик зондов в принципе могут быть определены параметры n и t.
Что касается аппаратуры высокочастотного индукционного каротажного изопараметрического зондирования (ВИКИЗ), то (в рамках модели однородной среды) анизотропия по с одинаково сказывается на показаниях всех пяти зондов этой аппаратуры в связи с их изопараметричностью.
Результаты измерений в ВДК зависят от е и среды. Поэтому для определения е требуется измерение как минимум двух разных характеристик компоненты Hl (Даев Д.С.). Расчеты показали, что (при высокой частоте f) в зависимости от угла электрические параметры анизотропной среды (t, t, n, n) сложным образом и по-разному влияют на различные характеристики компоненты Hl. Вследствие этого в анизотропной среде при значительно отличающихся от 0o углах определяемая по этим характеристикам кажущаяся диэлектрическая проницаемость к может иметь аномально- низкие (а в некоторых случаях - отрицательные) значения. Это ограничивает возможности применения ВДК с целью определении истинных диэлектрических проницаемостей t, n анизотропных горных пород.
Раздел 2.3. Приведен вывод выражений для поля элементарной тороидальной антенны, произвольно ориентированной относительно оси анизотропии. Они получены на основе дифференцирования по источнику (стороннему возбудителю поля) полученных в разделе 2.2 формул для векторного потенциала магнитного типа поля магнитного диполя с последующим вычислением компонент электромагнитного поля элементарной тороидальной антенны. На основе анализа полученного решения и результатов численных расчетов установлено, что в однородной анизотропной среде поле тороидальной антенны может испытывать значительное влияние n и .
При возбуждении поля тороидальной антенной T влияние параметров n и на результаты измерений - особенно велико, если в качестве датчика R также применить тороидальную антенну. В работе проанализированы особенности влияния параметров анизотропной среды на показания осевых и экваториальных зондов из тороидальных антенн.
На рис. 3 сопоставлены зависимости кажущихся сопротивлений к для осевых зондов из электрических диполей и осевых зондов из элементарных тороидальных антенн от угла б между осью зонда и осью анизотропии (в ближней зоне). При ориентации дипольно - осевого зонда по оси анизотропии (=0о) в соответствии с парадоксом анизотропии к=t, а для зонда из тороидальных антенн в этом случае значение к превосходит n в л2 раз.
Раздел 2.4. Приведено решение прямой задачи о поле линейной тороидальной антенны (в форме окружности), лежащей в плоскости, ортогональной оси анизотропии n однородной среды. Полученные выражения послужили при решении прямой задачи в разделе 3.5 для модели среды, соответствующей условиям каротажа в процессе бурения.
Раздел 2.5. Приведены решения прямых задач о поле линии AB переменного тока I и бесконечно-длинного кабеля (БДК), ориентированных по оси анизотропии.
Из полученных решений, в частности, следует, что если БДК лежит на оси Z цилиндрических координат, ориентированной по оси анизотропии n, то электрическую составляющую его поля определяют простые выражения: ,
Здесь I=I0exp(?it) - сторонний ток в кабеле, r - расстояние до кабеля, K0 - модифицированная функция Бесселя второго рода нулевого порядка аргумента . Согласно первому выражению в (2), компонента EzБДК поля БДК зависит от волнового числа kn и не зависит от kt. Из второго выражения в (2) следует, что в ближней зоне (|kn|r<<1) компонента то есть эта составляющая поля не зависит от электрических параметров среды и численно равна произведению ?210?7If. При |kn|r<<1 поле близко к однородному. Компоненту квазистационарного поля в ближней зоне определяет выражение: Следовательно, эта составляющая поля не зависит от t, но испытывает небольшое влияние n.
Глава 3. Поле при учете влияния цилиндрических границ (скважины, зоны проникновения, буровой колонны).
В этой главе приведены решения осесимметричных прямых задач электродинамики для моделей анизотропной среды с цилиндрической границей или двумя коаксиальными цилиндрическими границами и для различных сторонних возбудителей поля, а также результаты численных расчетов на основе полученных решений. Такие модели отвечают пройденному скважиной пласту неограниченной мощности. На рис. 4,а,б показаны 2-х слойная и 3-х слойная модели среды (последняя - для пласта с зоной проникновения), а также модель среды для условий каротажа в процессе бурения (рис. 4,в). Приведенные в этой главе решения задач послужили также для задания конечно-разностных аналогов краевых условий и условий у особых точек и линий при 2D- моделировании.
Разделы 3.1, 3.2. Приведены решения прямых задач о поле переменного электрического диполя и поле элементарной тороидальной антенны, расположенных на оси 2-х слойной модели среды (рис. 4,а).
Сравнение полученных решений показало, что в случае, когда эти два сторонних возбудителя поля расположены в изотропном цилиндре (модель скважины), электромагнитное поле, возбуждаемое переменным электрическим диполем, либо элементарной тороидальной антенной, испытывает идентичное влияние параметров окружающей скважину среды. Например, осевую компоненту Ez электрического поля электрического диполя и элементарной тороидальной антенны в скважине определяют выражения:
,
где dl - длина диполя (отрезка переменного тока I), - «момент» тороидальной антенны, k1- волновое число среды в скважине, kt и л - продольное волновое число и коэффициент анизотропии в окружающей скважину среде, =(1/1)?iще0е1 - комплексная удельная электропроводность скважины, a- радиус скважины, I0, I1, K0, K1 - модифицированные функции Бесселя, r, z - цилиндрические координаты точки наблюдения при расположении начала координат в центре стороннего возбудителя поля. ток анизотропный пласт скважина
Таким образом, если тороидальная антенна или электрический диполь находятся в изотропной скважине, то влияние параметров окружающей скважину среды, в частности ее анизотропии, на поле этих сторонних возбудителей - одинаково, что коренным образом отличается от того влияния, которое оказывают их на поле параметры однородной анизотропной среды (рис. 3). То, что расположенная в скважине элементарная тороидальная антенна аналогична электрическому диполю в отношении «чувствительности» поля к параметрам окружающей скважину среды подтвердили также результаты решений прямых задач и численных расчетов для моделей среды, более сложных, чем показанная на рис. 4,а.
Результаты расчета кривых зондирования для лежащего на оси скважины дипольно-осевого зонда (либо осевого зонда из элементарных тороидов) показали, что в присутствие скважины поперечное удельное сопротивление n окружающей скважину среды может оказывать лишь небольшое влияние на результаты измерений (до 20 - 40 % при =2)
Раздел 3.3. Получены выражения для электромагнитного поля, возбуждаемого лежащим на оси 2-х слойной модели среды (рис. 4,а) кабелем с точечным токовым электродом A. Электрический ток по кабелю I=I0?exp(?iщt). На основе полученного решения проведено моделирование, пример результатов которого показан на рис. 5. В представленном примере расчеты выполнены для составляющих осевой компоненты Ez поля в скважине и для частоты f=1кГц, значительно превышающей те частоты, при которых проводят измерения в каротаже КС.
На рис. 5,а,б показаны зависимости от координаты z (при расположении начала цилиндрических координат в электроде A) величины к, определенной по составляющей ReEz, и зависимости от z компоненты ImEz. При небольших расстояниях от точки наблюдения до электрода A кривые к(ReEz) соответствуют 2-х слойным кривым зондирования для предельного градиент-зонда каротажа КС.
Рис. 5,в количественно характеризует степень влияния n на характеристики поля (при л=2). Величина , где или соответствуют результатам расчетов для случаев, когда окружающая скважину среда - изотропна и ее удельное электрическое сопротивление равно сt, либо - анизотропна и ее характеризуют параметры сt и сn=л2•сt. Чем больше величина || (или другие аналогично вычисляемые величины), тем значительнее влияние, оказываемое поперечным удельным сопротивлением сn на характеристику поля. В рассматриваемом примере - это составляющие ReEz и ImEz.
На рис. 5,а,в видно, что величина n оказывает наибольшее влияние на к(ReEz) при тех расстояниях |z|, при которых наблюдаются максимумы на кривых зондирования к(z). Но влияние это - не превышает 20- 30% при л=2. Для модели пласта неограниченной мощности компонента ImEz испытывает лишь небольшое влияние среды (рис. 5,б,в). При достаточно больших расстояниях от электрода A до точки наблюдения, расположенной выше электрода, сторонний возбудитель поля можно рассматривать как БДК. В этом случае в соответствии с (2) ReEz не зависит от среды, а ImEz испытывает небольшое влияние n.
Раздел 3.4. Получены решения прямых задач для показанной на рис. 4,б 3-х слойной модели среды для случаев, когда сторонними возбудителями поля являются переменный электрический диполь, токовая линия AB и кабель с точечным электродом A.
Для случая, когда поле возбуждает кабель с электродом, представлены расчетные зависимости к(ReEz) от расстояния до электрода. При низкой частоте, при которой выполнены расчеты, эти зависимости можно рассматривать, как 3-х слойные кривые зондирования для предельных градиент-зондов каротажа КС. Показано, что наличие зоны проникновения в анизотропных пластах-коллекторах может привести к заметному изменению степени влияния n и t на результаты каротажа КС. При понижающем проникновении (чему при обозначениях на рис. 4,б отвечает неравенство 2t<3t) влияние поперечного удельного сопротивления 2n зоны проникновения проявляется значительно сильнее, чем при повышающем проникновении (2t>3t), что при учете особенностей преломления силовых линий электрического поля на границах и выражения (1) имеет достаточно простое физическое объяснение.
Раздел 3.5. Приведено решение прямой задачи для показанной на рис. 4,в модели среды при возбуждении поля линейной генераторной тороидальной антенной (в форме окружности) и проанализированы результаты численных расчетов, полученные на основе этого решения. Если при каротаже в процессе бурения в качестве датчика также пользоваться линейной тороидальной антенной, соосной с генераторной антенной, то измеряемое напряжение - пропорционально азимутальной компоненте H магнитного поля. При применении двух измерительных антенн могут быть измерены относительные или разностные характеристики компоненты H.
В этом разделе работы рассмотрены сначала результаты моделирования, полученные при удельном электрическом сопротивлении с скважины, типичном для условий, когда промывочной жидкостью является глинистый раствор. Показано, что при соответствующем выборе измеряемой характеристики и размеров зонда влияние поперечного сопротивления n на определяемое значение к может быть значительным (50% и более при л=2).
При использовании в качестве промывочной жидкости растворов на нефтяной основе удельное электрическое сопротивление с имеет высокие значения. Результаты расчетов показали, что при с, превышающем первые тысячи Омметров, и при соответствующем выборе измеряемых характеристик определяемая по ним величина к близка к n, но требования к точности измерений при этом - чрезвычайно высоки.
Глава 4. Поле наклонного магнитного диполя в присутствии анизотропного пласта и кривые для зондов электромагнитного каротажа против пластов.
Раздел 4.1. Представлено решение прямой задачи о поле переменного магнитного диполя с произвольно ориентированным моментом в присутствии анизотропного пласта (плоскопараллельного слоя) с осью анизотропии n по направлению нормали к границам пласта. Вмещающая пласт среда - изотропна.
Раздел 4.2. На основе полученного решения проведено моделирование с целью выяснения вопроса о том, как наклон осей зондов электромагнитного каротажа (ИК, ВКП, ВДК) относительно границ анизотропных и изотропных пластов сказывается на конфигурации кривых профилирования и характерных значениях против пластов.
При анализе результатов моделирования требовалось принять во внимание следующее. Для рассматриваемой модели пласта при направлении оси l зонда (и момента ) под ненулевым углом к оси анизотропии n и нормали к границам пласта (по сравнению со случаем, когда =0о) имеется компонента En электрической составляющей поля магнитного диполя по оси n и по нормали к границам. Поэтому есть и электрические заряды - поверхностные на границах пласта и объемные в пласте, если он анизотропен.
При тех частотах, которые применяют в ИК, кривые против пластов получены не только для осевой компоненты ImHl магнитного поля (и соответствующего ей значения к), но и для составляющих ImHN, ImEy, ортогональных оси зонда (рис. 1). Последняя из них представляет особый интерес, так как в однородной среде в ближней зоне она зависит лишь от одного электрического параметра среды - коэффициента . Данные моделирования показали на каких интервалах кривых против пластов и какие по величине искажения вносит влияние границ пласта (индуцированных на них зарядов) в значения этих компонент поля по сравнению моделью однородной среды.
Результаты расчетов для 3-х и 4-х элементных зондов ВКП (ВИКИЗ, MPR) и зонда ВДК против анизотропных пластов позволили установить некоторые особенности влияния границ анизотропных пластов как на кривые профилирования при измерении характеристик магнитного поля, так и на соответствующие кривые кажущихся значений к и кажущейся диэлектрической проницаемости к для зонда ВДК при различных углах б.
Глава 5. Численное решение осесимметричных прямых задач теории электромагнитного каротажа для анизотропных моделей среды методами конечных разностей.
Получены выражения, послужившие основой для создания алгоритмов численного решения осесимметричных (2D) прямых задач для соответствующих условиям ГИС моделей анизотропной среды при применении различных сторонних возбудителей поля. Алгоритмы разработаны для случаев, когда электрическая составляющая гармонически меняющегося поля имеет четную, а магнитная - нечетную цилиндрическую симметрию относительно оси симметрии Z модели среды. Тогда в цилиндрической системе координат r, , z с осью Z электрическое поле , а магнитное поле , где t - время. В этом случае сторонними возбудителями поля могут служить лежащие на оси Z переменный электрический диполь, линия AB переменного тока, кабель с точечным электродом или БДК, а также - элементарная или линейная (в форме окружности) тороидальные антенны с осью Z.
Сопоставлены результатов численных расчетов для 1D- моделей среды, полученные на основе аналитических решений прямых задач и разработанных алгоритмов конечно-разностных расчетов. Приведена оценка погрешности конечно-разностных расчетов для таких моделей среды.
Раздел 5.1. Представлены дифференциальное уравнение для комплексной амплитуды (r, z) компоненты H и выражения, позволяющие определить комплексные амплитуды , компонент электрического поля по величине . В однородных (по электрическим параметрам) областях вне сторонних токов величина удовлетворяет уравнению . Приведены условия краевой задачи (соответствующие данной модели среды и сторонним возбудителям поля), обеспечивающие единственность решения.
Подобные документы
Процессы разуплотнения горных пород. Электромагнитное поле в моделях разуплотненных структур трещиноватого типа. Зависимость электропроводности горных пород от доли трещин и их заполнения в процессе разуплотнения высокоомным или низкоомным флюидом.
курсовая работа [878,7 K], добавлен 18.04.2015Исторический образ, обзор первобытной обработки камня. Залегания горных пород и их внешний вид. Структура, текстура горных пород Южного Урала. Способы и оборудование для механической обработки природного камня. Физико-механические свойства горных пород.
курсовая работа [66,9 K], добавлен 26.03.2011Классификация горных пород по происхождению. Особенности строения и образования магматических, метаморфических и осадочных горных пород. Процесс диагенеза. Осадочная оболочка Земли. Известняки, доломиты и мергели. Текстура обломочных пород. Глины-пелиты.
презентация [949,2 K], добавлен 13.11.2011Сущность интрузивного магматизма. Формы залегания магматических и близких к ним метасоматических пород. Классификация хемогенных осадочных пород. Понятие о текстуре горных пород, примеры текстур метаморфических пород. Геологическая деятельность рек.
реферат [210,6 K], добавлен 09.04.2012Характеристика структуры, изучение строения и определение размеров пор горных пород. Исследование зависимости проницаемости и пористости горных пород. Расчет факторов проницаемости и методов определения содержания в пористой среде пор различного размера.
курсовая работа [730,4 K], добавлен 11.08.2012Проведение на электронных вычислительных машинах имитационных лабораторных испытаний горных пород и определение их механических свойств (пределов прочности, модуля упругости и коэффициента Пуассона). Теории определения прочности горных пород Кулона-Мора.
курсовая работа [3,8 M], добавлен 27.06.2014Методы определения возраста горных пород, слагающих Землю. Возраст пород слоя Базальт Карденас в восточной части Большого Каньона. Геологическая “блоковая" схема расположения пластов горных пород Большого Каньона. Ошибки радиологического датирования.
реферат [1,4 M], добавлен 03.06.2010Применяемое буровое оборудование и режимные параметры при разрушении горных пород. Характеристика термодинамических параметров зарядов промышленных взрывных веществ. Расчет параметров взрывных работ для рыхления пород при бурении в блоках на карьере.
курсовая работа [494,0 K], добавлен 02.06.2014Исследование особенностей осадочных и метафорических горных пород. Характеристика роли газов в образовании магмы. Изучение химического и минералогического состава магматических горных пород. Описания основных видов и текстур магматических горных пород.
лекция [15,3 K], добавлен 13.10.2013Разработка угольных месторождений. Факторы, влияющие на параметры процесса их сдвижения: вынимаемая мощность пласта, глубина горных разработок и угол падения пород, строение горного массива и физико-механические свойства пород, геологические нарушения.
контрольная работа [65,8 K], добавлен 15.12.2013
