共查询到20条相似文献，搜索用时 46 毫秒
1.
2.
Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain. 相似文献
3.
Myung Hoon Kim Yunseok Choi Young Ho Cha Changsoo Shin 《Pure and Applied Geophysics》2009,166(12):1967-1985
In order to correctly interpret marine exploration data, which contain many elastic signals such as S waves, surface waves and converted waves, we have developed both a frequency-domain modeling algorithm for acoustic-elastic coupled media with an irregular interface, and the corresponding waveform inversion algorithm. By applying the continuity condition between acoustic (fluid) and elastic (solid) media, wave propagation can be properly simulated throughout the coupled domain. The arbitrary interface is represented by tessellating square and triangular finite elements. Although the resulting complex impedance matrix generated by finite element methods for the acoustic-elastic coupled wave equation is asymmetric, we can exploit the usual back-propagation algorithm used in the frequency domain through modern sparse matrix technology. By running numerical experiments on a synthetic model, we demonstrate that our inversion algorithm can successfully recover P- and S-wave velocity and density models from marine exploration data (pressure data only). 相似文献
4.
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. 相似文献
5.
Ho Seuk Bae Changsoo Shin Young Ho Cha Yunseok Choi Dong‐Joo Min 《Geophysical Prospecting》2010,58(6):997-1010
Although waveform inversion has been intensively studied in an effort to properly delineate the Earth's structures since the early 1980s, most of the time‐ and frequency‐domain waveform inversion algorithms still have critical limitations in their applications to field data. This may be attributed to the highly non‐linear objective function and the unreliable low‐frequency components. To overcome the weaknesses of conventional waveform inversion algorithms, the acoustic Laplace‐domain waveform inversion has been proposed. The Laplace‐domain waveform inversion has been known to provide a long‐wavelength velocity model even for field data, which may be because it employs the zero‐frequency component of the damped wavefield and a well‐behaved logarithmic objective function. However, its applications have been confined to 2D acoustic media. We extend the Laplace‐domain waveform inversion algorithm to a 2D acoustic‐elastic coupled medium, which is encountered in marine exploration environments. In 2D acoustic‐elastic coupled media, the Laplace‐domain pressures behave differently from those of 2D acoustic media, although the overall features are similar to each other. The main differences are that the pressure wavefields for acoustic‐elastic coupled media show negative values even for simple geological structures unlike in acoustic media, when the Laplace damping constant is small and the water depth is shallow. The negative values may result from more complicated wave propagation in elastic media and at fluid‐solid interfaces. Our Laplace‐domain waveform inversion algorithm is also based on the finite‐element method and logarithmic wavefields. To compute gradient direction, we apply the back‐propagation technique. Under the assumption that density is fixed, P‐ and S‐wave velocity models are inverted from the pressure data. We applied our inversion algorithm to the SEG/EAGE salt model and the numerical results showed that the Laplace‐domain waveform inversion successfully recovers the long‐wavelength structures of the P‐ and S‐wave velocity models from the noise‐free data. The models inverted by the Laplace‐domain waveform inversion were able to be successfully used as initial models in the subsequent frequency‐domain waveform inversion, which is performed to describe the short‐wavelength structures of the true models. 相似文献
6.
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. 相似文献
7.
Alexey Stovas 《Geophysical Prospecting》2017,65(2):426-452
I derive the kinematic properties of single‐mode P, S1, and S2 waves as well as converted PS1, PS2, and S1S2 waves in elastic orthorhombic media including vertical velocity, two normal moveout velocities defined in vertical symmetry planes, and three anelliptic parameters (two of them are defined in vertical symmetry plane and one parameter is the cross‐term one). I show that the azimuthal dependence of normal moveout velocity and anellipticity is different in phase and group domains. The effects on‐vertical‐axis singularity and on‐vertical‐axis triplication are considered for pure‐mode S1 and S2 waves and converted‐mode S1S2 waves. The conditions and properties of on‐vertical‐axis triplication are defined in terms of kinematic parameters. The results are illustrated in four homogeneous orthorhombic models and one multilayered orthorhombic model with no variation in azimuthal orientation for all the layers. 相似文献
8.
9.
Anisotropic elastic full‐waveform inversion of walkaway vertical seismic profiling data from the Arabian Gulf 
下载免费PDF全文
下载免费PDF全文 John C Owusu Marwan Charara Scott Leaney Allan Campbell Shujaat Ali Igor Borodin Les Nutt Henry Menkiti 《Geophysical Prospecting》2016,64(1):38-53
Borehole seismic addresses the need for high‐resolution images and elastic parameters of the subsurface. Full‐waveform inversion of vertical seismic profile data is a promising technology with the potential to recover quantitative information about elastic properties of the medium. Full‐waveform inversion has the capability to process the entire wavefield and to address the wave propagation effects contained in the borehole data—multi‐component measurements; anisotropic effects; compressional and shear waves; and transmitted, converted, and reflected waves and multiples. Full‐waveform inversion, therefore, has the potential to provide a more accurate result compared with conventional processing methods. We present a feasibility study with results of the application of high‐frequency (up to 60 Hz) anisotropic elastic full‐waveform inversion to a walkaway vertical seismic profile data from the Arabian Gulf. Full‐waveform inversion has reproduced the majority of the wave events and recovered a geologically plausible layered model with physically meaningful values of the medium. 相似文献
10.
We measured in the laboratory ultrasonic compressional and shear‐wave velocity and attenuation (0.7–1.0 MHz) and low‐frequency (2 Hz) electrical resistivity on 63 sandstone samples with a wide range of petrophysical properties to study the influence of reservoir porosity, permeability and clay content on the joint elastic‐electrical properties of reservoir sandstones. P‐ and S‐wave velocities were found to be linearly correlated with apparent electrical formation factor on a semi‐logarithmic scale for both clean and clay‐rich sandstones; P‐ and S‐wave attenuations showed a bell‐shaped correlation (partial for S‐waves) with apparent electrical formation factor. The joint elastic‐electrical properties provide a way to discriminate between sandstones with similar porosities but with different clay contents. The laboratory results can be used to estimate sandstone reservoir permeability from seismic velocity and apparent formation factor obtained from co‐located seismic and controlled source electromagnetic surveys. 相似文献
11.
In order to account for the effects of elastic wave propagation in marine seismic data, we develop a waveform inversion algorithm for acoustic‐elastic media based on a frequency‐domain finite‐element modelling technique. In our algorithm we minimize residuals using the conjugate gradient method, which back‐propagates the errors using reverse time migration without directly computing the partial derivative wavefields. Unlike a purely acoustic or purely elastic inversion algorithm, the Green's function matrix for our acoustic‐elastic algorithm is asymmetric. We are nonetheless able to achieve computational efficiency using modern numerical methods. Numerical examples show that our coupled inversion algorithm produces better velocity models than a purely acoustic inversion algorithm in a wide variety of cases, including both single‐ and multi‐component data and low‐cut filtered data. We also show that our algorithm performs at least equally well on real field data gathered in the Korean continental shelf. 相似文献
12.
Leonardo Azevedo Rúben Nunes Amílcar Soares Guenther Schwedersky Neto Teresa S. Martins 《Geophysical Prospecting》2018,66(Z1):116-131
Seismic inversion plays an important role in reservoir modelling and characterisation due to its potential for assessing the spatial distribution of the sub‐surface petro‐elastic properties. Seismic amplitude‐versus‐angle inversion methodologies allow to retrieve P‐wave and S‐wave velocities and density individually allowing a better characterisation of existing litho‐fluid facies. We present an iterative geostatistical seismic amplitude‐versus‐angle inversion algorithm that inverts pre‐stack seismic data, sorted by angle gather, directly for: density; P‐wave; and S‐wave velocity models. The proposed iterative geostatistical inverse procedure is based on the use of stochastic sequential simulation and co‐simulation algorithms as the perturbation technique of the model parametre space; and the use of a genetic algorithm as a global optimiser to make the simulated elastic models converge from iteration to iteration. All the elastic models simulated during the iterative procedure honour the marginal prior distributions of P‐wave velocity, S‐wave velocity and density estimated from the available well‐log data, and the corresponding joint distributions between density versus P‐wave velocity and P‐wave versus S‐wave velocity. We successfully tested and implemented the proposed inversion procedure on a pre‐stack synthetic dataset, built from a real reservoir, and on a real pre‐stack seismic dataset acquired over a deep‐water gas reservoir. In both cases the results show a good convergence between real and synthetic seismic and reliable high‐resolution elastic sub‐surface Earth models. 相似文献
13.
Finite‐difference P‐SV simulations of seismic scattering characteristics of faulted coal‐seam models have been undertaken for near‐surface P‐ and S‐wave sources in an attempt to understand the efficiency of body‐wave to channel‐wave mode conversion and how it depends on the elastic parameters of the structure. The synthetic seismograms clearly show the groups of channel waves generated at the fault: one by the downgoing P‐wave and the other by the downgoing S‐wave. These modes travel horizontally in the seam at velocities less than the S‐wavespeed of the rock. A strong Airy phase is generated for the fundamental mode. The velocity contrast between the coal and the host rock is a more important parameter than the density contrast in controlling the amplitude of the channel waves. The optimal coupling from body‐wave energy to channel‐wave energy occurs at a velocity contrast of 1.5. Strong guided waves are produced by the incident S‐sources for source angles of 75° to 90° (close to the near‐side face of the fault). As the fault throw increases, the amplitude of the channel wave also increases. The presence of a lower‐velocity clay layer within the coal‐seam sequence affects the waveguiding characteristics. The displacement amplitude distribution is shifted more towards the lower‐wavespeed layer. The presence of a ‘washout’ zone or a brecciated zone surrounding the fault also results in greater forward scattering and channel‐wave capture by the coal seam. 相似文献
14.
Inversion of multicomponent seismic time shifts for reservoir pressure and length: a feasibility study 下载免费PDF全文
下载免费PDF全文 Pressure drops associated with reservoir production generate excess stress and strain that cause travel‐time shifts of reflected waves. Here, we invert time shifts of P‐, S‐, and PS‐waves measured between baseline and monitor surveys for pressure reduction and reservoir length. The inversion results can be used to estimate compaction‐induced stress and strain changes around the reservoir. We implement a hybrid inversion algorithm that incorporates elements of gradient, global/genetic, and nearest neighbour methods and permits exploration of the parameter space while simultaneously following local misfit gradients. Our synthetic examples indicate that optimal estimates of reservoir pressure from P‐wave data can be obtained using the reflections from the reservoir top. For S‐waves, time shifts from the top of the reservoir can be accurately inverted for pressure if the noise level is low. However, if noise contamination is significant, it is preferable to use S‐wave data (or combined shifts of all three modes) from reflectors beneath the reservoir. Joint wave type inversions demonstrate improvements over any single pure mode. Reservoir length can be estimated using the time shifts of any mode from the reservoir top or deeper reflectors. We also evaluate the differences between the actual strain field and those corresponding to the best‐case inversion results obtained using P‐ and S‐wave data. Another series of tests addresses the inversion of the time shifts for the pressure drops in two‐compartment reservoirs, as well as for the associated strain field. Numerical testing shows that a potentially serious source of error in the inversion is a distortion in the strain‐sensitivity coefficients, which govern the magnitude of stiffness changes. This feasibility study suggests which wave types and reflector locations may provide the most accurate estimates of reservoir parameters from compaction‐induced time shifts. 相似文献
15.
利用天然地震记录直达P波后续20s内的波形信息,研究了地壳、土地慢速度结构反演的方法、唯一性、精度及应注意的问题,并将求解有条件极值的惩罚数法和求解无条件极值的单纯形法引入到波形反演.数值计算证明,波形反演方法对记录的误差有压制作用,在深源、远震、各向同性水平层状近似合理的前提下,可用一个地震事件在射线平面内的二分量记录反演台站区地壳、上地幔P、S波速度结构.在新疆天山地区转换波流动台站观测中,选择出较好的记录进行波形反演,给出了该转换被测深剖面上6个台站下方的地壳、上地慢P、S波速度结构.通过塔里木盆地内一测点上波形反演和地震勘探结果的对比证明,波形反演方法具有满意的精度. 相似文献
16.
Analysis of different parameterisations of waveform inversion of compressional body waves in an elastic transverse isotropic Earth with a vertical axis of symmetry 下载免费PDF全文
下载免费PDF全文 In a multi‐parameter waveform inversion, the choice of the parameterisation influences the results and their interpretations because leakages and the tradeoff between parameters can cause artefacts. We review the parameterisation selection when the inversion focuses on the recovery of the intermediate‐to‐long wavenumbers of the compressional velocities from the compressional body (P) waves. Assuming a transverse isotropic medium with a vertical axis of symmetry and weak anisotropy, analytical formulas for the radiation patterns are developed to quantify the tradeoff between the shear velocity and the anisotropic parameters and the effects of setting to zero the shear velocity in the acoustic approach. Because, in an anisotropic medium, the radiation patterns depend on the angle of the incident wave with respect to the vertical axis, two particular patterns are discussed: a transmission pattern when the ingoing and outgoing slowness vectors are parallel and a reflection pattern when the ingoing and outgoing slowness vectors satisfy Snell's law. When the inversion aims at recovering the long‐to‐intermediate wavenumbers of the compressional velocities from the P‐waves, we propose to base the parameterisation choice on the transmission patterns. Since the P‐wave events in surface seismic data do not constrain the background (smooth) vertical velocity due to the velocity/depth ambiguity, the preferred parameterisation contains a parameter that has a transmission pattern concentrated along the vertical axis. This parameter can be fixed during the inversion which reduces the size of the model space. The review of several parameterisations shows that the vertical velocity, the Thomsen parameter δ, or the Thomsen parameter ε have a transmission pattern along the vertical axis depending on the parameterisation choice. The review of the reflection patterns of those selected parameterisations should be done in the elastic context. Indeed, when reflection data are also inverted, there are potential leakages of the shear parameter at intermediate angles when we carry out acoustic inversion. 相似文献
17.
We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls. 相似文献
18.
Philip Tillotson Angus Ian Best Jeremy Sothcott Clive McCann Wang Shangxu Xiang‐Yang Li 《Geophysical Prospecting》2011,59(1):111-119
P‐ and S‐wave velocity and attenuation coefficients (accurate to ±0.3% and ±0.2 dB/cm, respectively) were measured in synthetic porous rocks with aligned, penny‐shaped fractures using the laboratory ultrasonic pulse‐echo method. Shear‐wave splitting was observed by rotating the S‐wave transducer and noting the maximum and minimum velocities relative to the fracture direction. A block of synthetic porous rock of fracture density 0.0201 ± 0.0068 and fracture size 3.6 ± 0.38 mm (measured from image analysis of X‐ray CT scans) was sub‐sampled into three 20–30 mm long, 50 mm diameter core plugs oriented at 0°, 45° and 90° to the fracture normal (transversely isotropic symmetry axis). Full waveform data were collected over the frequency range 500–1000 kHz for both water and glycerin saturated cores to observe the effect of pore fluid viscosity at 1 cP and 100 cP, respectively. The shear‐wave splitting observed in the 90° core was 2.15 ± 0.02% for water saturated and 2.39 ± 0.02% for glycerin saturated, in agreement with the theory that suggests that the percentage splitting should be 100 times the fracture density and independent of the saturating fluid. In the 45° core, by contrast, splitting was 0.00 ± 0.02% for water saturation and ?0.77 ± 0.02% for glycerin saturation. This dependence on fracture orientation and pore fluid viscosity is consistent with the poro‐visco‐elastic theory for aligned, meso‐scale fractures in porous rocks. The results suggest the possible use of shear‐ or converted‐wave data to discriminate between fluids on the basis of viscosity variations. 相似文献
19.
Alireza Malehmir Magnus Andersson Maxim Lebedev Milovan Urosevic Vassili Mikhaltsevitch 《Geophysical Prospecting》2013,61(1):153-167
To provide a guide for future deep (<1.5 km) seismic mineral exploration and to better understand the nature of reflections imaged by surface reflection seismic data in two mining camps and a carbonatite complex of Sweden, more than 50 rock and ore samples were collected and measured for their seismic velocities. The samples are geographically from the northern and central parts of Sweden, ranging from metallic ore deposits, meta‐volcanic and meta‐intrusive rocks to deformed and metamorphosed rocks. First, ultrasonic measurements of P‐ and S‐wave velocities at both atmospheric and elevated pressures, using 0.5 MHz P‐ and S‐wave transducers were conducted. The ultrasonic measurements suggest that most of the measured velocities show positive correlation with the density of the samples with an exception of a massive sulphide ore sample that shows significant low P‐ and S‐wave velocities. The low P‐ and S‐wave velocities are attributed to the mineral texture of the sample and partly lower pyrite content in comparison with a similar type sample obtained from Norway, which shows significantly higher P‐ and S‐wave velocities. Later, an iron ore sample from the central part of Sweden was measured using a low‐frequency (0.1–50 Hz) apparatus to provide comparison with the ultrasonic velocity measurements. The low‐frequency measurements indicate that the iron ore sample has minimal dispersion and attenuation. The iron ore sample shows the highest acoustic impedance among our samples suggesting that these deposits are favourable targets for seismic methods. This is further demonstrated by a real seismic section acquired over an iron ore mine in the central part of Sweden. Finally, a laser‐interferometer device was used to analyse elastic anisotropy of five rock samples taken from a major deformation zone in order to provide insights into the nature of reflections observed from the deformation zone. Up to 10% velocity‐anisotropy is estimated and demonstrated to be present for the samples taken from the deformation zone using the laser‐interferometery measurements. However, the origin of the reflections from the major deformation zone is attributed to a combination of anisotropy and amphibolite lenses within the deformation zone. 相似文献
20.
Full waveform inversion algorithms are widely used in the construction of subsurface velocity models. In the following study, we propose a Laplace–Fourier-domain waveform inversion algorithm that uses both Laplace-domain and Fourier-domain wavefields to achieve the reconstruction of subsurface velocity models. Although research on the Laplace–Fourier-domain waveform inversion has been published recently that study is limited to fluid media. Because the geophysical targets of marine seismic exploration are usually located within solid media, waveform inversion that is approximated to acoustic media is limited to the treatment of properly identified submarine geophysical features. In this study, we propose a full waveform inversion algorithm for isotropic fluid–solid media with irregular submarine topography comparable to a real marine environment. From the fluid–solid system, we obtained P and S wave velocity models from the pressure data alone. We also suggested strategies for choosing complex frequency bands constructed of frequencies and Laplace coefficients to improve the resolution of the restored velocity structures. For verification, we applied our Laplace–Fourier-domain waveform inversion for fluid–solid media to synthetic data that were reconstructed for fluid–solid media. Through this inversion test, we successfully restored reasonable velocity structures. Furthermore, we successfully extended our algorithm to a field data set. 相似文献