弹性波被动源数据一次波估计方法研究

与声波被动源数据相比,弹性波被动源数据含有更加丰富的地下地质信息.弹性波被动源数据合成的虚拟炮记录中,同样存在表面相关多次波,影响数据的后续处理以及对地下地质构造的准确判断.对弹性波被动源数据合成的虚拟炮记录进行多次波的预测与匹配相减,将比声波数据更加困难.本文提出弹性波被动源数据一次波估计方法,针对各向同性介质条件下的弹性波被动源数据进行研究,依据声波被动源数据一次波估计方法的研究思路,推导弹性波被动源数据一次波估计方法的理论过程,充分考虑弹性波被动源数据所包含的波场成分,包括纵波、横波以及在传播过程中产生的转换波.在传播过程中,由于能量的损失,略去多次传播的高阶项,进而给出相应的目标函数.利用数据的稀疏性,结合L_1范数约束的最优化问题求解方法,避免对弹性波虚拟炮记录进行多次波的预测与匹配相减,解决多次波影响问题.本文通过合成的弹性波被动源数据,进行一次波估计的测试,并与主动源和合成的虚拟炮记录进行对比,验证本文提出方法的有效性与准确性.     

关于弹性应力波的三种传播模式

分析了小变形情况下变形体的平动变形和转动变形,利用连续介质的动量和动量矩守恒方程,应用各向同性线弹性体的本构关系和偶应力本构关系,提出无限弹性介质中存在体积波(压力波)、旋转波(偶应力波)和偏斜波(偏斜应力波).旋转波和偏斜波均满足形式一致的四阶波动方程,但引起不同的运动学行为和应力状态,四阶旋转波和偏斜波的传播速度不再是依赖材料参数的常值.不计旋转变形时,各向同性弹性固体中只有体积波和二阶偏斜波的传播,这时传播速度均为常数.在3种波动模式中,体积波与传统弹性应力波理论完全一致,四阶偏斜波可退化为二阶偏斜波,但后者不同于传统应力波理论中以位移作为波函数的二阶剪切波.从变形运动学与内力的关系以及能量传播的角度分析了传统应力波理论中关于二阶旋转波和二阶剪切波存在的问题.     

海洋单阶次表面多次波分离方法研究

海洋地震资料普遍发育强能量表面多次波,传统表面多次波压制技术(SRME)能够预测出所有阶次表面多次波,但是各阶次表面多次波相互混叠.为了能够单独利用不同阶次的表面多次波成像,降低干涉假象对多次波成像的影响,需要将不同阶次的表面多次波分离出来.本文提出一种基于扩展SRME的海洋单阶次表面多次波分离方法.首先,应用SRME...     

波路径偏移压制层间多次波的理论与应用

消除层间多次波是地震勘探资料处理研究领域的难题,尤其对于实际资料的处理,到目前为止还很难找到一种完全有效的方法. 本文给出了仅对一次波成像既波路径偏移方法压制层间多次波方法,在共炮道集和共检波点道集分别计算炮点射线的入射角和检波点射线的出射角,由此计算的角度作为射线追踪的初始角度,计算地震波射线的传播路径. 结合由程函方程计算的走时表,判断偏移范围是反射波还是多次波. 在前期偏移过程压制多次波的理论研究基础上,本文主要研究波路径偏移消除多次波的应用部分. 为了进一步说明效果的有效性,计算了在单炮和共成像点道集压制层间多次波,给出了实际资料的压制多次波的偏移结果.     

大气中的内船舶波

通过波动方程的分析解求得稳定层结大气中孤立地形激发的大气内船舶波,证明这类波是由非拦截波和拦截波两部分构成,非拦截波中的辐散波和横波的形态,王要由地形形状和走向决定,它们都可向大气上层传输,而拦截背风波的波形主要由层结、风速及风速切变等大气条件决定.     

依据浅层剪切波速度估计纵波速度

对大庆测井资料进行了曲线拟合,对深井摆记录地方震地表反射波进行了到时差计算和波速测定,获得了松散沉积建造地层的纵横波速比。结果表明松嫩盆地及其周边地区的浅层波速比约为3.3左右。根据此波速比,推算得泊松比=0.45。参考Peter Borman给出的数据,所反映的地质构造应为新生代,与松嫩盆地的地质结构非常相近。因此可根据浅层工程勘探测定的剪切波速度来估计浅层纵波速度。对于临时在松散沉积层组建的流动台网,采用该速度模型,定位结果可有一定改进。     

基于AR模型的相干函数有效频段范围的确定

目前人们在进行地震动相干函数模型参数的拟合时,选取的频段范围存在很大差异,而确定有效的频段范围对相干函数模型拟合结果的合理性很重要。本文利用AR模型(Autoregressive Model,自回归模型)计算地震动的功率谱,提出了根据功率谱"信噪比(SNR)"确定相干函数有效频段范围的方法。本文的信噪比中的"噪",专指AR模型中的激励白噪声,用于衡量AR模型中某一频率地震动功率谱与激励白噪声功率谱之间的比例。定义信噪比SNR=0时对应的频率为有效频段的最高频率,大于最高频率的频率成分,被当作是地震波中的随机成分,可以不考虑地震动的相关性。采用这一方法,对SMART-1第45次地震记录的不同台站间距下的地震波进行了功率谱估计,根据信噪比确定出了有效频段范围,并计算了相应的相干函数,结果表明此方法是可行的。     

Application of the Spatial Auto-Correlation Method for Shear-Wave Velocity Studies Using Ambient Noise

Ambient seismic noise or microtremor observations used in spatial auto-correlation (SPAC) array methods consist of a wide frequency range of surface waves from the frequency of about 0.1 Hz to several tens of Hz. The wavelengths (and hence depth sensitivity of such surface waves) allow determination of the site S-wave velocity model from a depth of 1 or 2 m down to a maximum of several kilometres; it is a passive seismic method using only ambient noise as the energy source. Application usually uses a 2D seismic array with a small number of seismometers (generally between 2 and 15) to estimate the phase velocity dispersion curve and hence the S-wave velocity depth profile for the site. A large number of methods have been proposed and used to estimate the dispersion curve; SPAC is the one of the oldest and the most commonly used methods due to its versatility and minimal instrumentation requirements. We show that direct fitting of observed and model SPAC spectra generally gives a superior bandwidth of useable data than does the more common approach of inversion after the intermediate step of constructing an observed dispersion curve. Current case histories demonstrate the method with a range of array types including two-station arrays, L-shaped multi-station arrays, triangular and circular arrays. Array sizes from a few metres to several-km in diameter have been successfully deployed in sites ranging from downtown urban settings to rural and remote desert sites. A fundamental requirement of the method is the ability to average wave propagation over a range of azimuths; this can be achieved with either or both of the wave sources being widely distributed in azimuth, and the use of a 2D array sampling the wave field over a range of azimuths. Several variants of the method extend its applicability to under-sampled data from sparse arrays, the complexity of multiple-mode propagation of energy, and the problem of precise estimation where array geometry departs from an ideal regular array. We find that sparse nested triangular arrays are generally sufficient, and the use of high-density circular arrays is unlikely to be cost-effective in routine applications. We recommend that passive seismic arrays should be the method of first choice when characterizing average S-wave velocity to a depth of 30 m (V_s30) and deeper, with active seismic methods such as multichannel analysis of surface waves (MASW) being a complementary method for use if and when conditions so require. The use of computer inversion methodology allows estimation of not only the S-wave velocity profile but also parameter uncertainties in terms of layer thickness and velocity. The coupling of SPAC methods with horizontal/vertical particle motion spectral ratio analysis generally allows use of lower frequency data, with consequent resolution of deeper layers than is possible with SPAC alone. Considering its non-invasive methodology, logistical flexibility, simplicity, applicability, and stability, the SPAC method and its various modified extensions will play an increasingly important role in site effect evaluation. The paper summarizes the fundamental theory of the SPAC method, reviews recent developments, and offers recommendations for future blind studies.     

3D S-wave velocity structure of the Ningdu basin in Jiangxi province inferred from ambient noise tomography with dense array

The Ningdu basin, located in southern Jiangxi province of southwest China, is one of the Mesozoic basin groups which has exploration prospects for geothermal energy. A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration. In this study, we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution. Based on the dense seismic array including 35 short-period (5 s-100 ?Hz) seismometers with an average interstation distance of ～5 ?km, Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion. Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm. The results revealed obvious low-velocity anomalies in the center of the basin, consistent with the low-velocity Cretaceous sedimentary rocks. The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies. The obvious seismic interface is about 2 ?km depth in the basin center and decreases to 700 ?m depth near the basin boundary, suggesting spatial thickness variations of the Cretaceous sediment. The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basin-controlling fault, which may provide possible upwelling channels for geothermal fluid. This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.     

地震动空间相干函数模型中拟合参数的修正

采用SMART-1台阵第5次和第45次地震的水平向分量加速度记录,首先计算了不同间距台站对的地震动空间相干函数;然后讨论了台站距离对相干函数拟合结果的影响,即某一特定距离的相干函数与所有不同距离的相干函数的拟合结果存在明显差异。为减小这一差异,提出了一种对各个不同距离的相干函数拟合参数进行二次回归的方法,并选用Loh相干函数模型进行了验证,最后给出了基于Loh相干函数模型的拟合参数的修正结果。结果表明本文提出的修正方法将大大提高相干函数模型中参数的拟合精度。      

地脉动空间自相关方法反演 浅层S波速度结构

在近地表地球物理领域, 基于地脉动(或称背景噪声)提取的面波频散曲线反演地下S波速度结构是一种简单经济的工程勘察方法. 本文基于地脉动的空间自相关方法对一个微型台阵观测的背景噪声记录进行处理, 介绍了一种简单易行的提取频散曲线的数据处理方法, 获得了6.7—23 Hz频段的可靠频散曲线; 通过对该观测频散曲线与预测模型的频散曲线进行拟合, 反演得到S波速度结构. 结果表明, 该速度结构与钻孔直接测试的结果相吻合.      

不同散射模式下S波包络合成

为揭示地震波在地壳小尺度非均匀介质中的散射过程,更准确地描述地震波的包络展宽现象,本文基于多次各向异性散射理论,采用离散波数法求解能量密度积分方程,选取高斯型自相关函数表征的散射模式,得到S波能量密度包络.基于此,本文首先分析了单次散射和多次散射在形成S波能量密度包络中的贡献规律;然后探讨了吸收系数和总散射系数对合成S...     

Observed Evolution of Linear and Nonlinear Effects at the Dahan Downhole Array,Taiwan: Analysis of the September 21, 1999 M 7.3 Chi-Chi Earthquake Sequence

— In this paper, the site characteristics of the Dahan downhole array are studied by analyzing the September 21, 1999 M 7.3 Chi-Chi earthquake sequence including the main shock and some aftershocks. The four-level array (0 m, 50 m, 100 m and 200 m) is located to the north of Hualien City in eastern Taiwan. Polarization analysis is used to check the orientation errors of the seismometers at different levels of depth. If the surface instrument is chosen as reference, the angle between the major polarization axes of the surface and any downhole records is the orientation error that must be corrected for the downhole accelerographs. The orientation errors at depths of 50 m, 100 m and 200 m are 32°, 120° and –84°. After the corrections, the coherency between the surface and downhole records is substantially improved. Spectral ratio analysis shows that the predominant frequency of the Chi-Chi main shock shifts to a lower frequency. We also simulate ground motions at different depths by using the Haskell method with a linear velocity structure model. The record at surface is chosen as the input motion. Compared with the observed data, ground acceleration can be well reproduced for the aftershocks (weak-motion events) of the September 21, 1999 M 7.3 Chi-Chi earthquake. However, for the Chi-Chi main shock, the synthetic waveform cannot match well with the observation neither in amplitude nor in phase. This indicates that large ground shaking probably induced the nonlinear site effect at that time, and the model used cannot support it.Acknowledgement. The authors would like to express their thanks to Dr. L.F. Bonilla and one anonymous reviewer for their valuable suggestions. This research was supported by the National Science Council under grant number NSC 89-2921-M-194-007. The Institute of Earth Sciences, Academia Sinica supplied the strong-motion data. The support of these organizations is gratefully acknowledged.     

Surface and downhole shear wave seismic methods for thick soil site investigations

Shear wave velocity–depth information is required for predicting the ground motion response to earthquakes in areas where significant soil cover exists over firm bedrock. Rather than estimating this critical parameter, it can be reliably measured using a suite of surface (non-invasive) and downhole (invasive) seismic methods. Shear wave velocities from surface measurements can be obtained using SH refraction techniques. Array lengths as large as 1000 m and depth of penetration to 250 m have been achieved in some areas. High resolution shear wave reflection techniques utilizing the common midpoint method can delineate the overburden-bedrock surface as well as reflecting boundaries within the overburden. Reflection data can also be used to obtain direct estimates of fundamental site periods from shear wave reflections without the requirement of measuring average shear wave velocity and total thickness of unconsolidated overburden above the bedrock surface. Accurate measurements of vertical shear wave velocities can be obtained using a seismic cone penetrometer in soft sediments, or with a well-locked geophone array in a borehole. Examples from thick soil sites in Canada demonstrate the type of shear wave velocity information that can be obtained with these geophysical techniques, and show how these data can be used to provide a first look at predicted ground motion response for thick soil sites.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.