[go: up one dir, main page]

CN111487477B - Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group - Google Patents

Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group Download PDF

Info

Publication number
CN111487477B
CN111487477B CN202010448671.XA CN202010448671A CN111487477B CN 111487477 B CN111487477 B CN 111487477B CN 202010448671 A CN202010448671 A CN 202010448671A CN 111487477 B CN111487477 B CN 111487477B
Authority
CN
China
Prior art keywords
axis
electric field
measured
cloud point
thunderstorm cloud
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.)
Active
Application number
CN202010448671.XA
Other languages
Chinese (zh)
Other versions
CN111487477A (en
Inventor
行鸿彦
杨旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Information Science and Technology
Original Assignee
Nanjing University of Information Science and Technology
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 Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN202010448671.XA priority Critical patent/CN111487477B/en
Publication of CN111487477A publication Critical patent/CN111487477A/en
Application granted granted Critical
Publication of CN111487477B publication Critical patent/CN111487477B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/08—Measuring electromagnetic field characteristics
    • G01R29/0807—Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814—Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0842—Measurements related to lightning, e.g. measuring electric disturbances, warning systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Electromagnetism (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，具体步骤包括：首先，建立了大气电场仪电场分量测量模型，定义了雷暴云点电荷的方位参数。根据镜像法原理，利用电位分布公式，得到了雷暴云点电荷坐标。定位参数除水平偏角和仰角外，还包括电场仪到雷暴云的距离。以主电场仪为基准，建立了大气电场仪阵列群，提出了一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，对各电场仪测得数据进行互补处理，重新得到雷暴云点电荷位置，有效解决了数据丢失或畸变会对雷暴云点电荷定位的准确性和稳定性产生不利影响的问题。结果表明，该方法能准确反映雷暴云点电荷所在方位，具有较好的定位效果。

Figure 202010448671

The invention discloses a method for complementing thunderstorm cloud point charge location data based on an atmospheric electric field instrument array group. The specific steps include: first, establishing an atmospheric electric field instrument electric field component measurement model, and defining azimuth parameters of thunderstorm cloud point charges. According to the principle of mirror image method, the point charge coordinates of thunderstorm cloud are obtained by using the formula of potential distribution. In addition to the horizontal declination and elevation, the positioning parameters also include the distance from the electric field meter to the thunderstorm cloud. Based on the main electric field meter, an atmospheric electric field instrument array group was established, and a method for complementary location data of thunderstorm cloud point charge based on the atmospheric electric field instrument array group was proposed. The point charge location effectively solves the problem that data loss or distortion will adversely affect the accuracy and stability of thunderstorm cloud point charge location. The results show that the method can accurately reflect the location of the thunderstorm cloud point charge, and has a good positioning effect.

Figure 202010448671

Description

基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group

技术领域technical field

本发明涉及雷暴云探测技术领域，特别是一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法。The invention relates to the technical field of thunderstorm cloud detection, in particular to a method for complementary location data of thunderstorm cloud point charges based on an atmospheric electric field instrument array group.

背景技术Background technique

雷暴云是在一定强度后产生的积雨云，而雷电则是雷暴云中电荷积累而产生的。随着社会现代化的发展，雷电引起的静电感应等效应将给社会经济造成不可估量的损失。雷电放电过程中产生的大电流和强电磁辐射会对地面森林、建筑物、电力电子设备、通信系统等造成破坏，甚至危及人身安全。因此，有效的雷暴云探测方法将有助于减少雷电危害。Thunderstorm clouds are cumulonimbus clouds produced after a certain intensity, while lightning is produced by the accumulation of electric charges in thunderstorm clouds. With the development of social modernization, the electrostatic induction and other effects caused by lightning will cause immeasurable losses to the social economy. The large current and strong electromagnetic radiation generated in the process of lightning discharge will cause damage to ground forests, buildings, power electronic equipment, communication systems, etc., and even endanger personal safety. Therefore, an effective thunderstorm cloud detection method will help reduce lightning hazards.

在气象领域，实时监测大气电场是研究雷暴云分布和雷电预警的重要手段之一。因此，利用三维大气电场仪同时观测大气电场的水平分量和垂直分量，可以得到更全面的雷暴云点电荷方位信息。有研究表明，科研人员已实现了电场分量实时测量，但如果考虑因实际环境的复杂性而可能导致的数据丢失问题，现有方法的实际应用效果还有待进一步验证。现有技术中存在雷暴云点电荷定位性能易受数据丢失影响的问题。In the field of meteorology, real-time monitoring of atmospheric electric fields is one of the important means to study the distribution of thunderstorm clouds and early warning of thunderstorms. Therefore, using the three-dimensional atmospheric electric field instrument to simultaneously observe the horizontal and vertical components of the atmospheric electric field, more comprehensive information on the azimuth of the thunderstorm cloud point charge can be obtained. Studies have shown that researchers have achieved real-time measurement of electric field components, but if the data loss problem that may be caused by the complexity of the actual environment is considered, the practical application effect of the existing method needs to be further verified. In the prior art, there is a problem that the localization performance of thunderstorm cloud point charges is easily affected by data loss.

发明内容SUMMARY OF THE INVENTION

本发明所要解决的技术问题是克服现有技术的不足而提供一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，本发明建立大气电场仪阵列群模型，以进一步提高雷暴云探测能力。The technical problem to be solved by the present invention is to overcome the deficiencies of the prior art and provide a complementary method for thunderstorm cloud point charge location data based on the atmospheric electric field instrument array group. The present invention establishes the atmospheric electric field instrument array group model to further improve the thunderstorm cloud detection. ability.

本发明为解决上述技术问题采用以下技术方案：The present invention adopts the following technical solutions for solving the above-mentioned technical problems:

根据本发明提出的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，包括以下步骤：According to a method for complementing thunderstorm cloud point charge location data based on an atmospheric electric field meter array group proposed in the present invention, the method includes the following steps:

步骤1、建立主大气电场仪模型，其中：主大气电场仪N1所在位置N1(0,0,0)为坐标原点，并以正南方向为X轴正半轴，正东方向为Y轴正半轴，建立三维直角坐标系；Step 1. Establish a model of the main atmospheric electric field meter, in which: the position N1(0,0,0) of the main atmospheric electric field meter N1 is the origin of the coordinates, and the due south direction is the positive semi-axis of the X axis, and the due east direction is the positive Y axis. Semi-axis to establish a three-dimensional Cartesian coordinate system;

由N1测得的雷暴云点电荷M1的坐标为M1(x1,y1,z1)，M1在X-Y平面上的投影为M1'(x1,y1,0)，M1到N1的X轴、Y轴、Z轴的距离分别为x1、y1、z1；主大气电场仪本身高度与其所处海拔之和为h1；N1测得的雷暴云点电荷的水平偏角和仰角分别表示为α1和β1；N1测得的M1到N1的距离为r1，由N1测得的雷暴云点电荷的电场强度为E1；The coordinates of the thunderstorm cloud point charge M1 measured by N1 are M1(x 1 , y 1 , z 1 ), the projection of M1 on the XY plane is M1'(x 1 , y 1 , 0), the X of M1 to N1 The distances of axis, Y axis and Z axis are x 1 , y 1 , z 1 respectively; the sum of the height of the main atmospheric electric field meter itself and its altitude is h1; the horizontal declination and elevation angles of the thunderstorm cloud point charge measured by N1 are respectively Expressed as α1 and β1; the distance from M1 to N1 measured by N1 is r1, and the electric field intensity of the thunderstorm cloud point charge measured by N1 is E1;

步骤2、基于步骤1建立的坐标系，将雷暴云看做一个点电荷q1，得到雷暴云点电荷M1在主大气电场仪N1处的电位分布

Figure BDA0002506728550000026
Step 2. Based on the coordinate system established in Step 1, consider the thunderstorm cloud as a point charge q1, and obtain the potential distribution of the thunderstorm cloud point charge M1 at the main atmospheric electric field meter N1
Figure BDA0002506728550000026

Figure BDA0002506728550000021
Figure BDA0002506728550000021

其中，q1'为点电荷q1的镜像电荷，ε1为空气介电常数，ε2为主大气电场仪N1所在地面的介电常数；Among them, q1' is the mirror charge of the point charge q1, ε 1 is the dielectric constant of air, and ε 2 is the dielectric constant of the surface where the atmospheric electric field meter N1 is located;

对电场强度E1进行正交分解，得到：Orthogonal decomposition of the electric field intensity E1, we get:

E1＝Ex1+Ey1+Ez1；E1=E x1 +E y1 +E z1 ;

其中，Ex1、Ey1、Ez1分别为N1测得X轴、Y轴、Z轴方向上的雷暴云点电荷的电场强度分量，且两两互相垂直；Among them, E x1 , E y1 , and E z1 are the electric field intensity components of the thunderstorm cloud point charges in the X-axis, Y-axis, and Z-axis directions measured by N1, respectively, and the two are perpendicular to each other;

对X，Y，Z轴方向的电位分布

Figure BDA0002506728550000027
进行求导，从而能够得到雷暴云点电荷M1的球坐标(r1,α1,β1)为：Potential distribution for X, Y, Z axis directions
Figure BDA0002506728550000027
By derivation, the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 can be obtained as:

Figure BDA0002506728550000022
Figure BDA0002506728550000022

中间变量

Figure BDA0002506728550000023
中间变量
Figure BDA0002506728550000024
Intermediate variables
Figure BDA0002506728550000023
Intermediate variables
Figure BDA0002506728550000024

结合雷暴云点电荷M1的球坐标(r1,α1,β1)，根据主大气电场仪模型的矢量关系，得到雷暴云点电荷M1的直角坐标(x1,y1,z1)为：Combined with the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1, and according to the vector relationship of the main atmospheric electric field meter model, the rectangular coordinates (x 1 , y 1 , z 1 ) of the thunderstorm cloud point charge M1 are obtained as:

Figure BDA0002506728550000025
Figure BDA0002506728550000025

步骤3、基于步骤1建立的坐标系，建立大气电场仪阵列群模型，且第一副大气电场仪N2和第二副大气电场仪N3的海拔高度与N1相同；N2所在位置为(xN2,yN2,0)，N2到N1的X轴、Y轴、Z轴的距离分别为xN2、yN2、0；N3所在位置为(xN3,yN3,0)，N3到N1的X轴、Y轴、Z轴的距离分别为xN3、yN3、0；Step 3. Based on the coordinate system established in step 1, establish an atmospheric electric field instrument array group model, and the altitudes of the first sub-atmospheric electric field instrument N2 and the second sub-atmospheric electric field instrument N3 are the same as N1; the location of N2 is (x N2 , y N2 ,0), the distances from N2 to N1 's X-axis, Y-axis, Z-axis are x N2 , y N2 , 0; The distances of , Y axis and Z axis are x N3 , y N3 , and 0 respectively;

此时，定义[Ex2,Ey2,Ez2]为N2测得的大气电场值，Ex2、Ey2、Ez2分别为N2测得X轴，Y轴，Z轴方向上的电场强度分量；[Ex3,Ey3,Ez3]为N3测得的大气电场值，Ex3、Ey3、Ez3分别为N3测得X轴，Y轴，Z轴方向上的电场强度分量；At this time, define [E x2 , E y2 , E z2 ] as the atmospheric electric field value measured by N2, and E x2 , E y2 , and E z2 are the electric field intensity components measured by N2 on the X-axis, Y-axis, and Z-axis directions respectively ; [E x3 , E y3 , E z3 ] is the atmospheric electric field value measured by N3, E x3 , E y3 , E z3 are the electric field strength components in the X-axis, Y-axis, and Z-axis directions measured by N3 respectively;

由N1测得雷暴云点电荷M1的坐标为(x1,y1,z1)，采用与步骤2同样的方法，利用[Ex2,Ey2,Ez2]、[Ex3,Ey3,Ez3]，得到由N2、N3直接测得的雷暴云点电荷M1的坐标分别为(X2,Y2,Z2)、(X3,Y3,Z3)，即：由N2测得的M1到N2的X轴、Y轴、Z轴的距离分别为X2、Y2、Z2；由N3测得的M1到N3的X轴、Y轴、Z轴的距离分别为X3、Y3、Z3；那么，从N1的观测角度看，由N2、N3间接测得的雷暴云点电荷M1的坐标分别为(x2,y2,z2)、(x3,y3,z3)，如下：The coordinates of the thunderstorm cloud point charge M1 measured by N1 are (x 1 , y 1 , z 1 ), using the same method as step 2, using [E x2 , E y2 , E z2 ], [E x3 , E y3 , E z3 ], the coordinates of the thunderstorm cloud point charge M1 directly measured by N2 and N3 are (X 2 , Y 2 , Z 2 ), (X 3 , Y 3 , Z 3 ), namely: measured by N2 The distances of the X-axis, Y-axis and Z-axis of M1 to N2 are X 2 , Y 2 , Z 2 respectively; the distances of X-axis, Y-axis and Z-axis of M1 to N3 measured by N3 are X 3 , Y 3 , Z 3 ; then, from the observation point of view of N1, the coordinates of the thunderstorm cloud point charge M1 indirectly measured by N2 and N3 are (x 2 , y 2 , z 2 ), (x 3 , y 3 , z 3 ), as follows:

Figure BDA0002506728550000031
Figure BDA0002506728550000031

式中，由N2测得的M1到N1的X轴、Y轴、Z轴的距离分别为x2、y2、z2；由N3测得的M1到N1的X轴、Y轴、Z轴的距离分别为x3、y3、z3；In the formula, the distances from M1 to N1's X-axis, Y-axis, and Z-axis measured by N2 are x 2 , y 2 , and z 2 respectively; The distances are x 3 , y 3 , z 3 ;

以步骤1建立的坐标系为基准，将N2、N3在X0Y平面上向外展开，联合定位雷暴云点电荷M1；预先设定偏差率阈值P％；对各电场仪测得数据进行互补处理后，重新得到雷暴云点电荷坐标(X,Y,Z)、水平偏角

Figure BDA0002506728550000039
和仰角θ的表达式如下；X、Y、Z分别为重新得到的M1到N1的X轴、Y轴、Z轴的距离；Taking the coordinate system established in step 1 as the benchmark, expand N2 and N3 outward on the X0Y plane to jointly locate the thunderstorm cloud point charge M1; pre-set the deviation rate threshold P%; perform complementary processing on the data measured by each electric field meter. , get the point charge coordinates (X, Y, Z), horizontal declination angle of thunderstorm cloud again
Figure BDA0002506728550000039
and the expression of the elevation angle θ are as follows; X, Y, and Z are the distances of the X-axis, Y-axis, and Z-axis from M1 to N1, respectively;

第一种情况：对于(xi,yi,zi)，i＝1,2,3，当N1与N2、N3测得的各轴数据偏差率均小于P％时，对各轴数据进行互补处理，重新得到的雷暴云点电荷位置为：The first case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured by N1, N2, and N3 is less than P%, the data of each axis are analyzed. Complementary processing, the re-obtained thunderstorm cloud point charge position is:

当

Figure BDA0002506728550000032
时，when
Figure BDA0002506728550000032
hour,

Figure BDA0002506728550000033
Figure BDA0002506728550000033

第二种情况：对于(xi,yi,zi)，i＝1,2,3：The second case: for (x i , y i , z i ), i=1,2,3:

①当N1与N2测得的各轴数据偏差率大于P％，且N1与N3测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：①When the data deviation rate of each axis measured by N1 and N2 is greater than P%, and the data deviation rate of each axis measured by N1 and N3 is less than P%, after processing by the data complementary method, the following expression is obtained:

当

Figure BDA0002506728550000034
或
Figure BDA0002506728550000035
或
Figure BDA0002506728550000036
且
Figure BDA0002506728550000037
Figure BDA0002506728550000038
时，when
Figure BDA0002506728550000034
or
Figure BDA0002506728550000035
or
Figure BDA0002506728550000036
and
Figure BDA0002506728550000037
Figure BDA0002506728550000038
hour,

Figure BDA0002506728550000041
Figure BDA0002506728550000041

②当N1与N3测得的各轴数据偏差率大于P％，且N1与N2测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：②When the data deviation rate of each axis measured by N1 and N3 is greater than P%, and the data deviation rate of each axis measured by N1 and N2 is less than P%, the following expression is obtained after processing by the data complementary method:

当

Figure BDA0002506728550000042
或
Figure BDA0002506728550000043
或
Figure BDA0002506728550000044
且
Figure BDA0002506728550000045
Figure BDA0002506728550000046
时，when
Figure BDA0002506728550000042
or
Figure BDA0002506728550000043
or
Figure BDA0002506728550000044
and
Figure BDA0002506728550000045
Figure BDA0002506728550000046
hour,

Figure BDA0002506728550000047
Figure BDA0002506728550000047

第三种情况：对于(xi,yi,zi)，i＝1,2,3，当N1、N2、N3中任意两个电场仪之间测得的任意轴数据的偏差率均大于P％时，重新将雷暴云点电荷定位数据置零，对应表达式为：The third case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of any axis data measured between any two electric field meters in N1, N2, and N3 is greater than When P%, reset the thunderstorm cloud point charge location data to zero again, and the corresponding expression is:

当

Figure BDA0002506728550000048
或
Figure BDA0002506728550000049
或
Figure BDA00025067285500000410
且
Figure BDA00025067285500000411
或
Figure BDA00025067285500000412
或
Figure BDA00025067285500000413
时，when
Figure BDA0002506728550000048
or
Figure BDA0002506728550000049
or
Figure BDA00025067285500000410
and
Figure BDA00025067285500000411
or
Figure BDA00025067285500000412
or
Figure BDA00025067285500000413
hour,

Figure BDA00025067285500000414
Figure BDA00025067285500000414

第四种情况：对于(xi,yi,zi)，i＝1,2,3，当N2与N3之间测得的各轴数据的偏差率均小于p％，且N2、N3分别与N1测得各轴数据的偏差率均大于P％时，利用数据互补方法处理后，得到如下表达式：The fourth case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured between N2 and N3 is less than p%, and N2 and N3 are respectively When the deviation rate of each axis data measured with N1 is greater than P%, the following expression is obtained after processing by the data complementary method:

当

Figure BDA00025067285500000415
且
Figure BDA00025067285500000416
或
Figure BDA00025067285500000417
或
Figure BDA00025067285500000418
时，when
Figure BDA00025067285500000415
and
Figure BDA00025067285500000416
or
Figure BDA00025067285500000417
or
Figure BDA00025067285500000418
hour,

Figure BDA0002506728550000051
Figure BDA0002506728550000051

作为本发明所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法进一步优化方案，步骤2中，得到雷暴云点电荷M1的球坐标(r1,α1,β1)的具体方法如下：As a further optimization scheme of the method for complementary location data of thunderstorm cloud point charges based on the atmospheric electric field meter array group of the present invention, in step 2, a specific method for obtaining the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 as follows:

电场强度E1，是一种三维矢量，对E1进行正交分解，得到：The electric field intensity E1 is a three-dimensional vector. Orthogonal decomposition of E1 can be obtained:

E1＝Ex1+Ey1+Ez1 (1)E1=E x1 +E y1 +E z1 (1)

其中，Ex1、Ey1、Ez1分别为N1测得X轴、Y轴、Z轴方向上的雷暴云点电荷的电场强度分量，且两两互相垂直。Among them, E x1 , E y1 , and E z1 are the electric field intensity components of the thunderstorm cloud point charges in the X-axis, Y-axis, and Z-axis directions measured by N1, respectively, and the two are perpendicular to each other.

对X，Y，Z轴方向的电位分布

Figure BDA0002506728550000052
进行求导：Potential distribution for X, Y, Z axis directions
Figure BDA0002506728550000052
Do a derivation:

Figure BDA0002506728550000053
Figure BDA0002506728550000053

z1比h1高出2个数量级，那么：z 1 is 2 orders of magnitude higher than h1, then:

z1≈z1-h1≈z1+h1 (3)z 1 ≈z 1 -h1≈z 1 +h1 (3)

基于主大气电场仪模型，M1到N1的距离r1为：Based on the main atmospheric electric field meter model, the distance r1 from M1 to N1 is:

Figure BDA0002506728550000054
Figure BDA0002506728550000054

利用式(3)、式(4)，将式(2)变为：Using formula (3) and formula (4), formula (2) becomes:

Figure BDA0002506728550000055
Figure BDA0002506728550000055

式(5)中，中间变量

Figure BDA0002506728550000056
中间变量
Figure BDA0002506728550000057
In formula (5), the intermediate variable
Figure BDA0002506728550000056
Intermediate variables
Figure BDA0002506728550000057

根据式(5)，得到雷暴云点电荷M1的球坐标(r1,α1,β1)为：According to formula (5), the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 are obtained as:

Figure BDA0002506728550000061
Figure BDA0002506728550000061

作为本发明所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法进一步优化方案，X2、Y2、Z2分别为As a further optimization scheme of the method for complementary location data of thunderstorm cloud point charges based on the atmospheric electric field meter array group, X 2 , Y 2 , and Z 2 are respectively

Figure BDA0002506728550000062
Figure BDA0002506728550000062

其中，r2为N2测得M1到N2的距离，α2为N2测得的雷暴云点电荷的水平偏角，β2为N2测得的雷暴云点电荷的仰角。Among them, r2 is the distance from M1 to N2 measured by N2, α2 is the horizontal declination angle of the thunderstorm cloud point charge measured by N2, and β2 is the elevation angle of the thunderstorm cloud point charge measured by N2.

作为本发明所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法进一步优化方案，X3、Y3、Z3分别为：As a further optimization scheme of the method for complementary location data of thunderstorm cloud point charge based on the atmospheric electric field meter array group described in the present invention, X 3 , Y 3 , and Z 3 are respectively:

Figure BDA0002506728550000063
Figure BDA0002506728550000063

其中，r3为N3测得M1到N3的距离，α3为N3测得的雷暴云点电荷的水平偏角，β3为N3测得的雷暴云点电荷的仰角。Among them, r3 is the distance from M1 to N3 measured by N3, α3 is the horizontal declination angle of the thunderstorm cloud point charge measured by N3, and β3 is the elevation angle of the thunderstorm cloud point charge measured by N3.

本发明采用以上技术方案与现有技术相比，具有以下技术效果：Compared with the prior art, the present invention adopts the above technical scheme, and has the following technical effects:

(1)本发明首先建立单个基于三维大气电场仪的大气电场测量模型，根据镜像法理论，推导了主电场仪的雷暴云点电荷坐标计算公式；在X0Y平面上围绕主电场仪建立电场仪阵列群模型，提出了一种雷暴云点电荷定位数据互补方法；(1) The present invention firstly establishes a single atmospheric electric field measurement model based on a three-dimensional atmospheric electric field instrument, and deduces the calculation formula of the coordinates of the thunderstorm cloud point charge of the main electric field instrument according to the mirror image method; establishes an electric field instrument array around the main electric field instrument on the X0Y plane A swarm model is proposed, and a complementary method of thunderstorm cloud point charge localization data is proposed;

(2)该方法不仅能利用阵列群测得雷暴云点电荷方位数据对其进行互补处理，重新获得点电荷位置，而且能有效解决定位数据丢失问题。(2) This method can not only use the azimuth data of thunderstorm cloud point charge measured by array group to perform complementary processing to obtain the point charge position again, but also can effectively solve the problem of location data loss.

附图说明Description of drawings

图1为主大气电场仪模型。Figure 1 is the main atmospheric electric field instrument model.

图2为大气电场仪阵列群模型。Figure 2 shows the model of the atmospheric electric field instrument array group.

图3为雷暴云点电荷到电场仪距离、电场分量测量误差与测距误差的关系曲线。Figure 3 shows the relationship between the distance from the thunderstorm cloud point charge to the electric field meter, the measurement error of the electric field component and the ranging error.

图4为雷暴云点电荷到电场仪距离、仰角与水平偏角测量误差的关系曲线。Figure 4 is the relationship curve between the distance from the thunderstorm cloud point charge to the electric field meter, the elevation angle and the horizontal declination angle measurement error.

图5为雷暴云点电荷到电场仪距离、仰角与仰角测量误差的关系曲线。Figure 5 shows the relationship between the distance from the thunderstorm cloud point charge to the electric field meter, the elevation angle and the measurement error of the elevation angle.

图6为实际大气电场仪阵列群站点分布图。Fig. 6 is the distribution map of the actual atmospheric electric field instrument array group site.

具体实施方式Detailed ways

为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图及具体实施例对本发明进行详细描述。In order to make the objectives, technical solutions and advantages of the present invention clearer, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

本发明提供一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，包括如下步骤：The present invention provides a method for complementing thunderstorm cloud point charge location data based on an atmospheric electric field meter array group, comprising the following steps:

①建立主大气电场仪模型，其具体为：①Establish the main atmospheric electric field instrument model, which is as follows:

如图1所示，主大气电场仪所在位置N1(0,0,0)为坐标原点，并以正南方向为X轴正半轴，正东方向为Y轴正半轴，建立三维直角坐标系；由N1测得的雷暴云点电荷的坐标为M1(x1,y1,z1)，M1在X-Y平面上的投影为M1'(x1,y1,0)，M1到N1的X轴、Y轴、Z轴的距离分别为x1、y1、z1。主大气电场仪本身高度与其所处海拔之和为h1。雷暴云点电荷的水平偏角和仰角分别表示为α1和β1。M1到N1的距离为r1，由N1测得的雷暴云点电荷的电场强度为E1。As shown in Figure 1, the position N1(0,0,0) of the main atmospheric electric field instrument is the origin of the coordinates, and the positive half-axis of the X-axis is the south direction, and the positive half-axis of the Y-axis is the positive east direction to establish a three-dimensional rectangular coordinate system; the coordinates of the thunderstorm cloud point charge measured by N1 are M1(x 1 , y 1 , z 1 ), the projection of M1 on the XY plane is M1'(x 1 , y 1 , 0), the distance from M1 to N1 The distances of the X-axis, Y-axis, and Z-axis are respectively x 1 , y 1 , and z 1 . The sum of the height of the main atmospheric electric field meter itself and its altitude is h1. The horizontal declination and elevation angles of thunderstorm cloud point charges are denoted as α1 and β1, respectively. The distance from M1 to N1 is r1, and the electric field strength of the thunderstorm cloud point charge measured by N1 is E1.

②将雷暴云看做一个点电荷q1，得到雷暴云点电荷M1在主大气电场仪N1处的电位分布

Figure BDA0002506728550000071
为：② Consider the thunderstorm cloud as a point charge q1, and obtain the potential distribution of the thunderstorm cloud point charge M1 at the main atmospheric electric field meter N1
Figure BDA0002506728550000071
for:

Figure BDA0002506728550000072
Figure BDA0002506728550000072

式(1)中，q1'为点电荷q1的镜像电荷，ε1为空气介电常数，ε2为大气电场仪N所在地面的介电常数。In formula (1), q1' is the image charge of the point charge q1, ε 1 is the dielectric constant of air, and ε 2 is the dielectric constant of the surface where the atmospheric electric field meter N is located.

电场强度E1，是一种三维矢量，对E1进行正交分解，得到：The electric field intensity E1 is a three-dimensional vector. Orthogonal decomposition of E1 can be obtained:

E1＝Ex1+Ey1+Ez1 (2)E1=E x1 +E y1 +E z1 (2)

式(2)中，Ex1、Ey1、Ez1分别为N1测得X轴，Y轴，Z轴方向上的雷暴云点电荷的电场强度分量，且两两互相垂直。In formula (2), E x1 , E y1 , and E z1 are the electric field intensity components of the thunderstorm cloud point charges in the X-axis, Y-axis, and Z-axis directions measured by N1, respectively, and they are perpendicular to each other.

对X，Y，Z轴方向的电位分布

Figure BDA0002506728550000073
进行求导：Potential distribution for X, Y, Z axis directions
Figure BDA0002506728550000073
Do a derivation:

Figure BDA0002506728550000081
Figure BDA0002506728550000081

通常，z1通常比h1高出2个数量级，那么：In general, z1 is usually 2 orders of magnitude higher than h1, then:

z1≈z1-h1≈z1+h1 (4)z 1 ≈z 1 -h1≈z 1 +h1 (4)

基于主大气电场仪模型，M1到N1的距离r1为：Based on the main atmospheric electric field meter model, the distance r1 from M1 to N1 is:

Figure BDA0002506728550000082
Figure BDA0002506728550000082

利用式(4)、式(5)，将式(3)变为：Using formula (4) and formula (5), formula (3) becomes:

Figure BDA0002506728550000083
Figure BDA0002506728550000083

式(6)中，中间变量

Figure BDA0002506728550000084
In formula (6), the intermediate variable
Figure BDA0002506728550000084

根据式(6)，得到雷暴云点电荷M1的球坐标(r1,α1,β1)为：According to formula (6), the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 are obtained as:

Figure BDA0002506728550000085
Figure BDA0002506728550000085

结合雷暴云点电荷M1的球坐标(r1,α1,β1)，根据主大气电场仪模型的矢量关系，得到雷暴云点电荷M1的直角坐标(x1,y1,z1)为：Combined with the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1, and according to the vector relationship of the main atmospheric electric field meter model, the rectangular coordinates (x 1 , y 1 , z 1 ) of the thunderstorm cloud point charge M1 are obtained as:

Figure BDA0002506728550000086
Figure BDA0002506728550000086

③提出一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，其具体为：(3) A complementary method of thunderstorm cloud point charge location data based on the atmospheric electric field meter array group is proposed, which is as follows:

建立图2所示的大气电场仪阵列群模型，基于步骤1建立的坐标系，建立大气电场仪阵列群模型，且第一副大气电场仪N2和第二副大气电场仪N3的海拔高度与N1相同；N2所在位置为(xN2,yN2,0)，N2到N1的X轴、Y轴、Z轴的距离分别为xN2、yN2、0；N3所在位置为(xN3,yN3,0)，N3到N1的X轴、Y轴、Z轴的距离分别为xN3、yN3、0；Establish the atmospheric electric field instrument array group model shown in Figure 2, based on the coordinate system established in step 1, establish the atmospheric electric field instrument array group model, and the altitudes of the first sub-atmospheric electric field instrument N2 and the second sub-atmospheric electric field instrument N3 are the same as N1. The same; the location of N2 is (x N2 , y N2 , 0), the distances from N2 to N1’s X-axis, Y-axis, Z-axis are x N2 , y N2 , 0 respectively; the location of N3 is (x N3 , y N3 ,0), the distances from N3 to N1 X axis, Y axis, Z axis are x N3 , y N3 , 0 respectively;

此时，定义[Ex2,Ey2,Ez2]为N2测得的大气电场值，Ex2、Ey2、Ez2分别为N2测得X轴，Y轴，Z轴方向上的电场强度分量；[Ex3,Ey3,Ez3]为N3测得的大气电场值，Ex3、Ey3、Ez3分别为N3测得X轴，Y轴，Z轴方向上的电场强度分量；At this time, define [E x2 , E y2 , E z2 ] as the atmospheric electric field value measured by N2, and E x2 , E y2 , and E z2 are the electric field intensity components measured by N2 on the X-axis, Y-axis, and Z-axis directions respectively ; [E x3 , E y3 , E z3 ] is the atmospheric electric field value measured by N3, E x3 , E y3 , E z3 are the electric field strength components in the X-axis, Y-axis, and Z-axis directions measured by N3 respectively;

由N1测得雷暴云点电荷M1的坐标为(x1,y1,z1)，采用与步骤2同样的方法，利用[Ex2,Ey2,Ez2]、[Ex3,Ey3,Ez3]，得到由N2、N3直接测得的雷暴云点电荷M1的坐标分别为(X2,Y2,Z2)、(X3,Y3,Z3)，即：由N2测得的M1到N2的X轴、Y轴、Z轴的距离分别为X2、Y2、Z2；由N3测得的M1到N3的X轴、Y轴、Z轴的距离分别为X3、Y3、Z3；那么，从N1的观测角度看，由N2、N3间接测得的雷暴云点电荷M1的坐标分别为(x2,y2,z2)、(x3,y3,z3)，如下：The coordinates of the thunderstorm cloud point charge M1 measured by N1 are (x 1 , y 1 , z 1 ), using the same method as step 2, using [E x2 , E y2 , E z2 ], [E x3 , E y3 , E z3 ], the coordinates of the thunderstorm cloud point charge M1 directly measured by N2 and N3 are (X 2 , Y 2 , Z 2 ), (X 3 , Y 3 , Z 3 ), namely: measured by N2 The distances of the X-axis, Y-axis and Z-axis of M1 to N2 are X 2 , Y 2 , Z 2 respectively; the distances of X-axis, Y-axis and Z-axis of M1 to N3 measured by N3 are X 3 , Y 3 , Z 3 ; then, from the observation point of view of N1, the coordinates of the thunderstorm cloud point charge M1 indirectly measured by N2 and N3 are (x 2 , y 2 , z 2 ), (x 3 , y 3 , z 3 ), as follows:

Figure BDA0002506728550000091
Figure BDA0002506728550000091

式中，由N2测得的M1到N1的X轴、Y轴、Z轴的距离分别为x2、y2、z2；由N3测得的M1到N1的X轴、Y轴、Z轴的距离分别为x3、y3、z3；In the formula, the distances from M1 to N1's X-axis, Y-axis, and Z-axis measured by N2 are x 2 , y 2 , and z 2 respectively; The distances are x 3 , y 3 , z 3 ;

Figure BDA0002506728550000092
Figure BDA0002506728550000092

其中，r2为N2测得M1到N2的距离，α2为N2测得的雷暴云点电荷的水平偏角，β2为N2测得的雷暴云点电荷的仰角。Among them, r2 is the distance from M1 to N2 measured by N2, α2 is the horizontal declination angle of the thunderstorm cloud point charge measured by N2, and β2 is the elevation angle of the thunderstorm cloud point charge measured by N2.

Figure BDA0002506728550000093
Figure BDA0002506728550000093

其中，r3为N3测得M1到N3的距离，α3为N3测得的雷暴云点电荷的水平偏角，β3为N3测得的雷暴云点电荷的仰角。Among them, r3 is the distance from M1 to N3 measured by N3, α3 is the horizontal declination angle of the thunderstorm cloud point charge measured by N3, and β3 is the elevation angle of the thunderstorm cloud point charge measured by N3.

基于图1所示坐标系，将N2、N3在X0Y平面上向外展开，联合定位雷暴云点电荷M1。各电场仪测得雷暴云点电荷定位数据可能会出现数据丢失或失真。针对这种可能性，预先设定偏差率阈值P％，以便更直观表示数据互补方法。Based on the coordinate system shown in Figure 1, N2 and N3 are spread out on the X0Y plane to jointly locate the thunderstorm cloud point charge M1. The location data of thunderstorm cloud point charge measured by each electric field meter may be lost or distorted. Aiming at this possibility, the deviation rate threshold P% is preset in order to express the data complementation method more intuitively.

对各电场仪测得数据进行互补处理后，重新得到雷暴云点电荷坐标(X,Y,Z)、水平偏角

Figure BDA00025067285500001014
和仰角θ的表达式如下：(重新得到的M1到N1的X轴、Y轴、Z轴的距离分别为X、Y、Z)After complementary processing of the data measured by each electric field meter, the point charge coordinates (X, Y, Z) and the horizontal declination angle of the thunderstorm cloud are re-obtained.
Figure BDA00025067285500001014
and the expression of the elevation angle θ are as follows: (The distances of the X-axis, Y-axis, and Z-axis of the re-obtained M1 to N1 are X, Y, and Z, respectively)

1)对于(xi,yi,zi)，i＝1,2,3，当N1与N2、N3测得的各轴数据偏差率均小于P％时，对各轴数据进行互补处理，重新得到的雷暴云点电荷位置为：1) For (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured by N1, N2, and N3 is less than P%, perform complementary processing on each axis data, The re-obtained thunderstorm cloud point charge positions are:

当

Figure BDA0002506728550000101
时，when
Figure BDA0002506728550000101
hour,

Figure BDA0002506728550000102
Figure BDA0002506728550000102

2)对于(xi,yi,zi)，i＝1,2,3：①当N1与N2测得的各轴数据偏差率大于P％，且N1与N3测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：2) For (x i , y i , z i ), i=1, 2, 3: ① When the data deviation rate of each axis measured by N1 and N2 is greater than P%, and the data deviation rate of each axis measured by N1 and N3 When it is less than P%, the following expression is obtained after processing by the data complementary method:

当

Figure BDA0002506728550000103
或
Figure BDA0002506728550000104
或
Figure BDA0002506728550000105
且
Figure BDA0002506728550000106
Figure BDA0002506728550000107
时，when
Figure BDA0002506728550000103
or
Figure BDA0002506728550000104
or
Figure BDA0002506728550000105
and
Figure BDA0002506728550000106
Figure BDA0002506728550000107
hour,

Figure BDA0002506728550000108
Figure BDA0002506728550000108

②当N1与N3测得的各轴数据偏差率大于P％，且N1与N2测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：②When the data deviation rate of each axis measured by N1 and N3 is greater than P%, and the data deviation rate of each axis measured by N1 and N2 is less than P%, the following expression is obtained after processing by the data complementary method:

当

Figure BDA0002506728550000109
或
Figure BDA00025067285500001010
或
Figure BDA00025067285500001011
且
Figure BDA00025067285500001012
Figure BDA00025067285500001013
时，when
Figure BDA0002506728550000109
or
Figure BDA00025067285500001010
or
Figure BDA00025067285500001011
and
Figure BDA00025067285500001012
Figure BDA00025067285500001013
hour,

Figure BDA0002506728550000111
Figure BDA0002506728550000111

3)对于(xi,yi,zi)，i＝1,2,3，当N1、N2、N3中任意两个电场仪之间测得的任意轴数据的偏差率均大于P％时，重新将雷暴云点电荷定位数据置零，对应表达式为：3) For (x i , y i , z i ), i=1, 2, 3, when the deviation rate of any axis data measured between any two electric field meters in N1, N2, N3 is greater than P% , reset the thunderstorm cloud point charge localization data to zero, the corresponding expression is:

当

Figure BDA0002506728550000112
或
Figure BDA0002506728550000113
或
Figure BDA0002506728550000114
且
Figure BDA0002506728550000115
或
Figure BDA0002506728550000116
或
Figure BDA0002506728550000117
时，when
Figure BDA0002506728550000112
or
Figure BDA0002506728550000113
or
Figure BDA0002506728550000114
and
Figure BDA0002506728550000115
or
Figure BDA0002506728550000116
or
Figure BDA0002506728550000117
hour,

Figure BDA0002506728550000118
Figure BDA0002506728550000118

4)对于(xi,yi,zi)，i＝1,2,3，当N2与N3之间测得的各轴数据的偏差率均小于p％，且N2、N3分别与N1测得各轴数据的偏差率均大于P％时，利用数据互补方法处理后，得到如下表达式：4) For (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured between N2 and N3 is less than p%, and N2 and N3 are measured with N1 respectively. When the deviation rate of the data of each axis is greater than P%, the following expression is obtained after processing by the data complementary method:

当

Figure BDA0002506728550000119
且
Figure BDA00025067285500001110
或
Figure BDA00025067285500001111
或
Figure BDA00025067285500001112
时，when
Figure BDA0002506728550000119
and
Figure BDA00025067285500001110
or
Figure BDA00025067285500001111
or
Figure BDA00025067285500001112
hour,

Figure BDA00025067285500001113
Figure BDA00025067285500001113

根据上述实施例方法，以下对基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法从测距测向方面进行性能分析：According to the method of the above-mentioned embodiment, the performance analysis of the complementary method of thunderstorm cloud point charge localization data based on the atmospheric electric field meter array group is performed from the aspect of ranging and direction finding:

结合空气电荷电场分布和雷暴云电荷结构原理，可将电场分量测量标准差设为

Figure BDA00025067285500001114
Combined with the electric field distribution of air charge and the charge structure principle of thunderstorm cloud, the standard deviation of electric field component measurement can be set as
Figure BDA00025067285500001114

1、雷暴云点电荷定位测距测向性能分析1. Performance analysis of thunderstorm cloud point charge location ranging and direction finding

基于间接测量误差理论，由电场分量测量误差

Figure BDA00025067285500001115
引起距离r1，水平偏角α1，仰角β1的测量误差
Figure BDA00025067285500001116
为：Based on the indirect measurement error theory, the measurement error is determined by the electric field component
Figure BDA00025067285500001115
Causes measurement errors of distance r1, horizontal declination angle α1, and elevation angle β1
Figure BDA00025067285500001116
for:

Figure BDA0002506728550000121
Figure BDA0002506728550000121

2、雷暴云点电荷定位测距性能分析2. Analysis of thunderstorm cloud point charge location and ranging performance

利用式(13)，研究距离r1，电场分量测量误差

Figure BDA0002506728550000122
与测距误差
Figure BDA0002506728550000123
的关系，仿真结果如图3所示。图3中，测距误差
Figure BDA0002506728550000124
受距离r1和电场分量测量误差
Figure BDA0002506728550000125
的影响，并且受前者影响较大。测距误差
Figure BDA0002506728550000126
均随着距离r1和电场测量误差
Figure BDA0002506728550000127
的增大而增大。当误差
Figure BDA0002506728550000128
处于0到1kV/m时，测距误差
Figure BDA0002506728550000129
几乎与距离r1的变化无关，误差小于0.06km。此外，当误差
Figure BDA00025067285500001210
大于1kV/m时，测距误差
Figure BDA00025067285500001211
随着距离r1的增大而急剧增大，符合三次函数特性。特别地，当距离r1大于1km时，测距误差
Figure BDA00025067285500001212
随距离r1的增大而线性增大，最大可达0.12km。综上所述，测距误差
Figure BDA00025067285500001213
小于0.12km，充分体现了数据互补方法测距性能的稳定性。Using equation (13), study the distance r1, the measurement error of the electric field component
Figure BDA0002506728550000122
and ranging error
Figure BDA0002506728550000123
The simulation results are shown in Figure 3. In Figure 3, the ranging error
Figure BDA0002506728550000124
Subject to distance r1 and electric field component measurement error
Figure BDA0002506728550000125
, and is greatly affected by the former. ranging error
Figure BDA0002506728550000126
Both with distance r1 and electric field measurement error
Figure BDA0002506728550000127
increases and increases. when the error
Figure BDA0002506728550000128
Ranging error at 0 to 1kV/m
Figure BDA0002506728550000129
Almost independent of the change in distance r1, the error is less than 0.06km. Furthermore, when the error
Figure BDA00025067285500001210
When greater than 1kV/m, ranging error
Figure BDA00025067285500001211
It increases sharply with the increase of the distance r1, which conforms to the characteristics of the cubic function. In particular, when the distance r1 is greater than 1km, the ranging error
Figure BDA00025067285500001212
It increases linearly with the increase of the distance r1, and the maximum can reach 0.12km. In summary, the ranging error
Figure BDA00025067285500001213
It is less than 0.12km, which fully reflects the stability of the ranging performance of the data complementary method.

3、雷暴云点电荷定位测向性能分析3. Analysis of thunderstorm cloud point charge location and direction finding performance

根据式(13)，得到如图4、5所示的测向性能仿真结果。在图3中，当距离r1和仰角β1增大时，水平偏角测量误差

Figure BDA00025067285500001214
会随之增大。特别地，当距离r1处于0到1km时，测量误差
Figure BDA00025067285500001215
几乎不受仰角β1变化的影响，误差
Figure BDA00025067285500001216
小于0.5度。但当距离r1处于1到2km时，误差
Figure BDA00025067285500001217
会随仰角β1的增加呈指数增长，误差
Figure BDA00025067285500001218
将达到1.9度。同样地，在图4中，仰角测量误差
Figure BDA00025067285500001219
随着距离r1和仰角β1的增大而缓慢增大。然而，当距离r1在1至2km范围内时，随着仰角β1的增加，误差
Figure BDA00025067285500001220
以抛物线形式缓慢升至0.16度。According to formula (13), the simulation results of direction finding performance as shown in Figures 4 and 5 are obtained. In Figure 3, when the distance r1 and the elevation angle β1 increase, the horizontal declination measurement error
Figure BDA00025067285500001214
will increase accordingly. In particular, when the distance r1 is between 0 and 1 km, the measurement error
Figure BDA00025067285500001215
Almost unaffected by changes in the elevation angle β1, the error
Figure BDA00025067285500001216
less than 0.5 degrees. But when the distance r1 is in 1 to 2km, the error
Figure BDA00025067285500001217
will increase exponentially with the increase of the elevation angle β1, the error
Figure BDA00025067285500001218
will reach 1.9 degrees. Likewise, in Figure 4, the elevation angle measurement error
Figure BDA00025067285500001219
It increases slowly with the increase of the distance r1 and the elevation angle β1. However, when the distance r1 is in the range of 1 to 2km, as the elevation angle β1 increases, the error
Figure BDA00025067285500001220
It rose slowly to 0.16 degrees in a parabolic fashion.

在实际测量实验中，南京信息工程大学测试站的主大气电场仪安装在电子与信息工程学院，电场仪距离平均海平面约28米，X轴和Y轴的正半轴分别指向南和东。同时，用于组成阵列群的第一副大气电场仪N2、第二副大气电场仪N3分别安装在滨江站和盘城站，如图6所示。In the actual measurement experiment, the main atmospheric electric field instrument of the test station of Nanjing University of Information Science and Technology was installed in the School of Electronics and Information Engineering. The electric field instrument was about 28 meters away from the mean sea level, and the positive semi-axes of the X-axis and Y-axis pointed to the south and east, respectively. At the same time, the first sub-atmospheric electric field instrument N2 and the second sub-atmospheric electric field instrument N3 used to form the array group are installed at Binjiang Station and Pancheng Station, respectively, as shown in Figure 6.

图6中，N2和N3的坐标分别为(0,-1,0)和(-1,1,0)，单位为km。将P％选为20％作为轴间数据的偏差率，并进行了以下两组实验。In Figure 6, the coordinates of N2 and N3 are (0,-1,0) and (-1,1,0), respectively, and the unit is km. P% was chosen as 20% as the deviation rate of the inter-axis data, and the following two sets of experiments were carried out.

1、阴天实验1. Cloudy day experiment

在2019年4月9日10时06分，三维大气电场仪阵列群测量数据(单位:kV/m)如表1所示。At 10:06 on April 9, 2019, the measurement data (unit: kV/m) of the three-dimensional atmospheric electric field meter array group are shown in Table 1.

表1 2019年4月9日10时06分的实验数据Table 1 Experimental data at 10:06 on April 9, 2019

Figure BDA0002506728550000131
Figure BDA0002506728550000131

表1中，从单个电场仪测得的三维电场分量中，不能直接看出它们之间的关系。仅根据水平分量和垂直分量的值，很难判断是否存在数据丢失。因此，进一步获得雷暴云点电荷的位置比较困难。值得注意的是，垂直电场分量较大于水平电场分量，因此，可以初步预测在主电场仪站点上空有雷暴云存在。In Table 1, from the three-dimensional electric field components measured by a single electric field meter, the relationship between them cannot be directly seen. Based on the values of the horizontal and vertical components alone, it is difficult to judge whether there is data loss. Therefore, it is difficult to further obtain the location of the thunderstorm cloud point charge. It is worth noting that the vertical electric field component is larger than the horizontal electric field component, so it can be preliminarily predicted that there will be thunderstorm clouds over the main electric field meter site.

如表1所示，阵列群实测数据所反映的问题符合式(8)及其对应的描述，表示各电场仪在10时06分测得的数据均是正常的。引入数据互补方法后，利用式(8)对实测数据进行处理。基于大气电场仪阵列群，雷暴云点电荷坐标为(0.094，-0.138,0.513)(单位:km)。这表明，点电荷位于南偏西55.74度，且距离主站约539米处。特别地，点电荷仰角较大，达到71.97度，说明主大气电场仪上空区域存在雷暴云，进一步验证了先前的假设。As shown in Table 1, the problem reflected by the measured data of the array group conforms to equation (8) and its corresponding description, indicating that the data measured by each electric field meter at 10:06 are normal. After the data complementation method is introduced, the measured data is processed by formula (8). Based on the atmospheric electric field instrument array group, the coordinates of the thunderstorm cloud point charge are (0.094, -0.138, 0.513) (unit: km). This shows that the point charge is located at 55.74 degrees west-south and about 539 meters from the main station. In particular, the point charge elevation angle is large, reaching 71.97 degrees, indicating that there is a thunderstorm cloud in the area above the main atmospheric electric field meter, which further verifies the previous hypothesis.

2、晴天实验2. Sunny day experiment

在2019年4月22日20时43分，三维大气电场仪阵列群测量数据(单位:kV/m)如表2所示。At 20:43 on April 22, 2019, the measurement data (unit: kV/m) of the three-dimensional atmospheric electric field meter array group are shown in Table 2.

表2 2019年4月22日20时43分的实验数据Table 2 Experimental data at 20:43 on April 22, 2019

Figure BDA0002506728550000132
Figure BDA0002506728550000132

同样地，在表2中，仅根据水平分量和垂直分量的值很难判断是否存在数据丢失。因此，进一步获得雷暴云点电荷的位置将较为困难。总体上，三维大气电场分量值较小。但电场仪测得垂直电场分量有超过1kv/m的较大值，在引入数据互补方法前，可对雷暴云的存在进行猜测。Likewise, in Table 2, it is difficult to judge whether or not there is data loss based only on the values of the horizontal and vertical components. Therefore, it will be more difficult to further obtain the location of the thunderstorm cloud point charge. In general, the three-dimensional atmospheric electric field component value is small. However, the vertical electric field component measured by the electric field meter has a large value of more than 1kv/m. Before the data complementary method is introduced, the existence of thunderstorm clouds can be guessed.

如表2所示，阵列群实测数据所反映的问题符合式(11)及其对应的描述，表示各台电场仪在20时43分测得的数据存在数据丢失。引入数据互补方法后，利用式(11)对实测数据进行处理，得到雷暴云点电荷坐标为(0,0,0)。此时，可判断出雷暴云没有出现在主站上空。对于先前的错误猜测，主要原因是单个电场仪难以判断雷暴云点电荷坐标的真实性。虽然数值相对较大，但在天气晴朗的情况下，这种偏差很可能是由于其他电场干扰，而不是天气本身。As shown in Table 2, the problem reflected by the measured data of the array group conforms to equation (11) and its corresponding description, indicating that the data measured by each electric field meter at 20:43 has data loss. After the data complementation method is introduced, the measured data is processed by formula (11), and the point charge coordinates of the thunderstorm cloud are obtained as (0, 0, 0). At this time, it can be judged that the thunderstorm cloud does not appear over the main station. For the previous wrong guesses, the main reason is that it is difficult for a single electric field meter to judge the authenticity of the coordinates of the point charge of the thunderstorm cloud. While the numbers are relatively large, in clear weather the deviation is likely due to other electric field disturbances rather than the weather itself.

根据上述实验结果，得到如下结论：According to the above experimental results, the following conclusions are drawn:

为减少雷暴云点电荷定位过程中数据丢失对其定位性能的负面影响，提出了一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法。将该方法应用于雷暴云探测中，减少了观测人员仅在单个电场仪正常工作时才能获取雷暴云点电荷方位的局限性。结果表明，与单台电场仪测得数据相比，该方法在点电荷定位方面有更好的效果，有效扩大了雷暴区域的监测范围。In order to reduce the negative impact of data loss on the localization performance of thunderstorm cloud point charges, a data complementary method for thunderstorm cloud point charge localization based on the atmospheric electric field meter array group was proposed. This method is applied to thunderstorm cloud detection, which reduces the limitation that observers can only obtain the point charge orientation of thunderstorm clouds when a single electric field meter is working normally. The results show that, compared with the data measured by a single electric field meter, the method has better effect in the localization of point charges, and effectively expands the monitoring range of the thunderstorm area.

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围内。The above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art who is familiar with the technical scope disclosed by the present invention can easily think of changes or substitutions. All should be covered within the protection scope of the present invention.

Claims (4)

1.一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，其特征在于，包括以下步骤：1. a thunderstorm cloud point charge positioning data complementary method based on atmospheric electric field instrument array group, is characterized in that, comprises the following steps: 步骤1、建立主大气电场仪模型，其中：主大气电场仪N1所在位置N1(0,0,0)为坐标原点，并以正南方向为X轴正半轴，正东方向为Y轴正半轴，建立三维直角坐标系；Step 1. Establish a model of the main atmospheric electric field meter, in which: the position N1(0,0,0) of the main atmospheric electric field meter N1 is the origin of the coordinates, and the due south direction is the positive semi-axis of the X axis, and the due east direction is the positive Y axis. Semi-axis to establish a three-dimensional Cartesian coordinate system; 由N1测得的雷暴云点电荷M1的坐标为M1(x1,y1,z1)，M1在X-Y平面上的投影为M1'(x1,y1,0)，M1到N1的X轴、Y轴、Z轴的距离分别为x1、y1、z1；主大气电场仪本身高度与其所处海拔之和为h1；N1测得的雷暴云点电荷的水平偏角和仰角分别表示为α1和β1；N1测得的M1到N1的距离为r1，由N1测得的雷暴云点电荷的电场强度为E1；The coordinates of the thunderstorm cloud point charge M1 measured by N1 are M1(x 1 , y 1 , z 1 ), the projection of M1 on the XY plane is M1'(x 1 , y 1 , 0), the X of M1 to N1 The distances of axis, Y axis and Z axis are x 1 , y 1 , z 1 respectively; the sum of the height of the main atmospheric electric field meter itself and its altitude is h1; the horizontal declination and elevation angles of the thunderstorm cloud point charge measured by N1 are respectively Expressed as α1 and β1; the distance from M1 to N1 measured by N1 is r1, and the electric field intensity of the thunderstorm cloud point charge measured by N1 is E1; 步骤2、基于步骤1建立的坐标系，将雷暴云看做一个点电荷q1，得到雷暴云点电荷M1在主大气电场仪N1处的电位分布
Figure FDA0002506728540000011
Step 2. Based on the coordinate system established in Step 1, consider the thunderstorm cloud as a point charge q1, and obtain the potential distribution of the thunderstorm cloud point charge M1 at the main atmospheric electric field meter N1
Figure FDA0002506728540000011
Figure FDA0002506728540000012
Figure FDA0002506728540000012
 其中，q1'为点电荷q1的镜像电荷，ε1为空气介电常数，ε2为主大气电场仪N1所在地面的介电常数；Among them, q1' is the mirror charge of the point charge q1, ε 1 is the dielectric constant of air, and ε 2 is the dielectric constant of the surface where the atmospheric electric field meter N1 is located; 对电场强度E1进行正交分解，得到：Orthogonal decomposition of the electric field intensity E1, we get: E1＝Ex1+Ey1+Ez1；E1=E x1 +E y1 +E z1 ; 其中，Ex1、Ey1、Ez1分别为N1测得X轴、Y轴、Z轴方向上的雷暴云点电荷的电场强度分量，且两两互相垂直；Among them, E x1 , E y1 , and E z1 are the electric field intensity components of the thunderstorm cloud point charges in the X-axis, Y-axis, and Z-axis directions measured by N1, respectively, and the two are perpendicular to each other; 对X，Y，Z轴方向的电位分布
Figure FDA0002506728540000013
进行求导，从而能够得到雷暴云点电荷M1的球坐标(r1,α1,β1)为：Potential distribution for X, Y, Z axis directions
Figure FDA0002506728540000013
By derivation, the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 can be obtained as:
Figure FDA0002506728540000014
Figure FDA0002506728540000014
 中间变量
Figure FDA0002506728540000015
中间变量
Figure FDA0002506728540000016
Intermediate variables
Figure FDA0002506728540000015
Intermediate variables
Figure FDA0002506728540000016
 结合雷暴云点电荷M1的球坐标(r1,α1,β1)，根据主大气电场仪模型的矢量关系，得到雷暴云点电荷M1的直角坐标(x1,y1,z1)为：Combined with the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1, and according to the vector relationship of the main atmospheric electric field meter model, the rectangular coordinates (x 1 , y 1 , z 1 ) of the thunderstorm cloud point charge M1 are obtained as:
Figure FDA0002506728540000021
Figure FDA0002506728540000021
 步骤3、基于步骤1建立的坐标系，建立大气电场仪阵列群模型，且第一副大气电场仪N2和第二副大气电场仪N3的海拔高度与N1相同；N2所在位置为(xN2,yN2,0)，N2到N1的X轴、Y轴、Z轴的距离分别为xN2、yN2、0；N3所在位置为(xN3,yN3,0)，N3到N1的X轴、Y轴、Z轴的距离分别为xN3、yN3、0；Step 3. Based on the coordinate system established in step 1, establish an atmospheric electric field instrument array group model, and the altitudes of the first sub-atmospheric electric field instrument N2 and the second sub-atmospheric electric field instrument N3 are the same as N1; the location of N2 is (x N2 , y N2 ,0), the distances from N2 to N1 's X-axis, Y-axis, Z-axis are x N2 , y N2 , 0; The distances of , Y axis and Z axis are x N3 , y N3 , and 0 respectively; 此时，定义[Ex2,Ey2,Ez2]为N2测得的大气电场值，Ex2、Ey2、Ez2分别为N2测得X轴，Y轴，Z轴方向上的电场强度分量；[Ex3,Ey3,Ez3]为N3测得的大气电场值，Ex3、Ey3、Ez3分别为N3测得X轴，Y轴，Z轴方向上的电场强度分量；At this time, define [E x2 , E y2 , E z2 ] as the atmospheric electric field value measured by N2, and E x2 , E y2 , and E z2 are the electric field intensity components measured by N2 on the X-axis, Y-axis, and Z-axis directions respectively ; [E x3 , E y3 , E z3 ] is the atmospheric electric field value measured by N3, E x3 , E y3 , E z3 are the electric field strength components in the X-axis, Y-axis, and Z-axis directions measured by N3 respectively; 由N1测得雷暴云点电荷M1的坐标为(x1,y1,z1)，采用与步骤2同样的方法，利用[Ex2,Ey2,Ez2]、[Ex3,Ey3,Ez3]，得到由N2、N3直接测得的雷暴云点电荷M1的坐标分别为(X2,Y2,Z2)、(X3,Y3,Z3)，即：由N2测得的M1到N2的X轴、Y轴、Z轴的距离分别为X2、Y2、Z2；由N3测得的M1到N3的X轴、Y轴、Z轴的距离分别为X3、Y3、Z3；那么，从N1的观测角度看，由N2、N3间接测得的雷暴云点电荷M1的坐标分别为(x2,y2,z2)、(x3,y3,z3)，如下：The coordinates of the thunderstorm cloud point charge M1 measured by N1 are (x 1 , y 1 , z 1 ), using the same method as step 2, using [E x2 , E y2 , E z2 ], [E x3 , E y3 , E z3 ], the coordinates of the thunderstorm cloud point charge M1 directly measured by N2 and N3 are (X 2 , Y 2 , Z 2 ), (X 3 , Y 3 , Z 3 ), namely: measured by N2 The distances of the X-axis, Y-axis and Z-axis of M1 to N2 are X 2 , Y 2 , Z 2 respectively; the distances of X-axis, Y-axis and Z-axis of M1 to N3 measured by N3 are X 3 , Y 3 , Z 3 ; then, from the observation point of view of N1, the coordinates of the thunderstorm cloud point charge M1 indirectly measured by N2 and N3 are (x 2 , y 2 , z 2 ), (x 3 , y 3 , z 3 ), as follows:
Figure FDA0002506728540000022
Figure FDA0002506728540000022
 式中，由N2测得的M1到N1的X轴、Y轴、Z轴的距离分别为x2、y2、z2；由N3测得的M1到N1的X轴、Y轴、Z轴的距离分别为x3、y3、z3；In the formula, the distances from M1 to N1's X-axis, Y-axis, and Z-axis measured by N2 are x 2 , y 2 , and z 2 respectively; The distances are x 3 , y 3 , z 3 ; 以步骤1建立的坐标系为基准，将N2、N3在X0Y平面上向外展开，联合定位雷暴云点电荷M1；预先设定偏差率阈值P％；对各电场仪测得数据进行互补处理后，重新得到雷暴云点电荷坐标(X,Y,Z)、水平偏角
Figure FDA0002506728540000023
和仰角θ的表达式如下；X、Y、Z分别为重新得到的M1到N1的X轴、Y轴、Z轴的距离；Taking the coordinate system established in step 1 as the benchmark, expand N2 and N3 outward on the X0Y plane to jointly locate the thunderstorm cloud point charge M1; pre-set the deviation rate threshold P%; perform complementary processing on the data measured by each electric field meter. , get the point charge coordinates (X, Y, Z), horizontal declination angle of thunderstorm cloud again
Figure FDA0002506728540000023
and the expression of the elevation angle θ are as follows; X, Y, and Z are the distances of the X-axis, Y-axis, and Z-axis from M1 to N1, respectively; 第一种情况：对于(xi,yi,zi)，i＝1,2,3，当N1与N2、N3测得的各轴数据偏差率均小于P％时，对各轴数据进行互补处理，重新得到的雷暴云点电荷位置为：The first case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured by N1, N2, and N3 is less than P%, the data of each axis are analyzed. Complementary processing, the re-obtained thunderstorm cloud point charge position is: 当
Figure FDA0002506728540000031
时，when
Figure FDA0002506728540000031
hour,
Figure FDA0002506728540000032
Figure FDA0002506728540000032
 第二种情况：对于(xi,yi,zi)，i＝1,2,3：The second case: for (x i , y i , z i ), i=1,2,3: ①当N1与N2测得的各轴数据偏差率大于P％，且N1与N3测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：①When the data deviation rate of each axis measured by N1 and N2 is greater than P%, and the data deviation rate of each axis measured by N1 and N3 is less than P%, after processing by the data complementary method, the following expression is obtained: 当
Figure FDA0002506728540000033
或
Figure FDA0002506728540000034
或
Figure FDA0002506728540000035
且
Figure FDA0002506728540000036
Figure FDA0002506728540000037
时，when
Figure FDA0002506728540000033
or
Figure FDA0002506728540000034
or
Figure FDA0002506728540000035
and
Figure FDA0002506728540000036
Figure FDA0002506728540000037
hour,
Figure FDA0002506728540000038
Figure FDA0002506728540000038
 ②当N1与N3测得的各轴数据偏差率大于P％，且N1与N2测得各轴数据偏差率小于P％时，利用数据互补方法处理后，得到如下表达式：②When the data deviation rate of each axis measured by N1 and N3 is greater than P%, and the data deviation rate of each axis measured by N1 and N2 is less than P%, the following expression is obtained after processing by the data complementary method: 当
Figure FDA0002506728540000039
或
Figure FDA00025067285400000310
或
Figure FDA00025067285400000311
且
Figure FDA00025067285400000312
Figure FDA00025067285400000313
时，when
Figure FDA0002506728540000039
or
Figure FDA00025067285400000310
or
Figure FDA00025067285400000311
and
Figure FDA00025067285400000312
Figure FDA00025067285400000313
hour,
Figure FDA00025067285400000314
Figure FDA00025067285400000314
 第三种情况：对于(xi,yi,zi)，i＝1,2,3，当N1、N2、N3中任意两个电场仪之间测得的任意轴数据的偏差率均大于P％时，重新将雷暴云点电荷定位数据置零，对应表达式为：The third case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of any axis data measured between any two electric field meters in N1, N2, and N3 is greater than When P%, reset the thunderstorm cloud point charge location data to zero again, and the corresponding expression is: 当
Figure FDA00025067285400000315
或
Figure FDA00025067285400000316
或
Figure FDA00025067285400000317
且
Figure FDA00025067285400000318
或
Figure FDA0002506728540000041
或
Figure FDA0002506728540000042
时，when
Figure FDA00025067285400000315
or
Figure FDA00025067285400000316
or
Figure FDA00025067285400000317
and
Figure FDA00025067285400000318
or
Figure FDA0002506728540000041
or
Figure FDA0002506728540000042
hour,
Figure FDA0002506728540000043
Figure FDA0002506728540000043
 第四种情况：对于(xi,yi,zi)，i＝1,2,3，当N2与N3之间测得的各轴数据的偏差率均小于p％，且N2、N3分别与N1测得各轴数据的偏差率均大于P％时，利用数据互补方法处理后，得到如下表达式：The fourth case: for (x i , y i , z i ), i=1, 2, 3, when the deviation rate of each axis data measured between N2 and N3 is less than p%, and N2 and N3 are respectively When the deviation rate of each axis data measured with N1 is greater than P%, the following expression is obtained after processing by the data complementary method: 当
Figure FDA0002506728540000044
且
Figure FDA0002506728540000045
或
Figure FDA0002506728540000046
或
Figure FDA0002506728540000047
时，when
Figure FDA0002506728540000044
and
Figure FDA0002506728540000045
or
Figure FDA0002506728540000046
or
Figure FDA0002506728540000047
hour,
Figure FDA0002506728540000048
Figure FDA0002506728540000048
 2.根据权利要求1所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，其特征在于，步骤2中，得到雷暴云点电荷M1的球坐标(r1,α1,β1)的具体方法如下：2. The method for locating data of thunderstorm cloud point charges based on atmospheric electric field meter array group according to claim 1, wherein in step 2, the spherical coordinates (r1, α1, β1 of the thunderstorm cloud point charge M1 are obtained) ) are as follows: 电场强度E1，是一种三维矢量，对E1进行正交分解，得到：The electric field intensity E1 is a three-dimensional vector. Orthogonal decomposition of E1 can be obtained: E1＝Ex1+Ey1+Ez1 (1)E1=E x1 +E y1 +E z1 (1) 其中，Ex1、Ey1、Ez1分别为N1测得X轴、Y轴、Z轴方向上的雷暴云点电荷的电场强度分量，且两两互相垂直；Among them, E x1 , E y1 , and E z1 are the electric field intensity components of the thunderstorm cloud point charges in the X-axis, Y-axis, and Z-axis directions measured by N1, respectively, and the two are perpendicular to each other; 对X，Y，Z轴方向的电位分布
Figure FDA0002506728540000049
进行求导：Potential distribution for X, Y, Z axis directions
Figure FDA0002506728540000049
Do a derivation:
Figure FDA0002506728540000051
Figure FDA0002506728540000051
 z1比h1高出2个数量级，那么：z 1 is 2 orders of magnitude higher than h1, then: z1≈z1-h1≈z1+h1 (3)z 1 ≈z 1 -h1≈z 1 +h1 (3) 基于主大气电场仪模型，M1到N1的距离r1为：Based on the main atmospheric electric field meter model, the distance r1 from M1 to N1 is:
Figure FDA0002506728540000052
Figure FDA0002506728540000052
 利用式(3)、式(4)，将式(2)变为：Using formula (3) and formula (4), formula (2) becomes:
Figure FDA0002506728540000053
Figure FDA0002506728540000053
 式(5)中，中间变量
Figure FDA0002506728540000054
中间变量
Figure FDA0002506728540000055
In formula (5), the intermediate variable
Figure FDA0002506728540000054
Intermediate variables
Figure FDA0002506728540000055
 根据式(5)，得到雷暴云点电荷M1的球坐标(r1,α1,β1)为：According to formula (5), the spherical coordinates (r1, α1, β1) of the thunderstorm cloud point charge M1 are obtained as:
Figure FDA0002506728540000056
Figure FDA0002506728540000056
 3.根据权利要求1所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，其特征在于，X2、Y2、Z2分别为3 . The method for complementing thunderstorm cloud point charge location data based on the atmospheric electric field meter array group according to claim 1 , wherein X 2 , Y 2 , and Z 2 are respectively
Figure FDA0002506728540000057
Figure FDA0002506728540000057
 其中，r2为N2测得M1到N2的距离，α2为N2测得的雷暴云点电荷的水平偏角，β2为N2测得的雷暴云点电荷的仰角。Among them, r2 is the distance from M1 to N2 measured by N2, α2 is the horizontal declination angle of the thunderstorm cloud point charge measured by N2, and β2 is the elevation angle of the thunderstorm cloud point charge measured by N2. 4.根据权利要求1所述的一种基于大气电场仪阵列群的雷暴云点电荷定位数据互补方法，其特征在于，X3、Y3、Z3分别为：4. A method for complementing thunderstorm cloud point charge location data based on an atmospheric electric field meter array group according to claim 1, wherein X 3 , Y 3 , and Z 3 are respectively:
Figure FDA0002506728540000061
Figure FDA0002506728540000061
 其中，r3为N3测得M1到N3的距离，α3为N3测得的雷暴云点电荷的水平偏角，β3为N3测得的雷暴云点电荷的仰角。Among them, r3 is the distance from M1 to N3 measured by N3, α3 is the horizontal declination angle of the thunderstorm cloud point charge measured by N3, and β3 is the elevation angle of the thunderstorm cloud point charge measured by N3.
CN202010448671.XA 2020-05-25 2020-05-25 Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group Active CN111487477B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010448671.XA CN111487477B (en) 2020-05-25 2020-05-25 Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010448671.XA CN111487477B (en) 2020-05-25 2020-05-25 Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group

Publications (2)

Publication Number Publication Date
CN111487477A CN111487477A (en) 2020-08-04
CN111487477B true CN111487477B (en) 2022-04-08

Family

ID=71795582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010448671.XA Active CN111487477B (en) 2020-05-25 2020-05-25 Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group

Country Status (1)

Country Link
CN (1) CN111487477B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112766127B (en) * 2021-01-12 2023-08-22 南京信息工程大学 Thundercloud Point Charge Location Method Based on Complementary Ensemble Mode Decomposition and SG Filter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107238826A (en) * 2017-06-09 2017-10-10 杨波 The method being distributed using Dual-Polarized Doppler Weather Radar echo inverting thunderstorm charge
CN109917196A (en) * 2019-01-11 2019-06-21 南京信息工程大学 A method for locating thunderstorm clouds based on the perspective of a three-dimensional atmospheric electric field meter
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
CN110297284A (en) * 2019-07-30 2019-10-01 南京信息工程大学 A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107238826A (en) * 2017-06-09 2017-10-10 杨波 The method being distributed using Dual-Polarized Doppler Weather Radar echo inverting thunderstorm charge
CN109917196A (en) * 2019-01-11 2019-06-21 南京信息工程大学 A method for locating thunderstorm clouds based on the perspective of a three-dimensional atmospheric electric field meter
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
CN110297284A (en) * 2019-07-30 2019-10-01 南京信息工程大学 A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Sound source omnidirectional positioning calibration method based on microphone observation angle;Xu Yang 等;《Complexity》;20181004;第1-15页 *
基于地面电场的雷暴云电荷结构混合反演;宋晨曦 等;《电瓷避雷器》;20151231(第6期);第73-79页 *

Also Published As

Publication number Publication date
CN111487477A (en) 2020-08-04

Similar Documents

Publication Publication Date Title
CN105279884B (en) A kind of meteorological transmission line forest fire sprawling method for early warning of combination public network
Gao et al. Three-dimensional propagation characteristics of the upward connecting leaders in six negative tall-object flashes in Guangzhou
CN109212475B (en) Dual-aircraft passive positioning method based on azimuth and pitch information
CN110174557B (en) Thunderstorm cloud positioning and calibrating method based on three-dimensional atmospheric electric field instrument observation visual angle
CN112305559A (en) Power transmission line distance measuring method, device and system based on ground fixed-point laser radar scanning and electronic equipment
CN110297284A (en) A kind of thunder cloud movement routine tracking based on three-dimensional atmospheric electric field instrument
CN107085167B (en) A kind of transmission line failure localization method based on big data
CN109917196A (en) A method for locating thunderstorm clouds based on the perspective of a three-dimensional atmospheric electric field meter
CN111487477B (en) Complementary method of thunderstorm cloud point charge location data based on atmospheric electric field array group
CN109581383A (en) Anisotropy sea radar raster-displaying analogy method and system
TWI745153B (en) Method, device and system of ultra-wideband localization
CN108073180A (en) Control method, the device and system of unmanned plane
CN112766127B (en) Thundercloud Point Charge Location Method Based on Complementary Ensemble Mode Decomposition and SG Filter
CN107883917A (en) High voltage transmission line tower height extracting method based on multidate shade difference
CN105988114A (en) Filling Method of Doppler Radar Wind Velocity Data
CN108362926A (en) A kind of voltage class recognition methods and device
CN207249147U (en) A kind of atmospheric electric field detector networking Lightning Warning system based on triangulation
JP2017156282A (en) Charge amount estimation method, estimation device, and estimation program
CN206388308U (en) A kind of automation foundation supervising device based on the Big Dipper
Sun et al. Mosaic research with 3D lidar point cloud of civil aircraft
CN115185300A (en) Unmanned aerial vehicle patrol route safety detection method, device, equipment and storage medium
CN112858979B (en) Thunderstorm cloud point charge positioning altitude correction method based on three-dimensional atmospheric electric field measurement
CN116086599B (en) Method, device, equipment and storage medium for acquiring low-altitude radiation sound field of aircraft
Yang et al. Analysis of shielding failure in transmission lines considering complex terrain
Liu et al. Research on Grid Information Model (GIM) Technology of Transmission Line Based on Full-area Feature Weight Template Matching Recognition

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
GR01 Patent grant
GR01 Patent grant