불확도 예산은 ‘오차를 크게 잡는 표’가 아니다

탄소감축 결과에 불확도를 붙이라고 하면 흔히 센서 사양서의 정확도, 예를 들어 ±2%를 가져와 최종 감축량에도 그대로 붙인다. 그러나 감축량은 센서 한 번의 판독값이 아니다. 기준기간 배출량에서 사업기간 배출량을 빼고, 경우에 따라 누출·프로젝트 배출·사육두수·생산량 변화까지 반영한 계산 결과다. 농도계가 매우 정확해도 환기량을 잘못 추정하거나 기준선이 비교 기간을 대표하지 못하면 감축량은 크게 빗나갈 수 있다.

국제측정용어집 VIM은 불확도 예산을 측정 불확도와 그 구성요소, 계산·결합 방법을 기술한 것으로 정의한다. 여기에는 측정모델, 입력량별 추정값과 불확도, 공분산, 확률분포, 자유도, 평가 유형, 필요할 경우 포함인자까지 들어간다. 따라서 불확도 예산은 막연한 ‘안전 마진’도, 결과를 약하게 만드는 면책문구도 아니다. 최종 숫자가 어떤 정보에 의존하며 어느 입력을 개선해야 결과가 실제로 좋아지는지를 보여주는 기술 문서다.

또 하나의 오해는 불확도와 오류를 같은 말로 쓰는 것이다. 잘못 입력한 사육두수, ppm과 %의 단위 혼동, 시계 불일치, 교정계수의 중복 적용은 발견해 수정해야 할 오류다. 반면 보정 후에도 남는 센서 교정의 한계, 개체 간 자연 변이, 표본이 모집단을 대표하는 정도처럼 현재 정보로 완전히 제거할 수 없는 분산은 불확도로 표현한다. 알려진 편향을 고치지 않고 단순히 불확도 구간만 넓히는 방식은 좋은 설계가 아니다.

첫 단계는 계산식보다 측정대상을 한 문장으로 고정하는 일이다

불확도 예산은 무엇을 측정하는지, 즉 측정량(measurand)을 먼저 정해야 한다. ‘메탄 감축량’만으로는 부족하다. 예를 들어 “A동 비육우 120두가 2026년 6월 1일부터 30일까지 배출한 메탄 질량을, 같은 생산단위로 정규화한 기준선과 비교한 월간 감축량(kg CH₄)”처럼 장소, 축군, 기간, 경계, 단위와 비교 규칙을 명시해야 한다. 장내발효만 포함하는지, 분뇨 저장에서 나온 메탄도 포함하는지, 외부 배경농도를 어떻게 빼는지도 이 문장에 연결돼야 한다.

그다음 최종 결과를 입력량의 함수로 쓴다. 단순화한 축사 배출률은 다음과 같이 표현할 수 있다.

E = Q × (x_out − x_in) × (M_CH4 / V_m) × t

여기서 Q는 건식 기준 체적유량, x_out과 x_in은 배출구와 유입구의 메탄 몰분율, M_CH4는 메탄 분자량, V_m은 해당 온도·압력의 몰부피, t는 적산 시간이다. ppm 농도는 10⁻⁶을 곱해 무차원 몰분율로 바꾸며 농도와 유량의 건식·습식 기준을 일치시켜야 한다. 감축량 R은 기준선 배출량 E_base와 사업 배출량 E_project의 차이에 경계 안의 추가 배출이나 누출 조정치를 반영해 계산한다. 실제 자연환기 축사는 하나의 배출구가 없고 기류가 역전되기도 하므로, 이 식을 그대로 쓰기보다 추적가스법, CO₂ 수지, 질량수지 또는 분산모델 등 채택한 방법의 완전한 모델을 기록해야 한다.

이 모델을 먼저 쓰면 빠진 항목이 드러난다. 농도 센서만 예산에 넣고 환기량을 빼놓았는지, 기준기간과 사업기간의 사육두수 보정이 별도 모델인지, 결측 구간을 메운 예측값이 관측값처럼 처리됐는지 확인할 수 있다. 측정모델의 버전과 적용 기간도 남겨야 한다. 계산식이 바뀌면 같은 원자료에서도 결과와 불확도가 달라지기 때문이다.

불확도 원인은 측정 사슬 전체에서 찾는다

실무에서는 원인을 여섯 묶음으로 나누면 누락을 줄일 수 있다.

  1. 계측기와 교정: 표준가스 인증값, 영점과 스팬, 선형성, 분해능, 반복성, 장기 드리프트, 교차감응, 온도·습도·압력 영향, 샘플 라인의 누설과 흡착이다.

  2. 공간 표본: 센서가 실제 유출 공기를 대표하는지, 급이구·분뇨구·팬 주변의 국소 농도에 치우쳤는지, 유입·배출 지점 구분이 맞는지다.

  3. 시간 표본: 반추·급이·분뇨 교반·환기 제어의 주기를 충분히 포착했는지, 측정하지 않은 야간·주말·계절을 어떻게 다뤘는지다.

  4. 활동자료: 일별 실제 사육두수, 체중, 섭취량, 사료 조성, 우유·증체량, 분뇨 관리기간과 장비 가동시간의 오차다.

  5. 모델과 변환: 환기량 추정, 배경농도 차감, 건식·습식 기준, 표준상태 환산, 시간 적분, 결측 대체와 외삽에서 생기는 불확도다.

  6. 기준선과 비교가능성: 기상, 축군 구성, 생산성, 사료 외 관리 변화가 기준선과 사업기간 사이에 남기는 차이다.

여기서 반복 측정의 통계로 구하는 성분은 보통 Type A 평가에 가깝고, 교정성적서·제조사 정보·과거 연구·전문적 판단에서 얻는 성분은 Type B 평가에 가깝다. Type A가 ‘좋은 불확도’, Type B가 ‘나쁜 추정’이라는 뜻은 아니다. 둘 다 표준불확도라는 공통 형태로 바꿔 결합해야 하며, 자료 출처와 가정을 투명하게 남기는 것이 핵심이다.

모든 퍼센트를 단순히 더하면 안 된다

각 입력의 불확도를 최종 결과 단위로 옮기려면 민감도계수가 필요하다. 이는 입력값이 조금 변할 때 결과가 얼마나 변하는지를 뜻한다. 독립적인 입력들이고 모델이 근사적으로 선형이라면 각 표준불확도에 민감도계수를 곱한 기여도를 제곱해 더한 뒤 제곱근을 취한다. NIST가 제시하는 전형적인 불확도 예산도 구성요소, 민감도계수, 표준편차와 자유도를 함께 표로 관리한다.

그러나 축산 데이터에는 독립성 가정이 자주 깨진다. 같은 센서가 기준기간과 사업기간을 모두 측정했다면 교정 편향은 두 값에 공통으로 작용한다. 배출률 계산에 쓰인 CO₂ 발생량과 사료섭취량이 같은 활동자료에서 유도됐다면 상관될 수 있다. 여러 센서가 하나의 온습도 보정값이나 동일 표준가스를 공유할 수도 있다. 공분산을 무시하면 불확도를 과대 또는 과소평가할 수 있다. ‘보수적일 것’이라며 무조건 더하는 것도, 같은 장비이니 상쇄될 것이라 가정하는 것도 근거가 필요하다.

비선형 모델, 비대칭 분포, 검출한계 부근 자료, 비율의 분모가 작은 경우에는 1차 오차전파가 불안정할 수 있다. 이때는 입력별 확률분포와 상관구조를 정의한 몬테카를로 시뮬레이션이 유용하다. 다만 시뮬레이션 횟수를 늘리는 것이 잘못된 분포 가정이나 대표성 부족을 해결하지는 않는다. 모델 검증과 입력 근거가 먼저다.

최종 보고에서 ±만 쓰면 해석할 수 없다. 결합표준불확도인지, 포함인자를 곱한 확장불확도인지, 포함확률이 얼마인지 밝혀야 한다. NIST는 확장불확도 U = k × u_c로 설명하며 정규분포와 충분한 조건 아래 k=2가 대략 95% 구간에 해당할 수 있다고 안내한다. 그러나 모든 데이터에서 자동으로 95%가 되는 것은 아니다. 분포, 유효 자유도와 채택한 방법을 함께 제시해야 한다.

농장 현장에서 작성하는 불확도 예산의 예

가령 저메탄 사료 전후의 한 달 배출량을 비교한다고 하자. 농도계 교정 불확도는 2%지만, 자연환기량 추정은 18%, 출입구 간 농도차의 공간 표본은 10%, 결측 대체는 6%, 기준선의 일간·주간 변동은 12%라고 평가됐다. 이 수치를 단순 합산해 48%라고 쓰면 안 된다. 각 항목이 어느 결과에 작용하는지, 서로 독립적인지, 전후 비교에서 공통편향이 상쇄되는지부터 확인해야 한다. 또한 예시의 값들은 같은 결과량에 대한 표준상대불확도로 환산된 뒤에만 비교·결합할 수 있으며, 실제 평가는 분포·포함인자·상관관계에 따라 달라진다.

민감도 분석 결과 환기량이 전체 분산의 절반 이상을 차지한다면 더 비싼 농도센서로 교체하는 것은 우선순위가 아니다. 추적가스 실험을 늘리거나 개구부 기류를 추가 측정하고, 환기 운전구간을 층화하는 편이 효과적일 수 있다. 반대로 기준선 변이가 지배적이면 측정기간 연장, 동시 대조군, 사육두수·섭취량 정규화가 먼저다. 불확도 예산의 가치는 최종 구간을 만드는 데 그치지 않고 다음 측정 투자의 순서를 정하는 데 있다.

감축량이 100 kg CH₄이고 확장불확도가 30 kg이라면 ‘100 kg 감축’만 전면에 쓰지 않는다. 100 kg CH₄, 확장불확도 30 kg, 적용 기간·경계·포함확률·방법 명시처럼 보고한다. 인정 기준이나 크레딧 발급 문턱과 비교할 때는 해당 제도의 보수성 규칙을 별도로 적용한다. 통계적 불확도 계산이 곧 제도상 할인율이나 검증 판정을 자동 결정하는 것은 아니다.

운영 가능한 예산으로 만드는 관리 기준

불확도 예산은 보고서 말미의 일회성 표가 아니라 데이터 파이프라인과 연결돼야 한다. 센서 교정성적서가 갱신되면 해당 기간의 Type B 성분도 바뀌어야 한다. 결측률이 사전 한계를 넘으면 정해 둔 대체모델을 적용하거나 그 기간을 제외하고, 변경 이유를 기록한다. 사육두수와 사료 배합이 바뀌면 기준선 적합성 검토를 다시 수행한다.

버전 관리에는 최소한 원자료 ID, 센서·표준가스 ID, 교정일, 계산코드 버전, 배출계수·보정계수 출처와 버전, 결측 처리 규칙, 담당자와 승인 이력이 필요하다. 원자료를 덮어쓰지 않고 보정값과 산출값을 분리하면 방법론이 바뀌어도 재계산할 수 있다. 불확도 항목마다 책임자와 재평가 주기를 배정해야 ‘표는 있는데 아무도 갱신하지 않는’ 상태를 피할 수 있다.

IPCC 국가 인벤토리 지침은 불확도 분석을 숫자 하나의 완벽성을 증명하는 절차가 아니라 주요 불확도원을 우선 개선하고 시간에 따른 인벤토리 품질을 높이는 과정으로 다룬다. 농장 단위 프로젝트의 목적과 국가 인벤토리는 같지 않지만, 모든 중요한 원천을 식별하고 일관된 시계열을 유지하며 QA/QC와 연결한다는 원칙은 그대로 유효하다.

실행 체크리스트

  • 측정량을 가스, 경계, 축군, 기간, 단위, 기준상태까지 한 문장으로 정의했는가?

  • 최종 감축량을 모든 입력과 보정의 함수로 표현하고 모델 버전을 고정했는가?

  • 센서뿐 아니라 공간·시간 표본, 환기량, 활동자료, 결측, 기준선 불확도를 포함했는가?

  • 알려진 오류와 수정 가능한 편향을 먼저 교정했는가?

  • 각 성분의 Type A·Type B 근거, 분포, 표준불확도와 자유도를 기록했는가?

  • 공통 센서·표준가스·활동자료에서 생기는 상관관계와 공분산을 검토했는가?

  • 민감도계수로 최종 결과에 대한 기여도를 계산했는가?

  • 비선형·비대칭·검출한계 문제에 맞는 전파 방법을 선택했는가?

  • 결합표준불확도와 확장불확도, 포함인자와 포함확률을 구분해 보고했는가?

  • 가장 큰 기여항목을 다음 측정 개선계획과 예산에 연결했는가?

결론: 좋은 불확도 예산은 숫자의 약점을 숨기지 않는다

탄소감축 측정에서 불확도는 실패의 표시가 아니다. 실제 세계의 변동과 정보의 한계를 정량적으로 드러내는 방식이다. 중요한 것은 임의로 넓은 구간을 붙이는 일이 아니라, 측정대상과 모델을 고정하고 각 불확도원이 결과에 들어오는 경로를 추적하는 것이다.

센서 정확도만 관리하면 정밀해 보이는 잘못된 감축량을 만들 수 있다. 반대로 공간 대표성, 환기량, 활동자료, 기준선과 상관관계까지 포함한 예산은 어떤 주장을 할 수 있고 무엇을 더 측정해야 하는지 알려준다. 결국 불확도 예산은 탄소성과를 약화하는 부록이 아니라, 현장 기록이 검증 가능한 주장으로 이동하기 위한 품질 설계도다.

출처