發(fā)布時間：2021-07-21 11:35

　　目前人們在進行地震動相干函數模型參數的擬合時,選取的頻段范圍存在很大差異,而確定有效的頻段范圍對相干函數模型擬合結果的合理性很重要。本文利用AR模型（Autoregressive Model,自回歸模型）計算地震動的功率譜,提出了根據功率譜"信噪比（SNR）"確定相干函數有效頻段范圍的方法。本文的信噪比中的"噪",專指AR模型中的激勵白噪聲,用于衡量AR模型中某一頻率地震動功率譜與激勵白噪聲功率譜之間的比例。定義信噪比SNR=0時對應的頻率為有效頻段的最高頻率,大于最高頻率的頻率成分,被當作是地震波中的隨機成分,可以不考慮地震動的相關性。采用這一方法,對SMART-1第45次地震記錄的不同臺站間距下的地震波進行了功率譜估計,根據信噪比確定出了有效頻段范圍,并計算了相應的相干函數,結果表明此方法是可行的。 

【文章來源】：地震工程與工程振動. 2020,40(01)北大核心CSCD

【文章頁數】：9 頁

【部分圖文】：

地震動記錄數據的組成

本文選取SMART-1臺陣的第45次地震的地震數據作為研究對象，并以第45次地震東西(EW)方向的內圈5個臺站(C00、I03、I06、I09、I12)的S波數據為例，簡述AR模型的建模過程。首先對加速度記錄積分得到速度記錄，如圖2。由速度記錄建立多維AR模型，階數(即公式(1)中的m)范圍選為2～15，采用FPE準則定階，如圖3(a)，取FPE值最小時的階次為適用模型的階次。因為2～5階的FPE值過大，圖中并未畫出，可直接忽略，在6～15階的范圍內，階次為10時的FPE值為最小，所以模型階次初選為10。除了考慮FPE值，還需要檢驗此階次下的模型的殘差自相關系數。如果殘差自相關系數落在以0為均值的99%置信區(qū)間內，則說明殘差之間是互不相關的，進而說明模型殘差符合白噪聲假定，模型滿足要求。圖3(b)為10階AR模型的殘差自相關系數，圖中灰色區(qū)域為99%置信區(qū)間，可見殘差自相關系數都在此置信區(qū)間內，此時模型滿足要求，模型的階次即定為10階。由最小二乘法求得在此階數下的模型參數Φ，參數的個數與模型階數相同，每個參數矩陣的維數與計算所使用的時程條數相同，因此，可以計算得到10個模型參數，每個參數為五維方陣，見式(11)。

由速度記錄建立多維AR模型，階數(即公式(1)中的m)范圍選為2～15，采用FPE準則定階，如圖3(a)，取FPE值最小時的階次為適用模型的階次。因為2～5階的FPE值過大，圖中并未畫出，可直接忽略，在6～15階的范圍內，階次為10時的FPE值為最小，所以模型階次初選為10。除了考慮FPE值，還需要檢驗此階次下的模型的殘差自相關系數。如果殘差自相關系數落在以0為均值的99%置信區(qū)間內，則說明殘差之間是互不相關的，進而說明模型殘差符合白噪聲假定，模型滿足要求。圖3(b)為10階AR模型的殘差自相關系數，圖中灰色區(qū)域為99%置信區(qū)間，可見殘差自相關系數都在此置信區(qū)間內，此時模型滿足要求，模型的階次即定為10階。由最小二乘法求得在此階數下的模型參數Φ，參數的個數與模型階數相同，每個參數矩陣的維數與計算所使用的時程條數相同，因此，可以計算得到10個模型參數，每個參數為五維方陣，見式(11)。2.2 功率譜估計與相干函數計算

【參考文獻】：
期刊論文
[1]基于能量的相干函數有效頻段范圍的選取[J]. 丁海平,朱越,于彥彥.  地震工程與工程振動. 2019(02)
[2]截止頻率的取值對地震動空間相干函數統(tǒng)計結果的影響[J]. 丁海平,羅翼,饒威波,朱越.  地震學報. 2018(05)
[3]窗函數和帶寬對地震動相干函數的影響[J]. 丁海平,袁莉莉,劉成浩.  自然災害學報. 2018(03)
[4]地震動的空間變化特性分析與修正相干模型[J]. 李英民,吳哲騫,陳輝國.  振動與沖擊. 2013(02)



本文編號：3294948