地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

基于卫-卫跟踪观测技术利用能量守恒法恢复地球重力场的数值模拟研究

根据卫－卫跟踪观测技术的测量原理,基于能量守恒法建立了一种新的双星相互跟踪和三星相互跟踪的卫星观测方程. 通过数值模拟,采用预处理共轭梯度法恢复120阶地球重力场. 模拟结果表明：第一,双星相互跟踪恢复地球重力场的精度和美国喷气动力实验室公布的EIGEN GRACE02S的结果相符合;第二,三星相互跟踪恢复地球重力场的精度较双星提高约2倍.     

弯曲射线追踪中Dijkstra算法的改进与实现

文章针对图论中寻找最短路径的Dijkstra算法内存占用量大,效率低的缺点,对该算法进行了改进,修改后的算法计算效率是原来的四倍,内存使用量和图中节点数呈线性关系.在此基础上,用新算法求出了激发点和接收点的最短走时路径,并由激发接收点的旅行时结合联合迭代法对理论模型和实际场地进行了反演.结果表明：和直射线追踪相比,弯曲射线路径能更好地反演出地质体内部的速度场分布.理论模型和实际探测结果证实改进后的算法是有效的.     

各向异性扩散在电磁测深曲线平滑中的应用

电磁类地球物理方法由于极易受到各类噪声的干扰,使得估计的视电阻率曲线或相位曲线发生畸变,严重地影响了反演解译精度,如何对这类观测曲线进行合理有效地平滑,是目前数据处理中的重点,也是难点.本文把各向异性扩散（anisotropic diffusion AD）引入到曲线平滑中,提出了各向异性扩散的曲线平滑方法,以平滑点梯度值的降函数作为扩散速度,在梯度较大的位置予以较小的扩散速度以保护特征,在梯度较小的位置予以较大的扩散速度平滑噪声引起的扰动.同时为了减小“飞点”(outliers)对平滑的干扰,采用局部单调扩散进行预处理.实验结果显示,本文提出的平滑方法有效地平滑了噪声与“飞点”的干扰,恢复了曲线的基本形态,有效地保持了曲线基本特征.     

基于APSO-SVR-GRU模型的白水河滑坡周期项位移预测

滑坡周期项位移的预测,是研究地质灾害中滑坡变形至关重要的一步。由于单一模型易受偶然因素影响,且无法充分利用有效信息,导致其预测精度不高,适用性不强。基于此,文中提出了一种结合自适应粒子群算法（APSO）、支持向量机回归算法（SVR）、门控神经网络算法(GRU)的组合模型。该模型通过自适应粒子群优化算法对支持向量机回归算...     

固体潮观测数据预处理中的常用方法及计算机实现

从地震台站产出的原始地形变前兆数据出发,提出实际操作性较强的预处理方法,从已编制好的软件中给出处理实例,从而提高地震台站观测人员的科研能力,填补资料预处理方面的某些空白,对于地震台站人员集观测、科研于一体有着十分重要的意义。从95-01-02项目厦门地震台前兆数据预处理的实践表明:这些方法取得了明显的实效,直到今天仍然具有很大的使用价值。     

固体潮汐理论值算法在VB 6.0中的应用与实现

太阳和月球对地球的力学影响在地球科学研究中占有重要地位,但目前对该物理量的计算还没有人机交互很好的计算软件。本文在前人引潮位计算基础上,使用了面向对象的程序设计方法,编制了界面友好的固体潮汐理论值计算软件。介绍了该软件的设计思路和基本功能。该软件采用窗口界面进行输入和输出,数据显示以波形图为主。软件的友好界面增强了数据分析的直观性,为固体潮理论的推广应用和教学奠定了一定的基础。     

基于耦合场分割算法的ANSYS二次开发及其应用

耦合场分割算法是当前研究的热点之一。本文首先对比研究了瞬态热传导方程和孔压增长消散方程,提出了一种适用于ANSYS二次开发的分割算法,并对2种不同的预测值计算模式对计算稳定性的影响进行了分析。然后,据之对ANSYS温度—结构耦合场模块进行了开发,形成了岩土工程有效应力计算模块。数值实验表明,这种方法具有较高的计算精度和良好的计算稳定性,并且能够充分利用现有程序资源,具有较高的工程应用价值。     

实时子结构实验Chang算法的稳定性和精度

与慢速拟动力子结构实验相比,实时子结构实验的优点在于它能真实地反映速度相关型试件的特性。实时子结构实验的逐步积分算法通常借用拟动力算法,但是目前液压伺服作动器很难实现速度反馈控制,因而试件速度不能实现原算法的假定值,这样一来算法的稳定性和计算精度将发生改变。台湾学者S.Y.Chang提出一种无条件稳定的显式拟动力算法,本文分析了这种方法应用于实时子结构实验时的稳定性和计算精度。研究发现在实时子结构实验中该方法由无条件稳定变成了有条件稳定的,精度也发生了改变。     

新一代雷达应用软件对富宁5.16等特大暴雨的检验

云南文山天气雷达站采用的新一代天气雷达是由四创电子公司生产的C INRAD-CC,随雷达配备的应用软件是由北京敏视达公司提供的,由美国W SR-88D雷达Bu ild 7(第七版本)软件移植而来的软件RPG与PUP,通过研究W SR-88D降水处理系统的质量控制步骤,对降水算法的可调参数进行重新设置与验证,发现降水等部分原厂设置并不适用于低纬高原的雷达站所在地区,通过调整后的设置经过文山天气雷达站几次大雨与暴雨的比较后得出:对一些参数进行重新设置后其输出的降水产品比较接近降水的真实情况。