剖面上存在明顯的地震放射同相軸下拉，河床底部與下覆地層反射 界面強。本發(fā)明提供的古河道河床侵蝕深度自動識別方法，對種子點和每個外推地震響應(yīng) 點，均采用9個正方形平面窗口來計算種子點和外推地震響應(yīng)點所在的反射層面的相干和 傾角（如圖2所示），每個正方形平面窗口代表一個種子點或外推地震響應(yīng)點的地質(zhì)狀況， 每個正方形平面窗口對應(yīng)圖5、圖6和圖7中的一個圓點。根據(jù)相干和傾角，驅(qū)動種子點和 每個外推地震響應(yīng)點向河道工區(qū)范圍內(nèi)的測線上，自動追蹤相位相同的外推地震響應(yīng)點， 在自動追蹤的過程中，計算每個外推地震響應(yīng)點處的相干平均值和方差值即為計算9個正 方形平面窗口的相干平均值和方差值，通過相干平均值和方差值判斷河道邊界，當(dāng)出現(xiàn)相 干平均值小，方差值大時，證明外推到了河床邊界，直至得到所有測線上的種子點，實現(xiàn)全 區(qū)河床地震響應(yīng)同相軸自動識別，最終實現(xiàn)多條多期河道的侵蝕深度預(yù)測。具體如下（下 文中所述的"外推"的含義均為自動追蹤相位相同的外推地震響應(yīng)點）：
[0062] 圖3是本發(fā)明實施例中古河道河床侵蝕深度自動識別方法的具體流程示意圖， 圖4是本發(fā)明實施例中在河床地震剖面反射波的最大波峰處設(shè)置種子點示意圖；圖5是 本發(fā)明實施例中基于河床種子點自動追蹤軌跡平面示意圖；如圖3、圖4和圖5所示，在 河床地震剖面反射波的最大波峰處選取一個地震響應(yīng)點作為種子點（如圖4所示），以 種子點U(i，為中心，其中，u為河床底界的地震響應(yīng)點，i為inline(主測線）， j為crossline (聯(lián)絡(luò)測線），為時間，采用9個包含9個采樣點的正方形平面窗口估 算種子點所在的反射層面的相干和傾角（如上段所述）。估算9個窗口的相干及視傾角 ck (Pk, Qk)，其中k = 1，2, ...，9。通過傾角驅(qū)動種子點外推，實現(xiàn)河床底部地震同相軸的自 動追蹤。
[0063] 如圖2所示，種子點具體向各測線方向外推（自動追蹤）的方法為：估算窗口 1、2、4、5、7、8的相干值，取其中最大的相干值及其對應(yīng)的視傾角，向位于同一 inline線 crossline減少的方向驅(qū)動河床底界層位自動追蹤，估算窗口 2、3、5、6、8、9相干值，取其中 最大的相干值及其對應(yīng)的視傾角，向同一 inline線crossline增大的方向驅(qū)動河床底界層 位自動追蹤，估算窗口 1、2、3、4、5、6相干值，取其中最大的相干值及其對應(yīng)的視傾角，向同 一 crossline道inline增大的方向驅(qū)動河床底界層位自動追蹤，估算窗口 4、5、6、7、8、9相 干值，取其中最大的相干值及其對應(yīng)的視傾角，向同一 crossline道inline減少的方向驅(qū) 動河床底界層位自動追蹤。每個外推地震響應(yīng)點具體向各測線方向外推時，也是如上所述， 在此不再贅述。
[0064] 同時結(jié)合多窗口旋轉(zhuǎn)掃描識別河床邊界，計算9個窗口的相干值的方差和平均值 分別為冰=Ilk-。在遠離河道側(cè)邊界，地震振幅、波形相似，相干平 均值大，方差值小。在靠近河道側(cè)邊界處的地震波形呈雜亂反射，其相似程度低，相干平均 值小，方差值大。因此，可以通過9個窗口相干值的方差及平均值的變化識別河道邊界。
[0065] 具體地，如圖3、圖4和圖5所示，首先，在河床地震反射同相軸中間部位的最大波 峰處選取一個采樣點作為種子點（圖4中圓點所示）：u (i，j，tu)，其中，u為河床底界的采 樣點，i為主測線inline，j為聯(lián)絡(luò)測線crossline，tu為河床位于線道（i, j)處的時間。 圖5中實心圓點為種子點在平面投影。
[0066] 然后，如圖3和5所示，在主測線Li線上向聯(lián)絡(luò)測線crossline增加和減少的兩個 方向外推（圖5中Li線上的橫向的箭頭所示），同時在每個外推地震響應(yīng)點進行旋轉(zhuǎn)窗口 掃描識別是否已外推到河道邊界處。在河床中間位置處，9個窗口的平均相干平均值大，方 差值小，種子點和鄰近的外推地震響應(yīng)點都位于同一河床上，平均相干值大，方差值小。在 兩側(cè)外推的過程中出現(xiàn)相干平均值小及方差值大的異常，表明已經(jīng)外推到了河床兩側(cè)邊界 處u(i, j-m tuO、u(i, j+n, ti，j+n)，此時，停止兩側(cè)的外推。
[0067] 將估算的u(i, j, tu) (j = j-m, j-m+1,--，j+n_l, j+n)作為新的種子點，向同 一聯(lián)絡(luò)測線crossline道，主測線inline增大和減少的方向外推（圖5中位于聯(lián)絡(luò)測線 Lj等多條聯(lián)絡(luò)測線的縱向箭頭所示）。以u(i, j, tu) (j = j-m, j-m+1,--,j+n-1, j+n) 外推鄰近線上的河床采樣點的過程中，首先從Li線上河床兩側(cè)邊界處u(i，j-m，ti+J、 u(i, j+n, ti，j+n)開始外推鄰近inline線上河床響應(yīng)點（即外推地震響應(yīng)點） u(i±l, j-m, ti±1，j_m)、u(i±l, j+n, ti±1，j+n)，對每個外推地震響應(yīng)點進行旋轉(zhuǎn)掃描進行河床 邊界識別。
[0068] 如圖6所示，如果以邊界點11(1，」_111，\」_"1)為種子點（圖6中的實心圓點）的 外推地震響應(yīng)點（圖6中位于聯(lián)絡(luò)測線Lj-m上實心圓點的上方、縱向箭頭所指空心圓） u(i+l, j-m, ti+1，j_m)處于河床內(nèi)部，則再以u(i+l, j-m, ti+1，j_m)為種子點，向同一 Li+Ι線 crossline減少的方向外推（圖6中Li+1線上橫向箭頭方向），一直外推到河道邊界終 止。而Li線上的非邊界點u(i, j-m+k, ，（k = 1，2, 3,...）（例如圖6中主測線Li 上的、位于聯(lián)絡(luò)測線Lj-m+l，Lj-m+2，Lj-m+3上的空心圓點），則根據(jù)已經(jīng)在Li線上外推得 到u(i，j-m+k，(k = 1，2, 3...)作為種子點，外推得到Li+Ι線上的外推地震響應(yīng)點 (圖6中Lj-m+1、Lj-m+2和Lj-m+3線上箭頭所指向的圓點）。
[0069] 如圖7所示，如果以邊界點u(i，為種子點（圖7中實心圓點）外推 (圖7中Lj-m線上箭頭方向）鄰線的河床響應(yīng)點（即外推地震響應(yīng)點）u(i+l, j-m, ti+1，j_m) 位于河床邊界及外部，則視u(i, j-m, tui)為無效種子點，而以u(i, j-m+k, I^jfk)，k = 1，2, 3, · ·，p-1, p, p+L ·.為種子點外推（圖7中Lj-m+l、Lj-m+2和Lj-m+3線上箭頭方向） 對應(yīng)的外推地震響應(yīng)點 u(i+l, j-m+k, ti+1，j_m+k)，k = I, 2, 3,. . , p-1, p, p+1...。直到當(dāng)外 推地震響應(yīng)u(i+l，j-m+p，ti+1，j_m+p)位于河道內(nèi)部，則將u(i，j-m+p，作為有效的種 子點外推其對應(yīng)的u(i+l, j-m+p, ti+1，j_m+p)。按上述同樣方法進行邊界點u(i, j-m, tui)、 u (i, j+n, ti>j+n)向鄰線進行逐次外推，最終實現(xiàn)全區(qū)河床地震響應(yīng)同相軸自動識別。
[0070] 本發(fā)明技術(shù)方案能夠自動識別古河床侵蝕深度，通過古河床的識別，間接指導(dǎo)暗 河巖溶儲層的預(yù)測。同時，古河床深度的識別有利于預(yù)測高效油藏的位置。本發(fā)明提出了 基于地震波相似性原理的古河道河床侵蝕深度自動識別方法。該方法在設(shè)置河床種子點的 基礎(chǔ)上，以地層傾角和地震波形相似性驅(qū)動種子點外推，同時采用多窗口旋轉(zhuǎn)掃描識別河 床邊界，河床邊界約束外推范圍。最終實現(xiàn)多條多期河道的侵蝕深度預(yù)測。識別深切河床 的空間位置可以間接指導(dǎo)暗河巖溶儲層的刻畫。
[0071] 另外，在哈拉哈塘