王敏超　李 欣　翁秀峰　張曉鵬　肖衛(wèi)國(guó)　沈旭峰　崔甲甲（西北核技術(shù)研究所　西安 710024）

?

譜估計(jì)方法在大型實(shí)驗(yàn)裝置地基振動(dòng)測(cè)量中的應(yīng)用

王敏超李 欣翁秀峰張曉鵬肖衛(wèi)國(guó)沈旭峰崔甲甲
（西北核技術(shù)研究所西安 710024）

摘要大型物理實(shí)驗(yàn)裝置為獲得高品質(zhì)的強(qiáng)流脈沖電子束流，通常對(duì)地基振動(dòng)都有較高的要求。詳盡的地質(zhì)振動(dòng)狀況分析是裝置選址的重要前提。本文針對(duì)某裝置待勘區(qū)域?qū)嶋H地質(zhì)條件和周圍人文環(huán)境，首先簡(jiǎn)要介紹了測(cè)點(diǎn)布設(shè)、振動(dòng)測(cè)試系統(tǒng)的基本組成，然后通過(guò)分階段測(cè)試比對(duì)試驗(yàn)對(duì)某裝置待勘區(qū)域進(jìn)行實(shí)際振動(dòng)測(cè)量，最后應(yīng)用譜分析的方法對(duì)區(qū)域內(nèi)不同時(shí)間段、不同方向的測(cè)量數(shù)據(jù)進(jìn)行分析，評(píng)估選址位置是否合理。數(shù)據(jù)結(jié)果表明，待勘區(qū)域在5μm·s-1以上速度被95%概率檢測(cè)，垂直方向位移振幅譜在185 nm，水平方向大于垂直方向約60 nm，物理實(shí)驗(yàn)表明，進(jìn)一步基建過(guò)程中應(yīng)采取適當(dāng)?shù)臏p震措施降低垂直位移幅度。

關(guān)鍵詞大型物理實(shí)驗(yàn)裝置，地基振動(dòng)，譜估計(jì)，位移譜，振幅譜

國(guó)家自然科學(xué)基金(No.11375145、No.41474035)資助

第一作者：王敏超，男，1981年出生，2013年于西安電子科技大學(xué)獲碩士學(xué)位，目前從事地震信號(hào)測(cè)量技術(shù)研究

Supported by National Natural Science Foundation of China(No.11375145,No.41474035)

First author:WANG Minchao,male,born in 1981,graduated from Xidian university with a master’s degree in 2013,focusing on the measurement technology of seismological signal

Application of ground vibration measurement in large experimental facilities based
on spectrum estimation methods

WANG MinchaoLI XinWENG XiufengZHANG XiaopengXIAO Weiguo SHEN XufengCUI Jiajia
(Northwest Institute of Nuclear Technology,Xi’an 710024,China)

AbstractBackground:In order to obtain high quality pulsed electron beam for large physical experimental facilities,it usually has higher requirement of the vibration of the foundation.Detailed geological vibration analysis is an important premise for the location of the experimental facilities.Purpose:Spectrum estimation method was applied to evaluate whether the location is reasonable through the actual vibration measurement for exploration area of different positions,different time-windows and different components.Methods:Firstly,the theories of traditional spectrum estimation were introduced.Then vibrations at different exploration positions,different time-windows,frequency bands and directions were analyzed according to the above methods,and amplitude spectrum and power spectrum of exploration area were given.Results &Conclusion:The data show that amplitude spectrum displacement of vertical component is 185 nm for exploration area which is 60 nm smaller than plane component,and give probability distribution of above result,providing the necessary basis for the mitigation measures to be taken in the process of further construction.

Key wordsLarge physical experimental facilities,Ground vibration,Spectrum estimation,Displacement spectrum,Amplitude spectrum

大型物理實(shí)驗(yàn)裝置（如大型加速器、激光器）為獲取高品質(zhì)的強(qiáng)流脈沖電子束流，通常對(duì)地基振動(dòng)都有較高的要求。它們?cè)诜燮陂g可能遭受到多種地基環(huán)境振動(dòng)，如重型卡車在附近駛過(guò)引起的地面振動(dòng)、靶室中靶的爆炸沖擊及其他爆轟試驗(yàn)產(chǎn)生的基礎(chǔ)振動(dòng)，以及高烈度地質(zhì)災(zāi)害引起的地基基礎(chǔ)振動(dòng)等，這些地基基礎(chǔ)環(huán)境振動(dòng)對(duì)裝置的正常運(yùn)行和束流品質(zhì)會(huì)造成重要影響[1-3]。

詳盡的地質(zhì)振動(dòng)狀況分析是裝置選址的重要前提。通過(guò)對(duì)裝置待勘區(qū)域的振動(dòng)測(cè)量，分析選址位置是否合理，以及進(jìn)一步對(duì)基建過(guò)程中應(yīng)采取的減震措施提供必要的依據(jù)?？睖y(cè)區(qū)域的場(chǎng)地條件包括巖土的物理力學(xué)性質(zhì)、動(dòng)力性能、土層結(jié)構(gòu)、覆蓋土的厚度，這些都會(huì)影響土層對(duì)于地震波的傳播。建筑物在場(chǎng)地地基之上，場(chǎng)地地基不僅承擔(dān)著建筑物的重量，又是傳遞地震能量的介質(zhì)，將地震波傳遞給上部結(jié)構(gòu)[4]。

目前國(guó)際上大型物理實(shí)驗(yàn)裝置在建設(shè)實(shí)驗(yàn)過(guò)程中，先期都開(kāi)展了一系列的振動(dòng)測(cè)量工作。位于西班牙巴塞羅那的新一代同步加速器輻射光源ALBA 在2005年項(xiàng)目實(shí)施前期，為評(píng)估場(chǎng)地的適應(yīng)性，曾在場(chǎng)地的不同區(qū)域采用寬頻帶地震計(jì)連續(xù)觀測(cè)一個(gè)月以上。位于德國(guó)漢堡的DESY(Deutsches Elektronen Synchrotron)機(jī)構(gòu)，在質(zhì)子加速器建設(shè)項(xiàng)目中，考慮振動(dòng)對(duì)線性對(duì)撞機(jī)的影響，先期充分測(cè)試了場(chǎng)地的振動(dòng)噪聲水平，在水冷和通風(fēng)系統(tǒng)設(shè)計(jì)和選址過(guò)程中利用了這一結(jié)果。在國(guó)內(nèi)，中國(guó)科學(xué)院上海應(yīng)用物理研究所在自由電子激光(Shanghai Free Electron Laser,SFEL)項(xiàng)目前期工程選址過(guò)程中，也曾在上海浦東張江園區(qū)內(nèi)的待勘區(qū)域中做了有關(guān)振動(dòng)測(cè)試方面的工作[5-7]。上海核工程研究設(shè)計(jì)院引進(jìn)的AP1000（國(guó)際標(biāo)準(zhǔn)化商業(yè)反應(yīng)堆）機(jī)組，它的結(jié)構(gòu)抗震要求很高，標(biāo)準(zhǔn)輸入值高達(dá)0.3 g（水平向及垂直向）[8]。文獻(xiàn)[9-10]是關(guān)于核電廠在選址過(guò)程中抗加速度等級(jí)和設(shè)計(jì)地震反應(yīng)譜的角度進(jìn)行的研究。

1　譜估計(jì)方法簡(jiǎn)介

信號(hào)的頻譜分析是研究信號(hào)特征的重要手段，對(duì)于確定性信號(hào)，可用Fourier變換來(lái)考察信號(hào)的頻譜特性，而對(duì)于廣義平穩(wěn)隨機(jī)信號(hào)而言，則需要求功率譜。功率譜不僅反映了隨機(jī)信號(hào)功率能量的分布特性，同時(shí)可以揭示信號(hào)中隱含的周期性以及相距很近的譜峰等有用信息，有著極為廣泛的應(yīng)用[11]。

譜估計(jì)方法通常分為兩類：非參數(shù)化方法和參數(shù)化方法。其中非參數(shù)化方法又稱經(jīng)典譜估計(jì)，如周期圖法、自相關(guān)法等；而參數(shù)化估計(jì)又稱現(xiàn)代譜估計(jì)，如AR模型法、平均移動(dòng)MA模型以及自回歸移動(dòng)平均模型法。

實(shí)際數(shù)字信號(hào)處理前，往往需要添加窗函數(shù)，常用的窗函數(shù)有矩形窗、Bartlett窗、Daniell窗、Hanming窗、Blackman窗等。不同窗函數(shù)的頻譜寬度和下降速度往往不盡相同，對(duì)于實(shí)際的譜估計(jì)會(huì)產(chǎn)生重要的影響。對(duì)于一定的數(shù)據(jù)長(zhǎng)度，周期圖估計(jì)首先取決于真實(shí)功率譜的形狀，其次是譜窗寬度的正確選擇。Barlett證明了平均周期圖等價(jià)于窗函數(shù)為三角形時(shí)的加窗周期圖（原始周期圖經(jīng)與譜窗卷積而得到的一種局部平均結(jié)果）?？紤]本文的實(shí)際應(yīng)用背景，為真實(shí)地反映待勘區(qū)實(shí)際的地基震動(dòng)情況，一般需要連續(xù)觀測(cè)數(shù)小時(shí)，涵蓋不同的工作時(shí)段，在基本的數(shù)據(jù)處理進(jìn)行之前將明顯具有人為振動(dòng)特征信號(hào)剔除掉（信噪比大于一定閾值），以便對(duì)實(shí)際的位移和速度盡可能做無(wú)偏估計(jì)。同時(shí)本文僅關(guān)注某一頻帶內(nèi)信號(hào)能量的平均值，對(duì)于具有明顯譜峰的信號(hào)往往不關(guān)心，一些現(xiàn)代譜估計(jì)方法在該應(yīng)用中不適用。因此本文采用經(jīng)典譜估計(jì)方法，窗函數(shù)選擇旁瓣較小的Hanming窗。

周期圖法是根據(jù)各態(tài)歷經(jīng)的隨機(jī)過(guò)程功率譜的定義進(jìn)行的譜估計(jì)。周期圖法是將隨機(jī)序列X(n)的N個(gè)觀測(cè)數(shù)據(jù)視為能量有限的序列，直接計(jì)算x(n)的離散傅里葉變換，得x(k)，然后再取其幅值的平方，并除以N，作為序列x(n)的真實(shí)功率譜的估計(jì)。

設(shè)平穩(wěn)隨機(jī)信號(hào)X(n)的有限個(gè)觀測(cè)值x(0)，x(1)，…，x(N-1)首先按式(1)求出其傅里葉變換，然后按式(2)進(jìn)行譜估計(jì)。

周期圖法應(yīng)用較為廣泛，主要是由于它與序列的頻譜有直接的對(duì)應(yīng)關(guān)系，且可以采用FFT快速算法來(lái)進(jìn)行計(jì)算。但該方法需要對(duì)無(wú)限長(zhǎng)的平穩(wěn)隨機(jī)序列進(jìn)行截?cái)?，相?dāng)于時(shí)域加矩形窗口，將會(huì)在頻域里使功率譜與窗函數(shù)卷積，從而產(chǎn)生頻譜泄露，容易使弱信號(hào)的主瓣被強(qiáng)信號(hào)的旁瓣所淹沒(méi)，造成頻譜的模糊和失真，使得譜分辨率降低，因此不滿足一致估計(jì)的條件。改進(jìn)的方法是將周期圖進(jìn)行平滑，使估計(jì)方差減小，從而得到一致譜估計(jì)。

Bartlett平均周期圖的方法是將N點(diǎn)的有限長(zhǎng)序列x(n)分段求周期圖再平均。將長(zhǎng)度為N的數(shù)據(jù)分為L(zhǎng)段，每段長(zhǎng)度為M。先對(duì)每段數(shù)據(jù)用周期圖法進(jìn)行譜估計(jì)，然后對(duì)L段求平均得到長(zhǎng)度為N的數(shù)據(jù)功率譜。Welch譜估計(jì)法是在Bartlett法基礎(chǔ)上進(jìn)行改進(jìn)，目的是在保持Bartlett法方差性能的同時(shí)，改善其分辨率。其基本原理是先對(duì)隨機(jī)序列分段，使每一段有部分重疊，然后對(duì)每一段數(shù)據(jù)用一個(gè)合適的窗函數(shù)進(jìn)行平滑處理，并對(duì)各段譜求平均，最后得到功率譜如式(3)，該估計(jì)的分辨率為式(4)所示[12]。

由于Welch譜估計(jì)法允許各段數(shù)據(jù)交疊，所以數(shù)據(jù)段L會(huì)增加，使方差得到很大的改善，但是數(shù)據(jù)的交疊又減少了每一段數(shù)據(jù)的不相關(guān)性，使方差的減少不會(huì)達(dá)到理論程度。同時(shí)由于選擇了合適的窗函數(shù)可以有效地改善信號(hào)的頻譜泄露，增加譜峰寬度，進(jìn)而提高分辨率。上述方法主要用于信號(hào)的功率譜估計(jì)。力傳感器按輸出信號(hào)類型的不同可分為速度型、加速度型和位移型，若x(n)為數(shù)字化后的速度信號(hào)，按式(5)即可得到對(duì)應(yīng)的速度振幅譜是S(ω)。其中運(yùn)算代表取模。由于速度是位移的導(dǎo)數(shù)，則位移譜只需將S(ω)除以jw因子即可。

2　地基振動(dòng)測(cè)量方法

通常完整的地基測(cè)量過(guò)程分為基本測(cè)點(diǎn)選擇和測(cè)試系統(tǒng)搭建。測(cè)點(diǎn)選擇首先考慮待勘區(qū)域內(nèi)的巖土性質(zhì)、周圍地理人文環(huán)境以及有關(guān)氣象資料。測(cè)試系統(tǒng)的搭建主要依據(jù)地震行業(yè)中振動(dòng)測(cè)量的典型方法，根據(jù)實(shí)際的測(cè)量要求，需要選擇相應(yīng)的傳感器。本文根據(jù)大型實(shí)驗(yàn)裝置的工作要求，僅關(guān)注地基振動(dòng)的位移和速度。

2.1振動(dòng)測(cè)試系統(tǒng)介紹

測(cè)試系統(tǒng)由地震數(shù)據(jù)采集器、地震傳感器、GPS時(shí)間系統(tǒng)和便攜式筆記本組成。其中地震計(jì)設(shè)備為德國(guó)的3Dlite地震計(jì)，該型地震計(jì)是德國(guó)Lennartz Electronic公司生產(chǎn)的一款輕便型負(fù)反饋式速率地震計(jì)，具有較高的靈敏度、較低的儀器噪聲和較為輕便的體積。其靈敏度為400 V·m-1·s，帶寬為1-80 Hz儀器噪聲水平小于3 nm·s-1(rms@1 Hz)。德國(guó)漢堡的DESY機(jī)構(gòu)也主要基于速度型檢波器進(jìn)行相應(yīng)的地基振動(dòng)測(cè)量。數(shù)據(jù)采集設(shè)備為美國(guó)REFTEK-130，該采集器的模數(shù)轉(zhuǎn)換位數(shù)為24位，并且在50 Hz采樣率條件下的噪聲水平RMS約為1 計(jì)數(shù)。該儀器設(shè)備是目前地震行業(yè)廣泛應(yīng)用的觀測(cè)設(shè)備，具有較高的采集精度和運(yùn)行穩(wěn)定性。外部電源采用直流免維護(hù)電池，避免引入工頻干擾。圖1(a)為速度型地震計(jì)原理示意圖，圖1(b)為振動(dòng)測(cè)試系統(tǒng)硬件連接配置圖[13-14]。

2.2振動(dòng)測(cè)點(diǎn)分布

測(cè)試過(guò)程主要分兩個(gè)階段：首先在西北核技術(shù)研究所反應(yīng)堆實(shí)驗(yàn)大廳內(nèi)部和待勘區(qū)域分別放置測(cè)點(diǎn)進(jìn)行觀測(cè)，該大廳在基建過(guò)程中曾采用了隔離減震措施；第二階段在大型實(shí)驗(yàn)裝置的待勘區(qū)域內(nèi)，選擇合適的位置沿地表按照?qǐng)D2斜下方進(jìn)行開(kāi)挖，分別在1 m、2 m、3 m、4 m處平整出1 m×1 m方形平臺(tái)分別放置測(cè)點(diǎn)進(jìn)行觀測(cè)。測(cè)點(diǎn)的儀器連接配置按§2.1所示。開(kāi)挖的俯視和側(cè)視圖如圖2所示。

圖2　待勘選區(qū)域測(cè)點(diǎn)開(kāi)挖剖面(a)與俯視圖(b)Fig.2　Section plane(a)and planform(b)of test station of reserved zone.

3　數(shù)據(jù)處理與結(jié)果分析

地震計(jì)輸出的電壓信號(hào)通過(guò)數(shù)據(jù)采集器模數(shù)轉(zhuǎn)換后成為24位數(shù)字信號(hào)，該轉(zhuǎn)換系數(shù)η為1.589mV·count-1，地震計(jì)靈敏度S為400 V·m-1·s。通過(guò)簡(jiǎn)單的線性關(guān)系可以推算出瞬時(shí)振動(dòng)速度量和位移量，其中速度為C/(S/η)，位移為C/(2πf·S/η)，代表實(shí)際的數(shù)字轉(zhuǎn)換量。累計(jì)采集一定時(shí)間內(nèi)的數(shù)據(jù)，即可以得到相應(yīng)時(shí)域變化v(t)和u(t)。在地基振動(dòng)測(cè)量當(dāng)中，我們主要關(guān)心是某時(shí)間段內(nèi)＞1 Hz的振動(dòng)位移值以及＞8 Hz的振動(dòng)速率值，均使用均方根值來(lái)描述：

式中：N代表實(shí)際采樣點(diǎn)數(shù)。所有測(cè)點(diǎn)的數(shù)據(jù)采集器采集參數(shù)配置為200 Hz采樣率、線性相位濾波器，無(wú)高通濾波器，前置放大器增益因子為1。濾波器采用巴特沃斯帶通濾波。

為全面了解在一天中不同時(shí)間段地基振動(dòng)的變化程度，各個(gè)測(cè)點(diǎn)連續(xù)測(cè)量72 h。數(shù)據(jù)回收后采用西北核技術(shù)研究所自研的地震波數(shù)據(jù)分析處理軟件SWA(Seismic Waveform Analysis)進(jìn)行處理，該軟件集成了地震數(shù)據(jù)功率譜估計(jì)（包括位移功率譜、速度功率譜、加速度功率譜）等多種譜估計(jì)的方法。通過(guò)對(duì)記錄后的數(shù)據(jù)頻域處理，得到不同時(shí)間段、不同方向的地基振動(dòng)情況。

圖3給出了不同方向各觀測(cè)點(diǎn)全時(shí)間段的位移功率譜統(tǒng)計(jì)結(jié)果?？梢钥吹椒磻?yīng)堆內(nèi)外兩個(gè)測(cè)點(diǎn)的本底位移基本一致，室外測(cè)點(diǎn)（待選區(qū)域）相對(duì)室內(nèi)兩個(gè)測(cè)點(diǎn)來(lái)說(shuō)本底位移偏高。三個(gè)方向在3 Hz存在平均位移峰值約為1μm。由圖3還可以看出，在20 Hz、40 Hz附近存在譜峰，可能是周圍建筑施工的影響。其中，CS2為反應(yīng)堆大廳減震槽內(nèi)測(cè)點(diǎn)，F(xiàn)Y3為反應(yīng)堆大廳減震槽外測(cè)點(diǎn)，F(xiàn)Y2為反應(yīng)堆大廳外測(cè)點(diǎn)。

圖3　不同方向各觀測(cè)點(diǎn)全時(shí)間段的位移功率統(tǒng)計(jì)結(jié)果Fig.3　Statistical results of the power spectra of different components to the full time period of station.

圖4、5按照不同時(shí)間段及不同方向分別統(tǒng)計(jì)了其振幅位移譜和速度譜，其反映了在不同頻率上的位移及速度分布情況；時(shí)域校正波形給出了地表的實(shí)際速度波動(dòng)及其最大值，另外給出了表示速度波動(dòng)的分布概率和概率分布曲線。不同時(shí)間段的統(tǒng)計(jì)結(jié)果顯示，室內(nèi)測(cè)點(diǎn)受人為環(huán)境噪聲影響較小，室外測(cè)點(diǎn)的本底位移受人為環(huán)境噪聲影響明顯。從振幅位移譜上可以看到，對(duì)于室內(nèi)測(cè)點(diǎn)在大于1 Hz的范圍內(nèi)位移峰值約在800 nm。垂向速度相對(duì)于水平向速度偏高，對(duì)于室內(nèi)測(cè)點(diǎn)在8 Hz以上小于4μm·s-1，水平向在8 Hz以上小于2μm·s-1。數(shù)據(jù)時(shí)間窗選擇凌晨零點(diǎn)至3點(diǎn)和早晨8點(diǎn)至11點(diǎn)，圖4、5分別給出了歸一化時(shí)域波形，及該時(shí)間段內(nèi)的速度分布概率。

圖4　北京時(shí)間凌晨0:00-3:00點(diǎn)CS2與FY3垂直方向比較(a)兩個(gè)臺(tái)站平均位移振幅譜，(b)兩個(gè)臺(tái)站平均速度振幅譜，(c)CS2臺(tái)站垂直通道原始波形，(d)CS2臺(tái)站去除儀器響應(yīng)后的波形圖，(e)CS2臺(tái)站速度分布概率，(f)CS2臺(tái)站速度累積概率曲線，(g)FY3臺(tái)站垂直通道原始波形，(h)FY3臺(tái)站去除儀器響應(yīng)后的波形圖，(i)FY3臺(tái)站速度分布概率，(j)FY3臺(tái)站速度累積概率曲線Fig.4　Vertical component at 0:00-3:00 a.m.between CS2 and FY3.(a)Mean displacement ASD,(b)Mean velocity ASD,(c)Original waveform of CS2,(d)Corrected velocity,(e)Velocity distribution probability of CS2,(f)Velocity cumulative distribution probability curve of CS2,(g)Original waveform of FY3,(h)Corrected velocity,(i)Velocity distribution probability of FY3,(j)Velocity cumulative distribution probability curve of FY3

圖5　北京時(shí)間上午8:00-11:00點(diǎn)CS2、FY2與FY3垂直方向比較Fig.5　Vertical component at 8:00-11:00 a.m.between CS2,FY2 and FY3.

圖5中，(a)三個(gè)臺(tái)站平均位移振幅譜，(b)三個(gè)臺(tái)站平均速度振幅譜，(c)CS2臺(tái)站垂直通道原始波形，(d)CS2臺(tái)站去除儀器響應(yīng)后的波形圖，(e)CS2臺(tái)站速度分布概率，(f)CS2臺(tái)站速度累積概率曲線，(g)FY2臺(tái)站垂直通道原始波形，(h)FY2臺(tái)站去除儀器響應(yīng)后的波形圖，(i)FY2臺(tái)站速度分布概率，(j)FY2臺(tái)站速度累積概率曲線，(k)FY3臺(tái)站垂直通道原始波形，(l)FY3臺(tái)站去除儀器響應(yīng)后的波形圖，(m)FY3臺(tái)站速度分布概率，(n)FY3臺(tái)站速度累積概率曲線。

在裝置待選位置當(dāng)前地基條件下（介質(zhì)為黃土），在工作時(shí)間8:00-18:00，大于1 Hz振動(dòng)位移在220 nm以內(nèi)，大于8 Hz振動(dòng)速率在4.1μm·s-1。反應(yīng)堆場(chǎng)房?jī)?nèi)地基為經(jīng)過(guò)處理的地基，測(cè)量結(jié)果比選址位置的位移噪聲小，大于1 Hz振動(dòng)位移小于183 nm，大于8 Hz振動(dòng)速率小于1.5μm·s-1。凌晨時(shí)，選址位置大于8 Hz振動(dòng)速率為1.2μm·s-1。

通過(guò)對(duì)測(cè)量數(shù)據(jù)的分析，得到的測(cè)點(diǎn)位置不同時(shí)間段（上午及凌晨）、不同方向振動(dòng)位移與振動(dòng)速率結(jié)果如表1所示。

表1　北京時(shí)間2014年6月12日08:00-11:00和00:00-03:00地震測(cè)量數(shù)據(jù)(RMS)對(duì)比Table 1　Comparison to seismic data(RMS)of different zones at June 12,2014 08:00-11:00 and 00:00-03:00.

圖6是按照第二階段測(cè)量計(jì)劃，在待選區(qū)域不同深度放置的測(cè)點(diǎn)，分別選擇凌晨和下午工作時(shí)間。地表附近測(cè)點(diǎn)T04與其他三個(gè)測(cè)點(diǎn)平均位移功率譜具有明顯差異，約為10倍。而其他三個(gè)測(cè)點(diǎn)在10 Hz以下基本一致，20 Hz以上平均位移幅度譜隨深度的變化有減小趨勢(shì)。凌晨T01、T02、T03測(cè)點(diǎn)在6 Hz附近存在位移幅度譜峰值約為14 nm2·Hz-1，而在白天（北京時(shí)間14:00-15:00）在6 Hz附近的峰值則會(huì)升高至63 nm2·Hz-1。

圖6　待勘區(qū)域不同測(cè)點(diǎn)相同時(shí)間段垂直方向位移幅度譜比較Fig.6　Comparison of displacement power spectra of different test zones at the same time in vertical component.

4　結(jié)語(yǔ)

本文對(duì)大型實(shí)驗(yàn)裝置中待勘區(qū)域的振動(dòng)進(jìn)行測(cè)試，對(duì)結(jié)果進(jìn)行不同位置、不同時(shí)刻、不同頻段、不同方向的分析。

1)待勘區(qū)域地表附近處在垂直方向1-80 Hz頻段內(nèi)，連續(xù)記錄24 h時(shí)間數(shù)據(jù)，其均方根位移約185nm。

2)人文背景噪聲振動(dòng)主要集中在低頻部分，一定的地基減震措施對(duì)8 Hz以上頻段起到消減左右，而該部分振動(dòng)卻是地基振動(dòng)的主要構(gòu)成部分。

3)垂直方向振幅大于水平方向。通常在安靜環(huán)境下，水平方向的振幅要大于垂直方向的振幅，本文提到的待勘區(qū)域由于距離高速公路較近，平時(shí)車流量較多，因此垂直方向是大型實(shí)驗(yàn)裝置后期工程施工需要考慮的主要問(wèn)題。

4)在未采取任何減震措施條件下不同深度20 Hz以下振動(dòng)的位移振幅譜影響不明顯。

參考文獻(xiàn)

1 姜巖松,王毅,蘇寶庫(kù).MLKF譜估計(jì)在地基噪聲特性分析中的應(yīng)用[J].電機(jī)與控制學(xué)報(bào),2011,15(2):103-104 JIANG Yansong,WANG Yi,SU Baoku.Application on ground noise analysis by MLKF spectral estimation[J].Electric Machines and Control,2011,15(2):103-104

2劉建波,李洪,陳楠,等.地基振動(dòng)對(duì)大型加速器對(duì)中的影響[J].激光與粒子束,2010,22(4):817-818 LIU Jianbo,LI Hong,CHEN Nan,et al.Alignment for large-scale accelerator after ground vibration[J].High Power Laser and Particle Beams,2010,22(4):817-818

3馮忠俊,王曉,杜涵文,等.上海自由電子激光工程地基振動(dòng)測(cè)試研究[J].核技術(shù),2008,31(4):245-249 FENG Zhongjun,WANG Xiao,DU Hanwen,et al.Ground vibration measurement on Shanghai free electron laser project[J].Nuclear Techniques,2008,31(4):245-249

4毛尚禮,余湘娟,張富有.地基隔震減震機(jī)理研究[J].華南地震,2010,30(3):75-76 MAO Shangli,YU Xiangjuan,ZHANG Fuyou.Studies on mechanism of ground’s seismic isolation and shock absorption[J].South China Journal of Seismology,2010,30(3):75-76

5Wilhelm Bialowons,HeiKo Ehrlichmann.Ground vibration measurements at ALBA-Barcelon[OL].Presentations in the XFEL LINAC Module meetings,ILC project meeting &other meetings:1-12[2005-09-06].http://vibration.desy.de/documents/presentations/[2-13]

6Wilhelm Bialowons,HeiKo Ehrlichmann.Ground vibration measurements with seismic sensors[OL].Workshop on Ambient Ground Motion and Civil Engineering for Low Emittance:1-4[2005-07-21].http://vibration.desy.de/documents/presentations/

7Juravlev V M.Investigations of power and spatial correlation characteristics of seismic vibrations in the CERN LEP tunnel for linear collider studies[R].CERN-SL/93-53,CLIC-Note 217,1993

8侯春林,李小軍,潘蓉,等.AP1000設(shè)計(jì)地震反應(yīng)譜在具體廠址評(píng)價(jià)中的應(yīng)用[J].原子能科學(xué)技術(shù),2013,47(7):1196-1197 HOU Chunlin,LI Xiaojun,PAN Rong,et al.Application of AP1000 certified seismic design response spectra in site specific evaluation[J].Atomic Energy Science and Technology,2013,47(7):1196-1197

9李琦,杜建勇,柳琳琳,等.AP1000主控室設(shè)備抗震鑒定試驗(yàn)研究[J].核技術(shù),2013,36(4):040614.DOI:10.11889/j.0253-3219.2013.hjs.36.040614 LI Qi,DU Jianyong,LIU Linlin,et al.Experimental study on seismic qualification test of AP1000 main control room equipment[J].Nuclear Techniques,2013,36(4):040614.DOI:10.11889/j.0253-3219.2013.hjs.36.040614

10 馬淵睿,朱翊洲,王赤虎,等.AP1000 核電站主控室盤臺(tái)抗震鑒定試驗(yàn)研究和工程實(shí)踐[J].核技術(shù),2013,36(4):040611.DOI:10.11889/j.0253-3219.2013.hjs.36.040611 MA Yuanrui,ZHU Yizhou,WANG Chihu,et al.Research and project practice in AP1000 nuclear plant MCR assemblies seismic test[J].Nuclear Techniques,2013,36(4):040611.DOI:10.11889/j.0253-3219.2013.hjs.36.040611

11 金連文,韋剛.現(xiàn)代數(shù)字信號(hào)處理簡(jiǎn)明教程[M].北京:清華大學(xué)出版社,2004:123-156 JIN Lianwen,WEI Gang.Advanced digital signal processing[M].Beijing:Tsinghua University Press,2004:123-156

12 王春興.基于MATLAB實(shí)現(xiàn)經(jīng)典功率譜估計(jì)[J].曲阜師范大學(xué)學(xué)報(bào),2011,37(2):59-61 WANG Chunxing.Achieve power spectrum estimation based on MATLAB[J].Journal of Qufu Normal University,2011,37(2):59-61

13 張曉鵬,李世維,王同東,等.基于檢波器的短周期地震計(jì)的改進(jìn)[J].地球物理學(xué)進(jìn)展,2014,29(5):2466-2471.DOI:10.6038/pg20140572 ZHANG Xiaopeng,LI Shiwei,WANG Tongdong,et al.Animprovedshortperiodgeophone-based seismometer[J].Progress in Geophysics,2014,29(5):2466-2471.DOI:10.6038/pg20140572

14 彼得·鮑曼.中國(guó)地震局監(jiān)測(cè)預(yù)報(bào)司(譯).新地震觀測(cè)實(shí)踐手冊(cè)(第1卷)[M].北京:地震出版社,2006:150-171 Peter Bormanm.Monitoring and Forecasting Department of China Earthquake Administration(translation).New manual of seismological observatory practice(Vol.1)[M].Beijing:Earthquake Press,2006:150-171

收稿日期：2015-08-19，修回日期：2015-11-04

DOI:10.11889/j.0253-3219.2016.hjs.39.030503

中圖分類號(hào)TL505，TU47