青海祁连克克里5.1级地震的主应变比分析

利用祁连县地应变观测资料，对１９９２年６月青海祁连克克里５．１级地震信息进行了综合分析，提出了主应变比异常理论，在该理论的基础上探讨了该次地震的发震机制，同时分析了祁连建站以来的地应变观测资料，认为肃南～祁连断裂段上及附近发生的地震具有相似性，是在变化的附加引张应力场作用下，载荷弱化，促使正断层或正走滑断层易于活动所导致的。     

祁连山构造带红外测距观测揭示的应变强度时空演变特征

利用截止2007年最新的、跨祁连山构造带红外测距观测资料，结合区域GPS观测所得形变背景，获取跨断层水平应变强度时空演变特征，分析其与2002年玉门Ms5．9级、2003年民乐Ms6．1．级等地震孕育一发生过程的关系。结果表明：①震前数月至1年左右，震中所在断层段（或整个构造带）应变加速或相对其它断层段显著，压性为主；其内一定数量场地最大剪应变强度明显“增强”或面应变强度“压性增强”（至绝对值2．0左右及以上），震后减弱或呈现张性转折。②2004年下半年以来祁连山构造带总体恢复压性变化为主的应力场状态，反映一定程度应力应变积累特性；2007年构造带西段压性应变明显增强。     

基于GNSS的青藏高原东北缘地壳运动场及强震趋势研究

利用“中国大陆构造环境监测网络”GNSS数据研究1998—2018年青藏高原东北缘排除同震影响等干扰后的速度场、主应变率场、最大剪切应变率场、面应变场等的变化,活动断裂滑动速率变化、跨活动断裂基线变化等。将研究区域内的二级块体再分区,获得各次级块体内部的应变率变化;获取研究区域地壳运动场的趋势性、动态特征。研究结果显示,阿尔金断裂带中东段、祁连块体和柴达木块体交界、巴颜喀拉块体与羌塘块体交界、祁连块体南边界中段、海原—六盘山断裂带和西秦岭北缘断裂带西段的逆冲运动,祁连块体北边界西段、庄浪河断裂的左旋走滑运动,祁连块体北边界东段、西秦岭北缘断裂带东段的左旋逆走滑运动,都属于造成一定程度地壳变形的持续性局部应变增强活动。阿尔金断裂带东段、东昆仑断裂带中西段、祁连块体北边界、庄浪河断裂北段、海原断裂南段、六盘山断裂北段、西秦岭北缘断裂带东段可能存在闭锁,未来十年可能发生M_S6.0以上地震。     

由GPS观测结果推导中国大陆现今水平应变场

以中国大陆及周边近400个GPS测站的水平运动速率为基础，给出了现今地壳水平应变场结果表明:①中国大陆水平应变为西强东弱，剪应变数值大于正应变数值(绝对值)，应变量级一般为10-8/a，局部区域达到10-7/a，但应变分布不均匀；②南北向应变最突出的部位为中国西南部西段的喜马拉雅条带、西北部的36N~42N段及柴达木断块的北缘；③东西向应变西边缘变化最大．此外，由西向东还具有正负交替的变化特征；④REN(东-北向剪切应变)与Rmax(最大剪切应变)数值较大的区域分别是喜马拉雅条带、西北部的36N～42N段、柴达木断块的西部、川滇菱形块体，以及阿拉善、祁连及塔里木断块的交界区；⑤青藏块体周边以面收缩为主，内部则以面膨胀为主．其以北的地区以面收缩为主．西界数值最大，东部数值最小(除燕山构造带外)；⑥西部西区主压应变为南北向，主张应变为东西向．西部东缘区主压应变为近东西向，主张应变为近南北向．川滇菱形块体主应变的方向发生了很大的变化，北部地区为东西压南北张，南部地区则恰好相反；⑦中国大陆的应变模式可能是断块模式与连续模式的组合．此外，小尺度优势应变可能是剪切应变．造成上述结果除与印度板块的碰撞及边界耦合有关外，还与深部物质的活动及地壳介质的物性有密切的关系．必须指出，由于GPS测站在空间上分布的不均匀性，那么，由此而来的应变场，其应变尺度也不一样．      

祁连山中段现代应力活动与浅层应力预报地震研究

本文在研究全新世断层活动的基础上，运用构造解析、水系应力场等多种方法，对祁连山中段的现代应力活动状态进行了综合研究，其结论是：区域主压应力方向为北东－南西向，平均Ｎ５０°Ｅ；但在大区范围内显示了应力场的空间变化，局部地区随时间亦有动态变化。此外，还结合祁连站的观测资料，研究了中强震和应变曲线的对应关系，并取得了较好的预报效果。     

基于GPS速度场采用RELSM模型分析青藏块体东北缘的形变—应变特征

利用青藏块体东北缘地区1999～2001年GPS观测获得的地壳水平运动速率场,通过对该地区进行块体划分,将该地区划分为9个块体,应用块体的整体旋转线性应变模型(RELSM)估计了各个块体的旋转与应变参数,以及计算了该地区内143个GPS站点的应变参数,以此分析了该地区的应变场的基本特征,结果表明:①阿拉善块体s较稳定,其旋转角为0.630×10-8,运动速率为0.688 mm/a,②相比其他块体,共和块体旋转角最大达到了6.589×10-8 ,运动速率达到了7.296 mm/a,③应变高值区主要集中在祁连山断裂,海原断裂等,在这些地区最大剪应变率达到了7.5×10-8、面膨胀率达到了-2.5×10-8、主压应变达到了-6×10-8.     

循环加载高压流变实验中塑性应变分析

在三轴应力下的流变实验中，通过循环加载方法测出轴向应变，轴向塑性应变和弹性模量随差应力的变化，再根据多孔介质有效弹模理论，计算出岩石变形中孔隙度的变化，在假定岩石骨架不可压缩下，得到体积膨胀应变，然后通过总塑性应变与扩容机制和位错机制应变的关系，分解出位错应变，结果表明：在低围压下（＜100MPa），扩容应变占优势，并伴随着部分位错应变，在高围压下（＞100MPa），位错应变占优势，只有较小的扩容应变。     

基于四分量钻孔应变的天山中段地区构造应变变化

根据钻孔应变观测理论，利用巴仑台、库米什和小泉沟分量钻孔应变观测数据定量计算测区附近构造应变变化。结果显示，三个台站的最大-最小主应变、面应变和剪应变的应变速率相对恒定，主方向大体不变；巴仑台附近区域受张-压应力相互作用，主压应变方向为N22°W,库米什台附近区域受拉张应力作用，主张应变方向为N8°E;小泉沟台附近区域受压应力作用，主压应变方向约为N46°W;精河M_S6.6地震前巴仑台和小泉沟的应变变化速率明显高于平均水平，均呈现在压缩背景下的应变加速变化异常，可以为应变资料同类异常的识别和判定提供参考。     

青藏块体东北缘地壳水平应变场数值模拟研究

基于青藏块体东北缘1999～2001年GPS结果,分别采用块体整体旋转与线性应变模型和弹性力学有限元法这两类地壳形变数值模拟方法,分析了该区地壳水平应变场特征.结果表明:(1)两类方法在研究地壳形变时各具一定的优势,前者对块体整体运动变形及块体与块体之间的相关性研究具有一定的优势,而后者则较强体现出了应变高值区与深大断层在空间分布的紧密结合性;(2)两类方法所获得的应变高值区具有良好的空间分布一致性,主要集中在阿尔金断裂、祁连山断裂中东段、东昆仑断裂和海原断裂这些深大断裂处及其附近,面压缩值达到了-3×10-8以上,最大剪应变值达到了10×10-8以上;(3)应变高值区的空间分布与中强地震的发生具有一定的对应关系.     

新疆钻孔应变异常与强震活动

在新疆钻孔应变18年的观测期间，发生了6级以上强震17组(次)，其中7级以上地震5次，8级以上地震1次。钻孔应变资料在64％强震前记录到了异常变化。其中属于天山和塔里木地块周边的15例强震中，在儿例地震前钻孔应变资料出现了异常，占全部强震的73％。利用钻孔应变资料正确预报的地震数仅4例，占全部强震的23％，占天山及塔里木地块周边地震数的26％。在这些强震中距离台站250km范围以内发生的地震有7例，正确预报的有4例，占57％。     

RESEARCH ON ACTIVITY OF ZHANGJIAKOU-BOHAI FAULT ZONE BASED ON GPS OBSERVATIONS

The northwestern section of the Zhangjiakou-Bohai fault zone starts in the west of Zhangjiakou, extending southeast through Huailai, Shunyi and Tianjin and entering into the Bohai Sea, with a width up to several tens of kilometers, narrow in the west and wide in the east. The Neogene-Quaternary has extended in the northwest and southeast direction, forming a large regional active structure. There are many earthquakes of magnitude 7 or above in the history on the Zhangjiakou-Bohai fault zone and it is also a strong earthquake activity zone in eastern China. Therefore, the modern tectonic activities of this fault zone have an important impact on regional seismic hazard, and are of great significance for earthquake prediction and disaster reduction. In this paper, using the mobile GPS station observation data of 1999, 2007, 2009, 2011, 2013 and 2015, and with the rigid-linear elastic block motion model equation proposed by LI Yan-xing, the horizontal deformation rate and strain rate of the Zhangjiakou-Bohai fault zone of the five adjacent periods of 1999-2007, 2007-2009, 2009-2011, 2011-2013 and 2013-2015 were calculated, the tectonic activity characteristics and evolution of the fault zone were studied. The results show that in the five periods, the average deformation rate of the Zhangjiakou-Bohai fault zone is 1. 74mm/a, the left-lateral strike-slip rate is 1.59mm/a, and the compression rate is -0.59mm/a. The Zhangjiakou-Bohai fault zone is characterized by left-lateral strike-slip and compression on the whole, and the left-lateral strike-slip rate is greater than the compression rate at each period. The strike-slip rate is significantly greater than the compression rate, indicating that the activity of Zhangjiakou-Bohai fault zone is dominated by left-lateral strike-slip faulting with compression. The minimum principal strain rate of the Zhangjiakou-Bohai fault zone in the five periods varies from -12.06×10^-9/a to -4.62×10^-9/a, and the average minimum principal strain axis direction is N63.9°E, with little change in direction. The maximum principal strain rate varies from 1.55×10^-9/a to 5.99×10^-9/a, and the average maximum principal strain axis direction is N333.9°W, the direction does not change much. The strike of the Zhangjiakou-Bohai fault zone is NWW(the overall strike is calculated by N300°W), and the normal strain rate of the fault zone is -5.87×10^-9/a(being compressional), and the shear strain rate is 12.70×10^-9/a. The shear strain rate on the fault zone is about twice the value of the normal strain rate, and the shear strain rate of the fault zone is greater than the normal strain rate, which indicates the shear stress of the 5 periods of 1999-2007, 2007-2009, 2009-2011, 2011-2013 and 2013-2015 is relatively significant, suggesting that the fault plain is dominated by left-lateral shear stress. This suggests that the Japan 3·11 earthquake has little effect on the deformation strain of the Zhangjiakou-Bohai fault zone, and it does not change the nature of activity of the fault zone. The tectonic activity is still inheriting. Since the tectonic activity of the Zhangjiakou-Bohai fault zone has gradually decreased after the Japan 3·11 earthquake, the deformation strain evolution trend has gradually returned to a unified consistent state. Therefore, the deformation strain state of the Zhangjiakou-Bohai fault zone does not have the condition for strong earthquakes.     

张家口地区应变参数变化及其特征研究

根据1988—2009年张家口伸缩仪3方向的线应变资料计算出最大主应变、最大剪应变以及主应变方向等平面应变参数,分析张家口地区的断层特征、张家口伸缩仪观测环境的变化,结合1988—2009年张家口地区Ms≥4.0级以上地震,研究张家口地区应变参数的变化特征。     

恭城—栗木断裂带的应变特征与变形环境探讨

通过系统的有限应变测量与差异应力估算,探讨了恭城—栗木断裂带的应变特征与变形环境,以其推动该断裂带的研究向更深层次发展。结果表明：①付林指数K的变化范围为0.32～1.67;应变主轴平均方位为276°∠21°,与断裂带的整体逆冲方向基本一致;断裂带的平均应变程度为2.26～3.00,属中等应变程度;断裂带的应变类型属压扁型应变与单剪型平面应变的一种复合叠加类型。②断裂带的差异应力值为40～80 MPa,变形深度小于10 km,属半脆性—半韧性的变形域。③断裂带的应变特征与变形环境及宏观构造样式基本吻合。     

地壳应变主体单元模型及其与中强地震关系研究

讨论了目前利用GPS观测结果研究块体及其边界带应变状态的研究进展,阐明了与地震相关的块体边界活动构造附近是地壳应变主体区域的认识. 利用Scholz断层能量动态平衡过程,建立了断层周边应变累积与断层面构造力的一般物理关系,据此提出以断裂带为主体的地壳应变主体单元模型. 并利用GPS基准站连续观测资料讨论了模型力学状态演化与2001年11月14日昆仑山口西8.1级地震过程的关系. 结果表明,强烈压剪背景中出现的相对松弛变化,可能有利于触发型中等强度地震的发生.      

Detection of Aseismic Slip on an Inland Fault by Crustal Movement and Groundwater Observations: A Case Study on the Yamasaki Fault, Japan

To understand the detailed process of fault activity, aseismic slip may play a crucial role. Aseismic slip of inland faults in Japan is not well known, except for that related to the Atotsugawa fault. To know whether aseismic slip does not occur, or is merely not detected, is an important question. The National Institute of Advanced Industrial Science and Technology constructed an observation site near Yasutomi fault, a part of the Yamasaki fault system, and has collected data on the crustal strain field, groundwater pressures, and crustal movement using GPS. In a departure from the long-term trend, a transient change of the crustal strain field lasting a few months was recorded. It indicated the possibility of an aseismic slip event. Furthermore, analyses of data from the extensometers at Yasutomi and Osawa observation vaults of Kyoto University, as well as GPS data from the Geographical Survey Institute (GEONET), revealed unsteady crustal strain changes. All data could be explained by local, left-lateral, aseismic slip of the order of 1 mm in the shallow part of the Yasutomi fault.     

利用GPS数据分析四川“Y”型构造区地壳运动状态

四川地区地质构造复杂,地壳活动剧烈,为了深入揭示该区“Y”型构造区地壳1999年以来近20年的动态演化规律,基于1999—2017年7期GPS数据,解算各周期网格速度场、应变率场,研究地壳应变场演化过程。结果表明：①2008年以前的3期GPS速度场相对稳定,汶川地震后,速度场变化最大的龙门山断裂带由4.0—5.0 mm/a增至8.0—10.0 mm/a;②汶川震后,“Y”型构造区最大剪应变高值区出现在汶川以东,由2.0×10^-8/a增到22.0×10^-8/a;龙门山断裂带以SE或SEE向主压应变为主,变化速率约5.0×10^-8/a—12.0×10^-8/a,鲜水河断裂由震前NS向主拉应变转为震后EW向主压应变,安宁河断裂东侧由震前SE向主压应变6.0×10^-8/a减至震后的2.0×10^-8/a;面膨胀结果显示,由震前低密度梯度带瞬间变为平行于龙门山断裂带走向的高密度变化区,且存在以金川至都江堰、北川至青川为条带的2个正负交替过渡区;③汶川地震发生压应力释放后,该区SEE向压性特征又逐渐增强,且持续至2017年,释放了龙门山断裂带地壳内部SEE向压应力多年累积能量,但汶川地震对鲜水河断裂与安宁河断裂的整体运动状态则无明显触发作用。     

福建省钻孔体应变同震响应分析

通过福建省4个钻孔体应变仪对印尼强震的同震响应波对比分析,研究了福建省体应变变化特征,得出体应变幅度与震级成正比、响应延迟时间与震中距成正比的结论,进一步探讨了福建省各台站对印尼远强震的体应变变化特征及其与断层性质的关系。     

基于GPS应变、地震应变率与震源应力场对帕米尔高原现今构造变形特征的分析

研究帕米尔高原的构造变形特征对于理解印度板块向北推挤过程中的应变分配方式以及应力转换模式具有重要的意义.本文利用区域GPS应变场、地震应变场与震源应力场分析帕米尔高原的构造形变特征.主要结论为:(1)该区域变形主要以NNW-SSE或近N-S向的挤压为主,在高原内部伴有明显的近ENE-WSW或E-W向拉张,应力方向在帕米尔高原与塔吉克盆地区域呈现逆时针旋转的趋势,而在塔里木盆地则显示几乎与帕米尔高原的一致的应力状态,这可能与两侧盆地块体的强度差异有关.(2)安德森断层参数A∅显示帕米尔高原北缘与西侧区域为逆断层应力状态,在高原内部为正断层应力状态,这与GPS应变的结果显示的应变主要集中在主帕米尔断裂与阿莱谷地附近而在高原内部应变较低是一致的,另外应力在喀喇昆仑断裂北段的方向基本平行于断层走向,以及断层北端较低的滑动速率,这说明了地壳挤压缩短可能是帕米尔高原主要的的构造变形特征,并不支持由于边界走滑断裂导致的应变分异或者块体挤出的模式.(3)综合考虑地震应变方向与SHmax从帕米尔北部NNW-SSE方向到天山北部的近N-S方向的转换,GPS应变方向在帕米尔高原两侧盆地都存在不同程度的旋转,应力场安德森参数A∅显示的应力状态以及SKS的结果显示的近ENE-WSW方向,我们认为印度板块向北推挤与天山造山带碰撞导致帕米尔高原不对称的径向逆冲是帕米尔高原现今构造变形的主要成因与构造模式.     

NEAR-FAULT DISPLACEMENT AND DEFORMATION OBTAINED FROM ONE-KILOMETER-LONG FAULT-CROSSING BASELINE MEASUREMENTS-A PRELIMINARY EXPERIMENT AT 2 SITES ON THE EASTERN BOUNDARY OF THE SICHUAN-YUNNAN BLOCK

The current and conventional fault-crossing short baseline measurement has a relatively high precision, but its measurement arrays usually fail to or cannot completely span major active fault zones due to the short length of the baselines, which are only tens to 100 meters. GNSS measurement has relatively low resolution on near-fault deformation and hence is not suitable for monitoring those faults with low motion and deformation rates, due to sparse stations and relatively low accuracy of the GNSS observation. We recently built up two experimental sites on the eastern boundary of the active Sichuan-Yunnan block, one crossing the Daqing section of the Zemuhe Fault and the other crossing the Longshu section of the Zhaotong Fault, aiming to test the measurement of near-fault motion and deformation by using fault-crossing arrays of one-kilometer-long baselines. In this paper, from a three-year-long data set we firstly introduce the selection of the sites and the methods of the measurement. We then calculate and analyze the near-field displacement and strain of the two sites by using three hypothetical models, the rigid body, elastic and composed models, proposed by previous researchers. In the rigid body model, we assume that an observed fault is located between two rigid blocks and the observed variances in baseline lengths result from the relative motion of the blocks. In the elastic model, we assume that a fault deforms uniformly within the fault zone over which a baseline array spans, and in the array baselines in different directions may play roles as strainmeters whose observations allow us to calculate three components of near-fault horizontal strain. In the composed model, we assume that both displacement and strain are accumulated within the fault zone that a baseline array spans, and both contribute to the observed variances in baseline lengths. Our results show that, from the rigid body model, variations in horizontal fault-parallel displacement component of the Zemuhe Fault at the Daqing site fluctuate within 3mm without obvious tendencies. The displacement variation in the fault-normal component keeps dropping in 2015 and 2016 with a cumulative decrease of 6mm, reflecting transverse horizontal compression, and it turns to rise slightly(suggesting extension)in 2017. From the elastic model, the variation in horizontal fault-normal strain component of the fault at Daqing shows mainly compression, with an annual variation close to 10^-5, and variations in the other two strain components are at the order of 10^-6. For the Longshu Fault, the rigid-body displacement of the fault varies totally within a few millimeters, but shows a dextral strike-slip tendency that is consistent with the fault motion known from geological investigation, and the observed dextral-slip rate is about 0.7mm/a on average. The fault-parallel strain component of the Longshu Fault is compressional within 2×10^-6, and the fault-normal strain component is mainly extensional. Restricted by the assumption of rigid-body model, we have to ignore homolateral deformation on either side of an observed fault and attribute such deformation to the fault displacement, resulting in an upper limit estimate of the fault displacement. The elastic model emphasizes more the deformation on an observed fault zone and may give us information about localizations of near-fault strain. The results of the two sites from the composed model suggest that it needs caution when using this model due to that big uncertainty would be introduced in solving relevant equations. Level surveying has also been carried out at the meantime at the two sites. The leveling series of the Daqing site fluctuates within 4mm and shows no tendency, meaning little vertical component of fault motion has been observed at this site; while, from the rigid-body model, the fault-normal motion shows transverse-horizontal compression of up to 6mm, indicating that the motion of the Zemuhe Fault at Daqing is dominantly horizontal. The leveling series of the Longshu site shows a variation with amplitude comparable with that observed from the baseline series here, suggesting a minor component of thrust faulting; while the baseline series of the same site do not present tendencies of fault-normal displacement. Since the steep-dip faults at the two sites are dominantly strike-slip in geological time scale, we ignore probable vertical movement temporarily. In addition, lengths of homolateral baselines on either side of the faults change somewhat over time, and this makes us consider the existence of minor faults on either side of the main faults. These probable minor faults may not reach to the surface and have not been identified through geological mapping; they might result in the observed variances in lengths of homolateral baselines, fortunately such variations are small relative to those in fault-crossing baselines. In summary, the fault-crossing measurement using arrays with one-kilometer-long baselines provides us information about near-fault movement and strain, and has a slightly higher resolution relative to current GNSS observation at similar time and space scales, and therefore this geodetic technology will be used until GNSS networks with dense near-fault stations are available in the future.     

华北地区断层形变异常与地震活动

根据断层形变求解的华北地区异常参数和应变累积率, 研究了华北主要断裂带的断层形变异常及其应震特征。 同一断裂带上的形变异常与该带上的地震有较好的对应性。 河套—张家口—蓬莱活动构造带上发生的强震, 北京地区的断层形变异常参数在总体上几乎都有明显的前兆性异常。 山西带北部的断层形变异常参数对山西断陷带及其延伸部位上的强震同样有较好的反映。 断层应变累积率反映了应力的积累程度, 其值相对较大时, 测点所在地区的地震较活跃, 反之亦然。