СПОСОБ РАДИОВОЛНОВОГО МЕЖСКВАЖИННОГО ПРОСВЕЧИВАНИЯ Российский патент 1997 года по МПК G01V3/12 

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

Известны способы обработки данных радиоволнового просвечивания, основанные на выделении аномалий в межскважинном пространстве по относительному уменьшению напряженности поля на наблюденной кривой ("Способ тени" [1, с. 150] ), путем сопоставления наблюденной кривой напряженности поля с кривыми напряженности нормального поля ("Способ сравнения" [2, с. 151]).

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

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

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

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

Наиболее близким к изобретению (его прототипом) является способ обработки данных радиоволнового просвечивания [3] заключающийся в получении полевых данных, оценивании уровня шума и его вычитании из полученных результатов по формуле
Е SQRT(Ue•Ue Un•Un) (1),
где Ue уровень сигнала,
Un уровень шума, учете траектории ствола скважины (коррекции данных по глубине), выделении лучей, удовлетворяющих условию малости радиальной компоненты зарегистрированного поля по формуле

где модуль волнового числа,
R расстояние передатчик-приемник,
Q1 угол между осью передающего диполя и лучом,
Q2 угол между осью приемного диполя и лучом, вычислении коэффициентов поглощения для этих лучей, построении зависимости коэффициентов поглощения от угла Q1 в вертикальной плоскости, оценке параметров анизотропии по этой зависимости для полярной компоненты электромагнитной волны, расчете величины нормального электромагнитного поля с учетом анизотропии по системе номограмм, расчете интерпретационного параметра (например, коэффициента экранирования) путем деления значения нормального поля на зарегистрированное значение при каждом измерении, вынесении полученных величин на обобщенную плоскость, построении изолиний интерпретационного параметра на обобщенной плоскости, выделении по ним аномалии и анализе ее формы, выделении прямолинейных участков изолиний в пределах аномалии, содержащих лучи, пересекающиеся в нижнем полупространстве, по критерию попадания угла наклона этих участков в диапазон углов О 90 градусов, переносе лучей с выделенных участков на межскважинный разрез так, чтобы оконтуриваемая ими область имела минимально возможную площадь (способ НЭП), и суждении на основе изучения указанной области о положении в межскважинном пространстве аномалиеобразующего объекта [3, с. 10-41] Как вспомогательный критерий для выделения лучей, пересекающихся в нижнем полупространстве, может быть также использовано попадание угла наклона касательной к изолинии в диапазон 0-90 градусов.

Недостатками известного способа является то, что процесс расчета нормального поля ведется с учетом влияния анизотропии только на одну полярную компоненту регистрируемой волны и только в вертикальной плоскости. Вычислительный алгоритм не предусматривает итерационной корректировки параметров анизотропии вмещающей среды, в связи с чем возникают существенные ошибки в интерпретации, а точная локализация аномалиеобразующего объекта в пространстве становится невозможной [3, с. 40]
Это обусловлено тем, что в силу слоистости горной среды на значение сигнала в точке приема существенным образом влияет угол встречи луча с плоскостью слоев. Однако подобная зависимость, в свою очередь, различна для различных расстояний от источника до приемника. На больших удалениях вклад от анизотропии в принятый сигнал может в несколько раз превышать вклад от объекта поиска.

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

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

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

К числу его недостатков следует отнести присутствующую всегда неоднозначность результатов, связанную с внутренними недостатками алгоритмов локализации и в частности способа НЭП.

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

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

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

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

Поиск и учет анизотропии в заявляемом способе разделяется на несколько этапов.

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

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

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

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

На фиг. 1 представлена блок-схема итерационного процесса учета и компенсации анизотропии и изменения коэффициента установки; на фиг. 2 - диаграмма углового распределения коэффициентов поглощения для полярной компоненты; на фиг. 3 инегральная диаграмма углового распределения коэффициентов поглощения; на фиг. 4 обобщенная плоскость коэффициентов экранирования; на фиг. 5 геоэлектрический разрез с результатами локализации по способу НЭП; на фиг. 6 геоэлектрический разрез с результатами локализации способом итерационных томографических алгоритмов (геотомограмма); на фиг. 7 - геотомограмма из эталонного каталога; на фиг. 8 кривые дифракционной составляющей электромагнитного поля для выбранной модели объекта; на фиг. 9 и 10 оценка электрического сопротивления аномального объекта: на фиг. 9 - график измерения сопротивления по линии L-M; на фиг. 10 диаграммы зависимости коэффициента экранирования от угла падения волны на объект.

Для существенного повышения точности интерпретации за счет раздельной обработки данных, описываемых оптико-геометрической и дифракционной волновыми моделями, в заявляемый способ вводится операция селекции исходных лучей по критерию угла наклона касательной к изолинии на обобщенной плоскости. Было теоретически и экспериментально показано, что лучи, несущие дифракционные искажения имеют углы наклона касательной к изолинии на обобщенной плоскости в пределах 90-180 и 270-360 градусов и, следовательно, могут быть выделены на основе этого критерия.

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

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

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

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

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

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

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

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

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

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

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

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

Противоречия могут быть разрешены лишь итерационно, многократной переобработкой полевого материала.

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

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

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

Отсутствие процесса сходимости означает, что примененная интерпретационная модель неверна и нуждается в замене.

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

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

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

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

6. Вычисляют нормальное поле как векторную сумму полярной и радиальной компонент с учетом анизотропии.

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

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

9. Выделяют аномалию и анализируют ее форму.

10. Рассчитывают угол наклона касательной к изолиниям и выделяют те участки изолиний, где значение угла попадает в диапазон О-90 и 180-270 градусов.

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

12. Судят на основании этого контура о границах аномального объекта в пространстве.

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

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

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

16. Строят набор зависимостей интерпретационного параметра от угла падения волны на объект для некоторого набора прямолинейных профилей.

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

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

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

20. Исключают из обработки интерференционные экстремумы по степени расхождения или закономерного изменения изолиний.

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

22. Многократно повторяя операции по пп. 19-21, строят ряды промежуточных результатов и определяют пределы их сходимости.

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

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

Расстояние между скважинами 1460 м глубина первой скважины 1240 м второй 980 м. Интервал перемещения излучателя в скважине 80-1000 м, приемника - 120-800 м. Шаг между стоянками излучателя 40 м, между стоянками приемника 20 м. Количество замеров 984.

Уровень шума был оценен в 0,23 мВ и вычтен по формуле (1).

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

По этим кривым были определены: для полярной компоненты угол наклона плоскости анизотропии Ао -7 градусов (точка F1), коэффициент анизотропии равный 0,15; для интегральной зависимости эти значения составили: Ао -10 градусов, (точка F2), коэффициент анизотропии 0,04.

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

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

Форма, размеры и положение на межскважинном разрезе аномалиеобразующего объекта первоначально были оценены по способу НЭП (фиг. 5). Интерпретация замкнутых контуров со значением коэффициентов экранирования "50", "100" и выше позволила сделать вывод о том, что аномалия образована двумя рудными телами, отстоящими друг от друга на значительное расстояние.

Для оценки такого предположения была рассчитана томограмма (фиг. 6). Ее сопоставление с эталоном (фиг. 7) позволило подтвердить истинность выбранной интерпретационной модели.

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

Определенные таким образом кромки объекта были вынесены на межскважинный разрез (фиг. 6, точка S1).

Сопротивление аномального объекта определено по линии L-М на межскважинном разрезе (фиг. 6). На фиг. 9 показаны расчетные графики изменения электрического сопротивления по этой линии, которое определено с помощью расчетных диаграмм зависимости коэффициента экранирования от угла падения электромагнитной волны на объекта (фиг. 10).

Удельное электрическое сопротивление выделенного объекта составило 60-80 Омм.

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

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

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

1. Способ радиоволнового межскважинного просвечивания, включающий возбуждение и регистрацию электромагнитного поля, расчет с учетом анизотропии среды нормального электромагнитного поля для полярной компоненты поля в вертикальной плоскости, расчет для полученного нормального поля интерпретационных параметров и суждение по ним о положении в межскважинном пространстве аномалиеобразующего объекта, отличающийся тем, что дополнительно при расчете нормального электромагнитного поля учитывают радиальную компоненту поля в вертикальной плоскости и полярную и радиальную компоненты в горизонтальной плоскости путем разграничения областей преимущественного действия полярной и радиальной компонент электромагнитного поля, построение отдельно в каждой выделенной области зависимости коэффициента поглощения от угла между лучом и плоскостью, перпендикулярной усредненной плоскости искривления стволов скважин и параллельной прямой, соединяющей их устья, построение на основе полученных зависимостей интегральной угловой зависимости коэффициентов поглощения для всей области исследований, по которой производят определение характеристик анизотропии среды отдельно для каждой компоненты электромагнитного поля, а также угловую зависимость коэффициента установки, вычисляют полное нормальное поле как векторную сумму полярной и радиальной компонент, сопоставляют его с зарегистрированным полем, оценивают степень их совпадения, корректируют характеристики анизотропии среды раздельно для полярной и радиальной компоненты и угловое распределение коэффициента установки, итерационно повторяют циклы расчетов, по результатам которых находят асимптотическую многомерную зависимость степени совпадения нормального и зарегистрированного электромагнитного поля от характеристик анизотропии среды и углового распределения коэффициента установки, определяют предел сходимости указанной зависимости, по достижении которого соответствующее ему нормальное электромагнитное поле принимается за итоговое, и на его основе производят расчет для всех зарегистрированных лучей интерпретационных параметров, по которым судят об истинной форме, размерах, электрических свойствах и положении в межскважинном пространстве аномалиеобразующего объекта. 2. Способ по п.1, отличающийся тем, что о границах аномалиеобразующего объекта в межскважинном пространстве судят по положению на межскважинном разрезе замкнутого контура наименьшей площади, который определяют путем расчета функции угла наклона касательной для каждой изолинии интерпретационного параметра на обобщенной плоскости, выделения участков изолиний, свободных от дифракционных и интерференционных искажений, по условию попадания угла наклона касательной к изолинии в диапазон 0 90^o и 180 - 270^o и переноса на межскважинный разрез с обобщенной плоскости лучей, попавших в выделенные участки изолиний, точки пересечения которых образуют на межскважинном разрезе замкнутый контур наименьшей площади. 3. Способ по п.1, отличающийся тем, что о положении в пространстве кромок аномалиеобразующего объекта судят по степени совпадения теоретической и зарегистрированной зависимостей амплитуды диффракционной составляющей зарегистрированного электромагнитного поля от взаимного расположения точек приемопередачи, которую определяют путем расчета функции угла наклона касательной для каждой изолинии интерпретационного параметра на обобщенной плоскости, выделения участков изолиний, несущих дифракционные и интерференционные искажения, по условию попадания угла наклона касательной в диапазон 90 180^o и 270 360^o, построения с использованием лучей из выделенных участков зависимостей амплитуды дифракционной составляющей зарегистрированного поля от взаимного расположения точек приемопередачи и построения теоретического каталога указанных зависимостей для различных положений в пространстве аномалиеобразующего объекта. 4. Способ по п.1, отличающийся тем, что о характерной форме и размерах объекта судят по степени совпадения томограммы интерпретационного параметра, построенной по данным зарегистрированного электромагнитного поля, с одной из расчетных томограмм из эталонного каталога. 5. Способ по п.1, отличающийся тем, что о численных значениях электрических свойств и мощности аномалиеобразующего объекта судят по степени совпадения теоретической и экспериментальной диаграмм зависимостей значения интерпретационного параметра от угла падения электромагнитной волны на данный объект, для чего рассчитывают каталог теоретических диаграмм этой зависимости для найденного положения в пространстве указанного объекта и набора различных значений его физических и геометрических характеристик.