[發(fā)明專利]微地震震級-頻度關系中的b值計算方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201811044723.6
• 申請日: 2018-09-07
• 公開（公告）號: CN110888161A
• 公開（公告）日: 2020-03-17
• 發(fā)明（設計）人: 楊小慧;丁圣
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創(chuàng)畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種微地震震級-頻度關系中的b值計算方法，其特征在于，包括：

通過壓裂段的每個事件的破裂能量，計算每個事件對應的震級數(shù)據(jù)；

根據(jù)所述每個事件對應的震級數(shù)據(jù)，獲得累計震級分布數(shù)據(jù)；

根據(jù)所述累計震級分布數(shù)據(jù)，通過最大曲率法計算最小震級；

根據(jù)大于所述最小震級的震級，通過線性最小二乘法計算b值。

2.根據(jù)權利要求1所述的微地震震級-頻度關系中的b值計算方法，其中，所述每個事件對應的震級數(shù)據(jù)為：

其中，M₀為每個事件的破裂能量，M_w為每個事件對應的震級數(shù)據(jù)。

3.根據(jù)權利要求1所述的微地震震級-頻度關系中的b值計算方法，其中，通過線性最小二乘法計算b值包括：

確定震級的經(jīng)驗分布函數(shù)；

根據(jù)所述震級的經(jīng)驗分布函數(shù)，確定震級目標函數(shù)；

根據(jù)所述震級目標函數(shù)，計算b值。

4.根據(jù)權利要求3所述的微地震震級-頻度關系中的b值計算方法，其中，所述震級的經(jīng)驗分布函數(shù)為：

其中，M₍₁₎≤M₍₂₎≤…≤M_(K)≤M_(K+1)≤…M_(N)為按大小排序后的震級，S_N(M)為震級M的經(jīng)驗分布函數(shù)。

5.根據(jù)權利要求4所述的微地震震級-頻度關系中的b值計算方法，其中，所述震級目標函數(shù)為：

其中，Q₁為震級目標函數(shù)，M_(i)為排序后的第i個震級，S_N(M_(i))為震級M_(i)的經(jīng)驗分布函數(shù)，b為待求b值。

6.根據(jù)權利要求4所述的微地震震級-頻度關系中的b值計算方法，其中，所述b值為：

其中，b為待求b值，M_(i)為排序后的第i個震級，Z_i＝-ln[1-S_N(M_(i))]。

7.一種微地震震級-頻度關系中的b值計算系統(tǒng)，其特征在于，該系統(tǒng)包括：

存儲器，存儲有計算機可執(zhí)行指令；

處理器，所述處理器運行所述存儲器中的計算機可執(zhí)行指令，執(zhí)行以下步驟：

通過壓裂段的每個事件的破裂能量，計算每個事件對應的震級數(shù)據(jù)；

根據(jù)所述每個事件對應的震級數(shù)據(jù)，獲得累計震級分布數(shù)據(jù)；

根據(jù)所述累計震級分布數(shù)據(jù)，通過最大曲率法計算最小震級；

根據(jù)大于所述最小震級的震級，通過線性最小二乘法計算b值。

8.根據(jù)權利要求7所述的微地震震級-頻度關系中的b值計算系統(tǒng)，其中，所述每個事件對應的震級數(shù)據(jù)為：

其中，M₀為每個事件的破裂能量，M_w為每個事件對應的震級數(shù)據(jù)。

9.根據(jù)權利要求7所述的微地震震級-頻度關系中的b值計算系統(tǒng)，其中，通過線性最小二乘法計算b值包括：

確定震級的經(jīng)驗分布函數(shù)；

根據(jù)所述震級的經(jīng)驗分布函數(shù)，確定震級目標函數(shù)；

根據(jù)所述震級目標函數(shù)，計算b值。

10.根據(jù)權利要求9所述的微地震震級-頻度關系中的b值計算系統(tǒng)，其中，所述b值為：

其中，b為待求b值，M_(i)為排序后的第i個震級，Z_i＝-ln[1-S_N(M_(i))]。