CN101802650A - 用于实时监测地球表面风暴活动的方法和装置 - Google Patents

用于实时监测地球表面风暴活动的方法和装置 Download PDF

Info

Publication number
CN101802650A
CN101802650A CN200880019493A CN200880019493A CN101802650A CN 101802650 A CN101802650 A CN 101802650A CN 200880019493 A CN200880019493 A CN 200880019493A CN 200880019493 A CN200880019493 A CN 200880019493A CN 101802650 A CN101802650 A CN 101802650A
Authority
CN
China
Prior art keywords
parameter
elf
earth
signal
spark discharge
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.)
Granted
Application number
CN200880019493A
Other languages
English (en)
Other versions
CN101802650B (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.)
Uniwersytet Jagiellonski
Original Assignee
Uniwersytet Jagiellonski
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 PL382624A external-priority patent/PL382624A1/pl
Priority claimed from PL385320A external-priority patent/PL385320A3/pl
Application filed by Uniwersytet Jagiellonski filed Critical Uniwersytet Jagiellonski
Publication of CN101802650A publication Critical patent/CN101802650A/zh
Application granted granted Critical
Publication of CN101802650B publication Critical patent/CN101802650B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/16Measuring atmospheric potential differences, e.g. due to electrical charges in clouds
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

一种包括记录ELF(超低频)电磁场信号的方法,ELF电磁场信号由电火花放电产生并在地球电离层谐振器中诱发,基于地球电离层腔中的ELF场谐振传播模型对ELF电磁场信号进行频谱分析,ELF电磁场信号的依赖于电火花放电的特征、距离和强度被分离并且然后与数据库模型参数比较。被选作最佳近似的模型数据能够重构电火花放电的位置和强度。一种能够重构电火花放电的位置和强度的装置具备两个感应磁天线。

Description

用于实时监测地球表面风暴活动的方法和装置
技术领域
本发明的主题是一种用于实时监测地球表面风暴活动的方法和装置,其用于远程测量风暴中心活动以及确定它们的位置。风暴活动的监测与航空和海洋运输的安全相关,并且在天气预报中扮演重要角色。
背景技术
当前,有多种方法来记录发生在雷雨云中、云之间和/或云和地面之间的大气中的电火花放电。
第一种监测风暴活动的方法是基于光学卫星对闪电的观测,闪电发生在前述的雷雨云中、云之间和/或云和地面之间。该方法提供由安装在卫星上的测量仪器直接观测的范围内的详细信息。所观察的区域覆盖ca.105km2,其只构成地球表面的一小部分。即使是具有卫星广阔扩展的网络的轨道观测系统也不能确保可靠的风暴活动全球监测系统。
追踪风暴的第二种方法与在大气中电火花放电期间产生的VLF、HF和VHF无线电频率内的电磁信号的探测有关。这种方法提供放电强度及其位置的评估。然而,其缺陷是在这些频率内的信号的范围,在HF内被限制到几百公里的距离,而在VLF内被限制到几千公里的距离。覆盖大范围以提供利用基站的观测需要密集的追踪基站网络。目前,最多开发的监测系统执行同时寄存在VLF和HF频率内的信号分析。这只在高度发达的国家中操作,因此其只覆盖地球表面一小部分。
第三种方法是基于只追踪VLF内的信号。单个基站的范围达到一千公里。当前,全球系统基于27个测量基站。
极低频(ELF)电磁场信号的传播从以“Studies of ELF propagation in the sphericalshell cavity using a field decomposition method based on asymmetry of Schumannresonance curves”为论文名发表在Journal of Geophysical Research,Vol.111,A 10304,doi:10.1029/2005JA011429,2006上的已公开的波兰克拉科夫的雅盖隆大学(JagiellonianUniversity)的科学家主导的研究成果得知。与该论文提出的论题相适应的,所观测的ELF谱中的曲线非对称和峰值谐振频率的变化由驻波场的叠加产生,这将与传出源头的波场一起产生谐振模型。所要解决的问题是是否能够分离谐振器中的两个场分量(component)并且单独地测量它们。关于这一问题的新的解决方式已经被提出并且其存在于频谱不对称性的测量或通过利用单个天线的电场或磁场分量的观测所获得的信号频谱。在这种方法中,假定谐振器的任意点中的信号频谱包括与谐振分量场相关的对称部分和与传播波场联系的非对称部分。场分量的功率谱α(θ,f)|2已经利用下述公式确定:
Figure GPA00001011160300021
该公式实现了谐振器的任意点到单个前进波源的近似距离θ的确定。
上述大气放电的观测方法,即使假定观测基线相当大的扩展,也不能实现覆盖地球整个表面的风暴活动的全球监测。同样,上述方法没有一个能够保证大气中的电火花放电的100%的测量效率。目前估计上述方法的效力介于60%和80%之间,这取决于所采用的方法以及所使用的信号分析算法。
发明内容
本发明的目的是提供一种横跨整个地球表面的风暴活动的全球监测的方法和装置。
根据本发明的监测风暴活动的方法是基于风暴单元中的大气的电火花放电诱发的电磁信号的分析,这类似于上述方法。本发明的想法是记录地球电离层谐振器内诱发的极低频(ELF)电磁场信号,然后基于涉及地球电离层谐振器中的ELF场的谐振传播模型进行频谱分析。接着,取决于电火花放电位置和天线之间的距离,即源头距离,和强度的特征被分开然后与模型数据库中的参数作比较。最接近的基准参数组确定这些源的位置和强度。分析可以实时处理,并且其结果可以以风暴活动图的形式图解。
由于ELF场内的波谐振传播的特殊性质,单个装置能够记录来自地球大气中的电火花放电的信号,并且寄存信号能够实现直到10000km距离的风暴活动图的明确预备。位于不同大陆的几个装置的使用允许在全球范围以地理坐标的形式轮廓清晰成像风暴活动。
优选地,利用两个彼此垂直设置,位于表面上或密切位于地表之下,并且沿着NS和EW方向放置的水平磁天线将ELF电磁场信号记录在没有局部电场源的区域中。该信号,被放大、滤波和例如在180Hz的取样频率变成数字观测信号的处理之后,通过无线电发射到电子数据处理单元,在那里进行频谱分析。所产生的信号功率谱然后与谐振曲线相配,这允许参数化这些频谱。接着,这些参数与存储在模型数据库中的基准参数相比较,该模型数据库已经利用地球电离层谐振器或腔或波导中的ELF波传播模型的知识产生。利用程序以最小化观测参数相对于那些基准参数之间的偏差,相对于一组基准参数作出选择,这组基准参数的值是最好的近似。一组基准参数重现了发生了大气中的电火花放电的风暴源头的位置和强度。根据这组基准参数,能够预备源图。
为了执行每个功率谱Sx,Sy的参数化,推荐利用下述公式匹配谐振曲线:
S ( ω ) = b + a ω α + Σ k = 1 K p k · [ 1 + e k · ( ω k - ω ) ] ( ω k - ω ) 2 + ( Γ k ) 2
其中:
S(ω)-匹配功率谱;
a-描述背景颜色噪声的参数;
b-描述背景宽带噪声的参数;
α-颜色噪声的频谱指数;
ω-比率(脉动);
pk-第k个谐振峰的最大功率值;
ek-第k个谐振峰的非对称参数;
ωk-第k个谐振峰的谐振比率,其等于2πfk
Γk-第k个谐振峰的半宽度。
本发明的另一个想法是实时监测地球表面的风暴活动的装置,该装置具有两个感应的、有源磁天线,该两个天线彼此沿着NS和EW方向垂直并放置在地表上或密切位于地表下,并且彼此通过屏蔽线连接到记录观测信号的块,利用两个相同的信号通道和放大器,滤波器,模数转换器和控制系统,其时针通过系统接收器与GPS时间同步,并且其中的数据电子处理单元是根据运用到ELF信号的频谱分析算法和取决于风暴放电源的距离和活动确定特征的算法编程。
优选地,记录观测信号的块连接到发射在模数转换器内处理的数据的发射天线,而计算单元连接到与连接到记录观测信号的块的发射天线通讯的接收天线。
附图说明
下面将通过例子并参考附图描述本发明,其中:
图1示出了产生ELF电磁波的风暴中心和地球;
图2示出了监测风暴活动的方法的流程图;
图3A示出了ELF电磁场的接收天线;
图3B和3C示出了实时监测地表上风暴活动的装置;以及
图4示出了三个选择的与第一三个连续传播极大值相关的谐振峰。
具体实施方式
附图中的实施例只是打算用于说明的目的并且不限制附加的权利要求所限定的本发明的范围。
图1-5示出了一种实时监测地球表面上风暴活动的装置,根据本发明的解决方法将作为用来远程测量风暴中心活动以及其在地球表面上的位置的风暴活动监测装置的例子呈现。装置的实现相同功能的部件和单元在所有的附图中以相同的编号标记或者仅以首数字不同的编号标记,其将分配给特殊的附图或数字。
图1示意性地示出了具有风暴中心101的地球100。风暴中心在空气腔104中,或者换句话说在地球电离层谐振器中,产生ELF电磁场信号,该信号是在地球电离层腔104内的所有方向传播的ELF电磁波102和103。在图1中,示出了单个大气电火花放电109,产生在电离层105和地面106之间产生的空气腔104中逆时针传播的第一电磁波102和顺时针传播的第二电磁波103。这些波环绕地球许多次,彼此干涉并且也到达以彼此相对垂直的有利位置布置的接收天线110、120。由于风暴中心101产生的波的干涉,特殊频率的波衰落并且以这种方式产生一系列谐振峰。
图2示意性地示出了实时监测地球表面上的风暴活动的方法的流程图。在步骤211,电火花放电产生的电磁场分量Hx,Hy产生信号,该信号由两个磁天线寄存。在步骤212,这些信号经历相同的频谱分析处理,结果产生功率谱Sx和Sy。由于地球电离层腔内的ELF信号传播的谐振特性,频谱遵循一系列与随后的传播极大值相关的谐振峰的形式。在步骤213,每个功率谱Sx和Sy与谐振曲线匹配,这将通过稍后详细描述的特殊公式描述。匹配程序的结果是频谱的参数化。每个频谱Sx和Sy被指派具有一组观测参数pk、fk、γk和ek。在步骤214,每个信号的参数组与基准参数pok、fok、γok和eok比较,基准参数pok、fok、γok和eok存储在模型参数数据库B中。比较系统的任务是一组基准参数pok、fok、γok和eok的选择,这组基准参数pok、fok、γok和eok的值尽可能接近一组观测参数pk、fk、γk和ek。为了这一目的,使用一种可用的最小化程序,用来最小化观测参数pk、fk、γk和ek和基准参数pok、fok、γok和eok之间的偏差。模型参数数据库B是,例如,在所观测的大气中的电火花放电的证实之后,在步骤215中利用地球电离层腔中ELF波传播模型知识制作的。由于对于观测者,即天线,和源头,即电火花放电,之间的每一个距离,以及对于每个放电强度都有特定的基准参数组,因此将观测参数组pk、fk、γk和ek指派给基准参数组pok、fok、γok和eok能够绘制源头的位置及其强度。在步骤216,准备源头图。所获得图的分辨能力严格地取决于信号分解的步骤213中获得的频谱观测参数组pk、fk、γk和ek的数目以及伴随测量的局部噪声等级。观测参数组pk、fk、γk和ek是取决于电火花放电的距离和强度的特征,其有助于重构电火花放电的位置和强度以及形成能够绘制风暴源头(M(x, y))图的基础。图的不确定性严格取决于场分布基准组的形式,即它们是假定传播模型的衍生物(derivative)。
图3A中示意性示出的在空气腔中产生的ELF电磁场301具有分量Hx,Hy和Hz,其中Hx和Hy由两个理想地彼此垂直的并且放置在没有局部电场源的地球区域的接收水平定向天线Ax321和Ay322接收。通常,天线Ax321和Ay322是感应的、有源的磁天线,例如,其是1m长,具有1cm2的芯截面。优选地,天线Ax321和Ay322是屏蔽外来电场(即使有的话)的,并且放置在地表106之上或密切地位于其下,通常沿着NS和EW方向。由于ELF电磁场301,在这些定向接收天线Ax321和Ay322中产生信号Ux331和Uy332。
在信号Ux331和Uy332经历基于地球电离层波导或地球电离层腔,简言之,空气腔中的ELF场的谐振传播模型的频谱分析,以及距离和由其决定的源头强度的特征分离之前,定向磁接收天线Ax321和Ay322中产生的信号Ux331和Uy332被发射到在两个相同的处理路径记录观测信号的组块(block),该两个相同的处理路径优选的是分开的,如图3B和3C所示,通过例如100m长度的屏蔽线发射到宽带低噪声放大器Gx341和Gy342的输入端,然后被放大,之后它们穿过例如通带1-60Hz的低通去阶梯效应滤波器BPFx351和BPFy352,接着被从模拟形式转换到数字形式,例如在与控制系统uP 365连接的16位模数转换器A/Cx361和A/Cy362中执行。信号的取样在具有例如180Hz频率的两个模数转换器A/Cx361和A/Cy362中同时执行。具有两个上述和下述组块的装置具有GPS接收器364,其用来使控制系统时钟与GPS时间同步。数字信号Ux331和Uy332然后被编码和引导到具有发射在模数转换器中处理的数据的天线T 367的低功率无线电发射系统366,或被编码和引导到波导线。另一种解决方法也是可能的,其中数字信号Ux331和Uy332可以被发射到控制系统以进一步分析,例如通过通讯线路发射到PC。然而,通过这种解决方法,引入噪声到测量路径的危险可能发生。在图3B所呈现的有利解决方法中,目前描述的所有装置块都远离电线放置并且都由自发电源供电,而到计算单元,例如上述的PC,的信号发射通过无线电或光发生。
在所呈现的例子中,数字信号Ux331和Uy332由远离天线几公里,例如3-10公里,并具有接收天线377的无线接收系统376接收,然后被编码和发送到计算单元375,例如由总线供电的PC。根据ELF信号频谱分析算法和与依赖电火花风暴放电源的距离和活动特征的确定相关的算法来编程作为电子处理系统的计算单元375。在计算单元375中,数字信号Ux331和Uy332经历相同的频谱分析处理以产生观测信号Ux231和Uy231的功率谱Sx和Sy并利用块DFTx371和DFTy372中的傅立叶变换使它们与谐振曲线匹配。由于地球电离层腔中的ELF观测信号传播的谐振特性,频谱是与表格1以及图4中的图表形式所呈现的连续传播极大值相联系的谐振峰系列或峰411、412和413的形式。
表格1
  Mod数k   Mod频率f   Mod功率p   Mod离心率e   Mod宽度Γ
  1234…   f1f2f3f4...   p1p2p3p4   e1e2e3e4   γ1γ2γ3γ4
  abα
在接下来的阶段,在匹配块Dx381和Dy382中,下述公式表示的谐振曲线被与每个功率谱Sx和Sy相匹配:
S ( ω ) = b + a ω α + Σ k = 1 K p k · [ 1 + e k · ( ω k - ω ) ] ( ω k - ω ) 2 + ( Γ k ) 2
其中:
S(ω)-匹配功率谱;
a-描述背景颜色噪声的参数;
b-描述背景宽带噪声的参数;
α-颜色噪声的频谱指数;
ω-比率(脉动);
pk-第k个谐振峰的最大功率值;
ek-第k个谐振峰的非对称参数;
ωk-第k个谐振峰的谐振比率,其等于2πfk
Γk-第k个谐振峰的半宽度。
上述公式是采用自发表在Journal of Geophysical Research,Vol.111,A10304,doi:10.1029/2005JA011429,2006上论文名为“Studies of ELF propagation in the sphericalshell cavity using a field decomposition method based on asymmetry of Schumannresonance curves”的文献,在利用描述背景颜色噪声的参数,描述背景宽带噪声的参数和颜色噪声的频谱指数完成其之后。增加描述背景颜色噪声的参数a,描述背景宽带噪声的参数b和颜色噪声的频谱指数α是更优选的,因为,如上所述,其改进了将曲线匹配到所观测背景的质量并增加了产生地图所需的频谱参数的确定的精度。需要重视的是参数a和b以及参数α是被确定的,但是不用于产生地图。
作为运用匹配程序的结果,频谱参数化随后发生在参数化块Paramx391和Paramy392中。对于每个频谱Sx和Sy,四个当前观测参数pk、fk、γk和ek被用于每个第k个模式,用于实现其参数化。在接下来的阶段,每个信号Ux和Uy的这四个观测参数pk、fk、γk和ek在比较器396中与基于地球电离层腔内的ELF波传播的知识产生的存储在模型参数数据库B 390中的基准参数组pok、fok、γok和eok比较。比较器396的任务是选择值最接近观测参数组pk、fk、γk和ek的一组基准参数pok、fok、γok和eok。为了这一目的,采用程序中最小化参数pk、fk、γk和ek与参数pok、fok、γok和eok间的偏差的一个程序。模型参数数据库B 390凭借地球电离层腔内的ELF波的传播知识的实力根据长期观测而创建。由于对于观测者和源头之间的每个距离以及每个放电强度都具有特定的参数组,所以将观测参数组pk、fk、γk和ek指派到基准参数组pok、fok、γok和eok能够重构源头的位置及其强度。换句话说,比较观测参数组pk、fk、γk、ek和基准参数组pok、fok、γok、eok,基准参数组pok、fok、γok、eok是模型参数数据库B的参数。基准参数组pok、fok、γok、eok被选择来作为观测参数组pk、fk、γk、ek的最佳近似之后,重构电火花放电或放电源和/或风暴中心的位置和强度。优选例子中的电子处理单元已经根据ELF信号频谱分析算法和借助于傅立叶变换确定取决于风暴放电的距离和活动的特征的算法编程,因为这在频谱分析的许多论文中都是已知的。
在接下来的阶段,源头图M(x,y)被展现,或者换句话说,在产生块395中,地球表面上风暴中心的分布和强度的源头图M(x,y)被展现。所获得的图的分辨率直接取决于在信号分解阶段获得的频谱分布参数的数目和在测量期间发生的噪声等级。图的不确定性(ambiguity)取决于场分布基准组的形式,即它们是假定传播模型的偏差。
已经被其他源头证实和/或通过其他监测地球表面上风暴活动的方法证实的观测参数组pk、fk、γk和ek能够被加入模型数据库B 390,并且用这种方法构成增补到用于连续重构放电源或风暴中心的位置和强度的基准参数组pok、fok、γok和eok
在当使用显著地增加了地球表面图上风暴中心的位置和强度的分辨能力的许多接收装置的情况下,能够通过卫星或因特网链接将数据发送到数据收集中心,在此该数据经历计算单元中的分析。
利用选择的例子呈现了根据本发明的解决方法。然而,这些例子并不限定本发明。显而易见的是在不改变该解决方法的本质特点下能够引入改进。所呈现的例子不是根据本发明的解决方法的应用的唯一的可能性。

Claims (8)

1.一种利用风暴单元中的大气的电火花放电诱发的电磁信号的分析实时监测地球表面上风暴活动的方法,其特征在于:
记录地球电离层谐振器中产生的ELF电磁场信号,基于所述地球电离层谐振器中的ELF场谐振传播模型对所述ELF电磁场信号进行频谱分析,其中取决于电火花放电的距离和强度的特征被分离并被与包含在模型参数数据库中的基准数据作比较,并且被选择作为取决于所述电火花放电的距离和强度的特征的最佳近似的所述基准参数用于开发绘制所述电火花放电的位置和强度图。
2.根据权利要求1所述的方法,其特征在于所述取决于所述电火花放电的距离和强度的特征是与存储在所述模型参数数据库中的基准参数(pok、fok、γok和eok)相比较的观测参数(pk、fk、γk和ek),所述基准参数(pok、fok、γok和eok)重构所述电火花放电的位置和强度并形成绘制风暴源头(M(x,y))图的基础。
3.根据权利要求1或2所述的方法,其特征在于通过两个有源磁性感应天线(Hx,Hy)将由所述电火花放电诱发的ELF电场信号记录在地球区域中,所述地球区域没有局部的电场源,所述两个天线(Hx,Hy)水平地放置并且优选地彼此垂直,并且沿着NS和EW方向位于地球表面上或密切地位于地表之下,然后这些已经被放大、滤波并在180Hz取样频率处被处理成数字观测信号(Ux,Uy)的信号被以无线电发射到电子数据处理单元,其中它们随后被频谱分析并产生匹配谐振曲线的功率谱(Sx,Sy),并且这些频谱被参数化以发现一组观测参数(pk、fk、γk和ek),并且接下来,所述观测参数(pk、fk、γk和ek)与利用地球电离层腔中的ELF波传播模型的知识产生的存储在所述模型参数数据库中的基准参数(pok、fok、γok和eok)相比较,在此之后,利用关于观测参数(pk、fk、γk和ek)和基准参数(pok、fok、γok和eok)之间的偏差的最小化程序,选择一组基准参数(pok、fok、γok和eok),该组基准参数(pok、fok、γok和eok)的值是观测参数(pk、fk、γk和ek)的最佳近似,其中所述一组基准参数(pok、fok、γok和eok)重构所述电火花放电的位置和强度并形成绘制风暴源(M(x,y))图的基础。
4.根据权利要求1所述的方法,其特征在于在记录所述地球电离层谐振器中诱发的ELF电磁场信号之后,基于所述地球电离层谐振腔中的ELF场谐振传播模型对所述ELF电磁场信号进行频谱分析,由此确定的源头距离以及强度特征被从所述ELF电磁场信号分离并与存储在所述模型参数数据库中的基准参数(pok、fok、γok和eok)相比较,并且其中被选作观测参数(pk、fk、γk和ek)的最佳近似的那些基准参数(pok、fok、γok和eok)重构所述电火花放电的位置及强度。
5.根据权利要求1或2所述的方法,其特征在于通过至少两个水平天线将所述ELF电场信号记录在地球区域中,所述地球区域没有局部的电场,所述天线优选地彼此垂直地放置并且位于地表上或密切地位于地表之下,其中所述ELF电磁场信号在被放大、滤波并在特定取样频率处被处理成数字观测信号之后,被以无线电发射到电子数据处理单元,在其中对它们进行频谱分析并产生它们的频谱以及匹配谐振曲线以及这些频谱的参数化以便确定每个第k个模式的观测参数组(pk、fk、γk和ek),为每个第k个模式作出参数化,并且随后所述观测参数组(pk、fk、γk和ek)与利用地球电离层腔中的ELF波传播模型的知识产生的存储在所述模型参数数据库中的基准参数组(pok、fok、γok和eok)相比较,然后,利用观测参数(pk、fk、γk和ek)和基准参数(pok、fok、γok和eok)之间偏差的最小化,选择其值尽可能接近所述观测参数(pk、fk、γk和ek)的基准参数(pok、fok、γok和eok),并且所选择的基准参数组(pok、fok、γok和eok)重构电火花放电的位置和强度并能够产生源头图。
6.根据权利要求1或2所述的方法,其特征在于通过下式所示匹配每个功率谱(Sx,Sy)和谐振曲线:
S ( ω ) = b + a ω α + Σ k = 1 K p k · [ 1 + e k · ( ω k - ω ) ] ( ω k - ω ) 2 + ( Γ k ) 2
其中:
S(ω)-匹配功率谱;
a-描述背景颜色噪声的参数;
b-描述背景宽带噪声的参数;
α-颜色噪声的频谱指数;
ω-比率(脉动);
pk-第k个谐振峰的最大功率值;
ek-第k个谐振峰的非对称参数;
ωk-第k个谐振蜂的谐振比率,其等于2πfk
Γk-第k个谐振峰的半宽度。
7.一种实时监测地球表面上风暴活动的装置,包括天线,记录观测信号的块,无线电发射系统和电子数据处理单元,其特征在于所述天线是两个有源磁感应天线(Hx,Hy),优选地彼此垂直,优选地沿着NS和EW方向,并且位于地表上或密切位于地表之下,它们通过屏蔽线连接到记录观测信号(Hx,Hy)的块,并且它们包含两个相同的具有放大器(Ax,Ay)、滤波器(BPFx,BPFy)、模数转换器(A/Cx,A/Cy)和控制系统(uP)的信号通道,所述控制系统(uP)的时钟是通过系统接收器与GPS的时间同步的,其中所述电子数据处理单元根据应用于ELF信号的频谱分析的算法和确定取决于风暴放电源的距离和强度的特征的算法而处理。
8.根据权利要求7所述的装置,其特征在于所述记录观测信号的块连接到发射在所述模数转换器中处理的数据的所述发射天线,而计算单元连接到与所述发射天线通讯的接收天线,所述发射天线连接到所述记录观测信号的块。
CN2008800194936A 2007-06-11 2008-06-11 用于实时监测地球表面风暴活动的方法和装置 Expired - Fee Related CN101802650B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
PLP382624 2007-06-11
PL382624A PL382624A1 (pl) 2007-06-11 2007-06-11 Sposób i urządzenie do monitorowania aktywności burzowej na powierzchni Ziemi w czasie rzeczywistym
PL385320A PL385320A3 (pl) 2008-05-30 2008-05-30 Sposób i urządzenie do monitorowania aktywności burzowej na powierzchni Ziemi w czasie rzeczywistym
PLP385320 2008-05-30
PCT/IB2008/052296 WO2008152587A2 (en) 2007-06-11 2008-06-11 Method and apparatus for monitoring storm activity on the earth's surface in real time

Publications (2)

Publication Number Publication Date
CN101802650A true CN101802650A (zh) 2010-08-11
CN101802650B CN101802650B (zh) 2013-06-05

Family

ID=39884341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008800194936A Expired - Fee Related CN101802650B (zh) 2007-06-11 2008-06-11 用于实时监测地球表面风暴活动的方法和装置

Country Status (12)

Country Link
US (1) US8332150B2 (zh)
EP (1) EP2165223B1 (zh)
JP (1) JP5562237B2 (zh)
KR (1) KR101357434B1 (zh)
CN (1) CN101802650B (zh)
BR (1) BRPI0813110B1 (zh)
CA (1) CA2690085C (zh)
ES (1) ES2400103T3 (zh)
MX (1) MX2009013584A (zh)
PL (1) PL2165223T3 (zh)
RU (1) RU2470332C2 (zh)
WO (1) WO2008152587A2 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107144883A (zh) * 2017-05-13 2017-09-08 中国地震局地质研究所 地震预测监测交变电磁场观测台网

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2559155C1 (ru) * 2014-01-30 2015-08-10 Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации Полевой индикатор естественного электромагнитного поля земли
KR101873347B1 (ko) 2017-04-26 2018-07-02 한국신발피혁연구원 폴리올레핀 난연필름 및 그 제조방법

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5153508A (en) * 1990-08-23 1992-10-06 United States Of America As Represented By The Administrator, National Aeronautics And Space Administration Method and apparatus for determining return stroke polarity of distant lightning
CN101833048A (zh) * 2010-03-23 2010-09-15 浙江大学 光电信号符合的闪电监测方法及其装置

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5305210A (en) * 1992-04-29 1994-04-19 B. F. Goodrich Flight Systems, Inc. Sampled data lightning strike detection and mapping system capable of early detection of an invalid strike from sampled data and quick resumption of monitoring an incoming signal
US5537318A (en) * 1994-07-13 1996-07-16 B. F. Goodrich Flightsystems, Inc. Lightning strike detection and mapping system
CA2159051A1 (en) * 1994-10-06 1996-04-07 John Howard Greenewald Thunderstorm cell detection and mapping system
US6125328A (en) * 1997-02-10 2000-09-26 Baron Services, Inc. System and method for projecting storms using NEXRAD attributes
JP3759831B2 (ja) * 1998-01-07 2006-03-29 株式会社サンコーシヤ ループアンテナ及びこれを用いた電磁波発生源位置標定システム
US6215294B1 (en) * 1998-12-30 2001-04-10 Ernest Wilson Coleman Storm warning method and apparatus
JP2000275358A (ja) * 1999-03-24 2000-10-06 Ikuo Watanabe 地球温暖化の計測方法、及び該方法に使用する装置
RU2188439C2 (ru) * 2000-10-04 2002-08-27 Центр международного сотрудничества по проблемам окружающей среды РАН Унифицированный генераторно-измерительный комплекс крайне низких и сверхнизких частот для геофизических исследований
JP2002258880A (ja) * 2001-03-06 2002-09-11 Sekiguchi Kikai Hanbai Kk 電波放出装置
RU2230336C2 (ru) * 2002-08-07 2004-06-10 Южно-Уральский государственный университет Однопунктовая система местоопределения гроз в ближней зоне
RU2253133C2 (ru) * 2003-07-31 2005-05-27 Государственное образовательное учреждение высшего профессионального образования "Южно-Уральский государственный университет" Способ определения местоположения молниевого разряда и многопунктовая система для его реализации
EP1711843B1 (de) * 2004-01-29 2014-05-28 nowcast GmbH System und verfahren zum erfassen, übermitteln und auswerten von durch elektromagnetische strahlung anfallenden daten und informationen
JP2005312755A (ja) * 2004-04-30 2005-11-10 Nippon Kenpan Seizo Kk 横臥式遠赤外線サウナ風呂
JP2006337329A (ja) * 2005-06-06 2006-12-14 Nippon Telegr & Teleph Corp <Ntt> 発雷位置推定システム及び方法
JP2007121127A (ja) * 2005-10-28 2007-05-17 Tokyo Electric Power Co Inc:The 雷放電位置標定システム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5153508A (en) * 1990-08-23 1992-10-06 United States Of America As Represented By The Administrator, National Aeronautics And Space Administration Method and apparatus for determining return stroke polarity of distant lightning
CN101833048A (zh) * 2010-03-23 2010-09-15 浙江大学 光电信号符合的闪电监测方法及其装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107144883A (zh) * 2017-05-13 2017-09-08 中国地震局地质研究所 地震预测监测交变电磁场观测台网

Also Published As

Publication number Publication date
ES2400103T3 (es) 2013-04-05
CA2690085C (en) 2014-10-28
US20100171485A1 (en) 2010-07-08
EP2165223A2 (en) 2010-03-24
EP2165223B1 (en) 2012-11-21
WO2008152587A2 (en) 2008-12-18
RU2009148713A (ru) 2011-07-20
WO2008152587A3 (en) 2009-02-05
RU2470332C2 (ru) 2012-12-20
KR101357434B1 (ko) 2014-02-03
JP5562237B2 (ja) 2014-07-30
PL2165223T3 (pl) 2013-04-30
BRPI0813110B1 (pt) 2019-02-05
CA2690085A1 (en) 2008-12-18
MX2009013584A (es) 2010-02-12
CN101802650B (zh) 2013-06-05
KR20100028579A (ko) 2010-03-12
JP2010529475A (ja) 2010-08-26
BRPI0813110A2 (pt) 2014-12-23
US8332150B2 (en) 2012-12-11

Similar Documents

Publication Publication Date Title
CN107085150B (zh) 一种短波发射天线三维立体方向图空中移动测量系统及方法
US9810804B2 (en) Collecting and transmitting control source electromagnetic signals
US20220128732A1 (en) System and method for determining characteristics of traveling ionospheric disturbances
CN104749640A (zh) 多源多频地空电磁探测发射方法
CN101802650B (zh) 用于实时监测地球表面风暴活动的方法和装置
Tsai et al. Three-dimensional numerical ray tracing on a phenomenological ionospheric model
CN102818941A (zh) 一种外场有扰环境下的电磁辐射发射测量方法
Belehaki et al. Ionospheric specification and forecasting based on observations from European ionosondes participating in DIAS project
CN1144065C (zh) 气象电磁测量系统
US9568632B2 (en) Practical electromagnetic gradiometer
Saraev et al. Foot, mobile and controlled source modifications of the radiomagnetotelluric method
Gurses et al. Ultra-sensitive broadband “AWESOME” electric field receiver for nanovolt low-frequency signals
CN116008722A (zh) 基于vmd变分模态分解的接触网故障测距系统的计算方法
Liatos et al. Characterization of 100-kHz noise in the lightning current derivative signals measured at the CN tower
Ganguly et al. New generation topside sounder
Yearby Magnetospheric VLF line radiation.
CN117074758A (zh) 局域电网地磁感应电流的风险评估方法
Pożoga et al. Ionospheric Scintillation Diagnostic on LOFAR ILT Network in Single Station Mode
CN113992287A (zh) 一种低频时码授时信号载波频率优选方法及系统
JPH11352207A (ja) 複数衛星による電磁波発生源観測方法
RU2379716C2 (ru) Система для определения предвестников землетрясения
Nickolaenko Impact of the ionosphere day-night non-uniformity on amplitude of ELF radio signal
CN118011479A (zh) 基于地震方法的道路地质扫描系统及数据处理成像方法
Lad Characterization of Atmospheric Noise and Precipitation Static in the Long Range Navigation (Loran-C) Band for Aircraft
Nanjo et al. Statistical relation between lightning and anomalous propagation in the line-of-sight VHF band

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130605

Termination date: 20160611

CF01 Termination of patent right due to non-payment of annual fee