基于选权拟合法的电离层电子密度层析重构

附加约束的电离层层析算法是解决电离层电子密度反演中不适定问题的主要方法,为避免此类方法中约束权阵的选取不当对电子密度分布重构产生的不良影响,本文将选权拟合法应用到电离层层析成像技术中.该方法特别设计了依据电子密度空间分布特性构造参数权矩阵的方案.新方法有明确的物理意义,挖掘了隐含的信息量,为解决电离层电子密度反演中由于观测数据的不足等因素引起的不适定问题提供了一种新途径,可以得到符合客观实际的结果.数值模拟实验和实测数据的反演结果证实了该算法有效性、可靠性和优越性.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

交通荷载作用下公路路基工作区深度研究

根据路面不平整实测资料,获得了汽车动荷载的计算模型和参数,并利用此荷载模式来分析公路结构的动力响应。考虑到公路结构的三维分层性状,利用Fourier变换技术和求解层状结构的精确刚度矩阵法,研究了公路结构在移动汽车荷载作用下的动力响应问题。为了方便工程应用,利用Odemark厚度和模量的当量转换公式,将公路结构简化为由路面、路基及地基组成的3层体系,并以此为基础分析了路基动应力的衰减规律,同时重点研究了汽车轴载大小、当量路面厚度及路基模量对路基工作区深度的影响,最后建议了路基工作区深度的定量表达式,研究结果可为公路路基设计提供参考。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

介质密度反演偏导矩阵的精确计算方法

实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法.     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

多肢剪力墙自振频率传递矩阵法简化计算

基于连续体系模型与传递矩阵法的结合,提出了一种简化方法来计算多肢剪力墙的动力特性。整个结构体系等效模拟为一个夹层结构梁;建立了相应的基本微分方程;通过求解微分方程得到了楼层间变形函数;利用边界条件及层间变形函数导出了层间传递矩阵;据此计算出了体系的振型与自振频率。算例分析表明本文方法计算结果与文献给出的解是比较吻合的。     

���α���ת����̽��

????????α???????????????????????????????????????????????λ???????????????????????????λ??????????????t???????1/2????????????????t?????????????????????????????????????????????????????????????????????λ?????????λ???????????????????????????????????????????????????????λ??????????????????????????????????????????λ??????????????λ??????????????????????????????λ?????????????????????????λ?????????????????????????????????????????????????????????Ч?????????????     

基于MRPs无陀螺卫星姿态确定的抗差估计算法

针对精确航天器姿态问题,采用修正罗德里格参数(MRPs)表示姿态,用动力学方程进行角速率的传播,分别基于等价协方差和观测残差统计特性的协方差调整方法,实现UKF算法的抗差效果,并通过仿真验证该方法的有效性。仿真结果表明：在含有粗差的情况下,两种抗差UKF算法均能得到稳定的状态估计值,而第二种方法计算效率更高。     

地基沉降预测模型的正则化算法

通过分析地基沉降预测模型,指出最小二乘的病态性会导致模型参数求解失败。应用正则化理论,基于矩阵求逆理论,提出了一种沉降预测模型参数的正则化无偏估计算法,说明了新算法的无偏性和方差最小性。在一定条件下,证明了新算法中正则参数的存在性,并给出了正则参数的计算公式。结合文献和工程实例进行的分析表明,新算法降低了矩阵条件数,减轻矩阵病态程度,能有效求得地基沉降预测模型参数。