平面单裂隙介质中核素的迁移模型

本文提出核素在平面单裂隙中的迁移模型 ,用于模拟放射性核素在裂隙介质中的迁移问题 ,并以实例说明该模型的可靠性     

双重介质中核素迁移模型的解析解及其应用

给出了具有单裂隙的双重介质中，裂隙一端有分时注入源情形下的迁移模型的解析解及一些简化条件下的解析解，并利用解析解研究了核素在双重介质中的迁移问题。     

双重介质中核素迁移模型的近似解析解及实验室使用条件

本文将双重介质中核过移模型的解析解进行了简化，得到了便于应用的近似解析解及其相关结论经推导给出了实验室利用近似解析解的条件以及求阻滞系数、固有扩散系数的公式。     

似圆状双重介质模型的上游迎风有限元法及其应用

针对似圆状双重介质模型给出了上游迎风有限元法，克服了一般有限单元法所导致的数值不稳定性，并以济南岩溶区为例证明了这一方法的可行性和有效性。     

双重介质中的一类核素迁移数学模型及其反演

探讨了孔隙与单裂隙双重介质中的一类核素迁移数学模型及其反演问题。该核素迁移模型是一个耦合的抛物型方程组定解问题。若已知排污点的核素浓度变化规律,利用Laplace变换及其逆变换方法,求得了核素迁移模型正问题的解析解;反之,由下游裂隙中某个点的实测核素浓度,利用偏微分方程的叠加原理和反问题的拟解法,反求出核素迁移模型反问题的解,即排污点的核素状态。最后,给出核素迁移模型的正问题和反问题的数值模拟。数值结果表明,正问题的解析解能够刻画核素的迁移规律,也显示出所提反问题方法能有效地反演核素污染源。     

利用算子分裂迎风均衡格式解对流为主溶质运移问题

水污染模拟问题是水流问题与溶质运移问题的耦合问题.各种常见的数值解法在以对流为主溶质运移问题的求解中都会遇到困难,如用有限单元法或有限差分法时,会产生数值弥散与过量这两类误差.引入算子分裂迎风均衡格式法求解对流为主的水污染模拟问题,较好地克服了数值弥散和数值解出现振荡问题,该格式具有良好的稳定性、单调性及守恒性特点.     

放射性核素在裂隙介质中迁移模型研究综述

指出了研究核素迁移问题的意义,归纳了裂隙介质的几种简化形式,概括总结了国内外对放射性核素在裂隙介质中迁移模式研究中的进展.提出了迁移数学模式的分类方法,并指出了目前迁移问题研究中存在的不足.     

似圆状双重介质模型的上游迎风有限元法及其应用

针对似圆状双重介质模型给出了上游迎风有限元法，克服了一般有限单元法所导致的数值不稳定性，并以济南岩溶区为例证明了这一方法的可行性和有效性。     

求解对流占优地下水水质模型的有限元高精度广义迎风单元均衡格式及应用

根据水质模型的具体特点，对不同的方程采用不同方法，水流问题用有限元法；对流弥散方程先用算子分裂的方法分解为两个方程，即对流方程和弥散方程，前者用高精度广义迎风格式求解，对弥散方程则采用多单元均衡格式法求解，最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟，可以看出在求解对流弥散定解问题时，广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题，与有限元相比有较大改进。     

页岩介质中放射性核素迁移的实验研究及其数学模型

在自行设计的实验装置中进行了放射性核素13 1I在页岩介质中的扩散、渗透 -弥散的实验。采用渗透法及扩散法获得了实验曲线并计算出弥散系数、阻滞系数、孔隙度等参数。实验表明13 1I在页岩介质中为弱吸附核素 ,其阻滞系数较小。在此实验基础上建立了放射性核素13 1I在页岩介质中的一维迁移物理模型及数学模型。结果表明13 1I随地下水迁移得比较慢 ,因此页岩对13 1I的迁移能起到一定的阻滞作用     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

三维溶质运移问题的分步广义迎风解法

对对流占优的三维溶质运移问题提出了分步广义迎风解法，首先利用Ｎ．Ｎ．Ｙａ－ｎｅｎｋｏ对水动力弥散方程分步求解的思想，将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题，对对流定解问题采用广义迎风对偶单元均衡法求解，对扩散定解问题采用一般的Ｇａｌｅｒｋｉｎ有限元法求解，不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题，而且避免了求节点速度这一步，节省运算步骤，对井点的浓度变化给出了更合适的求解方法。     

A TIME-STEPPING FINITE ELEMENT METHOD FOR ANALYSIS OF CONTAMINANT TRANSPORT IN FRACTURED POROUS MEDIA

This paper describes the development of a finite element method for analysing contaminant transport in double-porosity geomaterials using a time-stepping approach. In many cases, double-porosity models may be used to represent fractured rock formations and fissured soils. A distinctive feature of utilizing this kind of model is that it is not necessary to have an intimate knowledge of the nature, distribution and properties of individual fractures and fracture arrangement since the fracture geometry and details are considered only in an averaged or equivalent continuum sense. The flux exchange that occurs between the fluid in the fractures and in the solid matrix is represented by a linear heriditary process. This has the consequence that in order to carry the solution forward from time t to t+Δt, it is necessary to know and to store the complete contaminant history up to time t. This paper shows that all the hereditary information necessary to carry the solution forward is contained in the values of certain hereditary variables at time t so that it is not necessary to store the complete time history and consequently a more efficient numerical process can be developed.     

溶蚀岩体随机结构模型建立及其在岩体渗漏评价中的应用

针对中小规模裂隙和溶蚀洞穴造成的溶蚀岩体渗漏 ,本文提出了溶蚀岩体渗漏评价的一种新思路 :首先对现场开挖所揭露的有限的裂隙和溶蚀洞穴观测资料进行统计分析 ,利用统计分析结果建立溶蚀岩体随机结构模型 ;然后根据随机结构模型所反映出的裂隙和溶蚀洞穴有关参数计算渗透参数 ;最后结合区域水力学条件计算岩体渗漏量。多次模拟后可以得到渗漏量的分布规律。其中 ,建立合理的随机结构模型是关系本方法成败的关键 ,文中给出了建立溶蚀岩体随机结构模型的步骤 ,提供了随机模拟的编程思路和程序框图 ,并结合具体的工程实例对该方法的应用前景进行了讨论     

岩石复电阻率Dias模型及其参数求取方法

Dias模型是众多描述岩石复电阻率的频散特性模型中的一种,它不但参数少,且可由这些参数导出具有特定物理意义的参数。针对Dias模型参数的求取及反演解的不确定性评价,这里提出了一种模拟退火阻尼最小二乘联合反演方法。该方法首先将参数求解空间离散化,利用模拟退火找到全局最优解区间,然后利用阻尼最小二乘法在解空间进行求解。该方法有效地克服了原来两种方法的反演速度慢、反演结果依赖于初值的缺点,并实现了对反演结果多解性的评价。研究结果表明,对于岩石复电阻率数据的处理解释,Dias模型的使用比传统的Cole模型及其变形更具有潜在的优势。     

核废料贮库围岩介质THM耦合的定解问题及其加权积分方程

以核废料贮库裂隙岩体介质热-液-力耗散过程的定力解方程为基础，结合核废料地下贮存，分析了热、液、力三方面的边值及初始条件。根据计算力学加权残值法，导出了定解问题的加权积分方程，为实现核废料贮库围岩介质THM耦合有限元数值计算，打下了理论基础。     

烃类垂向微渗漏过程的积木块模型及差分格式

这里提出以代表质量守恒的反应对流扩散方程作为主控方程的烃类垂向微渗漏方程组的差分格式,即双向一维分裂校正差分格式,并建立地层积木块模型对该格式的边界进行讨论。差分格式是预估～校正差分格式的一种改进形式,它融合了Crank-Nicolson格式、交替方向隐格式、预估～校正差分格式的特点,具有二阶差分精度,且无条件稳定。由于差分格式将每一步都归结为求解三对角线方程组,因此适合并行运算。数值实验表明,应用差分格式的数值模拟结果符合烃类垂向微渗漏过程的理论模型,可作为烃类垂向微渗漏过程分析的计算方法。     

VTI介质转换波速度分析方法研究及应用

速度分析是转换波处理的核心问题之一,采用单平方根方程对具有垂直对称轴的横向各向同性(VTI)介质进行速度分析,使得转换波的处理不依赖于纵波的数据,既简化了转换波处理流程,又消除了由于纵波处理不精确所带来的传递误差。通过理论分析和数值计算,总结并设计了一套求取VTI介质转换波动校正参数的方法和处理流程,总结出了采用双参数扫描法获得较高分辨率和精度所需满足的必要条件。通过批处理和交互结合的方式对实际资料进行处理,取得了较好的处理效果,证明了该方法及软件的正确性。