Method and device for controlling axial power distribution of pressurized water reactor, and electronic equipment

Abstract

The application discloses a method and device for controlling axial power distribution of a pressurized water reactor and electronic equipment. The method comprises the following steps: generating an initial population comprising multiple individuals, the individuals being rod position sequences of temperature control rods, the rod position sequences comprising multiple time nodes and multiple rod positions, and the multiple time nodes corresponding to the multiple rod positions; taking the minimum axial power offset of the pressurized water reactor as an optimization target, and determining a fitness evaluation function; taking the initial population as a parent population, performing iterative processing for finding an optimal solution based on the fitness evaluation function, and obtaining a target rod position sequence associated with the initial population; and controlling the rod positions of the temperature control rods according to the target rod position sequence, so as to control the axial power distribution of the pressurized water reactor.

Description

Technical Field

[0001] This application relates to the fields of reactor core design and safety technology, and in particular to a method, apparatus and electronic equipment for controlling the axial power distribution of a pressurized water reactor. Background Technology

[0002] The reactor core is the central part of the reactor where the fission reaction occurs. During the power rise and fall of a pressurized water reactor, the axial power distribution of the core is prone to significant changes, even oscillations, due to effects such as xenon distribution and moderator density. Related technologies control the position of temperature control rods to regulate the axial power distribution of the core, thereby ensuring its stability. However, this method of controlling the position of temperature control rods suffers from poor control over the axial power distribution of the core. Summary of the Invention

[0003] This application aims to propose a control method, device, and electronic equipment for the axial power distribution of a pressurized water reactor, which can solve the problem that traditional core axial power distribution control technology has poor control effect on the core axial power distribution.

[0004] In a first aspect, embodiments of this application provide a method for controlling the axial power distribution of a pressurized water reactor, including:

[0005] An initial population is generated, comprising multiple individuals, wherein the individuals are a sequence of temperature control rod positions, and the sequence of rod positions includes multiple time points and multiple rod positions, wherein the multiple time points correspond to the multiple rod positions;

[0006] The fitness evaluation function is determined with the minimum axial power offset of the pressurized water reactor as the optimization objective.

[0007] Using the initial population as the parent population, an iterative process for finding the optimal solution is performed based on the fitness evaluation function to obtain the target rod sequence associated with the initial population;

[0008] The position of the temperature control rod is controlled according to the target rod position sequence to control the axial power distribution of the pressurized water reactor.

[0009] According to some embodiments of this application, the iterative process of finding the optimal solution based on the fitness evaluation function includes multiple iterative processes, and the k-th iterative process includes the following steps:

[0010] The fitness value of each individual in the target parent population is calculated based on the fitness evaluation function. The target parent population includes at least one of the individuals. If k is 1, the target parent population is the initial population. If k is an integer greater than 1, the target parent population is the parent population obtained after the previous iteration.

[0011] Based on the fitness values ​​of each individual in the target parent population, selection, crossover, and mutation operations are performed on the target parent population to obtain the iteratively processed parent population.

[0012] If the parent population after iterative processing does not meet the preset iteration stopping condition, then the target parent population is updated to the parent population after iterative processing, and the above steps are repeated until the parent population after iterative processing meets the preset iteration stopping condition.

[0013] The target bar sequence is the parent population that satisfies the preset iteration stopping condition.

[0014] According to some embodiments of this application, generating an initial population comprising multiple individuals includes:

[0015] Obtain a feasible solution space, which includes multiple rod position sequences, wherein the rod position sequences satisfy the insertion limit of the temperature control rod and the removal limit of the temperature control rod;

[0016] Select one of the bar positions from the feasible solution space as the transition sequence;

[0017] The initial population is obtained by applying a random number of steps to multiple bars in the transition sequence.

[0018] According to some embodiments of this application, determining the fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization objective includes:

[0019] An axial power offset function is established to calculate the difference in axial power offset at different time nodes. The difference in axial power offset is calculated from a first power offset and a second power offset corresponding to the time node. The first power offset is a reference offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset is the actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node.

[0020] A rod position fitness function is established based on the axial power offset function. The rod position fitness function is used to calculate the fitness value of the rod position based on the difference in axial power offset.

[0021] A fitness evaluation function is established based on the aforementioned stick position fitness function. The fitness evaluation function is used to calculate the fitness value of the stick position sequence based on the fitness values ​​of multiple stick positions in the stick position sequence.

[0022] According to some embodiments of this application, the axial power offset difference is calculated through the following steps:

[0023] Subtract the second power offset from the first power offset to obtain the transition value of the axial power offset;

[0024] The absolute value of the axial power offset transition value is taken to obtain the axial power offset difference value.

[0025] According to some embodiments of this application, the first power offset is calculated through the following steps:

[0026] The first power offset is obtained by multiplying the full power offset by the node relative power. The full power offset is the axial power offset when the pressurized water reactor is at full power, and the node relative power is the relative power of the pressurized water reactor at the time node.

[0027] According to some embodiments of this application, the fitness value of the rod position is obtained through the following steps:

[0028] If the difference in axial power offset corresponding to the rod position is less than a preset first offset difference, the fitness value of the rod position is confirmed to be the preset first fitness value.

[0029] If the difference in axial power offset corresponding to the rod position is greater than or equal to the preset first offset difference, and the difference in axial power offset corresponding to the rod position is less than or equal to the preset second offset difference, the fitness value of the rod position is confirmed to be the preset second fitness value.

[0030] If the difference in axial power offset corresponding to the rod position is greater than the preset second offset difference, the fitness value of the rod position is confirmed to be the preset third fitness value.

[0031] According to some embodiments of this application, the fitness value of the bar sequence is calculated through the following steps:

[0032] The fitness values ​​of each stick in the stick sequence are summed to obtain the overall fitness value of the stick.

[0033] Based on the axial power offset difference values ​​corresponding to multiple rod positions in the rod position sequence, the largest axial power offset difference value is selected.

[0034] The maximum difference in axial power offset is multiplied by a preset correction coefficient to obtain the correction value, wherein the preset correction coefficient is a positive number;

[0035] The fitness value of the rod sequence is obtained by subtracting the correction value from the overall fitness value of the rod position.

[0036] Secondly, embodiments of this application provide a control device, including:

[0037] A population generation module is used to generate an initial population consisting of multiple individuals, wherein the individuals are a sequence of rod positions of a temperature control rod, and the sequence of rod positions includes multiple time nodes and multiple rod positions, wherein the multiple time nodes correspond to the multiple rod positions.

[0038] The function determination module is used to determine the fitness evaluation function with the goal of minimizing the axial power offset of the pressurized water reactor.

[0039] The iterative processing module is used to perform iterative processing to find the optimal solution based on the fitness evaluation function, with the initial population as the target parent population, to obtain the target rod sequence associated with the initial population;

[0040] The rod position control module is used to control the position of the temperature control rod according to the target rod position sequence, so as to control the axial power distribution of the pressurized water reactor.

[0041] Thirdly, embodiments of this application provide an electronic device, including:

[0042] At least one processor;

[0043] At least one memory for storing at least one program;

[0044] When at least one of the programs is executed by at least one of the processors, the control method for the axial power distribution of a pressurized water reactor as described above is implemented.

[0045] Fourthly, embodiments of this application provide a computer-readable storage medium storing a processor-executable program, which, when executed by a processor, is used to implement the control method for axial power distribution of a pressurized water reactor as described above.

[0046] In this embodiment, by performing iterative processing on the initial population, the target rod position sequence can be obtained. Since the fitness evaluation function in the iterative processing aims to minimize the axial power offset of the pressurized water reactor, the optimal solution is found using the fitness evaluation function to obtain the target rod position sequence. Controlling the position of the temperature control rods according to the target rod position sequence can minimize the axial power offset of the pressurized water reactor. Compared with traditional core axial power distribution control technology, this method has better control effect on the core axial power distribution and improves the stability of the core axial power distribution.

[0047] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0048] The present application will be further described below with reference to the accompanying drawings and embodiments, wherein:

[0049] Figure 1 A flowchart illustrating an embodiment of the control method for axial power distribution in a pressurized water reactor provided in this application;

[0050] Figure 2 A flowchart of the iterative process in an embodiment of the pressurized water reactor axial power distribution control method provided in this application;

[0051] Figure 3 A flowchart illustrating the generation of the initial population in an embodiment of the pressurized water reactor axial power distribution control method provided in this application;

[0052] Figure 4 A flowchart for determining the fitness evaluation function in an embodiment of the control method for axial power distribution of pressurized water reactors provided in this application;

[0053] Figure 5 A flowchart for calculating the axial power offset difference in an embodiment of the control method for axial power distribution of pressurized water reactor provided in this application;

[0054] Figure 6 A flowchart illustrating the calculation of the fitness value of the rod position in an embodiment of the control method for axial power distribution of a pressurized water reactor provided in this application;

[0055] Figure 7 A flowchart illustrating the calculation of the fitness value of the rod position sequence in an embodiment of the control method for axial power distribution of a pressurized water reactor provided in this application;

[0056] Figure 8 A schematic diagram of the structure of an embodiment of the control device provided in this application;

[0057] Figure 9 A schematic diagram of the structure of the electronic device provided in this application.

[0058] Figure label:

[0059] Control device 100, population generation module 110, function determination module 120, iterative processing module 130, rod position control module 140

[0060] Electronic device 200, processor 210, memory 220. Detailed Implementation

[0061] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain this application, and should not be construed as limiting this application.

[0062] In the description of this application, it should be understood that the orientation descriptions, such as up, down, etc., are based on the orientation or positional relationship shown in the accompanying drawings, and are only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of this application.

[0063] In the description of this application, "multiple" refers to two or more. The use of "first" and "second" is for the purpose of distinguishing technical features only and should not be construed as indicating or implying relative importance, or implicitly indicating the number of technical features indicated, or the order in which the technical features are indicated.

[0064] In the description of this application, unless otherwise expressly defined, terms such as "setup," "installation," and "connection" should be interpreted broadly, and those skilled in the art can reasonably determine the specific meaning of the above terms in this application in conjunction with the specific content of the technical solution.

[0065] The following reference Figures 1 to 9 This application describes a method, apparatus, and electronic device for controlling the axial power distribution of a pressurized water reactor according to embodiments of the present application.

[0066] The control method for axial power distribution of pressurized water reactors according to the embodiments of this application is applied to the axial power distribution control of pressurized water reactors, such as... Figure 1 As shown, it includes the following steps:

[0067] Step S100: Generate an initial population consisting of multiple individuals, where each individual is a rod position sequence of a temperature control rod. The rod position sequence includes multiple time points and multiple rod positions, with multiple time points corresponding to multiple rod positions.

[0068] Step S200: Determine the fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization objective;

[0069] Step S300: Using the initial population as the parent population, perform iterative processing to find the optimal solution based on the fitness evaluation function to obtain the target rod sequence associated with the initial population;

[0070] Step S400: Control the position of the temperature control rods according to the target rod position sequence to control the axial power distribution of the pressurized water reactor.

[0071] In this embodiment, by performing iterative processing on the initial population, the target rod position sequence can be obtained. Since the fitness evaluation function in the iterative processing aims to minimize the axial power offset of the pressurized water reactor (PWR), the optimal solution is found using the fitness evaluation function to obtain the target rod position sequence. Controlling the position of the temperature control rods according to the target rod position sequence can minimize the axial power offset of the PWR. Compared with traditional core axial power distribution control technology, the control method for PWR axial power distribution in this embodiment has a better control effect on the core axial power distribution and improves the stability of the core axial power distribution.

[0072] In step S100 above, an initial population including multiple rod positions is generated so that subsequent steps can obtain the target rod position sequence based on the initial population.

[0073] For example, a subset is selected from the feasible solution space and used as the initial population. The feasible solution space includes multiple rod position sequences, which satisfy the insertion limit of the temperature control rod and the removal limit of the temperature control rod.

[0074] As can be understood, the rod position sequence includes multiple time nodes and multiple rod positions, with each time node corresponding one-to-one with a specific rod position. The rod position sequence is used to control the position of the temperature control rod; that is, according to the rod position sequence, the temperature control rod is moved to the corresponding position at each time node.

[0075] In step S200 above, the axial power offset of the pressurized water reactor (PWR) reflects the stability of its axial power distribution. The larger the axial power offset, the more unstable the axial power distribution. Taking the minimum axial power offset as the optimization objective, a fitness evaluation function is determined. This function can be used to evaluate the control effect of the rod position sequence on the core axial power distribution.

[0076] For example, the fitness value of the rod sequence can be calculated using a fitness evaluation function. The larger the fitness value of the rod sequence, the better the control effect of the rod sequence on the axial power distribution of the core, and the better the stability of the axial power distribution of the pressurized water reactor.

[0077] In step S300 above, the initial population is iteratively processed based on the fitness evaluation function to find the optimal solution. The optimal solution is the target rod position sequence associated with the initial population. The target rod position sequence can minimize the axial power offset of the pressurized water reactor.

[0078] In step S400 above, the temperature control rods are controlled at multiple time points according to the target rod position sequence, thereby controlling the axial power distribution of the pressurized water reactor.

[0079] One embodiment of this application, such as Figure 2As shown, the iterative process for finding the optimal solution based on the fitness evaluation function in step S300 includes multiple iterative processes, and the k-th iterative process includes the following steps:

[0080] Step S310: Calculate the fitness value of each individual in the target parent population based on the fitness evaluation function. The target parent population includes at least one individual. If k is 1, the target parent population is the initial population. If k is an integer greater than 1, the target parent population is the parent population obtained after the previous iteration.

[0081] Step S320: Based on the fitness values ​​of each individual in the target parent population, perform selection, crossover, and mutation operations on the target parent population to obtain the iteratively processed parent population;

[0082] Step S330: If the parent population after iterative processing does not meet the preset iteration stopping condition, update the target parent population to the parent population after iterative processing, and repeat the above steps until the parent population after iterative processing meets the preset iteration stopping condition.

[0083] The target bar sequence is the parent population that meets the preset iteration stopping condition.

[0084] In this embodiment, the initial population is iteratively processed. Each iteration includes: selecting multiple bar sequences with high fitness values ​​from the target parent population through a selection operation; performing a crossover operation on the selected bar sequences, i.e., pairing the multiple bar sequences; crossovering some bars in the paired bar sequences based on a certain probability; and changing the bar values ​​of the crossover bar sequences with a certain probability, thereby obtaining the iteratively processed parent population. If the iteratively processed parent population does not meet the preset iteration stopping condition, the target parent population is updated to the iteratively processed parent population, and steps S310 and S320 are re-executed until the iteratively processed parent population meets the preset iteration stopping condition. The bar sequences in the parent population that meet the preset iteration stopping condition are the target bar sequences.

[0085] It should be noted that for the selection operation, multiple target individuals can be selected from the target parent population according to a preset percentage or a preset number. For the crossover and mutation operations, there is no restriction on the order of the crossover and mutation operations; the target individual's position can be crossovered first, followed by mutation, or vice versa. It is not required that every target individual undergo both crossover and mutation operations; that is, a target individual can undergo only crossover or only mutation operations.

[0086] It is understandable that the preset iteration stopping condition can be that the iteration process stops after satisfying the preset number of iterations, or it can be that the fitness value of each individual in the parent population is greater than the preset individual fitness value.

[0087] In step S320 above, one of the following selection algorithms can be used for the selection operation: roulette wheel selection, direct truncation selection, or elite strategy selection. One of the following crossover algorithms can be used: single-point crossover, multi-point crossover, or uniform crossover. One of the following mutation algorithms can be used: uniform mutation, non-uniform mutation, or multi-point mutation.

[0088] One embodiment of this application further illustrates the generation of an initial population comprising multiple individuals in step S100, such as... Figure 3 As shown, step S100 includes the following steps:

[0089] Step S110: Obtain the feasible solution space, which includes multiple rod position sequences. The rod position sequences satisfy the insertion limit of the temperature control rod and the removal limit of the temperature control rod.

[0090] Step S120: Select a bar sequence from the feasible solution space as the transition sequence;

[0091] Step S130: Obtain the initial population by applying random steps to multiple bars in the transition sequence.

[0092] In this embodiment, a feasible solution space is obtained, in which multiple rod position sequences are all feasible solutions, and all multiple rod position sequences satisfy the insertion and removal limits of the temperature control rod. A rod position sequence is selected from the feasible solution space as a transition sequence, and multiple rod positions in the transition sequence are subjected to random steps, thereby generating multiple random rod position sequences, which is the randomly generated initial population.

[0093] One embodiment of this application further explains the step S200, "determining the fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization objective," as follows: Figure 4 As shown, step S200 includes the following steps:

[0094] Step S210: Establish an axial power offset function. The axial power offset function is used to calculate the difference in axial power offset at different time nodes. The difference in axial power offset is calculated from the first power offset and the second power offset corresponding to the time node. The first power offset is the reference offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset is the actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node.

[0095] Step S220: Establish a rod position fitness function based on the axial power offset function. The rod position fitness function is used to calculate the fitness value of the rod position based on the difference in axial power offset.

[0096] Step S230: Establish a fitness evaluation function based on the stick position fitness function. The fitness evaluation function is used to calculate the fitness value of the stick position sequence based on the fitness values ​​of multiple stick positions in the stick position sequence.

[0097] In this embodiment, the difference between the first power offset and the second power offset corresponding to the same time node is calculated using the axial power offset function to obtain the axial power offset difference corresponding to the time node. Further, the fitness value of the rod position is obtained using the rod position fitness function based on the axial power offset difference. Finally, the fitness value of the rod position sequence is calculated using the fitness evaluation function based on the fitness values ​​of multiple rod positions in the rod position sequence.

[0098] One embodiment of this application, such as Figure 5 As shown, the axial power offset difference in step S210 is calculated through the following steps:

[0099] Step S211: Subtract the second power offset from the first power offset to obtain the transition value of the axial power offset;

[0100] Step S212: Take the absolute value of the axial power offset transition value to obtain the axial power offset difference.

[0101] In this embodiment, the absolute value of the difference between the first power offset and the second power offset is taken to obtain the axial power offset difference.

[0102] It is understandable that the difference in axial power offset is calculated using the following formula:

[0103] ΔI diff,n =|ΔI ref -ΔI n |,

[0104] Where, ΔI diff，n The difference in axial power offset, ΔI ref The first power offset, ΔI n denoted as the second power offset, and n as the time node.

[0105] In one embodiment of this application, the first power offset in step S210 is calculated through the following steps:

[0106] Step S213: Multiply the full power offset by the node relative power to obtain the first power offset. The full power offset is the axial power offset when the pressurized water reactor is at full power, and the node relative power is the relative power of the pressurized water reactor at the time node.

[0107] In this embodiment, the first power offset can be obtained by multiplying the full power offset by the node relative power.

[0108] It should be noted that the relative power at the node is the relative power of the pressurized water reactor at the time node, which is preset. Therefore, the obtained first power offset can be used as a reference value for comparison with the second power offset.

[0109] Understandably, the first power offset is calculated using the following formula:

[0110] ΔI ref =AO0·P n ,

[0111] Where AO0 is the axial power offset of the pressurized water reactor at full power, P n This represents the relative power of the pressurized water reactor at a given time point.

[0112] Understandably, the second power offset can be obtained using core calculation software such as SCIENCE based on the rod position sequence.

[0113] One embodiment of this application, such as Figure 6 As shown, the fitness value of the rod position in step S220 is obtained through the following steps:

[0114] Step S221: If the difference in axial power offset corresponding to the rod position is less than the preset first offset difference, confirm that the fitness value of the rod position is the preset first fitness value.

[0115] Step S222: If the difference in axial power offset corresponding to the rod position is greater than or equal to the preset first offset difference, and the difference in axial power offset is less than or equal to the preset second offset difference, the fitness value of the rod position is confirmed to be the preset second fitness value.

[0116] Step S223: If the difference in axial power offset corresponding to the rod position is greater than the preset second offset difference, confirm that the fitness value of the rod position is the preset third fitness value.

[0117] In this embodiment, by comparing the axial power offset difference with a preset first offset difference and a preset second offset difference, if the axial power offset difference is less than the preset first offset difference, the fitness value of the rod position is confirmed as the preset first fitness value; if the axial power offset difference is greater than or equal to the preset first offset difference and less than or equal to the preset second offset difference, the fitness value of the rod position is confirmed as the preset second fitness value; if the axial power offset difference is greater than the preset second offset difference, the fitness value of the rod position is confirmed as the preset third fitness value.

[0118] Understandably, the fitness value of a rod position is calculated using the following formula:

[0119]

[0120] Where f(x) is the fitness value of the rod position, and x is the difference in axial power offset.

[0121] It should be noted that when calculating the fitness value of a certain rod position, the difference in axial power offset at the time node corresponding to that rod position is calculated based on the first power offset and the second power offset at the time node. The fitness value of the rod position is then calculated based on the difference in axial power offset.

[0122] One embodiment of this application, such as Figure 7 As shown, the fitness value of the bar sequence in step S230 is calculated through the following steps:

[0123] Step S231: Sum the fitness values ​​of each stick in the stick sequence to obtain the overall fitness value of the sticks;

[0124] Step S232: Select the largest axial power offset difference based on the axial power offset difference between multiple rod positions in the rod position sequence;

[0125] Step S233: Multiply the maximum axial power offset difference by a preset correction coefficient to obtain the correction value. The preset correction coefficient is a positive number.

[0126] Step S234: Subtract the correction value from the overall fitness value of the rod position to obtain the fitness value of the rod position sequence.

[0127] In this embodiment, the overall fitness value of the rod position is obtained by summing the fitness values ​​of each rod position in the rod position sequence. However, summing only the fitness values ​​of each rod position may result in multiple rod position sequences having the same fitness value, making it difficult to guarantee the optimization effect of the iterative processing. Therefore, a correction value is used to correct the fitness value of the rod position sequence. Even if multiple rod position sequences have the same fitness value, the rod position sequence with the larger maximum value among the multiple axial power offset differences has a larger correction value, making the fitness value of the rod position sequence smaller, thereby ensuring the optimization effect of the iterative processing.

[0128] Understandably, the fitness value of the bar sequence is calculated using the following formula:

[0129]

[0130] Where η is the fitness value of the bar sequence, n is the time node, and f(ΔI) diff，n ) represents the fitness value of the rod position, k is a correction coefficient and k is a positive number, ΔI diffIt is the largest difference in axial power offset among the differences in axial power offset among multiple positions in the bar sequence.

[0131] Understandably, the correction factor can be any small positive number, for example, a correction factor of 0.0125.

[0132] The control method for axial power distribution of a pressurized water reactor provided in this application embodiment can be executed by a control device 100 for axial power distribution of a pressurized water reactor. This application embodiment uses the control device 100 for axial power distribution of a pressurized water reactor executing the control method as an example to illustrate the control device 100 for axial power distribution of a pressurized water reactor provided in this application embodiment.

[0133] Please see Figure 8 This is a schematic diagram of the structure of a pressurized water reactor axial power distribution control device 100 provided in an embodiment of this application. This pressurized water reactor axial power distribution control device 100 is applied to an electronic device 200. For example... Figure 8 As shown, the control device 100 includes:

[0134] The population generation module 110 is used to generate an initial population consisting of multiple individuals, where each individual is a rod position sequence of a temperature control rod. The rod position sequence includes multiple time points and multiple rod positions, with multiple time points corresponding to multiple rod positions.

[0135] The function determination module 120 is used to determine the fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization objective.

[0136] The iterative processing module 130 is used to perform iterative processing to find the optimal solution based on the fitness evaluation function with the initial population as the target parent population, and obtain the target rod sequence associated with the initial population.

[0137] The rod position control module 140 is used to control the position of the temperature control rods according to the target rod position sequence in order to control the axial power distribution of the pressurized water reactor.

[0138] In one embodiment of this application, the iterative processing module 130 includes:

[0139] The individual fitness calculation module is used to calculate the fitness value of each individual in the target parent population based on the fitness evaluation function.

[0140] The target parent population generation module is used to perform selection, crossover, and mutation operations on the target parent population based on the fitness values ​​of each individual in the target parent population to obtain a new target parent population.

[0141] In one embodiment of this application, the population generation module 110 includes:

[0142] The feasible solution acquisition module is used to acquire the feasible solution space, which includes multiple rod position sequences. The rod position sequences satisfy the insertion limit of the temperature control rod and the removal limit of the temperature control rod.

[0143] The first selection module is used to select a bar sequence from the feasible solution space as a transition sequence;

[0144] The step-applying module is used to obtain the initial population by applying random steps to multiple bars in the transition sequence.

[0145] In one embodiment of this application, the function determination module 120 includes:

[0146] The first function establishment module is used to establish the axial power offset function. The axial power offset function is used to calculate the difference between the axial power offset based on the first power offset and the second power offset. The first power offset is the reference offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset is the actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node.

[0147] The second function establishment module is used to establish a rod position fitness function based on the axial power offset function. The rod position fitness function is used to calculate the fitness value of the rod position based on the difference in axial power offset.

[0148] The third function establishment module is used to establish a fitness evaluation function based on the stick position fitness function. The fitness evaluation function is used to calculate the fitness value of the stick position sequence based on the fitness values ​​of multiple stick positions in the stick position sequence.

[0149] In one embodiment of this application, the first function creation module includes:

[0150] The transition value calculation module is used to subtract the second power offset from the first power offset to obtain the axial power offset transition value;

[0151] The absolute value calculation module is used to take the absolute value of the axial power offset transition value to obtain the axial power offset difference.

[0152] In one embodiment of this application, the transition value calculation module includes:

[0153] The first power offset calculation module is used to multiply the full power offset by the node relative power to obtain the first power offset. The full power offset is the axial power offset when the pressurized water reactor is at full power, and the node relative power is the relative power of the pressurized water reactor at the time node.

[0154] In one embodiment of this application, the second function creation module includes:

[0155] The first judgment module is used to confirm the fitness value of the rod position as the preset first fitness value if the difference in axial power offset is less than the preset first offset difference.

[0156] The second judgment module is used to confirm that the fitness value of the rod position is the preset second fitness value if the difference in axial power offset is greater than or equal to the preset first offset difference and the difference in axial power offset is less than or equal to the preset second offset difference.

[0157] The third judgment module is used to confirm the fitness value of the rod position as the preset third fitness value if the difference in axial power offset is greater than the preset second offset difference.

[0158] In one embodiment of this application, the third function creation module includes:

[0159] The summation module is used to sum the fitness values ​​of each stick in the stick sequence to obtain the overall fitness value of the sticks.

[0160] The second selection module is used to select the largest axial power offset difference based on the axial power offset difference between multiple rod positions in the rod position sequence.

[0161] The correction module is used to multiply the maximum axial power offset difference by a preset correction coefficient to obtain the correction value. The preset correction coefficient is a positive number.

[0162] The difference module is used to subtract the correction value from the overall fitness value of the rod position to obtain the fitness value of the rod position sequence.

[0163] The control device 100 provided in this application embodiment can achieve... Figures 1 to 7 The various processes implemented in the method implementation examples achieve the same beneficial effects, and will not be described again here to avoid repetition.

[0164] In addition, one embodiment of this application also discloses an electronic device 200, such as... Figure 9 As shown, it includes:

[0165] At least one processor 210;

[0166] At least one memory 220 is used to store at least one program;

[0167] When at least one program is executed by at least one processor 210, the control method for the axial power distribution of the pressurized water reactor as described above is implemented.

[0168] The electronic device 200 provided in this application embodiment can achieve... Figures 1 to 7 The various processes implemented in the method implementation examples achieve the same beneficial effects, and will not be described again here to avoid repetition.

[0169] In addition, one embodiment of this application discloses a computer-readable storage medium storing a processor-executable program, which, when executed by a processor, is used to implement the control method for axial power distribution of a pressurized water reactor as described above.

[0170] The computer-readable storage medium provided in the embodiments of this application can achieve... Figures 1 to 7 The various processes implemented in the method implementation examples achieve the same beneficial effects, and will not be described again here to avoid repetition.

[0171] It will be understood by those skilled in the art that all or some of the steps and systems in the methods disclosed above can be implemented as software, firmware, hardware, and suitable combinations thereof. Some or all of the physical components can be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application-specific integrated circuit. Such software can be distributed on a computer-readable medium, which can include computer storage media (or non-transitory media) and communication media (or transient media). As is known to those skilled in the art, the term computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program modules, or other data). Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technologies, CD-ROM, digital versatile disc (DVD) or other optical disc storage, magnetic cartridges, magnetic tape, disk storage or other magnetic storage devices, or any other medium that can be used to store desired information and is accessible to a computer. Furthermore, as is known to those skilled in the art, communication media typically contain computer-readable instructions, data structures, program modules, or other data in modulated data signals such as carrier waves or other transmission mechanisms, and may include any information delivery medium.

[0172] The embodiments of this application have been described in detail above with reference to the accompanying drawings. However, this application is not limited to the above embodiments. Within the scope of knowledge possessed by those skilled in the art, various changes can be made without departing from the spirit of this application.

Claims

1. A method for controlling the axial power distribution of a pressurized water reactor, characterized in that, include: An initial population is generated, comprising multiple individuals, wherein the individuals are a sequence of temperature control rod positions, and the sequence of rod positions includes multiple time points and multiple rod positions, wherein the multiple time points correspond to the multiple rod positions; The fitness evaluation function is determined with the minimum axial power offset of the pressurized water reactor as the optimization objective. Using the initial population as the parent population, an iterative process for finding the optimal solution is performed based on the fitness evaluation function to obtain the target rod sequence associated with the initial population; According to the target rod position sequence, the position of the temperature control rod is controlled to control the axial power distribution of the pressurized water reactor; The fitness evaluation function is determined with the goal of minimizing the axial power offset of the pressurized water reactor, including: An axial power offset function is established to calculate the difference in axial power offset at different time nodes. The difference in axial power offset is calculated from a first power offset and a second power offset corresponding to the time node. The first power offset is a reference offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset is the actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node. A rod position fitness function is established based on the axial power offset function. The rod position fitness function is used to calculate the fitness value of the rod position based on the difference in axial power offset. A fitness evaluation function is established based on the aforementioned bar position fitness function. This fitness evaluation function is used to calculate the fitness value of the bar position sequence based on the fitness values ​​of multiple bars in the sequence. The fitness value of the rod position is obtained through the following steps: If the difference in axial power offset corresponding to the rod position is less than a preset first offset difference, the fitness value of the rod position is confirmed to be the preset first fitness value. If the difference in axial power offset corresponding to the rod position is greater than or equal to the preset first offset difference, and the difference in axial power offset corresponding to the rod position is less than or equal to the preset second offset difference, the fitness value of the rod position is confirmed to be the preset second fitness value. If the difference in axial power offset corresponding to the rod position is greater than the preset second offset difference, the fitness value of the rod position is confirmed to be the preset third fitness value; or, The fitness values ​​of each stick in the stick sequence are summed to obtain the overall fitness value of the stick. Based on the axial power offset difference values ​​corresponding to multiple rod positions in the rod position sequence, the largest axial power offset difference value is selected. The maximum difference in axial power offset is multiplied by a preset correction coefficient to obtain the correction value, wherein the preset correction coefficient is a positive number; The fitness value of the rod sequence is obtained by subtracting the correction value from the overall fitness value of the rod position.

2. The method for controlling the axial power distribution of a pressurized water reactor according to claim 1, characterized in that, The iterative process of finding the optimal solution based on the fitness evaluation function includes multiple iterative processes, and the k-th iterative process includes the following steps: The fitness value of each individual in the target parent population is calculated based on the fitness evaluation function. The target parent population includes at least one of the individuals. If k is 1, the target parent population is the initial population. If k is an integer greater than 1, the target parent population is the parent population obtained after the previous iteration. Based on the fitness values ​​of each individual in the target parent population, selection, crossover, and mutation operations are performed on the target parent population to obtain the iteratively processed parent population. If the parent population after iterative processing does not meet the preset iteration stopping condition, then the target parent population is updated to the parent population after iterative processing, and the above steps are repeated until the parent population after iterative processing meets the preset iteration stopping condition. The target bar sequence is the parent population that satisfies the preset iteration stopping condition.

3. The method for controlling the axial power distribution of a pressurized water reactor according to claim 1, characterized in that, The generation of an initial population comprising multiple individuals includes: Obtain a feasible solution space, which includes multiple rod position sequences, wherein the rod position sequences satisfy the insertion limit of the temperature control rod and the removal limit of the temperature control rod; Select one of the bar positions from the feasible solution space as the transition sequence; The initial population is obtained by applying a random number of steps to multiple bars in the transition sequence.

4. The method for controlling the axial power distribution of a pressurized water reactor according to claim 1, characterized in that, The axial power offset difference is calculated through the following steps: Subtract the second power offset from the first power offset to obtain the transition value of the axial power offset; The absolute value of the axial power offset transition value is taken to obtain the axial power offset difference value.

5. The method for controlling the axial power distribution of a pressurized water reactor according to claim 1, characterized in that, The first power offset is calculated through the following steps: The first power offset is obtained by multiplying the full power offset by the node relative power. The full power offset is the axial power offset when the pressurized water reactor is at full power, and the node relative power is the relative power of the pressurized water reactor at the time node.

6. A control device, characterized in that, include: A population generation module is used to generate an initial population consisting of multiple individuals, wherein the individuals are a sequence of rod positions of a temperature control rod, and the sequence of rod positions includes multiple time nodes and multiple rod positions, wherein the multiple time nodes correspond to the multiple rod positions. The function determination module is used to determine the fitness evaluation function with the goal of minimizing the axial power offset of the pressurized water reactor. The iterative processing module is used to perform iterative processing to find the optimal solution based on the fitness evaluation function, with the initial population as the target parent population, to obtain the target rod sequence associated with the initial population; The rod position control module is used to control the position of the temperature control rod according to the target rod position sequence, so as to control the axial power distribution of the pressurized water reactor. The fitness evaluation function is determined with the goal of minimizing the axial power offset of the pressurized water reactor, including: An axial power offset function is established to calculate the difference in axial power offset at different time nodes. The difference in axial power offset is calculated from a first power offset and a second power offset corresponding to the time node. The first power offset is a reference offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset is the actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node. A rod position fitness function is established based on the axial power offset function. The rod position fitness function is used to calculate the fitness value of the rod position based on the difference in axial power offset. A fitness evaluation function is established based on the aforementioned bar position fitness function. This fitness evaluation function is used to calculate the fitness value of the bar position sequence based on the fitness values ​​of multiple bars in the sequence. The fitness value of the rod position is obtained through the following steps: If the difference in axial power offset corresponding to the rod position is less than a preset first offset difference, the fitness value of the rod position is confirmed to be the preset first fitness value. If the difference in axial power offset corresponding to the rod position is greater than or equal to the preset first offset difference, and the difference in axial power offset corresponding to the rod position is less than or equal to the preset second offset difference, the fitness value of the rod position is confirmed to be the preset second fitness value. If the difference in axial power offset corresponding to the rod position is greater than the preset second offset difference, the fitness value of the rod position is confirmed to be the preset third fitness value; or, The fitness values ​​of each stick in the stick sequence are summed to obtain the overall fitness value of the stick. Based on the axial power offset difference values ​​corresponding to multiple rod positions in the rod position sequence, the largest axial power offset difference value is selected. The maximum difference in axial power offset is multiplied by a preset correction coefficient to obtain the correction value, wherein the preset correction coefficient is a positive number; The fitness value of the rod sequence is obtained by subtracting the correction value from the overall fitness value of the rod position.

7. An electronic device, characterized in that, include: At least one processor; At least one memory for storing at least one program; When at least one of the programs is executed by at least one of the processors, the control method for axial power distribution of a pressurized water reactor as described in any one of claims 1 to 5 is implemented.

Images