• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.57-57
• /
• 2011

지진해일은 진행속도가 빠르고 파장이 길며 파형의 변화 없이 먼 거리를 진행 할 수 있어 주변지역은 물론 멀리 떨어진 지역에도 심한 범람피해를 야기시킨다. 지진해일의 일반적인 특징으로 장파와 단파가 합성되어 있고 먼 거리를 전파할 경우 분산효과의 역할이 중요하게 된다. 특히 우리나라의 동해안에 영향을 주는 지진해일은 단주기파 성분이 강하고 파장에 비해 먼 거리를 전파하기에 분산을 고려하는 선형 Boussinesq 방정식을 지배방정식으로 사용하는 것이 바람직하다. 하지만 지금까지의 지진해일 전파모의를 위한 모형은 선형 Boussinesq 방정식의 복잡한 계산과 계산시간이 길다는 단점 때문에 선형 천수방정식을 지배방정식으로 사용하고 분산효과는 수치분산을 이용하여 고려해왔다. 지진해일 해석 시 일반적으로 사용되어 오던 기존의 leap-frog 유한차분 모형(Imamura et al., 1988; 조용식, 1996)은 지배방정식으로 선형 천수방정식을 사용하고 파의 분산효과는 수치분산을 이용하여 고려하므로 정해진 시간 간격에 대해 수심에 따라 격자 간격을 적절히 선택해야 하는데 수심이 복잡하게 변하는 경우 격자간격 조정이 불가능하여 분산효과를 정도 높게 고려할 수 없다. 이 문제점을 해결하기 위하여 윤성범 등(2004)은 파동방정식의 인위적인 분산항을 이용하여 Boussinesq 방정식의 분산효과를 고려할 수 있는 능동적인 분산보정기법을 제안하였고 Cho et al.(2007)는 일정한 수심에서 수치적인 분산오차가 Boussinesq 방정식의 물리적인 분산항을 대체하도록 수심, 격자 간격 및 계산 시간 간격 사이의 관계식을 유도하고 Boussinesq 방정식의 분산항과 일치하는 수치분산을 이용하여 실용적인 분산보정기법을 개발하였다. 이에 Ahn(2010)은 현재 컴퓨터의 계산 능력이 향상되어 선형 Boussinesq 방정식을 직접 차분하여 계산하는데 무리가 없다고 판단하여 선형 Boussinesq 방정식을 직접 차분한 모형을 개발하였다. 본 연구에서는 기존의 원해 지진해일 전파모의에 이용되어왔던 선형 천수방정식에 수치분산을 고려한 모형 대신 선형 Boussinesq 방정식의 유한차분 모형을 제안하였으며 기존의 선형 Boussinesq 방정식 모형의 격자와 수심간의 제약을 없애기 위해 차분 기법을 달리 한 2차 정확도의 유한차분 모형을 제안하였다. 검증을 위하여 선형 Boussinesq 방정식의 해석해(Carrier, 1991)와 비교하였다.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.107.2-107.2
• /
• 2010

본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 착용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 계산 효율이 개선되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.59-62
• /
• 2012

지진해일은 한 번 발생하면 매우 큰 피해를 주는 자연재해 중 하나이다. 특히 최근에 일어나고 있는 지진해일에서 주목할 점은 해저지진의 규모가 커짐에 따라 그로인해 생긴 지진해일이 발생 인접지역에만 피해를 주는 것이 아니라 2010년 칠레지진과 2011년 동일본 대지진으로 생긴 지진해일과 같이 태평양을 가로질러 먼 거리에까지도 영향을 준다는 것이다. 하지만 이러한 지진해일에 대하여 피해 경감 및 재해대책을 수립할 때는 지진해일의 발생시간에 대한 정확한 예측이 불가능하기 때문에 연구하는데 여러 가지 어려움이 있으므로 역사 및 가상지진해일의 수치모의실험 결과를 이용한다. 그래서 이렇게 활용성이 높은 지진해일 수치모의실험 결과의 정확도 검증을 위하여 본 연구를 시작하게 되었다. 본 연구에서 사용한 수치기법으로는 분산항을 제외한 선형 천수방정식을 이용하였고, 이를 차분하는 과정에 나타나는 수치분산을 조정하여 선형 Boussinesq 방정식의 분산항을 대체할 수 있도록 하였다. 그리고 실제 동일본 대지진으로 발생한 지진해일을 전파모의하여 분산항을 고려한 수치모의실험 결과와 이 당시 태평양에서 수면변화양상을 기록한 DART buoy 관측값과의 비교를 통해 정확도를 확인해보았다.

• Journal of Korea Water Resources Association
• /
• v.38 no.6 s.155
• /
• pp.461-469
• /
• 2005

A second-order upwind scheme is used to investigate the run-up heights of tsunamis in the East Sea and the predicted results are compared with the field data and results of a first-order upwind scheme. In the numerical model, the governing equations solved by the finite difference scheme are the linear shallow-water equations in deep water and nonlinear shallow-water equations in shallow water. The target events are 1983 Central East Sea Tsunami and 1993 Hokkaido Nansei Oki Tsunami. The predicted results represent reasonably well the run-up heights of tsunamis in the East Sea. And, the results of simulation are used for the design of inundation map.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.2105-2108
• /
• 2007

지진해일이 발생하여 먼 거리를 전파할 때 코리올리스 힘이 지진해일의 거동에 매우 큰 영향을 준다. 또한 지진해일은 파장에 비해 파고가 작은 장파이기 때문에 먼 거리를 이동하는 지진해일의 해석에는 코리올리스 효과가 반영된 선형 천수방정식을 사용하는 것이 적절하다. 본 연구에서는 구면좌표계에서 유한차분기법을 사용하여 지진해일의 거동을 모의하였다. 대상 지진해일은 1983년 동해 중부 지진해일이다. 관측된 결과를 본 연구에서의 모의 결과와 비교하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.52-56
• /
• 2010

본 연구에서는 지진해일의 전파 과정을 모의함에 있어 선형 천수방정식의 수치분산을 이용하는 기법이 아닌 선형 Boussinesq 방정식을 직접 차분하는 유한차분기법을 제안하였다. 지배방적식과 차분식의 일치성을 해석하기 위해 이산화 오차를 확인하고, 수치해의 안정적 수렴여부를 판단하기 위해 Von neumann 안정성 해석을 수행하였다. 또한 기법의 정확성을 검증하기 위하여 Gauss 분포의 초기 자유수면변위를 갖는 문제에 적용하여 선형 Boussinesq 방정식의 해석해와 비교하였다. 그 결과 기존의 선형 천수방정식을 차분화한 수치모형에 비하여 정확한 결과를 제공하였고 분산보정기법을 이용한 수치모형과 동일한 정확도를 보였으나 본 수치모형을 이용했을 때 비교적 넓은 범위의 조건에서 정확도 높은 결과를 제공하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.40-41
• /
• 2011

동해를 전파하는 지진해일은 세계적으로 다른 지역에서 발생하는 지진해일과 비교하였을 때 상대적으로 파장이 짧고, 이에 비해 먼 거리를 전파한다. 그러므로 동해에서 발생한 지진해일의 전파에 대한 해석을 수행할 때 물리적인 분산효과가 매우 중요하다. 따라서 지배방정식으로 분산 효과가 충분히 고려된 선형 Boussinesq 방정식을 사용한다. 기존의 연구에서는 leap-frog 기법을 사용하여 선형 천수방정식을 차분할 때 발생하는 수치분산항에 분산 보정계수를 이용하여 선형 Boussinesq 방정식의 물리적 분산항과 같은 형태로 나타나도록 유도하여 수치모의를 수행하였다. 그러나 기존에 사용한 지배방정식은 수심이 일정하다는 가정을 통하여 유도된 것으로, 수심에 변화가 있는 실제 지형을 통과하는 지진해일에 대한 수치모의를 수행한 결과의 정확도에 문제가 생길 수 있다. 본 연구에서는 기존의 연구에서 발생할 수 있는 수심 변화에 따른 오류를 개선하기 위하여 바닥 지형이 1차원으로 변한다는 가정을 이용하여 지배방정식을 유도하였으며, 이로 인해 발생하는 수심 변화가 고려된 항을 기존의 분산보정기법에 추가하였다. 그리고 적용성을 높이기 위하여 수치모의 기법의 제한을 최소화하는 연구를 진행하였다. 본 연구에서 제안한 수정 기법이 수심이 변화하는 지형을 전파하는 지진해일 수치모의 과정에서 경사에 대한 분산효과가 충분히 고려되는지 확인하기 위하여 Gaussian hump를 이용한 가상 지진해일을 원형 천퇴 지형에 통과시켰다. 본 연구에서 사용한 지형을 통과하는 Gaussian hump에 대한 해석해를 구하는 방법이 존재하지 않으므로, Boussinesq 방정식을 직접 차분하여 푸는 FUNWAVE를 사용하여 동일한 조건 하에서 수치모의를 수행하였다. 비교 결과를 통하여 본 연구에서 제안한 기법의 정확도 향상을 확인하게 되면, 실제 지형을 통과하는 지진해일의 수치모의에 대한 활용성을 높일 수 있을 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.2098-2101
• /
• 2007

본 연구에서는 최근 5년간 서해에서 발생한 해저 지진 중에서 가장 큰 지진인 2003년 3월 23일에 발생한 해저 지진의 진앙지(북위 $35^{\circ}$, 동경 $124.6^{\circ}$)에 우리 나라에 가장 큰 피해를 입힌 1983년 동해 중부 지진해일의 지진 매개변수(parameter)를 적용한 가상 지진해일을 발생시키고 새만금 지역의 최대 처오름 높이를 계산하였다. 가상지진해일에 대한 수치모의 결과 새만금 지역에서는 $0.4{\sim}0.9m$의 파고가 발생하겠으며, 이러한 결과는 지진해일로 인한 새만금 방조제 유실 및 방조제 내부의 재난방재를 위한 기초적 자료가 될 것이다.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.58-58
• /
• 2011

동해를 전파하는 지진해일은 다른 지역에서 발생하는 지진해일과 비교하였을 때 상대적으로 파장이 짧고, 파장에 비해 먼거리를 전파한다. 따라서 지진해일이 전파할 시 물리적인 분산효과가 매우 중요하다. 그러므로 동해에서 발생하는 지진해일을 수치모의 할 때는 분산효과가 충분히 고려될 수 있는 선형 Boussinesq 방정식을 사용한다. 그러나 이를 직접 풀 경우에는 상당히 많은 시간이 소비되며 효율적이지 못하다. 이와 같은 단점을 극복하기 위해 기존의 연구에서는 leap-frog 기법을 사용하여 선형 천수방정식을 차분할 때 발생하는 수치분산항을 선형 Boussinesq 방정식의 분산항과 같은 형태를 가질 수 있도록 분산보정계수를 사용하여 수치모의를 수행하였다. 하지만 이때 사용된 지배방정식은 수심이 일정하다는 가정을 이용하여 유도된 것이므로, 실제 경사가 있는 지형을 통과할 때의 수치모의 결과는 정확하다고 할 수 없다. 본 연구에서는 이를 극복하기 위하여 바닥 지형이 1차원으로 변한다는 가정으로 새로운 지배방정식을 유도하였으며, 수심변화로 인해 새로 발생하는 항을 기존의 분산보정기법에 추가하였다. 또한 수심이 변화는 지형을 통과하는 지진해일의 분산효과가 충분히 고려되는지 확인하기 위하여 Gaussian hump를 이용하여 가상 지진해일을 원형 천퇴지형에 통과시키는 수치모의를 수행하였다. 결과의 비교를 위한 정확해가 없으므로, 비선형 Boussinesq 방정식을 직접 차분하여 푸는 FUNWAVE를 이용하여 동일한 조건으로 수치모의를 수행하였다. 수치모의 시 중심선에 4개의 가상 gage를 설치하였으며, 이를 통해 각각의 수치모의 실험에 대한 자유수면 변위를 관찰하여 비교하였다. 수치모의 결과에 대한 비교를 통하여 기존의 분산보정기법에 비해 본 연구에서 제안한 새로운 수치기법이 분산효과를 비교적 잘 반영하는 것으로 나타났으며, 비교적으로 실제 지형에 적용하였을 때 정확도 향상의 가능성이 높다고 판단하였다.

• Journal of Korea Water Resources Association
• /
• v.35 no.4 s.129
• /
• pp.453-461
• /
• 2002

The run-up process of the 1983 Central East Sea Tsunami along the Eastern Coast is numerically investigated in this study. A finite difference numerical model based on the nonlinear shallow-water equations is employed. The maximum run-up height at Imwon is predicted and compared to field observation. A good agreement is observed. A maximum inundation map is made based on the maximum run-up heights to accentuate hazards of tsunami flooding.