地下爆炸地震波的数值模拟

合集下载

地震波波动方程数值模拟方法(严选优质)

地震波波动方程数值模拟方法地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。

克希霍夫积分法引入射线追踪过程，本质上是波动方程积分解的一个数值计算，在某种程度上相当于绕射叠加。

该方法计算速度较快，但由于射线追踪中存在着诸如焦散、多重路径等问题，故其一般只能适合于较简单的模型，难以模拟复杂地层的波场信息。

傅里叶变换法是利用空间的全部信息对波场函数进行三角函数插值，能更加精确地模拟地震波的传播规律，同时，利用快速傅里叶变换(FFT)进行计算，还可以提高运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型的适用性差，尤其是不适应于速度横向变化剧烈的模型.波动方程有限元法的做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最后求解总体矩阵得到波动方程的数值解；其主要优点是理论上可适宜于任意地质体形态的模型，保证复杂地层形态模拟的逼真性，达到很高的计算精度，但有限元法的主要问题是占用内存和运算量均较大，不适用于大规模模拟，因此该方法在地震波勘探中尚未得到广泛地应用。

相对于上述几种方法，有限差分法是一种更为快速有效的方法。

虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。

声波方程的有限差分法数值模拟对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：)()(2222222t S zu x u v t u +∂∂+∂∂=∂∂ （4-1） (,)v x z 是介质在点（x , z ）处的纵波速度，u 为描述速度位或者压力的波场，)(t s 为震源函数。

为求式（4-1）的数值解，必须将此式离散化，即用有限差分来逼近导数，用差商代替微商。

为此，先把空间模型网格化（如图4-1所示）。

设x 、z 方向的网格间隔长度为h ∆，t ∆为时间采样步长，则有：z∆,i j1,i j +2,i j+1,i j-h i x ∆= （i 为正整数）h j z ∆= (j 为正整数)t n t =∆ (n 为正整数)k j i u , 表示在(i,j)点，k 时刻的波场值。

数值分析在地震模拟中的应用研究

数值分析在地震模拟中的应用研究地震是一种常见的自然灾害，其给人们的生命安全和财产造成了严重的威胁。

为了更好地了解地震的性质和影响，科研人员一直在致力于地震模拟的研究。

数值分析作为一种有效的工具，已经在地震模拟中得到了广泛应用。

本文将探讨数值分析在地震模拟中的应用研究，并评估其在地震研究领域的影响。

一、地震模拟的背景和意义地震模拟是通过计算机模拟地震过程，以更好地了解和预测地震灾害。

传统的地震观测往往只能提供有限的信息，而地震模拟可以通过数值计算，提供更全面、详细的地震数据。

这对于地震研究和防灾减灾工作具有重要意义。

二、数值分析在地震模拟中的应用1. 地震波传播仿真地震波传播是地震模拟中的核心问题之一。

数值分析可以用来模拟地震波在地下介质中的传播过程。

通过建立复杂的波动方程和边界条件，数值模拟可以计算出地震波在各个点的传播路径和强度，从而得到地震波的传播规律。

2. 地震损伤预测地震模拟可以模拟地震对建筑物和结构的影响，帮助科研人员预测建筑物的损伤情况。

通过数值分析，可以模拟地震时建筑物的应力分布和变形情况，评估结构的强度和稳定性，为抗震设计提供依据。

3. 地震风险评估地震模拟可以通过模拟多个地震事件，分析某个地区地震灾害的可能性和潜在影响。

通过数值分析，可以得到地震活动的统计数据，包括频率、震级等，从而评估地震风险，制定相应的预警和防灾措施。

三、数值分析在地震研究领域的影响数值分析在地震研究领域的应用已经取得了显著的成果，并在某些方面取代了传统的试验方法。

通过数值模拟，研究人员可以得到更多的实验数据，进一步了解地震过程中的细节和规律。

同时，数值模拟可以模拟不同地震条件下的情景，提供更全面的研究结果。

然而，数值分析在地震研究中仍面临一些挑战。

首先，地震模拟需要大量的计算资源，包括高性能计算机和大规模存储设备。

其次，地震模拟需要建立复杂的数学模型和边界条件，对研究人员的能力和经验有一定要求。

最后，地震模拟结果的准确性和可靠性需要进一步验证和评估。

减震沟减震效果的数值模拟研究——打印

图 8 2种工况下水平峰值震速比较
第 23卷第 3期梁开水等减震沟减震效果的数值模拟研究
20
爆破 2006年 9月
的衰减规律 ,其中横坐标是地表质点的爆心距。 2. 3 不同抽槽长度的减震效果比较
在进行有减震沟情况下计算时 ,减震沟深 3 m , 宽 1. 2 m ,炮孔深 1. 5 m ,其中堵塞段为 0. 5 m ,最小抵抗线为 2 m;保护区长 40 m ,开挖侧长 35 m ,模型高为 10 m;抽槽长度分别是 10 m、15 m。典型计算模型见图 3,图 4、图 5 是保护区侧地表质点峰值震速随爆心距的衰减规律。
为了比较不同深度减震沟的减震效果 ,考虑了减震沟深 3 m 和深 6 m ,宽度为 1. 2 m 的工况 ,比较
图 7 2种工况下垂直峰值震速比较
由图 6、图 7 可以看出 ,减震沟越深 ,其隔震效果越好 ,对水平峰值的影响较大 ,对垂直峰值震速的影响相对较小。这是由于减震沟越深 ,爆破地震波需要绕射进入保护侧的路径越长 ,其损失的能量也越多。 2. 5 不同宽度减震沟减震效果比较
第 23卷第 3期 2006年 9月
爆破 BLAST ING
Vol. 23 No. 3 Sep. 2006
文章编号 : 1001 - 487X (2006) 03 - 0018 - 04
减震沟减震效果的数值模拟研究
梁开水 ,陈天珠 ,易长平
(武汉理工大学资源与环境工程学院 ,湖北武汉 430070)
Technology,W uhan 430070, China)
A bstract: V ibration2isolating effect of vibration2isolating slot is simulated in different conditions by an exp losive

大尺度地震波的数值模拟与预测

大尺度地震波的数值模拟与预测地震是地球内部能量释放的一种自然现象，会造成巨大的破坏和人员伤亡。

了解地震波的传播规律以及对其进行准确的模拟与预测对于地震灾害的预防和减轻具有重要意义。

本文将探讨大尺度地震波的数值模拟与预测方法，并分析其在地震灾害管理中的应用。

地震波是地震能量释放后在岩石、土壤和水体中传播所产生的波动，其传播速度和路径受到地质构造和介质性质的影响。

数值模拟地震波传播可以通过在计算机中解析地震方程来实现。

地震方程是描述地震波传播的基本方程，通常采用波动方程形式。

通过在三维空间中离散地震方程，可以得到地震波在不同地点的振幅和传播速度信息，从而实现对地震波传播的模拟。

为了进行大尺度地震波的数值模拟，需要获取大规模的地质结构模型和地震波速度模型。

地质结构模型可以通过地质勘探和地震资料分析得到，用于刻画地下介质的层状结构和性质。

地震波速度模型则是描述地震波在不同介质中传播速度的参数，可以通过地震资料和地震勘探技术获取。

利用这些模型，可以在计算机中建立相应的数值模型，在模拟地震波传播之前对其进行预测。

在进行地震波数值模拟之前，需要进行验证和校准。

验证是指将数值模拟结果与实测资料进行对比，以验证模拟的准确性和可靠性。

校准则是通过调整模拟参数，优化地震波模拟结果，使其与实测资料吻合程度更高。

验证和校准过程的完成可以提高地震波数值模拟的可信度，并为后续的预测工作打下基础。

大尺度地震波的数值模拟可以用于地震灾害管理的多个方面。

首先，通过模拟地震波在不同介质中的传播，可以预测地震造成的破坏范围和程度。

这对于城市规划、建筑设计以及灾害应急管理具有重要意义，可以提前采取相应的措施减轻地震灾害的影响。

其次，模拟地震波传播还可以用于评估地震烈度和地表运动速度，为地震灾害风险评估提供依据。

最后，地震波数值模拟还可以帮助科学家深入研究地震过程和地震发生机理，为地震灾害的原因和规律提供科学依据。

然而，大尺度地震波的数值模拟和预测也面临一些挑战。

有限差分法地震波传播数值模拟

2k
2 x
ω2
⎟⎟⎠⎞2
⎪⎫ ⎬ ⎪⎭
∂3P ∂z∂t 2
−
v2 4
∂3P ∂x2∂z
+
3v 4
∂3P ∂x2∂t
−
1 v
∂3P ∂t 3
=
0
Elastic Wave: (Bottom Boundary)
Utt = α 2U xx + β 2U zz + (α 2 − β 2 )Wxz Wtt = β 2Wxx + α 2Wzz + (α 2 − β 2 )U xz
2
MM 4
⎤ ⎥ ⎥
⎢⎢⎡CC12((MM
) )
⎤ ⎥ ⎥
M = M
6
⎥ ⎥ ⎥
⎢⎢C3(M
)
⎥ ⎥
⎢M⎥
2M ⎥
M⎦
⎢⎣CM(M ) ⎥⎦
⎡1⎤ ⎢⎢0⎥⎥ ⎢0⎥ ⎢⎢ M ⎥⎥ ⎢⎣0⎥⎦
⎡1
⎢ ⎢
13
⎢ 15
⎢ ⎢
M
⎢⎣12N −1
3 33 35
M 32 N −1
5 53 55
M 52 N −1
1992，1994年，Tessmer et al.在模拟二维以及三维不规则地表面波时，同样也是使用的如上吸收边界。
采取的措施
z 高阶差分解法--提高计算精度，减小数值频散 z 采用基于特征分析方法得到的吸收边界条件
数值频散问题------高阶差分解法
声波：
∑ ∂2 f
∂x 2
=1 Δx 2
M
Cm(M )[ f (x + mΔx) − 2 f (x) +
m=1
f (x − mΔx)] + o(Δx2M )

地下中深孔爆破地震波监测及衰减影响因素分析

地下中深孔爆破地震波监测及衰减影响因素分析
高峰;强小杰;杨潘磊;杨懿全;王巧莉
【期刊名称】《爆破》
【年(卷),期】2016(033)004
【摘要】消除及控制城市矿山井下爆破对地表环境及构筑物的振动影响具有重要意义.针对梅山铁矿爆破振动控制问题,在大量爆破振动监测的基础上,采用灰色关联度法对爆破振动影响因素开展分析.以振动合速度和主频作为系统特征指标,以单段最大药量、总药量、爆心距、高程差、分段数和炮孔数6个比较数列作为影响因素,通过灰关联度计算得到该6个指标的关联度排序.通过对实测数据进一步统计分析,将分段数控制在5段及以上对于控制地表爆破振动速度有显著效果.从最大限度控制地表振动主频的角度看,应优先使用逐孔起爆,一次爆破总药量控制在950 kg 为宜.
【总页数】7页(P27-33)
【作者】高峰;强小杰;杨潘磊;杨懿全;王巧莉
【作者单位】中南大学资源与安全工程学院,长沙410083;中南大学资源与安全工程学院,长沙410083;南京梅山冶金发展有限公司矿业分公司,南京210041;中南大学资源与安全工程学院,长沙410083;中南大学资源与安全工程学院,长沙410083【正文语种】中文
【中图分类】TD235
【相关文献】
1.在重复爆破作用下爆破地震波传播衰减规律数值模拟研究 [J], 赵禹
2.地铁车站开挖爆破地震波衰减特性监测分析 [J], 任仕国
3.延时起爆技术对中深孔爆破地震波强度的控制 [J], 付武;马玉磬;付博;王健
4.地下矿房中深孔爆破震动的监测与分析 [J], 郭学彬;肖正学
5.某地下矿山大直径中深孔爆破振动衰减规律分析 [J], 王雨波;汪为平;刘海林因版权原因，仅展示原文概要，查看原文内容请购买。

地震波传播模拟中的数值方法

地震波传播模拟中的数值方法一、引言对地球上发生的自然灾害进行研究和预测一直是人类所探究的课题之一。

其中，地震是一种造成极大灾害的自然现象，它的预测和探测对减轻地震对社会影响，提高人类对灾害的应对能力，具有重要意义。

地震波传播模拟是地震研究领域的重要课题，为了更好地预测地震和应对地震灾害，需要对地震波传播的数值模拟方法进行深入研究。

二、地震波传播数值模拟的方法1. 有限差分法（FDTD）有限差分法，英文全称为Finite Difference Time Domain，是一种常用的求解电磁场和声场传播问题的数值方法。

FDTD方法利用有限差分逼近微分算符，将偏微分方程离散化，然后通过差分方程组求解离散化问题。

FDTD方法的优点是较为简便和直观，对于一些基础场问题可以精确求解，但是FDTD方法在离散化问题域时会导致误差，对于具有复杂形状、边界不规则和含有多个介质的问题，其求解需要繁琐的预处理工作和较为复杂的网格划分，求解过程也较为复杂。

2. 有限元法（FEM）有限元法，英文全称为Finite Element Method，是一种广泛应用于工程和科学计算领域的数值方法。

它是通过将一个复杂的问题域分解成多个小问题域，用简单的数学公式在每个小问题域内求解，通过对这些小问题域的求解累加得到整个问题域的解。

FEM方法的特点是能够对不规则的计算域进行处理，求解过程较为直观和简单，对于多介质、弹性、非线性等问题也有很好的处理能力。

但FEM方法对于较为复杂的问题各向异性和自由面的处理比较困难。

3. 间接边界积分法（BEM）边界积分法，英文全称为Boundary Element Method，是近年来发展起来的一种求解偏微分方程的数值方法。

BEM方法将待求解的域分为界面和域外两部分，通过界面上的边界积分求解内部问题。

BEM方法对于不规则和异形问题的边界条件求解有很好的处理能力，并且具有较高的精度和较低的计算量。

但是对于非线性问题处理不够准确，对纯内部问题的求解效果不如其他方法。

地震学中的数值模拟技术研究

地震学中的数值模拟技术研究自古以来，地震一直是人类无法掌控的自然灾害之一。

虽然我们无法预测地震的发生，但是对于地震的研究和预防措施，可以减轻地震给人类带来的伤害。

而当今地震学领域中，数值模拟技术正得到日益重视的研究。

1.数值模拟技术数值模拟技术是一种利用计算机对实际问题进行数学模型化，并以数值计算为手段求得问题精确解的方法。

在地震学中，数值模拟技术能够模拟地震的过程，对地震的形成、发展以及危害进行研究。

通过数值模拟技术，我们可以预测地震对地表、建筑物和人员的破坏情况，为地震预防和人员疏散提供科学的依据。

2.地震数值模拟的基本步骤地震数值模拟的基本步骤分为三个部分：准备工作、模型建立和数值模拟。

准备工作包括搜集地震波数据、测量地震发生的条件和影响因素等。

模型建立需要确定模型的基本参数，如模型边界条件、材料特性、地震波输入等。

数值模拟则是利用计算机进行数值计算，得到地震波的传播、反射、衍射、能量传输和破坏情况。

3.数值模拟技术在地震学中的应用数值模拟技术在地震学中有着广泛应用，其中主要包括以下方面：3.1 地震波传播研究地震波传播是地震研究的基础，也是地震预测与震害评估的重要依据。

数值模拟技术可以对地震波的传播进行快速、准确的模拟研究，帮助我们理解地震波在地下介质中的传播规律、地震波在地面上的表现形式，以及地震波与建筑物、结构物的相互作用关系。

3.2 地震危害评估地震危害评估是对地震中各种影响因素进行评估的一项工作。

数值模拟技术不仅可以模拟地震波对建筑物和结构物的破坏情况，还可以研究地表水位的变化、坡面变形等地震带来的影响。

3.3 地震预测虽然目前没有任何一种方法可以完全准确地预测地震，但是数值模拟技术可以根据历史数据和地震破坏情况，对未来地震的可能发生地点、规模、引发危害等进行有限预测。

4. 数值模拟技术的发展前景数值模拟技术在地震学中的应用已经趋于成熟，但是仍有很多问题需要解决。

如模型精度、边界条件的处理、计算机运算速度等等都是需要进一步研究的问题。

地震勘探中的数值模拟

地震勘探中的数值模拟地震勘探是指通过部署地震探测仪器、采集地震波数据以及对数据进行处理和解释等一系列行动，来获取地下地貌、地质构造、岩层厚度、地层分布、矿产分布等信息的一种技术手段。

地震作为一种可计算的现象，因此我们可以通过数值模拟来更好地理解和预测地震勘探的效果和结果。

数值模拟是通过计算机模拟真实的物理现象，来寻求某些变量之间的关系以及对特定条件的反应，通常将物理规律以数学公式和算法的形式表达出来，然后利用计算机的计算能力求解。

在地震勘探中，数值模拟是一种重要的工具，可以对野外实测数据进行较为精确的模拟和计算，提高勘探效率和准确度。

一、数值模拟在地震勘探中的作用1. 预测地震波传播路径和覆盖面积在地震波传播过程中，地震波会在不同的介质中发生反射、折射等现象，因此地震波传播路径并不是单一的，而是呈现出复杂的波动形态。

通过数值模拟，可以在计算机中对地震波传播路径和覆盖面积进行预测，从而快速找到合适的实际测量点位，提高勘探效率和准确度。

2. 提高地震数据处理和解释的准确度在地震数据处理和解释过程中，我们通常需要更好地理解地震波在地下介质中传播的路径和规律，从而推断出地下介质的物理特征和构造情况。

而数值模拟可以提供较精确的地震波传播路径和规律预测，以及地下介质模型和物理特征参数的反演，在进行数据处理和解释时，可以帮助我们更好地理解数据意义，提高准确度。

3. 预测地震勘探效果和结果在进行地震勘探之前，我们通常需要事先进行预测，以确定合适的探测方法和点位，从而提高效率和准确度。

而数值模拟可以对采集数据进行模拟和计算，预测实际勘探的效果和结果，从而帮助我们快速找到合适的勘探目标，提高勘探效率和准确度。

二、数值模拟在地震勘探中的应用方法及技术1. 数值模拟方法数值模拟在地震勘探中的应用，通常采用声波传播方程模拟方法和有限差分法模拟方法等方法。

声波传播方程模拟方法是基于声波传播方程，根据研究对象和题目所需要的计算精度和运行速度，采用空间有限差分或积分算法进行数值求解。

地震波传播模型分析与数值模拟

地震波传播模型分析与数值模拟地震波是地震活动中的一种重要表现形式。

地震波是指地震时震源释放能量产生的波动，能以光、声、热和弹性等形式传播。

地震波传播的模型和数值模拟是研究地震学中的重要内容，也是应对地震灾害的重要手段。

地震波传播模型分析地震波传播模型是研究地震波传播规律的理论模型。

地震波传播模型的建立涉及到多学科的知识，包括地球物理学、地质学、数学等。

目前，地震波传播模型主要可以分为两种类型：解析方法和数值方法。

解析方法是利用物理学和数学知识分析和处理地震波传播的数学方程式，得出地震波传播的行为规律和传播特征，如椭球体、曲线波、双曲线波等特征。

这些解析方法主要包括：爆炸理论、时空系统理论、维克多立传感器理论、波形分析等。

数值方法是通过计算机模拟地震波的传播过程，利用差分和有限元等数值方法来计算地震波传播的各项性质。

数值方法是应对复杂地壳结构，更具有灵活性和适应性的一种方法。

这些数值方法主要包括：有限差分法、有限元法、模型元法等。

地震波传播数值模拟地震波传播数值模拟是利用数值方法计算地震波传播效果的一种方法。

数值模拟可模拟地震波在不同地质体系中的传播效果和地面运动情况。

通过数值模拟，可以得到地震波在地下的传播路径和传播速度，进而预测地震波对地表建筑物的影响。

地震波传播数值模拟主要有两种类型：二维数值模拟和三维数值模拟。

二维数值模拟通过计算地震波在水平方向上的传播情况，可以模拟较为简单的地质结构。

三维数值模拟则需要计算地震波在三维空间中的传播情况，比二维数值模拟更为复杂。

进行地震波传播数值模拟需要运用适当的计算机模型和软件。

常见的地震波模拟软件有ABAQUS、FLAC、COMSOL、MIDAS等。

这些软件可通过不同的数值方法、参数设定和建模操作，实现模拟地震波在不同地质体系中传播的效果，进而为地震防灾减灾提供参考依据。

总结地震波传播模型分析和数值模拟是研究地震学中的重要内容，能够为地震预测和防灾减灾提供依据。

三维地震数据处理中的数值模拟算法

三维地震数据处理中的数值模拟算法一、三维地震数据处理概述三维地震数据处理是地球物理学领域中的一项关键技术，它涉及到地震波在地下介质中的传播规律，以及如何通过地震数据来获取地下结构和性质的信息。

这项技术对于石油和天然气勘探、地质研究和工程勘察等领域具有极其重要的意义。

1.1 三维地震数据处理的重要性三维地震数据处理技术是勘探领域中不可或缺的工具，它能够提供地下结构的高分辨率图像，帮助地质学家和工程师更好地理解地下的地质构造、岩石类型以及流体分布等信息。

1.2 三维地震数据处理的流程三维地震数据处理包括多个步骤，从数据采集、预处理、地震波场模拟、速度建模、成像技术，到最终的解释和分析。

每一个步骤都对最终结果的准确性和可靠性起着至关重要的作用。

二、数值模拟算法在三维地震数据处理中的应用数值模拟算法是三维地震数据处理中的核心技术之一，它通过数学模型来模拟地震波在地下介质中的传播过程，从而预测地震数据。

2.1 数值模拟算法的基本原理数值模拟算法基于波动方程或弹性动力学方程，通过离散化方法将连续的地下介质转化为有限的网格系统。

然后，利用有限差分、有限元或谱方法等数值技术来求解这些方程，得到地震波在各个时间步长的波场分布。

2.2 数值模拟算法的关键技术- 波动方程求解：波动方程是描述地震波在地下介质中传播的基本方程，求解波动方程是模拟地震波传播的关键。

- 介质参数建模：介质参数如速度、密度和弹性模量等对地震波的传播特性有显著影响，准确的介质参数建模是数值模拟的基础。

- 边界条件和初始条件的设定：合理的边界条件和初始条件设定对于模拟结果的准确性至关重要。

- 并行计算技术：三维地震数据处理的数据量巨大，采用并行计算技术可以有效提高计算效率。

2.3 数值模拟算法的挑战- 计算复杂性：随着模型规模的增大，数值模拟的计算复杂性急剧增加，对计算资源的要求也越来越高。

- 多尺度问题：地下介质的多尺度特性给数值模拟带来了挑战，需要开发能够处理多尺度问题的算法。

数值模拟方法在地震预测中的应用

数值模拟方法在地震预测中的应用地震是地球上常见的自然灾害之一，给人们的生命财产安全带来了巨大威胁。

为了更好地预测和应对地震，科学家们通过数值模拟方法进行研究，以期提高地震预测的准确性和可靠性。

一、地震的成因和特征地震是地球内部能量释放的结果，通常由地壳中的断层活动引起。

地震的特征包括震源、震源深度、震级和震中等。

震源是地震发生的具体位置，震源深度是指地震发生的深度，震级是用来描述地震能量大小的指标，而震中则是地震发生地点的地理位置。

二、数值模拟方法的基本原理数值模拟方法是一种基于物理模型和数学计算的地震预测方法。

它通过建立地震波传播的数学方程，模拟地震波在地下的传播过程，从而预测地震的发生和影响。

数值模拟方法主要包括有限元法、有限差分法和边界元法等。

有限元法是一种将地震波传播过程离散化为有限个小区域进行计算的方法。

通过将地震波传播区域划分为许多小单元，建立各个小单元之间的相互作用关系，可以模拟地震波在地下的传播过程，并预测地震的发生和影响。

有限差分法是一种将地震波传播过程离散化为有限个时间步长进行计算的方法。

通过将地震波传播过程划分为一系列时间步长，分别计算每个时间步长内地震波的传播情况，可以模拟地震波在地下的传播过程，并预测地震的发生和影响。

边界元法是一种将地震波传播过程离散化为有限个边界进行计算的方法。

通过将地震波传播区域的边界划分为一系列小边界，计算每个小边界上地震波的传播情况，可以模拟地震波在地下的传播过程，并预测地震的发生和影响。

三、数值模拟方法在地震预测中具有重要的应用价值。

首先，数值模拟方法可以帮助科学家们更好地理解地震的产生机制和传播规律。

通过模拟地震波在地下的传播过程，科学家们可以研究地震波与地下结构的相互作用，揭示地震波传播的规律，为地震预测提供理论依据。

其次，数值模拟方法可以用于地震预测和灾害评估。

通过模拟地震波在地下的传播过程，科学家们可以预测地震的发生时间、震级和震中等，为地震预警和灾害应对提供重要依据。

分段计算方法在模拟爆炸应力波传播中的运用

ＡｂｔａｔＡｓｂｅｔｎｃｌｕａｉｇｎｍｅｉａｉｌｔｇｍｅｈｄｉｐｅｅｔｄａｅｎｔｅｄｅｅｔｙａｃｒｓｏｓｓｏｃｎｓｒｃ：ｕｓｃｉａｌｔｕｒｌｓｏｃｎｃｍｕａｉｔｏｒｓｎｅ．ｂｓｄｏｉｒｎｎｍｅｐｎｅｆｒｋｕ－ｎｓｈｆｄｉｏ
动力有限元方法是能有效模拟地下爆炸中菲弹性区冲击波和弹塑性波的传播方法之一且能考虑局部地质结构几何结构等对地震波的影响爆炸加载条件与真实情况比较一致但在有限的计算资源下由于岩石非弹性区动力有限元网格精度的要求单元的特征尺度过大就不能模拟冲击波的传播和计算时间步长的限制一般只能得到数百米或几千米范围内的地震波
国防科技大学学报
箜三鲞箜鱼Ｑ塑
Ｑ
Ｎ衄Ｑ
』Ｏ垦堡坐
：：Ｑ垒
文章编号：０１２８（ｏ８ｏ０８０１０ — ４６２ｏ）６— ０９— ５
分段计算方法在模拟爆炸应力波传播中的运用
严波曾新，
ｗｖ．ｈｂｅｔｎｃｃｌｔｇｍｔｄｄｉｓｈｍｕｔｇｒｇｎｉｔａｌｔｇｎｉｗｉｙａｉｆｔｅｍｎａｏｉｍａｅＴｅｓｓｉｕｉｅｏｉｄｅｓｌｉｉｏｎｅｓｃｒｏｈｈｄｎｍｃｎｅｌｅｔｌｒｈｓｕｃｏａａｎｌｈｖｅｔｉａｎｅｏｎａｉｅｉｎｃｉｅｇｔｉ
ＹＡＮ３ＩＺ１ｏＥＮＧＸｉ－２ｎｗｕ
，
（．ＣｌｇｆａｉＥｕａｏｏｆｃｒ，Ｎｔｎｌｎ．ｆＤｆｓｅｈｏｇ，ｌｎａ４０７，Ｃｉ；１ｏｅｅｏｓｄｃｔｎｆｌＢｃｉｒＯｉｓａｏａＵｉｏｅｎｅＴｃｎｌｙＣｍ￣ｈ１０３ｈｎｅｉｖｅｏａ

地下空洞地震波场的旋转交错网格有限差分数值模拟-精选文档

地下空洞地震波场的旋转交错网格有限差分数值模拟1 引言天然形成或人工挖凿形成的地下空洞，对城市基础建设具有潜在巨大隐患。

地震勘探在探测地下空洞方面取得了显著的进展，为了更好的理解地震波在地下存在空洞时传播的特征，本文采用旋转交错网格进行了高阶旋转交错网格空洞波场数值模拟，并分析了震源子波频率、地下空洞埋深及大小和形状对地震记录的影响。

在旋转交错网格一阶速度-应力波动方程波场数值模拟中，相同物理量的不同分量都被定义在单元网格中的一个相同位置，应力（质点速度）定义在单元网格的顶点，质点速度（应力）在相应的对角节点上。

质点密度和弹性常数定义在与应力相同的节点上，当模型存在物性分界面时，只需要对密度进行平均即可，因此在模拟裂隙介质和地形空洞等物性差异较大的非均匀介质时更具优势。

另外，为了减少人工边界产生的反射波的影响，在模型的人工边界外侧采用了完全匹配层。

经过分析不同模型参数下的地震记录及波场快照，加深了对地下空洞存在情况下地震波传播特征的理解。

2 旋转交错网格差分格式3 一阶应力速度弹性波方程二维各向同性介质一阶应力―速度弹性波方程可以表示为：（5）其中，为质点速度分量，为密度的倒数，为应力分量，为拉梅系数。

4 模型算例本文对埋深为2m和5m，边长为2m和5m的方形空洞模型分别在30Hz和80Hz雷克子波震源下进行了数值模拟。

空洞周围介质的纵横波速度及密度分别为866m/s、500m/s、2000Kg/m3，模型的左边界距震源40m，网格间距为0.1m，时间步长为0.00005s。

为了进一步研究震源频率、空洞大小及埋深对地震记录的影响，本文计算了含空洞和不含空洞地震记录之差，如图5所示。

模型的左边界距震源40m，最小偏移距为5m，即空洞的左边界位于检波器排列的第35道处。

5 结论通过将CZ真空法和高阶旋转交错网格相结合，较稳定而精确的模拟了地下存在空洞时的地震波传播。

通过模拟不同震源子波主频、不同深度和不同大小空洞情况下的地震记录及波场快照，分析了震源子波主频、深度及大小对地震波传播的影响。

下穿山岭隧道爆破振动效应的数值模拟分析

Science and Technology &Innovation ┃科技与创新2023年第24期·95·文章编号：2095-6835（2023）24-0095-03下穿山岭隧道爆破振动效应的数值模拟分析赵民强（中交路建交通科技有限公司，北京101300）摘要：钻爆法作为地下空间开挖的主要施工方法，爆破过程产生的地震波振动效应，是影响工程效率、危害地表构筑物的主要因素。

如何有效控制爆破产生的影响，正确处理减小爆破振动与提高施工效率之间的平衡，是隧道爆破安全、快速施工的关键问题。

依托新建崇礼铁路隧道下穿工程，通过现场实测辅以数值模拟的方法研究下穿隧道引起的环境问题，总结爆破振动的产生和传播机理。

利用数值模拟的方法，近似还原工程实际，结合上述拟得到的研究结果，分析地表建筑物不同测点在爆破振动激励荷载作用下的振动响应，以期控制隧道爆破对既有建筑物的不利影响，为类似工程提供参考意见。

关键词：隧道爆破振动；振动响应；现场监测；安全评价中图分类号：U452文献标志码：ADOI ：10.15913/ki.kjycx.2023.24.027目前，新建隧道下穿既有村庄等近接施工的实例工程日益增多[1-2]，钻爆法凭借经济高效的优势被广泛应用于隧道建设中。

作为施工工艺较复杂的隧道近接工程，新建隧道新奥法钻爆施工，开挖会打破围岩原有的平衡状态，对周围岩土体造成一定程度的扰动[3-4]，降低了地面既有建筑物的安全性，围岩变形及既有结构物沉降均难以控制，若不采取相应对策，将给生产建设和人们生命财产安全带来严重危害。

隧道掘进爆破施工诱发的负面效应，尤其是爆破振动会对周围建筑物造成不同程度的危害[5-6]。

当隧道临近区域内有建筑物存在时，为平衡新建隧道高效快速施工和周边建筑物安全问题，必须对隧道爆破振动波传播规律和建筑物受振动响应特性进行研究。

1工程概况以太锡铁路中某山岭隧道工程[7-8]为背景展开研究。

地震数值模拟技术2—1

二维一阶弹性波动方程速度X分量的差分方程
1 1 t 1 3 3 v1(m, n, l ) v1(m, n, l ) { [ xx (m , n, l ) xx (m , n, l )] 2 2 (m, n)x 24 2 2 9 1 1 t 1 3 [ xx (m , n, l ) xx (m , n, l )]} { [ xz (m, n , l ) 8 2 2 (m, n)z 24 2 3 9 1 1 t xz (m, n , l )] [ xz (m, n , l ) xz (m, n , l )]} fx 2 8 2 2 (m, n)
f ( x) x x { f ( x ) f (x )} / x x 2 2
f ( x) 1 3x 9 x 9 x 1 3x { f ( x ) f ( x ) f ( x ) f ( x )} / x x 24 2 8 2 8 2 24 2
波动方程原理
交错网格差分技术
波动方程原理
基于上式可以将各向同性介质弹性波动方程第一个方程差分为:
交错网格差分技术
提纲基于波动方程模拟技术提纲
波动方程基本理论基于波动方程地震模拟技术
声波方程模拟
速度模型设计
声波方程：
模型参数输入
S=0
T=1
空间计算波场值
波动方程原理
应力张量
波动方程原理
Z z
zy x xz yz xz O x X y dy zy z xy yx dz y Y
应力与应变分析
zx
yx
xy
y
yz
z
zx

爆炸冲击震动对砖墙破坏作用的数值模拟

好。研究表明：该种分析模型综合考虑了砖块和砂浆之间复杂的相互作用，并且对砂浆层进行了单独建模，保
证了砖墙在数值模拟上的真实性和正确性，因此可以准确地模拟出实验中砖块之间砂浆层的损伤积累破坏。
模型，即将砖块和砂浆作为独立的构成部分，考虑了二者之间的相互作用。在这种分离式模型中又可再分２种不同的模拟方法，一种是对砂浆不建模，将砂浆灰缝的厚度和质量折算到砖块上，通过砖块之间
＊
收稿日期：２０１３ — １２ — １７；修回日期：２０１４ — ０３ — ２２基金项目：国家科技支撑计划项目（２０１３ＢＡＫ０１Ｂ０１）
况，并与实验结果进行对比验证。
１砖墙的有限元模拟方法
目前，砖墙主要有２大类模拟方法。第１类是整体式模型，将其简化为一个整体来考虑，不考虑内
部砂浆与砖块复杂的相互作用，按照同一种连续化材料类型来进行分析，见图１（ａ）。第２类是分离式
关键词：爆炸力学；爆炸冲击震动；分离式共节点模型；ＬＳ－ＤＹＮＡ；数值模拟；砖墙中图分类号：０３８１国标学科代码ｉｌ３０３５文献标志码：Ａ

地震波场数值模拟方法

第42卷第2期2003年6月石　油　物　探GE OPHY SIC A L PROSPECTI NG FOR PETRO LE UMV ol.42,N o.2Jun.,2003文章编号:100021441(2003)022*******地震波场数值模拟方法张永刚(中国石油化工股份有限公司科技发展部,北京100029)摘要:简要总结了地震波场数值模拟的各种方法的基本原理及其主要特点,对最近在该领域出现的一些方法和研究结果做了简要的阐述,并对比了各种方法的优缺点。

在此基础上提出了运用波动方程数值模拟作为基础,结合射线方法辅助识别波场类型,用于分析异常波的产生机理和出现特点的基本思想,这对复杂条件下的地震勘探具有指导和借鉴意义。

关键词:地震波场;数值模拟;射线追踪;有限元;伪谱法;正演模拟中图分类号:P63114+1 文献标识码:AOn numerical simulations of seismic w avefieldZhang Y onggang(Department of Science and T echnology Development,SI NOPEC,Beijing100029,China)Abstract:This paper reviews the principles and characteristics of various numerical simulations of seismic wavefield,and com2 pares the merits and defects of the simulations.S ome newly emerged methods and results are briefly discussed.The author pro2 poses to study the generation mechanism and characteristics of abnormal waves based on wave equation numerical simulation supplemented by ray tracing.K ey w ords:seismic wavefield;numerical simulation;ray tracing;finite element;pseudo2spectrum;forward m odeling 地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段,地震波场数值模拟的主要方法包括2大类,即波动方程法和几何射线法。

爆炸中应力波理论分析及数值模拟

爆炸中应力波理论分析及数值模拟摘要：利用质量守恒定理以及动量守恒原理，对爆炸过程进行分析，推导出应力波在爆炸过程中的传播规律：应力波的幅值，波形和传播速度都会随着介质到重要中心的距离的变化而改变，并且呈现衰减趋势。

并用ANSYS模拟球形装药的应力波传播，对上述传播规律进行说明。

关键词：爆炸应力波数值模拟The Theoretical Analysis and Numerical Modeling ofExplosive Stress WaveAbstract:Analysis explosion process with the law of the law of conservation of energy and the law of conservation of mass.Propagation rule of stress wave in the explosion process is deduced.The rule suggests amplitude,waveform and wave velocity all change along with the change of media’s distance to the center of the explosion,and show a trend of attenuation. Simulate stress wave of spherical charge by ANSYS and prove the rule mentioned above.Keywords: Explosion , Stress Wave,Numerical modeling 爆炸时炸药会突然在物理和化学性质上发生巨大变化，同时伴随着巨大能量的释放，在爆炸冲击波向外传播是对周围介质进行作用，所以能够认为是应力波在介质中传播的过程。

随着介质中质点距离爆炸中心的距离的不同，应力波呈现出不同的特性，在炸药中传播的是爆轰波，附近介质中为冲击波，随距离增大变为塑性波和弹性波。