JP2012189387A - Lightning discharge position orientation system - Google Patents

Lightning discharge position orientation system Download PDF

Info

Publication number
JP2012189387A
JP2012189387A JP2011051799A JP2011051799A JP2012189387A JP 2012189387 A JP2012189387 A JP 2012189387A JP 2011051799 A JP2011051799 A JP 2011051799A JP 2011051799 A JP2011051799 A JP 2011051799A JP 2012189387 A JP2012189387 A JP 2012189387A
Authority
JP
Japan
Prior art keywords
lightning discharge
peak
time
electric field
differential signal
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.)
Withdrawn
Application number
JP2011051799A
Other languages
Japanese (ja)
Inventor
Narimitsu Okabe
成光 岡部
Sakae Taniguchi
栄 谷口
Nobumi Nagaraki
宣巳 長良木
Masaru Ishii
勝 石井
Motohisa Saito
幹久 齋藤
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.)
University of Tokyo NUC
Tokyo Electric Power Company Holdings Inc
Original Assignee
Tokyo Electric Power Co Inc
University of Tokyo NUC
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 Tokyo Electric Power Co Inc, University of Tokyo NUC filed Critical Tokyo Electric Power Co Inc
Priority to JP2011051799A priority Critical patent/JP2012189387A/en
Publication of JP2012189387A publication Critical patent/JP2012189387A/en
Withdrawn legal-status Critical Current

Links

Images

Abstract

PROBLEM TO BE SOLVED: To more improve position orientation accuracy in a lightning discharge position orientation system using an arrival time difference approach.SOLUTION: As a representative configuration of a lightning discharge position orientation system, there are provided an antenna 104 which receives electromagnetic waves including electromagnetic waves with lightning discharge and outputs a differentiation signal of an electric field, a plurality of receiving stations 102 including an integration circuit 114 which outputs an electric field signal and a processing section 112 which specifies an arrival time of electromagnetic waves with lightning discharge, and a central processing unit 134 which orients a position where lightning discharge occurs from the arrival time of the electromagnetic waves with lightning discharge in the plurality of receiving stations. The processing section 112 determines a zero level by averaging the differentiation signals for a predetermined time width sufficiently retroactively from the time of a peak of the electric field signals, determines a peak of differentiation signals at a time preceding to the peak of the electric field signals, and defines a time when the value of the differentiation signal becomes a predetermined value of 40 to 90% from the zero level to the peak of the differentiation signal at the time preceding to the peak of the differentiation signal, as an arrival time of electromagnetic waves with lightning discharge.

Description

本発明は、雷放電に伴う電磁波の到達時刻を用いて雷放電が発生した位置を高精度に標定する雷放電位置標定システムに関する。   The present invention relates to a lightning discharge position locating system for locating a position where a lightning discharge has occurred with high accuracy using the arrival time of an electromagnetic wave accompanying the lightning discharge.

落雷発生時には、被害を調べ、復旧作業を行うため、速やかに落雷位置を知る必要がある。そこで従来から、落雷位置の標定が行われている。現在用いられている落雷位置標定方法は交会法および到達時間差法と呼ばれる方法を組み合わせた方法で行われている。交会法は2箇所以上の受信局で雷放電に伴う電磁波の到来方位を検出し、その交点から落雷点の位置標定を行う。このため、雷放電に伴う電界変化波形が複雑、または複数の箇所でほぼ同時に落雷が発生しても位置標定を行う事が出来る。到達時間差法は、落雷位置標定(2次元の位置標定)では3ヶ所以上の受信局で、雷放電に伴う電磁波の各観測局への到達時刻を記録し、到達時間差から落雷点の位置標定を行う方法である。各局への到達時間差はGPS時計を用いて高精度で記録するため、交会法に比べて精度良く落雷位置の標定ができる。位置標定精度は到達時間差法の方が優れている為、現行のシステムでは到達時間差法で位置標定が難しい場合にのみ交会法で位置標定を行う。   When a lightning strike occurs, it is necessary to promptly know the location of the lightning strike in order to investigate damage and perform recovery work. Therefore, the location of lightning strikes has been conventionally performed. The lightning strike location method currently used is performed by a combination of the so-called meeting method and arrival time difference method. In the meeting method, the arrival direction of electromagnetic waves accompanying lightning discharge is detected at two or more receiving stations, and the location of the lightning strike is determined from the intersection. For this reason, even if the electric field change waveform accompanying lightning discharge is complicated or lightning strikes occur at a plurality of locations almost simultaneously, positioning can be performed. The arrival time difference method uses lightning strike location (two-dimensional location) to record the arrival time of electromagnetic waves associated with lightning discharge at three or more receiving stations and determine the location of the lightning strike from the arrival time difference. How to do it. The arrival time difference to each station is recorded with high accuracy using a GPS clock, so lightning strike location can be determined with higher accuracy than the meeting method. Since the position difference accuracy is superior to the arrival time difference method, the position determination is performed by the intersection method only when position determination is difficult by the arrival time difference method in the current system.

現行のシステムでは、各観測局は互いに100km程度以上離れて雷放電に伴う電磁波を観測している。しかし、伝搬距離が100kmを超える長距離だと、伝搬経路の地形や大地導電率の影響を受け、電磁界変化波形は変歪、減衰してしまい、結果として位置標定結果に誤差が発生する。これを防ぐ為には、伝搬距離が短くなるように観測局間の距離を短くすれば良い。しかし、現行の落雷位置標定システムでは、観測対象としている周波数帯が300kHz程度以下と低く、A/D変換のサンプリングレートも0.2μs程度と遅いため、伝搬距離が短くなるように受信局間の間隔を狭くしても位置標定精度は400m程度止まりであった。   In the current system, each observation station observes electromagnetic waves associated with lightning discharge at a distance of about 100 km or more. However, if the propagation distance is longer than 100 km, the electromagnetic field change waveform is distorted and attenuated due to the influence of the propagation path topography and the ground conductivity, resulting in an error in the position location result. In order to prevent this, the distance between observation stations may be shortened so that the propagation distance is shortened. However, in the current lightning location system, the frequency band to be observed is as low as about 300 kHz or less, and the sampling rate of A / D conversion is as low as about 0.2 μs, so that the interval between receiving stations is shortened so that the propagation distance is shortened. Even if it was narrowed, the positioning accuracy was only about 400 m.

これに対し特許文献1では、受信周波数及びサンプリングレートを高め、従来のシステムと比較して比較的近距離の所定半径円内(好ましくは50km)に3局または4局以上の受信局を設置することによって、位置標定精度を200m程度に向上させた。また、精度が向上したことにより、4局以上で観測する場合には三次元位置標定が可能となり、雲内のみの雷放電と落雷との分類も可能にしている。そのため従来の2次元での位置標定システムが一般的に「落雷位置標定システム」とされるのに対して、特許文献1に示す発明は「雷放電位置標定システム」と呼ぶ。   On the other hand, in Patent Document 1, the reception frequency and the sampling rate are increased, and three or four or more receiving stations are installed within a predetermined radius circle (preferably 50 km) that is relatively close compared to the conventional system. As a result, the positioning accuracy was improved to about 200 m. In addition, due to improved accuracy, three-dimensional positioning is possible when observing at four stations or more, and it is also possible to classify lightning discharges and lightning strikes only in the clouds. Therefore, the conventional two-dimensional positioning system is generally referred to as a “lightning strike positioning system”, whereas the invention disclosed in Patent Document 1 is referred to as a “lightning discharge positioning system”.

特開2007−121127号公報JP 2007-121127 A

落雷時に被害を確認する代表的な設備の1つに、送電鉄塔がある。雷放電位置標定システムによって落雷があったとされる位置に送電鉄塔がある場合には、これに登って送電線を確認する必要がある。しかし、特許文献1のシステムでは、位置標定精度が200m程度であり、隣接する鉄塔の間隔が300m程度の通常の場合でも、どの鉄塔に落雷があったか特定できない。この場合は、複数の鉄塔に登って確認作業を行う必要があり、費用と時間が掛かるという問題がある。   One of the typical facilities for checking damage during a lightning strike is a transmission tower. If there is a power transmission tower at a position where lightning strikes are detected by the lightning discharge location system, it is necessary to climb this to check the transmission line. However, in the system of Patent Document 1, even if the positioning accuracy is about 200 m and the interval between adjacent steel towers is about 300 m, it is not possible to specify which steel tower has a lightning strike. In this case, it is necessary to climb up a plurality of steel towers to perform confirmation work, and there is a problem that it takes cost and time.

発明者らが検討したところ、位置標定精度が低下する原因は、地形や大地導電率の影響で電磁界変化波形に変歪が生じ、特許文献1のシステムで電磁波到達時刻としている電界変化波形のピーク時刻に誤差が発生するためである。特許文献1のシステムを用いる事によって、現行のシステムに比べればこの影響は大幅に軽減する事が出来、位置標定精度は向上したが、前述の通り不十分であり、更なる改善が必要である。   As a result of investigations by the inventors, the cause of the decrease in positioning accuracy is that the electromagnetic field change waveform is distorted due to the influence of topography and ground conductivity, and the electric field change waveform that is the electromagnetic wave arrival time in the system of Patent Document 1 is used. This is because an error occurs at the peak time. By using the system of Patent Document 1, this effect can be greatly reduced compared to the current system and the positioning accuracy has been improved, but as described above, it is insufficient and further improvement is required. .

そこで本発明は、到達時間差法を用いた雷放電位置標定システムにおいて、位置標定精度をさらに高めることを目的としている。   Accordingly, an object of the present invention is to further improve the positioning accuracy in a lightning discharge positioning system using the arrival time difference method.

発明者らが検討したところ、数十km程度以内の比較的近距離を電磁波が伝搬した場合でも、落雷に伴う電磁界変化波形のわずかな変歪によるピーク時刻の変化はなお大きく、またピーク時刻の判定は、ラジオ放送波などのノイズの影響も受けやすい。これらのことから、従来のように電磁界変化波形のピークの時刻を電磁波到達時刻と判定すること自体に精度向上の限界があると考えるに至り、これに代わる電磁波到達時刻の判定方法について検討し、本発明を完成するに至った。   The inventors have examined that even when electromagnetic waves propagate in a relatively short distance of about several tens of kilometers, the change in peak time due to slight distortion of the electromagnetic field change waveform caused by lightning strike is still large, and the peak time This determination is easily affected by noise such as radio broadcast waves. From these facts, it has been considered that there is a limit to improving accuracy in determining the peak time of the electromagnetic field change waveform as the electromagnetic wave arrival time as in the past, and an alternative method for determining the electromagnetic wave arrival time has been studied. The present invention has been completed.

すなわち本発明に係る雷放電位置標定システムの代表的な構成は、雷放電に伴う電磁波を含む電磁波を受信し電界の微分信号を出力するアンテナと、微分信号を積分して電界信号を出力する積分回路と、雷放電が発生した時刻である雷放電に伴う電磁波の到達時刻を特定する処理部とを有する複数の受信局と、複数の受信局における雷放電に伴う電磁波の到達時刻から雷放電が発生した位置を標定する中央処理部とを備え、処理部は、電界信号のピークの時刻より十分にさかのぼった所定時間幅の微分信号の平均を取ってゼロレベルを決定し、電界信号のピークよりも前の時刻において微分信号のピークを決定し、微分信号のピークよりも前の時刻において微分信号の値がゼロレベルから微分信号のピークまでの40〜90%の間の所定値となった時刻を雷放電に伴う電磁波の到達時刻とすることを特徴とする。
That is, a typical configuration of the lightning discharge location system according to the present invention includes an antenna that receives an electromagnetic wave including an electromagnetic wave accompanying a lightning discharge and outputs an electric field differential signal, and an integration that integrates the differential signal and outputs an electric field signal. A plurality of receiving stations each having a circuit and a processing unit for identifying an arrival time of the electromagnetic wave accompanying the lightning discharge, which is a time when the lightning discharge occurs, and the lightning discharge from the arrival time of the electromagnetic wave accompanying the lightning discharge at the plurality of receiving stations. A central processing unit for locating the generated position, and the processing unit determines the zero level by taking an average of the differential signal having a predetermined time width that goes back sufficiently from the time of the peak of the electric field signal. The peak of the differential signal is determined at the previous time, and the value of the differential signal is a predetermined value between 40% and 90% from the zero level to the peak of the differential signal at the time before the peak of the differential signal. The time of Tsu, characterized in that the electromagnetic wave arrival time caused by lightning discharge.

上記構成によれば、電磁波の長距離伝搬による波形変歪の影響を極めて低減させることができ、雷放電による電磁波の到達時間差を精度良くかつ安定的に取得することができる。これは、特許文献1で雷放電に伴う電磁波の到達時刻としてきたピーク時刻より時間的に前の部分の波形から到達時刻を決定するため、先行する電界変化波形の変歪の重畳をさらに軽減出来るためである。これにより、雷放電の位置標定精度を飛躍的に向上させることが可能となった。   According to the above configuration, the influence of waveform distortion due to long-distance propagation of electromagnetic waves can be greatly reduced, and the arrival time difference of electromagnetic waves due to lightning discharge can be acquired accurately and stably. This is because the arrival time is determined from the waveform of the portion that is temporally earlier than the peak time that is the arrival time of the electromagnetic wave accompanying lightning discharge in Patent Document 1, so that the superposition of distortion of the preceding electric field change waveform can be further reduced. Because. This makes it possible to dramatically improve the positioning accuracy of lightning discharge.

上記所定値は、ゼロレベルから微分信号のピークまでの75〜85%の間であることが好ましい。上記所定値は電界の微分波形の立ち上がり部分の40〜90%の間で事案に応じて検討、設定すべきである。現在試用しているシステムで得られたデータでは、75%〜85%の間が、S/N(signal-noise ratio)良く到達時刻を定められる確率が高い。なお、40%以下の箇所は長距離伝搬による変歪は低減できるが、信号の変化が緩くS/Nが低下するため適さない。90%以上も変化が緩くS/Nが低下する上、時間的に後になり伝搬の影響による変歪も大きくなり適さない。   The predetermined value is preferably between 75% and 85% from the zero level to the peak of the differential signal. The predetermined value should be examined and set according to the case between 40 and 90% of the rising portion of the differential waveform of the electric field. In the data obtained with the currently used system, the probability of reaching the arrival time with a good S / N (signal-noise ratio) is high between 75% and 85%. Note that distortion of 40% or less can be reduced due to long-distance propagation, but is not suitable because the signal changes slowly and the S / N decreases. The change is not less than 90% and the S / N is lowered, and the distortion due to the influence of propagation is increased later, which is not suitable.

微分信号において、放送波の周波数成分を除去するデジタルフィルタを備えることが好ましい。これにより、ラジオ放送などの放送波によるノイズを低減し、位置標定精度をより高めることができる。   It is preferable to provide a digital filter that removes the frequency component of the broadcast wave in the differential signal. Thereby, the noise by broadcast waves, such as a radio broadcast, can be reduced and position location accuracy can be raised more.

本発明によれば、電磁波の伝搬経路の大地導電率及び地形による波形変歪の影響を低減させることができ、雷放電による電磁波の到達時間差を精度良くかつ安定的に取得することができる。これにより、雷放電の位置標定精度を飛躍的に向上させることができる。   ADVANTAGE OF THE INVENTION According to this invention, the influence of the waveform distortion by the earth conductivity of the propagation path of electromagnetic waves and topography can be reduced, and the arrival time difference of the electromagnetic waves by a lightning discharge can be acquired accurately and stably. Thereby, the positioning accuracy of lightning discharge can be dramatically improved.

実施形態に係る雷放電位置標定システムの原理を説明する図である。It is a figure explaining the principle of the lightning discharge position location system which concerns on embodiment. 受信局の回路構成を説明する図である。It is a figure explaining the circuit structure of a receiving station. 信号波形の例を示す図である。It is a figure which shows the example of a signal waveform. 雷放電時刻の特定を説明するフローチャートである。It is a flowchart explaining specification of the lightning discharge time. 微分信号にデジタルフィルタを適用した例を示す図である。It is a figure which shows the example which applied the digital filter to the differential signal. ゼロレベルの取得を説明する図である。It is a figure explaining acquisition of a zero level. ゼロレベルとゼロクロス点を説明する図である。It is a figure explaining a zero level and a zero crossing point. 雷放電に伴う電磁波の到達時刻とすべき40〜90%の間の所定値をどのように設定するかについて説明する図である。It is a figure explaining how the predetermined value between 40 to 90% which should be the arrival time of the electromagnetic waves accompanying a lightning discharge is set. 比較例と実施例を示す図である。It is a figure which shows a comparative example and an Example.

以下に添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。かかる実施形態に示す寸法、材料、その他具体的な数値などは、発明の理解を容易とするための例示に過ぎず、特に断る場合を除き、本発明を限定するものではない。なお、本明細書及び図面において、実質的に同一の機能、構成を有する要素については、同一の符号を付することにより重複説明を省略し、また本発明に直接関係のない要素は図示を省略する。   Hereinafter, preferred embodiments of the present invention will be described in detail with reference to the accompanying drawings. The dimensions, materials, and other specific numerical values shown in the embodiments are merely examples for facilitating understanding of the invention, and do not limit the present invention unless otherwise specified. In the present specification and drawings, elements having substantially the same function and configuration are denoted by the same reference numerals, and redundant description is omitted, and elements not directly related to the present invention are not illustrated. To do.

図1は実施形態に係る雷放電位置標定システムの原理を説明する図、図2は受信局の回路構成を説明する図である。本実施形態に係る雷放電位置標定システムは、到達時間差法を用いる。   FIG. 1 is a diagram illustrating the principle of a lightning discharge location system according to the embodiment, and FIG. 2 is a diagram illustrating a circuit configuration of a receiving station. The lightning discharge location system according to the present embodiment uses the arrival time difference method.

図1に示す雷放電位置標定システム100は、雷放電に伴う電磁波をアンテナで受信し、雷放電に伴う電磁波到達時刻を特定するための受信局102(102a、102b、102c、102d、…)を半径50kmの円内に4局以上設置する。各受信局102は電界信号を取得し、この電界信号に基づいて後述するように雷放電に伴う電磁波到達時刻を特定し、ホスト装置130(図2参照)に送信する。ホスト装置130は少なくとも4局の受信局102における雷放電に伴う電磁波の到達時刻を用いて、三次元の到達時間差法により雷放電が発生した三次元位置Pを標定する。   A lightning discharge position locating system 100 shown in FIG. 1 receives electromagnetic waves associated with lightning discharges with an antenna, and receives receiving stations 102 (102a, 102b, 102c, 102d,...) For specifying electromagnetic wave arrival times associated with lightning discharges. Install 4 stations or more in a circle with a radius of 50 km. Each receiving station 102 acquires an electric field signal, specifies an electromagnetic wave arrival time associated with lightning discharge based on the electric field signal, and transmits it to the host device 130 (see FIG. 2). The host device 130 locates the three-dimensional position P where the lightning discharge has occurred by the three-dimensional arrival time difference method using the arrival times of the electromagnetic waves accompanying the lightning discharge at at least four receiving stations 102.

ホスト装置130は、全ての受信局102の設置位置をx、y、z軸の三次元座標で記憶しており、各受信局102での雷放電に伴う電磁波の到達時刻を三次元の到達時間差法に適用することで、雷放電が発生した三次元位置Pを標定することができる。   The host device 130 stores the installation positions of all the receiving stations 102 in the three-dimensional coordinates of the x, y, and z axes, and the arrival time of the electromagnetic waves accompanying the lightning discharge at each receiving station 102 is a three-dimensional arrival time difference. By applying the method, the three-dimensional position P where the lightning discharge has occurred can be determined.

図2に示すように、受信局102は、電界アンテナ104、VLF−MF帯の電界アンプ106、A/Dコンバータ108、高精度時計110、処理部112、積分回路114、トリガ回路116、通信部118、デジタルフィルタ120を備える。   2, the receiving station 102 includes an electric field antenna 104, a VLF-MF band electric field amplifier 106, an A / D converter 108, a high-precision clock 110, a processing unit 112, an integration circuit 114, a trigger circuit 116, and a communication unit. 118, a digital filter 120 is provided.

電界アンテナ104は、VLF−MF帯の電磁波を受信し、電界Eの微分信号dE/dtを受信信号として出力する。ただし、電界アンテナ104の出力は微弱であるので、後段におけるサンプリングに適したダイナミックレンジにするために、次の電界アンプ106を使用する。   The electric field antenna 104 receives an electromagnetic wave in the VLF-MF band and outputs a differential signal dE / dt of the electric field E as a reception signal. However, since the output of the electric field antenna 104 is weak, the following electric field amplifier 106 is used in order to obtain a dynamic range suitable for sampling in the subsequent stage.

アンテナは磁界アンテナでもよい。その場合、磁界アンテナは電磁波の磁界の微分信号を受信信号として出力することになる。電界アンテナ104を用いるか、磁界アンテナを用いるかは受信局設置場所の条件により選択してもよい。なお、磁界アンテナの場合は、地面に平行で互いに直交する2成分を同時に観測するが、回路構成は図2に準じる。   The antenna may be a magnetic field antenna. In that case, the magnetic field antenna outputs a differential signal of the magnetic field of the electromagnetic wave as a reception signal. Whether to use the electric field antenna 104 or the magnetic field antenna may be selected according to the conditions of the receiving station installation location. In the case of a magnetic field antenna, two components that are parallel to the ground and orthogonal to each other are observed at the same time, but the circuit configuration conforms to FIG.

電界アンプ106は、電界アンテナ104の出力を増幅する。電界アンプ106は、少なくともVLF−MF帯(3kHz〜3MHz)の帯域で減衰なく増幅ができる周波数特性を有する。なお、VLF帯は3kHz〜30kHz、LF帯は30kHz〜300kHz、MF帯は300kHz〜3MHz、HF帯は3MHz〜30MHz、VHF帯は30MHz〜300MHzである。   The electric field amplifier 106 amplifies the output of the electric field antenna 104. The electric field amplifier 106 has frequency characteristics that can be amplified without attenuation in at least the VLF-MF band (3 kHz to 3 MHz). The VLF band is 3 kHz to 30 kHz, the LF band is 30 kHz to 300 kHz, the MF band is 300 kHz to 3 MHz, the HF band is 3 MHz to 30 MHz, and the VHF band is 30 MHz to 300 MHz.

A/Dコンバータ108は、電界アンプ106が出力した微分信号を時系列のデジタルデータに変換する。A/Dコンバータ108のサンプリング周期は、100ns以下である。   The A / D converter 108 converts the differential signal output from the electric field amplifier 106 into time-series digital data. The sampling period of the A / D converter 108 is 100 ns or less.

高精度時計110は、他の受信局に収容されている同様の装置と時刻同期を図るためのものであり、例えば、ルビジウム高精度時計からなる。これらの高精度時計110は、GPSシステム等を利用して100ns以内の精度で同期を取るようにしてもよい。   The high-precision timepiece 110 is for achieving time synchronization with a similar device accommodated in another receiving station, and is composed of, for example, a rubidium high-precision timepiece. These high-precision clocks 110 may be synchronized with accuracy within 100 ns using a GPS system or the like.

処理部112は、主として雷放電に伴う電磁波の各局への到達時刻の特定を行う。処理の詳細は後述する。処理部112は、具体例としては汎用のコンピュータ上で動作するソフトウェアによって実現することができる。処理部112は、トリガ回路116による雷放電イベント検出時にのみ、つまりトリガ信号が処理部112に到着したときのみ、雷放電による各局への電磁波の到達時刻の特定を実行する。   The processing unit 112 mainly specifies the arrival time of the electromagnetic wave accompanying the lightning discharge to each station. Details of the processing will be described later. As a specific example, the processing unit 112 can be realized by software operating on a general-purpose computer. The processing unit 112 specifies the arrival time of the electromagnetic wave to each station due to lightning discharge only when a lightning discharge event is detected by the trigger circuit 116, that is, only when the trigger signal arrives at the processing unit 112.

また本実施形態では、A/Dコンバータ108と処理部112の間にデジタルフィルタ120を備えている。デジタルフィルタ120は、ラジオ放送などの既知の放送波の成分を除去する。具体的には、デジタルフィルタ120によって放送波の周波数(搬送波の周波数)を減衰させる。このとき、各受信局102に到達する各放送波の電界強度をあらかじめ取得しておくことにより、適切な強度に減衰させることが可能となる。これにより、放送波によるノイズを低減し、より位置標定精度を高めることができる。   In this embodiment, a digital filter 120 is provided between the A / D converter 108 and the processing unit 112. The digital filter 120 removes known broadcast wave components such as radio broadcasts. Specifically, the frequency of the broadcast wave (carrier frequency) is attenuated by the digital filter 120. At this time, it is possible to attenuate to an appropriate intensity by acquiring in advance the electric field intensity of each broadcast wave that reaches each receiving station 102. Thereby, the noise by a broadcast wave can be reduced and position location accuracy can be improved more.

積分回路114は、電界アンプ106で増幅された受信信号を積分する。受信信号が電界の微分信号dE/dtであるから、これを積分した積分回路114の出力は電界Eの大きさ及び極性を表した電界信号である。   The integration circuit 114 integrates the reception signal amplified by the electric field amplifier 106. Since the received signal is the electric field differential signal dE / dt, the output of the integrating circuit 114 obtained by integrating the received signal is an electric field signal representing the magnitude and polarity of the electric field E.

トリガ回路116は、電界信号が雷放電の発生と関連する有意レベルに達したかどうかを判定する回路である。   The trigger circuit 116 is a circuit that determines whether the electric field signal has reached a significant level associated with the occurrence of lightning discharge.

図3は信号波形の例を示す図であって、図3(a)は微分信号の波形を示す図、図3(b)は電界信号の波形を示す図、図3(c)は微分信号の波形の拡大図、図3(d)は電界信号の波形の拡大図である。図3(a)に示されるように、受信信号dE/dtは、振幅や周期が様々なパルス状の波形が多く重畳している。そのため、dE/dtでは電界強度が雷放電の発生と関連する強度に達したか判断が難しい。そこで、積分回路114が受信信号を積分することで電界信号を作ると、図3(b)に示すような重畳するパルスの少ない電界信号の波形が得られる。そしてトリガ回路116は、電界信号の強度が雷放電の発生と認められるほど有意レベルに達したかどうかを判定する。つまり、トリガ回路116は、電界信号(アナログ信号)強度がしきい値を越えたことで雷放電があったというイベントを検出する。トリガ回路116の出力であるトリガ信号は、処理部112に入力される。   3A and 3B are diagrams illustrating examples of signal waveforms, in which FIG. 3A illustrates a waveform of a differential signal, FIG. 3B illustrates a waveform of an electric field signal, and FIG. 3C illustrates a differential signal. FIG. 3D is an enlarged view of the waveform of the electric field signal. As shown in FIG. 3A, the reception signal dE / dt has many superimposed pulse waveforms with various amplitudes and periods. Therefore, in dE / dt, it is difficult to determine whether the electric field intensity has reached the intensity associated with the occurrence of lightning discharge. Therefore, when the integration circuit 114 generates an electric field signal by integrating the received signal, a waveform of the electric field signal with few superimposed pulses as shown in FIG. 3B is obtained. The trigger circuit 116 determines whether or not the intensity of the electric field signal has reached a significant level such that lightning discharge is recognized. That is, the trigger circuit 116 detects an event that a lightning discharge has occurred due to the electric field signal (analog signal) intensity exceeding a threshold value. A trigger signal that is an output of the trigger circuit 116 is input to the processing unit 112.

なお、図3(a)に示す微分信号のピーク波形122と、図3(b)に示す電界信号のピーク波形124とは、おおむね同じ位置にある。しかし詳細には、図3(c)、図3(d)に示すように、微分信号のピークは電界信号の最も急峻な時刻に対応し、電界信号のピークは微分信号が減衰してゼロとなる時刻に対応する。   The peak waveform 122 of the differential signal shown in FIG. 3A and the peak waveform 124 of the electric field signal shown in FIG. 3B are substantially at the same position. However, in detail, as shown in FIGS. 3C and 3D, the peak of the differential signal corresponds to the steepest time of the electric field signal, and the peak of the electric field signal is zero as the differential signal attenuates. Corresponds to the time.

通信部118は、処理部112が出力した雷放電に伴う電磁波到達時刻および電界信号のピーク波高値を、インターネットやISDNなどのデータ通信網128を介してホスト装置130に送信する。   The communication unit 118 transmits the arrival time of the electromagnetic wave accompanying lightning discharge and the peak value of the electric field signal output from the processing unit 112 to the host device 130 via the data communication network 128 such as the Internet or ISDN.

ホスト装置130は、複数の受信局102からデータを受信する通信部132と、中央処理部134とを備えている。ホスト装置130は、具体例としては汎用のコンピュータを用いることができ、中央処理部134はそのコンピュータ上で動作するソフトウェアによって実現することができる。   The host device 130 includes a communication unit 132 that receives data from a plurality of receiving stations 102 and a central processing unit 134. The host device 130 can use a general-purpose computer as a specific example, and the central processing unit 134 can be realized by software operating on the computer.

中央処理部134は、図1の各受信局102(102a、102b、102c、102d)から、それぞれ雷放電に伴う電磁波の受信時刻t(t1、t2、t3、t4)を取得する。中央処理部134は、全ての受信局102の設置位置をx、y、z軸の三次元座標で記憶しており、各受信局102の設置位置を示した三次元空間において、各受信局102(102a、102b、102c、102d)における雷放電に伴う電磁波の受信時刻t(t1、t2、t3、t4)を三次元の到達時間差法に適用することで、雷放電の三次元位置Pを標定する。   The central processing unit 134 acquires the reception time t (t1, t2, t3, t4) of the electromagnetic wave accompanying the lightning discharge from each receiving station 102 (102a, 102b, 102c, 102d) in FIG. The central processing unit 134 stores the installation positions of all the reception stations 102 in three-dimensional coordinates of the x, y, and z axes, and in the three-dimensional space indicating the installation positions of the reception stations 102, each reception station 102. By applying the reception time t (t1, t2, t3, t4) of the electromagnetic wave accompanying the lightning discharge at (102a, 102b, 102c, 102d) to the three-dimensional arrival time difference method, the three-dimensional position P of the lightning discharge is determined. To do.

具体的には、中央処理部134は、受信局102a、102b間の到達時間差t1−t2をもたらす点の集合として回転双曲面S1(図示困難のため図示せず)を計算する。別の受信局の組からも同様に到達時間差t2−t3、t3−t4に基づいた回転双曲面S2、S3(図示せず)を計算する。これらの回転双曲面S1、S2、S3の交点が三次元位置Pであるから、ホスト装置は、回転双曲面S1、S2の交線K12を計算し、回転双曲面S2、S3の交線K23を計算し、交線K12と交線K23との交点である三次元位置Pを計算する。   Specifically, the central processing unit 134 calculates a rotating hyperboloid S1 (not shown due to difficulty in illustration) as a set of points that cause the arrival time difference t1-t2 between the receiving stations 102a and 102b. Similarly, rotation hyperboloids S2 and S3 (not shown) based on arrival time differences t2-t3 and t3-t4 are calculated from another set of receiving stations. Since the intersection of these rotation hyperboloids S1, S2, and S3 is the three-dimensional position P, the host device calculates the intersection line K12 of the rotation hyperboloids S1 and S2, and sets the intersection line K23 of the rotation hyperboloids S2 and S3. The three-dimensional position P that is the intersection of the intersection line K12 and the intersection line K23 is calculated.

回転双曲面S1、S2、S3が図示できないので、参考のため、回転双曲面S1と大地(XY平面)との交線K10、回転双曲面S2と大地との交線K20、回転双曲面S3と大地との交線K30が示してある。これらの交線K10、K20、K30はいずれも双曲線となる。   Since the rotation hyperboloids S1, S2, and S3 cannot be illustrated, for reference, the intersection line K10 of the rotation hyperboloid S1 and the ground (XY plane), the intersection line K20 of the rotation hyperboloid S2 and the ground, the rotation hyperboloid S3, and An intersection line K30 with the ground is shown. These intersecting lines K10, K20, and K30 are all hyperbolic curves.

実際の計算においては、観測領域全体をメッシュに区切り、各メッシュ点でのchi−square検定による誤差評価を行い、最適な地点を探し出すことが好ましい。chi−square検定の式は、以下の通りである。

Figure 2012189387

上記の式を用いて、各受信局102のchi−squareを合算し、最小となる位置、時刻を特定する。なお、最小となる位置の探索方法については、通常の最急降下法では最小解ではなく点在する極に入り込む可能性がある。極に入ることを防止するためには全てのメッシュ点でchi−squareを計算することが理想的だが、それでは計算量が膨大になるため、現実的ではない。そこで、探索区域をある程度粗いメッシュに区切り、そのメッシュでの最小点を探索し、その点を中心として更に細かいメッシュで探索する。これによりローカルな極に入り込むことを防止し、迅速に適切な解を得ることができる。 In actual calculation, it is preferable to divide the entire observation region into meshes, perform error evaluation by chi-square test at each mesh point, and find an optimum point. The formula of the chi-square test is as follows.
Figure 2012189387

Using the above equation, the chi-square of each receiving station 102 is summed to specify the minimum position and time. As for the search method of the minimum position, the normal steepest descent method may enter the scattered poles instead of the minimum solution. Although it is ideal to calculate chi-square at all mesh points in order to prevent entering the pole, this is not realistic because the calculation amount becomes enormous. Therefore, the search area is divided into meshes that are coarse to some extent, the minimum point in the mesh is searched, and the search is performed with a finer mesh centered on that point. As a result, it is possible to prevent entry into a local pole and obtain an appropriate solution quickly.

なお、本実施形態では受信局102の数を4つと説明している。受信局が4つあれば到達時間差法による三次元の雷放電位置標定をすることが可能である。しかし、受信局102の数は多いほど位置標定誤差の標準偏差は小さくできる。また、受信局102が5局以上あれば、chi−squareの計算を行ったときに値の大きくなってしまったデータ(他の落雷を観測、S/Nの悪いデータ)を除外して再計算することも可能である。これより、受信局102の数は5局以上の方が位置標定精度を高めることができる。   In the present embodiment, the number of receiving stations 102 is described as four. If there are four receiving stations, it is possible to determine the three-dimensional lightning discharge position by the arrival time difference method. However, the standard deviation of the positioning error can be reduced as the number of receiving stations 102 increases. Also, if there are 5 or more receiving stations 102, recalculate by excluding data that has become larger when chi-square is calculated (observation of other lightning strikes, data with poor S / N). It is also possible to do. As a result, the number of receiving stations 102 can be higher when the number of receiving stations is five or more.

上記構成の雷放電位置標定システム100において、本実施形態の特徴は、受信局102の処理部112における、雷放電に伴う電磁波の受信局102への到達時刻の特定にある。   In the lightning discharge position locating system 100 configured as described above, the feature of the present embodiment is that the processing unit 112 of the receiving station 102 specifies the arrival time of the electromagnetic wave accompanying the lightning discharge to the receiving station 102.

上述の特許文献1(特開2007−121127号公報)では、微分信号(電界微分波形)のゼロクロス値を到達時刻としている。これは、積分した電界信号の最大ピークに該当する。しかしながら、電界信号の最大ピーク時刻を用いると、長距離伝搬による変歪によって、各観測局への到達時間差に誤差が発生する。   In the above-mentioned Patent Document 1 (Japanese Patent Laid-Open No. 2007-121127), the zero cross value of the differential signal (electric field differential waveform) is set as the arrival time. This corresponds to the maximum peak of the integrated electric field signal. However, when the maximum peak time of the electric field signal is used, an error occurs in the arrival time difference to each observation station due to distortion due to long-distance propagation.

そこで本実施形態では、微分信号の立ち上がり部分の一点の時刻を到達時刻とする。なお、この立ち上がりの部分の一点はネットワークの構成、場所、季節により最適な位置が変化するため、微分信号の値がゼロレベルから微分信号のピークまでの40〜90%の間で逐次調整を行う。初期値はピークの80%の位置とする。   Therefore, in this embodiment, the time at one point of the rising portion of the differential signal is set as the arrival time. Note that the optimum position of this rising portion changes depending on the network configuration, location, and season, so that the differential signal value is sequentially adjusted between 40% and 90% from the zero level to the peak of the differential signal. . The initial value is 80% of the peak.

図4は雷放電時刻の特定を説明するフローチャートである。まずアンテナ104によって微分信号を受信すると(ステップ200)、積分回路114によって電界信号を取得する(ステップ202)。そしてトリガ回路116によって雷放電イベントを検出すると(ステップ204)、処理部112は雷放電に伴う電磁波の到達時刻の特定を開始する(ステップ206)。具体的には、処理部112は常時微分信号を蓄積しており、雷放電イベントが来ると、そのイベントの時刻から所定時間さかのぼった時点からの微分信号に対して処理を行う。   FIG. 4 is a flowchart for explaining the specification of the lightning discharge time. First, when a differential signal is received by the antenna 104 (step 200), an electric field signal is acquired by the integrating circuit 114 (step 202). When the lightning discharge event is detected by the trigger circuit 116 (step 204), the processing unit 112 starts specifying the arrival time of the electromagnetic wave accompanying the lightning discharge (step 206). Specifically, the processing unit 112 always accumulates a differential signal, and when a lightning discharge event occurs, the processing unit 112 performs processing on the differential signal from a time point that goes back a predetermined time from the time of the event.

処理部112は、まず微分信号のゼロレベルを取得する必要がある。しかし、受信した微分信号には多くのノイズが乗っている。ノイズの原因は様々であるが、最も影響が大きいと考えられる放送波によるノイズをデジタルフィルタ120によって除去する(ステップ208)。図5は微分信号にデジタルフィルタ120を適用した例を示す図である。図5に示すように、放送波によるノイズを低減することによって効果的にノイズを低減することができる。   The processing unit 112 first needs to acquire the zero level of the differential signal. However, a lot of noise is on the received differential signal. There are various causes of noise, but noise due to broadcast waves considered to have the greatest influence is removed by the digital filter 120 (step 208). FIG. 5 is a diagram showing an example in which the digital filter 120 is applied to the differential signal. As shown in FIG. 5, noise can be effectively reduced by reducing noise caused by broadcast waves.

次に処理部112は、微分信号のゼロレベルを取得する(ステップ210)。ゼロレベルは雷放電より前の時刻、すなわち電界信号のピークの時刻よりも十分にさかのぼった所定時間幅の微分信号の平均を取ってゼロレベルを決定する。これは雷放電による電界信号が観測されている間はゼロレベルが決められないためである。   Next, the processing unit 112 acquires the zero level of the differential signal (step 210). The zero level is determined by taking an average of the differential signal having a predetermined time width that goes back sufficiently before the time of lightning discharge, that is, the time of the peak of the electric field signal. This is because the zero level cannot be determined while the electric field signal due to lightning discharge is observed.

図6はゼロレベルの取得を説明する図である。図6では、微分信号の始め(雷放電イベントから所定時間さかのぼった時刻のデータ)から、100データ(4μs)、300データ(12μs)、500データ(20μs)で移動平均を取っている。1データは40nsである。図6から、100データや300データでは変動が大きく、500データでは変動が小さくなっていることから、これを用いてゼロレベルを決定することができる。   FIG. 6 is a diagram for explaining the acquisition of the zero level. In FIG. 6, the moving average is taken with 100 data (4 μs), 300 data (12 μs), and 500 data (20 μs) from the beginning of the differential signal (data at a time that goes back a predetermined time from the lightning discharge event). One data is 40 ns. From FIG. 6, since the fluctuation is large in 100 data and 300 data and the fluctuation is small in 500 data, the zero level can be determined using this.

次に処理部112は、微分信号とゼロレベルとの交点であるゼロクロス点を算出する(ステップ212)。図7は前述のゼロレベルとゼロクロス点を説明する図である。図7に示すように、ゼロクロス点となる時刻は、電界信号のピークの時刻となる。   Next, the processing unit 112 calculates a zero cross point that is an intersection of the differential signal and the zero level (step 212). FIG. 7 is a diagram for explaining the above-described zero level and zero cross point. As shown in FIG. 7, the time when the zero cross point is reached is the peak time of the electric field signal.

次に処理部112は、ゼロクロス点よりも前の時刻において微分信号のピークを検出する(ステップ214)。微分信号のピークは、電界信号が最も急峻となっていることを意味する。複数のピークが存在する場合があるが、ここでは単に最大のピークを検出する。   Next, the processing unit 112 detects the peak of the differential signal at a time before the zero cross point (step 214). The peak of the differential signal means that the electric field signal is the steepest. Although there may be a plurality of peaks, the maximum peak is simply detected here.

そして処理部112は、微分信号のピークよりも前の時刻において、微分信号の値がゼロレベルから微分信号のピークまでの40〜90%の間の所定値となった時刻を雷放電に伴う電磁波の到達時刻とする(ステップ216)。   Then, the processing unit 112 detects the time when the value of the differential signal becomes a predetermined value between 40% and 90% from the zero level to the peak of the differential signal at the time before the peak of the differential signal. (Step 216).

すなわち、電界信号のピークの時刻(微分信号のゼロクロス点)ではなく、微分信号のピークの時刻でもなく、微分信号の立ち上がり部分の一点の時刻を到達時刻とする。微分信号が急激に立ち上がる箇所ではS/Nが良い為、安定した到達時刻を定めることができる。これにより、複数の受信局102の間で雷放電に伴う電磁波の到達時刻の差を精度良く、かつ安定的に取得することが可能となる。   That is, the arrival time is not the peak time of the electric field signal (zero cross point of the differential signal) but the time of the rising portion of the differential signal, not the peak time of the differential signal. Since the S / N is good where the differential signal rises rapidly, a stable arrival time can be determined. Accordingly, it is possible to accurately and stably acquire the difference in arrival time of electromagnetic waves accompanying lightning discharge among the plurality of receiving stations 102.

図8は雷放電に伴う電磁波の到達時刻とすべき40〜90%の間の所定値をどのように設定するかについて説明する図である。図8に示すデータは本発明の検証のために試験的に稼動しているシステムで得られたデータであり、夏季に栃木県周辺で観測されたものである。図8(a)は所定値の比率に対する標定位置誤差を示す図、図8(b)は図8(a)に示すデータ2の微分信号の波形を示す図、図8(c)は推定電流値を示す図、図8(d)はchi−squareを示す図である。   FIG. 8 is a diagram for explaining how to set a predetermined value between 40% and 90% that should be the arrival time of the electromagnetic wave accompanying lightning discharge. The data shown in FIG. 8 is data obtained by a system operating on a trial basis for verification of the present invention, and was observed in the vicinity of Tochigi Prefecture in the summer. 8A is a diagram showing the orientation position error with respect to the ratio of the predetermined value, FIG. 8B is a diagram showing the waveform of the differential signal of the data 2 shown in FIG. 8A, and FIG. 8C is the estimated current. The figure which shows a value and FIG.8 (d) is a figure which shows chi-square.

図8(a)に示す表の右端にある合計の列を参照すると、40%のときの誤差が大きく、60%と80%の誤差は同等である。したがって、所定値を60%としても、80%としても同等の結果を得ることができることがわかる。ただし個別にデータを参照すると、図8(b)に示すように、データ2では60%の位置を複数回通過していて、7データ(280ns)のずれがある位置にも60%を通過するデータが存在していた。なお、図8(a)では1度目に通過した時刻で計算している。このようなデータがある場合には、80%の位置で計算した方が信頼性が高いと考えられる。   Referring to the total column at the right end of the table shown in FIG. 8A, the error at 40% is large, and the errors at 60% and 80% are equivalent. Therefore, it can be seen that the same result can be obtained even if the predetermined value is 60% or 80%. However, referring to the data individually, as shown in FIG. 8B, the data 2 passes through the position of 60% a plurality of times, and also passes through 60% to the position where there is a deviation of 7 data (280 ns). Data existed. In FIG. 8 (a), the calculation is performed at the time of passing the first time. When there is such data, it is considered that the reliability is higher when the calculation is performed at the position of 80%.

また図8(c)を参照すると、データ1の電流値が比較的大きめである。データ1は受信された電界変化波形の強度も強いため、S/Nが良く、ノイズに影響されず本発明の妥当性を検証するには図8のデータの中では最も適している。そのため、これに注目して図8(d)を参照する。すると、データ1では60%より80%のほうがchi−squareの値が小さい。このことからも、このときの一連のデータでは所定値として80%を用いることが適切である。   Further, referring to FIG. 8C, the current value of data 1 is relatively large. Since the intensity of the received electric field change waveform is strong, the S / N is good and the data 1 is most suitable among the data in FIG. 8 to verify the validity of the present invention without being affected by noise. Therefore, paying attention to this, reference is made to FIG. Then, in data 1, the value of chi-square is smaller in 80% than in 60%. For this reason as well, it is appropriate to use 80% as the predetermined value in the series of data at this time.

なお上述したように、この立ち上がりの部分の一点はネットワークの構成、場所、季節により最適な位置が変化するため、微分信号の値がゼロレベルから微分信号のピークまでの40〜90%の間で逐次調整を行うことが望ましい。ただし、特に75%〜85%の間ではS/Nが向上し、安定して所定値を設定することができると推定される。   As described above, since the optimal position of this rising portion changes depending on the network configuration, location, and season, the differential signal value is between 40% and 90% from the zero level to the peak of the differential signal. It is desirable to perform sequential adjustment. However, it is estimated that the S / N is improved particularly between 75% and 85%, and the predetermined value can be set stably.

(実施例と比較例)
図9は比較例と実施例を示す図である。比較例としては東京電力で運用されている従来からの落雷位置標定システム、実施例としては上記実施形態で説明した雷放電位置標定システムを用いた。同じ雷のデータ(いずれも落雷位置が確認されたもの)について比較例と実施例のシステムで雷放電位置を算出したところ、比較例においては平均誤差が367m、最大誤差では759mとなり、実施例では平均誤差が80m、最大誤差が230mであった。したがって、平均誤差で1/5程度、最大誤差で1/3程度と、大幅な精度向上が確認された。
(Examples and comparative examples)
FIG. 9 is a diagram showing a comparative example and an example. As a comparative example, a conventional lightning strike location system operated by TEPCO, and as an example, the lightning discharge location system described in the above embodiment was used. When the lightning discharge position was calculated with the system of the comparative example and the example for the same lightning data (both where the lightning strike position was confirmed), the average error in the comparative example was 367 m, and the maximum error was 759 m. The average error was 80 m and the maximum error was 230 m. Therefore, a significant improvement in accuracy was confirmed, with an average error of about 1/5 and a maximum error of about 1/3.

以上、添付図面を参照しながら本発明の好適な実施形態について説明したが、本発明は係る例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。   As mentioned above, although preferred embodiment of this invention was described referring an accompanying drawing, it cannot be overemphasized that this invention is not limited to the example which concerns. It will be apparent to those skilled in the art that various changes and modifications can be made within the scope of the claims, and these are naturally within the technical scope of the present invention. Understood.

本発明は、雷放電に伴う電磁波到達時刻を用いて雷放電が発生した位置を標定する雷放電位置標定システムとして利用することができる。   INDUSTRIAL APPLICABILITY The present invention can be used as a lightning discharge position locating system for locating a position where a lightning discharge has occurred using an electromagnetic wave arrival time associated with a lightning discharge.

100…雷放電位置標定システム、102…受信局、104…アンテナ、106…電界アンプ、108…A/Dコンバータ、110…高精度時計、112…処理部、114…積分回路、116…トリガ回路、118…通信部、120…デジタルフィルタ、122…ピーク波形、124…ピーク波形、128…データ通信網、130…ホスト装置、132…通信部、134…中央処理部、t…雷放電に伴う電磁波受信時刻 DESCRIPTION OF SYMBOLS 100 ... Lightning discharge location system, 102 ... Receiving station, 104 ... Antenna, 106 ... Electric field amplifier, 108 ... A / D converter, 110 ... High precision clock, 112 ... Processing part, 114 ... Integration circuit, 116 ... Trigger circuit, DESCRIPTION OF SYMBOLS 118 ... Communication part, 120 ... Digital filter, 122 ... Peak waveform, 124 ... Peak waveform, 128 ... Data communication network, 130 ... Host device, 132 ... Communication part, 134 ... Central processing part, t ... Electromagnetic wave reception accompanying lightning discharge Times of Day

Claims (3)

雷放電に伴う電磁波を含む電磁波を受信し電界の微分信号を出力するアンテナと、前記微分信号を積分して電界信号を出力する積分回路と、雷放電が発生した時刻である雷放電に伴う電磁波の到達時刻を特定する処理部とを有する複数の受信局と、
前記複数の受信局における雷放電に伴う電磁波の到達時刻から雷放電が発生した位置を標定する中央処理部とを備え、
前記処理部は、
前記電界信号のピークの時刻より十分にさかのぼった所定時間幅の前記微分信号の平均を取ってゼロレベルを決定し、
前記電界信号のピークよりも前の時刻において前記微分信号のピークを決定し、
前記微分信号のピークよりも前の時刻において微分信号の値が前記ゼロレベルから前記微分信号のピークまでの40〜90%の間の所定値となった時刻を雷放電に伴う電磁波の到達時刻とすることを特徴とする雷放電位置標定システム。
An antenna that receives an electromagnetic wave including an electromagnetic wave accompanying a lightning discharge and outputs a differential signal of an electric field, an integration circuit that integrates the differential signal and outputs an electric field signal, and an electromagnetic wave accompanying a lightning discharge that is the time when the lightning discharge occurs A plurality of receiving stations having a processing unit for identifying the arrival time of
A central processing unit for locating the position where the lightning discharge occurred from the arrival time of the electromagnetic wave accompanying the lightning discharge in the plurality of receiving stations,
The processor is
Determining the zero level by taking an average of the differential signal of a predetermined time width that goes back sufficiently from the peak time of the electric field signal;
Determining the peak of the differential signal at a time prior to the peak of the electric field signal;
The time when the value of the differential signal becomes a predetermined value between 40% and 90% from the zero level to the peak of the differential signal at the time before the peak of the differential signal is the arrival time of the electromagnetic wave accompanying lightning discharge Lightning discharge location system characterized by
前記所定値は、好ましくは前記ゼロレベルから前記微分信号のピークまでの75〜85%の間であることを特徴とする請求項1に記載の雷放電位置標定システム。   The lightning discharge position locating system according to claim 1, wherein the predetermined value is preferably between 75% and 85% from the zero level to the peak of the differential signal. 前記微分信号において、放送波の周波数成分を除去するデジタルフィルタを備えることを特徴とする請求項1に記載の雷放電位置標定システム。   The lightning discharge position locating system according to claim 1, further comprising a digital filter that removes a frequency component of a broadcast wave in the differential signal.
JP2011051799A 2011-03-09 2011-03-09 Lightning discharge position orientation system Withdrawn JP2012189387A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011051799A JP2012189387A (en) 2011-03-09 2011-03-09 Lightning discharge position orientation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011051799A JP2012189387A (en) 2011-03-09 2011-03-09 Lightning discharge position orientation system

Publications (1)

Publication Number Publication Date
JP2012189387A true JP2012189387A (en) 2012-10-04

Family

ID=47082737

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011051799A Withdrawn JP2012189387A (en) 2011-03-09 2011-03-09 Lightning discharge position orientation system

Country Status (1)

Country Link
JP (1) JP2012189387A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104569913A (en) * 2015-01-30 2015-04-29 武汉大学 High-precision full-lightning positioning method
CN105738712A (en) * 2016-03-18 2016-07-06 福克斯通(北京)科技有限公司 Fracturing truck set lightning early warning method, device and system
CN105866558A (en) * 2016-03-23 2016-08-17 清华大学深圳研究生院 Small-region lightning monitoring method and system
JP2017156158A (en) * 2016-02-29 2017-09-07 株式会社フランクリン・ジャパン Lightning position locating system, lightning position locating device and lightning position locating method
JP2017181039A (en) * 2016-03-28 2017-10-05 一般財団法人電力中央研究所 Estimation method, estimation device, and estimation program for lightening strike current characteristics value
JP2018054367A (en) * 2016-09-27 2018-04-05 株式会社フランクリン・ジャパン Electric charge existence altitude calculation system, and electric charge existence altitude calculation device
CN110174557A (en) * 2019-01-11 2019-08-27 南京信息工程大学 A kind of thunder cloud positioning calibration method based on three-dimensional atmospheric electric field instrument observation visual angle
CN112766127A (en) * 2021-01-12 2021-05-07 南京信息工程大学 Thundercloud point charge positioning method based on complementary set modal decomposition and SG filtering
CN113295933A (en) * 2021-05-25 2021-08-24 云南电网有限责任公司电力科学研究院 Group velocity extraction method for lightning electromagnetic wave propagation

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104569913A (en) * 2015-01-30 2015-04-29 武汉大学 High-precision full-lightning positioning method
JP2017156158A (en) * 2016-02-29 2017-09-07 株式会社フランクリン・ジャパン Lightning position locating system, lightning position locating device and lightning position locating method
CN105738712A (en) * 2016-03-18 2016-07-06 福克斯通(北京)科技有限公司 Fracturing truck set lightning early warning method, device and system
CN105866558B (en) * 2016-03-23 2018-06-12 清华大学深圳研究生院 A kind of zonule lightning monitoring localization method and system
CN105866558A (en) * 2016-03-23 2016-08-17 清华大学深圳研究生院 Small-region lightning monitoring method and system
WO2017161630A1 (en) * 2016-03-23 2017-09-28 清华大学深圳研究生院 Small-region lightning monitoring and positioning method and system
JP2017181039A (en) * 2016-03-28 2017-10-05 一般財団法人電力中央研究所 Estimation method, estimation device, and estimation program for lightening strike current characteristics value
JP2018054367A (en) * 2016-09-27 2018-04-05 株式会社フランクリン・ジャパン Electric charge existence altitude calculation system, and electric charge existence altitude calculation device
CN110174557A (en) * 2019-01-11 2019-08-27 南京信息工程大学 A kind of thunder cloud positioning calibration method based on three-dimensional atmospheric electric field instrument observation visual angle
CN110174557B (en) * 2019-01-11 2021-12-21 南京信息工程大学 Thunderstorm cloud positioning and calibrating method based on three-dimensional atmospheric electric field instrument observation visual angle
CN112766127A (en) * 2021-01-12 2021-05-07 南京信息工程大学 Thundercloud point charge positioning method based on complementary set modal decomposition and SG filtering
CN112766127B (en) * 2021-01-12 2023-08-22 南京信息工程大学 Lei Yundian charge positioning method based on complementary set modal decomposition and SG filtering
CN113295933A (en) * 2021-05-25 2021-08-24 云南电网有限责任公司电力科学研究院 Group velocity extraction method for lightning electromagnetic wave propagation

Similar Documents

Publication Publication Date Title
JP2012189387A (en) Lightning discharge position orientation system
Chung et al. Optimizing earthquake early warning performance: ElarmS‐3
CN101581745B (en) Location method of partial discharge fault source of transformer
CN101179283B (en) Method of detecting and suppressing pulse noise
KR101930354B1 (en) Apparatus and method for detecting deception signal in global navigation satellite receiver
JP4817665B2 (en) Lightning location method and system
Betz et al. LINET–An international VLF/LF lightning detection network in Europe
CN104808124A (en) Acoustic-electric combined positioning method of switch cabinet local discharge fault positions
JP6067741B2 (en) Lightning strike charge amount estimation system, lightning strike charge amount estimation method, and program
US6768946B2 (en) Lighting strike position locating method, apparatus, system and program
CN109188359B (en) Lightning positioning method and device
CN107015064A (en) Lightning Location Method based on thunder and lightning multivariate data auto-correlation Shicha algorithm
CN103278704B (en) A kind of three-dimensional lightning detection system based on the Big Dipper and method
CN107037272A (en) Lightning Location Method based on thunder and lightning multivariate data peak-seeking Shicha algorithm
CN105223481A (en) Based on the partial-discharge ultrahigh-frequency signal initial time defining method of difference energy function
CN102967868A (en) Positioning device and signal processing method thereof
CN204188711U (en) Lightning multiband detection research station and system
CN103399030A (en) System and method for detecting phase difference of answering signals under air traffic control 3/A mode
JP4721945B2 (en) Crustal deformation detection device and crustal deformation detection method
JP2010281667A (en) Gps satellite signal quality monitoring method having gps satellite signal quality monitoring function and gps satellite signal quality monitoring apparatus having gps satellite signal quality monitoring function
Akos et al. Signal quality monitoring: test results
CN211826249U (en) Thunder and lightning positioner based on very low frequency electromagnetic method
RU2392746C1 (en) Method of detecting electronic devices illegally installed on site
De et al. Studies on the anomalies in the behaviour of transmitted subionospheric VLF electromagnetic signals and the changes in the fourth Schumann resonance mode as signatures of two pending earthquakes
Kitov et al. Detection, estimation of magnitude, and relative location of weak aftershocks using waveform cross-correlation: The earthquake of August 7, 2016, in the town of Mariupol

Legal Events

Date Code Title Description
A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20140513