1. 引言
渗流是指流体渗入物质裂隙过程中形成的微观流动作用 [1]。裂隙在自然界普遍存在,大到可能发生于受到力学作用的岩土层,小到可以产生于增材制造等工艺下制备出的材料或元件。源于裂隙的渗流现象会改变物质原有的物理特性,对资源勘探开发、物质无损检测等技术成果的可解释性 [2] 具有广泛影响。
弹性波在物质中的传播机制是以上技术领域常常需要面对和解决的力学基本问题。当弹性波传播到物质上的裂隙时,受到渗流力学作用的影响,传播方向会发生偏离,波动性态也会有所改变 [3] [4]。为了准确描述上述现象,有必要聚焦于弹性波传播问题,构建一套行之有效的渗流力学数理方法。本文假设裂隙是作为渗流力学边界条件而存在的,首先在弹性波方程的基础上给出边界条件的若干具有直观意义的数理表达,然后绘制几种典型的渗流模型,采用频域有限元方法模拟弹性波传播到裂隙处时发生的奇特现象,揭示各种渗流力学作用对弹性波传播的影响。
2. 理论方法与基本假设
2.1. 带有完全匹配层的弹性波方程的等效积分弱形式
不失一般性,考虑二维直角坐标系(x, y)下无限大空间内的弹性波传播问题。在频域下引入完全匹配层(PML),以有限空间区域表达无限大空间内的弹性波传播。推导弹性波方程的等效积分弱形式(EIWF),以适用于频域有限元求解 [5]。EIWF的数理表达为
(1)
其中,S为一有限空间区域,L泛指区域内外边界,
为边界外法向,ρ为物质密度,ω为弹性波角频率,u、v分别为弹性波在x、y分量的位移,τxx、τyy分别为弹性波在x、y分量的法向应力,τyy为切向应力。另,sx、sy与PML有关,分别表征弹性波在x、y分量的吸收,具体设置方式详见文献 [6]。
2.2. 对裂隙的基本假设
如果物质表现为纯粹的固态,那么弹性波在物质上任一位置处的传播均满足法向位移连续、切向位移连续、法向应力连续、切向应力连续这四个边界条件;如果物质表现出固态和流态的接触作用,那么弹性波在接触处的传播满足法向位移连续、法向应力连续、切向应力为零这三个边界条件 [7]。据此,本文对物质上的裂隙性质作出如下假设:
① 裂隙宽度趋于零;
② 裂隙存在与否对弹性波传播具有实质影响,而裂隙内填充的物质对弹性波传播的影响不大;
③ 弹性波传播到裂隙处时,部分或全部法向连续性边界条件丧失;
④ 为确保弹性波穿透裂隙,切向连续性边界条件不作处理。
2.3. 边界条件的数学表达
根据章节2.2提出的假设,尤其是③号假设,给出渗流力学边界条件的三种具有直观意义的数理表达。
2.3.1. 弱渗流裂隙
弱渗流裂隙处满足的边界条件有两种情形:
① 法向位移为零、法向应力连续。此时有
(2)
其中,L0为裂隙所在位置。
② 法向位移连续、法向应力为零。此时有
(3)
其中,
为应力矩阵,
代表法向应力。
2.3.2. 强渗流裂隙
强渗流裂隙处同时满足法向位移为零、法向应力为零这两个边界条件。此时式(2)和式(3)共同成立。
3. 应用
3.1. 裂隙对弹性波传播的影响
3.1.1. 模型与参数
如图1所示绘制1号模型。物质长1.2 m,宽0.8 m,外围施加一厚度为0.2 m的PML。一条无限长的竖向裂隙将物质分割为左右对称的两部分。在1号模型下,渗流力学边界条件简化为:
① 弱渗流裂隙,法向位移为零、法向应力连续。此时根据式(2),有
(4)
② 弱渗流裂隙,法向位移连续、法向应力为零。此时根据式(3),有
(5)
在有限元求解框架下,式(5)化为
(6)
其中,vp、vs分别为弹性波在物质中传播的纵波、横波速度。
③ 强渗流裂隙,同时满足法向位移为零、法向应力为零。此时式(4)和式(6)共同成立。
在模型1中,点源和接收端处于同一水平面。点源到裂隙的距离为0.5 m。接收端一点在裂隙上,另一点在点源另一侧,到裂隙的距离也为0.5 m。
物质的弹性参数如表1所示。
Table 1. Elastic parameters of a material
表1. 物质的弹性参数
3.1.2. 结果与讨论
在点源处激发一个频率为10 kHz的u分量冲击函数,考察三种边界条件数理表达下1号模型内的u分量弹性波传播响应(以实部为例),分别如图2(b)、图2(c)、图2(d)所示。为便于比较,图2(a)展示了裂隙不发生作用的情形。同时,表2给出了以上四种情形下接收端处的u分量响应值。分析发现,法向应力为零的边界条件对弹性波传播响应的影响不大,图2(c)的波环相比图2(a)只是多了一些抖动,弹性波几乎毫无障碍地穿透到了裂隙另一侧,无法合理反映渗流力学作用;但是,法向位移为零的边界条件,会使得弹性波传播到边界条件处时产生的特性变化尤为明显,即图2(b)的波环形态(如幅度、波数等)相比图2(a)发生了更为明显的变化,这与渗流力学实际现象更加符合。当然值得注意的是,将裂隙处的法向位移强制取零,并不必然代表弹性波法向位移无法穿透裂隙传播,这可能是因为部分切向位移也在裂隙处发生了传播方向的偏转,“衍生”了法向位移,并在裂隙另一侧传播;此外,根据式(1),各个位移、应力分量的相互耦合作用也促进了这一穿透效应的产生。
(a) (b) (c) (d)
Figure 2. Effect of a crack on the elastic propagation response with a frequency of 10 kHz (In the model). (a) Without a crack; (b) With a weak permeable crack of type ①; (c) With a weak permeable crack of type ②; (d) With an enhanced permeable crack
图2. 1裂隙对10 kHz弹性波传播响应的影响(模型内)。(a) 无裂隙;(b) 弱渗流裂隙①;(c) 弱渗流裂隙②;(d) 强渗流裂隙
Table 2. Effect of a crack on the elastic propagation response with a frequency of 10 kHz (At the receiver)
表2. 裂隙对10 kHz弹性波传播响应的影响(接收端处)
3.2. 裂隙端点对弹性波传播的影响
3.2.1. 模型与参数
考虑一个更复杂些的模型。如图3所示绘制2号模型。物质长、宽、外围施加的PML厚度不变。模型内部有一条有限长的斜向45˚裂隙。鉴于章节3.1得到的重要结论,在2号模型下,仅考虑弱渗流裂隙①的边界条件,此时根据式(2),有
(7)
当
时,式(7)简化为
(8)
若执意要考虑法向应力为零的边界条件,则根据式(3),有
(9)
当
时,式(9)简化为
(10)
然后将各应力分量转换为位移即可,在此不再展开。
在模型2中,点源依旧与两个接收端处于同一水平面。裂隙一端在点源正上方,点源到该端的距离为0.5 m;裂隙另一端在最右侧接收端正下方,距离也为0.5 m。为了考察裂隙端点对弹性波传播的影响,以下将裂隙两端也视作接收端。
物质的弹性参数同如表1所示。
3.2.2. 结果与讨论
在点源处激发一个频率为10 kHz的u分量冲击函数,考察弱渗流裂隙①的数理表达下2号模型内的弹性波传播响应(以实部为例),分别如图4(a)~(b)所示。为便于比较,读者也可将图2(a)展示的裂隙不发生作用的情形拿来作为对比。分析发现,裂隙端点附近的波环发生了一定程度的扭曲,这就如同文献 [4] 所提到的,裂隙端点“似乎成为了一个新的‘波源’”,会“源源不断地在模型中传播某种波”。此外,当裂隙不发生作用时,根据点源处激发u分量冲击函数时的并矢格林函数解 [8],点源左右两侧、上下两侧的v分量本应呈现反对称特征,然而,倾斜裂隙的存在、裂隙端点的波源效应,共同使这一反对称特征受到了破坏。
(a)(b)
Figure 4. Effect of crack end point on the elastic propagation response with a frequency of 10 kHz (In the model). (a) u component; (b) v component
图4. 裂隙端点对10 kHz弹性波传播响应的影响(模型内)。(a) u分量;(b) v分量
表3给出了上述情形下各接收端处的u、v分量响应值。其中,最右侧接收端的u分量由于受到裂隙长度、位置、端点等因素的共同影响,其响应值−0.2483 − 0.3743i与1号模型中相应的0.2910 + 0.0973i
相差较大,且发生了极性翻转。顺带地,表3还给出了叠加分量
的响应值。其中,裂隙及其端
点处的叠加分量均为零,验证了模型求解的自圆性。
Table 3. Effect of crack end point on the elastic propagation response with a frequency of 10 kHz (At the receiver and the end point)
表3. 裂隙端点对10 kHz弹性波传播响应的影响(接收端、裂隙端点处)
4. 结论
本文聚焦于弹性波传播问题,对渗流力学边界条件的数理表达进行了探索,主要认识有:
(1) 将裂隙作为边界条件进行处理,不失为探究其渗流力学作用的有效策略;
(2) 对于弹性波传播问题,法向位移为零的渗流力学边界条件与实际现象更加符合;
(3) 裂隙处的法向位移为零,并不必然代表弹性波法向位移无法穿透裂隙传播。
同时,结合弹性波传播响应在不同渗流模型中的差异可知,本文提出的渗流力学边界条件的数理表达,能够合理反映裂隙长度、位置、端点等因素的共同影响。