[發(fā)明專利]一種三維VSP數(shù)據(jù)初至波拾取方法

有效

• 申請?zhí)枺?/td>: 200910241396.8
• 申請日: 2009-12-01
• 公開（公告）號: CN102081167A
• 公開（公告）日: 2011-06-01
• 發(fā)明（設(shè)計）人: 陳沅忠;李彥鵬
• 申請（專利權(quán)）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責(zé)任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京市中實友知識產(chǎn)權(quán)代理有限責(zé)任公司 11013
• 代理人: 劉天語

1.一種三維VSP數(shù)據(jù)初至波拾取方法，其特征在于：具體步驟包括：

1)用地震勘探方法取得三維三分量VSP數(shù)據(jù)集，進行觀測系統(tǒng)定義、去除單頻干擾、振幅恢復(fù)、靜校正處理；

2)抽取若干炮不同井源距的VSP數(shù)據(jù)，初步拾取其初至，形成每個觀測點深度的初至?xí)r間與井源距關(guān)系表；

3)對每一道數(shù)據(jù)，利用初至井源距關(guān)系表采用三次樣條插值出此道的先導(dǎo)初至t_p；

4)采用一個大于零的時窗Δt，截取(t_p-Δt)～(t_p+Δt)時間段的數(shù)據(jù)；

5)對截取的數(shù)據(jù)進行三分量合成模數(shù)道處理，公式為：

x_i、y_i、z_i表示三分量VSP數(shù)據(jù)第i個樣點振幅值，m_i表示模數(shù)道第i個樣點振幅值；

6)根據(jù)精度要求確定新的采樣率，采用三次樣條算法按新采樣率把模數(shù)道數(shù)據(jù)重采樣；

7)采用以下公式計算模數(shù)道數(shù)據(jù)每個樣點的相鄰振幅差的乘積值dm_i，dm_i＝(m_i-m_i-1)·(m_i+1-m_i)，m_i表示模數(shù)道第i個樣點振幅值，dm_i小于零的值為波峰，并找出其中的最大波峰值A(chǔ)；

8)用小于1的比例因子α，計算每一個波峰與A的比值，確定比值大于α的波峰，其中的第一個波峰是初至后的第一個波峰，這個波峰值計為B；

9)用小于1的比例因子β，計算B之前的數(shù)據(jù)與B的比值，確定比值小于β的距離B最近的一個樣點，對應(yīng)的時間為拾取初至；

10)重復(fù)(3)～(9)至全部三維VSP數(shù)據(jù)拾取；

11)在炮集內(nèi)利用相鄰深度初至?xí)r差輔助手動拾取初至；

12)對全部數(shù)據(jù)按井源距分組，在每個井源距內(nèi)計算相鄰深度初至?xí)r差dt_j，公式為：dt_j＝(t_j-t_j-1)/(d_j-d_j-1)，t_j表示第j道初至，d_j表示第j道深度值；對每一個深度的時差值序列按大小排序，找到中間值meddt，用允許時差范圍ΔT₂判斷，找出dt＞meddt+ΔT₂或dt＜meddt-ΔT₂的初至，返回到步驟11)重新拾取初至；

13)對全部數(shù)據(jù)按炮檢水平方位進行分組，在每個炮檢水平方位組內(nèi)把初至按井源距從小到大排序，對每一個深度初至計算后一炮初至與前一炮的初至差dt_ij，公式為：dt_ij＝(t_ij-t_(i-1)j)，t_ij表示第j深度道第i炮初至。找出dt_ij＜0的初至，返回到步驟11重新拾取初至；

14)對全部數(shù)據(jù)按接收點深度進行分組，在每個接收點深度內(nèi)有地理坐標(biāo)y、x及初至t三列數(shù)據(jù)，利用基于三角形的線性插值算法對(y，x，t)進行二元函數(shù)曲面擬合，形成規(guī)則網(wǎng)格數(shù)據(jù)(Y，X，T)，對(Y，X，T)進行多道二維中值濾波處理得到(Y，X，Tm)，計算dT＝|T-Tm|，用允許時差范圍ΔT₃判斷，找出dT＞ΔT₃的初至，返回到步驟11重新拾取初至。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于：步驟2)所述的井源距VSP數(shù)據(jù)包括最小井源距和最大井源距數(shù)據(jù)。

3.根據(jù)權(quán)利要求1所述的方法，其特征在于：步驟6)新的采樣率為0.01-0.1ms。