Method for quick simulation and complement of airborne data based on ground gamma-ray spectral data

Abstract

The present application relates to a kind of based on ground gamma energy spectrum data's flight data quick simulation complementing method, comprising: arranging existing aerial gamma energy spectrum data in research area, the average height h of statistics measurement is extracted data missing range;Collect ground gamma energy spectrum data information in data missing area, and obtain the potassium, uranium, thorium count rate of ground 0 m height;The count rate data of ground 0 m height is gridded, and artificial generation flight path data is extracted and established ground 0 m height count rate database;Obtain height attenuation coefficient, the aerial sensitivity coefficient of height h;According to ground 0 m height count rate database, according to the response contribution of measuring point data in certain area, integral obtains simulated aerial count rate data;According to simulated aerial count rate data, obtain aerial gamma energy spectrum data.The present application can quickly complement the data of missing area, uses ground gamma energy spectrum measurement to replace aerial gamma energy spectrum measurement, greatly saves cost.

Description

TECHNICAL FIELD

[0001] The application relates to a radioactive data processing method, in particular to a fast simulation and complementation method of airborne radioactive data based on ground gamma spectrum data. BACKGROUND

[0002] With the increasingly high requirement for radioactivity environment, the fusion and development of radioactive data are more and more urgent. In the existing work in a certain region, different regions have carried out airborne gamma spectrum measurement work with different purposes and ranges, and obtained the corresponding airborne gamma spectrum data of different regions. However, the ranges of airborne gamma spectrum measurement in different regions cannot be completely connected, resulting in the existence of a blank area in a local region, which leads to the incomplete utilization of regional airborne gamma spectrum data, and is not conducive to the development of regional radioactive analysis and radiation environment evaluation. The cost of carrying out airborne gamma spectrum measurement in the blank area is high, and it is urgent to explore a fast and effective airborne gamma spectrum data processing method and technology, and further study the method of simulating and complementing the airborne gamma spectrum data, which has important significance for the development and utilization of existing airborne gamma spectrum data, cost reduction and environmental damage reduction caused by exploration. SUMMARY

[0003] The purpose of the application is to provide a fast simulation and complementation method of airborne radioactive data based on ground gamma spectrum data, so as to solve the problem of high cost in the simulation and complementation of airborne gamma spectrum data in the blank area.

[0004] The application is implemented as follows: a fast simulation and complementation method of airborne radioactive data based on ground gamma spectrum data, comprising the following steps.

[0005] a. For the delineated research area, collate the existing airborne gamma spectrum data in the research area, count the measurement average height h of the existing airborne gamma spectrum data, and extract the range of data loss to obtain a data loss area.

[0006] b. Collect the ground gamma spectrum data of the data loss area, and calculate the potassium, uranium and thorium count rates at the ground 0m height.

[0007] c. Grid the count rate data at the ground 0m height obtained in step b, artificially generate flight path data in the gridded drawing, and extract and establish the count rate database at the ground 0m height according to the flight path.

[0008] d. Select a region where the ground gamma spectrum data and the airborne gamma spectrum data exist at the same time and are in a calm background field, calculate the height attenuation coefficient mu and the aerial sensitivity coefficient S of the height h according to the data of the region.

[0009] e. According to the count rate database of the ground 0m height obtained in step c, the response contribution of the data of the measuring points in a certain area is integrated to obtain the simulated aerial count rate data.

[0010] f. According to the simulated aerial count rate data, aerial gamma spectrum data is obtained, so that the aerial gamma spectrum data of the entire study area is supplemented.

[0011] Further, in step a, the existing aerial gamma spectrum data includes aerial gamma spectrum potassium, uranium, thorium count rate and content value with coordinates and measuring height.

[0012] Further, in step b, the ground gamma spectrum data includes potassium content, uranium content, thorium content and instrument sensitivity at the ground 0m height, and the potassium, uranium and thorium count rates at the ground 0m height are obtained according to the potassium content, uranium content and thorium content multiplied by the corresponding sensitivity.

[0013] Further, in step c, the flight path line spacing, the flight path of the aerial survey is referenced to the existing aerial gamma spectrum data, the flight path point spacing is 20m-50m, the ground gamma spectrum data including potassium content, uranium content, thorium content and coordinate information is extracted according to the generated flight path, and a count rate database at the ground 0m height is formed.

[0014] Further, in step d, the height attenuation coefficient μ = (lnN 地面 -lnN 空中 ) / h, wherein N 地面 is the ground count rate, N 空中 is the aerial survey count rate, and h is the simulated measuring height; the aerial sensitivity coefficient S = N 空中 / Q 地面 , wherein N 空中 is the aerial survey count rate, Q 地面 is the known content on the ground, and the aerial sensitivity coefficient S distribution includes SK, SU and STh.

[0015] Further, in step e, the ground measuring points with a horizontal distance of less than 180m to the simulated aerial measuring point participate in the calculation, and the calculation formula of the simulated aerial count rate data is as follows:

[0016]

[0017] , wherein N 地面 is the ground count rate; N 空中 is the aerial count rate; μ is the height attenuation coefficient; n is the number of known ground measuring points within the 180m radius of the ground projection of the simulated aerial measuring point; and r is the distance of the simulated measuring point.

[0018] Further, in step f, the calculation formula of the aerial gamma spectrum data is as follows:

[0019] Q=S*N 空中 .

[0020] Wherein: Q is the calculated content; N 空中 is the simulated air measurement point count rate; S is the air sensitivity coefficient.

[0021] In step b, the vehicle-mounted gamma energy spectrum measurement is carried out along the highway to obtain the ground gamma energy spectrum data of the data missing area.

[0022] The application carries out the measurement of the ground gamma energy spectrum in the data missing area, obtains the potassium, uranium and thorium count rates of the ground 0m height according to the ground gamma energy spectrum, carries out the gridding of the count rate data of the ground 0m height, artificially generates the flight path data, extracts and establishes the count rate database of the ground 0m height according to the flight path, integrates the simulated air count rate value according to the count rate database of the ground 0m height, and then obtains the aerial gamma spectrum data (including the potassium, uranium and thorium content) according to the simulated air count rate value, so as to supplement the aerial gamma spectrum data of the data missing area. At the same time, by selecting the area where the ground data and the aerial data exist simultaneously and is a quiet background field, the required height attenuation coefficient and the air sensitivity coefficient of the height h are calculated.

[0023] The application quickly supplements the blank area in the aerial gamma spectrum data by the ground gamma spectrum data, solves the problem that some small areas lack aerial gamma spectrum data, adopts the ground gamma spectrum measurement to replace the aerial gamma spectrum measurement, greatly saves the cost, and can reduce the environmental damage caused by exploration. BRIEF DESCRIPTION OF DRAWINGS

[0024] Figure 1 It is a schematic diagram for calculating the simulated air count rate data of the application.

[0025] Figure 2 It is the ground potassium content equivalent color block diagram of the example provided by the application.

[0026] Figure 3 It is the simulated aerial potassium count rate equivalent color block diagram of the example provided by the application. DETAILED DESCRIPTION

[0027] The application is a kind of based on the ground gamma spectrum data's aerial gamma spectrum data fast simulation method for supplementing, and the steps include the following contents.

[0028] a. For the delineated research area, collate the existing aerial gamma spectrum data in the research area, count the measurement average height h of the existing aerial gamma spectrum data, and obtain the data missing area by extracting the range of data missing.

[0029] Specifically, the existing aerial gamma-ray spectrometry database information of the study area is collected, the database including the count rate and content value of potassium, uranium and thorium with coordinates and measured height. The average height h of the measured area is calculated, the data is projected to a certain coordinate system (Mapgis map), and the coordinate information of the data missing area is extracted to obtain the data missing area.

[0030] b. Collecting the ground gamma-ray spectrometry data of the data missing area, and calculating the count rate of potassium, uranium and thorium at the height of 0 m.

[0031] Specifically, the ground gamma-ray spectrometry data includes the content of potassium, uranium and thorium at the height of 0 m and the instrument sensitivity. The count rate of potassium, uranium and thorium at the height of 0 m is calculated according to the content of potassium, uranium and thorium multiplied by the corresponding sensitivity.

[0032] The ground gamma-ray spectrometry data of the data missing area can be obtained by carrying out vehicle-mounted gamma-ray spectrometry along the road in the data missing area.

[0033] c. Griding the count rate data at the height of 0 m obtained in step b, manually generating flight path data in the gridded map, and extracting and establishing the count rate database at the height of 0 m according to the flight path.

[0034] Specifically, the count rate data at the height of 0 m is gridded, and the gridding interval is selected as 1 / 2 of the aerial gamma-ray spectrometry interval outside the data missing area.

[0035] Then the flight path data is manually generated, the flight path line interval and direction are referenced to the aerial flight path of the existing aerial gamma-ray spectrometry data, preferably the extension line of the existing aerial flight path of the existing aerial gamma-ray spectrometry data, and the flight path point interval on the flight path line is 20 m to 50 m.

[0036] The ground gamma-ray spectrometry data including the content of potassium, uranium and thorium and coordinate information is extracted from the count rate database at the height of 0 m according to the generated flight path, to form the count rate database at the height of 0 m. The database is similar to the existing aerial gamma-ray spectrometry data format, which is convenient for further data processing.

[0037] d. Selecting an area where the ground gamma-ray spectrometry data and the aerial gamma-ray spectrometry data exist and the background field is calm, and calculating the height attenuation coefficient μ and the aerial sensitivity coefficient S of the height h according to the data of the area.

[0038] The height attenuation coefficient and the aerial sensitivity coefficient of the height h are calculated using the ground gamma-ray spectrometry data and the aerial gamma-ray spectrometry data, the height attenuation coefficient μ = (lnN 地面 -lnN 空中 ) / h, wherein N 地面 is the ground count rate, and N 空中N is the aerial survey count rate, h is the simulated measurement height, which is the same as the average measurement height of the existing aerial gamma-ray spectrum data; the aerial sensitivity coefficient S = N 空中 / Q 地面 wherein: N 空中 is the aerial survey count rate, Q 地面 is the known content on the ground.

[0039] e. According to the count rate database of the ground 0m height obtained in step c, the simulated aerial count rate data is obtained by integrating the data response contribution of the measuring points within a certain area.

[0040] Specifically, the simulated aerial count rate value is integrated by using the count rate database of the ground 0m height, and the ground measuring points within a distance of less than 180m from the simulated aerial measuring point are involved in the calculation (the contribution of the surface source within a certain range on the ground), while satisfying the characteristics of ray exponential attenuation, i.e. N 空中 =N 地面 e -μh (wherein μ is the height attenuation coefficient, and h is the simulated measurement height).

[0041] As shown in Figure 1 , the ground measuring points within a horizontal distance of less than 180m from the simulated aerial measuring point are involved in the calculation, and the calculation formula of the simulated aerial count rate data is as follows:

[0042]

[0043] wherein: N 地面 is the ground count rate; N 空中 is the aerial count rate; μ is the height attenuation coefficient; n is the number of known ground measuring points within a 180m radius range of the ground projection of the simulated aerial measuring point; and r is the distance of the simulated measuring point.

[0044] f. According to the simulated aerial count rate data, the aerial gamma-ray spectrum data is obtained, so as to complete the aerial gamma-ray spectrum data of the entire study area.

[0045] The calculation formula of the aerial gamma-ray spectrum data is as follows:

[0046] Q = S * N 空中 .

[0047] wherein: Q is the calculated content; N 空中 is the count rate of the simulated aerial measuring point; and S is the aerial sensitivity coefficient.

[0048] As shown in Figure 2 and Figure 3 , an example is provided, in which the ground gamma-ray spectrum is measured in the study area to obtain data, and then the data is processed to obtain Figure 1The ground potassium content equivalent color block diagram is obtained, and then corresponding aerial gamma spectrum data is calculated by the method of the application, so that Figure 2 The simulated aerial potassium count rate equivalent color block diagram is obtained.

[0049] The application measures the ground gamma spectrum in the data missing area, obtains the potassium, uranium and thorium count rates at the ground 0m height according to the ground gamma spectrum, grids the count rate data at the ground 0m height, artificially generates the flight path data, extracts and establishes the count rate database at the ground 0m height according to the flight path, integrates the simulated aerial count rate value according to the count rate database at the ground 0m height, and then obtains the aerial gamma spectrum data (including the potassium, uranium and thorium content) according to the simulated aerial count rate value, so as to supplement the aerial gamma spectrum data in the data missing area. Meanwhile, the height attenuation coefficient and the aerial sensitivity coefficient of the height h are calculated by selecting the area where the ground data and the aerial data exist simultaneously and the background field is quiet.

[0050] The application quickly supplements the blank area in the aerial gamma spectrum data by the ground gamma spectrum data, adopts the ground gamma spectrum measurement to replace the aerial gamma spectrum measurement, greatly saves the cost, and can reduce the environmental damage caused by exploration.

Claims

1. A method for quick simulation and filling of airborne data based on ground gamma-ray spectral data, characterized in that, The method comprises the following steps: a. For the delineated study area, collate the existing aerial gamma-ray spectrum data in the study area, count the measurement average height h of the existing aerial gamma-ray spectrum data, and extract the range of missing data to obtain a data missing area; b. Collect ground gamma-ray spectrum data information in the data missing area, and calculate the potassium, uranium, and thorium count rates at a ground height of 0 m; c. Grid the count rate data at a ground height of 0 m obtained in step b, manually generate flight path data in the gridded map, and extract and establish a count rate database at a ground height of 0 m according to the flight path; d. Select an area where ground gamma-ray spectrum data and aerial gamma-ray spectrum data exist simultaneously and is a calm background field, calculate the height attenuation coefficient μ and the aerial sensitivity coefficient S of the height h according to the data of the area; e. According to the count rate database at a ground height of 0 m obtained in step c, integrate to obtain simulated aerial count rate data according to the response contribution of the data of measuring points in a certain area; f. Calculate aerial gamma-ray spectrum data according to the simulated aerial count rate data, so as to complete the aerial gamma-ray spectrum data of the entire study area.

2. The method according to claim 1, wherein, In step a, the existing aerial gamma-ray spectrum data includes aerial gamma-ray spectrum potassium, uranium, and thorium count rates and content values with coordinates and measurement heights.

3. The method of claim 1, wherein the method is characterized by: In step b, the ground gamma-ray spectrum data includes potassium content, uranium content, thorium content, and instrument sensitivity at a ground height of 0 m, and the potassium, uranium, and thorium count rates at a ground height of 0 m are calculated according to the potassium content, uranium content, and thorium content multiplied by the corresponding sensitivity.

4. The method of claim 1, wherein the method is characterized by, In step c, the flight path line spacing and direction refer to the aerial survey flight path of the existing aerial gamma-ray spectrum data, and the flight path point spacing is 20 m to 50 m. The ground gamma-ray spectrum data, including potassium count rate, uranium count rate, thorium count rate, and coordinate information, are extracted according to the generated flight path to form a count rate database at a ground height of 0 m.

5. The method of claim 1, wherein, In step d, the ground data and the air data section are repeated section, and the height attenuation coefficient μ = (lnN 地面 -lnN 空中 ) / h is calculated, wherein N 地面 is the ground count rate, N 空中 is the air count rate, and h is the simulation measurement height. The air sensitivity coefficient S = N 空中 / Q 地面 is calculated, wherein N 空中 is the air count rate, and Q 地面 is the known content of the ground.

6. The method of claim 1, wherein, In step e, the horizontal distance from the simulated aerial measuring point to the ground measuring point is less than 180 m, and the ground measuring point participates in the calculation. The calculation formula of the simulated aerial count rate data is as follows: Where: N 地面 is the ground count rate; N 空中 is the aerial count rate; μ is the height attenuation coefficient; n is the number of known ground measurement points within a 180 m radius of the ground projection of the simulated aerial measurement point; and r is the distance of the simulated measurement point.

7. The method of claim 1, wherein, In step f, the calculation formula of the aerial gamma-ray spectrum data is as follows: Q = S * N 空中； Where: Q is the calculated content; N 空中 is the simulated air measurement point count rate; S is the air sensitivity coefficient.

Images