WIND SPEED FORECASTING DEVICE, WIND SPEED FORECASTING METHOD AND RADAR DEVICE

Description

TECHNICAL FIELD

[0001] The present disclosure relates to a wind speed prediction device, and a wind speed prediction method.BACKGROUND ART

[0002] Some wind speed prediction devices that predict a wind speed in an observation region include a wind speed prediction device (hereinafter referred to as a "conventional wind speed prediction device") that predicts a wind speed in an observation region by simulating the wind speed using a weather model (see Patent Literature 1). The simulation of the wind speed is three-dimensional calculation processing.

[0003] Patent Literature 2 discloses a method for determining the flow of a fluid in a volume of interest, including steps of remotely measuring, at a plurality of measurement points distributed along at least three axes having different spatial orientations passing through the volume of interest, the radial velocity of the fluid in the vicinity of the measurement points, and for calculating the velocity of the fluid at a plurality of calculation points distributed in a grid in the volume of interest, wherein the calculation of the velocity of the fluid includes the use of a mechanical behavior model of the fluid.

[0004] Non-Patent Literature 1 briefly reviews procedures employed to retrieve vector wind information from single-Doppler radar observations. In particular, those procedures implying a linearity hypothesis for the wind field are shown to be particular cases of a method termed VVP (Volume Velocity Processing). This method, which makes full use of radar velocity data filling a volume, is first tested on simulated observations in order to assess its accuracy, then applied to actual data and shown to yield unbiased parameters of the horizontal vector wind field, as well as an estimate of the hydrometeor fall velocity and several control parameters.

[0005] In Non-Patent Literature 2, first, a Doppler velocity dealiasing technique is developed for difficult aliasing problems. A dealiasing algorithm is described that automatically processes one radial at a time by comparing the radial with a previous or surrounding radial. This technique has worked reliably on numerous Doppler radar datasets for severe storm environments observed by an X-band Doppler radar. Second, an analysis technique for the retrieval of the mesoscale (20 -100 km) wind field, such as the VVP ( Volume Velocity Processing ), is extended. This extended VVP analysis has adaptability for a nonlinear wind field better than the original method. This algorithm can compute both estimates and residual variances simultaneously for several analysis volumes, and select the optimal horizontal wind and its divergence which minimizes the residual variance. The horizontal wind field estimated has a larger scale of the observation resolution and this estimate aids in tracking the entire movement of rainfall distribution.

[0006] Patent Literature 3 discloses a system and method for determining wind velocity. The LIDAR 104 transmits a plurality of laser pulses. The paths of the laser pulses surround a center line that does not coincide with a vertical axis intersecting the LIDAR , i.e., an off-vertical scan. Each laser pulse reflects off particles in the air at a given location. This location is the intersection of the path of each LIDAR pulse and the target altitude. The target altitude is the altitude at which the wind velocity is to be determined. The reflected pulses are received by the LIDAR. The invention groups the received LIDAR pulses according to the altitude at which they were reflected, as opposed to the distance each pulse travels. Grouping the received pulses in this manner enables the invention to perform an off-vertical scan that accurately determines the wind velocity at a target altitude. The received laser pulses have an associated Doppler shift. The Doppler shift is a measure of the relative radial wind velocity at the target location. The invention measures the Doppler shift from a plurality of received laser pulses reflected from the target altitude. These measurements enable the invention to determine the wind velocity at the target location.CITATION LISTPATENT LITERATURES AND NON-PATENT LITERATURES

[0007] Patent Literature 1: Japanese Patent Application Laid-Open No. 2014-202190 Patent Literature 2: US 2013 / 282285 A1 Patent Literature 3: US 5 724 125 A Non-Patent Literature 1: Waldteufel PET AL: "On the Analysis of Single-Doppler Radar Data", Journal of Applied Meteorology and Climatology, 1 April 1979 (1979-04-01), pages 5382-542, XP093171238, DOI: https: / / doi.org / No-Patent Literature 2: Soichiro Sugimoto ET AL: "A rainfall prediction method using a meteorological radar: Estimation of horizontal wind and its divergence fields using Doppler radar data", Report of the Central Research Institute of the Electric Power Industry (Report of the Central Research Institute of the Electric Power Industry), 1 April 2001 (2001-04-01), pages 1-29, XP093171258, SUMMARY OF INVENTIONTECHNICAL PROBLEM

[0008] The conventional wind speed prediction device disclosed in Patent Literature 1 performs three-dimensional calculation processing, and thus has a problem that a calculation load may increase.

[0009] The present disclosure has been made to solve the above problems, and an object thereof is to obtain a wind speed prediction device and a wind speed prediction method capable of predicting a wind speed without performing three-dimensional calculation processing.SOLUTION TO PROBLEM

[0010] A wind speed prediction device according to the present disclosure is set forth in claim 1. A wind speed prediction method according to the present disclosure is set forth in claim 6.ADVANTAGEOUS EFFECTS OF INVENTION

[0011] According to the present disclosure, a wind speed can be predicted without performing three-dimensional calculation processing.BRIEF DESCRIPTION OF DRAWINGS

[0012] FIG. 1 is a configuration diagram illustrating a radar device 1 including a wind speed prediction device 3 according to a first embodiment. FIG. 2A is an explanatory diagram illustrating a plurality of beams in a case where the position of a flight area is directly above a beam transmitting and receiving unit 2 in a vertical direction, and FIG. 2B is an explanatory diagram illustrating a plurality of beams in a case where the position of the flight area is deviated from directly above the beam transmitting and receiving unit 2 in the vertical direction. FIG. 3 is a configuration diagram illustrating a wind speed prediction device 3 according to the first embodiment. FIG. 4 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the first embodiment. FIG. 5 is a hardware configuration diagram of a computer in a case where the wind speed prediction device 3 is implemented by software, firmware, or the like. FIG. 6 is a flowchart illustrating a wind speed prediction method which is a processing procedure performed by the wind speed prediction device 3. FIG. 7A is an explanatory diagram illustrating that altitudes of a plurality of observation points of the same range bin Rb m are different from each other when an elevation angle θ n is different, and FIG. 7B is an explanatory diagram illustrating that a plurality of observation points of the same range bin Rb m are at the same altitude even when the elevation angle θ n is changed. FIG. 8 is an explanatory diagram illustrating positions (x, y, z) of observation points of the wind speed. FIG. 9 is a configuration diagram illustrating a wind speed prediction device 3 according to a second embodiment. FIG. 10 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the second embodiment. FIG. 11 is a configuration diagram illustrating a wind speed prediction device 3 according to a third embodiment. FIG. 12 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the third embodiment. DESCRIPTION OF EMBODIMENTS

[0013] Hereinafter, in order to describe the present disclosure in more detail, modes for carrying out the present disclosure will be described with reference to the accompanying drawings.First Embodiment

[0014] FIG. 1 is a configuration diagram illustrating a radar device 1 including a wind speed prediction device 3 according to a first embodiment.

[0015] The radar device 1 illustrated in FIG. 1 includes a beam transmitting and receiving unit 2 and a wind speed prediction device 3.

[0016] FIG. 2 is an explanatory diagram illustrating a plurality of beams emitted to space from the beam transmitting and receiving unit 2. In the example of FIG. 2, the observation region is present in the space. The observation region is, for example, a flight area in which a flying object is scheduled to fly. The flying object is, for example, a drone or a helicopter.

[0017] FIG. 2A is an explanatory diagram illustrating a plurality of beams in a case where the position of the flight area is directly above the beam transmitting and receiving unit 2 in the vertical direction.

[0018] FIG. 2B is an explanatory diagram illustrating a plurality of beams in a case where the position of the flight area is deviated from directly above the beam transmitting and receiving unit 2 in the vertical direction.

[0019] The beam transmitting and receiving unit 2 includes, for example, a beam generating unit, a beam transmitter, a radiation direction switching unit, an antenna, and a beam receiver. In FIG. 1, notations of the beam generating unit, the beam transmitter, the radiation direction switching unit, the antenna, and the beam receiver are omitted.

[0020] As illustrated in FIG. 2, the beam transmitting and receiving unit 2 emits each of a plurality of beams having mutually different elevation angles θ n (n = 1, ..., N), which are angles formed by the line-of-sight direction and the horizontal direction, into space. N is an integer of 2 or more.

[0021] That is, as illustrated in FIG. 2A, when the flight area is present directly above in the vertical direction, the beam transmitting and receiving unit 2 emits beams in different directions at respective elevation angles θ n by performing beam scanning at the respective elevation angles θ n .

[0022] In FIG. 2A, N = 3, and each of the beam with the elevation angle θ 1 , the beam with the elevation angle θ 2 , and the beam with the elevation angle θ 3 is emitted from the beam transmitting and receiving unit 2 to space. In FIG. 2A, beams are emitted in four directions at the respective elevation angles θ n . In FIG. 2A, black circles indicate observation points of the wind speed.

[0023] Further, as illustrated in FIG. 2B, when the position of the flight area is deviated from directly above in the vertical direction, the beam transmitting and receiving unit 2 emits beams at respective elevation angles θ n while switching the elevation angles θ n .

[0024] In FIG. 2B, N = 12, and beams of elevation angles θ 1 to θ 12 are emitted from the beam transmitting and receiving unit 2 to space. For simplification of the drawing, only the elevation angles θ 1 to θ 6 are illustrated in FIG. 2B, and the illustration of the elevation angles θ 7 to θ 12 is omitted. In FIG. 2B, a beam is emitted in one direction at each elevation angle θ n . In FIG. 2B, black circles indicate observation points of the wind speed.

[0025] Each of FIGS. 2A and 2B illustrates an example of obtaining wind speed distributions in three two-dimensional planes. However, this is merely an example, and wind speed distributions in one or two two-dimensional planes may be obtained, or wind speed distributions in four or more two-dimensional planes may be obtained.

[0026] As the beam emitted from the beam transmitting and receiving unit 2, for example, in addition to continuous wave (CW) pulsed light, frequency-modulated pulsed light can be used.

[0027] The beam transmitting and receiving unit 2 receives each beam scattered in space as a scattering signal.

[0028] The beam transmitting and receiving unit 2 converts each scattering signal from an analog signal to a digital signal, and outputs the digital signal to the wind speed prediction device 3.

[0029] Further, the beam transmitting and receiving unit 2 outputs angle information indicating the elevation angle θ n of each beam emitted to space to the wind speed prediction device 3.

[0030] In the radar device 1 illustrated in FIG. 1, it is assumed that the beam transmitting and receiving unit 2 includes a beam generating unit, a beam transmitter, a radiation direction switching unit, an antenna, and a beam receiver. However, the beam transmitting and receiving unit 2 is only required to emit each of a plurality of beams having different elevation angles θ n into space and receive each beam scattered in the space as a scattering signal, and may have any configuration.

[0031] FIG. 3 is a configuration diagram illustrating the wind speed prediction device 3 according to the first embodiment.

[0032] FIG. 4 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the first embodiment.

[0033] The wind speed prediction device 3 illustrated in FIG. 3 includes a scattering signal acquiring unit 11, a Doppler frequency calculating unit 12, a first wind speed distribution estimating unit 13, a second wind speed distribution estimating unit 14, and a wind speed prediction unit 15.

[0034] The scattering signal acquiring unit 11 is implemented by, for example, a scattering signal acquiring circuit 21 illustrated in FIG. 4.

[0035] The scattering signal acquiring unit 11 acquires each of digital signals s θn,j (t) from the beam transmitting and receiving unit 2 as each of scattering signals after being scattered in space. j = 1, ..., J. j is the scan number of beam scanning at the elevation angle θ n . In the example of FIG. 2A, since beams are emitted in four directions at each elevation angle θ n , J = 4. In the example of FIG. 2B, since the beam is emitted in one direction at each elevation angle θ n , J = 1. t is a variable indicating time.

[0036] Further, the scattering signal acquiring unit 11 acquires the angle information indicating the elevation angle θ n of each beam emitted into space by the beam transmitting and receiving unit 2.

[0037] The scattering signal acquiring unit 11 outputs each of the digital signals s θn,j (t) and the angle information to the Doppler frequency calculating unit 12.

[0038] The Doppler frequency calculating unit 12 is implemented by, for example, a Doppler frequency calculating circuit 22 illustrated in FIG. 4.

[0039] The Doppler frequency calculating unit 12 acquires each of the digital signals s θn,j (t) and the angle information from the scattering signal acquiring unit 11.

[0040] Before performing fast Fourier transformation (FFT) on each digital signal s θn,j (t) in the hit direction, the Doppler frequency calculating unit 12 sets a range bin width as the distance resolution of each digital signal s θn,j (t) according to the elevation angle θ n indicated by the angle information. That is, the Doppler frequency calculating unit 12 sets an FFT length L θn corresponding to the range bin width.

[0041] The FFT length L θn is set shorter as the elevation angle θ n of the beam emitted to space is wider, and is set longer as the elevation angle θ n of the beam emitted to the space is narrower.

[0042] In the example of FIG. 2A, the FFT length L θ1 of the digital signal s θ1,j (t) corresponding to the beam at the elevation angle θ 1 is longer than the FFT length L θ2 of the digital signal s θ2,j (t) corresponding to the beam at the elevation angle θ 2 , and the FFT length L θ2 of the digital signal s θ2 (t) corresponding to the beam at the elevation angle θ 2 is longer than the FFT length L θ3 of the digital signal s θ3,j (t) corresponding to the beam at the elevation angle θ 3 .

[0043] The Doppler frequency calculating unit 12 converts each of the digital signals s θn,j (t) into frequency domain signals r θn,j (f) by performing FFT on the digital signals s θn,j (t) acquired by the scattering signal acquiring unit 11 in the hit direction. f is a variable indicating a frequency.

[0044] The Doppler frequency calculating unit 12 calculates the respective Doppler frequencies dpf θn,j,Rbm in a plurality of range bins Rb 1 to Rb M from the signals r θn,j (f) in the frequency domain. m = 1, ..., M, and M is an integer of 1 or more.

[0045] FIG. 2 illustrates an example of obtaining wind speed distributions in three two-dimensional planes. As illustrated in FIG. 2, in a case of obtaining wind speed distributions in three two-dimensional planes, the Doppler frequency calculating unit 12 calculates the Doppler frequencies dpf θn,j,Rbm in the range bins Rb 1 , Rb 2 , and Rb 3 from the signals r θn,j (f) in the frequency domain.

[0046] In the wind speed prediction device 3 illustrated in FIG. 3, it is assumed that M is an integer of 2 or more. Here, M is not limited to an integer of 2 or more, and M = 1 may be set.

[0047] When M = 1, the Doppler frequency calculating unit 12 calculates only the Doppler frequency dpf θn,j,Rb1 of one range bin Rb 1 including the observation region from the frequency domain signal r θn,j (f).

[0048] The Doppler frequency calculating unit 12 outputs the Doppler frequencies dpf θn,j,Rbm of the respective range bins Rb m to the first wind speed distribution estimating unit 13.

[0049] The first wind speed distribution estimating unit 13 is implemented by, for example, a first wind speed distribution estimating circuit 23 illustrated in FIG. 4.

[0050] The first wind speed distribution estimating unit 13 estimates a wind speed distribution u Rbm (x, y) in a two-dimensional plane corresponding to each range bin Rb m using a volume velocity processing (VVP) method from a plurality of Doppler frequencies dpf θ1,j,Rbm to dpf θN,j,Rbm in each range bin Rb m calculated by the Doppler frequency calculating unit 12.

[0051] When M is 1, the first wind speed distribution estimating unit 13 estimates the wind speed distribution u Rb1 (x, y) in the two-dimensional plane corresponding to the range bin Rb 1 as the two-dimensional plane including the observation region, and outputs the wind speed distribution u Rb1 (x, y) in the two-dimensional plane to the wind speed prediction unit 15.

[0052] When M is 2 or more, the first wind speed distribution estimating unit 13 outputs the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m to the second wind speed distribution estimating unit 14.

[0053] The second wind speed distribution estimating unit 14 is implemented by, for example, a second wind speed distribution estimating circuit 24 illustrated in FIG. 4.

[0054] The second wind speed distribution estimating unit 14 estimates a wind speed distribution u Rbm' (x, y) in a two-dimensional plane between a plurality of two-dimensional planes from the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes estimated by the first wind speed distribution estimating unit 13. For example, the wind speed distribution u Rb1' (x, y) in the two-dimensional plane is a wind speed distribution in the two-dimensional plane between the two-dimensional plane corresponding to the range bin Rb 1 and the two-dimensional plane corresponding to the range bin Rb 2 . The wind speed distribution u Rb2' (x, y) in the two-dimensional plane is a wind speed distribution in the two-dimensional plane between the two-dimensional plane corresponding to the range bin Rb 2 and the two-dimensional plane corresponding to the range bin Rb 3 .

[0055] The second wind speed distribution estimating unit 14 outputs the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional plane between the plurality of two-dimensional planes to the wind speed prediction unit 15.

[0056] The wind speed prediction unit 15 is implemented by, for example, a wind speed prediction circuit 25 illustrated in FIG. 4.

[0057] When M is 1, the wind speed prediction unit 15 predicts a wind speed u(t) in the observation region from the wind speed distribution u Rb1 (x, y) in the two-dimensional plane estimated by the first wind speed distribution estimating unit 13 using the two-dimensional Navier-Stokes equation (hereinafter referred to as a "two-dimensional N-S equation"). The observation region is, for example, an area including a flight area in which the flying object is scheduled to fly.

[0058] When M is 2 or more, the wind speed prediction unit 15 selects the wind speed distribution in the two-dimensional plane including the flight area as the observation region from the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes estimated by the first wind speed distribution estimating unit 13 and the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional plane estimated by the second wind speed distribution estimating unit 14.

[0059] The wind speed prediction unit 15 predicts the wind speed u(t) in the flight area from the wind speed distribution in the selected two-dimensional plane using the two-dimensional N-S equation.

[0060] In FIG. 3, it is assumed that each of the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, and the wind speed prediction unit 15, which are components of the wind speed prediction device 3, is implemented by dedicated hardware as illustrated in FIG. 4. That is, it is assumed that the wind speed prediction device 3 is implemented by the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, and the wind speed prediction circuit 25.

[0061] Each of the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, and the wind speed prediction circuit 25 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination thereof.

[0062] The components of the wind speed prediction device 3 are not limited to those implemented by dedicated hardware, and the wind speed prediction device 3 may be implemented by software, firmware, or a combination of software and firmware.

[0063] Software or firmware is stored in a memory of the computer as a program. The computer means hardware that executes a program, and corresponds to, for example, a central processing unit (CPU), a central processing device, a processing device, an arithmetic device, a microprocessor, a microcomputer, a processor, or a digital signal processor (DSP).

[0064] FIG. 5 is a hardware configuration diagram of a computer in a case where the wind speed prediction device 3 is implemented by software, firmware, or the like.

[0065] In a case where the wind speed prediction device 3 is implemented by software, firmware, or the like, a program for causing a computer to execute each processing procedure performed in the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, and the wind speed prediction unit 15 is stored in a memory 31. Then, a processor 32 of the computer executes the program stored in the memory 31.

[0066] Further, FIG. 4 illustrates an example in which each of the components of the wind speed prediction device 3 is implemented by dedicated hardware, and FIG. 5 illustrates an example in which the wind speed prediction device 3 is implemented by software, firmware, or the like. However, this is merely an example, and some components in the wind speed prediction device 3 may be implemented by dedicated hardware, and the remaining components may be implemented by software, firmware, or the like.

[0067] Next, the operation of the radar device 1 illustrated in FIG. 1 will be described.

[0068] FIG. 6 is a flowchart illustrating a wind speed prediction method which is a processing procedure performed by the wind speed prediction device 3.

[0069] As illustrated in FIG. 2, the beam transmitting and receiving unit 2 emits each of a plurality of beams having different elevation angles θ n (n = 1, ..., N) to space.

[0070] That is, as illustrated in FIG. 2A, when the flight area is present directly above in the vertical direction, the beam transmitting and receiving unit 2 emits beams in different directions at respective elevation angles θ n by performing beam scanning at the respective elevation angles θ n .

[0071] As illustrated in FIG. 2B, when the position of the flight area is deviated from directly above in the vertical direction, the beam transmitting and receiving unit 2 emits beams at respective elevation angles θ n while switching the elevation angles θ n . The beam emitted from the beam transmitting and receiving unit 2 is, for example, pulsed light.

[0072] Each beam emitted from the beam transmitting and receiving unit 2 is backscattered by fine particles such as aerosol floating in space. For example, when a plurality of fine particles is present in the line-of-sight direction of one beam, one beam is backscattered by each fine particle. The backscattered beam returns to the beam transmitting and receiving unit 2 as a scattering signal. Since the distances from the beam transmitting and receiving unit 2 to the plurality of fine particles are different from each other, the times at which the respective scattering signals return to the beam transmitting and receiving unit 2 are different from each other.

[0073] The beam transmitting and receiving unit 2 converts each scattering signal from an analog signal into a digital signal s θn,j (t) (n = 1, ..., N: j = 1, ..., J), and outputs the digital signal s θen,j to the wind speed prediction device 3.

[0074] Further, the beam transmitting and receiving unit 2 outputs the angle information indicating the elevation angle θ n of each beam emitted to space to the wind speed prediction device 3.

[0075] The scattering signal acquiring unit 11 acquires each of the digital signals s θn,j (t) from the beam transmitting and receiving unit 2 (step ST1 in FIG. 6).

[0076] Further, the scattering signal acquiring unit 11 acquires the angle information indicating the elevation angle θ n of each beam from the beam transmitting and receiving unit 2.

[0077] The scattering signal acquiring unit 11 outputs each of the digital signals s θn,j (t) and the angle information to the Doppler frequency calculating unit 12.

[0078] The Doppler frequency calculating unit 12 acquires each of the digital signals s θn,j (t) from the scattering signal acquiring unit 11, and acquires each angle information from the scattering signal acquiring unit 11.

[0079] Before FFT each of the digital signals s θn,j (t) in the hit direction, the Doppler frequency calculating unit 12 sets the FFT length L θn of each of the digital signals s θn,j (t) according to the elevation angle θ n indicated by the angle information (step ST2 in FIG. 6). The FFT length L θn corresponds to a range bin width of the scattering signal.

[0080] That is, the Doppler frequency calculating unit 12 sets the FFT length L θn of the digital signal s θn,j (t) in such a way that the range bin width corresponding to the FFT length L θn is proportional to 1 / sin (θ n ).

[0081] In a case where the FFT lengths L θn of the respective digital signals s θn,j (t) are the same, when the elevation angles θ n are different, as illustrated in FIG. 7A, the altitudes of a plurality of observation points of the same range bin Rb m are different from each other.

[0082] On the other hand, in a case where the FFT lengths L θn of the respective digital signals s θn,j (t) are set according to the elevation angle θ n , as illustrated in FIG. 7B, the altitudes of the plurality of observation points of the same range bin Rb m are the same altitude even if the elevation angles θ n are different.

[0083] In the radar device 1 illustrated in FIG. 1, the FFT length L θn is set shorter as the elevation angle θ n is wider, and the FFT length L θn is set longer as the elevation angle θ n of the scattering signal is narrower, so that a plurality of observation points in the same range bin Rb m have the same altitude.

[0084] FIG. 7 is an explanatory diagram illustrating an observation point of the Doppler frequency.

[0085] FIG. 7A is an explanatory diagram illustrating that the altitudes of a plurality of observation points of the same range bin Rb m are different from each other when the elevation angle θ n is different.

[0086] FIG. 7B is an explanatory diagram illustrating that a plurality of observation points in the same range bin Rb m have the same altitude even when the elevation angle θ n changes.

[0087] When the FFT length L θn is set to be short, the width of the range bin Rb m becomes narrow and the distance resolution becomes high. On the other hand, when the FFT length L θn is set to be long, the width of the range bin Rb m becomes wide and the distance resolution becomes low.

[0088] The Doppler frequency calculating unit 12 converts each digital signal s θn,j (t) into a frequency domain signal r θn,j (f) by performing FFT with an FFT length of L θn in the hit direction for each digital signal s θn,j (t) (step ST3 in FIG. 6).

[0089] The Doppler frequency calculating unit 12 calculates the respective Doppler frequencies dpf θn,j,Rbm in the plurality of range bins Rb 1 to Rb M from the signals r θn,j (f) in the frequency domain (step ST4 in FIG. 6).

[0090] The Doppler frequency calculating unit 12 outputs the Doppler frequencies dpf θn,j,Rbm of the respective range bins Rb m to the first wind speed distribution estimating unit 13.

[0091] Since the processing itself of calculating the Doppler frequencies dpf θn,j,Rbm of the range bin Rb m is a known technique, detailed description thereof will be omitted.

[0092] Note that calculation accuracy of the Doppler frequencies dpf θn,j,Rbm may be lowered due to a low signal to noise ratio (SNR) of the digital signal s θn,j (t). In such a case, the Doppler frequency calculating unit 12 may increase the calculation accuracy of the Doppler frequencies dpf θn,j,Rbm , for example, by performing the following processing.

[0093] For example, the Doppler frequency calculating unit 12 acquires digital signals s θn,j (t) related to a plurality of pulsed beams having the same elevation angle θ n and different hit numbers, and converts the digital signals s θn,j (t) related to the plurality of beams into frequency domain signals r θn,j (f). Then, the Doppler frequency calculating unit 12 coherently integrates the signals r θn,j (f) in the plurality of frequency domains to increase the calculation accuracy of the Doppler frequencies dpf θn,j,Rbm . Since the processing itself for enhancing the calculation accuracy of the Doppler frequencies dpf θn,j,Rbm is a known technique, detailed description thereof will be omitted.

[0094] The first wind speed distribution estimating unit 13 acquires the Doppler frequencies dpf θn,j,Rbm of the respective range bins Rb m from the Doppler frequency calculating unit 12.

[0095] The first wind speed distribution estimating unit 13 estimates the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m from the plurality of Doppler frequencies dpf θ1,j,Rbm to dpf θN,j,Rbm in each range bin Rb m using the VVP method (step ST5 in FIG. 6).

[0096] Hereinafter, estimation processing of the wind speed distribution u Rbm (x, y) in the two-dimensional plane by the first wind speed distribution estimating unit 13 will be specifically described.

[0097] FIG. 8 is an explanatory diagram illustrating positions (x, y, z) of observation points of the wind speed.

[0098] The position (x, y, z) of the observation point of the wind speed is expressed by polar coordinates centered on the position O of the beam transmitting and receiving unit 2 included in the radar device 1 as expressed by the following Expression (1). x y z = r cos θ sin ϕ r cos θ cos ϕ r sin θ

[0099] In Expression (1), θ is the elevation angle of a beam. φ is an angle formed by a line segment connecting the center point O' of the x-y plane where the observation point is present and the position (x, y, z) of the observation point and the y axis. r is a distance between the position O of the beam transmitting and receiving unit 2 and the position of the observation point (x, y, z).

[0100] For example, when the terrain of the ground surface vertically below the observation point is not a special terrain but a general terrain, there is no practical problem in observing the wind affecting the flying object as long as the wind speed prediction device 3 can observe the wind speed in a direction parallel to the two-dimensional plane even if the wind speed prediction device 3 cannot observe the wind speed in the vertical direction with respect to the two-dimensional plane where the observation point is present.

[0101] The general terrain corresponds to a flat terrain, an inclined terrain, a terrain having small irregularities, or the like. The special terrain corresponds to a terrain having a valley bottom, a terrain having a cliff with a large difference in elevation, or the like.

[0102] Accordingly, in the wind speed prediction device 3 illustrated in FIG. 3, the first wind speed distribution estimating unit 13 ignores the wind speed in the vertical direction, and a true wind speed vector u at the observation point is expressed as the following Expression (2). u = u x u y

[0103] In Expression (2), u x represents a component in a direction parallel to the x-axis of the wind speed vector u, and u y represents a component in a direction parallel to the y-axis of the wind speed vector u.

[0104] The first wind speed distribution estimating unit 13 performs Taylor expansion of the wind speed vector u centered on the center point O' with respect to x and y on the assumption that the wind speed u 0 at the center point O' of the x-y plane is expressed as the following Expression (3). For example, when the third and subsequent terms are ignored, the following approximate Expression (4) holds. u 0 = u x 0 , u y 0 T u x y ≃ u 0 + ∂ u ∂ x x + ∂ u ∂ y y + 1 2 ∂ 2 u ∂ x 2 x 2 + 2 ∂ 2 u ∂ x ∂ y xy + ∂ 2 u ∂ y 2 y 2

[0105] Since what is observed by the radar device 1 is the speed in the line-of-sight direction of the wind speed at the observation point, an observed wind speed u r (θ, φ) is expressed as the following Expression (5). u r θ ϕ = u x y ⋅ i x cos θ sin ϕ + i y cos θ cos ϕ + i z sin θ = a θ , ϕ T ⋅ p a θ , ϕ = cos θ sin ϕ , cos θ cos ϕ , rcos 2 θ sin 2 ϕ , rcos 2 θ sin ϕ cos ϕ , rcos 2 θ cos 2 ϕ , r 2 cos 3 θ sin 3 ϕ , r 2 cos 3 θ sin 2 ϕ cos ϕ , r 2 cos 3 θ sin ϕ cos 2 ϕ , r 2 cos 3 θ cos 3 ϕ p = u x 0 , u y 0 , ∂ u x ∂ x , ∂ u x ∂ y + ∂ u y ∂ x , ∂ u y ∂ y , ∂ 2 u x ∂ x 2 , ∂ 2 u x ∂ x ∂ y + ∂ 2 u y ∂ x 2 , ∂ 2 u x ∂ y 2 + ∂ 2 u y ∂ x ∂ y , ∂ 2 u y ∂ y 2

[0106] The first wind speed distribution estimating unit 13 acquires the observed wind speeds u r (θ 1 , φ 1 ) to u r (θ K , φ K ) of K beams, and substitutes the acquired observed wind speeds u r (θ 1 , φ 1 ) to u r (θ K , φ K ) into the following Expression (8) to estimate p illustrated in Expression (7). K is an integer of 2 or more. The K beams are K pieces of pulsed light in which sets of the elevation angle θ and the angle φ are different from each other. u r θ 1 ϕ 1 u r θ 2 ϕ 2 ⋮ u r θ K ϕ K = a θ 1 , ϕ 1 T a θ 2 , ϕ 2 T ⋮ a θ K , ϕ K T p ⇒ p ≃ a θ 1 , ϕ 1 T a θ 2 , ϕ 2 T ⋮ a θ K , ϕ K T † u r θ 1 ϕ 1 u r θ 2 ϕ 2 ⋮ u r θ K ϕ K

[0107] A symbol "†" denotes pseudo-inverse matrix in Expression 8 In Expression (8), † is a symbol representing a pseudo-inverse matrix.

[0108] Among the elements of p illustrated in Expression (7), the fourth, seventh, and eighth elements from the top are expressed as the sum of two parameters, and thus are uncertain elements.

[0109] The first wind speed distribution estimating unit 13 determines the uncertain elements by equally allocating the parameter of the element of p estimated by Expression (8).

[0110] When the fourth element from the top among the elements of p is p 4 , the first wind speed distribution estimating unit 13 determines p 4 as in the following Expression (9), for example.

[0111] When the seventh element from the top among the elements of p is p 7 , the first wind speed distribution estimating unit 13 determines p 7 as in the following Expression (10), for example.

[0112] When the eighth element from the top among the elements of p is p 8 , the first wind speed distribution estimating unit 13 determines p 8 as in the following Expression (11), for example. ∂ u x ∂ y = ∂ u y ∂ x = 1 2 p 4 ∂ 2 u x ∂ x ∂ y = ∂ 2 u y ∂ x 2 = 1 2 p 7 ∂ 2 u x ∂ y 2 = ∂ 2 u y ∂ x ∂ y = 1 2 p 8

[0113] As described above, since all the coefficients u 0 , ∂u / ∂x, ∂u / ∂y, ∂ 2< u / ∂x 2< , ∂ 2< u / ∂x∂y, and ∂ 2< u / ∂y 2< illustrated in the approximate Expression (4) are determined, the first wind speed distribution estimating unit 13 can obtain the wind speed vector u(x, y) at each observation point from the approximate Expression (4).

[0114] When the wind speed vector u(x, y) at each observation point can be obtained, the first wind speed distribution estimating unit 13 can obtain the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m .

[0115] Here, the first wind speed distribution estimating unit 13 approximates the Taylor expansion of the wind speed vector u(x, y) by 0th to 2nd order terms. However, this is merely an example, and the first wind speed distribution estimating unit 13 may approximate the Taylor expansion of the wind speed vector u(x, y) with terms of 0 degree to 3 degrees or more.

[0116] When M is 1, the first wind speed distribution estimating unit 13 outputs the wind speed distribution u Rb1 (x, y) in the two-dimensional plane corresponding to the range bin Rb 1 to the wind speed prediction unit 15.

[0117] When M is 2 or more, the first wind speed distribution estimating unit 13 outputs the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m to the second wind speed distribution estimating unit 14.

[0118] The second wind speed distribution estimating unit 14 acquires the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m from the first wind speed distribution estimating unit 13.

[0119] The second wind speed distribution estimating unit 14 estimates a wind speed distribution u Rbm' (x, y) in a two-dimensional plane between the plurality of two-dimensional planes from the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes (step ST6 in FIG. 6).

[0120] The second wind speed distribution estimating unit 14 outputs the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional planes to the wind speed prediction unit 15.

[0121] Hereinafter, estimation processing of the wind speed distribution u Rbm' (x, y) in the two-dimensional plane by the second wind speed distribution estimating unit 14 will be specifically described.

[0122] For example, when estimating the wind speed distribution u Rbm' (x, y) in the two-dimensional plane between the two-dimensional plane in which the wind speed distribution is u Rbm (x, y) and a two-dimensional plane in which the wind speed distribution is u Rbm+1 (x, y), the second wind speed distribution estimating unit 14 calculates the distance L m-m' in the vertical direction from the two-dimensional plane in which the wind speed distribution u Rbm' (x, y) is estimated to the two-dimensional plane in which the wind speed distribution is u Rbm (x, y).

[0123] Further, the second wind speed distribution estimating unit 14 calculates a distance L (m+1)-m' in the vertical direction from the two-dimensional plane in which the wind speed distribution u Rbm' (x, y) is estimated to the two-dimensional plane in which the wind speed distribution is u Rbm+1 (x, y).

[0124] The second wind speed distribution estimating unit 14 calculates a weighting coefficient w m for the wind speed distribution u Rbm (x, y) in the two-dimensional plane and a weighting coefficient w m+1 for the wind speed distribution u Rbm+1 (x, y) in the two-dimensional plane on the basis of the distance L m-m' and the distance L (m+1)-m' as expressed in the following Expressions (12) to (14). w m = L − L m − m ′ L w m + 1 = L − L m + 1 − m ′ L L m − m ′ + L m + 1 − m ′ = L

[0125] The second wind speed distribution estimating unit 14 estimates the wind speed distribution u Rbm' (x, y) in the two-dimensional plane between the two-dimensional plane in which the wind speed distribution is u Rbm (x, y) and the two-dimensional plane in which the wind speed distribution is u Rbm+1 (x, y) on the basis of the weighting coefficient w m and the weighting coefficient w m+1 as expressed in the following Expression (15). u Rbm ′ x y = w m × u Rbm x y + w m + 1 × u Rbm + 1 x y w m + w m + 1

[0126] Here, the second wind speed distribution estimating unit 14 estimates the wind speed distribution u Rbm' (x, y) in the two-dimensional plane as a weighted average based on the distance from the two-dimensional plane in which the wind speed distribution u Rbm' (x, y) is estimated. However, this is merely an example, and for example, the second wind speed distribution estimating unit 14 compares the distance L m-m' with the distance L (m+1)-m' . When the distance L m-m' is shorter than the distance L (m+1)-m' , the second wind speed distribution estimating unit 14 may use the wind speed distribution u Rbm (x, y) in the two-dimensional plane as the wind speed distribution u Rbm' (x, y) in the two-dimensional plane, and when the distance L m-m' is equal to or longer than the distance L (m+1)-m' , the second wind speed distribution estimating unit 14 may use the wind speed distribution u Rbm+1 (x, y) in the two-dimensional plane as the wind speed distribution u Rbm' (x, y) in the two-dimensional plane.

[0127] Further, here, the second wind speed distribution estimating unit 14 estimates the wind speed distribution u Rbm' (x, y) in one two-dimensional plane as the wind speed distribution in the two-dimensional plane between the two-dimensional plane in which the wind speed distribution is u Rbm (x, y) and the two-dimensional plane in which the wind speed distribution is u Rbm+1 (x, y). However, this is merely an example, and the second wind speed distribution estimating unit 14 may estimate wind speed distributions u Rbm' (x, y) of two or more two-dimensional planes having different positions in the vertical direction as the wind speed distribution in the two-dimensional plane between the two-dimensional plane in which the wind speed distribution is u Rbm (x, y) and the two-dimensional plane in which the wind speed distribution is u Rbm+1 (x, y).

[0128] When M is 1, the wind speed prediction unit 15 acquires the wind speed distribution u Rb1 (x, y) in the two-dimensional plane from the first wind speed distribution estimating unit 13.

[0129] When M is 2 or more, the wind speed prediction unit 15 acquires the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes from the first wind speed distribution estimating unit 13, and acquires wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional plane from the second wind speed distribution estimating unit 14.

[0130] When M is 1, the wind speed prediction unit 15 predicts the wind speed u(t) in the observation region in the two-dimensional plane from the wind speed distribution u Rb1 (x, y) in the two-dimensional plane using the two-dimensional N-S equation (step ST7 in FIG. 6).

[0131] When M is 2 or more, the wind speed prediction unit 15 selects the wind speed distribution in any two-dimensional plane from among the wind speed distributions u Rb1 (x, y) to u Rbm (x, y) in the plurality of two-dimensional planes and the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional planes.

[0132] When the observation point of the prediction target of the wind speed u(t) is present in the observation region in the two-dimensional plane corresponding to the range bin Rb m , the wind speed prediction unit 15 selects the wind speed distribution u Rbm (x, y). When the observation point of the prediction target of the wind speed u(t) is present in the observation region in the two-dimensional plane corresponding to the range bin Rb m' , the wind speed prediction unit 15 selects the wind speed distribution u Rbm' (x, y).

[0133] The wind speed prediction unit 15 predicts the wind speed u(t) in the observation region in the selected two-dimensional plane from the wind speed distribution in the selected two-dimensional plane using the two-dimensional N-S equation (step ST7 in FIG. 6).

[0134] Hereinafter, prediction processing of the wind speed u(t) performed by the wind speed prediction unit 15 will be specifically described.

[0135] Wind is generally considered an uncompressed fluid. The two-dimensional N-S equation regarding the uncompressed fluid is expressed by the following Expressions (16) and (17). ∂ u ∂ t = − u ⋅ ∇ u + ν Δ u − 1 ρ ∇ p + f ∇ ⋅ u = 0

[0136] In Expression (16), u is a wind speed vector, v is a viscosity coefficient, ρ is a density, p is a pressure, f is an external force vector, and t is a time variable.

[0137] Since the viscosity coefficient v is very small, the viscosity of the wind can be ignored. In addition, the term of the external force vector F can also be ignored. Therefore, under the condition that each of the viscosity of the wind and the external force is absent, Expression (16) is expressed as the following Expression (18), and Expression (17) is expressed as the following Expression (19). ∂ u ∂ t = − u ⋅ ∇ u − 1 ρ ∇ p ∇ ⋅ u = 0

[0138] The first term on the right side of Expression (18) is referred to as an advection term, and the second term on the right side is referred to as a pressure term.

[0139] When the wind speed vector u is a two-dimensional vector u = [u x , u y ] T< , Expression (18) is expressed as the following Expressions (20) and (21). ∂ u x ∂ t = − u x ∂ u x ∂ x − u y ∂ u x ∂ y − 1 ρ ∂ p ∂ x ∂ u y ∂ t = − u x ∂ u y ∂ x − u y ∂ u y ∂ y − 1 ρ ∂ p ∂ y ∂ u x ∂ x + ∂ u y ∂ y = 0

[0140] When Expression (22) is used, Expression (20) is expressed as the following Expression (23), and Expression (21) is expressed as the following Expression (24). ∂ u x ∂ t = ∂ u x 2 ∂ x − ∂ u x u y ∂ y − 1 ρ ∂ p ∂ x ∂ u y ∂ t = − ∂ u x u y ∂ x − ∂ u y 2 ∂ y − 1 ρ ∂ p ∂ y

[0141] When the wind speed vector u is discretized with respect to time by the difference method, Expression (23) is expressed as the following Expression (25), and Expression (24) is expressed as the following Expression (26). u x g + 1 − u x g dt = − ∂ u x 2 ∂ x − ∂ u x u y ∂ y − 1 ρ ∂ p ∂ x u y g + 1 − u y g dt = − ∂ u x u y ∂ x − ∂ u y 2 ∂ y − 1 ρ ∂ p ∂ y

[0142] In the Expressions (25) and (26), g and g + 1 are symbols indicating time steps.

[0143] By using an intermediate variable u x *, Expression (25) can be divided into two Expressions as the following Expressions (27) and (28).

[0144] In addition, by using an intermediate variable u y *, Expression (26) can be divided into two Expressions as the following Expressions (29) and (30). u x * − u x g dt = − ∂ u x 2 ∂ x − ∂ u x u y ∂ y u x g + 1 − u x * dt = − 1 ρ ∂ p ∂ x u y * − u y g dt = − ∂ u x u y ∂ x − ∂ u y 2 ∂ y u y g + 1 − u y * dt = − 1 ρ ∂ p ∂ y

[0145] Expression (27) is transformed as Expression (31), and Expression (29) is transformed as Expression (32). u x * = u x g − dt ∂ u x 2 ∂ x + ∂ u x u y ∂ y u y * = u y g − dt ∂ u x u y ∂ x + ∂ u y 2 ∂ y

[0146] Expression (28) is expressed as Expression (33) below when differentiated in the x direction, and Expression (30) is expressed as Expression (34) below when differentiated in the y direction. ∂ u x g + 1 ∂ x − ∂ u x * ∂ x dt = − 1 ρ ∂ 2 p ∂ x 2 ∂ u y g + 1 ∂ y − ∂ u y * ∂ y dt = − 1 ρ ∂ 2 p ∂ y 2

[0147] When the left sides of Expressions (33) and (34) are added to each other and the right sides of Expressions (33) and (34) are added to each other, the following Expression (35) is obtained. à 1 dt ∂ u x g + 1 ∂ x + ∂ u y g + 1 ∂ y + ∂ u x * ∂ x − ∂ u y * ∂ y = − 1 ρ ∂ 2 p ∂ x 2 + ∂ 2 p ∂ y 2

[0148] According to Expression (22), ∇·u g+1< = 0 needs to be satisfied for the estimated wind speed vector u at the next time step g + 1. Therefore, the following Expression (36) is obtained as the Poisson's equation regarding pressure from Expression (35). 1 ρ ∂ 2 p ∂ x 2 + ∂ 2 p ∂ y 2 = 1 dt ∂ u x * ∂ x + ∂ u y * ∂ y

[0149] The wind speed prediction unit 15 discretizes the wind speed vector u g< = [u x g< , u y g< ] T< at the current time step g in the spatial direction.

[0150] The wind speed prediction unit 15 calculates the intermediate variable u x * by substituting the discretized u x g< into Expression (31), and calculates the intermediate variable u y * by substituting the discretized u y g< into Expression (32).

[0151] The wind speed prediction unit 15 calculates the pressure p from Expression (36) using the intermediate variable u x * and the intermediate variable u y *.

[0152] The wind speed prediction unit 15 calculates a wind speed vector u x g+1< of the next time step g + 1 from Expression (28) using the intermediate variable u x * and the pressure p.

[0153] The wind speed prediction unit 15 calculates a wind speed vector u y g+1< of the next time step g + 1 from Expression (30) using the intermediate variable u y * and the pressure p.

[0154] The wind speed u(t) predicted by the wind speed prediction unit 15 is [u x g+1< , u y g+1< ] T< , and the wind speed u(t) is displayed on, for example, a display (not illustrated).

[0155] A conventional wind speed prediction device performs a simulation using a weather model in order to predict a wind speed in an observation region, and the simulation using the weather model is three-dimensional calculation processing. For example, when the number of observation points in the three-dimensional observation region is Cx × Cy × Cz, and C = Cx = Cy = Cz, the calculation amount order is approximately proportional to C 3< .

[0156] On the other hand, the calculation processing of the wind speed prediction device 3 illustrated in FIG. 3 is two-dimensional calculation processing. For example, when the number of two-dimensional planes is M, the number of observation points in the observation region present in the two-dimensional plane is Cx × Cy, and C = Cx = Cy, the calculation amount order is approximately proportional to C 2< × M. When M = C, the calculation amount order of the wind speed prediction device 3 illustrated in FIG. 3 is substantially the same as the calculation amount order of the conventional wind speed prediction device.

[0157] However, it is assumed that three-dimensional calculation processing uses a three-dimensional N-S equation that does not ignore vertical wind. In the three-dimensional N-S equation, it is difficult to greatly expand the interval between the observation points in the vertical direction, and the maximum value of the interval between the observation points at which the three-dimensional N-S equation can be solved is determined as a Courant condition in the vertical direction.

[0158] On the other hand, in the two-dimensional N-S equation ignoring the vertical wind, which is the two-dimensional calculation processing, since there is no Courant condition in the vertical direction, it is possible to satisfy Cz >> M. For example, when (M / Cz) = 1 / Q, the calculation amount order of the wind speed prediction device 3 illustrated in FIG. 3 is approximately 1 / Q of the calculation amount order of the conventional wind speed prediction device. Q is a value of about 10 to 1000.

[0159] As described above, for example, when the terrain of the ground surface vertically below the observation point is not a special terrain but a general terrain, there is no practical problem in observing the wind affecting the flying object as long as the wind speed prediction device 3 can observe the wind speed in the direction parallel to the two-dimensional plane even if the wind speed prediction device 3 cannot observe the wind speed in the vertical direction with respect to the two-dimensional plane where the observation point is present.

[0160] Even if the terrain of the ground surface vertically below the observation point is a special terrain, the wind speed in the direction parallel to the two-dimensional plane can be used as useful information for observing the wind affecting the flying object.

[0161] In a first embodiment described above, the wind speed prediction device 3 includes the scattering signal acquiring unit 11 to acquire a scattering signal that is each of a plurality of beams after being emitted to space and scattered in the space, the beams having mutually different elevation angles, which are angles formed by a line-of-sight direction and a horizontal direction, the Doppler frequency calculating unit 12 to set a range bin width as distance resolution of each of scattering signals acquired by the scattering signal acquiring unit 11 according to an elevation angle of each of the beams emitted to space, and calculate a Doppler frequency of a range bin corresponding to a two-dimensional plane including an observation region from the each of the scattering signals, the first wind speed distribution estimating unit 13 to estimate a wind speed distribution in the two-dimensional plane from a plurality of Doppler frequencies calculated by the Doppler frequency calculating unit 12 using a VVP method, and the wind speed prediction unit 15 to predict a wind speed in the observation region from the wind speed distribution in the two-dimensional plane estimated by the first wind speed distribution estimating unit 13 using a two-dimensional Navier-Stokes equation. Therefore, the wind speed prediction device 3 can predict the wind speed without performing three-dimensional calculation processing.

[0162] In the wind speed prediction device 3 illustrated in FIG. 3, the first wind speed distribution estimating unit 13 estimates the wind speed distribution u Rbm (x, y) in the two-dimensional plane corresponding to each range bin Rb m . However, the wind speed distribution u Rbm (x, y) estimated by the first wind speed distribution estimating unit 13 may slightly deviate from the actual wind speed distribution due to the influence of observation noise. The first wind speed distribution estimating unit 13 may smooth the wind speed distribution u Rbm (x, y) in the time direction in order to suppress the deviation due to the influence of the observation noise.

[0163] The smoothing of the wind speed distribution u Rbm (x, y) in the time direction can be implemented by fusing a fluid movement model such as a two-dimensional N-S equation and the estimated wind speed distribution u Rbm (x, y). The fusion of the fluid movement model and the wind speed distribution u Rbm (x, y) is called data assimilation.

[0164] For example, it is assumed that the true wind speed distribution u Rbm (x, y) at the previous time step g-1 is represented by u t=tg-1 , and the true wind speed distribution u Rbm (x, y) at the current time step g is represented by u t=tg . The two-dimensional N-S equation is a fluid model that can nonlinearly describe the relationship between u t=tg-1 and u t=tg , and if a nonlinear function is represented by f(·), the following Expression (37) holds. u t = t n = f u t = t n − 1

[0165] The data assimilation processing is processing performed using a data assimilation function As(·) as expressed in the following Expression (38). That is, the data assimilation processing calculates a smoothed u(x, y; t=tg ) bar of the current time step g by adding u(x, y; t=tg ) which is the wind speed distribution u Rbm (x, y) of the current time step g and f(u(x, y; t=tg-1 ) bar) obtained by time-developing the wind speed distribution u Rbm (x, y) of the previous time step g-1 to the current time step g. In the text of the specification, since it is not possible to attach a symbol "-" above a character due to the electronic application, for example, it is written as a u(x, y; t=tg ) bar. u x , y ; t = t g ¯ = As u x , y ; t = t g , f u x , y ; t = t g − 1 ¯

[0166] Here, the first wind speed distribution estimating unit 13 performs the data assimilation processing using the data assimilation function As(·). However, this is merely an example, and the first wind speed distribution estimating unit 13 may perform the data assimilation processing by, for example, a particle filter using a two-dimensional N-S equation or an ensemble Kalman filter using the two-dimensional N-S equation. In addition, the first wind speed distribution estimating unit 13 may perform the data assimilation processing by an extended Kalman filter in which the two-dimensional N-S equation is approximated to a linear function or a variational method using the two-dimensional N-S equation.Second Embodiment

[0167] In the wind speed prediction device 3 according to the first embodiment, the wind speed prediction unit 15 selects the wind speed distribution in any two-dimensional plane from among the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes and the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional planes.

[0168] In a second embodiment, a wind speed prediction device 3 including a wind speed prediction unit 16 will be described, the wind speed prediction unit 16 estimating a wind speed distribution in a two-dimensional plane including a flight area in which a flying object is scheduled to fly as an observation region from the wind speed distributions u Rb1 (x, y) to u Rbm (x, y) in the plurality of two-dimensional planes and the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional planes.

[0169] FIG. 9 is a configuration diagram illustrating the wind speed prediction device 3 according to the second embodiment. In FIG. 9, the same reference numerals as those in FIG. 3 denote the same or corresponding parts, and thus description thereof is omitted.

[0170] FIG. 10 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the second embodiment. In FIG. 10, the same reference numerals as those in FIG. 4 denote the same or corresponding parts, and thus description thereof is omitted.

[0171] The wind speed prediction device 3 illustrated in FIG. 9 includes a scattering signal acquiring unit 11, a Doppler frequency calculating unit 12, a first wind speed distribution estimating unit 13, a second wind speed distribution estimating unit 14, and the wind speed prediction unit 16.

[0172] The wind speed prediction unit 16 is implemented by, for example, a wind speed prediction circuit 26 illustrated in FIG. 10.

[0173] The wind speed prediction unit 16 estimates a wind speed distribution in a two-dimensional plane including a flight area as an observation region from the wind speed distributions u Rb1 (x, y) to u RbM (x, y) in the plurality of two-dimensional planes estimated by the first wind speed distribution estimating unit 13 and the wind speed distributions u Rb1' (x, y) to u RbM-1' (x, y) in the two-dimensional plane estimated by the second wind speed distribution estimating unit 14.

[0174] The wind speed prediction unit 16 predicts the wind speed u(t) in the flight area from the estimated wind speed distribution in the two-dimensional plane using the two-dimensional N-S equation.

[0175] In FIG. 9, it is assumed that each of the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, and the wind speed prediction unit 16, which are components of the wind speed prediction device 3, is implemented by dedicated hardware as illustrated in FIG. 10. That is, it is assumed that the wind speed prediction device 3 is implemented by the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, and the wind speed prediction circuit 26.

[0176] Each of the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, and the wind speed prediction circuit 26 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, ASIC, FPGA, or a combination thereof.

[0177] The components of the wind speed prediction device 3 are not limited to those implemented by dedicated hardware, and the wind speed prediction device 3 may be implemented by software, firmware, or a combination of software and firmware.

[0178] In a case where the wind speed prediction device 3 is implemented by software, firmware, or the like, a program for causing a computer to execute each processing procedure in the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, and the wind speed prediction unit 16 is stored in the memory 31 illustrated in FIG. 5. Then, the processor 32 illustrated in FIG. 5 executes the program stored in the memory 31.

[0179] Further, FIG. 10 illustrates an example in which each of the components of the wind speed prediction device 3 is implemented by dedicated hardware, and FIG. 5 illustrates an example in which the wind speed prediction device 3 is implemented by software, firmware, or the like. However, this is merely an example, and some components in the wind speed prediction device 3 may be implemented by dedicated hardware, and the remaining components may be implemented by software, firmware, or the like.

[0180] Next, the operation of the wind speed prediction device 3 illustrated in FIG. 9 will be described. The wind speed prediction device 3 is similar to the wind speed prediction device 3 illustrated in FIG. 3 except for the wind speed prediction unit 16, and thus only the operation of the wind speed prediction unit 16 will be described here.

[0181] Here, for convenience of description, it is assumed that the flight area is present in a two-dimensional plane between a two-dimensional plane in which the wind speed distribution is u Rbm (x, y) and a two-dimensional plane in which the wind speed distribution is u Rbm' (x, y).

[0182] The wind speed prediction unit 16 calculates a distance L p-m in the vertical direction from the two-dimensional plane in which the flight area is present to the two-dimensional plane in which the wind speed distribution is u Rbm (x, y).

[0183] Further, the wind speed prediction unit 16 calculates a distance L p-m' in the vertical direction from the two-dimensional plane in which the flight area is present to the two-dimensional plane in which the wind speed distribution is u Rbm' (x, y).

[0184] The wind speed prediction unit 16 calculates a weighting coefficient w p-m for the wind speed distribution u Rbm (x, y) in the two-dimensional plane and a weighting coefficient w p-m' for the wind speed distribution u Rbm' (x, y) in the two-dimensional plane on the basis of the distance L p-m and the distance L p-m' as expressed in the following Expressions (39) to (41). w p − m = L ′ − L p − m L ′ w p − m ′ = L ′ − L p − m ′ L ′ L p − m + L p − m ′ = L ′

[0185] The wind speed prediction unit 16 estimates the wind speed distribution u Rbp (x, y) of the two-dimensional plane where the flight area is present on the basis of the weighting coefficient w p-m and the weighting coefficient w p-m' as expressed in the following Expression (42). u Rbp x y = w p − m × u Rbm x y + w p − m ′ × u Rbm ′ x y w p − m + w p − m ′

[0186] The wind speed prediction unit 16 predicts the wind speed u(t) in the flight area from the estimated wind speed distribution u Rbp (x, y) in the two-dimensional plane using the two-dimensional N-S equation. The prediction processing of the wind speed u(t) performed by the wind speed prediction unit 16 is similar to the prediction processing of the wind speed u(t) performed by the wind speed prediction unit 15 illustrated in FIG. 3.Third Embodiment

[0187] In a third embodiment, a wind speed prediction device 3 including a time notification unit 17 that calculates a time t at which the wind speed u(t) predicted by the wind speed prediction unit 15 becomes equal to or higher than the threshold Th and gives a notification of the time t will be described.

[0188] FIG. 11 is a configuration diagram illustrating the wind speed prediction device 3 according to the third embodiment. In FIG. 11, the same reference numerals as those in FIGS. 3 and 9 denote the same or corresponding parts, and thus description thereof is omitted.

[0189] FIG. 12 is a hardware configuration diagram illustrating hardware of the wind speed prediction device 3 according to the third embodiment. In FIG. 12, the same reference numerals as those in FIGS. 4 and 10 denote the same or corresponding parts, and thus description thereof is omitted.

[0190] The wind speed prediction device 3 illustrated in FIG. 11 includes a scattering signal acquiring unit 11, a Doppler frequency calculating unit 12, a first wind speed distribution estimating unit 13, a second wind speed distribution estimating unit 14, the wind speed prediction unit 15, and the time notification unit 17.

[0191] The time notification unit 17 is implemented by, for example, a time notification circuit 27 illustrated in FIG. 12.

[0192] The time notification unit 17 specifies the time t at which the wind speed u(t) predicted by the wind speed prediction unit 15 becomes equal to or higher than the threshold Th.

[0193] The time notification unit 17 gives a notification of the specified time t to, for example, a flying object or a control tower monitoring a flight state of the flying object.

[0194] FIG. 11 illustrates an example in which the time notification unit 17 is applied to the wind speed prediction device 3 illustrated in FIG. 3. However, this is merely an example, and the time notification unit 17 may be applied to the wind speed prediction device 3 illustrated in FIG. 9.

[0195] In FIG. 11, it is assumed that each of the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, the wind speed prediction unit 15, and the time notification unit 17, which are components of the wind speed prediction device 3, is implemented by dedicated hardware as illustrated in FIG. 12. That is, it is assumed that the wind speed prediction device 3 is implemented by the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, the wind speed prediction circuit 25, and the time notification circuit 27.

[0196] Each of the scattering signal acquiring circuit 21, the Doppler frequency calculating circuit 22, the first wind speed distribution estimating circuit 23, the second wind speed distribution estimating circuit 24, the wind speed prediction circuit 25, and the time notification circuit 27 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, ASIC, FPGA, or a combination thereof.

[0197] The components of the wind speed prediction device 3 are not limited to those implemented by dedicated hardware, and the wind speed prediction device 3 may be implemented by software, firmware, or a combination of software and firmware.

[0198] In a case where the wind speed prediction device 3 is implemented by software, firmware, or the like, a program for causing a computer to execute each processing procedure in the scattering signal acquiring unit 11, the Doppler frequency calculating unit 12, the first wind speed distribution estimating unit 13, the second wind speed distribution estimating unit 14, the wind speed prediction unit 15, and the time notification unit 17 is stored in the memory 31 illustrated in FIG. 5. Then, the processor 32 illustrated in FIG. 5 executes the program stored in the memory 31.

[0199] Further, FIG. 12 illustrates an example in which each of the components of the wind speed prediction device 3 is implemented by dedicated hardware, and FIG. 5 illustrates an example in which the wind speed prediction device 3 is implemented by software, firmware, or the like. However, this is merely an example, and some components in the wind speed prediction device 3 may be implemented by dedicated hardware, and the remaining components may be implemented by software, firmware, or the like.

[0200] Next, the operation of the wind speed prediction device 3 illustrated in FIG. 11 will be described. The wind speed prediction device 3 is similar to the wind speed prediction device 3 illustrated in FIG. 3 except for the time notification unit 17, and thus only the operation of the time notification unit 17 will be described here.

[0201] The time notification unit 17 acquires the predicted wind speed u(t) from the wind speed prediction unit 15 every time the wind speed prediction unit 15 predicts the wind speed u(t).

[0202] Every time the predicted wind speed u(t) is acquired, the time notification unit 17 compares the wind speed u(t) with the threshold Th. The threshold Th is, for example, a lower limit value of the wind speed that may affect the flight of the flying object. The threshold Th may be stored in the internal memory of the time notification unit 17 or may be given from the outside.

[0203] The time notification unit 17 specifies the time t at which the wind speed u(t) is equal to or higher than the threshold Th on the basis of the comparison result between the wind speed u(t) and the threshold Th.

[0204] The time notification unit 17 gives a notification of the specified time t to, for example, a flying object or a control tower monitoring a flight state of the flying object.

[0205] In the third embodiment described above, the wind speed prediction device 3 illustrated in FIG. 11 includes the time notification unit 17 that specifies the time at which the wind speed predicted by the wind speed prediction unit 15 becomes equal to or higher than the threshold and gives a notification of the specified time t. Therefore, the wind speed prediction device 3 illustrated in FIG. 11 can predict the wind speed without performing three-dimensional calculation processing, as does the wind speed prediction device 3 illustrated in FIG. 3. In addition, it is possible to notify the time when wind that may affect the flight of the flying object blows.INDUSTRIAL APPLICABILITY

[0206] The present disclosure is suitable for a wind speed prediction device, a wind speed prediction method, and a radar device.REFERENCE SIGNS LIST

[0207] 1: radar device, 2: beam transmitting and receiving unit, 3: wind speed prediction device, 11: scattering signal acquiring unit, 12: Doppler frequency calculating unit, 13: first wind speed distribution estimating unit, 14: second wind speed distribution estimating unit, 15, 16: wind speed prediction unit, 17: time notification unit, 21: scattering signal acquiring circuit, 22: doppler frequency calculating circuit, 23: first wind speed distribution estimating circuit, 24: second wind speed distribution estimating circuit, 25, 26: wind speed prediction circuit, 27: time notification circuit, 31: memory, 32: processor

Claims

1. A wind speed prediction device (3), comprising: a scattering signal acquiring unit (11) to acquire a scattering signal that is each of a plurality of beams after being emitted to space and scattered in the space, the beams having mutually different elevation angles, which are angles formed by a line-of-sight direction and a horizontal direction; a Doppler frequency calculating unit (12) to set a range bin width as distance resolution of each of scattering signals acquired by the scattering signal acquiring unit (11) according to an elevation angle of each of the beams emitted to the space to enable monitoring of observation points of same range bins, the range bins being at the same altitude even for the plurality of beams at different elevation angles, and calculate a Doppler frequency of a range bin corresponding to a two-dimensional plane including an observation region from the each of the scattering signals; a first wind speed distribution estimating unit (13) to estimate a wind speed distribution in the two-dimensional plane from the plurality of Doppler frequencies calculated by the Doppler frequency calculating unit (12); and a wind speed prediction unit (15, 16) to predict a wind speed in the observation region from the wind speed distribution in the two-dimensional plane estimated by the first wind speed distribution estimating unit (13); characterised in that; the first wind speed distribution estimating unit (13) is configured to estimate the wind speed distribution in the two-dimensional plane using a volume velocity processing (VVP) method; and the wind speed predication unit (15, 16) is configured to predict the wind speed in the observation region using a two-dimensional Navier-Stokes equation.

2. The wind speed prediction device according to claim 1, wherein a plurality of range bins is provided, and the Doppler frequency calculating unit (12) is configured to calculate each of the Doppler frequencies in the plurality of range bins from the each of the scattering signals to express a three-dimensional region to be observed as a stack of a plurality of two-dimensional planes, and the first wind speed distribution estimating unit (13) is configured to estimate a wind speed distribution in each of the two-dimensional planes having a different altitude corresponding to each of the range bins corresponding to the two-dimensional planes having different altitudes using the VVP method from the plurality of the Doppler frequencies in each of the range bins calculated by the Doppler frequency calculating unit (12).

3. The wind speed prediction device according to claim 2, further comprising: a second wind speed distribution estimating unit (14) to estimate a wind speed distribution in a two-dimensional plane between two adjacent two-dimensional planes, selected from the plurality of two-dimensional planes, from the wind speed distributions in said two-dimensional planes having different altitudes estimated by the first wind speed distribution estimating unit (13), and to estimate a two-dimensional wind speed distribution having a designated altitude by linearly interpolating the wind speed distribution in the two-dimensional plane in an altitude direction, wherein the wind speed prediction unit (15,16) is configured to select the wind speed distribution in the two-dimensional plane including a flight area in which a flying object is scheduled to fly from among the wind speed distributions in the plurality of two-dimensional planes having different altitudes estimated by the first wind speed distribution estimating unit (13) and the second wind speed distribution estimating unit (14), and predict the wind speed in the flight area at which the flying object is scheduled to fly, from the selected wind speed distribution using a two-dimensional Navier-Stokes equation.

4. The wind speed prediction device according to claim 2, further comprising: a second wind speed distribution estimating unit (14) to estimate the wind speed distribution in a two-dimensional plane between two adjacent two-dimensional planes, selected from the plurality of two-dimensional planes, from the wind speed distributions in said two adjacent two-dimensional planes having different altitudes estimated by the first wind speed distribution estimating unit (13), and estimate a two-dimensional wind speed distribution having a designated altitude by linearly interpolating the wind speed distribution in the two-dimensional plane in an altitude direction, wherein the wind speed prediction unit (16) is configured to estimate the wind speed distribution in the two-dimensional plane at a flight altitude at which a flying object is scheduled to fly using the second wind speed distribution estimating unit, and predict a wind speed in the flight area at which the flying object is scheduled to fly, from the estimated wind speed distribution using a two-dimensional Navier-Stokes equation.

5. The wind speed prediction device according to claim 1, further comprising a time notification unit (17) to specify a time at which the wind speed predicted by the wind speed prediction unit (15, 16) is equal to or higher than a threshold and gives a notification of the time.

6. A wind speed prediction method comprising, when each of a plurality of beams having mutually different elevation angles, which are angles formed by a line-of-sight direction and a horizontal direction, is emitted to space: acquiring, by a scattering signal acquiring unit (11), a scattering signal that is each of the beams after being scattered in the space; setting, by a Doppler frequency calculating unit (12), a range bin width as distance resolution of each of scattering signals acquired by the scattering signal acquiring unit (11) according to an elevation angle of each of the beams emitted to the space to enable monitoring of observation points of same range bins, the range bins being at the same altitude even for the plurality of beams at different elevation angles; calculating, by the Doppler frequency calculating unit (12), a Doppler frequency of a range bin corresponding to a two-dimensional plane including an observation region from the each of the scattering signals; estimating, by a first wind speed distribution estimating unit (13), a wind speed distribution in the two-dimensional plane from a plurality of Doppler frequencies calculated by the Doppler frequency calculating unit (12); and predicting, by a wind speed prediction unit (15, 16), a wind speed in the observation region from the wind speed distribution in the two-dimensional plane estimated by the first wind speed distribution estimating unit (13); characterised in that; the first wind speed distribution estimating unit (13) estimates the wind speed distribution in the two-dimensional plane using a volume velocity processing (VVP) method; and the wind speed predication unit (15, 16) predicts the wind speed in the observation region using a two-dimensional Navier-Stokes equation.