岩石蠕变过程中声发射活动的实验研究

在常温下、单轴加压下对两种岩石完整样品进行蠕变破坏过程中声发射活动实验.根据实验获得的数据讨论了蠕变三个阶段中声发射活动特征.在瞬态蠕变阶段,随着加载,声发射活动逐渐增强;在稳态蠕变阶段,在较低的声发射活动背景上,可能出现较强的信号,声发射幅值-频度关系的 m 值可能稳定于某一低值,也可能下降;进入加速蠕变阶段,声发射活动水平增强,m 值下降加快,也可能表现为明显回升后再下降,并与岩石样品的性质有关.将蠕变破坏三个阶段中声发射活动与地震孕育过程中某些地震活动图象变化(如 b 值、前震活动)进行了联系,提出地壳岩石蠕变破坏机制形成前震——主震——(余震)型地震序列类型的可能性.      

前震或广义前震识别的频谱偏移法及其应用研究

在前期研究工作的基础上,本文提出了一种前震或广义前震识别的频谱偏移法。其判定指标为,正常背景的中等地震事件频谱较宽,而具有前震或广义前震性质的中等地震具有频带变窄且向低频端偏移集中的特征。究其根源可能与强震前震源区或主震相关区域亚稳态扩展破裂或区域介质破碎损伤、裂隙扩展成核等过程有关。对2010年4月14日青海玉树7.1级、2014年2月12日新疆于田7.3级等地震前震的计算结果显示,2次前震的频谱都发生了偏移,与非前震性同等震级地震频谱差异明显。     

The 2010 M w 7.2 El Mayor-Cucapah Earthquake Sequence, Baja California, Mexico and Southernmost California, USA: Active Seismotectonics along the Mexican Pacific Margin

The El Mayor-Cucapah earthquake sequence started with a few foreshocks in March 2010, and a second sequence of 15 foreshocks of M?>?2 (up to M4.4) that occurred during the 24?h preceding the mainshock. The foreshocks occurred along a north?Csouth trend near the mainshock epicenter. The M _w 7.2 mainshock on April 4 exhibited complex faulting, possibly starting with a ~M6 normal faulting event, followed ~15?s later by the main event, which included simultaneous normal and right-lateral strike-slip faulting. The aftershock zone extends for 120?km from the south end of the Elsinore fault zone north of the US?CMexico border almost to the northern tip of the Gulf of California. The waveform-relocated aftershocks form two abutting clusters, each about 50?km long, as well as a 10?km north?Csouth aftershock zone just north of the epicenter of the mainshock. Even though the Baja California data are included, the magnitude of completeness and the hypocentral errors increase gradually with distance south of the international border. The spatial distribution of large aftershocks is asymmetric with five M5+ aftershocks located to the south of the mainshock, and only one M5.7 aftershock, but numerous smaller aftershocks to the north. Further, the northwest aftershock cluster exhibits complex faulting on both northwest and northeast planes. Thus, the aftershocks also express a complex pattern of stress release along strike. The overall rate of decay of the aftershocks is similar to the rate of decay of a generic California aftershock sequence. In addition, some triggered seismicity was recorded along the Elsinore and San Jacinto faults to the north, but significant northward migration of aftershocks has not occurred. The synthesis of the El Mayor-Cucapah sequence reveals transtensional regional tectonics, including the westward growth of the Mexicali Valley and the transfer of Pacific?CNorth America plate motion from the Gulf of California in the south into the southernmost San Andreas fault system to the north. We propose that the location of the 2010 El Mayor-Cucapah, as well as the 1992 Landers and 1999 Hector Mine earthquakes, may have been controlled by the bends in the plate boundary.     

雁列式断层组合变形过程中的声发射活动特征

着重研究拉张或挤压型雁列式断层标本变形过程中声发射(简称AE，下同)的时空演化特征. 结果表明，预置构造对AE空间分布格局具有较强的控制作用，随着差应力的增加，AE首先在两个端点附近丛集，之后向两端点连线附近集中，出现明显的破裂局部化现象，较大事件还通常集中于某一端点附近反复发生. 前期微破裂丛集图象指示后期宏观破裂的扩展方向及扩展范围. 拉张和挤压型雁列区宏观破裂方向分别与轴向应力方向垂直和平行. 雁列式断层标本变形过程中，破裂前的弱化阶段相对较长、弱化过程明显. 微破裂事件累积频次指数增长可能是系统失稳前的典型征兆之一，而雁列区宏观破裂之后，AE数量逐渐减少、应变释放相对减弱. 摩擦滑动过程中，大的粘滑失稳前未见有AE活动前兆性的增强过程.雁列区的b值变化在失稳前显示趋势性降低——快速回升这一典型的变化过程，b值降低一般发生在差应力增强过程之中，并有可能延续至弱化阶段，而快速回升则一般发生在破裂失稳前的短时期内. 对比研究表明，构造差异所导致的b值差异远大于b值随差应力的增加而产生的变化量，而对同一构造标本，力学状态的改变会导致AE序列时序特征的急剧变化，较高的加载速率对应较高的应变能释放及明显的低b值. 先期破碎带由于较低的破坏强度，其对差应力的微小变化具有特殊敏感性，从而成为源、兆分离、窗口或敏感点效应等地震活动性前兆现象的一种可能的原因.      

Spatio-temporal characteristics of acoustic emission during the deformation of rock samples with compressional and extensional en-echelon faults

Introduction By the experimental techniques of acoustic emission (AE) and surface strain observation, a lot of achievements on the study of deformation and fracture process of samples with en-echelon faults has been obtained (Bomblakis, 1973; LIU, et al, 1986; DU, et al, 1989; LI, et al, 1989; ZHANG, FANG, 1989; ZHANG, LI, 1989; MA, et al, 1995a, b; LI, et al, 2000). However, there still exist many problems. For example, most of discussion on the spatio-temporal evolution of microf…     

云南地区强震的短临中强前震序列特征初探

对1965￣2003年云南地区强震前3个月内发生过的短临中强前震进行整理,将其按距离主震发生时间分为直接前震、临震前震与短期前震三类,然后通过地震序列动态分析地震活动,发现各类前震有以下显著特征:(1)直接前震:3、4级地震时间间隔临震时骤减;将地震序列从时间上分成三等分后,第二时段的地震频次没有明显衰减。(2)临震前震:震后第一天地震的频次较高而第二天相对过低,同时同一台站P波初动符号基本上一致或振幅比稳定在某一范围。(3)短期前震:震后10天前震频度衰减至零;20-30天左右地震频度骤增,同时h、p值也反应出类似的特征,30天h值迅速减小接近1或小于1,而b值基本稳定且均低于区域背景值。     

1975年2月4日海城地震(M=7.3)的前震系列/

我们研究了1975年2月4日海城地震(M=7.3)前震的定位及其辐射图象.用六个区域台的到时,相对于一次参考地震,将前震及主震进行了定位.这组前震开始相互很接近,然后随时间及其分布有一定的延伸.最大前震前,前震均位于直径约两公里的小体积内,而在最大前震后,其活动则向北西和南东方向扩展,形成六公里长的在北西方向上伸展的分布带.初动及 P 波 S 波振幅比表明,在前震系列中有两种不同的断裂机制.我们推测这两种辐射图象可能和前震处于分布带的不同部位有关.可能主震震源不处于前震震源所决定的断层上,而是位于这组前震南面6公里、且较这组前震浅几公里处.我们认为,在垂直于主震破裂面及前震分布带走向的方向上,前震和主震相距这么大的距离,可能是由于主震时产生滑动的断层是以雁行排列的.分析了在前震期间断层上滑动所引起的应力变化,认为由前震引起的主震断层上剪应力的增加是很小的.因之由前震直接触发主震的可能性不大.      

Use of Potential Foreshocks to Estimate the Short-term Probability of Large Earthquakes, Tohoku, Japan

This study seeks to construct a hazard function for earthquake probabilities based on potential foreshocks. Earthquakes of magnitude 6.5 and larger that occurred between 1976 and 2000 in an offshore area of the Tohoku region of northeast Japan were selected as events for estimating probabilities. Later occurrences of multiple events and aftershocks were omitted from targets. As a result, a total of 14 earthquakes were employed in the assessment of models. The study volume spans 300 km (East-West) × 660 km (North-South) × 60 km in depth. The probability of a target earthquake occurring at a certain point in time-space depends on the number of small earthquakes that occurred per unit volume in that vicinity. In this study, we assume that the hazard function increases geometrically with the number of potential foreshocks within a constrained space-time window. The parameters for defining potential foreshocks are magnitude, spatial extent and lead time to the point of assessment. The time parameter is studied in ranges of 1 to 5 days (1-day steps), and spatial parameters in 20 to 100 km (20-km steps). The model parameters of the hazard function are determined by the maximum likelihood method. The most effective hazard function examined was the following case: When an earthquake of magnitude 4.5 to 6.5 occurs, the hazard for a large event is increased significantly for one day within a 20 km radius surrounding the earthquake. If two or more such earthquakes are observed, the model expects a 20,000 times greater probability of an earthquake of magnitude 6.5 or greater than in the absence of such events.     

Simulation of Unstable Fault Slip in Granite Using a Bonded-particle Model

—?A bonded-particle model is used to simulate shear-type microseismic events induced by tunnel excavation in granite. The model represents a volume of granite by an assembly of 50,000 individual particles bonded together at points of contact. A plane of weakness is included in the model and this plane is subjected to increasing shear load while the normal load across the plane is held constant. As shear stress in the model increases, bonds begin to break and small acoustic emissions (AE) result. After enough bonds have broken, macro-slip occurs across the large portions of the fault in an unstable manner. Since the model is run dynamically, seismic source information can be calculated for the simulated AE and macro-slip events. This information is compared with actual results obtained from seismic monitoring around an underground excavation. Although the modelled events exhibit larger magnitudes than the actual recorded events, there are many similarities between the model and the actual results, namely the presence of foreshocks before the macro-slip events and the patterns of energy release during loading. In particular, the model provides the ability to examine the complexity of the slip events in detail.     

Effect of Injected Water on Hydraulic Fracturing Deduced from Acoustic Emission Monitoring

—In order to investigate the effects of injected water in hydraulic fracturing, experiments were conducted on cubic granite specimens, comparing fracturings induced by conventional water injection with those induced by pressurization of a urethane sleeve, thereby realizing "hydraulic fracturing" without the use of fracturing fluid. In both experiments, a shear type mechanism was found to be dominant in fault plane solutions of AE events. However, in the case of water injection, cracks extended rapidly with large drops in hole water pressure and bursts of AE, whereas in pressurization by the urethane sleeve, cracks extended stepwise with no such large drops in hole pressure and no bursts of AE. The difference in crack extension in the two experiments can be analyzed by comparing relations between crack length and stress intensity factor of mode I at a crack tip. The observation and analysis indicate that existence of fracturing fluid like water helps initiated cracks to extend rapidly and widely in hydraulic fracturing in actual HDR fields. Received September 12, 1996, accepted January 24, 1997     

海城地震前震系列与震群——兼论前震系列在地震预报中的一种功能

对比了海城前震系列与其周围同期震群的异同.两者的震源空间分布都比较集中,一般只有几公里的尺度.发震机制都比较稳定.但有些震群在其主要地震发生后,机制变化较显著,借此有可能识别一部份震群.但也有些震群其机制的稳定性始终不比海城前震系列差,因之借此无法识别前震系列.频次的分布两者也有非常相似的震例.总的印象是有些震群在晚期有可能识别,但未发现有确认前震(系列)的肯定判据.这一初步的结论仅仅是由于我们分辨技术的局限,抑是具有成因上的深刻含义是很值得思考的.最后引入了震级序列的缺震特征,借此可以判断系列的高潮是否已经渡过.这对(从其他证据认识到的)具有大震危险的地区,利用前震系列来改善大震临震时刻的机率性判断是有用的.注意到频次分布中有相隔12小时的低值点,对比就地的固体潮具有较好的同步性,说明固体潮可能对前震系列以及震群具有触发作用.      

基于深度学习到时拾取自动构建长宁地震前震目录

将深度学习到时拾取、震相关联技术与传统定位方法联系起来,构建一套连续波形自动化处理与地震目录自动构建流程,对于高效充分利用地震资料,提升微震检测能力具有十分重要的意义.我们应用最新发展的迁移学习震相识别技术、震相自动关联技术,对长宁M S6.0地震震中附近21个台站震前半个月(6月1日—6月17日)的连续记录波形进行P、S震相识别、震相自动关联和初步定位,并应用传统绝对定位和相对定位技术得到了长宁地震震前微震活动的绝对和相对定位目录.其中绝对定位目录能在较小的误差范围匹配85%的人工处理目录,其发震时刻平均误差为0.36±0.07 s,震级平均误差为0.15±0.024级,水平定位平均误差为1.45±0.028 km,其识别的1.0级以下微震数目是人工的8倍以上,将长宁地震震前微震目录的检测下限提升至M L-1左右,证明了基于深度学习到时识取和REAL(Rapid Earthquake Association and Location,快速震相关联和定位技术)震相自动关联来构建微震目录具有较好的实用性.我们的自动地震目录揭示了长宁M S6.0主震所发生的区域震前异常频繁的微震活动,以及与区域内盐矿注水井的关联性,更好地描绘了这些微震活动的时空演化特征,其空间活动性分布特征与长宁M S6.0余震序列的分布一致.     

The M<Subscript>W</Subscript> 4.5 Vallorcine (French Alps) earthquake of 8 September 2005 and its complex aftershock sequence

On 8 September 2005 a moderate M_W 4.5 earthquake occurred in the north-western Alps midway between Chamonix (France) and Martigny (Switzerland). The focal mechanism corresponds to a right-lateral strike-slip on a N60°E fault plane. The foreshock–mainshock–aftershock sequence is investigated on the basis of data recorded by a temporary network of 28 stations deployed for 1 month just after the mainshock, and data from permanent, regional seismic networks. Absolute and relative locations of more than 400 events are obtained with a mean uncertainty of approximately 0.2 km. Small foreshocks, the mainshock, and early and late aftershocks are located relative to the main aftershock set. The seismic sequence exhibits a surprisingly complex structure, with at least five clusters on distinct fault planes. The main elongated cluster agrees with the location of the mainshock, its hypocenter being 4.3 km below sea level. We discuss the relationship between the right-lateral fault beneath the Loriaz peak (the source of the Vallorcine event), the nearby normal Remuaz fault, and the regional seismotectonic stress field.     

华北地区中强以上地震与前兆震群的相关分析

应用“八五”攻关推广的有关震群参数计算的软件 ,对华北地区 (33°～ 43°N,1 1 4°～ 1 2 4°E)自 1 970年以来所发生的 40次震群事件的各项参数 (b,U,H,K,F,ρ等值 )进行了计算 ,旨在依据前兆震群的判定指标 ,寻找前兆震群与中强以上地震 (MS≥ 434)的对应关系。结果表明 ,无论是单项指标或综合指标的报准率均在 65%左右 ,仅有个别达到 75%以上 ,并且有约30 %具有前兆意义的震群属于漏报和虚报。所以 ,如何判定前兆震群和充分利用前兆震群预报地震和监视未来地震活动有待更深入的研究和探讨。     

Diversity of the 2014 Iquique’s foreshocks and aftershocks: clues about the complex rupture process of a Mw 8.1 earthquake

We study the foreshocks and aftershocks of the 1 April 2014 Iquique earthquake of Mw 8.1. Most of these events were recorded by a large digital seismic network that included the Northern Chile permanent network and up to 26 temporary broadband digital stations. We relocated and computed moment tensors for 151 events of magnitude Mw?≥?4.5. Most of the foreshocks and aftershocks of the Iquique earthquake are distributed to the southwest of the rupture zone. These events are located in a band of about 50 km from the trench, an area where few earthquakes occur elsewhere in Chile. Another important group of aftershocks is located above the plate interface, similar to those observed during the foreshock sequence. The depths of these events were constrained by regional moment tensor (RMT) solutions obtained using the records of the dense broad band network. The majority of the foreshocks and aftershocks were associated to the interplate contact, with dip and strike angles in good agreement with the characteristics of horst and graben structures (>2000 m offset) typical of the oceanic Nazca Plate at the trench and in the outer rise region. We propose that the spatial distribution of foreshocks and aftershocks, and its seismological characteristics were strongly controlled by the rheological and tectonics conditions of the extreme erosive margin of Northern Chile.     

Using the match-and-locate method to characterize foreshocks of the July 2019 MW6.4 Ridgecrest,California earthquake

The July 2019 M_W6.4 Ridgecrest, California earthquake and its distinct foreshocks were well recorded by local and regional stations, providing a great opportunity to characterize its foreshocks and investigate the nucleation mechanisms of the mainshock. In this study, we utilized the match-and-locate (M&L) method to build a high-precision foreshock catalog for this M_W6.4 earthquake. Compared with the sequential location methods (matched-filter + cross-correlation-based hypoDD), our new catalog contains more events with higher location accuracy. The M_W6.4 mainshock was preceded by 40 foreshocks within ～2 h (on July 4, 2019 from 15:35:29 to 17:32:52, UTC). Their spatiotemporal distribution revealed a complex seismogenic structure consisting of multiple fault strands, which were connected as a throughgoing fault by later foreshocks and eventually accommodated the 2019 M_W6.4 mainshock. To better understand the nucleation mechanism, we determined the rupture dimension of the largest M_L4.0 foreshock by calculating its initial rupture and centroid points using the M&L method. By estimating Coulomb stress change we suggested that the majority of foreshocks following the M_L4.0 event and M_W6.4 mainshock occurred within regions of increasing Coulomb stress, indicating that they were triggered by stress transfer. The nucleation process before the M_L4.0 event remains unclear due to the insufficient sampling rate of waveforms and small magnitude of events. Thus, our study demonstrates that the M&L method has superior detection and location ability, showing potential for studies that require high-precision location (e.g., earthquake nucleation).     

A summary of seismic activities in and around China in 2021

In this article, we review the general characteristics of seismicity in and around China and the overall statistics of earthquake damage in 2021, focusing on several significant events and related scientific topics. Among them, the largest event is the M_S 7.4 Madoi earthquake in Qinghai Province, northwest China. The event marks another M_S ?≥ ?7 earthquake occurring near the boundary of the Bayan Har Block that has ended a remarkable quiescence of the M_S ?≥ ?7 earthquakes within the Chinese mainland. In addition, the M_S 6.4 Yangbi earthquake in Yunnan Province, southwest China draws the most attention because of its abundant foreshocks, which are well recorded by the densely distributed seismic stations in the surrounding regions. Regarding this event, we review several recent publications focusing on the Gutenberg-Richter b-value change and the physical mechanism of foreshocks associated with this sequence. The M_S 6.0 Luxian earthquake in Sichuan Province, southwest China has caused serious damage with a relatively low magnitude, partly because the focal depth of the mainshock is relatively shallow (3.5 ?km). It is another strong earthquake occurring within the southeast Sichuan basin with low historical seismicity yet has increased significantly since 2015, probably due to shale gas development and associated hydraulic fracturing.     

Effect of Sliding Rate on the Activity of Acoustic Emission During Stable Sliding

--The mechanical and statistical characteristics of acoustic emission (AE) events during stable sliding are investigated through a laboratory experiment using a granite specimen with a pre-cut fault. Numerous AE events are found to be generated on the pre-cut fault, indicating that microscopically unstable fracture occurs during macroscopically stable sliding. The composite focal mechanism solution of AE events is determined from the first motion directions of P-waves. The determined mechanism is consistent with the double-couple one expected for the slip on the pre-cut fault. The source radii of large AE events are estimated to be about 10 mm from the widths of the first P-wave pulses. These indicate that the AE events are generated by shear fracture whose faulting area is a part of the pre-cut fault plane. The occurrence of AE events as a stochastic process approximately obeys the Poisson process, if the effect of mutually dependent events consti tuting clusters is corrected. The observed amplitude-frequency relation of AE events approximately follows a power law for a limited amplitude range. As the macroscopic sliding rate increases, the number of AE events per unit sliding distance decreases. This rate dependence of the AE activity is qualitatively consistent with the observation that the real area of contact between sliding surfaces decreases with an increase in the sliding rate as reported in the literature.     

由地震分布丛集性给出断层参数的一种新方法

由于大范围内地质构造的复杂性和介质的非均匀性,发震断层面的几何形态一般十分复杂.如果大地震的破裂过程涉及多个断层的活动,则发震断层并非是单一断层平面,而是多个断层面的组合.利用地震空间位置分布丛集性,即震源点成丛位于断层面附近的假设,结合稳健扩充算法和主成分分析给出一种可以重构活断层网络三维空间结构的新方法.该方法每次从震源点集中处开始,利用假设检验扩充子断层面,并得到多个子断层面.接着按震源点属于最近断层面的准则把各子断层面内的震源点进行竞争,并根据一定假设合并和删除一些子断层面,最后用主成分分析确定每个子断层面参数.于是可根据地震事件目录给出一组矩形区域来描述断层面网络结构,其中每个矩形断层面由其位置、走向和倾角确定.通过计算机模拟发现,新方法可成功地重建模拟地震目录的断层面,最后用于南加州1992年6月28日发生的Landers地震部分余震目录中,得到各个子断层面参数与已知地质破裂或隐伏断层相当一致.     

前震研究进展综述

简要介绍了当前国内外前震研究领域的一些进展。对比了不同前震定义条件下前震序列的共性特征,对几种主要的前震机理及主要的前震识别方法进行了简要综述,对其特点及存在的问题进行了评述和讨论。前震指主震之前在主震断层面上、紧邻主震破裂起始点发生的小地震。主震发生之前的一系列前震活动构成前震序列。在不同的前震定义条件下,具有“直接前震”的震例比例从10%至40%不等。理论上前震可用级联应力触发或预滑动模型进行解释。前震空间上主要集中分布在距离主震10~75km范围内,但其时间分布形式非常复杂,大多发生在主震前1~2天,部分前震序列的地震活动率显示明显的加速特征,但许多前震序列单独来看却往往显示主-余型序列的衰减特征。序列地震震源机制一致、序列b值偏低是前震序列最突出的特征。前震的发生与主震破裂形式和构造环境似有一定关系,在有限的前震震例中,逆冲型主震似乎具有相对更多的前震。部分震例的研究结果显示,随主震的临近,前震震源深度有逐渐下迁的特点。到目前为止,震前很难判定一次地震或一个地震序列是否前震或前震序列,所使用的前震识别主要有基于统计类比的方法、基于震源机制一致及衍生的相关方法以及基于对地震成核过程精细检测的方法。从现有不多的震例研究结果来看,尽管地震时空丛集和震源机制高度一致是前震序列的最显著特征,但却不是判定前震序列的充分条件。由于成核的破裂扩展速度和滑动位移有随时间较快增大的趋势,因而基于对地震成核过程精细检测的方法有望在前震识别中发挥更为重要的作用,但需更多震例进行验证。