[發明專利]基于模擬退火和多尺度脊波分析的縱橫波相關方法

在審

• 申請號: 201810140245.2
• 申請日: 2018-02-11
• 公開（公告）號: CN108363099A
• 公開（公告）日: 2018-08-03
• 發明（設計）人: 韓宏偉;程遠鋒;王興謀;羅霞;曲志鵬;曹麗萍;雷蕾;伊慧
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 濟南日新專利代理事務所 37224
• 代理人: 崔曉艷

本發明提供一種基于模擬退火和多尺度脊波分析的縱橫波相關方法，包括：步驟1，確定縱橫波原始剖面；步驟2，對縱波剖面進行多尺度脊波分析；步驟3，對橫波剖面進行多尺度脊波分析；步驟4，對多尺度脊波分析后的大尺度縱波剖面和橫波剖面進行模擬退火相關；步驟5，對多尺度脊波分析后的中尺度縱波剖面和橫波剖面進行模擬退火相關匹配；步驟6，對多尺度脊波分析后的小尺度縱波剖面和橫波剖面進行模擬退火相關。該基于模擬退火和多尺度脊波分析的縱橫波相關方法借助了多尺度脊波分析的優點和模擬退火算法對全局最優解的搜索能力，提高了縱橫波自動相關的精度，提高了對噪音影響的抵抗能力。

技術領域

本發明涉及地震資料處理領域，特別是涉及到一種基于模擬退火和多尺度脊波分析的縱橫波相關方法。

背景技術

近年來發展的多波地震勘探技術，充分利用了縱波和橫波在含油氣介質中傳播的不同特征及其豐富的波場信息，在解釋含油氣方面凸顯了良好的應用前景。多波勘探的最大的優點就是同時利用縱波和橫波的波場信息進行綜合解釋。其中解釋任務一個必不可少的環節就是要進行縱橫波的層位相關，然后應用縱波和橫波剖面的異同來進行地下未知巖層結構的解釋。縱觀國內外縱橫波相關技術的發展，大多數方法都是基于縱波和橫波同一地層反射同相軸相似性最大的原則來求取縱橫波的速度比，進而實現縱橫波在時間域的相關，其相關程度因而比較粗略。為此我們發明了一種新的基于模擬退火和多尺度脊波分析的縱橫波相關方法，解決了以上技術問題。

發明內容

本發明的目的是提供一種可提高縱橫波時間域精確自動相關的基于模擬退火和多尺度脊波分析的縱橫波相關方法。

本發明的目的可通過如下技術措施來實現：基于模擬退火和多尺度脊波分析的縱橫波相關方法，該基于模擬退火和多尺度脊波分析的縱橫波相關方法包括：步驟1，確定縱橫波原始剖面；步驟2，對縱波剖面進行多尺度脊波分析；步驟3，對橫波剖面進行多尺度脊波分析；步驟4，對多尺度脊波分析后的大尺度縱波剖面和橫波剖面進行模擬退火相關；步驟5，對多尺度脊波分析后的中尺度縱波剖面和橫波剖面進行模擬退火相關匹配；步驟6，對多尺度脊波分析后的小尺度縱波剖面和橫波剖面進行模擬退火相關。

本發明的目的還可通過如下技術措施來實現：

在步驟1中，取得多波多分量地震勘探資料中的原始縱波數據體和原始橫波數據體。

在步驟2中，對縱波剖面進行多尺度脊波分析，得到大尺度、中尺度和小尺度縱波剖面；在脊波分析中，首先對縱波剖面進行Radon變換，把圖像中的一維奇異性信號轉換成Radon域的一個點，然后用一維小波進行點的奇異性檢測。

在步驟2中，在將縱波剖面利用多尺度脊波分析分解成不同尺度的剖面步驟中，需要對圖像進行分塊，使每個分塊中的線條都近似于直線，然后再對每個分塊進行多尺度脊波分析，即首先對圖像進行子帶分解，然后將不同尺度的子帶圖像分成大小不同的子塊，分別對每個子塊進行多尺度脊波分析。

在步驟2中，定義雙變量函數f(x,y)的脊波變換為其中為脊波變換，a，b和θ分別表示尺度、位移和方向向量；通過Radon變換，把脊波變換與小波變換聯系起來，對于一個雙變量函數f(x,y)，其Radon變換為：式中R_f(a,b,θ)表示f(x,y)的Radon變換，而脊波變換系數可以通過對其Radon變換系數進行小波變換來計算，即式中t為積分變量；上式表明，脊波變換是對Radon變換切片的一維小波分析，其中方向角θ是固定的，而變量t是小波分析的對象。

在步驟3中，按照步驟2中的脊波變換方法，對橫波剖面進行多尺度脊波分析，得到大尺度、中尺度和小尺度橫波波剖面。

在步驟4中，首先設定初始解和初始控制參數，包括溫度、迭代次數以及終止條件，然后隨機產生新的解并計算轉移概率，判斷接受或放棄新解，最后判斷終止條件，滿足條件則停止迭代，否則重新隨機產生新解。

在步驟4中，模擬退火的求解步驟為：