A calculation method and system for the change amount of the reserves of geochemical elements

Abstract

This application relates to the field of geochemical technology, and specifically discloses a calculation method and system for the change amount of geochemical element reserves, which includes: calculating and obtaining a first set based on the first sample collection result of the sampling area; combining the first set and the spatial position relationship of the sampling points to obtain a first prediction model of the geochemical element density value in the sampling area, and obtaining the first reserve data of the geochemical elements; calculating and obtaining a second set based on the second sample collection result of the sampling area; combining the second set and the spatial position relationship of the sampling points to obtain a second prediction model of the geochemical element density value in the sampling area, and obtaining the second reserve data of the geochemical elements; combining the first reserve data and the second reserve data to obtain the change amount data of the geochemical elements in the sampling area. Through this application, the spatial heterogeneous factors of the sampling points can be taken into account in the calculation process of the geochemical element reserves, greatly improving the accuracy of the calculation results, and at the same time, the error range of the results can be given.

Description

Technical Field

[0001] This application relates to the field of geochemical technologies, and in particular, to a method and system for calculating the change amount of geochemical element reserves. Background Art

[0002] Existing methods for measuring geochemical element reserves obtain the content values of geochemical elements at sampling points through sampling at sampling points and laboratory analysis, and then calculate the density values of geochemical elements at the sampling points. Then, based on the density of geochemical elements at the sampling points, the geochemical element reserves in the region are estimated. However, due to the spatial heterogeneity of the relevant parameters of the sampling points, the calculation results of the geochemical element reserves obtained by the above method have relatively large errors at this time. Secondly, this method cannot evaluate the error of the final calculation result of the geochemical element reserves. In the case of a relatively large real error, it will also misjudge the change amount of the chemical element reserves. Summary of the Invention

[0003] In order to solve the problem that the spatial heterogeneity of the relevant parameters of the sampling points leads to relatively large errors in the calculation results of the change amount of existing geochemical element reserves, this application provides a method and system for calculating the change amount of geochemical element reserves.

[0004] According to one aspect of this application, a method for calculating the change amount of geochemical element reserves is provided, including: calculating and obtaining a first set based on the first sample collection result of the sampling area, where the first set includes density values of geochemical element reserves corresponding to multiple sampling points;

[0005] Combining the first set and the spatial position relationship of the sampling points to obtain a first prediction model of the geochemical element density value in the sampling area, and obtaining the first reserve data of the geochemical element based on the first prediction model and the sampling area;

[0006] Calculating and obtaining a second set based on the second sample collection result of the sampling area, where the second set includes density values of geochemical element reserves corresponding to multiple sampling points;

[0007] Combining the second set and the spatial position relationship of the sampling points to obtain a second prediction model of the geochemical element density value in the sampling area, and obtaining the second reserve data of the geochemical element based on the second prediction model and the sampling area;

[0008] Combining the first reserve data and the second reserve data to obtain the change amount data of the geochemical element reserves in the sampling area.

[0009] By adopting the above technical solution, first, a first set including density values of geochemical element reserves corresponding to multiple sampling points is obtained based on the first sample collection result of the sampling area. A first prediction model of the geochemical element density value of the sampling area is obtained by combining the first set and the spatial position relationship of the sampling points, and a first reserve data of the geochemical element is obtained based on the first prediction model and the sampling area. The second set is calculated and obtained based on the second sample collection result of the sampling area, and the second set includes density values of geochemical element reserves corresponding to multiple sampling points. A second prediction model of the geochemical element density value of the sampling area is obtained by combining the second set and the spatial position relationship of the sampling points, and a second reserve data of the geochemical element is obtained based on the second prediction model and the sampling area. The change amount of the reserves of the geochemical element in the sampling area is obtained by combining the first reserve data and the second reserve data. Compared with the classical statistical model that does not consider the spatial heterogeneity of the relevant parameters of the sampling points, the prediction model obtained by combining the spatial position relationship of the sampling points is more accurate.

[0010] Optionally, the obtaining of the first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points includes:

[0011] Calculating the empirical semivariogram of the density value of the geochemical element reserve based on the first set and the spatial position relationship of the sampling points;

[0012] Using an approximation criterion to screen the target model to fit with the empirical semivariogram to obtain the theoretical semivariogram;

[0013] Based on the theoretical semivariogram, the first prediction model of the geochemical element density value of the sampling area is obtained by using the block kriging method.

[0014] By adopting the above technical solution, due to the different spatial distributions of multiple sampling points, the empirical semivariogram of the density value of the geochemical element reserve is obtained by correlating the density value of the geochemical element reserve of the sampling points with the spatial position relationship of the sampling points. Then, the target model is screened to fit and generate the theoretical semivariogram. Finally, the block kriging method is used for analysis to obtain the first prediction model of the geochemical element density value of the sampling area with higher accuracy.

[0015] Preferably, the first reserve data includes the first reserve value of the geochemical element in the sampling area and the corresponding first standard error. The obtaining of the first reserve data of the geochemical element based on the first prediction model and the sampling area includes:

[0016] Obtain the first Kriging prediction value and the first variance of the geochemical element density value in the sampling area based on the first prediction model; obtain the first reserve value according to the Kriging prediction value of the geochemical element density value and the area of the sampling area, and obtain the first standard error by using the first variance of the geochemical element density value and the area of the sampling area.

[0017] By adopting the above technical solution, the spatial position relationship of the sampling area is incorporated into the calculation factors, and the calculation using the first prediction model can intuitively and accurately obtain the first reserve value and the first standard error of the geochemical elements in the sampling area, which is beneficial to the intuitive and accurate calculation results.

[0018] Preferably, the second prediction model for obtaining the geochemical element density value in the sampling area by combining the spatial position relationship of the second set and the sampling points includes:

[0019] Calculate the empirical semi-variogram of the geochemical element reserve density value based on the spatial position relationship of the second set and the sampling points;

[0020] Use the approximation criterion to screen the target model to fit with the empirical semi-variogram to obtain the theoretical semi-variogram;

[0021] Based on the theoretical semi-variogram, use the Kriging method to obtain the second prediction model of the geochemical element density value in the sampling area.

[0022] By adopting the above technical solution, due to the different spatial distributions of multiple sampling points, the empirical semi-variogram of the geochemical element reserve density value is obtained by correlating the reserve density value of the geochemical elements at the sampling points with the spatial position relationship of the sampling points, then the target model is screened to fit and generate the theoretical semi-variogram, and finally the Kriging method is used for analysis to obtain the second prediction model of the geochemical element density value in the sampling area with higher accuracy.

[0023] Preferably, the second reserve data includes the second reserve value of the geochemical elements in the sampling area and the corresponding second standard error. The obtaining of the second reserve data of the geochemical elements based on the second prediction model and the area of the sampling area includes:

[0024] Obtain the second Kriging prediction value and the second variance of the geochemical element density value in the sampling area based on the second prediction model; obtain the second reserve value according to the Kriging prediction value of the geochemical element density value and the area of the sampling area, and obtain the second standard error by using the second variance of the geochemical element density value and the area of the sampling area.

[0025] By adopting the above technical solution, the spatial position relationship of the sampling area is incorporated into the calculation factors, and through the calculation of the first prediction model, the second reserve value and the second standard error of the geochemical elements in the sampling area can be intuitively and accurately obtained, which is beneficial to the intuitive and accurate calculation results.

[0026] Preferably, the obtaining of the reserve change amount data of the geochemical elements in the sampling area by combining the first reserve data and the second reserve data includes:

[0027] Taking the difference between the first reserve value and the second reserve value to obtain the reserve change amount of the geochemical elements in the sampling area; calculating the first standard error and the second standard error to obtain the standard error of the reserve change amount;

[0028] Determining the confidence interval of the reserve change amount of the geochemical elements at a preset significance level according to the reserve change amount of the geochemical elements and the standard error of the reserve change amount.

[0029] By adopting the above technical solution, taking the difference between the first reserve value and the second reserve value to obtain the reserve change amount of the geochemical elements in the sampling area, calculating the first standard error and the second standard error to obtain the standard error of the reserve change amount, the reserve change amount and the corresponding standard error of the geochemical elements in the sampling area can be obtained. By determining the confidence interval of the reserve change amount of the geochemical elements at a preset significance level, the accuracy of the calculation result of the reserve change amount can be further improved.

[0030] Preferably, the first set includes a first training sample set, and the first prediction model is generated through the following training:

[0031] Generating the first training sample set, where the first training sample set includes a preset number of sampling points and the density values of the geochemical element reserves corresponding to the sampling points;

[0032] Training a preset model with the samples in the first training sample set, using the sampling points as inputs and the virtual density values of the geochemical element reserves corresponding to the sampling points as outputs. When the unification rate of the output virtual density values and the density values of the geochemical element reserves meets the preset threshold, the training of the preset model is completed and the first prediction model is generated.

[0033] By adopting the above technical solution, taking the density value of the geochemical element reserves obtained from the sampling area collection results as the training extrusion of the deep learning model, training to obtain the virtual density value of the geochemical element reserves corresponding to the sample points, and judging the authenticity of the virtual density value according to the density value of the geochemical element reserves, and then feeding back the judgment result to the deep learning model, so that the virtual density value generated by the deep learning model subsequently is closer to the density value of the geochemical element reserves. When the unification rate of the output virtual density value and the density value of the geochemical element reserves meets the preset threshold, the training of the preset model is completed and the first prediction model is generated.

[0034] Preferably, before obtaining the first reserve data of the geochemical element based on the first prediction model and the sampling area area, it includes:

[0035] Judging whether the sampling quantity of the sampling points exceeds a preset quantity threshold;

[0036] If the sampling quantity exceeds the quantity threshold, then obtain the mean value and standard error of the density value of the geochemical element reserves in the sampling area according to the first set;

[0037] If the sampling quantity does not exceed the quantity threshold, then process the first set by interpolation method to obtain a prediction curve model;

[0038] Obtain the mean value and standard error of the density value of the geochemical element reserves in the sampling area according to the prediction curve model.

[0039] By adopting the above technical solution, when the sampling quantity of the sampling points is small, it may lead to inaccurate prediction results. Therefore, the sampling quantity is judged by a preset quantity threshold. When the sampling quantity does not exceed the quantity threshold, the interpolation method can be used to process and obtain a prediction curve with higher prediction accuracy to calculate the mean value and standard error of the density value of the geochemical element reserves.

[0040] Preferably, after calculating and obtaining the first set based on the first sample collection result of the sampling area, it further includes:

[0041] Judging whether the density value of the geochemical element reserves in the first set is abnormal. The abnormality includes whether the density value of the geochemical element reserves corresponding to the sampling points exceeds the preset range. If the density value is abnormal, a reminder for re-sampling the corresponding sampling points will appear.

[0042] By adopting the above technical solution, since the density value of the geochemical element reserves is usually collected by artificial sampling, it is judged by presetting the range of the density value of the geochemical element reserves, so as to avoid the interference of obvious abnormal density values on the collection results.

[0043] According to another aspect of the present application, there is also provided a calculation system for the change amount of geochemical element reserves, including:

[0044] A first acquisition module, configured to calculate and obtain a first set based on the first sample acquisition result of the sampling area, where the first set includes density values of geochemical element reserves corresponding to multiple sampling points;

[0045] A first model processing module, configured to obtain a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain first reserve data of geochemical elements based on the first prediction model and the area of the sampling area;

[0046] A second acquisition module, configured to calculate and obtain a second set based on the second sample acquisition result of the sampling area, where the second set includes density values of geochemical element reserves corresponding to multiple sampling points;

[0047] A second model processing module, configured to obtain a second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain second reserve data of geochemical elements based on the second prediction model and the area of the sampling area;

[0048] A calculation module, configured to obtain the change amount data of the geochemical element reserves in the sampling area by combining the first reserve data and the second reserve data.

[0049] In summary, the present application includes the following beneficial technical effects:

[0050] Calculate and obtain a first set based on the first sample acquisition result of the sampling area; obtain a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain first reserve data of geochemical elements; calculate and obtain a second set based on the second sample acquisition result of the sampling area; obtain a second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain second reserve data of geochemical elements; obtain the change amount of the geochemical element reserves in the sampling area by combining the first reserve data and the second reserve data. Through the present application, the spatial heterogeneous factors of the sampling points can be considered in the calculation process of the geochemical element reserves, greatly improving the accuracy of the calculation results and giving the error range of the results. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 Shows a schematic diagram of the scenario of sampling points in a sampling area of the present application.

[0052] Figure 2 Shows a schematic flow diagram of a calculation method for the change amount of geochemical element reserves of the present application.

[0053] Figure 3 Shows a schematic flow chart of obtaining a first prediction model provided by an embodiment of the present application.

[0054] Figure 4 Shows a schematic flow chart of obtaining first reserve data provided by an embodiment of the present application.

[0055] Figure 5 Shows another schematic flow chart of obtaining a first prediction model provided by an embodiment of the present application.

[0056] Figure 6 Shows a schematic flow chart of obtaining the mean and standard error of the reserve density value of geochemical elements in the present application.

[0057] Figure 7 Shows a schematic structural diagram of a calculation system for the change amount of geochemical element reserves provided by an embodiment of the present application. Detailed implementation manners

[0058] The following makes the purpose, technical solutions and advantages of the present application clearer. The following further details the present application through the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0059] The following combines the attached Figure 1-7 Further details the present application.

[0060] Figure 1 Shows a schematic diagram of the scenario of sampling points in a sampling area. Usually, in the method for calculating the reserves of geochemical elements, multiple sampling points are set in the sampling area for sampling detection and analysis of geochemical elements, and the reserves of geochemical elements are calculated based on the sample sampling results of the sampling points.

[0061] Figure 2 Shows a schematic flow chart of a calculation method for the change amount of geochemical element reserves in the present application. It can be seen from Figure 2 that the calculation method for the change amount of geochemical element reserves in this embodiment includes the following steps:

[0062] S102, calculating and obtaining a first set based on the first sample collection result of the sampling area, where the first set includes density values of reserves of geochemical elements corresponding to multiple sampling points.

[0063] Further, the density value of the reserves of geochemical elements at the sampling point can be obtained through laboratory measurement and analysis of the samples obtained from the sampling points in the sampling area.

[0064] In one embodiment, the sampling area includes a soil sampling area, the geochemical element includes carbon element, and the carbon element reserve density value of the sampling points in the soil sampling area can be obtained by the above method. Here, the first sample collection result is the measurement and analysis result of the sampling points in the sampling area. For example, for sampling point A, the density value of the corresponding geochemical element reserve is B, then (A, B) is the sample data of a sampling point in the first sample sampling result.

[0065] S104. Obtain a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain the first reserve data of the geochemical element based on the first prediction model and the sampling area area.

[0066] It can be understood that by combining the first set and the position relationship of the sampling points to obtain the first prediction model of the geochemical element density value, compared with the classical statistical model based only on the geochemical element content density value of the sampling points, the prediction result obtained by combining the spatial position relationship factor of the sampling points will be more accurate.

[0067] Figure 3 The flowchart showing the process of obtaining the first prediction model provided by the embodiment of the present application is shown.

[0068] As Figure 3 shown, in some embodiments, the obtaining a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points includes:

[0069] S202. Calculate the empirical semivariogram of the geochemical element reserve density value based on the first set and the spatial position relationship of the sampling points.

[0070] Based on the first set and the spatial position relationship of the sampling points, according to the mathematical modeling method of the empirical semivariogram, it can well describe the spatial continuous variability of the sampling points in the process of geochemical element measurement, and reflect the change between the differences in soil properties and the carbon content density value. Among them, to create an empirical semivariogram, determine the squared difference of all sampling point position pairs, and plot these position pairs (the X-axis coordinate is the position spacing, and the Y-axis coordinate is half of the squared difference). After forming a semivariogram cloud, the semivariogram cloud can be used to explore and quantify the spatial dependence of the reserve density value of the geochemical element and the sampling points, and then quantify the hypothesis that things closer are more similar. Generally, the empirical semivariogram model is usually not affected by the random sampling of sampling points because all data are used to generate predicted values, which is beneficial to improving the accuracy of the prediction result.

[0071] S204. Use an approximation criterion to screen a target model to fit with the empirical semivariogram to obtain a theoretical semivariogram.

[0072] Among them, by using an approximation criterion to screen the target model, a theoretical semivariogram is obtained by fitting the empirical semivariogram, and then a model with higher prediction accuracy is obtained.

[0073] S206. Based on the theoretical semivariogram, a first prediction model of the geochemical element density value in the sampling area is obtained by using the block kriging method.

[0074] Among them, the kriging method is based on the theory of variograms and structural analysis, and is a method for linearly unbiased optimal estimation of the values of regionalized scalars in a limited area. It has a stable prediction effect on the reserve density value of geochemical elements. Based on the theoretical semivariogram, obtaining a first prediction model of the geochemical element density value in the sampling area by using the block kriging method can improve the prediction accuracy.

[0075] Figure 4 The flowchart showing the process of obtaining the first reserve data provided by the embodiment of the present application is shown.

[0076] Further, as Figure 4 shown, the first reserve data includes the first reserve value of the geochemical elements in the sampling area and the corresponding first standard error. Obtaining the first reserve data of the geochemical elements based on the first prediction model and the area of the sampling area includes:

[0077] S302. Based on the first prediction model, obtain the first block kriging prediction value and the first variance of the geochemical element density value in the sampling area.

[0078] S304. Obtain the first reserve value according to the block kriging prediction value of the geochemical element density value and the area of the sampling area, and obtain the first standard error by using the first variance of the geochemical element density value and the area of the sampling area.

[0079] Among them, the calculation formula for the reserve value of the geochemical element sampling area is:

[0080] SCS = Mean × S

[0081] In the formula for the reserve value of the sampling area, SCS is the reserve value of the sampling area, Mean is the kriging prediction value (i.e., the mean value of the density value of the geochemical element reserves), and S is the regional area of the sampling area.

[0082] The calculation formula for the standard error of the geochemical element is:

[0083]

[0084] Var is the variance of the geochemical element density value. The formulas for the reserve value and its standard error of the geochemical element sampling area are: SCS ± SE SCS。

[0085] S106. Calculate and obtain a second set based on the second sample collection result of the sampling area, where the second set includes density values of geochemical element reserves corresponding to multiple sampling points.

[0086] S108. Obtain a second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain second reserve data of the geochemical elements based on the second prediction model and the area of the sampling area.

[0087] Further, the obtaining of the second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points includes:

[0088] Calculate the empirical semivariogram of the geochemical element reserve density value based on the second set and the spatial position relationship of the sampling points;

[0089] Use an approximation criterion to screen a target model to fit with the empirical semivariogram to obtain a theoretical semivariogram;

[0090] Based on the theoretical semivariogram, obtain the second prediction model of the geochemical element density value of the sampling area by using the block kriging method.

[0091] Further, the second reserve data includes the second reserve value of the geochemical elements in the sampling area and the corresponding second standard error. The obtaining of the second reserve data of the geochemical elements based on the second prediction model and the area of the sampling area includes:

[0092] Obtain the second block kriging prediction value and the second variance of the geochemical element density value of the sampling area based on the second prediction model; obtain the second reserve value according to the block kriging prediction value of the geochemical element density value and the area of the sampling area, and obtain the second standard error by using the second variance of the geochemical element density value and the area of the sampling area.

[0093] Among them, the obtaining method and principle of the second reserve data are the same as those of the first reserve data. For details, refer to the embodiments of the first reserve data obtaining, which will not be elaborated here.

[0094] Further, since both the first prediction model and the second prediction model are prediction models for the reserve value of the geochemical element sampling area, the sampling points corresponding to the second sample collection result do not have to be the same as the sampling points corresponding to the first sample sampling result.

[0095] S110. Obtain the reserve change amount data of the geochemical elements in the sampling area by combining the first reserve data and the second reserve data.

[0096] Further, obtaining the reserve change amount data of the geochemical elements in the sampling area by combining the first reserve data and the second reserve data includes:

[0097] Taking the difference between the first reserve value and the second reserve value to obtain the reserve change amount of the geochemical elements in the sampling area; calculating the standard error of the first standard error and the second standard error to obtain the standard error of the reserve change amount;

[0098] Determining the confidence interval of the reserve change amount of the geochemical elements at a preset significance level according to the reserve change amount of the geochemical elements and the standard error of the reserve change amount.

[0099] Among them, since it is necessary to obtain the reserve change amount of the geochemical elements in the sampling area, the two sample collection results correspond to different periods of the same area, that is, the early stage corresponding to the first sample collection result and the later stage corresponding to the second sample collection result. The calculation formula for the reserve change amount of the geochemical elements is:

[0100] SCS d = SCS2 - SCS1

[0101] In the formula, SCS2 is the reserve value of the sampling area corresponding to the second sample collection result, and SCS1 is the reserve value of the sampling area corresponding to the first sample collection result.

[0102] The calculation formula for the standard error of the reserve change amount of the geochemical elements is:

[0103]

[0104] In the formula, SE SCS1 is the standard error corresponding to the first reserve data, and SE SCS2 is the standard error corresponding to the second reserve data.

[0105] Further, taking the significance level α = 0.05, the confidence interval of the reserve change amount of the geochemical elements is: [SCS d - 1.96SE SCSd , SCS d + 1.96SE SCSd . If 0 is within the above interval, it cannot be considered that there is a significant change in the reserve of the geochemical elements between the two periods. If both the upper and lower bounds of this interval are positive values, it is considered that the reserve of the geochemical elements in the sampling area has increased. If both the upper and lower bounds of this interval are negative values, it is considered that the reserve of the geochemical elements in the sampling area has decreased.

[0106] Figure 5 Fig. shows a schematic flowchart of another method for obtaining the first prediction model provided by the embodiment of the present application.

[0107] Further, as Figure 5 shown, the first set includes a first training sample set, and the first prediction model is trained and generated in the following manner:

[0108] S402, generate the first training sample set, where the first training sample set includes a preset number of sampling points and the density values of the geochemical element reserves corresponding to the sampling points;

[0109] S404, train a preset model with the samples in the first training sample set, use the sampling points as inputs, and use the virtual density values of the geochemical element reserves corresponding to the sampling points as outputs. When the unification rate between the output virtual density value and the density value of the geochemical element reserves meets a preset threshold, complete the training of the preset model and generate the first prediction model.

[0110] By adopting the above technical solution, using the density values of the geochemical element reserves obtained from the sampling area collection results as the training extrusion of the deep learning model, training to obtain the virtual density values of the geochemical element reserves corresponding to the sample points, and judging the authenticity of the virtual density values according to the density values of the geochemical element reserves, and then feeding the judgment results back into the deep learning model, so that the virtual density values generated by the deep learning model subsequently are closer to the density values of the geochemical element reserves. When the unification rate between the output virtual density value and the density value of the geochemical element reserves meets a preset threshold, thus complete the training of the preset model and generate the first prediction model.

[0111] Figure 6 shows a schematic flow diagram of obtaining the mean value and standard error of the density value of geochemical element reserves in this application. As Figure 6 shown, before obtaining the first reserve data of geochemical elements based on the first prediction model and the sampling area area, the following steps are included:

[0112] S502, judge whether the sampling quantity of the sampling points exceeds a preset quantity threshold;

[0113] S504, if the sampling quantity exceeds the quantity threshold, then obtain the mean value and standard error of the density value of the geochemical element reserves in the sampling area according to the first set;

[0114] S506, if the sampling quantity does not exceed the quantity threshold, then process the first set by interpolation method to obtain a prediction curve model;

[0115] S508, obtain the mean value and standard error of the density value of the geochemical element reserves in the sampling area according to the prediction curve model.

[0116] Among them, when the number of samples at the sampling points is small, it may lead to inaccurate prediction results. Therefore, the number of samples is judged by a preset quantity threshold. When the number of samples does not exceed the quantity threshold, interpolation can be used to obtain a prediction curve with higher prediction accuracy to calculate the mean and standard error of the geochemical element reserve density values.

[0117] Further, after calculating and obtaining the first set based on the first sample collection result of the sampling area, it further includes:

[0118] Judging whether the density values of the geochemical element reserves in the first set are abnormal. The abnormality includes whether the density values of the geochemical element reserves corresponding to the sampling points exceed the preset range. If the density values are abnormal, a resampling reminder for the corresponding sampling points will appear.

[0119] Among them, since the density values of the geochemical element reserves are usually collected by manual sampling, the range of the density values of the preset geochemical element reserves is used for judgment, so as to avoid the interference of obvious abnormal density values on the collection results.

[0120] The above is the introduction of the method embodiment. The following further illustrates the solution of the present application through the system embodiment.

[0121] Figure 7 The structural schematic diagram of a calculation system for the change amount of geochemical element reserves provided by an embodiment of the present application is shown. The calculation system includes:

[0122] The first collection module 61 is used to calculate and obtain a first set based on the first sample collection result of the sampling area. The first set includes multiple density values of geochemical element reserves corresponding to the sampling points;

[0123] The first model processing module 62 is used to obtain a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain the first reserve data of the geochemical element based on the first prediction model and the sampling area;

[0124] The second collection module 63 is used to calculate and obtain a second set based on the second sample collection result of the sampling area. The second set includes multiple density values of geochemical element reserves corresponding to the sampling points;

[0125] The second model processing module 64 is used to obtain a second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain the second reserve data of the geochemical element based on the second prediction model and the sampling area;

[0126] A calculation module 65, configured to obtain the data of the change in the reserves of geochemical elements in the sampling area by combining the first reserve data and the second reserve data.

[0127] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the above-described modules can refer to the corresponding processes in the foregoing method embodiments and will not be elaborated herein.

[0128] The implementation principle of this embodiment is as follows:

[0129] A first acquisition module 61, configured to calculate and obtain a first set based on the first sample acquisition result of the sampling area, where the first set includes density values of the reserves of geochemical elements corresponding to a plurality of sampling points; a first model processing module 62, configured to obtain a first prediction model of the density value of the geochemical elements in the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain the first reserve data of the geochemical elements based on the first prediction model and the area of the sampling area; a second acquisition module 63, configured to calculate and obtain a second set based on the second sample acquisition result of the sampling area, where the second set includes density values of the reserves of geochemical elements corresponding to a plurality of sampling points; a second model processing module 64, configured to obtain a second prediction model of the density value of the geochemical elements in the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain the second reserve data of the geochemical elements based on the second prediction model and the area of the sampling area; a calculation module 65, configured to obtain the data of the change in the reserves of geochemical elements in the sampling area by combining the first reserve data and the second reserve data. By means of the present application, the spatial heterogeneity factors of the sampling points can be considered in the calculation process of the reserves of geochemical elements, greatly improving the accuracy of the calculation results and being able to give the error range of the results.

[0130] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0131] The above are all the preferred embodiments of the present application, and the protection scope of the present application is not limited thereby. Therefore, all equivalent changes made according to the structure, shape, and principle of the present application should be covered within the protection scope of the present application.

Claims

1. A calculation method for the change amount of the reserves of geochemical elements, characterized in that Including: Calculating and obtaining a first set based on the first sample collection result of the sampling area, the first set including density values of geochemical element reserves corresponding to multiple sampling points; Combining the first set and the spatial position relationship of the sampling points to obtain a first prediction model of the geochemical element density value of the sampling area, and obtaining first reserve data of the geochemical element based on the first prediction model and the sampling area; Calculating and obtaining a second set based on the second sample collection result of the sampling area, the second set including density values of geochemical element reserves corresponding to multiple sampling points; Combining the second set and the spatial position relationship of the sampling points to obtain a second prediction model of the geochemical element density value of the sampling area, and obtaining second reserve data of the geochemical element based on the second prediction model and the sampling area; Combining the first reserve data and the second reserve data to obtain reserve change amount data of the geochemical element in the sampling area; both the first prediction model and the second prediction model are prediction models for the reserve value of the geochemical element sampling area, and the sampling points corresponding to the second sample collection result do not have to be the same as the sampling points corresponding to the first sample collection result; The two sample collection results correspond to different periods of the same area, i.e., the early stage corresponding to the first sample collection result and the later stage corresponding to the second sample collection result; The obtaining of the first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points includes: Calculating an empirical semi-variogram of the geochemical element reserve density value based on the first set and the spatial position relationship of the sampling points; Using an approximation criterion to screen a target model to fit with the empirical semi-variogram to obtain a theoretical semi-variogram; Based on the theoretical semi-variogram, obtaining the first prediction model of the geochemical element density value of the sampling area by using block kriging; the first reserve data includes the first reserve value of the geochemical element in the sampling area and the corresponding first standard error, and the obtaining of the first reserve data of the geochemical element based on the first prediction model and the sampling area includes: Obtaining a first block kriging predicted value and a first variance of the geochemical element density value of the sampling area based on the first prediction model; obtaining the first reserve value according to the block kriging predicted value of the geochemical element density value and the area of the sampling area, and obtaining the first standard error by using the first variance of the geochemical element density value and the area of the sampling area.

2. The calculation method according to claim 1, wherein The obtaining of the second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points includes: Calculating an empirical semi-variogram of the geochemical element reserve density value based on the second set and the spatial position relationship of the sampling points; Using an approximation criterion to screen a target model to fit with the empirical semi-variogram to obtain a theoretical semi-variogram; Based on the theoretical semi-variogram, obtaining the second prediction model of the geochemical element density value of the sampling area by using block kriging.

3. The calculation method according to claim 2, wherein The second reserve data includes the second reserve value of the geochemical elements in the sampling area and the corresponding second standard error. The obtaining of the second reserve data of the geochemical elements based on the second prediction model and the sampling area includes: Obtaining the second Kriging prediction value and the second variance of the geochemical element density value in the sampling area based on the second prediction model; obtaining the second reserve value according to the Kriging prediction value of the geochemical element density value and the area of the sampling area, and obtaining the second standard error by using the second variance of the geochemical element density value and the area of the sampling area.

4. The calculation method according to claim 3, characterized in that, The obtaining of the reserve change data of the geochemical elements in the sampling area by combining the first reserve data and the second reserve data includes: Obtaining the reserve change of the geochemical elements in the sampling area by taking the difference between the first reserve value and the second reserve value; Calculating the first standard error and the second standard error to obtain the standard error of the reserve change; Determining the confidence interval of the reserve change of the geochemical elements at a preset significance level according to the reserve change of the geochemical elements and the standard error of the reserve change.

5. The calculation method according to claim 1, characterized in that The first set includes a first training sample set, and the first prediction model is generated through the following method: Generating the first training sample set, where the first training sample set includes a preset number of sampling points and the density values of the geochemical element reserves corresponding to the sampling points; Training a preset model with the samples in the first training sample set, using the sampling points as inputs and the virtual density values of the geochemical element reserves corresponding to the sampling points as outputs. When the unification rate of the output virtual density values and the density values of the geochemical element reserves meets a preset threshold, the training of the preset model is completed and the first prediction model is generated.

6. The calculation method according to claim 1, characterized in that, The obtaining of the first reserve data of the geochemical elements based on the first prediction model and the sampling area includes: Judging whether the sampling quantity of the sampling points exceeds a preset quantity threshold; If the sampling quantity exceeds the quantity threshold, obtaining the mean value and the standard error of the geochemical element reserve density value in the sampling area according to the first set; If the sampling quantity does not exceed the quantity threshold, processing the first set by interpolation to obtain a prediction curve model; Obtaining the mean value and the standard error of the geochemical element reserve density value in the sampling area according to the prediction curve model.

7. The calculation method according to claim 1, characterized in that After calculating and obtaining the first set based on the first sample collection result of the sampling area, it further includes: Judging whether there is an abnormality in the density values of the geochemical element reserves in the first set. The abnormality includes whether the density values of the geochemical element reserves corresponding to the sampling points exceed a preset range. If the density value is abnormal, a reminder for re-sampling the corresponding sampling point will appear.

8. A calculation system for the change amount of the reserves of geochemical elements, characterized in that, Including: A first collection module, configured to calculate and obtain a first set based on the first sample collection result of the sampling area. The first set includes multiple density values of the geochemical element reserves corresponding to the sampling points; The first model processing module is used to obtain a first prediction model of the geochemical element density value of the sampling area by combining the first set and the spatial position relationship of the sampling points, and obtain the first reserve data of the geochemical elements based on the first prediction model and the sampling area; The second acquisition module is used to calculate and obtain a second set based on the second sample acquisition result of the sampling area, and the second set includes density values of geochemical element reserves corresponding to multiple sampling points; The second model processing module is used to obtain a second prediction model of the geochemical element density value of the sampling area by combining the second set and the spatial position relationship of the sampling points, and obtain the second reserve data of the geochemical elements based on the second prediction model and the sampling area; The calculation module is used to obtain the reserve change amount data of the geochemical elements in the sampling area by combining the first reserve data and the second reserve data; Both the first prediction model and the second prediction model are prediction models for the reserve value of the geochemical element sampling area, and the sampling points corresponding to the second sample acquisition result do not have to be the same as the sampling points corresponding to the first sample acquisition result; The two sample acquisition results correspond to different periods of the same area, namely the early stage corresponding to the first sample acquisition result and the later stage corresponding to the second sample acquisition result.

Images