A method for identifying the tension of pipe structures based on tension coefficient matrix correction

Abstract

This invention provides a pipeline structure tension identification method based on tension coefficient matrix correction. This method can identify the tension distribution during pipeline service, providing data support for real-time monitoring and evaluation of pipeline service status, thereby improving pipeline service safety and stability. The method includes: constructing a pipeline coordinate system; determining pipeline tension function parameters based on the relative positions between the pipeline's constrained ends; and then constructing a tension function model. Based on this, the tension identification equation for the pipeline under concentrated force is derived. A gravity constraint iterative method is introduced to solve the tension identification equation when the pipeline is not under concentrated force. Simultaneously, by extracting strain sensing information from fiber optic sensors integrated into the pipeline, a tension feature correction coefficient matrix is ​​constructed to significantly improve the tension identification accuracy.

Description

Technical Field

[0001] This invention belongs to the field of pipeline structure monitoring technology, and particularly relates to a pipeline structure tension identification method based on tension coefficient matrix correction. Background Technology

[0002] Modern industry, urban lifelines, and the aviation sector (such as aerial refueling hoses) heavily rely on various pipeline systems. These pipeline structures often operate under complex load environments, and pipeline tension, as a core safety and performance indicator, directly affects their buckling resistance, fatigue life, overall stability, and reliability. However, existing detection technologies primarily target pipeline structural defects and are ill-suited for long-distance, continuous, and real-time monitoring. Therefore, developing high-precision, intelligent tension identification methods is urgently needed. It provides a scientific basis for assessing the condition of pipelines in extreme environments, dynamic service conditions, and aging pipelines, optimizing maintenance / operation decisions, extending asset life, and improving mission success rates. It also meets the strategic requirements of increasingly stringent safety regulations and achieving lean pipeline integrity management. This research is of irreplaceable importance for building quality and safety, improving economic efficiency, and enhancing national defense capabilities.

[0003] In recent years, strain sensing technology based on fiber optic sensors has been increasingly applied to the field of structural health monitoring. Fiber optic sensors offer advantages such as high precision, resistance to electromagnetic interference, lightweight, and corrosion resistance, enabling real-time monitoring of strain distribution in structures. However, a mature, universally applicable, and effective method for tension identification that can utilize the abundant strain data acquired by fiber optic sensors is currently lacking, making it difficult to establish a mapping relationship between distributed strain and tension distribution in pipeline structures under complex operating conditions. Summary of the Invention

[0004] Objective of the Invention: The technical problem to be solved by this invention is to address the shortcomings of existing technologies by providing a method for identifying the tension of pipe structures based on tension coefficient matrix correction, comprising the following steps:

[0005] Step 1: Establish a rectangular coordinate system with the lower end of the pipe as the origin to determine the pipe length. l Height difference between the two constrained ends of the pipeline h 0, Angle between the lines connecting the two ends of the pipe θ The gravitational load p on the pipeline;

[0006] Step 2: Based on the tension function theory, derive the tension identification equation corresponding to the pipe being subjected to a vertically downward concentrated force;

[0007] Step 3: The gravity constraint iterative method is used to numerically solve the key parameters of the tension identification equation when the pipeline is not subjected to a vertically downward concentrated force.

[0008] Step 4: Attach fiber optic sensors to the pipe surface to acquire strain distribution data of the pipe surface during service. Based on the measured original strain data matrix and the first derivative matrix, construct a tension characteristic correction coefficient matrix. Substitute the tension characteristic correction coefficient matrix into the tension identification equation to correct the tension identification error caused by the lack of detailed material coefficient characteristics of the pipe.

[0009] In step 1, define the pipes. AB The horizontal length is l The gravitational load p on the pipe is calculated by the pipe density.

[0010] Step 2 includes:

[0011] Step 2.1, take any infinitesimal segment on the pipeline. ds The pipe is set to be under tension. T The horizontal component of the force is H The vertical component of the force is V Then, the static equilibrium equation of the pipeline during service can be obtained from static analysis:

[0012] (1),

[0013] in d y is the differential symbol; y represents the infinitesimal segment of the pipe. Y Axis coordinates; x For pipeline micro-segments X Axis coordinates;

[0014] Solve equation (1) by performing a second integral, and based on the boundary conditions: the x-coordinate of the pipeline starting point. x 1=0, ordinate y 1=0; x-coordinate of the pipeline endpoint x 2=l, ordinate y 2=l, thus obtaining the linear equation of the pipe in its free state:

[0015] (2),

[0016] in , Let be an arbitrary constant generated during the integration process, determined by equation (3):

[0017] (3),

[0018] in arsinh Represents the inverse hyperbolic sine function;

[0019] Step 2.2: Let point C be the lowest point of the pipeline. Differentiate equation (3) and take the extreme point to obtain the x-coordinate of the lowest point C. for:

[0020] (4),

[0021] Substituting equation (4) into equation (2), we obtain the pipeline alignment function under the condition of the known minimum point. y ( x ):

[0022] (5),

[0023] Assume any position of the pipeline Q ( x , y Point tension is T The vertical component of the tension is V The pipe tension identification equation is obtained from static analysis:

[0024] (6),

[0025] (7),

[0026] According to equation (7), the pipe is at its lowest point The tension value is minimum at the lowest point. The tension gradually increases towards both ends, reaching its maximum at point B, the higher end. Therefore, the maximum tension in the pipe is... and minimum value for:

[0027] (8),

[0028] in The vertical component of the tension at point B;

[0029] Step 2.3: When a pipe is subjected to a concentrated force acting vertically downward, the pipe shape can be classified into two types. One type of concentrated force is applied to the center of the pipe, which is called pipe type I. The other type of concentrated force is applied near both ends of the pipe, which is called pipe type II.

[0030] In step 2.3, when the pipe shape is linear type I, combining equations (2), (6), and (7) yields:

[0031] (9),

[0032] When the pipeline configuration is linear type II, combining equations (2), (6), and (7) yields:

[0033] (10)

[0034] In the formula:

[0035] (11),

[0036] in Location where concentrated force is applied X Axis coordinates for Y Axis coordinates The concentrated force load on the pipeline;

[0037] We obtain the results by numerically solving equations (9) and (10). The value will Substituting into equations (6) and (7), we can obtain the tension distribution characteristics of the pipeline when it is subjected to a vertically downward concentrated force.

[0038] Step 3 includes:

[0039] Step 3.1, based on the known x-coordinate of the lowest point of the pipeline y-axis To obtain information about the parameters The constraint equations are:

[0040] (12)

[0041] Step 3.2: Numerical solution is performed using the gravity constraint iterative method to construct the iterative function. :

[0042] (13)

[0043] Step 3.3: Solve for the parameters using the gravity constraint iterative method. When calculating the derivative of the iterative function... :

[0044] (14)

[0045] in:

[0046] (15)

[0047] Step 3.4, the gravity constraint iterative method format is as follows:

[0048] (16)

[0049] in, Represents the parameters obtained in the nth iteration. ; The convergence factor; Indicates parameter value The function, that is, the function M and N The convergence condition is: , The set convergence precision;

[0050] Step 3.5: Verify the calculation results using gravity constraint equations, assuming the pipeline... A The vertical component of the tension is V A ,pipeline B The vertical component of the tension is V B ,but:

[0051] (17)

[0052] According to the gravity constraint equation ps = V A + V B The calculation results must meet the following requirements:

[0053] (18)

[0054] in s Indicates the arc length of the pipe.

[0055] Step 4 includes:

[0056] Step 4.1: Arrange three optical fibers along the circumferential direction at 0°, 120°, and 240° on the surface of the pipe. Bragg The strain sensing path of the grating sensor; composed of optical fiber Bragg Actual tensile and compressive strain of the pipeline collected by the grating sensor ;

[0057] Step 4.2, the strains obtained at each discrete measurement point are as follows: Then the strain measurement matrix E is defined as:

[0058] (19)

[0059] in Indicates the first n The strain obtained from the measurement of discrete measuring points;

[0060] The first derivative matrix of the measured strain matrix is ​​calculated using the central difference method. :

[0061] (20)

[0062] in For the first i Strain values ​​at individual fiber optic measuring points; For the first i -1 strain value at fiber optic measuring point; For the first i fiber optic measurement points X Axis coordinate values; iValues ​​range from 1 to n ;

[0063] The first derivative matrix of the measured strain matrix was calculated. for:

[0064] (twenty one),

[0065] Step 4.3: Based on the measured strain matrix and its first derivative obtained from the fiber optic sensing system, the tension characteristic correction coefficient matrix is ​​obtained. K :

[0066] (twenty two),

[0067] Step 4.4: Establish the pipeline tension identification equation based on the tension correction coefficient matrix.

[0068] In step 4.1, the actual tensile and compressive strain of the pipeline The calculation formula is:

[0069] (twenty three),

[0070] in, , , Fiber optic cables of 0°, 120°, and 240° respectively. Bragg The strain is measured by a grating sensor.

[0071] In step 4.1, the pipe tension identification equation based on the tension correction coefficient matrix is:

[0072] (twenty four),

[0073] (25),

[0074] in Y ( x ) represents the pipe profile function based on the tension correction coefficient matrix; This represents the tension identification function for pipe structures based on the tension correction coefficient matrix.

[0075] The present invention also provides an electronic device, including a processor and a memory, the memory storing program code that, when executed by the processor, causes the processor to perform the steps of the method.

[0076] The present invention also provides a storage medium storing a computer program or instructions that, when the computer program or instructions are run on a computer, execute the steps of the method described.

[0077] This invention constructs a tension function model for the pipeline by determining its coordinate system and considering the spatial relationship between its two constrained ends. Furthermore, it proposes a numerical iteration method incorporating gravity constraints to solve for the tension distribution of the pipeline under its own weight, thus resolving the difficulty in solving the equations due to the transcendental nature of the tension function. Simultaneously, by extracting strain information from the fiber optic sensor integrated with the pipeline, a tension feature correction coefficient matrix is ​​constructed, significantly improving the accuracy of tension identification.

[0078] This invention aims to provide data support for real-time monitoring and evaluation of pipeline service status, thereby improving the safety and stability of pipeline service.

[0079] Beneficial Effects: This invention integrates fiber optic sensors with an innovative tension recognition model, significantly improving pipeline monitoring. It offers the following beneficial effects:

[0080] By associating the pipeline's shape function with the load it experiences, tension identification equations for the pipeline under different stress conditions are derived, overcoming the limitation that the pipeline tension function only includes the pipeline's gravity parameters. A numerical iteration method with gravity constraints is proposed to solve for the tension distribution of the pipeline under self-weight conditions, solving the problem of difficulty in solving the equations caused by the tension function being a transcendental equation. At the same time, by extracting strain information from the fiber optic sensor integrated with the pipeline, a tension feature correction coefficient matrix is ​​constructed, significantly improving the accuracy of tension identification. Attached Figure Description

[0081] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments, and the advantages of the present invention in the above and / or other aspects will become clearer.

[0082] Figure 1 This is a flowchart of the method of the present invention.

[0083] Figure 2 This is a schematic diagram of the simulation constraint application obtained from an embodiment of the present invention.

[0084] Figure 3 This is a comparison diagram of the tension distribution identified and simulated by the method of this invention.

[0085] Figure 4 This is a schematic diagram of the relative tension error identified and simulated by the method of the present invention. Detailed Implementation

[0086] like Figure 1 As shown, this embodiment provides a method for identifying the tension of a pipe structure based on tension coefficient matrix correction, including the following steps:

[0087] Step 1: Establish a rectangular coordinate system with the lower end of the pipe as the origin, and determine the horizontal length of the pipe. lHeight difference between the two constrained ends of the pipeline h 0, Angle between the lines connecting the two ends of the pipe θ The gravitational load p on the pipeline;

[0088] Step 2: Based on the tension function theory, derive the tension identification equation when the pipe is subjected to a vertically downward concentrated force;

[0089] Step 3: The key parameters of the tension identification equation when the pipeline is not subjected to a vertically downward concentrated force are numerically solved using the gravity constraint iterative method;

[0090] Step 4: Attach fiber optic sensors to the pipe surface to acquire strain distribution data during service. Based on the measured original strain data matrix and its first derivative matrix, construct a tension characteristic correction coefficient matrix. Substitute the tension characteristic correction coefficient matrix into the tension identification equation to correct the tension identification error caused by the lack of detailed material coefficient characteristics of the pipe.

[0091] Step 2 includes:

[0092] Step 2.1, take any infinitesimal segment on the pipeline. ds Assume the pipe is under tension. T The horizontal component of the force is H The vertical component of the force is V Then, the static equilibrium equation of the pipeline during the process can be obtained from static analysis:

[0093] (1),

[0094] In the formula, d y is the differential symbol; y represents the infinitesimal segment of the pipe. Y Axis coordinates; x For pipeline micro-segments X Axis coordinates.

[0095] Solve equation (1) by performing a second integral, and based on the boundary conditions: the x-coordinate of the pipeline starting point. x 1=0, ordinate y 1=0; x-coordinate of the pipeline endpoint x 2=l, ordinate y 2=l, thus obtaining the linear equation of the pipe in its free state:

[0096] (2),

[0097] in , Let be an arbitrary constant generated during the integration process, determined by equation (3):

[0098] (3),

[0099] in arsinh Represents the inverse hyperbolic sine function;

[0100] Step 2.2: Let point C be the lowest point of the pipeline. Differentiate equation (3) and take the extreme point to obtain the x-coordinate of the lowest point C. for:

[0101] (4),

[0102] Substituting equation (4) into equation (2), we obtain the pipeline alignment function under the condition of the known minimum point. y ( x ):

[0103] (5),

[0104] Assume any position of the pipeline Q ( x , y Point tension is T The vertical component of the tension is V The pipe tension identification equation is obtained from static analysis:

[0105] (6),

[0106] (7),

[0107] According to equation (7), the pipe is at its lowest point The tension value is minimum at the lowest point. The tension gradually increases towards both ends, reaching its maximum at point B, the higher end. Therefore, the maximum tension in the pipe is... and minimum value for:

[0108] (8),

[0109] in The vertical component of the tension at point B;

[0110] Step 2.3: When a pipe is subjected to a concentrated force acting vertically downward, the pipe shape can be classified into two types. One type of concentrated force is applied to the center of the pipe, which is called pipe type I. The other type of concentrated force is applied near both ends of the pipe, which is called pipe type II.

[0111] In step 2.3, when the pipe shape is linear type I, combining equations (2), (6), and (7) yields:

[0112] (9),

[0113] When the pipeline configuration is linear type II, combining equations (2), (6), and (7) yields:

[0114] (10)

[0115] In the formula:

[0116] (11),

[0117] in Location where concentrated force is applied X Axis coordinates for Y Axis coordinates The concentrated force load on the pipeline;

[0118] We obtain the results by numerically solving equations (9) and (10). The value will Substituting into equations (6) and (7), we can obtain the tension distribution characteristics of the pipeline when it is subjected to a vertically downward concentrated force.

[0119] Step 3 includes:

[0120] Step 3.1, based on the known x-coordinate of the lowest point of the pipeline y-axis To obtain information about the parameters The constraint equations are:

[0121] (12)

[0122] Step 3.2: Numerical solution is performed using the gravity constraint iterative method to construct the iterative function. :

[0123] (13)

[0124] Step 3.3: Solve for the parameters using the gravity constraint iterative method. When calculating the derivative of the iterative function... :

[0125] (14)

[0126] in:

[0127] (15)

[0128] Step 3.4, the gravity constraint iterative method format is as follows:

[0129] (16)

[0130] in, Represents the parameters obtained in the nth iteration. ; The convergence factor is used to ensure convergence. , It can be halved based on the previous iteration. Indicates parameter value The function, that is, the function M and N The convergence condition is: , The set convergence precision;

[0131] Step 3.5: Verify the calculation results using gravity constraint equations, assuming the pipeline... A The vertical component of the tension is V A ,pipeline B The vertical component of the tension is V B ,but:

[0132] (17)

[0133] According to the gravity constraint equation ps = V A + V B The calculation results must meet the following requirements:

[0134] (18)

[0135] in s Indicates the arc length of the pipe.

[0136] Step 4 includes:

[0137] Step 4.1: Arrange three optical fibers along the circumferential direction at 0°, 120°, and 240° on the surface of the pipe. Bragg The strain sensing path of the grating sensor; composed of optical fiber Bragg Actual tensile and compressive strain of the pipeline collected by the grating sensor ;

[0138] Step 4.2, the strains obtained at each discrete measurement point are as follows: Then the strain measurement matrix E is defined as:

[0139] (19)

[0140] in Indicates the first n The strain obtained from the measurement of discrete measuring points;

[0141] The first derivative matrix of the measured strain matrix is ​​calculated using the central difference method. :

[0142] (20)

[0143] in For the first i Strain values ​​at individual fiber optic measuring points; For the first i -1 strain value at fiber optic measuring point; For the first i fiber optic measurement points X Axis coordinate values; i Values ​​range from 1 to n ;

[0144] The first derivative matrix of the measured strain matrix was calculated. for:

[0145] (twenty one),

[0146] Step 4.3: Based on the measured strain matrix and its first derivative obtained from the fiber optic sensing system, the tension characteristic correction coefficient matrix is ​​obtained. K :

[0147] (twenty two),

[0148] Step 4.4: Establish the pipeline tension identification equation based on the tension correction coefficient matrix.

[0149] In step 4.1, the actual tensile and compressive strain of the pipeline The calculation formula is:

[0150] (twenty three),

[0151] in, , , Fiber optic cables of 0°, 120°, and 240° respectively. Bragg The strain is measured by a grating sensor.

[0152] In step 4.1, the pipe tension identification equation based on the tension correction coefficient matrix is:

[0153] (twenty four),

[0154] (25),

[0155] in Y ( x ) represents the pipe profile function based on the tension correction coefficient matrix; This represents the tension identification function for pipe structures based on the tension correction coefficient matrix.

[0156] In one specific embodiment of the present invention, a method for identifying the tension of a pipe structure based on tension coefficient matrix correction is provided, comprising:

[0157] First, based on the tension function theory, the tension identification equation for a pipe subjected to a vertically downward concentrated force is derived to solve the tension identification problem when a pipe is subjected to external forces.

[0158] Secondly, the gravity constraint iterative method is used to numerically solve the key parameters of the tension identification equation when the pipeline is not subjected to a vertically downward concentrated force, in order to solve the tension identification problem when the pipeline is only subjected to its own gravity.

[0159] Next, fiber optic sensors are attached to the pipe surface to acquire strain data during the process. A tension characteristic correction coefficient matrix is ​​constructed from the original strain data matrix and its first derivative matrix. This correction coefficient matrix is ​​then substituted into the tension identification equation to correct tension identification errors caused by the lack of detailed material coefficients of the pipe.

[0160] Finally, using Ansys Workbench simulation software, a pipe model was built, and fixed constraints were applied to both ends of the pipe. Simultaneously, its own gravity load was applied to the pipe, such as... Figure 2 As shown; extract the pipeline simulation tension data and strain data, substitute the strain data into the method of this invention to obtain the tension data identified by this invention, and compare the two tension data, as shown. Figure 3 As shown. The relative error between the tension data obtained by this method and the simulated tension data is calculated, as shown. Figure 4 As shown.

[0161] This invention provides a method for identifying pipeline structure tension based on tension coefficient matrix correction. Many methods and approaches exist for implementing this technical solution; the above description is merely a preferred embodiment of the invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of this invention, and these improvements and modifications should also be considered within the scope of protection of this invention. All components not explicitly stated in this embodiment can be implemented using existing technologies.

Claims

1. A method for identifying the tension of a pipe structure based on tension coefficient matrix correction, characterized in that, Includes the following steps: Step 1: Establish a rectangular coordinate system with the lower end of the pipe as the origin to determine the pipe length. l Height difference between the two constrained ends of the pipeline h 0, Angle between the lines connecting the two ends of the pipe θ The gravitational load p on the pipeline; Step 2: Based on the tension function theory, derive the tension identification equation corresponding to the pipe being subjected to a vertically downward concentrated force; Step 3: The gravity constraint iterative method is used to numerically solve the key parameters of the tension identification equation when the pipeline is not subjected to a vertically downward concentrated force. Step 4: Attach fiber optic sensors to the pipe surface to acquire strain distribution data on the pipe surface during service. Based on the measured original strain data matrix and the first derivative matrix, construct a tension characteristic correction coefficient matrix. Substitute the tension characteristic correction coefficient matrix into the tension identification equation to correct the tension identification error caused by the lack of detailed material coefficient characteristics of the pipe. Step 4 includes: Step 4.1: Arrange three optical fibers along the circumferential direction at 0°, 120°, and 240° on the surface of the pipe. Bragg The strain sensing path of the grating sensor; composed of optical fiber Bragg Actual tensile and compressive strain of the pipeline collected by the grating sensor ; Step 4.2, the strains obtained at each discrete measurement point are as follows: Then the strain measurement matrix E is defined as: （19）， in Indicates the first n The strain obtained from the measurement of discrete measuring points; The first derivative matrix of the measured strain matrix is ​​calculated using the central difference method. : （20）， in For the first i Strain values ​​at individual fiber optic measuring points; For the first i -1 strain value at fiber optic measuring point; For the first i fiber optic measurement points X Axis coordinate values; i Values ​​range from 1 to n ; The first derivative matrix of the measured strain matrix was calculated. for: （21）， Step 4.3: Based on the measured strain matrix and its first derivative obtained from the fiber optic sensing system, the tension characteristic correction coefficient matrix is ​​obtained. K : （22）， Step 4.4: Establish the pipeline tension identification equation based on the tension correction coefficient matrix.

2. The method according to claim 1, characterized in that, In step 1, the horizontal length of pipe AB is defined as... l The gravitational load p on the pipe is calculated by the pipe density.

3. The method according to claim 2, characterized in that, Step 2 includes: Step 2.1, take any infinitesimal segment on the pipeline. ds The pipe is set to be under tension. T The horizontal component of the force is H The vertical component of the force is V Then, the static equilibrium equation of the pipeline during service can be obtained from static analysis: （1）， in d y is the differential symbol; y represents the infinitesimal segment of the pipe. Y Axis coordinates; x represents the infinitesimal segment of the pipe. X Axis coordinates; Solve equation (1) by performing a second integral, and based on the boundary conditions: the x-coordinate of the pipeline starting point. x 1=0, ordinate y 1=0; x-coordinate of the pipeline endpoint x 2=l, ordinate y 2=l, thus obtaining the linear equation of the pipe in its free state: （2）， in , Let be an arbitrary constant generated during the integration process, determined by equation (3): （3）， in arsinh Represents the inverse hyperbolic sine function; Step 2.2: Let point C be the lowest point of the pipeline. Differentiate equation (3) and take the extreme point to obtain the x-coordinate of the lowest point C. for: （4）， Substituting equation (4) into equation (2), we obtain the pipeline alignment function under the condition of the known minimum point. y ( x ): （5）， Assume any position of the pipeline Q ( x , y Point tension is T The vertical component of the tension is V The pipe tension identification equation is obtained from static analysis: （6）， （7）， According to equation (7), the pipe is at its lowest point The tension value is minimum at the lowest point. The tension gradually increases towards both ends, reaching its maximum at point B, the higher end. Therefore, the maximum tension in the pipe is... and minimum value for: （8）， in The vertical component of the tension at point B; Step 2.3: When a pipe is subjected to a concentrated force acting vertically downward, the pipe shape can be classified into two types. One type of concentrated force is applied to the center of the pipe, which is called pipe type I. The other type of concentrated force is applied near both ends of the pipe, which is called pipe type II.

4. The method according to claim 3, characterized in that, In step 2.3, when the pipe shape is linear type I, combining equations (2), (6), and (7) yields: （9）， When the pipeline configuration is linear type II, combining equations (2), (6), and (7) yields: （10）， In the formula: （11）， in Location where concentrated force is applied X Axis coordinates for Y Axis coordinates The concentrated force load on the pipeline; We obtain the results by numerically solving equations (9) and (10). The value will Substituting into equations (6) and (7), we can obtain the tension distribution characteristics of the pipeline when it is subjected to a vertically downward concentrated force.

5. The method according to claim 4, characterized in that, Step 3 includes: Step 3.1, based on the known x-coordinate of the lowest point of the pipeline y-axis To obtain information about the parameters The constraint equations are: （12）， Step 3.2: Numerical solution is performed using the gravity constraint iterative method to construct the iterative function. : （13）， Step 3.3: Solve for the parameters using the gravity constraint iterative method. When calculating the derivative of the iterative function... : （14）， in: （15）， Step 3.4, the gravity constraint iterative method format is as follows: （16）， in, The parameter obtained in the nth iteration is represented by ; The convergence factor; Indicates parameter value The function, that is, the function M and N The convergence condition is: , The set convergence precision; Step 3.5: Verify the calculation results using gravity constraint equations, assuming the pipeline... A The vertical component of the tension is V A ,pipeline B The vertical component of the tension is V B ,but: （17）， According to the gravity constraint equation ps = V A + V B The calculation results must meet the following requirements: （18）， in s Indicates the arc length of the pipe.

6. The method according to claim 5, characterized in that, In step 4.1, the actual tensile and compressive strain of the pipeline The calculation formula is: （23）， in, , , Fiber optic cables of 0°, 120°, and 240° respectively. Bragg The strain is measured by a grating sensor.

7. The method according to claim 6, characterized in that, In step 4.1, the pipe tension identification equation based on the tension correction coefficient matrix is: （24）， （25）， in Y ( x ) represents the pipe profile function based on the tension correction coefficient matrix; This represents the tension identification function for pipe structures based on the tension correction coefficient matrix.

8. An electronic device, characterized in that, It includes a processor and a memory, the memory storing program code that, when executed by the processor, causes the processor to perform the steps of the method as described in any one of claims 1 to 7.

9. A storage medium, characterized in that, It stores a computer program or instructions that, when run on a computer, perform the steps of the method as described in any one of claims 1 to 7.

Images