共查询到20条相似文献，搜索用时 15 毫秒
1.
Margarita Luneva 《Pure and Applied Geophysics》1996,148(1-2):113-136
Numerical examples of high-frequency synthetic seismograms of body waves in a 2-D layered medium with complex interfaces (faults, wedges, curvilinear, corrugated) are presented. The wave field modeling algorithm combines the possibilities of the ray method and the edge wave superposition method. This approach preserves all advantages of the ray method and eliminates restrictions related to diffraction by boundary edges and to caustic effects in singular regions. The method does not require two-point ray tracing (source-to-receiver), and the position of the source, as well as the type of source, and the position of receivers can be chosen arbitrarily. The memory and the time required for synthetic seismogram computation are similar to ray synthetic seismograms. The computation of the volume of the medium (the Fresnel volume or Fresnel zones), which gives the essential contribution to the wave field, is included in the modeling program package. In the case of complicated irregular interface (or a layered medium with a regular ray field at the last interface), the method displays a high accuracy of wave field computation. Otherwise, the method can be considered a modification of the ray method with regularization by the superposition of edge waves. 相似文献
2.
Michel Bouchon 《Pure and Applied Geophysics》1996,148(1-2):3-20
We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering. 相似文献
3.
地震诱导电磁现象是国内外地学领域十分关注的前沿问题,前人对地震波和电磁场耦合波场的认识主要是基于规则模型获得的.为研究含起伏地表和地下界面的地层中震电波场激发、传播特性,本文采用有限元软件COMSOL Multiphysics模拟点震源激发的电磁场.首先给出频率域二维SHTE模式震电耦合方程组,然后利用COMSOL软件建立计算模型,并求解出点力源激发震电波场的频率域响应,最后利用FFT变换得到地震波场和电磁场的时间域波形.模拟结果表明,震电波场中存在三种类型的电磁信号,第一种是震源直接激发的电磁波;第二种是地震波在分界面处激发的电磁波(包括自由表面、地下不同介质分界面);第三种是伴随地震波的同震信号,前两种电磁波比地震波更早到达远处观测台站,对地震预警有重要意义.此外,研究还发现:当地震波传播至地表并沿着地表传播时,在地表附近空气层中同样记录到了伴随地震波传播的电磁扰动信号,该信号与相同水平源距条件下、地下观测点接收到的电磁信号相同,这与前人的一些观测结果相符.本文研究结果为今后地震电磁信号的解释提供了理论证据. 相似文献
4.
近源时域电磁场具有信号强、探测深度大和精度高等优点,但传统瞬变电磁场理论中偶极子近似在近源区会引起较大误差,推导瞬变电磁场直接时域解析式是解决这一问题的关键.本文在点电荷微元假设下通过时域格林函数,采用分离变量等方法推导出了上半空间一次有源波动场和反射波的时域解析式和下半空间二次无源波动场的时域解析式,结合均匀半空间瞬变电磁场的边界条件给出了均匀半空间瞬变电磁场的直接时域解析式,进而利用第一型曲线积分,通过沿回线源叠加推导出圆回线源在瞬变电磁场中的直接时域解析式.然后在半空间表面上,与传统的电偶极源假设下的表达式作了比较.数值结果表明两者在远源区的计算结果相差甚微,而近源区则存在很大误差.本文利用真正点元（点电荷）严密推导给出的均匀半空间表面上瞬变电磁场的直接时域解析式适用于全场区探测,克服了偶极子假设下只适用远场区的不足,为瞬变电磁法的进一步发展和实际勘探提供了新的理论基础. 相似文献
5.
本文研究声电效应测井波场的有限差分模拟算法.忽略井外地层中诱导电磁场对孔隙弹性波的影响,将求解动电耦合波方程组的问题解耦,先计算孔隙弹性波,再计算其诱导电磁场.基于轴对称柱坐标系下的速度-应力交错网格,采用时域有限差分计算井孔流体声波和井外地层孔隙弹性波.将电磁场近似看作似稳场,基于轴对称柱坐标系下的5点式有限差分网格,求解不同时刻的电位Poisson方程,计算诱导电场.结果表明:本文算法可准确模拟频率6.0 kHz的声电效应测井全波;在声波测井频率范围内,电导率、动电耦合系数和动态渗透率的低频近似对伴随电磁场的计算影响不大;地层水平界面导致伴随反射斯通利波的电场和显著的界面电磁波,后者对于探测地层界面具有潜在的应用价值. 相似文献
6.
Using reciprocal theorems for dynamic and static boundary value problems, boundary integral equations are presented for wave propagation in elastic, isotropic media and compressible, inviscid fluids in the time domain as well as in the frequency domain. For the analysis of fluid–soil and fluid–structure systems, suitable coupling conditions are prescribed along the interfaces. The numerical treatment of the boundary integral equations consists of a point collocation and of a discretization of the boundary, in which constant and linear approximation functions are assumed. Step-by-step integration is applied to the time-dependent equations, where again the states are taken to be linear and constant over each time interval. These boundary element procedures are used to analyse the response of dams due to horizontal and vertical ground motions considering dam–water interaction and absorption of hydrodynamic pressure waves at the reservoir bottom or at the far end into the soil medium. Both the frequency response and the impulse generated transient response are investigated. 相似文献
7.
The surface motion during an earthquake is different from point to point depending on the propagation properties of the seismic waves. Rocking and torsion are thus present in the free field, in proportion to the spatial derivatives of the surface motion with respect to a given direction. These derivatives are inversely proportional to the apparent wave velocity in that direction, so the smaller the wave apparent velocity, the more important its contribution to the rotations. In this respect, a marked contribution to surface rotations from surface waves is expected. A mathematical model is presented, based on a detailed representation of soil impedance, an approximate identification of surface waves and a deconvolution of body waves in P and SV contributions. Through this model the surface motion obtained from the records of strong-motion accelerometers can be expressed as a superposition of plane waves of known wavelengths. Rocking response spectra are computed and results are compared with previously published spectra. A sensitivity analysis is performed on some parameters of the model. 相似文献
8.
Munirdin Tohti Yibo Wang Evert Slob Yikang Zheng Xu Chang Yi Yao 《Geophysical Prospecting》2020,68(6):1927-1943
Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ. 相似文献
9.
用波函数展开法求解界面圆孔的SH波散射问题 总被引：2，自引：1，他引：1
研究了ＳＨ波绕界面圆孔散射的波函数展开理论。由入射,反射和透射波组成的自由波场与孔的散射场叠加成总波场。将定义于两个半平面的散射场按一定方式延拓到全平面。通过Ｈａｎｋｅｌ－Ｆｏｕｒｉｅｒ型积分来满足界面与孔缘处的边界条件,得到了确定特定系数的封闭代数方程组。给出了具体算例,计算了孔缘处的动应力系数分布。 相似文献
10.
The seismic wave field is considerably influenced by local structures close to the source and to the receiver. This applies to sources and receivers situated close to localized inhomogeneities, to structural interfaces, to the earth's surface, etc. In this paper we concentrate our attention mainly to the ray-theoretical radiation patterns of point sources situated close to the structural interfaces and to the earth's surface. In numerical modeling of high-frequency seismic wave fields by the ray method, the interaction of the source with the earth's surface has not usually been taken into account.The proposed procedure of the computation of the radiation patterns of point sources situated directly on structural interfaces and on the earth's surface is based on the zero-order approximation of the ray method, assuming that the length of the ray between the source and the receiver is long. The derived equations are extended to point sources located close to structural interface, to the earth's surface and to thin transition layers using the hybrid ray-reflectivity method, seeervený (1989). The thin layer need not be homogeneous; it may include an arbitrary inner layering (transition layers, laminas, etc.) The only requirement is for the layer to be thin. Roughly speaking, we require its thickness to be less than one quarter of the prevailing wavelength. The hybrid ray-reflectivity method describes well even certain non-ray effects (tunneling.S
* waves, etc.). Explicit analytical expressions for radiation patterns for all above listed point sources are found. These expression have a local character and may be easily implemented into computer codes designed for the routine computation of ray amplitudes and synthetic ray seismograms in 2-D and 3-D, laterally varying isotropic layered and block structures by the ray method.Numerical examples of radiation patterns ofP andS waves of point sources situated close to the earth's surface and to a thin low-velocity surface layer are presented and discussed. The explosive point source (center of dilatation) and the vertical and horizontal single force point sources are considered. It has been ascertained that the radiation patterns of point sources depend drastically on the depth of the source below the surface even if the depths vary within one quarter of the prevailing wavelength. 相似文献
11.
RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA 
下载免费PDF全文
下载免费PDF全文 Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection. 相似文献
12.
The mechanical model for plane strain, time-harmonic seismic wave propagation problems in cracked, multi-layered geological regions with surface topography and non-parallel interfaces was described in the first part of this work. Here, this model is used to investigate the response of such a region to the presence of traveling elastic waves generated by a seismic source. The computational methodology that was developed in the first part is based on a combination of both the regular (displacement-based) and the hypersingular (traction-based) Boundary Integral Equation Method (BIEM). First, the accuracy and convergence characteristics of this hybrid BIEM are studied. Then, a series of problems involving four different configurations of a reference geological deposit with both interface and internal cracks are solved, for a loading that is due to a seismically-induced pressure wave propagating upwards from the underlying rigid half-plane. The purpose of the numerical study is to investigate the influence of various key parameters of the problem, such as frequency and incidence angle of the incoming wave, size of the surface relief, location and size of the buried cracks, interaction effects between cracks and finally the presence of layers, on both the scattered displacement field and the stress concentration field. 相似文献
13.
Expansion of a plane wave into Gaussian beams 总被引：1，自引：0，他引：1
Summary An integral expansion which expresses a plane monochromatic wave as a superposition of Gaussian beams is found. The expansion can be used to solve many wave propagation problems in complicated structures, including laterally inhomogeneous media with curved interfaces. 相似文献
14.
The formal solutions of displacement field to the problem of elastic wave scattering and diffraction due to an infinitely long rigid cylinder embedded in an infinite elastic medium by an impulsive point source have been obtained in the integral form. The integrals for the reflected and the diffracted waves both in the shadow zone and in the illuminated zone are evaluated asymptotically for the early time motion by the Reisdue-Cagniard method and the Saddle-point-Cagniard method.Numerical results of the diffractedP, S andPS waves at a fixed circum-distance from the surface of the rigid cylinder show noticeably that (1) the energy partition for the diffractedS wave is small in comparison with that for the diffractedP wave, (2) the wave form of the diffractedS wave is broader and more diffused than that of the diffractedP wave, (3) the direction of the radial motions of the diffractedP andS waves varies as a function of the observational point, and (4) the energy partition for the diffractedP wave is much smaller than that for the direct or the reflectedP waves.This paper has been presented at the 46th Annual International Meeting of Society of Exploration Geophysicists in Houston, Texas, Oct. 28, 1976. 相似文献
15.
COMPUTATION OF ZERO-OFFSET VERTICAL SEISMIC PROFILES INCLUDING GEOMETRICAL SPREADING AND ABSORPTION*
Synthetic vertical seismic profiles (VSP) provide a useful tool in the interpretation of VSP data, allowing the interpreter to analyze the propagation of seismic waves in the different layers. A zero-offset VSP modeling program can also be used as part of an inversion program for estimating the parameters in a layered model of the subsurface. Proposed methods for computing synthetic VSP are mostly based on plane waves in a horizontally layered elastic or anelastic medium. In order to compare these synthetic VSP with real data a common method is to scale the data with the spherical spreading factor of the primary reflections. This will in most cases lead to artificial enhancement of multiple reflections. We apply the ray series method to the equations of motion for a linear viscoelastic medium after having done a Fourier transformation with respect to the time variable. This results in a complex eikonal equation which, in general, appears to be difficult to solve. For vertically traveling waves in a horizontally layered viscoelastic medium the solution is easily found to be the integral along the ray of the inverse of the complex propagation velocity. The spherical spreading due to a point source is also complex, and it is equal to the integral along the ray of the complex propagation velocity. Synthetic data examples illustrate the differences between spherical, cylindrical, and plane waves in elastic and viscoelastic layered media. 相似文献
16.
Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green’s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method. 相似文献
17.
岩石和岩体是具有复杂细微观结构的非均匀介质.弹性波在岩体中传播时,与岩体细微观缺陷相互作用表现出弹性波的频散效应.为研究岩体内部细观结构对弹性波频散效应的作用,本文采用双裂纹模型:在模型内部,考虑裂纹间的相互作用对弹性波的影响,以分析弹性波在双裂纹体系间的多次散射作用;在双裂纹体系间,采用线性叠加分析法,以考虑岩体缺陷影响的局部化.对波动方程应用Green函数基本解,利用边界积分方法,将双裂纹体系作为内边界处理,得到相应的频散方程,由此对比分析了双裂纹体系在上述两种分析方法下的区别,进一步探讨了双裂纹体系参数、孔隙流体压力和卸荷对岩体频散特性的影响. 相似文献
18.
19.