Method for extracting radioactivity spectrum peak area for area radiation monitor

Abstract

The present application relates to the technical field of energy spectrum peak analysis, and particularly relates to a method for extracting radioactivity spectrum peak area of a regional radiation monitor, determining a center peak channel position of energy spectrum data, evaluating the degree of subordination of a channel position to the center peak channel position by analyzing the similarity of energy characteristics at the local distribution channel position, and then preliminarily dividing the energy spectrum data into energy intervals, performing adaptive division and effective boundary identification based on the energy distribution characteristics of the energy spectrum in the energy interval, and obtaining a net peak area according to the effective boundary. The present application determines a more accurate range of the radioactivity performance peak by analyzing the change process of the energy characteristics from the peak center to each channel position in the energy interval, and improves the accuracy of positioning the radioactivity performance peak range in a complex energy environment.

Description

Technical Field

[0001] This invention relates to the field of energy spectrum peak analysis technology, and specifically to a method for extracting the area of ​​radioactive energy spectrum peaks for use in regional radiation monitoring instruments. Background Technology

[0002] With the widespread application of nuclear technology in energy, medical, and industrial fields, the need for continuous and accurate monitoring of environmental radioactivity levels has become increasingly important. Regional radiation monitoring instruments identify and quantify the radionuclides and their activities by collecting energy spectral data from the environment. In energy spectrum analysis, the area of ​​characteristic peaks is directly related to the activity of the corresponding radionuclides. Therefore, analyzing and identifying radioactive peaks from complex energy spectra and extracting their net peak areas is the core step in the entire regional radiation detection and analysis process. Existing technologies often use the Covell method to calculate the area of ​​energy spectral peaks.

[0003] The existing Covell method takes n channels on each side of the peak as the center and uses them as boundaries. After determining the peak and boundary channels, it uses a straight line to connect the count values ​​of the boundary points as the background line to obtain the background area. Then, it obtains the peak area by combining the total peak area and the background area. However, the Covell method is extremely sensitive to the choice of the n value when obtaining the peak area. It often uses 0.5 to 1 times the peak half-width as its n value. However, for regions with complex energy spectrum backgrounds, simply determining the n value by the half-width of the radiometric peak ignores the changes in other energies in the environment, which can easily introduce additional biases and lead to errors in the peak area calculation. Summary of the Invention

[0004] To address the technical problem in existing technologies where radioactive energy peaks are mixed with some energy peaks in the environment, easily introducing environmental background into the defined full-energy peak boundaries, thus leading to inaccurate peak areas determined by the Covell method, the present invention aims to provide a method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument. The specific technical solution adopted is as follows: This invention provides a method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument, comprising: Energy spectrum data is acquired at the monitoring point using a detector, and the central peak location is determined. Based on the change in energy intensity between each channel address and the nearest central peak channel address, the signal peak difference of each channel address is determined; based on the signal peak difference, the distance distribution and energy distribution characteristics between the two channels address and the nearest central peak channel address are analyzed to obtain the approximation of the energy characteristics of the two channels address. Based on the similarity of energy characteristics between each channel address and other corresponding channels address corresponding to the nearest central peak channel address, a subordinate evaluation index for each channel address relative to the nearest central peak channel address is obtained; based on the subordinate evaluation index, the channels address are clustered with the central peak channel address as the center to obtain the energy range. Within each energy range, the energy variation difference of each channel is determined based on the degree of deviation of the energy characteristics between each channel and other channels. The effective boundary of each energy range is determined based on the local distribution of the energy variation difference. The net area of ​​the energy spectrum peak is extracted based on the analysis of the net peak area of ​​the effective boundary.

[0005] Furthermore, the method for obtaining the signal peak difference includes: For any channel address, the difference in signal strength between the channel address and the nearest center peak channel address is taken as the center strength deviation of the channel address; the ratio of the center strength deviation of the channel address to the signal strength of the nearest center peak channel address is taken as the signal peak difference degree of the channel address.

[0006] Furthermore, the method for obtaining the approximation of the energy characteristics includes: For any central peak address, any two different addresses that are the nearest central peak address are taken as the analysis address group, and the address distance between each address in the analysis address group and the central peak address is obtained as the distribution distance of the address. The distribution distance difference between two road addresses in the road address group and the road address difference between two road addresses are combined and used as the distribution approximation of the road address group. Both the distribution distance difference and the road address difference are negatively correlated with the distribution approximation. By combining the difference in signal peaks and the similarity in distribution between two channels in the analysis channel group, as well as the distribution distance, the similarity in energy characteristics of the analysis channel group is obtained. The difference in signal peaks and the similarity in energy characteristics are negatively correlated.

[0007] Furthermore, the method for obtaining the subordinate evaluation indicators includes: For any given road address, the road address distance between that road address and the corresponding nearest central peak road address is used as the distance subordinate index for that road address. There is a negative correlation between the road address distance and the distance subordinate index. For any other channel address that is the closest to the central peak channel address, the channel distance between the channel address and the other channel address is negatively correlated and normalized to serve as the distance coefficient between the channel address and the other channel address; the product of the energy characteristic similarity between the channel address and the other channel address and the distance coefficient is used as the energy characteristic subordination index between the channel address and the other channel address. By combining the distance subordination index and energy characteristic subordination index of this tunnel site with all other tunnel sites, the subordination evaluation index of this tunnel site is obtained.

[0008] Furthermore, the method for obtaining the energy range includes: Using the central peak tunnel site as the cluster center, the subordinate evaluation index between the tunnel site and the center point is used as the metric for clustering. Each cluster result is used as an energy interval, and the total number of energy intervals is the same as the total number of central peak tunnel sites.

[0009] Furthermore, the method for obtaining the energy change difference includes: For any energy range, the mean of the energy characteristic approximation between all pairs of addresses in the energy range is negatively correlated and mapped as the background inconsistency mean of the energy range. For any address in the energy range, perform a negative correlation mapping between the address and the maximum energy characteristic similarity between the address and other addresses in the energy range to obtain the minimum inconsistency of the address. The difference between the minimum inconsistency of the channel address and the background inconsistency uniformity is taken as the energy variation difference of the channel address in the energy range.

[0010] Furthermore, the method for obtaining the effective boundary includes: For any side of the central peak address, starting from the central peak address, traverse each address one by one, and count the proportion of the number of positive energy change differences in the preset analysis window. When the proportion of the number is greater than the preset proportion threshold, stop traversing and take the address corresponding to the stop as the effective boundary; there are two effective boundaries in each energy interval.

[0011] Furthermore, the extraction of the net area of ​​the energy spectrum peak based on the net peak area analysis of the effective boundary includes: For any energy range, connect the signal intensity points corresponding to the two effective boundaries of the energy range, and obtain the background area based on the connection and the boundary; obtain the total area based on the energy spectrum data of all channels in the effective boundary; obtain the net area by the difference between the total area and the background area.

[0012] Furthermore, the method for obtaining the nearest central peak address includes: For any given channel address, calculate the channel address difference between that channel address and each center peak channel address, and take the center peak channel address with the smallest channel address difference as the nearest center peak channel address for that channel address.

[0013] Furthermore, the horizontal axis of the energy spectrum data represents the channel address, and the vertical axis represents the signal strength.

[0014] This invention offers the following advantages: By determining the location of the central peak and leveraging the similarity of energy characteristics at local distribution points, the entire energy spectrum is initially divided into energy intervals. Adaptive segmentation and boundary identification are performed based on the spectrum's inherent energy distribution characteristics, reducing the impact of single parameters. This allows for more accurate reduction of background contributions for peaks with different locations, resulting in a more realistic net peak area and lowering the risk of misjudgment due to environmental changes or fluctuations in spectrum quality. Furthermore, by analyzing the changes in energy characteristics from the peak center to each point within the energy interval, this invention determines a more accurate range for radioactive peaks, improving the accuracy of locating radioactive peak ranges in complex energy environments. Attached Figure Description

[0015] Figure 1 A flowchart of a method for extracting the peak area of ​​a radioactive energy spectrum for a regional radiation monitoring instrument, provided by the present invention; Figure 2 This invention provides a schematic diagram for obtaining the center peak address superimposed with high noise. Detailed Implementation

[0016] like Figure 1 As shown, an embodiment of the present invention provides a method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument, comprising the following steps: S1: At the monitoring point, energy spectrum data is obtained through a detector, and the central peak location is determined.

[0017] In this embodiment of the invention, a suitable detector is placed, and three standard sources with known energies (such as Cs 137662 keV, Co60 1173 / 1332 keV) are selected to calibrate the detector's energy scale. The calibrated detector is then used to monitor the target area. During monitoring, points are selected using a grid-based method, and the detector is placed in these points for a period of time to obtain radiation information from the environment. The signal is then classified according to amplitude using MCA to form an initial energy spectrum, which is then filtered and denoised using a moving average method to obtain the energy spectrum data. The monitoring placement time can be preset to half an hour, and the specific monitoring setup can be adjusted by the implementer according to the specific implementation scenario; no restrictions are imposed here.

[0018] It should be noted that the classification data is processed by a multichannel analyzer (MCA). The MCA divides the pulse amplitude into several "intervals," i.e., channel addresses, and counts the number of pulses within each interval, i.e., the signal intensity. The final output is two-dimensional data of channel address versus intensity. Therefore, the horizontal axis of the energy spectrum data represents the channel address, indicating the X-ray signal energy (keV). Different elements require different excitation energies for electronic transitions; therefore, the element type can be determined from the horizontal axis position of the spectral peak. The vertical axis of the energy spectrum data represents the signal intensity, reflecting the energy intensity. Based on the detected signal intensity, it is compared with stored standard values, and after correction, the elemental content can be obtained.

[0019] The energy spectrum contains several radiation peaks, reflecting energy intensity information when subjected to changes in energy intensity at corresponding frequencies. Due to the presence of background energy, complex environments can affect the determination of the area of ​​the radiation spectrum. To improve the determination of the range of the radiation peaks, the center position of the peak is first determined. When the energy intensity reaches its peak, it will change from rising to falling. Therefore, in this embodiment of the invention, the peak value can be obtained by the derivative method, and the channel address corresponding to the peak value is the center peak channel address. It should be noted that the derivative method for determining the peak value is a well-known technique in the art, that is, the point where the first derivative is zero and the second derivative is less than zero is taken as the peak value, which will not be further elaborated here.

[0020] Specifically, due to the possibility of spurious peaks caused by noise interference, when the distance between peak addresses is less than the preset interval, the average of the corresponding peak addresses is used as the center peak address. The preset interval can be set to 3, which can be adjusted by the implementer and is not limited here. Please refer to [link to relevant documentation]. Figure 2 It shows a schematic diagram of obtaining the center peak address superimposed with high noise, provided by the present invention. Figure 2 As shown, due to the influence of noise, there may be multiple false peaks on the same true peak (the bands between the boundaries of the true peaks are considered as one true peak). By averaging to find the center, the influence of false peaks can be eliminated as much as possible, and a more accurate center peak channel address can be obtained.

[0021] S2: Based on the change in energy intensity between each channel address and the nearest central peak channel address, determine the signal peak difference degree of each channel address; based on the signal peak difference degree, analyze the distance distribution and energy distribution characteristics between the two channels address and the nearest central peak channel address, and obtain the approximation degree of the energy characteristics of the two channels address.

[0022] The energy spectrum contains many peaks. Since the energy types corresponding to each peak are similar, that is, the emitted rays or energies are within the same small energy range, in order to accurately determine the energy distribution range of each peak, a preliminary analysis of the approximate situation between the channel sites is conducted for each central peak.

[0023] Preferably, in this embodiment of the invention, the method for obtaining the signal peak difference includes: For any given path address, the path address of the nearest center peak is taken as the nearest center peak path address of that path address. In this embodiment of the invention, the path address difference between the path address and each center peak path address is calculated, and the center peak path address corresponding to the path address with the smallest path address difference is taken as the nearest center peak path address of that path address, so as to facilitate subsequent measurement of the intensity deviation between the path address location and the target expected center peak.

[0024] The difference in signal strength between this channel address and the nearest center peak channel address is taken as the center strength deviation of this channel address, reflecting the degree of strength deviation between the channel address and the target center peak. The smaller the center strength deviation, the more likely the channel address belongs to the nearest center peak. Furthermore, the ratio of the center strength deviation of this channel address to the signal strength of the nearest center peak channel address is taken as the signal peak difference degree of this channel address, and the degree of deviation is uniformly quantified by the signal strength of the nearest center peak channel address.

[0025] Specifically, if multiple center peak addresses correspond to the smallest channel address difference, the center intensity deviation is calculated for each address, and the one with the smallest center intensity deviation is taken as the nearest center peak address. If there are still channels with consistent center intensity deviations, it indicates that the channel address has no deviation in its response to the center peak and has virtually no impact on peak affiliation determination. In this case, any center peak address can be arbitrarily selected as the nearest center peak address. In this embodiment of the invention, the difference is obtained through the absolute value of the numerical difference.

[0026] The evaluation of environmental energy characteristics is based on a combination of the location and energy intensity. When the energy intensity is similar and the frequency distribution is also similar, it is more likely that the two locations belong to the same energy peak. Therefore, the similarity of the energy characteristics of the two locations is analyzed by combining the energy distribution and the distance distribution between the locations.

[0027] Preferably, in this embodiment of the invention, the method for obtaining the approximation of energy characteristics includes: For any given central peak address, any two distinct addresses that are the closest central peak address are used as the analysis address group. The address distance between each address in the analysis address group and the central peak address is obtained as the distribution distance of the address. For the same nearest central peak, the closer the distance between the address and the central peak address, the more likely the trend state of the two addresses is consistent with the central peak, and the more consistent the similar state is.

[0028] Furthermore, the distribution distance difference between two road addresses in the analysis road address group and the road address difference between two road addresses are used as the distribution approximation of the road address group. Both the distribution distance difference and the road address difference are negatively correlated with the distribution approximation. The smaller the road address difference, the more likely the distribution is to be a continuous and concentrated state, and the higher the approximation. In a specific embodiment of the invention, the inverse ratio of the product of the distribution distance difference between two road addresses in the analysis road address group and the road address difference (specifically, the reciprocal in this embodiment) is used as the distribution approximation of the road address group. The smaller the difference, the higher the similarity of the features in the location distribution. If the product is zero, the distribution approximation can be directly assigned a value of 1.

[0029] Finally, by combining the differences in signal peak intensity and distribution similarity between two channels in the analysis group, as well as the distribution distance, the energy characteristic similarity of the analysis group was obtained. The differences in signal peak intensity and the energy characteristic similarity were negatively correlated. The smaller the differences in signal peak intensity, the more consistent the attenuation intensity relative to the peak center between the two channels, and the higher the similarity in energy intensity characteristics. Furthermore, the smaller the distribution distance, the closer the channel analysis is to the central peak, and the more significant the attribution characteristics.

[0030] In one specific embodiment of the present invention, the inverse ratio of the product of the distribution distances of two channel addresses in the analysis channel address group is used as the distance distribution index. If the product of the distribution distances of the two channel addresses is zero, the distance distribution index can be directly assigned a value of 1. The inverse ratio of the absolute value of the difference in signal peak difference between two channel addresses in the analysis channel address group is used as the deviation approximation. If the absolute value of the difference is zero, the deviation approximation can be directly assigned a value of 1. The product of the deviation approximation, the distribution approximation, and the distance distribution index is used as the energy characteristic approximation of the analysis channel address group. The larger the energy characteristic approximation, the more consistent the two channel addresses' responses to the same peak, and the more likely they are channel addresses within the energy peak.

[0031] S3: Based on the similarity of energy characteristics between each channel address and other corresponding channels address of the nearest central peak channel address, obtain the subordinate evaluation index of each channel address to the nearest central peak channel address; based on the subordinate evaluation index, cluster the channels address with the central peak channel address as the center to obtain the energy range.

[0032] For an energy peak, considering the peak's distribution characteristics, the energy characteristics of the channel address at the peak concentration point are highly similar. The closer the channel address is to the peak center, the higher the overall similarity of its energy characteristics with other channel addresses, which indicates that the channel address is more likely to belong to the corresponding central peak.

[0033] Therefore, considering the approximate energy characteristics and the possibility of subordination, in this embodiment of the invention, the method for obtaining the subordination evaluation index includes: For any non-central peak address, the distance between that address and the corresponding nearest central peak address is used as the distance subordination index. There is a negative correlation between the address distance and the distance subordination index. For distributed distances, the smaller the difference between addresses (i.e., the closer the distances), the higher the probability of address assignment. In one specific embodiment of the invention, the reciprocal of the address difference between that address and the corresponding nearest central peak address is used as the distance subordination index.

[0034] Furthermore, for any other channel that is the closest central peak channel to the current channel, the channel distance between the current channel and the other channel is negatively correlated and normalized, and used as the distance coefficient between the current channel and the other channel. The similarity of the energy characteristics between the current channel and the other channels that are the closest central peak channels is analyzed. Since the intensity change of the peak is continuous, the closer the distance between the current channel and other channels, the more important the similarity analysis should be, and the more relevant the attribution should be.

[0035] In one specific embodiment of the present invention, an exponential function with the natural constant e as its base is utilized. Perform negative correlation mapping and normalization, i.e. , The independent variable is the distance between the road and the road surface. This is also known as the distance coefficient.

[0036] Therefore, the product of the energy characteristic similarity between the address and other addresses and the distance coefficient is used as the energy characteristic subordination index between the address and other addresses. The greater the energy characteristic similarity, the higher the consistency of subordination, and the greater the probability that the address is within the range of the nearest central peak.

[0037] Finally, by comprehensively analyzing the relationship between the tunnel and all other tunnel sites, and combining the distance subordination index and energy characteristic subordination index of the tunnel site with all other tunnel sites, a subordination evaluation index for the tunnel site is obtained. In a specific embodiment of the present invention, after calculating the product of the distance subordination index and the energy characteristic subordination index of the tunnel site with each other tunnel site, the mean of all products is used as the subordination evaluation index of the tunnel site. The larger the subordination evaluation index, the higher the probability that the tunnel site belongs to the nearest central peak.

[0038] Based on the affiliation assessment, local clustering is performed on each central peak channel address to initially divide different energy peaks. The higher the affiliation probability, the more likely they are to be the same peak. In this embodiment of the invention, the central peak channel address is used as the cluster center point, and the affiliation assessment index between the channel address and the center point is used as a metric for clustering. The higher the affiliation assessment index, the higher the probability of them belonging to the same cluster. Each clustering result is used as an energy interval, and each cluster corresponds to an energy interval. Each energy interval contains one central peak channel address, and the total number of energy intervals is the same as the total number of central peak channel addresses.

[0039] It should be noted that clustering is a technique well-known to those skilled in the art, such as K-means clustering, and will not be limited or elaborated upon here.

[0040] S4: In each energy range, determine the degree of energy variation difference for each channel address based on the degree of deviation of the energy characteristic approximation between each channel address and other channel addresses; determine the effective boundary of each energy range based on the local distribution of the degree of energy variation difference in the energy range; and extract the net area of ​​the energy spectrum peaks based on the analysis of the net peak area according to the effective boundary.

[0041] For an energy range, the background energy is generated by many environmental factors, but the radioactive radiation energy is generated by the same factor. Therefore, the radiation energy characteristics in the energy peak are significantly different from other background energy characteristics. Subsequently, the boundary between this energy peak and the environmental energy can be determined based on the energy difference changes in the energy range.

[0042] Therefore, based on the overall energy characteristic approximation distribution among the trace addresses in the interval, the deviation state of each trace address in the interval is analyzed. In this embodiment of the invention, the method for obtaining the energy change difference of each trace address includes: For any given energy range, a negative correlation is applied to the mean of the energy characteristic approximation between all pairwise addresses within that range. Specifically, the reciprocal of the mean of the energy characteristic approximation between all pairwise addresses within that energy range is used as the background inconsistency mean for that energy range. Considering the inconsistency in the signal changes of the background energy for any two addresses, the energy characteristic approximation between the two addresses will be smaller. Through negative correlation, the overall mean reflects the baseline level of the background addresses. If the mean of the energy characteristic approximation is zero, the background inconsistency mean can be directly assigned a value of 1.

[0043] For any channel address within this energy range, a negative correlation mapping is performed between this channel address and the maximum energy characteristic approximation of other channels within the same energy range to obtain the minimum inconsistency of that channel address. Specifically, the reciprocal of the maximum energy characteristic approximation is taken as the minimum inconsistency. If there exists a maximum energy characteristic approximation of zero, the minimum inconsistency can be directly assigned a value of 1. For channels within the energy peak, due to the concentrated extremely high approximation, the inconsistency should theoretically be smaller than that of background channels. Therefore, the case with the minimum inconsistency is selected for comparison.

[0044] The difference between the minimum inconsistency of the channel address and the background inconsistency uniformity (i.e., the difference between the minimum inconsistency and the background inconsistency uniformity) is used as the energy variation difference of the channel address within that energy range. When the minimum inconsistency is lower than the background inconsistency uniformity, it indicates that the channel address has a very high degree of energy characteristic approximation and is more likely to be at the energy peak. Since the energy peak channel address is included in the background inconsistency uniformity analysis, the inconsistency of the benchmark is lower than that of a typical background channel address. Therefore, when the minimum inconsistency is higher than the background inconsistency uniformity, the channel address is more likely to be a background channel address.

[0045] Therefore, the range is divided by the distribution of local energy change differences within the energy interval. In this embodiment of the invention, for any side of the central peak address, the system traverses address by address, starting from the central peak address, and counts the percentage of addresses with positive energy change differences in the preset analysis window. When the percentage exceeds a preset percentage threshold, the traversal stops, and the address corresponding to the stop point is taken as the valid boundary. There are two valid boundaries in each energy interval. The preset analysis window is set to a side length of 5 addresses, and the preset percentage threshold is set to 60%. The specific values ​​can be adjusted by the implementer and are not limited here.

[0046] Specifically, when the traversal stops at the boundary of the energy range, the boundary of the energy range is taken as the effective boundary. The range of the energy spectrum peaks is taken as the range of the corresponding channel addresses between the effective boundaries, and the area of ​​the radioactive energy spectrum peaks is calculated by subtracting the background conditions in the environment within the range of the energy spectrum peaks.

[0047] In this embodiment of the invention, for any energy range, a line is drawn connecting the signal intensity points corresponding to the two effective boundaries of the energy range. The portion below this line represents the background area of ​​the corresponding energy spectrum peak. The background area is obtained based on the line and the boundary, which is the area enclosed by the effective boundary and the line. The total area is obtained based on the energy spectrum data of all channels within the effective boundary, which is the total area corresponding to the energy spectrum peak range. The net area is obtained by the difference between the total area and the background area, representing the area of ​​the final radioactive energy spectrum peak.

[0048] In summary, this invention, by determining the location of the central peak and leveraging the similarity of energy characteristics at local distribution points, initially divides the entire energy spectrum into energy intervals. Based on the inherent energy distribution characteristics of the spectrum, adaptive division and boundary identification are performed, reducing the influence of single parameters. This allows for more accurate reduction of background contributions for peaks with different locations, resulting in a more realistic net peak area and lowering the risk of misjudgment due to environmental changes or fluctuations in energy spectrum quality. Furthermore, by analyzing the changes in energy characteristics from the peak center to each point within the energy interval, this invention determines a more accurate range of the radioactive peak, improving the accuracy of radioactive peak range location in complex energy environments.

Claims

1. A method for extracting the peak area of ​​a radioactive energy spectrum for use in a regional radiation monitoring instrument, characterized in that, include: Energy spectrum data is acquired at the monitoring point using a detector, and the central peak location is determined. Based on the change in energy intensity between each channel address and the nearest center peak channel address, the signal peak difference degree of each channel address is determined, including: for any channel address, the signal intensity difference between the channel address and the nearest center peak channel address is taken as the center intensity deviation of the channel address; the ratio of the center intensity deviation of the channel address to the signal intensity of the nearest center peak channel address is taken as the signal peak difference degree of the channel address; based on the signal peak difference degree, the distance distribution and energy distribution characteristics between the two channels at the nearest center peak channel address are analyzed to obtain the approximation of the energy characteristics of the two channels, including: for any center peak channel address, taking the center peak channel address as... Any two distinct channels from the nearest central peak channel are taken as the analysis channel group. The channel distance between each channel in the analysis channel group and the central peak channel is obtained as the channel distribution distance. The difference in distribution distance between two channels in the analysis channel group and the channel difference between two channels are used as the distribution approximation of the analysis channel group. Both the difference in distribution distance and the channel difference are negatively correlated with the distribution approximation. The difference in signal peak difference between two channels in the analysis channel group and the distribution approximation, as well as the distribution distance, are used to obtain the energy characteristic approximation of the analysis channel group. The difference in signal peak difference is negatively correlated with the energy characteristic approximation. Based on the similarity of energy characteristics between each tunnel address and the corresponding other tunnel addresses of the nearest central peak tunnel address, a subordinate evaluation index for each tunnel address relative to the nearest central peak tunnel address is obtained. Based on the subordinate evaluation index, the tunnel addresses are clustered with the central peak tunnel address as the center to obtain energy intervals. This includes: using the central peak tunnel address as the cluster center point, clustering is performed using the subordinate evaluation index between the tunnel address and the center point as a metric, and each clustering result is used as each energy interval. The total number of energy intervals is the same as the total number of central peak tunnel addresses. Within each energy range, the energy variation difference of each channel is determined based on the degree of deviation of the energy characteristics between each channel and other channels. The effective boundary of each energy range is determined based on the local distribution of the energy variation difference. The net area of ​​the energy spectrum peak is extracted based on the analysis of the net peak area of ​​the effective boundary.

2. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The methods for obtaining the subordinate evaluation indicators include: For any given road address, the road address distance between that road address and the corresponding nearest central peak road address is used as the distance subordinate index for that road address. There is a negative correlation between the road address distance and the distance subordinate index. For any other channel address that is the closest to the central peak channel address, the channel distance between the channel address and the other channel address is negatively correlated and normalized to serve as the distance coefficient between the channel address and the other channel address; the product of the energy characteristic similarity between the channel address and the other channel address and the distance coefficient is used as the energy characteristic subordination index between the channel address and the other channel address. By combining the distance subordination index and energy characteristic subordination index of this tunnel site with all other tunnel sites, the subordination evaluation index of this tunnel site is obtained.

3. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The method for obtaining the energy change difference includes: For any energy range, the mean of the energy characteristic approximation between all pairs of addresses in the energy range is negatively correlated and mapped as the background inconsistency mean of the energy range. For any address in the energy range, perform a negative correlation mapping between the address and the maximum energy characteristic similarity between the address and other addresses in the energy range to obtain the minimum inconsistency of the address. The difference between the minimum inconsistency of the channel address and the background inconsistency uniformity is taken as the energy variation difference of the channel address in the energy range.

4. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The method for obtaining the effective boundary includes: For any side of the central peak address, starting from the central peak address, traverse each address one by one, and count the proportion of the number of positive energy change differences in the preset analysis window. When the proportion of the number is greater than the preset proportion threshold, stop traversing and take the address corresponding to the stop as the effective boundary; there are two effective boundaries in each energy interval.

5. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The extraction of the net area of ​​the energy spectrum peak based on the net peak area analysis of the effective boundary includes: For any energy range, connect the signal intensity points corresponding to the two effective boundaries of the energy range, and obtain the background area based on the connection and the boundary; obtain the total area based on the energy spectrum data of all channels in the effective boundary; obtain the net area by the difference between the total area and the background area.

6. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The method for obtaining the nearest central peak address includes: For any given channel address, calculate the channel address difference between that channel address and each center peak channel address, and take the center peak channel address with the smallest channel address difference as the nearest center peak channel address for that channel address.

7. The method for extracting the peak area of ​​a radioactive spectrum for a regional radiation monitoring instrument according to claim 1, characterized in that, The horizontal axis of the energy spectrum data represents the channel address, and the vertical axis represents the signal strength.

Images