最优化波数在点源二维有限元正演计算中的应用

文章结合地电模型给出了计算误差随最优化波数个数的变化曲线,揭示出该曲线上存在突降稳定点,这一规律使计算所需波数大大减少。研究还发现误差随网格加密参数曲线上也存在突降稳定点。上述规律的发现可以大大节省二维有限元正演的计算时间,使之在反演中的充分应用成为可能。     

在波数域计算一维重磁异常导数的Matlab语言算法

利用Matlab内建的快速傅氏变换函数可以方便地在波数域计算重磁异常导数.介绍了基于Matlab语言的波数域求导的算法,给出了程序源代码,讨论了一些有助于提高计算精度的编程技巧.通过模型试验和数据分析,发现在计算垂向导数时波数域求导算法的精度比傅氏级数的精度有明显改善;而水平导数的计算,2种方法的精度相当.在某区钾盐勘探中,用该方法处理高精度重力剖面数据,取得了较好效果.     

二维视电阻率断面的快速最小二乘反演

快速最小二乘反演是以平滑限定的最小二乘方法为基础,是对二维视电阻率断面进行反演的一种方法。反演过程不需要提供初始模型,在首次迭代时使用一均匀介质地下模型作为初始模型,该模型的视电阻率偏导数值可以用解析法得到。在后面的迭代中,使用了拟牛顿法去修改每一次迭代的偏导数矩阵,避免了偏导数矩阵的直接计算,从而减少了计算时间和存储空间。同时运用牛顿矩阵校正技术解最小二乘方程组也减少了大量的计算时间。总之,该方     

三维电阻率正演计算中的Lanczos迭代算法

在三维电阻率的正反演计算中,快速、准确的正演计算是反演的关键。而正演计算往往涉及到求解大型线性方程组Ax=b的问题,通过Lanczos迭代构造出对称三对角阵方程组,并采用正交分解法进行求解,与传统算法相比,此算法占用内存少、收敛速度快、且稳定;针对大型稀疏矩阵的特点,采用简单地记录矩阵的非零元素值及其所在行、列值的方法,来存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。通过理论分析和点电源三维地电场计算实例,阐述该法是地电三维正演计算的有效方法。     

古洪水流量的误差计算

古洪水研究已在大型工程和长江三峡，黄河小浪底等的设计洪水中得到动用取得了令人虚心成功，但古洪水流量的推求与实测量的计算往往有大的差别，根据随机误差的传播给出了一般情况下的古洪不流量推误差计算公式，以小流底２３６０ａＢＰ古洪水注同误差的计算为例，对最终计算出的古洪水流量成果给出了一个误差范围。     

改进的局部波数法及其在磁场数据解释中的应用

局部波数法是一种进行磁场数据解释的常用方法。现有的局部波数法在进行反演时往往需要计算局部波数的导数,会明显地增大噪声的干扰,为解释结果带来误差。提出3种利用磁异常在不同位置或不同高度上局部波数的简单组合来进行场源体深度及构造指数的计算方法,不需要计算局部波数的导数,降低了噪声的干扰,增强了反演结果的稳定性。通过理论模型试验,证明改过的局部波数法在有无噪声的情况下均能很好地完成异常的反演,其反演结果与理论值之间的差距小于理论值的5%。将其应用于四川某地区磁异常的反演中,其结果与解析信号的欧拉反褶积法的反演结果吻合。     

混合离散优化算法反演应用

作者给出了一种新的混合离散优化算法ＣＤＯＰ,并用该算法反演了储层物性参数,取得了满意效果。     

改进射线追踪算法的微震源反演

基于射线追踪理论,提出了一种改进的试射法路径追踪算法,可以较准确地模拟射线在层状介质中的传播,不但可计算出任意检波点能接收到的所有的反射波的旅行时及射线路径,而且具有快速、准确、适应复杂构造模型的优点。利用共扼梯度法,通过正演模型的修改迭代,研究微地震层状模型的反演算法。最后,对实际的检测资料进行反演,反演结果基本上和实际资料吻合。     

界面位场异常的快速正反演计算及其应用

文章推导了快速计算界面位场异常正、反演公式，其原理及性质与R.L.Parker1972－1974提出的基本一致。但符合通常习惯，便于应用。推导途径与Parker有所区别，公式有所不同。迭代反演采取另一种形式，改进了收敛效果。此方法经理论模型验证，效果良好，并有应用实例。     

直流激电测深二维反演的若干问题研究

目前,激电数据的二维反演主要针对高密度激电数据,而对于常规的激电测深则很少涉及.据此,针对常规激电测深二维反演所涉及到的一些技术难点进行探讨.在正演过程中,根据实测数据实现有限元网格自动剖分,并采用双网格系统,用相邻节点电导率双线性连续进行有限元模拟,并对模拟的误差进行校正.在反演过程中,利用互换原理和Broyden的方法计算偏导数矩阵;采用与奇异值分解算法精度相当的共轭梯度算法求解偏导数矩阵方程.最后通过对模拟数据和实测数据进行反演试算,验证本计算方法是可行的.     

高密度电阻率法比值参数基于阻尼最小二乘反演

三电位电极系中装置的探测精度和数据处理方法的研究一直是地球物理工作者研究的一个热点,同时也存在很大的争议.利用正演模拟结果讨论了岩溶地区几种可能存在的地质条件下三种装置的探测精度,发现β装置和γ装置的探测效果明显优于α装置;然后利用正演计算得到的数据合成比值参数(T),对合成数据T进行最小二乘反演,发现T值反演结果和视电阻率反演结果一致,并在噪声影响较大的区域,T值反演结果优于T值等值线图,可作为判断异常体特征的一个依据,也可验证视电阻率的反演结果,弥补由于噪声对某种单一装置探测效果的影响.以义马某地的水文地质勘察为例,T值最小二乘反演结果表明,在含水低阻区域T值也呈现小值异常,且显示的异常体边界准确,结构特征明显.利用T值反演对数据处理具有重要的意义,应予以重视.     

A new potential function for the calculation of contact forces in the combined finite–discrete element method

Although the potential contact force proposed by Munjiza overcomes the difficulties inherent in the traditional discrete element methods, the physical meaning of the potential is not clear and the contact force derived from the original potential function is strongly dependent on the mesh configuration. In this study, we redefine a potential function and propose a new contact force calculation method based on a unified standard. Moreover, the new potential function retains all the advantages of the original potential function but has less mesh dependency. Copyright © 2016 John Wiley & Sons, Ltd.     

地学自由曲面径向基函数网络重建的详细算法

径向基函数(RadialBasisFunction,简称RBF)神经网络是一种理想的地学离散数据网格化工具,能够适应各种不同分布形式和边界条件的数据,收敛速度较快,可以逼近任何复杂曲面。这里详细介绍了RBF神经网络的算法。适当径向基函数的形式和偏差系数是使用RBF神经网络作地学曲面重建的关键。大量的实际数据验证结果表明,当选用Gauss型径向基函数时,一般可获得比较理想的网格化效果,同时具备"曲面平滑"和"拟合度高"的特点。     

利用AVO正反演技术识别深水凹陷气层

针对深水凹陷莺二段地质条件复杂,地震资料信噪比较低、烃检困难等难点,利用琼东南盆地陵水凹陷地气勘探阶段的钻孔资料及测井资料,应用zoeppritz方程进行射线数字模拟,研究该区典型气层及水层段的反射特征和AVO特征,对深水区坡折带低信噪比资料进行了保幅去噪及道集拉平处理,对常规的碳氢指示因子进行了改进.经结果表明:改进的碳氢指示因子对第三类及第四类AVO异常检测卓有成效,第四类AVO异常对于深水勘探有重要意义.     

离散KL变换在航磁轴向异常计算机检测中的应用

目前航磁异常的计算机检测手段单一,现有检测方法停留在定性研究的水平上,无法适应实际需要。为了克服这些局限性,将航磁异常计算机检测结果从定性水平提升到定量水平,这里提出了将图像处理中的离散KL变换方法应用于航磁轴向异常检测的新思想。经实地实例检验,效果良好,这一思想对研究航磁异常计算机检测有一定的借鉴意义。     

基于遗传算法的磁异常反演

遗传算法是模拟自然界生物进化思想所得到的一种全局非线性最优化算法,该算法可有效地避免局部最优解。这里介绍遗传算法的原理以及基于其算法的磁异常反演算例,对不同迭代次数所得理论模型反演的磁性体空间位置与磁性等参数的相对误差进行了对比。据此,对四川会理黎溪地区某磁铁矿区外围8号剖面上实测△T磁异常进行遗传算法,并迭代拟合反演,求其反问题,获得了与后续工程揭露及磁性测定结果相吻合的结果。     

利用多重网格方法研究波动方程正反演问题进展

多重网格方法是应用于大型复杂科学运算的高效率的方法，精度较高，它也可以解决波动方程的正反演问题。利用多重网格法可以提高波动方程正反演问题的运算效率和精度，还可以利用广义差分法或以差分方程替代微分方程，再结合多重网格方法求解波动方程正反演问题。     

K28边坡变形破坏机制的离散元模拟

在工程地质勘察的基础上,对K28边坡在天然和饱水状态下的变形进行离散单元法数值模拟,分析了此类缓倾角顺层边坡滑移压致拉裂的变形破坏机制,以及地下水促进边坡变形破坏的作用机理,为边坡的工程治理设计提供合理的建议。