煤層陷落柱散射波數(shù)值模擬與成像

曹志勇，王 偉，王 赟

1 中國科學(xué)院地質(zhì)與地球物理研究所，100029 北京

2 資源與環(huán)境信息系統(tǒng)國家重點實驗室，中國科學(xué)院地理科學(xué)與資源研究所，100101 北京

3 礦床地球化學(xué)國家重點實驗室，中國科學(xué)院地球化學(xué)研究所，550002 貴陽

煤層陷落柱散射波數(shù)值模擬與成像

曹志勇1，王 偉2＊，王 赟3

1 中國科學(xué)院地質(zhì)與地球物理研究所，100029 北京

2 資源與環(huán)境信息系統(tǒng)國家重點實驗室，中國科學(xué)院地理科學(xué)與資源研究所，100101 北京

3 礦床地球化學(xué)國家重點實驗室，中國科學(xué)院地球化學(xué)研究所，550002 貴陽

煤層陷落柱是煤田勘探開發(fā)中常見的一種典型的非均勻地質(zhì)體.由于來自陷落柱的反射信號少、反射能量弱，使得基于反射波原理的常規(guī)地震成像方法難以有效識別陷落柱.本文以散射波理論為基礎(chǔ)，采用數(shù)值模擬方法，研究了陷落柱的散射波場特征，研究表明地面接收的波場中含有來自陷落柱陡傾角界面的散射波場.通過共散射點道集波場的模擬，可以清晰地識別散射波，獲得地下散射點和非均勻地質(zhì)體的信息，判斷散射點的位置，從而勾畫出不均勻地質(zhì)體的形態(tài).采用等效偏移距假設(shè)抽取共散射點道集，在此基礎(chǔ)上進行疊前偏移，對陷落柱成像；模擬與實際數(shù)據(jù)成像結(jié)果對比表明此方法能夠合理地提取散射點的散射波場信息，對陷落柱形態(tài)及內(nèi)部結(jié)構(gòu)準(zhǔn)確成像，是一種有效的煤田陷落柱成像方法.

陷落柱，散射波，共散射點道集，等效偏移距，成像

1 引 言

陷落柱是中國華北、華南等地區(qū)煤田開采中的一種典型地質(zhì)體，對煤炭與煤層氣資源的開采具有危害.陷落柱的形成主要與構(gòu)造破碎帶、裂隙、地下水的活動和巖溶等因素有關(guān)［1－2］.在含煤地層的下部由于大量地下水的常年活動，導(dǎo)致周圍可溶性巖石溶蝕形成空洞，在重力作用下，上部含煤巖層塌陷，形成陷落柱.通常陷落柱是由塊度大小不均、排列雜亂的上部地層塌陷物膠結(jié)、堆積而成.與圍巖相比，在地層的連續(xù)性、產(chǎn)狀、巖性等方面均有很大的差別，導(dǎo)致正常反射波組不連續(xù)或能量變?nèi)酰?－4］.在陷落柱模擬與識別中，一般基于射線理論對陷落柱斷陷點的繞射波進行模擬與識別.在斷陷點以外或者陷落柱面陡傾情況下，很難觀測到反射波［5－7］.盡管陷落柱與圍巖存在著很大的巖性變化，即具有明顯的波阻抗差異，但陷落柱一般為吸收介質(zhì)，經(jīng)過陷落柱界面反射的能量很弱.在實際生產(chǎn)中，由于陷落柱在空間分布上的復(fù)雜性和地震響應(yīng)的特殊性，在偏移成像剖面上，陷落柱的幾何邊界很難被清晰成像［8－10］.

目前，常規(guī)二維、三維反射地震勘探方法在煤田勘探上已經(jīng)取得了良好的效果，但從陷落柱勘探的已有成果來看，效果并不理想.陷落柱的解釋主要是在反射波地震資料的基礎(chǔ)上，利用相干技術(shù)，通過煤層的連續(xù)性或陷落柱的特征反射波識別陷落柱；煤層水平切片上解釋的陷落柱位置與實際位置一般都存在較大偏差，且通常解釋范圍比實際范圍要大［11－12］.煤田陷落柱的精確勘探不僅影響井下開采效率和資源的充分利用，同時對安全生產(chǎn)有著重要意義.

地震波散射理論的發(fā)展為研究煤田陷落柱的成像方法提供了新的思路［13－18］.陷落柱可以看作地下非均勻地質(zhì)體，其彈性波場主要表現(xiàn)為多種類型的波干涉疊加，信噪比較低，能量和頻率差異較大，難以形成有效的強反射；但可以觀測到不均勻體產(chǎn)生的廣義繞射波（或稱為散射波），因此可以考慮采用散射波方法來實現(xiàn)陷落柱成像.

本文在散射波理論的基礎(chǔ)上，通過陷落柱的散射波場正演，分析波場特征，然后通過共散射點抽道集，對陷落柱的波場進行等效偏移距成像，以說明陷落柱的散射成像的可行性與有效性.

本文采用Bancroft提出的等效偏移距方法（Equivalent Offset Method，EOM）形成共散射點道集，然后對共散射點（Common Scatter Point，CSP）道集進行積分求和，此方法是一種疊前Kirchhoff偏移方法.該方法在形成共散射點道集時，不僅將共中心點道集上的反射能量映射到散射雙曲線上，還充分利用了散射能量，從而形成高信噪比的共散射點道集；而且該方法在形成共散射點道集時，對觀測系統(tǒng)要求不高，更適用于山地等復(fù)雜地形勘探數(shù)據(jù)，應(yīng)用靈活；該方法在實際應(yīng)用中，除形成共散射點道集單獨實現(xiàn)，處理流程中的其他步驟均可以借用現(xiàn)有的商業(yè)軟件模塊實現(xiàn).

2 原理與方法

2.1 散射波基本原理

在空間—頻率域（r，ω），波場可以分為背景場和總場，并滿足標(biāo)量波動方程

其中ω為角頻率，量綱為LT－1，c0是背景速度，單位m·s－1，r＝（x，y，z）是觀測點坐標(biāo)，rs＝（xs，ys，zs）是震源位置坐標(biāo)，δ是delta函數(shù).觀測點在背景場和總場中坐標(biāo)一致.定義速度擾動為：α（r）＝1－／c2（r），為無量綱量，并令k＝ω／c0，量綱為L－1，用格林函數(shù)解波動方程［19－20］（1）、（2）得到Lippmann－Schwinger積分方程：

其中方程第一項表示由震源以背景速度向外傳播的發(fā)散波，第二項是總場中由于介質(zhì)速度相對背景速度的擾動而產(chǎn)生的波場.求解方程（3）可得到Born序列：U＝U0＋U1＋U2＋U3＋…，設(shè)觀測點的位置為ri，則有：

其中，U1的物理含義是在震源rs產(chǎn)生的波以背景速度c0傳播，入射到擾動區(qū)域r′處發(fā)生散射，散射強度為k2α（r′），之后仍以c0自由傳播到觀測點的球面波之和.U2的物理含義是在震源rs產(chǎn)生的波以背景速度c0傳播，入射到擾動區(qū)域r″處發(fā)生散射，散射強度為k2α（r″），散射后傳播到擾動區(qū)域r′處再次發(fā)生散射，散射強度為k2α（r′），之后仍以c0自由傳播到觀測點的球面波之和.同理Un的物理含義是在震源rs產(chǎn)生的波在擾動區(qū)域內(nèi)r′處發(fā)生n次散射后，在觀測點記錄的波［21－22］.

式（4）、（5）中背景場格林函數(shù)分別為：

其中：ki、ks和km分別是觀測點、震源點、擾動點處波場在水平面內(nèi)的波數(shù)量綱均為L－1.本文散射波正演只計算一次散射，即Born近似.

2.2 散射波偏移方法

將地下不均勻地質(zhì)體看成由散射點組成，地震波從震源傳播到散射點，被散射點散射后，逆向散射傳播到地面所有的接收點，因此，每個接收點均能接收到來自散射點的能量.按照Bancroft和Geiger提出的等效偏移距方法［23－25］，將符合雙平方根方程的散射能量在沒有時移的情況下，映射到符合單平方根方程的共散射點道集上，并且散射能量沿雙曲線形態(tài)分布；在共散射點道集上進行速度分析得到偏移速度，應(yīng)用Kirchhoff積分方程對散射能量進行求和，得到偏移剖面.關(guān)于這一方法已有詳細介紹，本文不予說明［26］.

3 算 例

3.1 模 型

陷落柱幾何形態(tài)一般頂部為圓錐狀，下部為圓柱狀.以山西潞安礦區(qū)常見的典型陷落柱為參照建立陷落柱模型如圖1所示.該地區(qū)上覆第四紀(jì)黃土，煤系地層主要分布在石炭－二疊紀(jì)，模型參數(shù)如表1所示.陷落柱頂距上覆黃土層50m，陷落柱高250m，陷落柱頂部傾角從上至下依次為18°（深250～255m），30°（深255～265m）和65°（深265～300m），陷落柱柱壁從深300m延伸到500m，柱壁傾角達87°，因網(wǎng)格離散，部分為90°，陷落柱內(nèi)充填物為隨機介質(zhì)，波速范圍為2200～2800m／s，低于圍巖速度.

正演模擬采樣相移法［20］，模擬參數(shù)為：模型長1800m，深600m；網(wǎng)格剖分大小為5m×5m；采樣頻率為1ms，記錄時間為700ms；子波主頻為45Hz，頻帶范圍為1～100Hz；全排列接收，道間距為5m，共360道；炮間距為30m，共60炮.

3.2 散射波場特征

合成的炮集地震記錄如圖2所示，炮點位置分別位于810m、900m和1020m的位置.從圖中可知散射波場在炮集上的特征為：

圖1 圓錐狀陷落柱模型Fig.1 Sunken pillar model

表1 陷落柱模型參數(shù)表Table 1 Elastic parameters of sunken pillar in coal seams

（1）可以接收到來自陷落柱陡傾界面的弱散射，不存在地震波散射盲區(qū)，散射波的能量主要集中于反射波射線方向；

（2）煤層斷陷點左右兩支繞射波存在明顯的干涉區(qū)；

（3）斷陷點散射波與陷落柱邊界散射波也存在相互干涉.

如圖2所示，在炮集shot28上，煤層1、2右側(cè)斷陷點散射波與陷落柱左側(cè)邊界弱散射波相互干涉明顯.在炮集shot30上，煤層1、2左側(cè)斷陷點散射波與陷落柱右側(cè)邊界弱散射波相互干涉明顯.在炮集shot31上，煤層1、2兩側(cè)斷陷點散射波與陷落柱左右邊界弱散射波相互干涉，由于在陷落柱的正上方，能量較弱；但由于陷落柱內(nèi)部相對低速度及煤層反射缺失，造成煤層反射同相軸下陷和強度減弱，也呈現(xiàn)明顯的波場異常.陷落柱內(nèi)隨機介質(zhì)散射波在炮集上不明顯.

從陷落柱模型散射波單炮記錄中抽取CSP道集記錄，如圖3所示.3個CSP道集位置分別為702.5m、877.5m、897.5m.對比模型可知，散射波在CSP道集上的波場特征如下：

（1）CSP道集具有很高的覆蓋次數(shù)（200次），遠比CMP道集覆蓋次數(shù)高（CMP道集覆蓋次數(shù)為60）；（2）CSP道集中包含其對應(yīng)位置下部很大范圍內(nèi)地質(zhì)模型的有效信息.雙曲線頂點橫向、縱向分布范圍與模型中傾斜界面橫向、縱向分布范圍一致；（3）CSP道集等效于零偏移距自激自收剖面.所以在每個CSP道集上都能辨認(rèn)出靠近散射點部位地下模型的大致構(gòu)造形態(tài)；（4）來自水平地層界面反射波是散射波定向排列、相干加強的結(jié)果，因此水平反射界面的散射能量最強.陡傾角地層或地質(zhì)體界面散射能量較弱、分布范圍較小，具弱散射波特征，且其雙曲線頂點連線位置為傾斜界面的位置.

圖3中紅點表示散射點道集位置；綠點表示反射界面散射點；黑點表示陷落柱界面散射點.綠色箭頭指示反射界面散射點形成的散射波雙曲線；黑色箭頭指示陷落柱界面散射點形成的散射波雙曲線.

3.3 散射偏移成像

對CSP道集進行速度分析，拾取散射點偏移速度，見圖4，可以形成偏移速度場.

散射波偏移成像結(jié)果如圖5所示，由基于CSP道集的等效偏移距偏移剖面與基于CMP道集的反射疊前偏移剖面對比可見，等效偏移距偏移剖面成像效果具有比較優(yōu)勢：陷落柱頂點、傾斜界面以及底邊界能夠清晰成像，沒有畫弧現(xiàn)象，陷落柱內(nèi)部隨機介質(zhì)成像清晰.

4 實際數(shù)據(jù)試驗

為了驗證散射成像方法識別煤田陷落柱相對于反射地震的優(yōu)勢，我們在山西晉城某煤礦找到一條過已知開挖驗證陷落柱的二維地震測線，通過常規(guī)縱波反射地震疊后偏移和本文提出方法成像效果的對比，說明了散射成像相比于反射地震的偏移是有優(yōu)勢的.如圖6所示，其中圖6a是常規(guī)反射地震疊后偏移剖面，圖6b是基于散射波的成像剖面.在陷落柱的位置處（圖中箭頭指示），反射地震信號煤層反射連續(xù)，只煤層上下各存在不連續(xù)反射痕跡；而散射成像不僅可以精確刻畫陷落柱所引起的不連續(xù)反射，而且兩側(cè)柱面的形態(tài)及縱向延伸可以清晰地識別.

圖6 已知陷落柱的反射偏移（a）與散射成像（b）對比Fig.6 Contrast between reflect migration（a）and scattered imaging（b）of a sunken pillar

5 結(jié) 論

本文從散射波的角度討論煤田陷落柱的成像問題.從波動方程出發(fā)對陷落柱模型進行了散射波正演模擬；并進一步對模擬和實際的數(shù)據(jù)應(yīng)用等效偏移距方法進行偏移成像，獲得以下認(rèn)識：

（1）在反射地震的盲區(qū)內(nèi)；地面上能觀測到陡傾角的散射波場，但散射波的能量主要集中于反射波射線方向.

（2）陷落柱的相鄰斷陷點之間、以及相鄰散射點之間的散射波存在相互干涉現(xiàn)象.

（3）CSP道集具有很高的覆蓋次數(shù)，能夠提高信噪比，有利于弱散射信息的提取和疊加.

（4）CSP道集等效于零偏移距自激自收剖面.通過CSP道集可以判斷陷落柱等不均勻地質(zhì)體的形態(tài).

（5）來自水平界面散射點的散射波能量強于陡傾界面散射點的散射波能量.

（6）等效偏移距方法能夠合理、準(zhǔn)確地提取散射波場信息，基于等效偏移距疊前時間偏移方法是對陷落柱等不均勻地質(zhì)體成像的有效方法之一.

（References）

［1］ 張茂林，尹尚先.華北型煤田陷落柱形成過程研究.煤田地質(zhì)與勘探，2007，35（6）：26－29.Zhang M L，Yin S X.Forming process of subsided column in coalfields of North China.Coal Geology ＆Exploration（in Chinese），2007，35（6）：26－29.

［2］ 方 沛，王海龍.巖溶陷落柱成因分布規(guī)律淺析.煤炭工程，2008，（11）：84－85.Fang P，Wang H L.Karst in the fall of distribution of causes.Coal Engineering（in Chinese），2008，（11）：84－85.

［3］ 尹尚先，武 強，王尚旭.華北煤礦區(qū)巖溶陷落柱特征及成因探討.巖石力學(xué)與工程學(xué)報，2004，23（1）：120－123.Yin S X，Wu Q，Wang S X.Studies on characters and forming mechanism of karstic collapse columns at mine area of north China.Chinese Journal of Rock Mechanics and Engineering（in Chinese），2004，23（1）：120－123.

［4］ Wang J M，Liu X W，Liu W S，et al.The inner circulation mechanism in sinking column formation in North China coal field—A case study in Fengfeng mining areas.Journal of Coal Science ＆Engineering（China），2007，13（1）：32－36.

［5］ 曹志勇，王 偉，楊德義等.煤田陷落柱波場模擬與分析.太原理工大學(xué)學(xué)報，2008，39（spec）：841－843.Cao Z Y，Wang W，Yang D Y，et al.Coalfield subsided column wavefield modeling and analysis.Journal of Taiyuan University of Technology（in Chinese），2008，39（spec）：841－843.

［6］ 曹志勇，楊德義.陷落柱的高斯射線束法模擬.中國煤炭地質(zhì)，2008，20（6）：63－65，69.Cao Z Y，Yang D Y.Modeling of subsided column with Gaussian beam method.Coal Geology of China（in Chinese），2008，20（6）：63－65，69.

［7］ 楊德義，王 赟，王 輝.陷落柱的繞射波.石油物探，2000，39（4）：82－86.Yang D Y，Wang Y，Wang H.Diffraction waves from fallen pillars.Geophysical Prospecting for Petroleum（in Chinese），2000，39（4）：82－86.

［8］ 吳守華，周國興，楊素霞等.陷落柱地震響應(yīng)特征分析.煤田地質(zhì)與勘探，2004，32（3）：52－54.Wu G H，Zhou G X，Yang S X，et al.Analysis for the seismic response characteristics of collapse column.Coal Geology ＆Exploration（in Chinese），2004，32（3）：52－54.

［9］ 孟憲民，蔣維平.波阻抗反演技術(shù)及其應(yīng)用效果.中國煤炭地質(zhì)，2009，21（1）：52－54.Meng X M，Jiang W P.Elastic impedance inversion technology and its application effect.Coal Geology of China（in Chinese），2009，21（1）：52－54.

［10］ 寧建宏，張廣忠.陷落柱的地震識別技術(shù)及其應(yīng)用.煤田地質(zhì)與勘探，2005，33（3）：64－67.Ning J H，Zhang G Z.Seismic identification technique and its application of collapse column.Coal Geology ＆Exploration（in Chinese），2005，33（3）：64－67.

［11］ 楊曉東.三維地震勘探在構(gòu)造復(fù)雜地區(qū)的應(yīng)用.中國煤炭地質(zhì)，2008，20（3）：62－64.Yang X D.Applied results of 3Dseismic prospecting in shallow coal seam buried areas.CoalGeologyofChina（in Chinese），2008，20（3）：62－64.

［12］ 程建遠，王璽瑞，郭曉山等.東龐礦突水陷落柱三維地震處理效果與對比.煤田地質(zhì)與勘探，2008，36（1）：62－65.Cheng J Y，Wang X R，Guo X S，et al.3Dseismic processing for subsidence column detection applied in Dongpang coal mine.Coal Geology ＆Exploration（in Chinese），2008，36（1）：62－65.

［13］ 黃聯(lián)捷，吳如山.垂向非均勻背景多頻背向散射層析成像.地球物理學(xué)報，1994，37（1）：87－100.Huang L J，Wu R S.Multifrequency backscattering tomography in a vertically varying background.Chinese J.Geophys.（Acta Geophysica Sinica）（in Chinese），1994，37（1）：87－100.

［14］ 董其成，趙鴻儒.非均質(zhì)礦體地震模型的衍射、轉(zhuǎn)換波及次生波研究.地球物理學(xué)報，1994，37（2）：228－238.Dong Q C，Zhao H R.Diffraction，converted waves and derivative wave research of anisotropic geologic body by seismic model.Chinese J.Geophys.（Acta Geophysica Sinica）（in Chinese），1994，37（2）：228－238.

［15］ 陶智勇.彈性波逆散射微擾論.地球物理學(xué)報，2005，48（4）：918－923.Tao Z Y.Elastic wave inverse scattering perturbation theory.Chinese J.Geophys.（in Chinese），2005，48（4）：918－923.

［16］ Wu R S，Aki K.Scattering characteristics of elastic waves by an elastic heterogeneity.Geophysics，1985，50（4）：582－595.

［17］ Wu R S，Huang L J.Scattered field calculation in heterogenous media using aphase－screen propagator：Expanded Abstracts of the Technical Program.SEG 62nd Annual Meeting，1992.1289－1292.

［18］ Wu R S，Huang L J，Xie X B.Backscattered Wave calculation using the DE Wolf approximation and a phasescreen propagator：Expanded Abstracts of the Technical Program.SEG 65th Annual Meeting，1995.1293－1296.

［19］ Weglein A B，Gssparotto F A，Carvalho P M，et al.An inverse－scattering series method for attenuating multiples in seismic reflection data.Geophysics，1997，62（6）：1975－1989.

［20］ Carcione J M，Herman G C，Ten Kroode A P E.Seismic modeling.Geophysics，2002，44（6）：1304－1087.

［21］ 杜世通.地震波動力學(xué)理論與方法.東營：中國石油大學(xué)出版社，2008.57－83.Du S T.Seismic Wave Dynamics Theory and Methods.Dongying：Petroleum University Press，2008.57－83.

［22］ Innanen K A.Born series forward modelling of seismic primary and multiple reflections：an inverse scattering shortcut.Geophys.J.Int.，2009，177：1197－1204.

［23］ Bancroft J C，Geiger H D.Equivalent offsets and CRP gathers for prestack migration.64th.Ann.Internat.Mtg.，Soc.Expl.Geophys.Expanded Abstracts，1994.672－675.

［24］ Bancroft J C，Wang S，F(xiàn)oltinek D S，et al.Common scatterpoint（CSP）prestack migration：CSEG National Convention.Expanded Abstracts，1995.47－48.

［25］ Bancroft J C，Geiger H D，Margrave G F.The equivalent offset method of prestack time migration.Geophysics，1998，63（6）：2042－2053.

［26］ 王 偉，尹軍杰，劉學(xué)偉等.等效偏移距方法及應(yīng)用.地球物理學(xué)報，2007，50（6）：1823－1830.Wang W，Yin J J，Liu X W，et al.Equivalent offset method and its application.Chinese J.Geophys.（in Chinese），2007，50（6）：1823－1830.

Numerical simulation and imaging of scattered wave of sunk pillar in coal seam

CAO Zhi－Yong1，WANG Wei2＊，WANG Yun3

1 Institute of Geology ＆Geophysics，Chinese Academy of Sciences，Beijing100029，China
2 State Key Laboratory of Resources and Environmental Information System，Institute of Geographic Sciences and Natural Resource Research，Chinese Academy of Sciences，Beijing100101，China
3 State Key Laboratory of Ore Deposit Geochemistry，Institute of Geochemistry，Chinese Academy of Sciences，Guiyang550002，China

The scattered wave theory has built the basis of scattered wave imaging of the sunken pillar in coal field seismic exploration.When seismic waves meet with the scatter point underground，the waves propagate around and generate scattered waves.And then reflection is formed by enhancing energy interference of scattered waves.Through the numerical modeling，we can receive the scattered waves from the steep dip interface of the sunken pillar in coal seam.After researching wave characteristic on common scatter point（CSP）gathers，we can recognize the scattered waves，and obtain the information of geologic structure underground，including thestructure position and the outline of heterogeneous geologic body.Tested with field seismic data，the result of scattered wave imaging indicates that the equivalent offset method（EOM）can reasonably and accurately use the information of scattered waves，and the pre－stack time migration based on EOM is an effective method for the imaging of sunken pillars and inhomogeneous geologic bodies.

Sunken pillar，Scattered wave，CSP，Equivalent offset，Imaging

10.6038／j.issn.0001－5733.2012.05.032

P631

2011－04－06，2011－09－28收修定稿

科技部973項目（2006CB202207）、國家重大科技專項（2008ZX05035－001－003、2008ZX05035－003－006HZ.2008ZX05008－006－004）和地下信息探測技術(shù)與儀器教育部重點實驗室（中國地質(zhì)大學(xué)，北京）開放課題項目（GDL0803）聯(lián)合資助.

曹志勇，男，1980年生，中國科學(xué)院地質(zhì)與地球物理研究所在讀博士生，主要從事地震各向異性理論方法與多分量地震技術(shù)研究.E－mail：caozhy09＠gmail.com

＊通訊作者王偉，男，1972年生，河北涿州人，博士，現(xiàn)在中國科學(xué)院地理科學(xué)與資源研究所從事勘探地震理論與方法的研究.E－mail：wang＿wei＠lreis.ac.cn

曹志勇，王偉，王赟.煤層陷落柱散射波數(shù)值模擬與成像.地球物理學(xué)報，2012，55（5）：1749－1756，

10.6038／j.issn.0001－5733.2012.05.032.

Cao Z Y，Wang W，Wang Y.Numerical simulation and imaging of scattered wave of sunk pillar in coal seam.Chinese J.Geophys.（in Chinese），2012，55（5）：1749－1756，doi：10.6038／j.issn.0001－5733.2012.05.032.

（本文編輯 劉少華）