基于BP神經(jīng)網(wǎng)絡(luò)的場(chǎng)地等效剪切波速變化預(yù)測(cè)研究

摘要：利用日本KiK-net臺(tái)網(wǎng)提供的407個(gè)臺(tái)站的30 952條地震動(dòng)記錄，提出了一種基于BP神經(jīng)網(wǎng)絡(luò)的場(chǎng)地等效剪切波速比變化預(yù)測(cè)模型。模型采用了均方誤差函數(shù)及Adam優(yōu)化算法，由3個(gè)輸入?yún)?shù)、5個(gè)隱藏神經(jīng)元及1個(gè)輸出參數(shù)組成。輸入?yún)?shù)為地面峰值加速度PGA、Arias烈度I_a及場(chǎng)地剪切波速V_S30，輸出為場(chǎng)地等效剪切波速比（V_Sr）。研究結(jié)果表明：該神經(jīng)網(wǎng)絡(luò)模型殘差對(duì)于各輸入變量整體呈現(xiàn)出無(wú)偏的特性，在大部分的軟硬場(chǎng)地中均有較好的預(yù)測(cè)性能，該模型預(yù)測(cè)得到的PGA關(guān)于V_S30的相關(guān)系數(shù)曲線與用傳統(tǒng)的最小二乘法回歸得到的函數(shù)曲線相比，其相關(guān)系數(shù)有更好的表現(xiàn)。該模型預(yù)測(cè)曲線顯示，B類場(chǎng)地在PGA達(dá)到175 cm/s²時(shí)，場(chǎng)地剪切波速下降5%，D、E類場(chǎng)地在PGA達(dá)到140 cm/s²時(shí)，場(chǎng)地剪切波速下降5%，多數(shù)場(chǎng)地的非線性閾值為50～100 cm/s²。PGA在該網(wǎng)絡(luò)模型中占據(jù)著較高的權(quán)重，為場(chǎng)地等效剪切波速變化的最主要控制參數(shù)。該網(wǎng)絡(luò)模型捕捉到場(chǎng)地等效剪切波速比隨PGA的增大有下降的趨勢(shì)，而較為松軟的D、E類場(chǎng)地受PGA影響更大，下降幅度更大。

關(guān)鍵詞：神經(jīng)網(wǎng)絡(luò)；等效剪切波速；場(chǎng)地非線性；參數(shù)預(yù)測(cè)；地表峰值加速度

中圖分類號(hào)：P315.7 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-0666（2024）02-0280-10

doi：10.20015/j.cnki.ISSN1000-0666.2024.0016

0 引言

場(chǎng)地非線性研究是地震工程領(lǐng)域的重要課題之一，它對(duì)近斷層地震動(dòng)模擬、地震動(dòng)預(yù)測(cè)方程的建立、震害評(píng)估等都具有重要的理論意義和工程應(yīng)用價(jià)值。早期的研究表明，場(chǎng)地非線性反應(yīng)閾值介于100～200 cm/s²（Beresnev，Wen，1996），近期的研究則表明場(chǎng)地非線性反應(yīng)閾值低于該值，如Rubinstein（2011）通過(guò)美國(guó)加州地震動(dòng)記錄研究發(fā)現(xiàn)地震動(dòng)強(qiáng)度達(dá)到35 cm/s²時(shí)就能觀察到場(chǎng)地非線性反應(yīng)；一些學(xué)者基于日本KiK-Net臺(tái)網(wǎng)的地震動(dòng)觀測(cè)數(shù)據(jù)證實(shí)了非線性場(chǎng)地反應(yīng)普遍存在于中等強(qiáng)度的地震動(dòng)中，如Ghofrani等（2013）估計(jì)的49個(gè)臺(tái)陣的場(chǎng)地非線性反應(yīng)閾值介于40～150 cm/s²；Wu等（2010）估計(jì)的6個(gè)KiK-Net臺(tái)陣的場(chǎng)地非線性反應(yīng)閾值介于20～80 cm/s²；Regnier等（2013）估計(jì)的36個(gè)KiK-Net臺(tái)陣的場(chǎng)地非線性反應(yīng)閾值介于20～500 cm/s²。

近年來(lái)，有大量的研究表明神經(jīng)網(wǎng)絡(luò)在處理復(fù)雜回歸預(yù)測(cè)問(wèn)題時(shí)有較好表現(xiàn)，且在處理地震動(dòng)預(yù)測(cè)問(wèn)題時(shí)具有優(yōu)勢(shì)（Derras et al，2014；萬(wàn)永革，李鴻吉，1995），基于數(shù)據(jù)驅(qū)動(dòng)的人工神經(jīng)網(wǎng)絡(luò)可以在不構(gòu)建具體函數(shù)關(guān)系式的前提下直接構(gòu)建端到端的預(yù)測(cè)模型，如Kerb和Ting（2005）基于多層前饋神經(jīng)網(wǎng)絡(luò)，利用30個(gè)臺(tái)站的歷史地震記錄，估算了高速鐵路系統(tǒng)車站的峰值加速度；Alavi和Gandomi（2011）將人工神經(jīng)網(wǎng)絡(luò)和模擬退火算法相結(jié)合，使用震級(jí)、震中距、場(chǎng)地等效剪切波速和斷層機(jī)制來(lái)預(yù)測(cè)地面峰值加速度、峰值速度和峰值位移；Dhanya和Raghukanth（2017）基于人工神經(jīng)網(wǎng)絡(luò)構(gòu)建了一個(gè)適用于淺源地震的地震動(dòng)參數(shù)預(yù)測(cè)模型，利用震級(jí)、震中距、場(chǎng)地剪切波速V_S30和震源機(jī)制這4個(gè)參數(shù)預(yù)測(cè)地面峰值速度、加速度和5％的阻尼譜加速度。

場(chǎng)地非線性所導(dǎo)致的場(chǎng)地等效剪切波速的下降是場(chǎng)地分析中的重要問(wèn)題，目前的研究中多使用傳統(tǒng)回歸的方法來(lái)回歸單個(gè)臺(tái)站的變化曲線，該方法僅能反映特定臺(tái)站的情況且無(wú)法應(yīng)用于部分?jǐn)?shù)據(jù)量較少的臺(tái)站，存在一定的局限性。本文建立了一個(gè)BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型，用來(lái)預(yù)測(cè)在地震動(dòng)作用下，場(chǎng)地的剪切波速下降情況并判斷產(chǎn)生非線性變化的大致區(qū)間。

1 數(shù)據(jù)庫(kù)創(chuàng)建

1.1 數(shù)據(jù)選取

日本KiK-net臺(tái)網(wǎng)具有較為完整的場(chǎng)地信息和大量的井上井下地震動(dòng)信息，站井的深度分布介于100～2 500 m。本文搜集整理了KiK-net臺(tái)網(wǎng)觀測(cè)到的地表與井下地震動(dòng)記錄。考慮到大震發(fā)生后場(chǎng)地的剪切波速會(huì)在短期內(nèi)快速恢復(fù)到90%以上并呈現(xiàn)出長(zhǎng)期緩慢恢復(fù)的性質(zhì)（苗雨等，2018），因此對(duì)PGAgt;500 cm/s²地震后100 d內(nèi)的數(shù)據(jù)進(jìn)行刪減，將各地震動(dòng)事件東西向及南北向作為兩次獨(dú)立的觀測(cè)結(jié)果（董凱月，2020），最終篩選出近20年來(lái)具有完整地質(zhì)資料的407個(gè)臺(tái)站中的30 952條地震動(dòng)記錄。對(duì)數(shù)據(jù)進(jìn)行基線校正和1～13 Hz濾波處理后，數(shù)據(jù)庫(kù)中PGAlt;20 cm/s²的記錄共17 584條，PGAgt;20 cm/s²的記錄部分組成如圖1所示。20 cm/s²lt;PGAlt;60 cm/s²的記錄共7 910條，PGAgt;60 cm/s²的記錄共有5 458條。本文所選地震動(dòng)記錄的震級(jí)與震中距散點(diǎn)分布如圖2所示。根據(jù)美國(guó)NEHRP規(guī)范（表1），統(tǒng)計(jì)出每類場(chǎng)地地震動(dòng)記錄所占比例，如圖3所示，圖中無(wú)A類場(chǎng)地，B類場(chǎng)地所占比例為7%，E類場(chǎng)地所占比例最低，為2.1%，本文選取的臺(tái)站以C類和D類為主。

1.2 等效剪切波速估計(jì)

本文通過(guò)井上井下解卷積求取場(chǎng)地等效剪切波速（張亞沛，2019；劉啟方，陳長(zhǎng)龍，2022）。設(shè)地表和井下加速度記錄的傅立葉譜分別為U₁（ω）和U₂（ω），應(yīng)用解卷積法首先計(jì)算系統(tǒng)的傳遞函數(shù)G（ω）：

式中：*表示復(fù)共軛，分子為地表與井下加速度在頻域內(nèi)的互相關(guān)函數(shù)，分母為井下加速度的功率譜。為了避免分母中譜值很小引起的不穩(wěn)定性，一般采用U′₂（ω）對(duì)U₂（ω）進(jìn)行校正：

對(duì)G（ω）進(jìn)行傅立葉反變換到時(shí)域g（t），即為井下脈沖輸入時(shí)地表的響應(yīng)。通過(guò)脈沖響應(yīng)（脈沖的峰值到時(shí)）獲得井下臺(tái)站到地表地震波的傳播時(shí)間，用該站井深度除以傳播時(shí)間，計(jì)算該場(chǎng)地土體等效剪切波速V_S。

場(chǎng)地等效剪切波速比V_Sr計(jì)算公式為：

式中：V_S0為該場(chǎng)地的原始等效剪切波速值，取該場(chǎng)地PGAlt;20 cm/s²的弱震動(dòng)作用下場(chǎng)地等效剪切波速的平均值（苗雨等，2018，2021）。

2 BP神經(jīng)網(wǎng)絡(luò)模型搭建

2.1 BP神經(jīng)網(wǎng)絡(luò)模型

BP神經(jīng)網(wǎng)絡(luò)是一種按照誤差逆向傳播算法訓(xùn)練的多層前饋神經(jīng)網(wǎng)絡(luò)。該網(wǎng)絡(luò)模型具有函數(shù)逼近的屬性，可以在不考慮輸入與輸出間具體函數(shù)關(guān)系的前提下，通過(guò)機(jī)器學(xué)習(xí)的方法，構(gòu)建出一種輸入與輸出的映射，擬合任意復(fù)雜的連續(xù)函數(shù)。BP神經(jīng)網(wǎng)絡(luò)包括一個(gè)輸入層、多個(gè)隱藏層和一個(gè)輸出層，每一層由不同數(shù)量的人工神經(jīng)元組成，網(wǎng)絡(luò)結(jié)構(gòu)如圖4所示。

隱藏層中第i個(gè)神經(jīng)元的輸出可以寫為：

y_ij=φ（∑w_ij+b_j）""" （4）

式中：w_ij表示神經(jīng)網(wǎng)絡(luò)中第i層與第i-1層神經(jīng)元之間的權(quán)重；b_j表示神經(jīng)網(wǎng)絡(luò)中第i層第j個(gè)神經(jīng)元線性層的偏置；φ表示激活函數(shù)。激活函數(shù)的類型有很多種，本文輸出層采用線性激活函數(shù)linear（式5），隱藏層采用tanh激活函數(shù)（式6）：

linear（x）=x """"（5）

2.2 網(wǎng)絡(luò)模型的開(kāi)發(fā)及評(píng)價(jià)指標(biāo)

目前，研究人員多以PGA來(lái)表征場(chǎng)地的非線性閾值（Chin，Aki，1991；Wu et al，2010；Wang et al，2019），由于該值僅能表征某次事件的峰值強(qiáng)度，因此本文還使用了Arias烈度（Ia）來(lái)表征某次事件整體能量大小（Benito，Herraiz，1997；Chousianitis et al，2014，2016），以及場(chǎng)地剪切波速V_S30來(lái)表征場(chǎng)地情況。經(jīng)過(guò)多次嘗試，建立了一個(gè)包含3個(gè)輸入?yún)?shù)（PGA、I_a、V_S30）、5個(gè)隱藏神經(jīng)元、1個(gè)輸出參數(shù)（V_Sr）的BP神經(jīng)網(wǎng)絡(luò)模型，其預(yù)測(cè)過(guò)程可表示為：

V_Sr=f ［PGA，ln（I_a），V_S30］""" （7）

式中：I_a為Arias烈度，由于其值分布范圍過(guò)廣，因此取其自然對(duì)數(shù)，以減小其分布范圍，從而加快神經(jīng)網(wǎng)絡(luò)梯度的更新，以有效提高神經(jīng)網(wǎng)絡(luò)的收斂速度，加速模型的訓(xùn)練。

本文采用Adam優(yōu)化算法進(jìn)行神經(jīng)網(wǎng)絡(luò)模型訓(xùn)練，它可以根據(jù)學(xué)習(xí)環(huán)境自適應(yīng)調(diào)整學(xué)習(xí)率，使得神經(jīng)網(wǎng)絡(luò)得以更快收斂，同時(shí)對(duì)初始的超參數(shù)設(shè)置不敏感，方便了神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)（Kingma，Ba，2014）。對(duì)模型預(yù)測(cè)效果的評(píng)價(jià)，本文選用MSE作為模型的損失函數(shù)（Dhanya et al，2017），使用MSE，其計(jì)算公式分別為及皮爾遜相關(guān)系數(shù)作為模型的評(píng)價(jià)指標(biāo)：

式中：x_i為目標(biāo)數(shù)據(jù)或記錄數(shù)據(jù)；y_i為預(yù)測(cè)數(shù)據(jù)；m為該批數(shù)據(jù)點(diǎn)總數(shù)；x、 y為該批次目標(biāo)數(shù)據(jù)和預(yù)測(cè)數(shù)據(jù)的均值。

3 網(wǎng)絡(luò)模型性能分析

3.1 模型性能

本文篩選出地震動(dòng)記錄較多、PGA分布較廣的6個(gè)包含不同場(chǎng)地類型的臺(tái)站作為測(cè)試臺(tái)站進(jìn)行圖像分析及模型評(píng)價(jià)，其場(chǎng)地情況見(jiàn)表2，測(cè)試臺(tái)站不參與訓(xùn)練。

本文使用場(chǎng)地非線性研究中常用的傳統(tǒng)雙曲線回歸模型作為模型1，表示為：

式中：x為PGA；y為場(chǎng)地等效剪切波速比；k為回歸參數(shù)。使用模型1，對(duì)6個(gè)臺(tái)站的V_Sr關(guān)于PGA的變化做了回歸分析，擬合參數(shù)見(jiàn)表3。通過(guò)訓(xùn)練所得BP神經(jīng)網(wǎng)絡(luò)模型對(duì)6個(gè)臺(tái)站的V_Sr進(jìn)行預(yù)測(cè)，模型1與神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果如圖5所示。由于V_Sr為場(chǎng)地等效剪切波速與弱震動(dòng)作用下場(chǎng)地等效剪切波速均值的比值，因此比弱震動(dòng)作用下場(chǎng)地等效剪切波速均值高的計(jì)算值會(huì)出現(xiàn)V_Sr大于1.0的結(jié)果。

由圖5可見(jiàn)，用BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)所得曲線與模型1擬合曲線整體上較為相似，兩者均呈現(xiàn)出V_Sr隨PGA的增大而逐漸減小的特性。在測(cè)試臺(tái)站的各場(chǎng)地類型中，當(dāng)PGA=50 cm/s²時(shí)，用BP神經(jīng)網(wǎng)絡(luò)模型得到的V_Sr預(yù)測(cè)值均未小于0.99；當(dāng)PGA=100 cm/s²時(shí)，各場(chǎng)地V_Sr下降幅度約為2%～3%，場(chǎng)地V_Sr開(kāi)始呈現(xiàn)非線性變化，多數(shù)場(chǎng)地的非線性閾值約為50～100 cm/s²；當(dāng)PGAgt;100 cm/s²后，場(chǎng)地V_Sr的非線性變化較為顯著。在測(cè)試臺(tái)站中，對(duì)于較硬的B類場(chǎng)地，其PGA達(dá)到175 cm/s²左右時(shí)，場(chǎng)地V_Sr下降約5%；對(duì)于較松軟的D、E類場(chǎng)地，PGA達(dá)到140 cm/s²左右時(shí)，V_Sr即下降約5%。圖中E類場(chǎng)地的BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)曲線波動(dòng)較大，推測(cè)是由于E類場(chǎng)地訓(xùn)練數(shù)據(jù)相對(duì)較少，訓(xùn)練不夠充分所導(dǎo)致。訓(xùn)練得到的BP神經(jīng)網(wǎng)絡(luò)模型各層參數(shù)見(jiàn)表4、5。

訓(xùn)練完成后，訓(xùn)練集的MSE為6×10^-4，測(cè)試臺(tái)站數(shù)據(jù)集的MSE為6.7×10^-4。訓(xùn)練集與測(cè)試臺(tái)站數(shù)據(jù)集的均方誤差相似，模型整體表現(xiàn)良好，未出現(xiàn)過(guò)擬合現(xiàn)象。

使用BP神經(jīng)網(wǎng)絡(luò)模型及模型1計(jì)算得到6個(gè)臺(tái)站的Person相關(guān)系數(shù)，結(jié)果見(jiàn)表6。從表中可見(jiàn)，BP神經(jīng)網(wǎng)絡(luò)模型對(duì)場(chǎng)地等效剪切波速比的預(yù)測(cè)相較于模型1，其相關(guān)系數(shù)相對(duì)較高，其預(yù)測(cè)表現(xiàn)均不弱于模型1，整體上有著更好的預(yù)測(cè)性能。

3.2 殘差分析

通常使用殘差分析作為主要手段分析神經(jīng)網(wǎng)絡(luò)的性能。殘差分析可以反映預(yù)測(cè)的輸出與輸入變量的整體性關(guān)系，可分析該網(wǎng)絡(luò)在整個(gè)數(shù)據(jù)區(qū)間上整體預(yù)測(cè)性能的好壞（Dhanya，Raghukanth，2020）。本文將殘差定義為預(yù)測(cè)值與目標(biāo)值之差，如下所示：

R_i=y_pred-y_target""" （11）

圖6展示了BP神經(jīng)網(wǎng)絡(luò)模型3個(gè)輸入變量PGA、ln（I_a）和V_S30對(duì)應(yīng)于輸出變量V_Sr的模型殘差，圖中橙色誤差棒為各區(qū)間上的所有殘差作均值及標(biāo)準(zhǔn)差處理。表7給出了3個(gè)輸入變量在各區(qū)間的殘差分布統(tǒng)計(jì)，由表可見(jiàn)，在各個(gè)區(qū)間上，該網(wǎng)絡(luò)模型均值殘差均在0附近，整體呈現(xiàn)出無(wú)偏特性，對(duì)于PGA及I_a較高的強(qiáng)震數(shù)據(jù)，可能是由于數(shù)據(jù)量較少且數(shù)據(jù)的離散度較大，因此其殘差有著相對(duì)更大的標(biāo)準(zhǔn)差。V_S30為120～1 270 m/s時(shí)，均有較低的均值殘差及標(biāo)準(zhǔn)差，表明其在多數(shù)軟硬場(chǎng)地中均有較好的預(yù)測(cè)性能。

3.3 敏感性分析

對(duì)于具有明確函數(shù)形式的回歸分析，每個(gè)變量的相對(duì)重要性可以從優(yōu)化后得到的系數(shù)中確定。但是在BP神經(jīng)網(wǎng)絡(luò)中，由于網(wǎng)絡(luò)本身是高度復(fù)雜的互聯(lián)節(jié)點(diǎn)，所以想要得到其相關(guān)性較為困難。許多研究者嘗試了各種技術(shù)來(lái)獲得輸入節(jié)點(diǎn)對(duì)輸出的相對(duì)影響（Gevrey et al，2003），其中最簡(jiǎn)單直接的方法之一是Garson（1991）提出的基于隱藏神經(jīng)元的連接權(quán)值與輸入節(jié)點(diǎn)相關(guān)權(quán)值的劃分，根據(jù)該方法，特定輸出變量j的每個(gè)輸入變量i的相關(guān)重要性RI為：

式中：W_i，m為輸入層-隱藏層各神經(jīng)元所對(duì)應(yīng)的權(quán)值；W_j，m為隱藏層-輸出層各神經(jīng)元權(quán)值；n和p分別為輸入?yún)?shù)個(gè)數(shù)及隱藏層神經(jīng)元數(shù)量。

通過(guò)表4、5中的權(quán)值，將式（12）應(yīng)用于BP神經(jīng)網(wǎng)絡(luò)模型上，所得到的輸入?yún)?shù)相對(duì)輸出參數(shù)V_Sr的相關(guān)權(quán)重如圖7所示。由圖7可見(jiàn)，PGA在該模型中所占權(quán)重較大，為該模型的主要控制參數(shù)。

控制部分輸入變量不變，繪制輸出變量隨某一輸入變量的變化關(guān)系曲線，是神經(jīng)網(wǎng)絡(luò)模型驗(yàn)證的常用手段。由于PGA為該網(wǎng)絡(luò)模型的主要控制參數(shù)，而Arias烈度參數(shù)的權(quán)重較低，僅占23.3%，因此本文選取4種不同的V_S30以涵蓋4類場(chǎng)地（圖8），基于模型的預(yù)測(cè)結(jié)果探討PGA和V_Sr變化的關(guān)系。

圖8為PGA關(guān)于V_Sr變化關(guān)系。由圖8可見(jiàn)，V_Sr隨PGA的增大有下降的趨勢(shì)，PGA越大，V_Sr下降越快。同時(shí)，V_S30越小的松軟場(chǎng)地受PGA影響越大，下降幅度也越大。在以往的研究中普遍認(rèn)為在堅(jiān)硬的場(chǎng)地中僅會(huì)產(chǎn)生很小的非線性變化，但在本文中B類場(chǎng)地在PGA較大時(shí)仍有一定的等效剪切波速比下降，推測(cè)可能是因?yàn)閿?shù)據(jù)庫(kù)中B類場(chǎng)地?cái)?shù)據(jù)所占比例較少，且V_S30與場(chǎng)地整體土層情況仍有一定差異，導(dǎo)致該網(wǎng)絡(luò)模型預(yù)測(cè)結(jié)果在V_S30較高的場(chǎng)地中仍有一定程度的等效剪切波速比下降。

4 結(jié)論

本文使用日本KiK-net臺(tái)網(wǎng)的407個(gè)臺(tái)站記錄到的30 952條記錄，提出了一種基于BP神經(jīng)網(wǎng)絡(luò)的場(chǎng)地等效剪切波速比變化預(yù)測(cè)模型。模型采用均方誤差函數(shù)以及Adam優(yōu)化算法，包含3個(gè)輸入?yún)?shù)、5個(gè)隱藏神經(jīng)元和1個(gè)輸出參數(shù)。輸入?yún)?shù)為加速度峰PGA、Arias烈度I_a及場(chǎng)地剪切波速度V_S30，輸出為場(chǎng)地等效剪切波速比V_Sr，得到以下結(jié)論：

（1）神經(jīng)網(wǎng)絡(luò)模型可以在不構(gòu)建復(fù)雜函數(shù)回歸關(guān)系的情況下，預(yù)測(cè)場(chǎng)地的剪切波速變化情況。本文選取的包含不同場(chǎng)地類型的6個(gè)臺(tái)站，通過(guò)使用BP神經(jīng)網(wǎng)絡(luò)對(duì)6個(gè)臺(tái)站的場(chǎng)地剪切波速比的預(yù)測(cè)曲線可知，BP網(wǎng)絡(luò)模型預(yù)測(cè)曲線相對(duì)光滑，未出現(xiàn)過(guò)擬合現(xiàn)象。相較于使用最小二乘法回歸所得的剪切波速比函數(shù)曲線，BP神經(jīng)網(wǎng)絡(luò)模型的相關(guān)系數(shù)較高，該神經(jīng)網(wǎng)絡(luò)模型殘差對(duì)于各輸入變量整體呈現(xiàn)出無(wú)偏特性，在大部分軟硬場(chǎng)地中均有較好的預(yù)測(cè)性能。

（2）由BP神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)曲線可知，PGA為50～100 cm/s²時(shí)，多數(shù)場(chǎng)地開(kāi)始出現(xiàn)非線性變化，剪切波速開(kāi)始發(fā)生變化；PGAgt;100 cm/s²時(shí)，場(chǎng)地發(fā)生顯著的非線性變化，剪切波速開(kāi)始發(fā)生較快變化。PGA達(dá)到175 cm/s²左右時(shí)，B類場(chǎng)地剪切波速下降5%，PGA達(dá)到約140 cm/s²時(shí)，D、E類場(chǎng)地剪切波速下降5%。

（3）模型敏感性分析表明，PGA在該神經(jīng)網(wǎng)絡(luò)模型中占據(jù)著較高的權(quán)重，為場(chǎng)地等效剪切波速比變化的最主要控制參數(shù)。該神經(jīng)網(wǎng)絡(luò)模型捕捉到場(chǎng)地等效剪切波速比隨PGA的增大有下降的趨勢(shì)，而較為松軟的D、E類場(chǎng)地受PGA影響更大，下降的幅度也更大。

（4）在數(shù)據(jù)較少的B、E類場(chǎng)地中，基于數(shù)據(jù)驅(qū)動(dòng)的神經(jīng)網(wǎng)絡(luò)模型存在一定的訓(xùn)練風(fēng)險(xiǎn)，是后續(xù)研究中需要討論的一個(gè)關(guān)鍵問(wèn)題。

本文使用了日本KiK-net臺(tái)網(wǎng)提供的地震動(dòng)數(shù)據(jù)，在此表示感謝。

參考文獻(xiàn)：

董凱月.2020.基于深度神經(jīng)網(wǎng)絡(luò)的地下地震動(dòng)參數(shù)預(yù)測(cè)研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué).

Dong K Y.2020.Study on the prediction of underground motion parameters based on deep neural network［D］.Harbin：Harbin Institute of Technology.（in Chinese）

劉啟方，陳長(zhǎng)龍.2022.基于移動(dòng)窗解卷積法的2011年日本東北M_W9.0大地震場(chǎng)地非線性時(shí)變識(shí)別［J］.地震學(xué)報(bào)，44（1）：96-110.

Liu Q F，Chen C L.2022.Identification of the temporal changes of site nonlinearity during 2011 M_W9.0 Tohoku earthquake by moving time window deconvolution method［J］.Acta Seismologica Sinica，44（1）：96-110.（in Chinese）

苗雨，施洋，王蘇陽(yáng)，等.2018.基于豎向臺(tái)陣地震記錄的非線性場(chǎng)地反應(yīng)研究［J］.自然災(zāi)害學(xué)報(bào)，27（6）：51-58.

Miao Y，Shi Y，Wang S Y，et al.2018.Assessing nonlinear soil behavior using vertical array data：A case at TCGH16 station from KiK-Net in Japan［J］.Journal of Natural Disasters，27（6）：51-58.（in Chinese）

苗雨，施洋，張昊，等.2021.基于地震觀測(cè)記錄的場(chǎng)地豎向非線性性質(zhì)［J］.北京工業(yè)大學(xué)學(xué)報(bào)，47（7）：669-679.

Miao Y，Shi Y，Zhang H，et al.2021.Investigating vertical nonlinearity of site using KiK-net seismic data［J］.Journal of Beijing University of Technology，47（7）：669-679.（in Chinese）

萬(wàn)永革，李鴻吉.1995.人工神經(jīng)網(wǎng)絡(luò)在地球物理中的應(yīng)用綜述［J］.國(guó)際地震動(dòng)態(tài)，（1）：9-14.

Wan Y G，Li H J.1995.Overview of the application of artificial neural networks in Geophysics［J］.Progress in Earthquake Sciences，（1）：9-14.（in Chinese）

張亞沛.2019.基于KIK-NET數(shù)據(jù)庫(kù)的土層特性識(shí)別及地下地震動(dòng)幅值特性研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué).

Zhang Y P.2019.Research on identification of soil characteristics and amplitude characteristics of underground motion based on KiK-Net database［D］.Harbin：Harbin Institute of Technology.（in Chinese）

Alavi A H，Gandomi A H.2011.Prediction of principal ground-motion parameters using a hybrid method coupling artificial neural networks and simulated annealing［J］.Computers amp; Structures，89（23-24）：2176-2194.

Benito B，Herraiz M.1997.An approach to the measurement of the potential structural damage of earthquake ground motions［J］.Earthquake Engineering and Structural Dynamics，26：79-92.

Beresnev I A，Wen K L.1996.Nonlinear soil response—A reality？［J］.Bulletin of the Seismological Society of America，86（6）：1964-1978.

Chin B H，Aki K.1991.Simultaneous study of the source path and site effects on strong ground motion during the 1989 Loma Prieta earthquake：a preliminary result on pervasive nonlinear site effects［J］.Bulletin of the Seismological Society of America，81（5）：1859-1884.

Chousianitis K，Del G V，Kalogeras I，et al.2014.Predictive model of Arias intensity and Newmark displacement for regional scale evaluation of earthquake-induced landslide hazard in Greece［J］.Soil Dynamics and Earthquake Engineering，65（4）：11-29.

Chousianitis K，Del G V，Sabatakakis N，et al.2016.Assessment of earthquake-induced landslide hazard in Greece：From Arias intensity to spatial distribution of slope resistance demand［J］.Bulletin of the Seismological Society of America，106（1）：174-188.

Derras B，Bard P Y，Cotton F.2014.Towards fully data driven ground-motion prediction models for Europe［J］.Bulletin Earthquake Engineering，12（1）：495-516.

Dhanya J，Raghukanth S T G.2017.Ground motion prediction model using artificial neural network［J］.Pure and Applied Geophysics，175（3）：1035-1064.

Dhanya J，Raghukanth S T G.2020.Neural network-based hybrid ground motion prediction equations for Western Himalayas and North-Eastern India［J］.Acta Geophysica，172（3）：762-781.

Garson G D.1991.Interpreting neural-network connection weights［J］.AI Expert，6：47-51.

Gevrey M，Dimopoulos I，Lek S.2003.Review and comparison of methods to study the contribution of variables in artifificial neural network models［J］.Ecological Modelling，160（3）：249-264.

Ghofrani H，Atkinson G M，Goda K.2013.Implications of the 2011 M9.0 Tohoku Japan earthquake for the treatment of site effects in large earthquakes［J］.Bulletin of Earthquake Engineering，11（1）：171-23.

Kerb T，Ting S B.2005.Neural network estimation of ground peak acceleration at stations along Taiwan high-speed rail system［J］.Engineering Applications of Artificial Intelligence，18（7）：857-866.

Kingma D P，Ba J.2014.Adam：A method for stochastic optimization［C］.ICLA，arXiv preprint arXiv：1412.6980.

Regnier J，Cadet H，Bonilla L F，et al.2013.Assessing nonlinear behavior of soils in seismic site response：Statistical analysis on Kik-net strong-motion data［J］.Bulletin of the Seismological Society of America，103（3）：1750-1770.

Rubinstein J L.2011.Nonlinear site response in medium magnitude earthquakes near Parkfield California［J］.Bulletin of the Seismological Society of America，101（1）：275-286.

Wang H Y，Jiang W P，Wang S Y，et al.2019.In situ assessment of soil dynamic parameters for characterizing nonlinear seismic site response using KiK-net vertical array data［J］.Bulletin of Earthquake Engineering，17（5）：2331-2360.

Wu C，Peng Z，Benzion Y.2010.Refined thresholds for non-linear ground motion and temporalchanges of siteresponse associated with medium-size earthquakes［J］.Geophysical JournalInternational，182（3）：1567-1576.

Research on Prediction of Site Equivalent Shear Wave Velocity

Change Based on BP Neural Network

SU Wenhao，LIU Qifang

（Suzhou University of Science and Technology，Key Laboratory of Structure Engineering

of Jiangsu Province，Suzhou 215009，Jiangsu，China）

Abstract

In this paper，a total of 30952 records from 407 stations of the Japanese KiK-net is used to propose a prediction model for the change of equivalent shear wave velocity ratio based on BP neural network.The model adopts the mean square error function and the Adam optimization algorithm，consists of three inputs，five hidden neurons and one output.The input parameters are Peak Ground Acceleration（PGA），Arias intensity（I_a）and site V_S30.The output parameter is site equivalent shear wave velocity ratio（V_Sr）.The research shows that the residual error of the network model is unbiased for each input variable，and has good prediction performance in many kinds of sites.Compared with the function curve of the traditional least-square method，the neural network model has a relatively better performance.In the prediction curve of the network model，the shear wave velocity of the site of Class B decreases by 5% when the PGA reaches about 175 cm/s²，and the shear wave velocity of the sites of Class D and E decreases by 5% when the PGA reaches about 140 cm/s².The nonlinear threshold of most sites is between 50～100 cm/s².PGA occupies a high weight in the network model and is the main controlling parameter of the site equivalent shear wave velocity change.The network model captures that the equivalent shear wave velocity ratio of the site has a downward trend with the increase of PGA.At the same time，it shows that the Class D and E sites are greatly affected by PGA，and the declineing range is larger.

Keywords：neural network；equivalent shear wave velocity；site nonlinearity；parameter prediction；peak ground acceleration