[發(fā)明專利]一種地震數(shù)據(jù)快速相位校正處理的方法

有效

• 申請?zhí)枺?/td>: 201210016370.5
• 申請日: 2012-01-18
• 公開（公告）號: CN103217716A
• 公開（公告）日: 2013-07-24
• 發(fā)明（設計）人: 高少武;趙波;祝樹云;祝寬海
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;任潔

技術領域

本發(fā)明涉及油田的勘探、開發(fā)、開采技術，具體是為反映地下地層層位、油藏描述提供高分辨率的地震圖形和數(shù)據(jù)的一種地震數(shù)據(jù)快速相位校正處理的方法，特別適用于地震數(shù)據(jù)處理反褶積之后的實際地震數(shù)據(jù)，通過剩余子波相位校正處理，達到消除地震數(shù)據(jù)剩余子波相位影響、提高地震數(shù)據(jù)分辨率目的。

背景技術

隨著地震勘探技術的發(fā)展，油氣勘探的難度和深度也越來越大，對地震資料分辨率的要求也越來越高。反褶積處理是提高地震數(shù)據(jù)分辨率最有效的途徑。地震褶積模型是地震勘探數(shù)據(jù)處理中最基本的模型之一。常規(guī)的地震數(shù)據(jù)反褶積處理都是基于地震褶積模型。反褶積處理的一個根本假設就是地震子波是最小相位的。這個假設使得我們可以由子波振幅譜得到子波的相位譜并繼而得到完整的子波。而實際數(shù)據(jù)的地震子波是混合相位的，這便使得以最小相位子波假設為前提的反褶積處理后地震子波不是一個脈沖，而存在剩余子波。實際地震數(shù)據(jù)處理時，即使將混合相位地震子波轉換為最小相位地震子波，由于子波有限長度，反褶積處理后地震子波也還會存在剩余子波。地震數(shù)據(jù)中剩余子波的存在，降低了地震數(shù)據(jù)的分辨率。

在地震數(shù)據(jù)子波中，零相位子波具有最大分辨率，因此實際數(shù)據(jù)處理為了提高分辨率，希望處理后的地震子波是零相位的。“相同頻帶范圍的子波中，以零相位子波的分辨率為最高”(李慶忠，走向精確勘探的道路，P14，石油工業(yè)出版社，1994)。“具有相同振幅譜的諸子波中，零相位子波的分辨率最高”(俞壽朋，高分辨率地震勘探，P17，石油工業(yè)出版社，1993)。“子波相位校正的目的是使子波零相位化”(俞壽朋，俞壽朋文集，P206，石油工業(yè)出版社，2001)。“子波波形由振幅譜和相位譜二者決定。振幅很小的頻率成分對子波波形的影響很小，也就是說這些頻率成分的相位譜重要性不大，有些誤差也關系不大。而振幅譜比較大的頻率成分是起主要作用的，這些頻率成分的相位譜是重要的。但這個頻帶的相位譜一般很接近于直線。通過相位譜坐標原點做一條平行于主頻帶相位譜的直線，則主頻帶相位譜與此直線的差近于常數(shù)。因此可近似認為子波是常相位的。”(俞壽朋，高分辨率地震勘探，P167，石油工業(yè)出版社，1993)。

常規(guī)的消除地震數(shù)據(jù)剩余子波相位處理，采用相位校正方法。相位校正就是將剩余子波相位作為一個常數(shù)，顯然這樣的假設僅僅是一階近似。當子波的振幅譜較窄時，此法會有一定效果。在相位校正處理方法中，通常相位校正因子的估算采用對目標函數(shù)方差模進行掃描的方法來確定(Levy?S.and?Oldenburg?D.W.，Automatic?phase?correction?of?common-midpoint?stacked?data，Geophysics，1987，VOL.52，NO.1，P51-P59；周興元，常相位校正，石油地球物理勘探，1989，VOL.24，NO.2；李合群、周興元，時差、常相位校正及加權疊加，石油地球物理勘探，2000，VOL.35，NO.4，P415-418)。掃描方法估算相位校正因子非常費時，計算效率低。

發(fā)明內容

本發(fā)明目的在于提供一種不需要進行掃描，快速消除地震數(shù)據(jù)剩余子波相位影響，提高地震數(shù)據(jù)分辨率的地震數(shù)據(jù)快速相位校正處理的方法。

本發(fā)明采用的技術方案，包括以下步驟：

1)用人工震源激發(fā)和采集地震數(shù)據(jù)并做預處理；

步驟1)所述的預處理是指對地震數(shù)據(jù)置標簽、定義觀測系統(tǒng)、速度分析、動校正、反褶積、疊加處理。

2)計算凱澤窗希爾伯特因子序列；

步驟2)所述的計算凱澤窗希爾伯特因子序列包括計算第一類零階修正貝賽爾函數(shù)和凱澤窗希爾伯特因子序列：

按照以下公式計算第一類零階修正貝賽爾函數(shù)I₀(x)：

I₀(x)＝a₀+a₁+a₂+...+a_k???????(1)

式中，I₀(x)為第一類零階修正貝賽爾函數(shù)；x為第一類零階修正貝賽爾函數(shù)自變量，其取值范圍為0≤x≤1；a₀是第一類零階修正貝賽爾函數(shù)初始值，且a₀＝1；k表示第一類零階修正貝賽爾函數(shù)的第k項；a_k是第一類零階修正貝賽爾函數(shù)第k項值；

按照以下公式計算第一類零階修正貝賽爾函數(shù)第k項值a_k：