A method for determining freezing force of tunnels in cold regions and related equipment

Abstract

The present invention discloses a freezing force determination method for a cold region tunnel and related equipment. The method comprises: obtaining measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface; calculating first pressure correction coefficients in different measuring directions according to the measured pressure information and theoretical pressure information; and calculating freezing force distribution information in the corresponding measuring direction according to the first pressure correction coefficients in different measuring directions and the theoretical pressure information.

Description

Technical Field

[0001] The present specification relates to the field of tunnel engineering, and more specifically, the present invention relates to a method for determining freezing force of tunnels in cold regions and related equipment. Background Art

[0002] Tunnels in cold regions are often subject to freeze-thaw cycles and frost heaving, resulting in lining cracks, water leakage and other diseases. Studies have shown that the periodic frost heaving caused by changes in the temperature field is the main cause of various diseases in cold region tunnels. Therefore, studying the frost heaving effect has become a breakthrough in solving the diseases of cold region tunnels. At present, there are many theoretical models and empirical formulas for the study of frost heaving force, but they are all derived based on the assumption that the geology of the tunnel is symmetrically distributed. However, the actual strata in the tunnel are not evenly distributed, and the current theoretical model is not accurate enough for the calculation of the actual tunnel freezing force. Summary of the invention

[0003] A series of simplified concepts are introduced in the Summary of the Invention, which will be further described in detail in the Detailed Description of the Invention. The Summary of the Invention does not mean to attempt to define the key features and essential technical features of the claimed technical solution, nor does it mean to attempt to determine the scope of protection of the claimed technical solution.

[0004] In order to propose a method for determining freezing force of a cold region tunnel that can be used in engineering, in a first aspect, the present invention proposes a method for determining freezing force of a cold region tunnel, the method comprising:

[0005] Obtaining measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0006] Calculate the first pressure correction coefficient in different measuring directions according to the above measured pressure information and theoretical pressure information;

[0007] The freezing force distribution information in the corresponding measuring direction is calculated according to the first pressure correction coefficient in different measuring directions and the theoretical pressure information.

[0008] Optionally, the above method further includes:

[0009] The above theoretical pressure information P is calculated according to the following formula b :

[0010]

[0011] m 1 =(b 2 +a 2 ) / (b 2 -a 2 )

[0012] m 2 =(c 2 +b 2 ) / (c 2 -b 2 )

[0013] p=2b 2 / (c 2 -b 2 )

[0014] q=2c 2 / (c 2 -b 2 )

[0015] e=E 2 -E 1

[0016] Where a is the inner diameter of the tunnel lining in the cold region, b is the inner diameter of the frozen zone, c is the inner diameter of the unfrozen zone, and E 1 , μ 1 are the elastic modulus and Poisson's ratio of the lining, E 2 , μ 2 are the elastic modulus and Poisson's ratio of the freezing circle, respectively.

[0017] Optionally, the above method further includes:

[0018] Obtain the measured pressure information of adjacent target excavated surfaces in the same measuring direction;

[0019] Based on the above measured pressure information in the same direction, a first isobaric surface in the same measuring direction is constructed;

[0020] The freezing force of the remaining cross sections of the adjacent excavated surface in the measuring direction is determined based on the first isobaric surface.

[0021] Optionally, the above method further includes:

[0022] In the case where there are temperature sensors in the same measuring direction on at least two other cross sections between the adjacent target excavated surfaces, constructing a first isotherm according to the temperature data of the temperature sensors;

[0023] Adjusting the first isobaric surface according to the first isotherm;

[0024] The freezing force of the remaining cross sections of the adjacent excavated surface in the above-mentioned measuring direction is determined based on the adjusted first isobaric surface.

[0025] Optionally, adjusting the first isobaric surface according to the first isotherm includes:

[0026] Obtain an isobaric projection of the first isobaric surface on the excavated surface corresponding to the first isotherm;

[0027] Obtain the first derivative values ​​of the above isobar projections in different measurement directions;

[0028] Obtaining second derivative values ​​of the first isotherm in different measurement directions;

[0029] When the difference between the first derivative value and the second derivative value is greater than a preset difference, obtaining an average value of the first derivative value and the second derivative value;

[0030] Adjust the above isobar projection based on the above average value;

[0031] The first isobaric surface is adjusted according to the adjusted isobaric line projection.

[0032] Optionally, the above method further includes:

[0033] Obtain the measured pressure information of the same target excavated surface in adjacent measurement directions;

[0034] Based on the measured pressure information on the same surface, a second isobaric surface is constructed on the same excavation surface;

[0035] The freezing forces in other directions between the adjacent measuring directions are determined according to the second isobaric surface.

[0036] Optionally, the second isobaric surface is obtained by fitting a third-order B-spline curve based on the measured pressure information on the same surface.

[0037] In a second aspect, the present invention further provides a device for determining freezing force of a tunnel in a cold region, comprising:

[0038] An acquisition unit, used to acquire measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is acquired by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0039] A first calculation unit, used for calculating a first pressure correction coefficient in different measurement directions according to the measured pressure information and the theoretical pressure information;

[0040] The second calculation unit is used to calculate the freezing force distribution information in the corresponding measuring direction according to the first pressure correction coefficient and the measured pressure information in different measuring directions.

[0041] In a third aspect, an electronic device comprises: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor is used to implement the steps of the method for determining the freezing force of a tunnel in a cold region as described in any one of the first aspects above when executing the computer program stored in the memory.

[0042] In a fourth aspect, the present invention further proposes a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method for determining the freezing force of a cold-region tunnel according to any one of the above items in the first aspect is implemented.

[0043] In summary, the freezing force calculation method for cold-region tunnels proposed in the present application includes: obtaining the measured pressure information of the target excavated surface in the target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface; calculating the first pressure correction coefficient in different measuring directions according to the measured pressure information and the theoretical pressure information; calculating the freezing force distribution information in the corresponding measuring direction according to the first pressure correction coefficient in different measuring directions and the theoretical pressure information. The freezing force determination method for cold-region tunnels proposed in the embodiment of the present application obtains the first pressure correction coefficient by comparing the measured pressure information of the pressure sensor pre-set on the excavated surface of the target tunnel with the theoretical pressure information of the theoretical calculation model, thereby correcting the theoretical pressure information according to the first pressure correction coefficient, and fitting according to the freezing force information to obtain the freezing force distribution information of the target excavated surface, which can overcome the limitation of the theoretical calculation model assuming that the rock and soil layers are evenly distributed, and the obtained freezing force is more in line with the distribution of the rock and soil characteristics of the actual tunnel.

[0044] The method for determining freezing force of tunnels in cold regions of the present invention, and other advantages, objectives and features of the present invention will be partially reflected through the following description, and will also be partially understood by technicians in this field through research and practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Various other advantages and benefits will become apparent to those of ordinary skill in the art by reading the detailed description of the preferred embodiments below. The accompanying drawings are only for the purpose of illustrating the preferred embodiments and are not to be considered as limiting the present specification. Also, the same reference symbols are used throughout the accompanying drawings to represent the same components. In the accompanying drawings:

[0046] Figure 1 A schematic flow chart of a method for determining freezing force of a tunnel in a cold region provided in an embodiment of the present application;

[0047] Figure 2 A schematic diagram of the principle of a tunnel freezing force calculation model provided in an embodiment of the present application;

[0048] Figure 3 A schematic diagram of the structure of a freezing force determination device for a cold region tunnel provided in an embodiment of the present application;

[0049] Figure 4 A schematic diagram of the structure of an electronic device for determining freezing force in cold-region tunnels provided in an embodiment of the present application. DETAILED DESCRIPTION

[0050] The method for determining the freezing force of a tunnel in a cold region proposed in the embodiment of the present application obtains a first pressure correction coefficient by comparing the measured pressure information of a pressure sensor pre-set on the excavated surface of the target tunnel with the theoretical pressure information of the theoretical calculation model, thereby correcting the theoretical pressure information according to the first pressure correction coefficient, and fitting the freezing force information to obtain the freezing force distribution information of the target excavated surface. This can overcome the limitation of the theoretical calculation model that assumes that the rock and soil layers are uniformly distributed, and the resulting freezing force is more in line with the actual distribution of rock and soil characteristics of the tunnel.

[0051] The terms "first", "second", "third", "fourth", etc. (if any) in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments described herein can be implemented in an order other than that illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices. The technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all of the embodiments.

[0052] See also Figure 1 , which is a flow chart of a method for determining freezing force of a tunnel in a cold region provided in an embodiment of the present application, which may specifically include:

[0053] S110, obtaining measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0054] Exemplarily, the measured pressure information of the excavated surface is obtained by a pressure sensor, and a hole is drilled on the normal line of the tunnel contour line, and a plurality of pressure sensors separated by a preset distance are placed in the borehole, and cement slurry is injected into the borehole after the sensors are arranged. After the cement slurry is solidified, the pressure sensor is consolidated with the rock and soil, and the first pressure correction coefficient in different measurement directions is calculated according to the above measured pressure information and theoretical pressure information through S120;

[0055] Exemplarily, each measurement direction is the drilling direction in which the pressure sensor is installed. The drilling direction is generally consistent with the normal direction of the outer contour of the excavated surface. The theoretical pressure information is the theoretical pressure information calculated by the existing calculation model. The calculation model may include: the frost heave model of the water-containing weathering layer, the frost heave model of the water accumulation behind the lining, the overall frost heave model of the freeze-thaw circle, etc. The measured pressure information and the theoretical pressure information are used to calculate the first pressure correction coefficient, which is to compensate for the deviation between the theoretical pressure information and the measured pressure information.

[0056] S130, calculating freezing force distribution information in corresponding measuring directions according to the first pressure correction coefficient in different measuring directions and the theoretical pressure information.

[0057] Exemplarily, after obtaining the first pressure correction coefficient, the freezing force information in different measuring directions can be obtained according to the product of the correction coefficient in different directions and the theoretical pressure information, and the freezing force information in different directions is fitted along the outer contour direction of the tunnel to obtain the freezing force distribution information.

[0058] In summary, the method for determining the freezing force of a tunnel in a cold region proposed in the embodiment of the present application obtains a first pressure correction coefficient by comparing the measured pressure information of a pressure sensor pre-set on the excavated surface of the target tunnel with the theoretical pressure information of the theoretical calculation model, thereby correcting the theoretical pressure information according to the first pressure correction coefficient, and fitting the freezing force information to obtain the freezing force distribution information of the target excavated surface. This can overcome the limitation of the theoretical calculation model that assumes that the rock and soil layers are uniformly distributed, and the resulting freezing force is more in line with the actual distribution of rock and soil characteristics of the tunnel.

[0059] In some examples, the method further includes:

[0060] The above theoretical pressure information P is calculated according to the following formula b :

[0061]

[0062] m 1 =(b 2 +a 2 ) / (b 2 -a 2 )

[0063] m 2 =(c 2 +b 2 ) / (c 2 -b 2 )

[0064] p=2b 2 / (c 2 -b 2 )

[0065] q=2c 2 / (c 2 -b 2 )

[0066] e=E 2 -E 1

[0067] Where a is the inner diameter of the tunnel lining in the cold region, b is the inner diameter of the frozen zone, c is the inner diameter of the unfrozen zone, and E 1 , μ 1 are the elastic modulus and Poisson's ratio of the lining, E 2 , μ 2 are the elastic modulus and Poisson's ratio of the freezing circle, respectively.

[0068] For example, Figure 2 The figure shows a frost heave model based on the overall freeze-thaw circle. By calculating the frost heave deformation caused by the freezing of the surrounding rock of the frozen circle, the elastic theory is used to calculate the frost heave force on the lining caused by the frost heave deformation of the frozen circle according to the displacement continuity condition of the multi-layer thick-walled cylinder. The inner diameter of the frozen circle remains unchanged, and the frost heave causes the outer diameter of the frozen circle to increase. Based on this assumption, the above formula can be obtained, and the theoretical pressure information P can be calculated according to the above formula. b .

[0069] In some examples, the method further includes:

[0070] Obtain the measured pressure information of adjacent target excavated surfaces in the same measuring direction;

[0071] Based on the above measured pressure information in the same direction, a first isobaric surface in the same measuring direction is constructed;

[0072] The freezing force of the remaining cross sections between the adjacent excavated surfaces in the measuring direction is determined based on the first isobaric surface.

[0073] For example, the target excavated surfaces where pressure sensors are installed in the tunnel are spaced apart and do not exist continuously at a small distance, and the geotechnical structures between adjacent target excavated surfaces should be roughly similar, with only some changes in thickness or direction. Therefore, the measured pressure information between two adjacent target excavated surfaces is comprehensively considered, and isobaric surfaces are established in the same measuring direction based on the measured pressure information. The general direction of the isobaric surface is basically similar to the structural distribution of the geotechnical surface. Then, based on the intersection of the isobaric surface and the remaining cross-sections between the adjacent excavated surfaces, the pressure of the remaining cross-sections at different depths can be obtained, and the correction coefficient of the cross-section in the measuring direction can be obtained based on the pressure and its pressure position, and the freezing force distribution of the cross-section in the direction can be obtained based on the correction coefficient and the theoretical pressure information.

[0074] In some examples, the method further includes:

[0075] In the case where there are temperature sensors in the same measuring direction on at least two other cross sections between the adjacent target excavated surfaces, constructing a first isotherm according to the temperature data of the temperature sensors;

[0076] Adjusting the first isobaric surface according to the first isotherm;

[0077] The freezing forces of the remaining cross sections between the adjacent excavated surfaces in the above-mentioned measuring direction are determined according to the adjusted first isobaric surface.

[0078] Exemplarily, in the case where there may be temperature sensors on at least two other sections in the same measuring direction between every two adjacent target excavated surfaces, the temperature distribution information is also correlated with the structural distribution of the rock and soil. If the temperature isotherm is constructed based on the temperature distribution information, and if the shape change pattern of the first isotherm is basically similar to that of the first isobaric surface, it is considered that the temperature and pressure characterize the structural characteristics of the rock and soil in a basically consistent manner, and the first isobaric surface does not need to be corrected; if there is a large anomaly in the change pattern of the first isotherm and the first isobaric surface, it is necessary to adjust the first isotherm according to the first isotherm to correct the calculation result of the freezing force.

[0079] In some examples, adjusting the first isobaric surface according to the first isotherm includes:

[0080] Obtain an isobaric projection of the first isobaric surface on the excavated surface corresponding to the first isotherm;

[0081] Obtain the first derivative values ​​of the above isobar projections in different measurement directions;

[0082] Obtaining second derivative values ​​of the first isotherm in different measurement directions;

[0083] When the difference between the first derivative value and the second derivative value is greater than a preset difference, obtaining an average value of the first derivative value and the second derivative value;

[0084] Adjust the above isobar projection based on the above average value;

[0085] The first isobaric surface is adjusted according to the adjusted isobaric line projection.

[0086] Exemplarily, the isobaric projection of the first isobaric surface on the excavated surface provided with the temperature sensor is obtained, and the similarity between the shape contours of the isobaric projection and the first isotherm is compared. The similarity between the shape contours is determined based on the derivative values ​​of the isobaric projection and the first isotherm in different measuring directions. If the first derivative value and the second derivative value in the same direction are greater than the preset difference, that is, the difference is large, it is considered that the isotherm and the isobaric contour in this direction have a large deviation. At this time, the average value between the isotherm and the isobar is calculated, the shape of the isobaric projection is adjusted according to the average value, and the first isobaric surface is adjusted according to the adjusted isobaric projection and refitted. The isobaric surface is adjusted by combining the temperature information and the pressure information to make the estimation result more accurate.

[0087] In some examples, the method further includes:

[0088] Obtain the measured pressure information of the same target excavated surface in adjacent measurement directions;

[0089] Based on the measured pressure information on the same surface, a second isobaric surface is constructed on the same excavation surface;

[0090] The freezing forces in other directions between the adjacent measuring directions are determined according to the second isobaric surface.

[0091] For example, the measurement directions of the excavated surface of the same target are limited. The continuity of the stratum changes in the local area can be taken into account by fitting the isobaric surface, and the freezing force in other directions between two adjacent measurement directions can be estimated based on the two adjacent measurement directions.

[0092] In some examples, the second isobaric surface is obtained by fitting a third-order B-spline curve based on the measured pressure information on the same surface.

[0093] For example, the third-order B-spline curve is a smooth curve representation method, which can generate a smooth curve through data points. Since the third-order B-spline curve is a locally controlled curve representation method, it can accurately fit the data points while maintaining the smoothness of the curve. The curve generated by the third-order B-spline curve has the characteristics of being able to conform to the distribution of the strata.

[0094] See also Figure 3 , an embodiment of the device for determining freezing force of a tunnel in a cold region in the embodiment of the present application may include:

[0095] An acquisition unit 21 is used to acquire measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is acquired by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0096] A first calculation unit 22, used to calculate a first pressure correction coefficient in different measurement directions according to the measured pressure information and the theoretical pressure information;

[0097] The second calculation unit 23 is used to calculate the freezing force distribution information in the corresponding measuring direction according to the first pressure correction coefficient and the measured pressure information in different measuring directions.

[0098] like Figure 4 As shown, an embodiment of the present application also provides an electronic device 300, including a memory 310, a processor 320, and a computer program 311 stored in the memory 320 and executable on the processor. When the processor 320 executes the computer program 311, the steps of any of the above-mentioned methods for determining the freezing force of a cold zone tunnel are implemented.

[0099] Since the electronic device introduced in this embodiment is a device used to implement a method for determining the freezing force of a cold zone tunnel in the embodiment of the present application, based on the method introduced in the embodiment of the present application, technical personnel in this field can understand the specific implementation mode of the electronic device of this embodiment and its various variations. Therefore, how the electronic device implements the method in the embodiment of the present application is not introduced in detail here. As long as the equipment used by technical personnel in this field to implement the method in the embodiment of the present application is within the scope of protection of this application.

[0100] In the specific implementation process, when the computer program 311 is executed by the processor, it can achieve Figure 1 Any implementation method in the corresponding embodiment includes:

[0101] Obtaining measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0102] Calculate the first pressure correction coefficient in different measuring directions according to the above measured pressure information and theoretical pressure information;

[0103] The freezing force distribution information in the corresponding measuring direction is calculated according to the first pressure correction coefficient and the measured pressure information in different measuring directions.

[0104] In a possible implementation, the above method further includes:

[0105] The above theoretical pressure information P is calculated according to the following formula b :

[0106]

[0107] m 1 =(b 2 +a 2 ) / (b2 -a 2 )

[0108] m 2 =(c 2 +b 2 ) / (c 2 -b 2 )

[0109] p=2b 2 / (c 2 -b 2 )

[0110] q=2c 2 / (c 2 -b 2 )

[0111] e=E 2 -E 1

[0112] Where a is the inner diameter of the tunnel lining in the cold region, b is the inner diameter of the frozen zone, c is the inner diameter of the unfrozen zone, and E 1 , μ 1 are the elastic modulus and Poisson's ratio of the lining, E 2 , μ 2 are the elastic modulus and Poisson's ratio of the freezing circle, respectively.

[0113] In a possible implementation, the above method further includes:

[0114] Obtain the measured pressure information of adjacent target excavated surfaces in the same measuring direction;

[0115] Based on the above measured pressure information in the same direction, a first isobaric surface in the same measuring direction is constructed;

[0116] The freezing force of the remaining cross sections between the adjacent excavated surfaces in the measuring direction is determined based on the first isobaric surface.

[0117] In a possible implementation, the above method further includes:

[0118] In the case where there are temperature sensors in the same measuring direction on at least two other cross sections between the adjacent target excavated surfaces, constructing a first isotherm according to the temperature data of the temperature sensors;

[0119] Adjusting the first isobaric surface according to the first isotherm;

[0120] The freezing forces of the remaining cross sections between the adjacent excavated surfaces in the above-mentioned measuring direction are determined according to the adjusted first isobaric surface.

[0121] In a possible implementation manner, adjusting the first isobaric surface according to the first isotherm includes:

[0122] Obtain an isobaric projection of the first isobaric surface on the excavated surface corresponding to the first isotherm;

[0123] Obtain the first derivative values ​​of the above isobar projections in different measurement directions;

[0124] Obtaining second derivative values ​​of the first isotherm in different measurement directions;

[0125] When the difference between the first derivative value and the second derivative value is greater than a preset difference, obtaining an average value of the first derivative value and the second derivative value;

[0126] Adjust the above isobar projection based on the above average value;

[0127] The first isobaric surface is adjusted according to the adjusted isobaric line projection.

[0128] In a possible implementation, the above method further includes:

[0129] Obtain the measured pressure information of the same target excavated surface in adjacent measurement directions;

[0130] Based on the measured pressure information on the same surface, a second isobaric surface is constructed on the same excavation surface;

[0131] The freezing forces in other directions between the adjacent measuring directions are determined according to the second isobaric surface.

[0132] In a possible implementation manner, the second isobaric surface is obtained by fitting a third-order B-spline curve based on the measured pressure information on the same surface.

[0133] It should be noted that in the above embodiments, the description of each embodiment has its own emphasis, and for parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0134] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0135] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded computer, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0136] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0137] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0138] The present application also provides a computer program product, which includes computer software instructions. When the computer software instructions are executed on a processing device, the processing device is caused to execute the following Figure 1 The process of determining the freezing force of a cold region tunnel in the corresponding embodiment includes:

[0139] Obtaining measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface;

[0140] Calculate the first pressure correction coefficient in different measuring directions according to the above measured pressure information and theoretical pressure information;

[0141] The freezing force distribution information in the corresponding measuring direction is calculated according to the first pressure correction coefficient and the measured pressure information in different measuring directions.

[0142] In a possible implementation, the above method further includes:

[0143] The above theoretical pressure information P is calculated according to the following formula b :

[0144]

[0145] m 1 =(b 2 +a 2 ) / (b 2 -a 2 )

[0146] m 2 =(c 2 +b 2 ) / (c 2 -b 2 )

[0147] p=2b 2 / (c 2 -b 2 )

[0148] q=2c 2 / (c 2 -b 2 )

[0149] e=E 2 -E 1

[0150] Where a is the inner diameter of the tunnel lining in the cold region, b is the inner diameter of the frozen zone, c is the inner diameter of the unfrozen zone, and E 1 , μ 1 are the elastic modulus and Poisson's ratio of the lining, E 2 , μ 2 are the elastic modulus and Poisson's ratio of the freezing circle, respectively.

[0151] In a possible implementation, the above method further includes:

[0152] Obtain the measured pressure information of adjacent target excavated surfaces in the same measuring direction;

[0153] Based on the above measured pressure information in the same direction, a first isobaric surface in the same measuring direction is constructed;

[0154] The freezing force of the remaining cross sections between the adjacent excavated surfaces in the measuring direction is determined based on the first isobaric surface.

[0155] In a possible implementation, the above method further includes:

[0156] In the case where there are temperature sensors in the same measuring direction on at least two other cross sections between the adjacent target excavated surfaces, constructing a first isotherm according to the temperature data of the temperature sensors;

[0157] Adjusting the first isobaric surface according to the first isotherm;

[0158] The freezing forces of the remaining cross sections between the adjacent excavated surfaces in the above-mentioned measuring direction are determined according to the adjusted first isobaric surface.

[0159] In a possible implementation manner, adjusting the first isobaric surface according to the first isotherm includes:

[0160] Obtain an isobaric projection of the first isobaric surface on the excavated surface corresponding to the first isotherm;

[0161] Obtain the first derivative values ​​of the above isobar projections in different measurement directions;

[0162] Obtaining second derivative values ​​of the first isotherm in different measurement directions;

[0163] When the difference between the first derivative value and the second derivative value is greater than a preset difference, obtaining an average value of the first derivative value and the second derivative value;

[0164] Adjust the above isobar projection based on the above average value;

[0165] The first isobaric surface is adjusted according to the adjusted isobaric line projection.

[0166] In a possible implementation, the above method further includes:

[0167] Obtain the measured pressure information of the same target excavated surface in adjacent measurement directions;

[0168] Based on the measured pressure information on the same surface, a second isobaric surface is constructed on the same excavation surface;

[0169] The freezing forces in other directions between the adjacent measuring directions are determined according to the second isobaric surface.

[0170] In a possible implementation manner, the second isobaric surface is obtained by fitting a third-order B-spline curve based on the measured pressure information on the same surface.

[0171] The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on the computer, the process or function according to the embodiment of the present application is generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions can be transmitted from a website site, a computer, a server or a data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (digital subscriber line, DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) mode to another website site, computer, server or data center. The computer-readable storage medium can be any available medium that a computer can store or a data storage device such as a server or a data center that includes one or more available media integration. The available medium can be a magnetic medium, (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid state drive (SSD)), etc.

[0172] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0173] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that 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, which can be electrical, mechanical or other forms.

[0174] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0175] In addition, each functional unit in each embodiment of the present application 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-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0176] 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 application 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 to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the various embodiments of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), disk or optical disk and other media that can store program codes.

[0177] The above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some of the technical features therein by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for determining freezing force of tunnels in cold regions, characterized in that: include: Acquire measured pressure information of a target excavated surface in a target tunnel, wherein the measured pressure information is obtained by installing a preset number of pressure sensors at a preset depth of the excavated surface; Calculating first pressure correction coefficients in different measuring directions according to the measured pressure information and the theoretical pressure information; Calculate the freezing force distribution information in the corresponding measuring direction according to the first pressure correction coefficient in different measuring directions and the theoretical pressure information; Also includes: Obtain the measured pressure information of adjacent target excavated surfaces in the same measuring direction; Constructing a first isobaric surface in the same measuring direction based on the measured pressure information in the same direction; In the case where there are temperature sensors in the same measuring direction on at least two other cross sections between the adjacent target excavated surfaces, constructing a first isotherm according to temperature data of the temperature sensors; adjusting the first isobaric surface according to the first isotherm; The freezing forces of the remaining cross sections between the adjacent excavated surfaces in the measuring direction are determined according to the adjusted first isobaric surface.

2. The method according to claim 1, characterized in that The adjusting the first isobaric surface according to the first isotherm comprises: Obtain an isobaric projection of the first isobaric surface on the excavated surface corresponding to the first isotherm; Obtaining the first derivative values ​​of the isobar projection in different measurement directions; Obtaining second derivative values ​​of the first isotherm in different measurement directions; When the difference between the first derivative value and the second derivative value is greater than a preset difference, obtaining an average value of the first derivative value and the second derivative value; adjusting the isobar projection based on the average value; The first isobaric surface is adjusted according to the adjusted isobaric line projection.

3. The method according to claim 1, characterized in that Also includes: Obtain the measured pressure information of the same target excavated surface in adjacent measurement directions; Constructing a second isobaric surface on the same excavation surface based on the measured pressure information on the same surface; The freezing forces in other directions between the adjacent measuring directions are determined according to the second isobaric surface.

4. The method according to claim 3, characterized in that The second isobaric surface is obtained by fitting a third-order B-spline curve based on the measured pressure information on the same surface.

5. An electronic device comprising: A memory and a processor, characterized in that the processor is used to implement the steps of the method for determining the freezing force of a cold region tunnel as described in any one of claims 1 to 4 when executing the computer program stored in the memory.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for determining freezing force of a cold region tunnel as described in any one of claims 1 to 4 is implemented.

Images