一種三維剩余靜校正方法
【專利摘要】本申請實(shí)施例公開了一種三維剩余靜校正方法，屬于地球物理勘探【技術(shù)領(lǐng)域】。該方法包括：將計(jì)算的基準(zhǔn)面靜校正施加到野外拾取的地震初至?xí)r間數(shù)據(jù)體上，得到一個新的地震初至?xí)r間數(shù)據(jù)體，然后在共炮點(diǎn)道集和共檢波點(diǎn)道集的正負(fù)偏移距上，分別選擇劃分折射層范圍，對該道集正負(fù)偏移距上的折射層的初至?xí)r間與偏移距進(jìn)行擬合，得到各地震道的擬合初至?xí)r間，然后求取二者的擬合時差，最后對所有檢波點(diǎn)和炮點(diǎn)的擬合時差各自進(jìn)行統(tǒng)計(jì)計(jì)算，其平均值即為各檢波點(diǎn)和炮點(diǎn)的剩余靜校正。本申請實(shí)施例的方法不依賴地震處理系統(tǒng)，計(jì)算的數(shù)據(jù)量小，對地震資料信噪比要求不高，便于在現(xiàn)場微機(jī)上計(jì)算剩余靜校正，能夠解決過大的剩余靜校正。
【專利說明】一種三維剩余靜校正方法

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

【背景技術(shù)】
[0002] 剩余靜校正是對地震資料所作的校正，用于補(bǔ)償由高程、風(fēng)化層厚度以及風(fēng)化層 速度產(chǎn)生的影響，把資料校到一個指定的基準(zhǔn)面上。它是一種地震勘探中消除近地表影響 從而提高地震反射波成像質(zhì)量的主要數(shù)據(jù)處理方法。
[0003] 現(xiàn)有技術(shù)中：
[0004] 2004年井西利等人公開了一種自適應(yīng)混合反演剩余靜校正的方法，針對剩余靜 校正問題的非線性特點(diǎn)，分析了已有的線性反演及非線性反演方法的優(yōu)缺點(diǎn)，并在此基礎(chǔ) 上提出了一種混合反演方法；
[0005] 2005年李全勝等人公開了一種基于模型道的地表一致性剩余靜校正方法，論述了 模型道的構(gòu)造方法，將剩余靜校正問題作為廣義線性反問題來求解炮點(diǎn)和檢波點(diǎn)剩余靜校 正量；
[0006] 2006年段云卿公開了一種折射波剩余靜校正方法，利用折射波信噪比高的特點(diǎn)， 將反射波剩余靜校正原理應(yīng)用于折射波處理，通過交互手段，逐段估算折射波的速度，用合 適的速度對地震記錄進(jìn)行線性動校正，在共炮點(diǎn)或共中心點(diǎn)道集上，用相關(guān)方法計(jì)算各道 與模型道時差，再用統(tǒng)計(jì)方法計(jì)算出各炮點(diǎn)及檢波點(diǎn)的剩余靜校正量；
[0007] 2011年潘樹林公開了一種共炮（檢）點(diǎn)剩余靜校正方法，基于經(jīng)過動校正后的共 炮點(diǎn)數(shù)據(jù)集和共檢波點(diǎn)數(shù)據(jù)集分別求取檢波點(diǎn)剩余靜校正量和炮點(diǎn)剩余靜校正量，并分 別將共炮、檢點(diǎn)道集動校正后的疊加道作為模型道與道集內(nèi)各個道進(jìn)行互相關(guān)，進(jìn)而求出 各個炮、檢波點(diǎn)的剩余靜校正量。
[0008] 針對三維觀測系統(tǒng)，上述剩余靜校正方法都依賴地震處理系統(tǒng)，即需要在對地震 資料進(jìn)行一定處理后才能應(yīng)用，所要計(jì)算的數(shù)據(jù)量大，因此對軟件和硬件都提出了較高要 求，如需要配置工作站和資料處理系統(tǒng)。由于資料處理中心一般擁有大、中型計(jì)算機(jī)，數(shù)據(jù) 處理能力強(qiáng)，故上述剩余靜校正方法一般可以在資料處理中心進(jìn)行。而一般情況下，野外現(xiàn) 場計(jì)算機(jī)的數(shù)據(jù)處理能力較弱，故針對三維觀測系統(tǒng)，上述剩余靜校正方法無法在野外現(xiàn) 場得到實(shí)際應(yīng)用。


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

[0009] 本申請實(shí)施例的目的在于提供一種三維剩余靜校正方法，以能夠在野外現(xiàn)場或資 料處理中心針對三維觀測系統(tǒng)進(jìn)行剩余靜校正。
[0010] 本申請實(shí)施例提供的一種三維剩余靜校正方法是這樣實(shí)現(xiàn)的：
[0011]-種三維剩余靜校正方法，包括：
[0012] 1)利用表層資料或者大炮初至資料建立表層模型，根據(jù)所述表層模型進(jìn)行基準(zhǔn)面 靜校正；
[0013] 2)對拾取的每個地震道的初至?xí)r間數(shù)據(jù)施加該地震道對應(yīng)的炮點(diǎn)和檢波點(diǎn)的基 準(zhǔn)面靜校正，生成新的地震道的初至?xí)r間Tu，其中，Ts，j表示第S炮點(diǎn)第j檢波點(diǎn)施加基準(zhǔn) 面靜校正后的地震道的初至?xí)r間；
[0014] 3)選取檢波點(diǎn)，確定每個檢波點(diǎn)所對應(yīng)的炮點(diǎn)線號范圍，根據(jù)每個檢波點(diǎn)所對應(yīng) 的炮點(diǎn)線號范圍內(nèi)的炮點(diǎn)生成該檢波點(diǎn)的共檢波點(diǎn)道集，確定每個共檢波點(diǎn)道集所對應(yīng)的 地震道的初至?xí)r間Ty，其中，Ty表示第r檢波點(diǎn)第i炮點(diǎn)地震道的初至?xí)r間；
[0015] 4)選取炮點(diǎn)，生成每個炮點(diǎn)的共炮點(diǎn)道集，確定每個共炮點(diǎn)道集折射層范圍內(nèi)各 檢波點(diǎn)處的初至擬合時差，并將該初至擬合時差作為各檢波點(diǎn)的單次剩余靜校正，其 中，Λ表示第s炮點(diǎn)第j檢波點(diǎn)的單次剩余靜校正；
[0016] 5)計(jì)算檢波點(diǎn)的剩余靜校正ARtj,其中，
[0017] Δ Rtj表示第j檢波點(diǎn)的剩余靜校正；

【權(quán)利要求】
1. 一種三維剩余靜校正方法，其特征在于，包括： 1) 利用表層資料或者大炮初至資料建立表層模型，根據(jù)所述表層模型進(jìn)行基準(zhǔn)面靜校 正； 2) 對拾取的每個地震道的初至?xí)r間數(shù)據(jù)施加該地震道對應(yīng)的炮點(diǎn)和檢波點(diǎn)的基準(zhǔn)面 靜校正，生成新的地震道的初至?xí)r間Ts，j，其中，Ts，j表示第s炮點(diǎn)第j檢波點(diǎn)施加基準(zhǔn)面靜 校正后的地震道的初至?xí)r間； 3) 選取檢波點(diǎn)，確定每個檢波點(diǎn)所對應(yīng)的炮點(diǎn)線號范圍，根據(jù)每個檢波點(diǎn)所對應(yīng)的炮 點(diǎn)線號范圍內(nèi)的炮點(diǎn)生成該檢波點(diǎn)的共檢波點(diǎn)道集，確定每個共檢波點(diǎn)道集所對應(yīng)的地震 道的初至?xí)r間TM，其中，Ty表示第r檢波點(diǎn)第i炮點(diǎn)地震道的初至?xí)r間； 4) 選取炮點(diǎn)，生成每個炮點(diǎn)的共炮點(diǎn)道集，確定每個共炮點(diǎn)道集折射層范圍內(nèi)各檢波 點(diǎn)處的初至擬合時差△ ，并將該初至擬合時差作為各檢波點(diǎn)的單次剩余靜校正，其中， Aty.表示第s炮點(diǎn)第j檢波點(diǎn)的單次剩余靜校正； 5) 計(jì)算檢波點(diǎn)的剩余靜校正ARtj，其中， ARtj表示第j檢波點(diǎn)的剩余靜校正；
n表示在共炮點(diǎn)道集上計(jì)算了該檢波點(diǎn)單次剩余靜校正的炮點(diǎn)個數(shù)； S表示三維測線內(nèi)炮點(diǎn)的位置編號，范圍為在該炮點(diǎn)的共炮點(diǎn)道集上計(jì)算了第j檢波 點(diǎn)的單次剩余靜校正的炮點(diǎn)個數(shù)n; j表示三維測線內(nèi)檢波點(diǎn)的位置編號，范圍為從第一條檢波點(diǎn)線的第一個檢波點(diǎn)到最 后一條檢波點(diǎn)線的最后一個檢波點(diǎn)； 6) 對于步驟3)中生成的共檢波點(diǎn)道集，確定每個共檢波點(diǎn)道集折射層范圍內(nèi)各炮點(diǎn) 處的初至擬合時差A(yù)ty，并將該初至擬合時差作為各炮點(diǎn)的單次剩余靜校正，其中，A 表示第r檢波點(diǎn)第i炮點(diǎn)的單次剩余靜校正； 7) 計(jì)算炮點(diǎn)的剩余靜校正量ASti，其中
ASt,表示第i炮點(diǎn)的剩余靜校正； m表示在共檢波點(diǎn)道集上計(jì)算了該炮點(diǎn)單次剩余靜校正的檢波點(diǎn)個數(shù)； i表示三維測線內(nèi)炮點(diǎn)的位置編號，范圍為第一條炮點(diǎn)線的第一個炮點(diǎn)到最后一條炮 點(diǎn)線的最后一個炮點(diǎn)； r表示三維測線內(nèi)檢波點(diǎn)的位置編號，范圍為在該檢波點(diǎn)的共檢波點(diǎn)道集上計(jì)算了第i炮點(diǎn)的單次剩余靜校正的檢波點(diǎn)個數(shù)m。
2. 如權(quán)利要求1所述的方法，其特征在于，步驟3)所述的選取檢波點(diǎn)，包括： 對于三維觀測系統(tǒng)，按照預(yù)設(shè)的檢波點(diǎn)線增量選取檢波點(diǎn)線，并在選取的每個檢波點(diǎn) 線上按照預(yù)設(shè)的檢波點(diǎn)增量選取檢波點(diǎn)。
3. 如權(quán)利要求2所述的方法，其特征在于，所述預(yù)設(shè)的檢波點(diǎn)線增量包括1-10中的任 一整數(shù)，所述預(yù)設(shè)的檢波點(diǎn)增量包括1-20中的任一整數(shù)。
4. 如權(quán)利要求1所述的方法，其特征在于，步驟3)所述的確定每個檢波點(diǎn)所對應(yīng)的炮 點(diǎn)線號范圍，根據(jù)每個檢波點(diǎn)所對應(yīng)的炮點(diǎn)線號范圍內(nèi)的炮點(diǎn)生成該檢波點(diǎn)的共檢波點(diǎn)道 集，包括： 確定每個檢波點(diǎn)所對應(yīng)的炮點(diǎn)線號范圍，所述炮點(diǎn)線號范圍的最小炮點(diǎn)線號為Slinel，所述炮點(diǎn)線號范圍的最大炮點(diǎn)線號為Sline2,其中，Slinel=Rline-S，Sline2 = Rline+S，Rline表示該檢波點(diǎn)的檢波點(diǎn)線號，S表示該檢波點(diǎn)所對應(yīng)的炮點(diǎn)線號偏移值； 根據(jù)所述炮點(diǎn)線號范圍內(nèi)的炮點(diǎn)生成該檢波點(diǎn)的共檢波點(diǎn)道集。
5. 如權(quán)利要求1所述的方法，其特征在于，步驟4)所述的選取炮點(diǎn)，包括：選取三維觀 測系統(tǒng)的全部炮點(diǎn)。
6. 如權(quán)利要求1所述的方法，其特征在于，步驟4)所述的選取炮點(diǎn)，包括： 對于三維觀測系統(tǒng)，按照預(yù)設(shè)的炮點(diǎn)線增量選取炮點(diǎn)線，并在選取的每個炮點(diǎn)線上按 照預(yù)設(shè)的炮點(diǎn)增量選取炮點(diǎn)。
7. 如權(quán)利要求6所述的方法，其特征在于，所述預(yù)設(shè)的炮點(diǎn)線增量包括1-10中的任一 整數(shù)，所述預(yù)設(shè)的炮點(diǎn)增量包括1-20中的任一整數(shù)。
8. 如權(quán)利要求1所述的方法，其特征在于，步驟4)所述的確定每個共炮點(diǎn)道集折射層 范圍內(nèi)各檢波點(diǎn)處的初至擬合時差A(yù)ty，包括： 將每個共炮點(diǎn)道集的數(shù)據(jù)按照檢波點(diǎn)線、正負(fù)偏移距進(jìn)行分類，每個分類分別生成一 個數(shù)據(jù)子集，在每個數(shù)據(jù)子集上根據(jù)折射層的變化劃分出一個折射層的范圍； 對該數(shù)據(jù)子集折射層范圍內(nèi)的各檢波點(diǎn)的偏移距與初至?xí)r間進(jìn)行擬合，得到每個炮點(diǎn) 數(shù)據(jù)子集指定折射層內(nèi)各檢波點(diǎn)的擬合初至?xí)r間，擬合前初至?xí)r間為Tsj 計(jì)算該炮點(diǎn)數(shù)據(jù)子集內(nèi)最小絕對值偏移距對應(yīng)的初至擬合時差，將最小絕對值偏 移距對應(yīng)的初至擬合時差作為該炮點(diǎn)數(shù)據(jù)子集指定折射層范圍內(nèi)的擬合系統(tǒng)校正時差 constSsk； 根據(jù)C〇nstSsk計(jì)算每個共炮點(diǎn)道集所有子集指定折射層范圍內(nèi)各檢波點(diǎn)處的初至擬 合時差
T' 表示第s炮點(diǎn)第j檢波點(diǎn)的擬合初至?xí)r間： Ts，j表示第s炮點(diǎn)第j檢波點(diǎn)施加基準(zhǔn)面靜校正后的地震道的初至?xí)r間。
9. 如權(quán)利要求1所述的方法，其特征在于，步驟6)所述的對于步驟3)中生成的共檢波 點(diǎn)道集，確定每個共檢波點(diǎn)道集折射層范圍內(nèi)各炮點(diǎn)處的初至擬合時差A(yù)ty，包括： 在由步驟3)生成的所有共檢波點(diǎn)道集上按照正負(fù)偏移距分別生成一個共檢波點(diǎn)數(shù)據(jù) 子集，在每個子集上根據(jù)折射速度的大小和范圍劃分出一個折射層的范圍； 對每個共檢波點(diǎn)數(shù)據(jù)子集上指定折射層范圍內(nèi)各炮點(diǎn)對應(yīng)的炮檢距與初至?xí)r間分別 進(jìn)行擬合，得到每個共檢波點(diǎn)道集數(shù)據(jù)子集內(nèi)各炮點(diǎn)處擬合后的初至?xí)r間擬合前初至 時間為I；，k; 計(jì)算該共檢波點(diǎn)道集數(shù)據(jù)子集最小絕對值正負(fù)偏移距對應(yīng)的初至擬合時差，將最小絕 對值正負(fù)偏移距對應(yīng)的初至擬合時差作為該共檢波點(diǎn)道集數(shù)據(jù)子集指定折射層范圍內(nèi)的 擬合系統(tǒng)校正時差constSA; 根據(jù)C〇nStSrt計(jì)算每個共檢波點(diǎn)道集所有子集指定折射層范圍內(nèi)各檢波點(diǎn)處的初至 擬合時差A(yù)tr，i;
T'u表示第r檢波點(diǎn)第i炮點(diǎn)的擬合初至?xí)r間： 表示第r檢波點(diǎn)第i炮點(diǎn)施加基準(zhǔn)面靜校正后的地震道的初至?xí)r間。
【文檔編號】G01V1/36GK104483705SQ201410852329
【公開日】2015年4月1日 申請日期:2014年12月31日 優(yōu)先權(quán)日:2014年12月31日
【發(fā)明者】馮發(fā)全 申請人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司