自适应非线性有限元的多重网格法求解及应用

描述了非线性自适应有限元的多重网格法求解过程。 通过算例说明了方法的有效性以及自适应多重网格有限元用于开挖施工过程模拟分析的优越性。     

在AutoCAD200O中应用VBA技术自动绘制等高线

应用AutoCAD进行地形图的计算机辅助制图能极大地提高工作效率，但因其没有根据数字地形分析（DTM模型）自动产生等高线的功能，用VBA编程技术，将外业采集的离散点通过建立三角网，在三角形边上内插等高点的方式，能在AutoCAD中实现等高线的自动绘制。     

伏特拉核函数模型在红岩煤矿矿井涌水量预测中的应用

本文运用伏特拉核函数基本理论,建立了南桐矿区红岩煤矿矿井涌水量的非线性核函数模型,介绍了计算核函数模型所用的递归算法。同时对模拟结果的残差进行了分析,并用ＡＲ模型进行了改进。     

线性方程组迭代解的随机模型测试研究

本文讨论大型线性方程组迭代解的随机模型测试评价问题。给出了常用迭代解法ＣＧ、ＬＳＱＲ、ＳＩＲＴ、ＳＡＲＴ、ＳＡＳＩＲＴ等的测试结果。结果表明：（１）方程组系数矩阵的特性（条件数）及解结构都对解精度有重要影响。解模型越粗糙,解的精度越低。（２）各种求解算法都有一定的平滑效应,同时各种算法也都会产生误差大于２００％的奇异解,奇异解元素数一般约占１０％。（３）数据的拟合残差一般不能真实反映解的精度。（４）对含误差数据的求解问题,较好的求解算法是ＤＬＳＱＲ与ＳＡＳＩＲＴ。     

基于L-BFGS算法和同时激发震源的频率多尺度全波形反演

全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。     

LBG算法对初始码书敏感的实验性能分析

介绍矢量量化的基本概念及其算法—LBG算法。指出LBG算法的缺陷，并针对LBG算法对初始码书敏感这一缺陷进行分析。同时还分析了LBG算法对几种不同初始码书的实验性能，得到了一种度量初始码书性能“优、劣”的方法，并提出了一种性能优良、操作简便的初始码书生成方法—等间隔法。最后，针对LBG算法对初始码书敏感这一缺陷给出了一些改进方向。     

土质边坡失稳判别的CA-AA耦合模型

解决边坡稳定性问题的关键是确定最危险滑动面(潜在最危险滑动面)和边坡的稳定系数。将蚁群算法用于求解土质的边坡稳定性问题。针对蚁群算法在演化过程中存在停滞和过早收敛的现象,对基本蚁群算法的结构形式和蚂蚁转移概率的计算进行了改进。构建了一种新的蚁群算法——元胞蚁群算法(CA-AA)。在文武坡边坡稳定性判别中,利用CA-AA法计算得到其稳定系数为1.37,略小于《岩土工程勘察规范》不平衡推力法计算值,弥补了不平衡推力法的不足,更接近实际。算例结果表明:元胞蚁群算法(CA-AA)总能搜索到问题的全局最优解,且搜索效率也有较大的提高。     

方差体技术在地震勘探中的应用

三维地震数据体反映了地下一个规则网格的反射情况。当遇到地下有断层或某个局部区域地层不连续变化时，一些地震道的反射特征就会与附近地震道出现差异而导致地震首局部的不连续性。三维方差体技术就是求取三维数据体所有样点的方差值来反映这种不连续性信息。通过该技术在济宁2号井某采区的实际应用，可看出，方差体技术在三维地震信息的自动拾取以及提高地震资料解释成果精度等方面有明显优势，对断层、陷落柱有良好的自动识别能力。     

基于耦合热传导-冰流模型的冰川钻孔温度重建的气候反演算法比较

冰川内部温度与冰面温度变化过程密切相关,因此可由实测冰温-深度剖面重建冰面温度变化。耦合热传导-冰流物理模型及其相关反演算法是基于冰温-深度剖面重建冰面温度变化的理论基础和关键。通过构造理想条件下的数值实验,分析比较了最小二乘法、Tikhonov正则化方法和蒙特卡罗算法重建百年尺度冰面温度变化的实验结果。同时,以我国藏北高原腹地的马兰冰川钻孔资料为例,结合气象资料构造真实钻孔数据模拟实验。分别给出了三种反演算法重建10 a和40 a尺度的冰面温度变化结果,并讨论了不同算法的优劣性和适用性。两组模拟实验结果均表明:Tikhonov正则化方法是目前求解该反问题的最优方法。与其他算法相比,Tikhonov正则化方法在各个时间尺度的重建结果与真实冰面温度变化吻合更好。另外,分别对重建百年和40 a尺度的冰温-深度剖面添加±0.02℃和±0.01℃的随机误差,研究发现Tikhonov正则化方法能够有效降低噪声干扰,相较其他算法得到的解更加稳定。这在一定程度上解决了该反问题中解的稳定性问题。     

L1-L2范数联合约束稀疏脉冲反演的应用

稀疏脉冲反演实际上就是利用反褶积原理, 从带有噪声的地震道中计算出具有稀疏分布特征的反射系数的振幅和时间.稀疏脉冲反演是非线性优化问题, 通常都是把非线性优化问题转化为线性优化问题, 然后用线性优化算法求解.以范数约束为基础, 提出L1-L2范数联合约束求解的方法.该算法采用了目前国际流行的内点算法, 与传统的优化算法相比, 这种算法具有精度高, 速度快的优点.通过研究人工模型和南海某油田实际数据, 表明该算法的计算结果和测井记录吻合好, 提高了地震分辨率, 目的层段分辨率优于8m.利用反射系数剖面预测的储层厚度和开发井吻合很好, 大大地降低了资源量计算的风险和油田开发的不确定性.