[发明专利]一种微地震双力偶震源的加载方法

权利要求书

1.一种微地震双力偶震源的加载方法，其特征在于：所述方法基于波动方程有限差分数值模拟实现微地震双力偶源的加载；

所述方法包括：

(1)确定断层面参数：走向Φ、倾角δ和滑动角λ；

(2)根据所述断层面参数建立地震矩张量；

(3)由所述地震矩张量的各个元素值获得双力偶源的加载方式；

(4)通过集中力组合加载方式完成双力偶源的加载；

(5)完成双力偶源波场模拟。

2.根据权利要求1所述的微地震双力偶震源的加载方法，其特征在于：所述步骤(1)中，走向Φ是由正北至断层线顺时针量取的角度，范围为0°至360°；倾角δ是断层面与地平面间的夹角，范围为0°至90°；滑动角λ是以断层走向为基础，向上为正，向下为负，范围为-180°至+180°。

3.根据权利要求2所述的微地震双力偶震源的加载方法，其特征在于：所述步骤(2)是这样实现的：

地震矩张量为3×3形式的矩阵，共有9个元素，表示如下：

M=M11M12M13M21M22M23M31M32M33]]>

其中，M_ij表示一个力偶，M₁₁表示作用于1轴、平行于1轴的力偶；M₁₂表示作用于2轴、平行于1轴的力偶；M₁₃表示作用于3轴、平行于1轴的力偶；其他元素以此类推；

所述地震矩张量与断层面参数的表达式如下：

M₁₁＝-M₀(sinδcosλsin2φ+sin2δsinλsin²φ)

M₂₂＝M₀(sinδcosλsin2φ-sin2δsinλcos²φ)

M₃₃＝M₀sin2δsinλ＝-(M₁₁+M₂₂)

M12=M0(sinδcosλcos2φ+12sin2δsinλsin2φ)]]>

M₁₃＝-M₀(cosδcosλcosφ+cos2δsinλsinφ)

M₂₃＝-M₀(cosδcosλsinφ-cos2δsinλcosφ)

其中，M₀表示地震力矩。