A civil air defense engineering opening site air tightness detection method, device and medium

Abstract

This invention discloses a method, equipment, and medium for on-site airtightness testing of the entrance of civil defense engineering projects, specifically relating to the field of airtightness testing. It addresses the problem in existing technologies that can only assess overall airtightness and struggle to identify specific leakage channels and propagation paths. By releasing charged tracer smoke into the entrance space and deploying electrostatic induction arrays on the inner and outer sides, the induced current signal generated during the leakage process is collected. The signal arrival time difference and amplitude characteristics are analyzed to extract the migration characteristics of different leakage channels. Combined with a three-dimensional geometric model of the entrance structure, the leakage propagation path is reconstructed, thereby marking the main leakage channel and the detour path. This achieves refined identification of the leakage location and path in complex structural entrances, improving the spatial representation capability and engineering application accuracy of the test results. Simultaneously, by comprehensively analyzing the migration speed and concentration decay characteristics of different channels, a basis is provided for subsequent sealing defect location and reinforcement treatment.

Description

Technical Field

[0001] This invention relates to the field of airtightness testing technology, and more specifically, to a method, equipment, and medium for on-site airtightness testing of the entrance of civil defense engineering. Background Technology

[0002] In civil defense projects, the entrances and exits are critical areas for personnel, equipment, and pipelines. These areas are typically equipped with protective doors, airtight doors, and various through-wall pipeline structures. This area is also the most sensitive and prone to leakage in terms of the overall airtightness of the project. During long-term use, factors such as aging of door seals, installation deviations in door frames, shrinkage and cracking of sealing materials around pipeline penetrations, and repeated opening and closing or vibration can easily create minute but continuous leakage paths at the entrances. When external shock waves or toxic gases are present, these leakage paths will directly affect the internal protective effectiveness of the project.

[0003] Current methods for detecting the airtightness of the inlet mostly employ the overall pressure difference method or the ordinary smoke tracing method. These methods determine the presence of a leak by observing the pressure decay or smoke escape. However, these methods often only provide an overall assessment of the sealing performance or a rough estimate of the leak location. In some scenarios where the leak is not visible to the naked eye, it is difficult to pinpoint the specific leak path and channel distribution under complex structural conditions. This is especially true for inlet areas with dense wall-penetrating pipelines and complex structural connections, where smoke diffuses and overlaps in multiple paths, making it impossible to distinguish between the main leak channel and secondary detours. This results in difficulties in subsequent maintenance and location, leading to low efficiency.

[0004] Meanwhile, traditional detection methods lack the ability to analyze the spatiotemporal characteristics of the leakage propagation process and cannot quantitatively characterize the flow velocity and attenuation characteristics of different channels, making it difficult to support the needs of refined operation and maintenance and targeted reinforcement. Therefore, it is necessary to propose an on-site airtightness detection method that can identify and reconstruct the leakage propagation path in complex opening structures. Summary of the Invention

[0005] In order to overcome the above-mentioned defects of the prior art, the embodiments of the present invention provide a method, equipment and medium for on-site airtightness testing of the entrance of civil defense engineering to solve the problems mentioned in the background art.

[0006] To achieve the above objectives, the present invention provides the following technical solution:

[0007] A method for on-site airtightness testing of the entrance of a civil defense project includes the following steps:

[0008] S1. Seal the entrance of the civil defense project and apply an electrostatic field to the tracer in the smoke generator to drive the smoke to carry a net charge dominated by a single polarity.

[0009] S2. At the inlet of the pipe connecting the inside and outside of the civil defense project, start the electrostatic induction array to monitor the induced current inside the inlet in real time.

[0010] S3. Gas is pumped into the mouth through the reserved tube, causing the internal air pressure to rise to the preset test threshold, triggering the charged smoke to flow in a direction to the potential leak under the pressure difference.

[0011] S4. When charged smoke escapes from the leak channel corresponding to the area where the pipe wall hole seal fails to be sealed, it drives the electrostatic induction array on the outside of the opening to monitor the induced current.

[0012] S5. Analyze the arrival time difference and intensity ratio of the signals captured by each electrode in the electrostatic induction array, and extract the migration velocity and concentration decay distribution of charged particles in different leakage channels within the limited electric field influence range.

[0013] S6. Based on the migration velocity and concentration decay distribution, combined with the mouth structure constraints, the trajectory of charged particles in three-dimensional space is inverted, and the main leakage channel and detour path are marked on the mouth structure model.

[0014] As a further aspect of the present invention, in S1, driving the smoke to carry a net charge dominated by a single polarity specifically includes:

[0015] Based on the volume of the internal space of the mouth to be tested and the preset test overpressure value, the required charge intensity range of the smoke particles is set according to the smoke diffusion coverage requirements. The output voltage of the electrostatic generator is adjusted to the initial set value according to the charge intensity range, so that the tracer aerosol particles flowing through the electrostatic ring capture charge in the corona discharge area.

[0016] Smoke samples were collected at the exit of the electrostatic ring, and the charge-to-mass ratio of the samples was analyzed. By measuring the deflection displacement of the particles in the electric field, the actual charge-to-mass ratio distribution range of the smoke particles was calculated based on the deflection displacement.

[0017] The actual charge-to-mass ratio distribution range is compared with the preset target range. If the median of the actual distribution range deviates from the target range, the inter-electrode voltage of the electrostatic ring or the discharge current of the corona electrode is adjusted step by step according to the preset step size until the actual charge-to-mass ratio distribution range falls into the target range. The adjusted voltage and current values ​​are then locked.

[0018] As a further aspect of the present invention, in step S2, activating the electrostatic induction array to monitor the induced current inside the opening in real time specifically includes:

[0019] The background noise current signal of each electrostatic induction array in a smoke-free state is obtained, and the noise baseline value of the corresponding electrode of the electrostatic induction array is recorded. The baseline value is used as the subtraction reference for the induced current signal processing.

[0020] Based on the structural gap characteristics and expected leakage intensity of each monitoring point located on the inner side of the entrance of the civil defense project, an induced current trigger level threshold is set for each electrostatic induction array.

[0021] The electrostatic induction array is activated to collect the current signal output from each monitoring point inside the opening. If the amplitude of the current signal exceeds the trigger level threshold, it is determined that charged smoke has passed through the monitoring point. The waveform segment of the induced current at the corresponding moment is recorded, and the peak amplitude of the waveform segment is extracted as the characteristic parameter of the smoke passage event.

[0022] As a further aspect of the present invention, in step S3, triggering the directional flow of charged smoke towards the potential leak under pressure differential specifically includes:

[0023] Gas is continuously injected into the mouth through a pre-reserved tube, and the real-time pressure value output by the air pressure measuring device installed inside the mouth is read at a fixed sampling frequency. The real-time pressure value is then compared with a preset test threshold.

[0024] When the real-time pressure value approaches the preset test threshold and the difference between the two enters the set approximation range, the injection mode of the pumping device is switched from continuous injection to intermittent pulse injection.

[0025] When the real-time pressure value reaches the preset test threshold, the pumping stops, the electric valve of the reserved pipe is closed to lock the internal pressure at the preset test threshold, and the charged smoke is driven to flow directionally along the leakage channel.

[0026] As a further aspect of the present invention, in S4, the monitoring of the induced current by the electrostatic induction array on the outer side of the driving port specifically includes:

[0027] The current signal output by the electrostatic induction array corresponding to each monitoring point on the outside of the mouth is collected in real time, and the current amplitude at each sampling moment is compared with the corresponding channel trigger level threshold.

[0028] When the current amplitude exceeds the trigger level threshold, the moment is recorded as the time when the smoke arrives. At the same time, the current waveform data within the set time window before and after the corresponding moment is extracted, and the peak amplitude of the waveform data is calculated as the signal characteristic parameter of the event.

[0029] As a further aspect of the present invention, step S5, extracting the migration velocity and concentration decay distribution of charged particles in different leakage channels, specifically includes:

[0030] From each recorded smoke passage event, the arrival time of the signal generated on the sensing electrodes of the electrostatic induction array at different spatial locations of the same leakage channel is extracted, the signal arrival time difference between the inner and outer electrodes is calculated, and the peak amplitude of the signal corresponding to the leakage event is read to calculate the ratio of the peak amplitude between the inner and outer electrodes.

[0031] Based on the actual layout coordinates of each electrode in the space inside and outside the mouth, the signal arrival time difference is correlated with the distance between the inner and outer electrodes to calculate the average migration speed of charged smoke in this section of the leakage channel. At the same time, based on the attenuation relationship between the peak amplitude ratio and the electrode distance, the concentration attenuation coefficient of charged smoke on the path from the leakage source to the electrode is inferred.

[0032] The average migration velocity and concentration decay coefficient calculated for each leakage channel are input into the leakage channel parameter inversion model constructed based on the preset leakage channel flow relationship. The migration velocity, concentration decay coefficient, and inner and outer electrode spatial coordinates are used as inputs, and the equivalent aperture and geometric orientation of the leakage channel are used as outputs, which are used as the flow characteristic parameters of the leakage channel. The characteristic parameters of other leakage channels are calculated simultaneously to construct the leakage channel characteristic parameter set.

[0033] As a further aspect of the present invention, in step S6, marking the main leakage channel and the detour path specifically includes:

[0034] Read the characteristic parameters of each leakage channel from the characteristic parameter set of the leakage channel, and establish a three-dimensional geometric model of the opening structure. The three-dimensional geometric model includes the spatial position and geometric constraint boundary of the door, wall, pipeline through-wall hole and adjacent connecting pipeline.

[0035] The spatial coordinates of each electrode are mapped to the three-dimensional geometric model of the mouth structure. The electrode coordinate points are used as the starting points of the trajectory. The direction of movement of charged particles in the channel is determined according to the geometric direction parameters of each leakage channel. At the same time, the equivalent aperture is used as the channel cross-sectional size constraint to generate flow path segments that meet the direction and cross-sectional constraints.

[0036] The generated flow path segments are spliced ​​together according to spatial connection relationship to form a continuous three-dimensional trajectory line covering the entire leakage propagation path. At the same time, the bifurcation points and intersection points in the trajectory line are identified. Multiple trajectory lines after the bifurcation point are marked as potential detour paths, and the mainstream trajectory line before the intersection point is marked as the main leakage channel.

[0037] The generated 3D trajectory line, main leakage channel, and detour path are superimposed in different colors on the 3D geometric model of the mouth structure, and the output is a 3D leakage trajectory map including path annotations.

[0038] The present invention also includes a computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the on-site airtightness testing method for the entrance of a civil defense project as described in any of the above claims.

[0039] The present invention also includes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the on-site airtightness testing method for the entrance of a civil defense project as described in any of the above claims.

[0040] The technical effects and advantages of the present invention regarding the on-site airtightness testing method, equipment, and medium for the entrance of civil defense projects are as follows:

[0041] By introducing charged tracer smoke into the air intake space and simultaneously detecting it using inner and outer electrostatic induction arrays, dynamic response acquisition of the leakage process is achieved. Compared to traditional methods that rely solely on pressure difference or visible smoke, this approach not only determines the presence of a leak but also acquires the temporal characteristics and spatial distribution information of the leak. Furthermore, by analyzing the arrival time difference and amplitude characteristics of the induction signals, differences in migration behavior in different leakage channels can be extracted. This is further combined with a three-dimensional geometric model of the air intake structure to reconstruct the leakage propagation path, thereby distinguishing between the main leakage channel and the detour path, elevating the detection results from a single judgment to a path-level expression. This method is suitable for air defense air intake environments with dense through-wall pipelines and complex structural connections. It can identify concealed leakage channels without damaging the structure, providing a clear location basis for subsequent sealing repairs. Simultaneously, by quantifying the propagation characteristics of different channels, it supports the development of differentiated reinforcement measures, improving the overall precision of airtightness detection and its field application value. Attached Figure Description

[0042] Figure 1 This is a schematic diagram of a method for on-site airtightness testing of the entrance of a civil defense project according to the present invention. Detailed Implementation

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

[0044] Example 1

[0045] Figure 1 This invention provides a method for on-site airtightness testing of the entrance of civil defense engineering, which includes the following steps:

[0046] S1. Seal the entrance of the civil defense project and apply an electrostatic field to the tracer in the smoke generator to drive the smoke to carry a net charge dominated by a single polarity.

[0047] S2. At the inlet of the pipe connecting the inside and outside of the civil defense project, start the electrostatic induction array to monitor the induced current inside the inlet in real time.

[0048] S3. Gas is pumped into the mouth through the reserved tube, causing the internal air pressure to rise to the preset test threshold, triggering the charged smoke to flow in a direction to the potential leak under the pressure difference.

[0049] S4. When charged smoke escapes from the leak channel corresponding to the area where the pipe wall hole seal fails to be sealed, it drives the electrostatic induction array on the outside of the opening to monitor the induced current.

[0050] S5. Analyze the arrival time difference and intensity ratio of the signals captured by each electrode in the electrostatic induction array, and extract the migration velocity and concentration decay distribution of charged particles in different leakage channels within the limited electric field influence range.

[0051] S6. Based on the migration velocity and concentration decay distribution, combined with the mouth structure constraints, the trajectory of charged particles in three-dimensional space is inverted, and the main leakage channel and detour path are marked on the mouth structure model.

[0052] In S1, the process of driving the smoke to carry a net charge dominated by a single polarity specifically includes:

[0053] By calculating the volume of the internal space of the test opening, the effective gas capacity under closed conditions is obtained. For example, the internal volume of a typical air-raid shelter opening is approximately 30 to 80 cubic meters. Then, the test overpressure value is determined based on the testing requirements. This overpressure value is set using a fixed configuration, for example, selecting 250 Pa as the standard test pressure to ensure that this value can form a stable pressure difference without causing additional damage to the structure. Based on this, and considering the diffusion characteristics of smoke within this volume, the concentration level at which smoke particles need to form a continuous coverage per unit time is determined, and the required charge intensity range of the particles is deduced accordingly. Specifically, the detection sensitivity of the electrostatic induction array is used as a reference; when the number of charged particles per cubic centimeter reaches 10... 4 Up to 10 5 At this level, the sensing signal can be stably triggered, thereby setting the charge intensity range of a single particle to 10. -15 Up to 10 -13 The Coulomb range is defined, where each unit in the array corresponds to a sensing electrode. After setting the charge intensity range, the output voltage of the electrostatic generator is adjusted to the corresponding initial value, for example, setting the corona electrode voltage within the range of 6kV to 10kV, and fixing the distance between the electrostatic ring electrodes at the preset structural dimensions. This allows aerosol particles to continuously undergo ionization and acquire charge as they pass through the corona region. During this process, the smoke generator continuously and stably outputs tracer aerosols with a particle size distribution concentrated in the range of 0.5 micrometers to 2 micrometers, giving the particles good charge acceptance and spatial dispersion when passing through the corona region, thereby forming a charged smoke stream dominated by unipolarity.

[0054] A sampling channel is set up at the exit of the electrostatic ring to introduce a portion of the charged smoke into the detection area. This detection area consists of a stable electric field environment formed by a pair of parallel electrode plates with a fixed spacing, for example, 5 cm. A constant voltage, for example, 2 kV, is applied to create a uniform electric field. A particle trajectory acquisition device is set up at the exit of the sampling channel to record the trajectory of smoke particles within the electric field region using high-speed imaging. The initial direction of motion of the particles before entering the electric field is used as a reference baseline, and the lateral offset distance of the particles under the influence of the electric field is recorded. Since the time it takes for a particle to pass through the electric field region is determined by the airflow velocity, which is kept stable by a fixed flow control device, for example, set to 0.5 m / s, the force state of the particles in the electric field is consistent. Through statistical analysis of the overall trajectory of multiple particles, the deflection displacement of the particles is statistically distributed and divided into multiple intervals based on the magnitude of the deflection displacement, for example, segmented statistically at intervals of 0.1 mm, thus forming a deflection displacement distribution curve. Furthermore, based on the electric field strength and the direction of force on the particles, a correspondence is established between the deflection displacement and the amount of charge carried by the particles. The statistically obtained deflection displacement interval is mapped to the charge-to-mass ratio interval, thereby obtaining the actual charge-to-mass ratio distribution interval of the current smoke particles. Through continuous sampling and statistics, this distribution interval is made stable and representative, avoiding the fluctuations caused by single sampling.

[0055] A preset target charge-to-mass ratio range is established, set according to the detection sensitivity and spatial diffusion requirements of the electrostatic induction array. For example, the target range can be set as a symmetrical range centered on a certain standard value, specifically, the median value corresponding to a deflection displacement of 0.3 mm, with an allowable fluctuation range of ±0.1 mm. The obtained actual charge-to-mass ratio distribution range is then compared with this target range, and the median of the actual distribution range is extracted as the judgment criterion. When the median deviates from the target range, the parameter adjustment process is initiated. The adjustment method adopts graded step control, adjusting the corona electrode voltage in fixed increments or decrements, such as 0.5 kV each time, while simultaneously fine-tuning the corona discharge current, such as 0.1 mA each time, to ensure stable changes in ionization intensity. After each adjustment, the smoke sampling and deflection displacement statistics process is repeated to obtain the updated charge-to-mass ratio distribution range, and the comparison is performed again. When the median values ​​of two consecutive sampling results fall within the target range, the current voltage and current parameters are considered to have reached a stable state, and this set of parameters is used as the locked parameters for this detection and remains unchanged. During the locking process, the output stability of the electrostatic generator is maintained.

[0056] In step S2, the electrostatic induction array is activated to monitor the induced current inside the opening in real time.

[0057] In an initial state with no smoke input, background noise current is collected from the electrostatic induction array deployed at each monitoring point inside the opening. These monitoring points include the inner and outer entrances / exits of each through-wall pipeline and the entrances / exits of connecting pipelines between adjacent protective units. The local area is manually sealed off to ensure isolation from external airflow disturbances and charged particle interference. Then, the induction electrodes at each monitoring point are sequentially energized and initialized to a stable operating state. Subsequently, the output current signal of each electrode is collected within a continuous fixed time window, for example, continuously for more than 60 seconds, maintaining a stable sampling frequency at a uniform level, such as 1000 samples per second, to ensure sufficient capture of weak current fluctuations. During the collection process, no tracer smoke is introduced to ensure that the recorded signals only reflect the background electrical noise of the environment, the thermal noise of the electrodes themselves, and weak spatial electric field disturbances. Statistical processing is performed on the current sequences acquired by each electrode to remove occasional spike interference signals. For example, isolated peaks that significantly outweigh the overall data distribution are filtered out. Then, the average value of the remaining data within a stable range is calculated, and the natural fluctuation range above and below this average value is recorded. For example, a stable fluctuation band is determined by the distribution range of the statistical data. This average value is determined as the noise baseline value corresponding to that electrode, and the fluctuation band is used as a reference range for subsequent judgments. For different monitoring points, due to differences in their spatial structure (e.g., near door gaps, pipe interfaces, or wall corners), the electric field environment varies. Therefore, the noise baseline values ​​of each electrode are recorded independently, and a uniform value is not used. After completing the baseline acquisition for all monitoring points, the baseline values ​​of each electrode are numbered and bound to form a correspondence. In subsequent detection processes, this baseline value is used as the subtraction benchmark for the induced current signal. That is, the baseline value of the corresponding electrode is subtracted from the real-time acquired current signal, so that subsequent analysis only processes the current change caused by charged smoke.

[0058] After determining the noise baseline, the trigger level of the induced current is set for each electrostatic induction array based on the structural gap characteristics and expected leakage intensity at each monitoring point location. First, on-site measurements are taken of the internal structure of the opening to obtain the gap width, length, and morphological characteristics near each monitoring point. For example, the contact width under compression is recorded for the door sealing strip area, and the gap size between the pipe diameter and the sleeve is recorded for pipe penetrations through walls. Based on the actual conditions of different structural locations, the monitoring points are divided into several types, such as door gaps, holes, and structural joints. Different leakage response levels are set for different types, and the trigger level setting method is determined accordingly. Specifically, for tightly sealed structural areas with gap widths within 1 mm, the trigger level is set to approximately twice the natural fluctuation range of the baseline value; for structural areas with gap widths between 1 mm and 3 mm, the trigger level is set to approximately three times the natural fluctuation range of the baseline value; and for areas with obvious structural gaps or aged seals, the trigger level is set to approximately four times the natural fluctuation range of the baseline value. This setup method is calibrated through on-site pre-experimentation. A small amount of charged smoke is introduced into a limited area, and the response at each monitoring point is observed. The multiplier is then fine-tuned based on the actual detected signal amplitude until a stable detection of real smoke signals is achieved without false triggering. After setup, the trigger level for each electrode is fixed, and its correspondence with the monitoring point location and structural type is recorded. This ensures that each channel has a targeted trigger sensitivity, thereby avoiding misjudgments caused by environmental noise and guaranteeing effective response to weak leakage signals during subsequent detection.

[0059] After setting the baseline value and trigger level, the electrostatic induction array is activated to continuously acquire current signals from each monitoring point inside the opening. A uniform sampling frequency is maintained during acquisition, and time information is recorded synchronously for each channel to ensure comparability of signals between different monitoring points. During real-time acquisition, baseline subtraction is performed on each current signal, subtracting the baseline value of the corresponding electrode from the current sampled value, causing the signal to fluctuate around zero. When a current signal exceeds the trigger level corresponding to that channel within a series of consecutive sampling points, an event of charged smoke passing through that monitoring point is determined, and this moment is immediately recorded as the event start time. Simultaneously, a set time window is extended forward and backward from this moment, for example, recording 0.2 seconds forward and 0.5 seconds backward, extracting complete waveform segments from the current data within this time period. Feature extraction is performed on the extracted waveform segments. First, the maximum current value in the segment is identified as the peak amplitude, and its occurrence time and position are recorded. Simultaneously, the overall waveform shape is smoothed to remove local minor fluctuations, ensuring the stability of peak identification. For cases where smoke passes through the same monitoring point consecutively within a short period, the time interval between adjacent events is determined. Multiple events with intervals less than a set merging time are grouped into a single continuous smoke passage event, and the overall peak amplitude of the merged waveform is re-extracted. Finally, the peak amplitude and occurrence time of each smoke passage event are recorded as characteristic parameters of the monitoring point and linked to its location.

[0060] In step S3, the charged smoke is triggered to flow directionally toward the potential leak under the drive of pressure difference.

[0061] After the opening is sealed and a stable space is established, gas is continuously injected into the interior through a pre-installed pipe to gradually increase the internal pressure. The gas injection device is sealed to the pre-installed pipe, and after the gas source is turned on, gas is delivered into the opening at a constant flow rate. This flow rate is fixedly set by a mechanical flow regulating valve, for example, within the range of 0.3 to 0.6 cubic meters per minute, to ensure a smooth and continuous pressure rise. During the gas injection process, a pressure measuring device deployed inside the opening is simultaneously activated. This pressure measuring device uses a micro-pressure sensor with a range covering 0 to 500 Pa and is positioned near the center of the structure inside the opening to avoid local flow causing deviations in the measurement results. The sampling frequency is set using a fixed configuration, for example, 10 samples per second, to ensure continuous recording of the pressure change curve. The preset test threshold is determined according to the airtightness testing specifications for civil defense projects and is fixedly set in conjunction with the structural bearing capacity, for example, 250 Pa is selected as the standard test threshold. During continuous gas injection, the real-time collected pressure value is compared with the threshold value one by one, and the change in the pressure rise rate is recorded simultaneously. When the pressure value is in the initial rising stage, the pressure increase shows a linear trend. As it approaches the target threshold, the growth rate gradually slows down. A complete pressure change curve is formed through continuous sampling, providing a basis for judgment for subsequent injection mode switching.

[0062] As the real-time pressure gradually approaches the preset test threshold, the injection method is controlled in stages. When the difference between the real-time pressure value and the preset test threshold enters a pre-set approximation range, an injection mode switching operation is performed. This approximation range is set using a fixed interval method, for example, 10% of the target threshold. That is, when the target threshold is 250 Pa, the range from 225 Pa to 250 Pa is used as the approximation range. After the pressure value enters this range, the gas injection method is immediately switched from continuous injection to intermittent pulse injection. By controlling the opening and closing cycle of the inlet valve, gas is intermittently introduced into the inlet in a short-term opening and closing manner, for example, a cycle mode with an opening time of 2 seconds and a closing time of 3 seconds is set. During each injection cycle, the real-time output value of the pressure measuring device is continuously read, and the natural pressure drop during the closing phase is observed. Through multiple cycles of repetition, the pressure gradually converges as it approaches the target value. This process limits the amount of gas injected at one time to prevent the pressure from exceeding the target threshold, and simultaneously uses multiple small-amplitude increments to stabilize the pressure and bring it close to the target range. During intermittent injection, the pressure variation amplitude within each cycle is recorded. When the pressure peaks of several consecutive cycles gradually approach the target threshold and the fluctuation range converges to a smaller interval, the current pressure control process is determined to have entered a stable approximation state. When the real-time pressure value reaches the preset test threshold, gas injection termination and pressure locking operations are immediately executed. In specific implementation, the gas supply channel of the gas injection device is closed the moment the pressure value is detected to first reach or slightly exceed the target threshold. At the same time, the electric valve on the reserved pipe is controlled to close rapidly, forming a closed space and maintaining the internal pressure of the inlet at the current level. After the valve is closed, the output value of the pressure measuring device continues to be read at a fixed sampling frequency to continuously monitor the pressure stability. For example, pressure changes are recorded over the next 30 to 60 seconds. When the pressure drop remains within the preset allowable range, such as not exceeding 5 Pa fluctuation, the pressure locking is determined to be successful. Under this stable pressure condition, a constant pressure difference is formed between the inside and outside of the inlet. This pressure difference drives the charged smoke to flow directionally along the existing leakage channel, which is the leakage channel formed by the sealing failure area around the wall penetration hole. Because the pressure is stable, the direction of smoke flow remains consistent, avoiding changes in flow direction caused by pressure fluctuations, thus ensuring that the flow behavior in each leak channel is repeatable and identifiable.

[0063] In S4, the electrostatic induction array on the outer side of the drive port monitors the induced current.

[0064] Electrodes at each outer monitoring point are independently connected according to their spatial location, ensuring that each signal has an independent channel and synchronously records time information. The sampling frequency is set using a fixed configuration, for example, uniformly set to 10 samples per second, to ensure complete capture of transient current changes. During acquisition, baseline subtraction is performed on each current signal, that is, the noise baseline value recorded by the corresponding electrode in a smoke-free state is subtracted from the real-time sample value, causing the signal to fluctuate around the zero value range, highlighting the effective changes caused by charged smoke. Subsequently, the current amplitude at each sampling moment is compared point-by-point with the pre-set trigger level of that channel. The trigger level is set according to the aforementioned method based on structural location differences; for example, for monitoring points located immediately adjacent to wall penetrations, the trigger level is set to three times the natural fluctuation range above the baseline value; for monitoring points located far from major structural gaps, the trigger level is set to four times the natural fluctuation range above the baseline value. During continuous sampling, when the current amplitude of a channel consistently exceeds its corresponding trigger level at multiple consecutive sampling points, it is determined that the channel has detected charged smoke, and this moment is immediately recorded as the smoke arrival time. Meanwhile, to avoid false triggering by single-point noise, a continuity constraint is added to the trigger determination. For example, it is required that at least 5 consecutive sampling points exceed the trigger level to be considered a valid event.

[0065] After confirming a smoke arrival event at a monitoring point, the corresponding current waveform is fully extracted and its features are calculated. Specifically, a time window is set forward and backward from the recorded smoke arrival time, extracting current waveform data for the corresponding time period from the original continuous sampling data. The time window is set in a fixed manner, for example, extending forward by 0.2 seconds and backward by 0.5 seconds, to ensure coverage of the baseline fluctuation phase before smoke arrival and the complete response process after arrival. After extracting the waveform, this data segment is processed uniformly. First, the waveform is smoothed using a moving average method, for example, averaging the local data using a 5-point moving average window to eliminate high-frequency noise interference and make the overall waveform trend clearer. Then, the maximum current value is searched in the smoothed waveform data, and this maximum value is determined as the peak amplitude of this smoke passage event, while the time position of this peak is recorded. In the case of multiple local peaks, the amplitudes of each peak are compared, and the one with the largest amplitude is selected as the final feature value, while other smaller fluctuations are ignored. For situations where multiple triggers occur consecutively at the same monitoring point within a short period, the time interval between adjacent events is determined. Events with intervals shorter than a set merging time are merged, for example, a merging time of 0.3 seconds. Multiple consecutive triggering events are merged into a single complete event, and the peak amplitude is re-extracted within the merged waveform range. The smoke arrival time and corresponding peak amplitude of each monitoring point are bound and recorded to form the event feature data for that monitoring point.

[0066] In step S5, the migration velocity and concentration decay distribution of charged particles in different leakage channels are extracted.

[0067] After collecting smoke passage events from the electrostatic induction arrays inside and outside the mouth, the recorded event data are organized and matched to extract response information from electrodes at different spatial locations along the same leakage channel. Based on the spatial relationships and relative structural connections of each monitoring point, the inner and outer electrodes are paired to form several candidate channel correspondences. For example, electrodes located on both sides of the same through-wall hole are defined as an initial matching pair, and adjacency matching relationships are established for electrodes in adjacent structural gap areas. Within each matching pair, event sequences that are temporally continuous and satisfy propagation logic are selected from the recorded smoke passage events. Specifically, the response time of the outer electrode must be later than that of the inner electrode, and the time interval between them must be within a reasonable range, such as 0.05 seconds to 2 seconds. For event pairs that meet the criteria, the corresponding smoke arrival time is read, and the time difference between the inner and outer electrodes is calculated. Simultaneously, the peak amplitudes corresponding to the events on both sides are read, and the ratio of the outer peak amplitude to the inner peak amplitude is calculated. To avoid the influence of occasional fluctuations on the results, multiple time difference and amplitude ratio data generated during multiple smoke passages through the same channel were centrally statistically processed. A fixed interval distribution method was used to screen stable interval data; for example, data within the middle 50% of the overall distribution were selected as valid data, and representative values ​​were calculated as the characteristic time difference and amplitude ratio of that channel. Through this processing, the response relationship of the same leakage channel at different spatial locations was uniformly expressed, thus forming stable time difference and amplitude ratio characteristics.

[0068] After obtaining the time difference and amplitude ratio characteristics of each leakage channel, the propagation characteristics of smoke within the leakage channel are quantitatively calculated by combining the actual coordinates of each electrode in the space inside and outside the opening. First, the actual coordinates of each electrode in three-dimensional space are read, and the shortest connection path length between the inner and outer electrodes is calculated based on the opening structure model. This path length is limited by structural constraints, meaning it is only measured along the passable structural area and does not cross solid walls. This path length is correlated with the corresponding time difference characteristic, and the average migration velocity within that segment of the leakage channel is obtained by dividing the path length by the time difference. This average migration velocity uses an equivalent characterization method, reflecting the overall propagation efficiency of charged smoke in the channel. Regarding concentration decay calculation, the peak amplitude ratio of the inner and outer electrodes is correlated with the spatial distance between the two electrodes. Based on the decreasing amplitude with distance, the amplitude ratios in different distance intervals are statistically segmented, forming an amplitude decay curve. After smoothing this curve, the overall trend is extracted, and the degree of amplitude decay with distance is converted into a concentration decay coefficient. The specific setting method is as follows: the percentage decrease in amplitude when the distance doubles is used as a reference indicator. For example, if the amplitude decreases by approximately 40% when the distance increases from 0.5 meters to 1 meter, the attenuation coefficient corresponding to that channel is defined within the corresponding level range. By comprehensively statistically analyzing multiple sets of data, the attenuation coefficient is made stable and serves as a unified parameter describing the concentration change characteristics of smoke during its propagation within the channel.

[0069] After obtaining the average migration velocity and concentration decay coefficient of each leakage channel, a parameter inversion model based on the flow relationship of the leakage channels is constructed to infer the equivalent pore size and geometric orientation of the leakage channels. In specific implementation, a leakage channel characteristic sample library is established based on existing engineering structures and historical detection data. This sample library contains flow behavior data under different pore sizes, different leakage channel lengths, and different structural morphologies. Each sample group records the corresponding migration velocity and decay characteristics. Subsequently, the sample library is hierarchically organized, and the leakage channels are divided into multiple intervals according to pore size, such as small pore size interval (less than 1 mm), medium pore size interval (1 mm to 3 mm), and large pore size interval (greater than 3 mm). Within each interval, the corresponding migration velocity range and decay coefficient range are statistically analyzed to form a standard comparison table. Based on this, an inversion model is constructed, consisting of input mapping and parameter matching. Input mapping receives the actual detected migration velocity, attenuation coefficient, and electrode spatial coordinates, converting them into a standardized feature sequence. Parameter matching compares the input features with features in each interval according to a standard lookup table, determining the corresponding aperture interval based on the closest match principle. Furthermore, it combines the relative positions of the electrode spatial coordinates to infer the geometric orientation of the leakage channel, such as determining whether the leakage channel is straight-through, bend-shaped, or laterally extended. During model operation, fixed matching rules are used for judgment. For example, when the migration velocity is in the medium range and the attenuation coefficient is low, the leakage channel is classified as a medium-aperture straight-through structure; when the migration velocity is low and the attenuation coefficient is high, it is determined to be a narrow and bend-shaped leakage channel structure. By performing the above inversion process on each leakage channel, the corresponding equivalent aperture and geometric orientation parameters are obtained and recorded as the flow characteristic parameters of that leakage channel. Finally, the characteristic parameters of all leakage channels are summarized to form a complete set of leakage channel characteristic parameters.

[0070] In S6, the main leakage channel and the detour path are marked.

[0071] Based on engineering design drawings and on-site measurement data, the spatial dimensions and relative positions of doors, walls, pipe penetrations, and adjacent connecting pipes were obtained. Each structural unit was expressed using a unified coordinate system, establishing a three-dimensional connectivity topology between the internal and external spaces of the opening, forming a three-dimensional geometric model containing geometric constraints and spatial connectivity. During modeling, solid structures such as walls and doors were defined as impassable areas, and their geometric boundaries were clearly defined. Areas such as wall penetrations, pipe interiors, and structural gaps were defined as potentially passable areas, and their spatial connectivity was recorded. For pipe penetrations, their axial direction and diameter range were further marked, and the connection between the penetrations and the internal and external spaces was clearly identified. After completing the basic structural modeling, the characteristic parameters of each leakage channel were associated and bound to the structural model, mapping the equivalent aperture, geometric orientation parameters, and corresponding internal and external electrode positions of each leakage channel to their corresponding spatial positions in the three-dimensional geometric model. This method ensures that the structural model includes not only geometric information but also channel feature information corresponding to the detection data, thus forming a three-dimensional geometric model that integrates physical constraints and detection results. In this model, the positional relationships between each structural boundary and channel are expressed in explicit spatial coordinates, guaranteeing that the subsequent path generation process proceeds within the established geometric constraints.

[0072] After constructing the 3D geometric model, the spatial coordinates of each electrode are mapped into the model, and flow path segments are generated based on the geometric orientation parameters of the leakage channel and the equivalent aperture constraint. Specifically, the spatial coordinate points of the inner and outer electrodes are used as the starting and ending constraint positions for path generation, and the corresponding positions are marked in the model. Then, based on the geometric orientation information given in the leakage channel characteristic parameters, a dominant direction is determined in the region where the electrodes are located; that is, a preferred extension direction is defined in 3D space, such as along the axial direction of the through-wall opening or along the structural gap. Guided by this dominant direction, the path gradually extends from the electrode coordinate points in the 3D geometric model. Simultaneously, constraints are applied to the path cross-section during the extension process; that is, the path is only allowed to pass through spatial regions with a cross-sectional size larger than the equivalent aperture, and is not allowed to pass through regions smaller than this size. During path generation, the space is divided into continuous small-scale units, and connectivity is judged between adjacent units. Only when there is a geometric connection between two units and the cross-sectional size requirement is the path extended from the current unit to the next unit. By gradually expanding the path, it continues to extend outward while satisfying directional and dimensional constraints until it reaches the corresponding outer electrode position or enters another known connected region. In cases with multiple possible extension directions, the path with the smaller angle to the dominant direction is preferentially selected based on geometric parameters, thus ensuring the continuity and directional consistency of the generated path segments. Finally, a set of flow path segments conforming to structural constraints is generated between each pair of inner and outer electrodes, and these are recorded as candidate leakage paths.

[0073] After obtaining the candidate leakage path segments, they are spatially stitched and node identified to form a complete three-dimensional leakage propagation trajectory. Based on the spatial endpoint positions of the path segments, it is determined whether there is a connection between different segments. When the spatial distance between the endpoints of two segments is less than a preset connection threshold, they are considered connectable. This connection threshold is set using a fixed distance method, for example, 5 mm, to accommodate the structural modeling error range. Under the condition of meeting the connection criteria, the corresponding segments are connected end-to-end to gradually form a continuous path. During the stitching process, the continuity of the path is checked to ensure that the path always lies within a passable area and does not cross any physical structures. After the path stitching is completed, node analysis is performed on the formed continuous trajectory. Specifically, all connecting nodes in the trajectory are traversed. When a node connects two or more subsequent path segments simultaneously, the node is marked as a branch point; when a node is formed by two or more path segments merging and continuing to extend forward, the node is marked as a junction point. For a branch point, its subsequent extended paths are marked as different candidate branch paths; for a junction point, its preceding path is defined as the main merging path. Based on this, the characteristic parameters of each path are compared. A comprehensive judgment is made based on the peak amplitude of the signal, the propagation time, and the concentration attenuation degree corresponding to each path. The peak amplitude is sorted from largest to smallest, the arrival time from smallest to largest, and the concentration attenuation coefficient from smallest to largest. The combined ranking position of each trajectory in these three rankings is calculated, and the trajectory ranking at the top in most rankings is selected as the main leakage channel and marked as the main leakage channel. The remaining paths extending from the bifurcation points are marked as detour paths. After completing the path classification, all trajectory lines are visualized according to category, distinguished by different colors in the 3D geometric model. For example, the main leakage channel is displayed in red, and detour paths are displayed in blue or green, while maintaining the spatial correspondence between each path and the structural model. Finally, a 3D leakage trajectory map containing complete path annotation information is output, intuitively presenting the spatial distribution relationship of each leakage path.

[0074] The present invention also includes a computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the on-site airtightness testing method for the entrance of a civil defense project as described in any of the above claims.

[0075] The present invention also includes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the on-site airtightness testing method for the entrance of a civil defense project as described in any of the above claims.

[0076] The above embodiments can be implemented, in whole or in part, by software, hardware, firmware, or any other combination thereof. When implemented using software, the above embodiments can be implemented, in whole or in part, as a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, all or part of the processes or functions described in the embodiments of this application are generated. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more sets of available media. The available medium can be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium. The semiconductor medium can be a solid-state drive.

[0077] Those skilled in the art will recognize that the modules and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

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

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

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

[0081] In addition, the functional modules in the various embodiments of this application can be integrated into one processing module, or each module can exist physically separately, or two or more modules can be integrated into one module.

[0082] If the aforementioned functions are implemented as software functional modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0083] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

[0084] In conclusion, the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for detecting the air tightness of a site of a civil air defense engineering entrance, characterized in that, Includes the following steps: S1. Seal the entrance of the civil defense project, apply an electrostatic field to the tracer in the smoke generator, and drive the smoke to carry a net charge dominated by a single polarity, specifically including: Based on the volume of the internal space of the mouth to be tested and the preset test overpressure value, the required charge intensity range of the smoke particles is set according to the smoke diffusion coverage requirements. The output voltage of the electrostatic generator is adjusted to the initial set value according to the charge intensity range, so that the tracer aerosol particles flowing through the electrostatic ring capture charge in the corona discharge area. Smoke samples were collected at the exit of the electrostatic ring, and the charge-to-mass ratio of the samples was analyzed. By measuring the deflection displacement of the particles in the electric field, the actual charge-to-mass ratio distribution range of the smoke particles was calculated based on the deflection displacement. The actual charge-to-mass ratio distribution range is compared with the preset target range. If the median of the actual distribution range deviates from the target range, the inter-electrode voltage of the electrostatic ring or the discharge current of the corona electrode is adjusted step by step according to the preset step size until the actual charge-to-mass ratio distribution range falls into the target range. The adjusted voltage and current values ​​are then locked. S2. At the inlet of the pipe connecting the inside and outside of the civil defense project, start the electrostatic induction array to monitor the induced current inside the inlet in real time. S3. Gas is pumped into the mouth through the reserved tube, causing the internal air pressure to rise to the preset test threshold, triggering the charged smoke to flow in a direction to the potential leak under the pressure difference. S4. When charged smoke escapes from the leak channel corresponding to the area where the pipe wall hole seal fails to be sealed, it drives the electrostatic induction array on the outside of the opening to monitor the induced current. S5. Analyze the arrival time difference and intensity ratio of the signals captured by each electrode in the electrostatic induction array, and extract the migration velocity and concentration decay distribution of charged particles in different leakage channels, specifically including: From the recorded smoke passage events, the arrival time of the signals generated on the induction electrodes of the electrostatic induction array at different spatial locations of the same leakage channel is extracted, the signal arrival time difference between the inner and outer electrodes is calculated, and the peak amplitude corresponding to the current waveform data is read to calculate the ratio of the peak amplitude between the inner and outer electrodes as the intensity ratio. Based on the actual layout coordinates of each electrode in the space inside and outside the mouth, the signal arrival time difference is correlated with the distance between the inner and outer electrodes to calculate the average migration speed of charged smoke in the leakage channel. At the same time, based on the attenuation relationship between the peak amplitude ratio and the electrode distance, the concentration attenuation coefficient of charged smoke on the path from the leakage source to the electrode is inferred. The average migration velocity and concentration decay coefficient calculated for each leakage channel are input into the leakage channel parameter inversion model constructed based on the preset leakage channel flow relationship. The migration velocity, concentration decay coefficient and inner and outer electrode spatial coordinates are used as inputs, and the equivalent aperture and geometric orientation of the leakage channel are used as outputs, which are used as the flow characteristic parameters of the leakage channel. The characteristic parameters of other leakage channels are calculated simultaneously to construct the leakage channel characteristic parameter set. S6. Based on the migration velocity and concentration decay distribution, combined with the mouth structure constraints, the trajectory of charged particles in three-dimensional space is inverted, and the main leakage channel and detour path are marked on the mouth structure model.

2. The method for on-site airtightness testing of the entrance of a civil defense project according to claim 1, characterized in that, In step S2, activating the electrostatic induction array to monitor the induced current inside the opening in real time specifically includes: The background noise current signal of each electrostatic induction array in a smoke-free state is obtained, and the noise baseline value of the corresponding electrode of the electrostatic induction array is recorded. The baseline value is used as the subtraction reference for the induced current signal processing. Based on the structural gap characteristics and expected leakage intensity of each monitoring point located on the inner side of the entrance of the civil defense project, an induced current trigger level threshold is set for each electrostatic induction array. The electrostatic induction array is activated to collect the current signal output from each monitoring point inside the opening. If the amplitude of the current signal exceeds the trigger level threshold, it is determined that charged smoke is passing through the monitoring point.

3. The method for on-site airtightness testing of the entrance of a civil defense project according to claim 1, characterized in that, In step S3, triggering the directional flow of charged smoke towards the potential leak under pressure differential specifically includes: Gas is continuously injected into the mouth through a pre-reserved tube, and the real-time pressure value output by the air pressure measuring device installed inside the mouth is read at a fixed sampling frequency. The real-time pressure value is then compared with a preset test threshold. When the real-time pressure value approaches the preset test threshold and the difference between the two enters the set approximation range, the injection mode of the pumping device is switched from continuous injection to intermittent pulse injection. When the real-time pressure value reaches the preset test threshold, the pumping stops, the electric valve of the reserved pipe is closed to lock the internal pressure at the preset test threshold, and the charged smoke is driven to flow directionally along the leakage channel.

4. The method for on-site airtightness testing of the entrance of a civil defense project according to claim 1, characterized in that, In step S4, the monitoring of the induced current by the electrostatic induction array on the outer side of the driving port specifically includes: The current signal output by the electrostatic induction array corresponding to each monitoring point on the outside of the mouth is collected in real time, and the current amplitude at each sampling moment is compared with the corresponding channel trigger level threshold. When the current amplitude exceeds the trigger level threshold, record that moment as the time when the smoke arrives, and at the same time extract the current waveform data within the set time window before and after the corresponding moment, and calculate the peak amplitude of the waveform data.

5. The method for on-site airtightness testing of the entrance of a civil defense project according to claim 1, characterized in that, In step S6, the marking of the main leakage channel and the detour path specifically includes: Read the characteristic parameters of each leakage channel from the characteristic parameter set of the leakage channel, and establish a three-dimensional geometric model of the opening structure. The three-dimensional geometric model includes the spatial position and geometric constraint boundary of the door, wall, pipeline through-wall hole and adjacent connecting pipeline. The spatial coordinates of each electrode are mapped to the three-dimensional geometric model of the mouth structure. The electrode coordinate points are used as the starting points of the trajectory. The direction of movement of charged particles in the channel is determined according to the geometric direction parameters of each leakage channel. At the same time, the equivalent aperture is used as the channel cross-sectional size constraint to generate flow path segments that meet the direction and cross-sectional constraints. The generated flow path segments are spliced ​​together according to spatial connection relationship to form a continuous three-dimensional trajectory line covering the entire leakage propagation path. At the same time, the bifurcation points and intersection points in the trajectory line are identified. Multiple trajectory lines after the bifurcation point are marked as potential detour paths, and the mainstream trajectory line before the intersection point is marked as the main leakage channel. The generated 3D trajectory line, main leakage channel, and detour path are superimposed in different colors on the 3D geometric model of the mouth structure, and the output is a 3D leakage trajectory map including path annotations.

6. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements a method for on-site airtightness testing of the entrance of a civil defense project as described in any one of claims 1-5.

7. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements a method for on-site airtightness testing of the entrance of a civil defense project as described in any one of claims 1-5.

Images