AIDT 온디바이스 STEM 화산 분화 초저주파 음파(Infrasound) 어레이 F-K 파수 분석 서보 엔진

1. 개요: 화산 폭발의 음향 지문 초저주파 음파(Infrasound)

화산 분화 시 수증기와 화산재가 음속 이상의 속도로 폭발 분출될 때 인간의 가청 주파수 이하 대역인 초저주파 음파(Infrasound, \(0.1 \sim 10\,\text{Hz}\))가 강력하게 발생합니다. 초저주파 음파는 대기 점성에 의한 흡수 감쇠가 극히 적어 수백~수천 킬로미터 떨어진 거리까지 전파되므로, 구름이나 화산재 연기에 가려진 마그마 분출구를 원격 감시하는 최고의 지구물리 도구입니다. 수백 미터 간격으로 설치된 차압 마이크로폰 어레이를 통해 도달 시간차(TDOA)와 주파수-파수(Frequency-Wavenumber, F-K) 스펙트럼을 실시간 합성함으로써 음원의 방위각(Back-azimuth)과 위상 속도를 정밀 추적할 수 있습니다.

2. 주파수-파수(F-K) 슬로우니스 빔포밍 수식 모델

평면파가 센서 어레이로 입사할 때 음파의 전파 방향과 속도를 나타내는 슬로우니스(Slowness) 벡터 \(\vec{s} = (s_x, s_y)\)는 다음과 같습니다:

\vec{s} = \frac{\hat{k}}{c} = \left( \frac{\sin \theta}{c}, \frac{\cos \theta}{c} \right), \quad |\vec{s}| = \frac{1}{c_{app}}

센서 위치 \(\vec{r}_j = (x_j, y_j)\)에서의 신호 \(x_j(t)\)의 푸리에 변환 \(X_j(f)\)에 대해 주파수-파수 공간 빔 출력 파워 \(P(f, \vec{s})\)는 다음과 같이 계산됩니다:

P(f, \vec{s}) = \left| \sum_{j=1}^M X_j(f) \exp\left( i 2\pi f \vec{s} \cdot \vec{r}_j \right) \right|^2

음원 방위각 \(\theta\)와 외견 위상 속도 \(c_{app}\)는 최대 빔 파워를 형성하는 슬로우니스 \(\vec{s}_{opt}\)로부터 도출되며, 분출구에서의 초과 압력 \(\Delta p(t)\)로부터 분출 체적 유량 \(Q(t)\)가 적분 역산됩니다:

\theta = \arctan\left( \frac{s_x}{s_y} \right), \quad Q(t) = \frac{2\pi r}{\rho_0} \int_0^t \Delta p\left( t' + \frac{r}{c} \right) dt'

3. 0-Alloc C# 화산 초저주파 어레이 서보 엔진 구현

4원 마이크로폰 어레이의 교차 상관 행렬을 10,000Hz 주기로 연산하여 방위각과 유량을 추출하는 ref struct 아키텍처입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.VolcanicAcoustics
{
    public readonly ref struct InfrasoundArrayState
    {
        public readonly float BackAzimuthDeg;
        public readonly float ApparentVelocityMs;
        public readonly float AcousticPowerPa2;
        public readonly float EruptionFlowRateM3S;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public InfrasoundArrayState(float azDeg, float vApp, float pPa2, float qM3s)
        {
            BackAzimuthDeg = azDeg;
            ApparentVelocityMs = vApp;
            AcousticPowerPa2 = pPa2;
            EruptionFlowRateM3S = qM3s;
        }
    }

    public ref struct VolcanicInfrasoundServo
    {
        private const float AirDensityKgM3 = 1.225f;
        private const float SoundSpeedNominal = 340.0f; // m/s
        private const float DistanceToVentM = 15000.0f;  // 15 km

        private float _estimatedAzimuthDeg;
        private float _smoothedExcessPressurePa;

        public VolcanicInfrasoundServo(float initialAzimuthDeg = 45.0f)
        {
            _estimatedAzimuthDeg = initialAzimuthDeg;
            _smoothedExcessPressurePa = 0.0f;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public InfrasoundArrayState UpdateServo(float tdoaXSec, float tdoaYSec, float peakPressurePa, float dt = 0.0001f)
        {
            // Array aperture ~ 200m: sx = tdoaX / 200, sy = tdoaY / 200
            float sx = tdoaXSec / 200.0f;
            float sy = tdoaYSec / 200.0f;
            float sMag = MathF.Sqrt(sx * sx + sy * sy);
            sMag = MathF.Max(1e-6f, sMag);

            float currentAzimuthDeg = MathF.Atan2(sx, sy) * (180.0f / MathF.PI);
            if (currentAzimuthDeg < 0.0f) currentAzimuthDeg += 360.0f;

            _estimatedAzimuthDeg += (currentAzimuthDeg - _estimatedAzimuthDeg) * (dt * 30.0f);
            float apparentSpeed = 1.0f / sMag;

            // Acoustic power & mass flow rate Q = (2 * pi * r / rho0) * int(p dt)
            _smoothedExcessPressurePa += (peakPressurePa - _smoothedExcessPressurePa) * (dt * 10.0f);
            float flowRateM3S = (2.0f * MathF.PI * DistanceToVentM / AirDensityKgM3) * (_smoothedExcessPressurePa / SoundSpeedNominal);

            float acousticPower = peakPressurePa * peakPressurePa;

            return new InfrasoundArrayState(_estimatedAzimuthDeg, apparentSpeed, acousticPower, flowRateM3S);
        }
    }
}

4. 실시간 화산 분화 초저주파 어레이 F-K 극좌표계 시뮬레이터

음원 방위각과 폭발 강도를 조절하며 2차원 슬로우니스 평면 상의 F-K 빔 피크 및 분출 가스 체적 유량을 관측하십시오.

화산 분화구 방위각 (Back-Azimuth): 65° (ENE) 분화구 과압 (Excess Pressure Δp): 18 Pa

5. 교육 연계 및 방재/국방 음향 감시 실무적 의의

2022 개정 교육과정 지구과학 I/II: 화산 활동과 판 구조론, 화산 쇄설물 및 분화 기둥 역학, 대기 파동과 원격 지구물리 탐사.
물리학 II: 파동의 도플러 효과, 위상 속도와 슬로우니스, 마이크로폰 어레이 신호 처리와 공간 빔포밍(Beamforming).
국방 음향 정찰 / 포병 탐지: 적 포병 사격 원점 역추적(음향 표적 탐지 대포병 레이더), 지하 핵실험 초저주파 대기 감시망.