Matlab 프로그램을 이용한 해저 탄성파 탐사 단면도 층후 산정 및 층후도 생성

Abstract

해저 탄성파 탐사 단면도에서 퇴적층의 층후 계산을 보다 효율적으로 하기 위하여 MATLAB을 활용한 스크립트(script) 프로그램을 작성하였다. 특히 효율성을 높이기 위하여 SEG-Y 등의 별도 데이터 형식이 아닌 해저 탄성파 탐사 단면 이미지만으로도 층후를 계산하여 층후도를 생성할 수 있도록 코딩하였다. 우선 탄성파 탐사 처리 소프트웨어(Seisee 등)를 이용하여 X축(Shot Number)과 Y축(TWTT)의 스케일을 고정한, 단면 이미지를 파일로 전환한다. 그리고 MATLAB 프로그램에서 단면이미지를 불러들여 1개의 픽셀에 해당하는 Shot Number와 TWTT(TWO-WAY Travel Time)의 스케일을 산출한다. 그리고 해저면(seafloor)과 음향기저면(acoustic basement)의 층(layer)을 RGB값으로 분류하여 층후를 계산하고, 층후도를 생성하는 알고리즘을 적용하였다. MATLAB 스크립트 프로그램을 실행하기 위해서는 처리할 탐사라인에 대한 정보 파일, 탄성파 탐사 이미지 파일 그리고 항해 파일(navigation file)의 입력 조건이 필요하다. 정보 파일의 경우, 이미지 이름, Shot Number Tick Mark 간격(X축), Shot Number Max data, TWTT Tick Mark 간격(Y축), TWTT Max data, Navigation File name, 해저면 및 음향기저면의 color RGB data, Output Shot Number 간격(default:1)대한 내용을 포함하고 있어야 한다. 탄성파 단면도 이미지 파일(*.bmp)은 X축과 Y축에 스케일을 나타내주는 tick mark가 있어야 하고, 해저면과 음향기저면 층(Layer)의 색깔선(color line)으로 표시되어 있어야 한다. 그리고 항해 파일은 단면도 내의 각 shot number에 해당하는 위치정보(Shot Number, X, Y)가 있어야 한다. 프로그램 처리 결과후를 계산하여 층후도를 생성할 수 있도록 코딩하였다. 우선 탄성파 탐사 처리 소프트웨어(Seisee 등)를 이용하여 X축(Shot Number)과 Y축(TWTT)의 스케일을 고정한, 단면 이미지를 파일로 전환한다. 그리고 MATLAB 프로그램에서 단면이미지를 불러들여 1개의 픽셀에 해당하는 Shot Number와 TWTT(TWO-WAY Travel Time)의 스케일을 산출한다. 그리고 해저면(seafloor)과 음향기저면(acoustic basement)의 층(layer)을 RGB값으로 분류하여 층후를 계산하고, 층후도를 생성하는 알고리즘을 적용하였다. MATLAB 스크립트 프로그램을 실행하기 위해서는 처리할 탐사라인에 대한 정보 파일, 탄성파 탐사 이미지 파일 그리고 항&