I. 서 론
천해 해양환경에서는 심해에 비해 해저면과 해수면의 영향을 많이 받으며 복잡한 양상으로 음파전달이 이루어진다. 다층구조 퇴적층으로 구성된 환경에서 발생하는 하부퇴적층에 의한 다중 경로 특성은 수신 신호간의 간섭현상을 더욱 복잡하게 만드는 요인이 될 수 있기 때문에 하부 퇴적층의 층 구조를 파악하는 것이 상당히 중요하다. 특히 중주파수 및 저주파수 대역의 소나를 운용할 때에는 고주파수 대역에 비해 해저면 하부로 투과되는 음파 에너지의 양이 증가하기 때문에 하부퇴적층에 대하여 충분하게 고려해야한다. 해양 퇴적층 구조를 파악하기 위해 퇴적 샘플을 직접 탐사하는 방법은 인적/물적 비용이 크게 소모되는 단점이 있기 때문에, 음향을 이용한 원격 탐사 방법으로 퇴적층 구조를 역 추정하는 연구가 필요하다. 따라서 천해에서 수중 음향을 이용한 퇴적층 구조 및 물리적 특성 역산 연구는 국내·외에서 다양한 방법으로 연구가 진행되었다. 퇴적층 구조에 대한 역산 연구로는 하부 퇴적층으로부터 반사되어 수신한 신호를 이용한 하부 퇴적층의 두께 및 음속, 감쇄계수 역산,[1] 다층 구조 퇴적층에서 수신된 반사 손실을 이용한 퇴적층의 두께 및 지음향 인자 역산,[2] 해양에 존재하는 주변 소음의 수신 위상 정보를 이용한 하부 퇴적층 구조 역산[3] 등이 있다. 특히 하부 퇴적층에 의해 복합적으로 나타나는 해저면 반사 신호를 이용한 음향학적 특성 분석을 통해 퇴적층 구조 및 물성에 대한 역산을 수행한 최근 연구 사례가 다수 존재한다.[4-6] Holland와 Osler[4]는 1 ms 이내의 매우 짧은 펄스길이의 광대역 신호를 이용하여 하부 퇴적층을 고려한 반사 손실을 측정하였으며, 이를 이용하여 지음향 인자를 역산하여 퇴적층 코어 정보와 비교하였다. 이와 유사한 방법으로 Dettmer와 Dosso[5]는 하부 퇴적층 반사 신호에서 나타나는 위상 차이를 이용하여 반사 계수를 계산하였으며, 이를 모델 결과와 비교하여 지음향 인자를 역산하였다. 이와 같이 하부 퇴적층에 의한 수층에서의 음파 전달 양상 차이를 분석하여 퇴적층 구조 및 지음향 인자를 역산하는 연구가 최근에 활발하게 수행되고 있다.
본 논문에서는 수평입사각에 따른 중주파수 대역(4~8 kHz) 해저면 반사 신호 측정 실험에서 수신한 표층 해저면 반사 신호와 하부퇴적층 반사 신호의 도달 시간차이를 이용하여 하부 퇴적층의 두께를 역으로 추정하고자 하였다. 또한 추정된 하부퇴적층 두께를 한국해양과학기술원에서 제공받은 실험 해역의 천부지층탐사기(sub-bottom profiler, CAP-6600, benthos) 탐사 결과와 비교하여 그 가능성을 확인하였다.
II. 해상실험 및 환경
수평입사각에 따른 해저면 반사 신호 측정 실험은 2015년 10월, 경남 거제시 장승포항 7 km 남동쪽 해역(34° 48.906‘ N, 128° 47.34’ E)에서 수행되었다(Fig. 1). 실험 해역의 수심은 약 58 m이며, 송신기와 수신기를 2대의 조사선에서 각각 분리하여 설치하였다. 수신 조사선에 총 2개의 수신기(Ch.1: TC-4014, Ch.2: TC-4032, Reson)를 약 35, 45 m 수심에 수직 배열하였으며, 송신 조사선의 음원(D/11, Neptune)은 수심 약 45 m에 위치하였다. 또한 각 수신기와 음원의 상단에 심도기록계(DR-1050, RBR)를 설치하여 실험하는 동안의 송·수신기의 수심 변화를 측정하였다. 송·수신기 간의 거리는 각각 120, 150 m로 고정하였으며, 두 조사선을 밧줄로 연결하여 거리 정확도를 향상시켰다. 궁극적으로 9 ~ 14°의 수평입사각 범위에서 해저면 반사 신호를 측정하였으며[Fig. 2(a)], 실험이 진행되는 동안의 수직 음속구조는 CTD(Conductivity, Temperature, Depth, AML Oceanographic)를 이용하여 측정되었다[Fig. 2(b)]. 송신 신호는 주파수 4, 6, 8 kHz, 펄스길이 1 ms의 정현파 신호를 30핑 반복 송·수신하였으며, 수신 신호의 샘플링 주파수는 500 kHz이었다. 실험 해역의 해저면 구성성분은 한국지질자원연구원에서 제공받은 표층 퇴적물의 평균 입도 정보[7]를 사용하였으며, 실험 해역의 표층 퇴적물 평균 입도 범위는 8 ~ 10 ϕ이다.
하부 퇴적층 두께 역산을 위한 실험으로는 송·수신기 거리 120 m, 송신기 수심 45 m, 수신기 수심 45 m에서 수행한 실험을 실험 세트 1, 송·수신기 거리 150 m, 송신기 수심 45 m, 수신기 수심 35 m에서 수행한 실험을 실험 세트 2, 송·수신기 거리 150 m, 송신기 수심 45 m, 수신기 수심 45 m에서 수행한 실험을 실험 세트 3이라 정의하고, 이후 설명에서 실험 세트 1, 실험 세트 2, 실험 세트 3으로 구분한다.
송·수신기의 수평거리(
)는 각 조사선에 설치된 DGPS(Differential Global Positioning System, GPS850, ASCEN) 자료를 이용하여 계산하였다. 그리고 DGPS의 오차 범위인 2.5 m로 인해 발생할 수 있는 수평 입사각의 오차를 줄이기 위해, 각 실험 세트마다 수신된 실측 신호의 다중경로별 도달시간과 음선 이론이 기반인 음파전달모델로 모의한 신호의 도달시간을 비교하여 수평거리 
및 총 수심 오차를 보정하는 과정을 수행하였다. Fig. 3은 실험 세트 1의 주파수 4 kHz 측정 신호[Fig. 3(a)]와 음파전달모델을 이용하여 측정한 해양환경 정보를 입력하여 모의한 수신 신호[Fig. 3(b)]이다. 이와 같이 실측 신호와 모의된 신호의 각 다중 경로별 도달 시간을 비교함으로써 총 수심과 송·수신기간의 수평거리 정확도를 향상시켰다. 또한 모의된 신호와 달리 측정 신호에서는 해저면 반사 신호와 해수면 반사 신호의 사이에 추가적인 다중경로가 존재하는 것을 확인할 수 있다. 이 다중경로는 해저면 하부에 존재하는 퇴적층에 의한 반사 신호로 판단이 되며, 사용한 음파전달모델에서는 수층에서의 음파전달만 모의 가능하기 때문에 모의 결과에서는 하부퇴적층 반사 신호를 확인할 수 없었다. 본 논문에서는 수신 신호로부터 관찰된 하부 퇴적층에 의한 반사 신호와 표층 해저면 반사 신호의 도달 시간 차이를 이용하여 하부 퇴적층의 두께를 추정하는 것을 목적으로 한다.
III. 퇴적층 특성 역산 과정 및 방법
하부 퇴적층의 두께 
를 역으로 추정하기 위해서는 두 가지 다중 경로 전파에 대한 이해가 필요하다. 첫 번째로 표층 해저면에서 반사되어 수신된 경로(Bottom bounce path, B-path), 두 번째는 표층 해저면에서 투과한 이후에 하부 퇴적층의 경계면에서 다시 반사되어 수신된 경로(Sub-bottom bounce path, Sub- B-path)이다. 두 가지 다중경로의 전파에 대한 모식도는 Fig. 4와 같으며, 서로 다른 경로를 통해 수신되는 두 신호는 도달시간의 차이를 갖는다. 따라서 표층 퇴적층의 두께는 이 신호들의 도달시간 차이와 표층 퇴적층의 음속을 사용하여 추정할 수 있다.[6] 본 논문에서는 각각의 실험세트별로 하부 퇴적층의 두께와 음속을 변화시켜가면서 두 다중경로의 도달시간 차이 
를 모의하였고, 실험적으로 측정한 도달시간 차이 
와 비교하여 퇴적층 두께를 역산하였다.
: density, 
: sound speed (subscripts 1 refer to water and subscripts 2 refer to sediment), 
, 
: distance from the surficial sediment to transducer and receiver, respectively, 
: thickness of surficial sediment, 
: incidence grazing angle, 
: transmission grazing angle, 
: source-receiver range].표층 해저면과 하부 퇴적층 반사경로의 도달시간 차이를 모의하기 위해서, 우선 각 다중경로의 총 이동거리를 
, 
로 정의하였으며 Eqs. (1)과 (2)와 같다.
, (1)
, (2)
여기서 
은 송신기로부터 해저면까지의 수직 거리, 
는 각 수신기로부터 해저면까지의 수직 거리, 
는 하부 퇴적층 두께, 
는 해저면 수평 입사각, 
는 해저면 수평 굴절각이다(Fig. 4). 다중경로의 총 이동거리를 수층과 퇴적층에 해당하는 음속으로 각각 나누어주면 각 다중경로의 도달시간을 계산할 수 있으며, 이는 Eqs. (3)과 (4)로 표현 가능하다. 여기서 
은 수층에서 측정된 수직 음속 구조의 평균 음속이며, 
는 퇴적층 음속이다. 본 논문에서는 퇴적층의 음속구조(
)는 한국지질자원연구원[7]에서 제공받은 평균입도자료를 이용하여 계산되었다.
, (3)
, (4)
궁극적으로 표층 해저면 반사 신호와 하부 퇴적층 반사 신호의 도달시간 차이 
는 Eq. (5)로 정리할 수 있다.
. (5)
실험 해역의 해저면 구성성분은 한국지질자원연구원에서 제공받은 표층 퇴적물의 평균 입도 정보[7]와 남해 인근 해역의 표층 퇴적물 음속 정보[8]를 사용하였다. 실험 해역의 표층 퇴적물 평균입도는 8 ~ 10 ϕ, 음속 범위는 약 1500 ~ 1550 m/s를 나타냈으며 연안으로 갈수록 음속이 낮아지는 것을 확인했다.
본 실험과 같이 낮은 입사각으로 해저 경계면에 입사된 신호가 퇴적층 내로 투과되기 위해서는 퇴적층의 음속이 수층에 비해 낮아야 한다. 따라서 도달시간 차이를 모의하기 위해 설정된 환경 변수는 수층의 음속보다 낮으면서 두께 추정 결과의 전체적인 변화 양상을 파악하기 위해 약 1450 ~ 1505 m/s의 표층 퇴적층 음속범위를 설정하였으며, 표층 퇴적층 반사 신호와 하부 퇴적층 반사 신호의 도달 시간 차이를 고려하여 하부 퇴적층 두께 추정 범위를 0 ~ 10 m로 적용하였다. 또한 총 이동거리 200 m 이내의 짧은 전파거리이기 때문에 음파의 굴절 효과를 무시하였다.
Fig. 5는 실험 세트 1, 주파수 4 kHz 경우에서 모의된 두 다중경로의 도달시간 차이(
)를 각기 다른 색으로 구분하여 나타내는 기법(의사 색채 기법: Pseudo color plot)으로 표현한 그림이며, 이 때 실험적으로 측정한 도달시간 차이 
에 해당되는 시간을 하얀색으로 구분하여 나타냈다(: 
, 1.8 ms) 그 결과, 한국지질자원연구원에서 제공된 실험 해역의 추정 평균 입도 범위인 8 ~ 10 ϕ에 해당하는 퇴적층 음속(검은색 점선: 1480 ~ 1505 m/s)에서는 4.4 ~ 5.5 m의 하부퇴적층 두께를 갖는 것으로 추정되었다.
). The solid line is the estimated sound speed of surficial sediment and the dotted lines indicate its error range.IV. 결과 및 토의
실험 세트 1, 2, 3의 주파수 4, 6, 8 kHz에서 각각 표층 해저면 반사 신호와 해수면 반사 신호의 사이에 하부 퇴적층 반사 신호가 수신되었다(Fig. 6). 같은 실험 세트 내에서 하부 퇴적층 반사 신호는 주파수에 따라 수신 준위의 차이를 보였지만, 표층 해저면 반사 신호와의 도달 시간 차이는 매우 유사한 경향을 보였다. Fig. 7은 앞에서 계산한 방법으로 각 실험 세트별 
를 사용하여 퇴적층 두께를 추정한 전체 결과이다. 퇴적층 두께를 추정한 결과는 각 실험 세트 별로 다소 차이는 있지만 약 5 m를 중심으로 4 ~ 7 m의 범위에 분포하였다(Table 1).
|
Fig. 6. Received intensities for experimental set of 1, 2, and 3 with frequencies of 4, 6, and 8 kHz. |
| |||
Fig. 7. Geoacoustic inversion result of the surficial sediment thickness for all experimental set (1, 2, and 3) and frequency (4, 6, and 8 kHz) combinations. |
본 논문에서 음향신호를 이용하여 추정된 하부 퇴적층의 두께 결과는 한국 해양과학기술원에서 제공받은 실험해역의 천부지층탐사 자료와 비교하였다. 2개의 정선에 대한 천부지층탐사 자료는 각각 8 km의 수평거리에 대한 측정 결과이다(Fig. 8). 정선 1 (Line 1)은 북에서 남쪽방향, 정선 2(Line 2)는 서에서 동쪽으로 이동하며 측정되었으며, 해저면 반사 신호를 측정한 실험 정점과 정선 1은 최단거리로 약 200 m, 정선 2와는 약 1.1 km 떨어져있었다. Fig. 8에서 정선 1, 2 옆의 화살표는 음향 실험이 수행된 정점과 최단거리에 해당하는 천부지층탐사 해역이다. 이 해역은 표층 해저면 하부에 얇은 하부 퇴적층이 분포하는 것을 확인할 수 있으며, 하부 퇴적층 두께는 약 10 m 이하로 구성되어 있는 것을 알 수 있다(Fig. 9). 보유하고 있는 천부지층탐사 결과가 그림 자료의 형태이기 때문에 정밀하게 퇴적층 수심을 파악할 수 없으며, 또한 실험 해역이 정선 1과는 약 200 m, 정선 2와는 약 1.1 km 떨어져 있기 때문에 실험 해역을 정확히 대표할 수 없다는 제한점이 있다. 하지만 천부지층탐사 결과(약 10 m 이하의 표층 퇴적층 구조)를 토대로 보았을 때, 본 연구에서 추정한 하부 퇴적층의 두께 역산 결과(약 4~7 m)는 타당한 것으로 판단할 수 있다.
V. 요약 및 결론
본 논문에서는 해저면 반사 신호 측정 실험에서 나타난 표층 해저면 반사 신호와 하부 퇴적층 반사 신호의 도달 시간차이를 이용하여 하부 퇴적층의 두께를 역으로 추정하였다. 수평 입사각에 따른 해저면 반사 신호 측정 실험은 경남 거제시 장승포항 주변 해역에서 수행하였으며, 측정된 수신 신호에서 나타난 표층 해저면 반사 신호와 하부 퇴적층 반사 신호의 도달 시간 차이 
를 계산하였다. 이는 각각의 실험세트 별로 하부 퇴적층 두께와 퇴적층 음속의 함수로 모의한 표층 해저면과 하부 퇴적층 반사 신호의 도달시간 차이 
와 비교하였으며, 해당되는 반사 경로 모의 환경을 확인함으로써 하부 퇴적층의 두께를 추정할 수 있었다. 실험해역의 퇴적물 특성은 한국지질자원연구원에서 제공받은 추정 평균 입도 범위인 8~10 ϕ에 해당되는 퇴적층 음속(1480~1505 m/s)을 이용하였으며, 이를 이용하여 추정한 하부 퇴적층의 두께는 약 4~7 m로 도출되었다. 또한 추정한 하부퇴적층 두께 역산 결과는 한국해양과학기술원에서 제공받은 거제 인근 해역의 천부지층탐사 자료와 비교하였으며, 본 논문의 역산 결과가 타당한 것으로 판단할 수 있었다.
본 논문에서는 해저면 하부퇴적층으로부터 반사되어 수신된 신호의 도달시간을 이용하여 하부퇴적층 두께만을 역산하였지만, 추후에는 다층구조 반사손실 모델을 이용하여 하부퇴적층 두께 외에도 퇴적층 음속, 밀도, 감쇄계수와 같은 지음향 인자를 역산하기 위한 연구를 수행할 계획이다.
