The 2017 Pohang earthquake afflicted more significant economic losses than the 2016 Gyeongju earthquake, even if these earthquakes had a similar moment magnitude. This phenomenon could be due to local site conditions that amplify ground motions. Local site effects could be estimated from methods using the horizontal-to-vertical spectral ratio, standard spectral ratio, and the generalized inversion technique. Since the generalized inversion method could estimate the site effect effectively, this study modeled the site effects in the Korean peninsula using the generalized inversion technique and the Fourier amplitude spectrum of ground motions. To validate the method, the site effects estimated for seismic stations were tested using recorded ground motions, and a ground motion prediction equation was developed without considering site effects.

전역 최적화 문제의 해를 유전 알고리즘을 사용하여 얻어 완전파형역산을 수행하고 층상 반무한체의 물성치를 추정하는 기법을 제안한다. 조화 수직 하중이 작용하는 층상 반무한체의 동적 응답을 측정하고, 이를 추정 물성치를 사용하여 계산된 응답과 비교한다. 응답의 추정치는 mid-point integrated finite element와 perfectly matched discrete layer를 사용하여 구성된 thin-layer model로부터 얻는다. 전역 최적화 문제의 목적 함수는 응답의 관측치와 추정치의 차이에 대한 L2-norm으로 계산된다. 유전 알고리즘을 사용하여 전역 최적화 문제의 해를 구하여 완전파형역산을 수행한다. 제안된 기법을 기본 진동 모드 뿐만이 아니라 고차 진동 모드도 우세한 다양한 층상 반무한 매질에 적용하여, 측정치가 잡음을 포함하지 않는 경우와 포함하는 경우 모두에 대해서 제안된 완전파형역산 기법은 층상 반무한체의 재료 특성을 추정하는데 적합함을 확인할 수 있다.

PURPOSES: The objective of this study is to evaluate the structural capacity of asphalt pavement in subsurface cavity sections using falling weight deflectometer (FWD) backcalculation method. METHODS: It is necessary to analyze the reduction of structural capacity in asphalt pavements due to the occurrence of subsurface cavities. The FWD testing was conducted on the cavity and intact asphalt pavement in the city of Seoul. The GAPAVE, backcalculation program for FWD deflections, was utilized to determine the layer moduli in asphalt pavements. The remaining life of asphalt pavements in cavity sections were predicted using the pavement performance model for fatigue cracking. The backcalculated layer moduli between cavity and intact sections were compared to determine the reduction of structural capacity due to subsurface cavity. The relationship between the reduction of layer modulus and cavity depth/length was analyzed to estimate the effect of cavity characteristics on the structural capacity degradation. RESULTS: According to the FWD backcalculation results, the modulus of asphalt layer, subbase, and subgrade in cavity sections are generally lower than those in intact sections. In the case of asphalt layers, the backcalculated modulus in cavity section was reduced by 50% compared to intact section. A study for the prediction of remaining life of cavity section shows that the occurrence of subsurface cavity induces the decrease of the pavement life significantly. It is found that there is no close relationship between the backcalculated modulus and cavity length. However, the reduction of asphalt layer modulus is highly correlated with the cavity depth and was found to increase with the decrease of cavity depth. CONCLUSIONS : This reduction of structural capacity due to the occurrence of cavities underneath asphalt pavements was determined using FWD backcalculation analysis. In the future, this approach will be utilized to establish the criteria of road collapse risk and predict the remaining life of cavity sections under numerous varied conditions.

PURPOSES: The objective of this study is to analyze the relationship between the FWD back-calculated modulus and dynamic modulus of asphalt layers for existing asphalt pavements. METHODS: To evaluate the dynamic modulus of the asphalt mixture in the existing and new asphalt layers, the uniaxial direct tension test was conducted on small asphalt specimens obtained from the existing asphalt-covered pavements. A dynamic modulus master curve was estimated by using the uniaxial direct tension test for each asphalt layer. The falling weight deflectometer (FWD) testing was conducted on the test sections, and the modulus values of pavement layers were back-calculated using the genetic algorithm and the finite element method based back-calculation program. The relationship between measured and back-calculated asphalt layer moduli was examined in this study. The normalized dynamic modulus was adopted to predict the stiffness characteristics of asphalt layers more accurately. RESULTS: From this study, we can conclude that there is no close relationship between dynamic modulus of first layer and back-calculated asphalt modulus. The dynamic moduli of second and third asphalt layers have some relation with asphalt stiffness. Test results also showed that the normalized dynamic modulus of the asphalt mixture is closely related to the FWD back-calculated modulus with 0.73 of R square value. CONCLUSIONS: The back-calculated modulus of asphalt layer can be used as an indicator of the stiffness characteristics of asphalt layers in the asphalt-covered pavements.

최근 도로포장 건설은 신설보다는 현재 공용중인 도로포장을 확장하거나 유지보수를 하는 추세이다. 이에 따라 유지보수공법을 적용한 도로포장체에 대한 정확한 지지력 및 공용성 평가가 향후 도로포장의 수명 예측 및 유지보수공법 적용에 매우 중요한 요소가 되고 있다. 본 논문의 목적은 주기적인 유지보수 에 따른 도로포장의 층별 구조적 지지력을 정확하게 평가하기 위한 소형시편 동탄성계수 시험결과와 FWD 비파괴 시험 결과와의 상관관계를 규명하고자 한다. 현장 아스팔트 혼합물의 역학적 물성을 측정하기 위하여 소형시편을 이용한 일축직접인장 시험법을 이 용하여 아스팔트 층별 동탄성계수를 측정하였다. 노스캐롤라이나 주의 6개구간에서 구간별 3~5개 층, 2 개 지점에서 코어하여 소형시편을 채취하였다<그림 1>. 채취한 시편들은 5, 20, 40도에서 25, 10, 5, 1, 0.5, 0.1 Hz 하중조건에서 동탄성계수 시험을 수행하여 동탄성계수 마스터커브를 결정하였다. 본 연구에서는 도로포장층 탄성계수를 역산하기 위하여 유한요소 해석 아스팔트 포장 구조해석 프로그 램을 엔진으로 유전자 알고리즘을 활용한 역산 프로그램을 이용하였다. 본 프로그램은 한국건설기술연구 원에서 개발된 프로그램으로 노스캐롤라이나의 현장 및 실내시험 결과에 적용하였다. 도로포장층 FWD 역산 탄성계수는 아스팔트층의 경우, 탄성계수 범위가 1.5~6GPa이며, 보조기층 탄성계수 범위는 150~850MPa이며 노상층 탄성계수 범위는 60~400MPa이다. 도로포장체의 구조적 지지력은 아스팔트층 의 탄성계수와 두께에 의하여 결정된다. 실내시험을 통하여 결정된 동탄성계수에 두께를 곱한 후 합산한 값을 전체 아스팔트층 두께로 나눈 탄성계수 (Normalized Dynamic Modulus) 산정한 후 FWD 역산탄성 계수와 비교하였다. 본 탄성계수를 사용할 경우 도로포장의 두께의 영향을 고려할 수 있으며 두께가 지지 력 산정에 반영될 수 있다. <그림 2>에 나타난 바와 같이 역산과 실측 탄성계수간의 R square값이 0.73 으로 상관도가 매우 높음을 알 수 있다. 향후에는 FWD 하중 적용시 발생하는 아스팔트 각 층의 주파수 를 산정하여 주파수별 동탄성계수와 FWD 역산 탄성계수간의 상관관계를 검토하고자 한다.

도로포장의 구조적 상태를 평가하기 위하여 FWD(Falling weight deflectometer) 시험이 널리 사용되고 있다. FWD 시험은 자유낙하 하는 추에 의하여 포장표면에 발생하는 처짐량을 측정하는 것으로, 이 표 면처짐량으로부터 포장체 각 층의 탄성계수를 역해석하여 추정한다. 역해석시 포장체의 두께는 코어채취를 통하여 측정한 실측값을 사용하거나 설계값을 사용한다. 실제 도로에서 포장체의 두께는 시공상태, 보수이력 등 다양한 원인에 따라 설계 두께와 상이하고, 두께 측정을 위해 코어를 매번 채취하는 것도 현실적으로 거의 불가능하다. 따라서 역해석시 정확한 포장층 두께를 고려하기가 어려운 실정이다. 이러한 부정확한 포장층 두께로 인하여 역산된 탄성계수의 신뢰성이 저하되고, 포장체 구조해석 및 공용수명 예측에 큰 오류가 발생할 수 있다. 포장체의 층두께를 현장에서 빠르고 정확하게 측정하기 위하여 지표투과레이더 (Ground penetrating radar, GPR)가 사용되고 있다. GPR 시험은 지표면으로 송출된 전자기파가 반사, 회절, 산란된 후 돌아오는 시간과 형상을 기초하여 층구조 및 매설물의 위치 등을 탐지할 수 있는 비파괴 시험이다. 본 연구에서는 최근 건설기술연구원의 일반국도 포장관리시스템(Pavement management system, PMS)에 도입된 접촉식(Ground-coupled) GPR을 이용하여 아스팔트 포장의 층두께를 측정하는 방법을 소개하고, 정확한 층두께가 역산된 포장체의 탄성계수의 신뢰성에 미치는 영향을 분석하였다. 조사의 효율성과 공간적 데이터의 신뢰성을 향상시키기 위하여 접촉식 GPR 안테나를 기존 FWD 장비 하부에 설치하였다. 그림 1은 GPR 시험으로 측정한 데이터로 한 지점에서의 시간이력 데이터(A-scan)와 거리 에 따라 연속적으로 나타낸 데이터(B-scan)를 나타낸다. 그림 1b에 나타난 바와 같이 포장층의 두께가 거리에 따라 달라지는 것을 알 수 있다. 이러한 실제 층두께를 역해석에 고려하여 추정된 탄성계수와 가정된 층두께를 고려한 탄성계수를 비교하여 신뢰성 분석을 하였다.

국내 비행장은 민간전용 비행장 7개와 공군전용 비행장 6개, 민과 군이 공동으로 운영하고 있는 8개 등 약 20개의 주요 비행장(소규모 비행장은 제외)이 운영되고 있다. 그리고 대한민국공군과 한국공항공사에서는 정기적으로 포장평가를 실시하여 포장의 잔존수명 및 항공기가 이착륙할 수 있는 허용하중을 결정하는 포장분류등급(PCN)을 결정하고 있으며, 이는 안전한 비행장 포장유지관리에 중요한 부분이다. 이러한 비행장 콘크리트포장시 탄성계수는 코아채취 후 실내실험을 통한 정적탄성계수를 적용함이 가장 정확한 포장물성치를 적용하는 방법이지만 코아채취 및 실내실험(공항별 약300여개소)으로 포장의 손상 및 시간적, 경제적으로 비효율적이다. 따라서 현재 세계적인 추세인 비파괴조사장비(HWD를 이용하여 처짐량을 측정한 후 역산프로그램을 이용하여 동적탄성계수를 추정하고 결과값을 포장평가 분석에 적용하고 있다. 이러한 동적탄성계수 추정방법은 신속하게 시행될 수 있으나 완전한 선형적인 탄성으로 거동하지 않아 정 적탄성계수와 차이가 발생하게 되며, 이는 결국 포장잔존수명 및 항공기허용하중 분석결과에 영향을 미치게 된다. 따라서 7개 비행장 73개소 슬래브중앙에서 HWD 조사시험을 실시하고 동일 지점에서 코아를 채취하여 동적탄성계수를 역산프로그램을 이용하여 추정하고 정적탄성계수 시험을 실시하였다. 그리고 비행장 콘크리트포장 동적탄성계수 추정에 일반적으로 사용되고 있는 역산프로그램들의 결과 값을 비교분석하였고 상관관계 분석을 통해 보정계수 산출연구를 실시하였다. 현장실험을 통하여 분석한 결과, 최근에 미국 FAA에서 개발된 BAKFAA, 미공병단에서 개발된 WESDEF, 공군에서 주로사용하고 있는 ELMOD, Westergard 식 중에서 BAKFAA프로그램과 Westergard 식이 WESDEF과 ELMOD에 비해 정적탄성계수에 근접하게 탄성계수가 추정되었으며 WESDEF과 ELMOD 결과값은 비슷한 경향을 나타내고 상대적으로 높게 추정되었다. 상관관계분석 결과도 BAKFAA프로그램 과 Westergard 식이 WESDEF과 ELMOD에 비해 높게 분석되었다.

본 연구에서는 석모도 지역의 지열에너지 개발에 중요한 심부 파쇄대 시스템에 대한 보다 광범위하고 정밀한 조사를 위하여 2005-2006년에 수행되었던 44 측점의 탐사 영역을 확대하여 추가로 36 측점에서 자기지전류 및 가청주파수 자기지전류 심부물리탐사 현장 탐사 자료를 획득하였다. 이전의 자료와 함께 모든 자료는 격자상으로 분포하도록 설계하여 관심 영역 내에서 2차원 및 3차원 역산을 통하여 지질 구조 해석을 실시하였다. 역산 해석 결과 지열수가 발견된 시추공 주변으로 파쇄대로 추정되는 최소 2 조의 큰 전기비저항 불연속면을 영상화하였는데, NNE-SWW(북북동-남서) 방향으로 발달한 불연속대는 지열수가 발견된 시추공 방향으로 발달되어있고, 이와 함께 인접한 바다로부터 EW(동-서) 방향으로 불연속면이 발달되어 있음을 볼 수 있었다.

제주도는 한반도 남부에 위치한 화산섬으로서, 화산활동 및 한반도와의 지체구조적 연결성과 관련된 중요성 때문에 오랫동안 지질학적/지구물리학적 연구의 대상이었다. 최근 제주도에서 심부지열자원 조사를 위한 다수의 저주파수 MT 자료가 획득됨에 따라 제주도 심부구조에 대한 접근이 가능해졌다. 그러나 제주도와 같은 섬환경에서는 주변 해양의 영향으로 저주파수 MT 자료가 왜곡되며, 이는 섬의 심부구조 규명을 어렵게 만드는 요인이다. 이 연구에서는 기존에 획득된 총 102 측점의 MT 자료에 대해 해양효과 보정을 실시하였다. 해양효과 보정은 1 Hz 이하에서 XY 모드와 YX 모드 전기비저항과 위상을 매우 유사하게 만들며, 전도성 하부지각의 존재를 명확하게 보여주었다. 제주도 남쪽 및 북쪽 측선의 보정된 MT 자료에 대해 2차원 역산을 수행한 결과, 평균적으로 깊이 20 km에서 고비저항 상부지각과 저비저항 하부지각의 경계를 명확히 확인할 수 있었다.

본 연구에서는 최소모델영역 연산자를 MT(magnetotelluric) 2차원 역산 알고리듬에 적용하여 역산 해의 대비를 향상시키고자 하였다. 이를 위하여 creeping법에 기초한 최소자승 역산에 최소모델영역 연산자를 수치적으로 유도하여 알고리듬을 구현하였으며, 공간함수로서의 평활화 상수를 도입한 ACB (Active Constraint Balancing) 법을 동시에 적용하여 최소모델영역 연산자를 이용할 때 단점으로 지적되었던 역산 해의 안정성을 향상시켰다. 고립된 단일 이상체 모델에 대한 수치실험을 통하여 MT 역산에 있어서 최소영역 연산자의 효과를 기존의 2차 미분연산자와 비교 분석하여 MT 역산에서의 특징을 고찰하였다. 또한 다중 이상체 모델에 대한 실험을 통하여 Occam 역산과 비교하여 최소모델영역 연산자를 이용한 역산 해의 특징을 비교 분석하였으며 현장 자료에의 적용을 통하여 그 적용성을 살펴보았다.

최근 들어 국내 광산개발에 대한 관심이 고조되면서 새로운 광상탐사 및 기존 광상의 연장성 확인 등의 목적으로 물리탐사를 적용하는 사례가 늘고 있다. 금속광상 탐사의 경우, 광화대가 주변에 비해 높은 전기전도도를 갖기 때문에 이를 탐지하기 위하여 전기비저항 탐사를 수행하는 것이 일반적이며, 현실 여건상 대부분 2차원 탐사가 수행된다. 그러나 국내 금속광상의 대부분이 맥상으로 분포하고 있으며, 광맥의 폭이 변한다거나 광맥이 단층에 의해 끊어져 있는 등 매우 복잡한 3차원 구조를 띤다. 따라서 3차원 광체구조에 대하여 2차원 탐사 및 2차원 해석을 수행할 경우 왜곡된 해석을 초래할 수 있다. 이에 이 연구에서는 이러한 3차원 광체구조에 대한 3차원 전기비저항 탐사의 적용성을 검토하기 위해 주향에 수직한 여러 측선에 대하여 2차원 쌍극자-쌍극자 탐사를 수행하여 얻은 자료와 3차원 단극자-단극자 배열을 이용하여 얻은 자료들을 각각 3차원 역산한 후 이들을 2차원 해석결과와 비교하였다. 3차원 맥상광체로는 맥폭이 변하는 모형과 맥이 단층에 의해 어긋난 모형 등을 가정하였다. 2차원 쌍극자-쌍극자 배열자료에 대해 3차원 역산을 수행하여 얻은 결과를 3차원 단극자-단극자 배열자료와 비교하면 단극자-단극자 탐사자료가 전체적인 구조는 잘 보여주나 배열의 특성상 쌍극자-쌍극자 탐사자료만큼 정확한 결과를 제시하지 못함을 확인할 수 있었다. 따라서 실제 탐사시 주향의 방향에 대한 정보를 알 수 있다면 주향에 수직한 2차원 측선들에 대해 쌍극자-쌍극자 탐사를 수행하고 이를 3차원 해석하는 것이 가장 바람직하다. 그러나 지표면에 광체가 드러나 있는 경우에 대해서는 이미 개발이 완료된 상태이며, 현재 남아있는 광상은 지하 깊은 곳에 매몰되어 있는 경우가 대부분이다. 이러한 경우 주향에 대한 정보를 알기 어려우므로 비록 해상도는 떨어지더라도 3차원 단극자-단극자 탐사를 수행하고 해석하는 것이 바람직할 것으로 보인다.

By Multiple-measurement method using Inverse Heat Conduction Method and parameter estimation method, thermo physical properties were measured, and its estimation accuracies were compared with other reference data. Values of the thermophysical properties determined from the measurement method are in the range of data reported in the literature. From the above, it should be pointed out that the inverse algorithm and multiple measurement method is particularly simple to code and can be run on the personal computer used for the data acquisition. Especially, the inverse algorithm for a heat flux determination can be practically used in every case where the surface measurement are not accessible and where it is impossible to correctly place some sensors with the body. It is found that the measurement method is more practical, convenient and time-saving than other steady-state methods.

최근에 수행되고 있는 중력 탐사는 고분해능의 중력계와 GPS(Global positioning system)를 통한 정밀한 측정과 측지가 이루어지고 있다. 중력탐사에서 모델링과 역산의 기술은 많은 발전이 있어왔지만, 중력자료처리는 거의 변화가 없었다. 통상적인 정밀한 중력 자료 보정을 통한 부우게이상은 측정점의 고도에서 기준면까지의 물질의 영향을 일정한 밀도를 이용해 제거해 버리기 때문에 측정점 바로 하부의 이상체에 의한 영향을 상당히 왜곡시키게 된다. 본 연구에서는 탐사 지역의 지형을 DEM(Digital Elevation Map) 자료와 Multiquadric equation을 이용하여 실제 지형과 유사한 Multiquadric surface를 자동적으로 구성하고, 이를 블록화 함으로써 보정의 대상이었던 기준면 상부에 대한 밀도를 탐사 지역의 지질 정보와 지형을 포함하는 역산을 통해 수치적으로 계산하였다. 이러한 지형을 포함한 역산 방법을 3차 원중력지형역산(3DGTI; 3-D Gravity Terrain Inversion)이라 한다. 이 연구의 효율성을 검증하기 위하여 주변암과 밀도차가 존재하는 관입지역에 대한 모델을 구성하고 적용한 결과 기존의 부게 보정 방법을 적용한 부게 이상도에 비해 자료의 왜곡이 감소하는 효과를 얻을 수 있었다. 이를 통하여 지형 역산을 통한 객관적인 부게 밀도의 결정과 부게 보정시 실제의 수평적인 밀도 변화를 반영함으로써 기존의 문제점을 보완하였다. 게다가, 3DGTI로부터 얻어진 밀도분포는 지형의 윤곽을 그대로 표현하고 있어서 보다 실질적인 지질을 보여준다고 하겠다. 이 방법을 화강암체가 관입하고 있는 마산·창원 일대에서의 중력 탐사 자료에 적용해본 결과 기존 방법보다 관입 화강암체의 위치와 그 규모를 알아내는데 더 효과적이었다. 따라서, 수평적인 밀도 변화가 뚜렷하게 존재하는 지역의 경우, 새로운 중력 자료 처리 방법이 기존의 부게 보정에서 발생하였던 문제점을 해결함으로써 천부의 분해능을 높이고, 심부의 밀도 분포도 좀더 정확하게 계산할 수 있으리라 생각된다.

토양이나 암반의 물성을 조사하기 위하여 시추공조사가 흔히 이루어진다. 그러나 시추조사의 결과는 불연속적이고 시추공과 시추공 사이의 물성은 두 시추공의 조사결과를 내삽하여 구할 수 밖 에 없다. 그러나 이러한 내삽법을 이용한 해석은 지반의 수평적 변화가 심하지 않은 경우에만 가능하다. 연약지반의 연속적인 2차원 S파 속도구조를 구하기 위하여 표면파 역산 방법을 사용하였다. 역산 결과를 해석하기 위하여 역산 결과의 해상도를 역산 결과와함께 제시하였다.

남한 중남부에 위치하는 한국교원대학교 단주기 지진 관측망에 1997∼1998년 사이 포착된 원거리 지진들의 종파 초동 시간의 상대적 도달시간차를 분석하고 향후 있을 한반도 전체규모의 종파속도구조대비의 연구를 위하여 예비적인 토모그라피(tomography) 역산을 실시하였다. 각 지진 관측점에서 측정된 종파의 상대적 도달시간차들은 최대 0.7초로 측정되었다. 경기육괴와 옥천대의 경계 근처의 관측점들에서 계산된 종파의 상대적 도달시간차들은 역방위각의 변화에 따라 비교적 큰 값의 변화를 나타내며 그 지역 지각 하부에 속도구조 이상체가 존재할 수 있음을 지시한다. 조사지역 상부 맨틀이 수평적으로 균질하다는 가정하에 지각 전체에 대한 토모그라피 역산 결과는 조사지역 북서쪽에 북동-남서 방향으로 상대적으로 낮은 종파 속도를 가지는 구조가 나타난다. 역산된 예비적 성격의 상대적 속도 구조이상체는 경기육괴와 옥천대의 경계를 따라 대비되는 듯한 경향을 나타낸다. 이번 연구에 사용된 블록의 크기나 위치 등에대한 입력 파라메터는 차후 자료가 충분히 집적될 경우 정밀분석을 위한 입력 파라메터로서 이용될 것이다.

금강하구 부근의 연약지반에서 획득한 탄성파 자료(25개 shot gather)를 역산하여 2차원 S파 속도구조를 구하였다. 탐사측선위에 위치한 2개 시추홀에서 지질조사를 실시하고 표준관입시험을 실시하였다. 2차원 S파 속도구조는 대상지역의 지층이 두께 1∼3m의 상부층(S파 속도 200∼700m/sec), 두께 5∼8m의 중간 저속도층(S파 속도 100m/sec∼400m/sec)과 그 아래 S파 속도 1000m/sec 이상의 하부층으로 이루어져 있음을 보인다. 저속도층은 탐사측선의 한쪽 끝에서 다른 쪽 끝으로 가면서 그 두께가 얇아지고, 기반암의 깊이도 얕아진다. S파 속도구조와 지층의 지질, 표준관입시험 값을 검토한 결과, 저속도층은 clay층과 밀접한 관련이 있는 것으로 사료된다. 이에 비해 Standard Penetrarion Test 값은 지층의 성분과는 연관성을 보이지 않고, 깊이에 따라 증가하는 것으로 나타난다. 이 연구는 표면파 역산이 연약지반의 S파 속도구조를 밝히는 데 효과적으로 사용될 수 있음을 보여준다. 또한 표면파 역산 방법은 연약지반에 흔히 존재하는 지표에 가까운 지하수면, 또는 저속도 층으로 인한 굴절파 탐사방법의 한계를 극복할 수 있는 방법을 제공한다.

In most marine seismic refraction studies using shots detonated near the sea surface, the location of receivers and shots are usually assumed to be located at the point of deployment and the shot times are accurately known. In the case of bottom shots and receivers, however, the above assumption can be led to big errors in the shot and receiver positions. The errors in the shot and receiver positions can cause severe distort in the interpretation of deep ocean refraction studies. In this study, an inversion algorithm has been developed to reduce the errors in the shot and receiver position. The inversion method has been used to obtain the shot and receiver positions and the shot origin times of a deep ocean refraction seismic experiment which was performed near the Juan de Fuca ridge in 1985. The comparison of two seismic record sections, one using the initial shot and receiver positions and the other using the shot and receiver positions obtained by the inversion algorithm, shows that the seismic signals emerges much more gradually in the latter. This indicated that the posotions of shots and receivers should be relocated using some inversion methods for better interpretation of deep marine seismic refraction studies.