芦山M_S7.0地震与Northridge M_W 6.7地震地表变形场对比研究

2013年四川芦山"4·20"M_S7.0地震和1994年美国Northridge M_W6.7地震的近场同震位移均表现出破裂未达地表的特征,利用1999~2007年GPS资料对比分析了2次地震的震前变形特征与同震位移的关联性。结果显示,芦山震源区存在较为显著的挤压应变积累(挤压应变率方位角约为284°),沿芦山地震破裂节面剪切变形不明显;芦山地震的同震位移场表现为逆冲兼少量左旋特征,反映了整个孕震周期的应变积累特性,可能表明震前10余年的观测并不能代表整个孕震周期特征。Northridge地震前,地表变形表现为右旋-挤压特征,同震以逆冲破裂为主,近场表现出类似"圆屋顶形"的东西两侧水平运动和垂直向上运动,GPS观测到的震前变形与同震释放的变形存在一定程度的不一致性。相同震中距测点结果表明:Northridge地震的同震地表位移幅度明显大于芦山地震,前者震前挤压应变积累速度显著大于1999~2007年芦山地震震源区的应变积累背景。     

数值模型显示的隐伏弯曲断层同震地表位移特征

余震精定位资料显示,芦山7.0级地震破裂面可能为弯曲程度较高的三维弯曲断层.相关研究显示,这种弯曲断层的位错方式和破裂面同震应力加卸载模式与普通平直断层有明显不同.文中采用无限半空间位错模型模拟显示,隐伏弯曲断层和平直逆断层引起的地表位移特征相似,但是弯曲断层引起的地表水平位移更接近区域整体的地壳缩短方向,缩短方向水平位移的量明显高于同等规模的平直逆断层,因此能更好地传递断层上盘大范围物质的水平运动.相对于平直断层,弯曲断层下盘水平位移随距离衰减十分明显.同等规模的弯曲断层导致的同震地表抬升小于平直逆断层或左旋逆断层引起的同震抬升,但能造成更明显的地表下降.由于地震规模较小,GPS等低密度空间分布的形变观测可能无法有效分辨芦山地震震源结构是否为弯曲断层.对震源结构的细节研究,还有待于利用高空间密度和高分辨率的形变观测资料.     

盲逆断层震后变形耦合机制研究——以2017年伊朗MW7.3地震为例

震后变形是地震周期中一个重要阶段,深入研究其产生机制有助于提高对地震演化过程、断层性质及地震危险性评估的认识.2017年伊朗M_W7.3地震InSAR观测资料丰富,为研究扎格罗斯造山带山前褶皱带盲逆冲断层震后变形机制提供了切入点.合理详细的震源模型及高精度的形变观测是开展震后余滑或黏弹松弛研究的首要前提,本文采用多视角InSAR资料联合远场波形数据反演该地震的破裂滑动分布,并解算得到伪三维地表位移,利用InSAR时序分析提取震源区主震后一年半内变形特征,显示震源区震后形变特征与同震位移场类似,LOS向形变速率最大值约8 cm·a^-1,震后变形预示震源区可能存在明显的震后应力调整现象.利用震后半年内形变资料约束的纯运动学余滑模型表明同时期释放能量为矩震级M_W6.7,进一步探讨应力驱动震后余滑及下地壳黏弹性松弛对震后变形的贡献,基于分层黏弹模型的模拟计算表明震后余滑和黏弹松弛效应的耦合模型可以更好地解释地表震后变形特征,其中震后余滑主要分布在同震破裂区的上方浅部,对主震后一年半内震后地表变形起主导作用,震源区壳下黏滞系数量级下...     

一种子像素偏移追踪算法的实例验证

作为一种中强震形变场获取技术,子像素偏移追踪算法可以与InSAR、GPS等构成重要互补。验证了一种以强度互相关、双线性插值和二维高斯回归三种技术为基础的、新发展的子像素偏移追踪算法,且利用该算法求取了2014年2月新疆于田地震和2013年9月巴基斯坦地震的同震形变位移场。2014年于田地震地表破裂以左旋走滑为主,最大左旋滑移量达到1.5 m。2013年巴基斯坦地震地表破裂主要表现为左旋走滑的弧形破裂,最大滑移量可达到10 m。经计算得到的巴基斯坦地震同震形变位移场与前人的研究结果基本一致,精度可达到百分之三个像素水平。结果表明新发展的子像素偏移追踪算法可以有效求取中强震的同震形变位移场。     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

昆仑山MS81地震震源参数的多破裂段模拟研究

昆仑山MS81地震的已有研究结果在破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数方面存在较大差异.本文采用D-InSAR技术首次获得昆仑山MS81地震干涉同震形变场，结合野外科学考察的实测值，进行了主破裂带InSAR视线向变化量的分解，通过对InSAR分解结果、野外科学考察、遥感解译等多源数据综合分析，重新划分了昆仑山地震的次级破裂段.进而通过对地震南北盘同震应变的分析，发现了昆仑山地震的南北两盘分别受挤压和拉张两种应力作用，研究表明多种岩石在拉张和压力作用下其最小主应力下的杨氏模量表现出非线性弹性特征，从而提出对昆仑山地震地表位移及震源特征参数分析时应考虑非线弹性介质导致的非线性弹性位移分布特征.基于上述原因，本文对Okada线弹性位错模型的算法进行了改进，提出了“多震源、非均一位错分量、多破裂段叠加”的线弹性模型，该模型模拟出的形变场干涉纹图较好地体现了地震形变场的分布特征，并由此获得了一套较为完整的地震发震断层的几何学特征参数，为破裂带长度、破裂面方向、破裂面大小等震源破裂特征参数研究提供了较好的解释.     

基于高分辨率地形数据的富蕴M8.0地震地表破裂带精细特征

同震地表破裂带的空间展布及形变特征是地球深部断层活动在地表的直观地貌表现,不但记录着地震破裂和断层运动的信息,还反映了区域应力和地壳运动状况.因此,开展震后地震地表破裂带调查对于了解发震断层的构造活动尤为重要.高精度地形观测技术可以获取前所未有的高时空分辨率的地球表面特征,为辨别历史地震地表破裂遗迹、提取地表同震位移、活动构造地质填图等提供高质量数据.文中选取富蕴1931年地震地表破裂带作为研究区,利用SfM(Structural from Motion)摄影测量技术生成分辨率为1 m的数字高程模型(DEM),详细识别地表破裂并测量冲沟的右旋位移.基于地表破裂的几何及构造地貌特征,将富蕴地震地表破裂带由北向南分为S1、S2、S3、S44段,其间以挤压隆起或拉分盆地相连接.沿破裂带共获得194组最新冲沟的右旋水平位移,得到1931年同震位移的平均值为(5.06±0.13)m.同震位移局部缺失或突变的区域与几何阶区的位置也有良好的对应关系.以上结果填补了对富蕴地震地表破裂精细形态研究的空白,也进一步展示高分辨率的地形数据在活动构造研究中良好的应用价值.     

亚像素相位相关法在获取汶川地震近场形变中的应用

震后地表实际破裂带的分布及其近场的形变特征,是理解块体运动学特性、断层破裂特征、地震发生机制等科学问题的十分重要的约束条件。基于InSAR获取的汶川地震同震形变场,由于发震断层附近同震形变梯度巨大,沿断层带出现了非相干条带,以致于无法获得断层附近的形变量。而基于亚像素级的光学影像偏移量法为获取断层附近大形变分布提供了可能。文中以SPOT卫星影像为数据源,采用光学影像偏移量法获得了什邡及茂县地区的水平位移形变场。结果显示龙门山断裂带上至少2条断裂同时发生破裂,形成了主要地表破裂带(龙门山镇-高川破裂带)和次级地表破裂带(汉旺破裂带),沿龙门山镇-高川破裂带平均位移量为4～6m,在高川附近伴随的平均右旋水平位移为1～3m; 汉旺破裂带因逆冲导致水平缩短,平均位移量一般为1～2m。汶川-茂县断裂带没有明显的地表破裂带。研究表明,利用光学影像相位相关法能够获得近断层位错量,可以成为InSAR手段的重要补充。     

GNSS观测的2021年5月22日玛多MS7.4地震同震位移及其约束反演的滑动破裂分布

2021年5月22日青海省果洛州玛多县发生M_W7.4地震,此次地震产生的地表破裂在空间上表现出明显的分段特征.本文基于不同来源的GNSS连续观测网数据获取了此次地震的精细三维同震形变场,结果显示:观测到的最大水平位移量达到280 mm,最大垂直形变量仅为25 mm,暗示此次地震的逆冲分量较小;此次地震具有较为明显的左旋走滑特征,同震形变基本对称,在NW-SE向的影响范围更广,该方向上水平同震形变大于3 mm的震中距范围超过500 km.进而,本文以余震精定位结果和GNSS观测到的三维同震形变场为约束,构建了地表破裂线为折线、倾角为85°、倾向西南的断层模型,反演了滑动破裂分布.结果显示:滑动破裂分布在震中两侧不均匀,均破裂到地表,破裂深度达到15 km左右,最大滑移量为4.73 m,计算的矩震级为M_W7.37.该结果与余震精定位结果具有很好的一致性,破裂的极值区正好位于早期余震空区,推测该余震空区未来的发震风险性较低.最后基于反演结果模拟计算了震中区域形变和应变场,结合应变值在断层地表迹线东南侧呈现挤压特征和已有的研究成果,推测此次地震增强了巴颜喀拉块体在东部地区挤压应力的积累特征,导致东部地区发震危险性增强,值得后续跟踪研究.     

2010年青海玉树Ms7.1级地震地表破裂带和形变特征分析

通过对SPOT卫星影像上地表破裂的目视解译,以及对ALOS PALSAR卫星数据进行D-InSAR形变提取和分析,结合地震活动性、震源机制、活动构造等资料,确定了发震断层空间分布、断层性质和同震形变场分布特征.结果显示,玉树地震发生在甘孜-玉树断裂带上,总体走向约为300°,断层近乎直立.根据相干性强弱将Ⅰ区地表破裂划...     

湖泊藻类动态模拟

以过程模拟研究为主线,从藻类生长与吸收,生态过程,尺度与空间异质性等方面阐述了湖泊藻类动物模拟研究的发展态势与存在问题,认为结构动力学模型、生态与水动力耦合模型、面向对象模型和神经网络模型等具有良好的发展前景。     

基于VTK的医学图像系统研究

医学三维数据场可视化是当前科学计算可视化应用的重点。本文论述了可视化工具包VisualizationToolkit(VTK)的主要特点,介绍了三维医学可视化系统实现的原理和方法。根据一系列人体断面图像,应用VTK实现了三维医学图像重构,对重构后的三维图像进行解剖分析;并实现对人体的二维图像逐帧浏览及柱状统计。实验结果表明,VTK具有使用灵活,功能强大的优点,是医学三维数据场可视化的有力工具。     

有关钢结构地震作用的讨论

简要介绍了6本国内标准(TJ11-74、TJ11-78、GB J11-89、GB50011-2001、JG J99-98、CECS160:2004)有关钢结构地震作用的计算方法,包括荷载组合方法和底部剪力法。通过归一化后的结果比较,评价了各本标准地震作用的大小。分析表明,除JG J99-98和GB50011-2001 2本标准由于考虑钢结构的阻尼比为0.02或0.035后,其地震作用有显著增大外,其它4本标准的地震作用处于接近的水平。考虑到钢筋混凝土和砖石结构等的地震作用水平在新旧抗震规范中并无多大变化,过分提高抗震性能较优的钢结构的地震作用水平显然是不合理的。     

基于Lg波识别判据的特性分析

Lg波是检测、识别地下核爆炸的重要震相.通常情况下,地下核爆炸的低频段Lg波能量相对较强,基于Lg波的识别判据在高频段表现较好.然而在低频范围内却表现出难以用球对称爆炸源模式解释的现象,如Lg/Pg判据在1 Hz附近失效等现象.本文利用理论地震图方法,分析了伴随层裂过程产生的重要辅助源——CLVD源对Lg波低频成分的调...     

断层错动对隧道的影响分析

隧道通过潜在地震地区时.存在活动断层的错动形成的重大工程地质问题的威胁.基于断裂儿何学和震源的运动学特征,以广州地区活动断裂为研究背景,通过在断层左右盘施加一定的速度荷载,模拟了断层错动达到0.28 m时对隧道的影响.得出隧道结构在断层错动影响下的位移和应力的变化规律.分别建立隧道结构距离断层50 m、100 m、150 m三个模型,比较了三种情况下的位移和应力状况.归纳出隧道与断层之间的安全距离.     

Occurrence of water ponding on soil surfaces depending on infiltration rates on Mongolian rangeland

The occurrence of water ponding on soil surfaces during and after heavy rainfall produces surface run‐off or surface water accumulation in low‐lying areas, which might reduce the water supply to soils and result in a reduction of the soil water that plants can use, especially in arid climates. On Mongolian rangeland, we observed ponded water on the surface of a specific soil condition subjected to a heavy rainfall of 30 mm/hr. By contrast, ponded water was not observed for the same type of soil where livestock grazing had been removed for 6–8 years via a fence or for nearby soil containing less clay. We measured the infiltration rate (the saturated hydraulic conductivity of the surface soil, K_s) of the three sites by applying ponded water on the soil surface (an intake rate test). The results showed that K_s in the rangeland was lower than the rainfall intensity in the site where water ponded on the soil surface; however, K_s of the soil inside of the fence has recovered to 3 times that of the soil outside of the fence to exceed the rainfall intensity. Heavy rainfall that exceeds the infiltration rate occurs several times a year at the livestock grazing site where we observed ponded water. Slight water repellency of the soil reduces rain infiltration to increase the possibility of surface ponding for the soil.     

基于地震动参数的建筑物震害研究

本文介绍了多步简化的push—over(MSP)方法,引导出求解单体建筑物易损性指数和群体建筑物易损性指数的办法,并通过建筑物平均易损性指数曲线,来评价城市的总体抗震性能。本文不采用过去以烈度为地震输入的易损性评定方法,采用地震动参数作为城市建筑物易损性的评定标准。给出的易损性指数曲线,可用来查找不同地震动参数下各类建筑物的易损性结果,比较城市各类建筑物的抗震性能。本文研究可望在城市建筑物震害评估中有好的应用前景。     

基于ActiveMQ平台的地震消息服务探讨

重点介绍了ActiveMQ消息平台的特点和功能,并对遥测台网中速报、台站参数共享、标定信息等部分应用进行了探讨和阐述,并给出了部分相关代码和应用实例流程图。     

New results on radiation effects on human health

Humans are exposed to ionizing radiation all the time, and it is known that it can induce a variety of harmful biological effects. Consequently, it is necessary to quantitatively assess the level of exposure to this radiation as the basis for estimating risks due to ionization radiation. During the Work Package 2400 of the COST-724/WG-2 action, a number of spacecraft and aircraft experiments have been performed with both active and passive detectors. A large data base has been created. In this contribution we would like to stress the results obtained and their importance in three particular directions: (i) Simultaneous investigation of galactic cosmic rays on aircraft and on the International Space Station (ISS); (ii) Onboard spacecraft neutron contribution as estimated on the basis of the comparison of results measured with MDU/Liulin equipment onboard ISS, foton capsule and a commercial aircraft flying at subsonic altitudes; (iii) Complex analysis of the results of long term measurements onboard a Czech Airlines aircraft. The results obtained are presented, analyzed, and discussed, and their complementary nature is underlined. The contribution represents a version of the Final Report of the Work Package 2400 of the COST-724/WG-2: Radiation Environment of the Earth.     

A note on edge waves on a flat shelf

Abstract

A variational approximation to the dispersion relation for trapped waves on a flat shelf of depth h ₁, bounded internally by a vertical coast and externally by a semi-infinite ocean of depth h ₂>h ₁, is obtained through an integral-equation formulation that accounts for all of the non-propagated modes that are excited at the discontinuity in depth (the conventional formulation of the edge-wave problem allows only for the propagated mode on the shelf and the dominant, non-propagated mode in the deep water). Coriolis effects are neglected. The exact result in the limit ω² h ₂/g↓0 (ω = angular frequency) is obtained by conformal mapping and compared with the variational approximation, which proves to be quite accurate over the entire range 1>h ₂/h ₁>x. The effects of the higher-order, non-propagated modes are found to be small for the long waves observed over the Southern California shelf by Snodgrass, Munk and Miller (1962).