Soil environment monitoring method applied to ecological restoration of refuse landfill

Abstract

The invention discloses a soil environment monitoring method applied to ecological restoration of a refuse landfill, and relates to the technical field of soil monitoring. The method comprises the following steps: acquiring environment monitoring data in each soil monitoring area of the refuse landfill; based on the correction coefficient and the initial heavy metal curing strength, obtaining the corrected heavy metal curing strength of each soil monitoring area; and based on the corrected heavy metal solidification strength, obtaining the heavy metal content of each soil monitoring area at different depths. According to the method, the collected heavy metal content is corrected and expressed by analyzing the characteristics (namely correcting the heavy metal solidification strength) of the diffusion influence of the substances such as heavy metal substance pollution when the humus which is filled at different depths in the refuse landfill continuously accumulatively covers and fills the refuse, so that the particularity of the original heavy metal content is reduced, and the content of the heavy metal in the refuse landfill is improved. And the accuracy of monitoring the soil environment in the refuse landfill is improved.

Description

Technical Field

[0001] This application relates to the field of soil monitoring technology, specifically to a soil environmental monitoring method applied to the ecological restoration of landfills. Background Technology

[0002] Landfill ecological restoration refers to the maintenance, remediation, and restoration of landfills that have been used or abandoned. This process aims to reduce the negative environmental impact of landfills, restore the site's natural ecological functions, and minimize the potential harm of pollutants to the surrounding environment and human health. Because landfills primarily treat waste by burying it in soil, the soil becomes contaminated with waste, forming humus. This humus contains high levels of various heavy metals, posing significant health risks to humans and the environment. Therefore, when carrying out ecological restoration of used or abandoned landfills, it is necessary to conduct soil environmental monitoring of the humus within the landfill to ensure that the restored site's impact on the environment and community is minimized and that long-term environmental sustainability is achieved.

[0003] Because the types and quantities of waste buried in different areas of a landfill vary, and the topography of the landfill itself also plays a role, the types of waste buried in the soil differ in different areas, leading to variations in the diffusion pathways of heavy metals in these areas. Current technologies typically select a specific area within the landfill and test the heavy metal content in samples of humus from that area to determine the quality of the humus. However, these methods do not consider the impact of landfill depth on heavy metal infiltration. These factors contribute to the technical problem that existing monitoring methods cannot accurately monitor the heavy metal content in landfill soil. Summary of the Invention

[0004] The purpose of this application is to provide a soil environmental monitoring method for ecological restoration of landfills, in order to solve the technical problem that existing monitoring methods cannot accurately monitor the heavy metal content in the soil of landfills.

[0005] To achieve the above objectives, this application provides the following technical solution: This application proposes a technical solution for a soil environmental monitoring method applied to the ecological restoration of landfills. This soil environmental monitoring method for landfill ecological restoration includes: Acquire environmental monitoring data for each soil monitoring area of ​​the landfill; the environmental monitoring data shall include at least the heavy metal content at different depths and the permeability coefficient under different confining pressures in each soil monitoring area. Based on various environmental monitoring data, the initial heavy metal solidification intensity of humus at different depths was obtained; the initial heavy metal solidification intensity was used to characterize the heavy metal solidification capacity of humus at the corresponding depth. Based on the initial heavy metal solidification intensity of humus at different depths, correction coefficients are obtained for each soil monitoring area; the correction coefficients are used at least to characterize the extent to which the heavy metal content of the corresponding soil monitoring area is affected by the heavy metal content of its adjacent soil monitoring areas. Based on the correction coefficient and the initial heavy metal solidification intensity, the corrected heavy metal solidification intensity of each soil monitoring area is obtained; Based on the modified heavy metal solidification intensity, the heavy metal content at different depths in each soil monitoring area was obtained.

[0006] As a specific solution in this application, the step of obtaining environmental monitoring data in each soil monitoring area of ​​the landfill includes: The landfill was divided into multiple first soil monitoring zones; Based on each first soil monitoring area, a second soil monitoring area is obtained; the second soil monitoring area is any area in each first soil monitoring area where no environmental monitoring data was obtained. Based on the second soil monitoring area, soil samples at different depths were obtained; Environmental monitoring data for the second soil monitoring area were obtained based on each soil sample.

[0007] As a specific solution in this application, the step of obtaining the initial heavy metal solidification strength of humus at different depths based on various environmental monitoring data includes: Based on the permeability coefficient under different confining pressures, the geological similarity under the same confining pressure is obtained; the geological similarity under the same confining pressure is at least used to characterize the similarity of the permeability coefficients of humus soils under the same confining pressure. Based on the permeability coefficient under different confining pressures, the support for heavy metal movement at the same depth is obtained; the support for heavy metal movement at the same depth is at least used to characterize the magnitude of the support capacity of a single-depth humus soil for heavy metal movement within its own confining pressure variation range. Based on the geological similarity under the same confining pressure and the support for heavy metal movement at the same depth, the comprehensive soil heavy metal diffusivity is obtained; the comprehensive soil heavy metal diffusivity is at least used to characterize the ability of a single depth of humus soil to diffuse heavy metals into adjacent depths of humus soil. Based on the comprehensive soil heavy metal diffusivity, the initial heavy metal solidification intensity is obtained; the initial heavy metal solidification intensity is used to characterize the ability of humus soil at different depths to fix heavy metals.

[0008] As a specific solution in this application, the method of obtaining geological similarity under the same confining pressure based on the permeability coefficient under different confining pressures includes: A two-dimensional coordinate system is established based on the permeability coefficient under different confining pressures, and multiple coordinate points are obtained; the horizontal axis of the two-dimensional coordinate system is the confining pressure, and the vertical axis is the permeability coefficient. Fit the data to each coordinate point to obtain the permeation curve corresponding to each depth; Based on each permeation curve, a first difference value and a second difference value are obtained; the first difference value is the difference in permeation slope at different depths under the same confining pressure; the second difference value is the difference in permeability coefficient at different depths under the same confining pressure. Based on the first difference value and the second difference value, the geological similarity under the same confining pressure is obtained.

[0009] As a specific solution in this application, the method of obtaining the support for heavy metal movement at the same depth based on the permeability coefficient under different confining pressures includes: For any depth within the same waste soil monitoring area, extract all permeability coefficients and corresponding confining pressure data within the permeability curve area corresponding to the empirical confining pressure range at that depth; the empirical confining pressure range is preset. Based on each permeability coefficient and each confining pressure, multiple sets of third and fourth difference values ​​are obtained; the third difference value is the difference between adjacent permeability coefficients among each permeability coefficient; the fourth difference value is the difference between adjacent confining pressures among each confining pressure. Based on the third and fourth difference values ​​of each group, the support for heavy metal movement at the same depth is obtained.

[0010] As a specific solution in this application, the step of obtaining comprehensive soil heavy metal diffusivity based on the geological similarity under the same confining pressure and the support for heavy metal movement at the same depth includes: For any depth within the same waste soil monitoring area, analyze the empirical confining pressure range of that depth and adjacent depths to obtain the overlapping range; Based on the overlapping range, multiple first confining pressures are obtained; the first confining pressure is any confining pressure within the overlapping range corresponding to the depth of the permeation curve. Based on the geological similarity of each first confining pressure and the support for heavy metal movement at the same depth, the comprehensive soil heavy metal diffusivity is obtained.

[0011] As a specific solution in this application, obtaining the initial heavy metal solidification strength based on the comprehensive soil heavy metal diffusivity includes: Based on the heavy metal content at different depths within each soil monitoring area, the current state of heavy metal content in humus at each depth is obtained; the current state of heavy metal content is used at least to characterize the magnitude of heavy metal content in humus at the corresponding depth. The content of each heavy metal is arranged in order of depth from shallow to deep to obtain a depth sequence; The initial heavy metal solidification intensity is obtained based on the depth sequence and the comprehensive soil heavy metal diffusivity.

[0012] As a specific solution in this application, the calculation formula for obtaining the initial heavy metal solidification intensity based on the depth sequence and the comprehensive soil heavy metal diffusivity is as follows: in, Indicates the initial heavy metal curing strength; This indicates the number of sequence values ​​in the depth sequence that precede the sequence value at the current depth. This indicates the overall soil heavy metal diffusivity at the current depth; This represents the overall soil heavy metal diffusivity at the m-th depth; This represents the current state of heavy metal content at the m-th depth; \mathrm{n}\mathrm{o}\mathrm{r}\mathrm{m}[\, \, ] This represents a normalization function used to normalize the sequence values ​​within the parentheses to the range [0, 1]; \mathrm{e}\mathrm{x}\mathrm{p}[\, \, ] This represents an exponential function with the natural constant e as its base.

[0013] As a specific solution in this application, the step of obtaining the correction coefficient corresponding to each soil monitoring area based on the initial heavy metal solidification intensity of humus at different depths includes: Based on each soil monitoring area, a third soil monitoring area is obtained; the third soil monitoring area is any soil monitoring area in each soil monitoring area for which no correction coefficient has been obtained. Based on the third soil monitoring area, a fourth soil monitoring area was obtained; the fourth soil monitoring area is adjacent to the third soil monitoring area, and no value was obtained for the degree of horizontal spread between the third soil monitoring area and the fourth soil monitoring area; Based on the third and fourth soil monitoring areas, multiple sampling distances and multiple fifth difference values ​​are obtained; the sampling distance is the distance between sampling points at the same depth in the third and fourth soil monitoring areas; the fifth difference value is the difference in the initial heavy metal solidification intensity at the same depth in the third and fourth soil monitoring areas. Based on each sampling distance and each fifth difference value, the horizontal spread impact value at different depths of the third soil monitoring area is obtained; the horizontal spread impact value is at least used to characterize the extent to which heavy metal substances in the third soil monitoring area can easily spread horizontally to the fourth soil monitoring area. Based on the horizontal spread impact value, the correction coefficient for the third soil monitoring area is obtained.

[0014] As a specific solution in this application, obtaining the correction coefficient for the third soil monitoring area based on the horizontal spread impact value includes: Based on the third soil monitoring area, multiple fifth soil monitoring areas are obtained; the fifth soil monitoring area is any soil monitoring area adjacent to the third soil monitoring area in each soil monitoring area, and the degree of horizontal spread influence between the third soil monitoring area and the fifth soil monitoring area has been obtained; Based on each fifth soil monitoring area, multiple horizontal spread impact values ​​were obtained; each horizontal spread impact value corresponds one-to-one with the fifth soil monitoring area. Based on the degree of spread of influence at each level, the correction coefficient for the third soil monitoring area is obtained.

[0015] Compared with the prior art, the beneficial effects of this application are: This application analyzes the characteristics of the diffusion effect of heavy metal pollution on humus at different depths in landfills as it is continuously accumulated and covered on top of the waste (i.e., correcting the solidification intensity of heavy metals). It then corrects the representation of the collected heavy metal content, reducing the specificity of the original heavy metal content (i.e., the randomness and randomness of the corresponding data caused by local sampling of humus samples), and improving the universality of the information representation of heavy metal content in landfills, thereby improving the accuracy of soil environmental monitoring in landfills. Attached Figure Description

[0016] Figure 1 This is a schematic flowchart of a soil environmental monitoring method for ecological restoration of landfills proposed in this application. Detailed Implementation

[0017] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0018] The terms "first," "second," etc., in the specification and accompanying drawings of this application are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. For example, the first soil monitoring area and the second soil monitoring area mentioned below belong to different soil monitoring areas. It should be understood that such names can be used interchangeably where appropriate so that the embodiments described herein can be implemented in an order other than that illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or modules is not necessarily limited to those steps or modules explicitly listed, but may include other steps or modules not explicitly listed or inherent to these processes, methods, products, or devices. The division of modules in the embodiments of this application is merely a logical division. In actual applications, there may be other division methods. For example, multiple modules may be combined into or integrated into another system, or some features may be ignored or not performed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be through some interface, and the indirect coupling or communication connection between modules may be electrical or other similar forms. None of these are limited in the embodiments of this application. Furthermore, the modules or sub-modules described as separate components may or may not be physically separated, may or may not be physical modules, or may be distributed among multiple circuit modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the embodiments of this application.

[0019] To address the technical problem mentioned in the background art that existing monitoring methods cannot accurately monitor the heavy metal content in landfill soil, this application proposes an embodiment of a soil environmental monitoring method applied to the ecological restoration of landfills. For example... Figure 1 As shown, the soil environmental monitoring method applied to the ecological restoration of landfills includes steps S100 to S500.

[0020] Step S100: Obtain environmental monitoring data for each soil monitoring area of ​​the landfill.

[0021] In this embodiment, the environmental monitoring data includes at least the heavy metal content at different depths and the permeability coefficient under different confining pressures within each soil monitoring area.

[0022] In the embodiments of this application, environmental monitoring data for each soil monitoring area of ​​the landfill can be obtained in any reasonable manner. For example, if the landfill has previously conducted soil monitoring and the historical data (e.g., heavy metal content at different depths, permeability coefficients under various confining pressures) for each area (e.g., the first soil monitoring area and the second soil monitoring area hereinafter) are stored on a local hard drive (which can be a monitoring terminal hard drive or a laboratory server hard drive, etc.) categorized by "area-depth-detection time", the historical files on the hard drive can be directly read using data management software (e.g., Excel or a professional environmental monitoring data system, etc.). For example, in the folder "2024 Landfill Monitoring Data" on the hard drive, the corresponding heavy metal detection reports and permeability coefficient records can be retrieved by subdirectories such as "first area-1m depth" and "second area-9m depth", quickly obtaining the archived monitoring data without the need for repeated sampling and testing. Alternatively, for massive monitoring data from multiple areas of a large landfill, a "regional data index table" can be created on the hard drive, marking the storage path of each soil monitoring area's data (e.g., "D drive\Landfill Monitoring\Third Area\202506 Sampling-5m Depth-Heavy Metal Data.xlsx"), data type (e.g., heavy metals or permeability coefficient), and detection parameters (e.g., sampling depth or confining pressure range). During data retrieval, the index table quickly locates the target area's data file, allowing direct access to the corresponding values. Simultaneously, backup data from adjacent areas can be retrieved to verify data consistency and improve data retrieval efficiency.

[0023] In a specific embodiment of this application, step S100 involves acquiring environmental monitoring data within each soil monitoring area of ​​the landfill, including steps S110 to S140.

[0024] Step S110: Divide the landfill into multiple first soil monitoring zones.

[0025] In the embodiments of this application, the landfill can be divided into multiple first soil monitoring areas in any reasonable manner. For example, if the landfill has regular terrain (e.g., rectangular or square) and the landfill material is evenly distributed, rectangular grids can be divided according to preset sizes (e.g., 50m×50m or 100m×100m, which can be adjusted according to the total area of ​​the landfill) based on the landfill boundary. Each grid is a first soil monitoring area. For example, a landfill with a total area of ​​1000m×800m can be divided into 100m×100m areas to obtain 80 equal-sized first soil monitoring areas, ensuring that the area of ​​each area is consistent and reducing sampling deviations caused by differences in area size. This is suitable for landfills with no obvious differences in pollution in the initial stage. Alternatively, if the landfill is circular or nearly circular, multiple concentric circles can be drawn based on the center point of the landfill (i.e., the center of the circle) with preset radii (e.g., 20m, 40m, or 60m). Radial lines can then be drawn from the center point at angles of 30°, 45°, or 60°, dividing the concentric circles into multiple sector-shaped areas. Each sector represents a primary soil monitoring area. For example, for a circular landfill with a radius of 100m, drawing five concentric circles with a 20m radius difference, and then dividing them with 45° radial lines, yields 40 sector-shaped primary soil monitoring areas. This approach takes into account the gradient differences from the center to the edge of the landfill and is suitable for applications requiring attention to the radial diffusion of pollutants. If the landfill is constructed in phases using a "unit-based" approach (e.g., a first-phase landfill unit and a second-phase landfill unit), then each completed landfill unit can be directly designated as a primary soil monitoring area. For example, if a landfill is constructed in three phases, with each phase containing five independent landfill units, then the landfill can be divided into 15 primary soil monitoring zones. If the landfill is divided into zones based on waste type (e.g., municipal solid waste zone, food waste zone, and hazardous waste storage zone), then each waste type's landfill zone can be designated as a primary soil monitoring zone. For instance, the municipal solid waste zone could be divided into 100m x 100m grids, the food waste zone (prone to high-concentration leachate) into 50m x 50m grids, and the hazardous waste storage zone (high pollution risk) into 20m x 20m grids; each grid would constitute a primary soil monitoring zone. For example, a landfill containing three municipal solid waste zones, two food waste zones, and one hazardous waste storage zone could be divided into 25 primary soil monitoring zones, enabling differentiated monitoring of zones with different pollution risks, with a focus on the heavy metal and organic matter content in high-risk areas.

[0026] Step S120: Based on each of the first soil monitoring areas, obtain the second soil monitoring area.

[0027] In this embodiment, the second soil monitoring area is any area in each of the first soil monitoring areas where no environmental monitoring data has been obtained. That is, in this embodiment, the method for obtaining environmental monitoring data for any soil monitoring area in each of the first soil monitoring areas is the same as the method for obtaining environmental monitoring data for the second soil monitoring area.

[0028] Step S130: Based on the second soil monitoring area, obtain soil samples at different depths.

[0029] In this embodiment, the depth of the soil sample can be preset. In subsequent embodiments, soil samples with depths of 1m, 5m, and 9m are used to illustrate various embodiments of this application, which does not mean that the soil environmental monitoring method for landfill ecological restoration proposed in this application is only applicable to the acquisition of environmental monitoring data from soil samples with depths of 1m, 5m, or 9m.

[0030] Step S140: Based on each soil sample, obtain environmental monitoring data for the second soil monitoring area.

[0031] In this embodiment, soil samples at the same depth corresponding to the second soil monitoring area are divided into two parts. One part is used to detect the content of various heavy metals (e.g., lead, cadmium, and mercury), and the other part is tested using a flexible wall permeability tester, starting from 0 and with a step size of 1, sequentially controlling the confining pressure to 80 kPa (or other preset values, such as 70 kPa or 60 kPa) and stopping the pressure application, and obtaining the permeability coefficient under each iteration of the confining pressure. The landfill contains multiple soil monitoring areas (i.e., various first soil monitoring areas), each soil monitoring area contains multiple depths, each depth contains multiple heavy metal content data and multiple permeability coefficients, and each permeability coefficient corresponds to a confining pressure.

[0032] In this embodiment, confining pressure represents the compressive strength of the humus soil, and permeability coefficient represents the strength of the diffusion of minute pollutants such as heavy metals outward from the humus soil under the corresponding compressive strength.

[0033] Step S200: Based on various environmental monitoring data, obtain the initial heavy metal solidification intensity of humus at different depths.

[0034] In this embodiment, the initial heavy metal solidification intensity is used at least to characterize the heavy metal solidification capacity of humus at the corresponding depth.

[0035] It is important to note that in actual soil monitoring, landfills typically cover large areas. Due to time constraints, existing methods often involve dividing the landfill into different zones (i.e., multiple primary soil monitoring zones mentioned above) for soil sampling to complete the monitoring and assessment. However, this approach has significant drawbacks. It ignores the technical problem of "locality" in the monitoring data caused by the migration of heavy metals within the landfill. That is, the sampling data only reflects the soil monitoring situation in a specific, small area and cannot comprehensively represent the overall soil monitoring situation of the landfill. The uneven distribution of waste types and volumes within a landfill, along with variations in the site's topography, leads to different degrees of soil compaction in different areas, resulting in significant differences in the heavy metal diffusion paths formed in each area. Therefore, to improve the accuracy of subsequent soil monitoring results, the first step is to scientifically divide the landfill into zones and clearly define different experimental analysis units. For the same waste soil monitoring area, the changes in the characteristics of the internal humus soil also need to be considered. The humus soil originally covering the waste, due to long-term contact and compression with the waste, allows harmful heavy metals contained in the waste to diffuse upwards from the lower soil layer in direct contact with the waste, and the degree of diffusion gradually changes over time. Simultaneously, the continuous accumulation pressure within the soil further contributes, ultimately causing differentiated geological structures in the humus soil at different depths within the same monitoring area. In real-world environments, the geological structure of humus soil at different depths directly affects the diffusion path of heavy metals. Based on this, in a specific embodiment of this application, step S200, based on various environmental monitoring data, obtains the initial heavy metal consolidation intensity of the humus soil at different depths, including steps S210 to S240.

[0036] Step S210: Based on the permeability coefficient under different confining pressures, obtain the geological similarity under the same confining pressure.

[0037] In this embodiment, the geological similarity under the same confining pressure is used at least to characterize the similarity of the permeability coefficients of humus soils under the same confining pressure.

[0038] It is important to note that humus at different depths experiences varying confining pressures (i.e., landfill pressure) during landfilling, resulting in significantly different pore structures and permeability channels (e.g., shallow humus has larger pores and higher permeability coefficients, while deeper humus is denser due to compression, resulting in smaller pores and lower permeability coefficients). Directly comparing the solidification effects at different depths would render the data meaningless due to the different geological foundations. In this embodiment, geological similarity under the same confining pressure essentially quantifies whether the pore distribution and diffusion paths of humus at different depths are similar under the same compressive strength. Higher similarity indicates stronger commonality in the geological structure between depths, reducing errors caused by structural differences when calculating the initial heavy metal solidification strength and providing a basis for direct comparison of solidification capabilities at different depths. Based on this, step S210, obtaining geological similarity under the same confining pressure based on the permeability coefficients under different confining pressures, includes steps S211 to S214.

[0039] Step S211: Establish a two-dimensional coordinate system based on the permeability coefficient under different confining pressures and obtain multiple coordinate points.

[0040] In this embodiment, the horizontal axis of the two-dimensional coordinate system represents the confining pressure, and the vertical axis represents the permeability coefficient.

[0041] Step S212: Fit the coordinates to obtain the permeation curves corresponding to each depth.

[0042] It is important to understand that fitting multiple coordinate points to obtain the corresponding curve (i.e., the penetration curve) is a mature technology in the field of computer science, which will not be elaborated on here.

[0043] Step S213: Based on each permeation curve, obtain the first difference value and the second difference value.

[0044] In this embodiment, the first difference value is the difference in permeability slope at different depths under the same confining pressure. The second difference value is the difference in permeability coefficient at different depths under the same confining pressure.

[0045] Step S214: Based on the first difference value and the second difference value, obtain the geological similarity under the same confining pressure.

[0046] In this embodiment, the geological similarity under the same confining pressure can be obtained based on the first difference value and the second difference value in any reasonable manner, as long as the geological similarity under the same confining pressure is negatively correlated with the magnitude of the first difference value and the second difference value. For example, in step S214, the calculation formula for obtaining the geological similarity under the same confining pressure based on the first difference value and the second difference value is as follows: in, Indicates geological similarity to confining pressure; Indicates the first difference value; 'a' represents the second difference value; 'a' represents the zero-prevention coefficient, which is any positive number close to 0. For example, the zero-prevention coefficient can be 0.01 or 0.001, etc. This indicates the absolute value; This represents a normalization function used to normalize the values ​​within the parentheses to the range [0, 1].

[0047] Alternatively, in step S214, based on the first difference value and the second difference value, the calculation formula for the geological similarity under the same confining pressure is obtained as follows: in, Indicates geological similarity to confining pressure; Indicates the first difference value; Indicates the second difference value; This represents an exponential function with the natural constant e as its base. This indicates that the absolute value is being calculated.

[0048] In this embodiment, the smaller the first and second difference values, the greater the geological similarity under the same confining pressure. This means that under the same confining pressure, the pore distribution and permeation channels (i.e., heavy metal diffusion paths) of humus at different depths are more similar, resulting in a more consistent response to external compression and a more uniform physical structure. This reduces the deviation of depth differences from the heavy metal diffusion pattern. Conversely, the larger the first and second difference values, the smaller the geological similarity under the same confining pressure. This means that under the same confining pressure, the differences in pore size and permeation channels between humus at different depths are more significant, resulting in greater differences in compression response and more pronounced differences in physical structure. This increases the deviation of the heavy metal diffusion pattern between depths. In other words, geological similarity under the same confining pressure reflects, to some extent, the similarity of the upward diffusion paths of heavy metals that can be reserved in humus at different depths under depositional pressure. A greater geological similarity under the same confining pressure indicates that the pore spacing within the corresponding humus is more similar, reflecting a more uniform geological structure. When these humus soils are piled up and landfilled according to the actual distribution depth in the landfill, the compression state of the humus soil within the soil is usually not significantly different due to the relatively uniform geological structure and the relatively uniform degree of response to external compression. Therefore, by analyzing the impact of the landfill pressure on the cumulatively landfilled humus soil on heavy metal diffusion based on the similarity of the confining geological conditions, the initial heavy metal solidification strength of the humus soil at different depths can be calculated.

[0049] Step S220: Based on the permeability coefficient under different confining pressures, obtain the support for heavy metal movement at the same depth.

[0050] In this embodiment, the same-depth heavy metal movement support is at least used to characterize the magnitude of the support capacity of a single-depth humus soil for heavy metal movement within its own confining pressure variation range.

[0051] It is important to note that heavy metals, as minute and difficult-to-observe harmful abiotic substances, do not diffuse and spread autonomously within humus. Their diffusion within humus is entirely dependent on external forces. In decommissioned or abandoned landfills, there are typically no vehicles passing through or plants with strong root absorption capabilities present. Therefore, in this scenario, the humus primarily relies on the internal pressure naturally generated during landfilling as its driving force, while the voids left by the geological conditions within the humus itself serve as diffusion channels supporting the movement of heavy metals. The changes in these two states together support the fixed retention of heavy metals in the humus at the same depth. Based on this, step S220, based on the permeability coefficient under different confining pressures, obtains the support for heavy metal movement at the same depth, including steps S221 to S223.

[0052] Step S221: For any depth within the same waste soil monitoring area, extract all permeability coefficients and corresponding confining pressure data within the permeability curve area corresponding to the empirical confining pressure range at that depth.

[0053] It should be noted that, under normal circumstances, in landfills that have been used up or abandoned, the compressive strength of humus at the same depth under the influence of the soil layers above and below generally fluctuates within a certain range. Therefore, the corresponding confining pressure range (i.e., the empirical confining pressure range) can be obtained first through empirical preset. For example, in one embodiment of this application, the empirical confining pressure ranges corresponding to humus at depths of 1m, 5m, and 9m can be [10kPa, 20kPa), [15kPa, 50kPa), and [40kPa, 80kPa], respectively.

[0054] Step S222: Based on each permeability coefficient and each confining pressure, obtain multiple sets of third and fourth difference values.

[0055] In this embodiment, the third difference value is the difference between adjacent permeability coefficients. The fourth difference value is the difference between adjacent confining pressures.

[0056] In this embodiment, the third difference value directly reflects the fluctuation range of the permeability coefficient of humus at the same depth when the confining pressure changes. The permeability coefficient is a "channel indicator" for heavy metal diffusion. If the third difference value is larger, it indicates that even a small change in confining pressure will lead to a significant change in the heavy metal diffusion channels (e.g., pores or fissures). The fourth difference value reflects the magnitude of the change in confining pressure and is used to anchor the "pressure background" of the permeability coefficient change. The combination of the two can eliminate the interference of "abnormal fluctuations in permeability coefficient caused by drastic changes in confining pressure" and accurately capture the influence of the geological characteristics of humus itself on its permeability.

[0057] Step S223: Based on the third and fourth difference values ​​of each group, obtain the support for the movement of heavy metals at the same depth.

[0058] In this embodiment, the support for heavy metal movement at the same depth can be obtained based on each group of third and fourth difference values ​​using any reasonable method. For example, the support for heavy metal movement at the same depth can be the sum of each group of third and fourth difference values, or, in step S223, the calculation formula for obtaining the support for heavy metal movement at the same depth based on each group of third and fourth difference values ​​can be as follows: in, Indicates the support for heavy metal movement at the same depth; Indicates the number of groups with the third and fourth difference values; This represents the third difference value of the i-th group; This represents the fourth difference value in the i-th group. In this embodiment, if the support for heavy metal movement at the same depth is greater, it indicates that within the corresponding confining pressure range, pressure changes in the humus soil at that depth more easily drive heavy metals to move through the infiltration channels, i.e., the stronger the support for heavy metal movement.

[0059] Step S230: Based on the geological similarity of the same confining pressure and the support of heavy metal movement at the same depth, obtain the comprehensive soil heavy metal diffusivity.

[0060] In this embodiment, the comprehensive soil heavy metal diffusivity is used at least to characterize the ability of a single depth of humus to diffuse heavy metals into adjacent depths of humus.

[0061] It is important to note that, ideally, as the depth of humus in a landfill increases, the landfill pressure (i.e., confining pressure) gradually increases, meaning there is no situation where adjacent depths have the same landfill pressure. However, in reality, due to the diffusion of heavy metals and other substances between humus at different depths, the weight of humus at corresponding depths may change, potentially leading to situations where humus at different depths has the same landfill pressure. Therefore, in areas where the confining pressure ranges overlap between different depths, based on the support for heavy metal movement at the same depth and combined with the geological similarity of the same confining pressure, we can analyze the ability of other depths to accept heavy metals flowing out of the humus at that depth when it diffuses to other depths (i.e., the comprehensive soil heavy metal diffusivity). Based on this, step S230, based on the geological similarity of the same confining pressure and the support for heavy metal movement at the same depth, obtains the comprehensive soil heavy metal diffusivity, including steps S231 to S233.

[0062] Step S231: For any depth within the same waste soil monitoring area, analyze the empirical confining pressure range of that depth and adjacent depths to obtain the overlapping range.

[0063] In this embodiment, the overlapping range refers to the confining pressure interval that overlaps between the empirical confining pressure range of a certain target depth and the empirical confining pressure range of its adjacent depth within the same waste soil monitoring area. For example, the empirical confining pressure range of humus soil at a depth of 1m is [10kPa, 20kPa); the empirical confining pressure range of humus soil at a depth of 5m is [15kPa, 50kPa]. Since the adjacent depth of 5m is 1m, the overlapping range between the two is [15kPa, 20kPa].

[0064] Step S232: Based on the overlapping range, obtain multiple first confining pressures.

[0065] In this embodiment, the first confining pressure is any confining pressure within the overlapping range where the permeation curve corresponding to the depth is located.

[0066] Step S233: Based on the geological similarity of each first confining pressure and the support for heavy metal movement at the same depth, obtain the comprehensive soil heavy metal diffusivity.

[0067] In this embodiment, the comprehensive soil heavy metal diffusivity can be obtained using any reasonable method based on the geological similarity of the same confining pressure corresponding to each first confining pressure and the support for heavy metal movement at the same depth. For example, in step S233, the calculation formula for obtaining the comprehensive soil heavy metal diffusivity based on the geological similarity of the same confining pressure corresponding to each first confining pressure and the support for heavy metal movement at the same depth can be as follows: in, Indicates the overall soil heavy metal diffusivity; Indicates the support for heavy metal movement at the same depth; Indicates the quantity of the first confining pressure; This indicates the geological similarity to the first confining pressure corresponding to the r-th confining pressure; This represents a normalization function used to normalize the values ​​within the parentheses to the range [0, 1].

[0068] Alternatively, in step S233, based on the geological similarity of each first confining pressure and the support for heavy metal movement at the same depth, the calculation formula for the comprehensive soil heavy metal diffusivity can be obtained as follows: in, Indicates the overall soil heavy metal diffusivity; Indicates the support for heavy metal movement at the same depth; Indicates the quantity of the first confining pressure; This indicates the geological similarity to the first confining pressure corresponding to the r-th confining pressure; This represents an exponential function with the natural constant e as its base.

[0069] In this embodiment, a higher overall soil heavy metal diffusivity indicates that the geological structure of the humus at that depth is more similar to that of the adjacent humus at the same depth within the overlapping confinement zone, and that the humus itself has a stronger ability to support the movement of heavy metals. The combined effect of these two factors makes it easier for heavy metals at that depth to break through the depth limit and diffuse into the humus at the adjacent depth. Conversely, a lower overall soil heavy metal diffusivity indicates a greater difference in geological structure between depths or a weaker ability to support the movement of heavy metals, meaning that the diffusion resistance of heavy metals in the depth direction is greater, making it difficult for them to spread to adjacent depths.

[0070] Step S240: Based on the comprehensive soil heavy metal diffusivity, obtain the initial heavy metal solidification intensity.

[0071] In this embodiment, the initial heavy metal consolidation strength is used to characterize the ability of humus soil at different depths to fix heavy metals. In this embodiment, the overall soil heavy metal diffusivity can be directly used as the initial heavy metal consolidation strength.

[0072] It is important to note that the comprehensive soil heavy metal diffusivity characterizes the ability of humus at different depths within the same landfill soil monitoring area to diffuse heavy metals. Although the comprehensive soil heavy metal diffusivity considers the impact of humus at adjacent depths on the heavy metal diffusion caused by the landfilling of humus at the current depth, in actual landfills that have been used up or abandoned, humus at different depths is continuously accumulated and landfilled. Therefore, it is necessary to further analyze the overall effect of the accumulated landfilling of humus at different depths on the diffusion of heavy metals in order to obtain the initial heavy metal solidification intensity of humus at different depths. Based on this, step S240, based on the comprehensive soil heavy metal diffusivity, obtains the initial heavy metal solidification intensity, including steps S241 to S243.

[0073] Step S241: Based on the heavy metal content at different depths within each soil monitoring area, obtain the current state of heavy metal content in the humus soil at each depth.

[0074] In this embodiment, the prevalence of heavy metal content is at least used to characterize the amount of heavy metals in the humus at the corresponding depth. For example, the prevalence of heavy metal content can be the average value of the heavy metal content in the humus at each depth. Alternatively, in step S241, based on the heavy metal content at different depths within each soil monitoring area, the calculation formula for the prevalence of heavy metal content in the humus at each depth is as follows: in, R represents the current level of heavy metal content; R represents the amount of heavy metal content data in humus at the current depth. This represents the heavy metal content value corresponding to the nth heavy metal content data; This represents a normalization function used to normalize the values ​​within the parentheses to the range [0, 1]. In this embodiment, a higher level of existing heavy metal content indicates that, under the condition that the movement of humus to other depths is not satisfied, more heavy metal substances are fixed in the humus at that depth; and under the condition that the movement of humus to other depths is satisfied, the amount of heavy metal change is greater.

[0075] Step S242: Arrange the current levels of each heavy metal content in order of depth from shallow to deep to obtain a depth sequence.

[0076] In this embodiment, taking any depth of humus soil as an example, the depths preceding that depth are counted, and these preceding depths are considered shallow humus soil. Each depth may correspond to multiple shallow humus soil layers, one shallow humus soil layer, or none at all. For example, assuming the depths corresponding to humus soil layers are 1m, 5m, and 9m, then for humus soil at a depth of 1m, there is no shallow humus soil layer; for humus soil at a depth of 5m, the humus soil at a depth of 1m is its shallow humus soil layer; and for humus soil at a depth of 9m, both the humus soil at a depth of 1m and the humus soil at a depth of 5m are its shallow humus soil layers.

[0077] Step S243: Based on the depth sequence and the comprehensive soil heavy metal diffusivity, obtain the initial heavy metal solidification intensity.

[0078] In this embodiment, the initial heavy metal consolidation intensity can be obtained based on the depth sequence and the comprehensive soil heavy metal diffusivity using any reasonable method. For example, the shallower the depth of the humus at a given depth (i.e., the smaller the sequence number in the depth sequence), the greater the heavy metal consolidation intensity corresponding to that depth of humus, meaning the heavy metals are less likely to be lost; conversely, the smaller the comprehensive soil heavy metal diffusivity corresponding to that depth of humus, the greater the heavy metal consolidation intensity, meaning the heavy metals are less likely to be lost. Based on this, in one embodiment of this application, the initial heavy metal consolidation intensity corresponding to a certain depth of humus can be negatively correlated with the sequence number in the depth sequence corresponding to that depth of humus and the comprehensive soil heavy metal diffusivity. It should be noted that the comprehensive soil heavy metal diffusivity characterizes the ability of humus at different depths within the same waste soil monitoring area to diffuse heavy metal substances. If the overall soil heavy metal diffusivity at a certain depth in a soil monitoring area is low (e.g., below 0.3), it indicates that the humus at that depth has a weak ability to diffuse heavy metals to adjacent depths, and the heavy metals do not flow significantly within the soil and remain stably. Conversely, if the overall soil heavy metal diffusivity at a certain depth in a soil monitoring area is high (e.g., greater than or equal to 0.3), it indicates that the humus at that depth has a strong ability to diffuse heavy metals to adjacent depths, and the heavy metals are easily mobile and remain unstable. Therefore, to avoid errors caused by a single formula and to improve the accuracy of the initial heavy metal solidification strength calculation, in another embodiment of this application, step S243, based on the depth sequence and the overall soil heavy metal diffusivity, obtains the following formula for calculating the initial heavy metal solidification strength: in, Indicates the initial heavy metal curing strength; This indicates the number of sequence values ​​in the depth sequence that precede the sequence value at the current depth. This indicates the overall soil heavy metal diffusivity at the current depth; This represents the overall soil heavy metal diffusivity at the m-th depth; This represents the current state of heavy metal content at the m-th depth; \mathrm{n}\mathrm{o}\mathrm{r}\mathrm{m}[\, \, ] This represents a normalization function used to normalize the sequence values ​​within the parentheses to the range [0, 1]; \mathrm{e}\mathrm{x}\mathrm{p}[\, \, ] This represents an exponential function with the natural constant e as its base.

[0079] To avoid the inability of a single formula to cover the differences between "low diffusion" and "high diffusion" scenarios (i.e., if a negative exponential formula is used for low diffusion scenarios (i.e., the overall soil heavy metal diffusivity is less than 0.3), the numerical compression may result in insufficient differentiation of the initial heavy metal solidification intensity; if normalized summation is used for high diffusion scenarios (i.e., the overall soil heavy metal diffusivity is greater than or equal to 0.3), the weakening effect of the overall soil heavy metal diffusivity on the initial heavy metal solidification intensity may be ignored), this embodiment calculates the initial heavy metal solidification intensity in segments. This ensures that the detailed differences in the initial heavy metal solidification intensity in low diffusion scenarios are identifiable, and that the weakening trend of the initial heavy metal solidification intensity in high diffusion scenarios is consistent with reality. Ultimately, the calculated result of the initial heavy metal solidification intensity is closer to the true fixation capacity of humus for heavy metals.

[0080] Step S300: Based on the initial heavy metal solidification intensity of humus at different depths, obtain the correction coefficients corresponding to each soil monitoring area.

[0081] In this embodiment, the correction coefficient is at least used to characterize the extent to which the heavy metal content of the corresponding soil monitoring area is affected by the heavy metal content of its adjacent soil monitoring areas.

[0082] It is important to note that the initial heavy metal consolidation intensity primarily considers the vertical impact of humus at different depths on heavy metal content within the same landfill monitoring area. In real-world environments, landfill sites that have been used or abandoned are often large in area. They are subject to the impact of heavy objects from garbage trucks and other external factors during previous use, and the inherent uneven geological composition of the landfill itself can lead to a certain influence on the diffusion of heavy metal pollution when different landfill monitoring areas meet. Therefore, step S300, based on the initial heavy metal consolidation intensity of humus at different depths, obtains the correction coefficients corresponding to each soil monitoring area, including steps S310 to S350.

[0083] Step S310: Based on each soil monitoring area, obtain the third soil monitoring area.

[0084] In this embodiment, the third soil monitoring area is any soil monitoring area among the various soil monitoring areas for which no correction coefficient has been obtained. That is, in this embodiment, the method for obtaining the correction coefficient for any soil monitoring area among the various soil monitoring areas is the same as the method for obtaining the correction coefficient for the third soil monitoring area.

[0085] Step S320: Based on the third soil monitoring area, obtain the fourth soil monitoring area.

[0086] In this embodiment, the fourth soil monitoring area is adjacent to the third soil monitoring area, and no horizontal spread impact value is obtained between the third soil monitoring area and the fourth soil monitoring area. That is, in this embodiment, the calculation method for the horizontal spread impact value between the third soil monitoring area and any adjacent soil monitoring area is the same as the calculation method for the horizontal spread impact value of the fourth soil monitoring area.

[0087] Step S330: Based on the third soil monitoring area and the fourth soil monitoring area, obtain multiple sampling distances and multiple fifth difference values.

[0088] In this embodiment, the sampling distance is the distance between sampling points at the same depth in the third soil monitoring area and the fourth soil monitoring area. The fifth difference value is the difference in the initial heavy metal solidification intensity at the same depth in the third soil monitoring area and the fourth soil monitoring area.

[0089] Step S340: Based on each sampling distance and each fifth difference value, obtain the horizontal spread impact value at different depths in the third soil monitoring area.

[0090] In this embodiment, the horizontal spread effect value is used to characterize at least the extent to which heavy metals in the third soil monitoring area can easily spread horizontally to the fourth soil monitoring area.

[0091] In this embodiment, the degree of horizontal spread impact at different depths in the third soil monitoring area can be obtained using any reasonable method based on each sampling distance and each fifth difference value. For example, the degree of horizontal spread impact at a certain depth in the third soil monitoring area is equal to the sum or product of the corresponding sampling distance and the fifth difference value at that depth.

[0092] In this embodiment, the smaller the horizontal spread influence value, the more consistent the humus state of the soil monitoring area (i.e., the third soil monitoring area) is with the corresponding adjacent soil monitoring area at the same depth, and the less likely the heavy metal substances contained therein are to spread horizontally to the corresponding adjacent soil monitoring area.

[0093] Step S350: Based on the horizontal spread impact value, obtain the correction coefficient for the third soil monitoring area.

[0094] In the embodiments of this application, any reasonable method can be used to obtain the correction coefficient of the third soil monitoring area based on the horizontal spread impact degree value. For example, the horizontal spread impact degree value can be directly used as the correction coefficient of the third soil monitoring area. Alternatively, step S350, obtaining the correction coefficient of the third soil monitoring area based on the horizontal spread impact degree value, includes steps S351 to S353.

[0095] Step S351: Based on the third soil monitoring area, obtain multiple fifth soil monitoring areas.

[0096] In this embodiment, the fifth soil monitoring area is any soil monitoring area adjacent to the third soil monitoring area among all soil monitoring areas, and the degree of horizontal spread influence between the third soil monitoring area and the fifth soil monitoring area has been obtained.

[0097] Step S352: Based on each fifth soil monitoring area, obtain multiple horizontal spread impact values.

[0098] In this embodiment, the values ​​of the degree of horizontal spread correspond one-to-one with the fifth soil monitoring area.

[0099] Step S353: Based on the degree of spread of influence at each level, obtain the correction coefficient for the third soil monitoring area.

[0100] In this embodiment, the greater the degree of spread influence at each level, the more necessary it is to correct the initial heavy metal solidification intensity corresponding to the third soil monitoring area; conversely, the smaller the degree of spread influence at each level, the less necessary it is to correct the initial heavy metal solidification intensity corresponding to the third soil monitoring area. In other words, in this embodiment, the correction coefficient is positively correlated with the degree of spread influence at each level.

[0101] In this embodiment, the correction coefficient for the third soil monitoring area can be obtained based on the horizontal spread impact values ​​using any reasonable method. For example, the correction coefficient can be the average of the horizontal spread impact values, or, in step S353, the calculation formula for the correction coefficient of the third soil monitoring area based on the horizontal spread impact values ​​is as follows: in, This represents the correction factor for the third soil monitoring area; This indicates the number of soil monitoring areas adjacent to the third soil monitoring area (i.e., the fifth soil monitoring area); This represents the degree of horizontal spread of influence between the third soil monitoring area and the y-th fifth soil monitoring area. This represents the normalization function, used to normalize the values ​​within the parentheses to the range [0, 1]. In this embodiment, a larger correction coefficient indicates a stronger influence of heavy metal spread from adjacent soil monitoring areas on the third soil monitoring area.

[0102] Step S400: Based on the correction coefficient and the initial heavy metal solidification intensity, obtain the corrected heavy metal solidification intensity for each soil monitoring area.

[0103] In this embodiment, the corrected heavy metal solidification intensity can be the product of the correction coefficient and the initial heavy metal solidification intensity, or, in step S400, based on the correction coefficient and the initial heavy metal solidification intensity, the calculation formula for the corrected heavy metal solidification intensity of each soil monitoring area is as follows: in, This indicates the correction for the curing strength of heavy metals; Indicates the correction factor; Indicates the initial heavy metal curing strength; This represents a normalization function used to normalize the values ​​within the parentheses to the range [0, 1].

[0104] Step S500: Based on the modified heavy metal solidification intensity, obtain the heavy metal content at different depths in each soil monitoring area.

[0105] In this embodiment, a higher correction for heavy metal solidification intensity indicates that the detected heavy metal content at the corresponding depth in the soil monitoring area is closer to its true heavy metal content. Based on this, in step S500, the calculation formula for the heavy metal content at different depths in each soil monitoring area, based on the correction for heavy metal solidification intensity, is as follows: in, This indicates the heavy metal content at any depth in any soil monitoring area. This indicates the heavy metal content detected at this depth in the soil monitoring area; This indicates the corrected heavy metal solidification intensity at this depth in the soil monitoring area.

[0106] The embodiments of the soil environmental monitoring method for landfill ecological restoration proposed in this application analyze the characteristics of the diffusion effect of heavy metal pollution on humus at different depths in the landfill as it is continuously accumulated and covered on the landfill (i.e., correcting the heavy metal solidification intensity). The collected heavy metal content is corrected and represented, reducing the specificity of the original heavy metal content (i.e., the randomness and randomness of the corresponding data caused by local sampling of humus samples). This improves the universality of the information representation of heavy metal content in landfills, thereby improving the accuracy of soil environmental monitoring in landfills.

[0107] In the above embodiments, the descriptions of each embodiment have different focuses. For parts not described in detail in a certain embodiment, please refer to the relevant descriptions in other embodiments.

[0108] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the methods, apparatuses, and devices described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0109] In the several embodiments provided in this application, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For instance, the division of modules is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple modules or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the displayed or discussed mutual couplings, direct couplings, or communication connections may be through some interfaces; indirect couplings or communication connections between devices or modules may be electrical, mechanical, or other forms.

[0110] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical modules; that is, they may be located in one place or distributed across multiple network modules. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0111] Furthermore, the functional modules in the various embodiments of this application can be integrated into one processing module, or each module can exist physically separately, or two or more modules can be integrated into one module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can be stored in a computer-readable storage medium.

[0112] In the above embodiments, implementation can be achieved, in whole or in part, through software, hardware, firmware, or any combination thereof. When implemented in software, it can be implemented, in whole or in part, as a computer program product.

[0113] The computer program product includes one or more computer instructions. When the computer program is loaded and executed on a computer, all or part of the processes or functions described in the embodiments of this application are generated. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, fiber optic, digital subscriber line) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that a computer can store or a data storage device such as a server or data center that integrates one or more available media. The available medium may be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., digital video optical disc), or a semiconductor medium (e.g., solid-state disk (SSD)).

[0114] Although embodiments of this application have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and variations can be made to these embodiments without departing from the principles of this application.

Claims

1. A soil environmental monitoring method applied to the ecological restoration of landfills, characterized in that, include: Obtain environmental monitoring data from each soil monitoring area of ​​the landfill; The environmental monitoring data includes at least the heavy metal content at different depths and the permeability coefficient under different confining pressures in each soil monitoring area; Based on various environmental monitoring data, the initial heavy metal solidification strength of humus at different depths was obtained. The initial heavy metal solidification strength is used at least to characterize the heavy metal solidification capacity of humus soil at the corresponding depth. Based on the initial heavy metal solidification intensity of humus at different depths, correction coefficients are obtained for each soil monitoring area; the correction coefficients are used at least to characterize the extent to which the heavy metal content of the corresponding soil monitoring area is affected by the heavy metal content of its adjacent soil monitoring areas. Based on the correction coefficient and the initial heavy metal solidification intensity, the corrected heavy metal solidification intensity of each soil monitoring area is obtained; Based on the modified heavy metal solidification intensity, the heavy metal content at different depths in each soil monitoring area was obtained.

2. The soil environmental monitoring method for ecological restoration of landfills according to claim 1, characterized in that, The acquisition of environmental monitoring data in each soil monitoring area of ​​the landfill includes: The landfill was divided into multiple first soil monitoring zones; Based on each first soil monitoring area, a second soil monitoring area is obtained; the second soil monitoring area is any area in each first soil monitoring area where no environmental monitoring data was obtained. Based on the second soil monitoring area, soil samples at different depths were obtained; Environmental monitoring data for the second soil monitoring area were obtained based on each soil sample.

3. The soil environmental monitoring method for ecological restoration of landfills according to claim 2, characterized in that, The method of obtaining the initial heavy metal solidification intensity of humus at different depths based on various environmental monitoring data includes: Based on the permeability coefficient under different confining pressures, the geological similarity under the same confining pressure is obtained; the geological similarity under the same confining pressure is at least used to characterize the similarity of the permeability coefficients of humus soils under the same confining pressure. Based on the permeability coefficient under different confining pressures, the support for heavy metal movement at the same depth is obtained; the support for heavy metal movement at the same depth is at least used to characterize the magnitude of the support capacity of a single-depth humus soil for heavy metal movement within its own confining pressure variation range. Based on the geological similarity under the same confining pressure and the support for heavy metal movement at the same depth, the comprehensive soil heavy metal diffusivity is obtained; the comprehensive soil heavy metal diffusivity is at least used to characterize the ability of a single depth of humus soil to diffuse heavy metals into adjacent depths of humus soil. Based on the comprehensive soil heavy metal diffusivity, the initial heavy metal solidification intensity is obtained; the initial heavy metal solidification intensity is used to characterize the ability of humus soil at different depths to fix heavy metals.

4. The soil environmental monitoring method for ecological restoration of landfills according to claim 3, characterized in that, The method of obtaining geological similarity under the same confining pressure based on the permeability coefficient under different confining pressures includes: A two-dimensional coordinate system is established based on the permeability coefficient under different confining pressures, and multiple coordinate points are obtained; the horizontal axis of the two-dimensional coordinate system is the confining pressure, and the vertical axis is the permeability coefficient. Fit the data to each coordinate point to obtain the permeation curve corresponding to each depth; Based on each permeation curve, a first difference value and a second difference value are obtained; the first difference value is the difference in permeation slope at different depths under the same confining pressure; the second difference value is the difference in permeability coefficient at different depths under the same confining pressure. Based on the first difference value and the second difference value, the geological similarity under the same confining pressure is obtained.

5. The soil environmental monitoring method for ecological restoration of landfills according to claim 4, characterized in that, The method of obtaining the support for heavy metal movement at the same depth based on the permeability coefficient under different confining pressures includes: For any depth within the same waste soil monitoring area, extract all permeability coefficients and corresponding confining pressure data within the permeability curve area corresponding to the empirical confining pressure range at that depth; the empirical confining pressure range is preset. Based on each permeability coefficient and each confining pressure, multiple sets of third and fourth difference values ​​are obtained; the third difference value is the difference between adjacent permeability coefficients among each permeability coefficient; the fourth difference value is the difference between adjacent confining pressures among each confining pressure. Based on the third and fourth difference values ​​of each group, the support for heavy metal movement at the same depth is obtained.

6. The soil environmental monitoring method for ecological restoration of landfills according to claim 5, characterized in that, The method of obtaining comprehensive soil heavy metal diffusivity based on the geological similarity under the same confining pressure and the support for heavy metal movement at the same depth includes: For any depth within the same waste soil monitoring area, analyze the empirical confining pressure range of that depth and adjacent depths to obtain the overlapping range; Based on the overlapping range, multiple first confining pressures are obtained; the first confining pressure is any confining pressure within the overlapping range corresponding to the depth of the permeation curve. Based on the geological similarity of each first confining pressure and the support for heavy metal movement at the same depth, the comprehensive soil heavy metal diffusivity is obtained.

7. The soil environmental monitoring method for ecological restoration of landfills according to claim 3, characterized in that, The method of obtaining the initial heavy metal solidification strength based on the comprehensive soil heavy metal diffusivity includes: Based on the heavy metal content at different depths within each soil monitoring area, the current state of heavy metal content in humus at each depth is obtained; the current state of heavy metal content is used at least to characterize the magnitude of heavy metal content in humus at the corresponding depth. The content of each heavy metal is arranged in order of depth from shallow to deep to obtain a depth sequence; The initial heavy metal solidification intensity is obtained based on the depth sequence and the comprehensive soil heavy metal diffusivity.

8. The soil environmental monitoring method for ecological restoration of landfills according to claim 7, characterized in that, The calculation formula for obtaining the initial heavy metal solidification intensity based on the depth sequence and the comprehensive soil heavy metal diffusivity is as follows: in, Indicates the initial heavy metal curing strength; This indicates the number of sequence values ​​in the depth sequence that precede the sequence value at the current depth. This indicates the overall soil heavy metal diffusivity at the current depth; This represents the overall soil heavy metal diffusivity at the m-th depth; This indicates the current state of heavy metal content at the m-th depth; This represents a normalization function used to normalize the sequence values ​​within the parentheses to the range [0, 1]. This represents an exponential function with the natural constant e as its base.

9. The soil environmental monitoring method for ecological restoration of landfills according to any one of claims 1 to 8, characterized in that, The correction coefficients for each soil monitoring area are obtained based on the initial heavy metal solidification intensity of humus at different depths, including: Based on each soil monitoring area, a third soil monitoring area is obtained; the third soil monitoring area is any soil monitoring area in each soil monitoring area for which no correction coefficient has been obtained. Based on the third soil monitoring area, a fourth soil monitoring area was obtained; the fourth soil monitoring area is adjacent to the third soil monitoring area, and no value was obtained for the degree of horizontal spread between the third soil monitoring area and the fourth soil monitoring area; Based on the third and fourth soil monitoring areas, multiple sampling distances and multiple fifth difference values ​​are obtained; the sampling distance is the distance between sampling points at the same depth in the third and fourth soil monitoring areas; the fifth difference value is the difference in the initial heavy metal solidification intensity at the same depth in the third and fourth soil monitoring areas. Based on each sampling distance and each fifth difference value, the horizontal spread impact value at different depths of the third soil monitoring area is obtained; the horizontal spread impact value is at least used to characterize the extent to which heavy metal substances in the third soil monitoring area can easily spread horizontally to the fourth soil monitoring area. Based on the horizontal spread impact value, the correction coefficient for the third soil monitoring area is obtained.

10. The soil environmental monitoring method for ecological restoration of landfills according to claim 9, characterized in that, The step of obtaining the correction coefficient for the third soil monitoring area based on the horizontal spread impact value includes: Based on the third soil monitoring area, multiple fifth soil monitoring areas are obtained; the fifth soil monitoring area is any soil monitoring area adjacent to the third soil monitoring area in each soil monitoring area, and the degree of horizontal spread influence between the third soil monitoring area and the fifth soil monitoring area has been obtained; Based on each fifth soil monitoring area, multiple horizontal spread impact values ​​were obtained; each horizontal spread impact value corresponds one-to-one with the fifth soil monitoring area. Based on the degree of spread of influence at each level, the correction coefficient for the third soil monitoring area is obtained.

Images