Method for Testing the Accuracy of Output Results of Transient Safety Analysis Model for Nuclear Reactors

Abstract

The present application proposes a method for verifying the accuracy of the output results of a transient safety analysis model for a nuclear reactor, including: based on the current operating condition information, obtaining a first discrete data set output by the nuclear reactor transient safety analysis model and a second discrete data set obtained through experiments, and determining the corresponding first equally-spaced time discrete function and second equally-spaced time discrete function respectively; based on the first discrete data set and the first equally-spaced time discrete function, determining a first Fourier function, and based on the second discrete data set and the second equally-spaced time discrete function, determining a second Fourier function; determining a differential amplitude value according to the first Fourier function and the second Fourier function; and determining the accuracy of the first discrete data set based on the differential amplitude value. This technical solution improves the comprehensiveness and accuracy of the evaluation of the output accuracy of the nuclear reactor transient safety analysis program.

Description

Technical Field

[0001] This application relates to the technical field of nuclear reactor safety control, and particularly to a method for testing the accuracy of output results of a transient safety analysis model for a nuclear reactor.

Background Art

[0002] The transient safety analysis of a nuclear reactor relies on the simulation calculations of computer programs. Due to the extremely complex structural design of the nuclear reactor itself and various thermohydraulic phenomena existing during the transient process, for this reason, many constitutive models developed based on separation effect experiments are included in the nuclear reactor transient safety analysis program. Most of these models are obtained by fitting experimental measurement data. Considering the relatively large measurement errors in the experimental data itself and the residuals existing during the fitting process, in short, each single constitutive model has relatively large calculation uncertainties. And there are numerous constitutive models in the nuclear reactor transient safety analysis program, and there are high-order interactions among their respective uncertainties, which makes it difficult for the nuclear reactor transient safety analysis program to obtain accurate calculation results.

[0003] In response to this, in related technologies, conservative correction factors are often introduced into many constitutive models, and the calculation uncertainties of the program are enveloped by using conservative calculation results. Generally, for a single constitutive model, the difference or quotient between the experimental measurement value and the program calculation value can be defined as the uncertainty characterization correction coefficient, and the distribution of scatter data can be evaluated by using parametric statistical methods such as hypothesis testing and frequency estimation, or non-parametric statistical methods such as orthogonal sequences and kernel density estimation.

[0004] However, such methods are used to separately evaluate each constitutive model within the nuclear reactor transient safety analysis program. In other words, they are used to evaluate the calculation accuracy of the local part of the nuclear reactor transient safety analysis program, and cannot conduct an overall analysis of the nuclear reactor transient safety. The obtained evaluation results still have deficiencies in terms of integrity.

[0005] Therefore, how to comprehensively and accurately evaluate the calculation result accuracy of the nuclear reactor transient safety analysis program has become an urgent technical problem to be solved currently.

Summary of the Invention

[0006] The embodiments of this application provide a method for testing the accuracy of output results of a transient safety analysis model for a nuclear reactor, aiming to solve the technical problem in related technologies that it is difficult to comprehensively evaluate the calculation result accuracy of the nuclear reactor transient safety analysis program and affect the evaluation accuracy.

[0007] In a first aspect, the embodiments of this application provide a method for testing the accuracy of output results of a transient safety analysis model for a nuclear reactor, including:

[0008] Obtain the current operating condition information of the nuclear reactor;

[0009] Based on the current operating condition information, output a first discrete data set of the safety performance evaluation indexes of the current operating condition through the nuclear reactor transient safety analysis model, and obtain a second discrete data set obtained from experiments of the safety performance evaluation indexes of the current operating condition;

[0010] Based on the first discrete data set and the second discrete data set, determine the corresponding first equally-spaced time discrete function and second equally-spaced time discrete function respectively, wherein the first equally-spaced time discrete function and the second equally-spaced time discrete function respectively reflect the changes of the operating conditions corresponding to the first discrete data set and the second discrete data set on the same uniform time axis;

[0011] Based on the first discrete data set and the first equally-spaced time discrete function, determine a first Fourier function, and based on the second discrete data set and the second equally-spaced time discrete function, determine a second Fourier function;

[0012] According to the first Fourier function and the second Fourier function, determine the difference amplitude value between the first discrete data set and the second discrete data set;

[0013] Based on the difference amplitude value, determine the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0014] In an embodiment of the present application, optionally, the determining the corresponding first equally-spaced time discrete function and second equally-spaced time discrete function based on the first discrete data set and the second discrete data set includes:

[0015] Perform fitting processing on the first discrete data set to obtain a first fitting function;

[0016] In the first fitting function, determine a third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor;

[0017] Based on the third discrete data set, construct a first warping function; and

[0018] Perform fitting processing on the second discrete data set to obtain a second fitting function;

[0019] In the second fitting function, determine a fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor;

[0020] Based on the fourth discrete data set, construct a second warping function;

[0021] Based on the first warping function and the second warping function, determine the corresponding first equidistant time discrete function and second equidistant time discrete function respectively.

[0022] In an embodiment of the present application, optionally, the determining the corresponding first equidistant time discrete function and second equidistant time discrete function based on the first warping function and the second warping function includes:

[0023] Perform time point registration processing on the first warping function and the second warping function;

[0024] Perform uniform equidistant time discretization processing on the first warping function and the second warping function after time point registration respectively to obtain the corresponding first equidistant time discrete function and second equidistant time discrete function.

[0025] In an embodiment of the present application, optionally, the determining the difference amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function includes:

[0026] Based on the first Fourier function, calculate the first output value corresponding to each first data point in the first discrete data set; and

[0027] Based on the second Fourier function, calculate the second output value corresponding to each second data point in the second discrete data set;

[0028] Based on the first output value corresponding to each first data point and the second output value corresponding to each second data point, determine the difference amplitude value between the first discrete data set and the second discrete data set, where

[0029]

[0030] AA is the difference amplitude value between the first discrete data set and the second discrete data set, k is the total number of the first data points in the first discrete data set or the total number of the second data points in the second discrete data set, where the number of the first data points is the same as and corresponds one-to-one to the number of the second data points, n is the serial number of the first data point and the corresponding second data point, F calc is the first output value corresponding to the nth first data point in the first discrete data set, F exp is the second output value corresponding to the nth second data point in the second discrete data set.

[0031] In an embodiment of the present application, optionally, the determining the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model based on the difference amplitude value includes:

[0032] Determine the preset amplitude value range to which the difference amplitude value belongs;

[0033] Set the preset precision corresponding to the preset amplitude value range as the precision of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0034] In a second aspect, an embodiment of the present application provides an output result precision inspection device for a nuclear reactor transient safety analysis model, including:

[0035] A current working condition acquisition unit, configured to acquire the current working condition information of the nuclear reactor;

[0036] A model output acquisition unit, configured to output a first discrete data set of the safety performance evaluation index of the current working condition through the nuclear reactor transient safety analysis model based on the current working condition information;

[0037] An experimental output acquisition unit, configured to acquire a second discrete data set obtained by experiments of the safety performance evaluation index of the current working condition;

[0038] An equidistant time discretization unit, configured to determine a first equidistant time discretization function and a second equidistant time discretization function corresponding to each of them based on the first discrete data set and the second discrete data set, where the first equidistant time discretization function and the second equidistant time discretization function respectively reflect the changes of the working conditions corresponding to the first discrete data set and the second discrete data set on the same uniform time axis;

[0039] A Fourier transform unit, configured to determine a first Fourier function based on the first discrete data set and the first equidistant time discretization function, and determine a second Fourier function based on the second discrete data set and the second equidistant time discretization function;

[0040] A difference amplitude value calculation unit, configured to determine the difference amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function;

[0041] A precision determination unit, configured to determine the precision of the first discrete data set output by the nuclear reactor transient safety analysis model based on the difference amplitude value.

[0042] In an embodiment of the present application, optionally, the equidistant time discretization unit includes:

[0043] A first fitting unit, configured to perform fitting processing on the first discrete data set to obtain a first fitting function;

[0044] The first characteristic phenomenon extraction unit is used to determine, in the first fitting function, a third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor;

[0045] The first warping function construction unit is used to construct a first warping function based on the third discrete data set;

[0046] The second fitting unit is used to perform fitting processing on the second discrete data set to obtain a second fitting function;

[0047] The second characteristic phenomenon extraction unit is used to determine, in the second fitting function, a fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor;

[0048] The second warping function construction unit is used to construct a second warping function based on the fourth discrete data set;

[0049] The first execution unit is used to determine a first equidistant time discrete function and a second equidistant time discrete function corresponding to each of them based on the first warping function and the second warping function.

[0050] In an embodiment of the present application, optionally, the first execution unit is specifically configured to: perform time point registration processing on the first warping function and the second warping function, and respectively perform uniform equidistant time discretization processing on the first warping function and the second warping function after time point registration to obtain a first equidistant time discrete function and a second equidistant time discrete function corresponding to each of them.

[0051] In an embodiment of the present application, optionally, the difference amplitude value calculation unit includes:

[0052] The first output value calculation unit is used to calculate a first output value corresponding to each first data point in the first discrete data set based on the first Fourier function;

[0053] The second output value calculation unit calculates a second output value corresponding to each second data point in the second discrete data set based on the second Fourier function;

[0054] The second execution unit is used to determine the difference amplitude value between the first discrete data set and the second discrete data set based on the first output value corresponding to each first data point and the second output value corresponding to each second data point, where

[0055]

[0056] AA is the difference amplitude value between the first discrete data set and the second discrete data set, k is the total number of the first data points in the first discrete data set or the total number of the second data points in the second discrete data set, where the number of the first data points is the same as that of the second data points and they are in one-to-one correspondence, n is the serial number of the first data point and the corresponding second data point, F calc is the first output value corresponding to the nth first data point in the first discrete data set, F exp is the second output value corresponding to the nth second data point in the second discrete data set.

[0057] In an embodiment of the present application, optionally, the accuracy determination unit is specifically configured to: determine the preset amplitude value range to which the difference amplitude value belongs, and set the preset accuracy corresponding to the preset amplitude value range as the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0058] In a third aspect, an embodiment of the present application provides a computer device, including: at least one processor; and a memory communicatively connected to the at least one processor; wherein, the memory stores instructions executable by the at least one processor, and the instructions are set to execute the method described in the first aspect above.

[0059] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium storing computer-executable instructions for executing the method described in the first aspect above.

[0060] In the above technical solution, in view of the technical problem that it is difficult to comprehensively evaluate the calculation result accuracy of the nuclear reactor transient safety analysis program in the related art, which affects the evaluation accuracy, the output result of the nuclear reactor transient safety analysis model and the actual experimental output result are fitted and then time-registered, and then the Fourier transform is performed to obtain the respective Fourier functions corresponding to the two. Finally, the difference amplitude value between the two is calculated through the respective Fourier functions corresponding to the two, as the ultimate basis for determining the accuracy of the output result of the nuclear reactor transient safety analysis model. Thus, it replaces the technical solution of separately evaluating each constitutive model in the nuclear reactor transient safety analysis program, and can directly evaluate the nuclear reactor transient safety analysis program as a whole based on the output result of the nuclear reactor transient safety analysis program. Compared with the original calculation accuracy evaluation method of separately evaluating each constitutive model, the calculation amount of evaluating the output accuracy of the nuclear reactor transient safety analysis program is greatly reduced, and the comprehensiveness and accuracy of evaluating the output accuracy of the nuclear reactor transient safety analysis program are improved, which helps to enhance the working safety of the nuclear reactor and reduce the risks in the production process.

Description of the Drawings

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

[0062] Figure 1 Shows a flowchart of a method for verifying the accuracy of the output results of a nuclear reactor transient safety analysis model according to an embodiment of the present application;

[0063] Figure 2 Shows a flowchart of a method for verifying the accuracy of the output results of a nuclear reactor transient safety analysis model according to another embodiment of the present application;

[0064] Figure 3 Shows a block diagram of a computer device according to an embodiment of the present application;

[0065] Figure 4 Shows a block diagram of a computer device according to another embodiment of the present application.

Detailed Implementation Modes

[0066] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are some, but not all, of the embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present application without creative efforts fall within the scope of protection of the present application.

[0067] Figure 1 Shows a flowchart of a method for verifying the accuracy of the output results of a nuclear reactor transient safety analysis model according to an embodiment of the present application.

[0068] As Figure 1 shown, a method for verifying the accuracy of the output results of a nuclear reactor transient safety analysis model according to an embodiment of the present application includes:

[0069] Step 102, obtaining the current operating condition information of the nuclear reactor.

[0070] The current operating condition information of the nuclear reactor reflects the performance of the nuclear reactor at the current moment in each dimension affecting its transient safety. Various current operating condition information can be understood as the information in each dimension affecting the transient safety of the nuclear reactor, including but not limited to the initial conditions of the nuclear reactor operation (such as initial power, initial pressure), boundary conditions (such as geometric structure, heat transfer boundary), and the values of state parameters (such as coolant temperature, cladding temperature), etc.

[0071] Step 104: Based on the current operating condition information, output a first discrete data set of the safety performance evaluation indexes corresponding to the current operating condition through a nuclear reactor transient safety analysis model, and obtain a second discrete data set of the safety performance evaluation indexes obtained from experiments for the current operating condition.

[0072] The nuclear reactor transient safety analysis model takes the current operating condition information of the nuclear reactor as input and calculates and outputs the evaluation result of the nuclear reactor transient safety. Specifically, it is the calculation result of the safety performance evaluation indexes of the nuclear reactor transient safety. Generally speaking, the calculation result of the safety performance evaluation indexes of the nuclear reactor transient safety is discrete data points corresponding to different times, that is, the first discrete data set. At the same time, obtain the calculation result of the safety performance evaluation indexes of the nuclear reactor transient safety obtained from experiments. The output result of the nuclear reactor transient safety analysis model and the result obtained from experiments respectively reflect the evaluation result of the nuclear reactor transient safety in terms of model output and experimental output.

[0073] Step 106: Based on the first discrete data set and the second discrete data set, determine the corresponding first equidistant time discrete function and second equidistant time discrete function respectively.

[0074] Among them, the first equidistant time discrete function and the second equidistant time discrete function respectively reflect the changes of the operating conditions corresponding to the first discrete data set and the operating conditions corresponding to the second discrete data set on the same uniform time axis.

[0075] Specifically, first, the first discrete data set can be fitted to obtain a first fitting function, and in the first fitting function, determine the third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor. Based on the third discrete data set, construct a first warping function.

[0076] Conventional registration techniques include translational registration and landmark registration. Among them, translational registration is applicable to the case where the time processes are consistent and cannot consider the influence of events accelerating or decelerating on the time scale, while landmark registration can regularize time according to the specified landmarks and has high applicability. Landmarks, that is, feature points, refer to the feature points with representative significance on the curve of parameter change with time, such as maximum values, extreme values, inflection points, etc. These feature points often have certain physical meanings, such as representing the occurrence or disappearance of a certain thermohydraulic phenomenon. In landmark registration technology, it is first necessary to determine each landmark on the output curve related to the target time, and the given landmark is a hyperparameter, which is often artificially determined by the analyst according to needs.

[0077] Marker registration relies on aligning the feature points of different curves using a warping function. The functional output of the program before registration is y(t). After transforming the time t through the warping function h(t), the registered function y*(t) can be obtained, that is:

[0078] y*(t) = y[h(t)]

[0079] The warping function h(t) has many properties, such as being strictly monotonic, the regularized times corresponding to the time of the marker points on different curves being equal, and the start time and end time being irregular, etc.

[0080] Among the curves fitted from the first discrete data set, select the occurrence time of the characteristic phenomenon of the nuclear reactor, and use the data points corresponding to the occurrence time as the new third discrete data set to construct the first warping function. This first warping function reflects the variation of the simulated working conditions of the characteristic phenomenon occurring in the nuclear reactor simulated by the nuclear reactor transient safety analysis model over time.

[0081] At the same time, the second discrete data set can be fitted to obtain a second fitting function, and in the second fitting function, determine the fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor, and then construct a second warping function based on the fourth discrete data set.

[0082] That is to say, among the curves fitted from the second discrete data set, select the occurrence time of the characteristic phenomenon of the nuclear reactor, and use the data points corresponding to the occurrence time as the new fourth discrete data set to construct the second warping function. This second warping function reflects the variation of the actual working conditions of the characteristic phenomenon occurring in the nuclear reactor in actual experiments over time.

[0083] Then, based on the first warping function and the second warping function, determine the corresponding first equally spaced time discrete function and second equally spaced time discrete function respectively. Since the time distributions of the first warping function and the second warping function are not uniform, in order to increase the calculation accuracy, the time lines of the two can be registered.

[0084] Specifically, first perform time point registration processing on the first warping function and the second warping function, and then perform uniform equally spaced time discretization processing on the first warping function and the second warping function after time point registration respectively to obtain the corresponding first equally spaced time discrete function and second equally spaced time discrete function. Thus, the time points on the obtained first equally spaced time discrete function and second equally spaced time discrete function correspond one by one, and the actual time represented by any time period on the curves corresponding to the two functions is also the same.

[0085] Step 108: Determine the first Fourier function based on the first discrete data set and the first equally-spaced time discretization function, and determine the second Fourier function based on the second discrete data set and the second equally-spaced time discretization function.

[0086] Perform Fourier transforms on the first discrete data set and the second discrete data set respectively. The obtained Fourier functions, compared with the original equally-spaced time discretization functions, decompose the working conditions represented by the original discrete data at a deeper level, and can more deeply reflect the performance levels of the working conditions represented by the first discrete data set and the second discrete data set.

[0087] Step 110: Determine the differential amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function.

[0088] So far, based on the Fourier functions of the first discrete data set and the second discrete data set respectively, the differential amplitude value reflecting the difference level between the two is obtained. It can be seen that to a certain extent, this differential amplitude value reflects the performance levels of the first discrete data set and the second discrete data set under the same working conditions. Further, it also reflects the specific differences between the two output methods, namely the model output to which the first discrete data set belongs and the experimental output to which the second discrete data set belongs, under the same working conditions. The specific method for calculating the differential amplitude value is as follows.

[0089] First, calculate the first output value corresponding to each first data point in the first discrete data set based on the first Fourier function; and calculate the second output value corresponding to each second data point in the second discrete data set based on the second Fourier function.

[0090] Next, determine the differential amplitude value between the first discrete data set and the second discrete data set based on the first output value corresponding to each first data point and the second output value corresponding to each second data point, where,

[0091]

[0092] AA is the differential amplitude value between the first discrete data set and the second discrete data set, k is the total number of the first data points in the first discrete data set or the total number of the second data points in the second discrete data set, where the number of the first data points is the same as that of the second data points and they correspond one by one, n is the serial number of the first data point and the corresponding second data point, F calc is the first output value corresponding to the nth first data point in the first discrete data set, F expis the second output value corresponding to the nth second data point in the second discrete data set.

[0093] It can be seen from this that the differential amplitude value actually reflects the proportion of the total difference between the Fourier function output values of the corresponding data points in the first discrete data set and the second discrete data set in the sum of the Fourier function output values of the corresponding data points in the second discrete data set.

[0094] Step 112, based on the differential amplitude value, determine the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0095] Finally, based on the differential amplitude value that reflects the difference between the model output and the experimental output under the same working conditions, the accuracy of the output content of the nuclear reactor transient safety analysis model compared with the actual experimental output information can be deduced.

[0096] Specifically, the method for calculating this accuracy is: determine the preset amplitude value range to which the differential amplitude value belongs, and set the preset accuracy corresponding to the preset amplitude value range as the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0097] Optionally, each preset amplitude value range can be set to correspond to an accuracy level. After determining the preset amplitude value range to which the differential amplitude value belongs, set the accuracy level corresponding to the preset amplitude value range to which the differential amplitude value belongs as the accuracy of the first discrete data set.

[0098] The above technical solution performs time registration after fitting the output result of the nuclear reactor transient safety analysis model and the actual experimental output result, then performs Fourier transform to obtain the respective corresponding Fourier functions of the two, and finally calculates the differential amplitude value of the two through the respective corresponding Fourier functions of the two as the final basis for determining the accuracy of the output result of the nuclear reactor transient safety analysis model. Thus, it replaces the technical solution of separately evaluating each constitutive model in the nuclear reactor transient safety analysis program in the related art, and can directly evaluate the nuclear reactor transient safety analysis program based on the output result of the nuclear reactor transient safety analysis program. Compared with the original calculation accuracy evaluation method for separately evaluating each constitutive model, it greatly reduces the calculation amount of the output accuracy evaluation of the nuclear reactor transient safety analysis program, improves the comprehensiveness and accuracy of the output accuracy evaluation of the nuclear reactor transient safety analysis program, helps to enhance the working safety of the nuclear reactor, and reduces the risks in the production process.

[0099] Figure 2 Shows a flowchart of a method for testing the accuracy of the output result of a nuclear reactor transient safety analysis model according to another embodiment of the present application.

[0100] As Figure 2 shown, the method for verifying the accuracy of the output result of the nuclear reactor transient safety analysis model according to another embodiment of the present application includes the following steps:

[0101] First, determine the target operating condition for applying the nuclear reactor transient safety analysis model, obtain the input conditions in this operating condition, including but not limited to initial conditions, boundary conditions, and state parameter values, and determine the performance indicators to be evaluated by the nuclear reactor transient safety analysis model.

[0102] Next, on the one hand, conduct experiments to obtain experimental measurement data - the time-series discrete data points of the experimental measurement of this performance indicator, and construct functional data based on the time-series discrete data points of the experimental measurement to obtain the functional object of the experimental measurement data.

[0103] On the other hand, simulate the input conditions of the target operating condition through the nuclear reactor transient safety analysis model to obtain program simulation data, and then construct functional data based on the program simulation data to obtain the functional object of the program simulation data.

[0104] So far, perform dynamic time warping on the functional object of the experimental measurement data and the functional object of the program simulation data to obtain the registered data.

[0105] Specifically, determine the time points when characteristic phenomena occur in the functional objects of the experimental measurement data and the program simulation data, construct the corresponding warping function based on this time point, and then register the time points when characteristic phenomena occur with the respective warping functions of the experimental measurement data and the program simulation data based on the dynamic time warping technology to make the time processes of the two consistent. Further, perform uniform equidistant time discretization on the registered result.

[0106] Next, use the fast Fourier transform to fit the data after uniform equidistant time discretization, construct the Fourier functions of the two, and calculate the difference amplitude value AA between the two based on the Fourier functions of the two. The calculation method of AA is the same as that in the embodiment Figure 1 shown, and will not be elaborated here.

[0107] Finally, based on this difference amplitude value, determine the accuracy of the program simulation data output by the nuclear reactor transient safety analysis model. The smaller the AA value, the closer the program simulation data output by the nuclear reactor transient safety analysis model is to the experimental measurement data, and the higher the calculation result accuracy; the larger the AA value, the greater the deviation between the program simulation data output by the nuclear reactor transient safety analysis model and the experimental measurement data, and the lower the calculation result accuracy.

[0108] In summary, this application designs a temporal data difference quantification and evaluation technique based on the fast Fourier transform and dynamic time warping to achieve the accuracy evaluation of the temporal calculation results of the nuclear reactor transient safety analysis model. Among them, first for a specific target condition to be evaluated, the initial input conditions of the target condition are determined, including initial conditions, boundary conditions, and state parameter values. Subsequently, through conducting experiments and establishing a program simulation model respectively, the experimental measurement temporal data results of specific performance indicators under the target condition and the corresponding program simulation temporal data are obtained. After that, functional data corresponding to the experimental measurement temporal data results and the program simulation temporal data are constructed through a basis function system, and the functional data objects corresponding to the experimental measurement temporal data results and the program simulation temporal data are registered through the dynamic time warping technique. Finally, the fast Fourier transform is used to compare the differences between the two functional data objects and quantify the differences to achieve the purpose of evaluating the accuracy of the temporal calculation results of the nuclear reactor transient safety analysis model.

[0109] The technical solution of this application can quantify the calculation ability and accuracy of the nuclear reactor transient analysis model, solve the possible calculation inaccuracies introduced when relying on subjective will to select programs during traditional safety analysis, and at the same time provide a numerical reference metric for the optimization, development, evaluation, and safety review certification of the thermal-hydraulic system program for nuclear safety analysis. It can clarify the direction for improving the program accuracy, reduce the program R & D cost while promoting the update and iteration of the program, and increase the reliability of the nuclear safety analysis results.

[0110] An embodiment of this application provides an output result accuracy inspection device for a nuclear reactor transient safety analysis model, including:

[0111] The current condition acquisition unit is used to acquire the current condition information of the nuclear reactor;

[0112] The model output acquisition unit is used to output the first discrete data set of the safety performance evaluation indicators related to the current condition through the nuclear reactor transient safety analysis model based on the current condition information;

[0113] The experimental output acquisition unit is used to acquire the second discrete data set obtained from the experiment of the safety performance evaluation indicators related to the current condition;

[0114] The equidistant time discretization unit is used to determine the corresponding first equidistant time discretization function and second equidistant time discretization function based on the first discrete data set and the second discrete data set, where the first equidistant time discretization function and the second equidistant time discretization function respectively reflect the changes of the conditions corresponding to the first discrete data set and the second discrete data set on the same uniform time axis;

[0115] A Fourier transform unit for determining a first Fourier function based on the first discrete data set and the first equally-spaced time discretization function, and determining a second Fourier function based on the second discrete data set and the second equally-spaced time discretization function;

[0116] A differential amplitude value calculation unit for determining the differential amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function;

[0117] An accuracy determination unit for determining the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model based on the differential amplitude value.

[0118] In an embodiment of the present application, optionally, the equally-spaced time discretization unit includes:

[0119] A first fitting unit for performing a fitting process on the first discrete data set to obtain a first fitting function;

[0120] A first characteristic phenomenon extraction unit for determining a third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor in the first fitting function;

[0121] A first warping function construction unit for constructing a first warping function based on the third discrete data set;

[0122] A second fitting unit for performing a fitting process on the second discrete data set to obtain a second fitting function;

[0123] A second characteristic phenomenon extraction unit for determining a fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor in the second fitting function;

[0124] A second warping function construction unit for constructing a second warping function based on the fourth discrete data set;

[0125] A first execution unit for determining a first equally-spaced time discretization function and a second equally-spaced time discretization function corresponding thereto based on the first warping function and the second warping function.

[0126] In an embodiment of the present application, optionally, the first execution unit is specifically configured to: perform time point registration processing on the first warping function and the second warping function, and perform uniform equally-spaced time discretization processing on the time point-registered first warping function and second warping function respectively to obtain a first equally-spaced time discretization function and a second equally-spaced time discretization function corresponding thereto.

[0127] In one embodiment of the present application, optionally, the differential amplitude value calculation unit includes:

[0128] A first output value calculation unit, configured to calculate a first output value corresponding to each first data point in the first discrete data set based on the first Fourier function;

[0129] A second output value calculation unit, configured to calculate a second output value corresponding to each second data point in the second discrete data set based on the second Fourier function;

[0130] A second execution unit, configured to determine a differential amplitude value between the first discrete data set and the second discrete data set based on the first output value corresponding to each first data point and the second output value corresponding to each second data point, where

[0131]

[0132] AA is the differential amplitude value between the first discrete data set and the second discrete data set, k is the total number of the first data points in the first discrete data set or the total number of the second data points in the second discrete data set, where the number of the first data points is the same as and in one-to-one correspondence with the number of the second data points, n is the serial number of the first data point and the corresponding second data point, F calc is the first output value corresponding to the nth first data point in the first discrete data set, and F exp is the second output value corresponding to the nth second data point in the second discrete data set.

[0133] In one embodiment of the present application, optionally, the accuracy determination unit is specifically configured to: determine a preset amplitude value range to which the differential amplitude value belongs, and set the preset accuracy corresponding to the preset amplitude value range as the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0134] The device uses the solution described in any one of the above embodiments, and thus has all the above technical effects, which will not be elaborated here.

[0135] In addition, in one embodiment, the present application provides a computer device, which may be a server, and its internal structure diagram may be as Figure 3As shown in the figure. The computer device includes a processor, a memory, a network interface, and a database connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile and / or volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external client through a network connection. When the computer program is executed by the processor, it can implement the method described in any of the above embodiments.

[0136] In one embodiment, the present application also provides a computer device, which can be a client, and its internal structure diagram can be as shown in Figure 4 As shown in the figure. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external server through a network connection. When the computer program is executed by the processor, it can implement the method described in any of the above embodiments.

[0137] Any of the above computer devices in the embodiments of the present application exists in various forms, including but not limited to:

[0138] (1) Mobile communication devices: The characteristic of this type of device is that it has mobile communication functions and mainly aims to provide voice and data communication. This type of terminal includes: smart phones (such as iPhone), multimedia phones, functional phones, and low-end phones, etc.

[0139] (2) Ultra-mobile personal computer devices: This type of device belongs to the category of personal computers, has computing and processing functions, and generally also has the characteristic of mobile Internet access. This type of terminal includes: PDA, MID, and UMPC devices, etc., such as iPad.

[0140] (3) Portable entertainment devices: This type of device can display and play multimedia content. This type of device includes: audio and video players (such as iPod), handheld game consoles, e-books, and smart toys, wearable devices, and portable in-vehicle navigation devices.

[0141] (4) Server: A device that provides computing services. The server consists of a processor, hard disk, memory, system bus, etc. The server is similar to a general computer architecture, but due to the need to provide highly reliable services, it has higher requirements in terms of processing power, stability, reliability, security, scalability, manageability, etc.

[0142] (5) Other electronic devices with data interaction functions.

[0143] In addition, an embodiment of the present application provides a computer-readable storage medium storing computer-executable instructions, and the computer-executable instructions are used to perform the following steps:

[0144] Obtain the current operating condition information of the nuclear reactor;

[0145] Based on the current operating condition information, output a first discrete data set of safety performance evaluation indicators for the current operating condition through a nuclear reactor transient safety analysis model, and obtain a second discrete data set obtained from experiments of safety performance evaluation indicators for the current operating condition;

[0146] Based on the first discrete data set and the second discrete data set, determine the corresponding first equally spaced time discrete function and second equally spaced time discrete function respectively, where the first equally spaced time discrete function and the second equally spaced time discrete function respectively reflect the changes of the operating conditions corresponding to the first discrete data set and the second discrete data set on the same uniform time axis;

[0147] Based on the first discrete data set and the first equally spaced time discrete function, determine a first Fourier function, and based on the second discrete data set and the second equally spaced time discrete function, determine a second Fourier function;

[0148] According to the first Fourier function and the second Fourier function, determine the difference amplitude value between the first discrete data set and the second discrete data set;

[0149] Based on the difference amplitude value, determine the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

[0150] It should be noted that the functions or steps that the above computer-readable storage medium or computer device can achieve can be referred to the relevant descriptions in the foregoing method embodiments. To avoid repetition, they will not be described one by one here.

[0151] The technical solution of the present application has been described in detail above in conjunction with the accompanying drawings. Through the technical solution of the present application, after fitting the output results of the nuclear reactor transient safety analysis model and the actual experimental output results, time registration is performed, and then Fourier transform is executed to obtain the respective Fourier functions of the two. Finally, the difference amplitude value between the two is calculated through the respective Fourier functions of the two, as the ultimate basis for determining the accuracy of the output results of the nuclear reactor transient safety analysis model. Thus, it replaces the technical solution of separately evaluating each constitutive model in the nuclear reactor transient safety analysis program in the related art, and can directly evaluate the nuclear reactor transient safety analysis program based on the output results of the nuclear reactor transient safety analysis program. Compared with the original calculation accuracy evaluation method of separately evaluating each constitutive model, it greatly reduces the calculation amount of the output accuracy evaluation of the nuclear reactor transient safety analysis program, improves the comprehensiveness and accuracy of the output accuracy evaluation of the nuclear reactor transient safety analysis program, helps to enhance the working safety of the nuclear reactor, and reduces the risks in the production process.

[0152] It should be understood that although the terms first, second, etc. may be used in the embodiments of the present application to describe discrete data sets, these discrete data sets should not be limited to these terms. These terms are only used to distinguish discrete data sets from each other. For example, without departing from the scope of the embodiments of the present application, the first discrete data set may also be referred to as the second discrete data set, and similarly, the second discrete data set may also be referred to as the first discrete data set.

[0153] Depending on the context, the word "if" as used herein may be interpreted as "when" or "while" or "in response to determining" or "in response to detecting". Similarly, depending on the context, the phrase "if determined" or "if detected (stated condition or event)" may be interpreted as "when determined" or "in response to determining" or "when detected (stated condition or event)" or "in response to detecting (stated condition or event)".

[0154] The terms used in the embodiments of the present application are for the purpose of describing specific embodiments only and are not intended to limit the present application. The singular forms "a", "the", and "said" used in the embodiments of the present application and the appended claims are also intended to include the plural forms unless the context clearly indicates otherwise.

[0155] In several embodiments provided in this application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling, direct coupling, or communication connection to each other can be through some interfaces. The indirect coupling or communication connection of the devices or units can be in electrical, mechanical, or other forms.

[0156] In addition, in each embodiment of this application, the functional units can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above integrated units can be implemented in the form of hardware or in the form of a combination of hardware and software functional units.

[0157] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, storage, database, or other medium used in each embodiment provided in this application can include non-volatile and / or volatile memories. Non-volatile memories can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memories can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0158] The above-described embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be included within the protection scope of the present application.

Claims

1. A method for checking the accuracy of output results of a nuclear reactor transient safety analysis model, characterized in that: include: Obtain current operating status information of nuclear reactors; Based on the current operating condition information, output a first discrete data set of safety performance evaluation indicators related to the current operating condition through a nuclear reactor transient safety analysis model, and obtain a second discrete data set of safety performance evaluation indicators related to the current operating condition obtained by experiment; Based on the first discrete data set and the second discrete data set, determining a first equidistant time discrete function and a second equidistant time discrete function corresponding to each other, wherein the first equidistant time discrete function and the second equidistant time discrete function respectively reflect changes in the operating conditions corresponding to the first discrete data set and the operating conditions corresponding to the second discrete data set on the same uniform time axis; Determine a first Fourier function based on the first discrete data set and the first equidistant time discrete function, and determine a second Fourier function based on the second discrete data set and the second equidistant time discrete function; Determine a difference amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function; Based on the difference amplitude value, the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model is determined.

2. The method according to claim 1, characterized in that The determining, based on the first discrete data set and the second discrete data set, respectively corresponding first equidistant time discrete function and second equidistant time discrete function comprises: Performing fitting processing on the first discrete data set to obtain a first fitting function; In the first fitting function, determining a third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor; constructing a first warping function based on the third discrete data set; and Performing fitting processing on the second discrete data set to obtain a second fitting function; In the second fitting function, a fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor is determined; constructing a second warping function based on the fourth discrete data set; Based on the first warping function and the second warping function, respectively corresponding first equidistant time-discrete functions and second equidistant time-discrete functions are determined.

3. The method according to claim 2, characterized in that The determining, based on the first warping function and the second warping function, respectively corresponding first equidistant time discrete function and second equidistant time discrete function comprises: Performing time point registration processing on the first warping function and the second warping function; Uniform equidistant time discretization processing is performed on the first warping function and the second warping function after the time point registration, respectively, to obtain the first equidistant time discretization function and the second equidistant time discretization function corresponding to each other.

4. The method according to any one of claims 1 to 3, characterized in that The determining, according to the first Fourier function and the second Fourier function, a difference amplitude value between the first discrete data set and the second discrete data set comprises: Calculating a first output value corresponding to each first data point in the first discrete data set based on the first Fourier function; and Calculate, based on the second Fourier function, a second output value corresponding to each second data point in the second discrete data set; Based on the first output value corresponding to each of the first data points and the second output value corresponding to each of the second data points, the difference amplitude value between the first discrete data set and the second discrete data set is determined, wherein: AA is the difference amplitude value between the first discrete data set and the second discrete data set, k is the total number of the first data points in the first discrete data set or the total number of the second data points in the second discrete data set, wherein the number of the first data points is the same as the number of the second data points and they correspond one to one, n is the sequence number of the first data point and the second data point corresponding to the first data point, F calc is the first output value corresponding to the nth first data point in the first discrete data set, F exp is the second output value corresponding to the nth second data point in the second discrete data set.

5. The method according to claim 4, characterized in that The step of determining the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model based on the difference amplitude value comprises: Determining a preset amplitude value range to which the difference amplitude value belongs; The preset accuracy corresponding to the preset amplitude value range is set as the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model.

6. A device for checking the accuracy of output results of a nuclear reactor transient safety analysis model, characterized in that: include: A current operating condition acquisition unit, used to acquire current operating condition information of the nuclear reactor; A model output acquisition unit, configured to output a first discrete data set of safety performance evaluation indicators of the current operating condition through a nuclear reactor transient safety analysis model based on the current operating condition information; An experiment output acquisition unit, used to acquire a second discrete data set obtained by an experiment of a safety performance evaluation index of a current working condition; an equidistant time discrete unit, used to determine, based on the first discrete data set and the second discrete data set, respectively corresponding first equidistant time discrete function and second equidistant time discrete function, wherein the first equidistant time discrete function and the second equidistant time discrete function respectively reflect changes in the operating conditions corresponding to the first discrete data set and the operating conditions corresponding to the second discrete data set on the same uniform time axis; A Fourier transform unit, configured to determine a first Fourier function based on the first discrete data set and the first equidistant time discrete function, and to determine a second Fourier function based on the second discrete data set and the second equidistant time discrete function; a difference amplitude value calculation unit, configured to determine a difference amplitude value between the first discrete data set and the second discrete data set according to the first Fourier function and the second Fourier function; The accuracy determination unit is used to determine the accuracy of the first discrete data set output by the nuclear reactor transient safety analysis model based on the difference amplitude value.

7. The device according to claim 6, characterized in that The equidistant time discrete unit comprises: A first fitting unit, used for performing fitting processing on the first discrete data set to obtain a first fitting function; A first characteristic phenomenon extraction unit, used for determining, in the first fitting function, a third discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor; A first warping function constructing unit, configured to construct a first warping function based on the third discrete data set; A second fitting unit, used for performing fitting processing on the second discrete data set to obtain a second fitting function; a second characteristic phenomenon extraction unit, configured to determine, in the second fitting function, a fourth discrete data set corresponding to the occurrence time of the characteristic phenomenon of the nuclear reactor; A second warping function constructing unit, configured to construct a second warping function based on the fourth discrete data set; The first execution unit is configured to determine, based on the first warping function and the second warping function, a first equidistant time discrete function and a second equidistant time discrete function that respectively correspond to each other.

8. The device according to claim 7, characterized in that The first execution unit is specifically used to: perform time point alignment processing on the first warping function and the second warping function, and respectively perform uniform equidistant time discretization processing on the first warping function and the second warping function after time point alignment to obtain the first equidistant time discretization function and the second equidistant time discretization function corresponding to each other.

9. A computer device, characterized in that: include: at least one processor; and, a memory communicatively coupled to the at least one processor; The memory stores instructions executable by the at least one processor, and the instructions are configured to execute the method according to any one of claims 1 to 5.

10. A computer-readable storage medium, characterized in that: Computer executable instructions are stored, and the computer executable instructions are used to execute the method according to any one of claims 1 to 5.

Images