有限区域变分四维同化系统的设计和试验

使用共轭码技术，构造了中尺度模式的伴随模式，研究了模式侧边界条件。基于伴随模式技术的原理，设计了伴随同化和Derber同化方法两套方案；下降算法采用共轭梯度法，使用准真实资料，进行子同化数值试验和预报试验。结果表明，两套同化方案都取得了一定的同化效果，Derber同化方法比伴随同化更令人满意。证明了变分同化是提高预报准确率的一种非常有效的手段。     

唐山地震应变-电阻率前兆及虚错动模式

研究了唐山地震周围240 km内16个地电台及大同6.1级地震周围50～60 km内2个台地电记录的全过程，排除了干扰，在确定可靠前兆异常的基础上，研究了地电前兆场的分布规律，证明了地电前兆场与震源应力场之间存在物理联系．对比了同震电阻率变化符号分布与震源机制结果，表明同震电阻率变化与地电前兆变化符号相反，分布相似．因而认为唐山地震前后的地电观测表明，唐山地震是一个弹性回跳过程．提出了唐山地电前兆的虚错动模式，即设想地电前兆是由一个与地震时产生的实际错动符号相反的虚错动所产生的应变积累而引起的．结合所求得的10－7～10－5应变范围内放大系数K=(/)/的非线性特性，根据断裂力学理论，选取合理虚错动参数，计算了唐山地电前兆的理论分布．与实际地电前兆对比表明，半定量的理论值与实测值基本吻合，从而证明虚错动模式对唐山地电前兆是适用的.      

井地直流电法三维数值模拟中若干问题研究

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全Cholesky分解稳定化的双共轭梯度算法(ICBG),通过与均匀半空间垂直线源解析解的对比,证明了该算法是准确可靠的。     

积分方程法三维模拟井地电法并行算法研究

在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。     

直流激电测深二维反演的若干问题研究

目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.     

1991年夏季持续性异常环流的共轭敏感性分析

在共轭敏感性分析的基础上发展了最敏感扰动的概念。用共轭模式对江淮流域大水1991年东亚夏季环流的持续异常特征,包括乌拉尔阻塞高压、鄂霍次克海偶极子型阻塞以及副热带高压的变动过程进行了分析,分别得到了对这些过程影响最大的初始扰动。结果从新的角度揭示了这些系统最重要的影响因子,也显示了共轭敏感性分析方法在大气环流物理过程诊断中的有效性和应用前景。     

利用不完全Cholesky共轭梯度法求解点源三维地电场

点源三维地电场的求解是一大型数值计算问题．本文用有限差分方法求解,最后形成一个线性方程组Ax＝b,这里A是大型稀疏的带状对称矩阵．解大型稀疏方程组的完全Cholesky分解直接算法,一般要求巨大的机器内存来存储系数矩阵A,而且计算速度极慢．因此引入不完全Cholesky共轭梯度（ICCG）算法及按行索引的稀疏存储模式,充分利用A的稀疏性,使得计算速度大大提高,而内存要求则大大减少,因此ICCG算法是地电三维正演的强有力工具．     

鲁西北西向断裂系与沂沭断裂带晚中生代演化关系及其动力学背景探讨

晚中生代是华北地块构造演化的转折时期。由出露于沂沭断裂带、鲁西地体、鲁东地体的中生代地层、岩浆,结合断裂活动年代学、区域地质等资料分析,可以将沂沭断裂带晚中生代构造演化划分出距今约160Ma、130～110Ma、90～80Ma等3个关键时期,并分别与左行压剪、左行张剪、右行压剪构造活动相对应。晚中生代沂沭断裂带与鲁西北西向断裂系间的几何学、运动学、年代学的共轭匹配关系表明,它们为特定构造动力学背景下形成的一组共轭断裂系。同时对鲁东地体晚中生代构造演化、鲁东北西向断裂系特征,以及沂沭断裂带、北西向断裂系晚中生代构造演化的动力学背景进行探讨。     

地形对不同反演的影响

[摘 要]根据观测数据和模型参数之间的关系,将求误差函数的最优问题的方法通常分为线性反演和非线性反演; 线性反演要求有相当准确的背景介质模型的条件,这和地球物理问题中的实际情况相矛盾。特殊情况下,线性反演甚至出现发散的情况。而非线性反演不用求目标函数的偏导数,全局随机搜索模型空间的最优解的办法,避免了问题的线性化,在山区复杂地形条件下,通过对比论证,非线性反演能够获取更合理的真实地电模型,从而说明由于非线性反演适应性较强、精度较高已被广泛的应用于地球物理的各个领域。     

BP神经网络在似大地水准面精化中的应用

利用分辨率为5'×5'的EGM2008模型计算高程异常,从水准联测点上移走模型中的高精度中长波部分后,通过BP神经网络模拟余下的短波部分,最后再对数据进行恢复,使实验达到较高的精度。通过实例对确定高程异常的多种方案进行比较分析,得出了一些有益的结论。