一種確定近地表速度模型的方法及裝置的制造方法
【技術領域】
[0001] 本申請設及地球物理勘探技術領域，特別設及一種確定近地表速度模型的方法及 裝置。
【背景技術】
[0002] 在地球物理勘探技術領域，近地表速度模型不僅影響地震勘探中地震波的激發(fā)和 接收，還會對獲取的原始資料的信噪比產(chǎn)生影響。因此，弄清近地表速度模型，對于地震采 集中的激發(fā)問題、資料處理中的靜校正問題W及波場延拓問題等均具有重要意義。
[0003] 目前，可W利用微測井數(shù)據(jù)或者小折射數(shù)據(jù)對近地表進行調查W構建近地表速度 模型。其中，微測井數(shù)據(jù)在縱向的采樣間隔一般都較小，初至的信噪比較高，雖然存在著固 定偏移距的井口距，但經(jīng)過垂直校正后，可W認為地震波在地層中垂直傳播，所W在微測井 位置處，其解釋精度還是相當高的；小折射數(shù)據(jù)的常規(guī)解釋方法雖要求地表比較平坦，并且 要求地層速度在橫向方向變化較小，但其道距較小，初至信噪比高，其解釋精度在小折射位 置處往往較高。但是運兩種方法僅僅可W得到調查點處有效深度范圍內縱向的速度變化， 是一個點數(shù)據(jù)，至于近地表速度的空間變化如何，該方法無能無力。實際工作中也不可能布 設高密度的表層調查控制點進行近地表結構的調查。
[0004] 另外，利用初至走時層析反演結果也可W構建近地表速度模型。利用初至走時層 析反演的方法構建的近地表速度模型可W很好地反映近地表速度的空間變化趨勢，然而受 到炮檢距、近道初至精度W及地震數(shù)據(jù)采集密度等因素的限制，構建的近地表速度模型的 精度往往不高，特別在淺地表處表現(xiàn)的尤為明顯。
[0005] 應該注意，上面對技術背景的介紹只是為了方便對本申請的技術方案進行清楚、 完整的說明，并方便本領域技術人員的理解而闡述的。不能僅僅因為運些方案在本申請的
【背景技術】部分進行了闡述而認為上述技術方案為本領域技術人員所公知。

【發(fā)明內容】

[0006] 本申請實施例的目的在于提供一種確定近地表速度模型的方法及裝置，W提高近 地表速度模型的精度。
[0007] 本申請實施例提供的一種確定近地表速度模型的方法及裝置是運樣實現(xiàn)的：
[0008] 一種確定近地表速度模型的方法，包括：
[0009] 對待測區(qū)域的表層調查數(shù)據(jù)進行處理，得到所述待測區(qū)域的近地表第一速度模 型；
[0010] 對所述待測區(qū)域的地震數(shù)據(jù)進行走時層析反演，得到所述待測區(qū)域的近地表第二 速度模型；
[0011] 計算所述第一速度模型與所述第二速度模型中相同位置處的數(shù)據(jù)的差值，確定當 計算的差值大于預設闊值時對應的位置并從所述第一速度模型中剔除確定的位置處的數(shù) 據(jù)；
[0012] 對所述第一速度模型中剔除后的數(shù)據(jù)進行插值計算，得到近地表的初始速度模 型；
[0013] 對所述第二速度模型與剔除數(shù)據(jù)后的第一速度模型中相同位置處的數(shù)據(jù)進行計 算，得到各個位置處對應的權重系數(shù)；
[0014] 基于計算出的權重系數(shù)，在所述待測區(qū)域的縱向方向W及橫向方向分別進行插值 計算，得到所述待測區(qū)域中不同位置處的權重因子；
[0015] 利用所述權重因子，對所述初始速度模型進行迭代校正，直至滿足預設條件為 止；
[0016] 將滿足所述預設條件的迭代校正后的速度模型確定為近地表速度模型。
[0017] -種確定近地表速度模型的裝置，包括：
[0018] 第一速度模型獲取單元，用來對待測區(qū)域的表層調查數(shù)據(jù)進行處理，得到所述待 測區(qū)域的近地表第一速度模型；
[0019] 第二速度模型獲取單元，用來對所述待測區(qū)域的地震數(shù)據(jù)進行走時層析反演，得 到所述待測區(qū)域的近地表第二速度模型；
[0020] 數(shù)據(jù)剔除單元，用來計算所述第一速度模型與所述第二速度模型中相同位置處的 數(shù)據(jù)的差值，確定當計算的差值大于預設闊值時對應的位置并從所述第一速度模型中剔除 確定的位置處的數(shù)據(jù)；
[0021] 權重系數(shù)計算單元，用來對所述第二速度模型與剔除數(shù)據(jù)后的第一速度模型中相 同位置處的數(shù)據(jù)進行計算，得到各個位置處對應的權重系數(shù)；
[0022] 初始速度模型獲取單元，用來對所述第一速度模型中剔除后的數(shù)據(jù)進行插值計 算，得到近地表的初始速度模型；
[0023] 權重因子獲取單元，用來基于計算出的權重系數(shù)，在所述待測區(qū)域的縱向方向W 及橫向方向分別進行插值計算，得到所述待測區(qū)域中不同位置處的權重因子；
[0024] 迭代校正單元，用來利用所述權重因子，對所述初始速度模型進行迭代校正，直至 滿足預設條件為止；
[00巧]近地表速度模型確定單元，用來將滿足所述預設條件的迭代校正后的速度模型確 定為近地表速度模型。
[00%] 本申請實施例提供的一種確定近地表速度模型的方法及裝置，分別對待測區(qū)域的 表層調查數(shù)據(jù)W及地震數(shù)據(jù)進行處理，得到對應的近地表第一速度模型W及第二速度模 型，然后可W利用所述第二速度模型對所述第一速度模型進行空間約束，去除所述第一速 度模型中不準確的數(shù)據(jù)并基于去除數(shù)據(jù)后的第一速度模型構建初始速度模型。然后可W利 用計算出的權重因子對所述初始速度模型進行迭代校正，從而得到精度較高的近地表速度 模型。
[0027] 參照后文的說明和附圖，詳細公開了本申請的特定實施方式，指明了本申請的原 理可W被采用的方式。應該理解，本申請的實施方式在范圍上并不因而受到限制。在所附 權利要求的精神和條款的范圍內，本申請的實施方式包括許多改變、修改和等同。
[0028] 針對一種實施方式描述和/或示出的特征可相同或類似的方式在一個或更 多個其它實施方式中使用，與其它實施方式中的特征相組合，或替代其它實施方式中的特 征。
[0029] 應該強調，術語"包括/包含"在本文使用時指特征、整件、步驟或組件的存在，但 并不排除一個或更多個其它特征、整件、步驟或組件的存在或附加。
【附圖說明】
[0030] 所包括的附圖用來提供對本申請實施例的進一步的理解，其構成了說明書的一部 分，用于例示本申請的實施方式，并與文字描述一起來闡釋本申請的原理。顯而易見地，下 面描述中的附圖僅僅是本申請的一些實施例，對于本領域普通技術人員來講，在不付出創(chuàng) 造性勞動性的前提下，還可W根據(jù)運些附圖獲得其他的附圖。在附圖中：
[0031] 圖1為本申請實施例提供的一種確定近地表速度模型的方法流程圖；
[0032] 圖2為本申請實施例對不同地層處的垂直時間進行擬合的示意圖；
[0033] 圖3為本申請實施例提供的一種確定近地表速度模型的裝置的功能模塊圖。
【具體實施方式】
[0034] 為了使本技術領域的人員更好地理解本申請中的技術方案，下面將結合本申請實 施例中的附圖，對本申請實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施 例僅僅是本申請一部分實施例，而不是全部的實施例?；诒旧暾堉械膶嵤├?，本領域普通 技術人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其它實施例，都應當屬于本申請保護 的范圍。
[0035] 圖1為本申請實施例提供的一種確定近地表速度模型的方法流程圖。雖然下文描 述流程包括W特定順序出現(xiàn)的多個操作，但是應該清楚了解，運些過程可W包括更多或更 少的操作，運些操作可W順序執(zhí)行或并行執(zhí)行（例如使用并行處理器或多線程環(huán)境）。如圖 1所示，所述方法可W包括：
[0036] S1 :對待測區(qū)域的表層調查數(shù)據(jù)進行處理，得到所述待測區(qū)域的近地表第一速度 模型。
[0037] 在本申請實施例中，所述待測區(qū)域的表層調查數(shù)據(jù)可W為微測井數(shù)據(jù)或者小折射 數(shù)據(jù)，利用所述待測區(qū)域的表層調查數(shù)據(jù)可W得到所述待測區(qū)域的近地表第一速度模型。 該第一速度模型在所述待測區(qū)域的縱向單點處的精度較高，但在橫向方向上卻容易受到表 層控制點密度的影響而導致精度不足。具體地可W表現(xiàn)為在橫向方向上的速度值可能出現(xiàn) 突變，形成奇異值，因此需要對所述第一速度模型進行后續(xù)的處理。
[0038] 在本申請實施例中，當所述表層調查數(shù)據(jù)為微測井數(shù)據(jù)時，可W從所述微測井數(shù) 據(jù)中拾取不同地層處的初至時間，然后可W根據(jù)下述公式將所述不同地層處的初至時間轉 化為垂直時間：
[0039]
[0040] 其中，t。為某一地層處的垂直時間，t為該地層處的初至時間，Η為該地層處激發(fā) 點的深度，D為該地層處接收點到井口的距離。
[0041] 得到不同地層處的垂直時間后，可W利用最小二乘法對所述垂直時間進行擬合， 從而可W得到不同地層處的速度值W及與所述速度值相對應的厚度值。具體地，圖2為本 申請實施例對不同地層處的垂直時間進行擬合的示意圖。如圖2所示，可W將不同地層處 的垂直時間繪制在同一坐標系內，該坐標系的橫坐標可w為垂直時間值，縱坐標可w為與 所述垂直時間值相對應的地層深度值。由于位于同一速度層內的垂直時間的點可W連成一 直線，從而可W根據(jù)垂直時間點的分布規(guī)律，劃分出各個速度層的對應的垂直時間點，然后 可W對每組垂直時間點利用最小二乘法擬合成直線，該直線的斜率的倒數(shù)便可W為該速度 層對應的速度值。而該直線上起始垂直時間點W及終止垂直時間點對應的深度差即可W為 該速度層的速度值對應的厚度值。從而可W將不同地層處確定出的速度值W及與所述速度 值相對應的厚度值確定為近地表第一速度模型。
[0042] 另外，在本申請實施例中，當所述表層調查數(shù)據(jù)為小折射數(shù)據(jù)時，同樣可W先從所 述表層調查數(shù)據(jù)中拾取不同地層處的初至時間，然后利用同樣的方法利用最小二乘法對所 述不同地層處的初至時間進行擬合，得到擬合直線，從而可W得到不同地層處對應的速度 值。接著，本申請實施例可W按照下述公式計算不同地層處的速度值對應的厚度值：
[0043]
[0044] 其中，m為大于或者等于2的整數(shù)，hm1為第m