Ｅａｒｔｈｑｕａｋｅ Ｒｅｓｅａｒｃｈ ｉｎ Ｃｈｉｎａ
Ｖｏｌｕｍｅ ３０， Ｎｕｍｂｅｒ ３， ２０１６
Ｕｐｐｅｒ Ｃｒｕｓｔａｌ Ｖｅｌｏｃｉｔｙ Ｓｔｒｕｃｔｕｒｅ
ａｌｏｎｇ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｆｒｏｍ
Ｍａ̓ａｎｓｈａｎ ｔｏ Ａｎｑｉｎｇ１
Ｔｉａｎ Ｘｉａｏｆｅｎｇ， Ｗａｎｇ Ｆｕｙｕｎ， Ｌｉｕ Ｂａｏｆｅｎｇ， Ｙａｎｇ Ｚｈｕｏｘｉｎ， Ｚｈｅｎｇ Ｃｈｅｎｇｌｏｎｇ， ａｎｄ Ｇａｏ Ｚｈａｎｙｏｎｇ
Ｇｅｏｐｈｙｓｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ Ｃｅｎｔｅｒ， Ｃｈｉｎａ Ｅａｒｔｈｑｕａｋｅ Ａｄｍｉｎｉｓｔｒａｔｉｏｎ， Ｚｈｅｎｇｚｈｏｕ ４５０００３， Ｃｈｉｎａ
Ｗｅ ａｐｐｌｉｅｄ ｔｈｅ ３Ｄ ｆｉｒｓｔ ａｒｒｉｖａｌ ｔｒａｖｅｌ ｔｉｍｅ ｔｏｍｏｇｒａｐｈｙ ｍｅｔｈｏｄ ｔｏ ｔｈｅ Ａｎｈｕｉ ａｃｔｉｖｅ ｓｅｉｓｍｉｃ
ｓｏｕｒｃｅ ｅｘｐｅｒｉｍｅｎｔ ｄａｔａ， ａｎｄ ｏｂｔａｉｎｅｄ ｔｈｅ ｉｍａｇｉｎｇ ｏｆ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ ｖｅｌｏｃｉｔｙ ｓｔｒｕｃｔｕｒｅ
ｂｅｎｅａｔｈ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｆｒｏｍ Ｍａ̓ａｎｓｈａｎ， Ｔｏｎｇｌｉｎｇ ｔｏ Ａｎｑｉｎｇ． Ｄａｔａ ｆｉｔｔｉｎｇ ｒｅｖｅａｌｓ ｔｈｅ
ｔｏｍｏｇｒａｐｈｉｃ ｍｏｄｅｌ ｆｉｔｓ ｔｈｅ ｄａｔａ ｗｉｔｈ ｕｎｃｅｒｔａｉｎｔｉｅｓ， ｗｉｔｈｏｕｔ ｏｖｅｒｆｉｔｔｉｎｇ， ａｎｄ ｗｉｔｈ ａ
ｍｉｎｉｍｕｍ ｏｆ ｃｏｍｐｌｅｘｉｔｙ． Ｔｈｅ ｔｏｍｏｇｒａｐｈｉｃ ｒｅｓｕｌｔ ｓｈｏｗｓ ａｎ ｏｂｖｉｏｕｓ ｈｅｔｅｒｏｇｅｎｅｏｕｓ ｕｐｐｅｒ
ｃｒｕｓｔ ｗｈｉｃｈ ｃｏｎｓｉｓｔｓ ｏｆ ａ ｓｅｒｉｅｓ ｏｆ ｕｐｌｉｆｔｓ ａｎｄ ｄｅｐｒｅｓｓｉｏｎ ｂａｓｉｎｓ． Ｔｈｅ ｖｅｌｏｃｉｔｙ ｍｏｄｅｌ ａｎｄ
ｒｅｇｉｏｎ ｉｍｐｌｙ ｔｈａｔ ｔｈｉｓ ｒｅｇｉｏｎ ｈａｓ ｅｘｐｅｒｉｅｎｃｅｄ ｃｒｕｓｔａｌ ｕｐｌｉｆｔ ａｎｄ ｅｘｔｅｎｓｉｏｎａｌ ｔｅｃｔｏｎｉｓｍ ｗｉｔｈ
ｃｏｎｃｏｍｉｔａｎｔ ｍａｇｍａｔｉｓｍ ｓｉｎｃｅ ｔｈｅ Ｃｅｎｏｚｏｉｃ．
Ｋｅｙ ｗｏｒｄｓ： Ｓｅｉｓｍｉｃ ａｃｔｉｖｅ ｓｏｕｒｃｅ　Ａｉｒｇｕｎ　３Ｄ ｓｅｉｓｍｉｃ ｔｏｍｏｇｒａｐｈｙ　Ｓｕｒｖｅｙ ｇｅｏｍｅｔｒｙ
ＩＮＴＲＯＤＵＣＴＩＯＮ
Ｗｉｔｈ ｔｈｅ ｕｎｉｆｉｅｄ ｐｌａｔｆｏｒｍ ｏｆ “ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ Ｇｅｏｓｃｉｅｎｃｅ Ｐｒｏｊｅｃｔ” ， ｔｈｅ Ａｎｈｕｉ ａｒｔｉｆｉｃｉａｌ ｓｏｕｒｃｅ
ｓｃｉｅｎｔｉｆｉｃ ｅｘｐｅｒｉｍｅｎｔ ｏｆ ｕｎｄｅｒｇｒｏｕｎｄ ｓｔｒｕｃｔｕｒｅ ｅｘｐｌｏｒａｔｉｏｎ ｓｅｔ ｔｈｅ Ａｎｈｕｉ ｓｅｃｔｉｏｎ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ
Ｒｉｖｅｒ ａｓ ｔｈｅ ｅａｒｌｙ ｅｘｐｅｒｉｍｅｎｔａｌ ｚｏｎｅ ｔｏ ｅｘｐｌｏｒｅ ｔｈｅ ｓｕｂｓｕｒｆａｃｅ ｓｔｒｕｃｔｕｒｅ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ａｎｄ
ｏｂｔａｉｎ ｔｈｅ ｔｈｒｅｅ⁃ｄｉｍｅｎｓｉｏｎａｌ ｓｕｂｓｕｒｆａｃｅ ｓｔｒｕｃｔｕｒｅ ｍｏｄｅｌ ｂｙ ａｉｒｇｕｎ ｓｅｉｓｍｉｃ ｓｏｕｒｃｅ ｇｅｎｅｒａｔｉｏｎ ｉｎ ｔｈｅ
Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｗａｔｅｒｗａｙ ｗｈｅｒｅ ｐｅｒｍａｎｅｎｔ ｓｅｉｓｍｉｃ ｓｔａｔｉｏｎｓ ａｎｄ ａ ｔｈｒｅｅ⁃ｄｉｍｅｎｓｉｏｎａｌ ｍｏｂｉｌｅ
ｏｂｓｅｒｖａｔｉｏｎ ｓｙｓｔｅｍ ｗｅｒｅ ｓｅｔ ｕｐ ｔｏ ｒｅｃｅｉｖｅ ａｒｔｉｆｉｃｉａｌ ｓｅｉｓｍｉｃ ｗａｖｅ ｓｉｇｎａｌ ｗｉｔｈｉｎ ｔｈｅ ｓｃｏｐｅ ｏｆ ａｒｏｕｎｄ
３００ｋｍ， ｗｈｉｃｈ ｈｅｌｐ ｔｏ ｓｔｕｄｙ ｔｈｅ ｉｓｓｕｅｓ ｒｅｌａｔｅｄ ｔｏ ｔｈｅ ｒｅｇｉｏｎａｌ ｓｐｅｃｉａｌ ｔｅｃｔｏｎｉｃ ｅｎｖｉｒｏｎｍｅｎｔ，
ｆｏｒｍａｔｉｏｎ ｍｅｃｈａｎｉｓｍ ｏｆ ｔｈｅ ｍｅｔａｌｌｏｇｅｎｉｃ ｂｅｌｔ ａｎｄ ｔｈｅ ｅｎｄ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｇｉａｎｔ ｓｔｒｉｋｅ ｓｌｉｐ ｆａｕｌｔ２ ．
Ａｎｈｕｉ Ｐｒｏｖｉｎｃｅ ｉｓ ｌｏｃａｔｅｄ ｉｎ ｔｈｅ ｃｏｕｐｌｉｎｇ ｚｏｎｅ ｂｅｔｗｅｅｎ ｔｈｅ Ｎｏｒｔｈ Ｃｈｉｎａ ｂｌｏｃｋ ａｎｄ ｔｈｅ
Ｙａｎｇｔｚｅ ｂｌｏｃｋ， ｗｈｉｃｈ ｉｓ ｔｈｅ ｉｎｔｅｒｓｅｃｔｉｏｎ ｏｆ ｔｈｅ Ｑｉｎｌｉｎｇ⁃Ｄａｂｉｅ ｏｒｏｇｅｎｉｃ ｂｅｌｔ ａｎｄ Ｙａｎｇｔｚｅ ｂｌｏｃｋ． Ｉｎ
Ｒｅｃｅｉｖｅｄ ｏｎ Ａｐｒｉｌ １， ２０１６； ｒｅｖｉｓｅｄ ｏｎ Ｊｕｎｅ １６， ２０１６． Ｔｈｉｓ ｐｒｏｊｅｃｔ ｗａｓ ｊｏｉｎｔｌｙ ｓｐｏｎｓｏｒｅｄ ｂｙ ｔｈｅ Ｎａｔｉｏｎａｌ Ｎａｔｕｒａｌ
Ｓｃｉｅｎｃｅ Ｆｏｕｎｄａｔｉｏｎ ｏｆ Ｃｈｉｎａ （４１５７４０８４） ， ａｎｄ ｔｈｅ Ｓｐａｒｋ Ｐｒｏｇｒａｍ ｏｆ Ｅａｒｔｈｑｕａｋｅ Ｓｃｉｅｎｃｅｓ （ ＸＨ１５０５９） ．
２
Ｃｈｅｎ Ｔａｏ， ２０１５． Ｉｎｔｒｏｄｕｃｔｉｏｎ ｔｏ ｔｈｅ Ａｎｈｕｉ Ａｒｔｉｆｉｃｉａｌ Ｓｅｉｓｍｉｃ Ｓｏｕｒｃｅ Ｓｃｉｅｎｔｉｆｉｃ Ｅｘｐｅｒｉｍｅｎｔ ｏｆ Ｕｎｄｅｒｇｒｏｕｎｄ
Ｓｔｒｕｃｔｕｒｅ Ｅｘｐｌｏｒａｔｉｏｎ Ｅｘｐｅｒｉｍｅｎｔｓ ｆｏｒ ｔｈｅ Ｙａｎｇｔｚｅ ｒｉｖｅｒ， ｈｔｔｐ： ／ ／ ｗｗｗ． ｃｅａ． ｇｏｖ． ｃｎ ／ ｐｕｂｌｉｓｈ ／ ｄｉｚｈｅｎｊ ／ １２１ ／ ３７９ ／
２０１５０９１７１５４６１３３９７７３７０９０ ／ ｉｎｄｅｘ． ｈｔｍｌ
１
Ｖｏｌｕｍｅ ３０， Ｎｕｍｂｅｒ ３
４３１
ｔｈｅ ｆｏｒｅｌａｎｄ ｏｆ ｔｈｅ ｉｎｔｅｒｓｅｃｔｉｏｎ ｂｅｔｗｅｅｎ Ｑｉｎｌｉｎｇ⁃Ｄａｂｉｅ ｏｒｏｇｅｎｉｃ ｂｅｌｔ ａｎｄ ｔｈｅ Ｔａｎ⁃Ｌｕ（ Ｔａｎｃｈｅｎｇ⁃
Ｌｕｊｉａｎｇ） ｆａｕｌｔ ｚｏｎｅ ｃｕｔｔｉｎｇ ｔｈｅ ｃｒｕｓｔ， ｔｈｅｒｅ ｉｓ ａ Ｖ⁃ｓｈａｐｅｄ ｍｅｔａｌｌｏｇｅｎｉｃ ｂｅｌｔ， ｎａｒｒｏｗ ｉｎ ｔｈｅ ＳＷ
ａｎｄ ｂｒｏａｄ ｉｎ ｔｈｅ ＮＥ， ｗｈｉｃｈ ｉｓ ｃｅｎｔｅｒｅｄ ｂｙ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｉｎ ｔｈｅ ｒａｎｇｅ ｏｆ ｔｈｅ ５０ｋｍ － １００ｋｍ
ｉｎ ｔｈｅ ｎｏｒｔｈ ａｎｄ ｓｏｕｔｈ ｒｅｓｐｅｃｔｉｖｅｌｙ （ Ｃｈａｎｇ Ｙｉｎｆｏ ｅｔ ａｌ． ， １９９１ ） （ Ｆｉｇ． １） ． Ｗｉｔｈ ｔｈｅ ｉｎｔｒａ⁃
ｃｏｎｔｉｎｅｎｔａｌ ｏｒｏｇｅｎｙ ｆｒｏｍ ｔｈｅ ｒｅｍｏｔｅ ｅｘｔｒｕｓｉｏｎ ｆｏｒｃｅ ｃａｕｓｅｄ ｂｙ ｔｈｅ ＮＷ ｌｏｗ⁃ａｎｇｌｅ ｓｕｂｄｕｃｔｉｏｎ ｏｆ ｔｈｅ
ｐａｌｅｏ⁃Ｐａｃｉｆｉｃ ｐｌａｔｅ ｄｕｒｉｎｇ ｔｈｅ ｍｉｄｄｌｅ ｔｏ ｌａｔｅ Ｊｕｒａｓｓｉｃ⁃ｅａｒｌｙ Ｃｒｅｔａｃｅｏｕｓ ａｎｄ ｔｈｅ ｓｕｂｓｅｑｕｅｎｔ
ｃｏｎｃｏｍｉｔａｎｔ ｍａｓｓｉｖｅ ｖｏｌｃａｎｉｃ⁃ｉｎｔｒｕｓｉｖｅ ｍａｇｍａｔｉｓｍ ｅｘｔｅｎｓｉｏｎ， ａｎ ａｌｔｅｒｎａｔｉｖｅ ｔｅｃｔｏｎｉｃ ｆｒａｍｅｗｏｒｋ ｏｆ
ｕｐｌｉｆｔｓ （ ｍｉｎｅｒａｌ ｇａｔｈｅｒｉｎｇ ａｒｅａｓ ａｓ Ｔｏｎｇｌｉｎｇ， Ｎｉｎｇｚｈｅｎ， Ｇｕｉｃｈｉ ） ａｎｄ ｄｅｐｒｅｓｓｉｏｎｓ （ ｍｉｎｅｒａｌ
ｇａｔｈｅｒｉｎｇ ａｒｅａｓ ａｓ Ｎｉｎｇｗｕ， Ｌｕｚｏｎｇ） ｗａｓ ｆｏｒｍｅｄ ｉｎ ｔｈｅ ｍｅｔａｌｌｏｇｅｎｉｃ ｂｅｌｔ ｏｆ ｔｈｅ Ｍｉｄｄｌｅ ａｎｄ Ｌｏｗｅｒ
Ｒｅａｃｈｅｓ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ （ Ｃｈａｎｇ Ｙｉｎｆｏ ｅｔ ａｌ． ， １９９１； Ｌｕ Ｑｉｎｇｔｉａｎ ｅｔ ａｌ． ， ２０１４ ） ． Ｔｈｅ
ｎｏｒｔｈｗｅｓｔ ｐａｒｔ ｏｆ ｔｈｅ ｍｅｔａｌｌｏｇｅｎｉｃ ｂｅｌｔ ｉｓ ａｄｊａｃｅｎｔ ｔｏ Ｋｏｎｇｌｉｎｇ⁃Ｄｏｎｇｌｉｎｇ ｔｅｒｒａｎｅ ｉｎ ｔｈｅ Ａｎｑｉｎｇ ａｒｅａ
ｗｈｉｌｅ ｉｔｓ ｓｏｕｔｈｅａｓｔ ｉｓ ｃｏｎｎｅｃｔｅｄ ｔｏ ｔｈｅ ｓｏｕｔｈｅｒｎ ｔｅｒｒａｎｅ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ． Ｔｈｅｒｅｆｏｒｅ， ｔｈｅ
“ Ｐｏｌｙ⁃Ｂａｓｅｍｅｎｔ ｗｉｔｈ ｏｎｅ ｃｏｖｅｒ” ｔｅｃｔｏｎｉｃ ｐａｔｔｅｒｎ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ ｂｌｏｃｋ， （ Ｃｈａｎｇ Ｙｉｎｆｏ ｅｔ ａｌ． ，
１９９６） ｉｓ ｅｍｂｏｄｉｅｄ ｉｎ ｔｈｅ ｃｒｕｓｔ ｏｆ ｔｈｉｓ ｒｅｇｉｏｎ．
Ｆｉｒｓｔ⁃ａｒｒｉｖａｌ ｔｒａｖｅｌ ｔｉｍｅ ｔｏｍｏｇｒａｐｈｙ ｉｓ ｔｈｅ ｍａｉｎ ｍｅｔｈｏｄ ｔｏ ｏｂｔａｉｎ ｈｉｇｈ⁃ｒｅｓｏｌｕｔｉｏｎ ｃｒｕｓｔａｌ
ｖｅｌｏｃｉｔｙ ｓｔｒｕｃｔｕｒｅｓ （ Ｚｅｌｔ ｅｔ ａｌ． ， １９９８； ＤｕａｎＹｏｎｇｈｏｎｇｅｔ ａｌ． ， ２００２； Ｘｕ Ｚｈａｏｆａｎ ｅｔ ａｌ． ， ２００６；
Ｐａｎ Ｊｉｓｈｕｎ ｅｔ ａｌ． ２００８； Ｊｉａ Ｙｕｐｅｎｇ ｅｔ ａｌ． ， ２０１２ ） ． Ｔｈｅ ｈｉｇｈ⁃ｒｅｓｏｌｕｔｉｏｎ ｐｒｏｆｉｌｉｎｇ ａｌｏｎｇ ｔｈｅ
Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｏｂｔａｉｎｅｄ ｆｒｏｍ ｔｈｅ Ａｎｈｕｉ ｅｘｐｅｒｉｍｅｎｔ ｏｆ ｔｈｅ “ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ Ｇｅｏｓｃｉｅｎｃｅ Ｐｒｏｊｅｃｔ” ｈａｓ
ｒｅａｌｉｚｅｄ ｄｅｎｓｅ ｏｂｓｅｒｖａｔｉｏｎ ｏｎ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ ｏｆ ａｂｏｖｅ⁃ｍｅｎｔｉｏｎｅｄ ｍｅｔａｌｌｏｇｅｎｉｃ ｂｅｌｔ ａｎｄ ｔｈｅ
Ｋｏｎｇｌｉｎｇ⁃Ｄｏｎｇｌｉｎｇ ｔｅｒｒａｎｅ ｉｎ ｔｈｅ Ａｎｑｉｎｇ ａｒｅａ， ｗｈｉｃｈ ｐｒｏｖｉｄｅｓ ａ ｇｏｏｄ ｄａｔａｂａｓｅ ｆｏｒ ｒｅｓｅａｒｃｈ ｏｎ ｔｈｅ
ｒｅｇｉｏｎａｌ ｃｒｕｓｔａｌ ｓｔｒｕｃｔｕｒｅ， ｔｅｃｔｏｎｉｃ ｅｖｏｌｕｔｉｏｎ ａｎｄ ｍｉｎｅｒａｌｉｚａｔｉｏｎ ｍｅｃｈａｎｉｓｍ．
１　ＯＢＳＥＲＶＡＴＩＯＮ ＡＮＤ ＤＡＴＡ
Ｔｈｅ Ａｎｈｕｉ ｅｘｐｅｒｉｍｅｎｔ ｏｆ ｔｈｅ “ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ Ｇｅｏｓｃｉｅｎｃｅ Ｐｒｏｊｅｃｔ ” ｉｓ ｌｏｃａｔｅｄ ｉｎ Ａｎｈｕｉ
Ｐｒｏｖｉｎｃｅ． Ｗｉｔｈ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ａｓ ｔｈｅ ｃｅｎｔｅｒ， ｗｅ ｌａｉｄ ｏｕｔ ９ ｓｕｒｖｅｙ ｌｉｎｅｓ ｐａｒａｌｌｅｌ ｔｏ ｔｈｅ Ｙａｎｇｔｚｅ
Ｒｉｖｅｒ， ２ ｓｕｒｖｅｙ ｌｉｎｅｓ ｎｅａｒｌｙ ｐｅｒｐｅｎｄｉｃｕｌａｒ ｔｏ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ａｎｄ ｍｏｒｅ ｔｈａｎ ７００ ｔｈｒｅｅ⁃
ｃｏｍｐｏｎｅｎｔ ｄｉｇｉｔａｌ ｓｅｉｓｍｏｇｒａｐｈｓ， ｉｎｃｌｕｄｉｎｇ ３５０ ＰＤＳ⁃２ ｄｉｇｉｔａｌ ｓｅｉｓｍｏｇｒａｐｈｓ， ３５０ ＤＺＳ⁃１ ｄｅｅｐ
ｄｉｇｉｔａｌ ｓｅｉｓｍｏｇｒａｐｈ ａｎｄ ５０ ｔｅｓｔ ｉｎｓｔｒｕｍｅｎｔｓ ｗｉｔｈ ａ ｓａｍｐｌｉｎｇ ｒａｔｅ ｏｆ ２００ ｍａｄｅ ｂｙ Ｃｈｏｎｇｑｉｎｇ
Ｇｅｏｌｏｇｉｃａｌ Ｉｎｓｔｒｕｍｅｎｔ Ｆａｃｔｏｒｙ． Ｔｈｅ ｏｂｓｅｒｖａｔｉｏｎ ｓｙｓｔｅｍ ｉｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ． １． Ｔｈｅ ｃｏｎｔｉｎｕｏｕｓ
ｏｂｓｅｒｖａｔｉｏｎ ｉｓ ｃｏｎｄｕｃｔｅｄ ｉｎ ｍｏｂｉｌｅ ｏｂｓｅｒｖａｔｉｏｎ ｔｏ ｒｅｃｅｉｖｅ ｓｉｇｎａｌｓ ｇｅｎｅｒａｔｅｄ ｂｙ ｔｈｅ ａｉｒｇｕｎ ｉｎ ｔｈｅ
Ｙａｎｇｔｚｅ Ｒｉｖｅｒ． Ｉｎ ｏｒｄｅｒ ｔｏ ｅｎｓｕｒｅ ｔｈｅ ｏｂｓｅｒｖａｔｉｏｎ ｑｕａｌｉｔｙ， ａｌｌ ｓｅｉｓｍｏｇｒａｐｈｓ ａｒｅ ｐｌａｃｅｄ ｉｎ ｐｉｔｓ ｗｉｔｈ
ｄｉａｍｅｔｅｒ ｏｆ ｍｏｒｅ ｔｈａｎ ３０ｃｍ ａｎｄ ｔｈｅ ｄｅｐｔｈ ｏｆ ｍｏｒｅ ｔｈａｎ ４０ｃｍ， ａｎｄ ｔｈｅ ｐｌａｔｆｏｒｍｓ ａｒｅ ｓｅｔ ｕｐ ａｎｄ
ｗａｔｅｒｐｒｏｏｆｅｄ ｉｎ ｐａｒｔ ｏｆ ｔｈｅ ｒｅｇｉｏｎ．
Ｉｎ ｔｈｉｓ ｐａｐｅｒ， ｔｈｅ ｓｔｕｄｙ ｏｂｊｅｃｔ ｉｓ ｔｈｅ ｄａｔａ ｐｒｏｄｕｃｅｄ ｂｙ ２０ ｆｉｘｅｄ⁃ｐｏｉｎｔ ｅｘｃｉｔａｔｉｏｎｓ ｏｂｓｅｒｖｅｄ ｉｎ
ｔｈｅ ｓｕｒｖｅｙ ｌｉｎｅ ａｌｏｎｇ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ｗｈｅｒｅ ｔｈｅｒｅ ａｒｅ ａ ｔｏｔａｌ ｏｆ １００ ｉｎｓｔｒｕｍｅｎｔｓ ｄｅｐｌｏｙｅｄ ｗｉｔｈ
ｔｈｅ ａｖｅｒａｇｅ ｓｕｒｖｅｙ ｐｏｉｎｔ ｄｉｓｔａｎｃｅ ｏｆ ２ｋｍ． Ｔｈｅ ｗｅａｋ ｓｉｇｎａｌ ｉｓ ｅｘｔｒａｃｔｅｄ ｂｙ ｌｉｎｅａｒ ｓｔａｃｋｉｎｇ （ Ｗａｎｇ
ｅｔ ａｌ． ， ２０１２ ） ｏｆ ｔｈｅ ｄａｔａ ｇｅｎｅｒａｔｅｄ ｂｙ ｔｈｅ ２０ ｆｉｘｅｄ⁃ｐｏｉｎｔ ｅｘｃｉｔａｔｉｏｎｓ ａｎｄ ｆｉｎａｌｌｙ ２０００ ｓｔａｃｋ
ｒｅｃｏｒｄｓ ａｒｅ ｏｂｔａｉｎｅｄ． Ａ ｔｙｐｉｃａｌ ｓｔａｃｋ ｒｅｃｏｒｄ ｉｓ ｔｈｅ １６ ｔｈ ｆｉｘｅｄ⁃ｐｏｉｎｔ ｅｘｃｉｔａｔｉｏｎ （ Ｆｉｇ． ２ ） ． Ａｌｌ
ｒｅｃｏｒｄｓ ａｒｅ ３ － ８Ｈｚ ｂａｎｄｐａｓｓ ｆｉｌｔｅｒｅｄ， ｔｈｅｎ ｔｈｅ ｆｉｒｓｔ ａｒｒｉｖａｌ Ｐｇ ｐｈａｓｅ ｉｓ ｐｉｃｋｅｄ ｕｐ ａｎｄ ｔｈｅ ｓｅｉｓｍｉｃ
ｐｈａｓｅ ｐｉｃｋｉｎｇ ｅｒｒｏｒ ｉｓ ｏｂｔａｉｎｅｄ， ｗｈｉｃｈ ｉｓ １００ｍｓ． Ｐｇ ｔｒａｖｅｌ ｔｉｍｅ ｃａｎ ｂｅ ｔｒａｃｅｄ ｔｏ ｔｈｅ ｄｉｓｔａｎｃｅ ｏｆ
６０ｋｍ － ８５ｋｍ．
Ｔｈｅ ｖｅｒｔｉｃａｌ ａｘｅｓ ａｒｅ ｔｈｅ ｒｅｄｕｃｅｄ ｔｒａｖｅｌ ｔｉｍｅ ｗｉｔｈ ６ ｋｍ／ ｓ ａｎｄ ｂａｎｄｐａｓｓ ｆｉｌｔｅｒｅｄ ｂｅｔｗｅｅｎ ２ －８Ｈｚ．
２　ＲＥＧＵＬＡＲＩＺＥＤ ＴＲＡＶＥＬ⁃ＴＩＭＥ ＩＭＡＧＩＮＧ
Ｔｈｅ ｒｅｇｕｌａｒｉｚａｔｉｏｎ ｍｅｔｈｏｄ ｉｓ ｔｏ ｕｎｄｅｒｓｔａｎｄ ｔｈｅ ｕｎｄｅｒｄｅｔｅｒｍｉｎｅｄ ｐａｒｔ ｏｒ ｔｏ ｐｒｅｖｅｎｔ ｔｈｅ ｏｖｅｒ
ｆｉｔｔｉｎｇ ｏｆ ｄａｔａ ｂｙ ｔｈｅ ｉｎｔｒｏｄｕｃｔｉｏｎ ｏｆ ｓｏｍｅ ｃｏｎｓｔｒａｉｎｔｓ ｉｎ ｓｏｌｖｉｎｇ ｉｌｌ⁃ｃｏｎｄｉｔｉｏｎｅｄ ｅｑｕａｔｉｏｎｓ． Ｔｈｅｓｅ
４３２
Ｅａｒｔｈｑｕａｋｅ Ｒｅｓｅａｒｃｈ ｉｎ Ｃｈｉｎａ
Ｆｉｇ． １
３Ｄ ｓｕｒｖｅｙ ｇｅｏｍｅｔｒｙ ｏｆ ｔｈｅ Ａｎｈｕｉ ｓｅｉｓｍｉｃ ｅｘｐｅｒｉｍｅｎｔ
ｃｏｎｓｔｒａｉｎｔｓ ａｒｅ ｕｓｕａｌｌｙ ｐｕｎｉｔｉｖｅ ｍｅａｓｕｒｅｓ ｏｎ ｔｈｅ ｃｏｍｐｌｅｘｉｔｙ ｏｆ ｔｈｅ ｓｏｌｕｔｉｏｎ， ｓｕｃｈ ａｓ ｔｈｅ
ｓｍｏｏｔｈｎｅｓｓ ｃｏｎｔｒｏｌ ｏｆ ｔｈｅ ｉｍａｇｉｎｇ． Ｔｈｅｒｅｆｏｒｅ， ｔｈｅ ｒｅｇｕｌａｒｉｚａｔｉｏｎ ｐｒｏｃｅｓｓ， ｔｏ ａ ｃｅｒｔａｉｎ ｄｅｇｒｅｅ， ｃａｎ
ｂｅ ｒｅｇａｒｄｅｄ ａｓ ｔｈｅ ｉｎｖｅｒｓｉｏｎ ｏｆ ｔｈｅ ｔｒａｖｅｌ ｔｉｍｅ ｃｕｒｖｅ ｒａｔｈｅｒ ｔｈａｎ ｔｈｅ ｓｉｎｇｌｅ⁃ｐｏｉｎｔ ｔｒａｖｅｌ ｔｉｍｅ ｆｉｔｔｉｎｇ
（ Ｚｈａｎｇ ｅｔ ａｌ． ， １９９８） ． Ｔｈｕｓ， ｗｅ ｃａｎ ｃｏｎｓｔｒｕｃｔ ａｎ ｏｂｊｅｃｔｉｖｅ ｆｕｎｃｔｉｏｎ ｔｈａｔ ｉｎｃｌｕｄｅｓ ｔｈｅ ｓｍｏｏｔｈｎｅｓｓ
ｏｆ ｔｈｅ ｖｅｌｏｃｉｔｙ ｍｏｄｅｌ ａｎｄ ｔｈｅ ｄａｔａ ｆｉｔｔｉｎｇ ｄｅｇｒｅｅ （ Ｚｅｌｔ ｅｔ ａｌ． ， １９９８） ．
Φ ( ｍ ) ＝ δ ｔ Ｔ Ｃ ｄ－１ δｔ ＋ λ [ ｍ Ｔ Ｃ ｈ－１ ｍ ＋ ｓ ｚ ｍ Ｔ Ｃ ｖ－１ ｍ ]
（１）
Ｗｈｅｒｅ， ｍ ｉｓ ｍｏｄｅｌ ｖｅｃｔｏｒ； δｔ ｉｓ ｄａｔａ ｒｅｓｉｄｕａｌ； Ｃ ｄ ｉｓ ｄａｔａ ｃｏｖａｒｉａｎｃｅ ｍａｔｒｉｘ； Ｃ ｈ ａｎｄ Ｃ ｖ ａｒｅ
ｈｏｒｉｚｏｎｔａｌ ａｎｄ ｖｅｒｔｉｃａｌ ｓｍｏｏｔｈｎｅｓｓ ｍａｔｒｉｘ； λ ｉｓ ｃｏｏｒｄｉｎａｔｉｏｎ ｆａｃｔｏｒ ｏｆ ｔｈｅ ｄａｔａ ｆｉｔｔｉｎｇ ａｎｄ ｓｍｏｏｔｈｉｎｇ
ｄｅｇｒｅｅ； ｓ ｚ ｉｓ ｗｅｉｇｈｔｉｎｇ ｆａｃｔｏｒ ｏｆ ｖｅｒｔｉｃａｌ ｓｍｏｏｔｈｎｅｓｓ ｗｈｉｃｈ ｉｓ ｃｏｏｒｄｉｎａｔｉｏｎ ｗｉｔｈ ｔｈｅ ｗｅｉｇｈｔ ｏｆ
ｈｏｒｉｚｏｎｔａｌ ｓｍｏｏｔｈｎｅｓｓ． Ｓｏ， ｅｖｅｒｙ ｌｉｎｅａｒ ｉｔｅｒａｔｉｏｎ ｃａｎ ｂｏｉｌ ｄｏｗｎ ｔｏ ｏｎｅ ｓｏｌｕｔｉｏｎ ｏｆ δｍ， ｓｏ ｔｈａｔ ｔｈｅ
ｏｂｊｅｃｔｉｖｅ ｆｕｎｃｔｉｏｎ ｉｓ ｔｈｅ ｍｉｎｉｍｕｍ （ Ｚｅｌｔ ｅｔ ａｌ． ， １９９８） ．
－１ ／ ２
－１ ／ ２
éê Ｃ ｄ Ｌ ùú
éê Ｃ ｄ δｔ ùú
（２）
êê λ Ｃ ｈ úúδｍ ＝ êê － λ Ｃ ｈ ｍ０ úú
ë ｓｚ λ Ｃｖ û
ë － ｓ ｚ λ Ｃ ｖ ｍ０ û
Ｗｈｅｒｅ Ｌ ｉｓ ｔｈｅ ｐａｒｔｉａｌ ｄｅｒｉｖａｔｉｖｅ ｍａｔｒｉｘ ｏｆ ｔｈｅ ｏｂｊｅｃｔｉｖｅ ｆｕｎｃｔｉｏｎ； ｍ０ ｉｓ ｔｈｅ ｃｕｒｒｅｎｔ ｍｏｄｅｌ； δｍ ｉｓ
ｔｈｅ ｍｏｄｅｌ ｔｏ ｂｅ ｄｉｓｔｕｒｂｅｄ； ｔｈｅ ｎｅｗ ｍｏｄｅｌ ｖｅｃｔｏｒ ｉｓ ｍ ＝ ｍ０ ＋ δｍ． Ｆｏｒｍｕｌａ （２） ｃａｎ ｂｅ ｓｏｌｖｅｄ ｂｙ
ｔｈｅ ＬＳＱＲ ａｌｇｏｒｉｔｈｍ．
Ｖｏｌｕｍｅ ３０， Ｎｕｍｂｅｒ ３
４３３
Ｆｉｇ． ２
Ｅｘａｍｐｌｅｓ ｏｆ ｉｎ⁃ｌｉｎｅ ｄａｔａ ｆｏｒ ｓｈｏｔ Ｎｏ． １， Ｎｏ． ８， Ｎｏ． １６，
ａｎｄ Ｎｏ． ２０ ｉｎ ｔｈｅ ｐｒｏｆｉｌｅ ａｌｏｎｇ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ
Ｄｕｅ ｔｏ ｔｈｅ ｉｎｔｒｏｄｕｃｔｉｏｎ ｏｆ ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ａｎｄ ｖｅｒｔｉｃａｌ ｄｉｓｔｕｒｂａｎｃｅ ｆａｃｔｏｒｓ ｗｈｅｎ ａｐｐｌｙｉｎｇ ｔｈｅ
ｒｅｇｕｌａｒｉｚａｔｉｏｎ ａｌｇｏｒｉｔｈｍ， ｖｅｒｔｉｃａｌ ｒｅｓｏｌｕｔｉｏｎ ｉｓ ｇｒｅａｔｌｙ ｉｍｐｒｏｖｅｄ． Ｗｉｔｈ ｇｒｉｄｓ ｉｎ ｄｉｆｆｅｒｅｎｔ ｈｏｒｉｚｏｎｔａｌ
ａｎｄ ｖｅｒｔｉｃａｌ ｗｉｄｔｈｓ ａｎｄ ｔｈｅ ｒｅａｓｏｎａｂｌｅ ｓｅｌｅｃｔｉｏｎ ｏｆ ａ ｒｅｇｕｌａｒｉｚａｔｉｏｎ ｆａｃｔｏｒ， ｗｅ ｃｏｕｌｄ ｇｉｖｅ ａ
ｒｅａｓｏｎａｂｌｅ ｓｍｏｏｔｈ ｓｏｌｕｔｉｏｎ ｏｎ ｐａｒｔ ｏｆ ｔｈｅ ｎｏ⁃ｒａｙ ｃｏｖｅｒａｇｅ ａｒｅａ ａｌｏｎｇ ｔｈｅ ｐｒｏｆｉｌｅ， ａｎｄ ｔｈｅｎ
ｒｅａｓｏｎａｂｌｙ ｏｂｔａｉｎ ｔｈｅ ｓｍｏｏｔｈｅｓｔ ｓｏｌｕｔｉｏｎ ｏｆ ｔｈｅ ｍｉｎｉｍｕｍ ｇｒｉｄ ｐｏｉｎｔ ｉｎ ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ａｎｄ ｖｅｒｔｉｃａｌ
ｄｉｒｅｃｔｉｏｎｓ ｔｈｒｏｕｇｈｏｕｔ ｔｈｅ ｓｔｕｄｙ ａｒｅａ （ Ｚｅｌｔ ｅｔ ａｌ， １９９８）
２􀆰 １　Ｓｅｌｅｃｔｉｏｎ ｏｆ Ｉｎｉｔｉａｌ Ｍｏｄｅｌ
Ｆｏｒ ｔｈｅ ｓｅｉｓｍｉｃ ｆｉｒｓｔ ａｒｒｉｖａｌ ｔｒａｖｅｌ ｔｉｍｅ ｉｍａｇｉｎｇ， ａ ｒｅａｓｏｎａｂｌｅ ａｎｄ ｏｂｊｅｃｔｉｖｅ ｉｎｉｔｉａｌ ｍｏｄｅｌ
ｓｈｏｕｌｄ ｂｅ ｔｈｅ ｏｎｅ ｄｉｍｅｎｓｉｏｎａｌ ｍｏｄｅｌ ｗｈｅｒｅ ｔｈｅ ｌｏｎｇｉｔｕｄｉｎａｌ ｖｅｌｏｃｉｔｙ ｇｒａｄｉｅｎｔ ｄｉｓｔｒｉｂｕｔｉｏｎ ｉｓ
ｃｏｎｓｉｓｔｅｎｔ ｗｉｔｈ ｔｈｅ ｃｒｕｓｔａｌ ｖｅｌｏｃｉｔｙ ｇｒａｄｉｅｎｔ ｉｎ ｔｈｅ ｓｔｕｄｙ ａｒｅａ （ Ｚｅｌｔ ｅｔ ａｌ． ， １９９８） ． Ｉｎ ｔｈｉｓ ｐａｐｅｒ，
ｔｈｅ ｂｅｓｔ ｏｎｅ ｄｉｍｅｎｓｉｏｎａｌ ａｖｅｒａｇｅ ｍｏｄｅｌ ｉｓ ｓｅｔ ａｓ ｔｈｅ ｉｎｉｔｉａｌ ｍｏｄｅｌ ｏｂｔａｉｎｅｄ ｆｒｏｍ ｍｕｌｔｉｐｌｅ ｉｍａｇｉｎｇ
ｍｅｔｈｏｄ （ Ｚｅｌｔ ｅｔ ａｌ． ， ２００３） ， ａｎｄ ｔｈｅ ｉｎｉｔｉａｌ ｍｏｄｅｌ ｉｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ． ３（ ａ） ． Ｆｉｇ． ３（ ｂ） ｉｓ ｔｈｅ ｔｒａｖｅｌ⁃
ｔｉｍｅ ｃｏｎｔｒａｓｔ ｏｆ ｔｈｅ ｉｎｉｔｉａｌ ｍｏｄｅｌ ｉｎ ｃａｌｃｕｌａｔｉｏｎ ａｎｄ ｉｎ ｏｂｓｅｒｖａｔｉｏｎ． Ｉｎ Ｆｉｇ． ３（ ｂ） ， ｔｈｅ ｔｒａｖｅｌ⁃ｔｉｍｅ
ｃｕｒｖｅ ｃａｌｃｕｌａｔｅｄ ｆｒｏｍ ｔｈｅ ｉｎｉｔｉａｌ ｍｏｄｅｌ ｉｓ ｉｎ ｔｈｅ ｍｉｄ ｏｂｓｅｒｖｅｄ ｔｒａｖｅｌ ｔｉｍｅｓ， ｓｏ ｔｈｅ ｍｏｄｅｌ ｉｓ
ｒｅｌａｔｉｖｅｌｙ ｒｅａｓｏｎａｂｌｅ．
４３４
Ｅａｒｔｈｑｕａｋｅ Ｒｅｓｅａｒｃｈ ｉｎ Ｃｈｉｎａ
Ｆｉｇ． ３
（ ａ） １⁃Ｄ ｉｎｉｔｉａｌ ｍｏｄｅｌｓ ｕｓｅｄ ｉｎ ｔｈｅ ｍｉｎｉｍｕｍ⁃ｓｔｒｕｃｔｕｒｅ ａｐｐｒｏａｃｈ，
（ ｂ） Ｃｏｍｐｌｅｔｅ ｓｅｔ ｏｆ ｆｉｒｓｔ ａｒｒｉｖａｌ ｐｉｃｋｓ ｒｅｄｕｃｅｄ ａｔ ６􀆰 ０ｋｍ ／ ｓ ｖｅｒｓｕｓ
ｓｏｕｒｃｅ⁃ｒｅｃｅｉｖｅｒ ｏｆｆｓｅｔ ｗｉｔｈｏｕｔ ｒｅｇａｒｄ ｔｏ ｐｏｓｉｔｉｏｎ ａｌｏｎｇ ｔｈｅ ｐｒｏｆｉｌｅ
２􀆰 ２　Ｍｏｄｅｌ Ｐａｒａｍｅｔｅｒｉｚａｔｉｏｎ ａｎｄ Ｓｅｌｅｃｔｉｏｎ ｏｆ Ｆｏｒｗａｒｄ ａｎｄ Ｉｎｖｅｒｓｉｏｎ Ｐａｒａｍｅｔｅｒｓ
Ｔｈｅ ２ ｎｄ ｏｒｄｅｒ ｆｉｎｉｔｅ ｄｉｆｆｅｒｅｎｃｅ ａｌｇｏｒｉｔｈｍ ｉｓ ｕｓｅｄ ｔｏ ｓｏｌｖｅ ｔｈｅ ｗａｖｅ⁃ｆｒｏｎｔ ａｌｇｏｒｉｔｈｍ ｏｆ ｔｈｅ
ｅｉｋｏｎａｌ ｅｑｕａｔｉｏｎ （ Ｈｏｌｅ ｅｔ ａｌ． ， １９９５； Ｚｅｌｔ ｅｔ ａｌ． ， １９９８） ． Ｔｈｅ ０􀆰 ２５ｋｍ × ０􀆰 ２５ｋｍ ｇｒｉｄ ｉｓ ａｄｏｐｔｅｄ
ｉｎ ｔｈｅ ｆｏｒｗａｒｄ ｍｏｄｅｌ ｏｆ ｔｈｅ ｓｔｕｄｙ ａｒｅａ， ａｎｄ ｆｉｎａｌｌｙ ｔｈｅ ｎｕｍｂｅｒ ｏｆ ｇｒｉｄｓ ｆｏｒ ｔｈｅ ｆｏｒｗａｒｄ ｍｏｄｅｌ ｉｓ
１０８１ ｉｎ ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ｄｉｒｅｃｔｉｏｎ ａｎｄ ４１ｉｎ ｔｈｅ ｖｅｒｔｉｃａｌ ｄｉｒｅｃｔｉｏｎ， ｗｉｔｈ ａ ｔｏｔａｌ ｏｆ ４４，３２１ ｎｏｄｅｓ．
１􀆰 ００ｋｍ × ０􀆰 ２５ｋｍ ｇｒｉｄ ｉｓ ｕｓｅｄ ｉｎ ｔｈｅ ｉｎｖｅｒｓｉｏｎ ｍｏｄｅｌ， ｎａｍｅｌｙ， ｇｒｉｄ ｐｏｉｎｔ ｓｐａｃｉｎｇ ｉｓ １０００ｍ
ｉｎ ｈｏｒｉｚｏｎｔａｌ ｄｉｒｅｃｔｉｏｎ ａｎｄ ２５０ｍ ｉｎ ｔｈｅ ｖｅｒｔｉｃａｌ ｄｉｒｅｃｔｉｏｎ， ｆｏｒｍｉｎｇ ２７０ ｈｏｒｉｚｏｎｔａｌ ｇｒｉｄｓ ａｎｄ ４０
Ｖｏｌｕｍｅ ３０， Ｎｕｍｂｅｒ ３
４３５
ｖｅｒｔｉｃａｌ ｇｒｉｄｓ， ｗｉｔｈ ａ ｔｏｔａｌ ｏｆ １０８００ ｉｎｖｅｒｓｉｏｎ ｇｒｉｄｓ． Ｈｏｒｉｚｏｎｔａｌ ａｎｄ ｖｅｒｔｉｃａｌ ｓｍｏｏｔｈｉｎｇ ｗｅｉｇｈｔ
ｆａｃｔｏｒ ｓ ｚ ｔａｋｅｓ ｔｈｅ ｖａｌｕｅ ｏｆ ０􀆰 １５．
２􀆰 ３　Ｄａｔａ Ｆｉｔｔｉｎｇ Ｅｒｒｏｒ ａｎｄ Ｓｔａｂｉｌｉｔｙ Ａｎａｌｙｓｉｓ
Ｔｈｅ ｒｏｏｔ ｍｅａｎ ｓｑｕａｒｅ ｅｒｒｏｒ ｏｆ ｉｎｉｔｉａｌ ｍｏｄｅｌ ｒａｙ ｔｒａｃｉｎｇ ｉｓ １９８ｍｓ， ａｎｄ ｔｈｅ χ２ ｖａｌｕｅ ｉｓ ６􀆰 １３．
Ｔｈｅ ｒｅｄｕｃｔｉｏｎ ｒａｔｅ ｏｆ ｄａｍｐｉｎｇ ｆａｃｔｏｒ ｉｓ １００ ｆｏｒ ｅｖｅｒｙ １０ ｉｔｅｒａｔｉｏｎｓ， ａｎｄ ｔｈｅ ｆｉｎａｌ ｒｏｏｔ ｍｅａｎ ｓｑｕａｒｅ
ｅｒｒｏｒ ｉｓ ８０ｍｓ ａｎｄ χ２ ｖａｌｕｅ ｄｅｃｒｅａｓｅｓ ｔｏ １ ａｆｔｅｒ ５８ ｉｔｅｒａｔｉｏｎｓ． Ｔｈｅ ｔｒａｖｅｌ⁃ｔｉｍｅ ｒｅｓｉｄｕａｌ ｄｉｓｔｒｉｂｕｔｉｏｎｓ
ｏｆ ｉｎｉｔｉａｌ ａｎｄ ｆｉｎａｌ ｍｏｄｅｌｓ ａｒｅ ｓｈｏｗｎ ｉｎ Ｆｉｇ． ４， ｆｒｏｍ ｗｈｉｃｈ ｔｈｅ ｔｒａｖｅｌ ｔｉｍｅ ｆｉｔｔｉｎｇ ｅｒｒｏｒ ｒａｎｇｅ ｏｆ ｔｈｅ
ｉｎｉｔｉａｌ ｍｏｄｅｌ ｉｓ ± ０􀆰 ６ ａｎｄ ｔｈａｔ ｏｆ ｔｈｅ ｆｉｎａｌ ｍｏｄｅｌ ｃｏｎｓｔｒａｉｎｔｓ ｉｓ ｗｉｔｈｉｎ ± ０􀆰 １ｓ， ｉｎ ａｃｃｏｒｄａｎｃｅ ｗｉｔｈ
Ｆｉｇ． ４
Ｔｒａｖｅｌｔｉｍｅ ｒｅｓｉｄｕａｌｓ ａｓ ａ ｆｕｎｃｔｉｏｎ ｏｆ ｆｒｅｑｕｅｎｃｙ ｆｏｒ （ ａ） ｐｒｅｆｅｒｒｅｄ ｉｎｉｔｉａｌ
ｍｏｄｅｌ ａｎｄ （ ｂ） ｐｒｅｆｅｒｒｅｄ ｆｉｎａｌ ｍｏｄｅｌ ｆｒｏｍ ｒｅｇｕｌａｒｉｚｅｄ ｉｎｖｅｒｓｉｏｎ
Ｅａｒｔｈｑｕａｋｅ Ｒｅｓｅａｒｃｈ ｉｎ Ｃｈｉｎａ
４３６
ｓｅｉｓｍｉｃ ｐｈａｓｅ ｐｉｃｋｉｎｇ ｅｒｒｏｒ χ２ ｖａｌｕｅ ｏｆ １． Ｔｈｉｓ ｉｎｖｅｒｓｅ ｆｉｔｔｉｎｇ ｄａｔａ ｉｓ ｎｅｉｔｈｅｒ ｏｖｅｒｆｉｔｔｉｎｇ ｎｏｒ
ｌａｃｋｉｎｇ．
Ｔｈｅ ｆｉｎａｌ ｍｏｄｅｌ ｒａｙ ｃｏｖｅｒ ｉｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ． ５， ｆｒｏｍ ｗｈｉｃｈ ｗｅ ｃａｎ ｓｅｅ ａ ｔｏｔａｌ ｏｆ ６２８ ｒａｙｓ ｔｈａｔ
ｆｏｒｍ ａ ｂｅｔｔｅｒ ｃｏｖｅｒａｇｅ ｏｆ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ ｉｎ ｔｈｅ ｓｔｕｄｙ ａｒｅａ． Ｉｔ ｉｎｄｉｃａｔｅｓ ｔｈａｔ ｔｈｅ ｉｍａｇｉｎｇ ｒｅｓｕｌｔｓ ａｒｅ
ｒｏｂｕｓｔ ａｎｄ ｒｅｌｉａｂｌｅ．
３　ＩＭＡＧＩＮＧ ＲＥＳＵＬＴＳ ＡＮＤ ＤＩＳＣＵＳＳＩＯＮ
Ｆｉｇ． ６ ｓｈｏｗｓ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔａｌ ｖｅｌｏｃｉｔｙ ｓｔｒｕｃｔｕｒｅ ｍｏｄｅｌ ｏｆ Ｔｏｎｇｌｉｎｇ⁃Ａｎｑｉｎｇ⁃Ｍａ̓ａｎｓｈａｎ， ｗｈｅｒｅ
ｔｈｅ ｃｒｙｓｔａｌｌｉｎｅ ｕｐｐｅｒ ｃｒｕｓｔａｌ ｓｔｒｕｃｔｕｒｅ ｈａｓ ａｎ ｏｂｖｉｏｕｓ ｈｏｒｉｚｏｎｔａｌ ｄｉｆｆｅｒｅｎｃｅ， ｗｈｉｃｈ ｆｉｔｓ ｔｈｅ ｚｏｎｅｄ
ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ．
Ｉｎ ｔｈｅ Ａｎｑｉｎｇ ａｒｅａ， ｗｅｓｔ ｏｆ ｔｈｅ Ｌｕｊｉａｎｇ⁃Ｚｏｎｇｙａｎｇ ｂａｓｉｎ， ｔｈｅ ｄｅｐｔｈ ｏｆ ｂａｓｅｍｅｎｔ ｉｓ ｄｅｅｐ，
ａｂｏｕｔ ４ｋｍ － ６ｋｍ； ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ ｐｒｅｓｅｎｔｓ ａ ｔｙｐｉｃａｌ ａｌｔｅｒｎａｔｉｎｇ ｔｅｃｔｏｎｉｃ ｃｈａｒａｃｔｅｒｉｓｔｉｃ ｏｆ ｔｗｏ
ｕｐｌｉｆｔｓ ａｎｄ ｔｗｏ ｄｅｐｒｅｓｓｉｏｎｓ． Ｔｈｅ ｍｅａｓｕｒｉｎｇ ｐｏｉｎｔ ｉｓ ｌｏｃａｔｅｄ ｏｎ ｔｈｅ ｎｏｒｔｈｅｒｎ ｓｈｏｒｅ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ
Ｒｉｖｅｒ， ｗｉｔｈ ｔｈｉｃｋｅｒ ｓｅｄｉｍｅｎｔａｒｙ ｃｏｖｅｒ ａｎｄ ａ ｓｕｒｆａｃｅ ｖｅｌｏｃｉｔｙ ｏｆ ３􀆰 ５ｋｍ ／ ｓ ｗｈｅｒｅ ｔｈｅ ｓｕｒｆａｃｅ ｉｓ
ｍａｎｉｆｅｓｔｅｄ ａｓ ａ ｌａｒｇｅ⁃ａｒｅａ ａｌｌｕｖｉａｌ ｐｌａｉｎ．
Ｆｉｇ． ５
Ｒａｙ ｃｏｖｅｒａｇｅ ｔｈｒｏｕｇｈ ｔｈｅ ｆｉｎａｌ ｍｏｄｅｌ
Ｔｈｅ ｂｕｒｉｅｄ ｄｅｐｔｈ ｏｆ ｔｈｅ ｃｒｙｓｔａｌｌｉｎｅ ｂａｓｅｍｅｎｔ ｉｎ ｔｈｅ Ｌｕｊｉａｎｇ⁃Ｚｏｎｇｙａｎｇ ｂａｓｉｎ ｉｓ ａｂｏｕｔ ４ｋｍ －
５ｋｍ ａｎｄ ｉｔｓ ｓｅｄｉｍｅｎｔａｒｙ ｃｏｖｅｒ ｓｈｏｗｓ ａｎ ｏｂｖｉｏｕｓ ｍｏｒｐｈｏｌｏｇｙ ｏｆ ｄｅｐｒｅｓｓｉｏｎ ｂａｓｉｎ， ｗｈｉｃｈ ｉｎｄｉｃａｔｅｓ
ｔｈｅ ｅｘｔｅｎｓｉｏｎ⁃ｄｅｐｒｅｓｓｉｏｎ ｐｒｏｃｅｓｓ ｓｉｎｃｅ ｔｈｅ Ｍｅｓｏｚｏｉｃ． Ｔｈｅ ｓｕｒｆａｃｅ ｖｅｌｏｃｉｔｙ ｉｓ ａｂｏｕｔ ３􀆰 ５ － ４􀆰 ０ｋｍ ／ ｓ，
ｗｈｉｌｅ ｔｈｅｒｅ ａｒｅ ｏｂｖｉｏｕｓ ｌｏｗ⁃ｖｅｌｏｃｉｔｙ ａｎｏｍａｌｉｅｓ ａｔ ａｂｏｕｔ ３ｋｍ ｄｅｐｔｈ ｂｅｎｅａｔｈ ｔｈｅ ｅｄｇｅ ｏｆ ｔｈｉｓ ｂａｓｉｎ
ｗｈｉｃｈ ｍａｙ ｂｅ ｒｅｌａｔｅｄ ｔｏ ｖｏｌｃａｎｉｃ⁃ｍａｇｍａｔｉｃ ａｃｔｉｖｉｔｉｅｓ ｓｉｎｃｅ ｔｈｅ Ｍｅｓｏｚｏｉｃ， ｓｕｇｇｅｓｔｉｎｇ ｔｈａｔ ｔｈｅ
ｂａｓｅｍｅｎｔ ｄｅｐｒｅｓｓｉｏｎ ｉｓ ｃｈａｒａｃｔｅｒｉｚｅｄ ｂｙ ｖｏｌｃａｎｉｃ ｂａｓｉｎ ｆｅａｔｕｒｅｓ． Ｃｏｍｂｉｎｅｄ ｗｉｔｈ ｔｈｅ ｄｅｅｐ ｓｅｉｓｍｉｃ
ｒｅｆｌｅｃｔｉｏｎ ａｎｄ ＭＴ ｉｍａｇｉｎｇ ｒｅｓｕｌｔｓ ｃａｒｒｉｅｄ ｏｕｔ ｉｎ Ｇｕｉｃｈｉ ａｒｅａ， ｔｈｅ ｌｏｗ⁃ｓｐｅｅｄ ｍｅｄｉｕｍ （ ａｔ ｔｈｅ ｄｅｐｔｈ
ｏｆ ３ｋｍ） ｉｎ ｔｈｅ ｂａｓｉｎ ｍａｙ ｃｏｒｒｅｓｐｏｎｄ ｔｏ ｔｈｅ ｃｈａｎｎｅｌｓ ｏｆ ｍａｇｍａｔｉｃ ａｃｔｉｖｉｔｙ ｉｎ ｔｈｅ ｒｅｇｉｏｎ （ Ｄｏｎｇ
Ｓｈｕｗｅｎ ｅｔ ａｌ． ， ２０１０） ．
Ｔｈｅ ｍｅａｓｕｒｅｍｅｎｔ ｐｏｉｎｔｓ ｔｏ ｔｈｅ ｅａｓｔ ｏｆ Ｔｏｎｇｌｉｎｇ （ ａｃｒｏｓｓ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ） ａｒｅ ｌｏｃａｔｅｄ ｏｎ
ｔｈｅ ｓｏｕｔｈｅｒｎ ｂａｎｋ ｏｆ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ． Ｔｈｅ ｎｅａｒ⁃ｓｕｒｆａｃｅ ｖｅｌｏｃｉｔｙ ｖａｒｉｅｓ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｉｎ ｔｈｅ
ｈｏｒｉｚｏｎｔａｌ ｄｉｒｅｃｔｉｏｎ． Ａｆｔｅｒ ｃｒｏｓｓｉｎｇ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ， ｔｈｅ ｓｕｒｆａｃｅ ｖｅｌｏｃｉｔｙ ｓｏａｒｓ ｔｏ ａｂｏｕｔ ４􀆰 ５ｋｍ ／ ｓ
ａｎｄ ｔｈｅ ｓｅｄｉｍｅｎｔａｒｙ ｃｏｖｅｒ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｕｐｌｉｆｔｓ ｗｉｔｈ ｔｈｅ ｂａｓｅｍｅｎｔ ｄｅｐｔｈ ｂｅｉｎｇ ａｂｏｕｔ ２ － ４ｋｍ．
Ｆｒｏｍ ｔｈｅ ｅａｓｔ ｏｆ Ｔｏｎｇｌｉｎｇ ｔｏ Ｗｕｈｕ， ｔｈｅ ｕｐｌｉｆｔｉｎｇ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｔｈｅ ｓｅｄｉｍｅｎｔａｒｙ ｃｏｖｅｒ
ｃｏｒｒｅｓｐｏｎｄ ｔｏ ｔｈｅ ｃｏｍｐｌｅｘ ｒｅｆｌｅｃｔｉｏｎ ｉｎ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ ｓｈｏｗｅｄ ｂｙ ｔｈｅ ｄｅｅｐ ｓｅｉｓｍｉｃ ｒｅｆｌｅｃｔｉｏｎ ｉｎ
ｔｈｉｓ ｒｅｇｉｏｎ （ Ｌｕ Ｑｉｎｇｔｉａｎ ｅｔ ａｌ． ， ２００２ ） ， ｗｈｉｃｈ ｃｏｕｌｄ ｂｅ ｉｎｆｅｒｒｅｄ ａｓ ｆｏｌｄ， ｔｈｒｕｓｔ ａｎｄ ｉｎｔｒｕｓｉｖｅ
ｓｔｒｕｃｔｕｒｅｓ， ｉｎｄｉｃａｔｉｎｇ ｔｈａｔ ｓｔｒｏｎｇ ｃｏｍｐｒｅｓｓｉｖｅ ｄｅｆｏｒｍａｔｉｏｎ ｈａｓ ｏｃｃｕｒｒｅｄ ｉｎ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔ．
Ｖｏｌｕｍｅ ３０， Ｎｕｍｂｅｒ ３
４３７
Ｆｉｇ． ６
Ｕｐｐｅｒ ｃｒｕｓｔａｌ ｖｅｌｏｃｉｔｙ ｍｏｄｅｌ ｆｏｒ ｐｒｏｆｉｌｅ ａｌｏｎｇ ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ
ｏｂｔａｉｎｅｄ ｂｙ ｆｉｒｓｔ⁃ａｒｒｉｖａｌ ｔｒａｖｅｌｔｉｍｅ ｔｏｍｏｇｒａｐｈｙ
ＡＣＫＮＯＷＬＥＤＧＥＭＥＮＴＳ
Ｗｅ ａｒｅ ｇｒａｔｅｆｕｌ ｔｏ Ｐｒｏｆｅｓｓｏｒ Ｚｈａｎｇ Ｘｉａｎｋａｎｇ ｆｏｒ ｔｈｅ ｉｎｓｔｒｕｃｔｉｖｅ ｓｕｇｇｅｓｔｉｏｎｓ ｏｎ ｄｅａｌｉｎｇ ｗｉｔｈ
ｔｈｅ ｅｘｐｅｒｉｍｅｎｔａｌ ｄａｔａ ａｎｄ ｔｏ Ｐｒｏｆｅｓｓｏｒ Ｗａｎｇ Ｂａｏｓｈａｎ ａｎｄ Ｐｒｏｆｅｓｓｏｒ Ｙａｏ Ｈｕａｊｉａｎ ｆｏｒ ｔｈｅｉｒ ｓｕｐｐｏｒｔ
ａｎｄ ｈｅｌｐ ｉｎ ａｉｒｇｕｎ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ．
ＲＥＦＥＲＥＮＣＥＳ
Ｃｈａｎｇ Ｙｉｎｆｏ， Ｄｏｎｇ Ｓｈｕｗｅｎ， Ｈｕａｎｇ Ｄｅｚｈｉ． Ｏｎ ｔｅｃｔｏｎｉｃｓ ｏｆ “ Ｐｏｌｙ⁃Ｂａｓｅｍｅｎｔ ｗｉｔｈ Ｏｎｅ Ｃｏｖｅｒ” ｉｎ Ｍｉｄｄｌｅ⁃Ｌｏｗｅｒ
Ｙａｎｇｔｚｅ Ｃｒａｔｏｎ Ｃｈｉｎａ ［ Ｊ］ ． Ｖｏｌｃａｎｏｌｏｇｙ ＆ Ｍｉｎｅｒａｌ Ｒｅｓｏｕｒｃｅｓ， １９９６， １７（ Ｉ） ：１ － １５（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ
ａｂｓｔｒａｃｔ） ．
Ｃｈａｎｇ Ｙｉｎｆｏ， Ｌｉｕ Ｘｉａｎｇｐｅｉ， Ｗｕ Ｙａｎｃｈａｎｇ． Ｃｏｐｐｅｒ ａｎｄ Ｉｒｏｎ Ｍｅｔａｌｌｏｇｅｎｉｃ Ｂｅｌｔ ｉｎ ｔｈｅ Ｍｉｄｄｌｅ ａｎｄ Ｌｏｗｅｒ Ｒｅａｃｈｅｓ ｏｆ
ｔｈｅ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ［ Ｍ］ Ｂｅｉｊｉｎｇ： Ｇｅｏｌｏｇｉｃａｌ Ｐｕｂｌｉｓｈｉｎｇ Ｈｏｕｓｅ， １９９７􀆰 １ － ２３９（ ｉｎ Ｃｈｉｎｅｓｅ）
Ｄｏｎｇ Ｓｈｕｗｅｎ， Ｘｉａｎｇ Ｈｕａｉｓｈｕｎ， Ｇａｏ Ｒｕｉ ｅｔ ａｌ． Ｄｅｅｐ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｏｒｅ ｆｏｒｍａｔｉｏｎ ｗｉｔｈｉｎ Ｌｕｊｉａｎｇ⁃Ｚｏｎｇｙａｎｇ
ｖｏｌｃａｎｉｃ ｏｒｅ ｃｏｎｃｅｎｔｒａｔｅｄ ａｒｅａ ｉｎ Ｍｉｄｄｌｅ ｔｏ Ｌｏｗｅｒ Ｒｅａｃｈｅｓ ｏｆ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ［ Ｊ］ ． Ａｃｔａ Ｐｅｔｒｏｌｏｇｉｃａ Ｓｉｎｉｃａ，
２０１０， ２６（９） ：２５２９ － ２５４２（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｄｏｎｇ Ｓｈｕｗｅｎ， Ｘｉａｎｇ Ｈｕａｉｓｈｕｎ， Ｇａｏ Ｒｕｉ ｅｔ ａｌ． Ｄｅｅｐ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｏｒｅ ｆｏｒｍａｔｉｏｎ ｗｉｔｈｉｎ Ｌｕｊｉａｎｇ⁃Ｚｏｎｇｙａｎｇ
ｖｏｌｃａｎｉｃ ｏｒｅ ｃｏｎｃｅｎｔｒａｔｅｄ ａｒｅａ ｉｎ ｍｉｄｄｌｅ ｔｏ ｌｏｗｅｒ ｒｅａｃｈｅｓ ｏｆ Ｙａｎｇｔｚｅ Ｒｉｖｅｒ ［ Ｊ］ ． Ａｃｔａ Ｐｅｔｒｏｌｏｇｉｃａ Ｓｉｎｉｃａ，
２０１０， ２６（９） ：２５２９ － ２５４２（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｄｕａｎ Ｙｏｎｇｈｏｎｇ， Ｚｈａｎｇ Ｘｉａｎｋａｎｇ， Ｆａｎｇ Ｓｈｅｎｇｍｉｎｇ． Ｔｈｒｅｅ⁃ｄｉｍｅｎｓｉｏｎａｌ ｆｉｎｉｔｅ⁃ｄｉｆｆｅｒｅｎｃｅ ｔｏｍｏｇｒａｐｈｙ ｏｆ ｖｅｌｏｃｉｔｙ
ｓｔｒｕｃｔｕｒｅ ｏｆ ｔｈｅ ｕｐｐｅｒ ｃｒｕｓｔａｌ ｉｎ Ｎｏｒｔｈ Ｃｈｉｎａ ［ Ｊ］ Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ， ２００２， ４５（３） ：３６２ － ３６９（ ｉｎ
Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｈｏｌｅ Ｊ． Ａ． ， Ｚｅｌｔ Ｂ． Ｃ． ３⁃Ｄ ｆｉｎｉｔｅ⁃ｄｉｆｆｅｒｅｎｃｅ ｒｅｆｌｅｃｔｉｏｎ ｔｒａｖｅｌ ｔｉｍｅｓ ［ Ｊ］ ． Ｇｅｏｐｈｙｓ Ｊ Ｉｎｔ， １９９５， １２１： ４２７ － ４３４．
Ｊｉａ Ｙｕｐｅｎｇ， Ｗａｎｇ Ｆｕｙｕｎ， Ｔｉａｎ Ｘｉａｏｆｅｎｇ ｅｔ ａｌ． Ｔｏｍｏｇｒａｐｈｙ ｏｆ ｈｉｇｈ ｒｅｓｏｌｕｔｉｏｎ ｓｅｉｓｍｉｃ ｒｅｆｒａｃｔｉｏｎ ｔｒａｖｅｌｔｉｍｅｓ ｉｎ
Ｔｉａｎｊｉｎ⁃Ｂｅｉｊｉｎｇ Ｐｒｏｆｉｌｅ ［ Ｊ］ ． Ｎｏｒｔｈｗｅｓｔｅｒｎ Ｓｅｉｓｍｏｌｏｇｉｃａｌ Ｊｏｕｒｎａｌ， ２０１２， ３４（４） ：３７５ － ３８２ ３４（４） ：３７５ － ３８２
（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
４３８
Ｅａｒｔｈｑｕａｋｅ Ｒｅｓｅａｒｃｈ ｉｎ Ｃｈｉｎａ
Ｌｖ Ｑｉｎｇｔｉａｎ， Ｄｏｎｇ Ｓｈｕｗｅｎ， Ｓｈｉ Ｄａｎｉａｎ ｅｔ ａｌ． Ｌｉｔｈｏｓｐｈｅｒｅ ａｒｃｈｉｔｅｃｔｕｒｅ ａｎｄ ｇｅｏｄｙｎａｍｉｃ ｍｏｄｅｌ ｏｆ ｍｉｄｄｌｅ ａｎｄ ｌｏｗｅｒ
ｒｅａｃｈｅｓ ｏｆ Ｙａｎｇｔｚｅ Ｍｅｔａｌｌｏｇｅｎｉｃ Ｂｅｌｔ： Ａ ｒｅｖｉｅｗ ｆｒｏｍ ＳｉｎｏＰｒｏｂｅ ［ Ｊ］ ． Ａｃｔａ Ｐｅｔｒｏｌｏｇｉｃａ Ｓｉｎｉｃａ， ２０１４， ３０（４） ：
８８９ － ９０６（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｌｖ Ｑｉｎｇｔｉａｎ， Ｈｕａｎｇ Ｄｏｎｇｄｉｎｇ， Ｈｏｕ Ｚｅｎｇｑｉａｎ ｅｔ ａｌ． Ｄｅｅｐ ｓｅｉｓｍｉｃ ｒｅｆｌｅｃｔｉｏｎ ｉｍａｇｅ ｏｆ ｃｒｕｓｔａｌ ｓｔｒｕｃｔｕｒｅ ｉｎ Ｔｏｎｇｌｉｎｇ
Ｏｒｅ Ｄｉｓｔｒｉｃｔ ［ Ｊ］ ． Ｍｉｎｅｒａｌ Ｄｅｐｏｓｉｔｓ， ２００２， （ Ｓｕｐｐｌ． ） ：１１７３ － １１７６（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｐａｎ Ｊｉｓｈｕｎ， Ｚｈａｎｇ Ｘｉａｎｋａｎｇ， Ｚｈａｏ Ｐｉｎｇ ｅｔ ａｌ． ２Ｄ ｍｕｌｔｉｓｃａｌｅ ｎｏｎ⁃ｌｉｎｅａｒ ｓｅｉｓｍｉｃ ｖｅｌｏｃｉｔｙ ｉｍａｇｉｎｇ ［ Ｊ］ ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ， ２００８， ５１（１） ：１９７ － ２０５（ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ ａｂｓｔｒａｃｔ） ．
Ｗａｎｇ Ｂ． ， Ｇｅ Ｈ． ， Ｗｅｉ Ｙ． ， ｅｔ ａｌ． Ｔｒａｎｓｍｉｔｔｉｎｇ ｓｅｉｓｍｉｃ ｓｔａｔｉｏｎ ｍｏｎｉｔｏｒｓ ｆａｕｌｔ ｚｏｎｅ ａｔ ｄｅｐｔｈ ［ Ｊ］ ． Ｅｏｓ Ｔｒａｎｓａｃｔｉｏｎｓ
Ａｍｅｒｉｃａｎ Ｇｅｏｐｈｙｓｉｃａｌ Ｕｎｉｏｎ， ２０１２， ９３（５） ：４９ － ５０．
Ｘｕ Ｃｈａｏｆａｎ， Ｚｈａｎｇ Ｘｉａｎｋａｎｇ， Ｚｈａｎｇ Ｊｉａｎｓｈｉ ｅｔ ａｌ． Ｒａｙ ｈｉｔ ａｎａｌｙｓｉｓ ｍｅｔｈｏｄ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎ ｔｏ ｃｏｍｐｌｅｘ ｕｐｐｅｒ
ｃｒｕｓｔａｌ ｓｔｒｕｃｔｕｒｅ ｓｕｒｖｅｙ ［ Ｊ］ ． Ａｃｔａ Ｓｅｉｓｍｏｌｏｇｉｃａ Ｓｉｎｉｃａ， ２００６， １９（２） ：１７３ － １８２ （ ｉｎ Ｃｈｉｎｅｓｅ ｗｉｔｈ Ｅｎｇｌｉｓｈ
ａｂｓｔｒａｃｔ） ．
Ｚｅｌｔ Ｃ． Ａ． ， Ｂａｒｔｏｎ Ｐ． Ｊ． ３Ｄ ｓｅｉｓｍｉｃ ｒｅｆｒａｃｔｉｏｎ ｔｏｍｏｇｒａｐｈｙ： Ａ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｔｗｏ ｍｅｔｈｏｄｓ ａｐｐｌｉｅｄ ｔｏ ｄａｔａ ｆｒｏｍ ｔｈｅ
Ｆａｅｒｏｅ Ｂａｓｉｎ ［ Ｊ］ ． Ｊ Ｇｅｏｐｈｙｓ Ｒｅｓ， １９９８， １０３： ７１８７ － ７２１０．
Ｚｅｌｔ Ｃ． Ａ． ， Ｓａｉｎ Ｋ． ， Ｎａｕｍｅｎｋｏ Ｊ． Ｖ． ， ｅｔ ａｌ． Ａｓｓｅｓｓｍｅｎｔ ｏｆ ｃｒｕｓｔａｌ ｖｅｌｏｃｉｔｙ ｍｏｄｅｌｓ ｕｓｉｎｇ ｓｅｉｓｍｉｃ ｒｅｆｒａｃｔｉｏｎ ａｎｄ
ｒｅｆｌｅｃｔｉｏｎ ｔｏｍｏｇｒａｐｈｙ ［ Ｊ］ ． Ｇｅｏｐｈｙｓ Ｊ Ｉｎｔ， ２００３， １５３：６０９ － ６２６．
Ｚｈａｎｇ Ｊ． ， Ｔｏｋｓｏｚ Ｍ． Ｎ． Ｎｏｎｌｉｎｅａｒ ｒｅｆｒａｃｔｉｏｎ ｔｒａｖｅｌｔｉｍｅ ｔｏｍｏｇｒａｐｈｙ ［ Ｊ］ ． Ｇｅｏｐｈｙｓｉｃｓ， １９９８， ６３ （ ５ ） ： １７２６ －
１７３７．
Ａｂｏｕｔ ｔｈｅ Ａｕｔｈｏｒ
Ｔｉａｎ Ｘｉａｏｆｅｎｇ， ｂｏｒｎ ｉｎ １９７９， ｉｓ ａ ｓｅｎｉｏｒ ｅｎｇｉｎｅｅｒ ａｔ ｔｈｅ Ｇｅｏｐｈｙｓｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ Ｃｅｎｔｅｒ，
Ｃｈｉｎａ Ｅａｒｔｈｑｕａｋｅ Ａｄｍｉｎｉｓｔｒａｔｏｒ． Ｈｉｓ ｍａｊｏｒ ｉｓ ｔｈｅ ｓｔｕｄｙ ｏｆ ｓｅｉｓｍｏｌｏｇｙ， ｔｒａｖｅｌ ｔｉｍｅ ｉｍａｇｉｎｇ， ａｎｄ
ｃｒｕｓｔａｌ ｌｉｔｈｏｓｐｈｅｒｅ ｅｖｏｌｕｔｉｏｎ． Ｅ⁃ｍａｉｌ：ｔｉａｎｄｌｅｒ＠ ｈｏｔｍａｉｌ． ｃｏｍ