KR101469860B1 - 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법 - Google Patents

기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법 Download PDF

Info

Publication number
KR101469860B1
KR101469860B1 KR20130119358A KR20130119358A KR101469860B1 KR 101469860 B1 KR101469860 B1 KR 101469860B1 KR 20130119358 A KR20130119358 A KR 20130119358A KR 20130119358 A KR20130119358 A KR 20130119358A KR 101469860 B1 KR101469860 B1 KR 101469860B1
Authority
KR
South Korea
Prior art keywords
rainfall
precipitation
peak
data
intensity
Prior art date
Application number
KR20130119358A
Other languages
English (en)
Inventor
이동인
김지현
김주연
Original Assignee
부경대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 부경대학교 산학협력단 filed Critical 부경대학교 산학협력단
Priority to KR20130119358A priority Critical patent/KR101469860B1/ko
Application granted granted Critical
Publication of KR101469860B1 publication Critical patent/KR101469860B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • 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
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • 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
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W2203/00—Real-time site-specific personalized weather information, e.g. nowcasting
    • 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

본 발명은 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법에 관한 것으로서, 보다 상세하게는 1분 단위의 기상 레이더 자료를 이용하여 수집된 각각의 강수샘플에 대한 강우피크개수, 피크에서의 강우강도, 누적강우량, 강수지속시간 등의 대류성 강우특성정보를 얻는 방법에 관한 것이다.
본 발명의 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법은, 기상레이더를 통해 레이더 강우강도를 산출하는 단계; 상기 강우강도에 대한 극 좌표계에서의 기상자료를 직교좌표계에 대한 기상자료로 변환하는 단계; Cressman 내삽법을 이용하여 각각의 지리학적 단위 격자마다의 합성 레이더 자료값을 산출하는 단계; 상기 지리학적 단위 격자마다의 각각의 강수샘플에 대한 누적강우량(Rac) 데이터를 생성하는 단계; 각각의 강수샘플에 대해서, 불필요한 노이즈 값 제거를 위해, 시간에 따른 강우강도를 FIR(Finite Impulse Response) 필터링하는 단계; 및각각의 강수샘플에 대해서, 상기 필터링된 시간에 따른 강우강도를 이용하여 강우피크개수(Npeak)를 포함한다.

Description

기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법{Method for gaining information about the convection precipitation in use of the weather radar}
본 발명은 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법에 관한 것으로서, 보다 상세하게는 1분 단위의 기상 레이더 자료를 이용하여 수집된 각각의 강수샘플에 대한 강우피크개수, 피크에서의 강우강도, 누적강우량, 강수지속시간 등의 대류성 강우특성정보를 얻는 방법에 관한 것이다.
기상 레이더는 전파를 이용하여 강수의 입자와 크기를 측정하는 원격 관측 장비로써, 시공간적 해상도를 가지고 유효반경 240 km 이내 10분 정도로 넓은 영역의 강수량을 산출하는 기상 장비 중 하나이다. 따라서 중규모 강수시스템인 태풍이나 장마, 홍수 등 넓은 범위에 대한 전체적인 기상상태를 실시간 파악할 수 있고, 비구름의 발달 구조 상태를 관측할 수 있다.
그러나 지금까지 기상 레이더를 이용한 대류성 강우에 대한 연구는 5분에서 10분 정도의 3차원 볼륨데이터를 이용하여 강수량의 변화에 대해 분석해오고 있으나, 이는 현재 빈번히 발생하고 있는 국지적 대류성 강우에 대한 시간적 주기 분석에 대한 해석이 어렵다. 이를 해결하기 위해서는 고해상도 자료인 강우 레이더가 필수적이다.
X 밴드 레이더는 S 밴드나 C 밴드에 비해 파장이 짧고, 주파수가 높기 때문에 목표물을 식별하고 변별하기 위한 고해상도 이미지 레이더를 만들 수 있는 장점이 있다. 특히 이중편파인 경우 수직편파와 수평편파를 교대로 발사하여 목표물을 관측함으로써 좁은 지역에서 돌발적으로 발생하는 대류성 호우를 정확히 예측할 수 있다.
현재 우리나라 국립기상연구소에서는 고해상도 자료를 산출하는 X 밴드 이중편파 레이더를 연구목적으로 운용하고 있으나, 강수 보정에 대한 기술적인 부분이 미흡하여 현업에는 쓰이지 못하는 실정이다.
아울러, 국지성 강우인 대류성 강우에 정확한 예측을 위해서, 대류성 강우에 대한 특성 정보의 수집 또한 미흡한 실정이다.
한국공개특허공보 공개번호 10-20120040025(2012.04.26공개)(발명의 명칭: 선박 레이더를 이용하여 국지 강수를 추정하는 시스템 및 그 추정방법)
본 발명은 상기의 문제점을 해결하기 위한 것으로서, 본 발명의 목적은 일본 국토교통성이 운용하고 현업에 쓰이고 있는 X 밴드 이중 편파 레이더(X-band polarimetric radar rainfall data)를 이용하여 1분마다 관측되는 강우강도를 통해 각 지리학적 격자마다 시간에 따른 강우강도 변화를 파악하여 대류성 강우 특성 정보를 얻는 방법을 제공하는 것이다.
본 발명의 목적은 이상에서 언급한 목적으로 제한되지 않으며, 언급되지 않은 또다른 목적들은 아래의 기재로부터 명확하게 이해될 수 있을 것이다.
상기의 목적을 달성하기 위해서, 본 발명의 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법은, 기상레이더를 통해 레이더 강우강도를 산출하는 단계; 상기 강우강도에 대한 극 좌표계에서의 기상자료를 직교좌표계에 대한 기상자료로 변환하는 단계; Cressman 내삽법을 이용하여 각각의 지리학적 단위 격자마다의 합성 레이더 자료값을 산출하는 단계; 상기 지리학적 단위 격자마다의 각각의 강수샘플에 대한 누적강우량(Rac) 데이터를 생성하는 단계; 각각의 강수샘플에 대해서, 불필요한 노이즈 값 제거를 위해, 시간에 따른 강우강도를 FIR(Finite Impulse Response) 필터링하는 단계; 및 각각의 강수샘플에 대해서, 상기 필터링된 시간에 따른 강우강도를 이용하여 강우피크개수(Npeak)를 포함하는 대류성 강우의 강우특성 정보를 얻는 단계;를 포함한다.
바람직한 실시예에 따르면, 상기 기상레이더에서 얻는 강우강도는 1분 단위의 강우강도이며, 상기 산출된 강우강도는 비차등 위상 변위(KDP)에 따른 강우강도이다.
바람직하게는 상기 누적강우량(Rac) 데이터를 생성하는 단계에서, 강우강도가 0.5mmhr-1 이상인 경우의 강수샘플에 대해서만 누적강우량 데이터를 생성하며,
바람직하게는 상기 FIR 필터링은 해밍창 함수법을 이용하여 낮은 주파수 영역을 통과시키고, 높은 주파수 영역을 차단하는 저역통과 필터가 사용된다.
바람직하게는 상기 대류성 강우의 강우특성 정보는 강우피크의 수(Npeak), 피크 강우강도(Rpeak), 누적강우량(Rac), 강수 지속시간(Train), 강수시작에서 피크 강우강도까지의 시간(Tpeak)를 포함한다.
이상에서 상술한 바와 같이, 본 발명은 1분 단위의 기상 레이더를 이용하고, 수직좌표계로 전환, 내삽법, 필터링 등을 통해 아주 정확한 대류성 강우의 특성정보를 얻을 수 있으며, 강우피크 1인 경우만을 다시 추출하여 이를 바탕으로 매우 높은 상관관계를 얻을 수 있어 대류성 강우를 보다 정확히 예측할 수 있다는 장점이 있다.
도 1은 본 발명의 바람직한 실시예에 따른, 기상 레이더 자료를 이용한 강수 특성 분류 방법에 대한 순서도.
도 2는 본 발명의 일실시예에 따른, 두 레이더의 낮은 두 고도각을 이용한 1분 강우강도 합성 PPI SCAN 사진(2010년 9월 28일 11시 40분).
도 3은 본 발명의 일실시예에 따른, 실험에서의 레이더 위치 및 분석영역을 표시하는 지도.
도 4는 본 발명의 일실시예에 따른, 11시간 누적강우량에 대한 피크별로 분류하여 나타낸 사진.
도 5는 본 발명의 일실시예에 따른, 4시간 누적강우량과, 한 지점(+ 지점)상에서의 시간에 따른 강우강도 변화 그래프.
도 6은 본 발명의 일실시예에 따른, 한 지점(+)상 위·경도에 따른 시계열 강우강도 변화를 나타내는 사진.
도 7은 본 발명의 바람직한 실시예에 따른, 강수 샘플의 모식도.
도 8은 본 발명의 일실시예에 따른 Rac 에 대한 Rpeak, Train, and Tpeak 의 상관분석 및 회귀분석 그래프.
도 9는 도 8에서의 상관계수와 선형회귀분석 그래프.
이하에서는 본 발명의 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법에 대해서 보다 상세히 설명하기로 한다.
도 1은 본 발명의 바람직한 실시예에 따른, 기상 레이더 자료를 이용한 강수 특성을 얻는 방법에 대한 순서도로서, 먼저, 도 1을 기준으로 상세히 설명하기로 한다.
본 발명은 1분 강우강도 레이더 자료 생성을 위하여 ① X 밴드 레이더의 1분 강우강도 산출, ② 레이더 자료 격자 전환, ③ Cressman 내삽법을 이용한 합성 레이더 자료를 산출하고, ④ 강수에 대한 시간적 주기를 산출하기 위해 누적 강수량 데이터 생성, ⑤ 불필요한 노이즈 값 제거를 위한 FIR(Finite Impulse Response) 필터 설계, ⑥ 각 지리학적 격자마다 시간에 따른 강우강도 피크 분류(대류성 강우특성정보를 얻는 단계)로 구성되어 진다.
1. X 밴드 이중 편파 레이더의 1분 강우강도 산출
먼저, X 밴드 이중 편파 기상레이더를 통해 1분 단위의 레이더 강우강도를 산출한다.
X 밴드 이중 편파 레이더의 매개변수 중 차등위상변화(differential phase change, ΦDP)는 실제로 전파경로(propagation path) 상의 주어진 지점에서 수평편파의 위상변이(Φhh)와 수직편파의 위상변이(Φvv)와의 차이를 말하는 것으로, 일반적으로 빗방울은 편구(oblate spheroid) 형상으로 근사할 수 있으며, 수평축이 연직축(수직축)보다 더 길다. 따라서 레이더파의 전파경로(propagation path)가 연직보다는 수평방향으로 더 길다. 이로 인해서 수평편파의 위상과 연직편파의 차이가 생기며, Φhh > Φvv 가 되는 것이다. 따라서 이를 보정해 주기 위해서, 비차등 위상 변위(Specific differential phase shift, KDP)에 따른 강우강도를 추정하고 있다.
KDP는 단위 거리에 대한 차등위상변위로 정의되며 다음과 같이 표현된다.
Figure 112013090712906-pat00001
(식 1)
여기서, r1과 r2는 각각 레이더로 부터 레인지 bin 1과 bin 2까지의 거리(km)를 의미하며, KDP의 단위는 deg/km 이다. 차등위상변위는 관측거리에 비례하여 증가하고, 분모의 2는 레이더파가 주어진 레인지까지의 왕복경로에 따른 ΦDP의 증가를 고려한 것이다.
기상학적으로 에코가 가지는 일반적인 KDP값의 범위는 -1 ~ 6 deg/km 이고, KDP는 등방성입자에 대해 민감하지 않다. 그 이유는 이들 입자들은 거의 구형이므로 수평편파와 수직편파의 위상속도에 차이가 없기 때문이다. 따라서 KDP는 비와 우박이 혼재되어 있는 경우 강우량 추정을 하는데 매우 효율적으로 이용될 수 있다. 또한 KDP를 이용한 강우량 추정은 레이더 반사도 인자를 이용하는 경우보다 좀 더 효율적인 면이 있다. 즉, 부분적인 빔차폐가 발생할 경우, 레이더반사도 인자는 감소하지만 KDP는 전파가 산란되고 되돌아오는 한 부분적인 빔차폐의 영향을 받지 않는다는 점이다. 일반적으로 X 밴드 편파 레이더는 비의 감쇠가 크다는 단점이 있다. 그러나 역으로 취하면 빗방울이 편평하게 되었을 때의 연직과 수평 편파의 관측 값의 차이가 쉽게 나타난다는 것을 알 수 있다.
Maesaka et al. (2011)은 강우에 의한 감쇠의 영향을 줄이고 정확한 강우 강도 산출이 가능한 KDP 와 일반적으로 강우강도를 추정하는데 사용되는 ZH를 이용하여 강우강도를 산출하였다. 특히 추정된 KDP는 지상 우량계와 비교하였을시 과소평가되어 1.3배 보정하여 계산된다.
표 1은 X 밴드 편파 레이더에서 사용된 R(KDP)와 R(ZH)의 관계식을 나타낸 것이다. 즉, 표 1은 레이더 안테나 고도각(θ)을 고려한 강우강도 추정식이다.
약한 비의 경우 편광 파도 위상차 변화율 측정이 매우 작아지므로, 이때는 일반적으로 강우강도를 추정하는 R(ZH)로 계산한다. 이때 R(ZH)는 사용되는 a1, a2는 사용되는 레이더의 시스템에 따라 식이 달라진다. 30 dBz의 반사도를 기준으로 약한 비(< 30 dBz)의 경우 R(ZH)관계식을, 강한 비(≥ 30 dBz)의 경우 R(KDP)관계식으로 강우강도를 추정한다.
Equation Factor
Figure 112013090712906-pat00002
Figure 112013090712906-pat00003

Figure 112013090712906-pat00004

Figure 112013090712906-pat00005
Figure 112013090712906-pat00006
 a2 and b2 depend on individually radar system
2. 레이더 자료 격자 변환
다음으로, 지리학적 단위 격자로서의 변환을 위해서, 상기 강우강도에 대한 극 좌표계에서의 기상자료를 직교좌표계에 대한 기상자료로 변환한다.
기상레이더 관측 네트워크에서 관측된 반사도(dBZ)자료는 QC(Quality Check) 과정을 거쳐 3차원 반사도장 CAPPI(Constant Altitude Plan Position Indicator)으로 생산한다. 이때, 3차원 반사도장 자료가 수치 모형의 배경장을 보정하기 위해서는 일정한 간격을 가지는 격자의 변환이 필요하다. 따라서 격자의 크기가 일정치 않은 극 좌표계(Polar coordinate)형식의 3차원 반사도장은 내삽을 통해 일정한 간격을 가지는 직교 좌표계(Cartesian coordinate)로 변환한다.
본 발명에서는 '격자'란 지리학적 격자를 의미하는 것으로, 예컨데 평면상 250 m × 250 m의 지리학적 격자를 의미한다.
3. Cressman 내삽법을 이용한 합성 레이더 자료 산출
다음으로, 근가중치 보정(가까운 레이더에서 수집된 레이더 자료에 가중치를 부여하는 보정) 및 레이더 자료 중 단위 격자 상에서 빠진 정보 보완을 위해서, Cressman 내삽법을 이용하여 합성 레이더 자료를 산출한다.
레이더 자료는 12개의 고도각을 가진 3차원 볼륨 데이터를 1분당 3번의 PPI(Plan Position Indicator)를 스캔하지만, 이 중 가장 낮은 고도각(θ1, θ2)이 매 1회씩 관측되므로 분당(per minute) 낮은 고도각의 PPI 관측 데이터를 얻을 수 있다. 따라서 관측된 낮은 고도각(θ1, θ2)의 분당 자료로 두 레이더의 PPI 자료를 합성하여 강우강도 값을 산출한다.
사용된 X 밴드 이중 편파 레이더는 일본 동경에 위치한 사이타마(Saitama, STM)와 신요코하마(Shimyokohama, SYK)에 위치한 두 대의 레이더로써, STM는 1.4 °와 2.4 °, SYK는 1.7 °과 2.6 °의 가장 낮은 두 고도각을 이용하여 평균하여 내삽 및 합성되어진다. 이때 낮은 두 고도각을 사용하여 평균되어지기 때문에 1분 강우강도 데이터는 최소한 2분 정도의 계산 소요시간이 필요하다.
각 레이더 사이트에서의 추정된 강우강도는 내삽 되어지며, 해상도인 경도 45/4 °, 위도 30/4 ° (대략 250 m)에 의한 지리학적 격자에 합성되어진다. 수정된 Cressman 내삽법은 다음과 같다. Sampling Radius = Rs (m)는 레이더로부터의 거리(r)의 선형방정식에 의해 정의되었다.
Figure 112013090712906-pat00007
(식 2)
레이더의 특성과 지리적인 조건을 고려하여 a3은 0.013, b3은 150 m로 정의하였다.
Cressman 내삽법에서의 가중함수는 다음과 같이 정의 되었다.
Figure 112013090712906-pat00008
(식 3)
w h 은 수평 가중값, d는 격자와 bin 범위 사이의 거리, C h 는 w h 의 통제·지배계수 0.5이며, bin 범위가 샘플링 범위에서 벗어났을 때 (d< Rs), bin은 샘플링하지 않는다 (w h =0).
w a 은 고도 가중값, h는 지면보다 높은 고도값, C a 는 w a 의 통제·지배계수 20이며, H는 Maximum 고도 내삽 5000 m이며, bin 고도의 h가 H보다 높으면 (H<h), bin은 샘플링하지 않는다 (w a =0).
예를 들어, 도 2는 2010년 9월 28일 11시 40분에 대한 1분 레이더 강우강도 합성 이미지이다. 이 때 θ1인 STM 1.4 °와 SYK 1.7 °은 11시 38분 PPI이며, θ2인 STM 2.4 °와 SYK 2.6 °은 11시 39분 PPI이다. 이를 평균 합성하여 11시 40분 강우강도가 계산되어진다. (a)와 (b)는 STM 레이더, (c)와 (d)는 SYK 레이더이며, 회색영역은 전파 감쇄영역이며, 빨간색은 레이더로부터 20 km씩 거리를 표시하였다. (e)는 1분 레이더 강우강도의 합성 이미지를 나타내며 점선은 레이더의 반경(80 km)를 나타낸다.
4. 누적 강수 데이터 생성
다음으로, 지리학적 단위 격자마다의 각각의 강수샘플에 대한 누적강우량(Rac) 데이터를 생성한다.
도 3은 실시예에서의 레이더 위치 및 분석영역을 표시하는 지도이다.
X 밴드 이중 편파 레이더는 300 m × 250 m의 공간해상도와 1분 강우강도를 산출하는 고해상도 레이더 자료이며, X 밴드 레이더의 강우감쇠를 고려하여 50 ㎞ × 90 ㎞ 범위(실선) 내에서 발생하는 대류성 강우를 분석하였다[도 3 참조]. 따라서 본 발명의 실시예에서 생성되는 격자수는 약 6만개이다.
분석 영역에서 강수가 시작된 시점부터 소멸까지 1분 레이더 강우강도 자료를 통해 각 지리학적 격자마다의 누적된 강수 데이터를 생성한다. 하나의 격자에 대한 강우강도 임계값은 0.5 mmhr-1로 하여 시간변화에 따른 강우강도 값을 저장한다.
5. FIR 필터 설계
다음으로, 각각의 강수샘플에 대해서, 불필요한 노이즈 값 제거를 위해, 시간에 따른 강우강도를 FIR(Finite Impulse Response) 필터링한다.
X 밴드 이중 편파 레이더의 강우강도는 불필요한 단주기 성분이 포함되어 있어 이를 제거하고자 유한 임펄스 응답 필터인 FIR 필터를 사용하여 단주기 성분을 제거한다. 따라서 해밍창 함수법을 이용하여 낮은 주파수 영역을 통과시키고, 높은 주파수 영역을 차단하는 저역통과 필터를 사용하며, 필터 길이 31, 차단 주파수는 1/180 Hz 로 강우강도의 단주기 성분을 제거한다.
하나의 강수 샘플은 5분 이상의 강수 지속시간을 보였을 때로 정의한다.
도 4는 11시간 누적강우량에 대한 피크별로 분류하여 나타낸 사진이다.
즉, 2010년 7월 5일 12시부터 22시까지 저장된 11시간 누적강우량(accumulated raninfall)과, 각 격자마다 FIR 필터를 통해 얻어진 대류성 강우의 분당 강우강도 변화를 피크가 하나(single peak)일 때의 강우 샘플(cell)의 수(number of single-peak cell)와, 두 개 이상(multi peak)일 때의 강우 샘플의 수(number of multi-peak cell)와, 그리고 피크가 하나 이상일 때의 총 강우 샘플의 수(즉, 피크가 하나일 때의 샘플의 수와 피크가 두개 이상일 때의 샘플의 수의 합(number of total cell))에 대해서 지도상에 표시하였다.
도 4를 통해 하나의 대류성 강우가 발생하였을 때 각 격자마다의 대류성 피크 분포를 파악할 수 있다.
도 5는 본 발명의 일실시예에 따른, 4시간 누적강우량과, 한 지점(+ 지점)상에서의 시간에 따른 강우강도 변화 그래프이다. 즉, 도 5는 도 2에서의 11시간 누적강우량 중 대류성 강우가 활발히 발생했던 17시 30분부터 21시 30분까지의 4시간 누적강우량과, 한 지점(+로 표시된 지점) 상에서의 시간에 따른 강우강도 변화를 나타낸 것이다.
시계열 그래프에서 검은선은 초기 레이더 자료값, 초록선은 FIR 필터에 의한 자료값, 노란선은 피크개수를 나타내고 있다. 특히 빨간 테두리에서의 강우강도 변화는 FIR 필터를 통해 불필요한 단주기 성분이 제거되는 것을 보여준다. 또한 도 6과 같이, 위도와 경도에 따른 시계열 분석도 가능하다.
6. 피크 분류 - 대류성 강우특성정보 획득
다음으로, 각각의 강수샘플에 대해서, 상기 필터링된 시간에 따른 강우강도를 이용하여 강우피크개수(Npeak)를 포함하는 대류성 강우의 강우특성 정보를 얻는다.
강수는 보통 대류성과 층류형으로 나누어지는데, 특히 자연재해 및 인간에게 직접적인 영향을 미치는 강한 강수가 대류성 강우이다. Short et al. (1993)과 Rosenfeld and Amitai (1995)는 대류성과 층류형 강수의 분류 기준을 10 mmhr-1의 임계값으로 정의하였다.
따라서 본 기술에서는 하나의 연속적인 강수가 발생하였을 때 최소한 최대 강우강도가 10 mmhr-1 이상이었을 때 하나의 대류성 강우 샘플로 정의한다.
최대 강우강도 값에 의해 분류된 대류성 강우 샘플은 피크의 개수에 따라 각각 분류된다. 피크가 하나 (single peak)일 때의 강우샘플의 수와, 피크가 두 개 이상(multi peak) 일 때로 분류하여 누적된 강수 분포와 이에 대한 피크 분포에 대해 비교할 수 있다.
하나의 강수 샘플은 도 7과 같이 5가지 변수로 강수의 특성을 파악할 수 있다.
도 7은 본 발명의 바람직한 실시예에 따른, 강수 샘플의 모식도이다.
도 7에서, 각각의 변수는 피크개수(Npeak), 피크강우강도(Rpeak), 누적강우량(Rac), 강수지속시간(Train), 강수시작에서 피크 강우강도까지의 시간(Tpeak)의 5가지 변수를 나타낸다. 검은선은 관측된 강우강도(mmhr-1), 빨간선은 FIR 필터를 통해 얻어진 강우강도(mmhr-1)이며, 초록선은 피크 수를 나타낸다.
또한 ASCII 파일로 각 지리학적 격자마다 발생된 대류성 강우 샘플에 대해서, 피크개수(Npeak), 강수 시작시간, 소멸(소강)시간, 경도(°E), 위도(°N), 강수 지속시간(Train), 강수시작에서 피크 강우강도까지의 시간(Tpeak), 누적강우량(Rac), 최대 강우강도(Rmax), 각 피크에서의 강우강도(Rpeak) 등 다양한 변수들을 저장하여 강수 피크별로 대류성 강우의 특성을 분류할 수 있다.
따라서, 각 지리학적 격자에 대해서 강우강도 임계값인 10 mmhr-1을 통해 대류성 강우에 대한 강우특성을 분류하였으며, 피크 개수에 따른 각 변수 값들이 ASCII 파일로 저장된다.
표 2는 도 3에서의 모든 지리학적 격자(약 6만개)에 대한 각각의 강수 샘플에 대한 강우특성을 분류한 표로서, 각각의 강수샘플에 대한 피크 개수에 따른 강우특성 정보가 저장된다.
① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨ ⑩
1 3 201007050401 201007050417 139.7422 35.97187 17 5 1.371 8.6
2 1 201007051041 201007051052 139.3328 35.44479 12 6 1.153 19.8
3 3 201007221441 201007221451 139.6078 35.63021 11 9 1.14 12.2
4 10 201007310416 201007310532 139.3297 35.79688 77 11 8.691 27.3
5 5 201008160831 201008160903 139.7734 35.87188 33 17 3.476 38.7
6 2 201008270853 201008270908 139.6266 36.04271 16 5 1.373 14.9
7 1 201009280211 201009280219 139.5641 35.33229 9 3 0.216 3.8
·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·		 ·
·
·
345,948 2 201009280357 201009280405 139.9078 36.11562 9 4 0.218 3.4
상기 표 2에서, ① 강수샘플 ID(강수 샘플 번호), ② 피크개수(Npeak), ③ 강수 시작시간, ④ 강수 소강시간, ⑤ 경도(°E), ⑥ 위도(°N), ⑦ 강수 지속시간(Train), ⑧ 강수시작에서 피크 강우강도까지의 시간(Tpeak), ⑨ 누적강우량(Rac), ⑩ 피크강우강도(Rpeak)을 나타낸다.
8. 강우피크의 개수가 1인 경우 추출
다음으로, 상기 강우특성정보 중 강우피크의 개수(Npeak)가 1개인 경우만의 강우특성정보를 추출한다. 즉, 상기와 같은 방법으로 얻어진 강수 샘플들에 대한 대류성 강우특성 정보 중에 Npeak가 단일 피크 (Single peak)일 때의 대류성 강우만을 추출한다.
단일 피크(Npeak=1)일 때의 강수샘플 중에서, 피크 강우강도 Rpeak의 최대치는 128.1 mmhr- 1였고, 누적강우량 Rac의 최대치는 40.4 mm였으며, 강수 지속시간 Train의 최대치는 42 min., 강수시작에서 피크 강우강도까지의 시간 Tpeak 의 최대치는 21 min.으로 나타났다.
강우피크의 개수가 1일 경우만을 추출하는 이유는 하기에서와 같이, 대류성 강우특성정보들의 상관성이 매우 높이 나타나기 때문이다.
9. 상관관계 도출
다음으로, 추출된 강우피크의 개수가 1인 경우의 강우특성정보를 이용하여, 누적강우량과 다른 강우특성정보와의 상관 관계를 도출한다.
단일피크(Npeak=1) 일 때의 강수 샘플 수는 전체 강수 샘플 수(345,948개) 중에서 약 25%인 88,036개였다. 즉, 88,036개의 단일피크(Npeak=1) 강수샘플을 이용하여, 피크 강우강도 Rpeak, 누적강우량 Rac, 강수 지속시간 Train, 강수시작에서 피크 강우강도까지의 시간 Tpeak 에 대한 상관관계 분석을 실시하였다.
단일 피크(Npeak=1)일 때의 강수 샘플에 대한 통계분석은 피어슨 상관관계 분석법을 이용하여 실시하였다. 표 3은 피어슨 상관관계 분석 결과를 나타낸다.
Rpeak Rac Train Tpeak
Rpeak - 0.93 0.62 0.49
Rac 0.93 - 0.67 0.52
Train 0.62 0.67 - 0.76
Tpeak 0.49 0.52 0.76 -
상기 표 3의 피어슨 상관관계 분석결과, Rac - Rpeak는 0.9 이상의 매우 높은 상관 값을 보였으며, 또한 Tpeak - Rrain도 0.75 이상의 높은 상관값을 나타내었다.
여기서, 주목할 점은 4개의 변수(Rpeak, Rac, Train, Tpeak) 중 Rac가 다른 세 변수(Rpeak, Rac, Train, Tpeak)와의 상관성이 다른 변수들간의 상관성보다 높다는 것을 확인할 수 있었다. 즉, 누적강우량 Rac는 피크 강우강도 Rpeak과의 상관값은 0.93, 강수 지속시간 Train과의 상관값은 0.67, 강수시작에서 피크 강우강도까지의 시간 Tpeak 와의 상관값은 0.52였음으로, 모두 0.5이상의 높은 상관관계를 나타내었다.
따라서, 누적강우량 Rac를 통해 다른 변수들(Rpeak, Train, Tpeak)을 예측가능하다는 사실을 확인할 수 있었다.
한편, 일반적인 대류성 강우와 강한 대류성 강우의 특성을 보다 구체적으로 살펴보기 위해, 표 2의 자료에 대해서 피크 강우강도를 10 mmhr-1과 50 mmhr-1의 이상으로 재분류하여 상관분석 및 회귀분석을 실시하였다
도 8은 본 발명의 일실시예에 따른 Rac 에 대한 Rpeak, Train, and Tpeak 의 상관분석 및 회귀분석 그래프이다.
도 8에서는 (a), (b), (c)는 일반적인 대류성 강우(피크강우강도 Rpeak ≥10 mmhr-1)에 대한 그래프로서, (a)는 누적강우량 Rac와 Rpeak 과의 상관분석, 회귀분석 그래프, (b)는 누적강우량 Rac와 Train과의 상관분석, 회귀분석 그래프, (c)는 누적강우량 Rac와 Tpeak 과의 상관분석, 회귀분석 그래프이며, (d), (e), (f)는 강한 대류성 강우(피크강우강도 Rpeak ≥ 50 mmhr-1)에 대한 그래프로서, (d)는 누적강우량 Rac와 Rpeak 과의 상관분석, 회귀분석 그래프, (e)는 누적강우량 Rac와 Train과의 상관분석, 회귀분석 그래프, (f)는 누적강우량 Rac와 Tpeak 과의 상관분석, 회귀분석 그래프이다.
도 8을 참조하면, 두 변수 Rac 와 Rpeak 에서는 피크 강우강도가 10 mmhr-1과 50 mmhr-1의 이상일 때 0.93, 0.77의 높은 상관성을 보였고, Rac 와 Train 에서도 0.80, 0.78의 높은 상관값을 나타냈다. 피크 강우강도가 50 mmhr-1 이상일 때의 Rac 와 Tpeak 에서는 0.49의 다소 낮은 상관성을 보였으나, 피크 강우강도가 10 mmhr-1 이상에서는 0.68의 다소 높은 상관성을 보였다.
한편, 표 4 및 표 5에서는 도 8에서 나타낸 Rac 에 대한 Rpeak, Train, and Tpeak 의 선형 회귀 분석식을 표로 나타내었다.
표 4는 Rpeak ≥ 10 mmhr- 1 에서의 최소, 최적, 최대의 선형 회귀 분석식이며, 표 5는 Rpeak ≥ 50 mmhr- 1 에서의 최소, 최적, 최대의 선형 회귀 분석식이다.
Min. Opt. Max.
Rpeak
Figure 112013090712906-pat00009
Figure 112013090712906-pat00010
Figure 112013090712906-pat00011
Train
Figure 112013090712906-pat00012
Figure 112013090712906-pat00013
Figure 112013090712906-pat00014
Tpeak
Figure 112013090712906-pat00015
Figure 112013090712906-pat00016
Figure 112013090712906-pat00017
Min. Opt. Max.
Rpeak
Figure 112013090712906-pat00018
Figure 112013090712906-pat00019
Figure 112013090712906-pat00020
Train
Figure 112013090712906-pat00021
Figure 112013090712906-pat00022
Figure 112013090712906-pat00023
Tpeak
Figure 112013090712906-pat00024
Figure 112013090712906-pat00025
Figure 112013090712906-pat00026
아울러, 도 9는 도 8에서의 상관계수와 선형회귀분석 그래프로서, (a)는 Rpeak ≥ 10 mmhr- 1 의 경우(일반 대류성 강우의 경우)를, (b)는 Rpeak ≥ 50 mmhr- 1 의 경우(강한 대류성 강우의 경우)를 각각 나타낸다.
이상과 같이, 대류성 강우에서의 단일피크의 특성을 활용하여 향후 현업에서의 국지성 호우에 대한 단시간 예보 및 nowcasting에 대한 능력을 향상시키는데 도움이 될 것으로 사료된다.
이상의 설명은 본 발명의 기술 사상을 예시적으로 설명한 것에 불과한 것으로, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자라면 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 다양한 수정 및 변형이 가능할 것이다. 따라서, 본 발명에 개시된 실시예들은 본 발명의 기술 사상을 한정하기 위한 것이 아니라, 설명하기 위한 것이고, 이러한 실시예에 의하여 본 발명의 기술 사상의 범위가 한정되는 것도 아니다. 본 발명의 보호 범위는 아래 특허청구범위에 의하여 해석되어야 하며, 그와 동등한 범위 내에 있는 모든 기술사상은 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.

Claims (11)

  1. 기상레이더를 통해 레이더 강우강도를 산출하는 단계;
    상기 강우강도에 대한 극 좌표계에서의 기상자료를 직교좌표계에 대한 기상자료로 변환하는 단계;
    Cressman 내삽법을 이용하여 각각의 지리학적 단위 격자마다의 합성 레이더 자료값을 산출하는 단계;
    상기 지리학적 단위 격자마다의 각각의 강수샘플에 대한 누적강우량(Rac) 데이터를 생성하는 단계;
    각각의 강수샘플에 대해서, 불필요한 노이즈 값 제거를 위해, 시간에 따른 강우강도를 FIR(Finite Impulse Response) 필터링하는 단계; 및
    각각의 강수샘플에 대해서, 상기 필터링된 시간에 따른 강우강도를 이용하여 강우피크개수(Npeak)를 포함하는 대류성 강우의 강우특성 정보를 얻는 단계;를 포함하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  2. 제 1항에 있어서,
    상기 기상레이더에서 얻는 강우강도는 1분 단위의 강우강도인 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  3. 제 1항에 있어서,
    상기 산출된 강우강도는 비차등 위상 변위(KDP)에 따른 강우강도인 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  4. 제 1항에 있어서,
    상기 누적강우량(Rac) 데이터를 생성하는 단계에서,
    강우강도가 0.5 mmhr-1 이상인 경우의 강수샘플에 대해서만 누적강우량 데이터를 생성하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  5. 제 1항에 있어서,
    상기 FIR 필터링은 해밍창 함수법을 이용하여 차단 주파수 1/180 Hz 로 강우강도의 단주기 성분을 제거하는 저역통과 필터가 사용되는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  6. 제 1항에 있어서,
    상기 대류성 강우의 강우특성 정보는 강우피크의 수(Npeak), 피크 강우강도(Rpeak), 누적강우량(Rac), 강수 지속시간(Train), 강수시작에서 피크 강우강도까지의 시간(Tpeak)를 포함하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  7. 제 1항에 있어서,
    상기 강우특성정보 중 강우피크의 개수(Npeak)가 1개인 경우만을 추출하는 단계; 및
    추출된 강우피크의 개수가 1인 경우의 강우특성정보를 이용하여, 누적강우량과 다른 강우특성정보와의 상관 관계를 도출하는 단계를 추가로 포함하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  8. 제 1항에 있어서,
    상기 강수샘플은 5분 이상의 강수 지속시간을 지닌 경우에 한정되는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  9. 기상레이더를 통해 1분 단위의 레이더 강우강도를 산출한 후, 산출한 강우강도를 이용하여 지리학적 단위격자마다의 합성 레이더 자료값을 산출하고, 이를 바탕으로 각각의 강수샘플에 대한 누적강우량(Rac) 데이터를 생성하여, 강우피크개수(Npeak)별로 강우특성정보를 분류하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  10. 제 9항에 있어서,
    상기 강우피크개수가 1인 경우만을 추출하여, 상기 강우특성정보 중 누적강우량과 다른 강우특성정보와의 상관 관계를 도출하는 단계를 추가로 포함하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
  11. 제 9항에 있어서,
    상기 강우특성정보는 강우피크의 수(Npeak), 피크 강우강도(Rpeak), 누적강우량(Rac), 강수 지속시간(Train), 강수시작에서 피크 강우강도까지의 시간(Tpeak)를 포함하는 것을 특징으로 하는 기상 레이더 자료를 이용한 대류성 강우특성정보를 얻는 방법.
KR20130119358A 2013-10-07 2013-10-07 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법 KR101469860B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR20130119358A KR101469860B1 (ko) 2013-10-07 2013-10-07 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR20130119358A KR101469860B1 (ko) 2013-10-07 2013-10-07 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법

Publications (1)

Publication Number Publication Date
KR101469860B1 true KR101469860B1 (ko) 2014-12-08

Family

ID=52677918

Family Applications (1)

Application Number Title Priority Date Filing Date
KR20130119358A KR101469860B1 (ko) 2013-10-07 2013-10-07 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법

Country Status (1)

Country Link
KR (1) KR101469860B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101760900B1 (ko) * 2016-01-08 2017-07-24 한국도로공사 기상 레이더 영상을 이용한 어는비 예측 방법
KR101771369B1 (ko) 2016-09-20 2017-08-24 대한민국(기상청장) 기상 데이터 처리 장치 및 이를 이용한 기상 신호 추출을 위한 필터 변수값 산출 방법
CN111398964A (zh) * 2020-04-10 2020-07-10 中国水利水电科学研究院 基于强降水识别与数值大气模式驱动的雷达临近预报方法
CN113558551A (zh) * 2021-09-26 2021-10-29 深圳市资福医疗技术有限公司 一种图像获取方法及应用其的胶囊内窥镜
CN115097459A (zh) * 2022-04-13 2022-09-23 泰州市气象局 S和x波段组网天气雷达反射率因子交互校验方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06214024A (ja) * 1993-01-13 1994-08-05 Toshiba Corp レーダ雨量計
JP2000502794A (ja) * 1995-12-26 2000-03-07 トムソン―セーエスエフ 二偏波レーダを使用する降水率を決定する方法およびそれを利用する気象レーダ
KR100922130B1 (ko) * 2008-11-13 2009-10-16 부경대학교 산학협력단 도플러 기상레이더의 이차 에코 제거 방법

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06214024A (ja) * 1993-01-13 1994-08-05 Toshiba Corp レーダ雨量計
JP2000502794A (ja) * 1995-12-26 2000-03-07 トムソン―セーエスエフ 二偏波レーダを使用する降水率を決定する方法およびそれを利用する気象レーダ
KR100922130B1 (ko) * 2008-11-13 2009-10-16 부경대학교 산학협력단 도플러 기상레이더의 이차 에코 제거 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
김주연 외 2인, 'X 밴드 이중 편파레이더 1분자료를 이용한 대류성 강수시스템의 통계적 분석', 2012년도 한국기상학회 가을학술대회 논문집, 2012.11. *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101760900B1 (ko) * 2016-01-08 2017-07-24 한국도로공사 기상 레이더 영상을 이용한 어는비 예측 방법
KR101771369B1 (ko) 2016-09-20 2017-08-24 대한민국(기상청장) 기상 데이터 처리 장치 및 이를 이용한 기상 신호 추출을 위한 필터 변수값 산출 방법
CN111398964A (zh) * 2020-04-10 2020-07-10 中国水利水电科学研究院 基于强降水识别与数值大气模式驱动的雷达临近预报方法
CN113558551A (zh) * 2021-09-26 2021-10-29 深圳市资福医疗技术有限公司 一种图像获取方法及应用其的胶囊内窥镜
CN115097459A (zh) * 2022-04-13 2022-09-23 泰州市气象局 S和x波段组网天气雷达反射率因子交互校验方法及系统

Similar Documents

Publication Publication Date Title
KR101512015B1 (ko) 기상 레이더 자료를 이용한 대류성 강우를 예측하는 방법
KR101469860B1 (ko) 기상 레이더 자료를 이용한 대류성 강우 특성 정보를 얻는 방법
Isoguchi et al. An L-band ocean geophysical model function derived from PALSAR
Zhang et al. Radar quality index (RQI)—A combined measure of beam blockage and VPR effects in a national network
US9310480B2 (en) Method and arrangement for removing ground clutter
KR101531224B1 (ko) 이중편파 레이더 기반의 강수 추정 시스템 및 그 방법
KR100963532B1 (ko) 기상레이더의 강수량 추정 방법
CN110082842A (zh) 一种降水估测方法及装置
Delobbe et al. Uncertainties in radar echo top heights used for hail detection
Cao et al. Spectrum-time estimation and processing (STEP) for improving weather radar data quality
Vulpiani et al. Characterization of Mediterranean hail-bearing storms using an operational polarimetric X-band radar
Zou et al. Improved algorithms for removing isolated non-meteorological echoes and ground clutters in CINRAD
JP2014102235A (ja) レーダ信号処理装置
Suk et al. Operational real-time adjustment of radar rainfall estimation over the South Korea region
Kim et al. Three-dimensional analysis of the initial stage of convective precipitation using an operational X-band polarimetric radar network
KR102326564B1 (ko) 연직바람 관측장비를 이용한 이중 도플러 레이더의 바람장 보정을 위한 장치 및 방법
Capozzi et al. X-band weather radar monitoring of precipitation fields in Naples urban areas: data quality, comparison and analysis
JP7110003B2 (ja) 処理装置、処理方法、およびプログラム
Kvasov et al. Weather radar data visualization using first-order interpolation
Chen et al. Comparison and Validation of Reflectivity of Dual Linear Polarization SA Radar and CINRAD SA Radar
Chandrasekar et al. Precipitation type determination from spaceborne radar observations
KR101892502B1 (ko) 윈드프로파일러를 이용한 ｓｒｈ 산출 시스템 및 이를 이용한 ｓｒｈ 산출 방법
Lian et al. Denoising algorithm based on local distance weighted statistics for photon counting LiDAR point data
Montopoli et al. Spatial characterization of raincell horizontal profiles from C-band radar measurements at mid-latitude
Chen et al. Evaluation of the rainfall nowcasting system for a dense radar network over Dallas-Fort Worth (DFW)

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20171201

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20181203

Year of fee payment: 5