基于ITRF2005框架的南美板块运动及形变分析

本文通过利用ITRF2005框架速度场以及IGS参考站近两年观测数据解算所得速度场对南美板块的运动及形变特征进行分析。结果表明,南美板块整体以约12.4mm/a的速率向北西方向运动,其欧拉矢量参数与NNR-NULVEL1A模型基本一致。形变特征表现为中南部区域西缘因纳兹卡板块俯冲向东凹陷,东侧受大西洋中脊海底扩张推力向西运动,东西向的挤压使其约以22mm/a的速度收缩并阻碍其北向运动。板块南端则受南极洲板块挤压向北西运动。     

“准拉格朗日法”和“欧拉法”的数学一致性与三次插值函数算法时间积分方案

从预报方程组通式和欧拉算符出发,用泰勒级数展开,推导出时空间微商余项为二阶、四阶完全预报方程组。可以证明,同阶时空间微商余项的准拉格朗日法和向前差分欧拉法时间积分方案,具有相同的物理意义和数学一致性,相比之下,传统中央差欧拉法时间积分方案是"简单格式"与"增大不可预测计算误差"并存。进而讨论用"三次插值函数"实现二阶时空间微商余项准拉格朗日法、或同阶向前差分欧拉法(可"二选一"),它们应该分别替代"双线性插值"准拉格朗日法和传统时空间中央差欧拉法,因前二者时空间微商余项及计算精度高于后二者。所以,三次插值函数算法可将准拉格朗日法和欧拉法时间积分方案、以及CFL判据统一起来。由于三次插值函数具有对原函数"变量场"的数学定律"收敛性"和二阶可导"最优性",且一次"三次插值函数"运算,即具对网格变量场二阶可导拟合"等价性":不仅拟合变量场斜率、还拟合其曲率和挠率。并因周期"三次插值函数",可作全球变量场"三次插值函数"二阶可导拟合,实现全球"三次"数值模式,并且可按变量场曲率判断,作变量场局域或单点平滑,保持"三次"模式时间积分的稳定性。     

铁矿床深部矿体定位预测的一种新途径 ———以西藏恰功铁矿为例

恰功铁矿是近几年新发现的矿山,矿体的形态、大小、位置、产状、边界等几何特征还没有清楚的认识,为了弄清这些问题,文章首先通过小渡多尺度分析方法将恰功铁矿磁异常资料进行了位场信号分离,然后应用欧拉反褶积方法反演了异常与其场源深度的对应定量关系,最后把位场分离后的异常场源深度叠加起来,就可以对矿山深部找矿提供一种全新的思路.     

基于方向tilt-Euler的三维磁数据快速反演

三维磁异常快速反演是磁法勘探进行地质解释的一项重要内容.本文在方向tilt梯度基础上进行导数处理,并进行合理改进,得到了6个改进的方向tilt梯度导数,丰富了不同方向上的磁异常信息.对常规欧拉反褶积方程进行三个方向求导,所构成的方程组与6个改进的方向tilt梯度导数进行结合,从而构建出了两个基于方向tilt梯度导数的快速反演方程组(方向tilt-Euler法),用于反演场源的位置参数.另外,将三个方向导数下的欧拉反褶积进行平方和及均方根处理,得到了场源构造指数的求解公式.单一模型试验表明,6个改进的方向tilt梯度导数异常基本不受磁化方向影响,可从不同方向上反映磁异常信息,方向tilt-Euler能够准确地反演出场源的位置及构造指数.叠加模型试验证实了方向tilt-Euler法相对于tilt-Euler法,具有反演解聚集度更强、连续性更好、反演准确度更高、受噪声影响更小及无虚假反演解等优点.将本文方法应用于内蒙古塔木素航磁资料处理中,获得了丰富的地下磁源分布信息,为了解该地区隐伏断裂构造及岩体分布等提供了有效依据.     

对欧拉矢量附加不等式约束的RELSM模型及其应用

借鉴选权拟合思想,对区域地壳运动的整体旋转和线性应变模型中欧拉矢量部分附加不等式约束,将区域地壳运动的欧拉参数约束到由多年观测资料确定的平均欧拉参数值附近,建立附有不等式约束的整体旋转和线性应变模型。对中国大陆构造环境监测网络在环渤海区域的近几期GPS观测速度场数据进行拟合分析,结果表明,本文方法通过不等式约束可适当增强块体整体运动趋势的平滑性,获得对速度场更高的拟合精度,可进一步增强模型的物理解释能力。     

BHA受压失稳的模型及临界钻压公式

在教学及实践中发现,欧拉压杆模型与底部钻具组合(BHA)受压失稳模型之间存在较大差别,使用欧拉模型计算钻具的临界失稳钻压有很大的局限性。为此,推导了钻具在自重作用下失稳的临界钻压及临界失稳长度公式,以便使公式更符合BHA自身的受力特点,提高计算准确性。      

塔里木地块上石炭统的古地磁研究

本文对塔里木地块西北缘柯坪地区的3条上石炭统康克林组剖面进行了古地磁研究。95个定向岩心样品的逐步热退磁和交变退磁揭示出3个磁性组分,第一个是具有低解阻温度谱或低矫顽力谱在现代地磁场作用下形成的粘滞剩磁A(D=10°;I=56°,α95=4°);第二个是红色石英砂岩的具有反极性的高温组分B(D=217°,I=-43°,λ=56°,Φ=184°E;A25=6°详细研究表明,该分量受附近晚二叠世脉岩侵入的重磁化;第三个是灰岩高矫顽力谱的反极性组分C（D=241°、I=-51°,λ=41°N,Φ=160°E,A（?）=4°）。这一分量通讨倾斜检验,在95%概率水平上明显不同于塔里木地块目前已有的石炭纪以后的古地磁方向,它代表了原生剩磁。研究结果表明:（1）晚石炭世塔里木地块已和哈萨克斯坦地块、苏联西天山地块以及西伯利亚地块发生碰撞:（2）苏联西天山地块相对塔里木地块曾逆时针旋转了70°（欧拉极位置λ=8°N,Φ=69°E）,而塔里木地块则相对西伯利亚地块逆时计旋转了21°（欧拉极位置λ=20°N,Φ=34°E）。从该区已有的晚古生代至早三叠世占地磁数据看,这些地块之间的相对旋转运动发生于早三叠世之后。     

多分辨率地震信号反褶积

基于二进小波变换提出了一种新的反褶积方法─-多分辨率地震信号反褶积．在地震信号二进小波变换域中的各尺度上分别进行其分辨率随小波尺度变化的反褶积,利用不同分辨率反褶积结果之间的相关性,以及测量噪声随尺度的衰减特性,从低分辨率反褶积结果逼近高分辨率反褶积结果．理论分析和实验表明,该方法有较高的精度,并且在较低信噪比情况下有好的效果．     

高分辨率三维地震资料处理技术及应用

高分辨率地震勘探技术贯穿于野外采集、室内资料处理和解释各个环节。文章从三维地震资料处理角度出发,探讨高分辨率、高信噪比和高保真度处理的关键技术措施,即波场净化处理,振幅、频率和相位补偿,叠前反褶积,静校正和动校正及叠后频谱整形。最后以实际资料说明经高分辨率处理后,频宽达到350Hz,给解释人员提供优质剖面,为分辨小间距薄煤层奠定基础。