一種三維剩余靜校正方法
【專利摘要】本申請實施例公開了一種三維剩余靜校正方法，屬于地球物理勘探【技術(shù)領(lǐng)域】。該方法包括：將計算的基準(zhǔn)面靜校正施加到野外拾取的地震初至?xí)r間數(shù)據(jù)體上，得到一個新的地震初至?xí)r間數(shù)據(jù)體，然后在共炮點道集和共檢波點道集的正負(fù)偏移距上，分別選擇劃分折射層范圍，對該道集正負(fù)偏移距上的折射層的初至?xí)r間與偏移距進(jìn)行擬合，得到各地震道的擬合初至?xí)r間，然后求取二者的擬合時差，最后對所有檢波點和炮點的擬合時差各自進(jìn)行統(tǒng)計計算，其平均值即為各檢波點和炮點的剩余靜校正。本申請實施例的方法不依賴地震處理系統(tǒng)，計算的數(shù)據(jù)量小，對地震資料信噪比要求不高，便于在現(xiàn)場微機(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)混合反演剩余靜校正的方法，針對剩余靜 校正問題的非線性特點，分析了已有的線性反演及非線性反演方法的優(yōu)缺點，并在此基礎(chǔ) 上提出了一種混合反演方法；
[0005] 2005年李全勝等人公開了一種基于模型道的地表一致性剩余靜校正方法，論述了 模型道的構(gòu)造方法，將剩余靜校正問題作為廣義線性反問題來求解炮點和檢波點剩余靜校 正量；
[0006] 2006年段云卿公開了一種折射波剩余靜校正方法，利用折射波信噪比高的特點， 將反射波剩余靜校正原理應(yīng)用于折射波處理，通過交互手段，逐段估算折射波的速度，用合 適的速度對地震記錄進(jìn)行線性動校正，在共炮點或共中心點道集上，用相關(guān)方法計算各道 與模型道時差，再用統(tǒng)計方法計算出各炮點及檢波點的剩余靜校正量；
[0007] 2011年潘樹林公開了一種共炮（檢）點剩余靜校正方法，基于經(jīng)過動校正后的共 炮點數(shù)據(jù)集和共檢波點數(shù)據(jù)集分別求取檢波點剩余靜校正量和炮點剩余靜校正量，并分 別將共炮、檢點道集動校正后的疊加道作為模型道與道集內(nèi)各個道進(jìn)行互相關(guān)，進(jìn)而求出 各個炮、檢波點的剩余靜校正量。
[0008] 針對三維觀測系統(tǒng)，上述剩余靜校正方法都依賴地震處理系統(tǒng)，即需要在對地震 資料進(jìn)行一定處理后才能應(yīng)用，所要計算的數(shù)據(jù)量大，因此對軟件和硬件都提出了較高要 求，如需要配置工作站和資料處理系統(tǒng)。由于資料處理中心一般擁有大、中型計算機(jī)，數(shù)據(jù) 處理能力強(qiáng)，故上述剩余靜校正方法一般可以在資料處理中心進(jìn)行。而一般情況下，野外現(xiàn) 場計算機(jī)的數(shù)據(jù)處理能力較弱，故針對三維觀測系統(tǒng)，上述剩余靜校正方法無法在野外現(xiàn) 場得到實際應(yīng)用。


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

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

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