基于PSR-WSVM模型的边坡位移预测

为建立高精度的边坡位移预测模型，采用相空间重构（PSR）将边坡位移时间序列数据转换为多维数据，同时构造小波核函数改进的支持向量机模型，建立PSR-WSVM模型并应用于边坡位移预测。将PSR-WSVM模型预测结果与传统支持向量机（SVM）模型、小波支持向量机（WSVM）模型和基于相空间重构的支持向量机（PSR-SVM）模型预测结果进行对比，通过平均绝对误差（MAE）、平均绝对误差百分比（MAPE）和均方根误差（RMSE）3个精度评价指标验证PSR-WSVM模型的可行性。工程实例结果表明，PSR-WSVM模型预测结果的3个精度评价指标都优于另外3种模型，边坡位移预测的精度明显提升。     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

由地面重力数据确定扰动位径向二阶梯度

提出利用地面重力异常数据计算地面扰动位径向二阶梯度，将该梯度的积分表达式转换为卷形式的谱表达式，便于应用ＦＦＴ／ＦＨＴ技术进行快速计算。     

昆仑山8.1级地震前中国大陆的构造应变背景

利用“网络工程”1998～2001年累积的1181个测站的GPS重复观测资料，采用双三次样条函数模型建立中国大陆水平运动模型速度场，用大地坐标在椭球面上计算各类应变场，详细分析了2001年昆仑山8．1级地震前中国大陆水平构造应变场空间分布特征。各类构造应变场的最高值都出现在喜马拉雅构造带与昆仑山地块内(地震断裂带南侧)，鲜水河—安宁河断裂带次之。分析表明，昆仑山8．1级地震正好发生在张性面膨胀应变率的高值区，第一、第二和最大剪应变率高值区边缘的突变区和最大、最小主应变率的高值区。     

关于Diophantine方程的解

对于正整数a ,设S(a)是Smarandache函数。证明了 :方程S(1·2 ) +S(2·3) +… +S(x(x +1) ) =S(x(x +1) (x +2 ) /3)仅有正整数解x =1。     

桐柏─大别山北麓金矿带地质物化探找矿模型

从桐柏─大别山成矿带区域地质、地球物理、地球化学特征的研究入手，分析与金成矿有关的地层、构造、岩浆岩等控矿因素的地球物理场及地球化学特征；建立区域金矿找矿模型。     

台湾海峡西岸重力场变化梯度场图像分析

讨论了应用曲面拟合技术和利用离散的重力观测点的重力变化对特定研究区的重力场变化进行数值模拟以获取该区重力变化空间梯度场图像的具体方法;并应用该方法获得了2004年3月福建安溪3.6级地震前的重力变化空间梯度场图像,发现空间梯度场图像能清晰地反映出地震前重力场群体性变化的差异性;利用此方法对2005年6月~2006年3月间台湾海峡西岸重力场变化进行数值模拟,获得其空间梯度场图像。     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

三峡重庆库区深部地球物理特征与断裂构造

为了深入研究三峡重庆库区岩石圈动力学特征及其对断裂构造活动的控制和影响．并为之提供基础资料，为三峡重庆库区地震、地质灾害的监测与防治提供基础依据，在已有地球物理资料的基础上，从综合地球物理研究角度出发，通过实测地震测线资料的再解释．采用新的处理技术方法，对本区的东西向主剖面和南北向支测线的地震测深资料进行二维射线追踪处理、Pg波成像；选用场分离技术、位移数字成像技术重新处理了重力和航磁资料，通过联合反演来建立深部二维构造剖面，对剖面所揭示的基底构造特征和地壳结构特征、主要断裂构造特征以及莫霍面的起伏特征进行了精细分析和细致研究．从地球物理平面场特征出发建立了岩石圈构造三维框架。研究结果表明，这样的研究思路准确、方法得当．结论可靠；沿该剖面，把可解译的断裂分为Ⅲ级：Ⅰ级为超岩石圈断裂；Ⅱ级为壳断裂；Ⅲ级为盖层断裂。依据地球物理特征，准确揭示了库区的断裂构造特征，达到了预期效果。     

黄土土性的空间自相关性研究

运用土层的随机场模型研究了绛帐、潼关和延安黄土主要物理力学指标沿深度的统计规律,计算出黄土土性的自相关距离,探讨了黄土土性的空间变异性,结果认为:黄土的物理力学指标沿深度均是统计非均匀的,在实际工程的设计中,如果把土层看成统计均匀,其土性参数的标准差会比实际情况大,从而会导致较大的误差。