重磁遗传算法三维反演中高速计算及有效存储方法技术

将地下场源区域规则划分成很多小长方体单元，并且通过反演确定这些单元的物性变 化，勾画出场源的分布图像，这种方式逐步成为重磁反演，特别是三维反演的重要方向；遗 传算法等非线性技术进行该类反演将逐步成为发展趋势. 本文指出，在应用遗传算法进行该 类反演过程中，隐含着数据量较大时超常规的计算量，它已成为制约该类反演充分发挥作用 的瓶颈问题；同时，本文提出了针对性的分离并存储几何格架的计算策略、以及独特的几何 格架等效压缩存储技术，可以从根本上提高非线性反演计算速度，为该类反演的有效应用奠 定了坚实的基础.     

二维波动方程约束反演的大范围收敛广义脉冲谱方法

针对二维波动方程，将大范围收敛的同伦 方法引入算子参数识别反演的过程，并结合求解不适定问题的Tikhonov正则化方法，设计了二维波动方程反演的大范围收敛广义脉冲谱方法，并利用测井资料进行了约束反演. 大量 的数值模拟结果以及抗噪实验均表明了方法的有效性.     

上地幔变黏度小尺度对流的数值研究

基于二维模型，利用有限元方法，研究上地幔-岩石圈系统的变黏度小尺度对流. 考虑该系统的黏度随温度以指数形式变化，数值结果表明，当黏度随温度变化较剧烈时，由于低温高黏度，系统的最上部物质不参与对流，系统发育形成一个类似于岩石圈的静止盖层. 计算表面热流、地形起伏及重力异常与对流格局有较好的相关性，高热流、上升地形对应于对流的上升区，反之低热流、下降地形与对流的下降区对应.     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

根据激光雷达观测结果研究中国武汉地区钠层的分布

通过大量钠层荧光激光雷达观测数据的分析研究，给出我国武汉地区钠层平均分布形态的基本特征. 结合钠层模型分析，指出大气微量元素分布的纬度变化可能是导致不同纬度地区钠层平均分布形态差异的主要原因,同时给出武汉地区随季节变化的钠层平均分布形态. 对武汉地区钠层分布形态时间变化的分析，发现其短期变化与重力波活动及夜间变化与潮汐波活动相关.     

时间域航空电磁法一维正演研究

基于电磁勘探理论，导出了层状大地条件下时间域航空电磁法（偶极-偶极装置）的正演计算公式和算法，编制了相应的计算机程序，对若干典型地电断面作了正演计算. 计算结果说明时间域航空电磁法的探测能力和探测条件，进而为时间域直升机航空电磁系统的设计方案提供了依据.     

热层风场的理论模拟与分析

通过求解中性大气Navier Stokes动量方程建立了一个时变的三维风场理论模式，利用目前新版的中性大气模式NRLMSISE 00及国际电离层参考模式IRI2000作为输入参数给出热层风场. 基于该模式，计算得到中等太阳活动年磁静日风场的变化形态及其受电场和离子曳力的影响. 同时，将Navier Stokes动量方程作不同形式的简化，并利用简化模式与本文的模式计算结果的对比，分析中性大气Navier Stokes动量方程中黏性项以及非线性项(U·Δ)U的作用. 结果表明，本文所建立伪三维风场模式给出的结果更为合理，而简化模式在某些地区尤其在低纬和赤道区不适用，黏性项及非线性项的作用不可忽略. 本文所建立的风场模式将对研究电离层动力学过程、电离层与热层的耦合过程以及空间天气学研究都有着重要意义.     

1950～1990年中国地磁剩余场冠谐分析

根据中国1950、1960、1970、1980年和1990年地磁三分量绝对测量资料，使用球冠谐和分析方法，分别计算1950～1990年各个年代中国地磁剩余场冠谐模型. 球冠极点位于36°N和104°E，球冠半角为30°，冠谐模型的截断阶数为8. 地磁剩余场冠谐模型的均方偏差分别为：对于X分量，1950年为93.1 nT，1960年为128.9 nT，1970年为107.2 nT，1980年为107.6 nT，1990年为95.2 nT；对于Y分量，1950～1990年依次为74.8 nT，98.1 nT，89.2 nT，89.9 nT和84.0 nT；对于Z分量，分别为122.2 nT(1950年)，135.0 nT(1960年)，137.7 nT(1970年)，110.1 nT(1980年)和107.5 nT(1990年). 根据中国地磁剩余场冠谐模型和全球地磁场DGRF模型，得到中国地磁场的冠谐模型，并对冠谐模型的边界效应进行了分析和讨论.     

利用地震层析成像数据计算地幔对流新模型的探讨

假设地幔地震层析成像数据对应的地幔横向不均匀结构是地幔热对流的结果. 将地震层析成像数据转化为地幔温度(或密度)不均匀分布，考虑热流体动力学的三个基本方程，顾及热输运方程中的非线性项，直接将地震层析成像转化的地幔温度不均匀分布作为内部荷载引入基本方程, 反演计算地幔对流. 本文在利用地震层析成像数据计算地幔对流模型的新理论和方法的基础上，用SH12WM13地震层析成像模型数据，计算了全球地幔对流格局. 结果表明，对流格局不仅依赖地震层析成像数据，而且在很大程度上受地幔动力学框架、热动力参数和边界条件的所确定的系统响应函数的影响. 显示了地幔中复杂的对流格局，特别是区域性层状对流以及多层对流环可能在地幔中存在的现象.     

长白山天池火山区上地壳三维速度层析成像

利用三维地震透射观测系统所获得的莫霍界面反射波对长白山天池火山区下方上地壳的三维速度细结构进行了反演计算，得到了研究区下方地壳内不同深度的P、S波速度层析成像和VP/VS扰动分布图像，结果表明，天池火山区及其周边速度结构所呈现出的非均匀性变化与该区的地质构造和岩浆系统的存在有着密切的关系，异常变化比较明显的低P波速度异常体分布和高VP/VS扰动体意味这些介质体处在“软化高温”的状态，这与该地区广泛分布的地热异常和低Q值结构的研究结果是一致的.