本發(fā)明提供了一種在復(fù)雜各向異性介質(zhì)中模擬地震波的方法及裝置，方法包括：給數(shù)值化的三維地質(zhì)模型的每個(gè)網(wǎng)格點(diǎn)的彈性常數(shù)和密度賦值；獲取三維地質(zhì)模型中網(wǎng)格點(diǎn)沿三個(gè)坐標(biāo)軸的速度分量的一階空間導(dǎo)數(shù)；將速度分量的一階空間導(dǎo)數(shù)代入到彈性介質(zhì)中彈性波方程組中的第一方程式，并進(jìn)行求解獲取應(yīng)力分量；獲取三維地質(zhì)模型中網(wǎng)格點(diǎn)沿三個(gè)坐標(biāo)軸的應(yīng)力分量的一階空間導(dǎo)數(shù)；將應(yīng)力分量的一階空間導(dǎo)數(shù)代入到彈性介質(zhì)中彈性波方程組中的第二方程式，并進(jìn)行求解獲取速度分量；根據(jù)通過第二方程式獲取的速度分量得到當(dāng)前時(shí)間第一時(shí)刻的速度分量的值，作為當(dāng)前時(shí)刻的地震波的記錄；獲取下一時(shí)刻的地震波的記錄，直至獲取時(shí)間長度t的地震波的波場值。