劉小洪，馮明友，楊午陽(yáng)，孫 輝，魏新建，劉 哲

（1.西南石油大學(xué)資源與環(huán)境學(xué)院；2.中國(guó)石油勘探開(kāi)發(fā)研究院西北分院；3.中國(guó)石油長(zhǎng)慶油田分公司勘探部）

利用Kohonen神經(jīng)網(wǎng)絡(luò)劃分二維地震相
——以柴達(dá)木盆地E區(qū)風(fēng)險(xiǎn)勘探為例

劉小洪1，馮明友1，楊午陽(yáng)2，孫 輝2，魏新建2，劉 哲3

（1.西南石油大學(xué)資源與環(huán)境學(xué)院；2.中國(guó)石油勘探開(kāi)發(fā)研究院西北分院；3.中國(guó)石油長(zhǎng)慶油田分公司勘探部）

為研究柴達(dá)木盆地E區(qū)大型背斜構(gòu)造沉積相及砂體展布規(guī)律，在二維地震數(shù)據(jù)閉合差校正、鄰區(qū)井標(biāo)定引入及精細(xì)層位解釋基礎(chǔ)上，利用改進(jìn)算法的Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)開(kāi)展二維地震相劃分研究，識(shí)別出三角洲前緣水下分流河道、分流間灣及灘壩等微相。本文研究認(rèn)為，研究區(qū)古流向?yàn)槟蠔|—北西向，儲(chǔ)集砂體較發(fā)育，主要富集于研究區(qū)中部，現(xiàn)今構(gòu)造東高點(diǎn)位于有利沉積相帶。改進(jìn)算法的Kohonen神經(jīng)網(wǎng)絡(luò)二維地震相劃分技術(shù)補(bǔ)充了沉積相研究成果，適合于西部二維地震資料覆蓋的風(fēng)險(xiǎn)探區(qū)，具較強(qiáng)的推廣價(jià)值。

Kohonen；神經(jīng)網(wǎng)絡(luò)；地震相；沉積相；柴達(dá)木盆地

0 引言

地震相指特定沉積相或地質(zhì)體在地震振幅、相位、連續(xù)性及反射特征組合中的綜合響應(yīng)。地震相分析即以實(shí)際地質(zhì)資料為基礎(chǔ)，利用地震參數(shù)并結(jié)合其它資料進(jìn)行綜合沉積環(huán)境解釋的過(guò)程。

地震相研究的核心為神經(jīng)網(wǎng)絡(luò)技術(shù)［1-4］。神經(jīng)網(wǎng)絡(luò)技術(shù)主要分為有監(jiān)督的神經(jīng)網(wǎng)絡(luò)技術(shù)（應(yīng)用于勘探中后期及開(kāi)發(fā)期）和無(wú)監(jiān)督的神經(jīng)網(wǎng)絡(luò)技術(shù)（多應(yīng)用于勘探前期），多名學(xué)者［5-10］曾利用神經(jīng)網(wǎng)絡(luò)技術(shù)開(kāi)展地球物理方面的研究。

研究區(qū)位于柴達(dá)木盆地北部風(fēng)險(xiǎn)勘探區(qū)，該研究區(qū)二維地震資料豐富，有監(jiān)督的神經(jīng)網(wǎng)絡(luò)技術(shù)受限。本文開(kāi)發(fā)出二維地震相分析及成圖軟件，利用Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)對(duì)地震道波形進(jìn)行分類(lèi)并將此分類(lèi)形成離散的“地震相”；再根據(jù)“擬合度”準(zhǔn)則對(duì)實(shí)際地震道進(jìn)行對(duì)比和分類(lèi)，細(xì)致刻畫(huà)出地震信號(hào)橫向變化規(guī)律及地震波形平面分類(lèi)特征；最后結(jié)合相關(guān)研究結(jié)果對(duì)地震資料進(jìn)行綜合地質(zhì)解釋?zhuān)M(jìn)而開(kāi)展儲(chǔ)層預(yù)測(cè)和含油氣性判別研究。

1 Kohonen神經(jīng)網(wǎng)絡(luò)原理

Kohonen神經(jīng)網(wǎng)絡(luò)是無(wú)監(jiān)督學(xué)習(xí)的人工神經(jīng)網(wǎng)絡(luò)，也稱(chēng)為自組織特征圖（Self-Organizing Map）。該網(wǎng)絡(luò)可實(shí)現(xiàn)對(duì)輸入模式特征進(jìn)行拓?fù)溥壿嬘成?，是一種高維輸入向量向低維輸出向量投影的方法，主要由輸入層、隱藏層（Kohonen層）和輸出層組成（圖 1）［11］。其原理是：模擬生物神經(jīng)元之間的競(jìng)爭(zhēng)作用進(jìn)行信息處理，通過(guò)比較輸入樣本間的空間最小距離來(lái)不斷調(diào)整神經(jīng)元間的權(quán)值，從而使相似的輸入樣本（即優(yōu)勝神經(jīng)元）聚集為一種分類(lèi)，經(jīng)過(guò)每一類(lèi)反復(fù)迭代，最終輸出一組已分類(lèi)的數(shù)據(jù)。

圖1 Kohonen神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)示意圖Fig.1 Structural diagram of Kohonen neural network

2 實(shí)現(xiàn)步驟

利用Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)劃分二維地震相的主要步驟為：①對(duì)二維地震資料進(jìn)行去噪處理和閉合差校正，并結(jié)合研究區(qū)或相鄰區(qū)域的測(cè)井資料進(jìn)行精細(xì)層位標(biāo)定，在標(biāo)定結(jié)果上解釋目的層位；②設(shè)置地震相分類(lèi)數(shù)，對(duì)權(quán)值、學(xué)習(xí)速率及鄰域進(jìn)行初始化，以地震道值作為輸入樣本；③采用公式（式中 wmn表示輸入樣本到輸出結(jié)點(diǎn)的權(quán)值；t表示時(shí)刻，為正整數(shù)）計(jì)算輸入樣本xm到每一個(gè)輸出結(jié)點(diǎn)n之間的距離dn，此距離被稱(chēng)為Euclidean距離，并選取具有最小距離的結(jié)點(diǎn)（n為輸出結(jié)點(diǎn)）；④修改n與鄰近結(jié)點(diǎn)之間的權(quán)值，采用公式 wmn（t+1）=wmn（t）+η（t）［xm（t）-wmn（t）］（式中 η（t）為一個(gè)增益項(xiàng)，0＜η（t）＜1，隨時(shí)間增加下降至0）修正權(quán)值；⑤若增益項(xiàng)的值減小至0，則結(jié)果滿意，輸出結(jié)果，否則轉(zhuǎn)換至第二步重新訓(xùn)練，直至神經(jīng)網(wǎng)絡(luò)收斂為止；⑥對(duì)輸出的地震相結(jié)果進(jìn)行分析，將地震相轉(zhuǎn)化為沉積相。

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

3.1 研究區(qū)概況

研究區(qū)構(gòu)造位于柴達(dá)木盆地北緣斷塊帶鄂博梁—葫蘆山—伊克雅烏汝背斜帶的中段，北鄰昆特依—伊北凹陷，南接一里坪凹陷，呈“兩凹夾一隆”的構(gòu)造格局。由圖2可看出，研究區(qū)現(xiàn)今地表構(gòu)造為北西—南東向延伸的長(zhǎng)軸背斜，背斜長(zhǎng)約46 km，寬約8 km，面積約368 km2。由周緣鉆探可知，該區(qū)主要目的層為上新統(tǒng)上油砂山組（）、上新統(tǒng)下油砂山組（）、中新統(tǒng)上干柴溝組（N1）3個(gè)層位。 多名學(xué)者結(jié)合頂部構(gòu)造特征（圖3）、鄰區(qū)構(gòu)造研究成果［12-14］、二維地震資料似“平點(diǎn)”含氣性預(yù)測(cè)結(jié)果［15］等認(rèn)為，研究區(qū)具大型氣藏潛力。但鑒于資料缺乏、沉積特征復(fù)雜及構(gòu)造含氣性不確定等問(wèn)題，至今該區(qū)的勘探仍未能突破。

圖2 研究區(qū)構(gòu)造位置圖Fig.2 Structural location of the study area

圖3 研究區(qū)頂部時(shí)間構(gòu)造圖Fig.3 Time-structure map for the top ofin the study area

研究區(qū)內(nèi)二維測(cè)線經(jīng)連片處理，測(cè)線交點(diǎn)處的時(shí)差、振幅、相位、頻率基本閉合，二維地震測(cè)網(wǎng)為2 km×3 km～2 km×4 km，分布均勻，可滿足平面研究要求。所以，本文利用改進(jìn)算法的Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)對(duì)目的層開(kāi)展二維地震相劃分研究，并結(jié)合沉積相平面研究成果進(jìn)行綜合分析，以期為沉積相及砂體展布規(guī)律研究提供技術(shù)支持。

3.2 Kohonen神經(jīng)網(wǎng)絡(luò)地震相劃分

本文利用Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)進(jìn)行二維地震相劃分的流程主要由二維地震資料預(yù)處理、目的層層位精細(xì)解釋、地震反射特征分類(lèi)、地震波形Kohonen神經(jīng)網(wǎng)絡(luò)分析、地震相—沉積相平面對(duì)比分析等組成。首先，對(duì)研究區(qū)34條二維地震測(cè)線進(jìn)行相位校正和波形、振幅及頻率一致性處理，并通過(guò)閉合差校正來(lái)解決二維地震資料因施工年度、施工因素及處理流程不統(tǒng)一等引起的問(wèn)題；其次，結(jié)合相鄰區(qū)域葫蘆山構(gòu)造井-震標(biāo)定結(jié)果，將標(biāo)定層位引入研究區(qū)，并對(duì)研究區(qū)目的層進(jìn)行精細(xì)層位解釋?zhuān)蛔詈蟾鶕?jù)目的層地震反射同相軸的振幅、頻率、連續(xù)性等特征，將地震剖面分為4種類(lèi)型（對(duì)應(yīng)平面分布如圖3中Ⅰ—Ⅳ所示）：Ⅰ型，同相軸反射呈雜亂型，橫向極不連續(xù)，中等振幅，中等頻率，為混合間互沉積（圖4a）；Ⅱ型，振幅為雙強(qiáng)谷夾弱峰型，中等頻率，同相軸橫向連續(xù)性中等，為三角洲前緣分流間灣沉積（圖4b）；Ⅲ型，振幅為強(qiáng)谷強(qiáng)峰型，頻率較高，同相軸橫向連續(xù)性好，為灘壩砂沉積（圖4c）；Ⅳ型，振幅為強(qiáng)谷強(qiáng)峰型，頻率較低，同相軸橫向連續(xù)性較差，同相軸略顯下拉現(xiàn)象，可能為三角洲前緣分流河道沉積（圖4d）。其它地震反射特征介于4種類(lèi)型之間，難以直接甄別。

圖4 研究區(qū)典型地震反射特征剖面Fig.4 Typical seismic sections of reflection characteristics in the study area

在地震剖面分類(lèi)研究基礎(chǔ)上，沿層提取地震振幅，將振幅波形中不同道的每個(gè)采樣點(diǎn)作為輸入數(shù)據(jù)，送入神經(jīng)網(wǎng)絡(luò)進(jìn)行學(xué)習(xí)，反復(fù)迭代，修改權(quán)值，直至網(wǎng)絡(luò)訓(xùn)練成熟，網(wǎng)絡(luò)拓?fù)溱呌诜€(wěn)定，再重新將所有提取的振幅波形送入網(wǎng)絡(luò)，從而得到最終的地震相分類(lèi)結(jié)果。在計(jì)算時(shí)為避免受地震資料品質(zhì)的影響，本文選擇大時(shí)窗進(jìn)行屬性提取，巖層上、下分別提取16 ms?？紤]到地震相分類(lèi)數(shù)過(guò)少可能會(huì)導(dǎo)致結(jié)果失真，首先選擇分類(lèi)數(shù)為8進(jìn)行分析（圖5a），為使結(jié)果突出優(yōu)勢(shì)相且與地質(zhì)意義相吻合，進(jìn)一步將分類(lèi)數(shù)合并為4（中間采用過(guò)渡色充填）。因此，Kohonen神經(jīng)網(wǎng)絡(luò)地震相最終分類(lèi)結(jié)果（圖5b）能整體反映沉積體系及砂體展布特征。

圖5 研究區(qū)頂部Kohonen分類(lèi)剖面（a）、地震相分類(lèi)結(jié)果解釋?zhuān)╞）及沉積相平面圖（c）Fig.5 Seismic section of Kohonen classification （a）,interpretation result of seismic facies classification （b） and sedimentary facies ichnography （c） of the top ofin the study area

區(qū)域沉積相研究成果表明：漸新世早期由于盆地邊緣斷裂活動(dòng)強(qiáng)烈，盆地整體沉降進(jìn)入大型坳陷期，研究區(qū)接受了一套三角洲相沉積，周邊主要發(fā)育濱淺湖相沉積；中新世的沉積環(huán)境基本繼承了漸新世沉積格局，研究區(qū)繼續(xù)接受三角洲相和濱淺湖相沉積，水體范圍有所擴(kuò)大，但水體淺，為廣泛的弱氧化環(huán)境；晚第三紀(jì)上新世早期，研究區(qū)水體逐漸變深，沉積物顏色變暗，物源豐富，沉積厚度較大。各期有利沉積相帶均為分選較好、儲(chǔ)滲性佳的濱淺湖灘壩及三角洲前緣水下分流河道微相（圖5c）。

3.3 應(yīng)用效果分析

基于Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)劃分的地震相平面圖細(xì)化了沉積相研究成果。本文研究認(rèn)為：研究區(qū)內(nèi)主要發(fā)育三角洲前緣水下分流河道及濱淺湖灘壩等微相，古流向?yàn)槟蠔|—北西向，砂體廣泛分布，尤其中部地區(qū)（今構(gòu)造高點(diǎn)）相對(duì)富集；有利區(qū)域（紅—黃色）為三角洲前緣水下分流河道微相（圖3b），現(xiàn)今構(gòu)造東高點(diǎn)位于有利相帶內(nèi)。

實(shí)鉆表明，構(gòu)造東高點(diǎn)在N21氣測(cè)顯示活躍、儲(chǔ)層發(fā)育、測(cè)井解釋儲(chǔ)層厚，表明預(yù)測(cè)結(jié)果可信。

4 結(jié)論及認(rèn)識(shí)

Kohonen神經(jīng)網(wǎng)絡(luò)技術(shù)在淺層氣預(yù)測(cè)、油氣評(píng)價(jià)、巖性識(shí)別、地震相劃分等方面具有優(yōu)勢(shì)，在運(yùn)用時(shí)需注意以下幾點(diǎn)：

（1）高信噪比的地震資料、規(guī)則的測(cè)網(wǎng)、準(zhǔn)確的層位解釋是進(jìn)行二維地震相劃分的前提，因此二維地震資料預(yù)處理和標(biāo)定極為重要。

（2）通常選擇較大時(shí)窗進(jìn)行劃分，分類(lèi)時(shí)可對(duì)整個(gè)數(shù)據(jù)進(jìn)行快速掃描以選取最佳分類(lèi)數(shù)。

（3）神經(jīng)元數(shù)量應(yīng)適中，過(guò)多會(huì)影響分類(lèi)速度，過(guò)少則導(dǎo)致分類(lèi)結(jié)果重復(fù)，影響預(yù)測(cè)精度。

（4）神經(jīng)網(wǎng)絡(luò)運(yùn)用前可利用二維地震剖面反射特征進(jìn)行分析，劃分時(shí)需在地質(zhì)思路的指導(dǎo)下結(jié)合其它資料綜合對(duì)比，才能全面、準(zhǔn)確地得出符合實(shí)際的劃分結(jié)果。

［1］Ediriwickrema J，Khorram S.Hierachical Maximum Likelihood Classification for Improved Accuracies［J］.IEEE Transon Geoscience and Remote Sensing，1997，35（4）：810-816.

［2］EssenreiterR，KarrenbachM，TreitelS.Identification and classification of multiple reflections with self-organizing maps［J］.Geophysical Prospecting，2001，49（3）：341-352.

［3］Baan M V D，Jutten C.Neural networks in geophysical applications［J］.Geophysics，2000，65（4）：1032-1047.

［4］Taner M T，Berge，Walls J D，et al.Well log calibration of kohonenclassified seismic attributes using bayesian logic ［J］.Journal of Petroleum Geology，2001，24（4）：405-416.

［5］鄧傳偉，李莉華，金銀姬，等.波形分類(lèi)技術(shù)在儲(chǔ)層沉積微相預(yù)測(cè)中的應(yīng)用［J］.石油物探，2008，47（3）：262-265.

［6］孫勤華，劉曉梅，劉建新，等.利用波形分析技術(shù)半定量預(yù)測(cè)塔中碳酸鹽巖儲(chǔ)層［J］.巖性油氣藏，2010，22（1）：101-103.

［7］王宏斌，張虎權(quán)，衛(wèi)平生，等.碳酸鹽巖地震儲(chǔ)層學(xué)在塔中地區(qū)生物礁灘復(fù)合體油氣勘探中的應(yīng)用［J］.巖性油氣藏，2010，22（2）：18-23.

［8］侯斌，桂志先，許輝群，等.應(yīng)用多屬性神經(jīng)網(wǎng)絡(luò)方法預(yù)測(cè)油氣［J］.巖性油氣藏，2010，22（3）：118-120.

［9］衛(wèi)平生，潘建國(guó)，張虎權(quán)，等.地震儲(chǔ)層學(xué)的概念、研究方法和關(guān)鍵技術(shù)［J］.巖性油氣藏，2010，22（2）：1-5.

［10］楊杰，衛(wèi)平生，李相博.石油地震地質(zhì)學(xué)的基本概念、內(nèi)容和研究方法［J］.巖性油氣藏，2010，22（1）：1-6.

［11］Kohonen T.The Self-organizing maps［J］.Proceedings of the IEEE，1990，78（9）：1464-1480.

［12］宋成鵬，張曉寶，汪立群，等.柴達(dá)木盆地北緣天然氣成因類(lèi)型及氣源判識(shí)［J］.石油與天然氣地質(zhì)，2009，31（1）：90-96.

［13］羅群，黃捍東，李玉喜.柴達(dá)木盆地北緣淺層滑脫斷裂下盤(pán)油氣富集規(guī)律及控藏模式［J］.中國(guó)石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，51（3）：34-38.

［14］沈亞，陳元忠，孫維軍，等.柴達(dá)木盆地地震地質(zhì)研究新進(jìn)展［J］.石油地球物理勘探，2008，43（S1）：146-150.

［15］陳迎賓，張壽庭，萬(wàn)傳治，等.柴達(dá)木盆地鄂博梁Ⅲ號(hào)構(gòu)造地震低速異常識(shí)別及成因探討［J］.地球物理學(xué)進(jìn)展，2008，23（6）：1894-1902.

Application of Kohonen neural network to 2D seismic facies division in E area of Qaidam Basin

LIU Xiao-hong1， FENG Ming-you1， YANG Wu-yang2， SUN Hui2， WEI Xin-jian2， LIU Zhe3
（1.School of Resources and Environment， Southwest Petroleum University， Chengdu 610500， China； 2.Research Institute of Petroleum Exploration&Development-Northwest， PetroChina， Lanzhou 730020， China； 3.Department of Exploration，Changqing Oilfield Company， PetroChina， Xi’an 710018， China）

The situation ofE area is hard toproceed deep research and risk assessment bythe absent ofprospectingwell data.Based on seismic mis-tie calibration,adjacent well calibration and fine horizon interpretation,Kohonen neural network method is applied to carry out two-dimensional seismic facies classification of target zone.Microfacies of delta front such as distributary channel,interdistributary bay and sand bar are recognized.The paleocurrent direction is suspected from southeast to northwest.Reservoir sand bodies developed well in the middle of the study area,preliminary prospecting well is located in favorable sedimentary facies.Sedimentary facies division are supplied and refined bythe result ofseismic facies，which can supplysignificant foundation for geometric arrangement ofrisk wells and regional breakthrough.

Kohonen； neuralnetwork； seismicfacies； sedimentaryfacies； QaidamBasin

TE112.221

A

2010-12-20；

2011-02-16

劉小洪，1980年生，女，講師，主要從事儲(chǔ)層沉積學(xué)研究。地址：（610500）四川省成都市新都區(qū)新都大道 8號(hào)西南石油大學(xué)資源與環(huán)境學(xué)院。E-mail：liuxiaohong_swpu@163.com

1673-8926（2011）04-0115-04

郭言青）