PGNAA detection data comprehensive analysis method

Abstract

The application provides a PGNAA detection data comprehensive analysis method, which relates to the technical field of detection and comprises the following steps: acquiring the first element content, neutron absorption cross section and environmental temperature variation amount and humidity of a sample to be detected; adjusting the ray self-absorption correction coefficient of the first element content, adjusting the flux distribution based on the neutron absorption cross section, calibrating the peak position based on the environmental temperature variation amount, and correcting the background based on the environmental humidity; and through case-by-case processing, when the relevant parameters exceed the threshold value, the synergistic calculation and simulation are adopted, otherwise the standard value is used, so that the accurate comprehensive analysis of the PGNAA detection data is realized, and the accuracy and adaptability of the detection result are improved.

Description

Technical Field

[0001] This invention relates to the field of detection technology, and more specifically, to a comprehensive analysis method for PGNAA detection data. Background Technology

[0002] As industrial production demands increasingly higher precision and efficiency in material composition detection, PGNAA technology has been widely adopted due to its non-destructive testing capabilities. However, in actual testing, dynamic changes in factors such as the elemental content of the sample, neutron absorption characteristics, and ambient temperature and humidity can lead to deviations in X-ray self-absorption, flux distribution, peak position, and background, affecting the reliability of the test results.

[0003] Traditional PGNAA detection data analysis methods do not handle different cases of element content, neutron absorption cross section and environmental parameters separately. They often use uniform standard values ​​or simplified models, which cannot take into account the synergistic effect of multiple elements, complex flux distribution and the impact of dynamic environmental changes. This results in insufficient absorption correction for high-content elements and signal distortion for low-content elements, making it difficult to meet the detection accuracy requirements.

[0004] Therefore, it is necessary to design a comprehensive analysis method for PGNAA detection data to solve the problem that the existing technology does not differentiate between element content, neutron absorption cross section and environmental parameters, resulting in large deviations and insufficient accuracy in the detection results, and is unable to adapt to complex detection scenarios. Summary of the Invention

[0005] In view of this, the present invention proposes a comprehensive analysis method for PGNAA detection data, which aims to solve the problem that the existing technology does not process element content, neutron absorption cross section and environmental parameters separately, resulting in large deviations and insufficient accuracy in the detection results.

[0006] In one aspect, this invention proposes a comprehensive analysis method for PGNAA detection data, comprising:

[0007] S1, obtain the first element content, neutron absorption cross section, ambient temperature change, and ambient humidity of the sample being tested;

[0008] S2, adjust the X-ray self-absorption correction coefficient of the first element based on the content of the first element;

[0009] S3, adjust the flux distribution based on the neutron absorption cross section;

[0010] S4, calibrate the peak position based on the ambient temperature change;

[0011] S5, correct the background based on the ambient humidity.

[0012] Furthermore, when adjusting the X-ray self-absorption correction coefficient of the first element based on its content, the following steps are included:

[0013] When the content of the first element is greater than the threshold of the first element content, the X-ray self-absorption correction coefficient of the first element is determined by multi-element collaborative calculation.

[0014] When the content of the first element is less than or equal to the content threshold of the first element, the X-ray self-absorption correction coefficient of the first element is equal to the standard X-ray self-absorption correction coefficient of the first element.

[0015] Furthermore, when determining the X-ray self-absorption correction coefficient of the first element using a multi-element collaborative calculation method, the following steps are included:

[0016] Extract the absorption interaction elements in the sample that have absorption interaction with the first element;

[0017] Obtain the standard X-ray self-absorption correction factor, atomic number, content, and X-ray energy of the first element, and the standard X-ray self-absorption correction factor, atomic number, content, and X-ray energy of each of the absorption interaction elements;

[0018] The first correction factor is obtained by dividing the product of the standard X-ray self-absorption correction factor, atomic number, and content of the first element by the X-ray energy of the first element.

[0019] The second correction factor is obtained by dividing the product of the standard X-ray self-absorption correction factor, atomic number, and content of the absorbing interaction element by the X-ray energy of the absorbing interaction element.

[0020] The third correction coefficient is obtained by adding the second correction coefficient of each of the absorbing interaction elements;

[0021] The X-ray self-absorption correction coefficient of the first element is obtained by summing the third correction coefficient and the first correction coefficient.

[0022] Furthermore, adjusting the flux distribution based on the neutron absorption cross section includes:

[0023] When the neutron absorption cross section is greater than the neutron absorption cross section threshold, the flux distribution is determined using Monte Carlo simulation.

[0024] When the neutron absorption cross section is less than or equal to the neutron absorption cross section threshold, the flux distribution is equal to the standard flux distribution.

[0025] Furthermore, when determining the flux distribution using Monte Carlo simulation, the following steps are included:

[0026] Extract neutron interaction elements from the sample to be tested that interact with neutrons;

[0027] Obtain the neutron scattering cross section, neutron absorption cross section, and content of the neutron-interacting element;

[0028] The first distribution parameter is obtained by dividing the product of the neutron scattering cross section, neutron absorption cross section, and content of the neutron interacting element by a preset neutron energy constant.

[0029] The second distribution parameter is obtained by adding the first distribution parameters of each of the neutron interaction elements;

[0030] The flux distribution is obtained by substituting the second distribution parameter into the Monte Carlo simulation model.

[0031] Furthermore, when calibrating the peak position based on the aforementioned change in ambient temperature, the following steps are included:

[0032] When the change in ambient temperature is greater than the temperature change threshold, the peak position calibration value is determined using a temperature segment calibration method.

[0033] When the change in ambient temperature is less than or equal to the temperature change threshold, the peak position calibration value is equal to the standard peak position calibration value.

[0034] Furthermore, when determining the peak calibration value using the temperature segmentation calibration method, the following steps are included:

[0035] The ambient temperature is divided into several continuous temperature ranges;

[0036] Determine the temperature range to which the current ambient temperature belongs;

[0037] Obtain the peak temperature coefficient and standard peak value corresponding to the temperature range;

[0038] The peak position calibration value is obtained by multiplying the peak position temperature coefficient, the difference between the current ambient temperature and the standard temperature, and the standard peak position value.

[0039] Further, when correcting the background based on the ambient humidity, the following steps are included:

[0040] When the ambient humidity is greater than the humidity threshold, the background correction amount is determined using a dynamic humidity correction method.

[0041] When the ambient humidity is less than or equal to the humidity threshold, the background correction amount is equal to the standard background correction amount.

[0042] Furthermore, when determining the background correction amount using a dynamic humidity correction method, the following steps are included:

[0043] Obtain the initial humidity and current humidity of the sample being tested;

[0044] The change in humidity is obtained by calculating the difference between the current humidity and the initial humidity.

[0045] Obtain the preset humidity correction factor and standard background spectrum;

[0046] The background correction amount is obtained by multiplying the humidity change, the humidity correction factor, and the standard background spectrum.

[0047] Compared with the prior art, the beneficial effects of the present invention are as follows: the comprehensive analysis method of PGNAA detection data of the present invention determines the X-ray self-absorption correction coefficient according to the content of the first element, adjusts the flux distribution according to the neutron absorption cross section, calibrates the peak position in combination with the change of ambient temperature, and corrects the background based on the ambient humidity. It realizes targeted processing of different element contents, neutron absorption characteristics and environmental parameters, avoids the limitations of uniform standard values ​​or simplified models, and can comprehensively consider the synergistic effect of multiple elements and the influence of environmental factors, thereby improving the accuracy and adaptability of comprehensive analysis of PGNAA detection data and meeting the needs of complex detection scenarios. Attached Figure Description

[0048] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings:

[0049] Figure 1 A flowchart of the PGNAA detection data comprehensive analysis method provided in this embodiment of the invention. Detailed Implementation

[0050] Exemplary embodiments of the present invention will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided to enable a more thorough understanding of the present invention and to fully convey the scope of the invention to those skilled in the art. It should be noted that, without conflict, the embodiments and features described herein can be combined with each other. The present invention will now be described in detail with reference to the accompanying drawings and embodiments.

[0051] Reference Figure 1 As shown in some embodiments of this application, a comprehensive analysis method for PGNAA detection data includes:

[0052] S1, obtain the first element content, neutron absorption cross section, ambient temperature change, and ambient humidity of the sample being tested;

[0053] S2, adjust the X-ray self-absorption correction coefficient of the first element based on the content of the first element;

[0054] S3, adjust the flux distribution based on the neutron absorption cross section;

[0055] S4, calibrate the peak position based on the ambient temperature change;

[0056] S5, correct the background based on the ambient humidity.

[0057] Specifically, the first element refers to elements with a high atomic number (atomic number greater than 20) in the periodic table. These elements have a higher number of protons and neutrons in their nuclei, resulting in stronger absorption of gamma rays. In PGNAA detection, their content and distribution in the sample significantly affect the self-absorption effect of gamma rays, thus impacting the accuracy of the detection data. Common examples include iron (Fe, atomic number 26) and lead (Pb, atomic number 82). The content of each element in the sample is obtained through the following method: using PGNAA technology, the sample is irradiated with a neutron source, activating the elements and generating characteristic gamma rays. These characteristic gamma rays are then received by a detector, and the content is analyzed according to... The energy and intensity of characteristic gamma rays, combined with known elemental gamma ray information and relevant calibration curves, are used to calculate the content of each element. The neutron absorption cross section is obtained by analyzing the main elemental composition of the sample, using the known neutron absorption cross section data of each element, and combining the element's content ratio in the sample, through weighted calculation. The known neutron absorption cross section data of each element comes from experimentally determined standard values ​​in the nuclear physics database. The gamma ray self-absorption correction coefficient is a parameter used to correct the absorption effect of the sample itself on the gamma rays. This coefficient can adjust the deviation of the detection data caused by the absorption of the gamma rays by the sample, making the analysis results more consistent with the actual situation.

[0058] Understandably, by determining the X-ray self-absorption correction coefficient based on the content of the first element, adjusting the flux distribution according to the neutron absorption cross section, calibrating the peak position in conjunction with changes in ambient temperature, and correcting the background based on ambient humidity, targeted processing of different element contents, neutron absorption characteristics, and environmental parameters is achieved. This avoids the limitations of uniform standard values ​​or simplified models, and can comprehensively consider the synergistic effects of multiple elements and the influence of environmental factors, thereby improving the accuracy and adaptability of PGNAA detection data comprehensive analysis and meeting the needs of complex detection scenarios.

[0059] In some embodiments of this application, adjusting the X-ray self-absorption correction coefficient of the first element based on its content includes:

[0060] When the content of the first element is greater than the threshold of the first element content, the X-ray self-absorption correction coefficient of the first element is determined by multi-element collaborative calculation.

[0061] When the content of the first element is less than or equal to the content threshold of the first element, the X-ray self-absorption correction coefficient of the first element is equal to the standard X-ray self-absorption correction coefficient of the first element.

[0062] Specifically, the first element content threshold is the element content value determined by statistically analyzing the critical value at which the first element has a significant impact on X-ray self-absorption; the standard X-ray self-absorption correction coefficient is a benchmark coefficient predetermined for the first element or absorbing interacting elements under specific standard conditions (such as standard content, standard environment, etc.) to correct for the intensity attenuation of X-rays in the tested sample due to the self-absorption effect. It is usually obtained through experimental measurement, i.e., by performing PGNAA testing on standard samples with known composition and content, and then analyzing the attenuation law of X-rays in the standard samples using theoretical calculation models (such as Monte Carlo simulations) to determine the coefficient; alternatively, it can be obtained by consulting relevant standard databases or literature for verified correction coefficients of the same type of element under standard conditions.

[0063] Understandably, different methods for determining the X-ray self-absorption correction coefficient are set for the content of the first element. When the content exceeds the threshold, collaborative calculation is used; otherwise, the standard value is used. This avoids the limitations of a single processing method and improves the adaptability of the correction coefficient.

[0064] In some embodiments of this application, when determining the X-ray self-absorption correction coefficient of the first element using a multi-element collaborative calculation method, the following steps are included:

[0065] Extract the absorption interaction elements in the sample that have absorption interaction with the first element;

[0066] Obtain the standard X-ray self-absorption correction factor, atomic number, content, and X-ray energy of the first element, and the standard X-ray self-absorption correction factor, atomic number, content, and X-ray energy of each of the absorption interaction elements;

[0067] The first correction factor is obtained by dividing the product of the standard X-ray self-absorption correction factor, atomic number, and content of the first element by the X-ray energy of the first element.

[0068] The second correction factor is obtained by dividing the product of the standard X-ray self-absorption correction factor, atomic number, and content of the absorbing interaction element by the X-ray energy of the absorbing interaction element.

[0069] The third correction coefficient is obtained by adding the second correction coefficient of each of the absorbing interaction elements;

[0070] The X-ray self-absorption correction coefficient of the first element is obtained by summing the third correction coefficient and the first correction coefficient.

[0071] Specifically, by analyzing the differences in atomic number between each element in the tested sample and the first element, the overlap of the X-ray energy absorption cross section, and the correlation of element content, elements that significantly affect the X-ray self-absorption of the first element are screened out as absorption interaction elements that have absorption interactions with the first element. X-ray energy refers to the specific energy value of the γ-rays released by the first element and absorption interaction elements in the tested sample after neutron activation. Different elements will release characteristic γ-rays after activation, and their energies are unique and can be used to identify element types. X-ray energy is obtained in the following way: the detector is calibrated using a standard radioactive source with known energy (such as cobalt-60, cesium-137, etc.) to establish the correspondence between the detector output signal and the X-ray energy. Then, the γ-rays released by the tested sample after neutron activation are measured. The energy of the measured ray is determined by the detector output signal and the established correspondence.

[0072] Specifically, when the content of the first element is greater than the threshold value, its absorption interaction with other elements in the sample is significant. Therefore, it is not possible to rely solely on the standard value. After extracting the absorption interaction elements, based on the characteristics that atomic number reflects the difference in the element's ability to absorb radiation, the influence of content on the strength of absorption, and the correlation between radiation energy and the degree of absorption attenuation, the product of the standard radiation self-absorption correction coefficient, atomic number, and content of the first element and each absorption interaction element is divided by their respective radiation energy. Then, by summing and integrating the absorption effects of the first element itself and each absorption interaction element, the actual absorption situation under the synergistic effect of multiple elements can be comprehensively reflected, thereby obtaining the accurate radiation self-absorption correction coefficient of the first element.

[0073] Understandably, the specific steps for calculating the X-ray self-absorption correction coefficient through multi-element collaborative calculation have been refined. By extracting and quantifying interactive elements, the influence of multi-element absorption is comprehensively considered, making the correction of the high-content first element more accurate.

[0074] In some embodiments of this application, adjusting the flux distribution based on the neutron absorption cross section includes:

[0075] When the neutron absorption cross section is greater than the neutron absorption cross section threshold, the flux distribution is determined using Monte Carlo simulation.

[0076] When the neutron absorption cross section is less than or equal to the neutron absorption cross section threshold, the flux distribution is equal to the standard flux distribution.

[0077] Specifically, when the neutron absorption cross section is greater than the neutron absorption cross section threshold, the absorption of neutrons by the sample is stronger, and the scattering, attenuation, and other interaction processes of neutrons in the sample are more complex. At this time, methods based on simplified assumptions such as homogeneous medium are difficult to accurately reflect the actual propagation and distribution of neutrons. Monte Carlo simulation, on the other hand, simulates the trajectory of a large number of neutrons by random sampling, which can accurately characterize the complex microscopic processes of scattering and absorption between neutrons and matter, thus more reliably determining the flux distribution. The neutron absorption cross section threshold is a cross section value determined by statistically analyzing the critical value that has a significant impact on the flux distribution. The standard flux distribution is the neutron flux distribution data obtained by measuring the neutron flux of standard samples with known composition under preset standard experimental conditions. The standard flux distribution is obtained by measuring the neutron flux of standard samples without neutron interacting elements or with neutron interacting element content below a preset threshold under preset experimental conditions, recording the flux values ​​of neutrons at different positions in the standard samples, and then processing and averaging the data.

[0078] Understandably, the flux distribution is determined based on the neutron absorption cross section. Monte Carlo simulation is used when the threshold is exceeded, and the standard distribution is used otherwise. This approach takes into account both complex and simple scenarios and improves the rationality of the flux distribution.

[0079] In some embodiments of this application, determining the flux distribution using Monte Carlo simulation includes:

[0080] Extract neutron interaction elements from the sample to be tested that interact with neutrons;

[0081] Obtain the neutron scattering cross section, neutron absorption cross section, and content of the neutron-interacting element;

[0082] The first distribution parameter is obtained by dividing the product of the neutron scattering cross section, neutron absorption cross section, and content of the neutron interacting element by a preset neutron energy constant.

[0083] The second distribution parameter is obtained by adding the first distribution parameters of each of the neutron interaction elements;

[0084] The flux distribution is obtained by substituting the second distribution parameter into the Monte Carlo simulation model.

[0085] Specifically, the preset neutron energy constant is a constant related to neutron energy that is experimentally determined and fixed with reference to the average energy of the neutron source actually used. When extracting neutron-interacting elements that interact with neutrons in the sample, elements with atomic numbers greater than or equal to a preset threshold (e.g., Z≥10) are first screened based on the periodic table. The atomic nuclei of these elements have sufficient mass numbers to cause significant scattering or absorption of neutrons. Then, the sample is scanned with an energy-dispersive X-ray fluorescence spectrometer (EDXRF) to identify major and minor elements with a relative content greater than 0.1% as candidate interacting elements. When obtaining the neutron scattering and absorption cross sections of neutron-interacting elements, standard cross section data of each element in the thermal neutron energy region (0.0253 eV) were extracted from the ENDF / B-VIII.0 database published by the International Atomic Energy Agency (IAEA). For isotopes not included in the database, theoretical calculation models (such as optical model potentials) were used to extrapolate the data from experimental data. The elemental content was quantitatively calculated by measuring the characteristic gamma-ray peak area using neutron activation analysis (NAA) and combining it with the calibration curves of known standard materials. For light elements (such as H, C, and N), elastic recoil detection analysis (ERDA) was used for supplementary measurements. Finally, the content of each element was normalized to a mass percentage.

[0086] Specifically, when the neutron absorption cross section is greater than the threshold, the interaction between neutrons and the sample is strong and complex. After extracting the neutron interaction elements, their scattering cross section, absorption cross section, and content jointly determine the scattering and absorption behavior of neutrons. The first distribution parameter obtained by multiplying the product of the three and dividing by the neutron energy constant can quantify the influence of a single element on the neutron flux. The second distribution parameter obtained by summing them up can comprehensively reflect the overall effect of all interaction elements. Substituting them into the Monte Carlo simulation model can accurately simulate the trajectory and distribution law of neutrons under complex interactions, thereby obtaining an accurate flux distribution.

[0087] Understandably, the operation of Monte Carlo simulation in determining flux distribution has been clarified. By extracting neutron interaction elements and calculating distribution parameters, flux distribution calculation in complex scenarios becomes more targeted and accurate.

[0088] In some embodiments of this application, calibrating the peak position based on the ambient temperature change includes:

[0089] When the change in ambient temperature is greater than the temperature change threshold, the peak position calibration value is determined using a temperature segment calibration method.

[0090] When the change in ambient temperature is less than or equal to the temperature change threshold, the peak position calibration value is equal to the standard peak position calibration value.

[0091] Specifically, the temperature change threshold is the minimum temperature change required for segmented calibration, determined by statistically analyzing the impact of different temperature change amplitudes on peak position calibration accuracy. The standard peak position calibration value is a fixed peak position calibration value determined based on the standard temperature and initial measurement conditions when the ambient temperature change does not exceed the temperature change threshold. The standard peak position calibration value is determined by measuring the detector with a standard radiation source of known energy under preset standard temperature conditions, obtaining the characteristic peak position data corresponding to the standard radiation source, and using this characteristic peak position data as a reference.

[0092] Understandably, peak position is calibrated based on the amount of change in ambient temperature. When the threshold is exceeded, segmented calibration is used; otherwise, the standard value is used. This reduces the impact of temperature fluctuations on peak position and improves the reliability of peak position calibration.

[0093] In some embodiments of this application, determining the peak calibration value using a temperature segmentation calibration method includes:

[0094] The ambient temperature is divided into several continuous temperature ranges;

[0095] Determine the temperature range to which the current ambient temperature belongs;

[0096] Obtain the peak temperature coefficient and standard peak value corresponding to the temperature range;

[0097] The peak position calibration value is obtained by multiplying the peak position temperature coefficient, the difference between the current ambient temperature and the standard temperature, and the standard peak position value.

[0098] Specifically, the peak position temperature coefficient is determined within a specific temperature range by measuring the ratio of the peak position shift of a standard sample at different temperatures to the temperature change, and is used to quantify the influence of temperature on the peak position. The standard peak position value is a reference value of the peak position obtained by measuring a standard source ray of known energy under standard temperature conditions. The standard temperature is a constant temperature value set experimentally as a reference for peak position calibration, usually the ambient temperature during the calibration or initial measurement of the detection equipment. When obtaining the peak position temperature coefficient and standard peak position value corresponding to the temperature range, a standard sample covering each temperature range is first selected, and the characteristic peak position of the standard sample is measured at multiple temperature points within each temperature range. By calculating the difference between the peak position at different temperature points and the peak position at the standard temperature, and combining the temperature change, the peak position temperature coefficient of the range is obtained. At the same time, the characteristic peak position measured at the standard temperature is determined as the standard peak position value corresponding to the range.

[0099] Specifically, when the change in ambient temperature exceeds the temperature change threshold, the effect of temperature on the detector exhibits nonlinear characteristics. Dividing the ambient temperature into several continuous intervals makes the relationship between temperature and peak position shift within each interval approximately linear. By determining the interval to which the current temperature belongs, the peak position temperature coefficient (reflecting the relative shift of the peak position caused by a unit temperature change within the interval) and the standard peak position value corresponding to that interval are obtained. Then, the difference between the current ambient temperature and the standard temperature is combined, and the three are multiplied to quantify the peak position shift caused by the temperature change within that interval, thereby obtaining the peak position calibration value.

[0100] Understandably, the steps for temperature-segmented peak calibration have been refined. By dividing the range and quantifying the calibration values, the peak calibration is made more accurate when the temperature changes significantly, thus adapting to different temperature environments.

[0101] In some embodiments of this application, background correction based on the ambient humidity includes:

[0102] When the ambient humidity is greater than the humidity threshold, the background correction amount is determined using a dynamic humidity correction method.

[0103] When the ambient humidity is less than or equal to the humidity threshold, the background correction amount is equal to the standard background correction amount.

[0104] Specifically, the humidity threshold is a humidity value determined by statistically analyzing the critical value at which ambient humidity has a significant impact on background noise; the standard background correction is a background correction benchmark value obtained by measuring the environment without a test sample under preset humidity conditions; the process of obtaining the standard background correction is as follows: under preset humidity conditions, ensuring that there is no test sample in the detection environment, using the same PGNAA equipment and parameter settings as the actual detection, measuring the background rays in the environment to obtain background spectrum data, and then processing the background spectrum data by baseline subtraction, noise filtering, etc., to obtain a benchmark value that can be used for subsequent background correction, and this benchmark value is the standard background correction.

[0105] Understandably, the background correction is based on environmental humidity, with dynamic correction used when the threshold is exceeded and standard values ​​used otherwise. This effectively addresses the impact of humidity changes on the background and improves the applicability of background correction.

[0106] In some embodiments of this application, determining the background correction amount using a dynamic humidity correction method includes:

[0107] Obtain the initial humidity and current humidity of the sample being tested;

[0108] The change in humidity is obtained by calculating the difference between the current humidity and the initial humidity.

[0109] Obtain the preset humidity correction factor and standard background spectrum;

[0110] The background correction amount is obtained by multiplying the humidity change, the humidity correction factor, and the standard background spectrum.

[0111] Specifically, the initial humidity is the initial value of the ambient humidity recorded at the start of the detection; the preset humidity correction coefficient is a correction ratio coefficient corresponding to the humidity change obtained by experimentally measuring the background noise change pattern under different humidity conditions; the standard background spectrum is the background X-ray energy spectrum data obtained by measuring the background noise change corresponding to different humidity conditions under preset humidity conditions; the preset humidity correction coefficient is obtained by performing multiple detections on the sample under different humidity conditions, measuring the background noise change corresponding to different humidity conditions, and then calculating it based on the ratio of the background noise change to the humidity change; the standard background spectrum is obtained by detecting the environment without the sample under preset standard humidity conditions, collecting and recording the background X-ray spectrum at this time.

[0112] Specifically, when the ambient humidity exceeds the humidity threshold, humidity changes significantly affect the physical state of the sample (e.g., changes in water content may alter sample density and scattering characteristics) and the scattering background of the environment around the detector (e.g., water molecules enhance the scattering of X-rays). The difference between the initial humidity and the current humidity (the amount of humidity change) directly reflects the degree of this influence. The preset humidity correction coefficient is a proportional coefficient used to quantify the impact of humidity changes on the background, based on a large amount of experimental data. The standard background spectrum represents the background radiation level under the baseline conditions. Therefore, multiplying the three together can comprehensively quantify the impact of humidity changes on the background and obtain an accurate background correction amount.

[0113] Understandably, the specific operation of dynamic humidity correction for the background has been clarified. By calculating and quantifying the change in humidity, the background correction in high humidity environments is more accurate and humidity interference is reduced.

[0114] It should be noted that:

[0115] Numerous specific details are set forth in the specification provided herein. However, it will be understood that embodiments of this application may be practiced without these specific details. In some instances, well-known structures and techniques have not been shown in detail so as not to obscure the understanding of this specification.

[0116] Furthermore, those skilled in the art will understand that although some embodiments described herein include certain features included in other embodiments but not others, combinations of features from different embodiments are meant to be within the scope of this application and form different embodiments.

[0117] The above description is merely a preferred embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method of comprehensive analysis of PGNAA detection data, characterized in that, The method comprises the following steps: S1, obtaining the first element content, neutron absorption cross section, and environmental temperature variation of the measured sample; S2, adjusting the first element's ray self-absorption correction coefficient based on the first element content; S3, adjusting the flux distribution based on the neutron absorption cross section; S4, calibrating the peak position based on the environmental temperature variation; S5, correcting the background based on the environmental humidity; When adjusting the first element's ray self-absorption correction coefficient based on the first element content, the method comprises the following steps: When the first element content is greater than the first element content threshold, the first element's ray self-absorption correction coefficient is determined by using a multi-element collaborative calculation method; When the first element content is less than or equal to the first element content threshold, the first element's ray self-absorption correction coefficient is equal to the first element's standard ray self-absorption correction coefficient; When the first element's ray self-absorption correction coefficient is determined by using the multi-element collaborative calculation method, the method comprises the following steps: Absorption interaction elements that have absorption interaction with the first element in the measured sample are extracted; The first element's standard ray self-absorption correction coefficient, atomic number, content, and ray energy, and the standard ray self-absorption correction coefficient, atomic number, content, and ray energy of each absorption interaction element are obtained; The product of the first element's standard ray self-absorption correction coefficient, atomic number, and content is divided by the first element's ray energy to obtain a first correction coefficient; The product of the absorption interaction element's standard ray self-absorption correction coefficient, atomic number, and content is divided by the absorption interaction element's ray energy to obtain a second correction coefficient; The second correction coefficients of each absorption interaction element are added to obtain a third correction coefficient; The first correction coefficient and the third correction coefficient are summed to obtain the first element's ray self-absorption correction coefficient; When adjusting the flux distribution based on the neutron absorption cross section, the method comprises the following steps: When the neutron absorption cross section is greater than the neutron absorption cross section threshold, the flux distribution is determined by using a Monte Carlo simulation method; When the neutron absorption cross section is less than or equal to the neutron absorption cross section threshold, the flux distribution is equal to a standard flux distribution; When the flux distribution is determined by using the Monte Carlo simulation method, the method comprises the following steps: Neutron interaction elements that have interaction with neutrons in the measured sample are extracted; The neutron scattering cross section, neutron absorption cross section, and content of the neutron interaction element are obtained; The product of the neutron scattering cross section, neutron absorption cross section, and content of the neutron interaction element is divided by a preset neutron energy constant to obtain a first distribution parameter; The first distribution parameters of each neutron interaction element are added to obtain a second distribution parameter; The second distribution parameter is substituted into a Monte Carlo simulation model to obtain the flux distribution; When calibrating the peak position based on the environmental temperature variation, the method comprises the following steps: When the environmental temperature variation is greater than the temperature variation threshold, the peak position calibration value is determined by using a temperature segmented calibration method; When the environmental temperature variation is less than or equal to the temperature variation threshold, the peak position calibration value is equal to a standard peak position calibration value; When the peak position calibration value is determined by using the temperature segmented calibration method, the method comprises the following steps: The environmental temperature is divided into a plurality of continuous temperature intervals; Determine a temperature interval to which a current environment temperature belongs; Obtain a peak temperature coefficient corresponding to the temperature interval and a standard peak value; Multiply the peak temperature coefficient, a difference between the current environment temperature and a standard temperature, and the standard peak value to obtain a peak calibration value; When the environment humidity is greater than a humidity threshold value, a background correction amount is determined by using a humidity dynamic correction mode; When the environment humidity is less than or equal to the humidity threshold value, the background correction amount is equal to a standard background correction amount; When the background correction amount is determined by using the humidity dynamic correction mode, the method comprises the following steps: Obtain an initial humidity and a current humidity of the measured sample; Calculate a difference between the current humidity and the initial humidity to obtain a humidity change amount; Obtain a preset humidity correction coefficient and a standard background spectrum; Multiply the humidity change amount, the humidity correction coefficient, and the standard background spectrum to obtain the background correction amount. ​

Images