葛利華，姜 弢，徐學(xué)純，賈海青，楊志超

1.吉林大學(xué)地球信息探測儀器教育部重點實驗室，長春 130026

2.吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院，長春 130026

3.吉林大學(xué)地球科學(xué)學(xué)院，長春 130061

0 引言

隨著地震勘探［1］工作的不斷深入，容易勘探的平原區(qū)域日益減少，主要勘探任務(wù)轉(zhuǎn)向了地表條件復(fù)雜的區(qū)域［2］。這些區(qū)域近地表不均勻性突出，特別是有表層高速層［3］、山地［4］或灘淺海［5］等區(qū)域，給地震的激發(fā)、接收和靜校正帶來較大的困難，表層結(jié)構(gòu)調(diào)查及速度建模是優(yōu)化井深設(shè)計及數(shù)據(jù)處理靜校正的有效方法之一。

近年來，國內(nèi)外很多地球物理學(xué)家針對不同地區(qū)的近地表進行了表層調(diào)查和速度建模研究。2000年，Chandola等［6］對印度Cauvery盆地砂巖出露區(qū)的近地表模型進行了研究，并指出了近地表模型對該區(qū)地震數(shù)據(jù)采集和處理的重要意義。2001年，Zhu等［7］提出首波層析成像建立近地表模型的方法，實驗數(shù)據(jù)和真實數(shù)據(jù)都表明該方法明顯改善了靜校正和深度成像。2003年，夏竹等［8］對中國西部復(fù)雜地表地區(qū)近地表特征與近地表結(jié)構(gòu)進行了研究，對山地、黃土源、沙漠和戈壁4種典型復(fù)雜地表地區(qū)的近地表地層的地質(zhì)、地球物理特性及對靜校正、速度分析、疊加成像等的影響進行了詳盡的分析。2008年，Enikanselu等［9］利用29個小折射數(shù)據(jù)和1口微測井數(shù)據(jù)，揭示了尼日利亞尼日爾三角洲Mono地區(qū)的風化層地球物理特征。2010年，孫維昭等［10］發(fā)展了基于微測井資料進行層析反演重建高分辮率近地表速度模型的方法，進行了相應(yīng)的高精度靜校正量的計算，最后在我國南方某地區(qū)進行了實際應(yīng)用。2011年，李衛(wèi)忠等［11］提出近地表相和近地表相分析的概念，利用可控震源地震采集、近地表調(diào)查等數(shù)據(jù)，進行了近地表黏彈性參數(shù)反演，并建立了高精度的極淺近地表層模型。2011年，Guevara等［12］用2口微測井數(shù)據(jù)在哥倫比亞某地區(qū)建立了橫波速度模型，得到隨深部高速變化的速度模型及相關(guān)的巖性特征，用于數(shù)據(jù)處理的靜校正。2012年，Yang等［13］提出基于 Monte-Carlo抽樣方法的近地表速度建模方法，用合成數(shù)據(jù)和實際數(shù)據(jù)驗證了該方法的有效性。2013年，Lu等［14］提出基于多模塊快速行進法（multi-stencils fast marching method，MSFM）的地震走時層析成像方法，用于重建復(fù)雜近地表速度模型，提高了建立速度模型的效率和準確性。

此次，為滿足國家深部探測項目對方法比對研究的需求，同時為遼西葫蘆島東部進行首次區(qū)域地震勘探井深設(shè)計提供依據(jù)，為后期數(shù)據(jù)處理提供準確的靜校正量，筆者基于地表相似系數(shù)的表層模型建立方法對遼西葫蘆島東部進行了表層速度建模。

1 遼西葫蘆島東部地表條件

遼西葫蘆島東部地區(qū)以山地為主，出露的地層為典型的華北型。該測區(qū)地表起伏較大，地勢總體呈西北高、東南低的趨勢，為松嶺山脈延續(xù)分布丘陵地帶。此次實驗測線長16.9km，為西北－東南走向，樁號西北小，東南大，海拔高程為4～191m；中部至南部海邊為農(nóng)田，土層厚度小于0.5m，海拔高程8～80m；交通發(fā)達、人口密度大、地面干擾源多；北部為山區(qū)，地表裸露巖石，起伏較大且通行困難。該測區(qū)發(fā)育地層比較齊全，包括太古宇［15］、元古宇［16－17］、古生界、中生界［18］及新生界。大部分地段出露地層為太古宇變質(zhì)巖或巖漿巖，局部小面積出露沉積巖。

2 表層結(jié)構(gòu)速度建模原理

2.1 微測井工作原理

微測井按照實現(xiàn)方式不同分為地面微測井［19］、井中微測井和雙井微測井［20］等。參照表層調(diào)查實驗點多種方法比對實驗結(jié)果，筆者采用井中微測井方法。井中微測井采用井中激發(fā)，地面接收的方式，得到透射波［21－22］記錄，根據(jù)透射波初至?xí)r間來研究近地表結(jié)構(gòu)。

井中微測井通過改變炸藥或雷管在井中的激發(fā)位置獲得不同的初至?xí)r間。激發(fā)時井中注水，確保與大地良好耦合［23］，以得到可靠、有效的微測井地震記錄。

得到微測井地震記錄后，把透射波的初至?xí)r間轉(zhuǎn)化為井口位置自激自收［24］的垂直時間T0，設(shè)i點的垂直時間為T0i，則

其中：ti為i點的初至?xí)r間；Hi為激發(fā)點i的深度；d為井檢距。

將轉(zhuǎn)換后的時間T0與對應(yīng)的深度置于時間－深度坐標系。不同深度的T0位于同一速度層內(nèi)時，點的分布為一直線，直線斜率的倒數(shù)即為該層的層速度。2條直線的交點為上下層介質(zhì)的分界面。

2.2 地表速度模型建立原理

微測井點的數(shù)量有限，得到的微測井解釋成果不足以代表整條測線的表層情況。為了獲得整條測線的低、降速層厚度和速度，需要建立近地表模型。近地表與界面、界面與界面之間存在近代沉積的連續(xù)性和繼承性，它們之間具有一定的相干關(guān)系。可以利用表層調(diào)查結(jié)果和層間的這種相干關(guān)系，建立整個工區(qū)表層模型［25］。

設(shè)A、B為2個微測井點，水平距離為D1，低速層的厚度和速度分別為ZA、ZB、vA、vB，如圖1所示。

圖1 表層模型建立示意圖Fig.1 Schematic of surface modeling

設(shè)層間關(guān)系系數(shù)為K，A、B之間地面上任意一點為G。由A、B2點低速層厚度直線內(nèi)插得到G點的厚度為ZAB，并把內(nèi)插直線與G點的豎線交點記為C，G與C點的高程差為H，則G點的低速層厚度ZG為

其中：EG為G點的地表高程；EC為C點的地表高程，由A、B2點的地表高程直線內(nèi)插得到。H＝EG－EC：當C點位于G點以下時，H為正值；當C點位于G點以上時，H為負值。

G點的速度vG為

其中，D2為A與G間的水平距離。

設(shè)N為一定范圍內(nèi)的控制點數(shù)，則層間關(guān)系系數(shù)K為

其中：E1，j為第一個界面第j點的高程；ˉE1為第一個界面一定范圍內(nèi)的平均高程；E2，j為第二個界面第j點的高程；ˉE2為第二個界面一定范圍內(nèi)的平均高程；Em，j為（E1，j－ˉE1）和（E2，j－ˉE2）中較大的一個。當K＝1時，高速層頂界面和地表高程面完全相似；當K＝0時，高速層頂界面和地表高程面無關(guān)，完全為搜索半徑內(nèi)控制點高速層頂高程線性平滑。

根據(jù)公式（2）、（3）、（4），可逐一計算測線內(nèi)其他物理點表層的層速度和厚度。

3 微測井采集方法設(shè)計及參數(shù)

3.1 微測井采集方法設(shè)計

該測區(qū)內(nèi)地表黃土層較薄，為保證每個層位不少于4個控制點，淺層激發(fā)點加密。該測區(qū)以山地為主，大多數(shù)井潛水面較深，為保證每次都在水中激發(fā)，激發(fā)過程中需要不斷補充水源。利用Google Earth、地形圖等工具，提前適當調(diào)整微測井點位，盡量選擇離道路近的井位進行微測井，保證供水便利。為防止炸藥在井中激發(fā)時對井壁造成破壞，激發(fā)改用雷管，激發(fā)前注水時往井中倒入適量膨潤土保護井壁。向井中注水過程中，容易沖掉井壁上的碎石，使碎石脫落，測井后測井電纜拔不出來造成生產(chǎn)井浪費，選擇測井點位時適當微調(diào)，選擇井況較好的井，且注水過程采取保護井壁措施，微測井電纜下井時，與井底保留一定距離，防止雷管激發(fā)后掉落物堆積在井底拔卡住微測井電纜。

3.2 微測井采集參數(shù)

經(jīng)過試驗，確定施工因素。采用 WZG-24A地震儀，采集時采用固定增益12dB，記錄長度0.5s，采樣率0.2ms。自上而下各激發(fā)點雷管數(shù)量均為1發(fā)，以雷管底部對齊爆破電纜的標記處為準。采取相應(yīng)的防浮措施，確保每一層位至少有4個控制點（個別低速層很薄的點不少于3個），井中激發(fā)點間距布設(shè)遵循上密下疏的原則，0.50～5.00m激發(fā)間距為0.25m，6.00～10.00m激發(fā)間距為0.50m，10.00m以下激發(fā)點間距為1.00m，具體實施時依據(jù)低降速層變化情況微調(diào)激發(fā)點距。采用型號為PS-38的單只檢波器，共12道，扇形擺放接收，為保證檢波器與井口高差，1－5道距井口1.00m，6－12道距井口2.00m，如圖2所示。

4 微測井解釋及速度建模結(jié)果

山地微測井工作存在困難，鄰區(qū)微測井工作開展較少。建昌遼河油田2008年和2012年的2次山地勘探，2008年生產(chǎn)物理點2 587個，只有6口微測井，2012年有所提高，5 506個物理點中也僅有46個微測井點，其他點位的表層工作用準確度和困難性都較小的小折射代替。

此次實驗，該區(qū)生產(chǎn)物理點234個，2013年1月到2013年4月完成井中微測井8口，平均每2 km 1個點。

圖2 井中微測井觀測系統(tǒng)示意圖Fig.2 Schematic of well-microlog observing system

4.1 微測井解釋結(jié)果

按照野外施工方法得到8口井的炮集記錄，用克朗軟件繪制時深曲線。以12304樁號為例，該樁號微測井井深18m：0.75～4.25m激發(fā)點間距為0.25m，4.5～9.5m 激發(fā)間距為0.5m，11～18m激發(fā)間距為1m，炮集記錄如圖3所示。根據(jù)得到的炮集記錄，得到微測井解釋成果（圖4）。由圖4可見，該微測井點近地表分為4層。其中：低速層的速度為798m/s，厚度為2.39m；降速層的速度為2 354m/s，厚度為4.61m；在低速層和降速層之間存在速度為5 112m/s、厚度為2.11m的較薄高速夾層；該微測井點處高速層1頂界面的埋深為9.11 m，速度為4 887m/s。

微測井的表層結(jié)構(gòu)計算結(jié)果如表1所示。由表1可見：高速層速度為2 878～4 887m/s，頂界面埋深為4.50～9.56m，平均為6.36m；低速層速度為349～993m/s，平均厚度為3.34m；工區(qū)東側(cè)12304號樁附近地下3m處有約2m的較薄高速層。同時，對留作井中微測井的8口生產(chǎn)井進行巖性錄井，結(jié)果為6口井處地表被較薄的黃土覆蓋，2口井花崗巖出露；高速層頂界面最厚處為7.50m左右，巖性以花崗巖為主，1口井為石英砂巖，2口井為凝灰?guī)r。

圖3 12304樁號井中微測井炮集記錄Fig.3 Shot gather of uphole in pile No.12304

圖4 12304樁號井中微測井時深曲線Fig.4 Time-depth curve of uphole in pile No.12304

表1 微測井解釋結(jié)果Table 1 Interpretation results of uphole

4.2 速度建模結(jié)果

根據(jù)式（4）和表1的微測井解釋成果，計算可得地表相似系數(shù)，如圖5所示。

圖5 地表相似系數(shù)Fig.5 Similarity coefficient of near-surface

由圖5可見：整個測區(qū)低速層的相似系數(shù)均較大，為0.9～1.0，說明該測區(qū)低速層的連續(xù)性較好；降速層和高速層1相似系數(shù)左側(cè)較小，中間較大，右側(cè)介于兩者之間，可見降速層和高速層1在測線西北端變化最劇烈，連續(xù)性也較差；中間區(qū)段連續(xù)性較好，東南段介于兩者之間。

由激發(fā)點坐標、高程，接收點坐標、高程，以及微測井解釋成果、相似系數(shù)，得到該工區(qū)厚度分布（圖6）和速度分布（圖7）。

由圖6可見：低速層厚度變化不大，為1～5m，平均厚度為2.9m；降速層厚度除個別區(qū)段厚度變化較大外，大部分區(qū)段位于10m之內(nèi)，平均厚度為8.7m；結(jié)合8口微測井結(jié)果，高速層頂界面平均厚度為9.7m。根據(jù)經(jīng)驗值，最佳激發(fā)井深在高速層頂面以下約4m為宜，本測區(qū)的藥量試驗確定激發(fā)藥量為8kg，藥柱長度為2m，考慮到部分區(qū)段高速層較深，設(shè)計激發(fā)井深13～15m。

由巖性錄井結(jié)果知，13～15m處激發(fā)巖性主要為致密的花崗巖。

圖6 表層厚度分布圖Fig.6 Distribution of surface thickness

圖7 表層速度分布圖Fig.7 Distribution of surface velocity

由圖7可見：低速層速度變化較小，在1 000 m/s內(nèi)；降速層速度測線兩端較中間低，測線兩端速度變化較平緩，總體為1 000～2 800m/s；高速層1速度變化大，在10500樁號附近較低，范圍為2 800～4 900m/s。

根據(jù)解釋成果及建立的模型，得到靜校正量數(shù)據(jù)，并繪制表層界面高程圖（圖8），為工區(qū)的靜校正提供了真實可靠的數(shù)據(jù)。

為消除地表高程差大對低、降速層變化的影響，更好地了解表層厚度和速度的變化，將圖8地表拉平，得到地表拉平后的表層剖面圖（圖9）。由圖9可見測線內(nèi)近地表的速度分布以及各速度層的厚度分布。在深度15m處，除10500樁號附近速度較低外，大部分區(qū)段速度高于2 800m/s，從而驗證了設(shè)計激發(fā)井深13～15m的合理性，激發(fā)進入高速層，由巖性錄井結(jié)果可知激發(fā)巖性主要為致密的花崗巖。

圖8 表層界面高程圖Fig.8 Elevation of surface interfaces

圖9 地表拉平后表層剖面圖Fig.9 Surface profiles after leveling the surface

5 結(jié)論

1）該工區(qū)地形起伏大，測線西北端表層速度和厚度變化大，中間段連續(xù)性較好，東南端變化較緩和。

2）該測區(qū)表層分2～3層：低速層厚度為1.68～4.33m，速度為349～993m/s；降速層厚度為1.82～5.23m，速度為989～24 481m/s；高速層1埋深為4.64～12.28m，速度為2 878～4 887m/s。

3）表層結(jié)構(gòu)建模不僅為靜校正提供了數(shù)據(jù)，同時也為激發(fā)井深選擇提供了依據(jù)。

（

）：

［1］林君.電磁驅(qū)動可控震源地震勘探原理及應(yīng)用［M］.北京：科學(xué)出版社，2004.

Lin Jun.Seismic Exploration and Application of Vibrator Driven by Electromagnetic［M］.Beijing：Science Press，2004.

［2］丁建榮，毛鳳鳴，郝天珧.表層火成巖覆蓋地區(qū)的表層結(jié)構(gòu)調(diào)查［J］.地球物理學(xué)進展，2007，22（3）：872-878.

Ding Jianrong，Mao Fengming，Hao Tianyao.Surface Structure Survey in the Area Where Covered with the Igneous Rock Near Surface ［J］. Progress in Geophysics，2007，22（3）：872-878.

［3］王建花，李慶忠，邱睿.淺層強反射界面的能量屏蔽作用［J］.石油地球物理勘探，2003，38（6）：589-596，602.

Wang Jianhua，Li Qingzhong， Qiu Rui.Energy Shielding Effect of Strong Reflection Interface in Near-Surface［J］.Oil Geophysical Prospecting，2003，38（6）：589-596，602.

［4］林依華，張中杰，尹成，等.復(fù)雜地形條件下靜校正的綜合尋優(yōu)［J］.地球物理學(xué)報，2003，46（1）：101-106.

Lin Yihua，Zhang Zhongjie，Yin Cheng，et al.Hybrid Optimization of Static Estimation in Complex Topography［J］.Chinese Journal of Geophysics，2003，46（1）：101-106.

［5］徐錦璽，邱燕，何京國，等.灘淺海地震勘探采集技術(shù)應(yīng)用［J］.地球物理學(xué)進展，2005，20（1）：66-70.

Xu Jinxi，Qiu Yan，He Jingguo，et al.Seismic Exploration Acquisition Technique Application of Beach Shallow Sea［J］.Progress in Geophysics，2005，20（1）：66-70.

［6］Chandola S K，Prakash A.Near-Surface Characteristics of a Riyalur-Pondicherry Sub-Basin and Their Relevance to Seismic Data Acquisition and Processing［C］//SEG Expanded Abstract of 70th SEG Ann International Meeting.Calgary：SEG，2000：1283-1286.

［7］Zhu Tianfei，Cheadle S，Petrella A，et al.First-Arrival Tomography for Near-Surface Model Building［C ］//EAGE 63rd Conference ＆ Technical Exhibition.Amsterdam：EAGE，2001：11-15.

［8］夏竹，張少華，王學(xué)軍.中國西部復(fù)雜地表地區(qū)的近地表特征與表層結(jié)構(gòu)探討［J］.石油地球物理勘探，2003，38（4）：414-424.

Xia Zhu，Zhang Shaohua，Wang Xuejun.Discussion on Near-Surface Characters and Structures in Complex Area of Western China ［J］. Oil Geophysical Prospecting，2003，38（4）：414-424.

［9］Enikanselu P A.Geophysical Seismic Refraction and Uphole Survey Analysis of Weathered Layer Characteristics in the“Mono”Field，North Western Niger Delta，Nigeria［J］.Pacific Journal of Science and Technology，2008，9（2）：537-545.

［10］孫維昭，谷躍民，徐剛，等.高分辨率近地表速度模型重建及在靜校正中應(yīng)用［J］.地球物理學(xué)進展，2010，25（5）：1757-1762.

Sun Weizhao， Gu Yuemin， Xu Gang，et al.Reconstruction of the High-Resolution Near-Surface Velocity Model and Its Application to the Static Correction Under Complex Surface Conditions［J］.Progress in Geophysics，2010，25（5）：1757-1762.

［11］李衛(wèi)忠，劉連升，范留明，等.基于近地表相分析的全局約束近地表建模方法［J］.地質(zhì)勘探，2011，31（1）：41-47.

Li Weizhong，Liu Liansheng，F(xiàn)an Liuming，et al.Global Constrained Near Surface Modeling Based on a Neap Surface Seismic Facies Analysis［J］.Geologic Prospecting，2011，31（1）：41-47.

［12］Guevara S，Margrave G，Agudelo W，et al.Near-Surface S-Wave Velocity Models from Two Uphole Surveys［C］//Twelfth International Congress of the Brazilian Geophysical Society.Rio de Janeiro：The Brazilian Geophysical Society，2011.

［13］Yang Kai，Zhou Xing，Li Hui，et al.Near-Surface Velocity Model Construction Based on a Monte-Carlo Scheme［J］.Applied Geophysics，2012，9（4）：475-482.

［14］盧回憶，劉伊克，常旭.基于 MSFM的復(fù)雜近地表模型走時計算［J］.地球物理學(xué)報，2013，56（9）：3100-3108.

Lu Huiyi，Liu Yike，Chang Xu. MSFM-Based Travel-Times Calculation in Complex Near-Surface Model［J］.Chinese Journal of Geophysics，2013，56（9）：3100-3108.

［15］王慶龍.遼西興城地區(qū)新太古代花崗質(zhì)雜巖特征與成因［D］.長春：吉林大學(xué)，2012.

Wang Qinglong.Characteristics and Genesis of Neoarchean Granitic Complex in Xingcheng Area of Western Liaoning［D］.Changchun：Jilin University，2012.

［16］楊明春，鄭常青，張梅生，等.遼寧興城地區(qū)中元古界長城系大紅峪組石英砂巖質(zhì)礫巖古沉積環(huán)境研究［J］.礦物巖石，2012，32（3）：103-109.

Yang Mingchun， Zheng Changqing， Zhang Meisheng， et al. Study on Paleo-Sedimentary Environment of Dahongyu Formation Sandstone Conglomerate of Changcheng System of Mesoproterozoic in Xingcheng Area，Liaoning Province［J］.Journal of Mineralogy and Peterology，2012，32（3）：103-109.

［17］王旖旎，張梅生，李曉波，等.遼寧興城夾山地區(qū)中元古界團山子組疊層石的發(fā)現(xiàn)及其地層學(xué)意義［J］.吉林大學(xué)學(xué)報：地球科學(xué)版，2012，42（增刊1）：336-343.

Wang Yini，Zhang Meisheng，Li Xiaobo，et al.Discovery of Mesoproterozoic Stromatolites from the Tuanshanzi Formation of Changcheng System and Its Stratigraphical Significance in Jiashan Area of Xingcheng，Liaoning Province［J］.Journal of Jilin University：Earth Science Edition，2012，42（Sup.1）：336-343.

［18］許敏，薛林福，王東坡.遼西地區(qū)中生代火山－沉積盆地群特征及成因機制［J］.世界地質(zhì)，1997，16（2）：35-39.

Xu Mei， Xue Linfu， Wang Dongpo. The Characteristics and Mechanism of Mesozoic Volcanic-Sedimentary Basin Group of Western Liaoning［J］.Global Geology，1997，16（2）：35-39.

［19］李明海.地面微測井在山地表層結(jié)構(gòu)調(diào)查中的應(yīng)用［J］.勘探地球物理進展，2008，31（5）：378-382.

Li Minghai.Application of Uphole Methods in the Investigation of Near Surface Structure［J］.Progress in Exploration Geophysics，2008，31（5）：378-382.

［20］李天樹，陳寶德，蘇德仁.雙井微測井技術(shù)在表層結(jié)構(gòu)調(diào)查中的應(yīng)用［J］.石油物探，2004，43（5）：471-474.

Li Tianshu，Chen Baode，Su Deren.Application of Twin-Well Microlog in Near Surface Investigation［J］.Geophysical Prospecting for Petroleum，2004，43（5）：471-474.

［21］唐小平，劉寬厚.透射波勘探在尋找金屬礦中的探討［J］.地球物理學(xué)進展，2012，27（4）：1671-1678.

Tang Xiaoping，Liu Kuanhou.Transmitted Wave Exploration Discussion on Search Metal Minerals［J］.Progress in Geophysics，2012，27（4）：1671-1678.

［22］孫章慶.起伏地表條件下的地震波走時與射線路徑計算［D］.長春：吉林大學(xué)，2011.

Sun Zhangqing.The Seismic Traveltimes and Raypath Computation Under Undulating Earth’s Surface Condition［D］.Changchun：Jilin University，2011.

［23］邸志欣，丁偉，王增明.復(fù)雜山前帶地震勘探采集技術(shù)的實踐與認識［J］.石油物探，2012，51（6）：548-561.

Di Zhixin，Ding Wei，Wang Zengming.Practice and Understanding of Seismic Acquisition Technology in Complicated Foothill Area ［J］. Geophysical Prospecting for Petroleum，2012，51（6）：548-561.

［24］陳猛，劉建紅，韓波，等.模型正演技術(shù)在庫車地震勘探中的應(yīng)用［J］.天然氣工業(yè)，2007，27（增刊1）：167-169.

Chen Meng， Liu Jianhong， Han Bo， et al.Application of Modeling Technology of Seismic Exploration in the Kuche［J］.Natural Gas Industry，2007，27（Sup.1）：167-169.

［25］熊定鈺.復(fù)雜地表靜校正方法研究及應(yīng)用［D］.成都：成都理工大學(xué)，2011.

Xiong Dingyu.Static Methods on Rugged Topography：Investigation and Application［D］.Chengdu：Chengdu University of Technology，2011.