Soil radon concentration data processing method and system

Abstract

The invention is suitable for the technical field of data processing, and provides a soil radon concentration data processing method and system, and the method comprises the following steps: receiving radon concentration collection information; a comprehensive correction factor Fi is calculated, F = f geology (Gi) * f climate (Ti, Pi, Hi) * f vegetation (Vi), f geology (Gi) is a geology correction factor, f climate (Ti, Pi, Hi) is a climate correction factor, and f vegetation (Vi) is a vegetation correction factor; determining a corrected radon concentration value Rxi according to the comprehensive correction factor Fi and the depth attenuation function; calculating a corrected radon concentration value of each coordinate point in the regional grid based on an inverse distance weighting method and the corrected radon concentration value Rxi; and calling a standard topographic map, and performing inversion to obtain a standard radon concentration value of each coordinate point in the regional grid. By introducing the comprehensive correction factor, the collected original radon concentration value is corrected, and the distribution rule and the change trend of the radon concentration of the regional soil can be known more comprehensively and accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of data processing, and in particular to a soil radon concentration data processing method and system. Background Art

[0002] Soil radon concentration is an important indicator reflecting the distribution and enrichment of the radioactive element radon in the soil. It is of great significance to geological exploration, environmental monitoring, construction project safety assessment and other fields. At present, the commonly used method for collecting soil radon concentration mainly uses professional radon concentration detection equipment to conduct field measurements at different locations and depths to obtain a series of radon concentration values. However, these raw collected data are affected by various environmental factors. Due to the complex and changeable actual measurement environmental conditions, it is difficult to directly compare and analyze the measurement data of different situations and environments. There is currently a lack of methods to convert actual measurement data to standard environmental conditions, which is not conducive to the subsequent processing and analysis of soil radon concentration. Therefore, it is necessary to provide a soil radon concentration data processing method and system to solve the above problems. Summary of the Invention

[0003] In view of the deficiencies in the prior art, the present invention aims to provide a soil radon concentration data processing method and system to solve the problems existing in the above-mentioned background technology.

[0004] The present invention is achieved by providing a method for processing soil radon concentration data, the method comprising the following steps:

[0005] Receive radon concentration collection information, the radon concentration collection information including radon concentration values, each radon concentration value Ri corresponding to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, pressure Pi, humidity Hi and covering vegetation type Vi;

[0006] Calculate the comprehensive correction factor Fi, F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor;

[0007] Determine the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function;

[0008] Calculate the corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi;

[0009] A standard topographic map is retrieved, which includes standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location, and the standard radon concentration value of each coordinate point in the regional grid is obtained by inversion.

[0010] As a further solution of the present invention, the step of calculating the comprehensive correction factor Fi specifically includes:

[0011] Calculate the geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively;

[0012] Calculate the climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient;

[0013] Calculate the vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), δ is the vegetation suppression coefficient;

[0014] The comprehensive correction factor Fi is calculated based on the geological correction factor, climate correction factor and vegetation correction factor.

[0015] As a further solution of the present invention, the step of determining the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function specifically includes:

[0016] Determine the depth attenuation function w(di)=exp(-λdi), where λ is the radon diffusion attenuation coefficient;

[0017] The corrected radon concentration value Rxi is calculated, Rxi=Ri×Fi×w(di).

[0018] As a further solution of the present invention, the step of calculating the corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi specifically includes:

[0019] Generate evenly distributed regional grid points according to the regional range and grid resolution, and determine each coordinate point in the regional grid;

[0020] Determine the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi;

[0021] According to the collected coordinates and Rxi within the search radius, the corrected radon concentration value of the corresponding coordinate point is determined. where d ij is the horizontal distance from the collection coordinate to the grid coordinate point, p is the power parameter, which is used to control the speed at which the weight decays with distance, and N represents the number of collection coordinates within the search radius.

[0022] As a further solution of the present invention, the step of inverting to obtain the standard radon concentration value of each coordinate point in the regional grid specifically includes:

[0023] The comprehensive correction factor of each coordinate point is obtained based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location;

[0024] Convert the regional grid into a three-dimensional grid based on the depth gradient, and determine the depth attenuation function of each coordinate point in the three-dimensional grid;

[0025] Determine the standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value and the depth attenuation function of each coordinate point;

[0026] Horizontal slices, vertical sections and three-dimensional renderings are generated based on the standard radon concentration values at each coordinate point in the three-dimensional grid.

[0027] Another object of the present invention is to provide a soil radon concentration data processing system, the system comprising:

[0028] The collected information uploading module is used to receive radon concentration collected information, wherein the radon concentration collected information includes radon concentration values, and each radon concentration value Ri corresponds to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, air pressure Pi, humidity Hi and covering vegetation type Vi;

[0029] Comprehensive correction factor module, used to calculate the comprehensive correction factor Fi, F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor;

[0030] A radon concentration correction module is used to determine a corrected radon concentration value Rxi based on a comprehensive correction factor Fi and a depth attenuation function;

[0031] A modified concentration calculation module is used to calculate the modified radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the modified radon concentration value Rxi;

[0032] The standard concentration calculation module is used to retrieve the standard topographic map, which contains standard temperature, standard pressure, standard humidity, and geological structure and covering vegetation type of each location, and invert to obtain the standard radon concentration value of each coordinate point in the regional grid.

[0033] As a further solution of the present invention, the comprehensive correction factor module includes:

[0034] Geological correction factor unit, used to calculate geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively;

[0035] Climate correction factor unit, used to calculate climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient;

[0036] Vegetation correction factor unit, used to calculate vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), δ is vegetation suppression coefficient;

[0037] The comprehensive correction factor unit is used to calculate the comprehensive correction factor Fi based on the geological correction factor, climate correction factor and vegetation correction factor.

[0038] As a further solution of the present invention, the radon concentration correction module includes:

[0039] A depth attenuation function unit is used to determine a depth attenuation function w(di)=exp(-λdi), where λ is a radon diffusion attenuation coefficient;

[0040] The radon concentration value calculation unit is used to calculate the corrected radon concentration value Rxi, Rxi = Ri×Fi×w(di).

[0041] As a further solution of the present invention, the corrected concentration calculation module includes:

[0042] A regional grid determination unit is used to generate evenly distributed regional grid points according to the regional range and grid resolution, and determine each coordinate point in the regional grid;

[0043] An acquisition coordinate search unit, configured to determine the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi;

[0044] The corrected concentration calculation unit is used to determine the corrected radon concentration value of the corresponding coordinate point according to the collected coordinates within the search radius and Rxi. where d ij is the horizontal distance from the collection coordinate to the grid coordinate point, p is the power parameter, which is used to control the speed at which the weight decays with distance, and N represents the number of collection coordinates within the search radius.

[0045] As a further solution of the present invention, the standard concentration calculation module includes:

[0046] The coordinate point correction factor unit is used to obtain the comprehensive correction factor of each coordinate point based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location;

[0047] A coordinate point depth attenuation unit is used to convert the regional grid into a three-dimensional grid based on the depth gradient and determine the depth attenuation function of each coordinate point in the three-dimensional grid;

[0048] A radon concentration inversion unit is used to determine the standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value and the depth attenuation function of each coordinate point;

[0049] The concentration visualization unit is used to generate horizontal slice diagrams, vertical section diagrams and three-dimensional rendering diagrams according to the standard radon concentration value of each coordinate point in the three-dimensional grid.

[0050] Compared with the prior art, the present invention has the following beneficial effects:

[0051] The present invention introduces a comprehensive correction factor, comprehensively considering the impact of multiple factors such as geological structure, climatic conditions, and covering vegetation type on radon concentration, and corrects the collected original radon concentration values, providing a reliable data basis for subsequent regional radon concentration analysis. The method also calculates the corrected radon concentration value for each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value, fully considering the spatial relationship between the measurement points and the accuracy of the corrected data, and can more accurately reflect the spatial distribution characteristics of the regional radon concentration. In addition, the standard radon concentration value for each coordinate point in the regional grid is finally inverted based on the corrected radon concentration value, which helps to more comprehensively and accurately understand the distribution pattern and change trend of regional soil radon concentration. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] Figure 1 The figure is a flow chart of a method for processing soil radon concentration data.

[0053] Figure 2 This is a flow chart for calculating the comprehensive correction factor in a soil radon concentration data processing method.

[0054] Figure 3 The present invention is a flow chart for determining the corrected radon concentration value in a soil radon concentration data processing method.

[0055] Figure 4 The present invention is a flowchart for calculating the corrected radon concentration value of each coordinate point in a soil radon concentration data processing method.

[0056] Figure 5 This is a flow chart for inverting the standard radon concentration value of each coordinate point in a soil radon concentration data processing method.

[0057] Figure 6This is a structural diagram of a soil radon concentration data processing system. DETAILED DESCRIPTION

[0058] In order to make the purpose, technical solutions and advantages of the present invention clearer, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.

[0059] The specific implementation of the present invention is described in detail below with reference to specific embodiments.

[0060] like Figure 1 As shown, an embodiment of the present invention provides a method for processing soil radon concentration data, the method comprising the following steps:

[0061] S100, receiving radon concentration collection information, the radon concentration collection information including radon concentration values, each radon concentration value Ri corresponding to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, pressure Pi, humidity Hi and covering vegetation type Vi;

[0062] S200, calculating the comprehensive correction factor Fi, F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), where fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor;

[0063] S300, determining a corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function;

[0064] S400, calculating a corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi;

[0065] S500: Retrieve a standard topographic map, which includes standard temperature, standard pressure, standard humidity, and geological structure and vegetation type at each location, and invert to obtain a standard radon concentration value at each coordinate point in the regional grid.

[0066] It should be noted that the original collected soil radon concentration data is affected by a variety of environmental factors, such as geological structure, temperature, air pressure, humidity, and covering vegetation type. Differences in the porosity, permeability and other characteristics of rocks in different geological structures will lead to different migration and enrichment patterns of radon; changes in temperature, air pressure and humidity will affect the diffusion and decay process of radon; and the covering vegetation type may indirectly affect the radon concentration by affecting the air permeability and moisture content of the soil. Due to the complex and changeable actual measurement environmental conditions, it is difficult to directly compare and analyze measurement data from different times and places. There is currently a lack of methods to convert actual measurement data to standard environmental conditions. In addition, since the distribution of measurement points in the actual measurement process is often uneven, it is difficult to directly obtain continuous radon concentration distribution information for the entire area. The embodiments of the present invention are intended to solve the above problems.

[0067] In an embodiment of the present invention, the obtained radon concentration collection information needs to be uploaded first. The radon concentration collection information includes a large number of radon concentration values, each of which corresponds to collection coordinates, collection depth, geological structure, temperature, air pressure, humidity, and covering vegetation type. The collection coordinates here refer to plane coordinates. Then, a comprehensive correction factor Fi corresponding to each radon concentration value is calculated. The comprehensive correction factor is obtained by integrating the geological correction factor, the climate correction factor, and the vegetation correction factor. By introducing the comprehensive correction factor, the embodiment of the present invention comprehensively considers the influence of multiple factors on radon concentration, such as geological structure, climate conditions (temperature, air pressure, humidity), and covering vegetation type, and corrects the collected original radon concentration value. This can effectively eliminate the interference of these factors, obtain more accurate corrected radon concentration values, and provide a reliable data basis for subsequent regional radon concentration analysis. Next, a depth attenuation function is constructed, and the corrected radon concentration value Rxi is determined in combination with the comprehensive correction factor Fi. Then, based on the inverse distance weighted method and the corrected radon concentration value Rxi, the corrected radon concentration value is calculated for each coordinate point in the regional grid. This fully considers the spatial relationship between the measurement points and the accuracy of the corrected data, can more accurately reflect the spatial distribution characteristics of the regional radon concentration, improve the accuracy of spatial interpolation, and provide more reliable results for drawing regional soil radon concentration distribution maps. In addition, the embodiment of the present invention also requires the prior determination of a standard topographic map that covers the entire study area and includes standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location. In this way, the standard radon concentration value for each coordinate point in the regional grid can be inverted based on the corrected radon concentration value, allowing soil radon concentration data measured under different environmental conditions and at different times to be compared and analyzed under a unified standard environment. This provides a scientific basis for the evaluation and comparison of regional soil radon concentrations and helps to more comprehensively and accurately understand the distribution pattern and changing trend of regional soil radon concentrations.

[0068] like Figure 2 As shown, as a preferred embodiment of the present invention, the step of calculating the comprehensive correction factor Fi specifically includes:

[0069] S201, calculate the geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively;

[0070] S202, calculate the climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient;

[0071] S203, calculating the vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), where δ is the vegetation suppression coefficient;

[0072] S204: Calculate a comprehensive correction factor Fi based on the geological correction factor, the climate correction factor, and the vegetation correction factor.

[0073] In the embodiment of the present invention, Gi needs to be digitized according to the geological structure type. Gi of granite is 1.5, Gi of sedimentary rock is 1.0, and Gi of clay layer is 0.8. Vi needs to be digitized according to the vegetation type. Vi of forest is 1.0, Vi of grassland is 0.9, and Vi of bare soil is 0.7. Then the geological correction factor fgeology(Gi) can be calculated. fgeology α is the geological influence weight coefficient, α=0.3, μ G and σ G are the mean and standard deviation of geological structures, for example, μ G and σ G 1 and 0.5 respectively. Then calculate the climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, set to 0.001, and γ is the humidity attenuation coefficient, set to 0.02. Rising temperature and falling pressure promote radon exhalation, while increasing humidity inhibits radon migration. Next, calculate the vegetation correction factor, fvegetation(Vi), as fvegetation(Vi) = 1-δ × ln(Vi + 1), where δ is the vegetation inhibition coefficient, set to 0.1. Finally, multiply the geological correction factor, climate correction factor, and vegetation correction factor to obtain the combined correction factor, Fi.

[0074] like Figure 3 As shown in FIG. 1 , as a preferred embodiment of the present invention, the step of determining the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function specifically includes:

[0075] S301, determining a depth attenuation function;

[0076] S302: Calculate and obtain the corrected radon concentration value Rxi.

[0077] In the embodiment of the present invention, the depth attenuation function is determined as w(di)=exp(-λdi), where λ is the radon diffusion attenuation coefficient, λ=0.05m -1 Then, the corrected radon concentration value Rxi can be calculated, Rxi = Ri × Fi × w(di). This allows soil radon concentration data measured under different environmental conditions to be mapped to a unified standard environment.

[0078] like Figure 4 As shown in FIG. 1 , as a preferred embodiment of the present invention, the step of calculating the corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi specifically includes:

[0079] S401, generating evenly distributed regional grid points according to the regional range and grid resolution, and determining each coordinate point in the regional grid;

[0080] S402, determining the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi;

[0081] S403, determining the corrected radon concentration value of the corresponding coordinate point according to the collected coordinates within the search radius and Rxi.

[0082] In the embodiment of the present invention, evenly distributed regional grid points are generated based on the regional range and grid resolution, with a spacing of 10 meters between grid points. Each coordinate point in the regional grid is determined, and the coordinate point is a plane coordinate. Then, the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi are determined. The search radius needs to be determined in advance. Finally, the corrected radon concentration value of the corresponding coordinate point can be determined based on the acquisition coordinates within the search radius and Rxi. where d ij is the horizontal distance from collection coordinate i to grid coordinate point j. p is a power parameter that controls how quickly the weight decays with distance. N represents the number of collection coordinates within the search radius. A larger p value increases the weight of neighboring points, resulting in a more localized interpolation result (potentially less smooth). A smaller p value results in a more uniform distribution, but may lead to oversmoothing. This can be optimized through cross-validation. A value of 2 or 3 is generally recommended.

[0083] like Figure 5 As shown in FIG. 1 , as a preferred embodiment of the present invention, the step of inverting to obtain the standard radon concentration value of each coordinate point in the regional grid specifically includes:

[0084] S501, obtaining a comprehensive correction factor for each coordinate point based on the standard temperature, standard pressure, standard humidity, and the geological structure and vegetation type of each location;

[0085] S502, converting the regional grid into a three-dimensional grid based on the depth gradient, and determining a depth attenuation function for each coordinate point in the three-dimensional grid;

[0086] S503, determining a standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value, and the depth attenuation function of each coordinate point;

[0087] S504: Generate a horizontal slice diagram, a vertical section diagram, and a three-dimensional volume rendering diagram according to the standard radon concentration value of each coordinate point in the three-dimensional grid.

[0088] In an embodiment of the present invention, when performing inversion, a comprehensive correction factor is obtained for each coordinate point based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type at each location. The regional grid is then converted into a three-dimensional grid based on the depth gradient, and the depth attenuation function is determined for each coordinate point in the three-dimensional grid. The depth gradient needs to be determined in advance, for example, the depth gradient includes 1 meter, 3 meters, 5 meters, and 7 meters. The standard radon concentration value can then be determined based on the comprehensive correction factor, the corrected radon concentration value, and the depth attenuation function for each coordinate point. Horizontal slices, vertical profiles, and three-dimensional renderings are generated based on these standard radon concentration values, resulting in better visualization and facilitating the study of the distribution patterns and changing trends of regional soil radon concentrations.

[0089] like Figure 6 As shown, an embodiment of the present invention further provides a soil radon concentration data processing system, the system comprising:

[0090] The collected information uploading module 100 is used to receive radon concentration collected information, wherein the radon concentration collected information includes radon concentration values, each radon concentration value Ri corresponds to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, pressure Pi, humidity Hi and covering vegetation type Vi;

[0091] The comprehensive correction factor module 200 is used to calculate the comprehensive correction factor Fi, where F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), where fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor;

[0092] The radon concentration correction module 300 is used to determine the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function;

[0093] The modified concentration calculation module 400 is used to calculate the modified radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the modified radon concentration value Rxi;

[0094] The standard concentration calculation module 500 is used to retrieve a standard topographic map containing standard temperature, standard pressure, standard humidity, and geological structure and vegetation type at each location, and invert to obtain the standard radon concentration value of each coordinate point in the regional grid.

[0095] As a preferred embodiment of the present invention, the comprehensive correction factor module 200 includes:

[0096] Geological correction factor unit, used to calculate geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively;

[0097] Climate correction factor unit, used to calculate climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient;

[0098] Vegetation correction factor unit, used to calculate vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), δ is vegetation suppression coefficient;

[0099] The comprehensive correction factor unit is used to calculate the comprehensive correction factor Fi based on the geological correction factor, climate correction factor and vegetation correction factor.

[0100] As a preferred embodiment of the present invention, the radon concentration correction module 300 includes:

[0101] A depth attenuation function unit is used to determine a depth attenuation function w(di)=exp(-λdi), where λ is a radon diffusion attenuation coefficient;

[0102] The radon concentration value calculation unit is used to calculate the corrected radon concentration value Rxi, Rxi = Ri×Fi×w(di).

[0103] As a preferred embodiment of the present invention, the corrected concentration calculation module 400 includes:

[0104] A regional grid determination unit is used to generate evenly distributed regional grid points according to the regional range and grid resolution, and determine each coordinate point in the regional grid;

[0105] An acquisition coordinate search unit, configured to determine the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi;

[0106] The corrected concentration calculation unit is used to determine the corrected radon concentration value of the corresponding coordinate point according to the collected coordinates within the search radius and Rxi. where d ij is the horizontal distance from the collection coordinate to the grid coordinate point, p is the power parameter, which is used to control the speed at which the weight decays with distance, and N represents the number of collection coordinates within the search radius.

[0107] As a preferred embodiment of the present invention, the standard concentration calculation module 500 includes:

[0108] The coordinate point correction factor unit is used to obtain the comprehensive correction factor of each coordinate point based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location;

[0109] A coordinate point depth attenuation unit is used to convert the regional grid into a three-dimensional grid based on the depth gradient and determine the depth attenuation function of each coordinate point in the three-dimensional grid;

[0110] A radon concentration inversion unit is used to determine the standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value and the depth attenuation function of each coordinate point;

[0111] The concentration visualization unit is used to generate horizontal slice diagrams, vertical section diagrams and three-dimensional rendering diagrams according to the standard radon concentration value of each coordinate point in the three-dimensional grid.

[0112] The above is only a detailed description of the preferred embodiments of the present invention, which is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

[0113] It should be understood that, although the various steps in the flow chart of each embodiment of the present invention are shown in sequence according to the indication of the arrows, these steps are not necessarily performed in sequence according to the order indicated by the arrows. Unless otherwise specified herein, the execution of these steps is not strictly limited in order, and these steps can be performed in other orders. Moreover, at least a portion of the steps in each embodiment may include a plurality of sub-steps or a plurality of stages, and these sub-steps or stages are not necessarily performed at the same time, but can be performed at different times, and the execution order of these sub-steps or stages is not necessarily performed in sequence, but can be performed in turn or alternately with at least a portion of other steps or sub-steps or stages of other steps.

[0114] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0115] Those skilled in the art will readily appreciate other embodiments of the present disclosure after considering the disclosure in the specification and examples. This application is intended to cover any variations, uses, or adaptations of the present disclosure that follow the general principles of the present disclosure and include common knowledge or customary techniques in the art not disclosed herein. The description and examples are to be considered merely as exemplary, and the true scope and spirit of the present disclosure are indicated by the claims.

Claims

1. A soil radon concentration data processing method, characterized in that: The method comprises the following steps: Receive radon concentration collection information, the radon concentration collection information including radon concentration values, each radon concentration value Ri corresponding to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, pressure Pi, humidity Hi and covering vegetation type Vi; Calculate the comprehensive correction factor Fi, F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor; Determine the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function; Calculate the corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi; A standard topographic map is retrieved, which includes standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location, and the standard radon concentration value of each coordinate point in the regional grid is obtained by inversion.

2. The soil radon concentration data processing method according to claim 1, characterized in that: The step of calculating the comprehensive correction factor Fi specifically includes: Calculate the geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively; Calculate the climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient; Calculate the vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), δ is the vegetation suppression coefficient; The comprehensive correction factor Fi is calculated based on the geological correction factor, climate correction factor and vegetation correction factor.

3. The soil radon concentration data processing method according to claim 1, characterized in that: The step of determining the corrected radon concentration value Rxi based on the comprehensive correction factor Fi and the depth attenuation function specifically includes: Determine the depth attenuation function w(di)=exp(-λdi), where λ is the radon diffusion attenuation coefficient; The corrected radon concentration value Rxi is calculated, Rxi=Ri×Fi×w(di).

4. The soil radon concentration data processing method according to claim 1, characterized in that: The step of calculating the corrected radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the corrected radon concentration value Rxi specifically includes: Generate evenly distributed regional grid points according to the regional range and grid resolution, and determine each coordinate point in the regional grid; Determine the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi; According to the collected coordinates and Rxi within the search radius, the corrected radon concentration value of the corresponding coordinate point is determined. where d ij is the horizontal distance from the collection coordinate to the grid coordinate point, p is the power parameter, which is used to control the speed at which the weight decays with distance, and N represents the number of collection coordinates within the search radius.

5. The soil radon concentration data processing method according to claim 1, characterized in that: The step of inverting to obtain the standard radon concentration value of each coordinate point in the regional grid specifically includes: The comprehensive correction factor of each coordinate point is obtained based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location; Convert the regional grid into a three-dimensional grid based on the depth gradient, and determine the depth attenuation function of each coordinate point in the three-dimensional grid; Determine the standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value and the depth attenuation function of each coordinate point; Horizontal slices, vertical sections and three-dimensional renderings are generated based on the standard radon concentration values at each coordinate point in the three-dimensional grid.

6. A soil radon concentration data processing system, characterized in that: The system comprises: The collected information uploading module is used to receive radon concentration collected information, wherein the radon concentration collected information includes radon concentration values, and each radon concentration value Ri corresponds to the collection coordinates, collection depth di, geological structure Gi, temperature Ti, air pressure Pi, humidity Hi and covering vegetation type Vi; Comprehensive correction factor module, used to calculate the comprehensive correction factor Fi, F = fgeology (Gi) × fclimate (Ti, Pi, Hi) × fvegetation (Vi), fgeology (Gi) is the geological correction factor, fclimate (Ti, Pi, Hi) is the climate correction factor, and fvegetation (Vi) is the vegetation correction factor; A radon concentration correction module is used to determine a corrected radon concentration value Rxi based on a comprehensive correction factor Fi and a depth attenuation function; A modified concentration calculation module is used to calculate the modified radon concentration value of each coordinate point in the regional grid based on the inverse distance weighted method and the modified radon concentration value Rxi; The standard concentration calculation module is used to retrieve the standard topographic map, which contains standard temperature, standard pressure, standard humidity, and geological structure and covering vegetation type of each location, and invert to obtain the standard radon concentration value of each coordinate point in the regional grid.

7. The soil radon concentration data processing system according to claim 6, characterized in that: The comprehensive correction factor module includes: Geological correction factor unit, used to calculate geological correction factor fgeology (Gi), fgeology α is the geological influence weight coefficient, μ G and σ G are the mean and standard deviation of geological structures, respectively; Climate correction factor unit, used to calculate climate correction factor fclimate (Ti, Pi, Hi), fclimate β is the temperature-pressure coupling coefficient, γ is the humidity attenuation coefficient; Vegetation correction factor unit, used to calculate vegetation correction factor fvegetation(Vi), fvegetation(Vi)=1-δ×ln(Vi+1), δ is vegetation suppression coefficient; The comprehensive correction factor unit is used to calculate the comprehensive correction factor Fi based on the geological correction factor, climate correction factor and vegetation correction factor.

8. The soil radon concentration data processing system according to claim 6, characterized in that: The radon concentration correction module includes: A depth attenuation function unit is used to determine a depth attenuation function w(di)=exp(-λdi), where λ is a radon diffusion attenuation coefficient; The radon concentration value calculation unit is used to calculate the corrected radon concentration value Rxi, Rxi = Ri×Fi×w(di).

9. The soil radon concentration data processing system according to claim 6, characterized in that: The corrected concentration calculation module includes: A regional grid determination unit is used to generate evenly distributed regional grid points according to the regional range and grid resolution, and determine each coordinate point in the regional grid; An acquisition coordinate search unit, configured to determine the acquisition coordinates of each coordinate point within the search radius and the corresponding corrected radon concentration value Rxi; The corrected concentration calculation unit is used to determine the corrected radon concentration value of the corresponding coordinate point according to the collected coordinates within the search radius and Rxi. where d ij is the horizontal distance from the collection coordinate to the grid coordinate point, p is the power parameter, which is used to control the speed at which the weight decays with distance, and N represents the number of collection coordinates within the search radius.

10. The soil radon concentration data processing system according to claim 6, characterized in that: The standard concentration calculation module includes: The coordinate point correction factor unit is used to obtain the comprehensive correction factor of each coordinate point based on the standard temperature, standard pressure, standard humidity, and the geological structure and covering vegetation type of each location; A coordinate point depth attenuation unit is used to convert the regional grid into a three-dimensional grid based on the depth gradient and determine the depth attenuation function of each coordinate point in the three-dimensional grid; A radon concentration inversion unit is used to determine the standard radon concentration value based on the comprehensive correction factor, the corrected radon concentration value and the depth attenuation function of each coordinate point; The concentration visualization unit is used to generate horizontal slice diagrams, vertical section diagrams and three-dimensional rendering diagrams according to the standard radon concentration value of each coordinate point in the three-dimensional grid.

Images