３７６
吉林大学学报（地球科学版）
［４］Ａｋｉ
Ｋ，ＲｉｅｈａｒｄｓＰ
２００３妊
精确内插公式，并将其与显式有限元的内节点波动计算公式相结合，发展了计算水平成层介质自由场波动的时域数值模拟方法。该方法不仅具有与散射场的数值模拟相匹配、计算量与介质分层情况无关、占用计算机存储空间小的优点，而且为近场波动数值模拟中自由场与散射场的同步计算提供了可能。
本文是在中国地震局工程力学研究所廖振鹏院士的指导下完成的，在此深表感谢。
参考文献：
［１］廖振鹏．近场波动的有限元模拟［Ｊ］．地震工程与工程
振动．１９８４，４（２）：ｌ—１４．
Ｌｉａｏｗａｖｅ
Ｇ．定量地震学．第１卷［Ｍ］．李钦
祖，邹其嘉，等译．北京：地震出版社．１９８７．１３５—１６８
ＡｋｉＫ，ＲｉｃｈａｒｄｓＰ
Ｇ．ＱｕａｎｔｉｔａｔｉｖｃＳｅｉｓｍｏｌｏｇｙ．Ｖｏｌ
１
［Ｍ］Ｂｅｒｉｎｇ：ＣｈｉｎａＳｅｉｓｍｏｌｏｇｉｃａｌ
１６８．
Ｐｒｅｓｓ，１９８７．１３５一
［５］ＣｈｅｎＨｏｕｑｕｎ，Ｄｕ
Ｘｉｕｌｉ．Ｈ０ｕＳｈｕｎｚａｉ
ｔｏ
Ａｐｐｌｉｃａｔｉｏｎ
ｏｆ
ｔｒａｎｓｍｉｔｔｉｎｇｂｏｕｎｄａｒｉｅｓ
ｏｆ
ａｎ
ｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃａｎａｌｙｓｉｓ
ａｒｃｈ
ｄａｍ－ｆｏｕｎｄａｔｉｏｎ－ｒｅｓｅｒｖｏｉｒｓｙｓｔｅｍ［Ａ］．Ｉｎ：
ｓｏｉｌ－ｓｔｒｕｃｔｕｒｅｉｎ
ＺｈａｎｇＣ
Ｈ，ＷｏｌｆＪＰ，ｅｄｓ．Ｄｙｎａｍｉｃ
ｔｅｒａｃｒｉｏｎ（ｃｕｒｒｅｎｔｒｅｓｅａｒｃｈｉｎＣｈｉｎａａｎｄＳｗｉｔｚｅｒｌａｎｄ）
［Ｃ］．Ｂｅｉｊｉｎｇ：Ｉｎｔｅｒｎａｔｉｏｎａｌ
１９９７．１１５—１２４．
ＡｃａｄｅｍｉｃＰｕｂｌｉｓｈｅｒｓ，
［６］李山有，廖振鹏．地震体波斜入射情形下台阶地形引起
的波型转换［Ｊ］．地震工程与工程振动，２００２，２２（４）：９—
１５．Ｌｉ
Ｓｈａｎｙｏｕ，Ｌｉａｏ
ａ
Ｚｈｅｎｐｅｎｇ．Ｆｉｎｉｔｅｅｌｃｍｅｎｔｓｉｍｕｌａｔｉｏｎｏｆｎｅａｒ—ｆｉｅｌｄ
ｍｏｔｉｏｎ［Ｊ］．ＥａｒｔｈｑｕａｋｅＥｎｇｉｎｅｅｒｉｎｇ
Ｖｉｂｒａｔｉｏｎ，１９８４，４（２）：１—１４．
ａｎｄ
Ｅｎｇｉ—
ｎｅｅｒｉｎｇ
Ｚｈｅｎｐｅｎｇ．Ｗａｖｅ－ｔｙｐｅ
ｔｏ
ｃｏｎｖｅｒｓｉｏｎ
［２３廖振鹏，黄孔亮，杨柏坡，等．暂态波分析的透射边界条
件口］．中国科学（Ａ辑）。１９８４，２７：ｌ
ＬｉａｏＺｈｅｎｐｅｎｇ，Ｈｕａｎｇ
０６３一１
０７６．
ｃａｕｓｅｄｂｙ
ｓｔｅｐ
ｔｏｐｏｇｒａｐｈｙｓｕｂｊｅｃｔｅｄ
ｉｎｃｌｉｎｅｄｓｅｉｓ—
ｍｉｃｂｏｄｙｗａｖｅ［Ｊ］．ＥａｒｔｈｑｕａｋｅＥｎｇｉｎｅｅｒｉｎｇａｎｄＥｎｇｉ—
ｎｅｅｒｉｎｇ
Ｋｏｎｇｌｉａｎｇ，ＹａｎｇＢａｉｐｏ．ｅｔａ１．
Ｖｉｂｒａｔｉｏｎ．２００２，２２（４）：９—１５．
Ａｔｒａｎｓｍｉｔｔｉｎｇｂｏｕｎｄａｒｙｆｏｒｔｒａｎｓｉｅｎｔ％ｖａｖｅａｎａｌｙｓｅｓ
［７］李山有．重大工程结构的设计地震输人［Ｄ］．哈尔滨：
中国地震局工程力学研究所，２０００．５７～７４．
ＬｉＳｈａｎｙｏｕ．Ｄｅｓｉｇｎｅａｒｔｈｑｕａｋｅｉｎｐｕｔｆｏｒｉｍｐｏｒｔａｎｔ
ｇｉｎｅｅｒｉｎｇｉｎｇ
ｅｎ—
［Ｊ］．Ｓｃｉｅｎｃｅ
１０７６．
ｉｎＣｈｉｎａ（ＳｅｒｉｅｓＡ），１９８４，２７：１０８３—
［３３廖振鹏．工程波动理论导论第二版［Ｍ］．北京：科学出
版社，２００２．１３６—２８５．
Ｌｉａｏｉｎ
ｗａｖｅ
ｓｔｒｕｃｔｕｒｅ［Ｄ］．Ｈａｒｂｉｎ：Ｉｎｓｔｉｔｕｔｅ
ｏｆＥｎｇｉｎｅｅｒ—
５７
Ｍｅｃｈａｎｉｃｓ ＣｈｉｎａＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕ，２０００
Ｚｈｅｎｐｅｎｇ．Ｉｎｔｒｏｄｕｃｔｉｏｎ
ｔｏ
ｍｏｔｉｏｎｔｈｅｏｒｉｅｓ
Ｐｒｅｓｓ，
—７４
ｅｎｇｉｎｅｅｒｉｎｇ［Ｍ］，Ｂｅｒｉｎｇ；ＣｈｉｎａＡｃａｄｅｍｉｃ
２００２．１３６—２８５．
ＴＨＥ
ＴＩＭＥ—ＳＴＥＰＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴｌ０Ｎ
ｏＦＦＲＥＥＦＩＥＬＤ
ＭｏＴＩｏＮ
ｏＦ
ＬＡＹＥＲＥＤＨＡＬＦ二ＳＰＡＣＥＦＯＲＩＮＣＬＩＮＥＤ
ＬＩ
（１
ＩｎｓｔｉｔｕｔｅｏＪＥｎｇｉｎｅｅＨｎｇ
ＳＥＩｓＭＩＣＷＡＶＥＳ
Ｓｈａｎ—ｙｏｕ‘，ＷＡＮＧ
Ｘｕｅｌｉａｎ９２，ＺＨＯＵＺｈｅｎｇ—ｈｕａ‘
１５００８０，Ｃｈｉｎａ；２
Ｃｈｕｘｉｏｎｇ
Ｍｅｃｈａｎｉｃｓ，ＣｈｉｎａＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕ，Ｈａｒｂｉｎ
ＮｏｒｍａｌＣｏｌｌｅｇｅｏＪＹｕｎｎａｎ
Ｐｒｏｖｉｎｃｅ．Ｃｈｕｘｉｏｎｇ
６７５０００，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｔｈｅ
ｉｎｎｅｒ－ｔｒａｎｓｍｉｔｔｉｎｇ
ｗａｖｅｍｏｔｉｏｎ，ｉｅ．，ｔｈｅｆｒｅｅ－ｆｉｅｌｄ
ｌＺｌＯＩｉｏｎｉｎ
ａ
ｌａｙｅｒｅｄ
ｈａｌｆｓｐａｃｅ
ｂａｃｋｇｒｏｕｎｄ
ｆａＵＳｔ
ｂｅｄｅｔｅｒｍｉｎｅｄｉｎｏｒｄｅｒｔｏｃａｌｃｕｌａｔｅｔｈｅｍｏｔｉｏｎｏｆｂｏｕｎｄａｒｙｎｏｄｅｓａｎｄｉｎｐｕｔｗａｖｅｍｏ—
ｓｃａｔｔｅｒｉｎｇｐｒｏｂｌｅｍ．Ａ
ｔｉｏｎｅｎｅｒｇｙ，ｗｈｅｎａｐｐｌｙｉｎｇｆｉｎｉｔｅｅｌｅｍｅｎｔｉｎｔｈｅｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｗａｖｅｍｏｔｉｏｎ
ｔｉｍｅｄｏｍａｉｎｅｑｕａｔｉｏｎ
ｓｅｔ
ｉｓｐｒｅｓｅｎｔｅｄｆｏｒｃａｌｃｕｌａｔｉｎｇｆｒｅｅ—ｆｉｅｌｄｍｏｔｉｏｎｏｆｌａｙｅｒｅｄｈａｌｆ－ｓｐａｃｅ，ｔａｋｉｎｇｓｉｄｅ
ｂｏｕｎｄａｒｙｎｏｄｅｓｆｏｒｔｈｅｅｘａｍｐｌｅ．Ｔｈｅｅｑｕａｔｉｏｎｓｅｔｔａｋｅｓｔｈｅｃｏｎｓｉｄｅｒａｔｉｏｎｏｆｄｅｃｏｕｐｌｉｎｇｍｏｔｉｏｎｅｑｕａｔｉｏｎ
ｏｆｆｉｎｉｔｅｅｌｅｎｌｅｎｔａｎｄｔｈｅｒｅｌａｔｉｏｎｓｌｆｏｒｍｕｌａｅｂｅｔｗｅｅｎｎｅｉｇｈｂｏｒｉｎｇｎｏｄｅｓｂａｓｅｄ
ｔｉｎｇｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆ
ｓｅｎｓｅ
ａ
ｏｎ
ｔｈｅｈｏｒｉｚｏｎｔａｌｐｒｏｐａｇａ—
ｅｘａｃｔ
ｆｒｅｅ—ｆｉｅｌｄｍｏｔｉｏｎｓｙｓｔｅｍ．Ｔｈｅｍｅｔｈｏｄｇｉｖｅｓ
ｃａｎ
ａ
ｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｉｎｔｈｅ
ｏｆｆｉｎｉｔｅｅｌｅｍｅｎｔ，ｗｈｉｃｈ
ｂｅａｐｐｌｉｅｄ
ｔｏ
ｓｙｎｃｈｒｏｎｉｚｅｄｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｉｎｎｅｒ—ｔｒａｎｓｍｉｔｔｉｎｇａｎｄ
ｓｃａｔｔｅｒｉｎｇｗａｖｅ－ｍｏｔｉｏｎｆｉｅｌｄｉｎｔｉｍｅ—ｄｏｍａｉｎ．
Ｋｅｙｗｏｒｄｓ：ｌａｙｅｒｅｄｈａｌｆｓｐａｃｅ；ｉｎｃｌｉｎｅｄｓｅｉｓｍｉｃｗａｖｅ；ｆｒｅｅ—ｆｉｅｌｄｍｏｔｉｏｎ；ｔｉｍｅ—ｄｏｍａｉｎａｌｇｏｒｉｔｈｍ