CN110823356B - 基于梅尔频谱的分布式光纤入侵检测方法 - Google Patents
基于梅尔频谱的分布式光纤入侵检测方法 Download PDFInfo
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 40
- 239000013307 optical fiber Substances 0.000 title claims abstract description 25
- 238000001514 detection method Methods 0.000 title claims abstract description 15
- 238000013507 mapping Methods 0.000 claims abstract description 21
- 239000013598 vector Substances 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 8
- 230000009545 invasion Effects 0.000 claims abstract description 4
- 238000013135 deep learning Methods 0.000 claims description 14
- 238000012549 training Methods 0.000 claims description 14
- 230000006399 behavior Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000009432 framing Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000011897 real-time detection Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000004044 response Effects 0.000 abstract description 7
- 238000012544 monitoring process Methods 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 230000004069 differentiation Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Classifications
-
- 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
-
- 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
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/18—Speech 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/21—Speech 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
- G10L25/24—Speech 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
- G10L25/51—Speech 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:将滤波后的信号沿时间方向求取能量值:
步骤1-3:遍历El中的点,对每个点按如下规则进行筛选:
将满足如上条件的E’和E’所在的坐标位置Loc进行记录,记录的这些位置坐标Loc即为疑似入侵事件发生的位置点。
本发明步骤二中,所述提取疑似入侵事件发生位置点的信号,是指根据步骤一中得到的位置坐标Loc将原始信号中对应位置的信号提取出来,所述生成梅尔频谱是指将提取的每个疑似入侵位置点的信号都生成对应的梅尔频谱,生成步骤为:
步骤2-1:信号分帧,即将每个位置点的一维信号进行分帧处理;
步骤2-2:计算每帧的傅里叶变换及其功率谱;
步骤2-3:使用梅尔滤波器组将每帧的傅里叶频谱变换到梅尔频谱域中,变换公式为:
M(f)=1125ln(1+f/700)
步骤2-4:每一帧信号都重复上述步骤3操作,所有帧的信号将形成多个梅尔频谱图,按照帧顺序将多个梅尔频谱图进行组合。
本发明步骤三中,所述将多个梅尔频谱图通过频率映射形成一维特征向量是指对每个梅尔频谱图都通过映射形成一个特征值,多个梅尔频谱图通过映射后形成的多个特征值组合在一起形成一维的特征向量,具体实现步骤为:
步骤3-1:构建频谱映射关系:
其中,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:将滤波后的信号沿时间方向求取能量值:
步骤1-3:遍历El中的点,对每个点按如下规则进行筛选:
将满足如上条件的E’和E’所在的坐标位置Loc进行记录,记录的这些位置坐标Loc即为疑似入侵事件发生的位置点;
步骤二中,所述提取疑似入侵事件发生位置点的信号,是指根据步骤一中得到的位置坐标Loc将原始信号中对应位置的信号提取出来,所述生成梅尔频谱是指将提取的每个疑似入侵位置点的信号都生成对应的梅尔频谱,生成步骤为:
步骤2-1:信号分帧,即将每个位置点的一维信号进行分帧处理;
步骤2-2:计算每帧的傅里叶变换及其功率谱;
步骤2-3:使用梅尔滤波器组将每帧的傅里叶频谱变换到梅尔频谱域中,变换公式为:
M(f)=1125ln(1+f/700)
步骤2-4:每一帧信号都重复上述步骤三 操作,所有帧的信号将形成多个梅尔频谱图,按照帧顺序将多个梅尔频谱图进行组合;
步骤三 中,所述将多个梅尔频谱图通过频率映射形成一维特征向量是指对每个梅尔频谱图都通过映射形成一个特征值,多个梅尔频谱图通过映射后形成的多个特征值组合在一起形成一维的特征向量,具体实现步骤为:
步骤3-1:构建频谱映射关系:
其中,value为计算的特征值,E为梅尔频率为melf时的能量值,μ为权重,melf0、melf1和melf2为频带截止频率;
步骤3-2:每个梅尔频谱都根据步骤3-1中的映射关系,将频谱图中的频率分布映射到一个特征值上;
步骤3-3:根据分帧顺序,将多个特征值进行组合形成该疑似入侵位置点的特征向量;
步骤四中,所述构建深度学习训练模型的步骤为:
步骤4-1:采集大量不同入侵行为和扰动行为的样本信号;
步骤4-2:根据上述步骤提取所有样本信号的特征向量;
步骤4-3:将所有样本的特征向量分成训练样本和测试样本,使用训练样本对深度学习网络进行训练,使用测试样本对算法性能进行评估,并不断优化参数和算法结构,直至算法性能达到最优;通过上述步骤训练的模型可应用于***中,实时检测破坏管道的入侵行为并发出报警。
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)
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 | 中国人民公安大学 | 一种基于物联网的地下异常振动安全监测*** |
CN114998837B (zh) * | 2022-06-14 | 2024-07-09 | 中电海康集团有限公司 | 入侵事件的实时识别方法、装置、***及介质 |
CN116186642B (zh) * | 2023-04-27 | 2023-09-08 | 山东汇英光电科技有限公司 | 一种基于多维特征融合的分布式光纤传感事件预警方法 |
Family Cites Families (4)
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 | 中国工程物理研究院激光聚变研究中心 | 一种真空设备泄漏检测及定位方法 |
-
2019
- 2019-10-09 CN CN201910954339.8A patent/CN110823356B/zh active Active
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) | 一种地震采集资料质量定量分析与评价方法 | |
CN105527650A (zh) | 一种工程尺度下微震信号及p波初至自动识别算法 | |
CN111582160A (zh) | 基于分布式光纤传感***的车辆轨迹提取方法及装置 | |
CN104266894A (zh) | 一种基于相关性分析的矿山微震信号初至波时刻提取方法 | |
CN108648763B (zh) | 基于声学信道的个人计算机使用行为监测方法及*** | |
CN111399056A (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
CN105626502A (zh) | 基于小波包和拉普拉斯特征映射的柱塞泵健康评估方法 | |
Akhouayri et al. | Automatic detection and picking of P-wave arrival in locally stationary noise using cross-correlation | |
CN101221249A (zh) | 地震资料品质分析技术 | |
Li et al. | Distributed acoustic sensing vertical seismic profile data denoising based on multistage denoising network | |
CN106019377B (zh) | 一种基于时空域降频模型的二维地震勘探噪声去除方法 | |
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 | |
CN116186642B (zh) | 一种基于多维特征融合的分布式光纤传感事件预警方法 | |
CN104570118B (zh) | 一种基于双因素的自动识别与去除工业干扰的方法 | |
CN114626420A (zh) | 一种基于CWT-Xception-RF的φ-OTDR振动事件分类方法 | |
CN112464777B (zh) | 一种光纤振动源垂直距离智能估算方法 | |
CN117908128B (zh) | 一种波数域优势频带识别岩性体断溶体边界的检测方法 | |
CN112099080B (zh) | 一种基于局部叠加检测地面微地震事件的方法 | |
Chen et al. | Pre-earthquake prediction method based on electromagnetic wave abnormal signal detection |
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 |