一種基于垂直地震衰減規(guī)律的中值濾波方法以及系統(tǒng)的制作方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明關(guān)于地球物理勘探技術(shù)領(lǐng)域，特別是關(guān)于地球物理勘探中垂直地震數(shù)據(jù)處 理領(lǐng)技術(shù)，具體的講是一種基于垂直地震衰減規(guī)律的中值濾波方法及系統(tǒng)。
【背景技術(shù)】
[0002] 波場信息豐富是垂直地震剖面VSP數(shù)據(jù)的一個(gè)重要特點(diǎn)，在剖面中可以看到多種 波場相互交織的現(xiàn)象。通常，VSP數(shù)據(jù)又不具備地面地震較高的覆蓋次數(shù)，故而很難通過疊 加的方法強(qiáng)化或衰減某一種波場。因此，波場分離方法在VSP數(shù)據(jù)的處理中的應(yīng)用就顯示 尤為重要。
[0003] 現(xiàn)有技術(shù)中，VSP數(shù)據(jù)的處理方案中并沒有針對(duì)性的波場分離方法，其波場分離手 段都是借用地面地震數(shù)據(jù)處理的方法，如頻率一波數(shù)濾波（FK濾波）、均值濾波和中值濾波 等。其中，F(xiàn)K濾波容易產(chǎn)生假頻，在覆蓋次數(shù)低、質(zhì)量要求高的VSP數(shù)據(jù)波場中的應(yīng)用十 分有限；均值濾波和中值濾波方法都屬于時(shí)間域的多道速度濾波，中值濾波的抗噪能力優(yōu) 于均值濾波，在實(shí)際生產(chǎn)中應(yīng)用較多，并且派生出自適應(yīng)去異常、邊界變道數(shù)等中值濾波方 法。在VSP數(shù)據(jù)處理中，雖然這些波場分離方法并沒有考慮到VSP數(shù)據(jù)自身的特點(diǎn)，但在傳 統(tǒng)的處理流程中，波場分離步驟是在振幅補(bǔ)償之后進(jìn)行的，此時(shí)剖面振幅能量相對(duì)統(tǒng)一，是 可以利用中值方法進(jìn)行波場分離的。
[0004] 近年來，VSP數(shù)據(jù)在井控技術(shù)方面得到了快速發(fā)展，其振幅衰減參數(shù)、球面擴(kuò)散補(bǔ) 償參數(shù)（TAR)、地層吸收補(bǔ)償參數(shù)（Q)等在地面地震處理中大量應(yīng)用，起到了很好的井控約 束作用。然而，這些井控參數(shù)的提取都需要在VSP數(shù)據(jù)振幅相對(duì)保真前提下完成，要求在振 幅隨VSP數(shù)據(jù)記錄深度呈指數(shù)變化的剖面中進(jìn)行波場分離。通常，一個(gè)同相軸的不同位置 存在著非常大的振幅能量差異，這將極大的影響中值濾波方法的可靠性。
[0005] 因此，如何提出一種新的濾波方法，進(jìn)而解決用于井控的VSP數(shù)據(jù)數(shù)據(jù)處理中沒 有可靠波場分離是本領(lǐng)域亟待解決的技術(shù)難題。

【發(fā)明內(nèi)容】

[0006] 為了解決現(xiàn)有技術(shù)中用于井控的VSP數(shù)據(jù)數(shù)據(jù)處理中沒有可靠波場分離方法的 問題，本發(fā)明在自適應(yīng)去異常的、邊界變道數(shù)的中值濾波基礎(chǔ)上，提出一種基于垂直地震衰 減規(guī)律的中值濾波方法及系統(tǒng)，是一種精確的基于垂直地震衰減規(guī)律的中值濾波的方案， 由VSP數(shù)據(jù)獲得單分量數(shù)據(jù)，對(duì)單分量數(shù)據(jù)進(jìn)行零相位化處理，結(jié)合時(shí)距曲線確定TAR因 子，并沿時(shí)距曲線拉齊剖面，進(jìn)行定向加權(quán)中值濾波，極大提高了中值濾波的可靠性。
[0007] 本發(fā)明的目的之一是，提供一種基于垂直地震衰減規(guī)律的中值濾波方法，包括：采 集垂直地震剖面VSP數(shù)據(jù)；根據(jù)所述的VSP數(shù)據(jù)分選出單分量數(shù)據(jù)；對(duì)所述的單分量數(shù)據(jù) 進(jìn)行零相位化處理，得到零相位化處理后的剖面；根據(jù)所述零相位化處理后的剖面得到時(shí) 距曲線；根據(jù)所述零相位化處理后的剖面以及所述的時(shí)距曲線確定所述VSP數(shù)據(jù)的球面擴(kuò) 散補(bǔ)償TAR因子；將零相位化處理后的剖面沿所述時(shí)距曲線拉齊，得到拉齊后的剖面；根據(jù) 所述VSP數(shù)據(jù)的TAR因子對(duì)所述拉齊后的剖面進(jìn)行定向加權(quán)中值濾波。
[0008] 本發(fā)明的目的之一是，提供了一種基于垂直地震衰減規(guī)律的中值濾波的系統(tǒng)，包 括：剖面采集裝置，用于采集垂直地震剖面VSP數(shù)據(jù)；單分量數(shù)據(jù)分選裝置，用于根據(jù)所述 的VSP數(shù)據(jù)分選出單分量數(shù)據(jù)；零相位化處理裝置，用于對(duì)所述的單分量數(shù)據(jù)進(jìn)行零相位 化處理，得到零相位化處理后的剖面；時(shí)距曲線確定裝置，用于根據(jù)所述零相位化處理后的 剖面得到時(shí)距曲線；TAR因子確定裝置，用于根據(jù)所述零相位化處理后的剖面以及所述的 時(shí)距曲線確定所述VSP數(shù)據(jù)的球面擴(kuò)散補(bǔ)償TAR因子；拉齊裝置，用于將零相位化處理后 的剖面沿所述時(shí)距曲線拉齊，得到拉齊后的剖面；中值濾波裝置，用于根據(jù)所述VSP數(shù)據(jù)的 TAR因子對(duì)所述的拉齊后的剖面進(jìn)行定向加權(quán)中值濾波。
[0009] 本發(fā)明的有益效果在于，在自適應(yīng)去異常的、邊界變道數(shù)的中值濾波基礎(chǔ)上，提供 了一種基于垂直地震衰減規(guī)律的中值濾波方法及系統(tǒng)，是一種精確的基于垂直地震衰減規(guī) 律的中值濾波的方案，由VSP數(shù)據(jù)分選出單分量數(shù)據(jù)，對(duì)單分量數(shù)據(jù)進(jìn)行零相位化處理，結(jié) 合時(shí)距曲線確定TAR因子，并沿時(shí)距曲線拉齊剖面，進(jìn)行定向加權(quán)中值濾波，極大提高了中 值濾波的可靠性，為用于井控的VSP數(shù)據(jù)數(shù)據(jù)處理提供了可靠的波場分離方案。
[0010] 為讓本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉較佳實(shí)施例， 并配合所附圖式，作詳細(xì)說明如下。
【附圖說明】
[0011] 為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可 以根據(jù)這些附圖獲得其他的附圖。
[0012] 圖1為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波方法的流程 圖；
[0013] 圖2為圖1中的步驟S104的具體流程圖；
[0014] 圖3為圖1中的步驟S105的具體流程圖；
[0015] 圖4為圖1中的步驟S107的具體流程圖；
[0016] 圖5為圖4中的步驟S403的具體流程圖；
[0017] 圖6為圖4中的步驟S405的具體流程圖；
[0018] 圖7為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)的結(jié)構(gòu) 框圖；
[0019] 圖8為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)中的時(shí) 距曲線確定裝置的具體結(jié)構(gòu)框圖；
[0020] 圖9為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)中的TAR 因子確定裝置的具體結(jié)構(gòu)框圖；
[0021] 圖10為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)中的中 值濾波裝置的具體結(jié)構(gòu)框圖；
[0022] 圖11為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)中的加 權(quán)數(shù)據(jù)序列確定模塊的具體結(jié)構(gòu)框圖；
[0023]圖12為本發(fā)明實(shí)施例提供的一種基于垂直地震衰減規(guī)律的中值濾波系統(tǒng)中的去 異常值處理模塊的具體結(jié)構(gòu)框圖；
[0024] 圖13為本發(fā)明提供的具體實(shí)施例中原始VSP數(shù)據(jù)Z分量示意圖；
[0025] 圖14為本發(fā)明提供的具體實(shí)施例中地震子波零相位化示意圖；
[0026] 圖15為本發(fā)明提供的具體實(shí)施例中VSP數(shù)據(jù)初至拾取示意圖；
[0027] 圖16為本發(fā)明提供的具體實(shí)施例中利用VSP數(shù)據(jù)計(jì)算的TAR因子的計(jì)算結(jié)果圖；
[0028] 圖17為本發(fā)明提供的具體實(shí)施例中TAR的計(jì)算差異曲線圖；
[0029] 圖18為本發(fā)明提供的具體實(shí)施例中沿下行波時(shí)距曲線拉平后的VSP數(shù)據(jù)示意 圖；
[0030] 圖19為本發(fā)明提供的具體實(shí)施例中一個(gè)窗內(nèi)的樣點(diǎn)序列不意圖；
[0031] 圖20為本發(fā)明提供的具體實(shí)施例中依VSP數(shù)據(jù)規(guī)律的加權(quán)系數(shù)序列示意圖；
[0032] 圖21為本發(fā)明提供的具體實(shí)施例中加權(quán)之后一個(gè)窗內(nèi)的樣點(diǎn)序列示意圖；
[0033] 圖22為本發(fā)明提供的具體實(shí)施例中自適應(yīng)去異常后的樣點(diǎn)序列示意圖；
[0034] 圖23為本發(fā)明提供的具體實(shí)施例中定向加權(quán)的、自適應(yīng)去異常的、邊界變道數(shù)的 中值濾波后VSP數(shù)據(jù)示意圖1 ;
[0035] 圖24為本發(fā)明提供的具體實(shí)施例中定向加權(quán)的、自適應(yīng)去異常的、邊界變道數(shù)的 中值濾波后VSP數(shù)據(jù)示意圖2。
【具體實(shí)施方式】
[0036] 下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完 整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；?本發(fā)明中的實(shí)施例，本領(lǐng)域普通