Full-slope erosion sediment transport simulation method integrating soil properties and land utilization

Abstract

The invention discloses a full-slope erosion sediment transport simulation method integrating soil properties and land utilization, which comprises the following steps: respectively collecting erosion sediment transport basic data corresponding to full slopes of a ridge and loess hills slope and a gully slope in a drainage basin, including soil data, topographic data, land utilization data, vegetation data and water and sediment data; automatically judging the erosion environment type of the whole slope based on erosion sediment transport basic data, topographic conditions and land utilization type boundary conditions, and obtaining initial to-be-optimized parameters influencing erosion based on the erosion environment type; updating the initial to-be-optimized parameters based on a cuckoo search algorithm until optimal parameters are obtained; the slope combination and the slope length in the topographic condition and the land utilization type and vegetation coverage rate combination in the land utilization type boundary condition are changed, multiple different scenes are obtained, the whole slope sediment production and transportation evolution process is simulated based on the different scenes and the optimal parameters, and the slope sediment transportation rates under different conditions are obtained.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of slope soil erosion prediction, and a full-slope erosion sediment transport simulation method integrating soil properties and land use. BACKGROUND

[0002] Slopes are important sources of river sediment, and severe slope erosion leads to land degradation and a series of environmental, social and economic problems. With the implementation of major ecological restoration projects in recent years, the relationship between water and sediment has changed significantly, posing greater challenges for water and sediment process simulation and prediction, especially at the sub-rainfall scale. It is an important requirement for the simulation and prediction of watershed sediment production and transport processes and the benefit assessment of soil and water conservation measures to establish a soil erosion model with a widely adaptable physical mechanism for full-slope erosion sediment production, considering factors such as soil properties, vegetation cover, topography and land use change.

[0003] In related technologies, numerous soil erosion models have been developed, which have undergone development from empirical models to physical models. Empirical models are simple in structure, easy to use, and accurate in predicting multi-year average soil loss, but they are difficult to apply at the sub-rainfall scale and below, do not consider the impact of key factors on soil separation and sediment transport processes, have poor portability, and are difficult to control in terms of precision when extending the boundary conditions or region, and the model structure is unstable. Soil erosion models based on physical processes are developed to try to break the shackles of application limitations, but due to the existence of a large number of complex or assumed parameters and processes, and the development under relatively ideal conditions, they lack consideration of the direct impact of changing underlying surfaces such as topography, vegetation and land use on erosion and sediment transport, making them still unable to have more extensive extension and still have a certain distance from practical application.

[0004] Therefore, how to realize the effective expression of the key influencing factors of soil erosion physical models on the slope erosion of the Loess Plateau, and continuously simulate the gully-slope differentiation erosion and sediment transport processes generated by the ridge and valley slopes has become a problem to be solved. SUMMARY

[0005] Therefore, the embodiment of the present application provides a full-slope erosion sediment transport simulation method integrating soil properties and land use, which solves the problem of insufficient consideration of key influencing factors of soil erosion physical models on slope erosion of the Loess Plateau in related technologies, and the inability to continuously simulate the gully-slope differentiation erosion and sediment transport processes generated by the ridge and valley slopes.

[0006] According to a first aspect of the embodiment of the present application, a full-slope erosion sediment transport simulation method integrating soil properties and land use is provided, comprising: Collecting corresponding erosion sediment transport basic data of full slope surface of the gully slope and the valley slope in the watershed respectively, the erosion sediment transport basic data including soil data, terrain data, land use data, vegetation data and water and sediment data; Automatically judging an erosion environment type of the erosion full slope surface based on the erosion sediment transport basic data, terrain conditions and land use type boundary conditions, and obtaining initial to-be-optimized parameters affecting erosion based on the erosion environment type; Updating the initial to-be-optimized parameters based on the cuckoo search algorithm until optimal parameters are obtained; Changing slope combinations and slope lengths in the terrain conditions, and land use type and vegetation coverage combinations in the land use type boundary conditions to obtain a plurality of different scenarios, and simulating the sediment transport evolution process of the full slope surface based on the plurality of different scenarios and the optimal parameters to obtain the sediment transport rate of the slope surface under the plurality of different scenarios.

[0007] According to a second aspect of the embodiment of the present application, an electronic device is provided, comprising a processor, a memory, a communication interface and a communication bus, the processor, the memory and the communication interface complete communication with each other through the communication bus; the memory is used for storing at least one executable instruction, and the executable instruction makes the processor execute the operation corresponding to the method according to the first aspect.

[0008] According to a third aspect of the embodiment of the present application, a computer storage medium is provided, and the computer storage medium stores a computer program, and the program is executed by a processor to realize the method according to the first aspect.

[0009] According to the scheme provided by the embodiment of the present application, the erosion sediment transport basic data corresponding to the full slope surface of the ridge and gully slope in the watershed is collected respectively, the erosion sediment transport basic data including soil data, terrain data, land use data, vegetation data and water and sediment data; the erosion environment type of the erosion full slope surface is automatically judged based on the erosion sediment transport basic data, terrain conditions and land use type boundary conditions, and the initial optimization parameters affecting erosion are obtained based on the erosion environment type; the initial optimization parameters are updated based on the cuckoo search algorithm until the optimal parameters are obtained; the slope combination and slope length in the terrain conditions, and the land use type and vegetation coverage combination in the land use type boundary conditions are changed to obtain a plurality of different scenes, and the full slope surface sediment yield evolution process is simulated based on the plurality of different scenes and the optimal parameters to obtain the sediment yield rate of the slope surface under the plurality of different scenes. In the method, the erosion sediment transport basic data corresponding to the full slope surface of the actual ridge and gully slope is collected, and the erosion environment type of the slope surface is judged through the erosion sediment transport basic data, and the parameters to be optimized are further initialized and the range, required terrain conditions and land use type boundary conditions are obtained for subsequent simulation, the sediment source-sink effect generated by the ridge and gully slope differentiation is considered, and the spatial heterogeneity of soil properties, terrain, land use and vegetation coverage along the slope surface is reflected, and the problem that the existing soil erosion physical model cannot continuously simulate the ridge and gully slope differentiation generated by the ridge and gully slope is solved. BRIEF DESCRIPTION OF DRAWINGS

[0010] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor based on these drawings. Figure 1 The flowchart of the full slope surface erosion sediment transport simulation method provided by the embodiment of the present application is shown. Figure 2 The effect diagram of the simulated sediment yield rate of the full slope surface under different scenes provided by the embodiment of the present application is shown. Figure 3 The effect diagram of the comparison between the simulated single-width sediment yield rate and the measured single-width sediment yield rate provided by the embodiment of the present application is shown. Figure 4 The structural diagram of an electronic device provided by the embodiment of the present application is shown. DETAILED DESCRIPTION

[0011] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present application with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. The following embodiments are used to describe the present application, but not to limit the scope of the present application. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0012] In the following description, “some embodiments” are related to a subset of all possible embodiments, but it can be understood that “some embodiments” can be the same subset or different subsets of all possible embodiments, and can be combined with each other without conflict.

[0013] It should be noted that the terms “first”, “second”, “third” involved in the embodiments of the present application are only to distinguish similar objects, and do not represent the specific order of the objects. It can be understood that “first”, “second”, “third” can be interchanged in specific order or sequence as allowed, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein.

[0014] Those skilled in the art can understand that, unless otherwise defined, all terms (including technical terms and scientific terms) used herein have the same meaning as generally understood by those skilled in the art to which the embodiments of the present application belong. It should also be understood that terms such as those defined in general dictionaries should be understood as having meanings consistent with those in the prior art, and should not be interpreted in an idealized or overly formal sense unless specifically defined as such herein.

[0015] Figure 1 A flowchart of a comprehensive soil property and land use full-slope erosion and sediment transport simulation method provided by the embodiments of the present application, the comprehensive soil property and land use full-slope erosion and sediment transport simulation method provided by the embodiments of the present application can be executed by an electronic device, for example, a computer, a server, etc.

[0016] As shown in Figure 1 The comprehensive soil property and land use full-slope erosion and sediment transport simulation method comprises: S101, collecting erosion and sediment transport basic data corresponding to the full-slope of beam and gully slopes in the watershed respectively, the erosion and sediment transport basic data including soil data, topographic data, land use data, vegetation data and water and sediment data.

[0017] In the embodiments of the present application, the full slope surface in the watershed can be beam and valley slopes, and the corresponding erosion sediment transport basic data of the beam and valley slopes in the watershed are collected, and the erosion sediment transport basic data includes soil data, terrain data, land use data, vegetation data and water and sediment data. The soil data includes soil organic matter (g / kg) and soil median particle size (mm), the terrain data includes slope (radian) and slope length (m), the land use data includes bare land (set to 0), farmland (set to 1), grassland (set to 2) and forest land (set to 3), the vegetation data includes vegetation coverage (%) and the water and sediment data includes single-width runoff (m / s) and single-width sediment transport rate (kg / (s m). The collected erosion sediment transport basic data is the actual data corresponding to the slope surface of the beam and valley slopes. 2

[0018] S102, automatically determine the erosion environment type of the full slope surface based on the erosion sediment transport basic data, terrain conditions and land use type boundary conditions, and obtain the initial optimization parameters affecting erosion based on the erosion environment type.

[0019] In the embodiments of the present application, the erosion sediment transport basic data, terrain conditions and land use type boundary conditions are input into a preset discrimination program to obtain the erosion environment type of the full slope surface, wherein the erosion environment type includes a single slope type, a composite slope type with the same land use in the upper and lower sections, and a composite slope type with different land use in the upper and lower sections. The single slope type, the composite slope type with the same land use in the upper and lower sections, and the composite slope type with different land use in the upper and lower sections are the slope types in the original natural form of the slope surface. Further, the initial optimization parameters affecting the subsequent erosion of the full slope surface are obtained according to the erosion environment type. The initial optimization parameters are soil separability parameters, sediment transportability parameters, soil erodibility decay index, vegetation coverage decay coefficient, soil separation capacity index of beam and valley slopes, and sediment transport capacity index of beam and valley slopes. The terrain conditions can include slope combinations and slope lengths, and the land use type boundary conditions can be farmland-grassland, forest land-grassland, bare land-bare land and grassland-grassland. The land use type boundary conditions corresponding to the composite slope type with the same land use in the upper and lower sections are bare land-bare land and grassland-grassland, and the land use type boundary conditions corresponding to the composite slope type with different land use in the upper and lower sections are farmland-grassland and forest land-grassland. The terrain conditions and land use type boundary conditions are measured data.

[0020] S103, change the slope combinations and slope lengths in the terrain conditions, and the land use types and vegetation coverage combinations in the land use type boundary conditions to obtain a plurality of different scenarios, and simulate the sediment transport evolution process of the full slope surface based on the plurality of different scenarios and the optimal parameters to obtain the sediment transport rate of the full slope surface under the plurality of different scenarios.

[0021] ​Exemplarily, the erosion environment type obtained through the step S102 is single slope type, and the original condition is that the slope combination is 15%, the slope length is 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%. Different scenarios are obtained after the change, one of which is that the slope is increased from 15% to 35%, the slope length remains unchanged at 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%. The sediment yield evolution process of the scenario and the optimal parameters is simulated to obtain the sediment yield rate of the slope surface in the scenario, and the other scenarios are similar. Through comparison of different scenarios, it can be seen that the changes of the slope, the slope length, the land use type and the vegetation coverage have a significant influence on the soil erosion and the sediment yield rate.

[0022] Exemplarily, the erosion environment type obtained through the step S102 is compound slope type with the same land use on the upper and lower sections, and the slope gradients of the upper and lower sections are inconsistent, the slope gradient of the lower section is steep, and the slope gradient of the upper section is gentle. It is assumed that a mountain slope with a total length of 400 meters is divided into two sections, each with a length of 200 meters, and the land use type is the same, or grassland, or bare land, or farmland. The original condition is that the slope gradient of the first 200 meters is 15%, the slope length is 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%, the slope gradient of the last 200 meters is 15%, the slope length is 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%, that is, the entire slope surface is a uniform slope combination, and the slope gradients of the upper and lower sections are the same, both being 15%. Different scenarios are obtained after the change, one of which is that the slope gradient of the first 200 meters is 15%, the slope length remains unchanged at 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%, the slope gradient of the last 200 meters is increased from 15% to 35%, the slope length remains unchanged at 200 meters, the land use type and the vegetation coverage are grassland with a coverage of 30%. The sediment yield evolution process of the scenario and the optimal parameters is simulated to obtain the sediment yield rate of the slope surface in the scenario, and the other scenarios are similar. Through comparison of different scenarios, it can be seen that the changes of the slope, the slope length, the land use type and the vegetation coverage have a significant influence on the soil erosion and the sediment yield rate.

[0023] As Figure 2 shown, Figure 2 the effect diagram of the simulated sediment yield rate of the entire slope surface in different scenarios provided by the embodiment of the present application. In Figure 2 , the erosion environment type obtained is compound slope type, the slope gradient of the upper section is 15°, the slope gradient of the lower section is 30°, and the simulated land use type boundary condition is that the upper section is bare land-bare land-farmland-grassland, and the lower section is grassland-grassland-woodland-grassland. The horizontal coordinate is different slope length, and the vertical coordinate is the sediment yield rate of the entire slope surface.

[0024] It can be understood that in the embodiments of the present application, the erosion sediment transport basic data corresponding to the full slope of the hillside slope and the valley slope in the watershed are collected respectively, the erosion sediment transport basic data including soil data, terrain data, land use data, vegetation data and water and sediment data; the erosion environment type of the full slope is automatically determined based on the erosion sediment transport basic data, terrain conditions and land use type boundary conditions, and the initial optimization parameters affecting erosion are obtained based on the erosion environment type; the optimization parameters are updated based on the cuckoo search algorithm until the optimal parameters are obtained; the slope combination and slope length in the terrain conditions, and the land use type and vegetation coverage combination in the land use type boundary conditions are changed to obtain a plurality of different scenarios, and the sediment transport rate of the full slope under a plurality of different scenarios is obtained based on a plurality of different scenarios and the optimal parameters. In the method, the erosion sediment transport basic data corresponding to the full slope of the hillside slope and the valley slope are collected, and the erosion environment type of the slope is determined based on the erosion sediment transport basic data, and the initial optimization parameters are further determined and the range thereof, the required terrain conditions and the land use type boundary conditions are obtained for subsequent simulation, the sediment source and sink effect of the hillside slope and the valley slope due to the gully slope differentiation is considered, and the spatial heterogeneity of soil properties, terrain, land use and vegetation coverage along the slope is reflected, and the problem that the existing soil erosion physical model cannot continuously simulate the gully slope differentiation of the hillside slope and the valley slope due to the insufficient consideration of the key influencing factors of the loess plateau slope erosion is solved.

[0025] In the embodiments of the present application, the initial optimization parameters affecting erosion based on the erosion environment type in S102 can be realized by S1021 to S1022, which are described as follows.

[0026] S1021, setting the simulation terrain conditions and the simulation land use type boundary conditions based on the erosion environment type.

[0027] S1022, obtaining the parameter type and the preset parameter range based on the simulation land use type boundary conditions, and randomly generating a group of parameters in the preset parameter range as the initial optimization parameters according to the parameter type.

[0028] In the embodiments of the present application, the simulation terrain conditions and the simulation land use type boundary conditions are set based on the erosion environment type, and the parameter type and the preset parameter range are obtained based on the simulation land use type boundary conditions, and a group of parameters in the preset parameter range are randomly generated as the initial optimization parameters of the subsequent erosion of the full slope, so as to ensure the effectiveness of the simulation.

[0029] In the embodiments of the present application, S103 can be realized by S1031 to S1032, which are described as follows.

[0030] S1031. Obtain the current simulated single-width sediment transport rate based on the initial parameters to be optimized and the single-width sediment transport rate simulation formula, and when the comparison result between the simulated single-width sediment transport rate and the measured single-width sediment transport rate meets the preset accuracy judgment index condition, use the initial parameters to be optimized as the optimal parameters.

[0031] In an embodiment of the present invention, before using the single-width sediment transport rate simulation formula, the slope of the Liangmao slope and the gully slope of the current basin is first judged by the erosion environment type to determine whether it is a single slope type or a compound slope type. When it is a single slope type, the current single-width sediment transport rate is obtained by the single-width sediment transport rate simulation formula corresponding to the parameters to be optimized and the single slope type. When it is a compound slope type, the current single-width sediment transport rate is obtained by the single-width sediment transport rate simulation formula corresponding to the parameters to be optimized and the compound slope type. The preset accuracy judgment index conditions include the determination coefficient , Nash effective coefficient and relative error The three indicators comprehensively evaluate the simulation accuracy, and the effectiveness evaluation criteria are: good quality , meeting the calculation accuracy and does not meet the basic calculation accuracy .

[0032] Furthermore, when the current simulated sediment transport rate per unit width is obtained, the parameter to be optimized is determined as the optimal parameter when the comparison result between the simulated sediment transport rate per unit width and the measured sediment transport rate meets the preset accuracy judgment index condition. Alternatively, when the current corresponding number of iterations reaches the preset upper limit of the number of iterations, the parameter to be optimized is determined as the optimal parameter.

[0033] like Figure 3 As shown, Figure 3 A schematic diagram showing the effect of comparing the simulated single-width sediment transport rate and the measured single-width sediment transport rate provided by an embodiment of the present invention. Figure 3 The simulated land use boundary conditions in the four figures are grassland, farmland, forestland, and bare land, respectively. The horizontal axis represents different dates, and the vertical axis represents the sediment transport rate per width. The sediment transport rate per width includes the simulated sediment transport rate per width (the dashed line in the figure) and the measured sediment transport rate per width (the solid line in the figure). Figure 3 It also includes the coefficient of determination R calculated based on the simulated single-width sediment transport rate and the measured single-width sediment transport rate. 2 , Nash effective coefficient NSE and relative error RE.

[0034] In an embodiment of the present invention, the formula of the coefficient of determination is as follows: ; In the above formula, is the coefficient of determination, is the measured single-width sediment transport rate of the i-th rainfall event, is the simulated sediment yield per unit width of the i-th rainfall event, is the average of the measured sediment yield per unit width; is the average of the simulated sediment yield per unit width.

[0035] The Nash effective coefficient is shown in the following formula: ; In the above formula, is the Nash effective coefficient.

[0036] The relative error is shown in the following formula: ; In the above formula, is the relative error.

[0037] In the embodiments of the present application, the full slope is simultaneously limited by the soil separation capacity and the sediment transport capacity during the erosion and sediment yield process, and the core control equation can be shown in the following formula: ; In the above formula, is the sediment yield per unit width, is the soil separation rate, is the soil separation capacity, is the sediment transport capacity, is the distance from the current position of the full slope to the slope top.

[0038] Wherein, considering the influence of soil properties, land use change and vegetation, the soil separation capacity of Liangmao slope and the sediment transport capacity of Liangmao slope are shown in the following formulas respectively: ; ; In the above formula, is the distance from the current position of the full slope to the slope top; is the soil separability parameter; the subscript represents the slope of Liangmao; is the soil organic matter content of Liangmao slope; is the median particle size of Liangmao slope soil; is the organic matter index; is the median particle size index of soil separation capacity. is the index of ; is the sediment transportability parameter; is the median particle size index of sediment transport capacity; is the index of ; is the soil erodibility decay index; is the vegetation coverage of the ridge-saddle slope, wherein the organic matter index is-0.778; the median particle size index of soil separation capacity is 1.318; the median particle size index of sediment transport capacity is-0.322, is the power of the whole slope runoff of the ridge-saddle slope.

[0039] The power of the whole slope runoff of the ridge-saddle slope is expressed as follows: ; In the above formula, is the average runoff per unit width of the ridge-saddle slope, is the water density, is the slope of the ridge-saddle slope, is the attenuation coefficient of vegetation coverage.

[0040] In the embodiment of the present application, the formula for simulating the sediment transport rate per unit width of a single slope type (ridge-saddle slope) is as follows: ; wherein, is the simulated sediment transport rate per unit width of a single slope type, is the length of the ridge-saddle slope, is the ridge-saddle slope.

[0041] In the embodiment of the present application, the effects of the upper catchment, the significant increase in slope and the change in land use type need to be considered for the gully slope, but the mechanism of the runoff erosion energy does not change, and the core control equation of the sediment production of the gully slope does not change. The formula for the sediment transport rate per unit width of the gully slope considering the soil properties and the change in land use type is derived as follows: ; In the above formula, the subscript represents the gully slope, represents the length of the gully slope; and represent the parameter combination term; is the soil separability parameter of the gully slope; is the transportability parameter of the sediment of the gully slope; is the attenuation index of the soil erodibility of the gully slope; is the attenuation coefficient of the vegetation coverage of the gully slope; is the index of the gully slope ; is the index of the gully slope ; is the vegetation coverage of the gully slope; ​the soil organic matter content of the gully slope; the gully slope gradient, the soil median particle diameter of the gully slope; the initial sediment yield rate.

[0042] wherein the gully slope single-width sediment yield rate and the beam and mound slope single-width sediment yield rate are integrated to obtain a calculation formula of the single-width sediment yield rate of the combined beam and mound slope and gully slope composite slope type, as follows: ; S1032、In the case that the comparison result of the simulated single-width sediment yield rate and the measured single-width sediment yield rate does not meet the preset precision discrimination index condition, the to-be-optimized parameter is iteratively updated again until the comparison result of the simulated single-width sediment yield rate and the measured single-width sediment yield rate corresponding to the updated to-be-optimized parameter meets the preset precision discrimination index condition, and the updated to-be-optimized parameter is taken as the optimal parameter.

[0043] In the embodiment of the present application, in the case that the comparison result between the current simulated single-width sediment yield rate and the measured single-width sediment yield rate does not meet the preset precision discrimination index condition, the to-be-optimized parameter is iteratively updated again through the cuckoo search algorithm, and the current simulated single-width sediment yield rate corresponding to the updated to-be-optimized parameter and the single-width sediment yield rate simulation formula is obtained again to determine whether the comparison result of the current simulated single-width sediment yield rate and the measured single-width sediment yield rate meets the preset precision discrimination index, and the process is repeated until the optimal parameter is obtained.

[0044] wherein, if the maximum number of iterations has been reached in the iteration process but the current simulated single-width sediment yield rate and the measured single-width sediment yield rate do not meet the preset precision discrimination index condition, a parameter precision warning is given, the maximum number of iterations is increased, and the above steps are repeated until the optimal parameter is obtained.

[0045] In the embodiment of the present application, the erosion environment type of the slope surface obtained based on the preset discrimination program and the erosion and sedimentation basic data in S102 can be realized through S102a, which is specifically described as follows.

[0046] S102a, the terrain data, land use data, vegetation data, water and sediment data, terrain conditions and land use type boundary conditions are sorted into the comma-separated value file under the same directory of the preset discrimination program, and the erosion environment type is obtained in combination with the soil data input in the preset discrimination program.

[0047] In the embodiment of the present application, the slope and slope length in the terrain data, the vegetation coverage in the vegetation data, the single-width runoff and single-width sediment rate data in the water and sediment data, the terrain conditions and the land use type boundary conditions are arranged in the preset discrimination program, such as the data.csv in the same directory of the python program, in the example format, and it is ensured that each kind of data has a corresponding column name. Further, the soil organic matter content and soil median particle size basic data in the soil data are provided to the Python program in the input mode. For example, the input can be performed through a command line parameter, a configuration file or interactive input, etc. The preset discrimination program reads the data in the dataa.csv file, and combines the soil organic matter content and soil median particle size data input from the outside to further analyze and process to obtain the erosion environment type.

[0048] In the embodiment of the present application, S104 can be implemented by S1041 to S1042, which are described as follows.

[0049] S1041, the optimal parameter and the corresponding parameter value under different scenarios are substituted into the single-width sediment rate simulation formula to obtain the single-width sediment rate under different scenarios.

[0050] S1042, the sediment rate of the full slope surface under different scenarios is obtained based on the single-width sediment rate under different scenarios.

[0051] In the embodiment of the present application, different scenarios include different slope combinations, slope lengths, land use types, and vegetation coverage rate combinations, each corresponding to a parameter value. The parameter values under different scenarios and the optimal parameter are substituted into the single-width sediment rate simulation formula to obtain the single-width sediment rate under different scenarios. The sediment rate of the full slope surface under different scenarios can be obtained by multiplying the single-width sediment rate and the width of the full slope surface.

[0052] Referring to Figure 4 , a structural schematic diagram of an electronic device according to an embodiment of the present application is shown, and the specific implementation of the electronic device is not limited in the specific embodiments of the present application.

[0053] As Figure 4 shown, the electronic device can include a processor 502, a communications interface 504, a memory 506, and a communications bus 508.

[0054] Among them: The processor 502, the communications interface 504, and the memory 506 complete mutual communication through the communications bus 508.

[0055] The communications interface 504 is configured to communicate with other electronic devices or servers.

[0056] The processor 502 is configured to execute the program 510, and in particular, execute the steps in the above method embodiments.

[0057] In particular, the program 510 can include program codes including computer operation instructions.

[0058] The processor 502 can be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application. The one or more processors included in the smart device can be the same type of processor, such as one or more CPUs; or can be different types of processors, such as one or more CPUs and one or more ASICs.

[0059] The memory 506 is configured to store the program 510. The memory 506 can include a high-speed RAM memory, and can also include a non-volatile memory such as at least one disk memory.

[0060] The program 510 can be specifically configured to cause the processor 502 to perform operations corresponding to the methods described in the above method embodiments.

[0061] The specific implementation of each step in the program 510 can refer to the corresponding description in the corresponding steps and units in the above method embodiments, which will not be described here. Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the devices and modules described above can refer to the corresponding process description in the above method embodiments, which will not be described here.

[0062] It should be noted that, according to the needs of implementation, each component / step described in the embodiments of the present application can be split into more components / steps, or two or more components / steps or part of the operations of the components / steps can be combined into a new component / step, to achieve the purpose of the embodiments of the present application.

[0063] The above-described methods according to embodiments of the application can be implemented in hardware, firmware, or software, or any combination thereof, and can be stored in or implemented with the aid of one or more data storage media (e.g., removable storage media like CD-ROM, floppy disks, hard disks, optical disks, or memory such as ROM, RAM, flash memory, etc.), or implemented with the aid of one or more processors (e.g., general-purpose, special-purpose or programmable hardware) of a processing system, or any combination thereof. The methods described herein can be stored as software or computer code on a non-transitory recording medium which can be any available media that can be read by a general- purpose or special-purpose processing system, such as a computer, processor, or microprocessor, or any programmable hardware. It is understood that the methods described herein can be implemented with the aid of software stored on a recording medium, which is read and executed by a processing system to perform the tasks specified in the software. Furthermore, a general- purpose computer, special-purpose computer, or programmable hardware can be used to implement the described methods, wherein software code is stored in a non-transitory recording medium which is accessed by the computer, processor, or hardware to cause the computer, processor, or hardware to perform the methods described herein.

[0064] Those skilled in the art can realize that the units and method steps of each example described in connection with the embodiments disclosed herein can be realized by electronic hardware, or a combination of computer software and electronic hardware. Whether the functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the embodiments of the present application.

[0065] The above embodiments are only used to illustrate but not limit the embodiments of the present application, and ordinary skilled in the art can make various changes and modifications without departing from the spirit and scope of the embodiments of the present application, therefore all equivalent technical solutions belong to the scope of the embodiments of the present application, the patent protection scope of the embodiments of the present application should be defined by the claims.

Claims

1. A full-slope erosion and sediment transport simulation method that integrates soil properties and land use, characterized by: include: Collect basic erosion and sediment transport data corresponding to the entire slope surface of ridge slopes and valley slopes in the basin, including soil data, topography data, land use data, vegetation data and water and sediment data; Automatically determine the erosion environment type of the entire slope based on erosion and sediment transport basic data, terrain conditions, and land use type boundary conditions, and obtain initial parameters to be optimized that affect erosion based on the erosion environment type; The initial parameters to be optimized are updated based on the cuckoo search algorithm until the optimal parameters are obtained; By changing the slope combination and slope length in the terrain conditions, as well as the land use type and vegetation coverage combination in the land use type boundary conditions, a variety of different scenarios are obtained, and based on the various different scenarios and the optimal parameters, the sediment production and transport evolution process of the entire slope is simulated to obtain the sediment transport rate of the entire slope under the various different scenarios.

2. The method according to claim 1, characterized in that The obtaining of initial parameters to be optimized that affect erosion based on the erosion environment type includes: Setting simulated terrain conditions and simulated land use type boundary conditions based on the erosion environment type; A parameter type and a preset parameter range are obtained based on the simulated land use type boundary condition, and a group of parameters are randomly generated in the preset parameter range according to the parameter type as the initial parameters to be optimized.

3. The method according to claim 1, characterized in that The updating of the initial parameters to be optimized based on the cuckoo search algorithm until the optimal parameters are obtained includes: Obtaining a current simulated single-width sediment transport rate based on the initial parameters to be optimized and a single-width sediment transport rate simulation formula, and using the initial parameters to be optimized as optimal parameters when a comparison result between the simulated single-width sediment transport rate and the measured single-width sediment transport rate meets a preset accuracy judgment index condition; When the comparison result of the simulated single-width sediment transport rate and the measured single-width sediment transport rate does not meet the preset accuracy judgment index condition, the initial parameter to be optimized is iteratively updated again until the comparison result of the simulated single-width sediment transport rate and the measured single-width sediment transport rate corresponding to the updated parameter to be optimized meets the preset accuracy judgment index condition, and the updated parameter to be optimized is used as the optimal parameter.

4. The method according to claim 1, wherein The automatic determination of the erosion environment type of the entire slope based on the erosion and sediment transport basic data, terrain conditions, and land use type boundary conditions includes: The terrain data, the land use data, the vegetation data, the water and sand data, the terrain conditions and the land use type boundary conditions are sorted into a comma-separated value file in the same directory as the preset discrimination program, and the erosion environment type is obtained in combination with the soil data input in the preset discrimination program.

5. The method according to claim 3, characterized in that The simulating the sediment production and transport evolution process of the entire slope based on the multiple different scenarios and the optimal parameters to obtain the sediment transport rate of the entire slope under the multiple different scenarios includes: Substituting the optimal parameters and the corresponding parameter values ​​under the multiple different scenarios into the single-width sediment transport rate simulation formula to obtain the single-width sediment transport rates under the multiple different scenarios; The sediment transport rates of the entire slope under the various scenarios are obtained based on the single-width sediment transport rates under the various scenarios.

6. The method according to any one of claims 1 to 5, characterized in that The erosion environment types include a single slope type, a composite slope type with the same land use in the upper and lower sections, and a composite slope type with different land use in the upper and lower sections. The composite slope type is composed of the ridge slope and the valley slope. The terrain conditions include the slope combination and the slope length. The land use type boundary conditions corresponding to the composite slope type with the same land use in the upper and lower sections include the upper section being a combination of bare land and bare land, and the lower section being a combination of grassland and grassland. The land use type boundary conditions corresponding to the composite slope type with different land use in the upper and lower sections include the upper section being a combination of farmland and grassland, and the lower section being a combination of forestland and grassland.

Images