Neutron energy spectrum unfolding method, device, equipment, medium and product

Abstract

The invention provides a spectrum unfolding method, device and equipment of a neutron spectrum, a medium and a product, and relates to the field of reactors. The method comprises the following steps: acquiring a theoretical neutron energy spectrum of a neutron source at a target position after passing through M media; obtaining neutron energy spectrum response functions corresponding to the N detectors respectively; extracting N feature basis functions based on the theoretical neutron spectrum; according to the reading of each detector, the neutron spectrum response function corresponding to each detector and the N feature basis functions, determining an expansion coefficient corresponding to each feature basis function; and determining the sum value of the products of each characteristic basis function and the corresponding expansion coefficient as a target neutron energy spectrum at the target position under the action of the neutron source. Thus, when the feature basis functions are extracted based on the theoretical neutron energy spectrum, the number of the feature basis functions is equal to the number of the detectors, the unique solution of the neutron energy spectrum can be obtained through solving, the underdetermined problem is avoided, and the accuracy of the spectrum obtaining result is improved.

Description

Technical Field

[0001] The present invention relates to the field of reactors, and more particularly to a method, apparatus, device, medium, and product for resolving neutron energy spectra. Background Art

[0002] The neutron energy spectrum, which is the distribution of neutron flux density with energy, is a key parameter in reactor core design. Therefore, it usually needs to be verified through the measurement results of reactor physics experiments. In general energy spectrum measurement methods, the neutron energy spectrum cannot be directly obtained and needs to be resolved from the measured values by combining the detection response function of the detector system.

[0003] Existing spectrum resolution methods such as the generalized least squares method are prone to overdetermined or underdetermined problems during spectrum resolution. On the one hand, the generalized least squares method has a low adaptability to matrix ill-conditioning problems, and the spectrum resolution results are prone to show a certain degree of oscillation. On the other hand, the accuracy of the solution of the generalized least squares method depends greatly on the initial value of the iteration. Therefore, the spectrum resolution results of existing spectrum resolution methods are not stable and accurate enough, and their applicability to problems needs to be expanded. Summary of the Invention

[0004] The purpose of the embodiments of the present application is to provide a method, apparatus, device, medium, and product for resolving neutron energy spectra, so that the spectrum resolution results are unique and the accuracy is improved.

[0005] In a first aspect, the embodiments of the present application provide a method for resolving neutron energy spectra, which is applied to a spectrum resolution device. The spectrum resolution device includes N detectors, and each detector is used to measure the neutron energy spectrum at a target position after a neutron source passes through M media; both N and M are positive integers;

[0006] The spectrum resolution method includes:

[0007] Obtain the theoretical neutron energy spectrum at the target position after the neutron source passes through M media;

[0008] Obtain the neutron energy spectrum response functions corresponding to the N detectors respectively;

[0009] Extract N characteristic basis functions based on the theoretical neutron energy spectrum;

[0010] Determine the expansion coefficients corresponding to the respective characteristic basis functions according to the readings of the respective detectors, the neutron energy spectrum response functions corresponding to the respective detectors, and the N characteristic basis functions;

[0011] Determine the sum of the products of each characteristic basis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source.

[0012] In some embodiments, obtaining the theoretical neutron energy spectrum at the target position after the neutron source passes through M media includes:

[0013] Obtain a first theoretical neutron energy spectrum and a second theoretical neutron energy spectrum; the first theoretical neutron energy spectrum is the theoretical neutron energy spectrum at a target position after the neutron source passes through M media under an optimal working condition, and the second theoretical neutron energy spectrum includes I theoretical neutron energy spectra at the target position after the neutron source passes through M media under I sets of perturbed working conditions, where I is a positive integer.

[0014] In some embodiments, obtaining the first theoretical neutron energy spectrum includes:

[0015] Obtain the optimal values of the parameters of the M media;

[0016] Determine the first theoretical neutron energy spectrum according to the optimal values.

[0017] In some embodiments, obtaining the second theoretical neutron energy spectrum includes:

[0018] Obtain the optimal values and error amounts of the parameters of the M media;

[0019] For each parameter in the M media, determine the upper limit value of the parameter by summing the optimal value of the parameter and the error amount, and determine the lower limit value of the parameter by subtracting the error amount of the parameter from the optimal value of the parameter;

[0020] Determine I parameter combinations corresponding to the M media; each parameter combination in the I parameter combinations includes the target values of the parameters in the M media; the target value of each parameter is any value between the upper limit value and the lower limit value of the parameter; at least one of the target values of the parameters in each parameter combination in the I parameter combinations is different;

[0021] Determine I second theoretical neutron energy spectra according to the I parameter combinations.

[0022] In some embodiments, extracting N characteristic basis functions based on the theoretical neutron energy spectrum includes:

[0023] Construct a matrix of the theoretical energy spectrum based on the theoretical neutron energy spectrum;

[0024] Decompose the matrix to obtain N characteristic basis functions.

[0025] In some embodiments, determining the expansion coefficients corresponding to the characteristic basis functions according to the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and the N characteristic basis functions includes:

[0026] Substitute the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and the N characteristic basis functions into the target formula to obtain the expansion coefficients corresponding to the characteristic basis functions;

[0027] Wherein, the target formula is:

[0028]

[0029] Among them, C n represents the reading of the nth detector, and R n (E g ) represents the neutron energy spectrum response function corresponding to the nth detector, and φ k (E g ) represents the kth eigenbasis function, and a k represents the expansion coefficient corresponding to the kth eigenbasis function. N represents the number of detectors, and G represents the number of energy groups of the neutron energy spectrum response function; the value range of n is from 1 to N, the value range of k is from 1 to N, and the value range of g is from 1 to G.

[0030] In a second aspect, an embodiment of the present application provides a spectrum unfolding device for neutron energy spectrum, which is applied to a spectrum unfolding device. The spectrum unfolding device includes N detectors, and each detector is used to measure the neutron energy spectrum at a target position after a neutron source passes through M media; both N and M are positive integers;

[0031] The spectrum unfolding device includes:

[0032] A first acquisition module, configured to acquire the theoretical neutron energy spectrum at the target position after the neutron source passes through M media;

[0033] A second acquisition module, configured to acquire the neutron energy spectrum response functions corresponding to the N detectors respectively;

[0034] An extraction module, configured to extract N eigenbasis functions based on the theoretical neutron energy spectrum;

[0035] A first determination module, configured to determine the expansion coefficients corresponding to the eigenbasis functions according to the readings of the detectors, the neutron energy spectrum response functions corresponding to the detectors, and the N eigenbasis functions;

[0036] A second determination module, configured to determine the sum of the products of each eigenbasis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source.

[0037] In a third aspect, an embodiment of the present application provides a spectrum unfolding device, including:

[0038] A memory, configured to store instructions; and

[0039] A processor, configured to call the instructions from the memory and be able to implement the spectrum unfolding method for neutron energy spectrum provided in the first aspect of the embodiments of the present application when executing the instructions.

[0040] In a fourth aspect, an embodiment of the present application provides a machine-readable storage medium, on which instructions are stored, and the instructions are used to enable the machine to execute the spectrum unfolding method for neutron energy spectrum according to the above.

[0041] Fifth aspect, an embodiment of the present application provides a computer program product. When the instructions in the computer program product are executed by a processor of a deconvolution device, the deconvolution device is caused to execute the deconvolution method of neutron energy spectrum as described above.

[0042] In an embodiment of the present application, a theoretical neutron energy spectrum at a target position after a neutron source passes through M media is obtained; neutron energy spectrum response functions corresponding to N detectors are obtained; N characteristic basis functions are extracted based on the theoretical neutron energy spectrum; according to the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and the N characteristic basis functions, expansion coefficients corresponding to each characteristic basis function are determined; the sum of the products of each characteristic basis function and its corresponding expansion coefficient is determined as the target neutron energy spectrum at the target position under the action of the neutron source. In this way, when extracting the characteristic basis functions based on the theoretical neutron energy spectrum, the number of characteristic basis functions is made equal to the number of detectors, and a unique solution of the neutron energy spectrum can be obtained by solving, avoiding the occurrence of underdetermined problems and improving the accuracy of the deconvolution result. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] Figure 1 is one of the flow schematic diagrams of the deconvolution method of neutron energy spectrum provided by an embodiment of the present application;

[0044] Figure 2 is a schematic diagram of obtaining a neutron energy spectrum under a perturbed state provided by an embodiment of the present application;

[0045] Figure 3 is another flow schematic diagram of the deconvolution method of neutron energy spectrum provided by an embodiment of the present application;

[0046] Figure 4 is a schematic structural diagram of the deconvolution device of neutron energy spectrum provided by an embodiment of the present application;

[0047] Figure 5 is a schematic structural diagram of the deconvolution device provided by an embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0048] The technical solutions in the embodiments of the present application will be clearly 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 in the present application belong to the scope of protection of the present application.

[0049] The terms "first", "second", etc. in the description and claims of this application are used to distinguish similar objects, rather than to describe a specific order or sequence. It should be understood that such terms can be interchanged under appropriate circumstances, so that the embodiments of this application can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first", "second", etc. are generally of the same type, and the number of objects is not limited. For example, the first object can be one or more. In addition, "and / or" in the description and claims means at least one of the connected objects, and the character " / ", generally represents an "or" relationship between the associated objects before and after.

[0050] The following will, with reference to the accompanying drawings, through specific embodiments and their application scenarios, provide a detailed description of the neutron energy spectrum deconvolution method, device, equipment, medium and product provided by the embodiments of this application.

[0051] Please refer to Figure 1 , which is one of the flow diagrams of the neutron energy spectrum deconvolution method provided by the embodiments of this application. This method is applied to a deconvolution device. Please refer to Figure 2 , which is a schematic diagram of obtaining the neutron energy spectrum under a perturbed state provided by the embodiments of this application. As Figure 2 shown, the deconvolution device includes N detectors, and each detector is used to measure the neutron energy spectrum at the target position after the neutron source passes through M media; where N and M are both positive integers. The detector measures the neutron energy spectrum at the target position after the neutron source passes through M media.

[0052] In the embodiments of this application, the number of detectors can be multiple.

[0053] As Figure 1 shown, the deconvolution method includes the following steps S100 to S500.

[0054] Step S100: Obtain the theoretical neutron energy spectrum at the target position after the neutron source passes through M media.

[0055] In the embodiments of this application, the neutron source can be understood as a device that generates neutrons, including a dedicated neutron generation device or neutrons spontaneously generated by fuel. The medium can be understood as the interaction medium between neutrons and the substances or materials they pass through. The medium can include, but is not limited to, gases, liquids, solid substances, and other materials, etc., aiming to affect and change the energy distribution, propagation mode, and interaction with the detector of neutrons, so as to affect the neutron energy spectrum at the target position. The neutron energy spectrum can be understood as the energy distribution of neutron flux density. The target position can be understood as the position where the neutron energy spectrum to be measured is located, and can be specifically set according to the actual situation.

[0056] In the embodiments of the present application, there may be one neutron source, or a combination of multiple neutron sources. The neutron source may not overlap with the position of the medium, or there may be an overlap.

[0057] In the embodiments of the present application, there may be M media between the detector and the neutron source, and the detector may be used to measure the neutron energy spectrum at the target position after the neutron source passes through the M media. The number of detectors may be N, where both N and M are positive integers.

[0058] Step S200: Obtain the neutron energy spectrum response functions corresponding to the N detectors respectively.

[0059] In the embodiments of the present application, for the N detectors, the neutron energy spectrum response function of each detector is R n (E g ), n = 1, 2,..., N, g = 1, 2,..., G, where G is the number of energy groups of the response function, and the typical value of the number of energy groups of the response function is several hundred groups. The neutron energy spectrum response function can characterize the conversion relationship between the neutron energy spectrum and the detector signal. E can be understood as the neutron energy, and E g can represent the energy value of the neutron corresponding to the number of energy groups of the response function.

[0060] In the embodiments of the present application, there are no restrictions on the type of the detector and the type of the neutron energy spectrum to be measured.

[0061] Step S300: Extract N characteristic basis functions based on the theoretical neutron energy spectrum.

[0062] In the embodiments of the present application, the theoretical neutron energy spectrum may include one type of neutron energy spectrum, or may include two types of neutron energy spectra. In one example, the theoretical neutron energy spectrum may include a first theoretical neutron energy spectrum and a second theoretical neutron energy spectrum; the first theoretical neutron energy spectrum φ 0 (E g ) can be understood as the theoretical neutron energy spectrum at the target position after the neutron source passes through the M media under the optimal working condition, and the parameters of the M media under the optimal working condition are all optimal values. The parameters of the medium may include, but are not limited to, the density, size, and temperature of the medium, etc. The optimal values of the parameters of the M media are preset values. Specifically, they can be determined according to parameters such as the energy response of the detector and the arrangement position of the detector.

[0063] In the embodiments of the present application, since there will be certain deviations between the true values and the optimal values of the parameters such as the density, size, and temperature of each medium, and these deviations may be caused by manufacturing processes or measurements, etc. Therefore, the neutron energy spectra corresponding to multiple working conditions can be obtained by combining the optimal values of each parameter with the error amounts of each parameter.

[0064] In the embodiments of the present application, the second theoretical neutron energy spectrum can be understood as including I theoretical neutron energy spectra of M media under I sets of perturbation conditions of the neutron source, where I is a positive integer. The parameter values of M parameters under the perturbation conditions are obtained by adding or subtracting an error amount based on the optimal values. The second theoretical neutron energy spectrum can be expressed as φ i (E g ), i = 1, 2,..., I, where I can be understood as the number of perturbation conditions, and I needs to satisfy being greater than or equal to the number N of detectors.

[0065] In one example, for the density parameter ρ 1 of the first medium, its error amount is Δρ 1 . Then the actual range of the density parameter of the first medium may cover from ρ 1 - Δρ 1 to ρ 1 + Δρ 1 .

[0066] Step S400: Determine the expansion coefficients corresponding to each eigenfunction according to the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and N eigenfunctions.

[0067] In the embodiments of the present application, a corresponding table of the readings of the pre-constructed detectors, the neutron energy spectrum response functions corresponding to the detectors, the eigenfunctions, and the expansion coefficients corresponding to each eigenfunction can be obtained. And the expansion coefficients corresponding to each eigenfunction are found according to the table.

[0068] In the embodiments of the present application, a relational expression between the readings of the detectors, the neutron energy spectrum response functions corresponding to the detectors, the eigenfunctions, and the expansion coefficients corresponding to each eigenfunction can be constructed. The relationship between the readings of the detectors, the neutron energy spectrum response functions corresponding to the detectors, the eigenfunctions, and the expansion coefficients corresponding to each eigenfunction can be expressed as:

[0069]

[0070] Among them, C n represents the reading of the nth detector, R n (E g ) represents the neutron energy spectrum response function corresponding to the nth detector, represents the kth eigenfunction, a k represents the expansion coefficient corresponding to the kth eigenfunction, N represents the number of detectors, G represents the energy group number of the neutron energy spectrum response function; the value range of n is from 1 to N, the value range of k is from 1 to N, and the value range of g is from 1 to G.

[0071] In an embodiment of the present application, after reading the readings of each detector, the expansion coefficients corresponding to each eigenbasis function can be solved according to the neutron energy spectrum response function corresponding to each detector, N eigenbasis functions, and the relational expression.

[0072] Step S500: Determine the sum of the products of each eigenbasis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source.

[0073] In an embodiment of the present application, after determining each eigenbasis function, the sum of the products of each eigenbasis function and its corresponding expansion coefficient can be determined as the target neutron energy spectrum at the target position under the action of the neutron source. The target neutron energy spectrum can be expressed as:

[0074]

[0075] wherein, φ(E g ) represents the target neutron energy spectrum, represents the eigenbasis function, and a k represents the expansion coefficient corresponding to the eigenbasis function. k = 1, 2,..., N, that is, the number of eigenbasis functions satisfies the number equal to the number of detectors.

[0076] Through the above steps S100 - step S500, obtain the theoretical neutron energy spectrum at the target position after the neutron source passes through M media; obtain the neutron energy spectrum response functions corresponding to N detectors respectively; extract N eigenbasis functions based on the theoretical neutron energy spectrum; determine the expansion coefficients corresponding to each eigenbasis function according to the readings of each detector, the neutron energy spectrum response functions corresponding to each detector, and N eigenbasis functions; determine the sum of the products of each eigenbasis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source. In this way, when extracting the eigenbasis function based on the theoretical neutron energy spectrum, the number of eigenbasis functions is made equal to the number of detectors, and a unique solution of the neutron energy spectrum can be obtained through solution, avoiding the occurrence of the underdetermined problem and improving the accuracy of the spectrum analysis result.

[0077] In some embodiments, obtaining the theoretical neutron energy spectrum at the target position after the neutron source passes through M media includes:

[0078] Obtain a first theoretical neutron energy spectrum and a second theoretical neutron energy spectrum; the first theoretical neutron energy spectrum is the theoretical neutron energy spectrum at the target position after the neutron source passes through M media under the optimal working condition, and the second theoretical neutron energy spectrum includes I theoretical neutron energy spectra at the target position after the neutron source passes through M media under I sets of perturbed working conditions, where I is a positive integer.

[0079] Specifically, the theoretical neutron energy spectrum at the target position after the neutron source passes through M media can include a first theoretical neutron energy spectrum and a second theoretical neutron energy spectrum. The first theoretical neutron energy spectrum φ0 (E g ) can be understood as the theoretical neutron energy spectrum at the target position after the neutron source passes through M media under the optimal working condition, and the parameters of the M media under the optimal working condition are all optimal values. When determining the first theoretical neutron energy spectrum, particle transport can be obtained through the Monte Carlo method.

[0080] In the embodiment of the present application, the second theoretical neutron energy spectrum can be understood as including I theoretical neutron energy spectra at the target position after the neutron source passes through M media under I sets of perturbed working conditions, where I is a positive integer. The perturbed working condition can be understood as the parameter value obtained by adding or subtracting an error amount on the basis of the optimal value of the M parameters. The second theoretical neutron energy spectrum can be expressed as φ i (E g ), i = 1, 2,..., I, where I can be understood as the number of perturbed working conditions, and I needs to satisfy being greater than or equal to the number N of detectors.

[0081] In this embodiment, by obtaining the first theoretical neutron energy spectrum and the second theoretical neutron energy spectrum, both the theoretical neutron energy spectrum under the optimal working condition and the theoretical neutron energy spectrum under the influence of parameter tolerances can be considered, improving the accuracy of determining the target neutron energy spectrum.

[0082] In some embodiments, obtaining the first theoretical neutron energy spectrum includes:

[0083] Obtaining the optimal values of the parameters of the M media;

[0084] Determining the first theoretical neutron energy spectrum according to the optimal values.

[0085] Specifically, the first theoretical neutron energy spectrum φ 0 (E g ) can be understood as the theoretical neutron energy spectrum at the target position after the neutron source passes through M media under the optimal working condition. The parameters of the M media under the optimal working condition are all optimal values. First, obtain the optimal values of the parameters of the M media, and determine the first theoretical neutron energy spectrum according to the optimal values of the parameters of the M media. In one example, when determining the first theoretical neutron energy spectrum, particle transport can be obtained by using the Monte Carlo method. The Monte Carlo Method is a computational method based on random sampling and statistical simulation, and its core idea is to approximately solve mathematical problems or simulate complex phenomena through a large number of random experiments.

[0086] In this embodiment, by determining the first theoretical neutron energy spectrum, the theoretical neutron energy spectrum under the optimal working condition can be considered, improving the accuracy of subsequent calculation of determining the target neutron energy spectrum.

[0087] In some embodiments, obtaining the second theoretical neutron energy spectrum includes:

[0088] Obtain the optimal values and error amounts of the parameters of M media;

[0089] For each parameter among the M media, determine the upper limit value of the parameter by summing the optimal value and the error amount of the parameter, and determine the lower limit value of the parameter by taking the difference between the optimal value of the parameter and the error amount of the parameter;

[0090] Determine I parameter combinations corresponding to the M media; each parameter combination among the I parameter combinations includes the target values of the parameters of the M media; the target value of each parameter is any value between the upper limit value and the lower limit value of the parameter; at least one of the target values of the parameters in each parameter combination among the I parameter combinations is different;

[0091] Determine I second theoretical neutron energy spectra according to the I parameter combinations.

[0092] Specifically, the second theoretical neutron energy spectrum can be understood as including I theoretical neutron energy spectra at the target position after the neutron source passes through M media under I sets of perturbation conditions, where I is a positive integer. The perturbation conditions can be understood as the value ranges between parameters obtained by adding or subtracting the error amount on the basis of the optimal values of the M parameters. First, obtain the optimal values and error amounts of the parameters of the M media. For each parameter among the M media, determine the upper limit value of the parameter by summing the optimal value and the error amount of the parameter, and determine the lower limit value of the parameter by taking the difference between the optimal value of the parameter and the error amount of the parameter.

[0093] In the embodiments of the present application, after determining the optimal values, upper limit values, and lower limit values of the parameters, determine I parameter combinations corresponding to the M media according to the optimal values, upper limit values, and lower limit values of the parameters. Each parameter combination among the I parameter combinations includes the target values of the parameters of the M media; the target value of each parameter is any value between the upper limit value and the lower limit value of the parameter; at least one of the target values of the parameters in each parameter combination among the I parameter combinations is different.

[0094] In the embodiments of the present application, the target value of each parameter can be selected as the optimal value, upper limit value, or lower limit value of the parameter.

[0095] In one example, when the parameters of the medium include density ρ, size L, and temperature T, each parameter has its own error amount Δρ, ΔL, and ΔT. For density ρ, ρ + Δρ, ρ, and ρ - Δρ can be taken; for size L, L + ΔL, L, and L - ΔL can be taken; for temperature T, T + ΔT, T, and T - ΔT can be taken. When performing perturbations, they can be carried out sequentially or jointly, and the perturbation process needs to include the working conditions of single-parameter perturbations of all parameters. The perturbation amount can be selected according to the deviation between each parameter and the nominal value. Combining the above values gives various perturbation calculations. There can be at most 27 combinations, or any partial combinations can be selected from the 27 combinations.

[0096] In one example, taking the error amounts as Δρ1, ΔL1, and ΔT1 respectively, there can be at most 27 combinations:

[0097] i = 1: density is ρ1 + Δρ1, size is L1 + ΔL1, temperature is T1 + ΔT1;

[0098] i = 2: density is ρ1 + Δρ1, size is L1 + ΔL1, temperature is T1;

[0099] i = 3: density is ρ1 + Δρ1, size is L1 + ΔL1, temperature is T1 - ΔT1;

[0100] i = 4: density is ρ1 + Δρ1, size is L1, temperature is T1 + ΔT1;

[0101] i = 5: density is ρ1 + Δρ1, size is L1, temperature is T1;

[0102] i = 6: density is ρ1 + Δρ1, size is L1, temperature is T1 - ΔT1;

[0103] …

[0104] i = I: density is ρ1 - Δρ1, size is L1 - ΔL1, temperature is T1 - ΔT1;

[0105] At this time, I = 27.

[0106] In another example, some combinations can also be arbitrarily selected from the above 27 methods. For example:

[0107] i = 1: density is ρ1, size is L1, temperature is T1;

[0108] i = 2: density is ρ1 + Δρ1, size is L1, temperature is T1;

[0109] i = 3: density is ρ1 - Δρ1, size is L1, temperature is T1;

[0110] i = 4: density is ρ1, size is L1 + ΔL1, temperature is T1;

[0111] i = 5: density is ρ1, size is L1 - ΔL1, temperature is T1;

[0112] i = 6: density is ρ1, size is L1, temperature is T1 - ΔT1;

[0113] i = 7: density is ρ1, size is L1, temperature is T1 + ΔT1;

[0114] At this time, I = 7. When making the selection, it is necessary to ensure that the range covered by these state points is relatively large.

[0115] In this embodiment, by considering the possible error amounts of various parameters, the second theoretical neutron energy spectrum under multiple perturbation conditions is determined, which can improve the accuracy of subsequent calculation of the target neutron energy spectrum determination.

[0116] In some embodiments, N eigenbasis functions are extracted based on the theoretical neutron energy spectrum, including:

[0117] Construct a matrix of the theoretical energy spectrum based on the theoretical neutron energy spectrum;

[0118] Decompose the matrix to obtain N eigenbasis functions.

[0119] Specifically, after obtaining the theoretical neutron energy spectrum, a matrix φ of the theoretical neutron energy spectrum can be constructed based on the theoretical neutron energy spectrum G×(I+1) . And decompose the matrix to obtain N eigenbasis functions

[0120] In one example, constructing a matrix of the theoretical energy spectrum is to arrange the first theoretical neutron energy spectrum and the second theoretical neutron energy spectrum in matrix form, which can be specifically expressed as:

[0121]

[0122] In the embodiment of the present application, after obtaining the matrix of the theoretical energy spectrum, the matrix can be decomposed. When decomposing, the proper orthogonal decomposition method can be used for decomposition. The proper orthogonal decomposition is a process of decomposing a matrix into eigenvectors. The eigenvectors are the eigenbasis functions in the embodiment of the present application. For each energy spectrum, it can be obtained by multiplying and adding the eigenbasis functions and the corresponding expansion coefficients.

[0123] In this embodiment, by decomposing the matrix, the complex matrix can be decomposed into a set of independent basis vectors and corresponding coefficients, making the calculation of the target neutron energy spectrum more convenient.

[0124] In some embodiments, according to the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and N eigenbasis functions, the expansion coefficients corresponding to each eigenbasis function are determined, including:

[0125] Substitute the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and N eigenbasis functions into the target formula to obtain the expansion coefficients corresponding to each eigenbasis function;

[0126] Among them, the target formula is:

[0127]

[0128] Among them, C n represents the reading of the nth detector, and R n (Eg ) represents the neutron energy spectrum response function corresponding to the nth detector, represents the kth eigenbasis function, a k represents the expansion coefficient corresponding to the kth eigenbasis function, N represents the number of detectors, and G represents the number of energy groups of the neutron energy spectrum response function; the value range of n is from 1 to N, and the value range of k is from 1 to N.

[0129] Specifically, after determining the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and N eigenbasis functions, the above parameters can be substituted into the target formula to obtain the expansion coefficients corresponding to each eigenbasis function.

[0130] In the embodiment of the present application, the target formula can be:

[0131] Among them, C n represents the reading of the nth detector and can be directly read. R n (E g ) represents the neutron energy spectrum response function corresponding to the nth detector, represents the kth eigenbasis function, a k represents the expansion coefficient corresponding to the kth eigenbasis function, N represents the number of detectors, and G represents the number of energy groups of the neutron energy spectrum response function; the value range of n is from 1 to N, and the value range of k is from 1 to N. Even if the number of eigenbasis functions is equal to the number of detectors, a unique solution can be obtained.

[0132] In one example, in the case of including three detectors, three target formulas can be constructed:

[0133]

[0134] By solving the system of linear equations composed of the three target formulas, the expansion coefficients corresponding to each eigenbasis function can be determined.

[0135] In this embodiment, by constructing the relationship between the readings of each detector, the neutron energy spectrum response function corresponding to each detector, and N eigenbasis functions and solving the expansion coefficients, the accuracy of solving the expansion coefficients can be improved.

[0136] Please refer to Figure 3 , which is the second schematic flow chart of the method for solving the neutron energy spectrum provided by the embodiment of the present application. As Figure 3 shown, the method for solving the spectrum includes the following steps:

[0137] Step S301, obtaining the theoretically best-estimated neutron energy spectrum based on the analysis object.

[0138] Step S302: Based on the analysis object, perturbations are made to obtain a series of theoretical neutron energy spectra.

[0139] Step S303: Extract characteristic basis functions from the theoretical neutron energy spectra.

[0140] Step S304: Use the detector readings to determine the expansion coefficients of the characteristic basis functions for the best-estimated neutron energy spectrum.

[0141] Step S305: Obtain the final neutron energy spectrum.

[0142] In the embodiment of the present application, the analysis object refers to a device capable of providing a neutron energy spectrum. This device is arranged by various materials, and the optimal input parameters are the nominal densities, nominal sizes, etc. of these materials and their arrangements.

[0143] In the embodiment of the present application, for the specific implementation, reference can be made to the foregoing description, which will not be elaborated here.

[0144] Please refer to Figure 4 , which is a schematic structural diagram of a neutron energy spectrum deconvolution device provided by an embodiment of the present application. The neutron energy spectrum deconvolution device is applied to a deconvolution device, and the deconvolution device includes N detectors. Each detector is used to measure the neutron energy spectrum at the target position after the neutron source passes through M media; both N and M are positive integers. In the second aspect of the embodiment of the present application, a neutron energy spectrum deconvolution device 40 is provided. The device 40 includes:

[0145] A first acquisition module 41, configured to acquire the theoretical neutron energy spectrum at the target position after the neutron source passes through M media;

[0146] A second acquisition module 42, configured to acquire the neutron energy spectrum response functions corresponding to the N detectors respectively;

[0147] An extraction module 43, configured to extract N characteristic basis functions based on the theoretical neutron energy spectrum;

[0148] A first determination module 44, configured to determine the expansion coefficients corresponding to the respective characteristic basis functions according to the readings of the respective detectors, the neutron energy spectrum response functions corresponding to the respective detectors, and the N characteristic basis functions;

[0149] A second determination module 45, configured to determine the sum value of the products of each characteristic basis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source.

[0150] The neutron energy spectrum deconvolution device 40 provided in the second aspect of the embodiment of the present application can implement each process implemented in the above method embodiment and achieve the same beneficial effects. To avoid repetition, it will not be elaborated here.

[0151] Please refer to Figure 5, which is a schematic structural diagram of the spectral analysis device provided by the embodiments of the present application. The embodiments of the present application also provide a spectral analysis device 5000, including a processor 5100 and a memory 5200. The memory 5200 stores machine-executable instructions that can be executed by the processor 5100, and the processor 5100 can execute the machine-executable instructions to implement the above-mentioned spectral analysis method of neutron energy spectrum.

[0152] In some embodiments, the embodiments of the present application also provide a machine-readable storage medium, on which instructions are stored. When the instructions are executed by a processor, the processor implements the above-mentioned spectral analysis method of neutron energy spectrum.

[0153] In some embodiments, the embodiments of the present application also provide a computer program product, including a computer program, which implements the spectral analysis method of neutron energy spectrum according to the above embodiments when executed by a processor.

[0154] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0155] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, and the combination of the processes and / or blocks in the flowchart and / or block diagram can also be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for implementing the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks. These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured product including an instruction device, and the instruction device implements the functions in Figure 1 one process or multiple processes and / or blocks Figure 1The functions specified in one or more boxes. These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process. Thus, the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes and / or boxes. Figure 1 One or more processes and / or boxes Figure 1 The steps of the functions specified in one or more boxes.

[0156] In a typical configuration, a computing device includes one or more processors (CPUs), an input / output interface, a network interface, and memory.

[0157] The memory may include non-permanent memory in the form of computer-readable media, random access memory (RAM), and / or non-volatile memory such as read-only memory (ROM) or flash memory (flash RAM). The memory is an example of computer-readable media.

[0158] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can store information by any method or technology. The information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette tapes, magnetic disk storage or other magnetic storage devices, or any other non-transmission media that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transitory media such as modulated data signals and carrier waves.

[0159] It should also be noted that the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, commodity or device comprising a series of elements not only includes those elements, but also includes other elements not expressly listed, or also includes elements inherent to such process, method, commodity or device. Without further limitation, an element defined by the statement "comprising one..." does not exclude the presence of another identical element in the process, method, commodity or device comprising the element.

[0160] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, various changes and modifications can be made to the present application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present application shall be included within the scope of the claims of the present application.

[0161] In addition, any combination can be made among various different embodiments of the present invention as long as it does not violate the idea of the present invention, and it should also be regarded as the content disclosed by the present invention.

Claims

1. A method for interpreting a neutron spectrum, characterized in that: Applied to a spectrum analysis device, the spectrum analysis device comprises N detectors, each of which is used to measure the neutron energy spectrum of a neutron source at a target position after passing through M media; N and M are both positive integers; The spectrum decomposition method comprises: Obtaining a theoretical neutron energy spectrum of the neutron source at a target position after the neutron source passes through the M media; Obtaining neutron energy spectrum response functions corresponding to the N detectors respectively; Extracting N characteristic basis functions based on the theoretical neutron energy spectrum; Determine an expansion coefficient corresponding to each characteristic basis function according to the reading of each detector, the neutron energy spectrum response function corresponding to each detector, and the N characteristic basis functions; The sum of the products of each characteristic basis function and its corresponding expansion coefficient is determined as the target neutron energy spectrum at the target position under the action of the neutron source.

2. The spectrum interpretation method according to claim 1, characterized in that: The obtaining of a theoretical neutron energy spectrum of the neutron source at a target position after the neutron source passes through the M media includes: A first theoretical neutron energy spectrum and a second theoretical neutron energy spectrum are obtained; the first theoretical neutron energy spectrum is the theoretical neutron energy spectrum at the target position after the neutron source passes through the M media under optimal conditions, and the second theoretical neutron energy spectrum includes I theoretical neutron energy spectra at the target position after the neutron source passes through the M media under I groups of disturbance conditions, where I is a positive integer.

3. The spectrum interpretation method according to claim 2, characterized in that: Obtain a first theoretical neutron spectrum, including: Obtaining optimal values ​​of various parameters of the M media; The first theoretical neutron energy spectrum is determined according to the optimal value.

4. The spectrum interpretation method according to claim 2, characterized in that: Obtain the second theoretical neutron spectrum, including: Obtaining the optimal value and error amount of each parameter of the M media; For each parameter in the M media, the sum of the optimal value of the parameter and the error amount is determined as the upper limit of the parameter, and the difference between the optimal value of the parameter and the error amount of the parameter is determined as the lower limit of the parameter; Determine I parameter combinations corresponding to the M media; each parameter combination in the I parameter combinations includes a target value of each parameter in the M media; the target value of each parameter is any value of the parameter between the upper limit value and the lower limit value; each parameter combination in the I parameter combinations has at least one parameter with a different target value; According to the I parameter combinations, the I second theoretical neutron energy spectra are determined.

5. The spectrum interpretation method according to claim 1, characterized in that: The extracting N characteristic basis functions based on the theoretical neutron energy spectrum includes: constructing a matrix of theoretical energy spectrum based on the theoretical neutron energy spectrum; The matrix is ​​decomposed to obtain the N characteristic basis functions.

6. The spectrum interpretation method according to claim 1, characterized in that: Determining the expansion coefficient corresponding to each characteristic basis function according to the reading of each detector, the neutron energy spectrum response function corresponding to each detector, and the N characteristic basis functions includes: Substituting the readings of each detector, the neutron energy spectrum response function corresponding to each detector and the N characteristic basis functions into the target formula to obtain the expansion coefficient corresponding to each characteristic basis function; Wherein, the target formula is: Among them, C n Represents the reading of the nth detector, R n (E g ) represents the neutron energy spectrum response function corresponding to the nth detector, represents the kth characteristic basis function, a k represents the expansion coefficient corresponding to the kth characteristic basis function, N represents the number of detectors, and G represents the number of energy groups of the neutron energy spectrum response function; the value range of n is 1 to N, the value range of k is 1 to N, and the value range of g is 1 to G.

7. A neutron spectrum analysis device, characterized in that: Applied to a spectrum analysis device, the spectrum analysis device comprises N detectors, each of which is used to measure the neutron energy spectrum of a neutron source at a target position after passing through M media; N and M are both positive integers; The spectrum analysis device comprises: A first acquisition module is used to acquire a theoretical neutron energy spectrum of the neutron source at a target position after the neutron source passes through the M media; A second acquisition module is used to acquire neutron energy spectrum response functions corresponding to the N detectors respectively; An extraction module, used for extracting N characteristic basis functions based on the theoretical neutron energy spectrum; A first determination module is used to determine the expansion coefficient corresponding to each characteristic basis function according to the reading of each detector, the neutron energy spectrum response function corresponding to each detector and the N characteristic basis functions; The second determination module is used to determine the sum of the products of each characteristic basis function and its corresponding expansion coefficient as the target neutron energy spectrum at the target position under the action of the neutron source.

8. A spectrum analysis device, characterized in that: include: a memory configured to store instructions; and a processor configured to call the instructions from the memory and implement the neutron energy spectrum decomposition method according to any one of claims 1 to 6 when executing the instructions.

9. A machine-readable storage medium, characterized in that: The machine-readable storage medium stores instructions for causing a machine to execute the neutron energy spectrum decomposition method according to any one of claims 1 to 6.

10. A computer program product, characterized in that When the instructions in the computer program product are executed by a processor of a spectrum decomposition device, the spectrum decomposition device executes the neutron energy spectrum decomposition method according to any one of claims 1 to 6.

Images