[發明專利]一種基于輔助坐標系的起伏地表波形反演方法

有效

• 申請號: 201610037642.8
• 申請日: 2016-01-20
• 公開（公告）號: CN105549080B
• 公開（公告）日: 2017-08-25
• 發明（設計）人: 曲英銘;李振春;李金麗;黃建平
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 濟南舜源專利事務所有限公司37205
• 代理人: 王連君

技術領域

本發明屬于石油地球物理勘探領域，具體涉及一種基于輔助坐標系的起伏地表波形反演方法。

背景技術

劇烈起伏地表對地震資料采集、處理帶來了嚴重的影響，地球物理研究人員發展了一系列方法來克服這個問題。目前，針對起伏地表的處理主要有兩種策略：一是對表層波場進行校正，二是基于起伏地表進行深度偏移成像。但復雜條件下靜校正量難以準確計算，且靜校正無法徹底消除地表起伏對地震波場造成的扭曲，因此直接基于起伏地表的深度偏移成像方法逐漸成為研究的熱點。

由于全波形反演高精度以及高分辨率的特點，使其成為速度建模的一種有力工具，逐漸成為研究的熱點。全波形反演是一個非線性數據擬合的過程，通過減少觀測數據與預測數據的之間的差值來更新參數模型，這個過程以迭代的方式重復下去，直到數據差值足夠小為止。

發明內容

針對現有技術中存在的上述技術問題，本發明提出了一種基于輔助坐標系的起伏地表波形反演方法，設計合理，具有良好的推廣價值。

為了實現上述目的，本發明采用如下技術方案：

一種基于輔助坐標系的起伏地表波形反演方法，按照如下步驟進行：

步驟1：輸入初始全局速度場、常規疊前炮記錄、起伏高程及震源子波，并建立觀測系統；

步驟2：根據初始全局速度場及起伏高程進行網格剖分，將近地表附近的網格剖分成曲網格，深層的網格剖分成矩形網格；

步驟3：將初始全局速度場變換到輔助坐標系下的矩形網格，采用下式所示的變換格式：

其中，x和z表示笛卡爾坐標系下的橫縱坐標；ξ和η表示輔助坐標系下的橫縱坐標；z_i-1(ξ)和z_i(ξ)是笛卡爾坐標系下第i層頂界面、底界面的高程，定義最深層高程為零；η_i-1(ξ)和η_i(ξ)是輔助坐標系第i層頂界面、底界面的縱向采樣點數，定義最深層縱向采樣點數為零；

步驟4：對常規疊前炮記錄劃分時窗，在輔助坐標系下對近地表附近的曲網格區域速度場應用早至波波形反演，更新近地表的速度場，早至波波形反演的梯度方向如下：

其中，g為梯度；x_s為炮點坐標；v為介質速度；p為聲壓；p^*是殘差波場的反向傳播；t為時間；T_e表示早至波時窗；

步驟5：判斷使用近地表速度場正演模擬的炮記錄與步驟4中劃分時窗后的常規疊前炮記錄之差是否滿足誤差條件；

若：判斷結果是使用近地表速度場正演模擬的炮記錄與步驟4中劃分時窗后的常規疊前炮記錄之差滿足誤差條件，則近地表速度場更新完成，然后執行步驟6；

或判斷結果是使用近地表速度場正演模擬的炮記錄與步驟4中劃分時窗后的常規疊前炮記錄之差不滿足誤差條件，則執行步驟4；

步驟6：將步驟1中輸入的常規疊前炮記錄分解成不同主頻的多尺度炮記錄，采用如下式所示的分解公式：

其中，F_w是維納濾波器；W_o表示初始炮記錄子波；W_t是生成的炮記錄子波；ω是角頻率；ε為一個小數；*表示共軛轉置；

步驟7：將更新過近地表速度的全局速度場作為初始速度場，在輔助坐標系下應用全波形反演從低頻到高頻更新全局速度場，全波形反演的梯度方向如下：

其中，f表示主頻；f₁和f_max為多尺度分解的最低頻率和最高頻率；T_max表示常規疊前炮記錄的最大記錄時間；表示主頻是f的殘差波場反傳；

步驟8：判斷使用全局速度場正演模擬的炮記錄與常規疊前炮記錄之差是否滿足誤差條件；

若：判斷結果是使用全局速度場正演模擬的炮記錄與常規疊前炮記錄之差滿足誤差條件，則全局速度場更新完成，然后執行步驟9；

或判斷結果是使用全局速度場正演模擬的炮記錄與常規疊前炮記錄之差不滿足誤差條件，則執行步驟7；

步驟9：將更新完成的全局速度場反變換到笛卡爾坐標系下，反變換公式如下：

步驟10：輸出反演的速度場。

優選地，在步驟4中，具體包括

步驟4.1：定義目標函數：