KR20010090027A - 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 - Google Patents
탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 Download PDFInfo
- Publication number
- KR20010090027A KR20010090027A KR1020010045228A KR20010045228A KR20010090027A KR 20010090027 A KR20010090027 A KR 20010090027A KR 1020010045228 A KR1020010045228 A KR 1020010045228A KR 20010045228 A KR20010045228 A KR 20010045228A KR 20010090027 A KR20010090027 A KR 20010090027A
- Authority
- KR
- South Korea
- Prior art keywords
- time
- window
- data
- acoustic wave
- wave
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 99
- 238000012545 processing Methods 0.000 title claims abstract description 57
- 238000003672 processing method Methods 0.000 claims description 23
- 239000000463 material Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 8
- 238000001028 reflection method Methods 0.000 claims description 6
- 230000000379 polymerizing effect Effects 0.000 claims description 5
- 238000005422 blasting Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 abstract description 24
- 230000008569 process Effects 0.000 abstract description 20
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000006116 polymerization reaction Methods 0.000 description 28
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 125000003473 lipid group Chemical group 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 239000013049 sediment Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000283726 Bison Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005486 microgravity Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/34—Transmitting data to recording or processing apparatus; Recording data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
본 발명은 지하 지질 구조를 영상화하는 반사법을 이용한 탄성파 자료 처리를 단순화시킴으로써 전문적인 지식 없이도 지하 영상을 신속하게 얻을 수 있으며, 특히 속도에 관한 정보가 전혀 없는 상황에서도 실제와 근사한 결과를 도출하는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치에 관한 것이다. 본 발명에 따른 방법, 시스템 및 장치는 자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력받고, 상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하며, 상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 것을 특징으로 한다.
Description
본 발명은 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치에 관한 것으로서, 보다 상세하게는 지하 지질 구조를 영상화하는 반사법을 이용한 탄성파 자료 처리에 관한 것으로, 복잡하고 전문적인 여러 단계의 처리과정을 간단하게 단순화시킴으로써 전문적인 지식이 없이도 지하영상에 관한 정보를 얻을 수 있으며, 특히 속도에 관한 정보가 전혀 없는 상황에서도 실제와 근사한 결과를 도출하는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치에 관한 것이다.
지하 지질 구조를 영상화하기 위해서는 해당 지역의 지하로 탄성파를 보내어 그 탄성파에 대한 기초 정보 자료를 측정한 후 그 정보를 바탕으로 지질 구조를 분석 및 파악하여야 한다. 여기서 기초 정보란 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도 등을 의미한다.
도 1은 지질 구조를 영상화하기 위한 장치의 개략적인 구성도이다.
도 1을 참조하여 설명하면, 지질 구조를 영상화하기 위해서는 기초 정보를 측정하는 자료 획득 장치(100)와 측정된 기초 정보 자료를 바탕으로 지질 구조를 분석 및 파악하여 영상화하는 자료 처리 장치(150)가 필요하다.
자료 획득 장치(100)를 사용하여 기초 정보 자료를 측정하는 일 실시예를 나타낸 도 2를 참조하여, 자료 획득 장치에 대해 설명한다.
자료 획득 장치(100)는 신호를 발생시키는 발생부(110), 수신부(120), 기록부(130)로 구성되어 있다.
발생부(110)는 파원을 발생시키는 역할을 수행한다. 예를 들어, 육상의 경우에 심부 탐사에는 폭약을 발생부(110)로 사용하며, 천부 탐사에는 해머, 낙하추와 같은 중력추와 진동기 등을 발생부(110)로 사용한다.
수신부(120)는 상기 발생부(110)로부터 얻은 탄성파에 관한 자료를 수신하는 역할을 하며, 여러 개의 수신기를 포함하고 있다. 수신기는 육상의 경우 대표적으로 자석과 코일로 구성된 지오폰(geophone)과 해상의 경우 대표적으로 압전체(piezoelectric material)로 구성된 하이드로폰(hydrophone)이 사용되고 있다. 이하에서는 수신기 중 지오폰을 예로 들어 설명하기로 한다.
기록부(130)는 상기 수신부(120)에서 수신한 자료를 기록하며 그 자료를 자료 처리 장치(150)의 프로세서(160)에 전송하는 역할을 한다. 자료를 전송하는 방법으로는 케이블을 통해 획득한 자료를 자료 처리 장치에 직접 송신하거나, 획득한 자료를 저장 수단, 예를 들어 디스옵(diskette), 씨디(CD), 마그넷 테이프(magnet tape)에 저장하여 자료 처리 장치에 전송하는 방법 등이 이용되고 있다. 상기에서 언급한 바와 같이 자료 획득 장치(100)에서 획득된 자료는 자료 처리 장치(150)의 프로세서(160)가 지하의 영상 정보를 얻기 위해 처리하게 된다. 그리고 처리된 결과는 자료 처리 장치(150)의 디스플레이부(170)에 표시되게 된다.
이 때, 프로세서가 탄성파 자료의 기초 정보를 처리하는 방법이 문제가 된다.
종래의 프로세서(160)가 기초 정보 자료를 처리하기 위한 통상적인 기초 정보 처리 방법으로는 석유 탐사 등에 많이 쓰이고 있는 공중점 중합 방법이 주로 사용되어 왔다. 공중점 중합 방법은 자료 처리 장치(150)로 통상 워크 스태이션급의 이상의 전문 컴퓨터를 사용하여야 하며 10여단계의 처리 과정을 거쳐야 한다. 여기서 공중점이란 여러 공발점 즉, 발파점의 모음 자료를 재분류하여 하나의 반사점에서 반사된 여러 트레이스들은 모은 자료이다.
대표적인 처리 단계는 디멀티플렉스(demultiplex), 포맷 변환, 파원/수신기 위치 정보 입력, 실진폭 회수(true amplitude recovery), 고도 보정(elevation correction), 트레이스 편집, 주파수-파수 필터(f-k filter), 뮤트(내부 및 외부), 역콘볼루션(deconvolution), 광대역 주파수 필터, 공중점 분류(common midpoint sort), 속도 분석(velocity analysis), 수직 경로 시차 보정(normal moveout correction), 잔여정보정(residual statics correction), 중합(stack), 구조 보정(migration) 등이다.
대표적인 단계 중 몇 가지 중요한 단계의 역할에 대해 간단히 언급하면 다음과 같다.포맷 변환 단계는 기록 포맷을 자료 처리 소프트웨어 내부 포맷으로 변환하는 단계이다.
또한, 파원/수신기 위치 정보 입력 단계는 각 기록별 파원과 지오폰들의 위치에 관한 정보를 컴퓨터에 입력하는 단계이다.
또한, 고도 보정 단계는 발파점 및 수신점의 고도 변화에 의한 시간차를 보정하는 단계이다.
또한, 주파수-파수 필터 단계는 신호의 주파수-파수 대역은 통과시키고 잡음의 주파수-파수 대역은 감쇠시킴으로써 신호/잡음비(S/N)를 향상시키는 단계이다.
또한, 트레이스 편집 및 뮤트 단계는 기록이 되지 않았거나 잡음이 강하게 기록된 경우 기록별, 트레이스별, 측점별, 수신기별로 제거하고 극성이 바뀐 트레이스의 극성을 수정하는 단계이다.
또한, 공중점 분류 단계는 모든 트레이스의 발파 및 수신 위치 정보를 이용하여 공통 발파점 자료(common shot gather)를 공중점(common midpoint)으로 분류하는 단계이다.
또한, 속도 분석 단계는 공중점 자료를 이용하여 중합 속도를 구하는 단계이다.또한, 수직 경로 시차 보정(이하, NMO보정이라 함) 단계는 속도 함수를 이용하여 주시 t에 기록된 자료들을 영거리 주시 t0로 옮기는 단계이다.
또한, 중합 단계는 NMO보정한 자료를 중첩시켜 합하는 단계이다.
또한, 구조 보정 단계는 중합 단면상의 반사면을 제자리로 옮기는 단계이다.
상기 공중점 중합 방법은 정확한 지하 영상 정보를 얻을 수 있다는 장점이 있으나, 각각의 단계는 모듈(module)로 구성되어 있고, 처리 순서 및 처리 변수 선택에 따라 상이한 결과가 나타난다. 즉 분석자는 각 단계별로 최적 처리 변수 결정을 위한 시험 처리를 하여야 하며, 이러한 처리 단계에 많은 시간과 노력이 필요하다. 또한 분석자는 전문적 지식도 갖추고 있어야 한다. 아울러 상기 방법에 의할 때 얕은 지층 경계면에서 반사된 파는 이완, 역전, 압축, 중첩 등의 현상으로 자료가 심하게 왜곡되는 문제가 있다.이러한 공중점 중합 방법의 문제점을 해소하기 위한 종래의 다른 방법으로는 최적 옵셋 처리 방법이 있다.
최적 옵셋 처리 방법이란 상기 공중점 중합 방법이 자료 처리에 전문성이 필요하고 처리 시간이 많이 소요되는 문제점을 극복하기 위하여 공통 발파점 기록에서 가장 적절한 옵셋 트레이스 하나만을 사용하여 단면을 만드는 방법이다.
이 방법은 간단하고 신속하게 결과 단면을 얻을 수 있으며 더욱이 깊이에 따른 속도 함수를 전혀 모르는 상황에서도 훌륭히 지층 단면을 만들 수 있는 장점이 있다. 그러나 이 방법은 지층들의 깊이가 실제 깊이보다 훨씬 깊은 곳에 나타나서 전체적인 단면을 왜곡시킨다. 또한 천부 지층구조 파악을 위한 고분해능 반사법에서는 그리 문제가 되지 않으나, 석유 탐사 등 심부 지층구조 파악을 위해 이 방법을 사용한 경우 자료의 신호/잡음비(S/N)가 낮으면 깊은 지층 경계면을 잘 나타내지 못하는 문제가 있다.
따라서, 본 발명은 종래 기술의 제반 문제점을 해결하기 위하여 안출한 것으로서, 최적 옵셋이 반사면 깊이(반사파 주시)에 따라 변하게 함으로써 공중점 중합 방법이 가지고 있던 복잡성과 최적 옵셋 방법의 깊이 왜곡을 최소화시키는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치를 제공함에 그 목적이 있다.
또한, 본 발명의 목적은 공기파 도달 속도는 거의 일정하므로, 공기파 도달 이전의 적절한 창길이를 설정하여, 이 창 안의 기록만을 시간-거리 영역에서 수평 중합하는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치를 제공함에 그 목적이 있다.또한, 본 발명의 또 다른 목적은 공중점 분류 및 속도 정보 없이도 최적 옵셋 방법보다 훨씬 사실에 가까운 단면을 생성시키는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치를 제공함에 그 목적이 있다.
또한, 본 발명의 또 다른 목적은 공중점 모음 자료(common midpoint gather)에서의 시행을 통해 반사 깊이 변화에 따른 수평 이동 현상이 없어져서 공발점 모음(common shot gather)에서 보다 양호한 단면을 얻을 수 있으며, 속도 정보를 1-2곳에서만 알더라고 이 정보를 이용하여 동보정할 경우, 정확한 깊이를 알 수 있는 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치를 제공함에 그 목적이 있다.
도 1은 지질 구조를 영상화하기 위한 장치의 개략적인 구성도.
도 2는 자료 획득 장치를 사용하여 기초 정보 자료를 측정하는 일 실시예.
도 3은 본 발명에 따른 탄성파 자료 시변 최적 옵셋 처리 방법에 대한 순서도.
도 4a는 탄성파 신호 및 위치가 정해진 트레이스 창을 나타낸 예시도.
도 4b는 탄성파 신호 및 위치가 정해진 트레이스 창의 길이가 시간 또는 옵셋에 따라 변하는 것을 나타내는 예시도.
도 5는 NMO보정 과정을 나타낸 그래프.
도 6a 내지 도 6e(이하, 도 6이라 칭함)는 본 발명의 시변 최적 옵셋 처리 방법을 구현하기 위하여 포트란(FORTRAN)을 사용하여 구현한 일 실시예.
도 7은 컴퓨터의 시뮬레이션(simulation)을 통해 얻어진 결과물인 지하 영상 단면에 대한 예시도.
도 8a는 종래의 통상적인 탄성파 처리 방법인 공중점 중합 방법에 의한 결과물인 지하 영상 단면도.
도 8b는 도 8c 내지 도 8e에서 설명할 단면에서 측점 60번 트레이스와 대응되는 기록자료.
도 8c는 종래의 최적 옵셋 처리 방법에 의한 결과물인 지하 영상 단면도.
도 8d는 본 발명에 의한 속도 정보를 이용하지 않고 얻어진 결과를 나타낸 지하 영상의 단면도.
도 8e는 본 발명에 의한 속도 정보를 이용하여 얻어진 결과를 나타낸 지하 영상 단면도.
<도면의 주요부분에 대한 부호의 설명>
100 : 자료 획득 장치 110 : 발생부
120 : 수신부 130 : 기록부
150 : 자료 처리 장치 160 : 프로세서
170 : 디스플레이부
상술한 목적들을 달성하기 위하여 본 발명에 따르면, 반사법을 이용해 얻은 탄성파 자료를 처리하는 방법에 있어서, 자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력하고, 상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하는 단계 및 상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 단계를 포함하는 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템과 장치 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다.
또한, 상기 수평 중합된 자료를 도시하여 지하 영상 단면을 작성하는 것을 더 포함하는 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다.
또한, 상기 기초 정보는 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도이고, 상기 창은 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 상기 빠른 파의 속도보다 일정 시간만큼 빠른 부분인 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템과 장치 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다.또한, 상기 탄성파의 속도 정보를 이용한 경우 상기 탄성파의 속도 정보를 이용하지 않은 경우보다 정확한 정보를 얻을 수 있고, 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다. 상기 창은 트레이스 중 유용한 부분만을 정하는 부분으로 창의 길이가 시간 또는 옵셋에 따라 변할 수 있는 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템과 장치 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다.
또한, 상기 창 양쪽의 가장자리에서 상기 탄성파 트레이스의 갑작스런 단절에 의한 신호 왜곡 현상을 완화하기 위해 상기 창 양쪽에 테이퍼(taper)를 사용하는 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법, 상기 방법에 상응하는 시스템과 장치 및 이러한 방법을 수행할 수 있는 프로그램이 수록되어 있는 기록 매체가 제공된다.
이하, 본 발명에 따른 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치의 바람직한 실시예를 첨부 도면을 참조하여 보다 상세히 설명한다.
도 3은 본 발명에 따른 탄성파 자료 시변 최적 옵셋 처리 방법에 대한 순서도이다.
도 3에 도시한 바와 같이, 본 발명에 의할 경우, 먼저 자료 처리 장치의 프로세서는 자료 획득 장치로부터 입력된 기초 정보를 바탕으로 트레이스 창을 결정한다(410). 여기서 트레이스란 하나의 지오폰, 즉 수신기에 수신된 지면의 상대 운동 기록이다. 트레이스 창은 유용한 신호부분만을 정하는 부분으로, 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 이보다 빠른 일정한 시간만큼의 트레이스가 창 안의 부분으로 선택된다.
트레이스 창을 결정하는 방법에 대해 도 4a 및 도 4b를 참조하여 살펴보면 다음과 같다.
도 4a는 탄성파 신호 및 위치가 정해진 트레이스 창을 나타낸 예시도이다.
도 4a를 참조하여 설명하면, 레일리파의 속도는 약 210m/s이고 공기파의 속도는 약 335m/s이므로, 레일리파보다 빠른 공기파의 속도인 335m/s 의 기울기를 갖도록 창을 설정하였다. 창은 공기파 도달 시간에서 안전 구간 15 샘플(3ms)을 빼고 시작하였고, 신호의 우세주기인 6ms를 창 길이로 선정하였다.
도 4b는 탄성파 신호 및 위치가 정해진 트레이스 창의 길이가 시간 또는 옵셋에 따라 변하는 것을 나타내는 예시도이다.도 4b를 참조하여 설명하면, 창의 길이는 신호/잡음비가 좋은 경우에 신호의 한 파장 정도가 들어갈 수 있도록 정해진다. 예를 들어, 신호가 고주파인 경우 저주파인 경우보다 창의 길이를 짧게 해 주어야 한다. 또한 도4b에서는 창 길이를 시간에 따라 변하는 선택을 두어, 창 길이가 하부로 내려갈수록 즉, 시간이 증가함에 상응하여 증가시킬 수 있다. 하부가 내려갈수록 신호/잡음비(S/N)가 낮아져 동일한 창 길이를 유지할 경우 중합 효과가 감소하기 때문에 창 길이를 증가시킴으로써 중합 효과를 유지 또는 증가시킬 수 있다.
이 때, 창 가장자리의 갑작스런 자료 단절에 의한 신호 왜곡 등의 악영향 즉, 신호 왜곡 현상을 완화하기 위해 15 샘플(3ms)의 테이퍼(taper)를 사용하였다. 여기서 테이퍼(taper)란, 창의 경계를 기준으로 자료를 자를 경우 잘린 부분의 자료가 갑자기 0이 됨을 방지하기 위하여 가장자리 부분에 창을 설정하고 그 창 안에서 자료의 가중치가 0에서 1로 변하게 하는 경사 함수를 말한다.
다시 도 3을 참조하면, 프로세서는 다음으로 속도 정보 유무를 판단한다(420). 속도 정보의 유무 판단 결과, 탄성파 자료의 기초 정보에 속도 정보가 포함된 경우, NMO보정을 한다(425).
도 5는 NMO보정 과정을 나타낸 그래프이다.
도 5를 참조하여 설명하면, 510은 파원(S)에서 발생한 탄성파의 공중점 자료를 수신기(R)가 수신하는 일 실시예와 수신한 자료를 바탕으로 탄성파의 트레이스를 그린 그림이다. 여기서 TNMO는 NMO보정해야 할 시간을 나타낸다. 520은 NMO보정전을 나타낸 그림이고, 530은 NMO보정 후를 나타낸 그림이다. NMO보정이란 상기에서 언급한 바와 같이 속도 함수를 이용하여 주시 t에 기록된 자료들을 영거리 주시 t0로 옮기는 것으로써 속도 정보를 이용할 경우 520과 같이 정확한 깊이를 알 수 있다. 다시 도 3을 참조하면, 속도 정보의 유무 판단 결과 속도 정보가 없거나, 속도 정보가 있어 NMO보정을 한 경우, 프로세서는 트레이스 중 상기 창 안에 포함된 자료를 수평 중합한다(430). 여기서, 수평 중합이란 동일 시간대의 여러 트레이스 자료를 시간별로 더하는 과정을 말한다. 수평 중합 결과 지하 영상에 대한 시변 최적 옵셋 자료가 완성되게 된다.
그 후, 프로세서는 수평 중합에 의해 얻어진 시변 최적 옵셋 자료를 디스플레이부에 전송한다(460).
물론 전송하기 전, 신호/잡음비(S/N)를 향상시키기 위한 필터링 과정(442)과, 미약한 신호를 되살리는 자동 이득 조절 과정(444)을 거칠 수도 있으나 본 발명에서는 당업자라면 누구나 쉽게 알 수 있는 바 설명을 생략하기로 한다.프로세서로부터 지하 영상 단면에 대한 시변 최적 옵셋 자료를 얻은 디스플레이부는 모니터, 프린터 기타 출력장치를 통해 상기 자료를 출력한다.
다시 도 1을 참조하면, 종래의 지하 영상 단면을 얻기 위한 장치는 일반적으로 자료 획득 장치(100)와 자료 처리 장치(150)가 독립적으로 구분되어졌다. 자료 획득 장치(100)로부터 획득한 자료를 자료 처리 장치(150)의 프로세서(160)가 통상적인 방법 즉, 공중점 중합 방법으로 처리하는 경우에는 상기에서 언급한 바와 같이 통상 워크스테이션급 이상의 전문 컴퓨터를 이용하여야 하며, 10여단계의 복잡한 과정을 거쳐야 하므로 처리 시간이 길기 때문이다. 예를 들어, 처리 시간은 보통 짧게는 며칠에서 길게는 몇 주일 정도 소요된다.그런데 각각의 단계를 하나의 단계로 통합하여 단순화한 본 발명에 의할 경우, 아주 짧은 시간 안(보통, 3초에서 10초 사이)에 지하 영상 단면에 대한 결과를 얻을 수 있기 때문에 자료 획득 장치와 자료 처리 장치가 하나로 결합될 수도 있다. 그 결과 현장에서 즉시 지하 영상 단면에 대한 결과를 얻을 수 있다.
도 6은 본 발명의 시변 최적 옵셋 처리 방법을 구현하기 위하여 포트란(FORTRAN)을 사용하여 구현한 일 실시예를 나타낸다.
설명하기에 앞서, 발명이 속하는 기술 분야에서 통상의 지식을 가진 자라면, 쉽게 이해할 수 있는 부분은 설명에서 제외하도록 한다.
도 6을 참조하면, 이 포트란 언어로 작성한 코드는 탄성파 반사법 자료 처리용으로 개발된 프로그램으로 0.2ms 샘플률(dt)로 192ms(nt/dt)동안 12채널(nx)로 기록한 180개(nr) 자료(input.dat)를 읽어들여 처리한 후 중합 결과를 out1.dat란 파일로 출력한다.
아울러 속도 정보가 있을 경우, 속도 정보를 tv.dat에서 읽어 들인 후, NMO보정을 실시한 결과를 out2.dat에 기록한다. 또한 종래의 최적 옵셋 방법으로 구한 단면 자료도 비교 목록으로 출력토록 하여 결과를 optos.dat에 기록한다.
먼저 605 부분에서는 함수의 대한 각각의 변수, 데이터 등을 지정한다.
610 부분에서는 소스(source) 및 타겟(target)을 지정한다.
615 부분에서는 창 길이를 지정한다. 여기서 wd는 창 길이를 나타낸다.
620 부분에서는 자료를 입력하고 처리한다. 여기서 nr은 레코드 개수(180개)을 나타내고, nl은 지층 개수(3개)를 나타낸다.
구체적으로 살펴보면, 621 부분에서 속도 정보를 읽고, 622 부분에서 반사면의 샘플 번호에 따른 속도 정보를 tvmax에 저장한다. 또한, 623 부분에서는 입력된 자료를 읽는다. 624 부분에서는 최적 시변창을 결정하고 수평 중합을 실시하기 위해 부프로그램을 호출한다. 624 부분에서 부프로그램을 호출하면 프로세서는 660 부분으로 가서 최적 시변창을 결정하고 수평 중합을 실시한다. 그리고 그 결과를 627 부분에서 e1에 저장한다.
660 부분에서 661 부분은 최적 시변창을 결정하는 부분이고, 662 부분은 수평 중합을 실시하는 부분이다. 구체적으로 살펴보면, iwd는 창 경계(wd)를 시간 단위에서 샘플 단위로 바꾸어 준 것이고, itap는 창 경계에서의 모서리 효과를 없애기 위한 테이퍼(taper)의 길이를 나타낸다. 또한, nx는 채널 수, 즉 12개의 채널을 나타내고, xx는 각 채널의 옵셋을 나타내고, iair는 xx를 da로 나눈 값을 시간 단위에서 샘플 단위로 바꾼 후 1을 더한 값, 즉 공기파의 도달 시각(샘플수)을 나타낸다. 또한 nt는 전체 트레이스의 길이를 나타내고, iw1은 시작 부분의 창경계를 나타내고, iw2는 끝 부분의 창 경계를 나타낸다. 또한, is1은 iw1을 기준으로 한 창 시작 부분의 테이퍼 시작을 나타내고, is2는 iw1을 기준으로 한 창 시작 부분의 테이퍼 끝을 나타낸다. 또한 is2는 테이퍼 없이 통과하는 부분의 시작을 나타내고, ie2는 테이퍼 없이 통과하는 부분의 끝을 나타낸다. 또한, is3는 iw2를 기준으로한 창 끝 부분의 테이퍼 시작을 나타내고, is3는 iw2를 기준으로 한 창 끝 부분의 테이퍼 끝을 나타낸다.625 부분은 속도 정보를 이용할 경우에 필요한 부분으로, 샘플마다의 평균 제곱근 속도(vrms)를 계산하기 위한 부프로그램을 호출하는 부분이다. 625 부분에서 부프로그램을 호출하면 프로세서는 650 부분으로 가서 샘플마다의 평균 제곱근 속도를 구한다.
구체적으로 살펴보면, 650 부분중 651 부분은 Dix 구간 속도를 계산하는 식이고, 652 부분은 각 샘플마다의 평균 제곱근 속도를 계산하는 식이다.
626 부분은 최적 시변창을 결정하고 평균 제곱근 속도를 이용하여 NMO 보정한 후, 수평중합을 실시하기 위해 부프로그램을 호출하는 부분이다. 626 부분에서 부프로그램을 호출하면 프로세서는 665 부분으로 가서 625 부분에서 계산된 평균 제곱근 속도를 이용하여 NMO보정을 한 후 수평 중합을 실시한다. 그리고 그 결과를 628 부분에서 e2에 저장한다.
665 부분에서 667 부분은 최적 시변창을 결정하는 부분이고, 668 부분은 NMO 보정을 한 후, 수평 중합을 실시하는 부분이다.
665 부분 중 660 부분에서 설명한 부분과 다른 부분에 대해서만 설명하면, 665 부분의 666 부분에서 ddt1, ddt2, ddt3, ddt4는 각각 테이퍼 모서리에서의 NMO 보정량을 나타낸다.
도 3을 다시 참조하여 설명하면, 프로세서가 621 부분에서 속도 정보가 있는 지 판단하여(410) 속도 정보가 없는 경우 624 부분에서, 부프로그램인 660부분을 호출하여 수평 중합(430)을 수행하고, 속도 정보가 있는 경우 625 부분에서 부프로그램인 650 부분을 호출하여 평균 제곱근 속도를 구하고 665 부분을 호출하여 NMO 보정(425)을 하고 수평 중합(430)하게 된다.도 6을 다시 참조하면, 629 부분에서는 종래의 최적 옵셋 방법으로 구한 결과를 e3에 저장하는 단계이다. 종래의 최적 옵셋 처리 방법에 의한 자료도 저장하고 출력하는 이유는 본 발명에 의한 방법과 비교하기 위해서이다. 따라서 실제적으로 본 프로그램이 실제로 사용될 때는 본 부분이 생략될 수도 있다.
645 부분은 헤더 정보를 넣기 위해 부프로그램인 655 부분을 호출하고 자료를 출력하는 부분이다.
본 프로그램으로 기본적으로 기초 자료 입력 단계, 연산 단계, 출력 단계로 구성된다. 프로그램을 구현하는 방법상의 문제로서 입력 단계는 대화식, 배치 모드 모두 가능하다. 아울러 베이직(BASIC), 포트란(FORTRAN), 씨(C) 등 어떠한 사용 언어로도 본 프로그램을 구현할 수 있다. 또한, 결과물도 종래의 소프트웨어를 사용하여도 무방하다. 단지 본 명세서에서, 일 실시예로 사용한 것은, 배치 모드로 자료를 입력하였고, FORTRAN을 사용하였으며, 결과물은 제3소프트웨어인 WINSEIS를 통해 구현하였다.
이하에서는 본 발명의 이해를 돕기 위해서, 컴퓨터 시뮬레이션(simulation)을 통해 얻어진 결과물와 실제 측정에 의해 얻어진 결과물인 지하 영상 단면도에 대해 상세히 설명하기로 한다.
도 7은 컴퓨터의 시뮬레이션(simulation)을 통해 얻어진 결과물인 지하 영상 단면에 대한 예시도이다.도 7을 참조하여 살펴보면, 710 부분은 지하 영상에 대한이상적 처리 결과를 나타낸 지상 영상의 단면도이다. 본 단면도는 우세주기 4.2ms의 영위상 파형(zero-phase wavelet)을 반사면에 콘볼류션(convolution)한 결과이다.
720 부분은 종래의 통상적인 탄성파 자료 처리 방법인 공중점 중합 방법의 결과를 나타낸 지상 영상의 단면도이다. 본 단면도는 얕은 지층 경계면으로 합성한 반사면 1과 반사면 2는 이완된 파형의 모습으로 나타나며, 이상적인 결과와 비교했을 때 이상적인 결과보다 수직 해상도가 낮다.
730 부분은 종래의 최적 옵셋 처리 방법에 의한 결과를 나타낸 지하 영상의 단면도이다. 본 반면도에서는 반사면들이 실제 시간보다 지연되어 나타남을 알 수 있다. 이러한 현상은 특히 천부 반사면에서 심하다. 또한 54ms 부근의 공기파와 85ms 부근부터 뚜렷이 나타나는 레일리파의 파형을 확인할 수 있다.
740 부분은 본 발명에 의한 속도 정보를 이용하지 않고 얻어진 결과를 나타낸 지하 영상의 단면도이다. 본 단면도에서는 반사면이 잘 인지되어 있고 원치 않는 반사면이 나타나는 효과가 없어졌음을 알 수 있다. 그러나 이상적인 결과와 비교할 때 NMO 보정을 하지 않은 결과 반사면이 적절한 원래 시간대에 위치하지 않고 1.7ms에서 4.5ms 정도의 반사 시간이 증가되었음을 알 수 있다.
750 부분은 본 발명에 의한 속도 정보를 이용하여 얻어진 결과를 나타낸 지상 영상의 단면도이다. 본 단면도에서는 NMO보정 결과 반사면이 적절한 시간대에 위치하며, 파형도 이상적인 결과와 가장 흡사함을 알 수 있다.도 8a 내지 8e는 실제 측정을 통해 얻어진 결과물인 지하 영상 단면에 대한 예시도이다.
이 자료는 경상북도 경주시와 안강읍 사이의 기복이 거의 없는 논에서 획득되었다. 양산단층 주향과 거의 직각을 이루는 농로상에 553m의 측선을 설정하고, 192ms 기록 시간 동안 0.2ms 샘플률로 총 180개의 자료를 획득하였다. 신호 발생은 5kg 해머를 지면에 설치한 알루미늄판에 수직으로 타격하는 방식을 택하였으며, 수신부는 100Hz 고주파수용 지오폰 12개를 사용하였다. 기록부는 미국 Bison Instruments사의 GeoPro 8012A 12채널 디지털 탄성과 기록계를 이용하였으며, 150Hz의 아날로그(analog) 고주파 통과 필터를 사용하였다, 그리고 획득된 자료는 통신용 소프트웨어를 이용하여 9,600 bps의 속도를 가진 RS-232C 케이블로 현장에서 노트북 PC에 전송 보관하였다.
도 8a는 종래의 통상적인 탄성파 처리 방법인 공중점 중합 방법에 의한 결과물인 지하 영상 단면도이다.
도 8a를 참조하여 구체적으로 살펴보면, 총 180개의 공발점 자료를 동일 공중점으로 취합한 것이므로 이 단면에서의 트레이스 수는 373개이다. 층 1은 최소 7.8m에서 최대 12.6m 두께로 존재한다. Dix (1955)의 구간속도 계산식을 사용하여 구한 층 1의 구간속도는 최저 640m/s에서 최고 1020m/s이며, 계산된 구간속도로부터 이 층은 포화된 미고결 퇴적층으로 해석된다. 구간속도는 단면의 우측으로 가면서 일반적으로 감소하는 경향을 보이며, 인접한 지역의 소음분석 자료에서 인지되는 비교적 느린 속도의 얇은 층은 근거리 옵셋이 6m인 중합단면도 상에서는 큰 무브아웃으로 인한 이완뮤트로 대부분 제거되어 나타나지 않는다.
층 2는 4.7m에서 12.0m 정도의 두께로 존재하고, 구간 속도로써 판단한 이층은 포화된 미고결 퇴적층으로 총 3개의 구간으로 구분된다.
공중점 1번부터 190번까지로 분류되는 첫 번째 구간은 1638m/s 가량의 평균 구간속도와 평균 6.9m 정도의 두께로 존재하며 다른 구간에 비해서 높은 구간 속도의 분포를 나타내고 층 1과의 구간속도 차가 타구간에 비해서 크게 구해지므로 비교적 뚜렷하고 연장성이 양호한 반사면을 나타낸다.
이 미고결 퇴적층의 두께는 미소중력 및 지자기 연구를 통하여 김기영·이광자가 계산한 8.5m의 자기 기반암 평균 깊이와 비슷하다. 구간속도를 고려할 때 각각의 구간은 암상 종류가 동일하지 않은 포화된 미고결 퇴적층으로 해석된다. 두 번째 구간은 공중점 190번부터 254번 사이로 분류되며 중심으로 약 95m 정도에 걸쳐 존재하며 1347m/s 가량의 평균 구간속도와 평균 10.3m 정도의 두께를 갖는다. 이 구간을 절단하는 단층의 연장성은 인지되지 않으나, 층 1과의 경계면 연장성이 다른 구간과 비교해서 양호하지 않고 구간 하부 경계면에서 많은 수의 단층들이 발견됨으로써 파쇄대로 추정되며 단층작용으로 인한 파쇄 정도가 타구간보다 심한 것으로 해석된다.
세 번째 구간은 공중점 254번 이후의 구간에 해당되며 평균 구간속도는 1170m/s, 평균 두께는 8.7m 정도이며 다른 구간에 비해서 구간 속도가 비교적 낮게 구해진다. 층 1과 층 2 사이의 반사면 기울기는 공중점 254번에서부터 점진적으로 감소하며, 회절파의 존재, 반사면의 단절 혹은 상실, 중합속도의 갑작스런 변화 등의 단층 인지 요인에 따라 인지한 단층은 측선을 따라 총 16개가 존재한다.이 구간에서는 총 6개의 단층이 인지되고 그 중 3개 가량의 단층들이 파쇄대로 해석되는구간에서 발견된다. 이로써 이들 단층의 활동 시기는 그리 오래되지 않았을 것으로 해석되며, 대부분의 단층들이 70도 이상의 고경사를 나타내고 주변 지질을 참고할 때 주향 이동 단층으로 판단된다.
또한 파쇄대의 음향 기반암 위에 존재하는 반사면은 울림 현상(ringing)이 매우 심한 특성을 보이고 있는데, 이는 미고결 퇴적층과 음향 기반암 사이의 반사계수가 크기 때문인 것으로 해석된다.
층 3은 최소 11.0m에서 최대 19.1m 두께로 존재하며, 이 층의 구간속도는 1418m/s에서 2191m/s 사이의 값을 갖는다. 층 2와 마찬가지로 공중점 213번 부근의 파쇄대를 중심으로 왼쪽과 오른쪽 단면에서 구간별 층의 구간 속도차가 크게 나타나며, 왼쪽의 높은 구간 속도 구간과 오른쪽의 낮은 속도 구간으로 나눌 수 있다.
음향 기반암 상부가 풍화된 것으로 해석되는 이 층은 다른 층에 비해 중합속도 범위가 넓은 특징을 갖는다. 공중점 약 200번 부근까지 관찰되던 음향 기반암과 층 3 사이의 반사면이 오른쪽으로 가면서 나타나지 않는 이유는 반사면이 존재하지 않거나 또는 파원 세기가 부족한 것에서 기인한 것으로 예상된다.
도 8b, 8c, 8e의 가로축은 도 8a처럼 공중점이 아닌 파원을 발생시킨 타격점 또는 지오폰 수신점을 의미하며, 이하의 내용에서 측점이라 명명하였다. 때문에 간격이 3m인 총 180개 측점의 트레이스를 도시하였다. 도 8b는 도 8c 내지 도 8e에서 설명할 단면에서 측점 60번 트레이스와 대응되는 기록자료이다. 세 개의 반사면 중 첫 번째 반사면은 트레이스 1번의 26ms 부근에, 두 번째 반사면과 세 번째 반사면은 트레이스 2번에서 31ms 부근과, 38ms 부근에 각각 나타남을 알 수 있다.
도 8c는 종래의 최적 옵셋 처리 방법에 의한 결과물인 지하 영상 단면도이다.
도 8c를 참조하여 설명하면, 도 8c는 공발점 입력자료의 5번 트레이스만을 모아 얻은 지층 단면이다. 측점 60번에서 약 35ms에 나타나는 첫 번째 반사면과 약 41ms에 나타나는 두 번째 반사면을 인지할 수 있으며, 약 55ms에 공기파가 나타난다.
그러나 층이 더 깊어지는 100번에서 110번 사이, 50옵60ms 범위의 트레이스 상에서 공기파와 두 번째 반사면 신호가 서로 간섭을 일으켜 정확한 도달 시각을 알 수 없게 만든다.
입력 자료에서 첫 번째, 혹은 두 번째 트레이스(비교적 작은 옵셋 트레이스) 상의 반사면 도달 시각과 비교해 볼 때 약간의 시간 지연이 일어난 것을 알 수 있다. 측점 60번과 120번 사이, 70ms에서 110ms의 시간 구간 사이에서는 강하게 레일리파가 도달했음을 보여 주고 있다.
이 방법에서 사용된 5번 트레이스는 첫 번째 반사면 신호가 아주 미약하게 기록되어 있어서 이 트레이스만으로 도시한 최적 옵셋 단면에서는 첫 번째 반사면이 인지되지 않는다.
도 8d는 본 발명에 의한 속도 정보를 이용하지 않고 얻어진 결과를 나타낸 지하 영상의 단면도이다.
도 8d를 참조하여 설명하면, 우선 지층단면에서 공기파나 레일리파 도달 등으로 인해 나타나는 신호/잡음비(S/N)를 낮추는 요소들이 사라졌다. 또한, 최적 옵셋 단면에서 인지하기 어려운 첫 번째 반사면의 모습이 비교적 잘 보인다(측점 60번에서 29ms의 시간대에 나타나는 신호). 두 번째 반사면과 세 번째 반사면도 동일한 측점에서 각각 35ms와 41ms의 시간에 존재한다. 그러나 각각의 반사면 신호는 최적 옵셋 단면과 마찬가지로 실제 도달시각보다 약간 지연되어 나타남을 볼 수 있다.
도 8e는 본 발명에 의한 속도 정보를 이용하여 얻어진 결과를 나타낸 지하 영상 단면도이다.
도 8e를 참조하여 설명하면, 측점 60번에서 약 26ms, 31ms, 38ms에 반사 면들이 나타남을 인지할 수 있으며, 이들 도달시각은 실제 반사면이 존재하는 시간대에 거의 정확히 일치함을 인지할 수 있다. 다만 측점 96번 트레이스는 현장 자료 획득시 오류로 인해 양쪽의 트레이스와 일관성을 보이지 않는다.
상기에서와 같이 본 발명은 속도 정보를 이용하지 않고도 종래의 방법에 비해 훨씬 정확한 지하 영상 단면을 구할 수 있으며, 속도 정보를 이용할 경우 거의 실제 단면과 유사하게 지하 영상 단면을 구할 수 있다.
본 발명의 다른 실시예에 의하면, 본 발명에 의한 프로세서는 메모리에 저장된 것들과 같이, 본 발명의 일 실시예에 따라서 실시되는 프로세서를 실행할 수 있도록 설계된 응용 주문형 집적 회로(ASIC)으로 구현될 수 있다.
상기에서는 본 발명의 바람직한 실시예를 참조하여 설명하였지만, 해당 기술 분야에서 통상의 지식을 가진 자라면 하기의 특허 청구의 범위에 기재된 본 발명의 사상 및 영역으로부터 벗어나지 않는 범위 내에서 본 발명을 다양하게 수정 및 변경시킬 수 있음을 이해할 수 있을 것이다.
본 발명에 따른 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치는 종래의 통상적 처리 방법인 공중점 중합 방법이 거쳐야 하는 10여 단계의 처리 단계와 각 단계별로 최적 변수 결정을 위한 시험이 필요 없어서 처리 시간이 매우 짧은(기존에 1주일 이상 소요되던 것이 수초에 해결됨) 효과가 있다.
또한, 본 발명은 자료 처리의 전문적 지식 없이도 최종 지상 영상 단면을 손쉽게 작성할 수 있는 효과도 있다.
또한, 본 발명은 지하의 속도 분포에 관한 정보가 전혀 없어도, 종래의 최적 옵셋 방법으로 구한 지질 단면보다 월등히 사실에 가까운 단면을 얻을 수 있는 효과도 있다.
또한, 본 발명은 기존 조사비 중의 처리 과정 소요 경비가 불필요하여 조사비의 대폭 절감이 가능하게 되는 효과도 있다.
이와 같이 본 발명은 최근 대규모 토목, 건축, 지하수, 환경 조사 등에서 지하 지질 구조를 정밀하게 규명할 목적으로 고분해능 탄성파 반사법의 적용 사례가 증가하고 있으며, 아울러 최근 3차원적 방법에 치중되고 있는 추세를 감안할 때, 처리 시간이 짧아 현장에서 곧바로 최종 결과를 얻을 수 있고, 경비를 줄일 수 있으며, 전문 지식이 크게 필요 없어서, 전문 인력이 거의 없는 업계에서 충분히 활용할 수 있으며, PC급 컴퓨터만으로도 신속히 처리가 가능한 본 발명의 특징을 감안할 때 장래 활용이 크게 기대된다.
Claims (29)
- 반사법을 이용해 얻은 탄성파 자료를 처리하는 방법에 있어서,자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력받는 단계;상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하는 단계; 및상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 단계를 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 수평 중합된 자료를 도시하여 지하 영상 단면을 작성하는 단계를 더 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 기초 정보는 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도인것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 창은 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 상기 빠른 파의 속도보다 일정 시간만큼 빠른 부분인것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 탄성파의 속도 정보를 이용한 경우 상기 탄성파의 속도 정보를 이용하지 않은 경우보다 정확한 정보를 얻을 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 창은 트레이스 중 유용한 부분만을 정하는 부분으로 창의 길이가 시간 또는 옵셋에 따라 변할 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 제1항에 있어서,상기 창 양쪽의 가장자리에서 상기 탄성파 트레이스의 갑작스런 단절에 의한 신호 왜곡 현상을 완화하기 위해 상기 창 양쪽에 테이퍼(taper)를 사용하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 방법.
- 프로그램이 저장되어 있는 메모리;상기 메모리에 결합되어 상기 프로그램을 실행하는 프로세서를 포함하되,상기 프로세서는 상기 프로그램에 의해 자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력받는 단계;상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하는 단계; 및상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 단계를 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제8항에 있어서,상기 수평 중합된 자료를 도시하여 지하 영상 단면을 작성하는 단계를 더 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제8항에 있어서,상기 기초 정보는 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도인것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제8항에 있어서,상기 창은 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 상기 빠른 파의 속도보다 일정 시간만큼 빠른 부분인것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제8항에 있어서,상기 탄성파의 속도 정보를 이용한 경우 상기 탄성파의 속도 정보를 이용하지 않는 경우보다 정확한 정보를 얻을 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제12항에 있어서,상기 창은 트레이스 중 유용한 부분만을 정하는 부분으로 창의 길이가 시간 또는 옵셋에 따라 변할 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 제8항에 있어서,상기 창 양쪽의 가장자리에서 상기 탄성파 트레이스의 갑작스런 단절에 의한 신호 왜곡 현상을 완화하기 위해 상기 창 양쪽에 테이퍼(taper)를 사용하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 8항에 있어서,상기 프로세서는 상기 프로세서를 실행할 수 있도록 설계된 응용 주문형 직접 회로(ASIC)로 구현될 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 시스템.
- 탄성파 자료 시변 최적 옵셋 처리 방법을 수행하기 위하여 디지털 처리 장치에 의해 실행될 수 있는 명령어들의 프로그램이 유형적으로 구현되어 있으며 디지털 처리장치에 의해 판독될 수 있는 기록 매체에 있어서,상기 탄성파 자료 시변 최적 옵셋 처리 방법이,상기 프로그램에 의해 자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력받는 단계;상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하는 단계; 및상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 단계를 포함하는포함하는 것을 특징으로 하는 기록 매체.
- 제16항에 있어서,상기 수평 중합된 자료를 도시하여 지하 영상 단면을 작성하는 단계를 더 포함하는것을 특징으로 하는 기록매체.
- 제16항에 있어서,상기 기초 정보는 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도인것을 특징으로 하는 기록매체.
- 제16항에 있어서,상기 창은 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 상기 빠른 파의 속도보다 일정 시간만큼 빠른 부분인것을 특징으로 하는 기록매체.
- 제16항에 있어서,상기 탄성파의 속도 정보를 이용한 경우 상기 탄성파의 속도 정보를 이용하지 않은 경우보다 정확한 정보를 얻을 수 있는것을 특징으로 하는 기록매체.
- 제16항에 있어서,상기 창은 트레이스 중 유용한 부분만을 정하는 부분으로 창의 길이가 시간또는 옵셋에 따라 변할 수 있는것을 특징으로 하는 기록매체.
- 제16항에 있어서,상기 창 양쪽의 가장자리에서 상기 탄성파 트레이스의 갑작스런 단절에 의한 신호 왜곡 현상을 완화하기 위해 상기 창 양쪽에 테이퍼(taper)를 사용하는것을 특징으로 하는 기록매체.
- 반사법을 이용해 얻은 탄성파 자료를 처리하는 장치에 있어서,자료 획득 장치로부터 획득된 탄성파 자료의 기초 정보를 입력받는 수단;상기 기초 정보에 의해 얻어진 탄성파 트레이스의 시간 또는 옵셋(offset)에 따라 창(window)의 위치를 정하는 단계; 및상기 탄성파 트레이스(trace) 중 상기 창 안에 포함된 자료만을 수평 중합하는 단계를 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 수평 중합된 자료를 도시하여 지하 영상 단면을 작성하는 수단을 더 포함하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 기초 정보는 파원/수신기의 수직 및 수평 위치, 발파 간격, 수신기 간격, 기록 길이, 샘플률, 공기파 및 레일리파의 속도인 것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 창은 공기파나 레일리파 중 빠른 파의 속도를 기준으로, 상기 빠른 파의 속도보다 일정 시간만큼 빠른 부분인것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 탄성파의 속도 정보를 이용한 경우 상기 탄성파의 속도 정보를 이용하지 않은 경우보다 정확한 정보를 얻을 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 창은 트레이스 중 유용한 부분만을 정하는 부분으로 창의 길이가 시간 또는 옵셋에 따라 변할 수 있는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
- 제23항에 있어서,상기 창 양쪽의 가장자리에서 상기 탄성파 트레이스의 갑작스런 단절에 의한 신호 왜곡 현상을 완화하기 위해 상기 창 양쪽에 테이퍼(taper)를 사용하는것을 특징으로 하는 탄성파 자료 시변 최적 옵셋 처리 장치.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR20010045228A KR100412097B1 (ko) | 2001-07-26 | 2001-07-26 | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 |
JP2002188865A JP2003057357A (ja) | 2001-07-26 | 2002-06-28 | 地震波データの時変最適オフセット処理方法、処理システム、処理装置および記録媒体 |
US10/187,322 US20030031091A1 (en) | 2001-07-26 | 2002-06-28 | Method, system, and apparatus for processing seismic data by means of the time-varying optimum offset concept |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR20010045228A KR100412097B1 (ko) | 2001-07-26 | 2001-07-26 | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20010090027A true KR20010090027A (ko) | 2001-10-18 |
KR100412097B1 KR100412097B1 (ko) | 2003-12-24 |
Family
ID=19712579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR20010045228A KR100412097B1 (ko) | 2001-07-26 | 2001-07-26 | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20030031091A1 (ko) |
JP (1) | JP2003057357A (ko) |
KR (1) | KR100412097B1 (ko) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10895653B2 (en) | 2016-11-16 | 2021-01-19 | Korea Institute Of Geoscience And Mineral Resources | Method of first arrival picking of multi-channel seismic survey data considering sound source-receiver array |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080253226A1 (en) * | 2007-04-11 | 2008-10-16 | Stig Rune Lennart Tenghamn | System and method for marine seismic surveying |
US8380440B2 (en) * | 2008-06-02 | 2013-02-19 | Westerngeco L.L.C. | 3D residual binning and flatness error correction |
EP2864820A4 (en) | 2012-06-22 | 2015-06-17 | Schlumberger Technology Bv | SEISMIC DATA PROCESSING BY NON-LINEAR STACKING |
US9726772B2 (en) | 2013-10-17 | 2017-08-08 | Westerngeco L.L.C. | Iterative stacking of seismic image partitions |
KR101520399B1 (ko) * | 2013-12-02 | 2015-05-19 | 주식회사 코어밸런스 | 미소 위치 배열 중합을 이용한 지진파의 잡음제거 장치 및 방법 |
KR101842656B1 (ko) | 2015-10-05 | 2018-03-27 | 주식회사 에프에스 | 지층구조의 해석력 향상 기능을 갖는 공통반사표면법과 공통회절표면법을 이용한 탄성파 데이터의 복합처리방법 |
CN105445793B (zh) * | 2015-11-30 | 2018-04-06 | 中国石油天然气集团公司 | 一种确定坏道数据的方法及装置 |
CN106019370B (zh) * | 2016-05-06 | 2018-01-02 | 中国石油天然气股份有限公司 | 岩溶层断珠体油藏的井点确定方法及装置 |
CN107780400A (zh) * | 2016-08-24 | 2018-03-09 | 北京市水利规划设计研究院 | 一种大口径输水隧洞地基检测系统及方法 |
US11327188B2 (en) * | 2018-08-22 | 2022-05-10 | Saudi Arabian Oil Company | Robust arrival picking of seismic vibratory waves |
JP6734612B1 (ja) * | 2020-04-07 | 2020-08-05 | 株式会社ウオールナット | 空洞厚推定方法及びその装置 |
CN113075735B (zh) * | 2021-03-31 | 2024-03-08 | 中石化石油工程技术服务有限公司 | 基于直斜缆的速度拾取方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3651451A (en) * | 1969-10-23 | 1972-03-21 | Mobil Oil Corp | Continuous velocity estimation |
US4476551A (en) * | 1980-12-31 | 1984-10-09 | Mobil Oil Corporation | Selecting offset as a function of time and velocity to stack seismograms with high signal to noise ratio |
-
2001
- 2001-07-26 KR KR20010045228A patent/KR100412097B1/ko not_active IP Right Cessation
-
2002
- 2002-06-28 US US10/187,322 patent/US20030031091A1/en not_active Abandoned
- 2002-06-28 JP JP2002188865A patent/JP2003057357A/ja active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10895653B2 (en) | 2016-11-16 | 2021-01-19 | Korea Institute Of Geoscience And Mineral Resources | Method of first arrival picking of multi-channel seismic survey data considering sound source-receiver array |
Also Published As
Publication number | Publication date |
---|---|
US20030031091A1 (en) | 2003-02-13 |
KR100412097B1 (ko) | 2003-12-24 |
JP2003057357A (ja) | 2003-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Maggi et al. | An automated time-window selection algorithm for seismic tomography | |
US8687463B2 (en) | Method for processing at least two sets of seismic data | |
KR100412097B1 (ko) | 탄성파 자료 시변 최적 옵셋 처리 방법, 시스템 및 장치 | |
US20080270033A1 (en) | Methods of hydrocarbon detection using spectral energy analysis | |
CN113759424A (zh) | 基于频谱分解和机器学习的岩溶储层充填分析方法和系统 | |
Kluesner et al. | Practical approaches to maximizing the resolution of sparker seismic reflection data | |
CN104280777B (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
CA1261458A (en) | Sonic logging system | |
CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
WO2012071233A2 (en) | Identifying invalid seismic data | |
US5307268A (en) | Method for using P-S converted waves as a subsalt imaging tool | |
CA2606201C (en) | Method of processing seismic data for avo or avoa characterisation | |
EP0354813B1 (en) | Seismic data processing | |
CN108594301B (zh) | 一种具有差异特征的地震数据融合的方法及处理终端 | |
Marson-Pidgeon et al. | Source depth and mechanism inversion at teleseismic distances using a neighborhood algorithm | |
CN109490962A (zh) | 一种浅地层剖面非零偏移距的消除方法 | |
Mahvelati et al. | Comparison of dispersion-based analysis of surface waves and full waveform inversion in characterizing unknown foundations | |
Burdick | Estimation of the frequency dependence of Q from ScP and ScS phases | |
US20090052279A1 (en) | Method and system for determining azimuth parameters for seismic data processing | |
Mi et al. | Prestack depth imaging from topography with a Fourier method | |
CN104459788B (zh) | 一种确定地震数据采集参数的方法及系统 | |
CN115291300B (zh) | 一种基于复杂地质条件下的三维勘探分析系统及评估方法 | |
CN112824939B (zh) | 地震数据的剩余静校正方法、装置、设备及可读存储介质 | |
Poletto et al. | Seismic acquisition and processing of onshore dual fields by a reciprocal experiment | |
Pipatprathanporn et al. | Waveform modelling of hydroacoustic teleseismic earthquake records from autonomous Mermaid floats |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
N231 | Notification of change of applicant | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20081202 Year of fee payment: 6 |
|
LAPS | Lapse due to unpaid annual fee |