A method for dynamic monitoring and control of the spatial environment of underground utility tunnels

Abstract

This invention discloses a method for dynamic monitoring and control of the spatial environment in underground utility tunnels, relating to the field of urban infrastructure operation and maintenance technology. After sorting the tunnel segments, instead of applying the same control method directly to all abnormal segments, it continues to read the environmental monitoring items that contribute the most to the deviation in each segment. Combining the function type, function area, and cross-segment boundaries of existing actuators, it calculates the control mapping relationship between the tunnel segment and candidate actuators, and then selects the actuator that best matches the current dominant deviation state for each segment. Fans, drainage pumps, access control linkage programs, and ventilation zone switching programs are no longer simply called in fixed patterns, but are instead linked separately based on the dominant environmental problems of different tunnel segments. For common issues in underground utility tunnels such as gas migration, moisture diffusion, and water accumulation, this method enables the actuator to match the abnormal attributes and main chain relationships of the current tunnel segment.

Description

Technical Field

[0001] This invention relates to the field of urban infrastructure operation and maintenance technology, specifically to a method for dynamic monitoring and control of the spatial environment of underground integrated pipe corridors. Background Technology

[0002] The operation and management of underground integrated utility tunnels falls under the technical field of urban infrastructure operation and maintenance, and further involves the technical field of environmental status monitoring and coordinated response in underground linear enclosed spaces. Specifically, it involves a dynamic monitoring and control technology for the spatial environment, which continuously collects, determines the status of, identifies the propagation of, and controls spatial environmental parameters such as temperature, humidity, oxygen concentration, harmful gas concentration, and water level within the underground integrated utility tunnel. This type of technology does not target a single monitoring device or a single environmental quantity, but rather a continuous corridor space composed of multiple compartments, connecting doors, ventilation zones, drainage slopes, and monitoring points. The monitoring results are not only related to the values ​​of each monitoring point itself, but also to the connectivity of the compartments, airflow transmission, slope orientation, and the distribution of maintenance access points.

[0003] Current environmental monitoring methods for underground utility tunnels mostly revolve around individual equipment points. Typically, data such as temperature, humidity, oxygen concentration, harmful gas concentration, and water level are collected separately at each monitoring point. The data from each monitoring point is then compared to its respective preset threshold, triggering alarms or coordinated responses when the threshold is reached. While this method can reflect the numerical changes at a specific monitoring point at a specific moment, in practical applications, environmental anomalies within underground utility tunnels do not always appear in isolated, single-point forms. Instead, they often propagate segment by segment along the structural connections between adjacent sections. Especially when ventilation zones are continuous, drainage slopes are consistent, or maintenance openings are interconnected, a localized environmental anomaly often extends from the initial segment to adjacent segments. Existing point-based judgment methods lack an organizational process for the spatial relationships between tunnel segments, and also lack a joint judgment process for the sequence of anomalies appearing between adjacent segments, the extent of their spread, and the areas where they linger. Therefore, it is difficult to identify the starting point, propagation path, and main lingering segments of the same anomaly chain from multiple discrete alarm points.

[0004] The aforementioned situation arises from two main factors. First, existing monitoring systems primarily use sensor placement as the basic organizational unit, directly mapping data to individual device locations. Second, current judgment logic often employs independent threshold comparisons, failing to adequately utilize spatial information such as structural diagrams, compartment numbers, connecting door locations, ventilation zone boundaries, and drainage slopes. This results in monitoring results appearing as numerous scattered numerical and alarm records at the data level. When issues such as moisture diffusion, localized oxygen deficiency, harmful gas migration, or water accumulation along the slope occur within the underground utility tunnel, the system may generate alarms at multiple points sequentially. However, it cannot pinpoint which section the anomaly originated from, its direction of propagation, or which sections it lingers in. This easily leads to misjudging continuously propagating anomalies as multiple independent events, resulting in inconsistencies between the control sequence and the actual propagation order of the anomaly. The relationship between the initial and subsequent responses becomes confused, and residual anomalies may remain in adjacent sections even after localized intervention. Consequently, the anomaly chain persists within the utility tunnel, impacting subsequent inspections, ventilation, drainage, and compartment operation management. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this invention provides a method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels, solving the problems mentioned in the background art.

[0006] To achieve the above objectives, the present invention provides the following technical solution: a method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels, comprising the following steps:

[0007] S1. Read the structural data of the underground utility tunnel, divide the tunnel space into segments, and construct a segment topology diagram consisting of multiple segment nodes and the connection relationship between adjacent segments.

[0008] S2. Based on the corridor topology diagram, read the existing environmental monitoring data in the corridor and perform unified time axis mapping and corridor-level summary processing to generate the structured monitoring sequence corresponding to each corridor and calculate the comprehensive deviation E(n) corresponding to each corridor.

[0009] S3. Based on the corridor topology diagram and combined with the comprehensive deviation E(n), calculate the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path, and obtain the path-level spatial anomaly value Mp and the anomaly propagation main chain.

[0010] S4. Based on the main chain of abnormal propagation, calculate the control sequence value Kh corresponding to each corridor segment, and sort each corridor segment in the main chain according to the control sequence value to generate the control execution sequence.

[0011] S5. After executing the control sequence, continue to read the environmental monitoring data corresponding to each corridor segment in the anomaly propagation main chain, obtain the structural regression value Yp, and perform feedback control.

[0012] Preferably, S1 includes S11 and S12;

[0013] S11. Read the corridor centerline, compartment outline, compartment number, connecting door position, ventilation zone boundary, drainage slope change position, environmental monitoring point position and actuator layout data from the pipe gallery structure diagram or model data, and sort each structural position in a chain according to the extension direction of the corridor centerline to form an ordered structural position sequence.

[0014] The boundary criteria are as follows: the continuity of numbering between adjacent structural locations, changes in connecting components, switching of ventilation zones, changes in slope aspect, and mixed affiliation of monitoring points. The boundary value of each adjacent location is calculated using these criteria. The formula for obtaining the boundary value is:

[0015]

[0016] In the formula, Ba represents the corridor boundary value between the a-th position and the (a+1)-th position in the ordered structural position sequence; Na represents the numbering break mark, used to characterize whether the compartment number of two adjacent structural positions changes; it is 1 when a numbering change occurs and 0 when no change occurs; Ga represents the number of connecting components, used to characterize the number of doors, inspection ports, partitions, or grille openings that appear between the current positions; Va represents the air zone switching mark, used to characterize whether two adjacent structural positions are in different ventilation zones; it is 1 when a switch occurs and 0 when no switch occurs; Pa represents the slope reversal mark, used to characterize whether the drainage slope changes at this position; it is 1 when the slope direction changes and 0 when no change occurs; La represents the axial distance between two adjacent structural positions along the centerline direction; and Ma represents the number of monitoring points that are simultaneously connected to different functional objects between two adjacent structural positions.

[0017] The corridor segment boundary value Ba is compared with the corridor segment dividing baseline BX to obtain the corridor segment node set; the specific method is as follows:

[0018] When the segment boundary value Ba is greater than or equal to the segment dividing baseline BX, the current position is determined as the segment boundary point, and the ordered structure sequence is segmented to form multiple segment nodes.

[0019] The baseline for dividing the corridor segment, BX, is obtained as follows: BX = Med(Ba) + Med(|Ba - Med(Ba)|).

[0020] Med(Ba) represents the median value of all corridor segment boundary values, and Med(|Ba-Med(Ba)| represents the median dispersion of all corridor segment boundary values ​​relative to the median value.

[0021] Then, the starting position, ending position, cabin type, internal monitoring point set, and internal actuator set corresponding to each corridor node are archived to generate a corridor node set.

[0022] Preferably, in step S12, based on the set of corridor nodes, the doorway connectivity, ventilation direction consistency, drainage slope bearing capacity, and centerline adjacency distance between each corridor node are read respectively.

[0023] Determine whether there is a direct adjacency relationship between any two corridor segment nodes; for corridor segment nodes with a direct adjacency relationship, further calculate the topological connectivity Cij, obtaining the formula as follows:

[0024]

[0025] In the formula, Cij represents the topological connectivity between corridor node i and corridor node j; Kij represents the interface connectivity number, which is used to characterize the number of exchangeable interfaces such as doorways, inspection ports, and grilles between two corridor nodes; Rij represents the wind direction collinearity symbol, which is used to characterize whether the dominant ventilation direction is continuously extended along the connection direction of the two corridor nodes, taking 1 when they are consistent and 0 when they are inconsistent; Sij represents the slope bearing symbol, which is used to characterize whether the drainage slope direction is from corridor node i to corridor node j, taking 1 when there is a merging relationship and 0 otherwise; Dij represents the centerline adjacency distance between two corridor nodes; Hij represents the number of barrier layer levels, which is used to characterize the number of partition doors, fireproof partition units, or other barrier components between two corridor nodes;

[0026] When the topological connectivity Cij value is large, it indicates that there is both an actual interface connection between the two corridor nodes and the conditions for wind flow or slope connection, which are high-association adjacent segments in the subsequent environmental anomaly propagation analysis; when the topological connectivity Cij value is small, it indicates that although they are close in space, they should not be directly regarded as high-priority propagation adjacent segments due to the presence of barrier components or discontinuous wind slope conditions.

[0027] After calculating the topological connectivity Cij of all adjacent corridor nodes, the corridor node pairs that meet the direct adjacency condition are written into the adjacency relationship set. The corridor nodes are used as graph nodes, the adjacency relationships are used as graph edges, and the topological connectivity Cij is used as the edge attribute to construct the corridor topological relationship graph.

[0028] Preferably, S2 includes S21 and S22;

[0029] S21. Based on the corridor topology diagram, read the temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each monitoring point inside each corridor, and simultaneously read the collection start time, collection end time, record holding time and data type identifier corresponding to the original collection records of each monitoring point; due to the differences in data collection cycle, reporting interval and holding time of different monitoring points.

[0030] A unified time-slice sequence is constructed within the current analysis period, dividing the entire analysis period into multiple consecutive time slices, and then the original collection records of each monitoring point are mapped to the corresponding unified time slice;

[0031] The mapping process is as follows: for each original record, extract its start time, end time and corresponding monitoring value, and determine whether there is a time overlap between the coverage time interval of the original record and each unified time slice;

[0032] When there is a time overlap, the overlap duration between the original record and the corresponding unified time slice is calculated, and the overlap duration is used as the coverage length of the original record to the unified time slice; then the monitoring value of the original record is multiplied by the corresponding overlap duration to obtain the duration contribution value of the original record to the unified time slice, and the duration contribution values ​​of all original records falling within the same unified time slice are accumulated.

[0033] Meanwhile, the overlapping duration of all original records within the unified time slice is accumulated; finally, the total accumulated duration contribution value is divided by the total accumulated overlapping duration to obtain the mapping result of the monitoring point in the unified time slice.

[0034] After completing the unified time axis mapping of each monitoring point, the mapping values ​​of the same type of monitoring points in the same corridor segment on the same time slice are then summarized at the corridor segment level to form the structured monitoring sequence S={(T(n),U(n),O(n),G(n),W(n))|n=1,2,...,N} corresponding to the corridor segment;

[0035] In the formula, T(n), U(n), O(n), G(n) and W(n) represent the representative values ​​of temperature, humidity, oxygen concentration, harmful gas concentration and water level of the corridor segment in the nth unified time slice, respectively, and N represents the total number of unified time slices obtained in the current analysis period.

[0036] Preferably, S22, based on the structured monitoring sequence S, samples of the same type of operating conditions consistent with the current operating conditions are extracted from the historical database, and the temperature reference value, humidity reference value, oxygen concentration reference value, harmful gas concentration reference value and water level reference value of each corridor section under the same type of operating conditions are obtained respectively.

[0037] Simultaneously, the historical swing width of various monitoring quantities around their respective benchmark values ​​is obtained; then, the obtained structured monitoring sequence of the corridor segment is compared with the historical benchmark data on a time-by-time basis, and the comprehensive deviation E(n) corresponding to each corridor segment is calculated.

[0038] The overall deviation E(n) is obtained using the following formula:

[0039]

[0040] In the formula, Ei(n) represents the comprehensive deviation of corridor segment i in the nth unified time slice, ln represents the logarithmic function, Sji(n) represents the j-th data in the structured monitoring sequence of corridor segment i in the nth unified time slice, pSji represents the historical baseline value of the j-th data in the structured monitoring sequence of corridor segment i, and bSji represents the historical swing width of the j-th data in the structured monitoring sequence of corridor segment i.

[0041] When the value of Ei(n) is small, it indicates that the current environmental state of the corridor is close to the normal distribution under similar operating conditions; when the value of Ei(n) is large, it indicates that the corridor has deviated from its historical normal state in the current time slice and should be the focus of subsequent abnormal path identification.

[0042] The comprehensive deviation E(n) is combined with the historical comprehensive deviation sample sequence {Eih(z)} of the corridor segment under the same operating conditions in the historical database to calculate and obtain the dynamic judgment line EiD corresponding to the corridor segment;

[0043] The formula for obtaining the dynamic decision line EiD is: EiD = Med{Eih(z)} + Med{|Eih(z)-Med{Eih(z)}|};

[0044] In the formula, Eih(z) represents the z-th historical comprehensive deviation sample of corridor segment i under the same operating conditions, and Med{Eih(z)} represents the median value of the historical comprehensive deviation sample.

[0045] The dynamic judgment line EiD is not a fixed threshold, but a judgment boundary that changes with the historical environmental fluctuations of the corridor segment itself.

[0046] When a certain corridor segment has historically experienced relatively stable fluctuations, the resulting dynamic judgment line EiD is lower, making it easier to identify current deviations.

[0047] When a certain corridor segment has historically experienced relatively active fluctuations, the resulting dynamic judgment line EiD is higher, and only deviations exceeding the historical fluctuation background will be identified as abnormal.

[0048] Therefore, the comprehensive deviation E(n) in the current time slice is compared with the dynamic decision line EiD to determine whether the corridor segment has entered a deviation state in the current time slice; the time slice that satisfies Ei(n)≥EiD is taken as the time slice of the corridor segment that exceeds the line, and the corresponding result is sent to the subsequent abnormal propagation path identification step for use;

[0049] Preferably, S3 includes S31 and S32; S31: For any corridor segment i, the time slice that satisfies Ei(n)≥EiD is taken as the time slice that crosses the line for that corridor segment, and the time slice number that first satisfies the condition is recorded as the first time ti that the corridor segment crosses the line.

[0050] The total duration covered by the longest consecutively satisfied segment is denoted as the continuous dwell time Qi of the corridor segment; simultaneously, the over-line amplitudes on all over-line time segments of the corridor segment are aggregated to obtain the average over-line amplitude A of the corridor segment, obtained by the following formula:

[0051]

[0052] In the formula, Ai represents the average overshoot amplitude of corridor segment i; SAi represents the set of overshoot time slices of corridor segment i; and Ni+ represents the number of time slices in the set of overshoot time slices.

[0053] The average deviation is used to represent the overall degree of elevation of the corridor segment after it enters the deviation state, which deviates from its historical normal boundary. When the Ai value is large, it indicates that the corridor segment exceeds the dynamic judgment line by a large margin and should be regarded as a key observation segment in the subsequent path identification process. When the Ai value is small, it indicates that although the corridor segment has entered the deviation state, its deviation is still in a relatively shallow range.

[0054] After obtaining the first crossing time ti, continuous dwell time Qi, and average crossing amplitude A of each corridor segment, corridor segments with direct adjacent relationships and all having crossing time slices are taken as path connection objects and connected in series along the adjacent connection direction in the corridor segment topology diagram to form several candidate propagation paths p=(h1, h2, ..., hm);

[0055] Where p represents a candidate propagation path; hr represents the r-th segment in path p arranged in connection order; and m represents the number of segments contained in the candidate propagation path.

[0056] For any two adjacent corridor segments hr and hr+1 in the candidate propagation path, first calculate the correction order interval Δtr. When the first crossing time of the latter corridor segment is not earlier than that of the former corridor segment, take the time difference between the two. When a time reversal occurs after the former corridor segment, add an equal amount of penalty to the time difference, specifically obtained by the following formula:

[0057] ;

[0058] Wherein, Δtr specifically represents the correction order interval between the r-th corridor segment and the (r+1)-th corridor segment in the candidate propagation path; t(hr) and t(hr+1) represent the first crossing time of two adjacent corridor segments, respectively;

[0059] Based on the correction order interval Δtr and the average superline amplitude A, the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path are calculated and obtained.

[0060] The method for obtaining the transmission order quantity Op is as follows: First, along the arrangement order of the candidate propagation path, take out two adjacent corridor segments in the path one after the other; for each pair of adjacent corridor segments, read the topological coupling degree between them, and at the same time read the first crossing time of the previous corridor segment and the next corridor segment; then, divide the topological coupling degree between the adjacent corridor segment pairs by "1 and the sum of the successive intervals" to obtain the local contribution value of the adjacent corridor segment pairs to the transmission order quantity of the entire candidate propagation path; then, according to the adjacent order in the path, perform the same processing on all adjacent corridor segment pairs respectively, and accumulate all the obtained local contribution values ​​one by one to obtain the transmission order quantity Op of the candidate propagation path;

[0061] The diffusion span measure Up is obtained as follows: First, according to the order of the candidate propagation paths, the axial length and the corresponding average overshoot of each corridor segment on the path are read sequentially; for each corridor segment, the axial length of the corridor segment is combined with its average overshoot to obtain the single-segment contribution value of the corridor segment to the diffusion range of the entire candidate propagation path; then, the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the diffusion span measure Up of the candidate propagation path.

[0062] The dwell intensity Jp is obtained by sequentially reading the continuous dwell time and the corresponding average overshoot of each corridor segment on the path; for each corridor segment, the continuous dwell time and the average overshoot are combined to obtain the single-segment contribution value of the corridor segment to the dwell state of the entire candidate propagation path; then the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the dwell intensity Jp of the candidate propagation path.

[0063] Preferably, in step S32, for each candidate propagation path obtained in step S31, the propagation order quantity Op, the diffusion span quantity Up, and the dwell intensity quantity Jp are jointly calculated to obtain the path-level spatial anomaly value Mp of the candidate propagation path.

[0064] Path-level spatial outliers Mp are obtained in the following way:

[0065] ;

[0066] When a path has a high transmission order, a large diffusion span, and a strong dwell intensity, the resulting Mp value is large, indicating that the path is closer to the actual propagation chain of the current spatial anomaly inside the utility tunnel; when a path has only one of these characteristics and the other two are weak, the resulting Mp value is small.

[0067] After obtaining the path-level spatial outliers of all candidate propagation paths, the path-level spatial outliers Mp of each path are compared, and the candidate propagation path with the largest value is taken as the current anomaly propagation main chain.

[0068] Subsequently, within the main chain of abnormal propagation, the initial crossing time, continuous dwell time, and average crossing amplitude of each corridor segment were read, and the dwell contribution value Dhr of each corridor segment was calculated:

[0069] The dwell contribution value Dhr is obtained by multiplying the continuous dwell time of the corridor segment by the average over-line amplitude of the corridor segment;

[0070] The dwell contribution value is used to describe the degree to which a corridor segment can withstand abnormal dwelling states within the main chain. When a corridor segment maintains both a long continuous overrun duration and a high overrun amplitude, the dwell contribution value Dhr is larger, and the corridor segment is closer to an abnormal dwelling segment. When a corridor segment has entered the main chain, but the dwell time is short or the overrun amplitude is shallow, the dwell contribution value Dhr is smaller.

[0071] Within the main chain, the corridor segment with the earliest first crossing time is designated as the starting corridor segment; the corridor segment with a dwell contribution value Dhr greater than or equal to the median dwell contribution value of all corridor segments in the main chain is designated as the dwell corridor segment; and the remaining main chain corridor segments located after the starting corridor segment and not belonging to the dwell corridor segment are designated as the diffusion corridor segments.

[0072] Preferably, S4 includes S41 and S42; S41, according to the arrangement order in the main chain of abnormal propagation, extract the continuous dwell time, average over-line amplitude, topological coupling degree between each corridor segment and the upstream adjacent corridor segment, the shortest span to the actionable device and the segment-level distance of the corridor segment relative to the starting corridor segment.

[0073] Among them, the upstream coupling relationship of the starting corridor section is handled according to the main chain entrance mark;

[0074] For any corridor segment h in the main chain, first determine its upstream adjacent corridor segment uh, and then calculate the control sequence value Kh of the corridor segment.

[0075] The control sequence value Kh is obtained as follows: the continuous dwell time of the corridor segment, the average overshoot amplitude, and the topological coupling degree with the upstream corridor segment are taken as the sequential lifting term. The shortest span number from the corridor segment to the existing actuator and the segment-level distance relative to the starting corridor segment are taken as the sequential pressure drop term. The result of the product of the sequential lifting terms is divided by the sum of the sequential pressure drop term and 1 to obtain the control sequence value Kh corresponding to the corridor segment.

[0076] When the control sequence value Kh is large, it indicates that the corridor segment is closer to the main abnormal segment in terms of dwell state, over-line depth, and propagation connection relationship, and is more suitable to prioritize calling the existing execution device; when the control sequence value Kh is small, it indicates that the corridor segment is closer to the subsequent linkage segment or the segment affected by delay.

[0077] Preferably, in step S42, based on the control sequence value Kh, all corridor segments in the main chain are sorted from largest to smallest to obtain the control sequence within the main chain.

[0078] When two or more corridor segments have the same control priority value, the corridor segment with the smaller segment distance is selected first, and then the corridor segment with the larger average over-line amplitude is selected as the target of the advance control.

[0079] After completing the corridor segment sorting, the environmental monitoring items with the largest deviation contribution in each corridor segment, as well as the action type, action area and cross-segment action boundary of each existing actuator, are read. For each sorted corridor segment, an actuator or control program that can act on the corridor segment and is consistent with the dominant deviation type of the corridor segment is selected, and the regulation matching value between the corridor segment and the candidate actuator is calculated.

[0080] The method for calculating the control matching value is as follows: take the consistent state of the actuator and the current corridor segment's dominant deviation type, and the continuous coverage state of the actuator's action section on the current corridor segment and its downstream adjacent corridor segments as the matching lifting term, take the number of segments spanned from the actuator to the current corridor segment as the matching pressure drop term, and then divide the result of the multiplication of the matching lifting term by the sum of the matching pressure drop term and 1 to obtain the control matching value between the current corridor segment and the corresponding actuator.

[0081] For each sorted corridor segment, the control program with the largest control matching value is selected, and the "corridor segment number - execution object - execution order" are sequentially connected to form a control execution sequence.

[0082] Preferably, S5 includes S51 and S52; S51: After the control execution sequence is completed, a retest time window covering the observation period after control is selected to obtain temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each corridor segment in the main chain of abnormal propagation.

[0083] Following the unified time axis mapping method and corridor segment-level aggregation method in step S2, the corridor segment representative temperature, corridor segment representative humidity, corridor segment representative oxygen concentration, corridor segment representative harmful gas concentration and corridor segment representative water level within the retest time window after regulation are obtained, and the comprehensive deviation of each corridor segment after regulation on each unified time slice is recalculated.

[0084] For any segment in the main chain of abnormal propagation, the time slice that still meets the condition that "the comprehensive deviation is not lower than the dynamic judgment line of the segment" after regulation is taken as the residual over-line time slice of the segment; the total duration corresponding to the longest continuous segment in the residual over-line time slice is recorded as the remaining dwell time Qhr of the segment; the over-line portion on the residual over-line time slice is accumulated and converted according to the number of residual over-line time slices to obtain the remaining over-line amplitude Ahr of the segment;

[0085] The process of obtaining the remaining over-line amplitude Ahr is as follows: First, read all the remaining over-line time slices of the corridor segment after regulation, and extract the difference between the comprehensive deviation on each remaining over-line time slice and the dynamic judgment line of the corridor segment; then accumulate all the differences; finally, divide the total accumulated difference by the number of remaining over-line time slices to obtain the average over-line depth that the corridor segment still retains after regulation.

[0086] Therefore, the remaining over-line amplitude Ahr is used to characterize the degree to which the corridor segment still deviates from the historical normal boundary of the segment after the current regulation; when the value of the remaining over-line amplitude Ahr is large, it indicates that although the corridor segment has been regulated, it still retains a relatively obvious abnormal state; when the value of the remaining over-line amplitude Ahr is small, it indicates that the deviation state of the corridor segment is close to exiting the over-line range.

[0087] After obtaining the remaining dwell time Qhr and the remaining overshoot amplitude Ahr of each corridor segment in the main chain, the order of the remaining time within the main chain after regulation is re-examined according to the original arrangement of each corridor segment in the main chain.

[0088] If the residual first overrun time of the downstream corridor segment is earlier than the residual first overrun time of the upstream corridor segment, then the adjacent corridor segments are recorded as a residual back jump segment.

[0089] After checking all adjacent corridor pairs, the number of residual bounce segments Bp in the current abnormal propagation main chain is obtained; at the same time, all corridor segments in the main chain are checked one by one to see if there are still residual time slices, and the number of corridor segments with residual time slices is recorded as the number of residual time slices Sp.

[0090] Among them, the residual bounce number Bp is used to characterize whether there is still an abnormal sequence reversal phenomenon inside the main chain after regulation; when the residual bounce number Bp is large, it indicates that there is still a discontinuous residual propagation relationship inside the main chain after regulation.

[0091] The number of residual superline corridor segments Sp is used to characterize the coverage area of ​​corridor segments that still retain the superline state in the current main chain; when the number of residual superline corridor segments Sp is large, it indicates that the current anomaly has not been reduced to a small number of local corridor segments;

[0092] By combining the continuous dwell time Qi and average over-line amplitude A of each corridor segment in the main chain before regulation, and the remaining dwell time Qhr, remaining over-line amplitude Ahr, remaining bounce segment number Bp and remaining over-line corridor segment number Sp obtained after regulation, the structural regression value Yp corresponding to the current abnormal propagation main chain is obtained.

[0093]

[0094] In the formula, m represents the number of segments contained in the main chain of abnormal propagation;

[0095] S52. First, extract completed control events from the historical event database that are identical to the current abnormal propagation main chain in terms of cabin type, dominant deviation type, main chain length, and calling execution device type, and read the structural regression value sample corresponding to the historical event at the end of the control; then, obtain the regression decision line YpD of the current abnormal propagation main chain based on the structural regression value sample.

[0096] The regression decision line YpD is obtained as follows: First, sort the structural regression value samples of similar historical regulation events and extract their median values; then calculate the dispersion of each sample relative to the median value; finally, subtract the dispersion from the median value to obtain the regression decision line YpD corresponding to the current abnormal propagation main chain.

[0097] The regression threshold is used to indicate the structural fallback boundary that the current main chain should reach under the background of similar regulatory events. When the structural regression value of a certain event is not lower than the regression threshold, it indicates that the current main chain's dwell state, over-limit depth, and residual chain range have entered the regression range of similar events. When the structural regression value of the event is lower than the regression threshold, it indicates that the current main chain still retains a lot of residual structure.

[0098] Compare the structural regression value Yp with the regression decision line YpD to determine the current state of the anomaly propagation main chain.

[0099] The judgment method is as follows:

[0100] When the structural regression value Yp ≥ the regression decision line YpD, the current abnormal propagation main chain is determined to enter the regression state, the current control result is output, and the control end time, regression corridor set, residual corridor set and final structural regression value corresponding to the current main chain are written into the historical event database.

[0101] When the structural regression value Yp < regression decision line YpD, it is determined that the current abnormal propagation main chain has not entered the regression state. Further, the corridor segments that still have residual overshoot time slices in the current main chain are read, and the remaining dwell time Qhr, remaining overshoot amplitude Ahr, residual bounce relationship and residual segment position of the corridor segments are sent back to step S4 as inputs for the next round of regulation order calculation and regulation execution sequence rearrangement.

[0102] This invention provides a method for dynamic monitoring and control of the spatial environment of underground utility tunnels, which has the following beneficial effects:

[0103] (1) Based on the structured monitoring sequence, the comprehensive deviation status of each corridor segment is continuously analyzed, and the transmission sequence, diffusion span, and dwell intensity of the candidate propagation path are further calculated. Then, the path-level spatial anomaly value and the main chain of anomaly propagation are obtained. Based on this processing, the system output is no longer a number of isolated alarm points, but a spatial anomaly chain that can distinguish the anomaly initiation segment, diffusion segment, and dwell segment. This is helpful to identify the propagation direction, coverage area, and main dwelling segment of the same environmental anomaly from multiple discrete alarms, and solves the problem that existing methods cannot explain where the anomaly starts, along which path it expands, and in which segments it persists.

[0104] (2) In the process of dividing the corridor into sections, not only the section number and spatial location are considered, but also the connecting doors, ventilation zones, drainage slopes, and the mixed situation of monitoring point affiliation are considered. This makes the determination of the corridor boundary no longer rely solely on a single location marker, but rather on a comprehensive judgment based on the structural connectivity and environmental transmission conditions. For common situations in underground integrated pipe corridors such as moisture migration, gas diffusion, local hypoxia, and water accumulation spreading along the structural path, this division method is closer to the actual propagation boundary. It can refine the originally difficult-to-distinguish continuous space into multiple analysis units with environmental significance, making it easier to identify where the anomaly starts, where it is segmented, and where it continues to the next section.

[0105] (3) By further organizing the historical comprehensive deviation samples, dynamic judgment lines corresponding to each corridor segment are obtained, so that the judgment boundaries of different corridor segments correspond to their respective historical fluctuations. For corridor segments with relatively stable historical fluctuations, their abnormal changes are easier to identify; for corridor segments with relatively active historical fluctuations, it is necessary to combine their own historical distribution to determine whether the current state is abnormal. This helps to solve the problem that the unified threshold in the existing method is difficult to adapt to different corridor segments and different working conditions.

[0106] (4) By constructing transmission sequence quantity, diffusion span quantity and residence intensity quantity respectively, the temporal sequence relationship, spatial extension range and abnormal residence state of candidate propagation paths are jointly judged, which can distinguish between two different situations: "only multiple corridor segments fluctuate at the same time" and "the anomaly is transmitted segment by segment along adjacent corridor segments". For continuous propagation problems such as moisture migration, gas diffusion, local hypoxia and water accumulation extension that are common in underground integrated pipe corridors, this approach is more in line with the actual change process and makes it easier to identify how the anomaly unfolds along the main chain. Attached Figure Description

[0107] Figure 1 This is a schematic diagram illustrating the steps of a method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to the present invention.

[0108] Figure 2 This is a flowchart of the segment boundary determination method for a spatial environment dynamic monitoring and control method for underground integrated utility tunnels according to the present invention.

[0109] Figure 3 This is a flowchart illustrating the over-line state judgment of a spatial environment dynamic monitoring and control method for underground integrated utility tunnels according to the present invention. Detailed Implementation

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

[0111] Example 1

[0112] This invention provides a method for dynamic monitoring and control of the spatial environment of underground utility tunnels. Please refer to [link / reference]. Figures 1 to 3 This includes the following steps:

[0113] S1. Read the structural data of the underground utility tunnel, divide the tunnel space into segments, and construct a segment topology diagram consisting of multiple segment nodes and the connection relationship between adjacent segments.

[0114] S2. Based on the corridor topology diagram, read the existing environmental monitoring data in the corridor and perform unified time axis mapping and corridor-level summary processing to generate the structured monitoring sequence corresponding to each corridor and calculate the comprehensive deviation E(n) corresponding to each corridor.

[0115] S3. Based on the corridor topology diagram and combined with the comprehensive deviation E(n), calculate the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path, and obtain the path-level spatial anomaly value Mp and the anomaly propagation main chain.

[0116] S4. Based on the main chain of abnormal propagation, calculate the control sequence value Kh corresponding to each corridor segment, and sort each corridor segment in the main chain according to the control sequence value to generate the control execution sequence.

[0117] S5. After executing the control sequence, continue to read the environmental monitoring data corresponding to each corridor segment in the anomaly propagation main chain, obtain the structural regression value Yp, and perform feedback control.

[0118] In this embodiment, instead of using the existing method of reading, judging, and alarming separately for each monitoring point, this solution first divides the utility tunnel space into segments and establishes a topological relationship diagram of the segments. Then, it organizes environmental monitoring data around the nodes of the segments. After this processing, the temperature, humidity, oxygen concentration, harmful gas concentration, and water level information, which were originally scattered at different locations, are organized into structured monitoring results consistent with the actual spatial connectivity. This transforms the monitoring object from "single-point numerical changes" to "segmental state changes," which is more in line with the actual situation where the internal environment of the underground integrated utility tunnel is transmitted segment by segment along the connectivity of the compartments, the direction of airflow, and the drainage slope.

[0119] This scheme, based on structured monitoring sequences, continuously analyzes the comprehensive deviation status of each corridor segment and further calculates the transmission sequence, diffusion span, and dwell intensity of candidate propagation paths. It then obtains path-level spatial anomalies and the main propagation chain of anomalies. Based on this process, the system outputs not just a few isolated alarm points, but a spatial anomaly chain that can distinguish the anomaly initiation segment, diffusion segment, and dwell segment. This facilitates the identification of the propagation direction, coverage area, and main dwelling segments of the same environmental anomaly from multiple discrete alarms, solving the problem that existing methods struggle to explain where the anomaly starts, along what path it spreads, and in which segments it persists.

[0120] After completing one round of regulation, this scheme continues to retest each segment in the main chain of anomaly propagation and judges the residual state after regulation by using the structural return value. Through this feedback process, it can identify whether there are still residual over-line segments, residual dwell states, and residual order reversals in the current main chain, and decide whether to rearrange the regulation execution sequence accordingly. Compared with the existing processing logic of "alarm when threshold is reached, end when threshold is released", this scheme not only focuses on whether anomalies occur, but also on whether the anomaly chain truly falls back to the return state after regulation, making it suitable for phased tracking and handling of continuously propagating environmental anomalies.

[0121] This solution is based on existing utility tunnel structure data, existing environmental monitoring point data, and existing actuator layout data. Its technical focus is on spatial organization, sequence rearrangement, path determination, and control sequence generation, without relying on new hardware. For underground utility tunnels that have already been built and put into operation, this solution can be directly adapted to existing monitoring systems and control devices, facilitating the transformation of the original point-based monitoring system into a spatially structured monitoring and control system that considers the propagation relationships within the tunnel segments.

[0122] Example 2

[0123] This embodiment is an explanation based on Embodiment 1. Please refer to it. Figure 2 Specifically: S1 includes S11 and S12;

[0124] S11. Read the corridor centerline, compartment outline, compartment number, connecting door position, ventilation zone boundary, drainage slope change position, environmental monitoring point position and actuator layout data from the pipe gallery structure diagram or model data, and sort each structural position in a chain according to the extension direction of the corridor centerline to form an ordered structural position sequence.

[0125] The boundary criteria are as follows: the continuity of numbering between adjacent structural locations, changes in connecting components, switching of ventilation zones, changes in slope aspect, and mixed affiliation of monitoring points. The boundary value of each adjacent location is calculated using these criteria. The formula for obtaining the boundary value is:

[0126]

[0127] In the formula, Ba represents the corridor boundary value between the a-th position and the (a+1)-th position in the ordered structural position sequence, Na represents the numbering break mark, Ga represents the number of connecting components, Va represents the wind zone switching mark, Pa represents the slope aspect reversal mark, La represents the axial distance between two adjacent structural positions along the centerline direction, and Ma represents the number of monitoring points that are mixed and hung, which is used to characterize the number of monitoring points that are simultaneously hung on different functional objects between two adjacent structural positions.

[0128] The corridor segment boundary value Ba is compared with the corridor segment dividing baseline BX to obtain the corridor segment node set; the specific method is as follows:

[0129] When the segment boundary value Ba is greater than or equal to the segment dividing baseline BX, the current position is determined as the segment boundary point, and the ordered structure sequence is segmented to form multiple segment nodes.

[0130] Then, the starting position, ending position, cabin type, internal monitoring point set, and internal actuator set corresponding to each corridor node are archived to generate a corridor node set.

[0131] S12. Based on the set of corridor nodes, read the doorway connectivity, ventilation direction consistency, drainage slope bearing capacity, and centerline adjacency distance between each corridor node.

[0132] Determine whether there is a direct adjacency relationship between any two corridor segment nodes; for corridor segment nodes with a direct adjacency relationship, further calculate the topological connectivity Cij, obtaining the formula as follows:

[0133]

[0134] In the formula, Cij represents the topological connectivity between corridor node i and corridor node j; Kij represents the interface connection number; Rij represents the wind direction collinearity symbol; Sij represents the slope aspect connection symbol; Dij represents the centerline adjacency distance between two corridor nodes; and Hij represents the barrier layer level.

[0135] After calculating the topological connectivity Cij of all adjacent corridor nodes, the corridor node pairs that meet the direct adjacency condition are written into the adjacency relationship set. The corridor nodes are used as graph nodes, the adjacency relationships are used as graph edges, and the topological connectivity Cij is used as the edge attribute to construct the corridor topological relationship graph.

[0136] In this embodiment, S11 first reads the corridor centerline, cabin outline, cabin number, connecting door position, ventilation zone boundary, drainage slope change position, environmental monitoring point position, and execution device arrangement position in the underground integrated pipe gallery in a unified manner, and forms an ordered structural position sequence according to the extension direction of the corridor centerline. Then, based on the continuous status of the number, the change status of the connecting components, the switching status of the ventilation zone, the change status of the slope, and the status of the monitoring point, it is determined whether the adjacent positions are suitable as the corridor segment boundary positions. This can divide the originally continuously unfolding pipe gallery space into corridor segment nodes with clear start and end ranges and clear internal belonging relationships.

[0137] In this embodiment, the process of dividing the corridor segments considers not only the segment number and spatial location, but also the interconnecting doors, ventilation zones, drainage slopes, and the mixed situation of monitoring point affiliation. This ensures that the determination of the corridor segment boundaries no longer relies solely on a single location marker, but rather on a comprehensive judgment combining the structural connectivity and environmental transmission conditions. For common issues within underground utility tunnels such as moisture migration, gas diffusion, localized oxygen deficiency, and water accumulation spreading along the structural path segment by segment, this division method more closely approximates the actual propagation boundary. It can refine what was originally difficult to distinguish into continuous spaces into multiple analytical units with environmental significance, facilitating the identification of where anomalies originate, where they are segmented, and where they continue to the next segment.

[0138] In this embodiment, step S12 continues to read the doorway connectivity, ventilation direction consistency, drainage slope compatibility, and centerline adjacency distance between each corridor segment node. The topological connectivity of directly adjacent corridor segment nodes is calculated. Then, a corridor segment topological relationship graph is constructed using corridor segment nodes as graph nodes and adjacency relationships as graph edges. This unifies the spatial relationships originally scattered across structural diagrams, ventilation zoning data, and drainage slope data into a single relationship graph. After this processing, the relationships between adjacent corridor segments within the utility tunnel are no longer simply physical adjacencies, but rather segment-level spatial relationships with varying connectivity strengths, wind direction compatibility, and slope compatibility. This facilitates the subsequent understanding of environmental anomalies as a continuous process unfolding along a topological path, rather than a series of unrelated single-point fluctuations.

[0139] In this embodiment, the constructed corridor topology diagram can express in advance which two corridors are directly adjacent, which two corridors are more likely to form environmental transmission, and which segment has a clearer connection with which segment. When reading environmental data such as temperature, humidity, oxygen concentration, harmful gas concentration and water level, the data can be directly reorganized according to the corridor segments and their adjacent relationships.

[0140] In this embodiment, each corridor segment node, upon generation, simultaneously archives its start and end positions, compartment type, internal monitoring point set, and internal actuator set. This establishes a correspondence between the spatial structure, monitoring objects, and control objects within the same corridor segment unit. Consequently, when an anomaly is subsequently identified in a corridor segment, the monitoring source and control target corresponding to that segment can be directly located. This facilitates the control process to revolve around the specific corridor segment, rather than repeatedly matching alarms at multiple points and multiple actuators. For common issues in underground utility tunnel operation and management, such as the disconnect between upstream anomalies and downstream handling, and the unclear correspondence between monitoring results and control targets, this organizational method clarifies the connection between segment-level monitoring and segment-level control.

[0141] In this embodiment, after combining S11 and S12, the segmentation of corridor nodes is completed first, followed by the expression of the adjacency relationships between corridor segments, ultimately forming a corridor topology diagram composed of multiple corridor nodes and the connection relationships between adjacent corridor segments. This process transforms the analysis object of environmental problems within underground integrated utility tunnels from "single-point monitoring values" to "corridor nodes and their connection relationships," which corresponds to the problem described in the background art where existing methods struggle to identify the starting location, propagation direction, and continuous sections of anomalies. Based on this topology diagram, subsequent monitoring, judgment, and control of environmental anomalies can be conducted around the continuity of the spatial structure, making it suitable for segmented identification and sequential handling of continuously propagating environmental problems.

[0142] Example 3

[0143] This embodiment is an explanation based on Embodiment 2. Please refer to it. Figure 1 Specifically: S2 includes S21 and S22;

[0144] S21. Based on the corridor topology diagram, read the temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each monitoring point in each corridor, and simultaneously read the collection start time, collection end time, record holding time and data type identifier corresponding to the original collection records of each monitoring point.

[0145] A unified time-slice sequence is constructed within the current analysis period, dividing the entire analysis period into multiple consecutive time slices, and then the original collection records of each monitoring point are mapped to the corresponding unified time slice;

[0146] After completing the unified time axis mapping of each monitoring point, the mapping values ​​of the same type of monitoring points in the same corridor segment on the same time slice are then summarized at the corridor segment level to form the structured monitoring sequence S={(T(n),U(n),O(n),G(n),W(n))|n=1,2,...,N} corresponding to the corridor segment;

[0147] In the formula, T(n), U(n), O(n), G(n) and W(n) represent the representative values ​​of temperature, humidity, oxygen concentration, harmful gas concentration and water level of the corridor segment in the nth unified time slice, respectively, and N represents the total number of unified time slices obtained in the current analysis period.

[0148] S22. Based on the structured monitoring sequence S, extract samples of the same type of operating conditions that are consistent with the current operating conditions from the historical database, and obtain the temperature reference value, humidity reference value, oxygen concentration reference value, harmful gas concentration reference value and water level reference value of each corridor section under the same type of operating conditions.

[0149] Simultaneously, the historical swing width of various monitoring quantities around their respective benchmark values ​​is obtained; then, the obtained structured monitoring sequence of the corridor segment is compared with the historical benchmark data on a time-by-time basis, and the comprehensive deviation E(n) corresponding to each corridor segment is calculated.

[0150] The overall deviation E(n) is obtained using the following formula:

[0151]

[0152] In the formula, Ei(n) represents the comprehensive deviation of corridor segment i in the nth unified time slice, ln represents the logarithmic function, Sji(n) represents the j-th data in the structured monitoring sequence of corridor segment i in the nth unified time slice, pSji represents the historical baseline value of the j-th data in the structured monitoring sequence of corridor segment i, and bSji represents the historical swing width of the j-th data in the structured monitoring sequence of corridor segment i.

[0153] The comprehensive deviation E(n) is combined with the historical comprehensive deviation sample sequence {Eih(z)} of the corridor segment under the same operating conditions in the historical database to calculate and obtain the dynamic judgment line EiD corresponding to the corridor segment;

[0154] The formula for obtaining the dynamic decision line EiD is: EiD = Med{Eih(z)} + Med{|Eih(z)-Med{Eih(z)}|};

[0155] In the formula, Eih(z) represents the z-th historical comprehensive deviation sample of corridor segment i under the same operating conditions, and Med{Eih(z)} represents the median value of the historical comprehensive deviation sample.

[0156] In this embodiment, in step S21, temperature, humidity, oxygen concentration, harmful gas concentration, and water level data from different monitoring points within the same corridor segment, with different sampling periods and recording durations, are first uniformly organized into the same time-slice sequence. Then, the data are summarized by corridor segment to obtain a structured monitoring sequence that unfolds continuously over time, with the corridor segment as the unit. In this way, environmental records that were originally scattered across various monitoring points are transformed into continuous state results corresponding to the spatial units of the corridor segment, facilitating subsequent assessment of environmental changes based on the corridor segment rather than individual points.

[0157] In this embodiment, S22 does not directly compare the current monitoring value with a fixed threshold. Instead, it first calls up a sample of similar operating conditions that are consistent with the current operating situation to obtain the historical benchmark values ​​and historical swing widths of various monitoring quantities for each corridor segment, and then calculates the comprehensive deviation based on these values. This allows for comparison of the current corridor segment status under the same operating background, which helps to distinguish between normal operating condition fluctuations and abnormal deviations, making the environmental judgment of the corridor segment more consistent with the actual changes in the underground integrated utility tunnel under different compartments and different operating stages.

[0158] In this embodiment, by further organizing the historical comprehensive deviation samples, dynamic judgment lines corresponding to each corridor segment are obtained, so that the judgment boundaries of different corridor segments correspond to their respective historical fluctuations. For corridor segments with relatively stable historical fluctuations, their abnormal changes are easier to identify; for corridor segments with relatively active historical fluctuations, it is necessary to combine their own historical distribution to determine whether the current state is abnormal. This helps to solve the problem that the uniform threshold in the existing method is difficult to adapt to different corridor segments and different operating conditions.

[0159] In this embodiment, after combining S21 and S22, the point data is first organized into a segment-level structured monitoring sequence, and then the sequence is transcribed into a comprehensive deviation and dynamic judgment result, laying the foundation for subsequent identification of abnormal propagation paths and determination of abnormal starting and stopping segments. Compared with the prior art's method of issuing alarms separately for each point, which makes it difficult to see the continuous propagation relationship of abnormalities, this embodiment is more suitable for spatial environment problems such as underground integrated pipe corridors where the connectivity, airflow direction, and slope relationship change segment by segment.

[0160] Example 4

[0161] This embodiment is an explanation based on Embodiment 3. Please refer to it. Figure 3 Specifically: S3 includes S31 and S32;

[0162] S31. For any corridor segment i, the time slice that satisfies Ei(n)≥EiD is taken as the time slice of the corridor segment crossing the line, and the time slice number that first satisfies the condition is recorded as the first time slice ti of the corridor segment crossing the line.

[0163] The total duration covered by the longest consecutively satisfied segment is denoted as the continuous dwell time Qi of the corridor segment; simultaneously, the over-line amplitudes on all over-line time segments of the corridor segment are aggregated to obtain the average over-line amplitude A of the corridor segment, obtained by the following formula:

[0164]

[0165] In the formula, Ai represents the average overshoot amplitude of corridor segment i; SAi represents the set of overshoot time slices of corridor segment i; and Ni+ represents the number of time slices in the set of overshoot time slices.

[0166] After obtaining the first crossing time ti, continuous dwell time Qi, and average crossing amplitude A of each corridor segment, corridor segments with direct adjacent relationships and all having crossing time slices are taken as path connection objects and connected in series along the adjacent connection direction in the corridor segment topology diagram to form several candidate propagation paths p=(h1, h2, ..., hm);

[0167] Where p represents a candidate propagation path; hr represents the r-th segment in path p arranged in connection order; and m represents the number of segments contained in the candidate propagation path.

[0168] For any two adjacent corridor segments hr and hr+1 in the candidate propagation path, first calculate the correction order interval Δtr. When the first crossing time of the latter corridor segment is not earlier than that of the former corridor segment, take the time difference between the two. When a time reversal occurs after the former corridor segment, add an equal amount of penalty to the time difference, specifically obtained by the following formula:

[0169] ;

[0170] Wherein, Δtr specifically represents the correction order interval between the r-th corridor segment and the (r+1)-th corridor segment in the candidate propagation path; t(hr) and t(hr+1) represent the first crossing time of two adjacent corridor segments, respectively;

[0171] Based on the correction order interval Δtr and the average superline amplitude A, the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path are calculated and obtained.

[0172] The method for obtaining the transmission order quantity Op is as follows: First, along the arrangement order of the candidate propagation path, take out two adjacent corridor segments in the path one after the other; for each pair of adjacent corridor segments, read the topological coupling degree between them, and at the same time read the first crossing time of the previous corridor segment and the next corridor segment; then, divide the topological coupling degree between the adjacent corridor segment pairs by "1 and the sum of the successive intervals" to obtain the local contribution value of the adjacent corridor segment pairs to the transmission order quantity of the entire candidate propagation path; then, according to the adjacent order in the path, perform the same processing on all adjacent corridor segment pairs respectively, and accumulate all the obtained local contribution values ​​one by one to obtain the transmission order quantity Op of the candidate propagation path;

[0173] The diffusion span measure Up is obtained as follows: First, according to the order of the candidate propagation paths, the axial length and the corresponding average overshoot of each corridor segment on the path are read sequentially; for each corridor segment, the axial length of the corridor segment is combined with its average overshoot to obtain the single-segment contribution value of the corridor segment to the diffusion range of the entire candidate propagation path; then, the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the diffusion span measure Up of the candidate propagation path.

[0174] The dwell intensity Jp is obtained by sequentially reading the continuous dwell time and the corresponding average overshoot of each corridor segment on the path; for each corridor segment, the continuous dwell time and the average overshoot are combined to obtain the single-segment contribution value of the corridor segment to the dwell state of the entire candidate propagation path; then the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the dwell intensity Jp of the candidate propagation path.

[0175] S32. For each candidate propagation path obtained in step S31, the propagation order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp are jointly calculated to obtain the path-level spatial anomaly value Mp of the candidate propagation path.

[0176] Path-level spatial outliers Mp are obtained in the following way:

[0177] ;

[0178] After obtaining the path-level spatial outliers of all candidate propagation paths, the path-level spatial outliers Mp of each path are compared, and the candidate propagation path with the largest value is taken as the current anomaly propagation main chain.

[0179] Subsequently, within the main chain of abnormal propagation, the initial crossing time, continuous dwell time, and average crossing amplitude of each corridor segment were read, and the dwell contribution value Dhr of each corridor segment was calculated:

[0180] The dwell contribution value Dhr is obtained by multiplying the continuous dwell time of the corridor segment by the average over-line amplitude of the corridor segment;

[0181] Within the main chain, the corridor segment with the earliest first crossing time is designated as the starting corridor segment; the corridor segment with a dwell contribution value Dhr greater than or equal to the median dwell contribution value of all corridor segments in the main chain is designated as the dwell corridor segment; and the remaining main chain corridor segments located after the starting corridor segment and not belonging to the dwell corridor segment are designated as the diffusion corridor segments.

[0182] In this embodiment, based on whether the tunnel segment has entered the over-line state, the first over-line time, continuous dwell time, and average over-line amplitude are extracted. Then, candidate propagation paths are formed along the tunnel segment topology map. This can organize the anomaly records that were originally scattered across multiple monitoring points and multiple time slices into a path result that unfolds continuously according to spatial connection relationships. In this way, environmental anomalies no longer appear as several isolated alarm points, but rather a continuous propagation chain can be identified along the actual connection direction of the tunnel, which corresponds to the problem in the background technology that single-point judgment is difficult to reflect the anomaly propagation process.

[0183] In this embodiment, by constructing transmission sequence, diffusion span, and dwell intensity quantities respectively, the temporal sequence, spatial extension range, and abnormal dwell state of candidate propagation paths are jointly judged. This allows for the differentiation between two different situations: "multiple corridor segments fluctuating simultaneously" and "anomalies propagating segment by segment along adjacent corridor segments." For continuous propagation problems commonly encountered in underground utility tunnels, such as moisture migration, gas diffusion, localized oxygen deficiency, and water accumulation, this approach better reflects the actual change process and facilitates the identification of how anomalies unfold along the main chain.

[0184] In this embodiment, after comparing each candidate propagation path, the path with the largest spatial outlier is selected as the current main propagation chain of the anomaly. This allows for the selection of the main anomaly chain that best matches the current anomaly state from multiple possible paths. Thus, subsequent analysis no longer focuses on "which point exceeds the line," but can further clarify "which segment first exhibits an anomaly, which segments belong to the diffusion zone, and which segments belong to the persistent dwelling zone," thereby providing a clearer spatial representation of the anomaly's starting location, propagation direction, and main dwelling location.

[0185] Example 5

[0186] Please refer to Figure 1 Specifically: S4 includes S41 and S42;

[0187] S41. According to the arrangement order in the main chain of abnormal propagation, extract the continuous dwell time, average over-line amplitude, topological coupling degree with the upstream adjacent corridor segment, the shortest span to the actionable device, and the segment-level distance of the corridor segment relative to the starting corridor segment in sequence.

[0188] Among them, the upstream coupling relationship of the starting corridor section is handled according to the main chain entrance mark;

[0189] For any corridor segment h in the main chain, first determine its upstream adjacent corridor segment uh, and then calculate the control sequence value Kh of the corridor segment.

[0190] The control sequence value Kh is obtained as follows: the continuous dwell time of the corridor segment, the average overshoot amplitude, and the topological coupling degree with the upstream corridor segment are taken as the sequential lifting term. The shortest span number from the corridor segment to the existing actuator and the segment-level distance relative to the starting corridor segment are taken as the sequential pressure drop term. The result of the product of the sequential lifting terms is divided by the sum of the sequential pressure drop term and 1 to obtain the control sequence value Kh corresponding to the corridor segment.

[0191] S42. Based on the control sequence value Kh, sort all corridor segments in the main chain from largest to smallest to obtain the control sequence within the main chain.

[0192] When two or more corridor segments have the same control priority value, the corridor segment with the smaller segment distance is selected first, and then the corridor segment with the larger average over-line amplitude is selected as the target of the advance control.

[0193] After completing the corridor segment sorting, the environmental monitoring items with the largest deviation contribution in each corridor segment, as well as the action type, action area and cross-segment action boundary of each existing actuator, are read. For each sorted corridor segment, an actuator or control program that can act on the corridor segment and is consistent with the dominant deviation type of the corridor segment is selected, and the regulation matching value between the corridor segment and the candidate actuator is calculated.

[0194] The method for calculating the control matching value is as follows: take the consistent state of the actuator and the current corridor segment's dominant deviation type, and the continuous coverage state of the actuator's action section on the current corridor segment and its downstream adjacent corridor segments as the matching lifting term, take the number of segments spanned from the actuator to the current corridor segment as the matching pressure drop term, and then divide the result of the multiplication of the matching lifting term by the sum of the matching pressure drop term and 1 to obtain the control matching value between the current corridor segment and the corresponding actuator.

[0195] For each sorted corridor segment, the control program with the largest control matching value is selected, and the "corridor segment number - execution object - execution order" are sequentially connected to form a control execution sequence.

[0196] S5 includes S51 and S52; S51: After the control execution sequence is completed, select the retest time window covering the observation period after control to obtain the temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each corridor segment in the abnormal propagation main chain.

[0197] Following the unified time axis mapping method and corridor segment-level aggregation method in step S2, the corridor segment representative temperature, corridor segment representative humidity, corridor segment representative oxygen concentration, corridor segment representative harmful gas concentration and corridor segment representative water level within the retest time window after regulation are obtained, and the comprehensive deviation of each corridor segment after regulation on each unified time slice is recalculated.

[0198] For any segment in the main chain of abnormal propagation, the time slice that still meets the condition of "comprehensive deviation not lower than the dynamic judgment line of the segment" after regulation is taken as the residual over-line time slice of the segment; the total duration corresponding to the longest continuous segment in the residual over-line time slice is recorded as the remaining dwell time Qhr of the segment; the over-line portion on the residual over-line time slice is accumulated and converted according to the number of residual over-line time slices to obtain the remaining over-line amplitude Ahr of the segment;

[0199] After obtaining the remaining dwell time Qhr and the remaining overshoot amplitude Ahr of each corridor segment in the main chain, the order of the remaining time within the main chain after regulation is re-examined according to the original arrangement of each corridor segment in the main chain.

[0200] If the residual first overrun time of the downstream corridor segment is earlier than the residual first overrun time of the upstream corridor segment, then the adjacent corridor segments are recorded as a residual back jump segment.

[0201] After checking all adjacent corridor pairs, the number of residual bounce segments Bp in the current abnormal propagation main chain is obtained; at the same time, all corridor segments in the main chain are checked one by one to see if there are still residual time slices, and the number of corridor segments with residual time slices is recorded as the number of residual time slices Sp.

[0202] By combining the continuous dwell time Qi and average over-line amplitude A of each corridor segment in the main chain before regulation, and the remaining dwell time Qhr, remaining over-line amplitude Ahr, remaining bounce segment number Bp and remaining over-line corridor segment number Sp obtained after regulation, the structural regression value Yp corresponding to the current abnormal propagation main chain is obtained.

[0203]

[0204] In the formula, m represents the number of segments contained in the main chain of abnormal propagation;

[0205] S52. First, extract completed control events from the historical event database that are identical to the current abnormal propagation main chain in terms of cabin type, dominant deviation type, main chain length, and calling execution device type, and read the structural regression value sample corresponding to the historical event at the end of the control; then, obtain the regression decision line YpD of the current abnormal propagation main chain based on the structural regression value sample.

[0206] The regression decision line YpD is obtained as follows: First, sort the structural regression value samples of similar historical regulation events and extract their median values; then calculate the dispersion of each sample relative to the median value; finally, subtract the dispersion from the median value to obtain the regression decision line YpD corresponding to the current abnormal propagation main chain.

[0207] Compare the structural regression value Yp with the regression decision line YpD to determine the current state of the anomaly propagation main chain.

[0208] The judgment method is as follows:

[0209] When the structural regression value Yp ≥ the regression decision line YpD, the current abnormal propagation main chain is determined to enter the regression state, the current control result is output, and the control end time, regression corridor set, residual corridor set and final structural regression value corresponding to the current main chain are written into the historical event database.

[0210] When the structural regression value Yp < regression decision line YpD, it is determined that the current abnormal propagation main chain has not entered the regression state. Further, the corridor segments that still have residual overshoot time slices in the current main chain are read, and the remaining dwell time Qhr, remaining overshoot amplitude Ahr, residual bounce relationship and residual segment position of the corridor segments are sent back to step S4 as inputs for the next round of regulation order calculation and regulation execution sequence rearrangement.

[0211] In this embodiment, S4 first calculates the control priority value of each corridor segment based on the continuous dwell time, average overshoot amplitude, topological connection relationship, cross-segment relationship of the execution device, and segment-level positional relationship of each corridor segment within the main chain of anomaly propagation. Then, the corridor segments in the main chain are sorted accordingly. This ensures that the control order no longer depends on the order of occurrence of individual alarm points, but rather on the position of each segment within the main chain within the anomaly chain. Compared with the parallel handling of multiple alarm points in the prior art, this method better reflects the actual situation of environmental anomalies propagating segment by segment within the underground integrated utility tunnel. It facilitates handling the initial segment and the continuously dwelling segment first, followed by the subsequent connecting segments, ensuring that the control order corresponds to the spatial unfolding state of the anomaly chain.

[0212] In this embodiment, after the corridor segments are sorted, the same control method is not directly applied to all abnormal corridor segments. Instead, the environmental monitoring items that contribute the most to the deviation of each corridor segment are read, and the control mapping relationship between the corridor segment and the candidate control devices is calculated by combining the function type, function area, and cross-segment boundary of the existing actuators. Then, the execution object that is more suitable for the current dominant deviation state is selected for each corridor segment. In this way, the fan, drainage pump, access control linkage program, and ventilation zone switching program are no longer called in a fixed mode, but are linked according to the dominant environmental problems of different corridor segments. For common situations in underground integrated pipe corridors such as gas migration, local hypoxia, moisture diffusion, and water accumulation, this processing method can match the execution object with the abnormal attributes and main chain connection relationship of the current corridor segment, which is convenient for forming a control execution sequence that unfolds segment by segment.

[0213] In this embodiment, after the control execution sequence ends, S5 continues to conduct retests around each corridor segment within the anomaly propagation main chain. This involves not only rereading various environmental monitoring data but also recalculating the remaining dwell state, remaining out-of-line state, residual bounce relationship, and residual out-of-line corridor range after control. Thus, the system's judgment no longer stops at "whether control has been performed" or "whether a certain point has returned to within the threshold," but further determines whether the original anomaly main chain still retains residual propagation relationships after control. Compared to the background technology's processing logic of ending the process upon clearing the alarm at a single point threshold, this embodiment is more suitable for identifying the residual state of continuously propagating anomalies after control.

[0214] In this embodiment, the main chain state before and after regulation is compared by using structural regression values. Then, regression results from similar historical regulation events are used as the basis for judgment. This allows for the determination of whether the current main chain has entered a regression state within the historical context of similar compartment types, similar dominant deviation types, similar main chain lengths, and similar actuator types. Thus, regression judgment no longer relies on a single fixed boundary but corresponds to the current anomaly type and regulation scenario. For tunnel environments with varying degrees of fluctuation under different compartments and operating conditions, this approach makes it easier to distinguish between the states of "the main chain has fallen back" and "the main chain still has residual anomalies."

[0215] Although embodiments of the invention 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 and spirit of the invention, the scope of which is defined by the appended claims and their equivalents.

Claims

1. A method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels, characterized in that: Includes the following steps: S1. Read the structural data of the underground utility tunnel, divide the tunnel space into segments, and construct a segment topology diagram consisting of multiple segment nodes and the connection relationship between adjacent segments. S2. Based on the corridor topology diagram, read the existing environmental monitoring data in the corridor and perform unified time axis mapping and corridor-level summary processing to generate the structured monitoring sequence corresponding to each corridor and calculate the comprehensive deviation E(n) corresponding to each corridor. S3. Based on the corridor topology diagram and combined with the comprehensive deviation E(n), calculate the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path, and obtain the path-level spatial anomaly value Mp and the anomaly propagation main chain. S4. Based on the main chain of abnormal propagation, calculate the control sequence value Kh corresponding to each corridor segment, and sort each corridor segment in the main chain according to the control sequence value to generate the control execution sequence. S5. After executing the control sequence, continue to read the environmental monitoring data corresponding to each corridor segment in the anomaly propagation main chain, obtain the structural regression value Yp, and perform feedback control. 2.The method of claim 1, wherein the method comprises: S1 includes S11 and S12; S11. Read the corridor centerline, compartment outline, compartment number, connecting door position, ventilation zone boundary, drainage slope change position, environmental monitoring point position and actuator layout data from the pipe gallery structure diagram or model data, and sort each structural position in a chain according to the extension direction of the corridor centerline to form an ordered structural position sequence. The boundary criteria are as follows: the continuity of numbering between adjacent structural locations, changes in connecting components, switching of ventilation zones, changes in slope aspect, and mixed affiliation of monitoring points. The boundary value of each adjacent location is calculated using these criteria. The formula for obtaining the boundary value is: In the formula, Ba represents the corridor boundary value between the a-th position and the (a+1)-th position in the ordered structural position sequence, Na represents the numbering break mark, Ga represents the number of connecting components, Va represents the wind zone switching mark, Pa represents the slope aspect reversal mark, La represents the axial distance between two adjacent structural positions along the centerline direction, and Ma represents the number of monitoring points that are mixed and hung, which is used to characterize the number of monitoring points that are simultaneously hung on different functional objects between two adjacent structural positions. The corridor segment boundary value Ba is compared with the corridor segment dividing baseline BX to obtain the corridor segment node set; the specific method is as follows: When the segment boundary value Ba is greater than or equal to the segment dividing baseline BX, the current position is determined as the segment boundary point, and the ordered structure sequence is segmented to form multiple segment nodes. Then, the starting position, ending position, cabin type, internal monitoring point set, and internal actuator set corresponding to each corridor node are archived to generate a corridor node set.

3. The method of claim 2, wherein the method comprises: S12. Based on the set of corridor nodes, read the doorway connectivity, ventilation direction consistency, drainage slope bearing capacity, and centerline adjacency distance between each corridor node. Determine whether there is a direct adjacency relationship between any two corridor segment nodes; for corridor segment nodes with a direct adjacency relationship, further calculate the topological connectivity Cij, obtaining the formula as follows: In the formula, Cij represents the topological connectivity between corridor node i and corridor node j; Kij represents the interface connection number; Rij represents the wind direction collinearity symbol; Sij represents the slope aspect connection symbol; and Dij represents the centerline adjacency distance between two corridor nodes. Hij represents the number of barrier layers; After calculating the topological connectivity Cij of all adjacent corridor nodes, the corridor node pairs that meet the direct adjacency condition are written into the adjacency relationship set. The corridor nodes are used as graph nodes, the adjacency relationships are used as graph edges, and the topological connectivity Cij is used as the edge attribute to construct the corridor topological relationship graph.

4. The method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 3, characterized in that: S2 includes S21 and S22; S21. Based on the corridor topology diagram, read the temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each monitoring point in each corridor, and simultaneously read the collection start time, collection end time, record holding time and data type identifier corresponding to the original collection records of each monitoring point. A unified time-slice sequence is constructed within the current analysis period, dividing the entire analysis period into multiple consecutive time slices, and then the original collection records of each monitoring point are mapped to the corresponding unified time slice; After completing the unified time axis mapping of each monitoring point, the mapping values ​​of the same type of monitoring points in the same corridor segment on the same time slice are then summarized at the corridor segment level to form the structured monitoring sequence S={(T(n),U(n),O(n),G(n),W(n))|n=1,2,...,N} corresponding to the corridor segment; In the formula, T(n), U(n), O(n), G(n) and W(n) represent the representative values ​​of temperature, humidity, oxygen concentration, harmful gas concentration and water level of the corridor segment in the nth unified time slice, respectively, and N represents the total number of unified time slices obtained in the current analysis period.

5. The method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 4, characterized in that: S22. Based on the structured monitoring sequence S, extract samples of the same type of operating conditions that are consistent with the current operating conditions from the historical database, and obtain the temperature reference value, humidity reference value, oxygen concentration reference value, harmful gas concentration reference value and water level reference value of each corridor section under the same type of operating conditions. Simultaneously, the historical swing width of various monitoring quantities around their respective benchmark values ​​is obtained; then, the obtained structured monitoring sequence of the corridor segment is compared with the historical benchmark data on a time-by-time basis, and the comprehensive deviation E(n) corresponding to each corridor segment is calculated. The overall deviation E(n) is obtained using the following formula: In the formula, Ei(n) represents the comprehensive deviation of corridor segment i in the nth unified time slice, ln represents the logarithmic function, Sji(n) represents the j-th data in the structured monitoring sequence of corridor segment i in the nth unified time slice, pSji represents the historical baseline value of the j-th data in the structured monitoring sequence of corridor segment i, and bSji represents the historical swing width of the j-th data in the structured monitoring sequence of corridor segment i. The comprehensive deviation E(n) is combined with the historical comprehensive deviation sample sequence {Eih(z)} of the corridor segment under the same operating conditions in the historical database to calculate and obtain the dynamic judgment line EiD corresponding to the corridor segment; The formula for obtaining the dynamic decision line EiD is: EiD = Med{Eih(z)} + Med{|Eih(z)-Med{Eih(z)}|}; In the formula, Eih(z) represents the z-th historical comprehensive deviation sample of corridor segment i under the same operating conditions, and Med{Eih(z)} represents the median value of the historical comprehensive deviation sample.

6. The method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 5, characterized in that: S3 includes S31 and S32; S31. For any corridor segment i, the time slice that satisfies Ei(n)≥EiD is taken as the time slice of the corridor segment crossing the line, and the time slice number that first satisfies the condition is recorded as the first time slice ti of the corridor segment crossing the line. The total duration covered by the longest consecutively satisfied segment is denoted as the continuous dwell time Qi of the corridor segment; simultaneously, the over-line amplitudes on all over-line time segments of the corridor segment are aggregated to obtain the average over-line amplitude A of the corridor segment, obtained by the following formula: In the formula, Ai represents the average overshoot amplitude of corridor segment i; SAi represents the set of overshoot time slices of corridor segment i; and Ni+ represents the number of time slices in the set of overshoot time slices. After obtaining the first crossing time ti, continuous dwell time Qi, and average crossing amplitude A of each corridor segment, corridor segments with direct adjacent relationships and all having crossing time slices are taken as path connection objects and connected in series along the adjacent connection direction in the corridor segment topology diagram to form several candidate propagation paths p=(h1, h2, ..., hm); Where p represents a candidate propagation path; hr represents the r-th segment in path p arranged in connection order; and m represents the number of segments contained in the candidate propagation path. For any two adjacent corridor segments hr and hr+1 in the candidate propagation path, first calculate the correction order interval Δtr. When the first crossing time of the latter corridor segment is not earlier than that of the former corridor segment, take the time difference between the two. When a time reversal occurs after the former corridor segment, add an equal amount of penalty to the time difference, specifically obtained by the following formula: ； Wherein, Δtr specifically represents the correction order interval between the r-th corridor segment and the (r+1)-th corridor segment in the candidate propagation path; t(hr) and t(hr+1) represent the first crossing time of two adjacent corridor segments, respectively; Based on the correction order interval Δtr and the average superline amplitude A, the transmission order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp corresponding to the candidate propagation path are calculated and obtained. The method for obtaining the transmission order quantity Op is as follows: First, along the arrangement order of the candidate propagation path, take out two adjacent corridor segments in the path one after the other; for each pair of adjacent corridor segments, read the topological coupling degree between them, and at the same time read the first crossing time of the previous corridor segment and the next corridor segment; then, divide the topological coupling degree between the adjacent corridor segment pairs by "1 and the sum of the successive intervals" to obtain the local contribution value of the adjacent corridor segment pairs to the transmission order quantity of the entire candidate propagation path; then, according to the adjacent order in the path, perform the same processing on all adjacent corridor segment pairs respectively, and accumulate all the obtained local contribution values ​​one by one to obtain the transmission order quantity Op of the candidate propagation path; The diffusion span measure Up is obtained as follows: First, according to the order of the candidate propagation paths, the axial length and the corresponding average overshoot of each corridor segment on the path are read sequentially; for each corridor segment, the axial length of the corridor segment is combined with its average overshoot to obtain the single-segment contribution value of the corridor segment to the diffusion range of the entire candidate propagation path; then, the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the diffusion span measure Up of the candidate propagation path. The dwell intensity Jp is obtained by sequentially reading the continuous dwell time and the corresponding average overshoot of each corridor segment on the path; for each corridor segment, the continuous dwell time and the average overshoot are combined to obtain the single-segment contribution value of the corridor segment to the dwell state of the entire candidate propagation path; then the single-segment contribution values ​​of all corridor segments in the path are accumulated sequentially to obtain the dwell intensity Jp of the candidate propagation path.

7. The method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 6, characterized in that: S32. For each candidate propagation path obtained in step S31, the propagation order quantity Op, diffusion span quantity Up, and dwell intensity quantity Jp are jointly calculated to obtain the path-level spatial anomaly value Mp of the candidate propagation path. Path-level spatial outliers Mp are obtained in the following way: ； After obtaining the path-level spatial outliers of all candidate propagation paths, the path-level spatial outliers Mp of each path are compared, and the candidate propagation path with the largest value is taken as the current anomaly propagation main chain. Subsequently, within the main chain of abnormal propagation, the initial crossing time, continuous dwell time, and average crossing amplitude of each corridor segment were read, and the dwell contribution value Dhr of each corridor segment was calculated: The dwell contribution value Dhr is obtained by multiplying the continuous dwell time of the corridor segment by the average over-line amplitude of the corridor segment; Within the main chain, the corridor segment with the earliest first crossing time is designated as the starting corridor segment; the corridor segment with a dwell contribution value Dhr greater than or equal to the median dwell contribution value of all corridor segments in the main chain is designated as the dwell corridor segment; and the remaining main chain corridor segments located after the starting corridor segment and not belonging to the dwell corridor segment are designated as the diffusion corridor segments.

8. A method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 7, characterized in that: S4 includes S41 and S42; S41. According to the arrangement order in the main chain of abnormal propagation, extract the continuous dwell time, average over-line amplitude, topological coupling degree with the upstream adjacent corridor segment, the shortest span to the actionable device, and the segment-level distance of the corridor segment relative to the starting corridor segment in sequence. Among them, the upstream coupling relationship of the starting corridor section is handled according to the main chain entrance mark; For any corridor segment h in the main chain, first determine its upstream adjacent corridor segment uh, and then calculate the control sequence value Kh of the corridor segment. The control sequence value Kh is obtained as follows: the continuous dwell time of the corridor segment, the average overshoot amplitude, and the topological coupling degree with the upstream corridor segment are taken as the sequential lifting term. The shortest span number from the corridor segment to the existing actuator and the segment-level distance relative to the starting corridor segment are taken as the sequential pressure drop term. The result of the product of the sequential lifting terms is divided by the sum of the sequential pressure drop term and 1 to obtain the control sequence value Kh corresponding to the corridor segment.

9. A method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 8, characterized in that: S42. Based on the control sequence value Kh, sort all corridor segments in the main chain from largest to smallest to obtain the control sequence within the main chain. When two or more corridor segments have the same control priority value, the corridor segment with the smaller segment distance is selected first, and then the corridor segment with the larger average over-line amplitude is selected as the target of the advance control. After completing the corridor segment sorting, the environmental monitoring items with the largest deviation contribution in each corridor segment, as well as the action type, action area and cross-segment action boundary of each existing actuator, are read. For each sorted corridor segment, an actuator or control program that can act on the corridor segment and is consistent with the dominant deviation type of the corridor segment is selected, and the regulation matching value between the corridor segment and the candidate actuator is calculated. The method for calculating the control matching value is as follows: take the consistent state of the actuator and the current corridor segment's dominant deviation type, and the continuous coverage state of the actuator's action section on the current corridor segment and its downstream adjacent corridor segments as the matching lifting term, take the number of segments spanned from the actuator to the current corridor segment as the matching pressure drop term, and then divide the result of the multiplication of the matching lifting term by the sum of the matching pressure drop term and 1 to obtain the control matching value between the current corridor segment and the corresponding actuator. For each sorted corridor segment, the control program with the largest control matching value is selected, and the "corridor segment number - execution object - execution order" are sequentially connected to form a control execution sequence.

10. A method for dynamic monitoring and control of the spatial environment of underground integrated utility tunnels according to claim 9, characterized in that: S5 includes S51 and S52; S51. After the control execution sequence is completed, select a retest time window covering the observation period after control, and obtain temperature data, humidity data, oxygen concentration data, harmful gas concentration data and water level data corresponding to each corridor segment in the abnormal propagation main chain. Following the unified time axis mapping method and corridor segment-level aggregation method in step S2, the corridor segment representative temperature, corridor segment representative humidity, corridor segment representative oxygen concentration, corridor segment representative harmful gas concentration and corridor segment representative water level within the retest time window after regulation are obtained, and the comprehensive deviation of each corridor segment after regulation on each unified time slice is recalculated. For any segment in the main chain of abnormal propagation, the time slice that still meets the condition of "comprehensive deviation not lower than the dynamic judgment line of the segment" after regulation is taken as the residual over-line time slice of the segment. The total duration corresponding to the longest continuous segment in the residual overline time slice is recorded as the remaining dwell time Qhr of the corridor segment; the overline portion on the residual overline time slice is accumulated and converted according to the number of residual overline time slices to obtain the remaining overline amplitude Ahr of the corridor segment; After obtaining the remaining dwell time Qhr and the remaining overshoot amplitude Ahr of each corridor segment in the main chain, the order of the remaining time within the main chain after regulation is re-examined according to the original arrangement of each corridor segment in the main chain. If the residual first overrun time of the downstream corridor segment is earlier than the residual first overrun time of the upstream corridor segment, then the adjacent corridor segments are recorded as a residual back jump segment. After checking all adjacent corridor pairs, the number of residual bounce segments Bp in the current abnormal propagation main chain is obtained; at the same time, all corridor segments in the main chain are checked one by one to see if there are still residual time slices, and the number of corridor segments with residual time slices is recorded as the number of residual time slices Sp. By combining the continuous dwell time Qi and average over-line amplitude A of each corridor segment in the main chain before regulation, and the remaining dwell time Qhr, remaining over-line amplitude Ahr, remaining bounce segment number Bp and remaining over-line corridor segment number Sp obtained after regulation, the structural regression value Yp corresponding to the current abnormal propagation main chain is obtained. In the formula, m represents the number of segments contained in the main chain of abnormal propagation; S52. First, extract completed control events from the historical event database that are identical to the current abnormal propagation main chain in terms of cabin type, dominant deviation type, main chain length, and calling execution device type, and read the structural regression value sample corresponding to the historical event at the end of the control; then, obtain the regression decision line YpD of the current abnormal propagation main chain based on the structural regression value sample. The regression decision line YpD is obtained as follows: First, sort the structural regression value samples of similar historical regulation events and extract their median values; then calculate the dispersion of each sample relative to the median value; finally, subtract the dispersion from the median value to obtain the regression decision line YpD corresponding to the current abnormal propagation main chain. Compare the structural regression value Yp with the regression decision line YpD to determine the current state of the anomaly propagation main chain. The judgment method is as follows: When the structural regression value Yp ≥ the regression decision line YpD, the current abnormal propagation main chain is determined to enter the regression state, the current control result is output, and the control end time, regression corridor set, residual corridor set and final structural regression value corresponding to the current main chain are written into the historical event database. When the structural regression value Yp < regression decision line YpD, it is determined that the current abnormal propagation main chain has not entered the regression state. Further, the corridor segments that still have residual overshoot time slices in the current main chain are read, and the remaining dwell time Qhr, remaining overshoot amplitude Ahr, residual bounce relationship and residual segment position of the corridor segments are sent back to step S4 as inputs for the next round of regulation order calculation and regulation execution sequence rearrangement.

Images