１　引　言

地层孔隙结构特征作为重要的地质学参数，影响整个石油勘探开发过程，地质工程师需要利用其对地层的储量做出准确的预测，钻井工程师需要利用其预测井壁是否会失稳、发生井下事故，油藏工程师需要利用其了解渗流机理为建立开发方案提供依据，因此，地层孔隙结构的研究对整个勘探开发过程都具有重要意义．声波在固体中传播时能携带出大量内部孔隙结构特征信息而被广泛应用于地层孔隙结构的研究．目前主要利用声波的特征属性来反演地层孔隙结构，这对常规地层具有较好的适用性．而在碳酸盐岩储层［１］的研究过程中，由于该类储层饱含裂缝和孔洞系统且空间分布上存在非规律性，在形态、大小上呈现不规则性和多尺度性，这给声波实验研究造成了很大困难，使传统的超声波实验结果不能客观反映碳酸盐岩储层的固有特征．目前国内外通过室内实验统计、超声波测试［２－７］等方法对含缝洞孔隙结构的岩石进行了大量研究．魏建新［８］利用环氧树脂作为基质的裂缝模型，用面积相同的硅橡胶薄片模拟裂缝，再使用实验室人工灌注控制裂缝密度，研究了裂缝密度和形态对声波的影响；Ｈｕｓｓｅｉｎｙ等［９］研究再度证明裂缝密度的增加将导致声波速度的减小．Ｘｕ等［１０］研究表明在孔隙度一定的情况下孔洞的存在可能导致４０％的误差在纵波速度检测上；魏建新［８］、曹均［１１］、季敏［１２］等以超声波为实验手段，分别用人工物理模型模拟地层中的缝、孔洞结构，研究了对应的声波速度和衰减特性，结果表明，对速度的影响小于对衰减的影响，得到了振幅衰减属性对孔洞规模变化反应最敏感的结论，因此可将声波的衰减特性作为研究岩石孔隙结构的重要切入点．李琼等［１３］根据物理模型相似性原理，从缝洞系统中抽取并用４５钢制作了单孔洞缝物理模型，探讨单个缝洞的大小、形状与声波属性的关系；吴俊峰［１４］和季敏［１２］等利用物理模型进行了孔洞的声波响应特征研究，分析孔洞分布形态和孔洞规模对声波属性的影响，得知在缝洞处出现的“珠串”现象跟缝洞大小、形态和内部充填物有关．贺振华等［１５］测试了不同材料和不同裂缝模型的声波特征响应，得到了振幅衰减和偏振方向和裂缝夹角之间的关系；毕权贵等［１６］采用类岩石材料模型试件模拟了不同裂隙长度、裂隙厚度、裂隙倾角、非贯通裂隙对声波特性的影响．综上所述，目前国内外学者主要利用类岩石材料制作不同缝洞孔隙结构实体物理模型，通过超声波实验总结得到声波特征参数随孔隙结构变化的定性认识．由于实验条件的限制，使复杂缝洞结构的定量设置较为困难，同时，由于模型孔隙结构的多样化和研究材料的大量耗损导致无法开展大规模复杂结构的物理模型实验，使得相应的研究结果也具有一定的局限性．针对上述问题，笔者拟以裂缝系统为研究对象，从波动理论［１７－１８］出发，对超声波实验过程进行数值模拟研究，计算了缝产状、孔隙度及密度变化条件下，能够反映出岩体的各向异性特征且灵敏度高、可靠性强的声波衰减系数，并对观测结果作综合分析，进而总结声波在裂缝储层中的衰减特征．

２　超声波衰减系数测试理论声波在介质中传播时存在衰减现象，而衰减系数是用来度量声波在介质中振动或波动能量的非弹性衰减率的物理量．实验室测定岩样的声波衰减系数有三种方法：

２．１　长短岩样对比法将同一种待测岩样切割成长短不同的两块，分别测量两块样品的超声波首波峰值幅度．测量过程中，保持压力、温度及耦合条件一致．按公式（１）计算岩样的声波衰减系数：）式中：α为岩样的衰减系数，ｄＢ／ｍ；Ａ１、Ａ２分别表示两块岩样的声波幅度，Ｖ；Ｌ１、Ｌ２分别表示两块岩样的长度，ｍ．

２．２　标准样品对比法取长度与被测岩样相同的铝块作为标准样品，分别测量并记录岩样和铝块的声波首波峰值幅度．按公式（２）计算被测岩样的声波衰减系数：α＝（ｌｎＡ０－ｌｎＡ）／Ｌ＋α０，（２）式中：α０为铝块的声波衰减系数（实际计算时可近似为０），ｄＢ／ｍ；Ａ０、Ａ分别表示铝块和岩样的声波幅度，Ｖ；Ｌ表示岩样的长度，ｍ．２．３　信号对比法分别测量并记录岩样和探头对接时的首波峰值幅度．按公式（３）计算被测岩样的声波衰减系数：α＝（ｌｎＡ０－ｌｎＡ）／Ｌ，（３）式中：Ａ０、Ａ分别为探头对接和岩样的声波幅度，Ｖ；Ｌ表示岩样的长度，ｍ．本研究采用信号对比法，该方法更简洁、直观，避免了在铝块或对比岩样的声波测试过程中引入的实验误差．

３　数值模拟理论基础

３．１　波动理论用Ｕ表示某一时刻ｔ二维空间上任一点（ｘ，ｙ）处的位移，二维声波方程如下：

３．２　波动方程有限差分求解

３．２．１　初始条件当ｔ＝０时，位移速度为零．即）．（１２）在波场初始值已知的条件下，就可以通过递推公式来计算波场任何时刻的值．

３．２．２　振源条件振源的初始条件为３．２．３　边界条件假定在差分计算中，在计算区域的边界（ｘ＝０，ｘ＝ａ）上质点的位移为零，则边界条件可以定义为Ｕ（０，ｙ，ｔ）＝Ｕ（ａ，ｙ，ｔ）＝０．（１４）　　在计算区域的左右边界与实际基本相符，可以按照常规反射边界处理，而上、下边界（ｙ＝０，ｙ＝ｂ）则采用吸收边界条件［１９］：上边界：Ｕｋ＋１３．２．４　稳定性和收敛性条件当差分方程解的误差不随计算时间的推进而增加时，则该解具有稳定性．本次模拟中，对于一定的Δｘ，Δｙ，Δｔ来说，式（１０）中的ｋ取无限大时，要使差分方程的解保持在有意义的范围内，就需要差分方程的解趋近微分方程，即差分方程收敛于微分方程．利用基于傅里叶展开的Ｖｏｎ　Ｎｅｕｍａｎｎ方法，当时间步长经过ｋ次之后，差分方程与微分方程解的误差为式中ｎ为虚数，ｍ＝（０，１，２，…），ｑ表示一个时间步长上差分与微分方程的解的误差．如果ｑ≤１，差分方程的解的误差将不随时间步长的增加而增加，差分方程的解将是稳定的．根据Ｍｉ式（１７）满足差分方程，代入式（１０），得到ｑ的解为ｑ．（１９）　　当差分方程的空间和时间步长与介质速度之间满足式（１９）时，差分方程的解稳定．实验室超声波测试是采用透射法进行测量．但由于岩芯的非均质性及岩芯内裂缝分布的不确定性，给研究裂缝对声波传播特性的影响带来了较大困难．数值模拟则以被测试样的纵向剖面为模型，以探头的激发信号为振源，基于波动理论和有限差分方法用Ｍａｔｌａｂ编程来实现超声波透射实验的模拟．

４　数值模拟及分析

本文采用的圆柱体岩芯尺寸为６６ｍｍ（长度）×２５ｍｍ（直径），实体岩芯速度为６２００ｍ／ｓ，裂缝中空气速度为３４０ｍ／ｓ，提取实验室探头入射波作为初始波，保持声波传播方向不变，主频为２５０ｋＨｚ，将岩芯纵向剖面区域划分成３３０×１２５的网格，空间网格精度为０．２ｍｍ，时间网格精度为２０ｎｓ．通过不断调整岩芯的裂缝模型参数来实现孔隙结构的变化，对声波实验的模拟，其具体步骤如下：首先，载入裂缝模型参数，通过数值计算画出波场快照图；接下来，从波场模拟结果中提取出接收端波形即模拟波形（如图１）；最后，处理模拟波形中的首波幅值并按照相应的公式计算出衰减系数．

４．１　裂缝孔隙度变化对声波衰减系数的影响裂缝产状一定，通过改变缝宽（张开度）、缝长、裂缝数都可以实现裂缝孔隙度的变化．

４．１．１　缝宽对声波衰减系数的影响水平缝长１１ｍｍ，缝宽分别设为０和２０ｍｍ，记为Ａ０—Ａ５．声波激发方向与裂缝延伸方向垂直．由于裂缝内空气介质和岩芯骨架的声阻抗差距大，从波场快照图２中可以明显看出：超声波在裂缝边缘产生反射、散射，导致声波衰减．随缝宽增大，越来越多的波信号被反射回岩芯上边界，声波透射能力逐渐减弱，同时，从岩芯右壁散射到接收端的信号因裂缝右边界的反射而受到阻碍，且随缝宽变大，第二次反射的界面范围增大，信号在两界面间不断发生反射，使得声波传播的路径也变长，也导致接收端信号的能量越来越弱．从图３可知：当声波激发方向与裂缝延伸方向垂直时，随缝宽Ｌ增大，声波衰减系数α呈α＝ａｅｄＬ函数规律递增，对主频为２５０ｋＨｚ的纵波，当缝宽大于１ｍｍ，衰减系数变化幅度逐渐减小．

４．１．２　缝长对声波衰减系数的影响裂缝延伸方向仍然与声波激发方向垂直，裂缝起点位置变化：（１）裂缝以岩芯侧面一端为起点向质心扩展．缝宽度１ｍｍ，缝长分别设为０、１、２、４、５和１０ｍｍ，记为Ｂ０—Ｂ５．（２）裂缝以岩芯的质心为起点向岩芯两侧扩展．缝宽度１ｍｍ，缝长分别设为０、１、２、４、５和１０ｍｍ，记为Ｂ１０—Ｂ１５．从图６、图７见，在两种缝分布状态下，随缝长变大，岩芯衰减系数都呈明显的上升趋势．并且缝分布对衰减系数也有较大影响，缝在边缘时，衰减系数随缝长增加呈指数函数递增；而缝位于岩芯质心处时，衰减系数随缝长呈线性趋势增加．从波场快照（图４、５）中可以看出：缝宽恒定时，缝长增加，阻碍传播路径的反射界面范围也越来越大，到达接收端的声波信号逐渐减少，导致了衰减系数的变大．

４．１．３　裂缝数量对声波衰减系数的影响缝长１１ｍｍ，缝宽１ｍｍ，水平缝条数分别设为１、２、４、５、１０、２０条，记为Ｃ１—Ｃ６．随裂缝数量增加，裂缝孔隙度变大，衰减系数线性增加（如图９）．单位长度岩芯骨架中气饱和裂缝条数增加，两相分布更不均匀，由于岩芯骨架与气相的动力模量相差很大，在骨架与气相之间产生一定的应变相位差，从而在两相界面处产生内摩擦而使机械能转化为热能，较大地消耗能量；同时，随着裂缝数量急剧增加，声波传播的路径增多（如图８），也将导致声波衰减系数增大．

４．２　裂缝密度变化对衰减系数的影响在孔隙度恒定的情况下，通过以下两种手段来模拟缝密度变化对衰减系数的影响：（１）裂缝延伸方向与声波激发方向垂直时，缝长）．　　（２）裂缝延伸方向与声波激发方向一致时，缝长同为２０ｍｍ，在孔隙度不变的条件下，缝宽分别为１２，６，４，３，２，１ｍｍ（记为Ｄ１１—Ｄ１６）．　　从图１２、图１３见，在孔隙度恒定的条件下，当组成单个孔隙空间的裂缝尺寸变小时，缝数目增大，缝密度也变大，两种缝产状的岩芯衰减系数都呈明显的幂函数递增．从波场快照（图１０、１１）中可以看出：当缝宽逐渐变小时，骨架与空气相之间的界面个数也相应地增加，产生更多的反射、散射，使声波穿透的能量减小，岩石声波的衰减系数增大；当开度逐渐变小时，传播路径上裂缝个数急剧增加，声波传播的路径相应地增多，而该声波传播主能量路径的宽度却越来越小，受到不同裂缝边缘之间的反射、折射干扰就会越来越明显，消耗了部分能量，导致了该衰减系数显著增大．

４．３　裂缝产状变化对衰减系数的影响裂缝实体模型见图１４．缝长１１ｍｍ，缝宽８ｍｍ从图１５可以看出，当声波传播方向与裂缝延伸方向的夹角在１０°～９０°范围内时，随着裂缝角度的逐渐增大，衰减系数线性递增．由于缝角度的增加其阻碍声波传播的裂缝横截面积增大，被反射回岩芯顶端的信号越来越强，因此，使超声波散射的能量增多，穿透的能量减小，导致了岩芯的衰减系数增大．

５　结　论

（１）随孔隙度增加，衰减系数总体呈上升趋势但孔隙度增长方式不同，衰减系数的递增函数规律亦不同．缝宽增大，声波衰减系数呈α＝ａｅｂＬ＋ｃｅｄＬ函数规律递增；衰减系数随边缘缝长增加呈指数函数递增；衰减系数随中间缝长呈线性趋势增加；裂缝数量增加，衰减系数线性增大．（２）在孔隙度不变的前提下，随裂缝密度增大，衰减系数以幂函数形式递增．（３）当孔隙度恒定时，裂缝产状对岩芯的声波衰减系数影响很大，随着声波激发方向与裂缝延伸方向夹角角度的增大，衰减系数呈线性趋势增加．（４）该方法在含裂缝岩芯的超声波透射实验模拟中取得较好效果，能避免在实验操作过程中出现的人为误差，节约取芯成本，提高实验效率，为缝洞结构的岩芯声波实验研究提供了新的思路．

论文写好了，要发表到哪里去，这可是件大事，可以在本站了解一下评职的具体要求哦，可能与你单位的评职要求相吻哈哦，可咨询我们哦。