可控源音频磁大地电流法一维正演及精度评价

研究发现，两种目前比较常用的可控源音频磁大地电流法（ＣＳＡＭＴ）的正演计算方法，即简化积分法和二次插值法之间计算结果差异较大。为了正确地评价这两种方法的计算精度，本文提出了利用直接积分进行，ＣＳＡＭＴ一维正演的算法，即通过将任意角度频率测深正演结果沿发射源方向直接进行数值积分，得到有限长度发射源的ＣＳＡＭＴ正演结果，大量计算和精度评价发现，在这些算法中，简化积分法和本文提出的直接积分法计算结果吻合     

MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用

在矩形有限元和三角形有限元方法大地电磁测深（ＭＴ）二维正演计算的基础上，设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序，使得ＭＴ二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深（ＣＳＡＭＴ）ＴＥ极化方式的二维正演计算。文章给出了在不同模型上计算的ＭＴ和线源ＣＳＡＭＴ的响应，分析了计算ＣＳＡＭＴ响应时，计算点离源的距离对视电阻率和相位值计算的影响，比较了ＭＴ和ＣＳＡＭＴ网格剖分的特点，提出了计算ＣＳＡＭＴ响应的改进的边界条件设计     

线电流源声频大地电磁测深的二维正演计算及响应特点

论述了如何在大地电磁测深二维正演计算的基础上，结合线电流源声频大地电磁测深（ＣＳＡＭＴ）的边界条件，利用三角形６节点有限元法计算ＴＥ极化方式ＣＳＡＭＴ二维响应的方法；给出了４个模型的计算结果     

张量CSAMT数据的二维反演

在一单行本中,我们提出了一种ＣＳＡＭＴ正演模拟算法并说明了这种算法在解释野外数据中的应用效果,在本文中,提出了一种用于二维电阻率模型的ＣＳＡＭＴ数据反演算法,这种算法包含一个三维源并可以反演（ＴＥ）模式及（ＴＭ）模式下的数据,用这种算法对假设的数据进行试算,结果表明,不用做平面波校正就能对远区场,过渡区场及近区场进行反演,最终的模型结果显示,原始电阻率模型可被恢复,而不存在任何由上源于效应而产生的     

线源频率电磁测深二维正演(二)

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

线源频率电磁测探二维正演（二）

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必须要进行三维情况下的研究。线源频率电磁一维天演模拟研究为最终实现三维正演打下了基础。     

可控源音频大地电磁法在矿体定位预测中的应用--以山东省乳山市蓬家夼金矿床为例

可控源音频大地电磁法（ＣＳＡＭＴ法）是一种直接探测地下岩性及矿产的电磁测深新技术,在山东省乳山市蓬家夼金矿床的应用研究表明,ＣＳＡＭＴ法可以成功地用于确定矿体的形态,产状和埋藏深度,为矿床地球物理找矿模型的建立提供了科学的依据,是了解矿化规律和进行隐状矿体预测的有效手段。     

二维CSAMT三角部分人机交互解释程序的实现

二维ＣＳＡＭＴ人机交互解释程序，以往都是针对规则模型来实现的。本文主要介绍二维ＣＳＡＭＴ三角剖分的人机交互解释程序的实现过程；程序实现的软硬件环境，运算使用的参数以及程度运行操作方法。     

微机网络并行计算及2.5维CSAMT正演的并行实现

本简要介绍了如何建立一个简单、易行的微机网络并行计算的开发环境。在该并行环境上完成了２．５维ＣＳＡＭＴ有限单元正演的并行程序，并且获得较好的速比。     

复电阻率法二维正演并行计算研究

复电阻率法二维正演采用二次场算法,结合有限单元法进行离散化,然后引入ColeCole模型,完成数值模拟。通过理论模型试算验证了二维正演程序的正确性。当模型剖分较细时,多发射源、多频率的正演计算效率有待提高。复电阻率法二维数值模拟是按发射源或频率计算,且各发射源之间、各频率之间计算是相互独立且互不影响的,因此通过引入MPI并行计算来提高其计算效率,从而减少其正演所需的时间。采用串行、并行复电阻率法二维正演程序来计算同一理论模型,结果验证了该并行程序的正确性。对并行算法的效率进行分析表明,该并行算法是可靠的、高效的。     

边坡框架梁文克尔简化计算方法及工程应用

为了改进边坡框架梁的计算方法,基于文克尔地基梁假定,考虑框架梁与土体的互相作用,通过关键点拟合的方法,把指数函数和三角函数简化为一元二次函数,推导出边坡框架梁的简化计算公式。该方法被应用于福建一处边坡支护工程,采用锚索框架梁支护,施工过程中对边坡进行位移监测。研究结果表明,传统的文克尔地基梁计算需要求解微分方程,计算过程比较复杂,不适合设计中使用,而简化后的计算方法容易掌握,可以方便计算出框架梁的端部弯矩和跨中弯矩。边坡位移监测结果表明该简化计算方法安全可行。研究结果对边坡框架梁的工程设计有一定的参考作用。      

竖向荷载下桩基础弹性分析的改进计算方法

将Randolph剪切位移方法中桩身位移与桩端位移的函数关系简化为一多项式,并与Poulos积分方程中土体柔度系数矩阵相结合,提出了一种竖向受荷单桩弹性分析的改进计算方法,从而避免了Poulos积分方程法中的差分运算以及由此带来的其他矩阵运算,同时比Randolph方法更能准确模拟桩身剪切应力的分布,并将单桩的改进计算方法应用于群桩分析。将改进的计算方法与Poulos、Randolph方法以及Chow混合方法的单桩和群桩计算结果进行了比较,结果表明：改进计算方法是可行的,可满足精度要求。     

边坡稳定性分析条分法的细胞自动机求解和应用

边坡稳定性分析的条分法如瑞典条分法、毕肖普条分法等,现仍是边坡分析和评估的主要方法。由于刚体条分的假定,如何简化和计算条间力,是各种方法里最为关键的问题。合理的简化成为不同方法的求解途径,也是影响其计算结果差别的主要和直接的原因。然而,带着这一简化的遗憾,各种条分法仍是我们目前分析边坡稳定性的主要手段。针对研究现状,尝试更合理求解条分法条间力的求解方法,对条块刚性接触的条间力,不做简化,而引入细胞自动机求解方法,通过迭代得到条间力的解,从而使超静定的边坡条分力学模型,得到了唯一收敛的静定解。编制了相应的程序软件CASlope。结果表明,这一求解途径是可行的,所求解的结果是合理的。     

TEM中心回线法计算考虑关断时间的全区视电阻率

详细推导了TEM法中心回线装置发射电流为斜阶跃波形条件下全区视电阻率迭代反演计算方法。对正演过程中对关断时间的积分采用了解析表示。对反演过程中求导运算,通过变量变换,将对电阻率的求导变换为对时间的求导,从而也给出简单的解析关系。避开了数值积分运算和差分替代导数的近似,提高了计算精度和速度。为克服迭代反演过程中要多次计算概率积分的耗时问题,采用标准数据库的方法一劳永逸地加以解决,也大大提高了计算精度和速度。对实测资料进行了计算,并通过算例对一些问题进行了讨论。     

边坡危险滑动面及稳定安全系数的力学解析方法

为了分析边坡的稳定性,提出了一种基于极限平衡理论的分析方法。假定坡体为只受重力作用的弹性体,通过弹性力学解析解求得应力分布,将滑动面用分段二次多项式函数表示,当分段间隔足够小时,这种表示方式可以描述任意曲面形状,根据坡体滑面所满足的Mohr-Coulomb准则,定义边坡稳定安全系数,其为多项式系数的函数。通过混和罚函数优化方法可以求出最危险滑动面对应的多项式系数,获得边坡的稳定安全系数,从而判定边坡的稳定性。所提出的方法与以往方法不同,不用将滑体划分成垂直条块,不用假定滑面上的法向应力分布形式,而直接利用已有的坡体弹性应力解。算例表明:本文方法所得到的结果与简化Bishop法和有限元强度折减法获得的结果非常接近,且基本介于有限元法和简化Bishop法之间。     

层状地基中单桩负摩擦问题积分方程解法

利用Biot固结理论和积分方程方法研究了表面有堆载的层状地基中单桩负摩擦问题。根据层状饱和土的圆形载荷基本解得出了单桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化。再结合传递和刚度矩阵传递到各个层中去,对变换域内的积分方程采用Schapery 逆变换方法得到时域内单桩的近似积分方程。求解积分方程组并进行相应的数值逆变换,就可得出层状地基中的单桩在表面圆形均布载荷作用下的位移、轴力、孔压和桩侧摩阻力随时间的变化情况。计算结果表明,桩侧剪力和孔压分层明显。     

加速型改进迭代法在非饱和土渗流中的应用研究

Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。     

土坡稳定分析简化Bishop法的数值解

建立了适于边坡稳定分析简单的平面直角坐标系,得到了该坐标系下的边坡稳定分析瑞典圆弧法的积分表达式,与张天宝法进行了验证和对比;提出了基于边坡稳定分析简化Bishop法原理的积分法,推导了其安全系数Fs积分表达式,编制了数值积分计算和图形处理程序,并进行了工程应用,实践证明该方法的正确性。     

基于节点计算的自适应数值积分及其程序实现

数值积分是基于节点计算的数值分析方法的关键技术之一。针对完全无网格法的计算要求，提出了一种能自动适应任意计算域上各种节点分布方式的数值积分算法。该算法能随计算点的位置不同，自动确定积分域及积分域内的求积点，计算出数值积分结果，给出了临空边坡滑移体和圆孔方板2个实例。结果表明，按文中算法编制的程序具有良好的边界再生能力和较高的计算精度，可作为基于节点数值分析方法的一个程序模块，用于实际计算。底层算法思想，使得程序可方便地转换为其他计算语言。     

深部走向长壁开采地表移动数值模拟

在对深部开采进行界定的基础上,运用数值模拟手段,基于关键层理论建立三维深部开采模型,对走向长壁式开采地表移动进行分析,数值计算所得三维地表移动图形能够定量、直观描述地表移动,结果与概率积分法预测符合良好,实现对采煤-上覆岩土层破坏-地表沉陷全过程分析,是对煤矿开采沉陷可视化工作的重要补充。通过数值计算,总结了深部开采地表移动特点,并对比深部开采与浅部开采,分析地表影响半径及相关变形指标。分别对充分采动和非充分采动条件下概率积分法下沉系数进行计算,并与传统计算方法进行比较,结果显示:两种开采条件下,下沉系数均小于传统计算结果,基于浅部开采的经验下沉系数不再适合深部采动情况,得出下沉系数随深度增大而变小的结论,并分析这种变化产生的原因。