一種微小斷層檢測(cè)方法及斷層檢測(cè)裝置制造方法
【專(zhuān)利摘要】本發(fā)明提供一種微小斷層的檢測(cè)方法及斷層檢測(cè)裝置。所述的微小斷層檢測(cè)方法包括以下處理步驟：將采集的地震數(shù)據(jù)記錄進(jìn)行高保真數(shù)據(jù)處理；對(duì)高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波；將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別；將預(yù)測(cè)前地震數(shù)據(jù)與經(jīng)過(guò)縮放的預(yù)測(cè)后的地震數(shù)據(jù)的振幅對(duì)應(yīng)相減；對(duì)預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析。利用本發(fā)明，可以識(shí)別出地震數(shù)據(jù)中微小斷層的特征，進(jìn)而識(shí)別出微小斷層。
【專(zhuān)利說(shuō)明】—種微小斷層檢測(cè)方法及斷層檢測(cè)裝置
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震勘探中地震數(shù)據(jù)解釋領(lǐng)域，特別涉及一種微小斷層檢測(cè)方法及斷層檢測(cè)裝置。
【背景技術(shù)】
[0002]為了勘測(cè)開(kāi)發(fā)出地層中更多的油氣資源，就需要尋找出更多更隱蔽、更細(xì)小的構(gòu)造油氣田。油氣資源多分布在斷層中。地殼巖層因受力達(dá)到一定強(qiáng)度而發(fā)生破裂，并沿破裂面有明顯相對(duì)移動(dòng)的構(gòu)造稱斷層。斷層既是油氣藏的邊界，也是油氣資源運(yùn)移、聚集的通道，因此斷層的識(shí)別是油氣勘探開(kāi)發(fā)中非常重要的工作。
[0003]識(shí)別斷層首要步驟是獲得高保真地震數(shù)據(jù)，包括在野外加大地震數(shù)據(jù)的采集密度、在室內(nèi)對(duì)采集的地震數(shù)據(jù)做高保真數(shù)據(jù)處理等過(guò)程。第一步野外加大地震數(shù)據(jù)采集密度，主要是指增加不同炮點(diǎn)的震源激發(fā)次數(shù)和增大接收地震波的檢波點(diǎn)的分布密度，目的是從采集上提高地震數(shù)據(jù)的采集密度，最大可能的采集和記錄到地下斷層的反射信號(hào)。第二步是對(duì)采集后的地震數(shù)據(jù)記錄做高保真數(shù)據(jù)處理。高保真數(shù)據(jù)處理是地震數(shù)據(jù)處理中指一系列提高地震數(shù)據(jù)信噪比、提高地震數(shù)據(jù)分辨率和地震數(shù)據(jù)保持振幅的處理過(guò)程。
[0004]其次，分析上述獲得的高保真地震數(shù)據(jù)，通過(guò)對(duì)地震數(shù)據(jù)振幅、相位、波形等信息的分析，盡可能在地震數(shù)據(jù)解釋中準(zhǔn)確識(shí)別出斷層。
[0005]目前在高保真地震數(shù)據(jù)中檢測(cè)斷層的主要方法為相干體斷層識(shí)別技術(shù)。這種方法是利用在出現(xiàn)斷層、地層巖性突變、特殊地質(zhì)體的地帶，地震道之間的波形特征發(fā)生變化，進(jìn)而導(dǎo)致局部的道與道之間相關(guān)性的發(fā)生突變，相干值較小，然后求得所有樣點(diǎn)與其周邊樣點(diǎn)的相關(guān)值形成一個(gè)相干數(shù)據(jù)體，通過(guò)對(duì)相干數(shù)據(jù)體作水平切片來(lái)揭示斷層的地質(zhì)現(xiàn)象。地震數(shù)據(jù)采集時(shí)，每個(gè)接收地震波信息的檢波點(diǎn)所記錄的地震數(shù)據(jù)稱為地震道，或道。相干體斷層識(shí)別技術(shù)的實(shí)施過(guò)程主要包括三個(gè)步驟，第一步是將獲得的高保真三維地震數(shù)據(jù)做偏移處理。偏移處理能夠提高地震剖面分辨率，恢復(fù)子波波形和振幅特征，有利于地震數(shù)據(jù)解釋。第二步是在經(jīng)過(guò)偏移獲得的三維數(shù)據(jù)體中，對(duì)每一個(gè)樣點(diǎn)求得與周?chē)鷺狱c(diǎn)數(shù)據(jù)的相干值，形成一個(gè)相干數(shù)據(jù)體。第三步是對(duì)第二步形成的相干數(shù)據(jù)體做切片分析，進(jìn)行構(gòu)造和巖性解釋。
[0006]上述所述的斷層中，微小斷層是指斷距很小(通常不超過(guò)10米)且即使在信噪比很高的情況下也難以識(shí)別的斷層。相干體斷層識(shí)別技術(shù)是通過(guò)識(shí)別出地震數(shù)據(jù)中因斷層引起的地震數(shù)據(jù)的振幅突變、波形突變等特征，進(jìn)一步通過(guò)切片分析識(shí)別出斷層。而在地震數(shù)據(jù)中微小斷層兩側(cè)的波形、振幅、相位等特征并沒(méi)有發(fā)生顯著的特征改變。在地震數(shù)據(jù)記錄中即使有微小斷層存在，微小斷層處地震波特性的相干值變化不大，找不到微小斷層的特征，也就無(wú)法識(shí)別微小斷層，影響了微小斷層所在的油氣藏的勘探和開(kāi)發(fā)。

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

[0007]本發(fā)明目的在于提供一種微小斷層檢測(cè)方法及斷層檢測(cè)裝置，將預(yù)測(cè)濾波后地震數(shù)據(jù)進(jìn)行縮放和預(yù)測(cè)前地震數(shù)據(jù)比較，用預(yù)測(cè)前地震數(shù)據(jù)與頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)樣點(diǎn)據(jù)的振幅對(duì)應(yīng)相減，獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)，該方法得到的預(yù)測(cè)濾波殘差數(shù)據(jù)能體現(xiàn)出微小斷層的特征，在切片分析時(shí)能夠檢測(cè)出微小斷層，解決了相干體斷層識(shí)別技術(shù)無(wú)法識(shí)別微小斷層特性的技術(shù)問(wèn)題。
[0008]本發(fā)明提供的一種微小斷層檢測(cè)方法，其實(shí)現(xiàn)過(guò)程包括以下處理步驟:
[0009]S1:將采集的地震數(shù)據(jù)記錄進(jìn)行高保真數(shù)據(jù)處理；
[0010]S2:對(duì)高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波，其所述在頻率空間域進(jìn)行預(yù)測(cè)濾波的過(guò)程包括以下處理步驟:
[0011]S201:將所述預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換從時(shí)間域變換到頻率空間域，形成空間域矢量數(shù)據(jù)；
[0012]S202:在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子；
[0013]S203:在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)；
[0014]S204:將所述頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域；
[0015]S3:將頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別，獲得預(yù)測(cè)后縮放數(shù)據(jù)，其所述將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別的過(guò)程，包括以下處理步驟:
[0016]所述的經(jīng)過(guò)頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng)的系數(shù)；
[0017]所述系數(shù)為預(yù)測(cè)前地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)后地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值；
[0018]S4:將預(yù)測(cè)前地震數(shù)據(jù)與所述預(yù)測(cè)后縮放數(shù)據(jù)的樣點(diǎn)的振幅對(duì)應(yīng)相減；
[0019]S5:對(duì)預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后縮放數(shù)據(jù)樣點(diǎn)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析。
[0020]上述所述的微小斷層檢測(cè)方法，其優(yōu)選方案為，還包括為S4中預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)定一個(gè)門(mén)檻值，將所述預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的數(shù)據(jù)用零代替。
[0021]上述所述的微小斷層檢測(cè)方法，其優(yōu)選方案為，S3中所述將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別的過(guò)程，包括以下處理步驟:
[0022]所述預(yù)測(cè)前地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng)的系數(shù)；
[0023]所述系數(shù)為預(yù)測(cè)后地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)前地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值。
[0024]本發(fā)明利用本發(fā)明檢查微小斷層的方法提供一種斷層檢測(cè)裝置，所述斷層檢測(cè)裝置包括高保真數(shù)據(jù)處理模塊，頻率空間域?yàn)V波模塊，振幅調(diào)整模塊，振幅相減模塊，切片分析模塊；其中，
[0025]高保真數(shù)據(jù)處理模塊，用于提高地震的數(shù)據(jù)信噪比和分辨率；
[0026]頻率空間域?yàn)V波模塊，用于將高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波；
[0027]振幅調(diào)整模塊，用于將高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)和頻率空間域?yàn)V波模塊獲得的預(yù)測(cè)后地震數(shù)據(jù)的振幅調(diào)整到同一振幅級(jí)別；
[0028]振幅相減模塊，用于將振幅調(diào)整模塊獲得的預(yù)測(cè)后縮放數(shù)據(jù)與高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)對(duì)應(yīng)地震道的振幅相減；
[0029]上述所述的頻率空間域?yàn)V波模塊包括傅里葉變換模塊、維納濾波模塊、預(yù)測(cè)數(shù)據(jù)模塊、傅里葉逆變換模塊；其中，
[0030]傅里葉變換模塊，用于將高保真數(shù)據(jù)模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換從時(shí)間域變換到頻率空間域，形成空間域矢量數(shù)據(jù)；
[0031]維納濾波模塊，用于在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子；
[0032]預(yù)測(cè)數(shù)據(jù)模塊，用于在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)；
[0033]傅里葉逆變換模塊，用于將預(yù)測(cè)數(shù)據(jù)模塊獲得的頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域；
[0034]上述所述振幅調(diào)整模塊包括系數(shù)計(jì)算模塊和數(shù)據(jù)縮放模塊；其中，
[0035]系數(shù)計(jì)算模塊，用于計(jì)算預(yù)測(cè)前地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)后地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值；
[0036]數(shù)據(jù)縮放模塊，用于預(yù)測(cè)后地震數(shù)據(jù)中每一個(gè)地震道的采樣點(diǎn)的振幅乘以由所述系數(shù)計(jì)算模塊計(jì)算得出的與其地震道相對(duì)應(yīng)的系數(shù)。
[0037]上述所述的微小斷層檢測(cè)裝置，其優(yōu)選方案為，在所述振幅相減模塊中，還包括門(mén)檻值處理模塊，用于濾除所述預(yù)測(cè)濾波殘差數(shù)據(jù)中的干擾信息，所述門(mén)檻值處理模塊包括門(mén)檻值設(shè)定模塊和門(mén)檻值處理模塊；其中，
[0038]門(mén)檻值設(shè)定模塊，用于對(duì)所述預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)置一個(gè)門(mén)檻值；
[0039]門(mén)檻值處理模塊，用于將所述預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的振幅
用零代替。
[0040]上述所述的微小斷層檢測(cè)裝置，其優(yōu)選方案可以為，所述系數(shù)計(jì)算模塊，用于計(jì)算預(yù)測(cè)后地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)前地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值；
[0041]所述數(shù)據(jù)縮放模塊，用于預(yù)測(cè)前地震數(shù)據(jù)中每一個(gè)地震道的采樣點(diǎn)的振幅乘以由所述系數(shù)計(jì)算模塊計(jì)算得出的與其地震道相對(duì)應(yīng)的系數(shù)。
[0042]本發(fā)明首先將地震數(shù)據(jù)在頻率空間域進(jìn)行濾波預(yù)測(cè)獲得預(yù)測(cè)后地震數(shù)據(jù)，然后將預(yù)測(cè)前后地震數(shù)據(jù)相減獲得的預(yù)測(cè)殘差數(shù)據(jù)，從預(yù)測(cè)殘差數(shù)據(jù)的切片上識(shí)別微小斷層。本發(fā)明基于在頻率空間域線性同相軸是可以預(yù)測(cè)的原理，在頻率空間域線性同相軸存在微小斷層的地方預(yù)測(cè)效果不理想，當(dāng)預(yù)測(cè)前后地震數(shù)據(jù)相減時(shí)，存在微小斷層的地方就會(huì)出現(xiàn)比較明顯的非零殘差值，而在沒(méi)有微小斷層的地方，預(yù)測(cè)效果較好，預(yù)測(cè)殘差的理論值為零。本發(fā)明方法將微小斷層的識(shí)別與頻率空間域地震數(shù)據(jù)的預(yù)測(cè)殘差值聯(lián)系起來(lái)，將微小斷層的特征在預(yù)測(cè)殘差數(shù)據(jù)中體現(xiàn)出來(lái)，再通過(guò)對(duì)預(yù)測(cè)殘差數(shù)據(jù)的切片分析結(jié)合地震解釋資料以及當(dāng)?shù)貙?shí)際地質(zhì)情況識(shí)別出微小斷層。
【專(zhuān)利附圖】

【附圖說(shuō)明】
[0043]圖1是本發(fā)明實(shí)施例1提供的微小斷層檢測(cè)方法流程圖；
[0044]圖2是本發(fā)明實(shí)施例1中的一段預(yù)測(cè)前地震數(shù)據(jù)記錄圖；
[0045]圖3是本發(fā)明實(shí)施例1中圖2在頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)記錄圖；
[0046]圖4是本發(fā)明實(shí)施例1中圖2與圖3相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)圖；
[0047]圖5是本發(fā)明實(shí)施例1中圖4經(jīng)門(mén)檻值處理后的預(yù)測(cè)濾波殘差數(shù)據(jù)圖；
[0048]圖6是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置的模塊結(jié)構(gòu)示意圖；
[0049]圖7是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中頻率空間域?yàn)V波模塊的模塊結(jié)構(gòu)示意圖；
[0050]圖8是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中振幅調(diào)整模塊的模塊結(jié)構(gòu)示意圖；
[0051]圖9是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中門(mén)檻值處理模塊的模塊結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0052]為了使本【技術(shù)領(lǐng)域】的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí)施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述。顯然，所描述的實(shí)施例僅僅是本申請(qǐng)中的一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其它實(shí)施例，都應(yīng)當(dāng)屬于本發(fā)明保護(hù)的范圍。
[0053]實(shí)施例1是本發(fā)明提供的一種在高保真地震數(shù)據(jù)背景下檢測(cè)微小斷層的方法，圖1是本發(fā)明實(shí)施例1提供的微小斷層檢測(cè)方法流程圖。如圖1所示，一種微小斷層檢測(cè)方法，其實(shí)施過(guò)程包括以下處理步驟:
[0054]S1:將采集的地震數(shù)據(jù)記錄進(jìn)行高保真數(shù)據(jù)處理。
[0055]野外加大地震數(shù)據(jù)的采集密度獲得地震數(shù)據(jù)記錄，對(duì)采集獲得的地震數(shù)據(jù)記錄做高保真數(shù)據(jù)處理，獲得高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)記錄。
[0056]上述所述的高保真數(shù)據(jù)處理是指一系列提高地震數(shù)據(jù)信噪比和分辨率的處理過(guò)程。實(shí)際野外采集的地震數(shù)據(jù)記錄中，由于激發(fā)和接收因素(震源激發(fā)深度不同、炮檢距不同等)、地層層狀結(jié)構(gòu)、各種干擾波等因素的影響，造成原始的地震數(shù)據(jù)記錄信噪比和分辨率較低，不適用于后續(xù)的地震數(shù)據(jù)處理和解釋。所述的炮檢距是檢波點(diǎn)到炮點(diǎn)的距離。檢波點(diǎn)是地震數(shù)據(jù)采集時(shí)接收數(shù)據(jù)的檢波器所在的點(diǎn)。炮點(diǎn)是地震數(shù)據(jù)采集時(shí)激發(fā)地震波的震源所在的點(diǎn)。高保真數(shù)據(jù)處理就是在地震數(shù)據(jù)處理時(shí)，為盡量排除各種非地質(zhì)因素引起的干擾和誤差而進(jìn)行的一系列靜校正、動(dòng)校正、反褶積、保持振幅等提高地震數(shù)據(jù)信噪比和分辨率的處理過(guò)程。例如，靜校正是按靜校正量所確定的移動(dòng)量，將地震道的振幅離散值進(jìn)行整體移動(dòng)，消除地面山谷、地表風(fēng)化層和低速帶等地表異常的影響。動(dòng)校正是通過(guò)動(dòng)校正公式計(jì)算出動(dòng)校正量來(lái)消除地震波到達(dá)各不同檢波點(diǎn)時(shí)產(chǎn)生的正常時(shí)差。檢波器接收到的地震信號(hào)可以看成是反射系數(shù)序列與地震子波的褶積，子波的存在會(huì)降低分辨率，反褶積是將接收到的地震信號(hào)的子波進(jìn)行壓縮，最后在地震道上僅保留反射系數(shù)，提高地震資料分辨率。地震波的能量在傳播過(guò)程中受吸收和衰減作用的影響，使得淺、中、深層能量差異較大，保持振幅是通過(guò)使用球面擴(kuò)散補(bǔ)償、剩余振幅補(bǔ)償?shù)燃夹g(shù)盡量消除地震數(shù)據(jù)記錄中因空間、時(shí)間因素而產(chǎn)生的能量差異。
[0057]對(duì)野外采集的地震數(shù)據(jù)記錄進(jìn)行高保真數(shù)據(jù)處理，獲得高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)記錄。
[0058]S2:對(duì)高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波。
[0059]地震數(shù)據(jù)在頻率空間域的濾波是對(duì)地震數(shù)據(jù)記錄中隨機(jī)噪聲壓制的過(guò)程，其濾波方法可以參考申請(qǐng)?zhí)枮?00410102646.7的中國(guó)專(zhuān)利。
[0060]地震數(shù)據(jù)記錄上各道振動(dòng)相位相同的振幅極值(俗稱波峰或波谷)的連線稱為同相軸。在時(shí)間空間域，同相軸的頻率相同或相近，當(dāng)變換到頻率空間域時(shí)，同相軸則近似為一條直線。當(dāng)已知所述頻率空間域的同相軸直線上的幾個(gè)點(diǎn)時(shí)，該直線上其它點(diǎn)的數(shù)值可以通過(guò)一些特定的算法預(yù)測(cè)出來(lái)，這種特性稱為地震數(shù)據(jù)在頻率空間域的線性同相軸是可以預(yù)測(cè)的。
[0061]如上述所述專(zhuān)利《壓制低信噪比地震記錄中隨機(jī)噪聲的方法》中壓制隨機(jī)噪聲的方法，是基于地震數(shù)據(jù)在頻率空間域線性同相軸是可以預(yù)測(cè)的原理。所述預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波的方法主要過(guò)程包括以下步驟:
[0062]S201:將所述預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換將其從時(shí)間域變換到頻率空間域，形成地震數(shù)據(jù)中振幅有大小和方向的空間域矢量數(shù)據(jù)；
[0063]S202:在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子；
[0064]S203:在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)；
[0065]S204:將所述頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域，形成一個(gè)經(jīng)過(guò)頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)。
[0066]經(jīng)高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行濾波后，形成經(jīng)過(guò)頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)。
[0067]S3:將頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別，獲得預(yù)測(cè)后縮放數(shù)據(jù)。
[0068]將S2所述的經(jīng)過(guò)頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng)的系數(shù)，該系數(shù)將每個(gè)地震道采樣點(diǎn)的均方根振幅縮放成與預(yù)測(cè)前對(duì)應(yīng)的原始地震道采樣點(diǎn)的均方根振幅一致，獲得頻率空間域預(yù)測(cè)濾波后的縮放數(shù)據(jù)。該過(guò)程處理的目的是將地震數(shù)據(jù)濾波前后的振幅調(diào)整到一個(gè)能量級(jí)別上，當(dāng)所需處理的預(yù)測(cè)前地震數(shù)據(jù)和預(yù)測(cè)后地震數(shù)據(jù)相對(duì)應(yīng)地震道的均方根振幅的平均值一致的時(shí)候，可以認(rèn)為所述預(yù)測(cè)前地震數(shù)據(jù)和預(yù)測(cè)后地震數(shù)據(jù)的振幅為同一個(gè)振幅級(jí)別。
[0069]預(yù)測(cè)后地震數(shù)據(jù)的每一個(gè)地震道都有一個(gè)與其相對(duì)應(yīng)的系數(shù)，該系數(shù)為預(yù)測(cè)前地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)后地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值。若頻率空間域?yàn)V波后的預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)對(duì)應(yīng)地震道的振幅差值變化不大，即振幅變化誤差在地震數(shù)據(jù)解釋人員可接受的范圍內(nèi)，此時(shí)認(rèn)為上述所述的系數(shù)的值為1，預(yù)測(cè)后地震數(shù)據(jù)可以實(shí)際上不再進(jìn)行縮放。
[0070]上述所述的將預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別，也可以是預(yù)測(cè)前地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng)的系數(shù)，此時(shí)的系數(shù)為預(yù)測(cè)后地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)前地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值。
[0071]將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)的振幅級(jí)別調(diào)整到與預(yù)測(cè)前地震數(shù)據(jù)的振幅級(jí)別一致，獲得預(yù)測(cè)后縮放數(shù)據(jù)。
[0072]S4:將預(yù)測(cè)前地震數(shù)據(jù)與所述預(yù)測(cè)后縮放數(shù)據(jù)的樣點(diǎn)的振幅對(duì)應(yīng)相減。
[0073]預(yù)測(cè)前地震數(shù)據(jù)與上述所述獲得預(yù)測(cè)后縮放數(shù)據(jù)的樣點(diǎn)的振幅對(duì)應(yīng)相減，獲得預(yù)測(cè)濾波殘差數(shù)據(jù)?；陬l率空間域線性同相軸是可以預(yù)測(cè)的原理，在頻率空間域線性同相軸存在微小斷層的地方預(yù)測(cè)效果不理想。當(dāng)預(yù)測(cè)前后地震數(shù)據(jù)相減時(shí)，在沒(méi)有微小斷層的地方，預(yù)測(cè)效果較好，預(yù)測(cè)殘差的理論值為零，而存在微小斷層的地方就會(huì)出現(xiàn)比較明顯的非零殘差值。本發(fā)明方法將微小斷層的特征的識(shí)別與頻率空間域地震數(shù)據(jù)的預(yù)測(cè)殘差值聯(lián)系起來(lái)。
[0074]所述的預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后縮放數(shù)據(jù)樣點(diǎn)的振幅對(duì)應(yīng)相減，可以是預(yù)測(cè)前地震數(shù)據(jù)減去預(yù)測(cè)后地震數(shù)據(jù)，也可以是預(yù)測(cè)后地震數(shù)據(jù)減去預(yù)測(cè)前地震數(shù)據(jù)。
[0075]在此可以為上述所述的預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)定一個(gè)門(mén)檻值，將預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的數(shù)據(jù)用零代替，形成門(mén)檻值過(guò)濾后的預(yù)測(cè)濾波殘差數(shù)據(jù)，使預(yù)測(cè)濾波殘差數(shù)據(jù)中微小斷層體現(xiàn)出來(lái)的特征更加明顯。所述門(mén)檻值為經(jīng)驗(yàn)值，可以由所述領(lǐng)域技術(shù)人員經(jīng)多次實(shí)驗(yàn)后選擇能較好體現(xiàn)出微小斷層的特征的數(shù)值。
[0076]將預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后縮放數(shù)據(jù)樣點(diǎn)的振幅對(duì)應(yīng)相減，獲得預(yù)測(cè)濾波殘差數(shù)據(jù)。
[0077]S5:對(duì)預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后縮放數(shù)據(jù)的樣點(diǎn)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析。
[0078]對(duì)上述所述的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析，進(jìn)而從切片上識(shí)別出微小斷層。切片是沿著三維地震數(shù)據(jù)體的一個(gè)方向，以平面或曲面抽取的具有地質(zhì)物理意義的二維空間的數(shù)據(jù)體。常用切片分析方法為水平切片。水平切片是一種等時(shí)平面圖，是地下不同層位的地質(zhì)信息在同一時(shí)間平面內(nèi)的反映。斷層在水平切片上的特征主要是振幅呈規(guī)律性分布，由微小斷層引起的預(yù)測(cè)濾波殘差數(shù)據(jù)在水平切片上體現(xiàn)出的特征主要為規(guī)律性出現(xiàn)的一段同相軸。通過(guò)對(duì)預(yù)測(cè)濾波殘差數(shù)據(jù)的切片分析，識(shí)別出地震數(shù)據(jù)中的微小斷層。
[0079]通過(guò)對(duì)所述預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析，從切片上找到微小斷層體現(xiàn)出的特征，識(shí)別出微小斷層。
[0080]圖2至圖5是使用本發(fā)明方法在高保真地震數(shù)據(jù)記錄中識(shí)別出微小斷層的過(guò)程。首先，對(duì)野外采集的地震數(shù)據(jù)記錄做進(jìn)行高保真數(shù)據(jù)處理后形成的預(yù)測(cè)前地震數(shù)據(jù)，如圖2所示。如圖2所示，橫軸為地震數(shù)據(jù)記錄的地震道號(hào)，如圖2中地震道號(hào)1980至2080段。對(duì)每一個(gè)地震道進(jìn)行的編號(hào)稱為地震道號(hào)?？v軸為該地震數(shù)據(jù)的記錄時(shí)間，單位為毫秒，如圖2中的600毫秒至1500毫秒。如圖2所示，該預(yù)測(cè)前地震數(shù)據(jù)包括兩個(gè)明顯的同相軸，一個(gè)是在1000毫米處的水平同相軸，另一個(gè)是從1020ms到1500ms的傾斜同相軸。該預(yù)測(cè)前地震數(shù)據(jù)記錄的水平同相軸上有兩個(gè)微小斷層，傾斜同相軸有兩個(gè)微小斷層，四個(gè)微小斷層從圖2中都很難直觀識(shí)別出，即使做進(jìn)一步分析也難以識(shí)別。其次，對(duì)圖2預(yù)測(cè)前地震數(shù)據(jù)用本發(fā)明方法中S2所述的濾波方法在頻率空間域進(jìn)行預(yù)測(cè)濾波處理，獲得預(yù)測(cè)后地震數(shù)據(jù)，如圖3所示。再次，將頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別。對(duì)比圖3與圖2的數(shù)據(jù)可以發(fā)現(xiàn)，圖3預(yù)測(cè)后地震數(shù)據(jù)整體上振幅沒(méi)有發(fā)生明顯變化，與預(yù)測(cè)前地震數(shù)據(jù)各個(gè)相對(duì)應(yīng)的地震道的振幅相比比較均衡，故在本實(shí)施例1中圖3的預(yù)測(cè)后數(shù)據(jù)振幅縮放系數(shù)取I。然后，用圖2預(yù)測(cè)前地震數(shù)據(jù)減去圖3預(yù)測(cè)后地震數(shù)據(jù)，形成預(yù)測(cè)濾波殘差數(shù)據(jù)，如圖4所示。在圖4上可以看見(jiàn)A、B、C、D出有明顯的預(yù)測(cè)濾波殘差數(shù)據(jù)，同時(shí)周邊有較多的預(yù)測(cè)濾波殘差干擾數(shù)據(jù)存在，這些干擾有可能影響解釋人員對(duì)微小斷層的辨別。因此，在本實(shí)施例1中為圖4設(shè)置了 6個(gè)門(mén)榲值，大小分別為50，100，200，300，400，500，用這6個(gè)門(mén)檻值分別對(duì)圖4的數(shù)據(jù)分別進(jìn)行過(guò)濾。對(duì)這6個(gè)門(mén)檻值過(guò)濾的預(yù)測(cè)濾波殘差數(shù)據(jù)的結(jié)果進(jìn)行對(duì)比分析，其中當(dāng)選擇門(mén)檻值為400時(shí)能清楚的顯示出微小斷層在預(yù)測(cè)濾波殘差數(shù)據(jù)中的特征。因此本實(shí)施例1選擇400作為門(mén)檻值。在圖4中，振幅絕對(duì)值大于門(mén)檻值400的振幅被保留下來(lái)，小于門(mén)檻值400的振幅值用零代替，形成門(mén)檻值400過(guò)濾后的預(yù)測(cè)濾波殘差數(shù)據(jù)，如圖5所示。從圖5中可以清晰的識(shí)別出A1、B1、C1、D1處微小斷層體現(xiàn)出的特征。最后，對(duì)圖5預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析，進(jìn)而識(shí)別出微小斷層，達(dá)到本發(fā)明方法實(shí)現(xiàn)檢測(cè)微小斷層的目的。
[0081]本發(fā)明的實(shí)施例2為利用本發(fā)明檢測(cè)微小斷層方法的一種斷層檢測(cè)裝置。圖6示是本發(fā)明實(shí)施例2提供的一種斷層檢測(cè)裝置的模塊結(jié)構(gòu)示意圖。如圖6所示，所述的斷層檢測(cè)裝置包括高保真數(shù)據(jù)處理模塊1，頻率空間域?yàn)V波模塊2，振幅調(diào)整模塊3，振幅相減模塊4，切片分析模塊5。其中，
[0082]高保真數(shù)據(jù)處理模塊1，可以用于對(duì)采集的地震數(shù)據(jù)進(jìn)行提高地震數(shù)據(jù)信噪比和分辨率的處理，獲得預(yù)測(cè)前地震數(shù)據(jù)；
[0083]頻率空間域?yàn)V波模塊2，可以用于將高保真數(shù)據(jù)處理模塊I獲得的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行濾波，獲得預(yù)測(cè)后地震數(shù)據(jù)；
[0084]振幅調(diào)整模塊3，可以用于將高保真數(shù)據(jù)處理模塊I獲得的預(yù)測(cè)前地震數(shù)據(jù)和頻率空間域?yàn)V波模塊2獲得的預(yù)測(cè)后地震數(shù)據(jù)的振幅調(diào)整到同一振幅級(jí)別，獲得預(yù)測(cè)后縮放數(shù)據(jù)；
[0085]振幅相減模塊4，可以用于將振幅調(diào)整模3塊獲得的預(yù)測(cè)后縮放數(shù)據(jù)與高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)對(duì)應(yīng)地震道的振幅相減，獲得預(yù)測(cè)濾波殘差數(shù)據(jù)；
[0086]切片分析模塊5，可以用于對(duì)振幅相減模塊獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析。
[0087]其中，上述所述的所述頻率空間域?yàn)V波模塊2包括傅里葉變換模塊201、維納濾波模塊202、預(yù)測(cè)數(shù)據(jù)模塊203、傅里葉逆變換模塊204。圖7是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中所述頻率空間域?yàn)V波模塊2的模塊結(jié)構(gòu)示意圖。其中，
[0088]傅里葉變換模塊201，可以用于將高保真數(shù)據(jù)模塊I獲得的預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換從時(shí)間域變換到頻率空間域，形成空間域矢量數(shù)據(jù)；
[0089]維納濾波模塊202，可以用于在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子；
[0090]預(yù)測(cè)數(shù)據(jù)模塊203，可以用于在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)；
[0091]傅里葉逆變換模塊204，可以用于將預(yù)測(cè)數(shù)據(jù)模塊203獲得的頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域。
[0092]上述所述振幅調(diào)整模塊3包括系數(shù)計(jì)算模塊301和數(shù)據(jù)縮放模塊302。圖8是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中所述振幅調(diào)整模塊3的模塊結(jié)構(gòu)示意圖。其中，
[0093]系數(shù)計(jì)算模塊301，可以用于計(jì)算所述預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)對(duì)應(yīng)地震道的均方根振幅的平均值的比值；
[0094]數(shù)據(jù)縮放模塊302，可以用于預(yù)測(cè)后地震數(shù)據(jù)中每一個(gè)地震道振幅乘以由所述系數(shù)計(jì)算模塊301計(jì)算得出的與其地震道相對(duì)應(yīng)的系數(shù)，獲得預(yù)測(cè)后縮放數(shù)據(jù)。
[0095]上述所述振幅相減模塊4中，還可以包括門(mén)檻值處理模塊401，用于濾除所述振幅相減模塊獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)中的干擾信息。圖9是本發(fā)明實(shí)施例2提供的斷層檢測(cè)裝置中所述門(mén)檻值處理模塊401的模塊結(jié)構(gòu)示意圖。如圖9所示，所述門(mén)檻值處理模塊401包括門(mén)檻值設(shè)定模塊4011和門(mén)檻值處理模塊4012。其中，
[0096]門(mén)檻值設(shè)定模塊4011，可以用于為所述振幅相減模塊4獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)置一個(gè)門(mén)檻值；
[0097]門(mén)檻值處理模塊4012，可以用于將所述振幅相減模塊4獲得預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的數(shù)據(jù)用零代替。
[0098]本發(fā)明提供的斷層檢測(cè)裝置，通過(guò)各個(gè)功能模塊實(shí)現(xiàn)地震數(shù)據(jù)頻率空間域?yàn)V波、振幅縮放、門(mén)檻值過(guò)濾等獲得精確的預(yù)測(cè)濾波殘差數(shù)據(jù)，然后將斷層的特征在預(yù)測(cè)濾波殘差數(shù)據(jù)中體現(xiàn)出來(lái)，能快速高效的檢測(cè)出斷層。
[0099]需要說(shuō)明的是，所述微小斷層檢測(cè)方法中，所述的預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后的縮放數(shù)據(jù)的振幅對(duì)應(yīng)相減，可以是預(yù)測(cè)前地震數(shù)據(jù)減去預(yù)測(cè)后地震數(shù)據(jù)，也可以是預(yù)測(cè)后地震數(shù)據(jù)減去預(yù)測(cè)前地震數(shù)據(jù)。本發(fā)明對(duì)此并不作出限定。
[0100]需要說(shuō)明的是，所述斷層檢測(cè)裝置中，所述振幅相減模塊4中預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后的縮放數(shù)據(jù)的振幅對(duì)應(yīng)相減，可以是預(yù)測(cè)前地震數(shù)據(jù)減去預(yù)測(cè)后地震數(shù)據(jù)，也可以是預(yù)測(cè)后地震數(shù)據(jù)減去預(yù)測(cè)前地震數(shù)據(jù)。本發(fā)明對(duì)此并不作出限定。
[0101]本發(fā)明提供的一種利用頻率空間域地震數(shù)據(jù)預(yù)測(cè)濾波殘差數(shù)據(jù)檢測(cè)微小斷層的檢測(cè)方法及斷層檢測(cè)裝置，通過(guò)特定處理方法計(jì)算得出預(yù)測(cè)濾波殘差數(shù)據(jù)，再對(duì)預(yù)測(cè)濾波殘差數(shù)據(jù)的切片分析，該方法得到的預(yù)測(cè)濾波殘差數(shù)據(jù)能高效地檢測(cè)出微小斷層，實(shí)現(xiàn)相干體斷層檢測(cè)技術(shù)等無(wú)法檢測(cè)微小斷層的目的，為尋找和挖潛微小斷層構(gòu)造油氣田奠定了理論上的基礎(chǔ)。
【權(quán)利要求】
1.一種微小斷層檢測(cè)方法，其特征在于，包括以下處理步驟: S1:將采集的地震數(shù)據(jù)記錄進(jìn)行高保真數(shù)據(jù)處理； 52:將高保真數(shù)據(jù)處理后的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波，其所述在頻率空間域進(jìn)行預(yù)測(cè)濾波的過(guò)程包括以下處理步驟: 5201:將所述預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換從時(shí)間域變換到頻率空間域，形成空間域矢量數(shù)據(jù)； 5202:在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子； S203:在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)； S204:將所述頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域； 53:將頻率空間域?yàn)V波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別，獲得預(yù)測(cè)后縮放數(shù)據(jù)；其所述將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別的過(guò)程，包括以下處理步驟: 所述的經(jīng)過(guò)頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng) 的系數(shù)； 所述系數(shù)為預(yù)測(cè)前地震數(shù)據(jù)中地震道的均方根振幅值的平均值，與預(yù)測(cè)后地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值； 54:將預(yù)測(cè)前地震數(shù)據(jù)與所述預(yù)測(cè)后縮放數(shù)據(jù)的樣點(diǎn)的振幅對(duì)應(yīng)相減； 55:對(duì)預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后縮放數(shù)據(jù)樣點(diǎn)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析。
2.如權(quán)利要求1所述的微小斷層檢測(cè)方法，其特征在于，還包括為S4中預(yù)測(cè)前地震數(shù)據(jù)與預(yù)測(cè)后地震數(shù)據(jù)的振幅對(duì)應(yīng)相減獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)定一個(gè)門(mén)檻值，將所述預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的數(shù)據(jù)用零代替。
3.如權(quán)利要求1所述的微小斷層檢測(cè)方法，其特征在于，S3中所述將頻率空間域預(yù)測(cè)濾波后的預(yù)測(cè)后地震數(shù)據(jù)與預(yù)測(cè)前地震數(shù)據(jù)的振幅調(diào)整到同一個(gè)振幅級(jí)別的過(guò)程，包括以下處理步驟: 所述預(yù)測(cè)前地震數(shù)據(jù)中的每一個(gè)地震道，乘以一個(gè)與其地震道相對(duì)應(yīng)的系數(shù)； 所述系數(shù)為預(yù)測(cè)后地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)前地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值。
4.一種斷層檢測(cè)裝置，其特征在于，所述斷層檢測(cè)裝置包括高保真數(shù)據(jù)處理模塊，頻率空間域?yàn)V波模塊，振幅調(diào)整模塊，振幅相減模塊，切片分析模塊；其中， 所述高保真數(shù)據(jù)處理模塊，用于提高地震的數(shù)據(jù)信噪比和分辨率； 所述頻率空間域?yàn)V波模塊，用于將高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)在頻率空間域進(jìn)行預(yù)測(cè)濾波； 所述振幅調(diào)整模塊，用于將高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)和頻率空間域?yàn)V波模塊獲得的預(yù)測(cè)后地震數(shù)據(jù)的振幅調(diào)整到同一振幅級(jí)別； 所述振幅相減模塊，用于將振幅調(diào)整模塊獲得的預(yù)測(cè)后縮放數(shù)據(jù)與高保真數(shù)據(jù)處理模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)對(duì)應(yīng)地震道的振幅相減；所述切片分析模塊，用于對(duì)振幅相減模塊獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)做切片分析；上述所述的頻率空間域?yàn)V波模塊包括傅里葉變換模塊、維納濾波模塊、預(yù)測(cè)數(shù)據(jù)模塊、傅里葉逆變換模塊；其中， 所述傅里葉變換模塊，用于將高保真數(shù)據(jù)模塊獲得的預(yù)測(cè)前地震數(shù)據(jù)通過(guò)傅里葉變換從時(shí)間域變換到頻率空間域，形成空間域矢量數(shù)據(jù)； 所述維納濾波模塊，用于在所述空間域矢量數(shù)據(jù)的單個(gè)頻率點(diǎn)上，用維納濾波方法求得預(yù)測(cè)算子； 所述預(yù)測(cè)數(shù)據(jù)模塊，用于在所述空間域矢量數(shù)據(jù)的單點(diǎn)頻率點(diǎn)上，利用被預(yù)測(cè)地震道前后總共一個(gè)預(yù)測(cè)算子長(zhǎng)度的地震道與其對(duì)應(yīng)的預(yù)測(cè)算子逐個(gè)相乘，再將相乘的結(jié)果相加作為被預(yù)測(cè)地震道的單點(diǎn)頻率，完成頻率域所有單點(diǎn)頻率的預(yù)測(cè)，形成頻率空間域的預(yù)測(cè)數(shù)據(jù)； 所述傅里葉逆變換模塊，用于將預(yù)測(cè)數(shù)據(jù)模塊獲得的頻率空間域的預(yù)測(cè)數(shù)據(jù)通過(guò)傅里葉逆變換從頻率域變換到時(shí)間域； 上述所述振幅調(diào)整模塊包括系數(shù)計(jì)算模塊和數(shù)據(jù)縮放模塊；其中， 所述系數(shù)計(jì)算模塊，用于計(jì)算預(yù)測(cè)前地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)后地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值； 所述數(shù)據(jù)縮放模塊，用于預(yù)測(cè)后地震數(shù)據(jù)中每一個(gè)地震道的采樣點(diǎn)的振幅乘以由所述系數(shù)計(jì)算模塊計(jì)算得出的與其地震道相對(duì)應(yīng)的系數(shù)。
5.如權(quán)利要求4所述的一種斷層檢測(cè)裝置，其特征在于，所述振幅相減模塊中還包括門(mén)檻值處理模塊，用于濾除所述振幅相減模塊獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)中的干擾信息，所述門(mén)檻值處理模塊包括門(mén)檻值設(shè)定模塊和門(mén)檻值處理模塊；其中， 所述門(mén)檻值設(shè)定模塊，用于為所述振幅相減模塊獲得的預(yù)測(cè)濾波殘差數(shù)據(jù)設(shè)置一個(gè)門(mén)檻值； 所述門(mén)檻值處理模塊，用于將所述振幅相減模塊獲得預(yù)測(cè)濾波殘差數(shù)據(jù)中振幅絕對(duì)值小于門(mén)檻值的數(shù)據(jù)用零代替。
6.如權(quán)利要求4所述的一種斷層檢測(cè)裝置，其特征在于: 所述系數(shù)計(jì)算模塊，用于計(jì)算預(yù)測(cè)后地震數(shù)據(jù)中地震道的均方根振幅值的平均值與預(yù)測(cè)前地震數(shù)據(jù)與其相對(duì)應(yīng)的地震道的均方根振幅值的平均值的比值； 所述數(shù)據(jù)縮放模塊，用于預(yù)測(cè)前地震數(shù)據(jù)中每一個(gè)地震道的采樣點(diǎn)的振幅乘以由所述系數(shù)計(jì)算模塊計(jì)算得出的與其地震道相對(duì)應(yīng)的系數(shù)。
【文檔編號(hào)】G01V1/28GK103760599SQ201410001435
【公開(kāi)日】2014年4月30日 申請(qǐng)日期:2014年1月2日 優(yōu)先權(quán)日:2014年1月2日
【發(fā)明者】俞壽朋, 熊定鈺, 錢(qián)忠平, 趙波, 崔京彬 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司