Next Article in Journal
Monitoring of Atmospheric Carbon Dioxide over a Desert Site Using Airborne and Ground Measurements
Next Article in Special Issue
Structure-Preserving Random Noise Attenuation Method for Seismic Data Based on a Flexible Attention CNN
Previous Article in Journal
An Efficient Method to Compensate Receiver Clock Jumps in Real-Time Precise Point Positioning
Previous Article in Special Issue
Accurately Stable Q-Compensated Reverse-Time Migration Scheme for Heterogeneous Viscoelastic Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Envelope Travel-Time Objective Function for Reducing Source–Velocity Trade-Offs in Wave-Equation Tomography

1
Key Laboratory of Petroleum Resource Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Innovation Academy of Earth Science, Chinese Academy of Sciences, Beijing 100029, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(20), 5223; https://doi.org/10.3390/rs14205223
Submission received: 5 September 2022 / Revised: 15 October 2022 / Accepted: 17 October 2022 / Published: 19 October 2022
(This article belongs to the Special Issue Geophysical Data Processing in Remote Sensing Imagery)

Abstract

:
In conventional cross-correlation (CC)-based wave-equation travel-time tomography, wrong source wavelets can result in inaccurate velocity inversion results, which is known as the source–velocity trade-off. In this study, an envelope travel-time objective function is developed for wave-equation tomography to alleviate the non-uniqueness and uncertainty due to wrong source wavelets. The envelope of a seismic signal helps reduce the waveform fluctuations/distortions caused by variations of the source time function. We show that for two seismic signals generated with different source wavelets, the travel-time shift calculated by cross-correlation of their envelopes is more accurate compared to that obtained by directly cross-correlating their waveforms. Then, the CC-based envelope travel-time (ET) objective function is introduced for wave-equation tomography. A new adjoint source has also been derived to calculate the sensitivity kernels. In the numerical inversion experiments, a synthetic example with cross-well survey is first given to show that compared to the traditional CC travel-time objective function, the ET objective function is relatively insensitive to source wavelet variations and can reconstruct the elastic velocity structures more reliably. Finally, the effectiveness and advantages of our method are verified by inversion of early arrivals in a practical seismic survey for recovering near-surface velocity structures.

Graphical Abstract

1. Introduction

In global- and exploration-scale seismology, compared to ray-tracing tomography methods, wave-equation-based approaches show stronger capability to recover subsurface velocity heterogeneities with higher accuracy and spatial-resolution [1,2,3,4,5,6,7,8,9]. Due to the lack of low frequencies in seismic recordings and poor initial models, the standard waveform-difference objective function in wave-equation tomography or full-waveform inversion (FWI) is prone to cycle-skipping [10,11,12,13,14,15]. Because travel-time is linearly related to the velocity perturbations, cross-correlation (CC)-based wave-equation travel-time tomography is widely used to reconstruct low-wavenumber components of the velocity models [16,17,18,19,20,21,22]. High-wavenumber model details can then be obtained using FWI with the waveform-difference objective function based on the macro background structures provided by wave-equation travel-time tomography.
In most practical applications of wave-equation based imaging/inversion methods, the source signatures (e.g., source time function, mechanism, and location) are unknown. In wave-equation travel-time tomography, if the source wavelet used in forward modeling is not accurate, the travel-time shifts calculated by cross-correlating the synthetic and observed seismic signals are also erroneous, producing unexpected uncertainties in the inverted velocity models. This issue is referred to as the source–velocity (or structure) trade-off, which has bothered geophysicists in seismic tomography for decades [23,24,25,26,27,28]. Fichtner [24] quantified source–structure trade-offs in ambient noise tomography with off-diagonal Hessian elements and mentioned that the inherent trade-offs may lead to mispositioning of structural heterogeneities. Joint inversion of the source wavelet and velocity structures represents one potential strategy to overcome this problem. However, simultaneous inversion of model parameters with different dimensionalities also suffers from the problem of interparameter trade-off [29,30,31], leading to inaccurate velocity and source inversion results. Truncated-Newton (or Gauss–Newton) optimization methods are promising to reduce these trade-off uncertainties by applying the inverse Hessian implicitly [32,33]. However, calculating Hessian-vector products in linear conjugate-gradient methods increases the algorithm complexities and computation cost significantly [32,33,34]. Yuan et al., [26] developed a double-difference objective function for wave-equation travel-time tomography with differential travel-time measurements, which can obtain higher resolution velocity structures independent of the source wavelet. However, within this approach, as the number of sources and receivers increases, the computation burden increases significantly.
The envelope measures the instantaneous amplitude of a seismic signal [35]. Wu et al., [11] recognized that the envelope fluctuation and decay of the seismic records carries ultra-low frequencies that are essentially valuable for recovering long-wavelength model structures. In FWI studies and practices, the objective functions minimizing envelope differences of synthetic and observed seismograms are widely used to suppress non-linearity of the inverse problem [11,12,13,36,37,38]. Pan and Innanen [39] revealed that the envelope-based objective functions help resolve the attenuation structures more stably. In this research, we find that envelope of the seismic signal can reduce the effects of waveform fluctuations/distortions caused by source wavelet changes. The travel-time shifts calculated by cross-correlating the seismic envelopes are relatively insensitive to wavelet errors. Thus, a CC-based envelope travel-time (ET) objective function is developed for wave-equation tomography. Derivation of the corresponding new adjoint source for calculating the velocity sensitivity kernels is provided. This ET objective function is expected to construct the velocity structures more reliably, suffering from fewer trade-off uncertainties due to inaccurate source wavelets. Compared to the Newton-based optimization methods and double-difference adjoint tomography approach, this method can be implemented more easily without increasing the computational requirements. Furthermore, this approach inherits the advantage of the envelope-difference objective function for mitigating cycle-skipping. Using the cross-correlation-based ET objective function for reducing source–velocity trade-offs has not been investigated by other researchers and represents the innovative contribution of this study.
In the numerical section, we first provide a synthetic example of elastic wave-equation travel-time tomography with a cross-well survey. Inversion experiments using the conventional CC and ET objective functions with correct- and wrong source wavelets are performed for comparison. The output results indicate that when wrong source wavelet is used in traditional CC-based wave-equation travel-time tomography, the inverted P-wave velocity V P and S-wave velocity V S perturbations are blurred with inaccurate values. However, when using the ET objective function, the elastic velocity perturbations can be recovered accurately with both correct and wrong source wavelets. This suggests that this method is less sensitive to errors of source wavelets and can reduce the trade-off uncertainties. Synthetic inversion experiments are also given to show that with large velocity contrasts, the ET objective function still works effectively. To examine the performances of the ET objective function further, we apply it to practical seismic data acquired from southeast Poland. Early arrivals are extracted from the recordings for imaging the near-surface V P model variations. It is observed that the ET objective function can resolve more detailed velocity structures with inaccurate source wavelets, suggesting that this approach works effectively and stably for noisy field data.
In this paper, the ET objective function is first introduced and the corresponding adjoint source is derived following the chain rule. The theory of conventional CC-based wave-equation travel-time tomography is given in Appendix A. Then, both synthetic inversion experiments and practical field data application are provided to verify the effectiveness and advantages of the presented method.

2. Methods

Envelope Travel-Time Wave-Equation Tomography

In wave-equation-based tomography/inversion methods, the standard objective function measures the direct waveform differences between the synthetic data u i and seismic observation d i . Because the relationship between model perturbations and waveform differences is highly non-linear, this waveform-difference objective function is easily trapped at the local minimum, referring to as the cycle-skipping problem. Thus, the cross-correlation-based travel-time objective function was developed for reducing the non-linearity of wave-equation tomography [16]. The theory of conventional wave-equation travel-time tomography is introduced in Appendix A. However, when the source wavelet is inaccurate, the inverted velocity structures suffer from source–velocity trade-offs. In this section, we introduce an envelope travel-time wave-equation tomography approach for reducing the influences of wrong source wavelets.
An analytic signal s i can be constructed from a real signal s i and its Hilbert transform:
s i t = s i t i H s i t ,
where i denotes the imaginary unit and H s i t means Hilbert transform of s i :
H s i t = 1 π + s i τ t τ d τ .
The analytic signal can be written in terms of instantaneous amplitude E i and instantaneous phase φ i :
s i t = E i t exp i φ i t ,
where E i = R s i t 2 + I s i t 2 ( R and I mean real and imaginary parts) represents the envelope of the analytic signal. In FWI studies and applications, the envelope-difference objective function is commonly used to reduce the cycle-skipping problem. In this research, we reveal that the envelope of the seismic signal can reduce the fluctuations/distortions of waveforms caused by inaccurate source wavelets.
Figure 1 shows a minimum phase wavelet s i with a dominant frequency of 10 Hz (solid-red line) and a wrong source wavelet s i w (solid-black line). The wrong source wavelet is generated by applying phase and amplitude variations to the correct minimum phase wavelet combining with another different wavelet. The dash-dotted-red and dash-dotted-black lines indicate their envelopes ( s i , e n v and s i , e n v w ), respectively. Δ s i (solid-blue line) is the difference between these two wavelets. It can be seen that the waveforms of the two wavelets are obviously different. However, their envelopes can reduce the phase and amplitude distortions. We generate observed trace d i and synthetic trace u i c using the same minimum phase wavelet, as shown in Figure 2a,b. d i , e n v and u i , e n v c are their envelopes, respectively. Travel-time difference between the observed and synthetic traces is exactly 0.2 s. In Figure 2c, the synthetic trace u i w with the wrong wavelet and its envelope u i , e n v w are plotted. Then, we calculate the cross-correlation c i between d i and u i c , as plotted in Figure 2d. An exact travel-time shift 0.2 s can be obtained. c i w = d i u i w means the cross-correlation between the observed trace and the synthetic trace with the wrong wavelet (⊗ means cross-correlation). The travel-time shift is obtained as 0.187 s with 6.5% error. c i , e n v w = d i , e n v u i , e n v w is the cross-correlation between the observed trace envelope and the synthetic trace envelope with the wrong wavelet. The obtained travel-time shift is 0.2004 s with 0.2% error. This means that when the source wavelet is wrong, cross-correlation of the envelopes can produce a more accurate travel-time shift.
Thus, in this study, we introduce the ET objective function for wave-equation tomography to mitigate the influence of inaccurate source wavelets in velocity inversion:
ϕ e t m = 1 2 i Δ t i e t 2 ,
where the travel-time difference Δ t i e t is obtained by correlating the envelopes of synthetic and observed seismic traces:
Δ t i e t = arg max τ 0 T E i t + τ E i o b s t d t ,
where E i and E i o b s are the envelopes of synthetic and observed seismic data. Variation of the ET objective function with respect to the model perturbation is given by:
Δ ϕ e t = Δ t i e t Δ Δ t i e t ,
where Δ Δ t i e t means the variation of travel-time shift corresponding to model perturbation:
Δ Δ t i e t = 1 N e t t E i t Δ E i t ,
where N e t in the denominator is:
N e t = 0 T t 2 E i t E i t d t ,
and Δ E i represents the variation of envelope:
Δ E i t = u i t Δ u i t + H u i t Δ H u i t E i t .
Substituting Equation (9) into Equation (7) gives [40]:
Δ Δ t i e t = t E i t u i t N e t E i t Δ u i t H t E i t H u i t N e t E i t Δ u i t .
Inserting Equation (10) into Equation (6), we obtain
Δ ϕ e t = Δ t i e t t E i t u i t N e t E i t Δ u i t + Δ t i e t H t E i t H u i t N e t E i t Δ u i t .
Thus, we can derive the adjoint source of the ET objective function for calculating the velocity sensitivity kernels as:
s i , e t * t = Δ t i e t P i T t u i T t + Δ t i e t H P i T t H u i T t ,
where P i is
P i t = t E i t N e t E i t .
This ET objective function can be implemented easily and is expected to recover the velocity structures more accurately by reducing the trade-off uncertainties caused by wrong source wavelets.

3. Numerical Results

In this section, we present synthetic inversion experiments and field data application to verify the advantages of ET objective function for recovering velocity structures when an inaccurate source wavelet is used in forward simulations and adjoint tomography. We carry out the inversion experiments using the software package SeisElastic2D [41].

3.1. Synthetic Inversion Example with Cross-Well Survey

In the first numerical example, we design synthetic elastic velocity models with cross-well survey. Figure 3a,b show the target V P and V S velocity models with V P / V S = 2 . Velocity anomalies at the center parts of the target models are generated with ±10% perturbations of the background values. The initial models are homogeneous without the velocity perturbations at the center parts ( V P = 2 V S = 3500 m/s). The density model is homogeneous with ρ = 1900 kg/m 3 . Because there is no density contrast in the elastic models, boundaries of the inversion results are not affected. A number of 14 sources and 180 receivers are arranged regularly on the left and right sides of the model, as shown in Figure 3a. The observed data are generated using the spectral-element forward modeling method with a minimum phase wavelet of 50 Hz dominant frequency, as indicated by the red line in Figure 3c. Another wrong source wavelet is also created by applying phase and amplitude variations to the correct minimum phase wavelet (black line in Figure 3c).
When applying the inversion, P- and S-waves are separated in the shot profiles to invert for V P and V S models, respectively. In Figure 4a,b, we plot the discrete traces of P- and S-waves extracted from the observed data (vertical component), synthetic data with the wrong source wavelet, and their envelopes for comparison. Envelopes of the seismic traces can cancel the effects of phase variations and reduce the waveform distortions effectively. In Figure 5, the absolute travel-time shifts calculated by cross-correlating the seismic traces and their envelopes are provided for comparison. It can be seen that when using wrong source wavelet, the travel-time shifts obtained from cross-correlating envelopes are closer to those obtained with the correct source wavelet, which means that the ET objective function is more robust for calculating the travel-time shifts.
Figure 6 presents the inverted V P model perturbations using traditional CC and ET objective functions with the correct and wrong source wavelets, respectively. In Figure 7, the inverted V S model perturbations by different objective functions and source wavelets are given for comparison. It can be seen that when using the wrong source wavelet, the inverted Δ V P c c , w and Δ V S c c , w structures (Figure 6b and Figure 7b) from the traditional CC objective function are blurred. The velocity anomalies are not clearly resolved. This observation is consistent with that presented in Yuan et al., [26]. However, when using the ET objective function, the inverted Δ V P e t , w and Δ V S e t , w model perturbations (Figure 6d and Figure 7d) with the wrong source wavelet are comparable to those obtained with the correct wavelet. Figure 8 shows the reductions of normalized travel-time data misfits when using different objective functions and source wavelets. With the wrong source wavelet, the inversion with the traditional CC objective function stops after fewer iterations, with smaller reductions of travel-time data misfits. These observations suggest that compared to the traditional CC objective function, this ET objective function algorithm is less sensitive to source wavelet errors and thus can produce more accurate velocity estimations.
To examine the effectiveness and robustness of the ET objective function for recovering velocity structures with large contrasts, we increase the velocity perturbations at the center parts of the models from ±10% of background values to ±30%. In this condition, the inverse problem becomes more non-linear. The inversion experiments are then performed using the ET objective function with the wrong source wavelet. The reconstructed model perturbations of P- and S-wave velocities are presented in Figure 9a,b. It can be seen that the elastic velocity anomalies can still be reliably inverted. This experiment suggests that the ET objective function can work effectively for large velocity contrasts with wrong source wavelets by reducing cycle-skipping and source–velocity trade-offs.

3.2. Field Data Application

To illustrate the effectiveness and stability of the ET objective function for practical seismic data, in this section, we apply wave-equation travel-time tomography with different source wavelets to image the near-surface velocity heterogeneities using early arrivals of the practical recordings. This dataset was acquired from Poland with Vibroseis seismic sources and is publicly available.
The initial 1D V P model is created with conventional migration velocity analysis, as shown in Figure 10. The initial V S and ρ models are generated with V S = V P / 2.5 and ρ = 350 × V P 0.25 (Gardner’s rule) and kept unchanged in the inversion. The model is 13.5 km wide and 1.0 km deep. We extract 27 sources with an interval of 200 m from the raw dataset for inversion. The receiver spacing is 25 m. Figure 11a shows one raw shot gather. It is observed that the raw shot is contaminated by strong surface waves, but contains clear early arrivals, which are valuable for recovering near-surface V P structures. The shot gathers are pre-processed to obtain the early arrivals with a series of operations including de-noising, band-pass filtering ([10 Hz, 20 Hz]), amplitude normalization, etc., as presented in Figure 11b. Because the Vibroseis source (8–95 Hz linear sweep over 15 s) was used to generate the field data, the zero-phase Klauder source wavelet s i k can be directly generated, as indicated by the red line in Figure 12. To examine the performance of the traditional CC and the ET objective functions, we create another wrong source wavelet s i w by applying phase and amplitude variations to the Klauder wavelet for inversion, as indicated by the black line in Figure 12.
We first perform inversion experiments using the traditional CC objective function with the correct Klauder wavelet and wrong source wavelet. The inversion results ( V P c c , k and V P c c , w ) are plotted in Figure 13a,b for comparison. The corresponding model perturbations ( Δ V P c c , k and Δ V P c c , w ) are presented in Figure 14a,b, respectively. Data misfits reduction histories are plotted in Figure 15. It can be seen that when using the correct Klauder wavelet, detailed velocity structures at the near-surface (0.0–1.0 km) can be resolved. A low-velocity zone at shallow parts (approximately 0–0.2 km) and a high-velocity formation at depths of approximately 0.3–0.7 km can be imaged clearly. Due to the limitation of source illumination, only the velocity structures at horizontal distance of 2.0–10.0 km are recovered. However, when using the wrong source wavelet, the inversion stops with fewer iterations, as indicated by the dash-dotted-black line in Figure 15. Only the low-velocity zone at shallow depth is recovered. The whole velocity model is blurred and inaccurate. Then, we carry out inversion experiments using the ET objective function with the correct and wrong source wavelets. Figure 13c,d show the inverted models ( V P e t , k and V P e t , w ), respectively. Figure 14c,d present the corresponding model perturbations ( Δ V P e t , k and Δ V P e t , w ). We observe that the inverted velocity structures by the ET objective function are very close to those obtained by the traditional CC objective function with the correct source (Figure 13a). Even though the wrong source wavelet is used, the ET objective function can still recover the low- and high-velocity zones clearly and reduce the travel-time shifts significantly, as indicated by the dash-dotted-red line in Figure 15. In Figure 16a, the observed traces, synthetic traces calculated using initial model and wrong source wavelet, and the synthetic traces calculated from the inverted model by the traditional CC objective function are plotted for comparison. After inversion, the travel-time shifts between the synthetic traces and observed traces are still very large. However, the envelopes of the synthetic traces calculated from the inverted model by the ET objective function match the envelopes of the observed traces closely, as indicated by the gray arrows in Figure 16b. These observations suggest that compared to the traditional CC objective function, the ET objective function is less sensitive to the variation of the source wavelet and thus can reduce source–velocity trade-offs. This field data application example also suggests that our new method can work effectively for other non-Vibroseis datasets with unknown source wavelets.
Finally, we design a synthetic inversion example to examine the credibility of the inversion results presented in Figure 13, which is equivalent to a resolution test. The target model is created by adding the velocity perturbation Δ V P (Figure 17) to the initial V P model (Figure 10). The initial model used for inversion is the same as that (Figure 10) used in the practical data applications. Both the traditional CC and ET objective functions are used to recover the velocity perturbations with the correct Klauder wavelet and wrong source wavelet (Figure 12). The inversion parameters and settings are the same as those in the practical applications described above. The recovered model perturbations are presented in Figure 18. Compared to the traditional CC objective function, the ET objective function can recover the velocity anomalies more accurately and reliably, which is consistent with the observations in the practical data application. This test means that the inversion results in Figure 13 by different objective functions and source wavelets are convincing.

4. Discussions

In this study, we discuss and investigate the problem of source–velocity trade-offs in wave-equation travel-time tomography. This issue also exists in FWI with the waveform-difference objective function. Different source-independent algorithms have been developed for FWI [42,43,44,45]. However, these strategies may not be effective for inversion of seismic travel-time measurements and also suffer from the cycle-skipping challenge. Compared to the double-difference travel-time objective function proposed by Yuan et al., [26], our method can be implemented more easily without increasing the computation cost.
In wave-equation travel-time tomography, the inversion results are also affected by other source signatures including source mechanism and location. Here, the ET objective function only works for reducing the trade-offs caused by inaccurate source time functions. In travel-time tomography applications, P- and S-waves are always separated to invert for P- and S-wave velocities, respectively. In this condition, the influences of the source mechanism are limited. In exploration seismology, the source locations are always known. However, in earthquake seismology, the effects caused by inaccurate source locations should be carefully addressed.
In the inversion practices with traditional CC and ET objective functions, we realize that velocity perturbations also result in amplitude fluctuations of the waveforms. These effects produce errors in the calculation of travel-time shifts and estimation of velocity structures [20]. Furthermore, it is noticed that the ET objective function is more robust to these amplitude distortions. However, in some cases, the ET objective function does not work well for reducing the influences of wrong source wavelets. For examples, when the source wavelets contain many irregular slide lobes or the correct and wrong source wavelets have different frequency bands, this approach cannot produce very accurate velocity inversion results. In the large velocity contrasts example, the ET objective function inverts for the models effectively. However, if the model perturbations are increased further, the ET objective function may fail due to cycle-skipping.

5. Conclusions

In traditional cross-correlation based wave-equation travel-time tomography, the inversion results suffer from uncertainties caused by source–velocity trade-offs. In this research, we have revealed that the envelope of a seismic signal helps in reducing the amplitude and phase distortions caused by inaccurate source wavelets. The travel-time shifts calculated by cross-correlating envelopes are less sensitive to the variations of source wavelets. Thus, an envelope travel-time objective function is designed for wave-equation tomography. Derivation of the adjoint source for this envelope travel-time objective function is also provided. Compared to the traditional approach, this method is expected to recover the velocity structures more accurately with wrong source wavelets by suppressing the effects of source–velocity trade-offs. In the synthetic inversion experiments and field data applications, we verify that when a wrong source wavelet with phase and amplitude changes is used in forward modeling and tomography, the envelope travel-time objective function can still calculate the travel-time shifts accurately and invert for the velocity models reliably.

Author Contributions

Conceptualization, W.P. and Y.W.; methodology, W.P. and N.M.; validation, W.P. and N.M.; resources, Y.W.; writing—original draft preparation, W.P.; writing—review and editing, N.M. and Y.W.; supervision, Y.W.; project administration, Y.W.; funding acquisition, Y.W. and W.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (grant no. 42004116, 12171455), the Original Innovation Research Program of the Chinese Academy of Sciences (CAS) under grant no. ZDBS-LY-DQC003, and IGGCAS Research Start-up Funds (grant no. 2020000061).

Data Availability Statement

The synthetic data can be obtained by contacting the corresponding author. The field seismic data are available at https://www.freeusp.org/RaceCarWebsite/TechTransfer/Tutorials/Processing_2D/ (accessed on 15 January 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Conventional Wave-Equation Travel-Time Tomography

In conventional wave-equation-based travel-time tomography, the velocity structures are estimated by iteratively minimizing the travel-time differences between synthetic and observed measurements:
ϕ c c m , s = 1 2 i Δ t i 2 ,
where m and s mean the structure and source models, respectively, and travel-time difference Δ t i is obtained by cross-correlating the synthetic data u i and observed data d i :
Δ t i = arg max τ 0 T u i t + τ d i t d t ,
where T is the maximum recording time. The synthetic traces are obtained by numerically solving the acoustic/elastic wave-equations. Thus, seismic wave-equation inversion/tomography can be regarded as a PDE (partial differential equation)-constrained inverse problem. Based on the framework of the adjoint-state theorem, the objective function corresponding to the relative model perturbations ( Δ m / m and Δ s / s ) can be written as [4]:
Δ ϕ c c m , s = x K m Δ m + K s Δ s d x ,
where ⊕ means the space of studied area in subsurface, K m and K s are the sensitivity kernels (or gradients) for model structures and source, respectively. Following the adjoint-state method with divergence theorem and integration by parts [46], the sensitivity kernels for P-wave velocity V P and S-wave velocity V S in isotropic-elastic media can be derived as [31]:
K V P = 2 ρ V P 2 i u i * k u k ,
K V S = 2 ρ V S 2 j u i * i u j + j u i 2 i u i * k u k .
where u i and u i * indicate the forward and adjoint (or back-propagated) displacement wavefields, respectively, and summation convention is being used for the subscripts of these displacement wavefields. Note that in the expressions of sensitivity kernels, the integrations over source, receiver, space, and time are ignored for the sake of compactness. The adjoint wavefields can be obtained by solving the adjoint-state equation with the following adjoint source:
s i , c c * t = 1 N Δ t i t u i T t , N = 0 T t 2 u i t u i t d t .
Travel-time is linearly related to velocity variations; thus, the CC-based travel-time objective function is advantageous to overcome the cycle-skipping difficulty.
In many practical seismic imaging/inversion applications, the source signatures (e.g., source mechanism and wavelet) are unknown. When an inaccurate source wavelet is used in forward modeling, cross-correlation of the synthetic and observed seismic traces produces inaccurate travel-time shifts, resulting in unconvincing velocity models. In the vicinity of the global minimum of ϕ c c , incorporating the second-order derivative Hessian H gives [25]:
ϕ c c m + Δ m , s + Δ s ϕ c c m , s + Δ m Δ s H m m H m s H s m H s s Δ m Δ s ,
where † means matrix transpose and the off-diagonal block matrices H m s and H s m contain mixed derivatives of ϕ c c with respect to the structure and source. The Hessian-vector product H m s Δ s measures the trade-off artifacts in the velocity sensitivity kernels due to source errors [31]. Theoretically, the trade-offs can be effectively reduced by applying the inverse Hessian implicitly with linear (preconditioned) conjugate-gradient algorithms in a truncated (Gauss)-Newton iteration [33]. However, for structure and source models with different dimensionalities and units, it is complex and expensive to calculate the Hessian-vector products and perform inner conjugate-gradient iterations. In this study, we adopt the quasi-Newton l-BFGS optimization method [47,48,49] for approximating the inverse Hessian and updating models.

References

  1. Tarantola, A. Inversion of seismic reflection data in the acoustic approximation. Geophysics 1984, 49, 1259–1266. [Google Scholar] [CrossRef]
  2. Woodward, M. Wave-equation tomography. Geophysics 1992, 57, 231–248. [Google Scholar] [CrossRef]
  3. Pratt, R.G.; Shin, C.; Hicks, G.J. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion. Geophys. J. Int. 1998, 133, 341–362. [Google Scholar] [CrossRef]
  4. Tromp, J.; Tape, C.; Liu, Q. Seismic tomography, adjoint methods, time reversal, and banana-doughnut kernels. Geophysics 2005, 160, 195–216. [Google Scholar] [CrossRef] [Green Version]
  5. Xie, X.; Jin, S.; Wu, R. Wave-equation-based seismic illumination analysis. Geophysics 2006, 71, S169–S177. [Google Scholar] [CrossRef]
  6. Tape, C.; Liu, Q.; Maggi, A.; Tromp, J. Adjoint tomography of the southern California crust. Science 2009, 325, 988–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wu, R.; Wang, B.; Hu, C. Renormalized nonlinear sensitivity kernel and inverse thin-slab propagator in T-matrix formalism for wave-equation tomography. Inverse Probl. 2015, 31, 115004. [Google Scholar] [CrossRef] [Green Version]
  8. Bozdağ, E.; Peter, D.; Lefebvre, M.; Komatitsch, D.; Tromp, J.; Hill, J.; Podhorszki, N.; Pugmire, D. Global adjoint tomography: First-generation model. Geophys. J. Int. 2016, 207, 1739–1766. [Google Scholar] [CrossRef] [Green Version]
  9. Operto, S.; Miniussi, A. On the role of density and attenuation in 3D multi-parameter visco-acoustic VTI frequency-domain FWI: An OBC case study from the North Sea. Geophys. J. Int. 2018, 213, 2037–2059. [Google Scholar] [CrossRef] [Green Version]
  10. Bunks, C.; Saleck, F.M.; Zaleski, S.; Chavent, G. Multiscale seismic waveform inversion. Geophysics 1995, 60, 1457–1473. [Google Scholar] [CrossRef]
  11. Wu, R.; Luo, J.; Wu, B. Seismic envelope inversion and modulation signal model. Geophysics 2014, 79, WA13–WA24. [Google Scholar] [CrossRef]
  12. Yuan, Y.O.; Simons, F.J.; Bozdaǧ, E. Multiscale adjoint waveform tomography for surface and body waves. Geophysics 2015, 80, R281–R302. [Google Scholar] [CrossRef] [Green Version]
  13. Luo, J.; Wu, R. Seismic envelope inversion: Reduction of local minima and noise resistance. Geophys. Prospect. 2015, 63, 597–614. [Google Scholar] [CrossRef]
  14. Métivier, L.; Brossier, R.; Mérigot, Q.; Oudet, E.; Virieux, J. Measuring the misfit between seismograms using an optimal transport distance: Application to full waveform inversion. Geophys. J. Int. 2016, 205, 345–377. [Google Scholar] [CrossRef]
  15. Yao, G.; Wu, D.; Wang, S. A review of reflection-waveform inversion. Pet. Sci. 2020, 17, 334–351. [Google Scholar] [CrossRef] [Green Version]
  16. Luo, Y.; Schuster, G.T. Wave-equation traveltime inversion. Geophysics 1991, 56, 645–653. [Google Scholar] [CrossRef]
  17. Zhou, C.; Schuster, G.T.; Hassanzadeh, S.; Harris, J.M. Elastic wave equation traveltime and waveform inversion of crosswell data. Geophysics 1997, 62, 853–868. [Google Scholar] [CrossRef]
  18. Leeuwen, T.; Mulder, W.A. A correlation-based misfit criterion for wave-equation traveltime tomography. Geophys. J. Int. 2010, 182, 1383–1394. [Google Scholar] [CrossRef] [Green Version]
  19. Chi, B.; Dong, L.; Liu, Y. Correlation-based reflection full-waveform inversion. Geophysics 2015, 80, R189–R202. [Google Scholar] [CrossRef]
  20. Luo, Y.; Ma, Y.; Wu, Y.; Liu, H.; Gao, L. Full-traveltime inversion. Geophysics 2016, 81, R261–R274. [Google Scholar] [CrossRef]
  21. Zheng, Y.; Wang, Y.; Luo, Q.; Chang, X.; Zeng, R.; Wang, B.; Zhao, X. Frequency-dependent reflection wave-equation traveltime inversion from walkaway vertical seismic profile data. Geophysics 2019, 84, R947–R961. [Google Scholar] [CrossRef]
  22. Feng, B.; Xu, W.; Wu, R.; Xie, X.; Wang, H. Finite-frequency traveltime tomography using the Generalized Rytov approximation. Geophys. J. Int. 2020, 221, 1412–1426. [Google Scholar] [CrossRef]
  23. Zhang, H.; Thurber, C.H. Double-difference tomography: The method and its application to the Hayward Fault, California. Bull. Seismol. Soc. Am. 2003, 3, 1875–1889. [Google Scholar] [CrossRef]
  24. Fichtner, A. Full Seismic Waveform Inversion for Structural and Source Parameters. Ph.D. Thesis, Ludwig Maximilian University, Munich, Germany, 2010. [Google Scholar]
  25. Fichtner, A. Source-structure trade-offs in ambient noise correlations. Geophys. J. Int. 2015, 202, 678–694. [Google Scholar] [CrossRef] [Green Version]
  26. Yuan, Y.O.; Simons, F.J.; Tromp, J. Double-difference adjoint seismic tomography. Geophys. J. Int. 2016, 206, 1599–1618. [Google Scholar] [CrossRef] [Green Version]
  27. Sager, K.; Ermert, L.; Boehm, C.; Fichtner, A. Towards full waveform ambient noise inversion. Geophys. J. Int. 2018, 212, 566–590. [Google Scholar] [CrossRef]
  28. Blom, N.; Hardalupas, P.S.; Rawlinson, N. Mitigating the effect of errors in source parameters on seismic (waveform) tomography. Geophys. J. Int. 2022, 232, 810–828. [Google Scholar] [CrossRef]
  29. Operto, S.; Gholami, Y.; Prieux, V.; Ribodetti, A.; Brossier, R.; Métivier, L.; Virieux, J. A guided tour of multiparameter full waveform inversion with multicomponent data: From theory to practice. Lead. Edge 2013, 32, 1040–1054. [Google Scholar] [CrossRef]
  30. Pan, W.; Innanen, K.A.; Margrave, G.F.; Fehler, M.C.; Fang, X.; Li, J. Estimation of elastic constants for HTI media using Gauss-Newton and full-Newton multiparameter full-waveform inversion. Geophysics 2016, 81, R275–R291. [Google Scholar] [CrossRef]
  31. Pan, W.; Geng, Y.; Innanen, K.A. Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: Synthetic experiments and Hussar data set application. Geophys. J. Int. 2018, 213, 1305–1333. [Google Scholar] [CrossRef]
  32. Métivier, L.; Bretaudeau, F.; Brossier, R.; Virieux, J.; Operto, S. Full waveform inversion and the truncated Newton method: Quantitative imaging of complex subsurface structures. Geophys. Prospect. 2014, 62, 1353–1375. [Google Scholar] [CrossRef] [Green Version]
  33. Pan, W.; Innanen, K.A.; Liao, W. Accelerating Hessian-free Gauss-Newton full-waveform inversion via l-BFGS preconditioned conjugate-gradient algorithm. Geophysics 2017, 32, R49–R64. [Google Scholar] [CrossRef]
  34. Epanomeritakis, I.; Akçelik, V.; Ghattas, O.; Bielak, J. A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion. Inverse Probl. 2008, 24, 034015. [Google Scholar] [CrossRef]
  35. Taner, M.T.; Koehler, F.; Sheriff, R.E. Complex seismic trace analysis. Geophysics 1979, 44, 1041–1063. [Google Scholar] [CrossRef]
  36. Chen, G.; Wu, R.; Wang, Y.; Chen, S. Multi-scale signed envelope inversion. J. Appl. Geophys. 2018, 153, 113–126. [Google Scholar] [CrossRef]
  37. Gao, Z.; Pan, Z.; Gao, J.; Wu, R. Frequency controllable envelope operator and its application in multiscale full-waveform inversion. IEEE Trans. Geosci. Remote Sens. 2019, 57, 683–699. [Google Scholar] [CrossRef]
  38. Hu, Y.; Wu, R.; Huang, X.; Long, Y.; Xu, Y.; Han, L. Phase-amplitude-based polarized direct envelope inversion in the time-frequency domain. Geophysics 2022, 87, R245–R260. [Google Scholar] [CrossRef]
  39. Pan, W.; Innanen, K.A. Amplitude-based misfit functions in viscoelastic full-waveform inversion applied to walk-away vertical seismic profile data. Geophysics 2019, 84, B335–B351. [Google Scholar] [CrossRef]
  40. Bozdağ, E.; Trampert, J.; Tromp, J. Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements. Geophys. J. Int. 2011, 185, 845–870. [Google Scholar] [CrossRef]
  41. Pan, W.; Innanen, K.A.; Wang, Y. SeisElastic2D: An open-source package for multiparameter full-waveform inversion in isotropic-, anisotropic- and visco-elastic media. Comput. Geosci. 2020, 145, 104586. [Google Scholar] [CrossRef]
  42. Kim, H.J.; Lee, K.H. Source-independent full-waveform inversion of seismic data. Geophysics 2003, 68, 2010–2015. [Google Scholar]
  43. Xu, K.; Greenhalgh, S.A.; Wang, M. Comparison of source-independent methods of elastic waveform inversion. Geophysics 2006, 71, R91–R100. [Google Scholar] [CrossRef]
  44. Alkhalifah, T.; Choi, Y. Source-independent time-domain waveform inversion using convolved wavefields: Application to the encoded multisource waveform inversion. Geophysics 2011, 76, R125–R134. [Google Scholar]
  45. Zhang, Q.; Zhou, H.; Li, Q.; Chen, H.; Wang, J. Robust source-independent elastic full-waveform inversion in the time domain. Geophysics 2016, 81, R29–R44. [Google Scholar] [CrossRef]
  46. Liu, Q.; Tromp, J. Finite-frequency kernels based on adjoint methods. Bull. Seismol. Soc. Am. 2006, 96, 2383–2397. [Google Scholar] [CrossRef] [Green Version]
  47. Byrd, R.H.; Lu, P.; Nocedal, J. A limited memory algorithm for bound constrained optimization. SIAM J. Sci. Comput. 1995, 16, 1190–1208. [Google Scholar] [CrossRef]
  48. Nocedal, J. Updating quasi-Newton matrices with limited storage. Math. Comput. 1980, 35, 773–782. [Google Scholar] [CrossRef]
  49. Nocedal, L.; Wright, S.J. Numerical Optimization, 2nd ed.; Springer: Cham, Switzerland, 2006; pp. 30–63. [Google Scholar]
Figure 1. Comparison of the source wavelets and their envelopes. The solid-red and dash-dotted-red lines indicate the original minimum phase wavelet ( s i ) and its envelope ( s i , e n v ). The solid-black and dash-dotted-black lines indicate the wrong source wavelet ( s i w ) and its envelope ( s i , e n v w ), respectively. The solid-blue line indicates the difference between the correct and wrong source wavelets ( Δ s i ).
Figure 1. Comparison of the source wavelets and their envelopes. The solid-red and dash-dotted-red lines indicate the original minimum phase wavelet ( s i ) and its envelope ( s i , e n v ). The solid-black and dash-dotted-black lines indicate the wrong source wavelet ( s i w ) and its envelope ( s i , e n v w ), respectively. The solid-blue line indicates the difference between the correct and wrong source wavelets ( Δ s i ).
Remotesensing 14 05223 g001
Figure 2. (a) Observed trace d i (solid-red line) and its envelope d i , e n v (dash-dotted-red line); (b) synthetic trace u i c (solid-red line) with correct wavelet and its envelope u i , e n v c (dash-dotted-red line); (c) synthetic trace u i w with wrong wavelet (solid-black line) and its envelope u i , e n v w (dash-dotted-black line); (d) c i is the cross-correlation between observed trace d i and the synthetic trace u i c (solid-red line). c i w is the cross-correlation between d i and u i w (solid-black line). c i , e n v w is the cross-correlation between d i , e n v and u i , e n v w (dash-dotted-black line).
Figure 2. (a) Observed trace d i (solid-red line) and its envelope d i , e n v (dash-dotted-red line); (b) synthetic trace u i c (solid-red line) with correct wavelet and its envelope u i , e n v c (dash-dotted-red line); (c) synthetic trace u i w with wrong wavelet (solid-black line) and its envelope u i , e n v w (dash-dotted-black line); (d) c i is the cross-correlation between observed trace d i and the synthetic trace u i c (solid-red line). c i w is the cross-correlation between d i and u i w (solid-black line). c i , e n v w is the cross-correlation between d i , e n v and u i , e n v w (dash-dotted-black line).
Remotesensing 14 05223 g002
Figure 3. (a,b) are the true V P and V S model structures; (c) The red and black lines indicate the correct and wrong source wavelets ( s i c and s i w ). The red stars and blue triangles in (a) indicate the locations of sources and receivers, respectively.
Figure 3. (a,b) are the true V P and V S model structures; (c) The red and black lines indicate the correct and wrong source wavelets ( s i c and s i w ). The red stars and blue triangles in (a) indicate the locations of sources and receivers, respectively.
Remotesensing 14 05223 g003
Figure 4. (a) Comparison of the traces (vertical component) of observed P-waves (black lines) and synthetic P-waves with the wrong source wavelet (red lines); The yellow and blue lines indicate their envelopes, respectively. (b) Comparison of the traces (vertical component) of observed S-waves (black lines) and synthetic S-waves with the wrong source wavelet (red lines); The yellow and blue lines are their envelopes.
Figure 4. (a) Comparison of the traces (vertical component) of observed P-waves (black lines) and synthetic P-waves with the wrong source wavelet (red lines); The yellow and blue lines indicate their envelopes, respectively. (b) Comparison of the traces (vertical component) of observed S-waves (black lines) and synthetic S-waves with the wrong source wavelet (red lines); The yellow and blue lines are their envelopes.
Remotesensing 14 05223 g004
Figure 5. (a,b) show the comparisons of absolute travel-time shifts for P- and S-waves, respectively. The blue lines indicate the travel-time shifts Δ t i c c , w obtained by cross-correlating the observed traces with synthetic traces using the wrong source wavelet. The yellow lines indicate the travel-time shifts Δ t i e t , w obtained by cross-correlating the envelopes of the observed traces with envelopes of the synthetic traces using the wrong source wavelet. The red bars indicate the travel-time shifts Δ t i c c , c obtained by cross-correlating the observed traces with the synthetic traces using the correct source wavelet.
Figure 5. (a,b) show the comparisons of absolute travel-time shifts for P- and S-waves, respectively. The blue lines indicate the travel-time shifts Δ t i c c , w obtained by cross-correlating the observed traces with synthetic traces using the wrong source wavelet. The yellow lines indicate the travel-time shifts Δ t i e t , w obtained by cross-correlating the envelopes of the observed traces with envelopes of the synthetic traces using the wrong source wavelet. The red bars indicate the travel-time shifts Δ t i c c , c obtained by cross-correlating the observed traces with the synthetic traces using the correct source wavelet.
Remotesensing 14 05223 g005
Figure 6. (a,b) are the inverted P-wave velocity model perturbations Δ V P c c , c and Δ V P c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the inverted model perturbations Δ V P e t , c and Δ V P e t , w by the ET objective function with correct and wrong source wavelets.
Figure 6. (a,b) are the inverted P-wave velocity model perturbations Δ V P c c , c and Δ V P c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the inverted model perturbations Δ V P e t , c and Δ V P e t , w by the ET objective function with correct and wrong source wavelets.
Remotesensing 14 05223 g006
Figure 7. (a,b) are the inverted S-wave velocity model perturbations Δ V S c c , c and Δ V S c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the inverted model perturbations Δ V S e t , c and Δ V S e t , w by ET objective function with correct and wrong source wavelets.
Figure 7. (a,b) are the inverted S-wave velocity model perturbations Δ V S c c , c and Δ V S c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the inverted model perturbations Δ V S e t , c and Δ V S e t , w by ET objective function with correct and wrong source wavelets.
Remotesensing 14 05223 g007
Figure 8. (a) Reductions of the normalized travel-time data misfits for V P inversion; (b) reductions of the normalized travel-time data misfits for V S inversion. The solid-black and dash-dotted-black lines indicate the traditional CC travel-time data misfits with correct and wrong source wavelets. The solid-red and dash-dotted-red lines indicate the new ET travel-time data misfits with correct and wrong source wavelets.
Figure 8. (a) Reductions of the normalized travel-time data misfits for V P inversion; (b) reductions of the normalized travel-time data misfits for V S inversion. The solid-black and dash-dotted-black lines indicate the traditional CC travel-time data misfits with correct and wrong source wavelets. The solid-red and dash-dotted-red lines indicate the new ET travel-time data misfits with correct and wrong source wavelets.
Remotesensing 14 05223 g008
Figure 9. (a,b) are the inverted perturbations of P- and S-wave velocities Δ V P e t , w and Δ V S e t , w b they ET objective function with wrong source wavelets when the velocity perturbations are ±30% of the background values.
Figure 9. (a,b) are the inverted perturbations of P- and S-wave velocities Δ V P e t , w and Δ V S e t , w b they ET objective function with wrong source wavelets when the velocity perturbations are ±30% of the background values.
Remotesensing 14 05223 g009
Figure 10. The initial V P model. The red stars and blue triangles indicate the schematic locations of sources and receivers.
Figure 10. The initial V P model. The red stars and blue triangles indicate the schematic locations of sources and receivers.
Remotesensing 14 05223 g010
Figure 11. (a) Raw shot gather of vertical component; (b) the shot gather after pre-processing.
Figure 11. (a) Raw shot gather of vertical component; (b) the shot gather after pre-processing.
Remotesensing 14 05223 g011
Figure 12. Comparison of the zero-phase Klauder source wavelet ( s i k , red line) with band-pass filtering of [10 Hz, 20 Hz] and wrong source wavelet ( s i w , black line).
Figure 12. Comparison of the zero-phase Klauder source wavelet ( s i k , red line) with band-pass filtering of [10 Hz, 20 Hz] and wrong source wavelet ( s i w , black line).
Remotesensing 14 05223 g012
Figure 13. (a,b) show the inverted velocity models V P c c , k and V P c c , w by the traditional CC objective function with the correct Klauder wavelet and wrong source wavelet, respectively; (c,d) are the inverted velocity models V P e t , k and V P e t , w by the ET objective function with the correct Klauder wavelet and wrong source wavelet.
Figure 13. (a,b) show the inverted velocity models V P c c , k and V P c c , w by the traditional CC objective function with the correct Klauder wavelet and wrong source wavelet, respectively; (c,d) are the inverted velocity models V P e t , k and V P e t , w by the ET objective function with the correct Klauder wavelet and wrong source wavelet.
Remotesensing 14 05223 g013
Figure 14. (a,b) show the inverted model perturbations Δ V P c c , k and Δ V P c c , w by the traditional CC objective function with the correct Klauder wavelet and wrong source wavelet; (c,d) are the inverted model perturbations Δ V P e t , k and Δ V P e t , w by the ET objective function with the correct Klauder wavelet and wrong source wavelet.
Figure 14. (a,b) show the inverted model perturbations Δ V P c c , k and Δ V P c c , w by the traditional CC objective function with the correct Klauder wavelet and wrong source wavelet; (c,d) are the inverted model perturbations Δ V P e t , k and Δ V P e t , w by the ET objective function with the correct Klauder wavelet and wrong source wavelet.
Remotesensing 14 05223 g014
Figure 15. Reductions of the normalized data misfits for travel-time inversion of early arrivals. The solid-black and dash-dotted-black lines indicate traditional CC objective function data misfits with correct and wrong source wavelets. The solid-red and dash-dotted-red lines indicate the ET objective function data misfits with correct and wrong source wavelets.
Figure 15. Reductions of the normalized data misfits for travel-time inversion of early arrivals. The solid-black and dash-dotted-black lines indicate traditional CC objective function data misfits with correct and wrong source wavelets. The solid-red and dash-dotted-red lines indicate the ET objective function data misfits with correct and wrong source wavelets.
Remotesensing 14 05223 g015
Figure 16. (a) Comparison of the observed data (black lines), synthetic data calculated with the initial model and wrong source (gray lines), and synthetic data calculated from the inverted model (Figure 13b) and wrong source (red lines); (b) comparison of the envelopes of observed data (black lines), synthetic data calculated with the initial model and wrong source (gray lines), and synthetic data calculated from the inverted model (Figure 13d) and wrong source (red lines).
Figure 16. (a) Comparison of the observed data (black lines), synthetic data calculated with the initial model and wrong source (gray lines), and synthetic data calculated from the inverted model (Figure 13b) and wrong source (red lines); (b) comparison of the envelopes of observed data (black lines), synthetic data calculated with the initial model and wrong source (gray lines), and synthetic data calculated from the inverted model (Figure 13d) and wrong source (red lines).
Remotesensing 14 05223 g016
Figure 17. Target velocity perturbation Δ V P .
Figure 17. Target velocity perturbation Δ V P .
Remotesensing 14 05223 g017
Figure 18. (a,b) are the recovered perturbations Δ V P c c , k and Δ V P c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the recovered perturbations Δ V P e t , k and Δ V P e t , w by the ET objective function with correct and wrong source wavelets.
Figure 18. (a,b) are the recovered perturbations Δ V P c c , k and Δ V P c c , w by the traditional CC objective function with correct and wrong source wavelets; (c,d) are the recovered perturbations Δ V P e t , k and Δ V P e t , w by the ET objective function with correct and wrong source wavelets.
Remotesensing 14 05223 g018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pan, W.; Ma, N.; Wang, Y. An Envelope Travel-Time Objective Function for Reducing Source–Velocity Trade-Offs in Wave-Equation Tomography. Remote Sens. 2022, 14, 5223. https://doi.org/10.3390/rs14205223

AMA Style

Pan W, Ma N, Wang Y. An Envelope Travel-Time Objective Function for Reducing Source–Velocity Trade-Offs in Wave-Equation Tomography. Remote Sensing. 2022; 14(20):5223. https://doi.org/10.3390/rs14205223

Chicago/Turabian Style

Pan, Wenyong, Ning Ma, and Yanfei Wang. 2022. "An Envelope Travel-Time Objective Function for Reducing Source–Velocity Trade-Offs in Wave-Equation Tomography" Remote Sensing 14, no. 20: 5223. https://doi.org/10.3390/rs14205223

APA Style

Pan, W., Ma, N., & Wang, Y. (2022). An Envelope Travel-Time Objective Function for Reducing Source–Velocity Trade-Offs in Wave-Equation Tomography. Remote Sensing, 14(20), 5223. https://doi.org/10.3390/rs14205223

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop