劉盛陽 朱煜華

①中石化河南油田分公司油氣開發(fā)管理部 ②中石化河南油田分公司勘探開發(fā)研究院

常規(guī)成圖軟件在誤差趨勢(shì)面校正方面往往采用反距離加權(quán)法，本次研究重點(diǎn)采用克里金算法提高了誤差趨勢(shì)面校正精度。在連片拼圖方面，通過利用在工區(qū)重疊處設(shè)置偽井、確定最小時(shí)差線方法，提高連片成圖精度，使多工區(qū)連片構(gòu)造圖與已鉆井構(gòu)造誤差小于5 m，滿足了科研生產(chǎn)的需求。

TH地區(qū)由STM、YQD、YQZ等15個(gè)三維地震工區(qū)組成，面積6600 km2。因各個(gè)工區(qū)地震資料采集、處理參數(shù)不一致，尤其是地震資料處理過程中應(yīng)用的基準(zhǔn)面、充填速度不一致，導(dǎo)致地震相位和地震速度場(chǎng)產(chǎn)生時(shí)差，是影響構(gòu)造連片成圖的主要因素[1]。通過對(duì)15個(gè)工區(qū)的重疊部分的時(shí)差分析，不同工區(qū)間的時(shí)差較大，時(shí)差范圍在-120 ms～130 ms，給連片構(gòu)造成圖精度帶來較大難度。

1 方法原理

1.1 方法原理

一般常用的誤差校正趨勢(shì)面是反距離加權(quán)法，該插值方法效率較高，但是該算法的加權(quán)系數(shù)僅為距離的函數(shù)，忽略了采樣點(diǎn)處的空間關(guān)系；而克里金插值方法中，權(quán)重取決于樣本點(diǎn)、未知點(diǎn)的距離及空間關(guān)系的擬合模型，考慮了數(shù)據(jù)中存在的空間相關(guān)距離和方向偏差，其插值結(jié)果往往具有更強(qiáng)的規(guī)律性和合理性，對(duì)待插值點(diǎn)的誤差值進(jìn)行線性無偏最優(yōu)估計(jì)[2]，公式為：

Z（x0）表示未知樣點(diǎn)的值；Z（xi）表示未知樣點(diǎn)周圍的已知樣本點(diǎn)的值；N為已知樣本點(diǎn)的個(gè)數(shù)；λi為第i個(gè)樣本點(diǎn)的權(quán)重。

同時(shí)滿足無偏估計(jì)的條件：

與反距離插值算法相比，克里金算法具有以下優(yōu)點(diǎn)：①模型同時(shí)具有局部和全局的統(tǒng)計(jì)特性，可以分析已知信息的趨勢(shì)和動(dòng)態(tài)，能進(jìn)行外推插值；②克里金算法作為線性回歸分析的一種改進(jìn)技術(shù)，包含線性回歸部分和非參數(shù)部分；非參數(shù)部分被視作隨機(jī)過程的實(shí)現(xiàn)，可使插值具有較高的精度和保真度。

1.2 克里金誤差趨勢(shì)面校正

利用平均速度得到的初始深度構(gòu)造圖與已鉆井分層深度存在一定誤差，需利用鉆井?dāng)?shù)據(jù)校正。常用的校正方法是用井點(diǎn)處的誤差生成一個(gè)誤差校正趨勢(shì)面，利用趨勢(shì)面校正構(gòu)造圖，得到與鉆井深度吻合的構(gòu)造圖。

2 應(yīng)用實(shí)例

首先對(duì)TH地區(qū)各工區(qū)開展層位的精細(xì)解釋，完成15個(gè)工區(qū)的等T0圖；然后利用沿層平均速度進(jìn)行時(shí)深轉(zhuǎn)換，得到各工區(qū)的深度構(gòu)造圖；利用已知井分層深度數(shù)據(jù)，統(tǒng)計(jì)構(gòu)造圖與鉆井深度的誤差，采用克里金算法擬合誤差趨勢(shì)面，對(duì)初步構(gòu)造圖進(jìn)行誤差校正；在相鄰工區(qū)的重疊處設(shè)置偽井，確定最小時(shí)差拼接線，完成相鄰工區(qū)的拼接，最終完成連片成圖。

2.1 克里金誤差趨勢(shì)面校正

由于TH地區(qū)的鉆井分布不均勻，在應(yīng)用誤差趨勢(shì)面進(jìn)行初步構(gòu)造圖校正時(shí)，利用反距離插值算法擬合的誤差趨勢(shì)面在已知井點(diǎn)附近會(huì)出現(xiàn)“牛眼”現(xiàn)象（圖1），造成井校后構(gòu)造圖的局部形態(tài)發(fā)生畸變[3]。

圖1 STM工區(qū)反距離加權(quán)插值誤差趨勢(shì)面圖

采用克里金插值方法對(duì)鉆井分層誤差進(jìn)行插值，得到的誤差趨勢(shì)面如圖3所示。該算法產(chǎn)生的誤差趨勢(shì)面能較好地反映數(shù)據(jù)中蘊(yùn)含的趨勢(shì)性規(guī)律，誤差趨勢(shì)保持較好，沒有明顯的異?，F(xiàn)象（圖2）。

圖2 STM工區(qū)克里金插值誤差趨勢(shì)面圖

從STM工區(qū)克里金誤差趨勢(shì)面校正后的構(gòu)造圖對(duì)井誤差分析結(jié)果來看，對(duì)井誤差小于5 m（圖3），井校后的構(gòu)造圖精度較高。因此，采用克里金插值的算法完成各個(gè)工區(qū)的初步構(gòu)造圖的井校工作，使校正后的構(gòu)造圖具有較高精度。

圖3 STM工區(qū)T46z 克里金校正后對(duì)井誤差分布圖

2.2 最小時(shí)差連片成圖

在TH地區(qū)北部，由于鉆井?dāng)?shù)較少，不能覆蓋工區(qū)重疊處。根據(jù)地震資料品質(zhì)高低依次劃分為Ⅰ、Ⅱ、Ⅲ、Ⅳ共4個(gè)等級(jí)[4]，在工區(qū)重疊處適當(dāng)添加偽井參與校正[5]，控制相鄰工區(qū)的閉合差，提高構(gòu)造成圖精度。設(shè)計(jì)偽井遵循的原則是：加入偽井后，工區(qū)參與校正的所有井在平面上的分布能覆蓋工區(qū)的大部分區(qū)域和主要構(gòu)造單元；相鄰工區(qū)的邊界線附近，可以適當(dāng)增加偽井的密度，以使相鄰工區(qū)的構(gòu)造圖經(jīng)過鉆井校正后，構(gòu)造等值線能夠較好地實(shí)現(xiàn)拼接；因?yàn)檫吔绺浇牡卣鹳Y料品質(zhì)往往會(huì)變差，偽井參數(shù)難以準(zhǔn)確求取，偽井的位置需要與工區(qū)邊界保持一定的距離[6]。

在完成時(shí)深轉(zhuǎn)換及鉆井、偽井校正后，得到各個(gè)工區(qū)深度構(gòu)造圖，然后對(duì)相鄰工區(qū)的等值線進(jìn)行了精細(xì)的拼接。首先將某個(gè)單層地震反射層深度構(gòu)造圖等值線抽稀至20 m等值線間隔，然后按兩個(gè)工區(qū)等值線誤差最小處確定拼接線，刪除掉每個(gè)工區(qū)拼接線外側(cè)等值線，最后將兩個(gè)工區(qū)的等值線無縫拼接一起。采用同樣的方法完成各個(gè)工區(qū)的連片成圖，構(gòu)造成圖精度對(duì)井誤差達(dá)到5 m。

3 結(jié)論和認(rèn)識(shí)

（1）相比反距離加權(quán)算法，克里金插算法不僅考慮了數(shù)據(jù)點(diǎn)的空間相關(guān)距離，而且考慮了數(shù)據(jù)間的方向偏差，采用克里金插值法校正構(gòu)造圖的精度更高。

（2）在無鉆井控制的工區(qū)重疊處設(shè)置偽井校正，是減小工區(qū)間閉合差，提高連片成圖精度的有效手段。

（3）根據(jù)地震資料品質(zhì)確定各工區(qū)優(yōu)先級(jí)別，開展時(shí)差校正，校正過程中采用最小時(shí)差原則，是多工區(qū)連片構(gòu)造成圖成功的關(guān)鍵。