張永峰，王　鵬，張亞兵，胥小萍，蘇　勤，徐興榮

（1.中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司，西安710065；2.中國(guó)石油勘探開(kāi)發(fā)研究院西北分院，蘭州730020；3.中國(guó)石油吉林油田分公司勘探開(kāi)發(fā)研究院，吉林松原131200）

基于變差函數(shù)的高精度靜校正融合技術(shù)及其應(yīng)用

張永峰1，王鵬2，張亞兵3，胥小萍3，蘇勤2，徐興榮2

（1.中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司，西安710065；2.中國(guó)石油勘探開(kāi)發(fā)研究院西北分院，蘭州730020；3.中國(guó)石油吉林油田分公司勘探開(kāi)發(fā)研究院，吉林松原131200）

地震反射波時(shí)距曲線受復(fù)雜地表地質(zhì)結(jié)構(gòu)的影響產(chǎn)生畸變而不能同相疊加，需要對(duì)地震資料做靜校正處理。不同的靜校正方法都有各自的適用條件：高程靜校正無(wú)法消除低降速帶的影響；模型靜校正受限于追蹤的深度和排列的長(zhǎng)度，難以全面反映低速帶底界的形態(tài)；折射靜校正僅適用于地表較平緩、表層速度橫向均勻性較好且有明顯折射界面的地區(qū)；層析靜校正適應(yīng)任意表層模型的反演，但是反演結(jié)果不穩(wěn)定。在地震勘探地區(qū)地貌日益復(fù)雜的今天，優(yōu)選的某一種靜校正方法無(wú)法有效解決所有的問(wèn)題。鑒于上述問(wèn)題，提出了一種適用于多種地表類型復(fù)雜區(qū)域地震資料處理的基于變差函數(shù)擬合重構(gòu)的高精度三維靜校正方法，該方法能夠?qū)崿F(xiàn)不同靜校正量的擬合和重構(gòu)，并很好地解決復(fù)雜地表?xiàng)l件下不同靜校正方法優(yōu)勢(shì)擬合的難題，從而可有效提高復(fù)雜地表地形條件的三維地震資料的成像質(zhì)量，保證構(gòu)造形態(tài)的可靠性。

高精度靜校正；變差函數(shù)；擬合重構(gòu)；復(fù)雜地表地區(qū)

0　引言

近年來(lái)，地震勘探向山地、戈壁、沙漠及黃土塬等地表?xiàng)l件復(fù)雜的地區(qū)發(fā)展［1］，而復(fù)雜的地表?xiàng)l件給地震資料處理帶來(lái)很大影響。由于地震波在傳播過(guò)程中受復(fù)雜地表地質(zhì)結(jié)構(gòu)的影響會(huì)產(chǎn)生畸變，進(jìn)而影響速度分析和疊加成像，因此需要對(duì)地震資料做靜校正處理［2］。中國(guó)西部?jī)?nèi)陸盆地的靜校正問(wèn)題是決定地震剖面成像效果的主要因素，也是地震勘探取得突破的技術(shù)瓶頸，地球物理學(xué)家Dix曾說(shuō)：“解決好靜校正就等于解決了地震勘探中幾乎一半的問(wèn)題”［3］。目前，成熟的靜校正方法有模型靜校正、高程靜校正、折射靜校正和層析反演靜校正等，這些方法均基于某種特定的近地表模型假設(shè)和野外采集數(shù)據(jù)特征，適用于不同的地表地質(zhì)條件［4］。在包括山地、戈壁、沙漠及黃土塬等多種地表類型的復(fù)雜地區(qū)的實(shí)際應(yīng)用過(guò)程中，只能優(yōu)選某一種靜校正方法，而無(wú)法實(shí)現(xiàn)不同靜校正方法優(yōu)勢(shì)的擬合，若只是簡(jiǎn)單拼接靜校正量，則在2種靜校正的邊界位置，會(huì)出現(xiàn)地震剖面同相軸錯(cuò)位或抖動(dòng)，即不能很好地同向疊加的短波長(zhǎng)靜校正問(wèn)題，或者出現(xiàn)與實(shí)際地質(zhì)高程變化不符的剖面構(gòu)造形態(tài)，剖面淺、中、深層構(gòu)造形態(tài)不一致及構(gòu)造畸變等長(zhǎng)波長(zhǎng)靜校正問(wèn)題［5］。鑒于此，筆者首先分析常用靜校正方法的適用性，針對(duì)多種地表類型復(fù)雜區(qū)域地震資料處理，提出應(yīng)用變差函數(shù)進(jìn)行擬合重構(gòu)的高精度三維靜校正方法，以期實(shí)現(xiàn)不同靜校正量的擬合和重構(gòu)，并很好地解決復(fù)雜地表?xiàng)l件下不同靜校正方法優(yōu)勢(shì)擬合的難題，進(jìn)而有效提高復(fù)雜地表?xiàng)l件的三維地震資料的成像質(zhì)量，保證構(gòu)造形態(tài)的可靠性。

1　不同靜校正方法的適應(yīng)性分析

西部探區(qū)是以山地、戈壁、沙漠及黃土塬為主的近地表結(jié)構(gòu)復(fù)雜區(qū)。黃土塬溝塬交錯(cuò)，低速帶不穩(wěn)定；戈壁、沙漠區(qū)砂礫巖層厚度變化大，低降速帶厚度及速度變化劇烈；山地地表起伏劇烈，山體陡峭，老地層出露，折射層不穩(wěn)定。復(fù)雜地表的靜校正問(wèn)題往往歸結(jié)為近地表層速度的求解問(wèn)題，地表越復(fù)雜，靜校正問(wèn)題解決的難度就越大。常用的靜校正方法主要為高程靜校正、模型靜校正、折射靜校正及層析靜校正等。高程靜校正無(wú)法消除低降速帶的影響，而模型靜校正利用小折射、微測(cè)井解釋結(jié)果建立模型，但該方法受限于追蹤的深度和排列的長(zhǎng)度，難以全面反映低速帶底界的形態(tài)，靜校正不徹底。在實(shí)際應(yīng)用中，折射靜校正和層析靜校正在符合其適用條件的地區(qū)均有較好的應(yīng)用效果，這是對(duì)高精度三維靜校正方法進(jìn)行優(yōu)勢(shì)擬合的基礎(chǔ)。

1.1層析靜校正

層析靜校正技術(shù)是一種利用單炮初至進(jìn)行近地表速度反演的方法［6-7］，該方法假設(shè)近地表速度由多個(gè)速度單元組成，每個(gè)單元的速度不同且為常數(shù)。地震波的旅行時(shí)間是對(duì)地下介質(zhì)慢度函數(shù)沿著波的傳播射線路徑進(jìn)行線性積分，可表示為

式中：s（x，z）為地下介質(zhì)的慢度函數(shù)；dl為波的射線路徑的微分；t為地震波從震源s到檢波點(diǎn)r的走時(shí)。

式（1）離散后的矩陣形式表示為

式中：T為所有炮點(diǎn)到檢波點(diǎn)的旅行時(shí)矩陣；S為地下介質(zhì)的慢度矩陣；A為與地震波傳播射線路徑有關(guān)的距離矩陣。

層析運(yùn)算包括一個(gè)正演過(guò)程（計(jì)算每個(gè)炮點(diǎn)-檢波點(diǎn)對(duì)的旅行時(shí)間）和一個(gè)反演過(guò)程（根據(jù)初至剩余時(shí)間交互更新速度模型）。完整的層析靜校正算法主要包括以下幾個(gè)步驟：①大炮初至拾取；②成像域網(wǎng)格化；③射線追蹤和分割；④剩余時(shí)間（誤差）計(jì)算；⑤更新速度，減小誤差，得到更為合理的速度模型。

基于初至信息的非線性層析反演靜校正方法具有很明顯的優(yōu)勢(shì)，其初至波不受折射波的局限，還可以是直達(dá)波、透射波及繞射波，反演所利用的地震波類型決定了對(duì)反演模型的適應(yīng)性。另外，該方法假設(shè)近地表速度由多個(gè)速度單元組成，使其可以適應(yīng)任意表層模型的反演［8-9］，而不受地形起伏、橫向速度變化和地下界面傾斜等因素的影響，特別是在存在高陡構(gòu)造以及老地層出露的地區(qū)，層析靜校正表現(xiàn)出了更大的優(yōu)勢(shì)。不過(guò)該方法在解決一些問(wèn)題的同時(shí)又帶來(lái)了另一些問(wèn)題，如由于介質(zhì)被網(wǎng)格化為一系列單元，引入了大量的未知量，使反演結(jié)果不穩(wěn)定，需要進(jìn)行間接的正則化約束［10-11］。

1.2折射靜校正

折射靜校正方法基于層狀介質(zhì)假設(shè)［12］，在假定折射界面穩(wěn)定的條件下，利用拾取的折射波初至建立地下折射面模型，然后求出低速帶的速度和厚度，再通過(guò)給定的替換速度和統(tǒng)一基準(zhǔn)面求出靜校正量，并最終消除因地形及低速帶速度和厚度變化等因素所引起的靜校正問(wèn)題。假設(shè)只存在一個(gè)折射層的近地表模型，tAB表示經(jīng)過(guò)炮點(diǎn)s激發(fā)檢波點(diǎn)r所接收的折射旅行時(shí)：

設(shè)A點(diǎn)的延遲時(shí)為tA，B點(diǎn)的延遲時(shí)為tB，則式（1）可寫(xiě)為

式（4）就是折射波初至方程。對(duì)于每一個(gè)折射波初至都可以建立這樣一個(gè)初至方程，那么應(yīng)用于整個(gè)工區(qū)就可以得到一個(gè)折射波初至方程組。在近地表存在n個(gè)折射層的情況下，為了最佳地確定每一個(gè)折射層，需要記錄來(lái)自不同炮檢距的數(shù)據(jù)，再根據(jù)不同的炮檢距范圍逐層反演計(jì)算。

由于其假設(shè)條件是水平層狀介質(zhì)，因此折射靜校正適用于地表較平緩、表層速度橫向均勻性較好且有明顯折射界面的地區(qū)，然而其缺點(diǎn)也很明顯，即基于假設(shè)，為了獲得來(lái)自每一層的折射波，速度必須隨深度增大，所以在地層存在速度反轉(zhuǎn)或者低降速帶過(guò)厚時(shí)，都將直接影響該方法的適用性，并進(jìn)一步影響反演模型和靜校正量的精度。多年的實(shí)踐證明，初至折射靜校正方法存在低頻成分處理欠佳以及出現(xiàn)假構(gòu)造等現(xiàn)象。針對(duì)以上問(wèn)題，在實(shí)際應(yīng)用時(shí)可通過(guò)計(jì)算得到包含高、低頻分量的初至折射波靜校正量，并實(shí)現(xiàn)高、低頻靜校正量的分離，進(jìn)而結(jié)合地質(zhì)上的判別標(biāo)準(zhǔn)應(yīng)用變差函數(shù)進(jìn)行擬合重構(gòu)。

2　變差函數(shù)的原理

2.1基本原理

數(shù)學(xué)家Matheron在1962年創(chuàng)立的變差函數(shù)理論兼顧了區(qū)域化變量的隨機(jī)性和數(shù)據(jù)的空間結(jié)構(gòu)特征，是一個(gè)表征區(qū)域化變量隨機(jī)性和結(jié)構(gòu)性的有效工具，可用于研究地質(zhì)變量的空間相關(guān)性［13］。區(qū)域化變量Z（x）在點(diǎn)x和x+h處的參數(shù)值Z（x）與Z（x+h）之差的方差的一半為Z（x）在h方向上的變差函數(shù)，可表述為

式中：h為滯后距離，m；E表示數(shù)學(xué)期望；Z（x）為在位置x處的變量值；Z（x+h）為在位置x偏離h處的變量值。

定義Z（x）在x和x+h兩點(diǎn)的兩階混合中心矩為協(xié)方差函數(shù)，可表示為

協(xié)方差函數(shù)C（h）表示任意兩點(diǎn)的相關(guān)程度。兩點(diǎn)間的距離為h，當(dāng)h越大，這種相關(guān)程度越小，當(dāng)h充分大時(shí)，則完全不相關(guān)，即存在一個(gè)定值a，當(dāng)h≥a時(shí)，有C（h）=0。這表明距離大于a的2個(gè)區(qū)域化變量是不相關(guān)的。由式（6）可得

式中：a稱為y（h）的變程，它反映了區(qū)域化變量影響的范圍。

假設(shè)N（h）是間距為h的所有點(diǎn)對(duì)的總數(shù)，則變差函數(shù)可通過(guò)下式計(jì)算：

式中：N（h）為間距為h的點(diǎn)對(duì)數(shù)；Z（xi）為某位置的變量值；y（h）為實(shí)驗(yàn)變差函數(shù)。

圖1　變差函數(shù)示意圖Fig.1　Schematic diagram of variogram and its parameters

以h為橫坐標(biāo)，y（h）為縱坐標(biāo)所得到的一組［h，y（h）］點(diǎn)稱為變差函數(shù)圖（圖1）。變差函數(shù)y（h）隨滯后距離h變化的各項(xiàng)特征表達(dá)了區(qū)域化變量的各種空間變異性質(zhì)，這些特征可通過(guò)變差圖的各項(xiàng)參數(shù)，如變程、塊金值及基臺(tái)值來(lái)表示。變程的地震學(xué)意義在于其不僅反映靜校正變量在某一方向上變化，還反映對(duì)應(yīng)的地表類型在某一方向上的延伸程度及地表的規(guī)模和平面非均質(zhì)性。變程長(zhǎng)，表明靜校正量空間分布的相關(guān)性好，即非均質(zhì)性相對(duì)較弱；反之，則非均質(zhì)性相對(duì)較強(qiáng)。塊金效應(yīng)用以描述在很短的距離內(nèi)兩點(diǎn)間變量值的變化。基臺(tái)值反映變量在研究范圍內(nèi)的變化幅度，基臺(tái)值越小說(shuō)明數(shù)據(jù)的波動(dòng)程度越小，參數(shù)變化的幅度就越小。

2.2變差函數(shù)分析

如何合理利用變差函數(shù)進(jìn)行數(shù)據(jù)分析，直接關(guān)系到下一步靜校正融合的可靠性［14］。變差函數(shù)參數(shù)的求取原則［15］為：點(diǎn)對(duì)的個(gè)數(shù)隨著滯后距離的增加而減少，估計(jì)的精度正比于數(shù)據(jù)對(duì)的個(gè)數(shù)，但用于計(jì)算變差函數(shù)的滯后距離的數(shù)目一般應(yīng)大于30個(gè)點(diǎn)對(duì)。塊金是在距離為0時(shí)的模型值，是測(cè)量不確定性的標(biāo)準(zhǔn)，表現(xiàn)為在很短的距離內(nèi)有較大的空間變異性，若研究目標(biāo)為區(qū)域上靜校正的變化情況，那么小的塊金效應(yīng)常數(shù)則反映所研究地區(qū)靜校正具有很好的連續(xù)性。

2.3變差函數(shù)求取

變差函數(shù)的模型主要有3種：指數(shù)模型、球狀模型及高斯模型［16-18］。由于區(qū)域靜校正相對(duì)穩(wěn)定，故本次研究選用隨機(jī)性和連續(xù)性均適中的球狀模型，因?yàn)樗梢钥焖龠_(dá)到基臺(tái)值。求取變差函數(shù)的關(guān)鍵是依次設(shè)置順物源方向的主方向、次方向和垂直方向的參數(shù)，從而得到地質(zhì)變量的空間分布規(guī)律，其基本步驟概括如下：①調(diào)節(jié)主變程，確定一個(gè)搜索方向，并計(jì)算這個(gè)搜索方向上的帶寬和搜索半徑等參數(shù)；②改變搜索半徑以及步長(zhǎng)數(shù)值，使得變差函數(shù)曲線與回歸曲線基本重合，此時(shí)也可以適當(dāng)?shù)馗淖兯阉鞣较蛑钡蕉咧睾?；③通過(guò)擬合選擇樣點(diǎn)間相關(guān)性最好的方向作為主方向，使次方向與主方向重合，便可得到主方向和主變程。以此類推，可以依次得出次方向和垂直方向的變程值。

3　基于變差函數(shù)的高精度靜校正擬合重構(gòu)技術(shù)

地震資料處理中的靜校正量既要確保構(gòu)造形態(tài)的可靠，又要保證準(zhǔn)確成像，因此使得滿足全區(qū)地震資料處理要求的靜校正量的求取難度加大。在地震資料處理過(guò)程中，靜校正可以分離出低頻分量和高頻分量2個(gè)部分［19］，其中，低頻靜校正量主要影響地震剖面的地下地質(zhì)構(gòu)造形態(tài)，高頻靜校正量主要影響反射同相軸的成像質(zhì)量。本文方法以準(zhǔn)確的非線性層析成像靜校正為基礎(chǔ)，先進(jìn)行靜校正的高、低頻分離，再對(duì)局部信噪比低的成像部分，按照成像的質(zhì)量，選取有成像優(yōu)勢(shì)的其他靜校正方法的高頻靜校正，如折射靜校正和層析靜校正等，并在所選2種高頻靜校正的邊界用變差函數(shù)進(jìn)行擬合，進(jìn)而與分離的低頻靜校正重構(gòu)，得到一套新的靜校正量。具體實(shí)現(xiàn)過(guò)程包括以下步驟：

（1）搜集地震采集工區(qū)的生產(chǎn)炮初至波資料。

（2）以準(zhǔn)確拾取的初至波時(shí)間為基礎(chǔ)，采用優(yōu)化的參數(shù)進(jìn)行層析反演靜校正量的計(jì)算。

（3）利用初至拾取的折射時(shí)間計(jì)算折射靜校正量，與步驟（2）計(jì)算得到的層析反演靜校正量（基礎(chǔ)靜校正量）進(jìn)行疊加成像效果對(duì)比，選出折射靜校正量成像效果好的區(qū)域范圍；或者將利用近地表信息獲取的其他靜校正量，與步驟（2）計(jì)算得到的基礎(chǔ)靜校正量進(jìn)行疊加成像效果對(duì)比，選出其他靜校正量成像效果好的區(qū)域范圍。

（4）根據(jù)野外采集施工的設(shè)計(jì)排列長(zhǎng)度及地下構(gòu)造的形態(tài)，確定低頻平滑半徑，在炮域和檢波點(diǎn)域?qū)游鲮o校正量以及折射或其他靜校正量分別分解為低頻分量和高頻分量。

分別用STA，STC和STR表示重構(gòu)后的靜校正量、層析靜校正量、折射或其他靜校正量，根據(jù)給定的排列長(zhǎng)度將層析靜校正量、折射或其他靜校正量分解為各自的高頻分量和低頻分量，即

（5）利用定量信噪比分析方法，在疊加數(shù)據(jù)體上選出折射或其他靜校正量成像效果好于層析靜校正量的區(qū)域。

（6）在步驟（5）選定的區(qū)域，利用變差函數(shù)對(duì)炮點(diǎn)靜校正高頻分量與檢波點(diǎn)靜校正高頻分量分別在炮點(diǎn)域和檢波點(diǎn)域進(jìn)行擬合。

（7）分別在炮點(diǎn)域和檢波點(diǎn)域?qū)游鲮o校正量的低頻分量和經(jīng)步驟（6）變差函數(shù)擬合得到的高頻分量進(jìn)行重構(gòu)，那么，經(jīng)過(guò)高、低頻分離與重構(gòu)后的靜校正量則為

4　應(yīng)用效果

吐哈盆地玉北構(gòu)造帶西起勝金口油田，東至魯克沁構(gòu)造帶，南起吐玉克油田，北至勝北生油洼陷，區(qū)域面積為300 km2。工區(qū)地形地貌為山地和戈壁，地表高程起伏大，其中覆蓋火焰山的面積占60%以上，地震地質(zhì)條件較為復(fù)雜，加之受火焰山逆沖推覆等復(fù)雜地質(zhì)條件的影響，給該區(qū)地震資料的處理造成了一定困難。從原始預(yù)疊加剖面也可以看出，原始未做靜校正的疊加剖面其靜校正問(wèn)題非常嚴(yán)重，特別是在山體主體構(gòu)造部位及推覆體構(gòu)造部位，基本識(shí)別不出構(gòu)造的有效形態(tài)（圖2）。

圖2　單炮靜校正問(wèn)題分析Fig.2　Statics analysis of shot

基于變差函數(shù)的高精度靜校正融合技術(shù)依據(jù)不同地表?xiàng)l件下靜校正方法的優(yōu)劣，結(jié)合地質(zhì)條件，合理建立符合研究區(qū)地表?xiàng)l件的近地表表層結(jié)構(gòu)模型，再根據(jù)模型的特征選擇填充速度及基準(zhǔn)面，并最終獲取精確的基準(zhǔn)面靜校正量。在該區(qū)將層析靜校正作為基礎(chǔ)靜校正方法進(jìn)行CMP面的高、低頻分離，獲得低頻分量（因?yàn)閷游鲮o校正的低頻部分較為準(zhǔn)確、可靠，同時(shí)在整個(gè)區(qū)域內(nèi)是經(jīng)過(guò)閉合的），再利用折射靜校正等其他靜校正方法所獲得的綜合靜校正結(jié)果進(jìn)行CMP面的高、低頻分離，獲得高頻分量。從圖3可以看出，層析靜校正量與折射靜校正量的整體趨勢(shì)非常接近，但在山體高陡構(gòu)造部位層析靜校正的效果明顯好于折射靜校正［圖3（b）紅色標(biāo)出區(qū)域］，而在地勢(shì)平坦且存在穩(wěn)定折射層的區(qū)域，折射靜校正又有較明顯的優(yōu)勢(shì)［圖3（c）藍(lán)色標(biāo)出區(qū)域］。為了獲取更合理、更高精度的靜校正結(jié)果，對(duì)二者的優(yōu)勢(shì)部位進(jìn)行了融合，具體融合過(guò)程如圖4所示。從圖3（d）可以看出，融合后的靜校正量在工區(qū)北部成像效果較好，在工區(qū)南部對(duì)火焰山構(gòu)造部位也刻畫(huà)得很清晰，整體疊加效果有所突破。

圖3　靜校正擬合重構(gòu)前、后疊加剖面Fig.3　Stack sections before and after static correction fitting reconstruction

圖4　靜校正擬合重構(gòu)過(guò)程Fig.4　The process of statics fitting reconstruction

5　結(jié)論

（1）針對(duì)復(fù)雜地表?xiàng)l件地區(qū)應(yīng)用單一的靜校正方法不能有效處理好整個(gè)工區(qū)的靜校正問(wèn)題，通過(guò)分析層析靜校正和折射靜校正的適用條件與優(yōu)缺點(diǎn)，并在此基礎(chǔ)上提出了基于變差函數(shù)擬合重構(gòu)的高精度三維靜校正方法，該方法可有效提高山地、戈壁、沙漠及黃土塬等復(fù)雜地表?xiàng)l件地區(qū)地震資料的成像質(zhì)量。

（2）在區(qū)域地質(zhì)認(rèn)識(shí)的基礎(chǔ)上，經(jīng)變差函數(shù)區(qū)域性約束所重構(gòu)的近地表表層結(jié)構(gòu)模型能最終獲取精確的基準(zhǔn)面靜校正量，并能更加合理地實(shí)現(xiàn)不同靜校正量的擬合和重構(gòu)，將其應(yīng)用在吐哈盆地火焰山復(fù)雜地表工區(qū)很好地解決了靜校正問(wèn)題，有效提高了該區(qū)三維地震資料的成像質(zhì)量。

（References）：

［1］吳杰，蘇勤，王建華.層析靜校正技術(shù)在柴北緣地區(qū)的應(yīng)用［J］.巖性油氣藏，2008，20（3）：79-82.

Wu Jie，Su Qin，Wang Jianhua.Application of tomographic static correction in north margin of Qaidam Basin［J］.Lithologic Reservoirs，2008，20（3）：79-82.

［2］熊翥.復(fù)雜地區(qū)地震資料處理思路［M］.北京：石油工業(yè)出版社，2002：31-41.

Xiong Zhu.Thoughts of seismic data processing in complex areas［M］.Beijing：Petroleum Industry Press，2002：31-41.

［3］車(chē)建英，張偉.初至折射波靜校正技術(shù)在復(fù)雜探區(qū)地震數(shù)據(jù)處理中的應(yīng)用［J］.中國(guó)煤田地質(zhì)，2006，18（1）：54-56.

Che Jianying，Zhang Wei.Application of preliminary refracted wave static corrction technique in complicated area seismic data processing［J］.Coal Geology of China，2006，18（1）：54-56.

［4］胡自多，賀振華，王宇超，等．綜合靜校正技術(shù)及其在川東高陡構(gòu)造區(qū)應(yīng)用效果［J］.天然氣地球科學(xué)，2010，21（6）：1046-1053.

Hu Ziduo，He Zhenhua，Wang Yuchao，et al.Comprehensive static correction and its application for high-steep structure in east SichuanBasin［J］.NaturalGasGeoscience，2010，21（6）：1046-1053.

［5］王孝，賀振華，王述江，等．多信息約束層析反演靜校正技術(shù)及其應(yīng)用［J］.天然氣地球科學(xué)，2010，21（2）：316-320.

Wang Xiao，He Zhenhua，Wang Shujiang，et al.Technology of static correction for multi-information constrained tomographic inversion and its appication［J］.Natural Gas Geoscience，2010，21（2）：316-320.

［6］馮澤元，李培明，唐海忠，等．利用層析反演技術(shù)解決山地復(fù)雜區(qū)靜校正問(wèn)題［J］.石油物探，2005，44（3）：284-287.

Feng Zeyuan，Li Peiming，Tang Haizhong，et al.Using tomographic inversion techniques to solve statics problem of complex mountainousareas［J］.GeophysicalProspectingFor Petrolum，2005，44（3）：284-287.

［7］韓曉麗，楊長(zhǎng)春，麻三懷，等．復(fù)雜山區(qū)初至波層析反演靜校正［J］.地球物理學(xué)進(jìn)展，2008，23（2）：475-483.

Han Xiaoli，Yang Changchun，Ma Sanhuai，et al.Static of tomographic inversion by first breaks in complex areas［J］.Progress in Geophysics，2008，23（2）：475-483.

［8］劉全新，李道善.關(guān)于山地靜校正和偏移基準(zhǔn)面的一些認(rèn)識(shí)［J］.天然氣地球科學(xué)，2003，14（1）：57-59.

Liu Quanxin，Li Daoshan.Points about land static correction and migration datum［J］.Natural Gas Geoscience，2003，14（1）：57-59.

［9］林伯香，孫晶梅，劉清林.層析成像低速帶速度反演和靜校正方法［J］.石油物探，2002，41（2）：136-140.

LinBoxiang，SunJingmei，Liu Qinglin.Tomography for LVZ velocity inversion and statics［J］.Geophysical Prospecting for Petroleum，2002，41（2）：136-140.

［10］何光明，賀振華，黃德濟(jì)，等.非線性層析靜校正在川西地區(qū)資料處理中的應(yīng)用［J］.石油物探，2006，45（1）：88-92.

He Guangming，He Zhenhua，Huang Deji，et al.The application of nonlinear tomography statics in Chuanxi data processing［J］. Geo-physical Prospecting for Petroleum，2006，45（1）：88-92.

［11］王西文，雍學(xué)善，王宇超，等．面對(duì)重點(diǎn)勘探領(lǐng)域的地震技術(shù)研究和應(yīng)用實(shí)效［J］.巖性油氣藏，2010，22（3）：83-90.

Wang Xiwen，Yong Xueshan，Wang Yuchao，et al.Study and application of seismic technologies for key exploration field［J］.Lithologic Reservoirs，2010，22（3）：83-90.

［12］李霞，張印堂.折射層析反演靜校正在伊朗卡山地區(qū)的應(yīng)用［J］.勘探地球物理進(jìn)展，2006，29（1）：38-41.

Li Xia，Zhang Yintang.The application of refraction and tomographicinversionstaticsinKashan area of Iran［J］.Progress in Exploration Geophysics，2006，29（1）：38-41.

［13］侯景儒，尹鎮(zhèn)南，李維明，等.實(shí)用地質(zhì)統(tǒng)計(jì)學(xué)［M］.北京：地質(zhì)出版社，1998：36-45.

Hou Jingru，Yin Zhennan，Li Weiming，et al.Practical geostatistics［M］.Beijing：Geology Press，1998：36-45.

［14］王家華，劉倩.儲(chǔ)層建模中對(duì)變差函數(shù)分析的幾點(diǎn)認(rèn)識(shí)［J］.石油化工應(yīng)用，2011，30（10）：5-7.

Wang Jiahua，Liu Qian.Variogram analysis in reservoir modeling［J］.Petrochemical Industry Application，2011，30（10）：5-7.

［15］周游，程時(shí)清，張敏.儲(chǔ)層建模中變差函數(shù)參數(shù)的設(shè)置［J］.西安石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，25（5）：25-27.

Zhou You，Cheng Shiqing，Zhang Min.Parameter setting of variogram in reservoir modeling［J］.Jounal of Xi'an Shiyou University：Natural Science Edition，2010，25（5）：25-27.

［16］劉國(guó)棟，晏寧平，黃文科，等．精細(xì)變差函數(shù)分析在靖邊氣田下古氣藏的應(yīng)用［J］.石油化工應(yīng)用，2009，28（7）：31-32.

Liu Guodong，Yan Ningping，Huang Wenke，et al.The application of meticulous variogram in low Palaeozoic of Jingbian gasfield［J］. Petrochemical Industry Application，2009，28（7）：31-32.

［17］段天向，劉曉梅，張亞軍，等．Petrel建模中的幾點(diǎn)認(rèn)識(shí)［J］.巖性油氣藏，2007，19（2）：102-107.

Duan Tianxiang，Liu Xiaomei，Zhang Yajun，et al.Discussion on geologic modeling with Petrel［J］.Lithologic Reservoirs，2007，19（2）：102-107.

［18］孫立雙，王恩德，劉國(guó)華，等．變差函數(shù)穩(wěn)健性研究［J］.東北大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，29（2）：277-279.

Sun Lishuang，Wang Ende，Liu Guohua，et al.Study on robust variogram［J］.Journal of Northeastern University：Natural Science，2008，29（2）：277-279.

［19］酆少英，何瑞珠，王巍，等．復(fù)雜地表?xiàng)l件下的靜校正方法研究［J］.河南石油，2002，16（3）：18-20.

Feng Shaoying，He Ruizhu，Wang Wei，et al.Study of statics method in complicated surface area［J］.Henan Petroleum，2002，16（3）：18-20.

（本文編輯：于惠宇）

Research and application of high-precision merged static correction technology based on variogram

ZHANG Yongfeng1，WANG Peng2，ZHANG Yabing3，XU Xiaoping3，SU Qin2，XU Xingrong2
（1.Xibei Engineering Corporation Limited，PowerChina，Xi'an 710065，China；2.PetroChina Research Institute of Exploration and Development-Northwest，Lanzhou 730020，China；3.Reserch Institute of Exploration and Development，PetroChina Jilin Oilfield Company，Songyuan 131200，Jilin，China）

The travel-time curve of seismic refraction waves is distorted by complex surface geological structures and cannot in-phase stack，so we need to carry out static correction of seismic data.Each static correction method has its own application condition：the elevation statics method cannot eliminate the effect of low velocity；the model statics method is limited by the length of alignment and the depth of tracking，so it is difficult to fully reflect the form of lowvelocity zone bottom boundary；the refraction statics method is only applied for areas with flat surface，uniform lateral velocity and significant refractive layers；the tomography statics method is adapted to any surface model inversion，but the inversion results are instability.A certain kind of static correction cannot solve all problems for increasingly complex exploration areas.This paper presented a high-precision three-dimensional static correction method based on variogram fitting reconstruction which is applicable to complex region with a variety of surface types，it can achieve fitting and reconstruction of different statics，solve the problem of different static correction advantage fitting inthe complex surface conditions，and effectively improve the image quality of three-dimensional seismic data to ensure the reliability of structural morphology.

high-precision static correction；variogram；fittingand reconstruction；complexsurface areas

P631.4

A

1673-8926（2015）03-0108-07

2014-12-08；

2015-02-03

國(guó)家重大科技專項(xiàng)“十二五”計(jì)劃課題“天然氣復(fù)雜儲(chǔ)層預(yù)測(cè)與烴類檢測(cè)地球物理技術(shù)研究及應(yīng)用”（編號(hào)：2011ZX05007-006）資助

張永峰（1982-），男，工程師，主要從事物探資料處理研究工作。地址：（710065）陜西省西安市中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司。E-mail：9924176@qq.com。