廣西喀斯特地貌會(huì)放大地震晃動(dòng)？且衰減更快？

廣西昨天凌晨的5.2級(jí)地震，造成了數(shù)棟房屋倒塌。

難以想象這是5.2級(jí)地震造成的破壞。

除了這次地震震源深度淺，這些房屋可能年久失修，且設(shè)防標(biāo)準(zhǔn)低外。

當(dāng)?shù)靥厥獾目λ固氐孛?，?huì)不會(huì)是幫兇？

從廣西全境的烈度分布看，所有臺(tái)站的儀器烈度非常?。ㄐ∮?度），不至于造成如此大的破壞吧？

由于測(cè)量地震烈度的儀器距離震中不可能太近，而震中房屋又出現(xiàn)了倒塌，這只能說(shuō)明，該地震的強(qiáng)地震動(dòng)隨著遠(yuǎn)離震中，迅速衰減。

因此，抽象出來(lái)的物理問(wèn)題是這樣的：

在近地表溶洞廣布的喀斯特地貌地區(qū)，由于大量空洞的存在，相比于沒(méi)有空洞的情況，是否會(huì)放大震中的地面晃動(dòng)，并且這一晃動(dòng)僅集中在震中附近，遠(yuǎn)離震中則迅速衰減？

于是，筆者翻出了學(xué)生時(shí)代寫的二維彈性波場(chǎng)模擬代碼，改了改，完成了本次數(shù)值模擬。

先說(shuō)結(jié)論：確實(shí)是這樣的，由于空洞對(duì)波場(chǎng)的散射，會(huì)導(dǎo)致震中附近的晃動(dòng)更大且持續(xù)時(shí)間更長(zhǎng)，并且這一晃動(dòng)，隨著遠(yuǎn)離震中（可能只需要幾百米），就迅速衰減。

我們先看沒(méi)有空洞的波場(chǎng)（Ux）：

整個(gè)波場(chǎng)（U）：

現(xiàn)在看看有空洞的情況(Ux) ，近地表白色圓形就是空洞 ：

整個(gè)波場(chǎng)（U）：

由上圖所示，由于空洞的存在，近地表波場(chǎng)散色非常嚴(yán)重，導(dǎo)致波場(chǎng)持續(xù)時(shí)間很長(zhǎng)。

現(xiàn)在，我們看看地表臺(tái)站記錄到的波形：

先看震中的：

明顯可見(jiàn)，當(dāng)?shù)叵掠卸磿r(shí)（紅線），地表位移的振幅大了10倍，持續(xù)時(shí)間長(zhǎng)了數(shù)倍，這足以說(shuō)明 喀斯特地貌是這次地震的重要幫兇。

我們?cè)倏纯?，距離震中250米的臺(tái)站記錄（持續(xù)時(shí)間增長(zhǎng)，振幅略大）：

再看看距震中500米的臺(tái)站（持續(xù)時(shí)間增長(zhǎng)，振幅反而更?。?/p>

再看看距震中3500米的臺(tái)站（持續(xù)時(shí)間增長(zhǎng)，振幅更小了）：

這些結(jié)果，足以說(shuō)明：

喀斯特地貌區(qū)域由于地下大量空洞的存在，會(huì)造成波場(chǎng)散射，導(dǎo)致震中附近非常近（可能就數(shù)百米）區(qū)域內(nèi)地震波的振幅會(huì)有數(shù)量級(jí)的增加，并且隨著遠(yuǎn)離震中，迅速衰減（即散射衰減）。

這和我們觀測(cè)到的現(xiàn)象完全一致。

喀斯特地貌是這次地震震中附近房屋倒塌的“幫兇”！

下面講一講模擬參數(shù)設(shè)置：

模型長(zhǎng)8千米，深度5 千米 , 三角形網(wǎng)格，網(wǎng)格大小50米。

P波速度：3500 米/秒 ，S波速度：2020米/秒，密度2000 kg/m^3。

模型除地表外為PML吸收邊界。

震源為雷克子波，主頻5 Hz，位于模型中間，即x=4 千米 ，深度2.5千米處。

計(jì)算方法為有限元，進(jìn)行的二維彈性波場(chǎng)模擬。

對(duì)于空洞模型，除了存在空洞外，其余設(shè)置和沒(méi)有空洞模型完全一致。

空洞設(shè)置：在近地表100-500米深度范圍內(nèi)隨時(shí)生成了23個(gè)圓形空洞，其直徑為100到400米不等。

模型時(shí)間步長(zhǎng)：0.00286秒，模擬時(shí)長(zhǎng)為9秒。

空洞模型網(wǎng)格剖分：

空洞位置及大?。?/p>

Void 1: x=6668 m, y=393 m, d=250 m

Void 2: x=6168 m, y=340 m, d=123 m

Void 3: x=1162 m, y=488 m, d=383 m

Void 4: x=2841 m, y=310 m, d=173 m

Void 5: x=5688 m, y=180 m, d=282 m

Void 6: x=1618 m, y=374 m, d=222 m

Void 7: x=6423 m, y=204 m, d=114 m

Void 8: x=2134 m, y=488 m, d=319 m

Void 9: x=5262 m, y=404 m, d=148 m

Void 10: x=4051 m, y=463 m, d=226 m

Void 11: x=3469 m, y=402 m, d=200 m

Void 12: x=2389 m, y=271 m, d=144 m

Void 13: x=4538 m, y=371 m, d=116 m

Void 14: x=4072 m, y=191 m, d=107 m

Void 15: x=3217 m, y=206 m, d=137 m

Void 16: x=4861 m, y=171 m, d=110 m

Void 17: x=3796 m, y=317 m, d=114 m

Void 18: x=1312 m, y=116 m, d=146 m

Void 19: x=5083 m, y=128 m, d=138 m

Void 20: x=1847 m, y=106 m, d=185 m

Void 21: x=2178 m, y=116 m, d=169 m

Void 22: x=5983 m, y=494 m, d=156 m

Void 23: x=6191 m, y=109 m, d=110 m

無(wú)空洞模型網(wǎng)格剖分：