Pipeline stress assessment methods, devices and electronic equipment

Abstract

This invention discloses a pipeline stress assessment method, apparatus, and electronic device, relating to the field of pipeline safety assessment technology. The method includes: acquiring the target internal pressure of the target pipeline and the target ultrasonic time-series signals along each target detection path; performing stress response analysis on the target pipeline based on the target internal pressure, the target ultrasonic time-series signals, a preset internal pressure sensitivity vector, and a preset path sensitivity matrix to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response; performing a three-dimensional stress simulation on the target pipeline based on the subsets of the target external load stress response and the subsets of the target internal pressure stress response to obtain a target three-dimensional total stress tensor field; and assessing the pipeline stress state based on the target three-dimensional total stress tensor field to generate a target pipeline stress assessment result. The technical solution of this invention can improve the accuracy of pipeline stress assessment.

Description

Technical Field

[0001] This invention relates to the field of pipeline safety assessment technology, and in particular to a pipeline stress assessment method, apparatus and electronic equipment. Background Technology

[0002] With the development of pipeline safety assessment technology, ultrasonic non-destructive testing technology has emerged. This technology has the characteristics of being non-contact, having strong penetration, and being sensitive to internal defects and stress changes in materials, which in turn has given rise to the traditional method of pipeline stress detection based on unidirectional or single-point ultrasonic waves.

[0003] Traditional techniques primarily rely on deploying one or a few ultrasonic probes on the pipe surface to measure the time-of-flight variation of ultrasonic waves propagating in a single direction (usually the pipe axis). The time-of-flight difference is then converted into an average stress value along the corresponding path using the classical acoustoelastic formula. This method is essentially a discrete-point measurement; to obtain stress data over a larger area, multiple measurement points need to be pre-set on the pipe surface for repeated point-by-point testing, followed by rough estimation of the overall stress distribution trend using interpolation or other methods.

[0004] However, the current traditional method based on unidirectional single-point ultrasonic testing has obvious limitations: its single-path measurement results are essentially a mixture of stresses under the combined action of multiple loads inside and outside the pipeline, and cannot effectively separate the stress components generated by external mechanical loads and internal fluid pressures; at the same time, the discrete point measurement method is difficult to achieve continuous and high-resolution stress field reconstruction across the entire pipeline, especially lacking the ability to accurately locate and quantify stress concentration areas, resulting in a lack of accuracy in stress state assessment. Summary of the Invention

[0005] This invention provides a pipeline stress assessment method, apparatus, and electronic device, which can improve the accuracy of pipeline stress assessment.

[0006] According to one aspect of the present invention, a method for assessing pipeline stress is provided, the method comprising: Acquire the internal pressure of the target pipeline and the ultrasonic timing signals of each target detection path; Based on the target internal pressure, the target ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, stress response analysis is performed on the target pipeline to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. Based on the target external load stress response subset and the target internal pressure stress response subset, a three-dimensional stress simulation is performed on the target pipeline to obtain the target three-dimensional total stress tensor field. Based on the target three-dimensional total stress tensor field, the target pipeline stress state is evaluated, and the target pipeline stress evaluation result is generated.

[0007] According to another aspect of the present invention, a pipe stress assessment device is provided, the device comprising: The data acquisition module is used to acquire the internal pressure of the target pipeline and the ultrasonic time-series signals of each target detection path; The data analysis module is used to perform stress response analysis on the target pipeline based on the target internal pressure, the target ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. The pipeline stress synthesis module is used to perform three-dimensional stress simulation on the target pipeline based on the target external load stress response subset and the target internal pressure stress response subset, and obtain the target three-dimensional total stress tensor field. The pipeline stress assessment module is used to assess the pipeline stress state of the target pipeline based on the target three-dimensional total stress tensor field and generate the target pipeline stress assessment result.

[0008] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor, which enables the at least one processor to perform the pipeline stress assessment method according to any embodiment of the present invention.

[0009] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions for causing a processor to execute and implement the pipe stress assessment method according to any embodiment of the present invention.

[0010] According to another aspect of the present invention, a computer program product is provided, the computer program product comprising a computer program that, when executed by a processor, implements the pipeline stress assessment method according to any embodiment of the present invention.

[0011] The technical solution of this invention acquires the internal pressure and multi-path ultrasonic time-series signals of the target pipeline simultaneously. By using a preset path sensitivity matrix and a preset internal pressure sensitivity vector, the physical source of the pipeline stress response under complex conditions is first decoupled to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. Then, through field reconstruction technology based on the acoustoelastic principle, a high-resolution and quantifiable target three-dimensional total stress tensor field is generated. The pipeline stress state is evaluated by the target pipeline through the target three-dimensional total stress tensor field. This achieves a deep integration of non-destructive testing, solid mechanics, data analysis, and engineering decision-making, changing the limitations of traditional single-point pipeline inspection and improving the accuracy, comprehensiveness, and engineering practicality of pipeline stress state evaluation.

[0012] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0013] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0014] Figure 1 This is a flowchart of a pipeline stress assessment method provided in Embodiment 1 of the present invention; Figure 2 This is a flowchart of a pipeline stress assessment method provided in Embodiment 2 of the present invention; Figure 3 This is a schematic diagram of the structure of a pipeline stress assessment device provided in Embodiment 3 of the present invention; Figure 4 This is a schematic diagram of the structure of an electronic device that implements the pipeline stress assessment method of this invention. Detailed Implementation

[0015] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0016] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0017] Example 1 Figure 1 This is a flowchart of a pipeline stress assessment method provided in Embodiment 1 of the present invention. This embodiment of the invention is applicable to situations where pipeline stress is assessed based on ultrasonic multipath methods. The method can be executed by a pipeline stress assessment device, which can be implemented in hardware and / or software. This device can be configured in an electronic device that carries pipeline stress assessment functions, such as a detection terminal and / or a server. It is understood that this method can be applied to a detection terminal, a server, or a system including both a detection terminal and a server, and is implemented through the interaction between the detection terminal and the server.

[0018] See Figure 1 The pipe stress assessment method shown includes: S101. Obtain the target internal pressure in the target pipeline and the target ultrasonic timing signals of each target detection path.

[0019] The target pipeline is the pipeline to be tested for pipeline stress assessment. The target internal pressure is the internal pressure collected in real time by a fluid pressure sensor inside the target pipeline. The target internal pressure is used to characterize the static pressure of the fluid inside the target pipeline. Optionally, the target internal pressure can be a scalar.

[0020] The target ultrasonic timing signal is the ultrasonic timing signal collected by an array of ultrasonic sensors deployed on the surface of the target pipe. The ultrasonic sensor array consists of multiple ultrasonic probes. The target detection path is a detection path consisting of two ultrasonic probes. Each pair of ultrasonic probes includes a transmitting probe and a receiving probe. Optionally, for any detection path, the original analog signals collected by the ultrasonic probes within a single target detection path are converted from analog to digital signals to form discrete target ultrasonic timing signals. For example, the mathematical form of the target ultrasonic timing signal can be expressed as: Where i is the index of the target detection path; t is the sampling time sequence number.

[0021] For example, the target ultrasonic time-series signal for each target detection path includes multi-dimensional physical information such as the time of flight, waveform phase, and amplitude of the ultrasonic wave propagating along the target detection path. Among these, the time of flight is most sensitive to stress changes. The time of flight is the direct basis for calculating changes in sound velocity and inverting the magnitude of stress. The waveform phase reflects the minute waveform distortions and frequency component changes caused by stress during ultrasonic wave propagation. Amplitude attenuation is related to the microstructure, defect distribution, and energy dissipation within the target pipe material. Amplitude can serve as a reference to assist in judging the material homogeneity or damage of the target pipe.

[0022] Specifically, the internal pressure of the target pipe is collected in real time by a fluid pressure sensor inside the target pipe. Ultrasonic time-series signals of the target along each detection path are collected by an array of ultrasonic sensors deployed on the surface of the target pipe.

[0023] S102. Based on the target's internal pressure, the target's ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, stress response analysis is performed on the target pipeline to obtain the target's external load stress response subset and the target's internal pressure stress response subset.

[0024] The preset path sensitivity matrix is ​​used to characterize the mapping relationship between the standard external load on the target pipeline and the ultrasonic timing signal. Specifically, the preset path sensitivity matrix characterizes the theoretical change in the ultrasonic propagation time of a certain target detection path when the target pipeline is subjected to a certain standard external load. For example, the preset path sensitivity matrix can be an M×N matrix, which can be denoted as... Where M is the total number of target detection paths; N is the total number of types of standard external payloads. Correspondingly, Let be the theoretical change in ultrasonic propagation time of the i-th target detection path when the target pipeline is subjected to the j-th standard external load.

[0025] The standard external load is a standard form of a certain type of external load. It is used to characterize the standard pressure value of a particular type of external load. The total number of standard external load types can be set and adjusted according to the environment of the target pipeline and the standard external load conditions experienced during actual operation. For example, the types of standard external loads may include, but are not limited to, unit axial tensile force, unit axial compressive force, unit circumferential compressive force, unit positive bending moment (causing tension on one side of the pipeline and compression on the other), unit torque, and unit uniform external pressure simulating lateral soil pressure.

[0026] In an optional embodiment of the present invention, the process of generating the preset path sensitivity matrix includes: acquiring the target pipeline structural parameters, a set of standard external loads, and the target detection path; for a single standard external load, performing stress field simulation on the target pipeline based on the target pipeline structural parameters to obtain the first theoretical stress distribution field of the target pipeline under the standard external load; for a single standard external load, using a preset set of acoustoelastic physics formulas, integrating the first theoretical stress distribution field along the target detection path to obtain the first theoretical time change of the target ultrasonic time-series signal under the standard external load; and generating the preset path sensitivity matrix based on the first theoretical time change of each target ultrasonic time-series signal under each standard external load.

[0027] Target pipe structural parameters are used to characterize the geometric and structural features of the target pipe. Optionally, the target pipe structural parameters can be preset and adjusted based on pipe characteristics detected in actual operation. For example, target pipe structural parameters include, but are not limited to, the outer diameter, wall thickness, material density, Young's modulus of elasticity under reference conditions, and Poisson's ratio of the target pipe.

[0028] The standard external load set is a predefined set of standardized external mechanical conditions used to quantify the sensitivity of a target pipeline to external loads. The standard external load set can be set according to the types of external loads present in the pipeline during actual operation. The standard external load set contains multiple standard external loads. For example, any standard external load in the standard external load set can be denoted as j, representing the pressure of the j-th type of standard external load.

[0029] The first theoretical stress distribution field is used to characterize the stress tensor of a certain standard external load on the surface of a target pipe. The first theoretical stress distribution field contains the stress tensor at each coordinate point on the surface of the target pipe. For example, the first theoretical stress distribution field can be denoted as... Correspondingly, the stress tensor of the j-th standard external load at any coordinate point k on the surface of the target pipe can be denoted as: .

[0030] The preset acoustoelasticity formula set is a collection of physical laws describing the quantitative relationship between stress and changes in ultrasonic velocity. Optionally, the preset acoustoelasticity formula set can be set and adjusted based on experience and knowledge of acoustoelasticity in practical work.

[0031] The first theoretical time variation is the theoretical time variation of the ultrasonic wave's flight time along the target detection path when the target pipe is subjected only to a standard external load. The first theoretical time variation quantifies the purely theoretical influence of the standard external load on the ultrasonic timing signal. For example, the first theoretical time variation can be denoted as... Where i is the i-th target detection path; j is the j-th standard external payload.

[0032] Specifically, it involves obtaining predetermined target pipeline structural parameters, standard external load sets, and target detection paths.

[0033] Specifically, finite element analysis software is used to establish a high-fidelity three-dimensional digital model of the pipe segment using the structural parameters of the target pipe. For each standard external load in the set of standard external loads, the corresponding standard external load is applied to the three-dimensional digital model of the pipe segment. Using finite element analysis software and based on the theory of elasticity, static finite element analysis is performed on the target pipe to obtain the first theoretical stress distribution field on the outer surface of the entire target pipe structure caused by the standard external load. The three-dimensional digital model of the pipe segment can accurately reflect the relationship between the geometry and material constitutive model of the target pipe.

[0034] Specifically, for a single standard external load, using a preset set of acoustoelastic physical formulas, the first theoretical stress distribution field is numerically integrated along the target detection path to obtain the first theoretical time change of the target ultrasonic time-series signal under the standard external load.

[0035] For example, the time-of-flight variation of an ultrasonic wave propagating in a specific direction is related to the strain tensor integral of the region traversed by the target detection path. The target detection path can restrict the propagation of the ultrasonic wave in a specific direction.

[0036] For example, the following formula can be used to represent the pre-defined physical formula for acoustic elasticity: ; In the formula, This is the change in flight time of an ultrasonic wave propagating in a specific direction, i.e., the first theoretical time change. The reference value for the flight time of an ultrasonic wave propagating in a specific direction is used to normalize the change in the flight time of the ultrasonic wave propagating in a specific direction, that is, to convert the measured absolute time offset into the theoretically required relative rate of change; where, the absolute time offset is the change in the flight time of the ultrasonic wave propagating in a specific direction; K is the acoustoelastic coefficient tensor; and L is the target detection path. Let be the strain tensor; where, Hooke's law can be used to determine the relationship between stress tensor and stress tensor. (That is, the first theoretical stress distribution field of the target pipeline under standard external load includes the stress tensor of each coordinate point on the surface of the target pipeline) and elastic constants are obtained.

[0037] Specifically, the first theoretical time variation of the ultrasonic time-series signals of each target under each standard external load is arranged to generate a preset path sensitivity matrix.

[0038] For example, for M target detection paths and N standard external payloads, a total of The first theoretical time change The rows of the matrix can be mapped to the target detection path, and the columns to the standard external load, for each first theoretical time variation. The paths are arranged to obtain a preset path sensitivity matrix. The path sensitivity matrix is ​​then filled with the [missing information - likely a specific element or component]. Line number The elements of the column are Each element of the preset path sensitivity matrix It can be viewed as a "feature response vector," which fully describes the response of the target pipeline when subjected to the first... Under standard external loads, the first theoretical time variation of all M target detection paths is calculated. Here, the i-th row of the preset path sensitivity matrix corresponds to the i-th target detection path; the j-th column corresponds to the j-th standard external load. Therefore, a... The preset path sensitivity matrix S.

[0039] This scheme establishes an accurate parameterized finite element model of the pipeline (i.e., a three-dimensional digital model of the pipeline segment) based on the target pipeline parameters. Combined with the laws of acoustoelastic physics, it systematically simulates and calculates a set of theoretical responses of ultrasonic multipath under standard mechanical conditions, thereby constructing a preset path sensitivity matrix. This enables the pre-calibration of the complex pipeline mechanical behavior and ultrasonic propagation effect in digital space, providing a reliable quantitative reference for the subsequent accurate separation and identification of stress responses from different sources from actual mixed signals.

[0040] The preset internal pressure sensitivity vector is used to characterize the theoretical change in the propagation time of the target detection path when the standard internal pressure of the target pipeline increases by 1 unit (e.g., 1 MPa). For example, the preset internal pressure sensitivity vector can be a... The column vector. The preset internal pressure sensitivity vector can be represented as... Among them, the elements in the preset internal pressure sensitivity vector This represents the theoretical impact of standard internal pressure on the i-th target detection path.

[0041] Optionally, the process of generating the preset internal pressure sensitivity vector includes: acquiring the target pipeline structural parameters, standard internal pressure, and target detection path; performing stress field simulation on the target pipeline based on the target pipeline structural parameters and standard internal pressure to obtain the second theoretical stress distribution field of the target pipeline under standard internal pressure; integrating the second theoretical stress distribution field along the target detection path using a preset set of acoustoelastic physics formulas to obtain the second theoretical time change of each target ultrasonic time-series signal under standard internal pressure; and generating the preset internal pressure sensitivity vector based on the second theoretical time change of each target ultrasonic time-series signal under each standard internal pressure.

[0042] Standard internal pressure is a predefined unit of internal pressure used to quantify the sensitivity of a target pipeline to internal pressure.

[0043] The second theoretical stress distribution field is used to characterize the stress tensor of the standard internal pressure at any coordinate point on the surface of the target pipe.

[0044] The second theoretical time variation is the theoretical time variation of the ultrasonic wave's flight time along the target detection path when the target pipe is only subjected to standard internal pressure. This second theoretical time variation quantifies the purely theoretical influence of standard internal pressure on the ultrasonic timing signal.

[0045] Specifically, the target pipeline structure parameters, standard internal pressure, and target detection path are obtained in advance.

[0046] Specifically, finite element analysis software is used to establish a high-fidelity three-dimensional digital model of the pipeline segment using the structural parameters of the target pipeline. A corresponding standard internal pressure is applied to the three-dimensional digital model of the pipeline segment. Based on the theory of elasticity, static finite element analysis is performed on the target pipeline using finite element analysis software to obtain the second theoretical stress distribution field on the outer surface of the entire target pipeline structure caused by the standard internal pressure.

[0047] Specifically, by using a set of preset acoustic elastic physical formulas, the second theoretical stress distribution field is numerically integrated along the target detection path to obtain the second theoretical time change of the ultrasonic time sequence signal of each target under standard internal pressure.

[0048] Specifically, the first theoretical time change of each target ultrasonic time signal under each standard external load is arranged to generate a preset internal pressure sensitivity vector.

[0049] By establishing an accurate parameterized finite element model of the pipeline using the target pipeline parameters, and combining it with the laws of acoustoelastic physics, a set of theoretical responses of ultrasonic multipath under standard mechanical conditions are systematically simulated and calculated. This allows for the construction of a preset internal pressure sensitivity vector, enabling the pre-calibration of the complex pipeline mechanical behavior and ultrasonic propagation effects in digital space. This provides a reliable quantitative reference for the subsequent accurate separation and identification of stress responses from different sources from actual mixed signals.

[0050] A subset of the target external load stress response is used to characterize the response of the target ultrasonic time-series signal to the external load on the target pipeline. For example, the target external load stress response subset can be denoted as... The target external load stress response subset can be a... A real-valued column vector; M is the total number of target detection paths; Used to characterize the The time-of-flight variation of ultrasonic waves along the target detection path, solely contributed by the combined effect of all external payloads; among which, It is a scalar.

[0051] A subset of the target's internal pressure stress response is used to characterize the response of the target's ultrasonic time-series signal to the internal pressure of the target pipeline. For example, the subset of the target's internal pressure stress response can be denoted as... Among them, the subset of internal pressure stress response of the target can be a A real-valued column vector; M is the total number of target detection paths; Used to characterize the The time-of-flight variation of ultrasonic waves along the target detection path, solely contributed by the internal pressure of the target; among which, It is a scalar.

[0052] Specifically, by using the target ultrasonic time-series signal and the target internal pressure, combined with the pre-determined preset path sensitivity matrix and preset internal pressure sensitivity vector, the overall stress response state of the target pipeline is analyzed to obtain the target external load stress response subset and the target internal pressure stress response subset.

[0053] In an optional embodiment of the present invention, stress response analysis is performed on the target pipeline based on the target internal pressure, the target ultrasonic time-series signal, a preset internal pressure sensitivity vector, and a preset path sensitivity matrix to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. This includes: calculating the target time change vector based on each target ultrasonic time-series signal and a reference signal; constructing a target stress response separation equation based on the target time change vector, the target internal pressure, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix; estimating the target external load coefficient vector based on the target stress response separation equation; calculating the target external load stress response subset based on the target external load coefficient vector and the preset path sensitivity matrix; and calculating the target internal pressure stress response subset based on the target internal pressure and the preset internal pressure sensitivity vector.

[0054] The reference reference signal is an ultrasonic timing signal acquired through various target detection paths under a known, low-stress reference state of the target pipeline. For example, the reference state includes a state of shutdown and depressurization or no external interference. Optionally, the reference reference signal can be obtained by testing the target pipeline during actual operation. The reference reference signal can be pre-stored in this device. For example, the reference reference signal can be denoted as... Where i is the index of any target detection path.

[0055] The target time variation vector is used to comprehensively characterize the observable total effect of the current complex stress state of the target pipeline on the ultrasonic waves emitted by all target detection paths. The number of elements in the target time variation vector is the same as the number of target ultrasonic time-series signals. That is, the number of elements in the target time variation vector is the same as the number of target detection paths. The target time variation vector contains the target time variation of each target detection path. Here, the target time variation is the target time variation of the flight time of a single target ultrasonic time-series signal relative to the flight time of a reference signal.

[0056] The target stress response separation equation is a mathematical model used to describe the source of the observed vector, the target time change vector. The target stress response separation equation can be constructed based on the principle of linear superposition of acoustoelastic effects. This principle holds under the assumption of small deformation.

[0057] The target time change vector can be considered as the superposition of two linear contributions: one part is the sum of contributions from various external loads, and the other part is the contribution from internal pressure.

[0058] The contribution of the external load can be expressed as the product of a preset path sensitivity matrix and an unknown target external load coefficient vector. The target external load coefficient vector represents the amplification factor of various types of standard external loads in practice. This can be understood as the actual external load magnitude being a multiple of the target external load coefficient of the standard external load. For example, the contribution of the external load can be denoted as... For example, the target external load coefficient vector It is The unknown column vector. The target external load coefficient vector can be denoted as... Where N is the total number of types of external loads. This represents the magnification factor of the j-th standard external load in reality, that is, the actual load size is a fraction of the standard external load. This can be understood as the magnitude of the j-th type of external load actually borne by the target pipeline being 1 / 2 times the standard external load. times. It is a dimensionless coefficient.

[0059] The contribution of internal pressure can be expressed as the product of a preset internal pressure sensitivity vector and the real-time measured target internal pressure. For example, the contribution of internal pressure can be denoted as... The preset internal pressure sensitivity vector can be denoted as: The internal pressure of the target can be denoted as p.

[0060] Therefore, the target stress response separation equation can be constructed as follows: ; In the formula, S is the target time change vector; S is the preset path sensitivity matrix; It is the target external load coefficient vector to be solved; It is the preset internal pressure sensitivity vector; p is the target internal pressure; It is a subset of the target's external load stress response; It is a subset of the internal pressure stress response of the target.

[0061] In the target stress response separation equation, , , and It is a known quantity, and These are the unknowns to be solved. The core objective of solving the target stress-response separation equation is to observe the target's time-varying vector. Estimate the target external load coefficient vector This allows us to separate the response to the external load.

[0062] To account for measurement noise and model error, a noise residual vector is introduced. This noise residual vector describes the observational portion of the target time-varying quantity vector—the source of the observation vector—that is not explained by the mathematical model. For example, the noise residual vector includes measurement noise, model error, and unmodeled physical effects.

[0063] Therefore, the constructed target stress response separation equation takes the form of: ; In the formula, S is the target time change vector; S is the preset path sensitivity matrix; It is the target external load coefficient vector to be solved; It is the preset internal pressure sensitivity vector; p is the target internal pressure; It is a subset of the target's external load stress response; among which, Depend on The calculation shows that this represents the portion of the target time change vector purely contributed by external loads. It is a subset of the internal pressure stress response of the target; among which, Depend on The calculation shows that represents the portion of the target time change vector purely contributed by internal pressure; e is the noise residual vector.

[0064] In the target stress response separation equation, , , and It is a known quantity, and Let and be the unknowns to be solved. The core objective of solving the target stress response separation equation is to observe the target time change vector. Estimate the target external load coefficient vector The response to the external load is separated by the noise residual vector e.

[0065] By clearly defining the mathematical form of the target stress response separation equation and the physical and mathematical meaning of each parameter, a rigorous and transparent theoretical framework is provided for the entire stress decoupling process, transforming complex engineering problems into computable optimization problems.

[0066] Specifically, for a single target ultrasonic time-series signal, a cross-correlation algorithm is used to compare the target ultrasonic time-series signal corresponding to a single target detection path with a reference signal to calculate the target time change. The time changes of each target are then sorted according to the order of the target detection paths to obtain a target time change vector.

[0067] For example, for a single target detection path, the target ultrasonic time-series signal corresponding to the i-th target detection path can be calculated. With reference signal Cross-correlation function between We can find the time shift corresponding to the peak position of this cross-correlation function. The time shift corresponding to this peak position. It is the change in the target time of flight of the target ultrasonic timing signal relative to the flight time of the reference signal, i.e. Perform the above operation on all M target detection paths to obtain M target time changes. According to the order of the target detection path, the time changes of M targets are... Arrange them into an M-dimensional column vector to obtain the target time change observation vector. .

[0068] Specifically, based on the target time change observation vector, target internal pressure, preset internal pressure sensitivity vector, and preset path sensitivity matrix, a target stress response separation equation is constructed based on the assumption of linear superposition principle.

[0069] Specifically, L1-norm regularized regression or numerical optimization algorithms are used to solve the target stress-response separation equation and estimate the target external load coefficient vector. For example, the numerical optimization algorithm could be the coordinate descent method.

[0070] For example, when using L1 norm regularized regression to solve for the target external load coefficient vector of the target stress response separation equation, it is not only required that the fitting error... Minimize, also for the target external load coefficient vector L1 norm Imposing penalties. This is based on a reasonable physical prior: in actual target pipelines, only a few types of target external loads are often active (i.e., (It is a sparse vector).

[0071] Therefore, the optimization problem of the target external load coefficient vector in the target stress response separation equation can be expressed as: ; In the formula, It is the solution to the optimization problem, that is, obtaining the objective external load coefficient vector that minimizes the objective function. That is, the optimal target external load coefficient vector ; S is the target time change observation vector; S is the preset path sensitivity matrix; It is the target external load coefficient vector to be solved; It is the preset internal pressure sensitivity vector; p is the target internal pressure; It is a preset regularization parameter used to control the strength of sparsity.

[0072] Specifically, the product between the target's external load coefficient vector and the preset path sensitivity matrix is ​​calculated to obtain a subset of the target's external load stress response. The product between the target's internal pressure and the preset internal pressure sensitivity vector is calculated to obtain a subset of the target's internal pressure stress response.

[0073] For example, the following formula can be used to calculate the subset of the target external load stress response: ; In the formula, It is a subset of the target's external load stress response, representing the portion of the target's time variation vector attributable to the external load; S is the preset path sensitivity matrix. It is the optimal target external load coefficient vector .

[0074] For example, the following formula can be used to calculate a subset of the target's internal pressure stress response: ; In the formula, It is a subset of the target's internal pressure stress response, representing the portion of the target's time-varying vector that can be attributed to internal pressure; is the preset internal pressure sensitivity vector; p is the target internal pressure.

[0075] Correspondingly, the optimal noise residual vector This can be considered as noise that is not explained by the mathematical model.

[0076] This scheme obtains a precise target time change vector by comparing the real-time detected target ultrasonic time-series signal with a reference signal. Using a pre-calibrated preset path sensitivity matrix and preset internal pressure sensitivity vector, a target stress response separation equation based on linear superposition and sparse prior is constructed. Finally, through regression solving, the scheme achieves high-precision separation of the target external load stress response subset dominated by external mechanical load and the target internal pressure stress response subset dominated by internal fluid pressure from the mixed ultrasonic total response. This successfully transforms the complex physical source problem of pipeline stress into a computable sparse signal decomposition problem, providing reliable input data for the subsequent independent reconstruction of the external stress field (i.e., the target external load stress distribution vector) and the internal stress field (i.e., the target internal pressure stress distribution vector).

[0077] S103. Based on the target external load stress response subset and the target internal pressure stress response subset, perform three-dimensional stress simulation on the target pipeline to obtain the target three-dimensional total stress tensor field.

[0078] The target three-dimensional total stress tensor field is used to comprehensively and quantitatively characterize the continuous stress distribution across the entire target pipe segment. The target three-dimensional total stress tensor field stores the complete planar stress state of each target acquisition location in the global coordinate system. For example, the target three-dimensional total stress tensor field is a spatial function. For instance, the target three-dimensional total stress tensor field can be denoted as: For any point on the surface of the target pipe The total three-dimensional stress tensor field of the target provides the complete plane stress state at each coordinate point. The stress tensor at a single coordinate point (i.e., the target acquisition location) can be denoted as... In the formula, Characterizes the normal stress components along the axial direction of the target pipe; where tensile stress is positive and compressive stress is negative; Characterizes the normal stress components along the circumferential direction of the target pipeline; Characterizes the shear stress components within the plane of the target pipe surface.

[0079] Specifically, a field reconstruction algorithm based on physical principles is used to perform three-dimensional stress simulation on the target's external load stress response subset and internal pressure stress response subset, generating a continuous stress distribution on the target pipe surface and obtaining the target's three-dimensional total stress tensor field.

[0080] For example, for each discretized coordinate point on the surface of the target pipeline, a large set of linear equations based on acoustoelastic physics is solved using subsets of the target's external load stress response and internal pressure stress response. This allows the calculation of the stress tensors at that coordinate point caused by the external load and internal pressure, respectively. By vector superimposing these two stress tensors at each coordinate point, the final three-dimensional total stress tensor field covering the entire target pipeline is obtained.

[0081] In an optional embodiment of the present invention, a three-dimensional stress simulation is performed on the target pipeline based on a subset of the target's external load stress response and a subset of the target's internal pressure stress response to obtain the target's three-dimensional total stress tensor field. This includes: for a single target detection path, constructing a target stress tensor equation for the target detection path based on the subset of the target's external load stress response, the subset of the target's internal pressure stress response, and preset acoustoelastic physical coefficients corresponding to each target acquisition location on the target detection path; solving the target stress tensor equation for each target detection path to obtain the target's external load stress distribution vector and the target's internal pressure stress distribution vector at each target acquisition location; and superimposing the target's external load stress distribution vector and the target's internal pressure stress distribution vector at each target acquisition location to obtain the target's three-dimensional total stress tensor field.

[0082] Each target ultrasonic timing signal along a target detection path is associated with a target acquisition location. The target acquisition location characterizes the spatial coordinates of each point along the target pipeline traversed by the target detection path. The target acquisition location uniquely identifies the spatial location information for sampling within the target acquisition path. For any target ultrasonic timing signal along a target detection path, the target acquisition location must at least include the position coordinates of the transmitting probe along the target detection path. and the position coordinates of the receiving probe Optionally, the target acquisition location can be its position coordinates in a three-dimensional coordinate system referenced to the pipe axis. Optionally, the target acquisition location can be pre-determined and stored. For example, the target acquisition location can be denoted as... Correspondingly, the target ultrasonic time-series signal can include the target acquisition location. The complete dataset of the target ultrasonic time-series signal can be represented as follows: Where M is the total number of target detection paths; i is the index of any target detection path; It is the first Ultrasonic timing signal of the target along the target detection path; Let be the target acquisition location of the ultrasonic timing signal of the target along the i-th target detection path.

[0083] The preset acoustoelastic physical coefficients are coefficient tensors related to the waveform and propagation direction of the target ultrasonic time-series signal. Optionally, the preset acoustoelastic physical coefficients can be calculated based on the acoustoelastic constants and elastic constants of the target pipeline material, combined with the specific spatial geometric relationship of each "path-acquisition location" pair, according to the acoustoelastic theory formula. The acoustoelastic constants can be third-order elastic constants. The specific spatial geometric relationship of each "path-acquisition location" pair can be determined by the target acquisition location and the three-dimensional pipeline segment digital model. Optionally, the preset acoustoelastic physical coefficients can be directly calibrated by applying a unit test load to the three-dimensional pipeline segment digital model and performing acoustic-mechanical coupling simulation.

[0084] The preset acoustic elastic physical coefficients include external acoustic elastic inversion coefficients and internal acoustic elastic inversion coefficients. For example, the mathematical form of the preset acoustic elastic physical coefficients can be expressed as follows: In the formula, K is the acoustic elastic physical coefficient; These are the external acoustic elasticity inversion coefficients; These are the internal acoustoelastic inversion coefficients; i represents the i-th target detection path; k represents the k-th target acquisition location, i.e., the k-th coordinate point; m represents the component index of the stress tensor ( They can correspond to each ).

[0085] The physical meaning of the external acoustoelastic inversion coefficient is the theoretical change in the flight time of a target detection path caused by a unit value (e.g., 1 Pa) of a certain component of the stress tensor at a target acquisition location, considering only external loads. The external acoustoelastic inversion coefficient is a scalar. For example, This represents the theoretical change in the flight time of the i-th target detection path caused by the m-th component of the stress tensor at the k-th coordinate point being a unit value (1 Pa).

[0086] The internal acoustoelastic inversion coefficients, when considering only internal pressure, represent the theoretical change in the flight time of a target detection path caused by a unit value (e.g., 1 Pa) of a component of the stress tensor at a target acquisition location. The internal acoustoelastic inversion coefficients are scalars. For example, This represents the theoretical change in the flight time of the i-th target detection path caused by considering only external loads and only internal pressures, when the m-th component of the stress tensor at the k-th coordinate point is a unit value (1 Pa).

[0087] For a single target detection path, based on a three-dimensional coordinate system with the target pipe axis as the reference axis, the stress tensor of each discretized target acquisition location on the outer surface of the target pipe is calculated. For example, the stress tensor of coordinate point k at each target acquisition location can be denoted as... According to the theory of acoustoelasticity, the change in flight time of a target detection path is proportional to the line integral of the strain tensor (obtained from the stress tensor using Hooke's law) of all nodes traversed by the target detection path.

[0088] The target stress tensor equation characterizes the relationship between the stress response and the stress tensor at the target acquisition location. For example, the stress response includes the target's external load stress response and the target's internal pressure stress response. Accordingly, the target stress tensor equation can include the external load target stress tensor quantum equation and the internal pressure target stress tensor quantum equation. The external load target stress tensor quantum equation characterizes the relationship between the target's external load stress response and the stress tensor at the target acquisition location caused by the external load. The internal pressure target stress tensor quantum equation characterizes the relationship between the target's internal pressure stress response and the stress tensor at the target acquisition location caused by the internal pressure.

[0089] The target external load stress distribution vector is the stress tensor caused by the external load at a certain target acquisition location. The target internal pressure stress distribution vector is the stress tensor caused by the internal pressure at a certain target acquisition location. The target external load stress distribution vector and the target internal pressure stress distribution vector together constitute the discretized stress field on the pipe surface caused by different physical factors.

[0090] Specifically, for a single target detection path, based on the target external load stress response subset and the external acoustoelastic inversion coefficients in the preset acoustoelastic physical coefficients corresponding to each target acquisition position on the target detection path, the target stress tension quantum equation of the target detection path is constructed.

[0091] For example, the quantum equation for the external load target stress tension along the target detection path can be expressed by the following formula: ; In the formula, Let be the stress response of the i-th target external load in the subset of target external load stress responses, representing the time-of-flight variation of the ultrasonic wave along the i-th target detection path, which is purely contributed by the combined effect of all external loads; k is the acquisition position of the k-th target; K is the total number of target acquisition positions; This represents the comprehensive result of the external acoustoelastic inversion coefficients corresponding to the stress tensor at the k-th target acquisition location; it can be understood as... This is the theoretical change in the flight time of the i-th target detection path caused by each component of the stress tensor induced by the external load at the k-th coordinate point being a unit value; It is obtained by synthesizing the external acoustoelastic inversion coefficients. ; This is the target external load stress distribution vector at the k-th target acquisition location.

[0092] Accordingly, for a single target detection path, based on the subset of internal pressure stress response of the target and the internal acoustoelastic inversion coefficients in the preset acoustoelastic physical coefficients corresponding to each target acquisition position on the target detection path, the internal pressure target stress tension quantum equation of the target detection path is constructed.

[0093] For example, the quantum equation for the internal pressure target stress tension along the target detection path can be expressed by the following formula: ; In the formula, Let i be the target internal pressure stress response in the subset of target internal pressure stress responses, representing the i-th target internal pressure stress response. The time-of-flight variation of the ultrasonic waves along the target detection path, solely contributed by the combined effect of all internal pressures; k is the acquisition location of the kth target; K is the total number of target acquisition locations; This represents the comprehensive result of the external acoustoelastic inversion coefficients corresponding to the stress tensor at the k-th target acquisition location; it can be understood as... This is the theoretical change in the flight time of the i-th target detection path caused by each component of the stress tensor caused by internal pressure at the k-th target acquisition location being a unit value; Obtained through internal acoustoelastic inversion coefficients. ; This is the target internal pressure stress distribution vector at the k-th target acquisition location.

[0094] Specifically, for a single target detection path, the target stress tensor equation of the target detection path is constructed based on the quantum equations of the external load target stress tension and the internal pressure target stress tension.

[0095] Specifically, for the external load target stress tensor quantum equation in the target stress tensor equation of all target detection paths, all unknown target external load stress distribution vectors are arranged into an external load long vector according to the order of target acquisition positions, forming the external load equation set.

[0096] For example, the external load equations can be represented by the following formulas: ; In the formula, It is the external load coefficient matrix; This is the long vector of the external load, which comprehensively represents the stress distribution vector of the external load on each target; It is a subset of the target's external load stress response.

[0097] Correspondingly, for the internal pressure target stress tensor quantum equation in the target stress tensor equation of all target detection paths, all unknown target internal pressure stress distribution vectors are arranged into an internal pressure long vector according to the order of target acquisition positions, forming the internal pressure equation set.

[0098] For example, the internal pressure equations can be represented by the following formulas: ; In the formula, It is the internal pressure coefficient matrix; This is the long vector of internal pressure, which comprehensively represents the internal pressure stress distribution vector of each target. It is a subset of the internal pressure stress response of the target.

[0099] Specifically, the ridge regression regularization method is used to solve for the long vector of the external load in the external load equation system, thus obtaining the long vector of the external load. The long vector of the external load is then restored according to the order of the target acquisition positions to obtain the target external load stress distribution vector at each target acquisition position.

[0100] For example, taking the external load equations as an example, the ridge regression regularization method (e.g., the conjugate gradient method) is used to transform the problem into minimizing the objective function. For example, the objective function corresponding to the external load equations can be expressed by the following formula: ; In the formula, the first term This is the data fitting term, which requires the solution to conform as closely as possible to the observed data; among them, It is the external load coefficient matrix; This is the long vector of the external load, which comprehensively represents the stress distribution vector of the external load on each target; It is a subset of the target's external load stress response; the second term It is a regularization term; among which, is the preset regularization parameter; L is a discrete Laplacian operator matrix. The regularization term penalizes drastic spatial fluctuations in the solution, thus favoring a smooth, continuous stress distribution to ensure the solution conforms to most engineering realities.

[0101] Therefore, an estimated solution for the long vector stability of the external load can be obtained. Following the order of target acquisition locations, the solved external load long vectors are restored to obtain the target external load stress distribution vector at each target acquisition location k. .

[0102] Accordingly, the ridge regression regularization method is used to solve for the long internal pressure vector in the internal pressure equation system, thus obtaining the long internal pressure vector. The long internal pressure vector is then restored according to the order of the target acquisition locations to obtain the target internal pressure stress distribution vector for each target acquisition location.

[0103] For example, the exact same procedure as the external load equations can be applied to the internal pressure equations to obtain the target internal pressure stress distribution vector at each target acquisition location k. .

[0104] Specifically, for a single target acquisition location, the external load stress distribution vector and the internal pressure stress distribution vector are superimposed to obtain the target's three-dimensional total stress tensor at that location. Based on the target's three-dimensional total stress tensor at each acquisition location, a target three-dimensional total stress tensor field is constructed.

[0105] For example, the following formula can be used to superimpose the target's external load stress distribution vector and the target's internal pressure stress distribution vector to obtain the target's three-dimensional total stress tensor at a single target acquisition location: ; In the formula, Let the total three-dimensional stress tensor of the target at the k-th target acquisition location be denoted as . This represents the target external load stress distribution vector at the k-th target acquisition location; The vector representing the internal pressure stress distribution of the target at the k-th target acquisition location.

[0106] This scheme maps the separated path response data—namely, the target external load stress response subset and the target internal pressure stress response subset—based on acoustoelastic physical relationships and establishes a target stress tensor equation for the target detection path with the target external load stress distribution vector and the target internal pressure stress distribution vector at the target acquisition location as unknowns. By solving for the target external load stress distribution vector and the target internal pressure stress distribution vector at each target acquisition location on the target pipeline surface, and then superimposing the stresses, the overall three-dimensional total stress tensor field of the target is obtained. This achieves the dimensionality upgrade and reconstruction from limited, one-dimensional path response data to comprehensive, two-dimensional surface stress field data (i.e., the target three-dimensional total stress tensor field). This ensures that the reconstructed target three-dimensional total stress tensor field is not only complete but also has clear physical authenticity and engineering interpretability.

[0107] S104. Based on the target three-dimensional total stress tensor field, evaluate the pipeline stress state of the target pipeline and generate the target pipeline stress evaluation results.

[0108] The target pipeline stress assessment result is the target pipeline risk assessment result. For example, the target pipeline stress assessment result may include whether the target pipeline is at risk or not. Correspondingly, when the target pipeline is at risk, the target pipeline risk assessment result also includes the risk sampling location of the target pipeline. Here, the risk sampling location refers to the target sampling location on the target pipeline where the risk exists.

[0109] Optionally, the target three-dimensional total stress tensor field can be input into a pre-trained pipeline stress assessment model to assess the pipeline stress state of the target pipeline and output the target pipeline stress assessment results.

[0110] Optionally, the equivalent stress at each target acquisition location on the target pipeline can be calculated based on the target's three-dimensional total stress tensor field. The equivalent stress at each acquisition location is then compared to a preset stress threshold. If a target acquisition location has an equivalent stress greater than the preset stress threshold, the target pipeline is determined to be at risk; otherwise, it is determined not to be at risk. Optionally, when a target pipeline is at risk, the corresponding acquisition location is designated as a risk acquisition location.

[0111] The technical solution of this invention acquires the internal pressure and multi-path ultrasonic time-series signals of the target pipeline simultaneously. By using a preset path sensitivity matrix and a preset internal pressure sensitivity vector, the physical source of the pipeline stress response under complex conditions is first decoupled to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. Then, through field reconstruction technology based on the acoustoelastic principle, a high-resolution and quantifiable target three-dimensional total stress tensor field is generated. The pipeline stress state is evaluated by the target pipeline through the target three-dimensional total stress tensor field. This achieves a deep integration of non-destructive testing, solid mechanics, data analysis, and engineering decision-making, changing the limitations of traditional single-point pipeline inspection and improving the accuracy, comprehensiveness, and engineering practicality of pipeline stress state evaluation.

[0112] Example 2 Figure 2 This is a flowchart of a pipeline stress assessment method provided in Embodiment 2 of the present invention. Based on the above embodiments, this embodiment of the present invention specifies the process of "assessing the pipeline stress state based on the target three-dimensional total stress tensor field and generating target pipeline stress assessment results" as follows: "Based on the target three-dimensional total stress tensor field, calculate the target equivalent stress, target maximum shear stress, and target hydrostatic pressure at each target acquisition location on the target pipeline; screen target acquisition locations where the target equivalent stress is greater than a preset stress threshold to form an abnormal stress region location set; within the abnormal stress region location set, determine the target acquisition locations where the target equivalent stress is a local extreme point of a local abnormal stress region as high-risk locations; for a single high-risk location, calculate the target risk value of the high-risk location based on the target equivalent stress, target maximum shear stress, and target hydrostatic pressure." This achieves rapid and accurate location of potential risk points in the target pipeline, improves the safety of pipeline operation, and enhances the intuitiveness and accuracy of the target pipeline stress assessment results. It should be noted that parts not detailed in this embodiment of the present invention can be found in the descriptions of other embodiments.

[0113] See Figure 2 The pipe stress assessment method shown includes: S201. Obtain the target internal pressure in the target pipeline and the target ultrasonic timing signals of each target detection path.

[0114] S202. Based on the target's internal pressure, the target's ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, stress response analysis is performed on the target pipeline to obtain the target's external load stress response subset and the target's internal pressure stress response subset.

[0115] S203. Based on the target external load stress response subset and the target internal pressure stress response subset, perform three-dimensional stress simulation on the target pipeline to obtain the target three-dimensional total stress tensor field.

[0116] S204. Based on the target's three-dimensional total stress tensor field, calculate the target's equivalent stress, maximum shear stress, and hydrostatic pressure at each target acquisition location on the target pipeline.

[0117] Target equivalent stress is used to characterize the overall intensity of the stress state at the target sampling location from a single numerical dimension. Target equivalent stress eliminates the directionality of the stress tensor. Target equivalent stress is a key parameter for determining whether pipeline materials have entered a plastic state or experienced fatigue failure. By characterizing the overall intensity of the stress state at the target sampling location with a single positive number, target equivalent stress facilitates direct comparison with a preset stress threshold.

[0118] The target maximum shear stress is the maximum shear intensity experienced at the target acquisition location. The target hydrostatic pressure is the average normal stress state experienced at the target acquisition location. The average normal stress includes both tensile and compressive forces.

[0119] Specifically, for a single target acquisition location, the equivalent stress calculation formula is used to synthesize the stress tensor of each component at each target acquisition location on the target pipeline, thereby obtaining the target equivalent stress at the target acquisition location on the target pipeline.

[0120] For example, the equivalent stress calculation formula can be represented by the following formula: ; In the formula, Equivalent stress to the target; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress components along the axial direction of the target pipe; where tensile stress is positive and compressive stress is negative; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress components along the circumferential direction of the target pipeline; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target Shear stress components within the plane of the target pipe surface.

[0121] Specifically, for a single target acquisition location, the maximum shear stress calculation formula is used to synthesize the stress tensor of each component at each target acquisition location on the target pipeline, thereby obtaining the maximum shear stress of the target at the target acquisition location on the target pipeline.

[0122] For example, the formula for calculating the maximum shear stress can be expressed as follows: ; In the formula, The target is the maximum shear stress; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress components along the axial direction of the target pipe; where tensile stress is positive and compressive stress is negative; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress components along the circumferential direction of the target pipeline; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target Shear stress components within the plane of the target pipe surface.

[0123] Specifically, for a single target acquisition location, the hydrostatic pressure calculation formula is used to synthesize the stress tensor of each component at each target acquisition location on the target pipeline to obtain the target hydrostatic pressure at the target acquisition location on the target pipeline.

[0124] For example, the following formula can be used to represent the formula for calculating hydrostatic pressure: ; In the formula, Target hydrostatic pressure; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress components along the axial direction of the target pipe; where tensile stress is positive and compressive stress is negative; The stress tensor at target acquisition location k in the three-dimensional total stress tensor field of the target The normal stress component along the circumferential direction of the target pipeline.

[0125] S205. Collect the locations of each target whose equivalent stress is greater than the preset stress threshold to form a set of abnormal stress area locations.

[0126] A preset stress threshold is a standard value used to characterize the equivalent stress of the target pipeline when it is in a safe state. Optionally, the preset stress threshold can be determined based on the yield strength, allowable stress, safety factor, and / or relevant engineering specifications of the pipeline material. For example, the preset stress threshold can be 0.8 times the yield strength of the material.

[0127] The set of abnormal stress regions is the collection of data acquisition locations of all targets where the equivalent stress exceeds a preset stress threshold. This set of abnormal stress regions can be denoted as... In the formula, Let be the set of locations of abnormal stress regions; k is the acquisition location of the kth target; It is the target equivalent stress at the k-th target acquisition location; It is a preset stress threshold.

[0128] Specifically, the equivalent stress of all target acquisition locations is traversed, and all target acquisition locations that meet the requirement of having an equivalent stress greater than a preset stress threshold are selected to form a set of abnormal stress region locations.

[0129] S206. In areas of concentrated abnormal stress, the target acquisition location of the local extreme point of the target equivalent stress in the local abnormal stress area is determined as a high-risk location.

[0130] The local abnormal stress region is a region defined based on the target acquisition locations within the abnormal stress region location set. The local abnormal stress region is used to screen high-risk locations among the target acquisition locations within the region. The local abnormal stress region is a preset spatial neighborhood of a single target acquisition location within the abnormal stress region location set. For example, it can be defined as circular regions centered on a single target acquisition location within the abnormal stress region location set, with a radius equal to a preset proportion of the pipe diameter. Circular regions with overlapping areas are merged to obtain the local abnormal stress region. The preset proportion is used to adjust the range of the local abnormal stress region. The preset proportion can be set and adjusted according to the actual requirements of pipeline risk detection in practical work.

[0131] A local extreme point is the maximum equivalent stress of a target within a local abnormal stress region. The target acquisition location corresponding to a local extreme point is usually the most severe location within the local abnormal stress region. For example, the target acquisition location corresponding to a local extreme point could be the location where a crack is most likely to initiate. A high-risk location is the target acquisition location corresponding to a local extreme point within a local abnormal stress region. For example, a high-risk location can be denoted as... .

[0132] Specifically, the target acquisition locations within the concentrated areas of abnormal stress are divided into local abnormal stress regions. The equivalent stress of the targets at each acquisition location within these local abnormal stress regions is compared to determine the local extreme points. The target acquisition locations corresponding to these local extreme points are then designated as high-risk locations.

[0133] S207. For a single high-risk location, the target risk value of the high-risk location is calculated based on the target equivalent stress, target maximum shear stress, and target hydrostatic pressure at the high-risk location.

[0134] The target risk value is used to comprehensively assess the risk level of high-risk locations along a target pipeline. The higher the target risk value, the greater the risk level of those high-risk locations.

[0135] Specifically, for a single high-risk location, the target risk value of the high-risk location is calculated using the target risk assessment formula based on the target equivalent stress, target maximum shear stress, and target hydrostatic pressure of the high-risk location.

[0136] For example, the following formula can be used to represent the target risk assessment formula: ; In the formula, It is the target risk value for a high-risk location. It is a scalar greater than 0; These are weighting coefficients with a sum of 1, which can be preset based on the importance of the failure mechanism (such as yielding, shearing, or brittle fracture) or actual working requirements; It is a high-risk location. The target equivalent force at the location; It is the preset yield strength of the pipe material, which can be obtained by consulting the technical specifications of the pipe material or a standard material performance database; It is a high-risk location. The target maximum shear stress at the location; It is a preset shear yield strength, used to characterize the yield limit of pipe materials under pure shear. It can be determined by the yield criterion in mechanics of materials through the following formula. Calculated; It is a high-risk location. The target hydrostatic pressure at the location.

[0137] The technical solution of this invention calculates the target equivalent stress through the target total stress tensor field and quickly filters the location set of abnormal stress regions by applying a preset stress threshold. It also accurately locates the local extreme points of the target equivalent stress within the local abnormal stress region, thus determining high-risk locations. This achieves rapid and accurate location of potential risk points in the target pipeline, improving the safety of pipeline operation. By combining multi-dimensional stress components to calculate the target risk value of high-risk locations, the intuitiveness and accuracy of the target pipeline stress assessment results are improved.

[0138] In an optional embodiment of the present invention, after calculating the target risk value of a single high-risk location based on the target equivalent stress, target maximum shear stress, and target hydrostatic pressure of the high-risk location, the method further includes: obtaining the historical risk value of the single high-risk location, and obtaining the comprehensive risk value of the high-risk location based on the target risk value and the historical risk value.

[0139] Historical risk values ​​characterize the probability of risk existing at high-risk locations identified through analysis of historical maintenance records. The comprehensive risk value is a quantitative indicator that combines the severity of stress conditions and the urgency of risk at high-risk locations. The comprehensive risk value is a scalar.

[0140] Optionally, historical risk values ​​for high-risk locations can be assigned according to preset historical risk value generation rules and historical maintenance records. These preset rules are used to assign historical risk values ​​to high-risk locations. For example, preset historical risk value generation rules include 1.0 for no historical records, 1.5 for records of no repairs, and 0.8 for locations repaired within the last 60 days. "No historical records" can be understood as no record of the high-risk location in historical maintenance records, meaning no problem was detected at the high-risk location during historical maintenance. "Records of no repairs" can be understood as a record of the high-risk location in historical maintenance records, but after review, no repair was required, meaning the high-risk location was a false detection during historical maintenance. "Repaired within the last 60 days" can be understood as a high-risk location being detected as having a problem and being repaired during historical maintenance within the last 60 days.

[0141] Correspondingly, for a single high-risk location, the historical risk value of the high-risk location is determined based on historical maintenance records according to the preset historical risk value generation rules.

[0142] Optionally, a preset pipeline historical maintenance database can store maintenance records for the entire pipeline lifecycle and historical risk values ​​for each high-risk location on the pipeline.

[0143] Correspondingly, for a single high-risk location, the historical maintenance database of the pipeline can be queried to determine the historical risk value of the high-risk location.

[0144] Specifically, a comprehensive risk assessment formula is used to combine the target risk value and historical risk value to obtain the comprehensive risk value of high-risk locations.

[0145] For example, the comprehensive risk assessment formula can be represented by the following formula: ; In the formula, It is the comprehensive risk value of the high-risk location k; It is the target risk value at the high-risk location k; It is the historical risk value of high-risk location k; where k is the index of any high-risk location.

[0146] This solution calculates a comprehensive risk value by integrating the target risk value obtained from real-time mechanical analysis with the historical risk value determined based on historical maintenance data. By taking into account the historical risk value in addition to the target risk value, the comprehensive risk value obtained can more comprehensively reflect the risk status of the pipeline.

[0147] In an optional embodiment of the present invention, after obtaining the comprehensive risk value of high-risk locations based on the target risk value and historical risk values, the method further includes: spatially clustering each high-risk location to obtain each intervention pipeline segment; and sorting each intervention pipeline segment based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment to obtain a target intervention pipeline segment sequence.

[0148] The target intervention pipeline segment sequence is an ordered sequence of descriptions of intervention pipeline segments. It characterizes the maintenance order of each intervention pipeline segment. In essence, the earlier an intervention pipeline segment appears in the target sequence, the higher its maintenance priority. The target intervention pipeline segment sequence provides maintenance personnel with clear, ordered, and spatially aggregated maintenance guidance. Optionally, the target intervention pipeline segment sequence can be generated by spatially clustering and prioritizing high-risk locations. For example, the mathematical form of the target intervention pipeline segment sequence can be an ordered sequence of segment descriptions, which can be denoted as... A single intervention pipeline segment can be denoted as... .

[0149] An intervention pipeline segment is a single pipeline segment within a sequence of target intervention pipeline segments. An intervention pipeline segment represents a continuous pipeline segment that requires intervention. A single intervention pipeline segment defines a continuous physical interval on the target pipeline that needs maintenance. For example, an intervention pipeline segment includes a start coordinate point and an end coordinate point.

[0150] Specifically, a density-based spatial clustering algorithm can be used to perform cluster analysis on each high-risk location. The density-based spatial clustering algorithm automatically divides adjacent high-risk locations into a "cluster" or "group" based on preset parameters such as the spatial distance between high-risk locations and the minimum number of points in the neighborhood. Based on the high-risk locations contained in a single "cluster" or "group", the starting and ending coordinate points are determined to obtain the intervention pipeline segment.

[0151] Specifically, for a single intervention pipeline segment, the comprehensive risk values ​​of all high-risk locations within that segment are compared to determine the highest comprehensive risk value for that segment. The highest comprehensive risk values ​​of each intervention pipeline segment can then be sorted from highest to lowest to generate a sequence of target intervention pipeline segments arranged by priority from highest to lowest.

[0152] This scheme intelligently categorizes discrete high-risk locations into an ordered sequence of target intervention pipeline segments, avoiding isolated treatment of individual risk points and thus utilizing maintenance resources more efficiently. By spatially clustering each high-risk location, intervention pipeline segments are obtained, linking local risk points with the overall pipeline structure and achieving pipeline risk assessment from a global perspective. By ranking each intervention pipeline segment based on the comprehensive risk value of each high-risk location contained within it, a target intervention pipeline segment sequence is obtained. Prioritizing intervention pipeline segments with higher risk levels can eliminate potential fault points in advance, reducing the probability of pipeline leaks and ruptures, and further improving the safety of pipeline operation.

[0153] In an optional embodiment of the present invention, after sorting each intervention pipeline segment based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment to obtain a target intervention pipeline segment sequence, the method further includes: generating a target high-risk point list based on each high-risk location and the corresponding comprehensive risk value; and generating a target pipeline stress assessment report based on the target three-dimensional total stress tensor field, the target high-risk point list, and the target intervention pipeline segment sequence.

[0154] The target high-risk point list is used to comprehensively represent each high-risk location and its corresponding comprehensive risk value. For example, the target high-risk point list can be an ordered list of tuples representing each high-risk location. The target high-risk point list can be denoted as... Each element in the list of high-risk target points corresponds to an identified high-risk location. These are the three-dimensional spatial coordinates of the high-risk location q; It is the comprehensive risk value of the high-risk position q.

[0155] The target pipeline stress assessment report is a structured language result comprising the target's three-dimensional total stress tensor field, a list of high-risk points, and a sequence of pipeline segments requiring intervention. The report is used to synthesize the target's three-dimensional total stress tensor field, the list of high-risk points, and the sequence of pipeline segments requiring intervention.

[0156] Specifically, the high-risk locations and their corresponding comprehensive risk values ​​are concatenated to obtain a list of target high-risk points. Based on preset language organization rules, the target three-dimensional total stress tensor field, the list of target high-risk points, and the target intervention pipeline segment sequence are converted into structured language to obtain a target pipeline stress assessment report. The preset language organization rules are used to convert the target three-dimensional total stress tensor field, the list of target high-risk points, and the target intervention pipeline segment sequence into structured language. Optionally, the preset language organization rules can be set and adjusted according to the format requirements of the pipeline stress assessment report and the various data formats in actual work. This solution does not further limit the setting of the preset language organization rules.

[0157] This solution ultimately generates a target pipeline stress assessment report that integrates visualization, quantification, and decision-making recommendations. This makes the target pipeline stress assessment results more structured and traceable, and realizes a closed-loop transformation from raw stress data to directly driven maintenance decisions. It deeply integrates detection technology, mechanical analysis, and big data decision support, thereby improving the scientific nature, predictability, and efficiency of pipeline risk management.

[0158] Example 3 Figure 3 This is a schematic diagram of a pipeline stress assessment device provided in Embodiment 3 of the present invention. This embodiment of the present invention is applicable to situations involving pipeline stress assessment. The device can execute pipeline stress assessment methods and can be implemented in hardware and / or software. The device can be configured in electronic devices that perform pipeline stress assessment functions, such as detection terminals and / or servers.

[0159] See Figure 3 The pipeline stress assessment device shown includes: a data acquisition module 301, a data analysis module 302, a pipeline stress synthesis module 303, and a pipeline stress assessment module 304. Specifically, the data acquisition module 301 acquires the target internal pressure of the target pipeline and the target ultrasonic time-series signals of each target detection path; the data analysis module 302 performs stress response analysis on the target pipeline based on the target internal pressure, the target ultrasonic time-series signals, a preset internal pressure sensitivity vector, and a preset path sensitivity matrix, obtaining a subset of target external load stress response and a subset of target internal pressure stress response; the pipeline stress synthesis module 303 performs three-dimensional stress simulation on the target pipeline based on the subsets of target external load stress response and the subsets of target internal pressure stress response, obtaining a target three-dimensional total stress tensor field; and the pipeline stress assessment module 304 assesses the pipeline stress state of the target pipeline based on the target three-dimensional total stress tensor field, generating a target pipeline stress assessment result.

[0160] The technical solution of this invention acquires the internal pressure and multi-path ultrasonic time-series signals of the target pipeline simultaneously. By using a preset path sensitivity matrix and a preset internal pressure sensitivity vector, the physical source of the pipeline stress response under complex conditions is first decoupled to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. Then, through field reconstruction technology based on the acoustoelastic principle, a high-resolution and quantifiable target three-dimensional total stress tensor field is generated. The pipeline stress state is evaluated by the target pipeline through the target three-dimensional total stress tensor field. This achieves a deep integration of non-destructive testing, solid mechanics, data analysis, and engineering decision-making, changing the limitations of traditional single-point pipeline inspection and improving the accuracy, comprehensiveness, and engineering practicality of pipeline stress state evaluation.

[0161] In an optional embodiment of the present invention, the apparatus further includes: a standard external load acquisition module, used to acquire target pipeline structural parameters, a set of standard external loads, and a target detection path; a first theoretical stress distribution field generation module, used to perform stress field simulation on the target pipeline based on the target pipeline structural parameters for a single standard external load, to obtain a first theoretical stress distribution field of the target pipeline under the standard external load; a first theoretical time change calculation module, used to integrate the first theoretical stress distribution field along the target detection path for a single standard external load using a preset set of acoustoelastic physics formulas, to obtain a first theoretical time change of the target ultrasonic time-series signal under the standard external load; and a preset path sensitivity matrix generation module, used to generate a preset path sensitivity matrix based on the first theoretical time change of each target ultrasonic time-series signal under each standard external load.

[0162] In an optional embodiment of the present invention, the data analysis module 302 includes: a target time change vector calculation unit, used to calculate the target time change vector based on the target ultrasonic time-series signal and the reference signal; a target external load coefficient vector estimation unit, used to construct a target stress response separation equation based on the target time change vector, the target internal pressure, a preset internal pressure sensitivity vector, and a preset path sensitivity matrix, and to estimate the target external load coefficient vector based on the target stress response separation equation; an external load stress response calculation unit, used to calculate a subset of the target external load stress response based on the target external load coefficient vector and the preset path sensitivity matrix; and an internal pressure stress response calculation unit, used to calculate a subset of the target internal pressure stress response based on the target internal pressure and the preset internal pressure sensitivity vector.

[0163] In an optional embodiment of the present invention, the pipeline stress synthesis module 303 includes: a target stress tensor equation construction unit, used to construct a target stress tensor equation for a single target detection path based on a subset of the target external load stress response, a subset of the target internal pressure stress response, and preset acoustoelastic physical coefficients corresponding to each target acquisition position of the target detection path; a target stress tensor equation solving unit, used to solve the target stress tensor equations of each target detection path to obtain the target external load stress distribution vector and the target internal pressure stress distribution vector at each target acquisition position; and a pressure stress distribution vector superposition unit, used to superimpose the target external load stress distribution vector and the target internal pressure stress distribution vector at each target acquisition position to obtain a target three-dimensional total stress tensor field.

[0164] In an optional embodiment of the present invention, the pipeline stress assessment module 304 includes: a multi-dimensional stress calculation unit, used to calculate the target equivalent stress, target maximum shear stress, and target hydrostatic pressure at each target acquisition location on the target pipeline based on the target three-dimensional total stress tensor field; an abnormal stress region location set construction unit, used to screen each of the target acquisition locations where the target equivalent stress is greater than a preset stress threshold to form an abnormal stress region location set; a high-risk location screening unit, used to identify the target acquisition locations where the target equivalent stress is a local extreme point of a local abnormal stress region as high-risk locations within the abnormal stress region location set; and a target risk value calculation unit, used to calculate the target risk value of a single high-risk location based on the target equivalent stress, the target maximum shear stress, and the target hydrostatic pressure of the high-risk location.

[0165] In an optional embodiment of the present invention, the pipeline stress assessment module 304 further includes: a comprehensive risk value calculation unit, used to calculate the target risk value of a high-risk location based on the target equivalent stress, the target maximum shear stress, and the target hydrostatic pressure of the high-risk location for a single high-risk location, and then obtain the historical risk value of the high-risk location for a single high-risk location, and obtain the comprehensive risk value of the high-risk location based on the target risk value and the historical risk value.

[0166] In an optional embodiment of the present invention, the pipeline stress assessment module 304 further includes: an intervention pipeline segment determination unit, configured to perform spatial clustering on each of the high-risk locations to obtain each intervention pipeline segment after obtaining the comprehensive risk value of the high-risk location based on the target risk value and the historical risk value; and an intervention pipeline segment sequence generation unit, configured to sort each intervention pipeline segment based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment to obtain a target intervention pipeline segment sequence.

[0167] In an optional embodiment of the present invention, the pipeline stress assessment module 304 further includes: a target high-risk point list generation unit, used to sort each intervention pipeline segment based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment to obtain a target intervention pipeline segment sequence, and then generate a target high-risk point list based on each high-risk location and the corresponding comprehensive risk value; and a target pipeline stress assessment report generation unit, used to generate a target pipeline stress assessment report based on the target three-dimensional total stress tensor field, the target high-risk point list, and the target intervention pipeline segment sequence.

[0168] The pipeline stress assessment device provided in this embodiment of the invention can execute the pipeline stress assessment method provided in any embodiment of the invention, and has the corresponding functional modules and beneficial effects of the method.

[0169] In the technical solutions of this invention, the acquisition, storage, and application of the target internal pressure of the target pipeline, the target ultrasonic timing signals of each target detection path, the target pipeline structural parameters, the standard external load set, and the historical risk values ​​of the target detection path and high-risk locations all comply with the provisions of relevant laws and regulations and do not violate public order and good morals.

[0170] Example 4 Figure 4 A schematic diagram of an electronic device 400 that can be used to implement embodiments of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices (such as helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0171] like Figure 4 As shown, the electronic device 400 includes at least one processor 401 and a memory, such as a read-only memory (ROM) 402 or a random access memory (RAM) 403, communicatively connected to the at least one processor 401. The memory stores computer programs executable by the at least one processor. The processor 401 can perform various appropriate actions and processes based on the computer program stored in the ROM 402 or loaded into the RAM 403 from storage unit 408. The RAM 403 may also store various programs and data required for the operation of the electronic device 400. The processor 401, ROM 402, and RAM 403 are interconnected via a bus 404. An input / output (I / O) interface 405 is also connected to the bus 404.

[0172] Multiple components in electronic device 400 are connected to I / O interface 405, including: input unit 406, such as keyboard, mouse, etc.; output unit 407, such as various types of displays, speakers, etc.; storage unit 408, such as disk, optical disk, etc.; and communication unit 409, such as network card, modem, wireless transceiver, etc. Communication unit 409 allows electronic device 400 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0173] Processor 401 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 401 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 401 performs the various methods and processes described above, such as pipe stress assessment methods.

[0174] In some embodiments, the pipe stress assessment method may be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 408. In some embodiments, part or all of the computer program may be loaded and / or installed on electronic device 400 via ROM 402 and / or communication unit 409. When the computer program is loaded into RAM 403 and executed by processor 401, one or more steps of the pipe stress assessment method described above may be performed. Alternatively, in other embodiments, processor 401 may be configured to perform the pipe stress assessment method by any other suitable means (e.g., by means of firmware).

[0175] Various implementations of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), complex programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various implementations may include: implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0176] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0177] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0178] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0179] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or middleware components (e.g., application servers), or frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0180] A computing system can include clients and servers. Clients and servers are generally geographically separated and typically interact via communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system. It addresses the shortcomings of traditional physical hosts and VPS (Virtual Private Server) services, such as high management difficulty and weak business scalability.

[0181] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0182] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for assessing pipeline stress, characterized in that, The method includes: Acquire the internal pressure of the target pipeline and the ultrasonic timing signals of each target detection path; Based on the target internal pressure, the target ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, stress response analysis is performed on the target pipeline to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. Based on the target external load stress response subset and the target internal pressure stress response subset, a three-dimensional stress simulation is performed on the target pipeline to obtain the target three-dimensional total stress tensor field. Based on the target three-dimensional total stress tensor field, the target pipeline stress state is evaluated, and the target pipeline stress evaluation result is generated.

2. The method according to claim 1, characterized in that, The process of generating the preset path sensitivity matrix includes: Obtain the target pipeline structural parameters, standard external load set, and target detection path; For a single standard external load, based on the target pipeline structural parameters, stress field simulation is performed on the target pipeline to obtain the first theoretical stress distribution field of the target pipeline under the standard external load; For a single standard external load, a preset set of acoustoelastic physical formulas is used to integrate the first theoretical stress distribution field along the target detection path to obtain the first theoretical time change of the target ultrasonic time sequence signal under the standard external load. Based on the first theoretical time change of each of the target ultrasonic time-series signals under each of the standard external loads, a preset path sensitivity matrix is ​​generated.

3. The method according to claim 1, characterized in that, The stress response analysis of the target pipeline is performed based on the target internal pressure, the target ultrasonic time-series signal, a preset internal pressure sensitivity vector, and a preset path sensitivity matrix to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response, including: Based on the target ultrasonic time-series signal and the reference reference signal, the target time change vector is calculated; Based on the target time change vector, the target internal pressure, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, a target stress response separation equation is constructed, and based on the target stress response separation equation, the target external load coefficient vector is estimated. Based on the target external load coefficient vector and the preset path sensitivity matrix, calculate the target external load stress response subset; Based on the target internal pressure and the preset internal pressure sensitivity vector, calculate the target internal pressure stress response subset.

4. The method according to claim 1, characterized in that, The step of performing a three-dimensional stress simulation on the target pipeline based on the target external load stress response subset and the target internal pressure stress response subset to obtain the target three-dimensional total stress tensor field includes: For a single target detection path, based on the target's external load stress response subset, the target's internal pressure stress response subset, and the preset acoustoelastic physical coefficients corresponding to each target acquisition position on the target detection path, the target stress tensor equation of the target detection path is constructed. Solve the target stress tensor equations for each of the target detection paths to obtain the target external load stress distribution vector and the target internal pressure stress distribution vector at each target acquisition location; The target external load stress distribution vector and the target internal pressure stress distribution vector at each target acquisition location are superimposed to obtain the target three-dimensional total stress tensor field.

5. The method according to claim 1, characterized in that, The step of assessing the pipeline stress state based on the target three-dimensional total stress tensor field and generating target pipeline stress assessment results includes: Based on the target's three-dimensional total stress tensor field, calculate the target's equivalent stress, maximum shear stress, and hydrostatic pressure at each target acquisition location on the target pipeline. The target acquisition locations that are greater than a preset stress threshold are selected to form an abnormal stress region location set; The locations of the targets whose equivalent stress is a local extreme point in the local abnormal stress region are identified as high-risk locations. For a single high-risk location, the target risk value of the high-risk location is calculated based on the target equivalent stress, the target maximum shear stress, and the target hydrostatic pressure at the high-risk location.

6. The method according to claim 5, characterized in that, After calculating the target risk value of a single high-risk location based on the target equivalent stress, the target maximum shear stress, and the target hydrostatic pressure at that location, the method further includes: For a single high-risk location, the historical risk value of the high-risk location is obtained, and based on the target risk value and the historical risk value, the comprehensive risk value of the high-risk location is obtained.

7. The method according to claim 6, characterized in that, After obtaining the comprehensive risk value of the high-risk location based on the target risk value and the historical risk value, the method further includes: Spatial clustering is performed on each of the high-risk locations to obtain each intervention pipeline segment; Based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment, the intervention pipeline segments are sorted to obtain the target intervention pipeline segment sequence.

8. The method according to claim 7, characterized in that, After sorting the intervention pipeline segments based on the comprehensive risk value of each high-risk location contained in each intervention pipeline segment to obtain the target intervention pipeline segment sequence, the method further includes: Based on each of the high-risk locations and their corresponding comprehensive risk values, a list of target high-risk points is generated. Based on the target's three-dimensional total stress tensor field, the target's list of high-risk points, and the target's intervention pipeline segment sequence, a target pipeline stress assessment report is generated.

9. A pipe stress assessment device, characterized in that, The device includes: The data acquisition module is used to acquire the internal pressure of the target pipeline and the ultrasonic time-series signals of each target detection path; The data analysis module is used to perform stress response analysis on the target pipeline based on the target internal pressure, the target ultrasonic time-series signal, the preset internal pressure sensitivity vector, and the preset path sensitivity matrix, to obtain a subset of the target external load stress response and a subset of the target internal pressure stress response. The pipeline stress synthesis module is used to perform three-dimensional stress simulation on the target pipeline based on the target external load stress response subset and the target internal pressure stress response subset to obtain the target three-dimensional total stress tensor field. The pipeline stress assessment module is used to assess the pipeline stress state of the target pipeline based on the target three-dimensional total stress tensor field and generate the target pipeline stress assessment result.

10. An electronic device, characterized in that, The electronic device includes: At least one processor; and A memory communicatively connected to the at least one processor; wherein, The memory stores a computer program that can be executed by the at least one processor to enable the at least one processor to perform the pipe stress assessment method according to any one of claims 1-8.

Images