WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

 

 

  На правах рукописи

Кеворкянц Сурен Сергеевич

Теоретические и методологические основы проектирования и интерпретации межскважинного радиопросвечивания  при поисках рудных тел в слоисто-анизотропных средах

Специальность:  25.00.10геофизика, геофизические методы

  поисков полезных ископаемых

А В Т О Р Е Ф Е Р А Т

диссертация на соискание ученой степени

доктора физико-математических наук

ТРОИЦК - 2007

Работа выполнена в  Центре геоэлектромагнитных исследований

Института Физики Земли Российской Академии Наук

Официальные оппоненты: 

  академик РАЕН,  доктор физико-математических наук, профессор  Дмитриев В.И.

  доктор физико-математических наук, профессор  Юдин М.Н.

  доктор физико-математических наук Губатенко В.П.

Ведущая организация:

Московский государственный университет имени М.В. Ломоносова (кафедра

геофизических методов исследований земной коры геологического факультета)

Защита диссертации состоится 15 ноября  2007 г.  на заседании диссертационного Совета Д.212.121.07 при Московском  государственном  геологоразведочном университете (МГГРУ).  Адрес: 117997, Москва, ул. Миклухо-Маклая, 23,  в ауд. 6-38  в 15 часов.

Телефон/Факс: +7(495)4380828

С диссертацией можно ознакомиться в библиотеке Московского государственного

геологоразведочного университета. 

Автореферат разослан  “ ____ “ _____________  2007 г.

Ученый секретарь Диссертационного Совета Г.Н.  Боганик

кандидат технических наук, профессор

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Подземные радиоволновые методы, называемые сокращенно РВМ, на практике начали применяться с 50-х годов [Тархов, 1955; Даев, 1959; Петровский, 1959 и др.], при этом наиболее широкое применение имел теневой метод (радиопросвечивание) [Даев, 1959; Петровский, 1961; Бондаренко,  Коваленко, Тархов, 1962; Грачев, 1966, 1972; Иванов, Петровский, Привезенцев, 1969; Савицкий, 1969; Мамаев 1972; Мамаев, 1978; Борисов, Гуревич, Чигирина, 1983  и др.].

Метод радиопросвечивания (РП) развивался в вариантах межскважинном, шахтно-скважинном (скважина-выработка и т.п.), скважина-поверхность и  выработка-поверхность.

До 70-ых годов в круг задач РВМ, как правило, входили поиски и разведка  контрастных тел (локальных или пластообразных неоднородностей) расположенных в относительно однородных массивах пород, где можно было пренебрегать влиянием земной поверхности (при достаточной глубине исследуемой области) и границ с другими геологическими структурами. Поэтому способы интерпретации РВМ основывались на формальном представлении нормального поля как составляющих полей дипольных источников в однородной среде. Применительно к интерпретации данных РП при поисках и разведке пластообразных залежей (сульфидных жил, угольных месторождений, пластообразных залежей субгоризонтального залегания) развивались способы интерпретации, основанные на квазилучевых или физико-оптических представлениях о распространении радиоволн [А.Д. Петровский, 1967, 1971; Соколов, 1975; В.Н. Мамаев, 1978 и др.]. 

В 70-ые под руководством и при участии А.Д. Петровского началось опробование, а затем и практическое применение межскважинного радиопросвечивания (РП-МС) с целью поисков погребенных кимберлитовых тел в слоистых осадочных толщах нижнего палеозоя и верхнего докембрия районов Западной Якутии [М.С. Бехтерева, 1978; М.С. Бехтерева, В. И. Привезенцев 1978]. Основными поисковыми районами в то время являлись  Далдыно-Алакитский и Малоботуобинский районы.

В процессе исследований были установлены рабочие частоты РП:  312 кГц – основная и 625 кГц – дополнительная (вспомогательная), оптимальные для сетей поискового бурения 500×500 м  в первом из вышеуказанных районов и 350×350 м –  во втором районе, где породы обладают более высоким поглощением. В задачу РП-МС, применяемых при поисках кимберлитовых трубок (КТ), входило выявление их верхней (приторцевой) части в перекрытой траппами осадочной толще мощностью 30÷50 м, сложенной относительно слабо поглощающими (высокоомными) плотными породами. Обнаружение поискового объекта путем применения РП на 30÷50- метровом интервале, составляющем не более 1/51/8 длин радиоволн, при среднем расстоянии между скважинами от 350 до 500 м (в зависимости от поискового района) представляет весьма необычную для традиционных условий применения данного метода и непростую задачу. Заметим, что в настоящее время при использовании более мощной современной аппаратуры (например, передатчика РПД-2С) и рабочей частоты 156 кГц РП проводится при расстояниях между скважинами от 500 до 800 м в зависимости от района работ. Для успешной интерпретации РП, проводимого на коротком глубинном интервале при относительно больших межскважинных расстояниях, требуется непременный и достаточно точный учет влияния выше- и нижележащих пород, граничащих с породами рабочего интервала (лучевой подход здесь не подходил из-за многократных отражений и наличия скользящих волн). В связи с этим, были введены такие понятия, как  радиоволновые геомодели (РВ-ГМ) и  типы РВ-ГМ  вмещающей среды, на основе которых можно было бы строить соответствующую систему интерпретации РП-МС и получать более точные оценки нормального поля, учитывающие особенности распространения радиоволн в том или ином типе РВ-ГМ. С целью более эффективного применения РП-МС и повышения надежности интерпретации результатов радиопросвечивания в контрастно-слоистых толщах в качестве вспомогательного метода был предложен односкважинный радиоволновой метод (РВМ-ОС), выполняемый с метровыми антеннами на более высоких частотах (1250, 2500 кГц), и были привлечены петрофизические исследования [Петровский А.Д., Бехтерева М.С.,  Бондаренко А.Т. и др., 1985; Бондаренко, Ковалев, 1988 и др.]. Результаты петрофизических исследований указали, в частности, на наличие существенной электрической анизотропии [Зинчук, Бондаренко, Гарат, 2002] пород рабочего интервала и граничащих с ними осадочных толщ. Таким образом, был разработан скважинный радиоволновой комплекс РВМ (односкважинные и межскважинные радиоволновые измерения и петрофизические исследования) и была доказана геологическая и экономическая эффективность его применения в комплексе с бурением при поисках кимберлитовых трубок в Западной Якутии [Петровский, Бехтерева, Привезенцев, 1982; Кеворкянц, 1982;  Кеворкянц, Бехтерева, 1984 и др.]

В 1990 г. при использовании разработанной автором методики интерпретации комплекса РВМ его коллегами – Ю.Д. Ковалевым  и А.С. Фроловым на участке Бобкова Якутской кимберлитовой провинции были обнаружены три кимберлитовые трубки «Радиоволновая», «40 лет АМГРЭ» и «ЦНИГРИ».

С внедрением РВМ в комплекс поисковых работ при поисках КТ  Западной Якутии и, в частности с обнаружением КТ, а также с расширением области применения РП (поиск нефти и газа в слоистых осадочных толщах, инженерные задачи и т.д.) [Lytle, Lager et al., 1981;  Johnson, Richardson et al., 1984; Фролов, Истратов и др. 2001; и др.] существенно повысилась и значимость РВМ среди других методов подземной геофизики. С ростом значимости РВМ и степени их вовлечения в практику широкомасштабных геологоразведочных работ (поиски и разведка КТ и месторождений нефти и газа в слоисто-анизотропных толщах), а, следовательно, и в решение проблемы повышения отечественной сырьевой базы существенно возросла и необходимость повышения эффективности РВМ, что невозможно без создания соответствующих теоретико-интерпретационных основ.

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

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

  В соответствии с поставленной целью определились следующие задачи:

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

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

- создание строгих (с контролируемой точностью) и достаточно простых приближенных (в пределах допустимой для практики погрешности) алгоритмов решения прямой задачи РП в основных типах поисковых радиоволновых геомоделей слоисто-анизотропных сред как однородных по простиранию, так и содержащих локальные неоднородности;

2. Уточнение понятия зоны, существенной для распространения радиоволн в поглощающих средах (радиогеоинформационной области – РГИО), и разработка на его основе методологических принципов построения и типизации упрощенных радиофизических моделей (радиоволновых геомоделей – РВ-ГМ) распространения радиоволн, путем выделения их из слоисто-анизотропной толщи, как областей обеспечивающих эффективное применение РП;

3. Разработка алгоритмов и формул для эффективной интерпретации межскважинного РП при поисках контрастных объектов в типовых радиоволновых геомоделях слоисто-анизотропных сред и при изучении пластовых тел – экранов и волноводов на стадии поисков и оценки запасов;

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

6.  Разработка основных принципов проектирования межскважинного радиопросвечивания на основе учета результатов математического моделирования, использования результатов РП-МС на эталонных участках, критериев подобия и количественных оценок диаметра РГИО.

Научную новизну  диссертационной работы определяют такие ее результаты, как

1) вывод и доказательство оптимальности условий калибровки Лоренца (лежащих в основе строгого решения прямых задач) для векторных и скалярных потенциалов электрического и магнитного типов в средах с одноосной электрической и магнитной анизотропией;

2) установление линейной дифференциальной связи между диагональными и недиагональными элементами тензорных функций (аффиноров) Грина электрического и магнитного типов в одномерно-неоднородных (плоскослоистых) одноосно-анизотропных средах с электрической и магнитной анизотропией; эта связь позволяет получить физически наглядные представления элементов указанных тензоров и компонент векторов электрического и магнитного полей через две скалярные функции, каждая из которых является фундаментальным решением одного из двух уравнений типа Гельмгольца, одного, учитывающего электрическую, и другого, учитывающего магнитную анизотропию; 

3) две формы векторных интегральных выражений типа Стрэттона-Чу  для слоистых одноосно-анизотропных сред, содержащих локальную неоднородность, приводящие к различным видам систем  объемных и граничных интегральных уравнений электродинамики;

4) система из двух векторных сингулярных граничных интегральных уравнений 2-го рода,  получаемая из выражений типа Стрэттона-Чу в результате тензорной нормировки  интегральных соотношений, которая при сшивании интегралов по внешней и внутренней сторонам границы неоднородности, обеспечивает взаимное исключение членов подынтегральных выражений, имеющих сингулярность 3-го порядка; эта сингулярность, возникая у вторых производных тензоров Грина при M M0 , приводит к расхождению интегралов ;

5) алгоритмы приближенных (аналитических и полуаналитических) вычислений полей дипольных источников в типовых радиоволновых геомоделях продуктивной толщи слоисто-анизотропного разреза, и основанное на них приближенное решение прямой задачи РП-МС в двух основных типах радиоволновых геомоделей (экранированное полупространство и слой-волновод), содержащих поперечный круговой цилиндр конечной проводимости;

6) совокупность интерпретационных формул, алгоритмов и поправок, обеспечивающих интерпретацию РП-МС в основных типах РВ-ГМ, а также в моделях со значительным влиянием обходных волн;

7) система совместной интерпретации результатов комплекса РВМ односкважинных, межскважинных измерений (на разных частотах) и петрорадиофизических исследований, основанная на учете частотной дисперсии и анизотропии вмещающих пород;

8) новая трактовка понятия зоны, существенной для распространения радиоволн (обобщения первой зоны Френеля для поглощающих сред), названного радиогеоинформационной областью (РГИО), основанная на интегральной формуле Кирхгофа для полей диполей.

Практическая ценность  диссертации определяется тем, что

- полученные в ней теоретический материал может быть использован как при решении теоретических задач в различных методах электроразведки, так и при решении практических задач прикладной электродинамики;

- разработанные диссертантом методика обработки РП-МС в основных поисковых радиоволновых геомоделях слоисто-анизотропной среды и ее модификация для моделей, где является существенным вклад обходных волн, позволяют значительно повысить достоверность результатов обработки интерпретации данных РП-МС;

- предложенная система совместной интерпретации  комплекса РВМ, включающего межскважинные, односкважинные измерения и петрорадиофизические исследования, позволяет существенно (не менее, чем в 2 ÷3 раза) повысить вероятность обнаружения кимберлитовых тел и надежность отбраковки площадей, относимых к категории неперспективных;

- развитая в работе методология выделения и типизации поисковых радиоволновых геомоделей слоистых сред, проектирования РВМ и решения обратной задачи радиопросвечивания была передана в 2003 году в виде методических рекомендаций Амакинской ГРЭ АЛРОСА, ведущей поиск КТ в Западной Якутии; она также может быть использована при поисках и разведке разных видов полезных ископаемых, залегающих в слоисто-анизотропных толщах и обладающих повышенным поглощением по отношению к вмещающим породам.

Апробация и публикации. Основные положения и результаты работы докладывались на XXX  Международном геофизическом симпозиуме, Москва, 23-28 сентября 1985 г., часть 1, Международном семинаре «Вопросы теории и практики геологической интерпретации гравитационных, магнитных и электрических полей», Москва, МГРИ 1-4 февраля 1993 г.,  Международной научной конференции  "Геофизика и современный Мир", 9-13 августа, 1993 г., Годичной сессии ученого совета «НТД-96-ЦНИГРИ», 1-4 апреля 1997.,  Второй Всероссийской школе-семинаре по электромагнитным зондированиям Земли, Москва, 28-30 ноября 2005 г.  и т. д.  Всего по теме диссертации подготовлена 31 публикация в виде статей, научных сообщений и тезисов к докладам на конференциях, в том числе 8 статей в рецензируемых журналах. 

защищаемые положения :

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

2. Полученные в работе системы объемных и граничных сингулярных интегральных и объемных интегро-дифференциальных уравнений второго рода для электромагнитных полей в слоисто-анизотропных средах с локальными включениями предоставляют выбор различных подходов к решению прямых и обратных задач радиоволновых методов в зависимости от  специфики  решаемой  задачи,  класса  изучаемой модели и т.п.

Система граничных сингулярных интегральных уравнений 2-го рода получена на основе способа тензорной нормировки ядер интегралов, представляющего новый подход к конструированию граничных интегральных уравнений с устранением сильной сингулярности в ядрах тензорного типа.

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

интерпретации и проектирования РВМ при поисках рудных тел в слоистых анизотропных средах и разведке рудных месторождений пластового типа – экранов и волноводов;

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

4.  Для поисков локальных экранирующих объектов в слоисто-анизотропных средах в работе развита общая методология проектирования РВМ и системной интерпретации данных односкважинных и межскважинных радиоволновых измерений и петрофизических исследований  с уточнением параметров радиоволновой геомодели для каждой отдельной пары скважин. Применение данной методологии позволяет существенно повысить точность оценок нормального поля, надежность и достоверность выделения экранирующих объектов и отбраковки неперспективных площадей, что продемонстрировано на примерах интерпретации РВМ при поисках кимберлитовых трубок Западной Якутии. 

Структура и объем работ.  Диссертация состоит из введения, 4 глав и заключения, а также списка литературы из 97 наименований. Работа содержит 229  страниц, в том числе 43 рисунка  и  10  таблиц. 

Автор выражает глубокую благодарность  академику РАЕН  д. г.-м.н. А.Д. Петровскому за научное руководство при работе над проблемами проектирования и интерпретации и плодотворные обсуждения по структуре диссертации, д. ф.-м. н. проф. Б.С. Светову за внимательное ознакомление с материалами диссертации, ценные замечания в процессе подготовки рукописи, а также организационную помощь при подготовке ее к защите, д. ф.-м. В.В. Спичаку за ценные замечания, д. ф.-м. Файнбергу Э.Б  и  д. ф.-м. Александрову за полезные советы, внимание и поддержку при апробации работы, коллегам по ЦГЭМИ: В.В. Агееву, О.А. Агеевой, С.Д. Каринскому – за полезные советы и замечания при подготовке и  Ю.И. Куксе за полезные советы и неоценимую помощь при оформлении диссертации.

Автор также глубоко признателен своим коллегам и товарищам по ЦНИГРИ  М.С. 

Бехтеревой, А.Т.Бондаренко и  Ю.Д. Ковалеву, В.Ф. Кухареву,  В.В. Мамаеву, А.С. Фролову, В.Ю. Абрамову и др. за многолетнее плодотворное сотрудничество.

С О Д Е Р Ж А Н И Е Р А Б О Т Ы

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

Глава I.  Основы теории радиоволновых методов при изучении слоисто-анизотропных сред. Глава посвящена вопросу оптимизации представлений электромагнитных полей с гармоническим возбуждением  в плоскослоистой одноосно-анизотропной среде, при этом допускается, что электромагнитные параметры среды могут изменяться вдоль главной оси анизотропии как кусочно-дифференцируемые функции (слоистую среду такого типа, следуя работе [Дмитриев, 1969], условимся называть градиентной). Для общности рассмотрена среда с одноосной анизотропией электрических параметров (удельной электропроводности и диэлектрической проницаемости) и магнитной проницаемости при совмещенных главных осях анизотропии указанных параметров.

Электромагнитные поля в слоистых одноосно анизотропных средах, как известно,  могут быть описаны через электродинамические (или электромагнитные) потенциалы [Тихонов, 1959; Четаев 1970, Ваньян, 1965; Светов, Губатенко1, 1988 и др.] или вектор Герца [Петровский, 1971]. Другим способом описания полей является их представление через поперечные (в направлении главной оси анизотропии) компоненты векторов электрического и магнитного полей [Kong,1972; Таборовский, Эпов, 1977;  Кеворкянц, 1991  и др.]. Выражение векторов поля через электромагнитные потенциалы  представляет собой удобный путь для расчетов полей в слоистых и слоисто-анизотропных средах (с плоско-параллельными границами раздела слоев), традиционно применяемый в электроразведке.

Изучение более сложных сред и, в частности, слоистых (в том числе градиентных) сред, содержащих локальные неоднородности, требует более универсальных способов описания полей, к числу которых можно отнести их представление через электродинамические тензорные функции Грина (ТФГ) [Дмитриев, 1969]. Применение ТФГ позволяет одинаково эффективно, как описывать первичные поля от произвольных локальных источников, так и определять вторичные поля, рассеянные неоднородностями, расположенными в слоистой (в том числе градиентной) среде [Дмитриев, 1969; Захаров, Ильин, 1970]. Фундаментальные тензорные функции Грина слоистой задачи вводятся через векторные потенциалы (точнее их интегральные представления через ТФГ), в силу чего в начале данной главы изучены основные соотношения для векторных потенциалов электрического и магнитного типов (порожденных сторонними, соответственно, электрическими и магнитными токами).

В первом параграфе рассмотрена проблема калибровки векторного и скалярного электродинамических потенциалов  электрического Aэ, э  и магнитного Aм, м типов и постановки граничной задачи для векторных потенциалов в одномерно-неоднородной среде с электрической и магнитной анизотропией. В прямоугольной системе координат x, y, z,  ось z которой совмещена с главной осью анизотропии среды, тензоры комплексной диэлектрической и магнитной проницаемостей данной среды имеют следующее представление:

  , (1)

где ,  n и t – величины удельной электропроводности в направлении главной оси анизотропии и в перпендикулярной к ней плоскости,  – круговая частота  гармонических электромагнитных колебаний.  Показано, что электродинамические потенциалы электрического и магнитного типов, вводимые  традиционным образом [Четаев, 1966], связаны калибровочными соотношениями вида

, ,  (2)

являющимися единственно оптимальными в классе обобщенных условий Лоренца.  Соотношения  (2)  представляют собой обобщения на градиентные среды с электрической и магнитной анизотропией соотношений калибровки, полученных в работах [Тихонов, 1959] и [Четаев, 1962, 1966], соответственно, для потенциалов магнитного и электрического типов в электрически анизотропных средах. Соотношения вида (2) позволяют свести граничные задачу определения векторов электрического и магнитного полей к решению граничных задач для компонент векторных потенциалов или векторов Герца, для решения которых существуют различные подходы. В работе [Петровский, 1971] посредством введения векторов Герца решена  задача о гармоническом возбуждении элементарных поперечного (направленного вдоль главной оси анизотропии) и продольного электрических диполей в однородном электрически анизотропном пространстве при учете токов смещения. Для определения компоненты вектора Герца, описывающего поле продольного диполя (первый нижний индекс  x  указывает направление оси диполя), использована линейно-дифференциальная связь между компонентами и и , полученная автором упомянутой работы из теоремы взаимности.

В работе [Ваньян, 1965] на основе результатов [Четаев, 1962], предложен удобный алгоритм расчета электромагнитных полей дипольных источников и, в частности, определения компоненты поля горизонтального (направленного вдоль оси x) диполя на поверхности слоисто-анизотропного полупространства, параметры которого  –  кусочно-постоянные функции глубины. Решение получено в квазистационарном приближении.

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

Во втором параграфе рассмотрена задача определения тензорных функций Грина (ТФГ) электрического и магнитного типов в слоисто-анизотропной среде.

Общая постановка граничных задач для скалярных компонент тензорных функций Грина слоистой (градиентной) изотропной среды и универсальный алгоритм их расчета впервые даны в работе [Дмитриев, 1969]. Для заданных в ограниченном объеме V плотностей сторонних электрических или магнитных токов векторы-потенциалы связаны с  ТФГ соотношениями следующего вида:

, (3)

где M, M0∈V – точки пространства с координатами, соответственно, x, y, x и x0, y0, x0; и – тензорные функции Грина соответственно электрического и магнитного типов, имеющие в системе координат x, y, z  следующий вид:

  .

Из уравнений, граничных условий и соотношений, приведенных в первом параграфе  для скалярных компонент векторов-потенциалов по аналогии с работой [Дмитриев, 1969], с учетом (3) получены следующие соотношения для диагональных элементов ТФГ слоисто-анизотропной среды:

    ,  (4)

где и представляют собой фундаментальные решения следующих задач:

  , ,  (5)

, , (6)

, - операторы класса Гельмгольца электрического и магнитного типов для градиентной анизотропной среды, переходящие в однородных изотропных областях  в классические операторы Гельмгольца [Кеворкянц, 1987]. Для функций , автором получены следующие соотношения линейной дифференциальной связи с фундаментальными функциями [Кеворкянц, 2000].:

  ,  .   (7)  

Соотношения (4) и (7) связывают  недиагональные элементы ТФГ  с диагональными.

Функции , ,  и ,   выражаются через преобразования Ханкеля вида

(8)

где  r =, из которых следует, что задача определения  всех элементов  ТФГ  обоих типов (электрического и  магнитного ) сводится к определению пары функций  –  ,  , представляющих решение задачи вида

    (9) 

где  ;  , .

В третьем параграфе для решения задачи (9) автор отдает предпочтение известному алгоритму [Дмитриев, 1969] и приводит достаточно детальное его описание для случая  слоисто-анизотропной модели среды. Суть этого алгоритм заключается в переходе от функций и  к  спектральным импедансам, решении для последних уравнения Риккати, затем определении с их помощью функций и  для различных вариантов расположения точек M0 и M относительно друг друга и границ раздела слоистой пачки с верхним и нижним полупространствами. 

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

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

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

Одним из наиболее эффективных и хорошо разработанных методов расчета  электромагнитных полей в слоистых средах, содержащих локальные неоднородности, является метод интегральных уравнений (ИУ), развитый в работах [Дмитриев, 1969; Захаров Е.В., Ильин И.В., 1970; Дмитриев, Захаров, 1987;  Васильев, 1977, Светов, 1984 и др.].  Как один из перспективных методов решения объемных ИУ, в работе [Зингер, Фейнберг 1995] был предложен итеративно-диссипативный метод, получивший дальнейшее развитие и численную реализацию в работах  [Панкратов, Авдеев, Кувшинов 1995;  Avdeev, Kuvshinov, Pankratov, Newman, 1997; Авдеев 2002  и  др.]. Метод ИУ обладает такими достоинствами, как экономичность по времени реализации решения на ПК и возможность получения приближенных аналитических решений таких, например, как приближения Кирхгофа в высокочастотной (коротковолновой) и Рытова-Борна в длинноволновой области. Ценность аналитических приближений решений ИУ заключается в возможности их применения как для достаточно точных оперативных оценок возможностей, так и при решении конкретных задач поисков, разведки и интерпретации полевых материалов РВМ. Один из способов построения систем ИУ – это использование векторных интегральных соотношений типа Стрэттона-Чу,  которые для изотропных моделей однородной области и слоистой среды с неоднородностью приведены соответственно в работах [Вычислительные методы .., 1978; Жданов, Спичак, 1983;  Светов, 1984 и др.] и [ Захаров, Ильин, 1970; Дмитриев, Захаров, 1987  и др.].

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

В первом параграфе данной главы автором на основе подходов, развитых в работах [Дмитриев, 1969; Захаров и Ильин, 1970; Дмитриев, Захаров, 1987], получены две формы векторных интегральных соотношений типа Стрэттона-ЧУ в слоисто-анизотропных средах, в первой из которых ядра интегралов представлены функциями Грина , и их первыми тензорными и векторными производными, во второй  – функциями Грина и их тензорными производными до второго порядка.

Во втором параграфе из первой формы представлений типа Стрэттона-Чу получены система сингулярных граничных ИУ первого рода и система объемных интегро-дифференциальных уравнений второго рода, а из  второй – системы объемных сингулярных ИУ второго рода [Кеворкянц, 1992, 2007].

В третьем параграфе из второй формы представлений (типа Стрэттона-Чу) для векторов поля E и H через поверхностные интегралы получена система граничных ИУ второго рода для слоисто-анизотропных моделей вмещающей среды [Кеворкянц, 1995, 2007].  Представления  первого из упомянутых векторов через граничные интегралы имеют вид

  . (10)

где вектор напряженности первичного электрического поля в области ;     функция Грина области  ,  , R  евклидово расстояние между точками M  и  M0, единичный тензор, ;

, , соответственно волновое число, магнитная и комплексная диэлектрическая проницаемости области .  Вторые (тензорные) производные  функций Грина в ядрах интегралов (12) при стремлении точки M к M0 имеют сингулярность порядка [1/R3(M, M0)]R→0,  что для M0 ∈ S приводит к расхождению интегралов.

В скалярной и векторной задачах рассеяния волн на однородном изотропном теле, помещенном в однородную изотропную среду [Mller, 1957; Купрадзе, 1967], сингулярность порядка O(1/R3)R → 0  ( кубическая сингулярность) в интегралах типа (12) устраняется путем взаимного исключения при сложении векторных произведений первого и второго интегралов вида (10) на векторы нормалей к поверхности S, соответственно, с внешней и внутренней сторон. Описанный выше способ устранения кубической сингулярности не усложняется и при переходе к моделям слоистых изотропных сред, где локальная неоднородность расположена в одном из слоев [Захаров, Ильин 1970]. Это объясняется тем, что тензор Грина слоистой изотропной среды вблизи точки своей сингулярности M=M0 совпадает со скалярной функцией Грина однородной среды с параметрами слоя, в котором находится точка M0.

При переходе к анизотропным средам задача устранения кубической сингулярности  усложняется тем, что в таких средах функции Грина остаются тензорами и в точке своей сингулярности. Один из способов ухода от сингулярности в ядрах интегралов при построении поверхностных ИУ второго рода предложен в работе [Табаровский, Дубман,1989]. Он заключается в замене граничных интегралов по внешней и внутренней сторонам границы неоднородности S интегралами по вспомогательным (фиктивным) замкнутым поверхностям, одна из которых охватывает неоднородность снаружи, а другая заключена внутри нее.  Полученная таким образом система ИУ при дискретизации при дискретизации ИУ образует систему линейных алгебраических уравнений (СЛАУ) с размерностью приблизительно вдвое большей, чем у СЛАУ, получаемой из системы граничных ИУ с интегралами по S.  Кроме того, надо полагать, результат применения данного способа решения в определенной мере будет зависеть от выбора вспомогательных поверхностей. 

В диссертации развивается способ устранения кубической сингулярности на основе специальной тензорной нормировки второго из интегралов (10). Показана возможность получения нормирующих тензоров для моделей, где параметры вмещающей среды представляют собой кусочно-дифференцируемые функции переменной z! , , которыена интервале оси z, где расположена неоднородность, удовлетворяют вместе со своими первыми производными условию Липшица, то есть!

   

Если положить, что   и  - нормирующие тензоры, ,   -

неизвестные векторные функции, , -

-  диагональные тензоры, то система поверхностных (граничных) ИУ 2-го рода для слоисто-анизотропной среды, полученная в результате тензорной нормировки имеет вид 

где

  ,  ,

, - планарные тензоры в плоскости касательной к поверхности S в точке M0, первый из которых единичный, а второй - тензор поворота  на  -90 . Нормирующие тензоры определяются из следующего предельного соотношения:

.,

где  – проекция тензора на плоскость, касательную к  S  в точке M0 ;  ,  – тензорные функции Грина для однородного анизотропного пространства с параметрами  ,  .

 

  Из содержания главы II вытекают следующие выводы:

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

Для получения системы векторных граничных интегральных уравнений второго рода с двумя неизвестными предложена тензорная нормировка интегральных выражений типа Стрэттона-Чу, обеспечивающая исключение сингулярности порядка O(1/R3)R → 0 в поверхностных (граничных) интегралах и описан способ определения нормирующих тензоров.

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

Система граничных ИУ 2-го рода позволяет уменьшить размерность задачи, и, следовательно, существенно сократить объем вычислений, необходимых для ее решения. Эта система представляется более предпочтительной в тех случаях, когда линейные  размеры неоднородности сопоставимы или существенно больше длины волны во вмещающей среде. В таких случаях для получения приближенных решений возможна скаляризация системы с понижением числа неизвестных до одной тангенциальной составляющей вектора E или H, благодаря использованию приближенных граничных условий Леонтовича или более универсальных условий Петровского-Фейнберга  [Басс, 1960]. 

Глава III.  Численные  и  аналитические  решения  прямых  задач радиопросве­чивания. Задачей данной главы является получение алгоритмов и формул для строгих (с контролируемой точностью) и приближенных расчетов электромагнитных полей диполь­ных источников в слоисто-анизотропных моделях, характерных  для  геолого- радиофизических условий применения  РВМ.  Традиционные методики интерпретации радиопросвечивания в слоистых средах, как правило, опираются на лучевые и физико-оптические представления о распространении радиоволн [Петровский, 1967, 1971; Мамаев, 1978]. На самом деле модели распространения полей в РВМ сочетают в себе законы геометрической и физической оптики и явления дифракции и индукции, учет которых позволит существенно повысить надежность интерпретации и достоверность ее результатов. Создание  алгоритмов и приближенных формул, обеспечивающих удовлетворительную точность (в пределах погрешности измерений) в рамках характерных для РВМ типов  радиофизических моделей, позволяет разрабатывать на их основе эффективные методики интерпретации, объединяемые в эффективную систему интерпретации РВМ.

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

Можно выделить два класса слоистых сред: первый класс составляют среды, являющиеся вмещающими для поисковых объектов (чаще это слоистые осадочные толщи), и второй класс – это среды, в которых объектом исследования (поиска, оценки запасов и т.п.) являются протяженные пластовые тела или образования (жилы, дайки, угольные пласты).

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

Радиогеоинформационная область (РГИО), называемая иначе областью, существенной для распространения радиоволн, представляет собой обобщение понятия первой зоны Френеля [Фейнберг 1961, Черный 1978 и др.] на поглощающие среды [Петровский, 1971].  Автором диссертации предложен алгоритм определения диаметра РГИО на основе использования интегральной формулы Гюйгенса-Френеля [Кеворкянц,  Коновалова, 2000] и показано, что РГИО также как и первая зона Френеля имеет форму эллипсоида вращения, фокусы которого совмещены  с точками  излучения и наблюдения.. Важным отличием РГИО от первой зоны Френеля является зависимость площади ее поперечного сечения от задаваемого коэффициента информативности Kи (практический интерес представляют  значения Kи  от 0,9  до  0,99).  В диссертации  предложены принципы классификации (на основе учета границ РГИО) и установлены основные типы РВ-ГМ слоистых сред, благоприятных для применения  в них РП с целью обнаружения локальных экранирующих объектов (базовые РВ-ГМ), на основе которых строится интерпретационная система РВМ. В слоисто-анизотропной толще, состоящей из чередующихся слоев с разными значениями удельного электрического сопротивления (УЭС), поисковый объект, представляющий собой хорошо проводящую локальную структуру (тело, блок пород и т.п.) можно обнаружить только на интервале, или в слое, где вмещающая среда  имеет относительно высокое УЭС (низкий коэффициент поглощения). В зависимости от положения границ РГИО по отношению к границам указанного интервала можно выделить следующие три основных типа РВ-ГМ:

квазиоднородная модель (условно называемая 0-ым типом РВ-ГМ) имеющая место, когда границы эллипсоида вращения, представляющего РГИО не пересекаются ни с одной из границ толщи пород, представляющей благоприятный для  РП  интервал;

модель «экранированное полупространство» (1-ый тип РВ-ГМ) в которой только одна из границ, верхняя или нижняя граница РГИО пересекается, соответственно, с верхней или нижней границей благоприятного для РП интервала;

модель «слой-волновод» имеющая место, когда верхняя и нижняя граница РГИО пересекаются, соответственно, с верхней и нижней границами благоприятного для РП интервала.

Квазиоднородная анизотропная модель вмещающей среды, для которой теория и методика интерпретации РВМ хорошо разработаны [Петровский, 1971;  Руководство, 1977; Борисов Б.Ф. и др., 1984 и т. д.], в диссертации не рассматривается. Ниже приводится пример георазреза  с выделенными на нем интервалами, представляющими базовые РВ-ГМ 1-го и 2-го типов (рис.1). Благоприятный для РП интервал или часть этого интервала исследуемые методом радиопросвечивания, называют рабочим интервалом (РИ) РП.

На стратиграфо-петрорадиофизической колонке поискового района (рис.1) выделены два интервала, в которых можно применять РВМ для поисков кимберлитовых трубок. Верхний интервал представляет РВ-ГМ 2-го типа «слой-волновод» (известняки нижнего силура мощностью до 50 м). Он может служить рабочим интервалом РП, если его мощность не меньше 30-50 м.  Нижний интервал, точнее его верхняя часть (плотные доломиты нижнего ордовика, перекрытые экранирующей толщей глинистых известняков мощностью до 30 м), представляет РВ-ГМ 1-го типа  «экранированное полупространство» и служит рабочим интервалом РП на площадях, где силурийская толща полностью или частично размыта.

Говоря о слоистых моделях второго класса, где определенная пластовая структура является объектом поиска или разведки, можно выделить также два основных типа базовых РВ-ГМ. Первый тип – слой-экран в однородной среде (золото-сульфидные, медно-колчеданные жилы, зоны трещинноватости, обводнения и т.п.) и второй тип – слой-волновод в однородной среде (золото-кварцевые жилы, угольные пласты и т.п.).

Для обоих классов РВ-ГМ, описанных выше, решение прямой задачи сводится к изучению и расчетам полей электрических диполей в общей модели – анизотропный слой (экран или волновод) в однородной анизотропной среде (изотропная среда рассматривается, как частный случай анизотропной). Это решение, в зависимости от положения и ориентации диполей относительно границ слоя, мощности слоя и соотношения электрических параметров слоя и граничащих с ним полупространств, приводит к различным приближенным асимптотическим (аналитическим и полуаналитическим) выражениям, которые составляют основу количественной интерпретации РП в перечисленных выше типах РВ-ГМ.

Второй параграф главы посвящен вопросам численного преобразования Ханкеля для расчетов компонент полей диполей и элементов трехмерной тензорной функции Грина в слоисто-анизотропных средах, выраженных через интегралы вида 

  (l = 0, 1, 2;  s = 0, 1), (11)

– одна из функций вида , , ,  – решения задачи (9), сформулированной в первой главе, , – символы производных по z, z0 . Функция при   стремится к 0, если и к константе, если  . Для эффективного вычисления интегралов вида (11), когда из    следует  (при ),  он представляется в следующем виде:

  (12)

где l = 1, 2,  , , – коэффициенты следующего разложения 

(λ → ∞),

определяемые по трем значениям для заданных трех достаточно больших значениях параметра интегрирования . В правой части (12) второе слагаемое состоит из табличных интегралов, а множитель при функции Бесселя под интегралом, который обозначим , стремится к 0 при  ., что позволяет применить к данному интегралу метод численного интегрирования на интервале , где Q – достаточно большая величина. Для этого указанный интервал разбивается на N подинтервалов, на каждом из которых функция заменяется интерполяционным полиномом степени M, в котором узлами интерполяции являются корни многочленов Чебышева степени M+1, после чего задача сводится к вычислению интегралов вида

, , 

для которых при при  m = 0, 1 имеют место аналитические представления, а для  m >1  - следующие рекуррентные:

  (13) 

Описанный выше способ представляет собой вариант реализации метода Филона применительно к интегралам, содержащим функции Бесселя, предложенный в работах [Дмитриев, 1965], где приведено второе из рекуррентных соотношений (13), и [Ваньян, 1965], где функция типа на каждом подинтервале приближается квадратной функцией. Предложенный выше выбор узлов интерполяции позволяет использовать полиномы достаточно высокой степени (M = 8÷12), увеличив длину и сократив количество подинтервалов интегрирования, на концах которых вычисляются функции Бесселя и интегральная функция Бесселя. По скорости счета он уступает быстрому преобразованию Ханкеля с применением фильтров [Андерсон 1978,  Sourenson 1977,  Nielson 2003], который, однако, резко теряет точность при расчетах слабых полей, когда величина интеграла по модулю меньше 10-6  [Кеворкянц, Виноградова, 1997]. 

  В третьем параграфе  описывается асимптотический способ вычисления интегралов вида (11) и приводятся полученные на его основе приближенные аналитические и полуаналитические решения. Асимптотические методы расчетов полей в слоистых средах занимают важное место в теории акустики и радиофизики [Бреховских, 1973], [Макаров, Новиков, Рыбачек, 1991], [Фелсен, Маркувиц, 1978], в частности, в задачах подземной радиосвязи и радиопросвечивания [Кухарев,1978; Кеворкянц 1979, 1981, 1983,  2004, 2005]. Современные вычислительные средства и существующие алгоритмы численного счета позволяют достаточно быстро и точно рассчитывать интегралы вида (16). Аналитические приближения на основе метода наискорейшего спуска (МНС) достаточно просты и физически интерпретируемы, что делает их незаменимыми для решения задач интерпретации РП и физического понимания особенностей распространения радиоволн в различных типах РВ-ГМ. Автором получены оригинальные алгоритмы реализации МНС и определены области применимости, а также ограничения различных аналитических приближений для полей дипольных источников в радиоволновых геомоделях слоисто-анизотропных сред. При реализации МНС (за исключением случаев, когда луч, соединяющий источник и точку наблюдения, направлен близко к нормали по отношению к границам слоя),  интегралы типа (11) приводятся к виду 

    (l = 0, 1, 2;  s = 0, 1),  (14)

где - функция Ханкеля первого рода, s-го порядка  от аргумента ρ . Для вычисления интегралов вида (14) путь интегрирования деформируется – отводится в верхнюю область комплексной плоскости переменной интегрирования λ,  после чего задача расчета (14) сводится к определению его дискретного и сплошного спектров. Дискретный спектр представляет собой сумму вычетов в полюсах функции , пересеченных при деформации контура интегрирования, сплошной спектр – один или два интеграла по линиям наискорейшего убывания (спуска) подынтегральной функции, огибающим точки ветвления  функции на комплексной плоскости переменной λ.

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

Для интерпретации РП в основных типах РВ-ГМ особое значение имеет изучение полей электрических диполей в частном типе трехслойной анизотропной модели – слоя конечной мощности, ограниченного полупространствами с одинаковыми параметрами (модель «слой в однородной среде»). Анализ и решение задачи расчетов полей диполей в данном типе моделей при различных соотношениях электрических параметров анизотропных слоя и граничащей с ним среды позволили сделать следующие выводы относительно структуры поля для основных частот РП и параметров базовых РВ-ГМ: 

- из счетного числа полюсов функции только вычет в первом полюсе (первая мода)  дает основной вклад в поле дискретного спектра;

-  для модели, в которой слой имеет повышенное УЭС по отношению к вмещающей среде (слой-волновод), практически значимым для поля вблизи слоя является вклад дискретного спектра, а точнее его первой моды; случаи, когда это утверждение не справедливо (очень тонкий и слабоконтрастный слой) особого интереса для РВМ не представляют; 

- для модели, в которой слой имеет пониженное УЭС  по отношению к вмещающей среде (слой-экран),  поле представлено сплошным спектром;

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

На основе изложенных выше выводов получены аналитические представления для компонент поля в базовых радиоволновых геомоделях, наиболее важные из которых приведем ниже, записав компоненты Ez и Ex полей электрических диполей, направленных, соответственно по осям  z  и  x,  в виде

,  .

где  I – ток в диполе, lz  и lx  - длины диполей. Пусть h – мощность слоя-экрана, а плоскости z = z* и z = z*- h  – его границы. Для ez  в РВ-ГМ экранированное полупространство  имеют место следующие представления: . а) ζ = z + z0 -2z* ;  б) ζ = 2(z*-h) - z - z0 ; 

  (15)

,  ; , ;  ;

  -  функция ослабления;

  - величины, зависящие от ζ,  электрических и радиоволновых параметров модели. 

Для  РВ-ГМ слой-волновод приведем компоненту  ez  внутри слоя (z*-h < z, z0 < z*)

  , (16)

где – первый полюс функции из выражения (11) для интеграла , являющийся решением трансцендентного уравнения  =0, которое для данной модели приводится к достаточно простому виду.

При использовании РП для определения мощности и удельной электропроводности крутопадающих экранирующих пластов направления лучей, соединяющих источник и точку наблюдения, близки к нормальному (под углами 70110), а генераторной и приемной антенн – к параллельному границе пласта. Применительно к такой геометрии наблюдений в диссертации приведены первые приближения (в элементарных функциях) для коэффициентов отражения и прозрачности слоя-экрана. Используя соотношения  , || ≈1/ Э, где  - нормальное поле (в отсутствие пласта),  Э  -  коэффициент экранирования, можно по результатам измерения Э оценивать мощность и удельную электропроводность пласта, жилы и т.п. [Соколов, 1975].

На рис. 2 приведены кривые изменения при синхронном перемещении источника и точки наблюдения в моделях слой-экран и слой-волновод в однородной среде в направлении, перпендикулярном к границам слоя. Кривые рассчитаны по приведенным выше приближенным формулам и методом численного интегрирования НТВЦ МГУ (В.И. Дмитриев, Э.А. Федорова, 1979). Сходимость в пределах требуемой точности обоих результатов внутри и вблизи слоя подтверждает пригодность полученных автором приближений для получения достаточно точных оценок нормального поля, интерпретации и проектирования РВМ в слоисто-анизотропных средах, где ведутся поиски и разведка локальных объектов и структур, а также при поисках и разведке пластовых тел.

Автором был выполнен большой объем расчетов поля в РВ-ГМ «экранированное полупространство» и «слой-волновод» для рабочих частот РП, применяемых при поисках кимберлитовых трубок Якутии  (0, 312 кГц – 1,25 МГц), и  достаточно  широкого  диапазона

изменения  параметров  моделей (УЭС, диэлектрической проницаемости  и  коэффициентов

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

приведены результаты расчетов полей в РВ-ГМ экранированное полупространство  и  слой-волновод.

В четвертом параграфе главы приводится приближенное решение задачи рассеяния поля вертикального электрического диполя на поперечном круговом цилиндре повышенной электропроводности, расположенном в одной из базовых радиоволновых геомоделях 1-го (экранированное полупространство) и 2-го (слой-волновод)  типов, характерных для районов поисков кимберлитовых трубок Западной Якутии. Кимберлит имеет слабую контрастность по отношению к толщам песчано-глинистых пород, граничащим непосредственно с породами рабочего интервала (с полупространством в РВ-ГМ типа 1 и слоем-волноводом  в РВ-ГМ типа 2). Поэтому в РВ-ГМ 1-го типа рассматривается полубесконечный цилиндр, примыкающий торцом к внутренней поверхности полупространства (рис. 4a), а в РВ-ГМ типа 2 – цилиндр конечной длины, примыкающими торцами к внутренним поверхностям стенок волновода (рис. 4b). Поскольку поле, прошедшее в экранирующие толщи

Рис. 2. Сравнение асимптотических и численных расчетов поля вертикального электрического диполя в присутствии плоского горизонтального слоя в однородной среде (a слой-волновод, b слой-экран) при перемещении точек излучения (x0, y0, z0) и наблюдения (x, y, z) в направлении оси Z (z=z0):

f - частота поля, r – расстояние между точками излучения и наблюдения, h – мощность слоя; 2, 2 и 1, 1 - удельное электрическое сопротивление и диэлектрическая проницаемость соответственно слоя и среды, в которой он расположен; 1 – результаты расчетов первой моды поля дискретного спектра (a) и поля сплошного спектра (b), выполненных автором; 2 – результаты строгих расчетов полного поля, выполненных одним из способов численного преобразования Ханкеля, в НИВЦ МГУ (В.И. Дмитриев, Э.А. Федорова).

Рис. 3. Типовые теоретические кривые межскважинного радиопросвечивания с синхронным перемещением передатчика и приемника в радиоволновых геомоделях экранированное полупространство  (a)  и «слой-волновод» (b ):

H – глубина относительно подошвы слоя-экрана;  r   расстояние между излучателем и точкой наблюдения; h – мощность слоя-волновода; Z/h – относительная глубина расположения излучателя и точки наблюдения по отношению к середине слоя-волновода.

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

для компоненты Ez, поля в рабочем интервале РП (вне цилиндра), измеряемой при межскважинном РП.

  Поле в присутствии вертикального  цилиндра в РВ-ГМ типа 1:

  ,

где ,  ,  ,

.

Поле в присутствии вертикального цилиндра  в  РВ-ГМ  типа 2:

 

  ; ,  ,

kв –  эффективное  волновое  число  слоя-волновода,

  .

–  спектральные функции, определяемые из граничных условий. Индекс «e» обозначает поле во внешней по отношению к цилиндру области полупространства (рис. 4a) или слоя-волновода (рис. 4b).

  Из полученных в данной главе результатов можно заключить.

Предложены новая физико-математическая трактовка и способ определения (с заданной точностью) границ зоны, существенной для распространения радиоволн (обобщения  первой зоны Френеля в поглощающей среде), названной  радиогеоинформационной областью (РГИО). Сформулированы принципы выделения на основе  учета РГИО интервалов эффективного применения РП в слоисто-анизотропном разрезе, которые по характеру распространения радиоволн сводятся к трем основным типам радиоволновых геомоделей (РВ-ГМ) анизотропных сред, вмещающих поисковые объекты:  РВ-ГМ типа 1 – экранированное полупространство, РВ-ГМ типа 2 – слой-волновод и РВ-ГМ – квазиоднородная среда.

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

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

Глава IV. Интерпретационная система РВМ  при поисках  экранирующих объектов в поисковых горизонтах слоисто-анизотропных сред. При поисках рудных объектов в слоистых средах методом межскважинного радиоволнового просвечивания (РП-МС) одной из основных задач интерпретации является достоверная оценка нормального поля в поисковом горизонте (рабочем интервале РП). Для решения этой задачи в условиях влияния выше и нижележащих толщ и геопомех метод РП-МС применяется в комплексе с односкважинным радиоволновым  методом  (РВМ-ОС), включающем односкважинное радиоволновое профилирование вдоль всей скважины и односкважинное радиоволновое зондирование (РВЗ-ОС) в пределах РИ [Петровский 1968, Бехтерева 1978]. Для детальности исследования околоскважинного пространства односкважинные радиоволновые измерения проводятся с короткими антеннами (длина плеча 1-1,5 м ), эффективность излучения которых обеспечивается при частотах в 4÷8 раз превышающих рабочие частоты. РП-МС. В комплекс РВМ также включены петрорадиофизические (ПРФ) исследования, результатами которых являются ПРФ-характеристики горных пород (УЭС, диэлектрическая проницаемость, коэффициенты поглощения, а также коэффициенты анизотропии и частотные характеристики перечисленных параметров), образующих поисковый горизонт и прилегающие к нему слои [Петровский и др.,1985; Бондаренко, Ковалев, 1988]. Вышеописанный измерительный комплекс РП-МС, РВМ-ОС и ПРФ-исследований назван скважинным комплексом РВМ  для поиска рудных объектов в слоисто-анизотропных средах.

Интерпретационная система комплекса РВМ (ИС-РВМ) представляет собой совместные, взаимно дополняющие и взаимосогласованные анализ и обработку данных односкважинных, межскважинных радиоволновых измерений и петрорадиофизических исследований, конечным результатом которых является заключение о степени перспективности поискового горизонта в пределах изученной площади [Петровский и др.,1985;. Кеворкянц, Бондаренко, Абрамов, 2003; Кеворкянц, Абрамов, Ковалев, 2005].

Теоретико-интерпретационные основы ИС-РВМ разработаны на основе применения комплекса РВМ при поисках кимберлитовых трубок в районах Якутской кимберлитовой провинции (Западная Якутия), где РП проводилось на частотах (0,312, 0,625, 1,25 МГц) при расстояниях между скважинами с передающей и приемной антеннами (не более 500÷560 м). Изучаемые при этом модели сред физически достаточно точно подходили под одни из двух основных типов РВ-ГМ – «экранированное полупространство» и «слой-волновод», охарактеризованных в разделе 3.1. Поэтому эти модели названы чистымии моделями, а разработанная на основе их изучения интерпретационная система названа базовой интерпретационной системой РВМ. 

Начиная с 90-ых годов поиски кимберлитовых трубок в Якутии (Амакинская ГРЭ) методом РП-МС проводились при рабочей частоте РП –  156 кГц и расстояниях между скважинами от 500 до 800 и более метров, что привело проблеме, связанной с обходными волнами. Проходя через слой-экран, граничащий с породами РИ, обходные волны распространяются по вышележащей толще траппов (интрузивных пород, обладающих наиболее высокими УЭС в верхней части разреза) и поэтому не могут нести в себе достаточной информации о наличии экранирующих объектов в поисковом горизонте. Модели собходными волнами являются неблагоприятными, если доля обходных волн в них соизмерима с полезным сигналом или превышает его, и не  очень благоприятными, или нечеткими РВ-ГМ, если она составляет от 10% до 30% от измеряемого поля.

В разделе 4.5 приведены критерии и рекомендации, следование которым дает возможность существенного снижения фактора «нечеткости» модели и перехода к более благоприятным условиям поисков. Для повышения эффективности интерпретации в условиях нечетких моделей в разделе 4.4 автором предложена модифицированная методика обработки данных РП-МС (для нечетких моделей), учитывающая вклад обходных волн. Модификация обработки РП-МС является лишь небольшим дополнением к базовой интерпретационной системе комплекса РВМ, для которой далее термин «базовая» опустим.

.Интерпретационная система комплекса РВМ, включает три основных этапа.

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

Подробное описание первого этапа интерпретации приводится в первом параграфе.  Основные задачи этого этапа: 1) типовая  ПРФ-идентификация модели рабочего интервала с прилегающими к нему толщами; 2) отбраковка результатов, не интерпретируемых с точки зрения количественной интерпретации; 3) выявление визуальных признаков наличия крупных экранирующих объектов в поисковом горизонте межскважинного пространства; 4) оценочная ПРФ-идентификация моделей РИ – оценка ПРФ характеристик моделей рабочего интервала (на основе данных ПРФ-исследований), глубинной отметки верхней границы пород РИ и мощности перекрывающего их экранирующего слоя для каждой скважины, пары скважин, каждого участка, поискового полигона. Первый этап интерпретации проводится путем сопоставления и анализа кривых односкважинного, межскважинного профилирования и данных бурения, а также по оценкам результатов петрорадиофизических исследований. Важной частью его результата является определение типа РВ-ГМ и оценка ее геометрических параметров (верхней границы пород РИ, мощности  слоя-волновода в РВ-ГМ типа 2 и т. п.) и средних по полигону (или по его отдельным участкам) ПРФ-характеристик (ПРФ-идентификация). Полезную роль при ПРФ-идентификации  разреза по данным РП-МС может играть моделирование теоретических кривых поглубинного изменения измеряемой амплитуды поля Ez  при синхронном перемещении передающего и приемного устройств. Подбор теоретической кривой, наиболее близкой к экспериментальной, позволяет получить приближенную оценку электрических параметров физико-геологической модели поискового горизонта для данной пары скважин. Примеры сопоставления теоретических и экспериментальных кривых РП-МС  приведены на рис. 5 a, b. 

ПРФ-идентификация РВ-ГМ с погрешностью до ±30-35%  позволяет решать прямую задачу РП-МС и РВЗ-ОС в базовых типах РВ-ГМ, и по результатам ее решения  получать в процессе обработки более достоверные оценки коэффициентов поглощения и нормального поля в породах РИ каждой пары скважин. Оценочной  ПРФ-идентификация пород РИ и прилегающих к ним слоев имеет большое значение для повышения эффективности совместной количественной обработки и интерпретации комплекса  РВМ.

Во втором параграфе дается описание полной количественной интерпретации комплекса РВМ, основными элементами (этапами) которой являются 1) количественная обработка односкважинного радиоволнового зондирования (РВЗ-ОС) в рабочем интервале, 2) количественная обработка результатов межскважинного радиопросвечивания (РП-МС) рабочего интервала, 3) совместная обработка данных РВЗ-ОС и РП-МС, 4) геолого-геофизическая интерпретация РВМ на основе анализа результатов совместной обработки данных односкважинных и межскважинных радиоволновых измерений. Все четыре составляющие интерпретации, перечисленные выше, предполагают использование материалов ПРФ-исследований.

Задачей количественной обработки односкважинных радиоволновых зондирований рабочего интервала является определение коэффициента поглощения (КП) пород  РИ

для каждой скважины поискового полигона. Для определения   по  данным  РВЗ-ОС традиционно  используется формула для амплитуды поля Ez в однородной среде на оси электрического диполя в волновой зоне источника [Руководство, 1977].

  При проведении РВЗ-ОС в слоистых средах использование формулы для поля в волновой зоне в однородной среде позволяет получить лишь эффективное значении коэффициента поглощения . Истинный же коэффициент поглощения пород РИ можно получить, используя теоретическую зависимость от при определенных наборах значений других, менее значимых, параметров. Эта зависимость может быть построена в математической или графической форме. Пример  такой теоретической зависимости представлен в виде номограммы для модели слой-волновод  на рис. 6.  При известной связи между теоретическими значениями от , можно по значениям  , полученным из РВЗ-ОС, определять истинный коэффициент поглощения пород РИ.. Величина определяется по заданному набору длин разносов Z и соответствующих им значений измеренного поля методом наименьших квадратов. Ранее отмечалось, что односкважинные измерения выполнятся на частотах в 28 раз превышающих частоты РП-МС, между тем хорошо известен такой феномен, как частотная дисперсия материалов и горных пород [Петровский, 1964; Богородицкий и др., 1965; Каринский, Даев, Светов, 1997; Светов, Агеев, Лебедева, 1996 и др.]. Наиболее отчетливо она выражена для пород с УЭС от 700-1000 Ом⋅м и выше на частотах от первых сотен кГц до первых десятков МГц [Мишенин, Зинчук, Бондаренко, 2001; Зинчук, Бондаренко, Гарат, 2002 и др.]. Поэтому следующим этапом обработки РВЗ-ОС является учет частотной дисперсии, то есть пересчет величины , для частоты РВЗ-ОС, обозначенной , в значение , соответствующее частоте РП-МС, с помощью специальных номограмм или эмпирических формул частотных зависимостей  УЭС, диэлектрической проницаемости и КП для данной породы или группы пород.

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

Рис. 7. Частотная зависимость электрических параметров (, /0) и коэффициента поглощения  k”  мерзлых пород  (T=-10C)  Далдыно-Алакитского района:

1 – известняки, доломиты, мраморизованные доломиты S1ln, песчанники O1ss, 2 – глинистые известняки, мергели,  3 – кимберлиты с трубок Восток, Бобкова, Айхал и др., 4 – кимберлиты трубки Юбилейная. 

одного из районов поисков кимберлитовых трубок Якутии приведены на рис. 7 [Зинчук, Бондаренко, Гарат, 2002]. По материалам ПРФ-исследований мерзлых пород основных районов поисков кимберлитовых трубок Далдыно-Алакитского и Малоботуобинского, получены обобщенные эмпирические номограммы и формулы частотных зависимостей и взаимосвязи УЭС, диэлектрической проницаемости и КП для пород РИ, позволяющие опреде-лять значения УЭС, ДЭП и КП с учетом частотной дисперсии [Кеворкянц, Бондаренко, Абрамов 2002].  Подводя итог, весь процесс обработки результатов РВЗ-ОС можно представить в виде следующей после­довательности  операций:  1) определение  в породах рабочего интервала для каждой скважины; 2) пересчет (или перевод с помощью специальных номограмм) указанного пара­метра в значение истинного коэффициента поглощения пород рабочего интервала ; 3) пересчет значений с частоты РВМ-ОС в значения истинного коэффициента поглощений на частоте РП-МС с помощью обобщенных формул или номограмм частотной зависимости  КП  для данного типа пород.

Задачей количественной обработки межскважинного радиопросвечивания является определение интерпретационных параметров, необходимых для достоверных оценок  нормального поля в РИ каждой пары скважин. В процессе количественной обработки РП-МС определяются следующие интерпретациионные параметры: 1) амплитудный коэффициент E0 для изучаемого участка;  2) коэффициент поправки за модель для амплитудного множителя, определяемый по номограмме или рассчитываемый с помощью компьютерной программы; 3) среднеплощадные  эффективный и истинный коэффициенты поглощения в рабочем интервале в направлении простирания слоев  –  соответственно и ;

4)  кажущийся коэффициент поглощения  для  каждой пары  скважин  (j1, j2= =1, 2, …, N – порядковые номера скважин данной пары,  j1 j2 ; N – общее число скважин).

Интерпретационные формулы для оценки нормального поля основаны на приведенных в третьей главе приближенных аналитических представлениях (15) и (16) компоненты Ez поля поперечного электрического диполя  в базовых РВ-ГМ слоисто-анизотропной среды. Расчет и анализ нормального поля для частот 312 и 625 кГц при широком диапазоне изменений параметров РВ-ГМ обоих типов показали возможность его представления в виде следующих интерпретационных формул:

  (для РВ-ГМ типа 1) (17)

    (для РВ-ГМ типа 1) (18)

где r – расстояние  между  передающей  и  приемной  антеннами,  H  –  расстояние  от луча

просвечивания до кровли пород рабочего интервала в РВ-ГМ типа 1, h – мощность слоя-волновода в РВ-ГМ типа 2; -  коэффициент  поглощения  нижнего  полупространства  в

РВ-ГМ типа 1; - эффективный коэффициент слоя волновода в РВ-ГМ типа 2; - амплитудный множитель, слабо зависящий от ПРФ-параметров модели (в РВ-ГМ типа 1 он зависит также от H) и в связи с этим принимаемый в пределах изучаемого участка (при фиксированных значениях H в РВ-ГМ типа 1) за постоянную величину. Номограмма  величины , как функции от  двух параметров  r и  H, приведена на рис. 8.  Значения УЭС, коэффициентов анизотропии пород РИ и перекрывающей экранирующей толщи, приведенные в номограмме, наиболее реально отражают границы изменения указанных параметров в пределах Далдыно-Алакитского района поисков кимберлитовых трубок.

Зависимость выражений в правых частях соотношений (17) и (18) от ПРФ-параметров модели (вне экспоненты), как бы включена в множитель , получаемый при обработке РП-МС с использованием МНК, а параметр H (в РВ-ГМ типа 1) обычно фиксируется для заданного набора межскважинных измерений. Поэтому указанные выражения представляют собой окончательные интерпретационные формулы, используемые для оценок нормального поля и расчетов коэффициентов экранирования, в которых множители и рассматриваются как поправочные коэффициенты. Они определяются из результатов теоретических расчетов, проведенных с использованием данных ПРФ-идентификации модели, или оцениваются приближенно из таблиц, приведенных в диссертации, либо (для РВ-ГМ типа 1) номограмм вида, приведенного на рис. 8.

Среднеплощадные амплитудный множитель и коэффициент поглощения (или

) определяются по методу наименьших квадратов (для линейной регрессии)с итеративной фильтрацией  из  совокупности пар значений { rj, } j =1, 2,…N  (N – общее число обрабатываемых пар скважин), где rj – расстояние между центрами приемной и передающей антенн,– поправленная амплитуда поля. Применение фильтров позволяетисключить из обработки данные с аномально низкими значениями измеренного поля. 

Особенность совместной обработки данных межскважинных и односкважинных  измерений в анизотропно-слоистых средах заключается в том, что для их увязки, прежде всего, требуется количественная оценка анизотропии коэффициента поглощения пород РИ, равной . Определенный таким образом коэффициент анизотропии не будет иметь ясного физического смысла, если в правую часть последнего соотношения подставить эффективные значения КП.  Определение истинных значений по односкважинным измерениям изложено выше.  Значения - по межскважинным измерениям в РВ-ГМ типа 1 можно рассчитать из (17), если все остальные параметры в ней известны. Для определения истинного в РВ-ГМ типа 2 необходимо использовать теоретические расчеты  [выражение (16)]. Установив взаимосвязь между теоретическими значениями  и , можно по значениям второго из параметров, полученным из РП-МС, определить и, наоборот, по известным определить значения , необходимые для расчетов нормального поля  в РИ по интерпретационной формуле (18). Примеры теоретической зависимости от  приведены  в  виде номограмм на рис. 9.

Рис. 9.  Номограммы зависимости эффективного коэффициента поглощения  от  истинного  в  анизотропном слое-волноводе мощностью h

1, 1t  и  1n – соответственно диэлектрическая проницаемость, продольное и поперечное УЭС среды, граничащей со слоем-волноводом;  2, 2t  и  2n– соответственно диэлектрическая проницаемость, продольное и поперечное  УЭС  слоя-волновода.

Оценку коэффициента анизотропии можно получить по совокупности пар значений КП, включающих  значение , полученного из результатов РП-МС в j-ой паре, и значение , полученного из односкважинных измерений в каждой из скважин j-ой пары. Отобранные пары значений наносятся в виде точек на плоскость прямоугольной системы координат  с абсциссой и ординатой , и по этим точкам строится кривая регрессии  ( рис.10). 

Связь между и  можно получить и вычислительным способом, применив МНК. 

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

Рис. 10.  Кривые эмпирической взаимосвязи между коэффициентами поглощения ( по данным  РВМ-ОС) и (по данным РП-МС)  анизотропной толщи  пород рабочего интервала для северного (кривая 1) и южного (кривая 2) участков полигона Кюэллеэх-Бестях  Мало-Ботуобинского района  Западной Якутии (РВ-ГМ типа 1).

Примеры результатов определения для каждой отдельной скважины, полученных путем совместной обработки РП-МС, РВЗ-ОС и результатов ПРФ-исследований демонстрирует табл. 1  приложения 1 (фрагмент результатов обработки РВМ по участоку ЦЕНТР III  – РВ-ГМ слой-волновод).  Среднеплощадной коэффициент поглощения , полученный по данным РП- МС и средняя от совокупности значений по отдельным скважинам  для данного участка равны, соответственно, 0,0133 и  0,0138  (расхождение меньше 4%). Аналогичные параметры для участка Кюэллеэх-Бестях (РВ-ГМ типа 1) равны, соответственно, 0,0346 и 0,0339 Нп/м  (расхождение составляет чуть больше 2%), при этом оценки нормального поля, выполненные по указанным значениям КП при расстояниях между скважинами  r =250÷400 м, различаются соответственно на 10÷15 %.

Наличие следующих трех параметров для каждой пары скважин: значения , полученного по  данным РП-МС  в  j-ой  паре скважин,  и значений , , полученные для каждой из скважин той же пары (в результате совместной обработки односкважинных и межскважинных измерений), позволяет рассчитывать коэффициент экранирования Э j  для данной пары скважин с помощью соотношений:

    (для РВ-ГМ 1), 

  (22)    (для РВ-ГМ 2).

где  i (i = j1 , j2) – некоторые весовые коэффициенты, .  Пример получения  Э j по указанным трем параметрам для каждой пары скважин демонстрирует таблица 2 приложения 1  (фрагмент).

Детализированная проверка  на  перспективность.  Детализированная проверка предполагает выполнение следующих операций: 1) оценку коэффициента экранирования с помощью одного из соотношений (22);  2) выделение аномальных значений коэффициента экранирования  с учетом максимального ожидаемого уровня  геопомех, выражаемого в виде коэффициента экранирования  , для каждой (j-ой) пары скважин. В зависимости от соотношений параметров и , может быть установлен один из четырех уровней перспективности данного участка: 1 – неперспективный,  2 - условно неперспективный,  3 – условно перспективный, 4 – перспективный. Уровень перспективности РИ каждой ( j–ой) пары скважин определяется по критериям-неравенствам, которым будет удовлетворять значение Э j. Примеры оценки уровня перспективности для отдельных пар скважин, по которым получены относительно большие значения  Э j,  приведены  в  табл. 3. (приложение 1).

 

Рис. 11. Пример  обнаружения  и  локализации  кимберлитовых  трубок  с  помощью  РП-МС  на  участке  Бобкова Далдыно-Алакитского района:

–  поисковые и разведочные  скважины; 2 – скважины под заверку  аномалий  РП-МС;  3 –  луч  просвечивания  со  значением  коэффициента  экранирования  Э ;  4 – предполагаемый согласно бурению  контур поперечного сечения кимберлитового тела; 5 – контуры поперечных сечений вблизи от торцевой части кимберлитовых трубок  «40 лет Амакинской ГРЭ» (РВО №1),  «Радиоволновая» (РВО №2) и «ЦНИГРИ» (РВО №3) по данным РП-МС.

В третьем параграфе приводятся примеры применения интерпретационной системы РВМ на данных комплекса РВМ по участкам Кюэллеэх-Бестях (РВ-ГМ типа 1) и ЦЕНТР III (РВ-ГМ типа 2) Западной Якутии. Они представлены в виде сводных таблиц результатов обработки и планов участков с нанесенными на них значениями коэффициентов экранирования для каждой пары скважин. Приведены примеры локализации кимберлитовых трубок..

Наглядным примером обнаружения и локализации кимберлитовых трубок в результате применения комплекса РВМ и его интерпретации по предложенной автором методике является участок Бобкова, где было установлено наличие трех аномальных радиоволновых объектов (РВО) (рис. 11 – РВО-№1, РВО-№2 и РВО-№3). Заверка аномалий бурением показала, что данные объекты представляют собой кимберлитовые трубки. Детальная интерпретация данных РВМ при учете результатов разведочного и заверочного бурения позволила построить контуры горизонтального сечения верхней (приторцевой) части вышеупомянутых трубок, названных соответственно «40 лет Амакинской ГРЭ», «Радиоволновая» и «ЦНИГРИ». 

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

Для обработки данных межскважинного радиопросвечивания, выполненного в условиях нечетких моделей,  предложена модифицированная методика. Она заключается в замене линейной модели (в экспоненте) затухания поля в РИ нелинейной, в которой выделена линейная часть, соответствующая чистой модели, как полезная, и нелинейная, связанная со вкладом обходных волн, как помеховая, исключаемая на заключительном этапе обработки. Результаты применением модифицированной методики обработки данных на одном из участков Далдыно-Алакитского района приведены на рис. 12 в виде карты  распределения коэффициента поглощения на плоскости поискового горизонта, полученной методом томографической реконструкции. На карте хорошо выделяются зоны трещин и разломов, как области с повышенными значениями коэффициента поглощения.

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

Рис. 12. Результаты томографической обработки данных межскважинного радиопросвечивания на юго-западном подучастке одного из  поисковых участков Далдыно-Алакитского района. Карта распределения эффективного коэффициента поглощения  ()  пород поискового горизонта (рабочего интервала).

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

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

Проектируемые параметры: 1) минимальный  средний поперечный  размер поискового объекта (, где  S – параметр РГИО), 2) средний шаг разбуривания, то есть сторона элементарного квадрата разбуриваемойсети скважин (L ≈0,88÷0,83;  – среднее расстояние  между  скважинами),  3) основная рабочая частота РП-МС  рассчитываются из эталонных параметров с помощью следующих соотношений:

а)  экстраполяции , 

б)  подобия  ,

в которых  индексами «э» и «п» обозначены, соответственно, параметры эталонного участка и участка проектируемых поисковых работ.

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

  Параметры

Участки

  Эталонный

Проектный

Средний размер объекта,  м

147

  200

Рабочая частота радиопросве-

чивания, кГц

312

  156

Коэффициент поглощения, Нп/м

0.017

  0.015

Средний шаг сети скважин, м

230

  270

Из содержания данной главы следует:  предложены

1) методика интерпретации, основанная на учете типов РВ-ГМ и использовании односкважинных радиоволновых измерений, которая позволяет

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

- существенно повысить гарантию обнаружения поисковых объектов при относительно небольших коэффициентах экранирования (Э ≥ 5÷10), в то время как без данных ОС-измерений о наличии трубки можно уверенно говорить лишь при Э ≥ 20-30, что подтверждается анализом результатов, полученных на участках опытно-методических и поисковых работ и, в частности, обнаружением и четкой локализацией на участке Бобкова (Далдыно-Алакитский район) трех кимберлитовых  трубок, названных впоследствии Радиоволновая,  ЦНИГРИ  и  40 лет  Амакинской ГРЭ ;

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

3) способ проектирования поисковых работ (бурения и РВМ), основанный на

- учете результатов математического моделирования,

- непременном учете соответствия участков, в которых проектируются поиски одному 

из типов базовых радиоволновых геомоделей

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

ЗАКЛЮЧЕНИЕ

В диссертационной работе автором получены следующие

основные результаты:

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

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

3.  Для решения прямой задачи радиопросвечивания в слоисто-анизотропных средах, содержащих локальные неоднородности, разработаны теоретические основы в виде объемных и граничных сингулярных интегральных уравнений, объемных интегро-дифференциальных уравнений 2-го рода и получены приближенные аналитические решения в виде рядов несобственных интегралов или явных функций для моделей «экранированное полупространство» и «слой-волновод», содержащих поперечный круговой цилиндр повышенной проводимости.

4. Дана новая трактовка области, существенной для распространения радиоволн в поглощающих средах, названная радиогеоинформационной областью (РГИО), и предложен способ определения ее границ. Впервые сформулированы принципы типовой классификации моделей сред, изучаемых методом радиопросвечивания (РП), так называемых радиоволновых геомоделей (РВ-ГМ), на основе учета границ РГИО. В соответствии с указанными принципами определены основные типы РВ-ГМ  слоистых сред, в которых гарантируется эффективное применение РП с целью обнаружения там рудных тел или их фрагментов.

5. На основе приближенных аналитических (асимптотических) выражений для поля электрического диполя в основных типах радиоволновых геомоделей получены удобные интерпретационные формулы для расчетов нормального поля в основных типах РВ-ГМ  при поисках рудных тел в слоисто-анизотропных средах и формулы для оценки параметров поисковых объектов пластового типа.

7. Для повышения эффективности РВМ при поисках рудных месторождений в слоисто-анизотропных средах предложена оригинальная интерпретационная система  РВМ. Она включает раздельную и совместную интерпретации межскважинных и односкважинных радиоволновых измерений в рабочем интервале при использовании результатов петрофизических исследований, учете частотной дисперсии и анизотропии электрических параметров горных пород, слагающих поисковую толщу. Её применение позволяет существенно повысить (по сравнению с существующей методикой)  достоверность оценки нормального поля и эталонных коэффициентов экранирования  ,  для каждой исследуемой пары скважин, обеспечивая тем самым значительное (не менее, чем в 23 раза) повышение гарантии обнаружения поисковых объектов и надежности выявленных аномалий в условиях латеральной изменчивости петрофизических характеристик разреза.

8.  Предложен способ проектирования поисковых работ (бурения и РВМ) на основе учета результатов численного моделирования полей в слоисто-анизотропных средах, непременного учета типа РВ-ГМ поискового горизонта, границ РГИО, а также применения экстраполяции и критериев подобия к эталонным моделям. Его применение дает возможность более точного и скорого определения оптимальных метрологических характеристик (оптимальных густоты сети скважин и глубины бурения, рабочих частот РВМ, и т.п.), обеспечивающих наиболее эффективное применение комплекса РВМ при поисках экранирующих объектов (рудных тел) в слоисто-анизотропных средах.

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

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

ЦИТИРУЕМАЯ  ЛИТЕРАТУРА

1. Авдеев Д.Б. Разработка и применение метода интегральных уравнений для решения трех- мерных задач электроразведки. Автореф. докт. дисс. Троицк. 2002. 38 с.

2. Басс Ф.Г. Граничные условия для электромагнитного поля на поверхности с произвольным значением диэлектрической проницаемости. Ж. Радиотехника  и  электроника, 1960,

№3.  С. 389-392.

3. Бехтерева М.С. Некоторые результаты опробования скважинных радиоволновых изме-

рений при поисках кимберлитовых трубок под траппами. Тр. ЦНИГРИ. Вып. 137.  М., 

1978.  С. 76-82.

4 Богородицкий Н.П., Волокбинский Ю.М., Воробьев А.А. и  др. Теория диэлектриков.

«Энергия».  1965.  344 с.

5. Бондаренко А.Т., Ковалев Ю.Д.  Устройство и методика экспрессных измерений электрических параметров горных пород на естественном мерзлом керне.  Тр. ЦНИГРИ. Вып. 222,

1988. С. 53-58.

6. Бондаренко А.Т., Пыстин А.Б. Петрорадиофизические свойства карбонатных и магматических пород и их использование при поисках кимберлитовых трубок. Руды и металлы. 1995, №1. С. 69-76. 

7. Бондаренко В.М., Коваленко Н.Д., Тархов А.Г. Геофизические исследования урановых месторождений методом радиоволнового просвечивания. – «Изв. вузов, Сер. Геология и разведка», 1962, №2, с. 71-82. 

8. Борисов Б.Ф., Гуревич Г.Ф., Чигирина И.И. Методические указания по обработке и интерпретации результатов радиопросвечивания в анизотропных средах по редкой сети скважин. – М.: ЦНИГРИ, 1984.

9. Бреховских Л.М. Волны в слоистых средах. Изд. АН СССР. М. 1973.  260 с. 

10. Бурсиан В.Р. Теория электромагнитных полей, применяемых в геофизике. – Л.: Недра, 1972

10. Ваньян Л.Л. Основы электромагнитных зондирований. М. «Недра»,  1965. 168 с.

11. Васильев Е.Н. Алгоритмизация задач дифракции на основе интегральных уравнений. – В сб.: Научно-методические статьи по прикладной электродинамике. Вып. 1. М., «Высшая школа», 1977.  С. 94-128.

12. Грачев А.А. К вопросу об обработке результатов радиоволнового просвечивания. – «Разведочная геофизика». Вып. 14. М., «Недра», 1966.  С. 42-44.

13. Даев Д.С. Об обработке и истолковании данных радиоволнового просвечивания. Изв. вузов, «Геология и разведка». 1959,  №6.  С. 104-109.

14. Дмитриев В.И. Электромагнитные поля в слоисто-неоднородных средах. М.: Изд-во МГУ.  1969. 311 с.

15. Дмитриев В.И. Расчет электромагнитного поля в методе частотного зондирования. Сб. Вычислительные методы и программирование. Изд. МГУ.  М. 1965. С. 386-397. 

16. Дмитриев В.И., Захаров Е.В. Интегральные уравнения в краевых задачах электродинамики. Изд. Московского университета. 1987.168 с.

17. Дмитриев В.И., Федорова Э.А. Учет конечных длин передающей и приемнолй антенн в методе радиопросвечивания. В сб. «Численные методы в геофизических исследованиях». Изд. МГУ. 1979. С. 51-63.

18. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. Т.1. «Финансы и статистика». М. 1986. 366 с.

19. Жданов М.С., Спичак В.В.  Интегралы типа Стрэттона-Чу  для  неоднородных  сред  и некоторые их приложения к задачам геоэлектрики. В сб.: Математическое моделирование электромагнитных полей. М., ИЗМИРАН, 1983, с.4-25

20. Захаров Е.В., Ильин И.В. Интегральные представления электромагнитных полей в неоднородной слоистой среде. Изв. АН СССР, сер. «Физика Земли». 1970,  №8. С. 62-72.

21. Зингер Б.Ш., Файнберг Э. Б. Обобщенный итеративно-диссипативный метод расчета электромагнитных полей в неоднородной среде с учетом токов смещения. Физика Земли. 1995, №2. С. 69-72. 

22. Зинчук Н.Н., Бондаренко А.Т., Гарат М.Н. и др. Петрофизика кимберлитов и вмещающих пород. М.  Недра.  2002.  695 с.

23. Иванов В.Н., Петровский А.Д., Привезенцев В.И. Об использовании данных радиоволнового просвечивания при разведке медно-никелевого месторождения. – «Разведка и охрана недр». 1969,  №11. С. 29-34.

24. Каринский А.Д., Даев Д.С., Светов Б.С. и др.  Определение диэлектрической проницаемости горных пород по частотным зависимостям удельного электрического сопротивления на основе преобразования Гильберта. Изв. вузов. Геология и разведка. 1997, №5. С. 101-109.

25. Макаров Г.И., Новиков В.В., Рыбачек С.Т. Распространение электромагнитных волн  над земной поверхностью. Результаты исследований по международным геофизическим проектам. М.,  АН СССР.  1991. 197 с.

26. Мамаев В.Н. Особенности методики скважинного радиоволнового просвечивания в неоднородных вмещающих средах. Автореферат канд. дисс. М., МГРИ, 1978. 16 с.

27. Мамаев В.Н. Методика интерпретации результатов радиопросвечивания при оконтуривании промышленных руд, залегающих в низкоомных породах..- «Разведочная геофизика», вып. 54. М., «Недра», 1972, с. 78-86.

28. Микоев И.И. Методика применения радиоволновых методов при поисках кимберлитовых трубок Западной Якутии. Изд. Воронежского ун-та 2001. 31 с.

29. Мишенин С.Г., Зинчук Н.Н., Бондаренко А.Т. Петрофизические параметры кимберлитов, траппов и осадочных пород. В кн. Проблемы алмазной геологии и некоторые пути их решений. Изд. Воронежского ун-та  2001. С. 297-323.

30. Панкратов О.В., Авдеев Д.Б., Кувшинов А.В. Рассеяние электромагнитного поля в неоднородной земле. Решение прямой задачи. // Физика Земли. 1995.  №3. С. 17-25.

31. Петровский А.Д. О зависимости электрических свойств горных пород от частоты электромагнитного поля. Тр. ЦНИГРИ. Вып.59, 1964. С. 93-97.

32. Петровский А.Д. Методика вычисления волновых полей вблизи проводящих пластов и плоских границ раздела двух сред. М., ЦНИГРИ. 1967. 17 с.

34. Петровский А.Д. Об одном способе интерпретации результатов радиоволнового просвечивания. Тр. ЦНИГРИ. Вып. 33, 1959.

35. Петровский А.Д. Радиоволновые методы в подземной геофизике. М., «Недра» 1971. 323 с. Изд. 2-ое,  М., ЦНИГРИ,  2001. 290 с. 

35. Петровский А.Д., Бехтерева М.С., Привезенцев В.И. О геологической и экономической эффективности скважинного радиопросвечивания при поисках трубообразных тел в контрастно-слоистых средах. Тезисы докл. На Всесоюзном научно-техническом семинаре: Геофизические исследования скважин при поисках и разведке глубокозалегающих рудных месторождений. Г. Октябрьский, 1982. С. 45-46.

36. Руководство по радиоволновым методам скважинной и шахтной геофизики. Под ред.  А.Д. Петровского и А.А. Попова. М., «Недра» 1977. 335 с.

37. Савицкий А.П. Метод радиоволнового просвечивания в условиях малого различия  по сопротивлению пород и руд (на примере месторождений цветных металлов Восточного Забайкалья). Автореферат канд. дисс. Л., ЛГУ, 1969. 19 с.

38. Светов Б.С. Электродинамические основы квазистационарной геоэлектрики. – М.: Недра, 1984. 184 с.

39. Светов Б.С.  О частотной дисперсии электрических свойств среды. Ж. Физика Земли. 

  1992,  №4. С. 62-70.

40. Светов Б.С., Агеев В.В., Лебедева Н.А. Поляризуемость горных пород и феномен высокоразрешающей электроразведки. Геофизика. Ж. ЕАГО.  1996.  С. 42-52.

41. Светов Б.С., Губатенко  В.П. Аналитические решения задач электродинамики. М.: Наука, 1988. 

42. Соколов Ю.Н. Аналитический способ раздельного определения мощности и проводимости пластовых рудных тел по данным радиопросвечивания. Тр. ЦНИГРИ. Вып 134, 1975. С. 67-73. 

43. Табаровский Л.А. Применение метода интегральных уравнений в задачах геоэлектрики. Новосибирск: Наука, 1975. 139 с.

44. Табаровский Л.А.,  Эпов М.И. Электромагнитные поля гармонических источников в слоистых анизотропных средах.  «Геология и геофизика». Изд. «Наука», Сибирское отделение. 1977,  №1.

45. Табаровский Л.А., Дубман А.Л. Метод разнесенных вспомогательных источников в математическом моделировании задач геоэлектрики. Ж. Геология и геофизика. 1989, №2. с.132-139.

46. Талалов А.Д,  Даев Д.С. О структурном механизме частотной дисперсии электрических

свойств гетерогенных горных пород. Ж. Физика Земли 1996,  №8. С. 56-66.

47. Тархов А.Г. Радиоволновой метод электроразведки. Тр. МГРИ, т. 28, 1955.

48. Тихонов А.Н. о распространении переменного электромагнитного поля  в  слоистой анизотропной среде. Докл. АН СССР. 1959. Т. 126,  №5. С. 967-970.

49. Фейнберг Е.Л. Распространение радиоволн вокруг земной поверхности. М.. изд. АН СССР. 1961. 495 с. 

50. Фелсен Л., Маркувиц М. Излучение и рассеяние волн.  Т. 1. М., МИР. 1978. 547 с.

60. Фельдман Т.А.  Поле диполя в слое с идеальными границами. Тр. ЦНИГРИ. Вып. 116, «Рудная геофизика». М. 1974. С. 136-148.

61. Фролов А.Д., Истратов В.А., Лысов М.Г., Остапчук С.И. Радиоволновая геоинтроскопия мерзлых пород. Материалы второй конференции геокриологов России, Т. 4,  Инженерная геокриология. – М.: Изд-во МГУ, 2001.  С. 290-299.

62. Хмелевской В.К. Электроразведка. М., изд. МГУ, 1984, 420 с.

63. Четаев Д.Н. Новый метод решения основной задачи теории дипольных электромагнитных зондирований// Геология и геофизика. 1962. .№2. С. 118-122.

64. Четаев Д.Н. О электромагнитных потенциалах в слоисто-анизотропных средах. Изв. АН СССР, сер. «Физика Земли». 1966,  №10.  С. 48-61.

65. Anderson W.L. Computer Program Numerical integration of related Hankel transforms of orders 0 and 1 by adaptive digital filtering.  Geophysics.  Vol.. 44,  No 7 JULY 1979. P. 1287-1305.

66. Avdeev D.B., Kuvshinov A.V., Pankratov O.V., Newman G.A. High-performance three-dimensional electromagnetic modeling using modified Neuman series/ Wide-band numerical solution and examples // J. Geomagn. Geoelectr. 1997. V. 49. P. 1519-1539 

67. Christensen N.B. Optimal Fast Hankel Transform filters. Geophys. Prospect., 38, 1990. Р. 545-568.

68. Kong  J. A. Electromagnetic fields due to dipole antennas over stratified anisotropic media. Geophysics, Vol. 37, N 6 (December, 1972). Р. 985-986.

69. Lytle, R.J., Lager, D.L., Laine, E. F., Salisbury, J.D. and Okada, J.T.. Fluid-flow monitoring using electromagnetic probing: Geopys. Prosp. 29, 1981. Р. 627-638.

70. Mller C. Grundprobleme der matematischen Theorie der electromagnetischen Schwingungen. Berlin/ 1957. 

71. Johnson, H., Richardson, B., Pizziara, J., and Salisbury, W.W. Radio-frequency geothomography for remotely probing the interiors of operating mini and commercial sized oil shale restores. Geophysics, 49, 1984. Р. 1288-1300. 

СПИСОК РАБОТ АВТОРА

1. Кеворкянц С.С.  К расчету поля поперечного гармонического электрического диполя в присутствии анизотропного слоя, расположенного в однородной анизотропной среде. Труды ЦНИГРИ. Вып. 145. М. 1979. С. 52-58.

2. Кеворкянц С.С. К расчету электромагнитных полей гармонических диполей в трехслойной анизотропной среде. Труды  ЦНИГРИ. Вып. 161. М. 1981. С. 72-73.

3. Кеворкянц С.С. О представлении электромагнитного поля в однородной анизотропной среде  через вектор Герца магнитного типа. Тезисы к докл. на конференции молодых ученых и специалистов ЦНИГРИ. Изд. ЦНИГРИ, М.,1982. С. 133-135.

4. Кеворкянц С.С. О поле вертикального гармонического электрического диполя в анизотропной среде в присутствии проводящего цилиндра. Там же. С. 138-140.

5. Кеворкянц С.С. О способах повышения эффективности интерпретации данных межскважинного радиопросвечивания. Там же. С. 149-151.

4. Кеворкянц С.С. Об условиях Лоренца для электромагнитных потенциалов в одноосно анизотропных средах. Изв. АН СССР, «Физика Земли», М., 1983,  №1,  с. 69-73.

5. Кеворкянц С.С. Расчет поля вертикального электрического диполя в присутствии поперечного кругового цилиндра в слое-волноводе и полупространстве с идеально отражающими границами. Труды ЦНИГРИ. Вып. 179. .М. 1983.

6. Кеворкянц С.С. К обоснованию способа учета волновой анизотропии горных пород путем односкважинных радиоволновых измерений. Тезисы докл. в сб. «Критерии прогноза и пути повышения эффективности геологоразведочных работ на месторождениях благородных, цветных металлов и алмазов». М., ЦНИГРИ, 1983

7. Кеворкянц С.С. Теория и способы межскважинного радиопросвечивания при поисках кимберлитовых трубок Западной Якутии. Автореф. к канд. дис. М., МГРИ, 1986. 18 с.

8. Кеворкянц С.С. К постановке граничной задачи электродинамики слоистой анизотропной среды. Физика Земли. 1987. №8. С. 62-67.

9. Кухарев В.Ф., Кеворкянц С.С. Расчет поля «Ez» в слоисто-анизотропной среде применительно к задачам межскважинного радиопросвечивания. Тр. ЦНИГРИ. Вып. 222. М., 1988, с. 34-42.

10. Кеворкянц С.С. К определению полей произвольных диполей в электромагнитно-анизотропной слоистой среде. Изв. АН СССР. Физика Земли. М. 1991,  №7. С. 52-59.

11. Кеворкянц С.С. Кухарев В.Ф. Об одном способе расчета полей погруженных дипольных источников в слоисто-анизотропных средах. Изв. РАН. Физика Земли. М. 1992,  №5. С. 71-78.

12. Кеворкянц С.С. Интегральные уравнения электродинамики в задачах дифракции на ограниченном теле в слоисто-неоднородной анизотропной среде. Изв. РАН. Физика Земли.  М. 1992, №5. С. 59-71.

13. Кеворкянц С.С. О тензорной нормировке функции Грина в граничных интегральных уравнениях электродинамики слоистых одноосно-анизотропных сред. Изв. РАН. Физика Земли. М. 1995,  №10. С. 52-60. 

14. Кеворкянц С.С. К определению тензорных функций Грина электродинамической задачи слоисто-анизотропной среды. Физика Земли. М. 2000, №10.  С. 84-92.

15. Кеворкянц С.С. Дифракция радиоволн на поперечном круговом цилиндре в слоисто-анизотропных моделях месторождений кимберлитовых трубок Якутии. Геофизика. 2004, №2. С. 58-69.

16. Кеворкянц С.С. Объемные и граничные интегральные уравнения второго рода в задачах дифракции на ограниченном теле, помещенном в слоисто-анизотропную среду. Физика Земли. М. 2007. № 3. С. 4-10.

17. Кеворкянц С.С., Абрамов В.Ю. Компьютерный метод выделения нормального поля при межскважинном (межвыработочном) просвечивании в локально неоднородных георазрезах. Отечественная геология. 1994, №4. С. 54-58.

18. Кеворкянц С.С., Абрамов В.Ю., Ковалев Ю.Д. Скважинный радиоволновой комплекс при поисках кимберлитовых трубок в Западной Якутии. Геофизика. 2005, №3, с. 56-64.

19. Кеворкянц С.С., Бехтерева М.С. Об интерпретации результатов радиопросвечивания в слабопоглощающей среде с экранирующими границами. Тез. докл. Конференции молодых ученых  и  специалистов  ЦНИГРИ:  Вопросы  геологии,  прогнозирования,  методики поисков  и  разведки  месторождений  золота, цветных  металлов  и  алмазов.  М.  1984.  С. 55-56.

20. Кеворкянц С.С., Бондаренко А.Т., Абрамов В.Ю. Способ учета частотной дисперсии электрических свойств мерзлых карбонатных пород в интерпретационном комплексе скважинной радиогеоразведки при поисках кимберлитовых трубок Западной Якутии. В кн. Проблемы прогнозирования, поисков и изучения месторождений полезных ископаемых на пороге XXI века.: Материалы региональной научно-практической конференции., г. Мирный, апр. 2003 г. Воронеж 2003. С. 573-576. 

21. Кеворкянц С.С.,. Виноградова М.Н.  Решение прямой задачи радиогеоразведки для слабых электромагнитных полей от дипольных источников при изучении слоисто-анизотропных сред. Тезисы докл. на годичной сессии ученого совета «НТД-96-ЦНИГРИ», 1-4 апреля 1997. М. 1997.

22. Кеворкянц С.С., Ковалев Ю.Д., Мамаев В.Н., Фролов А.С. Результаты радиопросвечивания на угольных месторождениях КНР. // НТД-92-ЦНИГРИ. Материалы научно-практической конференции. 12-14 мая 1993.  М., 1993. С. 68.

23. Кеворкянц С.С., Коновалова Т.М. Объемная локализация объекта на основе выделения информационных зон радиопросвечивания. Руды и металлы  2000, №5. С.77-81.

24. Кеворкянц С.С., Лаптев М.М., Абрамов В.Ю. Интерпретация радиопросвечивания при оконтуривании пластообразных рудных залежей. Геофизика. М., ЕАГО 2003,  №3. С. 58-63.

25. Петровский  А.Д. , Бехтерева М.С. , Бондаренко А.Т., Кеворкянц С.С. и др.  Основы применения скважинных РВМ для подземного картирования и обнаружения неоднородностей в контрастно-слоистой анизотропной среде.  Тр. XXX Международного геофизического симпозиума. Тезисы докл. Ч. 1. М., 1985.  С.  65-75..

26. Петровский А.Д., Мамаев В.Н., Кеворкянц С.С. и др. О возможностях оптимизированной системы технологий радиогеоразведки на месторождениях нефти и газа // Руды и металлы. - 1997. – №4. – С. 72. 


1 В отмеченной работе рассмотрены градиентные среды с  анизотропией в различных системах координат.




© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.