CN110297284A - A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument - Google Patents
A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument Download PDFInfo
- Publication number
- CN110297284A CN110297284A CN201910692270.6A CN201910692270A CN110297284A CN 110297284 A CN110297284 A CN 110297284A CN 201910692270 A CN201910692270 A CN 201910692270A CN 110297284 A CN110297284 A CN 110297284A
- Authority
- CN
- China
- Prior art keywords
- electric field
- thunder cloud
- cloud
- thunder
- point
- 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
- 230000005684 electric field Effects 0.000 title claims abstract description 114
- 230000033001 locomotion Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 claims description 19
- BULVZWIRKLYCBC-UHFFFAOYSA-N phorate Chemical compound CCOP(=S)(OCC)SCSCC BULVZWIRKLYCBC-UHFFFAOYSA-N 0.000 claims description 6
- 230000001939 inductive effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 30
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 description 17
- 238000002474 experimental method Methods 0.000 description 8
- 230000004807 localization Effects 0.000 description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/12—Measuring electrostatic fields or voltage-potential
- G01R29/14—Measuring field distribution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/16—Measuring atmospheric potential differences, e.g. due to electrical charges in clouds
Landscapes
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument is based on model firstly, establishing the spatial model of tracking thunder cloud movement routine according to mirror image law theory, defines the dynamic parameter in thunder cloud orientation;Then, the Potential distribution formula using thunder cloud at three-dimensional atmospheric electric field instrument, obtains three-dimensional electric field component;Finally, orientation where thunder cloud is associated with time point, the tracking of thunder cloud movement routine is realized.Electric field component measurement error, azimuth are provided, the relationship of pitch angle and path trace accuracy analyzes tracking performance.The present invention can accurately track the movement routine of thunder cloud, have preferable tracking effect.
Description
Technical field
The invention belongs to thunder cloud forecasting and warning fields, and in particular to a kind of thunder cloud shifting based on three-dimensional atmospheric electric field instrument
Dynamic path following method.
Background technique
Thunder cloud is the major sources of lightning, and lightning is a kind of great natural calamity, and Lightning Disaster is joined already
It closes state and is classified as " one of ten kinds of natural calamities of most serious " in relation to part.Thunder and lightning is frequently accompanied by strong heavy rain and fitful wind, sometimes
It can also be with cyclone and hail.Estimate according to relevant department, all over the world per minute about 1800 lightning accompanied by peals of thunder and carrying out
In.Their each seconds about issue 600 lightnings, wherein there is 100 attack earth.Therefore, effective thunder cloud monitors hand
Section is particularly important to lightning protection mitigation ability is improved.
Now, Sensitive Units layout, the scientific research personnel such as atmospheric electric field detector is widely used in aerospace, oil field is stored in a warehouse also accumulate
It joins in the research work of atmospheric electric field measurement pole.The study found that atmospheric electric field is made of vertical component and horizontal component
's.Row letter man of virtue and ability et al. detects space electric field using three-dimensional electric field sensor and dielectric constant sensor, proposes one
Kind is suitable for the three-dimensional electric field system of atmospheric electric field measurement.Xuwei et al. develops a kind of three-dimensional atmospheric electric field instrument, for simultaneously
The horizontal and vertical component of atmospheric electric field is measured, and introduces particle swarm optimization algorithm inverting thunder cloud orientation.
However, these are studied while realizing the measurement of atmospheric electric field three-dimensional component, it can not be to the real-time movement of thunder cloud
Path is tracked, and thunder cloud forecasting and warning truly is not achieved, poor, the mobile road for thunder cloud positioning real-time
Diameter tracking cannot achieve problem.
Summary of the invention
Goal of the invention: the present invention proposes a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument, energy
Solving the problems, such as thunder cloud positioning, real-time is poor, movement routine tracking cannot achieve.
Technical solution: a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument of the present invention,
The following steps are included:
(1) according to mirror image law theory, the spatial model of tracking thunder cloud movement routine is established, and defines thunder cloud orientation
Dynamic parameter;
(2) Potential distribution formula of the thunder cloud at three-dimensional atmospheric electric field instrument is utilized, three-dimensional electric field component is obtained, acquires thunder
Sudden and violent cloud point STThe spherical coordinates of T variation at any time;
(3) the set P of the thunder cloud dynamic coordinate of T variation at any time is defined, and by orientation where thunder cloud and time point
It is associated, realizes the tracking of thunder cloud movement routine.
Further, the realization process of the step (1) is as follows:
Using N point as coordinate origin, three-dimensional cartesian coordinate system is established, in which: N (0,0,0) is three-dimensional atmospheric electric field instrument place
Position, ST(xT,yT,zT) it is thunderstorm cloud point STT and the dynamic position that changes, S at any timeT'(xT,yT, 0) and it is STIn x0y plane
Projection;The sum of the height above sea level of h expression atmospheric electric field detector oneself height and its position;Thunderstorm cloud point STDynamic level drift angle
It is respectively α with pitch angleT, βT, and RTFor thunderstorm cloud point STTo the real-time range of electric field instrument N;At electric field instrument N, thunder cloud is measured
Point STReal-time electric field strength be ET。
Further, the step (2) the following steps are included:
(21) according to charge inducing in image method and image charge phase scheduling theory, regard thunder cloud region as charge
Q obtains thunderstorm cloud point STReal-time Potential distribution at three-dimensional atmospheric electric field instrument NAre as follows:
Wherein, q' is the image charge of point charge q, ε1、ε2Ground respectively where air and three-dimensional atmospheric electric field instrument N
Dielectric constant;
(22) due to On The Electric Field of Thunderstorm intensity E that atmospheric electric field detector N point measuresT, it is a kind of trivector, it can be to ETIt carries out
Orthogonal Decomposition obtains:
ET=EX,T+EY,T+EZ,T
Wherein, EX,T、EY,T、EZ,TRespectively electric field instrument N measures x-axis, y-axis, the real-time electric-field strength of thunder cloud on z-axis direction
Degree, and it is orthogonal two-by-two;
(23) thunderstorm cloud point S is soughtTSpherical coordinates (the R of T variation at any timeT',αT',βT') are as follows:
Wherein, RT' it is thunder cloud STTo the distance R of electric field instrument NTInverting value, αT', βT' it is respectively thunder cloud dynamic level
Drift angle αT, pitch angle βTInverting value;As quantity of electric charge q, permittivity ε1, ε2When known, real-time thunder cloud coordinate (RT',αT',
βT') only with three-dimensional electric field vector EX,T、EY,T、EZ,TIt is related.
Further, the set P of the dynamic coordinate that T changes at any time of thunder cloud described in step (3) are as follows:
Wherein, T1, T2To TnRespectively refer to for the 1st, 2...n time point, and with time point T1, T2...TnCorresponding thunder cloud
Coordinate is respectively
The utility model has the advantages that compared with prior art, beneficial effects of the present invention: the tracking can accurately track thunder cloud
Movement routine, make it possible orientation where real-time detection thunder cloud, there is preferable tracking effect.
Detailed description of the invention
Fig. 1 is flow chart of the invention;
Fig. 2 is that thunder cloud movement routine tracks spatial model;
Fig. 3 is thunder cloud to atmospheric electric field detector distance, the relationship of electric field measurement error and tracking range error;
Fig. 4 is thunder cloud to atmospheric electric field detector distance, the relationship of pitch angle and level angle measurement error;
Fig. 5 is thunder cloud to atmospheric electric field detector distance, the relationship of pitch angle and elevation measurement error;
Fig. 6 is real-time thunder cloud movement routine tracing figure on April 9th, 2019;
Fig. 7 is the Nanjing 10:06 radar station radar echo map on April 9th, 2019;
Fig. 8 is the Nanjing 10:33 radar station radar echo map on April 9th, 2019;
Fig. 9 is the Nanjing 10:50 radar station radar echo map on April 9th, 2019;
Figure 10 is the Nanjing 11:07 radar station radar echo map on April 9th, 2019;
Figure 11 is real-time thunder cloud movement routine tracing figure on April 22nd, 2019;
Figure 12 is the Nanjing 20:43 radar station radar echo map on April 22nd, 2019;
Figure 13 is the Nanjing 20:54 radar station radar echo map on April 22nd, 2019;
Figure 14 is the Nanjing 21:05 radar station radar echo map on April 22nd, 2019;
Figure 15 is the Nanjing 21:16 radar station radar echo map on April 22nd, 2019.
Specific embodiment
In the following with reference to the drawings and specific embodiments, the present invention is furture elucidated.
As shown in Figure 1, the present invention provides a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument, packet
Include following steps:
1, thunder cloud movement routine trace model is established, specifically:
As shown in Fig. 2, establishing three-dimensional cartesian coordinate system using N point as coordinate origin, in which: N (0,0,0) is three-dimensional atmospheric
Electric field instrument position, ST(xT,yT,zT) it is thunderstorm cloud point STT and the dynamic position that changes, S at any timeT'(xT,yT, 0) and it is ST
Projection in x0y plane;The sum of the height above sea level of h expression atmospheric electric field detector oneself height and its position;Thunderstorm cloud point ST's
Dynamic level drift angle and pitch angle are respectively αT, βT, and RTFor thunderstorm cloud point STTo the real-time range of electric field instrument N;In electric field instrument N
Place, measures thunderstorm cloud point STReal-time electric field strength be ET。
2, according to charge inducing in image method and image charge phase scheduling theory, regard thunder cloud region as a charge q,
Obtain thunderstorm cloud point STReal-time Potential distribution at three-dimensional atmospheric electric field instrument NAre as follows:
In formula (1), q' is the image charge of point charge q, ε1、ε2Ground respectively where air and three-dimensional atmospheric electric field instrument N
Dielectric constant.
The On The Electric Field of Thunderstorm intensity E measured due to atmospheric electric field detector N pointT, it is a kind of trivector, it can be to ETIt carries out orthogonal
It decomposes, obtains:
ET=EX,T+EY,T+EZ,T (2)
In formula (2), EX,T、EY,T、EZ,TRespectively electric field instrument N measures x-axis, y-axis, the real-time electric field of thunder cloud on z-axis direction
Intensity, and it is orthogonal two-by-two.
Respectively to x, y, the Potential distribution in z-axis directionDerivation is carried out, E can be obtainedX,T、EY,T、EZ,T:
In general, thunder cloud coordinate parameters zTBe higher by 2 orders of magnitude than height h, then:
zT≈zT-h≈zT+h (4)
Vector correlation according to Fig. 1, thunder cloud STTo the distance R of three-dimensional atmospheric electric field instrument NTAre as follows:
Using formula (4), (5), formula (3) is become:
In formula (6),
Using formula (6), thunderstorm cloud point S is obtainedTSpherical coordinates (the R of T variation at any timeT',αT',βT') are as follows:
In formula (7), RT' it is thunder cloud STTo the distance R of electric field instrument NTInverting value, αT', βT' it is respectively thunder cloud dynamic
Level angle αT, pitch angle βTInverting value.As quantity of electric charge q, permittivity ε1, ε2When known, real-time thunder cloud coordinate (RT',
αT',βT') only with three-dimensional electric field vector EX,T、EY,T、EZ,TIt is related.E has been determinedX,T、EY,T、EZ,TAfterwards, thunderstorm can further be detected
Cloud movement routine.
3, the set P of the thunder cloud dynamic coordinate that T changes at any time is defined are as follows:
In formula (8), T1, T2To TnRespectively refer to for the 1st, 2...n time point, and with time point T1, T2...TnCorresponding thunderstorm
Cloud coordinate is respectively
Orientation where thunder cloud is associated with time point, thunder cloud movement routine can be realized by convolution (7), (8)
Tracking.
According to above-described embodiment method, below to the thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument from
Performance evaluation is carried out in terms of ranging direction finding:
Thunder cloud movement routine tracking performance evaluation
According to aerial charge field distribution and thunder cloud charge structure principle, dielectric constant of air ε is taken1It is 1, atmospheric electric field
The permittivity ε on ground where instrument N2It is 5, quantity of electric charge q is 5C.Separately setting atmospheric electric field measurement component standard deviation isWith
Study the overall performance of thunder cloud movement routine tracking.
Thunder cloud coordinate formula is the basis of thunder cloud movement routine tracking.It, can according to indirect error theory
It is obtained using formula (7) by electric field component measurement errorCause distance RT', level angle αT', pitch angle βT' measurement error
σR, σα, σβAre as follows:
Formula (9) are further simplified, are obtained:
By formula (10) it is found that the range error σ of thunder cloud movement routine method for tracingROnly with distance RT' and electric field measurement mistake
DifferenceIt is related, and with thunder cloud level angle αT' and pitch angle βT' unrelated;Its angle measurement error σα, σβWith distance RT', pitch angle
βT' and electric field measurement errorIt is related, but both with level angle αT' unrelated.Therefore, in practical applications, the track side
The performance of method is not influenced by level angle variation;In the maximum measurement range of three-dimensional atmospheric electric field instrument, distance RT' and electric field
Measurement errorSmaller, the ranging direction finding precision of this method can be higher.
The analysis of thunder cloud movement routine tracking range performance
Using formula (10), distance R is studiedT', electric field measurement errorWith thunder cloud movement routine tracking range error
σRRelationship, simulation result is as shown in Figure 3.In Fig. 3, the range error σ of thunder cloud movement routine trackingRWith distance
RT' and electric field measurement errorIncrease and increase.As distance RT' at 0 to 0.55km, range error σRIt is missed by electric field measurement
DifferenceInfluence smaller, σR maxOnly reach 1m;And as distance RT' be greater than 0.55km when, range error σRWith errorIncrease
And rise sharply, σR maxReach 6.4m.Similarly, when electric field measurement errorWhen less than 0.15kV/m, range error σRHardly by
Distance RT' influence, σR maxFor 1m;And work as errorAt 0.15 to 1kV/m, with distance RT' increase, error σRRise
To 6.4m.To sum up, although the range error σ of this methodRUnavoidably, but still can receive range in, show preferably with
Track range performance.
The analysis of thunder cloud movement routine tracking Measure direction performance
Using formula (10), distance R is studiedT', thunder cloud pitch angle βT' and its level angle measurement error σαRelationship, imitate
True result is as shown in Figure 4.In Fig. 4, with distance RT' and pitch angle βT' increase, the water of thunder cloud movement routine tracking
Flat measurement of angle error σαAlso it increases with it.When distance r' is 0 to 0.42km, measurement error σαBy pitch angle βT' variation shadow
Ring smaller, error σαLess than 0.2 degree;And when distance r' is greater than 0.42km, with pitch angle βT' increase, measurement error σαIt is in
Exponential type trend rises, σαmaxReach 1.07 degree.As distance RT' and pitch angle βT' it is larger when, measurement error σαAlso larger.To the greatest extent
Pipe is in this way, more often than not, error σαHardly influenced by the two.
Formula (10) are utilized again, study distance RT', thunder cloud pitch angle βT' and its elevation measurement error σβRelationship,
Simulation result is as shown in Figure 5.In Fig. 5, the elevation measurement error σ of thunder cloud movement routine trackingβWith distance RT' and
Pitch angle βT' increase and increase.As distance RT' at 0 to 0.35km, elevation measurement error σβBy pitch angle βT' influence compared with
It is small, σβmaxReach 0.005 degree;Similarly, as pitch angle βT' less than 30 degree when, error σβHardly by distance RT' influence,
σβmaxReach 0.01 degree.And as distance RT' and pitch angle βT' it is larger when, measurement error σβAlso 0.0375 is only slowly risen to
Degree.
In actual measurement experiment, the three-dimensional atmospheric electric field instrument N at Nanjing Information engineering Univ (NUIST) test point is pacified
Mounted in Electronics and Information Engineering institute, atmospheric electric field detector is placed on samely with the microphone array for being used for thunder auditory localization
Point.Herein, about 28 meters of atmospheric electric field detector range averaging sea level, three-dimensional system of coordinate X-axis positive axis where defining electric field instrument are south,
Y-axis positive axis is east.Three-dimensional electric field component value (the E recorded using host computer interfaceX,T,EY,T,EZ,T) (units:kV/m), it obtains
To by thunderstorm cloud point charge coordinate (RT',αT',βT') the set P (units:km, degrees, degrees) of composition, and and thunder
Auditory localization, radar return diagram data are compared analysis.
For the measurement distance of three-dimensional atmospheric electric field instrument 10 between 20km, field minimum resolution ratio is 10V/m, data sampling
Rate is 250 times per second, and can guarantee three-dimensional electric field component mean data of output per minute.
Rainy weather experiment:
Assigned at 11 points 07 minute on 04 09th, 2,019 10: 06, the electric field component data that three-dimensional atmospheric electric field instrument measures
(unit: kV/m) is as shown in table 1.
Table is tested on April 9th, 1 2019 for the first time
In table 1, relative to EX,TAnd EY,T, EZ,TIt is larger, by data shown in observation table 1, although can primarily determine 10
Point 06, which assigns to the 11 points of station 07 minute this period NUIST upper areas, may lightening activity, but because lacking electric field component and its
The association of his physical quantity, can not detect thunder cloud movement routine, can not achieve thunder cloud forecasting and warning eventually.Therefore, the number of table 1 is utilized
According to introducing thunder cloud movement routine tracking passes through correlation time point TnWith orientation (R where thunder cloudT',αT',βT'), into
One step obtains real-time tracking figure, as shown in Figure 6.10: 06 assign to 11: 07 by stages between thunder cloud movement routine be 1-2-3-4,
It is 1'-2'-3'-4' that it, which is projected to the path in X0Y plane,.As can be seen from Figure 6,10 points 06 minute, three-dimensional atmospheric electric field instrument measures
Thunder cloud orientation be (0.558, -55.71,72.61), thunderstorm cloud point about at 55.71 degree of directions of south by west, 0.558 km,
At this point, pitch angle is larger, 72.61 degree are reached;Assign at 10 points 33 minutes from 10: 06, thunderstorm cloud point substantially from southwest southeastwards
57.65 degree of direction movements, and it is in the overhead away from 0.576 km of electric field instrument;Assign at 10 points 50 minutes from 10: 33, thunderstorm cloud point is several
It does not move, rests on the overhead of 0.591 km of NUIST website southeastern direction;10 points after 50 minutes soon, thunder cloud is in 11
Point 56.02 degree of directions movements northeastward in 07 minute, thunderstorm cloud point is away from about 0.618 km of electric field instrument, at this point, pitch angle is down to 51.04
Degree, illustrates thunder cloud just far from test point.To sum up, this experimental period is in spring thunder and takes place frequently season, in conjunction with Fig. 6, stands in NUIST
Empty region on atmospheric electric field detector has reproduced thunder cloud from the process having to nothing using this method, has preferably demonstrated having for method
Effect property.
In addition, setting array element spacing is 0.8 meter using the microphone array put up, thunder is carried out within the same period
Auditory localization is compared with thunder cloud movement routine tracking.The auditory localization algorithm being introduced into document [1], by catching
Thunder source signal is caught, after screening out the smaller signal of amplitude and data closer from array, obtains result as shown in table 2.
Table thunder source positioning experiment on April 9th, 2 2019
In table 2, it is (0.652,52.14,79.38) that microphone array, which is listed in the sound source coordinate that measure for 42 minutes at 10 points,.The sound source
In 52.14 degree of directions of south by east, from 0.652 km of array center.Especially its pitch angle is 79.38 degree, is approximately perpendicular to Z
Axis.Greatly there may be cross thunder for empty region on the microphone array of the station NUIST for this explanation.In addition, assigning at 10 points 10: 33
50/, the thunder cloud level angle and pitch angle that tracking measures respectively may be about 58 degree, 63 degree, and both shown in table 2
Respectively 52.14,79.38 degree, which has better accuracy relative to existing tracking.
For the validity for further proving the tracking, official's data comparison can be measured with Nanjing radar station.Choose with
1 time point of table identical radar echo map, as is seen in figs 7-10.In addition, the cross center marked in radar map is the station NUIST
Three-dimensional atmospheric electric field instrument position.Find out from Fig. 7-10, is assigned to 11:07 minutes from 10: 06, empty region in the test point of the station NUIST
Radar echo intensity be not less than 40dBZ, highest has been even up to 55dBZ, this shows to accumulate in cloud layer above the test point
A large amount of charges, the region truly have thunder cloud to exist.It is found by analysis, from 10 points of 06 minute thunder clouds from west close to test point
Start, radar echo intensity is larger always, which is completely covered by test point overhead;Then, thunder cloud is from southwest
Yearn for east movement, is continued for being moved at 11 points 07 minute.After thunder cloud region is moved far from test point toward northeastward, radar
Echo strength is just declined.Fig. 6 to 10 is made a general survey of, in rainy weather experiment, thunderstorm region is rapid from test point southwestward
Mobile to its northeastward, the movement of thunder cloud during this period is consistent with the result that tracking measures.
Fair weather experiment:
Assigned at 21 points 16 minutes on 04 22nd, 2,019 20: 43, the electric field component data that three-dimensional atmospheric electric field instrument measures
(unit: kV/m) is as shown in table 3.
Table second of experiment on April 22nd, 3 2019
In table 3, due to electric field component EX,T, EY,T, EZ,TIt is smaller, it can tentatively judge to assign at 21 points 16 minutes 20: 43
Period, NUIST test point upper area do not have lightening activity.Similarly, if introducing thunder cloud movement routine tracking, just
The real-time movement routine that thunder cloud can further be detected, is more intuitive to see the motion process of thunder cloud.It is obtained using 3 data of table
To the measured result of tracking, as shown in figure 11.In Figure 11,20: 43 assign to 21: 16 by stages between thunder cloud movement routine
For 1-2-3-4, projection to the path in X0Y plane is 1'-2'-3'-4'.As can be seen from Figure 11,20 points 43 minutes, three-dimensional atmospheric
It is (3.141,143.13,10.20) that electric field instrument, which measures thunder cloud orientation, and thunderstorm cloud point is about in 36.87 degree of north by east direction.This
When, thunder cloud farther out from electric field instrument, has reached 3.141 kms, but pitch angle is only 10.20 degree.This shows that this is engraved in the station NUIST
Upper empty region does not have thunder cloud presence.Assign at 20 points 54 minutes from 20: 43, thunder cloud is still in away from 3.557 km of electric field instrument
Northeastward, it is just gradually far from electric field instrument.Assign at 21 points 05 minute from 20: 54, thunder cloud is past eastern from 39.81 degree of northeast direction
Southern 51.34 degree of directions fast move, and rest on the overhead away from 4.042 km of NUIST website.At 21 points after 05 minute, thunder cloud
It is mobile to continue southeastwards direction, assigns in 21: 16 up to 45 degree of directions of south by east, the position away from 4.552 km of electric field instrument.Entirely
Process, farther out away from atmospheric electric field detector, and pitch angle is smaller for thunder cloud.This can illustrate that the station NUIST test point overhead has no
Thunder cloud.In addition, this method traces into thunder cloud in the frequent activities of the east region away from 4 km of electric field instrument or so, it can be pre-
Surveying thunder cloud can continue to move to later towards southeastern direction.
Microphone array is utilized again, thunder auditory localization is carried out within the same period, the results are shown in Table 4.
Table thunder source positioning experiment on April 22nd, 4 2019
In table 4, it is (0.143, -24.81,35.06) that microphone array, which is listed at 21 points and measures within 08 minute sound source coordinate, which exists
24.81 degree of directions of south by west, pitch angle are 35.06 degree, and particularly, sound source to array center is only 143 meters, this illustrates the sound source
It should be a certain interference source, rather than measured target thunder source, empty region did not generated thunder in NUIST test point.It utilizes
Thunder source location algorithm provides information feedback for movement routine tracking, helps to improve thunder cloud forecasting and warning accuracy.
The radar return diagram data with 3 same time point of table is chosen, further verifies tracking validity, radar map is such as
(being still the station NUIST electric field instrument position at cross mark in figure) shown in Figure 12-15.Find out from Figure 12-15, from 20: 43
21:16 points are assigned to, down to 5dBZ, only up to 20dBZ, this shows the radar echo intensity of empty region in the test point of the station NUIST
A small amount of charge is had accumulated in cloud layer above the test point, which has no thunder cloud presence.Although test point overhead does not have thunderstorm
Cloud, but still prediction of the tracking to result shown in Figure 11 can be verified by observing radar map 12-15.It has been observed that away from three
The farther away east of atmospheric electric field detector is tieed up, there are large range of thunder cloud region, radar echo intensity is up to 45dBZ;With when
Between passage, thunder cloud region is stood test point far from NUIST gradually, mobile toward South of Jiangsu Province city Yixing City direction.In conjunction with figure
12-15, it is consistent that the movement of thunder cloud during this period with tracking measures result.In fair weather experiment, even if thunder cloud
Region is not covered in empty region in test point, but it is still within following range, shows that thunder cloud movement routine tracks
The validity of method.
According to above-mentioned experimental result, obtain as drawn a conclusion:
This example demonstrates that the three-dimensional electric field component that thunder cloud generates is necessary being.Explore electric field component and thunder cloud
The relationship of orientation physical quantity can make it possible to detect thunder cloud movement routine.By establishing the tracking of thunder cloud movement routine
Spatial model has derived thunder cloud orientation calculation formula, and orientation where it is associated with time point, one kind is proposed and is based on
The thunder cloud movement routine tracking of three-dimensional atmospheric electric field instrument.According to indirect error theory, surveyed to tracking is influenced
Error away from direction finding precision is analyzed.In conjunction with thunder auditory localization and official's radar diagram data, reality has been carried out at the station NUIST
It tests, preferably demonstrates the validity of the tracking.These researchs, can be further while realizing thunder cloud fixing by gross bearings
Real-time tracking is carried out to thunder cloud movement routine, truly realizes thunder cloud forecasting and warning.
Claims (4)
1. a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument, which comprises the following steps:
(1) according to mirror image law theory, the spatial model of tracking thunder cloud movement routine is established, and defines the dynamic in thunder cloud orientation
Parameter;
(2) Potential distribution formula of the thunder cloud at three-dimensional atmospheric electric field instrument is utilized, three-dimensional electric field component is obtained, acquires thunder cloud
Point STThe spherical coordinates of T variation at any time;
(3) the set P of the thunder cloud dynamic coordinate that T changes at any time is defined, and orientation where thunder cloud and time point are carried out
The tracking of thunder cloud movement routine is realized in association.
2. a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument according to claim 1, special
Sign is that the realization process of the step (1) is as follows:
Using N point as coordinate origin, three-dimensional cartesian coordinate system is established, in which: N (0,0,0) is three-dimensional atmospheric electric field instrument position,
ST(xT,yT,zT) it is thunderstorm cloud point STT and the dynamic position that changes, S at any timeT'(xT,yT, 0) and it is STThrowing in x0y plane
Shadow;The sum of the height above sea level of h expression atmospheric electric field detector oneself height and its position;Thunderstorm cloud point STDynamic level drift angle and bow
The elevation angle is respectively αT, βT, and RTFor thunderstorm cloud point STTo the real-time range of electric field instrument N;At electric field instrument N, thunderstorm cloud point S is measuredT
Real-time electric field strength be ET。
3. a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument according to claim 1, special
Sign is, the step (2) the following steps are included:
(21) according to charge inducing in image method and image charge phase scheduling theory, regard thunder cloud region as a charge q, obtain
To thunderstorm cloud point STReal-time Potential distribution at three-dimensional atmospheric electric field instrument NAre as follows:
Wherein, q' is the image charge of point charge q, ε1、ε2The dielectric on ground respectively where air and three-dimensional atmospheric electric field instrument N
Constant;
(22) due to On The Electric Field of Thunderstorm intensity E that atmospheric electric field detector N point measuresT, it is a kind of trivector, it can be to ETIt carries out orthogonal
It decomposes, obtains:
ET=EX,T+EY,T+EZ,T
Wherein, EX,T、EY,T、EZ,TRespectively electric field instrument N measures x-axis, y-axis, the real-time electric field strength of thunder cloud on z-axis direction, and
It is orthogonal two-by-two;
(23) thunderstorm cloud point S is soughtTSpherical coordinates (the R of T variation at any timeT',αT',βT') are as follows:
Wherein, RT' it is thunder cloud STTo the distance R of electric field instrument NTInverting value, αT', βT' it is respectively thunder cloud dynamic level drift angle
αT, pitch angle βTInverting value;As quantity of electric charge q, permittivity ε1, ε2When known, real-time thunder cloud coordinate (RT',αT',βT') only
With three-dimensional electric field vector EX,T、EY,T、EZ,TIt is related.
4. a kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument according to claim 1, special
Sign is, the set P of the dynamic coordinate of T variation at any time of thunder cloud described in step (3) are as follows:
Wherein, T1, T2To TnRespectively refer to for the 1st, 2...n time point, and with time point T1, T2...TnCorresponding thunder cloud coordinate
Respectively
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910692270.6A CN110297284A (en) | 2019-07-30 | 2019-07-30 | A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910692270.6A CN110297284A (en) | 2019-07-30 | 2019-07-30 | A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110297284A true CN110297284A (en) | 2019-10-01 |
Family
ID=68032150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910692270.6A Pending CN110297284A (en) | 2019-07-30 | 2019-07-30 | A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110297284A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487477A (en) * | 2020-05-25 | 2020-08-04 | 南京信息工程大学 | Thunderstorm cloud point charge positioning data complementation method based on atmospheric electric field instrument array group |
CN112213703A (en) * | 2020-08-18 | 2021-01-12 | 成都信息工程大学 | Method and device for inverting cloud parameters by using millimeter wave cloud radar |
CN112766127A (en) * | 2021-01-12 | 2021-05-07 | 南京信息工程大学 | Thundercloud point charge positioning method based on complementary set modal decomposition and SG filtering |
CN113255510A (en) * | 2021-05-21 | 2021-08-13 | 南京信息工程大学 | Rainbow cloud point charge moving path imaging method based on multi-time scale DBSCAN and sample entropy |
CN114812698A (en) * | 2022-07-01 | 2022-07-29 | 江西信息应用职业技术学院 | Lightning protection monitoring method and device, readable storage medium and electronic equipment |
CN116741012A (en) * | 2023-08-14 | 2023-09-12 | 河北省气象行政技术服务中心 | Lightning simulation generating and processing method and device, electronic equipment and storage medium |
CN117709209A (en) * | 2024-02-06 | 2024-03-15 | 水利部南京水利水文自动化研究所 | Thunderstorm cloud charge structure inversion method based on QGA_fmincon |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4245190A (en) * | 1978-12-12 | 1981-01-13 | Lightning Location And Protection, Inc. | Lightning detection system utilizing triangulation and field amplitude comparison techniques |
JPH03222295A (en) * | 1990-01-26 | 1991-10-01 | Central Res Inst Of Electric Power Ind | Laser lightning guiding method |
CN1635721A (en) * | 2003-12-30 | 2005-07-06 | 中国科学院空间科学与应用研究中心 | Regional lightning prewarning method |
CN203054567U (en) * | 2013-01-16 | 2013-07-10 | 曾金全 | Lightning monitoring and early warning system |
CN104614597A (en) * | 2015-01-28 | 2015-05-13 | 南京信息工程大学 | Thunderstorm early-warning method |
CN106353604A (en) * | 2016-08-31 | 2017-01-25 | 国网电力科学研究院武汉南瑞有限责任公司 | Method for early warning thunder and lightning on basis of automatic tracking of atmospheric electric field background noise |
CN109142896A (en) * | 2018-07-25 | 2019-01-04 | 南京信息工程大学 | Lightning Warning method based on three-dimensional atmospheric electric field and MEMD |
CN109917196A (en) * | 2019-01-11 | 2019-06-21 | 南京信息工程大学 | A kind of thunder cloud localization method based on three-dimensional atmospheric electric field instrument visual angle |
-
2019
- 2019-07-30 CN CN201910692270.6A patent/CN110297284A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4245190A (en) * | 1978-12-12 | 1981-01-13 | Lightning Location And Protection, Inc. | Lightning detection system utilizing triangulation and field amplitude comparison techniques |
JPH03222295A (en) * | 1990-01-26 | 1991-10-01 | Central Res Inst Of Electric Power Ind | Laser lightning guiding method |
CN1635721A (en) * | 2003-12-30 | 2005-07-06 | 中国科学院空间科学与应用研究中心 | Regional lightning prewarning method |
CN203054567U (en) * | 2013-01-16 | 2013-07-10 | 曾金全 | Lightning monitoring and early warning system |
CN104614597A (en) * | 2015-01-28 | 2015-05-13 | 南京信息工程大学 | Thunderstorm early-warning method |
CN106353604A (en) * | 2016-08-31 | 2017-01-25 | 国网电力科学研究院武汉南瑞有限责任公司 | Method for early warning thunder and lightning on basis of automatic tracking of atmospheric electric field background noise |
CN109142896A (en) * | 2018-07-25 | 2019-01-04 | 南京信息工程大学 | Lightning Warning method based on three-dimensional atmospheric electric field and MEMD |
CN109917196A (en) * | 2019-01-11 | 2019-06-21 | 南京信息工程大学 | A kind of thunder cloud localization method based on three-dimensional atmospheric electric field instrument visual angle |
Non-Patent Citations (2)
Title |
---|
宋晨曦: "基于地面电场资料的雷暴云电荷结构反演研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
杨春明 等: "大气电场仪在雷电监测预警中的应用", 《现代建筑电气 工程与应用》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487477A (en) * | 2020-05-25 | 2020-08-04 | 南京信息工程大学 | Thunderstorm cloud point charge positioning data complementation method based on atmospheric electric field instrument array group |
CN111487477B (en) * | 2020-05-25 | 2022-04-08 | 南京信息工程大学 | Thunderstorm cloud point charge positioning data complementation method based on atmospheric electric field instrument array group |
CN112213703A (en) * | 2020-08-18 | 2021-01-12 | 成都信息工程大学 | Method and device for inverting cloud parameters by using millimeter wave cloud radar |
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 |
CN113255510A (en) * | 2021-05-21 | 2021-08-13 | 南京信息工程大学 | Rainbow cloud point charge moving path imaging method based on multi-time scale DBSCAN and sample entropy |
CN113255510B (en) * | 2021-05-21 | 2023-07-25 | 南京信息工程大学 | Thunderstorm cloud point charge moving path imaging method based on multi-time scale DBSCAN and sample entropy |
CN114812698A (en) * | 2022-07-01 | 2022-07-29 | 江西信息应用职业技术学院 | Lightning protection monitoring method and device, readable storage medium and electronic equipment |
CN114812698B (en) * | 2022-07-01 | 2022-09-30 | 江西信息应用职业技术学院 | Lightning protection monitoring method and device, readable storage medium and electronic equipment |
CN116741012A (en) * | 2023-08-14 | 2023-09-12 | 河北省气象行政技术服务中心 | Lightning simulation generating and processing method and device, electronic equipment and storage medium |
CN116741012B (en) * | 2023-08-14 | 2023-10-27 | 河北省气象行政技术服务中心 | Lightning simulation generating and processing method and device, electronic equipment and storage medium |
CN117709209A (en) * | 2024-02-06 | 2024-03-15 | 水利部南京水利水文自动化研究所 | Thunderstorm cloud charge structure inversion method based on QGA_fmincon |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110297284A (en) | A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument | |
CN106443830B (en) | A method of Typhoon Monitoring and evaluation monitoring accuracy based on multi-source satellite data | |
Ruohoniemi et al. | Mapping high‐latitude plasma convection with coherent HF radars | |
CN106990401B (en) | full-waveform airborne laser radar data-based two-class elevation error correction method | |
CN108181635A (en) | A kind of laser point cloud sorting technique for transmission line of electricity scissors crossing analysis | |
CN106199605B (en) | Wind field error correcting method | |
CN110646782B (en) | Satellite-borne laser on-orbit pointing calibration method based on waveform matching | |
Teng et al. | Mini-UAV LiDAR for power line inspection | |
CN107843895B (en) | A kind of Dual-Doppler weather radar dimensional wind inversion method | |
CN110487266A (en) | A kind of airborne photoelectric passive high-precision localization method suitable for sea-surface target | |
US20120120230A1 (en) | Apparatus and Method for Small Scale Wind Mapping | |
CN106052694B (en) | Based on the method that gravitational vectors and its gradient tensor carry out locating and tracking to single moving object | |
CN109358325A (en) | The landform inversion method of radar altimeter under rolling topography background | |
CN109917196A (en) | A kind of thunder cloud localization method based on three-dimensional atmospheric electric field instrument visual angle | |
CN107783157A (en) | External sort algorithm intelligence satellite selection method, system and radar platform based on aeronautical satellite | |
CN104251699A (en) | Indoor space positioning equipment and positioning method thereof | |
Coïsson et al. | Tsunami signature in the ionosphere: A simulation of OTH radar observations | |
Yang et al. | A Moving Path Tracking Method of the Thunderstorm Cloud Based on the Three‐Dimensional Atmospheric Electric Field Apparatus | |
Rogers et al. | Measurements and simulation of ionospheric scattering on VHF and UHF radar signals: Channel scattering function | |
CN113848393A (en) | Method for correcting inversion error of ground flashback intensity | |
Wood et al. | Localization of individual lightning discharges via directional and temporal triangulation of sferic measurements at two distant sites | |
Moran et al. | Observations of interplanetary scintillation during the 1998 Whole Sun Month: a comparison between EISCAT, ORT and Nagoya data | |
Wilson | Infrasound from auroral electrojet motions at I53US | |
Jing et al. | Monitoring capabilities of a mobile mapping system based on navigation qualities | |
Mikhailovsky et al. | The features of thunderstorm activity control by different radiophysical measuring instruments (radar, lightning detection systems, ground-based fluxmeters) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191001 |