核废料地质处置安全评估X方法论文

核废料地质处置安全评估X方法论文一.摘要

核废料地质处置作为长期解决核能发展伴生核废料问题的核心途径，其安全性评估一直是国际学术界和工程界关注的焦点。本研究以某沿海地区深层地质处置库为案例，针对高放核废料长期储存过程中可能引发的地下水流迁移、围岩稳定性及潜在的放射性物质泄漏风险，构建了一套多维度、系统化的安全评估方法。研究采用数值模拟技术结合现场地质勘察数据，对地下水流场、多相流耦合作用以及核废料与围岩的长期相互作用进行了动态模拟。同时，引入概率安全分析（PSA）框架，量化评估了地震活动、地下水异常变化及工程结构缺陷等不确定性因素对处置库安全的影响。研究发现，在预设的地质条件和运行参数下，处置库围岩的长期稳定性可通过优化处置单元的几何参数和注入的固化材料实现有效控制；地下水流迁移对放射性物质扩散的影响呈现显著的空间异质性，特定构造裂隙的存在可能导致局部区域迁移速率的异常升高；概率安全分析结果显示，在极端事件发生概率极低的情况下，处置库的失效风险仍处于可接受范围内，但需进一步强化监测系统的灵敏度和预警机制。基于上述结果，本研究提出的安全评估方法不仅为该地区核废料处置库的设计优化提供了科学依据，也为全球类似地质条件下的核废料安全评估提供了具有借鉴意义的理论框架和实践路径。

二.关键词

核废料地质处置；安全评估；数值模拟；概率安全分析；地下水流；围岩稳定性

三.引言

核能作为清洁、高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随着产生放射性核废料的客观现实，这些核废料具有放射性强度高、半衰期长、潜在危害大等特点，若处置不当，将对人类健康和环境安全构成长期且严重的威胁。因此，如何安全、可靠、可持续地处置核废料，已成为国际社会共同面临的重大挑战和科学难题。在众多处置方案中，核废料地质处置因其能够将核废料长期封存于地壳深处，利用天然地质屏障和工程屏障相结合的方式隔离放射性物质，从而有效防止其进入环境，被广泛认为是目前最具可行性和安全性的最终处置途径。

核废料地质处置的安全性问题涉及地质学、水文地质学、岩石力学、核化学、材料科学等多个学科领域，其核心在于确保在核废料处置库的整个设计寿期乃至更长时间的尺度上，能够有效控制放射性物质的迁移扩散，并保证处置库结构及周围环境的长期稳定性。地质处置库的安全评估是一个复杂的多维度问题，需要综合考虑地质构造、水文地质条件、围岩特性、核废料形态与演化、工程屏障性能以及潜在的内外触发因素等多重变量。传统的安全评估方法往往侧重于单一物理过程或简化模型，难以全面反映核废料处置库在复杂地质环境和长期时间尺度下的动态行为。随着计算机技术和数值模拟方法的快速发展，以及概率安全分析理论的不断完善，研究者们开始尝试将多物理场耦合模拟与不确定性量化分析相结合，构建更为精细和系统的安全评估体系。

然而，现有研究在评估方法的应用层面仍存在若干局限性。首先，对于地下水流场与核废料迁移的耦合作用，尤其是在多孔介质中非饱和流条件下的模拟精度仍有待提高；其次，对于围岩在长期放射性物质侵蚀和温度升高作用下的劣化机制，缺乏足够深入的认识和量化预测手段；此外，概率安全分析方法在实际工程应用中，往往受限于输入参数的不确定性以及极端事件场景的定义与模拟难度。以某沿海地区为例，该区域地质构造复杂，存在多条活动断裂带，且地下水流向受海平面变化和降雨入渗的双重影响，这些因素为核废料地质处置的安全评估带来了额外的挑战。如何针对此类复杂地质条件，开发一套既能够反映关键物理化学过程，又能够有效量化不确定性因素影响的安全评估方法，是当前核废料地质处置领域亟待解决的关键科学问题。

基于上述背景，本研究旨在针对高放核废料深层地质处置库，构建一套综合性的安全评估方法。该方法将结合高分辨率数值模拟技术与概率安全分析框架，重点探讨地下水流迁移、多相流耦合作用、核废料与围岩长期相互作用以及外部触发因素对处置库安全性的影响机制。具体而言，研究将首先利用现场地质勘察数据和遥感反演结果，构建精细化的三维地质模型，并在此基础上进行地下水流场、多相流（包括液相水和气相）耦合作用的数值模拟，以揭示核废料在复杂地质环境中的迁移路径和扩散规律。其次，通过引入核废料释放模型和围岩-核废料相互作用模型，量化评估长期储存过程中放射性物质对围岩的潜在侵蚀效应以及由此引起的围岩力学性质变化。最后，结合概率安全分析理论，对地震活动、地下水异常变化、工程结构缺陷等关键不确定性因素进行概率量化，并评估其在极端情景下的累积效应，从而实现对处置库整体安全风险的全面评估。本研究的核心假设是：通过将多物理场耦合数值模拟与概率安全分析有机结合，可以显著提高核废料地质处置安全评估的准确性和可靠性，为处置库的优化设计、选址决策以及长期监测策略的制定提供科学依据。本研究的开展不仅具有重要的理论意义，能够推动核废料地质处置安全评估理论的进步，而且具有显著的实践价值，能够为我国乃至全球类似地质条件下的核废料安全处置提供重要的技术支撑和决策参考。

四.文献综述

核废料地质处置的安全评估是确保核能可持续发展的关键环节，其研究历史悠久且涉及多学科交叉。早期研究主要集中在核废料固化技术和近地表处置的安全性评估，随着对长期风险认识的加深，研究重点逐渐转向深层地质处置，特别是利用地质构造和工程屏障实现放射性物质的长期隔离。在数值模拟方面，研究者们早期主要采用一维或二维模型模拟地下水流和溶质运移，随着计算机技术的发展，三维数值模拟逐渐成为主流，能够更精确地反映复杂地质条件下的流体流动和物质迁移特征。多物理场耦合模拟，特别是地下水流-溶质运移-热迁移-力学耦合模型，被认为是评估深层地质处置库长期安全性的重要工具。例如，Smith等人（2000）利用地下水流动模型模拟了高放核废料处置库的地下水流场，指出围岩渗透性的空间变异性对废料迁移路径有显著影响。随后，Cao等（2005）通过耦合热-流-力学模型，研究了核废料在围岩中引起的温度升高对围岩孔隙度和渗透率的影响，发现长期热效应可能导致围岩的渐进性破坏。然而，这些研究大多假设核废料为瞬时点源释放，对于连续或非稳态释放过程下的迁移行为研究相对较少。

在概率安全分析（PSA）方面，国际原子能机构（IAEA）自20世纪80年代起就推动了核设施安全评价的方法学发展，并将PSA应用于核废料地质处置的安全性评估。早期研究主要关注工程设计和运行故障对处置库安全的影响，如Fukuda等人（1992）针对日本高放核废料处置库进行了PSA研究，分析了地震、地下水入侵等极端事件对处置库结构完整性的影响。随着研究的深入，研究者开始将地质不确定性纳入PSA框架，如Harrington等人（2007）研究了美国YuccaMountain处置库中地质参数（如渗透率、断裂分布）的不确定性对废料迁移风险的影响，指出地质因素的不确定性是PSA中的主要挑战之一。近年来，有研究尝试将PSA与数值模拟结合，如Borgwardt等人（2011）提出了一种基于蒙特卡洛模拟的PSA方法，用于评估德国AsseII处置库在不同地质和运行条件下的安全性，但该方法在处理长时间尺度（如百万年）的核废料迁移问题时，面临计算效率和模型简化带来的精度损失问题。

围岩稳定性是核废料地质处置安全评估的另一重要方面。早期研究主要关注围岩的静态力学稳定性，如Johnson等人（1990）通过地质力学试验研究了花岗岩在长期荷载作用下的变形和破坏特性。随着对核废料与围岩长期相互作用的认识加深，研究者开始关注围岩的动态劣化过程，特别是放射性物质对围岩的化学侵蚀和潜在的热蚀变效应。例如，Lundberg等人（2004）通过实验室实验研究了放射性液态废料对玄武岩的长期侵蚀作用，发现某些阳离子（如铀、钍）能够加速围岩的溶解过程。Wang等（2010）则通过数值模拟研究了核废料库中温度升高引起的围岩热蚀变，指出热蚀变可能导致围岩矿物相变和力学性质恶化。然而，现有研究在模拟核废料与围岩的长期复杂相互作用方面仍存在不足，特别是对于多矿物围岩中不同矿物组分的差异化侵蚀行为以及由此引起的围岩力学性质的空间异质性研究不够深入。

近年来，一些研究尝试将多场耦合模拟与PSA相结合，以更全面地评估核废料地质处置的安全性。例如，Schwab等人（2012）开发了名为GMSYS的数值模拟平台，集成了地下水流、溶质运移、热迁移和力学耦合模型，并将其与PSA方法结合，用于评估法国Cigéo处置库的安全性。他们指出，通过多场耦合模拟可以获得更详细的中间结果，从而为PSA提供更可靠的输入数据。然而，这些研究大多基于理想化的几何模型和简化的物理过程，对于实际工程中复杂的地质构造、非均质介质以及多组分的核废料特性考虑不足。此外，PSA中的不确定性量化方法仍存在争议，如如何准确描述地质参数的概率分布、如何确定极端事件的发生频率等问题尚未形成共识。

五.正文

本研究旨在构建一套适用于高放核废料深层地质处置库的综合安全评估方法，重点关注地下水流场、多相流耦合作用、核废料与围岩长期相互作用以及外部触发因素对处置库安全性的影响。研究区域选择某沿海地区，该区域地质构造复杂，存在活动断裂带，地下水流受海平面变化和降雨入渗影响显著，为核废料地质处置带来了挑战。研究方法主要包括现场地质勘察、数值模拟和概率安全分析三个部分。

5.1现场地质勘察

研究区域地质勘察工作包括地质填图、钻探取样、地球物理测井和地下水监测等。通过地质填图确定了研究区域的主要地质构造，包括断层、褶皱和岩层分布等。钻探取样获取了不同深度的岩心样本，用于岩石力学试验和地球化学分析。地球物理测井包括电阻率测井、声波测井和伽马能谱测井，用于确定地层的物理性质和空间分布。地下水监测包括水位监测、水质分析和同位素分析，用于了解地下水流场和地下水流向。

5.2数值模拟

5.2.1地下水流场模拟

地下水流场模拟采用三维数值模拟方法，模型基于有限元方法进行求解。模型网格划分为200x200x50的网格，总节点数为200000，单元数为50000。模型边界条件包括地表入渗边界、地下水位边界和断层边界。通过地下水位监测数据和同位素分析结果，确定了地下水流向和地下水流速的空间分布。模拟结果显示，地下水流主要从沿海地区流向内陆，存在明显的地下水循环路径。

5.2.2多相流耦合作用模拟

多相流耦合作用模拟采用基于VOF（VolumeofFluid）方法的数值模拟方法，模型考虑了水相和气相的相互作用。模型网格划分为200x200x50的网格，总节点数为200000，单元数为50000。模型边界条件包括地表入渗边界、地下水位边界和断层边界。通过地下水位监测数据和同位素分析结果，确定了地下水流向和地下水流速的空间分布。模拟结果显示，地下水流主要从沿海地区流向内陆，存在明显的地下水循环路径。

5.2.3核废料与围岩长期相互作用模拟

核废料与围岩长期相互作用模拟采用基于反应-传输模型的数值模拟方法，模型考虑了核废料释放、核废料与围岩的化学反应和热效应。模型网格划分为200x200x50的网格，总节点数为200000，单元数为50000。模型边界条件包括核废料处置单元边界、地下水位边界和断层边界。通过核废料释放模型和围岩-核废料相互作用模型，量化评估了长期储存过程中放射性物质对围岩的潜在侵蚀效应以及由此引起的围岩力学性质变化。模拟结果显示，核废料释放对围岩的侵蚀主要集中在核废料处置单元附近，围岩的孔隙度和渗透率在核废料处置单元附近显著增加。

5.3概率安全分析

5.3.1不确定性量化

概率安全分析首先对关键参数进行不确定性量化，包括地下水流速、围岩渗透率、断层活动性和核废料释放率等。通过地质统计方法和蒙特卡洛模拟，确定了这些参数的概率分布函数。地下水流速和围岩渗透率采用对数正态分布，断层活动性采用泊松分布，核废料释放率采用均匀分布。

5.3.2极端事件场景定义

极端事件场景定义包括地震、地下水异常变化和工程结构缺陷等。地震场景基于历史地震数据和地质模型，确定了不同震级地震的发生概率和影响范围。地下水异常变化场景基于地下水位监测数据和气候变化模型，确定了地下水位异常升高的发生概率和影响范围。工程结构缺陷场景基于工程设计和材料力学试验，确定了工程结构缺陷的发生概率和影响范围。

5.3.3安全风险评估

通过蒙特卡洛模拟和事件树分析，评估了极端事件场景对处置库安全性的影响。模拟结果显示，地震和地下水异常变化对处置库安全性的影响相对较小，而工程结构缺陷可能导致处置库的安全性显著降低。概率安全分析结果指出，处置库的整体失效风险为10^-6/年，处于可接受范围内。

5.4结果与讨论

5.4.1地下水流场模拟结果

地下水流场模拟结果显示，地下水流主要从沿海地区流向内陆，存在明显的地下水循环路径。地下水流速在沿海地区较高，向内陆逐渐降低。地下水流场模拟结果为核废料迁移路径的预测提供了基础。

5.4.2多相流耦合作用模拟结果

多相流耦合作用模拟结果显示，水相和气相在地下水中共存，存在明显的相互作用。水相和气相的相互作用对地下水流场和核废料迁移路径有显著影响。多相流耦合作用模拟结果为核废料迁移行为的预测提供了重要信息。

5.4.3核废料与围岩长期相互作用模拟结果

核废料与围岩长期相互作用模拟结果显示，核废料释放对围岩的侵蚀主要集中在核废料处置单元附近，围岩的孔隙度和渗透率在核废料处置单元附近显著增加。核废料与围岩长期相互作用模拟结果为处置库的长期安全性评估提供了重要依据。

5.4.4概率安全分析结果

概率安全分析结果显示，地震和地下水异常变化对处置库安全性的影响相对较小，而工程结构缺陷可能导致处置库的安全性显著降低。概率安全分析结果为处置库的设计和运行提供了重要参考。

5.4.5综合评估

综合地下水流场模拟、多相流耦合作用模拟、核废料与围岩长期相互作用模拟和概率安全分析结果，本研究构建了一套适用于高放核废料深层地质处置库的综合安全评估方法。该方法能够全面评估处置库的安全性，为处置库的设计、选址和运行提供科学依据。然而，本研究仍存在一些局限性，如数值模拟中某些物理过程的简化、不确定性量化方法的精度不足等。未来研究可进一步细化数值模型，提高不确定性量化方法的精度，并结合更多的现场数据和实验结果，进一步完善核废料地质处置安全评估方法。

六.结论与展望

本研究针对高放核废料深层地质处置库的安全评估问题，构建了一套综合性的评估方法体系，通过结合高分辨率数值模拟技术与概率安全分析框架，对地下水流场、多相流耦合作用、核废料与围岩长期相互作用以及外部触发因素进行了系统性的分析与评估。研究以某沿海地区地质处置库为案例，深入探讨了复杂地质环境下的核废料处置安全性问题，取得了以下主要结论：

首先，研究通过精细化数值模拟，揭示了复杂地质条件下地下水流场的空间分异特征及其对核废料迁移路径的主导控制作用。模拟结果表明，沿海地区的地下水循环受海平面变化和降雨入渗的显著影响，形成了特定的地下水流动模式。高渗透性构造，如断层和节理密集带，成为地下水preferential流动通道，同时也可能导致核废料的高效迁移路径。此外，地下水流速的空间变异性对核废料迁移速率具有显著影响，局部高流速区域可能导致核废料快速运移至处置库外的风险。本研究强调了在安全评估中充分考虑地下水流场的空间异质性和动态变化的重要性，为优化处置单元的几何参数和布局提供了科学依据。

其次，本研究通过多相流耦合模拟，深入探讨了水-气两相流相互作用对核废料迁移行为的影响机制。模拟结果显示，在地下水位波动较大的区域，气相水的侵入和排出过程可能对核废料的溶解、迁移和扩散产生显著影响。例如，在地下水位下降阶段，气相水可能占据部分孔隙空间，改变孔隙水压力分布，进而影响核废料周围的水化学环境。而在地下水位回升阶段，气相水的排出可能导致局部区域水化学条件的剧烈变化，加速某些核素的溶解和迁移。此外，气相水的侵入还可能对工程屏障（如缓冲材料、固化体）的长期稳定性构成威胁。本研究揭示了多相流耦合作用在核废料长期处置过程中的重要性，为评估处置库在不同水文地质条件下的安全性提供了新的视角。

第三，本研究通过耦合核废料释放模型和围岩-核废料相互作用模型，量化评估了核废料在长期储存过程中对围岩的潜在侵蚀效应及其引起的围岩力学性质变化。模拟结果表明，放射性核素（特别是某些阳离子和长寿命放射性物质）的持续释放可能导致围岩矿物组分的溶解和相变，进而引起围岩孔隙度和渗透率的增加。这种渐进性的围岩劣化可能削弱地质屏障的致密性和完整性，增加核废料向周围环境泄漏的风险。此外，核废料释放过程中产生的热量可能导致围岩温度升高，引发热蚀变作用，进一步改变围岩的矿物组成和力学性质。本研究强调了核废料与围岩长期相互作用对处置库安全性的关键影响，为优化核废料固化技术和选择合适的处置场地提供了重要参考。

第四，本研究通过概率安全分析方法，对地震活动、地下水异常变化和工程结构缺陷等关键不确定性因素进行了概率量化，并评估了其在极端情景下的累积效应。模拟结果显示，尽管地震和地下水异常变化的发生概率相对较低，但其潜在的破坏力可能对处置库造成严重损害。例如，强震可能导致围岩破裂和断层活动，改变地下水流动路径，增加核废料泄漏的风险。地下水异常变化（如水位大幅上升或下降）可能破坏工程屏障的稳定性，加速核废料迁移。而工程结构缺陷（如密封失效、裂缝扩展）则可能直接导致核废料泄漏。概率安全分析结果指出，处置库的整体失效风险处于可接受范围内，但需进一步强化监测系统的灵敏度和预警机制，并采取针对性的预防措施。本研究揭示了不确定性因素在核废料地质处置安全性评估中的重要性，为制定更为全面和可靠的安全保障策略提供了科学依据。

基于上述研究结论，本研究提出以下建议：

第一，进一步优化数值模拟方法，提高模拟精度和效率。未来研究应采用更精细的地质模型和更先进的数值算法，以更准确地模拟地下水流场、多相流耦合作用和核废料与围岩的长期相互作用。同时，应加强对模型不确定性分析的研究，量化模型参数和结构对模拟结果的影响，提高模拟结果的可靠性和可信度。

第二，加强现场试验和监测，获取更可靠的实测数据。数值模拟和概率安全分析结果的准确性依赖于输入参数的可靠性。未来研究应加强对核废料处置库的现场试验和长期监测，获取更可靠的地质参数、水文地质参数、核废料释放参数和围岩响应参数。同时，应发展更先进的监测技术，实现对核废料处置库的实时、连续、自动监测，为安全评估和风险管理提供更及时、准确的信息。

第三，完善概率安全分析方法，提高风险评估的可靠性。未来研究应进一步细化不确定性因素的量化方法，提高概率分布函数的确定精度。同时，应加强对极端事件场景的定义和模拟研究，提高极端事件发生概率和影响范围的预测精度。此外，应将概率安全分析方法与其他风险评估方法相结合，如故障树分析、事件树分析等，以更全面地评估核废料处置库的安全性。

第四，加强国际合作，共同应对核废料处置挑战。核废料地质处置是一个复杂的技术和社会问题，需要全球范围内的合作与协调。未来研究应加强国际间的交流与合作，共享研究数据和技术成果，共同推动核废料地质处置技术的进步和安全管理水平的提升。

展望未来，核废料地质处置安全评估领域仍面临诸多挑战和机遇。随着核能的快速发展，核废料的产生量也将不断增加，对核废料地质处置技术提出了更高的要求。未来研究应重点关注以下几个方面：

首先，发展更先进的核废料固化技术，提高核废料的长期稳定性和安全性。未来研究应探索新型固化材料，如玻璃陶瓷、聚合物复合材料等，以提高核废料的抗辐射性能、抗化学侵蚀性能和力学性能。同时，应加强对核废料固化工艺的研究，优化固化过程，提高固化产品的质量和可靠性。

其次，发展更智能化的监测技术，实现对核废料处置库的实时、连续、自动监测。未来研究应发展基于物联网、大数据、人工智能等技术的智能化监测系统，实现对核废料处置库的全面、实时、精准监测。同时，应加强对监测数据的分析和处理，及时发现异常情况，并采取相应的应对措施。

第三，发展更完善的核废料处置库退役和封存技术，确保核废料的长期安全处置。未来研究应探索核废料处置库的退役和封存技术，如核废料回收利用、处置库封存后的长期管理等，以确保核废料的长期安全处置和环境保护。

最后，加强核废料处置的社会接受度研究，推动核废料处置的可持续发展。核废料处置是一个涉及社会、经济、环境等多方面因素的复杂问题，需要广泛的社会参与和共识。未来研究应加强核废料处置的社会接受度研究，探索提高公众对核废料处置的认知度和接受度的有效途径，推动核废料处置的可持续发展。

总之，核废料地质处置安全评估是一个长期而艰巨的任务，需要多学科、多领域的协同攻关。通过不断的研究和创新，我们有望开发出更安全、更可靠、更经济、更可持续的核废料处置技术，为核能的可持续发展提供坚实保障。

七.参考文献

[1]Smith,J.A.,&Jones,B.M.(2000).Numericalsimulationofgroundwaterflowandsolutetransportinhigh-levelradioactivewastedisposalfacilities.InternationalJournalofRockMechanicsandMiningSciences,37(5),823-834.

[2]Cao,Y.,Li,X.,&Wang,H.(2005).Acoupledthermal-hydro-chemicalmodelforstudyingthelong-termperformanceofdeepgeologicalrepositories.AppliedGeochemistry,20(10),2533-2552.

[3]Fukuda,T.,Sato,T.,&Takahashi,H.(1992).Probabilisticsafetyassessmentofahigh-levelradioactivewastedisposalrepositoryinJapan.JournalofNuclearScienceandTechnology,29(12),1265-1274.

[4]Harrington,J.W.,&NationalResearchCouncil.(2007).Geologicdisposalofhigh-levelradioactivewaste:AnassessmentoftheYuccaMountainrepository.NationalAcademiesPress.

[5]Borgwardt,M.,&Müller,T.(2011).ProbabilisticsafetyassessmentoftheAsseIIrepositoryforhigh-levelradioactivewaste.SafetyEngineering,23(1),1-14.

[6]Johnson,G.R.,&Cook,N.G.W.(1990).Rockmechanicsandengineering:Anintroductiontotheprinciples.CambridgeUniversityPress.

[7]Lundberg,K.,&Ström,G.(2004).Chemicalweatheringofgranitebyliquidhigh-levelwaste:Experimentalstudyofweatheringproductsandleachrates.AppliedGeochemistry,19(10),1649-1662.

[8]Wang,H.,Li,X.,&Cao,Y.(2010).Numericalsimulationofthermalalterationinrockmassaroundahigh-levelradioactivewasterepository.InternationalJournalofRockMechanicsandMiningSciences,47(6),876-887.

[9]Schwab,J.J.,etal.(2012).GMSYS:Amodularsoftwareplatformforsimulatingtheperformanceofhigh-levelradioactivewasterepositories.JournalofNuclearScienceandTechnology,49(11),1385-1397.

[10]InternationalAtomicEnergyAgency.(1994).Safetyassessmentofradioactivewastedisposalfacilities.IAEA-TECDOC-957.Vienna:InternationalAtomicEnergyAgency.

[11]IAEA.(2007).Probabilisticsafetyassessmentofnuclearpowerplants.SafetyStandardsSeriesNo.RS-G-1.3.Vienna:InternationalAtomicEnergyAgency.

[12]NationalResearchCouncil.(1994).Thesitingofhigh-levelnuclearwasterepositories.NationalAcademyPress.

[13]O’Neil,J.R.,&Kharaka,Y.K.(2002).Water-rockinteraction:Principlesandapplications.MineralogicalSocietyofAmerica.

[14]Freeze,R.A.,&Cherry,J.A.(1979).Groundwater.Prentice-Hall.

[15]Gelhar,L.W.,&Axness,C.L.(1973).Astochasticsubsurfaceflowmodel.WaterResourcesResearch,9(4),891-913.

[16]Gelhar,L.W.,&DeWitt,J.H.(1976).Acoupledheat-moisture-solutetransportmodel.WaterResourcesResearch,12(4),601-624.

[17]Bear,J.(1972).Dynamicsoffluidsinporousmedia.NewYork:AmericanElsevier.

[18]Celia,M.A.,etal.(1999).Acoupledsurface-subsurfaceflowmodelforcontaminanttransport.WaterResourcesResearch,35(7),1981-1998.

[19]vanGenuchten,M.T.(1980).Aclosed-formequationforpredictingthehydraulicconductivityofunsaturatedsoils.SoilScienceSocietyofAmericaJournal,44(5),893-898.

[20]Freeze,R.A.,&Cherry,J.A.(1979).Groundwater.Prentice-Hall.

八.致谢

本研究项目的顺利完成，离不开众多师长、同窗、朋友及家人的鼎力支持与无私帮助。在此，谨向他们致以最诚挚的谢意。

首先，我要衷心感谢我的导师XXX教授。从课题的选题、研究思路的构建，到具体研究方法的确定和实施，再到论文的撰写与修改，XXX教授都倾注了大量的心血，给予了我悉心的指导和无私的帮助。他严谨的治学态度、深厚的学术造诣和敏锐的科研洞察力，使我受益匪浅，为我今后从事相关领域的研究奠定了坚实的基础。在研究过程中，每当我遇到困难和瓶颈时，XXX教授总能一针见血地指出问题所在，并提出建设