钍基次临界氯盐堆：概念、设计与前景探索

钍基次临界氯盐堆：概念、设计与前景探索一、绪论1.1研究背景与意义随着全球经济的持续增长和人口的不断增加，能源需求呈现出迅猛增长的态势。传统化石能源如煤炭、石油和天然气，不仅储量有限，且在使用过程中会产生大量的温室气体和污染物，对环境造成严重的负面影响，加剧全球气候变化。在此背景下，开发清洁、高效、可持续的能源成为当务之急。核能作为一种低碳、高效的能源，在全球能源结构中占据着日益重要的地位。自20世纪50年代第一座核电站建成以来，核能发电技术得到了广泛的应用和发展。截至目前，全球已有30多个国家和地区拥有核电站，核能发电量占全球总发电量的10%左右。然而，当前主流的核能系统，如压水堆、沸水堆等，主要以铀为核燃料，存在着铀资源有限、核废料处理困难、核扩散风险等问题。随着铀资源的逐渐消耗，这些问题愈发凸显，严重制约了核能的可持续发展。钍基次临界氯盐堆作为一种新型的核能系统，近年来受到了国际核能界的广泛关注。它以钍为核燃料，相较于传统的铀基反应堆，具有诸多独特的优势。钍是一种比铀更为丰富的天然放射性元素，其在地壳中的含量约为铀的3-4倍。这使得钍基反应堆在燃料供应上具有更大的优势，能够有效缓解铀资源短缺的问题，为核能的长期可持续发展提供坚实的保障。钍基次临界氯盐堆在运行过程中产生的核废料量较少，且放射性较低，处理难度和成本相对较低。这在很大程度上减轻了核废料对环境和人类健康的潜在威胁，降低了核能发展的环境风险。钍基次临界氯盐堆采用次临界运行模式，通过外部中子源来维持链式反应，具有更高的固有安全性。这种运行模式可以有效避免反应堆失控和核事故的发生，如切尔诺贝利核事故和福岛核事故等，极大地提高了核能利用的安全性。该反应堆还具有热效率高、运行灵活性强等优点，能够更好地满足不同的能源需求场景，为核能的多元化应用提供了可能。在全球积极推动能源转型和可持续发展的大背景下，研究钍基次临界氯盐堆的概念设计具有极其重要的意义。从能源安全角度来看，开发钍基核能系统有助于减少对进口铀资源的依赖，增强国家的能源自主供应能力，保障能源安全。从环境保护角度而言，钍基次临界氯盐堆产生的核废料少，对环境的影响小，符合全球对低碳、清洁能源的发展需求，有助于应对气候变化和实现可持续发展目标。从科学技术发展角度出发，研究钍基次临界氯盐堆可以推动核能技术的创新和进步，拓展核能领域的研究边界，为未来核能的发展开辟新的道路，提升国家在国际核能领域的竞争力。1.2钍资源与钍铀循环钍（Thorium）是一种天然放射性金属元素，化学符号为Th，原子序数90，在元素周期表中属于锕系元素。其在室温下呈银白色，质地柔软且具有一定的延展性。钍的熔点为1750℃，沸点达4788℃，密度约为11.7g/cm³。它的化学性质较为活泼，能与氧、卤素单质、酸等发生化学反应，比如粉末状的金属钍在空气中燃烧会生成ThO₂并发出白光。自然界中，钍主要以²³²Th同位素的形式存在，半衰期长达1.40×10¹º年，这使得钍在长期内保持相对稳定的放射性状态。钍在地壳中的平均含量为0.001-0.0012%，丰度比铀高，是一种潜在的丰富能源来源。全球钍资源量约为640万吨，广泛分布于二十多个国家。印度是钍储量最大的国家，其钍资源量约为85万吨，占世界已知钍资源总量的13%；巴西和澳大利亚的钍储量分别为63万吨和60万吨，各占10%和9%。中国也是钍资源较为丰富的国家之一，内蒙古白云鄂博矿区的钍矿资源类占中国总钍资源量一半以上，该矿床是独居石和氟碳铈矿的混合型轻稀土矿床，伴生着钍和铀核素。钍铀燃料循环是利用钍资源的关键过程。在这个循环中，原本不易发生核裂变的钍-232，在反应堆内通过吸收一个中子转变为钍-233。由于钍-233不稳定，会进一步发生β衰变，从而转化为铀-233。而铀-233是一种良好的裂变材料，能够在反应堆中引发链式反应，释放出巨大的核能。与传统的铀钚燃料循环相比，钍铀燃料循环具有诸多显著优势。钍资源的储量丰富，这使得钍铀燃料循环在燃料供应上更具可持续性，能有效缓解铀资源短缺带来的能源压力。钍铀燃料循环产生的核废料相对较少，且放射性水平较低，大大降低了核废料处理的难度和对环境的潜在危害。在反应堆中，钍燃料可以有多种利用方式。在熔盐堆中，钍基燃料可以溶解在熔融的盐载体中，形成液态燃料。这种液态燃料具有良好的流动性，能在反应堆堆芯和热交换器组成的回路中循环流动，将裂变产生的热量高效地输送到堆外。液态燃料还能实现在线燃料装卸和核素分离等功能，提高反应堆的运行效率和灵活性。在固态燃料反应堆中，钍可以与其他材料混合制成燃料棒，放入反应堆堆芯进行辐照。这种方式与传统的铀基燃料棒类似，但由于钍的特性，能够带来更好的中子经济性和较低的放射性废物产生量。钍基反应堆利用钍燃料具有明显的优势。从资源可持续性角度看，丰富的钍资源确保了核能发展在燃料供应上的长期稳定性，减少了对有限铀资源的依赖，为核能的长远发展提供了坚实的物质基础。在环境保护方面，较少的核废料产生量和较低的放射性水平，降低了核能发展对环境的负面影响，符合可持续发展的理念。从核安全角度而言，钍基反应堆的物理特性使得其在运行过程中具有更好的稳定性和安全性，降低了核事故发生的风险，增强了公众对核能利用的信心。1.3熔盐堆钍利用概述熔盐堆（MoltenSaltReactor,MSR）是一种采用熔融态的混合盐作为主冷却剂和核燃料载体的核裂变反应堆。其基本原理是利用核燃料在熔盐中的溶解和循环，实现核裂变反应的持续进行，并通过熔盐将裂变产生的热量带出堆芯，用于发电或其他工业应用。与传统的轻水堆相比，熔盐堆具有独特的技术特点。熔盐堆可以在高温下运行，且熔盐的蒸汽压较低，这使得系统可以在相对较低的压力下工作，从而降低了机械应力，提高了安全性，减少了对厚重压力容器的需求，降低了反应堆的建设成本和复杂性。熔盐堆的核燃料可以溶解在熔盐中，实现了燃料的在线添加和后处理，无需像传统反应堆那样频繁停堆更换燃料，提高了反应堆的运行效率和灵活性。熔盐堆的概念最早于20世纪40年代提出，美国在20世纪50-60年代开展了一系列熔盐堆相关的研究项目。1954年，美国进行的飞行器反应堆实验（AircraftReactorExperiment,ARE）是熔盐堆研究的重要开端，该实验旨在开发一种可用于核动力轰炸机的高功率密度反应堆，采用了熔融氟盐（NaF-ZrF4-UF4）作为燃料，氧化铍（BeO）作为慢化剂，液态钠作为二次冷却剂。此后，橡树岭国家实验室（OakRidgeNationalLaboratory,ORNL）在熔盐堆研究方面取得了重要进展。1965年，ORNL建成了熔盐堆实验（Molten-SaltReactorExperiment,MSRE），这是一座7.4MW热功率的试验堆，采用LiF-BeF2-ZrF4-UF4作为燃料，石墨作为慢化剂，运行温度达到650℃，并成功运行了四年。MSRE验证了熔盐堆的可行性和许多关键技术，为后续研究奠定了基础。然而，由于当时的技术条件、成本、政策等多方面因素，美国在20世纪70年代后期逐渐放缓了熔盐堆的研究。除美国外，其他国家也在熔盐堆领域开展了研究工作。法国在20世纪60-70年代对熔盐堆进行了理论研究和概念设计，探索了熔盐堆在核能领域的应用潜力。日本、俄罗斯等国也在不同时期参与了熔盐堆相关的国际合作研究项目，共同推动熔盐堆技术的发展。近年来，随着全球对清洁能源需求的增加以及对核能安全性和可持续性的关注，熔盐堆作为第四代核能系统的候选堆型之一，再次成为国际核能研究的热点。许多国家和国际组织纷纷加大对熔盐堆研究的投入，开展了一系列的研究项目和实验堆建设计划。目前，国际上多个国家和地区正在积极推进熔盐堆的研究与开发。美国重启了熔盐堆相关研究项目，旨在开发先进的熔盐堆技术，提高核能利用效率和安全性。美国能源部支持的一些研究项目聚焦于新型熔盐堆的设计、材料研发以及燃料循环技术的创新。欧洲一些国家通过合作项目共同开展熔盐堆研究，如欧盟的一些框架计划中包含了熔盐堆相关的研究课题，涉及熔盐堆的物理设计、热工水力分析、材料性能研究等多个方面。中国在熔盐堆领域也取得了显著进展。2011年，中国科学院启动了“未来先进核裂变能——钍基熔盐堆核能系统”战略性先导科技专项，旨在掌握钍基熔盐堆核能系统的核心技术，实现钍基熔盐堆的应用。2017年，中国在甘肃武威开工建设2兆瓦钍基熔盐实验堆（TMSR-SFR），并于2021年成功实现首次临界，标志着中国在钍基熔盐堆技术方面取得了重要突破。该实验堆的成功运行，为后续的技术发展和工程应用奠定了坚实基础，也表明中国在熔盐堆研究领域已处于世界前列。尽管熔盐堆在钍利用方面具有巨大潜力，但目前仍面临一些问题和挑战。熔盐堆的材料问题是一大挑战。熔盐具有强腐蚀性，在高温和辐照环境下，与熔盐接触的结构材料和包壳材料会受到严重腐蚀，影响反应堆的安全运行和寿命。开发能够耐受高温、辐照和熔盐腐蚀的新型材料是当前研究的重点之一。熔盐堆的燃料循环和后处理技术还不够成熟。实现高效、安全的钍铀燃料循环以及熔盐燃料的在线后处理，需要进一步的技术研发和验证。熔盐堆的核安全分析和监管技术也需要进一步完善，以确保反应堆在各种工况下的安全运行。此外，熔盐堆从实验研究到商业化应用还需要大量的资金投入和较长的研发周期，涉及技术、经济、社会等多方面的因素，需要综合考虑和协调推进。1.4研究目标与内容本研究旨在深入开展钍基次临界氯盐堆的概念设计，为新型核能系统的开发提供坚实的理论基础和技术支持。通过对反应堆物理、热工水力、材料科学等多学科的综合研究，优化反应堆的设计参数，提高其安全性、经济性和可持续性，使其成为一种具有竞争力的未来核能解决方案。具体研究内容如下：开发相关计算程序：基于反应堆物理和热工水力的基本原理，自主开发适用于钍基次临界氯盐堆的计算程序。在反应堆物理计算方面，利用蒙特卡罗方法，开发能够精确模拟中子输运过程的程序，考虑钍基燃料的核反应截面、中子散射特性等因素，准确计算堆芯内的中子通量分布、反应性变化等关键物理参数。在热工水力计算方面，开发能够模拟熔盐在复杂流道中流动和传热的程序，考虑熔盐的热物理性质随温度和压力的变化，以及堆芯内的发热功率分布，精确计算熔盐的温度场、流速场等热工参数。通过实验数据对程序进行验证和校准，确保程序的准确性和可靠性，为后续的堆芯设计和分析提供有力的工具。深入分析散裂反应：利用开发的计算程序，对散裂反应靶进行全面深入的分析。研究不同靶材料（如铅铋合金、钨等）在高能质子束轰击下的散裂反应特性，包括散裂中子产额、能量分布、空间分布等。分析靶的结构设计（如靶的形状、尺寸、冷却方式等）对散裂反应的影响，通过优化靶的结构，提高散裂中子的产生效率，降低靶材料的辐照损伤和热应力。研究散裂反应产生的热量在靶内的分布和传递规律，为靶的冷却系统设计提供依据，确保靶在高功率运行条件下的安全性和稳定性。优化堆芯设计：进行堆芯的初步设计，确定堆芯的结构形式（如圆柱形、棱柱形等）、燃料布置方式（如均匀布置、分区布置等）以及慢化剂、反射层的材料和结构。基于反应堆物理计算，优化堆芯的燃料组成和装载量，提高反应堆的增殖性能和中子经济性，确保反应堆能够在次临界状态下稳定运行，并实现钍铀燃料的有效转换。通过热工水力计算，优化堆芯的冷却剂流量分配和流道设计，确保堆芯内的温度分布均匀，避免局部过热现象的发生，提高反应堆的热效率和安全性。综合考虑反应堆物理和热工水力的要求，对堆芯设计进行多目标优化，寻求最佳的设计方案。开展安全分析：对钍基次临界氯盐堆进行全面的安全分析，识别可能出现的各种事故工况，如冷却剂丧失事故、反应性引入事故、超功率事故等。针对不同的事故工况，利用开发的计算程序，分析事故的发展过程和后果，评估反应堆的安全性。研究反应堆的固有安全特性，如负反应性温度系数、功率自动调节能力等，以及安全系统（如紧急停堆系统、余热排出系统等）的设计和性能，确保反应堆在各种事故工况下能够自动停堆，并有效排出堆芯余热，防止事故的扩大，保障公众和环境的安全。评估经济性：对钍基次临界氯盐堆的经济性进行评估，分析反应堆的建设成本、运行成本、维护成本以及核燃料循环成本等。考虑钍资源的价格、燃料制备和后处理的成本、设备的使用寿命和更换周期等因素，建立经济模型，预测反应堆在其寿命期内的发电成本。与传统的核能系统和其他可再生能源系统进行经济对比分析，评估钍基次临界氯盐堆在能源市场中的竞争力，为其商业化应用提供经济可行性依据。二、次临界系统燃耗程序开发2.1燃耗计算程序MOBAT在反应堆物理研究领域，MCNP5和ORIGEN2程序是两款重要的工具，各自具备独特的功能与特点。MCNP5（MonteCarloN-Particleversion5）是由美国洛斯阿拉莫斯国家实验室（LANL）开发的一款蒙特卡罗粒子输运模拟软件，在核工程、辐射防护、医学物理等多个领域有着广泛应用。它能模拟中子、光子、电子等多种粒子在复杂三维几何结构中的输运行为和相互作用，拥有强大的几何建模能力，可处理各种复杂的反应堆堆芯结构，如不同形状的燃料组件、复杂的冷却剂流道以及多样化的屏蔽结构等。MCNP5还具备丰富的核数据库，涵盖了众多核素的反应截面数据，能精确模拟各类核反应过程，如中子的散射、吸收、裂变等反应。它采用蒙特卡罗方法，通过大量的随机抽样和模拟粒子的运动轨迹，来获取系统的平均特性，从而为反应堆物理分析提供准确的数据支持。ORIGEN2程序则是一款功能强大的点燃耗及放射性衰变计算程序，主要用于计算反应堆中核素的燃耗、放射性衰变以及相关的核反应过程。它包含了完整的衰变链、裂变产额、各种核反应截面及其释放能等数据，能够详细追踪核素在反应堆运行过程中的变化情况。在计算核燃料的燃耗时，ORIGEN2可以考虑多种因素，如不同核素的中子吸收和裂变反应、燃料的初始成分和富集度、反应堆的运行功率和时间等。通过输入这些参数，程序能够准确计算出不同燃耗时刻核素的种类和数量，以及放射性废物的产生情况，为反应堆的燃料管理和核废物处理提供重要依据。然而，在次临界系统的燃耗计算中，单一使用MCNP5或ORIGEN2程序存在一定的局限性。MCNP5虽然能精确模拟粒子输运，但缺乏对核素燃耗过程的直接计算能力；ORIGEN2虽擅长燃耗计算，但在处理复杂几何结构和粒子输运方面能力有限。为了满足次临界系统燃耗计算的需求，开发了MOBAT耦合程序，将MCNP5和ORIGEN2的优势结合起来。MOBAT耦合程序的工作原理基于两者的互补特性。在计算过程中，首先利用MCNP5模拟次临界系统中粒子的输运过程，获取堆芯内精确的中子通量分布。中子通量分布是核反应发生的关键参数，它决定了核素与中子相互作用的概率。MCNP5通过对大量粒子的随机抽样和追踪，考虑了粒子在复杂堆芯结构中的散射、吸收等各种行为，从而得到准确的中子通量分布结果。然后，将MCNP5计算得到的中子通量分布数据传递给ORIGEN2程序。ORIGEN2根据这些中子通量数据，结合自身所包含的核反应截面、衰变链等信息，计算核素的燃耗变化，包括核素的生成、消耗以及放射性衰变等过程。在计算过程中，ORIGEN2会根据不同核素的特性和中子通量，精确计算每个核素的反应率和变化量。完成燃耗计算后，ORIGEN2将更新后的核素成分信息反馈给MCNP5。MCNP5利用这些新的核素成分数据，重新进行粒子输运模拟，考虑新的核素分布对中子输运的影响。如此循环迭代，逐步模拟次临界系统在燃耗过程中的物理变化。MOBAT耦合程序具有显著的优势。它能够更准确地模拟次临界系统的燃耗过程。通过结合MCNP5在粒子输运模拟的精确性和ORIGEN2在燃耗计算方面的专业性，避免了单一程序在处理复杂物理过程时的局限性，提高了计算结果的准确性和可靠性。在模拟钍基次临界氯盐堆的燃耗过程中，MOBAT耦合程序可以精确考虑堆芯内复杂的几何结构、钍基燃料的独特核反应特性以及燃耗过程中核素的变化，为反应堆的设计和运行提供更精确的数据支持。MOBAT耦合程序还能提高计算效率。相比于分别使用MCNP5和ORIGEN2进行独立计算，耦合程序通过数据的自动传递和迭代计算，减少了人工干预和数据处理的工作量，缩短了计算时间，提高了工作效率。这对于大规模的反应堆物理分析和优化设计具有重要意义，能够在更短的时间内为研究人员提供所需的计算结果，促进研究工作的快速推进。2.2次临界装置燃耗计算程序MOADS在次临界系统中，外源中子效率和加速器束流是两个至关重要的概念，它们对反应堆的运行性能和能量输出有着深远的影响。外源中子效率是指在次临界系统中，由外部中子源引入的中子能够有效引发核裂变反应并维持链式反应的比例。它反映了外源中子在反应堆中的利用效率，是衡量次临界系统性能的关键指标之一。较高的外源中子效率意味着更多的外源中子能够参与到核裂变反应中，从而提高反应堆的功率输出和能量转换效率。加速器束流则是指加速器产生并注入到次临界系统中的带电粒子束，通常为质子束。加速器束流的参数，如束流强度、能量等，直接影响着散裂反应产生的外源中子的数量和能量分布，进而影响次临界系统的运行性能。较强的束流强度能够产生更多的外源中子，为反应堆提供更充足的中子源，有助于维持反应堆的稳定运行。为了准确计算次临界装置的燃耗过程，开发了专门的程序MOADS。MOADS程序的开发是一个复杂而系统的工程，它基于对次临界系统物理过程的深入理解和现有计算方法的优化。在开发过程中，首先对次临界系统的燃耗物理过程进行了详细的理论分析。考虑到次临界系统中中子通量分布的特殊性，以及核素燃耗过程中的各种复杂核反应，建立了精确的物理模型。该模型充分考虑了中子的散射、吸收、裂变等反应，以及核素的衰变、生成和消耗等过程。采用了先进的数值算法来求解物理模型中的方程。针对蒙特卡罗方法计算量大、计算时间长的问题，通过优化随机数生成算法、改进粒子跟踪策略等方式，提高了计算效率。在核数据库的选择和处理方面，精心挑选了包含丰富核素信息和精确反应截面数据的数据库，并对数据进行了预处理和验证，以确保计算结果的准确性。经过多轮的测试和优化，最终完成了MOADS程序的开发。为了验证MOADS程序的准确性和可靠性，选取了IAEA公布的IAEA-ADS基准题模型进行验证计算。该基准题模型包含了详细的堆芯结构、材料组成、中子源等信息，是国际上广泛认可的用于验证次临界系统计算程序的标准模型。在验证计算过程中，将MOADS程序计算得到的有效增殖因数、燃耗区核素质量、外源中子效率以及外中子源强度等结果与基准题模型的参考值进行了对比分析。计算结果显示，MOADS程序计算得到的有效增殖因数与参考值的偏差在允许范围内，表明程序能够准确计算次临界系统的临界状态。对于燃耗区核素质量的计算，MOADS程序的结果与参考值也具有良好的一致性，能够精确追踪核素在燃耗过程中的变化。在外源中子效率和外中子源强度的计算上，MOADS程序同样表现出了较高的准确性，计算结果与参考值相符。通过与IAEA-ADS基准题模型的对比验证，充分证明了MOADS程序在计算次临界装置燃耗时的准确性和可靠性。这为钍基次临界氯盐堆的设计和分析提供了有力的工具，能够为反应堆的物理设计、燃料管理、安全分析等提供准确的数据支持。2.3本章小结本章围绕次临界系统燃耗程序开发展开了深入研究，取得了一系列具有重要意义的成果。针对次临界系统燃耗计算的需求，开发了MOBAT耦合程序，将MCNP5强大的粒子输运模拟能力与ORIGEN2在燃耗计算方面的优势相结合。通过MCNP5精确模拟次临界系统中粒子的输运过程，获取堆芯内准确的中子通量分布；再将该数据传递给ORIGEN2，利用其丰富的核数据库和专业算法计算核素的燃耗变化。这种耦合方式有效克服了单一程序在处理次临界系统复杂物理过程时的局限性，大幅提高了计算结果的准确性和可靠性。为了进一步满足次临界装置燃耗计算的特殊要求，开发了专门的MOADS程序。在开发过程中，充分考虑了次临界系统中中子通量分布的独特性以及核素燃耗过程中的复杂核反应，建立了精确的物理模型。采用先进的数值算法和优化的核数据库，有效提高了计算效率和结果的准确性。通过与IAEA公布的IAEA-ADS基准题模型进行对比验证，MOADS程序在计算有效增殖因数、燃耗区核素质量、外源中子效率以及外中子源强度等关键参数时，与参考值具有良好的一致性，充分证明了其准确性和可靠性。MOADS程序的开发对钍基次临界氯盐堆的研究具有极其重要的意义。它为钍基次临界氯盐堆的设计和分析提供了强有力的工具，能够准确计算反应堆在不同运行工况下的燃耗情况，为反应堆的物理设计、燃料管理、安全分析等提供精确的数据支持。在堆芯设计过程中，利用MOADS程序可以优化燃料组成和装载量，提高反应堆的增殖性能和中子经济性。在安全分析方面，该程序能够帮助研究人员准确评估反应堆在各种事故工况下的安全性，为安全系统的设计和优化提供依据。MOADS程序还可以用于研究钍基次临界氯盐堆的燃料循环方案，评估不同燃料循环策略对反应堆性能和环境的影响，为实现高效、可持续的钍铀燃料循环提供技术支持。三、散裂反应原理及散裂靶中子学性能分析3.1散裂反应原理散裂反应是一个复杂的核物理过程，当高能粒子（如质子）与重原子核相互作用时，会引发一系列的核反应，产生散裂产物和散裂中子。这一过程在钍基次临界氯盐堆中起着关键作用，为反应堆提供外源中子，维持链式裂变反应。散裂反应的第一步是核内级联过程。当高能质子入射到重原子核时，由于质子的能量显著大于核内核子之间的相互作用能量，且其德布罗意波长小于核子之间的平均距离，质子每次只与核内的一个核子发生碰撞。这种碰撞是一种非弹性散射，质子将部分能量传递给被撞核子，使其获得较高的动能。被撞核子也会像入射质子一样，与核内其他核子继续发生碰撞，从而引发一系列的级联碰撞。在这个过程中，会有多个核子从原子核中被发射出来，这些核子携带一定的能量，离开靶核。整个核内级联过程发生的时间极短，大约在10⁻²²秒的量级。在这个极短的时间内，核内的能量和动量发生了快速的转移和重新分配，使得原子核的结构和状态发生了显著变化。经过核内级联过程后，靶核处于高激发态，此时进入去激发过程。去激发过程主要通过蒸发和裂变两种方式进行。蒸发过程类似于受热液体蒸发出分子的过程。处于高激发态的残留核，会以发射中子、质子、α粒子等轻粒子的形式释放能量，逐渐降低激发态。这些被发射出来的粒子就像从高温原子核中“蒸发”出来的一样。当残留核的激发能足够高时，还可能发生裂变。裂变过程中，残留核会分裂成两个或多个中等质量的碎片，同时释放出大量的能量和中子。蒸发和裂变这两种去激发方式存在一定的竞争关系，具体哪种方式占主导，取决于残留核的激发能、质量数以及核的结构等因素。当残留核的激发能较低时，蒸发过程相对更容易发生；而当激发能较高时，裂变的概率会增加。在散裂反应中，会产生各种散裂产物。这些散裂产物是由级联－蒸发过程或级联－蒸发－裂变过程产生的。散裂产物的质量数分布范围较广，从接近靶核质量数的核素到质量数较小的轻核素都有。以400兆电子伏的质子轰击铋-209原子核为例，散裂产物的质量数主要位于150-205之间。反应截面的峰位约在质量数为195处，峰的半高宽约为16个质量数。质量数低于峰位的核素的反应截面随质量数的下降而迅速下降。当入射质子能量升高时，峰位向质量数降低方向移动，峰高下降，而宽度增大。质量数较低的核素集中于β稳定线附近，质量数较高的产物则处于β稳定线的缺中子的一侧，这与蒸发过程密切相关，并与蒸发计算的结果基本相符。散裂中子是散裂反应的重要产物之一。每个与原子核相作用的质子能够轰击出20-30个中子。这些散裂中子包括级联中子、蒸发中子以及裂变中子。级联中子是在核内级联过程中发射出来的，它们具有较高的能量，通常在兆电子伏量级。蒸发中子是在蒸发过程中发射的，其能量相对较低。裂变中子则是在裂变过程中产生的，能量也较高。散裂中子的能谱分布较宽，从低能到高能都有分布。其能量分布和产额与入射质子的能量、靶核的性质以及散裂反应的具体过程等因素有关。随着入射质子能量的增加，散裂中子的产额和平均能量一般也会增加。不同靶核材料对散裂中子的产生也有显著影响，例如，液态铅和铅铋合金（LBE）靶由于其高密度和高原子序数，通常具有较高的中子产额。3.2散裂反应计算模型在散裂反应的研究中，为了准确描述和预测散裂反应的过程和结果，发展了多种计算模型，这些模型基于不同的物理假设和数学方法，对散裂反应的不同阶段和现象进行模拟。核内级联模型是描述散裂反应中核内级联过程的重要工具。在核内级联过程中，高能质子与靶核内的核子发生一系列的非弹性碰撞，导致多个核子从原子核中被发射出来。常见的核内级联模型有Bertini模型和INCL4模型。Bertini模型是较早发展起来的核内级联模型，它基于简单的物理图像，将核内级联过程视为一系列的二元碰撞。在这个模型中，入射质子与靶核内的核子发生碰撞，每次碰撞都遵循一定的能量和动量守恒规则。通过模拟这些碰撞过程，可以计算出在核内级联阶段发射出的核子的能量和角度分布。Bertini模型在早期的散裂反应研究中得到了广泛应用，对于理解核内级联过程的基本特征起到了重要作用。然而，该模型也存在一些局限性。它对核子间相互作用的描述相对简单，没有充分考虑核介质的影响，导致在某些情况下计算结果与实际情况存在偏差。在处理高能量的入射质子时，Bertini模型的计算精度会下降。INCL4模型则是在Bertini模型的基础上发展起来的，它对核内级联过程的描述更加精细。INCL4模型考虑了更多的物理因素，如核子在核介质中的多次散射、核子-核子相互作用的截面随能量的变化等。它采用了更复杂的数学方法来描述核内级联过程，能够更准确地计算出散裂反应中发射出的粒子的能量和角度分布。在模拟高能质子与重核的散裂反应时，INCL4模型能够更准确地预测级联中子的产额和能谱。INCL4模型也并非完美无缺。它的计算过程相对复杂，计算量较大，对计算机的性能要求较高。在处理一些特殊的核反应情况时，如非常低能量的入射质子或特殊的靶核结构，INCL4模型可能也会存在一定的局限性。蒸发模型主要用于描述散裂反应中靶核去激发过程中的蒸发阶段。在蒸发过程中，处于高激发态的残留核会以发射中子、质子、α粒子等轻粒子的形式释放能量，逐渐降低激发态。常见的蒸发模型有Dresner模型和ABLA模型。Dresner模型是一种基于统计理论的蒸发模型，它假设残留核处于热平衡状态，发射粒子的过程类似于热蒸发。在Dresner模型中，通过计算残留核的激发能、温度等参数，根据统计力学的原理来确定发射粒子的概率和能量分布。该模型在计算蒸发粒子的产额和能谱方面具有一定的准确性，能够较好地描述蒸发过程的一些基本特征。Dresner模型对残留核的结构和性质的描述相对简单，没有充分考虑核结构对蒸发过程的影响。在处理一些具有特殊核结构的靶核时，Dresner模型的计算结果可能与实际情况存在一定的偏差。ABLA模型则是一种更加先进的蒸发模型，它综合考虑了核结构、能级密度等因素对蒸发过程的影响。ABLA模型采用了更精确的能级密度公式和衰变概率计算方法，能够更准确地描述不同核素在蒸发过程中的行为。它不仅可以计算蒸发粒子的产额和能谱，还可以对蒸发过程中产生的各种核素的丰度进行预测。在研究散裂产物的分布时，ABLA模型能够提供更详细和准确的信息。ABLA模型的计算过程较为复杂，需要输入更多的核物理参数，这在一定程度上限制了它的应用范围。对于一些缺乏详细核数据的核素，ABLA模型的计算精度可能会受到影响。裂变模型用于描述散裂反应中靶核去激发过程中的裂变阶段。当残留核的激发能足够高时，裂变过程可能会发生，残留核会分裂成两个或多个中等质量的碎片，同时释放出大量的能量和中子。常见的裂变模型有基于液滴模型的裂变模型和考虑壳效应的裂变模型。基于液滴模型的裂变模型将原子核视为一个液滴，当液滴的能量超过一定阈值时，就会发生裂变。该模型通过计算液滴的表面能、库仑能等参数，来确定裂变的概率和裂变碎片的质量分布。这种模型在描述裂变过程的一些宏观特征方面具有一定的成功，但它对裂变过程中量子效应的考虑不足。考虑壳效应的裂变模型则在液滴模型的基础上，进一步考虑了原子核的壳结构对裂变的影响。原子核的壳结构会导致核的稳定性和裂变性质发生变化，考虑壳效应的裂变模型能够更准确地描述一些具有特殊壳结构的核素的裂变过程。该模型需要更精确的核结构数据，计算过程也更为复杂。在散裂反应中，蒸发和裂变这两种去激发方式存在竞争关系。当残留核的激发能较低时，蒸发过程相对更容易发生，因为发射轻粒子所需的能量较低。随着残留核激发能的增加，裂变的概率会逐渐增大，因为较高的激发能使得原子核更容易克服裂变位垒。蒸发和裂变的竞争还与靶核的质量数、中子数等因素有关。重核相对更容易发生裂变，因为它们的库仑能较高，裂变时释放的能量更多。而轻核则更倾向于通过蒸发来释放能量。准确描述蒸发和裂变的竞争关系是散裂反应计算模型中的一个关键问题，目前的模型在这方面还存在一定的改进空间，需要进一步深入研究核结构、激发能等因素对蒸发和裂变竞争的影响。3.3散裂反应计算模拟程序在散裂反应计算模拟领域，多种程序被广泛应用，它们各自具有独特的优势和适用范围，同时也存在一定的局限性。HETC（High-EnergyTransportCode）程序是早期用于高能粒子输运模拟的重要工具。它基于蒙特卡罗方法，能够模拟高能质子与原子核相互作用过程中的粒子输运行为。在散裂反应计算中，HETC程序能够考虑多种粒子，如质子、中子、π介子等在物质中的输运和相互作用。它在处理高能区的核反应时具有一定的优势，能够较为准确地计算出级联过程中粒子的能量和角度分布。由于其开发时间较早，在核反应模型和数据方面存在一定的局限性。HETC程序对核反应截面数据的描述相对简单，对于一些复杂的核反应过程，尤其是涉及到多种核素和复杂核结构的反应，其计算精度有待提高。它在处理低能区的核反应时，准确性不如专门针对低能核反应开发的程序。LAHET（LosAlamosHigh-EnergyTransport）程序也是一款常用的散裂反应计算程序。它同样基于蒙特卡罗方法，在散裂反应的模拟中具有较强的功能。LAHET程序包含了丰富的核反应模型，如核内级联模型、蒸发模型等，能够较为全面地描述散裂反应的过程。在计算散裂中子产额、能量沉积等关键参数方面，LAHET程序表现出较高的准确性。在研究高能质子轰击不同靶材料时，LAHET程序能够准确地计算出散裂中子的产额和能谱分布。它还可以与其他程序，如MCNP程序进行耦合，进一步提高计算的准确性和可靠性。LAHET程序的计算效率相对较低。由于散裂反应涉及到复杂的核物理过程，需要进行大量的粒子输运模拟，LAHET程序在处理大规模计算时，计算时间较长，对计算机资源的需求较大。它的输入参数和模型设置较为复杂，对于初学者来说，需要花费较多的时间和精力来掌握。MCNPX（MonteCarloN-ParticleeXtended）程序是在MCNP基础上发展起来的，具有更强大的功能和更广泛的应用范围。在散裂反应计算中，MCNPX程序不仅能够精确模拟中子、质子等粒子的输运过程，还能处理复杂的几何结构和多种材料的相互作用。它包含了最新的核反应数据和模型，能够更准确地描述散裂反应中的各种物理过程。MCNPX程序可以考虑材料的温度、密度等因素对核反应的影响，在模拟实际的散裂靶时，能够提供更符合实际情况的计算结果。它还具备强大的可视化功能，能够直观地展示粒子的输运轨迹和能量沉积分布。MCNPX程序的计算结果对输入数据的依赖性较强。如果输入的核数据不准确或不完整，可能会导致计算结果出现较大偏差。它在处理某些特殊的核反应或极端条件下的散裂反应时，可能需要进一步的模型改进和验证。在实际应用中，选择合适的散裂反应计算模拟程序需要综合考虑多种因素。对于研究高能区的散裂反应，且对计算效率要求不是特别高的情况，可以优先考虑LAHET程序或MCNPX程序，因为它们在处理复杂核反应过程和准确描述高能粒子输运方面具有优势。如果需要处理复杂的几何结构和多种材料的相互作用，MCNPX程序则是更好的选择。而对于一些对计算精度要求不是特别高，且主要关注散裂反应的基本特征和趋势的研究，HETC程序也可以作为一种快速计算和初步分析的工具。在进行散裂反应计算模拟时，还可以结合多种程序的优势，进行对比计算和验证，以提高计算结果的可靠性。3.4散裂靶中子学性能分析在散裂靶中子学性能分析中，计算方法与工具的选择至关重要。本研究采用蒙特卡罗方法，借助MCNPX程序进行模拟计算。蒙特卡罗方法基于概率统计理论，通过大量随机抽样来模拟粒子的输运过程，能够有效处理复杂的几何结构和物理过程。MCNPX程序则是一款功能强大的蒙特卡罗粒子输运模拟软件，它能够精确模拟中子、质子等粒子在物质中的输运和相互作用，具备丰富的核反应模型和详细的核数据库，为散裂靶中子学性能分析提供了有力的支持。为了准确模拟散裂靶的中子学性能，构建了合理的计算模型。散裂靶选用常见的材料，如液态铅、铅铋合金（LBE）以及熔盐等，并将其设计为圆柱形结构。以液态铅散裂靶为例，其半径设定为0.3m，高度为1.0m。在模型中，明确了质子束的入射参数，入射质子能量为1GeV，束流强度为1mA。质子束沿着靶的中心轴方向垂直入射，这种入射方式能够简化计算模型，同时也符合实际的加速器驱动系统中的常见情况。为了准确模拟粒子的输运过程，对计算区域进行了精细的网格划分，确保计算结果的准确性。通过MCNPX程序的模拟计算，得到了一系列关于散裂靶中子学性能的关键数据，并对这些结果进行了深入分析。在中子产额及中子能谱对比方面，不同靶材料展现出各异的特性。液态铅和LBE靶由于其高密度和高原子序数，具有较高的中子产额。在1GeV质子束轰击下，液态铅靶的中子产额约为25个/质子，LBE靶的中子产额约为23个/质子。熔盐靶的中子产额相对较低，例如氯盐靶的中子产额约为15个/质子。中子能谱分布也因靶材料的不同而有所差异。液态铅和LBE靶产生的中子能谱较宽，从低能到高能都有分布，且高能中子的比例相对较高。熔盐靶产生的中子能谱相对较窄，低能中子的比例相对较大。这种差异主要是由于不同靶材料的核反应特性和中子散射性质不同所导致的。液态铅和LBE靶的原子核质量较大，与质子相互作用时更容易产生高能中子；而熔盐靶的原子核质量相对较小，中子在其中的散射过程更为频繁，导致低能中子的比例增加。在质子入射高度选择上，研究发现其对散裂靶的中子学性能有着显著影响。通过改变质子入射高度，对比不同高度下的中子产额和能量沉积分布，发现当质子入射高度在靶中心位置附近时，中子产额达到最大值。当质子入射高度偏离中心位置时，中子产额逐渐降低。这是因为在中心位置入射时，质子能够与靶材料充分相互作用，引发更多的散裂反应，从而产生更多的中子。而偏离中心位置入射时，质子与靶材料的作用路径变短，部分质子可能未充分反应就离开了靶体，导致中子产额下降。能量沉积分布也会随着质子入射高度的变化而改变。在中心位置入射时，能量沉积相对较为均匀；偏离中心位置时，能量沉积会出现明显的不均匀分布，局部区域的能量沉积会显著增加，这可能会对靶材料的结构和性能产生不利影响。能量沉积计算结果显示，质子在散裂靶内的能量沉积主要集中在入射区域附近。以液态铅靶为例，在1GeV质子束轰击下，约70%的能量沉积在靶的前10cm区域内。随着深度的增加，能量沉积逐渐减少。能量沉积还与靶材料的密度和原子序数有关。高密度和高原子序数的材料，如液态铅和LBE靶，由于其原子核与质子的相互作用概率较大，能量沉积相对较多；而熔盐靶由于其密度和原子序数相对较低，能量沉积相对较少。能量沉积的不均匀分布可能会导致靶材料局部温度升高，产生热应力，进而影响靶的使用寿命和安全性。因此，在散裂靶的设计中，需要充分考虑能量沉积的分布情况，采取有效的冷却措施，以确保靶材料的温度在安全范围内。在散裂产物计算方面，通过MCNPX程序模拟得到了不同靶材料在散裂反应中产生的散裂产物的种类和丰度。液态铅靶和LBE靶产生的散裂产物种类繁多，涵盖了从质量数较小的轻核素到接近靶核质量数的重核素。熔盐靶产生的散裂产物相对较少，且主要集中在质量数较小的区域。散裂产物的分布与靶材料的核结构以及散裂反应过程密切相关。不同的靶材料具有不同的核结构和核反应特性，在散裂反应中，质子与靶核的相互作用方式和程度不同，导致产生的散裂产物也不同。散裂产物的放射性和半衰期也各不相同，这对于散裂靶的放射性废物处理和环境影响评估具有重要意义。一些半衰期较长的放射性散裂产物需要特殊的处理和处置方式，以确保环境安全。通过对散裂靶中子学性能的分析可知，液态铅和LBE靶在中子产额方面具有优势，但能量沉积不均匀，可能对靶的结构和冷却造成挑战。熔盐靶虽然中子产额相对较低，但能量沉积分布较为均匀，且在简化堆芯结构方面具有潜力。在实际的钍基次临界氯盐堆设计中，需要综合考虑中子学性能、工程可行性、安全性等多方面因素，选择合适的散裂靶材料和设计方案。3.5本章总结本章深入研究了散裂反应原理及散裂靶中子学性能，取得了一系列具有重要价值的成果。详细阐述了散裂反应的原理，当高能质子与重原子核相互作用时，会依次经历核内级联过程和去激发过程。在核内级联过程中，质子与靶核内的核子发生非弹性碰撞，引发一系列级联碰撞，导致多个核子从原子核中发射出来。随后的去激发过程则通过蒸发和裂变两种方式进行，残留核以发射中子、质子、α粒子等轻粒子的形式（蒸发）或分裂成中等质量碎片（裂变）来释放能量。在这个过程中，会产生各种散裂产物，其质量数分布范围较广，散裂中子也是重要产物之一，每个与原子核作用的质子能轰击出20-30个中子。对散裂反应计算模型进行了全面分析，包括核内级联模型（如Bertini模型和INCL4模型）、蒸发模型（如Dresner模型和ABLA模型）以及裂变模型。这些模型基于不同的物理假设和数学方法，分别描述了散裂反应中不同阶段的物理过程。核内级联模型用于描述核内级联过程中粒子的相互作用和发射，蒸发模型用于描述残留核去激发过程中的蒸发阶段，裂变模型则用于描述激发能足够高时残留核的裂变过程。不同模型在描述散裂反应的不同方面各有优势和局限性，Bertini模型对核内级联过程的描述相对简单，而INCL4模型考虑了更多物理因素，计算更精确，但计算量也更大。介绍了散裂反应计算模拟程序，如HETC、LAHET和MCNPX等。这些程序基于蒙特卡罗方法，能够模拟高能质子与原子核相互作用过程中的粒子输运行为。HETC程序在早期用于高能粒子输运模拟，具有一定的基础功能，但在核反应模型和数据方面存在局限性。LAHET程序包含丰富的核反应模型，计算准确性较高，但计算效率相对较低，输入参数和模型设置也较为复杂。MCNPX程序功能强大，能处理复杂几何结构和多种材料相互作用，包含最新核反应数据和模型，计算结果准确性高，且具备可视化功能，但对输入数据的依赖性较强。通过MCNPX程序对散裂靶中子学性能进行了深入分析。构建了合理的计算模型，明确了散裂靶材料（如液态铅、铅铋合金、熔盐等）的选择和结构设计，以及质子束的入射参数。在中子学性能分析中，发现不同靶材料的中子产额和中子能谱存在显著差异。液态铅和铅铋合金靶由于其高密度和高原子序数，中子产额较高，能谱较宽且高能中子比例相对较高；熔盐靶中子产额相对较低，能谱较窄且低能中子比例较大。质子入射高度对散裂靶的中子学性能有显著影响，在靶中心位置附近入射时，中子产额达到最大值，能量沉积也更为均匀，偏离中心位置入射则会导致中子产额下降和能量沉积不均匀。能量沉积主要集中在入射区域附近，且与靶材料的密度和原子序数有关，液态铅和铅铋合金靶能量沉积较多，熔盐靶相对较少。不同靶材料的散裂产物种类和丰度也各不相同，液态铅靶和铅铋合金靶产生的散裂产物种类繁多，熔盐靶产生的相对较少。这些研究成果对钍基次临界氯盐堆的设计具有重要的指导意义。在堆芯设计中，需要根据散裂靶的中子学性能选择合适的靶材料和结构，以提高中子产额，优化中子能谱，确保反应堆能够获得足够的外源中子，维持稳定的链式裂变反应。考虑能量沉积分布和散裂产物特性，合理设计散裂靶的冷却系统和放射性废物处理方案，以保障反应堆的安全运行和环境安全。通过对散裂反应和散裂靶中子学性能的深入研究，为钍基次临界氯盐堆的概念设计提供了坚实的理论基础和关键的数据支持，有助于推动钍基核能系统的发展和应用。四、钍基次临界氯盐堆中子学设计4.1次临界氯盐堆TMSFEA原理介绍能量放大器（EnergyAmplifier）的概念最早由诺贝尔物理学奖获得者卡罗・鲁比亚（CarloRubbia）提出。这一概念的核心是将加速器产生的高能质子束引入到次临界堆中，通过与堆内的散裂靶发生散裂反应，产生大量的外源中子，以此驱动次临界堆的链式裂变反应。传统的临界反应堆依靠堆芯内的核燃料自身的链式反应来维持运行，而能量放大器则通过外部中子源的引入，打破了这种依赖，为核能的利用提供了一种全新的思路。在传统临界反应堆中，核燃料的裂变反应是自持的，一旦达到临界状态，链式反应就会持续进行。但这种反应堆对核燃料的纯度和富集度要求较高，且在运行过程中存在一定的安全风险。能量放大器采用次临界运行模式，自身产生的中子不足以维持链式反应，需要依靠外部中子源。这种运行模式使得反应堆的反应性更容易控制，大大提高了反应堆的安全性。即使在发生意外情况时，只要切断外部中子源，反应堆的链式反应就会迅速停止，从而有效避免了核事故的发生。能量放大器概念的提出，为解决传统核能系统中存在的一些问题提供了可能。在燃料利用方面，能量放大器可以更有效地利用核燃料，特别是对于钍资源的利用。钍是一种储量丰富的核燃料，但在传统反应堆中难以直接利用。能量放大器通过散裂反应产生的中子，可以将钍转化为可裂变的铀-233，从而实现钍资源的高效利用。在核废料处理方面，能量放大器产生的核废料相对较少，且放射性水平较低，处理难度和成本也相应降低。这是因为能量放大器的运行模式和燃料利用方式，使得核废料中的长寿命放射性核素含量减少，降低了核废料对环境的长期影响。在次临界氯盐堆中，能量增益G是一个至关重要的参数，它反映了反应堆的能量产出与输入之间的关系。能量增益G的定义为反应堆输出的总能量与加速器输入的能量之比。在TMSFEA中，能量增益G可以通过以下公式计算：G=\frac{P_{th}}{P_{acc}}其中，P_{th}为反应堆的热功率，即反应堆输出的能量，单位为瓦特（W）；P_{acc}为加速器的束流功率，即加速器输入的能量，单位为瓦特（W）。能量增益G受到多种因素的影响。中子产额是一个关键因素。散裂反应产生的中子产额越高，反应堆内参与裂变反应的中子就越多，从而能够产生更多的能量，提高能量增益G。不同的散裂靶材料具有不同的中子产额，液态铅和铅铋合金（LBE）靶由于其高密度和高原子序数，通常具有较高的中子产额，能够为反应堆提供更多的中子，有利于提高能量增益G。而熔盐靶的中子产额相对较低，可能会对能量增益G产生一定的影响。反应堆的中子效率也会影响能量增益G。中子效率是指参与裂变反应的中子数与入射中子数之比。如果反应堆内存在较多的中子吸收体，或者中子泄漏较多，就会导致中子效率降低，从而减少反应堆的能量产出，降低能量增益G。堆芯的结构设计、燃料的布置方式以及反射层的使用等都会影响中子的分布和泄漏情况，进而影响中子效率和能量增益G。核燃料转换比CR是另一个重要的中子学参数，它对于评估反应堆的燃料增殖性能具有关键意义。核燃料转换比CR的定义为反应堆中产生的易裂变核素（如铀-233、钚-239等）的数量与消耗的易裂变核素数量之比。在TMSFEA中，核燃料转换比CR的计算公式为：CR=\frac{N_{new}}{N_{consumed}}其中，N_{new}为产生的易裂变核素的数量，N_{consumed}为消耗的易裂变核素的数量。核燃料转换比CR反映了反应堆将非易裂变核素（如钍-232、铀-238等）转化为易裂变核素的能力。当CR大于1时，意味着反应堆在运行过程中不仅能够消耗易裂变核素，还能够产生更多的易裂变核素，实现燃料的增殖。在钍基次临界氯盐堆中，通过合理的设计和运行，可以使CR大于1，从而实现钍资源的有效利用和增殖。这对于解决核燃料资源短缺问题具有重要意义，能够延长反应堆的运行寿命，减少对外部核燃料的依赖。如果CR小于1，反应堆则需要不断补充易裂变核素，以维持链式反应的进行，这会增加核燃料的成本和管理难度。TMSFEA的设计目标是在简化堆芯设计和燃料处理流程的前提下，实现钍资源的高效利用和稳定的能量输出。为了实现这一目标，TMSFEA采用了一系列独特的设计思路和关键参数。在堆芯设计方面，TMSFEA采用单流、无慢化体的简单堆芯设计。这种设计可以有效降低堆芯的复杂度，减少系统的不确定性和潜在风险。堆芯中的熔盐既作为燃料又作为散裂靶，这种一体化的设计避免了使用特殊的结构材料将靶与次临界堆芯分开，有助于简化整个堆芯的结构，降低系统复杂性。熔盐还具有良好的热物理性质，能够有效地传输热量，提高反应堆的热效率。在燃料选择上，TMSFEA选用具有钚-钍溶解度高、熔点合适以及散裂中子产额较高的三元氯盐NaCl-PuCl₃-ThCl₄作为散裂靶和燃料盐。这种燃料盐能够满足TMSFEA对燃料性能的多方面要求。较高的钚-钍溶解度可以保证燃料在熔盐中的均匀分布，有利于提高反应堆的中子经济性和运行稳定性。合适的熔点则确保了燃料盐在反应堆运行温度下能够保持液态，实现良好的流动性和传热性能。较高的散裂中子产额则有助于提高反应堆的能量增益和核燃料转换比。在能量输出方面，TMSFEA设定了明确的功率目标。在初始束流强度小于最大束流强度4mA的限定条件下，选定系统总功率为300MWth。这一功率水平既考虑了反应堆的技术可行性和安全性，又能够满足一定规模的能源需求。通过合理调整加速器的束流功率和堆芯的设计参数，可以实现反应堆在该功率水平下的稳定运行，持续输出能量。在钍资源利用方面，TMSFEA致力于实现高效的钍-铀循环。通过引入加速器提供的外中子源，补偿燃耗过程中裂变产物对中子的吸收，使得反应堆能够在无需燃料在线处理的条件下，实现钍资源的有效利用。在燃耗寿期结束时，有57%的钚和16%的钍被消耗，²³³U的裂变份额可以达到70.9%，表明TMSFEA在钍资源利用方面具有较高的效率。TMSFEA还通过多循环演化的燃料循环方案，进一步提高钍的利用效率，降低核废料的放射性毒性水平。在多循环演化过程中，系统的初始CR和钚的裂变率份额逐循环下降，而²³³U的裂变率份额逐循环上升。经过共四个循环的演化，TMSFEA的²³³U总净累积产量为2832.7kg，钍的利用效率基本随燃耗线性增长，四个循环寿期末的Th利用率分别可达11.38%、18.53%、23.42%和26.49%，钍的能量贡献从初循环寿期初的0%增加到第四循环寿期末的66.9%，而第四循环的乏燃料放射性毒性比初循环堆降低了63.7%。4.2TMSFEA的设计与优化4.2.1TMSFEA燃料盐在TMSFEA的设计中，燃料盐的选择是至关重要的环节，它直接影响着反应堆的性能、安全性和经济性。燃料盐需要具备多种优良特性，以满足反应堆复杂的运行需求。钚-钍溶解度高是燃料盐的重要特性之一。较高的溶解度能够确保钚和钍在熔盐中均匀分布，避免出现局部浓度过高或过低的情况。这对于维持反应堆内中子通量的均匀分布以及保证核反应的稳定进行至关重要。均匀的燃料分布可以提高中子的利用效率，减少中子的浪费，从而提高反应堆的中子经济性和运行稳定性。如果钚-钍在熔盐中溶解度低，可能会导致燃料沉淀或团聚，影响反应堆的正常运行。合适的熔点也是燃料盐的关键特性。TMSFEA运行在一定的温度范围内，燃料盐的熔点需要与之相匹配。合适的熔点能确保燃料盐在反应堆运行温度下保持液态，具有良好的流动性。液态的燃料盐可以在堆芯内自由流动，实现高效的热量传递，将核裂变产生的热量迅速带出堆芯，保证反应堆的热效率。如果燃料盐熔点过高，在反应堆运行温度下可能无法完全熔化，导致燃料流动不畅，影响热量传递，甚至可能引发局部过热，威胁反应堆的安全。而熔点过低，可能会增加燃料盐的挥发性和腐蚀性，对反应堆的材料和结构造成损害。散裂中子产额较高同样是燃料盐的重要考量因素。散裂中子是维持次临界堆链式裂变反应的关键，较高的散裂中子产额意味着反应堆能够获得更多的外源中子，有利于提高反应堆的能量增益和核燃料转换比。更多的散裂中子可以引发更多的核裂变反应，产生更多的能量，提高反应堆的输出功率。较高的散裂中子产额还可以促进钍-铀循环的进行，提高钍资源的利用效率。基于以上对燃料盐特性的要求，经过对多种熔盐的研究和分析，最终选定三元氯盐NaCl-PuCl₃-ThCl₄作为TMSFEA的散裂靶和燃料盐。这种三元氯盐在钚-钍溶解度、熔点以及散裂中子产额等方面都表现出了良好的性能。在钚-钍溶解度方面，它能够有效地溶解钚和钍，为反应堆提供均匀的燃料分布。其熔点经过测试和计算，符合TMSFEA的运行温度要求，在运行过程中能够保持稳定的液态，确保良好的流动性和传热性能。在散裂中子产额方面，三元氯盐NaCl-PuCl₃-ThCl₄展现出较高的数值，能够为反应堆提供充足的外源中子，满足反应堆对中子的需求。与其他常见的熔盐相比，如氟盐等，NaCl-PuCl₃-ThCl₄在钚-钍溶解度和散裂中子产额上具有明显的优势。氟盐虽然在某些方面也有良好的性能，但在钚-钍溶解度和散裂中子产额方面不如NaCl-PuCl₃-ThCl₄，这使得NaCl-PuCl₃-ThCl₄成为TMSFEA燃料盐的理想选择。4.2.2TMSFEA加速器质子能量与束流功率加速器质子能量和束流功率是影响TMSFEA性能的关键因素，它们与反应堆的能量增益和运行稳定性密切相关。为了确定合适的质子能量和束流功率，进行了一系列深入的研究。在研究质子能量与系统能量增益的关系时，通过改变质子能量进行模拟计算。当质子能量较低时，质子与散裂靶发生散裂反应产生的中子产额较低，导致反应堆内参与裂变反应的中子数量不足，能量增益较低。随着质子能量的增加，散裂反应更加剧烈，中子产额显著提高，反应堆内的裂变反应更加充分，能量增益逐渐增大。当质子能量过高时，虽然中子产额进一步增加，但同时也会带来一些负面效应，如散裂靶的能量沉积过高，可能导致靶材料的温度过高，影响靶的使用寿命和安全性；加速器的建设和运行成本也会大幅增加。经过综合分析和优化，发现1000MeV是一个较为合理、经济的入射质子能量。在这个能量下，既能保证产生足够的散裂中子，维持反应堆的能量增益，又能在一定程度上控制散裂靶的能量沉积和加速器的成本。对于质子入射位置对系统能量增益的影响，研究发现当质子入射位置在中心偏上位置时，系统的能量增益达到最大。这是因为在中心偏上位置入射时，质子能够与散裂靶充分相互作用，引发更多的散裂反应，产生更多的中子。质子在这个位置入射，其在散裂靶内的作用路径更长，能够与更多的靶核发生相互作用，从而提高中子产额。不同的入射位置会导致质子在散裂靶内的能量沉积分布不同，中心偏上位置入射时，能量沉积分布更加均匀，有利于热量的导出，提高反应堆的安全性和稳定性。而偏离中心偏上位置入射时，中子产额会逐渐降低，能量增益也会随之减小。在束流功率方面，它直接决定了加速器输入的能量，进而影响反应堆的输出功率。束流功率越大，加速器产生的质子数量越多，散裂反应产生的中子数量也相应增加，反应堆的输出功率就越高。束流功率的增加也会带来一些问题，如加速器的成本增加、散裂靶的辐照损伤加剧等。在设计TMSFEA时，需要在满足反应堆功率需求的前提下，合理选择束流功率。在初始束流强度小于最大束流强度4mA的限定条件下，选定系统总功率为300MWth。通过调整束流功率，使得反应堆能够在这个功率水平下稳定运行，同时兼顾加速器的成本和散裂靶的安全性。当束流功率过高时，散裂靶可能会因为辐照损伤而缩短使用寿命，增加维护成本；而束流功率过低，则无法满足反应堆的功率需求，影响能量输出。因此，在选择束流功率时，需要综合考虑反应堆的功率需求、加速器的成本以及散裂靶的安全性等多方面因素，找到一个最佳的平衡点。4.2.3TMSFEA反射层设计反射层在TMSFEA中起着至关重要的作用，它能够有效减少中子泄漏，提高中子的利用效率，进而对反应堆的性能产生显著影响。为了优化反射层的设计，对不同材料的反射层进行了详细的对比分析。液体铅作为一种常见的反射层材料，具有独特的性能特点。在中子能谱软化方面，液体铅反射层对堆芯中子能谱的软化作用较小。这意味着经过液体铅反射层反射回堆芯的中子，其能量分布更接近堆芯内中子的原始能量分布。这种特性有助于获得较大的转换比。在钍基次临界氯盐堆中，较大的转换比能够更有效地将钍转化为易裂变的铀-233，提高钍资源的利用效率。如果反射层对中子能谱的软化作用过大，会导致反射回堆芯的中子能量过低，无法有效地引发钍-232向铀-233的转化反应，从而降低转换比。液体铅还具有较高的密度和良好的中子反射性能，能够有效地反射中子，减少中子泄漏。其高密度使得中子在其中的散射概率增加，更多的中子能够被反射回堆芯，提高中子的利用效率。与液体铅反射层相比，其他材料的反射层在性能上存在一定的差异。例如，石墨反射层虽然在某些反应堆中也被广泛应用，但在TMSFEA中，其对中子能谱的软化作用相对较大。石墨的原子结构和中子散射特性导致反射回堆芯的中子能量有所降低，这可能会对转换比产生一定的负面影响。一些金属反射层，如铍反射层，虽然在中子反射性能上也有一定的优势，但铍具有毒性，在实际应用中需要考虑更多的安全防护措施，增加了系统的复杂性和成本。除了材料选择，反射层的厚度也是影响其性能的重要因素。随着反射层厚度的增加，中子泄漏逐渐减少。这是因为较厚的反射层能够提供更多的中子散射机会，使更多的中子被反射回堆芯。当反射层厚度增加到一定程度后，中子泄漏的减少趋势逐渐变缓。这是因为在达到一定厚度后，大部分能够被反射的中子已经被反射回堆芯，继续增加厚度对中子反射效果的提升作用有限。而且，过厚的反射层会增加反应堆的体积和重量，提高建设成本，同时也可能会对反应堆的热工水力性能产生一定的影响。在设计TMSFEA的反射层时，需要综合考虑反射层材料的性能和厚度，选择合适的反射层设计方案。通过优化反射层的设计，可以有效提高中子的利用效率，增加反应堆的能量增益和核燃料转换比，为TMSFEA的高效运行提供保障。4.2.4TMSFEA堆芯几何参数及熔盐配比堆芯几何参数和熔盐配比的优化对于TMSFEA的性能至关重要，它们直接影响着反应堆的中子学性能、热工水力性能以及运行稳定性。通过对不同堆芯尺寸和熔盐配比的遍历搜索，深入研究了它们对反应堆性能的影响。在堆芯尺寸方面，堆芯半径和高度是两个关键参数。当堆芯半径较小时，堆芯内的中子泄漏相对较大，因为中子更容易从较小的堆芯表面逃逸出去。这会导致反应堆内参与裂变反应的中子数量减少，影响能量增益和核燃料转换比。随着堆芯半径的增加，中子泄漏逐渐减少。较大的堆芯半径意味着中子在堆芯内有更多的机会与燃料发生相互作用，提高了中子的利用效率。堆芯半径过大也会带来一些问题。一方面，过大的堆芯半径会增加反应堆的体积和重量，提高建设成本。另一方面，堆芯半径过大可能会导致堆芯内的功率分布不均匀，局部区域的功率过高，对堆芯的热工水力性能和材料性能提出更高的要求。堆芯高度也对反应堆性能有重要影响。适当增加堆芯高度可以增加燃料的装载量，提高反应堆的输出功率。过高的堆芯高度可能会导致堆芯内的温度分布不均匀，增加冷却的难度，影响反应堆的安全运行。熔盐配比是另一个重要的参数。在三元氯盐NaCl-PuCl₃-ThCl₄中，PuCl₃和ThCl₄的比例对反应堆性能有显著影响。当PuCl₃的比例增加时，反应堆内的易裂变核素浓度增加，裂变反应更加剧烈，能量增益和核燃料转换比可能会提高。过高的PuCl₃比例也会带来一些问题。一方面，PuCl₃的增加会导致熔盐的放射性增强，对反应堆的防护和操作要求更高。另一方面，过高的PuCl₃比例可能会影响熔盐的物理性质，如熔点、黏度等，进而影响熔盐的流动性和传热性能。ThCl₄的比例增加则会增加钍的含量，有利于提高钍资源的利用效率。但如果ThCl₄比例过高，反应堆内的裂变反应可能会受到一定的抑制，因为钍-232需要吸收中子才能转化为易裂变的铀-233，过多的钍会稀释易裂变核素的浓度。通过对不同堆芯尺寸和熔盐配比的大量计算和分析，最终得到了TMSFEA的详细堆芯设计参数。在堆芯尺寸方面，确定了合适的堆芯半径和高度，使得中子泄漏和功率分布达到较好的平衡。在熔盐配比方面，找到了PuCl₃和ThCl₄的最佳比例，既能保证反应堆有较高的能量增益和核燃料转换比，又能兼顾熔盐的物理性质和放射性防护要求。这些优化后的堆芯几何参数和熔盐配比，为TMSFEA的高效稳定运行奠定了坚实的基础。4.2.5TMSFEA燃耗计算分析利用专门开发的次临界装置燃耗计算程序MOADS，对TMSFEA进行燃耗计算分析，深入研究反应堆在运行过程中的性能变化。在燃耗过程中，有效增殖因数（keff）是一个关键参数，它反映了反应堆内中子的平衡状态。对于TMSFEA这种次临界堆，keff始终小于1，需要依靠外源中子来维持链式裂变反应。随着燃耗的进行，keff会逐渐发生变化。在燃耗初期，由于堆芯内的燃料充足，中子吸收体相对较少，keff相对较高。随着燃耗的加深，燃料逐渐消耗，裂变产物等中子吸收体不断积累，keff会逐渐降低。当keff降低到一定程度时，反应堆的链式裂变反应将难以维持，需要采取相应的措施，如补充燃料或调整运行参数。核燃料转换比（CR）在燃耗过程中也会发生变化。在TMSFEA中，通过合理的设计，初始CR可以达到1.19。随着燃耗的进行，CR会逐渐下降。这是因为在燃耗过程中，易裂变核素不断消耗，而新产生的易裂变核素的速度逐渐减缓。在多循环演化过程中，通过离线燃料后处理等措施，可以在一定程度上维持CR的稳定，提高钍资源的利用效率。堆内主要核素的质量也会随着燃耗的进行而发生显著变化。钚和钍是TMSFEA中的重要核素。在燃耗过程中，钚不断参与裂变反应，其质量逐渐减少。到燃耗寿期结束时，有57%的钚被消耗。钍则通过吸收中子逐渐转化为易裂变的铀-233，其质量也会相应减少。在这个过程中，²³³U的裂变份额逐渐增加。在燃耗寿期结束时，²³³U的裂变份额可以达到70.9%，表明²³³U在反应堆的能量输出中起到了越来越重要的作用。这也说明TMSFEA在钍资源利用方面具有较高的效率，能够有效地将钍转化为易裂变核素，实现能量的持续输出。在整个燃耗寿期内，TMSFEA展现出了良好的性能。其寿期可以达到39年，在这期间能够稳定地输出能量。在无需燃料在线处理的条件下，TMSFEA能够高效利用钍资源，这为钍基核能系统的发展提供了有力的支持。通过燃耗计算分析，也为TMSFEA的运行管理和燃料循环方案的优化提供了重要依据。在实际运行中，可以根据燃耗计算的结果，合理调整反应堆的运行参数，优化燃料循环策略，进一步提高反应堆的性能和钍资源的利用效率。4.3本章小结本章围绕钍基次临界氯盐堆TMSFEA展开了全面且深入的设计与优化工作，取得了一系列关键成果，为该堆型的进一步发展和应用奠定了坚实基础。在设计原理方面，详细阐述了能量放大器概念以及TMSFEA的工作原理。TMSFEA利用加速器产生的高能质子与散裂靶发生散裂反应，产生外源中子驱动次临界熔盐堆，维持链式裂变反应。明确了能量增益G和核燃料转换比CR这两个重要参数的定义和计算公式，它们对于评估反应堆的性能至关重要。TMSFEA设定了明确的设计目标，在简化堆芯设计和燃料处理流程的前提下，实现钍资源的高效利用和稳定的能量输出。通过引入加速器提供的外中子源，补偿燃耗过程中裂变产物对中子的吸收，简化了燃料处理流程。在初始束流强度小于最大束流强度4mA的限定条件下，选定系统总功率为300MWth，初始转换比(CR)达到1.19，初始能量增益达到117，寿期可达39年。到燃耗寿期结束时，钚和钍的消耗比例以及²³³U的裂变份额表明TMSFEA在无需燃料在线处理的条件下，能够高效利用钍资源。在设计与优化过程中，对TMSFEA的各个关键部分进行了细致研究。在燃料盐选择上，确定了具有钚-钍溶解度高、熔点合适以及散裂中子产额较高的三元氯盐NaCl-PuCl₃-ThCl₄作为散裂靶和燃料盐。这种燃料盐能够满足TMSFEA对燃料性能的多方面要求，确保了反应堆的稳定运行和高效性能。在加速器质子能量与束流功率方面，通过研究发现1000MeV是更合理、更经济的入射质子能量，且当质子入射位置在中心偏上位置时，系统的能量增益达到最大。在束流功率的选择上，综合考虑了反应堆的功率需求、加速器的成本以及散裂靶的安全性等多方面因素，在初始束流强度小于最大束流强度4mA的限定条件下，选定系统总功率为300MWth。在反射层设计中，对比分析了不同材料的反射层，发现液体铅反射层对堆芯中子能谱的软化作用较小，有助于获得较大的转换比。通过对反射层厚度的研究，确定了合适的反射层厚度，有效减少了中子泄漏，提高了中子的利用效率。在堆芯几何参数及熔盐配比方面，通过对不同堆芯尺寸和熔盐配比的遍历搜索，确定了最佳的堆芯半径、高度以及PuCl₃和ThCl₄的比例。这些参数的优化使得中子泄漏和功率分布达到较好的平衡，兼顾了熔盐的物理性质和放射性防护要求，为反应堆的高效稳定运行提供了保障。利用MOADS程序对TMSFEA进行燃耗计算分析，研究了燃耗过程中有效增殖因数（keff）、核燃料转换比（CR）以及堆内主要核素质量的变化。结果表明，TMSFEA在整个燃耗寿期内能够稳定输出能量，且在钍资源利用方面具有较高的效率。综上所述，TMSFEA的设计方案在理论上具有较高的可行性和优势。通过优化设计，在钍资源利用效率、能量输出稳定性以及系统安全性等方面都展现出良好的性能。然而，该设计方案仍存在一些需要进一步改进的地方。在燃料盐的长期稳定性和腐蚀性方面，虽然目前的三元氯盐NaCl-PuCl₃-ThCl₄具有良好的性能，但在长期运行过程中，其稳定性和对反应堆结构材料的腐蚀性仍需进一步研究和验证。在加速器的可靠性和维护成本方面，作为TMSFEA的关键组成部分，加速器的高可靠性和低维护成本对于反应堆的长期稳定运行至关重要，需要进一步优化加