[發(fā)明專利]基于隨機傳遞函數(shù)法的泄洪誘發(fā)場地振動預(yù)測方法

在審

• 申請?zhí)枺?/td>: 201710548624.0
• 申請日: 2017-07-06
• 公開（公告）號: CN107103398A
• 公開（公告）日: 2017-08-29
• 發(fā)明（設(shè)計）人: 張龑;張國新;李松輝;劉毅;范哲
• 申請（專利權(quán)）人: 中國水利水電科學(xué)研究院
• 主分類號: G06Q10/04
• 分類號: G06Q10/04;G06K9/00
• 代理公司: 北京北新智誠知識產(chǎn)權(quán)代理有限公司11100
• 代理人: 陳英

1.基于隨機傳遞函數(shù)法的泄洪誘發(fā)場地振動預(yù)測方法，其特征在于，包括：

對每個單一振源的實測信號進行濾波處理后，組成多振源激勵源信號，作為輸入信號，

對輸出測點的實測信號進行濾波處理后，作為輸出響應(yīng)信號，

基于該輸入信號與輸出響應(yīng)信號，建立傳遞函數(shù)，基于傳遞函數(shù)的無偏估計法，計算預(yù)測結(jié)果信號，

在預(yù)測結(jié)果信號的基礎(chǔ)上疊加噪聲序列，作為最終的預(yù)測結(jié)果。

2.根據(jù)權(quán)利要求1所述的基于隨機傳遞函數(shù)法的泄洪誘發(fā)場地振動預(yù)測方法，其特征在于，所述對實測信號進行濾波處理的方法為：

S1：構(gòu)造與實測信號的信號長度相對應(yīng)的白噪聲信號，對其進行EEMD分解，計算噪聲指標(biāo)值η_j；

η_j＝σ_j/σ₁(3)

其中，σ_j為第j個IMF分量的標(biāo)準差，σ₁為第1個IMF分量的標(biāo)準差；

σj=Σk=1N[cj(k)-cj‾]2N-1---(1)]]>

其中，c_j(k)為第j個IMF分量，為c_j(k)的均值，N為信號長度；

S2：對實測信號進行EEMD分解，計算其各IMF分量的標(biāo)準差σ_j’；

S3：計算經(jīng)EEMD分解的實測信號的各IMF分量的噪聲標(biāo)準差λ_j；

噪聲標(biāo)準差λ_j的計算公式為：

λ_j＝η_jσ₁`(4)

其中，j＝2，3，4，...，n，n為整數(shù)；

S4：判別實測信號各IMF分量中包含的噪聲成分，當(dāng)λ_j等于或大于相應(yīng)的σ_j’時，將該第j個IMF分量直接濾除；當(dāng)λ_j小于相應(yīng)的σ_j’時，對該第j個IMF分量進行小波閾值濾波。

3.根據(jù)權(quán)利要求2所述的基于隨機傳遞函數(shù)法的泄洪誘發(fā)場地振動預(yù)測方法，其特征在于，所述小波閾值濾波的計算公式為：

T=λj2lnN---(6)]]>

4.根據(jù)權(quán)利要求1所述的基于隨機傳遞函數(shù)法的泄洪誘發(fā)場地振動預(yù)測方法，其特征在于，所述噪聲序列的方差為：

σx-y2=σy2e-SNR20---(8)]]>

其中，x為濾波處理前的輸出響應(yīng)信號，y為濾波處理后的輸出響應(yīng)信號，為濾波處理后時間序列的方差，為濾波處理后的噪聲序列的方差。