数值模式初值的敏感性程度对四维同化的影响

用著名的Ｌｏｒｅｎｚ系统作了共轭变分同化的数值试验。发现随着模式对初值敏感性程度的增加，用这种方法得到和模式相协调的初始场愈来愈困难，直到某些情况下的完全失败。这表明四维同化和可预报期限是联系在一起的。另一方面，随着方程不精确程度的增加，变分同化的效果愈来愈差，直到所做的预报无任何意义可言。如果在做变分同化的同时对模式参数也进行反演，就可使得基于Ｌｏｒｅｎｚ系统所做的预报效果大大提高。     

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

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

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

讨论了地下垂直线源分段计算和场叠加的方法,并实现了在套管上供直流电的三维数值模拟。讨论了大型容量矩阵的压缩存储方式,采用数组和结构体相结合的方法实现容量矩阵的一维链表式压缩存储。在求解超大型稀疏线形方程组时引入不完全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个关键时期,并分别与左行压剪、左行张剪、右行压剪构造活动相对应。晚中生代沂沭断裂带与鲁西北西向断裂系间的几何学、运动学、年代学的共轭匹配关系表明,它们为特定构造动力学背景下形成的一组共轭断裂系。同时对鲁东地体晚中生代构造演化、鲁东北西向断裂系特征,以及沂沭断裂带、北西向断裂系晚中生代构造演化的动力学背景进行探讨。