[發(fā)明專利]一種三維剩余靜校正方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201410852329.0
• 申請(qǐng)日: 2014-12-31
• 公開（公告）號(hào): CN104483705A
• 公開（公告）日: 2015-04-01
• 發(fā)明（設(shè)計(jì)）人: 馮發(fā)全
• 申請(qǐng)（專利權(quán)）人: 中國石油天然氣集團(tuán)公司;中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司
• 主分類號(hào): G01V1/36
• 分類號(hào): G01V1/36
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 黨曉林;李永強(qiáng)

技術(shù)領(lǐng)域

本發(fā)明涉及地球物理勘探技術(shù)領(lǐng)域，特別涉及一種三維剩余靜校正方法。

背景技術(shù)

剩余靜校正是對(duì)地震資料所作的校正，用于補(bǔ)償由高程、風(fēng)化層厚度以及風(fēng)化層速度產(chǎn)生的影響，把資料校到一個(gè)指定的基準(zhǔn)面上。它是一種地震勘探中消除近地表影響從而提高地震反射波成像質(zhì)量的主要數(shù)據(jù)處理方法。

現(xiàn)有技術(shù)中：

2004年井西利等人公開了一種自適應(yīng)混合反演剩余靜校正的方法,針對(duì)剩余靜校正問題的非線性特點(diǎn)，分析了已有的線性反演及非線性反演方法的優(yōu)缺點(diǎn)，并在此基礎(chǔ)上提出了一種混合反演方法；

2005年李全勝等人公開了一種基于模型道的地表一致性剩余靜校正方法，論述了模型道的構(gòu)造方法，將剩余靜校正問題作為廣義線性反問題來求解炮點(diǎn)和檢波點(diǎn)剩余靜校正量；

2006年段云卿公開了一種折射波剩余靜校正方法，利用折射波信噪比高的特點(diǎn)，將反射波剩余靜校正原理應(yīng)用于折射波處理，通過交互手段，逐段估算折射波的速度，用合適的速度對(duì)地震記錄進(jìn)行線性動(dòng)校正，在共炮點(diǎn)或共中心點(diǎn)道集上，用相關(guān)方法計(jì)算各道與模型道時(shí)差，再用統(tǒng)計(jì)方法計(jì)算出各炮點(diǎn)及檢波點(diǎn)的剩余靜校正量；

2011年潘樹林公開了一種共炮(檢)點(diǎn)剩余靜校正方法，基于經(jīng)過動(dòng)校正后的共炮點(diǎn)數(shù)據(jù)集和共檢波點(diǎn)數(shù)據(jù)集分別求取檢波點(diǎn)剩余靜校正量和炮點(diǎn)剩余靜校正量,并分別將共炮、檢點(diǎn)道集動(dòng)校正后的疊加道作為模型道與道集內(nèi)各個(gè)道進(jìn)行互相關(guān),進(jìn)而求出各個(gè)炮、檢波點(diǎn)的剩余靜校正量。

針對(duì)三維觀測(cè)系統(tǒng)，上述剩余靜校正方法都依賴地震處理系統(tǒng)，即需要在對(duì)地震資料進(jìn)行一定處理后才能應(yīng)用，所要計(jì)算的數(shù)據(jù)量大，因此對(duì)軟件和硬件都提出了較高要求，如需要配置工作站和資料處理系統(tǒng)。由于資料處理中心一般擁有大、中型計(jì)算機(jī)，數(shù)據(jù)處理能力強(qiáng)，故上述剩余靜校正方法一般可以在資料處理中心進(jìn)行。而一般情況下，野外現(xiàn)場(chǎng)計(jì)算機(jī)的數(shù)據(jù)處理能力較弱，故針對(duì)三維觀測(cè)系統(tǒng)，上述剩余靜校正方法無法在野外現(xiàn)場(chǎng)得到實(shí)際應(yīng)用。

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

本申請(qǐng)實(shí)施例的目的在于提供一種三維剩余靜校正方法，以能夠在野外現(xiàn)場(chǎng)或資料處理中心針對(duì)三維觀測(cè)系統(tǒng)進(jìn)行剩余靜校正。

本申請(qǐng)實(shí)施例提供的一種三維剩余靜校正方法是這樣實(shí)現(xiàn)的：

一種三維剩余靜校正方法，包括：

1)利用表層資料或者大炮初至資料建立表層模型，根據(jù)所述表層模型進(jìn)行基準(zhǔn)面靜校正；

2)對(duì)拾取的每個(gè)地震道的初至?xí)r間數(shù)據(jù)施加該地震道對(duì)應(yīng)的炮點(diǎn)和檢波點(diǎn)的基準(zhǔn)面靜校正，生成新的地震道的初至?xí)r間T_s,j，其中，T_s,j表示第s炮點(diǎn)第j檢波點(diǎn)施加基準(zhǔn)面靜校正后的地震道的初至?xí)r間；

3)選取檢波點(diǎn)，確定每個(gè)檢波點(diǎn)所對(duì)應(yīng)的炮點(diǎn)線號(hào)范圍，根據(jù)每個(gè)檢波點(diǎn)所對(duì)應(yīng)的炮點(diǎn)線號(hào)范圍內(nèi)的炮點(diǎn)生成該檢波點(diǎn)的共檢波點(diǎn)道集，確定每個(gè)共檢波點(diǎn)道集所對(duì)應(yīng)的地震道的初至?xí)r間T_r,i，其中，T_r,i表示第r檢波點(diǎn)第i炮點(diǎn)地震道的初至?xí)r間；

4)選取炮點(diǎn)，生成每個(gè)炮點(diǎn)的共炮點(diǎn)道集，確定每個(gè)共炮點(diǎn)道集折射層范圍內(nèi)各檢波點(diǎn)處的初至擬合時(shí)差Δt_s,j，并將該初至擬合時(shí)差作為各檢波點(diǎn)的單次剩余靜校正，其中，Δt_s,j表示第s炮點(diǎn)第j檢波點(diǎn)的單次剩余靜校正；

5)計(jì)算檢波點(diǎn)的剩余靜校正ΔRt_j，其中，

ΔRt_j表示第j檢波點(diǎn)的剩余靜校正；

n表示在共炮點(diǎn)道集上計(jì)算了該檢波點(diǎn)單次剩余靜校正的炮點(diǎn)個(gè)數(shù)；

s表示三維測(cè)線內(nèi)炮點(diǎn)的位置編號(hào)，范圍為在該炮點(diǎn)的共炮點(diǎn)道集上計(jì)算了第j檢波點(diǎn)的單次剩余靜校正的炮點(diǎn)個(gè)數(shù)n；

j表示三維測(cè)線內(nèi)檢波點(diǎn)的位置編號(hào)，范圍為從第一條檢波點(diǎn)線的第一個(gè)檢波點(diǎn)到最后一條檢波點(diǎn)線的最后一個(gè)檢波點(diǎn)；

6)對(duì)于步驟3)中生成的共檢波點(diǎn)道集，確定每個(gè)共檢波點(diǎn)道集折射層范圍內(nèi)各炮點(diǎn)處的初至擬合時(shí)差Δt_r,i，并將該初至擬合時(shí)差作為各炮點(diǎn)的單次剩余靜校正，其中，Δt_r,i表示第r檢波點(diǎn)第i炮點(diǎn)的單次剩余靜校正；

7)計(jì)算炮點(diǎn)的剩余靜校正量ΔSt_i，其中，

ΔSt_i表示第i炮點(diǎn)的剩余靜校正；

m表示在共檢波點(diǎn)道集上計(jì)算了該炮點(diǎn)單次剩余靜校正的檢波點(diǎn)個(gè)數(shù)；

i表示三維測(cè)線內(nèi)炮點(diǎn)的位置編號(hào)，范圍為第一條炮點(diǎn)線的第一個(gè)炮點(diǎn)到最后一條炮點(diǎn)線的最后一個(gè)炮點(diǎn)；