关于欧拉反褶积方法计算中的一点改进

在解欧拉方程的过程中,常用的方法是最小二乘法求解。但在求解过程中计算参数会根据需要做一些调整,每调整一次计算参数,就要重复计算一次解欧拉方程的过程,这样在很大程度上影响了计算速度。为了避免这种重复,这里找到了一组系数,该系数与场值及其梯度值有关,并直接建立起场源位置与构造指数N简单的对应关系,从而把解欧拉方程的最小二乘问题转化成了解线性方程的问题,避免了解欧拉方程组时的重复计算,大大地减少了计算工作量,有效地提高了工作效率。     

宁波城市地面沉降物理数学模型及沉降预测

根据宁波市地面沉降分层标历年观测数据,用不同的物理模型描述它们的压缩、回弹特征。将可能选用的沉降物理模型和分层标观测数据进行最优拟合,应用这种非线性规划(SUMT)的Powell最优化方法来求解宁波地面沉降的主要压缩层的变形数值,其结果表明宁波地面沉降处于主固结阶段。将宁波地面沉降数学总模型视为各土层变形之总和。各土层的计算主要的参数亦通过最优化方法进行反算求得。利用总模型所得的沉降计算值和同一地层结构区地面水准点实测数据进行比较获得满意的结果。将总模型外推则可预测同一地层结构区今后若干年内沉降数值。以上全部工作将由计算机完成。     

符号化计算法在地球物理应用上的初探——以计算大地电磁视电阻率及其偏导数为例

符号化计算法是基于符号运算的新型计算方法。分析了符号化计算法与数值计算法的区别,并以一具体示例实现为例,总结出其数值计算过程是按“创建需求母符号表达式”、“符号运算生成目标符号表达式”、“用真值替换符号计算出最终结果”三个步骤进行的。通过MATLAB编程对比了符号化法和怕特里克法数值法计算视电阻率及其偏导数的计算结果,证明了符号化计算法的可行性。实验也证明其数值计算速度较慢,但它具有为等效或近似数值计算提供标准结果、有助于发现新算法的优势,这种优势将使符号化计算法在地球物理领域中有广泛的应用前景。     

航空瞬变电磁数据一维Occam反演

分析了奥克姆(Occam)反演法的原理和计算方法,并将其应用于时间域航空电磁一维反演中,给出了一种简便的正则化参数搜索方法,对迭代过程进行了改进,大大提高了反演速度,5次迭代约20s。在反演初始模型中,电阻率数据由原始感应电动势数据通过自动迭代法转换而来,深度数据由初始电阻率计算得出。经算例表明,收敛过程中满足拟合差最小的拉格朗日乘子μ值随迭代次数增加而递减;反演过程非常稳定,一般5次迭代至10次迭代就能达到收敛,对低阻层的反演深度和电阻率值都接近真实值,但对高阻层的反演深度和电阻率值误差较大。     

利用数值法和解析法对祁南煤矿矿坑涌水量进行预测分析

以祁南煤矿为例,详细分析了其水文地质条件,利用两种计算方法(数值法和解析法)对矿坑涌水量进行了预测研究,并将观测数据与计算结果进行了对比分析,总体认为数值法得到的预测结果比较可靠。另外,在计算矿坑涌水量的过程中,解析法可以当做数值法的补充验证方法。两种预测方法相互结合,计算得到的矿坑涌水量更加可靠。     

非保守示踪剂条件下的河流水质参数的确定

将描述示踪剂为非保守、瞬时投源条件下的一维河流水团示踪试验的解析表达式进行2次适当的变形,分别得到一个二元和一元线性方程。对二元线性方程应用相关系数极值法,可以推导出计算河流横断面平均流速的公式;利用直线图解法或一元线性回归法能够计算出一元线性方程中的2个常数值,据此可以计算河流的纵向离散系数D_L与示踪剂的一级反应速率常数k。与现有的数据分析方法相比较,该方法具有:(1)适用于示踪剂为非保守的情况;(2)利用一组野外试验数据可以同时计算河流离散系数D_L、河流横断面平均流速V与示踪剂的一级反应速率k常数3个参数值;(3)全部数据分析过程可以程序化,由计算机完成等特点。     

纯球坐标系内各项重力校正值计算方案和计算过程

在国内,首次阐述了"纯球坐标系内各项重力校正值的计算方案和计算过程".该计算方案有如下优点:①以观测点P为中心,以球面距离am=am-1+Δa·bm-1(m=2,3,...,M)为半径将计算区划分为M个环带,再以Δλ0为由P点出发射线的角度间隔将所有环带等分为N (整数)份的计算区划分方案是最严密拼接,因而校正值的计算精度相当高.②在地形复杂和相对高差大的情况下,通过减小Δa、b和Δλ0 数值,可进一步提高计算精度.③可以直接采用带有经、纬度坐标的高程数据,进行大面积测区的各项重力校正值计算,而不受6°带和标准图幅限制.     

几种多元统计分析方法的综合通用程序

在用电子计算机处理地质数据时,往往需要对同一批数据采用各种方法进行计算。但在计算中多次输入程序和数据,既浪费机时,又增加了外部设备出错的可能性。这里介绍一种综合程序,它将常用的六种多元统计方法的程序编在一起,只要输入一次程序和数据,就可以反复计算。 这个程序是用TQ—16机BCY语言编写的,各种方法所用的共同计算步骤编成一些过程,计算时分别调用,篇幅不长,使用比较方便。     

切面截距加权平均网格化算法

本文介绍了切面截距加权平均网格化算法的方法原理.这种算法能较好地反映随机点的数据的趋势特征,适用广,计算速度也较快.文中通过实例将它与按距离平方倒数加权平均法和局域加权最小二乘拟合法作了比较。     

基于CUDA技术的矿产储量计算

插值计算是矿产储量计算中的重要环节,但由于计算量大,经常需要在插值密度和计算时间之间进行权衡.本项研究通过最新的CUDA技术,将插值运算转换到显卡上进行并行化计算,极大地提高了计算速度.以紫金山铜矿数据为研究对象,进行了方法技术的实验,与基于CPU的计算相比,基于GPU下CUDA的计算速度提高了67倍.     

快速排序法

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

GIS支持下计算格网自动生成技术

计算格网生成是流体力学数值模拟过程中的一项主要工作,以往普遍采用手工剖分方法,其结果对数值模拟产生很大影响。将地理信息系统(GIS)技术引入格网生成过程,利用计算机自动生成计算格网,并在此基础上讨论了以GIS为主控过程的集成数值模拟环境及其功能。     

基于CUDA技术的矿产储量计算

插值计算是矿产储量计算中的重要环节,但由于计算量大,经常需要在插值密度和计算时间之间进行权衡。本项研究通过最新的CUDA技术,将插值运算转换到显卡上进行并行化计算,极大地提高了计算速度。以紫金山铜矿数据为研究对象,进行了方法技术的实验,与基于CPU的计算相比,基于GPU下CUDA的计算速度提高了67倍。     

基于CUDA技术的矿产储量计算

插值计算是矿产储量计算中的重要环节,但由于计算量大,经常需要在插值密度和计算时间之间进行权衡。本项研究通过最新的CUDA技术,将插值运算转换到显卡上进行并行化计算,极大地提高了计算速度。以紫金山铜矿数据为研究对象,进行了方法技术的实验,与基于CPU的计算相比,基于GPU下CUDA的计算速度提高了67倍。     

P—Ⅲ型分布离均系数Φ值的高精度高效数值算法

本文提供了水文分析计算中常用Ｐ－Ⅲ型分布离均系数Φ值的高精度高效数值算法，并将该方法计算结果同其它方法成果进行比较与检验，该方法对偏态系数Ｃｓ和概率Ｐ适用范围之宽，计算速度之快及所能达到的计算精度足够满足实际应用的要求。     

饱水白云岩临界点、骨架和流体弹性参数的数值计算

基于临界孔隙度模型,提出了用岩石的整体弹性信息反演求取临界点、流体和骨架等局部弹性参数的数值计算方法和双向线性回归计算公式;结合饱水白云岩的样品测试数据,以孔隙度为自变量和因变量,对密度和密度与纵、横波速度平方的乘积分别进行了数值计算。以测试样品的整体信息求得其临界点、流体和骨架弹性参数值,并与实测数据做了相关性分析,其相关系数高达90%,充分表明数值计算公式的正确性和实现方法的有效性。岩石骨架、流体弹性参数的数值反演计算在油气勘探领域中具有巨大潜力,运用测井曲线和地震数据,可以反演求出岩石孔隙中流体弹性参数(密度、速度),对直接指示油层、气层起到重要作用。     

一种改进步长自调的非连续变形分析法

非连续变形分析（DDA）是一种针对块体系统变形和位移求解的数值计算方法。引入Newmark方法于结构动力学微分方程中,考虑惯性力和阻尼力作用,改进时间步长自动调节,并实现DDA求解程序;比较研究Newmark方法中的线性加速法、常加速法和平均加速法在DDA程序中计算的收敛速度,讨论块体系统动力学计算过程中DDA方法对惯性力和阻尼力的添加和删除,并提出根据计算精度要求的误差控制实现方案。将改进的DDA方法模拟一个典型的煤与瓦斯突出过程,取得了满意的计算结果,该改进算法为DDA方法处理动力学问题提供新的途径     

超大数据量的快速排序法

在超大数据量情况下,传统的排序方法及分级优化化法由于其自身的算法限制,难以更有效地提高排序效率。为此,作者在本文中提出了一种快速排序法,该方法的基本思路是用特殊设定的不同分界值把整个数据区域成尽可能多的小单元,每个数据通过与不同分界值的对比迅速归位于相应的单元内,然后对各单元内的数据进行排序,按序边接各单元,最终达到排序的目的,与分组优化法的计算工作量大小对比及试算结果,充分表明了该方法的算法可行,计算量少,排序速度更快等优点,特别适用于超大数据量情况的数据排序。     

降水空间不均匀性对径流过程模拟的影响

应用数字高程流域水系模型，通过子流域和网格2种空间离散方式，采用传统的子流域法、雨量插值子流域法、逐个网格法3种雨量数据输入处理方法，分别作日流量模拟与洪水过程模拟，模拟过程与实测水文过程相比较表明：日模中3种雨量处理方法计算精度相当，次模中逐个网格法优于传统子流域法和雨量插值子流域法；降雨空间分布不均匀性考虑得越充分，水文过程模拟精度越高。     

P—Ⅲ型分布Φ值数值计算方法比较

皮尔逊Ⅲ型频率曲线的数值计算方法是水文频率计算中的难点之一。本文对近年来提出的几种数值计算方法进行了分析比较，并推导出了各数值计算方法试算离均系数Φ值的迭代公式。最后推荐计算精度高、运算速度快、收敛范围广的级数展开法迭代公式。