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

1 Настоящее изобретение относится к способу обнаружения объектов в земле с использованием георадиолокаторов. Более конкретно,изобретение относится к установке для получения информации, указывающей на глубину и ориентацию объекта, находящегося под поверхностью земли, причем указанная установка предназначена для использования электромагнитного излучения, передаваемого и принимаемого антенной системой, связанной с установкой, содержащей передатчик и приемник соответственно для формирования электромагнитного излучения во взаимосвязи с упомянутой антенной системой и для приема электромагнитного излучения, отраженного от объекта,также во взаимосвязи с упомянутой антенной системой. Известны многочисленные технические решения и варианты осуществления таких наземных радиолокационных установок и различные варианты осуществления способов получения информации о глубине и положении указанного объекта путем его обнаружения с использованием наземных радиолокаторов. Примеры подобных технических решений описаны в следующих патентах США: 5 339 080,5 192 952, 5 499 029, 4 504 033, 5 352 095, 4 430 653,4 698 634, 4 062 010, 4 839 654 и 5 130 711. В частности, в патенте США 4 728 897 описан способ,использующий антенную систему, в которой набор диполей расположен симметрично относительно общего центра и выполнен с возможностью качания механическим или электронным способом для генерирования и приема электромагнитного излучения. Данная технология обеспечивает возможность сбора необходимой информации о положении и ориентации объекта рассматриваемого типа, однако данный способ имеет потенциальные возможности для усовершенствования и модификации для повышения эффективности генерирования информации упомянутого типа и, тем самым, для обеспечения более надежной идентификации глубины и ориентации обнаруживаемого объекта. Задачей настоящего изобретения является создание наземной радиолокационной установки вышеописанного типа, которая обеспечивает улучшенные функциональные возможности в отношении идентификации объектов рассматриваемого класса по сравнению с техническими решениями, известными из предшествующего уровня техники, и, следовательно, обеспечивает реализацию улучшенного способа измерений по сравнению с техническими решениями, известными из предшествующего уровня техники. Указанный результат, а также иные цели,преимущества, признаки настоящего изобретения, поясняемые в нижеследующем детальном описании, достигаются с помощью установки вышеупомянутого типа, характеризующейся тем, что она предназначена для использования электромагнитного излучения, передаваемого и принимаемого антенной системой, связанной с уста 002451 2 новкой, и содержит передатчик и приемник соответственно для генерирования электромагнитного излучения совместно с упомянутой антенной системой, а также для приема электромагнитного излучения, отраженного от объекта,во взаимосвязи с упомянутой антенной системой, отличающейся тем, что антенная система содержит множество индивидуальных дипольных антенн, которые размещены по отношению к геометрическому центру антенной системы так, что центры различных дипольных антенн смещены относительно геометрического центра антенной системы, и установка имеет средство для вращения механически или электронным способом антенной системы вокруг или относительно геометрического центра антенной системы. Реализация антенной системы, соответствующей настоящему изобретению, в виде множества дипольных антенн, смещенных относительно геометрического центра антенной системы, обеспечивает существенно более сложную характеристику излучаемого электромагнитного излучения и в то же время более точный и детальный прием электромагнитного излучения,что обеспечивает возможность достижения существенного усовершенствования способа измерения. Антенная система, соответствующая изобретению, может быть выполнена, согласно двум альтернативным вариантам осуществления, таким образом, что индивидуальные дипольные антенны расположены радиально от геометрического центра антенной системы,причем центры различных дипольных антенн расположены со смещением относительно геометрического центра антенной системы, в частности, центры различных дипольных антенн расположены на окружности или на различных окружностях с разными радиусами, причем геометрический центр антенной системы совпадает с центрами упомянутых окружностей, либо различные дипольные антенны могут, как вариант,формировать треугольную конфигурацию, расположенную симметрично относительно геометрического центра антенной системы, т.е. диполи расположены треугольником, а вершина треугольника расположена в геометрическом центре антенной системы, либо составную треугольную конфигурацию, расположенную симметрично относительно геометрического центра антенной системы. В качестве альтернативы эти два варианта выполнения дипольной антенны могут быть объединены с обеспечением другой диаграммы направленности, либо могут также использоваться другие конфигурации, отличные от треугольной, например конфигурации в форме квадрата, или многоугольника, или комбинации нескольких различающихся многоугольных конфигураций. Индивидуальные антенные элементы антенной системы установки, соответствующей настоящему изобретению, могут образовывать 3 параллельно ориентированные наборы передатчика и приемника, или, как вариант, индивидуальные антенные элементы антенной системы могут быть одинаково поляризованы или ортогонально поляризованы. Эти специальные признаки позволяют использовать широкое разнообразие вариантов осуществления. Установка, соответствующая настоящему изобретению, может содержать средство для обеспечения вращения антенной системы и, тем самым, обеспечивать вращения антенной системы с помощью механических или электронных средств. Предпочтительно в настоящем изобретении используется электронное вращение для обеспечения управления вращением с высокой точностью при существенном снижении полного веса установки. Электронное вращение может быть обеспечено введением сдвига между отдельными антенными элементами антенной системы, что обеспечивает поляризацию электромагнитного поля относительно геометрического центра с угловыми приращениями в типовом случае порядка 22,5, 30, 36 или долей или целых кратных величин этих углов. Антенная система установки, соответствующей настоящему изобретению, может в соответствии с одним из вариантов осуществления генерировать электромагнитное излучение и принимать электромагнитное излучение на нескольких отдельных частотах, предпочтительно в диапазоне от 100 МГц до 1 ГГц шагами по 5 МГц. При этом обеспечивается преимущество, состоящее в обеспечении широкого спектра отражений. Установка, соответствующая настоящему изобретению, может содержать средство обработки сигналов для измерения передаточной функции, например передаточной функции по напряжению, передаточной функции по току или их комбинаций, или передаточных функций по мощности для сигналов, соответствующих излучению, переданному антенной системой, по отношению к излучению, принятому антенной системой. В настоящем изобретении предпочтительно измеряется передаточная функция по напряжению для связанных значений углового изменения антенной системы и частоты сигнала,излученного антенной системой, и сигнала,принятого антенной системой. Средство обработки сигналов установки,соответствующей настоящему изобретению,может выполнять преобразование из частотной области во временную область с использованием преобразования Фурье или с применением математической экспоненциальной модели с соответствующей рациональной передаточной функцией для передаточной функции, связанной с конкретным углом, для генерирования непрерывной во времени функции для каждого угла, которая вычисляется в предварительно определенное число дискретных моментов времени при идентичных моментах вычисления для 4 каждого угла. Кроме того, средство обработки сигналов может выполнять преобразование из частотной области во временную область путем преобразования Фурье или с использованием математической экспоненциальной модели с соответствующей рациональной передаточной функцией для каждой угловой гармоники, зависящей от времени непрерывной во времени функции в предварительно определенное количество дискретных моментов времени при идентичных моментах вычислений для каждой угловой гармоники, генерирование набора чисел,соответствующих математической функции,имеющей дискретные угловые гармоники и дискретные моменты времени в качестве независимых переменных. Это обеспечивает надежное определение уровня и задержек по времени для основных отражателей даже при отсутствии идеальных условий. Средство для обработки сигналом установки, соответствующей настоящему изобретению, может выполнять анализ сигналов посредством преобразования Фурье над каждым из используемых моментов, связанных с каждым углом, непрерывно во времени, причем функцияиспользует время в качестве постоянной и угол в качестве независимой переменной для генерирования представления в области углов. Кроме того, средство для обработки сигналов может выполнять анализ сигналов с помощью преобразования Фурье для каждой частоты, причем частота используется как постоянная, а угол как независимая переменная для генерирования набора чисел с использованием угловых гармоник и частоты измерений в качестве независимых переменных. Назначением применения преобразования Фурье является выделение содержимого угловых гармоник сигнала при различных временных задержках (и глубинах). В установке, соответствующей настоящему изобретению, средство обработки сигналов может выполнять различные операции над обрабатываемыми сигналами. Средство обработки сигналов может предпочтительно выполнять сканирование локальных и глобальных максимумов, указывающих периодические по углу отражения от объектов, при представлении в области углов, причем пиковые значения указывают на временную задержку соответствующих отражений от объекта. Кроме того, средство обработки сигналов может выполнять вычисления углового положения относительно антенной системы для периодических по углу отражений,а также средство обработки сигналов может осуществлять использование собранных результатов измерений, выполненных в процессе перемещения по горизонтали, для подавления мешающих отражений. Настоящее изобретение также относится к способу генерирования информации, указывающей на глубину и ориентацию объекта, расположенного под поверхностью земли, причем 5 этот способ характеризуется тем, что включает использование электромагнитного излучения,передаваемого и принимаемого антенной системой посредством передатчика и приемника соответственно для формирования во взаимосвязи с антенной системой электромагнитного излучения и для приема электромагнитного излучения, отраженного от объекта, с помощью антенной системы, и отличается тем, что для передачи и приема электромагнитного излучения используется множество индивидуальных дипольных антенн, которые размещены по отношению к геометрическому центру антенной системы так, что центры антенных элементов смещены относительно геометрического центра антенной системы, и антенная система вращается, механически или электронным способом,вокруг или относительно геометрического центра антенной системы. Способ, соответствующий изобретению, предусматривает реализацию антенной системы в виде множества дипольных антенн, смещенных относительно геометрического центра антенной системы, что обеспечивает существенно более сложную характеристику излучаемого электромагнитного излучения и в то же время более точный и детальный прием электромагнитного излучения, что обеспечивает возможность достижения существенного усовершенствования способа измерения. В способе, соответствующем настоящему изобретению, антенная система может быть выполнена, согласно двум альтернативным вариантам осуществления, таким образом, что индивидуальные дипольные антенны расположены радиально от геометрического центра антенной системы, причем центры различных дипольных антенн расположены со смещением относительно геометрического центра антенной системы, в частности, центры различных дипольных антенн расположены на окружности или на различных окружностях с разными радиусами,причем геометрический центр антенной системы совпадает с центрами упомянутых окружностей, либо различные дипольные антенны могут, как вариант, формировать треугольную конфигурацию, расположенную симметрично относительно геометрического центра антенной системы, т.е. диполи расположены треугольником, а вершина треугольника расположена в геометрическом центре антенной системы, либо составную треугольную конфигурацию, расположенную симметрично относительно геометрического центра антенной системы. В качестве альтернативы эти два варианта выполнения дипольной антенны могут быть объединены с обеспечением другой диаграммы направленности, либо могут также использоваться другие конфигурации, отличные от треугольной, например конфигурации в форме квадрата, или многоугольника, или комбинации нескольких различающихся многоугольных конфигураций. 6 В способе, соответствующем настоящему изобретению, индивидуальные антенные элементы антенной системы установки, особенно дипольные антенны, образующие ортогонально ориентированные наборы для передатчика и приемника или, как вариант, ориентированные параллельно, могут использоваться для передачи и приема электромагнитного излучения. Кроме того, для передачи и приема электромагнитного излучения могут быть использованы одинаково поляризованные или ортогонально поляризованные антенны передатчика и приемника. Эти специальные признаки позволяют использовать широкое разнообразие вариантов осуществления. В способе, соответствующем изобретению,могут быть предусмотрены средства для вращения антенной системы, обеспечивающие вращение антенной системы с помощью механических или электронных средств. Предпочтительно в настоящем изобретении используется электронное вращение для обеспечения управления вращением с высокой точностью при существенном снижении полного веса установки. Электронное вращение может быть обеспечено введением сдвига между отдельными антенными элементами антенной системы, что обеспечивает поляризацию электромагнитного поля относительно геометрического центра с угловыми приращениями в типовом случае порядка 22,5, 30, 36 или долей или целых кратных величин этих углов. В способе, соответствующем настоящему изобретению, антенная система, соответствующая вышеупомянутому варианту выполнения,может использоваться для передачи и приема электромагнитного излучения для генерирования электромагнитного излучения и приема электромагнитного излучения на нескольких отдельных частотах, предпочтительно в диапазоне от 100 МГц до 1 ГГц шагами по 5 МГц. При этом обеспечивается преимущество, состоящее в обеспечении широкого спектра отражений. В способе, соответствующем изобретению,может использоваться средство обработки сигналов для измерения передаточной функции,например передаточной функции по напряжению, передаточной функции по току или их комбинаций, или передаточной функции по мощности для сигналов, соответствующих излучению, переданному антенной системой, по отношению к излучению, принятому антенной системой. В настоящем изобретении предпочтительно измеряется передаточная функция по напряжению для связанных значений изменения угла антенной системы и частоты сигнала, излученного антенной системой, и сигнала, принятого антенной системой. В способе, соответствующем изобретению,средство обработки сигналов может использоваться для выполнения преобразования из час 7 тотной области во временную область с использованием преобразования Фурье или с применением математической экспоненциальной модели с соответствующей рациональной передаточной функцией для передаточной функции, связанной с конкретным углом, для генерирования непрерывной во времени и зависящей от времени функции для каждого угла, который вычисляется в предварительно определенное количество дискретных моментов времени при идентичных моментах вычисления для всех углов. Кроме того, средство обработки сигналов может использоваться для выполнения преобразования из частотной области во временную область путем преобразования Фурье или с использованием математической экспоненциальной модели с соответствующей рациональной передаточной функцией для каждой угловой гармоники, зависящей от времени непрерывной во времени функции в предварительно определенное количество дискретных моментов времени при идентичных моментах вычислений для всех угловых гармоник, для генерирования набора чисел, соответствующих математической функции,имеющей дискретные угловые гармоники и дискретные моменты времени в качестве независимых переменных. Это обеспечивает надежное определение уровня и задержек по времени для основных отражателей даже при отсутствии идеальных условий. В способе, соответствующем настоящему изобретению, средство обработки сигналов может использоваться для выполнения анализа сигналов посредством преобразования Фурье для каждого из используемых моментов, связанных с каждым углом, непрерывной во времени функции, с использованием времени в качестве постоянной и угла в качестве независимой переменной для генерирования представления в области углов. Кроме того, средство для обработки сигналов может выполнять анализ сигналов с помощью преобразования Фурье для каждой частоты, причем частота используется как постоянная, а угол - как независимая переменная для генерирования набора чисел с использованием угловой гармоники и частоты измерений в качестве независимых переменных. Назначением применения преобразования Фурье является выделение содержимого угловых гармоник сигнала при различных временных задержках (и глубинах). В способе, соответствующем настоящему изобретению, средство обработки сигналов может использоваться для выполнения различных операций над обрабатываемыми сигналами. Средство обработки сигналов может использоваться для выполнения сканирования локальных и глобальных максимумов, указывающих периодические по углу отражения от объектов,при представлении в области углов, причем пиковые значения указывают на временную задержку соответствующих отражений от объек 002451 8 та. Кроме того, средство обработки сигналов может выполнять вычисления углового положения относительно антенной системы для периодических по углу отражений, а также средство обработки сигналов может использоваться для применения собранных результатов измерений,выполненных в процессе перемещения по горизонтали, для подавления мешающих отражений. Ниже изобретение описывается более детально со ссылками на чертежи, на которых представлено следующее: фиг. 1-12 относятся к описанию предшествующего уровня техники. Фиг. 1 иллюстрирует, каким образом антенны перемещаются вдоль произвольной траектории по поверхности земли в процессе вращения антенн или путем поворота поляризации посредством постоянного электрического сдвига между антенными элементами. Фиг. 2 иллюстрирует возможную конфигурацию антенн с двумя выводами, символизируемую тонкими электрическими диполями со смещенными центрами. За счет электрического сдвига между 8 антеннами может быть обеспечен поворот поляризации. Фиг. 3 а иллюстрирует одиночную пару коллинеарных антенн с двумя выводами, одна из которых используется в качестве передающей, а другая - в качестве приемной. На фиг. 3b показаны 4 такие пары. На фиг. 3 с показана одиночная пара параллельных антенн с двумя выводами, а на фиг. 3d - 4 такие пары. За счет электрического сдвига между антеннами по фиг. 3b и 3d может быть обеспечен поворот поляризации с приращениями по 45. Фиг. 4 - фиг. 4 а иллюстрирует одиночную пару антенн с двумя выводами, имеющих смещенные центры и поляризованных ортогонально, причем одна из них используется в качестве передающей, а другая - в качестве приемной. На фиг. 4b показаны 8 таких пар. За счет электрического сдвига между антеннами по фиг. 4b может быть обеспечен поворот поляризации с приращениями 45 в процессе измерений с использованием пар, показанных на фиг. 4 а. На фиг. 4 с показана одиночная пара антенн с двумя выводами, имеющих смещенные центры, ортогональные поляризации и асимметричную взаимную связь, причем одна из них используется в качестве передающей, а другая - в качестве приемной. На фиг. 4d показаны 4 такие пары. За счет электрического сдвига между антеннами по фиг 4d может быть обеспечен поворот поляризации приращениями по 45 в процессе измерений с использованием пар антенн, показанных на фиг. 4 с. Фиг. 5 - на фиг. 5 а показан одиночный диполь, который может быть использован как передающая, а также как приемная антенна соответственно измерению s-параметра S11, представляющего измерение отражения (в противоположность измерению передачи, которое используется во всех случаях, когда передающая и 9 приемная антенны разделены). На фиг. 5d показаны 8 таких диполей. За счет электрического сдвига между антеннами по фиг. 5b может быть обеспечен поворот поляризации приращениями по 45 в процессе измерений с использованием пар антенн, показанных на фиг. 5 а. Фиг. 6 - на фиг. 6 а представлены два ортогональных диполя, один из которых используется как передающий, а другой - как приемный. На фиг. 6b показаны две такие пары. Теоретическое применение конфигурации по фиг. 6b заключается в измерении отражения (S11) при использовании только одного диполя в каждый данный момент времени. Фиг. 7 - фиг. 7 а иллюстрирует угол протяженных объектов относительно набора антенн типа 2. Фиг. 7b иллюстрирует угол протяженных объектов относительно набора антенн типа 1. Фиг. 8 иллюстрирует набор антенн в двух положениях относительно объекта. Для ясности в набор включены только 4 диполя. На фиг. 8 а показан вид сбоку набора в двух положениях, а на фиг. 8b - соответствующий вид сверху. Фиг 9 иллюстрирует конфигурации диполей, в которых боковые стороны не параллельны. Такой вариант выполнения одиночного диполя может иметь преимущество в отношении ширины полосы частот диполя. На фиг. 10 представлен схематичный вид возможной конструкции диполей, в которой разъединенные металлические листы 1 электрически соединены с помощью резистивных элементов связи а. На фиг. 10 а показана половина такого диполя. Число элементов может меняться. На фиг. 10b показаны две такие половины с фидерным симметрирующим устройством 2. На фиг. 11 показан схематичный вид возможной конструкции диполей, в которой разъединенные металлические листы 1-5 электрически соединены с помощью резистивных элементов связи а. Боковые стороны не параллельны, и ширина диполей возрастает с увеличением расстояния от центра. На фиг. 11 а показана половина такого диполя. Число элементов может меняться. На фиг. 11b показаны две такие половины с фидерным симметрирующим устройством 2. На фиг. 12 показана фундаментальная конструкция антенны в поперечном сечении по линии, соответствующей диаметру на фиг. 9 а или 9b. Антенные элементы е лежат на опорной плате d из соответствующего диэлектрического материала. Вся конструкция удерживается самой верхней платой а, которая может быть выполнена из любого материала с учетом электрических характеристик. Ниже поддерживающей платы а находится экранирующая пленка, или металлическая плата b размещена так, чтобы предотвратить возможность излучения в направлении вверх; с - поглощающий материал,например ECCOSORB1s или 1s (торговая 10 марка Emerson Cuming) или их комбинации. Ниже активных элементов имеется плата из диэлектрического материала. Антенные элементы е состоят из плоских металлических элементов,соединенных элементами резистивной связи,как показано на фиг. 10 или на фиг. 11, f - диэлектрический материал, который защищает активные элементы механически, а также обеспечивает электрическую настройку. Относительная диэлектрическая проницаемость такого материала в типовом случае находится в пределах от 3 до 10. Фиг. 13-19 относятся к измерениям, осуществленным в соответствии с предпочтительным вариантом осуществления изобретения. На фиг. 13 представлен результат процедуры 1 вплоть до этапа 16 включительно. На фиг. 13 а показан уровень сигнала в функции от индекса n угла и времени (нc) в форме трехмерного графика. На фиг. 13b показан уровень сигнала в функции от индекса n угла и времени (нc) в форме контурного графика. Фиг. 13 с и 13d иллюстрируют результат описываемой процедуры 1 вплоть до этапа 15 включительно. На фиг. 13 с показан уровень измеренного сигнала в функции от индекса n угла и частоты (ГГц) в форме трехмерного графика. На фиг. 13d показан уровень сигнала в функции от индекса n угла и частоты (ГГц) в форме контурного графика. Фиг. 13 е и 13f иллюстрируют результат описываемой процедуры 2 вплоть до этапа 22 включительно. На фиг. 13 е показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде трехмерного графика. На фиг. 13f показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде контурного графика. На фиг. 13g и 13h иллюстрируются кривые обнаружения как результат процедуры 6 вплоть до этапа 6.3 включительно. Фиг. 13g иллюстрирует D в абсолютных величинах, а фиг 13h иллюстрирует D в дБ. На фиг. 13i и 13j иллюстрируются кривые обнаружения как результат процедуры 2 вплоть до этапа 2.3 включительно. Фиг. 13i иллюстрируетD в абсолютных величинах, а фиг. 13j иллюстрирует D в дБ. На фиг. 14 показан результат описываемой процедуры 1 вплоть до этапа 1.6 включительно. На фиг. 14 а показан уровень сигнала в функции от индекса n угла и времени (нc) в форме трехмерного графика. На фиг. 14b показан уровень сигнала в функции от индекса n угла и времени (нс) в форме контурного графика. Фиг 14 с и 14d иллюстрируют результат описываемой процедуры 1 вплоть до этапа 1.5 включительно. На фиг 14 с показан уровень измеренного сигнала в функции от индекса n угла и частоты (ГГц) в форме трехмерного графика. На фиг 14d показан уровень сигнала в функции от индекса n угла и частоты (ГГц) в форме контурного графика. Фиг. 14 е и 14f иллюстрируют 11 результат описываемой процедуры 2 вплоть до этапа 2.2 включительно. На фиг. 14 е показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде трехмерного графика. На фиг. 14f показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде контурного графика. На фиг. 14g и 14h иллюстрируются кривые обнаружения как результат процедуры 6 вплоть до этапа 6.3 включительно. Фиг. 14g иллюстрирует D в абсолютных величинах, а фиг 14h иллюстрирует D в дБ. На фиг. 14i и 14j иллюстрируются кривые обнаружения как результат процедуры 2 вплоть до этапа 2.3 включительно. Фиг. 14i иллюстрирует D в абсолютных величинах, а фиг. 14j иллюстрирует D в дБ. На фиг. 15 представлен результат описываемой процедуры 1 вплоть до этапа 1.6 включительно. На фиг. 15 а показан уровень сигнала в функции от индекса n угла и времени(нc) в форме трехмерного графика. На фиг. 15b показан уровень сигнала в функции от индекса n угла и времени (нc) в форме контурного графика. Фиг. 15 с и 15d иллюстрируют результат описываемой процедуры 1 вплоть до этапа 1.5 включительно. На фиг. 15 с показан уровень измеренного сигнала в функции от индекса n угла и частоты (ГГц) в форме трехмерного графика. На фиг. 15d показан уровень сигнала в функции от индекса n угла и частоты (ГГц) в форме контурного графика. Фиг. 15 е и 15f иллюстрируют результат описываемой процедуры 2 вплоть до этапа 2.2 включительно. На фиг. 15 е показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде трехмерного графика. На фиг. 15f показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде контурного графика. На фиг. 15g и 15h иллюстрируются кривые обнаружения как результат процедуры 6 вплоть до этапа 6.3 включительно. Фиг. 15g иллюстрирует D в абсолютных величинах, а фиг. 15h иллюстрируетD в дБ. На фиг. 15i-15 о иллюстрируются кривые обнаружения как результат процедуры 2 вплоть до этапа 2.3 включительно. Фиг. 15i иллюстрируетD в абсолютных величинах, а фиг. 15j иллюстрирует D в дБ, фиг. 15k и фиг. 15l иллюстрируют результаты вычислений, полученные в описываемой процедуре 3. На фиг. 16 иллюстрируется результат описываемой процедуры 1 вплоть до этапа 1.6 включительно. На фиг. 16 а показан уровень сигнала в функции от индекса n угла и времени(нc) в форме трехмерного графика. На фиг. 16b показан уровень сигнала в функции от индекса n угла и времени (нc) в форме контурного графика Фиг. 16 с и 16d иллюстрируют результат описываемой процедуры 1 вплоть до этапа 1.5 включительно. На фиг. 16 с показан уровень измеренного сигнала в функции от индекса n угла 12 и частоты (ГГц) в форме трехмерного графика. На фиг. 16d показан уровень сигнала в функции от индекса n угла и частоты (ГГц) в форме контурного графика. Фиг. 16 е и 16f иллюстрируют результат описываемой процедуры 2 вплоть до этапа 2.2 включительно. На фиг. 16 е показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде трехмерного графика. На фиг. 16f показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде контурного графика На фиг. 16g и 16h иллюстрируются кривые обнаружения как результат процедуры 6 вплоть до этапа 6.3 включительно. Фиг. 16g иллюстрирует D в абсолютных величинах, а фиг. 16h иллюстрирует D в дБ. На фиг. 16i и 16j иллюстрируются кривые обнаружения как результат процедуры 2 вплоть до этапа 2.3 включительно. Фиг 16i иллюстрируетD в абсолютных величинах, а фиг. 16j иллюстрирует D в дБ. На фиг. 17 представлен результат описываемой процедуры 1 вплоть до этапа 1.6 включительно. На фиг. 17 а показан уровень сигнала в функции от индекса n угла и времени (нс) в форме трехмерного графика. На фиг. 17b показан уровень сигнала в функции от индекса n угла и времени (нс) в форме контурного графика. Фиг. 17 с и 17d иллюстрируют результат описываемой процедуры 1 вплоть до этапа 1.5 включительно. На фиг. 17 с показан уровень измеренного сигнала в функции от индекса n угла и частоты (ГГц) в форме трехмерного графика. На фиг. 17d показан уровень сигнала в функции от индекса n угла и частоты (ГГц) в форме контурного графика. Фиг. 17 е и 17f иллюстрируют результат описываемой процедуры 2 вплоть до этапа 2.2 включительно. На фиг. 17 е показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде трехмерного графика. На фиг. 17f показан вычисленный уровень сигнала в функции от индекса k угловых гармоник и времени (нс) в виде контурного графика. На фиг. 17g и 17h иллюстрируются кривые обнаружения как результат процедуры 6 вплоть до этапа 6.3 включительно. Фиг. 17g иллюстрирует D в абсолютных величинах, а фиг. 17h иллюстрируетD в дБ. На фиг. 17i и 17j иллюстрируются кривые обнаружения как результат процедуры 2 вплоть до этапа 2.3 включительно. Фиг. 17i иллюстрирует D в абсолютных величинах, а фиг. 17j иллюстрирует D в дБ. Фиг. 18 иллюстрирует семейство кривых обнаружения, каждая из которых снята с использованием опытного образца и затем помещена на одном графике, при этом видно, как изменяется кривая обнаружения при прохождении через объект. Объект, соответствующий графикам на фиг. 18, 19 А и 19b, представляет собой трубу из поливинилхлорида, заполнен 13 ную водой, имеющую диаметр 110 мм и находящуюся на глубине примерно 140 см. На фиг. 18 временная задержка преобразована в оценку глубины, причем для скорости распространения радиоволн в земле принята оценка среднего значения. На фиг. 19 а показана типовая кривая обнаружения, полученная для изображения опытного образца. Четкая индикация обнаружения видна на глубине, соответствующей временной задержке примерно 31 нc. На фиг. 19b показана кривая обнаружения,на которой отсутствует индикация обнаруживаемого объекта, такого как кабель или труба. Фиг. 20-44 иллюстрируют представленный вариант осуществления изобретения. На фиг. 20 показана полностью установка,содержащая передающую и приемную антенны,закрепленные на круговом диске, корпус батареи, электронные блоки передатчика и приемника и установку обработки данных, включающую в себя монитор. На фиг. 21 показана блок-схема системы аппаратных средств для опытного образца. Система состоит из трех основных частей: антенной решетки, электронного блока векторного анализатора реального времени и персонального компьютера, снабженного модулем обработки сигналов. На фиг. 22 показана блок-схема электронного блока векторного анализатора реального времени (ВАРВ). На фиг. 23 показана блок-схема электронного блока приемника. На фиг. 24 показана принципиальная схема блока 'а' приемника. На фиг. 25 показана принципиальная схема блока 'b' приемника. На фиг. 26 показана принципиальная схема блока 'с' приемника. На фиг. 27 показана принципиальная схема блока 'd' приемника. На фиг. 28 показана принципиальная схема передатчика. На фиг. 29 показана принципиальная схема блока синтеза. На фиг. 30 показана принципиальная схема блока питания блока синтеза. На фиг. 31 показана принципиальная схема блока непрерывного (быстрого) синтеза. На фиг. 32 показана принципиальная схема блока переменного синтеза. На фиг. 33 показана принципиальная схема опорного блока для блока синтеза. На фиг. 34 показана принципиальная схема блока переключения для блока синтеза. На фиг. 35 показана принципиальная схема схемы обнаружения, содержащей логические цифровые входные/выходные линии для управления блоком синтеза, антенным блоком и блоками поддержки. 14 На фиг. 36 показана принципиальная схема схемы обнаружения, состоящей из аналогоцифровых преобразователей. На фиг. 37 показана принципиальная схема платы антенного переключателя. На фиг. 38 показана принципиальная схема переключателя декодирования/уровня антенного переключателя. На фиг. 39 показана принципиальная схема дипольной антенны. На фиг. 40 показан вид верхней части топологии двусторонней печатной платы схемы дипольной антенны. На фиг. 41 показан вид нижней части топологии двусторонней печатной платы схемы дипольной антенны. На фиг. 42 показан последовательный интерфейс (EZLITE). На фиг. 43-45 показаны альтернативные конфигурации антенных элементов. Фиг. 43 иллюстрирует конфигурацию из четырех пар антенн, все центры диполей смещены относительно геометрического центра. Элементы имеют форму бабочки, причем такая конструкция практически соответствует показанному на фиг. 11, а конфигурация в целом соответствует конфигурации по фиг. 12. Фиг. 44 иллюстрирует конфигурацию из пяти пар антенн, все центры диполей смещены относительно геометрического центра. Элементы имеют форму бабочки, причем такая конструкция практически соответствует показанному на фиг. 11, а конфигурация в целом соответствует конфигурации по фиг. 12. Фиг. 45 иллюстрирует конфигурацию из пяти пар антенн, все центры диполей смещены относительно геометрического центра. Элементы имеют форму бабочки, причем такая конструкция практически соответствует показанному на фиг. 11, а конфигурация в целом соответствует конфигурации по фиг. 12. Для более детального пояснения изобретения представлены следующие таблицы. Таблицы 1-8 содержат информацию относительно представленного варианта осуществления изобретения. Таблица 1 содержит указания выводов и перечень компонентов для электронного блока приемника. Таблица 2 содержит указания выводов и перечень компонентов для электронного блока передатчика. Таблица 3 содержит указания выводов и перечень компонентов для блока синтеза. Таблица 4 содержит указания выводов и перечень компонентов для блока обнаружения. Таблица 5 содержит указания выводов и перечень компонентов для блока антенного переключателя. Таблица 6 содержит указания выводов и перечень компонентов для переключателя декодирования/уровня антенного переключателя. 15 Таблица 7 содержит перечень компонентов для последовательных интерфейсов по фиг. 4244. Таблица 8 содержит перечень компонентов для элемента дипольной антенны. Таблицы 9-11 содержат информацию, касающуюся возможных вариантов осуществления изобретения. Таблица 9 содержит данные для ряда возможных антенных конфигураций. Таблица 10 содержит примеры числа поворотов антенной системы, причем- угол поворота между каждым измерением, N - число измерений. Таблица 11 содержит примеры индексов k,соответствующих двум периодам на поворот. Значения индексов k даны парами (k2, k-2),угол поворота между каждым измерением, N число измерений. Таблица 12 содержит информацию, касающуюся измерений согласно представленному варианту осуществления изобретения. В таблице 12 представлены итоговые данные для объектов, проиллюстрированных в примерах 1-5. На фиг. 20 представлен пример осуществления 10 настоящего изобретения, который содержит конструкцию, выполненную для обеспечения перемещения, включающую в себя набор колес 32, опору 26, U-образную ручку 30, на открытом конце которой закреплена опора 26,причем колеса 32 закреплены на валу, закрепляемом с помощью ряда выступов U-образной ручки 30 на открытом ее конце. Корпус 12 установлен на платформе 26, содержащей индивидуальные антенные элементы, выполненные как описано ниже. Индивидуальные антенные элементы закреплены на верхней стороне корпуса 12 с помощью фиксаторов 18 со струбцинами. Металлические винты 24 обеспечивают дополнительное крепление индивидуальных антенных элементов. Болты 22 из пластика скрепляют вместе верхнюю пластину и нижнюю пластину корпуса 12. Передаваемые сигналы подаются, и принимаемые сигналы принимаются, как описано ниже, соответственно на индивидуальные антенные элементы и от них через проводники 20, соединенные с блоками 16 переключателя,содержащими, совместно с электронным логическим блоком 14, средства управления переключением передачи-приема индивидуальных антенных элементов. Электронный логический блок 14 размещен под блоками 16 переключателя. Передача сигналов между индивидуальными антенными элементами и электронным блоком 34 установки, включающим в себя электронный блок синтеза, электронный блок измерения частоты и электронный блок передатчикаприемника, как описано ниже, обеспечивается с помощью набора проводников 28. Средство обработки сигналов содержит компьютер 38, свя 002451 16 занный с монитором 42 и клавиатурой 40. Для питания компьютера 38, монитора 42, клавиатуры 40 и электронного блока 34 установки предусмотрен блок питания 36, расположенный на противоположной стороне ручки 30 по отношению к компьютеру 38. Ниже представлено подробное описание изобретения. 1. Предпосылки изобретения Настоящее изобретение относится к способам и устройствам для обнаружения объектов в земле с помощью соответствующих радиолокаторов. Такие радиолокаторы, облучающие землю или иные материалы, давно используются и известны как георадары, радиолокаторы глубинного обнаружения, подповерхностные радиолокаторы. В течение длительного времени пытались использовать такие подповерхностные радиолокаторы для поиска объектов в земле, в том числе кабелей, труб и оптических волокон. Однако известные из предшествующего уровня техники радиолокаторы обнаружения объектов под землей, по существу, не пригодны для этой цели, что может считаться причиной того, почему отсутствовала коммерческая технология радиолокационных обнаружителей кабелей и труб, несмотря на настоятельную потребность и множество экспериментов. При благоприятных условиях коммерчески доступные радиолокационные системы обнаружения объектов под землей обеспечивают обнаружение кабелей и труб,но только в случае, когда радиолокационная антенна перемещается вперед по одной или более прямым линиям, поперечно обнаруживаемым объектам, и только путем интерпретации зарегистрированных измерений, которые требуют высокого уровня квалификации интерпретатора. Радиоволны имеют поляризацию, связанную с пространственной ориентацией соответствующих электромагнитных полей. Использование поляризации добавляет еще одну степень измерений радиолокаторам обнаружения объектов под землей, в частности, таким образом, что объекты могут быть селектированы на окружающем их фоне, в дополнение к традиционному разрешению по времени. Использование поляризации в качестве средства различения объектов, особенно протяженных объектов, не привело, однако, в ранее использовавшихся радиолокаторах к возможности обнаружения кабелей и труб. Действительно, известны патенты и публикации, посвященные поляризации, в частности работа авторовYoung et al. из Университета штата Огайо (см. пример [8]). В этом направлении работали и другие исследователи, включая группу Guntonet al. в Великобритании (см. примеры [1], [2],[10]). Кроме того, был проведен ряд испытаний по переносу методов, известных как синтезированная апертура или синтезированная антенная решетка, используемых в качестве бортовых и спутниковых радиолокационных систем, на ра 17 диолокаторы обнаружения объектов под землей. Применение поляризации также известно из этих систем, что часто определяется как поляриметрия. Общей особенностью методов поляриметрии и синтезированной апертуры является то, что они требуют довольно однородной среды, что облучаемые объекты находятся в дальней зоне антенны и что на взаимную фазу радиолокационных сигналов не оказывают существенного влияния неоднородности в среде распространения. Когда антенна находится вблизи земли или в условиях распространения волны в почве,почти всегда будет иметься сильное влияние на фазу и амплитуду волн со стороны неоднородностей в форме камней, металлических объектов, неровностей поверхности, влаги и т.д. Влияние таких неоднородностей может привести к значительным колебаниям при перемещении антенны на несколько сантиметров. Методы использования поляризации, как описано в [1],[2], оказываются неэффективными в условиях таких сигналов, как указано выше. Периодичность, присущая радиолокационным сигналам от тонких удлиненных объектов в земле при идеальных условиях, как описано в [1], [2], на практике оказывается маскированной случайными изменениями значительного уровня. Шум,обусловленный такими неоднородностями, делает на практике ненадежным обнаружение на основе простых алгебраических комбинаций измеренных радиолокационных сигналов. Вращение поляризации является полезным и эффективным средством для улучшения характеристик радиолокаторов обнаружения объектов под землей. В [1, Gunton et al.] представлена следующая математическая модель зависимости угла поляризации для длинных тонких объектов:F(t)=A(t)+B(t)sin(X)+C(t)cos(X)+D(t)sin(2X)+E(t)cos(2X),где Х - угол поворота антенной системы. Не совсем ясно, как фазы A(t)E(t) связаны в данном уравнении. В любом случае можно заключить, что математические выражения, подобные приведенному выше, которые могут рассматриваться как первая часть в ряду Фурье, не могут использоваться на практике для описания сигналов, измеренных в радиолокаторе обнаружения объектов под землей, с каким-либо полезным приближением. Изобретателями было найдено, что надлежащий анализ Фурье угловой зависимости совместно с другими математическими операциями является весьма полезным инструментом для обнаружения объектов. Проанализировав ортогональные антенны с двумя выводами с совпадающими электрическими центрами, использованные в патентах [8],[1], [2], изобретатели пришли к выводу, что они не являются наиболее пригодными для обнаружения объектов, подобных кабелям или трубам. Аргументом, обосновывающим использование взаимно перпендикулярных антенн с двумя вы 002451 18 водами с совмещенными центрами (т.е. скрещенных диполей), согласно описаниям указанных патентов, является то, что взаимная электромагнитная связь между ними равна теоретически нулю и что, таким образом, отсутствует прямая связь между сигналами от передатчика к приемнику (определяемая понятием проникающий сигнал). Изобретателями было найдено, что антенны со смещенными центрами предпочтительны в большинстве случаев. Причина этого состоит в том, что симметрия, которая необходима в конструкции со скрещенными диполями, на практике подвергается возмущениям со стороны окружающей среды, в частности неровностей поверхности земли. Поскольку ближние зоны расположенных рядом антенн в большой степени совпадают, малые отклонения от симметрии вносят большие искажения, в результате чего возникает существенная взаимная связь. Поэтому разнесение центров антенн является более предпочтительным. Предложенная предпочтительная антенная конфигурация представляет собой коллинеарные или параллельные диполи со смещенными центрами, в противоположность описаниями патентов [1], [8], где антенны, а также способы обнаружения основывались на ортогональности. Однако следует подчеркнуть, что предложенный способ обнаружения, включая обработку сигнала, чрезвычайно полезен и при использовании ортогональных антенн. Предложенная предпочтительная форма сигнала соответствует так называемому радиолокатору с непрерывным сигналом со ступенчато изменяющейся частотой, в котором излучаются сигналы синусоидальной волны, имеющие частоту, изменяющуюся ступенчато. Термин непрерывный сигнал означает, что синусоидальное колебание излучается на протяжении такого времени, что в процессе каждого измерения переходные процессы могут считаться закончившимися. В радиолокаторе обнаружения подземных объектов, предлагаемом в настоящем изобретении, излучаются частоты в интервале от 100 до 1000 МГц с приращениями частоты 5 МГц (однако способ, соответствующий изобретению, может быть использован и в любом другом диапазоне частот). Данная форма сигнала требует способа преобразования во временную область. Изобретателями было установлено, что электромагнитное затухание в земле является,по существу, экспоненциальной функцией частоты в соответствующем частотном диапазоне. На основе этого была создана экспоненциальная модель сигнала, предназначенная для использования в качестве основы для оценки временных задержек радиолокационного сигнала. Способ несколько сходен со способами Рrоnу (см. [7]). Способы, полученные из исходного способа Рrоnу, уже известны в анализе радиолокационных сигналов. Новым признаком настоящего изобретения является то, что математически 19 когерентным способом он учитывает дисперсионные характеристики материалов в земле. В принципе, этим способом можно определить интенсивность отражателя, независимо от того,что отражатель погружен в материал с сильным ослаблением. Кроме того, он обеспечивает лучшее временное разрешение, чем преобразование Фурье, которое часто выполняется с использованием алгоритмов быстрого преобразования Фурье (БПФ). Классическая взаимосвязь между шириной полосы и разрешением по времени,согласно которой наименьшая обнаруживаемая разность по времени обратно пропорциональна ширине полосы сигнала, применима, как известно, когда используются способы на основе преобразования Фурье. При определенных обстоятельствах лучшее разрешение может быть получено способом Рrony. Изобретателями было найдено, что обязательным условием удовлетворительного обнаружения объектов является то, что определяемая поляризацией угловая зависимость должна исследоваться отдельно при различных временных задержках или, иными словами, отражатели, находящиеся на различных расстояниях от антенн по отношению ко времени, анализируются по поляризации отдельно. Следовательно,необходим двумерный анализ сигнала, при котором упомянутые две размерности соответствуют углам и времени. Кроме того, изобретателями было обнаружено, что перемещение антенн в процессе вращения имеет полезный эффект нарушения периодичности в угловой зависимости очень близко расположенных объектов, в то время как периодическая зависимость более удаленных (и,вероятно, более крупных) объектов сохраняется. Поэтому анализ Фурье по углу поворота для сигнала, измеренного при таких условиях, покажет, что сигналы от объектов вблизи поверхности будут иметь энергию, распределенную в пределах широкого спектра гармоник, в то время как, например, кабели и трубы обусловят значительно более концентрированный спектр в окрестности гармоник, соответствующих 2 периодам на вращение. Таким образом, достигается ослабление нежелательных отражений от объектов и неровностей поверхности. Это очень хорошо соответствует практическому использованию радиолокатора обнаружения подземных объектов для сканирования, осуществляемого в направлении вперед вручную или при установке на транспортное средство или машину и использовании механического или электрического вращения поляризации. Анализ сигналов и, следовательно, способ обнаружения, как описано в [1] и [2], основываются на сравнении с опорными измерениями,в [2] в форме согласованных фильтров, что эквивалентно корреляции с известными кривыми. Проблема с использованием опорных измерений состоит в том, что не существует двух 20 объектов с соответствующей окружающей средой, дающих одинаковое отражение. Поэтому способ, основанный на опорных измерениях,трудно использовать практически. Настоящее изобретение направлено на решение этих проблем за счет использования существующих способов и устройств. Как очевидно из некоторых примеров, иллюстрирующих обработанные данные измерений, проведенных при реальных условиях, отношение сигнал/шум 20 дБ и более при обычных условиях может быть получено с помощью способов, как описано ниже. Для сравнения,путем сканирования известным радиолокатором обнаружения подземных объектов в поперечном направлении относительно кабеля или трубы может быть получен эхо-сигнал от этих объектов примерно того же уровня, что и от различных объектов окружающей среды, т.е. отношение сигнал/шум составляет примерно 0 дБ при тех же самых условиях. Использование способа Рrоnу для преобразования из частотной области во временную область вместе с анализом Фурье для углов и вращением поляризации обеспечивает надежность обнаружения, которая не могла быть достигнута ранее. 2. Введение Настоящее изобретение продемонстрировано в варианте осуществления, назначением которого является обнаружение и определение ориентации удлиненных объектов в земле, на практике представляющих собой трубы из различных материалов, кабели, содержащие металлические проводники и оптические волокна. Однако описываемый способ, включая математическое описание и применение электромагнитных волн с вращающейся поляризацией,может быть использован и для других целей,например для обнаружения противопехотных мин. Описываемый способ, что касается определения направления, т.е. ориентации, ограничен, однако, удлиненными объектами. В последующем описании это имеется в виду при упоминании об обнаружении объектов, находящихся в земле. Настоящее изобретение также пригодно для обнаружения объектов, скрытых в других материалах, например в бетоне, в других строительных материалах, в скальных породах,что может рассматриваться как материалы земли по отношению к электромагнитным волнам. Кроме того, объекты в воздушной среде также могут обнаруживаться описываемым способом. Изобретение основывается на использовании электромагнитных волн, которые передаются в материал среды, в котором скрыты один или более объектов. Таким образом, способ основан на принципе радиолокации. Электромагнитные (ЭМ) волны могут, в соответствии с общепринятым определением, называться радиоволнами. Изобретение основывается на излучении ЭМ волн с помощью антенны, измерении, 21 регистрации и последующей обработке электромагнитных откликов или отражений, обусловленных объектом. Радиолокатор, как известно, например, в условиях аэропортов или некоторых иных мест использования, зачастую располагается на относительно большом расстоянии от отражающих объектов по сравнению с длиной используемых ЭМ волн и/или длительности излучаемых импульсов. В непосредственной близости от антенны радиолокатора часто находятся объекты, поверхность земли, что в соответствии с общепринятыми терминами определяется как ближняя зона антенны. Электромагнитное влияние таких ближних объектов часто не толкуется как отражения, поскольку превалируют так называемые реактивные поля. В последующем изложении термин отражения используется, без детальных различий, применительно ко всем зарегистрированным откликам, обусловленным объектами вне антенны, как ближними, так и удаленными. Настоящее изобретение основывается на использовании комбинации временного разрешения электромагнитных откликов и их поляризационной зависимости. При этом не важно,вызваны ли временные задержки и поляризационная зависимость от элементов, определенных понятиями передача, отражение, ближние зоны,дальние зоны, индуктивная связь, емкостная связь или реальные радиолокационные отражения. Предпочтительный вариант осуществления изобретения выполнен так, что объекты,подлежащие обнаружению, находятся вне ближней зоны антенны, но существенным преимуществом изобретения является то, что оно обеспечивает подавление эффектов ближней зоны, даже если они характеризуются более высокими уровнями, чем полезные отражения. В радиолокационной системе, соответствующей изобретению, одна антенна используется для излучения ЭМ волн, а другая антенна используется для приема откликов на излучения, или используется множество таких антенн. В возможном варианте можно осуществлять выбор одной или множества антенн для передачи и приема, возможно, в одно и то же время. В последующем передаточная функция между передатчиком и приемником используется независимо от конструкции антенны. Важным моментом является то, что в системе проводится различие между излученным сигналом и принятым сигналом и производится их сравнение. Сравнение осуществляется косвенным образом на основании того, что излучаемый сигнал, или непосредственно переменное напряжение, определяется и является известным, а принимаемый сигнал измеряется и регистрируется. В последующем описании сигнал означает электромагнитный отклик или связь между из 002451 22 лученным и принятым сигналом. В последующем описании объекты, помимо хорошо определенных объектов, подобных трубам, кабелям и т.д., также означают различия в характеристиках материала, т.е. неоднородности в среде, в которой распространяются ЭМ волны. Под вращением антенн понимается вращение поляризации с помощью физического (механического) вращения, а также электрических средств,включая электрический сдвиг между рядом антенных элементов. Измерения означают осуществляемое вручную или автоматически определение физических размеров с помощью соответствующей аппаратуры. Регистрация и запоминание могут осуществляться также вручную или автоматически. Изобретение использует анализ сигналов электромагнитного отклика для определения поляризационной зависимости и временных задержек (или одного из двух) сигналов электромагнитных откликов. Анализ включает вычисления на базе измеренных величин. На практике возможные варианты осуществления настоящего изобретения наилучшим образом реализуются с помощью цифровой аппаратуры и некоторых физических элементов, например частоты излучаемых ЭМ волн, изменяемой с определенными интервалами. Поэтому множество чисел будет использоваться в последующем для представления физических элементов,которые по своей природе могут быть непрерывными. В последующем упоминаются двумерные множества чисел, которые являются концентрациями дискретных значений. Эти наборы чисел могут толковаться как элементы матрицы, а каждая строка и каждый столбец в ней могут толковаться как последовательности чисел. В качестве их синонимов используются кривые или функции или значения функций. Возможные варианты осуществления изобретения базируются на дискретных частотах, углах,угловых гармониках и моментах времени, что обосновано использованием цифровых блоков памяти и вычислительных блоков, хотя представленные примеры продемонстрированы как непрерывные изображения. Принципы изобретения применимы к непрерывным частотам, угловым поворотам, угловым гармоникам, а также к соответствующим дискретным признакам. Точное воплощение математических формул, приведенных ниже, не является единственно возможным. Например,существуют формулы для непрерывного и дискретного преобразования Фурье, но в технике хорошо известно, что преобразования Фурье и последовательности Фурье могут быть вычислены различными путями. Приведенные формулы использованы для того, чтобы продемонстрировать возможные варианты осуществления новых принципов изобретения на практике. Показанные примеры кривых и вычислений, основанных на измерениях физических элементов, 23 являются всего лишь возможными иллюстрациями принципов изобретения. Числа, указывающие ссылочные источники, приведены в квадратных скобках []. Номера уравнений, представленных в описании, показаны в круглых скобках . 3. Вращение поляризации Если в процессе распространения электромагнитная волна падает на объект, электрические характеристики которого отличаются от электрических характеристик окружающей среды, то часть энергии волны отклоняется или отражается при падении на объект. Хорошо известным фактом является то, что при определенных обстоятельствах это отражение (рассеяние) зависит от поляризации волны по отношению к пространственной ориентации объекта. Под поляризацией электромагнитной волны понимается направление вектора электрического поля. Эта поляризация может быть переменной, диффузной или может быть примерно линейной в заданном направлении, или может быть круговой, что означает, что она вращается вокруг оси. В возможном варианте осуществления настоящего изобретения, показанном ниже, используется линейная поляризация. Кроме того,другие типы поляризации, включая круговую и эллиптическую, могут генерироваться физически и могут анализироваться как линейные комбинации линейно поляризованных составляющих, поэтому линейная поляризация рассмотрена ниже отдельно. 3.1. Линейная поляризация Под линейно поляризованными волнами понимаются электромагнитные волны, которые сформированы как линейно поляризованные в горизонтальном (наиболее часто) направлении настолько, насколько это возможно практически. В способе предусматривается, что приближенно линейная горизонтальная (наиболее часто) поляризация волн вращается относительно перпендикулярной (наиболее часто) оси. Возможный вариант способа использует волны,которые в возможно большей степени должны характеризоваться следующим:- линейной поляризацией в различных направлениях, параллельной уровню, параллельному продольной оси удлиненного объекта(объектов), которые обнаруживаются и местоположение которых определяется,- вращением направления линейной поляризации относительно оси, перпендикулярной упомянутому уровню (и перпендикулярной продольной оси объекта), но не обязательно пересекающей продольную ось объекта. Поскольку обнаружение и определение местоположения кабелей. оптических волокон или труб, погруженных в землю, является основным назначением данного изобретения, направления линейной поляризации зачастую будут горизонтальными, а вращение будет произ 002451 24 водиться относительно вертикальной оси. То,что в данном контексте такие электромагнитные волны не описываются как волны с круговой поляризацией, объясняется тем, что в данном способе нет необходимости, чтобы линейная поляризация была на всех возможных углах или чтобы вращение было непрерывным. Дополнительно причиной этого является то, что в изобретении учитывается, что на практике очень трудно получить хорошо определенную круговую или эллиптическую поляризацию при распространении в земле ввиду весьма существенно неоднородного состава почвенной среды. Электромагнитные волны, упоминаемые в настоящем описании, генерируются приложением высокочастотного электрического напряжения к антенне с соответствующей конструкцией. Для получения ЭМ волны с практически линейной поляризацией часто используются электрические диполи или двухэлементные антенны. Затем один или более диполей располагаются, например, над поверхностью земли с горизонтальной ориентацией их продольной оси(см. фиг. 2, 3 и 4). В последующем описании для удобства используются двухэлементные дипольные антенны, однако и другие типы антенн могут также использоваться, например так называемые щелевые антенны. Для обнаружения кабелей и труб в земле часто используют горизонтально ориентированные дипольные антенны, поскольку кабели и трубы обычно расположены в земле горизонтально. Не является обязательным условием в изобретении, что объекты располагаются в земле точно горизонтально, и также не является необходимым, чтобы положение антенн было строго горизонтальным над землей. Хорошо известно, что если на удлиненный объект, например кабель или трубу, в земле падает электромагнитная волна, поляризация которой практически горизонтальна, то уровень отраженного от объекта сигнала будет периодической функцией угламежду направлением поляризации волны и продольным направлением объекта. Это описано в [2] и в других источниках, а также в [1]. Наконец, также хорошо известно, что если на линейно поляризованную антенну, например диполь, падает линейно поляризованная ЭМ волна, то величина напряжения, индуцированного на выводах антенны, пропорциональна косинусу угла между направлениями поляризации антенны и волны соответственно. Эти хорошо известные факты использованы в настоящем изобретении для обеспечения возможности обнаружения, определения местоположения и направления удлиненных объектов в земле без использования опорных измерений,подобных описанным в [1]. Возможный вариант осуществления предложенного способа с использованием линейной 25 поляризации не требует применения полностью линейно поляризованных ЭМ волн или применения антенн, предназначенных для излучения ЭМ волн, излучающих такие полностью горизонтально поляризованные волны, или применения антенн, предназначенных для приема отражений, чувствительных к волнам, имеющим определенную поляризацию. Также не является необходимым, чтобы обнаруживаемые объекты отражали только волны с определенной поляризацией. Для применения заявленного способа достаточно, чтобы эти условия выполнялись в некоторой степени. 3.2. Другие формы поляризации Круговая и эллиптическая поляризации могут быть сформированы физически и могут как таковые рассматриваться как линейные комбинации линейных поляризаций. Например,круговая поляризация может быть получена путем приложения напряжения переменного тока к ортогонально линейно поляризованным антеннам со взаимным смещением по фазе 90. 4.1. Антенные конфигурации для линейной поляризации Линейная вращающаяся поляризация применена в возможном варианте осуществления настоящего изобретения, который продемонстрирован изобретателями. Ниже некоторые возможные конфигурации антенн с линейной поляризацией будут описаны схематично. Антенны определяют поляризацию, которая должна быть использована. Представляет важность угол между поляризацией передающей и приемной антенн. Последующее описание основано на двух основных конфигурациях для линейной поляризации: измерение 1 типа: совпадающая поляризация: передающая и приемная антенны с одинаковой поляризацией; измерение 2 типа: кросс-поляризация: поляризации передающей и приемной антенн взаимно перпендикулярны. На фиг. 3 и 4 иллюстрируются основные конструкции возможных комбинаций антенн в форме очертаний комбинаций дипольных антенн типа 1 и типа 2 соответственно. Ниже ссылки будут даваться на Sпараметры, поскольку они являются хорошо определенными. S21 характеризует измерение прохождения между двумя наборами выводов(даже хотя отражения в среде передачи являются важными факторами), и S11 характеризует измерение отражения на одиночном соединении. Использование S-параметров является символическим, также могут использоваться и другие определения. В таблице 9 приведены характеристики некоторых возможных конфигураций антенн. Все комбинации на фиг. 3 и 4 составлены из пар диполей: один диполь на передачу, другой диполь на прием. На фиг. 5 иллюстрируется измерение с одиночной антенной, что означает 26 измерение отражения, которое может представлять параметр S11, так как излученный и принятый сигналы существуют на одних и тех же выводах. Конфигурация, показанная на фиг. 6 а и 6b,соответствует использованной в [1] и [8], причем на практике эта конфигурация оказалась менее полезной для обнаружения объектов в земле, в числе прочего вследствие того, что прямая электромагнитная связь между передающей и приемной антеннами имеет величину,сравнимую с ЭМ волнами, отраженными от объектов в земле. Электрические центры антенн совпадают на фиг. 6 а и 6b. Несмотря на практические ограничения, такая конфигурация пригодна для приближенного теоретического анализа соотношений, связанных с поляризационно-чувствительным обнаружением объектов в земле. Если, как показано на фиг. 7 а, два скрещенных диполя размещены так, что их взаимный центр находится на оси, пересекающей перпендикулярно удлиненный тонкий объект и диполи, лежащие в плоскости, параллельной упомянутому линейному объекту, то применимы приведенные ниже математические соотношения для амплитуды электрического напряжения переменного тока на выводах приемной антенны, причем угол между одним из диполей и объектом равен . Предполагаем, что высокочастотное напряжение переменного тока прикладывается к передающей антенне с постоянной амплитудой,так что прикладываемое напряжение может быть представлено математическим уравнениемVtx(t)=Atxcos(2f),(1) где Vtх(t) - прикладываемое зависящее от времени напряжение,t - время,Atх и f - амплитуда и частота приложенного сигнала соответственно. При условии, что этот электрический сигнал приложен к выводам передающей антенны,для напряжения на выводах приемной антенны справедливо следующее приближенное уравнение:Vr1(t, )=A1sin(2)Atxcos(2ft+1),(2) где A1 - постоянный амплитудный коэффициент, определяемый физическими условиями,1 - начальная фаза. Если два коллинеарных диполя расположены так, как показано на фиг. 7b, и длина диполей и расстояние между их центрами достаточно малы в сравнении с расстоянием между диполями и объектом, применимо следующее соотношение:Vr2(t, )=А 2(1+cos(2Atxcos(2ft+2), (3) где А 2 и 2 - постоянные. Следует отметить, что напряжение приемной антенны зависит от времени t и от угла поворотаантенны относительно объекта и что 27 изменения, обусловленные изменением ,обычно происходят намного медленнее, чем изменения по времени. Аналогичные математические уравнения применимы к типам антенн, иных, чем электрические диполи, если только они имеют характеристику линейной поляризации. Два математических соотношения (2) и (3) основаны на ряде допущений:- что отсутствует отражение от приемной антенны к передающей антенне,- что антенны, а также объекты весьма тонки по сравнению с длиной используемых ЭМ волн, что является условием для антенн и для объектов, рассматриваемых как линейно поляризованные,- что толщина антенн и объектов мала относительно расстояния между антеннами и объектом,- что вся среда, окружающая антенны и объект, является однородной во всех направлениях, что также может быть выражено как то,что среда, в которой распространяются электромагнитные поля и волны, является однородной. Настоящее изобретение обеспечивает возможность надежного обнаружения кабелей и труб в земле за счет использования поляризационного анализа без выполнения вышеуказанных идеальных предварительных условий. 4.2. Уравнения для угловой зависимости при линейной поляризации В последующем особое внимание будет уделено передаточной функции по напряжению между передатчиком и приемником, поскольку она содержит важную информацию, а не напряжению на выводах приемной антенны. Если удалить зависимые от времени части из уравнений (2) и (3), чтобы подчеркнуть важность углового поворота , получим следующее уравнение для угловой зависимости передаточной функции по напряжению:(4) для угловой зависимости в случае антенной конфигурации на кросс-поляризации (типа 2) иH2=A2(1+cos(2 для угловой зависимости в случае антенной конфигурации на совпадающей поляризации(типа 1). Однако вышеприведенное математическое уравнение с упомянутыми предварительными условиями никогда не применяется точно. На практике, при измерении соотношения между напряжением переменного тока, приложенным к передающей антенне, и имеющимся напряжением на приемной антенне возникает более сложная ситуация. 4.3. Уравнения для угловой зависимости при круговой поляризации Волны с круговой поляризацией могут быть получены с помощью линейной комбина 002451 28 ции двух линейно поляризованных волн. Антенна с круговой поляризацией может быть создана путем комбинации линейно поляризованных антенн. Например, могут быть использованы схематичные конфигурации, в которых каждая из показанных комбинаций удваивается, так что одна включается в передающую антенну, а другая - в приемную антенну. Если передающая и приемная антенны имеют круговую поляризацию и если, как упомянуто, предполагаем, что объект является тонким и удлиненным и, следовательно, является чувствительным к линейной поляризации, результатом является передаточная функция по напряжению, которая не изменяется по величине. При идеальных условиях, соответствующих упомянутым выше условиям, получим передаточную функцию вида(5a2) При определенных условиях отклик двух антенн с круговой поляризацией на линейный объект будет постоянным, т.е. Нc=Аc. Здесь j - оператор мнимого представления в комплексных числах,j2=-1,Аc - амплитудный коэффициент,Нc - теперь комплексная величина и выражает (идеальным образом) зависимый от угла вращения поворот фазы в передаточной функции по напряжению. Поскольку поляризация излученной волны при круговой (эллиптической) поляризации вращается непрерывно в процессе распространения,отсутствует, подобно случаю линейной поляризации, в уравнении для поляризационного состояния, которое сохраняется, пока антенная система находится в заданном состоянии или положении. В противоположность этому,присутствует в уравнении (5 а 1) и (5a2) для начального поворота поляризации волны. Такой начальный поворот может изменяться соответствующими средствами, подобно вышеупомянутой линейной поляризации. Угловая зависимость, как в уравнениях(5 а 1) и (5a2), может быть получена комбинациями ортогональных антенн со взаимным поворотом фазы на 90. Следует заметить, что нет необходимости в физическом существовании такой схемы поворота фазы или в том, чтобы электрические сигналы имелись в ортогональных антеннах в одно и то же время. Линейные комбинации линейных поляризаций, включая круговую, могут быть образованы путем возбуждения наборов, например диполей, в произвольном порядке, запоминания измеренных сигналов и затем суммирования их алгебраически. Поворот фазы сигналов получается путем 29 умножения на соответствующие комплексные постоянные перед суммированием. Если круговая или эллиптическая поляризация получена с помощью комбинации линейно поляризованных антенн, начальный поворот полученной в результате (приближенно) круговой поляризации отслеживает механический или электрический (электронный) поворот использованных линейно поляризованных элементов. Если используется излучение волн с эллиптической или круговой поляризацией, то измеренные сигналы могут быть использованы для обнаружения объектов с применением способов, раскрытых в патенте на использование линейной поляризации. Практическое использование круговой или эллиптической поляризации может быть основано на разложении измеренной комплексной угловой зависимости на действительную и мнимую части, соответственно тому, что часто называют квадратурным разложением на синфазную (I) и квадратурную (Q) составляющие. Например, уравнение (5 а 2) может быть записано в виде Нc=Acexp(-j2)=Accos(2)-jAcsin(2),откуда видно, что мнимая часть соответствует уравнению для линейной поляризации (4). 5. Анализ Фурье и измерения,основанные на вращении Представленные выше выражения H1 иH2 для угловой зависимости передаточной функции (при линейной поляризации) характеризуются поворотом за два периода, когда антенная система делает один поворот. Если разместить любую комбинацию передающих и приемных антенн над землей и поворачивать ее вокруг вертикальной оси, то в процессе вращения будет иметь место существенное изменение передаточной функции, особенно если антенна размещена вблизи земли, ввиду неоднородностей в земле, таких как камни, меняющиеся виды почв, влажности, материалов отходов и т.д. Кроме того, неровности поверхности земли часто оказывают заметное влияние. Помимо изменений по величине передаточной функции, будут иметь место изменения фазы более сложным образом, чем сдвиг знака полярности,имеющийся в вышеприведенных выражениях для H1 и H2. Разность фаз между напряжением переменного тока, приложенным к передающей антенне, и напряжением, измеренным на приемной антенне, является фазой передаточной функции. Если повернуть антенну точно на один оборот, т.е. на 360, то можно точно измерить одну и ту же передаточную функцию перед поворотом и после него, поэтому очевидно, что существует периодическая угловая зависимость. Но помимо функций, испытывающих изменение за два периода на один поворот, часто исследуются зависимости отс более высо 002451 30 кой периодичностью и зависимость, соответствующая одному периоду на поворот. В настоящем изобретении было обнаружено, что угловая зависимость передаточной функции по напряжению между двумя полностью или частично линейно поляризованными антеннами, соответственно излучающими и принимающими электромагнитные поля или волны, распространяющиеся в землю и полностью или частично отражающиеся от объектов в земле или на земле, с хорошей степенью приближения может быть описана посредством ряда Фурье и что это может быть использовано для существенного улучшения обнаружения удлиненных объектов по сравнению с известными способами. Такой ряд Фурье может быть выражен математически следующим образом: где n=0, 1, 2, N/2, что представляет собой ряд Фурье из N+1 элементов. Сo представляет среднее значение, которое постоянно в процессе вращения. Постоянные Фурье Сn в уравнении(6) являются комплексными и могут быть вычислены в форме интегралов с помощью следующего уравнения: где n=0, 1, 2, N/2.- поворот антенной системы, основанный на произвольном начальном положении, - переменная, обозначающая переменную угловую гармонику,j2=-1,j - оператор мнимого представления в комплексных числах. Уравнение (6) описывает изменение передаточной функции по напряжению в процессе поворота антенной системы как функцию угла поворота . Введенная переменнаяаналогична частотной переменной , которую часто используют для представления (угловой) частоты при Фурье-преобразовании временных гармонических функций. Заметим, что в (6) имеется конечное число элементов; на практике обнаружилось, что лишь малое число гармоник необходимо для эффективного представления Н. Н - комплексная функция, содержащая амплитудную и фазовую информацию передаточной функции. Математический способ записи комплексных чисел, используемый в данной связи, соответствует способу записи, используемому инженерами в области электронной техники и иными специалистами для представления периодически изменяющейся величины, и описан в ряде источников, например [3], [4]. Примечания относительно постоянных Фурье: 31 Со в вышеприведенном уравнении (6) описывает среднее значение, которое в типовом случае обусловлено отражениями в земле, не зависящими от поворота поляризации. С 2 и С-2 - две постоянные, соответствующие двум периодам вращения антенной системы. При линейно поляризованных антеннах C2 и С-2 при идеальных условиях являются взаимно комплексно сопряженными, т.е. С 2=С-2 (символозначает комплексное сопряжение). В случае антенн круговой поляризации С 2 или С-2 при идеальных условиях равно нулю.H - по своей сущности является непрерывной функцией углаи до сих пор рассматривалась как таковая, но зачастую большое практическое значение будет иметь измерение Н только при малом числе конечных значений , например 8 или 16 значений в течение одного поворота, что соответствует повороту антенной системы на 45 или на 22,5 соответственно между каждым измерением. Если Н измеряется, как описано, только при некоторых дискретных значениях , угол может быть записан математически в следующем виде:=n,где n=0N-1, - наименьший угловой поворот, например 45 или 22,5. Основываясь на таких измерениях Н,комплексные компоненты Фурье могут быть вычислены по формуле где Ck - периодическая функция по k. Ck образует дискретную функцию, независимая переменная которой может рассматриваться как квантованная копия вышеупомянутой угловой гармонической переменной , но для практического использования для обнаружения объектов полезно рассматривать ее как функцию по индексуN может равняться 8, 16, 32 или 64. Число измерений в (8) равно числу точек вычислений для функции Vf(k). Это уравнение (8) может быть эффективно вычислено посредством быстрого преобразования Фурье (БПФ) (см., например,[4]).Vf(k), как определено в уравнении (8), в настоящем изобретении определяется как угловая гармоническая функция, так как она возникает при преобразовании Фурье зависящей от угла передаточной функции и является, следовательно, решением для угловой зависимости в гармонических составляющих. Vf(k) имеет ин 002451 32 декс f, так как H в уравнении (8) представляет собой частотно-зависимую переменную. В радиолокаторе со ступенчатым изменением частоты, где H(n) измеряется на М частотах для каждого n, Vf(k) становится двумерной функцией с МхN точками, т.е. Vf(n, m), n=0 N-1 иk используется в качестве индекса для угловой гармонической переменной в дискретной форме. Если радиолокационная система измеряет сигналы как функцию времени - как в импульсном радиолокаторе, наиболее широко распространенном в коммерчески доступных радиолокационных системах обнаружения подземных объектов - угловое преобразование Фурье измеренных данных приведет к получению функции,которая зависит от угловой гармоники и от времени. Вычисление проводится согласно (8), но часто измеренный результат будет не комплексным, а действительным. Соответственно (8) можно определить следующее: где ht - измеренная функция времени, наиболее часто - импульсный отклик, вновь при различных углах n. Если измерения используют N углов и It временных точек, то Vt становится функцией с ItxN точками, т.е. Vt(n, i), n=0 N-1,i=1 It. Оказалось, что такое решение на основе уравнения (8), (8 а) или тому подобных уравнений с последующим выбором соответствующих гармонических компонентов является чрезвычайно эффективным в выделении электромагнитных отражений от объектов, таких как кабели, трубы, в земле, если они возникают на фоне других отражений, например от почвы и т.п. В предпочтительном варианте осуществления настоящего изобретения используются дискретные углы поворота поляризации (антенной системы), и предпочтительные значения N в (8) и (8 а) пригодны для вычисления БПФ. Здесь можно сослаться на таблицу 10, где показано число поворотов антенной системы при различных N и при двух имеющихся значениях .- угловое расстояние между каждым измерением. Дискретный индекс k для угловой гармоники, когда Vf(k) или Vt(k) вычисляется в соответствии с (8) или (8 а) или эквивалентного соотношения, зависит от использованных угловых приращенийи числа измерений N, которое включено в вычисления (и, следовательно, от числа физических поворотов поляризации/ антенн). В соотношениях (8) и (8 а) нет отрицательных значений индексов. Два значения k,соответствующие вышеупомянутым С 2 или С-2,соответственно двум периодам на поворот, для некоторых типовых значений приведены в таблице 11. Примеры индексов k соответствуют двум периодам на поворот. Значения k опреде 33 лены попарно (k2, k-2),- угловой поворот между каждым измерением, N - число измерений. 6. Вращение и поступательное перемещение Настоящее изобретение использует тот факт, что электромагнитные волны излучаются в землю в процессе вращения поляризации волн. Если некоторая область почвы должна исследоваться на предмет наличия в ней кабелей или труб, необходимо, чтобы антенная система перемещалась вперед над поверхностью земли. Известные георадиолокационные способы требуют, чтобы антенная система перемещалась в соответствии с предварительно определенным законом, например, по прямым, параллельным траекториям. Настоящее изобретение обеспечивает возможность перемещения вперед над землей антенной системы по произвольной траектории,причем представленный здесь способ обеспечивает повышение надежности правильного обнаружения и определения местоположения удлиненных объектов в земле вследствие такого перемещения. Если, как показано на фиг. 8 и 9, выполнить два поворота плоскости поляризации антенной системы, оба относительно центров, лежащих вблизи, но не обязательно точно над заглубленным удлиненным объектом, и если измерить передаточную функцию по напряжению в процессе вращения в обоих случаях, то можно видеть, что между двумя наборами измерений имеют место определенные сходства и что существуют некоторые различия. Различия могут быть вызваны отклонениями в составе верхних слоев почвы и неровностями поверхности. Отклонения в составе почвы, влажности, появление камней, отбросы, металлические объекты и т.п. часто приводят к значительным изменениям в измеренных результатах при перемещении даже на несколько сантиметров. Сходство зачастую вызвано отражениями от заглубленного объекта. Мельчайшие перемещения антенны,возможно, изменяют незначительно эти отражения и их угловую зависимость, но периодичность, в основном, сохраняется. В целом, периодическая последовательность зависящих от угла отражений от удлиненных объектов будет приближенно поддерживаться в процессе незначительного перемещения, в то время как соответствующая последовательность неоднородностей, находящихся над объектом, приводит к значительным изменениям. Если вычисление компонентов Фурье согласно уравнению (8) выполняется на основе повторяемых измерений в процессе комбинации поворота и соответствующего перемещения, тоV(k), соответствующее этим значениям k, будет выделяться по отношению к другим. То, что выделяется, будет относиться к удлиненным объектам, в то время как то, что ослабляется,будет относиться к вышеупомянутым другим 34 причинам отражений. Это условие, являющееся благоприятным для обнаружения, например,кабелей и труб, обусловлено тем обстоятельством, что вычисление компонентов Фурье в соответствии с (8), в основном, состоит из суммирования произведений, что означает определение среднего значения измеренных величин. Определение среднего, таким образом, важно в связи с объектами в земле, вызывающей зависимые от углов отражения. Например, горизонтальные слои почвы будут одинаково отражать электромагнитные волны независимо от направления горизонтальной поляризации, такие отражения при вычислении (8) будут суммироваться в компонент Фурье V(0), что концептуально аналогично составляющей постоянного тока или среднему значению решения по Фурье для переменного тока. Следует отметить, что отклонения в измеренном сигнале, которые вызваны ближней зоной антенн и которые, как общепринято, не обусловлены действительными радиолокационными отражениями, могут быть, таким образом,ослаблены или выделены соответственно посредством вращения, перемещения и решения по Фурье. Это является результатом того факта,что ближние зоны антенн часто имеют хорошо определенную ориентацию. Способ, описанный здесь и обеспечивающий выделение удлиненных объектов путем вращения поляризации электромагнитных волн,может быть использован для улучшения определения местоположения, например, кабелей посредством индуктивных методов, использующих низкочастотные магнитные поля, генерируемые катушками и токами. В противоположность другим способам обнаружения кабелей и труб в земле, как, например, описано в [11], заявленный способ не требует, чтобы антенная система перемещалась с пересечением относительно указанных объектов, перемещение вдоль такого объекта приведет к усилению отражений от этого объекта. Однако данный способ может реально использоваться при перемещении вдоль произвольной траектории, если только расстояние и угол относительно обнаруживаемого объекта не изменяются слишком значительно между каждым измерением. В возможном варианте осуществления настоящего изобретения поляризация поворачивается так, что полный поворот антенн достигается за каждые несколько сантиметров перемещения. Для каждого четвертого поворота, соответственно N=32 и =45, вычисляются 32 компонента Фурье. При этих условиях нежелательные отражатели и неровности, находящиеся между антенной и удлиненным объектом, сглаживаются до такой степени, что обеспечивается надежное обнаружение. 35 7. Измерение в частотной области,радиолокатор непрерывного излучения со ступенчатым изменением частоты Настоящее изобретение также относится к способу, обеспечивающему возможность использования электромагнитных полей и волн,которые имеют синусоидальную форму и которые сами по себе не содержат информации о распространении волн переходных процессов или отражений переходных процессов от объектов в земле. Способ позволяет использовать так называемый радиолокационный способ со ступенчатым изменением частоты, при котором генерируются и излучаются синусоидальные(или почти синусоидальные) волны, генерируются вторичные волны, интенсивность и фаза которых измеряются приемной антенной, которая может быть той же самой, что и передающая антенна. Радиолокационный способ со скачкообразным изменением частоты включает передачу ряда волн с различными частотами, часто с одинаковыми приращениями по частоте между отдельными передачами сигнала. Частоты, используемые в радиолокаторе со ступенчатым изменением частоты, математически могут быть выражены следующим образом:F(m)=fm+f0, m=0, 1, М-1,(9) где f(m) - используемая частота,f - частотный интервал, т.е. наименьшая разность по частоте,m - целое число,f0 - наименьшая частота. Для простоты в практических целях могут быть использованы, например, следующие числовые параметры: f=5 МГц, М=201, f0=100 МГц,что означает, что использована 201 различная частота в интервале от 100 до 1100 МГц включительно. Порядок, в котором излучаются частоты, определяемые соотношением (9), не важен. Все частоты прикладываются к передающей антенне в течение такого времени, чтобы переходной процесс, существующий одновременно с синусоидальным изменением во времени электромагнитных полей, можно было рассматривать как несущественный по сравнению с полной зависящей от времени напряженностью поля. На практике это может потребовать интервала времени 10 мкс или менее. Измерения, осуществляемые с использованием таких синусоидальных полей или волн с подавленными переходными процессами, приводят в результате к зависимому от частоты измерению передаточной функции по напряжению между передающей и приемной антеннами. Такое измерение в частотной области само по себе не обеспечивает информации о временных задержках отражений. Возможность различения между разными временными задержками принципиально важна в радиолокационной системе для определения местоположения объектов в земле, поскольку 36 такое различение является существенным средством для сортировки, например, отражающих объектов на земле или непосредственно под землей. В случае радиолокационной системы,спроектированной как система непрерывного излучения со ступенчатым изменением частоты,становится необходимым обеспечить измерение временных задержек радиолокационного сигнала на основе измеренных синусоидальных колебаний. Традиционным способом получения таких измерений является преобразование Фурье,часто реализуемое как так называемые алгоритмы быстрого преобразования Фурье (БПФ). С помощью этого способа измеренные элементы в частотной области могут быть преобразованы в вычисленную оценку соответствующих временных задержек. Однако применение методов Фурье для преобразования из частотной области во временную область сопровождается значительной неопределенностью, поскольку все радиолокационные измерения могут быть выполнены в ограниченном частотном интервале, в соответствии с фундаментальными физическими условиями. Радиоволны в земле на практике ограничены со стороны низких частот, поскольку антенны для низких частот становятся слишком громоздкими, чтобы обеспечивать достаточную чувствительность при приеме и достаточное излучение при передаче. Применяемые частоты электромагнитных волн в почве, кроме того, на практике ограничены и со стороны высоких частот, поскольку высокие частоты затухают при распространении. Затухание на практике часто увеличивается с изменением частоты с высокой скоростью на частотах, которые, согласно общей оценке, являются практически применимыми. Поэтому радиоволны в почве могут зачастую выбираться с учетом конфликтующих требований, принимая во внимание ограниченные габариты антенн, с одной стороны, и снижение проникающей мощности излучения, с другой стороны. Частоты в интервале от 100 МГц до 1 ГГц часто применяются по практическим соображениям. Затухание при распространении радиоволн в диапазоне от 100 МГц до 1 ГГц существенно возрастает с частотой. Это очевидно частично из измерений, на которых основывается настоящее изобретение, и, кроме того, известно из ряда источников, в числе которых [6,Hoekstra m.fl.], [5, Turner], [9, Freundhofer et al.]. Кроме того, частотная зависимость затухания на практически применимых частотах имеет примерно экспоненциальную форму. Если затухание проиллюстрировать на логарифмической шкале в дБ на метр в функции от частоты в МГц, то результат даст хорошее приближение к линейной функции. Затухание, измеренное в дБ/м, таким образом, относительно линейно зависит от частоты, что означает, что затухание 37 можно выразить математически как экспоненциальную функцию частоты. Это эмпирическое соотношение используется в настоящем изобретении для получения разрешения, которое ранее не было реализовано при вычислениях временных задержек с использованием сигналов непрерывного излучения со ступенчатым изменением частоты. 8. Математическая модель электромагнитных отражений от объектов в земле и других скрытых материалов Если отклик приемной антенны при передаче электромагнитных волн, характеризуемых синусоидальным изменением во времени, измерить с помощью электронной аппаратуры, то для каждой излученной частоты может быть получено значение относительной амплитуды и фазы принятого сигнала. Для каждой использованной частоты измеряется синусоидальный отклик, который математически может быть описан следующим образом:t - время. Каждое измерение дает комплексное числоHm(t)=Amехр(jm). На приемную антенну воздействует приходящая ЭМ волна, которая с хорошим приближением может толковаться как сумма волн,образующих задержанный во времени и ослабленный отклик на излученную волну. В радиолокаторе непрерывного излучения со ступенчатым изменением частоты такие отклики имеют,разумеется, ту же частоту, что и напряжение переменного тока передающей антенны. Задержанные во времени составляющие являются следствием наличия ряда отражателей. Следовательно, математически принятый отклик может быть выражен следующим образом: Здесь измеряемая величина Нm с амплитудой Аm и начальной фазой m выражена как сумма L составляющих, каждая с амплитудой ak,m и с фазой, пропорциональной частоте fm измерений,и с задержкой k,m для одиночного отражателя. Ввиду экспоненциальной частотной зависимости затухания, ak,m можно выразить следующим образом:ck - комплексные коэффициенты, представляющие частотно-зависимое затухание, а также поворот фазы, что, в числе прочего, обусловлено эффективным коэффициентом отра 002451 38 жения для одиночного отражателя. В (13) исключен индекс m, так как в математической модели предполагаем, что сk, rk, k не зависят от частоты. Если осуществлять измерения на ряде частот, которые эквидистантно расположены на интервале, как выражено в (9), то можно представить (13) в виде или В (15) ck=ckexp[(rk+jk)2fo], и, следовательно, получено уравнение, в котором комплексные экспоненты пропорциональны m, поскольку f - постоянная величина. Основываясь на этих уравнениях и в результате частотной характеристики передаточной функции, измеренной на ряде дискретных частот, как указано в (9), можно записать систему уравнений в следующем виде: НM=C1(P1)M+C2(P2)M+ CL(PL)M Данная система состоит из М уравнений с 2L неизвестными комплексными компонентами. Обычно будет обеспечено больше измерений,чем неизвестных компонентов, т.е. М 2L, и,следовательно, система уравнений становится переопределенной. Такая переопределенная система уравнений должна решаться способом наименьших квадратов. Система уравнений в(16), основанная на М измерениях, может быть решена посредством модифицированного так называемого Рrоnу-способа, использующего декомпозицию вырожденных значений (ДВЗ) и/или так называемую QR-декомпозицию. Хорошие результаты могут быть получены с помощью модифицированного Рrоnу-способа и с помощью Shanks-способа в комбинации с QR-декомпозицией или ДВЗ. Решение вышеупомянутой системы уравнений по методу наименьших квадратов приводит к получению 2L комплексных чисел: Сk=Вkехр(jk) или Pk=exp(-rk+jk)f, k=1L, (17) где Bk - амплитуда,k - фаза,rk - коэффициент затухания,k - временная задержка для k-го отражателя,f - использованный интервал дискретизации в частотной области. 39 Описанная здесь математическая модель для измеряемых радиолокационных сигналов сформулирована в соответствии с тем фактом,что затухание в почве проявляет примерно экспоненциальную частотную зависимость. Модель обеспечивает улучшенное решение для временных задержек по сравнению с методом Фурье при идентичных физических условиях. Для решения уравнений, использованных в модели, можно использовать практически различные аппроксимации и формы записи, часть которых очевидна из возможных вариантов осуществления, включая прилагаемые компьютерные программы. Практическая проблема в решении соотношения (17) состоит в том, что временная задержка для отражателей будет различной при различных углах поворота, что является результатом электрического шума и вышеупомянутых неизбежных отклонений характеристик почвы и т.п. Это является практической проблемой, так как необходимо обеспечить возможность выполнения углового Фурьеанализа при некоторых выбранных временных задержках. По этой причине и по другим причинам, связанным с вычислениями, может быть использована альтернативная математическая формулировка уравнений (15)-(17).P1 PL в (16) и (17) могут трактоваться как корни так называемого характеристического полинома, полученного из однородного дифференциального уравнения, описывающего контекст между значениями Нm для различных m. Уравнение (16) может быть сформулировано в виде так называемой модели рационального полинома видаB(z) числителя и с полиномом A(z) знаменателя по переменной z, представляющей оператор частотного сдвига для частотного сдвига f. Следует отметить, что Нm, как выражено в (16),может рассматриваться как разложение в ряд по корням полинома A(z). Представление (18) также определяется как функция с авторегрессивным смещением среднего. Коэффициенты в (18) и, следовательно, коэффициенты в (16) могут быть аппроксимированы посредством ряда методов. Возможный вариант осуществления изобретения использует так называемый Рronyспособ в модифицированной форме для определения полинома знаменателя и Shanks-способ для определения полинома числителя; они реализованы в прилагаемых компьютерных программах и использованы для вычисления кривых, иллюстрируемых в данном описании. Прилагаемые компьютерные программы иллюстрируют возможные варианты осуществления алгоритмов, которые позволили получить хоро 002451 40 шие результаты, причем их можно выполнить посредством общедоступных электронных интегральных схем, выполняющих вычислительные функции, причем за такое малое время, что изобретение можно использовать в устройстве,которое перемещается с довольно высокой скоростью над землей. Если функция H(z), описываемая уравнением (18), определена в форме коэффициентовb0, bq или a1, aL и тем самым получено функциональное уравнение для аппроксимирования частотно-зависимой передаточной функции (по одной на угол, при котором осуществлено измерение), то легко можно вычислить соответствующую временную функцию. Соответствующая временная функция вычисляется путем подстановки z в (18):(20) Следует подчеркнуть, что временная функция здесь не идентична временной функции в уравнениях (2) и (3). Зависящая от времени функция h(t), полученная из (19) и (20), не существует физически, но она оказалась полезной для обнаружения объектов в земле. Функция h(t) является непрерывной (и периодической по t), но часто на практике вычисляется для ряда дискретных значений во времени t. Математическая модель, представленная здесь, основывается на ряде предположений относительно измеренных частотных данных и,следовательно, о физических свойствах окружающей среды, включая следующее:- затухание электромагнитных волн при распространении имеет экспоненциальную частотную зависимость,- временные задержки являются частотно независимыми,- число отражателей конечно, здесь определено как L. В реальности ни одно из этих предварительных условий не выполняется полностью. Существенным преимуществом изобретения является то, что оно обеспечивает возможность надежного определения уровня и временной задержки сигналов существенных отражателей даже в ситуации, когда упомянутые идеальные условия не имеют места. Необходимо обратить внимание на следующее. Обсуждаемая промоделированная здесь функция Нm сформирована как измеряемый частотный отклик в радиолокаторе непрерывного излучения со ступенчато изменяемой частотой и поэтому является комплексной с фазовым и амплитудным значением. Как вариант,Нm может быть представлено действительной и мнимой частями. Помимо шума и иных неточностей, действительная часть и мнимая часть Нm содержат ту же самую информацию, и обе они 41 достаточны для выполнения вычислений коэффициентов в (18) или в (16). Вышеуказанная модель может быть успешно использована для моделирования действительных последовательностей. Соотношение(16 а) может быть модифицировано в следующем виде:Pk в вышеприведенных соотношениях будет существовать в виде комплексносопряженных пар, чтобы суммирование было действительным. Соответствующие соотношения применимы для корней в A(z) в соотношении (18). Могут иметься некоторые практические преимущества в использовании только действительной или только мнимой частей Нm. В возможном варианте заявленного способа используется Rе(Нm). Преимущество, главным образом, связано с измерением Нm, которое в представленном варианте осуществления реализуется в системе, использующей как аналоговые, так и цифровые методы. 9. Двумерный анализ сигналов Если используются сигналы радиолокатора непрерывного излучения, вышеописанное решение Фурье для зависимости измерений от угла поворота антенной системы, расположенной в произвольном угловом положении, зачастую оказывается недостаточным само по себе для достижения надежного обнаружения удлиненных объектов в земле. Во многих случаях в полном измеренном сигнале будут преобладать отражения, нежелательные в данном случае и обусловленные поверхностью земли и неоднородностями в земле. Поэтому непрерывный сигнал лишь в незначительной степени будет проявлять угловую периодичность, хотя он содержит отражение, например, от кабеля в земле. Для получения наилучшего возможного выделения угловой зависимости в поляризационночувствительных отражениях от объектов полезно разложить измеренный сигнал как по углу,так и по времени. Если сигнал измеряется во времени и результатом является импульс или иной переходный процесс, электрическое напряжение, приложенное к передающей антенне,будет достаточным для осуществления углового решения Фурье для различных сигнальных составляющих, соответствующих различным временным задержкам. Но если сигнал представляет собой непрерывный сигнал со ступенчатым изменением частоты, комбинация вышеописанных двух методов анализа может быть использована для углового разрешения Фурье и временного анализа на основе экспоненциального моделирования соответственно. Функция H, введенная в разделе 3 и 4,представляющая комплексную передаточную функцию, измеряется при различных углах по 002451 42 воротаи при различных частотах и может поэтому рассматриваться как функция двух переменных: углаи частоты f, т.е. Н(, f). Поскольку на практике измерение выполняется при дискретных углах и дискретных частотах, и если используетсяf(m)=fm+f0, m=0, 1, М-1 и =n, n=0, 1, N-1,то получим Н в виде: Н(n, fm) или просто Н(m, n), m=0, 1, М-1, n=0, 1, N-1 (21) Последняя компактная форма в (21) применяется, когда n и fm определены. Н(m, n) это набор чисел с NхМ комплексными элементами. Если производится преобразование частотной переменной во время, например, методом, описанным в разделе (8), для каждого из значений , т.е. для каждого n, то получим N зависящих от времени функций, в целом рассматриваемые как дискретная двумерная функция: пои по времени t: h(, t). На практикеVf(k), как определено в уравнении (8) в разделе 4, является преобразованием Фурье функции Н и, следовательно, не дает непосредственной информации о временной зависимости измеренного значения. Поэтому на практике является важным преимуществом осуществление преобразования Фурье функции h(, t) или h(n, i) в угол или в n соответственно, чтобы выделить содержимое угловых гармоник сигнала при различных временных задержках (и, следовательно, глубинах). Таким образом, на практике используется очень полезная альтернатива соотношения (8):(23), представляющая фиксированное значение временной задержки. Если (23) вычислить для ряда значений i,то получим двумерный набор чисел по k, i, где k соответствует угловой гармонике, а i соответствует квантованной временной задержке, т.е. Число точек во времени равно: It=i2-i1+1. Временной индекс проходит от i1 до i2, включая оба эти значения, соответственно подходящему временному интервалу.Vt(k, i) формируется путем преобразования Фурье It функций для каждых N точек. При этом используются не действительно двумерные функции, а комбинации одномерных функций(последовательности чисел). В (24) преобразование частоты во время выполняется при наличии результата h(n, i) перед вычислением угловой гармоники. Альтернативный способ вычис 43 ления функции, соответствующей Vt(k, i), состоит в следующем: сначала преобразование Фурье измеренных частотных данных в угол,затем преобразование частоты во время. Не приводится функциональное уравнение для процедуры по Рrоnу-способу, описанному в разделе 8, для преобразования частоты во время. Для удобства здесь используется следующая запись: x(t)=Prony[X(f)] для процедуры в целом. После этого альтернативное вычисление Vt(k, i) может быть записано в краткой форме: Следует обратить внимание, что m, n, k и i индексы для частоты, угла, угловой гармоники и времени соответственно. 10. Кривые обнаружения Угловая зависимость на различных глубинах в земле соответственно различным расстояниям, основанным на времени, будет наиболее часто переменной. Аналогичным образом, это применимо к окружающей обстановке для антенн в целом. Например, удлиненный отражающий объект будет демонстрировать периодичность, примерно сходную с описываемыми уравнениями (4) и (5), в то время как материал над объектом будет проявлять иную, возможно меньшую, степень периодичности. Путем вычисления уровня угловой гармоники, соответствующей двум периодам на вращение, при различных временных задержках может быть получена полезная кривая обнаружения, которая содержит в себе функцию времени (глубины) и которая концентрирует результат измерений и связанные с этим вычисления в наглядном виде. Кривые обнаружения вычисляются путем суммирования соответствующей гармоники вVt(k, i). Примерами формул для кривых обнаружения являются:D1(i)=Vt(k1, i)+Vt(K2, i), i=i1i2,(26) где k1 и k2 соответствуют выбранной угловой гармонике, например, соответствующей двум периодам на поворот. Значение D1 вычисляется для каждого выбранного значения i, соответствующего каждой анализируемой временной задержке. Необходимо отметить следующее. Два значения k используются потому, что вычисление коэффициента Фурье в соответствии с формулами (24) или (25) приводит к вычислению того, что соответствует положительной и отрицательной частоте для каждой гармоники, что является математически обусловленным следствием используемого определения коэффициентов Фурье и фундаментальных значений. Если разлагаемая в ряд по Фурье последовательность выполняет определенные требования, то одним из двух связанных коэффициентов можно пренебречь, но в общем случае это не имеет места. 44 При идеальных условиях угловая зависимость,соответствующая уравнениям (5 а 1) или (5 а 2),будет получена с использованием круговой поляризации. Такие угловые зависимости приводят к определенному значению k при разложении угловых гармоник по Фурье. Другими примерами кривых обнаружения,которые можно использовать, являются следующие:(28) В общем виде эти кривые обнаружения можно выразить следующим образом: где Wkпредставляет подходящую функцию выбранной гармоники. Примеры приведены в(26), (27) и (28). А - множество, содержащее выбранные гармоники. Пример. Необходимо обнаружить удлиненные объекты, причем предполагается, что они имеют относительно высокое содержание гармоник, соответствующих двум периодам на поворот антенн, которые имеют линейную поляризацию, т.е. типа 1. Измерение выполняется с угловыми интервалами 22,5 и рассматривается 4 поворота сразу, т.е.: =2/16 (радиан) иN=64. В этом случае релевантные значения k таковы: k1=8, k2=56. Можно сослаться на уравнения (24) и (25). Поскольку локальные неоднородности вблизи антенны приводят к распространению энергии сигнала по большому количеству угловых гармоник, то имеет смысл взвесить уравнение (29) по отношению к энергии гармоник за пределами множества А, что может быть реализовано следующим образом: Данное выражение представляет измерение относительного уровня выбранных гармоник при данной временной задержке. Кривые обнаружения могут быть в общем случае получены из функций Vt(k, i), характеризующих поляризационно-чувствительный временной отклик угловых гармоник от объектов. Изменения времени, связанные с местоположением объекта, в измеренных сигналах могут быть рассчитаны в среднем, путем определения кривых обнаружения в малом временном интервале, соответствующем большинству значений индекса i. Кривые обнаружения концентрируют большое число измерений в единую кривую,выражающую определение того, существуют ли объекты с заданной периодичностью или постоянной угловой зависимостью в сканируемых по углам областях на земле. 45 При вычислениях D(i) в типовом случае включаются только две гармоники (2 значенияk), что означает, что необходимое частотновременное преобразование с использованием способа на основе радиолокатора непрерывного излучения со ступенчато изменяемой частотой требуется выполнить только для двух релевантных гармоник, если перед частотно-временным преобразованием выполнено разложение по Фурье по угловым гармоникам. Это соответствует использованию уравнения (25), символизирующего процедуру, включающую в себя Рrоnуспособ. Рrоnу-способы связаны с довольно значительными вычислениями. Поэтому способ,начинающийся с разложения по Фурье по угловым гармоникам с последующим вычислением по Рrоnу-способу в типовом случае двух выбранных гармоник и с заключительным вычислением кривой обнаружения, представляет собой наиболее эффективный способ. Данный способ проиллюстрирован с помощью примеров, приведенных в приложении. В вышеизложенном описании речь шла о выделении удлиненных объектов, но подобные кривые обнаружения могут быть сформированы для любых типов объектов. Объекты с вращательной симметрией (относительно линии визирования антенн), например, проявляют относительно высокое содержание среднего значения в ряду Фурье соответственно k=0 в формулах (24) и (25). Поэтому А в (29) может быть в этом случае ограничено до А=0. 11. Определение направления удлиненных объектов Настоящее изобретение обеспечивает возможность определения направления удлиненных объектов без использования эталонных измерений. Определение направления выполняется в связи с известным физическим выполнением и ориентацией антенн. Определение направления производится с учетом множества кривых, представляющих уровень отражений, измеренных при соответствующих значениях углов поворота и временных задержек, например полученных с помощью методов, описанных выше. Используется математическая запись,приведенная выше, но когерентные значения угла и временные задержки, полученные иными способами, могут дополнительно использоваться в качестве исходной точки для определения направления. Набор чисел, символизируемыйht(n, i) в уравнениях (24) или (25), при наличии отражений от удлиненного объекта, будет показывать угловую периодичность в определенном временном интервале, где объект позиционируется по времени, как показано на графиках, иллюстрирующих периодическую последовательность на перпендикулярном поперечном сечении при 31 или при 22,5 нc. Индекс i в ht(n, i) соответствует времени, и если ht(n, i) представить графически как функцию от n для фиксированного значения i=ia вблизи связанного с 46 временем положения объекта, то получится частично периодическая последовательность по n,соответствующая углу. Графически представленная функция может быть математически записана как(30a) Как вариант, может быть использовано соотношение Аv(n, ia)=ht(n, ia)ht(n, ia), i=ia, n=0N-1(30b) Здесь символсоответствует комплексному сопряжению. В процессе радиолокационных измерений антенная система вращается относительно удлиненного объекта в земле, начиная с неизвестного начального угла 0. Если используемая антенная система работает на совпадающей поляризации (тип 1), локальные максимумы Av(n,ia) будут иметься для k значений, соответствующих 0+0, 180 и т.д. Если антенная система работает на кросс-поляризации (тип 2), то будут существовать максимумы при значениях углов 0+45, 90, 135 и т.д. Следует отметить, что частота для выражения Av(n, ia) по углу удвоена для антенны типа 2 по сравнению с антенной системой типа 1. Кривые для примеров, иллюстрируемых чертежами, соответствуют антенной системе типа 1. На практике угловая зависимость ht(n, ia) в(30) не будет соответствовать простым уравнениям, как в случае уравнений (4) и (5). Вследствие этого угловое определение не выполняется непосредственно на основе Av(n, ia), а основывается на преобразовании Фурье от Av(n, ia). Сканированные максимумы для Av(n, ia) совпадают с максимумами в синусоидальной кривой, соответствующей определенной гармонике в ряду Фурье Av(n, ia). Необходимая гармоника является одной соответствующей двум максимумам на поворот для антенн типа 1 и четырем максимумам на поворот для антенн типа 2. Термин поворот определен как вращение поляризации,возможно путем вращения антенной системы. Фаза релевантной комплексной составляющей Фурье прямо пропорциональна смещению максимумов, упомянутых в связи с положением, с которым совпадает первый максимум при n=0. Поэтому значение индекса n, соответствующее смещению, вычисляется из фазы для релевантной гармоники. Индекс гармоники k2, соответствующий двум максимумам на поворот, может быть получен следующим образом:k2=N/180 ( в градусах),(31 а) а соответствующий четырем максимумам на поворот - как показано ниже:k4=N/180 ( в градусах) Соответствующая составляющая Фурье вычисляется следующим образом:Xv(ia) - комплексное число, из фазы которого просто определяется начальный угол антенны по отношению к удлиненному объекту,как описано ниже. Антенны типа 1. Смещение соответствующих максимумов в связи с n=0 может быть определено относительно n=0 следующим образом:n0 - не является целым. Следовательно,угол объекта по отношению к линейной поляризации антенной системы типа 1 равен: 0=a/2 (все углы в градусах),a=Xv(ia)=фаза k2-гармоник(32b) Следует отметить, что начальный угол антенной системы относительно объекта, вычисляемый как 0 в (32), может быть получен с лучшим разрешением, чем соответствующий приращению угламежду каждым измерением. Можно использовать возможность интерполяции разрешения Фурье и 0 можно в принципе определить с произвольным разрешением, ограничиваемым только шумом, неточностями измерения и т.д. Смысл 0 поясняется в связи с фиг. 7b. Антенны типа 2. В данном случае число периодов в Xv(ia) на поворот удваивается по сравнению с антеннами типа 1. Смысл 0 поясняется в связи с фиг. 7 а,причем необходимо дополнительно ввести угол 45. Другие типы антенн. Антенны с круговой или эллиптической поляризацией могут рассматриваться как линейные комбинации антенн типа 2, а ориентация объектов может быть определена по процедуре, соответствующей представленной в настоящем описании. 12. Возможные способы обнаружения и определения местоположения кабелей и труб в почве Вышеописанные математические методы для анализа сигнала в области углов/угловых гармоник и в частотно-временной области соответственно могут быть эффективным образом объединены различными путями. В последующем будут описаны некоторые такие комбинации, применимые на практике. Ниже приводятся двумерные множества чисел, поскольку наиболее часто наиболее соответственными являются дискретные значения. Эти множества чисел могут трактоваться как элементы матрицы, и каждая из строк или столбцов в матрице может рассматриваться как последовательность. Кривые или функции и функциональные значения описываются синонимичным образом. Предпочтительные варианты осуществления, представленные изобретателями, основываются на дискретных частотах,углах, угловых гармониках и моментах времени, 002451 48 что обосновывается применением цифровой памяти и вычислительных блоков, несмотря на то, что представленные примеры реализуются посредством непрерывных функций. Принципы изобретения применимы также к непрерывным частотам, углам поворота, угловым гармоникам и моментам времени. Процедура 1. Содержит следующие этапы: 1.1. Антенное устройство типа 1 или 2, как описано в разделах 3 и 4, размещается на земле или над землей при произвольной ориентации,однако примерно горизонтально. 1.2. Передаточная функция по напряжению между передающей и приемной антеннами измеряется на ряде частот. Измеренные значения запоминаются. 1.3. Антенное устройство поворачивается механически или электрически на целое кратное угла . Типовое значениеможет составлять 22,5 или 45 соответственно 1/16 или 1/8 полного поворота. Одновременно с вращением антенна может быть незначительно перемещена поступательно, в данном случае это соответствует перемещению антенны в процессе ее вращения над землей. 1.4. Передаточная функция по напряжению между передающей и приемной антеннами измеряется на нескольких частотах. Измеренные значения запоминаются. 1.5. Этапы 1.3. и 1.4. повторяются до тех пор, пока не будет выполнено соответствующее число измерений при соответствующем числе угловых позиций. Это число зависит от конкретных обстоятельств. Практически приемлемыми количествами используемых углов могут быть, например, следующие: 8, 16, 32, 64. В результате получают множество комплексных чисел, представляющих измеренные значения передаточной функции по напряжению для когерентных значений углового поворота и частоты. Результат может быть сведен к множеству действительных чисел соответственно действительным частям упомянутого множества комплексных чисел. 1.6. Полученные результаты измерений преобразуются из частотной области во временную область. Для этой цели применяется математический метод, основанный на экспоненциальной модели, описанной в разделе (8). Результатом является ряд оценок отражателей, их уровень, фаза и временная задержка. Преобразование из частотной области во временную область выполняется для каждой частотно-зависимой передаточной функции, связанной с конкретным углом. Таким образом, для каждого угла получается зависящая от времени непрерывная во времени периодическая функция, которая вычисляется для соответствующего числа дискретных моментов времени. Одно и то же время вычислений используется для всех углов. Результатом является множество двумерных ком 49 плексных чисел. Последовательность соответствует математической функции, имеющей дискретизированный угол поворота и время, дискретизированное как независимая переменная. Комментарий: Примеры результата процедуры 1.1.-1.6. показаны на фиг. 13 а, 14 а, 15 а,16 а, 17 а. Во многих случаях при рассмотрении этих результатов могут быть выявлены периодические по углу объекты. На фиг. 13 а иллюстрируется объект, представляющий собой трубу из пластика длиной 160 мм, содержащую воду,при временной задержке около 31 нс на глубине 135 см в почве. Процедура 2. 2.1. Этапы 1.1.-1.6. соответствуют процедуре 1. 2.2. Для каждой точки во времени, использованной на этапе 2.1., выполняется преобразование Фурье с использованием угла в качестве независимой переменной. В результате получают двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, имеющей угловую гармонику и время в качестве независимых переменных. 2.3. На основе двумерного множества чисел, полученного на этапе 2.2., вычисляется одномерная функция по времени D(t) с использованием одного из методов, описанных выше в разделе 10. 2.4. Дискретная функция времени, полученная на этапе 2.3., сканируется для нахождения локальных и глобальных максимумов, которые демонстрируют значительный уровень по сравнению со средним значением. Такие максимумы представляют периодические по углу отражения от удлиненных объектов. Значение времени, соответствующее индивидуальному максимуму, устанавливает временную задержку отражений от соответствующего объекта. Сканирование для нахождения максимумов может выполняться механически или путем анализа кривой. Комментарий: Временные задержки, вычисленные таким образом, могут, с учетом известности скорости распространения радиоволн в земле, использоваться для вычисления глубины объектов. Примеры результата этапов 2.1.2.2. процедуры показаны на фиг. 13 с, 14 с, 15 с,16 с, 17 с. Примеры результата этапов 2.1.- 2.3. процедуры показаны на фиг. 13 е, 14 е, 15 е, 16 е, 17 е,где представлены кривые обнаружения, вычисленные в соответствии с (28). Процедура 3. Состоит из процедуры 1 или соответствующих разделов и дополнительных вычислений. 3.1. Этапы 1.1.-1.6. соответствуют аналогичным этапам процедуры 1. 3.2. На основе результата этапа 3.1. вычисляется угловое положение относительно антенной системы с использованием метода, описан 002451 50 ного в разделе 11. Угловое положение вычисляется для одного или нескольких периодических по углу объектов, которые могут быть обнаружены или для которых может быть определено время путем визуального анализа результата этапов 1.1.-1.6. по процедуре 1 или с использованием кривых обнаружения, как описано в процедуре 2. Комментарий: Если антенная система предусматривает механическое вращение, то при этом определяется угловое положение в начальном положении относительно удлиненных объектов. Если вращение выполняется путем сдвигов между антенными элементами при фиксированном угловом положении, то угловое положение определяется направлением поляризации,соответствующим первому измерительному положению, и тем самым угловое положение всего антенного устройства по отношению к удлиненному объекту (объектам). Примеры результата этой процедуры представлены на фиг. 15k,l. Процедура 4. Содержит следующие этапы: 4.1. Антенна типа 1 или 2, как описано в разделе 3 и 4, размещена на земле или над землей при произвольной, однако приближенно горизонтальной ориентации. 4.2. Передаточная функция по напряжению между передающей и приемной антеннами измеряется на ряде частот. Измеренные значения запоминаются. 4.3. Антенное устройство поворачивается механически или электрически на целое кратное угла . Типовое значениеможет составлять 22,5 или 45 соответственно 1/16 или 1/8 полного поворота. Одновременно с вращением антенна может быть незначительно перемещена поступательно, в данном случае это соответствует перемещению антенны в процессе ее вращения над землей. 4.4. Передаточная функция по напряжению между передающей и приемной антеннами измеряется на нескольких частотах. Измеренные значения запоминаются. 4.5. Этапы 4.3. и 4.4. повторяются до тех пор, пока не будет выполнено соответствующее число измерений при соответствующем числе угловых позиций. Это число зависит от конкретных обстоятельств. Практически приемлемыми количествами используемых углов могут быть, например, следующие: 8, 16, 32, 64. В результате получают множество комплексных чисел, представляющих измеренные значения передаточной функции по напряжению для когерентных значений углового поворота и частоты. Результат может быть сведен к множеству действительных чисел соответственно действительным частям упомянутого множества комплексных чисел. 4.6. Полученные на этапе 4.5. результаты измерений подвергаются преобразованию Фу 51 рье по углу. Для каждой из частот, использованных на этапе 4.5., выполняется преобразование Фурье с использованием угла в качестве независимой переменной. В результате получается двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, использующей угловую гармонику и частоту измерений в качестве независимой переменной. Комментарий: Полученное множество чисел является двумерным, что может трактоваться как ряд частотно-зависимых кривых, по одной на каждую угловую гармонику, или как ряд функций гармоники, по одной на каждую частоту измерений. Процедура 5. Содержит этапы процедуры 4 или соответствующие им этапы и дополнительные вычисления. 5.1. Этапы 4.1.-4.6. по процедуре 4 или соответствующие им этапы. 5.2. Множество чисел, полученных на этапе 5.1., преобразуются из частотной области во временную область. Ряды, соответствующие каждой из угловых гармоник, преобразуются посредством математического метода, основанного на экспоненциальной модели, или соответствующего метода, описанного в разделе 8. Результатом является ряд оценок отражателей, их уровень, фаза и временная задержка. Таким образом, для каждой угловой гармоники получают зависящую от времени непрерывную по времени и периодическую функцию, которая вычисляется для соответствующего числа дискретных моментов времени. Идентичные точки вычислений во времени используются для каждой из угловых гармоник. В результате получают двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, использующей дискретизированную угловую гармонику и дискретизированное время в качестве независимых переменных. Комментарий: Полученное множество чисел является двумерным, что может трактоваться как ряд зависимых от времени кривых, по одной на каждую угловую гармонику, или как ряд кривых или функций гармоники, по одной на каждый момент времени. Данное множество чисел может быть проиллюстрировано, и при анализе могут быть обнаружены объекты путем идентификации глобальных или локальных максимумов для выделенных гармоник, в типовом случае двух гармоник соответственно двум периодам на поворот поляризации. Процедура 6. Содержит этапы процедуры 4 или соответствующие им этапы и дополнительные вычисления. 6.1. Этапы 4.1.-4.6. по процедуре 4 или соответствующие им этапы. 6.2. Множество чисел, полученных на этапе 6.1., преобразуются из частотной области во 52 временную область. Ряды, соответствующие каждой из угловых гармоник, преобразуются посредством математического метода, основанного на экспоненциальной модели, или соответствующего метода, описанного в разделе 8. Результатом являются ряд оценок отражателей, их уровень, фаза и временная задержка. Таким образом, для каждой угловой гармоники получают зависящую от времени непрерывную по времени и периодическую функцию, которая вычисляется для соответствующего числа дискретных моментов времени. Идентичные точки вычислений во времени используются для каждой из угловых гармоник. В результате получают двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, использующей дискретизированную угловую гармонику и дискретизированное время в качестве независимых переменных, но только для выбранных угловых гармоник. 6.3. На основе двумерного множества чисел, полученного на этапе 6.2., вычисляется одномерная функция времени D(t) с использованием методов, описанных выше в разделе 10. 6.4. Дискретная функция времени, полученная на этапе 6.3., сканируется для нахождения локальных и глобальных максимумов, которые демонстрируют значительный уровень по сравнению со средним значением. Такие максимумы представляют периодические по углу отражения от удлиненных объектов. Значение времени, соответствующее индивидуальному максимуму, устанавливает временную задержку отражений от соответствующего объекта. Сканирование для нахождения максимумов может выполняться механически или путем анализа кривой. Комментарий: Как описано в разделе 10,две из угловых гармоник часто выбираются в соответствии с двумя периодами на поворот. Если антенна типа 1, как описано выше, часто используются две гармоники, которые соответствуют положительной или соответственно отрицательной частоте, соответствующей двум периодам на поворот. Поскольку только временные последовательности требуют вычисления, например, двух гармоник, и, следовательно,используются только два Рrоnу-преобразования,то процедура является весьма эффективной с точки зрения вычислений. Временные задержки,вычисленные таким способом, могут, с учетом известности скорости распространения радиоволн в земле, использоваться для вычисления глубины объектов. Примеры результата этапа 6.1. процедуры 6 представлены на фиг. 13D,14D, 15D, 16D, 17DВ, на которых кривые обнаружения вычислены в соответствии с уравнением (28). Процедура 7. Содержит следующие этапы: 7.1. Антенное устройство типа 1 или 2, как описано в разделе 3, размещается на земле или 53 над землей с произвольной, однако примерно горизонтальной ориентацией. 7.2. Импульс напряжения, который должен быть импульсом с коротким периодом, прикладывается к передающей антенне. Напряжение передающей антенны затем измеряется для ряда временных задержек. 7.3. Антенное устройство поворачивается механически или электрически на целое кратное угла . Типовое значениеможет составлять 22,5 или 45 соответственно 1/16 или 1/8 полного поворота. Одновременно с вращением антенна может быть незначительно перемещена поступательно, в данном случае это соответствует перемещению антенны в процессе ее вращения над землей. 7.4. Импульс напряжения, который должен быть импульсом с коротким периодом, прикладывается к передающей антенне. Напряжение приемной антенны затем измеряется для ряда временных задержек. 7.5. Этапы 7.3. и 7.4. повторяются до тех пор, пока не будет выполнено соответствующее число измерений при соответствующем числе угловых позиций. Это число зависит от конкретных обстоятельств. Практически приемлемыми количествами используемых углов могут быть, например, следующие: 8, 16, 32, 64. В результате получают множество чисел, представляющих измеренные значения передаточной функции по напряжению для когерентных значений углового поворота и частоты. Множество чисел часто будет действительным. Если антенное устройство выполнено так, чтобы измерять две временные функции,которые взаимно ортогональны, - часто определяются как I(t) и Q(t) - множество чисел будет комплексным. Если измеренное множество чисел является действительным, то можно сформировать соответствующее ему множество комплексных чисел - в форме так называемого преобразования Гильберта - либо путем электрической фильтрации, либо путем вычислений с использованием измеренных дискретных значений. Процедура 8. Содержит этапы процедуры 7 или соответствующие им этапы и дополнительные вычисления. 8.1. Этапы 7.1.-7.5. по процедуре 7 или соответствующие им этапы. 8.2. Для каждой точки во времени, использованной на этапе 8.1. (=временным задержкам),выполняется преобразование Фурье с использованием угла в качестве независимой переменной. В результате получают двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, имеющей угловую гармонику и время в качестве независимых переменных. 8.3. На основе двумерного множества чисел, полученного на этапе 8.2., вычисляется од 002451 54 номерная функция по времени D(t) с использованием одного из методов, описанных выше в разделе 10, или соответствующих методов. 8.4. Дискретная функция времени, полученная на этапе 8.3., сканируется для нахождения локальных и глобальных максимумов, которые демонстрируют значительный уровень по сравнению со средним значением. Такие максимумы представляют периодические по углу отражения от удлиненных объектов. Значение времени, соответствующее индивидуальному максимуму, устанавливает временную задержку отражений от соответствующего объекта. Сканирование для нахождения максимумов может выполняться механически или путем анализа кривой. Комментарий: Временные задержки, вычисленные таким образом, могут, с учетом известности скорости распространения радиоволн в земле, использоваться для вычисления глубины объектов. Процедура 9. Содержит этапы процедуры 7 или соответствующей ей процедуры и дополнительные вычисления. 9.1. Этапы 7.1.-7.5. по процедуре 7 или соответствующие им этапы. 9.2. На основе результата этапа 9.1. вычисляется угловое положение относительно антенной системы с использованием метода, описанного в разделе 11. Угловое положение вычисляется для одного или нескольких периодических по углу объектов, которые могут быть обнаружены или для которых может быть определено время путем визуального анализа результата этапа 9.1. или с использованием кривых обнаружения, как описано в процедуре 8. Комментарий: Если антенная система предусматривает механическое вращение, то при этом определяется угловое положение в начальном положении относительно удлиненных объектов. Если вращение выполняется путем сдвигов между антенными элементами при фиксированном угловом положении, то угловое положение определяется направлением поляризации,соответствующим первому измерительному положению, и тем самым угловое положение всего антенного устройства по отношению к удлиненному объекту (объектам). Процедура 10. Содержит этапы процедуры 7 или соответствующие им этапы и дополнительные вычисления. 10.1. Этапы 7.1.-7.5. по процедуре 7 или соответствующие им этапы. 10.2. Для каждой точки во времени, использованной на этапе 10.1. (=временным задержкам), выполняется преобразование Фурье с использованием угла в качестве независимой переменной. В результате получают двумерное множество комплексных чисел. Это множество чисел соответствует математической функции, 55 имеющей угловую гармонику и время в качестве независимых переменных. Комментарий: Результат преобразования Фурье измеренных радиолокационных сигналов в угол поляризации может посредством визуального анализа дать ценную информацию относительно объектов в земле. Процедуры для использования круговой или эллиптической поляризации Настоящее изобретение предусматривает также использование круговой или эллиптической поляризации. Процедуры, описанные выше, например процедуры 1 и 4, используют антенны с линейной поляризацией, вращаемой механически или электрически, например, с помощью электронного сдвига между элементами. Если использовать круговую или эллиптическую поляризацию, то процедуры будут осуществляться соответствующим способом. Волны с круговой или эллиптической поляризацией и соответствующие им сигналы антенн могут быть разложены на составляющие, соответствующие линейной поляризации. Если круговую или эллиптическую поляризацию получают с помощью линейно поляризованных антенных элементов и сигналы на их выводах можно измерить, то способы, описанные выше, можно применить непосредственно к индивидуальным линейным частям - с последующей их комбинацией или без нее. Это применимо независимо от того, возбуждаются ли линейно поляризованные элементы одновременно или по отдельности, и независимо оттого,существуют ли физически при измерениях круговая или эллиптическая поляризации, или соответствующие сигналы формируются в виде линейных комбинаций результатов измерений. Настоящее изобретение также предусматривает то, что круговая или эллиптическая поляризации могут быть получены посредством одной или нескольких антенн, самих по себе излучающих электромагнитные волны с круговой или эллиптической поляризацией, например спиральных антенн. Электрические сигналы на выводах таких антенн делятся на ортогональные составляющие, и эти составляющие обрабатываются посредством способов, описанных в настоящем описании. Ортогональные составляющие сигналов часто определяются как квадратурные составляющие I и Q, которые образуют в целом комплексное представление сигнала. Независимо от разделения сигнала на ортогональные составляющие, получаемого с помощью электрических схем или посредством алгебраических операций в цифровом компьютере(наиболее часто), обработка полученных сигналов производится согласно настоящему изобретению посредством описанных здесь способов. 56 Использование методов Фурье для преобразования из частотной области во временную область В процедурах 1, 5 и 6, описанных выше, используется способ, основанный на экспоненциальной модели, для определения временных задержек в радиолокационном сигнале,как описано в разделе 8. Настоящее изобретение также предусматривает процедуры, соответствующие описанным, где преобразование Фурье используется для преобразования из частотной области во временную область, включая использование решения Фурье по углу, как описано выше в комбинации с преобразованием Фурье частотно-зависимых данных измерений. В процедуре 1, этап 1.6., процедуре 5,этап 5.2., и в процедуре 6, этап 6.2., способ,основанный на экспоненциальной модели, может быть заменен на метод Фурье, включающий быстрое преобразование Фурье (БПФ). Использование других видов сигналов,отличных от сигнала радиолокатора с непрерывным излучением со ступенчатым изменением частоты Настоящее изобретение также предусматривает использование методов, описанных во взаимосвязи с другими формами сигналов, отличными от сигналов непрерывного излучения со ступенчато изменяемой частотой, которые соответствуют предпочтительному варианту осуществления. Импульсный радиолокатор используется в процедурах 7-10. Методы, описываемые здесь, в комбинации с вращающейся поляризацией с перемещением над землей или без такого перемещения,могут быть с выгодой использованы независимо от формы используемого сигнала. Могут использоваться следующие формы сигналов:- с качанием частоты, т.е. линейночастотно-модулированное непрерывное излучение, при котором частота сигнала изменяется непрерывно в течение интервала времени,- сигнал корреляционного радиолокатора,в котором используются частоты для произвольного времени, причем временные задержки отражений вычисляются посредством корреляционной обработки. Обнаружение нескольких объектов Обнаружение и определение направления описаны выше для индивидуальных объектов. Если несколько объектов присутствуют в одном и том же множестве радиолокационных измерений, методы, описанные здесь, дают возможность индивидуальной обработки для каждого из объектов, если связанные с ними временные задержки различны. Объекты на различных глубинах, в том числе пересекающиеся кабели и трубы, могут быть определены индивидуально посредством описанных способов. Модифицированные Рrоnу-способы часто позволяют селектировать близко расположенные объекты. 57 Примеры выполненных вычислений Представленные примеры основаны на измерениях, выполненных при обычных условиях в Дании. Таким образом, примеры являются полностью реалистичными и соответствуют обычно существующим условиям. Однако их не следует рассматривать как примеры вычислений, приведенные с целью характеризовать теоретические или принципиальные аспекты способов, соответствующих изобретению. Таблица 12 иллюстрирует схему объектов,представленных в примерах 1-5. Фиг. 13a-13b, 14a-14b, 15a-15b, 16a-16b и 17a-17b - это результаты вычислений типа, упомянутого в разделе 9, иллюстрирующие вычисленный уровень сигнала в функции индекса угла n и времени в нc. Результаты вышеописанной процедуры 1 на этапе 1.6. на фиг. 13 а-17 а представлены в форме трехмерных графиков и на фиг. 13b-17b - то же самое в виде контурных графиков. На фиг. 13c-13d, 14c-14d, 15c-15d, 16c-16d и 17c-17d измеренный уровень сигнала представлен как функция индекса угла n и частоты в ГГц. Результаты вышеописанной процедуры 1 на этапе 1.5. на фиг. 13 с-17 с представлены в форме трехмерных графиков и на фиг. 13d-17d - то же самое в виде контурных графиков. Фиг. 13e-13f, 14e-14f, 15e-15f, 16e-16f и 17e-17f - это результаты вычислений, подобных упомянутым в разделе 9, иллюстрирующие вычисленный уровень сигнала в функции индекса угловой гармоники и времени в нc. Результаты вышеописанной процедуры 2 на этапе 2.2. на фиг. 13 е-17 е представлены в форме трехмерных графиков и на фиг. 13f-17f - то же самое в виде контурных графиков. Фиг. 13g-13h, 14g-14h, 15g-15h, 16g-16h и 17g-17h - это результаты вычислений типа, упомянутого в разделе 10, иллюстрирующие кривые обнаружения. Представлены результаты вышеописанной процедуры 6 на этапе 6.3. Фиг. 13i-13j, 14i-14j, 15i-15j, 16i-16j и 17i17j - это результаты вычислений типа, упомянутого в разделе 10, иллюстрирующие кривые обнаружения. Представлены результаты вышеописанной процедуры 2 на этапе 2.3. Фиг. 15k-15o иллюстрируют пример определения направления с помощью вычислений,как описано в разделе 11, а также иллюстрируют результат вычислений, включенных в вышеописанную процедуру 3. Комментарий к примерам: все представленные примеры используют измерения, осуществленные с помощью антенной системы, в которой внутренняя временная задержка и временной интервал для прямой связи от передатчика к приемнику составляет примерно 11 нс. Эта включенная задержка не показана отдельно на приведенных кривых. Примечание. Показанные кривые в приведенных примерах получены путем соединения 58 дискретных значений прямыми линиями и, соответственно, не должны трактоваться как непрерывные функции или кривые. Пример 1 иллюстрирует на фиг. 13 а и 13b явную угловую гармонику примерно на 31 нc, проявляющуюся в форме явных локальных максимумов на фиг. 13 е и 13f и дающую на кривой обнаружения на фиг. 13h глобальный максимум, превышающий шумовой уровень примерно на 17 дБ. Следовательно, обеспечивается четкое обнаружение трубы, которая в случае георадиолокатора без вращения поляризации и соответствующего анализа сформировала бы сигнал на уровне многих других, без каких-либо характеристических признаков. На фиг. 13 е и 13f видно, что существенная периодичность имеет место в районе 15-20 нc, и что она селектируется процедурой в целом. Угловая гармоника, соответствующая 0 ому порядку, т.е. среднему значению, имеет различимый максимум в районе 12 нc, что также отфильтровывается. Измерения в примере 2 осуществлялись при нормальных, но предпочтительных условиях. Здесь надежное обнаружение получено посредством только одного поворота с 16 измерениями. Если число вращений уменьшается наполовину, то получаются соответствующие результаты. Видно, что исходные частотные измерения на фиг. 14 с и 14b демонстрируют значительную степень периодичности. Пример 3 иллюстрирует измерение с существенными мешающими отражениями на глубине, соответствующей временной задержке 14 нc. Однако видно, что эти мешающие отражения на окончательной кривой обнаружения,представленной на фиг. 15h, ослабляются на 20 дБ по сравнению с глобальным максимумом,расположенным примерно в районе 22 нc. Фиг. 15k-15o иллюстрируют пример определения направления. На фиг. 15m показана амплитуда функции, соответствующей функцииh(k, t) на фиг. 15 а и 15b для t=21,8 нc, что соответствует положению глобального максимума кривой обнаружения по фиг. 15g и 15h. Периодичность отчетливо видна. Угловой индекс n используется как единица по горизонтали. В середине фиг. 15n абсолютное значение вычислений БПФ показано над вышеупомянутой кривой при числе точек N=64. Нижняя часть фиг. 15 о иллюстрирует фазу двух угловых гармоник, соответствующих двум периодам на поворот (антенна типа 1) соответственно k1=8 и k2=56 при повороте антенн. Значение -0,1844 представляет собой вычисленное угловое положение, определенное в числе угловых шагов в соответствии с вышеупомянутой формулой (32b). Из верхней кривой 15m и при рассмотрении фиг. 15 а и 15b видно, что вычисленное положение объекта, почти совпадающее с положением для первого измерения (n=0) является корректным. 59 Пример 5 иллюстрирует, каким образом, несмотря на очень слабую периодичность при частотных измерениях (фиг. 17 с и 17d) и несмотря на сильные мешающие отражения на расстоянии, соответствующем примерно 12 нc(фиг. 17 а и 17b), может быть обеспечено эффективное обнаружение, иллюстрируемое на фиг. 17g и 17h, на которых глобальный максимум примерно при 12 нc может быть режектирован, если не принимается обнаружение объекта непосредственно под поверхностью, как на фиг. 17i и 17j. 14. Возможные применения Настоящее изобретение, в числе прочего,имеет следующие применения:- Бортовой обнаружитель ручного режима для определения местоположения кабелей, труб и оптических волокон. Такой обнаружитель может обеспечить оператору немедленную индикацию наличия, положения и ориентации таких объектов. Обнаружитель может перемещаться над землей по произвольному закону и отслеживать местоположение заглубленного объекта.- Соответствующий обнаружитель, устанавливаемый на транспортных средствах.- Автоматический обнаружитель, устанавливаемый на экскаваторах или иных машинах для перемещения почвы для исключения возможности повреждения объектов.- Картографирование кабелей, труб, иных подземных объектов. Обнаружитель может быть связан с аппаратурой позиционирования и аппаратурой для сбора данных. В типовом случае посредством картографирования будет выполняться сканирование по предварительно определенному закону.- Инструмент для сканирования структуры и состояния материалов, включающих бетон.- Инструмент для определения местоположения труб и кабелей в строительных конструкциях: полах, стенах, например, из бетона или иных материалов. 15. Конструкции антенн Ступенчатый поворот поляризации без использования механических средств осуществляется наилучшим образом путем сдвига между одинаково выполненными антеннами, расположенными с соответствующими углами и со смещенными центрами, как показано схематично на фиг. 3, 4 и 9. На фиг. 3 и 4 представлены линейные тонкие диполи. Для получения достаточной ширины полосы предпочтительно использовать ветви диполей конической формы, как показано на фиг. 9. Для исключения или, по меньшей мере, для ослабления отражений от концов диполей (соответствует стоячим волнам по оси диполя) часто используется некоторый тип ослабления. В предпочтительном варианте ослабление обеспечивается использованием комбинации поглощающего материала, расположенного вблизи с антеннами, над антеннами, как показано на фиг. 12, и использованием последовательных элементов резистивной связи, как показано на фиг. 60 10 и 11. Резистивные элементы могут состоять из резистора поверхностного монтажа, имеющего переменное значение, так что наивысшие значения используются на наибольшем удалении от центра диполя. Типовые значения находятся в пределах от 10 до 100 Ом. Резистивные элементы могут также состоять из последовательности, связывающей вышеупомянутый резистор и переключающий диод, например PIN-диод. Путем приложения напряжения постоянного тока между концевыми точками диполя и центром может быть обеспечена возможность перевода PIN-диодов в проводящее или непроводящее состояние. Один из эффектов такого устройства состоит в том, что диполь в режиме запирания будет иметь меньшую электромагнитную связь с рядом расположенными аналогичными диполями,чем в случае без диодов. Таким путем можно уменьшить габариты всей конструкции без нарушения связи между активными и пассивными диполями. Нежелательные связи между диполями могут обусловить искажение линейности поляризации и нежелательные помехи в сигнале. Только к тем антеннам, которые желательно перевести в активное состояние, прикладывается напряжение постоянного тока, часто сразу к двум, что приводит к переводу PIN-диодов в проводящее состояние. Коммутатор, описанный выше, содержащий последовательность PIN-диодов, обеспечивает перекрытие диполей друг с другом без существенной взаимной связи. Диполь с четырьмя металлическими поверхностями и тремя резистивными элементами связи с PIN-диодами или без них с общей длиной 24 см был продемонстрирован изобретателями и при измерениях проявил частотную характеристику относительно плоскую в области от 100 до 900 МГц. Остальная часть конструкции показана на фиг. 11 и 12. 16. Использованные источники 1. US-Patent 4,728,897, Gunton et al. 2. US-Patent 4,812,850, Gunton et al. 3. Balanis C.A., Advanced Engineering Electromagnetics, John WileySons, 1989. 4. Oppenheim A.V.Schafer R.W., Discrete-Time Signal Processing, Prentice Hall, 1989. 5. Turner G.Siggins A.F., Constant Q