СПОСОБ ТЕПЛОВОГО НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ ТЕПЛОТЕХНИЧЕСКОГО СОСТОЯНИЯ ПРОТЯЖЕННЫХ, СЛОЖНОПРОФИЛЬНЫХ И ТРУДНОДОСТУПНЫХ ОБЪЕКТОВ
METHOD FOR THERMAL NONDESTRUCTIVE INSPECTION OF THERMAL-TECHNICAL STATE OF LONG, NON-UNIFORM AND HARD-TO-REACH OBJECTS

Область техники
Изобретение относится к области измерительной техники, в частности к тепловому неразрушающему контролю объектов, и может быть использовано для теплового контроля (мониторинга) технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа на основе обратного рамановского рассеяния в стандартных одномодовых и многомодовых оптических волокнах.
Предшествующий уровень техники
Тепловой неразрушающий контроль включает две основные технологические операции:
1. Регистрацию температурного поля поверхности контролируемого объекта.
2. Обработку температурного поля с целью решения поставленных задач: обнаружение внутренних несплошностей, определения качества теплоизоляции и т.п.
В большинстве случаев для регистрации температурного поля поверхности используется тепловизионная техника (О.Н.Будадин, А.И.Потапов, В.И.Колганов и др. Тепловой неразрушающий контроль изделий. - М., Наука, 2002, 476 с.).
К способам теплового контроля технического состояния протяженных объектов, например трубопроводов, теплотрасс и т.п., использующим тепловизинную технику, относится способ, раскрытый в указанном источнике. Способ заключается в перемещении над земной поверхностью системы регистрации температурного поля - тепловизора, регистрации видеоизображения температурного поля и его компьютерной обработке. По результатам обработки определяют аномалии температурного поля, характерные для дефектных областей теплоизоляции. Имея определенные преимущества, данный способ регистрации температурного поля имеет существенные недостатки, которые значительно ограничивают его применение на практике, а также при решении задач долговременного периодического мониторинга технического состояния объектов. Эти недостатки следующие:
- большое влияние неравномерного качества поверхности, с которой осуществляется регистрация температурного поля, наличие помех, трудностей при перемещении аппаратуры регистрации (тепловизора) и др., которые снижают достоверность контроля,
- низкая производительность контроля вследствие необходимости перемещения тепловизора вдоль теплотрассы,
- большая стоимость контроля в случае применения авиационной техники для перемещения тепловизора,
- трудности дефектометрии трубопроводов, расположенных в полых коробах,
- затруднения постоянного в течение длительного времени контроля одного и того же объекта, особенно с небольшим периодом контроля.
За последние 30 лет оптические волокна (ОВ) стали важнейшим элементом систем связи благодаря очень малым потерям, которые испытывает свет при распространении, большой полосе пропускания для информационного сигнала и нечувствительности к воздействию внешних электромагнитных полей. Для телекоммуникационных применений разработаны различные конструкции оптических кабелей (ОК), способные работать в широком диапазоне внешних температур и механических воздействий, в условиях высокой влажности в химически активных средах.
Современные волоконно-оптические линии связи позволяют передавать широкополосные сигналы (10 ГГц и более) на расстояние в десятки километров без ретрансляторов. Эти же свойства делают ОВ и ОК крайне привлекательными для использования в устройствах контроля и измерения различных физических величин, таких как температура, механические напряжения, давление и т.п. Такие устройства называют волоконно-оптическим датчиками (ВОД).
В общем случае ВОД можно разделить на два больших класса [1]. К первому классу относятся внешние или гибридные датчики. В них ОВ служит средством доставки сигнала к различным чувствительным приборам, в которых пучки света преобразуются внешними воздействиями, подвергаясь модуляции по интенсивности или поляризации.
В другом типе датчиков само оптическое волокно выступает в качестве чувствительного элемента. Изменение параметров внешней среды приводит к изменению параметров сигнала, распространяющегося по ОВ. Такие датчики по самой своей сути являются распределенными, поскольку воздействию подвергаются различные участки ОВ (все ОВ или каждая его точка) и, тем самым, имеется возможность отслеживать изменения среды в большом количестве точек или непрерывно по всей длине волокна. В этом случае одно ОВ заменяет систему точечных датчиков (их число в некоторых случаях может составлять сотни и тысячи штук), которые должны быть установлены для решения конкретной задачи контроля.
В настоящее время самым распространенным способом измерения распределения характеристик ОВ по его длине является оптическая рефлектометрия, когда параметром является мощность сигнала обратного рассеяния, приходящего от различных точек ОВ. Импульсные оптические рефлектометры (OTDR - optical time domain reflectometer), использующие сигнал обратного рэлеевского рассеяния, широко применяются в волоконно-оптической связи для измерения затухания и длины ОВ. При рэлеевском рассеянии частота (длина волны) падающего и рассеянного излучения совпадают.
Измерение температуры ОВ осуществляется путем регистрации бриллюэновского или рамановского рассеяния. Спектр рассеяния в этих случаях состоит из двух компонент: стоксовой и антистоксовой (их длины волн соответственно больше и меньше длины волны падающего излучения).
Бриллюэновское рассеяние возникает в результате взаимодействия между падающим излучением и акустическими волнами, которые образуются вследствие механических и температурных напряжений в ОВ.
Бриллюэновское рассеяние имеет две спектральные компоненты, частота которых отличается от частоты падающего излучения на 10.5…11 ГГц [2]. Температура ОВ непосредственно влияет как раз на величину этого частотного сдвига. Метод бриллюэновской рефлектометрии применим только в одномодовых ОВ.
Рамановское рассеяние возникает в результате взаимодействия между падающим излучением и колебательными состояниями молекул вещества. Интенсивность антистоксовой компоненты рассеяния зависит от температуры, которая, с другой стороны, практически не влияет на интенсивность стоксовой компоненты. Для регистрации этих компонент используется рефлектометрический метод. Как и в стандартном случае, он позволяет получить данные о непрерывном распределении исследуемого параметра (в данном случае температуры) вдоль оптического волокна.
Метод рамановской рефлектометрии позволяет использовать для этой цели как стандартные телекоммуникационные ОВ, так и специализированные оптические волокна и кабели, предназначенные для работы в условиях высоких температур, давлений, в активных средах и т.п. ОВ одновременно выполняет функции сенсора и канала распространения измерительных сигналов, а рамановский рефлектометр - прибора, регистрирующего и обрабатывающего эти сигналы.
При прохождении электромагнитного излучения через физическую среду часть этого излучения рассеивается на неоднородностях и атомах (молекулах) вещества. При этом в рассеянном излучении присутствуют компоненты с той же длиной волны, что и падающее излучение, а так же и компоненты с другими длинами волн. На фиг.1 схематически представлен спектральный состав рассеянного излучения в кварцевом ОВ.
Основные параметры рэлеевского рассеяния - длина волны и интенсивность не зависят от температуры ОВ. С другой стороны, ее изменение приводит к сдвигу частоты бриллюэновских компонент и изменению интенсивности антистоксовой компоненты рамановского рассеяния. Эти эффекты схематически отображены на фиг.2.
В соответствии с указанными физическими явлениями существует два способа измерения температуры ОВ.
При бриллюэновском рассеянии происходит взаимодействие между распространяющимся оптическим сигналом и акустическими волнами, которые возникают вследствие механических и температурных напряжений в ОВ. Эти напряжения создают волны показателя преломления, которые движутся со звуковой скоростью. Вследствие эффекта Доплера у рассеянного излучения появляется сдвиг по длине волны (частоте). Скорость акустической волны непосредственно связана с плотностью среды, которая, в свою очередь, зависит от температуры и механического напряжения. Таким образом, бриллюэновский сдвиг частоты рассеянного сигнала несет в себе информацию и о температуре и о напряжении ОВ.
В этом заключается главная проблема при реализации измерителя температуры ОВ с помощью анализа бриллюэновского рассеяния - невозможность отделить механические и температурные влияния на сдвиг частоты. Кроме того, в измерительной аппаратуре необходимо применять лазерные излучатели с узким спектром излучения и высокодобротные волоконно-оптические фильтры, а для регистрации бриллюэновских компонент использовать либо гетеродинный, либо интерференционный метод приема. Из этого следует, что реализация метода регистрации бриллюэновского рассеяния для измерения температуры возможна только в одномодовых ОВ, а использование сложных высокоточных волоконно-оптических компонентов приводит к значительному удорожанию аппаратной части прибора.
При рамановском рассеянии изменение температуры ОВ приводит к изменению интенсивности (мощности) антистоксовой компоненты. При этом температура не влияет на характеристики стоксовой компоненты и рэлеевского рассеяния. Отношение интенсивностей компонент рамановского рассеяния определяется из выражения [3]:
AS
S
где I и I - интенсивности антистоксовой и стоксовой компонент соответственно;
AS
S
λ и λ - длины волн этих компонент;
h - постоянная Планка;
с - скорость света в вакууме;
Т - абсолютная температура;
ν - величина частотного сдвига между падающим излучением и рамановскими линиями.
-1
В кварцевом ОВ рамановские линии отстоят от частоты падающего излучения на 13.2 ТГц или 440 см [3]. Например, если длина волны падающего излучения 1064 нм, то пики стоксовой и антистоксовой компонент будут приходиться на 1116 нм и 1016 нм соответственно.
AS
S
Из соотношения (1) следует, что для измерения температуры ОВ достаточно знать отношение интенсивностей I/I. При этом расстояние по оси длин волн между компонентами довольно велико. Это позволяет:
- использовать для регистрации рамановского рассеяния такие же методы и устройства, как и в обычной рэлеевской рефлектометрии;
- применять стандартные телекоммуникационные оптические мультиплексоры и фильтры для разделения спектральных компонент рассеяния;
- измерять температуру как одномодовых, так и многомодовых ОВ.
Однако интенсивность антистоксовой компоненты сигнала обратного рассеяния будет зависеть не только от температуры ОВ, но и от затухания в нем. Поэтому аппаратура и методы (методики) измерения должны учитывать этот фактор.
Известен способ теплового контроля технического состояния протяженных объектов, например трубопроводов, теплотрасс и т.п. (О.Н.Будадин, А.И.Потапов, В.И.Колганов и др. Тепловой неразрушающий контроль изделий. - М., Наука, 2002, 476 с.). Он включает регистрацию температурного поля области поверхности земли в районе расположения трубопровода и анализ температурного поля. За счет повышенной температуры трубопровода (теплоноситель, протекающий по трубопроводу имеет повышенную температуру по сравнению с окружающей средой) и нарушенной теплоизоляции либо нарушенной сплошности самого трубопровода, на поверхности земли образуются аномалии температурного поля, характеристики которых содержат информацию о техническом состоянии трубопровода. Анализ температурного поля позволяет оценить характеристики нарушений трубопровода и объективно определить его техническое состояние.
Однако этот способ имеет существенные недостатки:
- необходим свободный доступ к поверхности для обследования тепловизионной системы,
- большие трудности определения толщины и качества технического состояния теплоизоляции из-за неравномерного качества поверхности, наличия помех, трудностей при перемещении аппаратуры регистрации (тепловизора) и т.п.,
- низкая производительность контроля вследствие необходимости перемещения тепловизора вдоль теплотрассы,
- высокая стоимость контроля в случае применения авиационной техники для перемещения тепловизора,
- трудности дефектометрии трубопроводов, расположенных в полых коробах,
- затруднения постоянного в течение длительного времени контроля одного и того же объекта, особенно с небольшим периодом контроля.
Сущность изобретения
В свете изложенного возникает необходимость создания достоверного метода теплового неразрушающего контроля (мониторинга) технического состояния протяженных сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа на основе обратного рамановского рассеяния в стандартных одномодовых и многомодовых оптических волокнах.
Это позволит решать следующие задачи:
1. Повышение объективности, достоверности результатов теплового контроля технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа, в т.ч. обеспечение достоверности результатов контроля не менее 0,98.
2. Расширение области применения для различных классов объектов и условий проведения контроля, в т.ч.
- обеспечение проведения контроля без условий свободного доступа к поверхности земли, где расположен контролируемый объект,
- повышения помехозащищенности за счет исключения влияния качества поверхности земли при подземном расположении трубопровода и т.п.,
- повышение производительности контроля, обеспечения производительности контроля не менее 10 км протяженности трубопровода за 2-3 мин в течение не менее 15-20 лет.
3. Снижение затрат на проведения контроля за счет исключения использования перемещающихся узлов и блоков системы контроля.
4. Возможность проведения теплового контроля качества и дефектометрии трубопроводов, расположенных в полых коробах.
Таким образом, появится возможность решения актуальных для практики задач повышения объективности, достоверности результатов теплового контроля технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа, расширение области его применения, снижение затрат на проведения контроля и снижение потерь энергии за счет своевременного ремонта объекта.
Техническая задача, на решение которой направлено настоящее изобретение, состоит в:
1. обеспечении достоверности результатов контроля, в т.ч. повышении объективности, достоверности результатов контроля технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа,
2. расширении области применения теплового контроля для различных классов объектов и условий проведения контроля, в т.ч.
- обеспечении проведения контроля без условий свободного доступа к поверхности земли, где расположен контролируемый объект,
- повышении помехозащищенности за счет исключения влияния качества поверхности земли при подземном расположении трубопровода и т.п.,
3. повышении производительности контроля, обеспечения производительности контроля не менее 10 км протяженности трубопровода за 2-3 мин в течение не менее 15-20 лет,
4. возможности проведения теплового контроля качества и дефектометрии трубопроводов, расположенных в полых коробах.
внутр
0
Технический результат достигается за счет того, что согласно способу теплового неразрушающего контроля теплотехнического состояния протяженных, сложнопрофильных и труднодоступных объектов прокладывают одномодовую или многомодовую оптоволоконную линию по поверхности теплоизоляции трубопровода, максимальная допустимая температура нагрева которой составляет 200-250°С, имеряют температуру теплопроводной жидкости, протекающей внутри трубопровода в течение времени Т(t), измеряют тепловой поток q характерного качественного участка трубопровода, имеющего качественную, принимаемую за эталонную теплоизоляцию, измеряют температуру T(x,t) в дискретных точках по длине оптоволоконной линии, формируют матрицу измерений, располагая в столбцах значения температуры в одной точке поверхности теплозащитного слоя объекта, измеренные в разное время, а в строках - значения температуры, полученные в разных точках по длине объекта одновременно:
определяют матрицу толщин теплозащитного слоя вдоль контролируемого объекта от времени
внутр
i
0
Элементы матрицы определяются по формуле: δ(j,i)=λ(Т(t)-T(j,i))/q,
по определенному значению толщины теплозащитного слоя определяют техническое состояние трубопровода и величину потерь передаваемой энергии, где λ - коэффициент теплопроводности материала теплоизоляции.
Технический результат достигается также за счет того, что температуру теплоизоляции измеряют прибором, включающим импульсный твердотельный лазер со схемой запуска, оптический демультиплексор, предназначенный для разделения сигналов рэлеевского и рамановского рассеяния; оптические приемники, предназначенные для регистрации сигналов рэлеевского и рамановского рассеяния и сигнала синхронизации, оптические разветвители, предназначенные для направления оптических импульсов в измеряемое оптическое волокно, на оптический приемник сигнала синхронизации и на оптический демультиплексор; катушку оптического волокна, предназначенного для создания опорного оптического сигнала обратного рассеяния; электронный датчик температуры оптического волокна; блок вычислений, предназначенный для первичной обработки измерительной информации; блок управления, предназначенный для управления процессом измерения и накопления информации; блок интерфейсный, предназначенный для связи прибора с персональным компьютером, и импульсный преобразователь напряжения, предназначенный для выработки напряжений для питания блоков прибора.
При этом используют импульсный твердотельный лазер, генерирующий зондирующие импульсы с длиной волны излучения 1064 нм.
Описание фигур чертежей
Сущность изобретения будет далее раскрыта более подробно со ссылками на сопровождающие чертежи, где
фиг.1. показывает спектральный состав рассеянного излучения в кварцевом ОВ,
фиг.2. - сдвиг частоты бриллюэновских компонент и изменения интенсивности антистоксовой компоненты рамановского рассеяния,
на фиг.3. приведена схема укладки опто-волоконной линии вдоль протяженного объекта - теплотрассы,
на фиг.4. показана укладка опто-волоконной линии вдоль трубопровода,
фиг.5. иллюстрирует характер систематической погрешности,
фиг.6 - характер скачка систематической погрешности,
фиг.7 показывает рефлектограммы стоксовой (верхняя линия) и антистоксовой компонент рамановского рассеяния,
фиг.8 - характерный вид динамического диапазона - разности между начальным уровнем обратного рассеяния и максимальным уровнем шума за пределами рефлектограммы,
на фиг.9. приведены графики зависимости амплитуды шума на рефлектограмме от уровня сигнала обратного рассеяния (уровня рефлектограммы) при различных значениях динамического диапазона рефлектометра,
на фиг.10 - зависимости амплитуды шума на рефлектограмме от погрешности измерения температуры от длины многомодового ОВ. Предполагается, что затухание ОВ на длине волны 1064 нм составляет 0,9 дБ/км и чувствительность интенсивности антистоксовой компоненты к изменению температуры - 0,013 дБ/°С,
на фиг.11 показан график зависимости амплитуды антистоксовой компоненты по отношению к стоксовой от температуры,
фиг.12 иллюстрирует шумовую схему оптического приемника,
на фиг.13 приведена структурная схема прибора,
на фиг.14 изображен внешний вид прибора,
на фиг.15 - блок-схема измерительного алгоритма,
фиг.16 иллюстрирует примеры результатов измерений, отображенных в виде графиков, полученных с помощью программного обеспечения (ПО) прибора,
фиг.17 показывает оптическое волокно и участки, температура которых подлежит анализу,
фиг.18 - распределение весовых коэффициентов,
на фиг.19 изображена схема алгоритма мониторинга заданных участков,
на фиг.20 приведена фотография прибора с катушкой опто-волоконного кабеля,
на фиг.21 - блок-схема измерительного алгоритма,
на фиг.22 - результаты мониторинга,
фиг.23 показывает интерфейс оператора, и
фиг.24 - фрагмент видеоизображения температурного поля поверхности земли в районе контролируемого трубопровода.
Описание предпочтительного варианта изобретения
Ниже рассмотрено использование способа для теплового неразрушающего контроля технического состояния протяженного сложнопрофильного и труднодоступного трубопровода в условиях ограниченного доступа на основе обратного рамановского рассеяния в стандартных одномодовых и многомодовых оптических волокнах.
Прокладывают одномодовую или многомодовую опто-волоконную линию по поверхности теплоизоляции трубопровода. Потенциальная температура нагрева опто-волоконной линии не должна превышать максимальной допустимой (200-250 град.С).
Схема укладки опто-волоконной линии вдоль трубопровода теплотрассы представлена на фиг.3.
На фиг.4 приведены фотографии укладки опто-волоконной линии вдоль трубопровода.
внутр
Измеряют температуру теплопроводной жидкости, протекающей внутри трубопровода в течение времени - Т(t).
Измерение температуры осуществляется либо непосредственно самой теплопроводной жидкости внутри трубопровода, путем введения термощупа через специальные технологические отверстия, либо измеряют температуру поверхности трубопровода на участке отсутствия теплоизоляции.
Измерение может осуществляться термодатчиком - Термометр контактный ТК-5.09 производитель ООО «ТехноАс» (г.Коломна) с погрешностью не более ±0,5 град.
0
Измеряют тепловой поток характерного качественного участка трубопровода (участка качественной эталонной теплоизоляции) - q.
Измерения осуществляют датчиком теплового потока - Электронный измеритель плотности тепловых потоков пятиканальный ИТП-МГ4.03 "ПОТОК" и Зонды. Производитель ООО «Стройприбор» (г.Челябинск).
Измеряют температуру в дискретных точках по длине оптоволоконной линии - T(x, t).
Формируют (располагают) результаты измерения по поверхности теплозащитного слоя, таким образом получают матрицу измерений.
Определяют матрицу толщин теплозащитного слоя вдоль контролируемого объекта от времени
внутр
i
0
Элементы матрицы определяют по формуле δ(j,i)=λ(Т(t)-T(j,i))/q.
По определенному значению толщины теплозащитного слоя определяют техническое состояние трубопровода и величину потерь передаваемой энергии, например, следующим образом:
min
Здесь δ - минимальное допустимое значение толщины теплоизоляции, определенное исходя из предельно допустимых потерь передаваемого по трубопроводу энергоносителя.
Значения матриц температуры Т(nxk) и толщины δ(nxk) теплоизоляции сохраняют в базе данных для ретроспективного анализа функционирования контролируемого трубопровода.
Определим погрешность способа.
Исследование величины сигналов обратного рамановского рассеяния
При определении температуры оптического волокна (ОВ) осуществляется регистрация компонент рамановского рассеяния и производится процедура сравнения сигналов. Поскольку предполагается, что характеристики рэлеевского и стоксового сигнала не зависят от температуры, то отличие в их поведении от антистоксового сигнала (изменение наклона антистоксовой рефлектограммы, появление на ней больших неоднородностей и т.п.) трактуется как изменение температуры на соответствующих участках ОВ.
Однако для более точного измерения температуры ОВ необходимо учитывать различия сигналов, которые имеются даже при постоянной температуре вдоль ОВ. На фиг.5 показаны рефлектограммы стоксовой и антистоксовой компонент линии, состоящей из двух ОВ, соединенных с помощью сварного соединения. Температура обоих ОВ одинаковая, однако видно, что наклон стоксовой (верхней) рефлектограммы, имеющей большую длину волны, меньше. Кроме того, затухания в соединении ОВ разные. Если для определения температуры просто вычислить разность между сигналами, то указанные различия приведут к систематическим погрешностям в виде наклона температурного графика и скачка в месте соединения ОВ, как показано на фиг.5, 6.
имп-S
имп-AS
имп
S
AS
При анализе сигналов комбинационного рассеяния необходимо учитывать, что длина волны рассеянного излучения отличается от длины волны зондирующего импульса. Здесь необходимо отметить, что коэффициент обратного рассеяния В (или В) определяют долю мощности оптического импульса (с длиной волны λ), которая преобразуется в соответствующую компоненту сигнала рамановского рассеяния рассеивается (с длиной волны λ или λ).
имп-AS
Параметры стоксового сигнала, как и рэлеевского, не зависят от температуры, а коэффициент антистоксового обратного рассеяния В зависит и от температуры, и от расстояния, поскольку температура может быть различна в разных точках ОВ. Это свойство можно выразить следующим образом:
имп-AS
имп-AS
0
имп-AS
В=В(Т)·В(ΔT(L)),
имп-AS
0
0
где В(Т) - значение коэффициента антистоксового обратного рассеяния, когда ОВ по всей длине имеет некоторую постоянную температуру Т,
имп-AS
В(ΔT(L)) - изменение этого коэффициента при изменении температуры.
В соответствии с (12), величина затухания в соединении ОВ, измеренная рефлектометром на антистоксовой компоненте, определяется выражением (см. (7) и (9)):
1
2
где T и Т - температура 1-го и 2-го ОВ соответственно;
1,2
0
1,2
Т=Т+ΔT.
Из вышеприведенных равенств следует, что измеренная величина затухания в соединении зависит от:
12,имп
- затухания на длине волны зондирующего оптического импульса (β),
12,AS
- затухания на длине волны антистоксовой компоненты (β),
- отношения коэффициентов антистоксового обратного рассеяния волокон при одинаковой температуре;
- изменения этих коэффициентов при изменении температуры.
Информация о температуре ОВ заложена именно в последнем слагаемом равенства, поэтому для ее точного определения необходимо иметь возможность учесть значения других составляющих. Это значит, что необходимо независимыми измерениями определить:
- коэффициенты затухания γ каждого ОВ на длинах волн оптического импульса и антистоксовой компоненты (не зависят от температуры);
- значения затухания β в соединениях ОВ на этих длинах волн (не зависят от температуры);
- составляющую коэффициента антистоксового обратного рассеяния каждого ОВ при некоторой постоянной вдоль ОВ температуре;
- и, наконец, долю мощности антистоксового сигнала в каждой точке ОВ, зависящую от температуры.
Факторы, определяющие точность измерения температуры оптического волокна при регистрации рамановского рассеяния
Мощность сигнала и амплитуда шума.
При определении температуры оптического волокна путем регистрации интенсивности компонент рамановского рассеяния используется метод импульсной рефлектометрии. Сигналы обратного рассеяния имеют малую амплитуду, и их анализ проводится при низком отношении сигнала к шуму (для больших расстояний это отношение может быть меньше 1). Поэтому основными факторами, определяющими погрешность измерения температуры оптического волокна, являются:
- мощность источника излучения, поскольку мощность сигнала обратного рассеяния (и рэлеевского и рамановского) прямо пропорциональна мощности зондирующего оптического импульса;
- длительность зондирующего оптического импульса, поскольку мощность сигнала обратного рассеяния (и рэлеевского и рамановского) прямо пропорциональна этой длительности;
- собственный шум оптического приемника;
- время измерения (накопления) сигнала.
Для численных оценок влияния указанных факторов используются формула (2) и выражения:
ш
1
2
1/2
ΔU=5·lg(t/t),
ш
1
2
где ΔU (дБ) - изменение амплитуды шума оптического приемника при изменении времени измерения в (t/t) раз.
Из формулы (2) следует, что при увеличении мощности импульса в 2 раза уровень рефлектограммы возрастает на
имп2
имп1
5·lg(Р/Р)=5·lg(2)=1,5 дБ
Аналогично из формулы (16) следует, что при увеличении времени измерения в 2 раза амплитуда шума уменьшается на 0,75 дБ.
интенсивности
Мощность сигнала обратного рассеяния и собственные шумы оптического приемника рефлектометра определяют отношение сигнала к шуму и, значит, амплитуду шума на рефлектограмме. Поскольку измерение температуры ОВ основано на измерении  антистоксовой компоненты рамановского рассеяния, то амплитуда шума на рефлектограмме определяет минимально возможную погрешность определения температуры.
На фиг.7 показаны рефлектограммы стоксовой (верхняя линия) и антистоксовой компонент рамановского рассеяния.
Амплитуда шума на антистоксовой рефлектограмме значительно выше; она составляет примерно 0,022 дБ. Температурная чувствительность интенсивности антистоксовой компоненты составляет 0,45…0,65%/°С [6] или 0,01…0,014 дБ/°С. Это значит, что для данного измерения минимальная погрешность измерения температуры может составлять 1,6…2,2°С.
Поэтому анализ способов повышения точности измерения температуры, прежде всего, должен основываться на анализе отношения сигнала к шуму в рефлектометре и путях его увеличения.
Влияние динамического диапазона рефлектометра на амплитуду шума на рефлектограмме
Теоретическую оценку амплитуды шума на рефлектограмме удобно провести пользуясь понятием динамического диапазона рефлектометра.
Динамический диапазон - это разность между начальным уровнем обратного рассеяния и максимальным уровнем шума за пределами рефлектограммы (см. фиг.8); указанная разность определяется при определенной длительности импульса и заданном времени измерения.
Определим амплитуду шума на рефлектограмме.
Обозначим
c
U - напряжение сигнала обратного рассеяния, соответствующего некоторой точке ОВ, на выходе оптического приемника рефлектометра, В;
ш
U - напряжение шума на выходе оптического приемника рефлектометра, В;
0
U - максимальное напряжение сигнала обратного рассеяния, соответствующего начальной точке ОВ, на выходе оптического приемника рефлектометра, В;
c
А - уровень на рефлектограмме, соответствующий напряжению сигнала U, дБ.
0
0
A - максимальный уровень рефлектограммы, соответствующий напряжению сигнала U, дБ.
Тогда исходя из определения, динамический диапазон
0
ш
D=5·lg(U/U).
Кроме того, поскольку максимальному уровню рефлектограммы в соответствии со стандартом [7] обычно приписывается значение 0 дБ, то все более низкие уровни будут иметь отрицательные значения, и можно записать
c
0
A=5·lg(U/U).
Амплитуду шума в точке на рефлектограмме, соответствующей уровню А, можно определить следующим образом:
Необходимо еще раз отметить, что уровень точки рефлектограммы А - отрицательная величина, а динамический диапазон D - положительная.
На фиг.9 показаны графики зависимости амплитуды шума на рефлектограмме от уровня сигнала обратного рассеяния (уровня рефлектограммы) при различных значениях динамического диапазона рефлектометра.
На фиг.10 эти данные преобразованы в зависимость погрешности измерения температуры от длины многомодового ОВ. При этом предполагается, что затухание ОВ на длине волны 1064 нм составляет 0,9 дБ/км и чувствительность интенсивности антистоксовой компоненты к изменению температуры - 0,013 дБ/°С.
Полученные оценки позволяют конкретизировать требования к величине динамического диапазона рефлектометра, который определяет минимальную достижимую погрешность измерения температуры ОВ.
Мощность сигналов обратного рамановского рассеяния
S
AS
Мощность стоксовой P и антистоксовой P компонент сигнала, рассеянного от участка ОВ длиной dL, находящегося на расстоянии L от начала ОВ, можно определить по формуле [8]:
S(AS)
S(AS)
S(AS)
dP=P(L)·П·Г·dL,
где P(L) - мощность падающего излучения в точке L;
S(AS)
П - коэффициенты Бозе-Эйнштейна распределения плотностей вероятностей фононов для стоксовой (S) и антистоксовой (AS) компонент;
S(AS)
Г - коэффициент преобразования падающей оптической мощности в сигнал обратного рассеяния
Согласно [8]
,
,
где ΔЕ - энергия собственных колебаний кристаллической решетки твердого вещества, для кварцевого стекла ΔЕ=50 мэВ;
k - постоянная Больцмана;
Т - абсолютная температура.
Длина участка рассеяния связана с длительностью зондирующего оптического импульса соотношением [4]:
имп
dL=(τ·с)/(2·n),
где с - скорость света в вакууме;
n - показатель преломления сердцевины ОВ.
S
AS
Используя данные, приведенные в [8], можно определить, что в многомодовом ОВ при комнатной температуре (T=300 K) при длине волны зондирующего излучения 1064 нм и длительности оптического импульса 1 нс мощности стоксовой и антистоксовой компонент, приведенные к входной мощности, составят P=4,67 нВт/Вт, P=0,85 нВт/Вт и стоксова компонента больше антистоксовой в 5,5 раз.
Из этих данных следует, что для длины волны 1064 нм и длительности импульса 1 нс коэффициент обратного рассеяния для стоксовой компоненты будет составлять -83.3 дБ; а для антистоксовой -90.7 дБ. Коэффициент обратного рэлеевского рассеяния для стандартного ММ ОВ составляет -68 дБ для длины волны 850 нм и -76 дБ для длины волны 1300 нм [9].
На фиг.11. приведен график зависимости амплитуды антистоксовой компоненты по отношению к стоксовой от температуры. Видно, что в диапазоне 100 градусов отношение интенсивностей изменяется примерно в 2 раза.
Метод корреляционного рефлектометра
Как указывалось выше, мощность сигнала обратного рассеяния, а значит, и динамический диапазон, прямо пропорциональны мощности и длительности зондирующего оптического импульса, т.е. его энергии. В работе [6] предложено для увеличения энергии полезного сигнала использовать не одиночные оптические импульсы, псевдослучайную последовательность, образованную с использованием комплиментарно-корреляционного кодирования. Принцип использования такого кодирования в оптической рефлектометрии был предложен в статье [10]. Главное отличие комплиментарных кодов от других состоит в том, что теоретически корреляционная функция последовательности таких импульсов не имеет боковых лепестков, вся ее энергия сосредоточена в центральном пике. В [10] получены выражения для отношения сигнала к шуму для обычного одноимпульсного метода и при использовании кодированной последовательности:
В этих формулах
N - число усреднений отсчета рефлектограммы;
М - длина (количество бит) псевдослучайной последовательности;
βh - энергия одиночного импульса с учетом формы импульсного отклика оптического приемника рефлектометра;
σ - амплитуда шума на выходе оптического приемника рефлектометра.
Из этих равенств видно, что кодирование улучшает отношение сигнала к шуму пропорционально квадратному корню из длины кодированной последовательности.
В статье [8] описаны экспериментальные результаты измерения температуры многомодового ОВ. В качестве источника излучения рефлектометра использовался полупроводниковый лазер с внешним резонатором. Длина волны лазера 1064 нм, мощность излучения ~20 мВт. С помощью такого прибора было получено разрешение около 0,3°С при измерении ОВ длиной 4 км в течение 10 минут.
Чувствительность оптического приемника
Оптимизация чувствительности оптического приемника рефлектометра является важнейшей задачей при проектировании прибора.
На фиг.12 изображена шумовая схема оптического приемника.
Здесь обозначены:
- А - усилитель;
- R и С - нагрузочный резистор и паразитная емкость входной цепи;
- is - источник тока сигнала, связан с оптической мощностью Р через чувствительность фотодиода S;
- in - источник шумов, обусловленный паразитными токами во входной цепи ОПр, дробовым шумом фотодиода и тепловым шумом нагрузочного резистора R;
- еа - шум активного усилительного элемента входной цепи.
2
Суммарная мощность шума St (А), приведенная к входу оптического приемника запишется [11]:
где z - комплексное сопротивление входной цепи оптического приемника:
При дальнейших расчетах будем предполагать, что все источники шумов являются "белыми".
Для достижения максимального динамического диапазона в оптическом приемнике рефлектометра используется лавинный фотодиод (ЛФД). Поэтому можно считать, что спектральная плотность мощности Sd обусловлена темновым током ЛФД id и дробовым шумом, вызываемым падающим на ЛФД сигналом с оптической мощностью Р:
2+x
Sd=2·e·(In+M(Im+S·P)),
где е - заряд электрона;
In - неумножаемый темновой ток ЛФД;
М - коэффициент лавинного усиления фотодиода;
х - коэффициент избыточного шума ЛФД (для германиевых ЛФД х=1, для ЛФД на основе InGaAs х=0.8);
Im - умножаемый темновой ток ЛФД.
После интегрирования получаем:
где Ia - входной ток усилителя;
В - полоса пропускания ОПр по уровню 0.5.
Проведем расчет порогового значения чувствительности Pn оптического приемника, разработанного для измерителя температуры со следующими параметрами:
1/2
R=10 кОм; С=3 пФ; In=10 нА; Im=2 нА; Ia=750 нА; М=30; ea=8n В/(Гц); х=0.8; B=15 МГц.
-16
2
При подстановке этих значений величина мощности шума, приведенного ко входу оптического приемника, получается St=1.95·10 A. Отсюда следует, что значение пороговой чувствительности Pn=2.8 нВт, т.е. при таком значении оптической мощности отношение сигнал/шум на выходе оптического приемника составит 1.
Расчеты, проведенные выше, позволяют определить оптимальные значения основных характеристик прибора - динамического диапазона и чувствительности оптического приемника, что даст возможность повысить точность измерения температуры ОВ.
Структурная схема и конструкция прибора
В результате проведенных исследований был разработан прибор для измерения температуры ММ ОВ.
Структурная схема прибора представлена на фиг.13. Прибор содержит:
- импульсный твердотельный лазер 1 со схемой запуска, предназначенный для генерации зондирующих импульсов с длиной волны излучения 1064 нм;
- оптические разветвители (ОР-1 и ОР-2) 2, 3, предназначенные для направления оптических импульсов в измеряемое оптическое волокно, на оптический приемник синхронизации и на оптический демультиплексор;
- оптический демультиплексор 4, предназначенный для разделения сигналов рэлеевского и рамановского рассеяния;
- оптические приемники 5, 6 и 7, предназначенные для регистрации сигналов рэлеевского и рамановского рассеяния (ОПр-Р, ОПр-С и ОПр-А) и сигнала синхронизации (ОПр-синхро);
- катушку оптического волокна ОВ-Т 8, предназначенного для создания опорного оптического сигнала обратного рассеяния;
- электронный датчик температуры 9 этого оптического волокна;
- блок вычислений 10, предназначенный для первичной обработки измерительной информации;
- блок управления 11, предназначенный для управления процессом измерения и накопления информации;
- блок интерфейсный 12, предназначенный для связи прибора с персональным компьютером;
- импульсный преобразователь напряжения (ИПН) 13, предназначенный для выработки необходимых напряжений для питания блоков прибора.
Прибор работает следующим образом.
Лазер 1 генерирует короткие оптические импульсы лазера с частотой, обеспечивающей измерение ОВ длиной не менее 10 км и определяемой параметрами схемы запуска. Основная мощность лазерных импульсов через оптический разветвитель ОР-1 2 с коэффициентом деления 99:1 поступает на выходной разветвитель ОР-2 3 с коэффициентом деления 50:50, а малая ее часть - на оптический приемник синхронизации ОПр-синхро 5. Выходной электрический сигнал этого приемника 5 служит для синхронизации работы остальных блоков прибора.
К выходу разветвителя ОР-2 3 присоединено опорное оптическое волокно ОВ-Т 14 длиной около 300 м. Оно помещено в теплоизолированный кожух, в котором также находится электронный датчик температуры 9. Датчик 9 измеряет температуру ОВ-Т 14, и это значение используется для привязки шкалы измерителя к абсолютным значениям температуры.
Выход ОВ-Т 14 является выходом прибора, к нему подключается измеряемое ОВ.
Сигналы обратного рассеяния измеряемого ОВ возвращаются в измеритель и через разветвитель ОР-2 3 поступают на вход оптического демультиплексора 4. В нем происходит разделение сигналов рэлеевского и рамановского рассеяния. Эти сигналы поступают на соответствующие оптические приемники 5-7, преобразуются в электрическую форму и усиливаются. Затем они переводятся в цифровую форму и дополнительно обрабатываются (усредняются, логарифмируются и т.п.) в блоке вычислительном 10.
Процессор (блок управления 11) управляет работой измерителя. Данные передаются в персональный компьютер, запоминаются в нем и отображаются на экране.
Внешний вид прибора приведен на фиг.14.
На передней панели измерителя расположены:
- выключатель 15 и индикаторный светодиод "Сеть" 16;
- оптический разъем типа FC/APC 17 для подключения измеряемого ОВ;
- индикаторный светодиод "Лазер" 18.
Оптический разъем закрыт защитным колпачком 19.
На задней панели измерителя расположены:
- разъем 20 для подключения сетевого кабеля питания "~230 В";
- разъем 21 USB для подключения кабеля интерфейсного для связи прибора с персональным компьютером;
- гнезда предохранителей 22;
- клемма заземления 23.
Для обеспечения возможности выполнения измерений с прибором было разработано программное обеспечение (ПО), позволяющее управлять измерительным модулем с персонального компьютера под управлением ОС Windows XP через USB интерфейс, основываясь на низкоуровневой системе команд. Система команд позволяет устанавливать необходимые параметры измерения и считывать результаты в режиме реального времени. В результате выполнения алгоритма согласно методике измерений (см. фиг.15) на экране отображается график зависимости температуры от расстояния (см. фиг.16), построенный по точкам с установленным пространственным разрешением.
Согласно своему назначению, разрабатываемый продукт должен иметь возможность диагностики большого количества потенциально опасных объектов, расположенных вдоль многомодового оптического волокна. После того как получены положительные результаты в измерении температуры непрерывно вдоль всего оптического волокна (в произвольном (точечном) месте, часть 1), необходимо было разработать методику и программный алгоритм, позволяющий контролировать параметры дискретных объектов (участков) (заранее известных и заданных для температурного анализа) - см. фиг.17.
Для того чтобы программа выводила значения температуры (максимальной, минимальной или некоторой средней расчетной) сразу в цифровом виде (не аналоговом, т.е. просто в виде графика), необходимо указать участки, на которые воздействуют температурные поля интересующих объектов. В этом случае непрерывный мониторинг объектов может осуществляться оптимальным образом.
Ввиду того, что воздействие температурного поля объекта по длине участка не равномерно для расчета весовых коэффициентов значений температуры для каждой точки участка, применяется распределение Гаусса (см. фиг.18)
Тогда распределение температуры по длине ОВ будет описываться выражением
Поэтому алгоритм мониторинга участков может иметь вид, показанный на фиг.19.
Основные технические характеристики прибора
Режимы работы:
- анализатор температуры;
- рефлектометр.
Длина волны зондирующего оптического импульса - (1064±10) нм.
Тип оптического разъема - FC/APC.
Максимальная длина измеряемого OB - 10 км.
Технические характеристики в режиме анализатора температуры.
Диапазон измеряемых температур - от 0°С до плюс 250°С (диапазон температур определяется собственной температурной прочностью оптических волокон и кабелей).
Пределы допускаемой абсолютной погрешности измерения температуры MM OB в рабочих условиях применения соответствуют таблице 1.
Таблица 1
Максимальное расстояние при измерении температуры, км
Погрешность измерения температуры, °С
1
0,5
3
0,8
10
1,5
Пределы допускаемой абсолютной погрешности определения координаты до точки измерения температуры - не более ±3,0 м.
Технические характеристики в режиме рефлектометра.
Диапазон измеряемых расстояний - 12 км.
Длительность зондирующих оптических импульсов - 3 нс.
Пределы допускаемой абсолютной погрешности измерения расстояния составляют:
-5
ΔL=±(dL+3·10·L),
где dL - разрешение (интервал дискретизации сигнала обратного рассеяния), равное 2,5 м;
L - измеряемое расстояние, м.
Динамический диапазон измерения затухания в дБ при отношении сигнал/шум, равном 1 (ОСШ=1) и числе усреднений 128 соответствует таблице 2.
Таблица 2.
Динамический диапазон, дБ
Сигнал рэлеевского рассеяния
не менее 27
Сигнал стоксовой компоненты рамановского рассеяния
не менее 23
Сигнал антистоксовой компоненты рамановского рассеяния
не менее 24
Нелинейность шкалы затухания должна быть не более ±0,05 дБ/дБ.
Питание измерителя осуществляется от сети переменного тока с напряжением (230±23) В и частотой (50,0±0,4) Гц.
Габариты измерителя - 478×362×107 мм.
Масса измерителя - 5 кг.
На фиг.20 приведена фотография прибора с катушкой опто-волоконного кабеля.
На фиг.21 приведена блок-схема измерительного алгоритма.
Технические и эксплуатационные характеристики разработанной аппаратуры и предлагаемого способа контроля приведены в таблице 3.
Таблица 3.
Эксплуатационно-технические характеристики
п/п
Наименование характеристики
Значение
2
3
Количество обслуживающего персонала
1 человек (оператор, предусмотрен режим автоматического мониторинга)
Общий вес технических средств
Вес электронного блока до 5 кг
Кабель-сенсор длиной до 15 км, весом не более 30 кг/км
Контролируемые и диагностируемые объекты
Шахты, тоннели, склады, трубопроводы, нефте- и газохранилища, энергетическое и электрическое оборудование, промышленное оборудование, объекты природопользования и экологии, строительные объекты, инженерные коммуникации и сооружения, потенциально опасное оборудование
Вид регистрируемой
- Температурное распределение по кабелю-сенсору,
п/п
Наименование характеристики
Значение
2
3
информации в полевых условиях (условиях реальной эксплуатации объекта без вывода его из режима эксплуатации)
температурные поля,
- Места утечки теплоносителей
- температурные параметры в контрольных (реперных) точках.
Результаты контроля
- реальное техническое состояние и остаточный ресурс потенциально опасных объектов
- качество (отсутствие дефектов сплошности и структуры) материалов,
- безопасность и надежность конструкций,
- энергоэффективность конструкций
Экспериментальные исследования способа контроля с использованием разработанной аппаратуры проводились на участке теплотрассы (трубопровода) протяженностью 8 км.
В соответствии с предлагаемым способом и разработанной аппаратурой проводился мониторинг технического состояния теплоизоляции трубопровода в течение 3 лет.
Одновременно один раз в год осуществлялся тепловой контроль и диагностика технического состояния трубопровода в соответствии со способом, принятым в качестве прототипа. В качестве измерительной аппаратуры использовался тепловизор фирмы FLIR, размещенный на вертолете. Регистрация результатов в виде видеоизображений температурных полей (в режиме видеофильма) осуществлялась в память ЭВМ.
Результаты контроля качества теплоизоляции по заявляемому способу в виде температурного профиля на участке трубопровода протяженностью 450 м приведены на фиг.22. Из приведенного чертежа ясно видны три дефекта нарушения качества теплоизоляции, которые вызвали локальное повышение температуры на 8-10 град.
На фиг.23 приведен интерфейс оператора в процессе проведения экспериментальных исследований по заявляемому способу, который в оперативном режиме позволяет провести экспресс-анализ дефектных участков.
На фиг.24, в качестве примера, приведен фрагмент видеоизображения температурного поля поверхности земли в районе контролируемого трубопровода, полученного в результате теплового контроля по способу, принятому в качестве прототипа.
По результатам контроля определяют ориентировочную толщину теплозащитного слоя вдоль контролируемого объекта от времени по формуле
внутр
0
δ(x,t)=X(Т(t)-T(x,t))/q
По определенному значению толщины теплозащитного слоя определяют техническое состояние трубопровода и величину потерь передаваемой энергии.
min
Здесь δ - минимальная допустимая толщина слоя теплоизоляции на поверхности трубопровода.
Величина повышенных (сверхнормативных) теплопотерь определяется следующим образом:
внутр
min
внутр
Δq=λ[(Т(t)-T(x,t))/δ-(Т(t)-T(x,t))/δ(х,t)].
Результаты сравнительного анализа результатов экспериментальных исследований предлагаемого способа и способа, принятого в качестве прототипа, на одном и том же объекте приведены в таблице 4.
Экспериментальные исследования показали, что предлагаемый способ теплового неразрушающего контроля (мониторинга) технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа на основе обратного рамановского рассеяния в стандартных одномодовых и многомодовых оптических волокнах обеспечивает:
1. Повышение объективности, достоверности результатов теплового контроля технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа, в т.ч. обеспечение достоверности результатов контроля не менее 0,98.
2. Расширение области применения для различных классов объектов и условий проведения контроля, в т.ч.
- обеспечение проведения контроля без условий свободного доступа к поверхности земли, где расположен контролируемый объект,
- повышения помехозащищенности за счет исключения влияния качества поверхности земли при подземном расположении трубопровода и т.п.,
- повышение производительности контроля, обеспечения производительности контроля не менее 10 км протяженности трубопровода за 2-3 мин в течение не менее 15-20 лет.
3. Снижение затрат на проведения контроля за счет исключения использования перемещающихся узлов и блоков системы контроля.
4. Возможность проведения теплового контроля качества и дефектометрии трубопроводов, расположенных в полых коробах.
5. Значительно большую оперативность и производительность контроля.
6. В десятки раз меньшую стоимость проведения контроля.
7. Возможность выявления меньших по размеру дефектов (ориентировочно, в 10 раз).
Изобретение имеет следующие преимущества:
1. Повышение объективности, достоверности результатов теплового контроля технического состояния протяженных, сложнопрофильных и труднодоступных объектов (например, трубопроводов и теплотрасс) в условиях ограниченного доступа, в т.ч. обеспечение достоверности результатов контроля не менее 0,98.
2. Расширение области применения для различных классов объектов и условий проведения контроля, в т.ч.
- обеспечение проведения контроля без условий свободного доступа к поверхности земли, где расположен контролируемый объект,
- повышения помехозащищенности за счет исключения влияния качества поверхности земли при подземном расположении трубопровода и т.п.,
- повышение производительности контроля, обеспечения производительности контроля не менее 10 км протяженности трубопровода за 2-3 мин в течение не менее 15-20 лет.
3. Снижение затрат на проведения контроля за счет исключения использования перемещающихся узлов и блоков системы контроля.
4. Возможность проведения теплового контроля качества и дефектометрии трубопроводов, расположенных в полых коробах.
5. Значительно большую оперативность и производительность контроля.
6. В десятки раз меньшую стоимость проведения контроля.
7. Возможность выявления меньших по размеру дефектов (ориентировочно, в 10 раз).
ИСТОЧНИКИ ИНФОРМАЦИИ
1. Волоконно-оптические датчики. Вводный курс для инженеров и научных работников: Пер.с англ. / Под. ред. Э.Удда. - М.: Техносфера, 2008. - 520 с.
2. Листвин А.В., Листвин В.Н. Рефлектометрия оптических волокон. - М.: ЛЕСАРарт, 2005. - 208 с.
3. Wait P.C., Souza K De, Newson T.P. A theoretical comparison of spontaneous Raman and Brillouin based fibre optic distributed temperature sensors Opt. Commun. 1997, V.144 (1-3), p.17-23.
4. Григорьянц В.В., Чаморовский Ю.К. Диагностика волоконных световодов и оптических кабелей методом обратного рассеяния // Итоги науки и техники, ВИНИТИ, серия "Радиотехника", 1982 г., т.29, с.47-79.
5. Hartog A.H., Gold M.P. On the theory of backscattering in single-mode optical fibers - Journal of Lightwave Technology, 1984, V.LT-2, №.2, p.76-82.
6. Soto M.A., Sahu P.K., Faralli S. et al. High performance and highly reliable Raman-based distributed temperature sensors based on correlation-coded OTDR and multimode graded-index fibers - Proc. of SPIE, 2007, V.6619, p.66193B-1-66193B-4.
7. GR-196-Core. Generic requirements for optical time domain reflectometer (OTDR) - type equipment.: Bellcore, December 1998.
8. Farahani M.A., Gogolla T. Spontaneous Raman Scattering in Optical Fibers with Modulated Probe Light for Distributed Temperature Raman Remote Sensing - Journal of Lightwave Technologies, 1999, V.17, №.8, p.1379-1391.
9. Corning InfiniCor 50 µm optical fibers. Product information. http://www.corning.com/docs/opticalfiber/pil457.pdf
10. Nazarathy M., Newton S.A., Giffard R.F. et al. Real-time long range complementary correlation optical time domain reflectometer - Journal of Lightwave Technology, 1989, V.7, №1, p.24-37.
11. Personic S.D. Receiver design for digital fiber optic communication systems - The Bell System Technical Journal, 1973, V.52, No.6.