鄂尔多斯块体及其东南缘剪切波速度结构与波速比研究

利用深源远震的纯S波波形拟合方法研究了鄂尔多斯块体及其东、南缘的剪切波速度结构。结果显示，鄂尔多斯块体内部榆林台站下方速度结构相当稳定，没有很明显的高低速层变化，随着深度的增加，剪切波速度稳步升高；在该台下方25km深处有一个厚度达12km的层位(ys=3．90km／s)，该层的波速度相对其东部同等深度处的波速度偏低(以≥4．00km／s)。鄂尔多斯地块东部边缘的地壳厚度都比榆林厚，而且呈现高、低速层相间结构，低速层明显增多，与块体内部比较一致的特征为下地壳是一个厚达10km的速度较低的下地壳(yS=3．8km／s)，但埋深在35km左右。另外反演得到了各个台站的P波速度结构，并对相应的波速比(Vp/Vs)进行了研究，发现位于鄂尔多斯块体内部的榆林台波速比平均值仅为1．68，上地壳波速比处在非常低的状态(1．60左右)，中下地壳一直到上地幔的波速比稳定在1．73左右，反映介质均匀稳定，处于一种刚性状态；而鄂尔多斯块体东、南边缘的台站，都有几个波速比较高(1．80左右)的层位，其中东部边缘苛岚和离石地区的波速比平均值系统超过一般的弹性体波速比平均值1．732，反映了边缘地区地壳活动相对比较活跃，介质层中可能有其他物质的影响。     

Velocity changes associated with generalized triaxial deformation of pyrophyllite

Compressional and shear-wave velocities (V _p andV _s ) were measured during the generalized triaxial deformation (i.e. _{1 2}=2 ₃) of pyrophyllite. Observed velocity changes could be ascribed to crack development during dilatancy. Velocity changes were very localized with respect to the ultimate failure plane. The orientation and development of the failure plane was continuously observed with laser holography. Velocity reverals, i.e. changes from a decreasing trend to an increasing trend, were documented in a wet sample in bothV _p andV _s . These changes in bothV _p andV _p are inconsistent with dialatancy-diffusion models. The reversals were interpreted as a reflection of local stress reorientation caused by a slowly propagating fault.     

用尾波干涉法监测介质波速变化研究进展

Multiply scattered waves are sensitive to media changes owing to the effect of repeated sampling,superposition and amplification. Based on this characteristic,small-medium changes could be detected by using coda wave interferometry. In recent years,coda wave interferometry has been widely used in estimating velocity variation with high precision in areas such as seismology and non-destructive testing. This paper systematically presents the principle and research status of coda wave interferometry,and especial focus is placed on the research of media velocity variations by using repeating earthquakes,artificial sources,and ambient noise. Applications of coda wave interferometry can contribute to the more subtle understanding of dynamic evolution process in the medium.     

The Effect of Velocity Inversions on H/V

We analyzed the phenomenology of microtremor H/V curves under inversions in the shear-wave velocity (V_s) profile in the subsoil. Under no V_s inversion the spectral signature of the H/V peaks is found to be ‘eye-shaped’ with the horizontal components higher than the vertical. Conversely, under negative velocity gradients, numerous of differences emerge. I) A H/V ratio below 1 is observed for a wide range of frequencies, due to the decrease of the horizontal components below the vertical one. II) In the presence of persistent H/V < 1, small bumps in the H/V ratio given by local minima in the vertical spectral component may represent the relics of the peaks indicating resonances and stratigraphic discontinuities. As a consequence, in the presence of velocity inversions the H/V > 2 SESAME (2004) criterion fails but a stratigraphic interpretation may still be possible. III) The H/V curves should always be interpreted together with the single component spectra. IV) Microtremor H/V measurements for stratigraphic/microzonation purposes on stiff artificial soils, (asphalt, concrete, cement, pavements) should always be avoided since the latter often produce velocity inversions. This may have consequences in the intermediate to high frequency domain ( > 1 Hz) also in the application of reference site methods, like H_site/H_bedrock, to microtremor. Theoretical modeling confirms these experimental findings.     

Velocity anomalies: An alternative explanation based on data from laboratory experiments

Locations and velocities were calculated for microseisms occurring in samples of rock subjected to triaxial loading and injection of pore fluid. This was accomplished by analyzing arrival times of acoustic emission using an automatic first arrival picker. Apparent velocity anomalies were observed prior to both failure of intact samples and violent slip in samples containing saw cuts. Further analysis revealed that these fluctuations in calculated velocity were not due to changes in the true seismie velocity. Instead, variations in calculated velocity are shown to be related to sampling errors in picking first arrivals. The systematic picking of late first arrivals for small magnitude events was found to be a persistent bias resulting in low calculated velocities. This has encouraged the reexamination of earthquake records to determine how important sampling biases are in contributing to reported velocity anomalies.     

北京地区主要活动断裂深部速度结构特征及强震构造分析

根据分别由人工地震测深速度结构剖面数字化和天然地震层析成像建立的北京地区三维地壳速度结构模型,本文绘制了剖面图。同时将剖面两侧5km的地震震中投影到剖面上,分析了北京地区主要活动断裂的深部速度结构特征,证实两种方法得到的地壳速度结构模型具有很大的相似性。同时,顺义凹陷、马池口凹陷下方存在近垂直的深断裂,此深断裂可能为未来的发震断裂,地壳浅部与此对应的活动断裂是顺义-前门-良乡断裂北段、黄庄-高丽营断裂北段以及清河断裂,地壳内部的深断裂与地壳浅部控制盆地发育的活动断裂之间可能处于"汇而未交"的状态,具有地震孕育的深部构造背景。     

Velocity Structure Beneath Active Faults and Seismic Tectonic Characteristics in the Beijing Area

Velocity structure beneath active faults in the Beijing area has been discussed,based on the digital crustal model of velocity from deep seismic sounding profiles and tomography imaging of P waves. We found that there exists nearly vertical deep faults beneath the Shunyi Depression and the Machiko Depression,which are very likely to be seismogenic faults in the future. In the superficial crust,the north segment of the Shunyi-QianmenLiangxiang fault,the north segment of the Huangzhuang-Gaoliying fault and th...     

2003年大姚地震震中区的速度和衰减结构

使用2003年7月21日和10月20日云南省大姚县M62、M61地震震后流动数字化地震监测台网所记录的地震走时和波形资料,反演了震中区的Vp、Vs和Qp的三维结构.前人根据余震分布推测该震区存在一条北西西走向右旋走滑断裂带.三维结果表明,Vp和Qp的低值异常可以佐证该断裂带的存在;但除此断裂带,可能还存在一条北北东向断裂带,倾向西,表现为低Vp、Vs以及低Qp异常,并且异常深度达到8～10 km.因此,大姚震中区可能存在两条相互交错的断层,这与整个南北带断层的交错特征非常相似.另外我们推断该北北东向断带在南北方向可能有较长的延伸.     

九江—瑞昌序列震源位置及震源区速度结构的联合反演

采用震源位置和速度结构的联合反演方法确定2005年九江一瑞昌5.7级地震序列的分布和震源区的速度结构.结果显示九江一瑞昌地震序列存在NW330.方向和NE550方向的两组共轭分布,其中NW330.方向的优势分布更为明显;震源深度分布在3～18km的范围.地震主体破裂发生在上地壳-中地壳过渡层的P波速度低速区内,主震发生后余震向东南与西北方向两头扩展,其中往西北方向破裂时可能遇到了障碍体,积累的能量转向西南-北东方向继续破裂,引发了4.8级强余震.     

Velocity and concentration profiles of particle movement in sheet flows

Sheet flows occur widely in natural free-surface flows including rivers in flood, tidal estuaries and coastal waters in storm conditions when bed shear stress becomes sufficiently high. Particle volumetric concentration in sheet flows normally follows a linear distribution with the Rouse [Rouse H. Modern conceptions of the mechanics of fluid turbulence. Trans ASCE, 1937;102:463–543] distribution applicable in the dilute water column above the sheet-flow layer. However, a well-verified formula for determining particle velocity distribution in sheet flows is still lacking. This paper presents formulas to describe the particle velocity profile in steady or oscillatory sheet flows. They compare well with measured data. In particular, the novel formula for determining the particle velocity at the top of bedload–sediment-dominated sublayer in sheet flows is also well verified with measured data.     

Simultaneous Inversion of Earthquake Relocation and Velocity Structure in the Shanxi-reservoir,Wenzhou

Using the data of P-wave travel time recorded at the Shanxi-reservoir seismological network and Zhejiang and Fujian local networks, we implemented a simultaneous inversion of earthquake relocation and velocity structure and determined the new locations of earthquakes in the Shanxi-reservoir. The results show that: (1) the overall epicenter distribution is NW directed, and the Shanxi reservoir induced seismicity has a close relationship to the Shuangxi-Jiaoxiyang fault; (2) the focal depth of the Shanxi rese...     

基于微动技术探测盆地浅部地层速度结构1

本文介绍了利用微动技术探测地表浅层速度结构的方法,并在江苏常州地区一个凹陷盆地内进行了一次野外台阵观测试验,采集场地数据的质量良好.利用空间自相关法得到了不同半径下的自相关函数曲线,该场地的自相关函数显示,对于半径较大的台阵,其自相关系数较好;根据自相关函数提取了场地的频散曲线.通过设定初始模型,利用邻域反演算法得到了该场地的基本速度结构和各层深度分布范围;最后通过与场地附近的人工地震剖面进行对比分析表明,试验结果所揭示的土层结构与浅层人工地震勘探解释结果接近.     

高斯波包反射走时速度反演方法

扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.     

用河北张家口台速度与位移记录测定近震震级的对比

对河北张家口台FBS-3宽频带大动态范围的数字化速度地震仪和DD—2短周期模拟位移地震仪所记录的河北张北地区的地震，利用速度和位移数据分别测定了近震震级并进行了对比分析。结果显示，在一定的震级范围内，这两种资料测定的震级基本一致，没有发现二者之间存在系统误差，标准误差也较小。可以认为，无论从其物理意义还是具体做法上，用速度测定的震级与由位移测定的震级之间能够保持较好的一致性和连续性。     

Velocity model building from seismic reflection data by full‐waveform inversion

Full‐waveform inversion is re‐emerging as a powerful data‐fitting procedure for quantitative seismic imaging of the subsurface from wide‐azimuth seismic data. This method is suitable to build high‐resolution velocity models provided that the targeted area is sampled by both diving waves and reflected waves. However, the conventional formulation of full‐waveform inversion prevents the reconstruction of the small wavenumber components of the velocity model when the subsurface is sampled by reflected waves only. This typically occurs as the depth becomes significant with respect to the length of the receiver array. This study first aims to highlight the limits of the conventional form of full‐waveform inversion when applied to seismic reflection data, through a simple canonical example of seismic imaging and to propose a new inversion workflow that overcomes these limitations. The governing idea is to decompose the subsurface model as a background part, which we seek to update and a singular part that corresponds to some prior knowledge of the reflectivity. Forcing this scale uncoupling in the full‐waveform inversion formalism brings out the transmitted wavepaths that connect the sources and receivers to the reflectors in the sensitivity kernel of the full‐waveform inversion, which is otherwise dominated by the migration impulse responses formed by the correlation of the downgoing direct wavefields coming from the shot and receiver positions. This transmission regime makes full‐waveform inversion amenable to the update of the long‐to‐intermediate wavelengths of the background model from the wide scattering‐angle information. However, we show that this prior knowledge of the reflectivity does not prevent the use of a suitable misfit measurement based on cross‐correlation, to avoid cycle‐skipping issues as well as a suitable inversion domain as the pseudo‐depth domain that allows us to preserve the invariant property of the zero‐offset time. This latter feature is useful to avoid updating the reflectivity information at each non‐linear iteration of the full‐waveform inversion, hence considerably reducing the computational cost of the entire workflow. Prior information of the reflectivity in the full‐waveform inversion formalism, a robust misfit function that prevents cycle‐skipping issues and a suitable inversion domain that preserves the seismic invariant are the three key ingredients that should ensure well‐posedness and computational efficiency of full‐waveform inversion algorithms for seismic reflection data.     

用接收函数研究川滇地区国家地震台下地壳厚度及波速比

本文利用远震接收函数的方法,对川滇地区的昆明、腾冲、成都和攀枝花等4个国家地震台的台基下方不同方向的莫霍面深度及波速比进行了研究和分析。结果表明:昆明地震台台基下方的莫霍面深度基本在50km左右,波速比为1.62～1.69,地壳厚度和波速比不因方向不同而发生明显的变化;腾冲地震台台基下方的地壳厚度有着比较明显的方向性,东北方向厚为40.7km,东南方向为49.7km,两个方向的波速比相差也很大,差值达到0.2;成都地震台台基下方莫霍面的深度在40km左右,但是东北和西南方向要加深8km,两个方向波速比相差0.13;攀枝花地震台台基下方的地壳厚度比较稳定,厚度在60km左右,波速比变化也不明显。     

Precise Measurement of Subsurface Seismic Velocity Variation by Coda Wave Interferometry

A filed experiment was conducted continuously for three days, and the velocity variation was measured using coda wave interferometry. The measurement error is estimated to be around 10^-4, which coincides well with the theoretical error. The velocity variation during this period is up to 10^-3. The relationship between velocity variation and changes in air temperature, barometric pressure and solid earth tide was analyzed with linear least square fitting. The velocity has no dependence on air temperature, but displayed change of the order of 10^-6 - 10^-7 when the barometer or earth tide changed one Pa.     

利用远震接收函数反演黔江地震台下方的速度结构

本文选取黔江地震台2007年至2011年记录到的60个远震宽频带数字地震记录,采用频率域反褶积法获得台站的接收函数,并用H-Kappa叠加法来反演台站下方的地壳厚度和泊松比,最终得到了黔江地震台下方的地壳速度结构。同时,通过H-Kappa法反演得到的台站下方的地壳厚度作为波速反演的约束条件,以减少反演的非唯一性。计算结果显示,黔江地震台下方的地壳厚度为44km,这与广泛认同的中国大陆中西部地区莫霍深度在38—45km基本一致。本文对增强该地区深部地质构造特征研究和孕震机制分析具有积极的意义。     

环境激励下桥梁结构模态频率与波速变化监测研究

基于环境激励的桥梁健康监测方法是未来的发展方向。我们采用宽频带地震计Guralp CMG-6TD,建立了在线桥梁监测系统,对北京市昌平区鲁疃西路大桥进行了连续监测。分别使用随机子空间识别获得桥梁模态参数和尾波干涉法监测桥梁波速变化。结果表明：（1）采用的高灵敏度宽频带速度型地震计,能够连续准确记录宽频带范围内的三分量振动信号,非常适用于长期桥梁动态监测;（2）基于连续记录的环境激励信号,随机子空间法能够稳定可靠地识别桥梁低阶模态参数,尾波干涉法能够连续精确地监测桥梁材料的细微波速变化;（3）温度变化会引起桥梁材料的弹性模量发生显著变化。波速变化呈现振幅为1%的日变化特征,与温度变化呈明显的反相关性,温度敏感系数约为。同时模态频率的变化不太明显,说明波速可能是比模态频率更为敏感的量化损伤指标;（4）模态频率是桥梁整体状态的一种度量,而波速变化可以监测桥梁局部材料弹性性质的变化。将这两者结合起来,利用环境激励信号,连续准确识别桥梁模态频率以及监测波速变化,对于桥梁健康监测具有潜在重要应用价值。     

滇西地区地壳浅部基底速度细结构的研究

本文论述滇西地区思茅－中甸剖面基底导速度细结构的研究结果。基底层项部Ｐｇ界面深度为０－３．５ｋｍ，而底部Ｐ１＾０界面深度为１１．０－１７．０ｋｍ。基底层速度在金河－洱海断裂以南为５．７０－６．３０ｋｍ／ｓ，断裂以北增至６．３０－６．５０ｋｍ／ｓ，其过渡带位于剑川附近，剖面上断裂附近，除界面深度变化外，Ｐｇ面速度横向变化也较明显。速度等值线较为稀疏的景云桥炮南侧、大仓炮附近与支梯炮南侧地区，估计发生