[发明专利]一种声辐射剪切波波速检测方法及系统

摘要

本申请公开了一种声辐射剪切波波速检测方法及系统，包括分别选择冲击位置、第一检测位和第二检测位，对冲击位置进行声辐射力冲击；沿深度方向分别对两组时间‑深度‑剪切波位移矩阵进行累加，得到两组由累加和组成的累加位移数组；分别对两组累加位移数组进行最大值搜索，将最大值作为剪切波波峰位置处的时间；根据所述第一检测位和所述第二检测位之间的距离及所述第一检测位和所述第二检测位的剪切波波峰位置时间，计算得出剪切波速度。本申请无需寻找位移矩阵中的焦点位置，而是先把位移矩阵沿深度方向累加，综合利用不同深度的剪切波波峰信息，找出在具有最多剪切波波峰通过的时间位置，求出最大值所在的横向位置，即可计算剪切波速度。

主权项

1.一种声辐射剪切波波速检测方法，其特征在于，包括：分别选择冲击位置、第一检测位和第二检测位，对所述位置进行声辐射力冲击；采集所述第一检测位和所述第二检测位的数据分别进行处理，得出两组时间‑深度‑剪切波位移矩阵；沿深度方向分别对两组所述时间‑深度‑剪切波位移矩阵进行累加，得到两组由累加和组成的累加位移数组；分别对两组所述累加位移数组进行最大值搜索，将最大值作为剪切波波峰位置处的时间；根据所述第一检测位和所述第二检测位之间的距离及所述第一检测位和所述第二检测位的剪切波波峰位置时间，计算得出剪切波速度；所述采集所述第一检测位和所述第二检测位的数据分别进行处理，具体包括：利用互相关分析分别对在所述第一检测位和所述第二检测位采集的数据进行处理，其中，利用互相关分析分别对在所述第一检测位和所述第二检测位采集的数据进行处理，包括对数据矩阵相邻列进行互相关位移估计计算，设步长为Step，窗长为WinLength，得出互相关位移矩阵大小为K*(N‑1),其中K＝(M‑WinLength)/Step+1，M为在第一检测位采集数据的纵向长度，N为采集次数，此矩阵即为时间‑深度‑剪切波位移矩阵。