周田宇，梁富源，馬完君，謝志南

（中國地震局工程力學(xué)研究所地震工程與工程振動重點(diǎn)實(shí)驗(yàn)室；地震災(zāi)害防治應(yīng)急管理部重點(diǎn)實(shí)驗(yàn)室，哈爾濱 150080）

引言

局部起伏地形對地震動，尤其是高頻地震動存在較為顯著的放大或抑制作用，簡稱“地形效應(yīng)”。若不考慮地形場地淺表土層的影響，一般凸起地形頂部對高頻地震動存在放大作用，凹陷地形底部對高頻地震動存在抑制作用［1-2］。其中：凸起地形對地震動放大進(jìn)而導(dǎo)致震害加劇的現(xiàn)象在多次強(qiáng)震震害調(diào)查中均有體現(xiàn)，如1971年San Fernando地震［3］；1975年海城地震［1］；1989年Loma Prieta地震［4］；1994年Northridge地震［5］；2008年汶川地震［6］；2010年Haiti地震［7］；2014年蘆山地震［8］。我國《建筑與市政工程抗震通用規(guī)范》（GB 55002-2021）亦明確指出：考慮局部突出地形對地震動參數(shù)（設(shè)防烈度）的放大作用，這對山區(qū)建筑的抗震計(jì)算十分必要。

地形效應(yīng)研究大致分兩類：形成機(jī)理研究和定量表征模型構(gòu)建。研究方法包括：參考臺站法，解析方法，數(shù)值方法和振動臺試驗(yàn)。形成機(jī)理研究結(jié)果表明：地形效應(yīng)影響因素復(fù)雜多樣，主要包括：地形的幾何形狀（山體、邊坡和河谷）［2］；地形場地入射地震動的類型（P波、SV波、SH波和Rayleigh波）、角度（垂直入射和斜入射）和頻譜分布（高頻成分占優(yōu)和低頻成分占優(yōu)）［9-11］；震源參數(shù)（震源與場點(diǎn)之間的方位、震源深度、震源機(jī)制）［12］；場地淺表土層波速與阻尼結(jié)構(gòu)等［13］。然而，量化地形效應(yīng)影響因素較為困難。原因有二：（1）不同地震中上述因素對地震動的放大和抑制作用不一致。（2）影響因素自身難以量化，比如場地淺表土層與地形的耦合作用。已有定量表征模型大多僅考慮了地形幾何形狀參數(shù)的影響。例如：利用NGA-West2數(shù)據(jù)庫，RAI等［14］構(gòu)建了適用于CY14地震動衰減關(guān)系（亦稱為地震動預(yù)測方程）的地形效應(yīng)修正模型。模型參數(shù)為場點(diǎn)高程與場點(diǎn)所在1 500 m范圍圓形區(qū)域平均高程的差值。模型給出了對0.2 s以上周期反應(yīng)譜的修正系數(shù)：就0.25 s至1.5 s周期反應(yīng)譜而言，高程差（場點(diǎn)高程減去場點(diǎn)及鄰域的平均高程）大于20 m場點(diǎn)的修正系數(shù)為正，最大修正系數(shù)為0.12，對應(yīng)于場點(diǎn)地形對反應(yīng)譜的放大倍數(shù)為1.13；就0.2 s至10 s周期反應(yīng)譜而言，高程差小于-20 m場點(diǎn)的修正系數(shù)為負(fù)，最小修正系數(shù)為-0.29，對應(yīng)于場點(diǎn)地形對反應(yīng)譜的放大倍數(shù)為0.75。但這一修正模型是否可以獨(dú)立應(yīng)用于表征地形效應(yīng)有待進(jìn)一步研究。利用二維數(shù)值模擬數(shù)據(jù)庫，郝明輝等［15］以坡角為參數(shù)，針對梯形狀凸起地形，建立了地形對PGA、0.2 s周期反應(yīng)譜和1.0 s周期反應(yīng)譜的放大倍數(shù)模型。放大倍數(shù)隨坡角的增加而增大。坡角大于60°情形，凸起地形對PGA的平均放大倍數(shù)為1.547。利用三維數(shù)值模擬數(shù)據(jù)庫，MAUFROY等［16］針對一般地形以頻率標(biāo)度曲率為參數(shù)，建立了地形對高頻段地震動傅里葉幅值譜的放大倍數(shù)模型，模型所給放大倍數(shù)與直接基于數(shù)值模擬記錄計(jì)算所得放大倍數(shù)誤差在±30%以內(nèi)；后續(xù)，WANG等［13］發(fā)展了考慮場地淺表土層厚度和阻尼比參數(shù)影響的地形效應(yīng)表征模型。這一模型是否獨(dú)立于場地效應(yīng)模型（淺表土層對地震動的放大或抑制，如VS30模型）有待深入分析且這兩類參數(shù)就一般山區(qū)地形場地而言難以獲取。

當(dāng)前，工程常用的寬頻帶地震動模擬與預(yù)測方法未考慮局域地形對高頻地震動的影響。如基于地震動參數(shù)預(yù)測方程的寬頻帶地震動合成方法，結(jié)合低頻地震動波動模擬和高頻地震動隨機(jī)模擬的寬頻帶地震動模擬與預(yù)測混合方法以及單獨(dú)基于隨機(jī)模擬的寬頻帶地震動模擬與預(yù)測方法?；谏鲜龇椒M和預(yù)測地震動開展高頻地震動起主導(dǎo)作用的地震災(zāi)害和地震次生災(zāi)害評估可能導(dǎo)致對震害的錯(cuò)誤估計(jì)。為此，本文利用相對獨(dú)立于場地效應(yīng)的Maufroy地形效應(yīng)表征模型，以中國地震科學(xué)實(shí)驗(yàn)場建設(shè)區(qū)域-川滇區(qū)域?yàn)槔?，?gòu)建區(qū)域高頻地震動地形效應(yīng)的計(jì)算方案，并在此基礎(chǔ)上探討計(jì)算結(jié)果在后續(xù)寬頻地震動模擬與預(yù)測方法中的應(yīng)用，為地震科學(xué)試驗(yàn)場“從地震破裂過程到工程結(jié)構(gòu)響應(yīng)”全鏈條地震動模擬與預(yù)測平臺的構(gòu)建提供一定參考。

1 Maufroy地形效應(yīng)量化表征模型

本文選用Maufory地形效應(yīng)量化表征模型開展高頻地震動地形效應(yīng)分析［16］，原因在于該模型構(gòu)建過程：（1）將地形效應(yīng)本質(zhì)歸結(jié)為復(fù)雜三維地震波傳播效應(yīng)。（2）將場點(diǎn)地形效應(yīng)與場地效應(yīng)分離。（3）考慮了地形效應(yīng)分析結(jié)果的不確定性。另外，該模型有望后續(xù)拓展至三分量高頻地震動地形效應(yīng)分析。該模型的詳細(xì)構(gòu)建過程見參考文獻(xiàn)［16］。Maufroy地形效應(yīng)量化表征模型為統(tǒng)計(jì)得到的地形對不同頻率傅里葉幅值譜的放大倍數(shù)在不同分位數(shù)處的取值為：

式中：AF16th(f)、AF50th(f)和AF84th(f)分別為16%、50%和84%分位數(shù)取值；λs為頻率f處對應(yīng)的地震波波長，λs=VSf，VS為彈性半空間剪切波波速；CS(LS)為頻率標(biāo)度曲率，其中特征長度為：式中：h為曲率計(jì)算地形面剖分尺寸；n為計(jì)算頻率標(biāo)度曲率所采用的平滑矩陣的維度。

式中：C為地形面上任意點(diǎn)的近似曲率，以(xi,yi)點(diǎn)為例：

式中：E(xi,yi)為(xi,yi)地形場點(diǎn)的高程。AF16th(f)、AF50th(f)和AF84th(f)計(jì)算流程可歸納如下：（1）確定場點(diǎn)VS。（2）確定地形面剖分尺寸h，一般直接取為數(shù)字地形的最高分辨率。（3）利用數(shù)字地形計(jì)算曲率C，并利用式（5）分別計(jì)算n=3,5,7,…情形下的LS處對應(yīng)的頻率標(biāo)度曲率CS(LS)。（4）利用f=VSλs確定計(jì)算頻率，λs由式（4）得到。（5）將λs和CS(LS)分別代入式（1）、式（2）和式（3）得到AF16th(f)、AF50th(f)和AF84th(f)。

2 區(qū)域地形效應(yīng)分析方案

本節(jié)闡述基于Maufory地形效應(yīng)量化表征模型的區(qū)域地形效應(yīng)計(jì)算方案［16］，具體包含兩部分：（1）基礎(chǔ)數(shù)據(jù)的獲取。（2）地形放大倍數(shù)的計(jì)算。地形效應(yīng)計(jì)算所需基礎(chǔ)數(shù)據(jù)包含高程數(shù)據(jù)以及場點(diǎn)VS。就全球范圍而言，30 m和90 m分辨率高程數(shù)據(jù)可從公開的數(shù)字高程Digital Elevation Model（DEM）模型中提取得到，如覆蓋了北緯83°到南緯83°之間的所有陸地區(qū)域的Advanced Spaceborne Thermal Emission and Reflection Radiometer（ASTER）Global Digital Elevation Model（Version 3）數(shù)字高程模型，其下載網(wǎng)址為https：//search.earthdata.nasa.gov/search。12.5m分辨率的數(shù)據(jù)也可以直接從網(wǎng)頁獲取，但一般需要付費(fèi)。場點(diǎn)VS可以從區(qū)域公共波速結(jié)構(gòu)模型中提取。伴隨能源物探工作的持續(xù)推進(jìn)和大中城市防震減災(zāi)工作的陸續(xù)開展，目前加州、日本和臺灣等地區(qū)已有若干公共波速結(jié)構(gòu)模型，如加州的Unified Community Velocity Model（UCVM）、日本的Japan Integrated Velocity Structure Model（JIVSM）和臺灣的Taiwan Velocity Model（TVM）。近年來，國內(nèi)公共波速模型也發(fā)展迅速，如段永紅等［17］利用人工地震寬角反射/折射探測剖面波速數(shù)據(jù)建立的華北克拉通中東部地殼三維速度結(jié)構(gòu)模型（HBCrust1.0）；姚華建［18］利用體波和面波走時(shí)數(shù)據(jù)反演成像得到的中國地震科學(xué)實(shí)驗(yàn)場地區(qū)地殼剪切波速模型，South West China Community Velocity Model version 1.0（SWChinaCVM-1.0）；HAN等［19］利用中國大陸數(shù)字化地震臺網(wǎng)記錄的區(qū)域地震走時(shí)數(shù)據(jù)成像得到的中國大陸巖石圈深達(dá)150 km的高分辨率三維地震縱波和橫波速度結(jié)構(gòu)模型，Unified Seismic Tomography Models for Continental China Lithosphere version 2.0（USTClitho2.0）。圖1給出了USTClitho2.0模型川滇地區(qū)（97.5°～105.5°E，21°～32°N）圖1（a）與SWChinaCVM-1.0模型圖1（b）表層VS分布圖及二者之間的差異圖1（c），其中：圖1（c）為USTClitho2.0模型的VS值減去SWChinaCVM-1.0的VS值后的差值圖像?？梢钥闯觯篣STClitho2.0模型的VS值普遍小于SWChinaCVM-1.0模型，兩類模型在圖示東南區(qū)域差異顯著，最大值差值為0.81 km/s。由于USTClitho2.0模型構(gòu)建過程中考慮了地形對成像結(jié)果的影響，本文選用USTClitho2.0確定場點(diǎn)VS。

圖1 川滇區(qū)域表層波速結(jié)構(gòu)Fig.1 Top shear wave velocity structure of different CVM model in Chuan-Dian region

結(jié)合數(shù)字高程和場點(diǎn)VS，即可利用Maufory模型計(jì)算區(qū)域任一點(diǎn)地形效應(yīng)，具體計(jì)算流程見上一節(jié)。然而，由圖1可知：受公共波速模型分辨率的影響，區(qū)域交接位置地形場點(diǎn)VS取值不唯一。就此類情形，參照交錯(cuò)網(wǎng)格有限差分法中介質(zhì)分界面處等效介質(zhì)參數(shù)的建立［20］，定義場點(diǎn)VS(xi,yi)為：

圖2為利用基礎(chǔ)資料和式（7）計(jì)算所得的川滇某一局部地區(qū)（3 km×3 km）對不同頻率地震動的放大倍數(shù)。經(jīng)緯度與大地坐標(biāo)系的轉(zhuǎn)換選用WGS84地理坐標(biāo)系。圖2（b）和圖2（d）分別給出了地形對1.15 Hz和8.06 Hz地震動傅里葉幅值譜的放大倍數(shù)AF50th(f)的空間分布圖。從圖中可以看出：Maufroy構(gòu)建的地形效應(yīng)定量表征模型，整體呈現(xiàn)出以下特征：凸起地形和凹陷地形分別對地震動存在放大和抑制作用，且放大倍數(shù)空間分布隨考察地震動的頻率降低而趨于平滑；放大區(qū)域集中在山頂與山脊的凸起地區(qū)，縮小區(qū)域集中在山底與山谷等凹陷地區(qū)；山頂與山脊地形對高頻地震動的平均放大系數(shù)在1.5左右，山底與山谷地形對高頻地震動的平均縮小系數(shù)在0.5左右。

圖2 利用Maufroy模型對局部地形效應(yīng)的估計(jì)Fig.2 Topographic amplification factors computed using Maufroy’s model

3 區(qū)域地形效應(yīng)分析方案在川滇區(qū)域的應(yīng)用

本節(jié)闡述區(qū)域地形效應(yīng)分析方案在川滇區(qū)域的應(yīng)用。川滇是中國首個(gè)研究“從地震破裂過程到工程結(jié)構(gòu)響應(yīng)”全鏈條的地震科學(xué)實(shí)驗(yàn)場的建設(shè)區(qū)域［21］。研究川滇地形對地震動的影響，對實(shí)驗(yàn)場全鏈條地震動模擬與預(yù)測平臺的構(gòu)建具有一定參考價(jià)值。本文根據(jù)川滇建設(shè)區(qū)域范圍97.5°～105.5°E，21°～32°N，按1°×1°分區(qū)下載地形圖。下面以30°～31°N，102°～103°E地區(qū)分塊為例，闡述川滇區(qū)域地形效應(yīng)分析。

首先利用USTClitho2.0和式（7）確定每個(gè)地形場點(diǎn)計(jì)算所需的VS，該地區(qū)VS為2.90 km/s。結(jié)合下載數(shù)字地形圖的最高分辨率確定h為30 m。在此基礎(chǔ)上，針對n的不同取值，確定地震動放大倍數(shù)的計(jì)算頻率點(diǎn)，結(jié)果見表1。

表1 不同n取值下對對應(yīng)的特征長度與計(jì)算頻率Table 1 Computational characteristic lengths and frequencies corresponding to different values of n

其次，計(jì)算地形場點(diǎn)對應(yīng)的頻率標(biāo)度曲率。用Global Mapper導(dǎo)入數(shù)字高程數(shù)據(jù)，利用轉(zhuǎn)換工具將場點(diǎn)經(jīng)緯度轉(zhuǎn)換為場點(diǎn)WGS84大地坐標(biāo)。由于經(jīng)過坐標(biāo)轉(zhuǎn)換后得到的地形圖為非規(guī)則矩形，為克服非規(guī)則矩形邊界點(diǎn)處數(shù)據(jù)急劇變化導(dǎo)致的邊界頻率標(biāo)度曲率計(jì)算結(jié)果不準(zhǔn)確的問題，本文通過截?cái)嗵幚韺⒎且?guī)則矩形轉(zhuǎn)化為規(guī)則矩形（見圖3）。重復(fù)上述操作直至川滇區(qū)域被所有矩形塊覆蓋。在此基礎(chǔ)上，計(jì)算表1中LS頻率標(biāo)度曲率CS。

圖3 下載的高程空間分布圖和經(jīng)過截?cái)嗵幚砗蟮母叱炭臻g分布圖Fig.3 Spatial distribution of the original elevations and the truncated distribution of elevations

然后，將CS代入式（1）、式（2）和式（3），即可得到不同頻率點(diǎn)處區(qū)域地形對不同頻率地震動的放大倍數(shù)。由圖4給出了頻率8.06 Hz、4.84 Hz、2.69 Hz和1.15 Hz處不同分位數(shù)處放大倍數(shù)空間分布，表2中給出了區(qū)域放大倍數(shù)計(jì)算結(jié)果的統(tǒng)計(jì)特征值。注意到實(shí)際計(jì)算所得放大倍數(shù)范圍中包含負(fù)值，這可能是由于一些場點(diǎn)位于極端凹陷的區(qū)域，這些場點(diǎn)曲率絕對值為較大負(fù)值，這一負(fù)值可能導(dǎo)致放大倍數(shù)小于0。為驗(yàn)證這一猜想，計(jì)算了圖4某一局部地區(qū)三維地形中紅色圓形區(qū)域內(nèi)（凹陷區(qū)）最底部場點(diǎn)在不同頻率下的AF50th：（1）8.06 Hz時(shí)為-1.39。（2）4.84 Hz時(shí)為-0.36。（3）2.69 Hz時(shí)為0.62。（4）1.15 Hz時(shí)為0.81。由此可知：當(dāng)使用該定量表征模型計(jì)算某一場點(diǎn)不同頻率下的放大倍數(shù)時(shí)，應(yīng)注意這個(gè)場點(diǎn)所位于的地形是否存在異常。反之，利用定量表征模型亦可以搜尋特殊的地形，例如尋找圖5所示的凹陷較深的地形，只需要找到放大倍數(shù)小于零的場點(diǎn)即可。

圖5 存在顯著凹陷的局部三維地形Fig.5 Local 3Dtopography with significant depressions

表2 30°～31°N，102°～103°E區(qū)域內(nèi)不同頻率地震動對應(yīng)的不同分位數(shù)地形放大倍數(shù)的統(tǒng)計(jì)特性Table 2 Statistics on topographic amplification factors for ground motion of different frequencies for different percentiles in 30°～31°N，102°～103°E region

圖4 不同頻率地震動對應(yīng)的不同分位數(shù)地形放大倍數(shù)Fig.4 Topographic amplification factors for ground motion of different frequencies for different percentiles

最后，在得到8.06 Hz至1.15 Hz的10個(gè)頻率下的放大倍數(shù)后進(jìn)行曲線擬合，即可得到1～8 Hz內(nèi)任意頻率下的三類放大倍數(shù)。將上述工作對剩余221塊區(qū)域重復(fù)進(jìn)行，即可得到川滇區(qū)域在30 m精度下的1～8 Hz的所有放大倍數(shù)值。表3給出了川滇地區(qū)（97.5°～105.5°E，21°～32°N）1～8 Hz的不同分位數(shù)地形放大倍數(shù)的統(tǒng)計(jì)特性。由表3可以看出：（1）同一分位數(shù)放大倍數(shù)隨頻率增加而變大，放大倍數(shù)分布范圍變寬，局部區(qū)域放大倍數(shù)最高可至5.69。（2）同一分位數(shù)倍數(shù)在不同頻率的均值是一樣的，且均值和中值近似相等。不同頻率下的AF16th的均值和中值為0.7，AF50th的均值和中值為1.0，AF84th的均值和中值為1.4。這三個(gè)值分別與公式（1）、式（2）和式（3）的截距一樣。對比AF16th，AF50th和AF84th與《建筑與市政工程抗震通用規(guī)范》（GB55002-2021）4.1.1條給出的不利地形地段處水平設(shè)計(jì)地震參數(shù)放大系數(shù)（1.1～1.6）可知：AF16th和AF50th的中值和均值低于規(guī)范給出的最小放大倍數(shù)1.1，雖然AF84th的中值和均值大于最小放大倍數(shù)1.1但卻低于規(guī)范給出的最大放大倍數(shù)1.6。因此就平均意義而言，規(guī)范取值略微偏于保守，規(guī)范取值將低估顯著凸起地形對高頻地震動的放大作用，如圖4所示。（3）S波頻率一定時(shí)，AF84th的方差大于AF50th，AF50th的方差大于AF16th。（4）由于統(tǒng)計(jì)所得偏度特征值均大于零，均值略大于中值可知川滇區(qū)域凸起地形覆蓋范圍略低于凹陷地形覆蓋范圍，但凸起地形的凸起程度高于凹陷地形的凹陷程度。值得注意的是：偏度大于0部分也有可能是文章未計(jì)入放大系數(shù)小于0場點(diǎn)所導(dǎo)致的。

表3 川滇區(qū)域不同頻率地震動對應(yīng)的不同分位數(shù)地形放大倍數(shù)的統(tǒng)計(jì)特性Table 3 Statistics on topographic amplification factors for ground motion of different frequencies for different percentiles in Chuan-Dian region

進(jìn)一步，本文結(jié)合漾濞地震震害考察結(jié)果分析川滇地形效應(yīng)放大倍數(shù)與震害的對應(yīng)。2021年5月21日漾濞縣發(fā)生6.4級地震。因?yàn)檠h城范圍較小，為了保證精度使用了分辨率較高的12.5 m的高程數(shù)字模型。計(jì)算漾濞縣在1～8Hz頻率下的AF50th。為了全面考慮各個(gè)頻率對地形放大的影響，漾濞地區(qū)每一點(diǎn)的AF50th取由8個(gè)頻率計(jì)算放大倍數(shù)中的最大值。計(jì)算結(jié)果如圖6（a）所示，并在圖6（b）給出漾濞地區(qū)建筑分布及對應(yīng)的震害分布圖。結(jié)合圖6（a）和圖6（b）可以觀察到：圖6（b）中山體突出部分（紅色矩形區(qū)域處）AF50th最高；漾濞縣城南部（橢圓區(qū)域）的放大倍數(shù)比縣城內(nèi)部的放大倍數(shù)高。這表明在一次地震動中，漾濞縣城南部可能遭受由地形引起的更大的地震動反應(yīng)，因此可能遭受更嚴(yán)重的破壞。結(jié)合圖6（a）和圖6（b）可以看出受災(zāi)情況與使用AF50th分布預(yù)測的災(zāi)害分布較為一致，即縣城南側(cè)災(zāi)害最為嚴(yán)重。這一結(jié)果初步驗(yàn)證了川滇地形效應(yīng)計(jì)算結(jié)果的可靠性；同時(shí)表明在強(qiáng)地震情景再現(xiàn)中考慮地形效應(yīng)的必要性。

圖6 漾濞地區(qū)AF50th圖像和漾濞地區(qū)震害分布Fig.6 AF50th of the Yangbi area and earthquake survey result in Yangbi area

4 川滇區(qū)域地形效應(yīng)分析結(jié)果的局限及應(yīng)用

本節(jié)結(jié)合地形效應(yīng)定量表征模型的構(gòu)建和震害調(diào)查結(jié)果，闡述區(qū)域地形效應(yīng)計(jì)算結(jié)果的局限性，并討論計(jì)算結(jié)果在工程常用寬頻帶地震動模擬及預(yù)測中的應(yīng)用。

本文采用的地形效應(yīng)量化表征模型存在如下局限：（1）模型僅考慮了S波段的放大，未計(jì)入地形對P波段以及S波段后續(xù)面波及尾波段的影響。（2）模型未考慮地震方位的影響，如多次震害表明面向震源方向山體坡面對地震動的放大顯著高于背向震源方向的山體坡面。（3）模型未明確地形場點(diǎn)VS的選用準(zhǔn)則，本文采用的低分辨率公共波速模型上層介質(zhì)的VS值，這一取值相當(dāng)于幾公里范圍深部地下波速的平均值，這可能與場地基巖的VS值不對應(yīng)，后續(xù)有必要深入研究VS選取對計(jì)算結(jié)果的影響。（4）模型假定地形效應(yīng)獨(dú)立于與震源、路徑、場地地震動效應(yīng)，即四者可以進(jìn)行解耦，但這一解耦是否成立有待進(jìn)一步論證。WANG等［13］認(rèn)為地形效應(yīng)模型構(gòu)建必須考慮淺表土層阻尼與厚度信息。（5）模型僅考慮了地形對高頻地震動傅里葉幅值譜的放大效應(yīng)，這與RAI等［14］為ASK14地震動衰減關(guān)系構(gòu)建的場地效應(yīng)修正模型不一致，RAI修正模型表明地形對低頻地震動反應(yīng)譜同樣存在顯著影響。當(dāng)然該不一致跟ASK14地震動衰減關(guān)系是否隱含納入了地形效應(yīng)的影響有關(guān)。（6）模型僅考慮了水平向地震動的地形效應(yīng)，未考慮地形對豎向地震動的影響。

忽略上述局限，下面給出區(qū)域地形效應(yīng)計(jì)算結(jié)果在工程常用寬頻帶地震動模擬與預(yù)測中應(yīng)用的初步建議：（1）計(jì)算模擬與預(yù)測區(qū)域范圍內(nèi)目標(biāo)場點(diǎn)處的（1～8 Hz）范圍地震動幅值譜的放大倍數(shù)AF50th，超過8 Hz的地震動幅值譜放大倍數(shù)直接采用8 Hz處的地震動幅值譜放大倍數(shù)，低于1 Hz地震動幅值譜放大倍數(shù)區(qū)為1。（2）利用工程常用寬頻帶地震動模擬與預(yù)測方法計(jì)算目標(biāo)場點(diǎn)地震動，根據(jù)計(jì)算所得地震動幅值譜放大倍數(shù)修正計(jì)算所得兩個(gè)水平方向地震動傅里葉幅值譜。（3）忽略地形對水平向地震動相位的影響，對修正后水平方向地震動傅里葉幅值譜做傅里葉逆變換得到考慮高頻地震動地形效應(yīng)的水平方向地震動。

5 結(jié)論

本文首先闡述了Maufroy表征模型以及地形放大因子的計(jì)算流程。在此基礎(chǔ)上，構(gòu)建基于USTClitho2.0波速模型和ASTGTMV003數(shù)字高程模型的區(qū)域地形效應(yīng)分析方案。以30°～31°N，102°～103°E區(qū)域?yàn)槔瑢?shí)現(xiàn)了對川滇區(qū)域地形效應(yīng)的整體分析。其結(jié)果表明：凸起地形對高頻地震動存在較為顯著的放大，而凹陷地形對高頻地震動存在抑制作用；川滇區(qū)域凸起地形覆蓋范圍略低于凹陷地形覆蓋范圍，但凸起地形的凸起程度高于凹陷地形的凹陷程度。此外，結(jié)合漾濞地震震害考察結(jié)果初步闡明了川滇地形效應(yīng)分析結(jié)果的可靠性以及山區(qū)建筑震害與地形效應(yīng)存在關(guān)聯(lián)，驗(yàn)證了凸起地形對高頻地震動的放大是導(dǎo)致山區(qū)建筑震害加劇的原因之一。最后，討論了文中地形效應(yīng)分析所用模型的局限性及模型在基于地震動衰減關(guān)系或混合寬頻帶地震動模擬方法開展的地震動模擬與預(yù)測工作中的應(yīng)用。本文工作還相當(dāng)粗略，后續(xù)擬針對量化效應(yīng)表征模型的局限性開展進(jìn)一步的研究與實(shí)踐工作。

致謝：

感謝中國地震局工程力學(xué)研究所林旭川研究員為本文提供的圖6b。感謝中國地震局工程力學(xué)研究所王輔臣博士提供的若干文獻(xiàn)。