해양 구조물 고장확률 효율 추정을 위한 순차 설계

초록

본 논문은 해양 구조물의 고장확률을 추정하기 위해 환경 조건의 조건부 밀도(CDE)를 효율적으로 추정하는 두 가지 방법을 제안한다. 첫 번째는 병렬 템퍼링 MCMC와 중요표본추출(IS)을 결합한 IS‑PT 방법이며, 두 번째는 적응형 가우시안 에뮬레이션과 베이지안 사분법을 결합한 AGE 방법이다. 두 방법을 다중모드 CDE를 갖는 합성 구조와 비선형 공진 응답을 보이는 모노파일에 적용해 비교했으며, λ 파라미터가 알려진 경우 AGE가 가장 적은 계산 비용을 보이고, λ가 미지인 경우 IS‑PT가 더 신뢰할 수 있음을 확인하였다.

상세 분석

논문은 해양 구조물의 장기 환경 변수 X와 단기 환경 변수 Y 사이의 확률적 관계를 명시적으로 모델링하고, 구조 응답 R = g_R(Y)와 조건부 밀도 f_{R|X}(r|x)를 통해 고장확률 p = P(R > r_{cr})를 정의한다. 핵심은 적분식 (3)에서 나타나는 “조건부 환경 밀도”(CDE) ˜f_X(x;r_{cr}) = P(R > r_{cr}|X = x)·f_X(x) 를 정확히 추정하는 것이다. CDE는 고장에 크게 기여하는 X 영역을 강조하므로, 샘플링 비용을 최소화하려면 이 분포에 비례하는 제안밀도 g_{pr}(x)를 사용해야 한다.

첫 번째 접근법인 IS‑PT는 CDE 자체를 목표 분포로 하는 병렬 템퍼링 MCMC를 수행한다. 병렬 템퍼링은 다중 모드 혹은 비타원형 형태의 CDE에서도 효율적인 탐색을 가능하게 하며, 체인 간 스와핑을 통해 낮은 온도 체인이 고밀도 영역을 오가게 한다. MCMC 샘플을 커널 밀도 추정으로 부드럽게 하여 g_{pr}(x)를 구성하고, 이를 기반으로 중요표본추출을 수행해 고장확률을 추정한다. 이 과정은 CDE를 직접 추정하므로 제안밀도와 목표밀도의 차이가 최소화되어 분산이 크게 감소한다.

두 번째 접근법인 AGE는 로그 CDE w(x)=log