난류 극한 현상의 예측 가능성 계층 구조와 기계학습

**1. 서론** 대기·해양 등 자연계의 고차원·다중스케일 시스템에서 극한 현상의 조기 경보는 사회·경제적 파급 효과가 크다. 기존 예측 가능성 연구는 Lyapunov 지수와 같은 전역적인 불안정성 지표에 의존하거나, 초기 조건을 미세하게 교란한 앙상블을 통해 지역적 예측 가능성을 추정한다. 그러나 이러한 방법은 (i) 정확한 지배 방정식이 필요하고, (ii) 계산 비용이 매우 높으며, (iii) 초기 교란 방식에 민감하다는 한계가 있다. 최근 데이터‑구동 접근법—극값 이론 기반 로컬 지표, 시간 지연 재현성, 정보 이론적 발산—이 제안되었지만, 이들은 시스템 전체의 평균적 특성을 다룰 뿐 개별 사건에 대한 예측 가능성 한계를 제공하지 못한다. **2. 방법론** 본 연구는 2차원 콜모고로프 흐름(Kolmogorov flow)을 실험 플랫폼으로 선택한다. 이 흐름은 128 × 128 격자, Re = 100, 강제 파수 n_f = 4 로 설정된 비압축성 Navier‑Stokes 방정식을 직접수치시뮬레이션(DNS)한다. enstrophy Ω = ⟨ω²⟩가 상위 1 %를 차지하는 순간을 극한 사건으로 정의하고, 총 211개의 사건을 추출한다. **2.1 자동회귀 조건부 확산 모델** 시간 t의 속도장 xₜ 을 입력으로, 확산 과정(노이즈 추가 → 역확산을 통한 복원)을 여러 단계에 걸쳐 수행한다. 각 단계는 U‑Net 구조로 구현되며, 손실 함수는 변분 베이즈 프레임워크에 기반한 ELBO를 최소화한다. 학습 후, 모델은 현재 상태 xₜ 으로부터 미래 상태 xₜ₊₁ 의 확률 분포 p(xₜ₊₁|xₜ) 를 샘플링한다. 이 과정을 재귀적으로 적용해 원하는 리드 타임 T 까지의 앙상블을 생성한다. **2.2 예측 기술 평가와 예측 가능성 한계 정의** 각 사건 e 에 대해 리드 타임 T (= tₑ − t)에서 모델이 생성한 enstrophy 분포 D 에 대해 CRPS_Te 를 계산한다. 기준점으로는 장기 enstrophy 분포 F_L 에 기반한 CRPS_ref를 사용한다. 사건별 기술 점수는 Sₑ(T) = 1 − CRPS_Te / CRPS_ref 으로 정의되며, Sₑ(T) > 0인 최대 T 을 예측 가능성 한계 T*ₑ 라 한다. **3. 결과** *예측 기술 히트맵* – 사건별 Sₑ(T)와 리드 타임을 시각화한 히트맵에서, 대부분의 사건은 t − tₑ ≈ −7 (시뮬레이션 시간 단위)에서 기술이 급격히 감소한다. 이는 약 1–2 Lyapunov 시간에 해당한다. 그러나 일부 사건은 t − tₑ ≈ −15 까지 높은 Sₑ(T)를 유지, 즉 > 4 Lyapunov 시간에 걸친 예측 가능성을 보인다. *스펙트럴 필터링 실험* – 초기 조건에 대해 isotropic spectral cut‑off k_c 을 적용해 작은 스케일을 차례로 제거하였다. λ_c = 0.30 L (전체 파장 대비 30 %) 이하의 대규모 모드만 남겨도 원본과 거의 동일한 T*ₑ 분포를 유지한다(R ≥ 0.84). 반대로 λ_c > 0.30 L 인 경우, 모든 사건군에서 T*ₑ가 크게 감소하고, Wasserstein 거리와 permutation test을 통해 통계적으로 유의미한 차이가 확인된다. 이는 예측 가능성이 주로 대규모 구조에 의해 지배된다는 증거이다. *코히어런트 구조 분석* – 고예측 사건을 역추적하면, 사건 직전 강한 변형(strain) 코어가 형성되고, 이 코어 주변에 사중극(quadrupolar) 소용돌이 패킷이 조직된다. 이 패킷은 평균 ≈ 10 Lyapunov 시간(시뮬레이션 시간 기준) 동안 지속되며, 변형 코어와 소용돌이의 수명이 길수록 T*ₑ가 크게 증가한다. 반면, 저예측 사건은 변형 코어가 짧게 존재하거나, 소용돌이 패킷이 빠르게 붕괴되는 특징을 보인다. **4. 논의** 본 연구는 (i) 데이터‑구동 확산 모델이 복잡한 난류 흐름에서도 사건별 예측 가능성 한계를 정량화할 수 있음을, (ii) 대규모 코히어런트 구조—특히 변형 코어와 사중극 소용돌이—가 극한 사건의 예측 가능성을 결정한다는 물리적 메커니즘을 제시한다. 기존의 Lyapunov 기반 전역 불안정성 지표와 달리, 여기서 제시한 사건‑별 기술 점수는 관측만으로도 적용 가능하므로, 실제 대기·해양 관측망에 바로 활용될 수 있다. 또한, 스펙트럴 필터링 결과는 저해상도(대규모 모드만 보존) 모델에서도 핵심 예측 정보를 유지할 수 있음을 시사한다. **5. 결론** 자동회귀 조건부 확산 모델과 CRPS 기반 사건‑별 기술 점수를 결합한 프레임워크는 난류 극한 현상의 예측 가능성 계층을 명확히 드러낸다. 대규모 변형 코어와 사중극 소용돌이 패킷의 지속성이 예측 가능성의 물리적 근거임을 확인함으로써, 향후 관측 기반 조기 경보 시스템 설계와 저비용 모델링 전략에 중요한 인사이트를 제공한다.