共查询到20条相似文献，搜索用时 15 毫秒
1.
The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly. 相似文献
2.
Rayleigh面波的频散特性可以用来研究地表浅层结构. 本文使用时域有限差分法来模拟复杂黏弹性介质中的Rayleigh面波,研究了Q值对面波频散特性的影响.文中采用旋转交错网格有限差分,以非分裂卷积形式的完全匹配层为吸收边界,推出了求解二阶位移-应力各向同性黏弹性波动方程的数值方法.为了检验数值解的精度,首先将简单模型的正演结果与解析解对比,验证了方法的正确性;然后模拟了横向缓变层状介质和含有洞穴的介质中的面波,对弹性和黏弹性介质中的面波的频散特性进行对比分析.模拟结果表明浅层Q值对面波的频散特性有显著的影响;强吸收情况下,高阶面波的能量相对低阶面波能量显著增强. 相似文献
3.
We derive the phase velocity dispersion and the scattering for wave vertically propagating in a periodically weak‐contrast horizontally layered medium with arbitrary number of layers in a period. Phase velocity dispersion is defined as the frequency dependence of vertical travel time, and scattering is defined as a reflection coefficient at the interface between the multilayered system and the corresponding Backus medium. Low‐frequency approximation is used to define a dynamic effective medium with frequency‐dependent phase velocity. The results are compared with those obtained earlier for a gradient medium. We show that the low‐frequency weak‐contrast approximation is valid for models with realistic contrasts in elastic properties. 相似文献
4.
The presence of fractures in fluid‐saturated porous rocks is usually associated with strong seismic P‐wave attenuation and velocity dispersion. This energy dissipation can be caused by oscillatory wave‐induced fluid pressure diffusion between the fractures and the host rock, an intrinsic attenuation mechanism generally referred to as wave‐induced fluid flow. Geological observations suggest that fracture surfaces are highly irregular at the millimetre and sub‐millimetre scale, which finds its expression in geometrical and mechanical complexities of the contact area between the fracture faces. It is well known that contact areas strongly affect the overall mechanical fracture properties. However, existing models for seismic attenuation and velocity dispersion in fractured rocks neglect this complexity. In this work, we explore the effects of fracture contact areas on seismic P‐wave attenuation and velocity dispersion using oscillatory relaxation simulations based on quasi‐static poroelastic equations. We verify that the geometrical and mechanical details of fracture contact areas have a strong impact on seismic signatures. In addition, our numerical approach allows us to quantify the vertical solid displacement jump across fractures, the key quantity in the linear slip theory. We find that the displacement jump is strongly affected by the geometrical details of the fracture contact area and, due to the oscillatory fluid pressure diffusion process, is complex‐valued and frequency‐dependent. By using laboratory measurements of stress‐induced changes in the fracture contact area, we relate seismic attenuation and dispersion to the effective stress. The corresponding results do indeed indicate that seismic attenuation and phase velocity may constitute useful attributes to constrain the effective stress. Alternatively, knowledge of the effective stress may help to identify the regions in which wave induced fluid flow is expected to be the dominant attenuation mechanism. 相似文献
5.
6.
7.
High-frequency surface-wave analysis methods have been effectively and widely used to determine near-surface shear (S) wave velocity. To image the dispersion energy and identify different dispersive modes of surface waves accurately is one of key steps of using surface-wave methods. We analyzed the dispersion energy characteristics of Rayleigh and Love waves in near-surface layered models based on numerical simulations. It has been found that if there is a low-velocity layer (LVL) in the half-space, the dispersion energy of Rayleigh or Love waves is discontinuous and ‘‘jumping’’ appears from the fundamental mode to higher modes on dispersive images. We introduce the guided waves generated in an LVL (LVL-guided waves, a trapped wave mode) to clarify the complexity of the dispersion energy. We confirm the LVL-guided waves by analyzing the snapshots of SH and P–SV wavefield and comparing the dispersive energy with theoretical values of phase velocities. Results demonstrate that LVL-guided waves possess energy on dispersive images, which can interfere with the normal dispersion energy of Rayleigh or Love waves. Each mode of LVL-guided waves having lack of energy at the free surface in some high frequency range causes the discontinuity of dispersive energy on dispersive images, which is because shorter wavelengths (generally with lower phase velocities and higher frequencies) of LVL-guided waves cannot penetrate to the free surface. If the S wave velocity of the LVL is higher than that of the surface layer, the energy of LVL-guided waves only contaminates higher mode energy of surface waves and there is no interlacement with the fundamental mode of surface waves, while if the S wave velocity of the LVL is lower than that of the surface layer, the energy of LVL-guided waves may interlace with the fundamental mode of surface waves. Both of the interlacements with the fundamental mode or higher mode energy may cause misidentification for the dispersion curves of surface waves. 相似文献
8.
S波速度结构能够反映地球介质的物性差异,是地壳内低速区结构特征判别的重要依据.本文尝试利用空间自相关法(SPAC法)从地震台站微动信号的垂直分量中提取瑞利波相速度频散曲线,通过对频散曲线的反演获得地下介质的S波速度结构.以国家数字测震台网8个宽频带地震台站的实测微动数据为例,采用SPAC方法获得了首都圈地区北京附近约30 km 深度范围内的一维S波速度结构.结果表明,该区结晶基底埋深较浅约2 km;分别在5~8 km 和12~16 km 深处发育S波低速层;8 km 和 20 km 处是S波速度差异较大的速度分界面.这一结果与以往地震学及人工地震探测结果较为吻合,表明SPAC法估算地壳S波速度结构是可行、有效的. 相似文献
9.
THREE-DIMENSIONAL S-WAVE VELOCITY DISTRIBUTION BASED ON AMBIENT NOISE ANALYSIS IN EASTERN NORTH 
下载免费PDF全文
下载免费PDF全文 We apply ambient noise tomography to continuous vertical component broadband seismic data between January 1, 2010 and December 31, 2011from the regional networks of 190 stations deployed by China Earthquake Administration in Hebei, Shanxi and Inner Mengolia. Ambient noise cross-correlations were performed to produce the Green's functions of each station-pair. Firstly, we used the multiple-filter analysis method to extract surface wave group and phase velocity dispersion curves from inter-station paths at periods from 7 to 40s. Then the study area was discretized into a 0.2°×0.2° grid to obtain the group and phase velocity distributions using O'ccam inversion method. After that, three dimensional (3-D) S-wave velocity structures from the surface down to 50km are inverted from group and phase velocities dispersion results. the results of S wave velocity distribution maps generally demonstrate good correlations with surface geological and tectonic features, and they also clearly revealed the lateral velocity variation in the crust. In the mid-upper crust, the basins are clearly resolved with low S wave velocity due to its thick sedimentary layer, and the Taihang and Yanshan uplifts show relative higher S wave velocity distribution. With the increase of depth (>30km), the S wave velocity distribution presents a contrary characteristic compared to that of the shallow layer, and the S wave velocity beneath the Taihang and Yanshan uplifts are much lower than basin areas, which is possibly correlated with the thickness of the crust. 3-D S wave velocity shows a low-velocity zone at~10~20km depth observed beneath the Tanshan-Hejian-Xintai-Cixian belt and Bohai Bay. the low-velocity zone at~20~30km depth beneath the Datong area may be associated with the thermal material in the crust-mantle. Our S wave velocity distribution maps clearly show that Taihang Mountains is not only the boundary of topography and tectonic zone, but also the transition zone of high and low S wave velocity. 相似文献
10.
Site effects in Mexico City: Constraints from surface wave inversion of shallow refraction data 总被引：1，自引：0，他引：1
J. Ramos-Martínez F.J. Chávez-García E. Romero-Jiménez J.L. Rodríguez-Zúñiga J.M. Gómez-González 《Journal of Applied Geophysics》1997,36(4):157-165
In order to understand and simulate site effects on strong ground motion records of recent earthquakes in Mexico City, it is fundamental to determine the in situ elastic and anelastic properties of the shallow stratigraphy of the basin. The main properties of interest are the shear wave velocities and Q-quality factors and their correlation with similar parameters in zones of the city. Despite population density and paved surfaces, it is feasible to gather shallow refraction data to obtain laterally homogeneous subsoil structures at some locations. We focused our analysis in the Texcoco Lake region of the northeastern Mexico City basin. This area consists of unconsolidated clay sediments, similar to those of the lake bed zone in Mexico City, where ground motion amplification and long duration disturbances are commonly observed. We recorded Rayleigh and Love waves using explosive and sledgehammer sources and 4.5 Hz vertical and horizontal geophones, respectively. Additionally, for the explosive source, we recorded three-component seismograms using 1 Hz seismometers. We obtained phase velocity dispersion curves from ray parameter-frequency domain analyses and inverted them for vertical distribution of S wave velocity. The initial model was obtained from a standard first-break refraction analysis. We also obtained an estimation of the QS shear wave quality factor for the uppermost stratigraphy. Results compare well with tilt and cone penetrometer resistance measurements at the same test site, emphasizing the importance of these studies for engineering purposes. 相似文献
11.
Acoustic transversely isotropic models are widely used in seismic exploration for P‐wave processing and analysis. In isotropic acoustic media only P‐wave can propagate, while in an acoustic transversely isotropic medium both P and S waves propagate. In this paper, we focus on kinematic properties of S‐wave in acoustic transversely isotropic media. We define new parameters better suited for S‐wave kinematics analysis. We also establish the travel time and relative geometrical spreading equations and analyse their properties. To illustrate the behaviour of the S‐wave in multi‐layered acoustic transversely isotropic media, we define the Dix‐type equations that are different from the ones widely used for the P‐wave propagation. 相似文献
12.
用适配滤波频时分析技术处理了锡龙（SHIO）、清迈（CHTO）、昆明（KMI）和拉萨LSA）台记录的长周期数字化面波记录,获取了穿过缅甸弧及周边地区的530条路径的Rayleigh波频散,这些频散的周期范围为10.45～105.03 s.在此基础上,以分格频散反演方法从混合路径频散中提取了1°×1°网格内的纯路径频散,并且由网格内的纯路径频散反演出深达200 km的S波速度结构,最后重建了缅甸弧及周边地区的S波速度三维结构.所得结果表明：大致以实皆断裂为分界,其东部地壳波速较低,其西部地壳波速较高.印度-缅甸地区岩石圈厚度为110～130 km,上地幔顶部S波速度为43～4.4 km/s;而缅甸弧东侧的滇缅泰地块下方为一低速地幔柱上涌区,其宽度为150～200 km左右,这里的岩石圈厚度为70～80 km,上地幔顶部S波速度为41～4.2 km/s.另外,S波速度结构还反映出这一构造格局呈南北向的空间展布,并且与该区地震震源分布、断裂走向、火山分布有很好的对应关系. 相似文献
13.
Rayleigh wave imaging is efficient in estimating the shear- (S) wave velocity in near-surface exploration. The key is to accurately extract the dispersion of Rayleigh wave. We propose a method to calculate the dispersion of the active-source Rayleigh wavefield by using the Aki formulation. The spectrum after the cross correlation of two-channel records in the frequency domain is expressed by the Bessel function. Using the corresponding relation between the zero point of the spectrum real part and the Bessel function root, the phase velocity at the discrete frequency point is obtained and the dispersion curve is extracted. First, the theoretical basis and calculation method used in the active-source Rayleigh wave data are introduced. Then, three sets of theoretical models are calculated by this method and the results are consistent with the theoretical dispersion. Finally, we process a group of real Rayleigh wave data and obtain the phase velocity profiles and compared them with the results obtained by the multichannel surface wave analysis method. The effectiveness and applicability of the Aki method in active-source data processing are verified. Compared with multichannel wave processing, the advantage of the Aki method lies in the use of two-channel data in a single-shot record. When the number of acquisition channels in a shot gathers is insufficient or there is a bad channel, the quality of the extracted dispersion is guaranteed. 相似文献
14.
A nearly analytic exponential time difference method for solving 2D seismic wave equations 下载免费PDF全文
下载免费PDF全文 In this paper, we propose a nearly analytic exponential time difference (NETD) method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approximate the high-order spatial differential operators and transform the seismic wave equations into semi-discrete ordinary differential equations (ODEs). Then, the converted ODE system is solved by the exponential time difference (ETD) method. We investigate the properties of NETD in detail, including the stability condition for 1-D and 2-D cases, the theoretical and relative errors, the numerical dispersion relation for the 2-D acoustic case, and the computational efficiency. In order to further validate the method, we apply it to simulating acoustic/elastic wave propagation in multilayer models which have strong contrasts and complex heterogeneous media, e.g., the SEG model and the Marmousi model. From our theoretical analyses and numerical results, the NETD can suppress numerical dispersion effectively by using the displacement and gradient to approximate the high-order spatial derivatives. In addition, because NETD is based on the structure of the Lie group method which preserves the quantitative properties of differential equations, it can achieve more accurate results than the classical methods. 相似文献
15.
在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息. 相似文献
16.
裂缝广泛分布于地壳岩石中并具有多尺度的特点,研究地震波在裂缝介质中的传播规律具有重要的科学意义和应用价值.含流体的裂缝中传播的Krauklis波(也曾被称为慢流体波、裂缝波、斯通利导波等)近年来逐渐成为地球物理学界的研究热点之一.本研究利用地震物理模型技术研究横波在单裂缝样品中传播时随传播方向变化的响应特征,在实验室用环氧树脂为基质构建模型样品,一个为不含裂缝的参考样品,一个为含直径30 mm、厚度0.2 mm的薄币状单裂缝的模型,利用0.5 MHz换能器测试了横波在模型中传播时透射信号的特征.实验结果显示在单裂缝样品中透射横波发生了波形的畸变和振幅降低,振幅谱显示透射信号的主频明显降低,在初至横波之后为具有一定能量的尾波.对尾波进行时频分析,发现在方位角较小时尾波能量较弱而方位角较大时尾波能量较强,尾波中信号的主频在0.1 MHz左右,与前人理论和实验研究结果一致. 相似文献
17.
针对软地层套管井中弯曲波频散移向高频,以致超出了现行偶极声波测井仪器的激发频带的问题,本文在井外为孔隙地层时,采用Biot模型对套管井多极源激发的声波场进行了理论推导,对套管井偶极弯曲波的频散特性进行了数值模拟,考察了高密度水泥(或快速水泥)和低密度水泥(或轻质水泥)情况下的弯曲波频散,分析了不同水泥环对弯曲波主频散区和激发谱偏移的影响,重点对轻质水泥套管井中偶极源激发的模式波频散与激发谱及临界折射P、S波激发特性进行了研究,并考察了渗透率、孔隙度对弯曲波频散及衰减的影响,利用实轴积分法计算了偶极声源激发的时域全波波形.分析对比结果表明,在快速水泥情况下弯曲波频散曲线随着地层特征横波速度的减小会迅速向高频移动,随着快速水泥向轻质水泥变化,弯曲波频散曲线向高频移动将会减缓,对于特征横波速度低于1400m·s~(-1)的软地层,偶极弯曲波基础模式主频散区(或截止频率)可以由原来13kHz向低频移动至4kHz;在轻质水泥套管井中,无论是硬地层还是软地层,在目前偶极声波测井仪器声源主频激发下,接收波形中弯曲波均占主导地位. 相似文献
18.
应用三维非均匀交错网格有限差分程序模拟了多极子随钻声波测井的波场,考察了硬地层中单极子随钻声波测井响应及软地层中声源频率对多极子随钻声波测井的影响;计算了接收阵列波形的时域相干谱及频散特征,结果表明随钻测井仪器居中时,在硬地层中可以直接利用单极子声系得到地层横波速度,但地层纵波受到钻铤波的极大干扰.而在软地层中可以利用单极子声系直接测得地层纵波速度,但地层横波速度则需要使用工作在较低频率下的四极子随钻声波仪来测得,此时测井频段又会落入钻井噪声频带而受到影响.为此本文提出数值模拟手段结合物理实验的方法得到钻铤波的影响,在实际测井中通过信号处理方法而不是用在仪器本体上刻槽的方式来更完全地消除钻铤的影响.为了避免钻井噪声的影响,本文提出六极子随钻声波仪更适合软地层的横波测量.本文用数值实验结果来说明这些方法的可行性. 相似文献
19.
四川盆地中部浅层地壳一维剪切波速度结构初步研究 总被引：2，自引：0，他引：2
2010年1月31日四川盆地中部的遂宁发生了M5.0地震,四川盆地内部的台站观测到了很强的短周期面波波形数据。本文对地震数据的Rayleigh波和Love波提取基阶群速度频散曲线并反演得到四川盆地中部浅层地壳一维剪切波波速度结构。研究结果表明,四川盆地中部近地表剪切波平均速度约为2km/s,并且随深度逐渐增加,地壳深度在10km左右时,剪切波速度达到了3.5km/s。此结果适用于四川盆地中部以西,从遂宁到龙泉山附近,而四川盆地东部,从遂宁到华蓥山断裂不适用。该结果可为龙门山断裂附近的三维结构研究提供参考,并对强地面运动计算和区域内地震定位的研究有一定参考价值。 相似文献
20.
多极子阵列声波测井仪器采集的单极和偶极数据受到地层、井孔、仪器测量系统的影响.在处理实际声波测井数据时,必须考虑多极子模式波的频散效应,以及测井仪器在其中的影响.根据仪器等效理论和相位匹配方法,本文提出了一种从多极子阵列声波测井数据中同时获得纵、横波慢度的联合反演方法.这种方法的关键在于利用相同仪器-地层模型计算多极子模式波频散曲线,以此来匹配频域内纵波与横波数据的相位.相对于将泄漏纵波和弯曲波频散效应分开处理的其他方法,该方法不仅可以减少纵横波速度反演的不确定性,而且还避免了从声波数据中提取频散数据的繁琐过程.通过理论分析和现场数据处理证明了本文联合反演方法的准确性和有效性. 相似文献