스트리머 방식 탐사 자료의 동시 송신원 전파형 역산을 위한 Global correlation 기반 목적함수 최적화 연구 A Study on Optimization of the Global-Correlation-Based Objective Function for the Simultaneous-Source Full Waveform Inversion with Streamer-Type Data원문보기
동시 송신원 전파형 역산 기법은 계산량을 획기적으로 줄여 전파형 역산의 적용성을 높여준다. 그러나 다수의 송신원 모음 자료를 동시에 모델링하여 사용하기 때문에 관측 자료의 수진기 위치가 송신원에 따라 다른 경우, 나머지(residual) 파동장에 불필요한 값을 생성하게 되고 이는 파형역산의 수렴성을 저해하게 된다. 특히, 제한된 벌림 거리(offset)를 갖는 스트리머 방식의 탐사자료는 동시 송신원 기법을 적용하기에 가장 어려운 자료 형태이다. 이러한 문제점을 극복하기 위해 최근에 global correlation에 기반한 목적함수가 제안되었고, 시간영역 전파형 역산에 성공적으로 적용되었다. 그러나 이 기법은 변형된 목적함수를 사용하기 때문에 나머지 파동장이 왜곡되고 경우에 따라 역산 결과에 부정적인 영향을 주기도 한다. 또한, 여러 가지 장점을 갖고 있는 주파수 영역 파형역산에 적용된 사례는 아직 보고된 적이 없다. 본 논문에서는 이러한 나머지 파동장의 왜곡을 최소화하기 위해 global correlation 계산 시 사용하는 자료에 진폭감쇠 기법을 적용한다. 진폭감쇠를 적용한 자료는 global correlation의 특성을 최적화하여 나머지 파동장의 왜곡을 줄이고 파형역산 결과를 향상시킨다. 시간 영역에서 구한 나머지 파동장을 주파수 영역에서 역전파시킴으로써 global correlation기법을 주파수 영역에서 구현한다. 스트리머 방식의 합성 탐사자료를 이용한 예제를 통해 본 논문에서 제안한 기법이 기존의 global correlation 목적함수에 기반한 동시 송신원 전파형 역산보다 향상된 결과를 얻을 수 있음을 보여준다.
동시 송신원 전파형 역산 기법은 계산량을 획기적으로 줄여 전파형 역산의 적용성을 높여준다. 그러나 다수의 송신원 모음 자료를 동시에 모델링하여 사용하기 때문에 관측 자료의 수진기 위치가 송신원에 따라 다른 경우, 나머지(residual) 파동장에 불필요한 값을 생성하게 되고 이는 파형역산의 수렴성을 저해하게 된다. 특히, 제한된 벌림 거리(offset)를 갖는 스트리머 방식의 탐사자료는 동시 송신원 기법을 적용하기에 가장 어려운 자료 형태이다. 이러한 문제점을 극복하기 위해 최근에 global correlation에 기반한 목적함수가 제안되었고, 시간영역 전파형 역산에 성공적으로 적용되었다. 그러나 이 기법은 변형된 목적함수를 사용하기 때문에 나머지 파동장이 왜곡되고 경우에 따라 역산 결과에 부정적인 영향을 주기도 한다. 또한, 여러 가지 장점을 갖고 있는 주파수 영역 파형역산에 적용된 사례는 아직 보고된 적이 없다. 본 논문에서는 이러한 나머지 파동장의 왜곡을 최소화하기 위해 global correlation 계산 시 사용하는 자료에 진폭감쇠 기법을 적용한다. 진폭감쇠를 적용한 자료는 global correlation의 특성을 최적화하여 나머지 파동장의 왜곡을 줄이고 파형역산 결과를 향상시킨다. 시간 영역에서 구한 나머지 파동장을 주파수 영역에서 역전파시킴으로써 global correlation기법을 주파수 영역에서 구현한다. 스트리머 방식의 합성 탐사자료를 이용한 예제를 통해 본 논문에서 제안한 기법이 기존의 global correlation 목적함수에 기반한 동시 송신원 전파형 역산보다 향상된 결과를 얻을 수 있음을 보여준다.
The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field...
The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.
The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.
* AI 자동 식별 결과로 적합하지 않은 문장이 있을 수 있으니, 이용에 유의하시기 바랍니다.
문제 정의
본 연구에서는 나머지 파동장의 왜곡 현상을 최소화하기 위해서 진폭감쇠를 수행한 후 global correlation을 적용하는 동시 송신원 전파형 역산 알고리듬을 제안한다. Global correlation을 적용한 기존의 동시 송신원 역산 기법이 시간 영역에서 수행된 것과 달리 본 논문에서는 주파수 영역에서 수행하고, global correlation의 특성을 최적화하기 위해 나머지 파동장을 시간 영역에서 분석한다.
본 연구에서는 스트리머 자료에 적용 가능한 효율적인 주파수 영역 동시 송신원 전파형 역산 기법을 제안하였다. 알고리듬의 효용성을 확인하기 위해 나머지 파동장을 분석하였고, 이를 통해 l2-노옴을 이용한 역산 기법은 모델링 자료와 스트리머 방식의 관측자료 간의 벌림 거리 차이로 인해 불필요한 이벤트가 발생됨을 확인하였다.
, 2009). 본 연구에서는 스트리머 탐사 자료에 대해 효과적으로 동시 송신원 기법을 적용할 수 있는 알고리듬을 구축하기 위해 나머지 (residual) 파동장을 분석하였다. Fig.
본 연구에서는 시간 영역에서 정의된 global correlation 기법을 주파수 영역 동시 송신원 파형역산에 적용하기 위해 주파수 영역에서 생성된 자료를 시간 영역으로 변환시켜 나머지 파동장을 계산하였다(Fig. 4). 이 때 주파수 영역의 자료를 시간 영역으로 변환 시 발생하는 두루마리 현상(wrap around)을 억제하기 위해 직접파 도달 이전의 신호에 대해 뮤팅(muting)을 적용하였다(Fig.
가설 설정
동시 송신원 전파형 역산 기법을 보다 명확하게 이해하기 위해 송신원이 2개인 경우를 가정하여 이론을 전개하기로 한다. 먼저, 기존의 주파수 영역 l2-노옴(norm) 목적함수(objective function)는 다음과 같이 정의된다.
제안 방법
본 연구에서는 나머지 파동장의 왜곡 현상을 최소화하기 위해서 진폭감쇠를 수행한 후 global correlation을 적용하는 동시 송신원 전파형 역산 알고리듬을 제안한다. Global correlation을 적용한 기존의 동시 송신원 역산 기법이 시간 영역에서 수행된 것과 달리 본 논문에서는 주파수 영역에서 수행하고, global correlation의 특성을 최적화하기 위해 나머지 파동장을 시간 영역에서 분석한다. 역산 알고리듬의 유효성을 검증하기 위해 제안된 동시 송신원 파형역산 결과를 진폭감쇠가 적용되지 않은 기법과 비교한다.
앞 절에서 언급한 바와 같이 시간 영역으로의 역 푸리에 변환 시 관측 자료 및 모델링 자료에 두루마리 현상(wrap around)을 억제하기 위한 전 처리로써 뮤팅을 적용하였다. 또한 global correlation 기법을 최적화하기 위해 진폭감쇠 작업을 수행하였다. 전 처리된 자료에 대해 global correlation 기법을 적용하여 나머지 파동장을 계산하고, 푸리에 변환을 수행한다.
4에 나타내었다. 먼저, 동시 송신원 기법을 적용하여 주파수 영역에서 자료를 생성하고, 역 푸리에(Fourier) 변환을 수행한다. 이 때 여러 개의 송신원을 동시에 가하여 하나로 합쳐진 자료를 역 푸리에 변환하기 때문에 추가적인 계산 시간이 많이 요구되지 않아, 기존의 주파수 영역 동시 송신원 전파형 역산과 비슷한 계산 효율성을 지닌다.
본 연구에서는 global correlation 기법을 적용한 동시 송신원 전파형 역산 기법을 주파수 영역에서 구현하는 방법을 제시하였고, 전체적인 역산 과정을 Fig. 4에 나타내었다. 먼저, 동시 송신원 기법을 적용하여 주파수 영역에서 자료를 생성하고, 역 푸리에(Fourier) 변환을 수행한다.
스트리머 탐사 자료에 대한 동시 송신원 기법의 적용 가능성을 확인하기 위해 제한된 벌림(limited offset) 거리를 갖는 관측 자료를 생성하여 나머지 파동장을 분석하였다. Fig.
이 때 여러 개의 송신원을 동시에 가하여 하나로 합쳐진 자료를 역 푸리에 변환하기 때문에 추가적인 계산 시간이 많이 요구되지 않아, 기존의 주파수 영역 동시 송신원 전파형 역산과 비슷한 계산 효율성을 지닌다. 앞 절에서 언급한 바와 같이 시간 영역으로의 역 푸리에 변환 시 관측 자료 및 모델링 자료에 두루마리 현상(wrap around)을 억제하기 위한 전 처리로써 뮤팅을 적용하였다. 또한 global correlation 기법을 최적화하기 위해 진폭감쇠 작업을 수행하였다.
4). 이 때 주파수 영역의 자료를 시간 영역으로 변환 시 발생하는 두루마리 현상(wrap around)을 억제하기 위해 직접파 도달 이전의 신호에 대해 뮤팅(muting)을 적용하였다(Fig. 3 (e)).
진폭감쇠는 신호에 감쇠함수를 곱해주어 직접파의 에너지를 강조시키는 효과가 있으므로, 직접파에 가중치를 주고 global correlation을 수행하면 효과적으로 탐사 자료 획득 구조에 맞는 나머지 파동장을 얻을 수 있었다. 주파수 영역에서 생성된 자료를 시간 영역으로 변환하여 나머지 파동장을 계산하는 과정에서 발생하는 두루마리 현상(wrap around)을 억제하기 위해 선처리 작업으로써 뮤팅을 적용하였다. 수치 예제를 통해 진폭감쇠와 global correlation 기법에 기반한 본 알고리듬은 기존의 global correlation 기법만을 적용하는 동시 송신원 전파형 역산에 비해 모델의 심부 및 암염 구조에서 향상된 역산 결과를 생성하며, 스트리머 방식의 탐사 자료에 대한 역산에 적합함을 확인하였다.
(2009)이 제안한 주파수 별 가중평균 기법을 적용하여 최종 경사 방향을 얻는다. 최종적으로 얻은 경사 방향에 따라 모델 변수를 반복적으로 개선하여 역산 속도모델을 구한다.
대상 데이터
송신원은 192 m 간격으로 48개를 설정하였고, 각 송신원에 대한 수진기 간격은 16 m이다. 스트리머 탐사 방식의 관측 자료를 생성하기 위해 근거리 벌림(offset)은 176 m, 원거리 벌림은 6.4 km로 설정하였다. 인공합성 관측 자료는 주파수 영역의 유한요소 모델링 기법을 적용하여 생성하였다.
데이터처리
Global correlation을 적용한 기존의 동시 송신원 역산 기법이 시간 영역에서 수행된 것과 달리 본 논문에서는 주파수 영역에서 수행하고, global correlation의 특성을 최적화하기 위해 나머지 파동장을 시간 영역에서 분석한다. 역산 알고리듬의 유효성을 검증하기 위해 제안된 동시 송신원 파형역산 결과를 진폭감쇠가 적용되지 않은 기법과 비교한다.
이론/모형
본 연구에서 제안한 역산 알고리듬의 효용성을 확인하기 위해 Fig. 1 (a)와 같이 상부에 물층을 추가한 Marmousi 모델을 사용하였다(Versteeg, 1994). 물층이 추가된 Marmousi 모델의 크기는 수평 거리가 9.
4. The flowchart of the simultaneous-source full waveform inversion algorithm using our attenuated global correlation technique.
, 2007)시켜 효율적으로 목적함수의 경사방향을 계산한다. 수렴 속도를 향상시키기 위해 유사-헤시안(pseudo-Hessian)행렬의 대각 성분을 이용하여 경사 방향을 선처리(pre-conditioning)하고, Ha et al. (2009)이 제안한 주파수 별 가중평균 기법을 적용하여 최종 경사 방향을 얻는다. 최종적으로 얻은 경사 방향에 따라 모델 변수를 반복적으로 개선하여 역산 속도모델을 구한다.
또한 기존의 global correlation 목적함수에 기반한 동시 송신원 기법은 이러한 이벤트를 제거하는 효과를 보이나, 일부 구간에서 왜곡된 신호가 생성되었다. 왜곡 현상을 최소화하기 위해 진폭감쇠를 적용한 후 global correlation에 기반한 동시 송신원 기법을 수행하였다. 진폭감쇠는 신호에 감쇠함수를 곱해주어 직접파의 에너지를 강조시키는 효과가 있으므로, 직접파에 가중치를 주고 global correlation을 수행하면 효과적으로 탐사 자료 획득 구조에 맞는 나머지 파동장을 얻을 수 있었다.
4 km로 설정하였다. 인공합성 관측 자료는 주파수 영역의 유한요소 모델링 기법을 적용하여 생성하였다. 최대 주파수가 18 Hz인 일차 미분 가우스 함수(first derivative Gaussian function)를 송신원 파형으로 이용하였다.
또한 global correlation 기법을 최적화하기 위해 진폭감쇠 작업을 수행하였다. 전 처리된 자료에 대해 global correlation 기법을 적용하여 나머지 파동장을 계산하고, 푸리에 변환을 수행한다. 주파수 영역으로 변환된 나머지 파동장을 역전파(Shin et al.
인공합성 관측 자료는 주파수 영역의 유한요소 모델링 기법을 적용하여 생성하였다. 최대 주파수가 18 Hz인 일차 미분 가우스 함수(first derivative Gaussian function)를 송신원 파형으로 이용하였다. 역산을 위한 초기 모델은 심도가 깊어짐에 따라 속도가 1.
성능/효과
3 (e))를 얻을 수 있기 때문에 기존의 global correlation을 적용한 동시 송신원 전파형 역산 기법(Fig. 5 (b))보다 모델의 배사 구조 및 암염 구조의 연속성이 향상된 역산 결과를 얻을 수 있었다.
알고리듬의 효용성을 확인하기 위해 나머지 파동장을 분석하였고, 이를 통해 l2-노옴을 이용한 역산 기법은 모델링 자료와 스트리머 방식의 관측자료 간의 벌림 거리 차이로 인해 불필요한 이벤트가 발생됨을 확인하였다. 또한 기존의 global correlation 목적함수에 기반한 동시 송신원 기법은 이러한 이벤트를 제거하는 효과를 보이나, 일부 구간에서 왜곡된 신호가 생성되었다. 왜곡 현상을 최소화하기 위해 진폭감쇠를 적용한 후 global correlation에 기반한 동시 송신원 기법을 수행하였다.
주파수 영역에서 생성된 자료를 시간 영역으로 변환하여 나머지 파동장을 계산하는 과정에서 발생하는 두루마리 현상(wrap around)을 억제하기 위해 선처리 작업으로써 뮤팅을 적용하였다. 수치 예제를 통해 진폭감쇠와 global correlation 기법에 기반한 본 알고리듬은 기존의 global correlation 기법만을 적용하는 동시 송신원 전파형 역산에 비해 모델의 심부 및 암염 구조에서 향상된 역산 결과를 생성하며, 스트리머 방식의 탐사 자료에 대한 역산에 적합함을 확인하였다. 진폭감쇠를 적용한 본 파형역산 기술은 global correlation의 특성을 최적화하여 나머지 파동장의 왜곡을 줄이고 파형역산 결과를 향상시키는 것으로 생각된다.
본 연구에서는 스트리머 자료에 적용 가능한 효율적인 주파수 영역 동시 송신원 전파형 역산 기법을 제안하였다. 알고리듬의 효용성을 확인하기 위해 나머지 파동장을 분석하였고, 이를 통해 l2-노옴을 이용한 역산 기법은 모델링 자료와 스트리머 방식의 관측자료 간의 벌림 거리 차이로 인해 불필요한 이벤트가 발생됨을 확인하였다. 또한 기존의 global correlation 목적함수에 기반한 동시 송신원 기법은 이러한 이벤트를 제거하는 효과를 보이나, 일부 구간에서 왜곡된 신호가 생성되었다.
5 (a)에 나타내었다. 이를 통해 스트리머 자료에 기존의 동시 송신원 역산 기법을 적용하여 역산할 경우 수렴하지 않음을 확인하였다. 이는 이론에서 설명한 바와 같이 관측 자료와 모델링 자료 사이의 벌림 차이로 인해 나머지(residual) 파동장에 원치 않는 이벤트가 발생하기 때문이다(Fig.
왜곡 현상을 최소화하기 위해 진폭감쇠를 적용한 후 global correlation에 기반한 동시 송신원 기법을 수행하였다. 진폭감쇠는 신호에 감쇠함수를 곱해주어 직접파의 에너지를 강조시키는 효과가 있으므로, 직접파에 가중치를 주고 global correlation을 수행하면 효과적으로 탐사 자료 획득 구조에 맞는 나머지 파동장을 얻을 수 있었다. 주파수 영역에서 생성된 자료를 시간 영역으로 변환하여 나머지 파동장을 계산하는 과정에서 발생하는 두루마리 현상(wrap around)을 억제하기 위해 선처리 작업으로써 뮤팅을 적용하였다.
수치 예제를 통해 진폭감쇠와 global correlation 기법에 기반한 본 알고리듬은 기존의 global correlation 기법만을 적용하는 동시 송신원 전파형 역산에 비해 모델의 심부 및 암염 구조에서 향상된 역산 결과를 생성하며, 스트리머 방식의 탐사 자료에 대한 역산에 적합함을 확인하였다. 진폭감쇠를 적용한 본 파형역산 기술은 global correlation의 특성을 최적화하여 나머지 파동장의 왜곡을 줄이고 파형역산 결과를 향상시키는 것으로 생각된다. 향후 현장 자료에 대한 적용성 연구를 통해 자동화된 최적 감쇠계수 선택 기법을 개발해야 될 것이다.
후속연구
진폭감쇠를 적용한 본 파형역산 기술은 global correlation의 특성을 최적화하여 나머지 파동장의 왜곡을 줄이고 파형역산 결과를 향상시키는 것으로 생각된다. 향후 현장 자료에 대한 적용성 연구를 통해 자동화된 최적 감쇠계수 선택 기법을 개발해야 될 것이다.
질의응답
핵심어
질문
논문에서 추출한 답변
동시 송신원 기법의 문제점은 무엇이며 어떻게 해결할 수 있는가?
이 기법은 다수의 송신원을 동시에 모델링함으로써 기존의 방법에 비해 송신원 개수에 비례하는 계산시간 단축효과가 있다. 물론 다중 송신원의 동시 모델링에 의한 혼선 잡음(crosstalk noise) 등의 문제점이 있지만, 이러한 혼선 잡음은 위상 부호화(Phase encoding) 기법 등으로 해결할 수 있다(Morton and Ober, 1998; Jing et al., 2000; Romero et al.
동시 송신원 기법은 기존의 방법에 비해 어떤 효과가 있는가?
, 2005). 이 기법은 다수의 송신원을 동시에 모델링함으로써 기존의 방법에 비해 송신원 개수에 비례하는 계산시간 단축효과가 있다. 물론 다중 송신원의 동시 모델링에 의한 혼선 잡음(crosstalk noise) 등의 문제점이 있지만, 이러한 혼선 잡음은 위상 부호화(Phase encoding) 기법 등으로 해결할 수 있다(Morton and Ober, 1998; Jing et al.
동시 송신원 기법은 어떤 경우에 추가적인 문제가 발생하는가?
무작위 위상 부호화를 통해 제거할 수 있는 혼선 잡음과는 별개로 동시 송신원 기법은 관측 자료와 모델링 자료 간의 자료 획득 배열(acquisition geometry)이 다를 경우 추가적인 문제가 발생한다. 관측 자료에서 벌림거리의 제한으로 자료가 기록되지 않는 부분에도 모델링 자료는 파동장이 기록되고 이는 나머지 파동장을 심각하게 왜곡시키게 된다.
참고문헌 (16)
Ben-Hadj-Ali, H., Operto, S., and Virieux, J., 2009, Three dimensional frequency-domain full waveform inversion with phase encoding, 79th Annual International Meeting, SEG, Expanded Abstracts, 2288-2292.
Ben-Hadj-Ali, H., Operto, S., and Virieux, J., 2011, An efficient frequency-domain full waveform inversion method using simultaneous encoded sources, Geophysics, 76, R109-R124.
Capdeville, Y., Gung, Y., and Romanowicz, B., 2005, Towards global earth tomography using the spectral element method: A technique based on source stacking, Geophysical Journal International, 162, 541-554.
Choi, Y., and Alkhalifah, T., 2011, Application of encoded multi-source waveform inversion to marine-streamer acquisition based on the global correlation, 73rd Conference and Exhibition, EAGE, Extended Abstracts, F026.
Fomel, S., 2007, Local seismic attributes, Geophysics, 72, A29-A33.
Jing, X., Finn, C. J., Dickens, T. A., and Willen, D. E., 2000, Encoding multiple shot gathers in prestack migration, 70th Annual International Meeting, SEG, Expanded Abstracts, 786-789.
Kim, Y., Cho, H., Min, D.-J., and Shin, C., 2011, Comparison of Frequency-Selection Strategies for 2D Frequency-Domain Acoustic Waveform Inversion, Pure and Applied Geophysics, 168, 1715-1727.
Krebs, J. R., Anderson, J. E., Hinkley, D., Neelamani, R., Lee, S., Baumstein, A., and Lacasse, M.-D., 2009, Fast fullwavefield seismic inversion using encoded sources, Geophysics, 74, WCC177-WCC188.
Morton, S. A., and Ober, C. C., 1998, Faster shot-record depth migration using phase encoding, 68th Annual International Meeting, SEG, Expanded Abstracts, 1131-1135.
Romero, L. A., Ghiglia, D. C., Ober, C. C., and Morton, S. A., 2000, Phase encoding of shot records in prestack migration, Geophysics, 65, 426-436.
Routh, P. S., Krebs, J. R., Lazaratos, S., Baumstein, A. I., Chikichev, I., Lee, S., Downey, N., Hinkley, D., and Andorson, J. E., 2011, Full-wavefield inversion of marine streamer data with the encoded simultaneous source method, 73rd Conference and Exhibition, EAGE, Extended Abstracts, F032.
Shin, C., Pyun, S., and Bednar, J. B., 2007, Comparison of waveform inversion, part 1: conventional wavefield vs logarithmic wavefield, Geophysical Prospecting, 55, 449-464.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.