岩石破裂孕育过程中的几种Q值变化

采用脉冲传播谱振幅比法、能量法和尾波法分别测定片麻岩破裂孕育过程中的Ｑ值变化．结果表明:谱振幅比法所得Ｑ值呈现相对平稳值——升高至甚高值(极大值)——回降——破裂的变化形式；能量法所得Ｑ值变化形态与上述相反，即平稳值——低值——回升——破裂；而尾波法所得Ｑ值，主频Ｑ值呈现相对高值——下降至低值(极小)——回升——破裂，与能量法Ｑ值变化大体一致．而对于高频(大于主频)Ｑ值则呈现平稳值——升高至极大值——回降——破裂．同时尾波振幅谱还表明:应力达70％破裂应力后，首先是高频Ｓ波迅速减少(QC)变大)，临近主破裂时主频Ｓ波则大幅度减少(QC回升)，这可能是两种尾波QC变化形态不同的原因所致．而P，S波(初至波)振幅谱也表明，随着裂隙出现，Ｓ波频带变窄，高频成分迅速减少．即S波谱变化不同．这使采用不同方法所得的QS的变化特征不同．      

负坡段加卸载响应比的物理意义及预测效能

根据实验结果，用应变能研究岩石在极值应力后主破裂前（负坡段）加卸载响应比可能出现奇异点或负值的变化特征，并且说明其普遍性．以应变（能）及Ｐ波走时为响应，分别得出加卸载响应比的倒数在岩石破裂孕育过程中的变化形态．结果表明:在逼近主破裂时，加卸载响应比的倒数急剧下降后出现极小值（甚至为负值）．此特征可以看作临近主破裂的标志．      

尾波振幅比和能量比异常特征研究

利用中国数字地震台网（CDSN）兰州台记录的短周期数字地震资料，测定了尾波最大振幅与S波最大振幅之比r1、尾波能量与S波能量之比的平方根值r2随时间的变化。结果表明：在永登主震（19950和天祝-古浪主震（1996）前，r1和r2均出现上升或下降异常变化，且有较大的偏离平均值的异常高值或异常低值出现。     

云南地区地震序列早期地震波动力学特性

研究了云南地区各类型地震序列的早期地震波动力学特征－地震波垂直分量Ｐ波初动和Ｐ、Ｓ波最大振幅比，结果表明：前震序列振幅比为过原点的直线，比值稳定、Ｐ波初动一致；双震型地震序列在第一主震后Ｐ波实动出现一组时，仍有一类波形的振幅比保持为过原点的直线，第一主震后Ｐ波初动出现两组时，仍有初动相同的一组地震振幅比保持过原点的直线，或者两组初动的地震振幅比分别、共同成过原点的直线；主余震地震序列在非大震孕育期     

岩石临破坏前波速奇异变化的再研究

从实验上确证：岩石在主破裂前一瞬间，纵波速度升高甚至超过原最大值的两倍以上，纵、横波速比Vp/Vs升高甚至达2.5～3.0，之所以造成这种现象：对前一个原因是主破裂前有些张裂缝闭合，另一个原因是主破裂前微破裂急剧增多，在某一时刻同时发生几个小破裂，各检波器所记讯息可能不是同一事件。对后，P、S波的走时在破裂孕育过程中的下降--回升（相对变化量）并不同步而是P波的变化先于或大于S波也是重要原因。     

脉冲超声波穿过单轴压缩岩样的动态传播特性

本文叙述脉冲超声穿过单轴压缩岩样时的瞬态波形变化特征。1.大约从破裂应力的80％起,脉冲超声波形的衰减系数K开始下降,并于主破裂前有一显著低值;2.脉冲超声第一波峰的幅值Ap_1从破坏应力的70％,走时t_(p1)从破坏应力的80％开始,随载荷应力σ的增加幅值减小,走时增加;3.最大振幅与其持续时间的对数比α和波形衰减系数K显示同步变化。据此认为,上述变化特征是单轴压缩岩样对脉冲超声的综合影响。     

利用祁连山主动源资料研究2016年门源6.4级地震前后波速变化

采用叠加、互相关、插值拟合等方法对祁连山气枪主动源台网数据进行处理,针对2016年1月21日青海门源6.4级地震前后气枪激发地震P波、S波震相的走时变化特征进行了分析,结果表明,地震前约6个月时,震中附近3个台站的相对走时出现下降变化（走时减少）,至震前约3个月时低值异常恢复正常,之后再次出现走时下降变化,地震即发生于走时变化恢复过程中。S波走时变化最大下降幅度达18ms,震后走时变化逐渐恢复正常,且3个台站变化趋势较为一致,其中,距震中最近的台站的S波走时变化最明显（ZDY38台）,较远台站的走时变化幅度较小,其变化特征与震源区位置有关。走时缩短意味着速度增加,可能与区域应力积累间存在一定的关系。     

云南武定地震余震近场尾波衰减的区域特征

根据地方震尾波散射理论，利用云南武定MS6.5主震后近场田心和云龙两个流动数字化台站记录的余震序列地震资料，以不同中心频率（1.5～20Ｈｚ）取样分频研究了近场尾波振幅衰减率（f ）的空间分布特征．结果表明，两个台站给出的尾波振幅衰减率（f ）存在明显的小区域范围非均匀性变化．位于主发震断层两侧的云龙台（东侧）和田心台（西侧）测得的尾波振幅衰减率（f ）值，当 f ＜6.0Ｈｚ时基本一致．相应的Qc（f ）值为66～120的低值；但在高频范围（f 6.0Ｈｚ），由于穿过主断层的近场尾波遭到断层破碎带强烈的吸收，因而云龙台测得的尾波振幅衰减率（f ）值，比不穿过主断层的田心台给出的（f ）降低30％，并讨论了产生近场尾波振幅衰减空间分布差异的原因．本文的研究结果还表明，波源因子A0（f）不仅与震源强度相关，也与频率犳相关，用1.5～20Ｈｚ求出的Ｓ波平均自由程犔在主断层东侧为30ｋｍ、西侧为40ｋｍ ，由此证实了地震波穿过主断层破碎带时有更强的散射．      

差应力条件下岩石Q值测定方法的改进

本文指出了脉冲传播和谱振幅比法测定岩石破裂孕育过程中Q值变化矛盾的原因,提出用能量法测定Q值,并以大理岩,大姚铜矿砂岩为例给出了真三轴压缩时岩石破裂孕育过程中P波总能量的变化及不同方法所测Q值的变化特征。其中能量法所测Q值与P波总能量变化特征一致,即主断裂前它们都出现了“下降—回升”变化及变化的各向异性。     

地震波形群体特征在强余震预测中的应用

研究了华东地区１５个地震序列的振幅比、初动方向和初动半周期资料。得到：大震（主余型序列中的主震或震群型序列中最大地震）前后的Ｐ波初动方向，除人们熟知的由优势分布转化为紊乱外，还有可能呈优势分布，优势方向改变的特征。余震序列的Ｐ波初动方向和振幅比特征可能有助于预测强余震。初动半周期特征可能用于区别的前震和余震。     

The Q value variations in the preparing pro－cess of rock rupture

ＴｈｅＱｖａｌｕｅｖａｒｉａｔｉｏｎｓｉｎｔｈｅｐｒｅｐａｒｉｎｇｐｒｏ┐ｃｅｓｏｆｒｏｃｋｒｕｐｔｕｒｅＢＩＮＷＡＮＧ１）（王彬），ＺＨＡＯ－ＹＯＮＧＸＵ１）（许昭永），ＪＩＮ－ＭＩＮＧＺＨＡＯ１）（赵晋明）ＹＩ－ＬＩＨＵ２）（胡毅力）１）Ｓｅｉ...     

TheQ value variations in the preparing process of rock rupture

In this paper, the "spectral amplitude ratio method" (SAR), "energy method" (EN) and "coda wave method" (CW) are used to calculate theQ value variations of gneiss in the preparing rupture process. The obtained results show that the variation state ofQ values by SAR features the shape of relative stability—gradual increment to the maximum—then decrement and final rupture. The variation state ofQ values by EN is just contrary to that by SAR, i. e. with the shape of stability—decrement—increment—and final rupture. The varation state ofQ values by CW is similar to that by EN, its main frequency features the shape of relatively high value—decrement to the minimum—increment—and final rupture. But to the high frequency (higher than the main frequency), the variation state ofQ values features the shape of the stable value-increment to the maximum-decrement-and final rupture. At the same time, the results by coda wave amplitude spectrum show that, when stress reaches 70% of rupture stress, the high frequency component of S wave rapidly reduces (Q _c increasing); at the time of impending the main rupture, the main frequency component reduces with a large scale (Q _c increasing again), this may be the reason which causes the different variation states of two codaQ values. The result of amplitude spectra of P, S (initial wave) waves also show that with the appearance of microcracks the frequency band of S wave turn to be narrow, the high frequency component is reduced quickly, i. e. the S wave spectra have different variation states with different frequency components. That is why theQ _s obtained by different methods have different variation characteristics.     

The Q value variations in the preparing process of rock rupture

In this paper, the "spectral amplitude ratio method" (SAR), "energy method" (EN) and "coda wave method" (CW) are used to calculate theQ value variations of gneiss in the preparing rupture process. The obtained results show that the variation state ofQ values by SAR features the shape of relative stability—gradual increment to the maximum—then decrement and final rupture. The variation state ofQ values by EN is just contrary to that by SAR, i. e. with the shape of stability—decrement—increment—and final rupture. The varation state ofQ values by CW is similar to that by EN, its main frequency features the shape of relatively high value—decrement to the minimum—increment—and final rupture. But to the high frequency (higher than the main frequency), the variation state ofQ values features the shape of the stable value-increment to the maximum-decrement-and final rupture. At the same time, the results by coda wave amplitude spectrum show that, when stress reaches 70% of rupture stress, the high frequency component of S wave rapidly reduces (Q _c increasing); at the time of impending the main rupture, the main frequency component reduces with a large scale (Q _c increasing again), this may be the reason which causes the different variation states of two codaQ values. The result of amplitude spectra of P, S (initial wave) waves also show that with the appearance of microcracks the frequency band of S wave turn to be narrow, the high frequency component is reduced quickly, i. e. the S wave spectra have different variation states with different frequency components. That is why theQ _s obtained by different methods have different variation characteristics.     

VARIATION OF THE TRAVEL TIME DIFFERENCE BETWEEN TELESEISMIC PS CONVERTED WAVE AND PP TRANSMITTED WAVE IN FOCAL REGION BEFORE AND AFTER WENCHUAN MS8.0 EARTHQUAKE

When P waves from distant earthquakes meet a velocity discontinuity in the earth's crust and upper mantle, they give rise to a series of converted PS waves besides PP refracted waves. It is possible to monitor the variation of the physical properties of the medium in the limited formation space above the transition zone of the seismogenic zone by measuring the time difference between the teleseismic PS converted wave and the first arrival P wave, that is, time-variation Δt_PS=t_PS-t_P. The advantage of this method is that the transition point of the teleseismic source with similar source is relatively stable at the transition interface, and the accuracy of the measured relative time is high, and the change of the medium in a small range of the seismogenic zone above the conversion interface can be monitored. This paper studies the variation of the travel time difference Δt_PS in focal region before and after Wenchuan M_S8.0 earthquake. We select 2001 to 2012 as the research period, use teleseismic waveforms which occurred in the southern region of Sumatra and Hindu Kush area recorded by Sichuan station YZP and JJS. These teleseisms satisfy 5.0 ≤ M ≤ 6.5, and their waveform signal-to-noise ratio is high with clear initial P-wave motion. The epicentral distance of teleseisms is less than 3 degrees. Then we obtain the variation of the travel time difference Δt_PS between teleseismic PS converted wave and PP transmitted wave recorded during the study period in the two stations. The results show that there is a slow increasing trend of Δt_PS before 2006, and an obvious low value process of Δt_PS appeared in the period about 2 years before the Wenchuan earthquake. The maximum decline was about 0.2~0.3s, more than 4~5 times the measurement error. The low value has a certain degree of return about 2~3 months before the earthquake. The change of arrival time difference indicates that the medium is in different states in different periods of seismogenic process. The sharp decrease in Δt_PS from 2006 to January 2008 may be due to the strong disturbance caused by the stress accumulation of the medium. At this stage, the velocity of P wave and S wave increases with the increase of stress, and the increase of S wave velocity will result in the decrease of Δt_PS. The change of Δt_PS is greatly affected by S wave velocity, so Δt_PS appears to decrease rapidly. Regarding the low value that has a certain degree of return about 2~3 months before the Wenchuan earthquake, the possible reason is that the release of stress is much higher than the accumulation of stress in meta-instability stage. At this stage, the velocity of S wave decreases and the decrease of S wave causes Δt_PS to increase. Then, the Wenchuan earthquake of magnitude 8.0 occurred. It is shown that the teleseismic converted wave method in this paper can monitor the variation of medium's wave velocity before large earthquakes, and it has a good prospect in seismic monitoring and worth further experimental study.     

长白山天池火山区上地壳三维速度层析成像

利用三维地震透射观测系统所获得的莫霍界面反射波对长白山天池火山区下方上地壳的三维速度细结构进行了反演计算，得到了研究区下方地壳内不同深度的P、S波速度层析成像和VP/VS扰动分布图像，结果表明，天池火山区及其周边速度结构所呈现出的非均匀性变化与该区的地质构造和岩浆系统的存在有着密切的关系，异常变化比较明显的低P波速度异常体分布和高VP/VS扰动体意味这些介质体处在“软化高温”的状态，这与该地区广泛分布的地热异常和低Q值结构的研究结果是一致的.     

复杂场地条件下地震动空间效应研究

产生地震动空间效应的主要原因是地质构造的变化会极大地改变复杂场地条件下到达不同观测点的地震波,因此地震波振幅和相位发生明显变化。本文利用有限元和边界元耦合的方法,通过建立雷克子波入射时的非线性响应模型求解复杂场地条件下的波动积分方程,进行复杂近地表构造中地震波场数值模拟,旨在分析复杂场地条件下入射方向不同的地震波的动力响应规律。研究结果显示近地表会极大地改变地震波传播特征,导致地震波传播趋向不一致性的因素是由于复杂场地条件重建了地震波的传播路径,尤其须指出的是非一致性地表位移响应源于软土层中低频子波的传播。在地震响应分析中导致系统不同动力响应的其他原因是主频和地震波速度结构的改变,同时地震波入射角度的改变会引起到达观测点的时间滞后,地表位移响应随着入射角度的增大而增加,地震动的空间效应随着入射角度的减小而愈趋明显。     

2012年印尼8.6级地震应变地震波的Hilbert-Huang时频分析

本文通过端点效应压制的Hilbert-Huang变换, 对大同及沁源台布置的四分量钻孔应变仪记录的印尼8.6级地震激发的应变地震波形进行时频分析, 结果显示印尼8.6级地震的主震和8.2级余震的应变地震波序列各个震相具有不同的时频特征: ① 地震波到达之前的所谓“环境噪声”部分, 瞬时频率低, 瞬时振幅小; ② P波初至时, 高频成分突然增加, 振幅也随即增强; ③ S波到达时, 频率有所降低而振幅剧烈上升; ④ 面波到达时, 振幅进一步剧烈上升达到整个序列的极大值; ⑤ 尾波部分振幅逐渐降低, 但与噪声部分相比频率依然偏高, 振幅依然偏大。 本文也将应变地震波与地震仪记录的地震波进行对比, 虽然应变地震波与地震波波形和Fourier谱具有极高的相关系数, 但从Hilbert-Huang变换得到的边际谱上看, 应变地震波与地震波有显著的区别, 应变地震波比地震波记录的低频成分相对更多。 通过Hilbert谱, 有助于更好地了解非平稳信号的局部特征, 对于突变信号的地震波, Hilbert-Huang变换是一个较好的时频分析工具。     

地震波动力学特征变化指标在短期地震预报中的应用

本文研究了地震前地震波动力学特征的变化，主要是波形和频谱的变化，并从中提取出可望在短期（半年以内）地震预报中应用的一些预报指标。这些指标包括：Ｐ波初始部分波形的时间线性度和空间线性度，初始段平均半周期，Ｓ，Ｐ波频谱的峰值频率，拐角频率，相对带宽度，相对频谱峰值，高频段的谱线斜率绝对值等。同时还研究了短周期地脉动的波形线性度及频谱特征量（峰值频率等）在大震前出现的短期（半年以内）和临震（９天以内）的     

POISSON'S RATIO VARIATIONS OF CRUSTAL MEDIA BEFORE AND AFTER XINYUAN-HEJING MS6.6 EARTHQUAKE IN 2012

We analyzed the variation characteristics of Poisson's ratio in crustal media from January 2009 to December 2012 at 11 fixed seismic stations(for station SCH, it is from January 2006 to December 2012)within an epicenter distance of 200km of the Xinyuan-Hejing M_S6.6 earthquake in Xinjiang on June 30, 2012 using the methods of P wave receiver functions, H-κ stacking of receiver functions, and time sliding window, and obtained the following conclusions: (1)The crustal media's Poisson ratio of five stations in an epicenter distance less than 130km showed a significant and long-lasting decline about 2~3 years before Xinyuan-Hejing M_S6.6 earthquake. Taking the crustal Poisson ratio mean value as reference, the decrease ranges between 0.003 and 0.014, the decrease in 4 stations are more than twice the mean error. The variations of the Poisson's ratio in crust are characterized by "V" shape or "double V" shape. Earthquakes occur at the end of the formation of "V" shape. After the occurrence of earthquakes, the Poisson's ratio continues to rise. The earliest initial fall appeared in July 2009 at WUS station which has the minimum epicentral distance(77km). The Poisson ratio of the crustal media of 6 stations with epicentral distance more than 150km fluctuated up and down around the mean value, and there is no significant decline or persistent low value. (2)We analyzed the arrival-time variations of the quasi-repetitive receiver functions Ps converted wave(t_Ps)of the 3 stations WUS, SCH and XNY and found that the travel times of Ps converted waves became smaller in the crust before the earthquake and increased after the earthquake. (3)Through the comprehensive analysis on the descending process, decline ranges, variations process, duration of Poisson' ratio, the Ps converted waves arrival time variations, the original time of earthquake, and the number of stations, it is inferred that the cause for Poisson's ratio anomalous variations is the change of physical properties of crustal media in the process of earthquake preparation and occurrence. Since the variation characteristics of crustal media may be related to the earthquake magnitude, the size of seismogenic area, the medium properties under stations, and the focal distance, whether the medium variation characteristics exist before and after Xinyuan-Hejing M_S6.6 earthquake will need more earthquake cases analyses. (4)The H-κ stacking of receiver functions is used to calculate the velocity ratio. Because P-wave velocity is given, this method can only be applied when the Ps converted wave velocity of Moho surface of receiver functions changes before an earthquake. With the application of receiver functions to the analysis of more earthquake cases, we can gain more insights into the variation of crustal medium parameters during the seismogenic process. This observation indicates that the receiver function method may become a new approach to detect the Poisson's ratio change of the crustal media before strong earthquake under the condition of high seismic network density.