[go: up one dir, main page]

KR101476111B1 - 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치 - Google Patents

단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치 Download PDF

Info

Publication number
KR101476111B1
KR101476111B1 KR20140131556A KR20140131556A KR101476111B1 KR 101476111 B1 KR101476111 B1 KR 101476111B1 KR 20140131556 A KR20140131556 A KR 20140131556A KR 20140131556 A KR20140131556 A KR 20140131556A KR 101476111 B1 KR101476111 B1 KR 101476111B1
Authority
KR
South Korea
Prior art keywords
data
frequency
time
seismic reflection
analyzing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
KR20140131556A
Other languages
English (en)
Inventor
신성렬
정우근
하지호
정승재
Original Assignee
한국해양대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한국해양대학교 산학협력단 filed Critical 한국해양대학교 산학협력단
Priority to KR20140131556A priority Critical patent/KR101476111B1/ko
Application granted granted Critical
Publication of KR101476111B1 publication Critical patent/KR101476111B1/ko
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Geology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치가 제시된다. 본 발명에서 제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 시계열 자료의 입력자료를 생성하는 단계, 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계, 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 단계, 상기 주파수의 개별 주파수 자료를 추출하는 단계, 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계를 포함할 수 있다.

Description

단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치{Method and Apparatus for Separation Frequency Analysis of Seismic Reflection Data using Short-Time Fourier Transform}
본 발명은 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치에 관한 것이다.
탄성파 반사법 탐사란, 지층의 음향 임피던스 차이에 의해 지층 경계면 에서 반사되는 탄성파를 수진기에서 기록하여 지하구조에 대한 규명 등을 목적으로 사용되는 물리탐사 기법이다. 이러한 기법에서 사용하는 정보는 지층 경계면에서 반사되어 수진기에 기록되는 신호이다. 이때, 실제 기록에서는 지층 경계면으로부터의 반사파 이외에도 외부로부터의 잡음이 직접적으로 기록되거나, 다중 반사파, 기계적 잡음 등이 한계 기록될 수 있다. 따라서 신호와 잡음을 구분하거나, 기록으로부터 신호만을 추출해 내기 위한 절차는 탄성과 반사법 탐사에 필수적으로 잡고 있으며, 이를 위한 연구는 아직도 계속 되고 있는 실정이다.
종래 기술에 따른 탄성파 반사법 탐사 기법의 연구 및 프로그램의 경우에는 창 함수(window function)의 크기(window length)나 이동범위(window step) 등에 대한 설정이 자유롭지 못했으며, 이에 따라 증가되는 연산량(computational cost)의 증가를 감당하기 위하여 높은 사양의 컴퓨터가 요구되거나, 오랜 연산시간을 요구하게 되었다.
본 발명이 이루고자 하는 기술적 과제는 시간영역으로 기록되는 자료를 단시간 푸리에 변환(sort time Fourier transform)을 이용하여 시간-주파수 영역으로 변환하고 원하는 주파수를 선택적으로 기록하여, 신호와 잡음을 분리하는 방법 및 장치를 제공하는데 있다.
일 측면에 있어서, 본 발명에서 제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 시계열 자료의 입력자료를 생성하는 단계, 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계, 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 단계, 상기 주파수의 개별 주파수 자료를 추출하는 단계, 3차 스플라인 보간법(Spline interpolation method)을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계를 포함할 수 있다.
상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계는 푸리에 변환(Fourier transform)을 이용한 주파수 분석을 실시하여 중심 주파수 정보를 계산할 수 있다.
상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계는 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수들 중 하나를 선택하여 메인 알고리즘을 수행할 수 있다.
상기 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계는 상기 시간-주파수 영역의 자료에 대하여 트레이스(trace) 별로 단시간 푸리에 변환을 적용할 수 있다. 그리고, 상기 트레이스 별로 단시간 푸리에 변환을 적용한 결과 및 보간 된 결과를 이진 포맷으로 저장할 수 있다. 또한, 특정 주파수를 선택하여 주파수 별 반향 특성을 확인할 수 있는 스펙트럼 분해 자료 및 보간 된 결과를 이진 포맷으로 저장할 수 있다.
또 다른 일 측면에 있어서, 본 발명에서 제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치는 시계열 자료의 입력자료를 생성하는 입력자료 생성부, 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 시간구간별 주파수 분석부, 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 시간-주파수 영역 자료 생성부, 상기 주파수의 개별 주파수 자료를 추출하는 개별 주파수 자료 추출부, 3차 스플라인 보간법(Spline interpolation method)을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 보정부를 포함할 수 있다.
본 발명의 실시예들에 따르면 시간영역으로 기록되는 자료를 단시간 푸리에 변환(sort time Fourier transform)을 이용하여 시간-주파수 영역으로 변환하고 원하는 주파수를 선택적으로 기록하여, 신호와 잡음을 분리할 수 있다.
이를 통하여 특정 주파수에서 발행 가능한 이벤트(spectral anomaly) 등을 쉽게 발견하고, 지하구조의 특징에 대하여 쉽게 강조할 수 있다. 또한, 단시간 푸리에 변환에 사용되는 변수를 자유롭게 변환할 수 있고, 기존의 방식과 같은 연산과정뿐만 아니라, 3차원 스플라인 보간법을 이용하여 빠른 시간 안에 결과를 도출해 냄으로써, 연산 결과에 대한 빠른 추정이 가능하다. 이러한 결과는 현장 자료 등에 대한 적용성을 높여줄 뿐만 아니라, 학술적인 의미 또한 가질 것으로 예상된다.
도 1은 본 발명의 일 실시예에 따른 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법을 설명하기 위한 흐름도이다.
도 2는 본 발명의 일 실시예에 따른 창 함수의 이동구간 및 창 함수의 길이를 나타내는 도면이다.
도 3은 본 발명의 일 실시예에 따른 창 함수를 선택하는 과정을 설명하기 위한 도면이다.
도 4는 본 발명의 일 실시예에 따른 입력자료 생성 결과를 나타내는 도면이다.
도 5는 본 발명의 일 실시예에 따른 시간-주파수 영역 자료 생성 결과를 나타내는 도면이다.
도 6은 본 발명의 일 실시예에 따른 개별 주파수 성분의 시계열 자료 생성 및 3차 스플라인 보간법을 적용한 결과를 나타낸 도면이다.
도 7은 본 발명의 일 실시예에 따른 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치를 나타내는 도면이다.
이하, 본 발명의 실시 예를 첨부된 도면을 참조하여 상세하게 설명한다.
도 1은 본 발명의 일 실시예에 따른 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법을 설명하기 위한 흐름도이다.
탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 시계열 자료의 입력자료를 생성하는 단계(110), 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계(120), 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 단계(130), 상기 주파수의 개별 주파수 자료를 추출하는 단계(140), 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계(150)를 포함할 수 있다.
단계(110)에서, 시계열 자료의 입력자료를 생성할 수 있다. 다시 말해, 시계열 자료의 입력자료를 생성하여 원하는 주파수를 선택하기 위한 시간-주파수 영역으로 변환할 수 있다.
단계(120)에서, 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행할 수 있다. 다시 말해, 푸리에 변환을 이용한 기존 주파수 분석을 실시하여 중심 주파수(peak frequency) 정보를 계산할 수 있다. 이때, 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수(window function)들 중 하나를 선택하여 메인 알고리즘을 수행할 수 있다. 다시 말해, 창 함수(window shape)를 선택하고, 상기 창 함수의 길이(window length) 및 상기 창 함수의 이동구간 변수(window step)를 설정할 수 있다.
단계(130)에서, 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하고, 단계(140)에서, 상기 주파수의 개별 주파수 자료를 추출할 수 있다.
단계(150)에서, 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정할 수 있다. 이때, 상기 시간-주파수 영역의 자료에 대하여 트레이스 별로 단시간 푸리에 변환을 적용할 수 있다. 그리고, 그 결과와 보간 된 결과를 이진 포맷(binary format) 으로 저장할 수 있다. 또한, 특정 주파수를 선택하여 주파수 별 반향 특성을 확인할 수 있는 스펙트럼 분해(spectral decomposition) 자료 및 그 보관된 결과를 이진 포맷(binary format)으로 저장할 수 있다.
도 2는 본 발명의 일 실시예에 따른 창 함수의 이동구간 및 창 함수의 길이를 나타내는 도면이다.
도 2A는 창 함수의 이동구간 변수(window step)에 따른 시간 샘플 수(number of time samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이다.
도 2A(a)는 창 함수의 이동구간 변수(window step)가 5일 때 시간 샘플 수(number of time samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이고, 도 2A(b)는 창 함수의 이동구간 변수(window step)가 10일 때 시간 샘플 수(number of time samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이고, 도 2A(c)는 창 함수의 이동구간 변수(window step)가 15일 때 시간 샘플 수(number of time samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이다.
도 2B는 창 함수의 길이(window length)에 따른 샘플 수(number of samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이다. 도 2B(a)는 창 함수의 길이(window length)가 10일 때 샘플 수(number of samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이고, 도 2B(b)는 창 함수의 길이(window length)가 50일 때 샘플 수(number of samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이고, 도 2B(c)는 창 함수의 길이(window length)가 100일 때 샘플 수(number of samples) 및 진폭(Amplitude)의 관계를 나타내는 도면이다.
도 3은 본 발명의 일 실시예에 따른 창 함수를 선택하는 과정을 설명하기 위한 도면이다.
입력 자료에 대한 시간구간별 주파수 분석을 수행할 때, 푸리에 변환을 이용한 기존 주파수 분석을 실시하여 중심 주파수(peak frequency) 정보를 계산할 수 있다. 이때 도 3에 나타낸 것과 같이, 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수(window function)들 중 하나를 선택하여 메인 알고리즘을 수행할 수 있다. 다시 말해, 창 함수(window shape)를 선택하고, 상기 창 함수의 길이(window length) 및 상기 창 함수의 이동구간 변수(window step)를 설정할 수 있다.
이와 같은 방법을 통해 단시간 푸리에 변환에 사용되는 창 함수의 길이(window length)와 창 함수의 이동구간(window step) 변수를 자유롭게 설정함으로써, 각 변수가 결과에 미치는 영향에 대해서 확인할 수도 있다.
또한, 단시간 푸리에 변환의 창 함수 창 함수의 이동구간(window step) 변수를 자유롭게 설정할 수 있게 함으로써, 사용자가 제안하는 알고리즘의 계산시간 및 정확도를 편의에 따라 조절할 수 있다.
도 4는 본 발명의 일 실시예에 따른 입력자료 생성 결과를 나타내는 도면이다.
제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 도 4와 같이, 시계열 자료의 입력자료를 생성할 수 있다. 다시 말해, 시계열 자료의 입력자료를 생성하여 원하는 주파수를 선택하기 위한 시간-주파수 영역으로 변환할 수 있다.
도 5는 본 발명의 일 실시예에 따른 시간-주파수 영역 자료 생성 결과를 나타내는 도면이다.
제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 입력 자료에 대한 시간구간별 주파수 분석을 수행할 수 있다. 다시 말해, 푸리에 변환을 이용한 기존 주파수 분석을 실시하여 중심 주파수(peak frequency) 정보를 계산할 수 있다. 이때, 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수(window function)들 중 하나를 선택하여 메인 알고리즘을 수행할 수 있다. 다시 말해, 창 함수(window shape)를 선택하고, 상기 창 함수의 길이(window length) 및 상기 창 함수의 이동구간 변수(window step)를 설정할 수 있다.
그리고, 도 5와 같이 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축할 수 있고, 상기 주파수의 개별 주파수 자료를 추출할 수 있다.
도 6은 본 발명의 일 실시예에 따른 개별 주파수 성분의 시계열 자료 생성 및 3차 스플라인 보간법을 적용한 결과를 나타낸 도면이다.
제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정할 수 있다. 이때, 상기 시간-주파수 영역의 자료에 대하여 트레이스 별로 단시간 푸리에 변환을 적용할 수 있다. 그리고, 그 결과와 보간 된 결과를 이진 포맷(binary format)으로 저장할 수 있다. 또한, 특정 주파수를 선택하여 주파수 별 반향 특성을 확인할 수 있는 스펙트럼 분해(spectral decomposition) 자료 및 그 보관된 결과를 이진 포맷(binary format)으로 저장할 수 있다.
이와 같은 방법을 통해 단시간 푸리에 변환에 사용되는 창 함수의 길이(window length)와 창 함수의 이동구간(window step) 변수를 자유롭게 설정함으로써, 각 변수가 결과에 미치는 영향에 대해서 확인할 수도 있다.
또한, 단시간 푸리에 변환의 창 함수 창 함수의 이동구간(window step) 변수를 자유롭게 설정할 수 있게 함으로써, 사용자가 제안하는 알고리즘의 계산시간 및 정확도를 편의에 따라 조절할 수 있다. 예를 들어, 창 함수의 이동구간(window step) 변수가 1이 아닐 경우, 단시간 푸리에 변환 결과를 3차원 스플라인 보간법을 이용 및 내삽하여 시간영역에 대한 해상도를 상향 조절할 수 있다.
제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 고속 푸리에 변환 및 3차원 스플라인 보간법을 이용하여 종래기술의 단시간 푸리에 변환 방식에 비하여 계산속도가 더욱 빠르다. 또한, 입력 자료를 통하여 나이퀴스트 주파수(Nyquist frequency), 주파수 영역의 최소 단위(df) 및 중심주파수(peak frequency)를 계산하여 프린트(print)함으로써 특정 주파수에서의 반향 효과를 알아보기 위한 스펙트럼 분해(spectral decomposition) 기법에 대한 설계를 위해 추가적인 주파수 분석이 요구되지 않는다.
그리고, 제안하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법은 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수(window function)들 중 하나를 선택하여 메인 알고리즘을 수행함으로써 여러 가지 목적의 자료처리에 대한 적용성이 우수하다. 또한, 알고리즘 수행 결과에 대하여 수행된 트레이스의 위치와 추출한 주파수를 함께 표기함으로써 자료의 보관 및 사용이 편리하다.
도 7은 본 발명의 일 실시예에 따른 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치를 나타내는 도면이다.
탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치(700)는 입력자료 생성부(710), 시간구간별 주파수 분석부(720), 시간-주파수 영역 자료 생성부(730), 개별 주파수 자료 추출부(740), 보정부(750)를 포함할 수 있다.
입력자료 생성부(710)는 시계열 자료의 입력자료를 생성할 수 있다. 다시 말해, 시계열 자료의 입력자료를 생성하여 원하는 주파수를 선택하기 위한 시간-주파수 영역으로 변환할 수 있다.
시간구간별 주파수 분석부(720)는 상기 입력 자료에 대한 시간구간별 주파수 분석을 수행할 수 있다. 다시 말해, 푸리에 변환을 이용한 기존 주파수 분석을 실시하여 중심 주파수(peak frequency) 정보를 계산할 수 있다. 이때, 스퀘어(squared), 가우시안(gaussian), 해밍(hamming), 해닝(hanning), 터키(turkey) 윈도우(window)를 포함하는 복수의 창 함수(window function)들 중 하나를 선택하여 메인 알고리즘을 수행할 수 있다. 다시 말해, 창 함수(window shape)를 선택하고, 상기 창 함수의 길이(window length) 및 상기 창 함수의 이동구간 변수(window step)를 설정할 수 있다.
시간-주파수 영역 자료 생성부(730)는 상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축할 수 있다.
개별 주파수 자료 추출부(740)는 상기 주파수의 개별 주파수 자료를 추출할 수 있다.
보정부(750)는 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정할 수 있다. 이때, 상기 시간-주파수 영역의 자료에 대하여 트레이스 별로 단시간 푸리에 변환을 적용할 수 있다. 그리고, 그 결과와 보간 된 결과를 이진 포맷(binary format)으로 저장할 수 있다. 또한, 특정 주파수를 선택하여 주파수 별 반향 특성을 확인할 수 있는 스펙트럼 분해(spectral decomposition) 자료 및 그 보관된 결과를 이진 포맷(binary format)으로 저장할 수 있다.
이상에서 설명된 장치는 하드웨어 구성요소, 소프트웨어 구성요소, 및/또는 하드웨어 구성요소 및 소프트웨어 구성요소의 조합으로 구현될 수 있다. 예를 들어, 실시예들에서 설명된 장치 및 구성요소는, 예를 들어, 프로세서, 콘트롤러, ALU(arithmetic logic unit), 디지털 신호 프로세서(digital signal processor), 마이크로컴퓨터, FPA(field programmable array), PLU(programmable logic unit), 마이크로프로세서, 또는 명령(instruction)을 실행하고 응답할 수 있는 다른 어떠한 장치와 같이, 하나 이상의 범용 컴퓨터 또는 특수 목적 컴퓨터를 이용하여 구현될 수 있다. 처리 장치는 운영 체제(OS) 및 상기 운영 체제 상에서 수행되는 하나 이상의 소프트웨어 애플리케이션을 수행할 수 있다. 또한, 처리 장치는 소프트웨어의 실행에 응답하여, 데이터를 접근, 저장, 조작, 처리 및 생성할 수도 있다. 이해의 편의를 위하여, 처리 장치는 하나가 사용되는 것으로 설명된 경우도 있지만, 해당 기술분야에서 통상의 지식을 가진 자는, 처리 장치가 복수 개의 처리 요소(processing element) 및/또는 복수 유형의 처리 요소를 포함할 수 있음을 알 수 있다. 예를 들어, 처리 장치는 복수 개의 프로세서 또는 하나의 프로세서 및 하나의 콘트롤러를 포함할 수 있다. 또한, 병렬 프로세서(parallel processor)와 같은, 다른 처리 구성(processing configuration)도 가능하다.
소프트웨어는 컴퓨터 프로그램(computer program), 코드(code), 명령(instruction), 또는 이들 중 하나 이상의 조합을 포함할 수 있으며, 원하는 대로 동작하도록 처리 장치를 구성하거나 독립적으로 또는 결합적으로(collectively) 처리 장치를 명령할 수 있다. 소프트웨어 및/또는 데이터는, 처리 장치에 의하여 해석되거나 처리 장치에 명령 또는 데이터를 제공하기 위하여, 어떤 유형의 기계, 구성요소(component), 물리적 장치, 가상 장치(virtual equipment), 컴퓨터 저장 매체 또는 장치, 또는 전송되는 신호 파(signal wave)에 영구적으로, 또는 일시적으로 구체화(embody)될 수 있다. 소프트웨어는 네트워크로 연결된 컴퓨터 시스템 상에 분산되어서, 분산된 방법으로 저장되거나 실행될 수도 있다. 소프트웨어 및 데이터는 하나 이상의 컴퓨터 판독 가능 기록 매체에 저장될 수 있다.
실시예에 따른 방법은 다양한 컴퓨터 수단을 통하여 수행될 수 있는 프로그램 명령 형태로 구현되어 컴퓨터 판독 가능 매체에 기록될 수 있다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체에 기록되는 프로그램 명령은 실시예를 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다. 상기된 하드웨어 장치는 실시예의 동작을 수행하기 위해 하나 이상의 소프트웨어 모듈로서 작동하도록 구성될 수 있으며, 그 역도 마찬가지이다.
이상과 같이 실시예들이 비록 한정된 실시예와 도면에 의해 설명되었으나, 해당 기술분야에서 통상의 지식을 가진 자라면 상기의 기재로부터 다양한 수정 및 변형이 가능하다. 예를 들어, 설명된 기술들이 설명된 방법과 다른 순서로 수행되거나, 및/또는 설명된 시스템, 구조, 장치, 회로 등의 구성요소들이 설명된 방법과 다른 형태로 결합 또는 조합되거나, 다른 구성요소 또는 균등물에 의하여 대치되거나 치환되더라도 적절한 결과가 달성될 수 있다.
그러므로, 다른 구현들, 다른 실시예들 및 특허청구범위와 균등한 것들도 후술하는 특허청구범위의 범위에 속한다.

Claims (7)

  1. 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법에 있어서,
    시계열 자료의 입력자료를 생성하는 단계;
    상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계;
    상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 단계;
    상기 주파수의 개별 주파수 자료를 추출하는 단계; 및
    3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계
    를 포함하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  2. 제1항에 있어서,
    상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계는,
    푸리에 변환을 이용한 주파수 분석을 실시하여 중심 주파수 정보를 계산하는
    탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  3. 제1항에 있어서,
    상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 단계는,
    스퀘어, 가우시안, 해밍, 해닝, 터키 윈도우를 포함하는 복수의 창 함수들 중 하나를 선택하여 메인 알고리즘을 수행하는
    탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  4. 제1항에 있어서,
    상기 3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 단계는,
    상기 시간-주파수 영역의 자료에 대하여 트레이스 별로 단시간 푸리에 변환을 적용하는
    탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  5. 제4항에 있어서,
    상기 트레이스 별로 단시간 푸리에 변환을 적용한 결과 및 보간 된 결과를 이진 포맷으로 저장하는
    탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  6. 제4항에 있어서,
    특정 주파수를 선택하여 주파수 별 반향 특성을 확인할 수 있는 스펙트럼 분해 자료 및 보간 된 결과를 이진 포맷으로 저장하는
    탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법.
  7. 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치에 있어서,
    시계열 자료의 입력자료를 생성하는 입력자료 생성부;
    상기 입력 자료에 대한 시간구간별 주파수 분석을 수행하는 시간구간별 주파수 분석부;
    상기 주파수의 전 주파수 범위 내 분석을 위한 시간-주파수 영역 자료를 구축하는 시간-주파수 영역 자료 생성부;
    상기 주파수의 개별 주파수 자료를 추출하는 개별 주파수 자료 추출부; 및
    3차 스플라인 보간법을 이용하여 상기 개별 주파수 자료의 샘플 수를 보정하는 보정부
    를 포함하는 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 장치.
KR20140131556A 2014-09-30 2014-09-30 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치 Expired - Fee Related KR101476111B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR20140131556A KR101476111B1 (ko) 2014-09-30 2014-09-30 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR20140131556A KR101476111B1 (ko) 2014-09-30 2014-09-30 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치

Publications (1)

Publication Number Publication Date
KR101476111B1 true KR101476111B1 (ko) 2014-12-24

Family

ID=52679837

Family Applications (1)

Application Number Title Priority Date Filing Date
KR20140131556A Expired - Fee Related KR101476111B1 (ko) 2014-09-30 2014-09-30 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치

Country Status (1)

Country Link
KR (1) KR101476111B1 (ko)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104849754A (zh) * 2015-05-20 2015-08-19 中国地质大学（武汉） 用于金属破裂声发射信号检测的归一化分析方法
CN105044772A (zh) * 2015-08-06 2015-11-11 中国石油天然气集团公司 一种实现时变谱模拟反褶积的方法和装置
KR20160095690A (ko) * 2015-02-03 2016-08-12 한국해양대학교 산학협력단 단일채널 탄성파 탐사자료로부터 지하 속도정보 도출을 위한 탄성파 파형역산 방법 및 장치
KR20160116459A (ko) * 2015-03-30 2016-10-10 한양대학교 산학협력단 탄성파 트레이스 내삽을 위한 방법 및 그 장치
CN106990449A (zh) * 2017-04-20 2017-07-28 中国地质大学（武汉） 一种基于多维数据的陆相层序地层格架构建表征方法
KR101891269B1 (ko) * 2018-03-05 2018-08-27 주식회사 이에스피 구조물 결함 및 내층 변화 탐지를 위한 uwb 내부 변화 감지 시스템 장치
CN111045079A (zh) * 2019-12-20 2020-04-21 核工业北京地质研究院 一种增强地震反射特征的数据处理方法
KR20230017702A (ko) 2021-07-27 2023-02-06 브레인소프트주식회사 해석적 방법에 기반한 디제이 변환 주파수 추출 방법

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0572347A (ja) * 1991-09-13 1993-03-26 Nec Corp 音波探査装置
KR20040005288A (ko) * 2002-07-09 2004-01-16 신용준 시간-주파수 영역 반사파 처리장치 및 방법
KR20070113517A (ko) * 2006-05-24 2007-11-29 단국대학교 산학협력단 수중 초음파 원격이탈 제어 시스템
KR100923685B1 (ko) * 2009-07-08 2009-10-28 삼성탈레스 주식회사 근사 고속 퓨리에 변환을 기반한 능동 소나 탐지기에서 표적 위치를 탐지하는 장치 및 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0572347A (ja) * 1991-09-13 1993-03-26 Nec Corp 音波探査装置
KR20040005288A (ko) * 2002-07-09 2004-01-16 신용준 시간-주파수 영역 반사파 처리장치 및 방법
KR20070113517A (ko) * 2006-05-24 2007-11-29 단국대학교 산학협력단 수중 초음파 원격이탈 제어 시스템
KR100923685B1 (ko) * 2009-07-08 2009-10-28 삼성탈레스 주식회사 근사 고속 퓨리에 변환을 기반한 능동 소나 탐지기에서 표적 위치를 탐지하는 장치 및 방법

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101695080B1 (ko) 2015-02-03 2017-01-11 한국해양대학교 산학협력단 단일채널 탄성파 탐사자료로부터 지하 속도정보 도출을 위한 탄성파 파형역산 방법 및 장치
KR20160095690A (ko) * 2015-02-03 2016-08-12 한국해양대학교 산학협력단 단일채널 탄성파 탐사자료로부터 지하 속도정보 도출을 위한 탄성파 파형역산 방법 및 장치
KR101667481B1 (ko) 2015-03-30 2016-10-18 한양대학교 산학협력단 탄성파 트레이스 내삽을 위한 방법 및 그 장치
KR20160116459A (ko) * 2015-03-30 2016-10-10 한양대학교 산학협력단 탄성파 트레이스 내삽을 위한 방법 및 그 장치
CN104849754B (zh) * 2015-05-20 2017-06-09 中国地质大学（武汉） 用于金属破裂声发射信号检测的归一化分析方法
CN104849754A (zh) * 2015-05-20 2015-08-19 中国地质大学（武汉） 用于金属破裂声发射信号检测的归一化分析方法
CN105044772A (zh) * 2015-08-06 2015-11-11 中国石油天然气集团公司 一种实现时变谱模拟反褶积的方法和装置
CN106990449A (zh) * 2017-04-20 2017-07-28 中国地质大学（武汉） 一种基于多维数据的陆相层序地层格架构建表征方法
CN106990449B (zh) * 2017-04-20 2019-12-17 中国地质大学（武汉） 一种基于多维数据的陆相层序地层格架构建表征方法
KR101891269B1 (ko) * 2018-03-05 2018-08-27 주식회사 이에스피 구조물 결함 및 내층 변화 탐지를 위한 uwb 내부 변화 감지 시스템 장치
CN111045079A (zh) * 2019-12-20 2020-04-21 核工业北京地质研究院 一种增强地震反射特征的数据处理方法
CN111045079B (zh) * 2019-12-20 2021-11-30 核工业北京地质研究院 一种增强地震反射特征的数据处理方法
KR20230017702A (ko) 2021-07-27 2023-02-06 브레인소프트주식회사 해석적 방법에 기반한 디제이 변환 주파수 추출 방법

Similar Documents

Publication Publication Date Title
KR101476111B1 (ko) 단시간 푸리에 변환기법을 이용한 탄성파 반사법 탐사자료의 개별 주파수 특성 분석 방법 및 장치
US8767978B2 (en) System and method for processing sound signals implementing a spectral motion transform
JP6236282B2 (ja) 異常検出装置、異常検出方法、及びコンピュータ読み取り可能な記憶媒体
JP2017090912A (ja) ニューラルネットワークトレーニング装置及び方法と、音声認識装置及び方法
JP2019535448A5 (ko)
CN107340536B (zh) 地震反演剥离煤层强反射的方法和装置
BR102013012369A2 (pt) Método e aparelho para remoção de fantasmas pré-empilhamento de dados sísmicos
CN109644304B (zh) 混响环境的源分离
JP2015138100A5 (ko)
KR101483513B1 (ko) 음원위치추적장치 및 음원위치추적방법
KR101695080B1 (ko) 단일채널 탄성파 탐사자료로부터 지하 속도정보 도출을 위한 탄성파 파형역산 방법 및 장치
KR20170101614A (ko) 분리 음원을 합성하는 장치 및 방법
CN108897036A (zh) 一种地震数据处理方法及装置
JP2008261720A (ja) アンビギュイティ処理装置
JP2016093492A5 (ja) 情報処理装置、情報処理方法、およびプログラム
WO2016197629A1 (en) System and method for frequency estimation
Sahnoun et al. Sparse modal estimation of 2-D NMR signals
KR102044520B1 (ko) 음성 존재 구간 판별 장치 및 방법
Candy Environmentally adaptive processing for shallow ocean applications: A sequential Bayesian approach
KR102084315B1 (ko) 주파수 검출 및 위상 차이 추정에 기반한 선박 기적소리의 방향 탐지 방법 및 그 장치
US20160170058A1 (en) Method and system of seismic data processing
JP6559579B2 (ja) 観測地震の距離パラメータを推定する方法ならびに距離パラメータ推定プログラムおよび距離パラメータ推定プログラムを記録したコンピュータ読み取り可能な記録媒体
KR101838408B1 (ko) 오디오 신호 처리 방법 및 장치
CN110673207B (zh) 一种高频重建方法、装置及计算机存储介质
Briani et al. Validated exponential analysis for harmonic sounds

Legal Events

Date Code Title Description
PA0109 Patent application

Patent event code: PA01091R01D

Comment text: Patent Application

Patent event date: 20140930

PA0201 Request for examination
PA0302 Request for accelerated examination

Patent event date: 20141008

Patent event code: PA03022R01D

Comment text: Request for Accelerated Examination

Patent event date: 20140930

Patent event code: PA03021R01I

Comment text: Patent Application

E701 Decision to grant or registration of patent right
PE0701 Decision of registration

Patent event code: PE07011S01D

Comment text: Decision to Grant Registration

Patent event date: 20141030

GRNT Written decision to grant
PR0701 Registration of establishment

Comment text: Registration of Establishment

Patent event date: 20141217

Patent event code: PR07011E01D

PR1002 Payment of registration fee

Payment date: 20141217

End annual number: 3

Start annual number: 1

PG1601 Publication of registration
FPAY Annual fee payment

Payment date: 20191209

Year of fee payment: 6

PR1001 Payment of annual fee

Payment date: 20191209

Start annual number: 6

End annual number: 6

PR1001 Payment of annual fee

Payment date: 20210617

Start annual number: 7

End annual number: 7

PC1903 Unpaid annual fee

Termination category: Default of registration fee

Termination date: 20220928