封閉爆炸近場地震縱波及橫波質(zhì)點(diǎn)速度模型*

李學(xué)政，張誠鎏，劉文學(xué)

（西北核技術(shù)研究所，陜西 西安 710024）

封閉爆炸近場地震縱波及橫波質(zhì)點(diǎn)速度模型*

李學(xué)政，張誠鎏，劉文學(xué)

（西北核技術(shù)研究所，陜西 西安 710024）

基于第四紀(jì)黃土層小爆炸當(dāng)量系列封閉爆炸的近場地震觀測數(shù)據(jù)，對縱波和橫波分別進(jìn)行了分析處理，獲得了縱波和橫波的質(zhì)點(diǎn)速度，確定了縱波和橫波質(zhì)點(diǎn)速度現(xiàn)有模型的參數(shù)。綜合分析了有代表性的3個(gè)模型，得出現(xiàn)存模型本質(zhì)上的差別在于爆炸當(dāng)量指數(shù)不同。對此提出了雙極模型，并介紹了該模型參數(shù)的計(jì)算方法。對這幾種模型在描述觀測和實(shí)測數(shù)據(jù)方面的殘差進(jìn)行了比較。結(jié)果表明，雙極模型反演數(shù)據(jù)與觀測數(shù)據(jù)之間的殘差平方和最小，而傳統(tǒng)模型中最好的效果僅是接近雙極模型。

爆炸力學(xué)；質(zhì)點(diǎn)速度；封閉爆炸；近場地震；雙極模型

工程爆破地震效應(yīng)中有關(guān)的地運(yùn)動參數(shù)普遍受到關(guān)注［1－2］，這些參數(shù)直接和巖土松動與破碎、礦山開采、坑道工程和地震探礦等眾多領(lǐng)域密切相關(guān)。一些學(xué)者［2－5］對近地表爆炸地震的地面豎向振動速度地表效應(yīng)和隨機(jī)特性進(jìn)行了研究，給出了地面豎向振動速度的理論模擬結(jié)果，在隨機(jī)統(tǒng)計(jì)特性上與實(shí)驗(yàn)結(jié)果具有較好的一致性，而用短時(shí)傅立葉變換和基于小波變換理論的時(shí)頻分析方法對地震效應(yīng)的時(shí)間－頻率特征的分析結(jié)果表明，小波變換方法可以更好地給出爆炸地震效應(yīng)的細(xì)節(jié)信息。然而，在具體應(yīng)用領(lǐng)域中，特別是在爆炸拆除、爆炸安全距離估計(jì)和地震探礦中，不僅關(guān)心爆炸地震幅值大小，而且關(guān)注爆炸地震引起的振動參數(shù)隨觀測距離和爆炸當(dāng)量的變化。

已建立了眾多地運(yùn)動參數(shù)與爆炸當(dāng)量及觀測距離的關(guān)系，這些數(shù)據(jù)大多是通過礦山開采和人類工程采集的，所涉及的爆炸當(dāng)量是數(shù)十噸甚至數(shù)百噸TNT，觀測距離為幾千米至數(shù)百千米，爆炸能量往往有相當(dāng)大的部分消耗到空氣中，而有關(guān)完全封閉爆炸，特別是幾千克、數(shù)十千克直至數(shù)百千克的封閉爆炸地震資料相對缺乏。爆炸源所處的空間位置不同，爆炸能量在空氣和介質(zhì)中的能量分布差別較大，如果炸藥中心位于地表水平面上，分別有65%和35%的能量傳到空氣和土中；如果炸藥的下底位于地表水平面上，分別有78%和22%的能量傳到空氣和土中。而且爆炸地震傳播對介質(zhì)的依賴性很強(qiáng)，對于不同的介質(zhì)條件，地震波傳播規(guī)律也存在顯著差異。由于考慮到地震破壞效應(yīng)，傳統(tǒng)上往往考慮的是最大速度模型，而爆炸最大速度往往不是同類波形，數(shù)據(jù)之間相關(guān)性較低。

本文中，基于黃土層爆炸當(dāng)量千克級至百千克級一系列封閉爆炸的地震觀測數(shù)據(jù)，分別對縱波和橫波進(jìn)行分析處理，計(jì)算爆炸地震縱波和橫波質(zhì)點(diǎn)速度隨爆炸當(dāng)量和觀測距離的變化，確定質(zhì)點(diǎn)速度模型。并提出雙極模型，該模型能夠同時(shí)實(shí)現(xiàn)相關(guān)系數(shù)為極大、實(shí)測與模型計(jì)算數(shù)據(jù)間殘差平方和極小。詳細(xì)介紹該模型參數(shù)的計(jì)算方法，并對這幾種模型在描述觀測數(shù)據(jù)方面的殘差進(jìn)行比較。

1 數(shù)學(xué)模型的建立

1.1 模型理論分析

物理上爆炸地震的質(zhì)點(diǎn)振動參數(shù)與爆炸當(dāng)量、爆炸方式、爆炸源的埋深、源區(qū)介質(zhì)、傳播路徑、觀測點(diǎn)地質(zhì)條件等多種因素有關(guān)，這些因素作用的表現(xiàn)形式為爆炸地震波的激發(fā)、耦合或衰減等效應(yīng)。爆炸當(dāng)量為激發(fā)爆炸地震的能量源，顯然，爆炸當(dāng)量越大，激發(fā)出的應(yīng)力波和地震波越強(qiáng)。爆炸方式、源區(qū)介質(zhì)和爆炸源的埋深影響到能量耦合［6－8］，對于相同的爆炸當(dāng)量，因爆炸源所在的介質(zhì)不同，爆炸能量耦合到地層中的能量差別很大。一般，硬巖耦合系數(shù)高，軟巖耦合系數(shù)低；封閉式耦合系數(shù)高，噴發(fā)式耦合系數(shù)低。由于介質(zhì)的吸收和球面擴(kuò)散作用，地震波參數(shù)隨觀測距離增大而發(fā)生非線性衰減。于是，質(zhì)點(diǎn)振動參數(shù)p可表示為

式中：質(zhì)點(diǎn)振動參數(shù)p包括位移、速度或加速度，W 為爆炸當(dāng)量，R為觀測距離；η為能量耦合系數(shù)，是變量c1、c2和c3的函數(shù)，c1、c2和c3分別表示爆炸方式、爆炸埋深和源區(qū)介質(zhì)等因素。

1.2 現(xiàn)存模型

根據(jù)上面分析可知，爆炸地震參數(shù)除與爆炸當(dāng)量、傳播距離明顯有關(guān)外，還和能量耦合系數(shù)有關(guān)，而爆炸方式、爆炸埋深和源區(qū)介質(zhì)直接影響能量耦合系數(shù)。針對具體問題，研究者們提出眾多數(shù)學(xué)模型或稱經(jīng)驗(yàn)公式時(shí)，對封閉爆炸與拋灑爆炸、硬巖與黃土等分別進(jìn)行研究，從而簡化了模型。但歸結(jié)起來，具有代表性的模型有以下3種。

P.B.Attwell等［9］提出了地震動最大速度模型Ⅰ為

式中：K 為場地系數(shù)，K＝0.013～0.148，α＝0.64～0.96。

郝保田［10］根據(jù)高能封閉爆炸地震動觀測數(shù)據(jù)，提出了最大速度模型Ⅱ?yàn)?/p>

式中：K＝8.2，α＝2.23。

謝毓壽等［11］基于大量的工業(yè)化爆數(shù)據(jù)，對經(jīng)驗(yàn)關(guān)系進(jìn)行了深入研究。根據(jù)炸藥量為50g的實(shí)驗(yàn)數(shù)據(jù)，得出質(zhì)點(diǎn)速度隨觀測距離的衰減因子為約－1.8，采用爆炸當(dāng)量0.1～2.0kg的爆炸地震觀測數(shù)據(jù)，觀測距離在35～404m之間，爆炸當(dāng)量指數(shù)為約0.6。提出地震動最大速度模型Ⅲ為

上式可改寫為v＝KWαRβ，場地為堅(jiān)硬基巖時(shí)，K1＝2.2～2.6，α＝0.60，β＝－1.8。

1.3 雙極模型

式（2）～（4）有一個(gè)共同特點(diǎn)，質(zhì)點(diǎn)速度都是爆炸當(dāng)量與觀測距離的函數(shù)，但又分為2類：（1）質(zhì)點(diǎn)速度是以爆炸當(dāng)量和觀測距離為變量的函數(shù)，如式（4）；（2）質(zhì)點(diǎn)速度是爆炸當(dāng)量和觀測距離的復(fù)合函數(shù)。其復(fù)合函數(shù)又有2種，詳見式（2）和（3）。根據(jù)上述公式及參數(shù)的分布范圍計(jì)算可知，爆炸當(dāng)量指數(shù)在0.33～0.96之間。上述3個(gè)模型中，由于模型Ⅲ是一個(gè)二元函數(shù)，當(dāng)爆炸當(dāng)量范圍較寬、實(shí)驗(yàn)次數(shù)較多時(shí)，統(tǒng)計(jì)處理結(jié)果相對可靠。模型Ⅰ和Ⅱ的復(fù)合變量形式上是不同的，但是，通過恒等變形，均可變?yōu)橥环N形式。由此看出，模型比例距離的爆炸當(dāng)量指數(shù)差別較大。換句話說，由于在比例距離計(jì)算中所選取的爆炸當(dāng)量指數(shù)不同，才出現(xiàn)了不同的模型，在預(yù)測數(shù)據(jù)時(shí)，精度也不同。模型建立的目的就是使模型預(yù)測的數(shù)據(jù)和觀測數(shù)據(jù)盡可能地接近，既然決定爆炸當(dāng)量模型的本質(zhì)是爆炸當(dāng)量指數(shù)的選取，在此對上述模型Ⅰ和Ⅱ進(jìn)行統(tǒng)一，選取最優(yōu)的爆炸當(dāng)量指數(shù)，使數(shù)據(jù)間相關(guān)系數(shù)最大，使模型反映的數(shù)據(jù)與觀測數(shù)據(jù)之間的殘差平方和最小，簡稱雙極模型Ⅳ，可表示為v＝K（Wα／R）β 。

圖1 爆炸源與測點(diǎn)分布示意圖Fig.1Layout of explosion sources and observation points

2 觀測數(shù)據(jù)

2.1 觀測概況

所用的地震計(jì)有2種，3D-lite和 CMG-40T。3D-lite為三分向速度地震計(jì)，速度頻響在1～80Hz范圍內(nèi)為一平直線，速度換能靈敏度為400V／（m／s），地震計(jì)噪聲小于3nm／s；CMG-40T速度頻響為0.033～50Hz，速度換能靈敏度為2kV／（m／s）。

在地層為第四紀(jì)黃土層地區(qū)10km2多的范圍內(nèi)，在爆炸源的東南－西北方向，距爆炸源0.60～2.50 km的范圍內(nèi)，布設(shè)了10多個(gè)地震測點(diǎn)，測點(diǎn)與爆炸源呈半包圍形式，爆炸源和測點(diǎn)分布見圖1。

爆炸源被埋放在地下20m處，上覆土層封閉，爆炸當(dāng)量為千克級至百千克級。為了研究地震波質(zhì)點(diǎn)速度隨爆炸當(dāng)量的變化，應(yīng)保證爆炸源處的介質(zhì)和位置“不變”，盡可能縮短爆炸源位置間的距離，所有爆炸源布設(shè)在半徑約50m的較小區(qū)域內(nèi)。

在野外黃土介質(zhì)中進(jìn)行了多次系列封閉小爆炸當(dāng)量化爆，獲得了一批地震數(shù)據(jù)，除個(gè)別測點(diǎn)因本底噪聲較大或震中距較遠(yuǎn)，記錄資料信噪比較小外，大多數(shù)觀測記錄波形完整、初至清晰。

2.2 觀測結(jié)果

近場地震記錄波形如圖2所示，宏觀上可以看出，所有的波形頻率均在10～30Hz范圍內(nèi)，觀測距離越長，波形的持續(xù)時(shí)間也越長，頻率也越低。每條記錄明顯發(fā)育2個(gè)震相，前段為縱波，后段為橫波。對記錄進(jìn)行了分析處理，獲得了縱波質(zhì)點(diǎn)速度和橫波質(zhì)點(diǎn)速度。在此列出了3次爆炸記錄波形處理結(jié)果，如表1所示，R為觀測點(diǎn)到爆炸源的距離，v為質(zhì)點(diǎn)速度，下標(biāo)“l(fā)”和“t”分別代表縱波和橫波，“50、8、100”分別代表50、8、100kg的爆炸當(dāng)量。

圖2 爆炸地震垂直向波形Fig.2 Seismic waves from an explosion earthquake

表1 黃土層中B2小爆炸當(dāng)量封閉爆炸地震動垂向速度部分?jǐn)?shù)據(jù)Table 1 Partial vertical velocities of low-yield explosion earthequake in soil for explosion B2

3 傳統(tǒng)模型

3.1 資料預(yù)處理

場地地震觀測點(diǎn)未變，只是爆炸源位置發(fā)生了少量變化，爆炸方式完全密封，場地的本底噪聲為約0.5μm／s，總體上觀測數(shù)據(jù)質(zhì)量是好的。宏觀上，測點(diǎn)質(zhì)點(diǎn)速度隨爆心距的增大而變小。但是，由于自然介質(zhì)的復(fù)雜性、爆炸源工程上實(shí)現(xiàn)的技術(shù)差異，地震波傳播過程中，數(shù)據(jù)還是具有一定的分散性。由于測點(diǎn)局部地質(zhì)條件所致，如2號爆炸，觀測點(diǎn)3和觀測點(diǎn)4，距離相差近一倍，而縱波質(zhì)點(diǎn)速度相差并不顯著，而個(gè)別測點(diǎn)的質(zhì)點(diǎn)速度出現(xiàn)異常，在模型建立時(shí)，將被視為奇點(diǎn)而剔除。

3.2 模型的建立

3.2.1 顯函數(shù)模型

模型Ⅲ為爆心距和爆炸當(dāng)量的顯函數(shù)模型，需要確定距離衰指數(shù)β和爆炸當(dāng)量指數(shù)α。只要所有觀測點(diǎn)在一條直線上，或在各向同性的介質(zhì)中，通過觀測一次或多次爆炸在不同距離上的地震數(shù)據(jù)，即可計(jì)算β。固定爆炸源位置，通過改變爆炸當(dāng)量，根據(jù)觀測資料，即可確定爆炸當(dāng)量指數(shù)α。由于爆炸沖擊效應(yīng)，爆炸源處介質(zhì)結(jié)構(gòu)受到破壞，工程技術(shù)上很難實(shí)現(xiàn)爆炸源的位置完全重復(fù)，因此，爆炸當(dāng)量指數(shù)α不易準(zhǔn)確確定。而表1中的各次爆炸源相對位置變化不大，僅為數(shù)10m，震中距離變化也不大，均在1%～2%之間，因此，可認(rèn)為爆炸源位置不變。

回歸結(jié)果如圖3、表2和表3所示，縱波質(zhì)點(diǎn)速度vl的距離衰減指數(shù)βl介于－2.13～－1.67之間，均值為約－2.0，相關(guān)系數(shù)γl在－0.97～－0.78之間；而橫波質(zhì)點(diǎn)速度vt的距離衰減指數(shù)βt介于－2.30～－1.97之間，均值為約－2.2，相關(guān)系數(shù)γt在－0.97～－0.92之間。比較發(fā)現(xiàn)，橫波質(zhì)點(diǎn)速度與觀測距離的相關(guān)性更顯著。

圖3 B2爆炸質(zhì)點(diǎn)速度與觀測距離的關(guān)系曲線Fig.3 Relation between particle velocities and observation distances for explosion B2

表2 距離衰減指數(shù)Table 2 Distance attenuation exponents

縱波和橫波質(zhì)點(diǎn)速度隨爆炸當(dāng)量的變化也不相同，縱波質(zhì)點(diǎn)速度隨爆炸當(dāng)量變化的指數(shù)αl在0.45～1.17之間，橫波波質(zhì)點(diǎn)速度隨爆炸當(dāng)量變化的指數(shù)αt在0.78～1.34之間。上述數(shù)據(jù)表明，質(zhì)點(diǎn)速度隨爆炸當(dāng)量變化的指數(shù)存在一定的分散性，橫波質(zhì)點(diǎn)速度隨爆炸當(dāng)量增加而增大較快。

根據(jù)3次爆炸地震觀測，計(jì)算了模型Ⅲ的參數(shù)，縱波Kl分別為1.17、1.32和1.23，均值Kl為約1.24，橫波Kt分別為0.44、0.36和0.42，均值為約0.41。計(jì)算數(shù)據(jù)表明，盡管爆炸當(dāng)量差別較大，K值仍具有較好的一致性。

于是模型Ⅲ的縱波和橫波質(zhì)點(diǎn)速度分別為

表3 爆炸當(dāng)量指數(shù)Table 3 Yield exponents

3.2.2 復(fù)合函數(shù)模型

復(fù)合函數(shù)模型以模型Ⅰ和模型Ⅱ?yàn)榇?，都是將爆炸?dāng)量和觀測距離的比例作為一個(gè)變量，模型的差別表現(xiàn)在計(jì)算爆炸當(dāng)量和觀測距離比例時(shí)，采用了不同的爆炸當(dāng)量指數(shù)。首先計(jì)算各個(gè)模型比例距離，然后再進(jìn)行數(shù)據(jù)擬和，模型中，比例距離的指數(shù)α和參數(shù)K的計(jì)算結(jié)果如表4和圖4所示，其中b為截距，σ為方差，下表“l(fā)”和“t”分別代表縱波和橫波。

表4 復(fù)合函數(shù)模型計(jì)算結(jié)果Table 4 Calculated results of compound function models

圖4 模型Ⅰ和模型Ⅱ縱波和橫波質(zhì)點(diǎn)速度的擬合曲線Fig.4 Fitted particle velocity curves of longitudinal and transversal waves based on modelsⅠandⅡ

4 雙極模型

雙極模型，稱為模型Ⅳ，有2個(gè)參數(shù)，分別是爆炸當(dāng)量指數(shù)和比例距離指數(shù)，通過2個(gè)步驟，可以計(jì)算數(shù)值的大小。首先設(shè)定爆炸當(dāng)量指數(shù)，對所有觀測數(shù)據(jù)進(jìn)行預(yù)處理，然后進(jìn)行回歸，計(jì)算出比例距離指數(shù)β、截距b、相關(guān)系數(shù)γ和質(zhì)點(diǎn)速度殘差平方和σ，然后改變爆炸當(dāng)量指數(shù)，再計(jì)算殘差和等。這樣反復(fù)循環(huán)計(jì)算，當(dāng)殘差和最小時(shí)所對應(yīng)的參數(shù)α和值β，即為所求的模型參數(shù)。

本模型具體計(jì)算時(shí)，設(shè)定指數(shù)α在1／1.0～1／3.0范圍內(nèi)按照一定間隔變化，并計(jì)算比例距離，然后進(jìn)行回歸。不同爆炸當(dāng)量指數(shù)所對應(yīng)模型的相關(guān)系數(shù)和數(shù)據(jù)殘差不同。相關(guān)系數(shù)是一條凹向下的具有極大值的曲線，殘差是一條凹向上的具有極小值的曲線，如圖5所示。圖5（c）和（d）分別為縱波相關(guān)系數(shù)γl和殘差σl曲線，爆炸當(dāng)量指數(shù)為1／2.4時(shí)，相關(guān)系數(shù)最大，殘差最小。圖5（e）和（f）分別為橫波相關(guān)系數(shù)γt和殘差σt曲線，最小殘差對應(yīng)的指數(shù)為1／2.0。圖5（a）和（b）分別第四紀(jì)黃土層縱波和橫波模型Ⅳ擬和曲線，詳細(xì)計(jì)算結(jié)果如表5所示，模型Ⅳ的縱波和橫波質(zhì)點(diǎn)速度分別為

圖5 關(guān)于模型Ⅳ的質(zhì)點(diǎn)速度、殘差和相關(guān)系數(shù)Fig.5 Particle velocities，remains and related coefficients on modelⅣ

表5 模型Ⅳ不同α計(jì)算結(jié)果部分?jǐn)?shù)據(jù)Table 5Partial calculation data on modelⅣin the case of different yield exponents

5 模型比較

在顯含模型數(shù)據(jù)處理中，爆炸當(dāng)量指數(shù)是通過固定爆炸當(dāng)量的位置不變得到的，這一假定，在物理上較難實(shí)現(xiàn)。其一，由于爆炸當(dāng)量不同，爆炸的埋深不同；其二，盡管在各次爆炸中，爆炸源位置間距很短，但位置不可能作到完全重復(fù)。該模型的參數(shù)與模型Ⅱ的參數(shù)接近。

表4中給出了模型Ⅰ和模型Ⅱ的計(jì)算參量，縱波的比例距離指數(shù)分別為0.866和2.162，而相關(guān)系數(shù)分別為0.943和0.942，方差分別為0.185和0.198。數(shù)據(jù)表明，這2個(gè)模型在描述觀測數(shù)據(jù)的變化規(guī)律方面存在一定的差異，模型Ⅰ描述各觀測數(shù)據(jù)效果較好。2個(gè)模型的橫波相關(guān)系數(shù)分別為0.975和0.951，方差分別為0.156和0.220，顯然，模型Ⅰ能較好地描述近場地震波的傳播規(guī)律。對于同一種模型而言，縱波和橫波相比，未表現(xiàn)出哪種震相的模型在反映質(zhì)點(diǎn)速度變化規(guī)律方面具有更強(qiáng)的優(yōu)勢。

由表5可知，模型Ⅳ中縱波爆炸當(dāng)量指數(shù)為1／2.4時(shí)，相關(guān)系數(shù)最大為0.948，殘差最小為0.178。與模型Ⅰ相比，該模型能夠更好地描述地震波質(zhì)點(diǎn)速度的變化規(guī)律。模型Ⅳ中當(dāng)橫波爆炸當(dāng)量指數(shù)為1／2.0時(shí)，殘差最小為0.157，而與模型Ⅰ相比，二者殘差基本一致。換句話說，無論是縱波還是橫波，在模型計(jì)算數(shù)據(jù)與觀測數(shù)據(jù)一致性方面，模型Ⅳ是最佳的，二者之間的殘差平方和最小，而傳統(tǒng)的有代表性的3個(gè)模型中最好的結(jié)果僅是接近雙極模型。

6 結(jié) 論

物理上爆炸地震的質(zhì)點(diǎn)振動參數(shù)與爆炸當(dāng)量、爆炸方式、源區(qū)介質(zhì)、傳播路徑、觀測點(diǎn)地質(zhì)條件等多種因素有關(guān)，這些因素作用的表現(xiàn)形式為彈性波激發(fā)、耦合和衰減。因此，質(zhì)點(diǎn)振動速度可表示為爆炸當(dāng)量、觀測距離和能量耦合的函數(shù)。

縱波和橫波質(zhì)點(diǎn)速度隨爆炸當(dāng)量呈指數(shù)變化，縱波當(dāng)量指數(shù)為0.80，橫波當(dāng)量指數(shù)為1.14；縱波和橫波質(zhì)點(diǎn)速度隨距離衰減指數(shù)不同，縱波衰減指數(shù)為2.0，橫波衰減指數(shù)為2.2。數(shù)據(jù)表明，橫波與縱波相比，質(zhì)點(diǎn)速度隨爆炸當(dāng)量增大而增大較快，隨距離增大而快速衰減。

綜合分析了有代表性的3個(gè)模型，得出模型本質(zhì)上的差別在于爆炸當(dāng)量指數(shù)不同。對此提出一種雙極模型，介紹了該模型參數(shù)的計(jì)算方法，并對這幾種模型在描述觀測和實(shí)測數(shù)據(jù)方面的殘差進(jìn)行了對比。無論是縱波還是橫波，雙極模型反演數(shù)據(jù)與觀測數(shù)據(jù)之間的殘差平方和最小，而傳統(tǒng)模型中最好的結(jié)果僅是接近雙極模型。

［1］劉愛文，俞言祥，付長華，等.50t科學(xué)爆破振動的衰減特征及地形效應(yīng)［J］.爆炸與沖擊，2010，30（1）：21-26.

LIU Ai-hua，YU Yan-xiang，F(xiàn)U Chang-hua，et al.Attenuation characteristics and topographic effect of a scientific explosion with 50texplosive［J］.Explosion and Shock Waves，2010，30（1）：21-26.

［2］劉晶波.局部不規(guī)則地形對地震地面運(yùn)動的影響［J］.地震學(xué)報(bào)，1996，18（2）：239-245.

［3］林大超，張奇，白春華.爆炸地震效應(yīng)的時(shí)頻分析［J］.爆炸與沖擊，2003，23（1）：32-36.

LIN Da-chao，ZHANG Qi，BAI Chun-hua.Time-frequency analysis of explosion seismic effects［J］.Explosion and Shock Waves，2003，23（1）：32-36.

［4］張雪亮，黃樹棠.爆破地震效應(yīng)［M］.北京：地震出版社，1981：159-161.

［5］Iyama J，Kuwamura H.Application of wavelets to analysis and simulation of earthquake motions［J］.Earthquake Engineering ＆ Structural Dynamics，1999，28（3）：255-272.

［6］李孝蘭.空腔解耦爆炸實(shí)驗(yàn)研究的基礎(chǔ)理論［J］.爆炸與沖擊，2000，20（2）：186-192.

LI Xiao-lan.Basic theory of decoupled explosions in cavities［J］.Explosion and Shock Waves，2000，20（2）：186-192.

［7］Haskell N A.A static theory of the seismic coupling of a contained underground explosion［J］.Journal of Geophysical Research，1961，66（9）：2937-2944.

［8］Glenn L A.Energy-density effects on seismic decoupling［J］.Journal of Geophysical Research，1993，98（2B）：1933-1942.

［9］Attewell P B，F(xiàn)armeer I W，Haslam D.Prediction of ground vibration parmeters from major quarry blasts［J］.Mining and Materials Engineering，1965：621-626.

［10］郝保田.地下核爆炸及其應(yīng)用［M］.北京：國防工業(yè)出版社，2002：227-229.

［11］謝毓壽，王耀文.工業(yè)爆破的地震效應(yīng)［J］.地球物理學(xué)報(bào)，1962，11（2）：154-161.

HSIEH Yü-show，WANG Yao-wen.Seismic effect of construction blasting［J］.Diqiu Wuli Xuebao，1962，11（2）：154-163.

Particle velocity models of longitudinal and transversal waves in the near field of sealed explosions＊

LI Xue-zheng，ZHANG Cheng-liu，LIU Wen-xue
（Northwest Institute of Nuclear Technology，Xi’an 710024，Shaanxi，China）

The seismic data of low-yield explosion earthquakes which took place in the Quaternary period soil were based on to calculate the particle velocities of the longitudinal and transversal waves and determine the parameters for the existent models.The three representative models were synthetically analyzed.It displays that the difference of the existent models lies in different yield exponents.Accordingly，a double-extreme model was proposed and the calculation method of the model parameters was introduced in detail.Comparison of these models shows that the data calculated by the double-extreme model is the most consistent with the real observation and the amplitude error between them is the smallest.The best results calculated by the traditional velocity models only approach these calculated by the double-extreme model.

mechanics of explosion；particle velocity；sealed explosion；near-field earthquake；double-extreme model

4February 2010；Revised 22July 2010

LI Xue-zheng，lixuezheng01＠yahoo.com.cn

（責(zé)任編輯 張凌云）

O384 國標(biāo)學(xué)科代碼：130·3599

A

1001-1455（2011）02-0196-08＊

2010-02-04；

2010-07-22

李學(xué)政（1957— ），男，博士，高級工程師。