본 논문에서는 LNG 벙커링 바지에 대한 예인안정성을 평가하였다. 친환경 에너지원인 LNG(액화천연가스)의 전환을 위한 인프라 로 LNG 벙커링 바지가 개발되고 있다. LNG 벙커링 바지는 예인줄에 연결되어 이동하는 부선의 형태이기 때문에 LNG 공급용 부선의 예인 안정성의 확보는 LNG 벙커링 바지 뿐 만아니라 주변 항해 선박의 안전을 위해 매우 중요하다. LNG 벙커링 바지의 예인안정성을 초기설계 단계에서 추정하기 위해 예인시뮬레이션을 위한 수치코드를 개발하였다. 초기설계 단계에서 활용 가능하도록 운동방정식에는 조종수학그룹 모델이 적용되었고 조종미계수에는 경험식이 적용되었다. 개발된 수치코드의 유효성을 확인하기 위해 발표된 계산 및 모형시험 결과와 비교 하였다. 개발된 수치코드를 이용하여 LNG 벙커링 바지의 선미 스케그 면적과 예인 위치의 변화에 따른 예인시뮬레이션을 수행하였다. 수행 된 시뮬레이션을 통해 설계된 선미 스케그 면적의 적정성을 확인하였다.