初至波層析影響因素分析

1、初至波層析影響因素分析劉玉柱 3董良國( 同濟(jì)大學(xué)海洋地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室)劉玉柱 ,董良國 . 初至波層析影響因素分析 . 石油地球物理勘探 ,2007 ,42 (5) :544553摘要 在山地地震勘探中 ,目前廣泛采用初至波走時(shí)層析方法反演起伏地表速度結(jié)構(gòu) ,但其最終反演結(jié)果受許多因素影響 。本文通過大量數(shù)值試驗(yàn)對影響表層模型層析反演精度的方程解法 、先驗(yàn)信息利用 、觀測系統(tǒng) 、初 始模型及平滑因子等主要因素進(jìn)行了詳細(xì)的分析對比 ,總結(jié)出了最佳參數(shù)的選擇原則和范圍 ,并建立了內(nèi)部約 束方程 。將本文方法應(yīng)用于實(shí)際資料處理 ,取得了較好的成像效果 。關(guān)鍵詞 層析反演 表層結(jié)構(gòu) 初至波 先驗(yàn)信

2、息 正則化影響因素1引言2參數(shù)分析在山地地震勘探中 ,近地表速度結(jié)構(gòu)反演是一個(gè)亟待卻又很難解決的問題 。盡管國內(nèi)外學(xué)者對該問題 進(jìn)行了大量的研究 18 ,但至今仍沒有很好地解決 , 從而嚴(yán)重影響表層校正效果及成像質(zhì)量 。通過細(xì)致調(diào)研 ,我們認(rèn)為應(yīng)充分利用地震波的走時(shí)和波形 (包 括振幅 、相位 、頻率等) 信息 ,結(jié)合線性反演 、擬線性反 演 、完全非線性反演等多種反演手段 ,在地表調(diào)查資 料和已知地質(zhì)背景知識等先驗(yàn)信息的約束下 ,采用多 信息 、多方法聯(lián)合反演的方法可望有效地解決此問題 。盡管聯(lián)合反演是最終目標(biāo) ,但它仍依賴于單一反 演技術(shù)的提高 。目前廣泛采用的單一反演方法包括折射層析 9

3、10 和初至波走時(shí)層析 2 ,47 ,11 ,12 兩種 。折射層析以層狀模型假設(shè)為基礎(chǔ) ,盡管折射層析有時(shí)可 以得到比較好的靜校正效果 ,但層狀假設(shè)的前提條件并不總能得到滿足 。初至波走時(shí)層析沒有層狀模型 假設(shè) ,它適用于任意表層速度結(jié)構(gòu) ,反演結(jié)果能較好 地反映表層速度的低頻趨勢 。鑒于初至層析仍存在 反演精度低 、穩(wěn)定性差 、反演結(jié)果受諸多因素影響等 問題 ,本文通過大量數(shù)值試驗(yàn) ,對影響反演精度的主要因素進(jìn)行詳細(xì)對比分析 ,旨在探討各參數(shù)影響層析 反演效果的規(guī)律 ,總結(jié)出各參數(shù)優(yōu)選的原則和范圍 , 達(dá)到提高初至層析反演精度及反演穩(wěn)定性的目的 。在線性假設(shè)下 ,初至波走時(shí)層析成像簡化為求

4、解大規(guī)模線性方程組ls = t( 1)其中 : l 為 m n 維矩陣 , 矩陣元素 l kj 代表第 k 條射線在第 j 個(gè)模型參數(shù)單元內(nèi)的長度 ;s 是長度為 n 的列向量 , 代表模型慢度參數(shù)修正量 ;t 為 m 維列 向量 , 代表觀測的初至波走時(shí)與理論計(jì)算的初至波 走時(shí)殘差 13 。不考慮拾取誤差與射線追蹤誤差 ,概括起來影響層析反演結(jié)果的因素包括 : 方程解法 ; 先驗(yàn)信息的利用 ; 觀測系統(tǒng) ; 初始模型 ; 平滑因子 。下面通過理論分析 、結(jié)合數(shù)值實(shí)驗(yàn)探討上述各因素對層析反演結(jié)果的影響 。2 . 1 方程解法求解方程組式 (1) 通常有反投影法 (b p t) 、代數(shù) 重建法

5、( a r t) 、同時(shí)迭代重建法 ( sir t) 、奇異值分 解法 ( svd) 和投影類迭代方法等 14 多種算法 。其 中反投影法是最傳統(tǒng)的一種方法 ,其基本思想是將 旅行時(shí)殘差按射線在模型參數(shù)單元內(nèi)的覆蓋權(quán)數(shù)進(jìn)行反投影 。這種方法雖然分辨率低于代數(shù)重建法和 投影法 ,但具有較高的計(jì)算效率 ,而且穩(wěn)定性好 。計(jì) 算公式為3 上海市同濟(jì)大學(xué)海洋與地球科學(xué)學(xué)院 ,200092本文于 2006 年 11 月 23 日收到 ,修改稿于 2007 年 6 月 8 日收到 。本課題獲得教育部“新世紀(jì)優(yōu)秀人才支持計(jì)劃”( n cep20520384) 與“國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目”( 2006

6、cb202402) 資助 。巨大的稀疏矩陣 ,而且矩陣的條件數(shù)有時(shí)很大 (病態(tài)問題比較嚴(yán)重) ,則求解方程組時(shí)收斂速度慢且不穩(wěn) 定 ;l sc g 法改變了對該矢量的計(jì)算方式 ,使其能夠 處理矩陣為病態(tài) 的情 況 ; l sq r 法 是利 用 l a nczo s 法求解最小二乘問題的一種投影法 ,它可壓制數(shù)據(jù) 誤差的傳遞 ,收斂速度較快 ,適合求解大型病態(tài)線性 代數(shù)方程組 14 。圖 1a 與圖 1 b 是基于實(shí)際數(shù)據(jù)的拾取初至采用 相同的初始模型與反演參數(shù) ,分別利用反投影法與 l sq r 算法得到的層析成像結(jié)果 。圖 2 是圖 1a 與 圖 1 b 地表以下 40 m 處的速度切片對

7、比圖 。不考慮 精度 ,可以看出 l sq r 算法確實(shí)比反投影法具有更高的分辨率 , 但 l sq r 法的計(jì)算效率僅為反投影方 法的 1/ 51/ 10 。tkl kjkksj =( 2)l kj +1k其中 : l kj 代表第 j 個(gè)網(wǎng)格單元內(nèi)的所有射線總長k度 ; tk 表示第 k 條射線的初至波走時(shí)殘差 ; gk =l kj +2 , 代表第 k 條射線的長度 ; 1 、2 為防止計(jì)j算溢出的阻尼因子 13 。投影類迭代方法的基本思路是 :尋找方程組的近似解 , 使得該近似解屬于由 n 維空間中 m ( m n) 個(gè)無關(guān)向量構(gòu)成的子空間 ,而殘 差量與子空間垂直 。投影類迭代方法還

8、可細(xì)分為 l a nczo s 法 、共軛梯度法 ( c g) 、最小二乘共軛梯度法 (l sc g) 及最小二乘 q r 分解法 ( l sq r) 。共軛梯度 法具有快速高效的特點(diǎn) ,其主要問題是要求存儲g2 . 2先驗(yàn)信息的利用地球物理反演是一個(gè)典型的混定問題 ,在走時(shí) 層析中表現(xiàn)為有些模型參數(shù)的射線覆蓋率高 ,有些 模型參數(shù)的射線覆蓋率低 ,甚至為零 。取該混定問 題的目標(biāo)函數(shù)為(s) = ( ls - t) t ( ls - t) +2s ts( 3)混定特性導(dǎo)致了模型參數(shù)的修改量與射線的覆蓋密度密切相關(guān) ,但根據(jù)第三類先驗(yàn)信息 ,反演結(jié)果應(yīng)該 比較平坦或平滑 ,因此必須在每次迭代后

9、對模型參 數(shù)進(jìn)行平滑 。為了更有效地利用先驗(yàn)信息 ,應(yīng)該采用正則化 方法將先驗(yàn)信息的約束納入層析方程組中 ,這里稱 之為內(nèi)部約束方法 ?；诨於▎栴}阻尼最小二乘解(式 (4) ) ,筆者推導(dǎo)出包含這三種先驗(yàn)信息的內(nèi)部約 束方法對應(yīng)的阻尼最小二乘解公式 ,即其中為阻尼因子或正則化因子 。令 5 = 0, 則有5s= ( l t l +2 i) - 1 l ttse st( 4)t2 i +2 t2 t2 te st( l l +1 w1 w1 +2 w2 w2 +3 w3 w3 )s式 ( 4) 即為式 ( 1) 的阻尼最小二乘解 517 。上述反演中 , 先驗(yàn)信息的利用非常重要 。本文 將初至

10、波走時(shí)層析中的先驗(yàn)信息分為三類 : 根據(jù) 近地表調(diào)查資料得到的某些參數(shù)的精確值 ; 從已 知背景地質(zhì)知識得到的某些參數(shù)的取值范圍 ; 模型空間的參數(shù)分布特點(diǎn) 。針對不同類型的先驗(yàn)信息 構(gòu)造出相應(yīng)的正則化公式 , 這三種正則化公式可用 一個(gè)統(tǒng)一式子表達(dá) 。對于一般的走時(shí)層析 , 先驗(yàn)信息的應(yīng)用主要是 約束更新后的速度模型 , 這里稱之為外部約束方法 。 根據(jù)先驗(yàn)信息的不同 , 外部約束方法分為三種 : 一種是緊約束 , 即根據(jù)第一類先驗(yàn)信息強(qiáng)行將某些參數(shù) 值修改為先驗(yàn)值 ; 第二種為寬約束 , 即根據(jù)第二類先 驗(yàn)信息約束模型的每個(gè)參數(shù)到預(yù)先設(shè)定的范圍之t 2 tl t + 1 w1 w1 ( s

11、- s0 ) +=2 t2 t+2 w2 k +3 w3 w3 ( - s0 )( 5)其中 :w1 ,w2 ,w3 依次為第一 、二 、三類正則化矩陣 ;為衡量數(shù) 據(jù)誤 差權(quán) 重與 模 型 誤 差 權(quán) 重 的 阻 尼 因 子 ; 1 ,2 ,3 分別為第一 、二 、三類正則化因子 , ,1 ,2 ,3 可以是由嘗試得到的較小量 ; s0 為當(dāng)前模 型的慢度向量 ; s為期望模型的慢度向量 。為了驗(yàn)證不同約束方法的層析反演效果 ,本文 設(shè)計(jì)了圖 3 所示的二維起伏地表理論模型 。模型中第一層 為強(qiáng) 烈 橫向 變速 ( 1000 3600 m/ s) 的 低 速 帶 ,第二 、三層為均勻介質(zhì) ,

12、速度分別為 3500 m/ s 和1900 m/ s 。本文特意設(shè)計(jì)了一個(gè)低速層 , 目的是檢 驗(yàn)層析方法對低速層的反演能力 ,后文將看到 ,由于 炮檢距小 ,初至波無法穿越低速層而不能將其反演出來 ,但這并不影響本文后續(xù)分析 。內(nèi) ( 即 若vmax , 則 vi=vmax ; 若 vivmin , 則vivi = vmin ) ;第三種為平滑處理 ,走時(shí)層析反演本身的圖 3 二維起伏地表理論模型該理論模型采用 608 炮觀測數(shù)據(jù) ,均勻分布在橫向 9950 39310 m 的 范 圍 內(nèi) , 觀 測 系 統(tǒng) 為 20002202022022000 m 。內(nèi)部約束方法對應(yīng)的阻尼最小二乘解為式

13、 ( 5) ,1 ,3 都等于 1 , w3 取一 階 差分 矩 陣 。外部約束方法對應(yīng)的阻尼最小二乘解為式 (4) 。兩種方法都采用 l sq r 算法求解線性代數(shù)方程組 。外部約束和內(nèi)部約束對應(yīng)的反演結(jié)果分別如圖 4a和圖 4 b 所示 ,目標(biāo)函數(shù)值 ( 拾取初至?xí)r間與最終層 析結(jié)果的理論計(jì)算初至?xí)r間之間殘差的均方根 ,單 位為 ms) 分別為 8 . 0928 ms 和 4 . 781 ms 。不難看出 , 內(nèi)部約束比外部約束具有更高的分辨率和反演精 度 ,特別是圈出的位置更加明 顯 。圖 4c 是 理論 模型 、外部約束方法 、內(nèi)部約束方法層析結(jié)果 ( 地表以下 40 m 速度切片)

14、對比圖 , 圖 4 d 是拾取初至?xí)r 間 (綠線) 與基于外部約束層析結(jié)果 ( 圖 4a ) 計(jì)算的初 至?xí)r間 (紅線) 對比圖 ,圖 4e 是拾取初至?xí)r間 ( 綠線)與基于內(nèi)部約束層析結(jié)果 ( 圖 4 b) 計(jì)算的初至?xí)r間(紅線) 對比圖 。從這些圖也可看出 ,內(nèi)部約束方法 的層析結(jié)果更接近理論模型 ,而且內(nèi)部約束層析成 像方法可較準(zhǔn)確地反演出更加豐富的高波數(shù)成分 。2 . 3觀測系統(tǒng)觀測系統(tǒng)參數(shù)包括炮間距 、道間距 、最大炮檢距 等 ,它關(guān)系到地面接收到的地震數(shù)據(jù)中信息的豐富程度 ,直接影響表層反演的效果 。理論上 ,炮點(diǎn) 、檢波點(diǎn)越密反演結(jié)果越好 ,但考慮到成本問題 ,道間距 不 宜 太

15、小 。圖 5 a 顯示 的是 2 0 m 道 距的 層 析結(jié) 果 ,圖 5 b 顯示的是 1 m 道距的層析結(jié)果 。與圖 3 的理論模型對比可以發(fā)現(xiàn) ,除了個(gè)別地 方存在差異外 ,總體上基本一致 ,這一點(diǎn)從圖 5c 兩 者的收斂曲線也可看出 , 因?yàn)閮烧咦罱K收斂到幾乎相同的目標(biāo)函數(shù)值 。圖 5d 是理論模型 、20 m 道距 、1 m 道距層析結(jié)果在地表以下 60 m 處的速度切片對比圖 ,從該圖也可以看出 ,20 m 道距和 1 m 道距的層 析結(jié)果基本相當(dāng) (這可能是由平滑因子造成的 ,筆者 盡量將平滑因子縮小再對比 20 m 道距和 1 m 道距的 層析結(jié)果 ,得到相同結(jié)論) , 所以

16、, 為了節(jié)約計(jì)算成 本 ,可選取相對較大的道間距 。圖 5c 20m (紅線) 和 1 m (綠線) 道距的層析收斂曲線圖 5d 理論模型 (紅線) 、20 m (藍(lán)線) 和 1 m (綠線) 道距的層析反演結(jié)果對比速度切片 ( 地表以下 60 m)另一方面 ,野外觀測系統(tǒng)提供了足夠大的炮檢距 ,實(shí)際應(yīng)用中可只用小炮檢距的初至波進(jìn)行反演 。 盡管炮檢距越大反演的深度越大 ,但由于層析的平 均效應(yīng) ,反演的表層速度的分辨率會降低 ,反演的精度也會降低 。所以如果只為進(jìn)行表層校正而只關(guān)心較淺的表層速度結(jié)構(gòu) ,則只需小炮檢距初至信息即 可 。圖 6a 和圖 6 b 是最大炮檢距分別為 1000 m

17、和4000 m 的層析反演結(jié)果 。兩者差異較大 ,圖 6a 反演深度小但精度高 , 圖 6 b 反演深度大但精度低 。從圖 6c的層析收斂曲線也可得出同樣結(jié)論 。2 . 4 初始模型初始模型的質(zhì)量直接影響最終反演結(jié)果的精 度 ?，F(xiàn)今實(shí)際資料處理的初始模型一般是通過近地 表調(diào)查資料結(jié)合先驗(yàn)地質(zhì)知識建立起來的 。在缺乏 先驗(yàn)知識的情況下 ,可采用常速模型或等梯度模型 。下面以這兩種初始模型的實(shí)際資料處理來說明不同 初始模型對最終層析結(jié)果的影響 。圖 7a 與圖 7c 分 別是采用常速初始模型與等梯度初始模型 (圖 7b) 得 到 的層析結(jié)果 , 可見淺部反演結(jié)果相當(dāng) , 深部較大差圖 6c 最大炮

18、檢距為 1000 m (紅線) 和 4000m (綠線)的層析收斂曲線圖 7a 實(shí)際資料常速初始模型對應(yīng)的層析結(jié)果別是由初始模型不同引起的 ,層析并沒有修改中深層的速度信息 ,圖 7d 表明兩者在地表以下 100m 處仍具 有很好的相似性 。從圖 7e 可見梯度初始模型的反演 結(jié)果好于常速初始模型 ,而且層析收斂速度快 。2 . 5 平滑因子盡管正則化方法中的一階差分算子或二階差 分算子已 經(jīng) 對 模 型 參 數(shù) 進(jìn) 行 了 平 均 或 光 滑 處 理 ,但仍有必要對更新過的模型再進(jìn)行平滑處理 。平滑過程中的平滑因子對反演的質(zhì)量影響很大 。平 滑窗過大 會 降 低 反 演 分 辨 率 , 平

19、滑 窗 過 小 雖 可 提高反演分 辨 率 , 但 模 型 變 化 劇 烈 可 能 導(dǎo) 致 射 線 追圖 7e 常速初始模型 (綠) 與梯度初始模型 (紅)對應(yīng)的層析收斂曲線蹤 失 敗 , 進(jìn) 而 使 下 一 輪 的 層 析 方 程 求 解 不 穩(wěn) 定 。因此需合 理 選 取 平 滑 因 子 , 以 兼 顧 反 演 結(jié) 果 的 分辨率和穩(wěn)定性 。下面的數(shù)值試驗(yàn)分別采用 100 m 和 50 m 的平滑 窗進(jìn)行處理 。從圖 8 可知 ,較小平滑窗得到的層析結(jié)果并不理想 ,平滑窗過小甚至?xí)?dǎo)致反演過程不穩(wěn)定 。如圖 8c 所示 ,50 m 平滑窗層析很快收斂 ,但 目標(biāo)函 數(shù) 值 仍 然 很 大

20、。多 次 試 驗(yàn) 表 明 , 平 滑 窗 取100200 m 較適宜 。炮 ,炮間距與檢波器間距均為 20 m 。根據(jù)前面理論與數(shù)值試驗(yàn)分析結(jié)果 ,本次處理 選用如下參數(shù)和設(shè)置 : 采用 l sq r 算法求解反演 方程 ; 采用正則化處理方法修正反演方程 ,最終阻尼最小二乘解對應(yīng)式 (5) ,其中第三類正則化矩陣采 用一 階 差 分 算 子 ; 拾 取 初 至 的 最 大 炮 檢 距 為1000 m ,道距為 20 m ,炮距為 20 m ; 初始模型的淺 層 50 m 左右范圍是根據(jù)地表調(diào)查資料建立起來的 (圖 9a) ,中深層采用常速 ( 為了驗(yàn)證層析處理實(shí)際問題的能力 ,也為了簡便)

21、; 平滑因子取 100 m 。圖 9 b 是該 測 線 近 地 表 速 度 的 最 終 層 析 反 演 結(jié)果 , 圖 9c 顯示了該測線拾取初至 ( 綠線) 與據(jù)最 終反演模型計(jì)算的初至 ( 紅線) 的對比 。從該圖可 見 , 拾取初 至 與 理 論 計(jì) 算 初 至 吻 合 較 好 , 說 明 層析反演結(jié) 果 基 本 反 映 了 真 實(shí) 的 近 地 表 速 度 結(jié) 構(gòu) 。 圖 10a 、圖 10 b 分別是在折射層析靜校正與初至層 析靜校正基礎(chǔ)上得到的水平疊加剖面 。兩張剖面 的主要構(gòu)造基 本相 同 , 但 圖 10a 比 圖 10 b 的 分 辨 率更高一些 。另外 , 兩者在細(xì)節(jié)上有些差別

22、 , 圖中圖 8c 平滑因子分別為 50m (紅) 和 100m (綠)對應(yīng)的層析收斂曲線3實(shí)際資料處理為了驗(yàn)證以上分析結(jié)果 ,本文選取了中國西部地區(qū)一條二維測線的地震資料進(jìn)行了實(shí)際地震資料 處理 。該 測 線 所 在 地 區(qū) 地 表 起 伏 劇 烈 , 海 拔 為11003500 m , 地表結(jié)構(gòu)復(fù)雜 , 橫向速度變化嚴(yán)重 。測線長約為 35000 m ,排列長度為 15000 m ,共計(jì) 556圖 10b 實(shí)際數(shù)據(jù)初至層析靜校正疊加剖面圈出部分圖 10 b 比圖 10a 反映得更加細(xì)致 。4結(jié)論與討論初至波走時(shí)層析是目前廣泛應(yīng)用的一種近地表速度結(jié)構(gòu)反演方法 ,但影響層析反演效果的因素很 多

23、 。本文通過大量的數(shù)值試驗(yàn)分析了走時(shí)層析中的五大主要影響因素 ,即方程解法 、先驗(yàn)信息的利用 、觀測系統(tǒng) 、初始模型及平滑因子 。文中從理論分析與數(shù)值試驗(yàn)結(jié)果得出以下結(jié)論 :(1) 投影類比反投影法和迭代類反演法更適合 求解初至波層析中的大型方程組 ,分辨率更高 ,但效mindo ng2j ooa ndshincha ngsoo . ref ractio nto mo grap hy using a wavefo r m2inver sio n back2p ropagatio n tech2nique . geo p h y s ics , 2006 , 71 : 2130ivano vj

24、.j oint a nalysi sof ref ractio ns wit h surf acewav2es : a n i nver se sol utio n to t he ref ractio n2t raveltime p ro blem. geo p h ys ics , 2006 , 71 : 131138zho u h ua2wei . multi scale t raveltime to mo grap hy .geo p h y s ics , 2003 , 68 : 16391649zho u h ua2wei . fir st 2 break vertical sei

25、 smic p rofili ng to mo grap hy fo r vinto n salt do me . geo p h y s ics ,2006 , 71 : 2936李錄明等. 初至波表層模型層析反演. 石油地球物理勘探 , 2000 , 35 (5) : 559564李錄明 , 羅省賢. 復(fù)雜三維表層模型層析反演與靜校正. 石油地球物理勘探 , 2003 , 38 (6) : 636641井西利等. 建立速度模型的層析成像方法研究. 石油物探 , 2002 , 41 (1) : 7275shtivel man v . kinematic i nver sio n of fir

26、 st a r rival s of ref racted wave s a co mbined app roach . geo p h y s ics ,1996 , 61 ( 2) : 509519楊海申等. 復(fù)雜區(qū)三維折射靜校正技術(shù)與應(yīng)用效果.石油地球物理勘探 , 2005 , 40 (2) : 219225sheng j ia n2mi ng. ea rly ar rival wavefo r m to mo grap hy o n nea r2surf ace ref ractio n data . geo p h y s ics , 2006 ,71 ( 4) : 4757zha

27、ng j ie a nd to k soz m n . no nlinea r ref ractio n t raveltime to mo grap hy. geo p h y s ics , 1998 , 63 ( 5 ) :17261737sto r k c. l inear a sp ect s of to mo grap hic velocit y a nal2ysi s. geo p h ys ics , 1991 , 56 ( 4) : 483495吳律著. 層析基礎(chǔ)及其在井間地震中的應(yīng)用. 北京 :石油工業(yè)出版社 , 1997ta rantola a and valet te

28、b . generalized no n2li nea r in2 ver se p ro blems solved using t he lea st2squa re s criteri2 o n. re v o f geo p h y s an d s p ace p h ys ics , 1982 , 20 :219232fo mel s. shaping regula rizatio n in geop hysical esti2matio n p ro blems. e x p an de d a bs t racts o f 75 t h s e g m t g , 2005 ,

29、16731676姚姚著. 地球物理反演基本理論與應(yīng)用方法. 武漢 :中國地質(zhì)大學(xué)出版社 , 2005sheng j ia n2ming. h i g h resol ut i on seis m ic t om o g ra p h y w i t h gr t a n d e a rl y a r ri v al w ave f o rm t om o g ra p h y . di ssertatio n fo r docto r degree of u ta h u niver stiy ,2004woo dwa r d m j . wave2equatio n to mo grap h

30、y. geo p h y s2ics , 1992 , 57 ( 1) : 1526va sco d w. beyo nd ray to mo grap hy : wavepat hs andfre snel vol ume s. geo p h y s ics , 1995 , 60 ( 6 ) : 1790 1804zho u bi ng. cro sshole inver sio n wit h no r malized f ull 2wavefo r m amplit ude data . geo p h y s ics , 2003 , 68 ( 4) :13201330(本文編輯

31、:朱漢東)率較低而且占用內(nèi)存較大 。(2) 為了保證走時(shí)層析的線性假設(shè)條件和更有 效地利用先驗(yàn)信息 ,應(yīng)該采用正則化方法將先驗(yàn)信 息融入反演方程組中 ,替代傳統(tǒng)的外部約束方法 。 (3) 觀測系統(tǒng)關(guān)系到能夠采集哪些地震信息 ,進(jìn) 而影響到能夠反演出哪些信息 。理論上講 ,觀測系 統(tǒng)的炮點(diǎn) 、檢波點(diǎn)越密集 (小道距 、小炮距) 反演結(jié)果越接近于真實(shí) ,但小道距的層析結(jié)果相對于大道距 的改進(jìn)并不明顯 ,出于成本考慮 ,可以使用較大道距 (20 m) 數(shù)據(jù)進(jìn)行層析反演 。另外 ,由于層析的平均 效應(yīng) ,炮檢距與反演的淺層速度的分辨率成反比 ,為 了提高淺層的分辨率需要舍棄大炮檢距初至波信息 ,而只用

32、小炮檢距 (小于 2000 m 甚至更小) 信息 。 (4) 初始模型制約著反演的精度 ,初始模型越精 確越好 。通過理論模型試驗(yàn)發(fā)現(xiàn) ,如果沒有任何初 始模型的信息 ,那么利用梯度初始模型比常速初始 模型能得到更好的反演效果 ,收斂速度也較快 。一般來講速度是隨深度增加而增加的 ,所以梯度模型 比常速模型更準(zhǔn)確 、有效 。(5) 無論是否進(jìn)行正則化處理 ,都有必要對更新 過的模型進(jìn)行平滑處理 。平滑窗過大會降低反演分 辨率 ,平滑窗過小會降低層析的穩(wěn)定性 。理論模型試驗(yàn)表明 ,100200 m 的平滑窗在確保層析穩(wěn)定性 的同時(shí)可得較好的反演結(jié)果 。綜上所述 ,文中給出了一個(gè)比較好的層析參數(shù)

33、組合 。最后 ,通過實(shí)際資料的處理驗(yàn)證了初至波走 時(shí)層析方法解決近地表速度模型的有效性和本文總結(jié)的最佳參數(shù)組合選擇的合理性 。另外 ,走時(shí)層析的分辨率限制問題普遍存在 ,這 是由子波的高頻假設(shè)與觀測數(shù)據(jù)的局限性所致 。為 了能夠得到精度更高 、分辨率更高的反演結(jié)果 ,需要 從地震波傳播的物理本質(zhì)出發(fā) ,并考慮地震數(shù)據(jù)的多頻特征 ,將地震走時(shí) 、振幅和波形等觀測信息結(jié)合 起來 。當(dāng)然 ,這些信息的利用有待菲涅耳體層析成 像與波動(dòng)方程層析成像方法的進(jìn)一步發(fā)展與實(shí)用化 。 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 參 考 文 獻(xiàn) 21 1 sc

34、huster g t. fer mat s interf ero met ric p rinciple fo rtar get2o riented t raveltime to mo grap hy. geo p h y s ics ,2005 , 70 : 4750作 者王建民 教授級高級工程師 ,1961 年生 ;1982 年畢業(yè)于華東 石油學(xué)院物探專業(yè) , 1999 年畢業(yè)于長春科技大學(xué)地球 探測與信息專業(yè) ,獲碩士學(xué)位 ; 2004 年畢業(yè)于吉林大學(xué) 地球探測與信息技術(shù)專業(yè) ,獲博士學(xué)位 ;現(xiàn)為清華大學(xué)在 站博士后 。主要從事地震資料采集 、處理 、解釋及石油地 質(zhì)綜合研究工作 。陳學(xué)

35、強(qiáng) 高級工程師 ,1973 年生 ;1995 年獲中國地質(zhì)大學(xué)學(xué)士學(xué)位 ,2005 年獲中國石油大學(xué)碩士學(xué)位 。一直從事陸上地震勘探數(shù)據(jù)采集技術(shù)與方法研究工作 。印興耀 教授 ,博士生導(dǎo)師 ,1962 年生 ;1982 年在華東石油學(xué)院物探專業(yè)畢業(yè)后留校任教 ,1989 年和 1998 年分別獲蘭州大學(xué)無線電專業(yè)碩士學(xué)位和中國科學(xué)院地球物理所固體地球物理專業(yè)博士學(xué)位 。參與或主持完成了多項(xiàng)國家級和省部級科研課題 ,發(fā)表論文 30 余篇 ?，F(xiàn)在中國石油大學(xué) (華東) 地球資源與信息學(xué)院從事科研 、教學(xué)與管理工作 。王 棣 高級工程師 ,1969 年生 ;2006 年獲同濟(jì)大學(xué)固體地球物理專業(yè)博士

36、學(xué)位 。一直從事物探方法技術(shù)研究和地震資料疊前成像處理工作 ,發(fā)表論文 10 余篇 。李振春 教 授 , 博 士 生 導(dǎo) 師 , 1963 年 生 ; 中 國 石 油 大 學(xué) ( 華東) 地球物理系主任 。主 要從事地震波傳播與正演模擬 、地震成像與偏移速度分析 、多尺度地震資料聯(lián)合反演 、cf p2a v p 分析理論與方法等領(lǐng)域的教學(xué)與研究工作 。梁 鍇 1982 年生 ;2004 年畢業(yè)于中國石油大學(xué) ( 華東) 勘查技術(shù)與工程專業(yè) ,獲學(xué)士學(xué)位 ,2006 年于該校獲地球探測與信息技術(shù)專業(yè)碩士學(xué)位 ?，F(xiàn)為中國石油大學(xué) ( 華東) 地球資源與信息學(xué)院博士研究生 ,研究方向?yàn)楦飨虍愋越橘|(zhì)地震波傳播與成像 。王立歆 高級工程師 ,1971 年生 ;現(xiàn)為同濟(jì)大學(xué)海洋與地球科學(xué)學(xué)院在讀博士研究生 ,研究方向?yàn)榈卣鹳Y料處理及疊前成像方法研究 。郭 棟 高級工程師 , 1967 年生 ; 1991 年畢業(yè)于成都地質(zhì)學(xué)院石油地質(zhì)專業(yè) , 2006 年獲中國科學(xué)院廣州地球化學(xué)研究所構(gòu)造地質(zhì)學(xué)專業(yè)博士學(xué)位 。長期從事油氣地震地質(zhì)綜合研究工作 ,發(fā)表論文 20 余篇 。孫偉家 1984 年生 ; 現(xiàn)在中國科學(xué)院地質(zhì)與地球物理研究所攻讀碩士學(xué)位 ,研究方向?yàn)榈卣鹂碧椒椒ê蛻?yīng)用 。劉玉柱 講師 , 1979 年生 ; 2001 年和 2004 年先