CN117594267A - 堆芯中子学模型的修正方法、装置、设备和存储介质 - Google Patents
堆芯中子学模型的修正方法、装置、设备和存储介质 Download PDFInfo
- Publication number
- CN117594267A CN117594267A CN202311458625.8A CN202311458625A CN117594267A CN 117594267 A CN117594267 A CN 117594267A CN 202311458625 A CN202311458625 A CN 202311458625A CN 117594267 A CN117594267 A CN 117594267A
- Authority
- CN
- China
- Prior art keywords
- neutron flux
- measured
- neutron
- section
- group
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 77
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 230000004907 flux Effects 0.000 claims abstract description 334
- 238000004590 computer program Methods 0.000 claims description 23
- 238000013178 mathematical model Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 abstract description 9
- 238000009826 distribution Methods 0.000 description 23
- 230000000694 effects Effects 0.000 description 13
- 238000012544 monitoring process Methods 0.000 description 11
- 238000010521 absorption reaction Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 9
- 239000000446 fuel Substances 0.000 description 7
- 239000003758 nuclear fuel Substances 0.000 description 7
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 5
- 229910052796 boron Inorganic materials 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 229910052724 xenon Inorganic materials 0.000 description 5
- FHNFHKCVQCLJFQ-UHFFFAOYSA-N xenon atom Chemical compound [Xe] FHNFHKCVQCLJFQ-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004992 fission Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 239000003607 modifier Substances 0.000 description 3
- 229910052778 Plutonium Inorganic materials 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- OYEHPCDNVJXUIW-UHFFFAOYSA-N plutonium atom Chemical compound [Pu] OYEHPCDNVJXUIW-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 229910052770 Uranium Inorganic materials 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
- G21C17/108—Measuring reactor flux
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21D—NUCLEAR POWER PLANT
- G21D3/00—Control of nuclear power plant
- G21D3/001—Computer implemented control
- G21D3/002—Core design; core simulations; core optimisation
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21D—NUCLEAR POWER PLANT
- G21D3/00—Control of nuclear power plant
- G21D3/08—Regulation of any parameters in the plant
- G21D3/10—Regulation of any parameters in the plant by a combination of a variable derived from neutron flux with other controlling variables, e.g. derived from temperature, cooling flow, pressure
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本申请涉及一种堆芯中子学模型的修正方法、装置、设备和存储介质。方法包括:根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;根据各节块的状态参数确定各节块的等效截面参数;根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。本申请实施例,可以提高堆芯中子学模型修正的准确度,以便于基于修正后的模型可以更准确地计算理论功率,有利于更加准确地监测堆芯内的核反应情况。
Description
技术领域
本申请涉及核反应堆技术领域,特别是涉及一种堆芯中子学模型的修正方法、装置、设备和存储介质。
背景技术
核反应在核电站中扮演着重要的角色,堆芯是核反应的核心。堆芯在线监测是保障核电站安全运行的重要手段之一。通常情况下,堆芯监测系统通过将基于中子学模型计算的理论功率分布与堆芯检测仪器测量的功率分布结合的方式,监测核反应堆芯内的核反应情况。若理论功率分布与测量功率分布之间的偏差超过预设阈值,则表明核反应的运行可能存在风险。因此,如何获取更加准确的理论功率分布对于核电站的监测十分重要。
相关技术中,通过中子等效截面偏差的全局优化算法对中子学模型进行调整,以便于基于调整后的中子学模型可以计算理论功率分布。但是,相关技术中模型调整方式的准确度不太高。
发明内容
基于此,有必要针对上述技术问题,提供一种能够提高中子学模型修正准确度的堆芯中子学模型的修正方法、装置、设备和存储介质。
第一方面,本申请提供了一种堆芯中子学模型的修正方法,包括:
根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;
根据各节块的状态参数确定各节块的等效截面参数;
根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。
在其中一个实施例中,根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型,包括:
对于各节块,根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定第一泄漏曲率参数;
根据第一泄漏曲率参数,确定修正后的堆芯中子学模型。
在其中一个实施例中,根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定第一泄漏曲率参数,包括:
根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和快群泄漏修正量,确定快群中子平衡方程,其中,快群泄漏修正量中包括第一泄漏曲率参数;
根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和热群泄漏修正量,确定热群中子平衡方程,其中,热群泄漏修正量中包括第一泄漏曲率参数;
根据快群中子平衡方程和热群中子平衡方程,确定第一泄漏曲率参数。
在其中一个实施例中,根据第一泄漏曲率参数,确定修正后的堆芯中子学模型,包括:
根据第一泄漏曲率参数分别对快群中子平衡方程和热群中子平衡方程进行修正,得到修正后的堆芯中子学模型。
在其中一个实施例中,根据各节块的状态参数确定各节块的等效截面参数,包括:
对于各节块,根据预设状态参数集合确定节块对应的多个状态参数;
根据多个状态参数按照预设差值算法,确定节块的等效截面参数。
在其中一个实施例中,目标堆芯的测量功率包括目标堆芯的各节块的测量功率,根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量,包括:
对于各节块,根据修正前的堆芯中子学模型确定节块的理论快群中子通量、理论热群中子通量、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面;
根据节块的理论快群中子通量、节块的理论热群中子通量、节块的测量功率、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面,确定节块的测量快群中子通量和测量热群中子通量。
在其中一个实施例中,上述方法还包括:
根据扰动堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率;
根据各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率,确定第二泄漏曲率参数;
根据第二泄漏曲率参数,确定修正后的扰动堆芯中子学模型;
基于各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定修正后的扰动堆芯中子学模型是否正确。
在其中一个实施例中,基于各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定修正后的扰动堆芯中子学模型是否正确,包括:
根据修正后的扰动堆芯中子学模型,确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率;
根据各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及各扰动节块的测量功率与对应的修正功率之间的差值,确定修正后的扰动堆芯中子学模型是否正确。
在其中一个实施例中,根据各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及各扰动节块的测量功率与对应的修正功率之间的差值,确定修正后的扰动堆芯中子学模型是否正确,包括:
若各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值小于第一预设阈值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值小于第二预设阈值,且各扰动节块的测量功率与对应的修正功率之间的差值小于第三预设阈值,则确定修正后的扰动堆芯中子模型正确。
在其中一个实施例中,上述方法还包括:
根据修正后的堆芯中子学模型确定目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率。
第二方面,本申请还提供了一种堆芯中子学模型的修正装置,包括:
第一确定模块,用于根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;
第二确定模块,用于根据各节块的状态参数确定各节块的等效截面参数;
第三确定模块,用于根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。
第三方面,本申请还提供了一种计算机设备,该计算机设备包括存储器和处理器,存储器存储有计算机程序,处理器执行计算机程序时以实现上述第一方面的方法的步骤。
第四方面,本申请还提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现上述第一方面的方法的步骤。
第五方面,本申请还提供了一种计算机程序产品,包括计算机程序,该计算机程序被处理器执行时实现上述第一方面的方法的步骤。
上述堆芯中子学模型的修正方法、装置、设备和存储介质,通过根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;根据各节块的状态参数确定各节块的等效截面参数;根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。相对于传统技术中通过中子等效截面偏差的全局优化算法对中子学模型进行调整的方式而言,本申请实施例根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定各节块的测量快群中子通量和热群中子通量,通过将实际测量数据与理论数据相结合的方式,可以准确地确定各个节块在快群中子通量和热群中子通量的分布情况。进一步地,根据各节块的状态参数确定各节块的等效截面参数的方式,由于考虑不同状态参数下的等效截面参数是不同的,因此,基于状态参数可以得到更准确的等效截面参数。进一步地,通过根据准确的各节块的测量快群中子通量、测量热群中子通量和等效截面参数,可以更加准确地确定修正后的堆芯中子学模型,从而提高了模型修正的准确度,以便于基于修正后的模型可以更准确地计算理论功率,有利于更加准确地监测堆芯内的核反应情况。
附图说明
为了更清楚地说明本申请实施例或相关技术中的技术方案,下面将对实施例或相关技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本申请一个实施例中堆芯中子学模型的修正方法的实施环境示意图;
图2为本申请一个实施例中堆芯中子学模型的修正方法的流程示意图;
图3为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图4为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图5为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图6为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图7为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图8为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图;
图9为本申请一个实施例中堆芯中子学模型的修正装置的结构示意图;
图10为本申请一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
图1为本申请一个实施例中堆芯中子学模型的修正方法的实施环境示意图,如图1所示,本申请实施例的实施环境可以包括终端10和服务器11;其中,终端10通过网络可以与服务器11进行通信。其中,终端10可以包括但不限于是各种个人计算机、笔记本电脑、智能手机和平板电脑。服务器11可以用独立的服务器或者是多个服务器组成的服务器集群来实现。数据存储系统可以集成在服务器11上,也可以放在云上或其他网络服务器上。
结合图1所示的实施环境,本申请实施例中,服务器11可以根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;并根据各节块的状态参数确定各节块的等效截面参数。进一步地,服务器11可以根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。进一步地,服务器11可以根据修正后的堆芯中子学模型确定目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率,并将目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率发送给终端10。
在实际应用中,通过对测量功率的监测和分析,可以及时发现堆芯内部的异常情况,并采取相应的措施,以确保核电站的稳定运行。当理论功率分布与测量功率分布的偏差超过预设阈值时,需要采取进一步的措施。一种方法是通过降低功率水平来减少理论-测量偏差大带来的影响,但该方法会影响核电站的经济性。另一种方法是调整中子学模型,进而调整理论功率分布,使理论-测量偏差逐渐缩小。相关技术中,通过中子等效截面偏差的全局优化算法对中子学模型进行调整,以便于基于调整后的中子学模型可以计算理论功率分布,但是相关技术中的方法需要进行大量的计算来生成中子等效截面,计算复杂度较高,需要更多的计算资源和时间。进一步地,该方法的搜索结果依赖于初始参数的设置和调整,因此需要进行一定的人工干预,会导致结果的不确定性。因此,这种中子学模型参数调整的方式可能会导致中子学模型修正的准确性较低。
本申请实施例提供的堆芯中子学模型的修正方法、装置、设备和存储介质,通过根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;根据各节块的状态参数确定各节块的等效截面参数;根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。相对于传统技术中通过中子等效截面偏差的全局优化算法对中子学模型进行调整的方式而言,本申请实施例根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定各节块的测量快群中子通量和热群中子通量,通过将实际测量数据与理论数据相结合的方式,可以准确地确定各个节块在快群中子通量和热群中子通量的分布情况。进一步地,根据各节块的状态参数确定各节块的等效截面参数的方式,由于考虑不同状态参数下的等效截面参数是不同的,因此,基于状态参数可以得到更准确的等效截面参数。进一步地,通过根据准确的各节块的测量快群中子通量、测量热群中子通量和等效截面参数,不需要预先设置初始条件和初始参数,从而减少了人为引入的不确定性,可以更加准确地确定修正后的堆芯中子学模型,从而提高了模型修正的准确度,以便于基于修正后的模型可以更准确地计算理论功率,有利于更加准确地监测堆芯内的核反应情况。
在一个实施例中,图2为本申请一个实施例中堆芯中子学模型的修正方法的流程示意图,本申请实施例中以该方法应用于服务器11为例进行说明。如图2所示,本申请实施例的方法可以包括以下步骤:
步骤S201,根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量。
示例性地,本申请实施例中涉及的堆芯中子学模型是指用于描述核反应堆中中子输运和相互作用的数学模型;修正前的堆芯中子学模型即未进行修正之前的模型。
示例性地,本申请实施例中涉及的目标堆芯的测量功率是指对目标反应堆芯(或者简称为目标堆芯)实际测量得到的功率值。在核反应堆中,功率是衡量反应堆能量输出的指标,可以通过测量热量释放或使用其他监测设备来测量。其中,监测设备可以是堆芯功率分布的在线监测系统,在线监测系统通过安置多个固定式中子探测器和信号采集系统在核电站内部进行数据采集和处理,得到堆芯的测量功率。当然,目标堆芯的测量功率还可以通过其他方式测量。
示例性地,本申请实施例中涉及的测量快群中子通量是指根据在线监测系统测量得到的快中子能区域内的单位时间内通过单位面积的中子数目。本申请实施例中涉及的各节块的测量热群中子通量是指在线监测系统测量得到的热中子能区域内的单位时间内通过单位面积的中子数目。
示例性地,本申请实施例中可以将目标堆芯划分为多个节块。
本申请实施例中,服务器可以根据修正前的堆芯中子学模型和监测得到的目标堆芯测量功率,确定目标堆芯中各节块的测量快群中子通量和测量热群中子通量。本申请实施例中,通过将实际测量数据与理论数据相结合的方式,可以准确确定各个节块在快群中子通量和热群中子通量的分布情况。
示例性地,服务器可以根据修正前的堆芯中子学模型确定各节块的理论参数,并根据各节块的理论参数和各节块的测量功率确定各节块的测量快群中子通量和测量热群中子通量。
步骤S202,根据各节块的状态参数确定各节块的等效截面参数。
示例性地,本申请实施例中涉及的状态参数可以包括但不限于以下至少一项:燃耗、硼浓度、慢化剂密度、燃料有效温度、氙浓度和Pu/U比。其中,燃耗是指核燃料在反应堆中燃烧的程度。硼是一种常用的中子吸收材料,在核反应堆中用于控制中子通量,通过控制堆芯中的硼浓度,可以调节中子吸收,从而影响到中子通量的分布情况。慢化剂是用于减慢中子速度的材料,慢化剂的密度会影响中子的慢化过程,进而影响到中子通量的分布情况。燃料有效温度是指核燃料的温度,它会影响到燃料中的原子核振动状态以及燃料材料的物理性质。氙是一种短寿命的放射性核素,是核燃料燃耗后的裂变产物之一,氙的存在会影响中子的吸收和散射行为,从而对中子通量的分布产生影响。Pu/U比是指核燃料中钚与铀的比例,钚是一种重要的核燃料裂变产物。Pu/U比的变化会影响燃料材料的截面积和吸收性质,进而对中子通量的分布产生影响。
本步骤中,服务器可以获取各节块的状态参数,并根据各节块中不同的状态参数,来确定各节块的等效截面参数。
示例性地,本申请实施例中涉及的等效截面参数可以包括但不限于以下至少一项:快群扩散系数、热群扩散系数、向下宏观散射截面、快群宏观吸收截面、热群宏观吸收截面、快群宏观中子产生截面、热群宏观中子产生截面、快群宏观能量产生截面和热群宏观能量产生截面。
本申请实施例中,通过根据各节块的状态参数确定各节块的等效截面参数的方式,由于考虑不同状态参数下的等效截面参数是不同的,因此,基于状态参数可以得到更准确的等效截面参数。
步骤S203,根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。
本步骤中,服务器可以根据步骤S201得到的各节块的测量快群中子通量、测量热群中子通量和步骤S202得到的等效截面参数,确定修正后的堆芯中子学模型。
示例性地,本申请实施例中涉及的修正后的堆芯中子学模型可以是指对修正前的堆芯中子学模型进行修正之后得到的堆芯中子学模型。
本申请实施例中,通过根据准确的各节块的测量快群中子通量、测量热群中子通量和等效截面参数,从而可以更加准确地确定修正后的堆芯中子学模型。
上述堆芯中子学模型的修正方法中,通过根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量;根据各节块的状态参数确定各节块的等效截面参数;根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。可见,本申请实施例中,根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定各节块的测量快群中子通量和热群中子通量,通过将实际测量数据与理论数据相结合的方式,可以准确地确定各个节块在快群中子通量和热群中子通量的分布情况。进一步地,根据各节块的状态参数确定各节块的等效截面参数的方式,由于考虑不同状态参数下的等效截面参数是不同的,因此,基于状态参数可以得到更准确的等效截面参数。进一步地,通过根据准确的各节块的测量快群中子通量、测量热群中子通量和等效截面参数,可以更加准确地确定修正后的堆芯中子学模型,从而提高了模型修正的准确度,以便于基于修正后的模型可以更准确地计算理论功率,有利于更加准确地监测堆芯内的核反应情况。
在一个实施例中,在上述实施例的基础上,本申请实施例中对上述实施例中涉及的步骤S203中根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型的相关内容作示例性的介绍说明。如图3所示,上述步骤S203可以包括步骤S2031以及步骤S2032。
步骤S2031,对于各节块,根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定第一泄漏曲率参数。
示例性地,本申请实施例中涉及的第一泄漏曲率参数用于指示快群中子和热群中子在节块之间的泄漏程度的参数。
本步骤中,服务器可以根据测量快群中子通量和热群中子通量、等效截面参数以及相应的泄漏修正量,确定快群中子平衡方程和热群中子平衡方程,从而确定第一泄漏曲率参数。
可选地,服务器可以根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和快群泄漏修正量,确定快群中子平衡方程,其中,快群泄漏修正量中包括第一泄漏曲率参数。
示例性地,本申请实施例中涉及的快群泄漏修正量可以表示为其中,为第n个节块的第一泄漏曲率参数;Dn,1为第n个节块的快群扩散系数的导数;/>为第n个节块的测量快群中子通量;n为大于0的整数。
当然,快群泄漏修正量还可以表示为上述公式的其他变形或等效公式。
示例性地,本申请实施例中涉及的快群中子平衡方程可以表示为:
其中,为第n个节块表面的快群中子通量;Σa,n,1为第n个节块的快群宏观吸收截面;Σr,n为第n个节块的向下宏观散射截面;vΣf,n,g为第n个节块的快群宏观中子产生截面;keff为中子系数;/>为第n个节块的测量热群中子通量。进一步地,/>表示由节块六个表面净流出的快群中子通量;/>表示由吸收或转移(由快群到热群慢化)而消失的快群中子通量;/>表示由于裂变产生的快群中子通量。
当然,快群中子平衡方程还可以表示为上述公式(1)的其他变形或等效公式。
进一步地,服务器可以根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和热群泄漏修正量,确定热群中子平衡方程,其中,热群泄漏修正量中包括第一泄漏曲率参数。
示例性地,本申请实施例中涉及的热群泄漏修正量可以表示为其中,Dn,2为第n个节块的热群扩散系数的导数。
示例性地,本申请实施例中涉及的热群中子平衡方程可以表示为:
其中,为第n个节块表面的热群中子通量;Σa,n,2为第n个节块的热群宏观吸收截面。进一步地,/>表示由节块六个表面净流出的热群中子通量;/>表示由吸收而消失的热群中子通量。
当然,热群中子平衡方程还可以表示为上述公式(2)的其他变形或等效公式。
进一步地,服务器可以根据快群中子平衡方程和热群中子平衡方程,确定第一泄漏曲率参数。
示例性地,服务器可以将快群中子平衡方程(1)和热群中子平衡方程(2)转换成一个矩阵求解的形式,可以表示为:
AnΦn=0 (3)
其中,Φn为相应的快群热群中子通量;An为根据快群中子平衡方程(1)和热群中子平衡方程(2)的所有系数得到的矩阵。
若矩阵An的行列式为零,则公式(3)存在一个解,即第n个节块上的快群中子通量和热群中子通量。因此,使det(An)=0,可以表示为:
需要说明的是,公式(4)有两个解,但只有一个具有物理意义,另一个解会导致负的群中子通量分布。因此,通过使矩阵An的行列式为零并求解矩阵的秩可以得到第一泄漏曲率参数Bn 2的值。
当然,公式(4)还可以表示为上述公式(4)的其他变形或等效公式。
可见,本申请实施例中,服务器可以通过增加快群泄漏修正量和热群泄漏修正量的方式,确定快群中子平衡方程和热群中子平衡方程,以便于可以更加准确地确定第一泄漏曲率参数。
步骤S2032,根据第一泄漏曲率参数,确定修正后的堆芯中子学模型。
本步骤中,服务器可以根据第一泄漏曲率参数对中子平衡方程进行修正,从而得到修正后的堆芯中子学模型。
示例性地,服务器可以根据第一泄漏曲率参数分别对快群中子平衡方程和热群中子平衡方程进行修正,得到修正后的堆芯中子学模型。
一种可能的实现方式中,如果修正前的堆芯中子学模型为不带泄漏修正量,服务器可以根据第一泄漏曲率参数确定快群泄漏修正量和热群泄漏修正量。进一步地,服务器可以通过在修正前的快群中子平衡方程中增加快群泄漏修正量得到修正后的快群中子平衡方程,以及通过在修改前的热群中子平衡方程中增加热群泄漏修正量得到修正后的热群中子平衡方程,从而得到了修正后的堆芯中子学模型。
另一种可能的实现方式中,如果修正前的堆芯中子学模型已经有泄漏修正量,服务器可以根据第一泄漏曲率参数分别对修正前的快群中子平衡方程中的快群泄漏修正量进行修正,以及对修正前的热群中子平衡方程中的热群泄漏修正量进行修正,从而得到修正后的堆芯中子学模型。
综上,对于所述各节块,服务器可以根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定第一泄漏曲率参数。进一步地,服务器可以根据第一泄漏曲率参数确定修正后的堆芯中子学模型。由于本申请实施例中的测量快群中子通量、测量热群中子通量和等效截面参数比较准确,因此,本申请实施例中服务器根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,可以准确地确定第一泄漏曲率参数,从而根据准确的第一泄漏曲率参数可以得到准确的修正后的堆芯中子学模型,以便于基于修正后的模型可以更准确地计算理论功率。
在一个实施例中,图4为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图,如图4所示,在上述实施例的基础上,本申请实施例中对上述实施例中涉及的步骤S202中根根据各节块的状态参数确定各节块的等效截面参数作示例性的介绍说明。如图4所示,上述步骤S202包括步骤S2021以及步骤S2022。
步骤S2021,对于各节块,根据预设状态参数集合确定节块对应的多个状态参数。
示例性地,本申请实施例中涉及的预设状态参数集合是指在服务器中预先生成的插值表。该状态参数集合中可以覆盖不同截面参数所对应的所有状态参数范围,可以包括但不限于燃耗、硼浓度、慢化剂密度、燃料有效温度、氙浓度和Pu/U比等参数范围。
本步骤中,对于各节块,考虑到等效截面参数包括一个或多个截面参数,服务器可以根据预设状态参数集合确定各截面参数分别对应的一个或多个状态参数。
例如,假设预设状态参数集合中包括截面参数A对应的状态参数A1和状态参数A2,以及截面参数B对应的状态参数B1、状态参数B2和状态参数B3,则服务器可以根据预设状态参数集合确定截面参数A对应的状态参数A1和状态参数A2,以及截面参数B对应的状态参数B1、状态参数B2和状态参数B3。
步骤S2022,根据多个状态参数按照预设插值算法,确定节块的等效截面参数。
示例性地,本申请实施例中涉及的预设的插值算法可以包括但不限于:线性插值、多项式插值、或者样条插值。其中,线性插值是假设给定的数据点之间的函数为线性关系。多项式插值是指假设给定的数据点之间的函数可以表示为一个多项式。样条插值是指假设给定的数据点之间的函数可以通过多个插值段拼接而成。
本步骤中,对于各节块的各截面参数,服务器在根据预设状态参数集合中确定该截面参数对应的一个或多个状态参数之后,可以按照预设的插值算法确定该节块的该截面参数。
例如,对于快群宏观吸收截面可能与预设状态参数集合中的燃耗、硼浓度、慢化剂密度相关;对于热群宏观能量产生截面可能与预设状态参数集合中的燃料有效温度、氙浓度和Pu/U相关。
可选的,对于各节块的各截面参数,服务器根据该截面参数对应的一个或多个状态参数,可以通过如下多项式插值公式(5)确定该节块的该截面参数。
其中,SP1、SP2和SP3为节块的三种状态参数;ci,j为展开系数;Σ为各节块的截面参数信息。
当然,服务器根据该截面参数对应的一个或多个状态参数,还可以通过上述多项式插值公式(5)的其他变形或等效公式确定该节块的该截面参数。
综上,对于各节块,服务器可以在预设状态参数集合中选择并确定不同截面对应的多个状态参数,并按照预设的插值算法,从而确定节块的等效截面参数。本申请实施例中通过使用预设的状态参数集合,可以选择适合不同截面的不同状态参数,这样可以考虑到不同的工况条件对截面参数的影响,以更准确地确定节块的等效截面参数。进一步地,服务器通过按照预设的插值算法可以准确地得到在实际运行中节块的等效截面参数,从而有利于提高模型修正的准确性。
在一个实施例中,图5为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图,如图5所示,在上述实施例的基础上,本申请实施例中对上述实施例中涉及的步骤S201中根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量作示例性的介绍说明。如图5所示,上述步骤S201包括步骤S2011以及步骤S2012。
步骤S2011,对于各节块,根据修正前的堆芯中子学模型确定节块的理论快群中子通量、理论热群中子通量、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面。
示例性地,本申请实施例中涉及的理论快群中子通量指根据修正前的堆芯中子学模型计算得到的快中子能区域内的单位时间内通过单位面积的中子数目;理论热群中子通量指根据修正前的堆芯中子学模型计算得到的热中子能区域内的单位时间内通过单位面积的中子数目。
示例性地,本申请实施例中涉及的快群宏观能量产生截面是指在高能区域内,中子与核反应堆燃料中的原子核相互作用并产生能量的过程;热群宏观能量产生截面是指在较低能区域内,中子与核反应堆燃料中的原子核相互作用并产生能量的过程。
本步骤中,对于各节块,服务器可以根据修正前的堆芯中子学模型计算得到节块的理论快群中子通量、理论热群中子通量、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面。
步骤S2012,根据节块的理论快群中子通量、节块的理论热群中子通量、节块的测量功率、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面,确定节块的测量快群中子通量和测量热群中子通量。
本步骤中,对于各节块,服务器可以根据节块的理论快群中子通量、理论热群中子通量、测量功率、快群宏观能量产生截面和热群宏观能量产生截面,从而确定节块的测量快群中子通量和测量热群中子通量。
示例性地,本申请实施例中假设修正前的堆芯中子学模型能够有效预测局部的中子能谱,即理论快群中子通量、理论热群中子通量和测量快群中子通量、测量热群中子通量之间的关系可以表示为:
其中,为理论快群中子通量;/>为理论热群中子通量。
当然,公式(6)还可以表示为上述公式(6)的其他变形或等效公式。
根据节块的测量功率,可以得到:
其中,κΣf,n,1为第n个节块的快群宏观能量产生截面;κΣf,n,2为第n个节块的热群宏观能量产生截面;为第n个节块的测量功率。
当然,公式(7)还可以表示为上述公式(7)的其他变形或等效公式。
对于各节块,服务器根据公式(6)和公式(7)可以计算得到节块的测量快群中子通量和测量热群中子通量。
综上,对于各节块,服务器可以根据修正前的堆芯中子学模型计算得到节块的理论快群中子通量、理论热群中子通量、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面。进一步地,服务器可以根据节块的理论快群中子通量、节块的理论热群中子通量、节块的测量功率、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面,从而确定节块的测量快群中子通量和测量热群中子通量。本申请实施例中,通过修正前的堆芯中子学模型可以准确地计算得到理论中子通量,并通过将理论中子通量与实际测量功率相结合,可以更准确地计算得到测量中子通量,本申请仅依赖于测量功率,而不需要其他高质量数据的支持,从而有利于提高中子学模型的准确性。
在一个实施例中,图6为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图,如图6所示,在上述实施例的基础上,本申请实施例中对上述堆芯中子学模型的修正方法的模型修正进行验证的相关内容作示例性的介绍说明。如图6所示,本申请实施例的方法还可以包括步骤S601、步骤S602、步骤S603以及步骤S604。
步骤S601,根据扰动堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率。
本步骤中,服务器可以通过预先建立扰动的堆芯中子学模型,并根据扰动的堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率。
示例性地,本申请实施例中涉及的扰动堆芯中子学模型用于指示基于对堆芯设计和运行参数进行扰动分析所建立的模型,通过建立扰动模型,可以模拟和评估在实际运行中可能发生的变化。各扰动节块的等效截面参数是指根据扰动堆芯中子学模型确定的等效截面参数;各扰动节块的测量快群中子通量是指根据扰动堆芯中子学模型确定的测量快群中子通量;各扰动节块的测量热群中子通量是指根据扰动堆芯中子学模型确定的测量热群中子通量;各扰动节块的测量功率是指根据扰动堆芯中子学模型确定的测量功率。
步骤S602,根据各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率,确定第二泄漏曲率参数。
本步骤中,服务器可以根据扰动的堆芯中子学模型得到的各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率,确定扰动的快群中子平衡方程和扰动的热群中子平衡方程。进一步地,服务器可以根据扰动的快群中子平衡方程和扰动的热群中子平衡方程,确定第二泄漏曲率参数。
应理解,服务器确定第二泄漏曲率参数的方式,可以参考上述确定第一泄漏曲率参数的方式,此处不再赘述。
示例性地,本申请实施例中涉及的第二泄漏曲率参数用于指示在扰动的堆芯中子学模型中,快群中子和热群中子在节块之间的泄漏程度的参数。
步骤S603,根据第二泄漏曲率参数,确定修正后的扰动堆芯中子学模型。
示例性地,本申请实施例中涉及的修正后的扰动堆芯中子学模型用于指示修正之后的扰动堆芯中子学模型。
本步骤中,服务器可以根据第二泄漏曲率参数,分别对扰动的快群中子平衡方程中的快群泄漏修正量和扰动的热群中子平衡方程中的热群泄漏修正量进行修正,从而得到修正后的扰动的堆芯中子学模型。
步骤S604,基于各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定修正后的扰动堆芯中子学模型是否正确。
本步骤中,服务器可以基于得到的修正后的扰动堆芯中子学模型,计算得到各扰动节块的修正快群中子通量、修正热群中子通量和修正功率,从而确定修正后的扰动堆芯中子学模型是否正确。
综上,本申请实施例中,服务器可以通过预先建立扰动的堆芯中子学模型,根据扰动的堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率。服务器可以确定扰动的快群中子平衡方程和扰动的热群中子平衡方程,从而确定第二泄漏曲率参数。进一步地,服务器可以根据第二泄漏曲率参数,分别对扰动的快群中子平衡方程和扰动的热群中子平衡方程进行修正,得到修正后的扰动的堆芯中子学模型。并基于得到的修正后的扰动堆芯中子学模型,计算得到各扰动节块的修正快群中子通量、修正热群中子通量和修正功率,从而确定修正后的扰动堆芯中子学模型是否正确。可见,本申请实施例通过建立扰动的堆芯中子学模型可以提供堆芯的响应和性能的预测和评估;通过确定扰动的快群中子平衡方程和扰动的热群中子平衡方程,可以更加准确地得到快群中子通量和热群中子通量。通过修正扰动的中子平衡方程,可以准确地得到修正后的扰动的堆芯中子学模型,从而提高了模型修正的准确度。进一步地,根据修正后的扰动中子学模型可以更准确地计算理论功率,更准确地反映堆芯内的实际情况和测量数据。
在一个实施例中,图7为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图,如图7所示,在上述实施例的基础上,本申请实施例中对上述实施例中涉及的步骤S604中基于各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定修正后的扰动堆芯中子学模型是否正确作示例性的介绍说明。如图7所示,上述步骤S604可以包括步骤S6041以及步骤S6042。
步骤S6041,根据修正后的扰动堆芯中子学模型,确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率。
示例性地,本申请实施例中涉及的各扰动节块的修正快群中子通量用于指示根据修正后的扰动堆芯中子学模型计算得到的快群中子通量;各扰动节块的修正热群中子通量用于指示根据修正后的扰动堆芯中子学模型计算得到的热群中子通量;各扰动节块的修正功率用于指示根据修正后的扰动堆芯中子学模型计算得到的修正功率。
本步骤中,服务器可以根据修正后的扰动堆芯中子学模型,从而确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率。
S6042,服务器根据各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及各扰动节块的测量功率与对应的修正功率之间的差值,确定修正后的扰动堆芯中子学模型是否正确。
一种可能的实现方式中,若各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值小于第一预设阈值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值小于第二预设阈值,且各扰动节块的测量功率与对应的修正功率之间的差值小于第三预设阈值,则确定修正后的扰动堆芯中子模型正确。
示例性地,本申请实施例中涉及的第一预设阈值可以表示各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值的允许范围。如果差值小于第一预设阈值,则可以认为修正后的模型在快群中子通量方面是正确的。
示例性地,本申请实施例中涉及的第二预设阈值可以表示各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值的允许范围。如果差值小于第二预设阈值,则可以认为修正后的模型在热群中子通量方面是正确的。
示例性地,本申请实施例中涉及的第三预设阈值可以表示各扰动节块的测量功率与对应的修正功率之间的差值的允许范围。如果差值小于第三预设阈值,则可以认为修正后的模型在功率方面是正确的。
另一种可能的实现方式中,若各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值不小于第一预设阈值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值不小于第二预设阈值,或者各扰动节块的测量功率与对应的修正功率之间的差值不小于第三预设阈值,则确定修正后的扰动堆芯中子学模型错误。
示例性地,各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值小于第一预设阈值,可以表示为:
其中,为各扰动节块的修正快群中子通量;/>为各扰动节块的测量快群中子通量;/>为第一预设阈值。
应理解,针对每个扰动节块,从测量快群中子通量和修正快群中子通量的差异中选择最大值作为扰动节块的测量快群中子通量与修正快群中子通量之间的差值。
当然,公式(8)还可以表示为上述公式(8)的其他变形或等效公式。
示例性地,各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值小于第二预设阈值,可以表示为:
其中,为各扰动节块的修正热群中子通量;/>为各扰动节块的测量热群中子通量;/>为第二预设阈值。
应理解,针对每个扰动节块,从测量热群中子通量和修正热群中子通量的差异中选择最大值作为扰动节块的测量热群中子通量与修正热群中子通量之间的差值。
当然,公式(9)还可以表示为上述公式(9)的其他变形或等效公式。
示例性地,各扰动节块的测量功率与对应的修正功率之间的差值小于第三预设阈值,可以表示为:
max(|PR-PM|)<εP (10)
其中,PR为各扰动节块的修正功率;PM为各扰动节块的测量功率;εP为第三预设阈值。
应理解,针对每个扰动节块,从扰动的测量功率与对应的修正功率之间的差值之中选择最大值作为各扰动节块的测量功率与对应的修正功率之间的差值。
当然,公式(10)还可以表示为上述公式(10)的其他变形或等效公式。
综上,服务器可以根据修正后的扰动堆芯中子学模型,从而确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率。进一步地,服务器可以根据各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及各扰动节块的测量功率与对应的修正功率之间的差值,确定修正后的扰动堆芯中子学模型是否正确。可见,本申请实施例中,通过使用修正后的扰动堆芯中子学模型,服务器可以准确地计算出每个扰动节块的修正快群中子通量和修正热群中子通量,修正后的中子通量将更准确地反映实际堆芯中的中子分布情况。进一步地,通过和预设的三个阈值进行对比,若差值均小于各自的预设阈值,则可以认为修正后的模型是正确的,从而有利于提高模型修正的准确性。
在一个实施例中,图8为本申请另一个实施例中堆芯中子学模型的修正方法的流程示意图,在上述实施例的基础上,本实施例堆芯中子学模型的修正方法还可以包括以下步骤。
S801、根据修正后的堆芯中子学模型确定目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率。
本申请实施例中,服务器可以根据修正后的堆芯中子学模型,重新计算得到目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率。
这样,通过上述实施方式,服务器基于修正后的堆芯中子学模型,可以更加准确得计算得到各节块的理论快群中子通量、理论热群中子通量和理论功率,有利于更加准确地监测堆芯内的核反应情况。
应该理解的是,虽然如上的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,如上的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段,这些步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
基于同样的发明构思,本申请实施例还提供了一种用于实现上述所涉及的堆芯中子学模型的修正方法的堆芯中子学模型的修正装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似,故下面所提供的一个或多个堆芯中子学模型的修正装置实施例中的具体限定可以参见上文中对于堆芯中子学模型的修正方法的限定,在此不再赘述。
在一个实施例中,图9为本申请一个实施例中堆芯中子学模型的修正装置的结构示意图,本申请实施例提供的堆芯中子学模型的修正装置可以应用于计算机设备中。如图9所示,本申请实施例的堆芯中子学模型的修正装置,包括:第一确定模块10、第二确定模块11和第三确定模块12,其中:
第一确定模块10,用于根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定目标堆芯的各节块的测量快群中子通量和测量热群中子通量。
第二确定模块11,用于根据各节块的状态参数确定各节块的等效截面参数。
第三确定模块12,用于根据各节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定修正后的堆芯中子学模型。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第三确定模块12,包括:第一确定单元和第二确定单元,其中:
第一确定单元,对于各节块,用于根据节块的测量快群中子通量、测量热群中子通量和等效截面参数,确定第一泄漏曲率参数。
第二确定单元,用于根据第一泄漏曲率参数,确定修正后的堆芯中子学模型。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第一确定单元具体用于:
根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和快群泄漏修正量,确定快群中子平衡方程,其中,快群泄漏修正量中包括第一泄漏曲率参数;
根据节块的测量快群中子通量、测量热群中子通量、等效截面参数和热群泄漏修正量,确定热群中子平衡方程,其中,热群泄漏修正量中包括第一泄漏曲率参数;
根据快群中子平衡方程和热群中子平衡方程,确定第一泄漏曲率参数。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第二确定单元具体用于:
根据第一泄漏曲率参数分别对快群中子平衡方程和热群中子平衡方程进行修正,得到修正后的堆芯中子学模型。
在一个实施例中,上述第二确定模块11具体用于:
对于各节块,根据预设状态参数集合确定节块对应的多个状态参数;根据多个状态参数按照预设插值算法,确定节块的等效截面参数。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第一确定模块10具体用于:
对于各节块,根据修正前的堆芯中子学模型确定节块的理论快群中子通量、理论热群中子通量、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面;
根据节块的理论快群中子通量、节块的理论热群中子通量、节块的测量功率、节块的快群宏观能量产生截面和节块的热群宏观能量产生截面,确定节块的测量快群中子通量和测量热群中子通量。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,堆芯中子学模型的修正装置还包括:第四确定模块、第五确定模块、第六确定模块和第七确定模块。
其中,第四确定模块,用于根据扰动堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率;
第五确定模块,用于根据各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率,确定第二泄漏曲率参数;
第六确定模块,用于根据第二泄漏曲率参数,确定修正后的扰动堆芯中子学模型;
第七确定模块,用于基于各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定修正后的扰动堆芯中子学模型是否正确。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第七确定模块包括:第三确定单元和第四确定单元。
其中,第三确定单元,用于根据修正后的扰动堆芯中子学模型,确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率;
第四确定单元,用于根据各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及各扰动节块的测量功率与对应的修正功率之间的差值,确定修正后的扰动堆芯中子学模型是否正确。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述第四确定单元具体用于:
若各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值小于第一预设阈值、各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值小于第二预设阈值,且各扰动节块的测量功率与对应的修正功率之间的差值小于第三预设阈值,则确定修正后的扰动堆芯中子模型正确。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
在一个实施例中,上述堆芯中子学模型的修正装置还包括:
第八确定模块,用于根据修正后的堆芯中子学模型确定目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率。
本实施例提供的堆芯中子学模型的修正装置,可以执行上述方法实施例,其实现原理和技术效果类似,在此不再赘述。
上述堆芯中子学模型的修正装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个示例性的实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图10所示。该计算机设备可以包括处理器、存储器、输入/输出接口(Input/Output,简称I/O)和通信接口。其中,处理器、存储器和输入/输出接口通过系统总线连接,通信接口通过输入/输出接口连接到系统总线。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质和内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的输入/输出接口用于处理器与外部设备之间交换信息。该计算机设备的通信接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种堆芯中子学模型的修正方法。
本领域技术人员可以理解,图10中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现本申请上述堆芯中子学模型的修正方法实施例中关于服务器的技术方案,其实现原理和技术效果类似,此处不再赘述。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现本申请上述堆芯中子学模型的修正方法实施例中关于服务器的技术方案,其实现原理和技术效果类似,此处不再赘述。
在一个实施例中,提供了一种计算机程序产品,包括计算机程序,该计算机程序被处理器执行时实现本申请上述堆芯中子学模型的修正方法实施例中关于服务器的技术方案,其实现原理和技术效果类似,此处不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、磁变存储器(Magnetoresistive Random Access Memory,MRAM)、铁电存储器(FerroelectricRandom Access Memory,FRAM)、相变存储器(Phase Change Memory,PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器等。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(StaticRandom Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory,DRAM)等。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等,不限于此。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利要求为准。
Claims (14)
1.一种堆芯中子学模型的修正方法,其特征在于,所述方法包括:
根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定所述目标堆芯的各节块的测量快群中子通量和测量热群中子通量;
根据所述各节块的状态参数确定所述各节块的等效截面参数;
根据所述各节块的测量快群中子通量、所述测量热群中子通量和所述等效截面参数,确定修正后的堆芯中子学模型。
2.根据权利要求1所述的方法,其特征在于,所述根据所述各节块的测量快群中子通量、所述测量热群中子通量和所述等效截面参数,确定修正后的堆芯中子学模型,包括:
对于所述各节块,根据所述节块的测量快群中子通量、所述测量热群中子通量和所述等效截面参数,确定第一泄漏曲率参数;
根据所述第一泄漏曲率参数,确定所述修正后的堆芯中子学模型。
3.根据权利要求2所述的方法,其特征在于,所述根据所述节块的测量快群中子通量、所述测量热群中子通量和所述等效截面参数,确定第一泄漏曲率参数,包括:
根据所述节块的测量快群中子通量、所述测量热群中子通量、所述等效截面参数和快群泄漏修正量,确定快群中子平衡方程,其中,所述快群泄漏修正量中包括所述第一泄漏曲率参数;
根据所述节块的测量快群中子通量、所述测量热群中子通量、所述等效截面参数和热群泄漏修正量,确定热群中子平衡方程,其中,所述热群泄漏修正量中包括所述第一泄漏曲率参数;
根据所述快群中子平衡方程和所述热群中子平衡方程,确定所述第一泄漏曲率参数。
4.根据权利要求2所述的方法,其特征在于,所述根据所述第一泄漏曲率参数,确定所述修正后的堆芯中子学模型,包括:
根据所述第一泄漏曲率参数分别对所述快群中子平衡方程和所述热群中子平衡方程进行修正,得到所述修正后的堆芯中子学模型。
5.根据权利要求1-4中任一项所述的方法,其特征在于,所述根据所述各节块的状态参数确定所述各节块的等效截面参数,包括:
对于所述各节块,根据预设状态参数集合确定所述节块对应的多个状态参数;
根据所述多个状态参数按照预设插值算法,确定所述节块的等效截面参数。
6.根据权利要求1-4中任一项所述的方法,其特征在于,所述目标堆芯的测量功率包括所述目标堆芯的各节块的测量功率,所述根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定所述目标堆芯的各节块的测量快群中子通量和测量热群中子通量,包括:
对于所述各节块,根据所述修正前的堆芯中子学模型确定所述节块的理论快群中子通量、理论热群中子通量、所述节块的快群宏观能量产生截面和所述节块的热群宏观能量产生截面;
根据所述节块的理论快群中子通量、所述节块的理论热群中子通量、所述节块的测量功率、所述节块的快群宏观能量产生截面和所述节块的热群宏观能量产生截面,确定所述节块的测量快群中子通量和测量热群中子通量。
7.根据权利要求1-4中任一项所述的方法,其特征在于,所述方法还包括:
根据扰动堆芯中子学模型确定各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率;
根据所述各扰动节块的等效截面参数、测量快群中子通量、测量热群中子通量和测量功率,确定第二泄漏曲率参数;
根据所述第二泄漏曲率参数,确定修正后的扰动堆芯中子学模型;
基于所述各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定所述修正后的扰动堆芯中子学模型是否正确。
8.根据权利要求7所述的方法,其特征在于,所述基于所述各扰动节块的测量快群中子通量、测量热群中子通量和测量功率确定所述修正后的扰动堆芯中子学模型是否正确,包括:
根据所述修正后的扰动堆芯中子学模型,确定各扰动节块的修正快群中子通量、修正热群中子通量和修正功率;
根据所述各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、所述各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及所述各扰动节块的测量功率与对应的修正功率之间的差值,确定所述修正后的扰动堆芯中子学模型是否正确。
9.根据权利要求8所述的方法,其特征在于,所述根据所述各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值、所述各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值,以及所述各扰动节块的测量功率与对应的修正功率之间的差值,确定所述修正后的扰动堆芯中子学模型是否正确,包括:
若所述各扰动节块的测量快群中子通量与对应的修正快群中子通量之间的差值小于第一预设阈值、所述各扰动节块的测量热群中子通量与对应的修正热群中子通量之间的差值小于第二预设阈值,且所述各扰动节块的测量功率与对应的修正功率之间的差值小于第三预设阈值,则确定所述修正后的扰动堆芯中子模型正确。
10.根据权利要求1-4中任一项所述的方法,其特征在于,所述方法还包括:
根据所述修正后的堆芯中子学模型确定所述目标堆芯的各节块的理论快群中子通量、理论热群中子通量和理论功率。
11.一种堆芯中子学模型的修正装置,其特征在于,所述装置包括:
第一确定模块,用于根据修正前的堆芯中子学模型和目标堆芯的测量功率,确定所述目标堆芯的各节块的测量快群中子通量和测量热群中子通量;
第二确定模块,用于根据所述各节块的状态参数确定所述各节块的等效截面参数;
第三确定模块,用于根据所述各节块的测量快群中子通量、所述测量热群中子通量和所述等效截面参数,确定修正后的堆芯中子学模型。
12.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至10任意一项所述方法的步骤。
13.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至10任意一项所述方法的步骤。
14.一种计算机程序产品,包括计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至10任意一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311458625.8A CN117594267A (zh) | 2023-11-03 | 2023-11-03 | 堆芯中子学模型的修正方法、装置、设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311458625.8A CN117594267A (zh) | 2023-11-03 | 2023-11-03 | 堆芯中子学模型的修正方法、装置、设备和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117594267A true CN117594267A (zh) | 2024-02-23 |
Family
ID=89910653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311458625.8A Pending CN117594267A (zh) | 2023-11-03 | 2023-11-03 | 堆芯中子学模型的修正方法、装置、设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117594267A (zh) |
-
2023
- 2023-11-03 CN CN202311458625.8A patent/CN117594267A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ryu et al. | Solution of the BEAVRS benchmark using the nTRACER direct whole core calculation code | |
Daeubler et al. | High-fidelity coupled Monte Carlo neutron transport and thermal-hydraulic simulations using Serpent 2/SUBCHANFLOW | |
Ivanov et al. | Challenges in coupled thermal–hydraulics and neutronics simulations for LWR safety analysis | |
Shimjith et al. | Space–time kinetics modeling of advanced heavy water reactor for control studies | |
Kerlin et al. | Dynamics and control of nuclear reactors | |
Gill et al. | Numerical methods in coupled Monte Carlo and thermal-hydraulic calculations | |
Zhang et al. | An assessment of coupling algorithms in HTR simulator TINTE | |
RU2550689C2 (ru) | Способ содействия работе ядерного реактора | |
Choi et al. | Development of high-fidelity neutron transport code STREAM | |
Ikonen et al. | Module for thermomechanical modeling of LWR fuel in multiphysics simulations | |
Magedanz et al. | High-fidelity multi-physics system TORT-TD/CTF/FRAPTRAN for light water reactor analysis | |
Mirvakili et al. | Developing a computational tool for predicting physical parameters of a typical VVER-1000 core based on artificial neural network | |
Valtavirta et al. | Coupled neutronics–fuel behavior calculations in steady state using the Serpent 2 Monte Carlo code | |
Yu et al. | Fuel performance analysis of BEAVRS benchmark Cycle 1 depletion with MCS/FRAPCON coupled system | |
Lei et al. | Development and validation of a deep learning‐based model for predicting burnup nuclide density | |
Valtavirta et al. | Validating the Serpent-Ants calculation chain using BEAVRS fresh core hot zero power data | |
Tuominen et al. | BEAVRS pin-by-pin calculations with Ants-SUBCHANFLOW-SuperFINIX code system | |
WO2024093077A1 (zh) | 核电厂导致芯块与包壳相互作用的事故瞬态分析方法 | |
Xia et al. | Performance evaluation of a 3-D kinetic model for CANDU reactors in a closed-loop environment | |
CN117594267A (zh) | 堆芯中子学模型的修正方法、装置、设备和存储介质 | |
Blakely et al. | Demonstration of LOTUS multiphysics BEPU analysis framework for LB-LOCA simulations | |
CN114582532A (zh) | 燃料组件的截面确定方法、装置、设备、介质和程序产品 | |
Duchnowski et al. | Neutronic and thermal-hydraulic calculations under steady state and transient conditions for a generic pebble bed fluoride salt cooled high-temperature reactor | |
Wysocki et al. | Coupling of CTF and RELAP5-3D Within an Enhanced-Fidelity Nuclear Power Plant Simulator | |
Jiménez et al. | Comparative analysis of neutronics/thermal-hydraulics multi-scale coupling for LWR analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |