趙志新 徐紀(jì)人

（中國(guó)地質(zhì)科學(xué)院地質(zhì)研究所，北京 100037）

地震波動(dòng)數(shù)值模擬方法研究非平勻介質(zhì)中地震波傳波問(wèn)題不僅是十分重要的理論課題，而且應(yīng)用也極為廣泛。

本文用錯(cuò)格實(shí)數(shù)傅立葉變換微分(SGRFFTD)擬譜法的數(shù)值模擬求解波動(dòng)方程，分析了地震波在沖積扇，盆地域中等多種不平勻地震構(gòu)造體的地球表層區(qū)域的傳播過(guò)程和地面運(yùn)動(dòng)分布。該方法的運(yùn)算速度比錯(cuò)格復(fù)數(shù)傅里葉變換微分法快1/3。而且在整個(gè)錯(cuò)格實(shí)數(shù)傅里葉變換微分法運(yùn)算過(guò)程中保留了奈奎斯特譜分量.使得新方法比普通分格的實(shí)數(shù)傅里葉變換微分法快精度高，穩(wěn)定性好，以新方法與De-hoop解析法在求解無(wú)限半空時(shí)地震波動(dòng)問(wèn)題中進(jìn)行比較表明，新數(shù)值解法的精度是令人滿(mǎn)意的。研究結(jié)果表明：

（1）當(dāng)入射波S1，在地球表面與土層底部形成多重反射。 由于每次上下反射波的反射與入射面的角度不盡相同，S1在不同地方振幅，頻率不同。而且由S1的衍生的反射波S2，S3……的振幅不斷發(fā)生變化，S2，S3在不同地方振幅、頻率不同，隨地而異。

（2）多重震蕩反射；例子中的多次反射多達(dá)5次。多次反射波S2，S3與S5的周期(p)一般要比初至波S1的周期大。因此在同速介質(zhì)體中，反射波的波長(zhǎng)也要長(zhǎng)些。有些地方反射波的振幅比入射波振幅大。

（3）空時(shí)各處遭受到反射次數(shù)的多少，反射傳播特征與地質(zhì)結(jié)構(gòu)有關(guān)。

這些特征對(duì)于識(shí)別多重反射波有積極意義。多次反射波的頻率和振幅變化，在強(qiáng)地面運(yùn)動(dòng)地震學(xué)研究中，對(duì)于研究不同頻率的波與不同建筑物造成不同破壞之時(shí)的關(guān)系、研究吸收介質(zhì)特征都有著積極的作用。再者不同譜特征將會(huì)應(yīng)用于勘探地球物中識(shí)別多重反射波。特別是對(duì)于在反射勘探中消除海洋地震反射勘探中的多次反射波干擾這一前沿學(xué)科研究有良好的應(yīng)用前景。

波的多重干涉與地下結(jié)構(gòu)有關(guān)。詳細(xì)調(diào)查，分析地質(zhì)結(jié)構(gòu)對(duì)于用地面運(yùn)動(dòng)數(shù)值模擬研究，對(duì)于地震波傳播和抗震研究都是十分重要的。