陳 猛，魏 巍，刁永波，屈 元，朱亞東，王 鋒

(1.塔里木油田分公司 勘探開發(fā)研究院，庫爾勒 841000；2.川慶鉆探工程公司 地球物理勘探公司，成都 610213)

西昆侖山前巨厚黃土區(qū)層析靜校正深化應(yīng)用及效果分析

陳 猛1，魏 巍1，刁永波2*，屈 元2，朱亞東2，王 鋒2

(1.塔里木油田分公司 勘探開發(fā)研究院，庫爾勒 841000；2.川慶鉆探工程公司 地球物理勘探公司，成都 610213)

塔西南西昆侖山前帶油氣資源豐富，是近年塔里油田分公司地震勘探主攻區(qū)之一。然而，長期以來，由于受西昆侖山前復(fù)雜地表地下的地震地質(zhì)條件影響，該區(qū)地震資料成像效果較差，難以滿足地質(zhì)勘探需要。在眾多影響地震資料成像效果的因素中，西昆侖山前普遍存在的巨厚黃土及劇烈的地表起伏使得該區(qū)靜校正問題尤為突出，為后期去噪、速度分析和偏移成像處理帶來不利影響，是導(dǎo)致該區(qū)地震資料成像較差的主要原因之一。這里基于波動方程正演模擬結(jié)合地震勘探實(shí)踐，分析了西昆侖山前帶靜校正問題產(chǎn)生的主要原因，從而有針對性地提出了以層析靜校正為主的綜合靜校正方法，達(dá)到了精細(xì)刻畫巨厚黃土層復(fù)雜近地表結(jié)構(gòu)的目的。

層析靜校正； 西昆侖山前； 巨厚黃土； 起伏地表； 模擬正演

0 引言

塔西南西昆侖山前區(qū)位于塔里木板塊西南，毗鄰青藏高原西北部的帕米爾-西昆侖造山帶，包括塔西南前陸坳陷盆地和山前帶[1]。

西昆侖山前帶油氣資源豐富，然而黃土塬地震勘探是一個世界性難題，國內(nèi)、外雖進(jìn)行過大量研究，但進(jìn)展有限[2]。西昆侖山前帶發(fā)育的巨厚黃土、礫石層以及復(fù)雜的地表地下結(jié)構(gòu)，導(dǎo)致了地震波散射、衰減嚴(yán)重，反射波場雜亂，資料信噪比低。造成地震資料品質(zhì)極端低下，給正確認(rèn)識山前沖斷帶結(jié)構(gòu)造成了困難[3]。尤其是鹽下層可靠的反射信息更少，難以準(zhǔn)確落實(shí)構(gòu)造形態(tài)，嚴(yán)重制約了西昆侖山前帶的油氣勘探進(jìn)程[4]。

近年來，塔里木油田分公司在該區(qū)進(jìn)行了地震勘探的持續(xù)攻關(guān)，采集方法從傳統(tǒng)的二維單線采集發(fā)展為寬線采集，再到高密度寬線采集，進(jìn)而發(fā)展到束狀三維以及寬方位三維采集。隨著采集參數(shù)地不斷強(qiáng)化以及針對于該區(qū)的處理解釋技術(shù)的不斷成熟，西昆侖山前帶的地震資料成像品質(zhì)逐步得到提升。然而，受西昆侖山前巨厚黃土層以及劇烈起伏的近地表所引起的靜校正問題，仍然是該區(qū)地震資料成像處理面臨的一大挑戰(zhàn)。

1 研究區(qū)靜校正問題分析及解決思路

靜校正的目的是，消除由于地形和低降速帶變化以及炮點(diǎn)激發(fā)深度不同等因素引起的地震波傳播時間的畸變[5]。多年來，靜校正方法的研究和應(yīng)用一直是地震勘探中需要重點(diǎn)解決的難題。地球物理學(xué)家迪克斯[6]曾指出：“解決好靜校正就等于解決了地震勘探中幾乎一半的問題”。近年來，勘探工作者們在西昆侖山前持續(xù)的地震成像處理攻關(guān)中，總結(jié)出了針對巨厚黃土礫石區(qū)的一套以層析靜校正為主的綜合靜校正研究方法和思路，并取得了較好的效果。

1.1 起伏地表和黃土礫石引起的靜校正問題

西昆侖山前帶地表高程起伏較大，近地表覆蓋的黃土速度及厚度在縱橫向上分布嚴(yán)重不均，是導(dǎo)致該區(qū)復(fù)雜靜校正問題的兩大主要因素。

圖1為區(qū)內(nèi)柯東三維區(qū)塊的數(shù)字高程模型圖，從圖1中可以看到，在同一探區(qū)內(nèi)，海拔范圍為1 800 m～2 800 m，局部高差可達(dá)800 m。如此劇烈的地表起伏在進(jìn)行基準(zhǔn)面校正時，勢必會引起較大的靜校正量值應(yīng)用，加上受替換速度影響，時間域處理容易出現(xiàn)成像畸變。

圖1 區(qū)內(nèi)某三維工區(qū)數(shù)字高程模型Fig.1 In a 3D area digital elevation model

圖2為區(qū)內(nèi)相鄰二維測線的近地表結(jié)構(gòu)模型圖，地表低速帶速度從400 m/s～800 m/s分布，地層厚度最大處達(dá)500 m以上，最薄只有10 m，復(fù)雜的近地表結(jié)構(gòu)，必然會導(dǎo)致低降速帶靜校正量值出現(xiàn)較大的中高頻抖動，資料處理過程中，容易產(chǎn)生中長波靜校正問題。

圖2 區(qū)內(nèi)相鄰二維寬線近地表模型Fig.2 Near surface model of a two dimensional wide line in the region

為近一步驗(yàn)證西昆侖山前帶近地表結(jié)構(gòu)對靜校正的影響，采用基于聲學(xué)方程的模型正演進(jìn)行對比分析。速度模型依據(jù)區(qū)內(nèi)柯東某測線建立，近地表低降速帶厚度以實(shí)際調(diào)查結(jié)果而來(圖3)，模型局部高差達(dá)到300 m，對比近地表高速填充(2 500 m/s)與低速填充(400 m/s～800 m/s)正演單炮效果。從圖4可以看出，近地表模型填充速度為低速時的正演單炮，較高速填充的正演單炮初至抖動明顯劇烈。不僅如此，由于受近地表巨厚低速黃土層的吸收、衰減作用，低速填充正演單炮資料信噪比大幅降低，幾乎難見有效反射，與該地區(qū)實(shí)際采集資料較為吻合。通過圖4(a)、圖4(b)對比可見，在討論區(qū)山前帶巨厚低降速黃土層靜校正影響因素中，近地表復(fù)雜的低速帶結(jié)構(gòu)對資料的影響要大于高程的局部起伏帶來的影響。因此，在處理塔西南巨厚黃土礫石區(qū)靜校正問題中，應(yīng)該更加注意對低降速層近地表結(jié)構(gòu)的刻畫和分析，從而合理地進(jìn)行靜校正量值的計算。

圖3 正演模型表層低降速帶結(jié)構(gòu)Fig.3 Forward model for the low down and down speed zone structure of the model

1.2 采集方法強(qiáng)化導(dǎo)致靜校正問題復(fù)雜化

近年在西昆侖山前進(jìn)行的地震勘探中，為壓制側(cè)面及隨機(jī)干擾，增加覆蓋次數(shù)，提高資料信噪比，常以寬線加大組合(多串檢波器較大組合基距埋置)的采集方式為主；為接收廣角反射信息而采用的大偏移距采集方式也應(yīng)運(yùn)而生[6]。不可否認(rèn)，上述采集方案的運(yùn)用，在高陡構(gòu)造的山地地震資料采集中取得了較好的效果。但隨之而來的卻對室內(nèi)精細(xì)的靜校正處理帶來新的挑戰(zhàn)，主要體現(xiàn)在：

圖4 近地表不同填充速度的正演單炮Fig.4 Forward modling shot record with different replacement velocity near the surface(a)高速填充；(b)低速填充

1)寬線采集使炮線與接收線距增加(最大間距達(dá)120 m以上)，在室內(nèi)資料處理過程中一般使用橫向大面元方式進(jìn)行疊加，位于同一面元內(nèi)的炮檢點(diǎn)分配相同的一致性靜校正量值，從而導(dǎo)致靜校正模糊化；大組合接收存在同樣的問題，即從面上獲得的地震信息使用點(diǎn)上的靜校正量值，這種模糊化的靜校正應(yīng)用在低降速帶巨厚及橫向分布不均的近地表?xiàng)l件下，是不合理的。

2)較高的覆蓋次數(shù)(500次以上)，有利于剖面信噪比的提高 ，但是在后期進(jìn)行地表一致性剩余靜校正計算時，卻會導(dǎo)致靜校正量分解困難，最終形成的地表一致性靜校正量值較小，難以進(jìn)行有效的疊加成像迭代。

3)不論是基于模型的基礎(chǔ)靜校正量計算，還是后期的地表一致性的剩余靜校正計算，都是基于地表一致性假設(shè)的，即每一個物理點(diǎn)享有唯一的靜校正量值。然而，如果同一物理點(diǎn)在近偏移距和遠(yuǎn)偏移距使用同一靜校正量值，會使真實(shí)的走時產(chǎn)生畸變，這種畸變隨偏移距的增加而凸顯，最終導(dǎo)致速度及成像誤差增大，即產(chǎn)生中長波長靜校正問題。排列過長的大偏移距接收在地形起伏較大的山地地表情況下，會使得其與靜校正的地表一致性假設(shè)矛盾加劇。

1.3 西昆侖山前靜校正問題帶來的啟示

綜上所述，面對近年來西昆侖山前采集參數(shù)的不斷強(qiáng)化，需要配套的靜校正手段進(jìn)行相應(yīng)處理，在以提高成像信噪比為第一前提下，還要充分考慮如何保證靜校正計算的合理性，從而保證成像誤差可控。

當(dāng)前折射靜校正與層析靜校正，是地震資料處理中最常用的兩種靜校正方法。折射靜校正適合近地表存在穩(wěn)定的折射層，表層速度和厚度縱橫向變化不太劇烈的地區(qū)，而層析靜校正在算法上較前者復(fù)雜許多，適合地形起伏較大，近地表結(jié)構(gòu)復(fù)雜，速度變化劇烈的地區(qū)。因而，它更能夠解決因速度橫向劇烈變化和速度倒轉(zhuǎn)引起的靜校正問題[7]，該方法能夠直接利用實(shí)際地震記錄的初至?xí)r間，傳播距離短，信噪比高、畸變小，易于識別[8]。西昆侖山前帶的原始地震資料信噪比極低，但是初至波旅行時卻可以相對較為準(zhǔn)確地拾取，根據(jù)地表第四系地層的沉積機(jī)理可知，區(qū)內(nèi)近地表低降速帶受緩慢壓實(shí)作用影響，縱、橫向上速度呈漸變趨勢，即近地表結(jié)構(gòu)屬于連續(xù)介質(zhì)，初至波符合回折波傳播原理，通過拾取初至波旅行時進(jìn)行層析反演，進(jìn)而刻畫近地表結(jié)構(gòu)，是較為實(shí)用可行的。折射靜校正由于自身的局限性，顯然不能很好適應(yīng)西昆侖山前帶的復(fù)雜近地表結(jié)構(gòu)。

針對西昆侖山前帶資料特點(diǎn)，總結(jié)出一套針對性的綜合靜校正處理思路，解決巨厚黃土礫石區(qū)的山前高精度靜校正處理顯得十分必要(圖5)，這里將著重對層析靜校正技術(shù)的深化應(yīng)用進(jìn)行論述。

圖5 塔西南巨厚黃土區(qū)綜合靜校正思路Fig.5 Comprehensive static correction of the Loess area in the southwest of Tarim basin

2 高精度層析反演近地表建模思路

2.1 初至波層析反演的基本原理

初至波表層模型層析反演的原理，是利用地震波射線的走時和路徑反演介質(zhì)速度結(jié)構(gòu)的一種高精度反演方法[9]。因此，在反演中所利用的地震波類型決定了對反演模型的適應(yīng)性。在任意表層模型反演中，初至波包括直達(dá)波、回折波、折射波以及幾種波組合后首先到達(dá)地表的波。由于直達(dá)波主要體現(xiàn)了均勻介質(zhì)模型，回折波主要體現(xiàn)連續(xù)介質(zhì)模型，而折射波主要體現(xiàn)層狀介質(zhì)模型，通過三者的組合以及層析法對介質(zhì)橫向變化的適應(yīng)性[9]，從而使得該方法適應(yīng)任意表層模型的反演問題。

初至波表層模型層析反演首先是將一個非線性問題簡化為一個線性方程組,如式(1)所示，然后用迭代逼近的方法求解。

AΔV=ΔT

(1)

式中:ΔT為m維初至波初至?xí)r間與模型正演初至波走時的差值列向量；ΔV為n維模型參數(shù)修改量列向量；A為m×n維與射線路徑及模型參數(shù)有關(guān)的矩陣,也稱為微分系數(shù)矩陣。式( 1)可寫成如下的分量形式

Δti=aijΔvj

(2)

式中微分系數(shù)矩陣為

(3)

初至波層析反演的基本步驟為：

1)初至波旅行時的精確拾取。

2)建立初始模型，進(jìn)行正演旅行時生成。

3)反演計算，計算計算旅行時與實(shí)際旅行時之差，進(jìn)行模型修正。

4)重復(fù)步驟2、3，對初始模型進(jìn)行迭代修正，直至擬合差符合要求。

5)計算模型靜校正量。

影響層析反演質(zhì)量的主要因素有：①野外空間采樣密度；②初始模型的合理性；③計算網(wǎng)格大小等。經(jīng)過多年的實(shí)踐，在該區(qū)采用了一套針對性的層析反演方法進(jìn)行精細(xì)的近地表建模及基礎(chǔ)靜校正量計算，取得了一定效果。

2.2 擬三維層析反演

擬三維層析反演主要針對寬線資料處理。圖6為該區(qū)典型的3炮2線的寬線采集觀測系統(tǒng)，兩條炮線最大相距達(dá)120 m。如果使用常規(guī)的二維方式進(jìn)行層析反演，將忽略低速帶及高程橫向變化帶來的影響，將處于橫向同一大面元的炮檢點(diǎn)歸入單獨(dú)一個點(diǎn)進(jìn)行計算，反演精度必然受到影響。使用三維層析反演方式建立針對寬線的初始層析模型進(jìn)行反演計算，可以較好地解決這個問題。具體做法是，縱向面元尺寸不變，定義一個小于最小線距的數(shù)值為橫向面元尺寸，建立矩陣方程進(jìn)行擬三維層析迭代計算。

圖6 區(qū)內(nèi)典型寬線觀測系統(tǒng)圖示Fig.6 A typical wide line geometry

圖7為常規(guī)二維方式進(jìn)行層析反演求取的近地表模型，圖8為三維小面元方式求取的近地表模型。對比圖7、圖8可以看出，三維層析反演模型精細(xì)地刻畫出了縱、橫向方向的低降速變化規(guī)律。理論上講，通過三維模型進(jìn)行低速帶剝離計算靜校正更為合理。從圖9(a)、圖9(b)對比看，三維反演計算的靜校正量應(yīng)用于水平疊加成像效果要優(yōu)于二維。

圖7 對寬線進(jìn)行常規(guī)二維方式層析反演的近地表模型Fig.7 A near surface model for the conventional two dimensional tomography inversion is carried out

2.3 區(qū)域連片層析建模

區(qū)域連片層析建模主要是針對二維測線的層析反演，主測線與聯(lián)絡(luò)線方向在擬合計算上存在差異，往往會導(dǎo)致反演模型在交結(jié)點(diǎn)位置出現(xiàn)厚度及速度的差異，最終應(yīng)用校正量會引起剖面不閉合的問題。

在對區(qū)內(nèi)多條二維測線進(jìn)行層析反演后，將反演模型統(tǒng)一測網(wǎng)進(jìn)行修正，使交點(diǎn)位置的速度模型結(jié)構(gòu)趨于一致，然后進(jìn)行連片的靜校正值計算，最終可以較好地解決二維測網(wǎng)存在閉合差的問題(圖10)。

圖8 對寬線進(jìn)行擬三維方式層析反演的近地表模型反演的近地表模型Fig.8 Near surface model for the quasi three dimensional model of the wide line(a)二維展示；(b)三維展示

圖9 寬線二、三維方式層析反演靜校正量應(yīng)用疊加對比反演的近地表模型Fig.9 Comparison between the static correction and the application of the static correction by the wide line two and three dimensional method(a)二維方式；(b)三維方式

圖10 區(qū)域連片層析建模Fig.10 Area contiguous tomographic modeling

2.4 低速測井資料約束的層析反演建模

由于受空間采樣的限制，近地表的低速特征通過傳統(tǒng)的層析反演難以精細(xì)刻畫，近地表速度普遍比實(shí)際偏高。研究區(qū)內(nèi)有較多的微測井、小折射等低降速帶調(diào)查成果資料，這些資料是控制地表以下0 m～50 m范圍內(nèi)低速帶特征的最直觀最可靠的資料。因此，在建立初至波層析反演的初始模型時，利用這些低速測井成果進(jìn)行約束建模，對淺近地表的速度及厚度特征進(jìn)行有效控制，能夠保證層析反演更為有效的收斂迭代。

圖11為直接使用大炮初至進(jìn)行層析反演后的表層模型(未約束)與使用低速測井成果進(jìn)行約束層析反演后的表層模型(井約束)示意圖。對比圖11(a)、圖11(b)可以看出，約束后的淺近地表的低降速帶信息更為豐富，與地層實(shí)際速度分布更吻合，所建反演模型更為準(zhǔn)確合理。通過對比圖12(a)、圖12(b)可以看出，有井約束的剖面成像效果更好。

2.5 高精度層析反演

由式(1)可以知道，建立層析反演初始模型時，初始模型網(wǎng)格的尺寸影響反演模型的細(xì)節(jié)豐富程度?？s小野外采集道距，從而增加初至旅行時空間采樣密度，可以豐富的旅行時參與層析反演的擬合迭代計算，便于快速收斂，同時更利于精細(xì)刻畫近地表的低速特征。

2.5.1 小道距層析調(diào)查

近年來，在該區(qū)試驗(yàn)性地采集了多條單點(diǎn)小道距(10 m道距)的層析調(diào)查資料，以期通過對高空間采樣密度的旅行時進(jìn)行層析反演獲得精細(xì)可靠的近地表模型，同常規(guī)大炮初至(30 m道距為主)層析反演對比，小道距層析模型刻畫出近地表更為豐富的低降速細(xì)節(jié)特征(圖13)，利于更精細(xì)的靜校正量值計算。

圖11 低速測井資料約束的層析反演對比圖Fig.11 Contrast chart of constrained tomography inversion of low velocity log data(a)井約束；(b)未約束

圖12 使用低速測井資料約束反演疊加效果對比Fig.12 Comparison of the effect of constrained inversion with low velocity log data(a)未約束；(b)井約束

圖13 不同道距的層析反演模型對比Fig.13 Comparison of different models of different trace spacing(a)10 m道距；(b)30 m道距

2.5.2 小網(wǎng)格層析反演

小道距層析要求在大炮采集線的同樣位置重新布設(shè)炮檢點(diǎn)進(jìn)行采集，增加了野外生產(chǎn)成本，特別是針對三維工區(qū)的表層調(diào)查，如果大面積使用小道距層析調(diào)查方式進(jìn)行采集，而獲得的成像效果改善有限，性價比不高，可操作性不強(qiáng)。

在實(shí)際生產(chǎn)應(yīng)用中，層析反演面元按經(jīng)驗(yàn)通常定為道距的一半。實(shí)際上，根據(jù)層析反演的算法可知，反演網(wǎng)格的尺寸定義越小，對連續(xù)介質(zhì)的速度刻畫越為精細(xì)，計算量也就越大。

圖14 50 m道距25 m面元反演層析模型Fig.14 50 m trace spacing inversion tomography model of 25 m cell size

圖15 50 m道距5 m面元反演層析模型Fig.15 50m trace spacing inversion tomography model of 5m cell size

圖14為通過50 m道距采集的初至旅行時反演獲得的層析模型，反演面元為25 m，圖15為50 m道距5 m面元反演模型，圖16為10 m道距5 m面元反演模型。對比圖14～圖16可以看出，即使是較大的空間采樣，使用較小的面元進(jìn)行反演計算，仍能獲得與較小道距采樣的較為近似的效果，甚至于反演迭代收斂度也基本吻合(圖17)，這種吻合表示兩者反演效果基本接近。計算的模型靜校正量應(yīng)用于單炮初至?xí)r(圖18)，小面元反演效果明顯要優(yōu)于大面元反演，初至校正后更為光滑。

圖16 10 m道距5 m面元反演層析模型Fig.16 10 m trace spacing inversion tomography model of 5m cell size

圖17 不同反演方式的迭代收斂度對比Fig.17 Iterative convergence of different inversion methods

隨著計算機(jī)運(yùn)算能力的增強(qiáng)，使得較小面元的層析反演計算能順利完成，結(jié)合本區(qū)大量的寬線采集資料(以30 m道距為主)，三維資料觀測面元也較大(25 m×25 m)，使用較小面元的初始模型進(jìn)行層析反演，既能較大程度地降低地震采集成本，又能充分挖掘現(xiàn)有資料的潛力，利于該區(qū)精細(xì)的近地表建模。

圖18 不同面元層析反演計算靜校正單炮應(yīng)用效果對比Fig.18 Comparison of the effect of different cell size tomographic inversion on the application of static correction(a)原始單炮；(b)25 m面元；(c)5 m面元

3 應(yīng)用效果分析

綜合以上方法將層析反演進(jìn)行優(yōu)化應(yīng)用，能夠更為精細(xì)地對近地表進(jìn)行刻畫，從而計算出更為合理的基礎(chǔ)靜校正量值。

從區(qū)內(nèi)單炮應(yīng)用效果看(圖19)，使用優(yōu)化層析反演后計算靜校正量的單炮初至更為平滑，局部的高頻抖動能夠得到一定控制，從而為后期各項(xiàng)預(yù)處理及成像打好資料基礎(chǔ)。

對比圖20、圖21可以看出，在巨厚黃土區(qū)內(nèi)，如果能夠?qū)乇斫Y(jié)構(gòu)進(jìn)行更為精細(xì)合理地刻畫，計算出靜校正量對于后期成像是有明顯的積極作用。

圖19 優(yōu)化層析反演方法精細(xì)建模求取靜校正量單炮應(yīng)用效果對比Fig.19 Comparison of the application effect of the seismic record between optimization and non optimization tomographic inversion(a)原始單炮；(b)常規(guī)層析靜校正量應(yīng)用后單炮；(c)優(yōu)化層析建模靜校正量應(yīng)用單炮

圖20 常規(guī)層析反演計算基礎(chǔ)靜校正量應(yīng)用后疊加剖面Fig.20 Profile of conventional tomographic inversion for the application of static correction

圖21 綜合優(yōu)化層析反演方法后計算基礎(chǔ)靜校正量后疊加剖面Fig.21 A comprehensive optimization method for the calculation of the superimposed profile after static correction

4 結(jié)論與認(rèn)識

西昆侖山前帶前期基礎(chǔ)靜校正的求取，主要依賴于對近地表結(jié)構(gòu)的精細(xì)刻畫，尤其是在巨厚黃土礫石層的靜校正處理中，速度縱、橫向變化較大的復(fù)雜近地表?xiàng)l件下,層析靜校正是切實(shí)可以行的靜校正方法。

1)針對厚黃土區(qū)近地表的精細(xì)建模，是該區(qū)解決靜校正問題的基礎(chǔ)。在資料信噪比不高情況下，利用配套的高精度層析反演技術(shù)，能較為準(zhǔn)確地對近地表結(jié)構(gòu)進(jìn)行刻畫。

2)在現(xiàn)有資料條件下，層析反演應(yīng)該綜合運(yùn)用微測井、小折射等低速測井資料，有助于更為精細(xì)的反演計算，從而獲得較為合理的靜校正量值。

3)小面元、小道距層析反演將是高精度層析反演發(fā)展的一個方向，在野外空間采樣密度有限的條件下，縮小初始模型的網(wǎng)格尺寸，同樣利于反演精度的提高。

[1]謝會文,朱亞東,曾昌明.西昆侖山前地區(qū)白堊系分布研究[J].沉積與特提斯地質(zhì),2014,34(3):57-58.XIE H W,ZHU Y D,ZENG C M.The piedmont area of West KunLun cretaceous distribution research[J].Sedimentary Geology And Tethyan Geology,2014,34(3):57-58.(In Chinese)

[2]羅仁澤,梁黎明,吳希光.寬線大組合理論及其在黃土塬地震采集中的應(yīng)用[J].天然氣工業(yè),2009,29(2):53-55.LUO R Z,LIANG L M,WU X G.The Application of the theory of wide line array and seismic acquisition in Loess Plateau[J].Natural Gas Industry,2009,29(2):53-55.(In Chinese)

[3]張瑋,漆家福,雷剛林.塔西南坳陷西昆侖山前沖斷帶的收縮構(gòu)造變形模式[J].新疆石油地質(zhì),2010,31(6):567-560.ZHANG W,QI J F,LEI G L.The Southwest Tarim Basin KunLun foreland thrust tectonic belt shrinkage deformation model[J].Xinjiang Petroleum Geology,2010,31(6):567-560.(In Chinese)

[4]劉依謀,梁向豪,孫建華.塔里木盆地西南坳陷巨厚黃土區(qū)地震采集技術(shù)[J].新疆石油地質(zhì),2009，30( 2):255-257.LIU Y M,LIANG X H,SUN J H.Seismic acquisition technology of giant thick loess area in southwest depression of Tarim Basin[J].Xinjiang Petroleum Geology,2009，30( 2):255-257.(In Chinese)

[5]孟慶利.基于大炮初至的層析靜校正算法研究與實(shí)現(xiàn)[D].北京:中國石油大學(xué),2011.MENG Q L.The research and implementation of the algorithm of tomography static correction based on the initial arrival of shot record[D].Beijing:China University of Petroleum,2011.(In Chinese)

[6]鐘本善，周熙襄.中國西部地區(qū)地震勘探的靜校正問題[J].物探化探計算技術(shù)，1999，21(4)：358-366.ZHONG B S,ZHOU X X.The static correction prollen of seismic exploration in westerm china[J].Computing Techniques for Geophysical and Geochemical Eexploration,1999,21(4):358-366.(In Chinese)

[7]莊道川,曾慶才,張亞斌.超大偏移距地震資料處理方法[J].西南石油學(xué)院報,2003,25(2):11-13.ZHUANG D C,ZENG Q C,ZHANG Y B.Large offset seismic data processing method[J].Journal of Southwest Petroleum Institute,2003,25(2):11-13.(In Chinese)

[8]于豪.折射波靜校正與層析靜校正技術(shù)適用性分析[J].地球物理學(xué)進(jìn)展,2012,27(6)：2577-2579.YU H.Applicability analysis of refraction statics and tomographic static correction technique[J].Progress in Geophysics,2012,27(6)：2577-2579.(In Chinese)

[9]周國婷，潘冬明，牛歡.層析靜校正方法研究與應(yīng)用[J].物探與化探,2012,36(5):802-805.ZHOU GT,PAN D M,NIU H.Research and application of tomography static correction method[J].Geophysical and Geochemical Exploration,2012,36(5):802-805.(In Chinese)

[10]李錄明,羅省賢,趙波.初至波表層模型層析反演[J].石油地球物理勘探,2000,35(5):559-564.LI L M,LUO S X,ZHAO B.Surface model tomography inversion by first wave[J].Oil Geophysical Prospecting,2000,35(5):559-564.(In Chinese)

West Kunlun piedmont thick loess area tomostatics deepen the application and effect analysis

CHEN Meng1,WEI Wei1,DIAO Yong-bo2*,QU Yuan2,ZHU Ya-dong2,WANG Feng2

(1.Tarim Oilfield Company exploration and development Research Institute,Korla 841000，China；2.CCDC Geophysical Company,Chengdu 610213，China)

Rich oil and gas resources in west Kunlun piedmont belt,southwest of Tarim basin,is one of the main seismic exploration area of Tarim oilfield company in recent years.However,due to geological conditions underground complex surface effects and poor seismic imaging area in the west Kunlun piedmont,it is difficult to meet the needs of geological exploration for a long time.Among the many factors affecting the seismic data imaging,the western Kunlun piedmont ubiquitous thick loess and dramatic ups and downs so that the area of surface statics problem is particularly prominent for late denoising,migration velocity analysis and imaging processing adversely affected area is one of the main causes of poor seismic imaging.Based on the wave equation forward modeling of seismic exploration practice combined analysis of the main reasons the west kunlun piedmont zone statics problems,which put forward a comprehensive static correction method based tomostatics,reach fine portray thick the purpose of the near-surface layer of loess complex structure.

tomostatics; west Kunlun piedmont; thick loess; undulating surface; forward modeling

2015-09-30 改回日期：2015-10-23

陳猛(1969-)，男，高級工程師，現(xiàn)從事地震資料處理、變速成圖、工程地質(zhì)等領(lǐng)域研究,E-mail：chengmen-tlm@petrochina.com.cn。

*通信作者：刁永波(1982-),男，工程師，現(xiàn)主要從事陸上復(fù)雜構(gòu)造低信噪比條件下的地震資料處理成像研究工作,E-mail：diaoyb_wt@cnpc.com.cn。

1001-1749(2016)06-0796-09

P 631.4