CN105929457A

CN105929457A - 一种地下径流动态信息的频谱识别方法

Info

Publication number: CN105929457A
Application number: CN201610235250.2A
Authority: CN
Inventors: 黄采伦; 陈超洋; 王靖; 吴亮红; 赵延明; 何斌华
Original assignee: Hunan University of Science and Technology
Current assignee: Hunan University of Science and Technology
Priority date: 2016-04-18
Filing date: 2016-04-18
Publication date: 2016-09-07
Anticipated expiration: 2036-04-18
Also published as: CN105929457B

Abstract

本发明公开一种地下径流动态信息的频谱识别方法，通过探测仪器接收天然电场信号并细化得到分析带宽f ₁～f ₂的分辨率为∆f的256条谱线，经最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据、动态信息识别判据后得出有无地下径流动态信息的结论并相应给出主频信息、动态信息的幅值、频率。本发明从研究地下径流动态信息的形成机理出发，运用信号调制原理与边频特性相结合的原则在细化频谱内自动识别地下径流动态信息，并充分考虑了干扰因素的影响，可使地下径流动态信息识别准确率达到90%以上，解决了现有物探方法中难以识别地下径流动态信息的问题。

Description

一种地下径流动态信息的频谱识别方法

技术领域

本发明属于物探成果分析技术，具体地说是一种用于地下径流动态信息分析的频谱识别方法。

背景技术

地下径流又称地下水径流、地下径流水，是重力作用下地下水在自然界陆地水循环过程中的流动。地下径流是岩石介质中由水头高处向水头低处流动的地下水流，它主要通过断裂带等连通的空隙流动。流动带走空隙中的物质，淘空断裂带，影响断裂带上建筑物的稳定、矿山的安全。掌握地下径流通道的类型、特征和快速、有效探测其位置是矿山防治水和地下水勘探过程中的关键问题之一。

地下水资源探测作为一个世界性的科技难题，就目前来说，仍未能得到很好的解决。目前国内外地面物探方法找水首推电法，其次有放射性探测法及测量地层磁分量的甚低频法。采用甚低频仪测地下磁场随空间变化的方法，由于受现代工业发展所带来的强电波干扰，很难在城市及郊区开展探测工作。放射性找水是利用天然放射性元素氟在岩石裂隙中富集造成放射性异常，仪器可测到异常带，但这个异常带可以是充填的、张开的，是否有水流动仍然不知，所以与地质目测断层找水差不多。自二十世纪30年代以来，我国开始用电法找矿，同时也用电法来找水；到目前为止，单就找水方面就有许多方法，如电法就有直流电法、音频电法等，其中直流电法又可分为电测深与电剖面法，分别用于测量纵深方向的地质变化和某一深度沿剖面方向的地质变化。上述方法的共同点是在人工场作用下进行测量，把寻找固体矿产的物探方法应用在找地下径流水上。地面仪器测量值反映的是地质体物性综合值，属于静态信息；这个物理量所显示的是地下的那种固体矿产或地下径流水全凭解释者的主观经验。所以上述方法的找水成功率只有40～50％，其根源在物探曲线的多解性。为解决了以往勘探中主要依据静态信息探测断裂带等地质问题精度较低的不足，发展了天然电场选频、地下磁流体探测等方法，从分析天然电场信号的时域、频域特征的角度出发来解决地下各种地质体探测的问题，特别是地下径流动态信息的发现，为地下径流水探测找到了一种准确、经济、快速的方法。

地下径流是一个动态系统，是指在有关因素影响下，地下水的水位、水量、水化学成分、水温等随时间的变化状况。地下水动态提供含水层或含水系统的系列信息。在验证所作出的水文地质结论或所采取的水文地质措施是否正确时，地下水动态是十分重要的。地下水动态受气候、水文、地质和人类活动等因素的影响，受气候、水文、地质等因素影响的，称“天然因素影响的地下水动态”，受人类活动影响的称“人类活动影响下的地下水动态”。研究地下水动态有助于解决一系列理论和实际问题。分析地下水动态可以帮助查明补给来源，查明含水层之间或含水层与地表水体之间的联系情况。确定供水井的深度时，需要了解最低水位，以保证干旱季节和干旱年份的水量供应。计算地下水资源，必须具备一定年限的地下水动态观测资料。监测人为活动影响下的地下水动态，可以及早发现不利变化（如咸水入侵淡含水层，地下水污染），不失时机地采取措施。地震前地应力的变化会引起地下水位乃至水质异常变化。因此，观测地下水动态可作为预报地震的一种辅助手段。现有物探方法中，虽然发现了地下径流动态信息的存在，但对其识别和判断的方法研究较少，目前主要依靠人工看图识别，存在效率低、受人为因素影响大等不足。能否发明一种将地下径流与地下固体矿产资源区分开来的分析方法，目前国内外科技界一直没有解决，尚无相关的研究成果和产品报道。

发明内容

针对现有物探方法在地下径流动态信息识别方面存在的不足，本发明为解决该问题提供了一种地下径流动态信息的频谱识别方法。

本发明的术方案是：一种地下径流动态信息的频谱识别方法，通过探测仪器接收天然电场信号并细化得到分析带宽f ₁～f ₂的分辨率为∆f的256条谱线，经最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据、动态信息识别判据后得出有无地下径流动态信息的结论并相应给出主频信息、动态信息的幅值、频率。

在本发明中，最大谱线离散判据是在256条谱线中找到最大幅值谱线F_MAX0并计算其及前后各3条谱线一起的平均幅值F_J0，若F_MAX0≥3F_J0则离散，反之不离散。

在本发明中，频谱能量集中判据是计算256条谱线的平均幅值F_J，若F_J0≥10F_J则频谱能量集中，否则频谱能量发散。

在本发明中，极点谱线搜索判据是缓存256条谱线到一指定存储区，在其中按从前往后两两比较幅值低者置0、从后往前两两比较幅值低者置0的原则找出频谱中的所有极值点，结合原频谱计算次大极值点谱线F_MAX1、F_MAX2及其前后各3条谱线一起的平均幅值F_J1、F_J2，若F_MAX1≥3F_J1则F_MAX1离散、F_MAX2≥3F_J2则F_MAX2离散，然后再将幅值小于F_MAX0/30的极值点置0。

在本发明中，动态信息识别判据是若频谱能量集中且3≤极值点个数≤30，则通过边频识别规则一、规则二、规则三判断是否存在关于F_MAX0、F_MAX1及F_MAX2的对称极值谱线来确定地下径流动态信息的存在；否则主频信息为F_MAX0、动态信息为0且退出识别。

所述的对称极值谱线是指两条极值谱线中间存在另一条极值谱线使两条极值谱线与其点距之差在±2条谱线范围内，称该对称极值谱线为边频，其频率为一阶边频频率f _b、幅值为两条谱线的幅值均值。

所述的边频识别规则一是F_MAX0离散且极值点分布于其左、右两边，若存在关于F_MAX0对称的边频极值谱线，则主频信息为F_MAX0，动态信息为一阶边频的信息；若不存在关于F_MAX0对称的边频极值谱线、但存在关于F_MAX1或F_MAX2对称的边频极值谱线且相应的F_MAX1或F_MAX2离散，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。

所述的边频识别规则二是F_MAX0离散但主要极值点分布在其一边，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。

所述的边频识别规则三是F_MAX0不离散，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；若F_MAX1和F_MAX2均不离散但存在包括F_MAX0、F_MAX1、F_MAX2中至少2个在内的4个以上的点距误差在±2条谱线内的极值谱线，则主频信息的频率为F_MAX0谱线频率、幅值为F_MAX0与F_MAX1的均值，动态信息的频率为一阶边频频率、幅值为F_MAX1与F_MAX2的均值；否则主频信息为F_MAX0、动态信息为0且退出识别。

本发明的有益效果是：从研究地下径流动态信息的形成机理出发，运用信号调制原理与边频特性相结合的原则在细化频谱内自动识别地下径流动态信息，并充分考虑了干扰因素的影响，可使地下径流动态信息识别准确率达到90%以上，解决了现有物探方法中难以识别地下径流动态信息的问题。

附图说明

图1是本发明的动态信息识别流程图；

图2是本发明的动态信息形成原理示意图；

图3是本发明的实施例一细化谱图；

图4是本发明的实施例二细化谱图；

图5是本发明的实施例三细化谱图；

图6是本发明的实施例四细化谱图；

图7是本发明的实施例五细化谱图；

图8是本发明的实施例六细化谱图；

图9是本发明的实施例七细化谱图；

图10是本发明的实施例八细化谱图。

具体实施方式

下面结合附图和实施例对本发明的方案作进一步详细说明。

参见附图1、附图2，一种地下径流动态信息的频谱识别方法，探测仪器接收天然电场信号并细化得到分析带宽f ₁～f ₂的分辨率为∆f的256条谱线，经最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据、动态信息识别判据后得出有无地下径流动态信息的结论并相应给出主频信息、动态信息的幅值、频率。由附图2可知：岩溶、裂隙中脉状流的地下水，按水文学规律流动在地底下会切割地磁场磁力线产生感应电动势f _b，但该感应电动势较微弱，无法穿透地层传送到地面；但它会对经过该处的天然电磁波f _c产生持续不断的干扰，最终附加在该电磁波信号上传输到地表，把这种附加有地下水流动信息的信号定义为地下径流动态信息。探测仪器通过在地表接收该天然电磁波的分布电场信号，并据此分析地下径流的赋存及状态。天然电磁波在穿透地层的过程中与不同结构的地质体发生折射、反射等耦合，以及与流动的地下径流的相互作用，都会产生天然的电磁耦合效应；其穿透到地面的电磁信号，是包含上述不同作用的信号，且不同的地质构造表现出的特征都不完全相同。目前物探找水方法主要是电法，它是通过探测地下低阻区域来区分地下水与矿物质；该方法从逻辑上看较合理，实际上这种方法由于其结果的多解性，成功率不高；因为低值异常可以由各种原因造成，所有金属矿物、煤、页岩、断层泥、软硬岩互层等都会反映低值异常；怎样才能使物探结果由多解上升到单一解、确定解，这是本发明方法要解决的问题。只要细心观察一下上述地质体与地下径流的共性和差异就可以找到解决问题的途径。以上金属矿等地质体都是固体，在地下处于静态；而地下径流与它们的共同点是低阻，但地下水是流动的，它不仅传递静态信息，而且还能传递天然动态信息；地下水具有低阻、良导体、流动、流量随时间变化等特性。地下水、天然气、石油，都是地下存在的流体，其中只有地下水是良导体，在地磁场中地下水在岩石裂缝、溶洞中运动就要切割磁力线，按法拉第定律，水流产生电流、环绕水流产生磁场，地下水与地磁场进行反复复杂的电磁感应，结果在地下形成不断发出瞬变电磁波的场源；这样在地磁场中，地下有地下水流动系统（流场）地面就可以测到一个相匹配的瞬变电磁平面波，二者是正投影关系。探测仪器的任务就是探测出瞬变电磁波场在地表的分布形状，在这个场内确定地下水源的位置及深度。地磁北极为S，南极为N，它与地理南北极有一个夹角，在我国大概偏西5度左右，我国的断裂方向以北东向为主导方向，使断裂带中的地下水流与磁力线成较大的交角，地下水是良导体，它在裂隙溶洞中按水文学方式运动就会切割磁力线，按法拉第原理，水流中产生电流，而电流周围又产生感应磁场，这个作用连续不断地进行地磁场干扰，结果就从地下水所处的深度向地面发射出瞬变电磁波，按照麦克斯韦定律，这种波的振动方向与前进方向垂直，属平面波，可直达地面。因此，在地面探测出瞬变电磁波的分布场，就知道地下是否存在地下水流动系统，通过改变测试频率可知瞬变电磁波场源的深度，这样就可以确定地下水流动系统实际分布情况。在图2中，假设地下径流动态信息为单一频率f _b的正弦波、天然电场信号为单一频率f _c的正弦波，地面探测仪器接收到的是f _b与f _c的耦合信号，即是经过f _b调制后的f _c信号；根据信号调制原理：在该信号的频域波形中，中间最高的谱线称为主频谱线f _c，分布于主频谱线两边并关于主频对称的谱线称为边频谱线f _b。主频谱线f _c表征地下某深度岩层的电性值，边频谱线f _b表征地下水的活动情况，边频谱线幅值表示地下水瞬时流量的大小，边频谱线距离主频谱线的远近称边频频率，边频频率表示地下水瞬时流动速度的快慢。地下径流动态信息即为附图2中所示的边频，边频具有以下特性：（1）边频关于主频对称，（2）边频幅值较主频幅值小，（3）边频、主频均是一定频谱范围的相对极值点，（4）边频存在1、2、…、n阶。

附图2所示为地下径流动态信息的理想状态，附图3所示是较为规则的地下径流动态信息。通常，地下径流在岩溶裂隙中流动的流速、流量都是随时间变化的，其产生的地下径流动态信息也是不规则的；作为地下径流动态信息载体的天然电磁波受场源、传输途径、天气等诸多因素的影响，也是变化的、不规则的；另探测仪器在现场采集数据的过程中，会受周边电力、通讯设施和厂矿大型设备运行产生的电磁干扰影响，导致探测数据的不稳定；这就给地下径流动态信息的准确识别带来很大的难度。本发明为克服以上问题，通过对大量的设施例分析，总结了以下的识别方法。

附图3～附图10均是本发明实施例的细化谱图，图中上半部分为天然瞬变电磁波的时域波形，下半部分为天然瞬变电磁波中指定分析带宽f ₁～f ₂的细化频谱图。频谱图中间的f _c谱线称为主频谱线，分布于主频谱线两边并关于主频对称的谱线f _b称为边频谱线。根据电磁场理论，天然瞬变电磁波在地层中的传导电流远远大于位移电流，则其穿透深度为

h＝1/(π f _c μ/ρ)^0.5 (单位：m)

式中，ρ为探测点的地表电阻率，一般地表大都是松散的土层，取其平均电阻率ρ≈10（Ω .m）；假设天然瞬变电磁波穿透的地层为无磁性介质，则磁导率取 μ＝4π×10^-7 H/m；据此得到地下水径流的埋藏深度与主频频率的关系如下

h＝1/(π f _c μ/ρ)^0.5＝1591.58(1/f _c)^0.5(单位：m)

流动的地下水会切割地磁场的磁力线而产生感应电动势e，根据电磁流量原理，感应电动势e（单位：mV）与流量的方程式为

e＝（4BK/L）Q

式中，B为探测地点的大地电磁场感应强度（单位：mV/m ²）、K为常数，对于一个已知的探测区域来说它们都是常数，但需要在探测区域附近找已知的水量的参考测试点进行整定；L为地下水流通道截面的周长（单位：m）。实际探测工作表明，地下水在岩溶裂隙中流动的流速、流量都是随时间变化的，其产生的感应电动势e是时间t的函数，计为e (t)＝0.5〔e ₁(t)＋e ₂(t)〕。e (t)与地下水储量Sw之间的关系为

Sw=∫₀ ^T[L.e (t)/4BK]dt=( L/4BK)∫₀ ^T e (t) dt (单位：m ³/h)

式中，T为地下水径流的变化周期，探测仪器能探测出地下水径流的变化周期T（边频频率f _b的倒数）和周期T内感应电动势e (t)的综合值，由上式可简单地估算出地下水储量Sw。但由于地下水径流的变化周期T受地质构造、补给情况、开采量等因素的影响而各不相同，且该周期也是随时间变化的，这给准确评估地下水储量或开采量带来一定困难。通常采用离散化方法来计算上式，具体做法是：由探测仪器对该固定深度连续重复探测1小时，得到M组探测数据，而这其中只有N组探测数据具有地下水径流特征信息，即边频综合幅值e (t ₀)、e (t ₁)、e (t ₂)……e (t _N-1)，边频频率f _b0、f _b1、f _b2……f _bN-1；因此，上式的离散化计算公式为：

(单位：m ³/h)

因此，本发明的目的是有效识别主频f _c的频率与幅值、边频f _b的频率与幅值，其具体做法及步骤参见附图1：探测仪器接收天然电场信号并细化得到分析带宽f ₁～f ₂的分辨率为∆f的256条谱线，经最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据、动态信息识别判据后得出有无地下径流动态信息的结论并相应给出主频信息、动态信息的幅值、频率。

最大谱线离散判据的理论依据是信号调制原理，在调制信号频谱中基频f _c的幅值是最大且离散的，据此判断探测信号中是否存在调制现象。本发明的具体做法是在探测仪器接收天然电场信号并细化得到的256条谱线中，首先找到最大幅值谱线F_MAX0并缓存其频率及幅值，然后计算其及前后各3条谱线一起的7条谱线的平均幅值F_J0，若幅值F_MAX0≥3F_J0则最大幅值谱线F_MAX0离散，反之则最大幅值谱线F_MAX0不离散。

频谱能量集中判据的理论依据是电磁波传输特性，在电磁波的传输过程中，高频率的电磁波信号可传输得更远且不易衰减，而低频率的电磁波信号则传输得不远且容易衰减；这样，调制信号频谱中的主要能量是高频率电磁波信号（载波信号或基波）贡献的，据此可评价调制信号频谱受干扰的程度。本发明判据是在最大谱线离散判据的基础上提出的，其具体做法是首先探测并细化所得的256条谱线的平均幅值F_J，然后与最大谱线及其前后各3条谱线的平均幅值F_J0比较，计算若F_J0≥10F_J则所分析频谱的能量集中，否则分析频谱的能量发散。

极点谱线搜索判据的理论依据是调制信号的频谱分配原则及边频特性，目的是去除分析频谱中的不相关谱线。本发明的具体做法是：（1）将原频谱的256条谱线值复制到一指定存储区缓存；（2）在缓存区中，按从前往后两两比较幅值低者置0（后向积分）、从后往前两两比较幅值低者置0（前向积分）的原则，使缓存区中只保留极值点的谱线值，其余无关的谱线值均被置0；（3）最大幅值谱线F_MAX0肯定是极值点，其离散性已在最大谱线离散判据中判断，接着在（2）得到的极值谱线中找到第二大极值点F_MAX1、第三大极值点F_MAX2（其中，F_MAX1幅值> F_MAX2幅值），并结合原频谱计算次大极值点谱线F_MAX1、F_MAX2及其前后各3条谱线一起的平均幅值F_J1、F_J2，若F_MAX1≥3F_J1则F_MAX1离散、F_MAX2≥3F_J2则F_MAX2离散，缓存F_MAX1、F_MAX2的离散性判断结果及其频率、幅值；（4）将幅值小于最大幅值谱线F_MAX0幅值的三十分之一的极值谱线置0，进一步减少无关谱线以提高后续识别搜索的效率。

经过最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据后，接下来由动态信息识别判据根据以上三判据的结果做出分析频谱中有无地下径流动态信息的结论；若存在动态信息，则相应给出主频信息、动态信息的幅值、频率；若不存在动态信息，只给出主频信息的幅值、频率。地下径流动态信息即为分析频谱中关于主频的边频，同样具有边频的特性：①边频关于主频对称，②边频幅值较主频幅值小，③边频、主频均是一定频谱范围的相对极值点，④边频存在1、2、…、n阶；因此，本发明的动态信息识别判据是根据边频特性提出的，其具体做法是：若分析频谱的频谱能量集中，且3≤极值点个数≤30，则通过边频识别规则一、规则二、规则三判断是否存在关于F_MAX0、F_MAX1及F_MAX2的对称极值谱线来确定地下径流动态信息的存在，否则主频信息为F_MAX0、动态信息为0且退出识别；所谓的对称极值谱线是指两条极值谱线中间存在另一条极值谱线使两条极值谱线与其点距之差在±2条谱线范围内，称该对称极值谱线为边频，其频率为一阶边频频率f _b、幅值为两条谱线的幅值均值。针对边频在频谱中的分布情况不同，本发明分三种方式进行识别，即边频识别规则一、规则二、规则三。

边频识别规则一：F_MAX0离散且极值点分布于其左、右两边，若存在关于F_MAX0对称的边频极值谱线，则主频信息为F_MAX0，动态信息为一阶边频的信息；若不存在关于F_MAX0对称的边频极值谱线、但存在关于F_MAX1或F_MAX2对称的边频极值谱线且相应的F_MAX1或F_MAX2离散，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。在附图3中，最大幅值谱线F_MAX0的左边和右边各存在一个关于F_MAX0对称的极值点，该极值点对应的谱线即为一阶边频谱线，与附图2中所示的理想情况非常相似；因此，主频（即载波信号）信息为F_MAX0谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.5×左边频幅值＋0.5×右边频幅值）。在附图4中，明显可见最大幅值谱线F_MAX0的左边存在二阶边频极值点、右边存在三阶边频极值点，且它们都关于F_MAX0谱线对称；因此，主频信息为F_MAX0谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.5×左一阶边频幅值＋0.5×右一阶边频幅值）。在附图5中，虽然在最大幅值谱线F_MAX0的左、右两边都存在极值点，但却不关于F_MAX0谱线对称，图中明显可见次大幅值谱线F_MAX1的左边存在一阶边频极值点、右边存在二阶边频极值点，且关于F_MAX1谱线对称，左边的二阶边频极值点重合在最大幅值谱线F_MAX0的频带内；因此，最大幅值谱线F_MAX0应是干扰或其他频率成分的谐波所致，取主频信息为F_MAX1谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.5×左一阶边频幅值＋0.5×右一阶边频幅值）。

边频识别规则二：F_MAX0离散但主要极值点分布在其一边，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。在附图6中，次大幅值谱线F_MAX1的左边是最大幅值谱线F_MAX0、右边存在一个边频极值点，且它们关于F_MAX1对称；因此，最大幅值谱线F_MAX0应是干扰或其他频率成分的谐波所致且频带内包含左一阶边频，取主频信息为F_MAX1谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.25×左一阶边频幅值＋0.75×右一阶边频幅值）。在附图7中，次大幅值谱线F_MAX1的左、右两边都存在二阶边频极值点，且它们关于F_MAX1谱线对称，而最大幅值谱线F_MAX0与它们均不相关；因此，最大幅值谱线F_MAX0应是干扰或其他频率成分的谐波所致，取主频信息为F_MAX1谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝＝0.5×左一阶边频幅值＋0.5×右一阶边频幅值）。在附图8中，次大幅值谱线F_MAX1、F_MAX2都分布在最大幅值谱线F_MAX0的左边，而次大幅值谱线F_MAX1不离散、F_MAX2离散，且F_MAX2的右边存在一个极值点与左边的F_MAX1关于F_MAX2谱线对称，如果包含F_MAX0、F_MAX1在内，则在F_MAX2的左、右两边都存在二阶边频极值点，且关于F_MAX2谱线对称；因此，最大幅值谱线F_MAX0应是干扰或其他频率成分的谐波所致且将F_MAX2的右二阶边频包含在其频带内，次大幅值谱线F_MAX1也是干扰或其他频率成分的谐波所致且将F_MAX2的左一阶边频包含在其频带内，取主频信息为F_MAX2谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.25×左一阶边频幅值＋0.75×右一阶边频幅值）。

边频识别规则三：F_MAX0不离散，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；若F_MAX1和F_MAX2均不离散但存在包括F_MAX0、F_MAX1、F_MAX2中至少2个在内的4个以上的点距误差在±2条谱线内的极值谱线，则主频信息的频率为F_MAX0谱线频率、幅值为F_MAX0与F_MAX1的均值，动态信息的频率为一阶边频频率、幅值为F_MAX1与F_MAX2的均值；否则主频信息为F_MAX0、动态信息为0且退出识别。在附图9中，最大幅值谱线F_MAX0不离散、但次大幅值谱线F_MAX1离散，在F_MAX1的左边存在一阶边频极值点、右边存在包含F_MAX0在内的二阶以上边频极值点，且关于F_MAX1谱线对称；因此，最大幅值谱线F_MAX0应是干扰或其他频率成分的谐波所致且将F_MAX1的右一阶边频包含在其频带内，取主频信息为F_MAX2谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（可从频谱中直接读取），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.75×左一阶边频幅值＋0.25×右一阶边频幅值）。在附图10中，最大幅值谱线F_MAX0与次大幅值谱线F_MAX1、F_MAX2均不离散，但存在包括F_MAX0、F_MAX1、F_MAX2在内的4个以上的点距误差在±2条谱线内的极值谱线，且关于最大幅值谱线F_MAX0对称；因此，频谱应是受到较严重的干扰所致，取主频信息的频率为F_MAX0谱线所对应的频率f _c（＝f ₁＋f ₁与f _c的点距×∆f）、幅值（＝0.5×F_MAX0幅值＋0.5×F_MAX1幅值），动态信息为一阶边频的频率f _b（＝0.5×两个一阶边频的点距×∆f）、幅值（＝0.5×F_MAX1幅值＋0.5×F_MAX2幅值）。

综上所述，本发明从研究地下径流动态信息的形成机理出发，运用信号调制原理与边频特性相结合的原则在细化频谱内自动识别地下径流动态信息，并充分考虑了干扰因素的影响，可使地下径流动态信息识别准确率达到90%以上，解决了现有物探方法中难以识别地下径流动态信息的问题。但地下径流是一个完整的地下水流动系统，其变化会因时间、地点的不同而不同；虽然其基本规律、识别原则是一致的，但还存在一定的差异；故还需更多的经验总结来完善本发明的识别方法。

Claims

1.通过探测仪器接收天然电场信号并细化得到分析带宽f ₁～f ₂的分辨率为∆f的256条谱线，经最大谱线离散判据、频谱能量集中判据、极点谱线搜索判据、动态信息识别判据后得出有无地下径流动态信息的结论并相应给出主频信息、动态信息的幅值、频率。

2.根据权利要求1所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的最大谱线离散判据是在256条谱线中找到最大幅值谱线F_MAX0并计算其及前后各3条谱线一起的平均幅值F_J0，若F_MAX0≥3F_J0则离散，反之不离散。

3.根据权利要求1所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的频谱能量集中判据是计算256条谱线的平均幅值F_J，若F_J0≥10F_J则频谱能量集中，否则频谱能量发散。

4.根据权利要求1所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的极点谱线搜索判据是缓存256条谱线到一指定存储区，在其中按从前往后两两比较幅值低者置0、从后往前两两比较幅值低者置0的原则找出频谱中的所有极值点，结合原频谱计算次大极值点谱线F_MAX1、F_MAX2及其前后各3条谱线一起的平均幅值F_J1、F_J2，若F_MAX1≥3F_J1则F_MAX1离散、F_MAX2≥3F_J2则F_MAX2离散，然后再将幅值小于F_MAX0/30的极值点置0。

5.根据权利要求1所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的动态信息识别判据是若频谱能量集中且3≤极值点个数≤30，则通过边频识别规则一、规则二、规则三判断是否存在关于F_MAX0、F_MAX1及F_MAX2的对称极值谱线来确定地下径流动态信息的存在；否则主频信息为F_MAX0、动态信息为0且退出识别。

6.根据权利要求5所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的对称极值谱线是指两条极值谱线中间存在另一条极值谱线使两条极值谱线与其点距之差在±2条谱线范围内，称该对称极值谱线为边频，其频率为一阶边频频率f _b、幅值为两条谱线的幅值均值。

7.根据权利要求5所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的边频识别规则一是F_MAX0离散且极值点分布于其左、右两边，若存在关于F_MAX0对称的边频极值谱线，则主频信息为F_MAX0，动态信息为一阶边频的信息；若不存在关于F_MAX0对称的边频极值谱线、但存在关于F_MAX1或F_MAX2对称的边频极值谱线且相应的F_MAX1或F_MAX2离散，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。

8.根据权利要求5所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的边频识别规则二是F_MAX0离散但主要极值点分布在其一边，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；否则主频信息为F_MAX0、动态信息为0且退出识别。

9.根据权利要求5所述的一种地下径流动态信息的频谱识别方法，其特征在于：所述的边频识别规则三是F_MAX0不离散，若F_MAX1或F_MAX2离散且存在关于F_MAX1或F_MAX2对称的边频极值谱线，则主频信息为F_MAX1或F_MAX2，动态信息为一阶边频的信息；若F_MAX1和F_MAX2均不离散但存在包括F_MAX0、F_MAX1、F_MAX2中至少2个在内的4个以上的点距误差在±2条谱线内的极值谱线，则主频信息的频率为F_MAX0谱线频率、幅值为F_MAX0与F_MAX1的均值，动态信息的频率为一阶边频频率、幅值为F_MAX1与F_MAX2的均值；否则主频信息为F_MAX0、动态信息为0且退出识别。