共查询到20条相似文献，搜索用时 46 毫秒
1.
2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D . 相似文献
2.
Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth. 相似文献
3.
Scattering theory, a form of perturbation theory, is a framework from within which time‐lapse seismic reflection methods can be derived and understood. It leads to expressions relating baseline and monitoring data and Earth properties, focusing on differences between these quantities as it does so. The baseline medium is, in the language of scattering theory, the reference medium and the monitoring medium is the perturbed medium. The general scattering relationship between monitoring data, baseline data, and time‐lapse Earth property changes is likely too complex to be tractable. However, there are special cases that can be analysed for physical insight. Two of these cases coincide with recognizable areas of applied reflection seismology: amplitude versus offset modelling/inversion, and imaging. The main result of this paper is a demonstration that time‐lapse difference amplitude versus offset modelling, and time‐lapse difference data imaging, emerge from a single theoretical framework. The time‐lapse amplitude versus offset case is considered first. We constrain the general time‐lapse scattering problem to correspond with a single immobile interface that separates a static overburden from a target medium whose properties undergo time‐lapse changes. The scattering solutions contain difference‐amplitude versus offset expressions that (although presently acoustic) resemble the expressions of Landro ( 2001 ). In addition, however, they contain non‐linear corrective terms whose importance becomes significant as the contrasts across the interface grow. The difference‐amplitude versus offset case is exemplified with two parameter acoustic (bulk modulus and density) and anacoustic (P‐wave velocity and quality factor Q) examples. The time‐lapse difference data imaging case is considered next. Instead of constraining the structure of the Earth volume as in the amplitude versus offset case, we instead make a small‐contrast assumption, namely that the time‐lapse variations are small enough that we may disregard contributions from beyond first order. An initial analysis, in which the case of a single mobile boundary is examined in 1D, justifies the use of a particular imaging algorithm applied directly to difference data shot records. This algorithm, a least‐squares, shot‐profile imaging method, is additionally capable of supporting a range of regularization techniques. Synthetic examples verify the applicability of linearized imaging methods of the difference image formation under ideal conditions. 相似文献
4.
Filters for migrated offset substacks are designed by partial coherence analysis to predict ‘normal’ amplitude variation with offset (AVO) in an anomaly free area. The same prediction filters generate localized prediction errors when applied in an AVO‐anomalous interval. These prediction errors are quantitatively related to the AVO gradient anomalies in a background that is related to the minimum AVO anomaly detectable from the data. The prediction‐error section is thus used to define a reliability threshold for the identification of AVO anomalies. Coherence analysis also enables quality control of AVO analysis and inversion. For example, predictions that are non‐localized and/or do not show structural conformity may indicate spatial variations in amplitude–offset scaling, seismic wavelet or signal‐to‐noise (S/N) ratio content. Scaling and waveform variations can be identified from inspection of the prediction filters and their frequency responses. S/N ratios can be estimated via multiple coherence analysis. AVO inversion of seismic data is unstable if not constrained. However, the use of a constraint on the estimated parameters has the undesirable effect of introducing biases into the inverted results: an additional bias‐correction step is then needed to retrieve unbiased results. An alternative form of AVO inversion that avoids additional corrections is proposed. This inversion is also fast as it inverts only AVO anomalies. A spectral coherence matching technique is employed to transform a zero‐offset extrapolation or near‐offset substack into P‐wave impedance. The same technique is applied to the prediction‐error section obtained by means of partial coherence, in order to estimate S‐wave velocity to P‐wave velocity (VS/VP) ratios. Both techniques assume that accurate well ties, reliable density measurements and P‐wave and S‐wave velocity logs are available, and that impedance contrasts are not too strong. A full Zoeppritz inversion is required when impedance contrasts that are too high are encountered. An added assumption is made for the inversion to the VS/VP ratio, i.e. the Gassmann fluid‐substitution theory is valid within the reservoir area. One synthetic example and one real North Sea in‐line survey illustrate the application of the two coherence methods. 相似文献
5.
Elastic full waveform inversion of seismic reflection data represents a data‐driven form of analysis leading to quantification of sub‐surface parameters in depth. In previous studies attention has been given to P‐wave data recorded in the marine environment, using either acoustic or elastic inversion schemes. In this paper we exploit both P‐waves and mode‐converted S‐waves in the marine environment in the inversion for both P‐ and S‐wave velocities by using wide‐angle, multi‐component, ocean‐bottom cable seismic data. An elastic waveform inversion scheme operating in the time domain was used, allowing accurate modelling of the full wavefield, including the elastic amplitude variation with offset response of reflected arrivals and mode‐converted events. A series of one‐ and two‐dimensional synthetic examples are presented, demonstrating the ability to invert for and thereby to quantify both P‐ and S‐wave velocities for different velocity models. In particular, for more realistic low velocity models, including a typically soft seabed, an effective strategy for inversion is proposed to exploit both P‐ and mode‐converted PS‐waves. Whilst P‐wave events are exploited for inversion for P‐wave velocity, examples show the contribution of both P‐ and PS‐waves to the successful recovery of S‐wave velocity. 相似文献
6.
—?The structural amplitude effect, associated with focusing and defocusing due to the reflector curvature, importantly contributes to reflection seismic amplitudes. This paper develops a conciliatory approach for estimating the structural amplitude effect and the attributes of amplitude variation versus offset (AVO). The AVO attributes are extracted from raw amplitudes, in which the structural effect is taken into account explicitly based on a structural model reconstructed from travel-time inversion. One of the goals is to conduct the AVO analysis not just locally (per CDP) but also horizontally to see the global variation along the reflection. The lateral variations of AVO attributes are decomposed by the Chebyshev expansion. The method is demonstrated with an example of weak shallow gas-water contact appearing on a 2-D seismic profile of a site survey in the North Sea. 相似文献
7.
Amir Abbas Babasafari Deva Prasad Ghosh Ahmed Mohammed Ahmed Salim Masoumeh Kordi 《Geophysical Prospecting》2020,68(8):2471-2493
Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification. 相似文献
8.
9.
Acoustic impedance is one of the best attributes for seismic interpretation and reservoir characterisation. We present an approach for estimating acoustic impedance accurately from a band‐limited and noisy seismic data. The approach is composed of two stages: inverting for reflectivity from seismic data and then estimating impedance from the reflectivity inverted in the first stage. For the first stage, we achieve a two‐step spectral inversion that locates the positions of reflection coefficients in the first step and determines the amplitudes of the reflection coefficients in the second step under the constraints of the positions located in the first step. For the second stage, we construct an iterative impedance estimation algorithm based on reflectivity. In each iteration, the iterative impedance estimation algorithm estimates the absolute acoustic impedance based on an initial acoustic impedance model that is given by summing the high‐frequency component of acoustic impedance estimated at the last iteration and a low‐frequency component determined in advance using other data. The known low‐frequency component is used to restrict the acoustic impedance variation tendency in each iteration. Examples using one‐ and two‐dimensional synthetic and field seismic data show that the approach is flexible and superior to the conventional spectral inversion and recursive inversion methods for generating more accurate acoustic impedance models. 相似文献
10.
Seismic velocity analysis in the scattering-angle/azimuth domain 总被引：2，自引：0，他引：2
Migration velocity analysis is carried out by analysing the residual moveout and amplitude variations in common image point gathers (CIGs) parametrized by scattering angle and azimuth. The misfit criterion in the analysis is of the differential-semblance type. By using angles to parametrize the imaging we are able to handle and exploit data with multiple arrivals, although artefacts may occur in the CIGs and need to be suppressed. The CIGs are generated by angle migration, an approach based on the generalized Radon transform (GRT) inversion, and they provide multiple images of reflectors in the subsurface for a range of scattering angles and azimuths. Within the differential semblance applied to these CIGs, we compensate for amplitude versus angle (AVA) effects. Thus, using a correct background velocity model, the CIGs should have no residual moveout nor amplitude variation with angles, and the differential semblance should vanish. If the velocity model is incorrect, however, the events in the CIGs will appear at different depths for different angles and the amplitude along the events will be non-uniform. A standard, gradient-based optimization scheme is employed to develop a velocity updating procedure. The model update is formed by backprojecting the differential semblance misfits through ray perturbation kernels, within a GRT inverse. The GRT inverse acts on the data, subject to a shift in accordance with ray perturbation theory. The performance of our algorithm is demonstrated with two synthetic data examples using isotropic elastic models. The first one allows velocity variation with depth only. In the second one, we reconstruct a low-velocity lens in the model that gives rise to multipathing. The velocity model parametrization is based upon the eigentensor decomposition of the stiffness tensor and makes use of B-splines. 相似文献
11.
Yanan Liu Shengchang Chen Guoxin Chen Feng Wang Zhiwei Gu 《Geophysical Prospecting》2020,68(8):2504-2517
Impedance is a physical parameter that plays an important role in seismic data processing and interpretation. A relative impedance perturbation (the ratio of the impedance perturbation and the impedance for the background models) imaging method in depth domain based on the reflection wave equation is proposed. Under the small perturbation assumption, primary wave and high-frequency approximation condition, a linear propagation equation of the primary reflection waves based on the relative impedance perturbation was first derived. On this basis, we further derived the imaging formula of the relative impedance perturbation using a linear inversion theory. Then, the source–receiver bidirectional illumination compensation was used to improve the image quality of the subsurface structures. The image result obtained by this method can be used to estimate the relative impedance perturbation. In the angle domain, the extracted near-angle-domain image gather with amplitude compensation can estimate the relative impedance perturbation, and the far-angle image gather provides the estimation of the relative velocity perturbation (the ratio of the velocity perturbation and the background velocity). Finally, several numerical tests demonstrate the effectiveness of the method. 相似文献
12.
In seismic waveform inversion, non‐linearity and non‐uniqueness require appropriate strategies. We formulate four types of L2 normed misfit functionals for Laplace‐Fourier domain waveform inversion: i) subtraction of complex‐valued observed data from complex‐valued predicted data (the ‘conventional phase‐amplitude’ residual), ii) a ‘conventional phase‐only’ residual in which amplitude variations are normalized, iii) a ‘logarithmic phase‐amplitude’ residual and finally iv) a ‘logarithmic phase‐only’ residual in which the only imaginary part of the logarithmic residual is used. We evaluate these misfit functionals by using a wide‐angle field Ocean Bottom Seismograph (OBS) data set with a maximum offset of 55 km. The conventional phase‐amplitude approach is restricted in illumination and delineates only shallow velocity structures. In contrast, the other three misfit functionals retrieve detailed velocity structures with clear lithological boundaries down to the deeper part of the model. We also test the performance of additional phase‐amplitude inversions starting from the logarithmic phase‐only inversion result. The resulting velocity updates are prominent only in the high‐wavenumber components, sharpening the lithological boundaries. We argue that the discrepancies in the behaviours of the misfit functionals are primarily caused by the sensitivities of the model gradient to strong amplitude variations in the data. As the observed data amplitudes are dominated by the near‐offset traces, the conventional phase‐amplitude inversion primarily updates the shallow structures as a result. In contrast, the other three misfit functionals eliminate the strong dependence on amplitude variation naturally and enhance the depth of illumination. We further suggest that the phase‐only inversions are sufficient to obtain robust and reliable velocity structures and the amplitude information is of secondary importance in constraining subsurface velocity models. 相似文献
13.
Guido Baeten Jan Willem de Maag René‐Edouard Plessix Rini Klaassen Tahira Qureshi Maren Kleemeyer Fons ten Kroode Zhang Rujie 《Geophysical Prospecting》2013,61(4):701-711
Velocity model building and impedance inversion generally suffer from a lack of intermediate wavenumber content in seismic data. Intermediate wavenumbers may be retrieved directly from seismic data sets if enough low frequencies are recorded. Over the past years, improvements in acquisition have allowed us to obtain seismic data with a broader frequency spectrum. To illustrate the benefits of broadband acquisition, notably the recording of low frequencies, we discuss the inversion of land seismic data acquired in Inner Mongolia, China. This data set contains frequencies from 1.5–80 Hz. We show that the velocity estimate based on an acoustic full‐waveform inversion approach is superior to one obtained from reflection traveltime inversion because after full‐waveform inversion the background velocity conforms to geology. We also illustrate the added value of low frequencies in an impedance estimate. 相似文献
14.
Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media 
下载免费PDF全文
下载免费PDF全文 Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process. 相似文献
15.
Sensitivity analysis and application of time‐lapse full‐waveform inversion: synthetic testing and field data example from the North Sea,Norway 下载免费PDF全文
下载免费PDF全文 Time‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (?30 m/s, or ?1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep). 相似文献
16.
Amplitude variation with offset (AVO) analysis and waveform inversion are techniques used to determine qualitative or quantitative information on gas hydrates and free gas in sediments. However, the quantitative contribution of gas hydrates to the acoustic impedance contrast observed at the bottom‐simulating reflector and the reliability of quantitative AVO analyses are still topics of discussion. In this study, common‐midpoint gathers from multichannel wide‐angle reflection seismic data, acquired offshore Costa Rica, have been processed to preserve true amplitude information at the bottom‐simulating reflector for a quantitative AVO analysis incorporating angles of incidence of up to 60°. Corrections were applied for effects that significantly alter the observed amplitudes, such as the source directivity. AVO and rock‐physics modelling indicate that free gas immediately beneath the gas‐hydrate stability zone can be detected and low concentrations can be quantified from AVO analysis, whereas the offset‐dependent reflectivity is not sensitive to gas‐hydrate concentrations of less than about 10% at the base of the gas‐hydrate stability zone. Bulk free‐gas saturations up to 5% have been determined from the reflection seismic data assuming a homogeneous distribution of free gas in the sediment. Assuming a patchy distribution of free gas increases the estimated concentrations up to 14%. There is a patchy occurrence of bottom‐simulating reflectors south‐east of the Nicoya Peninsula on the continental margin, offshore Costa Rica. AVO analysis indicates that this phenomenon is related to the local presence of free gas beneath the gas‐hydrate stability zone, probably related to a focused vertical fluid flow. In areas without bottom‐simulating reflectors, the results indicate that no free gas is present. 相似文献
17.
开发导致的时移地震差异代表了油藏的岩石物理性质变化,本文从差异波形角度出发,明确了声波阻抗界面的反射系数与油层声波阻抗变化之间的函数关系,进而分析薄油层时移地震差异波形特征.通过分析得出,开发前后油层阻抗变化越大,时移地震的差异振幅幅度越大;薄油层的差异波形特征仅与油层的声波阻抗变化有关,油层声波阻抗增大导致与地震子波微分形式相同的差异波形,油层声波阻抗减小产生的差异波形与180°相位的地震子波微分形式相同;油层声波阻抗变化方向导致差异波形振幅幅度变化的快慢. 相似文献
18.
常规的地震道反演方法建立在反射P波垂直入射假设 的基础上,而实际地震资料采集时多数是非零炮检距的，反射振幅是共中心点道集叠加的结 果 . 因此，利用常规地震道反演方法就不能得到可靠的波阻抗或其他岩性信息. 本文利用Patr ick Connolly弹性阻抗的思想，通过对Zoeppritz方程的进一步简化，推导出适合常规叠后 资料的、非零炮检距条件下纵波反射系数递推公式，提出了广义弹性阻抗的概念，解决了非 零炮检距条件下，常规叠后地震道正反演的关键问题. 广义弹性阻抗不仅包含波阻抗，还包 含了纵横波速度等岩性信息，具有很好的实用价值. 进行广义弹性阻抗的反演，能较常规地 震道反演获得更多、更可靠的流体、孔隙度、砂泥含量等信息，有助于解释常规地震道反演 和道积分剖面中的假象，降低反演的多解性，提高储层预测的精度. 相似文献
19.
《中国科学:地球科学(英文版)》2015,(8)
The existing expressions of elastic impedance,as the generalized form of acoustic impedance,represent the resistance of subsurface media to seismic waves of non-normal incidence,and thus include information on the shear-wave velocity.In this sense,conventional elastic impedance is an attribute of the seismic reflection and not an intrinsic physical property of the subsurface media.The derivation of these expressions shares the approximations made for reflectivity,such as weak impedance contrast andisotropic or weakly anisotropic media,which limits the accuracy of reflectivity reconstruction and seismic inversion.In this paper,we derive exact elastic impedance tensors of seismic P-and S-waves for isotropic media based on the stress-velocity law.Each componentof the impedance tensor represents a unique mechanical property of the medium.Approximations of P-wave elastic impedance tensor components are discussed for seismic inversion and interpretation.Application to synthetic data and real data shows the accuracy and robust interpretation capability of the derived elastic impedance in lithology characterizations. 相似文献
20.
Tariq Alkhalifah 《Geophysical Prospecting》2016,64(2):505-513
While velocity contrasts are responsible for most of the events recorded in our data, the long wavelength behavior of the velocity model is responsible for the geometrical shape of these events. For isotropic acoustic materials, the wave dependency on the long (wave propagation) and short (scattering) wavelength velocity components is stationary with the propagation angle. On the other hand, in representing a transversely isotropic with a vertical symmetry axis medium with the normal moveout velocity, the anellepticity parameter η, the vertical scaling parameter δ, and the sensitivity of waves vary with the polar angle for both the long and short wavelength features of the anisotropic dimensionless medium parameters (δ and η). For horizontal reflectors at reasonable depths, the long wavelength features of the η model is reasonably constrained by the long offsets, whereas the short wavelength features produce very week reflections at even reasonable offsets. Thus, for surface acquired seismic data, we could mainly invert for smooth η responsible for the geometrical shape of reflections. On the other hand, while the δ long wavelength components mildly affects the recorded data, its short wavelength variations can produce reflections at even zero offset, with a behavior pattern synonymous to density. The lack of the long wavelength δ information will mildly effect focusing but will cause misplacement of events in depth. With low enough frequencies (very low), we may be able to recover the long wavelength δ using full waveform inversion. However, unlike velocity, the frequencies needed for that should be ultra‐low to produce long‐wavelength scattering‐based model information as δ perturbations do not exert scattering at large offsets. For a combination given by the horizontal velocity, η, and ε, the diving wave influence of η is absorbed by the horizontal velocity, severely limiting the η influence on the data and full waveform inversion. As a result, with a good smooth η estimation, for example, from tomography, we can focus the full waveform inversion to invert for only the horizontal velocity and maybe ε as a parameter to fit the amplitude. This is possibly the most practical parametrization for inversion of surface seismic data in transversely isotropic with vertical symmetry axis media. 相似文献