首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 57 毫秒
1.
来自海底高速层径向波的理论地震图研究   总被引:1,自引:0,他引:1       下载免费PDF全文
本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波(简称径向波),这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.  相似文献   

2.
It is well established that the upper—cooler—part of the crust is brittle, while deeper zones present ductile behaviour. In some cases, this brittle–ductile transition is a single seismic reflector with an associated reflection coefficient. We first develop a stress–strain relation including the effects of crust anisotropy, seismic attenuation and ductility in which deformation takes place by shear plastic flow. Viscoelastic anisotropy is based on the eigenstrain model and the Zener and Burgers mechanical models are used to model the effects of seismic attenuation, velocity dispersion, and steady-state creep flow, respectively. The stiffness components of the brittle and ductile media depend on stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the reflection and transmission coefficients of a single brittle–ductile interface and of a ductile thin layer. The PP scattering coefficient has a Brewster angle (a sign change) in both cases, and there is substantial PS conversion at intermediate angles. The PP coefficient is sensitive to the layer thickness, unlike the SS coefficient. Thick layers have a well-defined Brewster angle and show higher reflection amplitudes. Finally, we compute synthetic seismograms in a homogeneous medium as a function of temperature.  相似文献   

3.
The construction of S-wave velocity models of marine sediments down to hundreds of meters below the seafloor is important in a number of disciplines. One of the most significant trends in marine geophysics is to use interface waves to estimate shallow shear velocities which play an important role in determining the shallow crustal structure. In marine settings, the waves trapped near the fluid–solid interface are called Scholte waves, and this is the subject of the study. In 1998, there were experiments on the Ninetyeast Ridge (Central Indian Ocean) to study the shallow seismic structure at the drilled site. The data were acquired by both ocean bottom seismometer and ocean bottom hydrophone. A new type of seafloor implosion sources has been used in this experiment, which successfully excited fast and high frequency (>500 Hz) body waves and slow, intermediate frequency (<20 Hz) Scholte waves. The fundamental and first higher mode Scholte waves have both been excited by the implosion source. Here, the Scholte waves are investigated with a full waveform modeling and a group velocity inversion approach. Shear wave velocities for the uppermost layers of the region are inferred and results from the different methods are compared. We find that the full waveform modeling is important to understand the intrinsic attenuation of the Scholte waves between 1 and 20 Hz. The modeling shows that the S-wave velocity varies from 195 to 350 m/s in the first 16 m of the uppermost layer. Depths levels of high S-wave impedance contrasts compare well to the layer depth derived from a P-wave analysis as well as from drilling data. As expected, the P- to S-wave velocity ratio is very high in the uppermost 16 m of the seafloor and the Poisson ratio is nearly 0.5. Depth levels of high S-wave impedance contrasts are comparable to the layer depth derived from drilling data.  相似文献   

4.
Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients.  相似文献   

5.
利用非均匀P-偏振电磁波(电场平行于入射面)在导电界面反射系数所产生的附加相角导出了非均匀电磁波在导电介质界面的横向偏移,并给出了一个低频电磁波算例,绘制了横向偏移随入射角的变化曲线. 计算结果表明,入射角在相移常数临界角、衰减常数临界角和90°处,横向偏移曲线存在三个间断点;当入射角等于这三个角时,电磁波将沿界面传播;当入射角在这三个角附近时,横向偏移效应影响较大;在其他角度入射时横向偏移与波长在一个数量级. 论文还讨论了横向偏移效应对电磁波测井测量结果所产生的误差,并给出了误差计算方法和算例.  相似文献   

6.
Synthetic seismograms are usually computed for reflections from vertical incidence of P waves for a horizontally layered medium. In actual practice the angle of incidence departs from the vertical, as receivers are usually located at some distance from the source. At angles other than the vertical, the conversion of P- to S-wave energy and changes in the reflection coefficient affect the shape of the synthetic seismograms. The effect of non-vertical incidence on synthetic seismograms is examined in this paper. Seismograms at non-vertical incidence have been computed using the plane-wave approach of Haskell (1953) for a layered medium. The use of plane waves is an approximation to the actual case of spherical wavefronts from a surface source. Using plane-wave theory, the expected wave forms as a function of angle of incidence were computed numerically for several simple models. The results indicate that the synthetic seismograms do not change significantly for angles of incidence between o and 25 degrees. For larger angles the changes in the wave forms may be severe. The effect is more pronounced for high-velocity layers than for low-velocity layers.  相似文献   

7.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大.  相似文献   

8.
Using Biot’s poroelasticity theory, we derive expressions for the reflection and transmission coefficients for a plane shear wave incident on an interface separating two different poroelastic solids. The coefficients are formulated as a function of the wave incidence angle, frequency and rock properties. Specific cases calculated include the boundary between water-saturated sand and water-saturated sandstone and the gas–water interface in sand. The results show a very different interface response to that of an incident P wave. Plane SV wave incidence does not significantly excite the Biot slow P wave if the frequency of the wave is below the transition frequency. Above this frequency, an incident plane SV wave can generate a mode-converted slow Biot P wave which is actually a normal propagating wave and not highly attenuating as in the usual (diffusive) case. For an incident SV wave onto a gas–water interface, even at very high frequency, there is no significant Biot second P wave produced. For small incident angles, the gas–water interface is essentially transparent. With increasing angles, there can arise an unusual "definitive angle" in the reflection/transmission coefficient curves which is related to the change of fluid viscosity on both sides of the interface and provides a possible new means for underground fluid assessment.  相似文献   

9.
利用电磁波在导电媒质界面反射系数的附加相角,导出了电场垂直于入射面的非均匀电磁波在导电介质界面的类全反射横向偏移,并对横向偏移进行了相关计算,绘出了横向偏移随入射角的变化曲线. 结果表明:入射角在相移常数临界角、衰减常数临界角和90°处,横向偏移曲线存在三个间断点. 当入射角等于这三个角时,横向偏移为无穷大,即电磁波将沿界面传播;当入射角在这三个角附近时,横向偏移变得非常大.  相似文献   

10.
In this paper, we derive an approximation of the SS-wave reflection coefficient and the expression of S-wave ray elastic impedance (SREI) in terms of the ray parameter. The SREI can be expressed by the S-wave incidence angle or P-wave reflection angle, referred to as SREIS and SREIP, respectively. Our study using elastic models derived from real log measurements shows that SREIP has better capability for lithology and fluid discrimination than SREIS and conventional S-wave elastic impedance (SEI). We evaluate the SREIP feasibility using 25 groups of samples from Castagna and Smith (1994). Each sample group is constructed by using shale, brine-sand, and gas-sand. Theoretical evaluation also indicates that SREIP at large incident angles is more sensitive to fluid than conventional fluid indicators. Real seismic data application also shows that SREIP at large angles calculated using P-wave and S-wave impedance can efficiently characterize tight gas-sand.  相似文献   

11.
线性波浪加载下海底斜坡失稳机制的数值分析   总被引:2,自引:0,他引:2       下载免费PDF全文
基于大型有限元软件ABAQUS中的荷载模块,添加一阶波浪力载荷模式,并结合强度折减技术,实现波浪力作用下海底斜坡稳定性与失稳机制的弹塑性有限元数值分析。引入典型算例,利用先前提出的波浪荷载下海底斜坡稳定性的极限分析上限方法开展数值解的对比验证;在此基础上,通过深入地变动参数比较分析,探讨不同波长、波高和水深等波浪参数对计算结果的影响以及波浪力影响下海底斜坡潜在滑动面的变化规律,获得波浪荷载下海底斜坡失稳滑动机制的初步认识。  相似文献   

12.
An exact analysis for two-dimensional dynamic interaction of monochromatic progressive plane compressional and shear seismic waves with a permeable circular tunnel lining of circumferentially varying wall thickness buried in a boundless porous elastic fluid-saturated formation is presented. The novel features of Biot dynamic theory of poroelasticity in conjunction with the translational addition theorems for cylindrical wave functions, along with the appropriate wave field expansions and the pertinent boundary conditions are employed to develop a closed-form solution in form of infinite series. The analytical results are illustrated with numerical examples in which an air-filled and water-saturated permeable tunnel lining of variable wall thickness, embedded within water-saturated surrounding formations of distinct frame properties (soft, stiff, and stiff viscoelastic soil), is insonified by fast compressional or shear waves at selected angles of incidence. The effects of liner eccentricity, interface permeability, formation material type, incident wave frequency, and angle of incidence on the hoop stress amplitude are evaluated and discussed for representative values of the parameters characterizing the system. Limiting cases are considered and good agreements with the solutions available in the literature are obtained.  相似文献   

13.
Imaging the PP- and PS-wave for the elastic vector wave reverse-time migration requires separating the P- and S-waves during the wave field extrapolation. The amplitude and phase of the P- and S-waves are distorted when divergence and curl operators are used to separate the P- and S-waves. We present a P- and S-wave amplitude-preserving separation algorithm for the elastic wavefield extrapolation. First, we add the P-wave pressure and P-wave vibration velocity equation to the conventional elastic wave equation to decompose the P- and S-wave vectors. Then, we synthesize the scalar P- and S-wave from the vector Pand S-wave to obtain the scalar P- and S-wave. The amplitude-preserved separated P- and S-waves are imaged based on the vector wave reverse-time migration (RTM). This method ensures that the amplitude and phase of the separated P- and S-wave remain unchanged compared with the divergence and curl operators. In addition, after decomposition, the P-wave pressure and vibration velocity can be used to suppress the interlayer reflection noise and to correct the S-wave polarity. This improves the image quality of P- and S-wave in multicomponent seismic data and the true-amplitude elastic reverse time migration used in prestack inversion.  相似文献   

14.
Minimization of seismic residuals does not guarantee uniqueness of the model, and this implies ambiguities in the inversion. Amplitude vs. offset (AVO) inversion does not lead to a unique solution of single elastic interface parameters unless converted and S-wave or critical angle reflections are available. Given the ambiguity of AVO inversion, this paper discusses the interaction between AVO and velocity estimation. The number of independent parameters necessary to describe an isolated reflection with AVO behaviour and residual velocity error is determined. Statistical analysis allows the establishment of an approximate equivalence of the effects of AVO and slight velocity variations; this equivalence cannot be solved without geological a priori information (kinematic equivalence). The data are then decomposed into compound events (i.e. sequences of N interfaces that follow each other at a fixed time lag). The decomposition is obtained by extrapolating the results of the analysis from narrowband to wideband data. Compound events decomposition demonstrates that AVO inversion is ambiguous, not only in the physical parameter space (P- and S-wave velocities, and density) but also kinematically. As an example of compound event decomposition, a medium is derived. This medium is geologically implausible but is kinematically equivalent.  相似文献   

15.
光滑处理使得单界面成为非均匀薄层,界面反射转变为层反射.为了探讨光滑处理的影响,以平面波作为入射波场,首先利用过渡层反射系数推导了反射信号的理论公式,然后就非均匀薄层下反射系数的计算问题,给出了具体的实现算法,并通过与经典Epstein过渡层反射系数解析结果的对比说明了算法的精度.最后在单界面及其被光滑后界面的对比分析中,得出了几点重要结论:随着光滑次数的增加,反射信号的到时基本保持不变,而反射信号的主频与能量呈减少趋势,其中信号能量在低光滑次数的衰减速率明显大于高光滑次数.  相似文献   

16.
It has been known since the beginning of reflection seismics that several disturbing events seen in seismic records are caused by waves with S-wave velocities instead of P-wave velocity. When using dynamite and recording with vertical geophones these events are primarily caused by converted waves. On the basis of known P- and S-wave velocities in a certain area a theoretical seismogram is calculated, displaying traveltime as well as energy relation for different wave configurations. By comparison with seismograms recorded in the same area it can be shown that converted wave events can be clearly recognized. These events can be described theoretically. Thus, either more effective computer programs can be applied to eliminate these disturbing events, or these events can be evaluated to get additional information about specific strata.  相似文献   

17.
套管井声波全波列中地层波的幅度特性   总被引:1,自引:0,他引:1       下载免费PDF全文
套管井声全渡中的地层波可以检测固井第二界面的水泥胶结质量状况.地层波不仅受地层衰减、套管波干扰等因素影响,更主要受第二界面两侧的声阻抗差值的影响.为此,基于固体介质的声势函数和界面的边界条件,推导了固-固界面纵横波的反射系数和透射系数,以及第一临界角入射情况下反射系数和透射系数的近似公式,并且分析了透射系数与声阻抗之间的关系.数值计算结果发现,地层波幅度随第二界面两侧声阻抗差增大而减小,且这种关系是非线性的.基于临界入射的透射系数,对经过衰减校正和干扰校正的地层波进行声阻抗校正,消除地层声阻抗对地层波幅度影响.使用地层波幅度评价第二界面的水泥胶结状况,并在实际资料处理中见到了好的效果.  相似文献   

18.
在地震勘探中,P波和S波入射到一固体液体分界面处时,在该分界处的水平检波器和垂直检波器将接受到相位的响应垂直的响应。垂直检波器的响应与在自由界面处的垂直检波器 响应大概相同,而水平检波器对P波的响应相对于对S波的 来说要强。  相似文献   

19.
Cracks play a very important role in many geotechnical issues and in a number of processes in the Earth’s crust. Elastic waves can be used as a remote sensing tool for determining crack density. The effect of varying crack density in crystalline rock on the P- and S-wave velocity and dynamic elastic properties under confining pressure has been quantified. The evolution of P- and S-wave velocity were monitored as a suite of dry Westerly granite samples were taken to 60, 70, 80 and 90 % of the unconfined uniaxial strength of the sample. The damaged samples were then subjected to hydrostatic confining pressure from 2 MPa to 200 MPa to quantify the effect of varying crack density on the P- and S-wave velocity and elastic properties under confining pressure. The opening and propagation of microcracks predominantly parallel to the loading direction during uniaxial loading caused a 0.5 and 6.3 % decrease in the P- and S-wave velocity, respectively. During hydrostatic loading, microcracks are closed at 130 MPa confining pressure. At lower pressures the amount of crack damage in the samples has a small but measureable effect. We observed a systematic 6 and 4 % reduction in P- and S-wave velocity, respectively, due to an increase in the fracture density at 2 MPa confining pressure. The overall reduction in the P- and S-wave velocity decreased to 2 and 1 %, respectively, at 50 MPa. The elastic wave velocities of samples that have a greater amount of microcrack damage are more sensitive to pressure. Effective medium modelling was used to invert elastic wave velocities and infer crack density evolution. Comparing the crack density results with experimental data on Westerly granite samples shows that the effective medium modelling used gave interpretable and reasonable results. Changes in crack density can be interpreted as closure or opening of cracks and crack growth.  相似文献   

20.
Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy.  相似文献   

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

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