[发明专利]一种自动切除地震波波至的方法

摘要

一种自动切除地震波波至的方法，适用于地震勘探资料处理领域，特点是计算第k层的折射波临界角αk‑1(0≦k≦n‑1)的余弦值；将折射波临界角αk‑1(0≦k≦n‑1)的余弦值代入第n(n>2)层的折射波时距曲线方程；对每一个折射层求解折射波时距曲线方程，得到偏移距X上所有折射波的波至时间，当X小于盲区半径时，折射波的波至时间取值为0；取其中的最大值tXmax，为偏移距X上的切除值；把所有偏移距对应的tXmax值连接成线后获得折射波自动切除线；取地层的层数n＝2，此时的地层模型为低降速带模型，采用浅层折射波法实现地震波波至的自动切除。本发明能准确剔除折射干扰能够保证取得较真实的分辨率。

主权项

一种自动切除地震波波至的方法，特点是通过以下步骤实现：1)根据由微测井或测井得到的n层水平层状介质的层速度和厚度，第k层的速度为vk‑1，地层厚度为hk‑1，用下式计算第k层的折射波临界角αk‑1的余弦值：cosαk‑1＝(1‑(vk‑1/vk)2)0.5  (1)式中，αk‑1是第k层的折射波临界角，vk‑1是第k层的速度，vk是第k+1层的速度，0≦k≦n‑1，n是地层的层数；2)将折射波临界角αk‑1的余弦值代入第n层的折射波时距曲线方程：t=Xvn-1+Σk=1n-12hk-1cosαk-1vk-1---(2)]]>式中，0≦k≦n‑1，n是地层的层数，n>2，X是偏移距，vn‑1是第n层地层的速度，hk‑1是第k层地层的厚度，αk‑1是第k层地层的折射波临界角，vk‑1是第k层的速度，t是折射波在偏移距X处的旅行时间；3)对每一个折射层求解折射波时距曲线方程，得到偏移距X上所有折射波的波至时间，当X小于盲区半径时，折射波的波至时间取值为0；取其中的最大值tXmax，作为偏移距X上的切除值；其中tXmax为同一偏移距X上所有折射波的波至时间中的最大值；4)把所有偏移距对应的tXmax值连接成线后获得折射波自动切除线；5)取地层的层数n＝2，此时的地层模型为低降速带模型，采用浅层折射波法实现地震波波至的自动切除。