岩体的弹性波频散特性及等效介质模型研究

岩体的弹性波频散特性,及等效介质模型研究,报告人：黄世强,2012,年,9,月,24,日,中国水力发电工程学会地质及勘探专业委员会,中 国 水 利 电 力 物 探 科 技 信 息 网,2012,年会交流论文,目录,二、岩体频散特性理论分析,三、现场测试与室内试验研究,一、前言,四、等效介质,模型,研究与分析,五、结语,一、 前 言,前 言,岩体是天然形成的地质体，属于一种多组分非均匀材料的复杂结构体。岩体中的宏观结构和细观结构对弹性波的传播均有较大的影响，使得岩体具有弹性波频散特性，即对于不同频率的弹性波，在岩体中的传播速度将有所差异。水电工程一般采用声波或地震波测试岩体的弹性波速度，弹性波速度是岩体分类与完整性评价的重要参数之一，由于岩体的频散特性使得不同测试方法得到不同的波速值，导致岩体质量评价结果的不一致性，这种情况在诸多工程中均有出现，为此，我院与中国科技大学合作开展相关研究。,前 言,岩体的结构体形式,目录,二、岩体频散特性理论分析,岩体频散特性理论分析,大量的研究成果资料表明，岩体的弹性波频散特性客观存在，但无法利用经典的连续介质弹性波理论进行解释。近年来国内外学者针对岩体的频散特性开展研究，相继提出基于复合材料细观力学的弹性波粘弹性理论：基于孔隙介质液体流动的,Biot,理论模型，考虑喷射机制的,BISQ,理论模型，考虑岩体裂纹的,Hudson,微扰理论模型和考虑岩体弛豫的热弛豫模型等，或研究岩体中弹性波的散射问题，试图利用这些理论或模型解释岩体的弹性波频散特性。,岩体频散特性理论分析,1,、基于复合材料细观力学的弹性波粘弹性理论：,当弹性波频率由,5001500Hz,（地震波频率）提高到,25kHz,（声波频率）时，根据粘弹性理论模型计算的波速增长不超过,20%,。,岩体频散特性理论分析,2,、基于孔隙介质液体流动的,Biot,理论模型：,基本假定：,波长远远大于所研究的宏观体积单元。,流体和固体相的位移均为小位移。,所研究的液相是连续的。,假设基质是各向同性的。,不考虑热力耦合等耦合现象。,低频时，孔隙中的流体为泊松流体，高频时需作与孔隙形状有关的校正。,岩体频散特性理论分析,3,、考虑喷射机制的,BISQ,理论模型：,认为：,Biot,机制是流体受固体骨架作用发生与波传播方向同向的运动引起的，同时考虑由于弹性波对细小孔隙的喷射压榨导致垂直波传播的方向的流体流动。,岩体频散特性理论分析,根据,Biot,模型和,BISQ,模型，在低频时，纵波速度随频率的变化幅度均不大。,Biot,模型频散效应只在临界频率附近狭窄的区域有变化，但幅度约在,1%,；,BISQ,模型频散效应在高于临界频率附近狭窄的区域较明显，幅度约在,10%,。,Log,（,f,）,Vp,岩体频散特性理论分析,4,、考虑岩体裂纹的,Hudson,微扰理论模型：,Hudson,用等效介质理论，考虑具有较长波长的波在破碎岩体中的传播问题，得到由于裂隙的存在波速和弹性常数的表达式：,根据,Hudson,的微扰理论模型，弹性波由,5001500Hz,（地震波频率）提高到,25kHz,（声波频率）时，纵波速度的增长幅度不超过,20%,。,岩体频散特性理论分析,5,、考虑岩体弛豫的热弛豫模型：,席道瑛等在,Biot,方程中引入唯象的复模量替换有关固体骨架模量的部分，在新的,Biot,方程中同时引入,Biot,机制和热弛豫机制。,E,1,h,E,0,热驰豫机制模型可以在全频率域表现出较明显的频散特性，最大增长幅度可达到,40%,左右。但该模型在考虑岩体局部细观结构时有困难。,岩体频散特性理论分析,根据现有各种弹性波频散理论模型计算，不同频率弹性波的速度差异并不大，不能完全解释岩体的频散现象。最新的研究成果表明，岩体还具有弹性波尺度效应特征。,Gettemy,等,(2004),对比研究了,San Gregorio,断裂带与岩块的波速，分析不同岩块尺度对波速的影响；,Mukerji(1995),、,Kinra(1982),等通过复合结构模型实验来研究弹性波的尺度效应，采用不同频率的弹性波测试其传播速度，试图揭示弹性波尺度效应的现象。研究结果表明，复合结构的弹性波速度和波长与模型尺度比值有较大的关联性，高低频弹性波速度的最大差异达到,70%,左右。,复合结构模型的尺度效应图,（,Mukerji,）,岩体频散特性理论分析,岩体频散特性理论分析,弹性波频散理论模型计算和尺度效应研究成果均表明，岩体的频散特性客观存在，岩体在不同频率的弹性波激励下，其传播速度是不同的。在工程实践中，岩体的声波速度往往高于地震波速度，对于新鲜、完整的均质性较好的岩体，其差异一般为,7%12%,；而对于风化、卸荷岩体，其差异要大得多，甚至达到一倍。若分别根据岩体的声波速度和地震波速度评价岩体的完整性、风化程度或围岩类别，将导致不同的评价结果。,目录,三、现场测试与室内试验研究,现场测试与室内试验研究,某水电工程坝区出露以玄武岩为主，主要有斜斑玄武岩、隐晶质玄武岩、柱状节理玄武岩、杏仁状玄武岩等，在前期勘察阶段对坝区,69,个平洞的洞壁进行钻孔声波测试和地震波测试，钻孔声波测试使用主频为,25kHz,的声波换能器，分别进行单孔声波测试和对穿声波测试，钻孔深度一般为,1.5m,，孔间距为,1m,。洞壁地震波测试则以锤击作震源，分别测试洞壁岩体的纵波速度和横波速度，地震波纵波的主频一般为,5001500Hz,，检波距为,1m,。,洞壁声波测试钻孔与地震波测线重合，检波点位于钻孔孔口。,1,、现场岩体测试,平洞洞壁地震波与声波测试,现场测试与室内试验研究,平洞洞壁地震波与声波对穿测试,现场测试与室内试验研究,将洞壁间孔深,1m,范围内的岩体,(,对穿、单孔,),平均声波速度与洞壁地震波纵波速度进行比较，并按岩体卸荷、风化分带进行统计。,根据,69,个平洞的洞壁岩体声波速度和地震波纵波速度统计分析，结果表明各类岩体的对穿声波速度和单孔声波速度均普遍高于同一部位的地震波纵波速度。,对穿声波速度：,弱卸荷和无卸荷岩体的平均比值分别为,1.37,和,1.16,，弱风化上段、弱风化下段和微新岩体的平均比值分别为,1.71,、,1.40,和,1.16,；,单孔声波速度：,弱卸荷和无卸荷岩体的平均比值分别为,1.45,和,1.43,，弱风化上段、弱风化下段和微新岩体的平均比值分别为,1.96,、,1.49,和,1.18,。,现场测试与室内试验研究,平洞洞壁岩体,(,对穿,),声波与地震波纵波速度比值统计表,岩性,弱卸荷,无卸荷,弱风化上段,弱风化下段,微新,柱状节理玄武岩,1.30,1.15,1.40,1.30,1.14,微晶隐晶玄武岩,1.50,1.15,1.95,1.46,1.15,斜斑玄武岩,1.36,1.16,2.14,1.42,1.17,杏仁,状玄武岩,1.31,1.19,1.35,1.40,1.20,平,均,1.37,1.16,1.71,1.40,1.16,现场测试与室内试验研究,平洞洞壁岩体,(,单孔,),声波与地震波纵波速度比值统计表,表,2-1,平洞洞壁岩体声波与地震波纵波速度比值统计表,岩性,弱卸荷,无卸荷,弱风化上段,弱风化下段,微新,柱状节理玄武岩,1.36,1.13,1.92,1.37,1.12,微晶隐晶玄武岩,1.61,1.15,2.19,1.56,1.14,斜斑玄武岩,1.44,1.11,1.92,1.53,1.15,杏仁,状玄武岩,1.39,1.18,1.79,1.49,1.18,平,均,1.45,1.43,1.96,1.49,1.18,现场测试与室内试验研究,由此可见，无论单孔声波或对穿声波，岩体声波速度与地震波速度的差异较大，且随着岩体风化、卸荷的加剧，岩体声波速度与地震波速度之间的差异愈加明显，弱风化、卸荷段岩体的声波速度比地震波速度高一倍左右。考虑声波的传播路径，洞壁对穿声波的传播路径与洞壁地震波传播路径更加相近，两者的波速更具有可比性。,对于微新岩体，对穿声波速度与地震波纵波速度的比值范围为,1.141.20,，平均值均为,1.16,，说明节理相对不发育的微新岩体的弹性波频散效应相对较弱。当然，现场测试结果所表现的岩体频散现象是岩体频散效应和尺度效应的综合反映。,现场测试与室内试验研究,为研究玄武岩岩块的频散特性，在现场拾取的柱状节理玄武岩岩块进行环境扫描电镜试验，在不同尺度观察玄武岩的细观和微观结构，在,50,m,以上尺度观察时，岩样外表比较致密；在,50,m,以下尺度观察时，可观察到原始结构面和微裂缝。,2,、岩块室内试验,现场测试与室内试验研究,不同尺度下破裂面的微观形貌,200,m,20,m,10,m,现场测试与室内试验研究,不同尺度下观察玄武岩微观形貌,50,m,尺度,10,m,尺度,现场测试与室内试验研究,为研究玄武岩岩块的弹性波频散特性，制作一块玄武岩试件，试件尺寸为,41mm41mm21mm,。分别采用,25 kHz,、,50 kHz,、,400 kHz,、,600 kHz,和,1000kHz,频率的声波测试其声波纵波速度，为,4700m/s,、,4750m/s,、,5200m/s,、,5500m/s,和,6000m/s,；当弹性波频率从,25kHz,增加到,1000kHz,时，玄武岩岩块的波速约提高,27.6%,，呈现明显的频散特征。,现场测试与室内试验研究,玄武岩岩块声波频率,声速关系图,现场测试与室内试验研究,目录,四、等效介质模型,研究,与分析,等效介质模型研究与分析,无论现场岩体测试，还是室内岩块试验，结果均表明，玄武岩岩体（或岩块）的弹性波速度与频率相关，频率越高，波速越高，频散特征愈加明显。,根据试验室和现场研究成果，考虑岩体的细观和宏观结构特性，以岩体的密度、泊松比、节理发育率等地质参数，建立岩体单节理等效介质模型，采用多尺度分析方法研究岩体的弹性波频散效应（包含尺度效应）。在波动方程的基础上，引入裂纹面自由边界条件，采用基本奇异格林函数，结合边界积分方法，研究弹性波入射岩体时内部诸如节理、裂隙等不连续结构面对弹性波的散射作用，利用频散方程计算不同频率的弹性波速度，提出整个频率域上的频散方程：,等效介质模型研究与分析,利用上述频散方程，根据岩体的地质参数，获得不同参数条件的频散特性曲线，用于计算岩体的不同频率的弹性波速度。,典型的单节理等效介质模型，弹性波频率从,50Hz,到,100kHz,，弹性波速度总体上有增大趋势，呈现明显的频散效应，而且节理密度越大，频散效应越明显；在弹性波频率从,1kHz,至,10kHz,曲线段，弹性波速度增长速率最大，表明此频率段的岩体频散效应愈加明显 。,等效介质模型研究与分析,典型单节理等效介质模型频散图,表,2-1,平洞洞壁岩体声波与地震波纵波速度比值统计表,（,图中,n,为节理密度，,u,为泊松比）,等效介质模型研究与分析,为验证单节理等效介质模型的符合性，在,PD133,、,PD921,、,PD923,、,PD931,、,PD933,平洞内选取,6,个测段进行验证分析。测段的岩性有斜斑玄武岩、隐晶质玄武岩和柱状节理玄武岩等，包括不同风化、卸荷带的岩体，实测岩体地震波纵波速度范围为,41005400m/s,，声波速度范围为,45255698m/s,，完整岩块的波速为,6450m/s,，根据平洞地质测绘资料，各测段的节理发育情况见下表。,等效介质模型研究与分析,各测段洞壁地质测绘成果表,洞号,洞深,(m),迹长大于,10cm节理,迹长小于,10cm节理,条,/m,长度,(m),张开率,张开,(mm),条,/m,平均长度,(cm),张开率,张开,(mm),PD921,95102,3,1,多闭合,/,10,10,多闭合,/,PD923,110113,5,1,多闭合,/,20,10,多闭合,/,P,D133,2933,13,1,多闭合,/,80,8,多闭合,/,8894,10,1,多闭合,/,50,8,多闭合,/,P,D933,174178,6,1,多闭合,/,20,10,多闭合,/,P,D931,104107,7,1,多闭合,/,60,10,多闭合,/,等效介质模型研究与分析,根据单节理等效介质模型，分别计算各测段的岩体声波速度（,25kHz,）和地震波纵波速度（,1000Hz,），与现场实测值进行对比，其中地震波速度的相对误差范围为,-6%+27%,，平均相对误差为,+13%,；声波速度的相对误差范围为,-7%+23%,；平均相对误差为,+11%,。由此表明，利用单节理等效介质模型计算的结果与实测值基本符合。,等效介质模型研究与分析,模型计算波速与实测波速比对分析表,等效介质模型研究与分析,目录,五、结 语,结 语,通过现场岩体测试和室内岩块试验研究玄武岩岩体的弹性波频散特性，采用多尺度分析方法，提出以岩体密度、泊松比、节理发育密度、节理张开度等地质参数为基础的单节理等效介质模型，探讨岩体频散效应对不同频率弹性波传播速度影响程度的计算方法，为在工程中利用不同频率弹性波速度对岩体进行一致性评价提供了可能。,通过对弹性波测试资料的地质分析发现，地震波速度与包含地质结构的岩体关联性较好，而声波速度与岩块质量关联性较好，两者从不同方面反映宏观岩体的宏观质量和岩块的细观质量。,结 语,由于岩体的弹性波频散特性和尺度效应影响因素多而复杂，见诸于公开发表的相关研究并不多。对于岩体的频散问题一直是困扰物探与地质工作者的现实问题，希望本文起到抛砖引玉的作用，推动有更多开创性的研究与成果应用。,感 谢！,