CN113514834B - Wind speed and direction inversion method, device, equipment and storage medium - Google Patents

Wind speed and direction inversion method, device, equipment and storage medium Download PDF

Info

Publication number
CN113514834B
CN113514834B CN202111071562.1A CN202111071562A CN113514834B CN 113514834 B CN113514834 B CN 113514834B CN 202111071562 A CN202111071562 A CN 202111071562A CN 113514834 B CN113514834 B CN 113514834B
Authority
CN
China
Prior art keywords
wind speed
wind
radial
real
wind direction
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
CN202111071562.1A
Other languages
Chinese (zh)
Other versions
CN113514834A (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.)
Aerospace Hongtu Information Technology Co Ltd
Original Assignee
Aerospace Hongtu Information Technology Co Ltd
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 Aerospace Hongtu Information Technology Co Ltd filed Critical Aerospace Hongtu Information Technology Co Ltd
Priority to CN202111071562.1A priority Critical patent/CN113514834B/en
Publication of CN113514834A publication Critical patent/CN113514834A/en
Application granted granted Critical
Publication of CN113514834B publication Critical patent/CN113514834B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00—Indicating or recording presence, absence, or direction, of movement
    • G01P13/02—Indicating direction only, e.g. by weather vane
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The application provides a wind speed and direction inversion method, a device, equipment and a storage medium, which relate to the technical field of radar wind measurement and specifically comprise the following steps: acquiring radial wind measurement data of at least two available meteorological radars around a measurement point, and combining the available meteorological radars two by two to form a plurality of combinations; based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring point, the real wind speed and the real wind direction of the multiple measuring points corresponding to each combination are obtained through inversion; and carrying out statistical processing on the real wind speed and the real wind direction of the plurality of measurement points to obtain the final real wind speed and the real wind direction of the measurement points. According to the method and the device, based on the radial wind data of the two meteorological radars at the measuring point, the real wind speed and the wind direction of the measuring point are obtained through inversion, and the accuracy of radar wind measurement is improved.

Description

Wind speed and direction inversion method, device, equipment and storage medium
Technical Field
The invention relates to the technical field of radar wind measurement, in particular to a wind speed and direction inversion method, device, equipment and storage medium.
Background
At present, the measurement of wind direction and wind speed distribution in the air by a meteorological radar is mainly based on the Doppler effect. Due to the limitation of the measurement mechanism of the method, the measured wind speed is the wind speed relative to the direction of the connecting line of the measurement point and the radar, and the measured wind direction is the direction deviating from or facing the radar on the connecting line, namely the wind direction and the wind speed measured by the radar are the wind speed projection of the real wind direction and the wind speed of the measurement point on the connecting line of the measurement point and the center of the radar, and the wind speed projection is only the component of the real wind direction and the wind speed, namely radial wind, and cannot represent the real wind direction and the wind speed of the measurement point. For the radar measurement data, misunderstanding of the air stroke distribution by a forecaster is often caused, and the applicability of the measurement data is reduced.
At present, many meteorological radars are networked in China, data are shared, but measurement of distribution of wind direction and wind speed in the air is still in measurement of radial wind, and the measurement is not true wind direction and wind speed.
Disclosure of Invention
In view of this, the present application provides a wind speed and direction inversion method, apparatus, device and medium, so as to improve the accuracy of wind speed and direction measurement of the existing meteorological radar based on the doppler effect.
On one hand, the embodiment of the application provides a wind speed and direction inversion method, which comprises the following steps:
acquiring radial wind measurement data of at least two available meteorological radars around a measurement point, and combining the available meteorological radars two by two to form a plurality of combinations;
based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring point, the real wind speed and the real wind direction of the multiple measuring points corresponding to each combination are obtained through inversion;
and carrying out statistical processing on the real wind speed and the real wind direction of the plurality of measurement points to obtain the final real wind speed and the real wind direction of the measurement points.
Further, the distance between two weather radars in a combination is not more than half of the maximum measurement distance of the weather radar.
Further, when the origin-point meteorological radar in any combination is the first radar and the non-origin-point meteorological radar is the second radar, the inversion is performed based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measurement point to obtain the real wind speed and the real wind direction of the multiple measurement points corresponding to each combination, and the inversion method includes:
obtaining v of measuring point at moment to be measured from measuring data of first weather radar1And a radial wind direction α;
obtaining the radial wind speed v of the measuring point at the same moment from the measuring data of the second meteorological radar2And a radial wind direction β;
according to radial wind speed v1And radial wind direction alpha, and radial wind speed v2And radial wind direction beta, calculating the real wind direction d and the real wind speed s of the measuring point:
Figure P_210910151732805_805495001
wherein the parameters
Figure P_210910151732852_852360001
Arctan () is an arctangent function.
Further, when the origin-point meteorological radar in any combination is the first radar and the non-origin-point meteorological radar is the second radar, the inversion is performed based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measurement point to obtain the real wind speed and the real wind direction of the multiple measurement points corresponding to each combination, and the inversion method includes:
determining the distribution range of the possible wind speed of the measuring point according to the real wind speed of the right-angle characteristic measuring point between the first meteorological radar and the second meteorological radar;
determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial wind speed of the measuring points;
sampling a possible wind speed distribution range according to a preset speed interval to obtain M sampling wind speeds, and sampling a possible wind direction distribution range according to a preset angle interval to obtain N sampling wind directions; for any combination between the M sampling wind speeds and the N sampling wind directions, calculating two corresponding radial wind speeds, and generating an inversion radial wind speed set;
calculating the error between each element in the inversion radial wind speed set and the radial wind speed of the measuring point to obtain an error set;
and taking the corresponding sampling wind direction and sampling wind speed which are smaller than the minimum value of the error threshold value in the error set as the real wind speed and the real wind direction of the measuring point.
Further, determining the distribution range of the possible wind speed of the measuring point according to the real wind speed of the right-angle characteristic measuring point between the first meteorological radar and the second meteorological radar; the method comprises the following steps:
determining a right-angle characteristic measurement point between the first meteorological radar and the second meteorological radar; the right-angle characteristic measuring point is positioned on a circle with the diameter of the connecting line of the first meteorological radar and the second meteorological radar;
radial wind speed v based on first meteorological radar and second meteorological radar at right-angle characteristic measurement point01And v02Calculating the true wind speed v of the right-angle feature measurement point0：
Figure P_210910151732885_885579001
Measuring the true velocity v of a point from a right-angle feature0Determining the distribution range of the possible wind speeds of the measuring points: when v is0When the wind speed is less than or equal to 10m/s, the distribution range of the possible wind speeds of the measuring points is (0, 2 v)0) Otherwise, the possible wind speed distribution range of the measuring points is as follows: (v)0/2，2 v0）。
Further, determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial wind speed of the measuring points; the method comprises the following steps:
taking the first radar as an origin and the due north as the Y-axis positive direction, when the radial wind speed of the first radar at the measuring point is greater than zero, the first possible wind direction distribution range of the measuring point is (alpha-90)o，α+90o) When the radial wind speed of the first radar at the measuring point is less than zero, the first possible wind direction distribution range of the measuring point is (alpha + 90)o，α+270o) (ii) a Wherein alpha is an included angle between a connecting line of the measuring point and the first meteorological radar and due north; if the value in the first possible wind direction distribution range is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting 360 to this valueoThen convert it into wind direction value interval [0 ]o, 360o）；
Taking the second meteorological radar as an origin, taking the true north as the Y-axis forward direction, and when the radial wind speed of the second meteorological radar at the measuring point is greater than zero, the possible wind direction distribution range of the measuring point is (beta-90)o，β+90o) When the radial wind speed of the second meteorological radar at the measuring point is less than zero, the second possible wind direction distribution range of the measuring point is (beta + 90)o，β+270o) (ii) a Wherein beta is a clamp between the connecting line of the measuring point and the second meteorological radar and the true northAn angle; if the value in the second possible wind direction distribution range is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting to this value
Figure M_210910151732932_932954001
Then convert it into wind direction value interval [0 ]o, 360o）；
And performing cross operation on the first possible wind direction distribution range and the second possible wind direction distribution range, wherein the obtained overlapped part is the possible wind direction distribution range of the measuring point.
Further, for any combination between the M sampling wind speeds and the N sampling wind directions, two corresponding radial wind speeds are calculated, and the two radial wind speeds are used as an element to generate an inversion radial wind speed set; the method comprises the following steps:
let the mth sampling wind speed be smThe nth sampling wind direction is dnThen the radial wind speed of the first weather radar at the measuring point
Figure P_210910151732979_979834001
And the radial wind speed of the second meteorological radar at the measurement point
Figure P_210910151733011_011066002
Comprises the following steps:
Figure P_210910151733042_042309001
where l = M × N, denotes the number in M × N combinations of wind speed and wind direction;
will be provided with
Figure P_210910151733092_092108001
And
Figure P_210910151733123_123386002
as an element of the set of inverted radial wind speeds, the total number of elements of the set is M N.
On the other hand, an embodiment of the present application provides a wind speed and direction inversion apparatus, including:
the acquisition unit is used for acquiring the measurement data of at least two available meteorological radars around the measurement point and combining the available meteorological radars two by two to form a plurality of combinations;
the inversion unit is used for obtaining the real wind speed and the real wind direction of a plurality of measuring points corresponding to each combination through inversion based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring points;
and the statistical unit is used for performing statistical processing on the real wind speed and the real wind direction of the plurality of measurement points to obtain the final real wind speed and the real wind direction of the measurement points.
On the other hand, an embodiment of the present application provides an electronic device, including: the wind speed and direction inversion method comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the computer program to realize the wind speed and direction inversion method of the embodiment of the application.
On the other hand, the embodiment of the present application provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the method for inverting wind speed and wind direction according to the embodiment of the present application is implemented.
Compared with the prior art, the beneficial effects of the embodiment of the application are as follows:
the embodiment of the application utilizes two meteorological radars to measure radial wind in the same region simultaneously, and utilizes the measured data to invert the true wind direction and wind speed of each measuring point in the region, thereby improving the accuracy of radar wind direction and wind speed measurement and making up the authenticity and integrity of Doppler meteorological radar wind direction and wind speed measurement.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flowchart of a wind speed and direction inversion method provided in an embodiment of the present application;
FIG. 2 is a schematic diagram of radial velocities of two orthogonal feature measurement points of a radar provided in an embodiment of the present application;
fig. 3 is a functional structure schematic diagram of a wind speed and direction inversion apparatus provided in an embodiment of the present application;
fig. 4 is a schematic diagram of a hardware structure of an electronic device according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
First, the design idea of the embodiment of the present application is briefly introduced.
Currently, when a single weather radar measures wind speed and wind direction based on the doppler effect, the measured wind speed is the wind speed in the direction of the connection line between the measurement point and the radar, and the measured wind direction is the direction departing from or facing the radar on the connection line, that is, the wind direction and the wind speed measured by the radar are radial wind, and cannot represent the true wind direction and the wind speed of the measurement point.
In order to solve the problem, the embodiment of the application is based on the radial wind measurement data of the two meteorological radars in the same region, the real wind direction and wind speed of each measurement point in the region are inverted, the applicable efficiency of the wind direction and wind speed of the meteorological radars is improved, and the authenticity and integrity of the measurement of the wind direction and wind speed of the Doppler meteorological radars are compensated.
For the radial wind measurement data of the same region based on the two meteorological radars, the two methods adopted by the embodiment of the application invert the real wind direction and speed of each measurement point in the region.
The first method adopts a direct geometric equation to solve, and directly solves the true wind speed and the true wind direction of a measuring point according to the mathematical relationship between the true wind speed and the true wind direction of the measuring point and the wind directions and the wind speeds of radial wind of two radars at the measuring point.
The second method adopts a searching method, firstly, a right-angle characteristic measuring point between two meteorological radars is found, and the distribution range of the possible wind speed of the measuring point is determined according to the real wind speed of the right-angle characteristic measuring point; then, determining the distribution range of the possible wind directions of the measuring points by utilizing the positive and negative radial speeds of the measuring points; sampling a possible wind speed distribution range according to a preset speed interval to obtain M sampling wind speeds, sampling the possible wind direction distribution range according to a preset angle interval to obtain N sampling wind directions, calculating two corresponding radial wind speeds for any combination of the M sampling wind speeds and the N sampling wind directions, and generating an inversion radial wind speed set by taking the two radial wind speeds as an element; calculating the error between each element in the inversion radial wind speed set and the radial wind speed of the measuring point to obtain an error set; and finally, taking the corresponding sampling wind direction and sampling wind speed of the minimum value in the error set as the real wind speed and real wind direction of the measuring point.
After introducing the application scenario and the design concept of the embodiment of the present application, the following describes a technical solution provided by the embodiment of the present application.
The first embodiment is as follows:
as shown in fig. 1, an embodiment of the present application provides a wind speed and direction inversion method, including:
step 101: acquiring measurement data of at least two available meteorological radars around a measurement point, and combining the available meteorological radars two by two to form a plurality of combinations;
the distance between the two meteorological radars A and B is preferably less than half of the maximum measurement distance of the meteorological radars, the meteorological radar A is used as an origin survey station, and the meteorological radar B is used as a non-origin survey station.
Step 102: based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring point, the real wind speed and the real wind direction of the multiple measuring points corresponding to each combination are obtained through inversion;
in step 102, for any combination of meteorological radar A and meteorological radar B, the inversion of the true wind speed and true wind direction of the measurement point can be achieved in two ways.
The first mode is a geometric formula method, which specifically comprises the following steps:
step 1021: acquiring the radial wind speed v of the meteorological radar A at a measuring point at a moment to be measured1The wind direction angle alpha is an included angle between a connecting line of the measuring point and the meteorological radar A and the true north;
step 1022: obtaining the radial wind speed v of the meteorological radar B at the measuring point at the same moment2And wind direction angle β; the wind direction angle beta is an included angle between a connecting line of the measuring point and the meteorological radar B and the true north;
in step 1022, since the measurement of the weather radar a and the weather radar B are not completely synchronized, the measurement data of the weather radar B can be regarded as two measurement data satisfying the simultaneity as long as the measurement time of the two measurement data is less than 1 minute.
Step 1023: calculating the real wind direction d and the real wind speed s of the measuring point at the time t;
in step 1023, according to the radial wind speed v1Wind direction angle alpha, and radial wind speed v2And a wind direction angle beta, calculating a real wind direction angle d and a real wind speed s according to the following formulas:
Figure P_210910151733154_154842001
wherein the parameters
Figure P_210910151733185_185865001
Arctan () is an arctangent function.
The second method is a search method, which specifically includes:
step 102-1: determining a right-angle characteristic measuring point between the meteorological radars A and B, and calculating the real speed of the right-angle characteristic measuring point, thereby determining the distribution range of the possible wind speed of the measuring point;
in step 102-1, for a group of measurement points in a certain area, a measurement point can always be found, and the measurement point and the right triangle at the position of the two measurement radars are formed, and the measurement point is located at the right corner point, as shown in fig. 2. The specific method comprises the following steps: and connecting the two radars, drawing a circle by taking the midpoint of the connecting line as a circular point and taking one half of the distance of the connecting line as a radius, wherein any point on the circle can be used as a right-angle characteristic measuring point. Based on the position characteristics of the right-angle feature points, the real speed v of the right-angle feature measurement points0The evolution of the sum of the squares of the radial velocities of the two radars at this measurement point is:
Figure P_210910151733248_248386001
wherein v is01And v02The radial wind speeds of the meteorological radars A and B at the right-angle characteristic measuring point are respectively.
According to the characteristics of the atmospheric continuity distribution, the possible wind speed distribution range of the measuring point is generally as follows: (v)0/2，2 v0) When v is0At 10m/s or less, the possible wind speed distribution of the measuring points is as follows: (0, 2 v)0）。
Step 102-2: determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial speed of the measuring points;
in step 102-2, with the meteorological radar a as the origin and the due north as the Y-axis forward direction, when the radial wind speed of the measurement point relative to the meteorological radar a is greater than zero, that is, the radial wind of the measurement point is toward the meteorological radar a, the possible wind direction distribution range of the measurement point is (α -90)o，α+90o) When the measuring point is radial to the meteorological radar AThe wind speed is less than zero, namely the radial wind of the measuring point deviates from the meteorological radar A, and the possible wind direction distribution range of the measuring point is (alpha + 90)o，α+270o) (ii) a Wherein alpha is an included angle between a connecting line of the measuring point and the meteorological radar A and true north;
if the value in the possible wind direction distribution range of the measuring point is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting 360 to this valueoThen convert it into wind direction value interval [0 ]o, 360o) (ii) a The method specifically comprises the following steps:
when alpha-90o<Wind direction range (alpha-90) at 0o，α+90o) To (360)o+α-90o，0o) And (0 °, α + 90)o) (ii) a When alpha +90o>360oTime, wind direction range (alpha-90)o，α+90o) Is converted into (alpha-90 DEG, 360 DEG)
Figure M_210910151733282_282023001
) And (0, alpha + 90)o-360o). When alpha +90o>360oIn the wind direction range (alpha + 90)o，α+270o) To (alpha + 90)o－360°，α+270o-360 °); when alpha +90o<360 DEG and alpha +270
Figure M_210910151733313_313788002
>360
Figure M_210910151733345_345063003
The range of the time wind direction is changed to (alpha + 90)o360 DEG and (0 DEG, alpha +270 DEG)o－360°）。
Taking the meteorological radar B as an origin, when the radial wind speed of a measuring point relative to the meteorological radar B is greater than zero, namely the radial wind of the measuring point faces the meteorological radar B, the possible wind direction distribution range of the measuring point is (beta-90)o，β+90o) When the radial wind speed of the measuring point relative to the meteorological radar A is less than zero, namely the radial wind of the measuring point deviates from the meteorological radar B, the true wind direction distribution range of the measuring pointIs (beta + 90)o，β+270o) (ii) a Wherein beta is an included angle between a connecting line of the measuring point and the meteorological radar B and due north.
If the value in the possible wind direction distribution range of the measuring point is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting 360 degrees to this value, it is converted into a wind direction value interval [0 ]o, 360o) (ii) a The method specifically comprises the following steps:
when beta-90 degree<Wind direction range (beta-90) at 0 deg.Co，β+90o) Is turned to (360 degrees + beta-90 degrees, 0
Figure M_210910151733376_376297001
) And (0, beta + 90)o) (ii) a When beta +90o>360 DEG time, wind direction range (beta-90)o，β+90o) Is converted into (beta-90 DEG, 360 DEG)
Figure M_210910151733423_423187002
) And (0 degrees, beta +90 degrees-360 degrees). When beta +90 DEG>Wind direction range (beta + 90) at 360 DEGo，β+270o) To (beta + 90)o－360°，β+270o-360 °); when beta +90 DEG<360 DEG and beta +270 DEG>The wind direction range is changed into (beta +90 degrees, 360 degrees) and (0 degrees, beta +270 degrees-360 degrees) at 360 degrees;
and performing cross operation on two possible wind direction distribution ranges determined by the two meteorological radars, wherein the overlapped part is the possible wind direction distribution range of the measuring point.
Step 102-3: sampling a real wind speed distribution range according to a preset speed interval to obtain M possible wind speeds, sampling a real wind direction distribution range according to a preset angle interval to obtain N possible wind directions, combining any one sampling wind speed and any one sampling wind direction, and calculating the radial speed of wind formed by the combination relative to the two radars;
in step 102-3, the speed interval is preferably 0.1 m/s and the angle interval is preferably 1 degree, which may be determined based on the calculation accuracy requirements.
If the mth sampling wind speed is smM is more than or equal to 1 and less than or equal to M, and the nth sampling wind direction angle is dnN is more than or equal to 1 and less than or equal to N, smAnd dnThe combination of the data to obtain a wind measurement data, the radial wind speed of the meteorological radar A at the measuring point
Figure P_210910151733454_454424001
And the radial wind speed of the meteorological radar B at the measuring point
Figure P_210910151733489_489601002
Comprises the following steps:
Figure P_210910151733536_536958001
where l = M × N, denotes the number in M × N combinations of wind speed and wind direction;
will be provided with
Figure P_210910151733568_568180001
And
Figure P_210910151733599_599435002
as an element of the set of inverted radial wind speeds, the total number of elements of the set is M N.
Step 102-4: calculating the error between each element in the inversion radial wind speed set and the measured radial wind speed to obtain an error set;
item I for the set of inverted radial wind speeds
Figure P_210910151733646_646370001
And
Figure P_210910151733680_680463002
error Err from two measured radial wind speedslComprises the following steps:
Figure P_210910151733712_712238001
or, error ErrlComprises the following steps:
Figure P_210910151733743_743508001
err islAs an element of the error set, the total number of elements of the set is M × N.
Step 102-5: and taking the corresponding sampling wind direction and sampling wind speed which are smaller than the minimum value of the error threshold value in the error set as the real wind direction and the real wind speed of the measuring point.
In step 102-5, the error threshold is preferably 0.2 m/s. The purpose of setting the error threshold is to ensure that the true wind speed and true wind direction of the inverted measurement point are reliable. If all the errors in the error set are larger than the error threshold value, the wind speed and the wind direction obtained by inversion do not meet the precision requirement, and the possible reason is that the distribution range is not properly selected, and at this time, the possible wind speed distribution range and the possible wind direction distribution range need to be determined again, and the searching process needs to be carried out again.
Step 103: and carrying out statistical processing on the real wind speed and the real wind direction of the plurality of measurement points to obtain the final real wind speed and the real wind direction of the measurement points.
In step 103, if K weather radars are distributed around the measurement point, theoretically, two or two of the K weather radars are combined to obtain
Figure P_210910151733790_790356001
A radar combination, then obtaining
Figure P_210910151733821_821598002
And calculating the final real wind speed and real wind direction of the measuring point by solving the mean value of the multiple values or finding the value of the maximum distribution probability in the multiple values.
Example two:
based on the foregoing embodiments, an embodiment of the present application provides an apparatus for inverting a wind speed and a wind direction, and referring to fig. 3, an apparatus 200 for inverting a wind speed and a wind direction provided by an embodiment of the present application at least includes:
an obtaining unit 201, configured to obtain measurement data of at least two available weather radars around a measurement point, and combine the available weather radars two by two to form multiple combinations;
the inversion unit 202 is configured to invert to obtain real wind speeds and real wind directions of a plurality of measurement points corresponding to each combination based on the radial wind speeds and the radial wind directions of the two available meteorological radars at the measurement points in each combination;
in one possible implementation, the inversion unit 202 is specifically configured to:
acquiring v of meteorological radar A at a measuring point at a moment to be measured1And wind direction angle α; the wind direction angle alpha is an included angle between a connecting line of the measuring point and the meteorological radar A and the true north;
obtaining the radial wind speed v of the meteorological radar B at the measuring point at the same moment2And wind direction angle β; the wind direction angle beta is an included angle between a connecting line of the measuring point and the meteorological radar B and the true north;
according to radial wind speed v1And radial wind direction alpha, and radial wind speed v2And radial wind direction beta, calculating the real wind direction d and the real wind speed s of the measuring point:
Figure P_210910151733868_868132001
wherein the parameters
Figure P_210910151733914_914927001
Arctan () is an arctangent function.
In another possible implementation, the inversion unit 202 is specifically configured to:
determining a right-angle characteristic measuring point between the meteorological radars A and B, and calculating the real speed of the right-angle characteristic measuring point, thereby determining the distribution range of the possible wind speed of the measuring point;
determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial speed of the measuring points;
sampling a real wind speed distribution range according to a preset speed interval to obtain M possible wind speeds, sampling a real wind direction distribution range according to a preset angle interval to obtain N possible wind directions, combining any one sampling wind speed and any one sampling wind direction, and calculating the radial speed of wind formed by the combination relative to the two radars;
calculating the error between each element in the inversion radial wind speed set and the measured radial wind speed to obtain an error set;
and taking the corresponding sampling wind direction and sampling wind speed of the minimum value in the error set as the real wind direction and the real wind speed of the measuring point.
The statistical unit 203 is configured to perform statistical processing on the real wind speed and the real wind direction of the multiple measurement points to obtain a final real wind speed and a final real wind direction of the measurement point.
It should be noted that, because the principle of the wind speed and direction inversion apparatus 200 provided in the embodiment of the present application for solving the technical problem is similar to that of the wind speed and direction inversion method provided in the embodiment of the present application, the implementation of the wind speed and direction inversion apparatus 200 provided in the embodiment of the present application can refer to the implementation of the wind speed and direction inversion method provided in the embodiment of the present application, and repeated details are not repeated.
Example three:
based on the foregoing embodiments, an embodiment of the present application further provides an electronic device, and referring to fig. 4, an electronic device 300 provided in an embodiment of the present application at least includes: the wind speed and direction inversion method comprises a processor 301, a memory 302 and a computer program stored on the memory 302 and capable of running on the processor 301, wherein the wind speed and direction inversion method provided by the embodiment of the application is realized when the processor 301 executes the computer program.
The electronic device 300 provided by the embodiment of the present application may further include a bus 303 connecting different components (including the processor 301 and the memory 302). Bus 303 represents one or more of any of several types of bus structures, including a memory bus, a peripheral bus, a local bus, and so forth.
The Memory 302 may include readable media in the form of volatile Memory, such as Random Access Memory (RAM) 3021 and/or cache Memory 3022, and may further include Read Only Memory (ROM) 3023.
The memory 302 may also include a program tool 3024 having a set (at least one) of program modules 3025, the program modules 3025 including, but not limited to: an operating subsystem, one or more application programs, other program modules, and program data, each of which, or some combination thereof, may comprise an implementation of a network environment.
Electronic device 300 may also communicate with one or more external devices 304 (e.g., keyboard, remote control, etc.), with one or more devices that enable a user to interact with electronic device 300 (e.g., cell phone, computer, etc.), and/or with any device that enables electronic device 300 to communicate with one or more other electronic devices 300 (e.g., router, modem, etc.). This communication may be through an Input/Output (I/O) interface 303. Also, the electronic device 300 may communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), and/or a public Network, such as the internet) via the Network adapter 306. As shown in FIG. 4, the network adapter 306 communicates with the other modules of the electronic device 300 via the bus 303. It should be understood that although not shown in FIG. 4, other hardware and/or software modules may be used in conjunction with electronic device 300, including but not limited to: microcode, device drivers, Redundant processors, external disk drive Arrays, disk array (RAID) subsystems, tape drives, and data backup storage subsystems, to name a few.
It should be noted that the electronic device 300 shown in fig. 4 is only an example, and should not bring any limitation to the functions and the scope of the application of the embodiments.
Example four:
the embodiment of the application also provides a computer-readable storage medium, which stores computer instructions, and the computer instructions are executed by a processor to implement the wind speed and direction inversion method provided by the embodiment of the application.
It should be noted that although several units or sub-units of the apparatus are mentioned in the above detailed description, such division is merely exemplary and not mandatory. Indeed, the features and functions of two or more units described above may be embodied in one unit, according to embodiments of the application. Conversely, the features and functions of one unit described above may be further divided into embodiments by a plurality of units.
Further, while the operations of the methods of the present application are depicted in the drawings in a particular order, this does not require or imply that these operations must be performed in this particular order, or that all of the illustrated operations must be performed, to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step execution, and/or one step broken down into multiple step executions.
While the preferred embodiments of the present application have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all alterations and modifications as fall within the scope of the application.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (11)

1. A wind speed and direction inversion method is characterized by comprising the following steps:
acquiring radial wind measurement data of at least two available meteorological radars around a measurement point, and combining the available meteorological radars two by two to form a plurality of combinations;
based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring point, the real wind speed and the real wind direction of the multiple measuring points corresponding to each combination are obtained through inversion;
carrying out statistical processing on the real wind speed and the real wind direction of the plurality of measuring points to obtain the final real wind speed and the real wind direction of the measuring points;
when the origin-point weather radar in any combination is the first weather radar and the non-origin-point weather radar is the second weather radar, the true wind speed and the true wind direction of the plurality of measurement points corresponding to each combination are obtained by inversion based on the radial wind speed and the radial wind direction of the two available weather radars in each combination at the measurement points, and the method comprises the following steps:
obtaining the radial wind speed v of the measuring point at the moment to be measured from the measuring data of the first weather radar1And a radial wind direction α;
obtaining the radial wind speed v of the measuring point at the same moment from the measuring data of the second meteorological radar2And a radial wind direction β;
according to radial wind speed v1And radial wind direction alpha, and radial wind speed v2And radial wind direction beta, calculating the real wind direction d and the real wind speed s of the measuring point:
Figure P_211202103214952_952687001
wherein the parameters
Figure P_211202103215015_015020001
Arctan () is an arctangent function.
2. The wind speed and direction inversion method of claim 1, wherein the distance between two meteorological radars in a combination is no more than half of the maximum measurement distance of the meteorological radars.
3. An apparatus for inverting a wind speed and a wind direction, comprising:
the acquisition unit is used for acquiring the measurement data of at least two available meteorological radars around the measurement point and combining the available meteorological radars two by two to form a plurality of combinations;
the inversion unit is used for obtaining the real wind speed and the real wind direction of a plurality of measuring points corresponding to each combination through inversion based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring points;
the statistical unit is used for performing statistical processing on the real wind speed and the real wind direction of the plurality of measuring points to obtain the final real wind speed and the real wind direction of the measuring points;
when the origin observation station weather radar in any combination is the first weather radar and the non-origin observation station weather radar is the second weather radar, the inversion unit is specifically configured to:
obtaining the radial wind speed v of the measuring point at the moment to be measured from the measuring data of the first weather radar1And a radial wind direction α;
obtaining the radial wind speed v of the measuring point at the same moment from the measuring data of the second meteorological radar2And a radial wind direction β;
according to radial wind speed v1And radial wind direction alpha, and radial wind speed v2And radial wind direction beta, calculating the real wind direction d and the real wind speed s of the measuring point:
Figure P_211202103215046_046509001
wherein the parameters
Figure P_211202103215093_093248001
Arctan () is an arctangent function.
4. A wind speed and direction inversion method is characterized by comprising the following steps:
acquiring radial wind measurement data of at least two available meteorological radars around a measurement point, and combining the available meteorological radars two by two to form a plurality of combinations;
based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring point, the real wind speed and the real wind direction of the multiple measuring points corresponding to each combination are obtained through inversion;
carrying out statistical processing on the real wind speed and the real wind direction of the plurality of measuring points to obtain the final real wind speed and the real wind direction of the measuring points;
when the origin-point weather radar in any combination is the first weather radar and the non-origin-point weather radar is the second weather radar, the true wind speed and the true wind direction of the plurality of measurement points corresponding to each combination are obtained by inversion based on the radial wind speed and the radial wind direction of the two available weather radars in each combination at the measurement points, and the method comprises the following steps:
determining the distribution range of the possible wind speed of the measuring point according to the real wind speed of the right-angle characteristic measuring point between the first meteorological radar and the second meteorological radar;
determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial wind speed of the measuring points;
sampling a possible wind speed distribution range according to a preset speed interval to obtain M sampling wind speeds, and sampling a possible wind direction distribution range according to a preset angle interval to obtain N sampling wind directions; for any combination between the M sampling wind speeds and the N sampling wind directions, calculating two corresponding radial wind speeds, and generating an inversion radial wind speed set;
calculating the error between each element in the inversion radial wind speed set and the radial wind speed of the measuring point to obtain an error set;
and taking the corresponding sampling wind direction and sampling wind speed which are smaller than the minimum value of the error threshold value in the error set as the real wind speed and the real wind direction of the measuring point.
5. The wind speed and direction inversion method of claim 4, wherein the distance between two meteorological radars in a combination is not more than half of the maximum measurement distance of the meteorological radars.
6. The wind speed and direction inversion method of claim 4, wherein the distribution range of the possible wind speed of the measuring point is determined according to the real wind speed of the right-angle characteristic measuring point between the first meteorological radar and the second meteorological radar; the method comprises the following steps:
determining a right-angle characteristic measurement point between the first meteorological radar and the second meteorological radar; the right-angle characteristic measuring point is positioned on a circle with the diameter of the connecting line of the first meteorological radar and the second meteorological radar;
radial wind speed v based on first meteorological radar and second meteorological radar at right-angle characteristic measurement point01And v02Calculating the true wind speed v of the right-angle feature measurement point0：
Figure P_211202103215125_125960001
Measuring the true velocity v of a point from a right-angle feature0Determining the distribution range of the possible wind speeds of the measuring points: when v is0When the wind speed is less than or equal to 10m/s, the distribution range of the possible wind speeds of the measuring points is (0, 2 v)0) Otherwise, the possible wind speed distribution range of the measuring points is as follows: (v)0/2，2 v0）。
7. The wind speed and direction inversion method of claim 4, wherein the distribution range of possible wind directions of the measuring points is determined according to the positive and negative of the radial wind speed of the measuring points; the method comprises the following steps:
taking the first meteorological radar as an origin and the due north as the Y-axis forward direction, and when the radial wind speed of the first meteorological radar at a measuring point is greater than zero, the first possible wind direction distribution range of the measuring point is (alpha-90)o，α+90o) When the radial wind speed of the first weather radar at the measuring point is less than zero, the first possible wind direction distribution range of the measuring point is (alpha + 90)o，α+270o) (ii) a Wherein alpha is an included angle between a connecting line of the measuring point and the first meteorological radar and due north; if the value in the first possible wind direction distribution range is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting 360 to this valueoThen convert it into wind direction value interval [0 ]o, 360o）；
Taking the second meteorological radar as an origin, taking the true north as the Y-axis forward direction, and when the radial wind speed of the second meteorological radar at the measuring point is greater than zero, the possible wind direction distribution range of the measuring point is (beta-90)o，β+90o) When the radial wind speed of the second meteorological radar at the measuring point is less than zero, the second possible wind direction distribution range of the measuring point is (beta + 90)o，β+270o) (ii) a Wherein beta is an included angle between a connecting line of the measuring point and the second meteorological radar and due north; if the value in the second possible wind direction distribution range is not in the wind direction angle value range [0 ]o, 360o) By adding or subtracting 360 to this valueoThen convert it into wind direction value interval [0 ]o, 360o）；
And performing cross operation on the first possible wind direction distribution range and the second possible wind direction distribution range, wherein the obtained overlapped part is the possible wind direction distribution range of the measuring point.
8. The wind speed and direction inversion method according to claim 7, wherein for any combination between M sampled wind speeds and N sampled wind directions, two corresponding radial wind speeds are calculated, and an inversion radial wind speed set is generated with the two radial wind speeds as an element; the method comprises the following steps:
let the mth sampling wind speed be smThe nth sampling wind direction is dnThen the radial wind speed of the first weather radar at the measuring point
Figure P_211202103215157_157747001
And the radial wind speed of the second meteorological radar at the measurement point
Figure P_211202103215189_189031002
Comprises the following steps:
Figure P_211202103215220_220263001
where l = M × N, denotes the number in M × N combinations of wind speed and wind direction;
will be provided with
Figure P_211202103215282_282832001
And
Figure P_211202103215316_316499002
as an element of the set of inverted radial wind speeds, the total number of elements of the set is M N.
9. An apparatus for inverting a wind speed and a wind direction, comprising:
the acquisition unit is used for acquiring the measurement data of at least two available meteorological radars around the measurement point and combining the available meteorological radars two by two to form a plurality of combinations;
the inversion unit is used for obtaining the real wind speed and the real wind direction of a plurality of measuring points corresponding to each combination through inversion based on the radial wind speed and the radial wind direction of the two available meteorological radars in each combination at the measuring points;
the statistical unit is used for performing statistical processing on the real wind speed and the real wind direction of the plurality of measuring points to obtain the final real wind speed and the real wind direction of the measuring points;
when the origin observation station weather radar in any combination is the first weather radar and the non-origin observation station weather radar is the second weather radar, the inversion unit is specifically configured to:
determining the distribution range of the possible wind speed of the measuring point according to the real wind speed of the right-angle characteristic measuring point between the first meteorological radar and the second meteorological radar;
determining the distribution range of possible wind directions of the measuring points according to the positive and negative of the radial wind speed of the measuring points;
sampling a possible wind speed distribution range according to a preset speed interval to obtain M sampling wind speeds, and sampling a possible wind direction distribution range according to a preset angle interval to obtain N sampling wind directions; for any combination between the M sampling wind speeds and the N sampling wind directions, calculating two corresponding radial wind speeds, and generating an inversion radial wind speed set;
calculating the error between each element in the inversion radial wind speed set and the radial wind speed of the measuring point to obtain an error set;
and taking the corresponding sampling wind direction and sampling wind speed which are smaller than the minimum value of the error threshold value in the error set as the real wind speed and the real wind direction of the measuring point.
10. An electronic device, comprising: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor when executing the computer program implementing a wind speed and direction inversion method according to any of claims 1-2 or a wind speed and direction inversion method according to any of claims 4-8.
11. A computer-readable storage medium, characterized in that a computer program is stored thereon, which when being executed by a processor implements the wind speed and direction inversion method according to any one of claims 1-2 or the wind speed and direction inversion method according to any one of claims 4-8.
CN202111071562.1A 2021-09-14 2021-09-14 Wind speed and direction inversion method, device, equipment and storage medium Active CN113514834B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111071562.1A CN113514834B (en) 2021-09-14 2021-09-14 Wind speed and direction inversion method, device, equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111071562.1A CN113514834B (en) 2021-09-14 2021-09-14 Wind speed and direction inversion method, device, equipment and storage medium

Publications (2)

Publication Number Publication Date
CN113514834A CN113514834A (en) 2021-10-19
CN113514834B true CN113514834B (en) 2022-03-08

Family

ID=78062993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111071562.1A Active CN113514834B (en) 2021-09-14 2021-09-14 Wind speed and direction inversion method, device, equipment and storage medium

Country Status (1)

Country Link
CN (1) CN113514834B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114324975B (en) * 2021-12-29 2024-01-09 航天科工火箭技术有限公司 Wind direction determining method, device, equipment and medium

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7365675B2 (en) * 2005-09-26 2008-04-29 The Boeing Company Measuring wind vectors remotely using airborne radar
DE102007058345A1 (en) * 2007-12-03 2009-06-04 Selex Sistemi Integrati Gmbh Method for determining composite data of weather radars in an overlapping region of the observation regions of at least two weather radars
CN107843895B (en) * 2017-10-20 2019-10-18 厦门市气象灾害防御技术中心(海峡气象开放实验室、厦门市避雷检测技术中心) A kind of Dual-Doppler weather radar dimensional wind inversion method
CN108107434B (en) * 2017-12-26 2020-06-09 厦门市气象灾害防御技术中心(海峡气象开放实验室厦门市避雷监测技术中心) Regional three-dimensional wind field picture splicing method based on double-Doppler radar inversion
CN112946657B (en) * 2021-02-03 2023-10-27 南京信息工程大学 Method for identifying ground wind field in strong convection weather

Also Published As

Publication number Publication date
CN113514834A (en) 2021-10-19

Similar Documents

Publication Publication Date Title
US10145933B2 (en) Angle determining system and method
US11815625B2 (en) Methods and devices for correcting underwater photon displacement and for depth sounding with single-photon Lidar
CN112083387B (en) Radar calibration method and device
CN113514834B (en) Wind speed and direction inversion method, device, equipment and storage medium
CN109195110B (en) Indoor positioning method based on hierarchical clustering technology and online extreme learning machine
CN112689775A (en) Radar point cloud clustering method and device
CN103116173A (en) Error test device for photoelectric tracking
CN111609849A (en) Multi-station rendezvous positioning method and system
JP2010146398A (en) Patrol support system, patrol support method and program
CN114879156A (en) Multi-target resolution method and device in radio frequency detection main beam and electronic equipment
CN115508775A (en) Using azimuth difference of incoming wave node positioning method for measurement
CN109041213A (en) Adaptive location modification method and device in a kind of room AP
Xie et al. A comparison of azimuth, misalignment, and orthogonality estimation methods in seismometer testing
CN113534130A (en) Multi-station radar multi-target data association method based on sight angle
CN114814779A (en) Buoy surge wave height observation data error evaluation method, system, equipment and medium
Sirola Exhaustive global grid search in computing receiver position from modular satellite range measurements
Wang et al. EMLR for multiple radars on a mobile platform
CN117075035B (en) Spin short baseline high-precision direction finding method, system, equipment, medium and terminal
CN115118366B (en) Multi-target resolution all-digital link modeling and checking method and device and electronic equipment
CN113960556B (en) Ship target RCS dynamic test data processing method and device
CN112763974B (en) Sensor deployment method, device and medium in positioning system
CN111880141B (en) Information source three-dimensional space positioning method of fully independent distributed linear array and related device
CN114966581A (en) Radar calibration method, device, equipment and storage medium
CN117518149A (en) Method, system, equipment and medium for calculating relative radial velocity
Ghassemi et al. A method for constructing the observer trajectory in bearings-only tracking of targets with a Markovian model

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