張 昭，陳中山，楊 彪

（1.河北省煤田地質(zhì)局 物測(cè)地質(zhì)隊(duì)，河北 邢臺(tái) 054000；2.河北省煤田地質(zhì)局，河北 石家莊 050085）

0 引 言

在內(nèi)蒙古塔然高勒礦區(qū)油房壕煤礦先期三維地震勘探工作的基礎(chǔ)上，對(duì)比以往成果報(bào)告，18 口鉆孔成果數(shù)據(jù)一致，12 口鉆孔成果數(shù)據(jù)不一致，成果報(bào)告中部分成果圖件對(duì)該區(qū)地層志丹群層位認(rèn)識(shí)有偏差。經(jīng)過(guò)數(shù)據(jù)篩選，以相同數(shù)據(jù)，利用測(cè)井約束波阻抗反演方法對(duì)地震數(shù)據(jù)進(jìn)行反演工作，對(duì)志丹群層位進(jìn)行了重新厘定。

1 概 況

白堊系下統(tǒng)志丹群（K1zh） 在勘探區(qū)西部和北部沖溝內(nèi)局部出露，其上部為深紅色泥巖、砂質(zhì)泥巖夾細(xì)粒砂巖，下部以灰綠色，淺紅色礫巖為主。志丹群與下伏地層侏羅系中統(tǒng)（J2） 不整合接觸，接觸面為紫紅色、雜色砂質(zhì)泥巖、泥巖與灰綠色、黃綠色粉砂巖[1-2]。該地層與圍巖波阻抗差異大，可形成反射波（TK1zh），該反射波能量強(qiáng)，連續(xù)性好，可全區(qū)追蹤。

2 資料整理

對(duì)勘探區(qū)30 口鉆孔（圖1） 資料進(jìn)行了統(tǒng)計(jì)整理。其中鉆孔柱狀數(shù)據(jù)由該區(qū)鉆孔柱狀成果圖統(tǒng)計(jì)而來(lái)，地質(zhì)剖面數(shù)據(jù)由以往成果報(bào)告中的勘探線(xiàn)剖面圖統(tǒng)計(jì)而來(lái)，地震反演解釋數(shù)據(jù)為工作站提取解釋層位數(shù)據(jù)結(jié)合鉆孔資料進(jìn)行反演工作的成果。依據(jù)統(tǒng)計(jì)數(shù)據(jù)，3 種統(tǒng)計(jì)方法鉆孔成果數(shù)據(jù)一致為18 口，不一致的為12 口。以往成果報(bào)告中勘探線(xiàn)剖面圖和鉆孔柱狀中志丹群地層層位存在偏差。

圖1 鉆孔分布Fig.1 Borehole distribution

3 常規(guī)地震解釋方法

地震解釋前選取鉆孔柱狀統(tǒng)計(jì)和地質(zhì)剖面線(xiàn)統(tǒng)計(jì)一致的鉆孔作為基準(zhǔn)孔。依據(jù)基準(zhǔn)鉆孔的測(cè)井曲線(xiàn)進(jìn)行人工合成記錄對(duì)志丹群反射波的地質(zhì)層位進(jìn)行標(biāo)定。通過(guò)聯(lián)井剖面的反射波進(jìn)行層位追蹤、網(wǎng)格對(duì)反射波層位進(jìn)行對(duì)比追蹤以及利用工作站的自動(dòng)追蹤拾取功能，由粗網(wǎng)格40 m×40 m 到細(xì)網(wǎng)格10 m×10 m，局部加密至5 m×5 m 網(wǎng)格解釋。利用三維數(shù)據(jù)體任意切取剖面，即切取縱向剖面、橫向剖面和水平切片。根據(jù)聯(lián)井剖面及任意方向剖面的特點(diǎn)，綜合利用各種資料，進(jìn)行檢查，由粗到細(xì)使解釋逐步深化，最終確定志丹群賦存形態(tài)。

4 反演工作

4.1 反演原理

地震勘探長(zhǎng)期以來(lái)是利用巖層的聲學(xué)特征確定巖性的分界面，這樣造成地震資料與煤田地質(zhì)的結(jié)合發(fā)生困難。為了使地震資料能與鉆井資料直接連接對(duì)比，需要將界面型的反射剖面轉(zhuǎn)換成巖層型的測(cè)井剖面，地震資料變成可與鉆井直接對(duì)比的形式。實(shí)現(xiàn)這種轉(zhuǎn)換的處理過(guò)程就是地震反演技術(shù)[3]。

反演方法主要分為疊前反演和疊后反演（圖2）。目前在生產(chǎn)中用的較多的是疊后振幅反演。較 為 流 行 的 軟 件 有 Jason、 STRATA、 ISIS、VanGurad 等。該方法以一維褶積模型為理論基礎(chǔ)，僅使用指定的地質(zhì)界面的振幅信息，而不是其反射時(shí)間。基本反演方法包括有遞歸反演、稀疏脈沖反演和基于模型的反演等。

圖2 地震反演方法流程Fig.2 Fiowchart of seismic inversion method

模型反演是基于地震正演模型技術(shù)發(fā)展而來(lái)的一種反演技術(shù)。首先要使用測(cè)井和地震解釋成果建立對(duì)礦藏地質(zhì)模型的初始猜測(cè)，即初始模型（用速度、密度表示）。利用正演算法合成地震剖面，即地質(zhì)模型[4-6]。其次將地震模型與實(shí)際地震剖面對(duì)比，反復(fù)修正初始模型，實(shí)現(xiàn)地震模型最優(yōu)化。即反演的最終結(jié)果為震模型與地震資料對(duì)比誤差足夠小且能夠滿(mǎn)足期望的邊界條件時(shí)的礦藏地質(zhì)模型。

此次反演工作采用的是測(cè)井約束反演方法。使用地震解釋層位的控制將測(cè)井聲阻抗曲線(xiàn)沿層位橫向逐道外推。與此同時(shí)，按層序的厚度變化對(duì)聲阻抗曲線(xiàn)進(jìn)行拉伸和壓縮[4]。多井約束時(shí)，井間則按井的距離加權(quán)內(nèi)插。初始模型此時(shí)保留了測(cè)井曲線(xiàn)的高分辨率。但井的外推和線(xiàn)性?xún)?nèi)插仍不能解決礦藏的橫向非均質(zhì)性問(wèn)題，需要地震資料來(lái)修正井外推和線(xiàn)性?xún)?nèi)插的聲阻抗曲線(xiàn)，實(shí)現(xiàn)非線(xiàn)性外推和內(nèi)插。具體做法為在給定的約束條件下不斷地在初始模型與地震道之間求解一個(gè)反射系數(shù)校正值|Δr(j)|，并用此來(lái)修正初始模型的阻抗曲線(xiàn)，制作合成地震道，使其盡可能地與地震道溫和[5-8]。最終修正的模型即為反演得到的波阻抗。

4.2 反演實(shí)現(xiàn)過(guò)程

測(cè)井約束波阻抗反演是一個(gè)系統(tǒng)化的處理、解釋過(guò)程，每個(gè)步驟都會(huì)影響最終的結(jié)果。多井約束波阻抗反演的處理過(guò)程中，測(cè)井資料的標(biāo)準(zhǔn)化處理、層位的精細(xì)標(biāo)定、地震地質(zhì)模型的建立和反演方法的選擇決定著反演的質(zhì)量和精度[9-10]。

4.2.1 資料準(zhǔn)備

收集研究區(qū)內(nèi)鉆孔柱狀信息、測(cè)井曲線(xiàn)（聲速、密度、自然伽馬、視電阻率等），以及主要目的層的鉆井分層數(shù)據(jù)，加載高保真地震數(shù)據(jù)體。

4.2.2 測(cè)井資料標(biāo)準(zhǔn)化處理

測(cè)井資料是波阻抗反演的最重要基礎(chǔ)資料之一。反演的關(guān)鍵之一是做好測(cè)井資料的標(biāo)準(zhǔn)化和歸一化處理[9-10]。

4.2.3 地震數(shù)據(jù)的應(yīng)用

測(cè)井約束波阻抗反演要求在地震資料處理過(guò)程中始終保持“三高”的處理流程，從而保證地震煤層反演結(jié)果的真實(shí)性。其中，“高保真”是獲得精確地下波阻抗資料的基礎(chǔ)，尤其對(duì)測(cè)井約束波阻抗反演這種基于地震數(shù)據(jù)的反演方法。地震資料分辨率在很大程度上決定了反演結(jié)果的真實(shí)分辨率，其原因就在于利用褶積模型實(shí)現(xiàn)的誤差判斷原則，可以使多個(gè)模型假設(shè)滿(mǎn)足誤差要求，在地震有限頻帶之外，往往不能實(shí)現(xiàn)通過(guò)不斷修改模型參數(shù)達(dá)到正確反演波阻抗的目的。

三維地震反演技術(shù)可以將地震資料的常規(guī)界面反射型剖面反演與測(cè)井資料直接對(duì)比的巖層型波阻抗物性剖面進(jìn)行結(jié)合，最終將地質(zhì)構(gòu)造模型、地震數(shù)據(jù)、井點(diǎn)測(cè)井資料有機(jī)地結(jié)合為一體。

5 應(yīng)用分析

綜合分析研究區(qū)地質(zhì)資料，以鉆孔綜合成果為基礎(chǔ)，通過(guò)常規(guī)地震解釋與反演工作，應(yīng)用地震綜合解釋成果和地震反演成果，對(duì)研究區(qū)志丹群地層層位進(jìn)行了重新厘定，重新掌握了該地層賦存形態(tài)的變化，并詳細(xì)分析了各鉆孔修正工作前變化情況（圖3 ～ 圖4）。

由修正后的成果可以看出，志丹群地層底界面的賦存形態(tài)埋藏深度較穩(wěn)定，起伏變化較修正前緩和平穩(wěn)，更接近真實(shí)的地質(zhì)形態(tài)，該成果已得礦方認(rèn)可，指導(dǎo)了區(qū)域地層新的認(rèn)識(shí)。

6 結(jié) 語(yǔ)

經(jīng)認(rèn)真分析研究該區(qū)地質(zhì)、鉆探、測(cè)井資料、地震資料和反演工作成果，并結(jié)合以往成果報(bào)告對(duì)志丹群地層層位的劃定，以鉆探綜合成果為主，應(yīng)用地震資料與地震反演成果，對(duì)研究區(qū)志丹群底界面重新修正、認(rèn)識(shí)。同時(shí)通過(guò)對(duì)波阻抗反演技術(shù)的運(yùn)用，更多地了解了研究區(qū)的地質(zhì)屬性特征，為全面認(rèn)識(shí)和分析該區(qū)的地質(zhì)情況提供了方便。

圖3 修正前志丹群底板等高線(xiàn)圖Fig.3 Base contour map of Zhidan Group before modification

圖4 修正后志丹群底板等高線(xiàn)圖Fig.4 Modified contour map of Zhidan Group floor