Quadrant power tilt control method, device and computer equipment for reactor core

By determining the quadrant power inclination change rate and target parameters in the pressurized water reactor nuclear power plant, the use of combustible poisons is solved, and the operational safety of the reactor is improved.

CN117747163BActive Publication Date: 2025-06-27CHINA NUCLEAR POWER TECH RES INST CO LTD
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202311517797.8
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2023-11-15
Publication Date
2025-06-27
Estimated Expiration
2043-11-15

AI Technical Summary

Technical Problem

How to control the increase in the quadrant power in the core in the pressurized water reactor nuclear power plant, affecting the operation and safety margin of the reactor.

Method used

The target parameters of the combustible poison are determined by determining the rate of change of the quadrant power inclination caused by fuel combustion in the fuel cycle containing the combustible poison, and the quadrant power inclination of the core is controlled according to these parameters.

Benefits of technology

Effectively control the quadrant power tilt of the core, reducing the negative impact on reactor operation and safety margin.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN117747163B_ABST
    Figure CN117747163B_ABST
Patent Text Reader

Abstract

The present application relates to a method, device and computer equipment for controlling the quadrant power tilt of a reactor core. The method includes: determining a first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poison according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core in the fuel cycle with burnable poison; then, according to the first expression of the first change rate of the quadrant power tilt, determining the target parameters of the burnable poison that affect the quadrant power tilt of the reactor core; and controlling the quadrant power tilt of the reactor core according to the target parameters of the burnable poison, so as to control the increase of the quadrant power tilt.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present application relates to the field of nuclear power plant control technology, and in particular to a method, device and computer equipment for controlling quadrant power tilt of a core. Background Art

[0002] Quadrant Power Tilt Ratio (QPTR) refers to the ratio between the power of the four quadrants of the reactor and their average value. It is an indicator to measure the symmetry of the radial power distribution of the core. Therefore, quadrant power tilt is an important safety factor for daily monitoring of pressurized water reactor nuclear power plants. Although the increase in quadrant power tilt will not directly affect the safe operation of the reactor, it will indirectly affect the operation and safety margin of the reactor.

[0003] Therefore, how to control the increase of quadrant power tilt becomes a technical problem to be solved urgently in this field. Summary of the invention

[0004] Based on this, it is necessary to provide a quadrant power tilt control method, device and computer equipment for a core that can control the increase of quadrant power tilt in order to solve the above technical problems.

[0005] In a first aspect, the present application provides a method for controlling quadrant power tilt of a core. The method comprises:

[0006] Determine a first expression for a first change rate of quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison according to a first change amount of quadrant power tilt caused by fuel combustion in each quadrant of the core in the fuel cycle containing burnable poison, a second change amount of quadrant power tilt caused by burning of burnable poison, and a burnup increment of the core;

[0007] Determining a target parameter of the burnable poison affecting the quadrant power tilt of the core according to a first expression for a first rate of change of the quadrant power tilt;

[0008] The quadrant power tilt of the core is controlled according to the target parameters of the burnable poison.

[0009] In one embodiment, the first expression for determining the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison according to the first change amount of the quadrant power tilt caused by the combustion of fuel in each quadrant of the core in the fuel cycle containing burnable poison, the second change amount of the quadrant power tilt caused by the combustion of burnable poison, and the burnup increment of the core includes:

[0010] Dividing the first change amount of the quadrant power tilt by the fuel consumption increment, a second expression for a second change rate of the quadrant power tilt is obtained;

[0011] Dividing a second change amount of the quadrant power tilt by the burnup increment to obtain a third expression of a third change rate of the quadrant power tilt;

[0012] Determining the first expression according to the second expression and the third expression.

[0013] In one embodiment, the determining the first expression according to the second expression and the third expression includes:

[0014] Adding the second expression and the third expression to obtain the first expression.

[0015] In one embodiment, the method further includes:

[0016] Determining a fourth expression of a fourth change rate of the quadrant power tilt corresponding to each quadrant in a fuel cycle without burnable poisons according to the relationship between the quadrant power of each quadrant and the burnup increment;

[0017] The determining a target parameter affecting the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt includes:

[0018] Determining the target parameter of the burnable poison affecting the quadrant power tilt of the core according to the fourth expression and the first expression.

[0019] In one embodiment, the determining the target parameter of the burnable poison affecting the quadrant power tilt of the core according to the fourth expression and the first expression includes:

[0020] Determining the similarity between the fourth expression and the second expression;

[0021] When the similarity is greater than a preset similarity threshold, determining the target parameter according to the third expression of the third change rate of the quadrant power tilt.

[0022] In one embodiment, the method further includes:

[0023] For each quadrant, determining a product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power of the burnable poison assembly with burnable poison in the quadrant, the burnup increment, and the slope corresponding to the burnable poison assembly with y roots; the slope is the slope of the curve segment where the reactivity increases with the increase of burnup in the K-infinity vs. burnup change curve, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant;

[0024] Determining a summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant;

[0025] Determine a second change amount of the quadrant power tilt based on the summation result and the quadrant power of this quadrant;

[0026] Determine the target parameter according to a third expression of a third change rate of the quadrant power tilt, including:

[0027] Determine the target parameter according to the slope in the third expression.

[0028] In one embodiment, controlling the quadrant power tilt of the core according to the target parameter of the burnable poison includes:

[0029] If the target parameter is greater than a first preset parameter, prohibit the use of the burnable poison to control the quadrant power tilt of the core;

[0030] If the target parameter is less than a second preset parameter, control the quantity of the burnable poison used to control the quadrant power tilt of the core.

[0031] In a second aspect, the present application also provides a device for controlling the quadrant power tilt of a core. The device includes:

[0032] A first determination module, configured to determine a first expression of a first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison based on a first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the core, a second change amount of the quadrant power tilt caused by burnable poison combustion, and a burnup increment of the core;

[0033] A second determination module, configured to determine a target parameter of the burnable poison that affects the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt;

[0034] A control module, configured to control the quadrant power tilt of the core according to the target parameter of the burnable poison.

[0035] In a third aspect, the present application also provides a computer device, including a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the steps of any of the above methods are implemented.

[0036] In a fourth aspect, the present application also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of any of the above methods are implemented.

[0037] The above-mentioned method, device, computer equipment, storage medium and program product for controlling the quadrant power tilt of the core determine the first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poison according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by the combustion of burnable poison, and the burnup increment of the core in the fuel cycle with burnable poison. Then, according to the first expression of the first change rate of the quadrant power tilt, the target parameter of the burnable poison affecting the quadrant power tilt of the core is determined, and the quadrant power tilt of the core is controlled according to the target parameter of the burnable poison, so as to control the increase of the quadrant power tilt. Description of the Drawings

[0038] Figure 1 is the internal structure diagram of a computer device provided by an embodiment of the present application;

[0039] Figure 2 is the schematic flowchart of a method for controlling the quadrant power tilt of a core provided by an embodiment of the present application;

[0040] Figure 3 are two possible quadrant division methods provided by an embodiment of the present application;

[0041] Figure 4 is the schematic diagram showing the change of the quadrant power tilt TILT with burnup in a fuel cycle provided by an embodiment of the present application;

[0042] Figure 5 is the schematic flowchart of a method for obtaining the first expression provided by an embodiment of the present application;

[0043] Figure 6 is the schematic flowchart of a method for obtaining the target parameter provided by an embodiment of the present application;

[0044] Figure 7 is the schematic diagram showing the change of the quadrant power tilt TILT with burnup in another fuel cycle provided by an embodiment of the present application;

[0045] Figure 8 is the schematic flowchart of another method for obtaining the target parameter provided by an embodiment of the present application;

[0046] Figure 9 is the schematic flowchart of yet another method for obtaining the target parameter provided by an embodiment of the present application;

[0047] Figure 10 is the schematic diagram of the curve of K-infinity changing with burnup provided by an embodiment of the present application;

[0048] Figure 11It is a flow chart of a core design method for controlling quadrant power tilt increase during operation of a pressurized water reactor nuclear power plant provided by an embodiment of the present application;

[0049] Figure 12 It is a structural block diagram of quadrant power tilt control of a core provided in an embodiment of the present application. DETAILED DESCRIPTION

[0050] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0051] Quadrant Power Tilt Ratio (QPTR) refers to the ratio between the power of the four quadrants of the reactor and their average value. It is an indicator to measure the symmetry of the radial power distribution of the core. Therefore, quadrant power tilt is an important safety factor for daily monitoring of pressurized water reactor nuclear power plants. Although the increase in quadrant power tilt will not directly affect the safe operation of the reactor, it will indirectly affect the operation and safety margin of the reactor.

[0052] Therefore, how to control the increase of quadrant power tilt becomes a technical problem to be solved urgently in this field.

[0053] The quadrant power tilt control method of the core provided in the embodiment of the present application can be applied to Figure 1 In the application environment shown. Figure 1 is an internal structure diagram of a computer device provided in an embodiment of the present application. The computer device may be a server, and its internal structure diagram may be as follows Figure 1 As shown. The computer device includes a processor, a memory and a network interface connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a quadrant power tilt control method of a core is implemented.

[0054] Those skilled in the art will understand that Figure 1 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0055] In one embodiment, as Figure 2 shown, Figure 2 FIG. is a schematic flowchart of a method for controlling the quadrant power tilt of a reactor core provided by an embodiment of the present application. This method can be applied to Figure 1 a computer device, and the method includes the following steps:

[0056] S201. Determine a first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poisons according to the first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the reactor core in the fuel cycle with burnable poisons, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core.

[0057] In the embodiment of the present application, the reactor core can be evenly divided into four sub-core regions, and the regions where the four sub-core regions are located are used as four quadrants. Exemplarily, as Figure 3 shown, Figure 3 FIG. are two possible quadrant division methods provided by the embodiment of the present application. The reactor core can be evenly divided into four sub-core regions according to the coordinate axis or diagonal, and the regions where the four sub-core regions are located are used as four quadrants.

[0058] In the embodiment of the present application, the fuel cycle with burnable poisons can be a fuel cycle in which burnable poisons are contained in the reactor core. Among them, the burnable poisons can be boron poison, gadolinium poison, erbium poison, etc. The burnup increment of the reactor core can be the burnup increase amount caused by fuel combustion or burnable poison combustion in the fuel cycle, and the quadrant power tilt change amount can be the change amount of the quadrant power tilt when the reactor core burnup increases in the fuel cycle.

[0059] In the embodiment of the present application, the quadrant power tilt can be represented by TILT, and the TILT of each quadrant of the reactor core can be defined as where Pi is the quadrant power of the i-th quadrant, which can be defined as TILTi is the quadrant power tilt of the i-th quadrant, Pj is the component power of the j-th component in the i-th quadrant, and the total number of components in the i-th quadrant is N.

[0060] Exemplarily, as Figure 4 shown, Figure 4 FIG. is a schematic diagram showing the change of the quadrant power tilt TILT with burnup in a fuel cycle provided by the embodiment of the present application. At the beginning of the cycle, in the fourth quadrant at the lower part of the reactor core, the TILT value is the largest, about 0.4%; in the second quadrant at the upper part of the reactor core, which is the opposite quadrant, the TILT value is the smallest, about -0.7%. As the reactor core burnup increases, the TILT value of the fourth quadrant becomes larger and starts to decrease rapidly at a burnup of 10000 MWd / tU; on the contrary, the TILT value of the second quadrant, which is the opposite quadrant, becomes smaller and starts to increase rapidly at a burnup of 10000 MWd / tU.

[0061] S202. Determine the target parameters of the burnable poison that affects the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt.

[0062] In the embodiment of the present application, the first change rate of the quadrant power tilt may be the ratio of the change amount of the quadrant power tilt to the burnup increment, that is, the first expression may be ΔTILT / ΔBu, where ΔTILT is the change amount of the quadrant power tilt when the core burnup increases from Bu1 to Bu2, and ΔBu is the burnup increment.

[0063] In the embodiment of the present application, a target parameter interval may be preset according to the target parameters of the burnable poison, and the type and dosage of the burnable poison are determined according to the target parameter interval, so as to control the quadrant power tilt of the core.

[0064] S203. Control the quadrant power tilt of the core according to the target parameters of the burnable poison.

[0065] When designing the core loading scheme, the quadrant power tilt of the core can be controlled according to the target parameters of the burnable poison. Exemplarily, if the target parameter is greater than the first preset parameter, the burnable poison corresponding to the target parameter is prohibited from being used to control the quadrant power tilt of the core. Or, if the target parameter is less than the second preset parameter, the quantity of the burnable poison corresponding to the target parameter is controlled to control the quadrant power tilt of the core.

[0066] In the embodiment of the present application, according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the core in the fuel cycle containing burnable poison in each quadrant of the core, the first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison is determined, and then according to the first expression of the first change rate of the quadrant power tilt, the target parameters that affect the quadrant power tilt of the core are determined, where the target parameters are used to control the quadrant power tilt of the core, so as to control the increase of the quadrant power tilt.

[0067] Refer to Figure 5 , Figure 5 is a schematic flowchart of a method for obtaining a first expression provided by an embodiment of the present application. This embodiment relates to a possible implementation manner of how to determine the first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the core in each quadrant of the core. On the basis of the above embodiment, the above S201 includes the following steps:

[0068] S501. Divide the first change amount of the quadrant power tilt by the burnup increment to obtain a second expression of the second change rate of the quadrant power tilt.

[0069] In an embodiment of the present application, the second expression may be ΔTILTBu / ΔBu, where ΔTILTBu is the first change amount of the quadrant power tilt caused by fuel combustion in a fuel cycle with burnable poisons when the core burnup increases from Bu1 to Bu2, and ΔBu is the burnup increment. Therefore, ΔBu = Bu2 - Bu1.

[0070] S502. Divide the second change amount of the quadrant power tilt by the burnup increment to obtain a third expression of the third change rate of the quadrant power tilt.

[0071] In an embodiment of the present application, the third expression may be ΔTILTGd / ΔBu, where ΔTILTGd is the first change amount of the quadrant power tilt caused by burnable poison combustion in a fuel cycle with burnable poisons when the core burnup increases from Bu1 to Bu2, and ΔBu is the burnup increment.

[0072] S503. Determine the first expression according to the second expression and the third expression.

[0073] In an embodiment of the present application, the first change rate of the quadrant power tilt corresponding to the first expression may be the ratio of the change amount of the quadrant power tilt to the burnup increment. Among them, the change amount of the quadrant power tilt and the burnup increment are jointly caused by fuel combustion and burnable poison combustion in a fuel cycle with burnable poisons. Therefore, the change amount of the quadrant power tilt and the burnup increment can be determined according to the first change amount of the quadrant power tilt and the burnup increment in the second expression, and the second change amount of the quadrant power tilt and the burnup increment in the third expression. Then, the first expression can be determined according to the change amount of the quadrant power tilt and the burnup increment.

[0074] In an embodiment of the present application, the first expression can be obtained by adding the second expression and the third expression. Or, the first expression can be obtained by adding the second expression and the third expression after multiplying them by a preset coefficient respectively.

[0075] In the embodiment of the present application, the first change amount of the quadrant power tilt is divided by the burnup increment to obtain a second expression of the second change rate of the quadrant power tilt. Then, the second change amount of the quadrant power tilt is divided by the burnup increment to obtain a third expression of the third change rate of the quadrant power tilt. Finally, the first expression is determined according to the second expression and the third expression. Since the change amount of the quadrant power tilt and the burnup increment are jointly caused by fuel combustion and burnable poison combustion in the fuel cycle with burnable poison, the first expression can be determined according to the second expression and the third expression. Then, the target parameter affecting the quadrant power tilt of the core is determined according to the first expression, where the target parameter is used to control the quadrant power tilt of the core, so as to control the increase of the quadrant power tilt.

[0076] Based on the above embodiment, for the above S503, determining the first expression according to the second expression and the third expression can be achieved in the following manner:

[0077] Add the second expression and the third expression to obtain the first expression.

[0078] In the embodiment of the present application, the first expression can be obtained by adding the second expression and the third expression. Alternatively, the first expression can be obtained by multiplying the second expression and the third expression by preset coefficients respectively and then adding them.

[0079] Exemplarily, the first expression can be ΔTILT / ΔBu = ΔTILTBu / ΔBu + ΔTILTGd / ΔBu.

[0080] In the embodiment of the present application, the first expression is obtained by adding the second expression and the third expression. Since the change amount of the quadrant power tilt and the burnup increment are jointly caused by fuel combustion and burnable poison combustion in the fuel cycle with burnable poison, the first expression is obtained by adding the second expression and the third expression. Then, the target parameter affecting the quadrant power tilt of the core is determined according to the first expression, where the target parameter is used to control the quadrant power tilt of the core, so as to control the increase of the quadrant power tilt.

[0081] Refer to Figure 6 , Figure 6 is a schematic flowchart of a method for obtaining a target parameter provided by an embodiment of the present application. Based on the above embodiment, the method includes the following steps:

[0082] S601, according to the relationship between the quadrant power of each quadrant and the burnup increment, determine a fourth expression of the fourth change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poison.

[0083] In the embodiments of the present application, the relationship between the quadrant power and the burnup increment in each quadrant can be ΔBu = P×Δt. Wherein, P is the quadrant power, and ΔBu is the burnup increment within the time increment Δt. The fourth expression for each quadrant can be the ratio of the quadrant power tilt change amount ΔTILT' to the burnup increment ΔBu in each quadrant, that is, ΔTILT' / ΔBu.

[0084] According to the relationship between the quadrant power and the burnup increment, it can be concluded that the burnup has a "self - inhibition" function on the power distribution. That is, when the power at a certain place is high, the amplitude of burnup deepening is also large, and the relatively large burnup will make the relative power at that place tend to decrease; on the contrary, when the power at a certain place is low, the amplitude of burnup deepening is also small, and the relatively small burnup will make the relative power at that place tend to increase. That is, the fourth expression can reflect the variation law of the "self - inhibition" function of burnup on the power distribution.

[0085] Exemplarily, as Figure 7 shown, Figure 7 is a schematic diagram showing the change of the quadrant power tilt TILT with burnup in another fuel cycle provided by the embodiments of the present application. At the beginning of the cycle, in the 7th quadrant at the lower part of the core, the TILT value is the largest, about 2%; in its reverse quadrant, the 5th quadrant at the upper part of the core, the TILT value is the smallest, about - 2%. As the core burnup increases, in the 7th quadrant, due to the relatively high power in this quadrant, the burnup increment is large, and the large burnup causes its TILT value to continuously decrease and finally tend to a stable value; on the contrary, in the reverse quadrant, the 5th quadrant, due to the relatively small power with the smallest TILT value in this quadrant, the burnup increment is small, and the small burnup causes its TILT value to continuously increase and finally tend to a stable value. That is, the above - mentioned burnup has a "self - inhibition" function on the power distribution.

[0086] S602. Determine the target parameters affecting the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt, including:

[0087] Determine the target parameters of the burnable poison affecting the quadrant power tilt of the core according to the fourth expression and the first expression.

[0088] In the embodiments of the present application, determine the target parameters of the burnable poison affecting the quadrant power tilt of the core according to the relationship between the quadrant power and the burnup increment corresponding to the fourth expression and the first change rate of the quadrant power tilt corresponding to the first expression.

[0089] In the embodiment of the present application, according to the relationship between the quadrant power and the burnup increment of each quadrant, a fourth expression of the fourth change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poisons is determined. According to the fourth expression and the first expression, the target parameter of the burnable poison affecting the quadrant power tilt of the reactor core is determined. Since the target parameter is used to control the quadrant power tilt of the reactor core, it is possible to control the increase in the quadrant power tilt.

[0090] Referring to Figure 8 , Figure 8 FIG. is a schematic flowchart of another method for obtaining a target parameter provided by an embodiment of the present application. This embodiment relates to a possible implementation manner of determining the target parameter affecting the quadrant power tilt of the reactor core according to the fourth expression and the first expression. On the basis of the above embodiment, the above S602 includes the following steps:

[0091] S801, determine the similarity between the fourth expression and the second expression.

[0092] S802, when the similarity is greater than a preset similarity threshold, determine the target parameter according to the third expression of the third change rate of the quadrant power tilt.

[0093] In the embodiment of the present application, the similarity between the fourth expression and the second expression is determined. According to the similarity and the preset similarity threshold, if the similarity is greater than the preset similarity threshold, it is determined that the fourth expression is similar to the second expression. Since the burnup has a "self-inhibition" function on the power distribution in the fourth expression, the second expression also has a "self-inhibition" function and will not increase the quadrant power tilt. Therefore, it is determined that the third expression is the main factor for increasing the quadrant power tilt, and the target parameter is determined according to the third expression of the third change rate of the quadrant power tilt.

[0094] In the embodiment of the present application, the similarity between the fourth expression and the second expression is determined. When the similarity is greater than the preset similarity threshold, the target parameter is determined according to the third expression of the third change rate of the quadrant power tilt. Since the target parameter is used to control the quadrant power tilt of the reactor core, it is possible to control the increase in the quadrant power tilt.

[0095] Referring to Figure 9 , Figure 9 FIG. is a schematic flowchart of yet another method for obtaining a target parameter provided by an embodiment of the present application. On the basis of the above embodiment, the method includes the following steps:

[0096] S901. For each quadrant, determine the product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power, burnup increment, and slope corresponding to the burnable poison assembly containing burnable poison in the quadrant; the slope is the slope of the curve segment where the reactivity increases with the increase of burnup in the K-infinity vs. burnup change curve graph, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant.

[0097] Exemplarily, the graph of the relationship between K-infinity (Kinf) and burnup can be as Figure 10 shown. Figure 10 It is a schematic diagram of the K-infinity vs. burnup change curve provided by an embodiment of the present application. Specifically, Figure 10 shown is the graph of the change of K-infinity with burnup of the assembly with a U-235 enrichment of 4.45% and containing 0, 8, 12, 16, and 24 gadolinium rod assemblies in each component respectively.

[0098] In the embodiment of the present application, exemplarily, the power of the m-th burnable poison assembly containing burnable poison in the quadrant can be P m , the burnup increment can be ΔBu, and the slope of the burnable poison assembly containing y burnable poisons can be X y . Therefore, the product result can be X y *ΔBu*Pm m .

[0099] S902. Determine the summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant.

[0100] In the embodiment of the present application, the product results corresponding to each burnable poison assembly in the quadrant can be added in sequence to obtain the summation result. Exemplarily, the summation result can be ∑ m X y *ΔBu*P m , where m = 1,..., M, y = 1,..., T. Here, M is the total number of burnable poison assemblies in the quadrant, and T is the total number of types of burnable poison assemblies in the core.

[0101] S903. Determine the second change amount of the quadrant power tilt according to the summation result and the quadrant power of the quadrant.

[0102] The quotient obtained by dividing the summation result by the quadrant power of the quadrant can be used as the second change amount of the quadrant power tilt. Alternatively, the result obtained by multiplying the quotient obtained by dividing the summation result by the quadrant power of the quadrant by a preset coefficient can be used as the second change amount of the quadrant power tilt.

[0103] Exemplarily, the second change amount of the quadrant power tilt can be Wherein, M is the total number of combustible poison assemblies in the quadrant, T is the total number of types of combustible poison assemblies in the core, and P i is the quadrant power of the i-th quadrant.

[0104] S904. Determine the target parameter according to the third expression of the third change rate of the quadrant power tilt, including:

[0105] Determine the target parameter according to the slope in the third expression.

[0106] In the embodiment of the present application, according to the first expression, the third expression, and the K-infinity vs. burnup change curve, it is determined that the slope in the third expression can be used to control the quadrant power tilt of the core. Therefore, the slope in the third expression is used as the target parameter.

[0107] In the embodiment of the present application, for each quadrant, determine the product result corresponding to the type of each combustible poison assembly. Wherein, the product result is the product of the power, burnup increment, and slope of the combustible poison assembly containing combustible poison in the quadrant. The slope is the curve slope of the section where the reactivity increases with the increase of burnup in the K-infinity vs. burnup change curve. Then, according to the product results corresponding to the quantities of combustible poison assemblies, determine the summation result. According to the summation result and the quadrant power of the quadrant, determine the third change rate of the quadrant power tilt. According to the slope in the third expression, determine the target parameter. Since the target parameter is used to control the quadrant power tilt of the core, it is possible to control the increase of the quadrant power tilt.

[0108] Based on the above embodiment, in S203, controlling the quadrant power tilt of the core according to the target parameter of the combustible poison can be achieved in the following manner:

[0109] If the target parameter is greater than the first preset parameter, prohibit the use of combustible poison to control the quadrant power tilt of the core.

[0110] If the target parameter is less than the second preset parameter, control the quantity of combustible poison used to control the quadrant power tilt of the core.

[0111] When designing the core loading scheme, if the target parameter is greater than the first preset parameter, prohibit the use of combustible poison to control the quadrant power tilt of the core. That is, when designing the core loading scheme, do not use the type of combustible poison with a large rising slope in the curve section where the reactivity increases with the increase of burnup. Taking the assembly with uranium-235 enrichment of 4.45% and gadolinium enrichment of 8% using uranium dioxide fuel as an example, do not use the type of assembly containing 24 gadolinium rods.

[0112] When designing the core loading scheme, if the target parameter is less than the second preset parameter, control the quantity of burnable poison used to control the quadrant power tilt of the core. That is, when designing the core loading scheme, try not to use burnable poison of the type with a smaller rising slope in the curve where the reactivity increases with burnup. If it is necessary to use the above-mentioned burnable poison of the type with a smaller rising slope in the curve where the reactivity increases with burnup, its total quantity should be controlled. Taking the assembly with uranium-235 enrichment of 4.45% and gadolinium enrichment of 8% using uranium dioxide fuel as an example, try not to use the type of assembly with 20 gadolinium rods. If it is necessary to use the above-mentioned type of assembly with 20 gadolinium rods, the number of assemblies of the type with 20 gadolinium rods in the whole core should be less than or equal to 8 groups.

[0113] In the embodiment of the present application, if the target parameter is greater than the first preset parameter, the use of burnable poison is prohibited to control the quadrant power tilt of the core; if the target parameter is less than the second preset parameter, control the quantity of burnable poison used to control the quadrant power tilt of the core, so as to control the increase in quadrant power tilt.

[0114] Refer to Figure 11 , Figure 11 is a schematic flowchart of a core design method for controlling the increase in quadrant power tilt during the operation of a pressurized water reactor nuclear power plant provided by the embodiment of the present application. The method includes the following steps:

[0115] S1101, Divide the first change amount of the quadrant power tilt by the burnup increment to obtain the second expression of the second change rate of the quadrant power tilt.

[0116] S1102, Divide the second change amount of the quadrant power tilt by the burnup increment to obtain the third expression of the third change rate of the quadrant power tilt.

[0117] S1103, Add the second expression and the third expression to obtain the first expression.

[0118] S1104, According to the relationship between the quadrant power of each quadrant and the burnup increment, determine the fourth expression of the fourth change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poison.

[0119] S1105, Determine the similarity between the fourth expression and the second expression.

[0120] S1106, When the similarity is greater than the preset similarity threshold, use the slope in the third expression as the target parameter.

[0121] S1107, Control the quadrant power tilt of the core based on the target parameter.

[0122] It should be understood that although the steps in the flowcharts involved in the above-described embodiments are sequentially shown in the direction of the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless there is a clear description in this article, there is no strict order limit for the execution of these steps, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-described embodiments may include multiple steps or multiple stages. These steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be executed alternately or in turn with at least a part of other steps or steps or stages in other steps.

[0123] Based on the same inventive concept, an embodiment of the present application also provides a device for controlling the quadrant power tilt of a reactor core for implementing the method for controlling the quadrant power tilt of the reactor core involved above. The solution provided by this device to solve the problem is similar to the solution described in the above method. Therefore, the specific limitations in one or more embodiments of the device for controlling the quadrant power tilt of the reactor core provided below can refer to the limitations on the method for controlling the quadrant power tilt of the reactor core in the above text, and will not be repeated here.

[0124] In one embodiment, as Figure 12 shown, Figure 12 FIG. is a structural block diagram of a device for controlling the quadrant power tilt of a reactor core provided by an embodiment of the present application. The device 1200 includes:

[0125] A first determination module 1201, configured to determine a first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poisons according to the first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the reactor core, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core.

[0126] A second determination module 1202, configured to determine the target parameter of the burnable poison that affects the quadrant power tilt of the reactor core according to the first expression of the first change rate of the quadrant power tilt.

[0127] A control module 1203, configured to control the quadrant power tilt of the reactor core according to the target parameter of the burnable poison.

[0128] In one of the embodiments, the first determination module 1201 includes:

[0129] A first calculation unit, configured to divide the first change amount of the quadrant power tilt by the burnup increment to obtain a second expression of the second change rate of the quadrant power tilt.

[0130] A second calculation unit, configured to divide a second change amount of the quadrant power tilt by a burnup increment to obtain a third expression of a third change rate of the quadrant power tilt.

[0131] A first determination unit, configured to determine a first expression according to the second expression and the third expression.

[0132] In one embodiment, the first determination unit is specifically configured to add the second expression and the third expression to obtain the first expression.

[0133] In one embodiment, the apparatus 1200 further includes:

[0134] A third determination module, configured to determine a fourth expression of a fourth change rate of the quadrant power tilt corresponding to each quadrant in a fuel cycle without burnable poisons according to the relationship between the quadrant power of each quadrant and the burnup increment.

[0135] Correspondingly, the second determination module 1202 is specifically configured to determine a target parameter of the burnable poison affecting the quadrant power tilt of the reactor core according to the fourth expression and the first expression.

[0136] In one embodiment, the second determination module 1202 includes:

[0137] A second determination unit, configured to determine the similarity between the fourth expression and the second expression;

[0138] A third determination unit, configured to determine the target parameter according to the third expression of the third change rate of the quadrant power tilt when the similarity is greater than a preset similarity threshold.

[0139] In one embodiment, the apparatus 1200 further includes:

[0140] A fourth determination module, configured to, for each quadrant, determine a product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power of the burnable poison assembly containing the burnable poison in the quadrant, the burnup increment, and the slope corresponding to the burnable poison assembly containing y roots; the slope is the slope of the curve segment where the reactivity increases as the burnup increases in the K-infinity versus burnup curve graph, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant.

[0141] A fifth determination module, configured to determine a summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant.

[0142] A sixth determination module, configured to determine a second change amount of the quadrant power tilt according to the summation result and the quadrant power of the quadrant.

[0143] Correspondingly, the third determination unit is specifically configured to determine the target parameter according to the slope in the third expression.

[0144] In one embodiment, the control module 1203 is specifically configured to prohibit the use of burnable poisons if the target parameter is greater than a first preset parameter, so as to control the quadrant power tilt of the reactor core; if the target parameter is less than a second preset parameter, control the quantity of burnable poisons used, so as to control the quadrant power tilt of the reactor core.

[0145] Each module in the above quadrant power tilt control device of the reactor core can be implemented in whole or in part by software, hardware, and their combination. The above modules can be embedded in the processor of the computer device in the form of hardware or be independent of it, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above respective modules.

[0146] In one embodiment, a computer device is provided, including a memory and a processor. A computer program is stored in the memory. When the processor executes the computer program, the following steps are implemented:

[0147] Determine a first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poisons according to the first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the reactor core, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core;

[0148] Determine the target parameter of the burnable poison that affects the quadrant power tilt of the reactor core according to the first expression of the first change rate of the quadrant power tilt;

[0149] Control the quadrant power tilt of the reactor core according to the target parameter of the burnable poison.

[0150] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0151] Divide the first change amount of the quadrant power tilt by the burnup increment to obtain a second expression of the second change rate of the quadrant power tilt;

[0152] Divide the second change amount of the quadrant power tilt by the burnup increment to obtain a third expression of the third change rate of the quadrant power tilt;

[0153] Determine the first expression according to the second expression and the third expression.

[0154] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0155] Add the second expression and the third expression to obtain the first expression.

[0156] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0157] Determine a fourth expression for the fourth rate of change of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poisons according to the relationship between the quadrant power and the burnup increment of each quadrant;

[0158] Determine the target parameters affecting the quadrant power tilt of the reactor core according to the first expression for the first rate of change of the quadrant power tilt, including:

[0159] Determine the target parameters of the burnable poisons affecting the quadrant power tilt of the reactor core according to the fourth expression and the first expression.

[0160] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0161] Determine the similarity between the fourth expression and the second expression;

[0162] In the case where the similarity is greater than the preset similarity threshold, determine the target parameters according to the third expression for the third rate of change of the quadrant power tilt.

[0163] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0164] For each quadrant, determine the product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power of the burnable poison assembly with burnable poisons in the quadrant, the burnup increment, and the slope corresponding to the burnable poison assembly with y roots; the slope is the slope of the curve segment where the reactivity increases with the increase of burnup in the K-infinity vs. burnup curve, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant;

[0165] Determine the summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant;

[0166] Determine the second change amount of the quadrant power tilt according to the summation result and the quadrant power of the quadrant;

[0167] Determine the target parameters according to the third expression for the third rate of change of the quadrant power tilt, including:

[0168] Determine the target parameters according to the slope in the third expression.

[0169] In one embodiment, when the processor executes the computer program, the following steps are further implemented:

[0170] If the target parameter is greater than the first preset parameter, prohibit the use of burnable poisons to control the quadrant power tilt of the reactor core;

[0171] If the target parameter is less than the second preset parameter, control the quantity of the burnable poison used to control the quadrant power tilt of the reactor core.

[0172] In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented:

[0173] According to the first change amount of the quadrant power tilt caused by fuel combustion in the fuel cycle with burnable poison in each quadrant of the reactor core, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core, determine the first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poison;

[0174] According to the first expression of the first change rate of the quadrant power tilt, determine the target parameter of the burnable poison that affects the quadrant power tilt of the reactor core;

[0175] According to the target parameter of the burnable poison, control the quadrant power tilt of the reactor core.

[0176] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0177] Divide the first change amount of the quadrant power tilt by the burnup increment to obtain the second expression of the second change rate of the quadrant power tilt;

[0178] Divide the second change amount of the quadrant power tilt by the burnup increment to obtain the third expression of the third change rate of the quadrant power tilt;

[0179] Determine the first expression according to the second expression and the third expression.

[0180] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0181] Add the second expression and the third expression to obtain the first expression.

[0182] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0183] According to the relationship between the quadrant power of each quadrant and the burnup increment, determine the fourth expression of the fourth change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poison;

[0184] According to the first expression of the first change rate of the quadrant power tilt, determine the target parameter that affects the quadrant power tilt of the reactor core, including:

[0185] According to the fourth expression and the first expression, determine the target parameter of the burnable poison that affects the quadrant power tilt of the reactor core.

[0186] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0187] Determine the similarity between the fourth expression and the second expression;

[0188] When the similarity is greater than a preset similarity threshold, determine the target parameter according to the third expression of the third change rate of the quadrant power tilt.

[0189] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0190] According to the summation result and the quadrant power of the quadrant, determine the third change rate of the quadrant power tilt. For each quadrant, determine the product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power of the burnable poison assembly containing burnable poison in the quadrant, the burnup increment, and the slope corresponding to the burnable poison assembly containing y roots; the slope is the slope of the curve segment where the reactivity increases with the increase of burnup in the K-infinity vs. burnup curve graph, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant;

[0191] Determine the summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant;

[0192] Determine the target parameter according to the third expression of the third change rate of the quadrant power tilt, including:

[0193] Determine the target parameter according to the slope in the third expression.

[0194] In one embodiment, when the computer program is executed by a processor, the following steps are further implemented:

[0195] If the target parameter is greater than the first preset parameter, prohibit the use of burnable poison to control the quadrant power tilt of the reactor core;

[0196] If the said target parameter is less than the second preset parameter, control the quantity of burnable poison used to control the quadrant power tilt of the reactor core.

[0197] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memory can include Read-Only Memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc. The databases involved in the embodiments provided in the present application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in the present application can be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, data processing logics based on quantum computing, etc., without limitation.

[0198] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope recorded in this specification.

[0199] The above-described embodiments only represent several implementation manners of the present application. Their descriptions are relatively specific and detailed, but they should not be construed as limiting the scope of the patent of the present application. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the appended claims.

Claims

1. A method for controlling the quadrant power tilt of a reactor core, characterized in that The method includes: Determining a first expression of a first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poison according to a first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the core, a second change amount of the quadrant power tilt caused by burnable poison combustion, and a burnup increment of the core; Determining a target parameter of the burnable poison affecting the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt; Controlling the quadrant power tilt of the core according to the target parameter of the burnable poison; Wherein, the determining a first expression of a first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle with burnable poison according to a first change amount of the quadrant power tilt caused by fuel combustion in each quadrant of the core, a second change amount of the quadrant power tilt caused by burnable poison combustion, and a burnup increment of the core includes: Dividing the first change amount of the quadrant power tilt by the burnup increment to obtain a second expression of a second change rate of the quadrant power tilt; Dividing the second change amount of the quadrant power tilt by the burnup increment to obtain a third expression of a third change rate of the quadrant power tilt; Determining the first expression according to the second expression and the third expression.

2. The method according to claim 1, characterized in that, The determining the first expression according to the second expression and the third expression includes: Adding the second expression and the third expression to obtain the first expression.

3. The method according to claim 1 or 2, characterized in that, The method further includes: Determining a fourth expression of a fourth change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle without burnable poison according to the relationship between the quadrant power of each quadrant and the burnup increment; the relationship between the quadrant power of the quadrant and the burnup increment is: the burnup increment is the product of the quadrant power of the quadrant and the time increment; the fourth change rate is the ratio between the change amount of the quadrant power tilt corresponding to the fuel cycle without burnable poison and the burnup increment; The determining a target parameter of the burnable poison affecting the quadrant power tilt of the core according to the first expression of the first change rate of the quadrant power tilt includes: Determining the target parameter of the burnable poison affecting the quadrant power tilt of the core according to the fourth expression and the first expression.

4. The method according to claim 3, wherein The determining the target parameter of the burnable poison affecting the quadrant power tilt of the core according to the fourth expression and the first expression includes: Determining the similarity between the fourth expression and the second expression; When the similarity is greater than a preset similarity threshold, determining the target parameter according to the third expression of the third change rate of the quadrant power tilt.

5. The method according to claim 4, wherein The method further includes: For each of the quadrants, determine the product result corresponding to the m-th burnable poison assembly in the quadrant; the product result is the product of the power of the burnable poison assemblies containing burnable poison in the quadrant, the burnup increment, and the slope corresponding to the burnable poison assembly containing y roots; the slope is the slope of the curve segment where the reactivity increases with the increase of burnup in the K-infinity versus burnup curve, and m is an integer greater than or equal to 1 and less than or equal to the total number of burnable poison assemblies in the quadrant; Determine the summation result according to the product result corresponding to the m-th burnable poison assembly in the quadrant; Determine the second change amount of the quadrant power tilt according to the summation result and the quadrant power of the quadrant; The determining the target parameter according to the third expression of the third change rate of the quadrant power tilt includes: Determine the target parameter according to the slope in the third expression; 6. The method according to any one of claims 1-2, characterized in that The controlling the quadrant power tilt of the reactor core according to the target parameter of the burnable poison includes: If the target parameter is greater than the first preset parameter, prohibit the use of the burnable poison to control the quadrant power tilt of the reactor core; If the target parameter is less than the second preset parameter, control the quantity of the burnable poison used to control the quadrant power tilt of the reactor core.

7. A quadrant power tilt control device for a reactor core, characterized in that, The device includes: A first determination module, configured to determine a first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core in each quadrant of the reactor core; A second determination module, configured to determine the target parameter of the burnable poison that affects the quadrant power tilt of the reactor core according to the first expression of the first change rate of the quadrant power tilt; A control module, configured to control the quadrant power tilt of the reactor core according to the target parameter of the burnable poison; Wherein, the determining the first expression of the first change rate of the quadrant power tilt corresponding to each quadrant in the fuel cycle containing burnable poison according to the first change amount of the quadrant power tilt caused by fuel combustion, the second change amount of the quadrant power tilt caused by burnable poison combustion, and the burnup increment of the reactor core in each quadrant of the reactor core includes: Dividing the first change amount of the quadrant power tilt by the burnup increment to obtain a second expression of the second change rate of the quadrant power tilt; Dividing the second change amount of the quadrant power tilt by the burnup increment to obtain a third expression of the third change rate of the quadrant power tilt; Determine the first expression according to the second expression and the third expression.

8. A computer device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Citation Information

Patent Citations

  • Pressurized water reactor core long-cycle refueling and loading method based on gadolinium enrichment

    CN112420223A

  • Method for inhibiting quadrantal power inclination of pressurized water reactor nuclear power station

    CN1783354A