位移观测值协方差矩阵模拟及其在整体解算地壳应变场中的应用

借助高斯型函数模拟分散分布的不同监测网点位移观测值的协方差矩阵,并将其应用于广东新丰江地区地壳应变场的计算中．结果表明：本文的方法不仅计算简便,模拟效果及计算稳定性较好,而且对联系分散分布的单元有一定的意义     

地壳应变速率与地震活动关系的研究——从地壳应变场探索强震活动的场源关系

通过研究应变积累速度（亦即地壳应变速率）与地震活动的关系,籍以探讨从地壳应变场中寻找近期强震危险区的方法。针对不同地区的地震活动水平特别是地震复发期的差异,丁国瑜（１９８４）指出,各地震区大震复发时间的差异主要取决于由大震孕育环境和构造物理条件所决定的构造活动速率（亦即应变累积速率）的差异。从近几年地震实况的初步（阶段性）检验看,应用大地测量中垂直形变速率梯度判定强震危险区,具有较好的效果和可喜的前景。作者从理论分析、实际资料对比及计算机模型研究等方面来讨论地壳应变速率与地震活动水平之间的关系,进而为场源关系研究中的“场中求源”的科学思路探索可能的技术途径,并讨论了从应变速率场中寻找强震危险区的方法。     

气侯平均风场作用下热带太平洋主流系和赤道行星波研究

将两层约化重力原始流体动力方程耦合气候月平均风场,数值计算流场基本能够正确反映热带太平洋上层主流系和温跃层的空间分布和季节变化．在气候平均条件下,东太平洋125°W附近经向风应力可激发出高阶混合Rossby重力波．海洋高阶赤道Kelvin波流速模态可从西太平洋边界传播到东太平洋边界,而高阶赤道Kelvin波温跃层模态从西太平洋边界东传后,在中太平洋受到高阶混合Rossby重力波诱发的西传温跃层扰动的阻挡．     

基于中小地震应变能密度的地震活动性图像分析

当前的地震活动图像分析方法未能充分利用大量的中小地震观测资料。本文提出中小地震能量密度值的定义及计算方法,以华北地区为研究区域,利用1970～2006年ML2.0～4.0地震资料,计算了中小地震能量密度图像,初步分析了中小地震活动的时空迁移特性,以及中小地震能量密度图像变化与强震的关系。该方法为地震危险区圈定、地震活动性参数确定、历史强震参数复核等提供了一种可能的工具。     

大同窗地震活动频次和应变能特征

通过对大同地震活动窗自开启以来(1989年10月~2016年10月)的地震频次和应变能进行统计分析,发现大同窗的月频次衰减拟合结果呈现分段交替活动特征,可分为3个活跃时段和2个背景活动时段。当处在活跃时段时,山西和张渤这2个地震活动带的地震密集发生,但在正常背景活动时,2个地震带发生的地震则较少。应变能统计结果表明,鄂尔多斯块体周缘强震对大同窗的地震活动有较直接的影响,显示出强震的震级越大、与大同窗的距离越近,则大同窗出现应变能异常就越早、幅度也越大;另外,大同窗受青藏块体强震的影响较强,而受西太平洋板块俯冲的影响较弱;大同窗应变能在周边块体强震后的反应,可能与其所处的地质构造环境有关,因此对应力场的变化较敏感。     

基于分量钻孔应变资料分析测区构造应变变

根据“双衬套”理论模型和应变理论, 以易县台分量钻孔应变观测资料为例解算测区构造应变变化。 由于易县台分量钻孔应变元件灵敏度不一致, 首先进行相对实地标定, 数据信度由0.7557升高到0.9344; 其次, 应用理论固体潮对易县台分量钻孔应变进行绝对标定, 得到耦合系数; 最后, 通过应变换算得到易县台观测点附近的应变参数, 面应变ε_a每年以2440 ns的速度变化, 最大剪应变ε_s以每年230 ns的速度变化, 应变主方向为138°。 综合结果表明, 易县台测区附近的构造应变呈稳定拉伸的状态。     

Research note: the effect of strain amplitude produced by ultrasonic waves on its velocity

The effect of the amplitude of ultrasonic waves propagating through a sample is not often taken into account in laboratory experiments. However, ultrasonic waves can produce relatively large strain inside the sample, and thus change the properties of the sample. To investigate the effect of strain amplitude on the P-wave velocity, a series of ultrasonic wave propagation experiments were carried out on three different media. All measurements were performed at 1 MHz central frequency and at the strain levels inside propagating waves of  ∼3.0 × 10⁻⁶ to 6.0 × 10⁻⁵ without applying confining pressure to the sample. Strains in the waves were measured by a laser Doppler interferometer upon wave arrival on a free surface of the sample. The ultrasonic velocities were measured by a pair of P-wave transducers located at the same measuring point as the laser beam of the LDI. The effect of strain on P-wave velocity varied for different material. The P-wave velocity was calculated using both a first arrival and a first maximum peak at different applied voltage. The P-wave velocity remained unchanged for a pure elastic medium (aluminium); however, the velocity increased continuously with the increasing of the strain for polymethylmethacrylate and Gosford sandstone. For Gosford sandstone, velocity increases up to 0.8% with strain increase from 7.0 × 10⁻⁶ to 2.0 × 10⁻⁵. This effect of velocity increase with the strain induced by an ultrasonic wave can be explained by the in-elasticity of both the polymethylmethacrylate and Gosford sandstone samples.     

利用GNSS连续观测资料分析川滇地块与断层活动的方法与尝试

本文利用GNSS连续站的资料分析2016、2017年川滇块体和巴彦喀拉块体速度场、应变场特征,认为道孚、西昌和马尔康附近应力积累相对集中。分析GNSS连续站点的相对速度发现,川滇块体内GNSS连续站的运动状态呈明显的顺时针旋转,巴彦喀拉块体也有顺时针旋转态势。此外,利用块体相对运动求解鲜水河断裂运动速率,鲜水河断裂西部的运动速率比东部稍大,约有2—3mm/a的差异。相比2016年,2017年鲜水河断裂运动速率变小,有利于缓解鲜水河断裂应变能的积累。     

用应变地震观测求解震源矩张量的基本原理

震源机制解,即对地震矩张量的推断,对于地震研究具有至关重要的意义.应变地震观测是张量观测,与摆式地震仪的位移矢量观测不同,可以为地震研究提供新的数据源.本文讨论用应变地震观测求解震源矩张量的基本原理.在距离震源足够远的地方,地震波可以看成平面波,其性质决定于震源矩张量.假设平面地震波的应变张量可以由震源矩张量通过坐标变换计算得到,就可以通过观测应变地震波求解震源机制.这个假设至少对于双力偶震源机制是成立的.由此可以证明,在理想的无限介质中,只要有两个以上不同地点的应变地震波观测,就可以解出震源矩张量.这为解决震源机制问题提供了新的方法.目前的地震矩张量求解方法需要两方面的条件:或者需要很多观测点(例如体波反演),或者需要长周期地震波资料(例如面波反演).这些方法只适用于分析比较大的地震.对于小震,因为通常其震中周围不会有足够多的摆式地震仪观测点观测到其地震波,而地震波周期又短,难以利用传统方法给出可靠的震源机制解,所以只需少数观测点就能求解震源矩张量的新方法特别有意义.用应变地震观测求解震源机制,可以给出更为精确的结果.     

GEOLOGY OF THE NORTHERN ARUN TECTONIC WINDOW

GEOLOGY OF THE NORTHERN ARUN TECTONIC WINDOW1　BordetP .Recherchesg啨ｏｌｏｇｉｑｕｅｓｄａｎｓｌ’ＨｉｍａｌａｙａｄｕＮ啨ｐａｌ,r啨ｇｉｏｎｄｕＭａｋａｌｕ[R].EditionsduCNRS ,Paris ,196 12 75 . 2　BordetP .G啨ｏｌｏｇｉｅｄｅｌａｄａｌｌｅｄｕＴｉｂｅｔ (Himalayacentral) [J].M啨ｍｏｉｒｅｓｈｏｒｓｓ啨ｒｉｅｄｅｌａＳｏｃｉｅｔ啨ｇ啨ｏｌｏｇｉｑｕｅｄｅＦｒａｎｃｅ,1977,8:2 35～ 2 5 0 . 3　BurcfielBC ,ChenZ ,HodgesKV ,etal.TheSouthTibetanDetachmentSystem ,Hima…