一種基于地形因子的炮點(diǎn)偏移方法
【專(zhuān)利摘要】本發(fā)明提供了一種基于地形因子的炮點(diǎn)偏移方法，包括：根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件；根據(jù)工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值；根據(jù)每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)；在需要偏移的炮點(diǎn)的縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移。本發(fā)明實(shí)施例的基于地形因子的炮點(diǎn)偏移方法，可以實(shí)現(xiàn)在地形起伏較大的山地地震采集施工時(shí)，通過(guò)地形因子坡度定量、準(zhǔn)確地計(jì)算出炮點(diǎn)最佳的偏移位置，通過(guò)炮點(diǎn)偏移，確保覆蓋次數(shù)分布滿足采集設(shè)計(jì)要求，得到淺中層資料，從而為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
【專(zhuān)利說(shuō)明】一種基于地形因子的炮點(diǎn)偏移方法【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地球物理勘探技術(shù)，尤其涉及一種基于地形因子的炮點(diǎn)偏移方法。
【背景技術(shù)】
[0002]在復(fù)雜山地區(qū)域進(jìn)行勘探施工時(shí)，由于地表地形復(fù)雜，起伏高差大，在一些山地高陡地形區(qū)無(wú)法布設(shè)炮點(diǎn)，導(dǎo)致覆蓋次數(shù)降低，淺層資料缺失。為提高覆蓋次數(shù)確保淺層資料，需要對(duì)高陡地區(qū)無(wú)法正常布設(shè)的炮點(diǎn)在緩坡區(qū)進(jìn)行偏移和加密。
[0003]目前，對(duì)于山地高陡地區(qū)的炮點(diǎn)偏移和加密都是通過(guò)施工人員的實(shí)地觀察或借助衛(wèi)片憑經(jīng)驗(yàn)來(lái)確定，往往誤差大、耗時(shí)耗力。因此，有必要根據(jù)地形因子坡度、坡向自動(dòng)確定最佳的炮點(diǎn)位置，通過(guò)地形因子坡度定量自動(dòng)對(duì)山地高陡地形區(qū)無(wú)法布設(shè)的炮點(diǎn)進(jìn)行合理偏移，以確保山地復(fù)雜區(qū)覆蓋次數(shù)分布滿足采集設(shè)計(jì)要求，確保得到好的淺中層資料。但是當(dāng)前，尚無(wú)根據(jù)地形因子自動(dòng)確定最佳炮點(diǎn)位置進(jìn)行偏移的方法。

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

[0004]本發(fā)明的目的是實(shí)現(xiàn)一種基于地形因子的炮點(diǎn)偏移方法，可以在地表復(fù)雜、高差起伏較大的地區(qū)無(wú)法正常布設(shè)炮點(diǎn)時(shí)，利用計(jì)算的地形因子坡度自動(dòng)定量確定最佳炮點(diǎn)位置進(jìn)行偏移的方法，為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
[0005]為了達(dá)到上述目的，本發(fā)明實(shí) 施例提供一種基于地形因子的炮點(diǎn)偏移方法，包括:根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件；根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值；根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)；在所述需要偏移的炮點(diǎn)的縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移。
[0006]進(jìn)一步地，在一實(shí)施例中，所述根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件，包括:利用遙感技術(shù)獲取所述工區(qū)數(shù)字高程模型；根據(jù)所述工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，所述地形因子坡度值S是地形曲面函數(shù)
Z=f(x, y)在東西、南北方向上高程變化率的函數(shù)，數(shù)學(xué)表示為:S =肌.吆」八:;其中:
fx為所述工區(qū)數(shù)字高程模型中X方向南北上的高程變化率；fy為所述工區(qū)數(shù)字高程模型中Y方向東西上的高程變化率；將計(jì)算得到的所述工區(qū)地形因子坡度值與所述工區(qū)數(shù)字高程模型相結(jié)合，得到所述工區(qū)坡度數(shù)字模型文件。
[0007]進(jìn)一步地，在一實(shí)施例中，所述根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值，包括:所述根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用4點(diǎn)雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值。
[0008]進(jìn)一步地，在一實(shí)施例中，所述根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)，包括:如果某炮點(diǎn)的坡度值超出設(shè)定的坡度限值，則該炮點(diǎn)需要偏移到平緩地方。[0009]進(jìn)一步地，在一實(shí)施例中，所述在需要偏移的炮點(diǎn)縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移，包括:將所述需要偏移的炮點(diǎn)沿接收線方向進(jìn)行縱向偏移，以最淺目的層深度為炮點(diǎn)縱向偏移的最大量，縱向偏移量為道距的整數(shù)倍，計(jì)算每一縱向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
[0010]進(jìn)一步地，在一實(shí)施例中，所述在需要偏移的炮點(diǎn)縱向和橫向偏移限制值范圍內(nèi)，選取滿足地形因子坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移，包括:將所述需要偏移的炮點(diǎn)沿在垂直接收線方向進(jìn)行橫向偏移，以最淺目的層深度為炮點(diǎn)橫向偏移的最大量，橫向偏移量是接收線距的整數(shù)倍，計(jì)算每一橫向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
[0011]進(jìn)一步地，在一實(shí)施例中，當(dāng)所述需要偏移的炮點(diǎn)偏移后的位置與正常炮點(diǎn)重合時(shí)，所述正常炮點(diǎn)也要依次偏移。
[0012]進(jìn)一步地，在一實(shí)施例中，所述炮點(diǎn)偏移方法還包括:在對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移前，計(jì)算所述需要偏移的炮點(diǎn)偏移前的覆蓋次數(shù)，并且對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移后，計(jì)算所述需要偏移的炮點(diǎn)偏移后的覆蓋次數(shù)，并進(jìn)行對(duì)比。
[0013]本發(fā)明實(shí)施例的基于地形因子的炮點(diǎn)偏移方法，可以實(shí)現(xiàn)在地形起伏較大的山地地震采集施工時(shí)，通過(guò)地形因子坡度定量、準(zhǔn)確地計(jì)算出炮點(diǎn)最佳的偏移位置，通過(guò)炮點(diǎn)偏移，確保覆蓋次數(shù)分布滿足采集設(shè)計(jì)要求，得到淺中層資料，從而為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
【專(zhuān)利附圖】

【附圖說(shuō)明】
[0014]為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。
[0015]圖1為本發(fā)明實(shí)施例的基于地形因子的炮點(diǎn)偏移方法的流程示意圖；
[0016]圖2為本發(fā)明另一實(shí)施例的基于地形因子的炮點(diǎn)偏移方法的流程示意圖；
[0017]圖3為一具體實(shí)施例中的工區(qū)數(shù)字高程模型圖；
[0018]圖4為一具體實(shí)施例中的工區(qū)地形因子坡度數(shù)字模型圖；
[0019]圖5為一具體實(shí)施例中的炮點(diǎn)偏移分布圖；
[0020]圖6為一具體實(shí)施例中的炮點(diǎn)偏移前的覆蓋次數(shù)分布圖；
[0021]圖7為一具體實(shí)施例中的炮點(diǎn)偏移后的覆蓋次數(shù)分布圖。
【具體實(shí)施方式】
[0022]下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0023]本發(fā)明實(shí)施例的主要目的是在地表復(fù)雜、高差起伏較大的地區(qū)無(wú)法正常布設(shè)炮點(diǎn)時(shí)，利用計(jì)算的地形因子坡度自動(dòng)定量確定最佳炮點(diǎn)位置進(jìn)行偏移，為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
[0024]圖1為本發(fā)明實(shí)施例的基于地形因子的炮點(diǎn)偏移方法的流程示意圖。如圖所示，本實(shí)施例的炮點(diǎn)偏移方法包括:步驟S101，根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件；步驟S102，根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值；步驟S103，根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)；步驟S104，在所述需要偏移的炮點(diǎn)的縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移。
[0025]在本實(shí)施例的步驟SlOl中，可以利用遙感技術(shù)獲取所述工區(qū)數(shù)字高程模型。并且，根據(jù)所述工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，可以利用一下公式實(shí)現(xiàn):所述地形因子坡度值S是地形曲面函數(shù)Z=f (X，y)在東西、南北方向上高程變化率的函數(shù)，數(shù)學(xué)表
示為:S = arc(g JJJTTt Ii中:fx為所述工區(qū)數(shù)字高程模型中X方向南北上的高程變化
率；fy為所述工區(qū)數(shù)字高程模型中Y方向東西上的高程變化率。將計(jì)算得到的所述工區(qū)地形因子坡度值與所述工區(qū)數(shù)字高程模型相結(jié)合，得到所述工區(qū)坡度數(shù)字模型文件，該工區(qū)坡度數(shù)字模型文件為網(wǎng)格數(shù)據(jù)文件。
[0026]在本實(shí)施例的步驟S102中，根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值，包括:所述根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用4點(diǎn)雙線性?xún)?nèi)插炮點(diǎn)位置所在的坡度網(wǎng)格，得到每一炮點(diǎn)的坡度值。
[0027]在本實(shí)施例的步驟S103中，根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)，包括:如果某炮點(diǎn)的坡度值超出設(shè)定的坡度限值，則該炮點(diǎn)需要偏移到平緩地方。例如，如果高陡山地作業(yè)炮點(diǎn)地形因子坡度值如果超出40度，施工難度與風(fēng)險(xiǎn)大，需要偏移到該炮點(diǎn)到平緩地方。
[0028]在本實(shí)施例的步驟S104中，在需要偏移的炮點(diǎn)縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移，包括如下三種情況:
[0029](1)、首先將所述需要偏移的炮點(diǎn)沿接收線方向進(jìn)行縱向偏移，以最淺目的層深度為炮點(diǎn)縱向偏移的最大量，縱向偏移量為道距的整數(shù)倍，計(jì)算每一縱向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
[0030]將需要偏移的炮點(diǎn)沿接收線方向(縱向)偏移，為了削弱炮點(diǎn)偏移對(duì)覆蓋次數(shù)的影響，以最淺目的層深度為炮點(diǎn)縱向偏移的最大量Xd，縱向偏移量為道距的整數(shù)倍，以確保覆蓋次數(shù)的均勻性。
[0031]將炮點(diǎn)沿接收線方向移動(dòng)幾個(gè)道距的距離，三個(gè)以上的需要偏移的炮點(diǎn)分布在高陡地形兩側(cè)。根據(jù)上述規(guī)則，計(jì)算可縱向偏移位置(n*Ds、(n-l)*Ds……3*Ds、2*Ds、
Ds、-Ds、-2*Ds、-3*Ds......- (n_l) *Ds、n*Ds 其中 Ds 為道距，n*Ds〈Xd)的地形因子坡度值，
在地形因子坡度限值范圍內(nèi)，選擇偏移最小的位置，連續(xù)點(diǎn)的偏移均勻分布在高陡地域兩偵||，依次進(jìn)行選擇，確定炮點(diǎn)的偏移位置。當(dāng)偏移的炮點(diǎn)與正常炮點(diǎn)重合時(shí)，正常炮點(diǎn)也需要依次偏移。[0032](2)、將所述需要偏移的炮點(diǎn)沿在垂直接收線方向進(jìn)行橫向偏移，以最淺目的層深度為炮點(diǎn)橫向偏移的最大量，橫向偏移量是接收線距的整數(shù)倍，計(jì)算每一橫向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
[0033]如果無(wú)效炮點(diǎn)無(wú)法沿接收線方向(縱向)偏移，則沿垂直接收線方向(橫向)偏移，使用地形因子坡度為標(biāo)準(zhǔn)，最大橫向偏移距離通過(guò)確定的最淺目的層的深度得到，以最淺目的層深度為炮點(diǎn)縱向偏移的最大量Yd，橫向偏移量是接收線距的整數(shù)倍，以確保覆蓋次數(shù)的均勻性。多個(gè)連續(xù)需要偏移的炮點(diǎn)按左右兩側(cè)等分，偏移距離為接收線整線距η倍。
[0034]依據(jù)上述規(guī)則，計(jì)算橫向偏移位置(n*Dl、(n_l)*Dl……2*D1、_D1、_2*D1、……-(n-l)*Dl、n*Dl，其中Dl為接收線距，n*Dl〈Yd)的坡度值，在地形因子坡度限值范圍內(nèi)，選擇偏移量最小的位置，連續(xù)點(diǎn)的偏移按左右兩側(cè)依次選擇。有效偏移范圍內(nèi)偏移整接收線距，計(jì)算炮點(diǎn)位置的地形因子坡度數(shù)值，滿足地形因子坡度限值，距離原位置的最近的為最佳炮點(diǎn)偏移位置。多個(gè)炮點(diǎn)偏移平均分布高陡地域兩側(cè)。
[0035](3)、不滿足步驟(1)、(2)要求的作為空炮點(diǎn)處理。
[0036]在另一實(shí)施例中，如圖2所示，該實(shí)施例的基于地形因子的炮點(diǎn)偏移方法包括:步驟S101，根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件；步驟S102，根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值；步驟S103，根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)；步驟S104’，在對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移前，計(jì)算所述需要偏移的炮點(diǎn)偏移前的覆蓋次數(shù)；步驟S104，在所述需要偏移的炮點(diǎn)的縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移；步驟S104’’，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移后，計(jì)算所述需要偏移的炮點(diǎn)偏移后的覆蓋次數(shù)，并與步驟S104’中的覆蓋次數(shù)進(jìn)行對(duì)比。觀測(cè)系統(tǒng)的覆蓋次數(shù)對(duì)比圖，使采集設(shè)計(jì)人員更能直觀清楚的看到炮點(diǎn)位置偏移后覆蓋次數(shù)得到了 提高。
[0037]上述實(shí)施例的基于地形因子的炮點(diǎn)偏移方法，可以實(shí)現(xiàn)在地形起伏較大的山地地震采集施工時(shí)，通過(guò)地形因子坡度定量、準(zhǔn)確地計(jì)算出炮點(diǎn)最佳的偏移位置，通過(guò)炮點(diǎn)偏移，確保覆蓋次數(shù)分布滿足采集設(shè)計(jì)要求，得到淺中層資料，從而為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
[0038]具體實(shí)施例:
[0039]本具體實(shí)施例可以應(yīng)用到山地施工的觀測(cè)系統(tǒng)設(shè)計(jì)，通過(guò)利用地形因子坡度進(jìn)行炮點(diǎn)偏移加密，彌補(bǔ)地形復(fù)雜空炮太多造成的覆蓋次數(shù)不足，同時(shí)也可以實(shí)現(xiàn)同一觀測(cè)系統(tǒng)確定炮點(diǎn)偏移前后的覆蓋次數(shù)對(duì)比，具體實(shí)現(xiàn)步驟如下:
[0040]( I)工區(qū)觀測(cè)系統(tǒng)與基礎(chǔ)數(shù)據(jù)
[0041]觀測(cè)系統(tǒng):32L3S180R，道距50米，炮點(diǎn)距100米，接收線距200米，炮線距200米；
工區(qū)南部山地起伏較大，部分地區(qū)施工困難。
[0042]要確定該觀測(cè)系統(tǒng)的炮點(diǎn)偏移位置，首先要通過(guò)工區(qū)的數(shù)字高程模型(如圖3所示)，計(jì)算工區(qū)地形因子坡度數(shù)字模型(如圖4所示)，計(jì)算各炮點(diǎn)的地形因子坡度，與地形因子坡度限值(坡度值小于40)對(duì)比，確定炮點(diǎn)偏移位置，如圖5所示的炮點(diǎn)偏移分布圖，從圖中可以看出需要偏移的炮點(diǎn)。[0043](2)沿接收線方向確定炮點(diǎn)偏移位置。如:
[0044]使用采集設(shè)計(jì)要求沿接收線方向進(jìn)行炮點(diǎn)偏移，炮點(diǎn)沿線方向偏移限值300米，偏移距離為整道距50米的整數(shù)倍。
[0045]例如單炮點(diǎn)進(jìn)行偏移，計(jì)算需要偏移的炮點(diǎn)在允許偏移范圍內(nèi)的地形因子坡度:允許偏移的范圍為300米、250米、200米、150米、100米、50米、-50米、-100米、-150米、-200米、-250米、-300米，分別計(jì)算出對(duì)應(yīng)整道距位置上地形因子坡度值為46、19、32、
18、27、57、48、52、32、36、31、26，則使用距離最近，滿足地形因子坡度限值最小的值27即為
炮點(diǎn)偏移位置。
[0046]例如連續(xù)炮點(diǎn)偏移，需要偏移3個(gè)炮點(diǎn)(a、b、C)，使用采集設(shè)計(jì)要求，炮點(diǎn)沿線方向偏移300米，偏移距離為整道距50米的整數(shù)倍。a、b點(diǎn)沿接收線方向上偏移，c點(diǎn)反向偏移。計(jì)算中間炮點(diǎn)b在允許偏移范圍內(nèi)的地形因子坡度(允許偏移范圍300米、250米、200米、150米、100米、50米，對(duì)應(yīng)整道距位置上地形因子坡度值為46、19、32、18、27、57)使用距離最近，并滿足地形因子坡度限值最小的值27正向偏移100米即為炮點(diǎn)b偏移位置。同向計(jì)算炮點(diǎn)a在允許偏移范圍內(nèi)的地形因子坡度(允許偏移范圍300米、250米、200米、150米、100米、50米，對(duì)于地形因子坡度值為32、46、19、32、18、27)偏移50米位置被b點(diǎn)占用，使用距離最近，并滿足地形因子坡度限值最小的值18正向偏移100米即為炮點(diǎn)a偏移位置。同理c點(diǎn)偏移位置反向偏移100米地形因子坡度值32。
[0047](3)沿垂直接收線方向進(jìn)行炮點(diǎn)偏移。
[0048]例如炮點(diǎn)沿垂直接收線方向上進(jìn)行炮點(diǎn)偏移，計(jì)算炮點(diǎn)在允許偏移范圍內(nèi)的地形因子坡度(允許偏移范圍600米、400米、200米、-200米、-400米、-600米，對(duì)應(yīng)整接收線距位置上，地形因子坡度值為16、18、27、51、34、17)使用距離最近，滿足地形因子坡度限值最小的值27正向偏移200米即為炮點(diǎn)偏移位置。
[0049](4)同一觀測(cè)系統(tǒng)，確定最佳炮點(diǎn)位置進(jìn)行偏移前后的覆蓋次數(shù)對(duì)比。
[0050]分別使用炮點(diǎn)位置偏移前后的觀測(cè)系統(tǒng)，計(jì)算覆蓋次數(shù)進(jìn)行對(duì)比，其中圖6為炮點(diǎn)偏移前的覆蓋次數(shù)分布圖，圖7為炮點(diǎn)偏移后的覆蓋次數(shù)分布圖，可以明確看出，確定最佳炮點(diǎn)位置進(jìn)行偏移后，彌補(bǔ)了覆蓋次數(shù)的不足。
[0051]本發(fā)明實(shí)施例的基于地形因子的炮點(diǎn)偏移方法，可以實(shí)現(xiàn)在地形起伏較大的山地地震采集施工時(shí)，通過(guò)地形因子坡度定量、準(zhǔn)確地計(jì)算出炮點(diǎn)最佳的偏移位置，通過(guò)炮點(diǎn)偏移，確保覆蓋次數(shù)分布滿足采集設(shè)計(jì)要求，得到淺中層資料，從而為三維地震采集觀測(cè)系統(tǒng)變觀設(shè)計(jì)提供參考。
[0052]本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實(shí)施例可提供為方法、系統(tǒng)、或計(jì)算機(jī)程序產(chǎn)品。因此，本發(fā)明可采用完全硬件實(shí)施例、完全軟件實(shí)施例、或結(jié)合軟件和硬件方面的實(shí)施例的形式。而且，本發(fā)明可采用在一個(gè)或多個(gè)其中包含有計(jì)算機(jī)可用程序代碼的計(jì)算機(jī)可用存儲(chǔ)介質(zhì)(包括但不限于磁盤(pán)存儲(chǔ)器、CD-ROM、光學(xué)存儲(chǔ)器等)上實(shí)施的計(jì)算機(jī)程序產(chǎn)品的形式。
[0053]本發(fā)明是參照根據(jù)本發(fā)明實(shí)施例的方法、設(shè)備(系統(tǒng))、和計(jì)算機(jī)程序產(chǎn)品的流程圖和/或方框圖來(lái)描述的。應(yīng)理解可由計(jì)算機(jī)程序指令實(shí)現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合?？商峁┻@些計(jì)算機(jī)程序指令到通用計(jì)算機(jī)、專(zhuān)用計(jì)算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個(gè)機(jī)器，使得通過(guò)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的裝置。
[0054]這些計(jì)算機(jī)程序指令也可存儲(chǔ)在能引導(dǎo)計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備以特定方式工作的計(jì)算機(jī)可讀存儲(chǔ)器中，使得存儲(chǔ)在該計(jì)算機(jī)可讀存儲(chǔ)器中的指令產(chǎn)生包括指令裝置的制造品，該指令裝置實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能。這些計(jì)算機(jī)程序指令也可裝載到計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備上，使得在計(jì)算機(jī)或其他可編程設(shè)備上執(zhí)行一系列操作步驟以產(chǎn)生計(jì)算機(jī)實(shí)現(xiàn)的處理，從而在計(jì)算機(jī)或其他可編程設(shè)備上執(zhí)行的指令提供用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的步驟。
[0055]本發(fā)明中應(yīng)用了具體實(shí)施例對(duì)本發(fā)明的原理及實(shí)施方式進(jìn)行了闡述，以上實(shí)施例的說(shuō)明只是用于幫助理解本發(fā)明的方法及其核心思想；同時(shí)，對(duì)于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在【具體實(shí)施方式】及應(yīng)用范圍上均會(huì)有改變之處，綜上所述，本說(shuō)明書(shū)內(nèi)容不應(yīng)理解為對(duì)本發(fā)明的限制。
【權(quán)利要求】
1.一種基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述炮點(diǎn)偏移方法包括: 根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件； 根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值； 根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)； 在所述需要偏移的炮點(diǎn)的縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移。
2.根據(jù)權(quán)利要求1所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述根據(jù)工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，得到工區(qū)坡度數(shù)字模型文件，包括: 利用遙感技術(shù)獲取所述工區(qū)數(shù)字高程模型； 根據(jù)所述工區(qū)數(shù)字高程模型計(jì)算工區(qū)地形因子坡度值，所述地形因子坡度值S是地形曲面函數(shù)Z=f(x，y)在東西、南北方向上高程變化率的函數(shù)，數(shù)學(xué)表示為:
S = arc(g」fx2+fv2 ； 其中:fx為所述工區(qū)數(shù)字高程模型中X方向南北上的高程變化率；fy為所述工區(qū)數(shù)字高程模型中Y方向東西上的高程變化率； 將計(jì)算得到的所述工區(qū)地形因子坡度值與所述工區(qū)數(shù)字高程模型相結(jié)合，得到所述工區(qū)坡度數(shù)字模型文件。
3.根據(jù)權(quán)利要求2所述 的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值，包括: 所述根據(jù)所述工區(qū)坡度數(shù)字模型文件，使用4點(diǎn)雙線性?xún)?nèi)插法得到每一炮點(diǎn)的坡度值。
4.根據(jù)權(quán)利要求3所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述根據(jù)所述每一炮點(diǎn)的坡度值和設(shè)定的坡度限值，確定需要偏移的炮點(diǎn)，包括: 如果某炮點(diǎn)的坡度值超出設(shè)定的坡度限值，則該炮點(diǎn)需要偏移到平緩地方。
5.根據(jù)權(quán)利要求4所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述在需要偏移的炮點(diǎn)縱向和橫向偏移限制值范圍內(nèi)，選取滿足所述設(shè)定的坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移，包括: 將所述需要偏移的炮點(diǎn)沿接收線方向進(jìn)行縱向偏移，以最淺目的層深度為炮點(diǎn)縱向偏移的最大量，縱向偏移量為道距的整數(shù)倍，計(jì)算每一縱向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
6.根據(jù)權(quán)利要求4所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述在需要偏移的炮點(diǎn)縱向和橫向偏移限制值范圍內(nèi)，選取滿足地形因子坡度限值的位置，對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移，包括: 將所述需要偏移的炮點(diǎn)沿在垂直接收線方向進(jìn)行橫向偏移，以最淺目的層深度為炮點(diǎn)橫向偏移的最大量，橫向偏移量是接收線距的整數(shù)倍，計(jì)算每一橫向偏移位置的地形因子坡度值，在滿足所述設(shè)定的坡度限值的位置中，選擇偏移最小的位置，作為所述需要偏移的炮點(diǎn)的偏移位置。
7.根據(jù)權(quán)利要求5或6所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，當(dāng)所述需要偏移的炮點(diǎn)偏移后的位置與正常炮點(diǎn)重合時(shí)，所述正常炮點(diǎn)也要依次偏移。
8.根據(jù)權(quán)利要求1所述的基于地形因子的炮點(diǎn)偏移方法，其特征在于，所述炮點(diǎn)偏移方法還包括: 在對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移前，計(jì)算所述需要偏移的炮點(diǎn)偏移前的覆蓋次數(shù)，并且對(duì)所述需要偏移的炮點(diǎn)進(jìn)行偏移后，計(jì)算所述需要偏移的炮點(diǎn)偏移后的覆蓋次數(shù)，并進(jìn)行對(duì)比。
【文檔編號(hào)】G01V1/04GK103592676SQ201310507505
【公開(kāi)日】2014年2月19日 申請(qǐng)日期:2013年10月24日 優(yōu)先權(quán)日:2013年10月24日
【發(fā)明者】胡永貴, 倪宇東, 劉麗, 施麗娟, 白志宏, 李偉波, 李平, 李海翔 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司