Способ и система для определения положения источника звука

Способ и система для определения положения источника звука

Иллюстрации

Показать все

Реферат

ОБЛАСТЬ ТЕХНИКИ

Настоящее изобретение относится к определению положения источников звука и предлагает способ и систему для определения положения источников звука;

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

ПРЕДШЕСТВУЮЩИЙ УРОВЕНЬ ТЕХНИКИ

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

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

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

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

СУЩНОСТЬ ИЗОБРЕТЕНИЯ

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

обеспечивают несколько акустических датчиков;

с помощью этих нескольких акустических датчиков принимают акустический сигнал от источника звука и

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Согласно третьей особенности настоящего изобретения предлагается способ обследования, содержащий способ согласно первого варианта. Предлагаемый способ может быть использован для мониторинга активов, например, нефте- и газопроводов.

Факультативно способ содержит размещение распределенного оптоволоконного датчика вблизи от актива или места, подлежащего мониторингу.

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

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

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

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

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

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

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

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

Согласно пятому варианту настоящего изобретения предлагается устройство для определения положения источника звука, содержащее:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

КРАТКОЕ ОПИСАНИЕ ГРАФИЧЕСКИХ МАТЕРИАЛОВ

Настоящее изобретение раскрывается ниже на примере со ссылками на сопроводительные чертежи, на которых:

На ФИГ.1 кратко представлены различные компоненты настоящего изобретения,

На ФИГ.1b схематически показан интерферометр для измерения оптической амплитуды, фазы и частоты оптического сигнала, из которых могут быть получены акустические амплитуда, фаза и частота; интерферометр может входить в детектор согласно одному из вариантов осуществления настоящего изобретения,

На ФИГ.2 показана связь источника звука и положений обнаружения детекторами,

На ФИГ.3 показано геометрическое место точек, образованное положениями обнаружения источника звука детекторами,

На ФИГ.4 показан разрез через источник звука на поверхности земли и детектор, заглубленный под поверхность,

На ФИГ.5 показан разрез через источник звука на поверхности земли и два детектора, заглубленные под поверхность,

На ФИГ.6 показана карта поверхности с источником звука на поверхности земли и двумя детекторами, заглубленными под поверхность,

На ФИГ.7 показана карта поверхности с источником звука на поверхности земли и детектором, заглубленным под поверхность; на карте показаны также различные зоны поверхности,

На ФИГ.8 показана карта поверхности, на которой детекторы не расположены вдоль прямой линии,

На ФИГ.9 показана карта поверхности, на которой детекторы охватывают некоторую область,

На ФИГ.10 показан пример дисперсионной кривой из уровня техники и

На ФИГ.11 показан пример дисперсионного анализа и профиля сдвиговой волны для известного из уровня техники решения.

ОПИСАНИЕ ПРЕДПОЧТИТЕЛЬНЫХ ВАРИАНТОВ ОСУЩЕСТВЛЕНИЯ ИЗОБРЕТЕНИЯ

Как показано на ФИГ.1, объект настоящего изобретения состоит из одного или нескольких детекторов, причем каждый из них состоит из нескольких акустических датчиков (101), от которых могут быть получены один или несколько сигналов (102), зависящих от расстояния и времени и являющихся результатом излучения одним или несколькими оптическими источниками (100) звуковой волны, уловленной детекторами (101); средств обработки (103), извлекающих из сигнала (102) одно или несколько геометрических мест точек (104), и средств обработки (105), получающих из одного или нескольких геометрических мест точек (104) положение одного или нескольких источников (100) звука. В предпочтительных вариантах осуществления настоящего изобретения количество акустических датчиков велико, они расположены тесно и, предпочтительно, реализованы в виде распределенного акустического датчика; более предпочтительно - в виде распределенного оптоволоконного акустического датчика и еще более предпочтительно - в виде датчика, регистрирующего как акустическую амплитуду, так и фазу. Пример такого предпочтительного детектора раскрыт в патентных заявках WO 2010/136809 A2 и WO 2010/136810 A2 и будет раскрыт подробнее ниже при обсуждении ФИГ.1b.

Время первого поступления акустической волны от источника (100) звука, определенное в сигнале (102), может быть использовано средством обработки (103) для получения требуемого геометрического места точек (104), но достичь более высокого отношения сигнал/шум и, в целом, более эффективной работы при получении геометрического места точек (104) можно в условиях значительной длительности поступления энергии звуковой волны. Это требует когерентности обнаружения звуковой волны детекторами (101) и предпочтительно достигается выполнением детекторов (101) в виде оптоволоконного акустического датчика, регистрирующего как акустическую амплитуду, так и фазу.

Геометрические места точек (104) состоят из наборов временных и пространственных координат точек, в которых в сигнале (102), выданном детекторами (101), обнаруживаются характеристики, общие для сигналов, зарегистрированных несколькими акустическими датчиками детекторов (101). Средство обработки (103) может использовать любые параметры звуковой волны, обнаруженной детекторами (101), для получения требуемых геометрических мест точек (102) взаимной корреляцией, если звуковая волна передается от источников (100) звука к детекторам (101) изотропно с малой дисперсией, так что сигналы, полученные всеми акустическими датчиками (101), имеют существенно одинаковую форму. Однако в одном из предпочтительных вариантов осуществления настоящего изобретения средство обработки (103) получает фазу или фазы одного или нескольких существенно несмешанных частотных компонентов акустического спектра сигнала каждого акустического датчика детектора (101). Таким образом, может быть снижено влияние дисперсии. Специалистам в области обработки сигналов известны способы, применяемые для получения фаз, например, такие, как узкополосный фильтр и преобразования Фурье, поэтому соответствующие подробности здесь не приводятся.

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

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

Свет, излученный лазером (21), модулируется импульсным сигналом (22). Оптический усилитель (25) используется для усиления импульсного лазерного пучка, а далее идет полосовой фильтр (26), чтобы отфильтровать шумы Усиленного Спонтанного Излучения, УСИ (АМЕ) усилителя. Далее оптический сигнал поступает в оптический циркулятор (27). На одном из входов циркулятора (27) может быть подключен дополнительный оптический фильтр (28). Свет направляется в волокно

(32) датчика, которое, например, является одномодовым или многомодовым волокном. Некоторая длина этого волокна может быть изолирована и использована в качестве опорного участка (30), например, для размещения в «тихом» месте или для подачи управляемого опорного сигнала. Этот опорный участок (30) может быть сформирован между отражателями - или сочетанием светоделительных элементов и отражателей (29) и (31). Отраженный и обратно рассеянный свет, возникший вдоль волокна (32) датчика, через циркулятор (27) направляется в интерферометр (33).

Поступивший в интерферометр свет усиливается в оптическом усилителе (1) и подается на оптический фильтр (2). Фильтр (2) отфильтровывает лежащие вне полосы шумы УСИ (АМЕ) усилителя (1). Затем свет поступает в оптический циркулятор (3), соединенный с 3×3 оптроном (4). Часть пучка направляется в фотодетектор (12) для отслеживания интенсивности входного пучка света. Другие части пучка направляются в первый и второй оптические тракты (5) и (6), отличающиеся длиной оптического пути. Зеркала Фарадеевского Вращения, ЗФВ (FRM) (7) и (8) отражают свет назад по первому и второму трактам (5) и (6), соответственно. Зеркала фарадеевского вращения обеспечивают компенсацию самопроизвольной поляризации вдоль оптических трактов (5) и (6), так что два частичных световых пучка эффективно интерферируют на каждом из входов 3×3 оптрона (4). Оптрон (4) вводит сдвиги фаз на 0 градусов, +120 градусов и -120 градусов относительно интерференционного сигнала, таким образом создаются первый, второй и третий компоненты интерференционного сигнала с относительными сдвигами по фазе.

Первый и второй компоненты интерференционного сигнала направляются оптроном (4) на фотодетекторы (13) и (14), а третий компонент интерференционного сигнала, падающий на оптический циркулятор (3), направляется на фотодетектор (15).

Фотодетекторы (12), (13), (14) и (15) преобразуют свет в электрические сигналы. Эти электрические сигналы преобразуются в цифровую форму, и затем быстродействующее устройство (34) обработки рассчитывает относительную модуляцию оптической фазы вдоль опорного участка (30) волокна и вдоль волокна (32) датчика. Это устройство обработки синхронизировано по времени с импульсным сигналом (22). Разница длин оптических путей (5) и (6) определяет пространственное разрешение.

Способы расчета относительных фаз и амплитуд трех сдвинутых по фазе компонентов интерференционного сигнала описаны в литературе. Например, в документах Zhiqiang Zhao et al. ("Improved Demodulation Scheme for Fiber Optic Interferometers Using an Asymmetric 3×3 Coupler", J. Lightwave Technology, Vol.13, No.11, November 1997, pp.2059 - 2068) и Huang et al. (US 5946429) раскрыты методики демодуляции выходов 3×3 оптронов в приложениях с мультиплексированием незатухающих волн.

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

Обратимся к ФИГ.3; средство обработки (103) может осуществлять рисование «от руки» геометрических мест точек (104) на графике интенсивности в координатах времени (301) и пространства (302), что может быть полезно для быстрой оценки конкретного акустического поведения источника (100) звука. Однако в одном из предпочтительных вариантов осуществления настоящего изобретения средство обработки (103) может быть реализовано с помощью компьютера и компьютерных программ, отрабатывающих математические алгоритмы.

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

Как показано на ФИГ.2, источник S (200) звука расположен на ортогональном расстоянии d (203) от точки О (201), обозначенной как положение z0. Звуковая волна, распространяющаяся с постоянной скоростью, выйдя из S (200) в момент времени ts, достигнет точки О (201) в момент t0 - до достижения какой-то другой части детектора в момент t, например, точки Q (202) в положении z. Если средняя скорость распространения звуковой волны обозначить как v, то по теореме Пифагора получим выражение:

связывающее положение и время, в которое произвольная отдельная точка фронта сферической звуковой волны, испущенной источником S (200) звука, достигнет произвольной отдельной точки Q (202), относительно точки первого обнаружения О (201). Но время, за которое звуковая волна достигнет точки первого обнаружения О (201), определяется выражением:

Таким образом, комбинируя выражения 1 и 2, получим:

Хотя при определенных условиях может быть возможно определить величины z0 и t0 или ts независимо, но обычно z0 и t0 в точке О (201) или ts в точке S (200), а также расстояние d (203) и скорость v должны считаться неизвестными, получаемыми в средстве обработки (105). Геометрические места точек (104) состоят из пар координат (z, t), так что, используя общеизвестные математические приемы, например, подбор кривой методом наименьших квадратов, специалист в области обработки сигналов может получить наилучшие оценки неизвестных в выражении 1 или выражении 3, получая величину расстояния d (203), что является предпочтительной целью настоящего изобретения. Выражение 1 и выражение 3 основаны на предположении, что скорость распространения звуковой волны постоянна, поэтому любое отклонение замеренных геометрических мест точек (104) от расчетов по выражению 1 или выражению 3 может служить указанием на вариации средней скорости распространения между различными частями геометрических мест точек (104) и источником (100) звука. Если эта средняя скорость вблизи детектора известна, например, из предшествующих измерений, то сравнение величины, полученной средством обработки (105) из выражения 1 или выражения 3, с ранее известным значением может дать поверку качества величины d, одновременно полученной из выражения 1 или выражения 3. Полученные таким образом малые расхождения величин v дадут уверенность в величине d, тогда как большие расхождения величин средней скорости заставят предположить неожиданное поведение и, следовательно, подвергнут сомнению полученную величину d.

Изобретатели обнаружили, что на практике при определенных условиях форма звуковой волны, достигшей акустического детектора (101) и регистрируемой в точке Q (202), может для величин z, больших, чем z0, отличаться от формы волны для величин z, меньших, чем z0. Это различие форм может быть достаточно велико, что позволяет использовать точку различия, в которой происходит изменение, как точку первого поступления О (201), получая таким образом значения t0 и z0. Кроме того, средняя скорость звука v на участке между детектором (101) и источником (100) звука может быть известна благодаря знанию материала или, более точно, из предшествующего измерения. Подставляя выражение 2 в выражение 1, получим:

Таким образом, при известных t0, z0 и v, выражение 4 дает оценку d для каждого измерения расстояния z до геометрического места точек (104) в момент времени t. В этом случае средство обработки (103) для получения геометрических мест точек и средство обработки (105) для получения положения источника звука предпочтительно могут быть скомбинированы, как будет сейчас объяснено.

Выражение 4 можно представить в следующем виде:

Где τ(z) это задержка между моментом, когда звуковая волна достигает точки первого поступления О (201), и положением Q (202). Сигнал (102), полученный от детектора (101), может быть представлен в виде двумерной матрицы замеров A[i, j], в которой строки (индекс i) содержат пространственные замеры (z[i]) в качестве первого размера матрицы, а столбцы (индекс j) содержат временные замеры (t[j]) в качестве второго размера. Вторая двумерная матрица замеров B[i, j, k] может быть получена из A[i, j] сдвигом каждой строки на соответствующую величину, полученную из временной задержки, которую дает соответствующее значение z[i], подставленное в выражение 5 вместе с оценочным значением d[k]. Получив набор матриц В для некоторого диапазона оценочных значений d[k], можно в пространственном измерении определить корреляцию, или функцию стоимости, следующим выражением:

которое равно единице, если для каждого значения j величины B[i, j, k] равны для всех i, и меньше единицы, если не равны. Так как форма обнаруженной звуковой волны принимается существенно неизменной, величина C[k], то есть ее максимум, указывает значение d[k], которое наиболее вероятно, т.е. наилучшим образом описывает измеренную конфигурацию звуковой волны, принятой детектором (101). В тех условиях, когда формы звуковой волны различаются для z>z0 и z<z0, необходимо оценить выражение 6 отдельно для каждой из двух пространственных зон z>z0 и z<z0, а затем два полученных значения наиболее вероятной величины d (203) могут быть скомбинированы для получения одного оптимального значения с использованием взвешенного среднего, где веса даются максимальными значениями функции стоимости C[k] для каждой из зон.

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

Конкретнее, в применениях, относящихся к сфере обеспечения безопасности, распределенный оптоволоконный акустический датчик обычно заглублен примерно на метр под поверхностью земли. В аналогичных устройствах, использующих источник волны сжатия, более двух третей всей излученной энергии преобразуется в волны Релея. Разные частотные компоненты поверхностной волны имеют разные скорости распространения. Это свойство называется дисперсией. Дисперсионные кривые могут быть получены анализом данных, записанных группой обычных сейсмоприемников (см. ФИГ.10, взятую из работы Park CB, Miller RD, Xia J (1999). Multichannel analysis of surface Waves (Masw); Geophys., 64, 800-8).

На ФИГ.11 (взятой из работы Gabriels, P., Snider, R., and Noletm G., 1987, In situ measurements of shear-wave velocity in sediments with higher mode Rayleigh waves. Geophys. Prospecting, 35, 187-196) показан пример первого задокументированного многоканального подхода к исследованию поверхностных волн, относящийся к 1980-м годам, когда голландские исследователи применили 24-канальную систему сбора данных для получения распределения скоростей сдвиговых волн приливно-отливной полосы посредством анализа зарегистрированных поверхностных волн.

Чтобы применить многоканальный анализ для определения положения источника, прежде всего получают две дисперсионных матрицы замеров A[i, j], как описано выше. Т.е. сигнал (102), полученный от детектора (101), может быть представлен в виде двумерной матрицы замеров A[i, j], в которой строки (индекс i) содержат пространственные замеры (z[i]) в качестве первого размера матрицы, а столбцы (индекс j) содержат временные замеры (t[j]) в качестве второго размера. Затем третья двумерная матрица замеров G[i, j, k] может быть выведена из A[i, j] получением частотного спектра временного размера (например, с использованием БПФ (FFT)) и фазовым сдвигом n-го ряда умножением его