A method for predicting the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes based on a deep operator network architecture

Abstract

The present invention provides a method for predicting the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes based on a deep operator network architecture. First, a simulated usage environment is constructed, and multi-batch ultimate stress tests are carried out on polyurethane directly buried pipes of different specifications under multiple environmental conditions. The corresponding specification and environmental data are recorded as inputs, and the maximum stress is used as the output to construct a data set. A feature screening module is designed to eliminate features with low relevance. Then, an ultimate state prediction model is designed and constructed based on the deep operator network architecture. A reinforcement learning strategy is used to design the optimization training for the model. The trained prediction model is deployed to the relevant environment. By inputting the specification of the thermal insulation pipe and environmental parameters on the computer terminal, the ultimate stress situation can be obtained. The present invention can eliminate the redundant features introduced in the initial features, further optimize the model architecture, reduce the complexity of model training, and ensure that the prediction model design can accurately predict the ultimate stress of polyurethane directly buried thermal insulation pipes under various working conditions.

Description

Technical Field

[0001] The present invention belongs to the technical field of predicting the ultimate state of pre-insulated pipes based on machine learning, and particularly relates to a method for predicting the ultimate state of prefabricated polyurethane directly buried insulation pipes based on a deep operator network architecture. Background Art

[0002] Polyurethane directly buried insulation pipes have excellent heat insulation performance, which can greatly reduce the heat loss during transportation, effectively save energy, and can be achieved whether it is for hot water, steam transportation or maintaining a low-temperature environment. Their characteristics of waterproofing, moisture-proofing and anti-corrosion protection can protect the pipes from external erosion, extend the service life and reduce the maintenance cost. At the same time, they are convenient for construction, adaptable to complex terrains, quick to install, and also have a certain structural support capacity to withstand various external forces. In addition, the flame retardant performance and alarm function of some products provide guarantees for safe operation, meet the requirements of the energy-saving and environmental protection era, and play an irreplaceable role in the pipeline transportation systems in many fields. However, in the actual operation of polyurethane directly buried insulation pipes, they will be affected by external forces such as soil pressure and vehicle loads, as well as thermal stresses caused by the transportation of hot media. Stress analysis can accurately calculate the stress distribution, avoid structural failures such as pipeline rupture, deformation and interface loosening caused by excessive stress, and ensure the structural safety of the pipeline; at the same time, it can also prevent the stress from being transmitted to the insulation layer, avoid damage to the insulation layer, and maintain the insulation performance and structural integrity. In addition, through stress analysis, the design parameters can be optimized, the fatigue life can be predicted, the service life of the pipeline can be extended, the maintenance and replacement costs can be reduced, and the requirements of engineering specifications and standards can be met, ensuring public safety and ensuring the smooth acceptance and commissioning of the project. Therefore, ultimate stress analysis is of great significance in the specific use process of polyurethane insulation pipes.

[0003] Currently, the analysis and prediction of the ultimate state of polyurethane directly buried pipes mainly focus on the following methods:

[0004] (1) Analytical method: The analytical method is a method based on theories such as elastic mechanics and material mechanics, which solves stresses by establishing accurate mathematical models and using mathematical formulas and derivations. For some polyurethane directly buried insulation pipe structures with simple geometric shapes and clear boundary conditions, this method can be used to directly obtain the analytical expression of stress. For example, for a straight pipe with a constant cross-section under uniform temperature change and axial force, according to the theory of thermoelasticity, by establishing equilibrium equations, geometric equations and physical equations, the stress distribution formula in the pipe can be derived, which can accurately give the stress values at each point of the pipe and provide a theoretical basis for stress analysis. However, for cases with complex shapes and boundary conditions, the analytical method is difficult to solve or even impossible to solve;

[0005] (2) Numerical analysis method: The numerical analysis method is to discretize the actual structure of the polyurethane direct-buried thermal insulation pipe into a finite number of units and nodes with the help of computer technology, perform mechanical analysis on each unit, and then solve the stress by assembling the overall equation group. Among them, the finite element method is one of the most commonly used numerical analysis methods. It divides the pipeline structure into a large number of small units, such as tetrahedral units, hexahedral units, etc., establishes a stiffness matrix for each unit, and then solves the displacement field of the entire structure according to the boundary conditions and load conditions, and then calculates the stress field through physical equations. This method can handle various complex geometric shapes, material properties and boundary conditions, and can accurately simulate the stress distribution of the pipeline under different working conditions. It is widely used in the stress analysis of polyurethane direct-buried thermal insulation pipes, but it requires professional software and certain computing resources, and has high requirements for model establishment and parameter setting.

[0006] (3) Experimental stress analysis method: The experimental stress analysis method is a method that directly measures stress or physical quantities related to stress by conducting experiments on actual polyurethane direct-buried insulated pipes or their models, and then determines the stress distribution. A common method is the resistance strain gauge method, in which a resistance strain gauge is attached to the surface of the pipe. When the pipe is deformed by force, the resistance value of the strain gauge changes. By measuring the resistance change and calculating the strain on the pipe surface based on the relationship between strain and resistance, the stress is then obtained based on the constitutive relationship of the material. There is also the photoelastic method, in which a pipe model is made of photoelastic material. Under loading conditions, the fringe pattern that appears in the photoelastic model is observed, and the stress distribution in the model is analyzed and calculated to infer the stress state of the actual pipe. The experimental stress analysis method can directly obtain the stress data of the pipe. The results are intuitive and reliable, and can be used to verify the results of numerical analysis and analytical methods. However, the experimental cost is high, and due to the limitations of experimental conditions and measurement technology, there are certain difficulties in measuring stress under some internal stresses and complex working conditions. Summary of the Invention

[0007] In response to the above problems, the first aspect of the present invention proposes a method for constructing a limit state prediction model for prefabricated polyurethane direct-buried insulated pipes based on a deep operator network architecture, comprising the following steps:

[0008] Step 1: Based on the use environment simulation, carry out extreme stress tests on direct buried pipes of different specifications to collect the maximum stress they can withstand under such conditions. The maximum stress and corresponding specification data and environmental data are recorded to construct an initial data set.

[0009] Step 2: Based on the feature screening module of the relief algorithm, the features with low correlation with the ultimate stress of the polyurethane direct-buried thermal insulation pipe are eliminated, and the specification data and environmental data with the highest feature importance are obtained. The training data set is constructed based on the corresponding maximum stress.

[0010] Step 3: Construct a limit state prediction model based on the deep operator network architecture, including a backbone network architecture, a branch network architecture, and a fusion output module. The backbone network is used to extract the physical property features of the directly buried thermal insulation pipe. The branch network is used to extract the environmental features of the directly buried thermal insulation pipe and perform weighted processing on the feature channels in the environmental features based on the self-attention mechanism to obtain self-attention environmental features. The fusion output module is used to adaptively fuse the two types of features output by the two network architectures and output the maximum stress prediction result based on the fused features;

[0011] Step 4: Optimize and train the constructed limit state prediction model based on the reinforcement learning strategy to obtain the final limit state prediction model.

[0012] Preferably, the specific data processing process of the feature screening module based on the relieff algorithm in Step 2 is as follows:

[0013] S21. Data selection and classification: Based on the initial dataset Data of the limit state of the polyurethane directly buried thermal insulation pipe, randomly select M sample data from it, and group the M sample data according to the data similarity SDU. The number of groups is L;

[0014] S22. Feature difference calculation: Calculate the feature difference TZC corresponding to each feature in the same group of samples and the feature difference BZC in different groups of samples respectively;

[0015] S23. Feature correlation calculation: Calculate the feature correlation Tmax based on the calculated difference TZC and BZC of each feature;

[0016] S24. Iterative calculation of feature correlation multiple times: Repeat the process from S21 to S23 for a total of B iterations. Perform average weighted processing on the feature correlation t∈[1, B], and define its result TO as the final feature importance. The higher the feature importance, the stronger the connection degree between this feature and the ultimate stress of the polyurethane directly buried thermal insulation pipe;

[0017] S25. Feature screening process: Based on the calculated feature importance of each feature, discard the redundant features with the feature importance in the last 40%.

[0018] Preferably, in S21,

[0019]

[0020] where Datayz j represents the j-th parameter value of a piece of data in the dataset Data, It represents the j-th parameter value of another piece of data in the dataset Data, where j ∈ [1, C]; here, C is the total number of environmental data, specification data, and maximum stress data.

[0021] The specific process of grouping is as follows: Randomly select one piece of data from the M pieces of sample data to start grouping. Then, calculate the SDU values of all the remaining data and this piece of data in turn. After the calculation, combine the piece of data with the largest SDU value with this piece of data to form a data group; After the grouped data is removed, continue to perform the above-mentioned data classification on the remaining data until the data grouping is completed;

[0022] In S22, calculate the feature difference degree TZC corresponding to each feature in the same-group samples and the feature difference degree BZC in different-group samples respectively. The formulas are as follows:

[0023]

[0024] Among them, It is the parameter value of the k-th feature of the w-th and m-th data in the same-group sample data, where T is the sum of environmental data and specification data; and respectively represent the parameter values of the k-th feature of the g-th and l-th data in different-group sample data, |*| represents the absolute value operation function.

[0025] Preferably, the backbone network architecture consists of three parts: an input layer, an intermediate layer, and an output layer; Select the specification characteristic parameters of the polyurethane directly buried insulation pipe as the input of this backbone network architecture. After feature screening, the number of retained polyurethane directly buried insulation pipe specification characteristic parameters is p. Therefore, the number of neurons in the input layer is p; The network architecture of the intermediate layer is mainly composed of three parts: a fully connected layer, a Relu activation function, and a Silu activation function. In the intermediate layer, there are first Wq fully connected layers, with Ws neurons in each fully connected layer. Then, it is activated by the Relu function to enhance the non-linear expression ability of the model. Finally, there are Wh fully connected layers, with Wf neurons in each layer of each fully connected layer and the Silu activation function; The designed output layer contains 4 different neurons. Then, through the spatial attention layer, the outputs of the 4 neurons in the hidden layer are adaptively fused to obtain the output Feaz of the backbone network architecture.

[0026] Preferably, a stacking structure is adopted in the branch network architecture; the stacking structure is composed of two parallel network paths, and each network path is composed of an input layer, an intermediate layer, and an output layer; the environmental characteristic parameters are used as the input of this branch network architecture. There are a total of q neurons in the input layer, and q is the total number of environmental characteristic data retained after screening; the structure of the intermediate layer is as follows: first, there are three sequentially connected fully connected layers, and then the data is non-linearly activated through the simgoid activation function. Subsequently, there are four sequentially connected fully connected layers. In these four fully connected layers, finally, non-linear normalization processing is performed through the Softmax activation function; the number of neurons in the output layer of the network path is set to 1; the output FZS of the output layer in each network path is sent to the self-attention mechanism unit for processing to obtain the output result Fe, and then a fully connected layer is used to further integrate the outputs of the two network paths to obtain the output Feaf of the branch network architecture.

[0027] Preferably, the self-attention mechanism unit maps the output FZS of each network path in the branch network architecture to multiple feature spaces. The self-attention mechanism includes L parallel self-attention taps, and then the channel attention layer is used to integrate the outputs of the L attention taps to obtain the output result Fe of the final self-attention mechanism; the feature processing process of each self-attention tap is as follows:

[0028]

[0029] Among them, Softmax(*) represents the Softmax linear normalization function, and Fea is the output result of the self-attention tap; the matrices Q, K, and V are the query, key, and value matrices mapped after the feature FZS is processed by three different matrices. Specifically:

[0030] Q = W Q *FZS

[0031] K = W K *FZS

[0032] V = W V *FZS

[0033] Among them, W Q ,W K ,W V are three different feature mapping matrices to process the feature FZS in three different feature spaces.

[0034] Preferably, the fusion output module further fuses the output Feaz obtained from the backbone network architecture and the output Feaf obtained from the branch network architecture to obtain the predicted result Yout of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe. First, the BN normalization layer is used to normalize Feaz and Feaf respectively to optimize the distribution of feature channels and enhance the robustness of the model when outputting the predicted result of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe, obtaining and Secondly, the bilinear interpolation upsampling layer is used to upsample to the same size; Thirdly, the dot product layer is used to perform a dot product operation on the upsampled and to obtain the fused feature Ffz; Thirdly, the fused feature is input into the Dropout layer to optimize the feature channel parameters in Ffz to obtain the feature Ffz* with frozen parameters, reducing the risk of model overfitting; Finally, Ffz* is input into the L2 regularization layer and then into the Relu activation function to obtain the final predicted result Yout of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe, where the L2 regularization layer constrains the model weights and enhances the robustness of the model during prediction.

[0035] Preferably, the reinforcement learning strategy comprehensively optimizes the number of fully connected layers in the intermediate layer of the model backbone network and the number of neurons Ws and Wf contained in each fully connected layer Wq and Wh to find the specific number of fully connected layers and the number of neurons in each layer corresponding to the best model performance;

[0036] The vector composed of the four quantity combinations of the number of fully connected layers Wq, Wh and the corresponding number of neurons Ws, Wf is used as the state space corresponding to this reinforcement learning; in its action space, it adaptively optimizes and changes the parameters Wq, Wh, Ws, Wf to find the model parameters corresponding to the best model performance; the design of its reward factor is realized by comparing the magnitudes of the loss functions. Using the constructed training dataset, based on the model parameters corresponding to the state space, the loss function value Lossq of the model at this time is obtained; when the state is updated, that is, an action to update the model parameters is taken, and then the loss function value Lossh corresponding to the model after the parameter update is calculated; if the value of Qjl is greater than 0, a positive reward value of 2*Qjl is given to this action, and if the value of Qjl is less than 0, a negative reward value of 2*Qjl is given to this action:

[0037] Qjl = Lossh - Lossq

[0038] Meanwhile, a reinforcement learning network is constructed. A policy network is constructed to optimize the parameters Wq, Wh, Ws, and Wf, and a value network is constructed to adaptively judge the value corresponding to this set of model parameters. The designed policy network mainly consists of 4 interconnected fully connected layers and a Relu activation function, and the designed value network mainly consists of 6 interconnected fully connected layers and a Silu activation function. The input of its network architecture is a set of four-dimensional vectors, representing a corresponding set of Wq, Wh, Ws, and Wf situations. The constructed dataset is used and the prediction model constructed with the input of this architecture at this time is trained, and the loss function after the training process is recorded to calculate the reward factor.

[0039] In the second aspect of the present invention, a method for predicting the ultimate state of prefabricated polyurethane directly buried insulation pipes based on a deep operator network architecture is proposed. The ultimate state prediction model constructed by the construction method described in the first aspect is deployed in a computer terminal, and the following process is included:

[0040] The specification data and environmental data of the polyurethane directly buried insulation pipe are obtained in real time and standardized into model input data;

[0041] The standardized specification data and environmental data are input into the ultimate state prediction model; the ultimate stress result of the polyurethane directly buried insulation pipe of this specification under this environmental condition is output.

[0042] Preferably, the specification characteristic parameters of the polyurethane directly buried insulation pipe retained for prediction after screening are: the steel pipe layer thickness GT of the insulation pipe, the insulation layer thickness WT of the insulation pipe, the diameter BD of the insulation pipe, the compressive strength BQ of the insulation pipe, the thermal conductivity BR of the insulation pipe, and the water absorption rate BH of the insulation pipe; the retained environmental characteristic parameters are the environmental temperature TM, the buried depth SD of the directly buried pipe, the surrounding soil density TD, and the soil friction coefficient TS.

[0043] Compared with the prior art, the present invention has the following beneficial effects:

[0044] (1) Feature screening module based on the relieff algorithm: A feature screening module is designed based on the relieff algorithm. By screening the importance of features in the dataset, only the features with a relatively high correlation with the ultimate stress are retained. According to the design of this feature screening module, the redundant features introduced in the initial features can be eliminated, thereby further optimizing the model architecture, reducing the complexity of model training, enhancing the deep fitting ability of the model, optimizing the complexity of feature selection when predicting the ultimate state of the polyurethane directly buried insulation pipe, and enhancing the prediction ability of the model;

[0045] (2) Establishing a limit state prediction model for polyurethane direct-buried thermal insulation pipes based on a deep operator network architecture: Based on the deep fitting capability of the deep operator network architecture, the present invention designs a limit state prediction model for polyurethane direct-buried thermal insulation pipes based on this architecture, and introduces a self-attention mechanism to further enhance the model's generalization and deep fitting capabilities. This ensures that the ultimate stress of polyurethane direct-buried thermal insulation pipes can be accurately predicted under various working conditions.

[0046] (3) Construction of a reinforcement learning-based model optimization module: A model parameter architecture optimization module was constructed based on a reinforcement learning strategy. The loss function, state, action, and reward factor were designed to find the model architecture parameters that achieve the best model performance. Based on this module design, the model fitting performance can be further optimized, enhancing the accuracy and robustness of the ultimate stress prediction of polyurethane direct-buried insulation pipes. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions of the present invention or the prior art, a brief introduction will be given below to the drawings required for use in the embodiments or the description of the prior art. Obviously, what is described below is only one embodiment of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0048] Figure 1 It is the logic block diagram of the overall process of the present invention;

[0049] Figure 2 This is a flow chart of the feature screening operation of the present invention;

[0050] Figure 3 This is a network structure diagram of the limit state prediction model of the present invention;

[0051] Figure 4 This is the network structure diagram of the fusion output module of the present invention;

[0052] Figure 5 2 is a comparison chart of the performance experimental results of the model of the present invention. DETAILED DESCRIPTION

[0053] The overall process of the present invention is as follows Figure 1 As shown:

[0054] Construction of the prediction data set for the ultimate state of polyurethane directly buried thermal insulation pipes: In the laboratory scenario, a simulated use environment for prefabricated polyurethane directly buried pipes is constructed. Ultimate stress tests are carried out on polyurethane directly buried pipes of different specifications under different environmental conditions to collect the maximum stress that the polyurethane directly buried pipes can withstand under such conditions. Record the corresponding specification data of the polyurethane directly buried pipes and the environmental data during the test, and use this as the input data of the model. Record the maximum stress it can withstand and use this as the output data of the model. To ensure the comprehensiveness of data collection, multiple batches of ultimate stress tests are carried out under multiple specifications and multiple environmental conditions to complete the construction of the prediction data set for the ultimate state of polyurethane directly buried thermal insulation pipes;

[0055] Construction of the feature screening module based on the relieff algorithm: In the data set construction stage, in order not to miss features with a high degree of correlation with the ultimate stress of polyurethane directly buried thermal insulation pipes, some redundant features will inevitably be introduced. These redundant features will increase the complexity of the model and also increase the difficulty of model training, resulting in the model being difficult to fit successfully. To eliminate the adverse effects of redundant features, the present invention designs a feature screening model based on the relieff algorithm to eliminate features with a relatively low degree of correlation with the ultimate stress of polyurethane directly buried thermal insulation pipes;

[0056] Construction of the prediction model for the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes based on the deep operator network architecture: Since the deep operator network architecture can approximate non-linear continuous operators with high accuracy and can maintain a small prediction error in the field of ultimate state prediction. Therefore, the present invention constructs a prediction model for the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes based on the deep operator network architecture. To further improve the performance of the model, a self-attention mechanism is designed inside the model;

[0057] Construction of the model architecture optimization module based on reinforcement learning: Based on the reinforcement learning strategy, a parameter optimization module for the model parameters of the prediction model for the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes is designed. The network architecture, state, action, and reward factor used in reinforcement learning are designed in detail to obtain the optimized prediction model for the ultimate state of polyurethane directly buried thermal insulation pipes;

[0058] Model deployment and use: Deploy the hardware devices required by the present invention and the corresponding prediction model for the ultimate state of polyurethane directly buried thermal insulation pipes to the relevant use environment. Input the specification parameter conditions of the thermal insulation pipe and the relevant environmental specification parameter conditions on the computer terminal, so as to obtain the ultimate stress situation of this prefabricated polyurethane directly buried thermal insulation pipe under such environmental conditions, achieving the effect of improving the engineering quality.

[0059] The method of the present invention will be further described below in conjunction with specific embodiments.

[0060] I. Data set construction

[0061] The present invention aims to predict the ultimate stress of directly buried polyurethane insulating pipes under real usage scenarios. Accurate training data is a prerequisite for ensuring the prediction accuracy. To this end, the present invention constructs an experimental scenario in the laboratory to simulate the actual usage of directly buried polyurethane insulating pipes, and conducts multiple batches of ultimate stress tests on directly buried polyurethane insulating pipes of different specifications under different experimental conditions to complete the construction of the ultimate state prediction data set of directly buried polyurethane insulating pipes.

[0062] Selection of directly buried polyurethane insulating pipe specifications: To ensure the comprehensiveness of data collection and not miss the features highly correlated with the ultimate stress of polyurethane insulating pipes. The present invention selects a variety of directly buried polyurethane insulating pipes of different specifications for experiments. The specification parameters selected in this embodiment include: the steel pipe layer thickness GT of the insulating pipe, the insulation layer thickness WT of the insulating pipe, the diameter BD of the insulating pipe, the width and length BW, BC of the insulating pipe, the average density BM of the insulating pipe, the compressive strength BQ of the insulating pipe, the thermal conductivity BR of the insulating pipe, and the water absorption rate BH of the insulating pipe;

[0063] Selection of experimental environment parameters: The ultimate stress of directly buried polyurethane insulating pipes is not only related to the self-specification parameters of the insulating pipes, but also greatly affected by their usage environment. Therefore, the experimental environment parameters selected in this embodiment are environmental temperature TM, buried depth SD of the directly buried pipe, surrounding soil density TD, soil friction coefficient TS, soil moisture TC, soil pH TJ, soil porosity TK, and soil bulk density TR.

[0064] Collection of ultimate stress data: Based on the described directly buried polyurethane insulating pipe specification parameters and experimental environment parameters. Conduct the ultimate stress test on directly buried polyurethane insulating pipes under a specific set of insulating pipe specification parameters and experimental environment parameters, and record the ultimate stress value JY. Record the directly buried polyurethane insulating pipe specification parameters GE and environmental parameters HC used in this set of tests. According to the above data records, obtain a complete set of data Datayz;

[0065] GE = [GT, WT, BD, BW, BC, BM, BQ, BR, BH]

[0066] HC = [TM, SD, TD, TS, TC, TJ, TK, TR]

[0067] Datayz = [GE, HC, JY]

[0068] Construction of the initial dataset for predicting the ultimate state of polyurethane directly buried insulation pipes: To ensure the comprehensiveness of data collection, a total of N groups of tests on the ultimate state of insulation pipes were carried out to cover all the specification parameters of polyurethane directly buried insulation pipes and the experimental environment parameters. The data collected from each test was combined to obtain the initial dataset Datay for predicting the ultimate state of polyurethane directly buried insulation pipes;

[0069] Datay = [Datayz1, Datayz2, Datayz3,..., Datayz N

[0070] Data standardization processing: To eliminate the dimensional difference and scale difference existing in the initial data, the data in the dataset Datay is unfolded and standardized. Take a piece of data Datayz i from the initial prediction dataset Datay of the ultimate state of polyurethane directly buried insulation pipes, where i ∈ [1, N]. Calculate the average value μ and standard deviation α of all specific parameter values α contained in this piece of data.

[0071]

[0072] Among them, α* represents the parameter value after the standardization operation, and α represents the initial parameter value in the data.

[0073] Construction of the standard prediction dataset for the ultimate state of polyurethane directly buried insulation pipes: Based on the described data standardization operation, each piece of data in the dataset Datay is standardized in turn. After the standardization operation is completed, the standard prediction dataset Data for the ultimate state of polyurethane directly buried insulation pipes is obtained.

[0074] II. Feature selection based on the relieff algorithm

[0075] In the dataset construction stage, to ensure the comprehensiveness of feature collection, some redundant features will inevitably be introduced into the dataset. These redundant features will increase the difficulty of model training and at the same time reduce the generalization ability of the model. To eliminate the influence of redundant features and enhance the model performance, the present invention constructs a feature selection module based on the relieff algorithm, which only retains the features with a high degree of correlation with the ultimate stress of polyurethane directly buried insulation pipes, so as to achieve the purpose of enhancing the generalization ability of the model. The process is as Figure 2 shown.

[0076] S1. Data selection and classification: Based on the obtained initial dataset Data of the ultimate state of polyurethane directly buried insulation pipes, randomly select M sample data from it. To facilitate the subsequent calculation of feature importance, after the data selection is completed, the M sample data are grouped according to the data similarity SD U.

[0077] ​

[0078] Among them, Datayz j represents the j-th parameter value of a piece of data in the data set Data, represents the j-th parameter value of another piece of data in the data set Data, where j ∈ [1, 18]. In this embodiment, the number of groups is 8, and the specific grouping process is as follows. Start grouping by randomly selecting a piece of data from M pieces of sample data, and sequentially calculate the SDU values of all the remaining data and this piece of data. After the calculation is completed, the piece of data with the largest SDU value and this piece of data are combined to form a data group. After the grouped data is removed, the remaining data continues to be classified as described above until the data grouping is completed;

[0079] S2, Feature difference degree calculation: Based on the described parameter selection, a total of 9 different specification characteristic parameters of polyurethane directly buried insulation pipes and 8 different environmental characteristic parameters are selected in this embodiment. Based on the sample grouping situation after the S1 process operation, calculate the feature difference degree TZC corresponding to each feature in the same group of samples and the feature difference degree BZC in different groups of samples respectively;

[0080]

[0081] Among them, the parameter values of the k-th feature of the w-th and m-th data in the same group of sample data, and respectively represent the parameter values of the k-th feature of the g-th and l-th data in different groups of sample data, |*| represents the absolute value operation function;[[ID=X]] [[ID=Y]]

[0082] S3, Feature correlation degree calculation: Based on the difference degrees TZC and BZC of each feature calculated in the S2 process, calculate the feature correlation degree Tmax:

[0083]

[0084] In the same group of samples, the smaller the feature difference degree, the stronger the correlation degree of this feature. While in different groups of samples, the larger the feature difference degree, the weaker the correlation degree of this feature;

[0085] S4, Multiple iteration feature correlation degree calculation: To avoid the calculation error caused by a single data sampling, in the present invention, the processes S1 to S3 are repeated, and the number of iterations is 6. The feature correlation degree calculated each time Perform average weighting processing for \(t\in[1,6]\), and define the result \(TO\) as the final feature importance. The higher the feature importance, the stronger the relationship between this feature and the ultimate stress of the polyurethane directly buried insulation pipe;

[0086]

[0087] Among them, represents the feature correlation degree calculated for the \(k\)-th feature in the \(t\)-th calculation;

[0088] S5. Feature screening process: Based on the calculated feature importance of each feature, it is found that the feature importance of the features in the last 40% among the 17 features has a significant decrease compared to the features in the first 60%. Therefore, in the present invention, the redundant features with feature importance in the last 40% will be discarded. After the feature screening is completed, the retained specification feature parameters of the polyurethane directly buried insulation pipe are: the steel pipe layer thickness \(GT\) of the insulation pipe, the insulation layer thickness \(WT\) of the insulation pipe, the diameter \(BD\) of the insulation pipe, the compressive strength \(BQ\) of the insulation pipe, the thermal conductivity \(BR\) of the insulation pipe, and the water absorption rate \(BH\) of the insulation pipe. The retained environmental feature parameters are the environmental temperature \(TM\), the buried depth \(SD\) of the directly buried pipe, the surrounding soil density \(TD\), and the soil friction coefficient \(TS\).

[0089] III. Construction of the ultimate state prediction model for prefabricated polyurethane directly buried insulation pipes based on the deep operator network architecture

[0090] The ultimate state prediction model for prefabricated polyurethane directly buried insulation pipes constructed in the present invention mainly consists of three parts, namely the backbone network structure, the branch network architecture, and the fusion output module, as Figure 3 shown.

[0091] Backbone network architecture design: The backbone network architecture consists of an input layer, an intermediate layer, and an output layer. In the present invention, the feature parameters of the polyurethane directly buried insulation pipe are selected as the input of this backbone network architecture. After feature screening, the number of retained feature parameters of the polyurethane directly buried insulation pipe is 6. Therefore, the number of neurons in the input layer is 6. The network architecture of the intermediate layer is mainly composed of a fully connected layer, a Relu activation function, and a Silu activation function. In the intermediate layer, first there are \(Wq\) fully connected layers, each fully connected layer has \(Ws\) neurons, and then it is activated by the Relu function to enhance the nonlinear expression ability of the model. Finally, there are \(Wh\) fully connected layers, each fully connected layer has \(Wf\) neurons in each layer and the Silu activation function. The designed output layer contains 4 different neurons, and then the outputs of the 4 neurons in the hidden layer are adaptively fused through the spatial attention layer to obtain the output \(Feaz\) of the backbone network architecture.

[0092] Branch Network Architecture Design: In the designed branch network architecture, a stacking structure is adopted to enhance the overall fitting ability of the model. The stacking structure consists of two parallel network paths, and each network path still consists of an input layer, an intermediate layer, and an output layer. In the branch network architecture, the environmental feature parameters are used as the input of this network architecture. Since the number of environmental feature parameters after screening is 4, a total of 4 neurons are included in the input layer. The structure design of the intermediate layer is as follows: First, there are 3 sequentially connected fully connected layers, and each fully connected layer contains 300 neurons. Subsequently, the data is non-linearly activated through the sigmoid activation function. Then, there are 4 sequentially connected fully connected layers, and in each of these 4 fully connected layers, there are 400 neurons. Finally, non-linear normalization processing is performed through the Softmax activation function. The number of neurons in the output layer of this network path is set to 1. The output FZS of the output layer in each network path is sent into the self-attention mechanism for further processing. Subsequently, a fully connected layer is used to further integrate the outputs of the two network paths to obtain the output Feaf of the branch network architecture.

[0093] Self-Attention Mechanism Design: This self-attention mechanism is used to map the output FZS of each network path in the branch network architecture into multiple feature spaces to further enhance the overall fitting ability of the model. The feature processing process of each self-attention head is as follows:

[0094]

[0095] Among them, Softmax(*) represents the Softmax linear normalization function, and Fea is the output result of the self-attention head. The matrices Q, K, and V are the query, key, and value matrices mapped after processing the feature FZS through three different matrices. Specifically:

[0096] Q = W Q *FZS

[0097] K = W K *FZS

[0098] V = W V *FZS

[0099] Among them, W Q ,W K ,W V are three different feature mapping matrices to send the feature FZS into three different feature spaces for processing and enhance the deep data fitting ability of the model.

[0100] In the self-attention mechanism designed in the present invention, a total of 8 parallel self-attention taps as described above are designed, and then the output of the 8 attention taps is integrated by the channel attention layer to obtain the output result Fe of the final self-attention mechanism:

[0101]

[0102] Among them, zyl(*) represents the operation function of the channel attention layer.

[0103] Fusion output module: The fusion output module further fuses the output Feaz obtained from the backbone network architecture and the output Feaf obtained from the branch network architecture to obtain the predicted result Yout of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe. First, the BN normalization layer is used to normalize Feaz and Feaf respectively to optimize the distribution of the feature channels to enhance the robustness of the model when outputting the predicted result of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe, and obtain and Secondly, the bilinear interpolation upsampling layer is used to upsample to the same size; again, the dot product layer is used to perform a dot product operation on the upsampled and to obtain the fused feature Ffz; again, the fused feature is input into the Dropout layer to optimize the feature channel parameters in Ffz to obtain the feature Ffz* with frozen parameters, reducing the risk of model overfitting; finally, Ffz* is input into the L2 regularization layer and then into the Relu activation function to obtain the final predicted result Your of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe, where the L2 regularization layer constrains the model weights and enhances the robustness of the model during prediction. The structure is as Figure 4 shown.

[0104] Yout = zrh(Feaf, Feaz)

[0105] Among them, zrh(*, *) represents the fusion output module.

[0106] IV. Model Optimization Training Based on Reinforcement Learning

[0107] In the present invention, in order to further improve the overall performance of the prefabricated polyurethane directly buried insulation pipe ultimate state prediction model and enhance the deep fitting ability of the model. After the overall architecture of the model is constructed, a model optimization module is constructed based on reinforcement learning technology to comprehensively optimize the number of fully connected layers in the intermediate layer of the model backbone network and the number of neurons Ws, Wf contained in each fully connected layer Wq, Wh, so as to find the specific number of fully connected layers and the number of neurons in each layer corresponding to the best model performance;

[0108] Loss function construction: It is clear that the loss function used by the model during the training phase is reinforcement learning to help construct the reward factor used in reinforcement learning. The present invention constructs the loss function Loss used in the model based on the mean square error and the absolute value error:

[0109]

[0110] where Yout o , o ∈ [1, N] represents the ultimate stress test result of the polyurethane directly buried thermal insulation pipe predicted by the model, and JY o represents the ultimate stress test result of the polyurethane directly buried thermal insulation pipe obtained from the actual test.

[0111] State, action, and reward factor construction: The reinforcement learning strategy used in the present invention is to optimize the model parameters in the ultimate state prediction model of the polyurethane directly buried thermal insulation pipe. Therefore, a vector composed of the number of fully connected layers Wq, Wh and the corresponding number of neurons Ws, Wf, a total of 4 numerical combinations, constitutes the state space corresponding to this reinforcement learning. In its action space, it is the adaptive optimization and change of the parameters Wq, Wh, Ws, Wf to find the model parameters corresponding to the best model performance. The design of its reward factor is achieved by comparing the magnitudes of the loss functions. Using the data set constructed in process S1, based on the model parameters corresponding in the state space, the loss function value Lossq of the model at this time is obtained. When the state is updated, that is, an action to update the model parameters is taken. Then calculate the loss function value Lossh corresponding to the model after the parameter update. If the value of Qjl is greater than 0, a positive reward value of 2 * Qjl is given to this action. If the value of Qjl is less than 0, a negative reward value of 2 * Qjl is given to this action;

[0112] Qjl = Lossh - Lossq

[0113] Reinforcement learning network construction: In the present invention, a policy network is constructed to optimize the parameters Wq, Wh, Ws, Wf to find the network model parameters that make the model performance the best. A value network is constructed to adaptively judge the value corresponding to taking this set of model parameters. In the present invention, the optimized parameters are relatively simple. Therefore, the designed policy network mainly consists of 4 interconnected fully connected layers and the Relu activation function. The designed value network mainly consists of 6 interconnected fully connected layers and the Silu activation function. The input of its network architecture is a four-dimensional vector, representing a corresponding situation of Wq, Wh, Ws, Wf. Using the data set constructed in process S1 and placing the prediction model constructed by the input of this architecture at this time for training, and recording the loss function after the training process to calculate the reward factor.

[0114] Reinforcement learning network training: Randomly initialize a set of network parameters \(W_q\), \(W_h\), \(W_s\), \(W_f\), and send them into the reinforcement learning network for optimization training. Adopt a stochastic gradient descent training strategy and combine the reward factor to adaptively update the parameters \(W_q\), \(W_h\), \(W_s\), \(W_f\) using the policy network and the value network. Continuously update the model parameters. When the loss function of the model is less than the set threshold \(T_h\), stop the training of the model. And set the obtained network parameters \(W_q\), \(W_h\), \(W_s\), \(W_f\) as the number of fully connected layers and the number of neurons used in the final fully connected layer in the middle layer of the model. Obtain the final optimized limit state prediction model of the polyurethane directly buried thermal insulation pipe.

[0115] V. Model Deployment and Application

[0116] Model deployment: Deploy the final optimized limit state prediction model of the polyurethane directly buried thermal insulation pipe and the parameter standardization program to the computer terminal to ensure that this prediction model can function smoothly and guarantee the smooth realization of its prediction function.

[0117] Development of parameter input interface: Develop a visual interface for inputting the specification parameters and usage environment parameters of the polyurethane directly buried thermal insulation pipe, and deploy this interface and the software in the backend to the computer terminal, and adapt it to the prediction model adaptively to ensure that its parameters can be smoothly passed into the limit state prediction model of the polyurethane directly buried thermal insulation pipe.

[0118] Ultimate stress prediction: Send various specification parameters input from the parameter input interface into the limit state prediction model of the polyurethane directly buried thermal insulation pipe after being processed by the data standardization program to obtain the ultimate stress result of this specification of polyurethane directly buried thermal insulation pipe under this environmental condition.

[0119] Data analysis: According to the analysis results of the predicted ultimate stress, help relevant staff understand the ultimate stress conditions of various specifications of polyurethane prefabricated directly buried thermal insulation pipes under different environmental conditions, so as to reasonably select the polyurethane prefabricated directly buried thermal insulation pipes used in the project according to the predicted ultimate stress conditions and the actual usage scenarios, so as to further improve the overall quality of the project.

[0120] VI. Experimental Results and Explanation

[0121] To verify that the method proposed in the present invention has a high detection accuracy, a performance comparison experiment is carried out with the existing limit state prediction algorithm for polyurethane pre-buried thermal insulation pipes. The results are as Figure 5As shown. The two existing algorithms are the TFRN algorithm and the BRTY algorithm respectively. Except for the environmental temperature parameter, the insulation pipe specification parameters and the remaining environmental parameters are kept the same, and comparative experiments are carried out at different operating temperatures (-20 °C, 0 °C, 10 °C, 30 °C). The quality of the algorithm performance is realized by comparing its prediction error value (MAPE value). The larger the MAPE value, the greater the prediction error of this algorithm, that is, the overall performance of the model is poor.

[0122]

[0123] Among them, yi represents the true ultimate stress value of the polyurethane pre-buried insulation pipe in the experiment, represents the predicted ultimate stress value of the polyurethane pre-buried insulation pipe. The results show that at the four temperature points, the algorithm proposed in the present invention has a smaller prediction error, that is, it shows that the algorithm proposed in the present invention has better performance than the existing algorithms.

[0124] The above are only the preferred embodiments of the present application and are not used to limit the present application. For those skilled in the art, the present application can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

[0125] Although the specific implementation manners of the present invention are described above, it is not a limitation on the protection scope of the present invention. Those skilled in the art should understand that based on the technical solution of the present invention, various modifications or deformations that can be made without creative efforts by those skilled in the art are still within the protection scope of the present invention.

Claims

1. A method for constructing a prediction model for the ultimate state of prefabricated polyurethane directly buried thermal insulation pipes, characterized in that, It includes the following processes: Step 1, based on the simulation of the use environment, carry out ultimate stress tests on polyurethane directly buried insulation pipes of different specifications to collect the maximum stress that can be withstood under such conditions, record the maximum stress, the corresponding specification data and environmental data, and construct an initial data set; Step 2, the feature screening module based on the relieff algorithm eliminates features with low relevance to the ultimate stress of polyurethane directly buried insulation pipes, obtains the specification data and environmental data with the top feature importance, and constructs a training data set in combination with the corresponding maximum stress; Step 3, construct an ultimate state prediction model based on the deep operator network architecture, including a backbone network architecture, a branch network architecture, and a fusion output module. The backbone network is used to extract the physical property features of polyurethane directly buried insulation pipes, and the branch network is used to extract the environmental features of polyurethane directly buried insulation pipes and perform weighted processing on the feature channels in the environmental features based on the self-attention mechanism to obtain self-attention environmental features. The fusion output module is used to adaptively fuse the two types of features output by the two network architectures and output the maximum stress prediction result based on the fused features; Step 4, optimize and train the constructed ultimate state prediction model based on the reinforcement learning strategy to obtain the final ultimate state prediction model.

2. The method for constructing a prediction model for the ultimate state of a prefabricated polyurethane directly buried thermal insulation pipe according to claim 1, wherein: The specific data processing process of the feature screening module based on the relieff algorithm in Step 2 is as follows: S21, data selection and classification: Based on the initial data set Data of the ultimate state of polyurethane directly buried insulation pipes, randomly select M sample data from it, and group the M sample data according to the data similarity SDU. The number of groups is L; S22, feature difference degree calculation: Calculate the feature difference degree TZC corresponding to each feature in the same group of samples and the feature difference degree BZC in different groups of samples respectively; S23, feature correlation calculation: Calculate the feature correlation Tmax based on the calculated difference degrees TZC and BZC of each feature; S24, Calculate the feature correlation multiple times: Repeat the process from S21 to S23 for a total of B iterations. Average and weight the feature correlation obtained in each calculation, where t ∈ [1, B], and define the result TO as the final feature importance. The higher the feature importance, the stronger the connection between this feature and the ultimate stress of polyurethane directly buried insulation pipes; S25, feature screening process: Based on the calculated feature importance of each feature, discard the redundant features with the feature importance in the last 40%.

3. A method for constructing a prediction model of the ultimate state of a prefabricated polyurethane directly buried insulation pipe according to claim 2, characterized in that: In S21, Among them, Datayz j represents the j-th parameter value of a piece of data in the data set Data, represents the j-th parameter value of another piece of data in the data set Data, where j ∈ [1, C]; C is the total number of environmental data, specification data plus maximum stress data; The specific process of grouping is as follows: randomly select one piece of data from M pieces of sample data to start grouping, and sequentially calculate the SDU values of all the other remaining data and this piece of data. After the calculation, combine the piece of data with the largest SDU value and this piece of data to form a data group; eliminate the grouped data, and continue to perform the above-mentioned data classification on the remaining data until the data grouping is completed; In S22, calculate the feature difference degree TZC corresponding to each feature in the same group of samples and the feature difference degree BZC in different groups of samples respectively. The formula is as follows: Among them, the parameter value of the k-th feature of the w-th and m-th data in the same group of sample data, k ∈ [1, T], where T is the sum of the environmental data and the specification data; and respectively represent the parameter values of the k-th feature of the g-th and l-th data in different groups of sample data, |*| represents the absolute value operation function.

4. The method for constructing a prediction model for the ultimate state of a prefabricated polyurethane directly buried insulation pipe according to claim 1, characterized in that: The backbone network architecture consists of three parts: an input layer, a middle layer, and an output layer; the specification characteristic parameters of the polyurethane directly buried thermal insulation pipe are selected as the input of this backbone network architecture. After feature screening, the number of retained specification characteristic parameters of the polyurethane directly buried thermal insulation pipe is p. Therefore, the number of neurons in the input layer is p; the network architecture of the middle layer is mainly composed of three parts: a fully connected layer, a Relu activation function, and a Silu activation function. In the middle layer, there are first Wq fully connected layers, with Ws neurons in each fully connected layer. Subsequently, it is activated by the Relu function to enhance the non-linear expression ability of the model. Finally, there are Wh fully connected layers, with Wf neurons in each layer of each fully connected layer and a Silu activation function; the designed output layer contains 4 different neurons. Subsequently, through the spatial attention layer, the outputs of the 4 neurons in the hidden layer are adaptively fused to obtain the output Feaz of the backbone network architecture.

5. A method for constructing a prediction model for the ultimate state of a prefabricated polyurethane directly buried thermal insulation pipe as described in claim 1, characterized in that: A stacked structure is adopted in the branch network architecture; the stacked structure is composed of 2 parallel network paths, and each network path consists of three parts: an input layer, a middle layer, and an output layer; the environmental characteristic parameters are used as the input of this branch network architecture. A total of q neurons are included in the input layer, where q is the total number of retained environmental characteristic data after screening; the structure of the middle layer is as follows: first, there are 3 sequentially connected fully connected layers, and then the data is non-linearly activated by the simgoid activation function. Subsequently, there are 4 sequentially connected fully connected layers. In these 4 fully connected layers, finally, non-linear normalization processing is performed through the Softmax activation function; the number of neurons in the output layer of the network path is set to 1; the output FZS of the output layer in each network path is sent to the self-attention mechanism unit for processing to obtain the output result Fe. Subsequently, 1 fully connected layer is used to further integrate the outputs of the 2 network paths to obtain the output Feaf of the branch network architecture.

6. The construction method of the limit state prediction model of a prefabricated polyurethane directly buried thermal insulation pipe as described in claim 5, characterized in that: The self-attention mechanism unit maps the output FZS of each network path in the branch network architecture into multiple feature spaces. The self-attention mechanism includes L parallel self-attention taps, and then the channel attention layer is used to integrate the outputs of the L attention taps to obtain the output result Fe of the final self-attention mechanism; the feature processing process of each self-attention tap is as follows: Among them, Softmax(*) represents the Softmax linear normalization function, and Fea is the output result of the self-attention tap; the matrices Q, K, and V are the query, key, and value matrices mapped after the feature FZS is processed by three different matrices. Specifically: Q = W Q *FZS K = W K *FZS V = W V *FZS Among them, W Q , W K , W V are three different feature mapping matrices to send the feature FZS into three different feature spaces for processing.

7. The method for constructing a prediction model for the ultimate state of a prefabricated polyurethane directly buried thermal insulation pipe according to claim 1, characterized in that: The fusion output module further fuses the output Feaz obtained from the backbone network architecture and the output Feaf obtained from the branch network architecture to obtain the predicted result Yout of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe; First, the BN normalization layer is used to normalize Feaz and Feaf respectively to optimize the distribution of feature channels, and obtain and Secondly, the bilinear interpolation upsampling layer is used to upsample to the same size; Thirdly, the dot product layer is used to perform a dot product operation on the upsampled and to obtain the fused feature Ffz; Subsequently, the fused feature is input into the Dropout layer to optimize the feature channel parameters in Ffz to obtain the feature Ffz * with frozen parameters; Finally, Ffz * is input into the L2 regularization layer and then into the Relu activation function to obtain the final predicted result Yout of the ultimate stress of the prefabricated polyurethane directly buried insulation pipe, where the L2 regularization layer constrains the model weights.

8. The construction method of the limit state prediction model of a prefabricated polyurethane directly buried thermal insulation pipe as described in claim 1, wherein: The reinforcement learning strategy comprehensively optimizes the number of fully connected layers in the middle layer of the model backbone network and the number of neurons Ws and Wf included in each fully connected layer Wq and Wh to find the specific number of fully connected layers and the number of neurons in each layer corresponding to the best model performance. The vector composed of the number of fully connected layers Wq, Wh and the corresponding number of neurons Ws, Wf is used as the state space corresponding to this reinforcement learning; in its action space, it is the adaptive optimization and change of the parameters Wq, Wh, Ws, Wf to find the model parameters corresponding to the best model performance; the design of its reward factor is realized by comparing the magnitudes of loss functions. Using the constructed training dataset, based on the model parameters corresponding to the state space, the loss function value Lossq of the model at this time is obtained. When the state is updated, that is, an action to update the model parameters is taken, and then calculate the loss function value Lossh corresponding to the model after the parameter update; if the value of Qil is greater than 0, give this action a positive reward value of 2*Qil, if the value of Qil is less than 0, give this action a negative reward value of 2*Qil: Qil = Lossh - Lossq At the same time, construct a reinforcement learning network, construct a policy network to optimize the parameters Wq, Wh, Ws, Wf, and construct a value network to adaptively judge the value corresponding to this set of model parameters; the designed policy network is mainly composed of 4 interconnected fully connected layers and a Relu activation function, and the designed value network is mainly composed of 6 interconnected fully connected layers and a Silu activation function; the input of its network architecture is a four-dimensional vector, representing a corresponding situation of Wq, Wh, Ws, Wf. Use the constructed dataset and place the prediction model constructed by the input of this architecture at this time for training, and record the loss function after the training process to calculate the reward factor.

9. A method for predicting the ultimate state of a prefabricated polyurethane directly buried thermal insulation pipe, characterized in that: Deploy the ultimate state prediction model constructed by the construction method described in any one of claims 1 to 8 to a computer terminal, and include the following process: Obtain the specification data and environmental data of the polyurethane directly buried insulation pipe in real time, and standardize them into model input data; Input the standardized specification data and environmental data into the ultimate state prediction model; output the ultimate stress result of the polyurethane directly buried insulation pipe of this specification under this environmental condition.

10. A method for predicting the ultimate state of a prefabricated polyurethane directly buried insulation pipe as described in claim 9, characterized in that: The polyurethane directly buried insulation pipe specification characteristic parameters retained after screening for prediction are: the steel pipe layer thickness GT of the insulation pipe, the insulation layer thickness WT of the insulation pipe, the diameter BD of the insulation pipe, the compressive strength BQ of the insulation pipe, the thermal conductivity BR of the insulation pipe, and the water absorption rate BH of the insulation pipe; the retained environmental characteristic parameters are the environmental temperature TM, the buried depth SD of the polyurethane directly buried insulation pipe, the surrounding soil density TD, and the soil friction coefficient TS.

Images