I. 서 론
저고도 항공기 탐지 및 추적 레이다는 기상 및 지상클러터가 존재하는 환경에서도 표적을 탐지 및 추적하기 위해서 거리 및 속도 분해능이 좋은 PT(pulse train) 파형을 사용한다[1]. PT 파형 사용 시, 시선속도 0 m/s인 클러터의 속도축 null to null 주엽폭은 송신 주파수와 처리 시간(processing time)에 반비례한다. 수신신호가 포화되지 않는 환경에서 클러터 오플롯의 속도는 클러터의 속도 주엽폭보다 작으며, 해당 속도의 오플롯은 제거된다[2]. 본 논문에서는 이 방법을 ZDR(zero Doppler rejection)이라 하며, 클러터의 속도축 주엽폭도 ZDR로 명시한다.
폭우가 내리는 환경에서 수신신호 포화 시, 속도 축에서 높은 파워의 부엽이 다수 발생하고 해당 부엽이 오플롯이 된다[3]. 본 논문에서 제안하는 알고리듬은 폭우가 내리는 환경에서 PT 파형 사용 시, 속도축 첫 번째 부엽까지 발생한 오플롯을 제거한다. 제안한 알고리듬은 기상클러터에 의해 수신신호가 포화되는 것을 탐지하고, 기상클러터의 시선속도를 추정하여, 수신신호 포화 시 속도축 첫 번째 부엽이 발생한 속도를 계산한 후, 해당 속도 범위 안의 오플롯을 제거한다. 해당 속도범위와 오플롯을 제거하는 알고리듬은 VDR(variable Doppler rejection)이라 한다. 제안된 알고리듬은 그림 1과 같이 CFAR(constant false alarm rate) 탐지 이후에 적용된다. 그림 2는 제안된 알고리듬의 순서도이며, 4단계로 구성된다. ZDR과 VDR로 기상클러터 오플롯이 제거되는 속도 범위를 비교하면 그림 3과 같다.
- Step 1. 기상클러터 clutter to noise ratio(CNR) 계산
- Step 2. 기상클러터 시선속도 추정
- Step 3. 기상클러터 포화 판단
- Step 4.1. 기상클러터 미포화 시, ZDR 적용
- Step 4.2. 기상클러터 포화 시, VDR 적용
II. 기상클러터 대응 알고리듬
PT 파형에 대한 펄스압축 및 도플러 처리 결과인 a(n, k)를 활용하여 속도 게이트별(k) 평균 CNR ψ(k)는 식 (1)과 같이 계산한다.
여기서, n은 거리 게이트, Nr는 거리 게이트의 개수, 그리고 σ2은 사전에 측정된 배경 잡음의 파워이다.
ψ(k)의 최댓값이 발생한 off-grid 속도 게이트 kc를 추정하여 기상클러터의 시선속도를 그림 4와 같이 계산하며, 이를 정리하면 식 (2)~식 (4)와 같다.
여기서, k0은 ψ(k)의 최댓값이 발생한 on-grid 속도 게이트, λ는 파장, vc는 기상클러터의 시선속도, fPRF는 PRF(pulse repetition frequency), 그리고 NFFT는 PT 파형의 도플러 처리 시 적용된 FFT(fast Fourier transform) 사이즈이다.
ψ(k)와 ADC(analog to digital converter)의 동적 범위 γ를 활용한 기상클러터의 포화 판단 조건은 식 (5)와 같다.
여기서, ρ는 기상클러터 신호 크기의 변화폭을 의미하며, 본 논문에서는 5 dB로 설정하고, 동적 범위 γ 보다 5 dB 작은 값보다 큰 값을 수신하면, 기상클러터 수신신호가 포화된 것으로 판단한다.
기상클러터에서 발생한 오플롯을 제거하는 속도 범위는 다음과 같다. 식 (6)은 기상클러터 미포화 시 사용되며, 식 (7)은 기상클러터 포화 시 사용된다. 식 (6) 및 식 (7)의 범위를 비교한 예시는 그림 3에서 확인할 수 있다.
여기서, kp는 탐지된 플롯의 off-grid 속도 게이트, ΔZDR은 PT 사용 시 클러터의 속도축 peak to null 속도 게이트 범위이다. 기상클러터 포화 시 속도축 첫 번 째 부엽의 속도 게이트는 기상클러터 시선속도의 약 3배이며, 식 (3)의 추정값으로 식 (8)의 ΔVDR 값을 계산하였다.
III. 실측 시험 결과
본 장에서는 레이다에 기상클러터 오플롯 제거 알고리듬을 적용한 결과를 소개한다. 맑은 날씨와 폭우가 쏟아지는 날씨에서 시험한 결과를 소개한다.
Variable | Value and description |
---|---|
Surveillance area | 135°~225° |
Frequency band | S |
Waveform | Pulse train |
Test date | 2023.06.01. - Sky clear, 2023.06.29. - Heavy rain. |
그림 5 ~ 그림 8은 맑은 날씨와 폭우의 환경에서 수신신호를 저장하여 분석한 시험 결과이다. 그림 7은 맑은 날씨에서의 Max of Mean CNR 값과 VDR/ZDR의 범위를 보여주고 있다. 맑은 날씨에서는 CNR 값이 대체적으로 낮게 분포되어 기상클러터 포화 조건을 만족하지 못하였고, 그에 따라 ZDR을 적용하여 오플롯 제거를 수행하였다. 수행한 결과는 그림 7과 같다.
그림 6은 강우 환경에서의 결과로 그림 5와 비교해 보면 Max of mean CNR 값이 상대적으로 높게 분포되어 있는 것을 확인할 수 있다. 그래서 앞서 언급한 기상클러터의 포화 판단조건 식 (5)를 만족하였고, 기상클러터에 의해 수신신호가 포화된 것으로 판단하여 식 (7)을 적용하여 탐지된 오플롯을 제거하였다. 그림 8은 VDR을 적용해서 오플롯을 제거한 결과이다. ZDR을 적용했을 때보다 도플러 제거 영역이 늘어났고, 늘어난 영역에 포함된 오플롯이 제거되었다.
IV. 결 론
본 논문에서는 강우 환경에서 기상클러터에 의한 레이다의 성능 저하를 개선하기 위해 기상클러터 포화로 발생하는 오플롯 제거 방법을 제안하였다. 먼저, 기상클러터 환경에서의 CNR 값 계산과 기상클러터의 포화 판단조건을 설명하고, 시선속도 추정을 통해 도플러 제거 영역을 넓히는 알고리듬을 제안하였다. 그리고 레이다에 제안하는 알고리듬을 적용하여 빔방사 후 수신신호를 저장하여 결과를 분석하였다. 분석결과 포화 기상클러터에 의해 발생한 오플롯을 제거하기 위한 방안으로 오플롯 제거영역을 고정하는 ZDR은 포화 기상클러터로 인해 부엽으로 확산된 오플롯들을 제거하지 못하였다. 하지만 본 논문에서 제안한 VDR은 도플러 제거 영역을 가변하여 약 90개 정도의 오플롯이 추가로 제거됨을 확인하였다. 오탐지 플롯들을 제거함으로써 실표적의 탐지 기회는 늘어나 포화 기상클러터의 대응으로써 효과적으로 판단된다.
알고리듬 성능개선을 위해 기상클러터 환경에서 더 많은 시험을 통한 데이터 축적과 분석을 하여 적절한 임계치를 설정하고, 축적한 데이터를 기반하여 VDR 적용 범위를 조절하면, 더 향상된 성능의 기상클러터 탐지와 대응 방안이 될 것으로 기대된다.