Control method for axial power distribution of pressurized water reactor, and apparatus

Abstract

A control method for axial power distribution of a pressurized water reactor, and an apparatus. The control method comprises: generating an initial population comprising a plurality of individuals, each individual being a rod position sequence of temperature control rods, the rod position sequence comprising a plurality of time nodes and a plurality of rod positions, and the plurality of time nodes corresponding to the plurality of rod positions (S100); taking the minimum axial power offset of the pressurized water reactor as an optimization objective, so as to determine a fitness evaluation function (S200); taking the initial population as a parent population, and, on the basis of the fitness evaluation function, performing iterative processing for finding an optimal solution, so as to obtain a target rod position sequence associated with the initial population (S300); and controlling the rod positions of the temperature control rods on the basis of the target rod position sequence, so as to control the axial power distribution of the pressurized water reactor (S400).

Description

Method and device for controlling axial power distribution of pressurized water reactor TECHNICAL FIELD

[0001] The present application relates to the technical field of reactor core nuclear design and safety, in particular to a method and device for controlling axial power distribution of pressurized water reactor. BACKGROUND

[0002] The reactor core is the core part of the reactor where the fission reaction occurs. In the process of power rise and fall of the pressurized water reactor, the axial power distribution of the reactor core is prone to large changes or even oscillation due to the influence of xenon distribution and moderator density effects. In the related art, the axial power distribution of the reactor core is controlled by controlling the rod position of the temperature control rod, so as to ensure the stability of the axial power distribution of the reactor core. However, the method of controlling the rod position of the temperature control rod has the problem of poor control effect on the axial power distribution of the reactor core. SUMMARY

[0003] The present application aims to provide a method and device for controlling axial power distribution of pressurized water reactor, which can solve the problem of poor control effect on the axial power distribution of the reactor core in the traditional reactor core axial power distribution control technology.

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

[0005] generating an initial population comprising a plurality of individuals, the individuals being rod position sequences of temperature control rods, the rod position sequence comprising a plurality of time nodes and a plurality of rod positions, the plurality of time nodes corresponding to the plurality of rod positions;

[0006] determining a fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization target;

[0007] performing iterative processing of finding optimal solution based on the fitness evaluation function with the initial population as the parent population, to obtain a target rod position sequence associated with the initial population;

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

[0009] According to some embodiments of the present application, the iterative processing of finding optimal solution based on the fitness evaluation function comprises a plurality of iterative processing processes, and the kth iterative processing process comprises the following steps:

[0010] calculating fitness values of each individual in a target parent population based on the fitness evaluation function, the target parent population comprising at least one of the individuals, and if the k is 1, the target parent population is the initial population; if the k is an integer greater than 1, the target parent population is a parent population obtained after a previous iteration process;

[0011] performing selection operation, crossover operation and mutation operation on the target parent population based on the fitness values of each individual in the target parent population, to obtain a parent population after iteration processing;

[0012] if the parent population after iteration processing does not satisfy a preset iteration stop condition, updating the target parent population to the parent population after iteration processing, and re-executing the above steps until the parent population after iteration processing satisfies the preset iteration stop condition,

[0013] wherein the target rod position sequence is the parent population satisfying the preset iteration stop condition.

[0014] According to some embodiments of the present application, the initial population comprising a plurality of individuals is generated, comprising:

[0015] obtaining a feasible solution space, the feasible solution space comprising a plurality of rod position sequences, the rod position sequences satisfying the insertion limit of the temperature control rod, and the rod position sequences satisfying the extraction limit of the temperature control rod;

[0016] selecting one of the rod position sequences in the feasible solution space as a transition sequence;

[0017] obtaining the initial population by applying a random number of steps to a plurality of the rod positions in the transition sequence.

[0018] According to some embodiments of the present application, the fitness evaluation function is determined with the minimum axial power offset of the pressurized water reactor as the optimization target, comprising:

[0019] establishing an axial power offset function, the axial power offset function being used to calculate axial power offset difference values at different time nodes, wherein the axial power offset difference values are calculated from a first power offset and a second power offset corresponding to the time nodes, the first power offset being 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 being an actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node;

[0020] establishing a rod position fitness function based on the axial power offset function, the rod position fitness function being used to calculate fitness values of the rod positions according to the axial power offset difference values.

[0021] A fitness evaluation function is established based on the rod position fitness function, and the fitness evaluation function is used to calculate a fitness value of the rod position sequence according to fitness values of a plurality of rod positions in the rod position sequence.

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

[0023] The first power offset is subtracted from the second power offset to obtain an axial power offset transition value;

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

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

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

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

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

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

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

[0031] According to some embodiments of the present application, the fitness value of the rod position sequence is calculated by the following steps:

[0032] The fitness values of each rod position in the rod position sequence are summed to obtain a rod position fitness overall value;

[0033] According to the axial power offset difference values corresponding to a plurality of rod positions in the rod position sequence, the largest axial power offset difference value is selected;

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

[0035] The total value of the rod position fitness is reduced by the correction value to obtain a fitness value of the rod position sequence.

[0036] In a second aspect, an embodiment of the present application provides a control device, comprising:

[0037] A population generation module is configured to generate an initial population comprising a plurality of individuals, the individuals being rod position sequences of temperature control rods, the rod position sequence comprising a plurality of time nodes and a plurality of rod positions, the plurality of time nodes corresponding to the plurality of rod positions;

[0038] A function determination module is configured to determine a fitness evaluation function with the minimum axial power offset of a pressurized water reactor as an optimization target;

[0039] An iteration processing module is configured to take the initial population as a target parent population, perform iteration processing for finding an optimal solution based on the fitness evaluation function, and obtain a target rod position sequence associated with the initial population;

[0040] A rod position control module is configured to control the rod positions of the temperature control rods according to the target rod position sequence to control the axial power distribution of the pressurized water reactor.

[0041] In a third aspect, an embodiment of the present application provides an electronic device, comprising:

[0042] At least one processor;

[0043] At least one memory configured to store at least one program;

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

[0045] In a fourth aspect, an embodiment of the present application provides a computer readable storage medium, wherein a processor executable program is stored, and the processor executable program is executed by a processor to implement the control method for the axial power distribution of the pressurized water reactor.

[0046] In the embodiment of the present application, the target rod position sequence can be obtained by performing the iteration processing on the initial population. Since the fitness evaluation function is used to find the optimal solution to obtain the target rod position sequence, and the fitness evaluation function is used to minimize the axial power deviation of the pressurized water reactor in the iteration processing, the axial power deviation of the pressurized water reactor can be minimized by controlling the rod position of the temperature control rod according to the target rod position sequence. Compared with the traditional core axial power distribution control technology, the control effect of the core axial power distribution is good, and the stability of the core axial power distribution is improved.

[0047] Additional aspects and advantages of the present application will be made apparent from the following description. BRIEF DESCRIPTION OF DRAWINGS

[0048] The present application will be further described below in conjunction with the drawings and embodiments, in which:

[0049] FIG. 1 is a flowchart of an embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0050] FIG. 2 is a flowchart of the iteration processing in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0051] FIG. 3 is a flowchart of generating the initial population in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0052] FIG. 4 is a flowchart of determining the fitness evaluation function in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0053] FIG. 5 is a flowchart of calculating the axial power deviation difference in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0054] FIG. 6 is a flowchart of calculating the fitness value of the rod position in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0055] FIG. 7 is a flowchart of calculating the fitness value of the rod position sequence in the embodiment of the control method of the axial power distribution of the pressurized water reactor provided by the present application;

[0056] FIG. 8 is a structural schematic diagram of an embodiment of the control device provided by the present application;

[0057] FIG. 9 is a structural schematic diagram of an electronic device provided by the present application.

[0058] Reference signs: control device 100, population generation module 110, function determination module 120, iteration processing module 130, rod position control module 140, electronic device 200, processor 210, memory 220. DETAILED DESCRIPTION

[0059] Embodiments of the present application are described below in detail with reference to the accompanying drawings, wherein the same or similar components or components having the same or similar functions are denoted by the same or similar reference numerals throughout the drawings. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present application and cannot be understood as a limitation of the present application.

[0060] In the description of the present application, it should be understood that the orientation description, such as the orientation or position relationship indicated by up, down, etc. is based on the orientation or position relationship shown in the drawings, and is only for the purpose of facilitating the description of the present application and simplifying the description, and does not indicate or imply that the device or component referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation of the present application.

[0061] In the description of the present application, multiple refers to more than two. If there is a description of first, second, it is only for the purpose of distinguishing technical features, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated or the sequence of technical features indicated.

[0062] In the description of the present application, unless otherwise explicitly limited, the words such as arrangement, installation, connection, etc. should be broadly understood, and those skilled in the art can reasonably determine the specific meaning of the above words in the present application in combination with the specific content of the technical solution.

[0063] A method, device and electronic equipment for controlling axial power distribution of a pressurized water reactor are described below with reference to FIGS. 1-9 according to embodiments of the present application.

[0064] The method for controlling axial power distribution of a pressurized water reactor according to embodiments of the present application is applied to control the axial power distribution of a pressurized water reactor, as shown in FIG. 1, and includes the following steps:

[0065] Step S100: generating an initial population including a plurality of individuals, the individuals being rod position sequences of temperature control rods, the rod position sequence including a plurality of time nodes and a plurality of rod positions, the plurality of time nodes corresponding to the plurality of rod positions;

[0066] Step S200: determining a fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization target;

[0067] Step S300: taking the initial population as the parent population, performing iterative processing to find the optimal solution based on the fitness evaluation function, and obtaining a target rod position sequence associated with the initial population;

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

[0069] In the embodiment, the target rod position sequence can be obtained by performing the iterative processing on the initial population. Since the fitness evaluation function is used to minimize the axial power deviation of the pressurized water reactor in the iterative processing, the target rod position sequence is obtained by using the fitness evaluation function to find the optimal solution, and the rod position of the temperature control rod is controlled according to the target rod position sequence, so that the axial power deviation of the pressurized water reactor is minimized. Compared with the traditional core axial power distribution control technology, the control method of the axial power distribution of the pressurized water reactor has good control effect on the core axial power distribution and improves the stability of the core axial power distribution.

[0070] In the step S100, the initial population including a plurality of rod position sequences is generated, so that the target rod position sequence is obtained from the initial population in the subsequent steps.

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

[0072] It can be understood that the rod position sequence includes a plurality of time nodes and a plurality of rod positions, and the plurality of time nodes correspond to the plurality of rod positions one by one. The rod position sequence is used to control the rod position of the temperature control rod, that is, according to the rod position sequence, the temperature control rod is controlled to move to the corresponding rod position at each time node.

[0073] In the step S200, the axial power deviation of the pressurized water reactor can reflect the stability of the axial power distribution of the pressurized water reactor. The larger the axial power deviation of the pressurized water reactor is, the more unstable the axial power distribution is. The fitness evaluation function is determined by taking the minimum axial power deviation of the pressurized water reactor as the optimization target, and the control effect of the rod position sequence on the core axial power distribution can be evaluated by the fitness evaluation function.

[0074] For example, the fitness value of the rod position sequence is calculated by the fitness evaluation function. The larger the fitness value of the rod position sequence is, the better the control effect of the rod position sequence on the core axial power distribution is, and the better the stability of the axial power distribution of the pressurized water reactor is.

[0075] In the step S300, the initial population is processed iteratively based on the fitness evaluation function, and the optimal solution is found by the iterative processing. The optimal solution is the target rod position sequence associated with the initial population, and the target rod position sequence can minimize the axial power deviation of the pressurized water reactor.

[0076] In the step S400, the rod position of the temperature control rod at the plurality of time nodes is controlled according to the target rod position sequence, so that the axial power distribution of the pressurized water reactor is controlled.

[0077] In an embodiment of the present application, as shown in FIG. 2, the iteration process of finding the optimal solution based on the fitness evaluation function in step S300 includes multiple iteration processes, and the kth iteration process includes the following steps:

[0078] 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, and 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 last iteration process;

[0079] Step S320: Perform selection operation, crossover operation and mutation operation on the target parent population based on the fitness value of each individual in the target parent population to obtain the parent population after iteration;

[0080] Step S330: If the parent population after iteration does not meet the preset iteration stop condition, update the target parent population to the parent population after iteration, and re-execute the above steps until the parent population after iteration meets the preset iteration stop condition,

[0081] Wherein, the target rod position sequence is the parent population that meets the preset iteration stop condition.

[0082] In this embodiment, the initial population is iteratively processed, and each iteration process includes: selecting multiple rod position sequences with larger fitness values from the target parent population through selection operation, then performing crossover operation on the selected multiple rod position sequences, i.e. pairing the multiple rod position sequences, crossing part of the rod positions in the paired rod position sequences based on a certain probability, and changing the value of the rod position of the crossed rod position sequence with a certain probability, thereby obtaining the parent population after iteration. If the parent population after iteration does not meet the preset iteration stop condition, update the target parent population to the parent population after iteration, and re-execute the above steps S310 and S320 until the parent population after iteration meets the preset iteration stop condition, and the rod position sequence in the parent population that meets the preset iteration stop condition is the target rod position sequence.

[0083] It should be noted that for the selection operation, a certain percentage or a certain number of target individuals can be selected from the target parent population through the selection operation. For the crossover operation and the mutation operation, the order of the crossover operation and the mutation operation is not limited, the rod position of the target individual can be first subjected to the crossover operation and then subjected to the mutation operation, or the rod position of the target individual can be first subjected to the mutation operation and then subjected to the crossover operation. Each target individual does not have to perform the crossover operation and the mutation operation, i.e. the target individual can only perform the crossover operation or only perform the mutation operation.

[0084] It can be understood that the preset iteration stopping condition can be that the iteration process stops after a preset iteration number is met, or that the fitness values of each individual in the parent population are all greater than a preset individual fitness value.

[0085] In step S320, one of roulette selection, direct cutting method, and elite strategy method can be used for the selection operation. One of single-point crossover method, multi-point crossover method, and uniform crossover method can be selected for the crossover operation. A uniform mutation method, a non-uniform mutation method, and a multi-point mutation method can be selected for the mutation operation.

[0086] In an embodiment of the present application, the generating of the initial population including multiple individuals in step S100 is further described. As shown in FIG. 3, step S100 includes the following steps:

[0087] In step S110, a feasible solution space is obtained, the feasible solution space including multiple rod position sequences, the rod position sequences satisfying the insertion limit of the temperature control rod and the extraction limit of the temperature control rod;

[0088] In step S120, one rod position sequence is selected from the feasible solution space as a transition sequence.

[0089] In step S130, a random number of steps is applied to multiple rod positions in the transition sequence to obtain the initial population.

[0090] In the embodiment, the feasible solution space is obtained, and the multiple rod position sequences in the feasible solution space are all feasible solutions, and the multiple rod position sequences all satisfy the insertion limit and the extraction limit of the temperature control rod. One rod position sequence is selected from the feasible solution space as a transition sequence, and a random number of steps is applied to multiple rod positions in the transition sequence, thereby generating multiple random rod position sequences, that is, a randomly generated initial population is obtained.

[0091] In an embodiment of the present application, the "determining of the fitness evaluation function with the minimum axial power offset of the pressurized water reactor as the optimization target" in step S200 is further described. As shown in FIG. 4, step S200 includes the following steps:

[0092] In step S210, an axial power offset function is established, the axial power offset function being used to calculate axial power offset difference values at different time nodes, wherein the axial power offset difference values are calculated from a first power offset and a second power offset corresponding to the time nodes, the first power offset being 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 being an actual offset of the axial power of the pressurized water reactor when the temperature control rod is at the time node.

[0093] Step S220: establishing a rod position fitness function based on the axial power offset function, the rod position fitness function being used to calculate a fitness value of a rod position according to the axial power offset difference value;

[0094] Step S230: establishing a fitness evaluation function based on the rod position fitness function, the fitness evaluation function being used to calculate a fitness value of a rod position sequence according to the fitness values of a plurality of rod positions in the rod position sequence.

[0095] In this embodiment, the difference value between the first power offset and the second power offset corresponding to the same time node is calculated by using the axial power offset function to obtain the axial power offset difference value corresponding to the time node. Then, the fitness value of the rod position is obtained from the axial power offset difference value by using the rod position fitness function. Finally, the fitness value of the rod position sequence is calculated from the fitness values of a plurality of rod positions in the rod position sequence by using the fitness evaluation function.

[0096] In an embodiment of the present application, as shown in FIG. 5, the axial power offset difference value in step S210 is calculated by the following steps:

[0097] Step S211: subtracting the second power offset from the first power offset to obtain an axial power offset transition value;

[0098] Step S212: taking an absolute value of the axial power offset transition value to obtain the axial power offset difference value.

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

[0100] It can be understood that the axial power offset difference value is calculated by the following formula: diff,n = |ΔI ref - ΔI n |,

[0101] wherein, ΔI diff，n is the axial power offset difference value, ΔI ref is the first power offset, ΔI n is the second power offset, and n is the time node.

[0102] In an embodiment of the present application, the first power offset in step S210 is calculated by the following steps:

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

[0104] In the embodiment, the full power offset is multiplied by the relative power of the node to obtain the first power offset.

[0105] It should be noted that the relative power of the node at the time node is pre-set for the pressurized water reactor, and thus the first power offset obtained can be used as a reference for comparison with the second power offset.

[0106] It can be understood that the first power offset is calculated by the following formula: ref = AO0·P n ,

[0107] wherein, AO0 is the axial power offset of the pressurized water reactor at full power, P n is the relative power of the pressurized water reactor at the time node.

[0108] It can be understood that the second power offset can be obtained according to the rod position sequence by using the SCIENCE software or other core calculation software.

[0109] In an embodiment of the present application, as shown in FIG. 6, the fitness value of the rod position in step S220 is obtained by the following steps:

[0110] Step S221: If the axial power offset difference corresponding to the rod position is less than a pre-set first offset difference, the fitness value of the rod position is confirmed as a pre-set first fitness value.

[0111] Step S222: If the axial power offset difference corresponding to the rod position is greater than or equal to the pre-set first offset difference and less than or equal to a pre-set second offset difference, the fitness value of the rod position is confirmed as a pre-set second fitness value.

[0112] Step S223: If the axial power offset difference corresponding to the rod position is greater than the pre-set second offset difference, the fitness value of the rod position is confirmed as a pre-set third fitness value.

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

[0114] It can be understood that the fitness value of the rod position is calculated by the following formula:

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

[0116] It should be noted that when calculating the fitness value of a rod position, the axial power offset difference value of the rod position at the time node is calculated according to the first power offset and the second power offset at the time node corresponding to the rod position, and the fitness value of the rod position is further calculated according to the axial power offset difference value.

[0117] In an embodiment of the present application, as shown in FIG. 7, the fitness value of the rod position sequence in step S230 is calculated by the following steps:

[0118] Step S231: summing the fitness values of each rod position in the rod position sequence to obtain a rod position fitness overall value;

[0119] Step S232: selecting the largest axial power offset difference value from the axial power offset difference values corresponding to the plurality of rod positions in the rod position sequence;

[0120] Step S233: multiplying the largest axial power offset difference value by a preset correction coefficient to obtain a correction value, wherein the preset correction coefficient is a positive number;

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

[0122] In the embodiment, the rod position fitness overall value 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 can result in the fitness values of a plurality of rod position sequences being the same, which makes it difficult to guarantee the optimization effect of the iterative processing. Therefore, the fitness value of the rod position sequence is corrected by the correction value. Even if the fitness values of a plurality of rod position sequences are the same, the larger the largest value of the plurality of axial power offset difference values in the rod position sequence, the larger the correction value, so that the fitness value of the rod position sequence is smaller, thereby guaranteeing the optimization effect of the iterative processing.

[0123] It can be understood that the fitness value of the rod position sequence is calculated by the following formula:

[0124] wherein η is the fitness value of the rod position sequence, n is the time node, f(ΔI diff,n ) is the fitness value of the rod position, k is the correction coefficient and k is a positive number, and ΔI diff is the largest axial power offset difference value in the plurality of axial power offset difference values of the rod position sequence.

[0125] It can be understood that the correction coefficient can be any small positive number, for example, the correction coefficient is 0.0125.

[0126] The method for controlling axial power distribution of a pressurized water reactor provided in the embodiments of the present application can be executed by a control device 100 for controlling axial power distribution of a pressurized water reactor. The embodiments of the present application take the control device 100 for controlling axial power distribution of a pressurized water reactor as an example to illustrate the control device 100 for controlling axial power distribution of a pressurized water reactor provided in the embodiments of the present application.

[0127] Please refer to FIG. 8, which is a structural schematic diagram of the control device 100 for controlling axial power distribution of a pressurized water reactor provided in the embodiments of the present application. The control device 100 for controlling axial power distribution of a pressurized water reactor is applied to an electronic device 200. As shown in FIG. 8, the control device 100 includes:

[0128] A population generation module 110 is configured to generate an initial population including a plurality of individuals. Each individual is a rod position sequence of a temperature control rod. The rod position sequence includes a plurality of time nodes and a plurality of rod positions. The plurality of time nodes correspond to the plurality of rod positions.

[0129] A function determination module 120 is configured to determine a fitness evaluation function with the minimum axial power offset of the pressurized water reactor as an optimization target.

[0130] An iteration processing module 130 is configured to take the initial population as a target parent population and perform an iteration processing for finding an optimal solution based on the fitness evaluation function, so as to obtain a target rod position sequence associated with the initial population.

[0131] A rod position control module 140 is configured to control the rod 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.

[0132] In an embodiment of the present application, the iteration processing module 130 includes:

[0133] An individual fitness calculation module is configured to calculate the fitness value of each individual in the target parent population based on the fitness evaluation function.

[0134] A target parent population generation module is configured to perform selection operation, crossover operation and mutation operation on the target parent population based on the fitness value of each individual in the target parent population, so as to obtain a new target parent population.

[0135] In an embodiment of the present application, the population generation module 110 includes:

[0136] A feasible solution acquisition module is configured to acquire a feasible solution space. The feasible solution space includes a plurality of rod position sequences. Each rod position sequence satisfies the insertion limit of the temperature control rod and satisfies the extraction limit of the temperature control rod.

[0137] A first selection module is configured to select one rod position sequence from the feasible solution space as a transition sequence.

[0138] A step number applying module is configured to obtain an initial population by applying random step numbers to a plurality of rod positions in a transition sequence.

[0139] In an embodiment of the present application, the function determining module 120 comprises:

[0140] A first function establishing module is configured to establish an axial power offset function, which is configured to calculate an axial power offset difference value according to a first power offset and a second power offset, the first power offset being a reference offset of axial power of the pressurized water reactor when the temperature control rod is at the time node, and the second power offset being an actual offset of axial power of the pressurized water reactor when the temperature control rod is at the time node.

[0141] A second function establishing module is configured to establish a rod position fitness function based on the axial power offset function, the rod position fitness function being configured to calculate a fitness value of a rod position according to the axial power offset difference value.

[0142] A third function establishing module is configured to establish a fitness evaluation function based on the rod position fitness function, the fitness evaluation function being configured to calculate a fitness value of a rod position sequence according to fitness values of a plurality of rod positions of the rod position sequence.

[0143] In an embodiment of the present application, the first function establishing module comprises:

[0144] A transition value calculating module is configured to subtract the second power offset from the first power offset to obtain an axial power offset transition value.

[0145] An absolute value calculating module is configured to take an absolute value of the axial power offset transition value to obtain the axial power offset difference value.

[0146] In an embodiment of the present application, the transition value calculating module comprises:

[0147] A first power offset calculating module is configured to multiply a full power offset by a node relative power to obtain the first power offset, the full power offset being an axial power offset when the pressurized water reactor is at full power, and the node relative power being a relative power of the pressurized water reactor at the time node.

[0148] In an embodiment of the present application, the second function establishing module comprises:

[0149] A first judging module is configured to, if the axial power offset difference value is less than a preset first offset difference value, confirm that the fitness value of the rod position is a preset first fitness value.

[0150] A second judging module is configured to, if the axial power offset difference value is greater than or equal to the preset first offset difference value and less than or equal to a preset second offset difference value, confirm that the fitness value of the rod position is a preset second fitness value.

[0151] The third judging module is configured to confirm that the fitness value of the rod position is a preset third fitness value if the axial power offset difference value is greater than a preset second offset difference value.

[0152] In an embodiment of the present application, the third function establishing module comprises:

[0153] The sum module is configured to sum the fitness values of the rod positions in the rod position sequence to obtain a total fitness value of the rod positions.

[0154] The second selecting module is configured to select the maximum axial power offset difference value according to the axial power offset difference values corresponding to the plurality of rod positions in the rod position sequence.

[0155] The correction module is configured to multiply the maximum axial power offset difference value by a preset correction coefficient to obtain a correction value, wherein the preset correction coefficient is a positive number.

[0156] The difference obtaining module is configured to subtract the correction value from the total fitness value of the rod positions to obtain the fitness value of the rod position sequence.

[0157] The control device 100 provided by the embodiment of the present application can realize each process of the method embodiments of FIGS. 1 to 7 and achieve the same beneficial effects. To avoid repetition, details are not described herein.

[0158] In addition, an embodiment of the present application further discloses an electronic device 200, as shown in FIG. 9, comprising:

[0159] at least one processor 210;

[0160] at least one memory 220 configured to store at least one program;

[0161] The at least one program is executed by the at least one processor 210 to implement the control method of the axial power distribution of the pressurized water reactor as described above.

[0162] The electronic device 200 provided by the embodiment of the present application can realize each process of the method embodiments of FIGS. 1 to 7 and achieve the same beneficial effects. To avoid repetition, details are not described herein.

[0163] In addition, an embodiment of the present application further discloses a computer readable storage medium, wherein the computer readable storage medium stores a program executable by a processor, and the program executable by the processor is executed by the processor to implement the control method of the axial power distribution of the pressurized water reactor as described above.

[0164] The computer readable storage medium provided by the embodiment of the present application can realize each process of the method embodiments of FIGS. 1 to 7 and achieve the same beneficial effects. To avoid repetition, details are not described herein.

[0165] As will be appreciated by one of ordinary skill in the art, all or some steps, systems of the above-disclosed methods can be implemented as software, firmware, hardware, or 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, a digital signal processor, or a microprocessor, or as hardware, or as an integrated circuit, such as an application-specific integrated circuit. Such software can be distributed on computer readable media, which can comprise computer storage media (or non-transitory media), and communication media (or transitory media). As is well known to those of ordinary skill in the art, computer storage media includes all computer-readable media in which data, computer executable instructions, or other computer readable data is / are publicized, embodied, or otherwise accessed. Computer storage media does not include communication media unless the communication media facilitates access to computer readable data by the individual or machine. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can be accessed by the computer. Further, as will be appreciated by one skilled in the art, communication media typically embodies computer readable instructions, data structures, program modules or other data in a modulated data signal such as a carrier wave or other transport mechanism and includes any information delivery media.

[0166] The above detailed description of the application has been given to understand the application better, but the application is not limited to the above embodiments, and various changes can be made within the knowledge of those skilled in the art without departing from the spirit of the 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; 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.

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 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 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.

5. The method for controlling the axial power distribution of a pressurized water reactor according to claim 4, 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.

6. The method for controlling the axial power distribution of a pressurized water reactor according to claim 4, 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.

7. The method for controlling the axial power distribution of a pressurized water reactor according to claim 4, characterized in that, 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.

8. The method for controlling the axial power distribution of a pressurized water reactor according to claim 4, characterized in that, The fitness value of the bar sequence is calculated through the following steps: 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.

9. The method for controlling the axial power distribution of a pressurized water reactor according to claim 2, characterized in that, The selection, crossover, and mutation operations performed on the target parent population include: Multiple bar sequences are selected from the target parent population through a selection operation. These multiple bar sequences are then paired, and some bars in the paired bar sequences are crossed. Finally, the values ​​of the bars in the crossed bar sequences are changed.

10. 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.

Images