首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
Damped Rayleigh-type waves at the interface of two solid media with strong density and wave velocity contrast have been studied theoretically and experimentally with a two-dimensional model; plexiglas and brass were used as model materials.  相似文献   

2.
为了研究弯曲界面曲率变化对分辨率的影响,首先推导了垂直入射下来自弯曲界面的反射波和透射波界面菲涅尔带近似的解析公式,证明了公式中曲率为零恰好对应已经被推导的平界面菲涅尔带的解析公式,然后给出了利用网格走时计算方法计算弯曲界面下反射波和透射波菲涅尔体的数值实现策略,这一实现策略同时保证了网格走时计算的精度和菲涅尔体计算的精度,最后对比了不同弯曲界面(不同曲率)下的菲涅尔体相对于平界面(曲率为零)下菲涅尔体的变化.研究结果表明,界面下高速时,向斜弯曲造成菲涅尔体在界面附近变宽,使得分辨率降低;背斜弯曲造成菲涅尔体在界面附近变窄,使得分辨率提高.并且向斜弯曲对分辨率的影响程度要明显大于背斜弯曲.而界面下低速时,结论正好相反.  相似文献   

3.
Summary An attempt was made to apply the Schlieren method for the investigation of kinematic and dynamic parameters of surface waves. The wave field generated by an exploding wire was studied in two-dimesional plexiglas models of a half space. Travel-time curves of the observed wave groups are given. Densitograms of surface waves obtained by microphotometric profiling parallel or perpendicular to the direction of the wave propagation enable to estimate the horizontal or vertical components of compressional stresses forming these waves.  相似文献   

4.
对于微地震正演模拟,本文以射线追踪的原理为基础,对两点间的射线追踪问题进行了研究,应用二分算法、改进二分算法和微变网格算法对水平层状匀速模型、弯曲层状匀速模型和复杂地质模型进行射线追踪,使得计算效率和适用范围都得到了很大的改善.文中对每种算法误差范围和计算效率进行了对比验证,对于不同的地质模型,选用合适的算法才能在计算速度和精度上得到双重保证,最后正演模拟了多波三分量记录.在模型建立上引入了超薄层概念,并在前人模拟的直达波、透射波、反射波基础上拟了折射波,使正演模拟的多波场信息更丰富.文中的应用实例及模型结果表明:与二分法相比,改进二分法能够对弯曲界面进行射线追踪,并能保证结果的精度.弯曲层状模型中,改进二分法与微变网格法相比计算速度有显著提高,能够应用到资料的反演中.  相似文献   

5.
6.
Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.  相似文献   

7.
In exploration geophysics, the efforts to extract subsurface information from wave characteristics exceedingly depend on the construction of suitable rock physics model. Analysis of different rock physics models reveals that the strength and magnitude of attenuation and dispersion of propagating wave exceedingly depend on wave-induced fluid flow at multiple scales. In current work, a comprehensive analysis of wave attenuation and velocity dispersion is carried out at broad frequency range. Our methodology is based on Biot's poroelastic relations, by which variations in wave characteristics associated with wave-induced fluid flow due to the coexistence of three fluid phases in the pore volume is estimated. In contrast to the results of previous research, our results indicate the occurrence of two-time pore pressure relaxation phenomenon at the interface between fluids of disparate nature, that is, different bulk modulus, viscosity and density. Also, the obtained results are compatible with numerical results for the same 1D model which are accounted using Biot's poroelastic and quasi-static equation in frequency domain. Moreover, the effects of change in saturation of three-phase fluids were also computed which is the key task for geophysicist. The outcomes of our research reveal that pore pressure relaxation phenomenon significantly depends on the saturation of distinct fluids and the order of saturating fluids. It is also concluded that the change in the saturation of three-phase fluid significantly influences the characteristics of the seismic wave. The analysis of obtained results indicates that our proposed approach is a useful tool for quantification, identification and discrimination of different fluid phases. Moreover, our proposed approach improves the accuracy to predict dispersive behaviour of propagating wave at sub-seismic and seismic frequencies.  相似文献   

8.
When a seismic source is placed in the water at a height less than a wavelength from the water–solid interface, a prominent S‐wave arrival can be observed. It travels kinematically as if it was excited at the projection point of the source on the interface. This non‐geometric S‐wave has been investigated before, mainly for a free‐surface configuration. However, as was shown in a field experiment, the non‐geometric S‐wave can also be excited at a fluid–solid configuration if the S‐wave speed in the solid is less than the sound speed in the water. The amplitude of this wave exponentially decreases when the source is moved away from the interface revealing its evanescent character in the fluid. In the solid, this particular converted mode is propagating as an ordinary S‐wave and can be transmitted and reflected as such. There is a specific region of horizontal slownesses where this non‐geometric wave exists, depending on the ratio of the S‐wave velocity and the sound speed of water. Only for ratios smaller than 1, this wave appears. Lower ratios result in a wider region of appearance. Due to this property, this particular P‐S converted mode can be identified and filtered from other events in the Radon domain.  相似文献   

9.
Kirchhoff-Helmholtz积分方法推广到远震转换波的合成地震图的计算,其优点是,能够计算复杂界面的地震波。通过与反射率法及动力学射线追踪的对比,表明KH积分方法能很好地模拟远震转换波震相,且精度较高。KH积分方法能够计算任意复杂界面的地震波,是一种研究地壳上地幔结构的有效方法。  相似文献   

10.
慢地震慢前兆的机制研究   总被引:12,自引:1,他引:12  
通过岩石高温高压破裂实验与有机玻璃破裂试验的应力—应变曲线对比、地震波形记录(观测资料)与试样破裂波形(实验资料)对比、地震前后定点应变与主破裂前后应变形态变化的对比,认为慢地震是临界或亚临界破裂或预先滑动所致,是低频应变波动。它是材料屈服、弱化或塑性变化的必然结果。而慢前兆则是临近大地震(大破裂)前出现的诸如形变、低频地震波(破裂弹性波)等的短期及临震前兆现象,它比较可靠。但由于许多岩石主破裂前并无明显的临界破裂或预滑动现象,也并非所有岩石都一定出现明显屈服,所以也并非每个地震前都有慢地震,因此,也就不一定都出现低频波动。从而,也并不一定都出现相同的短期、临震前兆现象。  相似文献   

11.
We consider the partial derivatives of travel time with respect to both spatial coordinates and perturbation parameters. These derivatives are very important in studying wave propagation and have already found various applications in smooth media without interfaces. In order to extend the applications to media composed of layers and blocks, we derive the explicit equations for transforming these travel–time derivatives of arbitrary orders at a general smooth curved interface between two arbitrary media. The equations are applicable to both real–valued and complex–valued travel time. The equations are expressed in terms of a general Hamiltonian function and are applicable to the transformation of travel–time derivatives in both isotropic and anisotropic media. The interface is specified by an implicit equation. No local coordinates are needed for the transformation.  相似文献   

12.
We report on an experimental study conducted to investigate the influence of small-scale wind waves on the airflow structure in the immediate vicinity of the air–water interface. PIV technique was used to measure the two-dimensional velocity fields at wind speeds of 3.7 and 4.4 m?s?1 and at a fetch of 2.1 m. The flow structure was analyzed as a function of wave phase. In the near-surface region, significant variations were observed in the flow structure over the waveform. The phase-averaged profiles of velocity, vorticity, and Reynolds stress showed different behavior on the windward and leeward sides of the wave in the near-surface region. The influence of wave-induced velocity was restricted within a distance of three significant wave heights from the surface, which also showed opposite trends on the windward and leeward sides of the crest. The results also show that the turbulent Reynolds stress mainly supports downward momentum transfer whereas the wave-induced Reynolds stress is responsible for the upward momentum transfer from wave to wind. In the immediate vicinity of the air–water interface, the momentum is transferred from waves to wind along the windward side, whereas, the momentum transfer is from wind to waves along the leeward side.  相似文献   

13.
We report on an experimental study conducted to investigate the influence of small-scale wind waves on the airflow structure in the immediate vicinity of the air–water interface. PIV technique was used to measure the two-dimensional velocity fields at wind speeds of 3.7 and 4.4 m s−1 and at a fetch of 2.1 m. The flow structure was analyzed as a function of wave phase. In the near-surface region, significant variations were observed in the flow structure over the waveform. The phase-averaged profiles of velocity, vorticity, and Reynolds stress showed different behavior on the windward and leeward sides of the wave in the near-surface region. The influence of wave-induced velocity was restricted within a distance of three significant wave heights from the surface, which also showed opposite trends on the windward and leeward sides of the crest. The results also show that the turbulent Reynolds stress mainly supports downward momentum transfer whereas the wave-induced Reynolds stress is responsible for the upward momentum transfer from wave to wind. In the immediate vicinity of the air–water interface, the momentum is transferred from waves to wind along the windward side, whereas, the momentum transfer is from wind to waves along the leeward side.  相似文献   

14.
When applying the conventional Fourier pseudospectral method (FSM) on a Cartesian grid that has a sufficient size to propagate a pulse, spurious diffractions from the staircase representation of the curved interfaces appear in the wavefield. It is demonstrated that these non-physical diffractions can be eliminated by using curved grids that conform to all the interfaces of the subsurface. Methods for solving the 2D acoustic wave equation using such curved grids have been published previously by the authors. Here the extensions to the full 2D elastic wave equations are presented. The curved grids are generated by using the so-called multiblock strategy which is a well-known concept in computational fluid dynamics. In principle the sub-surface is divided into a number of contiguous subdomains. A separate grid is generated for each subdomain patching the grid lines across domain boundaries to obtain a globally continuous grid. Using this approach, even configurations with pinch outs can be handled. The curved grid is taken to constitute a generalized curvilinear coordinate system. Thus, the elastic equations have to be written in a curvilinear frame before applying the numerical scheme. The method implies that twice the number of spatial derivatives have to be evaluated compared to the conventional FSM on a Cartesian grid. However, it is demonstrated that the extra terms are more than compensated for by the fewer grid points needed in the curved approach.  相似文献   

15.
We use Kirchhoff prestack depth migration to calculate migrated sections in 3-D simple anisotropic homogeneous velocity models in order to demonstrate the impact of anisotropy on migrated images. The recorded wave field is generated in models composed of two homogeneous layers separated by one either non-inclined or inclined curved interface. The anisotropy in the upper layer is triclinic. We apply Kirchhoff prestack depth migration to velocity models with different types of anisotropy: a triclinic anisotropic medium, an isotropic medium, transversely isotropic media with a horizontal (HTI) and vertical (VTI) symmetry axis. We observe asymmetry in migration caused by triclinic anisotropy and we show the errors of the migrated interface caused by inaccurate velocity models used for migration. The study is limited to P-waves.  相似文献   

16.
Equipment for simulation in laboratory conditions of hydrate-containing artificial samples and measuring their acoustic properties (wave velocities, absorption and attenuation) at different temperature and pressures is designed and constructed. The plant consists of a high-pressure chamber (up to 45 MPa), a measuring system intended for the excitation and reception of acoustic waves, systems for temperature and pressure control (axial and lateral) and for gas/liquid delivery into the sample. The measurements are performed on cylindrical samples with a 30-mm diameter and height of 10–50 mm. A set of successful test experiments was performed, including measurements of acoustic velocities of consolidated (plexiglas, sandstone, and frozen sand) and unconsolidated (dry and wet quartz sand) samples and formation of methane-hydrate bearing samples.  相似文献   

17.
A major complication caused by anisotropy in velocity analysis and imaging is the uncertainty in estimating the vertical velocity and depth scale of the model from surface data. For laterally homogeneous VTI (transversely isotropic with a vertical symmetry axis) media above the target reflector, P‐wave moveout has to be combined with other information (e.g. borehole data or converted waves) to build velocity models for depth imaging. The presence of lateral heterogeneity in the overburden creates the dependence of P‐wave reflection data on all three relevant parameters (the vertical velocity VP0 and the Thomsen coefficients ε and δ) and, therefore, may help to determine the depth scale of the velocity field. Here, we propose a tomographic algorithm designed to invert NMO ellipses (obtained from azimuthally varying stacking velocities) and zero‐offset traveltimes of P‐waves for the parameters of homogeneous VTI layers separated by either plane dipping or curved interfaces. For plane non‐intersecting layer boundaries, the interval parameters cannot be recovered from P‐wave moveout in a unique way. Nonetheless, if the reflectors have sufficiently different azimuths, a priori knowledge of any single interval parameter makes it possible to reconstruct the whole model in depth. For example, the parameter estimation becomes unique if the subsurface layer is known to be isotropic. In the case of 2D inversion on the dip line of co‐orientated reflectors, it is necessary to specify one parameter (e.g. the vertical velocity) per layer. Despite the higher complexity of models with curved interfaces, the increased angle coverage of reflected rays helps to resolve the trade‐offs between the medium parameters. Singular value decomposition (SVD) shows that in the presence of sufficient interface curvature all parameters needed for anisotropic depth processing can be obtained solely from conventional‐spread P‐wave moveout. By performing tests on noise‐contaminated data we demonstrate that the tomographic inversion procedure reconstructs both the interfaces and the VTI parameters with high accuracy. Both SVD analysis and moveout inversion are implemented using an efficient modelling technique based on the theory of NMO‐velocity surfaces generalized for wave propagation through curved interfaces.  相似文献   

18.
一、引言超声测井是一种新的地震測井法,它是利用超声波来研究小間隔內(1-2米)鈷井剖面岩层的物理机械性貭。近年来,由于石油探測事业发展的需要,这种方法在地球物理勘探和地球物理测井中受到了广泛的注意。根据超声測井的試驗和应用的結果看来,这种方法可以簡捷、經济、精确、詳細地測定出岩石的物理机械性貭(纵波和横波的波速,波的吸收系数等等)及其随深度的变化。  相似文献   

19.
Sandbars, submerged ridges of sand parallel to the shoreline, affect surfzone circulation, beach topography and beach width. Under time‐varying wave forcing, sandbars may migrate onshore and offshore, referred to as two‐dimensional (2D) behaviour, and vary in planshape from alongshore uniform ridges to alongshore non‐uniform ridges through the growth and decay of three‐dimensional (3D) patterns, referred to as 3D behaviour. Although 2D and 3D sandbar behaviour is reasonably well understood along straight coasts, this is not the case for curved coasts, where the curvature can invoke spatial variability in wave forcing. Here, we analyse sandbar behaviour along the ~3000 m man‐made curved coastline of the Sand Engine, Netherlands, and determine the wave conditions governing this behaviour. 2D and 3D behaviour was quantified within a box north and west of the Sand Engine's tip, respectively, using a 2.4‐year dataset of daily low‐tide video images and a sparser bathymetric dataset. The northern and western sides behaved similarly in terms of 2D behaviour, with seasonal onshore and offshore migration, resulting in a stable position on inter‐annual timescales. However, both sandbar geometry and 3D behaviour differed substantially between both sides. The geometric differences (bar shape, bar crest depth and wavelength of 3D patterns) are consistent with computed alongshore differences in breaker height due to refraction. The differences in the timing in growth, decay and morphological coupling of 3D patterns in the sandbar and shoreline are likely related to differences in the local wave angle, imposed by the curved coast. Similar dependency of bar behaviour on local wave height and angle may be expected elsewhere along curved coasts, e.g. shoreline sandwaves, cuspate forelands or embayed beaches. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

20.
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).  相似文献   

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

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