Calculation method of multi-nuclide gamma energy spectrum characteristic peak area

Abstract

The invention provides a multi-nuclide gamma energy spectrum characteristic peak area calculation method, which comprises the following steps: acquiring a multi-nuclide gamma energy spectrum, and determining a characteristic peak range of the multi-nuclide gamma energy spectrum; gaussian fitting processing is carried out on the determined characteristic peak range; determining the position and the area of a corresponding characteristic peak according to a Gaussian fitting processing result; and calculating the area of the characteristic peak according to the determined position and area of the characteristic peak. By applying the method, the area of the characteristic peak of the multi-nuclide gamma energy spectrum can be quickly and accurately calculated.

Description

Technical Field

[0001] This application relates to the technical field of gamma spectroscopy analysis, and particularly to a method for calculating the area of characteristic peaks of multi-nuclide gamma spectra. Background Art

[0002] In gamma spectroscopy analysis, the calculation of the area of characteristic peaks is an important step, usually used to determine the activity or concentration of radioactive isotopes. Commonly used methods for calculating the area of characteristic peaks in the prior art include the integration method, baseline correction method, peak fitting method, peak separation method, and region segmentation method, etc.

[0003] Among them, the integration method is the most direct calculation method. This method obtains the characteristic peak area by numerically integrating the area under the characteristic peak (for example, the trapezoidal method or Simpson's method, etc.). This method is simple and easy to use, but its calculation accuracy may be affected when the background noise is large or the data resolution is not high. The baseline correction method can obtain a more accurate peak area by selecting a suitable baseline area for correction during the measurement process. This method can effectively reduce the influence of background noise on the calculation of the characteristic peak area. However, certain experience may be required to select the appropriate baseline position and method. The peak fitting method fits the characteristic peak with a mathematical model (for example, Gaussian or Lorentz function, etc.), and calculates the characteristic peak area through the fitting parameters. This method has stronger adaptability when dealing with data with high noise or irregular peak shapes, but certain skills in model selection and parameter adjustment are required. The peak separation method is suitable for analyzing overlapping peaks. It separates multiple overlapping characteristic peaks through a computer algorithm and calculates their areas respectively. This method can improve the interpretation accuracy of multi-peak spectrum analysis, especially in complex samples, but it often requires strong computing power and appropriate algorithm support during implementation. The region segmentation method divides the entire energy spectrum into multiple regions and only calculates the area of the region occupied by the characteristic peak, which can effectively avoid the interference of background noise. This method can simplify data analysis, but strong professional knowledge may be required when selecting regions. Summary of the Invention

[0004] In view of this, this application provides a method for calculating the area of characteristic peaks of multi-nuclide gamma spectra, so that the area of characteristic peaks of multi-nuclide gamma spectra can be calculated quickly and accurately.

[0005] The technical solution of this application is specifically implemented as follows:

[0006] A method for calculating the area of characteristic peaks of multi-nuclide gamma spectra, the method includes:

[0007] Obtain a multi-nuclide gamma spectrum and determine the range of characteristic peaks of the multi-nuclide gamma spectrum;

[0008] Perform Gaussian fitting processing on the determined range of characteristic peaks;

[0009] Determine the positions and regions of the corresponding characteristic peaks according to the results of Gaussian fitting processing;

[0010] Calculate the areas of the characteristic peaks based on the determined positions and regions of the characteristic peaks.

[0011] Preferably, determine the characteristic peak range of the multi-nuclide γ energy spectrum according to the empirical relationship between the target energy and the energy resolution of the scintillation detector for measuring the multi-nuclide γ energy spectrum.

[0012] Preferably, extract the relevant parameters of the characteristic peaks from the determined characteristic peak range through the Gaussian fitting processing.

[0013] Preferably, extract the relevant parameters of multiple characteristic peaks through least squares fitting.

[0014] Preferably, when performing Gaussian fitting processing, the expression of the Gaussian function is:

[0015]

[0016] Wherein, A is the maximum height of the characteristic peak, μ is the central position of the characteristic peak, and σ is the standard deviation, which reflects the width of the peak.

[0017] Preferably, determine the positions and regions of the corresponding characteristic peaks according to the central position and standard deviation of the characteristic peaks.

[0018] Preferably, calculate the area of the characteristic peak by performing numerical integration on the fitting curve.

[0019] Preferably, calculate the area of the characteristic peak through the following formula:

[0020]

[0021] Wherein, S is the area of the characteristic peak, A is the maximum height of the characteristic peak, and σ is the standard deviation.

[0022] As can be seen above, in the method for calculating the area of the characteristic peak of the multi-nuclide γ energy spectrum in this application, since the characteristic peak range of the multi-nuclide γ energy spectrum is first determined and then Gaussian fitting processing is performed on the determined characteristic peak range, the positions and regions of the corresponding characteristic peaks can be determined according to the results of the Gaussian fitting processing, and then the area of the characteristic peak can be calculated based on the determined positions and regions of the characteristic peaks, so that the area of the characteristic peak of the multi-nuclide γ energy spectrum can be calculated quickly and accurately. Brief Description of the Drawings

[0023] Figure 1 It is a schematic flow chart of the method for calculating the area of the characteristic peak of the multi-nuclide γ energy spectrum in the specific embodiment of this application. Detailed Description of the Invention

[0024] To make the technical solutions and advantages of this application clearer and more understandable, the following further elaborates on this application in conjunction with the accompanying drawings and specific embodiments.

[0025] Figure 1 It is a schematic flowchart of the calculation method for the characteristic peak area of the multi - nuclide γ energy spectrum in the embodiments of this application. As Figure 1 shown, the calculation method for the characteristic peak area of the multi - nuclide γ energy spectrum in the embodiments of this application includes the following steps:

[0026] Step 101, obtain the multi - nuclide γ energy spectrum and determine the characteristic peak range of the multi - nuclide γ energy spectrum.

[0027] In the technical solution of this application, since the characteristic peak usually represents the presence of certain specific substances or characteristics in the sample and reflects the presence of certain specific nuclides in the sample, the characteristic peak range of the multi - nuclide γ energy spectrum will be determined first in this step. This step is the starting point of the entire process and is also the basis for ensuring the accuracy of the final result.

[0028] For example, as an example, in a specific embodiment of this application, the characteristic peak range of the multi - nuclide γ energy spectrum can be determined according to the empirical relationship [E t - E t ·η, E t + E t ·η] between the target energy E t and the energy resolution η of the scintillation detector for measuring the multi - nuclide γ energy spectrum.

[0029] Among them, the energy resolution η reflects the ability of the detector to distinguish different energy peaks and is one of the important factors affecting the correctness of spectral analysis. Therefore, by combining the energy resolution η, the boundaries of the characteristic peaks of the multi - nuclide γ energy spectrum can be identified and defined, ensuring that the analysis focuses on the real signal while reducing the interference of background noise, and thus effectively improving the accuracy and reliability of the data. Especially when dealing with multi - nuclide samples, accurately identifying and separating different characteristic peaks is even more crucial because it will directly affect the subsequent quantitative analysis.

[0030] Step 102, perform Gaussian fitting processing on the determined characteristic peak range.

[0031] In the technical solution of this application, after determining the characteristic peak range of the multi - nuclide γ energy spectrum, Gaussian fitting processing can be performed on the determined characteristic peak range.

[0032] The above - mentioned Gaussian fitting processing is a process of fitting the signal within the determined characteristic peak range. In the γ energy spectrum, the characteristic peak can usually be well described by the Gaussian function.

[0033] For example, as an illustration, in a specific embodiment of the present application, relevant parameters of the characteristic peak (such as parameters like the height, position, and full width at half maximum of the characteristic peak) can be extracted from the determined characteristic peak range through Gaussian fitting processing.

[0034] For another example, as an illustration, in a specific embodiment of the present application, relevant parameters of multiple characteristic peaks (such as parameters like the height, position, and full width at half maximum of the characteristic peak) can be extracted through least squares fitting.

[0035] Based on the above parameters, the shape and intensity of the relevant characteristic peak can be understood, providing an analysis basis for quantitative and qualitative analysis, and non - peak situations can also be fed back.

[0036] For another example, as an illustration, in a specific embodiment of the present application, when performing Gaussian fitting processing, the expression of the Gaussian function is: where A is the maximum height of the characteristic peak, μ is the center position of the characteristic peak, and σ is the standard deviation, which reflects the width of the peak.

[0037] Through the above - mentioned Gaussian fitting processing, the difference between the theoretical curve and the experimental data can be minimized, ensuring that the fitting model can truly reflect the physical characteristics of the peak. The result of the Gaussian fitting processing can not only help the analyst distinguish accurate characteristic signals from complex backgrounds but also enhance the quantitative analysis of γ - spectrum information.

[0038] Step 103: Determine the position and region of the corresponding characteristic peak according to the result of the Gaussian fitting processing.

[0039] In the technical solution of the present application, after performing Gaussian fitting processing, the position and region of the corresponding characteristic peak can be determined according to the result of the Gaussian fitting processing.

[0040] For example, as an illustration, in a specific embodiment of the present application, the position and region of the corresponding characteristic peak can be determined according to the center position x0 and the standard deviation σ of the characteristic peak.

[0041] Step 104: Calculate the area of the characteristic peak according to the determined position and region of the characteristic peak.

[0042] In the technical solution of the present application, after determining the position and region of the characteristic peak, the area of the characteristic peak can be calculated according to the determined position and region of the characteristic peak.

[0043] This step is an important step for quantifying the intensity of the characteristic peak and can reflect the concentration or abundance of nuclides in the sample.

[0044] For example, as an illustration, in a specific embodiment of the present application, the area of the characteristic peak can be calculated by performing numerical integration on the fitting curve.

[0045] For example, as an illustration, in a specific embodiment of the present application, the integrated area of the characteristic peak can be determined by the following empirical formula:

[0046]

[0047] where FWHM is the integrated area of the characteristic peak and σ is the standard deviation;

[0048] The range of the integrated area of the obtained characteristic peak is:

[0049] The area of the characteristic peak is calculated by the following formula:

[0050]

[0051] where S is the area of the characteristic peak, A is the maximum height of the characteristic peak, and σ is the standard deviation.

[0052] For example, as an illustration, in a specific embodiment of the present application, the definite integral of equation (1) over the range of the integrated area of the characteristic peak can be calculated to obtain the area S:

[0053]

[0054] where the last term in the above formula is the error function, and its value can be obtained by looking up a table or calculated by substituting into the corresponding approximate calculation formula.

[0055] Therefore, through the above steps 101 to 104, the area of the characteristic peak of the multi-nuclide γ energy spectrum can be obtained.

[0056] The area of the characteristic peak is directly related to the concentration or abundance of the target substance (e.g., nuclide) in the sample. Therefore, after calculating the area of the characteristic peak, accurate quantitative analysis can be carried out, providing basic support for subsequent research and applications. For example, in technical fields such as radiation protection, nuclear energy development, and environmental monitoring, the area of the characteristic peak is an indispensable and important data.

[0057] In summary, in the technical solution of the present application, since the range of the characteristic peak of the multi-nuclide γ energy spectrum is first determined and then Gaussian fitting processing is performed on the determined range of the characteristic peak, the position and region of the corresponding characteristic peak can be determined according to the result of the Gaussian fitting processing. Subsequently, according to the position and region of the determined characteristic peak, the area of the characteristic peak can be calculated, so that the area of the characteristic peak of the multi-nuclide γ energy spectrum can be calculated quickly and accurately.

[0058] Compared with the prior art, the technical solution in the present application has a faster processing speed, and there is no need for manual discrimination of the detector performance, and it can automatically determine whether a characteristic peak is detected within the data range. In addition, the calculation area result in the present application has better effects in subsequent processing.

[0059] The above are only the preferred embodiments of the present application and are not intended to limit 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 protection of the present application.

Claims

1. A method for calculating the area of a characteristic peak in the multi - nuclide γ energy spectrum, characterized in that, The method includes: Obtaining a multi - nuclide γ energy spectrum and determining the characteristic peak range of the multi - nuclide γ energy spectrum; Performing Gaussian fitting processing on the determined characteristic peak range; Determining the position and region of the corresponding characteristic peak according to the result of the Gaussian fitting processing; Calculating the area of the characteristic peak according to the determined position and region of the characteristic peak.

2. The method according to claim 1, wherein: The characteristic peak range of the multi - nuclide γ energy spectrum is determined according to the empirical relationship between the target energy and the energy resolution of the scintillation detector for measuring the multi - nuclide γ energy spectrum.

3. The method according to claim 1, wherein: Through the Gaussian fitting processing, the relevant parameters of the characteristic peak are extracted from the determined characteristic peak range.

4. The method according to claim 3, wherein: The relevant parameters of multiple characteristic peaks are extracted by least - squares fitting.

5. The method according to claim 4, characterized in that When performing Gaussian fitting processing, the expression of the Gaussian function is: where A is the maximum height of the characteristic peak, μ is the center position of the characteristic peak, and σ is the standard deviation, which reflects the width of the peak.

6. The method according to claim 5, wherein: The position and region of the corresponding characteristic peak are determined according to the center position and standard deviation of the characteristic peak.

7. The method according to claim 6, wherein: The area of the characteristic peak is calculated by performing numerical integration on the fitting curve.

8. The method according to claim 7, wherein The area of the characteristic peak is calculated by the following formula: where S is the area of the characteristic peak, A is the maximum height of the characteristic peak, and σ is the standard deviation.

Images