近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.     

中国大陆西部重力场变化与强震危险性关系

本文基于中国大陆西部流动重力测量获得的区域重力场变化数据, 尝试对其进行多种位场变换, 在一定空间尺度约束下定义了6种区域重力场变化指标量. 在此基础上, 以中国大陆2002年、 2005年和2008年等3期数据为研究对象, 计算并统计了中国大陆西部M_S≥6.0地震震中位置处震前重力场变化及各种导出参数值, 尝试初步构建以地震预报为目的的地震重力学科指标体系参数. 研究结果表明, 震前区域重力场变化与地震发生位置没有显著的统计关系, 但是重力场变化的垂直梯度和解析信号模量具有一定的地震预报意义, 适合作为地震前兆异常使用.      

利用PPP反演海潮负荷位移参数

选取高纬地区的3个CORS站以及中纬地区的SHAO测站,利用PPP浮点解计算测站未经海潮负荷位移改正的坐标时间序列,采用傅里叶变换方法提取分潮的周期信息,进而反演其中8个主要分潮的海潮负荷位移参数,将其与5个全球海潮模型DTU10、EOT11a、GOT4.8、TPXO.7.2、FES2004计算得到的海洋负荷位移参数比较。结果表明,PPP反演得到的8个主要分潮的海潮负荷位移参数与海潮模型计算结果具有较好的一致性,两者均方根差异为mm级,说明利用PPP反演海潮负荷位移参数是可行的。     

Locata单点定位精度研究

介绍Locata伪距和载波相位单点定位的方法。通过仿真数据和试验,分析基准站个数对定位精度的影响,比较伪距定位和观测值组合定位的精度,验证了载波相位定位结果的稳定性。分析表明,Locata静态单点定位的精度可以达到cm～mm级。     

基于HHT变换瞬时能量谱的长周期结构地震反应分析

基于HHT变换得到瞬时能量谱,并以简谐地震反应为理论基础进行多自由度体系在地震动作用下的瞬态反应计算。通过对一个长周期(十五层)结构工程实例的地震瞬态反应计算分析,揭示长周期结构有别于短周期结构动力响应的独特特性:长周期结构位移响应的瞬时能量峰值的时间滞后于地震动能量峰值的时间,且主要受到瞬态振动的控制,并以低阶振型的控制为主,第一阶振型峰值位移滞后的时间最长,滞后的时间主要取决于结构的第一阶自振周期。研究结果表明,长周期结构在强烈地震动作用下,地震动强度出现的最大时刻并不总是对应于结构破环与倒塌的时间。     

3D Hough Transform在激光点云特征提取中的应用

近年来,激光扫描测量技术得到了迅速的发展和广泛的应用,激光点云数据处理成为当前计算机科学和测绘科学等领域的研究热点。特征识别在图像数据处理中占有重要地位,在三维激光点云数据处理中,特征可以用于多站数据拼接、大数据量的简化压缩及建模,因此对特征的提取同样具有重要的研究价值。Hough变换最初用于二维图像的直线提取,后来被广泛应用于二维图像的特征识别,但其运算量极大,并且直接运用到三维可能会导致错误识别等问题。本文根据三维点云特征提取的需要,通过对二维和三维标准Hough变换的介绍和分析,结合激光点云数据处理特点,提出了一种依据法向计算结果进行Hough变换的新方法,该方法不仅能识别平面特征,还可将类似墙角等特征突变点进行有效分离,以便进行轮廓提取等操作。对模拟数据和激光扫描仪实测数据分别进行处理,平面提取的有效率和正确率达到要求,从而证明该方法能够满足点云中特征提取的需要。     

Matlab平台下遥感影像的北方岸线提取研究——以大连长兴岛为例

在Matlab平台上, 利用卫星遥感影像数据, 结合RS与GIS技术, 选取大连长兴岛作为研究区域, 使用阈值分割法、边缘检测法及小波变换法进行了海岸线的提取, 并对结果进行对比分析。试验结果表明, 边缘检测法中Canny算子与小波变换法对海岸线的提取结果与真实地物最为接近, 其中Canny算子可以较为准确地提取围填海区域的海岸线。     

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。     

关于频率电磁测深几个问题的探讨(八)——频率电磁测深与瞬变电磁测深的关系

从探测原理出发,说明人工源频率电磁测深与瞬变电磁测深的关系密切,这是一个问题的两个方面,可以互相转换。两种探测方法实质是一样的,具有共性;但毕竟从频率域和时间域两个侧面反映问题,又具有各自特性。两种方法的原理所依据的麦克斯韦方程,其正演公式可通过傅立叶变换相互转换;从实际探测所发收的信号特征,场区划分与探测深度等主要方面阐明两种探测方法异同,以加深对探测理论与方法的认识,结合实际择优应用。      

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.