本發明公開了一種基于進退法的橫波分裂分析方法，主要步驟：首先、根據快慢橫波相關性構造目標函數；然后，通過確定的真實值(或實際數據)得到R分量數據和T分量數據，進一步計算得到S1(t)和S2(t)；根據區間大小，將區間劃分為4個子區間，取每個子區間的中心點作為初始點；在每個區間采用進退法確定該區間的結果值，得到4個結果值，比較4個結果值，取最大COV值對應的時差和角度作為輸出值。本發明以快慢波的相關性作為目標函數，通過縮小極值區間來確定極值點，并且通過分區的方式獲得真實的裂縫參數，為非常規油氣儲層提供準確的裂縫信息。

技術領域

本發明涉及多分量地震資料處理技術領域，尤其是一種針對裂縫型儲層預測的橫波分裂分析方法。

背景技術

裂縫作為檢測非常規油氣含量的重要參數，它不僅可以連接孤立的孔隙，增加儲層的有效孔隙度，還可以作為重要的油氣運移空間，提高儲層的滲透率。20世紀90年代開始，隨著多波多分量勘探技術開始登上舞臺，橫波分裂技術開始應用轉換橫波，近年來，隨著3D3C多波多分量勘探技術的飛速發展，橫波分裂技術的發展又到了一個新的高潮。當橫波的偏振方向以斜交地方式穿過含裂縫的地層時，會發生橫波分裂現象(SWS)。初始橫波分裂為一個快橫波和一個慢橫波，分裂的快橫波偏振方向與裂縫走向平行，快慢波時差與裂縫的強度密切相關，因此，橫波分裂裂縫檢測技術成為研究裂縫方向和其發育程度最直接最可靠的方法之一。當今常規的橫波分裂分析方法如旋轉相關法、切向能量法、協方差矩陣法都是一種枚舉式的網格掃描算法。旋轉相關法通過快慢橫波相關性能很好的尋找最佳的裂縫參數(快橫波偏振方向和快慢橫波時差)。但是，在面對龐大的三分量采集數據時，枚舉式的計算特點必定會造成計算效率低下。所以，如何快速地提高計算效率，得到準確的裂縫參數信息，對裂縫型油氣藏儲層預測十分重要。

發明內容

針對現有的橫波分裂分析方法存在的面對龐大的三分量采集數據時計算效率低下的問題，本發明提供一種計算效率快的新的基于進退法的橫波分裂分析方法，該方法以快慢波的相關性作為目標函數，通過縮小極值區間來確定極值點，快速而準確的找到真實的裂縫參數。從而，得到計算速率最快，精度最高的橫波分裂分析方法，為非常規油氣儲層提供準確的裂縫信息。

本發明提供的基于進退法的橫波分裂分析方法，步驟如下：

S1、根據快慢橫波相關性構造目標函數如下：

式中，θ和Δt為SWS(橫波分裂)參數，S1(t)和S2(t)分別為一個時窗內的快橫波和慢橫波的時間序列；和分別是S1(t)和S2(t)在一個時窗內的平均值。

S2、通過確定的真實值或實際數據得到此時的R分量數據R(t)和T分量數據T(t)，由R(t)和T(t)根據以下公式計算得到S1(t)和S2(t)：

S3、將測試區間劃分為4個子區間，取每個子區間的中心點作為初始點。例如總區間為180×60，180為方位角范圍，60為時差范圍，分為四個子區間后，每個子區域的大小為90×30。

S4、根據給定的初始點確定區間，采用二維進退法確定該區間的結果值。具體方法步驟如下：

(1)對Δt方向零點的搜索；