邢廷棟，薛詩桂，黎小偉，索重輝，王輝明，焦艷艷

黃土塬區(qū)延安組煤層地震響應(yīng)特征物理模擬研究

邢廷棟1,2，薛詩桂1,2，黎小偉3，索重輝1,2，王輝明1,2，焦艷艷1,2

(1. 中國石化地球物理重點(diǎn)實(shí)驗(yàn)室，江蘇 南京 211103；2. 中國石化石油物探技術(shù)研究院，江蘇 南京 211103；3. 中國石化華北油氣分公司勘探開發(fā)研究院，河南 鄭州 450006)

針對黃土塬區(qū)延安組煤層地震響應(yīng)復(fù)雜問題進(jìn)行地震物理模擬研究。首先，研發(fā)適用于黃土塬地表干黃土層和煤層的模型材料，經(jīng)過多次試驗(yàn)，最終選用在硅橡膠中添加硅氣凝膠粉末的混合材料模擬疏松黃土塬地表層，選用在硅橡膠中添加超細(xì)碳粉的混合材料模擬低速低密度煤層；通過模具控制層位、逐層澆筑、三維雕刻起伏地層等方法，制作黃土塬區(qū)典型地質(zhì)結(jié)構(gòu)的三維地震物理模型，并開展地震物理模擬及地震成像分析。結(jié)果表明，延安組煤層與圍巖較大的波阻抗差異形成較強(qiáng)的反射振幅，對下部地層的成像有較強(qiáng)的屏蔽作用；較厚的多組煤層之間會形成層間多次波，影響下伏地層的成像。用地震振幅屬性對煤層進(jìn)行刻畫時(shí)，計(jì)算時(shí)窗大于40 ms更有利于煤層識別，但由于煤層存在調(diào)諧效應(yīng)，用地震屬性預(yù)測煤層厚度存在一定陷阱。

地震物理模擬；延安組煤層；黃土塬地表；振幅屬性

我國西部黃土塬地區(qū)煤炭資源十分富集，隨著國內(nèi)煤炭消耗量的不斷增長和東部煤炭資源的減少，煤炭資源勘探與開發(fā)的重心轉(zhuǎn)向西部地區(qū)[1-4]。而西部地表?xiàng)l件十分復(fù)雜，最主要的特征就是黃土塬，黃土塬地區(qū)土質(zhì)疏松、地表起伏大、彈性差、地震波傳播速度小、衰減大，對地震勘探造成很大的影響，使得黃土塬地區(qū)長期被地球物理學(xué)界視為“地震勘探禁區(qū)”[5-10]。地震正演模擬是理解地震波在地下介質(zhì)中的傳播特點(diǎn)、幫助解釋觀測數(shù)據(jù)的有效手段，在地震勘探中得到了廣泛的應(yīng)用[11]。地震正演模擬分為數(shù)值模型正演和物理模型正演，物理模型正演早在20世紀(jì)20年代就被用于地震勘探領(lǐng)域，但直到1974年著名的French模型誕生，地震勘探領(lǐng)域?qū)W者才真正認(rèn)識到物理模擬技術(shù)在三維地震勘探中的可行性及優(yōu)越性[12]。目前，對于煤層的地震模擬大多數(shù)是采用數(shù)值模擬方法，在物理模擬方面僅有少許學(xué)者進(jìn)行了研究，但并沒有考慮上覆黃土塬地層對煤層的影響[13-15]。筆者嘗試?yán)玫卣鹞锢砟M技術(shù)研究黃土塬區(qū)延安組煤層的地震響應(yīng)特征，從三維地質(zhì)建模到黃土塬區(qū)特性材料研究，再到三維地震物理模型的制作、數(shù)據(jù)采集與處理，最后對數(shù)據(jù)進(jìn)行解釋成像，系統(tǒng)分析延安組煤層成像特點(diǎn)，以及黃土塬區(qū)對煤層刻畫的影響，以期指導(dǎo)實(shí)際煤層地震勘探。

1 研究區(qū)域地質(zhì)條件

本次研究區(qū)主要位于鄂爾多斯盆地西南、天環(huán)坳陷南部，地層由下至上發(fā)育為中生代的三疊系、侏羅系、白堊系，新生代的古近系、新近系和第四系。其中，主要含煤地層位于侏羅系的延安組[16]，是此次模擬勘探的目的層。

1.1 地表?xiàng)l件

研究區(qū)地表為典型的黃土山地溝壑區(qū)，地勢總體走向?yàn)槲鞅备邧|南低，地形非常復(fù)雜，地表起伏大，地表高程為900～1 600 m。表層覆蓋有幾米到幾百米的黃土，且其厚度縱橫向隨機(jī)變化，總體呈現(xiàn)塬上厚，溝底薄的特點(diǎn)。表層黃土層土質(zhì)松散，從上到下大致可以分為干黃土(縱波波速為300～500 m/s)、濕黃土(縱波波速為800～1 000 m/s)和潛水面以下含水黃土(縱波波速為1 700～1 800 m/s)，表明煤層上覆層的速度較低。

1.2 延安組煤層地質(zhì)條件

通過測井資料分析可知，該區(qū)域在延8(KT_y8)、延9(KT_y9)等多組地層發(fā)育有煤層，且煤層厚度不均勻，變化較大[17]。其中延9煤層較厚，為6～15 m，延8煤層厚度較薄，為2～4 m。部分地區(qū)存在煤層分叉現(xiàn)象，發(fā)育多組煤層，而部分地區(qū)存在空白或弱反射，表明煤層并不是全區(qū)分布。煤層的速度和密度都較低，縱波速度約為2 500 m/s，密度約為1.6 g/cm3。

2 模擬材料研發(fā)

地震物理模型試驗(yàn)的基礎(chǔ)是地震物理模型的制作，其困難在于彈性參數(shù)的匹配與形態(tài)的可控，研究發(fā)現(xiàn)，環(huán)氧樹脂和硅橡膠的混合材料是最常用、最合適的模型材料[18-19]。然而這種混合材料有其自身的局限性，一方面是其速度變化范圍較小，一般為1 000～2 600 m/s，難以模擬低速層與高速層；另一方面是其均質(zhì)性很好，難以滿足非均質(zhì)性地層的要求。

2.1 煤層模擬材料

在地震剖面上，煤層常常表現(xiàn)為強(qiáng)振幅反射，延安組煤層具有速度低、密度小、厚度薄且分布不均勻等特點(diǎn)，厚度變化大，從幾米到幾十米不等，這就要求模擬材料有黏稠性、可塑性好、低速低密度等特性，而環(huán)氧樹脂和硅橡膠的混合材料無法達(dá)到這樣的要求，需要探尋新的材料或材料組合。

經(jīng)過多次試驗(yàn)，最終選用將1 000目(13 μm)的碳粉作為添加劑，以硅橡膠作為基質(zhì)按照1∶1的質(zhì)量比進(jìn)行混合，用該混合材料來模擬延安組煤層。首先碳粉是一種低密度的有機(jī)質(zhì)，而煤層也是一種非常好的有機(jī)質(zhì)，二者在特性上有相似之處；其次，純的硅橡膠速度僅為1 000 m/s左右，是一種非常好的低速材料，又是一種膠結(jié)劑，可以和碳粉混合后形成固體而不需要其他試劑；最后，加入碳粉末的硅橡膠材料變得較為黏稠，流動性差，有利于分布不均煤層的制作和成型。圖1a為1 000目(13 μm)的碳粉末。

2.2 黃土塬表層模擬材料

對于濕黃土和含水黃土層而言，常規(guī)的物理模型材料便可滿足模擬需求。但對于干黃土層而言，其主要的特點(diǎn)有：①較為疏松，縱橫波速度低，孔隙率大；②衰減大，對地震波具有很強(qiáng)的吸收作用；③厚度差異大，分布不均勻，地貌特征復(fù)雜。故而常規(guī)的模擬材料無法滿足這些特性，需要探索新的材料來模擬干黃土層。經(jīng)過大量的試驗(yàn)探索，本次最終選擇用一種名為硅氣凝膠材料的粉末作為添加劑，以硅橡膠為基質(zhì)的混合材料作為上覆低速高衰減黃土層的模擬材料。圖1b為塊狀硅氣凝膠。氣凝膠是一種固體物質(zhì)，是世界上密度最小的固體之一，常見的氣凝膠便是硅氣凝膠，其密度為 15 kg/m3，比空氣重十幾倍，此外，硅氣凝膠孔隙率非常高，約有95%，孔隙均為納米級，孔徑為20～70 nm，因此，其對于超聲波有較強(qiáng)的吸收作用。

制作干黃土地層物理模型的具體過程如下：①按比例稱取一定量的硅氣凝膠和硅橡膠液體(具體比例要根據(jù)黃土層的厚度和速度確定)，將其放置在溫度為40℃的恒溫箱中保溫2～3 h；②將硅氣凝膠取出，用小型粉碎機(jī)打磨約10 s(時(shí)間過長硅氣凝膠會變?yōu)榧?xì)粉，模擬效果不佳，過短則存在大顆粒，形成較大的非均質(zhì)體)；③取出硅橡膠液體，將打碎的硅氣凝膠粉混入，充分?jǐn)嚢枋苟呋旌暇鶆?；④依?jù)干黃土層的分布規(guī)律，將混合材料涂抹至物理模型表面。應(yīng)用這種材料制作的物理模型，其超聲波縱波速度變化范圍為600～1 000 m/s，拓展了原有地震物理模型材料在低速地層的限制。

3 物理模型構(gòu)建與實(shí)驗(yàn)采集

地震物理模擬是在實(shí)驗(yàn)室將野外的地質(zhì)構(gòu)造和地質(zhì)體按照一定的相似比模擬制作成物理模型，并用超聲波模擬野外地震波的一種模擬方法。具體來說，地震物理模擬就是依照地質(zhì)結(jié)構(gòu)中原體，應(yīng)用特定的相似準(zhǔn)則，縮制成模型，根據(jù)其速度、密度及結(jié)構(gòu)，在模型中復(fù)演與原體相似的天然狀況，進(jìn)行超聲地震實(shí)驗(yàn)，通過觀測獲取數(shù)據(jù)，然后再按照相似準(zhǔn)則將結(jié)果引申到原體，用以指導(dǎo)實(shí)踐。在這一系列實(shí)驗(yàn)過程中，首先要求模型和原型相比是相似的，其次要對實(shí)驗(yàn)數(shù)據(jù)做出合理的分析，從而得到真正可以反映原體實(shí)際情況的數(shù)據(jù)，原型是模型的基礎(chǔ)[20-23]。

3.1 物理模型設(shè)計(jì)與制作

首先，通過研究區(qū)域測井?dāng)?shù)據(jù)，獲取研究區(qū)內(nèi)各層位的層位數(shù)據(jù)，在此基礎(chǔ)上建立研究區(qū)的層位模型。模型以地表黃土層為頂面，以長9層(KT_c91)向下延伸至高程–1 100 m為底面，共16個(gè)層位(表1)。通過測井分層數(shù)據(jù)和區(qū)內(nèi)聲波測井?dāng)?shù)據(jù)，對研究區(qū)的各層位速度進(jìn)行了統(tǒng)計(jì)，并以此為各分層速度建模的基礎(chǔ)。由微測井和小折射分析，得到黃土塬上、斜坡、溝底的地層速度結(jié)構(gòu)變化和速度模型，從而完成對黃土層的地質(zhì)建模。黃土塬地貌大致有三層速度結(jié)構(gòu)，第一層和第二層為低(降)速度帶，速度范圍為400～1 000 m/s，深部黃土速度較為穩(wěn)定，速度為1 800 m/s。采用自然伽馬和聲波時(shí)差曲線交匯的方法識別煤層位置及厚度，分析發(fā)現(xiàn)研究區(qū)發(fā)育多組煤層。為簡化模型的制作，在模型設(shè)計(jì)時(shí)只考慮延8(KT_y8)和延9(KT_y9)兩組煤層，其中，延9煤層存在分叉現(xiàn)象，煤層分布不均勻，厚度變化范圍為6～15 m，縱波速度均為2 500 m/s。圖2為最終設(shè)計(jì)的地質(zhì)模型圖，為制作模型方便，縱向指標(biāo)由高程改為深度。模型尺寸與實(shí)際地質(zhì)體的比例為1∶20 000，模型的速度與實(shí)際速度比例為1∶2，采集數(shù)據(jù)的時(shí)間與野外采集的時(shí)間比例為1∶10 000，相應(yīng)的，實(shí)驗(yàn)室采集頻率與野外地震波頻率的比例為10 000∶1，即實(shí)驗(yàn)室主頻為170 kHz的超聲波代表野外主頻為17 Hz的地震波。物理模型各層參數(shù)設(shè)計(jì)見表1。

根據(jù)物理模型的設(shè)計(jì)，開始對三維物理模型進(jìn)行制作，模型的制作主要應(yīng)用三維雕刻起伏層位、精細(xì)模具控制煤層制作、雕刻反模澆筑黃土層等手段。圖3a為延9煤層制作，圖3b為延8煤層制作。依據(jù)煤層的厚度和分布范圍，通過精細(xì)模具控制2組煤層的制作，接著，開始對其上覆地層進(jìn)行澆筑，利用模具控制每個(gè)層位界面的制作，完成黃土層之下的所有地層澆筑。由于地表黃土層起伏較大，采用雕刻反模的方法進(jìn)行澆筑。最后是對干黃土層的澆筑，在澆筑過程中要保證干黃土層在塬上較厚，溝中較薄的特點(diǎn)，分多次澆筑。完成的模型如圖4所示，圖4a是煤層及其覆蓋層的三維物理模型，圖4b是最終的三維地震物理模型。

圖2 地質(zhì)模型原型

表1 三維地質(zhì)模型設(shè)計(jì)參數(shù)與測量參數(shù)

圖3 2組煤層物理模型

圖4 三維地震物理模型

3.2 地震物理模擬數(shù)據(jù)采集

完成物理模型制作后便可進(jìn)行地震物理模擬實(shí)驗(yàn)，將制作好的三維地震物理模型放入水槽中，應(yīng)用地震物理模擬采集儀器對模型進(jìn)行地震數(shù)據(jù)采集，在采集之前需要先設(shè)計(jì)好觀測系統(tǒng)，并由數(shù)據(jù)采集程序讀取和執(zhí)行。

為了更好地說明地表高衰減、低波速的干黃土層對煤層地震波場特征的影響，物理模擬實(shí)驗(yàn)過程中進(jìn)行了2次數(shù)據(jù)采集，完成目的層的制作后(圖4a)即進(jìn)行第一次三維地震數(shù)據(jù)采集，完成整個(gè)三維地震模型制作后(圖4b)再進(jìn)行一次三維地震數(shù)據(jù)采集。2次采集設(shè)計(jì)的觀測系統(tǒng)完全一致，采用40束30排炮4炮(共4 352炮)22線左135道右140道，其他參數(shù)見表2，依照觀測系統(tǒng)進(jìn)行地震數(shù)據(jù)的采集，最終獲得2套三維地震數(shù)據(jù)體。

表2 物理模型實(shí)驗(yàn)的觀測系統(tǒng)參數(shù)

4 數(shù)據(jù)處理、解釋及煤層響應(yīng)特征分析

4.1 物理模型數(shù)據(jù)處理與解釋

對2組地震數(shù)據(jù)分別進(jìn)行處理，可獲得2套偏移剖面。選取位于模型中間位置目的層的偏移剖面顯示，如圖5所示。從圖5可以看出，煤層有著非常強(qiáng)的波阻抗反射界面，圖中區(qū)域①為延9煤層的反射界面，該區(qū)域僅發(fā)育延9一組煤層，界面反射較強(qiáng)，其下除頁巖層尚可識別外，其他各層同相軸模糊且不連續(xù)；圖中區(qū)域②發(fā)育2組煤層，分別為延9煤層和延8煤層，反射界面振幅強(qiáng)且層間多次波發(fā)育，其下頁巖層成像雜亂，其余地層受到嚴(yán)重的屏蔽作用，基本無法識別。

圖5 目的層疊后偏移剖面(Inline)

圖6為2套數(shù)據(jù)在同一個(gè)Crossline方向上相應(yīng)的疊前時(shí)間偏移剖面，其中，圖6a為目的層疊前時(shí)間偏移剖面，圖6b為最終模型疊前時(shí)間偏移剖面。由于最終模型比目的層模型多了上覆地層及黃土塬地表層，故第二次采集的模擬數(shù)據(jù)對應(yīng)煤層的旅行時(shí)間較第一次多出約500 ms。由圖6可以看出，在沒有上覆層時(shí)，煤層地震響應(yīng)特征清晰，具有強(qiáng)振幅異常的特征，2組煤層之間發(fā)育層間多次波。當(dāng)有上覆地層及黃土塬地表層時(shí)，煤層的強(qiáng)振幅特性依然存在，但煤層較不連續(xù)，且受地表影響，形成傾斜界面，2組煤層之間的層間多次波消失不見。此外，煤層的下伏頁巖層和其他地層均無法識別?？梢婞S土塬地表層與煤層的疊加效應(yīng)對煤層下伏地層有著很強(qiáng)的干擾，嚴(yán)重影響其成像。

4.2 煤層地震響應(yīng)特征分析

地震屬性分析技術(shù)是地震儲層預(yù)測的一種主要手段[16]，本次使用屬性分析技術(shù)對延安組煤層進(jìn)行屬性刻畫。由于煤層存在異常振幅值，較其他方面的屬性而言，其對振幅有著更好的敏感性，因此，主要選取了最大振幅、最小振幅和均方根振幅3種屬性，并在不同時(shí)窗下進(jìn)行計(jì)算，以此識別和刻畫延安組煤層的邊界和厚度。共選取了4組時(shí)窗進(jìn)行振幅屬性計(jì)算，時(shí)窗分別為10、40、100、200 ms，圖7為目的層模型的振幅屬性計(jì)算結(jié)果，圖8為最終模型的振幅屬性計(jì)算結(jié)果。

從圖7可以看出，2套數(shù)據(jù)中，當(dāng)計(jì)算時(shí)窗較小(≤40 ms)時(shí)，各振幅屬性都難以準(zhǔn)確識別煤層的分布，也難以刻畫煤層的厚度；當(dāng)計(jì)算時(shí)窗較大(>40 ms)時(shí)，可以較好地刻畫煤層的分布。對比2套數(shù)據(jù)的振幅屬性，可以發(fā)現(xiàn)，黃土塬地層對煤層的識別有影響，主要表現(xiàn)在降低了地震數(shù)據(jù)的信噪比，使得煤層的邊界刻畫不準(zhǔn)確，但總體上煤層還是可識別的。對比3種不同振幅屬性，可以發(fā)現(xiàn)，在目的層模型中，三者識別的煤層差異較小，而且在大時(shí)窗條件下，基本可以準(zhǔn)確識別出煤層，而在有黃土塬地表層模型中，不同屬性識別出的煤層存在差異，相比而言，均方根振幅屬性刻畫的煤層更接近實(shí)際煤層的分布。

圖8 不同時(shí)窗計(jì)算下最終模型振幅屬性

5 結(jié)論

a.延安組煤層的速度和密度與圍巖相差較大，波阻抗差異明顯，在地震剖面上表現(xiàn)為強(qiáng)反射振幅特性，容易識別。

b. 煤層對下部地層的成像具有較強(qiáng)的屏蔽作用，尤其是發(fā)育多組煤層時(shí)，會產(chǎn)生層間多次波，使得其下伏地層成像模糊甚至無法成像。

c.由于煤層的調(diào)諧效應(yīng)，應(yīng)用振幅屬性對煤層進(jìn)行刻畫時(shí)，不同屬性和不同時(shí)窗刻畫的煤層厚度會有差別，故用地震屬性預(yù)測煤層厚度存在一定的陷阱；總體而言，大時(shí)窗(大于40 ms)計(jì)算得到的屬性更有利于煤層的刻畫。

d.黃土塬地表層對煤層的刻畫有著較大的影響，地表溝壑的存在使得煤層同相軸連續(xù)性變差，低速高衰減的干黃土層降低了煤層信號的信噪比，使得煤層的刻畫不準(zhǔn)確。

[1] 劉占勇，江濤，宋洪柱，等. 中國煤炭資源勘查開發(fā)程度分析[J]. 煤田地質(zhì)與勘探，2013，41(5)：1–5.

LIU Zhanyong，JIANG Tao，SONG Hongzhu，et al. Analysis of explorative and exploitative degree of China coal resources[J]. Coal Geology & Exploration，2013，41(5)：1–5.

[2] 程建遠(yuǎn)，聶愛蘭，張鵬. 煤炭物探技術(shù)的主要進(jìn)展及發(fā)展趨勢[J].煤田地質(zhì)與勘探，2016，44(6)：136–141.

CHENG Jianyuan，NIE Ailan，ZHANG Peng. Outstanding progress and development trend of coal geophysics[J]. Coal Geology & Exploration，2016，44(6)：136–141.

[3] 滕吉文，司薌，王玉辰. 我國化石能源勘探、開發(fā)潛能與未來[J]. 石油物探，2021，60(1)：1–12.

TENG Jiwen，SI Xiang，WANG Yuchen. Potential and future of fossil fuel exploration and development in China[J]. Geophysical Prospecting for Petroleum，2021，60(1)：1–12.

[4] 程建遠(yuǎn)，王壽全，宋國龍. 地震勘探技術(shù)的新進(jìn)展與前景展望[J].煤田地質(zhì)與勘探，2009，37(2)：55–58.

CHENG Jianyuan，WANG Shouquan，SONG Guolong. The new development and foreground expectation of seismic exploration[J]. Coal Geology & Exploration，2009，37(2)：55–58.

[5] 滕吉文，喬勇虎，宋鵬漢. 我國煤炭需求、探查潛力與高效利用分析[J]. 地球物理學(xué)報(bào)，2016，59(12)：4633–4653.

TENG Jiwen，QIAO Yonghu，SONG Penghan. Analysis of exploration，potential reserves and high efficient utilization of coal in China[J]. Chinese Journal of Geophysics(in Chinese)，2016，59(12)：4633–4653.

[6] 程建遠(yuǎn)，李寧，侯世寧，等. 黃土塬區(qū)地震勘探技術(shù)發(fā)展現(xiàn)狀綜述[J]. 中國煤炭地質(zhì)，2009，21(12)：72–76.

CHENG Jianyuan，LI Ning，HOU Shining，et al. Development status overview of seismic prospecting technology in loess tableland[J]. Coal Geology of China，2009，21(12)：72–76.

[7] 陳超群，高秦，何爭光，等. 鄂爾多斯盆地西南部巨厚黃土塬區(qū)非縱地震資料處理技術(shù)[J]. 煤田地質(zhì)與勘探，2017，45(1)：143–151.

CEHN Chaoqun，GAO Qin，HE Zhengguang，et al. The off-line seismic exploration and its application in huge-thick loess area in Ordos Basin[J]. Coal Geology & Exploration，2017，45(1)：143–151.

[8] 白萬山，劉田田，李紅桃. 黃土塬地區(qū)煤田地震勘探資料處理技術(shù)[J]. 煤田地質(zhì)與勘探，2014，42(4)：82–85.

BAI Wanshan，LIU Tiantian，LI Hongtao. Processing technologies of seismic data from the coalfield in loess plateau area[J]. Coal Geology & Exploration，2014，42(4)：82–85.

[9] 周俊杰，王雨，侯瑋. 黃土塬地區(qū)煤田三維地震綜合處理技術(shù)[J].地球物理學(xué)進(jìn)展，2016，31(5)：2299–2305.

ZHOU Junjie，WANG Yu，HOU Wei. 3D seismic comprehensive processing technology of coalfield in loess tableland[J]. Progress in Geophysics(in Chinese)，2016，31(5)：2299–2305.

[10] 陳超群，田媛媛，高秦，等. 基于隨機(jī)函數(shù)數(shù)據(jù)重構(gòu)的分頻異常振幅衰減技術(shù)在巨厚黃土塬區(qū)的應(yīng)用[J]. 石油物探，2019，58(5)：741–749.

CHEN Chaoqun，TIAN Yuanyuan，GAO Qin，et al. Frequency-division abnormal amplitude attenuation after data reconstruction based on random function and its application in the very thick loess tableland area，Ordos Basin[J]. Geophysical Prospecting for Petroleum，2019，58(5)：741–749.

[11] 熊曉軍，賀振華，黃德濟(jì). 三維波動方程正演及模型應(yīng)用研究[J].石油物探，2005，44(6)：554–556.

XIONG Xiaojun，HE Zhenhua，HUANG Deji. The application of 3D wave equation forward and modeling[J]. Geophysical Prospecting for Petroleum，2005，44(6)：554–556.

[12] 魏建新，牟永光，狄?guī)妥? 三維地震物理模型的研究[J]. 石油地球物理勘探，2002，37(6)：556–561.

WEI Jianxin，MOU Yongguang，DI Bangrang. Study of 3D seismic physical model[J]. Oil Geophysical Prospecting，2002， 37(6)：556–561.

[13] 戴世鑫. 基于物理模型的煤田地震屬性響應(yīng)特征的關(guān)鍵技術(shù)研究[D]. 北京：中國礦業(yè)大學(xué)(北京)，2012.

DAI Shixin. Research on key technology of response characteristics of seismic attributes based on the physical model[D]. Beijing：China University of Mining and Technology，2012.

[14] 韓堂惠，戴世鑫，李小華，等. 淮南煤系地層地震物理模型研究[J]. 煤炭學(xué)報(bào)，2011，36(4)：588–592.

HAN Tanghui，DAI Shixin，LI Xiaohua，et al. Seismic physical modeling research on coal measure strata in Huainan[J]. Journal of China Coal Society，2011，36(4)：588–592.

[15] 胡朝元. 薄煤層三維地震勘探技術(shù)：以淮南張集礦區(qū)11-2煤為例[D]. 北京：中國礦業(yè)大學(xué)(北京)，2011.

HAN Chaoyuan. Three-Dimensional seismic exploration technology for thin coal seam：Taking the 11-2 coal seam in Zhangji coalming for example[D]. Beijing：China University of Mining & Technology(Beijing)，2011.

[16] 陳曉智，湯達(dá)禎，許浩，等. 彬長礦區(qū)延安組煤層發(fā)育特征及其控制因素分析[J]. 中國礦業(yè)，2011，20(2)：110–113.

CHEN Xiaozhi，TANG Dazhen，XU Hao，et al. Development characteristics of coal seam and their controlling factors in Yan’an Formation in Binchang area[J]. China Mining Magazine，2011，20(2)：110–113.

[17] 王雪秋. 復(fù)雜近地表地震波響應(yīng)特征研究：以中國西部地區(qū)為例[D]. 長春：吉林大學(xué)，2009.

WANG Xueqiu. Study on seismic wave field under complex surface conditions of special physiognomy in West China[D]. Changchun：Jilin University，2009.

[18] 李智宏，朱海龍，趙群，等. 地震物理模型材料研制與應(yīng)用研究[J]. 地球物理學(xué)進(jìn)展，2009，24(2)：408–417.

LI Zhihong，ZHU Hailong，ZHAO Qun，et al. Study and materialization of new seismic physical model building materials[J]. Progress in Geophysics，2009，24(2)：408–417.

[19] 魏建新，狄?guī)妥? 地震物理模型中三維地質(zhì)模型材料特性研究[J]. 石油物探，2006，45(6)：586–590.

WEI Jianxin，DI Bangrang. Properties of materials forming the 3D geological model in seismic physical model[J]. Geophysical Prospecting for Petroleum，2006，45(6)：586–590.

[20] 趙群，馬國慶，宗遐齡. 超聲地震物理模型連續(xù)數(shù)據(jù)采集系統(tǒng)[J]. 地球物理學(xué)進(jìn)展，2004，19(4)：786–788.

ZHAO Qun，MA Guoqing，ZONG Xialing. Continuance data acquisition system of ultrasonic seismic physical modeling[J]. Progress in Geophysics，2004，19(4)：786–788.

[21] 狄?guī)妥?，魏建新，夏永? 三維地震物理模型技術(shù)的效果與精度研究[J]. 石油地球物理勘探，2002，37(6)：562–568.

DI Bangrang，WEI Jianxin，XIA Yongge. Study on effects and precision of 3D seismic physical model technique[J]. Oil Geophysical Prospecting，2002，37(6)：562–568.

[22] 趙鴻儒，王鐵男，唐文榜. 中國地球物理模型試驗(yàn)的發(fā)展[J]. 地球物理學(xué)報(bào)，1994，37(增刊1)：269–276.

ZHAO Hongru，WANG Tienan，TANG Wenbang. The developments of geophysical modeling[J]. Acta Geophysica Sinica，1994，37(Sup.1)：269–276.

[23] 王開燕，徐清彥，張桂芳，等. 地震屬性分析技術(shù)綜述[J]. 地球物理學(xué)進(jìn)展，2013，28(2)：815–823.

WANG Kaiyan，XU Qingyan，ZHANG Guifang，et al. Summary of seismic attribute analysis[J]. Progress in Geophysics，2013，28(2)：815–823.

Physical modeling of seismic response for the coal seams of Yan’an Formation in loess tableland of North China

XING Tingdong1,2, XUE Shigui1,2, LI Xiaowei3, SUO Chonghui1,2, WANG Huiming1,2, JIAO Yanyan1,2

(1. SINOPEC Key Laboratory of Geophysics, Nanjing 211103, China; 2. SINOPEC Geophysical Research Institute, Nanjing 211103, China; 3. Exploration and Development Research Institute, SINOPEC North China Company, Zhengzhou 450006, China)

A seismic physical simulation is carried out for the complex seismic response of the coal seam of Yan’an Formation in the loess tableland. First, the physical model materials for the dry layer near the surface of the loess tableland and the coal seam are developed. After several tests, a mixture of silicone rubber with silicon aerogel powder is selected to simulate the surface of the loess tableland, and the mixture with ultra-fine carbon powder added to silicone rubber is used to simulate the low-speed and low-density coal seam. Through the methods of molds control layer, layer by layer pouring, three-dimensional carving of undulating strata, a three-dimensional seismic physical model of the typical geological structure in the loess tableland is constructed, and seismic physical simulation and seismic imaging analysis are carried out. The results show that there is a strong reflection amplitude because of the great difference in impedance between the coal seam and surrounding rock of Yan’an Formation. Therefore, it has a strong shielding effect on the imaging of the lower strata. There are internal multiple waves in the multiple sets of coal seams, affecting the imaging of the underlying layer. The calculation time window greater than 40 ms is more favorable to coal seam identification when the coal seam is described by the seismic amplitude attribute. But there are some traps in predicting coal seam thickness by seismic attributes because of the coal seams tuning effect.

seismic physical modeling; Yan’an Formation coal seam; loess tableland surface; amplitude attribute

語音講解

P631

A

1001-1986(2021)06-0087-07

2021-08-20；

2021-10-22

國家自然科學(xué)基金企業(yè)創(chuàng)新發(fā)展聯(lián)合基金項(xiàng)目(U19B6003-004)

邢廷棟，1991年生，男，陜西榆林人，碩士，工程師，從事地震物理模擬實(shí)驗(yàn)及研究工作. E-mail：xingtd.swty@sinopec.com

薛詩桂，1980年生，男，河南南陽人，研究員，從事地震物理模擬及地震數(shù)據(jù)分析方法技術(shù)研究. E-mail：38473815@qq.com

邢廷棟，薛詩桂，黎小偉，等. 黃土塬區(qū)延安組煤層地震響應(yīng)特征物理模擬研究[J]. 煤田地質(zhì)與勘探，2021，49(6)：87–94. doi: 10.3969/j.issn.1001-1986.2021.06.010

XING Tingdong，XUE Shigui，LI Xiaowei，et al. Physical modeling of seismic response for the coal seams of Yan’an Formation in loess tableland of North China[J]. Coal Geology & Exploration，2021，49(6)：87–94. doi: 10.3969/j. issn.1001-1986.2021.06.010

移動閱讀

(責(zé)任編輯 聶愛蘭)