阅读说明：本技术 基于多普勒气象雷达的风场获取方法、装置、设备和介质 (Wind field acquisition method, device, equipment and medium based on Doppler meteorological radar ) 是由 周康明 申影影 于 2020-06-12 设计创作，主要内容包括：本申请涉及一种基于多普勒气象雷达的风场获取方法、装置、设备和介质。所述方法包括：获取低空下多普勒气象雷达的目标径向速度；根据所述目标径向速度,对雷达资料数据进行变分反演处理,得到低空下的目标反演风场；其中,所述雷达资料数据用于表征多普勒气象雷达探测区域的背景场、雨水含量、径向速度和噪音。采用本方法能够提高确定低空下的目标反演风场的精确性和可靠性,不仅能够用于天气过程的精细预报、边界层动力结构研究、风环境检测和邻近预报,而且也能够为机场附近天气过程的精细化预报,临近预报等提供依据,从而肯定了对航空、大型体育活动等的气象保障的重要现实意义。(The application relates to a wind field acquisition method, a device, equipment and a medium based on a Doppler meteorological radar. The method comprises the following steps: acquiring the target radial velocity of the Doppler meteorological radar under the low altitude; according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area. The method can improve the accuracy and reliability of determining the target inversion wind field at low altitude, can be used for fine forecasting of the weather process, boundary layer dynamic structure research, wind environment detection and proximity forecasting, and can provide basis for fine forecasting of the weather process near an airport, proximity forecasting and the like, thereby confirming the important practical significance of weather guarantee for aviation, large-scale sports activities and the like.)

1. A wind field acquisition method based on a Doppler meteorological radar is characterized by comprising the following steps:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

2. The method of claim 1, wherein the performing variational inversion processing on the radar data to obtain a target inversion wind field at low altitude comprises:

and performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of the low-altitude wind field, and taking the first target longitude and latitude wind field and the first target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

3. The method of claim 1, wherein the performing variational inversion processing on the radar data to obtain a target inversion wind field at low altitude comprises:

and performing four-dimensional variation inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

4. The method of claim 1, wherein the performing variational inversion processing on the radar data to obtain a target inversion wind field at low altitude comprises:

performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of a low-altitude wind field;

and performing four-dimensional variation inversion processing on radar data according to the first target longitude and latitude wind field and the first target height wind field to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field as the low-altitude target inversion wind field.

5. The method according to claim 2 or 4, wherein the three-dimensional variation inversion processing is performed on the radar data to obtain a first target longitude and latitude wind field and a first target altitude wind field of a wind field under low altitude, and the method comprises the following steps:

substituting the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation and a smoothness constraint equation to construct a target function;

minimizing the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of a wind field under low altitude;

and performing inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target height wind field.

6. The method of claim 4, wherein the obtaining of the second target longitude and latitude wind field and the second target altitude wind field of the wind field under low altitude by performing four-dimensional variational inversion processing on radar data according to the first target longitude and latitude wind field and the first target altitude wind field comprises:

judging whether the first target longitude and latitude wind field and the first target height wind field meet a preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment;

and when the first target longitude and latitude wind field and the first target height wind field are determined to meet the time discontinuity condition, performing four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

7. The method according to claim 4, wherein when the radar data includes merged data of radar networking, the four-dimensional variation inversion processing is performed on the radar data according to the first target longitude and latitude wind field and the first target altitude wind field to obtain a second target longitude and latitude wind field and a second target altitude wind field of a wind field under low altitude, and the method further comprises:

acquiring a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and acquiring reference wind speeds of a buoy in a Doppler meteorological radar detection range at different heights and a reference average root mean square error or a reference average absolute error of the reference wind speeds;

judging whether a first error difference value between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range or not;

and determining the first average root mean square error or the average absolute error, and performing four-dimensional variation inversion processing on radar data when an error difference value between the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

8. The method according to claim 3 or 4, wherein the four-dimensional variational inversion processing is performed on the radar data to obtain a second target longitude and latitude wind field and a second target altitude wind field of the wind field under the low altitude, and the method comprises the following steps:

constructing a value function according to a background field change constraint equation and a space-time smooth term constraint equation in the radar data;

and carrying out minimization processing on the value function to obtain a second target longitude and latitude wind field and a second target height wind field.

9. The method of claim 7, further comprising:

judging whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target height wind field and a second error difference value between the reference average root mean square error or the reference average absolute error of the reference wind speed are within a preset target difference value range or not;

when the second error difference value is determined to be within the target difference value range, outputting first prompt information with higher consistency level between the target inversion wind field and the reference wind speed;

and when the second error difference is determined not to be in the target difference range, outputting second prompt information that the target inversion wind field is invalid.

10. The method of claim 7, wherein the obtaining the reference wind speed of the buoy at different heights within the detection range of the Doppler meteorological radar comprises calculating the reference wind speed at the height z of the buoy by using the following formula:

wherein, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Denotes the coefficient of friction, k denotes the von Karman constant,represents the stable correction function, and l represents the Monin-Obukhov length.

11. The method according to claim 1, wherein before the step of performing variational inversion processing on the radar data to obtain the target inversion wind field at low altitude, the method comprises the following steps:

judging the magnitude relation between the target radial speed and a preset speed threshold;

and when the target radial velocity is determined to be greater than or equal to the velocity threshold, executing the variational inversion processing on the radar data to obtain a target inversion wind field under low altitude.

12. The method according to claim 1, wherein the number of the Doppler weather radar is at least two, and when a radar network is formed based on at least two Doppler weather radars, the obtaining of the target radial velocity of the Doppler weather radar under low altitude comprises:

acquiring the target radial velocity of each Doppler meteorological radar;

and performing composite operation on the target radial speed by adopting a quadrilateral rule to obtain the target composite radial speed of the low-altitude radar network.

13. A wind field acquisition device based on doppler meteorological radar, the device comprising:

the acquisition module is used for acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

the processing module is used for carrying out variation inversion processing on the radar data according to the target radial velocity to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

14. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the method of any of claims 1 to 12.

15. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 12.

Technical Field

The application relates to the technical field of aeronautical weather, in particular to a wind field acquisition method, a device, equipment and a medium based on a Doppler weather radar.

Background

The triggering and development of the wind and the precipitation process are closely related, and particularly, a low-layer amplitude line, a tangent line and the like are not only closely related to the generation of the convection cloud, but also one of main meteorological factors for forecasting the development of the convection cloud in advance. Therefore, in the field of meteorological guarantees for large activities like the olympic games and airports, the demand for fine observation of wind fields is becoming stronger and stronger.

The traditional technology comprises three parts of radar image pretreatment, wind direction measurement and wind speed measurement, wherein in a wind direction measurement index, image gradient, gray level and smooth items are organically combined, and the proportion of the three items is adjusted through a proportionality coefficient to establish a model suitable for sea surface wind field characteristics; in the wind speed measurement indexes, when the radar measures independently, the NRCS, the actually measured wind direction and the SNR are used as BP network input; in the wind speed measurement index, the sea surface wind speed of the marine radar is measured by taking the sea air boundary layer parameter as the additional input of the BP network.

However, the wind speed measurement in the conventional technology requires radar image processing and scaling factor adjustment, resulting in poor accuracy of the measured wind speed.

Disclosure of Invention

In view of the above, it is necessary to provide a doppler meteorological radar-based wind field acquisition method, apparatus, computer device and storage medium capable of measuring the accuracy of a low-altitude wind field.

A method of doppler meteorological radar based wind field acquisition, the method comprising:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

In one embodiment, the method further comprises the following steps:

and performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of the low-altitude wind field, and taking the first target longitude and latitude wind field and the first target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the method further comprises the following steps:

and performing four-dimensional variation inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the method further comprises the following steps:

performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of a low-altitude wind field;

and performing four-dimensional variation inversion processing on radar data according to the first target longitude and latitude wind field and the first target height wind field to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field as the low-altitude target inversion wind field.

In one embodiment, the method further comprises the following steps:

substituting the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation and a smoothness constraint equation to construct a target function;

minimizing the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of a wind field under low altitude;

and performing inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target height wind field.

In one embodiment, the method further comprises the following steps:

judging whether the first target longitude and latitude wind field and the first target height wind field meet a preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment;

and when the first target longitude and latitude wind field and the first target height wind field are determined to meet the time discontinuity condition, performing four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the method further comprises the following steps:

acquiring a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and acquiring reference wind speeds of a buoy in a Doppler meteorological radar detection range at different heights and a reference average root mean square error or a reference average absolute error of the reference wind speeds;

judging whether a first error difference value between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range or not;

and determining the first average root mean square error or the average absolute error, and performing four-dimensional variation inversion processing on radar data when an error difference value between the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the method further comprises the following steps:

constructing a value function according to a background field change constraint equation and a space-time smooth term constraint equation in the radar data;

and carrying out minimization processing on the value function to obtain a second target longitude and latitude wind field and a second target height wind field.

In one embodiment, the method further comprises the following steps:

judging whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target height wind field and a second error difference value between the reference average root mean square error or the reference average absolute error of the reference wind speed are within a preset target difference value range or not;

when the second error difference value is determined to be within the target difference value range, outputting first prompt information with higher consistency level between the target inversion wind field and the reference wind speed;

and when the second error difference is determined not to be in the target difference range, outputting second prompt information that the target inversion wind field is invalid.

In one embodiment, the method further comprises the following steps:

the acquiring of the reference wind speeds of the buoy in the detection range of the Doppler meteorological radar at different heights comprises calculating the reference wind speed when the height of the buoy is z by adopting the following formula:

wherein, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Expressing the coefficient of friction, k expressing the von Karman constant,

represents the stability correction function, l represents the length of the Morin-obuHough (Monin-Obukhov).

In one embodiment, the method further comprises the following steps:

judging the magnitude relation between the target radial speed and a preset speed threshold;

and when the target radial velocity is determined to be greater than or equal to the velocity threshold, executing the variational inversion processing on the radar data to obtain a target inversion wind field under low altitude.

In one embodiment, the number of the doppler meteorological radars is at least two, and when a radar network is formed based on at least two doppler meteorological radars, the acquiring the target radial velocity of the doppler meteorological radar under the low altitude includes:

acquiring the target radial velocity of each Doppler meteorological radar;

and performing composite operation on the target radial speed by adopting a quadrilateral rule to obtain the target composite radial speed of the low-altitude radar network.

A doppler meteorological radar based wind field acquisition apparatus, the apparatus comprising:

the acquisition module is used for acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

the processing module is used for carrying out variation inversion processing on the radar data according to the target radial velocity to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

According to the wind field obtaining method and device based on the Doppler meteorological radar, the computer equipment and the storage medium, the target radial velocity of the Doppler meteorological radar under the low altitude is obtained firstly, so that the variation inversion processing is determined to be carried out on radar data according to the target radial velocity, and the accuracy and reliability of determining the target inversion wind field under the low altitude are improved; moreover, the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area, so that a target inversion wind field with a small detection range can be determined when the radar data of a single Doppler meteorological radar are obtained, and a target inversion wind field with a large detection range can be determined when the radar data of a radar networking are obtained; furthermore, the target inversion wind field is obtained by performing the variation inversion processing on the radar data, so that the target inversion wind field can be obtained by using the three-dimensional variation inversion processing and the four-dimensional inversion processing on the radar data, the detailed information of the three-dimensional wind field and the fine structure of the three-dimensional wind field can be provided according to the obtained target inversion wind field, the target inversion wind field can be used for fine forecasting of the weather process, boundary layer dynamic structure research, wind environment detection and proximity forecasting, and the basis can be provided for fine forecasting of the weather process near an airport, proximity forecasting and the like, so that the important practical significance of weather guarantee of aviation, large-scale sports activities and the like is determined.

Drawings

FIG. 1 is a schematic flow diagram of a wind farm method based on Doppler weather radar in one embodiment;

FIG. 2 is a schematic diagram illustrating consistency verification of results obtained after radar data is processed by a variational inversion wind field in one embodiment;

FIG. 3A is a schematic flow chart illustrating a wind farm method based on Doppler weather radar in accordance with yet another embodiment;

FIG. 3B is a schematic flow chart illustrating a wind field method based on Doppler weather radar in accordance with yet another embodiment;

FIG. 3C is a schematic flow chart illustrating a wind field method based on Doppler weather radar in accordance with yet another embodiment;

FIG. 4 is a schematic flow chart of a wind field method based on Doppler meteorological radar in yet another embodiment;

FIG. 5 is a schematic flow chart of a wind field method based on Doppler meteorological radar in yet another embodiment;

FIG. 6 is a schematic flow chart of a wind field method based on Doppler meteorological radar in yet another embodiment;

FIG. 7 is a schematic flow chart of a wind field method based on Doppler meteorological radar in yet another embodiment;

FIG. 8 is a block diagram of a wind farm apparatus based on Doppler weather radar in a further embodiment;

FIG. 9 is a diagram illustrating an internal structure of a computer device according to an embodiment.

Detailed Description

In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.

According to the method for acquiring the low altitude wind shear region, an execution subject can be an acquiring device of the low altitude wind shear region, and the acquiring device of the low altitude wind shear region can be implemented as part or all of computer equipment in a software, hardware or combination of the software and the hardware. Optionally, the computer device may be an electronic device with a camera function, such as a Personal Computer (PC), a portable device, a notebook computer, a smart phone, a tablet computer, a portable wearable device, and the like, for example, a tablet computer, a mobile phone, and the like, and the specific form of the computer device is not limited in the embodiment of the present application.

It should be noted that the execution subject of the method embodiments described below may be part or all of the computer device described above. The following method embodiments are described by taking the execution subject as the computer device as an example.

In one embodiment, as shown in fig. 1, there is provided a wind field acquisition method based on a doppler meteorological radar, comprising the following steps:

and step S11, acquiring the target radial velocity of the Doppler meteorological radar under the low altitude.

Wherein the low altitude may include a flight area of 1000 meters or less. Alternatively, the low altitude may be 600 meters or less, and the target radial velocity includes one of a radial velocity of a single doppler meteorological radar and a target composite radial velocity corresponding to at least two doppler meteorological radars.

Specifically, the computer device may obtain the target radial velocity according to the number of the doppler meteorological radars, that is, when it is determined that the number of the doppler meteorological radars is one, the radial velocity of the doppler meteorological radar may be used as the target radial velocity of the doppler meteorological radar in the low altitude; when it is determined that the number of the Doppler meteorological radars is at least two and a radar networking is formed based on the at least two Doppler meteorological radars, the target radial speed of each Doppler meteorological radar can be obtained first, and then the quadrangle rule is adopted to perform composite operation on the obtained at least two target radial speeds, so that the target composite radial speed of the low-altitude radar networking is obtained. Optionally, the quadrilateral rules may comprise parallelogram rules.

Step S12, according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

Specifically, the variation inversion process may include at least one of a three-dimensional variation inversion process and a four-dimensional variation inversion process, the radar profile data may include profile data detected by a doppler meteorological radar in an X-band, and when the number of the doppler meteorological radars is one, the radar profile data may include profile data of the doppler meteorological radar; when the number of the doppler meteorological radars is at least two, the radar data may include combined data obtained by combining the data of each doppler meteorological radar.

In the actual processing process, when the computer equipment acquires the target radial velocity of the Doppler meteorological radar under the low altitude, three-dimensional variational inversion processing can be performed on radar data by adopting a three-dimensional variational (3DVAR) inversion wind field algorithm to obtain a first target longitude and latitude wind field and a first target height wind field under the low altitude; and a four-dimensional variation (4DVAR) inversion wind field algorithm can be adopted to perform four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field under low altitude.

Then, after consistency verification is performed on the first target longitude and latitude wind field and the first target altitude wind field, and the second target longitude and latitude wind field and the second target altitude wind field, a consistency effect graph as shown in fig. 2 is obtained. As can be seen from fig. 2, the refined wind field inversion can be realized by performing the 3DVAR inversion wind field algorithm and the 4DVAR inversion wind field algorithm on the radar data, and the consistency of the inversion wind field is better than that of weather with smaller wind speed. Therefore, the target inversion wind field under the low altitude may include a first target longitude and latitude wind field and a first target altitude wind field of the low altitude wind field, and may also include a second target longitude and latitude wind field and a second target altitude wind field of the low altitude wind field. Alternatively, the weather in which the wind speed is small may be a weather in which the wind speed is lower than 8.2 m/s.

In the actual processing process, the computer equipment can also perform quantitative analysis on parameters such as correlation coefficients, root mean square errors, average absolute errors and the like of the first target longitude and latitude wind field and the first target altitude wind field, and the second target longitude and latitude wind field and the second target altitude wind field to obtain 3DVAR and 4DVAR wind field inversion radar data, so that refined wind field inversion can be realized, and wind direction changes can be reflected.

According to the wind field acquisition method based on the Doppler meteorological radar, the target radial velocity of the Doppler meteorological radar under the low altitude is firstly acquired, so that the variation inversion processing is conveniently determined to be carried out on radar data according to the target radial velocity, and the accuracy and reliability of determining the target inversion wind field under the low altitude are improved; moreover, the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area, so that a target inversion wind field with a small detection range can be determined when the radar data of a single Doppler meteorological radar are obtained, and a target inversion wind field with a large detection range can be determined when the radar data of a radar networking are obtained; furthermore, the target inversion wind field is obtained by performing the variation inversion processing on the radar data, so that the target inversion wind field can be obtained by using the three-dimensional variation inversion processing and the four-dimensional inversion processing on the radar data, the detailed information of the three-dimensional wind field and the fine structure of the three-dimensional wind field can be provided according to the obtained target inversion wind field, the target inversion wind field can be used for fine forecasting of the weather process, boundary layer dynamic structure research, wind environment detection and proximity forecasting, and the basis can be provided for fine forecasting of the weather process near an airport, proximity forecasting and the like, so that the important practical significance of weather guarantee of aviation, large-scale sports activities and the like is determined.

In one embodiment, as shown in fig. 3A, step S12 includes:

step S121, three-dimensional variation inversion processing is carried out on the radar data, and a first target longitude and latitude wind field and a first target height wind field of the low-altitude wind field are obtained.

Specifically, as shown in fig. 3B, the computer device may obtain the first target longitude and latitude wind field and the first target altitude wind field by using the following sub-steps:

and S1211, substituting the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation and a smoothness constraint equation to construct an objective function.

Wherein the objective function is: j (μ, ν, ω) ═ J_B+J_o+J_E+J_C+J_P(1)

In the above formulas (1) to (6), a three-dimensional wind field coordinate system is established with the central point in the radar effective detection area of the doppler meteorological radar as the origin and with the longitude and latitude directions and the height direction, mu, v, and omega represent the longitude and latitude wind field and the height wind field to be inverted, and J_BRepresenting the ambient field constraint equation (J when the ambient field constraint is not considered)_B＝0)，J_oRepresenting an observation fitting constraint equation, J_EExpressing the conservation equation of rainwater content, J_CRepresenting the continuous wind field constraint equation, J_PRepresenting the smoothness constraint equation, u^xyRepresenting the wind field longitude, v, of the effective detection area of the radar in the xy plane of a Cartesian coordinate system^xyRepresenting the wind field latitude, w, of the effective detection area of the radar in the xy plane of a Cartesian coordinate system^xyRepresenting the wind field height u of the effective detection area of the radar in the xy plane of a Cartesian coordinate system^BRepresenting the wind field longitude, v, of the background field in the radar effective detection area in the three-dimensional wind field coordinate system^BRepresenting the wind field latitude, w, of the background field in the radar effective detection area in the three-dimensional wind field coordinate system^BRepresenting the wind field height W of the background field in the radar effective detection area in the three-dimensional wind field coordinate system_uBA wind field longitude weight representing a background field within the radar active detection area; w_vBA wind field latitude weight representing a background field within the radar active detection area; w_wBA wind field height weight, W, representing the background field in the effective detection area of the radar_ERepresenting a weight of rain water content in the radar active detection area, E representing an increment of rain water content in the radar active detection area, W_rA weight matrix representing Doppler meteorological radar, when the number of Doppler meteorological radars is at least two, W_rThe value size of the Doppler weather radar is determined by the relative position of each Doppler weather radar and hardware information; v. of_rTo representA wind field latitude weight for each grid in the radar active detection area,representing a wind field latitude weight of each grid point in the radar effective detection area; w_cRepresenting a continuous wind field weight factor, (i, j, k) representing a grid point number divided into a plurality of three-dimensional grid points in the radar effective detection area, d representing a horizontal grid distance, d²Representing the horizontal Laplacian, W_puRepresents the horizontal smoothness constraint weight, W_pvRepresents the vertical smoothness constraint weight, W_pwRepresenting the height direction smoothness constraint weight. Optionally, the radar effective detection area includes one of a detection range of a single doppler weather radar and a combined detection range obtained by combining detection ranges of at least two doppler weather radars.

In the actual process, the wind field longitude u^xyLatitude v of wind field^xyHeight w of wind field^xyWind field latitude weight v_rThe rain water content increment E can be calculated using the following formula:

u^xy＝0.5u_i,j+0.125(u_i+1,j+u_i-1,j+u_i,j+1+u_i,j-1) (7)

v^xy＝0.5v_i,j+0.125(v_i+1,j+v_i-1,j+v_i,j+1+v_i,j-1) (8)

w^xy＝0.5w_i,j+0.125(w_i+1,j+w_i-1,j+w_i,j+1+w_i,j-1) (9)

M＝0.01Z^0.5，

in the above formulae (6) to (11), u_i,jRepresenting the wind field longitude, u of the (i, j) th grid point in the radar active detection area_i+1,jRepresenting the wind field longitude, u of the (i +1, j) th grid point in the radar active detection area_i-1,jRepresenting the wind field longitude, u, of the (i-1, j) th grid point in the radar active detection area_i,j+1Representing a wind field longitude, u, of an (i, j +1) th grid point in the radar active detection area_i,j-1Representing a wind field longitude of an (i, j-1) th grid point in the radar active detection area; v. of_i,jRepresenting the wind field latitude, v, of the (i, j) th grid point in the radar active detection area_i+1,jRepresenting the wind field latitude, v, of the (i +1, j) th grid point in the radar active detection area_i-1,jRepresenting the wind field latitude, v, of the (i-1, j) th grid point in the radar active detection area_i,j+1Representing the wind field latitude, v, of the (i, j +1) th grid point in the radar active detection area_i,j-1Representing the wind field latitude of the (i, j-1) th grid point in the radar effective detection area; w is a_i,jRepresenting the wind field height, w, of the (i, j) th grid point in the radar active detection area_i+1,jRepresenting the wind field height, w, of the (i +1, j) th grid point in the radar active detection area_i-1,jRepresenting the wind field height, w, of the (i-1, j) th grid point in the radar active detection area_i,j+1Representing the wind field height, w, of the (i, j +1) th grid point in the radar active detection area_i,j-1Representing the wind field height of an (i, j-1) th grid point in the radar effective detection area; (x, y, z) represents the coordinate position of each grid point in the radar active detection area, (x)_r,y_r,z_r) Coordinates, w, of the antenna position of the Doppler weather radar_TRepresents an initial height of each grid point in the radar effective detection area, t represents a determination time, Z represents a height under a geodetic coordinate system,indicating a determined altitudeRain water content, rho₀Indicating the amount of rainfall at sea level.

Step S1212, performing minimization processing on the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of the wind field in the low altitude.

Specifically, when the computer device performs minimization processing on the objective function, the objective function may be solved by adopting a gradient descent method or a derivation method, and the obtained solution is expressed as a second-order lux polynomial expansion, so as to obtain the three-dimensional target longitude and latitude wind field u (x, y, z), v (x, y, z) and the three-dimensional target height wind field w (x, y, z); wherein:

in the above formula, P_nx(x) Legendre basis function, P, representing the x-direction at time n_ny(y) Legendre basis function in the y-direction at time n, P_nz(z) Legendre basis function in z direction at time n, a_nx,ny,nzThe expansion coefficients of the second order lux polynomial expansion are shown. Alternatively, n is 1,2 when radar profile data of a doppler meteorological radar is used three times per inversion.

In the actual processing process, the computer device may determine a minimization processing mode according to the order of the objective function, that is, when the order of the objective function is lower than the second order, the objective function may be differentiated; when the order of the target function is larger than the second order, the target function can be solved through a gradient descent method so as to obtain the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field.

Step S1213, performing inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target altitude wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target altitude wind field.

Specifically, when obtaining the three-dimensional target longitude and latitude wind field u (x, y, z), v (x, y, z) and the three-dimensional target height wind field w (x, y, z), the computer device may further perform inversion processing on the three-dimensional target longitude and latitude wind field u (x, y, z), v (x, y, z) and the three-dimensional target height wind field w (x, y, z) in a rectangular coordinate system to obtain the first target longitude and latitude wind field u (x, y), v (x, y) and the first target height wind field w (x, y).

In an actual processing process, after performing inversion processing on the three-dimensional target longitude and latitude wind field u (x, y, z), v (x, y, z) and the three-dimensional target height wind field w (x, y, z) in a rectangular coordinate system, the obtained first target longitude and latitude wind field may be u (x, z) and v (x, z), and the first target height wind field may be w (x, z); the first target longitude and latitude wind field can also be u (y, z) and v (y, z), and the first target height wind field can also be w (x, z).

And S122, performing four-dimensional variation inversion processing on radar data according to the first target longitude and latitude wind field and the first target height wind field to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field as target inversion wind fields under low altitude.

Specifically, the computer device may calculate a quantized value of consistency between the first target longitude and latitude wind field and the first target altitude wind field and a reference wind speed and a reference wind direction of a buoy in the radar effective detection area according to the first target longitude and latitude wind field and the first target altitude wind field, and further determine that the quantized value of consistency is less than a preset quantized threshold value of consistency, determine that the first target longitude and latitude wind field and the first target altitude wind field have low consistency with the reference wind speed and the reference wind direction, and then the computer may perform four-dimensional variation inversion processing on radar data; optionally, the consistency quantization threshold may take a value of 8.2 m/s.

As shown in fig. 3C, the computer device may obtain the second target longitude and latitude wind field and the second target altitude wind field by using the following sub-steps:

and step S1221, constructing a cost function according to a background field change constraint equation and a space-time smooth term constraint equation in the radar data.

Specifically, the computer device may substitute the radar data into a continuity equation, a momentum equation, an energy equation, and a thermodynamic equation, and then solve the data to obtain a background field change constraint equation, and then construct a cost function according to the background field change constraint equation, where:

in the above formula, J_B'Representing a constraint equation of the variation of the background field, and the constraint equation of the variation of the background field is used for measuring the difference between the mode analysis value of the initial moment of the assimilation window in the effective detection area of the radar and the background field η_vWeighting factor representing radial velocity, η calculated for convenience_v＝1；V_riRepresents the radial velocity of the doppler meteorological radar forecast,represents the radial velocity of the Doppler meteorological radar observation; j. the design is a square_PRepresenting a noise equation, also referred to as a penalty term or a spatio-temporal smoothing term, which may reduce noise and may also speed up the rate of minimization convergence of the cost function.

Step S1222, minimize the cost function to obtain the second target longitude and latitude wind field and the second target altitude wind field.

Specifically, after obtaining the cost function, the computer device may solve the cost function, that is, perform minimization on the cost function and then interpolate the value function into a cartesian coordinate system, so as to obtain the second target longitude and latitude wind field and the second target height wind field; and in the process of solving the longitude and latitude wind field and the altitude wind field of the second target, the background field added into the effective radar detection area is an atmospheric average field.

In the actual processing process, the computer device may determine a minimization processing mode according to the order of the cost function, that is, when the order of the cost function is lower than the second order, the derivation may be performed on the cost function; when the order of the cost function is larger than the second order, the cost function can be solved through a gradient descent method so as to obtain the second target longitude and latitude wind field and the second target altitude wind field.

In this embodiment, the computer device determines to perform four-dimensional variation inversion processing on radar data through the first target longitude and latitude wind field and the first target altitude wind field obtained after three-dimensional variation inversion processing, so as to obtain a target inversion wind field in low altitude, and not only can the first target longitude and latitude wind field and the first target altitude wind field be used as a refined wind field inversion result to reflect wind direction changes on the premise that the requirement for consistency index is low, but also the second target longitude and latitude wind field and the second target altitude wind field obtained after four-dimensional variation inversion processing can be used as the target inversion wind field on the premise that the requirement for consistency index is high, so that the flexibility and reliability of the target inversion wind field can be obtained according to different requirements for consistency index.

In one embodiment, as shown in fig. 4, the performing, in step S122, four-dimensional variational inversion processing on radar data according to the first target longitude and latitude wind field and the first target altitude wind field to obtain a second target longitude and latitude wind field and a second target altitude wind field of the wind field under low altitude includes:

step S21, judging whether the first target longitude and latitude wind field and the first target altitude wind field meet the preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment.

Specifically, when the first target longitude and latitude wind field and the first target altitude wind field are obtained through three-dimensional variation inversion processing, the computer device may further obtain a first corresponding relationship between the target longitude and latitude wind field at each time in the first target longitude and latitude wind field and the radar data at each time, and a second corresponding relationship between the radar data at each time and the target altitude wind field in the first target altitude wind field, and then, when it is determined that the current target longitude and latitude at the current time is only obtained for the current radar data at the current time according to the first corresponding relationship and the second corresponding relationship, the process proceeds to step S22.

And step S22, when the first target longitude and latitude wind field and the first target altitude wind field are determined to meet the time discontinuity condition, performing four-dimensional variational inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target altitude wind field of the low-altitude wind field.

Specifically, the computer device determines that the first target longitude and latitude wind field and the first target altitude wind field meet the time discontinuity condition, which may include that the current target longitude and latitude at the current time is only obtained for the current radar data at the current time, that is, the first target longitude and latitude wind field and the first target altitude wind field are independent in time, and at this time, four-dimensional variation inversion processing may be performed on the radar data to obtain the second target longitude and latitude wind field and the second target altitude wind field which are continuous in time.

In this embodiment, the computer device determines that the radar data is subjected to four-dimensional variation inversion processing by that the first target longitude and latitude wind field and the first target altitude wind field meet a preset time discontinuity condition, so as to obtain the second target longitude and latitude wind field and the second target altitude wind field which are continuous in time, thereby embodying necessity and effectiveness of obtaining the second target longitude and latitude wind field and the second target altitude wind field, and improving accuracy and reliability of the second target longitude and latitude wind field and the second target altitude wind field.

In an embodiment, when the radar data includes merged data of a radar network, as shown in fig. 5, in step S122, performing four-dimensional variation inversion processing on the radar data according to the first target longitude and latitude wind field and the first target altitude wind field to obtain a second target longitude and latitude wind field and a second target altitude wind field of a wind field under low altitude, includes:

step S31, obtaining a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and obtaining reference wind speeds of buoys in the detection range of the Doppler meteorological radar at different heights and a reference average root mean square error or a reference average absolute error of the reference wind speeds.

Specifically, when the computer device obtains the reference wind speeds of the buoy in the detection range of the doppler meteorological radar at different heights, the reference wind speed when the height of the buoy is z can be calculated by adopting the following formula:

in the above formula, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Denotes the coefficient of friction, k denotes the von Karman constant,represents the stable correction function, and l represents the Monin-Obukhov length.

In an actual processing process, when obtaining the first target longitude and latitude wind field and the first target altitude wind field, the computer device may further determine a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field by using an existing average root mean square error or average absolute error calculation method, so as to determine necessity of performing four-dimensional variation inversion processing on radar data according to a magnitude relation between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error.

And when the number of the Doppler meteorological radars is at least two, the acquired radar data can avoid the defect of incomplete data caused by the mirror image blind area of a single Doppler meteorological radar, so that the integrity and the accuracy of radar combined data obtained by combining the radar data of at least two Doppler meteorological radars are ensured.

Step S32, it is determined whether the first average root mean square error or the first average absolute error, and the first error difference between the reference average root mean square error or the reference average absolute error are within a preset difference range.

The preset difference range can be set according to the accuracy requirement expected by the user and fed back by the client. Alternatively, when the first error difference value is a difference between the first average root mean square error and the reference average root mean square error, the preset difference range may be (1, 3); when the first error difference is a difference between the first average absolute error and the reference average absolute error, the preset difference range may be (1, 2.4).

Specifically, when the computer device obtains the first average root mean square error or the first average absolute error, it may further determine whether a first error difference between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference range, so as to embody a disadvantage and a disadvantage of a three-dimensional variation processing algorithm for processing radar data, and also embody a necessity of performing four-dimensional variation inversion processing on the radar data.

And step S33, when the first average root mean square error or the average absolute error is determined and the error difference between the reference average root mean square error or the reference average absolute error is within a preset difference range, performing four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

Specifically, when the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error are within a preset difference range, the computer device may determine that the accuracy of the first target longitude and latitude wind field and the first target altitude wind field obtained by processing radar data by the three-dimensional variation processing algorithm does not meet the user expected accuracy requirement fed back by the client, and further perform four-dimensional variation inversion processing on the radar data, so as to obtain the second target longitude and latitude wind field and the second target altitude wind field meeting the user expected accuracy requirement.

In this embodiment, the computer device obtains the first mean square root error or the first mean absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and the difference between the reference mean square error or the reference mean absolute error of the reference wind speed of the buoy in the detection range of the Doppler meteorological radar at different heights, so as to realize the purpose of carrying out four-dimensional variation inversion processing on radar data when the error difference value is within the preset difference value range, therefore, the purpose of executing four-dimensional variational inversion processing radar data when the accuracy of the first target longitude and latitude wind field and the first target altitude wind field does not meet the accuracy requirement of user expectation fed back by a client can be achieved, and the accuracy and reliability of the second target longitude and latitude wind field and the second target altitude wind field are greatly improved.

In one embodiment, as shown in fig. 6, the method further comprises:

step S41, determining whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target altitude wind field, and a second error difference between the reference average root mean square error or the reference average absolute error of the reference wind speed, are within a preset target difference range.

The target difference range can be set according to the accuracy requirement expected by the user and fed back by the client. Alternatively, when the second error difference value is a difference between the second average root mean square error and the reference average root mean square error, the target difference value range may be (0, 2); when the second error difference is a difference between the second average absolute error and the reference average absolute error, the target difference range may be (0, 1.4).

Specifically, when obtaining the second average root mean square error or the second average absolute error, the computer device may further determine whether a second error difference between the second average root mean square error or the second average absolute error and the reference average root mean square error or the reference average absolute error is within a target difference range, so as to determine whether the second target longitude and latitude wind field and the second target altitude wind field are target inversion wind fields meeting the accuracy desired by the user.

And step S42, when the second error difference is determined to be within the target difference range, outputting first prompt information with higher consistency level between the target inversion wind field and the reference wind speed.

Specifically, when the second error difference is within the target difference range, the computer device may determine that the second target longitude and latitude wind field and the second target altitude wind field are target inversion wind fields meeting the accuracy desired by the user, and may further determine that the second target longitude and latitude wind field and the second target altitude wind field have higher consistency with a floating mark true value within a detection range of the doppler meteorological radar, and may output first prompt information having a higher consistency level between the target inversion wind field and the reference wind speed at this time, so as to embody the advantage of the four-dimensional variational inversion algorithm over the three-dimensional variational inversion algorithm in processing radar data.

And step S43, when the second error difference is determined not to be within the target difference range, outputting second prompt information that the target inversion wind field is invalid.

Specifically, when the second error difference is not within the target difference range, the computer device may determine that the second target longitude and latitude wind field and the second target altitude wind field acquired this time are not target inversion wind fields meeting the accuracy desired by the user, so that second prompt information indicating that the target inversion wind fields are invalid may be output to prompt the client to check the integrity of the radar data and/or the accuracy of the four-dimensional variation inversion processing algorithm corresponding to the user, thereby achieving the result validity and accuracy after the radar data is subsequently processed by the four-dimensional variation inversion processing algorithm.

In this embodiment, the computer device outputs first prompt information that the consistency level between the target inversion wind field and the reference wind speed is higher when the second error difference is within the target difference range, and outputs second prompt information that the target inversion wind field is invalid when the second error difference is not within the target difference range, by obtaining a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target altitude wind field, and a second error difference between reference average root mean square errors or reference average absolute errors of the buoy within the doppler meteorological radar detection range, so that flexibility and reliability of processing radar data by a four-dimensional variable division inversion algorithm can be realized.

In one embodiment, as shown in fig. 7, before the step of performing variation inversion processing on the radar data to obtain the target inversion wind field at low altitude in step S12, the method further includes:

step S51, determining a magnitude relationship between the target radial velocity and a predetermined velocity threshold.

In particular, the speed threshold may be used to characterize the wind field strength within the radar active detection area of the doppler meteorological radar. Alternatively, the speed threshold may be 8.2 m/s.

Specifically, when acquiring the target radial velocity of the doppler meteorological radar in the low altitude, the computer device may further determine a magnitude relationship between the target radial velocity and the velocity threshold value, so as to determine whether the wind field intensity in the radar effective detection area is a large wind field intensity.

And step S52, when the target radial velocity is determined to be greater than or equal to the velocity threshold value, executing the variational inversion processing on the radar data to obtain a target inversion wind field under low altitude.

Specifically, when the target radial velocity is greater than or equal to the velocity threshold, the computer device may determine that the wind field intensity in the radar effective detection area is a large wind field intensity, and at this time, may further perform a step of performing variation inversion processing on radar data, so as to obtain a target inversion wind field at a low altitude.

In the actual processing process, when the target radial velocity is smaller than the velocity threshold, it can be determined that the wind field intensity in the radar effective detection area is smaller wind field intensity, and correlation between low-altitude generated low-layer amplitude resultant lines, shear lines and the like and convection clouds is small under the smaller intensity, so that wind field inversion processing operation is not performed on the radar effective detection area corresponding to the smaller wind field intensity.

In this embodiment, the computer device executes the purpose of performing the variational inversion processing on the radar data to obtain the target inversion wind field at low altitude based on the condition that the determined target radial velocity is greater than or equal to the preset velocity threshold, so as to implement the purpose of performing the variational inversion processing on the radar data when the wind field intensity in the radar effective detection area of the doppler meteorological radar is large, thereby improving the pertinence and reliability of the variational inversion processing on the radar data.

It should be understood that although the various steps in the flowcharts of fig. 1, 3A, 3B, 3C, 4-7 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 1, 3A, 3B, 3C, and 4-7 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed in turn or alternately with other steps or at least some of the other steps.

In one embodiment, as shown in fig. 8, there is provided a wind field acquisition apparatus based on a doppler meteorological radar, including: an acquisition module 11 and a processing module 12, wherein:

and the obtaining module 11 is configured to obtain a target radial velocity of the doppler meteorological radar in the low altitude.

The processing module 12 is configured to perform variation inversion processing on the radar data according to the target radial velocity to obtain a target inversion wind field at low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

The processing module 12 may also be configured to perform three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target altitude wind field of the low-altitude wind field, and use the first target longitude and latitude wind field and the first target altitude wind field of the low-altitude wind field as target inversion wind fields in low altitudes.

The processing module 12 may also be configured to perform four-dimensional variation inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target altitude wind field of the low-altitude wind field, and use the second target longitude and latitude wind field and the second target altitude wind field of the low-altitude wind field as target inversion wind fields in the low altitude.

The processing module 12 may further specifically include: a first processing unit and a second processing unit.

Specifically, the first processing unit may be configured to perform three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target altitude wind field of the low-altitude wind field; the second processing unit may be configured to perform four-dimensional variation inversion processing on the radar data according to the first target longitude and latitude wind field and the first target altitude wind field to obtain a second target longitude and latitude wind field and a second target altitude wind field of the low altitude wind field, and use the second target longitude and latitude wind field and the second target altitude wind field as target inversion wind fields in the low altitude.

The first processing unit may further specifically include: a first building subunit, a first processing subunit and an inversion subunit.

Specifically, the first constructing subunit may be configured to substitute the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation, and a smoothness constraint equation to construct an objective function; the first processing subunit is configured to perform minimization processing on the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of a low-altitude wind field; and the inversion subunit is configured to perform inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target height wind field.

The second processing unit may further specifically include: the device comprises a first judgment subunit and a second processing subunit.

Specifically, the first determining subunit may be configured to determine whether the first target longitude and latitude wind field and the first target altitude wind field satisfy a preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment; the second processing subunit may be configured to, when it is determined that the first target longitude and latitude wind field and the first target altitude wind field satisfy the time discontinuity condition, perform four-dimensional variational inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target altitude wind field of the low-altitude wind field.

The second processing unit may further specifically include: the device comprises an acquisition subunit, a second judgment subunit and a third processing subunit. Specifically, the obtaining subunit may be configured to obtain a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and obtain a reference wind speed of a buoy in a detection range of the doppler meteorological radar at different altitudes and a reference average root mean square error or a reference average absolute error of the reference wind speed; the second judging subunit is configured to judge whether a first error difference between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference range; the third processing subunit may be configured to determine the first average root mean square error or the average absolute error, and perform four-dimensional variation inversion processing on the radar data when an error difference between the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference range, to obtain a second target longitude and latitude wind field and a second target altitude wind field of the low-altitude wind field.

The second processing unit may further specifically include: a second building subunit and a fourth processing subunit.

Specifically, the second constructing subunit may be configured to construct a cost function according to a background field change constraint equation and a space-time smoothing term constraint equation in the radar data; the fourth processing subunit may be configured to perform minimization processing on the cost function to obtain the second target longitude and latitude wind field and the second target altitude wind field.

The processing module 12 may further include: the device comprises a judging unit, a first output unit and a second output unit.

Specifically, the determining unit may be configured to determine whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target altitude wind field, and a second error difference between the reference average root mean square error or the reference average absolute error of the reference wind speed are within a preset target difference range; the first output unit may be configured to output first prompt information that the consistency level between the target inversion wind field and the reference wind speed is higher when it is determined that the second error difference is within the target difference range; and the second output unit may be configured to output second prompt information that the target inversion wind field is invalid when it is determined that the second error difference is not within the target difference range.

The acquiring subunit may be further configured to calculate a reference wind speed at the buoy height z using the following equation:

wherein, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Denotes the coefficient of friction, k denotes the von Karman constant,

represents the stable correction function, and l represents the Monin-Obukhov length.

The wind field acquisition device based on the doppler meteorological radar can further specifically include: the device comprises a judging module and an executing module.

Specifically, the determining module may be configured to determine a magnitude relationship between the target radial velocity and a preset velocity threshold; and the execution module can be used for executing the step of performing variation inversion processing on the radar data to obtain a target inversion wind field under low altitude when the target radial velocity is determined to be greater than or equal to the velocity threshold.

The obtaining module 11 may further specifically include: an acquisition unit and a composition unit.

Specifically, the acquiring unit may be configured to acquire a target radial velocity of each doppler meteorological radar; and the composite unit can be used for carrying out composite operation on the target radial speed by adopting a quadrilateral rule to obtain the target composite radial speed of the low-altitude radar network.

For specific limitations of the wind field acquisition device based on the doppler meteorological radar, reference may be made to the above limitations of the wind field acquisition method based on the doppler meteorological radar, and details are not repeated here. All or part of the modules in the wind field acquisition device based on the Doppler meteorological radar can be realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.

In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 9. The computer device includes a processor, a memory, a communication interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless communication can be realized through WIFI, an operator network, NFC (near field communication) or other technologies. The computer program is executed by a processor to implement a method for wind field acquisition based on a doppler meteorological radar. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.

Those skilled in the art will appreciate that the architecture shown in fig. 9 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.

In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

and performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of the low-altitude wind field, and taking the first target longitude and latitude wind field and the first target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

and performing four-dimensional variation inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of a low-altitude wind field; and performing four-dimensional variation inversion processing on radar data according to the first target longitude and latitude wind field and the first target height wind field to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field as the low-altitude target inversion wind field.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

substituting the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation and a smoothness constraint equation to construct a target function; minimizing the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of a wind field under low altitude; and performing inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target height wind field.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

judging whether the first target longitude and latitude wind field and the first target height wind field meet a preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment; and when the first target longitude and latitude wind field and the first target height wind field are determined to meet the time discontinuity condition, performing four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

acquiring a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and acquiring reference wind speeds of a buoy in a Doppler meteorological radar detection range at different heights and a reference average root mean square error or a reference average absolute error of the reference wind speeds; judging whether a first error difference value between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range or not; and determining the first average root mean square error or the average absolute error, and performing four-dimensional variation inversion processing on radar data when an error difference value between the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

constructing a value function according to a background field change constraint equation and a space-time smooth term constraint equation in the radar data; and carrying out minimization processing on the value function to obtain a second target longitude and latitude wind field and a second target height wind field.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

judging whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target height wind field and a second error difference value between the reference average root mean square error or the reference average absolute error of the reference wind speed are within a preset target difference value range or not; when the second error difference value is determined to be within the target difference value range, outputting first prompt information with higher consistency level between the target inversion wind field and the reference wind speed; and when the second error difference is determined not to be in the target difference range, outputting second prompt information that the target inversion wind field is invalid.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

the acquiring of the reference wind speeds of the buoy in the detection range of the Doppler meteorological radar at different heights comprises calculating the reference wind speed when the height of the buoy is z by adopting the following formula:

wherein, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Denotes the coefficient of friction, k denotes the von Karman constant,

represents the stable correction function, and l represents the Monin-Obukhov length.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

judging the magnitude relation between the target radial speed and a preset speed threshold; and when the target radial velocity is determined to be greater than or equal to the velocity threshold, executing the variational inversion processing on the radar data to obtain a target inversion wind field under low altitude.

In one embodiment, the processor, when executing the computer program, further performs the steps of:

the number of Doppler meteorological radar is at least two, and when constituting the radar network deployment based on at least two Doppler meteorological radars, acquire Doppler meteorological radar's target radial velocity under the low latitude, include: acquiring the target radial velocity of each Doppler meteorological radar; and performing composite operation on the target radial speed by adopting a quadrilateral rule to obtain the target composite radial speed of the low-altitude radar network.

It should be clear that, in the embodiments of the present application, the process of executing the computer program by the processor is consistent with the process of executing the steps in the above method, and specific reference may be made to the description above.

In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of:

acquiring the target radial velocity of the Doppler meteorological radar under the low altitude;

according to the target radial velocity, performing variation inversion processing on radar data to obtain a target inversion wind field under low altitude; the radar data are used for representing background fields, rainwater content, radial speed and noise of a Doppler meteorological radar detection area.

In one embodiment, the computer program when executed by the processor further performs the steps of:

and performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of the low-altitude wind field, and taking the first target longitude and latitude wind field and the first target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the computer program when executed by the processor further performs the steps of:

and performing four-dimensional variation inversion processing on the radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field of the low-altitude wind field as target inversion wind fields under low altitudes.

In one embodiment, the computer program when executed by the processor further performs the steps of:

performing three-dimensional variation inversion processing on the radar data to obtain a first target longitude and latitude wind field and a first target height wind field of a low-altitude wind field; and performing four-dimensional variation inversion processing on radar data according to the first target longitude and latitude wind field and the first target height wind field to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field, and taking the second target longitude and latitude wind field and the second target height wind field as the low-altitude target inversion wind field.

In one embodiment, the computer program when executed by the processor further performs the steps of:

substituting the radar data into a background field constraint equation, an observation fitting constraint equation, a rainwater content conservation equation, a continuous wind field constraint equation and a smoothness constraint equation to construct a target function; minimizing the target function to obtain a three-dimensional target longitude and latitude wind field and a three-dimensional target height wind field of a wind field under low altitude; and performing inversion processing on the three-dimensional target longitude and latitude wind field and the three-dimensional target height wind field in a rectangular coordinate system to obtain the first target longitude and latitude wind field and the first target height wind field.

In one embodiment, the computer program when executed by the processor further performs the steps of:

judging whether the first target longitude and latitude wind field and the first target height wind field meet a preset time discontinuity condition; the time discontinuity condition is used for representing radar data at the current moment and is used for acquiring a target longitude and latitude wind field and a target height wind field at the current moment; and when the first target longitude and latitude wind field and the first target height wind field are determined to meet the time discontinuity condition, performing four-dimensional variation inversion processing on radar data to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the computer program when executed by the processor further performs the steps of:

acquiring a first average root mean square error or a first average absolute error of the first target longitude and latitude wind field and the first target altitude wind field, and acquiring reference wind speeds of a buoy in a Doppler meteorological radar detection range at different heights and a reference average root mean square error or a reference average absolute error of the reference wind speeds; judging whether a first error difference value between the first average root mean square error or the first average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range or not; and determining the first average root mean square error or the average absolute error, and performing four-dimensional variation inversion processing on radar data when an error difference value between the first average root mean square error or the average absolute error and the reference average root mean square error or the reference average absolute error is within a preset difference value range to obtain a second target longitude and latitude wind field and a second target height wind field of the low-altitude wind field.

In one embodiment, the computer program when executed by the processor further performs the steps of:

constructing a value function according to a background field change constraint equation and a space-time smooth term constraint equation in the radar data; and carrying out minimization processing on the value function to obtain a second target longitude and latitude wind field and a second target height wind field.

In one embodiment, the computer program when executed by the processor further performs the steps of:

judging whether a second average root mean square error or a second average absolute error of the second target longitude and latitude wind field and the second target height wind field and a second error difference value between the reference average root mean square error or the reference average absolute error of the reference wind speed are within a preset target difference value range or not; when the second error difference value is determined to be within the target difference value range, outputting first prompt information with higher consistency level between the target inversion wind field and the reference wind speed; and when the second error difference is determined not to be in the target difference range, outputting second prompt information that the target inversion wind field is invalid.

In one embodiment, the computer program when executed by the processor further performs the steps of:

the acquiring of the reference wind speeds of the buoy in the detection range of the Doppler meteorological radar at different heights comprises calculating the reference wind speed when the height of the buoy is z by adopting the following formula:

wherein, V_zRepresents the reference wind speed u at the buoy height z in the detection range of the Doppler meteorological radar_*Denotes the coefficient of friction, k denotes the von Karman constant,

represents the stable correction function, and l represents the Monin-Obukhov length.

In one embodiment, the computer program when executed by the processor further performs the steps of:

judging the magnitude relation between the target radial speed and a preset speed threshold; and when the target radial velocity is determined to be greater than or equal to the velocity threshold, executing the variational inversion processing on the radar data to obtain a target inversion wind field under low altitude.

In one embodiment, the computer program when executed by the processor further performs the steps of:

the number of Doppler meteorological radar is at least two, and when constituting the radar network deployment based on at least two Doppler meteorological radars, acquire Doppler meteorological radar's target radial velocity under the low latitude, include: acquiring the target radial velocity of each Doppler meteorological radar; and performing composite operation on the target radial speed by adopting a quadrilateral rule to obtain the target composite radial speed of the low-altitude radar network.

It should be clear that, in the embodiments of the present application, the process of executing the computer program by the processor is consistent with the process of executing the steps in the above method, and specific reference may be made to the description above.

It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.

The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.

The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.