測井在露天礦邊坡軟弱夾層識別上的研究

趙耀

（煤科集團沈陽研究院有限公司，遼寧 撫順 113122）

測井在露天礦邊坡軟弱夾層識別上的研究

趙耀

（煤科集團沈陽研究院有限公司，遼寧 撫順 113122）

針對軟弱夾層賦存規(guī)律以及地球物理特征，運用費歇爾判別分析方法結(jié)合測井?dāng)?shù)據(jù)對鉆孔控制地層進行軟弱層識別。將該方法運用于F露天煤礦，以密度、聲波、中子、伽馬測井?dāng)?shù)據(jù)作為基礎(chǔ)樣本建立判別函數(shù)，結(jié)果表明，采用該方法能夠有效的對鉆孔周圍地層軟弱層位的識別，為邊坡安全性評價提供依據(jù)。

露天礦邊坡；軟弱夾層；測井；費歇爾判別分析；自動識別

0 引 言

露天礦邊坡高達(dá)300～500 m，邊坡的變形、破壞極為普遍，邊坡失穩(wěn)已成為危害人民生命財產(chǎn)安全及工程正常運營的災(zāi)害。露天礦邊坡通常為巖體邊坡，巖體內(nèi)的結(jié)構(gòu)面，尤其是軟弱結(jié)構(gòu)面的存在，常常是巖坡不穩(wěn)定的主要因素。因此在進行邊坡穩(wěn)定分析時，應(yīng)當(dāng)特別注意結(jié)構(gòu)面和軟弱層的影響。

目前我國在露天礦邊坡工程地質(zhì)勘察中，盡管沒有統(tǒng)一的軟弱層定義，但軟弱層具有強度低、天然含水量大、密度小、變形量大的特點是不可否認(rèn)的。主要由粘土頁巖、軟質(zhì)凝灰?guī)r、泥灰?guī)r、滑石片巖以及含有巖鹽或石膏成分的構(gòu)成。這類巖層遇水浸泡后易軟化，大大降低了巖體強度[1-2]。

通常在鉆探過程中，巖層中的薄層、軟弱結(jié)構(gòu)面難以發(fā)現(xiàn)。地球物理測井技術(shù)能夠探測到巖體內(nèi)破碎、斷裂帶以及軟弱夾層，并且可以準(zhǔn)確定位[3-5]，給邊坡失穩(wěn)預(yù)防提供了可靠的技術(shù)支持。

1 測井方法選用及軟弱層識別

地球物理測井常用方法有普通電阻率測井、感應(yīng)測井、自然伽馬測井、密度測井、中子孔隙度測井、井徑測井以及近年發(fā)展流行的成像測井和聲波陣列測井等。在實際工作中，根據(jù)探測的目的和要求可以選擇幾種方法進行組合，互相對比，驗證。

1.1 邊坡軟弱層測井響應(yīng)特征及原因

1）側(cè)向測井。側(cè)向測井又稱聚焦電阻率測井，是具有一個測量電極和多個屏蔽電極的一種測井方法。受地層厚度、相鄰層影響小，是劃分不同地層的有效方法[6]。

與普通電阻率測井相比，三側(cè)向測井更能反映原狀地層電性特征：由于屏流作用，使得視電阻率曲線受井眼影響小，主電極尺寸小，縱向分辨率高，有利于薄層劃分；其曲線具有對地層中心上、下對稱的特點，在巖層上、下界面處，曲線陡度大，劃分地層分界面更為精確；軟弱夾層由于常處與含水狀態(tài)，導(dǎo)致電阻率降低，結(jié)合深、淺側(cè)向重合方法能夠較好的識別軟弱層位[7]。

2）聲速測井。聲波測井是利用聲波在巖石等介質(zhì)中傳播時，幅度的衰減、速度或頻率的變化等聲學(xué)特性來研究鉆井地質(zhì)剖面、判斷固井質(zhì)量等問題的測井方法。目前邊坡工程地質(zhì)中常用聲速測井方法，又稱聲波時差測井。

聲速測井通過測量地層聲波速度，用來提供地層孔隙度資料，工程中常采用縱波聲速測井，通過測試縱波波速來判斷層位。當(dāng)探測地層有破碎帶或隱伏斷裂時，聲波時差曲線常出現(xiàn)較大波動，即周波跳躍現(xiàn)象。因此聲波測井在識別邊坡巖體中斷裂破碎弱層具有十分重要的作用[8]。

3）自然伽馬測井。自然伽瑪測井測量地層的天然放射性，測量自然伽瑪射線的強度。在沉積巖地層中，除了在地層中夾有火山灰、花崗巖質(zhì)沖積物等放射性物質(zhì)，或地層水中溶有鉀鹽外，通常純地層的放射性是很微弱的，而放射性元素主要聚集在粘土（泥質(zhì)）中。因此，自然伽瑪測井一般反應(yīng)地層中的粘土（泥質(zhì)）含量。它與自然電位測井同屬于粘土指示測井系列。由于不同地層具有不同的自然放射性強度，因此可根據(jù)自然伽瑪測井曲線判斷巖性。泥巖層或含有放射性物質(zhì)的地層呈高自然伽馬特征；砂巖層、致密地層及純灰?guī)r地層呈低自然伽馬特征。

4）井徑測井。井徑測井是測量井筒直徑大小的一種機械測井方法。在能夠利用測井中所測得的多條井徑曲線，計算出最大井徑、最小井徑及平均井徑。利用記錄的多條井徑數(shù)據(jù)，把井徑數(shù)據(jù)按井周順序插值計算，求出其井周的均勻的井徑數(shù)據(jù)列同時計算出每一深度點的最大值、最小值及平均值。

5）密度測井。密度測井屬于核測井的一種，通過探測地層反射的伽馬射線能量，經(jīng)標(biāo)定換算給出不同地層的密度值。通常情況下，斷層、破碎帶以及軟弱泥化夾層以及內(nèi)部含水的因素，常會導(dǎo)致密度曲線的變化。

6）中子孔隙度測井。中子測井利用中子源所放射出具有一定能量的中子與井周圍巖石起作用,從而實現(xiàn)發(fā)射中子的彈性散射和俘獲輻射等核反應(yīng)的測井方法。地層中的含氫量決定中子在地層中的減速過程，而含氫量與飽含水或油的孔隙體積相關(guān)?？梢岳脺y井所獲得的中子、熱中子參數(shù)來反映巖層天然含水量（濕度）的大小。

綜上所述，通過單一的測井曲線來判別弱層必然存在較大誤差，同時考慮若層機械強度低、天然含水量大、密度較小且孔隙度略大的特征，需綜合多種測井方法來識別。

其中：24 h連續(xù)運行設(shè)備的日平均耗電量為三鑒/微波探測器日平均耗電量、光口交換機日平均耗電量、網(wǎng)絡(luò)繼電器日平均耗電量三者之和。

1.2 判別分析方法識別軟弱層

根據(jù)鉆孔中軟弱層位測井?dāng)?shù)據(jù)的明顯差異，將異常數(shù)據(jù)作為樣本，經(jīng)過處理建立軟弱層識別標(biāo)準(zhǔn)，進而對井段各個層位進行歸類識別。主要步驟包括測井曲線選擇、建立判別函數(shù)確定判定標(biāo)準(zhǔn)、歸類識別。

1）測井?dāng)?shù)據(jù)特征分析。軟弱泥化夾層在物理性質(zhì)、化學(xué)性質(zhì)、放射性等方面具有明顯差異，測井曲線會有相應(yīng)的異常響應(yīng)，總體來說呈現(xiàn)高伽馬、高中子、低聲波和低密度的曲線特征，為多參數(shù)判別軟弱層提供了基礎(chǔ)條件，因此在選取曲線分析樣本數(shù)據(jù)時常會使用以上幾種曲線[10]。

2）費歇爾判別分析識別。費歇爾判別法是根據(jù)方差分析建立起來的一種區(qū)分各個總體的線性判別法，其基本思想就是投影，將n組k維數(shù)據(jù)投影到某一個方向上，使其投影的組和組之間盡可能的分開。具體步驟（圖1）首先確定訓(xùn)練樣本集：

其中，第一類、二類分別為：

定義各類類內(nèi)的離散度矩陣，即各類的協(xié)方差矩陣：

式中mj表示各類內(nèi)均值。

類內(nèi)離散度、總類內(nèi)離散度和類間離散度分別為：

為使得向需求的方向投影后兩類能夠盡可能的分離，而類內(nèi)部又盡可能的聚集，表示成如下準(zhǔn)則，及Fisher準(zhǔn)則：

將公式帶入并通過拉格朗日求極值方法得到所需投影方向和判別閾值：

最后將待確定樣本代入

通過判斷上式的符號和哪個類別相同，確定其屬于哪個類別[11-13]。

2 工程實例應(yīng)用

以F露天煤礦為例，該礦區(qū)地層由老到新有凝灰?guī)r、煤層、泥頁巖、砂巖。其中軟弱泥化夾層多發(fā)育在泥頁巖層中，于構(gòu)造作用下層間流動形成，遇水極易泥化。

測井采用重慶地質(zhì)儀器廠綜合數(shù)字測井系統(tǒng)JGS3，進行了深淺側(cè)向電阻率、自然伽馬、聲波、密度、中子、井徑共7項測井項目。在對原始數(shù)據(jù)的拼接校正和標(biāo)準(zhǔn)化基礎(chǔ)上，繪制了相應(yīng)的測井曲線并劃分了巖性段。

圖1 費歇爾判別程序圖

以ZL-1鉆井650～690 m井段為例（ZL-1井650～690 m曲線如圖2），根據(jù)綜合測井曲線，結(jié)合實際鉆孔資料進行對比煤層底板泥頁巖層位厚度約為1.5 m，具有明顯的高伽馬、高中子、低聲波和低密度的曲線特征。因此選擇GR、AC、DEN、CNL這4個最主要的影響因素作為判別因子，另選取5組不同深度軟弱層及多組其它巖層段數(shù)據(jù)為判別樣本，統(tǒng)計測井特征值，統(tǒng)計結(jié)果見表1。

采用MATLB計算軟件，編程實現(xiàn)測井樣本數(shù)據(jù)的分析，通過建立費歇爾判別分析的函數(shù)[14]，由此得到的投影方向多維向量方程和閾值分別為：

式中：Y為綜合判別值，X1為自然伽馬，X2為密度；X3為聲波；X4為中子。

對該鉆孔192.0～249.5 m井段285個數(shù)據(jù)進行了測試，對比鉆孔巖心，成功判別了該孔位內(nèi)3處軟弱層位，判別結(jié)果見表2及成果如圖3。

圖2 ZL-1井650～690 m曲線

表1 測井?dāng)?shù)據(jù)判別分析樣本

表2 費歇爾判別結(jié)果

通過建立判別函數(shù)對其余多個鉆孔識別了巖層中的軟弱夾層，其深度、厚度情況與鉆孔、錄井情況吻合度較高，說明了該方法的可行性。但該判別方法只能在軟弱層和其它巖性進行識別，無法進一步識別巖性，主要受測井項目數(shù)量的限制。

圖3 測井自動判別成果圖

3 結(jié) 論

利用測井方法結(jié)合費歇爾判別，對F露天煤礦其它施工的鉆孔進行驗證，均取得了較好的地質(zhì)效果。在同一鉆孔內(nèi)，通過此方法能夠快速識別軟弱夾層，確定潛在滑動面位置，對于邊坡穩(wěn)定的預(yù)防和評價解釋有著一定的輔助作用。今后應(yīng)在測井項目組合上、資料處理的程序上再做一些研究工作，相信測井技術(shù)在邊坡工程地質(zhì)勘察中的應(yīng)用前景將會更加樂觀。由于地質(zhì)條件差別大，導(dǎo)致判別函數(shù)不甚適用，需要另外建立判別函數(shù)。

［1］ 謝宗球，肖宏斌.邊坡軟弱層評價中全波列聲波測井的理論與實驗研究［J］.中南林業(yè)科技大學(xué)學(xué)報，2010，30（8）：109-111.

［2］ 丁立明，才慶祥.軟弱夾層對露天礦邊坡穩(wěn)定性的影響［J］．金屬礦山，2012，430（4）：40-43.

［3］ 桑黎明.地球物理測井在公路工程勘察中的應(yīng)用[J].黑龍江科技信息，2014，（25）：14-15.

［4］ 王好龍.基于測井?dāng)?shù)據(jù)判別分析法識別構(gòu)造煤的應(yīng)用［J］.煤炭技術(shù)，2014，33（12）：317-319.

［5］ 劉大野，陳立云.露天煤礦邊坡工程地質(zhì)勘察中幾個常見問題的芻議——以勝利煤田東二號露天煤礦、白音華煤田一號露天煤礦為例［J］.中國煤炭地質(zhì)，2010，430（4）：52-55.

［6］ 尉中良.地球物理測井［M］.北京：地質(zhì)出版社，2005.

［7］ 何公民.三側(cè)向測井在煤田地質(zhì)中的應(yīng)用［J］.煤炭技術(shù)，2009，28（6）：145-146.

［8］ 傅唯一，羅上強.聲波測井在勘察滑坡中的應(yīng)用［J］.勘查科學(xué)技術(shù)，2010（1）：52-54.

［9］ 丁國輝.煤系地層巖石物理力學(xué)性質(zhì)與測井曲線參數(shù)的相關(guān)性研究［D］.合肥:安徽理工大學(xué)，2010（6）.

［10］ 建林青.淺談測井方法在工程勘察中的應(yīng)用［J］.中國煤炭地質(zhì)，2008，20（12）：40-41.

［11］ 王鵬，王志章.核Fisher判別分析在火山巖巖性識別中的應(yīng)用［J］.測井技術(shù)，2015，39（3）：115-117.

［12］ 李偉偉，楊成.基于Fisher判別分析法的流體識別［J］.石油化工應(yīng)用，2013，32（12）：13-18.

［13］ 楊建，楊靜宇，王化鋒，等.Fisher線性鑒別分析的理論研究及其應(yīng)用［J］.自動化學(xué)報，2003，29（4）：482-492.

［14］ 謝中華.MATLAB統(tǒng)計分析與應(yīng)用：40個案例分析［M］.北京：北京航空航天大學(xué)出版社，2010.

【責(zé)任編輯：解連江】

Research of well logging for weak layer identification in Open-pit Mine slope

ZHAO Yao
(China Coal Technology and Engineering Group Shenyang Research Institute,Fushun 113122,China)

According to the law of occurrence and geophysical characteristics of weak interlayer,using Fischer discriminant analysis method,the article recognizes drilling control strata combining with logging data.This method is applied for open-pit mine,using density,sonic,neutron,gamma logging data as a basis sample discriminant function.The results show that this method can effectively identify weak layer of surrounding formation,which provide evidence for safety evaluation of open-pit slope

open-pit mine slope;weak layer;well logging;fisher discriminant-analysis;automatic recognition

TD824.7

B

1671－9816（2017）06－0052－04

10.13235/j.cnki.ltcm.2017.06.015

趙耀.測井在露天礦邊坡軟弱夾層識別上的研究［J］.露天采礦技術(shù)，2017，32（6）：52-55.

2017-03-17

國家自然科學(xué)基金面上項目（51274122）；國家自然科學(xué)基金煤炭聯(lián)合基金資助項目（U1361211）

趙 耀（1989—），男，黑龍江綏化人，助理工程師，碩士，2015年畢業(yè)于中國地質(zhì)大學(xué)（北京）地球探測與信息技術(shù)專業(yè)，主要從事露天礦邊坡穩(wěn)定性研究工作。