一種直接生成共轉換點道集的方法
【專利摘要】本發明提供一種直接生成共轉換點道集的方法，假設轉換波地震數據已預處理好，已獲得準確的動校正參數和計算轉換點位置的速度等參數，按炮檢距分組，分時窗疊加獲得時間和空間變化的共轉換點道集。本發明實現簡便、適應大數據量處理。
【專利說明】一種直接生成共轉換點道集的方法
【技術領域】
[0001]本發明涉及轉換波地震資料處理【技術領域】，具體地，涉及一種直接時空變生成共轉換點道集的方法。
【背景技術】
[0002]在多分量轉換波地震資料處理中，共轉換點道集的抽取方法一直是一項關鍵技術。因為轉換波的下行波是縱波、上行波是橫波，造成傳播路徑的不對稱性，轉換點位于從中心點偏向接收點的一側，且隨深度及界面上下縱橫波速度的變化而變化。這使得基于常規共中心點(CMP)道集的速度分析、動校正和水平疊加技術不能適用于轉換波，轉換波疊加成像必須按共轉換點進行道集分選。目前，有關轉換點位置的計算已有很多方法，主要包括漸近線法(ACP)、單層各向同性介質迭代解、多層介質指定深度各向同性近似解(Thomsen，1999)、多層介質指定旅行時各向同性近似解(Thomsen, 1999)和多層VTI介質指定旅行時近似解(Yuan&Li，2001)等。另外，周竹生等(1993)導出了一種快速、高精度轉換點的迭代算法；許士勇等(2002)給出了四次方程的近似表達式；姚陳等(2005)給出了三維傾角轉換波共轉換點的解析表達式；苑春方等(2005，2006，2006)推導出了單層各向同性介質與深度有關的唯一的解析解和含有旅行時間唯一解析解，傾斜界面上的轉換點的精確解；Chang等(2005)利用傾斜界面上轉換波射線路徑的三角幾何關系，推導出傾斜界面上的轉換點的近似公式；畢麗飛等(2005)用水平界面推導轉換點坐標迭代解類似的方法。但目前，共轉換點道集的選排方法仍然是以整道選排為主，主要是漸近逼近道分類法、依賴深度的選排方法。漸近逼近道分類法是通過一個工區的固定速度比，計算轉換點的漸近線，然后將整個工區的轉換波數據抽取到漸近線道集中，即ACP道集，這種方法能一定程度上能夠滿足深部地層的成像要求，但是淺層的成像效果較差。依賴深度的選排方法只對某一深度最為正確，都是能在目標層聚焦成像，更適合目的層較淺的轉換波處理，更適合更新速度分析，并且指定目的層計算提供了多種轉換點位置的計算方法，但在所選反射層的上面或下面，反射界面的疊后成像可能會變得模糊。另外，分層選排方法能夠使不同深度的轉換波均能達到較好的聚焦效果，但需要在不同層之間進行道集的拼接，特別是在大數據量的情況下實際操作困難。
[0003]上述的共轉換點道集的選排方法都有一定的應用條件，并不能很好的滿足轉換波對轉換點在時間上和空間上變化的要求。
【發明內容】
[0004]本發明目的是提供一種能適應時間上和空間上變化的共轉換點道集的選排方法，一種直接時空變生成共轉換點道集的方法。
[0005]本發明通過以下技術方案實現:
[0006]一種直接生成共轉換點道集的方法，具體步驟包括:
[0007]I)采集原始轉換波地震數據，經預處理和速度分析，得到轉換波動校正和共轉換點位置計算需要的參數，并對預處理后轉換波地震數據進行動校正，且通過切除消除動校拉伸的影響，
[0008]步驟1)所述的預處理包括加載觀測系統、靜校正、波場分離、疊前去噪、振幅恢復、預置共轉換點CCP號；
[0009]2)步驟1)處理后的轉換波地震數據按炮檢距大小進行分組，每一組對應一個炮檢距范圍；
[0010]3)輸入一組炮檢距范圍內的數據；
[0011]4)按道頭中的共轉換點CCP號讀取轉換波用來計算轉換點位置的參數；
[0012]5)對地震道進行分時窗處理，
[0013]步驟5)所述的分時窗處理是指將地震道分成多個時窗，時窗長度不變，時窗重疊，重疊時窗的重疊長度為時窗內樣點數的百分比perc, perc取值為0-100 ；
[0014]6)從第一個時窗開始，計算時窗中點的共轉換點CCP號，并由這個共轉換點號重新獲得該時窗用來計算轉換點位置的參數，之后重新計算時窗中點相對準確的共轉換點CCP號，直接把時窗內的樣點值與窗函數相乘疊在共轉換點地震道上，即轉換波共反射面元的寬度上，相應共轉換點地震道每個樣點的覆蓋次數加上窗函數的數值，
[0015]步驟6)所述的計算轉換點位置的參數是預先速度分析得到的，是轉換波均方根速度ve2、縱波均方根速度各向異性參數Xrff、縱橫波垂直速度比Y0和縱橫波有效速度
比 reff，
[0016]步驟6)所述的共轉換點CCP號，是指按工區面元劃分的，與常規處理中縱波共中心點CMP號一致，
[0017]步驟6)所述的共轉換點CCP號是由下述方法計算的:
[0018]xc = xm+D (1)
[0019]式中xM、xc為共中心點CMP和共轉換點CCP點的水平坐標，D為轉換點偏離中心點的水平距離，
[0020]D = xp -x/2(2)
[0021]式中Xp為共轉換點與炮點之間的水平距離，
[0022]已知共中心點CMP號As，計算相應的共轉換點CCP號jxc:
[0023]
【權利要求】
1.一種直接生成共轉換點道集的方法，具體步驟包括: 1)采集原始轉換波地震數據，經預處理和速度分析，得到轉換波動校正和共轉換點位置計算需要的參數，并對預處理后轉換波地震數據進行動校正，且通過切除消除動校拉伸的影響， 步驟I)所述的預處理包括加載觀測系統、靜校正、波場分離、疊前去噪、振幅恢復、預置共轉換點CCP號； 2)步驟I)處理后的轉換波地震數據按炮檢距大小進行分組，每一組對應一個炮檢距范圍； 3)輸入一組炮檢距范圍內的數據； 4)按道頭中的共轉換點CCP號讀取轉換波用來計算轉換點位置的參數； 5)對地震道進行分時窗處理， 步驟5)所述的分時窗處理是指將地震道分成多個時窗，時窗長度不變，時窗重疊，重疊時窗的重疊長度為時窗內樣點數的百分比perc, perc取值為0-100 ； 6)從第一個時窗開始，計算時窗中點的共轉換點CCP號，并由這個共轉換點號重新獲得該時窗用來計算轉換點位置的參數，之后重新計算時窗中點相對準確的共轉換點CCP號，直接把時窗內的樣點值與窗函數相乘疊在共轉換點地震道上，即轉換波共反射面元的寬度上，相應共轉換點地震道每個樣點的覆蓋次數加上窗函數的數值， 步驟6)所述的計算轉換點位置的參數是預先速度分析得到的，是轉換波均方根速度D縱波均方根速度％、各向異性參數Xeff、縱橫波垂直速度比Ytl和縱橫波有效速度比reff， 步驟6)所述的共轉換點CCP號，是指按工區面元劃分的，與常規處理中縱波共中心點CMP號一致， 步驟6)所述的共轉換點CCP號是由下述方法計算的:
xC = xm+D 式中xM、x。為共中心點CMP和共轉換點CCP點的水平坐標，D為轉換點偏離中心點的水平距離， 式中Xp為共轉換點與炮點之間的水平距離， 已知共中心點CMP號4*，計算相應的共轉換點CCP號Lc'
【文檔編號】G01V1/28GK103645498SQ201310553398
【公開日】2014年3月19日 申請日期:2013年11月8日 優先權日:2013年11月8日
【發明者】陳海峰, 李向陽, 狄幫讓, 魏建新, 蔡志光, 陳雙全 申請人:中國石油大學(北京)
專業數控銑床產品供應商
20年生產經驗，實力工廠，精工品質，質量有保障
