[發明專利]一種利用凸函數優化的微震事件檢測方法和系統

在審

• 申請號: 202011068043.5
• 申請日: 2020-10-08
• 公開（公告）號: CN112180439A
• 公開（公告）日: 2021-01-05
• 發明（設計）人: 翟明岳
• 申請（專利權）人: 廣東石油化工學院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種利用凸函數優化的微震事件檢測方法，其特征在于，包括：

步驟101獲取按時間順序采集的信號序列S；

步驟102通過實際調查，獲得微震事件類型矢量p，具體為：

其中：

N_i：微震事件類型總數

i＝1,2,···,N_i：微震事件類型序號

p_i：第i個類型微震事件的均方差，且按升序排列，即：

p₁p₂···

步驟103求取候選微震事件矢量e和候選狀態矩陣α^*，具體為：

其中：

如果ΔP_j≥p₁，則所述候選微震事件矢量e的第m個元素e_m＝j

所述候選狀態矩陣α^*的第o行第m列元素

m＝1,2,···,N_E：候選微震事件序號

N_E：候選微震事件總數

第e_m點RMS值

第j點RMS值

偏移值

對*上取整

*：表示任一自變量

N：所述信號序列S的長度

s_l：所述信號序列S的第l個元素

l：求和參量

如果所述元素序號l≤0，則所述元素序號置為1；

如果所述元素序號l≥N，則所述元素序號置為N；

m₀：所述信號序列S的均值

o＝1,2,···,N_i：微震事件類型序號；

步驟104求取狀態矩陣α，具體為：

其中，

所述狀態矩陣α的第o行第u列元素

中間參數

第e_u點RMS值

類型o微震事件的置信度

p_o：類型o微震事件RMS值

p_q：類型q微震事件RMS值

q＝1,2,···,N_i：微震事件類型序號；

步驟105判斷微震事件，具體為：對于狀態矩陣α中第i行第j列元素α_ij，如果α_ij＝1，則表明在數據點j處存在類型i的微震事件；如果α_ij＝0，則表明在數據點j處，不存在類型i的微震事件。

2.一種利用凸函數優化的微震事件檢測系統，其特征在于，包括：

模塊201獲取按時間順序采集的信號序列S；

模塊202通過實際調查，獲得微震事件類型矢量p，具體為：

其中：

N_i：微震事件類型總數

i＝1,2,···,N_i：微震事件類型序號

p_i：第i個類型微震事件的均方差，且按升序排列，即：

模塊203求取候選微震事件矢量e和候選狀態矩陣α^*，具體為：

其中：

如果ΔP_j≥p₁，則所述候選微震事件矢量e的第m個元素e_m＝j

所述候選狀態矩陣α^*的第o行第m列元素

m＝1,2,···,N_E：候選微震事件序號

N_E：候選微震事件總數

第e_m點RMS值

第j點RMS值

偏移值

對*上取整

*：表示任一自變量

N：所述信號序列S的長度

s_l：所述信號序列S的第l個元素

l：求和參量

如果所述元素序號l≤0，則所述元素序號置為1；

如果所述元素序號l≥N，則所述元素序號置為N；

m₀：所述信號序列S的均值

o＝1,2,···,N_i：微震事件類型序號；

模塊204求取狀態矩陣α，具體為：

其中，

所述狀態矩陣α的第o行第u列元素

中間參數

第e_u點RMS值

類型o微震事件的置信度

p_o：類型o微震事件RMS值

p_q：類型q微震事件RMS值

q＝1,2,···,N_i：微震事件類型序號；

模塊205判斷微震事件，具體為：對于狀態矩陣α中第i行第j列元素α_ij，如果α_ij＝1，則表明在數據點j處存在類型i的微震事件；如果α_ij＝0，則表明在數據點j處，不存在類型i的微震事件。