A Safety and Environmental Early Warning Method for Large-Scale Industrial Solid Waste Storage Sites Based on Digital Twins

Abstract

This disclosure provides a method for early warning of safety and environmental protection for large-scale industrial solid waste stockpiles based on digital twins. The method includes: acquiring the stockpile state variables within the target stockpile and constructing a comprehensive stockpile state set with a unified time reference based on these variables. Based on the comprehensive stockpile state set, a dual-state mapping structure is constructed in the digital twin space to obtain a dual-state digital twin stockpile model. The dual-state digital twin stockpile model is used to analyze the target stockpile area to identify its mechanical state characteristics and determine the critical instability pre-state based on these characteristics. The critical instability pre-state is used to determine whether the target stockpile area is in a critical state of mechanical degradation. When the target stockpile area is in a critical state of mechanical degradation, safety and environmental protection early warning information is output, enabling early warning of instability and pollution risks in high-moisture-content fine-particle solid waste stockpiles, thus improving the reliability and timeliness of safety and environmental protection risk assessment.

Description

Technical Field

[0001] This invention relates to the field of intelligent environmental protection technology, and in particular to a method for early warning of safety and environmental protection in large-scale industrial solid waste storage sites based on digital twins. Background Technology

[0002] In high-moisture-content fine-particle solid waste (such as fly ash and red mud) stockpiles, when the local medium is in the critical moisture content range transitioning from subsaturation to saturation for a long period, the stockpiles typically exhibit a typical phenomenon of strength degradation and geometric response decoupling. Due to the continuous decay of matrix suction and the reduction of effective contact stress between particles, the shear strength and bearing path inside the stockpiles have undergone substantial deterioration. However, this deterioration process is not necessarily accompanied by observable changes in surface morphology or significant displacement response. Thus, digital twin state mapping, which is constrained by geometric consistency and displacement continuity, can only maintain synchronization at the morphological level, but it is difficult to identify at the mechanical state level that the stockpiles have entered the pre-critical instability state.

[0003] For example, under the combined conditions of continuous light rain and leachate recharge at a red mud dump, the surface elevation and slope profile remain relatively stable over several weeks. However, the fine-grained internal zone, due to the gradual accumulation of pore water pressure and weakened structural cementation, approaches shear failure conditions. A short period of heavy rainfall or minor external disturbance could easily trigger localized slippage or fluidization failure. However, within the existing digital twin framework, the lack of engineering representation capabilities to separate the apparent stability from the internal instability makes it difficult for current digital twin models to identify dump risks in a timely and reliable manner.

[0004] Therefore, for high-moisture-content fine-particle solid waste stockpiles in the critical moisture range transitioning from subsaturation to saturation, where the internal shear strength and load-bearing structure of the stockpile have undergone substantial degradation, but the macroscopic geometry and surface displacement remain stable, existing digital twin methods can only achieve morphological synchronization and cannot identify the pre-critical instability state. This invention proposes a safety and environmental protection early warning method for large-scale industrial solid waste stockpiles based on digital twins. By explicitly constructing a separable expression and discrimination mechanism for the geometrically stable state and the internal mechanical degradation state within the digital twin framework, it achieves early identification and safety and environmental protection warnings for the transition of fine-particle solid waste stockpiles from apparent stability to critical instability. Summary of the Invention

[0005] Therefore, it is necessary to propose a safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins to address the aforementioned technical issues.

[0006] The present invention adopts the following technical solution.

[0007] The first aspect of this invention discloses a safety and environmental early warning method for large-scale industrial solid waste storage sites based on digital twins, the method comprising:

[0008] Obtain the state variables of the stack within the target stack, and construct a unified time-based integrated state set of the stack based on the stack state variables;

[0009] Based on the aforementioned integrated state set of the heap, a dual-state mapping structure of the heap is constructed in the digital twin space to obtain a dual-state digital twin heap model.

[0010] The dual-state digital twin reactor model is invoked to analyze the target reactor region in order to identify the reactor mechanical state characteristics, and the reactor critical instability pre-state is determined based on the reactor mechanical state characteristics;

[0011] Based on the pre-critical instability state of the reactor body, it is determined whether the target reactor body region is in a critical state of mechanical degradation, and when the target reactor body region is in a critical state of mechanical degradation, safety and environmental protection early warning information is output.

[0012] Furthermore, the step of acquiring the stack state quantities within the target stack yard and constructing a unified time-based integrated stack state set based on the stack state quantities includes:

[0013] Multiple monitoring points are selected within the target stockpile and numbered to obtain the water content, pore pressure, and apparent morphology of the stockpile at each monitoring point according to a set sampling period.

[0014] Based on the water content, pore pressure, and apparent morphology of the pile body, and combined with the preset critical reference values, the comprehensive state quantity of each monitoring point at each time is calculated.

[0015] The stack state quantities include the stack water content state quantity, pore pressure state quantity, and apparent morphology state quantity.

[0016] Furthermore, the step of acquiring the stack state quantities within the target stack yard and constructing a unified time-based integrated stack state set based on the stack state quantities also includes:

[0017] The comprehensive state quantities corresponding to each monitoring point at each time are summarized into the comprehensive state quantity of the reactor body, and a subcritical criterion threshold is set to determine the subcritical coverage ratio when the comprehensive state quantity is not lower than the subcritical criterion threshold.

[0018] Based on the overall state quantity of the reactor body and the subcritical coverage ratio, a reactor body strength threshold and a coverage ratio threshold are set, and a transition interval state identifier is determined when the overall state quantity of the reactor body is not lower than the reactor body strength threshold and the subcritical coverage ratio is not lower than the coverage ratio threshold.

[0019] The integrated state set of the reactor body is composed of the integrated state quantity of the reactor body, the subcritical coverage ratio, and the transition interval state identifier.

[0020] Furthermore, based on the comprehensive state set of the heap, a two-state mapping structure for the heap is constructed in the digital twin space to obtain a two-state digital twin heap model, including:

[0021] A geometric mapping channel is established in the digital twin space, and a geometric mapping reference and a geometric mapping effective area are established for the integrated state set of the heap through the geometric mapping channel;

[0022] Based on the geometric mapping datum and the effective area of ​​geometric mapping, the surface of the digital twin stack is divided into multiple partition units, and the elevation of the digital twin surface of each partition unit is determined.

[0023] Each monitoring point in the effective area of ​​the geometric mapping is projected to the corresponding partition unit, and the equivalent shape update amount of the partition unit is calculated. The elevation of the digital twin surface is updated by the equivalent shape update amount to obtain the geometric state mapping amount and the geometric twin surface state.

[0024] Furthermore, the step of constructing a dual-state mapping structure for the heap in the digital twin space based on the heap's integrated state set to obtain a dual-state digital twin heap model further includes:

[0025] The target pile body region is divided into multiple material partitions, and the reference shear strength parameters of each material partition are obtained. At the same time, the mechanical degradation coefficient of each material partition is defined according to the comprehensive state set of the pile body.

[0026] The twin mechanical parameters are obtained by equivalently reducing the reference shear strength parameters using the mechanical degradation coefficient.

[0027] Obtain the equivalent values ​​of normal stress on the slope of the pile body, void pressure, and weighted average value of friction force, and calculate the equivalent shear strength and equivalent driving shear stress based on the equivalent values ​​of normal stress on the slope of the pile body, void pressure, and weighted average value of friction force.

[0028] The mechanical state mapping quantity is calculated based on the equivalent shear strength and equivalent driving shear stress, and the twin mechanical parameters, equivalent shear strength and equivalent driving shear stress are integrated into a mechanical twin bearing state.

[0029] The dual-state digital twin stack model is composed of a stack dual-state mapping structure, which is composed of the time-synchronized mechanical state mapping quantity and the mechanical twin bearing state.

[0030] Furthermore, the step of calling the dual-state digital twin reactor model to analyze the target reactor region to identify the reactor mechanical state characteristics, and determining the reactor's critical instability pre-state based on the reactor mechanical state characteristics, includes:

[0031] Based on the aforementioned dual-state mapping structure and integrated state set of the stack, a transition interval identifier for subcritical effective discrimination is constructed, and the geometric stability constraint is calculated according to the set geometric stability threshold.

[0032] Obtain the start time of the transition interval identifier, define the shear capacity attenuation characteristic based on the geometric stability constraint and the equivalent shear strength at the start time of the interval, and calculate the shear capacity attenuation rate based on the shear capacity attenuation characteristic and the set sampling period.

[0033] Furthermore, the step of calling the dual-state digital twin reactor model to analyze the target reactor region to identify reactor mechanical state characteristics, and determining the reactor's critical instability pre-state based on these characteristics, also includes:

[0034] The difference between the equivalent shear strength and the equivalent driving shear stress is normalized to obtain the critical equilibrium characteristic quantity, and the critical approximation index is calculated based on the critical equilibrium characteristic quantity and the rate of change of shear capacity decay.

[0035] Based on the critical equilibrium characteristic quantity and the critical approximation index, the partition strength reduction ratio corresponding to each material partition is defined, and the partition strength reduction ratio is weighted and summarized according to the geometric stability constraint quantity to obtain the bearing path weakening characteristic quantity.

[0036] When the bearing path weakening characteristic, critical approximation index, and shear capacity attenuation characteristic are all not lower than their respective set thresholds, the pre-critical instability state of the reactor body is output.

[0037] Furthermore, the step of determining whether the target reactor region is in a critical state of mechanical degradation based on the pre-critical instability state of the reactor body, and outputting safety and environmental protection early warning information when the target reactor region is in a critical state of mechanical degradation, includes:

[0038] Based on the pre-critical instability state of the reactor and the time-synchronized dual-state mapping structure of the reactor, a warning trigger identifier is defined, and the warning information is classified and processed in response to the warning trigger identifier to obtain multiple warning levels and the warning confidence strength corresponding to each level.

[0039] Based on the warning level, warning confidence strength, and warning trigger identifier, safety warning types and environmental warning types are divided, and priorities are set for the safety warning types and environmental warning types so that the safety and environmental warning information is output according to the priorities.

[0040] The second aspect of this invention discloses a safety and environmental protection early warning device for large-scale industrial solid waste storage sites based on digital twins, used to implement the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins as described in any one of the first aspects, the device comprising:

[0041] The stack state construction module is used to obtain the stack state quantities in the target stack yard and construct a unified time-based integrated stack state set based on the stack state quantities.

[0042] The model building module is used to construct a dual-state mapping structure of the heap in the digital twin space based on the heap's comprehensive state set, so as to obtain a dual-state digital twin heap model.

[0043] The mechanical state analysis module is used to call the dual-state digital twin reactor model to analyze the target reactor region, identify the mechanical state characteristics of the reactor, and determine the critical instability pre-state of the reactor based on the mechanical state characteristics of the reactor.

[0044] The early warning information output module is used to determine whether the target reactor area is in a critical state of mechanical degradation based on the pre-critical instability state of the reactor body, and to output safety and environmental protection early warning information when the target reactor area is in a critical state of mechanical degradation.

[0045] A third aspect of the present invention discloses a terminal, including a processor and a storage medium;

[0046] The storage medium is used to store instructions;

[0047] The processor is configured to operate according to the instructions to perform the steps of the method described in the first aspect.

[0048] A fourth aspect of the present invention discloses a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method described in the first aspect.

[0049] The present invention has the following advantages:

[0050] (1) This invention collects state quantities reflecting the water content, pore pressure, and apparent morphology of a stockpile of fine-particle solid waste during its operation in the stockpile. A comprehensive set of state quantities is constructed under a unified time reference to describe the current water content and stress state range of the stockpile. Subsequently, a dual-state mapping structure is constructed in a digital twin space to describe the geometric state mapping of the stockpile's external stability and the mechanical state mapping of its internal bearing capacity and shear resistance. These two types of states are updated synchronously in the time dimension but are independent in the discrimination logic. This ensures that subsequent judgments do not use inconsistent state outputs and maintains the independence of the discrimination logic, thus achieving the engineering expression capability of separating the externally stable and internally unstable states.

[0051] (2) Based on the dual-state digital twin model of the stack, this invention performs continuous analysis on the mechanical state mapping of the stack region in the transition zone from subsaturation to saturation. This analysis identifies the characteristics of weakened bearing path, reduced shear capacity, or critical equilibrium reflected in the mechanical state mapping under the condition that the geometric state remains stable. Based on this, it determines whether the stack has entered the pre-critical instability state. Finally, the pre-critical instability state judgment result is used as a direct triggering basis. When it is determined that the stack is in a critical state of geometric stability but significant mechanical degradation, the corresponding safety and environmental protection early warning information is output to prompt the stack site manager to take targeted measures before any abnormalities occur in appearance. This achieves early warning of instability and pollution risks of high-moisture fine particulate solid waste stack sites, improving the reliability and timeliness of safety and environmental protection risk judgment. Attached Figure Description

[0052] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0053] Figure 1 This is a flowchart illustrating the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by this invention.

[0054] Figure 2 This is a schematic diagram of the structure of the safety and environmental protection early warning device for large-scale industrial solid waste storage sites based on digital twins provided by the present invention. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0056] like Figure 1 As shown in one embodiment, a method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins includes the following steps:

[0057] Step S110: Obtain the stack state variables within the target stack yard, and construct a unified time-based integrated stack state set based on the stack state variables.

[0058] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention includes the following steps in step S110:

[0059] Step S111: Select multiple monitoring points within the target stockpile and number them to obtain the water content, pore pressure, and apparent morphology of each monitoring point according to the set sampling period.

[0060] Step S112: Based on the water content state quantity, pore pressure state quantity, and apparent morphology state quantity of the pile body, and combined with the preset critical reference value, calculate the comprehensive state quantity corresponding to each monitoring point at each time.

[0061] Among them, the state parameters of the pile body include the water content state parameter, the pore pressure state parameter, and the apparent morphology state parameter. The apparent morphology state parameter is a state parameter used to describe the changes in the appearance of the pile body, specifically the change in the local elevation of the slope at the monitoring point over time.

[0062] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention further includes the following steps in step S110:

[0063] Step S113: The comprehensive state quantities corresponding to each monitoring point at each time are summarized into the comprehensive state quantity of the reactor body, and a subcritical criterion threshold is set to determine the subcritical coverage ratio when the comprehensive state quantity is not lower than the subcritical criterion threshold.

[0064] Step S114: Based on the overall state quantity of the reactor body and the subcritical coverage ratio, set the reactor body strength threshold and the coverage ratio threshold, and determine the transition interval state identifier when the overall state quantity of the reactor body is not lower than the reactor body strength threshold and the subcritical coverage ratio is not lower than the coverage ratio threshold.

[0065] The integrated state set of the reactor body consists of the integrated state quantity of the reactor body, the subcritical coverage ratio, and the state identifier of the transition zone.

[0066] In a specific embodiment, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention includes steps 1 to 4:

[0067] Step 1: Constructing the unified state variables for the subsaturated and saturated transition states of the heap.

[0068] During the operation of the stockpile, for high-moisture-content fine-particle solid waste stockpiles, state variables reflecting the water content, pore pressure, and apparent morphology of the stockpile are collected. A comprehensive set of stockpile state variables is then constructed under a unified time reference to describe the current water content and stress state range of the stockpile. This includes the following sub-steps:

[0069] Sub-step 1.1: Construct the original sampling records under a unified time reference.

[0070] Specifically, firstly, monitoring points are selected within the stockpile area. The engineering values ​​for these monitoring points are integers ranging from 20 to 500, covering the top of the stockpile, the slope surface, the toe of the slope, and the leachate collection area. Then, a uniform sampling period is set, ranging from 1 to 60 minutes. At each sampling time, water content, pore pressure, and apparent morphology parameters are collected for each monitoring point and written to the original sampling record with the same timestamp. Furthermore, to avoid directly mixing different units of measurement, sampling and archiving are only performed under the same time reference; no fusion calculations are performed.

[0071] Sub-step 1.2: Calculate the comprehensive state variables at the monitoring point level.

[0072] Specifically, for each monitoring point, the corresponding point-level comprehensive state quantity is calculated at each time point to compress the water content, pore pressure, and apparent morphology into a single engineering quantity that can be used to determine the subcritical range. In calculating the point-level comprehensive state quantity, the water content is first defined as volumetric water content, ranging from 0.05 to 0.70; the pore pressure ranges from -50 to 300 kPa; and the apparent morphology is defined as the local elevation change of the slope, ranging from -0.5 to 0.5 m. Subsequently, three benchmark constants were defined as critical reference values ​​for these three state variables. The critical reference value for volumetric water content ranged from 0.20 to 0.45, the critical reference value for pore pressure ranged from 0 to 50 kPa, and the critical reference value for local elevation change on the slope ranged from 0.01 to 0.05 m. Finally, based on the weighting coefficients assigned to these three state variables (ranging from 0.2 to 2.0, and the sum of the three weighting coefficients being 1) and their corresponding critical reference values, the comprehensive state variable at the monitoring point level was calculated. The numerical value of this comprehensive state variable at the monitoring point level is equal to the weighted sum of the products of the ratios of the three state variables to their respective critical reference values ​​and the corresponding weighting coefficients. Since elevation change can be negative, its absolute value must be taken.

[0073] Sub-step 1.3: Construct the heap's comprehensive state variable set.

[0074] Specifically, to describe the overall water content and stress state range of the reactor body at a certain moment, it is necessary to aggregate the comprehensive state quantities at each monitoring point into a reactor-level comprehensive state quantity, and simultaneously output the subcritical coverage ratio to characterize how many areas of the reactor body are in the critical transition phase. First, a representative area is defined for each monitoring point, ranging from 50 to 5000 square meters. The reactor-level comprehensive state quantity is calculated using an area-weighted average. Next, a subcritical criterion threshold is defined, ranging from 1.0 to 2.5, to mark dangerous areas that are close to saturation but have not undergone significant deformation. The ratio of the weighted sum of the areas of these dangerous areas to the total area is the subcritical coverage ratio. Areas with a monitoring point-level comprehensive state quantity greater than or equal to the subcritical criterion threshold are considered dangerous areas. This dangerous area can be calculated by assigning a value of 1 to the monitoring point-level comprehensive state quantity when it is greater than or equal to the subcritical criterion threshold using an indicator function. Finally, the reactor-level comprehensive state quantity set, consisting of the reactor-level comprehensive state quantity and the subcritical coverage ratio, is output.

[0075] Sub-step 1.4: Determine the state indicators of the subsaturated and saturated transition intervals.

[0076] Specifically, to determine whether the reactor is currently in the critical water-bearing range transitioning from subsaturation to saturation, a state indicator is constructed at each time point. This indicator must reflect both the overall degradation intensity and the proportion of the critical region to avoid misjudgments such as localized high danger but low overall average, or high overall average but small danger zone. First, the reactor-level strength threshold is set to a range of 1.0-2.0, and the coverage ratio threshold is set to a range of 0.10-0.40. If the reactor-level comprehensive state quantity is greater than or equal to the reactor-level strength threshold and the subcritical coverage ratio is greater than or equal to the coverage ratio threshold, it is identified as a state indicator for the subsaturation to saturation transition range and assigned a value of 1; otherwise, it is assigned a value of 0. Finally, instead of relying on a single indicator to trigger the subcritical transition state, a comprehensive reactor state set consisting of the reactor-level comprehensive state quantity, the subcritical coverage ratio, and the transition range state indicator is output.

[0077] Step S120: Based on the integrated state set of the heap, construct a dual-state mapping structure of the heap in the digital twin space to obtain a dual-state digital twin heap model.

[0078] Understandably, the dual-state mapping structure of the reactor core is achieved by establishing geometric mapping channels and mechanical mapping channels separately within the digital twin space. The geometric mapping channel, constrained by the geometric mapping reference and the effective geometric mapping area, projects the apparent morphological state quantities of the monitoring points onto the partitioned units, forming geometric state mapping quantities and geometric twin surface states. The mechanical mapping channel, based on material partitions and reference shear strength parameters, incorporates a mechanical degradation coefficient into the reactor core's comprehensive state set to generate twin mechanical parameters, equivalent shear strength, and equivalent driving shear stress, and outputs the mechanical state mapping quantities. The two types of mappings are updated synchronously in the time dimension and remain independent in their discrimination logic.

[0079] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention includes the following steps in step S120:

[0080] Step S121: Establish a geometric mapping channel in the digital twin space, and establish a geometric mapping benchmark and a geometric mapping effective area for the integrated state set of the heap through the geometric mapping channel.

[0081] Step S122: Based on the geometric mapping datum and the effective area of ​​geometric mapping, the surface of the digital twin stack is divided into multiple partition units, and the elevation of the digital twin surface of each partition unit is determined.

[0082] In this context, the digital twin surface elevation refers to the elevation representation of the stack surface corresponding to each partition unit in the digital twin space, used to reflect the geometric shape of the stack. It is obtained by updating the initial twin surface elevation in conjunction with the equivalent shape update amount every moment.

[0083] Step S123: Project each monitoring point in the effective area of ​​geometric mapping to the corresponding partition unit, and calculate the equivalent morphological update amount of the partition unit, so as to update the elevation of the digital twin surface through the equivalent morphological update amount, and obtain the geometric state mapping amount and the geometric twin surface state.

[0084] Among them, the geometric twin surface state is the overall geometric state expression composed of the updated results of the digital twin surface elevation on each partition unit. It is used to uniformly reflect the stability of the pile appearance and is derived from the partition convergence results of the morphological changes of monitoring points within the effective area of ​​geometric mapping.

[0085] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention further includes the following steps in step S120:

[0086] Step S124: Divide the target stack area into multiple material partitions, obtain the baseline shear strength parameters of each material partition, and define the mechanical degradation coefficient of each material partition according to the stack's comprehensive state set.

[0087] Step S125: The reference shear strength parameters are equivalently reduced by the mechanical degradation coefficient to obtain the twin mechanical parameters.

[0088] Step S126: Obtain the equivalent values ​​of normal stress on the slope of the accumulator, void pressure, and weighted average value of frictional force, and calculate the equivalent shear strength and equivalent driving shear stress based on the equivalent values ​​of normal stress on the slope of the accumulator, void pressure, and weighted average value of frictional force.

[0089] Step S127: Calculate the mechanical state mapping quantity based on the equivalent shear strength and equivalent driving shear stress, and integrate the twin mechanical parameters, equivalent shear strength and equivalent driving shear stress into a mechanical twin bearing state.

[0090] The dual-state digital twin stack model consists of a stack dual-state mapping structure, which is composed of time-synchronized mechanical state mapping quantities and mechanical twin bearing states.

[0091] In a specific embodiment, the present invention provides a method for early warning of safety and environmental protection of large-scale industrial solid waste dumps based on digital twins. Step 2 involves constructing a dual-state mapping of the digital twin dump body. Based on the comprehensive state set of the dump body output in step 1, a dual-state mapping structure of the dump body is constructed in the digital twin space. This structure includes a geometric state mapping to describe the stability of the dump body's appearance and a mechanical state mapping to describe the load-bearing and shear resistance of the dump body's interior. The two types of states are updated synchronously in the time dimension but are independent in their discrimination logic. This includes the following sub-steps:

[0092] Sub-step 2.1: Construct the geometric mapping reference quantity and the set of geometric mapping effective regions of the twin space.

[0093] Specifically, a geometric mapping channel for appearance stability is first established in the digital twin space. This geometric mapping is not equivalent to stuffing all sensor data into the 3D model, but rather establishing a reproducible mapping benchmark only for visible geometric quantities. First, the slope elevation change of the pile body at each time point is determined, and then the geometric stability contribution of each monitoring point is defined. This contribution is used to compress the elevation change into stability intensity. Next, the geometric reference threshold is set to a range of 0.01-0.05 m to filter out minor undulations and construction disturbances. The geometric stability contribution is equal to the ratio of the absolute value of the elevation change to its critical reference value plus the reciprocal of the result (1). This represents the closer the contribution is to 1 for more stable geometric shapes, and the closer the contribution is to 0 for greater geometric fluctuations.

[0094] Subsequently, the geometric stability contributions of each monitoring point were weighted according to the area represented by the monitoring point and aggregated into a stack-level geometric mapping benchmark. This stack-level geometric mapping benchmark ranges from 0 to 1, with values ​​closer to 1 indicating greater geometric stability. Finally, to form a set of effective geometric mapping zones, monitoring points were screened. Points whose absolute elevation changes were less than or equal to their critical reference values ​​were considered to be in the geometrically stable zone and included in the set of effective geometric mapping zones. This was to constrain subsequent geometric twin surface updates to maintain morphological synchronization only within the stable zone, thus avoiding misinterpreting local construction disturbances as overall deformation.

[0095] Sub-step 2.2: Calculate the geometric state mapping quantity and the geometric twin surface state.

[0096] Specifically, a geometric state mapping is constructed in the twin space, not aiming for full-field fine 3D reconstruction, but rather using geometric states that can be used for stability assessment as output. First, the surface of the twin stack is defined as being represented by partitioned cells, with the number of partitions ranging from 50 to 500. Simultaneously, the twin surface elevation of each partitioned cell is determined. Each monitoring point in the set of effective geometric mapping regions is projected onto its corresponding partitioned cell. Within the geometrically stable region, area weighting is used to map discrete elevation changes to equivalent morphological updates of the partitioned cells, yielding the equivalent morphological update amount of the partitioned cell. This is the weighted sum of the elevation changes of the partitioned cells belonging to the set of effective geometric mapping regions and their representative areas, used as the numerator, and the ratio of this sum to the total area of ​​the partitioned cells belonging to the set of effective geometric mapping regions.

[0097] Next, the surface elevation at the current moment minus the previous sampling period is summed with the equivalent morphological update amount, which is the updated twin surface elevation. Then, a geometric state mapping quantity is further defined to output a unified representation of appearance stability in twin space. Its value ranges from 0 to 1, considering both the overall geometric stability benchmark and the total fluctuation intensity of twin surface partition updates, and is used to output the geometric mapping as a discriminable state quantity. Finally, the updated twin surface elevation is used as the twin surface state as the geometric state mapping output.

[0098] Sub-step 2.3: Construct the mechanical state mapping quantity and the mechanical twin bearing state.

[0099] Specifically, a mechanical state mapping of internal bearing capacity and shear resistance is constructed in the twin space. This mechanical mapping does not rely on the stability conclusions of the geometric mapping, but rather on the expression of water content and pore pressure degradation in the comprehensive state quantity set output from step 1, and allows for the output of high-risk mechanical states even when geometrically stable. First, material zones of the reactor body are defined, with a suggested range of 3-20 zones to distinguish between fly ash layers, red mud layers, and overburden layers. Two benchmark shear strength parameters are determined for each material zone, one ranging from 5-50 kPa and the other from 10-35 kPa. Next, a mechanical degradation coefficient for each zone is defined, driven by the comprehensive state quantities of the reactor body from step 1. When the transition interval state indicates a subcritical transition state, the mechanical degradation coefficient increases with the increase of the comprehensive state quantities of the reactor body and the proportion of subcritical overburden. Therefore, the mechanical degradation coefficient is numerically equal to the product of the reactor-level integrated state quantity minus 1 and the subcritical coverage ratio, plus a set constant, which ranges from 0.5 to 2.0. The mechanical degradation coefficient is a coupled amplification of degradation intensity and subcritical coverage ratio, used to explicitly inject degradation effects in the critical range into the mechanical mapping.

[0100] Subsequently, the zoned strength parameters are equivalently reduced based on the mechanical degradation coefficient to obtain twin mechanical parameters. These twin mechanical parameters are also divided into two, and are numerically equal to the ratio of the numerator of the reference shear strength parameter to the product of the reduction sensitivity coefficient and the mechanical degradation coefficient, plus the value of 1 as the denominator. This ratio is used to characterize that as the degradation coefficient increases, the shear parameter decreases in the form of the denominator, so as to avoid negative values ​​and maintain engineering feasibility.

[0101] In this embodiment, it is necessary to further define the mechanical state mapping quantity to output a scalar representation of the overall load-bearing safety margin of the reactor body. First, the representative normal stress of the reactor body slope (range 10-200 kPa), the equivalent value of void pressure (range 0-300 kPa), and the weighted average value of friction force by region are determined. Then, the equivalent shear strength and the equivalent driving shear stress are calculated. The equivalent shear strength is numerically equal to the product of the difference between the representative normal stress of the reactor body slope and the equivalent value of void pressure, multiplied by the weighted average value of friction force by region, plus the result of the weighted average value of twin mechanical parameters of each region by volume fraction (range 1-50 kPa). The equivalent driving shear stress is numerically equal to the product of the representative normal stress of the reactor body slope and 0.3 (where 0.3 is the empirical proportionality coefficient of slope shear stress to normal stress).

[0102] Finally, using the equivalent shear strength as the numerator and the equivalent driving shear stress plus a value of 1 (to prevent the denominator from being 0) as the denominator, the calculated ratio is the mechanical state mapping quantity. The smaller this mechanical state mapping quantity, the closer it is to instability. Simultaneously, the mechanical twin bearing state, composed of the two twin mechanical parameters, the equivalent shear strength, and the equivalent driving shear stress, is output.

[0103] Sub-step 2.4 generates a dual-state synchronous update package.

[0104] Specifically, to ensure that the two types of states are updated synchronously in the time dimension but independent in the discrimination logic, only time synchronization and packaging are performed, without fusion judgment. First, a synchronization consistency check is set to verify that the two types of states are updated at the same time, and the geometric mapping update timestamp and the mechanical mapping update timestamp are determined. The timestamps are allowed to have a small delay (range 0-5 minutes). Then, the synchronization consistency check is numerically equal to the ratio of the absolute value of the difference between the geometric mapping update timestamp and the mechanical mapping update timestamp to the small delay, plus the reciprocal of the result of 1. When the synchronization consistency check is close to 1, it indicates good synchronization. When the synchronization consistency check is significantly less than 1, it indicates that waiting or resampling is required to ensure that the two states at the same time are comparable. Finally, the output is a dual-state synchronous update package containing the geometric state mapping, the geometric twin surface state, the mechanical state mapping, and the mechanical twin bearing state.

[0105] Step S130: Call the dual-state digital twin reactor model to analyze the target reactor region in order to identify the reactor mechanical state characteristics and determine the critical instability pre-state of the reactor based on the reactor mechanical state characteristics.

[0106] Understandably, the mechanical state characteristics of a reactor core refer to the set of states used to characterize the internal bearing capacity and shear safety level of the reactor core in a dual-state digital twin reactor core model. Specifically, it consists of mechanical state mapping quantities, twin mechanical parameters, equivalent shear strength, equivalent driving shear stress, and derived shear capacity attenuation characteristics, safety margin, and bearing path weakening characteristics. The pre-critical instability state of the reactor core, on the other hand, is the judgment result formed when the above mechanical state characteristics are simultaneously satisfied in terms of temporal continuity and threshold conditions, under the constraint of maintaining geometric stability. The difference between the two is that the former is a continuously updated state description, while the latter is a critical identification of the state evolution result.

[0107] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention includes the following steps in step S130:

[0108] Step S131: Based on the dual-state mapping structure of the stack and the integrated state set of the stack, construct the transition interval identifier for subcritical effective discrimination, and calculate the geometric stability constraint quantity according to the set geometric stability threshold.

[0109] The transition interval identifier is used to mark the state of the reactor body in the critical interval of evolution from subsaturation to saturation. It is determined by combining the comprehensive state variables of the reactor body with the subcritical coverage ratio. The transition interval identifier for effective subcritical discrimination is based on the establishment of the state identifier, and further superimposed with the constraint that the geometric state mapping quantity meets the geometric stability threshold, which is used to limit the effective triggering interval of mechanical degradation analysis.

[0110] Step S132: Obtain the start time of the transition interval identifier, define the shear capacity attenuation characteristic quantity based on the geometric stability constraint quantity and the equivalent shear strength at the start time of the interval, and calculate the shear capacity attenuation change rate based on the shear capacity attenuation characteristic quantity and the set sampling period.

[0111] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention further includes the following steps in step S130:

[0112] Step S133: Normalize the difference between the equivalent shear strength and the equivalent driving shear stress to obtain the critical equilibrium characteristic quantity, and calculate the critical approximation index based on the critical equilibrium characteristic quantity and the rate of change of shear capacity decay.

[0113] Step S134: Based on the critical equilibrium characteristic quantity and the critical approximation index, define the partition strength reduction ratio corresponding to each material partition, and weight and summarize the partition strength reduction ratio according to the known weights based on the geometric stability constraint quantity to obtain the bearing path weakening characteristic quantity.

[0114] Step S135: When the load-bearing path weakening characteristic, critical approximation index and shear capacity attenuation characteristic are all not lower than their respective set thresholds, output the pre-critical instability state of the reactor body.

[0115] In a specific embodiment, the present invention provides a method for early warning of safety and environmental protection of large-scale industrial solid waste dumps based on digital twins. Step 3 involves identifying the mechanical degradation state within the subcritical water-bearing range. Based on the dual-state digital twin dump model constructed in step 2, continuous analysis is performed on the mechanical state mapping of the dump region in the transition range from subsaturation to saturation. This analysis identifies the weakening of the bearing path, the attenuation of shear capacity, or the critical equilibrium characteristics reflected in the mechanical state mapping under the condition that the geometric state remains stable. Based on this, it determines whether the dump has entered the pre-critical instability state, including the following sub-steps:

[0116] Sub-step 3.1: Construct the subcritical effective discrimination interval identifier and geometric stability constraint.

[0117] Specifically, mechanical degradation discrimination is initiated only in the critical range between subsaturation and saturation transition, under conditions of stable appearance. This avoids mixing significant deformation periods and construction disturbance periods into the discrimination logic between appearance stability and internal degradation. First, a subcritical effective discrimination interval identifier is constructed. When the reactor body is in the subcritical transition range, geometric stability is constrained by a geometric stability threshold (range 0.80-0.95). Therefore, when the reactor body is in the subcritical transition range (the transition range state identifier is assigned a value of 1) and the geometric state mapping is greater than or equal to the geometric stability threshold, the subcritical effective discrimination interval identifier is determined and assigned a value of 1; otherwise, the identifier is assigned a value of 0. This discrimination process first locks the subcritical water-bearing range and then the appearance stability state, ensuring that the subsequently derived degradation characteristics correspond to the key defect scenarios you are addressing.

[0118] Next, the constraint strength of appearance stability on mechanical judgment is quantified, and a geometric stability constraint is constructed to weight subsequent index calculations. This makes the system emphasize the consideration of intrinsic degradation when the appearance is more stable. Therefore, the geometric stability constraint is numerically equal to the 1st to 3rd power of the geometric state mapping. The closer the geometric state mapping is to the value of 1, the closer the geometric stability constraint is to the value of 1.

[0119] Sub-step 3.2: Calculate the shear capacity attenuation characteristic quantity and the shear capacity attenuation rate.

[0120] Specifically, the core of shear capacity attenuation is the decrease in the upper limit of the shear capacity that can be carried over time. In step 2, the equivalent shear capacity has already been characterized by the equivalent shear strength. Therefore, based on the equivalent shear strength, a directly identifiable attenuation characteristic quantity and rate of change are constructed. First, let the reference time be the first time entering the subcritical discrimination interval, that is, when the subcritical effective discrimination interval identifier at the current time is assigned a value of 1 and the subcritical effective discrimination interval identifier at the time before the sampling period is assigned a value of 0. The equivalent shear strength at the reference time is limited to a reasonable range of 5-300 kPa. The shear capacity attenuation characteristic quantity is defined to be calculated only when the subcritical effective discrimination interval identifier at the current time is assigned a value of 1, and the expression is: ; In the formula, for The characteristic quantity of shear capacity decay at time; As the reference time; The equivalent shear strength at the reference time; for The equivalent shear strength at a given moment; for Geometric stability constraints at time t.

[0121] As the equivalent shear strength decreases, the shear capacity decay characteristic increases, and the geometric stability constraint is used to emphasize the sensitivity to decay under geometrically stable conditions. Subsequently, a decay rate of change is further defined to identify whether decay is accelerating. Numerically, it is equal to the difference between the current shear capacity decay characteristic and the shear capacity decay characteristic at the time corresponding to the time after subtracting the sampling period. This difference is used as the numerator, and the ratio of the numerator to the sampling period plus 1 as the denominator. This ratio characterizes the decay increment per unit time; the rate of change increases as decay continues to accelerate.

[0122] Sub-step 3.3: Calculate the critical equilibrium characteristic and the critical approximation index.

[0123] Specifically, the engineering meaning of the critical equilibrium characteristic is that the safety margin between the shear capacity and the actual driving shear stress approaches zero. Since equivalent shear strength and equivalent driving shear stress already exist in step 2, this sub-step constructs two complementary quantities: one is the safety margin, and the other is the speed at which the safety margin approaches zero. First, the safety margin is defined as the normalized result of the difference between the shear capacity and the driving shear stress, expressed as: ; In the formula, for Safety margin at all times; for Geometric stability constraints at time t; for The equivalent shear strength at a given moment; for The equivalent driving shear stress at any given moment; It is a constant with a value range of 1-20 kPa, used to suppress false amplification under low stress conditions.

[0124] When the equivalent shear strength approaches the equivalent driving shear stress, the safety margin approaches zero. When the equivalent shear strength is less than the equivalent driving shear stress, the safety margin is negative, indicating the onset of instability. To identify whether the system is approaching the critical point, a critical approach index is defined. This index considers both the decrease in safety margin and the accelerated shear attenuation, and its expression is: ; In the formula, for The critical approach exponent at any given moment; , Both are weighting coefficients, with values ​​ranging from 0.2 to 0.8, and their sum equals the value 1; This serves as a margin reference threshold, with a value ranging from 0.05 to 0.30. for The rate of change of shear capacity at time t.

[0125] In this embodiment, when the safety margin is lower than the margin reference threshold and the rate of change of shear capacity decay is positive and increasing, the critical approximation index increases significantly, indicating a critical approximation trend. The safety margin and critical approximation index are then output for the final critical instability pre-state determination.

[0126] Sub-step 3.4: Output the results of the weakening characteristic quantity of the bearing path and the pre-critical instability determination.

[0127] Specifically, the engineering implications of weakened load-bearing paths are not simply increased stress, but rather a shift in the load-bearing mechanism upon which the reactor body maintains equilibrium from friction-cementation control to near-slip control. This manifests as a synchronous decrease in strength parameters over time, leading to a continuous reduction in the margin of load-bearing capacity relative to the driving stress. Therefore, a characteristic quantity for weakened load-bearing paths is constructed at the material zoning level. First, the number of material zoning zones remains constant, and the zoning weights are set as volume fractions or volume weights, ranging from 0 to 1, with all weights summed to 1. Simultaneously, the reference time for entering the discrimination interval remains unchanged, and the zoning strength reduction ratio is defined as follows: ; In the formula, , These are two twin mechanical parameters at the reference time; , They are respectively Two twin mechanical parameters at time; for The partition intensity reduction ratio at time.

[0128] Subsequently, the partition strength reduction ratios are weighted and aggregated to obtain the bearing path weakening characteristic quantity, which is used to characterize the overall weakening degree of the bearing mechanism of the entire reactor. When multiple critical partitions degrade simultaneously, the partition strength reduction ratio increases significantly. Then, based on the set thresholds (divided into three set thresholds, the first ranging from 0.05 to 0.30, the second from 0.02 to 0.30, and the third from 2.0 to 6.0), the critical instability pre-state judgment result is given. That is, when the subcritical effective discrimination interval identifier is assigned a value of 1, the shear capacity decay characteristic quantity is greater than or equal to the set threshold ranging from 0.05 to 0.30, the critical approach index is greater than or equal to the set threshold ranging from 0.02 to 0.30, and the bearing path weakening characteristic quantity is greater than or equal to the set threshold ranging from 2.0 to 6.0, it is determined that the shear capacity is indeed decaying, the critical margin is approaching and the decay is accelerating, and the bearing path is weakening as a whole at the partition level. Then, the final critical instability pre-state is output.

[0129] Step S140: Determine whether the target reactor body region is in a critical state of mechanical degradation based on the pre-critical instability state of the reactor body, and output safety and environmental protection early warning information when the target reactor body region is in a critical state of mechanical degradation.

[0130] In some embodiments, the safety and environmental protection early warning method for large-scale industrial solid waste storage sites based on digital twins provided by the present invention includes the following steps in step S140:

[0131] Step S141: Based on the pre-critical instability state of the reactor body and the time-synchronized dual-state mapping structure of the reactor body, define the early warning trigger identifier, and in response to the early warning trigger identifier, perform hierarchical classification processing on the early warning information to obtain multiple early warning levels and the early warning confidence strength corresponding to each level.

[0132] Step S142: Based on the warning level, warning confidence strength, and warning trigger identifier, classify the safety warning type and environmental warning type, and set priorities for the safety warning type and environmental warning type so as to output safety and environmental warning information according to priority.

[0133] In a specific embodiment, the present invention provides a method for early warning of safety and environmental protection for large-scale industrial solid waste stockpiles based on digital twins. Step 4 involves outputting a safety and environmental protection early warning based on the pre-critical instability state. The pre-critical instability state judgment result output in step 3 is used as a direct triggering basis. When it is determined that the stockpile is in a critical state of geometric stability but significant mechanical degradation, corresponding safety and environmental protection early warning information is output to prompt the stockpile manager to take targeted measures before any abnormalities occur in appearance. This achieves early warning of instability and pollution risks in high-moisture fine-particle solid waste stockpiles, including the following sub-steps:

[0134] Sub-step 4.1: Construct the effective identifier for early warning trigger and the feature set required for early warning calculation.

[0135] Specifically, engineering constraints are first applied to determine the effectiveness of the triggering mechanism to avoid false triggering under geometrically unstable or non-target conditions. First, a geometric stability threshold is set to remain constant, ranging from 0.80 to 0.95; a mechanical margin threshold is set, ranging from 0.80 to 1.20, to limit the triggering conditions for significant mechanical degradation. A smaller mechanical state mapping value indicates greater danger. Then, a valid warning trigger identifier is defined as follows: under the final critical instability pre-state, the geometric state mapping value is greater than or equal to the geometric stability threshold, and the mechanical state mapping value is less than or equal to the mechanical margin threshold. The valid warning trigger identifier is assigned the value 1 (0 in other cases). This indicates that the three conditions of a true pre-state determination, stable appearance, and significant mechanical degradation are simultaneously met, ensuring that the output warning truly corresponds to the scenario to be addressed.

[0136] Subsequently, to support the subsequent classification and type output, a feature set required for early warning calculation is constructed, including shear capacity attenuation feature quantity, safety margin, bearing path weakening feature quantity, equivalent shear strength, and equivalent driving shear stress.

[0137] Sub-step 4.2: Early warning level classification and early warning confidence strength generation.

[0138] Specifically, upon effective triggering, the warning is output in a tiered manner. This tiering does not rely on algorithmic classification but rather on three physically meaningful indicators: margin approximation degree, shear attenuation degree, and load-bearing path weakening degree. First, based on the comparison results of the shear capacity attenuation characteristic, safety margin, and load-bearing path weakening characteristic with set thresholds (with values ​​ranging from 0.05-0.15 and 0.15-0.30 respectively), margin thresholds (with values ​​ranging from 0.10-0.25 and 0.00-0.10 respectively), and weakening thresholds (with values ​​ranging from 2.0-3.5 and 3.5-6.0 respectively), the warning level is categorized into three levels: 1-3, as shown below:

[0139] Level 1 warning: The effective warning trigger is a value of 1, the shear capacity attenuation characteristic is greater than or equal to the set threshold in the range of 0.05-0.15, the safety margin is less than or equal to the margin threshold in the range of 0.10-0.25, and the bearing path weakening characteristic is greater than or equal to the weakening threshold in the range of 2.0-3.5.

[0140] Level 2 warning: The effective warning trigger is a value of 1, the shear capacity attenuation characteristic is greater than or equal to the set threshold in the range of 0.15-0.30, the safety margin is less than or equal to the margin threshold in the range of 0.00-0.10, and the bearing path weakening characteristic is greater than or equal to the weakening threshold in the range of 3.5-6.0.

[0141] Level 3 warning level: The effective warning trigger indicator is a value of 1 and the difference between the equivalent shear strength and the equivalent driving shear stress is less than or equal to 0;

[0142] Level 0 warning: The effective indicator for triggering a warning is a value of 0.

[0143] Among them, the Level 3 warning level indicates that the shear capacity has reached the limit state where the driving shear stress is not met. This classification rule is based on the physical margin, and the closer to zero margin, the higher the level.

[0144] Next, the confidence level of the warning is output to quantify the urgency of the warning rather than its accuracy. The expression is as follows: ; In the formula, for The confidence level of the warning at any given moment; for The characteristic quantity of shear capacity decay at time; for The weakened feature quantity of the carrying path at any time; for Safety margin at all times; It is a constant, with a value range of 0.5-2.0.

[0145] Sub-step 4.3: Divide the types of safety warnings into environmental warnings.

[0146] Specifically, to simultaneously output safety and environmental protection early warnings, without introducing algorithms, type attribution is performed based on features with clear physical meanings. The safety side focuses on instability trends, while the environmental side focuses on the critical triggering conditions for leachate overflow and dust risk. First, the water content of the pile related to leachate is determined. If the features required for early warning calculation are not included in the set, they are added according to the same timestamp. Dust risk is assumed to be the pile average of wind speed (range 0-20 m / s) and surface wet-dry ratio (range 0.05-0.70). Then, the safety early warning types are defined as follows:

[0147] When the warning level is greater than or equal to Level 1 and the warning confidence level is greater than or equal to 1.0, it is determined to be a critical instability pre-state warning, indicating that priority should be given to dealing with the slope toe and the local seepage sensitive area;

[0148] When the warning level is Level III, it is determined to be an extreme instability warning, indicating that priority should be given to dealing with the slope toe and the local seepage sensitive area.

[0149] Subsequently, the environmental early warning types were defined as two categories: leachate spill risk and dust risk. The leachate trigger threshold ranges were set to 0.10-0.40 and 1.0-2.0, respectively. The dust trigger threshold range was set to 6-12 m / s, and the average threshold range was set to 0.08-0.20.

[0150] When the subcritical coverage ratio is greater than or equal to the set threshold ranging from 0.10 to 0.40 and the overall state quantity at the stack level is greater than or equal to the set threshold ranging from 1.0 to 2.0, the environmental warning type is determined to be leachate overflow risk, and priority is given to checking the drainage system and the overflow path of the enclosure; when the wind speed is greater than or equal to the dust trigger threshold and the stack average value is less than or equal to the stack average threshold, the environmental warning type is determined to be dust risk, and priority is given to carrying out dust suppression and work surface control.

[0151] Sub-step 4.4: Output safety and environmental protection early warning information package.

[0152] Specifically, the warning level, urgency intensity, and warning type are organized into actionable warning information packages, and suggested handling priorities are given to prompt yard managers to initiate corresponding measures before any abnormalities are observed in the appearance. The priority value is set to a range of 1-5, with a higher priority value indicating a more urgent risk.

[0153] The following describes the safety and environmental protection early warning device for bulk industrial solid waste storage sites based on digital twins provided by this invention. The safety and environmental protection early warning device for bulk industrial solid waste storage sites based on digital twins described below can be referred to in correspondence with the safety and environmental protection early warning method for bulk industrial solid waste storage sites based on digital twins described above.

[0154] like Figure 2 As shown in one embodiment, a safety and environmental protection early warning device for a bulk industrial solid waste storage site based on digital twins includes a stack state construction module, a model construction module, a mechanical state analysis module, and an early warning information output module.

[0155] The stack state construction module is used to obtain the stack state quantities within the target stack yard and construct a unified time-based stack state set based on the stack state quantities.

[0156] The model building module is used to construct a dual-state mapping structure of the heap in the digital twin space based on the heap's comprehensive state set, thereby obtaining a dual-state digital twin heap model.

[0157] The mechanical state analysis module is used to call the dual-state digital twin reactor model to analyze the target reactor region, identify the mechanical state characteristics of the reactor, and determine the critical instability pre-state of the reactor based on the mechanical state characteristics.

[0158] The early warning information output module is used to determine whether the target reactor area is in a critical state of mechanical degradation based on the pre-critical instability state of the reactor body, and outputs safety and environmental protection early warning information when the target reactor area is in a critical state of mechanical degradation.

[0159] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0160] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of this patent should be determined by the appended claims.

Claims

1. A method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins, characterized in that, The method includes: Obtain the state variables of the stack within the target stack, and construct a unified time-based integrated state set of the stack based on the stack state variables; Based on the aforementioned integrated state set of the heap, a dual-state mapping structure of the heap is constructed in the digital twin space to obtain a dual-state digital twin heap model. The dual-state digital twin reactor model is invoked to analyze the target reactor region in order to identify the reactor mechanical state characteristics, and the reactor critical instability pre-state is determined based on the reactor mechanical state characteristics; Based on the pre-critical instability state of the reactor body, it is determined whether the target reactor body region is in a critical state of mechanical degradation, and when the target reactor body region is in a critical state of mechanical degradation, safety and environmental protection early warning information is output.

2. The method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 1, characterized in that, The process of acquiring the stack state quantities within the target stack yard and constructing a unified time-based integrated stack state set based on these stack state quantities includes: Multiple monitoring points are selected within the target stockpile and numbered to obtain the water content, pore pressure, and apparent morphology of the stockpile at each monitoring point according to a set sampling period. Based on the water content, pore pressure, and apparent morphology of the pile body, and combined with the preset critical reference values, the comprehensive state quantity of each monitoring point at each time is calculated. The stack state quantities include the stack water content state quantity, pore pressure state quantity, and apparent morphology state quantity.

3. The method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 2, characterized in that, The step of acquiring the stack state quantities within the target stack yard and constructing a unified time-based integrated stack state set based on the stack state quantities further includes: The comprehensive state quantities corresponding to each monitoring point at each time are summarized into the comprehensive state quantity of the reactor body, and a subcritical criterion threshold is set to determine the subcritical coverage ratio when the comprehensive state quantity is not lower than the subcritical criterion threshold. Based on the overall state quantity of the reactor body and the subcritical coverage ratio, a reactor body strength threshold and a coverage ratio threshold are set, and a transition interval state identifier is determined when the overall state quantity of the reactor body is not lower than the reactor body strength threshold and the subcritical coverage ratio is not lower than the coverage ratio threshold. The integrated state set of the reactor body is composed of the integrated state quantity of the reactor body, the subcritical coverage ratio, and the transition interval state identifier.

4. The method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 1, characterized in that, Based on the comprehensive state set of the heap, a dual-state mapping structure for the heap is constructed in the digital twin space to obtain a dual-state digital twin heap model, including: A geometric mapping channel is established in the digital twin space, and a geometric mapping reference and a geometric mapping effective area are established for the integrated state set of the heap through the geometric mapping channel; Based on the geometric mapping datum and the effective area of ​​geometric mapping, the surface of the digital twin stack is divided into multiple partition units, and the elevation of the digital twin surface of each partition unit is determined. Each monitoring point in the effective area of ​​the geometric mapping is projected to the corresponding partition unit, and the equivalent shape update amount of the partition unit is calculated. The elevation of the digital twin surface is updated by the equivalent shape update amount to obtain the geometric state mapping amount and the geometric twin surface state.

5. The method for safety and environmental protection early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 4, characterized in that, The step of constructing a dual-state mapping structure for the heap in the digital twin space based on the heap's comprehensive state set to obtain a dual-state digital twin heap model further includes: The target pile body region is divided into multiple material partitions, and the reference shear strength parameters of each material partition are obtained. At the same time, the mechanical degradation coefficient of each material partition is defined according to the comprehensive state set of the pile body. The twin mechanical parameters are obtained by equivalently reducing the reference shear strength parameters using the mechanical degradation coefficient. Obtain the equivalent values ​​of normal stress on the slope of the pile body, void pressure, and weighted average value of friction force, and calculate the equivalent shear strength and equivalent driving shear stress based on the equivalent values ​​of normal stress on the slope of the pile body, void pressure, and weighted average value of friction force. The mechanical state mapping quantity is calculated based on the equivalent shear strength and equivalent driving shear stress, and the twin mechanical parameters, equivalent shear strength and equivalent driving shear stress are integrated into a mechanical twin bearing state; The dual-state digital twin stack model is composed of a stack dual-state mapping structure, which is composed of the time-synchronized mechanical state mapping quantity and the mechanical twin bearing state.

6. The method for safety and environmental early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 5, characterized in that, The step of calling the dual-state digital twin reactor model to analyze the target reactor region to identify the reactor's mechanical state characteristics, and determining the reactor's critical pre-instability state based on these characteristics, includes: Based on the aforementioned dual-state mapping structure and integrated state set of the stack, a transition interval identifier for subcritical effective discrimination is constructed, and the geometric stability constraint is calculated according to the set geometric stability threshold. Obtain the start time of the transition interval identifier, define the shear capacity attenuation characteristic based on the geometric stability constraint and the equivalent shear strength at the start time of the interval, and calculate the shear capacity attenuation rate based on the shear capacity attenuation characteristic and the set sampling period.

7. The method for safety and environmental protection early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 6, characterized in that, The step of calling the dual-state digital twin reactor model to analyze the target reactor region to identify the reactor mechanical state characteristics, and determining the pre-critical instability state of the reactor based on the reactor mechanical state characteristics, further includes: The difference between the equivalent shear strength and the equivalent driving shear stress is normalized to obtain the critical equilibrium characteristic quantity, and the critical approximation index is calculated based on the critical equilibrium characteristic quantity and the rate of change of shear capacity decay. Based on the critical equilibrium characteristic quantity and the critical approximation index, the partition strength reduction ratio corresponding to each material partition is defined, and the partition strength reduction ratio is weighted and summarized according to the geometric stability constraint quantity to obtain the bearing path weakening characteristic quantity. When the bearing path weakening characteristic, critical approximation index, and shear capacity attenuation characteristic are all not lower than their respective set thresholds, the pre-critical instability state of the reactor body is output.

8. The method for safety and environmental protection early warning of large-scale industrial solid waste storage sites based on digital twins according to claim 1, characterized in that, The step of determining whether the target reactor region is in a critical state of mechanical degradation based on the pre-critical instability state of the reactor body, and outputting safety and environmental protection early warning information when the target reactor region is in a critical state of mechanical degradation, includes: Based on the pre-critical instability state of the reactor and the time-synchronized dual-state mapping structure of the reactor, a warning trigger identifier is defined, and the warning information is classified and processed in response to the warning trigger identifier to obtain multiple warning levels and the warning confidence strength corresponding to each level. Based on the warning level, warning confidence strength, and warning trigger identifier, safety warning types and environmental warning types are divided, and priorities are set for the safety warning types and environmental warning types so that the safety and environmental warning information is output according to the priorities.

9. A terminal, comprising a processor and a storage medium; characterized in that: The storage medium is used to store instructions; The processor is configured to operate according to the instructions to perform the steps of the method according to any one of claims 1-8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the program implements the steps of the method according to any one of claims 1-8.

Images