CN115687854B - 高精度土样参数测量方法及其系统 - Google Patents

高精度土样参数测量方法及其系统 Download PDF

Info

Publication number
CN115687854B
CN115687854B CN202310005060.1A CN202310005060A CN115687854B CN 115687854 B CN115687854 B CN 115687854B CN 202310005060 A CN202310005060 A CN 202310005060A CN 115687854 B CN115687854 B CN 115687854B
Authority
CN
China
Prior art keywords
matrix
scale
decoding
soil sample
feature
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.)
Active
Application number
CN202310005060.1A
Other languages
English (en)
Other versions
CN115687854A (zh
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.)
Shandong Institute Of Earthquake Engineering Co ltd
Original Assignee
Shandong Institute Of Earthquake Engineering Co ltd
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 Shandong Institute Of Earthquake Engineering Co ltd filed Critical Shandong Institute Of Earthquake Engineering Co ltd
Priority to CN202310005060.1A priority Critical patent/CN115687854B/zh
Publication of CN115687854A publication Critical patent/CN115687854A/zh
Application granted granted Critical
Publication of CN115687854B publication Critical patent/CN115687854B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

公开了一种高精度土样参数测量方法及其系统，其采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。

Description

高精度土样参数测量方法及其系统
技术领域
本申请涉及工程建设领域，且更为具体地，涉及一种高精度土样参数测量方法及其系统。
背景技术
对于土层反应分析中，土的动剪切模量和阻尼是重要的计算参数，为获取较为准确的数值需要根据相关理论对离散数据进行计算。
现有技术中，对剪切模量和阻尼的计算较为粗糙，计算过程较为繁琐，计算结果不甚理想。如阻尼的计算需要将数据模糊成一个标准图形，这样计算简单但结果粗糙；或者将数据在cad中制图后求取面积，这样结果虽然较为准确但计算过程繁琐。同时，由于室内动三轴试验的原理，反算的最大动剪切模量数值较为离散，需要利用现场较容易获得的剪切波速或者密度、含水率等进行修正。这样，不仅过程较为繁琐，而且计算出的数值精准度也难以保证。
因此，期望一种高精度土样参数测量方案。
发明内容
为了解决上述技术问题，提出了本申请。本申请的实施例提供了一种高精度土样参数测量方法及其系统，其采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
根据本申请的一个方面，提供了一种高精度土样参数测量方法，其包括：获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
在上述高精度土样参数测量方法中，所述计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵，包括：计算所述第一输入向量的转置向量与所述第二输入向量之间的向量乘积以得到所述关联参数矩阵。
在上述高精度土样参数测量方法中，所述将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵，包括：使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在上述高精度土样参数测量方法中，所述将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵，包括：使用所述第二卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第二卷积神经网络模型的最后一层的输出为所述第二尺度关联参数特征矩阵，所述第二卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在上述高精度土样参数测量方法中，所述融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵，包括：以如下公式来融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；其中，所述公式为：
Figure 58178DEST_PATH_IMAGE001
其中，
Figure 95405DEST_PATH_IMAGE002
为所述解码特征矩阵，
Figure 577201DEST_PATH_IMAGE003
为所述第一尺度关联参数特征矩阵，
Figure 787603DEST_PATH_IMAGE004
为所述第二尺度关联参数特征矩阵，“
Figure 405666DEST_PATH_IMAGE005
”表示所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵相对应位置处的元素相加，
Figure 676110DEST_PATH_IMAGE006
和
Figure 379624DEST_PATH_IMAGE007
为用于控制所述解码特征矩阵中所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵之间的平衡的加权参数。
在上述高精度土样参数测量方法中，所述对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵，包括：以如下公式对所述解码特征矩阵进行特征分布校正以得到所述校正后解码特征矩阵；其中，所述公式为：
Figure 862558DEST_PATH_IMAGE008
其中
Figure 663024DEST_PATH_IMAGE009
是所述解码特征矩阵的第
Figure 42053DEST_PATH_IMAGE010
位置的特征值，
Figure 498442DEST_PATH_IMAGE011
和
Figure 41857DEST_PATH_IMAGE012
是所述解码特征矩阵的宽度和高度，且
Figure 368933DEST_PATH_IMAGE013
表示以2为底的对数函数值，
Figure 981180DEST_PATH_IMAGE014
是所述校正后解码特征矩阵的第
Figure 456024DEST_PATH_IMAGE010
位置的特征值。
在上述高精度土样参数测量方法中，所述将所述解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量，包括：使用解码器以如下公式对所述解码特征矩阵进行解码回归以获得所述解码值；其中，所述公式为：
Figure 546339DEST_PATH_IMAGE015
，其中
Figure 993501DEST_PATH_IMAGE016
是所述解码特征矩阵，
Figure 511070DEST_PATH_IMAGE017
是所述解码值，
Figure 676472DEST_PATH_IMAGE018
是权重矩阵，
Figure 570479DEST_PATH_IMAGE019
表示矩阵乘。
根据本申请的另一方面，提供了一种高精度土样参数测量系统，包括：数据获取模块，用于获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；数据结构化模块，用于将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；关联模块，用于计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；第一尺度编码模块，用于将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；第二尺度编码模块，用于将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；融合模块，用于融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；特征分布校正模块，用于对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和解码模块，用于将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
在上述高精度土样参数测量系统中，所述关联模块，用于计算所述第一输入向量的转置向量与所述第二输入向量之间的向量乘积以得到所述关联参数矩阵。
在上述高精度土样参数测量系统中，所述第一尺度编码模块，进一步用于：使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在上述高精度土样参数测量系统中，所述第二尺度编码模块，进一步用于：使用所述第二卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第二卷积神经网络模型的最后一层的输出为所述第二尺度关联参数特征矩阵，所述第二卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在上述高精度土样参数测量系统中，所述融合模块，进一步用于：以如下公式来融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；其中，所述公式为：
Figure 137726DEST_PATH_IMAGE020
其中，
Figure 563547DEST_PATH_IMAGE021
为所述解码特征矩阵，
Figure 278562DEST_PATH_IMAGE022
为所述第一尺度关联参数特征矩阵，
Figure 913943DEST_PATH_IMAGE023
为所述第二尺度关联参数特征矩阵，“
Figure 335697DEST_PATH_IMAGE024
”表示所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵相对应位置处的元素相加，
Figure 460648DEST_PATH_IMAGE025
和
Figure 335063DEST_PATH_IMAGE026
为用于控制所述解码特征矩阵中所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵之间的平衡的加权参数。
在上述高精度土样参数测量系统中，所述特征分布校正模块，进一步用于：以如下公式对所述解码特征矩阵进行特征分布校正以得到所述校正后解码特征矩阵；其中，所述公式为：
Figure 570872DEST_PATH_IMAGE027
其中
Figure 643871DEST_PATH_IMAGE028
是所述解码特征矩阵的第
Figure 142985DEST_PATH_IMAGE029
位置的特征值，
Figure 567013DEST_PATH_IMAGE030
和
Figure 547126DEST_PATH_IMAGE031
是所述解码特征矩阵的宽度和高度，且
Figure 943472DEST_PATH_IMAGE032
表示以2为底的对数函数值，
Figure 613488DEST_PATH_IMAGE033
是所述校正后解码特征矩阵的第
Figure 259233DEST_PATH_IMAGE029
位置的特征值。
在上述高精度土样参数测量系统中，所述解码模块，进一步用于：使用解码器以如下公式对所述解码特征矩阵进行解码回归以获得所述解码值；其中，所述公式为：
Figure 368003DEST_PATH_IMAGE034
，其中
Figure 884435DEST_PATH_IMAGE035
是所述解码特征矩阵，
Figure 459773DEST_PATH_IMAGE036
是所述解码值，
Figure 61656DEST_PATH_IMAGE037
是权重矩阵，
Figure 442959DEST_PATH_IMAGE038
表示矩阵乘。
与现有技术相比，本申请提供的高精度土样参数测量方法及其系统，其采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
附图说明
通过结合附图对本申请实施例进行更详细的描述，本申请的上述以及其他目的、特征和优势将变得更加明显。附图用来提供对本申请实施例的进一步理解，并且构成说明书的一部分，与本申请实施例一起用于解释本申请，并不构成对本申请的限制。在附图中，相同的参考标号通常代表相同部件或步骤。
图1为根据本申请实施例的高精度土样参数测量方法的流程图。
图2为根据本申请实施例的高精度土样参数测量方法的架构图。
图3为根据本申请实施例的高精度土样参数测量系统的框图。
图4为根据本申请实施例的应力应变滞回曲线图。
图5为根据本申请实施例的最大动剪切模量对比图。
图6为根据本申请实施例的剪切模量比结果对比图。
图7为根据本申请实施例的动阻尼结果对比图。
具体实施方式
下面，将参考附图详细地描述根据本申请的示例实施例。显然，所描述的实施例仅仅是本申请的一部分实施例，而不是本申请的全部实施例，应理解，本申请不受这里描述的示例实施例的限制。
申请概述:如上述背景技术所言，对于土层反应分析中，土的动剪切模量和阻尼是重要的计算参数，为获取较为准确的数值需要根据相关理论对离散数据进行计算。
现有技术中，对剪切模量和阻尼的计算较为粗糙，计算过程较为繁琐，计算结果不甚理想。如阻尼的计算需要将数据模糊成一个标准图形，这样计算简单但结果粗糙；或者将数据在cad中制图后求取面积，这样结果虽然较为准确但计算过程繁琐。同时，由于室内动三轴试验的原理，反算的最大动剪切模量数值较为离散，需要利用现场较容易获得的剪切波速或者密度、含水率等进行修正。这样，不仅过程较为繁琐，而且计算出的数值精准度也难以保证。因此，期望一种高精度土样参数测量方案。
相应地，由于在实际工程试验中，通常根据土样在各级荷载下轴向应力和轴向应变时程曲线绘出应力应变滞回曲线，由滞回曲线的应力应变值可求得试样的轴向模量。也就是，考虑到在剪切模量的测量计算过程中，土样在各级荷载下轴向应力值和轴向应变时程值对于土样的剪切模量测量具有着关联性的影响，因此，可以基于所述土样在各级荷载下轴向应力值和轴向应变时程值在时序上的关联性特征分布信息来进行所述待测试土样的剪切模量测量。具体地，在本申请的技术方案中，采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
具体地，在本申请的技术方案中，首先，获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值。然后，将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量以整合所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序上的分布信息。接着，进一步再计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵，以构建所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的关联性分布表示。
进一步地，使用在隐含特征提取方面具有优异表现的卷积神经网络模型来进行所述关联参数矩阵中关于所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性隐含特征提取。特别地，考虑到在所述关联参数矩阵中关于所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性分布信息在时间维度和参数样本维度上具有不同的尺度关联，也就是说，具有所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性分布信息的所述关联参数矩阵中的各个特征值间在高维特征空间中具有多尺度的隐含关联特征。因此，在本申请的技术方案中，为了提高对于所述待测试土样的剪切模量测量的准确性，使用具有不同尺度的卷积核的卷积神经网络模型来进行特征挖掘，来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序上的关联性分布信息的多尺度隐含关联特征。具体地，将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵，并将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵。
然后，再进一步融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵，以融合所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序维度上的关联性特征分布表示的多尺度隐含关联特征信息，并以此作为解码特征矩阵通过解码器中进行解码回归以得到用于表示待测试土样的剪切模量的解码值。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
特别地，在本申请的技术方案中，考虑到所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵分别表达不同尺度下的所述轴向应力值和所述轴向应变时程值的时序关联参数在高维特征空间下的关联模式表达特征，所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵在高维特征空间内的特征分布既有同相重合又有异相差异，这使得当通过点加的方式融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵时，由于同相增强和异相削弱，所述解码特征矩阵的整体特征分布的离散程度高，在通过解码器进行解码时导致解码回归的归纳发散，从而影响通过解码器进行解码时的收敛性，影响解码值的准确性。
因此，在通过解码器进行解码之前，首先对所述解码特征矩阵进行面向回归有界闭域的分布转移优化，表示为：
Figure 813897DEST_PATH_IMAGE039
Figure 560136DEST_PATH_IMAGE040
是所述解码特征矩阵的第
Figure 917824DEST_PATH_IMAGE041
位置的特征值，
Figure 102818DEST_PATH_IMAGE042
和
Figure 328263DEST_PATH_IMAGE043
是所述解码特征矩阵的宽度和高度，且
Figure 245403DEST_PATH_IMAGE044
表示以2为底的对数。
也就是，所述面向回归有界闭域的分布转移优化针对所述解码特征矩阵所代表的高维特征分布在转移到解码回归问题的目标域时可能存在的归纳发散，通过基于条件化解码回归的结构化信息约束来向着特征集合的有界闭域进行特征分布的收敛，以使得所述解码特征矩阵的特征分布转移到目标域下的具有稳定的可结构化边界的范围内，提升了回归解的归纳迭代的稳定性，即，提高了所述解码特征矩阵通过解码器进行解码时的收敛性，改进了解码值的准确性。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
基于此，本申请提出了一种高精度土样参数测量方法，其包括：获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和，将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
在介绍了本申请的基本原理之后，下面将参考附图来具体介绍本申请的各种非限制性实施例。
实施例1:图1为根据本申请实施例的高精度土样参数测量方法的流程图。如图1所示，根据本申请实施例的高精度土样参数测量方法，包括：S110，获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；S120，将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；S130，计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；S140，将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；S150，将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；S160，融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；S170，对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和S180，将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
图2为根据本申请实施例的高精度土样参数测量方法的架构图。如图2所示，在该架构中，首先，获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值。接着，将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量。然后，计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵。进而，将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵，同时，将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵。接着，融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵。然后，对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵。进而，将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
在步骤S110中，获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值。如上述背景技术所言，对于土层反应分析中，土的动剪切模量和阻尼是重要的计算参数，为获取较为准确的数值需要根据相关理论对离散数据进行计算。现有技术中，对剪切模量和阻尼的计算较为粗糙，计算过程较为繁琐，计算结果不甚理想。如阻尼的计算需要将数据模糊成一个标准图形，这样计算简单但结果粗糙；或者将数据在cad中制图后求取面积，这样结果虽然较为准确但计算过程繁琐。同时，由于室内动三轴试验的原理，反算的最大动剪切模量数值较为离散，需要利用现场较容易获得的剪切波速或者密度、含水率等进行修正。这样，不仅过程较为繁琐，而且计算出的数值精准度也难以保证。因此，期望一种高精度土样参数测量方案。
相应地，由于在实际工程试验中，通常根据土样在各级荷载下轴向应力和轴向应变时程曲线绘出应力应变滞回曲线，由滞回曲线的应力应变值可求得试样的轴向模量。也就是，考虑到在剪切模量的测量计算过程中，土样在各级荷载下轴向应力值和轴向应变时程值对于土样的剪切模量测量具有着关联性的影响，因此，可以基于所述土样在各级荷载下轴向应力值和轴向应变时程值在时序上的关联性特征分布信息来进行所述待测试土样的剪切模量测量。具体地，在本申请的技术方案中，采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。具体地，在本申请的技术方案中，首先，获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值。其中，所述轴向应力值和所述轴向应变时程值可以由应力传感器和应变传感器来采集获取。
在步骤S120中，将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量。也就是，将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量以整合所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序上的分布信息。此外，将所述轴向应力值和所述轴向应变时程值进行数据结构化也利于后续模型的计算。
在步骤S130中，计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵。也就是，计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵，以构建所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的关联性分布表示。具体地，在本申请的一个实施例中，计算所述第一输入向量的转置向量与所述第二输入向量之间的向量乘积以得到所述关联参数矩阵。
在步骤S140中，将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵。也就是，使用在隐含特征提取方面具有优异表现的卷积神经网络模型来进行所述关联参数矩阵中关于所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性隐含特征提取。
特别地，考虑到在所述关联参数矩阵中关于所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性分布信息在时间维度和参数样本维度上具有不同的尺度关联，也就是说，具有所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值的关联性分布信息的所述关联参数矩阵中的各个特征值间在高维特征空间中具有多尺度的隐含关联特征。因此，在本申请的技术方案中，为了提高对于所述待测试土样的剪切模量测量的准确性，使用具有不同尺度的卷积核的卷积神经网络模型来进行特征挖掘，来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序上的关联性分布信息的多尺度隐含关联特征。具体地，首先，将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵。
具体地，在本申请实施例中，所述将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵，包括：使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在步骤S150中，将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵。这里，在将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵的同时，将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵。应可以理解，步骤S140和步骤S150中所述第一卷积神经网络模型和所述第二卷积神经网络模型是并行结构。
具体地，在本申请实施例中，所述将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵，包括：使用所述第二卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第二卷积神经网络模型的最后一层的输出为所述第二尺度关联参数特征矩阵，所述第二卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在步骤S160中，融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵。也就是，融合所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时序维度上的关联性特征分布表示的多尺度隐含关联特征信息，并以此作为解码特征矩阵。
具体地，在本申请实施例中，所述融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵，包括：以如下公式来融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；其中，所述公式为：
Figure 821878DEST_PATH_IMAGE045
其中，
Figure 544983DEST_PATH_IMAGE046
为所述解码特征矩阵，
Figure 624935DEST_PATH_IMAGE047
为所述第一尺度关联参数特征矩阵，
Figure 775293DEST_PATH_IMAGE048
为所述第二尺度关联参数特征矩阵，“
Figure 839064DEST_PATH_IMAGE049
”表示所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵相对应位置处的元素相加，
Figure 365860DEST_PATH_IMAGE050
和
Figure 34739DEST_PATH_IMAGE051
为用于控制所述解码特征矩阵中所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵之间的平衡的加权参数。
在步骤S170中，对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵。特别地，在本申请的技术方案中，考虑到所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵分别表达不同尺度下的所述轴向应力值和所述轴向应变时程值的时序关联参数在高维特征空间下的关联模式表达特征，所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵在高维特征空间内的特征分布既有同相重合又有异相差异，这使得当通过点加的方式融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵时，由于同相增强和异相削弱，所述解码特征矩阵的整体特征分布的离散程度高，在通过解码器进行解码时导致解码回归的归纳发散，从而影响通过解码器进行解码时的收敛性，影响解码值的准确性。因此，在通过解码器进行解码之前，首先对所述解码特征矩阵进行面向回归有界闭域的分布转移优化。
具体地，在本申请实施例中，所述对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵，包括：以如下公式对所述解码特征矩阵进行特征分布校正以得到所述校正后解码特征矩阵；其中，所述公式为：
Figure 636227DEST_PATH_IMAGE052
其中
Figure 921715DEST_PATH_IMAGE053
是所述解码特征矩阵的第
Figure 721043DEST_PATH_IMAGE054
位置的特征值，
Figure 41166DEST_PATH_IMAGE055
和
Figure 533327DEST_PATH_IMAGE056
是所述解码特征矩阵的宽度和高度，且
Figure 571691DEST_PATH_IMAGE057
表示以2为底的对数函数值，
Figure 440289DEST_PATH_IMAGE058
是所述校正后解码特征矩阵的第
Figure 614919DEST_PATH_IMAGE054
位置的特征值。
也就是，所述面向回归有界闭域的分布转移优化针对所述解码特征矩阵所代表的高维特征分布在转移到解码回归问题的目标域时可能存在的归纳发散，通过基于条件化解码回归的结构化信息约束来向着特征集合的有界闭域进行特征分布的收敛，以使得所述解码特征矩阵的特征分布转移到目标域下的具有稳定的可结构化边界的范围内，提升了回归解的归纳迭代的稳定性，即，提高了所述解码特征矩阵通过解码器进行解码时的收敛性，改进了解码值的准确性。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
在步骤S180中，将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。也就是，通过所述解码器对所述解码特征矩阵进行解码回归以得到用于表示待测试土样的剪切模量的解码值。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
具体地，在本申请实施例中，所述将所述解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量，包括：使用解码器以如下公式对所述解码特征矩阵进行解码回归以获得所述解码值；其中，所述公式为：
Figure 77649DEST_PATH_IMAGE059
，其中
Figure 72150DEST_PATH_IMAGE060
是所述解码特征矩阵，
Figure 478860DEST_PATH_IMAGE061
是所述解码值，
Figure 39154DEST_PATH_IMAGE062
是权重矩阵，
Figure 607539DEST_PATH_IMAGE063
表示矩阵乘。
综上，基于本申请实施例的高精度土样参数测量方法被阐明，其采用基于深度学习的人工智能算法来提取出所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值在时间维度上的多尺度关联性特征分布表示，并以此来进行解码回归得到所述待测试土样的剪切模量。这样，能够提高土样剪切模量测量的精准度，以智能准确地对于土样的剪切模量进行测量。
实施例2:图3为根据本申请实施例的高精度土样参数测量系统的框图。如图3所示，根据本申请实施例的高精度土样参数测量系统100，包括：数据获取模块110，用于获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；数据结构化模块120，用于将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；关联模块130，用于计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；第一尺度编码模块140，用于将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；第二尺度编码模块150，用于将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；融合模块160，用于融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；特征分布校正模块170，用于对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和，解码模块180，用于将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
在一个示例中，在上述高精度土样参数测量系统100中，所述第一尺度编码模块140，进一步用于：使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在一个示例中，在上述高精度土样参数测量系统100中，所述第二尺度编码模块150，进一步用于：使用所述第二卷积神经网络模型的各层在层的正向传递中分别进行：对输入数据进行卷积处理以得到卷积特征图；对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及，对所述池化特征图进行非线性激活以得到激活特征图；其中，所述第二卷积神经网络模型的最后一层的输出为所述第二尺度关联参数特征矩阵，所述第二卷积神经网络模型的第一层的输入为所述关联参数矩阵。
在一个示例中，在上述高精度土样参数测量系统100中，所述融合模块160，进一步用于：以如下公式来融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；其中，所述公式为：
Figure 886074DEST_PATH_IMAGE064
其中，
Figure 565317DEST_PATH_IMAGE065
为所述解码特征矩阵，
Figure 245697DEST_PATH_IMAGE066
为所述第一尺度关联参数特征矩阵，
Figure 453824DEST_PATH_IMAGE067
为所述第二尺度关联参数特征矩阵，“
Figure 222585DEST_PATH_IMAGE068
”表示所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵相对应位置处的元素相加，
Figure 908781DEST_PATH_IMAGE069
和
Figure 709247DEST_PATH_IMAGE070
为用于控制所述解码特征矩阵中所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵之间的平衡的加权参数。
在一个示例中，在上述高精度土样参数测量系统100中，所述特征分布校正模块170，进一步用于：以如下公式对所述解码特征矩阵进行特征分布校正以得到所述校正后解码特征矩阵；其中，所述公式为：
Figure 619434DEST_PATH_IMAGE071
其中
Figure 606981DEST_PATH_IMAGE072
是所述解码特征矩阵的第
Figure 362448DEST_PATH_IMAGE073
位置的特征值，
Figure 955103DEST_PATH_IMAGE074
和
Figure 567350DEST_PATH_IMAGE075
是所述解码特征矩阵的宽度和高度，且
Figure 511035DEST_PATH_IMAGE076
表示以2为底的对数函数值，
Figure 335772DEST_PATH_IMAGE077
是所述校正后解码特征矩阵的第
Figure 317022DEST_PATH_IMAGE073
位置的特征值。
在一个示例中，在上述高精度土样参数测量系统100中，所述解码模块180，进一步用于：使用解码器以如下公式对所述解码特征矩阵进行解码回归以获得所述解码值；其中，所述公式为：
Figure 631328DEST_PATH_IMAGE078
，其中
Figure 62310DEST_PATH_IMAGE079
是所述解码特征矩阵，
Figure 894000DEST_PATH_IMAGE080
是所述解码值，
Figure 257985DEST_PATH_IMAGE081
是权重矩阵，
Figure 149717DEST_PATH_IMAGE082
表示矩阵乘。
这里，本领域技术人员可以理解，上述高精度土样参数测量系统 100中的各个单元和模块的具体功能和操作已经在上面参考图1到图2的高精度土样参数测量方法的描述中得到了详细介绍，并因此，将省略其重复描述。
如上所述，根据本申请实施例的高精度土样参数测量系统 100可以实现在各种终端设备中，例如用于高精度土样参数测量的服务器等。在一个示例中，根据本申请实施例的高精度土样参数测量系统 100可以作为一个软件模块和/或硬件模块而集成到终端设备中。例如，该高精度土样参数测量系统100可以是该终端设备的操作系统中的一个软件模块，或者可以是针对于该终端设备所开发的一个应用程序；当然，该高精度土样参数测量系统 100同样可以是该终端设备的众多硬件模块之一。
替换地，在另一示例中，该高精度土样参数测量系统 100与该终端设备也可以是分立的设备，并且该高精度土样参数测量系统 100可以通过有线和/或无线网络连接到该终端设备，并且按照约定的数据格式来传输交互信息。
实施例3:对于土层反应分析中，
Figure 395891DEST_PATH_IMAGE083
，土的动剪切模量和阻尼是重要的计算参数，为获取较为准确的数值需要根据相关理论对离散数据进行计算。
图4为根据本申请实施例的应力应变滞回曲线图。在实际工程试验中，剪切模量的基本实现步骤为，首先，根据各级荷载下轴向应力值和轴向应变时程值的数据，绘出应力应变滞回曲线（如图4所示）。然后，以如下公式由滞回曲线的应力应变值求得土样的轴向模量E，以及，土样在该级荷载下剪切模量
Figure 765692DEST_PATH_IMAGE084
，
Figure 249763DEST_PATH_IMAGE085
，
Figure 303608DEST_PATH_IMAGE086
，式中，
Figure 240340DEST_PATH_IMAGE087
表示同一次循环荷载中最大轴向动应力；
Figure 944991DEST_PATH_IMAGE088
表示同一次循环荷载中最大轴向动应变；
Figure 221252DEST_PATH_IMAGE089
为泊松比，对于饱和试样取值为0.5；
Figure 48262DEST_PATH_IMAGE090
为剪应变幅值；
Figure 409973DEST_PATH_IMAGE091
为轴向应变幅值；
Figure 387157DEST_PATH_IMAGE092
为轴向模量。
一般而言，阻尼
Figure 845820DEST_PATH_IMAGE093
的计算公式为：
Figure 250256DEST_PATH_IMAGE094
，式中，△W为一个循环消耗的能量，用应力应变滞回曲线所围的面积表示，W为三角形的面积即弹性应变能。
下面介绍△W应力应变滞回曲线所围面积的计算方式。由于滞回曲线不光滑性、不确定性，如何计算一个循环消耗的能量△W成为动阻尼比计算的一个难点。在现有技术中，将滞回曲线简化成一个个三角形或者梯形，或者将点导入到cad里面描图计算面积，前者计算结果较为粗糙，后者较为繁琐。本申请的改进思路为：利用有限元计算的思路，将试验结果划分为50个单元格，将应力加一个数值后均变为大于0，然后计算。
而在传统计算方法中，通过公式
Figure 364843DEST_PATH_IMAGE095
可以得出每一级应变下的剪切模量、阻尼，将其所求得各级荷载下不同幅值的模量，采用Hardin-Drnevich模型表示，公式为：
Figure 210964DEST_PATH_IMAGE096
。 当
Figure 258554DEST_PATH_IMAGE097
时，可以得到系数A等于最大轴向模量的倒数，即1/Emax，将最大轴向模量值带入公式
Figure 99471DEST_PATH_IMAGE098
中，即可得到最大剪切模量值。
但对比现场波速结果、Hardin公式推导结果和动三轴反算结果（图5为根据本申请实施例的最大动剪切模量对比图）可以看出，Hardin公式推导结果、场波速结果规律性较好，动三轴反算结果较为离散。在数值上，Hardin公式推导结果、场波速结果随深度增大有明显的增大趋势，且两者大小区别不大，而动三轴结果在深度较小时与前两者区别不是很大，随着深度增加，差距越来越大。其主要原因在于在较小的应变范围内（应变小于10-5），Hardin公式能够较好地模拟黏性土的剪切模量随有效围压及孔隙比的变化规律。
《土工试验方法标准》（GBT 50123-2019）在振动三轴试验一节中提出，最大动弹性模量需要在微小应变下测得。同时李广信等认为一般的动三轴仪在动应变很小时测量精度很差，不适用于测定Gmax值。本申请采用地震工程现场取得的现场剪切波速值及室内测得的含水率、密度等参数对室内动三轴结果进行现场波速测试修正或者Hardin公式结果修正。
根据室内动三轴结果，以1/G为纵坐标，
Figure 498092DEST_PATH_IMAGE099
为横坐标，将试验结果标记后，近似用一条直线来表示它们的关系，即可得到B值。将B值及现场波速测试取得的A’和 Hardin公式推导取得的A”分别带入公式
Figure 817077DEST_PATH_IMAGE100
，即可获得不同结果下修正的剪切模量比，图6为根据本申请实施例的剪切模量比结果对比图。
根据Hardin-Drnevich模型，阻尼比D与剪切模量比的关系可写为：
Figure 188016DEST_PATH_IMAGE101
根据室内动三轴实测阻尼比λd，计算试样应变较大时的最大阻尼比λdmax，本申请取剪应变为5%时结果为最大阻尼比λdmax。图7为根据本申请实施例的动阻尼结果对比图。
以上结合具体实施例描述了本申请的基本原理，但是，需要指出的是，在本申请中提及的优点、优势、效果等仅是示例而非限制，不能认为这些优点、优势、效果等是本申请的各个实施例必须具备的。另外，上述公开的具体细节仅是为了示例的作用和便于理解的作用，而非限制，上述细节并不限制本申请为必须采用上述具体的细节来实现。
本申请中涉及的器件、装置、设备、系统的方框图仅作为例示性的例子并且不意图要求或暗示必须按照方框图示出的方式进行连接、布置、配置。如本领域技术人员将认识到的，可以按任意方式连接、布置、配置这些器件、装置、设备、系统。诸如“包括”、“包含”、“具有”等等的词语是开放性词汇，指“包括但不限于”，且可与其互换使用。这里所使用的词汇“或”和“和”指词汇“和/或”，且可与其互换使用，除非上下文明确指示不是如此。这里所使用的词汇“诸如”指词组“诸如但不限于”，且可与其互换使用。
还需要指出的是，在本申请的装置、设备和方法中，各部件或各步骤是可以分解和/或重新组合的。这些分解和/或重新组合应视为本申请的等效方案。
提供所公开的方面的以上描述以使本领域的任何技术人员能够做出或者使用本申请。对这些方面的各种修改对于本领域技术人员而言是非常显而易见的，并且在此定义的一般原理可以应用于其他方面而不脱离本申请的范围。因此，本申请不意图被限制到在此示出的方面，而是按照与在此公开的原理和新颖的特征一致的最宽范围。
为了例示和描述的目的已经给出了以上描述。此外，此描述不意图将本申请的实施例限制到在此公开的形式。尽管以上已经讨论了多个示例方面和实施例，但是本领域技术人员将认识到其某些变型、修改、改变、添加和子组合。

Claims (10)

1.一种高精度土样参数测量方法，其特征在于，包括：
获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；
将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；
计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；
将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；
将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；
融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；
对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和
将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
2.根据权利要求1所述的高精度土样参数测量方法，其特征在于，所述计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵，包括：计算所述第一输入向量的转置向量与所述第二输入向量之间的向量乘积以得到所述关联参数矩阵。
3.根据权利要求2所述的高精度土样参数测量方法，其特征在于，所述将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵，包括：
使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：
对输入数据进行卷积处理以得到卷积特征图；
对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及
对所述池化特征图进行非线性激活以得到激活特征图；
其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
4.根据权利要求3所述的高精度土样参数测量方法，其特征在于，所述将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵，包括：
使用所述第二卷积神经网络模型的各层在层的正向传递中分别进行：
对输入数据进行卷积处理以得到卷积特征图；
对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及
对所述池化特征图进行非线性激活以得到激活特征图；
其中，所述第二卷积神经网络模型的最后一层的输出为所述第二尺度关联参数特征矩阵，所述第二卷积神经网络模型的第一层的输入为所述关联参数矩阵。
5.根据权利要求4所述的高精度土样参数测量方法，其特征在于，所述融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵，包括：
以如下公式来融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；
其中，所述公式为：
Figure FDA0004081246790000021
其中，Md为所述解码特征矩阵，Ma为所述第一尺度关联参数特征矩阵，Mb为所述第二尺度关联参数特征矩阵，
Figure FDA0004081246790000022
表示所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵相对应位置处的元素相加，α和β为用于控制所述解码特征矩阵中所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵之间的平衡的加权参数。
6.根据权利要求5所述的高精度土样参数测量方法，其特征在于，所述对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵，包括：
以如下公式对所述解码特征矩阵进行特征分布校正以得到所述校正后解码特征矩阵；
其中，所述公式为：
Figure FDA0004081246790000023
其中mi,j是所述解码特征矩阵的第(i,j)位置的特征值，W和H是所述解码特征矩阵的宽度和高度，且log表示以2为底的对数函数值，mi,j'是所述校正后解码特征矩阵的第(i,j)位置的特征值。
7.根据权利要求6所述的高精度土样参数测量方法，其特征在于，所述将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量，包括：
使用解码器以如下公式对所述校正后解码特征矩阵进行解码回归以获得所述解码值；其中，所述公式为：
Figure FDA0004081246790000031
其中X是所述校正后解码特征矩阵，Y是所述解码值，W是权重矩阵，
Figure FDA0004081246790000032
表示矩阵乘。
8.一种高精度土样参数测量系统，其特征在于，包括：
数据获取模块，用于获取待测试土样在各级荷载下的轴向应力值和轴向应变时程值；
数据结构化模块，用于将所述待测试土样在各级荷载下的轴向应力值和轴向应变时程值分别按照时间维度排列为第一输入向量和第二输入向量；
关联模块，用于计算所述第一输入向量和所述第二输入向量之间的关联参数矩阵；
第一尺度编码模块，用于将所述关联参数矩阵通过使用具有第一尺度的卷积核的第一卷积神经网络模型以得到第一尺度关联参数特征矩阵；
第二尺度编码模块，用于将所述关联参数矩阵通过使用具有第二尺度的卷积核的第二卷积神经网络模型以得到第二尺度关联参数特征矩阵；
融合模块，用于融合所述第一尺度关联参数特征矩阵和所述第二尺度关联参数特征矩阵以得到解码特征矩阵；
特征分布校正模块，用于对所述解码特征矩阵进行特征分布校正以得到校正后解码特征矩阵；和
解码模块，用于将所述校正后解码特征矩阵通过解码器以得到解码值，所述解码值为待测试土样的剪切模量。
9.根据权利要求8所述的高精度土样参数测量系统，其特征在于，所述关联模块，用于计算所述第一输入向量的转置向量与所述第二输入向量之间的向量乘积以得到所述关联参数矩阵。
10.根据权利要求9所述的高精度土样参数测量系统，其特征在于，所述第一尺度编码模块，进一步用于：
使用所述第一卷积神经网络模型的各层在层的正向传递中分别进行：
对输入数据进行卷积处理以得到卷积特征图；
对所述卷积特征图进行基于局部特征矩阵的均值池化以得到池化特征图；以及
对所述池化特征图进行非线性激活以得到激活特征图；
其中，所述第一卷积神经网络模型的最后一层的输出为所述第一尺度关联参数特征矩阵，所述第一卷积神经网络模型的第一层的输入为所述关联参数矩阵。
CN202310005060.1A 2023-01-04 2023-01-04 高精度土样参数测量方法及其系统 Active CN115687854B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310005060.1A CN115687854B (zh) 2023-01-04 2023-01-04 高精度土样参数测量方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310005060.1A CN115687854B (zh) 2023-01-04 2023-01-04 高精度土样参数测量方法及其系统

Publications (2)

Publication Number Publication Date
CN115687854A CN115687854A (zh) 2023-02-03
CN115687854B true CN115687854B (zh) 2023-03-21

Family

ID=85057642

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310005060.1A Active CN115687854B (zh) 2023-01-04 2023-01-04 高精度土样参数测量方法及其系统

Country Status (1)

Country Link
CN (1) CN115687854B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116340703B (zh) * 2023-04-04 2023-08-29 山东省地震工程研究院 基于数据的土样参数测量方法及其系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112255095A (zh) * 2020-09-25 2021-01-22 汕头大学 一种土体应力应变关系确定方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016064476A1 (en) * 2014-10-24 2016-04-28 Halliburton Energy Services, Inc. Shale fracturing characterization and optimization using three-dimensional fracture modeling and neural network
JP2020030084A (ja) * 2018-08-22 2020-02-27 敏寛 松倉 応力ひずみ曲線作成装置および応力ひずみ曲線作成方法
KR102023701B1 (ko) * 2018-11-26 2019-11-04 한국건설기술연구원 이종 광섬유 센서의 측정결과 융합을 이용한 구조물의 변형률 분포 측정방법 및 장치
CN111666671A (zh) * 2020-05-29 2020-09-15 中南大学 一种围岩体蠕变参数的实时反演方法
CN114254416B (zh) * 2020-09-25 2024-06-25 汕头大学 基于长短期记忆深度学习的土体应力应变关系确定方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112255095A (zh) * 2020-09-25 2021-01-22 汕头大学 一种土体应力应变关系确定方法

Also Published As

Publication number Publication date
CN115687854A (zh) 2023-02-03

Similar Documents

Publication Publication Date Title
US10439594B2 (en) Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation
Dupuis Exceedances over high thresholds: A guide to threshold selection
CN108875178B (zh) 用于减小结构模态识别不确定性的传感器布置方法
CN115687854B (zh) 高精度土样参数测量方法及其系统
CN108563837A (zh) 一种冲积河流水沙模型的模型参数实时校正方法和系统
Achatz et al. A two-layer model with empirical linear corrections and reduced order for studies of internal climate variability
CN112329349B (zh) 一种边坡可靠度评估方法
CN107862170A (zh) 一种基于动态缩聚的有限元模型修正方法
CN104102836A (zh) 一种电力系统快速抗差状态估计方法
Bujurke et al. Numerical solution of stiff systems from nonlinear dynamics using single-term Haar wavelet series
CN112347537B (zh) 工程结构数值模型的校准方法、装置、电子设备及介质
CN113176022A (zh) 一种分段式神经网络压力传感器压力检测方法及系统
CN111428419A (zh) 悬浮泥沙浓度预测方法、装置、计算机设备和存储介质
CN106599610A (zh) 预测长链非编码rna和蛋白质联系的方法及系统
CN106909738A (zh) 一种模型参数辨识方法
CN111859303B (zh) 基于动态贝叶斯平均的土壤湿度融合方法及系统
Akritas et al. 7 Statistical analysis of censored environmental data
CN102830430B (zh) 一种层位速度建模方法
Hui et al. Application of grey system theory to forecast the growth of larch
CN103983332A (zh) 一种基于hgsa-bp算法的传感器误差补偿方法
CN109187898B (zh) 水产养殖环境中水质氨氮含量的软测量方法及装置
Arab et al. Enhancing weighted uniform simulation for structural reliability analysis
CN116340384A (zh) 基于规则演化的核递归最大相关熵时间序列在线预测方法
CN116108745A (zh) 水环境模型多参数率定方法、终端设备及存储介质
Bishop et al. Optimization of the fixed global observing network in a simple model

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