地震各向异性——多组裂隙对横波偏振的影响

通过对多分量地震资料的分析,我们发现随着频率的增加横波分裂时差减小.对于深部接收的VSP数据来说快横波的偏振方向保持不变,而对于浅层接收的VSP数据来说偏振方向却存在一个最大可以达到20°的旋转.尽管多尺度随机分布微裂隙岩石物理模型已经成功地模拟并解释了横波分裂时差随频率变化的现象,却不能解释与频率相关的横波分裂.据推测,如果微裂隙的排列方向和大裂隙的排列方向不同,利用低频信息获得的偏振方向将指示裂隙主方向,而利用高频信息获得的偏振方向则指示微裂隙方向.在背景多孔隙介质中存在多组裂隙的情况下,推导出垂直入射条件下横波偏振方向的解析式,给出了系统研究横波在介质中传播的方法.研究结果表明,横波偏振方向会随着频率的变化而变化,并且在入射方位、角度一定的条件下,是裂隙方位和密度的函数,这些认识可能有助于揭示观测到的、依赖频率变化的横波偏振现象.     

昆明周边地区活动断层滑动与现代构造应力场

昆明盆地是一个受南北向活动断裂控制的新生代断陷盆地. 昆明盆地及其周边地区共有8条主要断裂带, 野外地震地质调查和断裂滑动观测表明, 南北向和北东向活动断裂以左旋走滑为主, 北西向活动断裂以右旋走滑为主. 采用野外观测取得的22个测点共706条活断层擦痕数据, 利用断层滑动资料反演构造应力张量的方法, 计算获得了研究区现代构造应力场主要特征: 昆明盆地及其周边地区现代构造应力场以北北西——南南东向挤压、 北东东——南西西向拉张为特征, 其最大主压应力（sigma;1）方向为333deg;——2deg;, 平均倾角21deg;; 最小主压应力（sigma;3）方向为44deg;——93deg;, 平均倾角14deg;; 中间主应力（sigma;2）倾角较陡, 近垂直. 应力结构均为走滑型. 断层滑动反演结果与震源机制解资料分析结果基本一致.      

地震照明分析及其在地震采集设计中的应用

地震照明分析是在给定地下背景结构的情况下用模拟方法研究采集系统对地下结构的探测能力.对目标照明的影响主要来自观测系统的设置,上覆地质结构的复杂性以及目标倾角等三项因素.我们将上述因素的影响统一纳入到对照明的定量计算中.本文阐述了关于地震照明的一系列基本概念,并基于波动理论给出对不同类型照明的计算方法.首先将入射和散射波分解到局部角度域并导出局部照明矩阵,它包括了在目标附近所有可能的入射和散射方向对照明的贡献.从照明矩阵出发,研究入射和反射波在角度域中与反射面的相互作用,并由此获得对地下结构的不同照明描述.作为例子,具体给出了对局部照明矩阵,空间体照明,反射面定向照明,反射面照明角度覆盖,以及对反射面可视性的计算方法,并给出了相应的数值计算结果. 最后,讨论了照明分析方法在地震采集设计中的一些可能应用.     

基于平面波模型重访地震背景噪声互相关及空间自相关(SPAC)

自Aki(1957)提出微震的空间自相关(SPatial AutoCorelation,SPAC)技术以来,SPAC技术一直独立发展,并在工程地震领域获得了广泛应用.近20年来,地震干涉(Seisimic Interferometry,SI)在多个领域引起人们的关注,该技术的核心思想是连续地震背景噪声的互相关函数(Noise Crosscorrelation Function,NCF)可以重建系统的格林函数(Green's Function,GF),对该技术的回溯性研究建立了SPAC和NCF的关系:它们是对同一物理现象的不同描述,SPAC在频率域中描述随机平稳噪声的空间相干,NCF在时间域中描述扩散场的互相关.理论上SAPC和NCF技术要求背景噪声源均匀分布,这样的噪声场可以用平面波叠加来模拟.本文基于平面波模型重访地震背景噪声的互相关和空间自相关技术,从单色平面波的互相关表示出发,对地震背景噪声互相关及空间自相关技术进行评述,试图使这些概念更易于理解.与之前众多研究地震干涉技术的理论相比,本文尤其关注以下几点:(1)基于简单的平面波模型,给出不同维度下,源或台站对方位均匀分布时,平面波互相关对入射波的方位平均和台阵对的方位平均结果,并对格林函数GF和时域互相关函数NCF的关系进行总结.(2)给出声源和(或)交叉台站方位分布不均匀时的互相关表示,指出这种非均匀性对方位的依赖关系,与弱各向异性介质中面波速度的方位依赖关系类似,因此,非均匀源的影响在反演时可能会映射到面波方位各向异性结果中.(3)互相关运算中,哪一个台站是虚拟源.NCF包含因果性和非因果性两部分,NCF的非对称性通常用于研究噪声源的方位分布,但由于源和接收的互易关系,及对互相关运算的不同定义和不同的傅里叶变换习惯,哪一个台站是虚拟源在目前的文献中并不明确.(4)方位平均和时间平均的关系.在SPAC处理中,需要对不同方位分布的台站对进行方位平均,本文从理论上说明,单个平面波入射时,交叉台站互相关系数对台站对的方位平均,等价于单个台站对互相关系数对入射波的时间平均.(5)几种特定分布非均匀噪声源的SPAC表示.包括单独的因果性噪声源和非因果性噪声源给出的互相关函数表示,及由此带来的相移问题.(6)利用SPAC、NCF和面波GF之间的关系,给出交叉分量的空间自相关系数表示.(7)衰减介质的空间相干表示.虽然利用地震干涉技术研究介质衰减在理论上仍然存在一些争议,但人们正试图研究从连续背景噪声记录中提取介质衰减的可能性.本文基于平面波模型,给出了不同坐标选择下,衰减介质的空间相干表示,这种表达的不同,指示了由地震干涉技术提取介质衰减的困难.与众多研究地震干涉的理论相比,比如稳相近似理论、互易定理、时间反转声学等,本文主要考虑均匀介质,不涉及非均匀介质的散射,从最简单的平面波模型,理解背景噪声重建系统格林函数这一地震干涉的核心思想和相应的基本概念.     

主动源和被动源面波浅勘方法综述

面波勘探在浅层勘探中获得了广泛应用.随着地震干涉理论和被动源噪声互相关技术的发展,发现传统的基于微动技术的空间自相关技术(SPAC)方法和背景噪声互相关函数(NCF)的物理本质是一致的,因此,将地震干涉理论用于被动源面波勘探中的技术和方法越来越受到关注.本文基于地震干涉理论和应用的发展,对面波勘探的原理进行了简要介绍,阐述了主动源方法中稳态法、瞬态法等各种常规方法.并根据最新的研究进展,对被动源SPAC和NCF方法在浅层勘探中的应用进行了评述,分析了各种勘探方法的原理和特点,叙述了主动源与被动源SPAC方法及NCF方法三者联合勘探的新思路,有望将传统的面波勘探技术拓展到浅层三维速度结构的精细探测.叙述了面波反演过程中多模式问题的一些重要概念和处理方法,最后讨论了面波方法的发展趋势和几个值得关注的研究方向.     

伪随机震源面波法频谱特征及其勘探深度分析

多道面波分析勘探法在工程场地安全评价中应用广泛,主要有人工源和天然源两种方法,但人工源法具有勘探深度相对较小的缺陷;天然源法其数据采集周期长现场操作便捷性差.文章通过引入伪随机震源作为面波勘探震源,希望能够在较短的采样周期内获取质量较高的频散能量谱.对获得的面波数据记录分别使用适合主动源法的F-K法以及适用于被动源法的SPAC法获取相应频散谱,并将伪随机震源的频散谱与同场地、同排列条件下人工源和天然源面波数据结果对比,认为:(1)运用SPAC法进行伪随机震源面波法频散分析所得到的频散谱收敛性较高;(2)伪随机震源面波法能够弥补人工源面波法勘探深度较小的缺点,相对于天然源面波法能够兼顾浅层高频区域,且现场操作简单、工作效率更高,适宜于不同地质条件的场地应用评价.     

推断勒夫波相速度的圆形台阵微震新技术

介绍了简单而又新颖的公式,使用微震的二分量,即圆形台阵记录的水平运动直接推断勒夫波相速度(cL)时,这些公式具有帮助作用。我们的公式与一般用于推断瑞雷波相速度(cR)的微震探测技术的空间自相关(SPAC)方法的公式类似。虽然SPAC方法现有的理论确实提供了估计cL的可能性,但这只能通过求解非线性的方程组才能实现,而其中待求解的未知数也包括cR以及瑞雷波?勒夫波的幂分配比。相反,在我们公式中cL是出现的唯一未知数。基于我们推荐公式——我们称为SPAC+L、SPAC-L和CCA-L(其中CCA表示无中心圆形台阵)法——的cL估算法的现场适用范围通过2个测试点的分析结果说明。在所提出的3种技术中,SPAC+L法由于表现最佳而出名。3种方法的有效波长范围下限是台阵半径r的2～5倍区域,上限是10～25r的区域。近年来,提出了类似的仅涉及cL的圆形台阵微震技术,但与它们相比,我们在此给出的方法,不是在数理计算上繁冗更少,就是在数学上更为简单。     

基于高通滤波的频率-空间域经验模态分解压制高频噪声

高频噪声压制是高分辨率地震数据处理中提高信噪比的关键性问题.本文针对f-x(频率-空间)反褶积空间预测滤波器无法处理非平稳、非线性信号的缺点,提出了一种基于高通滤波的频率-空间域经验模态分解(Empirical Mode Decomposition in the frequency-space domain,f-xEMD)压制地震剖面中高频噪声的方法.该方法采用全域高通滤波从原始数据中分离出含有部分有效信号的高频数据,将其变换到f-x域,然后在滑动的短窗口内提取每一个频率的空变数据序列进行EMD分解得到高频复本征模态函数(Intrinsic Mode Function,IMF)IMF1,将所有频率的IMF1序列反Fourier变换到时间域得到噪声剖面,将其与原始数据相减,达到高频噪声压制的目的.该方法可克服传统EMD分解方法中的模态混叠现象,保护陡倾角反射同相轴;压制后的噪声剖面中不包含有效信号能量,地震剖面的信噪比得到了提高.模拟数据和实际数据处理结果充分证明了该方法的有效性.     

SV波斜入射下凸起地形对地震动特性的影响分析

针对10种不同坡角的凸起地形,采用基于透射边界的有限元-有限差分计算方法,研究了局部凸起在SV波斜入射时,地震动峰值加速度放大倍数和反应谱谱比随入射角度的变化规律。结果表明:①地震波以一定的角度斜入射时,地表大部分观测点的峰值加速度放大倍数大于垂直入射的情况;②在计算模型宽高比一定的情况下,最不利入射角度与坡角有一定关系,且地表最不利位置随坡角增大由凸起台地边缘向中心移动;③入射角度对凸起地形地震反应谱特性的影响十分显著,不但影响谱比的幅值,也影响谱比曲线的形状,斜入射时各关键节点大部分周期点的谱比值大于垂直入射,入射角度对反应谱中的高频成分影响较为显著,而对长周期成分影响不大。     

地下衬砌洞室对Rayleigh波的放大作用

本文利用波函数展开法,给出了半空间中圆形衬砌洞室引起的Rayle igh波散射问题的一个级数解,并通过对求解技术的改进,获得了高频解答。论文通过入射波长、洞室直径、衬砌刚度和埋置深度等参数分析,研究了地下洞室对Rayle igh波的放大作用。数值结果表明,地下洞室情况下地表位移幅值可以达到无洞室情况的2.6倍以上;多数情况下,柔性衬砌洞室情况地表位移幅值最大,无衬砌洞室情况次之,刚性衬砌情况最小。     

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

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

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.     

Analysis of the spatial variation of seismic waves and ground movements from smart-1 array data

Specially designed arrays of strong motion seismographs located near earthquake sources are required for engineering studies of the near-source properties and the spatial variation of seismic waves. The SMART-1 array in Taiwan provides good records for this type of study. Careful study of the observed strong motion data permits the identification of wave types, directions and apparent wave velocities. In this paper, a principal direction ratio R (f,α) is defined; this indicates the principal direction of the motion (along a nearly straight line) within the range 0 < R < 1. Vertical motion of the ground is also included in this study. Orbit spectrum analysis is used to verify the identification of wave directions and wave types. The spatial variation of seismic waves along the principal direction is studied. From frequency-domain analysis, mathematical models of the spatial variation of ground displacement are developed using a wave-number spectrum and the cross-spectral density function between two spatial coordinates; these models in turn can provide two alternative models for the random vibration analysis of extensive structures subject to multiple point seismic excitation. The SMART-1 array data gathered during the January 29, 1981 earthquake also are used to demonstrate calculation of the ground strains and differential movements of the array site. From time-domain analysis, the spatial variation of seismic waves is defined for ground motion along the identified principal direction. The time variation of evolutionary spectra characterized by frequency-dependent parameters is used for this formulation. The SMART-1 array data again form the basis for discussion of the spatial variation of model parameters.     

Effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves

With the development of the dense array, the surface wave velocity and azimuthal anisotropy under the array can be directly obtained by beamforming the noise cross-correlation functions (NCFs). However, the retrieval of the Green's function by cross-correlating the seismic noise requires that the noise source has a uniform distribution. For the case with uneven noise source, the azimuthal dependence on the sources in the expression for the spatial coherence function, which corresponds to the NCF in the time domain, has the same form as the azimuthal dependence of the surface wave velocity in weakly anisotropic media. Therefore, the uneven noise source will affect the surface wave anisotropy extraction. In this study, three passive seismic methods, i.e., beamforming, SPAC (spatial autocorrelation), and NCF, are compared to demonstrate that an uneven source distribution and uneven station distribution have equivalent effects on the outcome from each method. A beamforming method is proposed to directly extract the velocity and azimuthal anisotropy of surface waves. The effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves was investigated using data from the ChinArray Phase II. A method for correcting the apparent anisotropy in beamforming results caused by an uneven station distribution is suggested.     

四川宜宾地区S波分裂特征

本文采用纵横比与偏振分析相结合的方法测定了2013年4月25日~2015年12月31日四川宜宾地区10个台站S波分裂参数,即快波偏振方向和慢波延迟时间。结果表明,华蓥山断裂两侧台站呈现不同的快波偏振优势方向,断裂带以西的台站偏振优势方向为NW向,与区域应力场方向一致;位于断裂带以东的台站优势偏振方向为NE向,与断裂走向一致。在地震密集分布区域内的CNI台的优势偏振方向为NE向,与台站附近的断裂带走向基本一致。研究区域南段的3个台站(JLI、YAJ、XWE)优势偏振方向近NS向。各个台站平均慢波延迟时间在3.07~11.95ms/km范围内,慢波延迟时间最大的台站是CNI台,距离2013年4月25日06时10分M_L5.2地震震中位置最近,这反映出震源区地震波各向异性程度较强。CNI台站的慢波延迟时间显示,在2015年2月7日M_L4.8地震前观测到慢波延迟时间有明显的上升趋势。     

Propagation characteristics of blast-induced shock waves in a jointed rock mass

The propagation characteristics of blast-induced shock waves in a jointed rock mass have been monitored and studied. Accelerometers were set up on a rock surface along three lines, at 0°, 45° and 90° with respect to the orientation of the predominant joint strikes. Cylindrical charges were detonated in a charge hole, and ground accelerations in both vertical and radial directions at various points on the rock surface were recorded. Results show that rock joints have significant effects on the propagation characteristics of blast-induced shock waves. The amplitude and principal frequency of shock waves attenuate with the increase of distance from the charge centre, and the increase of incident angle between the joint strike and the wave propagation path. The measured data were compared with the empirical equations of shock wave attenuation proposed by other authors. The mechanism of rock joint effect, the attenuation of shock waves in relation to the propagation distance, the charge weight and the incident angle, are discussed in this paper.     

ICEARRAY: the first small-aperture, strong-motion array in Iceland

A small-aperture, strong-motion array, the ICEARRAY, has been deployed in South Iceland, a region with a history of destructive earthquakes, some exceeding magnitude 7. The array’s purpose is: (1) monitoring future significant events in the region, (2) quantifying spatial variability of strong-motion over short distances and (3) shedding light on earthquake source processes. The number of array stations and their arrangement were based on an optimisation of the shape of the corresponding array transfer function (ATF). The optimal ICEARRAY configuration comprises 14 stations, has an aperture of ~1.9 km and a minimum interelement distance of ~50 m and possesses a near-azimuthally independent ATF with a sharp main lobe, negligible sidelobes and a wavenumber range of 1.5–24 rad/km. Accordingly, the ICEARRAY has the intended capabilities of capturing seismic waves in the frequency range of 1–20 Hz, which is of main interest to earthquake engineering and engineering seismology applications.     

微动探测方法软件的实现

利用微动进行勘探无需振动源,方便快捷、对环境没有干扰,是一种很有发展前景的勘探方法。针对微动探测数据处理的需求,开发研制了一套可视化微动探测处理分析软件系统MSMSPAC。该软件采用C~#语言开发,主要功能包括台阵设计、台阵方位角分布分析、数据格式转换、微动有效数据段挑选、数据滤波、重采样、参数计算及频散曲线提取、编辑、合并、项目管理等功能。软件系统界面友好,操作简便。     

Paleomagnetism and Cyclostratigraphy of the Middle Ordovician Krivolutsky Suite,Krivaya Luka Section,Southern Siberian Platform: Record of Non-Synchronous NRM-Components or a Non-Axial Geomagnetic Field?

The Middle Ordovician Volginsky and Kirensky fossil zones were sampled in the Krivaya Luka section (Krivolutsky suite) that outcrops along the Lena river in Siberia. The Volginsky and Kirensky zones are coeval to the Llandeilo in the global geologic time scale. The Krivaya Luka section consists of siltstones, clays, sandstones, and limestones, and displays a remarkably distinct sedimentary cyclicity, especially in its reddish middle part.Stepwise thermal demagnetization yields three NRM components. Component A, isolated in the 100—250°C interval can be either normal or reversed. The normal A-component has a direction close to recent local magnetic field. The reversed A-component directions are scattered around a direction close to that of the lower Triassic Siberian traps. Component B has unblocking temperatures that range from 400 to 500°C and is represented mainly by normal polarity directions. The B-component, isolated from rocks of the middle part of the section is of a normal polarity with D = 176.5°, I = 30.0° and a North pole position at 16.2°S, 111.3°E. The other parts of the section are characterized by intermediate B-directions, which resulted possibly by partially overlapping A- and C-components. The highest temperature dual-polarity component C was isolated in the 550—670°C interval, resulting in the detection of two complete polarity zones and three magnetic reversals. The C-component is characterized by the following mean directions: for the reversed component D = 335.7°, I = 6.9°, and for the normal component D = 188.6°, I = 28.0°, which is very close to the normal polarity directions of the B-component. The corresponding paleomagnetic North pole for reversed polarity rocks is 32.6°S, 137°E, which is typical of Middle Ordovician rocks from Siberia – the mean pole for Llanvirn-Llandeilo is 30°S, 136°E (cf. Smethurst et al., 1998) – whereas for normal polarity rocks the pole position 17.2°S, 99.1°E is markedly different. Nevertheless, we assume that the C-component records the ancient geomagnetic field of Ordovician times, even though it does not pass the reversals test. This could be explained by overlapping NRM unblocking temperature spectra for the B and C components. In this case, the paleomagnetic pole positions should be interpreted with some caution.In addition, the section was logged and sampled in detail for cyclostratigraphic purposes. Spectral analysis in the depth domain using the high-field susceptibility as input parameter showed that the observed cyclicity is most likely orbitally forced. Detected spectral peaks (significant at the 95% confidence level) were close to the expected positions of the periodicities of precession, obliquity and eccentricity for the Ordovician. Consequently, the average sediment accumulation rate is estimated at 3.5 cm/kyr. Extrapolating this sedimentation rate yields a total duration of at least 1 Myr for the Volginsky fossil zone and 1.2 Myr for the entire Krivaya Luka section.