General soil erodibility K value calculation method based on feature structure operator
Through soil sample point investigation and characteristic structure operator calculation, the regional limitation problem of soil erosion K value estimation in the existing technology is solved, and the accurate calculation and regional applicability of soil erosion K value are achieved, providing key parameters for soil erosion monitoring.
Patent Information
- Application Number
- CN202510288287.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-05-26
- Publication Date
- 2025-07-11
AI Technical Summary
In the prior art, foreign soil erosion K value estimation methods cannot be directly applied to my country, and the existing methods have large errors in the soil types in my country. It is impossible to establish a general soil erosion K value estimation method, and cannot accurately reflect the soil erosion characteristics, and there are regional limitations.
Through soil sample point survey, soil physical and chemical attribute characteristic data were obtained, characteristic structure operators such as soil pH, structure, texture and hydraulic power were calculated, and soil erosion K value was calculated based on operator coefficients, and spatial prediction of the target area was carried out.
A general soil erosion K value calculation method based on characteristic structure operator is provided, which can accurately reflect the characteristics of soil erosion and provide key parameters for regional soil erosion monitoring and soil and water conservation. The prediction results have good applicability and accuracy in different types of areas.
Smart Images
Figure CN120294291A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of soil erosion and water and soil conservation, and particularly to a general soil erodibility K value calculation method based on a characteristic structure operator. Background Art
[0002] The soil erodibility K value is a comprehensive index to measure the ease of soil erosion, reflecting the magnitude of the internal erosion force of the soil, and is closely related to the physical and chemical properties of the soil itself, that is, the material composition and structural state of the soil. The relational expressions for obtaining the soil erodibility K value established by domestic and foreign scholars all have obvious regional characteristics and are only suitable for estimating the soil erodibility K value in their respective research areas. The experimental observation results of the soil erodibility K value of subtropical soils in China in the prior art show that most soil types in this region cannot directly use the nomograph method to estimate the soil erodibility factor K value. The actual observation results also show that the existing foreign soil erodibility estimation methods cannot be directly applied to the estimation of the soil erodibility K value in China. In the prior art, the K values of typical cultivated soils in Guizhou calculated by various commonly used international methods have a large gap with the measured values; it is found that the K value estimation methods of the EPIC and USLE models both "overestimate" the soil erodibility in the thin-layer black soil area, and a correction method is proposed; a correction formula is also proposed for the "nomograph" results in the black soil area: Ks = 0.9823Kg + 0.0269 (Ks is the measured value, Kg is the estimated value). Therefore, it is particularly necessary to establish a general soil erodibility K value estimation method in China.
[0003] In the prior art, it is considered that the soil erodibility K value is variable and varies greatly with the rainfall intensity, especially in the small rainfall intensity range, and the variation range is the most significant. Research shows that the soil erodibility K value increases significantly with the increase of the previous soil water content; there are large differences in the soil erodibility K value between years, and the difference between the highest and lowest values of the same soil erodibility K value between years is nearly 8 times, and it is high in spring and summer and low in autumn and winter, and the difference between the seasonal highest and lowest values is as high as 4 times; the seasonal average erodibility K value of the soil over the years is closely related to the seasonal average rainfall, temperature, soil water content and bulk density.
[0004] Therefore, how to construct a structural operator for the soil erodibility K value through these influencing factors, establish a general soil erodibility K value estimation method, accurately and sensitively reflect the soil erodibility characteristics, and overcome the regional limitations of the existing algorithms is the key technical problem to be solved by the method of the present invention. Summary of the Invention
[0005] The object of the present invention is to provide a general soil erodibility K value calculation method based on a characteristic structure operator, which can provide key parameters of soil erodibility for regional soil and water loss monitoring and soil and water conservation benefit evaluation.
[0006] To achieve the above object, the present invention provides the following solutions:
[0007] A general soil erodibility K value calculation method based on a feature structure operator, comprising the following steps:
[0008] S1. Obtain the surface soil physical and chemical property characteristic data of the surveyed sample points through soil sample point surveys, soil sample collection, experimental analysis, and testing;
[0009] S2. Based on the surface soil physical and chemical property characteristic data of the surveyed sample points, calculate the feature structure operators of soil erodibility respectively, where the feature structure operators include the soil pH feature structure operator, the soil structure feature structure operator, the soil texture feature structure operator, and the soil hydraulic feature structure operator;
[0010] S3. Based on the feature structure operators calculated in S2, calculate the soil erodibility K value of the surveyed sample points;
[0011] S4. Conduct spatial prediction of the target area based on the calculated soil erodibility K value.
[0012] Preferably, in S1, the obtained surface soil physical and chemical property characteristic data of the surveyed sample points include soil pH, soil organic matter content, soil hydraulic characteristic parameters, soil bulk density, and soil mechanical composition.
[0013] Preferably, in S2, the formula for calculating the soil pH feature structure operator is as follows:
[0014] I = (1.03 - log10(pH)) -0.9946
[0015] Wherein, I represents the soil pH feature structure operator, and pH represents the soil pH.
[0016] Preferably, in S2, the formula for calculating the soil structure feature structure operator is as follows:
[0017] S = (1 - bulk / 10) 4.3488 ×(1 - omc / 100) -5.02474
[0018] Wherein, S is the soil structure feature structure operator, bulk is the soil bulk density, and omc is the soil organic matter content.
[0019] Preferably, in S2, the formula for calculating the soil texture feature structure operator is as follows:
[0020] T = (1 - clay / 100) 5.8882 ×(1 - sand / 100) 3.8301 ×(1 - silt / 100) 1.9064 ×(mm / 10)-0.7835
[0021] Among them, T is the soil texture characteristic structure operator, clay is the clay content, sand is the sand content, silt is the silt content, and mm is the ratio of the silt content to the clay content.
[0022] Preferably, in S2, the formula for calculating the soil hydraulic characteristic structure operator is as follows:
[0023] O = (1 - fc / 100) 7.4736 ×(1 - wc / 100) -6.9427
[0024] Among them, O is the soil hydraulic characteristic structure operator, fc is the field capacity, and wc is the wilting coefficient. Preferably, in S3, the calculation formula for the soil erodibility K value of the survey sample point is as follows:
[0025] K = §×I×S×T×O
[0026] Among them, § is the operator structure coefficient.
[0027] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements a general soil erodibility K value calculation method based on a characteristic structure operator as described in any one of the above.
[0028] According to the specific embodiments provided by the present invention, the following technical effects are disclosed by the present invention:
[0029] The present invention provides a general soil erodibility K value calculation method based on a characteristic structure operator. This method calculates and determines the characteristic structure operators of soil acidity, soil structure, soil texture, and soil hydraulics respectively by obtaining and establishing a physical and chemical property characteristic database of soil sample points in the target area, using the algorithms of 4 newly invented soil erodibility characteristic structure operators, and then comprehensively calculates the soil erodibility K value based on this, which can provide key parameters of soil erodibility for regional soil erosion monitoring and soil and water conservation benefit evaluation. BRIEF DESCRIPTION OF THE DRAWINGS
[0030] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required to be used in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.
[0031] Figure 1 It is a flowchart of a general soil erodibility K value calculation method based on a characteristic structure operator provided by the present invention;
[0032] Figure 2 Relationship diagram between the predicted residual probability distribution of soil erodibility K value and the normal mathematical expectation in the embodiment of the present invention (n = 19);
[0033] Figure 3 Relationship diagram between the predicted value and the observed value of soil erodibility K in the embodiment of the present invention (n = 19);
[0034] Figure 4 Relationship diagram between the predicted value and the observed value of soil erodibility K in the embodiment of the present invention (n = 5);
[0035] Figure 5 Distribution diagram of the predicted value and the prediction residual of soil erodibility K in the test area in the embodiment of the present invention. Specific implementation manner
[0036] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
[0037] To make the above objects, features, and advantages of the present invention more obvious and understandable, the present invention will be further described in detail below in conjunction with the accompanying drawings and specific implementation manners.
[0038] As Figure 1 shown, a general soil erodibility K value calculation method based on a feature structure operator provided by the present invention includes the following steps:
[0039] S1. Obtain the surface soil physical and chemical property characteristic data of the survey sample points through soil sample point investigation, soil sample collection, experimental analysis, and testing;
[0040] S2. Based on the surface soil physical and chemical property characteristic data of the survey sample points, calculate the feature structure operators of soil erodibility respectively, where the feature structure operators include the soil pH feature structure operator, the soil structure feature structure operator, the soil texture feature structure operator, and the soil hydraulic feature structure operator;
[0041] S3. Based on the feature structure operators calculated in S2, calculate the soil erodibility K value of the survey sample points;
[0042] S4. Conduct spatial prediction of the target area based on the calculated soil erodibility K value.
[0043] Further, in S1, the obtained surface soil physical and chemical property characteristic data of the survey sample points include soil pH, soil organic matter content, soil hydraulic characteristic parameters, soil bulk density, and soil mechanical composition.
[0044] Furthermore, in S2, the formula for calculating the soil pH characteristic structure operator is as follows:
[0045] I = (1.03 - log10(pH)) -0.9946
[0046] where I represents the soil pH characteristic structure operator and pH represents the soil pH value.
[0047] Further, in the said S2, the formula for calculating the soil structure characteristic structure operator is as follows:
[0048] S = (1 - bulk / 10) 4.3488 ×(1 - omc / 100) -5.02474
[0049] where S is the soil structure characteristic structure operator, bulk is the soil bulk density (g / cm3), and omc is the soil organic matter content (%).
[0050] Furthermore, in S2, the formula for calculating the soil texture characteristic structure operator is as follows:
[0051] T = (1 - clay / 100) 5.8882 ×(1 - sand / 100) 3.8301 ×(1 - silt / 100) 1.9064 ×(mm / 10) -0.7835
[0052] where T is the soil texture characteristic structure operator, clay is the clay (<0.002mm) content (%), sand is the sand (2 - 0.05mm) content (%), silt is the silt (0.05 - 0.002mm) content (%), and mm is the ratio of silt content to clay content (silt / clay) (dimensionless).
[0053] Still further, in S2, the formula for calculating the soil hydraulic characteristic structure operator is as follows:
[0054] O = (1 - fc / 100) 7.4736 ×(1 - wc / 100) -6.9427
[0055] where O is the soil hydraulic characteristic structure operator, fc is the field capacity (V / V, %), and wc is the wilting coefficient (V / V, %).
[0056] Further, in S3, the calculation formula for the soil erodibility K value of the surveyed sample points is as follows:
[0057] K = § × I × S × T × O
[0058] Among them, § is the operator structure coefficient, and usually can take the value of 11.9031.
[0059] In a specific embodiment, based on the data of 7,279 typical soil profile sample points in China during the second national soil census in the 1980s, the soil erodibility K value was calculated by using the method of the present invention. The specific implementation methods and steps are as follows:
[0060] (1) Establish a sample point soil property database. Based on the surface soil physical and chemical property data of 7,279 typical soil profile sample points across the country, such as soil pH, soil bulk density, soil organic matter content, soil mechanical composition (sand, silt, clay content), etc., a sample point soil property database is established.
[0061] (2) Calculate the numerical values of the soil erodibility characteristic structure operators. Based on the characteristic data of the sample point soil physical and chemical properties, the numerical values of the soil erodibility characteristic structure operators are calculated respectively by using the algorithms of the characteristic structure operators for soil acidity and alkalinity, soil structure, soil texture, soil hydraulics, etc.
[0062] (3) Calculate the soil erodibility K value. By using the soil erodibility characteristic structure operator and through a comprehensive calculation method (Equation 10), the surface soil erodibility K value (0.1317 t·hm²·h / hm²·MJ·mm) of each soil profile sample point is calculated.
[0063] (4) Based on the soil types of the profile sample points and their surface soil erodibility K values, the statistical distribution characteristics of the soil erodibility K values of each type (soil type) are obtained by using mathematical statistics, as shown in Table 1 below; by using the soil type GIS connection method, the soil profile sample points are connected with the soil type vector graphic patches, and then superimposed with the basin vector distribution map. Based on the soil type distribution map patches, the average soil erodibility K value of each basin is statistically calculated, as shown in Table 2 below.
[0064] Table 1
[0065]
[0066]
[0067]
[0068]
[0069]
[0070] Table 2
[0071]
[0072]
[0073]
[0074] Establish long-term experimental observation datasets for 19 soil erosion regions, such as the red soil region in southern China, the Loess Plateau region, the black soil region in Northeast China, and the purple soil region in Southwest China. The correlation matrix analysis between the soil erodibility K value and soil physical and chemical properties shows that there is a significant correlation between the soil erodibility K value and the soil organic matter content, while the correlation with other soil physical and chemical property indicators is weak. However, for the two principal components (F1, F2) formed by the principal component analysis of the characteristic indicators of soil physical and chemical properties in the erosion regions, the interpretation rates for the soil erodibility K value are 44.3% and 24.7% respectively, and the cumulative interpretation rate can reach 69.0%. And except for the soil gravel content and organic matter content, the contents of soil clay, silt, sand, as well as the soil silt-clay ratio, bulk density, field water holding capacity, wilting coefficient, and soil pH value and other 8 soil physical and chemical property characteristic indicators are significantly correlated with each other or with the F1 and F2 principal components. Therefore, the contents of soil clay, silt, and sand, soil bulk density, silt-clay ratio, soil pH, field water holding capacity, wilting coefficient, and organic matter content are used as characteristic variables affecting the soil erodibility K value. According to the soil physical and chemical property indicators and their influencing characteristic types, the 9 variable factors are divided into 4 types of structural factors, namely I soil acidity factor (soil pH), S soil structure factor (soil bulk density, soil organic matter), T soil texture factor (contents of soil clay, silt, sand and soil silt-clay ratio), and O soil hydraulic factor (field water holding capacity, wilting coefficient).
[0075] For the soil erodibility K value prediction method constructed based on the long-term experimental observation datasets of 19 soil erosion regions, through the analysis of the proportion of the absolute residual of the predicted value in the annual observation range of the 13 erosion test regions independently carried out, it is found that only 1 region is 26%, and the other regions are all less than 15%; 7 regions are less than 10%, 4 regions are less than 5%, and 1 region is only 0.4%. The smaller the proportion of the absolute residual of the predicted value in the annual observation range, the greater the probability that the predicted value is within the range of annual observation values; for all regions participating in the model parameter estimation, this prediction residual is much smaller than the annual observation range, indicating the accuracy of the prediction results. The probability of the prediction residual of the soil erodibility K shows a significant normal distribution law as Figure 2 shown, indicating that the probability of the prediction residual of the soil erodibility K being extremely large or extremely small is small, that is, the probability of the prediction result of the soil erodibility K being extremely accurate or extremely inaccurate is small, and most predictions are within a reasonable range, reflecting the marginal condition effect of the model prediction and the stability of the model prediction.
[0076] Meanwhile, the predicted values of soil erodibility K for the 19 erosion test areas used for modeling showed a significant linear relationship with the observed values, and the prediction interpretation rate or certainty of soil erodibility K could reach 85% (R2 = 0.8550, Adjusted R2 = 0.8465), as Figure 3 shown.
[0077] Six soil erosion test areas, such as Huairou cinnamon soil SS2 in the northern rocky soil and mountainous area, Miyun skeletal cinnamon soil SS1, Anxi red soil RS1 in the southern red soil area, Yujiang red soil P18, Yuexi brownish red soil RS2, and Suining purple soil PS1 in the southwestern purple soil area, can be used for model testing. Except for the SS1 area where the soil physical and chemical property characteristics are not detailed enough, the predicted values of soil erodibility K and the observed values in the remaining areas have a significant correlation, as Figure 4 shown. The average proportion of the absolute residuals of the predicted values of soil erodibility K in the 19 modeling areas to the annual observed range is 5.6%. Calculated according to this proportion, the predicted absolute residuals in the remaining 5 test areas except SS1 are less than the annual observed range, and the predicted values of soil erodibility K are also within the range of the annual observed values, as Figure 5 shown. Generally speaking, the constructed method for soil erodibility K value has good applicability in different types of areas across the country, and the predicted results of soil erodibility K value can basically reflect the soil erodibility characteristics of each type of area.
[0078] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it implements a general soil erodibility K value calculation method based on a feature structure operator as described in any one of the above.
[0079] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the above technical solution, in essence, or the part that contributes to the prior art, can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., including several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.
[0080] In this article, specific examples are used to elaborate on the principle and implementation manner of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific implementation manner and application scope. In summary, the content of this specification should not be construed as a limitation to the present invention.
Claims
1. A general method for calculating the soil erodibility K value based on a feature structure operator, characterized in that It includes the following steps: S1. Obtain the characteristic data of the physical and chemical properties of the surface soil at the investigation sampling points through soil sampling point surveys, soil sample collection, experimental analysis, and testing; S2. Based on the characteristic data of the physical and chemical properties of the surface soil at the investigation sampling points, calculate the characteristic structure operators of soil erodibility respectively. Among them, the characteristic structure operators include the characteristic structure operator of soil pH, the characteristic structure operator of soil structure, the characteristic structure operator of soil texture, and the characteristic structure operator of soil hydraulic characteristics; S3. Based on the characteristic structure operators calculated in S2, calculate the soil erodibility K value of the investigation sampling points; S4. Conduct spatial prediction of the target area based on the calculated soil erodibility K value.
2. The general soil erodibility K value calculation method based on the feature structure operator according to claim 1, wherein, In the above S1, the characteristic data of the physical and chemical properties of the surface soil obtained includes soil pH, soil organic matter content, soil hydraulic characteristic parameters, soil bulk density, and soil mechanical composition.
3. A general soil erodibility K value calculation method based on a feature structure operator according to claim 1, characterized in that In the above S2, the formula for calculating the characteristic structure operator of soil pH is as follows: I = (1.03 - log10(pH)) -0.9946 Where I represents the characteristic structure operator of soil pH, and pH represents soil pH.
4. A general soil erodibility K value calculation method based on a feature structure operator according to claim 3, characterized in that In the above S2, the formula for calculating the characteristic structure operator of soil structure is as follows: S=(1 - bulk / 10) 4.3488 ×(1 - omc / 100) -5.02474 Where S is the characteristic structure operator of soil structure, bulk is soil bulk density, and omc is soil organic matter content.
5. A general soil erodibility K value calculation method based on a feature structure operator according to claim 4, characterized in that, In the above S2, the formula for calculating the characteristic structure operator of soil texture is as follows: T = (1 - clay / 100) 5.8882 ×(1 - sand / 100) 3.8301 ×(1 - silt / 100) 1.9064 ×(mm / 10) -0.7835 Where T is the characteristic structure operator of soil texture, clay is clay content, sand is sand content, silt is silt content, and mm is the ratio of silt content to clay content.
6. The general soil erodibility K value calculation method based on a feature structure operator according to claim 5, characterized in that, In the above S2, the formula for calculating the characteristic structure operator of soil hydraulic characteristics is as follows: O = (1 - fc / 100) 7.4736 ×(1 - wc / 100) -6.9427 Where O is the characteristic structure operator of soil hydraulic characteristics, fc is field capacity, and wc is wilting coefficient.
7. A general soil erodibility K value calculation method based on a feature structure operator according to claim 1, characterized in that, In the above S3, the formula for calculating the soil erodibility K value of the investigation sampling points is as follows: K = § × I × S × T × O Where § is the operator structure coefficient.
8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements a general soil erodibility K value calculation method based on characteristic structure operators as described in any one of claims 1 to 7.