• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2005년도 춘계학술대회
• /
• pp.638-642
• /
• 2005

For gas hydrate exploration, long offset multichannel seismic data acquired using by the 4km streamer length in Ulleung basin of the East Sea. The dataset was processed to define the BSRs (Bottom Simulating Reflectors) and to estimate the amount of gas hydrates. Confirmation of the presence of Bottom Simulating reflectors (BSR) and investigation of its physical properties from seismic section are important for gas hydrate detection. Specially, faster interval velocity overlying slower interval velocity indicates the likely presences of gas hydrate above BSR and free gas underneath BSR. In consequence, estimation of correct interval velocities and analysis of their spatial variations are critical processes for gas hydrate detection using seismic reflection data. Using Dix's equation, Root Mean Square (RMS) velocities can be converted into interval velocities. However, it is not a proper way to investigate interval velocities above and below BSR considering the fact that RMS velocities have poor resolution and correctness and the assumption that interval velocities increase along the depth. Therefore, we incorporated Migration Velocity Analysis (MVA) software produced by Landmark CO. to estimate correct interval velocities in detail. MVA is a process to yield velocities of sediments between layers using Common Mid Point (CMP) gathered seismic data. The CMP gathered data for MVA should be produced after basic processing steps to enhance the signal to noise ratio of the first reflections. Prestack depth migrated section is produced using interval velocities and interval velocities are key parameters governing qualities of prestack depth migration section. Correctness of interval velocities can be examined by the presence of Residual Move Out (RMO) on CMP gathered data. If there is no RMO, peaks of primary reflection events are flat in horizontal direction for all offsets of Common Reflection Point (CRP) gathers and it proves that prestack depth migration is done with correct velocity field. Used method in this study, Tomographic inversion needs two initial input data. One is the dataset obtained from the results of preprocessing by removing multiples and noise and stacked partially. The other is the depth domain velocity model build by smoothing and editing the interval velocity converted from RMS velocity. After the three times iteration of tomography inversion, Optimum interval velocity field can be fixed. The conclusion of this study as follow, the final Interval velocity around the BSR decreased to 1400 m/s from 2500 m/s abruptly. BSR is showed about 200m depth under the seabottom

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2007년도 공동학술대회 논문집
• /
• pp.151-156
• /
• 2007

In December 2006, 2D surface streamer and Ocean Bottom Seismometer (OBS) data were acquired in the Ulleung basin in Korea where strong Bottom Simulating Reflectors (BSR) were shown as a result of 2D and 3D multichannel (MCS) reflection survey. The aim of this study is to provide another reliable source for estimating P wave velocity around BSR depth using OBS data in addition to velocity information from 2D surface seismic data. Four OBSs were deployed and four 20-km shot lines which pass two OBSs respectively were designed. To derive P wave velocity profile, interactive interval velocity analysis using ${\tau}$-p trajectory matching method (Kumar, 2005) was used for OBS data and semblance analysis was used for surface data. The seismic profiles cross the OBS instruments in two different directions yield recordings for four different azimuths. This raised the confidence for the results. All velocity profiles in the vicinity of BSR depth of four OBS sites show almost definite velocity changes which we could consider as upper BSR and free gas layer. Making comparison between velocity from OBS and that from 2D seismic semblance velocity analysis gives consistency in result.

• 지구물리와물리탐사
• /
• 제2권4호
• /
• pp.184-190
• /
• 1999

가스 하이드레이트는 전세계적으로 새로운 에너지 자원으로 활용 가능성을 포함하고 있어 연구가 활발하게 진행되고 있다. 한국자원연구소에서는 1997년 부터 동해에서 메탄 하이드레이트 부존 잠재력 규명을 위한 탄성파 탐사를 하고 있다. 탄성파 자료에서 하이드레이트 부존을 의미하는 일반적인 특성은 해저면과 평행하게 나타나는 BSR(Bottom Simulating Reflection)과 BSR 상부에서 보이는 진폭감소 그리고 BSR 하부에서 보이는 진폭증가와 구간속도의 감소 그리고 BSR에서 반사파의 역전현상 등이 있다. 따라서 위와 같은 하이드레이트 부존특성을 탐지하기 위한 목적으로 실시되는 자료처리는 진 진폭을 유지하는 자료처리, 정밀 속도분석 및 AVO분석 등이 요구된다. 본 연구는 1998년 동해에서 취득된 탄성파 탐사자료를 처리하여 하이드레이트 부존 가능성을 확인하고자 하였다. 적용된 자료처리 공정은 구형확산 보정과 주파수 필터링, 공심점 분류, 정밀 속도분석 공정 등이다. AVO분석은 이용된 현장자료가 AVO를 분석할 정도의 입사각을 유지하고 있지 않아 제외하였다. 정밀 속도분석은 반복적으로 속도 스펙트럼을 구하는 방법으로 정확한 중합속도 결정이 가능한 XYA를 이용하였으며 자료처리의 모든 공정은 국내고유의 탄성파 자료처리 소프트웨어인 Geobit 2.9.5 를 이용하였다. 자료처리 결과 음원위치 $1650\~1900$에서 해저면으로 부터 약 $367\~477m$ 깊이(왕복주시 약 1800ms)에 해저면과 평행하게 발달한 BSR을 확인할 수 있었으며, BSR부근에서 구간속도 감소 뿐만 아니라 해저면 반사파의 위상과 반대인 반사파 역전현상도 확인할 수 있었다.

• 자원환경지질
• /
• 제29권6호
• /
• pp.739-744
• /
• 1996

Kite는 조사지역을 한번 횡단함으로써 3차원 지하영상을 얻기 위한, 새로이 개발된 단일채널(single-channel) 탄성파 취득시스템이다. Kite시스템은 진행방향에 수직으로 전개되는 하이드로폰(hydrophone)과 하이드로폰의 한쪽 끝에 설치하는 에너지원(source)인 에어건(air-gun)으로 구성된다. 3차원 지하영상은 기존 다중채널 반사자료 처리방법을 사용하며, 공통 소스점 모음(CSS gather)을 딕스식을 이용하여 지하층들의 평균 구간속도를 구하게 된다. 이 구간속도와 정규 고결응력으로 부터 물리상수들을 선택적으로 이용하여 전단계수, 공극율, 횡파속도 등을 구하였다. 본 시스템은 북대서양의 대륙붕에서 실시되었으며 좋은 결과를 얻을 수 있었다.

• 지구물리와물리탐사
• /
• 제5권3호
• /
• pp.141-144
• /
• 2002

1년 시차를 두고 획득한 3차원 탄성파 자료를 분석한 결과 매질의 탄성파 속도 변화를 관찰할 수 있었다. 이 1년 기간 중에 연구지역 내에 있는 철로 지반에 대한 지반보강 공사가 실시되었다. 두 탄성파 자료중합 단면으로부터 첫번째 반사면의 시간을 비교한 결과 연구지역의 상당 부분에서 매질의 속도가 감소된 것을 관찰할 수 있었으나 철로지반 부근에서는 속도감소를 인지할 수 없었다.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2006년도 춘계학술대회
• /
• pp.377-381
• /
• 2006

A seismic survey for gas hydrate have performed over the East sea by the KIGAM since 1997. General indicator of gas hydrate in seismic data is commonly inferred from the BSR(Bottom Simulating Reflector) that occurred parallel to the sea floor, amplitude decrease at the top of the BSR, amplitude blanking at the bottom of the BSR, decrease of the interval velocity and the reflection phase reversal at the BSR. In this paper we had analyzed optimum parameters of the field data to detect the 9as hydrate. Shot delay correction is applied 95ms, spherical divergence correction is applied velocity library 3, bandpass filter is applied 25-30-115-120Hz deconvolution operator length is applied 60ms, lag is 6ms and accurate velocity analysis NMO correction, stack is performed. Geobit 2.11.0 developed by the KIGAM was used for all data processing. Processing results say that the BSR occurred parallel to the sea floor were shown at 3,150m/s of two way travel time from the sea floor through shot point 5,000-5,610, and identified the interval velocity decrease around BSR and the reflection phase reversal corresponding to the reflection at the sea floor.

• 한국지반공학회지:지반
• /
• 제14권3호
• /
• pp.95-108
• /
• 1998

지반의 S파 및 P파의 깊이에 따른 변화를 원위치에서 측정하기 위하여 다운흘 시험 (downhole testing)과 SCPT (seismic CPT) 등이 널리 사용되어 왔다. 다운홀 시험과 SCPT는 경제성, 운용의 용이성, 발진원의 단순성 등의 측면에서 효율적이기 때문에, 현재 지반조사에서 그 사용빈도가 더욱 증가하고 있는 추세이다. 특히 최근에는 다운흘과 SCPT의 자료 분석을 자동화하기 위한 노력의 일환으로 interval measurements의 기법이 활용되고 있는데, 현재 이에 대한 적절한 역산해석 (inversion analysis) 기법이 없는 형편이다. 따라서, 본 논문에서는 다운홀이나 SCPT의 interval measurements를 분석하기 위한 새로운 역산해석 기법을 제안하였다. 제안된 역산해석 기법의 정모델링(forward modeling)에서는 탄성파의 전파를 Snell의 법칙에 의거하여 굴절.반사되는 현상을 고려하였곡, 역산해석을 위해서는 최대공산법 (maximum likelihood method)을 적용하였다. 그리고, 본 논문에서 제안한 역산해석 기법의 검증을 위하여, 하나의 S파 주상도를 가정하고 이에 대하여 다운흘 시험을 모사하였다. 이론적으로 수행한 다운홀 시험 결과에 대하여 기존의 비 역산해석 방법과 본 논문에서 제안한 역산해 석 기법에 의해서 S파 주상도를 추정하였는데, 그 결과 본 논문에서 제시한 역산기법이 가장 정확한 결과를 도출하였으며, 다운홀 시험과 SCPT을 자동화하는데 효율적으로 적용이 될 수 있음을 입증하였다.

• 지구물리와물리탐사
• /
• 제12권3호
• /
• pp.255-262
• /
• 2009

가스하이드레이트 개발 시 중요한 요소 중의 하나인 가스하이드레이트 부존층의 정확한 심도를 얻기 위해서는 일반적으로 고해상도 정보를 주는 VSP 탐사가 많이 수행되어왔다. 본 연구에서는 가스하이드레이트 부존지역인 울릉분지 지역의 UBGH09 시추공에서 국내 최초로 삼성분 지오폰을 사용하여 획득한 제로오프셋 VSP 탐사자료로부터 구간속도를 추출하고 corridor 중합단면도를 생성하여 지표탄성파탐사 자료와 비교해 보았다. 먼저 Phase rotation을 통하여 삼성분 지오폰으로 기록된 신호들이 실제 수직 및 두 개의 수평성분 자료들이 되도록 변환 시키고, 파동장 분리과정을 통하여 직접파와 반사파를 분리하였다. 제로오프셋 VSP 자료의 직접파로부터 추출한 구간속도를 음파검층의 구간속도와 비교한 결과 두 구간속도 모두 비슷한 양상을 나타냈다. BSR 상부의 하이드레이트 층에서는 빠른 속도를 보이고 BSR 하부에서 속도가 급격히 떨어지는 것으로 보아 BSR 하부에 유리가스가 집적되어 있는 것으로 보인다. 제로오프셋 VSP 반사파 자료로부터 corridor 중합단면도를 생성하여 지표탄성파탐사의 CMP 중합단면도와 비교한 결과 대부분의 반사이벤트들이 잘 부합하였으며, phase rotation을 적용한 결과가 phase rotation을 적용하지 않았을 경우보다 더 잘 부합하였다.또한 transverse 성분을 이용하여 만든 corridor 중합단면도와 지표탄성파 CMP 중합단면도를 비교한 결과 지표탄성파 CMP 중합 단면도 내의 모오드 변환된 PS 반사파들을 확인할 수 있었다.

• 한국석유지질학회지
• /
• 제1권1호
• /
• pp.29-35
• /
• 1993

탄성파 트레이스 자료로부터 지구내부의 속도구조를 추정해내는 과정을 속도역산이라고 한다. 그러나 탄성파 트레이스자료가 담고 있는 정보의 주파수 대역은 제한되어 있으므로, 이것으로부터 광대역의 반사계수 수열을 구해내는 것은 쉽지가 않다. 본 연구에서는 선형 프로그래밍 처리기법을 이용하여 탄성파 트레이스의 진폭과 속도분석 자료를 적절히 결합함으로써 각 자료와 동일한 주파수 성분을 갖는 가장 단순한 형태의 속도구조를 구하고, 이것을 의속도 단면도로 도시하였다. 먼저 가상으로 설정한 모델로부터 프로그램의 효용성을 검증해 보았고 이를 토대로 우리나라 황해 제 2광구의 반사법 탄성파 자료에 적용하였다. 그리고 의속도 단면도와 이 지역에 위치한 시추공의 해석결과를 연관시켜 분석함으로써 백악기 화산암과 백악기 변성화산암의 경계를 구분하고 암상을 결정할 수 있었다.

• 한국지반공학회지:지반
• /
• 제10권4호
• /
• pp.29-38
• /
• 1994

공내하향 탄성파시험(downhole seismic tests) 결과를 해석하는 새로운 방법을 소개하였다. 이 방법은 역산이론(inverse thery)에 근거 하며, 직접산출법(direct measurements)과 간접산 출법(indirect measurements)과 같은 기존 방법으로는 불가능한 정도의 자세한 탄성파 속도 프로파일(profile)을 얻을 수 있다. 더욱이, 이 방법은 속도 프로파일의 신뢰도를 평가 할 수 있는 장점이 있다. 지반구조가 서로 다른 5개소의 현장에서 이 방법의 적용성을 검증하였다.