AIDT 온디바이스 STEM 사냐크 간섭계(Sagnac) 광섬유 자이로 초정밀 회전 지진계 서보 엔진

1. 개요

지진학에서 기존 지진계는 지면의 병진 운동(Translational Motion)만을 계측하여 지진파의 회전 성분(Rotational Motion)으로 인한 구조물 비틀림 붕괴 위험을 정확히 예측하기 어렵습니다. 사냐크 효과(Sagnac Effect)에 기반한 폐쇄 루프 간섭형 광섬유 자이로스코프(Interferometric Fiber-Optic Gyroscope, IFOG)는 1 km 이상의 편광유지(PM) 광섬유 코일 내부를 시계 방향(CW)과 반시계 방향(CCW)으로 동시 전파하는 두 광파 간의 상대적 위상차(\(\Delta \Phi_s\))를 측정함으로써 \(10^{-9}\text{ rad/s}\) 미만의 극초정밀 회전 지진파를 지면 가속도 간섭 없이 순수하게 계측합니다. 그러나 슈페(Shupe) 열 비가역 효과 및 광섬유 내 비선형 커(Kerr) 효과는 장기 드리프트를 유발합니다. 본 고에서는 디지털 위상 램프(Digital Phase Ramp) 피드백 서보를 10,000Hz 0-Alloc C# 엔진으로 구현합니다.

2. 수학적 모델링 및 지배 방정식

유효 반경 \(R\), 광섬유 길이 \(L\), 광파장 \(\lambda_0\)인 광섬유 코일이 각속도 \(\Omega(t)\)로 회전할 때 발생하는 사냐크 위상 천이 \(\Delta \Phi_s\)는 다음과 같습니다:

\Delta \Phi_s = \frac{4\pi L R}{\lambda_0 c} \Omega(t) = \frac{8\pi A N}{\lambda_0 c} \Omega(t)

광검출기에 도달하는 간섭 신호 세기는 변조 위상 \(\phi_m(t)\)와 위상 램프 피드백 \(\Delta \phi_{ramp}(t)\)에 의해 결정됩니다:

I(t) = \frac{I_0}{2} \left[ 1 + \cos\left( \Delta \Phi_s + \phi_m(t) - \phi_m(t - \tau) - \Delta \phi_{ramp} \right) \right]

여기서 \(\tau = \frac{n L}{c}\)는 광섬유 루프 내 광 지연 시간입니다. 10,000Hz 폐루프 서보 루프에서 \(\Delta \Phi_s - \Delta \phi_{ramp} = 0\)이 되도록 집적광학 위상변조기(MIOC)의 계단파 램프 전압 기울기(\(\frac{d\phi}{dt} = \frac{\Delta \Phi_s}{\tau}\))를 정밀 보정합니다.

3. 실시간 피드백 제어 아키텍처

+-----------------------------------------------------------------------------+
|    AIDT Sagnac Fiber-Optic Gyro Seismometer 10,000Hz Real-Time Servo Engine |
+-----------------------------------------------------------------------------+
|  [InGaAs Photodetector & Anti-Aliasing ADC] -> [Lock-In Demodulator]        |
|                                       |                                     |
|                                       v                                     |
|         [Closed-Loop Sagnac Phase Error & Shupe Thermal Estimator]          |
|                                       |                                     |
|                                       v                                     |
|      [Zero-Alloc Discrete Digital Phase Ramp & 2\pi Reset Modulator]        |
|                                       |                                     |
|                                       v                                     |
|   [16-Bit High-Speed DAC] -> [LiNbO3 Multi-Function Integrated Optical Chip]|
+-----------------------------------------------------------------------------+

4. Zero-Allocation C# 서보 엔진 구현

사냐크 간섭계의 널링(Nulling) 동작점 유지를 위한 C# 고속 폐루프 서보 클래스입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Seismology.Sagnac
{
    public readonly struct SagnacConfig
    {
        public readonly double ScaleFactorRadPerDps;
        public readonly double LoopTransitTimeSec;
        public readonly double Kp;
        public readonly double Ki;

        public SagnacConfig(double scaleFactorRadPerDps, double loopTransitTimeSec, double kp, double ki)
        {
            ScaleFactorRadPerDps = scaleFactorRadPerDps;
            LoopTransitTimeSec = loopTransitTimeSec;
            Kp = kp;
            Ki = ki;
        }
    }

    public ref struct SagnacSeismometerServo
    {
        private readonly SagnacConfig _config;
        private double _phaseIntegral;
        private double _rampPhaseRad;

        public SagnacSeismometerServo(in SagnacConfig config)
        {
            _config = config;
            _phaseIntegral = 0.0;
            _rampPhaseRad = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public (double estimatedRotationRateDps, double feedbackRampRad) Update(double demodulatedPhaseErrorRad, double dtSec)
        {
            _phaseIntegral += demodulatedPhaseErrorRad * dtSec;
            _phaseIntegral = Math.Clamp(_phaseIntegral, -100.0, 100.0);

            // 널링을 위한 피드백 램프 기울기 산출
            double rampSlope = (_config.Kp * demodulatedPhaseErrorRad) + (_config.Ki * _phaseIntegral);
            _rampPhaseRad += rampSlope * dtSec;

            // 2pi 위상 리셋 연산
            while (_rampPhaseRad >= 2.0 * Math.PI) _rampPhaseRad -= 2.0 * Math.PI;
            while (_rampPhaseRad < 0.0) _rampPhaseRad += 2.0 * Math.PI;

            double estimatedRate = (rampSlope * _config.LoopTransitTimeSec) / _config.ScaleFactorRadPerDps;

            return (estimatedRate, _rampPhaseRad);
        }
    }
}

5. 학술적 의의 및 교육적 활용 가치

본 설계는 2022 개정 물리 및 지구과학 교과 내 '특수 상대성 이론과 사냐크 효과', '지진파의 P/S파 및 회전 운동', '간섭계의 위상 제어 원리'를 다루는 융합 교육 플랫폼으로 활용됩니다. 광섬유 속을 회전 전파하는 빛의 결맞음 간섭을 활용하여 지진동의 6자유도 완벽 계측을 실감나게 학습할 수 있습니다.