三维层状黏弹性半空间中球面SH、P和SV波源自由场

采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数.     

基于瞬变电磁矩变换的快速三维反演方法

瞬变电磁法的严格三维反演计算复杂、占用资源多,在普通计算机上难以实现.本文引入瞬变电磁矩变换的概念,提出一种快速三维反演方法.该方法基于阻性限制(resistive limit)特性,建立包含异常体的三维大地的一阶矩响应正演算法,根据不同约束条件,选择优化的最速下降法实现瞬变电磁快速三维反演.文中通过含异常体的三维大地正演一阶矩与仿真数据一阶矩的对比,验证了快速三维正演算法的有效性,之后在不同约束条件下,利用优化的最速下降法实现了对含噪声的仿真瞬变电磁数据的快速三维反演.结果表明,该方法能够在普通计算机上短时间内较为准确地反演出地下异常体的体积和位置,在瞬变电磁数据的实时解释工作中具有良好的应用前景.     

融合多元遥感信息的新右江次级活动断裂解译

提出了融合多光谱影像、全色影像、DEM等多元遥感信息的断裂特征增强方法,建立了以地形和水体的空间特征变化为基础的断裂解译标志,解译并发现桂西北巴马-马山断裂带和右江断裂带之间存在一条北西向展布、长约40km的新的右江次级活动断裂,结合野外地震地质调查进行了验证。多元遥感信息融合处理中,将KL变换、HIS变换、乘积变换、比值变换、小波变换等融合方法进行了全面比较,对融合后影像进行主分量分析,确定了乘积变换融合的影像增强方法和解译底图,显著突出了活动断裂特征。     

近断层应力变化快速计算方法及对强余震空间分布的指示意义:以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反演海潮负荷位移参数是可行的。     

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

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

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

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

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

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

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

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