KR20230006193A - Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor - Google Patents
Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor Download PDFInfo
- Publication number
- KR20230006193A KR20230006193A KR1020210087075A KR20210087075A KR20230006193A KR 20230006193 A KR20230006193 A KR 20230006193A KR 1020210087075 A KR1020210087075 A KR 1020210087075A KR 20210087075 A KR20210087075 A KR 20210087075A KR 20230006193 A KR20230006193 A KR 20230006193A
- Authority
- KR
- South Korea
- Prior art keywords
- attenuation
- microwave
- band
- reflectivity
- radar reflectivity
- Prior art date
Links
- 238000002310 reflectometry Methods 0.000 title claims abstract description 185
- 238000012937 correction Methods 0.000 title claims abstract description 131
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000009977 dual effect Effects 0.000 title claims description 26
- 238000004891 communication Methods 0.000 claims abstract description 24
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 238000013480 data collection Methods 0.000 claims abstract description 18
- 230000010287 polarization Effects 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 24
- 238000005516 engineering process Methods 0.000 description 9
- 238000001556 precipitation Methods 0.000 description 9
- 238000012360 testing method Methods 0.000 description 9
- 230000003247 decreasing effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 5
- 238000000926 separation method Methods 0.000 description 4
- 238000013016 damping Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 241000533950 Leucojum Species 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
- G01S7/412—Identification of targets based on measurements of radar reflectivity based on a comparison between measured values and known or stored values
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/024—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
Description
본 발명은 엑스밴드 이중 편파 레이더 반사도 보정 장치 및 그 방법에 대한 것으로서, 특히 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치 및 그 방법에 관한 것이다.The present invention relates to an X-band dual polarization radar reflectivity correction device and method, and more particularly, to an X-band dual polarization radar reflectivity correction device and method using microwaves of multiple terrestrial communication base stations.
기상 레이더(Radio Detection And Ranging, RADAR)는 전자기파를 송신하여 목표물(강수입자)에 부딪혀서 후방 산란되어 되돌아오는 전자기파를 수신하여 강수입자의 강도와 거리를 측정하는 원격관측장비이다. 에스(S)밴드 및 씨(C)밴드 레이더에 비해 상대적으로 파장이 짧은 엑스(X)밴드 레이더는 분해능이 좋고 작은 크기의 안테나로 운용할 수 있어 강우 관측 비용을 절감할 수 있으며, 차량에 탑재가 가능하여 강우 관측 사각 지역에서 관측을 수행할 수 있는 장점이 있다. 에스밴드와 씨밴드와 같은 대형 레이더는 산 정상 부근에 대부분 구축되어 있기 때문에 도심 지역 및 저층 지역의 강우를 측정하기에는 빔 차폐 및 레이더 고도각의 한계가 있다. 미국과 일본 등에서는 CASA(Collaborative Adaptive Sensing of the Atmosphere) 프로그램이 시작된 이후로 강우 관측 공백 지역에 X밴드 레이더 네트워크를 구축하여 에스밴드 및 씨밴드 레이더의 강우 관측 한계를 극복하고자 노력하고 있다. 하지만 파장(3cm)이 짧은 엑스밴드 레이더 반사도는 강한 강수에코에 의해 감쇠가 크게 발생하여 해당 강수의 방위각의 후면으로 강수에코 탐지가 불가능하여 반사도 크기를 보정해야 한다. 엑스밴드 이중 편파 레이더가 도입된 이래로 차등위상차(Differential propagation phase, ) 변수를 활용하여 레이더 감쇠 보정 계수 를 산출하여 엑스밴드 레이더 감쇠 보정을 할 수 있는 방법이 개발되었다. 본 발명에서는 차등위상차 변수와 방송·통신·데이터 중계용으로 사용되고 있는 지상 통신 기지국(이하, 지상 마이크로파 링크) 마이크로파 감쇠 자료를 기준 값으로 이용하여 시공간적으로 변동하는 강수를 고려할 수 있는 레이더 감쇠 보정 계수 를 산출하여 감쇠가 발생한 엑스밴드 이중 편파 레이더의 반사도를 보정할 수 있는 방법을 새롭게 제안한다.Weather radar (Radio Detection And Ranging, RADAR) is a remote observation device that transmits electromagnetic waves, hits a target (precipitation particles), receives back-scattered electromagnetic waves, and measures the intensity and distance of precipitation particles. X-band radar, which has a relatively short wavelength compared to S-band and C-band radar, has good resolution and can be operated with a small-sized antenna, which can reduce rainfall observation costs and is mounted on vehicles. It has the advantage of being able to perform observations in the blind area of rainfall observation. Large radars such as S-Band and C-Band are mostly built near mountain peaks, so there are limitations in beam shielding and radar elevation angle to measure rainfall in downtown areas and low-rise areas. Since the beginning of the CASA (Collaborative Adaptive Sensing of the Atmosphere) program in the United States and Japan, X-band radar networks have been established in areas where rainfall observation is not available, and efforts have been made to overcome the limitations of S-band and C-band radar in rainfall observation. However, X-band radar reflectivity with a short wavelength (3 cm) is greatly attenuated by the strong precipitation echo, so it is impossible to detect the precipitation echo to the rear of the azimuth angle of the corresponding precipitation, so the size of the reflectivity must be corrected. Since the introduction of x-band dual polarization radar, differential propagation phase ) radar attenuation correction factor using the variable A method for calculating X-band radar attenuation correction has been developed. In the present invention, a radar attenuation correction coefficient capable of considering temporally and spatially varying precipitation by using differential phase difference variables and microwave attenuation data of a terrestrial communication base station (hereinafter, terrestrial microwave link) used for broadcasting, communication, and data relaying as reference values. We propose a new method for correcting the reflectivity of the X-band dual polarization radar in which attenuation has occurred by calculating .
본 발명의 목적은 강수에 의한 전파 감쇠가 큰 이중 편파 엑스밴드 레이더 반사도를 보정하는 보정 장치 및 그 방법을 제공하는 것이다.An object of the present invention is to provide a correction device and method for correcting dual polarization X-band radar reflectivity of which radio wave attenuation due to precipitation is large.
본 발명의 목적은 이상에서 언급한 목적으로 제한되지 않으며, 언급되지 않은 본 발명의 다른 목적 및 장점들은 하기의 설명에 의해서 이해될 수 있고, 본 발명의 실시예에 의해 보다 분명하게 이해될 것이다. 또한, 본 발명의 목적 및 장점들은 특허 청구 범위에 나타낸 수단 및 그 조합에 의해 실현될 수 있음을 쉽게 알 수 있을 것이다.The object of the present invention is not limited to the above-mentioned object, and other objects and advantages of the present invention not mentioned above can be understood by the following description and will be more clearly understood by the examples of the present invention. It will also be readily apparent that the objects and advantages of the present invention may be realized by means of the instrumentalities and combinations indicated in the claims.
본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치는, 마이크로파 링크의 송신 전력과 수신 전력을 수집하고, 마이크로파 링크 주변에 설치되어 있는 적어도 하나 이상의 기상청 자동기상관측장비 (Automatic Weather System, AWS)의 강우 감지기(rain detector)와 우량계에 일정 간격으로 기록된 누적 강우량 자료와 강우 시간 자료를 수집하는 관측 자료 수집 모듈과, 마이크로파 링크의 송신 전력과 수신 전력간의 차이를 이용하여 마이크로파 감쇠를 산출하는 마이크로파 감쇠 산출 모듈, 마이크로파 감쇠의 자기상관계수를 기반으로 강우 유무를 탐지하는 강우 유무 탐지 모듈, 및 강우 유무 탐지 모듈에서 강우 탐지 시 엑스밴드 레이더 게이트(r) 별 감쇠 크기()와, 엑스밴드 레이더의 레이더 반사도()를 이용하여 산출된 감쇠 보정된 레이더 반사도()로 엑스밴드 레이더 반사도를 보정하는 엑스밴드 레이더 반사도 보정 모듈을 가지는 반사도 보정 모듈을 포함한다.An X-band dual polarization radar reflectivity correction device using microwaves of a multi-terrestrial communication base station according to the present invention collects transmission power and reception power of a microwave link, and at least one or more automatic weather observation equipment installed around the microwave link ( Using the rain detector of the Automatic Weather System (AWS) and the observation data collection module that collects the cumulative rainfall data and rainfall time data recorded at regular intervals in the rain gauge, and the difference between the transmit power and receive power of the microwave link, A microwave attenuation calculation module that calculates microwave attenuation, a rain presence detection module that detects whether or not there is rain based on the autocorrelation coefficient of microwave attenuation, and attenuation size for each X-band radar gate (r) when the rain presence detection module detects rain ( ) and the radar reflectivity of the X-band radar ( ) attenuation-corrected radar reflectivity calculated using ) and a reflectivity correction module having an X-band radar reflectivity correction module for correcting the X-band radar reflectivity.
또한, 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법은, 관측 자료 수집 모듈이 마이크로파 링크의 송신 전력과 수신 전력을 수집하고, 마이크로파 링크 주변에 설치되어 있는 적어도 하나 이상의 기상청 자동기상관측장비 (Automatic Weather System, AWS)의 강우 감지기(rain detector)와 우량계에 일정 간격으로 기록된 누적 강우량 자료와 강우 시간 자료를 수집하는 단계와, 마이크로파 링크의 송신 전력과 수신 전력간의 차이를 이용하여 마이크로파 감쇠 산출 모듈이 마이크로파 감쇠를 산출하는 단계, 마이크로파 감쇠의 자기상관계수를 기반으로 강우 유무 탐지 모듈이 강우 유무를 탐지하는 단계, 및 반사도 보정 모듈이 비용함수를 판단하여 감쇠 보정 계수를 산출하고, 레이더 반사도반사도 와 감쇠 보정 계수 , b를 적용하여 계산한 게이트(r)별 엑스밴드 레이더 감쇠 크기로 감쇠 보정한 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정하는 단계를 포함한다.In addition, in the X-band dual polarization radar reflectivity correction method using multi-terrestrial communication base station microwaves according to the present invention, the observation data collection module collects the transmission power and reception power of the microwave link, and at least one installed around the microwave link Between the step of collecting the cumulative rainfall data and rainfall time data recorded at regular intervals in the rain detector and rain gauge of the Automatic Weather System (AWS) of the Korea Meteorological Administration and the transmission power and reception power of the microwave link Using the difference, the microwave attenuation calculation module calculates the microwave attenuation, the rain presence detection module detects whether or not there is rain based on the autocorrelation coefficient of the microwave attenuation, and the reflectance correction module determines the cost function and the attenuation correction coefficient Calculate , radar reflectivity and reflectivity and the attenuation correction factor X-band radar attenuation for each gate (r) calculated by applying , b Attenuation-corrected radar reflectivity with and correcting the X-band radar reflectivity by calculating .
여기서, 감쇠 보정된 레이더 반사도()는, 이며, 은 엑스밴드 레이더 감쇠 크기로서, 이고, , , , b는 엑스밴드 레이더의 감쇠 보정 계수이다.Here, the attenuation corrected radar reflectivity ( )Is, is, is the magnitude of the X-band radar attenuation, ego, , , , b is the attenuation correction coefficient of the X-band radar.
본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치 및 그 방법은 강우의 시공간적인 변동을 고려할 수 있는 최적의 감쇠 보정 계수를 결정할 수 있다.An apparatus and method for correcting X-band dual polarization radar reflectivity using microwaves of a multi-terrestrial communication base station according to the present invention can determine an optimal attenuation correction coefficient that can consider temporal and spatial variations in rainfall.
또한, 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치 및 그 방법은 특정 마이크로파 링크가 기계적인 원인으로 감쇠 자료를 제공할 수 없을 때에도 안정적으로 감쇠 보정 계수를 계산할 수 있는 장점이 있다.In addition, the X-band dual polarization radar reflectivity correction device and method using multi-terrestrial communication base station microwaves according to the present invention can stably calculate the attenuation correction coefficient even when a specific microwave link cannot provide attenuation data due to mechanical reasons. There are advantages to being
상술한 효과와 더불어 본 발명의 구체적인 효과는 이하 발명을 실시하기 위한 구체적인 사항을 설명하면서 함께 기술한다.In addition to the effects described above, specific effects of the present invention will be described together while explaining specific details for carrying out the present invention.
도 1은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치의 개념도이다.
도 2는 2016년 7월 5일 7시의 대류운 사례에서 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(Plan Position Indicator, PPI)를 나타낸 것이다.
도 3은 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다.
도 4는 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수 와 b를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다.
도 5는 산출한 엑스밴드 레이더 게이트 별 감쇠 크기를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다.
도 6은 기상청이 서울시에서 운영중인 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 7은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 8은 층상운 사례인 2016년 5월 2일 22시의 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(PPI)를 나타낸 것이다.
도 9는 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다.
도 10은 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다.
도 11은 산출한 엑스밴드 레이더 게이트 별 감쇠 크기를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다.
도 12는 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 13은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 14는 2016년 7월 5일의 대류운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 15는 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 16은 2016년 5월 2일의 층상운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 17은 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.
도 18은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정의 순서도이다.1 is a conceptual diagram of an X-band dual polarization radar reflectivity correction device using microwaves of a multi-terrestrial communication base station according to the present invention.
2 shows the X-band radar reflectivity (Plan Position Indicator, PPI) before attenuation correction of the Korea Institute of Civil Engineering and Building Technology in the case of convective cloud at 7:00 on July 5, 2016.
Figure 3 shows the S-band radar reflectivity (PPI) of the Korea Meteorological Administration.
4 is an optimal attenuation correction coefficient calculated by applying seven microwave link attenuations from Equation 1 to Equation 9 It shows the X-band radar attenuation field calculated by applying and b to Equation 6.
5 is the magnitude of attenuation for each calculated X-band radar gate It shows the X-band radar reflectivity obtained by applying attenuation correction to
6 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction at the AWS location operated by the Korea Meteorological Administration in Seoul.
7 shows a scatter diagram between attenuation-corrected X-band radar reflectivity and S-band radar reflectivity.
8 shows the X-band radar reflectivity (PPI) before attenuation correction of the Korea Institute of Civil Engineering and Building Technology at 22:00 on May 2, 2016, which is a case of stratified clouds.
9 shows S-band radar reflectivity (PPI) of the Korea Meteorological Administration.
10 shows an X-band radar attenuation field calculated by applying an optimal attenuation correction coefficient calculated by applying Equation 1 to Equation 9 to seven microwave link attenuations and applying Equation 6.
11 is the magnitude of attenuation for each calculated X-band radar gate It shows the X-band radar reflectivity obtained by applying attenuation correction to
12 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction at the AWS location.
13 shows a scatter diagram between attenuation-corrected X-band radar reflectivity and S-band radar reflectivity.
14 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of convective cloud on July 5, 2016.
15 shows a scatter diagram between X-band radar reflectivity and S-band radar reflectivity after attenuation correction.
16 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of stratified clouds on May 2, 2016.
17 shows a scatter diagram between X-band radar reflectivity and S-band radar reflectivity after attenuation correction.
18 is a flowchart of X-band dual polarization radar reflectivity correction using microwaves of a multi-terrestrial communication base station according to the present invention.
전술한 목적, 특징 및 장점은 첨부된 도면을 참조하여 상세하게 후술되며, 이에 따라 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 본 발명의 기술적 사상을 용이하게 실시할 수 있을 것이다. 본 발명을 설명함에 있어서 본 발명과 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 상세한 설명을 생략한다. 이하, 첨부된 도면을 참조하여 본 발명에 따른 바람직한 실시예를 상세히 설명하기로 한다. 도면에서 동일한 참조부호는 동일 또는 유사한 구성요소를 가리키는 것으로 사용된다.The above objects, features and advantages will be described later in detail with reference to the accompanying drawings, and accordingly, those skilled in the art to which the present invention belongs will be able to easily implement the technical spirit of the present invention. In describing the present invention, if it is determined that the detailed description of the known technology related to the present invention may unnecessarily obscure the subject matter of the present invention, the detailed description will be omitted. Hereinafter, preferred embodiments according to the present invention will be described in detail with reference to the accompanying drawings. In the drawings, the same reference numerals are used to indicate the same or similar components.
도 1은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치의 개념도이다.1 is a conceptual diagram of an X-band dual polarization radar reflectivity correction device using microwaves of a multi-terrestrial communication base station according to the present invention.
본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치는 도 1에 도시된 바와 같이, 관측 자료 수집 모듈(100)과, 마이크로파 링크의 송신 전력과 수신 전력의 차이를 이용하여 마이크로파 감쇠(rainfall-caused attenuation, )를 산출하는 마이크로파 감쇠 산출 모듈(200), 마이크로파 감쇠 자료를 기반으로 강우 유무를 탐지하는 강우 유무 탐지 모듈(300), 및 감쇠 보정 계수로 엑스밴드 레이더의 반사도를 보정하는 반사도 보정 모듈(400)을 포함한다.As shown in FIG. 1, the X-band dual polarization radar reflectivity correction apparatus using microwaves of a multi-terrestrial communication base station according to the present invention uses the difference between the transmission power and the reception power of the observation
표 1은 본 발명에 따른 관측 자료 수집 모듈에 입력될 7개의 마이크로파 링크의 특성과 각각의 마이크로파 링크의 수신 안테나에서 가장 근접하게 설치된 자동기상관측장비(AWS) 지점명을 나타낸 것이다.Table 1 shows the characteristics of seven microwave links to be input to the observation data collection module according to the present invention and the name of an automatic weather station (AWS) point installed closest to the receiving antenna of each microwave link.
name Link
name
AWSNearest
AWS
(GHz)Frequency
(GHz)
(dBm)Power
(dBm)
본 발명은 0.01㏈의 해상도로 15초 간격으로 7개의 마이크로파 링크의 송신 전력과 수신 전력을 집록하고, 상대적으로 파장이 긴 10㎓ 이하의 마이크로파를 활용하여 엑스밴드 레이더 반사도 보정을 위한 기준값으로 활용한다.The present invention collects the transmit power and receive power of 7 microwave links at 15 second intervals with a resolution of 0.01 dB, and utilizes microwaves with relatively long wavelengths of 10 GHz or less to use them as reference values for X-band radar reflectivity correction. .
표 2는 한국건설기술연구원의 엑스밴드 이중 편파 레이더와 기상청의 에스밴드 단일 편파 레이더의 특성을 나타낸 것이다.Table 2 shows the characteristics of the Korea Institute of Civil Engineering and Building Technology's X-band dual polarization radar and the Korea Meteorological Administration's S-band single polarization radar.
관측 자료 수집 모듈(100)은 엑스밴드 레이더 반사도 보정을 위해 관측 자료를 입력받아 수집한다. 여기서, 관측 자료 수집 모듈(100)이 수집하는 관측 자료는 마이크로파 링크의 송신 전력과 수신 전력, 마이크로파 링크 주변에 설치된 기상청의 자동기상관측장비(Automatic Weather System, AWS)의 강우 감지기에 기록된 강우 시간과 우량계에 기록된 누적 강우량과 강우 시간 자료, 엑스밴드 레이더의 반사도 및 에스밴드 레이더의 반사도를 포함한다. 또한, 이를 수집하기 위해서, 관측 자료 수집 모듈(100)은 마이크로파 링크의 송신 전력과 수신 전력을 수집하는 마이크로파 링크 자료 수집 모듈과, 마이크로파 링크 주변에 설치된 기상청 자동기상관측장비의 강우 감지기에 기록된 강우 지속 시간 자료를 수집하는 강우 감지기 자료 수집 모듈, 한국건설기술연구원의 엑스밴드 레이더의 반사도와 차등위상차 변수를 수집하는 엑스밴드 레이더 자료 수집 모듈, 및 기상청의 광덕산 에스밴드 레이더의 반사도를 수집하는 에스밴드 레이더 자료 수집 모듈을 포함한다.The observation
마이크로파 감쇠 산출 모듈(200)은 마이크로파 링크의 송신 전력과 수신 전력의 차이를 이용하여 마이크로파 감쇠(rainfall-caused attenuation, )를 산출한다. 여기서, 마이크로파 감쇠 산출 모듈(200)은 15초 해상도의 마이크로파 링크의 송신 전력과 수신 전력 간의 차이를 이용하여 마이크로파 감쇠를 산출한다.The microwave
강우 유무 탐지 모듈(300)은 마이크로파 감쇠 자료를 기반으로 강우 유무를 탐지한다. 이를 위해서, 강우 유무 탐지 모듈(300)은 강우 특성에 따른 마이크로파 감쇠의 자기상관계수가 0.65 이상을 보인 최소 지연 시간을 창 크기(moving time window, Wt)로 결정하는 창 크기 산출 모듈(310)과, 결정된 창 크기 동안의 평균 마이크로파 감쇠()를 산출하는 평균 마이크로파 감쇠 산출 모듈(320), 결정된 창 크기 동안의 마이크로파 감쇠와 평균 마이크로파 감쇠()간의 차이를 이용하여 마이크로파 감쇠의 검정 통계량을 산출하는 마이크로파 감쇠의 검정 통계량 산출 모듈(330), 마이크로파 감쇠를 이용한 강우 유무를 1차적으로 판단하기 위해 필요한 임계 값을 결정하는 임계 값 산출 모듈, 마이크로파 감쇠의 검정 통계량과 임계 값 간의 초과 여부를 통해 1차적으로 강우 유무를 탐지하는 1차 강우 유무 탐지 모듈(350), 1차적으로 결정된 강우 유무 탐지 결과와 특정 창 크기 동안의 평균 마이크로파 감쇠()를 이용하여 마이크로파 기준 감쇠(Ad)를 실시간으로 산출하는 마이크로파 기준 감쇠 산출 모듈(360), 마이크로파 감쇠와 마이크로파 기준 감쇠 간의 초과 여부를 통해 2차적으로 강우 유무를 탐지하는 2차 강우 유무 탐지 모듈(370)을 포함한다.The rain presence/absence detection module 300 detects the presence/absence of rain based on the microwave attenuation data. To this end, the rain presence detection module 300 includes a window
창 크기 산출 모듈(310)은 마이크로파 감쇠 산출 모듈(200)에서 산출된 마이크로파 감쇠 자료를 기반으로 마이크로파 감쇠의 자기상관계수가 통계학적으로 유의미한 수치, 예를 들어, 0.65를 보이는 최소 지연 시간을 창 크기(moving time window, Wt)로 결정한다.The window
여기서, 자기상관계수는 아래의 수학식 1과 같다.Here, the autocorrelation coefficient is equal to Equation 1 below.
수학식 1에서, 이고, 는 지연시간을 나타낸다.In Equation 1, ego, represents the delay time.
평균 마이크로파 감쇠 산출 모듈(320)은 창 크기 산출 모듈(310)에서 결정된 창 크기 동안의 마이크로파 감쇠를 평균하여 1분 간격으로 평균 마이크로파 감쇠()를 산출한다.The average microwave
여기서, 평균 마이크로파 감쇠()는 창 크기(Wt)내의 마이크로파 감쇠 크기의 평균이며, 아래의 수학식 2와 같다.Here, the average microwave attenuation ( ) is the average of the microwave attenuation in the window size (Wt), and is shown in Equation 2 below.
수학식 2에서, 는 실시간 측정되는 마이크로파 감쇠 크기(k) 이며, 는 창 크기(Wt)내의 마이크로파 감쇠 자료 개수이다.In Equation 2, is the microwave attenuation magnitude (k) measured in real time, is the number of microwave attenuation data within the window size (Wt).
마이크로파 감쇠의 검정 통계량 산출 모듈(330)은 마이크로파 감쇠 산출 모듈(200)에서 산출된 마이크로파 감쇠()와, 평균 마이크로파 감쇠 산출 모듈(320)에서 산출된 평균 마이크로파 감쇠()의 차이를 이용하여 마이크로파 감쇠의 검정 통계량을 산출한다. 여기서, 마이크로파 감쇠의 검정 통계량()은 아래의 수학식 3과 같다.The microwave attenuation test
수학식 3에서, 마이크로파 감쇠는 비정상성(nonstationary)이기 때문에 마이크로파 감쇠의 검정 통계량()과 마이크로파 감쇠의 표준 편차는 상호 수렴하지 않는다. 는 창 크기(Wt)내의 마이크로파 감쇠의 평균이다.In Equation 3, since microwave attenuation is nonstationary, the test statistic of microwave attenuation ( ) and the standard deviation of microwave attenuation do not converge. is the average of the microwave attenuation within the window size (Wt).
임계 값 산출 모듈(340)은 마이크로파 감쇠를 이용한 강우 유무를 1차 판단하기 위해 필요한 임계값()을 산출한다. 본 실시예는 강우 감지기 자료와 비교하여 90% 이상의 일치를 보이는 임계값을 최종 임계값()으로 결정하는 것을 예시하며, 이는 0.09인 것을 예시한다.The threshold
1차 강우 유무 탐지 모듈(350)은 마이크로파 감쇠의 검정 통계량()과, 임계 값 산출 모듈(340)에서 산출된 임계값 간의 초과 여부를 통해 1차적으로 강우 유무를 탐지한다. 즉, 창 크기(Wt)가 결정되면 임계값()이 결정되고, 아래의 수학식 4에 의해 1차적으로 마이크로파 감쇠를 이용한 강우 유무, 즉, 1차 강우 유무가 결정된다.The primary rainfall presence/
마이크로파 기준 감쇠 산출 모듈(360)은 평균 마이크로파 감쇠 산출 모듈(320)에서 산출된 평균 마이크로파 감쇠()와, 1차 강우 유무 탐지 모듈(350)에서 탐지된 1차 강우 유무 탐지 결과를 이용하여 마이크로파 기준 감쇠(Ad)를 실시간으로 산출한다. 여기서, 마이크로파 기준 감쇠(Ad)는 아래의 수학식 5와 같이 산출된다.The microwave reference
마이크로파 감쇠 크기가 실시간으로 산출된 마이크로파 기준 감쇠(Ad)를 초과할 때 강우에 의해 야기된 것으로 결정한다. 마이크로파 감쇠와 마이크로파 기준 감쇠(Ad)간의 초과 여부에 따라 1차적으로 결정된 강우 유무 결과를 보정하는 2차적으로 강우 유무를 최종 결정한다. 강우에 의해 야기된 마이크로파 감쇠 산출 모듈(200)은 마이크로파 감쇠 산출 모듈(200)에서 산출되는 마이크로파 감쇠()가 마이크로파 기준 감쇠 산출 모듈(360)에서 산출된 마이크로파 기준 감쇠와의 차이가 마이크로파 감쇠 크기의 2.5%보다 클 때 강우에 의해 야기된 마이크로파 감쇠()로 결정한다.When the magnitude of the microwave attenuation exceeds the microwave reference attenuation (Ad) calculated in real time, it is determined to be caused by rainfall. The presence or absence of rainfall is finally determined secondarily by correcting the result of whether or not there is rain primarily determined according to whether the microwave attenuation and the microwave standard attenuation (Ad) are exceeded. The microwave
2차 강우 유무 탐지 모듈(370)은 마이크로파 감쇠()와 마이크로파 기준 감쇠(Ad)간의 초과 여부를 통해 결정된 1차 강우 유무 탐지 결과를 보정하는 2차 강우 유무를 탐지한다. 여기서, 마이크로파 기준 감쇠(Ad)는 창 크기 내의 마이크로파 감쇠의 평균이며, 창 크기 내의 마이크로파 기준 감쇠와 큰 차이를 보이는 마이크로파 감쇠()는 강우가 내렸을 때의 마이크로파 감쇠()로 결정하기 위해 사용된다. 최종적으로 2차 강우 유무 탐지 모듈(370)에서 마이크로파 감쇠()를 이용한 강우 유무 탐지를 결정한다. 또한, 마이크로웨이브 링크 송신 및 수신 안테나로부터 가장 근접하게 설치된 자동기상관측장비(AWS)의 강우 감지기에 강우로 기록된 시간 자료와 비교하여 이를 검증한다.The secondary rainfall
반사도 보정 모듈(400)은 감쇠 보정 계수를 산출하여 엑스밴드 레이더의 반사도를 보정한다. 이를 위해서, 반사도 보정 모듈(400)은 마이크로파 감쇠 결정 모듈(410)과, 엑스밴드 레이더 게이트(r) 별 감쇠 크기()를 연산하는 엑스밴드 레이더 게이트 별 감쇠 크기 연산 모듈(420)과, 엑스밴드 레이더 평균 감쇠 크기()를 산출하는 엑스밴드 레이더 평균 감쇠 크기 산출 모듈(430), 레이더 주파수로 변환된 마이크로파 감쇠 크기()를 산출하는 마이크로파 감쇠 크기 산출 모듈(440), 비용함수()를 산출하는 비용함수 산출 모듈(450), 및 엑스밴드 레이더 반사도를 보정하는 엑스밴드 레이더 반사도 보정 모듈(460)을 포함한다.The
마이크로파 감쇠 결정 모듈(410)은 마이크로파 기준 감쇠를 초과하는 마이크로파 감쇠를 강우에 의해 야기된 마이크로파 감쇠()로 결정한다.The microwave
엑스밴드 레이더 게이트 별 감쇠 크기 연산 모듈(420)은 엑스밴드 레이더 게이트(r) 별 감쇠 크기()를 연산한다. 여기서, 엑스밴드 레이더 게이트 별 감쇠 크기는 아래의 수학식 6과 같다. 이중 편파 레이더가 도입된 이래로 엑스밴드 레이더 감쇠 크기()는 차등 위상차()의 합이 전파 경로를 따라 축적된 값으로 가정하고 계산한다. 차등 위상차는 수평 편파 위상과 수직 편파 위상의 차이를 의미하는데 전파 경로 누적 감쇠는 차등 위상차의 총합이라는 가정으로 감쇠를 추정한다.The attenuation size per X-band radar
수학식 6에서, 와 은 특정 레이더 방위각에서 반사도가 30dBZ 이상을 보이는 첫 번째 게이트와 마지막 게이트를 의미한다. r은 해당되는 게이트이다. 은 레이더의 게이트 간격(km)이며, 레이더 전자기파의 왕복 거리를 고려하기 위해서 2를 곱한다. 엑스밴드 레이더 게이트(r)별 감쇠 크기()는 아래의 수학식 7과 같이 계산한다. In Equation 6, Wow denotes the first and last gates whose reflectivity is greater than 30 dBZ at a specific radar azimuth. r is the corresponding gate. is the gate distance (km) of the radar, multiplied by 2 to consider the round-trip distance of radar electromagnetic waves. Attenuation size for each X-band radar gate (r) ( ) is calculated as in Equation 7 below.
수학식 7에서, 이며, 이다. 또한, 는 측정된 엑스밴드 레이더의 반사도 자료이며, , b는 경험적으로 결정되는 레이더 감쇠 보정 계수이다. 는 0.13에서 0.34까지 변화하며, 본 실시예는 로 0.14를 예시한다.In Equation 7, is, to be. also, is the reflectivity data of the measured X-band radar, , b is the radar attenuation correction factor determined empirically. changes from 0.13 to 0.34, and this embodiment 0.14 is exemplified by
엑스밴드 레이더 평균 감쇠 크기 산출 모듈(430)은 최적화된 , b 계수를 결정하기 위해서 아래의 수학식 8과 같이 마이크로파 링크 경로 길이에 해당하는 엑스밴드 레이더 평균 감쇠 크기()를 계산한다.The X-band radar average attenuation
수학식 8에서, 은 마이크로파 링크의 경로 길이에 해당하는 가장 인접한 레이더 게이트의 전체 개수이다.In
마이크로파 감쇠 크기 산출 모듈(440)은 레이더 주파수로 변환된 마이크로파 감쇠 크기()를 산출한다. 마이크로파 링크의 주파수와 엑스밴드 레이더 주파수가 다를 경우, 마이크로파 링크의 감쇠 크기(에 마이크로파 링크 주파수와 엑스밴드 레이더 주파수의 비를 곱하여 마이크로파 감쇠 크기()를 아래의 수학식 9와 같이 변환하고, 레이더 감쇠 보정 계수 를 결정하기 위한 기준 값으로 사용한다.The microwave attenuation
비용함수 산출 모듈(450)은 엑스밴드 레이더의 평균 감쇠 크기()와 변환된 마이크로파 감쇠 크기()간의 차이를 계산하여 아래의 수학식 10과 같이 비용함수()가 최소일 때의 감쇠 보정 계수를 최적의 감쇠 보정 계수로 결정한다.The cost
최적의 감쇠 보정 계수 를 산출하기 위해서 황금분할탐색방법을 적용하였다. 황금분할탐색방법으로 최적의 감쇠 보정 계수 를 산출하는 과정은 먼저 두 감쇠 보정 계수 (, ) 구간을 선택하고 구간 내의 감쇠 보정 계수(, )를 결정하기 위해 의 차이에 황금 비율인 0.618을 곱한 분리 거리(d)를 계산하고, 는 와 분리거리(d)의 차이 값으로 정의하고, 은 과 분리거리(d)의 합으로 정의하였다. 만약 를 레이더의 평균 감쇠 크기()계산에 필요한 감쇠 보정 계수 값으로 사용하고 수학식 10에 적용하여 산출한 비용함수()가 을 사용하여 산출한 비용함수()보다 크다면, 을 값으로 대체하고, 반대이면 를 값으로 대체하면 두 감쇠 보정 계수 (, )의 차이는 분리 거리(d)만큼 좁아진다. 반복 과정을 몇 번 거치면 두 감쇠 보정 계수 (, )의 거리는 0에 근접하게 되는데 그 때의 감쇠 보정 계수를 최적의 감쇠 보종 계수로 결정하였다. 감쇠 보정 계수 의 초기값으로 과 를 각각 0.1, 0.5으로 설정하였고 b는 1.0으로 상수로 사용하였다.Optimal attenuation correction factor To calculate , the golden partition search method was applied. Optimal damping correction coefficient by golden partition search method The process of calculating is first two attenuation correction coefficients ( , ) section and select the attenuation correction factor ( , ) to determine Calculate the separation distance (d) by multiplying the difference of 0.618, which is the golden ratio, Is Defined as the difference between and the separation distance (d), silver It was defined as the sum of and the separation distance (d). if the average attenuation magnitude of the radar ( ) attenuation correction factor required to calculate The cost function calculated by using it as a value and applying Equation 10 ( )go The cost function calculated using ( ) is greater than second Substitute the value, and vice versa cast Substituting the values for the two attenuation correction factors ( , ) is narrowed by the separation distance (d). After a few iterations, the two attenuation correction factors ( , ) is close to 0, and the attenuation correction coefficient at that time was determined as the optimal attenuation correction coefficient. Attenuation correction factor with an initial value of class were set to 0.1 and 0.5, respectively, and b was used as a constant at 1.0.
엑스밴드 레이더 반사도 보정 모듈(460)은 레이더 반사도 와 감쇠 보정 계수 , b를 적용하여 계산한 게이트(r)별 엑스밴드 레이더 감쇠 크기를 이용하여 엑스밴드 레이더 반사도를 보정한다. 이는 아래의 수학식 11과 같이 감쇠 보정한 레이더 반사도를 산출하여 수행할 수 있다.X-band radar
도 2는 2016년 7월 5일 7시의 대류운 사례에서 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(Plan Position Indicator, PPI)를 나타낸 것이다. 도 4는 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다. 도 5는 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수 와 b를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다. 도 5는 산출한 엑스밴드 레이더 게이트 별 감쇠 크기 를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다. 도 6은 기상청이 서울시에서 운영중인 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 도 7은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 도 2는 2016년 7월 5일 7시의 대류운 사례에서 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(Plan Position Indicator, PPI)를 나타낸 것이다. 층상운 사례에 비해 40 dBZ이상의 국지적인 강우 셀의 반사도를 보였다. 도 3은 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다. 도 4는 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수 와 b를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다. 최대 15 dB이상의 감쇠 크기를 보였다. 도 5는 산출한 엑스밴드 레이더 게이트 별 감쇠 크기를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다. 감쇠 보정 후 에스밴드 레이더 반사도와 공간 분포가 상당히 일치하는 것을 확인할 수 있다. 도 6은 기상청이 서울시에서 운영중인 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 절대 평균 오차(Mean Absolute Error, MAE)는 16.2 dB, RMSE는 17.7 dB 그리고 상관계수는 0.7를 보여 에스밴드 레이더 반사도에 비해 강우 감쇠로 인해 과소 측정된 것을 알 수 있다. 도 7은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것인데, 절대 평균 오차가 16.2 dB에서 6.7 dB로 크게 감소하였고 RMSE는 17.7 dB에서 8.2 dB로 감소하였으며 상관 계수는 0.7에서 0.8로 증가한 것을 확인하였다.2 shows the X-band radar reflectivity (Plan Position Indicator, PPI) before attenuation correction of the Korea Institute of Civil Engineering and Building Technology in the case of convective cloud at 7:00 on July 5, 2016. Figure 4 shows the S-band radar reflectivity (PPI) of the Korea Meteorological Administration. 5 is an optimal attenuation correction coefficient calculated by applying seven microwave link attenuations from Equation 1 to Equation 9 It shows the X-band radar attenuation field calculated by applying and b to Equation 6. 5 is the magnitude of attenuation for each calculated X-band radar gate It shows the X-band radar reflectivity obtained by applying attenuation correction to
도 8은 층상운 사례인 2016년 5월 2일 22시의 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(PPI)를 나타낸 것이고, 도 9는 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다. 도 10은 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다. 도 11은 산출한 엑스밴드 레이더 게이트 별 감쇠 크기 를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다. 도 12는 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 도 13은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.8 shows the X-band radar reflectivity (PPI) before attenuation correction of the Korea Institute of Civil Engineering and Building Technology at 22:00 on May 2, 2016, which is a case of layered cloud, and FIG. 9 shows the S-band radar reflectivity (PPI) of the Korea Meteorological Administration. . 10 shows an X-band radar attenuation field calculated by applying an optimal attenuation correction coefficient calculated by applying Equation 1 to Equation 9 to seven microwave link attenuations and applying Equation 6. 11 is the magnitude of attenuation for each calculated X-band radar gate It shows the X-band radar reflectivity obtained by applying attenuation correction to
도 8은 층상운 사례인 2016년 5월 2일 22시의 한국건설기술연구원의 감쇠 보정 전의 엑스밴드 레이더 반사도(PPI)를 나타낸 것이다. 도 9는 기상청의 에스밴드 레이더 반사도(PPI)를 나타낸 것이다. 에스밴드 레이더 전자기파는 강수에 의해 감쇠가 크지 않기 때문에 엑스밴드 레이더 반사도에 비해 그 크기가 상대적으로 큰 것을 확인할 수 있다. 도 10은 7개의 마이크로파 링크 감쇠를 수학식 1에서 수학식 9까지 적용하여 산출한 최적의 감쇠 보정 계수를 수학식 6에 적용하여 산출한 엑스밴드 레이더 감쇠장을 나타낸 것이다. 최대 12 dB이하의 감쇠 크기를 보였다. 도 11은 산출한 엑스밴드 레이더 게이트 별 감쇠 크기 를 수학식 10에 적용하여 감쇠 보정한 엑스밴드 레이더 반사도를 나타낸 것이다. 감쇠 보정 후 에스밴드 레이더 반사도와 공간 분포가 상당히 일치하는 것을 확인할 수 있다. 도 12는 AWS 위치에서의 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 절대 평균 오차는 8.9 dB, RMSE는 10.3 Db 그리고 상관 계수는 0.3을 보여 과소 측정된 것을 알 수 있다. 도 13은 감쇠 보정한 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것인데, 절대 평균 오차가 8.9 dB에서 3.1 dB 감소하였고 RMSE는 10.3 dB에서 3.7 dB로 감소하였으며 상관 계수는 0.3에서 0.7로 증가하였다.8 shows the X-band radar reflectivity (PPI) before attenuation correction of the Korea Institute of Civil Engineering and Building Technology at 22:00 on May 2, 2016, which is a case of stratified clouds. 9 shows S-band radar reflectivity (PPI) of the Korea Meteorological Administration. Since the S-band radar electromagnetic wave is not greatly attenuated by precipitation, it can be seen that the size is relatively large compared to the X-band radar reflectivity. 10 shows an X-band radar attenuation field calculated by applying an optimal attenuation correction coefficient calculated by applying Equation 1 to Equation 9 to seven microwave link attenuations and applying Equation 6. The maximum attenuation was less than 12 dB. 11 is the magnitude of attenuation for each calculated X-band radar gate It shows the X-band radar reflectivity obtained by applying attenuation correction to
도 14는 2016년 7월 5일의 대류운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이며, 도 15는 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 도 16은 2016년 5월 2일의 층상운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 도 17은 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다.14 shows the scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of convective cloud on July 5, 2016, and FIG. 15 shows the scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity after attenuation correction it is shown 16 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of stratified clouds on May 2, 2016. 17 shows a scatter diagram between X-band radar reflectivity and S-band radar reflectivity after attenuation correction.
도 14는 2016년 7월 5일의 대류운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 검은색 기호는 엑스밴드 레이더 중심으로부터 반경 10 km 이내에 위치한 AWS에서의 반사도 간의 산포도를 나타낸 것이고 노랑색 기호는 반경 20 km 이내의 반사도 간의 산포도, 파란색 기호는 반경 30 km, 빨강색 기호는 40 km 이내의 반사도 간의 산포도를 나타낸 것이다. 감쇠 보정 전의 엑스밴드 레이더와 에스밴드 레이더 반사도 간의 절대 평균 오차는 17.9 dB, RMSE는 20.6 dB 그리고 상관 계수는 0.4로 강우에 의한 감쇠로 엑스밴드 레이더 반사도가 낮은 것을 알 수 있다. 도 15는 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 감쇠 보정 후의 엑스밴드 레이더와 에스밴드 레이더 반사도 간의 절대 오차는 17.9 dB에서8.4 dB로 크게 감소하였고, RMSE는 20.6 dB 에서 10.8 dB로 감소하였으며 상관 계수는 0.4에서 0.7로 증가한 것을 확인할 수 있다. 도 16은 2016년 5월 2일의 층상운 사례에서 감쇠 보정 전의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 감쇠 보정 전의 엑스밴드 레이더와 에스밴드 레이더 반사도 간의 절대 평균 오차는 6.4 dB, RMSE는 8.1 dB이었고 상관 계수는 0.6이었다. 도 17은 감쇠 보정 후의 엑스밴드 레이더 반사도와 에스밴드 레이더 반사도 간의 산포도를 나타낸 것이다. 감쇠 보정 후의 엑스밴드 레이더와 에스밴드 레이더 반사도 간의 절대 평균 오차는 감쇠 보정 전의 6.4 dB 에서 6.3 dB 로 감소하였고, RMSE는 8.1 dB에서 7.9 dB로 감소하였다. 층상운 사례에 비해 대류운 사례에서 감쇠 보정 효과가 더 큰 것을 확인할 수 있다.14 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of convective cloud on July 5, 2016. The black symbol shows the scatter plot between reflectivities in AWS located within a radius of 10 km from the center of the X-band radar, the yellow symbol shows the scatter plot between reflectivities within a radius of 20 km, the blue symbol shows the scatter plot between the reflectivities within a radius of 30 km, and the red symbol shows the scatter plot within a radius of 40 km. It shows the scatter plot between reflectances. The absolute average error between the X-band radar and the S-band radar reflectivity before attenuation correction was 17.9 dB, the RMSE was 20.6 dB, and the correlation coefficient was 0.4. 15 shows a scatter diagram between X-band radar reflectivity and S-band radar reflectivity after attenuation correction. After attenuation correction, the absolute error between the X-band radar and the S-band radar reflectivity significantly decreased from 17.9 dB to 8.4 dB, the RMSE decreased from 20.6 dB to 10.8 dB, and the correlation coefficient increased from 0.4 to 0.7. 16 shows a scatter diagram between the X-band radar reflectivity and the S-band radar reflectivity before attenuation correction in the case of stratified clouds on May 2, 2016. The absolute average error between X-band radar and S-band radar reflectivity before attenuation correction was 6.4 dB, RMSE was 8.1 dB, and the correlation coefficient was 0.6. 17 shows a scatter diagram between X-band radar reflectivity and S-band radar reflectivity after attenuation correction. The absolute average error between the X-band radar and S-band radar reflectivity after attenuation correction decreased from 6.4 dB before attenuation correction to 6.3 dB, and the RMSE decreased from 8.1 dB to 7.9 dB. It can be seen that the attenuation correction effect is greater in the convective cloud case than in the stratified cloud case.
다음은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치를 이용하여 감쇠 보정한 엑스밴드 레이더 반사도를 검증하였다.Next, the attenuation-corrected X-band radar reflectivity was verified using the X-band dual polarization radar reflectivity correction device using microwaves of a multi-terrestrial communication base station according to the present invention.
감쇠 보정한 엑스밴드 레이더 반사도를 검증하기 위해서 에스밴드 레이더 반사도 자료와 비교 검증을 하였는데, 레이더 빔의 높이 차이와 엑스밴드 레이더와 에스밴드 레이더의 파장의 차이를 고려하였고 그 외의 차이는 무시할 수 있는 수준으로 가정하였다. 먼저 두 레이더 간의 수평 거리 차이가 수 백 km이기 때문에 지구 곡면 효과로 빔 높이의 차이가 발생한다. 따라서 본 발명에서는 표 3과 같이 에스밴드 레이더 고도각(0~1.7°) 빔의 높이를 조정하여 엑스밴드 레이더 고도각(5°)의 빔의 차이를 최소화하였다.In order to verify the attenuation-corrected X-band radar reflectivity, comparison and verification were performed with the S-band radar reflectivity data. was assumed. First, since the difference in horizontal distance between the two radars is several hundred km, a difference in beam height occurs due to the earth curvature effect. Therefore, in the present invention, as shown in Table 3, the height of the S-band radar elevation angle (0 to 1.7°) beam was adjusted to minimize the difference between the beams of the X-band radar elevation angle (5°).
SiteAWS
Site
RBH (km)KICT
RBH (km)
RBH (km) /Elevation angle(°)GDK
RBH (km) / Elevation angle (°)
SiteAWS
Site
RBH (km)KICT
RBH (km)
RBH (km) /Elevation angl(°)GDK
RBH (km) / Elevation angle (°)
엑스밴드 레이더 전자기파의 파장은 2.5~4cm, 에스밴드 레이더 전자기파의 파장은 8~15cm이고 강수 입자의 직경은 수 mm이고 눈송이나 우박은 수 mm에서 수 cm이기 때문에, 입자의 크기가 파장보다 작거나 비슷하다. 전파의 파장이 수 cm인 기상 레이더로 관측할 때 수 mm인 빗방울에서는 레일리 산란이, 직경이 수 mm에서 수 cm인 눈송이나 우박에서는 레일리 산란 또는 미 산란이 일어난다. 아래의 수학식 12를 이용하여 에스밴드 레이더 반사도를 엑스밴드 레이더 반사도로 변환하였다.The wavelength of the X-band radar electromagnetic wave is 2.5 to 4 cm, the wavelength of the S-band radar electromagnetic wave is 8 to 15 cm, and the diameter of precipitation particles is several mm, and snowflakes and hail are several mm to several cm. Similar. When observed by a weather radar with a wavelength of several cm, Rayleigh scattering occurs in raindrops with a wavelength of several mm, and Rayleigh scattering or Mie scattering occurs in snowflakes or hailstones with a diameter of several mm to several cm. The S-band radar reflectivity was converted to the X-band radar
수학식 12에서, 는 에스밴드 레이더 반사도이며, 는 에스밴드 레이더 반사도()가 엑스밴드 레이더 파장으로 변환되었을 때의 반사도이다. 에스밴드 레이더 반사도가 25dBZ 이하일 때는 식에 적용하여 반사도 변환을 하였고, 에스밴드 레이더 반사도가 25dBZ 에서 45dBZ 사이일 경우에는 식을 적용하였고, 45dBZ 이상일 경우에는 식을 적용하여 변환하였다.In
레이더 빔의 높이 차이와 엑스밴드와 에스밴드 레이더 간의 파장 차이를 고려하였고 그 외의 차이는 무시할 수 있는 수준으로 가정하였다. 정량적 검증은 아래와 같은 수식을 이용하여 통계적으로 수행되었다.The height difference of the radar beam and the wavelength difference between the X-band and S-band radars were considered, and other differences were assumed to be negligible. Quantitative verification was performed statistically using the following formula.
본 발명에서는 강우의 시공간적인 변동성을 고려한 최적의 엑스밴드 레이더 감쇠 보정 계수 를 결정하기 위해 레이더 반경 내에 위치한 7개의 마이크로파 링크 감쇠 자료를 모두 이용하는 방법을 제안하였다. 먼저, 각각의 마이크로파 링크 경로 길이에 해당하는 강우에 의한 엑스밴드 레이더의 평균 감쇠 크기(, , , , , , )를 엑스밴드 레이더 주파수로 변환한 7개의 마이크로파 링크 감쇠( )와 링크 길이에 해당하는 엑스밴드 레이더 반사도()를 이용하여 계산하였고, 각각의 비용함수()가 최소(=0.25)가 될 때의 감쇠 보정 계수(1, 2, 3, 4, 5, 6, 7)를 각각 산출하였다. 시공간적으로 변동성이 적은 층상운 강우 사례에서는 7개의 감쇠 보정 계수들의 평균값을 사용하였고, 대류운 사례에서는 7개의 감쇠 보정 계수들의 최대값을 최적의 감쇠 보정 계수로 결정하였다. 본 발명에서는 강우의 시공간적인 변동을 고려할 수 있는 최적의 감쇠 보정 계수를 결정할 수 있고, 특정 마이크로파 링크가 기계적인 원인으로 감쇠 자료를 제공할 수 없을 때에도 안정적으로 감쇠 보정 계수를 계산할 수 있는 장점이 있다.In the present invention, the optimal X-band radar attenuation correction coefficient considering the temporal and spatial variability of rainfall To determine , a method using all 7 microwave link attenuation data located within the radar radius is proposed. First, the average attenuation of the X-band radar by rainfall corresponding to each microwave link path length ( , , , , , , ) to the X-band radar frequency attenuating the seven microwave links ( ) and the X-band radar reflectivity corresponding to the link length ( ), and each cost function ( ) becomes the minimum (=0.25) attenuation correction factor ( One, 2, 3, 4, 5, 6, 7) were calculated respectively. In the case of stratiform rainfall with little temporal and spatial variability, the average value of 7 attenuation correction coefficients was used, and in the case of convective cloud, the maximum value of the 7 attenuation correction coefficients was determined as the optimal attenuation correction coefficient. In the present invention, it is possible to determine the optimal attenuation correction coefficient that can consider the temporal and spatial fluctuations of rainfall, and to stably calculate the attenuation correction coefficient even when a specific microwave link cannot provide attenuation data due to mechanical reasons. .
다음은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법에 대해 도면을 참조하여 설명한다. 후술될 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법의 각 단계는 전술된 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치를 이용한다.Next, a method for correcting X-band dual polarization radar reflectivity using microwaves of a multi-terrestrial communication base station according to the present invention will be described with reference to the drawings. Each step of the X-band dual polarization radar reflectivity correction method using microwaves of a multi-terrestrial communication base station according to the present invention, which will be described later, uses the above-described X-band dual polarization radar reflectivity correction device using microwaves of a multi-terrestrial communication base station.
도 18은 본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법의 순서도이다.18 is a flowchart of an X-band dual polarization radar reflectivity correction method using microwaves of a multi-terrestrial communication base station according to the present invention.
본 발명에 따른 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법은 도 18을 참조하면, 관측 자료를 수집하는 단계(S1, 미도시)와, 마이크로파 감쇠 자료를 기반으로 강우 유무를 탐지하는 단계(S2, 미도시), 엑스밴드 레이더의 반사도를 보정하는 단계(S3)를 포함한다.Referring to FIG. 18, the X-band dual polarization radar reflectivity correction method using microwaves of a multi-terrestrial communication base station according to the present invention includes collecting observation data (S1, not shown), and determining whether or not there is rainfall based on microwave attenuation data. A step of detecting (S2, not shown), and a step of correcting the reflectivity of the X-band radar (S3).
관측 자료를 수집하는 단계(S1)는 관측 자료 수집 모듈로 엑스밴드 레이더 반사도 보정을 위해 관측 자료를 수집한다. 여기서, 관측 자료를 수집하는 단계(S1)에서 수집하는 관측 자료는 마이크로파 링크의 송신 전력과 수신 전력, 마이크로파 링크 주변에 설치된 기상청의 자동기상관측장비(Automatic Weather System, AWS)의 강우 감지기에 기록된 강우 시간과 우량계에 기록된 누적 강우량과 강우 시간 자료, 엑스밴드 레이더의 반사도 및 에스밴드 레이더의 반사도를 포함한다.In the step of collecting observation data (S1), the observation data collection module collects observation data for X-band radar reflectivity correction. Here, the observation data collected in the step of collecting the observation data (S1) is the transmission power and reception power of the microwave link, recorded in the rainfall detector of the Automatic Weather System (AWS) of the Korea Meteorological Administration installed around the microwave link. It includes the rainfall time and cumulative rainfall and rainfall time data recorded in the rain gage, the reflectivity of the X-band radar, and the reflectivity of the S-band radar.
마이크로파 감쇠 자료를 기반으로 강우 유무를 탐지하는 단계(S2)는 강우 유무 탐지 모듈이 마이크로파 감쇠 자료를 기반으로 강우 유무를 탐지한다.In step S2 of detecting whether or not there is rain based on the microwave attenuation data, the rain presence detection module detects whether or not there is rain based on the microwave attenuation data.
엑스밴드 레이더의 반사도를 보정하는 단계(S3)는 비용함수를 판단하여 감쇠 보정 계수를 산출하고, 레이더 반사도 와 감쇠 보정 계수 , b를 적용하여 계산한 게이트(r)별 엑스밴드 레이더 감쇠 크기로 감쇠 보정한 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정한다. 이러한 엑스밴드 레이더의 반사도를 보정하는 단계(S3)는 비용함수를 판단하는 단계(S3-1)와, 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정하는 단계(S3-2)를 포함한다.In the step of correcting the reflectivity of the X-band radar (S3), the cost function is determined to calculate the attenuation correction coefficient, and the radar reflectivity and the attenuation correction factor X-band radar attenuation for each gate (r) calculated by applying , b Attenuation-corrected radar reflectivity with Calculate and correct the X-band radar reflectivity. The step of correcting the reflectivity of the X-band radar (S3) includes the step of determining the cost function (S3-1), the radar reflectivity Compensating the X-band radar reflectivity by calculating (S3-2).
비용함수를 판단하는 단계(S3-1)는 관측 자료를 수집하는 단계(S1)에서 수집된 관측 자료를 기반으로 비용함수를 연산하고, 연산된 비용함수가 최소인지 판단한다. 이를 위해서, 비용함수를 판단하는 단계(S3-1)는 엑스밴드 레이더 게이트 별 감쇠 크기를 연산하는 단계(S3-1-1)와, 링크 특성을 이용하여 강우에 의해 야기된 마이크로파 감쇠를 산출하는 단계(S3-1-2), 마이크로파 감쇠 크기를 연산하는 단계(S3-1-3), 엑스밴드 레이더 평균 감쇠 크기를 산출하는 단계(S3-1-4), 및 비용함수를 연산하는 단계(S3-1-5)를 포함한다.In the step of determining the cost function (S3-1), the cost function is calculated based on the observation data collected in the step of collecting observation data (S1), and it is determined whether the calculated cost function is minimum. To this end, the step of determining the cost function (S3-1) includes the step of calculating the attenuation size for each X-band radar gate (S3-1-1) and calculating the microwave attenuation caused by the rain using the link characteristics. Step (S3-1-2), calculating microwave attenuation (S3-1-3), calculating X-band radar average attenuation (S3-1-4), and calculating a cost function ( S3-1-5).
엑스밴드 레이더 게이트 별 감쇠 크기를 연산하는 단계(S3-1-1)는 엑스밴드 레이더 게이트 별 감쇠 크기 연산 모듈이 엑스밴드 레이더 게이트(r) 별 감쇠 크기()를 연산한다. 여기서, 엑스밴드 레이더 게이트(r) 별 감쇠 크기()는 전술된 수학식 7과 같다.In the step of calculating the attenuation size for each X-band radar gate (S3-1-1), the attenuation size calculation module for each X-band radar gate determines the attenuation size for each X-band radar gate r ( ) is computed. Here, the attenuation size for each X-band radar gate (r) ( ) is the same as Equation 7 described above.
링크 특성을 이용하여 강우에 의해 야기된 마이크로파 감쇠를 산출하는 단계(S3-1-2)는 마이크로파 감쇠 산출 모듈이 마이크로파 링크의 송신 전력과 수신 전력의 차이를 이용하여 마이크로파 감쇠(rainfall-caused attenuation, )를 산출한다.In step S3-1-2 of calculating the microwave attenuation caused by rainfall using the link characteristics, the microwave attenuation calculation module uses the difference between the transmission power and the reception power of the microwave link to calculate the microwave attenuation (rainfall-caused attenuation, ) is calculated.
마이크로파 감쇠 크기를 연산하는 단계(S3-1-3)는 마이크로파 감쇠 크기 산출 모듈이 레이더 주파수로 변환된 마이크로파 감쇠 크기()를 산출한다. 이는 전술된 수학식 9와 같이 구할 수 있다.In the step of calculating the microwave attenuation size (S3-1-3), the microwave attenuation size calculation module converts the microwave attenuation size to the radar frequency ( ) is calculated. This can be obtained as in Equation 9 above.
엑스밴드 레이더 평균 감쇠 크기를 산출하는 단계(S3-1-4)는 엑스밴드 레이더 평균 감쇠 크기 산출 모듈이 전술된 수학식 8과 같이 마이크로파 링크 경로 길이에 해당하는 엑스밴드 레이더 평균 감쇠 크기()를 계산한다.In the step of calculating the X-band radar average attenuation size (S3-1-4), the X-band radar average attenuation size calculation module corresponds to the microwave link path length as shown in
비용함수를 연산하는 단계(S3-1-4)는 비용함수 산출 모듈이 엑스밴드 레이더의 평균 감쇠 크기()와 변환된 마이크로파 감쇠 크기()간의 차이를 계산하여 전술된 수학식 10과 같이 비용함수()를 추정한다.In the step of calculating the cost function (S3-1-4), the average attenuation size of the X-band radar ( ) and the converted microwave attenuation magnitude ( ) by calculating the difference between the cost function (as shown in
레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정하는 단계(S3-2)는 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정한다. 이를 위해서, 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정하는 단계(S3-2)는 감쇠 보정 계수를 연산하는 단계(S3-2-1)와, 감쇠 보정 계수를 선택하는 단계(S3-2-2), 및 반사도를 보정하는 단계(S3-2-3)를 포함한다.radar reflectivity The step of correcting the X-band radar reflectivity by calculating (S3-2) is the radar reflectivity Calculate and correct the X-band radar reflectivity. For this purpose, the radar reflectivity Compensating the X-band radar reflectivity by calculating (S3-2) includes calculating an attenuation correction coefficient (S3-2-1), selecting an attenuation correction coefficient (S3-2-2), and reflectivity. and a step of correcting (S3-2-3).
감쇠 보정 계수를 연산하는 단계(S3-2-1)는 마이크로파 링크 경로 길이에 해당하는 강우에 의한 7개의 엑스밴드 레이더의 평균 감쇠 크기를 7개의 엑스밴드 레이더 주파수로 변환한 마이크로파 링크 감쇠와 각각의 비용함수가 최소가 될 때의 감쇠 보정 계수를 산출한다. 물론, 본 실시예는 관측 자료 수집 모듈에 7개의 마이크로파 링크의 특성이 입력되는 것이므로 7개의 엑스밴드 레이더의 평균 감쇠 크기를 예시하였으며, 관측 자료 수집 모듈에 입력되는 마이크로파 링크의 특성 개수에 따라 이는 달라질 수 있다.The step of calculating the attenuation correction coefficient (S3-2-1) is the microwave link attenuation obtained by converting the average attenuation of seven X-band radars by rainfall corresponding to the microwave link path length into seven X-band radar frequencies and each Calculate the damping correction coefficient when the cost function becomes the minimum. Of course, in this embodiment, since the characteristics of seven microwave links are input to the observation data collection module, the average attenuation of seven X-band radars is exemplified, and this will vary depending on the number of characteristics of microwave links input to the observation data collection module. can
감쇠 보정 계수를 선택하는 단계(S3-2-2)는 감쇠 보정 계수를 연산하는 단계(S3-2-1)에서 산출된 감쇠 보정 계수 중, 시공간적으로 변동성이 적은 층상운 강우 사례에서는 7개의 감쇠 보정 계수들의 평균값을 선택하고, 대류운 사례에서는 7개의 감쇠 보정 계수들의 최대값을 최적의 감쇠 보정 계수로 선택한다.In the step of selecting the attenuation correction coefficient (S3-2-2), among the attenuation correction coefficients calculated in the step of calculating the attenuation correction coefficient (S3-2-1), in the case of stratiform rainfall with little temporal and spatial variability, seven attenuations are selected. The average value of the correction coefficients is selected, and in the case of convective clouds, the maximum value of the seven attenuation correction coefficients is selected as the optimal attenuation correction coefficient.
반사도를 보정하는 단계(S3-2-3)는 엑스밴드 레이더 반사도 보정 모듈이 레이더 반사도 와 감쇠 보정 계수 , b를 적용하여 계산한 게이트(r)별 엑스밴드 레이더 감쇠 크기를 이용하여, 전술된 수학식 11과 같이 감쇠 보정한 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정한다.In the step of correcting the reflectivity (S3-2-3), the X-band radar reflectivity correction module adjusts the radar reflectivity and the attenuation correction factor X-band radar attenuation for each gate (r) calculated by applying , b Using
이상과 같이 본 발명에 대해서 예시한 도면을 참조로 하여 설명하였으나, 본 명세서에 개시된 실시 예와 도면에 의해 본 발명이 한정되는 것은 아니며, 본 발명의 기술사상의 범위 내에서 통상의 기술자에 의해 다양한 변형이 이루어질 수 있음은 자명하다. 아울러 앞서 본 발명의 실시 예를 설명하면서 본 발명의 구성에 따른 작용 효과를 명시적으로 기재하여 설명하지 않았을 지라도, 해당 구성에 의해 예측 가능한 효과 또한 인정되어야 함은 당연하다. 즉, 본 발명은 전술된 의료요청 상황 이외의 다양한 상황, 예를 들어, 지역 탐색, 맵 서비스, 특정 지역의 베이스 자료 구축 등에 적용될 수도 있다. 또한, 전술된 제어 방법은 프로그램으로 구현되어 저장 매체에 저장될 수 있다.As described above, the present invention has been described with reference to the drawings illustrated, but the present invention is not limited by the embodiments and drawings disclosed in this specification, and various modifications are made by those skilled in the art within the scope of the technical idea of the present invention. It is obvious that variations can be made. In addition, although the operational effects according to the configuration of the present invention have not been explicitly described and described while describing the embodiments of the present invention, it is natural that the effects predictable by the corresponding configuration should also be recognized. That is, the present invention may be applied to various situations other than the medical request situation described above, for example, area search, map service, base data construction of a specific area, and the like. In addition, the above-described control method may be implemented as a program and stored in a storage medium.
100: 관측 자료 수집 모듈
200: 마이크로파 감쇠 산출 모듈
300: 강우 유무 탐지 모듈
310: 창 크기 산출 모듈
320: 평균 마이크로파 감쇠 산출 모듈
330: 마이크로파 감쇠의 검정 통계량 산출 모듈
340: 임계 값 산출 모듈
350: 1차 강우 유무 탐지 모듈
360: 마이크로파 기준 감쇠 산출 모듈
370: 2차 강우 유무 탐지 모듈
400: 반사도 보정 모듈
410: 마이크로파 감쇠 결정 모듈
420: 엑스밴드 레이더 게이트 별 감쇠 크기 연산 모듈
430: 엑스밴드 레이더 평균 감쇠 크기 산출 모듈
440: 마이크로파 감쇠 크기 산출 모듈
450: 비용함수 산출 모듈
460: 엑스밴드 레이더 반사도 보정 모듈100: observation data collection module
200: microwave attenuation calculation module
300: rain presence detection module
310: window size calculation module
320 Average microwave attenuation calculation module
330: test statistic calculation module of microwave attenuation
340: threshold value calculation module
350: primary rain presence detection module
360: microwave reference attenuation calculation module
370: secondary rain presence detection module
400: reflectivity correction module
410: microwave attenuation determination module
420: attenuation size calculation module for each X-band radar gate
430: X-band radar average attenuation magnitude calculation module
440: microwave attenuation size calculation module
450: cost function calculation module
460: X-band radar reflectivity correction module
Claims (4)
상기 마이크로파 링크의 송신 전력과 수신 전력간의 차이를 이용하여 마이크로파 감쇠를 산출하는 마이크로파 감쇠 산출 모듈,
상기 마이크로파 감쇠의 자기상관계수를 기반으로 강우 유무를 탐지하는 강우 유무 탐지 모듈, 및
상기 강우 유무 탐지 모듈에서 강우 탐지 시 엑스밴드 레이더 게이트(r) 별 감쇠 크기()와, 엑스밴드 레이더의 레이더 반사도()를 이용하여 산출된 감쇠 보정된 레이더 반사도()로 엑스밴드 레이더 반사도를 보정하는 엑스밴드 레이더 반사도 보정 모듈을 가지는 반사도 보정 모듈을 포함하는 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치.
The transmission power and reception power of the microwave link are collected, and recorded at regular intervals in the rain detector and rain gauge of at least one Automatic Weather System (AWS) installed around the microwave link. An observation data collection module for collecting cumulative rainfall data and rainfall time data;
A microwave attenuation calculation module for calculating microwave attenuation using a difference between transmit power and received power of the microwave link;
A rain presence detection module for detecting whether or not there is rain based on the autocorrelation coefficient of the microwave attenuation, and
Attenuation size for each X-band radar gate (r) when rain is detected by the rain presence detection module ( ) and the radar reflectivity of the X-band radar ( ) attenuation-corrected radar reflectivity calculated using ) X-band dual polarization radar reflectivity correction device using multi-terrestrial communication base station microwaves including a reflectivity correction module having an X-band radar reflectivity correction module for correcting X-band radar reflectivity with ).
상기 감쇠 보정된 레이더 반사도()는,
이며,
상기 은 엑스밴드 레이더 감쇠 크기로서,
이고,
,
,
상기 , b는 엑스밴드 레이더의 감쇠 보정 계수인 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 장치.
According to claim 1,
The attenuation corrected radar reflectivity ( )Is,
is,
remind is the magnitude of the X-band radar attenuation,
ego,
,
,
remind X-band dual polarization radar reflectivity correction device using multi-terrestrial communication base station microwaves, where b is the attenuation correction coefficient of X-band radar.
상기 마이크로파 링크의 송신 전력과 수신 전력간의 차이를 이용하여 마이크로파 감쇠 산출 모듈이 마이크로파 감쇠를 산출하는 단계,
상기 마이크로파 감쇠의 자기상관계수를 기반으로 강우 유무 탐지 모듈이 강우 유무를 탐지하는 단계, 및
반사도 보정 모듈이 비용함수를 판단하여 감쇠 보정 계수를 산출하고, 레이더 반사도 와 감쇠 보정 계수 , b를 적용하여 계산한 게이트(r)별 엑스밴드 레이더 감쇠 크기로 감쇠 보정한 레이더 반사도를 산출하여 엑스밴드 레이더 반사도를 보정하는 단계를 포함하는 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법.
An observation data collection module collects transmission power and reception power of a microwave link, and transmits power to a rain detector and a rain gauge of at least one Automatic Weather System (AWS) installed around the microwave link. Collecting accumulated rainfall data and rainfall time data recorded at regular intervals;
Calculating, by a microwave attenuation calculation module, microwave attenuation using a difference between transmit power and received power of the microwave link;
Detecting presence or absence of rain by a rain presence detection module based on the autocorrelation coefficient of the microwave attenuation; and
The reflectivity correction module determines the cost function to calculate the attenuation correction coefficient, and the radar reflectivity and the attenuation correction factor X-band radar attenuation for each gate (r) calculated by applying , b Attenuation-corrected radar reflectivity with A method for correcting X-band dual polarization radar reflectivity using multi-terrestrial communication base station microwaves, comprising calculating and correcting X-band radar reflectivity.
상기 감쇠 보정된 레이더 반사도()는,
이며,
상기 은 엑스밴드 레이더 감쇠 크기로서,
이고,
,
,
상기 , b는 엑스밴드 레이더의 감쇠 보정 계수인 다중 지상 통신 기지국 마이크로파를 활용한 엑스밴드 이중 편파 레이더 반사도 보정 방법.According to claim 3,
The attenuation corrected radar reflectivity ( )Is,
is,
remind is the magnitude of the X-band radar attenuation,
ego,
,
,
remind X-band dual polarization radar reflectivity correction method using multi-terrestrial communication base station microwaves, where b is the attenuation correction coefficient of X-band radar.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020210087075A KR102632006B1 (en) | 2021-07-02 | 2021-07-02 | Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020210087075A KR102632006B1 (en) | 2021-07-02 | 2021-07-02 | Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20230006193A true KR20230006193A (en) | 2023-01-10 |
KR102632006B1 KR102632006B1 (en) | 2024-02-01 |
Family
ID=84893950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020210087075A KR102632006B1 (en) | 2021-07-02 | 2021-07-02 | Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR102632006B1 (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19990076841A (en) * | 1995-12-26 | 1999-10-25 | 게네비브 부 탄 | METHOD FOR DETERMINING RATE YEAR RATE USING VAPOR RADAR AND METHODIC RADAR FOR IMPLEMENTING THE METHOD |
JP2011047704A (en) * | 2009-08-25 | 2011-03-10 | Toshiba Corp | Weather radar system and rainfall rate calculation method thereof, and program |
KR20120125900A (en) * | 2011-05-09 | 2012-11-19 | 경북대학교 산학협력단 | Method of classify meteorological and non-meteorological echoes using dual polarization radars |
KR101291980B1 (en) | 2012-12-20 | 2013-08-09 | 경북대학교 산학협력단 | Method for making total quality index for radar reflectivity measurement |
KR101431707B1 (en) * | 2013-11-26 | 2014-09-22 | 한국건설기술연구원 | method of classification and quantification using data of X band dual polarization radar |
KR101915665B1 (en) * | 2018-01-22 | 2018-11-06 | 부경대학교 산학협력단 | Monitoring system for rainy and monitoring method using same by microwave attenuation character due to rainfall |
-
2021
- 2021-07-02 KR KR1020210087075A patent/KR102632006B1/en active IP Right Grant
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR19990076841A (en) * | 1995-12-26 | 1999-10-25 | 게네비브 부 탄 | METHOD FOR DETERMINING RATE YEAR RATE USING VAPOR RADAR AND METHODIC RADAR FOR IMPLEMENTING THE METHOD |
JP2011047704A (en) * | 2009-08-25 | 2011-03-10 | Toshiba Corp | Weather radar system and rainfall rate calculation method thereof, and program |
KR20120125900A (en) * | 2011-05-09 | 2012-11-19 | 경북대학교 산학협력단 | Method of classify meteorological and non-meteorological echoes using dual polarization radars |
KR101291980B1 (en) | 2012-12-20 | 2013-08-09 | 경북대학교 산학협력단 | Method for making total quality index for radar reflectivity measurement |
KR101431707B1 (en) * | 2013-11-26 | 2014-09-22 | 한국건설기술연구원 | method of classification and quantification using data of X band dual polarization radar |
KR101915665B1 (en) * | 2018-01-22 | 2018-11-06 | 부경대학교 산학협력단 | Monitoring system for rainy and monitoring method using same by microwave attenuation character due to rainfall |
Non-Patent Citations (2)
Title |
---|
강미영 外 3인, '연구용 X-band 이중 편파레이더를 이용한 강수량 추정 사례분석', 2011년도 한국기상학회 가을학술대회 논문집.* * |
강미영 外 6인, '이중 편파레이더를 이용한 정량적 강우량 추정', 2005년도 한국기상학회 가을학술대회 논문집.* * |
Also Published As
Publication number | Publication date |
---|---|
KR102632006B1 (en) | 2024-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shepherd | Radio wave loss deviation and shadow loss at 900 MHz | |
US8988275B2 (en) | Low energy radar system | |
Seidel et al. | Propagation measurements at 28 GHz to investigate the performance of local multipoint distribution service (LMDS) | |
KR101915665B1 (en) | Monitoring system for rainy and monitoring method using same by microwave attenuation character due to rainfall | |
CN110596710A (en) | Fire monitoring and decision-making system and method based on X-band dual-polarization weather radar network | |
CN109936855B (en) | Method, device, equipment and medium for inverting precipitation distribution by using base station signals | |
Vaccarono et al. | An integrated approach to monitoring the calibration stability of operational dual-polarization radars | |
US20210382170A1 (en) | Apparatus and method for removing noise for observation information of weather radar | |
Borga et al. | Rainfall estimation by weather radar | |
KR102632006B1 (en) | Device for attenuation correction of x-band dual polar radar reflectivity using adjacent multiple terrestrial microwave links and method therefor | |
Barrios et al. | Low altitude propagation effects—A validation study of the advanced propagation model (APM) for mobile radio applications | |
KR102632005B1 (en) | Device for attenuation correction of x-band single polarization radar reflectivity using adjacent multiple terrestrial microwave links and method therefor | |
Thurai et al. | Rain height variability in the tropics | |
Bringi et al. | Measurements and inferences of raindrop canting angles | |
KR101810405B1 (en) | System for detecting fog and fine particle using multiple frequencies | |
Rakshit et al. | Effective path length for estimating rain attenuation over an earth-space path using features of stratiform and convective precipitation at a tropical location | |
Sinatra et al. | Calibration of spatial rain scanner using rainfall depth of rain gauges | |
Adirosi et al. | NEFOCAST project for real-time precipitation estimation from Ku satellite links: Preliminary results of the validation field campaign | |
Toledo et al. | Absolute Calibration method for FMCW Cloud Radars | |
Bharadwaj et al. | Calibration system for ARM radars | |
Yu et al. | Application of X-band Polarimetric Phased-array Radars in Quantitative Precipitation Estimation | |
Thurai et al. | Precipitation scatter measurements from a transhorizon experiment at 11.2 GHz | |
Norland | Temporal variation of the refractive index in coastal waters | |
Lakanchanh et al. | Rainfall rate and rain attenuation in Ku-band satellite signal in Thailand and Laos | |
Ghaddar et al. | Experimental analysis of human body effects on NLOS 60 GHz propagation channel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant |