首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
  相似文献   

2.
Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.  相似文献   

3.
We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.  相似文献   

4.
Our objective was to discover the effect of variations in fluid properties and fracture geometry on the velocity of seismic wave propagation in fluid‐saturated media with parallel planar fractures. We used numerical models calculated by analytical solutions to examine the behaviour of P‐wave phase velocity dispersion in the normal direction to layering, in non‐porous and porous media with planar fractures. We also examined the anisotropy of low frequency phase and group velocities of fast and slow P‐waves and angular‐dependent reflection coefficients in media with planar fractures, under conditions of saturation by fluids with varying bulk moduli, densities, and fracture apertures. We defined several parametre , , and characterising dispersion, characterising anisotropy, characterising the difference between fast and slow modes, and R0 and characterising reflection. Our results show that the behaviour of dispersion shows wider stopbands in the case of gas saturation. Concavity indicator of dispersion for gas saturation was greater than that for liquid saturation and is usually greater than one. Anisotropy is more sensitive to bulk modulus contrast than to density contrast between the solid and the fluid, and is more sensitive to density contrast than to bulk modulus contrast. The case of gas saturation usually had a greater negative R0 and a greater value of compared with those of brine and heavy and light oil saturations. Our results are helpful in distinguishing fluid types saturating geophysical fractures and estimating the aperture and spacing of planar fractures. In seismic exploration, bulk modulus and fluid density can provide useful information in distinguishing among brine, oil, and gas; fracture geometry is important to estimate the permeability of reservoirs.  相似文献   

5.
Under certain circumstances, seismic propagation within porous media may be associated to the conversion of mechanical energy to electromagnetic energy, which is known as a seismo‐electromagnetic phenomenon. The propagation of fast compressional P‐waves is more specifically associated to the manifestations of a seismoelectric field linked to the fluid flows within the pores. The analysis of seismoelectric phenomena, which requires the combination of the theory of electrokinetics and Biot's theory of poroelasticity, provides us with transfer function that links the coseismic seismoelectric field E to the seismic acceleration . To measure the transfer function, we have developed an experimental setup enabling seismoelectric laboratory observation in unconsolidated quartz sand within the kilohertz range. The investigation focused on the impact of fluid conductivity and water saturation over the coseismic seismoelectric field. During the experiment, special attention was given to the accuracy of electric field measurements. We concluded that, to obtain a reliable estimate of the electric field amplitude, the dipole from which the potential differences are measured should be of much smaller length than the wavelength of the propagating seismic field. Time‐lapse monitoring of the seismic velocities and seismoelectric transfer functions were performed during imbibition and drainage experiments. In all cases, the quantitative analysis of the seismoelectric transfer function was in good agreement with theoretical predictions. While investigating saturation variations from full to residual water saturation, we showed that the ratio undergoes a switch in polarity at a particular saturation , which also implies a sign change of the filtration, traducing a reversal of the relative fluid displacement with respect to the frame. This sign change at critical saturation stresses a particular behaviour of the poroelastic medium: the dropping of the coseismic electric field to zero traduces the absence of relative pore/fluid displacements representative of a Biot dynamically compatible medium. We concluded from our experimental study in loose sand that the measurements of the coseismic seismoelectric coupling may provide information on fluid distribution within the pores and that the reversal of the seismoelectric field may be used as an indicator of the dynamically compatible state of the medium.  相似文献   

6.
We investigate fracture‐induced attenuation anisotropy in a cluster of events from a microseismic dataset acquired during hydraulic fracture stimulation. The dataset contains 888 events of magnitude ?3.0 to 0.0. We use a log‐spectral‐amplitude‐ratio method to estimate change in over a half‐hour time period where fluid is being injected and an increase in fracturing from S‐wave splitting analysis has been previously inferred. A Pearson's correlation analysis is used to assess whether or not changes in attenuation with time are statistically significant. P‐waves show no systematic change in during this time. In contrast, S‐waves polarised perpendicular to the fractures show a clear and statistically significant increase with time, whereas S‐waves polarised parallel to the fractures show a weak negative trend. We also compare between the two S‐waves, finding an increase in with time. A poroelastic rock physics model of fracture‐induced attenuation anisotropy is used to interpret the results. This model suggests that the observed changes in t* are related to an increase in fracture density of up to 0.04. This is much higher than previous estimates of 0.025 ± 0.002 based on S‐wave velocity anisotropy, but there is considerably more scatter in the attenuation measurements. This could be due to the added sensitivity of attenuation measurement to non‐aligned fractures, fracture shape, and fluid properties. Nevertheless, this pilot study shows that attenuation measurements are sensitive to fracture properties such as fracture density and aspect ratio.  相似文献   

7.
We have derived a convergent scattering series solution for the frequency-domain wave equation in acoustic media with variable density and velocity. The convergent scattering series solution is based on the homotopy analysis of a vectorial integral equation of the Lippmann–Schwinger type. By using the Green's function and partial integration, we have derived the vectorial integral equation of the Lippmann–Schwinger type that involves the pressure gradient field as well as the pressure field from the wave equation. The vectorial Lippmann–Schwinger equation can in principle be solved via matrix inversion, but the computational cost of matrix inversion scales like N 3 , where N is the number of grid blocks. The computational cost can be significantly reduced if one solves the vectorial Lippmann–Schwinger equation iteratively. A simple iterative solution is the Born series, but it is only convergent when the scattering potential is sufficiently small. In this study, we have used the so-called homotopy analysis method to derive an iterative solution for the vectorial Lippmann–Schwinger equation which can be made convergent even in strongly scattering media. The computational cost of our convergent scattering series scales as N 2 . Our algorithm, which is based on the homotopy analysis method, involves a convergence control operator that we select using hierarchical matrices. We use a three-layer model and a resampled version of the SEG/EAGE salt model to show the performance of the developed convergent scattering series.  相似文献   

8.
Logistic regression, neural networks and support vector machines are tested for their effectiveness in isolating surface waves in seismic shot records. To distinguish surface waves from other arrivals, we train the algorithms on three distinguishing features of surface-wave dispersion curves in the domain: spectrum coherency of the trace's magnitude spectrum, local dip and the frequency range for a fixed wavenumber k in the spectrum. Numerical tests on synthetic data show that the kernel-based support vector machines algorithm gives the highest accuracy in predicting the surface-wave window in the domain compared to neural networks and logistic regression. This window is also used to automatically pick the fundamental dispersion curve. The other two methods correctly pick the low-frequency part of the dispersion curve but fail at higher frequencies where there is interference with higher-order modes.  相似文献   

9.
Linear prediction filters are an effective tool for reducing random noise from seismic records. Unfortunately, the ability of prediction filters to enhance seismic records deteriorates when the data are contaminated by erratic noise. Erratic noise in this article designates non‐Gaussian noise that consists of large isolated events with known or unknown distribution. We propose a robust f‐x projection filtering scheme for simultaneous erratic noise and Gaussian random noise attenuation. Instead of adopting the ?2‐norm, as commonly used in the conventional design of f‐x filters, we utilize the hybrid ‐norm to penalize the energy of the additive noise. The estimation of the prediction error filter and the additive noise sequence are performed in an alternating fashion. First, the additive noise sequence is fixed, and the prediction error filter is estimated via the least‐squares solution of a system of linear equations. Then, the prediction error filter is fixed, and the additive noise sequence is estimated through a cost function containing a hybrid ‐norm that prevents erratic noise to influence the final solution. In other words, we proposed and designed a robust M‐estimate of a special autoregressive moving‐average model in the f‐x domain. Synthetic and field data examples are used to evaluate the performance of the proposed algorithm.  相似文献   

10.
While seismic reflection amplitudes are generally determined by real acoustical impedance contrasts, there has been recent interest in reflections due to contrasts in seismic‐Q. Herein we compare theoretical and modelled seismic reflection amplitudes for two different cases of material contrasts. In case A, we examine reflections from material interfaces that have a large contrast in real‐valued impedance () with virtually no contrast in seismic‐Q. In case B, we examine reflections from material interfaces that have virtually no contrast in but that have very large seismic‐Q contrasts. The complex‐valued reflection coefficient formula predicts non‐zero seismic reflection amplitudes for both cases. We choose physical materials that typify the physics of both case A and case B. Physical modelling experiments show significantly large reflections for both cases – with the reflections in the two cases being phase shifted with respect to each other, as predicted theoretically. While these modelling experiments show the existence of reflections that are predicted by theory, there are still intriguing questions regarding the size of the Q‐contrast reflections, the existence of large Q‐contrast reflections in reservoir rocks and the possible application of Q‐reflection analysis to viscosity estimation in heavy oilfields.  相似文献   

11.
When modelling the propagation of 3D non‐axisymmetric elastic and viscoelastic waves in cylindrical coordinates using the finite‐difference time‐domain method, a mathematical singularity occurs due to the presence of terms in the elastic and viscoelastic wave equations. For many years, this issue has been impeding the accurate numerical solution near the axis. In this work, we propose a simple but effective method for the treatment of this numerical singularity problem. By rotating the Cartesian coordinate system around the z‐axis in cylindrical coordinates, the numerical singularity problems in both 2D and 3D cylindrical coordinates can be removed. This algorithm has three advantages over the conventional treatment techniques: (i) the excitation source can be directly loaded at , (ii) the central difference scheme with second‐order accuracy is maintained, and (iii) the stability condition at the axis is consistent with the finite‐difference time‐domain in Cartesian coordinates. This method is verified by several 3D numerical examples. Results show that the rotating the Cartesian coordinate method is accurate and stable at the singularity axis. The improved finite‐difference time‐domain algorithm is also applied to sonic logging simulations in non‐axisymmetric formations and sources.  相似文献   

12.
The reflectivity method plays an important role in seismic modelling. It has been used to model different types of waves propagating in elastic and anelastic media. The diffusive–viscous wave equation was proposed to investigate the relationship between frequency dependence of reflections and fluid saturation. It is also used to describe the attenuation property of seismic wave in a fluid‐saturated medium. The attenuation of diffusive–viscous wave is mainly characterised by the effective attenuation parameters in the equation. Thus, it is essential to obtain those parameters and further characterise the features of the diffusive–viscous wave. In this work, we use inversion method to obtain the effective attenuation parameters through quality factor to investigate the characteristics of diffusive–viscous wave by comparing with those of the viscoacoustic wave. Then, the reflection/transmission coefficients in a dip plane‐layered medium are studied through coordinate transform and plane‐wave theory. Consequently, the reflectivity method is extended to compute seismograms of diffusive–viscous wave in a dip plane multi‐layered medium. Finally, we present two models to simulate the propagation of diffusive–viscous wave in a dip plane multi‐layered medium by comparing the results with those in a viscoacoustic medium. The numerical results demonstrate the validity of our extension of reflectivity method to the diffusive–viscous medium. The numerical examples in both time domain and time–frequency domain show that the reflections from a dip plane interface have significant phase shift and amplitude change compared with the results of horizontal plane interface due to the differences in reflection/transmission coefficients. Moreover, the modelling results show strong attenuation and phase shift in the diffusive–viscous wave compared to those of the viscoacoustic wave.  相似文献   

13.
We present a parsimonious wave‐equation travel‐time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave‐equation travel‐time inversion. Numerical results show that the parsimonious wave‐equation travel‐time tomogram has about the same accuracy as the tomogram computed by standard wave‐equation travel‐time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location.  相似文献   

14.
Ill-posedness is one of the most common and intractable issues that arise when solving geophysical inverse problems. Ill-posedness could be induced by various factors such as noise, band-limited intrinsic property of seismic data and inappropriate forward operators. Regularization has been proven to be an effective method widely accepted for mitigating the adverse effects of ill-posedness. Aiming to improve the stability and fidelity of the pre-stack seismic inversion process, we implement the inversion in a Bayesian framework, with a logarithmic absolute criterion taken as a likelihood function, and an -norm metric as a priori constraint. Here, we exploit the linear approximation as the forward operator, and optimize the regularized misfit function by the alternating direction method of multipliers. Applications of the method to synthetic and real data sets yielded improved inversion results in terms of accuracy and resolution, and demonstrated the robustness of the method to noise.  相似文献   

15.
A transmission + reflection wave-equation traveltime and waveform inversion method is presented that inverts the seismic data for the anisotropic parameters in a vertical transverse isotropic medium. The simultaneous inversion of anisotropic parameters and ε is initially performed using transmission wave-equation traveltime inversion method. Transmission wave-equation traveltime only provides the low-intermediate wavenumbers for the shallow part of the anisotropic model; in contrast, reflection wave-equation traveltime estimates the anisotropic parameters in the deeper section of the model. By incorporating a layer-stripping method with reflection wave-equation traveltime, the ambiguity between the background-velocity model and the depths of reflectors can be greatly mitigated. In the final step, multi-scale full-waveform inversion is performed to recover the high-wavenumber component of the model.  We use a synthetic model to illustrate the local minima problem of full-waveform inversion and how transmission and reflection wave-equation traveltime can mitigate this problem. We demonstrate the efficacy of our new method using field data from the Gulf of Mexico.  相似文献   

16.
Tensor algebra provides a robust framework for multi-dimensional seismic data processing. A low-rank tensor can represent a noise-free seismic data volume. Additive random noise will increase the rank of the tensor. Hence, tensor rank-reduction techniques can be used to filter random noise. Our filtering method adopts the Candecomp/Parafac decomposition to approximates a N-dimensional seismic data volume via the superposition of rank-one tensors. Similar to the singular value decomposition for matrices, a low-rank Candecomp/Parafac decomposition can capture the signal and exclude random noise in situations where a low-rank tensor can represent the ideal noise-free seismic volume. The alternating least squares method is adopted to compute the Candecomp/Parafac decomposition with a provided target rank. This method involves solving a series of highly over-determined linear least-squares subproblems. To improve the efficiency of the alternating least squares algorithm, we uniformly randomly sample equations of the linear least-squares subproblems to reduce the size of the problem significantly. The computational overhead is further reduced by avoiding unfolding and folding large dense tensors. We investigate the applicability of the randomized Candecomp/Parafac decomposition for incoherent noise attenuation via experiments conducted on a synthetic dataset and field data seismic volumes. We also compare the proposed algorithm (randomized Candecomp/Parafac decomposition) against multi-dimensional singular spectrum analysis and classical prediction filtering. We conclude the proposed approach can achieve slightly better denoising performance in terms of signal-to-noise ratio enhancement than traditional methods, but with a less computational cost.  相似文献   

17.
The conventional velocity scan can be computationally expensive for large‐scale seismic data sets, particularly when the presence of anisotropy requires multiparameter scanning. We introduce a fast algorithm for 3D azimuthally anisotropic velocity scan by generalizing the previously proposed 2D butterfly algorithm for hyperbolic Radon transforms. To compute semblance in a two‐parameter residual moveout domain, the numerical complexity of our algorithm is roughly as opposed to of the straightforward velocity scan, with N being the representative of the number of points in a particular dimension of either data space or parameter space. Synthetic and field data examples demonstrate the superior efficiency of the proposed algorithm.  相似文献   

18.
The nonlinearity of the seismic amplitude‐variation‐with‐offset response is investigated with physical modelling data. Nonlinearity in amplitude‐variation‐with‐offset becomes important in the presence of large relative changes in acoustic and elastic medium properties. A procedure for pre‐processing physical modelling reflection data is enacted on the reflection from a water‐plexiglas boundary. The resulting picked and processed amplitudes are compared with the exact solutions of the plane‐wave Zoeppritz equations, as well as approximations that are first, second, and third order in , , and . In the low angle range of 0°–20°, the third‐order plane‐wave approximation is sufficient to capture the nonlinearity of the amplitude‐variation‐with‐offset response of a liquid‐solid boundary with , , and ρ contrasts of 1485–2745 m/s, 0–1380 m/s, and 1.00–1.19 gm/cc respectively, to an accuracy value of roughly 1%. This is in contrast to the linear Aki–Richards approximation, which is in error by as much as 25% in the same angle range. Even‐order nonlinear corrective terms are observed to be primarily involved in correcting the angle dependence of , whereas the odd‐order nonlinear terms are involved in determining the absolute amplitude‐variation‐with‐offset magnitudes.  相似文献   

19.
Surface waves are often used to estimate a near‐surface shear‐velocity profile. The inverse problem is solved for the locally one‐dimensional problem of a set of homogeneous horizontal elastic layers. The result is a set of shear velocities, one for each layer. To obtain a P‐wave velocity profile, the P‐guided waves should be included in the inversion scheme. As an alternative to a multi‐layered model, we consider a simple smooth acoustic constant‐density velocity model, which has a negative constant vertical depth gradient of the squared P‐wave slowness and is bounded by a free surface at the top and a homogeneous half‐space at the bottom. The exact solution involves Airy functions and provides an analytical expression for the dispersion equation. If the ratio is sufficiently small, the dispersion curves can be picked from the seismic data and inverted for the continuous P‐wave velocity profile. The potential advantages of our model are its low computational cost and the fact that the result can serve as a smooth starting model for full‐waveform inversion. For the latter, a smooth initial model is often preferred over a rough one. We test the inversion approach on synthetic elastic data computed for a single‐layer P‐wave model and on field data, both with a small ratio. We find that a single‐layer model can recover either the shallow or deeper part of the profile but not both, when compared with the result of a multi‐layer inversion that we use as a reference. An extension of our analytic model to two layers above a homogeneous half‐space, each with a constant vertical gradient of the squared P‐wave slowness and connected in a continuous manner, improves the fit of the picked dispersion curves. The resulting profile resembles a smooth approximation of the multi‐layered one but contains, of course, less detail. As it turns out, our method does not degrade as gracefully as, for instance, diving‐wave tomography, and we can only hope to fit a subset of the dispersion curves. Therefore, the applicability of the method is limited to cases where the ratio is small and the profile is sufficiently simple. A further extension of the two‐layer model to more layers, each with a constant depth gradient of the squared slowness, might improve the fit of the modal structure but at an increased cost.  相似文献   

20.
Pre‐stack seismic data are indicative of subsurface elastic properties within the amplitude versus offset characteristic and can be used to detect elastic rock property changes caused by injection. We perform time‐lapse pre‐stack 3‐D seismic data analysis for monitoring sequestration at Cranfield. The time‐lapse amplitude differences of Cranfield datasets are found entangled with time‐shifts. To disentangle these two characters, we apply a local‐correlation‐based warping method to register the time‐lapse pre‐stack datasets, which can effectively separate the time‐shift from the time‐lapse seismic amplitude difference without changing the original amplitudes. We demonstrate the effectiveness of our registration method by evaluating the inverted elastic properties. These inverted time‐lapse elastic properties can be reliably used for monitoring plumes.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号