快速排序法

在物化探数据处理中,经常需要对原始数据按从小到大(或从大到小)的顺序进行排序,这里将介绍一种新型的排序方法。快速排序法的基本原则是:通过分部排序来完成全部排序。作为子过程调用,快速排序法的算法流程如图1、算法PAD图如图2。     

叠前非局部平均滤波压制随机噪音

非局部平均滤波方法的去噪性能优异,但其在地震资料处理中的应用刚刚起步。该方法利用数据具有的结构冗余,以包含局部结构的小窗口或邻域为单元,利用局部结构相似性进行加权运算,增强有效信号,压制随机噪音。针对叠前地震资料数据量大、噪音背景强、局部结构简单;原始非局部平均算法对每一点滤波,需要对数据体内所有点计算权系数后进行加权计算,计算量大,对强噪音背景适用性差等不足,对原始非局部平均算法进行了改进,主要包括:基于速度谱的搜索窗口分割;基于梯度域奇异值分解的局部结构相似集选择方法;基于相似集大小的自适应滤波参数选择方法。试验结果表明,该方法改进后对于叠前地震数据的随机噪声具有较好的压制作用。      

大量数据排序算法的优化设计

本文结合工作实际需要，对大量数据的排序问题设计了优化算法，并给出了排序优化算法程序。     

结构面抗剪强度参数线性拟合方法适用性研究

针对岩体结构面力学参数的重要性,开展最小二乘法、点群中心法、随机-模糊法、可靠度分析法4种结构面抗剪强度参数线性回归方法的理论对比与分析。在传统点群中心法的基础上,提出求解过程严密、计算精度较高的改进优化算法。为保证可靠度分析法中验算点法参数初始赋值的准确性,提出基于可靠度理论的结构面抗剪强度参数c, 的均值、标准差优化赋值方法。以武隆县鸡尾山滑坡为例,比较不同岩体结构面数据条件下上述4种方法的适用性。结果表明,各方法的适用性与试验数据的离散性及数据量有关;当数据线性理想时,4种方法均可直接用于结构面抗剪强度参数的线性回归分析;当试验数据较离散时,点群中心法已不再适用,可根据工程要求、异常值情况和计算收敛情况选取合适的数据处理方法。     

裂隙岩体不稳定温度场的复合单元算法研究

基于复合单元法,结合三维热传导-对流方程和“充填模型”,提出了裂隙岩体不稳定温度场的复合单元模型。该模型前处理简便快捷,计算网格生成时无需考虑裂隙的存在,网格剖分不受限制,随后利用复合单元前处理程序,依据裂隙的位置和方位将其自动离散在单元内。对常规热传导-对流方程进行自伴随性调整,应用变分原理,推导出裂隙岩体不稳定温度场的复合单元算法,该算法可分别计算出岩块子单元和裂隙的温度值,且可真实反映裂隙中水流与相邻岩块间的热能量交换规律。将复合单元数值模型计算的不稳定温度场结果与相应的实测数据进行对比分析可知,数值计算结果与实测数据基本一致,验证了裂隙岩体不稳定温度场复合单元算法的可靠性与有效性。算例分析表明,裂隙中水流与相邻岩块间有明显的热传导和热对流作用。     

二维有限元网格的自适应剖分及程序实现

有限单元法在计算岩土力学领域已得到广泛应用。在有限元的分析计算中，前处理工作即准备计算数据耗时费力。如何简化前处理工作、实现网络自动生成对于推广有限元分析至关重要。作者在综合考虑多种算法原理的基础上，对非规则几何图形、复连通域的二维模型提出了一套新的自适应剖分算法。该算法可以同时生成三角形单元和四边形单元，还能对单元形状进行规则性调整。给出了详细的实施步骤，并编写了相应的应用程序，应用程序同时提供了对单元节点进行优化排序的功能。实际应用表明，能够高效快速地给出良好的剖分网格。     

数据不协调时独居石EPMA CHIME定年计算方法的改进

Th-U-Pb系统数据不协调是独居石电子探针化学定年(EPMA CHIME Dating)中一种很常见的问题。独居石矿物产生数据不协调的主要原因包括:1)蚀变或重结晶造成的铅丢失;2)不同年龄域在空间上的重叠或者存在于很小颗粒上的小年龄域。独居石EPMA年龄必大于U等于0时的极端情况给出的值,即当U为0时,EPMA CHIME年龄给出的是~(208)Pb/~(232)Th年龄,这是测量区域内最老年龄的下限。当Th为0时,EPMA CHIME年龄值介于~(206)Pb/~(238)U和~(207)Pb/~(235)U年龄值之间,这是EPMA法所能得到的最老年龄的上限。分析表明,当独居石EPMA数据出现不协调时,传统等时线方法计算的年龄值误差较大。本文提出了一种处理数据不协调情况下的优化算法。该算法考虑了测量误差,并根据剩余铀的总量剔出大的离散数据。利用已公开的数据进行算法对比的结果表明,本文提出的优化算法计算结果可靠。     

曲线拟合在核探测器信号幅度提取中的应用

研究数字化核探测器输出脉冲信号的幅度提取问题。根据闪烁探测器的电路特性，推导了其脉冲信号理想曲线方程。综合利用三和值法、最小二乘法、试探法等方法，设计了该理想曲线拟合算法。由脉冲信号采样点数据，通过曲线拟合，计算出曲线的极大值，该值即为闪烁探测器输出脉冲信号的幅度。正反演数据结果对比表明，采用曲线拟合算法提取闪烁探测器输出脉冲信号幅度，其方案可行，计算结果更准确。相较于数字滤波和数字脉冲成形方法提取核探测器输出脉冲幅度，该算法更简单，对仪器的硬件处理速度要求更低。该方法开辟了核探测器脉冲幅度提取的新途径。     

改进"BT"法在平原区潜水系统水位预报中的应用研究

针对潜水系统水位预报有限单元计算方法--"BT"法存在的问题,提出了改进方法--改进"BT"法,该法是在第一类边界水位的修正中,考虑水位下降与潜水蒸发减少的相互作用;该方法的关键是计算潜水蒸发减少引起的第一类边界水位下降减少值Se.用数值离散及解析公式相结合计算Se值的方法、原理、过程,并给出了"BT"法、改进"BT"法的计算实例及结果对比.     

一种岩性提取新方法--光谱排序编码法

光谱排序编码法(SSE:Spectra Sort Encoding Algorithm)是在对目标光谱与地物光谱按光谱强度大小进行排序的基础上,计算两条光谱曲线的相似度,然后,根据相似度与识别误差进行编码,从而达到识别目标的目的.文章详细介绍了算法提出的背景、原理,并对算法在岩性提取中的精度进行对比验证,结果表明在识别误差从0到0.3变化的过程中,算法最高精度的平均值为1,平均精度的平均值为0.825.高精度是其进行岩性提取的基础,岩性端元提取对比试验的结果也证明它是一种有效的岩性提取算法.