[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN110715929B - Distributed strain micro crack detection system and method based on stacking self-encoder - Google Patents

Distributed strain micro crack detection system and method based on stacking self-encoder Download PDF

Info

Publication number
CN110715929B
CN110715929B CN201910973652.6A CN201910973652A CN110715929B CN 110715929 B CN110715929 B CN 110715929B CN 201910973652 A CN201910973652 A CN 201910973652A CN 110715929 B CN110715929 B CN 110715929B
Authority
CN
China
Prior art keywords
strain
encoder
subsequences
subsequence
distributed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201910973652.6A
Other languages
Chinese (zh)
Other versions
CN110715929A (en
Inventor
宋青松
武金睿
陈禹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changan University
Original Assignee
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN201910973652.6A priority Critical patent/CN110715929B/en
Publication of CN110715929A publication Critical patent/CN110715929A/en
Application granted granted Critical
Publication of CN110715929B publication Critical patent/CN110715929B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于堆叠自编码器的分布式应变裂缝检测系统及方法,利用深度神经网络良好的特征表征能力,将裂缝检测视为一个二分类问题,构建一个基于堆叠自编码器的深度神经网络,实现结构体应变子序列的裂缝与非裂缝分类。本发明方法在实验室钢梁上可以准确并且无遗漏地检测张口宽度为32μm的微小裂缝,为结构体表面分布式应变裂缝检测提供了一种具有良好噪声鲁棒性的解决方案。

Figure 201910973652

The invention discloses a distributed strain crack detection system and method based on a stacked self-encoder, which utilizes the good feature representation ability of a deep neural network, regards crack detection as a two-class problem, and constructs a depth based on the stacked self-encoder. A neural network for crack and non-crack classification of structural strain subsequences. The method of the invention can accurately and without omission detect tiny cracks with an opening width of 32 μm on a laboratory steel beam, and provides a solution with good noise robustness for the detection of distributed strain cracks on the surface of the structure.

Figure 201910973652

Description

一种基于堆叠自编码器的分布式应变微小裂缝检测系统及 方法A distributed strain micro-crack detection system and method based on stacked autoencoders

技术领域technical field

本发明属于模式识别领域,具体涉及一种基于堆叠自编码器的分布式应变微小裂缝检测系统及方法。The invention belongs to the field of pattern recognition, and in particular relates to a distributed strain micro-crack detection system and method based on a stacked self-encoder.

背景技术Background technique

裂缝检测一直是结构健康监测领域中的重要课题。裂缝检测方法包含了人工观测的方法和无损检测方法。人工观测的方法需要专门的维护人员使用专业的工具进行定期检查,该方法效率低、主观性强。无损检测方法主要通过超声波、X射线、探地雷达以及摄像机等获得的数据对结构体裂缝进行检测。这些传感器都是点对点传感器,无法针对结构体的整体数据进行测量,容易遗漏裂缝。Crack detection has always been an important topic in the field of structural health monitoring. Crack detection methods include manual observation methods and non-destructive testing methods. The manual observation method requires specialized maintenance personnel to use professional tools to conduct regular inspections, which is inefficient and highly subjective. Non-destructive testing methods mainly detect structural cracks through the data obtained by ultrasonic, X-ray, ground penetrating radar and cameras. These sensors are all point-to-point sensors, which cannot measure the overall data of the structure and are prone to missing cracks.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于提供一种基于堆叠自编码器的分布式应变微小裂缝检测系统及方法,以克服上述现有技术存在的缺陷,本发明能够对结构体的整体应变数据进行有效的采集,结构体整体应变数据进行裂缝检测,显著改善了裂缝检测的正确率,提升裂缝检测的检测效果,为结构健康监测提供了一种高效的裂缝检测方案。The purpose of the present invention is to provide a distributed strain micro-crack detection system and method based on stacked self-encoders, so as to overcome the above-mentioned defects in the prior art, and the present invention can effectively collect the overall strain data of the structure. The overall strain data is used for crack detection, which significantly improves the accuracy of crack detection, improves the detection effect of crack detection, and provides an efficient crack detection solution for structural health monitoring.

为达到上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:

一种基于堆叠自编码器的分布式应变微小裂缝检测系统,包括:应变序列采集模块,用于对结构体表面的分布式应变进行采集;应变序列预处理模块:用于将采集所得的应变序列进行z-score标准化并截取为应变子序列;基于堆叠自编码器的特征自学习与表征模块:用于提取所划分的应变子序列的特征;Softmax分类识别模块,用于对提取到的子序列特征进行二分类,判别每一个子序列属于裂缝子序列和非裂缝子序列的概率;A distributed strain micro-crack detection system based on stacked self-encoders, comprising: a strain sequence acquisition module for acquiring distributed strain on the surface of a structure; a strain sequence preprocessing module: for collecting the acquired strain sequence Standardize z-score and truncate into strain subsequences; feature self-learning and characterization module based on stacked autoencoder: used to extract the features of the divided strain subsequences; Softmax classification and recognition module, used to extract the subsequences The features are classified into two categories to determine the probability that each subsequence belongs to a cracked subsequence and a non-cracked subsequence;

进一步地,应变序列采集模块:将光纤传感器敷设于结构体表面,使用基于BOTDA的分布式光纤传感系统对结构体表面的分布式应变进行采集;应变序列预处理模块包括:Further, the strain sequence acquisition module: lays the optical fiber sensor on the surface of the structure, and uses the distributed optical fiber sensing system based on BOTDA to collect the distributed strain on the surface of the structure; the strain sequence preprocessing module includes:

进一步地,z-score标准化模块和滑窗截取模块,z-score标准化模块将应变序列标准化为0均值1标准差的数据。滑窗模块通过长度为21,步长为1的滑动窗口将标准化后的应变序列截取了一组长度都是21的应变子序列。基于堆叠自编码器的特征自学习与表征模块:由3个自动编码器模块构成,将3个自动编码器模块的编码部分用作特征表征。Further, the z-score normalization module and the sliding window interception module, the z-score normalization module normalizes the strain series to data with 0 mean and 1 standard deviation. The sliding window module intercepts a set of strain subsequences with a length of 21 from the normalized strain sequence through a sliding window with a length of 21 and a step size of 1. Feature Self-Learning and Representation Module Based on Stacked Auto-Encoder: It consists of 3 auto-encoder modules, and the encoded part of the 3 auto-encoder modules is used as feature representation.

进一步地,基于堆叠自编码器的特征自学习与表征模块:由3个自动编码器模块构成,用于提取所划分的应变子序列的特征,自动编码器模块对于输入数据x,特征h与输出

Figure BDA0002232923400000021
之间的关系可以表示为fθ(·)和gθ'(·)两个函数,具体如下所示:Further, the feature self-learning and characterization module based on stacked autoencoders: composed of 3 autoencoder modules, used to extract the features of the divided strain subsequences, the autoencoder module for input data x, feature h and output
Figure BDA0002232923400000021
The relationship between can be expressed as two functions f θ ( ) and g θ' ( ), as follows:

h=fθ(x)=sf(Wx+bh)h=f θ (x)=s f (Wx+b h )

Figure BDA0002232923400000022
Figure BDA0002232923400000022

其中,W和bh分别为输入数据x和特征h之间的连接矩阵和偏置向量,WT和bv分别为特征h和输出

Figure BDA0002232923400000023
之间的连接矩阵和偏置向量,WT为W的转置矩阵。sf和sg是激活函数。Among them, W and b h are the connection matrix and bias vector between the input data x and feature h, respectively, W T and b v are the feature h and output, respectively
Figure BDA0002232923400000023
Between the connection matrix and the bias vector, W T is the transpose matrix of W. s f and s g are activation functions.

基于堆叠自编码器的特征自学习与表征模块的过程具体如下所示:The process of feature self-learning and representation module based on stacked autoencoder is as follows:

hi=sf(Wixi+bh,i)h i =s f (W i x i +b h,i )

Figure BDA0002232923400000024
Figure BDA0002232923400000024

其中,xi,hi

Figure BDA0002232923400000025
分别为第i(i=1,2,…,n)个自动编码器模块的输入数据x,特征h与输出
Figure BDA0002232923400000031
Wi和Wi T为其连接权重矩阵,bh,i和bv,i为其偏置向量,xi=hi-1,sf和sg为激活函数。Among them, x i , h i ,
Figure BDA0002232923400000025
are the input data x, feature h and output of the i-th (i=1,2,...,n) autoencoder module, respectively
Figure BDA0002232923400000031
Wi and Wi T are their connection weight matrices, b h,i and b v,i are their bias vectors, xi =h i-1 , s f and s g are activation functions.

一种基于堆叠自编码器的分布式应变裂缝检测方法,包括以下步骤:A distributed strain crack detection method based on stacked autoencoders, comprising the following steps:

步骤1:应变序列采集;Step 1: Strain sequence acquisition;

步骤2:使用z-score标准化将采集到的应变序列进行标准化,使用长度为21,步长为1的滑动窗口截取应变序列,得到应变子序列,将应变子序列按照截取位置进行标记;Step 2: Use z-score normalization to standardize the collected strain sequence, use a sliding window with a length of 21 and a step size of 1 to intercept the strain sequence to obtain a strain subsequence, and mark the strain subsequence according to the interception position;

步骤3:使用基于堆叠自编码器的神经网络自动学习表征应变子序列的特征;Step 3: Automatically learn features representing strain subsequences using a stacked autoencoder-based neural network;

步骤4:采用Softmax分类器实现对提取到的应变子序列的特征进行二分类,完成裂缝检测;Step 4: Use the Softmax classifier to implement binary classification of the features of the extracted strain subsequences to complete crack detection;

进一步地,步骤1中应变序列采集具体过程为:将光纤传感器通过环氧树脂粘附于结构体表面,光纤的两端接于基于BOTDA的分布式光纤传感系统,基于BOTDA的分布式光纤传感系统通过两个光源——泵浦光和探测光测量光纤的布里渊频移,通过布里渊频移与应变的线性关系得到结构体表面的分布式应变。Further, the specific process of strain sequence acquisition in step 1 is as follows: the optical fiber sensor is adhered to the surface of the structure through epoxy resin, and the two ends of the optical fiber are connected to the BOTDA-based distributed optical fiber sensing system, and the BOTDA-based distributed optical fiber transmission system. The sensing system measures the Brillouin frequency shift of the fiber through two light sources—pump light and probe light, and obtains the distributed strain on the surface of the structure through the linear relationship between the Brillouin frequency shift and strain.

进一步地,步骤2中应变序列处理的具体过程为:Further, the specific process of the strain sequence processing in step 2 is:

步骤2.1:将采集到的应变序列减去其均值,除以其标准差,得到0均值1标准差的数据。Step 2.1: Subtract the mean value of the acquired strain series and divide it by its standard deviation to obtain data with 0 mean and 1 standard deviation.

步骤2.2:使用长度为21,步长为1的滑动窗口沿着采集到的应变序列进行滑动,将应变序列截取为一组长度为21的应变子序列。Step 2.2: Use a sliding window with a length of 21 and a step of 1 to slide along the acquired strain sequence, and cut the strain sequence into a set of strain subsequences with a length of 21.

步骤2.3:将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列,其余标记为非裂缝子序列。Step 2.3: Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the left 3 strain subsequences and the right 4 strain subsequences as crack subsequences sequences, and the rest are marked as non-cracked subsequences.

进一步地,步骤3中使用基于堆叠自编码器的神经网络自动学习表征应变子序列特征的具体过程如下:Further, the specific process of using the stacked autoencoder-based neural network to automatically learn the characteristics of the strain subsequence in step 3 is as follows:

步骤3.1:模型初始化,确定模型的层数与神经元个数。随机初始化模型中的连接权重矩阵和偏置向量。输入层的神经元个数等于21,为应变子序列的长度。Step 3.1: Model initialization, determine the number of layers and neurons of the model. Randomly initialize the connection weight matrix and bias vector in the model. The number of neurons in the input layer is equal to 21, which is the length of the strain subsequence.

步骤3.2:预训练堆叠自编码器,堆叠自编码器由3个自动编码器构成,用得到的应变子序列预训练每个自动编码器。预训练自动编码器的损失函数为输入与输出之间的均方误差,具体如下:Step 3.2: Pretrain the stacked autoencoder, which consists of 3 autoencoders, and pretrain each autoencoder with the resulting strain subsequences. The loss function of the pretrained autoencoder is the mean squared error between the input and the output, as follows:

Figure BDA0002232923400000041
Figure BDA0002232923400000041

其中,x为输入的应变子序列,

Figure BDA0002232923400000042
为自动编码器输出的重构数据,M为所有输入的应变子序列的数量,Xm
Figure BDA0002232923400000043
分别是输入模型的第m条应变子序列和对应的输出重构的第m条子序列。where x is the input strain subsequence,
Figure BDA0002232923400000042
is the reconstructed data output by the autoencoder, M is the number of all input strain subsequences, X m ,
Figure BDA0002232923400000043
are the mth strain subsequence of the input model and the mth subsequence of the corresponding output reconstruction, respectively.

进一步地,步骤4中采用Softmax分类器实现对应变子序列的分类,具体方法为:Further, in step 4, the Softmax classifier is used to realize the classification of the strain subsequences, and the specific method is:

步骤4.1:构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:Step 4.1: Build a Softmax classifier. For a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, suppose The function h δ (z) outputs a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:

Figure BDA0002232923400000044
Figure BDA0002232923400000044

其中,δ12是Softmax分类器的全部参数,

Figure BDA0002232923400000045
z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier,
Figure BDA0002232923400000045
With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:

Figure BDA0002232923400000051
Figure BDA0002232923400000051

其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;

步骤4.2:预训练Sofmax,将应变子序列输入预训练后的堆叠自编码器,得到输出的特征z(i),以z(i)及其标签类别y(i)预训练Softmax,损失函数为交叉熵函数,具体如下:Step 4.2: Pre-train Sofmax, input the strain subsequence into the pre-trained stacked autoencoder, get the output feature z (i) , pre-train Softmax with z (i) and its label category y (i) , the loss function is Cross entropy function, as follows:

Figure BDA0002232923400000052
Figure BDA0002232923400000052

其中,

Figure BDA0002232923400000053
为Softmax分类器的全部参数,
Figure BDA0002232923400000054
为输出的类别概率,λ1为Softmax中连接权重矩阵和偏置向量正则项的权重系数,M为输入应变子序列的总个数,K为类别数,为2。in,
Figure BDA0002232923400000053
are all parameters of the Softmax classifier,
Figure BDA0002232923400000054
is the output category probability, λ 1 is the weight coefficient of the regular term connecting the weight matrix and the bias vector in Softmax, M is the total number of input strain subsequences, and K is the number of categories, which is 2.

步骤4.3:微调,堆叠自编码器的编码部分后接Softmax分类器,可以使其具有分类功能。利用预训练得到的应变子序列微调堆叠自编码器的编码部分与Softmax分类器整体结构的连接权重矩阵和偏置向量。微调时的损失函数为交叉损失函数,具体如下:Step 4.3: Fine-tuning, the encoding part of the stacked autoencoder is followed by a Softmax classifier, which can make it have a classification function. Using the pretrained strain subsequences to fine-tune the connection weight matrix and bias vector between the encoding part of the stacked autoencoder and the overall structure of the Softmax classifier. The loss function during fine-tuning is a cross loss function, as follows:

Figure BDA0002232923400000055
Figure BDA0002232923400000055

其中,ω为堆叠自编码器中的连接权重矩阵和偏置向量,Θ为ω和δ,λ2为堆叠自编码器中连接权重矩阵和偏置向量正则项的权重系数。where ω is the connection weight matrix and bias vector in the stacked autoencoder , Θ is ω and δ, and λ2 is the weight coefficient of the regular term of the connection weight matrix and bias vector in the stacked autoencoder.

步骤4.4:Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。Step 4.4: The Softmax classifier receives the features output by the stacked autoencoder as its input, and outputs the category 0 or 1 of the strain subsequence, where 0 means non-slit and 1 means crack; for the feature z (i) output by the stacked autoencoder, Select the category l with the largest probability p(y (i) =l|z (i) ; δ) as the category corresponding to the feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.

与现有技术相比,本发明具有以下有益的技术效果:Compared with the prior art, the present invention has the following beneficial technical effects:

本发明实现以分布式光纤传感器采集数据,改变以往点对点传感的方式。通过标准化缩小了不同数据之间的差异。同时,以基于堆叠自编码器的方法克服了分布式光纤传感器高空间分辨率与低信噪比的矛盾。堆叠自编码器能够在低信噪比的数据中提取高鲁棒性、可辨识的特征用以分类。在裂缝检测中效果显著,能够检测微小裂缝,在微小裂缝的检测效果上得提升。The invention realizes data collection by distributed optical fiber sensors, and changes the previous point-to-point sensing method. Standardization narrows the differences between different data. At the same time, the method based on stacked autoencoder overcomes the contradiction between high spatial resolution and low signal-to-noise ratio of distributed optical fiber sensor. Stacked autoencoders can extract highly robust and identifiable features for classification in low signal-to-noise ratio data. It has a significant effect in crack detection, can detect tiny cracks, and has improved the detection effect of tiny cracks.

附图说明Description of drawings

图1是本发明系统的流程示意图;Fig. 1 is the schematic flow chart of the system of the present invention;

图2是本发明中的自动编码器的示意图;Fig. 2 is the schematic diagram of the autoencoder in the present invention;

图3是本发明中的堆叠自编码器的示意图;Fig. 3 is the schematic diagram of the stacked autoencoder in the present invention;

图4是本发明方法的过程示意图;Fig. 4 is the process schematic diagram of the inventive method;

图5是本发明中具体的预训练与微调示意图。FIG. 5 is a schematic diagram of specific pre-training and fine-tuning in the present invention.

具体实施方式Detailed ways

下面结合附图对本发明作进一步详细描述:Below in conjunction with accompanying drawing, the present invention is described in further detail:

参见图1至图5,一种基于堆叠自编码器的分布式应变裂缝检测系统,包括应变序列采集模块;应变序列预处理模块;基于堆叠自编码器的特征自学习与表征模块;Softmax分类识别模块(具体流程如图1所示)。Referring to Figures 1 to 5, a distributed strain crack detection system based on stacked autoencoders includes a strain sequence acquisition module; a strain sequence preprocessing module; a feature self-learning and characterization module based on stacked autoencoders; Softmax classification and recognition module (the specific process is shown in Figure 1).

应变序列采集模块,用于采集结构体的分布式应变,采集到的结构体分布式应变为一个一维序列;The strain sequence acquisition module is used to collect the distributed strain of the structure, and the collected distributed strain of the structure is a one-dimensional sequence;

应变序列预处理模块包括:z-score标准化模块和滑窗模块,z-score标准化模块将应变序列标准化为0均值1标准差的数据。滑窗模块通过长度为21,步长为1的滑动窗口将标准化后的应变序列截取了一组长度都是21的应变子序列。将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列,其余标记为非裂缝子序列。The strain sequence preprocessing module includes: z-score normalization module and sliding window module. The z-score normalization module normalizes the strain sequence to data with 0 mean and 1 standard deviation. The sliding window module intercepts a set of strain subsequences with a length of 21 from the normalized strain sequence through a sliding window with a length of 21 and a step size of 1. Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the left 3 strain subsequences and the right four strain subsequences as the crack subsequence, and the rest marked as non-cracked subsequences.

基于堆叠自编码器的特征自学习与表征模块包括:3个自动编码器模块如附图2,自动编码器模块对于输入数据x,特征h与输出

Figure BDA0002232923400000071
之间的关系可以表示为fθ(·)和gθ'(·)两个函数,具体如下所示:The feature self-learning and characterization module based on stacked autoencoder includes: 3 autoencoder modules as shown in Figure 2, the autoencoder module for input data x, feature h and output
Figure BDA0002232923400000071
The relationship between can be expressed as two functions f θ ( ) and g θ' ( ), as follows:

h=fθ(x)=sf(Wx+bh)h=f θ (x)=s f (Wx+b h )

Figure BDA0002232923400000072
Figure BDA0002232923400000072

其中,W和bh分别为输入数据x和特征h之间的连接矩阵和偏置向量,WT和bv分别为特征h和输出

Figure BDA0002232923400000073
之间的连接矩阵和偏置向量,WT为W的转置矩阵。sf和sg是激活函数。Among them, W and b h are the connection matrix and bias vector between the input data x and feature h, respectively, W T and b v are the feature h and output, respectively
Figure BDA0002232923400000073
Between the connection matrix and the bias vector, W T is the transpose matrix of W. s f and s g are activation functions.

基于堆叠自编码器的特征自学习与表征模块的过程具体如下所示:The process of feature self-learning and representation module based on stacked autoencoder is as follows:

hi=sf(Wixi+bh,i)h i =s f (W i x i +b h,i )

Figure BDA0002232923400000074
Figure BDA0002232923400000074

其中,xi,hi

Figure BDA0002232923400000075
分别为第i(i=1,2,…,n)个自动编码器模块的输入数据x,特征h与输出
Figure BDA0002232923400000076
Wi和Wi T为其连接权重矩阵,bh,i和bv,i为其偏置向量,xi=hi-1,sf和sg为激活函数。Among them, x i , h i ,
Figure BDA0002232923400000075
are the input data x, feature h and output of the i-th (i=1,2,...,n) autoencoder module, respectively
Figure BDA0002232923400000076
Wi and Wi T are their connection weight matrices, b h,i and b v,i are their bias vectors, xi =h i-1 , s f and s g are activation functions.

采用Softmax分类器实现对应变子序列的分类,具体方法为:The classification of strain subsequences is realized by Softmax classifier, and the specific method is as follows:

构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:Construct a Softmax classifier. For a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, assuming the function h δ (z) Output a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:

Figure BDA0002232923400000077
Figure BDA0002232923400000077

其中,δ12是Softmax分类器的全部参数,

Figure BDA0002232923400000081
z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier,
Figure BDA0002232923400000081
With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:

Figure BDA0002232923400000082
Figure BDA0002232923400000082

其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;

Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。The Softmax classifier receives the features output by the stacked self-encoder as its input, and outputs the category 0 or 1 of the strain subsequence, where 0 means no crack and 1 means crack; for the feature z (i) output by the stacked self-encoder, the probability p is chosen (y (i) =l|z (i) ; δ) The largest category l is taken as the category corresponding to this feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.

一种基于堆叠自编码器的分布式应变裂缝检测方法,具体步骤如图5所示:A distributed strain crack detection method based on stacked autoencoders, the specific steps are shown in Figure 5:

1)、应变序列采集;1), strain sequence collection;

2)、使用z-score标准化将采集到的应变序列进行标准化,使用长度为21,步长为1的滑动窗口截取应变序列,得到应变子序列,将应变子序列按照截取位置进行标记;2), use z-score standardization to standardize the acquired strain sequence, use a sliding window with a length of 21 and a step size of 1 to intercept the strain sequence, obtain a strain subsequence, and mark the strain subsequence according to the interception position;

2.1:将采集到的应变序列减去其均值,除以其方差,得到0均值1标准差的数据。2.1: Subtract the mean value of the acquired strain series and divide it by its variance to obtain data with 0 mean and 1 standard deviation.

2.2:使用长度为21,步长为1的滑动窗口沿着采集到的应变序列进行滑动,将应变序列截取为一组长度为21的应变子序列。2.2: Use a sliding window with a length of 21 and a step size of 1 to slide along the acquired strain sequence, and cut the strain sequence into a set of strain subsequences with a length of 21.

2.3:将所得应变子序列按照截取位置标记标签,将以裂缝处为中心截取的应变子序列标记为裂缝子序列,并将其左侧3个和右侧4个应变子序列标记为裂缝子序列。2.3: Label the obtained strain subsequence according to the interception position, mark the strain subsequence intercepted at the crack as the crack subsequence, and mark the 3 left and 4 right strain subsequences as the crack subsequence .

3)、用基于堆叠自编码器的神经网络自动学习表征应变子序列的特征;3), use the neural network based on the stacked autoencoder to automatically learn the characteristics of the strain subsequence;

3.1:模型初始化,确定模型的层数与神经元个数。随机初始化模型中的连接权重矩阵和偏置向量。输入层的神经元个数等于21,为应变子序列的长度。3.1: Model initialization, determine the number of layers and neurons of the model. Randomly initialize the connection weight matrix and bias vector in the model. The number of neurons in the input layer is equal to 21, which is the length of the strain subsequence.

3.2:预训练堆叠自编码器,堆叠自编码器由3个自动编码器构成,用得到的应变子序列预训练每个自动编码器。预训练自动编码器的损失函数为输入与输出之间的均方误差,具体如下:3.2: Pre-training stacked autoencoders. The stacked autoencoders consist of 3 autoencoders, and each autoencoder is pretrained with the resulting strain subsequences. The loss function of the pretrained autoencoder is the mean squared error between the input and the output, as follows:

Figure BDA0002232923400000091
Figure BDA0002232923400000091

其中,x为输入的应变子序列,

Figure BDA0002232923400000092
为自动编码器输出的重构数据,M为所有输入的应变子序列的数量,Xm
Figure BDA0002232923400000093
分别是输入模型的第m条应变子序列和对应的输出重构的第m条子序列。where x is the input strain subsequence,
Figure BDA0002232923400000092
is the reconstructed data output by the autoencoder, M is the number of all input strain subsequences, X m ,
Figure BDA0002232923400000093
are the mth strain subsequence of the input model and the mth subsequence of the corresponding output reconstruction, respectively.

4)、采用Softmax分类器实现对提取到的应变子序列的特征进行二分类,完成裂缝检测;4), adopt the Softmax classifier to realize the two-classification of the features of the extracted strain subsequence, and complete the crack detection;

4.1:构建Softmax分类器,对于给定的输入z,用假设函数hδ(z)针对每一个类别l估算出概率值p(y=l|z),l∈{0,1},假设函数hδ(z)输出一个t维的向量表示这t个估计的概率值,t=2,假设函数hδ(z)如下:4.1: Build a Softmax classifier, for a given input z, use the hypothesis function h δ (z) to estimate the probability value p(y=l|z) for each category l, l∈{0,1}, the hypothesis function h δ (z) outputs a t-dimensional vector representing the t estimated probability values, t=2, assuming that the function h δ (z) is as follows:

Figure BDA0002232923400000094
Figure BDA0002232923400000094

其中,δ12是Softmax分类器的全部参数,

Figure BDA0002232923400000095
z(i)为输入,y(i)为输出,Softmax分类器将z分为类别l的概率为:Among them, δ 1 , δ 2 are all parameters of the Softmax classifier,
Figure BDA0002232923400000095
With z (i) as input and y (i) as output, the probability that the Softmax classifier classifies z into class l is:

Figure BDA0002232923400000101
Figure BDA0002232923400000101

其中,z(i)为输入,y(i)为输出;where z (i) is the input and y (i) is the output;

4.2:预训练Sofmax,将应变子序列输入预训练后的堆叠自编码器,得到输出的特征z(i),以z(i)及其标签类别y(i)预训练Softmax,损失函数为交叉熵函数,具体如下:4.2: Pre-training Sofmax, input the strain subsequence into the pre-trained stacked autoencoder, get the output feature z (i) , pre-train Softmax with z (i) and its label category y (i) , and the loss function is crossover The entropy function, as follows:

Figure BDA0002232923400000102
Figure BDA0002232923400000102

其中,

Figure BDA0002232923400000103
为Softmax分类器的全部参数,
Figure BDA0002232923400000104
为输出的类别概率,λ1为Softmax中连接权重矩阵和偏置向量正则项的权重系数,M为输入应变子序列的总个数,K为类别数,为2。in,
Figure BDA0002232923400000103
are all parameters of the Softmax classifier,
Figure BDA0002232923400000104
is the output category probability, λ 1 is the weight coefficient of the regular term connecting the weight matrix and the bias vector in Softmax, M is the total number of input strain subsequences, and K is the number of categories, which is 2.

4.3:微调,堆叠自编码器的编码部分后接Softmax分类器,可以使其具有分类功能。利用预训练得到的应变子序列微调堆叠自编码器的编码部分与Softmax分类器整体结构的连接权重矩阵和偏置向量。微调时的损失函数为交叉损失函数,具体如下:4.3: Fine-tuning, the encoding part of the stacked autoencoder is followed by a Softmax classifier, which can make it have a classification function. Using the pretrained strain subsequences to fine-tune the connection weight matrix and bias vector between the encoding part of the stacked autoencoder and the overall structure of the Softmax classifier. The loss function during fine-tuning is a cross loss function, as follows:

Figure BDA0002232923400000105
Figure BDA0002232923400000105

其中,ω为堆叠自编码器中的连接权重矩阵和偏置向量,Θ为ω和δ,λ2为堆叠自编码器中连接权重矩阵和偏置向量正则项的权重系数。where ω is the connection weight matrix and bias vector in the stacked autoencoder , Θ is ω and δ, and λ2 is the weight coefficient of the regular term of the connection weight matrix and bias vector in the stacked autoencoder.

4.4:Softmax分类器接收堆叠自编码器输出的特征作为其输入,输出应变子序列的类别0或1,0表示非裂缝,1表示裂缝;对于堆叠自编码器输出的特征z(i),选择概率p(y(i)=l|z(i);δ)最大的类别l作为该特征对应的类别。将分布式应变子序列还原至截取的位置且将位置相邻的裂缝子序列合并为一个裂缝。4.4: The Softmax classifier receives the features output by the stacked self-encoder as its input, and outputs the category 0 or 1 of the strain subsequence, 0 means non-crack and 1 means crack; for the feature z (i) output by the stacked self-encoder, choose The category l with the largest probability p(y (i) =l|z (i) ; δ) is used as the category corresponding to the feature. The distributed strain subsequence is restored to the truncated position and the adjacent fracture subsequences are merged into one fracture.

实施效果Implementation Effect

首先将光纤传感器预张紧,然后通过环氧树脂粘附于钢结构体表面。光纤传感器的两端接于基于BOTDA的分布式光纤传感系统,得到结构体表面沿光纤传感器径向分布的分布式应变数据。采用本发明基于堆叠自编码器的微小裂缝检测方法,基于采集的分布式应变数据可以准确无遗漏地检出张口宽度为32μm的微小裂缝,是一种可以用于钢结构表面微小裂缝分布式检测的有效方法。The fiber optic sensor is first pre-tensioned and then adhered to the surface of the steel structure by epoxy. The two ends of the optical fiber sensor are connected to the distributed optical fiber sensing system based on BOTDA, and the distributed strain data distributed along the radial direction of the optical fiber sensor on the surface of the structure body are obtained. By adopting the micro-crack detection method based on the stacked self-encoder of the present invention, the micro-cracks with an opening width of 32 μm can be accurately detected based on the collected distributed strain data, which is a distributed detection method for micro-cracks on the surface of steel structures. effective method.

Claims (5)

1. A distributed strain microcrack detection system based on stacked self-encoders, comprising:
strain sequence acquisition module: the system is used for acquiring distributed strain of the surface of the structure; the strain sequence acquisition module specifically comprises: laying an optical fiber sensor on the surface of a structure, and collecting distributed strain on the surface of the structure by using a distributed optical fiber sensing system based on BOTDA;
a strain sequence preprocessing module: the strain acquisition device is used for performing z-score standardization on the acquired distributed strain and intercepting the distributed strain into a strain subsequence; the strain sequence pre-processing module comprises a z-score normalization module and a sliding window module: the z-score normalization module normalizes the strain sequences to 0-mean 1 standard deviation data; the sliding window module intercepts a group of strain subsequences with the lengths of 21 from the normalized strain sequence through a sliding window with the length of 21 and the step length of 1;
the self-learning and characterization module of the characteristics based on the stacking self-encoder comprises the following modules: the system comprises 3 automatic encoder modules, a data processing module and a data processing module, wherein the automatic encoder modules are used for extracting the characteristics of the divided strain subsequences, inputting the characteristics into the strain subsequences and outputting the characteristics into the strain subsequences;
and the Softmax classification and identification module is used for judging the probability that each strain subsequence belongs to the crack subsequence and the non-crack subsequence so as to carry out secondary classification on the extracted characteristics of the strain subsequences.
2. The distributed strain microcrack detection system according to claim 1, wherein the self-learning and characterization module based on the self-encoder comprises 3 autoencoder modules for extracting the characteristics of the divided strain subsequences, and the autoencoder modules input data x, characteristics h and output data x, h and output data h
Figure FDA0003436420230000011
The relationship between is represented as fθ(. and g)θ'Two functions, specifically as follows:
h=fθ(x)=sf(Wx+bh)
Figure FDA0003436420230000012
wherein W and bhRespectively, a connection matrix and an offset vector, W, between the input data x and the feature hTAnd bvRespectively characteristic h and output
Figure FDA0003436420230000021
A connection matrix and an offset vector between, WTIs a transposed matrix of W, sfAnd sgIs an activation function;
the process of the feature self-learning and characterization module based on the stacked self-encoder is specifically as follows:
hi=sf(Wixi+bh,i)
Figure FDA0003436420230000022
wherein x isi,hi
Figure FDA0003436420230000023
Input data x, characteristic h and output of the ith automatic encoder module respectively
Figure FDA0003436420230000024
And i is 1,2, L, n, WiAnd Wi TTo which a weight matrix is connected, bh,iAnd bv,iIs its offset vector, xi=hi-1,sfAnd sgIs an activation function.
3. A distributed strain micro crack detection method based on a stacked self-encoder is characterized by comprising the following steps:
step 1: collecting distributed strain on the surface of the structure; the strain sequence acquisition specific process comprises the following steps: the method comprises the following steps that an optical fiber sensor is adhered to the surface of a structure body through epoxy resin, two ends of an optical fiber are connected to a BOTDA-based distributed optical fiber sensing system, the BOTDA-based distributed optical fiber sensing system measures Brillouin frequency shift of the optical fiber through two light source pumping light and detection light, and distributed strain of the surface of the structure body is obtained through the linear relation between the Brillouin frequency shift and strain;
step 2: carrying out z-score standardization on the acquired distributed strain and intercepting the distributed strain into a strain subsequence; the method specifically comprises the following steps:
step 2.1: subtracting the mean value of the acquired strain sequence, and dividing the mean value by the standard deviation to obtain data of 0 mean value 1 standard deviation;
step 2.2: intercepting the normalized strain sequence into a group of strain subsequences with the length of 21 by using a sliding window with the length of 21 and the step size of 1;
step 2.3: marking the obtained strain subsequences as cut position mark labels, marking the strain subsequences which are cut by taking the crack as the center as crack subsequences, marking 3 strain subsequences on the left side and 4 strain subsequences on the right side as crack subsequences, and marking the rest as non-crack subsequences;
and step 3: automatically learning features characterizing the strain subsequence using a stacked self-encoder based neural network;
and 4, step 4: and (4) performing secondary classification on the extracted characteristics of the strain subsequence by adopting a Softmax classifier, and completing crack detection.
4. The distributed strain microcrack detection method based on the stacked self-encoder according to claim 3, wherein the specific process of using the neural network based on the stacked self-encoder to automatically learn and characterize the strain subsequence features in step 3 is as follows:
step 3.1: initializing a model, determining the number of layers and the number of neurons of the model, randomly initializing a connection weight matrix and a bias vector in the model, and inputting the number of the neurons of the layer to be equal to the length of a strain subsequence;
step 3.2: pre-training a stacked self-encoder, the stacked self-encoder consisting of 3 auto-encoders, pre-training each auto-encoder with the obtained strain sub-sequence, the loss function of the pre-trained auto-encoder being the mean square error L between input and output1The method comprises the following steps:
Figure FDA0003436420230000031
wherein x is an input strain subsequence,
Figure FDA0003436420230000032
for the reconstructed data output from the autoencoder, M is the number of all the input strain subsequences, Xm
Figure FDA0003436420230000033
Respectively the mth strain subsequence of the input model and the corresponding mth subsequence of the output reconstruction.
5. The distributed strain micro crack detection method based on the stacked self-encoder as claimed in claim 3, wherein a Softmax classifier is adopted in the step 4 to realize classification of the strain subsequences, and the specific method is as follows:
step 4.1: constructing a Softmax classifier using a hypothesis function h for a given input zδ(z) for each class l, a probability value p (y ═ l | z) is estimated, l ∈ {0,1}, assuming a function hδ(z) outputting a vector of dimensions t representing the probability values of the t estimates, t being 2, assuming the function hδ(z) is as follows:
Figure FDA0003436420230000041
wherein,
Figure FDA0003436420230000042
δ12is all the parameters of the Softmax classifier, z(i)To input, y(i)For output, the probability of the Softmax classifier classifying z into class l is:
Figure FDA0003436420230000043
wherein z is(i)To input, y(i)Is an output;
step 4.2: pre-training Sofmax, inputting the strain subsequence into the pre-trained stacked self-encoder to obtain the output characteristic z(i)In z is(i)And its label category y(i)Pre-training Softmax, wherein a loss function is a cross entropy function, and the method comprises the following specific steps:
Figure FDA0003436420230000044
wherein,
Figure FDA0003436420230000045
is the class probability of the output, λ1The weight coefficient is a weight coefficient connecting a weight matrix and a bias vector regular term in Softmax, M is the total number of input strain subsequences, and K is the category number which is 2;
step 4.3: utilizing a pre-trained strain subsequence to fine-tune a connection weight matrix and a bias vector of an encoding part of the stacked self-encoder and an overall structure of the Softmax classifier, wherein a loss function during fine tuning is a cross loss function, and the method specifically comprises the following steps:
Figure FDA0003436420230000046
where ω is in a stacked self-encoderConnecting the weight matrix and the offset vector, theta is omega and delta, lambda2The weight coefficient is a weight coefficient which is connected with the weight matrix and the bias vector regular term in the stacked self-encoder;
step 4.4: the Softmax classifier receives as its input the features stacked from the encoder output, outputs class 0 or 1 of the strain subsequence, 0 representing non-crack, 1 representing crack; for feature z of stacked self-encoder output(i)Selecting the probability p (y)(i)=l|z(i)(ii) a δ) the largest class i as the class to which the feature corresponds.
CN201910973652.6A 2019-10-14 2019-10-14 Distributed strain micro crack detection system and method based on stacking self-encoder Expired - Fee Related CN110715929B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910973652.6A CN110715929B (en) 2019-10-14 2019-10-14 Distributed strain micro crack detection system and method based on stacking self-encoder

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910973652.6A CN110715929B (en) 2019-10-14 2019-10-14 Distributed strain micro crack detection system and method based on stacking self-encoder

Publications (2)

Publication Number Publication Date
CN110715929A CN110715929A (en) 2020-01-21
CN110715929B true CN110715929B (en) 2022-03-25

Family

ID=69211557

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910973652.6A Expired - Fee Related CN110715929B (en) 2019-10-14 2019-10-14 Distributed strain micro crack detection system and method based on stacking self-encoder

Country Status (1)

Country Link
CN (1) CN110715929B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111750795A (en) * 2020-06-18 2020-10-09 哈尔滨工程大学 A distributed creep measurement system and measurement method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102322808B (en) * 2011-08-09 2013-03-27 中国计量学院 Very long range pulse coding distribution type Fiber Raman and Brillouin photon sensor
CN202195825U (en) * 2011-08-09 2012-04-18 中国计量学院 Extra-long distance pulse-coding distributed optical fiber Raman and Brillouin photon sensor
US10846447B2 (en) * 2016-04-29 2020-11-24 Exxonmobil Upstream Research Company Method and system for stacking fracture prediction

Also Published As

Publication number Publication date
CN110715929A (en) 2020-01-21

Similar Documents

Publication Publication Date Title
CN109765053B (en) Rolling bearing fault diagnosis method using convolutional neural network and kurtosis index
Iyer et al. Ultrasonic signal processing methods for detection of defects in concrete pipes
CN107144643B (en) A kind of damnification recognition method of Lamb wave monitoring signals statistical parameter
CN102928435A (en) Aircraft skin damage identification method and device based on image and ultrasound information fusion
CN114782753A (en) Lung cancer histopathological whole section classification method based on weakly supervised learning and converter
CN119023677A (en) A highly accurate intelligent detection system for wiring harness defects
CN108537790A (en) Heterologous image change detection method based on coupling translation network
CN110715929B (en) Distributed strain micro crack detection system and method based on stacking self-encoder
CN110738168A (en) A distributed strain micro-crack detection system and method based on stacked convolutional autoencoders
Chen et al. A mixed samples-driven methodology based on denoising diffusion probabilistic model for identifying damage in carbon fiber composite structures
CN114358189A (en) Fault diagnosis method of hydraulic waterproof valve based on multimodal deep residual shrinkage network
CN119086721B (en) Plastic product strength detection method and device based on non-contact measurement
TW202144763A (en) Tomography method, system and apparatus based on time-domain spectroscopy
Banerjee et al. Damage detection and localization by learning deep features of elastic waves in piezoelectric ceramic using point contact method
CN118447000A (en) Welding defect detection method and system
Masaki et al. Surface undulation detection system using wearable artificial skin layer with strain gauge
CN112446612B (en) Assessment method of damage assessment system of soft rigid arm mooring system connection structure
CN116818574A (en) Material mechanical properties measurement method based on terahertz spectrum characteristic parameter characterization
Liu et al. Automatic terahertz recognition of hidden defects in layered polymer composites based on a deep residual network with transfer learning
Lin et al. Automatic elimination of invalid impact-echo signals for detecting delamination in concrete bridge decks based on deep learning
CN119150120A (en) Impact damage identification method based on multi-mode data fusion and ensemble learning
CN114104328A (en) Aircraft state monitoring method based on deep migration learning
CN106679886A (en) Nonlinear fault detecting and identifying method of self-confirming air data system
CN118777562B (en) A steel quality analysis and feedback system
CN117274691B (en) Hyperspectral image classification method based on space pooling transducer

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220325

CF01 Termination of patent right due to non-payment of annual fee