CN110823356B - 基于梅尔频谱的分布式光纤入侵检测方法 - Google Patents

基于梅尔频谱的分布式光纤入侵检测方法 Download PDF

Info

Publication number
CN110823356B
CN110823356B CN201910954339.8A CN201910954339A CN110823356B CN 110823356 B CN110823356 B CN 110823356B CN 201910954339 A CN201910954339 A CN 201910954339A CN 110823356 B CN110823356 B CN 110823356B
Authority
CN
China
Prior art keywords
mel
frequency
intrusion
mel frequency
optical fiber
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.)
Active
Application number
CN201910954339.8A
Other languages
English (en)
Other versions
CN110823356A (zh
Inventor
邢陆雁
张妮娜
沙文广
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.)
Weihai Beiyang Photoelectric Information Technology Co ltd
Original Assignee
Weihai Beiyang Photoelectric Information Technology Co 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
Application filed by Weihai Beiyang Photoelectric Information Technology Co ltd filed Critical Weihai Beiyang Photoelectric Information Technology Co ltd
Priority to CN201910954339.8A priority Critical patent/CN110823356B/zh
Publication of CN110823356A publication Critical patent/CN110823356A/zh
Application granted granted Critical
Publication of CN110823356B publication Critical patent/CN110823356B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/004Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B13/00Burglar, theft or intruder alarms
    • G08B13/18Actuation by interference with heat, light, or radiation of shorter wavelength; Actuation by intruding sources of heat, light, or radiation of shorter wavelength
    • G08B13/181Actuation 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/183Actuation 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/186Actuation 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
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/18Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/21Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being power information
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/24Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being the cepstrum
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及输油输气管线中分布式光纤振动监测安防技术领域,具体地说是一种基于梅尔频谱的分布式光纤入侵检测方法,其特征在于包含以下步骤:读取光纤振动传感器采集的原始信号,定位出疑似入侵事件发生的位置点;提取疑似入侵事件发生位置点的信号,生成多个梅尔频谱图;多个梅尔频谱图通过频率映射形成一维特征向量;将一维特征向量输入到深度学习网络中,构建训练模型,并将模型用于光纤入侵事件的实时检测,本发明对低频区域信号响应灵敏,能够将不同事件对信号的低频响应差异化表征出来,能够很好的区分不同入侵事件和干扰事件,提高系统的识别准确率。

Description

基于梅尔频谱的分布式光纤入侵检测方法
技术领域:
本发明涉及输油输气管线中分布式光纤振动监测安防技术领域,具体地说是一种基于梅尔频谱的分布式光纤入侵检测方法。
背景技术
分布式光纤振动入侵检测技术在石油石化输油输气管线的安全监测领域上受到越来越多的关注,相比传统的监测方法,分布式光纤入侵检测技术可以实现长距离安全监测,并且对破坏入侵事件的响应速度快、定位准确,监测覆盖面广,这些优点是传统人工巡线和管道压力检测等手段所无法实现的,所以分布式光纤入侵检测技术的发展变得重要和迫切。
分布式光纤入侵检测方法是分布式入侵检测技术应用发展的一个重要因素。目前大多数方法都是将时域信号变换到频域中进行分析和识别,采用最多的时频变换方法是傅里叶变换和小波变换,这些时频分析方法可以将信号所覆盖的所有频率范围都提取出来,从整个频率范围分析信号特征,但是在实际中采集的入侵信号,由于土质对振动的吸收,信号的频率范围主要集中在低频区域,所以信号分析的重点是低频范围的信号特征。现阶段需要在保证计算速度的前提下,进一步提高系统对不同事件信号的响应灵敏度,以区分不同的入侵事件和干扰事件。
发明内容:
本发明针对现有技术中存在的缺点和不足,提出了一种对低频区域信号响应灵敏,能够将不同事件对信号的低频响应差异化表征出来,能够很好的区分不同入侵事件和干扰事件,提高系统的识别准确率的基于梅尔频谱的分布式光纤入侵检测方法。
本发明可以通过以下措施达到:
一种基于梅尔频谱的分布式光纤入侵检测方法,其特征在于包含以下步骤:
步骤一:读取光纤振动传感器采集的原始信号,定位出疑似入侵事件发生的位置点;
步骤二:提取疑似入侵事件发生位置点的信号,生成多个梅尔频谱图;
步骤三:多个梅尔频谱图通过频率映射形成一维特征向量;
步骤四:将一维特征向量输入到深度学习网络中,构建训练模型,并将模型用于光纤入侵事件的实时检测。
本发明所述步骤一中,光纤振动传感器采集的原始信号是包含时域信息和空间信息的二维数据矩阵,可表示为Dt×l=(dij)t×l,(i=1,2...,t;j=1,2...,l),其中t表示时间维度,l表示空间维度,定位出疑似入侵事件发生的位置点是将采集到的原始信号在时间方向上提取出能量信息,然后根据能量信息在空间位置上的相互关系确定疑似入侵事件发生的位置点,具体的实现步骤为:
步骤1-1:原始信号在时间方向采用滤波的方式将直流分量和高频噪声滤除;步骤1-2:将滤波后的信号沿时间方向求取能量值:
Figure BDA0002226776760000011
步骤1-3:遍历El中的点,对每个点按如下规则进行筛选:
Figure BDA0002226776760000021
将满足如上条件的E’和E’所在的坐标位置Loc进行记录,记录的这些位置坐标Loc即为疑似入侵事件发生的位置点。
本发明步骤二中,所述提取疑似入侵事件发生位置点的信号,是指根据步骤一中得到的位置坐标Loc将原始信号中对应位置的信号提取出来,所述生成梅尔频谱是指将提取的每个疑似入侵位置点的信号都生成对应的梅尔频谱,生成步骤为:
步骤2-1:信号分帧,即将每个位置点的一维信号进行分帧处理;
步骤2-2:计算每帧的傅里叶变换及其功率谱;
步骤2-3:使用梅尔滤波器组将每帧的傅里叶频谱变换到梅尔频谱域中,变换公式为:
M(f)=1125ln(1+f/700)
步骤2-4:每一帧信号都重复上述步骤3操作,所有帧的信号将形成多个梅尔频谱图,按照帧顺序将多个梅尔频谱图进行组合。
本发明步骤三中,所述将多个梅尔频谱图通过频率映射形成一维特征向量是指对每个梅尔频谱图都通过映射形成一个特征值,多个梅尔频谱图通过映射后形成的多个特征值组合在一起形成一维的特征向量,具体实现步骤为:
步骤3-1:构建频谱映射关系:
Figure BDA0002226776760000022
其中,value为计算的特征值,E为梅尔频率为melf时的能量值,μ为权重,melf0、melf1和melf2为频带截止频率;
步骤3-2:每个梅尔频谱都根据步骤1中的映射关系,将频谱图中的频率分布映射到一个特征值上;
步骤3-3:根据分帧顺序,将多个特征值进行组合形成该疑似入侵位置点的特征向量。
本发明上述步骤四中,所述构建深度学习训练模型的步骤为:
步骤4-1:采集大量不同入侵行为和扰动行为的样本信号;
步骤4-2:根据上述步骤提取所有样本信号的特征向量;
步骤4-3:将所有样本的特征向量分成训练样本和测试样本,使用训练样本对深度学习网络进行训练,使用测试样本对算法性能进行评估,并不断优化参数和算法结构,直至算法性能达到最优;通过上述步骤训练的模型可以应用于系统中,实时检测破坏管道的入侵行为并发出报警。
本发明与现有技术相比,采用了预判策略,在对入侵事件检测前先对光纤各个位置点进行筛选,保留疑似入侵事件发生的位置点并对这些位置点的信号进行进一步识别,这种策略可以提高计算速度,保证系统的实时计算性能。
附图说明
附图1为本发明实施例中能量曲线图。
附图2为本发明实施例中疑似入侵事件发生位置点的信号波形图。
附图3为本发明实施例中梅尔频谱图。
附图4为本发明实施例中特征向量曲线
具体实施方式
为了使本发明的技术方案更加清楚,以下我们结合附图进一步对本发明进行描述。
本发明提出了一种基于梅尔频谱的分布式光纤入侵检测方法,本方法在对疑似入侵事件发生位置预筛选的基础上,基于梅尔频谱提取特征,并通过深度学习网络进行最终的入侵事件判别,本方法的计算步骤如下:1)读取光纤振动传感器采集的原始信号,定位出疑似入侵事件发生的位置点;2)将疑似入侵事件发生位置点的信号生成梅尔频谱图;3)对梅尔频谱图进行映射,形成一维特征向量;4)将一维特征向量输入到深度学习网络中,构建训练模型,并将模型用于光纤入侵事件的实时检测。
其中,本发明所述定位出疑似入侵事件发生的位置点是指采用信号能量生成能量曲线,根据能量曲线定位出疑似入侵事件发生的位置点,具体为:原始信号在时间方向采用滤波的方式将直流分量和高频噪声滤除;将滤波后的信号沿时间方向求取能量值,生成能量曲线;遍历能量曲线中点,对每个点按照一定规则进行筛选,将满足规则的点进行保留。
所述将疑似入侵事件发生位置点的信号生成梅尔频谱图是指对筛选出来的信号进行分帧,对每帧信号提取梅尔频谱,最终形成多个梅尔频谱,具体实现步骤为:1)对筛选出来的每个位置点的一维信号进行分帧处理;2)计算每帧的傅里叶变换及其功率谱;3)使用梅尔滤波器组将每帧的傅里叶频谱变换到梅尔频谱域中,生成梅尔频谱图。
所述对梅尔频谱图进行映射,形成一维特征向量是指:将每个二维频谱图都通过映射关系映射为一个特征值,多个频谱图生成的特征值组合起来形成一维特征向量,其中一维特征向量中可以包含一个时间段的特征值。
所述深度学习网络模型训练和识别是指对一维特征向量组成的训练样本采用自定义深度学习网络进行训练,得到最优的深度学习网络模型。在实际应用中,采用最优深度学习网络模型实时的对各种入侵事件和干扰事件进行识别,准确区分出入侵事件,过滤掉干扰事件。
实施例1:
首先,读取原始信号:原始信号的大小为2000*5000,其中2000是每个位置点的采样频率,5000是空间位置点的个数。
其次,定位出疑似入侵事件发生的位置点:原始信号采用带通滤波器沿垂直方向进行滤波,得到去噪声和去直流的信号,然后对每个位置点的信号求能量得到如图1的能量曲线,曲线中的每个点表示在某一位置的信号能量值,最后根据该能量曲线,找到多个疑似入侵事件发生的位置点,如图中红色圈所标注。
然后,提取疑似入侵事件发生位置点的信号,生成多个梅尔频谱图:图1中得出位置1600是一个疑似入侵事件发生的位置点,将该位置点的原始信号提取出来,如图2所示,对该信号进行分帧,共分为5帧,每帧都生成梅尔频谱,梅尔滤波器组个数为40,所有梅尔频谱生成的梅尔声谱图如图3所示。
多个梅尔频谱图通过频谱映射形成一维特征向量:每个梅尔频谱图都通过频率映射关系生成一个特征值,多个梅尔频谱生成多个特征值,多个特征值组合形成特征向量,其中本实例中共包含了20秒的共100个特征的组合,图4是生成特征向量曲线。
构建深度学习网络模型:
首先,建立数据样本库,样本库中包括挖掘机挖掘、人工挖掘、车辆经过三种行为共5000个样本,其中3000条用于模型训练,2000条用于测试;然后搭建深度学习算法网络,并通过多次迭代优化得到最终的网络结构,最后将训练得到的模型用于系统中,对实时出现的入侵事件进行报警。
本发明提出的基于梅尔频谱的分布式光纤入侵检测方法根据实际入侵事件发生的频率响应规律,结合梅尔频谱对低频区域响应灵敏的特点,能很好的将入侵事件所引起的低频振动差异化的提取出来,使系统对不同事件的识别准确率提高,并且在入侵事件判断前通过位置筛选,先确定疑似入侵位置点,节省了计算时间,保证了实时计算的性能要求。

Claims (1)

1.一种基于梅尔频谱的分布式光纤入侵检测方法,其特征在于由以下步骤组成:
步骤一:读取光纤振动传感器采集的原始信号,定位出疑似入侵事件发生的位置点;
步骤二:提取疑似入侵事件发生位置点的信号,生成多个梅尔频谱图;
步骤三:多个梅尔频谱图通过频率映射形成一维特征向量;
步骤四:将一维特征向量输入到深度学习网络中,构建训练模型,并将模型用于光纤入侵事件的实时检测;
步骤一中,光纤振动传感器采集的原始信号是包含时域信息和空间信息的二维数据矩阵,可表示为Dt×l=(dij)t×l,i=1,2...,t;j=1,2...,l,其中t表示时间维度,l表示空间维度,定位出疑似入侵事件发生的位置点是将采集到的原始信号在时间方向上提取出能量信息,然后根据能量信息在空间位置上的相互关系确定疑似入侵事件发生的位置点,具体的实现步骤为:
步骤1-1:原始信号在时间方向采用滤波的方式将直流分量和高频噪声滤除;
步骤1-2:将滤波后的信号沿时间方向求取能量值:
Figure FDA0003526010130000011
步骤1-3:遍历El中的点,对每个点按如下规则进行筛选:
Figure FDA0003526010130000012
将满足如上条件的E’和E’所在的坐标位置Loc进行记录,记录的这些位置坐标Loc即为疑似入侵事件发生的位置点;
步骤二中,所述提取疑似入侵事件发生位置点的信号,是指根据步骤一中得到的位置坐标Loc将原始信号中对应位置的信号提取出来,所述生成梅尔频谱是指将提取的每个疑似入侵位置点的信号都生成对应的梅尔频谱,生成步骤为:
步骤2-1:信号分帧,即将每个位置点的一维信号进行分帧处理;
步骤2-2:计算每帧的傅里叶变换及其功率谱;
步骤2-3:使用梅尔滤波器组将每帧的傅里叶频谱变换到梅尔频谱域中,变换公式为:
M(f)=1125ln(1+f/700)
步骤2-4:每一帧信号都重复上述步骤三 操作,所有帧的信号将形成多个梅尔频谱图,按照帧顺序将多个梅尔频谱图进行组合;
步骤三 中,所述将多个梅尔频谱图通过频率映射形成一维特征向量是指对每个梅尔频谱图都通过映射形成一个特征值,多个梅尔频谱图通过映射后形成的多个特征值组合在一起形成一维的特征向量,具体实现步骤为:
步骤3-1:构建频谱映射关系:
Figure FDA0003526010130000021
其中,value为计算的特征值,E为梅尔频率为melf时的能量值,μ为权重,melf0、melf1和melf2为频带截止频率;
步骤3-2:每个梅尔频谱都根据步骤3-1中的映射关系,将频谱图中的频率分布映射到一个特征值上;
步骤3-3:根据分帧顺序,将多个特征值进行组合形成该疑似入侵位置点的特征向量;
步骤四中,所述构建深度学习训练模型的步骤为:
步骤4-1:采集大量不同入侵行为和扰动行为的样本信号;
步骤4-2:根据上述步骤提取所有样本信号的特征向量;
步骤4-3:将所有样本的特征向量分成训练样本和测试样本,使用训练样本对深度学习网络进行训练,使用测试样本对算法性能进行评估,并不断优化参数和算法结构,直至算法性能达到最优;通过上述步骤训练的模型可应用于系统中,实时检测破坏管道的入侵行为并发出报警。
CN201910954339.8A 2019-10-09 2019-10-09 基于梅尔频谱的分布式光纤入侵检测方法 Active CN110823356B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910954339.8A CN110823356B (zh) 2019-10-09 2019-10-09 基于梅尔频谱的分布式光纤入侵检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910954339.8A CN110823356B (zh) 2019-10-09 2019-10-09 基于梅尔频谱的分布式光纤入侵检测方法

Publications (2)

Publication Number Publication Date
CN110823356A CN110823356A (zh) 2020-02-21
CN110823356B true CN110823356B (zh) 2022-04-19

Family

ID=69548813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910954339.8A Active CN110823356B (zh) 2019-10-09 2019-10-09 基于梅尔频谱的分布式光纤入侵检测方法

Country Status (1)

Country Link
CN (1) CN110823356B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111238552B (zh) * 2020-02-27 2021-06-22 南昌航空大学 基于深度学习的分布式光纤传感系统扰动定位方法
CN112504428A (zh) * 2020-10-19 2021-03-16 威海北洋光电信息技术股份公司 低功耗嵌入式高速分布式光纤振动传感系统及其应用
CN114354744B (zh) * 2021-12-03 2023-08-29 国网河南省电力公司商丘供电公司 一种电力电缆防外破监测方法、系统及装置
GB2614421A (en) * 2021-12-30 2023-07-05 Aselsan Elektronik Sanayi Ve Ticaret Anonim Sirketi Frequency response estimation method to compensate for channel differences in distributed acoustic sensing systems
CN114495393B (zh) * 2022-02-28 2023-06-16 中国人民公安大学 一种基于物联网的地下异常振动安全监测系统
CN114998837A (zh) * 2022-06-14 2022-09-02 中电海康集团有限公司 入侵事件的实时识别方法、装置、系统及介质
CN116186642B (zh) * 2023-04-27 2023-09-08 山东汇英光电科技有限公司 一种基于多维特征融合的分布式光纤传感事件预警方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9476760B2 (en) * 2010-06-25 2016-10-25 Schlumberger Technology Corporation Precision measurements in a fiber optic distributed sensor system
CN103514877A (zh) * 2013-10-12 2014-01-15 新疆美特智能安全工程股份有限公司 振动信号特征参数提取方法
CN105095624B (zh) * 2014-05-15 2017-08-01 中国电子科技集团公司第三十四研究所 一种光纤传感振动信号的识别方法
CN106482912B (zh) * 2016-11-29 2018-12-28 中国工程物理研究院激光聚变研究中心 一种真空设备泄漏检测及定位方法

Also Published As

Publication number Publication date
CN110823356A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN110823356B (zh) 基于梅尔频谱的分布式光纤入侵检测方法
CN109212597B (zh) 基于深度信念网络和扫描叠加的微震源自动定位方法
Yu et al. A new statistical modeling and detection method for rolling element bearing faults based on alpha–stable distribution
CN102004264A (zh) 一种地震采集资料质量定量分析与评价方法
CN103064111A (zh) 一种基于形态滤波的微地震事件识别方法
CN105527650A (zh) 一种工程尺度下微震信号及p波初至自动识别算法
CN111582160A (zh) 基于分布式光纤传感系统的车辆轨迹提取方法及装置
CN108648763B (zh) 基于声学信道的个人计算机使用行为监测方法及系统
CN111399056B (zh) 一种基于划分方位角滤波预测裂缝强度的方法
CN105626502A (zh) 基于小波包和拉普拉斯特征映射的柱塞泵健康评估方法
CN101221249A (zh) 地震资料品质分析技术
CN106019377B (zh) 一种基于时空域降频模型的二维地震勘探噪声去除方法
Li et al. Distributed acoustic sensing vertical seismic profile data denoising based on multistage denoising network
CN108344795A (zh) 油气管道缺陷识别方法、装置及电子设备
Zhou et al. Unsupervised machine learning for waveform extraction in microseismic denoising
CN111694053A (zh) 初至拾取方法及装置
Lara et al. A real-time microearthquakes-detector based on voice activity detection and endpoint detection: An approach to Cotopaxi Volcano
CN112801033A (zh) 基于AlexNet网络的长输油气管线沿线施工扰动及泄漏识别方法
CN104570118B (zh) 一种基于双因素的自动识别与去除工业干扰的方法
CN116186642A (zh) 一种基于多维特征融合的分布式光纤传感事件预警方法
CN110727025A (zh) 一种隐蔽断层的识别方法
CN114626420A (zh) 一种基于CWT-Xception-RF的φ-OTDR振动事件分类方法
Liu et al. Application of convolutional neural network to GIS and physics
CN112464777B (zh) 一种光纤振动源垂直距离智能估算方法
CN112099080B (zh) 一种基于局部叠加检测地面微地震事件的方法

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
GR01 Patent grant
GR01 Patent grant