岩体中弹性波传播尺度效应的初步分析
岩石物理学-03岩石中波的传播与衰减
D en s ity o f co n stitu en t c ry s t a l (g / c m 3) 2 . 6 5 qu a rt z 2 . 6 5 qu a rt z 2 . 6 5 qu a rt z 2 .7 1 calcite 2 .7 1 calcite 2 .1 h alite (C a, M g ) C O 32 . 8 -2 . 9 -
D en s ity (g / c m 3)
S c ree , v e g e t a l s o i l D ry s a n ds W e t s a n ds S a t u rat ed s h a l e s a n d c l a y s M a rl s S a t u rat ed s h a l e an d s a n d s ec t i o n s P o ro u s a n d s a t u ra t e d s a n ds t o n e s L im es to n es C h alk S alt A n h y dri t e D o lo m ite G ran i t e B asalt G n eis s C o al W ater Ice O il
(1)波速和密度的关系 (1)波速和密度的关系 用途: 用途:为合成地震记录提供缺失的声波或密度 曲线; 曲线;定量分析孔隙度变化引起的地震振幅变 AVO分析估算横波速度 分析估算横波速度; 化;为AVO分析估算横波速度;提供远景区的异 常岩性信息。 常岩性信息。 a.Gardner公式（1974），主要适用于沉积岩 a.Gardner公式（1974），主要适用于沉积岩 公式 ）， 中的碎屑岩。 中的碎屑岩。
成都理工大学岩石物性测试系统MTS（ System） 成都理工大学岩石物性测试系统MTS（Mechanic Test System）的岩石波 MTS 速测量方法
纵_横波在岩石中的传播速度比及弹性模量与岩石所含流体的关系
根据石西地区 50 块岩样在不同温压下取平均 值按孔隙度小于 10%及孔隙度大于 10%分别所作 直方图(图 3)。从图中可以看出 , 孔隙度增加 , 模量 常数降低 (泊松比例外), 虽然饱和水样与饱和油样 模量有重叠部分 , 但二者仍可明显区分 。尤以体积 模量 、杨氏模量 、拉梅系数区别最明显 。玛北地区 60 块样品所做的相应关系图 , 也得出了同样的结论 。
different fluid saturation in Shixi area (据 50 块样品 , 实验温度为 90℃, 压力为 90 MPa)
图 4 是岩石在不同油饱和度时 , 弹性模量随压 力的变化曲线 。可以看出 , 体积模量 、杨氏模量 、剪 切模量随油饱和度的增高而增大 , 尤其体积模量表
摘 要 据对准噶尔盆地西部井下岩心进行的实验室模拟测试 , 饱和水的岩心样品 , 其纵横波速度比(vp/vs)、 泊松比(σ)、杨氏模量(E)、和拉梅常数(λ), 均大于饱和油的岩样 , 而其体积模量(K)和剪切模量(μ),则小于饱和油 的岩样 。岩石的弹性模量还与孔隙度有关 , 孔隙度增大则模量降低 ,孔隙度减小则模量提高 。当实验压力提高时 ,则 弹性模量增大 , 实验压力降低时 , 弹性模量减小 。岩石的上述特性是有规律的 , 这为利用地球物理勘探资料检测油 、 气 、水在岩石中的分布 ,提供了依据和可能 。
O c t .1998
1.3 1.2 1.1
量
模 切
弹性波检测法在岩土测试中的应用
弹性波检测法在岩土测试中的应用随着我国经济水平的不断提高，岩土工程领域也得到一定的发展，随之日益增多的岩土测试方法也得到广泛的应用，弹性波检测法是其中一个检测技术。
本文主要研究分析弹性波检测法在岩土测试中的应用，从而介绍该方法的优越性。
标签：弹性波检测法;岩土测试;应用分析一、前言近些年来，岩土工程测试技术取得卓越的发展，弹性波理论广泛应用于地震、地质勘探、采矿、材料的无损探伤、工程结构的抗震抗爆、岩土动力学等方面，弹性波技术作为一种物探技术，在岩土测试中应用广泛，是一种高效的探测技术。
下面将进一步介绍弹性波检测法的方法以及在岩土测试中的发展和应用。
二、弹性波技术在岩土领域的发展弹性波不仅对各类岩石都有一定的穿透力和分辨力而且它在介质中传播时与介质相互作用使接收波中携带了与岩石物理力学性质相关的各种信息，所以作为一种信息载体，弹性波是最理想的。
20世纪50年代初期，日本在对岩体的调查、评价和分类中，引入了弹性波检测技术。
约在同一时期，原苏联把岩石纵波速度与岩石密度的乘积，称之为波阻抗，作为岩石爆破性的分级指标。
到7O年代，又把波阻抗和岩石的裂隙性相结合补充了岩石爆破性分级表1959年。
西方学者Maurer曾建议把波阻抗作为岩石可钻性指标。
在日本的隧道开拓工程中.广泛使用一种所谓抗爆强度为准则的岩石分级，抗爆强度就是岩石抗剪强度和岩石纵波速度的特定函数值。
7O年代以来，我国和国外水电、铁道和工程地质界相关人员，选用纵波速度进行岩体分类比较普遍。
1997年，王让甲对运用弹性波技术进行岩石可钻性分级进行了较为系统、全面的阐述和研究。
论述了运用弹性波技术进行岩石可钻性分级的可行性，对以纵波速度划分岩石可钻性分级标准的制订进行了研究，得出可喜的成果。
并主张采用动弹性力学参数综合评定岩石可钻性，其中引入了横渡的应用。
指出纵波速度和横波速度之比是一种非常实用的物理概念，也是一个评价岩石质量的有效指标，每一种岩石具有其相应的波速比值，比值的增大或减少均说明岩性的变化。
岩石中裂纹对弹性波速度的影响
岩石中裂纹对弹性波速度的影响
方华;伍向阳;杨伟
【期刊名称】《地球物理学进展》
【年(卷),期】1998(13)4
【摘要】本文在常温常压下，对岩石中的层理、裂纹或裂缝引起的弹性波速度变化特征进行了实验研究，并通过简化的人工裂缝模型，研究了裂缝密度和相对位置的变化对波速的影响，取得了一些有意义的结果．
【总页数】3页(P79-81)
【关键词】岩石;裂纹;弹性波速度;地震勘探;油气勘探
【作者】方华;伍向阳;杨伟
【作者单位】中国科学院地球物理研究所
【正文语种】中文
【中图分类】P584;P618.130.8
【相关文献】
1.岩石受力状态对弹性波传播速度的影响初探 [J], 刘盛东;杨玉岷;胡绍龙
2.岩石裂隙对岩石的弹性性质及速度-孔隙率关系的影响 [J], 周巍;杨红霞
3.弹性各向异性介质岩石体波相速度计算 [J], 蔡晓刚
4.不同温压条件下弹性波在岩石中传播速度的实验研究 [J], 马瑞;卢民杰
5.岩石受力状态对弹性波传播速度影响的初步研究 [J], 刘盛东;杨玉岷
因版权原因，仅展示原文概要，查看原文内容请购买。
结构面刚度对岩体弹性纵波传播特性影响的数值模拟
结构面刚度对岩体弹性纵波传播特性影响的数值模拟黄真萍;曾焕接;曹洋兵;陈俊熙【摘要】针对结构面刚度对岩体弹性纵波传播特性的影响规律及机理不明的问题,基于3DEC软件及自编程序,以含单个水平、光滑结构面的岩体为研究对象,分析结构面法向刚度、切向刚度对岩体弹性纵波传播特性的影响.结果表明,初至纵波振幅和纵波波速仅与结构面法向刚度kn有关,而与切向刚度ks无关;初至纵波振幅和纵波波速皆随法向刚度kn的增加而增加,增长速率逐渐减小.考虑结构面间距S和岩块杨氏模量E的影响,引入E/(S?kn)作为综合影响参量,固定kn/ks为2.5时,初至纵波振幅和纵波波速都随着(S?kn)的增加而减小,以E/(S?kn)=0.1、1、10为界,衰减速率经历极慢-变快-变慢-极慢四个阶段.【期刊名称】《福州大学学报（自然科学版）》【年(卷),期】2019(047)001【总页数】6页(P107-112)【关键词】弹性波传播;岩体结构面;法向刚度;切向刚度;离散元模拟【作者】黄真萍;曾焕接;曹洋兵;陈俊熙【作者单位】福州大学环境与资源学院;地质工程福建省高校工程研究中心,福建福州 350108;福州大学环境与资源学院;福州大学环境与资源学院;地质工程福建省高校工程研究中心,福建福州 350108;福州大学环境与资源学院【正文语种】中文【中图分类】TU459+.30 引言岩体中存在大量性质不同、规模各异的结构面，弹性波传到这些结构面时，将发生反射、透射现象，并由此导致能量、振幅、波速等发生变化.鉴于弹性波在岩体中的传播特性问题涉及岩体工程爆破、声波检测及微震监测等领域，众多研究人员对此开展大量卓有成效的研究，取得重要进展.起初，研究人员将岩体看作连续、均匀、各向同性的介质，采用等效连续介质力学方法研究弹性波传播问题.考虑到这种模型过于简化，与实际情况差距甚大，因此位移不连续模型[1-4]被提出，此模型认为弹性应力波穿过节理面时，应力场是连续的，而位移是不连续的，从而较好地模拟应力波在节理岩体中的传播特性及过程.除此之外，结构面对于岩体弹性波传播的振幅衰减、高频滤波及信号延迟等效应[4-8]也被实验揭示，研究人员基于试验结果开展相关的理论建模及力学描述研究. 由于室内试验及解析计算等研究手段的局限性，研究人员尝试利用离散元程序进行结构面对岩体弹性波传播特性影响的研究[9-15]，离散元模拟方法逐渐成为一种重要的研究方法.茹忠亮等[9]得出节理面越粗糙，法向刚度越小，弹性波衰减程度越大，弹性波的波速降低越快.卢文波[10]指出节理刚度越大，其透射系数也越大；石崇等[11]认为节理面法向(或切向)刚度越大，各点透射率越高； Myer等[1]发现节理刚度越大，透射波形越接近于入射波，刚度越小，透射波波幅减小、高频滤波.王卫华等[12-13]发现节理面粗糙度相同时，透射系数随法向刚度的增大而增大.陈勇等[14]认为一维应力波垂直节理传播时，随着节理法向刚度的增大，反射波波幅减小、频率降低，透射波的波幅增大、频率降低.综上可知，目前的研究集中在岩体结构面对应力波的透射系数、反射系数的影响及波的振幅衰减、高频滤波、信号延迟等方面，而关于岩体结构面法向刚度、切向刚度如何定量影响岩体弹性纵波波速及初至纵波振幅等问题则鲜有研究报道.本文基于3DEC离散元分析平台，通过连续介质模型检验软件及自编程序的准确性，基于单裂隙岩体模型，定量研究结构面法向刚度、切向刚度对岩体弹性纵波传播特性的影响规律，从而为相关岩体工程动力响应、安全性评价及数值模拟参数设置等提供理论支撑.1 3DEC软件及自编程序的可靠性基于3DEC软件和FISH语言，自编可精确监测模型边界动力响应的离散元计算程序，研究裂隙岩体弹性纵波传播特性.以等效连续介质模型的解析解作为标准参考值，检验3DEC软件的适用性及自编程序的准确性.图1 岩块数值模型Fig.1 Numerical model of rock block对于连续、均匀、各向同性的弹性材料，其弹性纵波传播速度Cp的理论解如下：(1)式中：ρ为材料的密度； E为材料的杨氏模量；μ为材料的泊松比.1.1 数值模拟概况构建1 m×1 m×10 m的长方体模型如图1所示，将模型的顶面及四个立面设为自由边界条件(由于底面施加纵波，侧面的自由边界条件不会有反射和折射，能达到粘性边界条件的目的)，在模型底面施加速度型正弦纵波如下式 .vz=A sin(ωt)(2)式中： vz为底面z方向的速度(m·s-1)； A为振幅，本次取0.1 m·s-1；ω为角频率，本次取200 π，即频率f为100 Hz.表1 岩块物理力学参数Tab.1 Physical and mechanical parameters of rock block材料E / GPaμρ/ kg·m-3岩块300.252700根据岩体动力响应特征，选择瑞利阻尼模型，设置最小临界阻尼比为0.05、最小临界中心频率为100 Hz，设置的岩块物理力学参数如表1所示.在模拟弹性波传播时，网格单元的尺寸对计算结果的准确性有很大影响，网格尺寸过大容易导致波形失真及波峰、波谷的位置偏离真实值，进而导致波速计算值误差过大.钱七虎院士[15]建议，网格尺寸与波长之比不应大于1/10～1/8.根据以上设置及式(1)，该岩块的纵波速度理论解为3 651 m·s-1，故输入波的波长为36.51 m，由此设置网格尺寸为0.5 m，可以较精确地模拟波在模型中的传播.1.2 可靠性论证经计算发现，同样的程序每次计算的结果都有轻微差别，这可能是数值截断误差等导致的离散性.经过大量的计算尝试，发现取6次数值计算值的平均值作为最终计算结果是稳定、可靠的.因此，本文中所有列出的数值模拟结果数据都是取6次数值计算结果的平均值.图2 纵波波速理论值与模拟值随杨氏模量E的变化 Fig.2 Theoretic calculations and numerical simulation results of P-wave velocity as a functionof the Young's modulus E按照表1的参数设置，保持泊松比μ、密度ρ不变，依次改变杨氏模量E，得出的数值计算结果与理论解对比情况见图2.由图可知，数值模拟值与理论值相当接近，两者随杨氏模量E变化的曲线也相当吻合.同样地，分别改变泊松比μ和密度ρ进行数值模拟，结果也显示弹性纵波波速理论值与模拟值较为接近，两者的相对误差皆小于2%.从而证实了程序的可行性和可靠性，并间接证明了上述数值模型尺寸、网格尺寸等程序设置的合理性，可应用于后续结构面岩体研究.2 结构面刚度影响的数值模拟基于上述被论证为合理、可靠的3DEC数值模型及自编程序，为研究结构面刚度对岩体弹性纵波传播特性的影响，在连续介质模型中部施加一条水平光滑的结构面(为接触面单元)，如图3所示.岩块的物理力学参数、边界条件、阻尼模型及参数、网格尺寸等设置情况与上节相同.图3 单结构面岩体数值模型 Fig.3 Numerical model of rock mass2.1 研究方案结构面刚度是反映结构面几何构成的函数，与结构面的风化蚀变特征、张开度、粗糙度及其吻合情况等有关.因此，平直结构面和结构面刚度可等效模拟复杂天然结构面.为研究岩体结构面法向刚度、切向刚度对弹性波传播特性的影响，需要首先确定法向刚度、切向刚度的合理取值范围.需要说明的是，3DEC软件中结构面刚度的基本单位为Pa·m-1，并非N·m-1.由结构面的泊松效应可知，结构面法向刚度与切向刚度的比值应为[1，10].一般来说，对于离散元计算，法向刚度和切向刚度应小于与该节理邻接块体的等效刚度的10倍，公式如下：(3)式中： K为块体的体积模量； G为块体的剪切模量；Δzmin为与该节理邻接块体的最小棱边长度.此要求一般仅为提高数值计算效率，对本文的简单模型而言，结构面刚度超出此范围不会显著增加计算复杂度，反而能核实相关规定的正确性. 根据数值模拟基本情况，确定出的结构面法向刚度总体取值范围为7.5～7.5×106 MPa·m-1，切向刚度总体取值范围为3.0～3.0×106 MPa·m-1，并通过以下基本思路研究结构面刚度对岩体弹性纵波传播特性的影响：1) 保持结构面切向刚度ks不变，将kn/ks由1逐渐增加至10.2) 保持结构面法向刚度kn不变，将kn/ks由1逐渐增加至10.3) 保持结构面法向刚度与切向刚度之比为2.5不变，同时改变法向刚度与切向刚度的大小.2.2 数值模拟结果分析2.2.1 结构面法向刚度的影响分析固定结构面切向刚度为最低级(3 GPa·m-1)，将kn/ks由1逐渐增加至10，即逐渐增大结构面法向刚度的大小；再将结构面切向刚度调至下一级，逐渐增加法向刚度；最终将结构面切向刚度调至最高级(30 GPa·m-1)，并同样地逐渐增加法向刚度.从而在大跨度的切向刚度ks、法向刚度kn与切向刚度ks之比kn/ks范围内，研究法向刚度对弹性纵波传播特性的影响.不同级的切向刚度、不同kn/ks比值条件下，初至纵波振幅随法向刚度kn的变化如图4(a)所示，纵波波速随法向刚度kn的变化如图4(b)所示.给定激发纵波的振幅为100 mm·s-1，若无此结构面，岩块的弹性纵波波速理论值为3 651 m·s-1.图4 不同切向刚度下初至纵波振幅和纵波波速随kn/ks的变化Fig.4 Amplitudeof first break and P-wave velocity as a function of kn/ks at variousdiscontinuity tangential stiffness由图4(a)可知，当切向刚度不变时，初至纵波振幅随kn/ks比值的增加而增大，即初至纵波振幅随法向刚度的增加而增大，但增长速率逐渐减小，最终初至纵波振幅趋近于激发纵波振幅；当切向刚度固定为3～12 GPa·m-1时，初至纵波振幅在kn/ks比值为1～3.5之间为陡增段，且切向刚度越小，增长速率越大，在kn/ks比值大于3.5后为缓增段，增长速率逐渐减小；当切向刚度固定为12～30 GPa·m-1时，初至纵波振幅在kn/ks比值为1～10之间皆为缓增段，增长速率逐渐减小.由图4(b)可知，当切向刚度不变时，弹性纵波波速随kn/ks比值的增加而增大，即弹性纵波波速随结构面法向刚度的增加而增大，但增长速率逐渐减小，最后趋于稳定，最终稳定值介于3 566～3 670 m·s-1之间，与连续介质弹性纵波波速理论值基本相当.2.2.2 结构面切向刚度的影响分析由上节可知，初至纵波振幅和纵波波速与结构面法向刚度密切相关，但无法判定其与结构面切向刚度有关.为此，在本节中，固定结构面法向刚度为最低级(1.5 GPa·m-1)，将kn/ks比值由1逐渐增加至10，即逐渐减小结构面切向刚度的大小；再将结构面法向刚度调至下一级，逐渐减小切向刚度；最终将结构面法向刚度调至最高级(30 GPa·m-1)，并同样地逐渐减小切向刚度.从而在大跨度的法向刚度kn、法向刚度kn与切向刚度ks之比kn/ks范围内，研究切向刚度对岩体弹性纵波传播特性的影响.不同级的法向刚度条件下，初至纵波振幅随kn/ks比值的变化如图5(a)所示，纵波波速随kn/ks比值的变化如图5(b)所示.由图5可知，当法向刚度不变时，初至纵波振幅和纵波波速随kn/ks的增加都几乎不变(曲线上局部微小的波动是由于数值误差引起).由此表明，结构面切向刚度对岩体初至纵波振幅和纵波波速几乎没有影响.图5 不同法向刚度下初至纵波振幅和纵波波速随kn/ks的变化Fig.5 Amplitudeof first break and P-wave velocity as a function of kn/ks at variousdiscontinuity normal stiffness2.2.3 考虑结构面间距和岩块杨氏模量的结构面刚度的影响分析由上两节可知，结构面法向刚度对岩体纵波传播特性有重要作用，而切向刚度对纵波传播几乎没有影响.但是，上述研究以岩块纵波波速为基准值，仅孤立地研究法向刚度绝对值的影响规律，忽视了结构面法向刚度与岩块杨氏模量之间的紧密关联，因而难以将岩体弹性纵波传播特性进行整体考量.并且，岩体中的结构面数量庞大，仅一条结构面无法反映工程岩体实际情况.为此，引入比值E/(S·kn)，用以衡量岩块杨氏模量E(本文设置为30 GPa)、结构面间距S(对于本文的单结构面岩体模型，S为岩体模型高度，即10 m)和结构面法向刚度kn的组合关系对岩体弹性纵波传播特性的综合影响，从而使得研究结论具有更好的普适性.结构面法向刚度kn与切向刚度ks之比为2.5时，是硬岩工程的常见工况，故本次保持kn/ks比值为2.5不变，同时改变结构面法向刚度kn与切向刚度ks的大小，研究E/(S·kn)对岩体弹性纵波传播特性的影响.初至纵波振幅Ap/激发纵波振幅A0随E/(S·kn)的变化如图6(a)所示，纵波波速vp/连续介质纵波波速理论值Cp随E/(S·kn)的变化如图6(b)所示.图6 Ap/A0和vp/Cp随E/(S·kn)的变化Fig.6 Ap/A0 and vp/Cp as a functionof E/(S·kn)1) 由图6(a)可知，初至纵波振幅Ap/激发纵波振幅A0随E/(S·kn)的增大而减小，具体的减小规律为：① 当E/(S·kn)≤0.1时，初至纵波振幅衰减速率极小，衰减幅度不超过2%，Ap/A0接近于1；② 当0.1<E/(S·kn)≤1时，初至纵波振幅衰减速率迅速增大，有陡降现象，Ap/A0介于0.69～0.98；③ 当1<E/(S·kn)≤10时，初至纵波振幅衰减速率由大变小，Ap/A0介于0.12～0.69；④ 当E/(S·kn)>10时，初至纵波振幅衰减速率越来越小，但总体衰减量越来越大，最终Ap/A0减小至0.01.用三阶多项式对图6(a)的曲线进行拟合，获得以下拟合公式：y1=-0.000 011 2x3+0.005x2-0.205x+0.994(4)式中： y1为初至纵波振幅Ap/激发纵波振幅A0比值； x为E/(S·kn).2) 由图6(b)可知，弹性纵波波速vp/连续介质纵波波速理论值Cp随E/(S·kn)的增大而减小，具体减小规律为：① 当E/(S·kn)≤0.1时，纵波波速衰减速率极小，衰减幅度不超过3%，vp/Cp接近于1；② 当0.1<E/(S·kn)≤1时，纵波波速衰减速率迅速增大，有陡降现象，vp/Cp介于0.80～0.97；③ 当1<E/(S·kn)≤10时，弹性纵波波速衰减速率由大变小，vp/Cp介于0.69～0.80；④ 当E/(S·kn)>10时，衰减速率越来越小，但总体衰减量越来越大，最终vp/Cp减小至0.67.用四阶多项式对图6(b)的曲线进行拟合，获得以下拟合公式：y2=0.000 003 16x4-0.001 4x3+0.058x2-0.28x+1.004 2(5)式中： y2为弹性纵波波速vp/连续介质纵波波速理论值Cp比值； x为E/(S·kn).3 结语1) 结构面法向刚度对岩体弹性纵波传播特性有重要作用，而结构面切向刚度对岩体弹性纵波传播特性几乎没有影响.2) 依次固定结构面切向刚度为不同级大小(3 GPa·m-1到30 GPa·m-1)，在每级取值下，逐渐增大结构面法向刚度，当结构面切向刚度不变时，初至纵波振幅和纵波波速皆随法向刚度的增加而增大，但增长速率逐渐减小，最终初至纵波振幅趋近于激发纵波振幅，纵波波速趋近于连续介质纵波波速理论值.3) 考虑岩块杨氏模量E、结构面间距S与结构面法向刚度kn的综合影响，引入比值E/(S·kn)，在保持kn/ks为2.5，同时改变结构面法向刚度kn与切向刚度ks时，初至纵波振幅和纵波波速都随E/(S·kn)的增加而减小，衰减速率经历极慢-变快-变慢-极慢四个阶段，以E/(S·kn)=0.1、1、10为界，最终Ap/A0减小至0.01，vp/Cp减小至0.67.参考文献：【相关文献】[1] MYER L R, PYRAKNOLTE L J, COOK N G W.Effects of single fractures on seismic wave propagation[C]//The International Conference on Rock Joints.Berkeley: Elsevier Press, 1990.[2] 李业学, 谢和平, 朱哲明, 等.应力波穿越分形节理时的透反射规律研究[J].岩石力学与工程学报, 2009, 28(1): 120-129.[3] ZHAO J, GAI J G.Transsmission of elastic P-waves across single fractures with a nonlinear normal deformational behavior[J].Rock Mechanics and Rock Engineering, 2001, 34(1): 3-22.[4] 王志亮, 陈强, 张宇.黏弹性节理岩体的位移不连续模型研究[J].岩土力学, 2015, 36(8): 2177-2183.[5] KING M S, MYER L R, REZOWALLI J J.Experimental studies of elastic-wave propagation in a columnar-jointed rock mass[J].Geophysical Prospecting, 1986, 34(8): 1185-1199. [6] 孙冰, 郭闪闪, 曾晟, 等.应力波在层状节理岩体中的透反射与衰减规律[J].防灾减灾工程学报, 2015, 35(6): 828-832.[7] PYRAK-NOLTE L J.The seismic response of fractures and the interrelations among fracture properties[J].International Journal of Rock Mechanics & Mining Sciences & Geomechanics Abstracts, 1996, 33(8): 787-802.[8] 赵明阶.二维应力场作用下岩体弹性波速与衰减特性研究[J].岩石力学与工程学报, 2005, 26(1): 123-130.[9] 茹忠亮, 蒋宇静.弹性纵波入射粗糙节理面透射性能研究[J].岩石力学与工程学报, 2008, 27(12): 2535-2539.[10] 卢文波.应力波与可滑移岩石界面间的相互作用研究[J].岩土力学, 1996, 17(3): 70-75.[11] 石崇, 徐卫亚, 周家文, 等.二维波穿过非线性节理面的透射性能研究[J].岩石力学与工程学报, 2007, 26(8): 1645-1652.[12] 王卫华, 严哲, 姜海涛, 等.不同形貌节理对应力波传播影响的数值模拟[J].爆破, 2017, 34(2): 34-39.[13] 王卫华, 李夕兵, 胡盛斌.模型参数对3DEC动态建模的影响[J].岩石力学与工程学报, 2005,24(增刊1): 4790-4797 .[14] 陈勇, 王志亮.节理岩体中纵波传播特性数值研究[J].佳木斯大学学报(自然科学版), 2014, 32(6): 815-817.[15] 钱七虎, 黄小平, 唐德高, 等.冲击荷载作用下有限元方法求解波动过程的精度研究[J].爆炸与冲击, 1995, 15(1): 1-10.。
岩体的动力学性质.
岩体的动力学性质岩体的动力学性质是岩体在动荷载作用下所表现出来的性质，包括岩体中弹性波的传播规律及岩体动力变形与强度性质。
一、岩体中弹性波的传播规律1、弹性波在介质中的传播速度仅与介质密度ρ及其动力变形参数Ed ，μd有关。
因此可以通过测定岩体中的弹性波速来确定岩体的动力变形参数。
2、影响弹性波在岩体中的传播速度的因素:（1）岩性：不同岩性岩体中弹性波速度不同，岩体愈致密坚硬，波速愈大，反之，则愈小。
（2）结构面：沿结构面传播的速度大于垂直结构面传播的速度。
（3）应力：在压应力作用下，波速随应力增加而增加，波幅衰减少；反之，在拉应力作用下，则波速降低，衰减增大。
（4）含水量：随岩体中含水量的增加导致弹性波速增加。
（5）温度：岩体处于正温时，波速随温度增高而降低，处于负温时则相反。
二、岩体中弹性波速度的测定可以采用地震法、声波法来测试弹性波速，下面就介绍常用的声波法。
声波法测试步骤：（1）选择代表性测线，布置测点和安装声波仪，见下图。
（2）发生正弦脉冲，向岩体内发射声波。
声波法测弹性波原理图1.发射换能器；2.接收换能器；3.放大器；4.声波发射仪；5.计时装置（3）记录纵、横波在岩体中传播的时间。
（4）根据下面的公式计算波速。
三、岩体的动力变形与强度参数1、动力变形参数动力变形参数有：动弹性模量和动泊松比及动剪切模量。
可通过声波测试确定。
优点：不扰动被测岩体的天然结构和应力状态；测定方法简便，省时省力；能在岩体中各个部位广泛进行。
计算公式：岩体与岩块的动弹性模量都普遍大于静弹性模量。
坚硬完整岩体E d/E me约为1.2～2.0 ，风化、裂隙发育的岩体和软弱岩体E d/E me约为1.5～10.0左右，大者可超过20.0。
原因如下：①静力法采用的最大应力大部分在1.0～10.0MPa，少数则更大，变形量常以mm计，而动力法的作用应力约为10-4MPa量级，引起的变形量很微小。
因此静力法会测得较大的不可逆变形，而动力法则测不到这种变形。
含孔洞缺陷的岩体力学性能研究现状
含孔洞缺陷的岩体力学性能研究现状摘要：岩体中存在各种缺陷，其中孔洞缺陷是当前正在关注的重点之一。
学者将孔洞缺陷分为两类：一类是将开挖巷道、洞室等看作孔洞；一类是将岩体内的初始损伤视为孔洞。
众多学者对含孔洞缺陷的岩体进行了研究，得出孔洞的几何特征会影响岩体力学特性的结论。
实际上，天然岩体中存在各种组合缺陷，目前对含组合缺陷的岩体研究较少，因此，为了真实评价天然岩体的工程稳定性，需展开对含组合缺陷岩体力学的试验研究。
关键词：孔洞、组合缺陷、变形、力学特性随着基础设施建设和矿产资源的迅速发展，高速铁路修建、煤矿开采、能源储存及废弃物地下处置工程规模不断扩大。
岩体作为工程组成的主要介质，很大程度决定了工程的稳定性和安全性。
自然状态下，岩体是由岩石和各种结构面，如节理、裂隙、层理、软弱夹层、孔洞等组成的复杂地质体。
由于天然岩体内部结构的复杂性，岩石的力学性能很大程度上受到这些不连续结构的影响，边坡、洞室和地基的开挖必然会使地质体发生位移、变形，严重时会导致岩体的失稳破坏[1]。
隧道、洞室等地下工程在施工建设过程中充满了大量不确定因素，这些风险因素使得地下工程在施工和运营中可能发生各种事故，且事故一旦发生，不仅会造成重大人员伤亡和灾难性经济损失，而且会对周围环境和社造成恶劣影响[2]。
因此，借助理论与试验去分析含缺陷岩体变形破坏的过程并总结这些现象发生的规律，研究发生条件、过程及机理是非常重要的。
研究缺陷岩体是国际岩石力学正在关注的重点之一。
学者在研究岩体中存在孔洞缺陷时，通常把其归为两类：一类是在岩体中进行工程开挖，把巷道、硐室、隧道等看作孔洞[3-5]，研究孔洞周围裂纹演化规律；一类是把孔洞看作岩体内部初始损伤[6-9]，结合损伤、断裂力学研究孔洞岩体中裂纹萌发与演化过程，孔洞对岩体强度及变形特征的影响规律。
不论是地质作用产生的原生缺陷，抑或是人为开挖形成的孔洞，研究内容本质上都是探究荷载作用下孔洞岩体的裂纹萌发与扩展演化机制、变形破坏特征。
岩体的弹性波频散特性及等效介质模型研究
现场测试与室内试验研究
平洞洞壁地震波与声波测试
现场测试与室内试验研究
平洞洞壁地震波与声波对穿测试
现场测试与室内试验研究 将洞壁间孔深1m范围内的岩体(对穿、单孔)平均 声波速度与洞壁地震波纵波速度进行比较，并按岩体 卸荷、风化分带进行统计。 根据69个平洞的洞壁岩体声波速度和地震波纵波 速度统计分析，结果表明各类岩体的对穿声波速度和 单孔声波速度均普遍高于同一部位的地震波纵波速度。 对穿声波速度：弱卸荷和无卸荷岩体的平均比值分别 为1.37和1.16，弱风化上段、弱风化下段和微新岩体 的平均比值分别为1.71、1.40和1.16； 单孔声波速度：弱卸荷和无卸荷岩体的平均比值分别 为1.45和1.43，弱风化上段、弱风化下段和微新岩体 的平均比值分别为1.96、1.49和1.18。
根据Hudson的微扰理论模型，弹性波由 500~1500Hz（地震波频率）提高到25kHz（声波 频率）时，纵波速度的增长幅度不超过20%。
岩体频散特性理论分析 5、考虑岩体弛豫的热 弛豫模型：
E1 E0
席道瑛等在Biot方 程中引入唯象的 复模量替换有关 固体骨架模量的 部分，在新的Biot 方程中同时引入 Biot机制和热弛豫 机制。
目录
一、前言
二、岩体频散特性理论分析
三、现场测试与室内试验研究
四、等效介质模型研究与分析
五、结语
一、 前
言
前 言 岩体是天然形成的地质体，属于一种多组分非均 匀材料的复杂结构体。岩体中的宏观结构和细观结构 对弹性波的传播均有较大的影响，使得岩体具有弹性 波频散特性，即对于不同频率的弹性波，在岩体中的 传播速度将有所差异。水电工程一般采用声波或地震 波测试岩体的弹性波速度，弹性波速度是岩体分类与 完整性评价的重要参数之一，由于岩体的频散特性使 得不同测试方法得到不同的波速值，导致岩体质量评 价结果的不一致性，这种情况在诸多工程中均有出现， 为此，我院与中国科技大学合作开展相关研究。
弹性波的传播规律及其应用研究
弹性波的传播规律及其应用研究弹性波作为一种物理现象，在地震学、声学、无损检测等领域得到了广泛应用。
本文将探讨弹性波的传播规律及其应用研究。
一、弹性波的传播规律弹性波分为纵波和横波两种。
纵波是沿着传播方向的轴向变形所引起的波动，也称为压缩波；横波是垂直于传播方向的横向变形所引起的波动，也称为剪切波。
1. 纵波的传播规律纵波的速度与材料的弹性常数和密度有关。
在同一材料中，纵波的速度越大，材料的刚度就越大。
当纵波传播到两种密度不同的介质之间时，会发生折射和反射现象。
当入射角为0度时，不存在折射现象，只有反射现象。
2. 横波的传播规律横波的速度与材料的切变模量和密度有关。
在同一材料中，横波的速度越大，材料的脆性就越大。
当横波传播到两种密度不同的介质之间时，会发生折射和反射现象。
当入射角为90度时，不存在反射现象，只有折射现象。
二、弹性波在地震学中的应用弹性波在地震学中的应用主要体现在地震勘探和地震预警领域。
1. 地震勘探地震勘探是利用地震波在地下速度差异的传播原理获取地下结构信息的一种地球物理勘探技术。
在勘探中，利用人工震源产生的弹性波，通过地下不同介质的反射、折射和衍射等效应，获取地下结构的信息。
该技术被广泛应用于石油、天然气、煤炭等资源的勘探中。
2. 地震预警地震预警是在地震发生前通过地震波速度的差异来预测地震发生的时间、地点和震级等信息。
地震预警系统通常由三个部分组成：地震监测网络、地震数据分析和地震预警系统。
其中，地震监测网络负责监测地震波，地震数据分析将监测到的数据进行处理，地震预警系统根据处理后的数据来发出警报。
该技术被广泛应用于地震灾害的预防和减轻工作中。
三、弹性波在声学和无损检测中的应用弹性波在声学和无损检测中的应用主要体现在材料缺陷检测和医疗领域中。
1. 材料缺陷检测材料缺陷检测是利用弹性波在材料中传播的敏感性来检测材料中的缺陷和损伤。
该技术能够检测材料中的裂纹、缺陷、气泡等缺陷，并获取其位置和大小等信息。
弹性波在岩石与地质工程中的应用研究
弹性波在岩石与地质工程中的应用研究引言岩石与地质工程中，弹性波被广泛应用在地下勘探、地震监测和地质灾害预警等领域。
弹性波是指在介质中传播时呈现出弹性性质的波动，能够提供关于地下结构和地质特征的重要信息。
本文将讨论弹性波在岩石与地质工程中的应用研究，并探讨其对地质工程领域的影响。
弹性波在地震监测中的应用地震是地球内部能量释放的表现，对人类生活和建筑设施造成了严重威胁。
弹性波在地震监测中起着关键作用，能够提供地震早期预警和地震波传播路径的信息。
利用弹性波监测地震活动，可以预测地震发生的时间和地点，从而采取相应的防范措施，减少损失。
此外，弹性波还可以通过分析地震波的传播路径，推断出地下构造和岩层的性质，对地质工程设计和土壤条件评估提供重要帮助。
弹性波在地下勘探中的应用弹性波在地下勘探中发挥着不可替代的作用。
地下勘探是为了获取地下信息，评估勘探区域的地质条件，以指导地质工程设计和施工。
传统的地下勘探方法如钻探和采样，受制于采样点有限和局部性，容易导致勘探结果不准确。
而利用弹性波进行地下勘探，则能够获取地下结构和地质特征的全局信息。
利用弹性波测距方法，可以通过分析弹性波在介质中的传播时间和速度，准确测定地下岩层和构造的深度和分布。
此外，利用弹性波还可以判断勘探区域中存在的地下空洞和裂隙等地质隐患，为工程施工提供重要参考依据。
弹性波在地质灾害预警中的应用地质灾害如地面塌陷、滑坡、地震等对人类生命财产造成了巨大威胁。
弹性波在地质灾害预警中的应用，有助于提前预警可能发生的灾害并采取相应措施。
地质灾害往往伴随着异常的地下运动和变形，这些变化会引起地下介质中弹性波的传播特性的改变。
通过分析监测到的弹性波信号，可以推断地下介质的变化情况，识别出潜在的地质灾害风险区域。
此外，弹性波还可以用于监测地下水位和水质变化，提前预警可能导致地质灾害的水文因素，为防灾减灾提供重要依据。
结论弹性波在岩石与地质工程中的应用研究已取得了重要成果，对地下勘探、地震监测和地质灾害预警等领域具有重要意义。
岩石弹性波传播与地震监测技术
岩石弹性波传播与地震监测技术地震是地球上最常见的自然灾害之一，也是人类生活和安全的巨大威胁。
为了更好地理解地震的本质、预测和减轻地震带来的危害，科学家们开发出了岩石弹性波传播和地震监测技术。
这些技术在地震学和地质学等领域的研究中起着至关重要的作用，本文将对此进行探讨。
一. 岩石弹性波传播岩石弹性波传播是指地震波在岩石中传播的过程。
在地震发生时，能量会以波动形式传播，即弹性波。
岩石的弹性波传播速度与岩石的物理特性有关，不同类型的岩石对弹性波的传播速度有所差异。
通过研究弹性波在岩石中传播的特性，我们可以更好地了解地震波在地壳中的行为和性质。
在岩石中，弹性波分为纵波和横波两种类型。
纵波是指沿着波的传播方向进行振动的波，而横波则是指垂直于波的传播方向进行振动的波。
这两种波在岩石中的传播速度和传播方式不同，这些差异可以被地震监测设备所感知和记录。
二. 地震监测技术地震监测技术是指用于监测地震活动和预测地震的技术手段。
在过去的几十年里，地震监测技术取得了巨大的进步。
从最早的地震仪到现在的全球定位系统和高精度测震网络，这些技术的发展为我们提供了宝贵的地震数据，帮助科学家们更好地理解地震及其发生的原理。
地震监测技术的核心是测震仪器，它们被广泛安装在地震活动频繁的地区。
测震仪器可以感知到地震波的传播，并将数据传输到地震监测中心进行记录和分析。
通过这些数据，地震学家可以研究地震活动的规律，预测可能的地震发生时间和震级，并提供相关信息给政府和民众，以便他们进行适当的防范和应对。
除了测震仪器，地震监测技术还包括地震监测网络、地震目录、震源机制研究等。
地震监测网络是由多个测震站点组成的网络，它能够提供更准确的地震数据和更全面的地震监测覆盖范围。
地震目录则是将历史地震事件的发生时间、震级和震源位置进行记录和整理，用于研究地震活动的规律和趋势。
而震源机制研究则通过对地震波的传播和振动特性进行分析，来研究地震发生的原因和地震带来的破坏模式。
岩体内应力波传播的研究进展与展望
( College of Architecture and Civil Engineering, Beijing University of Technology, Beijing 100124, China)
Abstract: The study of rock mass dynamic mechanical behavior and the properties of stress wave propagation has become one of the key scientific issues. However, discontinuities of various scales ranging from micro鄄defects to macro鄄joints widely exist in natural rock mass, which affect the propagation of stress waves and lead to the attenuation of stress wave energy. The effects of micro鄄defects and macro鄄 joints in rock mass on stress wave propagation are different in mechanism, fully understanding of the stress wave propagation in different scale discontinuities is of great significance to the development of rock mass engineering. First, the attenuation mechanism of stress wave propagation through micro鄄defected rock mass was analyzed, the wave scattering method, equivalent elastic modulus method and equivalent viscoelastic continuous medium method of stress wave propagation in micro鄄defected rock mass were described in detail, and the applicability, advantages and disadvantages of the methods were compared. Second, the transmission and reflection mechanism of stress wave propagation across macro鄄jointed rock mass was further analyzed, the displacement discontinuity method and equivalent continuum method were
两种不混溶流体饱和岩石中弹性波的传播
岩石力学与工程学报 Chinese Journal of Rock Mechanics and Engineering
Vol.25 No.10 Oct.，2006
两种不混溶流体饱和岩石中弹性波的传播
蔡袁强，李保忠，徐长节
(浙江大学 土木工程系，浙江 杭州 310027)
时，用下标“w”和“n”分别表示湿润流体和非湿 润流体。设 ui ，vi ， wi 分别表示固体颗粒、湿润流 体和非湿润流体位移， Sw 和 Sn 分别表示湿润流体 和非湿润流体的饱和度，n 表示多孔介质的孔隙率。 引入渗流力学中的 Dupuit-Forchheimer 关系式，可得 湿润流体和非湿润流体相对于固体骨架的渗流速 度[30]：
=
dρn ρn
=
dpn Kn
⎪ ⎪⎭
式中： Kw ， Kn 分别为湿润、非湿润流体的体积压 缩模量。采用体积平均的方法可获得固体颗粒所受
到的围压 p = Sw pw + Sn pn [31]。 若假设固体 骨架承受的 平均有效应 力(mean
大小可取 2.0[32]。另外，两种不混溶流体的压力通
常不连续，两者的压力差值 pc 称为毛细管力，通常 表示为饱和度 Sw 的函数：
pc = pn − pw = pc (Sw )
(6)
2.2 物性方程
孔隙流体的体积应变可分别表示为
− dε w
=
dρ w ρw
=
dpw ⎫
Kw
⎪ ⎪
⎬
(7)
− dεn
ηπ = ηπ 0 F (θπ )
(3)
式中：ηπ 0 为流体的黏性系数；F (θπ ) 为 M. A. Biot[2] 引进的流体在高频下动力黏性修正系数，且
弹性波增宽技术在评价平洞岩体中的应用
弹性波增宽技术在评价平洞岩体中的应用杜平安【摘要】论述了地震弹性波的衰减机理及衰减参数的测定方法.并结合工程实例,利用脉冲增宽技术获取的衰减参数——质量系数,评价测区岩体质量及进行围岩分类.【期刊名称】《水利规划与设计》【年(卷),期】2012(000)006【总页数】3页(P47-49)【关键词】脉冲增宽;波能衰减;质量系数;岩体评价【作者】杜平安【作者单位】陕西省水利电力勘测设计研究院西安 710001【正文语种】中文【中图分类】TV221.21 前言地震波衰减参量是颇有价值的可以用于评价岩土工程地质现状的指标之一。
实践证明，脉冲增宽技术不失为简单、实用地测定地震波衰减的手段之一。
波能衰减致使脉冲增宽，而与之相关的波的质量系数（Q）则与岩体的弹性波速、构造、风化破碎程度等存在着密切的关系，是评价岩土工程地质现状的十分重要的参数。
本文以脉冲增宽现象在测区的地震资料中所呈现出的不同变化规律及特征为主要研究对象，同时结合生产实践，进行了求取引致脉冲增宽的有关参量的计算工作。
并以此为依据，对所测区域的岩土工程地质现状作了较为客观的揭示。
2 地震波的衰减机理地震波在地下介质中传播时会逐渐衰减，其衰减机理主要有2种，分述如下。
2.1 几何扩散波在介质中传播时，波前面会越来越大，由震源所产生的相同能量，散布在面积不断增加的波前面上，会使能量密度不断减小。
因此，波的振幅就会随着传播距离的增加而逐渐衰减。
2.2 介质吸收在波动过程中，介质的不同部位之间会产生某种摩擦力，亦称内摩擦力或黏滞力。
这种因摩擦而引起的波的内损耗，必然导致波动原具有的机械能量逐渐转化为热能。
这就是介质的非弹性性质所引起的弹性波的衰减现象，亦称之为介质（地层）对弹性波的吸收。
一般来说，弹性波衰减值的大小，是随介质的性质而变化的。
所以，研究岩土对地震波的吸收现象，可从中分析岩土性质、评价岩土质量，如岩土类型、密度、结构、裂隙、孔隙度、黏滞性、非均匀性等。
隧道围岩分级及围岩压力
隧道围岩分级及围岩压力隧道所穿过的地层是千变方化的，可能遇到各种工程性质不同的围岩。
隧道围岩分级是评价隧道围岩稳定性的重要参数，也是隧道支护方案设计和施工工艺确定的主要依据。
分级的正确与否直接影响着隧道施工和运营安全，因此，正确划分隧道围岩分级就显得尤为重要。
在围岩分级确定的情况下，如何确定支护结构上的作用力（即围岩压力）就成为正确、合理设计隧道结构的关键。
4.1 围岩岩性与初始应力4.1.1 围岩岩性隧道工程围岩是指地壳中受开挖活动影响的那一部分岩土体。
这个范围在横断面上约为6～10倍的洞径。
围岩的工程性质，一般包括三个方面：物理性质、水理性质和力学性质。
而对围岩稳定性最有影响的是力学性质，即围岩抵抗变形和破坏的性能。
围岩既可以是岩体，也可以是土体。
本书仅涉及岩体的力学性质。
岩体是在漫长的地质历史中形成的地质体，被许许多多不同方向、不同规模的断层面、层理面、节理面和裂隙面等各种地质界面切割为大小不等、形状各异的各种块体。
这些地质界面称为结构面或不连续面，这些块体称为结构体，岩体可以看作由结构面和结构体组合而成的具有结构特征的地质体。
所以，岩体的力学性质主要取决于岩体的结构特征、结构体岩石的特性及结构面的特性。
环境因素，尤其地下水和地应力对岩体的力学性质影响也很大。
在软弱围岩中，节理和裂隙比较发育，岩体被切割破碎，结构面对岩体的变形和破坏都不起主导作用，所以岩体的特性与结构体岩石的特性并无本质区别。
在完整而连续的岩体中亦是如此。
反之，在坚硬的块状岩体中，由于受软弱结构面切割，块体之间的联系减弱，此时，岩体的力学性质主要受结构面的性质及其在空间的组合所控制。
由此可见，岩体的力学性质必然是诸因素综合作用的结果。
岩体与岩石相比，两者有着很大的区别：与工程总体尺度相比，岩石几乎可以被认为是均质、连续和各向同性的介质；而岩体则具有明显的非均质性、不连续性和各向异性。
岩体抗拉变形能力差，因此，岩体受拉后很容易沿结构面发生断裂。
岩体的弹性波频散特性及等效模型研究
【 摘 要 】 通 过 分析 岩 体存 在 的弹 性 波 频 散现 象 ，研 究 影 响岩 体 弹性 波 频散 特 性 的 因素 ， 并根 据 试验
室 和现 场 研 究成 果 ，采用 多 尺度 分 析方 法 ，提 出以岩 体 密度 、 泊松 比、节 理发 育 密度 、节理 张 开度 等地 质 参数 为 基础 的 单节 理等 效 介质 模 型及 频 散方 程 ，探 讨岩 体 频散 效 应对 不 同频 率 弹性 波传 播 速 度影 响程 度 的计 算 方法 ， 为工 程 中利用 不 同频 率 波速 对岩 体 进行 一 致性 评价 提 供 了可 能 。 【 关 键词 】 弹 性波 频 散特 性 等 效介 质模 型 【 ＤｏＩ 编码 】 １ ０ ． ３ ９ ６ ９ ￣ ． ｉ ｓ ｓ ｎ ． １ ６ ７ ２ ． ２ ４ ６ ９ ． ２ ０ １ ４ ． ０ ２ ． ０ ０ ４
物探 方法 与技 术
水利 规划 与 设计
２ ０ １ ４年第 ２期
岩体的弹性波频散特性及等效模型研 究
黄世 强 李广 场 徐 松 林 ２
（ １ ．浙 江华东工程安全技术 有限公 司 浙江杭州 ３ １ ０ ０ １ ４ ； ２ ．中国科 学技术大 学 北京 １ ０ ０ ０ ２ ４）
驰豫机制 的模型计算 的结果也仅为 ４ ０ ％。 显 然 ， 利 用 这 些 理 论 仍 无 法 圆满 解 释 工 程 实 际 问题 。
的矿物颗粒组 成 ， 颗粒 间的裂纹 和微结构普遍存 在 ，而且在后 期风化 、卸荷 、剥蚀 及构造运动等 地质作用 下，岩体 内部形成大量 的节理 、劈理 、 褶 皱 、断层等 ，无论细观结 构或宏观结构均属于 种 多组 分非均匀材料 ， 具 有多种微结构 的复杂 结 构 体 。岩 体 中 的宏 观 和 细 观 结 构 对 弹 性 波 的传
应力波在节理岩体中传播模型实验研究【土木工程岩土工程专业优秀论文】
广西大掌硕士掌位论文矗‘力波在节理者体中传播掴‘型实验研究１．３预研究问题的研究现状１．３．１岩石节理图１．１不连续介质Ｆｉｇ．１．１Ｄｉｓｃｏｎｔｉｎｕｉｔｙｇｅｏｍｅｔｒｙ图１．２典型解理的岩体Ｆｉｇ．１．２Ｔｙｐｉｃａｌｊｏｉｎｔｅｄｒｏｃｋｍａｓｓ播规律进行了系统的研究，运用节理非线性位移不连续模型，建立了应力波在张开节理处传播的简化模型。
此外，有效介质原理“””通常用于研究应力波在微小而密集的裂隙当中的传播，在有效介质原理体系中，将空隙以及裂隙视为有效介质，总体考虑它们对应力波的波速和振幅的衰减作用。
有效介质原理的不足之处在于忽略了各个节理对振幅衰减的独特作用且没有考虑与频率有关的振幅及振动速度。
总而言之，要得到应力波在节理面传播的正确解，关键是要适当地考虑节理面在应力波传播过程中的应变状态，因为这直接影响到节理本身的应力一应变关系。
本文将着重探讨在室内动力荷载下，压缩波在包含填充物的单一节理处的传播规律，其中应力波垂直节理面入射，节理填充物的厚度和含水量皆可调，且相对入射波长来说，节理长度很大但厚度很小。
此外，实验主要是在改装的霍普金森压杆装置上进行的。
１．４研究的目的与途径图１．３新加坡地质结构图Ｆｉｇ．１．３ＭａｐｏｆＳｉｎｇａｐｏｒｅＧｅｏｌｏｇｉｃａｌＦｏｒｍａｔｉｏｎｓ从图１．１所示的地质图可看出，新加坡国土有三分之一的面积为花岗岩所覆盖（主要集中在ＢｕｋｉｔＴｉｍａｈ和Ｇｏｍｂａｋ两个区域），调查发现““，花岗岩天然结构中存在大第二章新加坡的地质概况２．１新加坡的地质概况新加坡位于赤道以北大约８０公里处的马来西亚半岛最南端，国土面积将近６５０平方公里，受热带海洋气候影响，新加坡常年炎热、多雨、潮湿，平均气温２６．７℃，每月平均降水量约１９６ｍｍ，湿度达８４％。
新加坡的地质构成材料只要有以下四种类型㈣：１）岩浆岩，包括北部ＢｕｋｉｔＴｉｍａｈ的花岗岩和中北部Ｇｏｍｂａｋ的辉长岩；２）沉积岩，主要分布在西部和西南部的Ｊｕｒｏｎｇ地区：３）第四纪冲击层沉积土，主要分布在东部；４）新沉积土，除大部分集中分布于ｋａｌｌａｎｇ一带，还有一部分同海洋冲积残积土一样，广泛分布于整个新加坡岛。
介观尺度流体流动作用对孔隙介质弹性波传播特征的影响研究的开题报告
介观尺度流体流动作用对孔隙介质弹性波传播特征的影响研究的开题报告一、研究背景岩石是一种多孔介质，具有很高的弹性波传播特性，因此具有广泛的应用前景。
岩石中的孔隙结构影响岩石的物理、力学和化学特性，因此对孔隙介质中弹性波的传播特性进行研究具有重要的理论和应用价值。
介观尺度流体流动是影响孔隙介质中弹性波传播的重要因素之一，因此研究介观尺度流体流动作用对孔隙介质弹性波传播特性的影响具有重要意义。
二、研究目的和意义本文旨在研究介观尺度流体流动作用对孔隙介质弹性波传播特征的影响。
具体来说，本文将对孔隙介质中介观尺度流体流动对介质中弹性波传播速度、衰减和呈现出的谱特征进行定量分析和描述。
对介观尺度流体流动对孔隙介质中弹性波传播的影响进行深入的理论和实验研究，将为地球物理学、石油勘探、地震预测等相关领域的研究提供新的思路和方法。
三、研究内容1. 建立介观尺度流体流动对孔隙介质弹性波传播特征的模型和数学描述。
2. 研究介观尺度流体流动对孔隙介质中弹性波传播速度、衰减和呈现出的谱特征的影响规律。
3. 通过数值模拟和实验研究验证模型的准确性，并分析实验结果与理论分析的差异。
四、研究方法和步骤1. 建立介观尺度流体流动对孔隙介质弹性波传播特征的模型和数学描述。
2. 运用物理学和数学分析方法研究介观尺度流体流动对孔隙介质中弹性波传播特征的影响规律。
3. 进行数值模拟，研究不同参数下模型呈现出的弹性波传播特征，验证模型的准确性，并对比实验和数值模拟结果。
4. 进行实验研究，利用岩石样本进行弹性波传播实验，通过对不同参数下的实验结果进行分析，验证理论分析的准确性。
五、预期成果1. 建立介观尺度流体流动对孔隙介质弹性波传播特征的模型和数学描述。
2. 描述介观尺度流体流动对孔隙介质中弹性波传播速度、衰减和呈现出的谱特征的影响规律。
3. 利用数值模拟和实验研究验证模型的准确性，并分析实验结果与理论分析的差异。
六、研究计划和进度安排1. 第一年：建立模型和数学描述，研究介观尺度流体流动对孔隙介质中弹性波传播速度、衰减和呈现出的谱特征的影响规律。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第33卷 第9期 岩 土 工 程 学 报 Vol.33 No.9 2011年9月 Chinese Journal of Geotechnical Engineering Sep. 2011 岩体中弹性波传播尺度效应的初步分析徐松林1，郑 文1，刘永贵1，席道瑛2，李广场3(1. 中国科学技术大学中国科学院材料力学行为和设计重点实验室，安徽 合肥 230027；2. 中国科学技术大学地球及空间科学系，安徽 合肥 230026；3. 浙江华东工程安全技术有限公司，浙江 杭州 310014)摘要：含缺陷岩体具有尺度效应，此类岩体中传播的弹性波也有尺度效应。
对现场测点EC37-201-06，在3.0×3.2 m2的范围内采用动态有限元方法进行了15种尺度的弹性波传播规律的分析研究。
对现场测点EC37-101-06，在1.2×1.2 m2的范围内采用准静态有限元方法进行了60种尺度的弹性波波速与围压及计算尺度的关系的计算分析。
前者采用了射线理论分析思想，而后者采用等效介质分析思想，得到了相应的弹性波的尺度效应，但二者规律有差异。
为建立二者间的联系，也为了工程应用，基于量纲理论分析方法，给出了一个半理论的波速与入射波频率的计算公式。
与现场声波和地震波测试结果，以及考虑随机分布单节理散射模型的计算结果进行比较，初步分析结果表明，此公式基本可行。
关键词：岩石动力学；弹性波；尺度效应；节理岩体；量纲分析中图分类号：TU45 文献标识码：A 文章编号：1000–4548(2011)09–1348–09作者简介：徐松林(1971–)，男，湖北人，博士，副教授，从事材料冲击作用下响应的研究。
E-mail: slxu99@。
A preliminary analysis of scale effect of elastic wave propagation in rock massXU Song-lin1, ZHENG Wen1, LIU Yong-gui1, XI Dao-ying2, LI Guang-chang3(1. CAS Key Laboratory for Mechanical Behavior and Design of Materials, University of Science and Technology of China, Hefei 230027,China; 2. Department of Earth and Space Science, University of Science and Technology of China, Hefei 230026, China; 3. Zhejiang EastChina Engineering Safety Technology Corporation Ltd., Zhejiang 310014, China)Abstract:The propagation rules of elastic wave in rock mass with defects take on scale effect, just like the rock mass. The dynamic finite element method (DFEM) is employed to investigate the propagation rules of elastic waves at site-EC37-201-06.The whole computation area is 3.0×3.2 m2 and 15 kinds of computation scales are applied. A static finite element method (FEM) is used to study the relations of elastic wave velocities to the confined pressure and computation scales at site-EC37-101-06.The whole computation area is 1.2×1.2 m2 and 60 kinds of computation scales are applied. The ray theory is used in the former method, and the effective media theory is used in the later. The scale effect of elastic waves is obtained, but there are differences for the two methods. To establish their relations and provide a simple model for engineering computation, a semi-theoretical phase velocity equation is proposed based on the dimensionless method. Compared with the in-situ sonic velocities, seismic velocities and velocities computed by the theoretical model with randomly distributed joints, the proposed equation can be well used in rock mass.Key words: rock dynamics; elastic wave; scale effect; joint rock; dimensional analysis0 引 言作为天然的地质体，原位岩体含有大量的节理、裂隙等缺陷，岩体具有较强的尺度效应。
岩石研究一般有4种尺度[1]，即矿物颗粒尺度、岩石尺度、岩体尺度和地质尺度。
在工程应用和研究中主要涉及3种尺度：微观(micro-scale)、细观(meso-scale)和宏观尺度(macro-scale)，与上述的前3个尺度相当。
不同尺度的作用机制和研究方法不同，如微观尺度研究的是矿物间的相互作用，而宏观尺度研究的是岩体作为等效介质的响应，存在较大的差异，但是不同尺度研究之间的联系尚无定论。
本文进行弹性波传播的尺度效应的研究，拟将弹性波的波长作为不同尺度间的联系进行初步探索，另外，也可通过此研究解释现场声波和地震波数据的离散性和差异性。
岩体中弹性波传播的尺度效应研究目前尚不系统。
实验研究方面主要关注岩石颗粒、随机分布裂纹等缺陷的影响[2-7]。
只有Gettemy等[7]对比分析了San ───────基金项目：国家自然科学基金项目（40874093）；中央高校基本科研业务费专项资金项目收稿日期：2010–06–23第9期徐松林，等. 岩体中弹性波传播尺度效应的初步分析 1349Gregorio断裂带与岩块的波速，探讨了岩体尺度效应对波速的影响。
理论研究方面主要集中在两类[8]：等效介质理论和射线理论。
等效介质理论是宏观尺度的分析方法，将岩体按一定方式等效为一种连续介质。
主要有Eshelby等效介质模型、Hashin-shtrikman和V oigt-Reuss上下限理论、自洽模型和广义自洽模型、微分等效介质模型等。
此类理论模型的优点是：概念很清晰，比较适于工程分析，是一种宏观唯象的描述。
但不能很好地描述复杂岩体的局部性态。
射线理论是细微观尺度的分析方法，主要考虑波与岩体中不同形状的异质体相互作用。
Ying CF等[9]、Achenbach等[10]、Erikson等[11]分别研究了波与球状、圆柱状等散射体的作用问题。
形成了诸如波函数展开法、积分方程法（后来发展到应用Born近似）、转换矩阵法等分析方法[12]，提供了很好的理论分析平台。
此方法的优点是：可以充分考虑岩体局部的细观形态和局部的荷载作用。
但其工程使用比较困难。
两类理论分别对应“较大”的尺度（介质可以均匀化）和“很小”的尺度（介质完全无法均匀化）两种极端情况，但两者间尚无法建立有机的联系。
Gubernatis等[13]以弹性波波长与平均颗粒尺寸的比值，而马沃克等[8]在Mukerji[14]工作的基础上，以入射波长与非均匀介质散射直径的比值为尺度变量研究岩体波速随岩体尺度变化的规律，得到了两种尺度间过渡区域的特性。
这些研究揭示的岩体波速的尺度效应实际上与频散效应是耦合在一起的，尺度效应在岩体中弹性波传播的地位很难认识清楚。
由于问题的复杂性，本文首先将结合金沙江上游含节理裂隙的玄武岩岩体进行从cm量级到m量级不同尺度弹性波传播的数值分析，分别研究在“射线理论”假定和“等效介质”假定下，含复杂节理分布的岩体中波速和衰减系数随统计尺度的变化规律，初步揭示弹性波传播的尺度效应。
然后应用量纲理论的分析方法，探讨弹性波波速与岩体内部特征尺度间的关系，提出一个可用于工程分析的弹性波波速随频率变化的模型，初步实现细微观尺度和宏观尺度的联接。
1 工程背景弹性波波速的尺度效应的研究除了有理论上的意义外，也可以用以解释现场声波和地震波测试中常遇到的一些现象。
金沙江上游某水电站坝区主要赋存节理玄武岩、隐晶玄武岩、杏仁玄武岩、角砾玄武岩及含斑玄武岩等五种玄武岩[15]，对其进行声波与地震波测试发现：坝区岩体的声波速度与地震波速度之间存在较大的差异，相关关系差。
坝区强卸荷、弱卸荷和无卸荷岩体的声波速度与地震波速度比值的平均值分别为2.23，1.39和1.16，弱风化上段、弱风化下段和微新岩体声波速度与地震波速度比值的平均值分别为1.73，1.44和1.16。
随着岩体风化、卸荷程度的加剧，两者差异明显。
本文图10中的现场声波和地震波测试结果与此相对应。
坝区各种岩性均存在类似的现象与规律。
岩体声波速度和地震波速度是工程岩体质量评价的指标，二者的巨大差异给工程岩体的分级和水工建筑物的设计带来了很大的困难，影响到了工程的施工。
此处现场弹性波测试中平硐与钻孔岩体声波测试换能器主频为25 kHz，测试间距20 cm；平洞地震波测试以锤击作震源，地震波纵波主频一般为500＃1500 Hz，测试间距100 cm。
两种弹性波测试方法所采用的入射波频率不同，测试覆盖的岩体尺度有差异。
两种弹性波属性也不一样：地震波是体积波，而声波为面波。
这是测试技术上的差异，工程岩体质量评价标准中在有限的范围（<30%）进行了考虑。