变异函数模型对渗透系数克里格插值的影响研究

在克里格计算或其它地质统计学研究中,变异函数模型选取的适当性直接影响渗透系数的结构分析与插值结果的合理性。分别基于高斯模型、球状模型、指数模型对同一沉积相介质的渗透系数进行结构分析,并采用普通克里格法进行插值,结果分析表明:变异函数模型对渗透系数克里格插值的影响主要体现在变程值的不同;变程值过小,搜索邻域内各数据点之间的相关性降低甚至消失,克里格法退化为简单的算术平均;随着变程值增大,参与克里格计算的实测点增多,插值结果趋近于平稳;当变程值与实际值较为吻合时,克里格插值精度达到最大。     

基于三次插值的大地电磁自适应有限元二维正演模拟

大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。     

基于改进BP网络算法的隧洞围岩分类

围岩分类对指导地下工程的设计和施工具有非常重要的意义.引入人工神经网络的方法, 进行隧洞围岩分类, 在传统BP算法的基础上, 通过改进学习算法、优化传递函数和网络结构进行神经网络方法优化.采用附加动量法和学习速率自适应调整的策略改进学习算法, 使得当误差大于上临界值时, 则降低学习率, 当误差小于下临界值时, 则适当提高学习率, 这样可加快网络的训练速度, 确保网络的稳定性; 通过引入调整学习率参数, 使得传递过程更加敏感, 加快了传递函数的收敛速度, 提高了训练函数的计算精度; 通过给定隐含层节点模型的取值范围, 对网络结构进行优化, 提高了泛化精度.将改进的BP网络模型应用于广东省东深供水改造工程的隧洞围岩分类中, 分类结果与根据《水工隧洞设计规范(SL279-2002) 》的分类结果完全一致, 表明该方法具有良好的工程实用性.      

无单元法布点优化的一种新方法

研究了无单元法中单个影响域内的布点形式对误差的影响,提出了半径权值的概念。在经过充分的数值计算后,给出在半径权值对无单元法插值结果的精度的影响值的大小。研究对优化无单元法的节点布置、减小模拟误差有很大的作用,并且使得影响域的大小不再依靠点数的多少来确定,避免了影响域过大的情况和为了单个影响域内有足够多的插值点而使得整体插值点数目过多的弊端。运用单个区域内的插值点健康度的理论,研究了整个区域布点对无单元法模拟函数--滑动最小二乘法模拟精度的影响,并运用遗传算法对整个布点区域上的插值点进行了优化,给出了区域上的优化值。通过一个实例,得到一个具有普遍意义的插值点优化值与插值点密度的关系,并给出无单元法布点的一般性原则。     

无压渗流问题分析的多节点有限元方法

采用基于平均值插值的多节点有限元方法分析有自由面渗流问题。在自由面附近采用多节点单元逼近自由面,利用平均值插值建立多节点单元的形函数,在远离自由面的区域采用四边形/三角形单元剖分。给定一个初始自由面位置,通过对渗流控制方程的多节点有限元求解,根据自由面上节点水头值判断自由面节点的调整方向和大小,最终迭代求出自由面的位置。土坝渗流问题的数值分析表明了所提方法的有效性和足够的计算精度。     

点源三维地形视电阻率的计算

利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。     

基于改进的双线性旅行时插值的三维射线追踪

旅行时射线追踪的精度和效率一直是影响地震层析成像质量和效果的关键因素。这里在三维近地表速度结构层析成像应用中,结合旅行时插值算法与最短路径算法,提出了一种三维初至波射线追踪算法。通过对双线性旅行时插值算法进行改进,并运用判定条件与简化插值计算公式进行快速计算,减少了插值次数,降低了运算量,同时也保证了较高的射线精度,有效地解决了三维射线追踪算法的计算低效率,射线精度不高的问题。     

地球物理勘探中几种二维插值方法的误差分析

为了研究采样和网格化方法对地球物理数据成图精度的影响,为野外数据采集布设提供一定的依据,采用数值模拟确定重力异常场场值,通过不同采样间距和不同插值方法计算重力异常绝对误差均方根值和节点处的绝对误差值,对比不同插值方法的误差,得到了如下认识：1）对于同一插值方法而言,存在小间距绝对误差均方根值小于大间距绝对误差均方根值的关系。2）对不同的插值方法而言:当采样间距小于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、反距离加权插值法、最近邻点法、最小曲率法,并且线性插值三角网法与自然邻点法具有几乎相同的数值;当采样间距大于最小异常地质体尺度时,绝对误差均方根值由小到大的顺序是径向基函数法、改进的谢别德法、克里金插值法、自然邻点法、最小曲率法、最近邻点法、反距离加权插值法,并且线性插值三角网法和自然邻点法具有几乎相同的数值。3）从绝对误差均方值看,径向基函数方法、改进的谢别德方法和克里金方法数值较小,其中径向基函数值绝对误差均方根值最小。4）从节点处绝对误差值来看,径向基函数方法、克里金方法、改进的谢别德方法相对其他插值方法具有更小的误差,不存在局部误差较小或较大的情况,是相对较好的插值方法,并且径向基函数方法是最好的。     

一种快速建立层析反演初始模型的方法

建立近地表层析反演模型时,通常利用射线追踪方法来计算模型初至旅行时,再与观测的初至实际旅行时进行比较,反复修改模型直到计算值与观测值的误差达到最小。因此,初始模型建立的质量直接影响到层析反演的精度、收敛速度和迭代次数。通过模拟近地表模型,利用时间项反演方法快速建立初始模型,提高了层析反演的精度,加快了收敛速度,减少了迭代次数。      

旅行时线性插值射线追踪提高计算精度和效率的改进方法

旅行时线性插值（LTI）射线追踪算法是基于线性假设的,在向前处理过程中仅用按行或按列扫描的方法来计算节点旅行时没有考虑逆向传播射线,导致其计算精度与网格剖分大小有关,在处理复杂介质时会使得追踪出来的射线路径不一定满足最短旅行时。因此,笔者提出了两点改进措施：在向前处理时需采用全方位循环的方法来计算节点最小旅行时;在网格边界加入次生节点。模型试算结果表明：采用全方位循环的LTI法考虑了来自各个方向的射线,可提高其对复杂模型的适应能力;在节点间距相同的情况下,网格边界插入次生节点的LTI法较传统的LTI法计算精度至少可以提高一个数量级,同时,计算速度也更快;随着节点间距剖分的越精细,计算耗时下降也越明显,计算速度较传统的方法可提高n~10n倍。     

独龙江公路下段Ⅰ号滑坡稳定性分析

独龙江地区岩体破碎，暴雨强度大，边坡稳定性差。独龙江公路下段Ⅰ号滑坡体位于贡山县普拉河旁，为确保公路安全，我们采用抗滑力／下滑力的方法来分析滑坡稳定性。     

梯度法解释复杂二维断裂重力异常

对以往的二维断裂重力异常的梯度解释方法作了几点改进:①采用最小二乘拟合方法求取台阶倾角α;②直接利用台阶倾角与台阶的水平边界位置来计算台阶的埋深;③将重力二维梯度解释方法推广到复杂断裂模型的解释。文中详细论述了这种方法的原理, 给出了其具体实现步骤。该方法比起以往重力梯度解释方法, 应用更加简便而且解释结果更加稳定。用该法成功地解释了某油气盆地的实测重力资料。     

在油气初次运移研究中压实曲线的应用

作者提出解释压实曲线的动态分析法。用它研究压实曲线,能够获得许多不同于静态分析的结论,如排烃深度、排烃方向以及对生油岩的评价等等。本文利用这种方法研究了板桥凹陷初次运移的问题,获得了与地化和构造研究一致的结论。说明这种方法在实践中得到了初步证实。     

区域稳定性的神经网络分区评价

以黄河黑山峡河段区域稳定性研究为例 ,本文系统地提出了区域稳定性的神经网络分区评价的程序与方法。结果表明 ,这一评价方法比其他方法更简洁 ,也更准确     

火山岩地区勘探地下水过程中自然电位异常的解释

当利用自然电位法（ＳＰ法）在火山岩地区勘探地下水时，发现负的自电异常大小与地下非饱和区厚度之间存在一个线性关系，从理论上对该经验关系进行了探讨，并提出对这一水文地质问题的自电测量解释方法，最后给出法国中央高原ＬａＣｈａｉｎｅｄｅｓＰｕｙｓ地区的勘探结果。     

面内云纹分析应变符号判断的两种新方法

提出了面内云纹法应变符号判断的两种新方法－零级法、同级法。零级法是根据零级云纹所处云纹图的特殊位置，如一角、一端来判断应变符号的方法；同级法是根据相同级数云纹间所夹云纹即为αｕ／αｘ，αｕ／αｙ（αｖ／αｘ，αｖ／αｙ）都为零的线来判断应变符号的一种方法。这两种方法分别在纯剪切和山字型构造应变分析中得到成功的应用。     

微机整谱拟合方法在矿物红外光谱多组分定性和定量分析中的应用

本文提出了一种利用微机整谱拟合技术对多组分物相的红外光谱进行定性鉴定和定量分析的方法,阐述了该方法的原理和计算方法。通过与其它红外光谱定性和定量分析方法的对比,该方法不需选择特征峰,不需制备更多的标样。每种物相只需一个含量(重量)已知的标样即可分析。它能够一次性地对多种组分物相做出鉴定并提供其中每种组分的含量。该方法完全实现了计算机程序化,减少了很多人为的因素影响,提高了分析精度,并达到了方便实用的要求。     

钢筋混凝土桩网复合地基沉降计算方法优选

钢筋混凝土桩网复合地基是一种新的地基处理技术。本文对适用于钢筋混凝土桩网复合地基沉降计算的方法进行探讨。     

有序地质量的简单分段法

近年来出现了用数学方法处理地质工作中的分层、分带或分段问题。费歇分段法(W.D.Fisher,1958),即“最优分割法”被普遍采用。应当指出,这个方法有一定的局限性。笔者曾对这个方法作了一些分析讨论¹⁾。为了克服费歇法的不足,笔者在本文中提出了另一种算法,称为简单分段法。     

岩石非弹性应变恢复(ASR)地应力测试方法中柔度研究进展评述

本文通过分析国内外非弹性应变恢复(ASR)法地应力测量研究现状,并对ASR法中重要的柔度参数作了系统整理,认为该测量方法在理论上逐渐趋于成熟,在技术应用上有广阔的前景,但该地应力测量原理中,针对温度、应力环境(作用应力)对ASR岩石柔度影响的研究相对较少,还存在一些缺陷,为此提出了该地应力测试方法有待改进或补充的建议,为完善ASR法地应力测量原理指明了研究方向,对今后ASR地应力测量技术的完善和发展具有重要的意义。