共成像點反射角角道集的生成方法及裝置的制造方法
【技術領域】
[0001] 本發明涉及地震勘探中反射地震信息處理技術領域，具體而言，涉及一種共成像 點反射角角道集的生成方法及裝置。
【背景技術】
[0002] 振幅隨偏移距變化（Amplitude Versus Offset，簡稱為AV0)技術廣泛應用于儲層 反演、預測巖石物性，有力地提高了儲層預測的準確性。與疊后反演相比，疊前地震反演可 以提供更豐富的反演信息，除縱波阻抗外，可以獲得橫波阻抗、總橫波密度和密度等，由此 可以計算各種彈性參數，為巖性和流體識別與預測提供豐富的信息。
[0003] AVO技術的理論基礎是振幅隨入射角變化（Amplitude Versus Angle，簡稱為 AVA)，二者在地下是水平層狀介質才具有簡單的幾何影射關系，可以由偏移距經過簡單的 幾何影射關系獲得入射角。國內外工業界俗稱的角道集是把正常時間校正（Normal Move Out，簡稱為ΝΜ0)校正后共中心點道集或者疊前偏移后偏移距道集影射到角度域。對于地 下復雜地質體，炮檢距向入射角采用上述簡單轉化關系會帶來入射角與振幅的誤差，從而 產生AVA假象。因此需要發展直接求取地下成像點處入射角的計算方法，進而獲得角道集， 其核心是獲得"真振幅"的成像結果，才能很好的反映地下的物性參數。
[0004] 疊前時間偏移方法（Pre-Stack Time Migration，簡稱為PSTM)可以對一類傾角、 斷層較為復雜但速度橫向變化不是很劇烈的構造較好成像，具有較高的計算效率。只是需 要疊加速度，即可簡單地通過速度掃描等方式得到恰當的速度模型，回避了疊前深度偏移 方法速度建模困難的問題，因此，采用疊前時間偏移技術輸出的共成像點反射角角道集，是 用于AVA分析較好的選擇。
[0005] 在相關技術中，疊前時間偏移生成角道集時，通常采用疊加速度或者均方根速度 來進行入射角計算，其小于層速度，因此計算的入射角小于實際角度，從而生成共成像點反 射角角道集不準確。
[0006] 針對相關技術中疊前時間偏移方法中生成共成像點反射角角道集不準確的問題， 目前尚未提出有效的解決方案。
【發明內容】
[0007] 本發明的主要目的在于提供一種共成像點反射角角道集的生成方法及裝置，以解 決相關技術中疊前時間偏移方法中生成共成像點反射角角道集不準確的問題。
[0008] 為了實現上述目的，根據本發明的一個方面，提供了一種共成像點反射角角道集 的生成方法。該生成方法通過疊前時間偏移方法生成共成像點反射角角道集，該生成方法 包括：獲取層速度，其中，層速度為地層中地震波傳播的速度；根據層速度計算入射角，其 中，入射角為地震射線的入射角；以及根據入射角生成共成像點反射角角道集。
[0009] 進一步地，根據層速度計算入射角包括：獲取均方根速度、走時和半偏移距，其中， 均方根速度為把地層中地震波的高次曲線看作二次曲線計算出地震波的速度，走時是地層 中地震波從震源傳到觀測點所經過的時間，半偏移距為一半偏移距的距離；以及根據層速 度、均方根速度、走時和半偏移距計算入射角。
[0010] 進一步地，根據層速度、均方根速度、走時和半偏移距計算入射角包括：根據第一 關系式計算入射角，其中，第一關系式為：
[0011]
[0012] ??η? rms
[0013] V是層速度，均方根速度，γ是入射角，T是走時，h是半偏移距，X s是炮點 橫坐標，ys是炮點縱坐標，X 1^是接收點橫坐標，y 1^是接收點縱坐標，X是成像點橫坐標，y是 成像點縱坐標。
[0014] 進一步地，在根據入射角生成共成像點反射角角道集之后，該生成方法還包括：對 共成像點反射角角道集進行時間-角度域剩余動校正；去除校正后的共成像點反射角角道 集中拉伸和噪音部分；以及根據去除后的共成像點反射角角道集形成用于疊前反演的角道 集。
[0015] 進一步地，在對切除后的共成像點反射角角道集應用常規噪音壓制技術形成用于 疊前反演的角道集之后，該生成方法還包括：將用于疊前反演的角道集中不同反射角進行 疊加，形成偏移疊加剖面；將偏移疊加剖面轉換為地下反射構造的剖面圖像；以及顯示剖 面圖像。
[0016] 進一步地，該第一關系式由以下方法得到：
[0017] rms \ rms： v .5
^
[0018] 其中，？：!是χ方向的射線參數，ρ ,是y方向的射線參數，χ。是成像點橫坐標，y。是 成像點縱坐標，Θ z是入射波與z軸的夾角；
[0019] 通過關系式（1)得到：
[0026] 其中
[0027] 根據關系式（I)、關系式（2)、關系式（3)、關系式（4)和關系式（5)得到第一關系 式。
[0028] 為了實現上述目的，根據本發明的另一方面，提供了一種共成像點反射角角道集 的生成裝置。該生成裝置通過疊前時間偏移裝置生成共成像點反射角角道集，該生成裝置 包括：獲取單元，用于獲取層速度，其中，層速度為地層中地震波傳播的速度；計算單元，用 于根據層速度計算入射角，其中，入射角為地震射線的入射角；以及生成單元，用于根據入 射角生成共成像點反射角角道集。
[0029] 進一步地，獲取單元包括：獲取模塊，用于獲取均方根速度、走時和半偏移距，其 中，均方根速度為把地層中地震波的高次曲線看作二次曲線計算出地震波的速度，走時是 地層中地震波從震源傳到觀測點所經過的時間，半偏移距為一半偏移距的距離；以及計算 模塊，用于根據層速度、均方根速度、走時和半偏移距計算入射角。
[0030] 進一步地，該生成裝置還包括：校正單元，用于對共成像點反射角角道集進行時 間-角度域剩余動校正；去除單元，用于去除校正后的共成像點反射角角道集中拉伸和噪 音部分；以及形成單元，用于根據去除后的共成像點反射角角道集形成用于疊前反演的角 道集。
[0031] 進一步地，該生成裝置還包括：疊加單元，用于將用于疊前反演的角道集中不同反 射角進行疊加，形成偏移疊加剖面；轉換單元，用于將偏移疊加剖面轉換為地下反射構造的 剖面圖像；以及顯示單元，用于顯示剖面圖像。
[0032] 通過本發明，由于采用層速度計算入射角，提升了入射角的精度，根據計算得到的 入射角生成共成像點反射角角道集，從而提升了生成共成像點反射角角道集的精度，解決 了相關技術中疊前時間偏移方法中生成共成像點反射角角道集不準確的問題，提高了生成 共反射點入射角角道集的精度，進而達到提升儲層預測的準確性。
【附圖說明】
[0033] 構成本申請的一部分的附圖用來提供對本發明的進一步理解，本發明的示意性實 施例及其說明用于解釋本發明，并不構成對本發明的不當限定。在附圖中：
[0034] 圖1是根據本發明實施例的共成像點反射角角道集的生成方法的流程圖；
[0035] 圖2是根據本發明實施例的共成像點反射角角道集的生成方法中偏移成像結果 的不意圖；
[0036] 圖3是根據本發明實施例的共成像點反射角角道集的生成方法中抽取的7-14Km 處的入射角10度偏移疊加剖面的示意圖；
[0037] 圖4是根據本發明實施例的共成像點反射角角道集的生成方法中IOKm處的偏移 距道集的示意圖；
[0038] 圖5是根據本發明實施例的共成像點反射角角道集的生成方法中IOKm處的角道 集的不意圖；
[0039] 圖6是根據本發明實施例的共成像點反射角角道集的生成方法中IOKm處的角道 集的局部放大的示意圖；
[0040] 圖7是根據本發明實施例的共成像點反射角角道集的生成方法中角道集中幅值 隨角度變化的曲線；以及
[0041] 圖8是根據本發明實施例的共成像點反射角角道集的生成裝置的示意圖。
【具體實施方式】
[0042] 需要說明的是，在不沖突的情況下，本申請中的實施例及實施例中的特征可以相 互組合。下面將參考附圖并結合實施例來詳細說明本發明。
[0043] 為了使本技術領域的人員更好地理解本申請方案，下面將結合本申請實施例中的 附圖，對本申請實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是 本申請一部分的實施例，而不是全部的實施例。基于本申請中的實施例，本領域普通技術 人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都應當屬于本申請保護的范 圍。
[0044] 需要說明的是，本申請的說明書和權利要求書及上述附圖中的術語"第一"、"第 二"等是用于區別類似的對象，而不必用于描述特定的順序或先后次序。應該理解這樣使 用的數