A tunnel construction management method and system

Abstract

The present invention relates to the technical field of data processing, and specifically relates to a tunnel construction management method and system, including: obtaining a plurality of initial seismic sources and a plurality of current seismic sources from all positions according to the time series of vibration data at all positions; obtaining a risk degree quantization value at the current moment according to the number of current seismic sources and initial seismic sources and the vibration data of the vibration data time series; obtaining an average increase rate for each position according to the difference in vibration data at each moment of each position; obtaining the number of focusing behaviors of the seismic source according to the number size of the current seismic source and the initial seismic source; obtaining a risk degree corrected quantization value at the current moment according to the risk degree quantization value at the current moment, the number of focusing behaviors of the seismic source, and the average increase rate of each position; and judging whether there is a risk in the tunnel according to the risk degree corrected quantization value at the current moment. The present invention improves the safety of tunnel construction management.

Description

Technical Field

[0001] The present invention relates to the field of data processing, and in particular to a tunnel construction management method and system. Background Art

[0002] In areas with complex terrains such as mountains or canyons, these terrains are often difficult to cross. If highways are built around mountains and other terrains, huge human and financial resources will be consumed, and it will not be convenient for traffic. Opening a tunnel inside the mountains to build a highway can effectively improve traffic efficiency, reduce the cost of crossing mountains, further promote economic and cultural exchanges between various places, and improve the happiness of the people. Therefore, it is necessary to build tunnels in specific terrains. However, the difficulty and risk of building tunnels are much greater than building ordinary highways. In the process of building tunnels, many construction operations such as blasting, excavation, and support operations will cause different degrees of vibration in the tunnel. Among these vibrations, some vibrations are normal and safe, while some vibrations may cause the risk of tunnel collapse, endangering the personal safety of tunnel construction workers. Therefore, it is very necessary to conduct a risk assessment of the vibrations generated in the tunnel during tunnel construction. In the prior art, when evaluating whether the tunnel vibrations generated during tunnel construction may cause tunnel collapse, the collected vibration data are often compared with the safety standard range, and the collapse risk of vibration is evaluated by whether the vibration data is within the safety range. However, this approach is somewhat one-sided. The safety standards are not universally applicable and do not analyze the specific location of vibrations in the tunnel and the changes in vibrations. Some smaller vibrations may slowly trigger huge vibrations and cause collapse, reducing the safety of tunnel construction management. Summary of the invention

[0003] The present invention provides a tunnel construction management method and system to solve the problem that in the prior art, when evaluating whether the tunnel vibration generated during tunnel construction may cause tunnel collapse, the collected vibration data is often compared with the safety standard range, and the collapse risk of vibration is evaluated by whether the vibration data is within the safety range. However, this is somewhat one-sided, the safety standard is not universal, and the specific location of the vibration in the tunnel and the change of the vibration are not analyzed. Some smaller vibrations may also slowly cause huge vibrations and cause collapse, reducing the safety of tunnel construction management.

[0004] A tunnel construction management method and system of the present invention adopts the following technical solutions:

[0005] An embodiment of the present invention provides a tunnel construction management method, the method comprising the following steps:

[0006] Real-time acquisition of the time series of vibration data at several different locations in the tunnel; among which there is an adjacent relationship between different locations;

[0007] Obtain a number of initial seismic sources and a number of current seismic sources from all positions according to the time series of vibration data at all positions;

[0008] Obtain the quantified risk level value at the current moment according to the number of current seismic sources and initial seismic sources and the vibration data in the time series of vibration data;

[0009] Obtain the average increase rate of each position according to the difference in vibration data at each moment of each position; obtain the number of focusing behaviors of the seismic source according to the number of current seismic sources and initial seismic sources; obtain the corrected quantified risk level value at the current moment according to the quantified risk level value at the current moment, the number of focusing behaviors of the seismic source, and the average increase rate of each position;

[0010] Judge whether there is a risk in the tunnel according to the corrected quantified risk level value at the current moment.

[0011] Preferably, the step of obtaining a number of initial seismic sources and a number of current seismic sources from all positions according to the time series of vibration data at all positions includes the following specific steps:

[0012] In the first vibration data of the time series of vibration data at all positions, when the first vibration data at the z-th position is greater than the first vibration data at all adjacent positions of the z-th position at the same time, mark the z-th position as an initial seismic source;

[0013] In the last vibration data of the time series of vibration data at all positions, when the last vibration data at the z-th position is greater than the last vibration data at all adjacent positions of the z-th position at the same time, mark the z-th position as a current seismic source.

[0014] Preferably, the step of obtaining the quantified risk level value at the current moment according to the number of current seismic sources and initial seismic sources and the vibration data in the time series of vibration data includes the following specific steps:

[0015] Obtain the seismic source characteristic value according to the number of current seismic sources and initial seismic sources;

[0016] Obtain the quantified risk level value at the current moment according to the vibration data in the time series of vibration data of the current seismic source, the vibration data in the time series of vibration data of the initial seismic source, and the seismic source characteristic value.

[0017] Preferably, the formula for obtaining the quantified risk level value at the current moment according to the vibration data in the time series of vibration data of the current seismic source, the vibration data in the time series of vibration data of the initial seismic source, and the seismic source characteristic value is as follows:

[0018]

[0019] In the formula, F represents the quantified value of the risk level at the current moment; represents the mean value of the last vibration data in the time series of vibration data of all current seismic sources; g 2 represents the number of current seismic sources; represents the mean value of the first vibration data in the time series of vibration data of all initial seismic sources; g 1 represents the number of initial seismic sources; x represents the number of positions that are both initial seismic sources and current seismic sources; is the seismic source eigenvalue; th() is the hyperbolic tangent function.

[0020] Preferably, obtaining the average increase rate of each position according to the differences in the vibration data of each moment of each position includes the following specific steps:

[0021] Obtain the first-order difference sequence of the time series of vibration data at each position, and record the latter vibration data in the two vibration data corresponding to the last negative number in the first-order difference sequence as the reference vibration data;

[0022] Obtain the average increase rate of each position according to the difference between the vibration data at the current moment and the reference vibration data of each position.

[0023] Preferably, obtaining the average increase rate of each position according to the difference between the vibration data at the current moment and the reference vibration data of each position includes the following specific steps:

[0024] Calculate the difference between the vibration data at the current moment of the i-th position and the reference vibration data of the i-th position, and record the ratio of the difference to the time interval between the current moment of the i-th position and the moment where the reference vibration data is located as the average increase rate of the i-th position.

[0025] Preferably, obtaining the number of focusing behaviors of the seismic source according to the quantity relationship between the current seismic source and the initial seismic source includes the following specific steps:

[0026] If g 2 -g 1 ≥0, then set the number of focusing behaviors of the seismic source to the preset initial value T;

[0027] Among them, g 2 represents the number of current seismic sources, g 1 represents the number of initial seismic sources;

[0028] In the time series of vibration data at all positions, when the last vibration data of the y-th position is greater than the first vibration data of all positions adjacent to the y-th position at the same time, mark the y-th position as the target position;

[0029] If g 2 -g 1 <0, then set the number of focusing behaviors of the seismic source to the number of target positions.

[0030] Preferably, obtaining the risk degree correction quantization value at the current moment according to the risk degree quantization value at the current moment, the number of focusing behaviors of the seismic source, and the average increase rate of each position includes the following specific steps:

[0031] Record the result after normalizing the product of the risk degree quantization value at the current moment, the number of focusing behaviors of the seismic source, and the maximum value of the average increase rate of all positions as the risk degree correction quantization value at the current moment.

[0032] Preferably, judging whether there is a risk in the tunnel according to the risk degree correction quantization value at the current moment includes the following specific steps:

[0033] If the risk degree correction quantization value at the current moment is less than or equal to the preset threshold a, it is determined that the tunnel is in the first risk level at the current moment;

[0034] If the risk degree correction quantization value at the current moment is greater than the preset threshold a and less than or equal to the preset threshold b, it is determined that the tunnel is in the second risk level at the current moment;

[0035] If the risk degree correction quantization value at the current moment is greater than the preset threshold b, it is determined that the tunnel is in the third risk level at the current moment.

[0036] The present invention also proposes a tunnel construction management system, including a memory, a processor, and a computer program stored on the memory and executable on the processor. The processor executes the computer program stored in the memory to implement the steps of the foregoing tunnel construction management method.

[0037] The beneficial effects of the technical solution of the present invention are as follows: on the basis of only judging through the safety standard range in the prior art, the position where the vibration is located and the degree of change of the vibration at this position in the near future are further considered, the vibration risk is comprehensively evaluated and predicted, and it is judged whether it may cause a landslide risk. Risk assessment can be carried out more accurately and there is a certain function of predicting risks in advance. The vibration in the tunnel generated during tunnel construction is risk-assessed to judge whether it may cause a landslide risk in the tunnel, and the risk is warned in time to avoid accidents. According to the difference in vibration data at each moment at each position, the average increase rate at each position is obtained; according to the quantity size of the current seismic source and the initial seismic source, the number of focusing behaviors of the seismic source is obtained; according to the risk degree quantization value at the current moment, the number of focusing behaviors of the seismic source, and the average increase rate at each position, the risk degree correction quantization value at the current moment is obtained; according to the risk degree correction quantization value at the current moment, it is judged whether there is a landslide risk in the tunnel; the risk of tunnel construction is warned in advance, and the safety of tunnel construction management is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0039] Figure 1 It is a step flow chart of a tunnel construction management method of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0040] In order to further elaborate on the technical means and effects adopted by the present invention to achieve the intended invention purpose, the following will, in combination with the drawings and preferred embodiments, describe in detail the specific implementation manners, structures, features, and effects of a tunnel construction management method and system according to the present invention. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures, or characteristics in one or more embodiments can be combined in any suitable form.

[0041] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the technical field to which the present invention belongs.

[0042] The following will specifically describe the specific solutions of a tunnel construction management method and system provided by the present invention in conjunction with the drawings.

[0043] Please refer to Figure 1, which shows a step flow chart of a tunnel construction management method provided by an embodiment of the present invention. The method includes the following steps:

[0044] Step S001: Obtain the vibration data time series at several different positions in the tunnel in real time; among them, there is an adjacent relationship between different positions.

[0045] Obtain the vibration data time series at several different positions in the tunnel in real time; among them, there is an adjacent relationship between different positions.

[0046] It should be noted that: Different positions in the tunnel, such as the support points in the tunnel, the ore loosening points found during construction, etc., are called risk points. A vibrating wire sensor signal converter is installed at each position to detect the tunnel vibration information within a certain range of the position. In this embodiment, the position data is 50, the vibration data acquisition frequency is once per second, and the duration is from the current moment to the previous 5 minutes. This is used as an example for description. In areas without risk points, they are evenly arranged according to the sensing range of the vibrating wire sensor signal converter to capture all the vibration information in the tunnel. Among them, each vibration data time series undergoes Gaussian filtering denoising processing. The specified values in the industry's highway tunnel construction technical specifications are known. If the maximum vibration data in the vibration data time series exceeds the specified range, a high risk is directly reported, the alarm is immediately turned on, and the vibration data time series is not analyzed subsequently. Otherwise, those that meet the safety standard range continue with the subsequent risk assessment steps. By analyzing the vibration data, it is possible to judge in real time whether risk warning is needed to drive the construction team to evacuate, or to change the construction team's operation content according to the vibration intensity.

[0047] Among them, Gaussian filtering is a well-known technology, and the specific method is not introduced here.

[0048] It should be further noted that: If the monitoring areas of two vibrating wire sensor signal converters overlap, that is, due to the uniform arrangement of sensors, in order to detect the vibrations in all areas, there will be a small overlapping area between the detection areas of the sensors, then the positions of the two sensors are considered adjacent.

[0049] Step S002: According to the vibration data time series at all positions, obtain several initial vibration sources and several current vibration sources from all positions; obtain the risk degree quantization value at the current moment according to the number of current vibration sources and initial vibration sources and the vibration data of the vibration data time series.

[0050] It should be noted that during the tunnel construction process, the vibrations generated by the construction are within the safety standards. Logically, the vibration intensity should remain stable and gradually decrease as the vibration propagates. However, for some vibrations that may cause collapses inside the tunnel, even if they are within the safety standards, these vibrations may exhibit the following behaviors. Along a certain vibration propagation path, the vibration intensity on this path may increase rather than decrease during propagation, gradually growing larger. Then, a certain position with the largest vibration intensity under this path may be on the verge of risk. For a certain position, if the vibration intensity at this point or in this area shows a trend of gradually increasing and large vibrations converging towards the interior of the area, then this area may be on the verge of risk. Therefore, the risk can be evaluated based on the changing trends of vibrations over time and along the propagation path.

[0051] It should be noted that in the time series of vibration data at all positions inside the tunnel, find the current seismic source and the seismic sources within the previous 5 minutes, and mark the positions of the seismic sources; and record the positions simultaneously marked by the current seismic source and the seismic sources within the previous 5 minutes.

[0052] In the first vibration data of the time series of vibration data at all positions, when the first vibration data at the z-th position is simultaneously greater than the first vibration data at all adjacent positions of the z-th position, mark the z-th position as the initial seismic source;

[0053] In the last vibration data of the time series of vibration data at all positions, when the last vibration data at the z-th position is simultaneously greater than the last vibration data at all adjacent positions of the z-th position, mark the z-th position as the current seismic source.

[0054] Thus, a number of initial seismic sources and a number of current seismic sources are obtained.

[0055] It should be noted that during normal construction vibrations, the seismic sources usually do not change much, and their vibration intensities are also in a relatively stable state. If the construction vibrations cause vibrations in some structurally sparse areas inside the tunnel, or if the tunnel instability leads to risk vibrations in some areas, and the vibrations propagate to some risk areas and reach the critical value of the risk areas, it may cause changes in some seismic sources or even increase the number of seismic sources, and the vibration intensities of the seismic sources will also change.

[0056] The calculation method for the quantified value of the risk level at the current moment is as follows:

[0057]

[0058] In the formula, F represents the quantified value of the risk level at the current moment; represents the mean value of the last vibration data in the time series of vibration data of all current seismic sources; g 2 represents the number of current seismic sources; represents the mean value of the first vibration data in the time series of vibration data of all initial seismic sources; g 1 represents the number of initial seismic sources; x represents the number of positions that are the same for both the initial seismic source and the current seismic source; is the eigenvalue of the seismic source; th() is the hyperbolic tangent function, which is used to normalize the data value between 0 and 1.

[0059] It should be noted that when the denominator in the formula is 0, let the denominator be 1, and this is used as an example for description to ensure the formula holds.

[0060] It should be noted that represents the difference between the mean value of the last vibration data in the time series of vibration data of all current seismic sources and the mean value of the first vibration data in the time series of vibration data of all initial seismic sources. The larger this value is, the greater the average vibration intensity at the current moment relative to five minutes ago, and the more likely there is a risk of collapse at the current moment; th(g 2 -g 1 ) represents the difference between the number of previous seismic sources and the number of initial seismic sources. The larger this value is, the more seismic sources have increased at the current moment relative to five minutes ago, and the more likely there is a risk of collapse at the current moment; x represents the number of positions that are the same for both the initial seismic source and the current seismic source. The smaller this value is, the more seismic sources have changed at the current moment relative to five minutes ago, and the more likely there is a risk of collapse at the current moment.

[0061] It should be noted that the quantization value comprehensively represents the risk level at the current moment, but the intensity change of some points within these five minutes and whether there is an aggregation behavior of seismic sources have not been considered yet, and further comprehensive calculation is still required.

[0062] Thus, the quantization value of the risk level at the current moment is obtained.

[0063] Step S003: According to the differences in the vibration data of each moment at each position, obtain the average increase rate of each position; according to the number of the current seismic source and the initial seismic source, obtain the number of times of the focusing behavior of the seismic source; according to the quantization value of the risk level at the current moment, the number of times of the focusing behavior of the seismic source, and the average increase rate of each position, obtain the corrected quantization value of the risk level at the current moment.

[0064] It should be noted that if the data obtained from some sensors shows a continuous sharp increase, it is necessary to consider whether this point is a high-risk point; and even if the number of seismic sources decreases at the current moment relative to five minutes ago, it may be a behavior of the aggregation of seismic sources, and it is necessary to further analyze this possibility.

[0065] It should be noted that first, the data in the time series of vibration data at each position is judged to observe whether there is a significant increase.

[0066] Obtain the first-order difference sequence of the vibration data time series at each position, and record the latter vibration data of the two vibration data corresponding to the last negative number in the first-order difference sequence as the reference vibration data.

[0067] It should be noted that the data change after the reference vibration data in the vibration data time series at each position is an increasing change.

[0068] Calculate the difference between the vibration data at the current moment at the i-th position and the reference vibration data at the i-th position, and record the ratio of the difference to the time interval between the current moment at the i-th position and the moment where the reference vibration data is located as the average increase rate at the i-th position;

[0069] The calculation method of the average increase rate at the i-th position is as follows:

[0070]

[0071] In the formula, v i represents the average increase rate at the i-th position; d i represents the vibration data at the current moment at the i-th position; d i ' represents the reference vibration data at the i-th position; t i represents the time interval between the current moment at the i-th position and the moment where the reference vibration data is located.

[0072] Record the maximum value of the average increase rates at all positions as MAX{v};

[0073] It should be noted that determine whether there is an aggregation behavior of tunnel seismic sources at the current moment in the tunnel; if g 2 -g 1 ≥0, it means that the current number of seismic sources has not decreased compared with the previous one, so the number of aggregation behaviors of seismic sources remains unchanged (if there is an aggregation behavior of seismic sources, that is, the case of the merger of seismic source points, it will lead to a decrease in the number of seismic sources); otherwise, it is necessary to further judge whether there is a situation of aggregation of vibration source points.

[0074] Preset the initial value T to 1, and take this as an example for description.

[0075] If g 2 -g 1 ≥0, then set the number of focusing behaviors of seismic sources to the preset initial value T;

[0076] Among them, g 2 represents the current number of seismic sources, and g 1 represents the initial number of seismic sources;

[0077] In the time series of vibration data at all positions, when the last vibration data at the y-th position is greater than the first vibration data at the y-th position and all adjacent positions to the y-th position, mark the y-th position as the target position;

[0078] If g 2 -g 1 < 0, then set the number of focusing behaviors of the seismic source to the number of target positions.

[0079] Record the result after normalizing the product of the risk degree quantization value at the current moment, the number of focusing behaviors of the seismic source, and the maximum value of the average increase rate at all positions as the risk degree correction quantization value at the current moment;

[0080] The calculation method of the risk degree correction quantization value at the current moment is as follows:

[0081] FX = sigmoid{F * T * MAX{v}}

[0082] In the formula, FX represents the risk degree correction quantization value at the current moment; F represents the risk degree quantization value at the current moment; T represents the number of focusing behaviors of the seismic source; MAX{v} represents the maximum value of the average increase rate at all positions; sigmoid{} is the normalization function.

[0083] It should be noted that T represents the number of focusing behaviors of the seismic source, and this value is a correction to the F value; MAX{v} represents the maximum value of the average increase rate at all positions. The larger this value is, the more likely it is to have a vibration intensity exceeding the safety standard, and the more necessary it is to give an early warning. This value and the F value jointly represent the current risk degree and the early warning of the subsequent risk degree.

[0084] Thus, the risk degree correction quantization value at the current moment is obtained.

[0085] Step S004: Determine whether there is a risk in the tunnel according to the risk degree correction quantization value at the current moment.

[0086] Evaluate the current risk degree of the tunnel according to the risk degree correction quantization value at the current moment obtained in the above steps;

[0087] The preset threshold a is 0.65, and the preset threshold b is 0.7;

[0088] If the risk degree correction quantization value at the current moment is less than or equal to the preset threshold a, it is determined that the tunnel is in the first risk level at the current moment; it is necessary for the staff to be more vigilant;

[0089] When the corrected quantization value of the risk level at the current moment is greater than the preset threshold a and less than or equal to the preset threshold b, it is determined that the tunnel is at the second risk level at the current moment; that is, there is a certain risk in the tunnel at present, and an emergency plan for preventing risks needs to be set up.

[0090] When the corrected quantization value of the risk level at the current moment is greater than the preset threshold b, it is determined that the tunnel is at the third risk level at the current moment; that is, there is a risk in the tunnel at present, and it is necessary to arrange for the staff to retreat in an orderly manner in a timely manner.

[0091] It should be noted that according to the above steps, the current risk assessment level is obtained. If there is a certain risk, it may be that the current construction behavior is inappropriate. It is necessary to give a warning to the relevant staff, remind them to change the operation strategy, let the construction personnel be more vigilant, reduce the work intensity during the process of evacuating to the safe area, and observe the change of the subsequent risk level; if it is prompted that there is a risk in the tunnel at present, immediately notify the relevant personnel to arrange for the construction personnel to retreat reasonably and orderly to the safe area outside the tunnel to prevent accidents from occurring.

[0092] The embodiment of the present invention also provides a tunnel construction management system, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the above steps S001 to S004 are implemented.

[0093] So far, this embodiment is completed.

[0094] The above are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A tunnel construction management method, characterized in that: The method comprises the following steps: Real-time acquisition of the time series of vibration data at several different locations in the tunnel; among which there is an adjacent relationship between different locations; According to the time series of the vibration data at all the locations, a number of initial earthquake sources and a number of current earthquake sources are obtained from all the locations; According to the number of current earthquake sources and initial earthquake sources and the vibration data of the vibration data time series, a quantitative value of the risk level at the current moment is obtained; According to the difference of vibration data at each position at each moment, the average increase rate of each position is obtained; according to the number of current and initial seismic sources, the number of focusing behaviors of the seismic sources is obtained; according to the quantitative value of the risk level at the current moment, the number of focusing behaviors of the seismic sources and the average increase rate of each position, the corrected quantitative value of the risk level at the current moment is obtained; Correct the quantitative value according to the risk level at the current moment to determine whether the tunnel is at risk; The specific steps for obtaining the focusing behavior times of the earthquake source are as follows: If g2-g1≥0, the number of focusing behaviors of the source is set to the preset initialization value T; Among them, g2 represents the number of current earthquake sources, and g1 represents the number of initial earthquake sources; In the vibration data time series at all positions, when the last vibration data of the yth position is greater than the first vibration data of all positions adjacent to the yth position, the yth position is marked as the target position; If g2-g1<0, let the number of focusing actions of the source be the number of target positions.

2. A tunnel construction management method according to claim 1, characterized in that: The method of obtaining a plurality of initial earthquake sources and a plurality of current earthquake sources from all locations according to the time series of the earthquake data at all locations includes the following specific steps: In the first vibration data of the vibration data time series at all positions, when the first vibration data at the z-th position is simultaneously greater than the first vibration data of all positions adjacent to the z-th position, the z-th position is marked as the initial earthquake source; In the last vibration data of the vibration data time series at all positions, when the last vibration data at the z-th position is simultaneously greater than the last vibration data of all positions adjacent to the z-th position, the z-th position is marked as the current earthquake source.

3. A tunnel construction management method according to claim 2, characterized in that: The step of obtaining the quantitative value of the risk level at the current moment according to the number of the current earthquake source and the initial earthquake source and the vibration data of the vibration data time series sequence includes the following specific steps: According to the number of current earthquake sources and initial earthquake sources, the earthquake source characteristic value is obtained; According to the vibration data of the vibration data time series sequence of the current earthquake source, the vibration data of the vibration data time series sequence of the initial earthquake source and the earthquake source characteristic value, a quantitative value of the risk degree at the current moment is obtained.

4. A tunnel construction management method according to claim 3, characterized in that: The specific formula for obtaining the quantitative value of the risk level at the current moment according to the vibration data of the current earthquake source vibration data time series sequence, the vibration data of the initial earthquake source vibration data time series sequence and the earthquake source characteristic value is as follows: In the formula, F represents the quantitative value of the risk level at the current moment; represents the mean value of the last vibration data in the vibration data time series of all current earthquake sources; g2 represents the number of current earthquake sources; represents the mean of the first vibration data of the vibration data time series of all initial seismic sources; g1 represents the number of initial seismic sources; x represents the number of locations that are both initial seismic sources and current seismic sources; is the source characteristic value; th() is the hyperbolic tangent function.

5. A tunnel construction management method according to claim 1, characterized in that: The method of obtaining the average increase rate of each position according to the difference of the vibration data of each position at each moment includes the following specific steps: Obtain a first-order difference sequence of the vibration data time series at each position, and record the latter vibration data of the two vibration data corresponding to the last negative number in the first-order difference sequence as the reference vibration data; According to the difference between the vibration data of each position at the current moment and the reference vibration data, the average increase rate of each position is obtained.

6. A tunnel construction management method according to claim 5, characterized in that: The method of obtaining the average increase rate of each position according to the difference between the vibration data at the current moment of each position and the reference vibration data includes the following specific steps: The difference between the vibration data at the current moment of the i-th position and the reference vibration data at the i-th position is calculated, and the ratio of the difference to the time interval between the current moment of the i-th position and the moment of the reference vibration data is recorded as the average increase rate of the i-th position.

7. A tunnel construction management method according to claim 1, characterized in that: The step of obtaining the corrected quantitative value of the risk level at the current moment according to the quantitative value of the risk level at the current moment, the number of focusing behaviors of the earthquake source, and the average increase rate of each position includes the following specific steps: The result of normalizing the product of the current risk level quantization value, the number of focusing behaviors of the earthquake source, and the maximum value of the average increase rate of all positions is recorded as the current risk level modified quantization value.

8. A tunnel construction management method according to claim 1, characterized in that: The specific steps of correcting the quantitative value according to the risk level at the current moment and judging whether the tunnel has risk are as follows: If the corrected quantized value of the risk level at the current moment is less than or equal to the preset threshold a, it is determined that the tunnel is at the first risk level at the current moment; If the corrected quantified value of the risk level at the current moment is greater than the preset threshold a and less than or equal to the preset threshold b, it is determined that the tunnel is at the second risk level at the current moment; If the corrected quantized value of the risk level at the current moment is greater than the preset threshold b, it is determined that the tunnel is at the third risk level at the current moment.

9. A tunnel construction management system, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the computer program is executed by a processor, the steps of a tunnel construction management method as described in any one of claims 1 to 8 are implemented.

Images