KR101892502B1 - calculation system for storm relative helicity using wind profiler and calculation method for storm relative helicity using the same - Google Patents

Abstract

The present invention relates to a Storm Relative Helicity (SRH) calculation system and an SRH calculation method using the same. More specifically, the present invention relates to the system and the method for calculating SRH using a wind profiler. The SRH calculation system using a wind profiler comprises: a wind vector calculation unit for removing non-Gaussian propagation signals from a source of the wind profiler and reassigning wind vectors; and an SRH calculation unit for calculating the SRH based on the wind vector calculated in the wind vector calculation unit. Accordingly, the present invention can provide the system and the method which can easily and accurately calculate the SRH.

Description

[0001] The present invention relates to an SRH calculation system using a wind profiler and a SRH calculation method using the same,

The present invention relates to an SRH calculating system and a SRH calculating method using the SRH calculating system, and more particularly, to a system and a method for calculating an SRH using a wind profiler.

)를 가지는 구름계를 의미한다. MCS에서 관측되는 가장 일반적인 현상은 강한 대류성 에코들이 선형으로 형성되고, 진행방향을 기준으로 그 후면에 층운형 에코가 따르는 스콜선(Squall Line)이다. 일반적으로 레이더 에코가 선형으로 길게 늘어선 현상이다. MCS는 매우 복잡한 열역학적인 과정을 통해서 형성되기 때문에 모델에 의한 예측이 어렵지만, 기상레이더, 윈드프로파일러, 라디오미터와 같은 원격탐사 장비를 이용한 관측으로 구조를 파악할 수 있다. 장마전선에 의해 자주 발생하는 스콜선은 연직 바람 쉬어를 이용하여 진단이 가능하다. 연직 바람 쉬어는 스콜선이 얼마나 강하고 오랫동안 유지될 것인지를 결정한다. 일반적으로 한반도에서 장마기간 동안 발생하는 집중호우는 역학적 불안정에 의해서 발생한다. 따라서, 쉽고 정확하게 SRH(Storm Relative Helicity)를 산출할 수 있는 방법이 필요하다.Muscular faults occurring on the rainy season are typically heavy rains caused by mesoscale convective systems (MCS). The MCS is associated with a combination of thunderstorms and areas of continuous precipitation, with a horizontal scale of at least 100 km in one direction (Meso-

). The most common phenomenon observed in MCS is a squall line in which strong convective echoes are formed linearly and followed by a stratified echo on the back of the direction of travel. Generally, the radar echo is linearly long. MCS is difficult to predict by model because it is formed through a very complex thermodynamic process, but it can be grasped by observations using remote sensing equipment such as weather radar, wind profilers, and radio meters. Scorelines, which are frequently generated by the power line, can be diagnosed using a vertical wind shear. A vertical wind shear determines how strong and long the squad line will remain. Generally, heavy rainfall during the rainy season on the peninsula is caused by mechanical instability. Therefore, there is a need for a method that can easily and accurately calculate SRH (Storm Relative Helicity).

Korean Patent Laid-Open Publication No. 2017-0002741 (published on April 1, 2017)

It is an object of the present invention to provide a system and method that can easily and accurately calculate SRH.

In order to achieve the above-mentioned object, the present invention provides a wind turbine generator comprising: a wind vector calculating unit for removing a non-vapor wave signal from raw data of a wind profiler and outputting a wind vector; and a wind vector calculating unit And an SRH calculating unit for calculating SRH (Storm Relative Helicity).

Wherein the wind vector calculator comprises: a raw data loading module for loading raw data of the wind profiler; a noise power level determining module for determining a noise power level of the Doppler spectrum of the raw data; A non-meteorological echo cancellation module that removes the radio signal from the altitude-dependent Doppler spectrum, a standby to determine the peak due to the atmospheric signal in the highly doped Doppler spectrum from which the non-terrestrial propagation signal exceeds the noise power level in the non- A motion peak signal selection module, a spectrum moment calculation module for determining a minimum point and a maximum point of a fast Fourier transform analysis point of the Doppler spectrum around the peak, and calculating a spectrum moment through the determined Doppler spectrum parameter, And in the spectrum moment calculation module A wind vector calculation module for calculating the wind vector of the wind profiler in the east-west direction, the north-south direction and the vertical direction through the first-order spectrum moments of the four slope beams and the vertical beams with respect to all the high- .

The SRH calculator includes an altitude average wind vector calculating module for calculating an average wind vector up to a predetermined altitude based on the wind vector calculated by the wind vector calculator, A vertical sheer calculation module for calculating the vertical sheath of the wind vector by using the upper-level elevation wind vector immediately for each elevation, a vector deviation calculation module for calculating the average wind vector after calculating the vertical sheath, And an accumulation module for calculating the SRH by accumulating the difference between the average wind vector at the altitude and the stom motion vector and the product of the vertical sheath of the wind vector up to a predetermined altitude.

The present invention further includes a step of removing the non-Gaussian propagation signal from the source data of the wind profiler and restoring the wind vector, and a step of generating the SRH output unit SRH based on the wind vector calculated by the wind vector calculation unit And calculating a Storm Relative Helicity of the SRH.

Wherein the step of removing the non-Gaussian propagation signal from the wind data of the wind vector profiler and restoring the wind vector comprises the steps of: loading the wind data of the wind profiler with the source data loading module; Determining a noise power level of the Doppler spectrum by a noise power level determination module; removing the non-Gaussian echo cancellation module in the Doppler spectrum of the altitude by the non-Gaussian propagation signal exceeding the noise power level; Wherein the atmospheric motion peak signal selection module determines the peak due to the atmospheric signal in a highly Doppler spectrum in which the non-Gaussian propagation signal exceeding the noise power level is removed from the noise power level, and a fast Fourier transform of the Doppler spectrum Transform < / RTI > analysis points, Calculating a spectrum moment by a spectrum moment calculating module through a spectrum parameter and calculating a spectrum moment by a spectrum analyzer based on the wind profiler through four slope beams and a first order spectral moment of a vertical beam for all the high- The wind vector of the high altitude wind vector is calculated in the east-west direction, the north-south direction and the vertical direction.

The step of calculating SRH (Storm Relative Helicity) based on the wind vector calculated by the wind vector calculating unit may include calculating an SRH (Storm Relative Helicity) based on the wind vector calculated by the wind vector calculating unit, Estimating a storm motion vector by a storm motion vector estimation module based on the altitude average wind vector, calculating a steep motion vector using a wind vector at an upper altitude level, Calculating a vertical shear of the wind vector, calculating a vertical shear, and a vector deviation calculation module calculating an average wind vector, And the cumulative module calculates the SRH by accumulating the product of the vertical sheer of the SRH to a predetermined altitude.

이며, 상기

는 수평속도벡터(Horizontal velocity vector)이고, 상기

는 스톰의 운동벡터(Storm motion vector)이다.In the SRH,

, And

Is a horizontal velocity vector,

Is a Storm motion vector of Storm.

The present invention provides a system and method for easily and accurately calculating SRH.

1 is a block diagram of an SRH calculation system using a wind profiler according to the present invention;
2 is a map showing an installation point of a wind profiler in an SRH calculation system using a wind profiler according to the present invention;
3 is a result of calculating the SRH of the wind profile using the SRH calculation system using the wind profiler according to the present invention.
4 is a flowchart of a method of calculating SRH using a wind profiler according to the present invention.

Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.

It will be apparent to those skilled in the art that the present invention may be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, It is provided to let you know. Like reference numerals refer to like elements throughout.

1 is a block diagram of an SRH calculation system using a wind profiler according to the present invention.

As shown in FIG. 1, the SRH calculating system using the wind profiler according to the present invention comprises a wind vector calculating unit 100 for calculating a wind vector, a wind vector calculating unit 100 for calculating a wind vector based on the wind vector calculated by the wind vector calculating unit 100 And an SRH calculation unit 200 for calculating SRH.

The wind vector calculator 100 calculates a wind vector using a wind profiler. Wind Profiler is a device used to observe the vertical wind vector of the atmosphere and is well suited for diagnosing the mechanical instability of the phenomenon of the wind. The Wind Profiler uses the Doppler beam method to determine the vertical profile of the wind vector and transmits the electromagnetic pulse wave according to the fixed azimuth and elevation angles. Due to the characteristics of the atmospheric scattering process, the received signal has a phase difference due to the Doppler delay and determines the atmospheric wind vector component through the Doppler spectrum calculated by Fast Fourier Transform (FFT). Wind profiler is generally used for meteorological studies. It can provide vertical structure of wind with time resolution of water and spatial resolution of several tens of meters. Therefore, it is necessary to analyze the risk weather or to improve prediction accuracy of numerical forecasting model It is widely used as basic data. However, since the wind profiler estimates the wind vector through the signal processing of the radar wave, an error due to various non-vapor propagation signals occurs. Therefore, the present invention treats the raw data of the wind profiler before calculating the SRH to remove the non-Gaussian propagation signal and recalculate the wind vector. For this purpose, the wind vector calculator 100 includes a raw data loading module 110 for loading raw data of the wind profiler, a noise power level determining module 120, a non-geophone echo cancellation module 130, A signal selection module 140, a spectrum moment calculation module 150, and an altitude dependent wind vector calculation module 160. [

The noise power level determination module 120 determines the noise power level of the Doppler spectrum of the wind profiler raw data. This determines the noise power level assuming the most widely used white noise in recent years. It also sets the level of spectral density lower than the spectral peak and collects the spectral data below the predetermined level to construct a new spectrum. This process is repeated until the newly constructed spectrum satisfies the white noise condition.

The non-geophone echo cancellation module 130 determines the noise power level in the noise power level determination module 120 and then removes the non-Gaussian propagation signal in the Doppler spectrum according to the altitude.

The atmospheric motion peak signal selection module 140 determines the peak due to the atmospheric signal in the Doppler spectrum at the altitude where the non-vapor phase signal is removed from the non-phase echo cancellation module 130.

The spectrum moment calculation module 150 determines a minimum point and a maximum point of the FFT analysis point of the Doppler spectrum around the peak due to the waiting signal determined by the waiting motion peak signal selection module 140. [ Further, the spectrum moment is calculated through the determined Doppler spectrum parameter.

,

,

)를 산출한다.The altitude-dependent wind vector calculation module 160 calculates the wind vector of the wind profiler (i.e., the altitude-based wind vector) through the first-order spectral moment of the inclined beam and the vertical beam for all altitude sampling points calculated by the spectrum-

,

,

).

The SRH calculating unit 200 calculates the SRH based on the wind vector calculated by the wind vector calculating unit 100. For this purpose, the SRH calculator 200 includes an altitude-dependent average wind vector calculation module 210, a storm motion vector estimation module 200, a vertical sheath calculation module 230, a vector deviation calculation module 240, 250).

The altitude-based average wind vector calculation module 210 calculates an average wind vector up to a predetermined altitude. Here, the predetermined altitude is 0 to 3 km.

The storm motion vector estimation module 200 estimates a storm motion vector based on the altitude-based average wind vector calculated by the altitude-based average wind vector calculation module 210. Here, the storm motion vector is 75% of the average wind vector per altitude, and the direction is +30 degrees.

The vertical shear output module calculates the vertical sheath of the wind vector using the upper-level elevation wind vector at each altitude.

The vector deviation calculation module 240 calculates a vertical sheer in the vertical sheer calculation module and then calculates an average wind vector.

The accumulation module 250 accumulates the product of the difference between the average wind vector and the storm motion vector for each altitude and the vertical sheath of the wind vector up to a predetermined altitude. Here, the predetermined altitude is 0 to 3 km as illustrated in the altitude-based average wind vector calculation module 210, and can be calculated as shown in Equation 1 below. That is, as shown in Equation (1) below, SRH is a value obtained by integrating the change of the wind turbine from the ground to the height of 3 km along the direction of the storm, which means dynamic instability.

는 수평속도벡터(Horizontal velocity vector),

는 스톰의 운동벡터(Storm motion vector)이다.In Equation (1)

A horizontal velocity vector,

Is a Storm motion vector of Storm.

2 is a map showing an installation point of a wind profiler in an SRH calculation system using a wind profiler according to the present invention. 3 is a result of calculating the SRH of the wind profile using the SRH calculation system using the wind profiler according to the present invention.

Referring to FIGS. 2 and 3, the SRH calculation system using the wind profiler according to the present invention is applied, and the SRH of the rainfall occurrence period is calculated as shown in FIG. 3 in the data observed at the wind profiler point in FIG. Here, as shown in FIG. 2, it can be seen that the value of SRH is high due to the mechanical instability caused by the storm.

Next, an SRH calculation method using a wind profiler according to the present invention will be described with reference to the drawings. The description of the SRH calculating system using the wind profiler according to the present invention, which will be described later, will be omitted or briefly explained.

4 is a flowchart of a method for calculating an SRH using a wind profiler according to the present invention.

As shown in FIG. 4, the SRH calculation method using a wind profiler according to the present invention includes a step S1 of calculating a wind vector and a step S2 of calculating an SRH based on the calculated wind vector .

The step S1 of calculating the wind vector calculates the wind vector using the wind vector calculator. To this end, the step S1 of calculating the wind vector includes a step S1-1 of loading the raw data, a step S1-2 of determining the noise power level, a step S1-3 A step S1-4 of selecting an atmospheric motion peak signal, a step S1-5 of calculating a spectrum moment, and a step S1-6 of calculating an altitude wind vector.

In the step (S1-1) of loading the raw data, the raw data loading module loads the raw data of the wind profiler.

Step S1-2 of determining the noise power level, in the raw data of the Wind Profiler loaded in step S1-1 of loading the raw data, determines the noise power level of the Doppler spectrum by the noise power level determination module .

The step of removing the non-vapor echo (S1-3) comprises the steps of (S1-2) determining the noise power level, and then the noise power level is determined, and then the non-vapor echo canceling module removes the non-vapor signal from the Doppler spectrum do. Here, the non-gaseous wave propagation signal is removed when it is higher than the determined noise power level.

In the step S1-4 of selecting the atmospheric motion peak signal, the atmospheric motion peak signal selection module determines the peak due to the atmospheric signal in the Doppler spectrum in which the non-vapor echo is removed.

The step S1-5 of calculating the spectrum moment includes a step of selecting a peak of the atmospheric motion peak signal by a spectral moment calculation module based on a peak due to the waiting signal determined in the step S1-4, , FFT) Determine the minimum and maximum points of the analysis point. Further, the spectrum moment is calculated through the determined Doppler spectrum parameter.

,

,

)를 산출한다.The step S1-6 of calculating the altitude-dependent wind vector is carried out by calculating the spectral moments through the first-order spectral moments of the four oblique beams and the vertical beams for all the high-order sampled points calculated in the step S1-5. The altitude-dependent wind vector output module calculates the wind vector of the wind profiler (

,

,

).

The step S2 of calculating the SRH calculates the SRH calculating part SRH based on the wind vector calculated in the step S1 of calculating the wind vector. The step S2 of calculating the SRH includes a step S2-1 of calculating an average wind vector per altitude, a step S2-2 of estimating a stomata motion vector, a step S2-3 calculating a vertical sheath, A step S2-4 of calculating a vector deviation, and an accumulating step S2-5.

In step S2-1 of calculating the average wind vector per altitude, the average wind vector calculation module for each altitude calculates an average wind vector up to a predetermined altitude, for example, from 0 to 3 km.

The step S2-2 of estimating the stomach motion vector estimates the stomach motion vector based on the altitude-based average wind vector calculated in the step S2-1 of calculating the average wind vector per altitude .

In the step S2-3 of calculating the vertical sheer, the vertical sheer calculating module calculates the vertical sheath of the wind vector using the upper-level elevation wind vector for each elevation.

In step S2-4 of calculating the vector deviation, the average shear is calculated after calculating the vertical shear in step S2-3, and then the average wind vector is calculated.

The accumulation step S2-5 is performed in the step S2-2 of estimating the average wind vector and the storm motion vector for each altitude calculated in the step S2-1 of calculating the average wind vector per altitude The SRH is calculated by accumulating the product of the difference of the estimated storm motion vector and the product of the vertical shear of the wind vector calculated in the step S2-3 of calculating the vertical shear to a predetermined altitude, that is, 0 to 3 km. Here, the calculated SRH is expressed by Equation (1).

It will be apparent to those skilled in the art that various modifications and variations can be made in the present invention without departing from the spirit and scope of the invention as defined in the appended claims. You will understand.

100: Wind vector calculation part 110: Raw data loading module
120: noise power level determination module 130: non-meteoric echo cancellation module
140: Atmospheric motion peak signal selection module
150: spectrum moment calculation module 160: altitude-dependent wind vector calculation module
200: SRH calculation unit
210: average wind vector calculation module by altitude
220: Storm motion vector estimation module 230: Vertical shear calculation module
240: vector deviation calculation module 250: cumulative module

Claims (8)

A wind vector calculator for removing the non-Gaussian propagation signal from the raw data of the wind profiler and reassigning the wind vector,
And an SRH calculating unit for calculating SRH (Storm Relative Helicity) based on the wind vector calculated by the wind vector calculating unit,
Wherein the wind vector calculator comprises: a raw data loading module for loading raw data of the wind profiler; a noise power level determining module for determining a noise power level of the Doppler spectrum of the raw data; A non-meteorological echo cancellation module that removes the radio signal from the altitude-dependent Doppler spectrum, a standby to determine the peak due to the atmospheric signal in the highly doped Doppler spectrum from which the non-terrestrial propagation signal exceeds the noise power level in the non- A motion peak signal selection module, a spectrum moment calculation module for determining a minimum point and a maximum point of a fast Fourier transform analysis point of the Doppler spectrum around the peak, and calculating a spectrum moment through the determined Doppler spectrum parameter, And in the spectrum moment calculation module And an altitude-dependent wind vector calculation module for calculating the wind vector of the wind profiler in the east-west direction, the north-south direction, and the vertical direction through the first beam spectrum of the four beams and the vertical beam for all the high- ,
The SRH calculator includes an altitude average wind vector calculating module for calculating an average wind vector up to a predetermined altitude based on the wind vector calculated by the wind vector calculator, A vertical sheer calculation module for calculating the vertical sheath of the wind vector by using the upper-level elevation wind vector immediately for each elevation, a vector deviation calculation module for calculating the average wind vector after calculating the vertical sheath, And an accumulation module for calculating an SRH by accumulating the difference between the average wind vector and the storm motion vector at each elevation and the product of the vertical sheer of the wind vector up to a predetermined altitude. delete delete The method according to claim 1,
In the SRH,

Lt;
remind

Is a horizontal velocity vector,
remind

Is an SRH calculation system using Wind Profiler which is a Storm motion vector of Storm. Removing the non-Gaussian propagation signal from the raw data of the wind profiler and reassigning the wind vector,
And calculating an SRH output unit SRH based on the wind vector calculated by the wind vector calculator,
Wherein the step of removing the non-Gaussian propagation signal from the wind data of the wind vector profiler and restoring the wind vector comprises the steps of: loading the wind data of the wind profiler with the source data loading module; Determining a noise power level of the Doppler spectrum by a noise power level determination module; removing the non-Gaussian echo cancellation module in the Doppler spectrum of the altitude by the non-Gaussian propagation signal exceeding the noise power level; Wherein the atmospheric motion peak signal selection module determines the peak due to the atmospheric signal in a highly Doppler spectrum in which the non-Gaussian propagation signal exceeding the noise power level is removed from the noise power level, and a fast Fourier transform of the Doppler spectrum Transform < / RTI > analysis points, Calculating a spectrum moment by a spectrum moment calculating module through a spectrum parameter and calculating a spectrum moment by a spectrum analyzer based on the wind profiler through four slope beams and a first order spectral moment of a vertical beam for all the high- The wind vector of the altitude wind vector is calculated in the east-west direction, the north-south direction and the vertical direction,
The step of calculating SRH (Storm Relative Helicity) based on the wind vector calculated by the wind vector calculating unit may include calculating an SRH (Storm Relative Helicity) based on the wind vector calculated by the wind vector calculating unit, Estimating a storm motion vector by a storm motion vector estimation module based on the altitude average wind vector, calculating a steep motion vector using a wind vector at an upper altitude level, Calculating a vertical shear of the wind vector, calculating a vertical shear, and a vector deviation calculation module calculating an average wind vector, And a cumulative module calculating the SRH by accumulating the product of the vertical sheer of the wind shear How to calculate RH. delete delete 6. The method of claim 5,
In the SRH,

Lt;
remind

Is a horizontal velocity vector,
remind

A method of SRH calculation using a wind profiler which is a storm motion vector of a storm.

Images