불확실성 정량화

시뮬레이션 통계 수치해석 마지막 수정: 2026-07-18 04:11:07

1. 개요[편집]

불확실성 정량화(Uncertainty Quantification, UQ)는 시뮬레이션의 입력·모델·수치오차에 들어있는 불확실성이 관심 있는 출력에 얼마나 전파되는지를 정량적으로 다루는 분야다. 한마디로 “이 해석 결과, 숫자 하나로 믿어도 되나요?”라는 질문에 확률로 답하는 학문이다. 결정론적 해석이 점 하나( $y = f(x)$ )를 내놓는다면, UQ는 그 점 대신 분포를 내놓는다.

엔지니어링에서 이게 중요한 이유는 단순하다. 재료 물성, 경계 조건, 하중, 제조 공차 중 정확히 아는 값은 하나도 없다. 그런데도 우리는 대표값 하나를 넣고 유한요소법이나 전산유체역학을 돌린 뒤 그 결과에 도장을 찍는다. UQ는 그 관행에 “입력이 흔들리는데 출력만 안 흔들릴 리가 없잖아”라고 태클을 건다.¹

2. 불확실성의 두 종류[편집]

UQ의 출발점은 불확실성을 성질에 따라 둘로 쪼개는 것이다. 이 구분을 놓치면 대화가 계속 어긋난다.

우연적 불확실성(aleatory) — 세상 자체의 본질적 무작위성. 난류의 순간 변동, 재료 결정립의 산포, 풍하중의 시간 변동처럼 줄일 수 없는(irreducible) 변동이다. 데이터를 아무리 모아도 사라지지 않고, 확률분포로 기술한다.
인식적 불확실성(epistemic) — 우리가 몰라서 생기는 불확실성. 모델을 잘못 골랐거나, 파라미터를 덜 측정했거나, 격자가 거칠어서 생긴다. 이건 데이터·연구·격자 세분화로 줄일 수 있다(reducible). 확률로 다룰지 구간(interval)으로 다룰지부터 논쟁거리다.

라틴어까지 동원해 굳이 나누는 이유가 있다. 우연적 불확실성은 안전계수로 감내할 대상이지만, 인식적 불확실성은 돈과 시간을 쓰면 줄어드는 대상이라 의사결정이 완전히 다르기 때문이다. 격자를 촘촘히 해서 사라지는 오차에 안전계수를 붙이는 건 낭비고, 재료의 본질적 산포를 연구비로 없애려는 건 삽질이다.

3. 전방 문제와 역문제[편집]

UQ가 푸는 문제는 방향에 따라 둘로 갈린다.

**전방 UQ(forward propagation)**는 입력 분포가 주어졌을 때 출력 분포를 구한다. 입력 확률변수 $\mathbf{X}$ 의 분포와 모델 $Y = f(\mathbf{X})$ 를 알 때, $Y$ 의 평균·분산·분위수·파손확률 $P(Y > y_{\text{crit}})$ 를 추정하는 것이 목표다. UQ라고 하면 보통 이쪽을 떠올린다.

**역 UQ(inverse problem)**는 반대로, 관측된 출력으로부터 입력 분포를 추정한다. 실험 데이터에 모델을 맞춰 미지 파라미터의 사후분포를 구하는 **베이지안 보정(Bayesian calibration)**이 대표적이다. 이건 사실상 역문제이자 베이지안 최적화와 같은 통계적 뿌리를 공유한다. 관측이 불완전하고 노이즈가 섞여 있으므로 결과 역시 분포로 나온다.

전방이든 역이든 최종적으로 궁금한 건 대개 하나로 수렴한다. 이 설계가 규정을 넘길 확률이 몇 퍼센트인가.

4. 전파 기법[편집]

입력 분포를 출력 분포로 옮기는 실전 도구들이다.

몬테카를로 방법 — 입력 분포에서 표본을 뽑아 모델을 반복 평가한다. 차원에 무관하게 수렴률이 $O(N^{-1/2})$ 로 일정한 게 최대 장점이자, 유일한 단점이다. 오차를 10배 줄이려면 표본을 100배 늘려야 하는데, 해석 한 번에 하루 걸리면 게임이 끝난다.
준몬테카를로(QMC) — 소볼·할톤 같은 저불일치 수열로 표본을 균등하게 깔아 수렴을 $O(N^{-1})$ 근처까지 끌어올린다. 매끄러운 저차원 문제에서 공짜로 얻는 이득.
다항식 카오스 전개(gPC) — 출력을 입력 확률변수의 직교다항식 급수로 전개한다. 가우시안 입력엔 에르미트, 균등 입력엔 르장드르 다항식을 쓰는 식으로 입력 분포에 맞춰 기저를 고른다(Wiener–Askey 체계).

Y = f(\mathbf{X}) \approx \sum_{i=0}^{P} c_i \, \Psi_i(\mathbf{X})

계수 $c_i$ 만 구하면 평균·분산은 물론 소볼 지수까지 급수에서 해석적으로 뽑힌다. 매끄러운 응답에서는 몬테카를로보다 수백 배 빠르지만, 입력 차원이 늘면 항의 수가 폭발하는 차원의 저주에 걸린다.

대리모델 기반 UQ — 대리 모델(크리깅 등)을 먼저 학습시켜 놓고 그 위에서 몬테카를로를 10만 번 돌린다. 원본 해석 대신 값싼 근사를 두들기는 것이라, 오늘날 비싼 해석에 대한 UQ는 사실상 전부 이 방식이다. gPC도 넓게 보면 대리모델의 일종이다.

5. 검증·확인과의 관계[편집]

UQ는 검증 및 확인(V&V)의 자연스러운 종착지다. V&V가 “코드를 맞게 풀었나(verification), 맞는 방정식을 풀었나(validation)“를 따진다면, UQ는 거기서 남는 수치오차와 모델 불확실성에 **정량적 크기표(error bar)**를 붙인다. 그래서 요즘은 아예 붙여서 VVUQ라 부른다.

특히 중요한 게 **모델 형태 불확실성(model-form uncertainty)**이다. 파라미터가 아니라 지배방정식이나 난류 모델 선택 자체가 틀렸을 가능성인데, 이건 확률분포로 깔끔하게 담기지 않아 UQ에서 가장 다루기 까다로운 짐승이다. 여러 모델을 앙상블로 돌려 그 산포를 보는 정도가 현실적 타협이다.

미국 원자력규제위원회(NRC)나 항공·국방 인증에서 UQ가 규제 요건으로 자리 잡은 것도 이 맥락이다. “우리 시뮬레이션이 안전합니다”가 아니라 “95% 신뢰수준에서 파손확률이 $10^{-6}$ 이하입니다”라고 말해야 서류가 통과된다.²

6. 현업에서의 현실[편집]

입력 분포를 정하는 게 UQ의 8할이다. “이 물성의 표준편차가 얼마죠?”라고 물으면 대부분 모른다. 결국 전문가 추정이나 빈약한 데이터로 분포를 지어내는데, 쓰레기를 넣으면 쓰레기 분포가 나온다(GIGO).
몬테카를로 표본 수를 물으면 “많이”라고밖에 답할 수 없어 보이지만, 추정하려는 게 평균이냐 $10^{-4}$ 짜리 꼬리확률이냐에 따라 필요 표본이 수천 배 차이 난다. 희소사건은 중요도 표본추출(importance sampling)로 따로 공략한다.
민감도 해석과 UQ는 붙어 다닌다. 민감도 해석으로 영향 없는 입력을 먼저 쳐내면 UQ 차원이 줄어 계산이 산다. 실무에선 스크리닝 → 차원 축소 → 전파 순서가 국룰이다.³
결과를 평균±표준편차 하나로 보고하면 절반만 전달된 것이다. 의사결정자가 진짜 봐야 하는 건 분포의 꼬리, 즉 최악의 경우가 얼마나 나쁜가다.

7. 관련 문서[편집]

그래서 나온 유명한 격언이 “모든 모델은 틀렸지만 일부는 유용하다”(조지 박스). UQ는 여기에 “얼마나 틀렸는지를 숫자로 말하자”를 덧붙인 셈이다. ↩
이 꼬리확률 게임의 무서운 점은, 관심 대상이 $10^{-6}$ 인데 몬테카를로를 백만 번 돌려봐야 사건이 평균 한 번 일어난다는 것이다. 그래서 희소사건 UQ는 그 자체로 하나의 전문 분야가 됐다. ↩
우연적/인식적 구분을 실무에서 진지하게 지키는 조직은 생각보다 적다. 다 섞어서 몬테카를로 한 방에 돌려버리는 게 흔한데, 그러면 “격자를 촘촘히 하면 줄어들 오차”와 “본질적 산포”가 한 분포에 뭉개져 의사결정의 근거가 흐려진다. ↩