[發(fā)明專利]一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法

有效

• 申請?zhí)枺?/td>: 201310553398.7
• 申請日: 2013-11-08
• 公開（公告）號(hào): CN103645498A
• 公開（公告）日: 2014-03-19
• 發(fā)明（設(shè)計(jì)）人: 陳海峰;李向陽;狄?guī)妥?/a>;魏建新;蔡志光;陳雙全
• 申請（專利權(quán)）人: 中國石油大學(xué)（北京）
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 北京匯澤知識(shí)產(chǎn)權(quán)代理有限公司 11228
• 代理人: 張瑾

1.一種直接生成共轉(zhuǎn)換點(diǎn)道集的方法，具體步驟包括：

1）采集原始轉(zhuǎn)換波地震數(shù)據(jù)，經(jīng)預(yù)處理和速度分析，得到轉(zhuǎn)換波動(dòng)校正和共轉(zhuǎn)換點(diǎn)位置計(jì)算需要的參數(shù)，并對預(yù)處理后轉(zhuǎn)換波地震數(shù)據(jù)進(jìn)行動(dòng)校正，且通過切除消除動(dòng)校拉伸的影響，

步驟1）所述的預(yù)處理包括加載觀測系統(tǒng)、靜校正、波場分離、疊前去噪、振幅恢復(fù)、預(yù)置共轉(zhuǎn)換點(diǎn)CCP號(hào)；

2）步驟1）處理后的轉(zhuǎn)換波地震數(shù)據(jù)按炮檢距大小進(jìn)行分組，每一組對應(yīng)一個(gè)炮檢距范圍；

3）輸入一組炮檢距范圍內(nèi)的數(shù)據(jù)；

4）按道頭中的共轉(zhuǎn)換點(diǎn)CCP號(hào)讀取轉(zhuǎn)換波用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù)；

5）對地震道進(jìn)行分時(shí)窗處理，

步驟5）所述的分時(shí)窗處理是指將地震道分成多個(gè)時(shí)窗，時(shí)窗長度不變，時(shí)窗重疊，重疊時(shí)窗的重疊長度為時(shí)窗內(nèi)樣點(diǎn)數(shù)的百分比perc，perc取值為0-100；

6）從第一個(gè)時(shí)窗開始，計(jì)算時(shí)窗中點(diǎn)的共轉(zhuǎn)換點(diǎn)CCP號(hào)，并由這個(gè)共轉(zhuǎn)換點(diǎn)號(hào)重新獲得該時(shí)窗用來計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù)，之后重新計(jì)算時(shí)窗中點(diǎn)相對準(zhǔn)確的共轉(zhuǎn)換點(diǎn)CCP號(hào)，直接把時(shí)窗內(nèi)的樣點(diǎn)值與窗函數(shù)相乘疊在共轉(zhuǎn)換點(diǎn)地震道上，即轉(zhuǎn)換波共反射面元的寬度上，相應(yīng)共轉(zhuǎn)換點(diǎn)地震道每個(gè)樣點(diǎn)的覆蓋次數(shù)加上窗函數(shù)的數(shù)值，

步驟6）所述的計(jì)算轉(zhuǎn)換點(diǎn)位置的參數(shù)是預(yù)先速度分析得到的，是轉(zhuǎn)換波均方根速度縱波均方根速度各向異性參數(shù)χ_eff、縱橫波垂直速度比γ₀和縱橫波有效速度比r_eff，

步驟6）所述的共轉(zhuǎn)換點(diǎn)CCP號(hào)，是指按工區(qū)面元?jiǎng)澐值模c常規(guī)處理中縱波共中心點(diǎn)CMP號(hào)一致，

步驟6）所述的共轉(zhuǎn)換點(diǎn)CCP號(hào)是由下述方法計(jì)算的：

x_C＝x_M+D

式中x_M、x_C為共中心點(diǎn)CMP和共轉(zhuǎn)換點(diǎn)CCP點(diǎn)的水平坐標(biāo)，D為轉(zhuǎn)換點(diǎn)偏離中心點(diǎn)的水平距離，

D=xp-x2]]>

式中x_p為共轉(zhuǎn)換點(diǎn)與炮點(diǎn)之間的水平距離，

已知共中心點(diǎn)CMP號(hào)計(jì)算相應(yīng)的共轉(zhuǎn)換點(diǎn)CCP號(hào)

jxC=INTn(xCΔxm)=INTn(xM+DΔxm)=jxM+INTn(DΔxm)]]>

其中：INT_n表示取最接近的整數(shù)，取整隱含做面元化處理,Δx_m為工區(qū)面元大小，

所述的共轉(zhuǎn)換點(diǎn)與炮點(diǎn)之間的水平距離，采用單層各向同性介質(zhì)迭代解、解析解、多層介質(zhì)指定旅行時(shí)各向同性近似解、多層VTI介質(zhì)指定旅行時(shí)近似解、傾斜界面迭代解、解析解計(jì)算得到，

步驟6）所述的直接把時(shí)窗內(nèi)的樣點(diǎn)值與窗函數(shù)相乘疊在共轉(zhuǎn)換點(diǎn)地震道上，即轉(zhuǎn)換波共反射面元的寬度上，是把這個(gè)時(shí)窗內(nèi)的樣點(diǎn)疊加在準(zhǔn)確的共轉(zhuǎn)換點(diǎn)CCP號(hào)上以及相鄰的I個(gè)共轉(zhuǎn)換點(diǎn)CCP號(hào)上，I=0,1,2,…,INT(2·m/Δx+0.5)，m為轉(zhuǎn)換波共反射面元的寬度，INT表示取整，

步驟6）所述的轉(zhuǎn)換波共反射面元的寬度采用以下關(guān)系式確定：

m=1Δm=0[1/(1-Δm)+0.5]Δm≥0.5[1/Δm+0.5]Δm<0.5]]>

Δm=(1+r)·Δsr·Δx-INT[(1+r)·Δsr·Δx]]]>

其中：Δx為道間距，Δs為炮間距，r為縱橫波速度比，INT表示取整；

7）重復(fù)步驟6）直到該道的每個(gè)時(shí)窗完成；

8）重復(fù)步驟4)至7)直到炮檢距分組內(nèi)所有道輸入完畢，輸出該組的共轉(zhuǎn)換點(diǎn)地震道，

步驟8）所述的輸出該組的共轉(zhuǎn)換點(diǎn)地震道為共轉(zhuǎn)換點(diǎn)地震道以及共轉(zhuǎn)換點(diǎn)地震道的覆蓋次數(shù)，并在道頭中置上該組炮檢距的大小；

9）重復(fù)步驟3)至8)直到所有的炮檢距分組輸入完畢，就得到了時(shí)間和空間變化的共轉(zhuǎn)換點(diǎn)疊前道集。