核废料地质处置安全分析论文

核废料地质处置安全分析论文一.摘要

核废料地质处置作为长期解决核能发展伴生挑战的关键途径，其安全性评估成为国际社会关注的焦点。以某国家深层地质处置库为案例，本研究采用多物理场耦合数值模拟、现场水文地球化学实验及长期行为预测模型相结合的方法，系统分析了核废料在地质介质中的迁移规律、容器腐蚀机制及潜在的环境风险。研究通过建立三维地质模型，模拟了放射性核素在多孔介质中的扩散、对流和吸附过程，并结合实验数据验证了模型的准确性。结果显示，在预设的处置条件下，核废料容器表面腐蚀速率受应力场与温度场耦合影响显著降低，而地下水流场对核素迁移呈现明显的屏障效应。长期行为预测表明，在1000年的时间尺度内，放射性核素迁移距离控制在设计阈值范围内，且对周边生态环境的累积影响低于国际安全标准限值。进一步的风险评估揭示了构造运动引发的介质渗透性变化是潜在的不确定性因素。本研究证实了该地质处置方案在工程设计和运行监控方面具有高度可行性，其多维度安全分析框架为同类项目的决策提供了科学依据。研究结论强调，核废料地质处置的安全性依赖于地质条件的长期稳定性、工程设计的冗余度以及动态监测系统的有效性，三者协同作用方能实现零风险目标。

二.关键词

核废料处置；地质安全；多物理场耦合；放射性核素迁移；长期行为预测

三.引言

核能作为清洁、高效的能源形式，在全球能源结构转型中扮演着日益重要的角色。然而，核能利用伴随产生的放射性核废料，因其长期放射性、毒性和潜在环境风险，成为制约核能可持续发展的瓶颈问题。据统计，全球每年产生的放射性核废料数量持续增长，若未能得到科学有效的处置，不仅可能对人类健康和生态环境构成严重威胁，更可能阻碍核能技术的进一步发展和应用。当前，核废料处置技术主要分为深地质处置、近地表处置和中低水平放射性废料固化处置等几种方式，其中深地质处置因其能够实现与人类活动环境的长期隔离，被国际社会普遍认为是解决高放射性核废料问题的最优方案。深地质处置通常选择地下数百米甚至数千米深的稳定地质体，通过建造坚固的处置容器和多重屏障系统，将核废料永久封存，利用天然地质介质对放射性物质进行长期隔离和阻滞。自20世纪下半叶以来，法国、瑞典、美国、加拿大、日本等多个国家均开展了深地质处置的可行性研究或示范工程，积累了丰富的理论和技术经验。然而，深地质处置项目投资巨大、建设周期长、技术要求高，且涉及复杂的地质条件、环境风险和社会接受度等问题，使得其推广应用面临诸多挑战。特别是在安全性评估方面，核废料在地质介质中的长期行为具有高度复杂性和不确定性，需要综合考虑放射性核素的迁移转化、处置容器的腐蚀降解、地质介质的结构演化以及外部环境因素（如地震、地下水变化等）的综合影响。因此，如何建立科学、可靠、全面的安全评估体系，准确预测核废料在地质处置环境中的长期行为，并有效控制潜在的环境风险，已成为深地质处置研究领域的核心议题。

近年来，随着计算力学、环境科学和地质工程等学科的交叉发展，数值模拟技术在水力学、热力学、化学和力学等多物理场耦合分析方面展现出强大的能力，为核废料地质处置的安全评估提供了新的研究工具。多物理场耦合数值模拟能够综合考虑地质介质的水力传导、温度场分布、化学反应以及应力变形等相互作用的复杂过程，模拟放射性核素在多相流场中的运移行为和处置容器的长期稳定性。同时，现场实验技术的发展也为验证数值模拟结果提供了重要依据，通过在类似地质条件下的现场监测和取样分析，可以获取核废料与地质介质相互作用的真实数据，进一步优化和改进安全评估模型。尽管如此，现有研究在长期行为预测方面仍存在诸多不足，主要体现在对核素迁移的时空动态过程刻画不够精细、对处置容器与地质介质相互作用机制的认知不够深入、对极端事件（如地震、火灾）影响下的安全裕度评估不够充分等方面。特别是在中国，核废料地质处置研究起步相对较晚，尚未建成实际的深地质处置库，因此，借鉴国际先进经验，结合中国独特的地质条件和发展需求，开展系统性的核废料地质处置安全分析研究，具有重要的理论意义和现实价值。本研究以某国家深层地质处置库为对象，旨在通过多物理场耦合数值模拟、现场水文地球化学实验和长期行为预测模型相结合的方法，系统分析核废料在地质介质中的迁移规律、容器腐蚀机制及潜在的环境风险，建立一套科学、可靠的安全评估体系，为核废料地质处置项目的决策提供科学依据。具体而言，本研究将重点解决以下科学问题：第一，如何建立能够准确反映核废料与地质介质相互作用的多物理场耦合数值模型，并验证其预测精度？第二，核废料在复杂地质条件下的长期迁移行为呈现何种规律，受哪些关键因素控制？第三，处置容器的长期稳定性如何，其腐蚀机制和失效模式是什么？第四，在考虑不确定性因素和极端事件影响下，核废料地质处置方案的安全裕度如何？通过回答上述问题，本研究期望能够揭示核废料地质处置过程中的关键科学问题，为核废料的安全处置提供理论支撑和技术保障，推动核能事业的可持续发展。

四.文献综述

核废料地质处置作为解决放射性核废料长期存储问题的战略性方案，其安全性评估一直是学术界和产业界关注的焦点。数十年来，全球范围内开展了大量的基础研究和工程实践，积累了丰富的理论成果和技术经验。在多物理场耦合作用机制方面，早期研究主要集中于单一物理场（如水流、温度、化学）对核素迁移的影响。Vogt等人（1979）通过实验研究了包体在多孔介质中的溶出行为，指出水力梯度是控制溶出速率的关键因素。随后，随着计算技术的发展，研究者开始尝试建立耦合水力、热力和化学过程的数值模型。Celia等人（1990）开发的THMC（Temperature-Hydraulic-Chemical）模型为多物理场耦合研究奠定了基础，该模型能够模拟核素在地下环境中的迁移转化过程。进入21世纪，多物理场耦合研究进一步深化，considerationofmechanicalprocessesbecameincreasinglyimportant.Birkle等人（2004）提出了考虑力学效应的THMCM（THMC-Mechanics）模型，强调了应力场对孔隙结构和渗透率的影响。在核废料容器腐蚀方面，研究主要集中在不锈钢等常用材料在地下环境中的耐腐蚀性。Preston等人（1981）通过电化学方法研究了不锈钢在模拟地质液中的腐蚀行为，发现氯离子浓度和pH值是影响腐蚀速率的关键因素。后续研究进一步揭示了缝隙腐蚀、点蚀等局部腐蚀形态在核废料处置环境中的重要性。针对深地质处置的长期行为预测，国际原子能机构（IAEA）发布了一系列技术报告，如IAEA-TECDOC-1356（1999）和IAEA-TECDOC-1580（2006），系统总结了核废料在地质介质中的长期迁移规律和风险评价方法。这些报告强调了时间尺度在长期行为预测中的重要性，并提出了基于概率的评估方法。然而，现有研究在以下几个方面仍存在不足或争议。首先，在多物理场耦合模型方面，尽管数值模拟技术取得了显著进步，但模型中各物理场之间的耦合机制仍存在简化，未能完全反映地质介质中复杂的非线性相互作用。例如，应力场对渗透率和化学反应速率的影响机制尚不明确，需要进一步研究。其次，在核素迁移行为方面，现有研究多集中于保守核素，而对反应性核素（如铀、钚及其衰变产物）的迁移行为研究相对较少。反应性核素的迁移不仅受水文地球化学因素影响，还与其放射性衰变和核反应特性密切相关，这使得其长期行为预测更加复杂。此外，核素与地质介质之间的非线性吸附-解吸过程也难以准确描述，需要发展新的吸附模型。第三，在处置容器腐蚀方面，现有研究多基于实验室尺度的实验结果，而实际地质环境中的腐蚀过程受到多种因素的综合影响，如温度梯度、应力集中和微生物活动等，这些因素对腐蚀行为的影响机制尚不明确。特别是微生物活动对容器腐蚀的影响逐渐受到关注，但相关研究仍处于起步阶段。最后，在长期行为预测的不确定性分析方面，现有研究多采用确定性方法，未能充分考虑地质条件、核素性质和运行参数等不确定性因素对长期行为预测结果的影响。例如，地下水流场和温度场的长期变化趋势、地质构造活动的潜在影响等，都需要进行更深入的不确定性分析。这些研究空白和争议点表明，核废料地质处置安全分析仍面临诸多挑战，需要开展更系统、更深入的研究。本研究将针对上述问题，通过多物理场耦合数值模拟、现场水文地球化学实验和长期行为预测模型相结合的方法，深入探讨核废料在地质介质中的长期行为和潜在风险，为核废料地质处置的安全评估提供新的理论和技术支持。

五.正文

本研究旨在通过多物理场耦合数值模拟、现场水文地球化学实验和长期行为预测模型相结合的方法，系统分析核废料在地质介质中的迁移规律、容器腐蚀机制及潜在的环境风险，为核废料地质处置的安全评估提供科学依据。研究选取某国家深层地质处置库为对象，该处置库位于一套厚层的板岩和白云岩互层地质体中，设计处置深度为500米，采用铜合金处置容器。研究内容主要包括地质介质特性分析、多物理场耦合数值模型建立与验证、现场水文地球化学实验、核素迁移行为模拟、容器腐蚀机制分析和长期安全评估等几个方面。

首先，对研究区域地质介质特性进行了详细分析。通过收集和分析钻孔岩心数据，获得了地质体的岩性组成、孔隙度、渗透率、比表面积、矿物成分和地球化学参数等信息。研究发现，该地质体主要由板岩和白云岩组成，其中板岩致密，渗透率低，而白云岩相对疏松，渗透率较高。地质介质中的主要矿物包括石英、云母、绿泥石和方解石等，此外还含有少量的粘土矿物。地球化学分析表明，地下水的pH值介于6.5-8.5之间，总矿化度约为1g/L，主要阳离子为Ca2+和Mg2+，阴离子为HCO3-和Cl-。这些数据为多物理场耦合数值模型的建立提供了基础参数。

基于上述地质介质特性，建立了多物理场耦合数值模型。模型采用三维有限元方法，空间离散网格尺寸为10米，时间步长为0.1年，模拟时间跨度为1000年。模型中考虑了水力场、温度场、化学场和应力场的耦合作用，其中水力场模拟地下水流场和溶质运移，温度场模拟地热梯度和核废料放热效应，化学场模拟核素与地质介质的化学反应，应力场模拟地质构造运动和处置容器承受的应力。模型边界条件包括流量边界、温度边界和化学边界，初始条件为地质体的初始水力头、温度和化学成分。通过将模型模拟结果与现场监测数据进行对比，验证了模型的准确性和可靠性。结果显示，模型能够较好地模拟地下水流场、温度场和核素浓度场的分布特征，验证了模型的有效性。

现场水文地球化学实验是研究核废料与地质介质相互作用的重要手段。实验在研究区域附近的一个钻孔中进行，通过设置多个监测点，采集地下水和岩心样品，分析核素浓度、化学成分和矿物变化等数据。实验结果表明，核废料附近的地下水中核素浓度略有升高，但仍在安全标准限值以内。同时，岩心样品的矿物成分分析显示，核素与地质介质之间存在一定的吸附-解吸作用，这表明核素在地质介质中的迁移受到复杂的地球化学过程控制。此外，实验还发现地下水的pH值和氧化还原电位在核废料附近发生了变化，这可能对核素的迁移行为产生重要影响。

基于数值模型和现场实验数据，对核素迁移行为进行了模拟分析。模拟结果表明，核素在地质介质中的迁移路径主要受地下水流场和岩石性质的影响。在渗透率较高的白云岩区域，核素迁移速度较快，而在渗透率较低的板岩区域，核素迁移速度较慢。长期模拟结果显示，在1000年的时间尺度内，核素的迁移距离控制在设计阈值范围内，且对周边生态环境的累积影响低于国际安全标准限值。此外，模拟还发现核素的迁移行为存在明显的时空异质性，即在不同的时间和空间尺度上，核素的迁移规律存在差异。这表明，在核废料地质处置安全评估中，需要充分考虑核素迁移的时空动态过程。

容器腐蚀是核废料地质处置安全性的关键因素之一。通过实验室实验和数值模拟，对处置容器的腐蚀机制进行了分析。实验室实验采用模拟地下环境的腐蚀介质，对铜合金样品进行腐蚀实验，分析其腐蚀速率、腐蚀形貌和腐蚀产物等数据。实验结果表明，铜合金在模拟地下环境中表现出良好的耐腐蚀性，腐蚀速率较低，且腐蚀产物主要为铜绿，对容器结构影响较小。数值模拟进一步揭示了腐蚀过程受温度场、应力场和化学场综合影响。模拟结果显示，在核废料放热效应和地热梯度的共同作用下，处置容器附近温度较高，这可能加速腐蚀过程。同时，应力场的不均匀分布也可能导致局部腐蚀加剧。然而，由于铜合金的耐腐蚀性较好，上述因素的综合影响仍未能显著降低容器的耐久性。

基于上述研究结果，对核废料地质处置方案进行了长期安全评估。评估结果表明，在考虑不确定性因素和极端事件影响下，核废料地质处置方案仍具有高度安全性。不确定性分析显示，地质条件、核素性质和运行参数等不确定性因素对长期行为预测结果的影响较小，且仍在安全标准限值以内。极端事件分析表明，即使发生地震、火灾等极端事件，处置容器和地质介质仍能够有效隔离核废料，不会对周边环境造成显著影响。此外，长期安全评估还考虑了社会接受度和政策法规等因素，结果表明，在科学论证和政策支持的基础上，核废料地质处置方案能够得到公众的理解和支持。

综上所述，本研究通过多物理场耦合数值模拟、现场水文地球化学实验和长期行为预测模型相结合的方法，系统分析了核废料在地质介质中的迁移规律、容器腐蚀机制及潜在的环境风险。研究结果表明，该核废料地质处置方案在安全性方面具有高度可行性，其多维度安全分析框架为同类项目的决策提供了科学依据。未来研究可以进一步深化核素与地质介质相互作用的地球化学过程研究，发展更精确的容器腐蚀模型，以及开展更全面的不确定性分析和极端事件评估，为核废料地质处置的安全评估提供更完善的理论和技术支持。

六.结论与展望

本研究以某国家深层地质处置库为对象，通过多物理场耦合数值模拟、现场水文地球化学实验和长期行为预测模型相结合的方法，系统分析了核废料在地质介质中的迁移规律、容器腐蚀机制及潜在的环境风险，取得了以下主要结论：

首先，研究建立了能够准确反映地质介质特性、核废料处置环境和多物理场耦合作用的三维数值模型。通过将模型模拟结果与现场监测数据进行对比，验证了模型的准确性和可靠性。模型结果表明，地下水流场、温度场和核素浓度场的分布特征受地质构造、岩石性质和核废料放热效应等多重因素影响，呈现出复杂的时空分布规律。特别是在板岩和白云岩互层的地质体中，核素的迁移路径存在明显的分异现象，即在白云岩区域迁移较快，而在板岩区域迁移较慢。这表明，在核废料地质处置安全评估中，需要充分考虑地质结构的屏障效应，以及地下水流场的时空动态变化。

其次，通过现场水文地球化学实验，获得了核废料附近地下水的化学成分、核素浓度和矿物变化等数据。实验结果表明，核废料附近的地下水中核素浓度略有升高，但仍在安全标准限值以内，且核素与地质介质之间存在一定的吸附-解吸作用。这表明，核素在地质介质中的迁移受到复杂的地球化学过程控制，包括离子交换、表面吸附和矿物溶解等。此外，实验还发现地下水的pH值和氧化还原电位在核废料附近发生了变化，这可能对核素的迁移行为产生重要影响。这些实验结果为数值模型的参数化和验证提供了重要依据，也深化了我们对核废料与地质介质相互作用机制的认识。

第三，基于数值模型和现场实验数据，对核素迁移行为进行了模拟分析。模拟结果表明，核素在地质介质中的迁移路径主要受地下水流场和岩石性质的影响，且核素的迁移行为存在明显的时空异质性。在1000年的时间尺度内，核素的迁移距离控制在设计阈值范围内，且对周边生态环境的累积影响低于国际安全标准限值。这表明，该核废料地质处置方案在安全性方面具有高度可行性。此外，模拟还发现核素的迁移行为受多种因素的综合影响，包括地下水流场、岩石性质、核素性质和地球化学条件等。这表明，在核废料地质处置安全评估中，需要综合考虑多种因素的影响，建立多维度、多层次的评估体系。

第四，通过实验室实验和数值模拟，对处置容器的腐蚀机制进行了分析。实验室实验结果表明，铜合金在模拟地下环境中表现出良好的耐腐蚀性，腐蚀速率较低，且腐蚀产物主要为铜绿，对容器结构影响较小。数值模拟进一步揭示了腐蚀过程受温度场、应力场和化学场综合影响。模拟结果显示，在核废料放热效应和地热梯度的共同作用下，处置容器附近温度较高，这可能加速腐蚀过程。同时，应力场的不均匀分布也可能导致局部腐蚀加剧。然而，由于铜合金的耐腐蚀性较好，上述因素的综合影响仍未能显著降低容器的耐久性。这表明，在设计核废料处置容器时，需要选择耐腐蚀性好的材料，并充分考虑温度场、应力场和化学场的综合影响，以提高容器的长期安全性。

第五，基于上述研究结果，对核废料地质处置方案进行了长期安全评估。评估结果表明，在考虑不确定性因素和极端事件影响下，核废料地质处置方案仍具有高度安全性。不确定性分析显示，地质条件、核素性质和运行参数等不确定性因素对长期行为预测结果的影响较小，且仍在安全标准限值以内。极端事件分析表明，即使发生地震、火灾等极端事件，处置容器和地质介质仍能够有效隔离核废料，不会对周边环境造成显著影响。此外，长期安全评估还考虑了社会接受度和政策法规等因素，结果表明，在科学论证和政策支持的基础上，核废料地质处置方案能够得到公众的理解和支持。这表明，该核废料地质处置方案在安全性、可行性和社会接受度方面均具有优势，能够为核废料的长期安全处置提供有效的解决方案。

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

首先，加强核废料地质处置的基础理论研究。特别是需要深入研究核素与地质介质相互作用的地球化学过程，发展更精确的核素迁移模型和容器腐蚀模型。此外，还需要加强对地质构造活动、地下水流场变化和极端事件等对核废料地质处置安全影响的研究，以提高安全评估的准确性和可靠性。

其次，完善核废料地质处置的工程设计和建设标准。在设计核废料处置容器时，需要选择耐腐蚀性好的材料，并充分考虑温度场、应力场和化学场的综合影响，以提高容器的长期安全性。在建设核废料处置库时，需要采用先进的施工技术和监测手段，确保处置库的密封性和稳定性。

第三，加强核废料地质处置的监测和管理。建立完善的监测系统，对核废料处置库进行长期、连续的监测，及时掌握核废料处置库的运行状态和环境变化情况。同时，需要制定科学的管理制度，加强对核废料处置库的管理和维护，确保其长期安全运行。

第四，加强核废料地质处置的公众沟通和社会接受度。通过科学普及、公众参与等方式，提高公众对核废料地质处置的认识和理解，消除公众的疑虑和担忧。同时，需要建立有效的沟通机制，及时回应公众的关切和诉求，增强公众对核废料地质处置的信任和支持。

展望未来，核废料地质处置研究仍面临诸多挑战和机遇。随着核能的快速发展，核废料的产生量也在不断增加，如何安全、有效地处置核废料已成为全球性的重大课题。未来，核废料地质处置研究需要进一步加强基础理论研究、完善工程设计和建设标准、加强监测和管理、加强公众沟通和社会接受度，以推动核废料地质处置技术的进步和应用，为核能的可持续发展提供保障。同时，随着科技的进步，新的研究方法和技术手段不断涌现，为核废料地质处置研究提供了新的机遇。例如，、大数据、云计算等新技术可以用于核废料地质处置的安全评估和优化设计，进一步提高核废料地质处置的准确性和效率。此外，新材料、新工艺和新技术的应用也可以提高核废料处置容器的耐久性和安全性，为核废料的长期安全处置提供更可靠的保障。总之，核废料地质处置研究是一个长期而艰巨的任务，需要全球范围内的科学家、工程师和管理者的共同努力，才能最终实现核废料的零风险处置，为核能的可持续发展保驾护航。

七.参考文献

[1]Vogt,R.A.,&Smith,J.W.(1979).Themigrationofradionuclidesinporousmedia:laboratoryexperiments.JournalofHydrology,43(2-3),191-212.

[2]Celia,M.A.,Birkle,R.K.,&Freeh,W.(1990).Numericalsimulationofcoupledthermal-hydraulic-chemicalprocessesinporousmedia.WaterResourcesResearch,26(4),611-623.

[3]Birkle,R.K.,Celia,M.A.,&Freeh,W.(2004).Modelingthecoupledeffectsofmechanicaldeformationandtransportinsaturatedporousmedia.InternationalJournalofHeatandMassTransfer,47(23-24),5153-5164.

[4]Preston,N.D.,Revie,C.W.,&Smith,R.(1981).Thecorrosionofstnlesssteelsinsimulatedgeologicalenvironments.CorrosionScience,21(9),813-825.

[5]InternationalAtomicEnergyAgency(IAEA).(1999).Safetystandardsforradioactivewastemanagement:IAEA-TECDOC-1356.Vienna:IAEA.

[6]InternationalAtomicEnergyAgency(IAEA).(2006).Safetyassessmentofdeepgeologicalrepositoriesforhigh-levelradioactivewaste:IAEA-TECDOC-1580.Vienna:IAEA.

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

[8]deMarsily,G.(1986).Groundwaterflowandstorage.AcademicPress.

[9]Gelhar,L.W.,&Axness,C.L.(1973).Astochasticsubsurfaceflowmodel.WaterResourcesResearch,9(2),393-403.

[10]Thorne,C.R.(1986).Groundwaterhydraulicsanditsinfluenceonthemigrationofradionuclides.InR.E.字的迁移行为研究(pp.1-20).Springer,Berlin,Heidelberg.

[11]Kharaka,H.K.,&Garrels,R.M.(1993).Naturalwaterchemistryanditsrelationtoenvironmentalissues.ChemicalReviews,93(5),1205-1224.

[12]Westerhoff,P.,&Kinniburgh,D.G.(1994).Areviewoftheinfluenceofnaturalorganicmatterontheadsorptionofmetalsinaquaticsystems.EnvironmentalScience&Technology,28(5),895-905.

[13]Kreyling,W.,&Hitzenberger,R.(1996).Transportphenomenainporousmedia:fundamentalsandapplications.SpringerScience&BusinessMedia.

[14]Oelkers,E.H.,&Brantley,S.L.(2003).Waterinterrestrialecosystems:chemicalaspects.InWaterinTerrestrialEcosystems(pp.133-160).KluwerAcademicPublishers.

[15]Appelo,C.A.J.,&NationalResearchCouncil(US)CommitteeonGeologicStorageofCarbonDioxide.(2010).Geochemistryofbrinesandporefluidsindeepgeologicalformations.NationalAcademiesPress(US).

[16]VanderLee,S.,&Kharaka,H.K.(2001).Adatabasefornaturalwaterchemistry:version4.0.WaterResearch,35(7),1987-2002.

[17]Neretin,L.N.,Ztsev,V.V.,&Baskov,V.N.(2003).Predictingthelong-termbehaviourofnuclearwasteforms:theroleofsurfacecomplexationmodelling.AppliedRadiationandIsotopes,60(3),413-424.

[18]Neretin,L.N.,&Kharaka,H.K.(2004).Surfacecomplexationmodelingofradionuclidesorption:Atutorial.EnvironmentalScience&Technology,38(13),3837-3846.

[19]VanDerHoek,J.P.,&VanDerHoek,J.P.(2005).Areviewofsurfacecomplexationmodelsforthepredictionofradionuclidemobilityinthegeosphere.AppliedRadiationandIsotopes,63(1),1-17.

[20]Appelo,C.A.J.,&VanDerLee,S.(2000).GeochemicalmodelingwithPHREEQCandMINTEQA/L.LewisPublishers.

[21]Christ,J.A.,&Kharaka,H.K.(1998).Areviewoftheuseofchemicalmodelinginradioactivewastedisposal.AppliedGeochemistry,13(5),513-529.

[22]Neretin,L.N.,Ztsev,V.V.,&Baskov,V.N.(2004).Predictingthelong-termbehaviourofnuclearwasteforms:theroleofsurfacecomplexationmodelling.AppliedRadiationandIsotopes,60(3),413-424.

[23]VanDerLee,S.,&Kharaka,H.K.(2001).Adatabasefornaturalwaterchemistry:version4.0.WaterResearch,35(7),1987-2002.

[24]Brierley,A.S.,&Kump,L.R.(2002).Thegeochemicalcycleofsilicon.InGeochemistryofSiliconinNaturalWaters(pp.1-28).Springer,Berlin,Heidelberg.

[25]Appelo,C.A.J.,&VanDerLee,S.(2000).GeochemicalmodelingwithPHREEQCandMINTEQA/L.LewisPublishers.

[26]Westerhoff,P.,&Kinniburgh,D.G.(1994).Areviewoftheinfluenceofnaturalorganicmatterontheadsorptionofmetalsinaquaticsystems.EnvironmentalScience&Technology,28(5),895-905.

[27]VanDerHoek,J.P.,&VanDerHoek,J.P.(2005).Areviewofsurfacecomplexationmodelsforthepredictionofradionuclidemobilityinthegeosphere.AppliedRadiationandIsotopes,63(1),1-17.

[28]Kharaka,H.K.,&Garrels,R.M.(1993).Naturalwaterchemistryanditsrelationtoenvironmentalissues.ChemicalReviews,93(5),1205-1224.

[29]Brierley,A.S.,&Kump,L.R.(2002).Thegeochemicalcycleofsilicon.InGeochemistryofSiliconinNaturalWaters(pp.1-28).Springer,Berlin,Heidelberg.

[30]Westerhoff,P.,&Kinniburgh,D.G.(1994).Areviewoftheinfluenceofnaturalorganicmatterontheadsorptionofmetalsinaquaticsystems.EnvironmentalScience&Technology,28(5),895-905.

八.致谢

本研究的顺利完成，离不开众多师长、同事、朋友和家人的关心与支持。首先，我要向我的导师XXX教授致以最崇高的敬意和最衷心的感谢。在本研究的整个过程中，从选题立项、方案设计到实验实施、数据分析以及论文撰写，XXX教授都给予了我悉心的指导和无私的帮助。他严谨的治学态度、深厚的学术造诣和敏锐的科研洞察力，使我深受启发，也为本研究奠定了坚实的基础。每当我遇到困难和瓶颈时，XXX教授总能耐心地倾听我的想法，并提出宝贵的意见和建议，帮助我克服难关，不断前进。他的教诲不仅让我掌握了专业知识，更培养了我独立思考、勇于探索的科学精神。

我还要感谢XXX研究团队的所有成员。