CN111487596A - Wind field detection data fusion method and device and electronic equipment - Google Patents

Wind field detection data fusion method and device and electronic equipment Download PDF

Info

Publication number
CN111487596A
CN111487596A CN202010314980.8A CN202010314980A CN111487596A CN 111487596 A CN111487596 A CN 111487596A CN 202010314980 A CN202010314980 A CN 202010314980A CN 111487596 A CN111487596 A CN 111487596A
Authority
CN
China
Prior art keywords
wind
detection
height
data
detection data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010314980.8A
Other languages
Chinese (zh)
Other versions
CN111487596B (en
Inventor
郑秀明
聂杨
王志锐
夏一凡
陈华彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aerospace New Weather Technology Co ltd
Original Assignee
Aerospace New Weather Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Aerospace New Weather Technology Co ltd filed Critical Aerospace New Weather Technology Co ltd
Priority to CN202010314980.8A priority Critical patent/CN111487596B/en
Publication of CN111487596A publication Critical patent/CN111487596A/en
Application granted granted Critical
Publication of CN111487596B publication Critical patent/CN111487596B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

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

Abstract

The invention relates to the technical field of data processing, in particular to a wind field detection data fusion method and device and electronic equipment. The method comprises the following steps: determining a first splicing height set of different detection modes of each wind profile radar; determining a second set of stitching heights between a plurality of wind profile radars having different detection ranges; acquiring detection data of a plurality of wind profile radars in the same wind field; selecting a first splicing height, and splicing and fusing detection data of different detection modes of each wind profile radar; and selecting a second splicing height, and splicing and fusing the detection data of the wind profile radar in different detection ranges. According to the method, the detection data of different detection modes of each wind profile radar and the detection data of the wind profile radars in different detection ranges are spliced and fused by determining the splicing heights of the different detection modes and the different wind profile radars of the wind profile radars, and compared with manual fusion, the method is good in fusion effect and the spliced data are not easy to mutate.

Description

Wind field detection data fusion method and device and electronic equipment
Technical Field
The invention relates to the technical field of data processing, in particular to a wind field detection data fusion method and device and electronic equipment.
Background
The wind profile radar is an atmospheric remote sensing device using turbulence as a tracer and is limited by a wind measuring principle, and the traditional wind profile radar cannot provide an accurate atmospheric wind field which is close to the ground and is about 100 meters below, so that a complete wind measuring system is usually provided with a wind speed and direction sensor on a ground station or a wind tower to complete wind measurement in cooperation with the wind profile. There is no clear method for effectively integrating the wind speed and direction sensor data into the wind profile radar data and calibrating the wind profile radar data.
In order to take into account the height resolution and the maximum detection power, a wind profile radar usually adopts multi-mode detection, the radar can obtain finer height resolution but smaller maximum detection power when detecting in a low mode, the radar can obtain larger maximum detection power but coarser height resolution when detecting in a high mode, and some wind profile radars are also designed with medium-mode detection, and the height resolution and the maximum detection power of the wind profile radars are between the low-mode detection and the high-mode detection. The detection data of different modes are crossed in height, mode splicing is mainly performed through subjective manual setting at present, and no specific optimal splicing method exists, so that sudden change of the data near the splicing height can occur.
In order to detect the atmospheric wind fields with different heights, two or more than two wind profile radars with different detection ranges are often adopted for networking detection, for example, a boundary layer wind profile radar with detection power of 3-5 kilometers and a troposphere wind profile radar with detection power of 8-12 kilometers are jointly used for detection. Data of different radars are crossed in height, and for how to perform data fusion, subjective manual setting is mainly used at present, so that the data fusion effect is poor, and data near the splicing height does not have a smooth transition process and is easy to mutate.
Disclosure of Invention
In view of this, embodiments of the present invention provide a method and an apparatus for fusing wind field detection data, and an electronic device, so as to solve the problem in the prior art that the detected wind field data is mainly fused by manual setting, and the fusion effect is poor.
According to a first aspect, an embodiment of the present invention provides a wind field detection data fusion method, including:
determining a first splicing height set of different detection modes of each wind profile radar; wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar;
determining a second set of stitching heights between a plurality of wind profile radars having different detection ranges; the second splicing height set comprises splicing heights for splicing detection data of a plurality of wind profile radars in different detection ranges;
acquiring detection data of the plurality of wind profile radars in the same wind field;
selecting a corresponding first splicing height from the first splicing height set, splicing and fusing detection data of different detection modes of each wind profile radar, wherein the detection data after splicing and fusing take the detection data at the corresponding first splicing height as boundary points;
and selecting a corresponding second splicing height from the second splicing set, splicing and fusing the detection data of the wind profile radars in different detection ranges, and obtaining the wind field detection data of the wind profile radars from the minimum detection height to the maximum detection height by using the detection data at the corresponding second splicing height as boundary points of the spliced and fused detection data.
The wind field detection data fusion method provided by the embodiment of the invention determines the first splicing height set of each wind profile radar in different detection modes, and a second splicing height set among a plurality of wind profile radars with different detection ranges, and a first splicing height for splicing and fusing detection data of different detection modes of each wind profile radar and a second splicing height for splicing and fusing detection data of the plurality of wind profile radars with different detection ranges are selected from the first splicing height set and the second splicing height set, the wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height are obtained, the wind profile radars with different detection modes and different detection ranges are objectively fused with the detection data of the same wind field, and randomness and uncertainty of artificially and subjectively setting the splicing height are avoided.
With reference to the first aspect, in a first implementation manner of the first aspect, the method further includes:
acquiring detection data of a wind speed and wind direction sensor;
preprocessing the detection data of the wind speed and wind direction sensor;
and performing linear interpolation on the preprocessed detection data by using the detection data of the wind profile radar with the lowest detection height, and performing fusion calibration on the detection data of the wind profile radar with the lowest detection height to obtain the detection data obtained by fusing the detection data of the wind speed and direction sensor and the detection data of the wind profile radar with the lowest detection height.
According to the wind field detection data fusion method provided by the embodiment of the invention, the detection data of the wind speed and wind direction sensor and the detection data of the wind profile radar with the lowest detection height are fused to obtain detection data. When the same wind field is detected, because the data of the ground cannot be detected due to the limitation of the lowest detection height of the wind profile radar, a wind speed and direction sensor is generally combined with the wind profile radar to obtain the overall detection data of the wind field from the ground to the highest detection height of the wind profile radar, and the defect of low-altitude detection performance of the wind profile radar is overcome.
With reference to the first aspect, in a second implementation manner of the first aspect, the determining a first splice height set of different detection modes of each wind profile radar includes:
acquiring a data acquisition rate of each height layer of each detection mode of each wind profile radar, wherein each wind profile radar has n detection modes, and n is greater than 1;
selecting the data acquisition within a detection range shared by the ith detection mode and the (i + 1) th detection modeRate above threshold η1Height set A ofi(ii) a Taking i from 1 to n-1;
calculating the height set AiThe wind measurement error is the square sum of the average error and the standard deviation of the detection data of the ith detection mode and the detection data of the (i + 1) th detection mode;
selecting the maximum height h in the M height layers with the minimum wind measurement errori；
Adding 1 to the value of i, and returning to execute that the data acquisition rate is higher than the threshold η in the detection range shared by the ith detection mode and the (i + 1) th detection mode1Height set A ofiUntil h is obtained1,…,hn-1Forming the first set of splice heights.
According to the wind field data fusion method provided by the embodiment of the invention, the splicing height determined through error analysis is an objective evaluation of the wind measuring performance of each wind profile radar in different detection modes, and compared with the traditional artificial subjective setting, the optimal splicing height of the wind profile radar in the same wind field and different detection modes can be reflected, so that the sudden change of the data of the same wind profile radar in different detection modes after fusion is avoided; furthermore, the maximum height h in the M height layers with the minimum wind measurement error is selectediThe detection mode with higher height resolution can provide effective detection data in a wider detection height range.
With reference to the first aspect, in a third implementation manner of the first aspect, the determining a second set of splice heights between a plurality of wind profile radars with different detection ranges includes:
arranging the plurality of wind profile radars according to the sequence of the maximum detection height from small to large, and selecting half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]J is 1 to m-1, and m is the total number of the wind profile radar;
according to historical data, the fact that the data acquisition rate in the detection height range shared by the jth wind profile radar and the jth +1 wind profile radar is higher than a threshold η is counted2Height set B ofj；
Calculating the height set BjWind measurement error of each height layer; the wind measurement error is the square sum of the average error and the standard deviation of the horizontal wind speed detected by the jth wind profile radar and the jth +1 wind profile radar and the standard deviation of the horizontal wind direction;
selecting the maximum height H in the N height layers with the minimum wind measurement errorj；
Adding 1 to the value of j, and returning to execute the selection of half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]Until H is obtained1,…,Hm-1And forming the second splicing height set.
According to the wind field data fusion method provided by the embodiment of the invention, the splicing height selected by the wind measurement error is utilized, the horizontal wind speed error and the horizontal wind direction error are optimized, the high spatial resolution of the wind profile radar with the minimum detection height is fully utilized, the square sum of the standard deviation is selected as the error of the horizontal wind direction, the influence of the difference of the selection of azimuth angles on the error of the horizontal wind direction of different wind profile radars is avoided, and the data of different wind profile radars are prevented from being mutated after being fused.
With reference to the first embodiment of the first aspect, in a fourth embodiment of the first aspect, the preprocessing the detection data of the wind speed and wind direction sensor includes:
carrying out primary smoothing on the detection data of the wind speed and direction sensor according to a vector averaging method to obtain detection data after primary smoothing;
and performing secondary smoothing on the detection data subjected to the primary smoothing by taking the time interval generated by the detection data of the wind profile radar as a smoothing window according to a vector averaging method to obtain the preprocessed detection data.
According to the wind field data fusion method provided by the embodiment of the invention, the detection data of the wind speed and wind direction sensor is preprocessed, so that the preprocessed detection data is relatively stable and representative, the time scale of the preprocessed detection data is close to that of a wind profile radar, so that data fusion is convenient, and meanwhile, the detection data of the wind speed and wind direction sensor is smoothly processed by adopting a vector average method instead of a scalar average method, so that the vector characteristic of a wind field is more objectively embodied.
With reference to the first embodiment of the first aspect, in a fifth embodiment of the first aspect, the performing linear interpolation on the preprocessed detection data by using the detection data of the wind profile radar with the lowest detection height, and performing fusion calibration on the detection data of the wind profile radar with the lowest detection height includes:
decomposing the preprocessed detection data into u with orthogonal directionsd、vdTwo sets of horizontal wind velocity components;
decomposing the detection data of the wind profile radar with the lowest detection height into U, V two groups of horizontal wind speed components with orthogonal directions;
and u is obtained by decomposing the data of the minimum height layer of which the data acquisition rate of the wind profile radar with the lowest detection height is higher than a preset threshold valuec、vcTwo sets of horizontal wind velocity components and ud、vdTaking the two groups of horizontal wind speed components as a reference, carrying out linear difference according to the layering height of the wind profile radar with the lowest detection height to obtain two groups of data of U 'and V', and synthesizing the two groups of data of U 'and V' to obtain wind speed and wind direction data of each height layer after linear interpolation;
comparing the wind speed and wind direction data of each height layer with the original wind speed and wind direction data detected by the wind profile radar with the lowest detection height;
and if the absolute value of the difference of the wind speeds or the absolute value of the difference of the wind directions is larger than a preset threshold value, replacing the original wind speed and wind direction data of the height layer of the wind profile radar with the lowest detection height with the wind speed and wind direction data of the height layer after linear interpolation.
With reference to the fifth embodiment of the first aspect, in the sixth embodiment of the first aspect, the preset threshold is calculated by the following formula:
Figure BDA0002458086220000051
wherein Th is a reference value of the preset threshold value, ThiA predetermined threshold value, h, representing the ith hierarchy leveliDenotes the ith hierarchy height, hdownIndicating the height, h, of the anemometric sensorupRadar u representing said wind profilec、vcThe corresponding stratification height, η, represents the adjustment factor.
According to a second aspect, an embodiment of the present invention provides a wind field detection data fusion apparatus, including:
the first determining module is used for determining a first splicing height set of different detection modes of each wind profile radar; wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar;
a second determination module to determine a second set of stitching heights between a plurality of wind profile radars having different detection ranges; the second splicing height comprises a splicing height for splicing detection data of a plurality of wind profile radars in different detection ranges;
the acquisition module is used for acquiring detection data of the plurality of wind profile radars in the same wind field;
the first fusion module is used for selecting a corresponding first splicing height from the first splicing height set, splicing and fusing the detection data of each wind profile radar in different detection modes, and the spliced and fused detection data takes the detection data at the corresponding first splicing height as a boundary point;
and the second fusion module is used for selecting a corresponding second splicing height from the second splicing set, splicing and fusing the detection data of the wind profile radars in different detection ranges, and obtaining the wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height by using the detection data at the corresponding second splicing height as boundary points in the spliced and fused detection data.
The wind field detection data fusion device provided by the embodiment of the invention determines the first splicing height set of each wind profile radar in different detection modes, and a second splicing height set among a plurality of wind profile radars with different detection ranges, and a first splicing height for splicing and fusing detection data of different detection modes of each wind profile radar and a second splicing height for splicing and fusing detection data of the plurality of wind profile radars with different detection ranges are selected from the first splicing height set and the second splicing height set, the wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height are obtained, the wind profile radars with different detection modes and different detection ranges are objectively fused with the detection data of the same wind field, and randomness and uncertainty of artificially and subjectively setting the splicing height are avoided.
According to a third aspect, an embodiment of the present invention further provides an electronic device, including:
the wind field detection data fusion method comprises a memory and a processor, wherein the memory and the processor are connected with each other in a communication mode, the memory stores computer instructions, and the processor executes the computer instructions so as to execute the wind field detection data fusion method according to the first aspect of the invention or any embodiment of the first aspect of the invention.
According to a fourth aspect, the embodiment of the present invention provides a computer-readable storage medium, which stores computer instructions for causing a computer to execute the wind farm detection data fusion method according to the first aspect of the present invention or any implementation manner of the first aspect.
The invention achieves the following technical effects by adopting the technical method:
1. by determining the splicing heights of different detection modes and different wind profile radars of the wind profile radars, the detection data of the different detection modes of each wind profile radar and the detection data of the wind profile radars in different detection ranges are spliced and fused, and compared with manual fusion, the fusion effect is good and the splicing data are not easy to mutate.
2. The splicing height is selected from the minimum error by calculating the error of the detection data of the common detection range of the radars with different wind profiles and the error of the detection data of the common detection range of different detection modes of the same wind profile, so that the sudden change after data fusion is avoided.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method of wind farm detection data fusion according to an embodiment of the present invention;
FIG. 2 is a flow chart of a data fusion method of a wind speed and direction sensor and a wind profile radar according to an embodiment of the invention;
FIG. 3 is a complete flow chart of a method for data fusion of a anemometry sensor and a wind profile radar according to an embodiment of the present invention;
FIG. 4 is a comparison of before and after pretreatment;
FIG. 5 is a flowchart of the steps for obtaining a first set of splice heights in accordance with an embodiment of the present invention;
FIG. 6 is a flowchart of the steps for obtaining a second set of splice heights in accordance with an embodiment of the present invention;
FIG. 7 is a schematic diagram of an application scenario of wind farm full data fusion according to an embodiment of the present invention;
FIG. 8 is a block diagram of a wind farm detection data fusion apparatus according to an embodiment of the present invention;
fig. 9 is a schematic diagram of a hardware structure of an electronic device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
There is provided a wind farm detection data fusion method according to an embodiment of the present invention, it should be noted that the steps illustrated in the flowchart of the drawings may be performed in a computer system such as a set of computer executable instructions, and that while a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different than here.
The method of the embodiment of the invention is described by taking a ground station wind speed and direction sensor, a low-altitude wind profile radar with the maximum detection height of 600m, a boundary layer wind profile radar with the maximum detection height of 3-5 Km and a troposphere wind profile radar with the maximum detection height of 8-10 Km as examples, and the specific embodiments are only used for explanation and are not used for limiting the scope of the invention.
In the embodiment of the invention, the speed of acquiring the detection data by the wind speed and wind direction sensor is in millisecond level; the different detection modes, detection ranges and height resolutions of the low-altitude wind profile radar, the boundary layer wind profile radar and the troposphere wind profile radar are shown in table 1:
TABLE 1 detection Performance of three different types of wind profile radar in different detection modes
Figure BDA0002458086220000081
Example 1
In this embodiment, a wind field data fusion method is provided, which can be used in the above electronic device, and fig. 1 is a flowchart of a wind field detection data fusion method according to an embodiment of the present invention, as shown in fig. 1, the flowchart includes the following steps:
s11, a first set of splice heights for different detection modes of each wind profile radar is determined.
Wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar.
Table 1 shows the detection performance of the different detection modes of the three different types of wind profile radars, each wind profile radar has a coincidence of detection ranges in the different detection modes, and when detecting the same wind field, the detection data corresponding to the coincident detection ranges also coincide, which is not beneficial to embodying the overall detection information of the wind field, so that the three wind profile radars are spliced and fused by using the height as a reference.
Firstly, the detection data of different detection modes of each of the above-mentioned wind profile radars are spliced and fused, specifically, taking the tropospheric wind profile radars in table 1 as an example, the coincidence detection range of the low mode and the medium mode is 2Km to 4Km, and the coincidence detection range of the medium mode and the high mode is 4Km to 7Km, so that when the tropospheric wind profile radar detects the same wind field, the detection data within 2Km to 4Km and 4Km to 7Km coincide, and further, the corresponding splice heights, for example, 3Km and 6Km, are respectively selected from the two coincidence detection ranges to constitute a first splice height set of different detection modes of the tropospheric wind profile radar.
It should be noted that the type and detection mode of the wind profile radar may also have other types and detection modes different from those shown in table 1, and table 1 only serves to explain the embodiment of the present invention.
S12, a second set of stitching heights between a plurality of wind profile radars having different detection ranges is determined.
Wherein the second set of stitching heights comprises stitching heights for stitching detection data of a plurality of wind profile radars of different detection ranges.
In combination with the above, the detection height ranges of the three wind profile radars shown in table 1 are 20m to 600m, 60m to 5Km and 120m to 10Km, respectively, and the coincidence detection ranges of the low-altitude wind profile radar and the boundary layer wind profile radar are as follows: 60 m-600 m, the coincidence detection range of the boundary layer wind profile radar and the troposphere wind profile radar is as follows: 120m to 5Km, and corresponding splicing heights, for example, 400m and 3Km, can be selected from two coincident detection ranges respectively by using the pareto optimal principle to form a second splicing height set among the three wind profile radars.
And S13, acquiring the detection data of the plurality of wind profile radars in the same wind field.
The detection data of the plurality of wind profile radars in the same wind field acquired by the electronic device can be data detected in real time from the wind profile radar installed in the wind field; or may be detection data of the wind field stored in the electronic device; or the electronic equipment acquires the detection data of the wind field from the outside in other modes. No matter what way the electronic equipment acquires the detection data of the wind field, the electronic equipment is only required to be ensured to acquire the detection data of the wind field.
And S14, splicing and fusing the detection data of different detection modes of each wind profile radar.
Specifically, a corresponding first splicing height is selected from the first splicing height set, splicing and fusing are performed on detection data of different detection modes of each wind profile radar, and the detection data after splicing and fusing take the detection data at the corresponding first splicing height as boundary points.
And combining S11, when carrying out data splicing on the troposphere wind profile radar, selecting low-mode detection data within the range of 120 m-3 Km, selecting medium-mode detection data within the range of 3 Km-6 Km, selecting high-mode detection data within the range of 6 Km-10 Km, and then splicing according to the height sequence.
And S15, splicing and fusing the detection data of the wind profile radar in different detection ranges.
Specifically, the corresponding second splicing heights are selected from the second splicing set, splicing and fusing the detection data of the wind profile radars in different detection ranges, and the spliced and fused detection data take the detection data at the corresponding second splicing heights as boundary points to obtain the wind field detection data of the wind profile radars from the minimum detection height to the maximum detection height.
And combining S12, when data splicing is carried out, selecting the data of the low-altitude wind profile radar as the detection data within the range of 20-400 m, selecting the data of the boundary layer wind profile radar as the detection data within the range of 400-3 Km, selecting the data of the troposphere wind profile radar as the detection data within the range of 3 Km-10 Km, and splicing according to the height sequence to obtain the wind field detection data from the minimum detection height to the maximum detection height of the three wind profile radars.
According to the wind field detection data fusion method provided by the embodiment of the invention, the wind field detection data from the minimum detection height to the maximum detection height of the plurality of wind profile radars are obtained by determining the first splicing height set of each wind profile radar in different detection modes and the second splicing height set of the plurality of wind profile radars with different detection ranges, selecting the first splicing height for splicing and fusing the detection data of each wind profile radar in different detection modes and the second splicing height for splicing and fusing the detection data of the plurality of wind profile radars with different detection ranges from the first splicing height set and the second splicing height set, and fusing the detection data of the plurality of wind profile radars and the detection data of the wind profile radars in different modes. By determining the splicing heights of different detection modes and different wind profile radars of the wind profile radars, the detection data of the different detection modes of each wind profile radar and the detection data of the wind profile radars in different detection ranges are spliced and fused, and compared with manual fusion, the fusion effect is good and the splicing data are not easy to mutate.
Optionally, fig. 2 is a flowchart of a data fusion method of a wind speed and direction sensor and a wind profile radar according to an embodiment of the present invention, and as shown in fig. 2, the method includes the following steps:
and S21, acquiring the detection data of the wind speed and wind direction sensor.
The detection data of the wind speed and wind direction sensor acquired by the electronic equipment can be data detected in real time by the wind speed and wind direction sensor installed in the wind field; or the detection data of a wind speed and wind direction sensor stored in the electronic equipment; or the detection data of the wind speed and wind direction sensor acquired by the electronic equipment from the outside in other ways. No matter what way the electronic equipment acquires the detection data of the wind speed and the wind direction, the electronic equipment is only required to be ensured to acquire the detection data of the wind speed and the wind direction sensor.
And S22, preprocessing the detection data of the wind speed and wind direction sensor.
Because the wind profile radar generally generates a group of detection data within 3-10 minutes, the wind speed and direction sensor generates 1-10 groups of detection data per second, and the time scales of the detection data of the wind speed and direction sensor are different, the detection data of the wind speed and direction sensor needs to be preprocessed, so that the time scales of the detection data of the wind profile radar and the detection data of the wind profile radar are consistent. In addition, the preprocessed wind speed and wind direction sensor detection data are relatively stable and representative.
As an alternative implementation manner of this embodiment, as shown in fig. 3, the S22 may include:
and S221, performing primary smoothing on the detection data of the wind speed and direction sensor according to a vector average method to obtain the detection data after the primary smoothing.
Specifically, the wind speed and direction sensor generally collects 1-10 groups of detection data every second, original wind speed and direction data with standard deviation larger than a preset value collected by the wind speed and direction sensor are decomposed to obtain two orthogonal components, the two components are respectively averaged according to a second-level time window and then synthesized to obtain an average value of the wind speed and the wind direction, and detection data with small jitter relative to the original data are formed.
S222, performing a second smoothing process according to a vector averaging method.
And performing secondary smoothing on the detection data subjected to the primary smoothing by using a time interval generated by the detection data of the wind profile radar as a smoothing window according to a vector averaging method to obtain the preprocessed detection data, wherein a comparison graph before and after preprocessing is shown in fig. 4.
And S23, performing linear interpolation on the preprocessed detection data by using the detection data of the wind profile radar with the lowest detection height, and performing fusion calibration on the detection data of the wind profile radar with the lowest detection height.
With reference to table 1, the lowest detection height of the wind profile radar is 20m (low-altitude wind profile radar), and the wind speed and direction sensor can be used for detecting within the height range of 0-20 m, so that when the detection data of the whole wind field are fused, only linear interpolation and fusion calibration are needed to be performed on the detection data preprocessed by the wind speed and direction sensor and the wind profile radar (low-altitude wind profile radar) with the lowest detection height, and after the detection data fused by the detection data of the wind speed and direction sensor and the detection data of the wind profile radar with the lowest detection height are obtained, the whole detection data of the wind field can be obtained.
According to the wind field detection data fusion method provided by the embodiment of the invention, the detection data of the wind speed and wind direction sensor and the detection data of the wind profile radar with the lowest detection height are fused to obtain detection data. When the same wind field is detected, because the data of the ground cannot be detected due to the limitation of the lowest detection height of the wind profile radar, a wind speed and direction sensor is generally combined with the wind profile radar to obtain the overall detection data of the wind field from the ground to the highest detection height of the wind profile radar, and the defect of low-altitude detection performance of the wind profile radar is overcome.
As an optional implementation manner of the embodiment of the present invention, the S23 may include:
(1) decomposing the preprocessed detection data into u with orthogonal directionsd、vdTwo sets of horizontal wind velocity components.
Specifically, the height of the preprocessed detection data is set as hdown。
(2) And decomposing the detection data of the wind profile radar with the lowest detection height into U, V two groups of horizontal wind speed components with orthogonal directions.
(3) And u obtained by decomposing the data of the minimum height layer with the data acquisition rate of the wind profile radar with the lowest detection height being higher than a preset threshold valuec、vcTwo sets of horizontal wind velocity components and ud、vdAnd taking the two groups of horizontal wind speed components as a reference, carrying out linear difference according to the layering height of the wind profile radar with the lowest detection height to obtain two groups of data of U 'and V', and synthesizing the two groups of data of U 'and V' to obtain wind speed and wind direction data of each height layer after linear interpolation.
In particular, said preset threshold is 90%, for a height h in U', ViHorizontal wind velocity component u ofi、viThe method is obtained by interpolation according to the following formula:
Figure BDA0002458086220000121
Figure BDA0002458086220000131
(4) and comparing the wind speed and wind direction data of each height layer with the original wind speed and wind direction data detected by the wind profile radar with the lowest detection height.
(5) And if the absolute value of the difference of the wind speeds or the absolute value of the difference of the wind directions is larger than a preset threshold value, replacing the original wind speed and wind direction data of the height layer of the wind profile radar with the lowest detection height with the wind speed and wind direction data of the height layer after linear interpolation.
Otherwise, if the absolute value of the difference of the wind speeds or the difference of the wind directions is smaller than or equal to the preset threshold value, maintaining the original detection data of the wind profile radar with the lowest detection height of the height layer.
As an optional implementation manner of the embodiment of the present invention, the preset threshold may be calculated by the following formula:
Figure BDA0002458086220000132
wherein Th is a reference value of the preset threshold value, ThiA predetermined threshold value, h, representing the ith hierarchy leveliDenotes the ith hierarchy height, hdownIndicating the height, h, of the anemometric sensorupRadar u representing said wind profilec、vcThe corresponding stratification height, η, represents the adjustment factor.
In particular, the predetermined threshold value may be optimized as a variable varying with height, considering that the strictness of the quality control should vary with the height level. For the height layer to be checked, if it is at the height h of the ground station (the height of the wind speed and direction sensor)downOr height hupIf yes, the check threshold is set Th; the farther away from hdownOr height hupThe higher the two height levels, i.e. the less stringent the quality control, the higher the verification threshold is at hdownOr height hupThe two height layers reach the maximum and are amplified by η times, the quality control strength of the detection data can be set by controlling the size of η, for example, when echoes are weak, clutter is strong, and identification is difficult, η can be reduced, quality control is enhanced, otherwise, in the weather change process or when echoes are strong, η can be increased, quality control is weakened, and the influence of a reference estimation value on the actual wind detection effect is reduced.
As an alternative to this embodiment, as shown in fig. 5, the first splicing height set is obtained by the following steps:
and S31, acquiring the data acquisition rate of each height layer of each detection mode of each wind profile radar. Each wind profile radar has n detection modes, n being greater than 1.
The data acquisition rate of each height layer of each detection mode of each wind profile radar acquired by the electronic device may be a data acquisition rate sent by the wind profile radar during or after detection; or may be a data acquisition rate stored in the electronic device; or the data acquisition rate acquired by the electronic device from the outside in other ways. No matter how the electronic equipment obtains the data obtaining rate, the electronic equipment only needs to be ensured to obtain the data obtaining rate.
The data acquisition rate refers to the percentage of the number of times of quality control of detection data to the total detection number of times in detection of the wind profile radar in a period of time, and since the echo signal of the wind profile radar is generally gradually weakened along with the rise of the height, the data acquisition rate is generally counted respectively according to different height layers.
Taking the tropospheric profile radar in table 1 as an example, n is 3, that is, the tropospheric profile radar has three detection modes, each detection mode has a different detection altitude range, and the electronic device acquires the data acquisition rate of each altitude layer in the detection altitude range.
S32, selecting the data acquisition rate in the shared detection range to be higher than the threshold η1Height set A ofi。
Specifically, the data acquisition rate is selected to be higher than the threshold η in the detection range shared by the ith detection mode and the (i + 1) th detection mode1Height set A ofi(ii) a The i is 1 to n-1.
Taking the tropospheric wind profile radar in table 1 as an example, where n is 3, the data acquisition rate in the common detection range of the 1 st detection mode (low mode) and the 2 nd detection mode (medium mode) is selected to be higher than η1The height layers of (A) form the height set AiSimilarly, the data acquisition rate in the 2 nd detection mode and the 3 rd detection mode (high mode) shared detection range is higher than η1The height layers of (A) form the height set AiSaid threshold η1And taking 90 percent.
S33, calculating the height set AiThe wind measurement error of each height level.
The wind measurement error is the sum of the squares of the mean error and the standard deviation of the detection data of the ith detection mode and the detection data of the (i + 1) th detection mode.
Taking tropospheric wind profile radar in table 1 as an example, the common detection range of the 1 st detection mode and the 2 nd detection mode is 2Km to 4Km, and the wind measurement error of each height level in 2Km to 4Km of the 1 st detection mode and the 2 nd detection mode is calculated, wherein the wind measurement error is the sum of the mean error of the detection data and the square of the standard deviation, and the specific calculation formula is as follows:
average error:
Figure BDA0002458086220000151
standard deviation:
Figure BDA0002458086220000152
wherein x is the detection data of different detection modes in the same height layer, a is the sample number of the detection data, and xkMeans that the kth detection data in the sample. In the calculation, the deviation average error 3 sigma is used as a mark of a coarse error, and the coarse error is eliminated.
S34, selecting the maximum height h in the M height layers with the minimum wind measurement errori。
Comparing the wind measurement errors to obtain M height layers with the minimum wind measurement errors, and then comparing the M height layers to obtain the maximum height hi。
S35, repeating the above S32-S34 until h is obtained1,…,hn-1Forming the first set of splice heights.
Adding 1 to the value of i, and returning to execute that the data acquisition rate is higher than the threshold η in the detection range shared by the ith detection mode and the (i + 1) th detection mode1Height set A ofiUntil h is obtained1,…,hn-1Forming the first set of splice heights.
The wind field data fusion method provided by the embodiment of the invention determines the fusion result through error analysisThe splicing height is an objective evaluation of the wind measuring performance of each wind profile radar in different detection modes, compared with the traditional artificial subjective setting, the optimal splicing height of the wind profile radar in the same wind field and different detection modes can be reflected, and in addition, the maximum height h in the M height layers with the minimum wind measuring error is selectediThe detection mode with higher height resolution can provide effective detection data in a wider detection height range.
As an alternative to this embodiment, as shown in fig. 6, the second splicing height set is obtained by the following steps:
and S41, determining the splicing height range among different wind profile radars.
Arranging the plurality of wind profile radars according to the sequence of the maximum detection height from small to large, and selecting half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]J is 1 to m-1, and m is the total number of the wind profile radar.
For 3 wind profile radars participating in data fusion, arranging the maximum detection heights from small to large in sequence, selecting half of the maximum detection heights of two adjacent wind profile radars as the upper limit and the lower limit of the splicing height range, namely, the splicing height range of the low-altitude wind profile radar and the boundary layer wind profile radar shown in table 1 is [300m,2.5km ], the splicing height range of the boundary layer wind profile radar and the convection layer wind profile radar is [2.5km,5km ], and when the data fusion is performed on the splicing height ranges, inputting the height splicing ranges into electronic equipment in advance to limit the splicing height limit value of the detection data between the adjacent wind profile radars and avoid the detection data exceeding the splicing height limit value from being fused by mistake.
S42, counting that the data acquisition rate in the common detection range is higher than a threshold η2Height set B ofj。
Specifically, the statistics are counted according to historical dataThe data acquisition rate in the detection high range shared by the jth wind profile radar and the jth +1 wind profile radar is higher than a threshold η2Height set B ofj。
Taking 3 wind profile radars shown in table 1 as an example, a data acquisition rate threshold η is set2Taking 80% of the height as a credibility constraint condition for the fusion of the detection data of the 3 wind profile radars, counting the height of a first wind profile radar (low-altitude wind profile radar) and a second wind profile radar (boundary layer wind profile radar) in each height layer within a common detection height range of 60 m-600 m according to historical data, wherein the data acquisition rate is higher than 80%, and forming a height set B1(ii) a Counting the height of the second wind profile radar and the third wind profile radar (troposphere wind profile radar) within the common detection height range of 4-5 Km, wherein the data acquisition rate is higher than 80%, and forming the height set B2。
S43, calculating the height set BjThe wind measurement error of each height level.
And the wind measurement error is the sum of the square of the average error and the standard deviation of the horizontal wind speed detected by the jth wind profile radar and the jth +1 wind profile radar and the standard deviation of the horizontal wind direction.
Specifically, the above B is calculated separately1The average error of the horizontal wind speed of the first wind profile radar and the second wind profile radar and the square sum of the standard deviation are used as wind speed errors, and B is calculated2And taking the average error and the square sum of the standard deviation of the horizontal wind speeds of the second wind profile radar and the third wind profile radar as the wind measurement error.
And, calculating the above B separately1Calculating the standard deviation of the horizontal wind directions of the first wind profile radar and the second wind profile radar as wind direction errors, and respectively calculating the B2And the standard deviation of the horizontal wind directions of the second wind profile radar and the third wind profile radar is used as a wind measurement error.
S44, selecting the maximum height H in the N height layers with the minimum wind measurement errorj。
From the wind speed errorAngle of difference, for said B1、B2Each of the height layers is graded, the height layer with the wind speed error arranged in the front 30% is judged as the preferred height layer, and the height layer arranged in the back 70% is judged as the normal height layer; and, from the angle of the wind direction error, for B1、B2The height levels in (b) are classified into classes, and the height level with the wind direction error arranged at the front 20% is judged as the preferred height level, and the height level arranged at the rear 80% is judged as the normal height level. Sorting N height layers which are judged to be optimal height layers simultaneously according to the wind speed error and the wind direction error, and selecting the maximum height Hj。
If the heights of the layers which are judged as the preferred height layers do not appear at the same time, selecting the height with the largest wind speed error as the height of the preferred height as Hj。
S45, repeating the steps S42-S44 until H is obtained1,…,Hm-1And forming the second splicing height set.
Adding 1 to the value of j, and returning to execute the selection of half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]Until H is obtained1,…,Hm-1And forming the second splicing height set.
According to the wind field data fusion method provided by the embodiment of the invention, the splicing height selected by the wind measurement error is utilized, the horizontal wind speed error and the horizontal wind direction error are optimized, the high spatial resolution of the wind profile radar with the minimum detection height is fully utilized, the square sum of the standard deviation is selected as the error of the horizontal wind direction, and the influence of the difference of the selection of azimuth angles on the error of the horizontal wind direction of different wind profile radars is avoided.
An application scenario of wind field full-data fusion according to an embodiment of the present invention is described below with reference to fig. 7, and as shown in fig. 7, the wind field detection data fusion method according to the embodiment of the present invention includes:
s51, determining the splicing height of each wind profile radar in different modes through error analysis;
selecting a height layer of which the data acquisition rate is higher than a threshold in a coincidence detection range of a single radar in an adjacent detection mode, calculating the wind measurement error of historical data corresponding to the height layer, and selecting the height corresponding to the historical data with the minimum wind measurement error as the splicing height of the single radar in different modes.
S52, determining data splicing heights among the wind profile radars of the different wind profile radars based on a pareto optimal principle;
and selecting a height layer with a data acquisition rate higher than a threshold in the coincidence detection range of the different wind profile radars, calculating the wind measurement error of historical data corresponding to the height layer, and selecting the height layer corresponding to the historical data with the minimum wind measurement error as the data splicing height of the different wind profile radars by utilizing the pareto optimal principle.
S53, acquiring and updating detection data of the wind speed and direction sensor and the wind profile radar;
because the data generation time of the wind speed and direction sensor is inconsistent with the data generation time of the wind profile radar, the data fusion can not be directly carried out, so that the data of the wind speed and direction sensor is preprocessed, the generation frequency of the data of the wind speed and direction sensor is kept consistent with the generation frequency of the data of the wind profile radar, and a foundation is provided for the subsequent data fusion.
S54, combining S51, completing data splicing of different modes of each wind profile radar;
and splicing the data of different modes of each wind profile radar according to the splicing height obtained in the step S51.
S55, carrying out fusion calibration on the detection data of the wind profile radar by using the detection data of the wind speed and direction sensor;
and performing linear interpolation on the detection data of the wind speed and wind direction sensor and the detection data of the wind profile radar with the lowest detection height to obtain fused detection data, comparing the fused detection data of each height layer with the original detection data of each height layer of the wind profile radar, and replacing the original detection data of the wind profile radar of the height layer with the fused detection data of the height layer if the error is greater than a preset threshold value to finish calibration.
S56, combining S52, completing data splicing of radars with different wind profiles;
and splicing the detection data of the different wind profile radars according to the splicing height obtained in the step S52.
S57, waiting for a new round of data generation, repeating the above S53-S56.
The invention achieves the following technical effects by adopting the technical method:
1. by determining the splicing heights of different detection modes and different wind profile radars of the wind profile radars, the detection data of the different detection modes of each wind profile radar and the detection data of the wind profile radars in different detection ranges are spliced and fused, and compared with manual fusion, the fusion effect is good and the splicing data are not easy to mutate.
2. The splicing height is selected from the minimum error by calculating the error of the detection data of the common detection range of the radars with different wind profiles and the error of the detection data of the common detection range of different detection modes of the same wind profile, so that the sudden change after data fusion is avoided.
Example 2
The present embodiment provides a wind field detection data fusion device, which is used for implementing the above embodiments and implementation manners, and the description of the device is omitted. As used below, the term "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
The present embodiment provides a wind field detection data fusion apparatus, as shown in fig. 8, including:
a first determining module 61, configured to determine a first stitching height set of different detection modes of each wind profile radar; wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar;
a second determination module 62 for determining a second set of splice heights between a plurality of wind profile radars having different detection ranges; the second splicing height comprises a splicing height for splicing detection data of a plurality of wind profile radars in different detection ranges;
an obtaining module 63, configured to obtain detection data of the multiple wind profile radars in the same wind field;
a first fusion module 64, configured to select a corresponding first splicing height from the first splicing height set, splice and fuse detection data of different detection modes of each wind profile radar, where the detection data at the corresponding first splicing height is used as a boundary point in the spliced and fused detection data;
and a second fusion module 65, configured to select a corresponding second splicing height from the second splicing set, splice and fuse the detection data of the wind profile radars in different detection ranges, and obtain wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height by using the detection data at the corresponding second splicing height as a boundary point in the spliced and fused detection data.
Further functional descriptions of the modules are the same as those of the corresponding embodiments, and are not repeated herein.
The wind field detection data fusion device provided by the embodiment of the invention determines the first splicing height set of each wind profile radar in different detection modes, and a second splicing height set among a plurality of wind profile radars with different detection ranges, and a first splicing height for splicing and fusing detection data of different detection modes of each wind profile radar and a second splicing height for splicing and fusing detection data of the plurality of wind profile radars with different detection ranges are selected from the first splicing height set and the second splicing height set, the wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height are obtained, the wind profile radars with different detection modes and different detection ranges are objectively fused with the detection data of the same wind field, and randomness and uncertainty of artificially and subjectively setting the splicing height are avoided.
Example 3
The present embodiment provides an electronic device having a wind field detection data fusion apparatus shown in fig. 8.
Referring to fig. 9, fig. 9 is a schematic structural diagram of an electronic device according to an alternative embodiment of the present invention, and as shown in fig. 9, the electronic device may include: at least one processor 71, such as a CPU (Central Processing Unit), at least one communication interface 73, memory 74, at least one communication bus 72. Wherein a communication bus 72 is used to enable the connection communication between these components. The communication interface 73 may include a Display (Display) and a Keyboard (Keyboard), and the optional communication interface 73 may also include a standard wired interface and a standard wireless interface. The Memory 74 may be a high-speed RAM Memory (volatile Random Access Memory) or a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The memory 74 may alternatively be at least one memory device located remotely from the processor 71. Wherein the processor 71 may be in connection with the apparatus described in fig. 8, an application program is stored in the memory 74, and the processor 71 calls the program code stored in the memory 74 for performing any of the above-mentioned method steps.
The communication bus 72 may be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The communication bus 72 may be divided into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one thick line is shown in FIG. 9, but this does not indicate only one bus or one type of bus.
The memory 74 may include a volatile memory (RAM), such as a random-access memory (RAM); the memory may also include a non-volatile memory (english: non-volatile memory), such as a flash memory (english: flash memory), a hard disk (english: hard-drive, abbreviated: HDD) or a solid-state drive (english: SSD); the memory 74 may also comprise a combination of memories of the kind described above.
The processor 71 may be a Central Processing Unit (CPU), a Network Processor (NP), or a combination of CPU and NP.
The processor 71 may further include a hardware chip, which may be an application-specific integrated circuit (ASIC), a programmable logic device (CP L D), or a combination thereof, and the P L D may be a complex programmable logic device (CP L D), a field-programmable gate array (FPGA), a general-purpose array logic (GA L), or any combination thereof.
Optionally, the memory 74 is also used for storing program instructions. The processor 71 may call a program instruction to implement the wind field detection data fusion method as described in embodiment 1 of the present application.
Example 4
The present embodiments provide a computer-readable storage medium having stored thereon computer-executable instructions that may perform the wind farm detection data fusion method in any of the above method embodiments. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD), a Solid State Drive (SSD), or the like; the storage medium may also comprise a combination of memories of the kind described above.
Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope defined by the appended claims.

Claims (10)

1. A wind field detection data fusion method is characterized by comprising the following steps:
determining a first splicing height set of different detection modes of each wind profile radar; wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar;
determining a second set of stitching heights between a plurality of wind profile radars having different detection ranges; the second splicing height set comprises splicing heights for splicing detection data of a plurality of wind profile radars in different detection ranges;
acquiring detection data of the plurality of wind profile radars in the same wind field;
selecting a corresponding first splicing height from the first splicing height set, splicing and fusing detection data of different detection modes of each wind profile radar, wherein the detection data after splicing and fusing take the detection data at the corresponding first splicing height as boundary points;
and selecting a corresponding second splicing height from the second splicing set, splicing and fusing the detection data of the wind profile radars in different detection ranges, and obtaining the wind field detection data of the wind profile radars from the minimum detection height to the maximum detection height by using the detection data at the corresponding second splicing height as boundary points of the spliced and fused detection data.
2. The method of claim 1, further comprising:
acquiring detection data of a wind speed and wind direction sensor;
preprocessing the detection data of the wind speed and wind direction sensor;
and performing linear interpolation on the preprocessed detection data by using the detection data of the wind profile radar with the lowest detection height, and performing fusion calibration on the detection data of the wind profile radar with the lowest detection height to obtain the detection data obtained by fusing the detection data of the wind speed and direction sensor and the detection data of the wind profile radar with the lowest detection height.
3. The method of claim 1, wherein determining the first set of splice heights for each of the wind profile radar different detection modes comprises:
acquiring a data acquisition rate of each height layer of each detection mode of each wind profile radar, wherein each wind profile radar has n detection modes, and n is greater than 1;
selecting η the data acquisition rate is above the threshold in the detection range shared by the ith detection mode and the (i + 1) th detection mode1Height set A ofi(ii) a Taking i from 1 to n-1;
calculating the height set AiThe wind measurement error is the square sum of the average error and the standard deviation of the detection data of the ith detection mode and the detection data of the (i + 1) th detection mode;
selecting the maximum height h in the M height layers with the minimum wind measurement errori；
Adding 1 to the value of i, and returning to execute that the data acquisition rate is higher than the threshold η in the detection range shared by the ith detection mode and the (i + 1) th detection mode1Height set A ofiUntil h is obtained1,…,hn-1Forming the first set of splice heights.
4. The method of claim 1, wherein determining a second set of splice heights between a plurality of wind profile radars with different detection ranges comprises:
arranging the plurality of wind profile radars according to the sequence of the maximum detection height from small to large, and selecting half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]J is 1 to m-1, and m is the total number of the wind profile radar;
according to historical data, the fact that the data acquisition rate in the detection height range shared by the jth wind profile radar and the jth +1 wind profile radar is higher than a threshold η is counted2Height set B ofj；
Calculating the height set BjWind measurement error of each height layer; the wind measurement error is the square sum of the average error and the standard deviation of the horizontal wind speed detected by the jth wind profile radar and the jth +1 wind profile radar and the standard deviation of the horizontal wind direction;
selecting the maximum height H in the N height layers with the minimum wind measurement errorj；
Adding 1 to the value of j, and returning to execute the selection of half of the maximum detection height of the jth wind profile radar as Hj1And half the maximum detection height of the j +1 th wind profile radar is denoted as Hj2Forming a splicing height range [ H ] of the jth wind profile radar and the jth +1 wind profile radarj1,Hj2]Until H is obtained1,…,Hm-1And forming the second splicing height set.
5. The method according to claim 2, wherein the preprocessing of the detection data of the wind speed and wind direction sensor comprises:
carrying out primary smoothing on the detection data of the wind speed and direction sensor according to a vector averaging method to obtain detection data after primary smoothing;
and performing secondary smoothing on the detection data subjected to the primary smoothing by taking the time interval generated by the detection data of the wind profile radar as a smoothing window according to a vector averaging method to obtain the preprocessed detection data.
6. The method according to claim 2, wherein the performing linear interpolation on the preprocessed detection data by using the detection data of the wind profile radar with the lowest detection height and performing fusion calibration on the detection data of the wind profile radar with the lowest detection height comprises:
decomposing the preprocessed detection data into u with orthogonal directionsd、vdTwo sets of horizontal wind velocity components;
decomposing the detection data of the wind profile radar with the lowest detection height into U, V two groups of horizontal wind speed components with orthogonal directions;
and u is obtained by decomposing the data of the minimum height layer of which the data acquisition rate of the wind profile radar with the lowest detection height is higher than a preset threshold valuec、vcTwo sets of horizontal wind velocity components and ud、vdTaking the two groups of horizontal wind speed components as a reference, carrying out linear difference according to the layering height of the wind profile radar with the lowest detection height to obtain two groups of data of U 'and V', and synthesizing the two groups of data of U 'and V' to obtain wind speed and wind direction data of each height layer after linear interpolation;
comparing the wind speed and wind direction data of each height layer with the original wind speed and wind direction data detected by the wind profile radar with the lowest detection height;
and if the absolute value of the difference of the wind speeds or the absolute value of the difference of the wind directions is larger than a preset threshold value, replacing the original wind speed and wind direction data of the height layer of the wind profile radar with the lowest detection height with the wind speed and wind direction data of the height layer after linear interpolation.
7. The method of claim 6, wherein the predetermined threshold is calculated by the following formula:
Figure FDA0002458086210000031
wherein Th is a reference value of the preset threshold value, ThiA predetermined threshold value, h, representing the ith hierarchy leveliDenotes the ith hierarchy height, hdownIndicating the height, h, of the anemometric sensorupRepresenting said windProfile radar uc、vcThe corresponding stratification height, η, represents the adjustment factor.
8. A wind field detection data fusion device, comprising:
the first determining module is used for determining a first splicing height set of different detection modes of each wind profile radar; wherein each wind profile radar has a plurality of detection modes, and the first set of stitching heights comprises stitching heights for stitching detection data of different detection modes of the same wind profile radar;
a second determination module to determine a second set of stitching heights between a plurality of wind profile radars having different detection ranges; the second splicing height comprises a splicing height for splicing detection data of a plurality of wind profile radars in different detection ranges;
the acquisition module is used for acquiring detection data of the plurality of wind profile radars in the same wind field;
the first fusion module is used for selecting a corresponding first splicing height from the first splicing height set, splicing and fusing the detection data of each wind profile radar in different detection modes, and the spliced and fused detection data takes the detection data at the corresponding first splicing height as a boundary point;
and the second fusion module is used for selecting a corresponding second splicing height from the second splicing set, splicing and fusing the detection data of the wind profile radars in different detection ranges, and obtaining the wind field detection data of the plurality of wind profile radars from the minimum detection height to the maximum detection height by using the detection data at the corresponding second splicing height as boundary points in the spliced and fused detection data.
9. An electronic device, comprising:
a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing therein computer instructions, and the processor executing the computer instructions to perform the wind farm detection data fusion method according to any one of claims 1 to 7.
10. A computer-readable storage medium having stored thereon computer instructions for causing a computer to perform the wind farm detection data fusion method of any one of claims 1-7.
CN202010314980.8A 2020-04-20 2020-04-20 Wind field detection data fusion method and device and electronic equipment Active CN111487596B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010314980.8A CN111487596B (en) 2020-04-20 2020-04-20 Wind field detection data fusion method and device and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010314980.8A CN111487596B (en) 2020-04-20 2020-04-20 Wind field detection data fusion method and device and electronic equipment

Publications (2)

Publication Number Publication Date
CN111487596A true CN111487596A (en) 2020-08-04
CN111487596B CN111487596B (en) 2022-06-21

Family

ID=71813697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010314980.8A Active CN111487596B (en) 2020-04-20 2020-04-20 Wind field detection data fusion method and device and electronic equipment

Country Status (1)

Country Link
CN (1) CN111487596B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006133203A (en) * 2004-11-09 2006-05-25 Nec Engineering Ltd Wind profiler system
US20080169975A1 (en) * 2007-01-12 2008-07-17 Young Paul Yee Process for generating spatially continuous wind profiles from wind profiler measurements
JP2014121993A (en) * 2012-12-21 2014-07-03 Jtekt Corp Electric power steering system
CN104345312A (en) * 2014-10-27 2015-02-11 北京无线电测量研究所 Data fusion method and system of millimeter-wave cloud radar
CN105607063A (en) * 2016-01-05 2016-05-25 北京无线电测量研究所 Detection method and system for low-altitude wind shear at airport
US9535158B1 (en) * 2013-11-21 2017-01-03 Rockwell Collins, Inc. Weather radar system and method with fusion of multiple weather information sources
CN109870692A (en) * 2019-04-16 2019-06-11 浙江力邦合信智能制动系统股份有限公司 A kind of radar viewing system and data processing method
US10408939B1 (en) * 2019-01-31 2019-09-10 StradVision, Inc. Learning method and learning device for integrating image acquired by camera and point-cloud map acquired by radar or LiDAR corresponding to image at each of convolution stages in neural network and testing method and testing device using the same

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006133203A (en) * 2004-11-09 2006-05-25 Nec Engineering Ltd Wind profiler system
US20080169975A1 (en) * 2007-01-12 2008-07-17 Young Paul Yee Process for generating spatially continuous wind profiles from wind profiler measurements
JP2014121993A (en) * 2012-12-21 2014-07-03 Jtekt Corp Electric power steering system
US9535158B1 (en) * 2013-11-21 2017-01-03 Rockwell Collins, Inc. Weather radar system and method with fusion of multiple weather information sources
CN104345312A (en) * 2014-10-27 2015-02-11 北京无线电测量研究所 Data fusion method and system of millimeter-wave cloud radar
CN105607063A (en) * 2016-01-05 2016-05-25 北京无线电测量研究所 Detection method and system for low-altitude wind shear at airport
US10408939B1 (en) * 2019-01-31 2019-09-10 StradVision, Inc. Learning method and learning device for integrating image acquired by camera and point-cloud map acquired by radar or LiDAR corresponding to image at each of convolution stages in neural network and testing method and testing device using the same
CN109870692A (en) * 2019-04-16 2019-06-11 浙江力邦合信智能制动系统股份有限公司 A kind of radar viewing system and data processing method

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
DIETER NAGEL: "MPRF waveform for ground based radars to suppress returns from windmills and rain clutter", 《2015 SENSOR DATA FUSION: TRENDS, SOLUTIONS, APPLICATIONS (SDF)》 *
NATIVI, S (NATIVI, S): "SSM/I and the Montagnana radar data fusion to analyse clouds over Tuscany", 《REMOTE SENSING OF CLOUDS AND THE ATMOSPHERE V》 *
吴心筱: "《视频中人的动作分析与识别》", 30 September 2019, 北京理工大学出版社 *
汪学渊: "丘陵地区边界层风廓线雷达数据统计特性分析", 《气象科学》 *
王柳柳: "毫米波云雷达冻雨-降雪微物理和动力特征分析", 《气象》 *
郭忠立: "双流机场风切变预警系统方案设计及LLWAS算法研究", 《气象水文海洋仪器》 *
黄立峰等: "大气遥感双视场米散射激光雷达信号拼接", 《遥感学报》 *

Also Published As

Publication number Publication date
CN111487596B (en) 2022-06-21

Similar Documents

Publication Publication Date Title
CN106950614B (en) A kind of region automatic weather station hour rainfall data method of quality control
CN105203023B (en) A kind of one-stop scaling method of vehicle-mounted three-dimensional laser scanning system placement parameter
CN105954543B (en) A kind of current meter high-speed straight-line formula integral checking method and its calibration equipment
CN114114198B (en) Precipitation data quality control method and device, storage medium and equipment
CN110071692A (en) Photovoltaic module failure judgment method, device and controller
CN110261857A (en) A kind of weather radar spatial interpolation methods
CN101231167A (en) Method for detecting and regulating sea survey line net systematical error
US11145082B2 (en) Method for measuring antenna downtilt angle based on deep instance segmentation network
CN115854999A (en) H-ADCP section average flow velocity self-correction method based on scene self-adaptation
CN116222670B (en) Ecological landscape slope monitoring method for urban green land planning
CN102426352B (en) Wind profiling radar based wind calculation method
CN109740684A (en) Shared bicycle lairage detection method and device
CN107831516B (en) Method for acquiring real-time high-precision displacement of dam by fusing GNSS and ground monitoring network
CN100478650C (en) Method and apparatus for improved fault detection in power generation equipment
CN109254291A (en) The Doppler power spectra data fusion method and device of cloud radar
CN111487596B (en) Wind field detection data fusion method and device and electronic equipment
CN110426717A (en) It is a kind of based on the co-located method and system of three-dimensional map assisted GNSS, positioning device, storage medium
CN106526651B (en) The method for building up and system of a kind of detector crystal position table
CN106501815A (en) A kind of extraterrestrial target orbit maneuver fusion detection method of only Space-based Angle Measured tracking
Wardah et al. Radar rainfall estimates comparison with kriging interpolation of gauged rain
EP4307172A1 (en) Data interpolation method and apparatus and storage medium
CN115015904A (en) Speed measurement method, system and equipment based on fusion of radar and camera and storage medium
CN110703012B (en) Distributed fault diagnosis method for power transmission line
CN109709555B (en) Method and system for identifying difference of adjacent scan data of weather radar
Helmert et al. An operational tool to quality control 2D radar reflectivity data for assimilation in COSMO-DE

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant