Programmed processing determination method and system for G9 curve

Abstract

The invention relates to a programmed processing determination method and system for a G9 curve, and the method comprises the steps: obtaining an original G9 curve which is a relation curve of a rod position and power and comprises N1 data points, each data point comprises the values of two parameter variables, and the two parameter variables comprise the rod position and the power; n2 data points with linearity meeting a preset condition are screened from N1 data points of the original G9 curve, and N2 is smaller than N1; and determining the new G9 curve according to the N2 data points. According to the method, the new G9 curve is determined by screening the data points meeting the linear conditions through programmed processing, the problems of low efficiency and large deviation caused by subjective judgment of personnel are effectively solved, the accuracy and efficiency of determining the G9 curve are improved, manual errors are reduced, and the automation degree is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of nuclear energy safety and control, and particularly to a method and system for determining the programmed processing of the G9 curve. Background Art

[0002] The present invention mainly focuses on solving the problem of selecting the optimal data points of the G9 curve in the rod control system of a nuclear power plant. During the operation of a nuclear power plant, the G9 curve plays a crucial role. The G9 curve is a two-dimensional array, one dimension is the rod position of the power rod, and the other dimension is the load of the turbogenerator. Usually, for a more vivid display, it is presented in the form of a curve, so it is called the G9 curve. The rod position of the power rod is interpolated and calculated according to the actual load of the turbogenerator set and transmitted to the downstream rod control system for automatic rod lifting and insertion operations. The G9 curve is used to automatically adjust the position of the control rod according to the load of the turbogenerator to maintain the dynamic balance between the core power and the power of the secondary loop steam turbine. Especially in the operation modes of load tracking or rapid load increase and decrease, this balance is particularly important to ensure the safe and stable operation of the reactor. The traditional DCS (Distributed Control System) has certain limitations in the setting of the G9 curve and usually only supports the setting of a limited number of data points. However, the G9 curve obtained through theoretical calculation or on-site test often contains a large number of data points, such as 51 groups or 26 groups of data points, and these data points are evenly spaced to cover the entire power range. However, due to system limitations and risk considerations, it is impossible to set all these data points into the DCS. Therefore, it is necessary to select the most representative data points from the original data for the setting of the new G9 curve to restore the shape and characteristics of the original G9 curve to the greatest extent. In the traditional method, the obvious inflection points on the curve are mainly identified by the naked eye of personnel, which actually involves a large number of combination schemes because a set of data that can best represent the original curve needs to be found from numerous possible selections. The accuracy of the G9 curve directly affects the compensation ability of the control rod for the core reactivity. If there is a deviation in the calculated value of the rod position of the power rod, it may cause the control rod to be unable to fully compensate for the change in the core reactivity, which may lead to a series of serious problems. For example, overheating of the primary loop may cause the unexpected opening of the GCT-A or VDA valves. Among them, the GCT-A valve is part of the steam turbine bypass discharge system, that is, these valves open under unexpected circumstances, thus triggering an unexpected license operation event; and subcooling may trigger the C22 signal, which may lead to a rapid load rejection of the unit. These events may pose a serious threat to the safe operation of the nuclear power plant. Summary of the Invention

[0003] The present invention provides a method and system for determining the programmed processing of the G9 curve to solve the problems of low efficiency and large deviation caused by relying on subjective judgment of personnel.

[0004] The technical solution adopted by the present invention to solve its technical problems is: to provide a method for determining the programmed processing of the G9 curve, which is applied to a processor, and the method includes:

[0005] Step S1: Obtain the original G9 curve, where the original G9 curve is a relationship curve between rod position and power and includes N1 data points, and the data points include the values of two parameter variables, and the two parameter variables include rod position and power;

[0006] Step S2: Select N2 data points with linearity meeting the preset conditions from the N1 data points of the original G9 curve, where N2 is less than N1;

[0007] Step S3: Determine the new G9 curve according to the N2 data points.

[0008] In one embodiment, the step S2 includes:

[0009] Step 2-1: Select two data points from the N1 data points of the original G9 curve as the data points of the new G9 curve;

[0010] Step 2-2: Screen out the data points of the new G9 curve from the remaining N1-2 data points according to the following method:

[0011] Respectively take the two determined data points on the new G9 curve as the head boundary point and the tail boundary point, and determine a linear curve according to the head boundary point and the tail boundary point;

[0012] Determine the data points between the head boundary point and the tail boundary point from the original G9 curve and use them as the intermediate data points;

[0013] Determine the intermediate data point with the largest deviation between the linear curve and the original G9 curve from the intermediate data points and use it as the data point of the new G9 curve.

[0014] In one embodiment, determining the intermediate data point with the largest deviation between the linear curve and the original G9 curve from the intermediate data points includes:

[0015] For each intermediate data point, according to the value of one of the parameter variables of the intermediate data point and the determined linear curve, use an interpolation algorithm to calculate the value of the other parameter variable of the intermediate data point, and calculate the difference between it and the value of the other parameter variable of the intermediate data point in the original G9 curve;

[0016] Determine the intermediate data point with the largest difference from the differences corresponding to the multiple intermediate data points respectively.

[0017] In one embodiment, the step S2-1 includes:

[0018] Select the first and last data points from the original G9 curve as the first and last data points of the new G9 curve.

[0019] In one embodiment, for each intermediate data point, according to the value of one of the variables of the intermediate data point and the determined linear curve, using an interpolation algorithm to calculate the value of the other variable of the intermediate data point, and calculating the difference between it and the value of the other variable of the intermediate data point in the original G9 curve includes:

[0020] For each intermediate data point, according to the rod position value of the intermediate data point and the determined linear curve, using a power interpolation algorithm to calculate the interpolated power value of the intermediate data point, and calculating the difference between it and the power value of the intermediate data point in the original G9 curve;

[0021] In one embodiment, for each intermediate data point, according to the value of one of the variables of the intermediate data point and the determined linear curve, using an interpolation algorithm to calculate the value of the other variable of the intermediate data point, and calculating the difference between it and the value of the other variable of the intermediate data point in the original G9 curve further includes:

[0022] For each intermediate data point, according to the power value of the intermediate data point and the determined linear curve, using a rod position interpolation algorithm to calculate the interpolated rod position value of the intermediate data point, and calculating the difference between it and the rod position value of the intermediate data point in the original G9 curve.

[0023] In one embodiment, the power interpolation algorithm includes:

[0024]

[0025] where P n is the interpolated power value of the nth point of the intermediate data point;

[0026] R 1，n is the rod position value of the nth point in the intermediate data point;

[0027] P 2,1 is the power value of the first boundary data point of the new G9 curve;

[0028] R 2,1 is the rod position value of the first boundary data point of the new G9 curve;

[0029] P 2,2 is the power value of the last boundary data point of the new G9 curve;

[0030] R 2,2 is the rod position value of the last boundary data point of the new G9 curve.

[0031] In one embodiment, step S3 includes:

[0032] Determine whether the number of data points of the new G9 curve reaches N2;

[0033] If it reaches N2, connect the N2 data points to obtain the new G9 curve.

[0034] In one embodiment, N1 is 51 or 26 groups, and N2 is 10 groups.

[0035] The present invention also provides a programmed processing determination system for the G9 curve, including a memory, a processor, and a computer program stored on the memory. The processor executes the computer program to implement the steps of the method described in any one of the above.

[0036] Implementing the present invention has the following beneficial effects: The present invention provides a programmed processing determination method and system for the G9 curve. Among them, the programmed processing determination method for the G9 curve includes obtaining an original G9 curve, where the original G9 curve is a relationship curve between rod position and power and includes N1 data points. The data points include values of two parameter variables, and the two parameter variables include rod position and power; screening out N2 data points with linearity meeting a preset condition from the N1 data points of the original G9 curve, where N2 is less than N1; determining the new G9 curve according to the N2 data points. This application determines the new G9 curve by programmatically processing and screening data points that meet the linear conditions, effectively solving the problems of low efficiency and large deviation caused by relying on subjective judgment of personnel, improving the accuracy and efficiency of G9 curve determination, reducing manual errors, and enhancing the degree of automation. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the technical solutions of the present invention, the drawings required for the description of the present invention will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0038] Figure 1 is a flowchart of the programmed processing determination method for the G9 curve of this application;

[0039] Figure 2 is a schematic diagram of the original G9 curve of this application;

[0040] Figure 3 is a schematic diagram of deviation calculation of this application;

[0041] Figure 4 is a schematic diagram of the selection result of the third point of the new G9 curve of this application;

[0042] Figure 5 It is a schematic diagram of the selection result of the fourth point of the new G9 curve in this application;

[0043] Figure 6 It is a schematic diagram of the final curve selection result of the new G9 curve in this application. Specific Embodiments

[0044] To make the objectives, features, and advantages of the present invention more obvious and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the scope of protection of the present invention.

[0045] It can be understood that the above embodiments only represent the preferred embodiments of the present invention, and the description thereof is relatively specific and detailed, but it should not be construed as a limitation on the scope of the patent of the present invention; it should be noted that for those of ordinary skill in the art, without departing from the concept of the present invention, the above technical features can be freely combined, and several deformations and improvements can also be made, all of which fall within the scope of protection of the present invention; therefore, all equivalent transformations and modifications made to the scope of the claims of the present invention shall fall within the scope covered by the claims of the present invention.

[0046] As Figure 1 and Figure 2 shown, Figure 1 It is a schematic flowchart of the method for determining the programmed processing of the G9 curve;

[0047] This application provides a method for determining the programmed processing of a G9 curve, including the following steps: 1. A method for determining the programmed processing of a G9 curve, which is applied to a processor, and the method includes:

[0048] Step S1: Obtain an original G9 curve, where the original G9 curve is a relationship curve between rod position and power, and includes N1 data points. The data points include values of two parameter variables, and the two parameter variables include rod position and power;

[0049] In this step, it should be noted that the original G9 curve is obtained through a large amount of experimental and actual operation data accumulation, and it depicts the relationship between the control rod position and power. This original curve is obtained from the database of the nuclear power plant and contains N1 data points. Each data point corresponds to a specific rod position and the corresponding power value, and these data points are collected under different working conditions.

[0050] Step S2: Select N2 data points from the N1 data points of the original G9 curve whose linearity meets the preset conditions, where N2 is less than N1;

[0051] Step S3: Determine a new G9 curve based on N2 data points.

[0052] A method for determining the programmed processing of a G9 curve proposed in this embodiment significantly improves the technical effect compared with the traditional manual adjustment and empirical judgment methods. Through programmed processing and data screening, it effectively avoids errors and deviations caused by human factors, thereby being able to more accurately determine the G9 curve and provide a more reliable basis for the safe and stable operation of nuclear power plants. At the same time, the automated method reduces the workload and time cost of manual operations, improves the efficiency of determining the G9 curve, enables nuclear power plants to respond more quickly to changes in operating conditions, and timely adjusts control strategies. In addition, this method can automatically screen and fit curves according to actual data, has strong adaptability, can be applied to control rod systems of nuclear power plants with different models and different operating states, and has a wide range of application prospects. Finally, the newly determined G9 curve has better linearity and continuity as a whole, helps the stable operation of the control rod control system, reduces control fluctuations and instability phenomena caused by uneven curves, and further improves the safety and reliability of nuclear power plant operation.

[0053] As Figure 2 and Figure 3 shown, further, step S2 includes:

[0054] Step 2-1: Select two data points from the N1 data points of the original G9 curve as the data points of the new G9 curve;

[0055] Step 2-2: Screen the data points of the new G9 curve from the remaining N1 - 2 data points according to the following method;

[0056] Respectively take the two determined data points on the new G9 curve as the head boundary point and the tail boundary point, and determine a linear curve according to the head boundary point and the tail boundary point;

[0057] Determine the data points between the head boundary point and the tail boundary point from the original G9 curve and take them as intermediate data points;

[0058] Determine the intermediate data point with the largest deviation between the linear curve and the original G9 curve from the intermediate data points and take it as the data point of the new G9 curve.

[0059] It should be noted that in this step, first, two data points are selected from the N1 data points of the original G9 curve as the data points of the new G9 curve. The first and last data points will be selected as the first and last data points of the new G9 curve. Next, the data points of the new G9 curve are screened from the remaining N1 - 2 data points according to the following method: The two determined data points on the new G9 curve are respectively used as the first boundary point and the last boundary point, and a linear curve is determined based on these two boundary points. The data points between the first boundary point and the last boundary point are determined from the original G9 curve and used as the intermediate data points. Among these intermediate data points, the intermediate data point with the largest deviation between the linear curve and the original G9 curve is determined and used as the data point of the new G9 curve.

[0060] As Figure 4 and Figure 5 shown, specifically, the determination of the intermediate data includes: Assume that the original G9 curve includes (R1, P1), (R2, P2), (R3, P3), (R4, P4), (R5, P5), (R6, P6), (R7, P7), (R8, P8), (R9, P9), (R10, P10), (R11, P11). The first and last data points, the first boundary point and the last boundary point, are selected for the new G9 curve: (R1, P1) and (R11, P11). At this time, the intermediate data are (R2, P2), (R3, P3), (R4, P4), (R5, P5), (R6, P6), (R7, P7), (R8, P8), (R9, P9), (R10, P10). Assume that (R8, P8) is selected as the next data point of the new G9 curve after interpolation calculation and difference calculation. When calculating the next data point of the new G9 curve, the first intermediate data are (R2, P2), (R3, P3), (R4, P4), (R5, P5), (R6, P6), (R7, P7), and the second intermediate data are (R9, P9), (R10, P10). The first intermediate data uses (R1, P1) as the first boundary point and (R8, P8) as the last boundary point, and the second intermediate data uses (R8, P8) as the first boundary point and (R11, P11) as the last boundary point. This process is repeated until the data points satisfying the preset quantity N2 are screened out.

[0061] Specifically, the selected data points are connected in sequence to obtain a new curve. The deviation between the corresponding power value of the original curve at different rod positions and the power value obtained by interpolation calculation of the new curve is repeatedly calculated, and the data point corresponding to the maximum deviation is obtained as the fourth data point in the new curve (marked with a red dot), as specifically Figure 5 shown.

[0062] Furthermore, step S2 - 1 includes:

[0063] Select the first and last data points from the original G9 curve as the first and last data points of the new G9 curve.

[0064] Further, for each intermediate data point, according to the value of one of the variables of the intermediate data point and the determined linear curve, use the interpolation algorithm to calculate the value of the other variable of the intermediate data point, and calculate the difference between it and the value of the other variable of the intermediate data point in the original G9 curve, including:

[0065] For each intermediate data point, according to the rod position value of the intermediate data point and the determined linear curve, use the power interpolation algorithm to calculate the interpolated power value of the intermediate data point, and calculate the difference between it and the power value of the intermediate data point in the original G9 curve;

[0066] And / or, for each intermediate data point, according to the power value of the intermediate data point and the determined linear curve, use the rod position interpolation algorithm to calculate the interpolated rod position value of the intermediate data point, and calculate the difference between it and the rod position value of the intermediate data point in the original G9 curve.

[0067] It should be noted that power interpolation or rod position interpolation can be used, or both power interpolation and rod position interpolation can be used simultaneously. By calculating the difference between the interpolated value and the original value of each intermediate data point, the data point with the largest deviation is found and used as the next data point of the new G9 curve. Specifically, for power interpolation, using the rod position and power values of the known head and tail boundary points, a linear relationship is constructed. According to the rod position value of the intermediate data point, the corresponding interpolated power value is calculated, and then subtracted from the actual power value at this rod position in the original G9 curve to obtain the power deviation; similarly, for rod position interpolation, the interpolated rod position value is calculated based on the power value and compared with the original value. By comparing the deviation values of all intermediate data points, the data point with the largest deviation is selected, and this point is the next data point of the new G9 curve.

[0068] When comparing the deviation distances between each data point and the new curve, one can choose the minimum deviation on the vertical axis, the minimum deviation on the horizontal axis, and the minimum perpendicular distance deviation. There are some differences in the final calculation results after adopting different selection schemes. If the minimum deviation on the vertical axis is adopted, it means that the final G9 curve has the minimum power deviation from the original G9 curve, and during the normal power operation of the nuclear power unit, the power coefficient changes less with the power level. A smaller power deviation means a smaller core reactivity deviation. If the minimum deviation on the horizontal axis is adopted, it means that the control rod position deviation is the smallest, and the differential worth of the control rod varies greatly with different rod positions. A smaller rod position deviation does not represent a smaller reactivity deviation. In this scheme, the minimum deviation on the vertical axis scheme (i.e., power deviation) is selected. Through this scheme, it can be more scientific and reasonable. At the same time, by implementing this algorithm in a programming manner, the selection of the G9 curve can be made faster, and the calculation can be completed within seconds. Compared with the original visual recognition by the naked eye, it shortens the time and improves the accuracy of the G9 curve, ensuring the safe operation of the unit during load tracking and peak shaving operation.

[0069] Furthermore, the power interpolation algorithm includes:

[0070]

[0071] where P n is the interpolation power value of the nth point of the intermediate data points;

[0072] R 1，n is the rod position value of the nth point in the intermediate data points;

[0073] P 2,1 is the power value of the first boundary data point of the new G9 curve;

[0074] R 2,1 is the rod position value of the first boundary data point of the new G9 curve;

[0075] P 2,2 is the power value of the last boundary data point of the new G9 curve;

[0076] R 2,2 is the rod position value of the last boundary data point of the new G9 curve.

[0077] The deviation calculation method is as follows:

[0078] where P 1,n is the power value of the nth point in the original curve.

[0079] In a specific embodiment, the left table is Table 1, which is the original G9 curve provided for the design; the right table is Table 2, which is the final new G9 curve. Visualize the original G9 curve on the coordinate axis as Figure 2(The data is shown in Table 1). The abscissa of the G9 curve is the position of the power rod (Step), and the ordinate is the power of the turbogenerator (WMe).

[0080]

[0081]

[0082] Table 3 shows the calculation process for the selection of the third point.

[0083] Table 3

[0084]

[0085]

[0086] After obtaining the interpolated power values of different data points through the interpolation calculation formula, calculate the power deviation of each data point, and obtain the data point with the largest power deviation as the next data point of the new G9 curve. Through 3, it can be obtained that the data point with a power of 353 and a rod position of 478 is the next data point of the new G9 curve. Table 4 shows the calculation process for the selection of the fourth point.

[0087] Table 4

[0088]

[0089]

[0090] Further, step S3 includes: determining whether the number of data points of the new G9 curve reaches N2;

[0091] If it reaches N2, connect the N2 data points to obtain the new G9 curve. Connecting the N2 data points to obtain the new G9 curve as Figure 6 shown.

[0092] Specifically, N1 is 51 or 26 groups, and N2 is 10 groups.

[0093] This application also provides a programmed processing determination system for the G9 curve, including a memory, a processor, and a computer program stored on the memory. The processor executes the computer program to implement the steps of any one of the above methods.

[0094] This application can be more scientific and reasonable. At the same time, by implementing this algorithm in a programming way, the selection of the G9 curve can be made faster, and the calculation can be completed within seconds. Compared with the original visual recognition, the time is shortened, and the accuracy of the G9 curve is also improved, ensuring the safe operation of the unit during load tracking and peak shaving operation. It can quickly select the most representative 10 groups of data from the original G9 curve data as the final curve, ensuring the minimum deviation between the final curve and the original curve, and thus obtaining the optimal G9 curve.

[0095] It can be understood that the above embodiments only represent the preferred embodiments of the present invention, and the description is relatively specific and detailed, but it cannot be understood as a limitation on the scope of the patent of the present invention; it should be noted that for those of ordinary skill in the art, without departing from the concept of the present invention, the above technical features can be freely combined, and several deformations and improvements can also be made, which all belong to the protection scope of the present invention; therefore, all equivalent transformations and modifications made to the scope of the claims of the present invention shall fall within the scope covered by the claims of the present invention.

Claims

1. A programmed processing determination method for a G9 curve, applied to a processor, characterized in that The method includes: Step S1: Obtain the original G9 curve, which is a curve showing the relationship between rod position and power and includes N1 data points. Each data point includes values of two parameter variables, and the two parameter variables are rod position and power; Step S2: Select N2 data points with linearity meeting a preset condition from the N1 data points of the original G9 curve, where N2 is less than N1; Step S3: Determine the new G9 curve according to the N2 data points.

2. The method for determining the programmed processing of the G9 curve according to claim 1, characterized in that The step S2 includes: Step 2-1: Select two data points from the N1 data points of the original G9 curve as data points of the new G9 curve; Step 2-2: Select data points of the new G9 curve from the remaining N1 - 2 data points in the following manner: Take the two determined data points on the new G9 curve as the head boundary point and the tail boundary point respectively, and determine a linear curve according to the head boundary point and the tail boundary point; Determine the data points between the head boundary point and the tail boundary point from the original G9 curve and take them as intermediate data points; Determine the intermediate data point with the largest deviation between the linear curve and the original G9 curve from the intermediate data points and take it as a data point of the new G9 curve.

3. The method for determining the programmed processing of the G9 curve according to claim 2, characterized in that, Determining the intermediate data point with the largest deviation between the linear curve and the original G9 curve from the intermediate data points includes: For each intermediate data point, according to the value of one of the parameter variables of this intermediate data point and the determined linear curve, use an interpolation algorithm to calculate the value of the other parameter variable of this intermediate data point, and calculate the difference between it and the value of the other parameter variable of this intermediate data point in the original G9 curve; Determine the intermediate data point with the largest difference from the differences corresponding to multiple intermediate data points respectively.

4. The method for determining the programmed processing of the G9 curve according to claim 2, characterized in that, The step S2-1 includes: Select the head and tail data points from the original G9 curve as the head data point and the tail data point of the new G9 curve.

5. The method for determining the programmed processing of the G9 curve according to claim 4, characterized in that, For each intermediate data point, according to the value of one of the variables of this intermediate data point and the determined linear curve, using an interpolation algorithm to calculate the value of the other variable of this intermediate data point, and calculating the difference between it and the value of the other variable of this intermediate data point in the original G9 curve includes: For each intermediate data point, according to the rod position value of this intermediate data point and the determined linear curve, use a power interpolation algorithm to calculate the interpolated power value of this intermediate data point, and calculate the difference between it and the power value of this intermediate data point in the original G9 curve.

6. The method for determining the programmed processing of the G9 curve according to claim 5, characterized in that For each intermediate data point, according to the value of one of the variables of this intermediate data point and the determined linear curve, using an interpolation algorithm to calculate the value of the other variable of this intermediate data point, and calculating the difference between it and the value of the other variable of this intermediate data point in the original G9 curve further includes: For each intermediate data point, according to the power value of this intermediate data point and the determined linear curve, use a rod position interpolation algorithm to calculate the interpolated rod position value of this intermediate data point, and calculate the difference between it and the rod position value of this intermediate data point in the original G9 curve.

7. The method for determining the programmed processing of the G9 curve according to claim 6, wherein The power interpolation algorithm includes: where P n is the interpolation power value of the nth point of the intermediate data point; R 1，n is the rod position value of the nth point among the intermediate data points; P 2,1 is the power value of the first boundary data point of the new G9 curve; R 2,1 is the rod position value of the first boundary data point of the new G9 curve; P 2,2 is the power value of the tail boundary data point of the new G9 curve; R 2,2 is the rod position value of the tail boundary data point of the new G9 curve.

8. The method for determining the programmed processing of the G9 curve according to claim 7, characterized in that The step S3 includes: Determine whether the number of data points of the new G9 curve reaches N2; If it reaches N2, connect the N2 data points to obtain the new G9 curve.

9. The method for determining the programmed processing of the G9 curve according to claim 8, characterized in that, The N1 is 51 or 26 groups, and the N2 is 10 groups.

10. A programmed processing determination system for a G9 curve, comprising a memory, a processor, and a computer program stored on the memory, characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1-9.

Images