CN112219361A - 入侵检测的事件统计生成方法及装置 - Google Patents

入侵检测的事件统计生成方法及装置 Download PDF

Info

Publication number
CN112219361A
CN112219361A CN201980037983.7A CN201980037983A CN112219361A CN 112219361 A CN112219361 A CN 112219361A CN 201980037983 A CN201980037983 A CN 201980037983A CN 112219361 A CN112219361 A CN 112219361A
Authority
CN
China
Prior art keywords
power region
signal power
event
signal
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201980037983.7A
Other languages
English (en)
Inventor
S·穆罕默德
A·米汉
B·博尔克
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ava Risk Group Ltd
Original Assignee
Ava Risk Group Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from AU2018901143A external-priority patent/AU2018901143A0/en
Application filed by Ava Risk Group Ltd filed Critical Ava Risk Group Ltd
Publication of CN112219361A publication Critical patent/CN112219361A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00—Burglar, theft or intruder alarms
    • G08B13/02—Mechanical actuation
    • G08B13/14—Mechanical actuation by lifting or attempted removal of hand-portable articles
    • G08B13/1481—Mechanical actuation by lifting or attempted removal of hand-portable articles with optical detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
    • G01D5/35354—Sensor working in reflection
    • G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
    • G01D5/35354—Sensor working in reflection
    • G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
    • G01D5/35361—Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/30—Testing of optical devices, constituted by fibre optics or optical waveguides
    • G01M11/31—Testing of optical devices, constituted by fibre optics or optical waveguides with a light emitter and a light receiver being disposed at the same side of a fibre or waveguide end-face, e.g. reflectometers
    • G01M11/3109—Reflectometers detecting the back-scattered light in the time-domain, e.g. OTDR
    • G01M11/3118—Reflectometers detecting the back-scattered light in the time-domain, e.g. OTDR using coded light-pulse sequences
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00—Burglar, theft or intruder alarms
    • G08B13/16—Actuation by interference with mechanical vibrations in air or other fluid
    • G08B13/1654—Actuation by interference with mechanical vibrations in air or other fluid using passive vibration detection systems
    • G08B13/1663—Actuation by interference with mechanical vibrations in air or other fluid using passive vibration detection systems using seismic sensing means
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00—Burglar, theft or intruder alarms
    • G08B13/18—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength
    • G08B13/181—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems
    • G08B13/183—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems by interruption of a radiation beam or barrier
    • G08B13/186—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems by interruption of a radiation beam or barrier using light guides, e.g. optical fibres
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00—Burglar, theft or intruder alarms
    • G08B13/18—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength
    • G08B13/181—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems
    • G08B13/187—Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength using active radiation detection systems by interference of a radiation field
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B29/00—Checking or monitoring of signalling or alarm systems; Prevention or correction of operating errors, e.g. preventing unauthorised operation
    • G08B29/18—Prevention or correction of operating errors
    • G08B29/185—Signal analysis techniques for reducing or preventing false alarms or for enhancing the reliability of the system
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/071—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using a reflected signal, e.g. using optical time domain reflectometers [OTDR]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
    • H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
    • H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
    • H04B10/0799—Monitoring line transmitter or line receiver equipment
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Computer Security & Cryptography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Burglar Alarm Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种用于入侵检测的计算机实施的事件统计生成方法，该方法包括针对多个传感器仓中的每一者，将来自相干光学时域反射仪的多个返回信号处理成时域信号，所述多个返回信号对应于在一个时间段内注入光学传感器光纤中的多个刺激脉冲。对于每个传感器仓：将所述各个时域信号变换成相应的频域信号；根据所述各个频域信号计算预期包含系统噪声的第一频带的第一信号功率区域和预期包含与至少第一事件相关的任何能量的第二频带的第二信号功率区域；以及通过将所述第二信号功率区域除以所述第一信号功率区域，至少部分地产生与所述第二信号功率区域与所述第一信号功率区域之比成比例的事件统计。

Description

入侵检测的事件统计生成方法及装置
技术领域
本发明涉及一种用于入侵检测的计算机实施的事件统计生成方法和事件统计生成装置。
背景技术
所有感测系统面临的挑战之一是要能在一些恶劣的环境中运行。通常安装在室外环境中的入侵检测系统也不例外。这些系统通常需要在大风或大雨期间运行，或靠近附近的交通车道或靠近其他与入侵无关的干扰。
一些入侵检测系统使用相干光时域反射仪(COTDR)技术从埋入式光纤获取信号。大多数埋入式COTDR光纤入侵检测系统在靠近滋扰警报源的位置运行，这些滋扰警报源通常包括来自公路或铁路交叉口的交通以及附近的挖掘设备。这种入侵检测系统的有效性取决于其抑制这些滋扰事件引起的任何警报的能力。
发明内容
本发明的实施方式提供一种用于入侵检测的计算机实施的事件统计生成方法，该方法包括：
针对多个传感器仓中的每一者，将来自相干光学时域反射仪的多个返回信号处理成时域信号，所述多个返回信号对应于在一个时间段内注入到光学传感器光纤中的多个刺激脉冲；以及
对于每个传感器仓：
将各个时域信号变换成相应的频域信号；
从所述各个频域信号计算预期包含系统噪声的第一频带的第一信号功率区域；
从频域信号计算预期包含与至少第一事件相关的任何能量的第二频带的第二信号功率区域；以及
通过将所述第二信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第二信号功率区域与所述第一信号功率区域之比成比例的事件统计。
在一个实施方式中，该方法包括针对多个时间段中的每一个时间段处理多个返回信号，以针对每个传感器仓获得多个事件统计。
在一个实施方式中，该方法包括将每个事件统计转换成指示所述事件统计相对大小的视觉表示，并在显示器上输出所述各个视觉表示的矩阵。
在一个实施方式中，每个视觉表示是灰度值。
在一个实施方式中，每个视觉表示是色阶值。
在一个实施方式中，生成事件统计包括应用从所述第一和第二频带中的频率仓的数量导出的缩放因子。
在一个实施方式中，所述方法包括从所述频域信号计算预期包含与至少第二事件相关的任何能量的第三频带的第三信号功率区域，并通过将所述第三信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第三信号功率区域与所述第一信号功率区域之比成比例的进一步的事件统计。
在另一实施方式中，提供一种事件统计生成装置，该装置包括：
处理器；和
存储器，其存储指令，当所述处理器执行所述指令时，所述指令使所述处理器：
针对多个传感器仓中的每一者，将来自相干光学时域反射仪的多个返回信号处理成时域信号，所述多个返回信号对应于在一个时间段内注入到光学传感器光纤中的多个刺激脉冲；以及
对于每个传感器仓：
将各个时域信号变换成相应的频域信号；
从各个频域信号计算预期包含系统噪声的第一频带的第一信号功率区域；
从所述频域信号计算预期包含与至少第一事件相关的任何能量的第二频带的第二信号功率区域；以及
通过将所述第二信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第二信号功率区域与所述第一信号功率区域之比成比例的事件统计。
在一个实施方式中，当所述处理器执行所述指令时，所述指令使所述处理器针对多个时间段中的每一个时间段处理多个返回信号，以针对每个传感器仓获得多个事件统计。
在一个实施方式中，当所述处理器执行所述指令时，所述指令使所述处理器将每个事件统计转换成指示所述事件统计相对大小的视觉表示，并在显示器上输出所述各个视觉表示的矩阵。
在一个实施方式中，每个视觉表示是灰度值。
在一个实施方式中，每个视觉表示是色阶值。
在一个实施方式中，所述处理器生成事件统计包括应用从所述第一和第二频带中的频率仓的数量导出的缩放因子。
在一个实施方式中，当所述处理器执行所述指令时，所述指令使所述处理器从所述频域信号计算预期包含与至少第二事件相关的任何能量的第三频带的第三信号功率区域，并通过将所述第三信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第三信号功率区域与所述第一信号功率区域之比成比例的进一步的事件统计。
在一个实施方式中，该装置包括光纤传感器，并且所述处理器构成实施实施相干光时域反射仪的控制器的一部分。
另一个实施方式提供了一种采用上述事件统计生成方法的计算机实施的入侵检测方法。
附图的简要说明
图示将结合以下图示描述本发明的实施方式，在图示中：
图1为示意图，其示出需要区分滋扰与由隐蔽埋入的传感器检测到的入侵振动；
图2示出了相干光学时域反射仪感测布置的实施例；
图3是来自图2的感测布置的原始拍摄数据的实施例；
图4是示出图3的信号如何被过滤的示意图；
图5示出了在隐蔽埋入的传感器上垂直行走的实施例；
图6示出了图5中的信号的功率谱密度；
图7示出了在隐蔽埋入的传感器上垂直行走的另一实施例；
图8示出了图7中的信号的功率谱密度；
图9示出了图7中的时间时域数据在2秒到3秒之间的一部分；
图10示出了系统噪声的实施例；
图11示出了图10中的信号的功率谱密度；
图12是能用于生成埋入事件检测统计的过程的示意图；
图13是用于生成埋入事件检测统计的过程的流程图；
图14示出了受关注的频率的实施例；
图15示出了在隐蔽埋入的传感器上垂直行走的位置块表示的实施例；
图16是针对在隐蔽埋入的传感器上的垂直行走使用小波分解的实施例；
图17是针对在隐蔽埋入的传感器上的垂直行走使用小波包分解的实施例；
图18是可用于生成埋入事件检测统计的流程图；
图19A和图19B是示出有人走过传感器产生的两个不同频带的信号；
图20A和图20B示出了有人在传感器上奔跑产生的两个不同频带的信号；
图21给出了不同频带的衰减表示；
图22是在低频带产生的事件的实施例；以及
图23是中频带产生事件的实施例。
具体实施方式
本发明的实施方式提供了一种埋入式入侵检测系统的事件检测方法。在本发明的实施方式中，入侵检测系统基于相干光学时域反射仪(COTDR)技术。示例性实施方式使用快速傅里叶变换(FFT)或离散小波变换(DWT)技术将每个位置仓(沿传感器)信号从时域转换到频域，并对频域信号进行处理以自动分离事件和系统噪声。
图1示出了可以采用本发明的实施方式的一个示例性场景。为了保护埋入的管道110，在管线旁边埋入光纤电缆形式的传感器120。在一个实施例中，电缆的埋入深度为100mm至150mm。在另一实施例中，电缆埋入的深度为100mm至300mm。
如图2中所示，控制器210(例如申请人的Aura Ai控制器)用于向光纤中注入120一系列脉冲220，并监测由Raleigh反向散射240引起的返回光信号250。通过对返回信号的处理，能够检测到当干扰作用于光纤120时的反向散射的变化。
在一个实施例中，控制器210将以每秒2000至5000脉冲的速率输出脉冲宽度为100ns并且峰值功率为125mW的光脉冲。当每个脉冲沿光纤传播时，控制器210以每秒2亿次(200MHz)的频率对反向散射进行采样，以得到电缆每0.5米长度的采样信号250。
该信号250被检测并被数字化以构建原始信号。当脉冲穿越整个光纤长度并接收到相应的反向散射信号时，这被称为"拍摄"(shot)。拍摄是单脉冲对光纤传感器的询问。每一拍摄被划分为0.5米的区段或传感器"仓"。图3中示出了原始拍摄的实施例，其中横轴310是以米为单位的距离，纵轴是以ADC为单位的振幅(ADC＝模数转换器)。
图4是示意图，其示出了在通过沿光纤120发送脉冲220并测量返回信号，为拍摄410的每个仓420创建原始信号后，然后用带通滤波器对其进行滤波以为每个仓产生带通滤波信号430，从而去除COTDR信号中存在的高频噪声。滤波后的信号430看起来类似于干涉信号。如图4中所示，滤波信号430的横轴450为时间(ms)，纵轴为电压(V)。
然而，虽然现有的控制器能够检测一些事件，但其他事件可能被噪声掩盖。例如，可能难以区分滋扰振动140(例如来自经过的火车)和入侵振动130。为了能够辨别诸如上述之类的不同的入侵和滋扰事件，本发明的实施方式旨在提供一种能够将入侵信号与系统噪声信号分开的高灵敏度的系统。
图5示出了由诸如图2中所示的COTDR系统从隐蔽的埋入传感器通过传感器上的垂直行走生成的行走入侵数据的实施例。该实施例数据是通过让"入侵者"从距离传感器位置15米的地方接近埋藏的传感器而生成的。在本实施例中，COTDR系统210参数被设定为200ns脉冲宽度，并且采样频率为2500Hz。图5中的数据是来自COTDR数据表示域上时空滤波数据(15至1200Hz之间的带通滤波(BPF))的单个仓。该仓的位置在沿35千米传感器的594米，并且与"入侵者"进行行走的物理位置对应。如能从图5中看到的，随着时间510的推移，振幅520变更。从时域数据500可以清楚地看出，在信号的区段530中，当入侵者朝传感器行走时，入侵者的步伐被埋在系统噪声内，而在区段540中，入侵者的存在从背景噪声中突显出来。本发明的实施例旨在提供一种即使事件存在于噪声内(例如在区段530中)，也能检测到这些事件的稳健方法。在一些实施方式中，这可以使得能够在入侵者离埋入的电缆更远时检测到入侵者的存在。
图6以线性620和对数标度630示出了图5中信号的功率谱密度610。即，它是由COTDR系统生成的以线性620和dB 630标度表示的行走入侵数据(图5中的信号)的频域。
图7示出了由在隐蔽埋入的传感器上垂直行走产生的行走入侵数据的另一实施例。同样，入侵者从距离传感器位置15米的地方接近埋入的传感器。COTDR系统210参数被设定为200ns脉冲宽度，并且采样频率为2500Hz。图7中的数据来自COTDR数据表示域上的时空滤波数据(在15至1200Hz之间的BPF)的单个仓。该仓的位置沿35千米传感器的32999米。同样，根据时域数据700可以清楚地看出，在信号的区段730中，当入侵者朝传感器行走时，入侵者的步伐被埋在系统噪声内，而在区段740中，入侵者的存在从背景噪声中突显出来。图8以线性820和对数830标度示出了图7中信号的功率谱密度810。
图9示出了图7中的时间时域数据在2秒到3秒之间的部分900。该图示出了入侵者步伐910被埋在噪声里面。所提出的方法具有从系统噪声中解析这些步伐的能力。
图10示出了沿35千米隐蔽埋入的传感器的682米的仓位置处的系统噪声1000的实施例。COTDR系统参数被设定为200ns脉冲宽度，并且采样频率为2500Hz。该图中噪声信号的强度(模数转换(ADC)值1010)与距离传感器隔开一定距离的入侵者的步伐产生的强度相当。图11以线性1120和对数1130标度示出图10中的信号的功率谱密度1110。
在本发明的一些实施方式中，使用诸如快速傅里叶变换(FFT)之类频率变换方法或诸如Welch之类的功率谱估计方法将每个位置仓从时域变换到频域，使得每个仓中的拍摄的每个块由该舱的频域表示的单个像素表示。在COTDR的仓-拍摄(时空)空间上，长度为N的有限长度拍摄的离散傅里叶变换(DFT)定义为：
Figure BDA0002818003060000071
图12是示出用于产生埋入事件检测统计的实施方式的光谱分析过程1300的示意图，图13是该过程1300的流程图。参照图12和图13，图示了由COTDR收集在仓-拍摄空间1200中的数据。即，横轴1210对应于仓，纵轴1220对应于唯一的拍摄。在步骤1310处，对于仓-1至仓-N的每个仓，从多次拍摄收集数据以形成块，因此将有多个块，块-1 1230，块-21240，...，块-N 1250。在一个实施例中，每个块由400次拍摄形成。在一个实施方式中，一旦收集了块的数据，就对数据进行处理。在其他实施方式中，在处理数据之前，收集多个块的数据。
每个仓的数据是上文结合图4所述的带通滤波后的数据(例如)。相应地，傅里叶变换将被应用于滤波后的数据块，例如，对于bin-1而言，来自拍摄块1280A的数据，对于bin-N而言，来自拍摄块1281A的数据。
在步骤1320处，将离散傅里叶变换过程1262(DFT)应用于每个块，以从高通滤波器频率f1到fs/2获得每个频率仓的幅度表示。
然后，方法1300在根据块的频域表示计算事件统计的过程1263中采用两个频带。第一频带是包括[f1，f2]之间的频率仓的信号带。第二频带是远离信号频带的系统噪声频带，并且由[f3，fs/2]之间的频率界定。在其他实施方式中，可以计算多个信号带的事件统计。
在步骤1330处，该方法涉及计算(或推断，compute)由[f1，f2]界定的信号功率区域As，并且在步骤1340处，该方法涉及计算由[f3，fs/2]限界定的噪声功率区域An。图14中示出了区域As和An。
在步骤1350处，该方法涉及使用等式2计算每个块的事件检测统计量，以产生如图11中所示的COTDR数据的仓-块表示1270。
Figure BDA0002818003060000081
其中，
Figure BDA0002818003060000082
Ns是信号频带[f1，f2]中的频率仓的数量，Nn是噪声频段[f3，fs/2]中的频率仓的数量。在这方面，可以理解的是，如果每个频带的频率仓的数量相同，则不需要比例系数。
也就是说，仓-块表示1270的横轴1271表示仓-1至仓-N，而纵轴对应于块数。例如，仓-块值1280B是从1280A的拍摄块计算出来的，并且仓-块值1281B是从1281A的拍摄块计算出来的。每个值在颜色空间中表示，在颜色空间中，颜色的相对强度根据功率比值得出。该颜色空间可以是灰度颜色空间或RGB或CYMK颜色空间。
图15示出了使用图13的方法，在电缆距离约11.5千米处，在隐蔽埋入的传感器上垂直行走的位置块表示(也称为"瀑布"图)。表1示出了图15瀑布图中所示的噪声和入侵步伐数据的功率比统计。统计数据是针对单个仓的块测量的。COTDR系统210参数被设定为200ns脉冲宽度，并采样频率为2500Hz。在该实施例中，功率比的标量值映射到不同的灰度值。在其他实施例中，功率比值映射到不同的颜色。也就是说，横轴1510代表以米为单位的位置，纵轴1520代表以秒为单位的时间。入侵者从距离传感器埋入的地方15米处接近埋入的传感器。此时的传感器长度为11.48千米。COTDR系统参数被设定为200ns脉冲宽度，并且采样频率为2500Hz。
表1
信号数据类型 功率比
纯噪声数据1531 1.24
噪声内的步伐1533 3.76
埋入的传感器顶上的步伐1532 36.16
显然，系统噪声产生的功率比接近于统一。这是由于以下事实：发现该系统中的系统噪声是白噪声，其强度与使As等价于An的不同频率大致相似。此外，随着入侵的增加，传感器顶部上的功率比开始增加到最大值。在埋入式应用中，入侵数据在小于200Hz的频带内有显著的频率，以便进行早期检测。对于诸如挖掘之类的粗暴入侵，当挖掘活动接近埋入的传感器电缆时，频带开始接近奈奎斯特(Nyquist)频率。然而，即使在这种情况下，低频带的能量也比高频带高得多。
在一个实施例中，位置块表示可以由用户在显示器上监测，以便确定是否存在入侵者。在另一个实施例中，可以对位置-块-表示的数据进行进一步的处理，以便例如将用户的注意力引导到位置-块表示的特定部分或产生警报。一个实施例是将阈值应用于统计，并将超过阈值的值作为事件应对。在这样的实施例中，同一仓和/或相邻仓中的事件可以被累积，直到有足够的事件表明各个事件很有可能不是来自滋扰事件，此时，事件的集合被看作"触发器"，并且可以提醒用户或触发警报。
在另选的实施方式中，使用小波分解或小波包变换来测量入侵信号和系统噪声之间的功率比。小波分解进行多级一维小波分析。在小波分解中，使用不同截止频率的滤波器来分析不同尺度(频率)的信号。信号通过一系列高通滤波器(细节系数)以分析高频，并且通过一系列低通滤波器(近似系数)以分析低频。图16示出了关于在隐蔽的埋入传感器上垂直行走的小波分解。入侵者从距离传感器埋入位置15米处接近埋入的传感器。入侵信号被分解为4级。小波分解与小波包的区别在于，图16中的小波分解的高通滤波带(细节系数)被分为高通和低通，以形成图17中的小波包。
图18示出了利用小波产生埋入检测统计的方法1800。与图13的方法中一样，在步骤1810处，对于仓-1至仓-N的每个仓，从多次拍摄中收集数据以形成块。
在步骤1820处，该方法涉及在级别J(在这些实施例中，当采样频率为2500Hz时，将其应用于级别4)在每个块上应用小波分解(图16)或小波包分解(图17)，并从低频率的节点(在图16的实施例中，入侵节点(4,0)1610)获得能量1611，从较高频率的系统噪声节点(在该实施例中是(1,1)1620)获得能量1621。(图17的实施例使用来自低频率节点(节点(4,0)1710)的能量1711和具有较高频率的系统噪声节点(在该实施例例中为(0,0)1720)的能量1721。)
参照图16的实施例，小波分解树上的这两个频带(节点)的使用方式与上述FFT方法的频带类似。在该实施例中，第一频带是来自节点(4,0)1610的信号频带。第二频带是远离信号频带的系统噪声频带，并且可以是节点(1,1)1620。
在步骤1830处，该方法涉及计算信号节点的信号功率区域As，并且在步骤1840处，该方法涉及计算噪声节点的噪声功率区域An。在这两种情况下，信号功率区域都是通过帕西瓦尔定理获得的。
Figure BDA0002818003060000101
其中，C是各频带的系数，N是该节点中信号的总采样数。
步骤1850与图13的方法中的步骤1350相同，即使用等式2计算每个块的事件检测统计，以产生COTDR数据的仓-块表示。
表2示出了使用小波分解方法对来自隐蔽埋入系统的噪声和入侵步伐数据的功率比统计的实施例。针对单个仓数据的块测量统计数据。COTDR系统参数被设定为200ns脉冲宽度，并且采样频率为2500Hz。
表2：使用小波分解法对隐蔽埋入系统的噪声和入侵步伐的功率比统计结果
信号数据类型 功率比
纯噪声数据(如图9中的数据) 1.5944
噪声内的步伐(如图8) 4.6032
埋入的传感器顶部上的步伐 67.2359
一些实施方式使用多个不同的频带，其中使用上述实施方式之一分别计算两个信号频带中每个频带的单独事件统计。由于不同的频率在经由地面传播时的衰减不同。这可以指示这些事件距离传感器发生的相对距离。因此，监测多个频带并采用作为使用来自一个频带的事件的触发器的补充或代替而使用来自两个频带的事件的触发器是有利的。
图19A和图19B示出了有人走过传感器产生的两个不同频带的信号。横轴1910是距离，纵轴1920是时间。在这种情形下，在较低频率(2至5Hz)下观察到的信号(如图19A中所示)比在20至40Hz之间的频率下观察到的信号强得多(如图19B中所示)。这就导致步伐更能可观察，并使发生入侵的置信度更大。20至40Hz频带提供了更清晰的时间定位，使得更容易分离出各个步伐。
不同类型的事件也会产生不同的频率和强度，例如奔跑的人产生较强的冲击力，因此更可能出现较高的频率。在地面上缓慢行走的人会踩得比较轻柔，很可能只会产生较低的频率。
图20A和图20B示出了有人在传感器上奔跑时采集的数据的位置块表示。横轴2010是距离，纵轴2020是时间。在这种情形下，20至40Hz的频带(图20B)更有利于识别更多的步伐。由于较低的频带(图20A)只能拾取传感器几米范围内的脚步声，因此怀疑该频带通过土壤的衰减较大。衰减的差异也可能归因于产生这些频率的声波类型。
图21给出了从1米高处落下一个8公斤重物所产生的不同频带的衰减的表示。该图示出了2至5Hz的第一频带2010、20至40Hz的第二频带2020以及70至150Hz的第三频带2030的衰减。该图示出了20至40Hz频带2020的衰减是显示的3个频带中最低的。
如何使用多个频带来改进入侵检测的实施例是标记来自对更接近事件(例如图22中所示的实施例)更敏感的频带的潜在警报。在这种情况下，2至5Hz频带被用于查找在区域2210中概括的触发事件。
图23示出了由中间频带(20至40Hz)产生的事件(例如2310)，在这个频带中，离传感器更远的地方的更多的事件可能被检测到。这个频率可以更好地用于检测与接近埋入传感器的潜在入侵者相关的事件。较高的频率也能给出更清晰的事件，能够更容易地分离出脚步声。
虽然已经针对图示描述了本发明，但将理解，本领域的技术人员可以在不脱离本发明精神的情况下进行许多变型和变化。根据以上描述和图示的任何变更和派生都包括在权利要求书所限定的本发明的范围内。
本发明的所附权利要求书和前面的描述中，除非由于明确的语言或必要的暗示，上下文另有要求，否则，"包括"一词或诸如"包括"之类的变体是在范围广泛的意义上使用的，即指明所述特征的存在，但不排除在本发明的各种实施方式中存在或增加进一步的特征。

Claims (17)

1.一种用于入侵检测的计算机实施的事件统计生成方法，该方法包括：
针对多个传感器仓中的每一者，将来自相干光学时域反射仪的多个返回信号处理成时域信号，所述多个返回信号对应于在一个时间段内注入到光学传感器光纤中的多个刺激脉冲；以及
对于每个传感器仓：
将各个时域信号变换成相应的频域信号；
从各个频域信号计算预期包含系统噪声的第一频带的第一信号功率区域；
从各个频域信号计算预期包含与至少第一事件相关的任何能量的第二频带的第二信号功率区域；以及
通过将所述第二信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第二信号功率区域与所述第一信号功率区域之比成比例的事件统计。
2.根据权利要求1所述的方法，该方法包括针对多个时间段中的每一个时间段处理多个返回信号，以针对每个传感器仓获得多个事件统计。
3.根据权利要求1或权利要求2所述的方法，该方法包括将每个事件统计转换成指示所述事件统计相对大小的视觉表示，并在显示器上输出各个视觉表示的矩阵。
4.根据权利要求3所述的方法，其中每个视觉表示是灰度值。
5.根据权利要求3所述的方法，其中每个视觉表示是色阶值。
6.根据权利要求1至5中任一项所述的方法，其中生成事件统计包括应用从所述第一频带和第二频带中的频率仓的数量导出的缩放因子。
7.根据权利要求1至6中任一项所述的方法，包括从所述频域信号计算预期包含与至少第二事件相关的任何能量的第三频带的第三信号功率区域，并通过将所述第三信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第三信号功率区域与所述第一信号功率区域之比成比例的进一步的事件统计。
8.根据权利要求1至7中任一项所述的方法，该方法还包括注入所述多个刺激脉冲并接收所述返回信号。
9.一种计算机实施的入侵检测方法，该方法包括从由权利要求1至8中任一项的方法产生的一个或多个事件统计来检测入侵。
10.一种事件统计生成装置，该装置包括：
处理器；和
存储器，其存储指令，当所述处理器执行所述指令时，所述指令使所述处理器：
针对多个传感器仓中的每一者，将来自相干光学时域反射仪的多个返回信号处理成时域信号，所述多个返回信号对应于在一个时间段内注入到光学传感器光纤中的多个刺激脉冲；以及
对于每个传感器仓：
将各个时域信号变换成相应的频域信号；
从各个频域信号计算预期包含系统噪声的第一频带的第一信号功率区域；
从所述频域信号计算预期包含与至少第一事件相关的任何能量的第二频带的第二信号功率区域；以及
通过将所述第二信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第二信号功率区域与所述第一信号功率区域之比成比例的事件统计。
11.根据权利要求10所述的装置，其中当所述处理器执行所述指令时，所述指令使所述处理器针对多个时间段中的每一个时间段处理多个返回信号，以针对每个传感器仓获得多个事件统计。
12.根据权利要求10或权利要求11所述的装置，其中当所述处理器执行所述指令时，所述指令使所述处理器将每个事件统计转换成指示所述事件统计相对大小的视觉表示，并在显示器上输出各个视觉表示的矩阵。
13.根据权利要求12所述的装置，其中每个视觉表示是灰度值。
14.根据权利要求12所述的装置，其中每个视觉表示是色阶值。
15.根据权利要求10至14中任一项所述的装置，其中所述处理器生成事件统计包括应用从所述第一频带和第二频带中的频率仓的数量导出的缩放因子。
16.根据权利要求10至15中任一项所述的装置，其中当所述处理器执行所述指令时，所述指令使所述处理器从所述频域信号计算预期包含与至少第二事件相关的任何能量的第三频带的第三信号功率区域，并通过将所述第三信号功率区域除以所述第一信号功率区域，至少部分地生成与所述第三信号功率区域与所述第一信号功率区域之比成比例的进一步的事件统计。
17.根据权利要求10至16中任一项所述的装置，该装置进一步包括光纤传感器，并且其中所述处理器构成实施相干光时域反射仪的控制器的一部分。
CN201980037983.7A 2018-04-06 2019-04-05 入侵检测的事件统计生成方法及装置 Pending CN112219361A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
AU2018901143A AU2018901143A0 (en) 2018-04-06 A computer-implemented event statistic generation method for intrusion detection
AU2018901143 2018-04-06
PCT/AU2019/050303 WO2019191815A1 (en) 2018-04-06 2019-04-05 Event statistic generation method and apparatus for intrusion detection

Publications (1)

Publication Number Publication Date
CN112219361A true CN112219361A (zh) 2021-01-12

Family

ID=68099665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980037983.7A Pending CN112219361A (zh) 2018-04-06 2019-04-05 入侵检测的事件统计生成方法及装置

Country Status (6)

Country Link
US (2) US11373493B2 (zh)
EP (1) EP3776915A4 (zh)
CN (1) CN112219361A (zh)
AU (2) AU2019248019B2 (zh)
IL (1) IL277826B2 (zh)
WO (1) WO2019191815A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115394024A (zh) * 2022-08-10 2022-11-25 武汉烽理光电技术有限公司 一种基于光栅阵列的步行监测和预测方法及装置

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111104891B (zh) * 2019-12-13 2022-03-08 天津大学 基于BiLSTM的复合特征光纤传感扰动信号模式识别方法
CN111369769B (zh) * 2020-02-17 2020-12-29 柳州达迪通信技术股份有限公司 一种基于数据分析的光缆监测预警系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070008123A1 (en) * 2005-07-06 2007-01-11 The Penn State Research Foundation A networked multiband waveguide intrusion detection and localization sensor
US20100073163A1 (en) * 2007-04-02 2010-03-25 Future Fibre Technologies Pty Ltd Method and apparatus for monitoring a structure
WO2010090752A1 (en) * 2009-02-09 2010-08-12 Bbn Technologies, Corp. Subsurface intrusion detection system
US20120230629A1 (en) * 2009-11-13 2012-09-13 Optasense Holdings Limited Fibre Optic Distributed Sensing
WO2016033199A1 (en) * 2014-08-28 2016-03-03 Adelos, Inc. Real-time fiber optic interferometry controller
CN105931402A (zh) * 2016-06-27 2016-09-07 上海波汇科技股份有限公司 基于图像识别的光纤周界入侵监测方法
US20170260839A1 (en) * 2016-03-09 2017-09-14 Conocophillips Company Das for well ranging
CN107369282A (zh) * 2017-08-23 2017-11-21 武汉理工光科股份有限公司 一种地埋式cotdr周界入侵报警系统及信号处理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2445364B (en) 2006-12-29 2010-02-17 Schlumberger Holdings Fault-tolerant distributed fiber optic intrusion detection
US20090201146A1 (en) * 2007-09-10 2009-08-13 Wayne Lundeberg Remote activity detection or intrusion monitoring system
US20090080898A1 (en) * 2007-09-24 2009-03-26 Fiber Sensys Llc Method and apparatus for reducing noise in a fiber-optic sensor
US8354929B2 (en) * 2008-09-09 2013-01-15 The United States of America, as respresented by the Secretary of the Navy Intruder detector and classifier
WO2011006210A1 (en) * 2009-07-17 2011-01-20 Future Fibre Technologies Pty Ltd Intrusion detection
GB2539254A (en) 2015-06-12 2016-12-14 Pimon Gmbh Method and apparatus for monitoring pipeline

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070008123A1 (en) * 2005-07-06 2007-01-11 The Penn State Research Foundation A networked multiband waveguide intrusion detection and localization sensor
US20100073163A1 (en) * 2007-04-02 2010-03-25 Future Fibre Technologies Pty Ltd Method and apparatus for monitoring a structure
WO2010090752A1 (en) * 2009-02-09 2010-08-12 Bbn Technologies, Corp. Subsurface intrusion detection system
US20120230629A1 (en) * 2009-11-13 2012-09-13 Optasense Holdings Limited Fibre Optic Distributed Sensing
WO2016033199A1 (en) * 2014-08-28 2016-03-03 Adelos, Inc. Real-time fiber optic interferometry controller
US20170260839A1 (en) * 2016-03-09 2017-09-14 Conocophillips Company Das for well ranging
CN105931402A (zh) * 2016-06-27 2016-09-07 上海波汇科技股份有限公司 基于图像识别的光纤周界入侵监测方法
CN107369282A (zh) * 2017-08-23 2017-11-21 武汉理工光科股份有限公司 一种地埋式cotdr周界入侵报警系统及信号处理方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115394024A (zh) * 2022-08-10 2022-11-25 武汉烽理光电技术有限公司 一种基于光栅阵列的步行监测和预测方法及装置
CN115394024B (zh) * 2022-08-10 2024-02-23 武汉烽理光电技术有限公司 一种基于光栅阵列的步行监测和预测方法及装置

Also Published As

Publication number Publication date
EP3776915A1 (en) 2021-02-17
EP3776915A4 (en) 2021-06-02
IL277826A (en) 2020-11-30
US11373493B2 (en) 2022-06-28
IL277826B2 (en) 2023-10-01
US11710390B2 (en) 2023-07-25
AU2019248019B2 (en) 2023-09-07
US20210027589A1 (en) 2021-01-28
IL277826B1 (en) 2023-06-01
WO2019191815A1 (en) 2019-10-10
AU2019248019A1 (en) 2020-11-19
AU2023278003A1 (en) 2023-12-21
US20220309888A1 (en) 2022-09-29

Similar Documents

Publication Publication Date Title
US11710390B2 (en) Event statistic generation method and apparatus for intrusion detection
CN105931402B (zh) 基于图像识别的光纤周界入侵监测方法
KR101688246B1 (ko) 지하 침투 탐지 시스템
Wang et al. Event identification based on random forest classifier for Φ-OTDR fiber-optic distributed disturbance sensor
WO2011006210A1 (en) Intrusion detection
Coviello et al. Detecting torrential processes from a distance with a seismic monitoring network
EP3332272A1 (en) Method and apparatus for monitoring pipeline using an optical fiber sensor system
Ghosh et al. Time–frequency analysis based robust vehicle detection using seismic sensor
CN103994817A (zh) 一种基于长距离光纤多发事件的振动振源的识别方法
CN114613116B (zh) 防外破预警方法、装置、设备及存储介质
Zhu et al. Low frequency acoustic signals associated with rock falls, thunderstorms, and wind turbulences in field environment
CN103994815A (zh) 一种识别光纤振源的方法
JP5507903B2 (ja) 震度推定方法及び装置
CN110657879B (zh) 一种基于fft的分布式光纤振动传感定位方法及装置
Mahmoud et al. Elimination of rain-induced nuisance alarms in distributed fiber optic perimeter intrusion detection systems
CN103994816A (zh) 一种基于光纤多发事件的识别方法
Wu et al. An effective signal separation and extraction method using multi-scale wavelet decomposition for phase-sensitive OTDR system
CA3010122C (en) Collecting apparatus, system and method for gravel transport pressure and transport audio
Arattano et al. Methods of data processing for debris flow seismic warning
Okal et al. Quantification of hydrophone records of the 2004 Sumatra tsunami
Hübl et al. Infrasound measurements of debris flow
RU2768227C1 (ru) Способ оперативно-технической охраны объектов и границ
RU2674809C9 (ru) Способ оперативно-технической охраны рубежей объектов и границ
WO2008038289A2 (en) A system and a method for detecting and classifying damage in a pipeline
Qu et al. Method for identifying mechanical vibration source based on detected signals by optical fiber

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination