Open access peer-reviewed chapter

Inverse Problem for an Electrical Dipole and the Lightning Location Passive Monitoring System

By Anatoly V. Panyukov

Submitted: July 1st 2017Reviewed: January 22nd 2018Published: February 22nd 2018

DOI: 10.5772/intechopen.74264

Downloaded: 221

Abstract

We solve the problem of the locating parameters, identifying equivalent dipole electromagnetic radiation source through measured horizontal magnetic and vertical electric components at some point of the infinite conducting ground. Methods based on analysis of measured signals are suggested. The problem under consideration, like any inverse problem of mathematical physics, is ill-conditioned. The consequences of this are the high sensitivity of the algorithm to the errors in the source data and calculation errors. All these circumstances do not allow to estimate the accuracy and reliability of the results obtained with the help of single-scale algorithms. The considered problem is contained in a complex of mathematical models of the practically important problem of forecasting the development of thunderstorm foci. Lightning meteorology focuses on investigating the lightning activities in different types of convective weather systems and the relationship of lightning to the dynamic and microphysical processes in thunderstorms. With the development and application of advanced lightning detection and location technologies, lightning meteorology has been developed into an important interdisciplinary between atmospheric electricity and meteorology. This paper reviews (1) methods to identify the dipole location and (2) possibilities to analyze the pre-radiation of thunderstorm clouds by the passive methods.

Keywords

  • algorithm
  • location finding
  • mathematical modeling
  • passive monitoring of the Earth’s electromagnetic field
  • electric dipole

1. Introduction

The problem of identifying the position parameters of an arbitrary-oriented electric dipole over a plane with infinite conductivity from its electromagnetic field induced at the observation point is considered. The considered problem is contained in a complex of mathematical models of the practically important problem of forecasting the development of thunderstorm foci. To estimate the location of the electromagnetic radiation source (EMR) based on the results of a single-point observation of the electromagnetic field induced by it, up to 1990, a number of devices were developed, based on the use of a vertical dipole as a model of an EMR source and physically realizable analog algorithms. In the next decade, the development of the use of an arbitrarily oriented dipole as a model of the EMR source and digital processing of observed signals was developed [1, 2, 3].

The effective algorithm for the single-point distance determination to a pulsed EMR source was proposed and investigated in the papers [4, 5], and pointed out to the emergence of irremovable uncertainty in the location of the EMR source by a one-point method caused by the difference in the orientation of the dipole from the vertical. The use of two or more observation points that do not belong to the same straight line makes it possible to determine the position parameters of the equivalent dipole [4, 6, 7].

A review of the status of passive storm monitoring systems by the end of 2003 and the demonstration of the use of lightning-position detection systems for the passive radar of hazardous meteorological phenomena are presented in [1, 2]. Modern methods for analyzing the field, allowing to determine the parameters of the source of EMP, characterizing its location and orientation are presented in the works [4, 5, 8] within the ISTC project #1822 was developed. As a result of the conducted field tests of this sample, during May–August 2004 more than 2.5 million atmospherics were recorded. Of these, not more than 10% were classified as radiation from a lightning discharge, the rest were classified as pre-threat radiation. Thus, the proportion of inter-cloud and intra-cloud discharges relative to the cloud-ground discharges turned out to be much higher than those noted for the work [9].

The registration of ominous cloud radiation (i.e., before the first lightning flash) by a single-point lightning protection system, for the purposes of forecasting thunderstorm development, was provided by the expansion of the dynamic range of receiving equipment and the further development of mathematical and software. At present, active radar facilities using a comparison of sounding and reflected signals are used to analyze the pre-threat state [9, 10, 11]. The obvious drawbacks of this approach are: (1) the high cost and the presence of the human factor; (2) the low probability of detecting individual discharges; (3) the lack of ecological compatibility due to the application of microwave radiation and its effect on the clouds. Passive methods of analysis of the pre-threat state, based on the analysis of the intrinsic radiation of clouds, exclude the noted shortcomings of active methods.

Previously, in passive thunderstorm monitoring systems, pre-threat radiation was either filtered out, or it was worked out incorrectly by single-point systems of thunderstorm location [12]. The database of ISTC project #1822 [8] field trials provides a large experimental material for testing the adequacy of thunderstorm models and testing of the developed software and devices. In general, the results of the project demonstrated the possibility and necessity of creating a new generation of thunderstorm detection systems and expanding the range of tasks they solve. This leads not only to a significant revision of the requirements for their technical characteristics, but also to the development of new mathematical models and algorithms for analyzing thunderstorm phenomena, their tracing and display, as well as archiving and their use by specialists from different subject areas.

2. Statement of the dipole location problem

2.1. Direct problem

Born and Wolf [13] give the expressions for the electromagnetic field of a dipole P0=ptn0δr0(moment p, direction n0, location r0, δ(*) is the Dirac delta function) in vacuum

EP0=14πε03ptr3+3ptcr2+ptc2rn0e0e014πε0ptr3+ptcr2+ptc2rn0,E1
HP0=14πptr2+ptrn0×e0E2

here e0is the direction from the dipole to the observer and the single and double primes denote time derivatives. For distances up to 300 km, a conducting plane surface is used to model the ground. The dipole in the half-space bounded by the infinitely conducting plane results in a field (E, H) consisting of the field (E(P0), H(P0)) of the source P0, and the reflected field; the latter can be represented by the field (E(P1), H(P1)) of an imaginary dipole P1=ptn1δr1, which is the mirror image of the dipole P1 therefore we have E=EP0+EP1,H=HP0+HP1(see Figure 1).

Figure 1.

Geometrical parameters of the mathematical problem.

Using the Cartesian system of coordinates with origin at the observation point О and the Oz axis being the normal of the bounding plane, we present the electric (Е) and magnetic (Н) components of the field in the coordinate form [4]:

Ext=Eyt=Hzt=0,Ezt=qt+qt3u2vα+qtuα2,E3
Hxt=sinφqt1α+qt1α2,Hyt=cosφqt1α+qt1α2,E4

here α=c/ris the inverse of the time of wave propagation from the source to the observer and the variables u,v,qtare determined from the equalities

wsinφ=sinθcosθ0sinψcosθsinθ0sinψ0,E5
wcosφ=sinθcosθ0cosψcosθsinθ0cosψ0,E6
u=sinθcosφψ,v=sinφψ0sinθsinψψ0,qt=wα2pt/2πc2r.E7

Other parameters are given in Figure 1.

2.2. Inverse problem

Components Ez, Hx, and Hycan be changed with the help of the antenna system consisting of a vertical electrical antenna and a pair of mutually orthogonal frame magnetic antennas, and hence used as initial data for the task of evaluating coordinates rψθfor the dipole position and ψ0θ0for the dipole orientation.

The inverse problem has a few specific aspects that differs it from the direct problem. If the source model is known, one often solves the inverse problem by searching the values of the parameters of the model, which give the best explanation of the measured data. Such a search is usually performed by an optimization algorithm which varies the values of the parameters until a necessary agreement between the measured and the model predicted fields is achieved. The search algorithms may get stuck in a local minimum and the solution will not be determined correctly. The use of a global optimization technique usually diminishes the problem of local minima but increase computational complexity.

For our case, it is possible to directly invert the values of the measured field into the values of the parameters describing the source location. In the present article, we use an explicit identification to infer the parameters of an electric dipole source situated over a perfectly conducting plane.

3. Calculation of the parameter φ. Exclusion of frames with interference of signals from different sources

Parameter φ, hereinafter called pseudo-bearing. For a vertical dipole (polar angle θ0=0we have φ=ψ,w=sinθ).For an inclined dipole, the pseudo-bearing φis differs from actual direction ψ, this difference depends on the orientation angles θ0and ψ0. For a horizontal dipole, we have φ=ψ0+πand w=cosθ, i.e., the pseudo-bearing is independent of the actual direction ψand is determined by only the dipole orientation ψ0. That is why the well-known direction finders [1, 2, 3] correctly record only cloud-ground discharges and give false estimates of bearing on inter-cloud and intra-cloud discharges. In fact, these methods evaluate pseudo-bearing φ.

To calculate the pseudo-bearing we use the following formulas:

cos2φ=hxxhyyhxxhyy2+4hxy2,sin2φ=hxxhyyhxxhyy2+4hxy2,E8
φ¯=atan1cos2φ1+cos2φsignhxy,E9

here hxx=T1T2Hxt2dt,hyy=T1T2Hyt2dt,hxy=T1T2HxtHytdt, [T1,T2] is signal observation interval.

Obviously, the signals observed in the antenna system can represent interference of signals from several sources. For the correct work of the algorithms described in the work, it is necessary to exclude frames with the presence of interference. Let us consider a possible way of detecting interference.

Interference leads to the appearance of a trend in the array of instantaneous error values

φ¯ti=Hxticosφ¯Hytisinφ¯,tiT1T2,i=0,1,2,,T2T1t=N.E10

Really let be Hxt=sinφ1S1t+sinφ2S2t+Nxt,Hyt=cosφ1S1t+cosφ2S2t+Nyt,Nx, Ny are uncorrelated white noise in the antenna channels having the same intensity N. Then φ¯t=S1tsinφ¯φ1+S2tsinφ¯φ2+NxtsinφNytcosφ.

It is easy to see that for φ1φ2function φ¯tfor any φ¯has a trend due to the difference in the functions S1t,S2tof the radiation sources. On the contrary, for φ¯=φ1=φ2, the function φ¯thas no such trend. Effect of interference make autocorrelation in function φ¯t. Thus, the fact of autocorrelation in the error function φ¯tis criterion for the presence of interference.

Durbin-Watson test is often used to determine the presence of autocorrelation. This criterion is usually used to establish the fact of the presence of an autocorrelation dependence of the first order in the error series, i.e., between its neighboring values φ¯tiand φ¯ti+1, ti,ti+1T1T2,i=0,1,2,,T2T1t=N,tis sampling step. Usually the neighboring error values are related by a stronger dependence than other values. Therefore, the absence of an autocorrelation dependence of the first order makes it possible to state quite confidently that there are no autocorrelation relationships in the errors.

Durbin-Watson statistics is calculated by the following formula:

d=i=0N1φ¯ti+1φ¯ti2i=0N1φ¯ti2.E11

We have 0d4. The values of d=0and d=4correspond to the cases when there is a strict positive or negative linear relationship between the shifted series φ¯tiand φ¯ti+1, respectively. If d=2, then the series is independent.

The method has been tested based on data generated as a result of three-month field testing of the lightning detection system in 2004 within the framework of the ISTC project #1822 [8] and containing 796,112 atmospheric spheres suitable for computation [8]. After calculating the value of the statistics dfor the entire set of signals from the signals (over the whole implementation), obtained because of the field tests of the lightning detection system and analysis of the results, the intermediate point values were chosen:

d1=0.4,d2=0.7,d3=3.3,d4=3.6.

If d0d1d44, then the hypothesis of the existence of first order autocorrelation in a series φ¯tiand φ¯ti+1is assumed. If dd2d3, then this hypothesis is rejected. If the value of the statistics dis in other intervals, then an unambiguous answer about the existence of autocorrelation cannot be given.

The analysis of the detected signals using the Durbin-Watson test showed the following results (see Table 1).

Value Range of Durbin-Watson statisticsNumber of signals
[0, 0.4]572,005
[0.4, 0.7]79,412
[0.7, 3.3]144,695
[3.3, 4]0

Table 1.

The relationship between the Durbin-Watson statistics interval and the number of corresponding signals.

Thus, we can confidently determine the coordinates of only 144,695 of the 796,112, based on the data of the entire implementation of each of the signals.

To eliminate this influence, you can select a part of the signal length (the implementation interval) at which interference appears less, or where there is none (there is no trend). Of course, the calculation of discharge characteristics over this implementation interval will be more accurate. In addition, the approach to isolating a “clean” signal can help in filtering the transients at the beginning of the signal reception.

4. Calculation of the parameters u,v,α

Further, we assume the parameter φ to be computed, and interference is absent in the recorded signals Ez, Hx, and Hy. Let us further use the function Ht=Hxtsinφ+Hytcosφ.

4.1. Primal algorithm

An analytical solution of the problem to identify parameters u,v,αis given in [4]

α=h1e2h2e1/h2e0h1e1,u=e0e2e12/h2e0h1e1,v=3h0uαe0αg2h0α,E12

here

ek=0EzktHktdt,hk=0Hkt2dt,g=0EztHtdt,k=0,1,2.E13

Calculated parameters u,v,α, φ can be used to determine the distance to the dipole r=c/α, as well as guaranteed interval estimates and probabilistic characteristics of angular coordinates of the dipole position.

Primal algorithm calculates parameters using formulas (3). It is assumed that the signals are represented by counts of instantaneous values with a given sampling rate. To calculate the coefficients (4), we use the known quadrature formulas.

An analysis of the sensitivity of the direct algorithm to the noise in the measured signals is given in [5]. It is shown that to achieve acceptable accuracy, pre-processing of input signals is necessary. It is found the dependence of mathematical expectation and variance of the distance determination error on the source spectrum and transfer function of the input filter. The estimates for the case corresponding to the averaged source of lightning discharges indicate the possibility of using the algorithm under study in systems for locating the lightning centers and allow us to develop the requirements for hardware of such systems.

4.2. Extremal algorithm

The primal algorithm was obtained in the assumption that the radiation source is the dipole, and the underlying surface—a plane with the infinite conductivity. Due to non-ideal models and noises in observed signals, we are forced to allow discrepancies in the resultant equations. Hence, intuition shows that most preferable is the least-squares principle. According to this principle, parameters of the model are determined from

Ix=T1T2x1Ezt+x2Ezt+x3Ht+x4Ht+x5Ht2dt=xAxmaxxE14
A=T1T2SSTdt,S=kEztEztkHtHtk1HtT,

k is the time constant for differentiating and integrating sections. This constant is introduced to make scales of coordinates for vector x consistent.

The extreme algorithm [14] for determining parameters u, ν, and α is calculated as follows:

Step 1. Let λ0λ1λ2λ34be eigenvalues of matrix А, xi,i=0,1,2,3,4are corresponding eigenvectors.

Step 2. Let mnbe the determinant of the matrix derived from the matrix x2x3x4through deleting lines having numbers m and n.

Step 3. For i=1,2do

βi=2252415+1i24+15241425if2514,24151i24+15241425214otherwise,

ui=βi13+2312,wi=45βi24+25
.

Step 4. Let i=argmini=1,21βiuiβiwiβi2wi2A1βiuiβiwiβi2wi2T,

Step 5. Let α=βik,u=ui,v=wi3u2.

A calculation experiment was performed to analyze accuracy and stability of the extreme algorithm, as well as to compare its characteristics with those of the primal algorithm. Software of the algorithm is implemented with the observation of agreements for the instrumental environment developed under the ISTC project 1822 [8]. Signals Ezt,Htare presented as integral counts of rapid values in discrete time moments, step of time discretization is t=2μs, number of quantification levels is 216. Input signals Ezt,Htare preliminarily processed by a smoothing filter. Usual difference and quadrature formulas with the accuracy of OΔt2were used for numerical differentiation and integration. The Givens rotation method is used to solve the complete symmetric problem of eigenvalues.

Extreme algorithm gives better accuracy if compared with the primal algorithms. Most significant difference in characteristics is observed for distances less than 50 km: error of the primal method sharply increases with the distance decrease, and the error of the extreme method is decreasing. Since with the decreasing distance, one has increased error in presentation of the real source as an electrical dipole, the experimental results demonstrate the extreme algorithm to be expedient at the distance r50km.

4.3. Parametrization of the algorithms

To solve the problem to identify the location parameters because of its poor condition is proposed to use a parameterized set of algorithms, and the final decision to accept the results of statistical analysis. Parameterizations of the primal and extremal methods are considered in paper [15].

The spectrum of dipole moment ptlies in a rather narrow frequency range, so the use of a bandpass filter with suitable lower and upper cutoff frequencies will maximize the use of all useful information that the signal carries and weaken the effect of noise present in the signal. Since the position parameters and the spectrum of the radiation source are not known a set of estimates E˜=αiuivi:i=12Nfor a family of bandpass filters with an amplitude-frequency characteristic

Fiω=1,ifωΩ̂iΩˇi,0,otherwise,i=12NE15

is constructed. The elements of the set E˜are considered as an implementation of a vector random variable, and for estimating the true position parameters u,v,and αto apply methods of statistical robust estimation.

It would be well to exploit more methods for increasing of statistical significance of result estimation. Put computational experiments’ results confirms efficiency the approach, but а wide variety of the measurement results do not lets to form up the well-provided measuring estimation of the dipole location parameters adequately. It is suggested that the filtration of the measurement results based on the cleaner functionals and combination of them for improving of the measuring estimation of the dipole location parameters. For construction of the cleaner functionals is the projection method [16]. This method is based on the projection of the inverse image of the vector-valued function EztHtonto the linear manifold Lof solutions of the system of differential Eqs. (1) and (2). There is a solution q(t) satisfying the system of differential Eqs. (1) and (2) under error-free values of the estimated parameters u,v,α. Therefore, the inverse image of the pair of signals EztHtbelongs to L, and the length of the projection is maximal. This allows you for filtering to use the length of the projection as a utility function. In terms of Fourier transforms, the utility function has the form [16]

Uuwα=Eωuω2+jωwα+wα2+Hωω2+jωαuω2+jωwα+wα22+ω2+jωα2.E16

Function Uuwαis continuous for uwα0. Algorithm for calculating Uuwαis completely stable for signals represented of 4096 readouts with the sampling time 2–10−6 s, i.e., containing 2047 harmonics with the first-harmonic discretization frequency 122 Hz.

Set EU0=uwαE˜:UuwαU0}contains most reliable estimates. Statistical analysis of the set of estimates makes it possible to estimate uvαprobability density.

4.4. Spectral statistical method for location parameters identifying of a dipole electromagnetic radiation source

The spectral statistical method makes it possible to obtain more stable solutions at a lower computation cost compared with the previously developed parametric extremum method. The spectral statistical method algorithm can be naturally parallelized. The proposed method based on the analysis of the measured-signal spectra, allows one to get many estimates of the source location, choose the final estimate of the results of analysis of the entire totality of these estimates, and therefore, reach stability in determining the source location. To reduce the computational complexity, it is preferable to analyze individual harmonics rather than the entire frequency band. In this case, the algorithms for estimating the dipole location are substantially simplified [17].

As before, we will consider signals represented as a fast Fourier transformation Eω,Hωcontaining 2047 harmonics with the first-harmonic discretization frequency 122 Hz. Let us divide the set Ω=ω1ω2ω2047of all harmonics into subsets

Ω1=ω1ω2ω131,Ω2=ω132ω133ω400andΩ3=ω401ω402ω2047.

Spectral statistical algorithm to find estimates û,v̂,α̂of the parameters u,v,αand their variance σû2,σv̂2,σα̂2is calculated as follows:

Step 1. For each ωΩcalculate

fω=ReEωHω¯,gω=ImEωHω¯,hω=Hω2.

Step 2. Calculate the noise intensity

σN2=ωΩ3Hω2Ω3.

Step 3. For each ωΩ2calculate

uω̂=1+σN2hωfωhω.

Step 4. Calculate

û=ωΩ2hωuω̂ωΩ2hω,σû2=ωΩ2fωhωûωΩ2hω.

Step 5. For each ωΩ1calculate

Aω=ωhωûfω1,Sω=σN2hω,αω̂=A2ωSωω2Sω+0.5+SωA2ωω2+0.25.

Step 6. Calculate

α̂=ωΩ1hωω2Aωαω̂ωΩ1hωω2Aω,σα̂2=ωΩ1hωω2Aωα̂αω̂2ωΩ1hωω2Aω.

Step 7. For each ωΩ1calculate

Wω=fωAωωgωhωAω.

Step 8. Calculate

ŵ=ωΩ1hωWωωΩ1hω,σŵ2=ωΩ1hωWωv̂2ωΩ1hω.

The rationale for the effectiveness of the algorithm and computational experiments are given in [17].

5. Secondary processing of electromagnetic field monitoring results

The most important problem in work of the single-point systems of storm location is ineradicable error in the bearing determination of each source of radiation taken separately. Amount of this error depends on orientation of the lightning equivalent dipole, and the “pseudo-bearing” determined by the device for horizontal dipoles does not depend from the real bearing at all and is defined only by projection of orientation of the dipoles to the ground plane.

The given problem can be solved only by the analysis of the information from all set of discharges of the storm. At International conference ETC’2006 [18], it was suggested to display the stream of signals as a map of lightning discharge hit probability density in this or that point of terrestrial surface. The representation proposed increases probability of definition of true thunderstorm discharges location, and hence, raises accuracy of their fixing.

Let us describe the technique of this approach realization, the difficulties arising, and the variants of their decision in detail.

Let values (r,ψ,θ)be the polar coordinates of the lightning equivalent dipole for the measuring device coordinate system. As it is known [4, 18], Cartesian coordinates x,yof the dipole meet the condition

xcosφ+ysinφ=rsinθcosφψ=ru.E17

This equation includes known parameters: ris distance from the dipole to the observation point, u,vare the variables of the source dipole model (3)—and unknown Cartesian coordinates x,yof the dipole location that is a feature of an uncertainty. The coordinates of the points of possible location of a dipole most remote from origin are the solutions of a set of equations

xcosφ+ysinφ=ru,x2+y2=r2.E18

They are equal

x1,2=rucosφ1u2sinφ,y1,2=rusinφ±1u2cosφ.E19

The coordinates x,yof all probable points of a source location introduce on a plane XOYa section of a straight line (10) between points x1y1and x2y2. The set of probable location of a dipole have one-parameter representation

Lruv=xλ=rucosφ+r2λ11u2sinφyλ=rusinφ+r12λ1u2cosφ,zλ=r2x2λy2λ=2rλ1λ1u20λ1.E20

A measure of the coordinates set Lruvequal

L=01dxλ2+dyλ2+dzλ2=4.844r1u2.

The value of L represents a measure of the uncertainty in the estimation of the Cartesian coordinates of the dipole location. Uncertainty is absent when L = 0. It is reached for a vertical dipole (u = ±1). It is possible value u = 0 for a horizontal dipole when L reaches the maximum value.

A more practical measure of uncertainty is probability

Pijk{[(xλ,yλ,zλ]Cijk}

of the dipole location in space cell

Cijk=xyzxi+1Δyj+1Δzk+1Δ,

here Δis cell size.

It is obvious that Pijk>0if CijkLruv. To find cells in which the dipole is located with a nonzero probability, it is sufficient to find the minimal covering CXOYruvof the projection of Lruvonto the plane XOY by flat cells Cij:

minrucosφarccosuΔrucosφ+arccosuΔimaxrucosφarccosuΔrucosφ+arccosuΔ,
minrucosφarcsinuΔrucosφ+arcsinuΔjmaxrucosφarcsinuΔrucosφ+arcsinuΔ.

An example of coverage CXOYruvis shown in the Figure 2.

Figure 2.

An example of coverage C XOY r u v .

Cell CijCXOYruvin accordance with equality (12) uniquely defines a space cell Cijk,k=rΔi2j2with a nonzero probability of the dipole location.

To find this probability, we accept a completely plausible hypothesis about the uniform distribution of the probabilities of all possible values of the angle ψ0(see Figure 1). In this case the probability density for angle ψis [4]

fψη=s2ψ̂s2cos2φη+sin2φη,arccosuηφarccosuE21

here

s=1uvuv,ψ̂=arctan1stanarccosu.

Thus, for CijCXOYruvone can take

Pijk=fψη,k=rΔi2j2,η=arctanji.

Let RTbe the set of registered signals in the time interval T.The integral estimate grade of radiation locus membership to cell Cijkover period Tlet be

gradeijk=dRTPijkdmaxlmndRTPlmnd.

It is suggested to display the stream of signals as a grademap. The representation proposed increases in probability of definition of true thunderstorm discharges location, and hence, raises accuracy of their fixing.

6. Using a network of sensors

Currently, mainly for monitoring thunderstorm activity, multi-point systems based on the application of the difference-ranging method and its variations to the results of monitoring the Earth’s electromagnetic field in the VLF and VHF bands are used. Some lacks such an approach are: (1) insufficient reliability, which results from the need communication networks for the exchange of information between observation stations spaced over significant distances, and the need for time synchronization at the sub-microsecond level; (2) the inability to analyze the cloud radiation before the thunderstorm. Previously, this was the reason for the development of single-point systems.

Now, the rapid development of computer technologies and communication systems have made it possible to combine single-point lightning location finders into a system and to organize joint processing of signals from individual lightning range finder. Automating the collection of information from observation points helps different services to respond more exactly to changes in thunderstorm conditions.

The uncertainty in the dipole location is fundamentally unavoidable when single-point systems are using, however these errors can be eliminated by determining its generalized parameters at two or more points [4, 6, 7]. Thus, by integrating single-point systems of passive monitoring of thunderstorm activity into a single computer network, it is possible to increase the probability of detecting a lightning discharge and the accuracy of determining its coordinates. In addition, an increase in the degree of automation of information collection from observation posts makes it possible to receive data promptly and make adequate decisions in accordance with the current thunderstorm situation.

6.1. Direction-range method

Let xyzbe the Cartesian coordinates of the dipole location point. Let xiyibe Cartesian coordinates of i-th observation post, iJ.Let riφiuibe the dipole location parameters of i-th observation post. Then the system of equations

xixcosφi+yiysinφi=riui,iJE22

holds. If the observation points do not lie on one line, then system (22) has full rank and can be solved by the least squares method

xy=1ΔiJaikJsinφksinφkφi1ΔiJaikJcosφksinφkφi,E23
Δ=1k<lJsin2φkφi,ai=xicosφi+yisinφi,iJ.

The method uses for its calculations the values ri=cαiuiφimeasured by i-th autonomous lightning detection finder, and not the bearing estimation ψ, as a classical direction-finding method. This makes it insensitive to the presence of an anomalous component in the observed magnetic field. There is no requirement for high-precision synchronization for the direction and distance method, which is a characteristic of the difference-range method, and the possibility of estimating the position of the lightning discharge at each point eliminates collisions in identifying the correspondence of the recorded signals to specific lightning discharges. Finally, this method makes it possible to determine the three coordinates of the location of an equivalent dipole source.

The development of this method, which makes it possible to find the angular coordinates θ0ψ0of the EMR dipole source, is given in [6, 7].

It should be noted that the use of the direction-ranging method involves knowledge of the ratio of the effective heights of magnetic and electric antennas and high requirements for the accuracy of the orientation of magnetic antennas. The determination of the ratio of effective antenna heights is a complex technical problem requiring a reference meter or a long time for the collection and processing of statistical data. Therefore, it is of interest to develop other methods for determining the location of a lightning discharge, which do not require knowledge of this relationship.

6.2. Range method

Range method for determining the coordinates of a dipole EMR source [19], in contrast to a direction- range finder, does not require a prior knowledge of the ratio of the effective antenna heights, and high accuracy of the orientation of magnetic antennas. The essence of the method is that the Cartesian coordinates xiyiof the i-th point of observation, iJ. Cartesian coordinates xyzof the dipole EMR source and ranges rifrom the i-th observation points up to the dipole EMR source satisfy the equations system

ri2=xxi2+yyi2+z2,iJE24

reducible by the introduction of the variable

m=xi2+yi2+zi2C

to a system of linear equations

mC2xix2yiy=ri2xi2yi2,iJ,E25

here constant Cis the scaling factor for reducing the condition number of the equation system.

If the observation points do not lie on one line then the resulting system will have full rank, and its least-squares method, the solution is:

mxy=ATA1ATB,whereA=C2x12y1C2x22y2C2xJ2yJ,B=r12x12y22r22x22y22rJ2xJ2yJ2.E26

At fixed coordinates of all observation points, the matrix A is constant. Therefore, even at the design stage of the system it is easy to construct an optimal algorithm for calculating unknowns, since ATA1ATcan be counted once for the implementation of the system.

For example, let observation points be

x1y1=RO,x2y2=0R,x3y3=RO,x4y4=0R.

Then minimum value of the condition number equal 1 and is reached for C=2R2.The estimates of the dipole location coordinates are the following:

x=r12r322R,y=r22r422R.E27

It is easy to determine the coordinates both the dipole orientation and its location with Cartesian coordinates (27) and the distances rito the observation pointsiJ.

7. Conclusion

New generation of thunderstorm passive monitoring systems expand the range of problems they solve: development of new mathematical models and algorithms for analyzing thunderstorm phenomena, their tracing and display [20].

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Anatoly V. Panyukov (February 22nd 2018). Inverse Problem for an Electrical Dipole and the Lightning Location Passive Monitoring System, Electric Field, Mohsen Sheikholeslami Kandelousi, IntechOpen, DOI: 10.5772/intechopen.74264. Available from:

chapter statistics

221total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Maxwell-Fredholm Equations

By Juan Manuel Velázquez-Arcos, Alejandro Pérez-Ricardez, Ricardo Teodoro Páez-Hernández, Samuel Alcántara Montes and Jaime Granados-Samaniego

Related Book

First chapter

Introductory Chapter: Nano-Enhanced Phase-Change Material

By Mohsen Sheikholeslami Kandelousi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us