用双三次样条函数和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.     

一种从经纬度网格资料获得正方形网格资料的插值方案

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

圆柱螺线法测斜计算中的数值方法

圆柱螺线法是井眼轨迹计算中最常使用的方法之一,在实钻轨迹监控、中靶分析和预测、井身质量评价等实际问题中都有重要的应用.但是当井段2个端点的井斜角之差或者方位角之差等于零时,圆柱螺线法坐标计算公式需要进行特殊处理才能正常计算;在井斜角之差或者方位角之差很小时,坐标计算会出现较大误差.对圆柱螺线法坐标计算公式进行了改进,用一组计算公式统一处理井斜角之差或者方位角之差等于零和不等于零的情况,提高了坐标的数值计算的精度和稳定性.对圆柱螺线型井眼轨迹的平均井眼曲率的数值计算进行了研究,提出了使用高斯型数值积分公式计     

确定任意磁性体磁化方向的样条函数法

由磁异常Z_α和重力异常g,利用样条函数的微、积分性质,从泊松方程出发,直接解得磁性体的磁化方向。该方法适用于任意形状磁性体的磁异常。     

用样条函数计算结果研究冲绳海槽热源机制

对冲绳海槽１９８４－１９９０年的实测热流值，利用数字滤波方法计算，得到反映该区深部地壳热状态的区域热流值为８０－１６０ｍＷ／ｍ＾２，用Ｂ样条函数法对其作数值模拟，计算出海槽地温场模式，并以此研究冲绳海槽海底高热流的形成机制，得出冲绳海槽热流值的区域性变化与綦中地质构造格局呈正相关，高热流是该区现代构造活动所致，其明显特征是地壳减薄，地幔上拱等结论。     

任意曲线曲面上重磁位场转换的B样条函数法

提出用B样条函数求解曲线、曲面上重磁位场的向上延拓,水平、垂向导数计算,磁异常分量互换的方法。该方法的特点是:原理简明,程序通用性强,计算精度高。