本申請涉及地球物理勘探技術(shù)領(lǐng)域，特別涉及一種確定近地表速度場的方法及裝置。

背景技術(shù)：

在地球物理勘探中，初至波走時層析是建立近地表速度場的主要方法之一，其建立的近地表速度場可以用于計算層析反演靜校正量，也可以作為疊前深度偏移或全波形反演過程中所采用的初始近地表速度場。在初至波走時層析過程中，可以將目的工區(qū)的近地表區(qū)域的地質(zhì)體模擬為三維網(wǎng)格模型。該網(wǎng)格模型可以包括：主測線方向、聯(lián)絡(luò)測線方向和垂直方向。這三個方向兩兩垂直。其中，主測線方向和聯(lián)絡(luò)測線方向均與地表平行，垂直方向與地表垂直。該網(wǎng)格模型是由三個方向上大小均相同的三維網(wǎng)格構(gòu)成。橫向邊界區(qū)域通常是指該網(wǎng)格模型中在主測線方向或聯(lián)絡(luò)測線方向上初至波射線覆蓋次數(shù)較少的區(qū)域。射線稀疏區(qū)域通常是指非橫向邊界區(qū)域中在垂直方向上穿過網(wǎng)格的初至波射線的條數(shù)較少的區(qū)域。垂直邊界區(qū)域通常是指該網(wǎng)格模型中在垂直方向上穿過每個網(wǎng)格的初至波射線的條數(shù)較少的區(qū)域?；谠摼W(wǎng)格模型，通過初至波走時層析所建立的近地表速度場實際上是網(wǎng)格化的，也就是說，近地表速度場中每一個地層位置處的速度與該網(wǎng)格模型中每一個網(wǎng)格的速度相對應(yīng)。

但是，在實際地震勘探施工過程中，目的工區(qū)中炮檢點的分布，以及從目的工區(qū)的地震數(shù)據(jù)中拾取的初至分布都可能存在一定程度的不均勻，可能導致初至波走時層析的網(wǎng)格模型中穿過每個網(wǎng)格的初至波射線的條數(shù)分布也不均勻。這些實際情況可能導致初至波走時層析的網(wǎng)格模型中較多網(wǎng)格的初至波射線不交叉或初至波射線比較稀疏，進而可能導致初至波走時層析反演過程中得到的基于網(wǎng)格模型的近地表速度場中橫向邊界區(qū)域的網(wǎng)格的速度、射線稀疏區(qū)域的網(wǎng)格的速度以及垂直邊界區(qū)域的網(wǎng)格的速度不合理，準確度較低，不符合目的工區(qū)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。因此，網(wǎng)格模型中邊界區(qū)域的網(wǎng)格的速度以及射線稀疏區(qū)域的網(wǎng)格的速度不合理的問題亟待解決。

技術(shù)實現(xiàn)要素：

本申請實施例的目的是提供一種確定近地表速度場的方法及裝置，以提高初至波走時層析的網(wǎng)格模型中網(wǎng)格的速度的準確度。

為解決上述技術(shù)問題，本申請實施例提供一種確定近地表速度場的方法及裝置是這樣實現(xiàn)的：

一種確定近地表速度場的方法，包括：

獲取目的工區(qū)的初至波數(shù)據(jù)，根據(jù)所述初至波數(shù)據(jù)，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述目的工區(qū)的初始近地表速度場；所述初至波數(shù)據(jù)包括：所述預設(shè)網(wǎng)格模型的初至波射線數(shù)據(jù)和有效炮檢對數(shù)據(jù)；

根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域；

基于所述初始近地表速度場，對所述橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第一修正處理；

基于所述初始近地表速度場，對所述射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度進行第二修正處理；

基于所述初始近地表速度場，對所述垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第三修正處理；

根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，確定所述目的工區(qū)的目標近地表速度場。

優(yōu)選方案中，所述根據(jù)初至波數(shù)據(jù)確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域，包括：

根據(jù)所述初至波數(shù)據(jù)，獲取所述預設(shè)網(wǎng)格模型中有效炮檢對位置處的網(wǎng)格的射線覆蓋次數(shù)；

根據(jù)所述有效炮檢對位置處的網(wǎng)格的射線覆蓋次數(shù)，確定所述預設(shè)網(wǎng)格模型中有效炮檢對的非零覆蓋區(qū)域；

獲取所述預設(shè)網(wǎng)格模型的頂表面區(qū)域；

根據(jù)所述預設(shè)網(wǎng)格模型的頂表面區(qū)域和所述有效炮檢對的非零覆蓋區(qū)域，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。

優(yōu)選方案中，所述根據(jù)預設(shè)網(wǎng)格模型的頂表面區(qū)域和預設(shè)有效炮檢對的非零覆蓋區(qū)域，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域，包括：

沿所述預設(shè)網(wǎng)格模型的主測線方向或聯(lián)絡(luò)測線方向，從所述有效炮檢對的非零覆蓋區(qū)域的兩端向內(nèi)減去所述預設(shè)有效炮檢距的一半，得到所述預設(shè)網(wǎng)格模型的標準非零覆蓋區(qū)域；

將所述預設(shè)網(wǎng)格模型的頂表面區(qū)域減去所述標準非零覆蓋區(qū)域，得到所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。

優(yōu)選方案中，所述根據(jù)初至波數(shù)據(jù)確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域，還包括：

根據(jù)所述初至波數(shù)據(jù)中初至波射線數(shù)據(jù)，獲取所述預設(shè)網(wǎng)格模型的最大覆蓋次數(shù)；

根據(jù)所述預設(shè)網(wǎng)格模型的最大覆蓋次數(shù)，確定所述預設(shè)網(wǎng)格模型中橫向邊界網(wǎng)格；

將所述橫向邊界網(wǎng)格構(gòu)成的區(qū)域作為所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域；所述橫向邊界網(wǎng)格表示所述預設(shè)網(wǎng)格模型的頂表面位置處射線覆蓋次數(shù)小于預設(shè)覆蓋閾值的網(wǎng)格。

優(yōu)選方案中，所述根據(jù)初至波數(shù)據(jù)確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域，包括：

獲取所述預設(shè)網(wǎng)格模型的非橫向邊界區(qū)域；

根據(jù)所述初至波數(shù)據(jù)中初至波射線數(shù)據(jù)，確定所述非橫向邊界區(qū)域中第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值；

根據(jù)預設(shè)分布閾值和所述第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域。

優(yōu)選方案中，所述根據(jù)預設(shè)分布閾值和所述第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域，包括：當所述預設(shè)標準網(wǎng)格的分布值小于所述預設(shè)分布閾值時，將所述第一非橫向邊界位置構(gòu)成的區(qū)域作為射線稀疏區(qū)域。

優(yōu)選方案中，所述根據(jù)預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及初至波數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的目標垂直邊界區(qū)域，包括：

根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置；

根據(jù)所述第一垂直邊界起始位置確定所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置；

將所述第一垂直邊界起始位置和所述第二垂直邊界起始位置作為所述目標垂直邊界起始位置；

將所述目標垂直邊界起始位置至所述預設(shè)網(wǎng)格模型的底表面位置所構(gòu)成的區(qū)域作為所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域。

優(yōu)選方案中，所述根據(jù)預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置，包括：

根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，確定所述預設(shè)網(wǎng)格模型的第一有效區(qū)域；

根據(jù)所述初至波數(shù)據(jù)中初至波射線數(shù)據(jù)，獲取所述第一有效區(qū)域的初至波射線在垂直上穿透的最大有效深度位置；

在所述預設(shè)網(wǎng)格模型中，從所述第一有效區(qū)域的初至波射線在垂直方向上穿透的最大有效深度位置處的網(wǎng)格開始，沿垂直方向向上推至第一個預設(shè)標準網(wǎng)格所處的位置；

將所述第一個預設(shè)標準網(wǎng)格所處的位置作為所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置。

優(yōu)選方案中，所述根據(jù)第一垂直邊界起始位置確定所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置，包括：

根據(jù)所述目的工區(qū)的地震數(shù)據(jù)，獲取所述第一垂直邊界起始位置處的高程值；

獲取所述第一垂直邊界起始位置的水平方向坐標值；

根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，確定所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值；

根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，確定所述射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值；

將所述橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值對應(yīng)的垂直位置，以及所述射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值對應(yīng)的垂直位置作為所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置。

優(yōu)選方案中，所述基于初始近地表速度場，對所述預設(shè)網(wǎng)格模型中橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第一修正處理，包括：

基于所述初始近地表速度場，獲取所述預設(shè)網(wǎng)格模型中第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度；

基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第二重疊區(qū)域的待處理網(wǎng)格的速度進行第一修正處理。

優(yōu)選方案中，所述基于第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第二重疊區(qū)域的待處理網(wǎng)格的速度進行第一修正處理，包括：

獲取所述預設(shè)網(wǎng)格模型中第一長方體區(qū)域；

基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，獲取所述第一長方體區(qū)域內(nèi)第一標準網(wǎng)格的個數(shù)以及第一標準網(wǎng)格的速度；

當所述第一標準網(wǎng)格的個數(shù)大于或等于預設(shè)第一網(wǎng)格個數(shù)時，基于所述第一標準網(wǎng)格的速度，采用反距離加權(quán)的方式替換所述第二重疊區(qū)域內(nèi)待處理網(wǎng)格的速度；或，當所述第一標準網(wǎng)格的個數(shù)小于預設(shè)第一網(wǎng)格個數(shù)時，分別將所述第一長方體區(qū)域沿主測線方向和聯(lián)絡(luò)測線方向上的長度增加一倍，直至滿足所述第一標準網(wǎng)格的個數(shù)大于或等于預設(shè)第一網(wǎng)格個數(shù)的條件即可。

優(yōu)選方案中，基于初始近地表速度場，對所述射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度進行第二修正處理，包括：

基于所述初始近地表速度場，獲取所述預設(shè)網(wǎng)格模型中第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度；

基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第三重疊區(qū)域的待處理網(wǎng)格的速度進行第二修正處理。

優(yōu)選方案中，所述基于第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第三重疊區(qū)域的待處理網(wǎng)格的速度進行第二修正處理，包括：

獲取所述預設(shè)網(wǎng)格模型中第二長方體區(qū)域；

基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，獲取所述第二長方體區(qū)域內(nèi)第二標準網(wǎng)格的個數(shù)以及第二標準網(wǎng)格的速度；

當所述第二標準網(wǎng)格的個數(shù)大于或等于預設(shè)第二網(wǎng)格個數(shù)時，基于所述第二標準網(wǎng)格的速度，采用反距離加權(quán)的方式替換所述第三重疊區(qū)域中待處理網(wǎng)格的速度；或，當所述第二標準網(wǎng)格的個數(shù)小于預設(shè)第二網(wǎng)格個數(shù)時，分別將所述第二長方體區(qū)域沿主測線方向和聯(lián)絡(luò)測線方向上的長度增加一倍，直至滿足所述第二標準網(wǎng)格的個數(shù)大于或等于預設(shè)第二網(wǎng)格個數(shù)的條件即可。

優(yōu)選方案中，所述基于初始近地表速度場，對所述垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第三修正處理，包括：

根據(jù)所述初始近地表速度場和所述垂直邊界區(qū)域中目標垂直邊界起始位置，確定所述目標垂直邊界起始位置處的網(wǎng)格的速度；

對所述目標垂直邊界初始位置處的網(wǎng)格的速度進行第一平滑處理，得到所述目標垂直邊界初始位置處的網(wǎng)格的目標速度；

基于所述目標垂直邊界初始位置處的網(wǎng)格的目標速度，確定所述垂直邊界區(qū)域中目標垂直邊界位置處的網(wǎng)格的目標速度；

將所述近地表速度場中目標垂直邊界位置處的網(wǎng)格的速度替換為所述目標垂直邊界位置處的網(wǎng)格的目標速度。

優(yōu)選方案中，所述根據(jù)初始近地表速度場和垂直邊界區(qū)域中目標垂直邊界起始位置，確定所述目標縱向邊界起始位置處的網(wǎng)格的速度，包括：

根據(jù)所述目標垂直邊界起始位置中第一垂直邊界起始位置和所述近地表速度場，確定所述第一垂直邊界起始位置處的網(wǎng)格的速度；

根據(jù)所述第一垂直邊界起始位置處的網(wǎng)格的速度，采用反距離加權(quán)的方式對所述目標垂直邊界起始位置中第二垂直邊界起始位置進行插值處理，確定所述第二垂直邊界起始位置處的網(wǎng)格的速度；

將所述第一垂直邊界起始位置的速度和所述第二垂直邊界起始位置的速度作為所述目標垂直邊界起始位置的速度。

優(yōu)選方案中，所述根據(jù)第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，確定所述目的工區(qū)的目標近地表速度場，包括：

根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，得到三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度；

對所述三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度進行第二平滑處理，得到所述預設(shè)網(wǎng)格模型的網(wǎng)格的目標速度；

根據(jù)所述預設(shè)網(wǎng)格模型的網(wǎng)格的目標速度，確定所述目的工區(qū)的目標近地表速度場。

優(yōu)選方案中，采用下述公式對所述三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度進行第二平滑處理：

公式中，vgnew表示所述預設(shè)網(wǎng)格模型中某一網(wǎng)格的目標速度，vgi表示所述第二平滑處理前以該網(wǎng)格為中心、以所述預設(shè)網(wǎng)格模型在主測線方向上ni為半徑的區(qū)域內(nèi)任一網(wǎng)格的速度，vgx表示所述第二平滑處理前以該網(wǎng)格為中心、以所述預設(shè)網(wǎng)格模型在聯(lián)絡(luò)測線方向上nx為半徑的區(qū)域內(nèi)任一網(wǎng)格的速度，vgz表示所述第二平滑處理前以該網(wǎng)格為中心、所述預設(shè)網(wǎng)格模型在垂直方向上長度為三個網(wǎng)格所構(gòu)成的區(qū)域內(nèi)任一網(wǎng)格的速度，ni表示所述預設(shè)網(wǎng)格模型中主測線方向上的平滑半徑，nx表示所述預設(shè)網(wǎng)格模型中聯(lián)絡(luò)測線方向上的平滑半徑。

一種確定近地表速度場的裝置，所述裝置包括：橫向邊界區(qū)域確定模塊、垂直邊界區(qū)域確定模塊、第一修正處理模塊、第二修正處理模塊、第三修正處理模塊和目標近地表速度場確定模塊；其中，

所述橫向邊界區(qū)域確定模塊，用于獲取目的工區(qū)的初至波數(shù)據(jù)，根據(jù)所述初至波數(shù)據(jù)，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述目的工區(qū)的初始近地表速度場；所述初至波數(shù)據(jù)包括：所述預設(shè)網(wǎng)格模型的初至波射線數(shù)據(jù)和有效炮檢對數(shù)據(jù)；

所述垂直邊界區(qū)域確定模塊，用于根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域；

所述第一修正處理模塊，用于基于所述初始近地表速度場，對所述橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第一修正處理；

所述第二修正處理模塊，用于基于所述初始近地表速度場，對所述射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度進行第二修正處理；

所述第三修正處理模塊，用于基于所述初始近地表速度場，對所述垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第三修正處理；

所述目標近地表速度場確定模塊，用于根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，確定所述目的工區(qū)的目標近地表速度場。

本申請實施例提供一種確定近地表速度場的方法，基于預設(shè)網(wǎng)格模型的初至波數(shù)據(jù)中初至波射線數(shù)據(jù)和有效炮檢對位置，合理地確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域、射線稀疏區(qū)域和縱向邊界區(qū)域。針對所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域、射線稀疏區(qū)域和縱向邊界區(qū)域的網(wǎng)格的速度，分別進行合理地修正處理，從而可以得到較為合理的近地表速度場，其準確度較高，且更加符合目的工區(qū)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

附圖說明

為了更清楚地說明本申請實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本申請中記載的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1是本申請一種確定近地表速度場的方法實施例的流程圖；

圖2是本申請實施例中橫向邊界區(qū)域劃分后的預設(shè)網(wǎng)格模型的俯視平面示意圖；

圖3是本申請實施例中采用本申請方法進行處理前的近地表速度場的剖面示意圖；

圖4是本申請實施例中采用本申請方法進行處理后的近地表速度場的剖面示意圖；

圖5是本申請實施例中基于處理前的近地表速度場進行正演處理得到的初至波射線密度的剖面示意圖；

圖6是本申請實施例中基于處理后的近地表速度場進行正演處理得到的初至波射線密度的剖面示意圖；

圖7是本申請確定近地表速度場的裝置實施例的組成結(jié)構(gòu)圖；

圖8是本申請確定近地表速度場的裝置實施例中垂直邊界區(qū)域確定模塊的組成結(jié)構(gòu)圖；

圖9是本申請確定近地表速度場的裝置實施例中第三修正處理模塊的組成結(jié)構(gòu)圖。

具體實施方式

本申請實施例提供一種確定近地表速度場的方法及裝置。

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請中的技術(shù)方案，下面將結(jié)合本申請實施例中的附圖，對本申請實施例中的技術(shù)方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本申請一部分實施例，而不是全部的實施例?；诒旧暾堉械膶嵤├?，本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都應(yīng)當屬于本申請保護的范圍。

圖1是本申請一種確定近地表速度場的方法實施例的流程圖。如圖1所示，所述確定近地表速度場的方法，包括以下步驟。

步驟s101：獲取目的工區(qū)的初至波數(shù)據(jù)，根據(jù)所述初至波數(shù)據(jù)，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述目的工區(qū)的初始近地表速度場。

初至波是指在所述目的工區(qū)的地震數(shù)據(jù)中，從目的工區(qū)的地表上的炮點出發(fā)，最先到達目的工區(qū)的地表上的檢波點的地震波。

在初至波走時層析過程中，可以將所述目的工區(qū)的近地表區(qū)域的地質(zhì)體模擬為三維網(wǎng)格模型，即所述預設(shè)網(wǎng)格模型。所述近地表區(qū)域可以表示所述目的工區(qū)的地震數(shù)據(jù)中0～200毫秒之間的采樣時間對應(yīng)的地質(zhì)區(qū)域。所述預設(shè)網(wǎng)格模型可以包括：主測線方向、聯(lián)絡(luò)測線方向和垂直方向。所述三個方向兩兩垂直。其中，所述主測線方向和所述聯(lián)絡(luò)測線方向均與地表平行，所述垂直方向與地表垂直。所述預設(shè)網(wǎng)格模型是由三個方向上大小均相同的三維網(wǎng)格構(gòu)成。所述預設(shè)網(wǎng)格模型的頂表面與所述目的工區(qū)的地表面相對應(yīng)。所述預設(shè)網(wǎng)格模型的底表面位置與所述目的工區(qū)的近地表區(qū)域的最深地層位置相對應(yīng)。

所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域是指所述預設(shè)網(wǎng)格模型的頂表面的主測線方向或聯(lián)絡(luò)測線方向上的區(qū)域。所述預設(shè)網(wǎng)格模型的射線稀疏區(qū)域也是指所述預設(shè)網(wǎng)格模型的頂表面的主測線方向或聯(lián)絡(luò)測線方向上的區(qū)域。

所述初至波數(shù)據(jù)可以包括：所述預設(shè)網(wǎng)格模型的初至波射線數(shù)據(jù)和有效炮檢對位置。所述有效炮檢對位置可以表示預設(shè)有效炮檢對的中點在所述預設(shè)網(wǎng)格模型中的網(wǎng)格位置。所述預設(shè)有效炮檢對可以表示所述目的工區(qū)中預設(shè)有效炮檢距的炮檢對。所述預設(shè)有效炮檢距可以為2000米。所述初至波射線數(shù)據(jù)可以包括：初至波射線條數(shù)、初至波射線的最大有效深度位置，以及初至波射線覆蓋次數(shù)。所述初至波射線條數(shù)可以表示所述預設(shè)網(wǎng)格模型中在垂直方向上穿過每個網(wǎng)格的初至波射線的條數(shù)。所述初至波射線的最大有效深度位置可以表示所述預設(shè)網(wǎng)格模型的初至波射線在垂直方向上穿透的最大有效深度位置。所述初至波射線覆蓋次數(shù)可以表示所述預設(shè)網(wǎng)格模型的初至波射線在主測線方向或聯(lián)絡(luò)測線方向上覆蓋所述預設(shè)網(wǎng)格模型的頂表面位置處的網(wǎng)格的次數(shù)。

所述根據(jù)初至波數(shù)據(jù)確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。具體地，根據(jù)所述初至波數(shù)據(jù)，可以獲取所述預設(shè)網(wǎng)格模型中有效炮檢對位置處的網(wǎng)格的射線覆蓋次數(shù)。根據(jù)所述有效炮檢對位置處的網(wǎng)格的初至波射線覆蓋次數(shù)，可以確定所述預設(shè)網(wǎng)格模型中有效炮檢對的非零覆蓋區(qū)域。所述有效炮檢對的非零覆蓋區(qū)域可以表示所述有效炮檢對位置處的初至波射線覆蓋次數(shù)大于零的網(wǎng)格構(gòu)成的區(qū)域?？梢垣@取所述預設(shè)網(wǎng)格模型的頂表面區(qū)域。根據(jù)所述預設(shè)網(wǎng)格模型的頂表面區(qū)域和所述有效炮檢對的非零覆蓋區(qū)域，可以確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。

進一步地，所述根據(jù)預設(shè)網(wǎng)格模型的頂表面區(qū)域和預設(shè)有效炮檢對的非零覆蓋區(qū)域，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。具體地，沿所述預設(shè)網(wǎng)格模型的頂表面的主測線方向或聯(lián)絡(luò)測線方向，從所述有效炮檢對的非零覆蓋區(qū)域的兩端向內(nèi)減去所述預設(shè)有效炮檢距的一半，可以得到所述預設(shè)網(wǎng)格模型的標準非零覆蓋區(qū)域。將所述預設(shè)網(wǎng)格模型的頂表面區(qū)域減去所述標準非零覆蓋區(qū)域，可以得到所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。

所述根據(jù)初至波數(shù)據(jù)，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。在另一種實施方式中，具體地，根據(jù)所述初至波數(shù)據(jù)中初至波射線數(shù)據(jù)，可以獲取所述預設(shè)網(wǎng)格模型的最大覆蓋次數(shù)。所述預設(shè)網(wǎng)格模型的最大覆蓋次數(shù)可以表示所述預設(shè)網(wǎng)格模型的頂表面位置處的網(wǎng)格的初至波射線覆蓋次數(shù)的最大值。根據(jù)所述預設(shè)網(wǎng)格模型的最大覆蓋次數(shù)，可以確定所述預設(shè)網(wǎng)格模型中橫向邊界網(wǎng)格。所述橫向邊界網(wǎng)格可以表示所述預設(shè)網(wǎng)格模型的頂表面位置處的初至波射線覆蓋次數(shù)小于預設(shè)覆蓋閾值的網(wǎng)格。所述預設(shè)覆蓋閾值可以為所述最大覆蓋次數(shù)的百分之五十?？梢詫⑺鰴M向邊界網(wǎng)格構(gòu)成的區(qū)域作為所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域。

所述根據(jù)初至波數(shù)據(jù)確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域。具體地，可以獲取所述預設(shè)網(wǎng)格模型的非橫向邊界區(qū)域。所述非橫向邊界區(qū)域可以為所述預設(shè)網(wǎng)格模型的頂表面上的橫向邊界區(qū)域以外的區(qū)域。根據(jù)所述初至波數(shù)據(jù)中初至波射線數(shù)據(jù)，可以確定所述非橫向邊界區(qū)域中第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值。所述第一非橫向邊界位置可以表示所述非橫向邊界區(qū)域中任一非橫向邊界位置。所述預設(shè)標準網(wǎng)格可以為所述預設(shè)網(wǎng)絡(luò)模型中初至波射線的條數(shù)大于預設(shè)射線閾值的網(wǎng)格。所述預設(shè)射線閾值可以為0.15乘以所述預設(shè)網(wǎng)格模型中在垂直方向上穿過每個網(wǎng)格的初至波射線的條數(shù)的最大值。所述預設(shè)標準網(wǎng)格的分布值可以表示所述非橫向邊界區(qū)域中每一非橫向邊界位置在垂直方向上的預設(shè)標準網(wǎng)格的個數(shù)。根據(jù)預設(shè)分布閾值和所述第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值，可以確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域。所述預設(shè)分布閾值可以為所述預設(shè)標準網(wǎng)格分布值的平均值的三分之一。

進一步地，所述根據(jù)預設(shè)分布閾值和所述第一非橫向邊界位置處的預設(shè)標準網(wǎng)格的分布值，確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的射線稀疏區(qū)域，可以包括：當所述預設(shè)標準網(wǎng)格的分布值小于所述預設(shè)分布閾值時，可以將所述第一非橫向邊界位置構(gòu)成的區(qū)域作為射線稀疏區(qū)域。

例如，圖2是本申請實施例中橫向邊界區(qū)域劃分后的預設(shè)網(wǎng)格模型的俯視平面示意圖。如圖2所示，圖2中灰色區(qū)域為所述橫向邊界區(qū)域，圖2中淺灰色點填充區(qū)域為所述非橫向邊界區(qū)域且非射線稀疏區(qū)域，圖2中深灰色斜線填充區(qū)域為所述射線稀疏區(qū)域。

所述根據(jù)初至波數(shù)據(jù)，確定所述目的工區(qū)的初始近地表速度場。具體地，可以獲取所述目的工區(qū)的地震數(shù)據(jù)。根據(jù)所述初至波數(shù)據(jù)，對所述地震數(shù)據(jù)進行層析反演靜校正處理，可以確定所述目的工區(qū)的初始近地表速度場。

對所述目的工區(qū)進行地震勘探和數(shù)據(jù)采集，可以獲取所述目的工區(qū)的地震數(shù)據(jù)和初至波數(shù)據(jù)。

步驟s102：根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，可以確定所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域。

具體地，根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，可以確定所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置。根據(jù)所述第一垂直邊界起始位置可以確定所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置?？梢詫⑺龅谝淮怪边吔缙鹗嘉恢煤退龅诙怪边吔缙鹗嘉恢米鳛樗瞿繕舜怪边吔缙鹗嘉恢谩？梢詫⑺瞿繕舜怪边吔缙鹗嘉恢弥了鲱A設(shè)網(wǎng)格模型的底表面位置所構(gòu)成的區(qū)域作為所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域。

所述根據(jù)預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置。具體地，根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，可以確定所述預設(shè)網(wǎng)格模型的第一有效區(qū)域。所述第一有效區(qū)域可以表示所述預設(shè)網(wǎng)格模型的頂表面上的非橫向邊界區(qū)域且非射線稀疏區(qū)域。所述非射線稀疏區(qū)域可以表示所述非橫向邊界區(qū)域中除射線稀疏區(qū)域以外的區(qū)域。根據(jù)所述初至波射線數(shù)據(jù)，可以獲取所述第一有效區(qū)域的初至波射線在垂直方向上穿透的最大有效深度位置。在所述預設(shè)網(wǎng)格模型中，可以從所述第一有效區(qū)域的初至波射線在垂直方向上穿透的最大有效深度位置處的網(wǎng)格開始，沿垂直方向向上推至第一個預設(shè)標準網(wǎng)格所處的位置?？梢詫⑺龅谝粋€預設(shè)標準網(wǎng)格所處的位置作為所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置。

所述根據(jù)第一垂直邊界起始位置確定所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置。具體地，根據(jù)所述目的工區(qū)的地震數(shù)據(jù)，可以獲取所述第一垂直邊界起始位置處的高程值。可以獲取所述第一垂直邊界起始位置的水平方向坐標值。所述水平方向坐標值可以表示一組主測線方向坐標值和聯(lián)絡(luò)測線方向坐標值。一組水平方向坐標值可以與一個橫向位置相對應(yīng)。根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，可以確定所述橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值。根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，可以確定所述射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值。可以將所述橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值對應(yīng)的垂直位置，以及所述射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值對應(yīng)的垂直位置作為所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置。所述高程值可以與所述預設(shè)網(wǎng)格模型中垂直位置相對應(yīng)。

進一步地，根據(jù)第一垂直邊界起始位置處的高程值，以及第一垂直邊界起始位置的水平方向坐標值，確定所述橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值，可以包括：根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，可以采用反距離加權(quán)的方式對所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域內(nèi)每一個橫向位置進行插值處理，可以得到所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域內(nèi)每一個橫向位置處的高程值。

進一步地，根據(jù)第一垂直邊界起始位置處的高程值，以及第一垂直邊界起始位置的水平方向坐標值，可以確定所述射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值，可以包括：根據(jù)所述第一垂直邊界起始位置處的高程值，以及所述第一垂直邊界起始位置的水平方向坐標值，可以采用反距離加權(quán)的方式對所述預設(shè)網(wǎng)格模型的射線稀疏區(qū)域內(nèi)每一個橫向位置進行插值處理，可以得到所述預設(shè)網(wǎng)格模型的射線稀疏區(qū)域內(nèi)每一個橫向位置處的高程值。

步驟s103：基于所述初始近地表速度場，對所述橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第一修正處理。

具體地，基于所述初始近地表速度場，可以獲取所述預設(shè)網(wǎng)格模型中第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度。所述第一重疊區(qū)域可以表示從所述預設(shè)網(wǎng)格模型的頂表面至所述目標垂直邊界起始位置所構(gòu)成的區(qū)域與所述第一有效區(qū)域重疊的區(qū)域。基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，可以對所述預設(shè)網(wǎng)格模型中第二重疊區(qū)域的待處理網(wǎng)格的速度進行第一修正處理。所述第二重疊區(qū)域可以表示從所述預設(shè)網(wǎng)格模型的頂表面至所述目標垂直邊界起始位置所構(gòu)成的區(qū)域與所述橫向邊界區(qū)域重疊的區(qū)域。所述第二重疊區(qū)域的待處理網(wǎng)格可以表示所述第二重疊區(qū)域內(nèi)射線條數(shù)小于所述預設(shè)射線閾值的網(wǎng)格。

進一步地，所述基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第二重疊區(qū)域的待處理網(wǎng)格的速度進行第一修正處理。具體地，可以獲取所述預設(shè)網(wǎng)格模型中第一長方體區(qū)域?；谒龅谝恢丿B區(qū)域的預設(shè)標準網(wǎng)格的速度，可以獲取所述第一長方體區(qū)域內(nèi)第一標準網(wǎng)格的個數(shù)以及第一標準網(wǎng)格的速度。當所述第一標準網(wǎng)格的個數(shù)大于或等于預設(shè)第一網(wǎng)格個數(shù)時，基于所述第一標準網(wǎng)格的速度，可以采用反距離加權(quán)的方式替換所述第二重疊區(qū)域中待處理網(wǎng)格的速度；或，當所述第一標準網(wǎng)格的個數(shù)小于預設(shè)第一網(wǎng)格個數(shù)時，可以分別將所述第一長方體區(qū)域沿主測線方向和聯(lián)絡(luò)測線方向上的長度增加一倍，直至滿足所述第一標準網(wǎng)格的個數(shù)大于或等于預設(shè)第一網(wǎng)格個數(shù)的條件即可。所述第一標準網(wǎng)格可以表示在所述第一長方體區(qū)域內(nèi)且在所述第一重疊區(qū)域內(nèi)的預設(shè)標準網(wǎng)格。所述預設(shè)第一網(wǎng)格個數(shù)可以為11。所述第一長方體區(qū)域可以表示以所述第二重疊區(qū)域內(nèi)的待處理網(wǎng)格為中心、主測線方向上的長度的一半為預設(shè)有效炮檢距、聯(lián)絡(luò)測線方向上的長度的一半為預設(shè)有效炮檢距、以及垂直方向上長度為連續(xù)三個網(wǎng)格所構(gòu)成的區(qū)域。

經(jīng)過第一修正處理后，所述預設(shè)網(wǎng)格模型中橫向邊界區(qū)域的網(wǎng)格的速度較為合理，準確度較高，較符合該區(qū)域?qū)?yīng)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

步驟s104：基于所述初始近地表速度場，對所述射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度進行第二修正處理。

具體地，基于所述初始近地表速度場，可以獲取所述預設(shè)網(wǎng)格模型中第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度。所述第一重疊區(qū)域可以表示從所述預設(shè)網(wǎng)格模型的頂表面至所述目標垂直邊界起始位置所構(gòu)成的區(qū)域與所述第一有效區(qū)域重疊的區(qū)域?；谒龅谝恢丿B區(qū)域的預設(shè)標準網(wǎng)格的速度，可以對所述預設(shè)網(wǎng)格模型中第三重疊區(qū)域的待處理網(wǎng)格的速度進行第二修正處理。所述第三重疊區(qū)域可以表示從所述預設(shè)網(wǎng)格模型的頂表面至所述目標垂直邊界起始位置所構(gòu)成的區(qū)域與所述射線稀疏區(qū)域重疊的區(qū)域。所述第三重疊區(qū)域內(nèi)的待處理網(wǎng)格可以表示所述第三重疊區(qū)域內(nèi)射線條數(shù)小于所述預設(shè)射線閾值的網(wǎng)格。

進一步地，所述基于所述第一重疊區(qū)域的預設(shè)標準網(wǎng)格的速度，對所述預設(shè)網(wǎng)格模型中第三重疊區(qū)域的待處理網(wǎng)格的速度進行第二修正處理。具體地，可以獲取所述預設(shè)網(wǎng)格模型中第二長方體區(qū)域?；谒龅谝恢丿B區(qū)域的預設(shè)標準網(wǎng)格的速度，可以獲取所述第二長方體區(qū)域內(nèi)第二標準網(wǎng)格的個數(shù)以及第二標準網(wǎng)格的速度。當所述第二標準網(wǎng)格的個數(shù)大于或等于預設(shè)第二網(wǎng)格個數(shù)時，基于所述第二標準網(wǎng)格的速度，可以采用反距離加權(quán)的方式替換所述第三重疊區(qū)域內(nèi)待處理網(wǎng)格的速度；或，當所述第二標準網(wǎng)格的個數(shù)小于預設(shè)第二網(wǎng)格個數(shù)時，可以分別將所述第二長方體區(qū)域沿主測線方向和聯(lián)絡(luò)測線方向上的長度增加一倍，直至滿足所述第二標準網(wǎng)格的個數(shù)大于或等于預設(shè)第二網(wǎng)格個數(shù)的條件即可。所述第二標準網(wǎng)格可以表示在所述第二長方體區(qū)域內(nèi)且在所述第一重疊區(qū)域內(nèi)的預設(shè)標準網(wǎng)格。所述預設(shè)第二網(wǎng)格個數(shù)可以為9。

所述第二長方體區(qū)域可以表示以所述第三重疊區(qū)域內(nèi)的待處理網(wǎng)格為中心、主測線方向上的長度為預設(shè)有效炮檢距、聯(lián)絡(luò)測線方向上的長度為預設(shè)有效炮檢距、以及垂直方向上長度為連續(xù)三個網(wǎng)格所構(gòu)成的區(qū)域。

經(jīng)過第二修正處理后，所述預設(shè)網(wǎng)格模型中射線稀疏區(qū)域的網(wǎng)格的速度較為合理，準確度較高，較符合該區(qū)域?qū)?yīng)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

需要說明的是，步驟s104可以在步驟s103之前或之后，本申請對此并不作出限定。

步驟s105：基于所述初始近地表速度場，對所述垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第三修正處理。

具體地，根據(jù)初始近地表速度場和所述垂直邊界區(qū)域中目標垂直邊界起始位置，可以確定所述目標垂直邊界起始位置處的網(wǎng)格的速度?？梢詫λ瞿繕舜怪边吔绯跏嘉恢锰幍木W(wǎng)格的速度進行第一平滑處理，得到所述目標垂直邊界初始位置處的網(wǎng)格的目標速度。基于所述目標垂直邊界初始位置處的網(wǎng)格的目標速度，可以確定所述垂直邊界區(qū)域中目標垂直邊界位置處的網(wǎng)格的目標速度。所述目標垂直邊界位置處的網(wǎng)格可以表示所述預設(shè)網(wǎng)格模型中所述目標垂直邊界起始位置以下至所述預設(shè)網(wǎng)格模型的底表面位置處的網(wǎng)格?？梢詫⑺鼋乇硭俣葓鲋心繕舜怪边吔缥恢锰幍木W(wǎng)格的速度替換為所述目標垂直邊界位置處的網(wǎng)格的目標速度。

所述根據(jù)初始近地表速度場和垂直邊界區(qū)域中目標垂直邊界起始位置，確定所述目標垂直邊界起始位置處的網(wǎng)格的速度。具體地，根據(jù)所述目標垂直邊界起始位置中第一垂直邊界起始位置和所述近地表速度場，可以確定所述第一垂直邊界起始位置處的網(wǎng)格的速度。根據(jù)所述第一垂直邊界起始位置處的網(wǎng)格的速度，采用反距離加權(quán)的方式對所述目標垂直邊界起始位置中第二垂直邊界起始位置進行插值處理，可以確定所述第二垂直邊界起始位置處的網(wǎng)格的速度?？梢詫⑺龅谝淮怪边吔缙鹗嘉恢玫乃俣群退龅诙怪边吔缙鹗嘉恢玫乃俣茸鳛樗瞿繕舜怪边吔缙鹗嘉恢玫乃俣?。

所述根據(jù)第一垂直邊界起始位置和近地表速度場，可以確定所述第一垂直邊界起始位置處的網(wǎng)格的速度。具體地，根據(jù)所述近地表速度場，可以確定所述預設(shè)網(wǎng)格模型中每一個網(wǎng)格的速度?；谒鲱A設(shè)網(wǎng)格模型中每一個網(wǎng)格的速度，可以獲取第一垂直邊界起始位置處的網(wǎng)格的速度。

進一步地，可以采用下述公式對所述目標垂直邊界初始位置處的網(wǎng)格的速度進行第一平滑處理：

公式(1)中，vgnew表示所述目標垂直邊界初始位置處的網(wǎng)格的目標速度，vgi表示所述目標垂直邊界初始位置在主測線方向上的網(wǎng)格的速度，vgx表示所述目標垂直邊界初始位置在聯(lián)絡(luò)測線方向上的網(wǎng)格的速度，ni表示所述預設(shè)網(wǎng)格模型中主測線方向上的平滑半徑，nx表示所述預設(shè)網(wǎng)格模型中聯(lián)絡(luò)測線方向上的平滑半徑。所述主測線方向上的平滑半徑和聯(lián)絡(luò)測線方向上的平滑半徑均可以所述預設(shè)有效炮檢距。

進一步地，可以采用下述公式確定所述垂直邊界區(qū)域中目標垂直邊界位置處的網(wǎng)格的速度：

vki＝vk+(ki-k)*dz*factor(2)

公式(2)中，vki表示所述目標垂直邊界位置處的第ki個網(wǎng)格的目標速度，vk表示所述垂直邊界初始位置處的第k個網(wǎng)格的目標速度，dz表示所述預設(shè)網(wǎng)格模型中垂直方向上相鄰網(wǎng)格的距離，factor表示梯度因子。所述梯度因子的取值范圍可以為0.001～0.1。所述梯度因子的物理意義可以是所述預設(shè)網(wǎng)格模型中在垂直方向上深度每增加一米，該深度位置處的網(wǎng)格的速度可以增加一定值。

經(jīng)過第三修正處理后，所述預設(shè)網(wǎng)格模型中垂直邊界區(qū)域的網(wǎng)格的速度較為合理，準確度較高，較符合該區(qū)域?qū)?yīng)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

需要說明的是，步驟s105可以在步驟s103之前或之后，也可以在步驟s104之前或之后，本申請對此并不作出限定。

步驟s106：根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，確定所述目的工區(qū)的目標近地表速度場。

具體地，根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，可以得到三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度。對所述三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度進行第二平滑處理，可以得到所述預設(shè)網(wǎng)格模型的網(wǎng)格的目標速度。根據(jù)所述預設(shè)網(wǎng)格模型的網(wǎng)格的目標速度，可以確定所述目的工區(qū)的目標近地表速度場。

進一步地，可以采用下述公式對所述三次修正處理后的預設(shè)網(wǎng)格模型的網(wǎng)格的速度進行第二平滑處理：

公式(3)中，vgnew表示所述預設(shè)網(wǎng)格模型中某一網(wǎng)格的目標速度，vgi表示所述第二平滑處理前以該網(wǎng)格為中心、以所述預設(shè)網(wǎng)格模型在主測線方向上ni為半徑的區(qū)域內(nèi)任一網(wǎng)格的速度，vgx表示所述第二平滑處理前以該網(wǎng)格為中心、以所述預設(shè)網(wǎng)格模型在聯(lián)絡(luò)測線方向上nx為半徑的區(qū)域內(nèi)任一網(wǎng)格的速度，vgz表示所述第二平滑處理前以該網(wǎng)格為中心、所述預設(shè)網(wǎng)格模型在垂直方向上長度為三個網(wǎng)格所構(gòu)成的區(qū)域內(nèi)任一網(wǎng)格的速度，ni表示所述預設(shè)網(wǎng)格模型中主測線方向上的平滑半徑，nx表示所述預設(shè)網(wǎng)格模型中聯(lián)絡(luò)測線方向上的平滑半徑。

在所述第二平滑處理過程中，當所處的網(wǎng)格為所述橫向邊界區(qū)域或所述垂直邊界區(qū)域內(nèi)的網(wǎng)格時，所述ni和nx的取值可以均為所述預設(shè)有效炮檢距。當所處的網(wǎng)格為非橫向邊界區(qū)域且非垂直邊界區(qū)域內(nèi)的網(wǎng)格時，所述ni和nx的取值可以均為所述預設(shè)有效炮檢距的一半。所述非橫向邊界區(qū)域且非垂直邊界區(qū)域可以為所述預設(shè)網(wǎng)格模型中除橫向邊界區(qū)域和垂直邊界區(qū)域以外的區(qū)域。

圖3是本申請實施例中采用本申請方法進行處理前的近地表速度場的剖面示意圖。圖4是本申請實施例中采用本申請方法進行處理后的近地表速度場的剖面示意圖。圖3和圖4的橫坐標均表示所述目的工區(qū)的主測線方向的位置，圖3和圖4的縱坐標均表示所述目的工區(qū)的地層海拔高度，圖3和圖4中的灰度值分別表示采用本申請方法進行處理前和處理后的近地表速度場。如圖3和圖4所示，相比所述處理前的近地表速度場，所述處理后的近地表速度場的兩端邊界位置的速度的連續(xù)性較好，更符合所述目的工區(qū)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

圖5是本申請實施例中基于處理前的近地表速度場進行正演處理得到的初至波射線密度的剖面示意圖。圖6是本申請實施例中基于處理后的近地表速度場進行正演處理得到的初至波射線密度的剖面示意圖。圖5和圖6的橫坐標均表示所述目的工區(qū)的主測線方向的位置，圖5和圖6的縱坐標均表示所述目的工區(qū)的地層海拔高度，圖5中的灰度值表示基于處理前近地表速度場進行正演處理得到的初至波射線密度。圖6中的灰度值表示基于處理后近地表速度場進行正演處理得到的初至波射線密度。如圖5和圖6所示，在所述預設(shè)網(wǎng)格模型的深層位置處，相比基于處理前的近地表速度場進行正演處理得到的初至波射線密度，基于處理后的近地表速度場進行正演處理得到的初至波射線密度更加收斂和穩(wěn)定，異常的射線密度明顯減少，表明邊界處理后的近地表速度場的準確度較高。

在另一種實施方式中，針對所述目的工區(qū)的地震數(shù)據(jù)在初至波走時層析反演迭代過程中所述得到的近地表速度場，采用本申請的方法均可以其進行處理，可以提高該近地表速度場的準確度，從而可以符合目的工區(qū)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

所述確定近地表速度場的方法實施例，基于預設(shè)網(wǎng)格模型的初至波數(shù)據(jù)中初至波射線數(shù)據(jù)和有效炮檢對位置，合理地確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域、射線稀疏區(qū)域和垂直邊界區(qū)域。針對所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域、射線稀疏區(qū)域和垂直邊界區(qū)域的網(wǎng)格的速度，分別進行合理地修正處理，從而可以得到較為合理的近地表速度場，其準確度較高，且更加符合目的工區(qū)的真實近地表的地質(zhì)結(jié)構(gòu)變化特征。

圖7是本申請確定近地表速度場的裝置實施例的組成結(jié)構(gòu)圖。如圖7所示，所述確定近地表速度場的裝置，可以包括：橫向邊界區(qū)域確定模塊100、垂直邊界區(qū)域確定模塊200、第一修正處理模塊300、第二修正處理模塊400、第三修正處理模塊500和目標近地表速度場確定模塊600。

所述橫向邊界區(qū)域確定模塊100，可以用于獲取目的工區(qū)的初至波數(shù)據(jù)，根據(jù)所述初至波數(shù)據(jù)，可以確定所述目的工區(qū)的預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述目的工區(qū)的初始近地表速度場。所述初至波數(shù)據(jù)可以包括：所述預設(shè)網(wǎng)格模型的初至波射線數(shù)據(jù)和有效炮檢對數(shù)據(jù)。

所述垂直邊界區(qū)域確定模塊200，可以用于根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域。

所述第一修正處理模塊300，可以用于基于所述初始近地表速度場，對所述橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第一修正處理。

所述第二修正處理模塊400，可以用于基于所述初始近地表速度場，對所述射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度進行第二修正處理。

所述第三修正處理模塊500，可以用于基于所述初始近地表速度場，對所述垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度進行第三修正處理。

所述目標近地表速度場確定模塊600，可以用于根據(jù)所述第一修正處理后的橫向邊界區(qū)域內(nèi)的網(wǎng)格的速度、所述第二修正處理后的射線稀疏區(qū)域內(nèi)的網(wǎng)格的速度，以及所述第三修正處理后的垂直邊界區(qū)域內(nèi)的網(wǎng)格的速度，確定所述目的工區(qū)的目標近地表速度場。

圖8是本申請確定近地表速度場的裝置實施例中垂直邊界區(qū)域確定模塊的組成結(jié)構(gòu)圖。如圖8所示，圖7中垂直邊界區(qū)域確定模塊200，可以包括：第一起始位置確定模塊210、第二起始位置確定模塊220、目標起始位置確定模塊230和區(qū)域確定模塊240。

所述第一起始位置確定模塊210，可以用于根據(jù)所述預設(shè)網(wǎng)格模型的橫向邊界區(qū)域和射線稀疏區(qū)域，以及所述初至波射線數(shù)據(jù)，確定所述預設(shè)網(wǎng)格模型的第一垂直邊界起始位置。

所述第二起始位置確定模塊220，可以用于根據(jù)所述第一垂直邊界起始位置可以確定所述預設(shè)網(wǎng)格模型中第二垂直邊界起始位置。

所述目標起始位置確定模塊230，可以用于將所述第一垂直邊界起始位置和所述第二垂直邊界起始位置作為所述目標垂直邊界起始位置。

所述區(qū)域確定模塊240，可以用于將所述目標垂直邊界起始位置至所述預設(shè)網(wǎng)格模型的底表面位置所構(gòu)成的區(qū)域作為所述預設(shè)網(wǎng)格模型的垂直邊界區(qū)域。

圖9是本申請確定近地表速度場的裝置實施例中第三修正處理模塊的組成結(jié)構(gòu)圖。如圖9所示，圖7第三修正處理模塊500，可以包括：起始位置速度確定模塊510、起始位置目標速度確定模塊520、邊界位置目標速度確定模塊530和替換模塊540。

所述起始位置速度確定模塊510，可以用于根據(jù)所述初始近地表速度場和所述垂直邊界區(qū)域中目標垂直邊界起始位置，確定所述目標垂直邊界起始位置處的網(wǎng)格的速度。

所述起始位置目標速度確定模塊520，可以用于對所述目標垂直邊界初始位置處的網(wǎng)格的速度進行第一平滑處理，得到所述目標垂直邊界初始位置處的網(wǎng)格的目標速度。

所述邊界位置目標速度確定模塊530，可以用于基于所述目標垂直邊界初始位置處的網(wǎng)格的目標速度，確定所述垂直邊界區(qū)域中目標垂直邊界位置處的網(wǎng)格的目標速度。

所述替換模塊540，可以用于將所述近地表速度場中目標垂直邊界位置處的網(wǎng)格的速度替換為所述目標垂直邊界位置處的網(wǎng)格的目標速度。

所述確定近地表速度場的裝置實施例與所述確定近地表速度場的方法實施例相對應(yīng)，可以實現(xiàn)所述確定近地表速度場的方法實施例，并取得方法實施例的技術(shù)效果。

在20世紀90年代，對于一個技術(shù)的改進可以很明顯地區(qū)分是硬件上的改進(例如，對二極管、晶體管、開關(guān)等電路結(jié)構(gòu)的改進)還是軟件上的改進(對于方法流程的改進)。然而，隨著技術(shù)的發(fā)展，當今的很多方法流程的改進已經(jīng)可以視為硬件電路結(jié)構(gòu)的直接改進。設(shè)計人員幾乎都通過將改進的方法流程編程到硬件電路中來得到相應(yīng)的硬件電路結(jié)構(gòu)。因此，不能說一個方法流程的改進就不能用硬件實體模塊來實現(xiàn)。例如，可編程邏輯器件(programmablelogicdevice,pld)(例如現(xiàn)場可編程門陣列(fieldprogrammablegatearray，fpga))就是這樣一種集成電路，其邏輯功能由用戶對器件編程來確定。由設(shè)計人員自行編程來把一個數(shù)字系統(tǒng)“集成”在一片pld上，而不需要請芯片制造廠商來設(shè)計和制作專用的集成電路芯片2。而且，如今，取代手工地制作集成電路芯片，這種編程也多半改用“邏輯編譯器(logiccompiler)”軟件來實現(xiàn)，它與程序開發(fā)撰寫時所用的軟件編譯器相類似，而要編譯之前的原始代碼也得用特定的編程語言來撰寫，此稱之為硬件描述語言(hardwaredescriptionlanguage，hdl)，而hdl也并非僅有一種，而是有許多種，如abel(advancedbooleanexpressionlanguage)、ahdl(alterahardwaredescriptionlanguage)、confluence、cupl(cornelluniversityprogramminglanguage)、hdcal、jhdl(javahardwaredescriptionlanguage)、lava、lola、myhdl、palasm、rhdl(rubyhardwaredescriptionlanguage)等，目前最普遍使用的是vhdl(very-high-speedintegratedcircuithardwaredescriptionlanguage)與verilog2。本領(lǐng)域技術(shù)人員也應(yīng)該清楚，只需要將方法流程用上述幾種硬件描述語言稍作邏輯編程并編程到集成電路中，就可以很容易得到實現(xiàn)該邏輯方法流程的硬件電路。

控制器可以按任何適當?shù)姆绞綄崿F(xiàn)，例如，控制器可以采取例如微處理器或處理器以及存儲可由該(微)處理器執(zhí)行的計算機可讀程序代碼(例如軟件或固件)的計算機可讀介質(zhì)、邏輯門、開關(guān)、專用集成電路(applicationspecificintegratedcircuit，asic)、可編程邏輯控制器和嵌入微控制器的形式，控制器的例子包括但不限于以下微控制器：arc625d、atmelat91sam、microchippic18f26k20以及siliconelabsc8051f320，存儲器控制器還可以被實現(xiàn)為存儲器的控制邏輯的一部分。

本領(lǐng)域技術(shù)人員也知道，除了以純計算機可讀程序代碼方式實現(xiàn)控制器以外，完全可以通過將方法步驟進行邏輯編程來使得控制器以邏輯門、開關(guān)、專用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來實現(xiàn)相同功能。因此這種控制器可以被認為是一種硬件部件，而對其內(nèi)包括的用于實現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)?；蛘呱踔粒梢詫⒂糜趯崿F(xiàn)各種功能的裝置視為既可以是實現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

上述實施例闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計算機芯片或?qū)嶓w實現(xiàn)，或者由具有某種功能的產(chǎn)品來實現(xiàn)。

為了描述的方便，描述以上裝置時以功能分為各種單元分別描述。當然，在實施本申請時可以把各單元的功能在同一個或多個軟件和/或硬件中實現(xiàn)。

通過以上的實施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請可借助軟件加必需的通用硬件平臺的方式來實現(xiàn)?；谶@樣的理解，本申請的技術(shù)方案本質(zhì)上或者說對現(xiàn)有技術(shù)做出貢獻的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來，在一個典型的配置中，計算設(shè)備包括一個或多個處理器(cpu)、輸入/輸出接口、網(wǎng)絡(luò)接口和內(nèi)存。該計算機軟件產(chǎn)品可以包括若干指令用以使得一臺計算機設(shè)備(可以是個人計算機，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請各個實施例或者實施例的某些部分所述的方法。該計算機軟件產(chǎn)品可以存儲在內(nèi)存中，內(nèi)存可能包括計算機可讀介質(zhì)中的非永久性存儲器，隨機存取存儲器(ram)和/或非易失性內(nèi)存等形式，如只讀存儲器(rom)或閃存(flashram)。內(nèi)存是計算機可讀介質(zhì)的示例。計算機可讀介質(zhì)包括永久性和非永久性、可移動和非可移動媒體可以由任何方法或技術(shù)來實現(xiàn)信息存儲。信息可以是計算機可讀指令、數(shù)據(jù)結(jié)構(gòu)、程序的模塊或其他數(shù)據(jù)。計算機的存儲介質(zhì)的例子包括，但不限于相變內(nèi)存(pram)、靜態(tài)隨機存取存儲器(sram)、動態(tài)隨機存取存儲器(dram)、其他類型的隨機存取存儲器(ram)、只讀存儲器(rom)、電可擦除可編程只讀存儲器(eeprom)、快閃記憶體或其他內(nèi)存技術(shù)、只讀光盤只讀存儲器(cd-rom)、數(shù)字多功能光盤(dvd)或其他光學存儲、磁盒式磁帶，磁帶磁磁盤存儲或其他磁性存儲設(shè)備或任何其他非傳輸介質(zhì)，可用于存儲可以被計算設(shè)備訪問的信息。按照本文中的界定，計算機可讀介質(zhì)不包括短暫電腦可讀媒體(transitorymedia)，如調(diào)制的數(shù)據(jù)信號和載波。

本說明書中的各個實施例均采用遞進的方式描述，各個實施例之間相同相似的部分互相參見即可，每個實施例重點說明的都是與其他實施例的不同之處。尤其，對于系統(tǒng)實施例而言，由于其基本相似于方法實施例，所以描述的比較簡單，相關(guān)之處參見方法實施例的部分說明即可。

本申請可用于眾多通用或?qū)Ｓ玫挠嬎銠C系統(tǒng)環(huán)境或配置中。例如：個人計算機、服務(wù)器計算機、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的消費電子設(shè)備、網(wǎng)絡(luò)pc、小型計算機、大型計算機、包括以上任何系統(tǒng)或設(shè)備的分布式計算環(huán)境等等。

本申請可以在由計算機執(zhí)行的計算機可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類型的例程、程序、對象、組件、數(shù)據(jù)結(jié)構(gòu)等等。也可以在分布式計算環(huán)境中實踐本申請，在這些分布式計算環(huán)境中，由通過通信網(wǎng)絡(luò)而被連接的遠程處理設(shè)備來執(zhí)行任務(wù)。在分布式計算環(huán)境中，程序模塊可以位于包括存儲設(shè)備在內(nèi)的本地和遠程計算機存儲介質(zhì)中。

雖然通過實施例描繪了本申請，本領(lǐng)域普通技術(shù)人員知道，本申請有許多變形和變化而不脫離本申請的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請的精神。