[發明專利]一種震電聯合獲取裂縫信息的定位方法

在審

• 申請號: 201811568064.6
• 申請日: 2018-12-21
• 公開（公告）號: CN109856677A
• 公開（公告）日: 2019-06-07
• 發明（設計）人: 王俊秋;馬立偉;田入運;孫佳;賓康成;羅世豪;許皓
• 申請（專利權）人: 吉林大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 沈陽銘揚聯創知識產權代理事務所(普通合伙) 21241
• 代理人: 屈芳

1.一種震電聯合獲取裂縫信息的定位方法，其特征在于，該方法包括：

利用微地震方法拾取的震源點陣列進行震源定位；

利用井地電法的視差分電阻率差值的變化對裂縫進行成像；

利用微地震方法拾取的震源點陣列和井地電法測量獲得的裂縫成像結果進行融合得到異常信息。

2.按照權利要求1所述的方法,其特征在于，所述利用微地震方法拾取的震源點陣列進行震源定位包括：

S1對微地震信號數據預處理；

S2建立速度模型，獲取界面的深度和各層的初始速度，求取到理論到時，通過正演所求取的到時與理論到時進行逼近，直到滿足精度要求，確定實際的各層的初始速度；

S3基于所述微地震信號，采用STA/LTA方法對微地震事件進行拾取，對微地震事件進行拾取，確定有效的地震時間；

S4基于微地震事件、實際的各層的初始速度和地質參數，利用震幅疊加法進行震源定位。

3.按照權利要求2所述的方法,其特征在于，所述步驟S2具體包括：

根據壓裂地層的實際資料建立速度模型，獲取界面的深度和各層的初始速度；

根據射孔位置和檢波器的位置利用射線追蹤方法求取每個檢波器位置和理論到時；

利用正演求取的到時差與實際的到時差進行逼近，求取偏離方差和，若不滿足探測精度要求，對速度值進行修正，修正后重新對實際到時差進行逼近，直到滿足精度為止。

4.按照權利要求2所述的方法，其特征在于，所述步驟S4具體包括：

步驟S41：獲取三維目標區域，根據震源定位精度將所述三維目標區域劃分成N個體元；

步驟S42：對體元i到檢波器j的距離進行計算，并根據地質參數和地震波的傳播速度計算得到體元i到檢波器j的走時δ_ij，i從1到N，j從1取到M；

步驟S43：將檢波器j的數據按照對應體元i的走時沿逆時間軸方向進行移動，得到新的移動數據X_j(t-δ_ij)，i從1取到N，j從1取到M；

步驟S44：對于體元i，將所有移動后的采集數據進行疊加，得到數據Y_i(t)；

步驟S45：將i從1取到N，得到所有體元的疊加數據，所述疊加數據為所有體元到所有檢波器的振幅疊加數據；

步驟S46：得到震幅疊加數據后，需對其進行掃描，獲取四維數組(x,y,z,δ)，得到數組極大值，為能量最高值點，確定為微地震震源點。

5.按照權利要求4所述的方法，其特征在于，所述步驟S46中利用洪水填充發對極大值進行確認，判斷真正震源點位置。

6.按照權利要求1所述的方法,其特征在于，利用井地電法的視差分電阻率差值的變化對裂縫進行成像包括：

S6對電法數據的濾波處理；

S7對電法數據進行數據地形校正，計算出純地形和水平地形下均勻半空間的視差分電阻率；

S8基于所述視差分電阻率計算校正因子，基于所述校正因子分別對電阻率和電位梯度進行校正，計算得到電位梯度；

S9基于所述電位梯度計算出每個采集點的視差分電阻率值，基于所述視差分電阻率值獲取對應的裂縫成像信息。

7.按照權利要求6所述的方法,其特征在于，所述步驟S6對電法數據的濾波處理包括：

S61采用低通濾波方法濾除諧波干擾；

S62采用向后差分的方法對對數據中的過沖或尖峰進行去除；

S63對數據的每個周期求取平均值，利用平均值構成一組一維數據并根據數據的特點進行最小二乘多項式擬合，從而獲取偏置或基線漂移擬合曲線，然后將原始數據減去漂移擬合曲線獲得校正后的電法數據；

S64采用多個周期疊加的方法，去除電法數據中的隨機噪聲干擾。

8.按照權利要求6所述的方法,其特征在于，電法數據通過布設為內圈、中圈、外圈的儀器采集得到，裂縫描繪有以下幾種情況：

1)當某一方向上內圈、中圈、外圈的視差分電阻率變化量大時，說明此方向上裂縫的長度長，涉及到外圈；

2)當某一方向上內圈視差分電阻率變化大時，中圈視差分電阻率變化不大時，說明此方向上裂縫長度短，未涉及到中圈；若中圈視差分電阻率變化大，而外圈視差分電阻率變化不大時，說明裂縫涉及到中圈，但未涉及到外圈；

3)當某一方向上內圈、中圈視差分電阻率變化不大時，無論外圈視差分電阻率變化是否明顯，該方向均沒有裂縫產生。