CN101799560A - 多个声音传感器雷电定位仪及其定位方法 - Google Patents
多个声音传感器雷电定位仪及其定位方法 Download PDFInfo
- Publication number
- CN101799560A CN101799560A CN200910310395A CN200910310395A CN101799560A CN 101799560 A CN101799560 A CN 101799560A CN 200910310395 A CN200910310395 A CN 200910310395A CN 200910310395 A CN200910310395 A CN 200910310395A CN 101799560 A CN101799560 A CN 101799560A
- Authority
- CN
- China
- Prior art keywords
- thunder
- lightning
- sound
- data
- coordinate
- 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
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种多个声音传感器雷电定位仪及其定位方法,其定位仪包括在四个观测点处布置的四个传感器、一个四通道数据采集卡和数据处理***。在雷电发生时,四个空间观测点的四个声音传感器分别记录各观测点雷声到达的时刻,四组时间数据经四通道数据采集卡A/D转换后,被传送到数据处理***,数据处理***通过逻辑运算,求取三个空间观测点雷声到达的时间差,根据声音在空气中传播的速度,计算出三个空间观测点的距离差;根据空间三个双曲面相交可以确定一个点的原理,根据方程组确定雷电产生位置的坐标,本发明在保证测量高精度的同时,简化了定位测量方法,降低了硬件要求和设备成本,实现了对雷电的准确自动定位,可在气象探测领域广泛应用。
Description
技术领域
本发明涉及气象探测技术领域,特别是一种雷电的定位方法。
背景技术
目前被国际上普遍认可的雷电定位方法大致可分为五种。第一种是改进了门控磁脉冲定向法,如广泛应用的地闪定向仪(DF),它能识别典型的对地闪电的LF(低频)频段辐射电磁场波形,测定放电波形的峰值和方位。多站DF布网则可确定闪电发生的位置。第二种是工作在LF频段的长基线TOA(时间到达法)技术,如闪电定位及跟踪***(LPATS)。第三种是工作在VHF(甚高频)频段的干涉仪方法,其特点是可以同时探测云闪和地闪,并可了解放电的大致过程,但探测距离较短。第四种是工作在VHF频段的TOA技术,如LDAR***。第五种是利用雷声差探测闪电通道的技术。从性能看,前两种定位技术可以在较大的范围内布网,后三种比较适用于特殊用途及闪电物理的研究。目前闪电定位仪的定位原理主要是根据电磁场的方法进行定位。
发明内容
本发明的目的是提供一种多个声音传感器雷电定位仪及其定位方法,要解决现有的雷电测量方法原理复杂,设备投资大,功能单一的技术问题。
为实现所述目的,本发明采用如下技术方案:
技术方案一:一种多个声音传感器雷电定位方法,其特征在于步骤如下:
步骤1,未知的雷电坐标为P点,在四个空间观测点O、A、B、C各布置一个声音传感器,四个声音传感器与一个四通道数据采集卡的四个信号输入通道连接,四通道数据采集卡的输出端与数据处理***的信号输入接口连接;
四个空间观测点的坐标分别为O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO;
步骤2,四个空间观测点的四个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻to、ta、tb、tc,四组时间数据经四通道数据采集卡A/D转换后,被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,将每个通道的雷声到达时刻,通过逻辑运算,求取A、B、C三个空间观测点分别相对于O点的雷声到达时间差toa、tob、toc;
步骤4,数据处理***根据声音在空气中传播的速度v,距离D=时间t×速度v的公式,计算出A、B、C三个空间观测点分别相对于P点和O点的距离差为:Doa、Dob、Doc;
步骤5,数据处理***根据空间三个双曲面相交可以确定一个点的原理,则四个不在同一平面的声音传感器就可以确定雷电产生位置P的坐标,根据方程组:
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xa)2+(y-ya)2+(z-za)2]1/2=Doa
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xb)2+(y-yb)2+(z-zb)2]1/2=Dob
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xc)2+(y-yc)2+(z-zc)2]1/2=Doc
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
所述声音传感器可以是麦克风或录音机,采样速度≥30微秒。
技术方案二:一种多个声音传感器雷电定位方法,其特征在于步骤如下:
步骤1,未知的雷电坐标为P点,在四个观测点O、A、B、C布置三个声音传感器和一个光传感器,三个声音传感器和一个光传感器与一个多通道数据采集卡的三个信号输入通道连接,多通道数据采集卡的输出端与数据处理***的信号输入接口连接;
四个观测点的坐标分别为,O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO;
步骤2,三个空间观测点的三个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻,三组时间数据to、ta、tb经三通道数据采集卡A/D转换后,被传送到数据处理***,广角摄像头记录闪电发生时刻,也被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,根据光和声音传播速度之差可以求出声音的相对距离再根据声音到达地面三点的时间差定位闪电的方位角和仰角,即利用摄像头拍摄闪电发生时刻,再获取到达地面的雷声,求取闪电到地面得距离。假设闪电坐标为P(x,y,z),地面三点分别为A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),闪电到三点的距离分别为:Da,Db,Dc。则方程组为:
(x-xa)2+(y-ya)2+(z-za)2=Da 2
(x-xb)2+(y-yb)2+(z-zb)2=Db 2
(x-xc)2+(y-yc)2+(z-zc)2=Dc 2
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
所述观测点A为坐标原点,坐标分别为:B(10,0,0),C(0,10,0)。由上三式可以解出点P1(x1,y1,z1),P2(x2,y2,z2)的坐标,令Da=A,Db=B,Dc=C则:
x1=x2=5+0.5A2+0.5B2
y1=y2=5+0.5A2-0.5C2
z1=0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)
z2=-0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)。
所述三个声音传感器位于同一海拔高度,并且三地连线可以构成一个直角三角形,光传感器位于直角三角形中任意一地。
所述声音传感器可以是麦克风或录音机,采样速度≥30微秒,光传感器是广角摄像头。
步骤6后利用公式求取雷电发生位置P的仰角。
一种多个声音传感器雷电定位仪,包括在四个观测点处布置的四个传感器,其特征在于:所述传感器是四个声音传感器或者三个声音传感器和一个光传感器的组合,四个传感器与一个四通道数据采集卡的四个信号输入端子连接,四通道数据采集卡采集每个传感器声光到达的时刻,其四个信号输出端子与数据处理***中的中央处理单元输入接口连接,中央处理单元内包括有控制器、运算器和寄存器,中央处理单元与存储器、显示器和输入设备连接,存储器内存入操作***、数据处理程序、计算公式和观测数据,计算公式包括求取相应两观测点声光到达时刻的时间差公式、根据声音在空气中传播的速度计算出雷电到达两地的距离差的公式、求得雷电发生位置的坐标、方位角和仰角的公式。
与现有技术相比本发明具有以下特点和有益效果:
本发明是利用闪电的声光特性进行定位,根据闪电光谱分析在可见光及近红外光谱范围内辐射能量最为集中,闪电的声音信号一般集中在频率低于100HZ的可闻区,和次声区,声强的峰值集中在40HZ左右。通过测量闪电的声光发生时刻利用其速度差可以定位。
本发明的方案一不需要摄像头只需四个麦克风即可定位,定位方法简单,对硬件要求较高,要求数据采集卡的响应速度快,在相隔10米情况下,要求能够精确到1厘米,即采样速度要能达到1/34000秒=30微秒。要求麦克风的一致性,稳定性好,但对其灵敏度没有要求。
方案二,在方案一基础上加上摄像头,降低对硬件的要求,提高定位精度,但定位方法要求有闪电识别和闪电同步。
本发明在保证测量高精度的同时,简化了定位测量方法,降低了硬件的要求,降低了设备成本,实现了对雷电的准确自动定位,可在气象探测领域广泛应用。
附图说明
下面结合附图对本发明做进一步详细的说明。
图1是实施例一的传感器布置示意图。
图2是实施例一的定位原理示意图。
图3是实施例一的设备结构框图。
图4是实施例一的定位观测流程图。
图5是实施例二的传感器布置示意图。
图6是实施例二的定位观测流程图。
附图标记:1-声音传感器、2-光传感器。
具体实施方式
一种多个声音传感器雷电定位仪,包括在四个观测点处布置的四个传感器,参见图1,所述传感器是声音传感器,声音传感器是四个麦克风。参见图5,所述传感器是三个声音传感器和一个光传感器的组合,例如:三个麦克风和一个广角摄像头的组合。参见图3,四个传感器与一个四通道数据采集卡的四个信号输入端子连接,四通道数据采集卡采集每个传感器声光到达的时刻,其四个信号输出端子与数据处理***中的中央处理单元输入接口连接,中央处理单元内包括有控制器、运算器和寄存器,中央处理单元与存储器、显示器和输入设备连接,存储器内存入操作***、数据处理程序、计算公式和观测数据,计算公式包括求取相应两观测点声光到达时刻的时间差公式、根据声音在空气中传播的速度计算出雷电到达两地的距离差的公式、求得雷电发生位置的坐标、方位角和仰角的公式。
四通道数据采集卡采集每个传感器声光到达的时刻,其四个信号输出端子与数据处理***中的中央处理单元输入接口连接,中央处理单元内有控制器、计算器、寄存器、计数器、缓冲器和内部总线,并与外部的显示器、输入设备和存储器连接,存储器内存入操作***、数据处理程序、计算公式和观测数据,计算公式包括求取相应两观测点声光到达时刻的时间差公式、根据声音在空气中传播的速度计算出雷电到达两地的距离差的公式、求得雷电发生位置的坐标、方位角和仰角的公式。
实施例一参见图1-图4所示,这种多个声音传感器雷电定位方法,其特征在于:
步骤1,未知的雷电坐标为P点,在四个空间观测点O、A、B、C各布置一个声音传感器,声音传感器是麦克风或录音机,采样速度≥30微秒,声音传感器1是麦克风或录音机,采样速度≥30微秒,光传感器2是广角摄像头,四个声音传感器与一个四通道数据采集卡的四个信号输入通道连接,四通道数据采集卡的输出端与数据处理***的信号输入接口连接;
四个空间观测点的坐标分别O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO;
步骤2,四个空间观测点的四个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻to、ta、tb、tc,四组时间数据经四通道数据采集卡A/D转换后,被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,将每个通道的雷声到达时刻,通过逻辑运算,求取A、B、C三个空间观测点分别相对于O点的雷声到达时间差toa、tob、toc;
步骤4,数据处理***根据声音在空气中传播的速度v,距离D=时间t×速度v的公式,计算出A、B、C三个空间观测点分别相对于P点和O点的距离差为:Doa、Dob、Doc;
步骤5,数据处理***根据空间三个双曲面相交可以确定一个点的原理,则四个不在同一平面的声音传感器就可以确定雷电产生位置P的坐标,根据方程组:
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xa)2+(y-ya)2+(z-za)2]1/2=Doa
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xb)2+(y-yb)2+(z-zb)2]1/2=Dob
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xc)2+(y-yc)2+(z-zc)2]1/2=Doc
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
利用公式求取雷电发生位置P的方位角θ。方位角:从观测点的指北方向线起,依顺时针方向到雷电P目标方向线之间的水平夹角,叫方位角。方位角的取值范围为0-360度。
实施例一的原理如下:根据几何学原理已知空间任一点到相对俩定点的距离之差为常数的轨迹为双曲面方程,空间俩双曲面相交为一曲线,空间三双曲面相交可以确定空间一点。为了能构成三双曲面需要四个空间点。
假设空间任一点坐标为P(x,y,z),相对的四个空间定点坐标分别为:O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc)。相对距离差为:Doa,Dob,Doc。
其中Doc一定大于零。建立方程组为:
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xa)2+(y-ya)2+(z-za)2]1/2=Doa
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xb)2+(y-yb)2+(z-zb)2]1/2=Dob
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xc)2+(y-yc)2+(z-zc)2]1/2=Doc
为了计算和布局方便使点O为坐标原点A,B,C坐标分别为:A(0,10,0),B(10,0,0),C(0,0,2)。由上三式可以解出点P1(x1,y1,z1),P2(x2,y2,z2)的坐标。令Doa=A,Dob=B,Doc=C,则:
x1=-0.5(25C2A2+100CA-100A2+10000-B3A+100AB+B2A2-25C3A-100B2-2500C2-(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C2+400A4B2+100B4A2-25C2A5+50A5C3-25A4C4-B4A4+20000BA3-200B3A3+2B3A5-A5B2-200A5B+100A5+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)1/2))(-100+B2+A2+25C2)
x2=-0.5(25C2A2+100CA-100A2+10000-B3A+100AB+B2A2-25C3A-100B2-2500C2+(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C2+400A4B2+100B4A2-25C2A6+50A5C3-25A4C4-B4A4+20000BA3-200B3A3+2B3A5-A6B2-200A5B+100A6+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)1/2))/(-100+B2+A2+25C2)
y1=0.5*(100A-B*(25C2A2+100CA-100A2+10000-B3A+100AB+B2A2-25C3A-100B2-2500C2-(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C2+400A4B2+100B4A2-25C2A6+50A5C2-25A4C4-B4A4+20000BA3-200B3A3+2B3A5-A5B2-200A5B+100A5+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)1/2)/(-100+B2+A2+25C2)-100B+BA2-B2A)/A
y2=0.5*(100A-B*(25C2A2+100CA-100A2+10000-B3A+100AB+B2A2-25C3A-100B2-2500C2+(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C2+400A4B2+100B4A2-25C2A5+50A5C3-25A4C4-B4A4+20000BA3-200B3A3+2B3A5-A5B2-200A5B+100A5+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)1/2)/(-100+B2+A2+25C2)-100B+BA2-B2A)/A
z1=0.25(4A-C(25C2A2+100CA-100A2+10000-B3A+100AB+B2A2-25C3A-100B2-2500C2-(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C4+400A4B2+100B4A2-25C2A4+50A5C3-25A4C4-B4A4+20000BA3-200B3A3+2B3A3-A4B3-200A5B+100A3+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)(1/2))/(-100+B2+A2+25C2)-100C-AC4+A2C)/A
z2=0.25(4A-C(25C2A2+100CA-100A2+10000-B3A+100AB+B2A4-25C5A-100B2-2500C2+(5200C2B2A2-5000C3BA2-25C4B2A2+20000CBA2-200CB3A2-30400A4-520000C2A2+20000CA3+2040000A2-25B4A2C2+400A4B2+100B4A2-25C2A5+50A5C3-25A4C4-B4A4+20000BA3-200B3A3+2B3A5-A5B2-200A5B+100A5+50C3B3A2-30400B2A2+5200C2A4+2500C4A2-200CA5-5000C3A3)(1/2))/(-100+B2+A2+25C2)-100C-AC2+A2C)/A
实施例二参见图5、图6所示,这种多个声音传感器雷电定位方法,
步骤1,未知的雷电坐标为P点,在四个观测点O、A、B、C布置三个声音传感器和一个光传感器,三个声音传感器位于同一海拔高度,并且三地连线可以构成一个直角三角形,光传感器位于直角三角形中任意一地,三个声音传感器和一个光传感器与一个多通道数据采集卡的三个信号输入通道连接,多通道数据采集卡的输出端与数据处理***的信号输入接口连接;
四个观测点的坐标分别为O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO。
步骤2,三个空间观测点的三个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻,三组时间数据to、ta、tb经三通道数据采集卡A/D转换后,被传送到数据处理***;广角摄像头记录闪电发生时刻,也被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,根据光和声音传播速度之差可以求出声音的相对距离再根据声音到达地面三点的时间差定位闪电的方位角和仰角,即利用摄像头拍摄闪电发生时刻,再获取到达地面的雷声,求取闪电到地面得距离。假设闪电坐标为P(x,y,z),地面三点分别为,闪电到三点的距离分别为:A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),闪电到三点的距离分别为:Da,Db,Dc。则方程组为:
(x-xa)2+(y-ya)2+(z-za)2=Da 2
(x-xb)2+(y-yb)2+(z-zb)2=Db 2
(x-xc)2+(y-yc)2+(z-zc)2=Dc 2
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
所述观测点A为坐标原点,坐标分别为:B(10,0,0),C(0,10,0)。由上三式可以解出点P1(x1,y1,z1),P2(x2,y2,z2)的坐标,令Da=A,Db=B,Dc=C则:
x1=x2=5+0.5A2+0.5B2
y1=y2=5+0.5A2-0.5C2
z1=0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)
z2=-0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)。
实施例二的原理如下:
根据光和声音传播速度之差可以求出声音的相对距离再根据声音到达地面三点的时间差定位闪电的方位角和仰角,即利用摄像头拍摄闪电发生时刻,再获取到达地面的雷声,求取闪电到地面得距离。假设闪电坐标为P(x,y,z),地面三点分别为A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),闪电到三点的距离分别为:Da,Db,Dc。则方程组为:
(x-xa)2+(y-ya)2+(z-za)2=Da 2
(x-xb)2+(y-yb)2+(z-zb)2=Db 2
(x-xc)2+(y-yc)2+(z-zc)2=Dc 2
为了计算和布局方便使点为A坐标原点坐标分别为:B(10,0,0),C(0,10,0)。由上三式可以解出点P1(x1,y1,z1),P2(x2,y2,z2)的坐标。令Da=A,Db=B,Dc=C则:
x1=x2=5+0.5A2+0.5B2
y1=y2=5+0.5A2-0.5C2
z1=0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)
z2=-0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)(舍)
Claims (10)
1.一种多个声音传感器雷电定位方法,其特征在于步骤如下:
步骤1,未知的雷电坐标为P点,在四个空间观测点O、A、B、C各布置一个声音传感器,四个声音传感器与一个四通道数据采集卡的四个信号输入通道连接,四通道数据采集卡的输出端与数据处理***的信号输入接口连接;
四个空间观测点的坐标分别为,O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO;
步骤2,四个空间观测点的四个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻totatbtc,四组时间数据经四通道数据采集卡A/D转换后,被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,将每个通道的雷声到达时刻,通过逻辑运算,求取A、B、C三个空间观测点分别相对于O点的雷声到达时间差toatobtoc;
步骤4,数据处理***根据声音在空气中传播的速度v,距离D=时间t×速度v的公式,计算出A、B、C三个空间观测点分别相对于P点和O点的距离差为:DoaDobDoc;
步骤5,数据处理***根据空间三个双曲面相交可以确定一个点的原理,则四个不在同一平面的声音传感器就可以确定雷电产生位置P的坐标,根据方程组:
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xa)2+(y-ya)2+(z-za)2]1/2=Doa
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xb)2+(y-yb)2+(z-zb)2]1/2=Dob
[(x-xo)2+(y-yo)2+(z-zo)2]1/2-[(x-xc)2+(y-yc)2+(z-zc)2]1/2=Dx
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
4.根据权利要求1所述多个声音传感器雷电定位方法,其特征在于:所述声音传感器是麦克风或录音机,采样速度≥30微秒。
5.一种多个声音传感器雷电定位方法,其特征在于步骤如下:
步骤1,未知的雷电坐标为P点,在四个观测点O、A、B、C布置三个声音传感器和一个光传感器,三个声音传感器和一个光传感器与一个多通道数据采集卡的三个信号输入通道连接,多通道数据采集卡的输出端与数据处理***的信号输入接口连接;四个观测点的坐标分别为O(xo,yo,zo),A(xa,ya,za),B(xb,yb,zb),C(xc,yc,zc),则A、B、C三个空间观测点分别相对于O点的距离为定值:AO、BO、CO;
步骤2,三个空间观测点的三个声音传感器在雷电发生时,分别记录各观测点雷声到达的时刻,三组时间数据totatb经三通道数据采集卡A/D转换后,被传送到数据处理***;广角摄像头记录闪电发生时刻,也被传送到数据处理***;
步骤3,数据处理***内存储有操作***、数据处理程序、计算公式和观测数据,根据光和声音传播速度之差可以求出声音的相对距离再根据声音到达地面三点的时间差定位闪电的方位角和仰角,即利用摄像头拍摄闪电发生时刻,再获取到达地面的雷声,求取闪电到地面得距离。假设闪电坐标为P(x,y,z),地面三点分别为A(xa,ya,za),B(xb,yb,zb),C (xc,yc,zc),闪电到三点的距离分别为:Doa,DobDoc。则方程组为:
(x-xa)2+(y-ya)2+(z-za)2=Da 2
(x-xb)2+(y-yb)2+(z-zb)2=Db 2
(x-xc)2+(y-yc)2+(z-zc)2=Dc 2
最终解出雷电发生位置P的坐标:P(x,y,z),通过显示设备显示数据处理结果。
6.根据权利要求5所述多个声音传感器雷电定位方法,其特征在于:所述观测点A为坐标原点,坐标分别为:B(10,0,0),C(0,10,0)。由上三式可以解出点P1(x1,y1,z1),P1(x2,y2,z2)的坐标,令Da=A,Db=B,Dc=C,则:
x1=x2=5+0.5A2+0.5B2
y1=y2=5+0.5A2-0.5C2
z1=0.5(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)
z2=-05(-20000+200B2-2A4+2A2B2-B4+200C2+2A2C2-C4)(1/2)。
7.根据权利要求5所述多个声音传感器雷电定位方法,其特征在于:所述三个声音传感器位于同一海拔高度,并且三地连线可以构成一个直角三角形,光传感器位于直角三角形中任意一地。
8.根据权利要求5所述多个声音传感器雷电定位方法,其特征在于:所述声音传感器是麦克风或录音机,采样速度≥30微秒,光传感器是广角摄像头。
10.一种根据权利要求1-9所述方法的多个声音传感器雷电定位仪,包括在四个观测点处布置的四个传感器,其特征在于:所述传感器是四个声音传感器或者三个声音传感器和一个光传感器的组合,四个传感器与一个四通道数据采集卡的四个信号输入端子连接,四通道数据采集卡采集每个传感器声光到达的时刻,其四个信号输出端子与数据处理***中的中央处理单元输入接口连接,中央处理单元内包括有控制器、运算器和寄存器,中央处理单元与存储器、显示器和输入设备连接,存储器内存入操作***、数据处理程序、计算公式和观测数据,计算公式包括求取相应两观测点声光到达时刻的时间差公式、根据声音在空气中传播的速度计算出雷电到达两地的距离差的公式、求得雷电发生位置的坐标、方位角和仰角的公式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910310395A CN101799560A (zh) | 2009-11-25 | 2009-11-25 | 多个声音传感器雷电定位仪及其定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910310395A CN101799560A (zh) | 2009-11-25 | 2009-11-25 | 多个声音传感器雷电定位仪及其定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101799560A true CN101799560A (zh) | 2010-08-11 |
Family
ID=42595313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910310395A Pending CN101799560A (zh) | 2009-11-25 | 2009-11-25 | 多个声音传感器雷电定位仪及其定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101799560A (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102413563A (zh) * | 2011-11-09 | 2012-04-11 | 中国科学院深圳先进技术研究院 | 无线信号源定位方法及*** |
CN102680939A (zh) * | 2012-04-06 | 2012-09-19 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种基于声光电同步观测的雷电定位方法及装置 |
CN103209030A (zh) * | 2012-01-12 | 2013-07-17 | 宏碁股份有限公司 | 电子装置及其数据传输方法 |
CN104616687A (zh) * | 2015-01-14 | 2015-05-13 | 小米科技有限责任公司 | 录音的方法及装置 |
CN104614709A (zh) * | 2015-01-19 | 2015-05-13 | 成都信息工程学院 | 一种基于声学和电磁学的雷声定位***及方法 |
CN105068150A (zh) * | 2015-07-24 | 2015-11-18 | 中国石油化工股份有限公司 | 一种基于多数据运算的雷云定位方法 |
CN105242790A (zh) * | 2015-11-04 | 2016-01-13 | 钱雄 | 一种智能设备控制方法及控制装置 |
CN105954653A (zh) * | 2016-04-26 | 2016-09-21 | 西安交通大学 | 一种gis现场耐压试验击穿点定位方法及*** |
CN106405253A (zh) * | 2016-08-24 | 2017-02-15 | 中国气象科学研究院 | 一种目标闪电辐射源定位方法及装置 |
CN107884675A (zh) * | 2016-09-29 | 2018-04-06 | 山东信通电子股份有限公司 | 具有输电线路雷击故障定位的监测装置 |
CN109001678A (zh) * | 2017-06-06 | 2018-12-14 | 中国科学技术大学 | 一种基于三维麦克风阵列的雷声检测与定位方法 |
CN109035709A (zh) * | 2018-09-13 | 2018-12-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种山洪泥石流的预警方法、装置及*** |
CN109599852A (zh) * | 2018-12-05 | 2019-04-09 | 北京雷布斯雷电安全科技有限公司 | 防雷*** |
CN109799482A (zh) * | 2017-11-16 | 2019-05-24 | 航天科工惯性技术有限公司 | 基于次声波三角定位的泥石流监测***及方法 |
CN110850509A (zh) * | 2019-11-13 | 2020-02-28 | 山东省地震局 | 一种可定向的综合雷电预警设备 |
CN111624412A (zh) * | 2020-04-21 | 2020-09-04 | 北京信息科技大学 | 一种闪电连接点定位方法、***、设备及可读存储介质 |
CN112684487A (zh) * | 2020-12-22 | 2021-04-20 | 西安交通大学 | 一种三维空间环境内快速准确定位放射源的定位方法 |
CN113421457A (zh) * | 2021-06-08 | 2021-09-21 | 北部湾大学 | 一种基于手机的雷声检测和定位方法 |
CN114325569A (zh) * | 2022-03-02 | 2022-04-12 | 广州国睿电子科技有限公司 | 一种雷电探测装置 |
-
2009
- 2009-11-25 CN CN200910310395A patent/CN101799560A/zh active Pending
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102413563A (zh) * | 2011-11-09 | 2012-04-11 | 中国科学院深圳先进技术研究院 | 无线信号源定位方法及*** |
CN102413563B (zh) * | 2011-11-09 | 2014-01-08 | 中国科学院深圳先进技术研究院 | 无线信号源定位方法及*** |
CN103209030A (zh) * | 2012-01-12 | 2013-07-17 | 宏碁股份有限公司 | 电子装置及其数据传输方法 |
CN103209030B (zh) * | 2012-01-12 | 2015-05-13 | 宏碁股份有限公司 | 电子装置及其数据传输方法 |
CN102680939A (zh) * | 2012-04-06 | 2012-09-19 | 国网电力科学研究院武汉南瑞有限责任公司 | 一种基于声光电同步观测的雷电定位方法及装置 |
CN104616687A (zh) * | 2015-01-14 | 2015-05-13 | 小米科技有限责任公司 | 录音的方法及装置 |
CN104614709A (zh) * | 2015-01-19 | 2015-05-13 | 成都信息工程学院 | 一种基于声学和电磁学的雷声定位***及方法 |
CN105068150A (zh) * | 2015-07-24 | 2015-11-18 | 中国石油化工股份有限公司 | 一种基于多数据运算的雷云定位方法 |
CN105242790A (zh) * | 2015-11-04 | 2016-01-13 | 钱雄 | 一种智能设备控制方法及控制装置 |
CN105242790B (zh) * | 2015-11-04 | 2019-06-21 | 钱雄 | 一种智能设备控制方法及控制装置 |
CN105954653A (zh) * | 2016-04-26 | 2016-09-21 | 西安交通大学 | 一种gis现场耐压试验击穿点定位方法及*** |
CN105954653B (zh) * | 2016-04-26 | 2020-03-31 | 西安交通大学 | 一种gis现场耐压试验击穿点定位方法及*** |
CN106405253A (zh) * | 2016-08-24 | 2017-02-15 | 中国气象科学研究院 | 一种目标闪电辐射源定位方法及装置 |
CN107884675A (zh) * | 2016-09-29 | 2018-04-06 | 山东信通电子股份有限公司 | 具有输电线路雷击故障定位的监测装置 |
CN109001678A (zh) * | 2017-06-06 | 2018-12-14 | 中国科学技术大学 | 一种基于三维麦克风阵列的雷声检测与定位方法 |
CN109001678B (zh) * | 2017-06-06 | 2021-10-01 | 中国科学技术大学 | 一种基于三维麦克风阵列的雷声检测与定位方法 |
CN109799482A (zh) * | 2017-11-16 | 2019-05-24 | 航天科工惯性技术有限公司 | 基于次声波三角定位的泥石流监测***及方法 |
CN109035709A (zh) * | 2018-09-13 | 2018-12-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种山洪泥石流的预警方法、装置及*** |
CN109599852A (zh) * | 2018-12-05 | 2019-04-09 | 北京雷布斯雷电安全科技有限公司 | 防雷*** |
CN109599852B (zh) * | 2018-12-05 | 2019-12-17 | 北京雷布斯雷电安全科技有限公司 | 防雷*** |
CN110850509A (zh) * | 2019-11-13 | 2020-02-28 | 山东省地震局 | 一种可定向的综合雷电预警设备 |
CN111624412A (zh) * | 2020-04-21 | 2020-09-04 | 北京信息科技大学 | 一种闪电连接点定位方法、***、设备及可读存储介质 |
CN111624412B (zh) * | 2020-04-21 | 2022-04-08 | 北京信息科技大学 | 一种闪电连接点定位方法、***、设备及可读存储介质 |
CN112684487A (zh) * | 2020-12-22 | 2021-04-20 | 西安交通大学 | 一种三维空间环境内快速准确定位放射源的定位方法 |
CN113421457A (zh) * | 2021-06-08 | 2021-09-21 | 北部湾大学 | 一种基于手机的雷声检测和定位方法 |
CN114325569A (zh) * | 2022-03-02 | 2022-04-12 | 广州国睿电子科技有限公司 | 一种雷电探测装置 |
CN114325569B (zh) * | 2022-03-02 | 2022-06-17 | 广州国睿电子科技有限公司 | 一种雷电探测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101799560A (zh) | 多个声音传感器雷电定位仪及其定位方法 | |
CN101957442B (zh) | 一种声源定位装置 | |
CN104898091B (zh) | 基于迭代优化算法的麦克风阵列自校准声源定位*** | |
CN108614268B (zh) | 低空高速飞行目标的声学跟踪方法 | |
KR20180063263A (ko) | 3차원 공간 검출 시스템, 측위 방법 및 시스템 | |
CN103454619A (zh) | 一种星载微波跟瞄雷达的电轴光学标定***及其标定方法 | |
CN102413563B (zh) | 无线信号源定位方法及*** | |
CN103941223A (zh) | 声源定位***及其方法 | |
CN104035065A (zh) | 基于主动旋转的声源定向装置及其应用方法 | |
US20200138408A1 (en) | Handheld three-dimensional ultrasound imaging system and method | |
CN103134387B (zh) | 一种低空慢速小目标探测与拦截***标定方法 | |
CN102680939A (zh) | 一种基于声光电同步观测的雷电定位方法及装置 | |
CN107861096A (zh) | 基于声音信号到达时间差的最小二乘测向方法 | |
EP0634881A1 (en) | Determination of position | |
CN110779616A (zh) | 一种飞行试验声爆信号测量方法 | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
CN107884743A (zh) | 适用于任意结构声音阵列的波达方向智能估计方法 | |
CN104655134A (zh) | 一种基于gps时标的多传感器数据采集*** | |
CN109061723B (zh) | 一种隧洞岩爆孕育过程的微震源高精度定位方法及*** | |
CN101035386A (zh) | 一种实现移动终端精确定位的装置及方法 | |
CN112162239B (zh) | 一种基于水平门型阵列的弹着点定位方法 | |
CN206638818U (zh) | 一种大深度长距离高精度水下超短基线定位控制*** | |
CN2667505Y (zh) | 双向激光测宽及定中仪 | |
CN111624552B (zh) | 一种基于声波渡越时间测量的地下管道定位***及方法 | |
CN206209409U (zh) | 无人机控制***和无人机 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20100811 |