• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.441-448
• /
• 2003

Refraction tomography requires an algorithm for efficiently computing the traveltimes and their $Fr\'{e}chet$ derivatives. We have attempted to solve the damped wave field using the frequency domain finite element model ing and then invoked the reciprocity theorem to calculate the $Fr\'{e}chet$ derivative of the traveltime with respect to the subsurface parameter. Then, we used a damped least square method to invert the traveltimes of the Marmousi 2 model. Numerical tests demonstrate that the refraction tomography with large aperture data can be used to estimate the smooth velocity model for the prestack depth migration.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.402-409
• /
• 2003

Inversion of traveltime requires an efficient algorithm for computing the traveltime as well as its $Frech\hat{e}t$ derivative. We compute the traveltime of the head waves using the damped wave solution in the Laplace domain and then present a new algorithm for calculating the $Frech\hat{e}t$ derivative of the head wave traveltimes by exploiting the numerical structure of the finite element method, the modem sparse matrix technology, and SWEET algorithm developed recently. Then, we use a properly regularized steepest descent method to invert the traveltime of the Marmousi-2 model. Through our numerical tests, we will demonstrate that the refraction tomography with large aperture data can be used to construct the initial velocity model for the prestack depth migration.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2000년도 가을 학술발표회 논문집
• /
• pp.577-584
• /
• 2000

Seismic refraction survey was peformed for 10 lines along NE-SW and NW-SE directions above Nampoong gallery at Makyo-ri, Dogye, Samcheok, Kangwon-do. 48 geophones were laid in line with the interval of 1m, and a 5Kg hammer was used as a source at 5 points for each line. Data processing was done using reciprocal time method, GRM, and traveltime tomography which utilizes wavefront expansion method for forward process and SIRT for inversion. The result shows that the first layer has its lower boundary between 3.49m and 8.88m. The P-wave velocity of the first and the second layer were estimated as 270 360m/s and 1550 1940m/s respectively. When the boundary of the first and second layer is smooth enough and the velocity difference is large enough, GRM has little advantage over reciprocal time method. The result of reciprocal method and traveltime tomography shows consistency. The northeast part of the boundary has syncline structure, which is similar to the topography above. This implies that the collapse of the cavities of Nampoong gallery result in the subsidence of the ground surface. The subsidence is in progress across the Youngdong railroad, therefore a proper reinforcement work is required.

• 지구물리와물리탐사
• /
• 제20권3호
• /
• pp.163-175
• /
• 2017

육상 탄성파 탐사에서 불규칙한 지표 고도차와 천부 지층의 풍화대는 자료의 반사파 신호를 왜곡시킨다. 그러므로 일반적인 육상 탄성파 탐사 자료는 정적 보정(static correction)을 통해 이러한 왜곡을 보정해야 한다. 정적 보정을 하기 위해 천부 지층의 속도가 필요하며, 이 속도는 굴절법 탐사를 통해 구할 수 있다. 그러나 육상 탄성파 탐사자료는 탐사 현장에 대한 허가나 현장 여건 그리고 장비 보유 현황에 따라 제한된 형태로 취득되는 경우가 많다. 이러한 상황에서는 송수신기 배열이 제한되어 굴절법 탐사에 적합한 자료를 얻을 수 없기 때문에 반사법 자료로 굴절법 해석을 수행할 수 밖에 없다. 따라서 본 연구에서는 불규칙 지표 모델에서 제한된 송수신기 배열을 이용하여 얻은 반사법 자료로 굴절 주시 역산을 수행하고 결과의 신뢰성을 분석하였다. 불규칙한 지표와 평탄한 지표를 갖는 모델에 대한 역산 결과를 비교하여 지표 고도 차이가 역산 결과의 신뢰성에 약간의 영향을 주는 것을 확인하였다. 또한, 송신원 개수에 따른 역산을 통해 극히 적은 수의 송신원이 아니라면 역산 결과의 신뢰성에 크게 영향을 주지 않음을 보였다. 반면에 제한된 배열의 수신기를 사용할 경우에는 수신기 전개영역의 중첩된 부분에서 속도왜곡이 발생함을 관찰하였고, 속도왜곡이 발생하지 않기 위한 최소 중첩영역의 크기를 제안하였다. 마지막으로 현장 자료의 지형 조건과 송수신기 배열 정보를 이용하여 수치 모형 실험을 수행하였고, 신뢰성 분석 결과를 검증하였다. 신뢰성 분석 결과를 토대로 현장 자료 역산 결과에 신뢰할 수 있는 영역과 해석에 유의해야 할 영역을 구분하였다.

• 지구물리와물리탐사
• /
• 제2권1호
• /
• pp.17-25
• /
• 1999

고해상의 지하 전기전도도 영상을 요하는 환경, 공학적인 적용을 위하여, 주파수영 역 전자탐사 자료를 이용한 주시 토모그래피를 수행하였다. 이를 위하여 우선 변환된 파동장을 파선급수(ray series)의 합으로 근사하여 자기장으로부터 직접 파동장의 초동을 구해낼 수 있는 방법을 제시하고 그 정확성 및 적용성을 검토하였다. 균질한 무한공간에서의 자기장을 이용한 발췌결과, 잡음이 없는 자료의 경우 주파수 2개를 사용하여도 아주 정확한 발췌가 이루어져 그 타당성이 입증되었으며, 기존의 파동장으로 직접변환하는 방법과 비교한 결과 더 적은 주파수 자료를 이용하여도 더 정확한 발췌가 이루어졌다. 층서구조 및 경사진 파쇄대 구조에 대하여 초동발췌 및 반복적 비선형 토모그래피를 적용하여 만족할만한 영상을 얻었다. 그러나 전기전도도의 비율이 큰 경우는 설정한 가정에 부합되지 않아 토모그래피 영상에서 전기비저항이 낮은 층이 확대되어 나타나는 결과를 보였다. 초동발췌를 위한 시간은 하나의 송, 수신 배열에 대하여 IBM PC로 10초 내외로, 현장에서 탐사를 수행하는 도중 실시간으로 발췌가 가능하다. 또한 토모그래피 역산을 위한 시간도 약 3분 이내로 전체 송, 수신 배열에 대한 측정이 끝남과 동시에 지하의 단면 영상을 확인할 수 있어 전자탐사 토모그래피의 현장 적용성을 한층 높일 수 있을 것으로 기대된다.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2000년도 정기총회 및 특별심포지움
• /
• pp.132-146
• /
• 2000

강원도 삼척시 도계읍 마교리의 남풍갱 상부 농경지에서 북동-남서 및 북서-남동 방향으로 총 10개의 측선에 대한 탄성파 굴절법 탐사를 실시하였다. 매 측선마다 지오폰은 일렬로 1m 간격으로 48 개를 설치하였으며 파원은 5Kg해머를 사용하여 5개의 위치에 타격을 가하였다. 굴절법 탐사의 자료처리는 역행주시법 , GRM과 함께 파면확장법을 통한 초동계산과 SIRT를 이용한 역산을 하는 주시 토모그래피를 사용하여 행해 졌다. 계산 결과 조사지역 의 상부 매립 및 퇴적층의 하부 경계면은 3.49m에서 8.88m의 심도로 분포하고 있으며 P파의 속도는 270${\~}$360m/s를 나타내었다. 하부 파쇄암반의 P파 속도는 1550${\~}$1940m/s의 분포를 보였다. 자료처리 결과 이처럼 상부와 하부층의 탄성파 속도 차이가 크게 나타나고 경계면의 굴곡이 완만할 경우에는 GRM이 역행주시법에 대해 갖는 이점이 거의 없음을 발견하였다. 역행주시법과 주시 토모그래피의 결과는 서로 잘 일치하였으며, 조사지역의 북동 방향으로는 상하부층의 경계면이 지표면이 겪은 변화와 동일한 굴곡을 보이고 있다. 이는 남풍갱 폐탄광 지역의 지하 채굴적에 의한 지반이완이 넓은 지역에 걸쳐 상부로 전이되어 나타난trough형 지반침하의 전형적인 양상으로 판단된다.

• 지구물리
• /
• 제9권4호
• /
• pp.343-349
• /
• 2006

소양댐 상류의 미고결 퇴적물 분포를 파악하기 위하여, 길이 94 m의 5개 측선을 따라 2 m 간격으로 204 ms동안 탄성파 자료를 기록하였다. 기록된 자료를 초동주시 토모그래피 역산법과 지연시간법으로 처리하여 조사지의 퇴적물 분포를 종합적으로 구명하였다. 기반암은 평균 14 m 깊이로 비교적 평탄하게 분포하며, 산사면이 있는 남동쪽으로 상승한 형태를 보인다. 이 기반암을 댐 건설 이전에 형성된 기존 퇴적물이 덮고 있으며, 주로 호수 생성 이전에 형성된 토양이나 사질 퇴적물로 구성되어 있다. 최상부층인 최신 퇴적물은 호수 생성 이후 유입된 부유 퇴적물 기원의 이질 퇴적층으로, 하류지역보다 두꺼운 평균 1.6 m 두께로 쌓여 있다.

• 지구물리와물리탐사
• /
• 제10권4호
• /
• pp.314-321
• /
• 2007

최근에 들어서 물성이 서로 다른 두 종류의 탐사자료의 복합역산에 대한 연구가 활발하게 이루어지고 있다. 이는 복합역산에 의하여 훨씬 더 정확한 지하구조 영상을 계산할 수 있을 뿐만 아니라 물리탐사 변수가 아닌 다른 물성 분포의 유도가 물리탐사로서 가능해지기 때문이다. 이 연구에서는 (1) cross-gradient로 정의되는 두 지하구조의 유사성의 최대화, (2) 두 물성간의 상관관계의 최대화, (3) 지하 물성 분포에 대한 선험적 정보의 3 종류 제한을 채택한 탄성파 굴절법 주시 토모그래피와 전기비저항 탐사 자료의 복합역산법을 개발하였다. 지표 전기비저항과 탄성파 굴절법 탐사의 수치실험을 통하여, 제안한 복합역산법의 효용성과 각종 제한조건의 효과를 분석하였다. 특히 제한조건을 적절히 이용할 경우, 탄성파 탐사의 저속도층에 의한 숨은 층 문제를 복합탐사 및 역산으로 해결할 수 있음을 알 수 있었다.

• 지구물리
• /
• 제9권1호
• /
• pp.37-45
• /
• 2006

풍암분지 내에 위치한 140 m 깊이 시험시추공 주변에서 지표 굴절법 및 원거리 수직탄성파 자료의 초동주시를 토모그래피방식으로 동시에 역산하여 시추공 주변의 속도구조를 구하였다 . 지진총으로 발생시킨 지진파는 시추공을 중심으로 N20°E 방향과 이에 수직인 측선을 따라 3 m 간격의 42개 지표 측점과 공내 71 m 깊이에서 수신하였다 . 는 시추공으로부터의 수평거리 -19.5∼+19.5 m인 지표 측점에서 5 kg 해머로 발생시켰으며 , 1∼2 m 간격의 깊이 9∼99 m 구간에서 3성분으로 204 ms 동안 기록하였다 . 지연시간 보정 후 , 이 자료를 지진총으로 기록한 지표 기록자료 초동주시와 동시에 SIRT 방식으로 역산하여 속도 토모그램을 작성하였다 . 속도 토모그램은 평균 1.5 m 두께의 속도 500 m/s 표토층이 속도 3,067∼5,717 m/s의 속도를 갖는 기반암 위에 분포하고 있음을 보여 준다 . 기반암 내부에는 깊이 30∼44 m 구간과 71 m 하부에 속도가 5,130 m/s 이상 되는 매우 단단한 암층이 존재하며 , 이는 코아분석 자료와 비교할 때 사암 및 이암층에 협재하는 역암층의 영향으로 판단된다 . 시추공 주변 조사범위 내에서 대규모 파쇄대나 단층 존재의 증거는 발견되지 않는다 .

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

To investigate the shear-wave velocity structures of the Korean peninsula, exploded seismic signals were recorded for 120 s along a 294-km WNW-ESE line and 150 s along a 335-km NNW-SSE line in 2002 and 2004, respectively. First arrival times of shear wave were inverted to derive the velocity tomograms. Initial shear-wave 1-D models were built using the initial P-wave velocity models used by Kim et al. and $V_p/V_s$ ratios of the IASP91 model. The raypaths indicate existence of mid-crust interfaces at the depth of 2-3 km and 16 km. The deepest significant interface corresponding to the Moho discontinuity varies in depth from 32 km to 36 km. The refraction velocity along the interface varies from 4.4 km/s to 4.6 km/s. The velocity tomograms also indicate existence of a low-velocity zone at the depth of 7.8 km under the Okchon fold belt.