[發明專利]一種聲輻射剪切波波速檢測方法及系統

有效

• 申請號: 201510514501.6
• 申請日: 2015-08-20
• 公開（公告）號: CN105212961B
• 公開（公告）日: 2018-08-31
• 發明（設計）人: 翁嘉淳;羅元亮;涂希;周偉峰;李強
• 申請（專利權）人: 深圳市紅源資產管理有限公司
• 主分類號: A61B8/00
• 分類號: A61B8/00
• 代理公司: 深圳鼎合誠知識產權代理有限公司 44281
• 代理人: 任葵;胥強

1.一種聲輻射剪切波波速檢測方法，其特征在于，包括：

分別選擇沖擊位置、第一檢測位和第二檢測位，對所述位置進行聲輻射力沖擊；

采集所述第一檢測位和所述第二檢測位的數據分別進行處理，得出兩組時間-深度-剪切波位移矩陣；

沿深度方向分別對兩組所述時間-深度-剪切波位移矩陣進行累加，得到兩組由累加和組成的累加位移數組；

分別對兩組所述累加位移數組進行最大值搜索，將最大值作為剪切波波峰位置處的時間；

根據所述第一檢測位和所述第二檢測位之間的距離及所述第一檢測位和所述第二檢測位的剪切波波峰位置時間，計算得出剪切波速度；

所述采集所述第一檢測位和所述第二檢測位的數據分別進行處理，具體包括：

利用互相關分析分別對在所述第一檢測位和所述第二檢測位采集的數據進行處理，其中，利用互相關分析分別對在所述第一檢測位和所述第二檢測位采集的數據進行處理，包括對數據矩陣相鄰列進行互相關位移估計計算，設步長為Step，窗長為WinLength，得出互相關位移矩陣大小為K*(N-1),其中K＝(M-WinLength)/Step+1，M為在第一檢測位采集數據的縱向長度，N為采集次數，此矩陣即為時間-深度-剪切波位移矩陣。

2.如權利要求1所述的聲輻射剪切波波速檢測方法，其特征在于，所述沿深度方向分別對兩組所述時間-深度-剪切波位移矩陣進行累加，包括：

選取局部數據進行累加。

3.如權利要求1所述的聲輻射剪切波波速檢測方法，其特征在于，所述沿深度方向分別對兩組所述時間-深度-剪切波位移矩陣進行累加前，還包括：

對所述時間-深度-剪切波位移矩陣進行橫向插值。

4.如權利要求1所述的聲輻射剪切波波速檢測方法，其特征在于，所述沿深度方向分別對兩組所述時間-深度-剪切波位移矩陣進行累加前，還包括：

對所述時間-深度-剪切波位移矩陣進行橫向濾波。

5.一種聲輻射剪切波波速檢測系統，其特征在于，包括：

選擇模塊，用于分別選擇沖擊位置、第一檢測位和第二檢測位，對所述沖擊位置進行聲輻射力沖擊；

采集模塊，用于采集所述第一檢測位和所述第二檢測位的數據分別進行處理，得出兩組時間-深度-剪切波位移矩陣；

累加模塊，用于沿深度方向分別對兩組所述時間-深度-剪切波位移矩陣進行累加，得到兩組由累加和組成的累加位移數組；

查找模塊，用于分別對兩組所述累加位移數組進行最大值搜索，將最大值作為剪切波波峰位置處的時間；

處理模塊，用于根據所述第一檢測位和所述第二檢測位之間的距離及所述第一檢測位和所述第二檢測位的剪切波波峰位置時間，計算得出剪切波速度；

所述采集模塊，還用于利用互相關分析分別對在所述第一檢測位和所述第二檢測位采集的數據進行處理，其中，利用互相關分析分別對在所述第一檢測位和所述第二檢測位采集的數據進行處理，包括對數據矩陣相鄰列進行互相關位移估計計算，設步長為Step，窗長為WinLength，得出互相關位移矩陣大小為K*(N-1),其中K＝(M-WinLength)/Step+1，此矩陣即為時間-深度-剪切波位移矩陣。

6.如權利要求5所述的聲輻射剪切波波速檢測系統，其特征在于，所述累加模塊，還用于選取局部數據進行累加。

7.如權利要求5所述的聲輻射剪切波波速檢測系統，其特征在于，所述累加模塊，還用于對所述時間-深度-剪切波位移矩陣進行橫向插值。

8.如權利要求5所述的聲輻射剪切波波速檢測系統，其特征在于，所述累加模塊，還用于對所述時間-深度-剪切波位移矩陣進行橫向濾波。