一種獲得初至波射線的方法和裝置的制造方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地球物理勘探技術(shù)領(lǐng)域，尤其涉及的一種獲得初至波射線的方法和裝 置。
【背景技術(shù)】
[0002] 地震勘探中的靜校正處理技術(shù)可以用來(lái)解決地震波在近地表介質(zhì)旅行時(shí)間變化 影響反射波疊加成像的問題。在靜校正處理技術(shù)中，初至是地震資料的一個(gè)重要屬性。在 地震資料處理過(guò)程中，通?？梢岳贸踔恋膶傩越鉀Q地震數(shù)據(jù)的靜校正問題。當(dāng)?shù)卣鸢l(fā)生 后，由于各種地震波的傳播速度不同，傳播到觀測(cè)點(diǎn)的時(shí)間也就有先后。地震觀測(cè)點(diǎn)最先接 收到的波就可以稱為初至波。
[0003] 野外近地表模型大多都是層狀的，而目前的商業(yè)軟件中使用的模型都是規(guī)則網(wǎng)格 模型來(lái)等效層狀模型，并且從淺到深，速度呈連續(xù)狀，這種射線追蹤的方法只適應(yīng)規(guī)則網(wǎng)格 的速度連續(xù)變化的模型，而無(wú)法適應(yīng)速度突變的層狀模型，這與實(shí)際地表情況不符。因此， 不能準(zhǔn)確的獲得初至波射線，也無(wú)法反演出相對(duì)精確的近地表模型。

【發(fā)明內(nèi)容】

[0004] 本申請(qǐng)的目的是提供一種獲得初至波射線的方法和裝置，以獲得準(zhǔn)確的層狀模型 的初至波射線。
[0005] 為了實(shí)現(xiàn)上述目的，本申請(qǐng)?zhí)峁┝?一種獲得初至波射線的方法，該方法包括：
[0006] 采集地震數(shù)據(jù)；
[0007] 根據(jù)所述地震數(shù)據(jù)建立層狀的速度模型；
[0008] 獲取所述速度模型中的層厚度和層速度，根據(jù)所述層厚度和層速度獲取激發(fā)點(diǎn)和 接收點(diǎn)到目的層的延遲時(shí)；
[0009] 根據(jù)所述延遲時(shí)、所述層速度和每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離獲取激發(fā)點(diǎn)通過(guò)每個(gè) 層的界面到達(dá)接收點(diǎn)的旅行時(shí)間；
[0010] 比較所述旅行時(shí)間，取所述旅行時(shí)間的最小值的射線，所述射線是初至波射線。
[0011] 在一個(gè)優(yōu)選的實(shí)施例中，該方法還包括：對(duì)所述初至波射線進(jìn)行網(wǎng)格剖分，獲取射 線穿過(guò)網(wǎng)格的索引和所述網(wǎng)格中射線旅行的距離。
[0012] 在一個(gè)優(yōu)選的實(shí)施例中，所述根據(jù)所述地震數(shù)據(jù)建立層狀的速度模型包括：
[0013] 根據(jù)所述地震數(shù)據(jù)以接收點(diǎn)和激發(fā)點(diǎn)的高程建立地表面；
[0014] 從所述地震數(shù)據(jù)獲取最小高程值，根據(jù)所述最小高程值獲得模型的水平底界面；
[0015] 從所述地震數(shù)據(jù)獲取地震道道距，根據(jù)所述道距設(shè)置模型水平方向的網(wǎng)格間距；
[0016] 根據(jù)地震數(shù)據(jù)獲取垂直方向的層數(shù)；
[0017] 根據(jù)地震數(shù)據(jù)線性內(nèi)插每層面上網(wǎng)格節(jié)點(diǎn)的高程值，層間用常速填充。
[0018] 在一個(gè)優(yōu)選的實(shí)施例中，所述根據(jù)所述層厚度和層速度獲取激發(fā)點(diǎn)和接收點(diǎn)到目 的層的延遲時(shí)包括：
[0019] 根據(jù)所述層厚度和層速度并利用公式獲取激發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)，公 式如下：
[0020]
【主權(quán)項(xiàng)】
1. 一種獲得初至波射線的方法，其特征在于，該方法包括： 采集地震數(shù)據(jù)； 根據(jù)所述地震數(shù)據(jù)建立層狀的速度模型； 獲取所述速度模型中的層厚度和層速度，根據(jù)所述層厚度和層速度獲取激發(fā)點(diǎn)和接收 點(diǎn)到目的層的延遲時(shí)； 根據(jù)所述延遲時(shí)、所述層速度和每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的 界面到達(dá)接收點(diǎn)的旅行時(shí)間； 比較所述旅行時(shí)間，取所述旅行時(shí)間的最小值的射線，所述射線是初至波射線。
2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，該方法還包括：對(duì)所述初至波射線進(jìn)行網(wǎng) 格剖分，獲取射線穿過(guò)網(wǎng)格的索引和所述網(wǎng)格中射線旅行的距離。
3. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述根據(jù)所述地震數(shù)據(jù)建立層狀的速 度模型包括： 根據(jù)所述地震數(shù)據(jù)以接收點(diǎn)和激發(fā)點(diǎn)的高程建立地表面； 從所述地震數(shù)據(jù)獲取最小高程值，根據(jù)所述最小高程值獲得模型的水平底界面； 從所述地震數(shù)據(jù)獲取地震道道距，根據(jù)所述道距設(shè)置模型水平方向的網(wǎng)格間距； 根據(jù)地震數(shù)據(jù)獲取垂直方向的層數(shù)； 根據(jù)地震數(shù)據(jù)線性內(nèi)插每層面上網(wǎng)格節(jié)點(diǎn)的高程值，層間用常速填充。
4. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述根據(jù)所述層厚度和層速度獲取激 發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)包括： 根據(jù)所述層厚度和層速度并利用公式獲取激發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)，公式如 下：
其中，
，Dn代表激發(fā)點(diǎn)或接收點(diǎn)到第n層目的層的延遲時(shí)（l〈n〈 = 6)，H玳 表第i層在激發(fā)點(diǎn)或接收點(diǎn)位置處的厚度，Vn代表第n層目的層的速度。
5. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，所述根據(jù)所述延遲時(shí)、所述層速度和 每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的界面到達(dá)接收點(diǎn)的旅行時(shí)間包括： 根據(jù)觀測(cè)系統(tǒng)，獲取每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離，利用公式獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的 界面到達(dá)接收點(diǎn)的旅行時(shí)間，公式如下：
其中，1；為激發(fā)點(diǎn)通過(guò)第n層的界面到達(dá)接收點(diǎn)的旅行時(shí)間，Dsn代表激發(fā)點(diǎn)到第n層 目的層的延遲時(shí)，DmR表接收點(diǎn)到第n層目的層的延遲時(shí)，Vn代表第n層目的層的速度，X 代表每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離。
6. -種獲得初至波射線的裝置，其特征在于，該裝置包括： 地震數(shù)據(jù)采集單元，用于采集地震數(shù)據(jù)； 速度模型建立單元，用于根據(jù)所述地震數(shù)據(jù)建立層狀的速度模型； 延遲時(shí)獲取單元，用于獲取所述速度模型中的層厚度和層速度，根據(jù)所述層厚度和層 速度獲取激發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)； 旅行時(shí)間獲取單元，用于根據(jù)所述延遲時(shí)、所述層速度和每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離 獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的界面到達(dá)接收點(diǎn)的旅行時(shí)間； 初至波射線獲取單元，用于比較所述旅行時(shí)間，取所述旅行時(shí)間的最小值的射線，所述 射線是初至波射線。
7. 根據(jù)權(quán)利要求6所述的裝置，其特征在于，該裝置還包括：網(wǎng)格剖分單元，用于對(duì)所 述初至波射線進(jìn)行網(wǎng)格剖分，獲取射線穿過(guò)網(wǎng)格的索引和所述網(wǎng)格中射線旅行的距離。
8. 根據(jù)權(quán)利要求6或7所述的裝置，其特征在于，所述速度模型建立單元包括： 地表面建立模塊，用于根據(jù)所述地震數(shù)據(jù)以接收點(diǎn)和激發(fā)點(diǎn)的高程建立地表面； 水平底界面獲取模塊，用于從所述地震數(shù)據(jù)獲取最小高程值，根據(jù)所述最小高程值獲 得模型的水平底界面； 網(wǎng)格間距獲取模塊，用于從所述地震數(shù)據(jù)獲取地震道道距，根據(jù)所述道距設(shè)置模型水 平方向的網(wǎng)格間距； 層數(shù)獲取模塊，用于根據(jù)地震數(shù)據(jù)獲取垂直方向的層數(shù)； 高程值內(nèi)插模塊，用于線性內(nèi)插每層面上網(wǎng)格節(jié)點(diǎn)的高程值，層間用常速填充。
9. 根據(jù)權(quán)利要求6或7所述的裝置，其特征在于，所述根據(jù)所述層厚度和層速度獲取激 發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)包括： 根據(jù)所述層厚度和層速度并利用公式獲取激發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)，公式如 下：
其中，
，Dn代表激發(fā)點(diǎn)或接收點(diǎn)到第n層目的層的延遲時(shí)（l〈n〈 = 6)，H玳 表第i層在激發(fā)點(diǎn)或接收點(diǎn)位置處的厚度，Vn代表第n層目的層的速度。
10. 根據(jù)權(quán)利要求6或7所述的裝置，其特征在于，所述根據(jù)所述延遲時(shí)、所述層速度和 每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的界面到達(dá)接收點(diǎn)的旅行時(shí)間包括： 根據(jù)觀測(cè)系統(tǒng)，獲取每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離，利用公式獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的 界面到達(dá)接收點(diǎn)的旅行時(shí)間，公式如下：
其中，1；為激發(fā)點(diǎn)通過(guò)第n層的界面到達(dá)接收點(diǎn)的旅行時(shí)間，Dsn代表激發(fā)點(diǎn)到第n層 目的層的延遲時(shí)，DmR表接收點(diǎn)到第n層目的層的延遲時(shí)，Vn代表第n層目的層的速度，X 代表每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離。
【專利摘要】本申請(qǐng)實(shí)施例提供的一種獲得初至波射線的方法和裝置。該方法包括：采集地震數(shù)據(jù)；根據(jù)所述地震數(shù)據(jù)建立層狀的速度模型；獲取所述速度模型中的層厚度和層速度，根據(jù)所述層厚度和層速度獲取激發(fā)點(diǎn)和接收點(diǎn)到目的層的延遲時(shí)；根據(jù)所述延遲時(shí)、所述層速度和每個(gè)激發(fā)點(diǎn)到接收點(diǎn)的距離獲取激發(fā)點(diǎn)通過(guò)每個(gè)層的界面到達(dá)接收點(diǎn)的旅行時(shí)間；比較所述旅行時(shí)間，取所述旅行時(shí)間的最小值的射線，所述射線是初至波射線。利用本申請(qǐng)實(shí)施例提供的技術(shù)方案可以準(zhǔn)確的獲得初至波射線。
【IPC分類】G01V1-36, G01V1-28
【公開號(hào)】CN104570091
【申請(qǐng)?zhí)枴緾N201410777016
【發(fā)明人】閆智慧, 呂哲健, 白志宏
【申請(qǐng)人】中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司
【公開日】2015年4月29日
【申請(qǐng)日】2014年12月15日