KR101548667B1 - 지자기교란요소 제거를 통한 해양자력탐사 총지자장 측정방법, 시스템 및 이를 위한 기록매체 - Google Patents

Abstract

본 발명은 해양자력탐사 측정방법 및 시스템에 관한 것이다.
본 발명에 따른 해양자력탐사 측정방법은, 조사 대상 영역에서 중심점을 기준으로 원형의 항로를 그리며 탐사선을 진행시키면서 자력계를 예인하며, 진북을 기준으로 시계방향으로 측정된 상기 탐사선의 진행 각도에 따라 상기 자력계를 이용하여 지자장값을 측정하는 측정단계; 탐사선의 각도에 대한 데이터와, 각도에서 측정된 지자장값에 대한 데이터를 모두 연산유닛에서 입력받는 데이터 입력단계; 연산프로그램이 수록되어 있는 연산유닛에서 측정된 지자장값에 포함되어 있는 탐사선에 의해 유도된 교란자기장값을 설정하는 연산단계; 및 측정된 자기장값으로부터 교란자기장값을 제하여 실제 지자장값을 산출하는 단계;를 포함하여 이루어진다.

Description

지자기교란요소 제거를 통한 해양자력탐사 총지자장 측정방법, 시스템 및 이를 위한 기록매체{Method and System for getting marine magnetic data by elimanating magnetic field disturbance from ship's heading effect and Recording media thereof}

본 발명은 해저의 금속 부존자원, 대규모 지질구조 및 파이프라인 등을 탐사하기 위한 해양자력탐사에 관한 것으로서, 더욱 상세하게는 해양자력탐사시 금속으로 이루어진 탐사선에 의해 자력탐사 측정값이 왜곡되는 것을 보정함으로써 정확하게 총지자장값을 취득하기 위한 방법 및 시스템, 그리고 이 방법을 컴퓨터에서 수행할 수 있도록 컴퓨터에서 읽기가능한 기록매체에 관한 것이다.

도 1은 해양자력탐사를 수행하는 모습을 나타낸 모식도이다.

도 1을 참고하면, 일반적으로 해양자력탐사는 탐사선의 후미에 해상자력계(marine magnetmeter)를 예인하면서 해상자력계에서 총지구자기장을 측정하는 방법을 말한다.

지구자기장(earth marnetic field)의 구성은, 도 2의 지구 핵 내부의 액체금속 운동에 의하여 야기되는 내부 요인이 대략 98%를차지하고 나머지 2% 정도는 태양에 의한 외부요인이다. 지구중력장과 달리 지구자기장은 지속적으로 변화하며 이러한 변화는 하루를 주기로 변화하는 일변화(diurnal variation), 수십 내지 수백년 주기로 천천히 변화하는 영년변화(secular variation) 그리고 급격히 변화하는 자기폭풍(magnetic storm)이 있다.

도 3에 도시된 바와 같이, 지구자기장은 크기와 방향을 가지는 벡터로 구성되어 있으며 지구자기장의 크기인 총자기 F, 지자기장의 수평성분 H와 진분 사이의 각인 편각(diclination) D, 그리고 H와 F 사이의 복각(inclination) I로 나타내며 이를 지자기의 3요소라 한다. 지구자기장의 단위는 nT(nano-Tesla)로 표시하며 적도와 극에서 각각 30,000nT 및 65,000nT 정도의 값을 나타낸다.

자력탐사는 위의 표 1에 나타낸 바와 같이 지표부근의 파이프라인이나 케이블을 찾는 소규모 조사부터 석유가스 자원탐사와 같은 대규모 탐사까지 적용 범위가 넓다는 장점이 있다. 자력탐사법이 많이 이용되면서 암석 및 광물의 자성에 관한 연구도 활발해 졌으며 특히 암석이 띄고 있는 잔류자기에 대한 연구와 응용은 1950년대 이후 지구물리학에서 중요한 위치를 차지하는 고지자기학으로 발전하여 대륙이동성이나 해양저 확장설과 같은 현대 지질학의 주요 학설을 입증하는 정량적인 증거로 제시되었다. 그리고 산업의 측면에서도 특히 석유가스자원 탐사시 기반암 등의 지하구조 규명에 이용되어 중력 자료 및 탄성파 자료와 함께 상호 보완자료로 유용하게 활용되고 있다.

즉, 해양자력탐사를 통해 조사대상영역의 총지자장을 측정하면서, 지자장의 변화가 급격하게 일어나는 것을 파악하여 금속 소재의 케이블이나 파이프라인 또는 금속 광맥 등 자기장에 영향을 미치는 요소를 탐사할 수 있다. 이에 해양자력탐사에서는 총지자장의 값을 정확하게 측정하는 것이 무엇보다도 중요하다. 그러나 해상자력계에서 측정된 값에는 탐사선(강자성을 띠는 철 소재로 제조)에 의해 유도된 자기장이 포함되어 있어 총지자장의 순수한 값을 취득하기 어렵다는 문제가 있었다.

본 발명을 상기한 문제점을 해결하기 위한 것으로서, 해양 자력탐사시 금속 소재로 이루어진 탐사선에 의해 유도된 자기장의 영향을 제거하여 탐사 대상 영역의 총지장을 정확하게 측정할 수 있는 지자기교란요소 제거를 통한 해양자력탐사 총지자장 측정방법, 시스템 및 컴퓨터 기록매체를 제공하는데 그 목적이 있다.

상기 목적을 달성하기 위한 본 발명에 따른 해양자력탐사 총지자장 측정방법은, 조사 대상 영역에서 중심점을 기준으로 원형의 항로를 그리며 탐사선을 진행시키면서 자력계를 예인하며, 진북을 기준으로 시계방향(또는 반시계방향)으로 측정된 탐사선의 진행 각도(진북방향일 때 0°)에 따라 상기 자력계를 이용하여 지자장값을 측정하는 측정단계; 상기 탐사선의 진행 각도에 대한 데이터와, 상기 각도에서 측정된 지자장값에 대한 데이터를 모두 연산유닛에서 입력받는 데이터 입력단계; 연산프로그램이 수록되어 있는 상기 연산유닛에서 상기 측정된 지자장값에 포함되어 있는 상기 탐사선에 의해 유도된 °°교란자기장값을 설정하는 연산단계; 및 상기 측정된 자기장값으로부터 상기 교란자기장값을 제하여 실제 지자장값을 산출하는 단계;를 포함하여 이루어진 것에 특징이 있다.

본 발명에서는 상기 연산단계에서 상기 교란자기장값은 하기의 식1에 의해서 산정되며,

교란자기장값 = C₀ + C₁cosθ + C₂cos2θ ... 식1

(여기서, θ는 탐사선의 방향(heading)으로 진북을 기준으로 (반)시계방향으로 측정한 탐사선의 진행 각도이며, C₀, C₁, C₂는 탐사선의 자성특성, 탐사지역 및 자력계와 탐사선의 거리에 따른 상수임)

상기 연산프로그램에서는 상기 식1의 상수를 결정할 수 있다.

또한, 본 발명의 일 예에서, 상기 연산프로그램에서는 최초시점(상기 탐사선과 진행 각도 0°)에서 측정한 지자장값을 기준으로 하여 상기 탐사선이 진행 각도를 변경하며 진행하면서 연속적으로 측정하는 지자장값들의 변화량과, 상기 탐사선이 진행 각도를 변경하며 진행하면서 유도되는 각 교란자기장값들과의 차이가 최소화되도록, 상기 교란자기장값에 대한 식1의 상기 상수 C₀, C₁, C₂를 결정한다.

본 발명에 따른 해양자력탐사 총지자장 측정방법 및 시스템에 의하여, 탐사 대상지역에서 선박을 이용한 자력탐사를 수행할 때, 선박에 의하여 유도된 자기장의 영향을 제거하여 탐사 대상지역의 총지자장 값을 정확하게 측정할 수 있다는 이점이 있다.

또한 본 발명에서는 이미 측정되어 있는 데이터들을 이용하여 탐사가 완료된 후에도 자력탐사 데이터에 대한 후처리 보정이 가능하다는 이점이 있다.

또한 본 발명에 따른 시스템 및 방법은 자력탐사에 사용되는 자력계의 종류와 상관없이 모두 적용가능하다는 이점이 있다.

도 1은 해양자력탐사를 설명하기 위한 개략적 모식도이다.
도 2는 지구자기장의 구성을 설명하기 위한 도면이다.
도 3은 지구자기장의 3요소를 설명하기 위한 도면이다.
도 4는 탐사선의 방향(heading)에 따라 측정된 총지자장값과, 이 지자장값과 가장 근사한 곡선(best fit curve)을 나타낸 그래프이다.
도 5는 본 발명에 따른 총지자장 측정방법의 개략적 흐름도이다.
도 6은 원형 탐사를 설명하기 위한 도면이다.
도 7은 최소자승법의 개념을 나타낸 그래프이다.
도 8은 Joides Resolution 탐사선박이 해양탐사를 하는 모습을 나타낸 모식도이며, 도 9는 Joides Resolution 탐사선박에 의해 원형탐사가 수행된 "Ori Massif" 지역의 사이트맵이다.
도 10은 실험에서 자력계에 의하여 측정된 자기장값과, 이 자기장값들 및 본 연산 프로그램을 이용하여 형성된 근사 곡선(best fit curve)을 함께 나타낸 그래프이다.

앞에서 설명한 바와 같이, 해양자력탐사에서 탐사선에 의해서 유도된 자기장의 영향으로 인하여 조사대상영역의 총지자장에 대한 정확한 값을 파악하기가 어려운 문제점이 있었다. 이에 'Bullard and Mason'은 1961년 자력측정자료에서 탐사선에 의한 영향을 제거하는 방법(하기의 식1)을 아래와 같이 제안하였으며, 본 발명은 'Bullard and Mason'의 방법으로부터 시작되었다.

F_Q = F + C₀ + C₁cosθ + C₂cos2θ + S₁sinθ + S₂sin2θ ... 식1

여기서, F_Q는 해상자력계에서 측정된 값, F는 총지자장, θ는 탐사선의 방향(heading)으로 북에서부터 시계방향으로 측정된 방향, C₀, C₁, C₂, S₁, S₂는 탐사선의 자성특성, 탐사지역 및 자력계와 탐사선의 거리에 따른 상수이다. 여기서, 좌우가 동형으로 대칭인 탐사선의 경우에는 sine값의 영향이 cosine값에 비해 매우 미미하여 S₁과 S₂는 0으로 설정할 수 있다(Bullard and Mason, 1961). 이에 따라 식1을 다시 정리하면 아래의 식2와 같이 된다.

F_Q = F + C ₀ + C ₁ cosθ + C ₂ cos2θ ... 식2

위 식2에서 밑줄친 부분은 탐사선의 영향으로 유도된 자기장에 의한 것으로서 측정된 지자장(F_Q) 값을 왜곡시키는 요소이다. 따라서 C₀, C₁, C₂ 값만 알면 측정된 측정된 자기장(F_Q)에서 밑줄친 요소를 제거하여 탐사지역의 순수 총지자장(F)을 구할 수 있다.

'Bullard and Mason'은 탐사선의 영향을 제거하는 방법을 제시하였다.

도 4는 탐사선의 방향(heading)에 따라 측정된 총지자장값과, 이 지자장값과 가장 근사한 곡선(best fit curve)을 나타낸 그래프이다.

도 4를 참고하면, 탐사선의 진행방향(수평축)은 0~360°로 나타나 있다. 즉 탐사선은 특정 지점으로부터 원형을 그리며 다시 제자리로 돌아오는 경로로 진행하였다. 그리고 각 경로의 지점(탐사선의 헤딩 각도)에서 지자장값을 측정하여 그래프의 수직축에 점으로 찍어 나타내었다. 그리고 'Bullard and Mason'은 측정된 값들과 가장 근사한 커브(best fit curve)를 구함으로써 위의 상수값들(C₀, C₁, C₂)을 파악할 수 있다고 제시하였다. 그러나 아쉽게도 'Bullard and Mason'은 측정값과 가장 유사한 커브를 구하는 방법에 대해서는 구체적으로 제시하지 않음으로써 위 방법을 사용하는데는 한계가 있었다.

본 발명에서는 측정된 값들에 대하여 행렬(matrix)과 최소자승법(least squares method) 기법을 도입하여 상수값(C₀, C₁, C₂)을 구하며, 이를 통해 해양탐사에서 총지자장을 정확하게 측정할 수 있는 방법을 제시하고자 한다.

이하, 본 발명의 일 실시예에 따른 지자기교란요소 제거를 통한 해양자력탐사 총지자장 측정방법(이하, '총지자장 측정방법', '총지자장 측정시스템')에 대하여 더욱 상세히 설명한다. 본 발명에 따른 총지자장 측정방법은 시계열적인 방법에 관한 것이며, 컴퓨터로 읽을 수 있는 프로그램에 의하여 보다 구체적으로 구현된다. 따라서, 본 발명은 프로그램이 설치되어 있는 컴퓨터 또는 이 프로그램이 수록되어 있는 기록매체를 이용하여 컴퓨터에서 프로그램을 가동하여 수행한다.

먼저 본 발명에 따른 총지자장 측정시스템은 탐사선, 해상자력계 및 컴퓨터 프로그램이 수록되어 있는 연산유닛으로 이루어진다. 앞에서 설명한 바와 같이, 탐사선은 금속 소재로 이루어지며, 해상자력계는 와이어라인에 의해 탐사선에 연결된 상태로 바다 내에 배치된다. 탐사선이 이동함에 따라 해상자력계는 예인되면서 각 지점의 지자장 값을 측정한다. 그리고 해상자력계에 의하여 획득된 지자장 데이터는 전기적 연결수단에 의하여 연산유닛으로 송신되며, 연산유닛에서는 뒤에서 설명할 연산 알고리즘에 의하여 탐사선에 의한 지자장 왜곡 요소를 제거하여 탐사지역의 순수한 총지자장 값을 산출한다.

상기한 시스템을 이용하여 본 발명에 따른 총지자장 측정방법에 대하여 설명한다.

도 5는 본 발명에 따른 총지자장 측정방법의 개략적 흐름도이며, 도 6은 원형 탐사를 설명하기 위한 도면이다.

도 5 및 도 6을 참고하면, 본 발명에 따른 총지자장 측정방법은 먼저 조사 영역에서 탐사선이 원형을 그리도록 진행시키면서 각 지점(또는 탐사선의 헤딩 각도)에서 지자장값을 측정하는 것으로부터 시작한다. 탐사선의 시작 방향은 상관없으며, 원형을 그리며 탐사선을 진행시키면서 자력값을 측정한다. 측정된 데이터들은 와이어라인에 포함되어 있거나 또는 와이어라인과 별도의 유무선통신라인을 통해 연산유닛으로 전송된다.

연산유닛에서는 측정된 데이터에 대하여 본 발명에 의한 연산 알고리즘을 적용하여 탐사선에 의하여 유도된 자기장의 영향을 결정하는 상수값(C₀, C₁, C₂)을 산출해 낸다.

그리고 이 상수값을 상기의 식2에 입력하여 측정된 총지자장값(F_Q)으로부터 탐사선에 의한 영향을 제거하여 탐사 영역의 순수한 총지자장값(F)를 획득할 수 있다.

이하에서는 본 발명의 연구진에 의하여 개발된 연산 프로그램에서 핵심이 되는 수식은 행렬 A와 행렬 C의 곱에서 행렬 Y를 제한 후, 이 값에 대하여 최소자승법(min과 자승으로 표현된 부분)와 2-Norm(2개의 절대값에 하첨자 2로 표현된 부분)을 도입하여 아래의 식3과 같이 표현된다.

min||AC-Y||² ₂ ... 식3

그리고 위 식3을 풀어서 우리가 구하고자 하는 자기교란상수값(C₀, C₁, C₂)에 대한 행렬 C를 아래의 식4에 의하여 구할 수 있다(식3을 푸는 과정은 뒤에서 후술함).

C = (A^TA)^-1·(A^TY) ... 식4

위 식4에서 상첨자 T는 각 행렬의 전치행렬을 의미한다. 예컨대, A^T는 행렬 A의 전치행렬이다.

먼저 연산 프로그램에 도입한 최소자승법과 2-Norm 개념은 대수학에서 주지의 개념이므로 간략하게만 소개한다. 도 7의 그래프를 참고하면 측정값들(빨간색 점)과, 이 측정값들의 경향에 가장 근사한 커브(파란색)를 함수 y로 나타내었다. 그리고 실제 측정값과 함수값의 차이(error)를 표시하였다. 최소자승법을 이용하여 측정된 값들과 가장 근사한 커브를 구할 수 있다. 이에 더하여 본 연구진에서는 최소자승법에 Norm이라는 수학적 도구를 접목한다. Norm이란 방정식의 수치적인 해의 오차를 효율적으로 분석하기 위한 것이다.

즉, X=(x₁, x₂, x₃ .... x_n)일 때 1-Norm과 2-Norm은 아래와 같이 정의된다.

1-Norm은 각각의 값에 절대값을 취하여 합한 것이고, 2-Norm은 각 값의 제곱하고 이들을 모두 합한 후 1/2승한 것으로서 본 연산 프로그램에서는 2-Norm을 이용한다.

본 발명에서 채용하는 연산 프로그램의 핵심이 되는 위의 식3은 측정된 자기장값과 실제 자기장값 사이의 차이, 즉 탐사선에 의하여 유도된 자기장값(이하, 교란자기장값)을 정의하면서 도출된다. 즉,위 식2를 이항정리하여, 실측된 자력값(F_Q)에서 탐사지역의 순수 자력값(F)을 제하여 아래와 같이 나타낼 수 있다.

F_Q - F = C ₀ + C ₁ cosθ + C ₂ cos2θ

그리고 탐사선의 진행방향(θ)에 따라 값들이 각각 생성되므로, 이들을 함께 표현하기 위하여 행렬을 도입하면 아래와 같다.

즉 A행렬과 C행렬의 곱(A·C)은 원형 탐사의 각 지점에서의 교란자기장값을 함께 표현한 것이다.

그리고 Y는 탐사선의 진행방향 0°(시작지점)에서 측정된 총자기장 값을 기준으로 설정하였을 때 탐사선 방향에 따라 나타나는 총자기장의 변화량으로 나타낼 수 있다. 즉, 탐사선 방향이 0°일 때에는 탐사선이 진행하고 있지 않은 상태이므로 자력계에서 실측된 값(F_Q)은 탐사선에 의한 교란자기장값의 영향을 받지 않은 순수한 지자장값(F₀)으로 볼 수 있다. 그리고 탐사선이 방향을 변화시키며 이동함에 따라 실측된 지자기의 변화값(F_Q-F₀)을 탐사선에 의해 발생한 교란자기장값(y_n)으로 볼 수 있다.

y_n = F_Q - F₀ = C₀ + C₁cosθ + C₂cos2θ

따라서 A행렬과 C행렬의 곱(AC)에서 Y행렬을 제하는 행렬 함수(AC-Y)를 설정한다. 그리고 이 값에 대하여 최소자승법과 2-Norm을 도입하여 위와 같이 식3 [min||AC-Y||² ₂ ]을 설정하고, 이 식3의 값이 최소값을 보이는 행렬C를 구하는 것이다. 위 식3이 최소가되는 행렬C에 대한 해를 구한다는 것은, 실측된 값들과 가장 근사한 커브를 구할 수 있다는 것을 의미한다.

이제 위의 식3을 풀어서 행렬C에 대한 해를 구하는 과정을 설명한다.

위의 식3의 행렬에 대하여 2-Norm과 제곱을 풀면 아래와 같이 나타낼 수 있다.

min||AC-Y||² ₂ = min(AC-Y)^T(AC-Y)

= min[(AC)^T(AC)-(AC)^TY-Y^T(AC)+Y^TY]

= min(C^TA^TAC-2C^TA^TY+Y^TY) ... 중간식1

위에서 상첨자 T는 행렬의 전치행렬(transpose matrix)을 나타낸다.

따라서 중간식1에서 C^TA^TAC-2C^TA^TY+Y^TY 값은 C²A²-2CAY+Y² 으로 바꿀 수 있다. 이는 전치행렬은 단지 열과 행을 바꾼 것일 뿐 절대값은 동일하기 때문이다.

행렬 C를 구하기 위하여 위의 중간식1에 대하여 편미분을 통해 아래의 중간식2를 도출한다.

∇(C^TA^TAC-2C^TA^TY+Y^TY)C = 2A^TAC-2A^TY = 0 ... 중간식2

즉, C²A²-2CAY+Y²을 C에 대하여 편미분을 하면 2CA²-2AY 가 된다. 그러나 행렬의 연산에는 순서를 바꿀 수 없기 때문에 위와 같이 순서대로 배열한 것이다.

따라서 위 중간식2의 2A^TAC-2A^TY=0에서 C에 대하여 풀면 위의 식4를 유도할 수 있다. 식 4를 다시 한 번 기재한다.

C = (A^TA)^-1·(A^TY) ... 식4

여기서, (A^TA)^-1는 A^TA의 역행렬이며, Y는 탐사선 방향 0°에서 측정된 총자기장 값을 기준으로 설정하였을 때 탐사선 방향에 따라 측정된 총자기장의 변화량(delta)으로 나타낼 수 있다. 따라서 위의 식4를 풀면 행렬C의 해(C₀,C₁,C₂)를 구할 수 있다. 본 발명의 일 구성요소인 연산유닛에서는 위의 식4가 기록되어 있어서, 데이터가 입력되면 행렬C의 해를 구하는 연산을 하게 된다.

행렬의 해를 구하면 탐사선에 의해 유도된 교란자기장값(C₀ + C₁cosθ + C₂cos2θ)을 알 수 있으므로, 탐사선의 각 방향에서 실측된 총자기장값으로부터 교란자기장값을 제하여 탐사지역의 순수한 총자기장값을 산출할 수 있게 된다.

즉, 본 발명에서는 원형탐사의 최초 시작시점(탐사선 헤딩각도 0°)에서 측정한 총자기장값이 실제 자기장값으로 선정함으로써 탐사선에 의해 유도된 교란자기장의 영향이 없는 것으로 설정(AC-Y=0)한다. 이후 탐사선이 진행하면서 측정되는 총자기장값(Y, 시작시점에 AC-Y=0 으로 설정했으므로 추후의 Y값은 실제로는 시작시점으로부터 변화량을 의미)의 차이를 최소화시키는 행렬C에 대한 해를 구한다. 행렬C의 해를 구하면 탐사선에 의해 유도된 교란자기장값을 산출할 수 있으므로, 실제 측정된 총자기장값으로부터 교란자기장값을 제하여 탐사영역의 실제 총자기장값을 정확하게 산출할 수 있다.

이하에서는 식4의 행렬C를 수학적으로 구하는 과정의 예에 대하여 설명한다. 수학적 복잡함을 피하기 위하여 4개의 각도(θ₁, θ₂, θ₃, θ₄)에서 측정한 값으로 간소화시켰다. 아래와 같이 STEP 1~4의 4개의 단계로 연산이 이루어진다.

상기한 바와 같이, 연산유닛에서는 탐사선의 진행각도(θ)에 대한 데이터와, 각 각도에서 측정된 자기장값에 대한 데이터를 입력받아 상기한 바와 같은 행렬을 풀어내어 자기교란상수 C(C₀, C₁, C₂)를 구하여 탐사선에 의해 유도된 교란자기장값을 획득할 수 있다. 그리고 이 교란자기장값을 실측된 자기장값으로부터 제하여 탐사영역의 실제 자기장값을 정확하게 파악할 수 있다.

본 연구진은 본 발명에 따른 총자기장 측정방법과 연산 프로그램의 효과를 알아보기 위하여 실험을 진행하였다. 도 8은 "Joides Resolution" 탐사선박이 해양탐사를 하는 모습을 나타낸 모식도이며, 도 9는 "Joides Resolution" 탐사선박에 의해 원형탐사가 수행된 "Ori Massif" 지역의 사이트맵이다.

도 8 및 도 9에 도시된 바와 같이, 탐사선박인 "Joides Resolution"을 이용하여 "Ori Massif" 지역(도 9 참고)에서 원형 탐사를 수행하였다. 실험에서 자력계에 의하여 측정된 자기장값과 본 연산 프로그램에 의하여 형성된 근사 곡선(best fit curve)을 도 10에 함께 나타내었다.

위 실험을 통해 "Ori Massif" 지역에서 "Joides Resolution" 선박을 이용하여 자력탐사를 수행하고, 연산 프로그램을 이용하여 아래와 같이 자기교란상수C(C₀, C₁, C₂)를 구하였고 탐사선박에 의한 교란자기장값(F_H)을 구하였다.

F_H = C₀ + C₁cosθ + C₂cos2θ = 40.91 - 45.89cosθ + 3.66cos2θ

도 10을 참고하면, 연산 프로그램에 의하여 상수값 C를 결정하여 best fit curve를 도출한 결과 실제 측정된 값(그래프 상의 점들)과 best fit curve가 거의 일치하는 것으로 나타났다.

아래의 표2 및 표3에는 위의 실험에서 본 연산 프로그램에 의하여 탐사선에 의한 자기교란을 보정하기 전과 보정한 후의 교차점(cross-over) 오차 비교를 수행를 수행한 결과가 나타나 있다.

[표 2] - 자기교란 보정전 cross-over 오차

자기교란 보정전에는 평균오차 20.2nT, 최대오차 48.2nT, 최소오차 5.7nT 로 나타났다.

[표 3] - 자기교란 보정후 cross-over 오차

자기교란 보정후에는 평균오차 6.2nT, 최대오차 12.1nT, 최소오차 -0.1nT 로 나타났다. 즉, 본 연산프로그램에 의하여 보정을 수행한 결과 오차 범위가 현저하게 줄어들었음을 알 수 있다. 이는 본 발명에 의하여 원형탐사를 수행하고 연산 프로그램을 작동시키면 해상 자력탐사에서 탐사선에 의해 유도된 교란자기장의 영향을 최소화시켜 신뢰성 있는 자력탐사를 수행할 수 있다는 것을 의미한다.

특히, 본 발명은 특정한 자력계에만 적용되는 것이 아니고 자력계의 종류에 관계없이 활용가능하다. 더욱이 탐사가 이미 수행된 영역이라고 하더라도 탐사선의 헤딩각도에 따른 자기장 측정값만 있다면 본 연산프로그램을 수행하여 데이터를 후처리를 통해 보정이 가능하다는 이점이 있다.

Claims (7)

1. 조사 대상 영역에서 중심점을 기준으로 원형의 항로를 그리며 탐사선을 진행시키면서 자력계를 예인하면서, 진북을 기준으로 시계방향으로 측정된 상기 탐사선의 진행 각도에 따라 상기 자력계를 이용하여 지자장값을 측정하는 측정단계;
   상기 탐사선의 각도에 대한 데이터와, 상기 각도에서 측정된 지자장값에 대한 데이터를 모두 연산유닛에서 입력받는 데이터 입력단계;
   연산프로그램이 수록되어 있는 상기 연산유닛에서 상기 측정된 지자장값에 포함되어 있는 상기 탐사선에 의해 유도된 교란자기장값을 설정하는 연산단계; 및
   상기 측정된 자기장값으로부터 상기 교란자기장값을 제하여 실제 지자장값을 산출하는 단계;를 포함하여 이루어지며,
   상기 연산단계에서 상기 교란자기장값은 하기의 식1에 의해서 산정되며,
   교란자기장값 = C₀ + C₁cosθ + C₂cos2θ ... 식1
   (여기서, θ는 탐사선의 방향(heading)으로 진북을 기준으로 시계방향으로 측정한 탐사선의 진행 각도이며, C₀, C₁, C₂는 탐사선의 자성특성, 탐사지역 및 자력계와 탐사선의 거리에 따른 상수임)
   상기 연산프로그램에서는 상기 식1의 상수를 결정하는 것을 특징으로 하는 해양자력탐사 총지자장 측정방법.
2. 삭제
3. 제1항에 있어서,
   상기 연산프로그램에서는 최초시점(상기 탐사선의 진행각도 0°진북방향)에서 측정한 지자장값을 기준으로 하여 상기 탐사선이 각도를 변경하며 진행하면서 연속적으로 측정하는 지자장값들의 변화량과, 상기 탐사선이 각도를 변경하며 진행하면서 유도되는 각 교란자기장값들과의 차이가 최소화되도록,
   상기 교란자기장값에 대한 식1의 상기 상수 C₀, C₁, C₂를 결정하는 것을 특징으로 하는 해양자력탐사 총지자장 측정방법.
4. 제1항에 있어서,
   상기 연산프로그램에서는 아래의 식2의 행렬을 연산하여 상기 상수 C₀, C₁, C₂를 결정하는 것을 특징으로 하는 해양자력탐사 총지자장 측정방법
   C = (A^TA)^-1·(A^TY) ... 식2
   (여기서, 아래에 표시한 바와 같이 A행렬은 상기 탐사선의 헤딩 각도를 변수로 하는 행렬이며, C행렬은 위 식1의 상수값들에 대한 행렬이고, Y는 상기 탐사선이 진행하며 실측한 지자장값 또는 최초시점(상기 탐사선의 진행 각도 0°)에서 측정한 지자장값을 기준으로 하여 상기 탐사선이 각도를 변경하며 진행하면서 연속적으로 측정하는 지자장값들의 변화량이며, 상첨자 T는 전치행렬이고, 상첨자 -1은 역행렬을 의미함)
   

   
5. 제4항에 있어서,
   상기 식2는 최소자승법과 2-Norm을 도입한 하기의 식3으로부터 유도된 것을 특징으로 하는 해양자력탐사 총지자장 측정방법.
   min||AC-Y||² ₂ ... 식3
6. 수중에서 운행되는 탐사선;
   와이어라인에 의하여 상기 탐사선에 연결되어 상기 탐사선에 의해 예인되면서 조사대상영역의 지자장을 측정하는 자력계; 및
   유무선 라인을 통해 상기 자력계로부터 측정된 지자장에 대한 데이터 및 상기 탐사선의 진행각도를 입력받으며, 상기 탐사선에 의하여 유도된 교란자기장값을 산출하는 컴퓨터로 읽기가능한 연산프로그램이 수록되어 있는 연산유닛;을 구비하며,
   상기 연산프로그램에서는 청구항 1, 청구항 3 및 청구항 4 중 어느 하나에 기재되어 있는 연산단계를 수행하는 것을 특징으로 하는 해양자력탐사 총지자장 측정 시스템.
7. 청구항 1, 청구항 3 및 청구항 4 중 어느 하나에 기재되어 있는 연산프로그램이 수록되어 있는 컴퓨터로 판독가능한 기록매체.

Images