2022 개정 교육과정 연계 온디바이스 성층권대기물리·광학계측 STEM 시리즈 #58

브루어 분광광도계 성층권 오존 전량(DU) 서보 엔진

5파장 자외선 차분 흡수비(Beer-Lambert) 및 0-Alloc C# 태양 직사광 천정각 추적 슬릿 마스크 제어 루프

ZeliDesk 대기광학·성층권환경팀 오존홀원격탐사·자외선분광 분과

2026년 9월 12일

1. 다파장 베어-람베르트 법칙 및 가중 차분 오존 산출 수리 모델

각 파장 $\lambda_i$에서 지표에 도달하는 태양 직사광 세기 $I_i$는 대기 상단 강도 $I_{0, i}$, 오존 흡수 계수 $\alpha_i$, 레일리 산란 계수 $\beta_i$, 에어로졸 산란 계수 $\delta_i$에 의해 베어-람베르트 법칙을 따릅니다:

여기서 $\Omega$는 총 오존량($\text{cm STP}$), $\mu$는 성층권 오존 광로(Air Mass Factor), $m_R, m_A$는 대기 질량 인자입니다. 레일리 산란과 에어로졸의 파장 의존성이 선형적으로 상쇄되도록 4개 파장에 가중치 $w_i = [1.0, -0.5, -2.2, 1.7]$을 부여하여 합산 함수 $F$를 구성합니다:

이에 따라 최종 오존 전량 $\Omega$ ($\text{DU}$)는 기상 간섭 요인이 완전 배제된 다음 공식으로 추출됩니다:

2. 0-Alloc C# 스택 기반 브루어 오존 모니터링 서보

광전증폭관(PMT) 파장별 펄스 계수로부터 차분 가중합을 계산하고 방위각 추적 오차를 1,000Hz로 보정하는 ref struct 루프입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.StratosphericPhysics
{
    public readonly ref struct BrewerSpectrometerTelemetry
    {
        public readonly double Count3101Nm;
        public readonly double Count3135Nm;
        public readonly double Count3168Nm;
        public readonly double Count3201Nm;
        public readonly double SolarZenithAngleDeg;
        public readonly double StationPressureMbar;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public BrewerSpectrometerTelemetry(double c1, double c2, double c3, double c4, double szaDeg, double pMbar)
        {
            Count3101Nm = c1;
            Count3135Nm = c2;
            Count3168Nm = c3;
            Count3201Nm = c4;
            SolarZenithAngleDeg = szaDeg;
            StationPressureMbar = pMbar;
        }
    }

    public readonly ref struct OzoneMeasurementResult
    {
        public readonly double TotalOzoneDobsonUnits;
        public readonly double AirMassFactorMu;
        public readonly double WeightedRatioF;
        public readonly bool IsOzoneHoleCondition;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OzoneMeasurementResult(double o3Du, double mu, double ratioF, bool holeAlert)
        {
            TotalOzoneDobsonUnits = o3Du;
            AirMassFactorMu = mu;
            WeightedRatioF = ratioF;
            IsOzoneHoleCondition = holeAlert;
        }
    }

    public ref struct BrewerOzoneServo
    {
        private const double F0Calibration = 3120.5; // 외기권 기기 교정 상수
        private const double DeltaAlpha = 0.342;     // 유효 차분 흡수 계수
        private const double EarthRadiusKm = 6371.0;
        private const double OzoneLayerHeightKm = 22.0;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OzoneMeasurementResult Step(in BrewerSpectrometerTelemetry telemetry)
        {
            // 1. 성층권 오존 대기 질량 인자 (mu) 기하학적 굴절 모델: mu = (1 + h/R) / sqrt(cos^2(SZA) + 2h/R)
            double szaRad = telemetry.SolarZenithAngleDeg * (Math.PI / 180.0);
            double cosSza = Math.Cos(szaRad);
            double ratioH = OzoneLayerHeightKm / EarthRadiusKm;
            double mu = (1.0 + ratioH) / Math.Sqrt(cosSza * cosSza + 2.0 * ratioH);

            // 2. 파장별 로그 광량 계산 (log10 Counts * 10000)
            double logI1 = Math.Log10(Math.Max(1.0, telemetry.Count3101Nm));
            double logI2 = Math.Log10(Math.Max(1.0, telemetry.Count3135Nm));
            double logI3 = Math.Log10(Math.Max(1.0, telemetry.Count3168Nm));
            double logI4 = Math.Log10(Math.Max(1.0, telemetry.Count3201Nm));

            // 3. 4파장 가중 차분합 F: w = [1.0, -0.5, -2.2, 1.7]
            double fRatio = 1.0 * logI1 - 0.5 * logI2 - 2.2 * logI3 + 1.7 * logI4;
            fRatio *= 10000.0;

            // 4. 오존 전량 (DU) 산출: Omega = (F0 - F) / (DeltaAlpha * mu)
            double ozoneDu = (F0Calibration - fRatio) / Math.Max(0.1, DeltaAlpha * mu);

            // 5. 오존홀 경보 판별 (220 DU 이하 시 오존홀 조건)
            bool isOzoneHole = ozoneDu < 220.0;

            return new OzoneMeasurementResult(ozoneDu, mu, fRatio, isOzoneHole);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용