基于f-x EMD波场分离的绕射多次波压制方法

自由表面多次波压制（SRME）方法是海洋地震数据处理中的关键步骤。当存在崎岖海底条件时,绕射波十分发育,常规SRME算法很难将复杂绕射多次波有效消除。对压制效果较差的多次波进行成因分析,首先基于f-x EMD波场分离算法,将全数据分解为反射部分和绕射部分,再对其进行多次波预测,获得绕射预测多次波子集：反射-绕射预测多次波、绕射-反射预测多次波和绕射-绕射预测多次波。然后,改进常规SRME算法的迭代过程,在SRME的迭代过程中使用与原数据中多次波更加匹配的绕射多次波子集进行多次波压制。与常规SRME算法的对比结果表明,改进的SRME算法在理论数据和实际数据处理中对多次波的压制都取得更好的效果。     

Use of P coda for determination of yield on nuclear explosions

Summary. Spectral analysis of P and P coda from NTS explosions recorded at NORSAR shows that magnitudes based on spectral integration of P coda provide a significantly more precise measure of yield than those based on the initial P . Coda amplitudes and spectra are considerably less sensitive to near-source variations than the direct P . A plot of P -coda magnitude, after correcting for the effects of attenuation in the upper mantle and the source spectrum based on knowledge of approximate yield or m_b , versus log yield appears to have a slope of unity. This implies that P coda from several stations (with different but known t * values) can be combined to improve the yield determinations further. Single station P coda from Soviet explosions can also be used to determine relative yields that appear to be at least as precise as those from Lg based on a network of stations.     

Wavefield downward extrapolation for migration velocity analysis on Marmousi data-set

The authors present a method for estimation of interval velocities using the downward continuation of the wavefield to perform layer-stripping migration velocity analysis. The generalized, phase-shift migration MG(F-K) in wavenumber-frequency domain was used for fulltime downward extrapolation of the wavefield. Such downward depth extrapolation accounts for strong changes of velocity in lateral and vertical directions and helps in correct positioning of the wavefield image in complex structures. Determination of velocity is the recursive process which means that the wavefield on depth level z _n−1 (n = 0, 1, ...) is an input data-set for determination of velocity on level z _n . The velocity ν [x, z _n − z _n−1] can be thus treated as interval velocity in Δz _n = z _n − z _n−1 step. This method was tested on synthetic Marmousi data-set and showed satisfactory results for complex, inhomogeneous media.     

Teleseismic waveform modelling with a one-way wave equation

It is now widely accepted that elastic properties of the continental lithosphere and the underlying sublithospheric mantle are both anisotropic and laterally heterogeneous at a range of scales. To fully exploit modern three-component broad-band array data sets requires the use of comprehensive modelling tools. In this work, we investigate the use of a wide-angle, one-way wave equation to model variations in teleseismic 3-D waveforms due to 2-D elastic heterogeneity and anisotropy. The one-way operators are derived based on a high-frequency approximation of the square-root operator and include the effects of wave propagation as well as multiple scattering. Computational cost is reduced through a number of physically motivated approximations. We present synthetic results from simple 1-D (layer over a half-space) and 2-D (subduction zone) models that are compared with reference solutions. The algorithm is then used to model data from an array of broad-band seismograph stations deployed in northwestern Canada as part of the IRIS-PASSCAL/LITHOPROBE CANOE experiment. In this region radial-component receiver functions show a clear continental Moho and the presence of crustal material dipping into the mantle at the suture of two Palaeo-Proterozoic terranes. The geometry of the suture is better defined on the transverse component where subduction is associated with a ∼10 km thick layer exhibiting strong elastic anisotropy. The modelling reproduces the main features of the receiver functions, including the effects of anisotropy, heterogeneity and finite-frequency scattering.     

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

Scattering parameters of the lithosphere below the Massif Central, France, from teleseismic wavefield records

Teleseismic P -wave recordings are analysed in the frequency range 0.3–6  Hz to derive structural (statistical) parameters of the lithosphere underneath the French Massif Central. For this we analyse differences in frequency-dependent intensities of the mean wavefield and the fluctuation wavefield. It is possible to discriminate a weak fluctuation regime of the wavefield in the frequency range below 1  Hz and a strong fluctuation regime starting above 1  Hz and continuing to higher frequencies. The observed wavefield fluctuations in the frequency range 0.3–3  Hz can be explained by scattering of the teleseismic P wave front at elastic inhomogeneities in the lithosphere. A statistical distribution of the inhomogeneities is assumed and the concept of random media is applied. The lithospheric structure under the Massif Central can be described as a 70  km thick heterogeneous layer with velocity fluctuations of 3–7 per cent and correlation lengths of the heterogeneities of 1–16  km.