本發(fā)明涉及一種成像方法及系統(tǒng)，屬于地理探測，具體是涉及一種地震反射探測時間剖面成像方法及系統(tǒng)。

背景技術：

1、對地震探測來說，數(shù)據(jù)處理是其中一個關鍵環(huán)節(jié)，處理參數(shù)和處理流程選取的優(yōu)劣直接決定了地震數(shù)據(jù)處理的成敗，選取合理的地震數(shù)據(jù)處理參數(shù)對深地震反射探測高信噪比地震時間剖面的獲取以及后續(xù)地震地質(zhì)解譯具有重要作用?，F(xiàn)有的深地震反射探測地震數(shù)據(jù)處理流程通常是基于成熟的石油地震勘探數(shù)據(jù)處理流程，但是深地震反射探測探測深度更大，深部區(qū)域地層介質(zhì)為變質(zhì)巖和巖漿巖，地震反射波信號多以弱反射為主；另外受復雜地形影響，激發(fā)炮點分布不規(guī)律，各個單炮激發(fā)能量差異大，這些因素使得原本反射信號較弱的深地震反射數(shù)據(jù)更不易進行精細數(shù)據(jù)處理，利用常規(guī)的處理流程難以對深達莫霍面的包括淺部、中部以及深部地殼區(qū)域高信噪比高精度成像。

2、目前，深地震反射探測已經(jīng)在多個造山帶區(qū)域開展，目的是探測和揭露地殼結構與構造、殼-幔交互關系以及主要深大斷裂的深淺構造關系等，對地殼結構變形演化、地震深部孕震構造以及礦產(chǎn)資源形成環(huán)境等研究意義重大。由于其探測深度大，在進行地震數(shù)據(jù)采集時常采用“大、中、小”炮等不同激發(fā)藥量模式進行采集，獲得的原始單炮品質(zhì)差異大，另由于造山帶區(qū)域地表地震地質(zhì)條件復雜，激發(fā)炮點不能按照理論設計布置，位置分散雜亂，這些因素導致地震數(shù)據(jù)空間采樣不均勻，多道疊加后的數(shù)據(jù)為不規(guī)則數(shù)據(jù)，在進行偏移處理時容易產(chǎn)生空間假頻，使地震反射波成像質(zhì)量下降，造成諸如成像腳印、成像噪音以及偏移畫弧等問題。另外，不規(guī)則的數(shù)據(jù)使得高精度偏移處理方法的優(yōu)勢發(fā)揮不出來，導致反射波成像振幅不保真，成像精度不高，處理獲得的地震時間剖面不能真實反映地下地質(zhì)構造形態(tài)和物性變化規(guī)律。

3、因此，在對深地震反射探測資料進行數(shù)據(jù)處理時，除了精細的靜校正和去噪處理，以及針對地殼淺、中、深部速度差異大而采取的速度模型分析外，還需要進行數(shù)據(jù)規(guī)則化和插值重建處理。該步處理方法將空間上不規(guī)則的樣點插值成規(guī)則均勻樣點，在一定范圍內(nèi)補全缺失的不規(guī)則地震道，使內(nèi)插后的地震數(shù)據(jù)具有完整的空間和偏移距信息，消除地震數(shù)據(jù)空間采樣不規(guī)則對成像效果和成像精度造成的影響，為高精度的疊前偏移處理提供了優(yōu)良的前置數(shù)據(jù)。

技術實現(xiàn)思路

1、以下給出一個或多個方面的簡要概述以提供對這些方面的基本理解。此概述不是所有構想到的方面的詳盡綜覽，并且既非旨在指認出所有方面的關鍵性或決定性要素亦非試圖界定任何或所有方面的范圍。其唯一的目的是要以簡化形式給出一個或多個方面的一些概念以為稍后給出的更加詳細的描述之序。

2、本發(fā)明主要的目的是解決現(xiàn)有技術中所存在的技術問題，提供了一種種地震反射探測時間剖面成像方法及系統(tǒng)。該系統(tǒng)及方法采用回轉波層析靜校正、疊前去噪、地表一致性振幅補償、速度分析和模擬退火-蒙特卡洛地表一致性剩余靜校正以及基于omp算法的多維數(shù)據(jù)規(guī)則化插值的組合處理流程，對深地震反射探測數(shù)據(jù)進行疊前處理，最后采用克?；舴虔B前時間和疊前深度偏移處理，達到深地震反射探測數(shù)據(jù)處理保真保幅精準成像的目的。

3、為解決上述問題，本發(fā)明的方案是：

4、一種地震反射探測時間剖面成像方法，包括：

5、步驟1，原始地震數(shù)據(jù)解編，建立炮點和檢測點的空間屬性；

6、步驟2，計算初至波走時得到近地表速度結構，求取使初至光滑連續(xù)、疊加剖面上有效反射波分辨率和信噪比高的靜校正量，消除地震反射波到時誤差；

7、步驟3，對深地震反射數(shù)據(jù)中不同種類、不同頻率以及不同深度區(qū)域的噪音進行壓制和去除，提高有效反射信號的信噪比和分辨率；

8、步驟4，依據(jù)深地震反射探測區(qū)由淺到深的速度變化規(guī)律，建立各標志層反射波能量團突出、層次清楚的速度譜；再運用模擬退火算法和蒙特卡洛剩余靜校正與精細速度分析迭代方式，得到標志層反射波能量團位置固定的速度函數(shù)和靜校正量；

9、步驟5，采用基于omp算法的多維數(shù)據(jù)規(guī)則化插值處理技術對疊加時間剖面的信噪比和分辨率進行處理；

10、步驟6，采用克?；舴虔B前時間和疊前深度偏移處理，獲得深地震反射探測偏移地震時間剖面。

11、優(yōu)選的，上述的一種地震反射探測時間剖面成像方法，步驟2中，采用回轉波層析靜校正處理方式，通過射線追蹤算法，計算初至波走時得到近地表速度結構，求取合適的靜校正量。

12、優(yōu)選的，上述的一種地震反射探測時間剖面成像方法，所述步驟3中，采用“分類、分步、分域、分區(qū)、分頻、分時”疊前去噪處理方式，針對深地震反射數(shù)據(jù)中不同種類、不同頻率以及不同深度區(qū)域的噪音進行壓制和去除，提高有效反射信號的信噪比和分辨率。

13、優(yōu)選的，上述的一種地震反射探測時間剖面成像方法，所述步驟4中，使剩余靜校正最終收斂到一個樣點范圍內(nèi)，提高成像精度。

14、一種地震反射探測時間剖面成像系統(tǒng)，包括：

15、觀測系統(tǒng)定義模塊，原始地震數(shù)據(jù)解編，建立炮點和檢測點的空間屬性；

16、靜校正量獲取模塊，計算初至波走時得到近地表速度結構，求取使初至光滑連續(xù)、疊加剖面上有效反射波分辨率和信噪比高的的靜校正量，消除地震反射波到時誤差；

17、疊加前處理模塊，對深地震反射數(shù)據(jù)中不同種類、不同頻率以及不同深度區(qū)域的噪音進行壓制和去除，提高有效反射信號的信噪比和分辨率；

18、迭代處理模塊，依據(jù)深地震反射探測區(qū)由淺到深的速度變化規(guī)律，建立各標志層反射波能量團突出、層次清楚的速度譜；再運用模擬退火算法和蒙特卡洛剩余靜校正與精細速度分析迭代方式，得到標志層反射波能量團位置固定的速度函數(shù)和靜校正量；

19、插值處理模塊，采用基于omp算法的多維數(shù)據(jù)規(guī)則化插值處理技術對疊加時間剖面的信噪比和分辨率進行處理；

20、時間剖面處理模塊，采用克?；舴虔B前時間和疊前深度偏移處理，獲得深地震反射探測偏移地震時間剖面。

21、優(yōu)選的，上述的一種地震反射探測時間剖面成像系統(tǒng)，靜校正量獲取模塊中，采用回轉波層析靜校正處理方式，通過射線追蹤算法，計算初至波走時得到近地表速度結構，求取合適的靜校正量。

22、優(yōu)選的，上述的一種地震反射探測時間剖面成像系統(tǒng)，疊加前處理模塊中，采用“分類、分步、分域、分區(qū)、分頻、分時”疊前去噪處理方式，針對深地震反射數(shù)據(jù)中不同種類、不同頻率以及不同深度區(qū)域的噪音進行壓制和去除，提高有效反射信號的信噪比和分辨率。

23、優(yōu)選的，上述的一種地震反射探測時間剖面成像系統(tǒng)，所述迭代處理模塊中，使剩余靜校正最終收斂到一個樣點范圍內(nèi)，提高成像精度。

24、與現(xiàn)有深地震反射探測數(shù)據(jù)處理技術相比，本發(fā)明帶來了以下有益技術效果：

25、本發(fā)明所用的深地震反射數(shù)據(jù)處理流程基于常規(guī)的地震勘探數(shù)據(jù)處理流程，不需要另做處理流程測試。

26、本發(fā)明針對深地震反射探測提出了數(shù)據(jù)處理中的關鍵處理方法組合步驟及流程，采用回轉波層析靜校正、疊前去噪、地表一致性振幅補償、速度分析和模擬退火-蒙特卡洛地表一致性剩余靜校正以及基于omp算法的多維數(shù)據(jù)規(guī)則化插值的組合處理流程，最后進行克希霍夫疊前時間和疊前深度偏移處理，消除了地形復雜區(qū)地震數(shù)據(jù)空間采樣不規(guī)則對成像效果和成像精度造成的影響，提高了疊加時間剖面的信噪比，獲得了高精度的疊前偏移地震時間剖面。