电偶源频率电磁测深三维地电模型有限元正演

在有源变频测深(包括可控源音频大地电磁测深CSAMT)的三维有限元正演中,应用了吸收边界条件和边元有限元算法,计算精度为3%左右,基本上可满足电偶极源三维地电模型情况下频率电磁测深响应特征的分析。同时由于广义变分概念的引入,可使泛函的变分原理应用于有耗媒质电磁场问题。      

电导率分块线性变化二维地电断面电阻率测深有限元数值模拟

在徐世浙等（１９９５）工作的基础上,给出了二维地电断面电阻率测深的变分问题以及电导率分块线性变化时的有限元方法数值解,并计算了电导率随深度变化时的地电断面的测深曲线．结果与解析解相比,精度很高．     

井地电位2.5维有限元数值模拟异常分析

从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序.然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性.在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图.结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置.     

二维地电影像——一种高分辨率岩性预测技术

在吸收传统的电阻率法勘探技术精华和方法理论的基础上,应用二维的观点和数字模拟重构技术,论述了二维地电影像系统的理论、方法、技术,并通过实测案例展示了该技术在研究工程与环境地质调查问题中的优势及可性行。     

井地电位2.5维有限元数值模拟异常分析CSCD

从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序。然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性。在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图。结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置。     

基于非结构混合网格求解二维浅水方程的一种数值方法

复杂边界下的流场数值模拟常基于非结构化网格进行求解,建立一种在非结构混合网格上求解水深平均的二维浅水方程的模型,以便精确模拟复杂边界、提高计算效率。该模型时间项离散采用隐格式使得模型具有较好的稳定性,对流项和扩散项分别采用总变差减小(Total Variation Diminishing,TVD)格式和构造辅助点的方法来离散,同时采用水深平均的标准k-ε模型来封闭湍流模型。选用两个经典验证算例检验模型,计算结果表明,基于非结构混合网格开发的模型具有较高的精度,且收敛性能较好。     

电导率分块线变化二维地电断面电阻率测深有限元？…

在徐世浙等（１９９５）工作的基础上,给出了二维地电断面电阻率测深的变分问题以及电导率分块线变化时的有限元方法数值解,并计算电导率随深度变化的地电断面的测深曲线。结果与解析解相比,精度很高。     

地下水水质模型及一种改进的特征有限元解法

本文讨论了地下水水质模型的表达形式,给出了一种适合于特征有限元法求解的表示方法,并对现有的特征有限元法进行了改进,提高了计算效率。     

二维电法数值模拟中的网格部分方法

本提出一种新的有限元及有限差分等方法中的网格剖分及单元电阻率设置算法，由此算法设计的软件，使用只需输入地电体的边界顶点坐标及其电阻率值，程序就可自动判断及设置各单元电阻率值，对程序使用来说极为简单方便。     

点源二维地电模型反演方法的研究

为了提高电法勘探数据解释的能力,寻找较好的二维反演方法,本文在研究改进的阻尼最小二乘法在二维反演中应用的同时,成功地引进和发展了Box-Kanemasu反演法。在微机上实现了电法勘探数据的二维反演解释。理论和实际地电模型的二维反演算例表明,该方法的理论是正确的,方法是切实可行的。     

三维地电断面激发极化法有限元数值模拟

用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.      

求解非稳定对流占优水质模型的非标准Galerkin有限元法

在求解非稳定地下水溶质运移模型时,若对流项占优,则模型表现出双曲方程的特性。针对这种特性,采用非标准Galerkin有限元方法进行求解是解决这类问题的有效途径。分别采用Wavelet-Galer-kin有限元方法、迎风有限元方法和特征有限元方法对强对流溶质运移模型进行了求解,并将其结果与标准Galerkin有限元和解析解进行对比。结果表明:标准Galerkin有限元方法会产生强烈的数值振荡;Wavelet-Galerkin有限元方法的时空定位效果好;迎风有限元方法能够有效降低数值振荡现象,但迎风因子对解的影响较大,而且会带来时间延迟;特征有限元方法能够提高解的精度,故可以认为特征有限元方法是求解强对流地下水溶质运移模型的首选方法。     

非线性有限单元法(续二)

(二)几何非线性问题--大应变与大位移1.司徐克林(Stricklin)等的方法在几何非线性问题中,总刚度矩阵可区分为两部分,即线性(弹性)刚度矩阵与非线性刚度矩阵     

时间序列模型和有限元模型在地下水资源评价和开发中的应用——以徐州市裂隙岩溶水为例

针对徐州市地下水长期无节制开采和管理不善而引起的多种环境地质问题,建立了与研究区裂隙岩溶含水层相适应的地下水流模型,并运用有限元方法进行求解。考虑到大气降水的多寡及其分配状况是影响地下水资源的最终决定因素,同时建立了研究区大气降水的时间序列模型,并应用于评价该区地下水资源,从而在宏观上为徐州市合理开采裂隙岩溶水资源提供科学依据。其结果对我国其它类似地区地下水资源的合理开采具有指导意义。     

求解地下水问题的三角形有限元法

从拉普拉斯方程开始,介绍了三角形有限元法求解地下水问题的基本原理,逐步推导了其数值解的公式、边界条件以及求解方法.提供了两个计算示例,并对其计算结果进行了讨论.     

Double Porosity Finite Element Method for Borehole Modeling

Summary. This paper considers the mechanical and hydraulic response around an arbitrary oriented borehole drilled in a naturally fractured formation. The formation is treated as a double porosity medium consisting of the primary rock matrix as well as the fractured systems, which are each distinctly different in porosity and permeability. The poro-mechanical formulations that couple matrix and fracture deformations as well as fluid flow aspects are presented. A double porosity and double permeability finite element solution for any directional borehole drilled in the fractured porous medium is given. Compared with conventional single-porosity analyses, the proposed double-porosity solution has a larger pore pressure in the matrix and a smaller tensile stress in the near-wellbore region. The effects of time, fracture, mud weight, and borehole inclination in the double-porosity solution are parametrically studied to develop a better understanding of physical characteristics governing borehole problems.     

岩体粘弹塑性-损伤本构模型及其有限元分析

本文将损伤力学的原理引入岩体流变学的研究,建议了一个损伤演化方程以描述岩体力学性质随时间不断弱化的现象。系统地推导了有限元计算公式。地下巷道位移的实测验算表明,本文的本构模型是合理而有效的。     

A multiresolution strategy for solving landslides using the Particle Finite Element Method

We present an approach for the simulation of landslides using the Particle Finite Element Method of the second generation. In this work, the multiphase nature (granular phase and water) of the phenomenon is considered in a staggered fashion using a single, indeformable Finite Element mesh. A fractional step and a monolithic strategy are used for the water flow and granular phase, respectively. In this way, the maximum accuracy with minimal computational times is reached. The method is completed by adding the interaction terms due to drag and pressure forces, together with a moving mesh strategy to reduce the size of the computational domain.     

关于构造带变形能的有限元法数学模拟

本文是在同一外力和边界及岩石力学性质条件下,用有限元数学模拟的方法研究不同构造变形带内部变形能的差别和分布问题。形状改变能（或称歪变能）ＵＷ比体积改变能（或称体变能）ＵＴ大将近一个数量级,ＵＷ、ＵＴ和构造变形能ＵＢ在挤压、剪切和引张变形岩带中总体依次逐步变小。ＵＷ和ＵＢ大致呈较大→大→小的变化,而ＵＴ却呈大→较大→小的趋势性变化。利用构造变形能这种规律,结合构造附加静水压力研究,进而探讨构造变形热及其影响流体运移等问题。