초임계 CO2 심부 지열 발전(EGS) 자연 순환 써모사이폰 서보 엔진

1. 지구물리학적 배경: 초임계 CO2 EGS와 써모사이폰 원리

지하 3~5km 깊이의 고온 건조 암반(Hot Dry Rock, $T > 200^\circ\text{C}$)에 $\text{sCO}_2$를 주입하면 다음과 같은 밀도 반전 사이클이 성립합니다:

1. 주입관(Cold Downcomer): 지상 냉각기 통과 후 고밀도($\rho_{\text{cold}} \approx 850\text{--}950\text{ kg/m}^3$) 상태의 액체형 $\text{sCO}_2$가 중력에 의해 자발적으로 강하합니다.
2. 파쇄대 열교환(Reservoir): 고온 암반 균열망을 통과하며 급격히 가열되어 기체형 초임계 상태로 팽창합니다.
3. 생산관(Hot Riser): 열팽창에 의해 밀도가 급감한 저밀도($\rho_{\text{hot}} \approx 300\text{--}450\text{ kg/m}^3$) $\text{sCO}_2$가 지상 터빈을 향해 맹렬히 부력 상승합니다.

2. 수학적 모델링: Darcy 마찰 손실과 자발 순환 유량 역산

수 킬로미터 관로의 수두차에 의한 자연 순환 구동 압력은 수십 bar($> 3\text{--}5\text{ MPa}$)에 달하여 순환 펌프 없이도 파쇄대의 Darcy 마찰 압력 강하를 완전히 상쇄합니다:

지상 팽창 배압 밸브 개도율을 10,000Hz 주기로 정밀 변조하여 핀치 포인트(Pinch Point) 온도 교차를 차단하고 최대 발전 효율 작동점을 유지합니다:

3. 10,000Hz 0-Alloc C# 써모사이폰 웰헤드 배압 서보

유정 두부(Wellhead) 압력과 온도를 실시간 감시하여 0.1ms 내에 배압 조절 밸브를 제어하는 C# 무할당 엔진입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.GeothermalEnergy
{
    public readonly ref struct GeothermalWellheadFrame
    {
        public readonly double InjectionPressureMpa;
        public readonly double ProductionPressureMpa;
        public readonly double ProductionTempCelsius;
        public readonly double MassFlowRateKgPerSec;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public GeothermalWellheadFrame(double pInj, double pProd, double tProd, double mDot)
        {
            InjectionPressureMpa = pInj;
            ProductionPressureMpa = pProd;
            ProductionTempCelsius = tProd;
            MassFlowRateKgPerSec = mDot;
        }
    }

    public ref struct Sco2ThermosiphonServo
    {
        private const double TargetBuoyancyHeadMpa = 4.5;
        private double _backPressureValvePercent;

        public Sco2ThermosiphonServo(double initValvePercent = 45.0)
        {
            _backPressureValvePercent = initValvePercent;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void UpdateWellheadControl(in GeothermalWellheadFrame frame, Span driveOut)
        {
            // 유효 열사이폰 구동 차압 산출
            double netBuoyancyMpa = frame.ProductionPressureMpa - frame.InjectionPressureMpa;
            double buoyancyDeficit = TargetBuoyancyHeadMpa - netBuoyancyMpa;

            // 10,000Hz 고속 배압 조절: 유량이 급감하면 밸브를 열어 순환 유량을 유도
            if (frame.MassFlowRateKgPerSec < 50.0)
            {
                _backPressureValvePercent = Math.Clamp(_backPressureValvePercent + 0.05, 10.0, 90.0);
            }
            else
            {
                _backPressureValvePercent = Math.Clamp(_backPressureValvePercent - buoyancyDeficit * 0.02, 10.0, 90.0);
            }

            driveOut[0] = _backPressureValvePercent; // 지상 배압 밸브 개도율 (%)
            driveOut[1] = netBuoyancyMpa;             // 순 써모사이폰 차압 (MPa)
            driveOut[2] = frame.MassFlowRateKgPerSec; // 질량 유량 (kg/s)
        }
    }
}

4. 교육 및 탄소중립 신에너지 응용 가치