共查询到20条相似文献，搜索用时 15 毫秒
1.
本文采用复变函数方法研究了稳态水平剪切波(SH波)在各向异性弹性半空间中任意形状孔洞上的散射及其对地面运动的影响.在变换空间中构造出自动满足各向异性半空间水平表面上应力为零的散射波函数,并根据孔洞表面应力为零的边界条件,采用最小二乘法求解散射波函数的系数.用介质的各向异性性质来模拟地质条件,给出了SH波作用下含圆孔、椭圆孔和方孔的各向异性半空间表面位移幅值的数值结果,并分析了介质特性、孔洞的形状、埋深、入射波波数及入射角度等因素对地面运动的影响规律.数值结果表明:介质的各向异性对含有孔洞的半空间表面的地表位移具有显著的影响;沿一定角度的入射波在某一频段内所引起的地表位移幅值比各向同性介质的可能要大,且随着孔洞埋深的增加,地表位移的幅值逐渐减小. 相似文献
2.
Based on the spherical cavity expansion theory in the elastic half space,the ground surface movement characteristics of shallowly buried explosions are analyzed.The results show that the induced seismic wave is a longitudinal wave in the near zone and a Rayleigh wave in the far zone.The maximum displacement(velocity) of the longitudinal wave and the Rayleigh wave are inversely proportional to the scaled distance,and can be described by exponential function with exponents equal to 1.4 and 0.5,respectively.The vibration frequencies of the waves have almost no change.The vibration frequency of the longitudinal wave approximates the natural vibration frequency of the cavity in the broken area,and the vibration frequency of the Rayleigh wave is about half that of the longitudinal wave.On the same reduced buried depth and reduced distance,the particle displacement is directly proportional to the product of the boundary loading and cavity radius,and is inversely proportional to the transversal wave velocity.Meanwhile,the particle velocity is directly proportional to the boundary loading and inversely proportional to the wave velocity ratio.In the far zone,the buried depth of the explosive only has a slight effect on the longitudinal wave,but has a larger effect on the Rayleigh wave. 相似文献
3.
地下衬砌洞室对Rayleigh波的放大作用 总被引：8，自引：3，他引：5
本文利用波函数展开法,给出了半空间中圆形衬砌洞室引起的Rayle igh波散射问题的一个级数解,并通过对求解技术的改进,获得了高频解答。论文通过入射波长、洞室直径、衬砌刚度和埋置深度等参数分析,研究了地下洞室对Rayle igh波的放大作用。数值结果表明,地下洞室情况下地表位移幅值可以达到无洞室情况的2.6倍以上;多数情况下,柔性衬砌洞室情况地表位移幅值最大,无衬砌洞室情况次之,刚性衬砌情况最小。 相似文献
4.
5.
SH波冲击下浅埋任意形孔洞的动力分析 总被引：5，自引：0，他引：5
求解了稳态SH波垂直弹性半空间水平表面入射时,浅埋任意形孔洞的动力响应。采用复变函数和多极坐标方法构造了一个能够自动满足水平表面上应力自由边界条件的散射波函数。应用这一波函数,将半空间中的问题转化为求解一个全空间中任意形孔洞的散射问题,最终将问题归结为对一组无穷代数方程组的求解。作为对抗爆问题的研究,给出了浅埋椭圆孔和方孔附近的动应力集中系数的数值结果,并对算例进行了讨论。 相似文献
6.
This paper presents a relationship between the focal depth in terms of Rayleigh-wave wavelength and the dominant frequency
of Rayleigh waves generated in a homogeneous half-space. Rayleigh waves were simulated using a (2, 4) staggered grid P-SV wave finite difference algorithm with VGR-stress imaging technique as a free surface boundary condition. VGR is an acronym
for vertical grid-size reduction. The simulated seismic responses using P-wave and SV-wave sources at different focal depths revealed Rayleigh-wave generation up to certain focal depth only for the considered
frequency bandwidth. A shift of normalized spectral shape of Rayleigh wave towards lower frequency with increasing focal depth
was inferred. Largest spectral amplitude was obtained in the wavelength for which the ratio of focal depth to the wavelength
of Rayleigh wave was around 0.17 in the case of P-wave source and 0.9 in the case of SV-wave source. An exponential decrease of spectral amplitude of Rayleigh wave with the departure of the ratio of focal depth
to Rayleigh wave wavelength from the above mentioned values was obtained. 相似文献
7.
The scattering and diffraction of plane SV waves underground, circular, cylindrical cavities at various depths in an elastic half space is studied in this paper. The cavities, studied here, are at depths of two to five cavity radii, measured from the surface to the center of the cavity. Fourier-Bessel series are used to satisfy the wave equation and the boundary conditions. When the angle of incidence of the plane SV wave exceeds the critical angle, surface waves are generated, which are expanded in terms of Fourier series, which also involve Bessel functions. The surface displacement amplitudes and phases that are presented show that the results depend on the following parameters: (1) The angle of incidence, θβ; (2) the ratio cavity depth to the cavity radius, h/a; (3) the dimensionless frequency of the incident SV wave, η; and (4) Poisson's ratio, v. The presence of the cavity in the half space results in significant deviation of both the displacement amplitudes and phases on the nearby half space surface from that of a uniform half space. 相似文献
8.
We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated. Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain. 相似文献
9.
扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果. 相似文献
10.
At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. In order to give our ‘measurements in the field’ a more quantitative interpretation, a series of experiments has been started under well-defined conditions. A cylindrical vessel contains water, simulating wet soil. A pulsed transmitting antenna (TA) is mounted above the water-surface irradiating horizontal underwater pipes. The reflected pulses are detected by a horizontal receiving dipole (AP, i.e. ‘air-probe’) in the air. The reflecting objects used in the experiments are steel pipes, PVC rods and hollow PVC pipes filled with air or water. The depth of the pipes varies from 0.25 to 2 m. The strength of the reflected pulse depends on the type of pipe, its diameter, its depth, the electromagnetic properties of the water and also on the strength and polarization of the incident E-field. The latter is (mainly) parallel to the axis of the pipe in the present experiments. The experimental results have been compared with calculated results using the theory of a continuous plane wave, incident on an infinitely long cylinder in a homogeneous dielectric medium with the same dielectric constant as water. In a previous paper an experiment was described in which a movable receiving dipole in the water measured the transverse, mutually orthogonal Eø- and Eθ-components of the pulses emitted by the TA. The amplitudes of Eθ versus depth, measured in that paper, are used in the calculations as amplitudes of the incident field. The attenuation of the scattered field is accounted for by assuming exponential damping. The relative amplitudes of the scattered field for different pipes lying at the same depth are reproduced very well by this simple theory. Also the relative values of the scattered amplitudes for a given pipe lying at different depths are described neatly, provided the refraction of the scattered field at the water-air interface is accounted for. 相似文献
11.
12.
The scattering of first mode linear baroclinic Rossby waves by a top-hat ridge in a continuously stratified ocean, with Brunt-Väisälä frequency that decays exponentially with depth below a surface mixed layer, is the subject of this study. A numerical mode matching technique is used to calculate the transmission coefficients for the propagating modes over the ridge. It is found that the scattered field depends crucially upon the stratification. For example, when the majority of the density variation is confined to a thin thermocline, corresponding to a small e-folding scale, gamma ?1, for the Brunt-Väisälä frequency, a large amount of the incident wave energy is reflected by a small amplitude ridge. Appreciable energy conversion between the propagating barotropic and baroclinic modes takes place in this case. An asymptotic analysis for a small amplitude ridge is presented that confirms these numerical results. In the limit gamma ?1→ 0, it is demonstrated that the scattered field in the continuously stratified ocean model differs markedly from the two-layer solution. The latter does not exhibit appreciable reflection of the incident wave energy for a small amplitude ridge. In conclusion, the application of a two-layer ocean model to describe Rossby wave scattering by ridges in place of a continuously stratified model cannot be recommended. 相似文献
13.
Numerical study on transient harbor oscillations induced by successive solitary waves 总被引：1，自引：1，他引：0
Tsunamis are traveling waves which are characterized by long wavelengths and large amplitudes close to the shore. Due to the transformation of tsunamis, undular bores have been frequently observed in the coastal zone and can be viewed as a sequence of solitary waves with different wave heights and different separation distances among them. In this article, transient harbor oscillations induced by incident successive solitary waves are first investigated. The transient oscillations are simulated by a fully nonlinear Boussinesq model, FUNWAVE-TVD. The incident successive solitary waves include double solitary waves and triple solitary waves. This paper mainly focuses on the effects of different waveform parameters of the incident successive solitary waves on the relative wave energy distribution inside the harbor. These wave parameters include the incident wave height, the relative separation distance between adjacent crests, and the number of elementary solitary waves in the incident wave train. The relative separation distance between adjacent crests is defined as the ratio of the distance between adjacent crests in the incident wave train to the effective wavelength of the single solitary wave. Maximum oscillations inside the harbor excited by various incident waves are also discussed. For comparison, the transient oscillation excited by the single solitary wave is also considered. The harbor used in this paper is assumed to be long and narrow and has constant depth; the free surface movement inside the harbor is essentially one-dimensional. This study reveals that, for the given harbor and for the variation ranges of all the waveform parameters of the incident successive solitary waves studied in this paper, the larger incident wave heights and the smaller number of elementary solitary waves in the incident tsunami lead to a more uniform relative wave energy distribution inside the harbor. For the successive solitary waves, the larger relative separation distance between adjacent crests can cause more obvious fluctuations of the relative wave energy distribution over different resonant modes. When the wave height of the elementary solitary wave in the successive solitary waves equals to that of the single solitary wave and the relative separation distance between adjacent crests is equal to or greater than 0.6, the maximum oscillation inside the harbor induced by the successive solitary waves is almost identical to that excited by the single solitary wave. 相似文献
14.
Analytical expressions for the borehole seismic signatures excited by an explosive source buried in an elastic half‐space 
下载免费PDF全文
下载免费PDF全文 Andrei M. Ionov 《Geophysical Prospecting》2018,66(Z1):161-174
Borehole guided waves that are excited by explosive sources outside of the borehole are important for interpreting borehole seismic surveys and for rock property inversion workflows. Borehole seismograms are typically modelled using numerical methods of wave propagation. In order to benchmark such numerical algorithms and partially to interpret the results of modelling, an analytical methodology is presented here to compute synthetic seismograms. The specific setup is a wavefield emanating from a monopole point source embedded within a homogeneous elastic medium that interacts with a fluid‐filled borehole and a free surface. The methodology assumes that the wavelength of the seismic signal is much larger than the borehole radius. In this paper, it is supposed that there is no poroelastic coupling between the formation and the borehole. The total wavefield solution consists of P, PP, and PS body waves; the surface Rayleigh wave; and the low‐frequency guided Stoneley wave (often referred as the tube wave) within the borehole. In its turn, the tube wave consists of the partial responses generated by the incident P‐wave and the reflected PP and PS body waves at the borehole mouth and by the Rayleigh wave, as well as the Stoneley wave eigenmode. The Mach tube wave, which is a conic tube wave, additionally appears in the Mach cone in a slow formation with the tube‐wave velocity greater than the shear one. The conditions of appearance of the Mach wave in a slow formation are formulated. It is shown that the amplitude of the Mach tube wave strongly depends on Poisson's ratio of the slow surrounding formation. The amplitude of the Mach tube wave exponentially decreases when the source depth grows for weakly compressible elastic media with Poisson's ratio close to 0.5 (i.e., saturated clays and saturated clay soils). Asymptotic expressions are also provided to compute the wavefield amplitudes for different combinations of source depth and source‐well offset. These expressions allow an approximate solution of the wavefield to be computed much faster (within several seconds) than directly computing the implicit integrals arising from the analytical formulation. 相似文献
15.
C.P. Sonett 《Physics of the Earth and Planetary Interiors》1975,10(3):313-322
A summary review of electromagnetic induction driven in the Moon by the interplanetary magnetic field is given. The point of view developed centers on inversion of Fourier transforms of the magnetic field in the free-stream solar wind (forcing function) and the response on the lunar surface measured by Lunar Surface Magnetometers. Conductivity profiles are shown to depend upon the central angle between the magnetometer given by local time and the incident wave-normal direction. The induction excites at least magnetic dipole and quadrupole “radiation”, but any scattered field is confined to the Moon's interior, save for propagation down the cavity where a TE-mode surface wave is generated. Confinement of the induced field on the sunward hemisphere and near the subsolar point is nearly complete, decreasing to the limbs, while in the diamagnetic cavity downstream of the Moon, partial confinement takes place. Both time and spatial multipoles of the induced field are present in the lunar interior complicating inversion into conductivity profiles. Profiles are reviewed and resolution limits are given and compared to those obtained from transient analysis. Finally a qualitative comparison to conductivity in the Earth is given. 相似文献
16.
— We develop a new scheme to compute 2-D SH seismograms for media with many flat cracks, based on the boundary integral method. A dry or traction-free boundary condition is applied to crack surfaces although other kinds of cracks such as wet or fluid-saturated cracks can be treated simply by assigning different boundary conditions. While body forces are distributed for cavities or inclusions to express scattered wave, dislocations (or displacement discontinuities between the top and the bottom surfaces of each crack) are used as fictitious sources along crack surfaces. With these dislocations as unknown coefficients, the scattered wave is expressed by the normal derivative of Green's function along the crack surface, which is called “double-layer potentials” in the boundary integral method, while we used “single-layer potentials” for cavities or inclusions. These unknowns are determined so that boundary conditions or crack surfaces are satisfied in the least-squared sense, for example, traction-free for dry cracks. Seismograms with plane-wave incidence are synthesized for homogeneous media with many cracks. First, we check the accuracy of our scheme for a medium with one long crack. All the predicted phases such as reflected wave, diffraction from a crack tip and shadow behind the crack are simulated quite accurately, under the same criterion as in the case for cavities or inclusions. Next, we compute seismograms for 50 randomly distributed cracks and compare them with those for circular cavities. When cracks are randomly oriented, waveforms and the strength of scattering attenuation are similar to the cavity case in a frequency range higher than k d $\simeq$ 2 where the size of scatterers d (i.e., crack length or cavity diameter) is comparable with the wavelength considered (k is the wavenumber). On the other hand, the scattering attenuation for cracks becomes much smaller in a lower frequency range (k d<2) because only the volume but not detail geometry of scatterers becomes important with wavelength much longer than each scatterer. When all the cracks are oriented in a fixed direction, the scattering attenuation depends strongly on the incident angle to the crack surface as frequency increases (k d>2): scattering becomes weak for cracks oriented parallel to the direction of the incident wave, while it gets close to the cavity case for cracks aligned perpendicular to the incident wave. 相似文献
17.
The full waveform synthetic seismogram of multiple scatteredSH waves by many cylindrical cavities in two-dimensional homogeneous elastic media is computed. We used the so-called single-layer potential integral representation of the scattered field and a discretization scheme with line source distribution for each cavity. The total field is the sum of the incident wave plus the field radiated from all sources, each multiplied by an unknown complex constant representing its strength. These constants are determined by imposing the appropriate boundary conditions in the least-squares sense. Here we solve scattering problems involving one, two, four, twelve and fifty cavities regularly distributed in a half-space. The seismograms computed along the free-surface show regions where the incident wave is strongly attenuated, as well as the arrivals of all multiple scattered phases. The accuray of the method is estimated from the degree of agreement of our solution for one cavity with the corresponding analytical solution, and also from the magnitude of the residual tractions along the boundaries of two cavities separated at various distances. Finally we apply the method to compute the case of fifty cylindrical cavities, each of radiusa, randomly distributed in a region 80a wide by 30a deep in a half-space. The value of scattering loss is obtained from the amplitude decay of the primary wave with distance for wavelengths in the range from 1.7a to 13.3a, using the synthetic seismogram calculated for the same distribution of 50 cavities as above, but in full-space. 相似文献
18.
Hydrodynamic constraints and storm wave characteristics on a sub‐horizontal shore platform 下载免费PDF全文
下载免费PDF全文 Few studies of wave processes on shore platforms have addressed the hydrodynamic thresholds that control wave transformation and energy dissipation, especially under storm conditions. We present results of a field experiment conducted during a storm on a sub‐horizontal shore platform on the east coast of Auckland, New Zealand. Small (<0.5 m) locally generated waves typically occur at the field site, whereas during the experiment the offshore wave height reached 2.3 m. Our results illustrate the important control that platform morphology has on wave characteristics. At the seaward edge of the platform a scarp abruptly descends beneath low tide level. Wave height immediately seaward of the platform was controlled by the incident conditions, but near the cliff toe wave height on the platform was independent of incident conditions. Results show that a depth threshold at the seaward platform edge > 2.5 times the gravity wave height (0.05–0.33 Hz) is necessary for waves to propagate onto the platform without breaking. On the platform surface the wave height is a direct function of water depth, with limiting maximum wave height to water depth ratios of 0.55 and 0.78 at the centre of the platform and cliff toe, respectively. A relative ‘platform edge submergence’ (water depth/water height ratio) threshold of 1.1 is identified, below which infragravity (<0.05 Hz) wave energy dominates the platform energy spectra, and above which gravity waves are dominant. Infragravity wave height transformation across the platform is governed by the relative platform edge submergence. Finally, the paper describes the first observations of wave setup on a shore platform. During the peak of the storm, wave setup on the platform at low tide (0.21 m) is consistent with measurements from planar sandy beaches, but at higher tidal stages the ratio between incident wave height and maximum setup was lower than expected. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
在Biot饱和土波动理论的基础上,用单排孔列作为屏障,对全空间饱和土中人射平面快纵波进行了隔离分析.采用波函数展开法,将人射波和散射波的势函数展开成Fourier-Bessel函数的无穷级数的形式,利用一组可以求解多散射问题的圆柱坐标系统和界面处位移连续的条件,得到问题的理论解.对影响隔离效果的孔间距和饱和土的渗透性等... 相似文献
20.
弹性层状半空间中无限长洞室对斜入射 平面SH波的三维散射(Ⅱ) 数值结果与分析 总被引：1，自引：1，他引：0
以基岩上单一土层场地为例, 计算分析了在斜入射平面SH波作用下弹性层状半空间中无限长洞室附近的地表位移. 研究表明, 层状半空间中地下洞室对波的散射与均匀半空间情况存在显著差别. 层状场地由于考虑了场地自身的动力特性, 使得洞室附近地表位移幅值的空间变化更为复杂, 基岩与土层刚度比、 土层厚度对散射效应均有着重要影响. 随着基岩与土层刚度比的增大, 地表位移幅值整体上逐渐增大; 随着土层厚度的增大, 土层对地表位移幅值的影响逐渐减小. 在频域解答的基础上, 给出了层状半空间中洞室对斜入射SH波散射的时域解答, 并以Ricker波为例进行了数值计算. 相似文献