共查询到20条相似文献，搜索用时 98 毫秒
1.
Stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter χ eff are correlated in the PS-converted-wave (PS-wave) anisotropic prestack Kirchhoff time migration (PKTM) velocity model and are thus difficult to independently determine. We extended the simplified two-parameter (stacking velocity V C2 and anisotropic parameter k eff) moveout equation from stacking velocity analysis to PKTM velocity model updating and formed a new four-parameter (stacking velocity V C2, vertical velocity ratio γ 0, effective velocity ratio γ eff, and anisotropic parameter k eff) PS-wave anisotropic PKTM velocity model updating and process flow based on the simplified two-parameter moveout equation. In the proposed method, first, the PS-wave two-parameter stacking velocity is analyzed to obtain the anisotropic PKTM initial velocity and anisotropic parameters; then, the velocity and anisotropic parameters are corrected by analyzing the residual moveout on common imaging point gathers after prestack time migration. The vertical velocity ratio γ 0 of the prestack time migration velocity model is obtained with an appropriate method utilizing the P- and PS-wave stacked sections after level calibration. The initial effective velocity ratio γ eff is calculated using the Thomsen (1999) equation in combination with the P-wave velocity analysis; ultimately, the final velocity model of the effective velocity ratio γ eff is obtained by percentage scanning migration. This method simplifies the PS-wave parameter estimation in high-quality imaging, reduces the uncertainty of multiparameter estimations, and obtains good imaging results in practice. 相似文献
2.
In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far
offsets is determined by four parameters. These are the C-wave stacking velocity V
C2, the vertical and effective velocity ratios γ
0and γ
eff, and the anisotropic parameter X
eff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine
this stacking velocity model.
The C-wave stacking velocity model V
C2, γ
0, γ
geff, and X
eff can be determined from P- and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout
inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy
to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z
from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield V
C2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes
it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose:
the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield
results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow.
This work is funded by the Edinburgh Anisotropy Project of the British Geological Survey.
First Author
Li Xiangyang, he is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project
in the British Geological Survey. He also holds a honorary professorship multicomponent seismology at the School of Geosciences,
University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984)
in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992)
in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute.
Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent
seismology. 相似文献
3.
在具有垂直对称轴横向各向同性介质中，利用四种参数来确定中间至远偏移距转换波（C-波）动校正。它们是C-波叠加速度Vc2，垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2，γ0，γeff，和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得，因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略，而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2，若数据不含燥音，其误差小于1％。方法稳健，允许γ启始假定值的误差达20％。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程：双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似，但是单扫描流程比双扫描更有效。 相似文献
4.
5.
The azimuthally varying non‐hyperbolic moveout of P‐waves in orthorhombic media can provide valuable information for characterization of fractured reservoirs and seismic processing. Here, we present a technique to invert long‐spread, wide‐azimuth P‐wave data for the orientation of the vertical symmetry planes and five key moveout parameters: the symmetry‐plane NMO velocities, V(1)nmo and V(2)nmo , and the anellipticity parameters, η(1), η(2) and η(3) . The inversion algorithm is based on a coherence operator that computes the semblance for the full range of offsets and azimuths using a generalized version of the Alkhalifah–Tsvankin non‐hyperbolic moveout equation. The moveout equation provides a close approximation to the reflection traveltimes in layered anisotropic media with a uniform orientation of the vertical symmetry planes. Numerical tests on noise‐contaminated data for a single orthorhombic layer show that the best‐constrained parameters are the azimuth ? of one of the symmetry planes and the velocities V(1)nmo and V(2)nmo , while the resolution in η(1) and η(2) is somewhat compromised by the trade‐off between the quadratic and quartic moveout terms. The largest uncertainty is observed in the parameter η(3) , which influences only long‐spread moveout in off‐symmetry directions. For stratified orthorhombic models with depth‐dependent symmetry‐plane azimuths, the moveout equation has to be modified by allowing the orientation of the effective NMO ellipse to differ from the principal azimuthal direction of the effective quartic moveout term. The algorithm was successfully tested on wide‐azimuth P‐wave reflections recorded at the Weyburn Field in Canada. Taking azimuthal anisotropy into account increased the semblance values for most long‐offset reflection events in the overburden, which indicates that fracturing is not limited to the reservoir level. The inverted symmetry‐plane directions are close to the azimuths of the off‐trend fracture sets determined from borehole data and shear‐wave splitting analysis. The effective moveout parameters estimated by our algorithm provide input for P‐wave time imaging and geometrical‐spreading correction in layered orthorhombic media. 相似文献
6.
Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t2(x2) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (VP0 and VS0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (VP0, VS0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves. 相似文献
7.
Elastic imaging from ocean bottom cable (OBC) data can be challenging because it requires the prior estimation of both compressional‐wave (P‐wave) and shear‐wave (S‐wave) velocity fields. Seismic interferometry is an attractive technique for processing OBC data because it performs model‐independent redatuming; retrieving ‘pseudo‐sources’ at positions of the receivers. The purpose of this study is to investigate multicomponent applications of interferometry for processing OBC data. This translates into using interferometry to retrieve pseudo‐source data on the sea‐bed not only for multiple suppression but for obtaining P‐, converted P to S‐wave (PS‐wave) and possibly pure mode S‐waves. We discuss scattering‐based, elastic interferometry with synthetic and field OBC datasets. Conventional and scattering‐based interferometry integrands computed from a synthetic are compared to show that the latter yields little anti‐causal response. A four‐component (4C) pseudo‐source response retrieves pure‐mode S‐reflections as well at P‐ and PS‐reflections. Pseudo‐source responses observed in OBC data are related to P‐wave conversions at the seabed rather than to true horizontal or vertical point forces. From a Gulf of Mexico OBC data set, diagonal components from a nine‐component pseudo‐source response demonstrate that the P‐wave to S‐wave velocity ratio (VP/VS) at the sea‐bed is an important factor in the conversion of P to S for obtaining the pure‐mode S‐wave reflections. 相似文献
8.
A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity VP0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces. 相似文献
9.
Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey 总被引：2，自引：0，他引：2
Roger A. Clark Philip M. Benson rew J. Carter Carlos A. Guerrero Moreno 《Geophysical Prospecting》2009,57(5):835-845
A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q−1eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q−1eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q−1eff= 52 ± 15) and slowest stacking velocity. The variation of Q−1eff with azimuth φ is well fitted by 1000Q−1eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 114 combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10−5, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here. 相似文献
10.
A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks 总被引：4，自引：0，他引：4
The laboratory ultrasonic pulse‐echo method was used to collect accurate P‐ and S‐wave velocity (±0.3%) and attenuation (±10%) data at differential pressures of 5–50 MPa on water‐saturated core samples of sandstone, limestone and siltstone that were cut parallel and perpendicular to the vertical borehole axis. The results, when expressed in terms of the P‐ and S‐wave velocity and attenuation anisotropy parameters for weakly transversely isotropic media (ɛ, γ, ɛQ, γQ) show complex variations with pressure and lithology. In general, attenuation anisotropy is stronger and more sensitive to pressure changes than velocity anisotropy, regardless of lithology. Anisotropy is greatest (over 20% for velocity, over 70% for attenuation) in rocks with visible clay/organic matter laminations in hand specimens. Pressure sensitivities are attributed to the opening of microcracks with decreasing pressure. Changes in magnitude of velocity and attenuation anisotropy with effective pressure show similar trends, although they can show different signs (positive or negative values of ɛ, ɛQ, γ, γQ). We conclude that attenuation anisotropy in particular could prove useful to seismic monitoring of reservoir pressure changes if frequency‐dependent effects can be quantified and modelled. 相似文献
11.
Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold‐and‐thrust belts) and in subsalt exploration. Here, we introduce a methodology for P‐wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters ε and δ and linearly varying symmetry‐direction velocity VP0 controlled by the vertical (kz) and lateral (kx) gradients. Since determination of tilt from P‐wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity VP0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients kz and kx and the anisotropy parameters ε and δ in the layer‐stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry‐axis direction is fixed and VP0 is known, the parameters kz, kx, ε and δ can be resolved from reflection data. It should be emphasized that estimation of ε in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas. 相似文献
12.
Anisotropy is often observed due to the thin layering or aligned micro‐structures, like small fractures. At the scale of cross‐well tomography, the anisotropic effects cannot be neglected. In this paper, we propose a method of full‐wave inversion for transversely isotropic media and we test its robustness against structured noisy data. Optimization inversion techniques based on a least‐square formalism are used. In this framework, analytical expressions of the misfit function gradient, based on the adjoint technique in the time domain, allow one to solve the inverse problem with a high number of parameters and for a completely heterogeneous medium. The wave propagation equation for transversely isotropic media with vertical symmetry axis is solved using the finite difference method on the cylindrical system of coordinates. This system allows one to model the 3D propagation in a 2D medium with a revolution symmetry. In case of approximately horizontal layering, this approximation is sufficient. The full‐wave inversion method is applied to a crosswell synthetic 2‐component (radial and vertical) dataset generated using a 2D model with three different anisotropic regions. Complex noise has been added to these synthetic observed data. This noise is Gaussian and has the same amplitude f?k spectrum as the data. Part of the noise is localized as a coda of arrivals, the other part is not localized. Five parameter fields are estimated, (vertical) P‐wave velocity, (vertical) S‐wave velocity, volumetric mass and the Thomsen anisotropic parameters epsilon and delta. Horizontal exponential correlations have been used. The results show that the full‐wave inversion of cross‐well data is relatively robust for high‐level noise even for second‐order parameters such as Thomsen epsilon and delta anisotropic parameters. 相似文献
13.
Several parameters are needed to describe the converted-wave （C-wave） moveout in processing multi-component seismic data, because of asymmetric raypaths and anisotropy. As the number of parameters increases, the converted wave data processing and analysis becomes more complex. This paper develops a new moveout equation with two parameters for C-waves in vertical transverse isotropy （VTI） media. The two parameters are the C-wave stacking velocity （Vc2） and the squared velocity ratio （7v,i） between the horizontal P-wave velocity and C-wave stacking velocity. The new equation has fewer parameters, but retains the same applicability as previous ones. The applicability of the new equation and the accuracy of the parameter estimation are checked using model and real data. The form of the new equation is the same as that for layered isotropic media. The new equation can simplify the procedure for C-wave processing and parameter estimation in VTI media, and can be applied to real C-wave processing and interpretation. Accurate Vc2 and Yvti can be deduced from C-wave data alone using the double-scanning method, and the velocity ratio model is suitable for event matching between P- and C-wave data. 相似文献
14.
Although it is widely recognized that anisotropy can have a significant influence on the focusing and positioning of migrated reflection events, conventional depth imaging methods still operate with isotropic velocity fields. Here, we present an application of a 2D migration velocity analysis (MVA) algorithm, designed for factorized v(x, z) VTI (transversely isotropic with a vertical symmetry axis) media, to an offshore data set from West Africa. By approximating the subsurface with factorized VTI blocks, it is possible to decouple the spatial variations in the vertical velocity from the anisotropic parameters with minimal a priori information. Since our method accounts for lateral velocity variation, it produces more accurate estimates of the anisotropic parameters than those previously obtained with time‐domain techniques. The values of the anellipticity parameter η found for the massive shales exceed 0.2, which confirms that ignoring anisotropy in the study area can lead to substantial imaging distortions, such as mis‐stacking and mispositioning of dipping events. While some of these distortions can be removed by using anisotropic time processing, further marked improvement in image quality is achieved by prestack depth migration with the estimated factorized VTI model. In particular, many fault planes, including antithetic faults in the shallow part of the section, are better focused by the anisotropic depth‐migration algorithm and appear more continuous. Anisotropic depth migration facilitates structural interpretation by eliminating false dips at the bottom of the section and improving the images of a number of gently dipping features. One of the main difficulties in anisotropic MVA is the need to use a priori information for constraining the vertical velocity. In this case study, we successfully reconstructed the time–depth curve from reflection data by assuming that the vertical velocity is a continuous function of depth and estimating the vertical and lateral velocity gradients in each factorized block. If the subsurface contains strong boundaries with jumps in velocity, knowledge of the vertical velocity at a single point in a layer is sufficient for our algorithm to determine all relevant layer parameters. 相似文献
15.
Tariq Alkhalifah 《Geophysical Prospecting》2004,52(1):39-48
The azimuth moveout (AMO) operator in homogeneous transversely isotropic media with a vertical symmetry axis (VTI), as in isotropic media, has an overall skewed saddle shape. However, the AMO operator in anisotropic media is complicated; it includes, among other things, triplications at low angles. Even in weaker anisotropies, with the anisotropy parameter η= 0.1 (10% anisotropy), the AMO operator is considerably different from the isotropic operator, although free of triplications. The structure of the operator in VTI media (positive η) is stretched (has a wider aperture) compared with operators in isotropic media, with the amount of stretch being dependent on the strength of anisotropy. If the medium is both vertically inhomogeneous, i.e. the vertical velocity is a function of depth (v(z)), and anisotropic, which is a common combination in practical problems, the shape of the operator again differs from that for isotropic media. However, the difference in the AMO operator between the homogeneous and the v(z) cases, even for anisotropic media, is small. Stated simply, anisotropy influences the shape and aperture of the AMO operator far more than vertical inhomogeneity does. 相似文献
16.
A practical method is presented for determining three‐dimensional S‐wave velocity (VS) profile from microtremor measurements. Frequency–wave number (f–k) spectral analyses of microtremor array records are combined, for this purpose, with microtremor horizontal‐to‐vertical (H/V) spectral ratio techniques. To demonstrate the effectiveness of the proposed method, microtremor measurements using arrays of sensors were conducted at six sites in the city of Kushiro, Japan. The spectral analyses of the array records yield dispersion characteristics of Rayleigh waves and H/V spectra of surface waves, and joint inversion of these data results in VS profiles down to bedrock at the sites. Conventional microtremor measurements were performed at 230 stations within Kushiro city, resulting in the H/V spectra within the city. Three‐dimensional VS structure is then estimated from inversion of the H/V spectra with the VS values determined from the microtremor array data. This reveals three‐dimensional VS profile of Kushiro city, together with an unknown hidden valley that crosses the central part of the city. The estimated VS profile is consistent with available velocity logs and results of subsequent borings, indicating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
17.
A multi‐azimuth seismic refraction study for a horizontal transverse isotropic medium: physical modelling results 
下载免费PDF全文
下载免费PDF全文 To investigate the characteristics of the anisotropic stratum, a multi‐azimuth seismic refraction technique is proposed in this study since the travel time anomaly of the refraction wave induced by this anisotropic stratum will be large for a far offset receiver. To simplify the problem, a two‐layer (isotropy–horizontal transverse isotropy) model is considered. A new travel time equation of the refracted P‐wave propagation in this two‐layer model is derived, which is the function of the phase and group velocities of the horizontal transverse isotropic stratum. In addition, the measured refraction wave velocity in the physical model experiment is the group velocity. The isotropic intercept time equation of a refraction wave can be directly used to estimate the thickness of the top (isotropic) layer of the two‐layer model because the contrast between the phase and group velocities of the horizontal transverse isotropic medium is seldom greater than 10% in the Earth. If the contrast between the phase and group velocities of an anisotropic medium is small, the approximated travel time equation of a refraction wave is obtained. This equation is only dependent on the group velocity of the horizontal transverse isotropic stratum. The elastic constants A11, A13, and A33 and the Thomsen anisotropic parameter ε of the horizontal transverse isotropic stratum can be estimated using this multi‐azimuth seismic refraction technique. Furthermore, under a condition of weak anisotropy, the Thomsen anisotropic parameter δ of the horizontal transverse isotropic stratum can be estimated by this technique as well. 相似文献
18.
A nonlinear method for multiparameter inversion of pre‐stack seismic data based on anisotropic Markov random field 下载免费PDF全文
下载免费PDF全文 Hongbing Zhang Qiang Guo Lifeng Liang Chenghao Cao Zuoping Shang 《Geophysical Prospecting》2018,66(3):461-477
Multiparameter inversion for pre‐stack seismic data plays a significant role in quantitative estimation of subsurface petrophysical properties. However, it remains a complicated problem due to the non‐unique results and unstable nature of the processing; the pre‐stack seismic inversion problem is ill‐posed and band‐limited. Combining the full Zoeppritz equation and additional assumptions with edge‐preserving regularisation can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function that includes edge‐preserving regularisation and soft constraints based on anisotropic Markov random fields and is intended especially for layered formations. We applied a fast simulated annealing algorithm to solve the nonlinear optimisation problem. The method directly obtains reflectivity RPP values using the full Zoeppritz equation instead of its approximations and effectively controls the stability of the multiparameter inversion by assuming a sectionally constant S‐ and P‐wave velocity ratio and using the generalised Gardner equation. We substituted the inverted parameters, i.e., the P‐wave velocity, the fitting deviation of S‐wave velocity, and the density were inverted instead of the P‐wave velocity, the S‐wave velocity, and the density, and the generalised Gardner equation was applied as a constraint. Test results on two‐dimensional synthetic data indicated that our substitution obtained improved results for multiparameter inversion. The inverted results could be improved by utilising high‐order anisotropic Markov random field neighbourhoods at early stages and low‐order anisotropic Markov random field neighbourhoods in the later stages. Moreover, for layered formations, using a large horizontal weighting coefficient can preserve the lateral continuity of layers, and using a small vertical weighting coefficient allows for large longitudinal gradients of the interlayers. The inverted results of the field data revealed more detailed information about the layers and matched the logging curves at the wells acceptably over most parts of the curves. 相似文献
19.
Simultaneous estimation of velocity gradients and anisotropic parameters from seismic reflection data is one of the main challenges in transversely isotropic media with a vertical symmetry axis migration velocity analysis. In migration velocity analysis, we usually construct the objective function using the l2 norm along with a linear conjugate gradient scheme to solve the inversion problem. Nevertheless, for seismic data this inversion scheme is not stable and may not converge in finite time. In order to ensure the uniform convergence of parameter inversion and improve the efficiency of migration velocity analysis, this paper develops a double parameterized regularization model and gives the corresponding algorithms. The model is based on the combination of the l2 norm and the non‐smooth l1 norm. For solving such an inversion problem, the quasi‐Newton method is utilized to make the iterative process stable, which can ensure the positive definiteness of the Hessian matrix. Numerical simulation indicates that this method allows fast convergence to the true model and simultaneously generates inversion results with a higher accuracy. Therefore, our proposed method is very promising for practical migration velocity analysis in anisotropic media. 相似文献
20.
土体剪切波速是进行土层地震反应分析的动力学参数，对场地地震动参数确定具有重要意义。基于地质地貌分析，将大同盆地划分为5类典型地质单元。对盆地1429个钻孔剪切波速资料进行分析，探讨VS30与VS20的相关性，研究土体埋深、岩性、地质单元、标贯击数及密实度等地质特征对VS的影响，并基于地质单元、剪切波速比、密实度系数及第四系上部覆盖层厚度相关性分析给出土体VS30预测模型。研究结果表明，基于典型地质特征的VS30预测模型拟合优度R2>0.90，预测精度很高，对于离散性较大、直接拟合估算较差及无剪切波速场地来说，以区分地质单元及土体类型的方式进行VS30分解预测是良好的研究思路。首次在区分地质单元及土体类型的前提下提出剪切波速比及密实度系数，并将其与第四系上部覆盖层厚度综合应用于VS30预测研究。研究结果可为大同盆地城市防震减灾规划、震害预测、区域性地震安全评价提供重要技术支撑。 相似文献