一种大体积混凝土结构隐患病害的定量识别方法及系统

1.本发明涉及一种大体积混凝土结构隐患病害的定量识别方法及系统，属于工程结构的健康监测技术领域。


背景技术：

2.随着大体积混凝土结构年龄的增长，在设计、施工及运行管理过程中存在的诸多不合理因素以及混凝土材料自身性能劣化等影响下，大体积混凝土结构不可避免地会暴露出先天和后天的缺陷或病害隐患。
3.基于波传播理论的常规无损检测方法(如：瑞利波检测技术、超声波检测技术、ct检测技术)已广泛应用于大体积混凝土结构的无损检测中，波在混凝土介质中传播时的观测信号携带有关混凝土材料性质、内部结构及组成的信息，然而受制于测试手段和技术，观测信号的解释是异常困难的，且检测结果仅能粗略判断缺陷是否存在，并不能较为准确地给出缺陷的具体位置和大小，期望通过对这些信号的解译定量地识别混凝土中的缺陷情况难度较大。
4.近年来发展的基于计算力学方法和智能优化算法的结合为结构内部缺陷的无损检测提供了一条新的途径，然而，该方法需要进行正分析和目标函数最小化迭代过程，在目标函数最小化迭代过程中易于陷入局部最优解，且迭代优化过程计算代价昂贵，此外，对于大体积混凝土结构由于体积庞大、结构形式异常复杂，对结构全域进行数值离散时在网格剖分及边界条件处理等方面难度极大。


技术实现要素：

5.本发明所要解决的技术问题是克服现有技术的缺陷，提供一种大体积混凝土结构隐患病害的定量识别方法及系统。
6.为解决上述技术问题，本发明提供一种大体积混凝土结构隐患病害的定量识别方法，包括：
7.获取实测的接收点信号xm(t)输入到预先构建的非线性映射关系模型中，输出对应的ym，ym即为表征真实结构内隐患病害的几何信息；
8.所述非线性映射关系模型的构建，包括：
9.对待检测大体积混凝土结构域进行识别，确定结构域中缺陷的存在方位；
10.根据结构域的几何维度以及所述缺陷的存在方位确定检测域的几何维度；
11.根据检测域范围，建立含边界吸收层的检测域的离散数值模型，并依据现场环境条件确定检测域中信号发射点和信号接收点的位置；
12.根据所述检测域中信号发射点和信号接收点的位置确定真实结构中的信号发射点和信号接收点的位置；
13.在真实结构中的信号发射点处输入波动信号f(t)，在真实结构中的信号接收点处采集接收点信号xn(t)；
14.在检测域中随机生成能够表征结构域内隐患病害的几何信息yi，根据所述波动信号f(t)模拟含隐患病害的检测域中的波传播过程，通过所述离散数值模型进行数值计算获得接收点的输出响应信号xi(t)，生成一对数据集(xi(t),yi)；通过反复执行生成数据集(xi(t),yi)的过程，直至生成满足条件的n对数据集(xi(t),yi),i＝1,2,
…
,n；
15.建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系模型。
16.进一步的，所述检测域几何维度至少有一条边界为结构域边界。
17.进一步的，所述边界吸收层为可模拟波在无限域中传播的吸收边界，且边界吸收层厚度aw为l
min
为检测域中边界的最小值。
18.进一步的，所述输入信号f(t)包括纵波主频fs和横波主频f
p
，所述纵波主频fs和横波主频f
p
的确定过程包括：
19.根据待检测大体积混凝土结构域的混凝土材料的质量密度ρ、弹性模量e和泊松比ν，通过下式计算混凝土介质中纵波波速cs和横波波速c
p
[0020][0021]
通过下式计算输入信号f(t)的两个主频fs和f
p
，
[0022]fs
＝cs/aw，f
p
＝c
p
/aw。
[0023]
进一步的，所述n对数据集中n的取值依据待检测的隐患病害的几何参数个数n确定，n＝η
×2n
，η取50到200之间的任意整数。
[0024]
进一步的，所述建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系模型，包括：
[0025]
采用数学模型建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系模型。
[0026]
进一步的，所述获取实测的接收点信号xn(t)输入到预先构建的非线性映射关系模型中，输出对应的yn，包括：
[0027]
反复将实测的接收点信号xm(t)输入到所述预先构建的非线性映射关系模型中，输入次数n≥2000，获取输出的n个yn值；
[0028]
计算置信水平为95％的n个yn样本中每个几何量的置信区间，取置信区间中值为表征真实结构内隐患病害的几何信息。
[0029]
进一步的，所述几何信息yi，包括：
[0030]
若结构内隐患病害为空洞，则采用几何形体拟合实际结构内的空洞；若结构内隐患病害为裂纹，则采用线段拟合实际结构内的裂纹。
[0031]
一种大体积混凝土结构隐患病害的定量识别系统，包括：
[0032]
获取模块，用于获取实测的接收点信号xm(t)；
[0033]
处理模块，用于将接收点信号xm(t)输入到预先构建的非线性映射关系模型中，输出对应的ym，ym即为表征真实结构内隐患病害的几何信息；
[0034]
所述处理模块包括：
[0035]
模型构建单元，用于对待检测大体积混凝土结构域进行识别，确定结构域中缺陷的存在方位；根据结构域的几何维度以及所述缺陷的存在方位确定检测域的几何维度；根
据检测域范围，建立含边界吸收层的检测域的离散数值模型，并依据现场环境条件确定检测域中信号发射点和信号接收点的位置；根据所述检测域中信号发射点和信号接收点的位置确定真实结构中的信号发射点和信号接收点的位置；在真实结构中的信号发射点处输入波动信号d(t)，在真实结构中的信号接收点处采集接收点信号xn(t)；在检测域中随机生成能够表征结构域内隐患病害的几何信息yi，根据所述波动信号f(t)模拟含隐患病害的检测域中的波传播过程，通过所述离散数值模型进行数值计算获得接收点的输出响应信号xi(t)，生成一对数据集(xi(t),yii；通过反复执行生成数据集(xi(t),yi)的过程，直至生成满足条件的n对数据集(xi(t),yi),i＝1,2,
…
,n；建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系模型。
[0036]
本发明所达到的有益效果：
[0037]
本发明通过结合现场测试和数值模拟可定量地检测出大体积混凝土结构中的隐患病害，本方法对结构本身无任何破坏作用，通过检测结果可诊断大体积混凝土结构中隐患病害的危险性等级，为进一步工程病险隐患处理措施的决策提供依据。
附图说明
[0038]
图1为本发明的定量识别方法的流程示意图；
[0039]
图2为含缺陷重力坝几何尺寸及检测域示意图；
[0040]
图3为检测域中边界吸收层、信号发射点及接收点示意图；
[0041]
图4为检测域有限元计算网格；
[0042]
图5为输入波动信号f(t)；
[0043]
图6为深度学习网络结构图；
[0044]
图7为根据置信区间中值确定的识别缺陷。
具体实施方式
[0045]
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0046]
一种大体积混凝土结构隐患病害的定量识别方法，包括以下步骤：
[0047]
步骤1.1)进行首轮现场测试，通过现场测试技术进行待检测大体积混凝土结构(结构域)的系统识别，确定结构域中缺陷的存在及存在缺陷的方位；
[0048]
步骤1.2)根据系统识别结果及结构域几何维度，确定检测域范围；
[0049]
步骤1.3)根据检测域的几何维度，建立含边界吸收层的检测域的离散数值模型，并依据现场环境条件确定检测域中信号发射点和接收点的位置；
[0050]
步骤1.4)进行二轮现场测试，在真实结构中，对应数值模拟中检测域的信号发射点在现场该点处输入波动信号f(t)，通过在接收点的位置粘贴的传感器记录测试得到的接收点信号xm(t)；
[0051]
步骤1.5)在检测域中随机生成能够表征结构域内隐患病害的几何信息yi(比如：若结构内隐患病害为空洞，可采用圆形、椭圆形、多边形等几何形体拟合实际结构内的空洞，圆形可通过圆心及半径即可唯一确定其几何信息，椭圆形可通过椭圆中心、长半轴、短半轴、方位角唯一确定其几何信息，多边形可通过多边形顶点位置唯一确定其几何信息；再
比如：若结构内隐患病害为裂纹，可通过裂纹的两个端点位置唯一确定其几何信息)，通过在检测域中信号发射点输入步骤1.4)所述波动信号f(t)，采用数值方法(如：有限元法、扩展有限元法、比例边界有限元法、无网格法等)模拟含隐患病害的检测域中的波传播过程，获得接收点的输出响应信号xi(t)，此时生成一对数据集(xi(t),yi)。通过反复执行上述步骤，在检测域中随机生成能够表征结构域内隐患病害的几何信息yi，并获得接收点的输出响应信号xi(t)，直至生成满足条件的n对数据集(xi(t),yi),i＝1,2,
…
,n；
[0052]
步骤1.6)采用数学模型建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系，即给定xi(t)，可输出对应的yi；
[0053]
步骤1.7)将步骤1.4)中实测的接收点信号xm(t)输入步骤1.6)的数学模型中，输出对应的ym，ym即为表征真实结构内隐患病害的几何信息。
[0054]
步骤1.1)中的现场测试技术可为超声波检测技术、瑞利波检测技术、ct检测技术、雷达探测等。
[0055]
步骤1.2)中的确定检测域时应至少有一条边界为结构域边界。
[0056]
步骤1.3)中的边界吸收层为可模拟波在无限域中传播的吸收边界且边界吸收层厚度aw为(l
min
为检测域中边界的最小值)。
[0057]
步骤1.4)中的现场测试技术为基于波传播理论的检测技术，如瑞利波检测技术、ct检测技术等。
[0058]
步骤1.4)中的输入信号f(t)的两个主频fs和f
p
确定方法为
[0059]
根据待检测大体积混凝土结构的混凝土材料的质量密度ρ、弹性模量e和泊松比ν，通过下式计算混凝土介质中纵波波速cs和横波波速c
p
[0060][0061]
输入信号f(t)的两个主频fs和f
p
可通过下式计算
[0062]fs
＝cs/aw，f
p
＝c
p
/aw[0063]
步骤1.5)中所述的n对数据集，n的确定与待检测的隐患病害的几何参数个数n有关，n＝η
×2n
，η可取50到200之间的任意整数。
[0064]
步骤1.6)中所述的数学模型应为具有统计意义的不确定性的数学模型，可为回归模型、基于人工神经网络的机器学习模型等。
[0065]
步骤1.7)中ym输出值为满足检测要求的置信区间，具体获取方式包括以下步骤：
[0066]
反复将步骤1.4)中实测的接收点信号xm(t)输入步骤1.6)所述的数学模型中，输入次数n≥2000，获取输出的n个ym值；
[0067]
计算置信水平为95％的n个ym样本中每个几何量的置信区间。
[0068]
实施例
[0069]
如图1所示为一种大体积混凝土结构隐患病害的定量识别方法的流程图，以大体积混凝土重力坝为例，其特征在于，包括以下步骤：
[0070]
步骤1.1)进行首轮现场测试，通过现场测试技术进行待检测大体积混凝土结构(结构域)的系统识别，确定结构域中缺陷的存在及存在缺陷的方位，如图1所示，通过探地雷达现场测试可初步判定在距重力坝坝顶30m～40m范围有一孔洞状缺陷的存在；
[0071]
步骤1.2)如图1，根据系统识别结果及结构域几何维度，确定检测域范围，检测域范围为abcde所围区域，ab边界为重力坝上游面部分区段，cde边界为重力坝下游面部分区段，ab边界和cde边界均为结构域边界；
[0072]
步骤1.3)根据检测域的几何维度，建立含边界吸收层的检测域的离散数值模型，检测域最小边界ab长度为l
min
＝10m，边界吸收层厚度取＝10m，边界吸收层厚度取并依据现场环境条件确定检测域中信号发射点和接收点的位置位于cde边界上，如图2所示，图3为数值计算时检测域的有限元网格；
[0073]
步骤1.4)进行二轮现场测试，在真实结构中，对应数值模拟中检测域的信号发射点在现场该点处输入波动信号f(t)，通过在接收点的位置粘贴的传感器记录测试得到的接收点信号xn(t)；
[0074]
混凝土的杨氏模量为e＝31.027gpa，泊松比ν＝0.15，质量密度ρ＝2643kg/m3，通过下式计算混凝土介质中纵波波速cs和横波波速c
p
[0075][0076]
输入信号f(t)的两个主频fs和f
p
可通过下式计算
[0077]fs
＝cs/aw＝1129.61hz，f
p
＝c
p
/aw＝1760.37hz
[0078]
如图5所示，输入信号f(t)由两主频为fs和f
p
的经汉宁窗过滤的正弦信号，可表达为
[0079][0080]
步骤1.5)用椭圆形拟合实际结构内的空洞，在检测域中随机生成能够表征结构域内隐患病害的几何信息yi＝{x
ei
,y
ei
,ai,bi,θi}，参数个数n＝5，(x
ei
,y
ei
)为椭圆中心点坐标，ai为椭圆长半轴，bi为椭圆短半轴，θi为椭圆方位角。通过在检测域中信号发射点输入图5所示波动信号f(t)，采用扩展有限元法仿真模拟含隐患病害的检测域中的波传播过程，获得接收点的输出响应信号xi(t)，此时生成一对数据集(xi(t),yi)。反复执行上述步骤，在检测域中随机生成能够表征结构域内隐患病害的几何信息yi，并获得接收点的输出响应信号xi(t)，生成满足条件的n＝3200对数据集(xi(t),yi),＝1,2,
…
,；满足条件，n＝η
×2n
，η取100。
[0081]
步骤1.6)如图6所示，建立深度学习网络结构模型，采用深度学习技术建立n＝3200对数据集中xi(t)和yi之间隐式的非线性映射关系，即通过将xi(t)输入深度学习模型，可输出对应的yi；
[0082]
步骤1.7)将实测的接收点信号xm(t)多次重复输入到深度学习模型中，输入2000次，输出对应的2000个不同的ym＝{x
em
,y
em
,am,bm,θm}，求出这2000个样本置信水平为95％的置信区间，如表1所示。取置信区间中值，即ym＝{14.86m,57.73m,1.95m,0.74m,12.75
°
}，识别的缺陷如图7所示。
[0083]
表1用于表征缺陷几何特征的物理量的置信区间
[0084][0085][0086]
相应的本发明还提供一种大体积混凝土结构隐患病害的定量识别系统，包括：
[0087]
获取模块，用于获取实测的接收点信号xm(t)；
[0088]
处理模块，用于将接收点信号xm(t)输入到预先构建的非线性映射关系模型中，输出对应的ym，ym即为表征真实结构内隐患病害的几何信息；
[0089]
所述处理模块包括：
[0090]
模型构建单元，用于对待检测大体积混凝土结构域进行识别，确定结构域中缺陷的存在方位；根据结构域的几何维度以及所述缺陷的存在方位确定检测域的几何维度；根据检测域范围，建立含边界吸收层的检测域的离散数值模型，并依据现场环境条件确定检测域中信号发射点和信号接收点的位置；根据所述检测域中信号发射点和信号接收点的位置确定真实结构中的信号发射点和信号接收点的位置；在真实结构中的信号发射点处输入波动信号f(t)，在真实结构中的信号接收点处采集接收点信号xn(t)；在检测域中随机生成能够表征结构域内隐患病害的几何信息yi，根据所述波动信号f(t)模拟含隐患病害的检测域中的波传播过程，通过所述离散数值模型进行数值计算获得接收点的输出响应信号xi(t)，生成一对数据集(xi(t),yi)；通过反复执行生成数据集(xi(t),yi)的过程，直至生成满足条件的n对数据集(xi(t),yi),i＝1,2,
…
,n；建立n对数据集中xi(t)和yi之间的显式或隐式的非线性映射关系模型。
[0091]
本发明的无损检测方法可应用于土木、水利工程中大体积混凝土结构的检测。本发明通过结合现场测试和数值模拟可定量地检测出大体积混凝土结构中的隐患病害，本方法对结构本身无任何破坏作用，通过检测结果可诊断大体积混凝土结构中隐患病害的危险性等级，为进一步工程病险隐患处理措施的决策提供依据。
[0092]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0093]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产
生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0094]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0095]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0096]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。