JS Journal of Sensors 1687-7268 1687-725X Hindawi Publishing Corporation 10.1155/2016/8126214 8126214 Research Article Modeling and Calculation of Dent Based on Pipeline Bending Strain http://orcid.org/0000-0001-9423-9706 Feng Qingshan 1,2 Li Rui 2,3 http://orcid.org/0000-0003-1010-4260 Zhang Hong 1 Tian Guiyun 1 China University of Petroleum Beijing 102249 China cup.edu.cn 2 Petrochina Pipeline Company Langfang 065000 China petrochina.com.cn 3 School of Automation Science and Electrical Engineering Beihang University Beijing 100191 China buaa.edu.cn 2016 16 8 2016 2016 27 01 2016 10 03 2016 03 04 2016 2016 Copyright © 2016 Qingshan Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The bending strain of long-distance oil and gas pipelines can be calculated by the in-line inspection tool which used inertial measurement unit (IMU). The bending strain is used to evaluate the strain and displacement of the pipeline. During the bending strain inspection, the dent existing in the pipeline can affect the bending strain data as well. This paper presents a novel method to model and calculate the pipeline dent based on the bending strain. The technique takes inertial mapping data from in-line inspection and calculates depth of dent in the pipeline using Bayesian statistical theory and neural network. To verify accuracy of the proposed method, an in-line inspection tool is used to inspect pipeline to gather data. The calculation of dent shows the method is accurate for the dent, and the mean relative error is 2.44%. The new method provides not only strain of the pipeline dent but also the depth of dent. It is more benefit for integrity management of pipeline for the safety of the pipeline.

Petrochina Pipeline Company
1. Introduction

With the development of the oil and gas production, the long-distance buried pipeline is used to the transport the production of oil and gas [1, 2]. Due to the reasons of time accumulated or construction, defects such as corrosion, gouge, dent, and displacement seriously threat the safety operation of the pipeline. In order to reduce the risk of these defects, the pipeline company usually selects in-line inspection tool to inspect the defects for body of pipeline. The magnetic flux leakage (MFL) and ultrasonic tool are used to inspect metal loss of pipeline such as corrosion, crack, and gouge.

Dent is another common defect in long-distance buried oil and gas pipeline. A lot of dents have been found in the pipeline due to the construction of mechanical damage. Dents can deform the pipeline and affect the integrity of body for the pipeline . Specifically, the dent located in the weld of pipeline or the dent that is very sharp can lead to the serious leak consequence. The information about the severity of the damage to a pipe due to dents is very important for the pipeline industry. The defect of dent is the radial deformations of the pipe wall. Generally, these dents can be measured by in-line inspection tools named geometry tool which equipped with instrumentation and mechanical fingers is able to provide data on the pipe wall geometry deformation. However, this tool only inspects the size (length, depth, and width) of the dent. In order to assess the hazard for the dent, the technical requires a complex modeling and computes the strain or other information. Another method for dent testing can be using the profile gauge to measure the depth, width, and length of the gouge or scratch damage. This method also can be used to determine the extent of damage and determine if any sharp edges are present. Recently, 3D laser scanner is used to map and measure the dent. This technology can more accurately measure the dimensions of the dent by the laser scan. However, the method for dent testing of profile gauge and 3D laser scan only use dents which have been dug. They are the verification for dents which are inspected by geometry tool.

Nowadays, another Pipeline In-Line Inspection tool which installed inertial measurement unit (IMU) is used to inspect the pipeline . As shown in Figure 1, this tool, which consists of a Data Acquisition System (DAS), an IMU, some weld detectors, and odometers, is driven to move through the pipeline and collect data regarding dents, bends, and navigation. Petrochina Pipeline Company has developed an IMU system based on a navigation system, which is suitable for a geometry PIG. The IMU can be used to record the attitude data for the tool during the inspection. A method of calculation for pipeline centerline is proposed by fusion for multisensors. Czyz et al. [11, 12] developed a pure geometry-based algorithm for the bending strain of a pipeline. Another pipeline strain testing method is using the strain gauge to be installed on the surface of pipeline. Two types of strain gauge, resistance and optical fiber, are used to measure the strain of pipeline. The strain gauge of optical fiber is more precise for measuring the strain. These two methods can be used to monitor the strain changes of the pipeline. But the limit for these two methods is only used to measure the part sections of pipeline. They cannot inspect the whole strain for long-distance buried pipeline.

Measurement system of in-line inspection of pipeline centerline.

The bending strain can be reported for the whole pipeline when the tool runs for once. Generally, the technical just uses the bending strain to detect the position of pipeline enduring the extra strain larger than the operation requirement and assessing the displacement for local pipeline. However, it is found that the bending strain can be affected by the dent . In this paper, a method for modeling and computation of dent which is based on bending strain is proposed. The technique takes Bayesian statistical theory and neural network to improve a new method to compute the depth of the dent. It is more useful to assess the dent for the pipeline integrity.

2. The Calculation Method of the Pipeline Bending Strain

In this section, the calculation method of the pipeline bending strain is described. Several factors such as internal pressure, temperature differential, external loads, and boundary conditions are affecting the strain. The strain can be consisted of two primary components: longitudinal and hoop strain . The bending strain is one of the components for longitudinal strain. It can be described by the following formula: (1) ε = D 2 k , ε v = D 2 k v , ε h = D 2 k h , where ε is the total bending strain, k is the total curvature, k v is vertical curvature, and k h is horizontal curvature.

The total curvature of the centerline of a pipe is described at each point along the pipeline by the curvature vector. In order to calculate the pipeline curvature, the centerline of a pipe is considered as a 3D parametric curve described in a Cartesian system by a vector v ( s ) , which is a function of a distance ( s ) along the curve [20, 21]: (2) v s = x s , y s , z s .

Assume that the vector t is tangent of v ( s ) , separating the vertical and horizontal curvature components as shown in Figure 2. The calculation of the pipeline bending strain can be given as (3) t x = cos P sin A , t y = cos P cos A , t z = sin P , where the pitch ( P ) and azimuth ( A ) of the pipeline centerline can be measured by the PIG.

Pitch ( P ) and azimuth ( A ) for the pipeline centerline.

Assume that the vector k is the curvature vector of a 3D curve at a given point, and k consists of the vertical curvature k v and the horizontal curvature k h , which can be given as (4) k s = d t d s , k = k v 2 + k h 2 .

The above equation can be written separately for each component of the curvature vector in the Cartesian system: (5) k x = d t x d s , k y = d t y d s , k z = d t z d s .

Based on (1)–(4), the components of the defined curvature vector can be calculated as follows: (6) k x = - sin P d P d s sin A + cos P cos A d A d s , k y = - sin P d P d s cos A - cos P sin A d A d s , k z = cos P d P d s .

The vertical curvature k v and the horizontal curvature k h can be given as (7) k v = - d P d s , k h = - d A d s cos P .

From (7), it can be seen that the pipeline bending strain can be calculated with the attitude data, which can be acquired with a PIG.

3. Modeling and Computation of Dent

In Section 2, a computation of the pipeline bending strain is discussed. The total of the features such as displacement and dent based on bending strain can be computed and identified. However, it is difficult to obtain the depth of dent by the bending strain measurement. Because of the lots of factors, it is hard to find a simple relationship between bending strain and size of dent. Neural network is a good solution for this practical project. It is based on complex connection of large number of neurons and deal with difficult language of the modeling information by self-learning, self-organized, and nonlinear dynamics.

3.1. BP Neural Network

The BP neural network consists of an input layer, a hidden layer, and a linear output layer. Each layer is composed of several neurons . The most basic three-layer BP neural network structure is shown in Figure 3. It is assumed that each layer consisted of “ N ” processing elements; the training set included “ M ” samples mode for ( X k , Y k ). For the “ p ” training sample, the sum of total input is N p j for j element. If the output is O p j , it follows (8) N P j = i - 0 N w j i O P j , O p j = f N P j , where w j i is the weight for neurons between i and j . f is the function as (9) f x = 1 1 + e - x .

A typical signal of dent feature from bending strain.

The error of neural network is as follows: (10) E = p E P , E P = 1 2 i d P j - O p j 2 , where d p j is expect output of the j output for input of p . If the learning rule used a method for gradient descent, the weights can be changed by error function. The modified function is as follows: (11) w i j t + 1 = w i j t + η δ p j O p j , δ p j = f N p j d p j - O p j , δ p j = f N p j k δ p k w k j , where t is the times of learning, η is the learning factor, and p is the value of modified error. δ p j is the modify error for output and δ p j is the modified error for hidden layer.

3.2. Improving the Neural Network Based on Bayesian Theory

The traditional neural network has difficulty in controlling its model complexity, leading to network overfitting, long time training, and network model low stability [25, 26]. However, the Bayesian neural network based on Bayesian reasoning can solve these problems effectively through correcting network training performance function.

In the actual complex system, the input quantity includes given input quantity and unknown disturbance quantity such as (noise signals), and normally the unknown disturbance quantity has an effect on the system output [27, 28]. Let x = [ x 1 , x 2 , , x n ] as an observable given input; then the relationship between system output quantity y and input quantity x can be described by the following formula: (12) y = f x + ε , where ε represents the effect of input on the output and the random quantity submitting to some distribution. The training performance function of neural network adopts mean square error function. Assume error function is as follows: (13) E 0 = 1 2 k - 1 K n - 1 N y n k - y n k 2 , where N is the sample numbers, K is the neural network output numbers, y n k is the expectation output, and y n k is the network actual output. Generally the not unique E 0 solution may lead the neural network training to fall into partial minimum value. To solve this problem, a constraint item can be introduced to make function f ( x ) have interpolation ability, and then the solution for E 0 is stable. The condition that f ( x ) has interpolation ability is that f ( x ) is smooth. While the square sum of network weight and threshold is small, its output is smoother. Therefore constraint item should represent smooth constraint, and the objective function is expressed as follows: (14) E ALL = α 0 1 2 i = 1 W w i 2 + β 0 1 2 k = 1 K n = 1 N y n k - y n k 2 = α E n e t w + β E n e t y .

W is the neural network weight and threshold numbers, w is network weight and threshold, E n e t ( w ) is the square sum of weight and threshold, E n e t ( y ) is the square sum of network actual output and objective expectation value residual, and α and β are hyperparameter which decide the value size of the neural network objective performance function and control the distribution form of weight and threshold. Through using new objective performance function, the network weight and threshold are as little as possible, on condition that the network training error is as small as possible. According to Bayesian, posterior distribution probability of α and β meets the following formula: (15) P α , β D , N h = P D α , β , N h P α , β N h P D N h .

D = { x ( N ) , y ( N ) } is a sample data set which consists of N samples, N h is neural network hidden layer number, P ( α , β N h ) is hyperparameter prior probability, P ( D N h ) is normalization factor, and P ( D α , β , N h ) is likelihood function. After partial derivatives to α , β , the most salient hyperparameter can be carried out: (16) α M P = γ 2 E n e t w M P , β M P = N - γ 2 E n e t w M P .

γ is the number that can reduce parameters numbers (0–W) in performance index function parameters in network. w M P is the weight and threshold when E A L L is minimum. M P is the subscript for w . w M P is the weight and threshold when E A L L is minimum. After working out the most salient hyperparameter α , β determine whether E A L L is convergence again. Through iterative judgment, the best and most salient neural network model can be obtained.

Relative to traditional neural network, the Bayesian neural network is focused on the probability distribution of the whole parameter space, and the predicted results are based on the statistical average of parameters posterior distribution. A single model is corresponding to one point in the space; all models are corresponding to the whole parameter space. Therefore the Bayesian neural network guarantees network stronger generalization in theory.

3.3. Modeling and Computation for Dent by Bending Strain

A typical signal of dent feature from bending strain is shown in Figure 3. The horizontal axis represents the dent length along the pipeline and the vertical axis is the bending strain for dent, respectively. The bending strain can be represented for the severity for the dent. To the identity of the dent feature, the input can be related to the bending strain information for the dent. The output selected the actual depth of the dent from the site measurement.

According to the relationship for the bending strain and dent as shown in Figure 3, pattern in a plot of bending strain against distance shows the characteristic of environmentally induced deformation for dent. The pattern is caused by a primary action in the center, with a reaction (bending strain in the opposite direction) on both sides. In the case of a dent the central deformation is an underbend with overbends at either end. The main influencing parameters are including the beginning bending strain S 1 and ending bending strain S 3 which is related to the reaction area of dent along the pipeline, respectively. The largest bending strain S 2 is related to the deepest of the dent. X 1 is the distance for reaction of the beginning bending strain for dent. X 2 is the distance for primary action for dent. X 3 is the distance for another reaction of the ending bending strain for dent. The process of modeling and computation for the dent based on Bayesian theory neural network is shown in Figure 4.

Dent modeling and computation based on Bayesian theory neural network.

4. Field Test and Data Analysis 4.1. Equipment and Performance

To test the proposed method experimentally, an in-line inspection tool (shown in Figure 5) with the proposed method is used to inspect oil pipeline, which is about 150 kilometers. About 150 dents have been dug to measure and verify depth gauge in past five years. These actual data can be selected to be used for sample for the proposed algorithm.

IMU Pipeline Inspection Tool.

The IMU pipeline ILI tool mainly consists of the following sensors to be used to calculate the bending strain and position of the pipeline:

Inertial measurement unit is the main device to gather the ILI data. Three gyroscopes and three accelerometers orthogonally mounted on the IMU to measure the attitude and positon for ILI tool. These gathered data can be used to calculate the pipeline bending strain with other sensors mounted on the ILI tool.

Two odometers mounted on the ILI to measure the distance and instantaneous and average speed, which can be used to modify the inertial system errors.

Two weld detector sensors are used to measure and record time of passing each girth weld for tool. They can use the alignment for every spool for repeat inspection.

As is known, most domestic oil pipelines are heating transportation and pass through mountains, hills, rivers, and other complex environmental areas ; the technical and safety requirements of electrical equipment are extremely strict. To safely inspect the actual pipeline, the IMU should meet not only the technical performance, but also the actual situation of the pipe and the external environment should be considered. Technical performances of inertial devices, which are used in the field test, are shown in Table 1.

Characteristics of sensors.

Sensor Characteristics Magnitude
Gyroscope Bias <0.01°/h
Random walk 0.002 ° / h

Accelerometer Bias stability <50  μ g
Scaling factor <50 ppm

Odometer Scaling factor <0.3%
White noise <0.1 m/s

Landmark White noise < ± 1 m
4.2. Data Analysis

According to Section 2, a typical bending strain for dent is shown in Figure 6, whose depth is 52.79 mm for 6.5% of outer diameter measured with gauge in field. The bending strain for this dent is −0.4125% which is located in bottom of the pipeline.

Typical bending strain of a dent.

For the method of Section 3, 120 data sets for bending strain of dent are selected for the input of the new neural network. Part of dents of training data are shown in Table 2. The input of neural network is consisted of S 1 , S 2 , S 3 , X 1 , X 2 , and X 3 introduced in Section 3.3. The output of neural network is the depth of the dent.

Part of dents of training data.

Depth ratio (%) Depth (mm) S 1 (%) S 2 (%) S 3 (%) X 1 (m) X 2 (m) X 3 (m)
2.02 16.45 0.067 0.067 0.079 5 18 9
1.76 14.29 0.111 0.111 0.023 9 16 2
1.39 11.28 0.016 0.016 0.046 1 23 6
1.49 12.13 0.128 0.128 0.149 12 17 13
1.99 16.15 0.068 0.068 0.101 6 15 9
1.88 15.3 0.055 0.055 0.080 3 14 5

The measured data input into the BP neural network to calculate the result of the dent. Parts of sample data input into the network for training and the other sample data to verify the accuracy. The calculation is shown in Table 3 and Figure 7. The mean relative error for BP calculation and actual depth is 12.26%. Although the BP neural network can be used to calculate the depth from the bending strain, the accuracy degree is low to assess the dent.

Calculation and error for BP neural network.

Actual depth (mm) Calculation of BP (mm) Error for calculation of BP and actual depth (mm) Relative error (%)
14.7153 13.38198 1.33332 9.06
15.8535 12.73971 3.11379 19.64
15.1218 12.85353 2.26827 15.00
15.1218 13.10556 2.01624 13.33
11.3007 13.32507 2.02437 17.91
9.6747 10.19502 0.52032 5.38
15.2844 14.43888 0.84552 5.53

Calculation and error for BP neural network.

To verify the accuracy of the proposed method, all of data input into the modified Bayesian neural network. The neural network can be trained by the sample data. Mean square error (MSE) is selected to assess the quality of the neural network. To adjust the parameter of neural network, the result of training is shown in Figure 8. The best validation performance of MSE for the method is 0.0040885 at the epoch 95.

Training for modified neural network based on Bayesian.

The calculation of modified neural network based on Bayesian is shown in Figure 9 and Table 4. It is clear to see the method can be used to calculate the depth of dent from bending strain. From Table 4, the mean relative error for calculation and actual depth is 2.44%. The accuracy of modified neural network is raised more than BP neural network. The calculation of modified Bayesian neural network can be used to calculate the depth of dent from dent. With this method, the pipeline company runs the IMU tool for once not only to obtain the pipeline bending strain but also to calculate all of the depths of dent. The proposed method offers a useful method for pipeline integrity evaluation.

Calculation and error for modified Bayesian neural network.

Actual depth (mm) Calculation of Bayesian (mm) Error for calculation of Bayesian and actual depth (mm) Relative error (%)
14.7153 14.98349 0.268193 1.82
15.8535 15.7455 0.108 0.68
15.1218 14.17717 0.94463 6.25
15.1218 14.96672 0.15508 1.03
11.3007 11.07863 0.22207 1.97
9.6747 10.10614 0.431435 4.46
15.2844 15.41688 0.132484 0.87

Calculation and error for modified Bayesian neural network.

5. Conclusions

The in-line inspection tool, which loads the IMU, used to inspect the centerline of the pipeline. The attitude information can be used to compute the bending strain of the pipeline. However, there is no paper or report to research the relationship between the bending strain and dent. In this paper, based on the analysis of the calculation method of pipeline bending strain, we propose a method based on modified Bayesian neural network to calculate dent depths from bending strain. To test the proposed method experimentally, a PIG with the proposed method is used to inspect a 150 km pipeline. It can be obtained that

A new method is proposed to verify the relationship between pipeline bending strain and dent based on the calculation of bending strain.

According to the characteristic of signal between dent and pipeline bending strain, a new model is presented in detail to be used for proposed calculation of algorithm.

The calculation of modified Bayesian neural network is proposed to compute the depth dent with pipeline bending strain to compare with actual data. And traditional BP neural network is used to calculate the depth of dent. According to the result of calculation, the mean accuracy of calculation for modified Bayesian neural network is better than BP neural network for 9.82%.

According to calculation, the proposed method is more accurate and suitable for the calculation of depth based on pipeline bending strain. The mean relative error for calculation of dent depth based the modified Bayesian neural network is 2.44%.

This paper provides a novel method for calculating dent depths from the pipeline bending strain. The dent can be evaluated by the bending strain not used to dig or use another tool to reinspect. The bending strain and calculation of dent depth is also useful to the evaluation of pipeline integrity.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the project of Petrochina Pipeline Company “the research of safety service of buried pipeline in permafrost region.”

Aihua Q. Chinese oil and gas pipeline transport development status and analysis of associated problems International Petroleum Economics 2009 17 12 17 12 Yang L.-J. Li R.-Q. Gao S.-W. Yang Y. Navigating and positioning technique for inner detection of pipeline Journal of Shenyang University of Technology 2012 34 4 427 432 2-s2.0-84865977628 Sha S. Feng Q. Excavated verification technology of pit defects of buried pipeline Oil and Gas Storage and Transportation 2014 8 834 838 Ghaednia H. Das S. Wang R. Kania R. Safe burst strength of a pipeline with dent-crack defect: effect of crack depth and operating pressure Engineering Failure Analysis 2015 55 288 299 10.1016/j.engfailanal.2015.06.005 2-s2.0-84934267422 Pinheiro B. Pasqualino I. Cunha S. Fatigue life assessment of damaged pipelines under cyclic internal pressure: pipelines with longitudinal and transverse plain dents International Journal of Fatigue 2014 68 38 47 10.1016/j.ijfatigue.2014.06.003 2-s2.0-84905369618 Allouti M. Schmitt C. Pluvinage G. Assessment of a gouge and dent defect in a pipeline by a combined criterion Engineering Failure Analysis 2014 36 1 13 10.1016/j.engfailanal.2013.10.002 2-s2.0-84886702955 Paeper S. Brown B. Beuker T. Inline inspection of dents and corrosion using ‘high quality’ multi-purpose smart-pig inspection data Proceedings of the 6th International Pipeline Conference (IPC '06) September 2006 243 248 10.1115/ipc2006-10157 2-s2.0-34547513613 Wang X. Song H. The inertial technology based 3-dimensional information measurement system for underground pipeline Measurement 2012 45 3 604 614 10.1016/j.measurement.2011.08.016 2-s2.0-84856758321 Mengjie L. The study of accurate in-line inspection technology to offshore pipeline route China Offshore Platform 2004 19 6 46 49 Yu J. Lee J. G. Park C. G. Han H. S. An off-line navigation of a geometry PIG using a modified nonlinear fixed-interval smoothing filter Control Engineering Practice 2005 13 11 1403 1411 10.1016/j.conengprac.2004.12.016 2-s2.0-21844453907 Czyz J. A. Fraccaroli C. Sergeant A. P. Measuring pipeline movement in geotechnically unstable areas using an inertial geometry pipeline inspection pig Proceedings of the ASME 1st International Pipeline Conference June 1996 Calgary, Canada Czyz J. A. Falk J. Use of geopig for prevention of pipeline failures in environmentally sensitive areas Proceedings of the Pipeline Pigging, Integrity Assessment and Repair Conference February 2000 Houston, Tex, USA Iflefel I. B. Moffat D. G. Mistry J. The interaction of pressure and bending on a dented pipe International Journal of Pressure Vessels and Piping 2005 82 10 761 769 10.1016/j.ijpvp.2005.06.002 2-s2.0-23844462252 Błachut J. Iflefel I. B. Collapse of pipes with plain or gouged dents by bending moment International Journal of Pressure Vessels and Piping 2007 84 9 560 571 10.1016/j.ijpvp.2007.04.007 2-s2.0-34548321739 Baek J.-H. Kim Y.-P. Kim W.-S. Koo J.-M. Seok C.-S. Load bearing capacity of API X65 pipe with dent defect under internal pressure and in-plane bending Materials Science and Engineering: A 2012 540 70 82 10.1016/j.msea.2012.01.078 2-s2.0-84862778705 Limam A. Lee L.-H. Kyriakides S. On the collapse of dented tubes under combined bending and internal pressure International Journal of Mechanical Sciences 2012 55 1 1 12 10.1016/j.ijmecsci.2011.10.005 2-s2.0-84858342115 Kim D. K. Cho S. H. Park S. S. Yoo H. R. Rho Y. W. Design and implementation of 30 geometry PIG KSME International Journal 2003 17 5 629 636 10.1007/BF02983858 Czyz J. A. Multi-pipeline geographical information system based on high accuracy inertial surveys Proceedings of the ASME 3rd International Pipeline Conference 2000 Calgary, Canada Hart J. D. Zulfiqar N. Moore D. H. Swank G. R. ‘Digital Pigging’ as a basis for improved pipeline structural integrity evaluations 2 Proceedings of the International Pipeline Conference September 2006 Calgary, Canada Integrity Management; Poster Session; Student Paper Competition 10.1115/IPC2006-10349 Czyz J. A. Adams J. R. Computations of pipeline—bending strains based on geopig measurements Proceedings of the ipeline Pigging and Integrity Monitoring Conference February 1994 Houston, Tex, USA 14 17 Hart J. D. Powel G. H. Geometry Monitoring of the Trans-Alaska Pipeline 2005 Trans-Alaska Pipeline System (TAPS) Yu D. Zhang L. Liang W. Ye Y. Wang Z. Noise reduction of signal and condition recognition of long-distance pipeline Acta Petrolei Sinica 2009 30 6 937 941 2-s2.0-72449162527 Zhou J. Gui B. Li M. Lin W. An application of the artificial neural net dominated by lithology to permeability prediction Acta Petrolei Sinica 2010 31 6 985 988 2-s2.0-78650506472 Li D. Lu D. Kong X. Du Y. Processing of well log data based on backpropagation neural network implicit approximation Acta Petrolei Sinica 2007 28 3 105 108 2-s2.0-34347268422 Peng K.-X. Yu S.-Z. Prediction and control of strip thickness based on Bayesian neural networks Journal of University of Science and Technology Beijing 2010 32 2 256 262 2-s2.0-77950567926 Yang H. Fu H.-Z. Stock index forecast based on bayesian regularization BP neural network Science Techno Logy and Engineering 2009 6 3306 3310 Wang H. Shi L. Zhao H. Yue Y. Temperature prediction model for sunlight greenhouse based on bayesian regularization BP neural network Hubei Agricultural Sciences 2015 9 4300 4303 Shenghua Y. Juan D. GDP prediction based on principal component analysis and bayesian regularization BP neural network Journal of Hunan University (Social Sciences) 2011 11 42 45 Li R. Feng Q. Cai M. Li H. Zhang H. Liu C. Zhao X. Measurement of long-distance buried pipeline centerline based on multi-sensor data fusion Acta Petrolei Sinica 2014 35 5 987 992 10.7623/syxb201405021 2-s2.0-84908018860