Damage analysis method and system for tunnel invert in gypsum-salt strata

Abstract

The present invention discloses an invert arch damage analysis method and system for a gypsum-salt stratum tunnel, the method comprising: obtaining displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert arch of a target gypsum-salt tunnel; constructing an invert arch deformation curve of the target section invert arch according to the displacement monitoring data; forming multiple calculation conditions; performing calculation analysis on multiple calculation conditions to obtain a displacement curve to be located; selecting a selected displacement curve; and using the selected displacement curve as a damage condition. The present invention simplifies the relatively complex erosion and expansion process of the base of a gypsum-salt stratum tunnel, so that the distribution of base defects of a gypsum-salt stratum tunnel can be analyzed through conventional invert arch displacement monitoring means, and the location of the defects can be found more accurately to perform corresponding repair operations, effectively improving the operational safety of the gypsum-salt tunnel, and at the same time providing corresponding guidance for the construction of the gypsum-salt tunnel.

Description

Technical Field

[0001] The present invention relates to the field of intelligent analysis technology, and in particular to a method and system for analyzing damage of an invert arch of a tunnel in a gypsum-salt stratum. Background Art

[0002] The geological conditions and geological environment of the tunnel in the gypsum-salt formation are extremely complex, and the lithology of the formation is changeable. Due to the strong solubility and expansion of gypsum-salt, the excavation of the tunnel in the gypsum-salt formation forms an underground dynamic water environment. When the gypsum-salt meets water, it will first dissolve and form cavities. Secondly, the gypsum-salt will also expand when it meets water, causing expansion, extrusion and cracking of the tunnel in the gypsum-salt formation, especially at the location of the tunnel invert, which is prone to dissolution cavities and expansion and extrusion. Tunnel construction inevitably passes through the gypsum-salt formation, which leads to major disasters and diseases in the construction and later operation of the gypsum-salt tunnel. Therefore, it is particularly important and practical to construct a method and system for analyzing the damage of the invert of the gypsum-salt formation tunnel. The base erosion of the gypsum-salt tunnel during operation is relatively complex. If the dissolution cavities or expansion crystallization of the base can be discovered in time and corresponding measures can be taken for maintenance, the operational safety of the tunnel in the gypsum-salt formation can be greatly improved. Summary of the invention

[0003] In order to at least overcome the above-mentioned deficiencies in the prior art, the purpose of the present application is to provide a method and system for analyzing damage of an invert arch of a tunnel in a gypsum-salt stratum.

[0004] In a first aspect, an embodiment of the present application provides a method for analyzing damage of an invert arch of a tunnel in a gypsum-salt stratum, comprising:

[0005] Obtaining displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert;

[0006] Constructing an invert deformation curve of the target section invert according to the displacement monitoring data; the vertical axis of the invert deformation curve is the vertical displacement value, and the horizontal axis of the invert deformation curve is the central angle corresponding to the monitoring point; extracting at least one characteristic point from the invert deformation curve, and combining the parts of the target section invert corresponding to the characteristic point that may have dissolution cavities or expansion stress to form a plurality of calculation conditions; the expansion stress is calculated according to the surrounding rock data;

[0007] Performing calculation analysis on a plurality of the calculation conditions, obtaining vertical displacement curves of the target section invert corresponding to the plurality of the calculation conditions to form a displacement curve to be positioned;

[0008] Comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve;

[0009] The calculation condition corresponding to the selected displacement curve is used as the damage condition of the target section invert of the target gypsum-salt tunnel.

[0010] When implementing the embodiments of the present application, it is necessary to obtain displacement monitoring data through monitoring means. On the same section, each monitoring point will collect corresponding displacement data, and each monitoring point will correspond to a different central angle; and for the lining stiffness, it can be calculated by the thickness and elastic modulus of the lining; the surrounding rock data needs to include the relevant parameters of the surrounding surrounding rock, as well as the stress conditions generated after the expansion and crystallization state may occur. These data can be obtained from the design-related data of the target gypsum salt tunnel, which belongs to the prior art and is not limited by the embodiments of the present application. For the deformation curve of the inverted arch that has been obtained, the main purpose of the embodiments of the present application is to analyze the possible base defect area from the deformation curve of the inverted arch for subsequent tunnel maintenance.

[0011] In the embodiment of the present application, different base defect conditions produce different inverted arch deformation curves. Since the mechanism of action is relatively complex, it is necessary to first perform feature analysis on the inverted arch deformation curve in the embodiment of the present application to find out the suspected defective parts, i.e., feature points, so as to find out the most correct set of combined working conditions by exhaustive calculation, wherein the defects mainly include dissolution cavities and expansion crystallization. After performing calculation and analysis on each calculated working condition, the displacement curves to be positioned corresponding to different working conditions can be obtained, which is also a curve with the vertical axis as the vertical displacement value and the horizontal axis as the central angle; after comparing the different displacement curves to be positioned and the inverted arch deformation curve, the working condition corresponding to the most similar displacement curve to be positioned is the defect distribution that is most likely to occur in the current section; at this time, the defects can be repaired according to this defect distribution, thereby improving the operational safety of the gypsum-salt tunnel. The embodiment of the present application simplifies the relatively complex dissolution and expansion process of the base of the gypsum-salt stratum tunnel through the above-mentioned technical solution, so that the distribution of base defects of the gypsum-salt stratum tunnel can be analyzed through conventional arch displacement monitoring means, and the location of the defects can be found more accurately and then the corresponding repair work can be carried out, which effectively improves the operational safety of the gypsum-salt tunnel and can provide corresponding guidance for the construction of the gypsum-salt tunnel.

[0012] In a possible implementation manner, extracting at least one characteristic point from the invert deformation curve includes:

[0013] When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points;

[0014] When there is only one maximum point and one minimum point in the invert deformation curve, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points.

[0015] In a possible implementation, multiple calculation conditions are formed by combining the location of the target section invert corresponding to the characteristic point with possible dissolution cavities or expansion stress, including:

[0016] Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality;

[0017] The state information of all the characteristic points is combined and exhaustively enumerated to form a plurality of the calculation conditions.

[0018] In a possible implementation, performing calculation analysis on the plurality of calculation conditions includes:

[0019] Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel;

[0020] Obtaining state information corresponding to different characteristic points in the calculation condition;

[0021] For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0;

[0022] For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point;

[0023] Perform mechanical calculations on the mechanical model with updated invert arch load information.

[0024] In a possible implementation, comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve includes:

[0025] Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve;

[0026] By calculating the cosine similarity between the proposed slope curve and the reference slope curve, the proposed slope curve with the highest cosine similarity is used as the selected displacement curve. In a second aspect, the embodiment of the present application also provides a gypsum-salt stratum tunnel invert damage analysis system, comprising:

[0027] An acquisition unit is configured to acquire displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert;

[0028] A construction unit is configured to construct an invert deformation curve of the target section invert according to the displacement monitoring data; the vertical axis of the invert deformation curve is the vertical displacement value, and the horizontal axis of the invert deformation curve is the central angle corresponding to the monitoring point;

[0029] An extraction unit is configured to extract at least one characteristic point from the deformation curve of the invert, and form a plurality of calculation conditions by combining the parts of the target section invert corresponding to the characteristic point with possible dissolution cavities or expansion stress; the expansion stress is calculated according to the surrounding rock data;

[0030] A calculation unit is configured to perform calculation analysis on a plurality of the calculation conditions, and obtain vertical displacement curves of the target section invert corresponding to the plurality of the calculation conditions to form a displacement curve to be determined;

[0031] a comparison unit configured to compare the plurality of displacement curves to be positioned with the invert deformation curve, and select the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve;

[0032] The damage unit is configured to use the calculation condition corresponding to the selected displacement curve as the damage condition of the target section invert of the target gypsum-salt tunnel.

[0033] In a possible implementation, the extracting unit is further configured to:

[0034] When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points;

[0035] When there is only one maximum point and one minimum point in the invert deformation curve, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points.

[0036] In a possible implementation, the extracting unit is further configured to:

[0037] Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality;

[0038] The state information of all the characteristic points is combined and exhaustively enumerated to form a plurality of the calculation conditions.

[0039] In a possible implementation manner, the computing unit is further configured to:

[0040] Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel;

[0041] Obtaining state information corresponding to different characteristic points in the calculation condition;

[0042] For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0;

[0043] For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point;

[0044] Perform mechanical calculations on the mechanical model with updated invert arch load information.

[0045] In a possible implementation, the comparison unit is further configured to:

[0046] Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve;

[0047] The cosine similarity between the pending slope curve and the reference slope curve is calculated, and the pending slope curve with the highest cosine similarity is used as the selected displacement curve.

[0048] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0049] The method and system for analyzing the damage of the invert arch of a tunnel in gypsum-salt stratum of the present invention simplifies the relatively complicated dissolution and expansion process of the base of the tunnel in gypsum-salt stratum through the above-mentioned technical scheme, so that the distribution of base defects of the tunnel in gypsum-salt stratum can be analyzed through conventional invert arch displacement monitoring means, and the location of the defects can be found more accurately and then corresponding repair operations can be carried out, which effectively improves the operational safety of the gypsum-salt tunnel and can provide corresponding guidance for the construction of the gypsum-salt tunnel. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] The drawings described herein are used to provide a further understanding of the embodiments of the present invention, constitute a part of this application, and do not constitute a limitation of the embodiments of the present invention. In the drawings:

[0051] Figure 1 This is a schematic diagram of the steps of the method of the embodiment of the present application;

[0052] Figure 2 This is a deformation curve diagram of the inverted arch in the embodiment of the present application;

[0053] Figure 3 This is a deformation curve diagram of the inverted arch according to an embodiment of the present application. DETAILED DESCRIPTION

[0054] To make the purpose, technical scheme and advantages of the embodiment of the present application clearer, the technical scheme in the embodiment of the present application will be clearly and completely described below in conjunction with the drawings in the embodiment of the present application. It should be understood that the drawings in the present application only serve the purpose of explanation and description and are not used to limit the scope of protection of the present application. In addition, it should be understood that the schematic drawings are not drawn in real proportion. The flowchart used in this application shows the operations implemented according to some embodiments of the embodiment of the present application. It should be understood that the operations of the flowchart can be implemented out of sequence, and the steps without logical context can be reversed in order or implemented simultaneously. In addition, those skilled in the art can add one or more other operations to the flowchart under the guidance of the content of the present application, or remove one or more operations from the flowchart.

[0055] In addition, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application described and shown in the drawings here can be arranged and designed in various configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the application claimed for protection, but merely represents the selected embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without making creative work belong to the scope of protection of the present application.

[0056] Please refer to Figure 1 , which is a flow chart of a method for analyzing damage of an invert arch of a tunnel in a gypsum-salt stratum provided by an embodiment of the present invention. Further, the method for analyzing damage of an invert arch of a tunnel in a gypsum-salt stratum may specifically include the contents described in the following steps S1 to S6.

[0057] S1: Acquire displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert;

[0058] S2: constructing an invert deformation curve of the target section invert according to the displacement monitoring data; the vertical axis of the invert deformation curve is the vertical displacement value, and the horizontal axis of the invert deformation curve is the central angle corresponding to the monitoring point;

[0059] S3: extracting at least one characteristic point from the deformation curve of the invert, and combining the parts of the target section invert corresponding to the characteristic point with possible dissolution cavities or expansion stress to form multiple calculation conditions; the expansion stress is calculated according to the surrounding rock data;

[0060] S4: performing calculation analysis on the plurality of calculation conditions, obtaining vertical displacement curves of the target section invert corresponding to the plurality of calculation conditions to form a displacement curve to be positioned;

[0061] S5: comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve;

[0062] S6: Using the calculation condition corresponding to the selected displacement curve as the damage condition of the target section invert of the target gypsum-salt tunnel.

[0063] When implementing the embodiments of the present application, it is necessary to obtain displacement monitoring data through monitoring means. On the same section, each monitoring point will collect corresponding displacement data, and each monitoring point will correspond to a different central angle; and for the lining stiffness, it can be calculated by the thickness and elastic modulus of the lining; the surrounding rock data needs to include the relevant parameters of the surrounding surrounding rock, as well as the stress conditions generated after the expansion and crystallization state may occur. These data can be obtained from the design-related data of the target gypsum salt tunnel, which belongs to the prior art and is not limited by the embodiments of the present application. For the deformation curve of the inverted arch that has been obtained, the main purpose of the embodiments of the present application is to analyze the possible base defect area from the deformation curve of the inverted arch for subsequent tunnel maintenance.

[0064] In the embodiment of the present application, different base defect conditions produce different inverted arch deformation curves. Since the mechanism of action is relatively complex, it is necessary to first perform feature analysis on the inverted arch deformation curve in the embodiment of the present application to find out the suspected defective parts, i.e., feature points, so as to find out the most correct set of combined working conditions by exhaustive calculation, wherein the defects mainly include dissolution cavities and expansion crystallization. After performing calculation and analysis on each calculated working condition, the displacement curves to be positioned corresponding to different working conditions can be obtained, which is also a curve with the vertical axis as the vertical displacement value and the horizontal axis as the central angle; after comparing the different displacement curves to be positioned and the inverted arch deformation curve, the working condition corresponding to the most similar displacement curve to be positioned is the defect distribution that is most likely to occur in the current section; at this time, the defects can be repaired according to this defect distribution, thereby improving the operational safety of the gypsum-salt tunnel. The embodiment of the present application simplifies the relatively complex dissolution and expansion process of the base of the gypsum-salt stratum tunnel through the above-mentioned technical solution, so that the distribution of base defects of the gypsum-salt stratum tunnel can be analyzed through conventional arch displacement monitoring means, and the location of the defects can be found more accurately and then the corresponding repair work can be carried out, which effectively improves the operational safety of the gypsum-salt tunnel and can provide corresponding guidance for the construction of the gypsum-salt tunnel.

[0065] In a possible implementation manner, extracting at least one characteristic point from the invert deformation curve includes:

[0066] When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points;

[0067] When there is only one maximum point and one minimum point in the invert deformation curve, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points.

[0068] When implementing the embodiments of the present application, the inventors found that generally speaking, the locations where defects such as dissolution voids or expansion crystals appear generally correspond to the maximum points, minimum points and 0 points of the deformation curve of the invert arch. For the working conditions with at least two maximum points or at least two minimum points, the cause of their occurrence is generally due to at least one defect appearing on one side of the invert arch, causing asymmetric deformation of the invert arch, so the 0 point can be ignored at this time to reduce the number of working conditions; and for the working conditions with only one maximum point and one minimum point, there may be a defect at the 0 point, which is equivalent to the local application of a radial force away from the center of the circle, thereby changing the original uplift deformation into a settlement deformation, so the 0 point needs to be considered.

[0069] In a possible implementation, multiple calculation conditions are formed by combining the location of the target section invert corresponding to the characteristic point with possible dissolution cavities or expansion stress, including:

[0070] Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality;

[0071] The state information of all the characteristic points is combined and exhaustively enumerated to form a plurality of the calculation conditions.

[0072] In a possible implementation, performing calculation analysis on the plurality of calculation conditions includes:

[0073] Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel;

[0074] Obtaining state information corresponding to different characteristic points in the calculation condition;

[0075] For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0;

[0076] For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point;

[0077] Perform mechanical calculations on the mechanical model with updated invert arch load information.

[0078] When the embodiment of the present application is implemented, when performing calculation and analysis, it is necessary to first construct a mechanical model of the target gypsum-salt tunnel. The calculation process can be divided into two ways, one of which can be calculated by the inverted arch mechanical model, and the other can be calculated by the finite element model. The mechanical model described in the embodiment of the present application generally adopts at least one of the inverted arch mechanical model and the finite element model. In the case where the state information corresponding to the characteristic point in the calculation condition has been known, it is necessary to perform different processing for the characteristic points of different state information; for the characteristic point of the erosion cavity, it is considered that the bottom of the inverted arch at this position has lost support, so the corresponding surrounding rock load is set to 0, and for the characteristic point whose state information is expansion stress, there will be expansion stress at its bottom, so the expansion stress and surrounding rock load need to be superimposed to perform the next step of calculation.

[0079] In a possible implementation, comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve includes:

[0080] Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve;

[0081] The cosine similarity between the pending slope curve and the reference slope curve is calculated, and the pending slope curve with the highest cosine similarity is used as the selected displacement curve.

[0082] When the embodiment of the present application is implemented, in the process of calculating the displacement curve to be positioned, the actual calculated curve amplitude will be slightly different from the monitored curve amplitude, which is caused by the complex formation conditions of the gypsum-salt tunnel itself. Therefore, in the embodiment of the present application, the curve is converted into a slope curve by means of a first-order derivative, which can reduce the impact of the change in the curve amplitude and pay more attention to the trend of the curve. In the embodiment of the present application, the method of calculating the curve similarity can adopt cosine similarity, which also pays more attention to the curve trend. The calculation process belongs to the prior art, and the embodiment of the present application does not make much limitation.

[0083] In a specific implementation, see Figure 2 and Figure 3 , respectively show the displacement curves of the bottom of the invert under two different working conditions, where Figure 2 The figure shows the deformation curve of the inverted arch calculated by simulation in the embodiment of the present application when the cavities appear at the -40° and 40° positions of the inverted arch at the same time, wherein the curves of different colors represent the deformation curves of the inverted arch under different lining elastic moduli. Figure 3The figure shows the deformation curve of the inverted arch calculated by simulation in the embodiment of the present application when expansion force occurs at both the -40° and 40° positions of the inverted arch, wherein curves of different colors represent the deformation curves of the inverted arch under different lining elastic moduli. Figure 2 and Figure 3 It can be seen that in these two working conditions, the defect position is at the maximum value of the curve. Similarly, with the changes in the stiffness and expansion force of the surrounding rock and lining, the defect will also appear at the position of 0 and the minimum value.

[0084] Based on the same inventive concept, the embodiment of the present application also provides a gypsum-salt stratum tunnel invert damage analysis system, comprising:

[0085] An acquisition unit is configured to acquire displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert;

[0086] A construction unit is configured to construct an invert deformation curve of the target section invert according to the displacement monitoring data; the vertical axis of the invert deformation curve is the vertical displacement value, and the horizontal axis of the invert deformation curve is the central angle corresponding to the monitoring point;

[0087] An extraction unit is configured to extract at least one characteristic point from the deformation curve of the invert, and form a plurality of calculation conditions by combining the parts of the target section invert corresponding to the characteristic point with possible dissolution cavities or expansion stress; the expansion stress is calculated according to the surrounding rock data;

[0088] A calculation unit is configured to perform calculation analysis on a plurality of the calculation conditions, and obtain vertical displacement curves of the target section invert corresponding to the plurality of the calculation conditions to form a displacement curve to be determined;

[0089] a comparison unit configured to compare the plurality of displacement curves to be positioned with the invert deformation curve, and select the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve;

[0090] The damage unit is configured to use the calculation condition corresponding to the selected displacement curve as the damage condition of the target section invert of the target gypsum-salt tunnel.

[0091] In a possible implementation, the extracting unit is further configured to:

[0092] When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points;

[0093] When there is only one maximum point and one minimum point in the invert deformation curve, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points.

[0094] In a possible implementation, the extracting unit is further configured to:

[0095] Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality;

[0096] The state information of all the characteristic points is combined and exhaustively enumerated to form a plurality of the calculation conditions.

[0097] In a possible implementation manner, the computing unit is further configured to:

[0098] Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel;

[0099] Obtaining state information corresponding to different characteristic points in the calculation condition;

[0100] For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0;

[0101] For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point;

[0102] Perform mechanical calculations on the mechanical model with updated invert arch load information.

[0103] In a possible implementation, the comparison unit is further configured to:

[0104] Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve;

[0105] The cosine similarity between the pending slope curve and the reference slope curve is calculated, and the pending slope curve with the highest cosine similarity is used as the selected displacement curve.

[0106] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0107] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, or it can be an electrical, mechanical or other form of connection.

[0108] The units described as separate components may or may not be physically separated. As units, it is obvious that a person of ordinary skill in the art can realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0109] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0110] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a grid device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0111] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for analyzing damage of tunnel inverts in gypsum-salt strata, characterized in that: include: Obtaining displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert; Constructing an inverted arch deformation curve of the target section inverted arch according to the displacement monitoring data; the vertical axis of the inverted arch deformation curve is the vertical displacement value, and the horizontal axis of the inverted arch deformation curve is the central angle corresponding to the monitoring point; Extracting at least one characteristic point from the deformation curve of the invert, and combining the corrosion cavities or expansion stresses existing in the part of the target section invert corresponding to the characteristic point to form a plurality of calculation conditions; The expansion stress is calculated based on the surrounding rock data; Performing calculation analysis on a plurality of the calculation conditions, obtaining vertical displacement curves of the target section invert corresponding to the plurality of the calculation conditions to form a displacement curve to be positioned; Comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve; The calculation condition corresponding to the selected displacement curve is used as the damage condition of the target section invert of the target gypsum-salt tunnel; Extracting at least one characteristic point from the invert deformation curve comprises: When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points; When there is only one maximum point and one minimum point in the deformation curve of the invert, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points; The dissolution cavities or expansion stresses existing in the position of the target section invert corresponding to the characteristic point are combined to form multiple calculation conditions including: Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality; Combining and exhaustively enumerating the state information of all the characteristic points to form a plurality of the calculation conditions; The calculation analysis of the plurality of calculation conditions comprises: Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel; Obtaining state information corresponding to different characteristic points in the calculation condition; For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0; For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point; Perform mechanical calculations on the mechanical model with updated invert arch load information.

2. The method for analyzing damage of an invert arch of a tunnel in gypsum-salt stratum according to claim 1, characterized in that: Comparing the plurality of displacement curves to be positioned with the invert deformation curve, and selecting the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve comprises: Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve; The cosine similarity between the pending slope curve and the reference slope curve is calculated, and the pending slope curve with the highest cosine similarity is used as the selected displacement curve.

3. The damage analysis system of tunnel invert in gypsum-salt stratum is characterized by: include: An acquisition unit is configured to acquire displacement monitoring data, lining stiffness data and surrounding rock data of a target section invert of a target gypsum-salt tunnel; the displacement monitoring data includes vertical displacement data corresponding to a plurality of different central angles on the target section invert; A construction unit is configured to construct an invert deformation curve of the target section invert according to the displacement monitoring data; the vertical axis of the invert deformation curve is the vertical displacement value, and the horizontal axis of the invert deformation curve is the central angle corresponding to the monitoring point; An extraction unit is configured to extract at least one characteristic point from the deformation curve of the invert, and to form a plurality of calculation conditions by combining the dissolution cavities or expansion stresses existing at the position of the target section invert corresponding to the characteristic point; The expansion stress is calculated based on the surrounding rock data; A calculation unit is configured to perform calculation analysis on a plurality of the calculation conditions, and obtain vertical displacement curves of the target section invert corresponding to the plurality of the calculation conditions to form a displacement curve to be determined; a comparison unit configured to compare the plurality of displacement curves to be positioned with the invert deformation curve, and select the displacement curve to be positioned that is most similar to the invert deformation curve as the selected displacement curve; A damage unit is configured to use the calculation condition corresponding to the selected displacement curve as the damage condition of the target section invert of the target gypsum-salt tunnel; The extraction unit is further configured to: When there are at least two maximum value points or at least two minimum value points in the invert deformation curve, the maximum value points and the minimum value points are used as the characteristic points; When there is only one maximum point and one minimum point in the deformation curve of the invert, the maximum point, the minimum point and the point with a displacement value of 0 are taken as the characteristic points; The extraction unit is further configured to: Randomly assigning each of the feature points with status information of existence of dissolution cavity, existence of expansion stress or no abnormality; Combining and exhaustively enumerating the state information of all the characteristic points to form a plurality of the calculation conditions; The computing unit is further configured to: Constructing a mechanical model of the target gypsum-salt tunnel, calculating the surrounding rock load of the target gypsum-salt tunnel according to the surrounding rock data and applying the load to the target gypsum-salt tunnel; Obtaining state information corresponding to different characteristic points in the calculation condition; For the feature point whose status information is a dissolution cavity, the surrounding rock load at the bottom of the invert area corresponding to the feature point is set to 0; For the characteristic point whose state information is expansion stress, superimposing the expansion stress on the surrounding rock load at the bottom of the invert area corresponding to the characteristic point; Perform mechanical calculations on the mechanical model with updated invert arch load information.

4. The gypsum-salt stratum tunnel invert damage analysis system according to claim 3 is characterized in that: The comparison unit is further configured to: Calculating the first-order derivative of the displacement curve to be determined to form a slope curve to be determined, and calculating the first-order derivative of the invert deformation curve to form a reference slope curve; The cosine similarity between the pending slope curve and the reference slope curve is calculated, and the pending slope curve with the highest cosine similarity is used as the selected displacement curve.

Images