本站小編為你精心準(zhǔn)備了軟件開發(fā)及強震數(shù)據(jù)分析參考范文，愿這些范文能點燃您思維的火花，激發(fā)您的寫作靈感。歡迎深入閱讀并收藏。

摘要

利用Matlab平臺，設(shè)計一款用于強震數(shù)據(jù)處理的Wave-site軟件，與在用軟件相比，增加臺站場地響應(yīng)計算等功能，并提供多種基線校準(zhǔn)和場地響應(yīng)計算方法作比較分析。應(yīng)用該軟件研究了三山強震臺場地響應(yīng)與震級的關(guān)系，數(shù)據(jù)表明，地震S波放大因子大于地脈動，且隨震級升高而增大。

關(guān)鍵詞

Matlab；Wave-site軟件；強震動記錄；基線校準(zhǔn)；場地響應(yīng)

強震研究不僅為地震區(qū)劃、抗震設(shè)防和應(yīng)急救災(zāi)提供科學(xué)依據(jù)，也在地震動衰減關(guān)系、反應(yīng)譜特征參數(shù)統(tǒng)計等地震工程研究領(lǐng)域發(fā)揮著重要作用（詹小艷，2014）。理論研究和震害調(diào)查表明，場地條件對于地震動特性和震害評估具有重要意義（Hiroshi，2004），利用強震觀測資料分析強震臺站場地的放大效應(yīng)十分必要（Nakamura，1989）。安徽強震動臺網(wǎng)（包含9個臺站）自2007年1月運行以來，獲取多條強震記錄。隨著“十二五”建設(shè)，安徽強震臺站數(shù)量將劇增至61個，均采用實時網(wǎng)絡(luò)通信。然而，強震分析軟件功能比較簡單，無法實現(xiàn)各種情況下記錄曲線的基線校準(zhǔn)（鄭水明，2010），也無法進行場地響應(yīng)計算（David，2002）。目前，已有使用Matlab平臺進行強震數(shù)據(jù)處理的做法（曹加良，2009），但可視化操作不強，難以推廣應(yīng)用。利用Matlab平臺開發(fā)設(shè)計Wave-site軟件，用于強震記錄的基線校準(zhǔn)和臺站的場地響應(yīng)研究，具備傅氏譜計算和數(shù)據(jù)輸出功能，為強震數(shù)據(jù)的廣泛應(yīng)用做準(zhǔn)備。

1系統(tǒng)設(shè)計

Matlab平臺擁有大量工具箱和豐富的庫函數(shù)，適于數(shù)據(jù)處理。Matlab可視化功能強大，方便制作各種繪圖。同時，MatlabGUI界面控件較多，程序編寫簡單，軟件可編譯為獨立于Matlab環(huán)境運行的C程序（安裝MCR環(huán)境）。Wave-site軟件通過GUI界面進行設(shè)計（圖1），采用Matlab內(nèi)部算法進行數(shù)據(jù)處理，軟件參數(shù)配置見表1。

2系統(tǒng)功能

2.1基線校準(zhǔn)強震加速度記錄基線校準(zhǔn)方法有基線法、Iwan法和低頻濾波法。基線法是從加速度記錄中減去震前地脈動平均值，無震前記錄則減去整個記錄的平均值。對復(fù)雜的基線漂移進行校準(zhǔn)，其結(jié)果可作為進一步處理的基礎(chǔ)。當(dāng)峰值加速度超過50cm/s2時，可以采用Iwan法。該方法基本思路是，將加速度時程分為3段，中間一段代表由磁滯效應(yīng)引起的基線偏移發(fā)生在強震部分，且認(rèn)為在同一段內(nèi)基線偏移是一個常量，用速度時程基線斜率表示加速度基線偏移量。使用濾波方法，能有效清除低頻噪音對加速度基線的影響，但加速度信息會被引入新的噪音信號，部分長周期信息和地面永久位移信息也可能被濾掉。實際工作中，為盡可能消除記錄中的低頻噪音，可以將3種方法相互結(jié)合使用。圖2為分別采用基線法（左圖）、Iwan法（右圖）得到的加速度曲線圖。

2.2傅氏譜計算應(yīng)用Matlab內(nèi)部離散FFT計算方法，對校準(zhǔn)后的加速度曲線進行傅氏譜計算，并在軟件主界面的坐標(biāo)圖內(nèi)顯示出來。

2.3場地響應(yīng)計算場地響應(yīng)計算方法有H/V譜比法、Atkinson法和Moya法，其中強震研究常用H/V譜比法。Nakamur（a1989）利用井下記錄數(shù)據(jù)證實HB/VB≈1，其中HB、VB分別代表基巖場地（下角標(biāo)為Bedrock縮寫）觀測資料的水平向和垂直向傅氏頻譜。而Lermo（1993）利用基巖為參考點也證實此值近似于1。因此，場地響應(yīng)可表示為此為H/V譜比法經(jīng)驗關(guān)系式，其中HS和VS為土層場地（下角標(biāo)為Soil縮寫）觀測資料的水平向和垂直向傅氏頻譜。公式（1）表明，只需在地表測量3個分量就可計算出場地傳遞函數(shù)。在實際研究中，通常采用H/V譜比法計算臺站場地響應(yīng)。

2.4其他功能Wave-site軟件支持?jǐn)?shù)據(jù)輸出，包括傅氏譜數(shù)據(jù)和場地響應(yīng)數(shù)據(jù)，且使用者可以根據(jù)需要不斷完善和擴充該軟件的功能。

3應(yīng)用實例

應(yīng)用Wave-site軟件，分別計算三山臺兩次強震記錄及震前地脈動的場地響應(yīng)。具體做法是：截取基線校準(zhǔn)后的強震記錄S波數(shù)據(jù)及地脈動數(shù)據(jù)（時間長度均為15s），應(yīng)用Matlab內(nèi)置FFT算法得到傅氏譜。為避免譜比失穩(wěn)，對各分量傅氏譜幅值采用0.1Hz平滑窗，即200個采樣點的平均進行平滑。采用origin7.0作圖，計算結(jié)果見圖3。在此將譜比值描述為場地的放大因子，從圖3可以看出，三山臺的譜比曲線具有較好的相似性。結(jié)合三山臺場地條件，對比地脈動與S波譜比曲線，可以得出以下結(jié)論。（1）三山臺地脈動和S波譜比曲線一致性好，在0.5—5s頻段，場地響應(yīng)相對較平坦。（2）由于覆蓋層較厚，地脈動和S波譜比曲線相似性較好，曲線趨向一致，但場地的放大因子差距較大（表2）。地震S波放大因子大于地脈動，且隨著震級升高而增大，可見土層場地對于地震的放大效應(yīng)較明顯。

4結(jié)束語

利用Matlab平臺開發(fā)的Wave-site軟件，具有功能全面、可持續(xù)升級等特點，適用于強震數(shù)據(jù)處理，具有推廣價值。利用該軟件計算安徽地區(qū)強震臺站場地響應(yīng)數(shù)據(jù)，并作為一項基礎(chǔ)資料加入強震臺站建臺報告。近期考慮將此軟件計算結(jié)果與手動計算結(jié)果進行對比，對軟件計算過程進行優(yōu)化，以期獲得更準(zhǔn)確的場地響應(yīng)數(shù)據(jù)。

參考文獻

曹加良.用Matlab軟件內(nèi)部函數(shù)校正強震記錄的探討[J].華南地震，2009，29（2）：88-94.

詹小艷，等.2012年7月20日江蘇寶應(yīng)M4.9地震加速度場的計算[J].地震地磁觀測與研究，2014，35（3/4）：47-50.

鄭水明，等.強震加速度記錄基線校正問題探討[J].大地測量與地球動力學(xué)，2010，30（3）：47-50.

作者：趙希磊 張學(xué)應(yīng) 龍劍鋒 李發(fā) 單位：安徽省地震局