심해에서 발생하여 천해로 진입하는 파는 다양한 요인들에 의해서 변형되며 이 중에서 파의 진행 방향에 큰 영향을 끼치는 것들이 회절과 굴절이다. 회절은 에너지 차이에 의해 발생하는 현상으로 섬, 부유 구조물 또는 방파제 주변에서만 발생한다. 반면, 굴절이란 파봉선 방향의 위상속도 차이에 의해 파의 진행 방향이 바뀌는 현상으로 전 영역에 걸쳐서 발생한다. 바닥이 불투수층인 경우에 에너지 손실을 무시하면 수심 또는 흐름의 공간적인 변화로 위상속도의 차이가 발생하여 굴절이 발생한다. 파랑이 투수층을 지나는 경우, 흐름의 공간적인 변화뿐만 아니라 항력 저항, 관성 저항 등의 에너지 손실로 위상속도의 변화가 발생하면서 굴절이 발생한다. 이러한 요인들로 인하여 굴절이 발생하여 파가 특정 부분에 집중되면 높은 파 에너지로 인하여 쇄파나 쇄파유도류가 발생하여 연안환경에 영향을 끼치며 또는 연안구조물의 붕괴를 초래하기도 한다.

파의 굴절을 계산하는 방법으로 등수심선이 나란한 경우 Snell의 법칙을 이용할 수 있으며, 수심이 복잡한 경우에는 위상속도를 이용하여 진행 방향을 예측하는 파향선법을 사용할 수 있다(Munk and Arthur, 1952; Wilson, 1966; Noda, 1974). 파향선법은 파의 진행거리 및 방향이 미분식으로 표현되어서 수치차분을 이용하여 파의 이동거리와 방향을 구할 수 있다. 이 방법을 이용하면 주어진 파랑 및 지형조건에서 파랑의 진행상황을 예측할 수 있지만 입사파의 주기, 항력계수 또는 관성계수 등의 조건이 변하였을 경우 굴절이 어떻게 바뀔지는 미리 계산해보기 전에 알 수가 없다.

본 연구에서는 굴절이 파봉선 방향으로의 위상속도 차이에 의해서 발생한다는 사실에 착안하여 다양한 조건에서 발생하는 파의 굴절 정도를 해석적으로 표현하는 방법을 제시하였다. 유도된 식은 불투수층 및 투수층에 대하여 적용하였으며 불투수층의 경우 기존에 알려진 사실과 잘 부합하는 것을 확인할 수 있었다.

파의 굴절은 나란히 진행하는 파와의 진행속도 차이에 의해 발생한다. 따라서 굴절의 정도를 나타내는 굴절수(index of refraction degree)는 다음과 같이 파봉선 방향으로의 위상속도차로 정의할 수 있다.

여기서, c는 위상속도, k는 파수, ω는 각 주파수이며 n은 Fig. 1과 같이 파 진행 방향의 법선(즉, 파봉선 방향) 벡터를 의미한다. 맨 앞에 붙은 1/ω는 ∂c/∂n이 [T⁻¹]의 차원을 갖기 때문에 무차원화시키기 위하여 사용되었다.

본 연구에서 제안한 해석적인 방법의 타당성을 검토하기 위하여 먼저 바닥이 불투수층인 지형에 적용해 보았다.

Fig. 5는 수심값에 변화를 주면서 식(6)과 (7)을 계산하여 도시한 결과이다. 식과 그림을 통해 알 수 있듯이 불투수층일 경우에는 굴절의 정도가 오직 상대수심에 의해서만 결정이 되며 상대수심이 작을수록 굴절이 크게 발생하는 것을 예측할 수 있다. 식(7)로 표현되는 근사식은 천해역에서는 정해와 잘 일치하지만 심해역에서는 다소 큰 값을 보인다. 이와 같은 결론은 Snell의 법칙이나 파향선법을 이용하여 파의 진행경로를 계산함으로써 확인할 수 있다. Fig. 6은 서로 다른 상대수심을 가지는 파들이 동일한 지점에서 동일한 입사각 θ = 30°으로 입사하는 경우에 대하여 Snell의 법칙을 이용하여 파의 진행 방향을 그린 것이다. 얕은 수심으로 갈수록 상대수심의 효과로 인하여 파의 굴절이 크게 발생하는 것을 확인할 수 있다. 이는 Lee et al.(2003)이 발견한 것과 동일하다.

Fig. 7은 서로 다른 공극률(λ)을 갖는 투수층에 대하여 굴절 정도를 계산한 그림이다. 투수층에서의 변수값은 a_l = 150, ν = 1.15 × 10⁻⁶(m²/s), d = 0.01 (m) 및 κ = 0.7으로 고정하고 공극률은 λ = 0.35, 0.45, 0.55의 값을 사용하였다. 투수층에서의 굴절 정도는 전체적으로 물 층인 경우와 비슷한 경향을 보였다. 상대수심이 작은 경우에는 굴절이 크게 발생하였으나 상대수심이 증가할수록 굴절의 정도가 약해지면서 상대 수심이 일정한 값 이상이 되면 굴절이 거의 발생하지 않았다. 공극률이 클수록 굴절이 크게 발생하였으나 그 차이가 크지는 않았다.

Fig. 8은 Fig. 7과 동일한 수심 조건에서 파의 진행을 계산한 결과이다. 입사파의 각도는 30°이며, 굴절이 크게 발생하도록 장주기(29.02초) 조건을 사용하였다. Fig. 7의 결과에서 보여지듯이 공극률에 따른 굴절 정도의 차이가 미미하기 때문에 실제 파의 진행경로도 거의 동일하게 나타났다.

마지막으로 투수층 위에 물 층이 있는 경우에 대하여 굴절 정도를 계산해 보았다. Fig. 9는 식(20) 및 (21)의 ∂c/∂h₁/ω과 ∂c/∂h₂/ω를 계산한 결과이다. 이는 투수층의 수심을 일정하게 유지하면서 물 층의 수심이 변하는 경우와 그 반대의 경우를 의미하는 것으로써 큰 값을 가질수록 굴절이 크게 발생한다는 것을 의미한다. 투수층에서의 조건은 Fig. 8에서 사용된 조건과 동일하게 사용하였다.

수심 변화가 위의 물 층이나 아래 투수층에서 동일하게 발생할 경우 물 층에서의 수심변화가 굴절에 더 큰 영향을 주었다. 이는 Fig. 10과 같은 지형에 파향선을 그려봄으로써 확인할 수 있다. Fig. 10(a)는 아래 투수층의 두께는 h₂으로 고정한 후, 위 물 층의 두께를 h에서 0까지 감소하는 지형이며 Fig. 10(b)는 반대되는 경우에 대한 지형이다.

Fig. 11은 h = 20 m, θ₀ = 30°의 조건에서 Fig. 10의 지형에 파향선을 그린 결과이다. 앞서 기술한대로, 물 층에서 수심변화가 발생한 경우에 투수층에서 수심변화가 발생한 경우보다 굴절이 더 크게 발생하였으며 이러한 경향은 상대수심이 작을수록 더 크게 발생하였다.

본 연구에서는 파의 굴절 정도를 해석적으로 예측할 수 있는 방법을 제시하였다. 파봉선 방향으로의 위상속도차에 의해 굴절이 발생한다는 원리에 착안하여 식을 유도하였다. 본 연구에서 제시한 방법은 파의 분산관계식을 이용하기 때문에 불투수층이나 투수층 또는 진흙층의 경우에도 분산관계식만 구할 수 있으면 적용할 수 있다. 먼저 불투수층 지형에 적용하여 기존의 결과와 잘 일치하는 것을 확인할 수 있었으며 투수층 및 2개 투수층에 대해서도 굴절정도와 파향선의 결과가 일관되게 나왔다.