Method and device for generating reactor core loading scheme of nuclear reactor

Abstract

The invention discloses a nuclear reactor core loading scheme generation method and device, and the method comprises the steps: obtaining a first assembly position set of a reactor core of a nuclear reactor, determining a first assembly position corresponding to a shutdown rod assembly of the reactor core from the first assembly position set, and for the remaining assembly positions except the first assembly position, determining the first assembly position corresponding to the shutdown rod assembly of the reactor core; and determining at least one symmetric component position set corresponding to the remaining component positions, determining component types corresponding to the symmetric component position sets from the at least one preset component type, and generating a reactor core loading scheme of the reactor core according to the component types corresponding to the component positions in the first component position set. Therefore, the reactor core loading scheme can be generated without using a machine learning model or a manual arrangement mode, so that the condition that the reactor core loading scheme is determined by manpower or by acquiring a large number of training samples can be avoided, and the time cost of generating the reactor core loading scheme of the reactor core can be reduced.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of core design of nuclear reactors, and particularly relates to a method and device for generating a core loading scheme of a nuclear reactor. BACKGROUND

[0002] At present, for determination of a core loading scheme of a nuclear reactor, an artificial arrangement method or a machine learning method is generally used to achieve the determination. However, the artificial arrangement method may combine more inherent experience to determine the core loading scheme, so that the determination of the core loading scheme is limited by the inherent experience. In addition, the artificial arrangement method needs to be operated manually, which may result in high time cost of the core loading scheme.

[0003] Secondly, the machine learning method needs a large number of existing core loading schemes to train a machine learning model. Since a large number of training samples with diversity need to be obtained, the time cost of the machine learning method for determining the core loading scheme is increased, so that the time cost of the core loading scheme is also high.

[0004] It can be seen that the existing method for determining the core loading scheme has low efficiency. Therefore, there is an urgent need for a method to improve the determination efficiency of the core loading scheme. SUMMARY

[0005] Embodiments of the application provide a method and device for generating a core loading scheme of a nuclear reactor, which can improve the determination efficiency of the core loading scheme.

[0006] In a first aspect, a method for generating a core loading scheme of a nuclear reactor is provided, which includes: obtaining a first component position set of a core of the nuclear reactor; the first component position set includes a plurality of component positions; determining, from the first component position set, a first component position corresponding to a shutdown rod component of the core according to an argument of each component position in the first component position set under a preset coordinate system; the preset coordinate system is a rectangular coordinate system with a midpoint corresponding to the first component position set as an origin; for a remaining component position in the first component position set except the first component position, determining at least one symmetric component position set corresponding to the remaining component position; each symmetric component position set in the at least one symmetric component position set includes component positions having a symmetric relationship in the remaining component position; for the at least one symmetric component position set, determining a component type corresponding to the symmetric component position set from at least one preset component type; and generating a core loading scheme of the core according to the component type corresponding to each component position in the first component position set.

[0007] In a second aspect, a device for generating a core loading scheme of a nuclear reactor is provided, which includes:

[0008] An acquisition module, configured to acquire a first component position set of a core of a nuclear reactor; the first component position set comprising a plurality of component positions;

[0009] an assembly position determination module, configured to determine, from the first assembly position set, a first assembly position corresponding to a shutdown rod assembly of the reactor core based on the arguments of the assembly positions in the first assembly position set in a preset coordinate system; the preset coordinate system being a rectangular coordinate system with a midpoint corresponding to the first assembly position set as its origin;

[0010] a component position set determination module configured to determine, for the remaining component positions other than the first component position in the first component position set, at least one symmetric component position set corresponding to the remaining component positions; each symmetric component position set in the at least one symmetric component position set includes component positions having a symmetric relationship among the remaining component positions;

[0011] a component type determination module, configured to determine, for at least one symmetrical component position set, a component type corresponding to the symmetrical component position set from at least one preset component type;

[0012] The plan generating module is used to generate a core loading plan for the core according to the component type corresponding to each component position in the first component position set.

[0013] An embodiment of the present application provides a method for generating a core loading plan for a nuclear reactor. The method obtains a first component position set of the core of the nuclear reactor, determines a first component position corresponding to a shutdown rod assembly of the core from the first component position set based on the argument of each component position in the first component position set in a preset coordinate system, determines at least one symmetric component position set corresponding to each remaining component position for each component position other than the first component position in the first component position set, each symmetric component position set in the at least one symmetric component position set including component positions having a symmetric relationship among the remaining component positions, determines a component type corresponding to each symmetric component position set from at least one preset component type for each symmetric component position set, and generates a core loading plan for the core based on the component type corresponding to each component position in the first component position set. In this manner, the core loading plan can be generated without using a machine learning model or manual arrangement, thus avoiding the need to manually determine the core loading plan or obtain a large number of training samples. Therefore, the present application can reduce the time cost of generating the core loading plan for the core and effectively improve the efficiency of determining the core loading plan. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following is a brief introduction to the drawings required for use in the embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0015] Figure 1 1 is a flow chart of a method for generating a core loading plan for a nuclear reactor provided in an embodiment of the present application;

[0016] Figure 2 This is a schematic diagram of a complete embodiment provided in the embodiments of the present application;

[0017] Figure 3 is a schematic diagram of a component distribution array provided in an embodiment of the present application;

[0018] Figure 4 This is a schematic diagram of isolated points provided in an embodiment of the present application;

[0019] Figure 5 This is a schematic diagram of the core range provided by an embodiment of the present application;

[0020] Figure 6 This is a schematic diagram of component positions along a 270° line provided in an embodiment of the present application;

[0021] Figure 7 1 is a schematic diagram of the position of the first component of the shutdown rod assembly provided in an embodiment of the present application;

[0022] Figure 8 Schematic diagram of the symmetric relationship between A2 and A2sym provided in the embodiments of the present application;

[0023] Figure 9 is a schematic diagram of the symmetry relationship Hex provided in an embodiment of the present application;

[0024] Figure 10 Schematic diagram of the distribution of component types at the A2 component position provided in an embodiment of the present application;

[0025] Figure 11 Schematic diagram of the component type distribution of the A2sym component location provided in an embodiment of the present application;

[0026] Figure 12 is a schematic diagram of a core loading solution provided in an embodiment of the present application;

[0027] Figure 13 This is a schematic diagram of the control rod screening of the internal component positions provided by the embodiment of the present application;

[0028] Figure 14 is a schematic diagram of the control rods in the position of the external components provided in an embodiment of the present application;

[0029] Figure 15 is a schematic diagram of the position of the second component of the control rod provided in an embodiment of the present application;

[0030] Figure 16Schematic diagram of the control rod ring provided in an embodiment of the present application;

[0031] Figure 17 is a schematic diagram of control rod grouping provided in an embodiment of the present application;

[0032] Figure 18 1 is a schematic structural diagram of a core loading plan generating device for a nuclear reactor provided in an embodiment of the present application;

[0033] Figure 19 It is a structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0034] The features and exemplary embodiments of various aspects of the present application will be described in detail below. In order to make the purpose, technical solutions and advantages of the present application clearer, the present application will be further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only intended to explain the present application, rather than to limit the present application. For those skilled in the art, the present application can be implemented without the need for some of these specific details. The following description of the embodiments is merely to provide a better understanding of the present application by illustrating the examples of the present application.

[0035] It should be noted that, in this document, relational terms such as first and second, etc., are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "comprises," "comprising," or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, the elements defined by the phrase "comprising..." do not exclude the presence of other identical elements in the process, method, article, or device comprising the elements.

[0036] To address related technical issues, embodiments of the present application provide a method, apparatus, device, computer storage medium, and computer program product for generating a core loading plan for a nuclear reactor. The core loading plan generation method of embodiments of the present application is primarily used to generate a core loading plan for a reactor core with circular cassette components, a triangular arrangement, and externally positioned control rods.

[0037] The following first introduces the method for generating a core loading plan for a nuclear reactor provided in an embodiment of the present application.

[0038] Figure 1FIG. 1 is a flow chart showing a method for generating a core loading plan for a nuclear reactor provided by an embodiment of the present application. Figure 1 As shown, the core loading plan generating method specifically includes the following steps S101 to S105.

[0039] Step S101: obtaining a first component position set of a core of a nuclear reactor.

[0040] Step S102 : determining a first assembly position corresponding to a shutdown rod assembly of a reactor core from the first assembly position set according to the arguments of the assembly positions in a preset coordinate system.

[0041] Step S103 : for the remaining component positions except the first component position in the first component position set, determine at least one symmetric component position set corresponding to the remaining component positions.

[0042] Step S104: for at least one symmetrical component position set, determine a component type corresponding to the symmetrical component position set from at least one preset component type.

[0043] Step S105 : generating a core loading plan for the reactor core according to the component type corresponding to each component position in the first component position set.

[0044] The first component location set includes a plurality of component locations.

[0045] In the above step S101, in one embodiment, the position of each component in the above first component position set can be obtained by pre-setting. In another embodiment, the position of each component in the above first component position set can be calculated based on the core parameters of the nuclear reactor.

[0046] In one embodiment, each component position may have a corresponding serial number, and the first component position set may be obtained by integrating the serial numbers corresponding to multiple component positions, where each serial number is used to represent the corresponding component position.

[0047] The above-mentioned preset coordinate system is a rectangular coordinate system with the midpoint corresponding to the first component position set as the origin. It can be understood that the area formed by the various component positions in the first component position set can cover the origin of the preset coordinate system, and the midpoint of the area formed by the various component positions is the origin of the preset coordinate system.

[0048] In the above step S102, for each component position in the first component position set, the polar coordinate radius of the component position in the preset coordinate system can be calculated based on the coordinates of the component position in the preset coordinate system, and then the corresponding angular radius can be calculated based on the polar coordinate radius of the component position to obtain the angular radius of each component position in the first component position set.

[0049] In one embodiment, the calculation formula of the polar coordinate radius rr_FA can be shown as the following formula (1):

[0050]

[0051] Among them, rr_FA represents the polar coordinate radius, x_core represents the horizontal coordinate, and y_core represents the vertical coordinate.

[0052] In one embodiment, the calculation formula of the argument sita_FA can be shown as the following formula (2):

[0053]

[0054] Among them, sita_FA represents the argument, x_core represents the horizontal coordinate, and y_core represents the vertical coordinate.

[0055] There are many types of components in the core. Therefore, different component positions can be used to place different types of components. In this embodiment, it is necessary to determine the component type at each component position to generate a core loading plan. Component types can be divided into shutdown rods, control rods, fuel assemblies, etc.

[0056] In this embodiment, first, based on the arguments of the assembly positions, a first assembly position that is most suitable for placing the shutdown rod assembly of the reactor core is determined from the first assembly position set.

[0057] In the above-mentioned step S103, the component positions in the above-mentioned first component position set are symmetrical. In one embodiment, for the remaining component positions, a symmetrical component position set can be obtained by obtaining multiple symmetrical component positions in the remaining component positions. The above-mentioned symmetrical component position set includes the above-mentioned multiple symmetrical component positions, and the angular spacing between the two closest component positions in the above-mentioned multiple symmetrical component positions is a preset angle, and the polar coordinate radius of each component position in the multiple symmetrical component positions is the same.

[0058] In another embodiment, a symmetric relationship between a plurality of symmetric component positions of the remaining component positions may be established first, and a symmetric component position set may be determined based on the component positions having the symmetric relationship.

[0059] In this embodiment, the component positions in the first component position set have axial symmetry. In this embodiment, for the remaining component positions, multiple symmetrical component positions can be obtained by obtaining component positions symmetrical about a preset axis. Specifically, the preset axis can include multiple preset axes.

[0060] In one embodiment, the angles between the preset axes are the same, the axes of the preset axes pass through the origin of the preset coordinate system, and the angles between the preset axes can be preset angles.

[0061] For example, the preset axes may include a first preset axis, a second preset axis, a third preset axis, a fourth preset axis, a fifth preset axis and a sixth preset axis, and the preset angle between each preset axis is 60°. The axis with an angle of 0° is the first preset axis, the axis with an angle of 60° is the second preset axis, the axis with an angle of 120° is the third preset axis, the axis with an angle of 180° is the fourth preset axis, the axis with an angle of 240° is the fifth preset axis, and the axis with an angle of 300° is the sixth preset axis. If component position A and component position B are symmetrical about the first preset axis, component position B is symmetrical about the first preset axis, component position A is symmetrical about the first preset axis, component position B ... Position B and component position C are symmetrical about the second preset axis, component position C and component position D are symmetrical about the third preset axis, component position D and component position E are symmetrical about the fourth preset axis, component position E and component position F are symmetrical about the fifth preset axis, and component position F and component position A are symmetrical about the sixth preset axis. Then, the multiple symmetrical component positions include component position A, component position B, component position C, component position D, component position E and component position F, and a symmetrical relationship is established based on component position A, component position B, component position C, component position D, component position E and component position F.

[0062] It can also be understood that the component positions in the first component position set have rotational symmetry, that is, among the two closest component positions in the plurality of symmetrical component positions, one component position can be determined by rotating the other component position by a preset angle.

[0063] Each symmetric component position set in the at least one symmetric component position set includes component positions having a symmetric relationship among the remaining component positions.

[0064] In step S103, a symmetric component position set may be determined based on component positions having a symmetric relationship. For example, if a symmetric relationship is established between component position A, component position B, component position C, component position D, component position E, and component position F, the symmetric component position set includes component position A, component position B, component position C, component position D, component position E, and component position F.

[0065] In the above step S104, the component type corresponding to each component position in each symmetrical component position set may be randomly determined from at least one preset component type.

[0066] In one embodiment, for each component position in each symmetrical component position set, a target preset component type is randomly selected from at least one preset component type, and the target preset component type is determined as the component type corresponding to the component position.

[0067] In an embodiment, each preset component type can be numbered, and a preset component type set can be obtained according to the numbers of the preset component types. For each component position in each set of symmetric component positions, a target number can be randomly selected from the preset component type set, and a target preset component type corresponding to the target number can be determined as a component type corresponding to the component position.

[0068] In another embodiment, in order to ensure that the distribution of the component types is symmetrical, the component types corresponding to each set of symmetric component positions can be randomly determined from at least one preset component type, and the component types corresponding to the set of symmetric component positions can be determined as the component types corresponding to each component position in the set of symmetric component positions.

[0069] In the above step S105, the component types corresponding to each component position in the first set of component positions are taken as the content of the core loading scheme of the core, and the core loading scheme of the core is output.

[0070] In some embodiments, the above step S102 can include, but is not limited to, the following steps:

[0071] According to the azimuth angle of each component position in the first set of component positions, at least one reference component position of the shutdown rod component is determined from the first set of component positions.

[0072] The azimuth angle of each component position in the first set of component positions can be calculated according to the above formulas (1) and (2).

[0073] From the first set of component positions, a first symmetric position corresponding to each reference component position in the at least one reference component position is determined.

[0074] The above first symmetric position can be a component position symmetrical to the reference component position about a preset axis, that is, the reference component position and the first symmetric position are symmetrical about the preset axis, and the preset axis includes a plurality of preset axes.

[0075] For example, the preset axes may include a first preset axis, a second preset axis, a third preset axis, a fourth preset axis, a fifth preset axis, and a sixth preset axis, with the axis with an angle of 0° being the first preset axis, the axis with an angle of 60° being the second preset axis, the axis with an angle of 120° being the third preset axis, the axis with an angle of 180° being the fourth preset axis, the axis with an angle of 240° being the fifth preset axis, and the axis with an angle of 300° being the fifth preset axis. The reference component position is component position A. If component position A and component position B are symmetrical about the first preset axis, component position B and component position C are symmetrical about the second preset axis, component position C and component position D are symmetrical about the third preset axis, component position D and component position E are symmetrical about the fourth preset axis, component position E and component position F are symmetrical about the fifth preset axis, and component position F and component position A are symmetrical about the sixth preset axis, then component position B, component position C, component position D, component position E, and component position F are all first symmetrical positions corresponding to reference component position A.

[0076] In another embodiment, the first symmetrical position can be obtained by rotating the reference component position by a preset fixed angle, and the preset fixed angle can include multiple different angles, and the difference between two angles with the smallest difference among the multiple different angles is a preset angle.

[0077] For example, the preset fixed angles include 60°, 120°, 180°, 240°, and 300°. The five component positions determined by rotating the reference component position by 60°, 120°, 180°, 240°, and 300° respectively are determined as the first symmetrical positions.

[0078] In another embodiment, the dth first symmetric position can be obtained by rotating the reference component position by a preset angle. When d is not equal to D, d=d+1 is set, and the dth first symmetric position is used as the reference component position. The step of rotating the reference component position by the preset angle to obtain the dth first symmetric position is returned until d=D, and D=360 / preset angle.

[0079] At least one reference component position and a first symmetrical position are determined as a first component position.

[0080] In this embodiment, by first determining the reference component position, and then quickly determining the first symmetrical position based on the reference component position through symmetry, the efficiency of determining the first component position can be improved.

[0081] In some embodiments, the step of determining at least one reference assembly position of the shutdown rod assembly from the first assembly position set may include, but is not limited to, the following steps:

[0082] A second component position set is determined from the first component position set according to the arguments of the component positions in the first component position set.

[0083] The second component position set includes component positions with an angle equal to the first preset angle.

[0084] In one embodiment, the first preset angle can be preset according to demand. To facilitate subsequent calculations, the first preset angle can be an odd multiple of 30, for example, the first preset angle is 270°.

[0085] In one embodiment, component positions having an angle equal to a first preset angle are selected from the first component position set, and the second component position set is obtained based on the component positions having an angle equal to the first preset angle.

[0086] For example, the first preset angle is 270°, the first component position set includes component position A, component position B, component position C, component position D, component position E and component position F, and the angles of component position A, component position B, component position C, component position D, component position E and component position F are 0°, 30°, 60°, 90°, 270° and 270° respectively, then the second component position set includes component position E and component position F.

[0087] The component position closest to the circumference of the preset circle in the second component position set is determined as the first selected position.

[0088] The above-mentioned preset circle takes the origin of the preset coordinate system as its center.

[0089] In one embodiment, the preset circle can be set as required.

[0090] In one embodiment, the first selected position can be obtained by respectively obtaining the vertical distance between each component position in the second component position set and the circumference, and determining the component position with the smallest vertical distance in the second component position set as the component position closest to the circumference of the preset circle.

[0091] The first selected position and a second selected position corresponding to the first selected position are determined as reference component positions.

[0092] The second selected position is a component position adjacent to the first selected position in the first component position set in a direction close to the origin of the preset coordinate system.

[0093] In one embodiment, at least one adjacent component position adjacent to the first selected position can be obtained from the first component position set, the polar coordinate radius of each component position in the at least one adjacent component position can be obtained respectively, and the component position with the smallest polar coordinate radius in the at least one adjacent component position can be determined as the second selected position.

[0094] In this embodiment, the most appropriate reference component position can be determined by the distance between the component position and the circumference of the preset circle, thereby improving the effectiveness of determining the reference component position.

[0095] In some embodiments, before the step of determining the component position closest to the circumference of the preset circle in the second component position set as the first selected position, the following steps may be included but are not limited to:

[0096] The component position farthest from the origin of the preset coordinate system in the second component position set is determined as the third selected position.

[0097] In one embodiment, the third selected position may be obtained by obtaining the polar coordinate radius of each component position in the second component position set and determining the component position with the largest polar coordinate radius as the component position farthest from the origin of the preset coordinate system.

[0098] For example, the second component position set includes component position A, component position B, and component position C. The polar coordinate radii corresponding to component position A, component position B, and component position C are rr_FA1, rr_FA2, and rr_FA3, where rr_FA1>rr_FA2>rr_FA3, then component position A is determined as the third selected position.

[0099] Obtain a first polar coordinate radius corresponding to the third selected position in the preset coordinate system.

[0100] In one embodiment, according to the coordinates of the third selected position in the preset coordinate system, the polar coordinate radius corresponding to the third selected position is calculated based on the above formula (1) to obtain the first polar coordinate radius.

[0101] According to the principle of bisection of volume, the target radius of the preset circle is determined according to the first polar coordinate radius.

[0102] In one embodiment, the target radius may be calculated as shown in the following formula (3):

[0103]

[0104] Wherein, rSFA represents the target radius, and R1 represents the first polar coordinate radius.

[0105] A circle is drawn with the origin of the preset coordinate system as the center and the target radius as the radius to obtain the preset circle.

[0106] In this embodiment, a preset circle that is more suitable for determining the position of the reference component can be obtained by combining the first polar coordinate radius of the third selected position with the volume bisection principle.

[0107] In some embodiments, step S103 may include but is not limited to the following steps:

[0108] A third component position set is determined from the remaining component positions according to the argument of each component position in the preset coordinate system.

[0109] The third component position set includes component positions whose angles are greater than or equal to the second preset angle and less than or equal to the third preset angle.

[0110] In one embodiment, the third preset angle and the second preset angle may be multiples of 30, the second preset angle may be 0°, and the third preset angle may be 60°.

[0111] The third component position set is formed by obtaining the component positions whose angles in the remaining component positions are within the range formed by the second preset angle and the third preset angle. For example, the second preset angle is 0°, the third preset angle is 60°, and the range formed by the second preset angle and the third preset angle is [0°, 60°]. The remaining component positions include component position A, component position B, component position C, and component position D. The angles of component position A, component position B, component position C, and component position D are 15°, 20°, 45°, and 75°, respectively. Then the third component position set includes component position A, component position B, and component position C.

[0112] For each component position in the third component position set, at least one second symmetrical position having a symmetrical relationship with the component position is determined from the remaining component positions.

[0113] The second symmetrical position is symmetrical to the component position about a preset axis. The preset axis may include multiple preset axes. The two closest preset axes among the multiple preset axes are spaced by a preset angle. The preset angle may be equal to the angle between the second preset angle and the third preset angle.

[0114] For example, the third set of component positions includes component position A and component position B, the preset angle is 60°, the preset axes include an axis with a 0° argument as a first preset axis, an axis with a 60° argument as a second preset axis, an axis with a 120° argument as a third preset axis, an axis with a 180° argument as a fourth preset axis, an axis with a 240° argument as a fifth preset axis, and an axis with a 300° argument as a sixth preset axis, for the component position A, the component position A1 in the remaining component positions is symmetrical to the component position A about the first preset axis, the component position A2 in the remaining component positions is symmetrical to the component position A1 about the second preset axis, the component position A3 in the remaining component positions is symmetrical to the component position A2 about the third preset axis, the component position A4 in the remaining component positions is symmetrical to the component position A3 about the fourth preset axis, the component position A5 in the remaining component positions is symmetrical to the component position A4 about the fifth preset axis, and the component position A5 in the remaining component positions is symmetrical to the component position A about the sixth preset axis, the second symmetrical positions corresponding to the component position A include the component position A1, the component position A2, the component position A3, the component position A4, and the component position A5, for the component position B, the component position B1 in the remaining component positions is symmetrical to the component position B about the first preset axis, the component position B2 in the remaining component positions is symmetrical to the component position B1 about the second preset axis, the component position B3 in the remaining component positions is symmetrical to the component position B2 about the third preset axis, the component position B4 in the remaining component positions is symmetrical to the component position B3 about the fourth preset axis, the component position B5 in the remaining component positions is symmetrical to the component position B4 about the fifth preset axis, and the component position B5 in the remaining component positions is symmetrical to the component position B about the sixth preset axis, the second symmetrical positions corresponding to the component position B include the component position B1, the component position B2, the component position B3, the component position B4, and the component position B5.

[0115] In an implementation, for each component position in the third set of component positions, a component position corresponding to a rotation of the component position by a preset fixed angle is determined as a corresponding second symmetrical position, the fixed angle includes a plurality of different angles, and a difference between two angles with the smallest difference in the plurality of different angles is a preset angle.

[0116] For example, the preset angle is 60°, and the fixed angle can include 60°, 120°, 180°, 240°, and 300°.

[0117] Since for the component position with the second preset angle, the component position corresponding to the component position rotated by the preset angle is the component position with the third preset angle, but the component position with the third preset angle is the component position in the third component position set, therefore, for each component position in the third component position set, the component position corresponding to the component position rotated by the preset fixed angle is determined as the second symmetrical position corresponding to the component position. It can be that for each component position in the target third component position set, the component position corresponding to the component position rotated by the preset fixed angle is determined as the second symmetrical position corresponding to the component position. The target third component position set is the component position set after screening out the component position with the third preset angle from the third component position set.

[0118] For each component position in the third component position set, a symmetric component position set is obtained according to the component position and at least one second symmetric position corresponding to the component position.

[0119] The plurality of symmetrical component positions may include a component position and at least one second symmetrical position corresponding to the component position.

[0120] For example, the third component position set includes component position A and component position B, the second symmetrical positions corresponding to component position A include component position A1, component position A2, component position A3, component position A4 and component position A5, and the second symmetrical positions corresponding to component position B include component position B1, component position B2, component position B3, component position B4 and component position B5, then the first symmetrical component position set includes A, A1, A2, A3, A4 and A5, and the second symmetrical component position set includes B, B1, B2, B3, B4 and B5.

[0121] In one embodiment, a symmetric relationship among component positions in the remaining component positions may be established first based on at least one second symmetric position corresponding to each component position in the third component position set.

[0122] For example, the third component position set includes component position A and component position B, the second symmetrical positions corresponding to component position A include component position A1, component position A2, component position A3, component position A4 and component position A5, and the second symmetrical positions corresponding to component position B include component position B1, component position B2, component position B3, component position B4 and component position B5, then the symmetric relationship is A-A1-A2-A3-A4-A5, B-B1-B2-B3-B4-B5.

[0123] After the symmetric relationship is obtained, for each component position in the third component position set, at least one second symmetric position corresponding to the component position is determined from the symmetric relationship to obtain a symmetric component position set.

[0124] In some embodiments, the method of determining the third component position set from the remaining component positions based on the arguments of the respective component positions in the preset coordinate system may include, but is not limited to, the following steps:

[0125] A fourth component position set is determined from the remaining component positions according to the arguments of each component position in the preset coordinate system.

[0126] The above-mentioned fourth component position set includes component positions whose radial angles are greater than or equal to the second preset radial angle and less than or equal to the fourth preset radial angle, and the angle between the third preset radial angle and the second preset radial angle is greater than the angle between the fourth preset radial angle and the second preset radial angle, wherein the angle between the third preset radial angle and the second preset radial angle can be twice the angle between the fourth preset radial angle and the second preset radial angle, and the fourth preset radial angle can be a multiple of 30. For example, if the third preset radial angle is 60° and the second preset radial angle is 0°, then the fourth preset radial angle can be 30°.

[0127] The fourth component position set is formed by obtaining the component positions in the remaining component positions whose angles are within the range formed by the second preset angle and the fourth preset angle. For example, the second preset angle is 0°, the fourth preset angle is 30°, and the range formed by the second preset angle and the fourth preset angle is [0°, 30°]. The remaining component positions include component position A, component position B, and component position C. The angles of component position A, component position B, and component position C are 15°, 20°, and 45°, respectively. Then the fourth component position set includes component position A and component position B.

[0128] A fifth component position set is determined from the remaining component positions according to the respective component positions in the fourth component position set.

[0129] The component positions in the fifth component position set and the component positions in the fourth component position set are axially symmetrical about the target axis, and the target axis is an axis passing through the origin of the preset coordinate system and having an angle of the fourth preset angle.

[0130] The third component location set includes a fourth component location set and a fifth component location set.

[0131] For example, the fourth preset angle is 30°, and the fourth component position set includes component position A and component position B, wherein among the remaining component positions, component position C and component position A are axially symmetrical about the axis with an angle of 30°, and among the remaining component positions, component position D and component position B are axially symmetrical about the axis with an angle of 30°. Then the fifth component position set includes component position C and component position D, and the third component position set includes component position A, component position B, component position C and component position D.

[0132] In some embodiments, step S104 may include but is not limited to the following steps:

[0133] For each component position in the fourth component position set and the fifth component position set, a component type corresponding to the component position is determined from at least one preset component type.

[0134] For each component position in the fourth component position set and the fifth component position set, obtain the symmetric component position set corresponding to the component position, and determine the component type corresponding to the component position as the component type corresponding to each component position in the symmetric component position set corresponding to the component position.

[0135] In this embodiment, by first determining the component type of each component position in the fourth component position set and the fifth component position set, the component type corresponding to the symmetric component position set corresponding to the component position is quickly determined based on the symmetry relationship.

[0136] In some embodiments, for each component location in the fourth component location set and the fifth component location set, the step of determining, from at least one preset component type, a component type corresponding to the component location may include, but is not limited to, the following steps:

[0137] For each component position in the fourth component position set, a component type corresponding to the component position is determined from at least one preset component type.

[0138] For each component position in the fourth component position set, a target symmetric component position symmetric to the component position is obtained from the fifth component position set, and the component type corresponding to the component position is determined as the component type corresponding to the target symmetric component position.

[0139] For example, the fourth component position set includes component position A1 and component position B1, and the fifth component position set includes component position A2 and component position B2. There is a symmetrical relationship between component position A1 and component position A2, and there is a symmetrical relationship between component position B1 and component position B2. If the component types corresponding to component position A1 and component position B1 are component type A and component type B respectively, then the component type corresponding to component position A2 is determined to be component type A, and the component type corresponding to component position B2 is determined to be component type B.

[0140] In some embodiments, step S101 may include but is not limited to the following steps:

[0141] According to the core parameters of the nuclear reactor, the component distribution dimensions of the core of the nuclear reactor are determined.

[0142] The above-mentioned component distribution dimension refers to the number of rows and columns in which the components are arranged in the core.

[0143] In one embodiment, the component distribution dimensions of the core may be calculated based on the core parameters using a preset algorithm.

[0144] The above core parameters may include the core radius Rc and the component center distance b1. The preset algorithm may determine the row spacing of the components in the core according to the component center distance, and determine the component distribution dimension of the core according to the core radius and the row spacing.

[0145] The preset algorithm can be expressed as follows:

[0146]

[0147] Where xs represents the component distribution dimension, Rc represents the core radius, b1 represents the component center distance, and b1×sin(60°) represents the row spacing.

[0148] In the preset coordinate system, the component distribution array of the core is constructed according to the component distribution dimensions.

[0149] The above-mentioned component distribution array is composed of a plurality of component position arrangements.

[0150] The number of rows and columns of the above-mentioned component distribution array is the same as the component distribution dimension, and the center of the component distribution array is the origin of the preset coordinate system.

[0151] At least one fourth selected position is screened from the component distribution array based on the coordinates of each component position in the component distribution array in a preset coordinate system.

[0152] The fourth selected position is within a preset core radius.

[0153] The first component position set is composed using the at least one fourth selected position.

[0154] In some embodiments, the core parameters include the component center distance b1, the component box outer radius rFA, and the core radius Rc. The step of selecting at least one fourth selected position from the component distribution array may also include, but is not limited to, the following steps:

[0155] For each component position, the coordinates of the component position are determined based on the component center distance, the component distribution dimension, and the target coordinates of the component position.

[0156] The above-mentioned target coordinates represent the row and column number of the component position in the component distribution array.

[0157] In a possible implementation, the number of rows of the component distribution array may be arranged from bottom to top, and the number of columns of the component distribution array may be arranged from left to right.

[0158] The coordinate calculation formulas for the component position are shown in the following formulas (5) and (6):

[0159]

[0160] Among them, (i, j) represents the target coordinates, i represents the number of rows, j represents the number of columns, x_core represents the horizontal coordinate, y_core represents the vertical coordinate, b1 represents the component center distance, and xs represents the component distribution dimension.

[0161] For each component position, the target distance corresponding to the component position is determined according to the coordinates of the component position and the outer radius of the component box.

[0162] The target distance is the distance between the component box boundary of the component position and the origin of the preset coordinate system.

[0163] In one embodiment, the target distance is calculated as shown in the following formula (7):

[0164]

[0165] Among them, rr_core represents the target distance, x_core and y_core represent the coordinates, and rFA represents the component box boundary.

[0166] For each assembly position, if the target distance of the assembly position is less than the core radius, the assembly position is determined as the first assembly position.

[0167] In some embodiments, after determining the first component location set, isolated component locations may be screened from the first component location set, and the first component location set after removing the isolated component locations may be determined as the final first component location set. Therefore, after determining the first component location set, the following steps may also be included but are not limited to:

[0168] An isolated component position reference set is obtained from the first component position set.

[0169] The first sequence number corresponding to the component position in the above-mentioned isolated component position reference set is equal to the sum of the target first sequence number and the target difference value, the target difference value is the second sequence number corresponding to the component position minus one, and the above-mentioned target first sequence number is the minimum first sequence number corresponding to the first component position set.

[0170] The first serial number is the serial number of the component position in the component distribution array of the core, and the second serial number is the serial number of the component position in the first component position set.

[0171] The formula for obtaining the isolated component position reference set n1 is shown in the following formula (8):

[0172] n1={Uno(i)|Uno(i)=Uno(1)+i-1,i∈[1,card(Uno)]} (8);

[0173] Among them, Uno(1) represents the target first serial number, that is, the first serial number corresponding to the minimum second serial number, card(Uno) represents the number of component positions in the first component position set, i represents the second serial number, Uno(i) represents the first serial number corresponding to the component position with the second serial number i, and n1 includes the component position corresponding to Uno(i).

[0174] For each component position in the isolated component position reference set, an isolated component position corresponding to the component position is determined from the first component position set.

[0175] There is a symmetrical relationship between the above-mentioned isolated component positions and component positions.

[0176] The isolated component position and the component position may be symmetrical about a preset axis, and the preset axis may include multiple preset axes.

[0177] Each component position in the isolated component position reference set and the isolated component position corresponding to each component position are deleted from the first component position set to obtain a filtered first component position set.

[0178] In some embodiments, step S105 may include but is not limited to the following steps:

[0179] Get the first sequence number and the second sequence number corresponding to each component position in the first component position set.

[0180] The first serial number corresponding to each component position is the serial number of the component position in the component distribution array of the core, and the second serial number corresponding to each component position is the serial number of the component position in the first component position set.

[0181] Before step S105, the component positions in the component distribution array may be arranged sequentially from left to right starting from the minimum row number to obtain the first serial number of each component position in the component distribution array. In one embodiment, the first serial number of the component position may also be calculated based on the target coordinates of the component position and the component distribution dimension. The calculation formula may be specifically as shown in the following formula (9):

[0182] nn=xs×(i-1)+j,i,j∈[1,xs] (9);

[0183] Among them, nn represents the first serial number of the component position, i and j represent the target coordinates of the component position, and xs represents the component distribution dimension.

[0184] The component positions in the first component position set may be re-sorted from small to large according to the first serial numbers to obtain the second serial numbers of the component positions in the first component position set.

[0185] An external component location reference set is determined from the first component location set according to the first sequence number and the second sequence number.

[0186] The first serial number corresponding to the component position in the external component position reference set is equal to the sum of the minimum first serial number and the target difference value, and the target difference value is the second serial number corresponding to the component position minus one.

[0187] The formula for determining the component position in the external component position reference set ex1 is shown in the following formula (10):

[0188] ex1={Uno(z)|Uno(z)=Uno(1)+z-1,z∈[1,card(Uno)]} (10);

[0189] Among them, Uno(1) represents the minimum first serial number corresponding to the first component position set, card(Uno) represents the number of component positions in the first component position set, z represents the second serial number, Uno(z) represents the first serial number corresponding to the component position with the second serial number z, and ex1 includes the component position corresponding to Uno(z).

[0190] The symmetric component position set corresponding to each component position in the external component position reference set is determined as the external component position set.

[0191] An internal component position set is determined based on component positions in the first component position set except for the external component position set.

[0192] According to the coordinates of each component position in the target component position set in a preset coordinate system, the second component position of the control rod around each component position is determined.

[0193] The target component position set is an internal component position set or an external component position set, and the number of second component positions around the component position of the internal component position set is greater than the number of second component positions around the component position of the external component position set.

[0194] A core loading plan for the core is generated according to the component type corresponding to each component position in the first component position set and the second component position.

[0195] In some embodiments, the step of determining the second assembly position of the control rods around each assembly position may include, but is not limited to, the following steps:

[0196] If the target component position set is an internal component position set, then for each component position in the internal component position set, the fifth selected position of the control rod around the component position is determined based on the coordinates of the component position in the preset coordinate system, the third serial number of the component position in the internal component position set, and the preset serial number.

[0197] The above-mentioned preset sequence number ranges from 1 to a preset number of control rods M, where the preset number of control rods M is a preset maximum number of control rods around the component position.

[0198] The calculation method of the fifth selected position can be shown in the following formulas (11) and (12):

[0199] xcr(M×(i-1)+j)=x_core(IFA(i))+b1 / 2 / sin(60°)×cos(30°+θ×(j-1)),

[0200] i∈[1,card(IFA)],j∈[1,M] (11);

[0201] ycr(M×(i-1)+j)=y_core(IFA(i))+b1 / 2 / sin(60°)×sin(30°+θ×(j-1)),

[0202] i∈[1,card(IFA)],j∈[1,M] (12);

[0203] Where M represents the preset number of control rods, card(IFA) represents the number of component positions in the internal component position set, i represents the third index of the component position in the internal component position set, j represents the preset index, and θ represents the quotient of 360° and M. For example, if M is 6, then θ is 60°. xcr(M×(i-1)+j) represents the abscissa of the j-th control rod in the component position corresponding to the third index i, and ycr(M×(i-1)+j) represents the ordinate of the j-th control rod in the component position corresponding to the third index i.

[0204] According to the above formulas (11) and (12), M fifth selected positions corresponding to each component position in the internal component position set can be obtained.

[0205] In one embodiment, the maximum number of control rods may be 6, and the preset serial numbers are from 1 to 6. The component positions in the internal component position set may be reordered from smallest to largest according to the first serial numbers to obtain third serial numbers for the component positions in the internal component position set.

[0206] For each component position in the internal component position set, the polar coordinate radius of the fifth selected position corresponding to the component position is obtained. If the polar coordinate radius of the fifth selected position is smaller than the second polar coordinate radius, the fifth selected position is determined as the second component position.

[0207] The second polar coordinate radius is the polar coordinate radius of the component position corresponding to the largest first serial number in the external component position reference set.

[0208] In some embodiments, the step of determining the position of the second component may include, but is not limited to, the following steps:

[0209] If the target component position set is an external component position set, then for each component position in the external component position set, the sixth selected position of the control rod around the component position is determined based on the coordinates of the component position in the preset coordinate system, the fourth serial number of the component position in the external component position set, and the preset serial number.

[0210] The sixth selected position can be calculated according to the above formulas (11) and (12).

[0211] For each component position in the first target position set in the external component position set, the polar coordinate radius of the sixth selected position corresponding to the component position is obtained, and the first preset number of seventh selected positions in the sixth selected position corresponding to the component position are determined as the second component position.

[0212] The first target position set includes a symmetric component position set corresponding to the component positions corresponding to the smallest first serial number and the largest first serial number in the external component position reference set.

[0213] The polar coordinate radius of the seventh selected position is smaller than the polar coordinate radius of the component positions in the sixth selected position that are not the seventh selected position.

[0214] In one embodiment, the sixth selected positions corresponding to the component positions are arranged from smallest to largest polar coordinate radius, and the sixth selected position located in the first preset number of positions is determined as the seventh selected position. Alternatively, the sixth selected positions corresponding to the component positions are arranged from largest to smallest polar coordinate radius, and the sixth selected position located in the last preset number of positions is determined as the seventh selected position. The first preset number can be 4.

[0215] For each component position in the second target position set in the external component position set, the polar coordinate radius of the sixth selected position corresponding to the component position is obtained, and the second preset number of eighth selected positions in the sixth selected positions corresponding to the component position are determined as the second component position.

[0216] The above-mentioned second target position set includes component positions in the external component position set except the first target position set, the polar coordinate radius of the eighth selected position is smaller than the polar coordinate radius of the component position in the sixth selected position that is not the eighth selected position, and the second preset number is greater than the first preset number.

[0217] In one embodiment, the sixth selected positions corresponding to the component positions are arranged from smallest to largest polar coordinate radius, and the sixth selected position located at the second-preset number before is determined as the eighth selected position. Alternatively, the sixth selected positions corresponding to the component positions are arranged from largest to smallest polar coordinate radius, and the sixth selected position located at the second-preset number after is determined as the eighth selected position. The second preset number can be 5.

[0218] In another embodiment, the step of determining the second assembly position of the control rods around each assembly position may include, but is not limited to, the following steps:

[0219] If the target component position set is an external component position set, then for each component position in the external component position reference set, the sixth selected position of the control rod around the component position is determined based on the coordinates of the component position in the preset coordinate system, the fourth serial number of the component position in the external component position set, and the preset serial number.

[0220] For the first target position in the external component position reference set, the polar coordinate radius of the sixth selected position corresponding to the first target position is obtained, and a first preset number of seventh selected positions in the sixth selected position corresponding to the first target position are obtained.

[0221] The first target position includes the component position corresponding to the smallest first serial number and the largest first serial number in the external component position reference set.

[0222] For the second target position in the external component position set, the polar coordinate radius of the sixth selected position corresponding to the component position is obtained, and a second preset number of eighth selected positions in the sixth selected position corresponding to the component position are obtained.

[0223] The second target position set includes component positions in the external component position set except the first target positions.

[0224] After the target selected position is rotated for the dth time by the target angle, the corresponding rotation component position is obtained, the target selected position includes the seventh selected position and the eighth selected position, and the target angle is the preset angle×d.

[0225] If d is not equal to D, the process returns to the step of rotating the target selected position by the target angle for the dth time to obtain the corresponding rotation component position until d is equal to D, where D=360 / preset angle.

[0226] The target selected position and the rotated component position are determined as a second component position.

[0227] All the determined second component positions are deduplicated to obtain a final plurality of second component positions.

[0228] In some embodiments, the control rod assemblies need to be grouped. Therefore, the step of generating a core loading plan for the reactor core based on the assembly type corresponding to each assembly position in the first assembly position set and the second assembly position may include, but is not limited to, the following steps:

[0229] Get the polar coordinate radius corresponding to each second component position.

[0230] The second component positions are divided into G sixth component position sets corresponding to different ranges according to the maximum polar coordinate radius among the polar coordinate radii corresponding to the second component positions and a preset ratio.

[0231] G is a positive integer greater than 1.

[0232] In one embodiment, the second component position representation area is divided into G circles according to the maximum polar coordinate radius and the preset ratio, wherein the radius of each circle can be calculated as follows:

[0233]

[0234] Among them, r_cir(u) represents the radius of the u-th cycle, and r_cir(G_cycle) represents the maximum polar coordinate radius.

[0235] The preset ratio may be 1:3:3:...:3, wherein the number of 3s in the preset ratio is G-1.

[0236] Taking the center circle as a range, and the range formed by each adjacent circle of other circles except the center circle as a range, G different ranges can be obtained.

[0237] Obtain G component positions within different ranges respectively to obtain G sixth component position sets.

[0238] The target sixth component position set among the G sixth component position sets is divided into a preset number of seventh component position sets.

[0239] The above-mentioned target sixth component position set does not include the component positions corresponding to the origin of the preset coordinate system, that is, the target sixth component position set is G sixth component position sets, excluding the G-1 sixth component position sets corresponding to the center circle.

[0240] The preset number of groups can be 3.

[0241] The seventh component position set in the wth target sixth component position set among the G sixth component position sets is separated from the adjacent seventh component position set in the w+1th target sixth component position set by a preset angle, where w takes a value from 1 to G-2, and the area corresponding to the wth target sixth component position set is adjacent to the area corresponding to the w+1th target sixth component position set.

[0242] For example, the preset angle is 60, and the second component position representation area is divided into 3 circles, wherein the 3 circles are sorted in sequence from the origin outward, and the 3 circles correspond to the first sixth component position set, the second sixth component position set, and the third sixth component position set, respectively. Among them, the first sixth component position set includes the component position corresponding to the origin of the preset coordinate system, then the second sixth component position set and the third sixth component position set are the first target sixth component position set and the second target sixth component position set, respectively, and the target sixth component position set is divided into 3 seventh component position sets, and the angle ranges corresponding to the 3 seventh component position sets of the first target sixth component position set are [0, 120], [120, 240], and [240, 360], respectively. Then the angle ranges corresponding to the 3 seventh component position sets of the second target sixth component position set are [60, 180], [180, 300], and [300, 60], respectively.

[0243] A core loading plan for the core is generated based on the component type corresponding to each component position in the first component position set, the second component position, the seventh component position set of the target sixth component position set, and the sixth component position sets other than the target sixth component position set in the G sixth component position sets.

[0244] In this embodiment, the second component position is divided into G sixth component position sets corresponding to different ranges through the maximum polar coordinate radius and the preset ratio, and the target sixth component position set in the G sixth component position sets is divided into a preset number of seventh component position sets, so that the control rods can be automatically grouped.

[0245] In order to better understand the above method, an embodiment of the present application provides a complete embodiment of a method for generating a core loading plan for a nuclear reactor.

[0246] This complete embodiment includes six steps, namely, core range determination, shutdown rod assembly determination, core symmetry relationship determination, core randomization scheme generation, control rod position determination, and control rod grouping. The core range determination step is used for a first component position set of a core of a nuclear reactor, the shutdown rod assembly determination step is used to determine a first component position corresponding to the shutdown rod assembly, the core symmetry relationship determination step is used to determine the symmetry relationship of the component positions, the core randomization scheme generation step is used to determine the component types of component positions other than the second component position in the first component position set, the control rod position determination step is used to determine the second component position corresponding to the control rod, and the control rod grouping step is used to group the determined second component positions to obtain a seventh component position set of G-1 target sixth component position sets and a sixth component position set other than the target sixth component position set in the G sixth component position sets.

[0247] Reference Figure 2 , Figure 2 A schematic diagram of the complete process of the embodiment provided in this application.

[0248] First, the steps for determining the core range are explained. Figure 2 The first point in the core range determines the process content.

[0249] The input core parameters are obtained. The core parameters include the core radius Rc, the outer radius of the circular assembly box rFA, the outer radius of the control rod guide tube rCR, the assembly center distance b1, the total number of radial control rod cycles G_cycle, and the total number of fuel assembly types (including shutdown rod assemblies) n_type.

[0250] Calculate the dimension of the component distribution array covering the center (0, 0) of the preset coordinate system and the core radius Rc (i.e., the component distribution dimension). The calculation formula is as follows:

[0251]

[0252] Wherein, b1×sin(60°) is the row spacing of the component distribution array in the y direction.

[0253] Calculate the array coordinates x_core and y_core covering the core range. The array coordinates are all xs×xs matrices. The row and column numbers (i, j) of the component distribution array (i.e., the target coordinates) are uniquely corresponding to the position number (first number) nn of the array coordinates in the component distribution array, as shown in the following example: Figure 3 As shown, Figure 3 A schematic diagram of a component distribution array shown in an embodiment of the present application.

[0254] The calculation formulas for array coordinates x_core and y_core are as follows:

[0255]

[0256] nn=xs×(i-1)+j,i,j∈[1,xs];

[0257] Calculate the target distance rr_core of the component box boundary relative to the center of the circle (i.e., the origin of the preset coordinate system) at each component position in the component distribution array. rr_core is an xs×xs matrix. The calculation formula for rr_core is as follows:

[0258]

[0259] Obtain the component serial number Uno within the core range according to rr_core, and obtain the first component position set Uno. Uno is a positive integer set and is calculated as follows:

[0260] Uno={nn|rr_core(nn) <Rc,nn∈[1,xs×xs]};

[0261] That is, the first serial number corresponding to the component position in the component distribution array whose target distance is less than the core radius is obtained, and Uno is composed according to the obtained first serial number.

[0262] Based on the continuity of Uno data and coordinate symmetry, the outermost isolated points are removed (for example, if the number of components in the first row of Uno is 1 or 2, then that row and its symmetrical positions are called isolated points) to determine the final Uno. Based on the continuity of Uno data, the one or two isolated points n1 to be removed are determined, where n1 is a set of one or two positive integers, which are the first one or two elements in Uno.

[0263] Reorder the component positions corresponding to the first serial numbers in Uno from small to large according to the first serial numbers to obtain the second serial numbers of the component positions in Uno.

[0264] The formula for obtaining the isolated component position reference set n1 is as follows:

[0265] n1={Uno(i)|Uno(i)=Uno(1)+i-1,i∈[1,card(Uno)]};

[0266] Among them, Uno(1) represents the target first sequence number, that is, the first sequence number corresponding to the minimum second sequence number, card(Uno) represents the number of first sequence numbers in Uno, i represents the second sequence number, Uno(i) represents the first sequence number corresponding to the component position with the second sequence number i, and n1 includes Uno(i).

[0267] For the isolated component position reference set n1, the first serial number that is symmetrical to each first serial number in the isolated component position reference set is obtained from Uno. In this embodiment, the coordinates of the component position corresponding to each first serial number (isolated point) in n1 are rotated 60 degrees five times to obtain the remaining five symmetrical isolated point sets n in Uno. k (k∈{2, 3, 4, 5, 6}), n k Including the first serial number determined by rotating the coordinates of the component position corresponding to each first serial number (isolated point) in n1 by 60°×(k-1), that is, obtaining the first serial number symmetrical with the isolated component position reference set n1 about the preset axis, the preset axis includes the first preset axis, the second preset axis, the third preset axis, the fourth preset axis, the fifth preset axis and the sixth preset axis, the first preset axis, the second preset axis, the third preset axis, the fourth preset axis, the fifth preset axis and the sixth preset axis are sequentially spaced by preset angles, and the preset angle is 60°.

[0268] The updated Uno is 6 sets of isolated point sets n kThe complement in the original Uno.

[0269] The formula of coordinate rotation of the isolated point n1 is as follows:

[0270]

[0271] The isolated point set n k The acquisition formula is as follows:

[0272]

[0273] Uno (new) = C Uno u_del.

[0274] Then, the stop rod assembly determination step is entered, and the flow content of the second point stop rod assembly determination in the Figure 2 can be referred to.

[0275] The polar coordinate radius rr_FA and the azimuth sita_FA (angle system) of the component position corresponding to each first serial number in Uno are calculated.

[0276] The polar coordinate radius rr_FA is calculated as follows:

[0277]

[0278] The azimuth sita_FA is calculated as follows:

[0279]

[0280] The first preset azimuth is 270°, and the first serial number set A1 (i.e., the second component position set) on the 270° line is obtained according to sita_FA, A1 is a positive integer set, and the calculation method is as follows:

[0281] A1 = {nn | sita_FA (nn) = 270°, nn ∈ Uno};

[0282] According to the first polar coordinate radius R1 corresponding to the outermost component position (the third selected position) in A1, the target radius rSFA is calculated according to the volume bisection principle. Since the target coordinate (i, j) is uniquely corresponding to the first serial number nn, the outermost component position in A1 is the component position corresponding to the smallest first serial number in the A1 set.

[0283] The target radius calculation formula is as follows:

[0284] R1 = rr_FA (min nn∈A1 nn);

[0285]

[0286] The component position corresponding to the first serial number in A1 that is closest to the circumference of a preset circle with the target radius rSFA as the radius and the origin as the center (the first selected position) is taken as the first component position of the first group of shutdown rods. The component positions corresponding to the two adjacent first serial numbers in a row of component positions toward the center of the circle are taken as the second and third component positions of the shutdown rods (the second selected positions).

[0287] The shutdown rod assembly sequence number is denoted as SFA, which is a 6×3 matrix. Each row represents the same group of shutdown rods. The calculation method is:

[0288] SFA(1,1)={nn|rr_FA(nn)=min ii∈A1 |rr_FA(ii)-rSFA|,nn∈A1};

[0289] SFA(1,2)=SFA(1)+xs-1;

[0290] SFA(1,3)=SFA(1)+xs;

[0291] That is, three reference component positions are obtained, including one first selected position SFA(1,1) and two second selected positions SFA(1,2) and SFA(1,3).

[0292] Based on coordinate symmetry, the serial numbers of the remaining five groups of shutdown rod assemblies are determined. The coordinates of the first group of shutdown rod positions (three reference assembly positions) are rotated 60° five times to obtain the corresponding first symmetric positions. This allows the serial numbers of the remaining five symmetric groups of shutdown rods to be obtained. The calculation formula is as follows:

[0293]

[0294] That is, the component position determined by rotating the reference component position by 60°×(i-1) is determined as the first symmetric position. According to the first symmetric position determined by rotating by 60°×(i-1), the i-th group of shutdown rod positions (i.e., the i-th group of first component positions) is formed, i∈[2,6].

[0295] Then, enter the step of determining the core symmetry relationship, which can be referred to Figure 2 The process content is determined by the core symmetry relationship in the third point.

[0296] The second preset angle of argument is 0°, and the fourth preset angle of argument is 30°. According to the angle sita_FA of the component position corresponding to each first serial number in Uno, the first serial number set A2 (fourth component position set) corresponding to the angle of argument in the range of [0°, 30°) is obtained from the remaining first serial numbers except the first serial number of the first component position in Uno, where A2 is a set of positive integers.

[0297] The formula obtained in A2 is as follows:

[0298] A2={nn|0≤sita_FA(nn)≤30°,nn∈C Uno SFA};

[0299] According to the axially symmetrical relationship between the coordinates, the first serial number set A2sym (fifth component position set) of the symmetrical component positions of A2 is determined. A2sym is a set of positive integers, each element of which uniquely corresponds to an element in A2. The component position corresponding to the first serial number in A2sym is symmetrical with the component position corresponding to the first serial number in A2 about the axis with a fourth preset radial angle of 30° (i.e., the target axis).

[0300] The formula for obtaining A2sym is as follows:

[0301]

[0302] Five 60° coordinate rotations are performed on the component positions within the range [0°, 60°) (i.e., the fourth component position set and the component positions in the fifth component position set excluding the component positions with an angle equal to the sum of the preset angle and the second preset angle) to obtain the corresponding second symmetric positions. This determines the symmetric relationship Hex between the component positions within the core. Hex is a 6-column matrix, the first column of which is the first sequence number set A3 (excluding the shutdown rod assembly) within [0°, 60°) (i.e., the fourth component position set and the component positions in the fifth component position set excluding the component positions with an angle equal to the sum of the preset angle and the second preset angle). In addition, because A3 includes the component position at the center (0, 0), five repeated calculations are performed during the symmetry operation. Therefore, the number of elements in the A3 set card(A3) satisfies card(Uno) = card(A3)*6-5+18. The number of rows in Hex is (card(Uno)-18+5) / 6. The remaining five columns in Hex are the first sequence numbers corresponding to the second symmetric positions after the rotation.

[0303] A3 obtains the formula as follows:

[0304] A3={nn|0≤sita_FA(nn)<60°,nn∈C Uno SFA};

[0305] The formula for determining Hex is as follows:

[0306] Hex(i,1)=A3(i),i∈[1,(card(Uno)-18+5) / 6];

[0307]

[0308] Then, enter the core random solution generation step. Figure 2The fourth point is the process content of core random scheme generation.

[0309] Set n_type-1 preset component types, number each preset component type, and the number can be a positive integer from 1 to n_type-1. The shutdown rod component type is numbered n_type.

[0310] A set of random numbers with the same length as A2 is generated according to the number of the preset component type, and the set of random numbers corresponds to the component type of the corresponding component position of A2.

[0311] The core loading plan is denoted as LP, which is an xs×xs matrix that uniquely corresponds to the x_core and y_core positions. The elements in LP are integers between 0 and n_type, representing the component type at that position. LP 0 means it is outside the core Rc range. Initializing LP to Oxs and denoting the random number as temp, we have:

[0312] LP(SFA(i,j))=n_type,i∈[1,6],j∈[1,3];

[0313] LP(A2(i))=temp(i),i∈[1,card(A2)];

[0314] Here, card(A2) indicates the number of first serial numbers in A2.

[0315] That is, for each first serial number in the fourth component position set A2, the component type corresponding to each first serial number is determined from the preset component types.

[0316] A2sym corresponds one-to-one with the first serial number in A2, thereby determining the component type at the component position corresponding to A2sym. At this point, the loading plan for the six shutdown rod assemblies and conventional components within the range of [0°, 60°) has been determined.

[0317] The formula for determining the component type corresponding to each first sequence number in A2sym is as follows:

[0318] LP(A2sym(t))=LP(A2(t)),t∈[1,card(A2)];

[0319] That is, obtain the first serial number corresponding to the target symmetric component position symmetrical to the component position with the first serial number in A2 from the fifth component position set A2sym, and determine the component type corresponding to the component position with the first serial number in A2 as the component type of the corresponding target symmetric component position.

[0320] According to the symmetry Hex of the whole stack components, the whole stack loading plan is determined. The specific formula is as follows:

[0321] LP(Hex(i,j))=LP(Hex(1,j)),i∈[1,(card(Uno)-18+5) / 6],j∈[2,6];

[0322] That is, for each first serial number in the fourth component position set and the fifth component position set, according to the symmetry relationship Hex, the symmetric component position set corresponding to the first serial number is obtained, and the component type corresponding to the first serial number is determined as the component type corresponding to each symmetric first serial number in the corresponding symmetric component position set, so as to obtain the component type of each component position.

[0323] Then, enter the control rod position determination step, you can refer to Figure 2 The fifth point is the process content for determining the control rod position.

[0324] According to Uno continuity and Hex symmetry, component positions are divided into two categories: internal and external component positions, which are denoted as IFA (internal component position set) and EFA (external component position set), respectively, and are both sets of positive integers.

[0325] First, obtain the external component position reference set nn_ex5, and the calculation formula is as follows:

[0326] nn_ex5={Uno(i)|Uno(i)=Uno(1)+i-1,i∈[1,card(Uno)]};

[0327] According to the symmetry relationship Hex, the position set of symmetric components that have a symmetric relationship with nn_ex5 is determined. The determination formula is as follows:

[0328] ii5={s|Hex(s,5)=nn,nn∈nn_ex5};

[0329] nn_ex k ={Hex(s,k)|s∈ii5},k∈[1, 2, 3, 4, 6];

[0330] The symmetrical component position set corresponding to each component position in the external component position reference set is determined as the external component position set EFA, and the determination formula is as follows:

[0331]

[0332] According to the component positions in the first component position set Uno except the external component position set, the internal component position set is determined. The determination formula is as follows:

[0333] IFA=C Uno EFA;

[0334] Set the preset number of control rods M to 6, retain 6 control rods outside the internal component position, and calculate the coordinate matrices xcr and ycr of the control rod positions (fifth selected positions) around the internal component position. The formula is as follows:

[0335] xcr(6×(i-1)+j)=x_core(IFA(i))+b1 / 2 / sin(60°)×cos(30°+60°×(j-1)),i∈

[0336] [1,card(IFA)],j∈[1,6];

[0337] ycr(6×(i-1)+j)=y_core(IFA(i))+b1 / 2 / sin(60°)×sin(30°+60°×(j-1)),i∈

[0338] [1,card(IFA)],j∈[1,6];

[0339] The preset serial numbers in the above formula are 1 to 6, and i in the above formula is the third serial number.

[0340] And calculate the polar coordinate radius rc of the fifth selected position of the internal component position, the formula is as follows:

[0341]

[0342] The polar coordinate radius of the outer edge component position is taken as r0 (the second polar coordinate radius), that is, the polar coordinate radius of the component position corresponding to the maximum value in nn_ex5.

[0343] According to rc, the fifth selected position of the control rod within the range of r0 is retained to obtain the position of the second component, and the xcr and ycr of the second component position are updated:

[0344] r0=rr_FA(max nn∈nn_ex1 nn);

[0345] xcr={xcr(ii)|rc(ii)≤r0,ii∈[1,card(IFA)×6]};

[0346] ycr={ycr(ii)|rc(ii)≤r0,ii∈[1,card(IFA)×6]};

[0347] For the outer assembly position, the first preset number is 4, the second preset number is 5, the edge assembly position (the first target position) of the outer assembly position reserves 4 control rods, and the middle assembly position (the second target position set) reserves 5 control rods. The coordinates of the seventh selected position of the control rods around the edge assembly position and the coordinates of the eighth selected position of the control rods around the middle assembly position are calculated xcrbis, ycrbis:

[0348]

[0349] In the above formula, j is a preset serial number, i in the above formula is a fourth serial number, and xcrbis k , ycrbis k is the coordinates corresponding to the rotating assembly position, xcrbis1(6×(i-1)+j) represents the horizontal coordinate of the jth control rod of the assembly position corresponding to the fourth serial number i, and ycrbis1(6×(i-1)+j) represents the vertical coordinate of the jth control rod of the assembly position corresponding to the fourth serial number i.

[0350] The seventh selected position of the control rods around the edge assembly position and the eighth selected position of the control rods around the middle assembly position are the second assembly positions.

[0351] The second coordinate positions of the control rods corresponding to the inner and outer assembly positions are merged and de-duplicated to obtain the final second assembly position coordinates xx and yy of the control rods.

[0352] Then, the control rod grouping step is entered, and the flow content of the seventh point control rod grouping in the above can be referred to. Figure 2

[0353] The polar coordinate radius r_CR and the argument sita_CR of the second assembly position of the control rods are calculated, and the calculation formula is as follows:

[0354]

[0355] The maximum value of r_CR is taken as the maximum radius (maximum polar coordinate radius) of control rod grouping, and is divided into G_cycle circles according to the volume 1:3:3:... principle. The radius r_cir corresponding to each G_cycle circle is calculated, and r_cir contains G_cycle radius data. The calculation formula is as follows:

[0356] r_cir(G_cycle)=max ii∈[1,card(xx)] rr_CR(ii);

[0357]

[0358] ​r_cir(i) to r_cir(i+1) is a range, and according to the component positions in each range, G sixth component position sets are obtained. The sixth component position set includes the first serial number corresponding to the component position of the corresponding circle.

[0359] The preset angle is 60°, the preset number of groups is 3, and the range outside the center circle corresponds to the sixth component position set (target sixth component position set). According to sita_CR, the same range is divided into 3 groups, and the two adjacent groups are staggered by 60° to obtain the corresponding seventh component position set. The xx and yy coordinates of the control rods included in each group of the seventh component position set are numbered Gi, where Gi is a positive integer set. The control rods are divided into 3×(G_cycle-1)+1 groups in total. The calculation formula is as follows:

[0360] G1={ii|rr_CR(ii) <r_cir(1),ii∈[1,card(xx)]};

[0361]

[0362] G 1+3×(g-2)+j =C {ii|r_cir(g-1)≤rr_cR(ii)＜r_cir(g)} (G 1+3×(g-2)+1 ∪G 1+3×(g-2)+2 ),g∈[2,G_cycle],j=3.

[0363] A core loading plan of the core is generated according to the component type corresponding to each component position in the first component position set, the second component position, the seventh component position set, and the sixth component position set in the sixth component position set except the target component position set.

[0364] Taking the input parameters in Table 1 as an example, the method of the present invention is used to determine the core range, generate a plan, determine and group control rods.

[0365] Table 1

[0366]

[0367] First, the core range is determined as follows:

[0368] Calculate the component distribution dimension xs=29 covering the circle center (0,0) and radius Rc.

[0369] Calculate the array coordinates x_core, y_core, x_core, y_core, row and column numbers (i, j) (target coordinates), and the first sequence number nn covering the core range.

[0370] Calculated based on the target distance rr_core, the initial Uno is as follows:

[0371] One={79 104 105 106 107 108 109 110 111 131 132 133 134 135 136 137 138 139 140 141 158 159 160 161 162 163 164 165 166 167 168 169 170 171 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 497 498 499 500 501 502 503 504505 506 507 508 509 510 511 512 513 514 515516 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543544 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 671 672 673 674 675 676 677 678 679 680 681 682 683 684 701 702 703 704 705 706 707 708 709 710 711 731732 733 734 735 736 737 738 763};

[0372] Among them, {79} in Uno is the isolated point n1.

[0373] According to the continuity, we can get the isolated points of each group n1 = {79}, n2 = {259}, n3 = {601}, n4 = {763}, n5 = {583}, n6 = {241}, u_del = {79 259 601 763 583 241}, as shown in the following example: Figure 4 As shown, Figure 4 The isolated point diagram provided in the embodiment of the present application is shown in FIG. 1 , in which the red-marked positions are isolated points. Figure 4 As shown, rr_core is Figure 4 The distance of the shown component's bounds relative to the circle's center.

[0374] After the update, Uno is as follows:

[0375] One={104 105 106 107 108 109 110 111 131 132 133 134 135 136 137 138 139 140 141 158 159 160 161 162 163 164 165 166 167 168 169 170 171 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 497 498 499 500 501 502 503 504 505 506 507508 509 510 511 512 513 514 515 516 526 527528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 555 556557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 584 585 586587 588 589 590 591 592 593 594 595 596 597 598 599 600 613 614 615 616 617618 619 620 621 622 623 624 625 626 627 628 642 643 644 645 646 647 648 649650 651 652 653 654 655 656 671 672 673 674 675 676 677 678 679 680 681 682683 684 701 702 703 704 705 706 707 708 709 710 711 731 732 733 734 735 736737 738}.

[0376] Finally, the core range is as Figure 5 As shown, card (Uno) = 379, Figure 5 A schematic diagram of the core range provided in an embodiment of the present application.

[0377] Then, determine the shutdown rod assembly as follows:

[0378] rr_FA and sita_FA are the polar coordinate radius and argument of the component position in the core, respectively.

[0379] The first serial number of the component position on the 270° line is A1 = {136 193 250 307 364}, refer to Figure 6 , Figure 6 Schematic diagram of component positions on the 270° line provided in an embodiment of the present application.

[0380] The first serial number of the outermost component position of A1 (the third selected position) is 136, the polar coordinate radius R1 = rr_FA(136) = 69.28 cm, and the target radius rSFA = 48.99 cm.

[0381] The first serial number of the component position closest to rSFA in A1 is 193, SFA(1, 1) = 193, SFA(1, 2) = 221, and SFA(1, 3) = 222.

[0382] After coordinate rotation, the serial numbers of the 6 groups of shutdown rod assemblies are:

[0383] SFA = [193 221 222 313 312 341 541 512 540 649 621 620 529 530 501 301 330 302].

[0384] Referring to Figure 7 , Figure 7 A schematic diagram of the first component position of the shutdown rod assembly provided in the embodiment of the present application.

[0385] Then, the core symmetry relationship is determined as follows:

[0386] The first serial number A2 of the component position in the range of [0°, 30°] (not including the shutdown rod) is A2 = {421 422 423 424 425 426 427 428 429 430 431 451 452 453 454 455 456 457 458 459 481 482 483 484 485 486 487 488 511 513 514 515 516 542 543 544 571 572}.

[0387] The first serial number of the symmetric component position of the component position corresponding to A2 about the 30° line is A2sym = {421 450 479 508 537 566 595 624 653 682 711 451 480 509 538 567 596 625 654 683 481 510 539 568 597 626 655 684 511 569 598 627 656 570 599 628 571 600}, and the component and corresponding relationship of A2 and A2sym are as shown in Figure 8 Figure 8 A schematic diagram of the symmetry relationship of A2 and A2sym provided in the embodiment of the present application.

[0388] ​42 423 424 425 426 427 428 429 430 431 451 452 453 454 455 456 457 458 459 480 481 482 483 484 485 486 487 488 509 510 511 513 514 515 516 538 539 542 543 544 567 568 569 570 571 572 596 597 598 599 600 625 626 627 628 654 655 656 683 684}, card(A3)=61, the symmetry relationship Hex is as follows, the symmetry relationship of the core component position is as follows Figure 9 As shown, Figure 9 A schematic diagram of the symmetry relationship Hex provided in an embodiment of the present application.

[0389] Hex=[421 421 421 421 421 421 422 450 449 420 392 393 423 479 477 419363 365 424 508 505 418 334 337 425 537 533 417 305 309 426 566 561 416 276281 427 595 589 415 247 253 428 624 617 414 218 225 429 653 645 413 189 197430 682 673 412 160 169 431 711 701 411 131 141 451 478 448 391 364 394 452507 476 390 335 366 453 536 504 389 306 338 454 565 532 388 277 310 455 594560 387 248 282 456 623 588 386 219 254 457 652 616 385 190 226 458 681 644384 161 198 459 710 672 383 132 170 480 506 447 362 336 395 481 535 475 361307 367 482 564 503 360 278 339 483 593 531 359 249 311 484 622 559 358 220283 485 651 587 357 191 255 486 680 615 356 162 227 487 709 643 355 133 199488 738 671 354 104 171 509 534 446 333 308 396 510 563 474 332 279 368 511592 502 331 250 340 513 650 558 329 192 284 514 679 586 328 163 256 515 708614 327 134 228 516 737 642 326 105 200 538 562 445 304 280 397 539 591 473303 251 369 542 678 557 300 164 285 543 707 585 299 135 257 544 736 613 298106 229 567 590 444 275 252 398568 619 472 274 223 370 569 648 500 273 194342 570 677 528 272 165 314 571 706 556 271 136 286 572 735 584 270 107 258596 618 443 246 224 399 597 647 471 245 195 371 598 676 499 244 166 343 599705 527 243 137 315 ​​600 734 555 242 108 287 625 646 442 217 196 400 626 675470 216 167 372 627 704 498 215 138 344 628 733 526 214 109 316 654 674 441188 168 401 655 703 469 187 139 373 656 732 497 186 110 345 683 702 440 159140 402 684 731 468 158 111 374].

[0390] Then, the core randomization scheme is generated as follows:

[0391] Generate a random number temp corresponding to the length of A2 = [5 5 1 5 4 1 2 3 5 5 1 5 5 3 5 1 3 5 4 54 1 5 5 4 4 4 2 4 1 4 1 2 1 1 5 4 2]. After the loading plan LP is initialized and the first assembly position and A2 position of the shutdown rod are confirmed, the component type distribution is as follows: Figure 10 As shown, Figure 10 A schematic diagram of the distribution of component types at the A2 component position provided in an embodiment of the present application.

[0392] According to the symmetry, the component type distribution of the A2sym component position is determined as follows Figure 11 As shown, Figure 11 A schematic diagram of the component type distribution of the A2sym component location provided in an embodiment of the present application.

[0393] According to coordinate symmetry, the full core loading scheme is:

[0394] LP=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 5 2 2 5 22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 5 4 1 1 4 1 1 4 5 1 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 2 5 5 4 4 4 1 1 4 4 4 5 5 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 4 45 5 4 1 6 1 4 5 5 4 4 2 0 0 0 0 0 0 0 0 0 0 0 0 0 5 1 4 5 3 3 5 6 6 5 3 3 5 41 5 0 0 0 0 0 0 0 0 0 0 0 0 2 1 4 4 3 2 1 5 4 5 1 2 3 4 4 1 2 0 0 0 0 0 0 0 00 0 0 2 4 1 1 5 1 1 5 1 1 5 1 1 5 1 1 4 2 0 0 0 0 0 0 0 0 0 0 5 1 1 6 6 5 5 43 4 3 4 5 5 6 6 1 1 5 0 0 0 0 0 0 0 0 0 2 1 4 1 6 4 1 3 5 5 5 5 3 1 4 6 1 4 12 0 0 0 0 0 0 0 0 2 4 4 4 5 5 1 4 5 1 5 1 5 4 1 5 5 4 4 4 2 0 0 0 0 0 0 0 0 54 5 3 1 5 3 5 5 5 5 5 5 3 5 1 3 5 4 5 0 0 0 0 0 0 0 0 1 5 5 3 2 1 4 5 1 5 5 51 5 4 1 2 3 5 5 1 0 0 0 0 0 0 0 0 5 4 5 3 1 5 3 5 5 5 5 5 5 3 5 1 3 5 4 5 0 00 0 0 0 0 0 2 4 4 4 5 5 1 4 5 1 5 1 5 4 1 5 5 4 4 4 2 0 0 0 0 0 0 0 0 2 1 4 16 4 1 3 55 5 5 3 1 4 6 1 4 1 2 0 0 0 0 0 0 0 0 0 5 1 1 6 6 5 5 4 3 4 3 4 5 56 6 1 1 5 0 0 0 0 0 0 0 0 0 0 2 4 1 1 5 1 1 5 1 1 5 1 1 5 1 1 4 2 0 0 0 0 0 0 0 0 0 0 0 2 1 4 4 3 2 1 5 4 5 1 2 3 4 4 1 2 0 0 0 0 0 0 0 0 0 0 0 0 5 1 4 5 33 5 6 6 5 3 3 5 4 1 5 0 0 0 0 0 0 0 0 0 0 0 0 0 2 4 4 5 5 4 1 6 1 4 5 5 4 4 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 5 5 4 4 4 1 1 4 4 4 5 5 2 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 1 5 4 1 1 4 1 1 4 5 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 2 5 22 5 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...

[0395] Loading plan such as Figure 12 As shown, Figure 12 A schematic diagram of the core loading scheme provided in an embodiment of the present application.

[0396] The control rod positions are then determined as follows:

[0397] According to Uno continuity, nn_ex5 = {104 105 106 107 108 109 110 111}, according to core symmetry, nn_ex6 = {171 200 229 258 287 316 345 374}, nn_ex1 = {488 516 544 572 600 628 656 684}, nn_ex2 = {738 737 736 735 734 733 732 731}, nn_ex3 = {671 642 613 584 555 526 497 468}, nn_ex4 = {354 326 298 270 242 214 186 158}, the first serial number of the peripheral component position is EFA = {104 105 106 107 108 109 110 111 171 200 229 258 287 316 345 374 488 516 544 572 600 628 656 684 738 737 736 735 734 733 732 731 671 642 613 584 555 526 497 468 354 326 298 270 242 214 186 158}, and the first serial number of the internal component position is IFA = {131 132 133 134 135 136 137 138 139 140 141 159 160 161 162 163 164 165 166 167 168 169 170 187 188 189 190 191 192 193 194 195 196 197 198 199 215 216 217 218 219 220 221 222 223 224 225 226 227 228 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 355 356 357 358 359 360 361 362 363 364 365 366 367368369 370 371 372 373 383 384 385 386 387 388 389 390 391 392 393 394 395 396397 398 399 400 401 402 411 412 413 414 415 416 417 418 419 420 421 422 423424 425 426 427 428 429 430 431 440 441 442 443 444 445 446 447 448 449 450451 452 453 454 455 456 457 458 459 469 470 471 472 473 474 475 476 477 478479 480 481 482 483 484 485 486 487 498 499 500 501 502 503 504 505 506 507508 509 510 511 512 513 514 515 527 528 529 530 531 532 533 534 535 536 537538 539 540 541 542 543 556 557 558 559 560 561 562 563 564 565 566 567 568569 570 571 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 614615 616 617 618 619 620 621 622 623 624 625 626 627 643 644 645 646 647 648649 650 651 652 653 654 655 672 673 674 675 676 677 678 679 680 681 682 683701 702 703 704 705 706 707 708 709 710 711}.

[0398] Calculate the coordinates xcr, ycr and polar radius rc of the six control rods around the internal component position IFA.

[0399] Take the polar coordinate radius corresponding to the outer edge component position (i.e. the component position corresponding to the maximum value 111 in nn_ex5), r0 = 81.19 cm, and the control rod screening of the inner component position is as follows: Figure 13 As shown, Figure 13A schematic diagram of the control rod screening of the internal component positions provided in an embodiment of the present application.

[0400] In the external component positions, 4 control rods are reserved at the edge component positions and 5 control rods are reserved at the middle component positions. According to the control rod coordinates corresponding to nn_ex5, the coordinates of the remaining 5 groups of corresponding position control rods are obtained after rotation, such as Figure 14 As shown, Figure 14 Schematic diagram of the control rods in the position of the external components provided in an embodiment of the present application.

[0401] The control rods corresponding to the positions of the inner and outer components are merged and deduplicated to obtain the final coordinates xx and yy of the control rods, such as Figure 15 As shown, Figure 15 A schematic diagram of the position of the second component of the control rod provided in an embodiment of the present application.

[0402] Finally, the control rods are grouped as follows:

[0403] Calculate the polar coordinate radius r_CR and argument sita_CR of the second assembly position of the control rod.

[0404] The maximum value of r_CR is taken as the maximum radius of the control rod group, r_cir(4) = 82.11 cm. It is divided into 4 circles according to the volume 1:3:3:3 principle. The corresponding radius is calculated as r_cir = [25.97 51.93 68.70 82.11]. The range from r_cir(i) to r_cir(i+1) is as follows: Figure 16 As shown, Figure 16 A schematic diagram of the control rod circles provided in an embodiment of the present application.

[0405] Except for the center circle, the same range is divided into 3 groups according to sita_CR, and the adjacent groups are staggered by 60°, such as Figure 17 As shown, Figure 17 This is a schematic diagram of the control rod grouping provided in an embodiment of the present application, wherein the number of control rods in each group is 84, 72, 72, 72, 80, 80, 80, 72, 72, 72, and 72, respectively, and the number of control rods in each group is approximately the same.

[0406] Based on the above method, the core loading plan can be generated without using a machine learning model or manual arrangement method. This can avoid determining the core loading plan manually or by obtaining a large number of training samples. Therefore, the present application can reduce the time cost of generating the core loading plan of the core and effectively improve the efficiency of determining the core loading plan.

[0407] In addition, the component type corresponding to the component position can be determined randomly, which can avoid excessive dependence on manual experience to determine the component type and avoid being limited by manual experience, thereby improving the diversity of core loading schemes.

[0408] Finally, the control rods are grouped efficiently and quickly through calculation.

[0409] The embodiments of the present application can quickly determine the core range and generate random solutions based on differentiated user needs, thereby increasing the diversity of core loading solutions, expanding the search range of core loading solutions, and providing more possibilities for optimizing core loading.

[0410] In order to better implement the above method, the embodiment of the present application provides a core loading plan generating device for a nuclear reactor, referring to Figure 18 , Figure 18 This is a structural block diagram of a core loading plan generating device for a nuclear reactor provided in an embodiment of the present application. The core loading plan generating device 20 for a nuclear reactor includes:

[0411] The acquisition module 201 is configured to acquire a first component position set of a core of a nuclear reactor; the first component position set includes a plurality of component positions.

[0412] The assembly position determination module 202 is configured to determine a first assembly position corresponding to a shutdown rod assembly of the core from the first assembly position set based on the arguments of the assembly positions in the first assembly position set in a preset coordinate system; the preset coordinate system is a rectangular coordinate system with the midpoint corresponding to the first assembly position set as its origin.

[0413] The component position set determination module 203 is used to determine at least one symmetric component position set corresponding to the remaining component positions other than the first component position in the first component position set; each symmetric component position set in at least one symmetric component position set includes component positions with symmetric relationships among the remaining component positions.

[0414] The component type determination module 204 is configured to determine, for at least one symmetrical component position set, a component type corresponding to the symmetrical component position set from at least one preset component type.

[0415] The plan generating module 205 is configured to generate a core loading plan for the reactor core according to the component type corresponding to each component position in the first component position set.

[0416] In some embodiments, the assembly position determination module 202 is specifically configured to: determine at least one reference assembly position of the shutdown rod assembly from the first assembly position set based on the argument of each assembly position in the first assembly position set; determine a first symmetric position corresponding to each reference assembly position in the at least one reference assembly position from the first assembly position set; and determine the at least one reference assembly position and the first symmetric position as the first assembly position.

[0417] In some embodiments, the component position determination module 202 is specifically used to: determine a second component position set from the first component position set based on the angular distribution of each component position in the first component position set; the second component position set includes component positions with an angular distribution of a first preset angular distribution; the component position in the second component position set that is closest to the circumference of a preset circle is determined as the first selected position; the preset circle has the origin of a preset coordinate system as the center; the first selected position and the second selected position corresponding to the first selected position are determined as reference component positions; the second selected position is the component position in the first component position set that is adjacent to the first selected position in the direction close to the origin of the preset coordinate system.

[0418] In some embodiments, the component position determination module 202 is specifically used to: concentrate the second component positions, and determine the component position farthest from the origin of the preset coordinate system as the third selected position; obtain the first polar coordinate radius corresponding to the third selected position in the preset coordinate system; determine the target radius of the preset circle according to the first polar coordinate radius according to the principle of bisection of the volume; draw a circle with the origin of the preset coordinate system as the center and the target radius as the radius to obtain the preset circle.

[0419] In some embodiments, the component position set determination module 203 is specifically used to: determine a fourth component position set from the remaining component positions based on the angular distribution of each component position in the preset coordinate system, the fourth component position set including component positions whose corresponding angular distribution is greater than or equal to the second preset angular distribution and less than or equal to the fourth preset angular distribution; determine a fifth component position set from the remaining component positions based on each component position in the fourth component position set; the component positions in the fifth component position set and the component positions in the fourth component position set have an axially symmetrical relationship about the target axis, and the target axis is an axis passing through the origin of the preset coordinate system and having an angle of the fourth preset angular distribution; for each component position in the fourth component position set and the fifth component position set, determine at least one second symmetrical position that has a symmetrical relationship with the component position from the remaining component positions; for each component position in the fourth component position set and the fifth component position set, obtain a symmetrical component position set based on the component position and the at least one second symmetrical position corresponding to the component position.

[0420] In some embodiments, the acquisition module 201 is specifically used to: determine the component distribution dimension of the core of the nuclear reactor based on the core parameters of the nuclear reactor, the component distribution dimension being the number of rows and columns of the component arrangement distribution in the core; construct a component distribution array of the core based on the component distribution dimension in a preset coordinate system; the component distribution array is composed of a plurality of component position arrangements, the number of rows and columns of the component distribution array is the same as the component distribution dimension, and the center of the component distribution array is the origin of the preset coordinate system; based on the coordinates of each component position in the component distribution array in the preset coordinate system, screen at least one fourth selected position from the component distribution array; the fourth selected position is within the preset core radius; based on at least one fourth selected position, a first component position set is formed.

[0421] The core parameters include component center distance, component box outer radius and core radius. In some embodiments, the acquisition module 201 is specifically used to: for each component position, determine the coordinates of the component position based on the component center distance, component distribution dimension and target coordinates of the component position; the target coordinates represent the row and column number of the component position in the component distribution array; for each component position, determine the target distance corresponding to the component position based on the coordinates of the component position and the component box outer radius; the target distance is the distance between the component box boundary of the component position and the origin of the preset coordinate system; for each component position, if the target distance of the component position is less than the core radius, the component position is determined as the first component position.

[0422] In some embodiments, the above-mentioned scheme generation module 205 is specifically used to: obtain the first serial number and the second serial number corresponding to each component position in the first component position set; the first serial number corresponding to each component position is the serial number of the component position in the component distribution array of the core, and the second serial number corresponding to each component position is the serial number of the component position in the first component position set; according to the first serial number and the second serial number, an external component position reference set is determined from the first component position set; the first serial number corresponding to the component position in the external component position reference set is equal to the sum of the minimum first serial number and the target difference, and the target difference is the second serial number corresponding to the component position minus one; the components of the external component position reference set are summed. The corresponding symmetrical component position set is determined as the external component position set; the internal component position set is determined based on the component positions in the first component position set except the external component position set; the second component position of the control rod around each component position is determined based on the coordinates of each component position in the target component position set in a preset coordinate system; the target component position set is the internal component position set or the external component position set, and the number of second component positions around the component position of the internal component position set is greater than the number of second component positions around the component position of the external component position set; according to the component type corresponding to each component position in the first component position set and the second component position, a core loading plan of the core is generated.

[0423] In some embodiments, the above-mentioned scheme generation module 205 is specifically used to: if the target component position set is an internal component position set, then for each component position in the internal component position set, according to the coordinates of the component position in the preset coordinate system, the third serial number of the component position in the internal component position set and the preset serial number, determine the fifth selected position of the control rod around the component position; the preset serial number ranges from 1 to the preset number of control rods, and the preset number of control rods is the preset maximum number of control rods around the component position; for each component position in the internal component position set, obtain the polar coordinate radius of the fifth selected position corresponding to the component position, if the polar coordinate radius of the fifth selected position is smaller than the second polar coordinate radius, determine the fifth selected position as the second component position; the second polar coordinate radius is the polar coordinate radius of the component position corresponding to the largest first serial number in the external component position reference set.

[0424] In some embodiments, the solution generation module 205 is specifically used to: if the target component position set is an external component position set, then for each component position in the external component position reference set, according to the coordinates of the component position in the preset coordinate system, the fourth serial number of the component position in the external component position set, and the preset serial number, determine the sixth selected position of the control rod around the component position; for each component position in the first target position set in the external component position set, obtain the polar coordinate radius of the sixth selected position corresponding to the component position, and determine the first preset number of seventh selected positions in the sixth selected position corresponding to the component position as the second component position; the first target position set includes the smallest first serial number and the largest seventh selected position in the external component position reference set. For the symmetrical component position set corresponding to the component position corresponding to the first serial number, the polar coordinate radius of the seventh selected position is smaller than the polar coordinate radius of the component position in the sixth selected position that is not the seventh selected position; for each component position in the second target position set in the external component position set, the polar coordinate radius of the sixth selected position corresponding to the component position is obtained, and the second preset number of eighth selected positions in the sixth selected positions corresponding to the component position are determined as the second component position; the second target position set includes the component positions in the external component position set except the first target position set, the polar coordinate radius of the eighth selected position is smaller than the polar coordinate radius of the component position in the sixth selected position that is not the eighth selected position, and the second preset number is greater than the first preset number.

[0425] In some embodiments, the above-mentioned scheme generation module 205 is specifically used to: obtain the polar coordinate radius corresponding to each second component position; divide the second component position into G sixth component position sets corresponding to different ranges according to the maximum polar coordinate radius and a preset ratio among the polar coordinate radii corresponding to each second component position, where G is a positive integer greater than 1; divide the target sixth component position set in the G sixth component position sets into a preset number of seventh component position sets; the target sixth component position set is the sixth component position set that does not include the component position corresponding to the origin of the preset coordinate system, and the seventh component position set in the wth target sixth component position set in the G sixth component position sets is separated from the adjacent seventh component position set in the w+1th target sixth component position set by a preset angle, and w takes a value from 1 to G-2; generate a core loading scheme for the core according to the component type corresponding to each component position in the first component position set, the second component position, the seventh component position set of the target component position set, and the sixth component position sets other than the target sixth component position set in the G sixth component position sets.

[0426] Based on the above-mentioned device, the core loading plan can be generated without using a machine learning model or manual arrangement method. This can avoid determining the core loading plan manually or by obtaining a large number of training samples. Therefore, the present application can reduce the time cost of generating the core loading plan of the core and effectively improve the efficiency of determining the core loading plan.

[0427] Figure 19 A schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present application is shown.

[0428] The electronic device may include a processor 301 and a memory 302 storing computer program instructions.

[0429] Specifically, the processor 301 may include a central processing unit (CPU), or an application-specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application.

[0430] The memory 302 may include a large capacity memory for data or instructions. By way of example and not limitation, the memory 302 may include a hard disk drive (HDD), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a magnetic tape, or a universal serial bus (USB) drive, or a combination of two or more of these. Where appropriate, the memory 302 may include removable or non-removable (or fixed) media. Where appropriate, the memory 302 may be inside or outside the integrated gateway disaster recovery device. In a specific embodiment, the memory 302 is a non-volatile solid-state memory.

[0431] In some embodiments, the memory 302 may include read-only memory (ROM), random access memory (RAM), magnetic disk storage media devices, optical storage media devices, flash memory devices, electrical, optical, or other physical / tangible memory storage devices. Thus, generally, the memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., memory devices) encoded with software including computer-executable instructions, and when the software is executed (e.g., by one or more processors), it is operable to perform the operations described with reference to the method according to an aspect of the present disclosure.

[0432] The processor 301 reads and executes computer program instructions stored in the memory 302 to implement any one of the methods for generating a core loading plan for a nuclear reactor in the above embodiments.

[0433] In one example, the electronic device may further include a communication interface 303 and a bus 310. Figure 19 As shown, the processor 301 , the memory 302 , and the communication interface 303 are connected via a bus 310 and communicate with each other.

[0434] The communication interface 303 is mainly used to implement communication between various modules, devices, units and / or equipment in the embodiments of the present application.

[0435] Bus 310 comprises hardware, software or both, couples the parts of electronic equipment to each other.For example, and not limitation, bus can comprise accelerated graphics port (AGP) or other graphics bus, enhanced industry standard architecture (EISA) bus, front side bus (FSB), hypertransport (HT) interconnection, industry standard architecture (ISA) bus, infinite bandwidth interconnection, low pin count (LPC) bus, memory bus, micro channel architecture (MCA) bus, peripheral component interconnection (PCI) bus, PCI-Express (PCI-X) bus, serial advanced technology attachment (SATA) bus, video electronics standard association local (VLB) bus or other suitable bus or two or more of these combinations.In suitable cases, bus 310 can comprise one or more buses.Although the present application embodiment describes and shows specific bus, the application considers any suitable bus or interconnection.

[0436] In addition, an embodiment of the present application provides a computer storage medium, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the method for generating a core loading plan for the nuclear reactor is implemented.

[0437] An embodiment of the present application provides a computer program product. When instructions in the computer program product are executed by a processor of an electronic device, the electronic device executes the above-mentioned method for generating a core loading plan for a nuclear reactor.

[0438] It should be understood that the present application is not limited to the specific configurations and processes described above and illustrated in the figures. For the sake of brevity, a detailed description of known methods is omitted here. In the above embodiments, several specific steps are described and illustrated as examples. However, the method process of the present application is not limited to the specific steps described and illustrated. Those skilled in the art can make various changes, modifications, and additions, or change the order of the steps after understanding the spirit of the present application.

[0439] The functional blocks shown in the above-described block diagram can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, they can be, for example, electronic circuits, application specific integrated circuits (ASICs), appropriate firmware, plug-ins, function cards, etc. When implemented in software, the elements of the present application are programs or code segments that are used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or a communication link by a data signal carried in a carrier wave. "Machine-readable medium" can include any medium that can store or transmit information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROMs, flash memories, erasable ROMs (EROMs), floppy disks, CD-ROMs, optical disks, hard disks, optical fiber media, radio frequency (RF) links, etc. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

[0440] It should also be noted that the exemplary embodiments mentioned in this application describe some methods or systems based on a series of steps or devices. However, this application is not limited to the order of the above steps. In other words, the steps can be performed in the order mentioned in the embodiments, or in a different order, or several steps can be performed simultaneously.

[0441] Aspects of the present disclosure have been described above with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer or other programmable data processing device to produce a machine so that these instructions executed by the processor of the computer or other programmable data processing device enable the implementation of the function / action specified in one or more boxes of the flowchart and / or block diagram. Such a processor can be, but is not limited to, a general-purpose processor, a special-purpose processor, a special application processor or a field programmable logic circuit. It is also understood that each box in the block diagram and / or flowchart and the combination of the boxes in the block diagram and / or flowchart can also be implemented by dedicated hardware that performs the specified function or action, or can be implemented by a combination of dedicated hardware and computer instructions.

[0442] The above description is only a specific embodiment of the present application. Those skilled in the art will clearly understand that for the convenience and brevity of description, the specific working processes of the systems, modules and units described above can refer to the corresponding processes in the aforementioned method embodiments, and will not be repeated here. It should be understood that the scope of protection of the present application is not limited thereto. Any person skilled in the art can easily think of various equivalent modifications or replacements within the technical scope disclosed in the present application, and these modifications or replacements should be included in the scope of protection of the present application.

Claims

1. A method for generating a core loading plan for a nuclear reactor, characterized in that: The method comprises: Obtaining a first component position set of a core of a nuclear reactor; the first component position set comprising a plurality of component positions; determining a first assembly position corresponding to the shutdown rod assembly of the core from the first assembly position set based on the arguments of the assembly positions in a preset coordinate system; the preset coordinate system is a rectangular coordinate system with the midpoint corresponding to the first assembly position set as the origin; For the remaining component positions other than the first component position in the first component position set, determining at least one symmetric component position set corresponding to the remaining component positions; each of the symmetric component position sets in the at least one symmetric component position set includes component positions having a symmetric relationship among the remaining component positions; For the at least one symmetrical component position set, determining a component type corresponding to the symmetrical component position set from at least one preset component type; A core loading plan for the core is generated according to the component type corresponding to each component position in the first component position set.

2. The method according to claim 1, characterized in that Determining a first assembly position corresponding to the shutdown rod assembly of the core from the first assembly position set according to the arguments of the assembly positions in a preset coordinate system includes: determining at least one reference assembly position of the shutdown rod assembly from the first assembly position set based on the arguments of the assembly positions in the first assembly position set; Determining, from the first component position set, a first symmetrical position corresponding to each of the at least one reference component position; The at least one reference component position and the first symmetrical position are determined as the first component position.

3. The method according to claim 2, characterized in that Determining at least one reference assembly position of the shutdown rod assembly from the first assembly position set based on the argument of each assembly position in the first assembly position set includes: determining a second component position set from the first component position set according to the arguments of the component positions in the first component position set, wherein the second component position set includes component positions having an argument of a first preset argument; Determine the component position closest to the circumference of a preset circle in the second component position set as the first selected position; the preset circle has the origin of the preset coordinate system as its center; The first selected position and the second selected position corresponding to the first selected position are determined as the reference component position; the second selected position is the component position adjacent to the first selected position in the first component position concentration in the direction close to the origin of the preset coordinate system.

4. The method according to claim 3, characterized in that Before determining the component position closest to the circumference of the preset circle in the second component position concentration as the first selected position, the method further includes: In the second component position set, the component position farthest from the origin of the preset coordinate system is determined as the third selected position; Obtaining a first polar coordinate radius corresponding to the third selected position in the preset coordinate system; According to the principle of bisection of volume, determining the target radius of the preset circle according to the first polar coordinate radius; A circle is drawn with the origin of the preset coordinate system as the center and the target radius as the radius to obtain a preset circle.

5. The method according to claim 1, wherein The determining, for the remaining component positions other than the first component position in the first component position set, at least one symmetric component position set corresponding to the remaining component positions includes: determining a third component position set from the remaining component positions based on the arguments of each component position in the preset coordinate system, the third component position set including component positions having an argument greater than or equal to the second preset argument and less than or equal to the third preset argument; For each component position in the third component position set, determining at least one second symmetrical position having a symmetrical relationship with the component position from the remaining component positions; For each component position in the third component position set, the symmetric component position set is obtained according to the component position and the at least one second symmetric position corresponding to the component position.

6. The method according to claim 5, characterized in that The determining a third component position set from the remaining component positions according to the arguments of each component position in the preset coordinate system includes: determining, based on the arguments of each of the remaining component positions in the preset coordinate system, a fourth component position set from the remaining component positions, the fourth component position set including component positions having corresponding arguments greater than or equal to the second preset argument and less than or equal to the fourth preset argument, and an angle between the third preset argument and the second preset argument being greater than an angle between the fourth preset argument and the second preset argument; A fifth component position set is determined from the remaining component positions based on the component positions in the fourth component position set; the component positions in the fifth component position set and the component positions in the fourth component position set are axially symmetric about a target axis, where the target axis is an axis passing through the origin of the preset coordinate system and having an angle equal to the fourth preset argument. The third component location set includes the fourth component location set and the fifth component location set.

7. The method according to claim 1, characterized in that The obtaining of a first component position set of a core of a nuclear reactor includes: Determining, based on core parameters of the nuclear reactor, component distribution dimensions of the core of the nuclear reactor, wherein the component distribution dimensions are the number of rows and columns in which components within the core are arranged; In a preset coordinate system, constructing a component distribution array of the core according to the component distribution dimension; the component distribution array is composed of a plurality of component positions, the number of rows and columns of the component distribution array is the same as the component distribution dimension, and the center of the component distribution array is the origin of the preset coordinate system; screening at least one fourth selected position from the component distribution array according to the coordinates of each component position in the component distribution array in the preset coordinate system; The fourth selected position is within a preset core radius; The first component position set is composed using the at least one fourth selected position.

8. The method according to claim 7, characterized in that The core parameters include component center distance, component box outer radius, and core radius. The selecting at least one fourth selected position from the component distribution array based on the coordinates of each component position in the component distribution array in the preset coordinate system includes: For each component position, determining the coordinates of the component position based on the component center distance, the component distribution dimension, and the target coordinates of the component position; the target coordinates represent the row and column number of the component position in the component distribution array; For each component position, determining a target distance corresponding to the component position based on the coordinates of the component position and the outer radius of the component box; the target distance is the distance between the boundary of the component box at the component position and the origin of the preset coordinate system; For each of the assembly positions, if the target distance of the assembly position is smaller than the core radius, the assembly position is determined as the first assembly position.

9. The method according to claim 1, characterized in that Generating the core loading plan of the core according to the component type corresponding to each component position in the first component position set includes: Obtaining a first serial number and a second serial number corresponding to each component position in the first component position set; the first serial number corresponding to each component position is the serial number of the component position in the component distribution array of the core, and the second serial number corresponding to each component position is the serial number of the component position in the first component position set; determining an external component position reference set from the first component position set based on the first sequence number and the second sequence number; wherein the first sequence number corresponding to the component position in the external component position reference set is equal to the sum of the minimum first sequence number and a target difference, where the target difference is one less than the second sequence number corresponding to the component position; Determine a symmetric component position set corresponding to each component position in the external component position reference set as an external component position set; determining an internal component location set based on component locations other than the external component location set in the first component location set; determining, based on coordinates of each component position in a target component position set in a preset coordinate system, second component positions of the control rod around each component position; the target component position set is the internal component position set or the external component position set, and the number of second component positions around the component position in the internal component position set is greater than the number of second component positions around the component position in the external component position set; A core loading plan for the core is generated according to the component type corresponding to each component position in the first component position set and the second component position.

10. The method according to claim 9, characterized in that Determining the second assembly position of the control rod around each assembly position according to the coordinates of each assembly position in the target assembly position set in a preset coordinate system includes: If the target component position set is the internal component position set, then for each component position in the internal component position set, determining a fifth selected position of the control rod around the component position based on the coordinates of the component position in a preset coordinate system, a third sequence number of the component position in the internal component position set, and a preset sequence number; the preset sequence number ranges from 1 to a preset number of control rods, where the preset number of control rods is a preset maximum number of control rods around the component position; For each component position in the internal component position set, obtain the polar coordinate radius of the fifth selected position corresponding to the component position; if the polar coordinate radius of the fifth selected position is smaller than the second polar coordinate radius, determine the fifth selected position as the second component position; the second polar coordinate radius is the polar coordinate radius of the component position corresponding to the largest first serial number in the external component position reference set.

11. The method according to claim 9, characterized in that Determining the second assembly positions of the control rods around each assembly position according to the coordinates of each assembly position in the first assembly position set in a preset coordinate system includes: If the target assembly position set is the external assembly position set, then for each assembly position in the external assembly position reference set, determining a sixth selected position of the control rod around the assembly position according to a coordinate of the assembly position in a preset coordinate system, a fourth sequence number of the assembly position in the external assembly position set, and a preset sequence number; For each component position in the first target position set of the external component position set, obtaining the polar coordinate radius of the sixth selected position corresponding to the component position, and determining a first preset number of seventh selected positions in the sixth selected positions corresponding to the component position as the second component position; the first target position set includes a symmetric component position set corresponding to the component positions corresponding to the smallest first sequence number and the largest first sequence number in the external component position reference set, and the polar coordinate radius of the seventh selected position is smaller than the polar coordinate radius of the component position in the sixth selected position that is not the seventh selected position; For each component position in the second target position set in the external component position set, obtain the polar coordinate radius of the sixth selected position corresponding to the component position, and determine the second preset number of eighth selected positions in the sixth selected positions corresponding to the component position as the second component position; the second target position set includes the component positions in the external component position set except the first target position set, the polar coordinate radius of the eighth selected position is smaller than the polar coordinate radius of the component positions in the sixth selected positions that are not the eighth selected positions, and the second preset number is larger than the first preset number.

12. The method according to claim 9, characterized in that Generating a core loading plan for the core according to the component type corresponding to each component position in the first component position set and the second component position includes: Obtain the polar coordinate radius corresponding to each position of the second component; Divide the second component positions into sixth component position sets corresponding to G different ranges according to the maximum polar coordinate radius among the polar coordinate radii corresponding to the second component positions and a preset ratio, where G is a positive integer greater than 1; Dividing the target sixth component position set in the G sixth component position sets into a preset number of seventh component position sets; the target sixth component position set is a sixth component position set that does not include the component position corresponding to the origin of the preset coordinate system, and the seventh component position set in the wth target sixth component position set among the G sixth component position sets is separated from the adjacent seventh component position set in the w+1th target sixth component position set by a preset angle, where w ranges from 1 to G-2; A core loading plan for the core is generated based on the component type corresponding to each component position in the first component position set, the second component position, the seventh component position set of the target sixth component position set, and the sixth component position sets other than the target sixth component position set in the G sixth component position sets.

13. A device for generating a core loading plan for a nuclear reactor, characterized in that: The device comprises: An acquisition module, configured to acquire a first component position set of a core of a nuclear reactor; the first component position set comprising a plurality of component positions; an assembly position determination module, configured to determine, from the first assembly position set, a first assembly position corresponding to the shutdown rod assembly of the core based on the arguments of the assembly positions in the first assembly position set in a preset coordinate system; the preset coordinate system being a rectangular coordinate system with a midpoint corresponding to the first assembly position set as its origin; a component position set determining module, configured to determine, for the remaining component positions other than the first component position in the first component position set, at least one symmetric component position set corresponding to the remaining component positions; each of the symmetric component position sets in the at least one symmetric component position set includes component positions having a symmetric relationship among the remaining component positions; a component type determination module, configured to determine, for the at least one symmetrical component position set, a component type corresponding to the symmetrical component position set from at least one preset component type; A plan generating module is used to generate a core loading plan for the core according to the component type corresponding to each component position in the first component position set.

Images