Simulation study of pressure wave fluctuation law along fracture under water hammer effect
-
摘要:
为了分析水锤效应对裂隙内脉动水压力传播规律的影响，基于一维瞬变流模型与水锤运动方程和连续性方程，通过数值模拟研究流体运动状态瞬时变化时水锤压力的跃升及波动规律。模拟结果表明：脉动水流在裂隙中以一定流速流动时，当遇到裂隙末端闭合时，流速受阻降低，受惯性和压力波叠加反射的作用会发生水锤效应；在水锤作用下流体压力会瞬时跃升并且以矩形压力波的形式冲击煤体裂隙；当压力波交替升降作用于煤体裂隙内部时，激增的流体压力有利于解堵及裂隙扩展，形成水锤正反馈；而因压力波的叠加反射，不同流段的流体会相互阻碍导致压力降低甚至变为负压，形成水锤负反馈；负反馈导致冲击裂隙的压力降低，但形成的周期性交变压力可以对煤体形成疲劳损伤，并衍生出更多新裂隙，提高压裂效果。
Abstract:In order to analyze the effect of water hammer effect on the propagation law of pulsating water pressure in the fissure, based on the one-dimensional transient flow model with water hammer equation of motion and continuity equation, the jump and fluctuation law of water hammer pressure when the fluid motion state changes instantaneously is studied by numerical simulation. The simulation results show that when the pulsating water flows at a certain flow rate in the fissure, the flow rate is hindered and reduced when it encounters the closure of the fissure end, and the water hammer effect occurs due to the superimposed reflection of inertia and pressure waves. Under the action of water hammer, the fluid pressure will jump instantaneously and hit the coal fissure in the form of rectangular pressure wave. When the pressure wave alternately rises and falls inside the coal fissure, the surge fluid pressure is beneficial to unblock and fissure expansion, forming a positive water hammer feedback; while the superimposed reflection of the pressure wave, the fluid in different flow sections will obstruct each other, resulting in pressure reduction or even negative pressure, forming a negative water hammer feedback. The negative feedback causes the pressure of impact fracture to decrease, but the formed periodic alternating pressure can form fatigue damage to the coal body and derive more new fractures to improve the fracturing effect.
-
在煤矿井下，煤层瓦斯压力测定（比如胶囊-黏液封孔等）和钻孔瓦斯涌出初速度测定均需要向钻孔内送入封孔器。当封孔器在钻孔内移动时，封孔器端头与钻孔残留煤渣相互作用；一方面，钻孔残留煤渣会阻碍封孔器移动，另一方面，受封孔器端头力学作用，煤渣出现堆积现象；如果堆积严重，封孔器难以快速到达钻孔预定位置，参数测定可能失败（比如钻孔瓦斯涌出初速度测定，要求2 min内将封孔器送入预定位置并密封钻孔）。
为减轻钻孔残留煤渣对封孔器移动的影响，学者们从排渣和减阻2个方面进行了研究。在排渣方面：刘涛[1]提出了环缝式引射器输送钻孔煤渣技术；王永龙等[2-3]提出了双动力低螺旋钻杆排渣技术和刻槽钻杆排渣技术；杜长胜等[4]提出了脉动通风与钻杆分风联合作用排渣技术，李国强等[5]提出了旋风分离除渣技术；李栋等[6]提出了射流排水排渣技术；张辉等[7]提出了泵吸反循环钻进排渣技术；潘竞涛等[8]提出了下行定向钻孔氮气泡沫幂律多相流携渣解堵技术；汪义龙等[9]提出了煤层定向钻孔气举排渣技术；李成成[10]提出了“逐级升压切缝、低压旋转排渣”的排渣技术；张宏杰等[11]提出了“高螺旋机械排渣为主、液动排渣为辅”的复合排渣技术；童碧等[12]提出了基于整体式螺旋钻杆的复合排渣工艺技术；何将福等[13]提出了三通道反循环排渣技术；张宏钧等[14]提出了螺旋钻杆与压风复合排渣技术。在减阻方面：齐黎明等[15]提出了在封孔器前端安设锥面端头的减阻技术；周福宝等[16]提出了让高压气流从钻孔孔壁与抽采管之间流出，从而减轻煤渣阻碍抽采管移动的技术方案。
如果排渣技术过度使用，会扰动钻孔孔壁并产生新的煤渣甚至于塌孔；由此可见，排渣技术虽然可降低钻孔残留煤渣数量，但不能将残留煤渣清除；因此，仅采用排渣技术是不够的，需补充应用减阻技术。在封孔器前端安设锥面端头可以有效降低残留煤渣阻力[15]，是一种简单、实用的减阻技术；封孔器锥面端头（简称锥面端头）与钻孔残留煤渣的力学作用及其减阻机理是该技术推广应用中亟待解决的基础理论问题，并且锥面端头设计也需要进一步优化，在公开的文献资料中，未见相关研究成果。为此，从理论上分析锥面端头与钻孔残留煤渣的力学作用，探索锥面端头减阻机理，并进行工程应用。
1. 锥面端头前方钻孔残留煤渣堆积状况分类
锥面端头前方残留煤渣堆积状况如图1所示。
由图1可以看出，随着锥面端头移动，锥面端头前方堆积的残留煤渣（简称堆积煤渣）越来越厚。根据堆积煤渣是否超过锥面端头锥顶，锥面端头前方钻孔残留煤渣堆积状况可分为2类：堆积煤渣未超过锥顶和堆积煤渣超过锥顶。
2. 锥面端头移动最小推力理论分析
2.1 堆积煤渣未超过锥顶时的最小推力
假定锥面端头在钻孔内呈水平状态，以锥面端头锥底最低点O为原点，锥面端头移动方向（锥面端头轴线方向）为x轴正方向，垂直向上为z轴正方向，y轴正方向垂直于平面xOz，并指向右侧，建立直角坐标系，锥面端头受力分析如图2所示。
在图2中，锥面CAOB为锥面端头与煤渣的接触面（简称接触面），O1点为锥底圆心，2θ为锥面角度；C点为锥面上接触面与非接触面分界线的顶点，位于x轴正上方；h为堆积煤渣厚度。
锥面端头受力如下：封孔器推力F1整体作用于x轴正方向；重力G整体作用于z轴负方向；煤渣压力F2（挤压应力σ的总称）垂直于接触面并指向锥面端头内部；煤渣摩擦阻力F3与接触面相切并垂直于z轴（锥面端头以近水平形式移动，并不断破坏堆积煤渣的稳定性，与滑块在固定斜面的移动有着本质区别）。
为便于受力分析，需对煤渣压力与煤渣摩擦阻力沿x轴、y轴、z轴3个方向进行分解，并求出各分力计算方法。
y轴、z轴方向分力微元挤压应力分析如图3所示。由于图2中锥面CAO和锥面CBO沿平面CDO对称分布，因此，取锥面CBO进行分析即可。先作2个平行于yOz平面的截面（间距为dx），再从圆锥的中心轴线作2个截面（夹角为dβ），在锥面CBO上取微元（面积为ds），并进行受力分析。
![]()
图 3 y轴、z轴方向分力微元挤压应力分析Figure 3. Extrusion stress analysis of differential element for y-axis and z-axis component of force根据图3（a），微元面积可表示为：
$$ \mathrm{d}s=\frac{\mathrm{d}x}{\mathrm{c}\mathrm{o}\mathrm{s}\;\varphi }（H-x）\frac{r}{H}\mathrm{d}\beta $$ (1) 式中：s为锥面CBO的面积，m2；x为x轴坐标，m；φ为沿x轴方向的锥面角度，（°）；H为圆锥高度，m；r为锥面端头锥底半径，m；β为夹角，（°）。
根据图3（b），微元承受的挤压应力可分解为平行于x轴的分应力dF2x和沿着截面指向圆心的分应力dF2yz，dF2yz可表示为：
$$ {\mathrm{d}}{F}_{2{yz}}=\sigma （H-x）\frac{r}{H}{\mathrm{d}}x{\mathrm{d}}\beta $$ (2) 式中：F2yz为锥面CBO的煤渣压力沿着截面指向圆心分力，N；σ为接触面挤压应力，Pa。
对F2yz进行分解，力分解示意图如图4所示。
根据图4，F2yz平行于y轴的分力F2y可表示为：
$$ F_{2{{{y}}}}=\int_0^{\frac{h}{\mathrm{t}\mathrm{a}\mathrm{n}\; \theta}}\; \sigma（H-x）\frac{r}{H}（1-\mathrm{cos}\; \alpha）\mathrm{d}x $$ (3) 式中：F2y为F2yz平行于y轴的分力（锥面CBO的煤渣压力沿y轴方向分力），N；h为堆积煤渣厚度，m；θ为锥面端头锥面角度的1/2值，（°）；α为弧线对应夹角，（°）。
根据三角函数关系，有：
$$ （H-x）\frac{r}{H}（1-\mathrm{c}\mathrm{o}\mathrm{s}\;\alpha ）=h-x\mathrm{t}\mathrm{a}\mathrm{n}\;\theta $$ (4) 将式（4）代入式（3）并积分，有：
$$ {F}_{2{y}}=\sigma \frac{{h}^{2}}{2\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }=\sigma {S}_{\vartriangle CDO} $$ (5) 式中：${S}_{\vartriangle CDO} $为$\vartriangle $CDO的面积，m2。
$\vartriangle $CDO为锥面CBO沿y轴的投影，因此，因此，锥面CBO的煤渣压力沿y轴方向分力等于挤压应力σ与投影面积SCDO之积。
同理，F2yz平行于z轴的分力F2z可表示为：
$$ {F}_{2{z}}={\int }_{0}^{\frac{h}{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }}\sigma （H-x）\frac{r}{H}\mathrm{s}\mathrm{i}\mathrm{n}\;\mathrm{\alpha }\mathrm{d}x $$ (6) 式中：F2z为F2yz平行于z轴的分力（锥面CBO的煤渣压力沿z轴方向分力），N。
根据三角函数关系，有：
$$ （H-x）\frac{r}{H}\mathrm{s}\mathrm{i}\mathrm{n}\;\mathrm{\alpha }=\sqrt{2Rh-{h}^{2}}（1-\frac{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }{h}x） $$ (7) 将式（7）代入式（6）并积分，有：
$$ {F}_{2{z}}=\sigma \frac{h}{2\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }\;\sqrt{2Rh-{h}^{2}}=\sigma {S}_{\vartriangle CBD} $$ (8) 式中：${S}_{\vartriangle CBD} $为三角形CBD的面积，m2。
$\vartriangle $CBD为锥面CBO沿z轴的投影，因此，锥面CBO的煤渣压力沿z轴方向分力等于挤压应力σ与投影面积SCBD之积。
x轴方向分力微元挤压应力分析图如图5所示。图3中平行于x轴的微元分应力dF2x是φ的函数，由于φ非固定值，不便求解；为此，先作2个平行于xOy平面的截面（间距为dz），再作2个平行于xOz平面的截面（间距为dy），在锥面CBO上取微元（面积为ds），并进行受力分析。
![]()
图 5 x轴方向分力微元挤压应力分析图Figure 5. Extrusion stress analysis of differential element for x-axis component of force
同理，图5中微元面积可表示为：
$$ \mathrm{d}s=\frac{\mathrm{d}y}{\mathrm{s}\mathrm{i}\mathrm{n}\;\varphi }\mathrm{d}z $$ (9) 式中：y为y轴坐标，m；z为z轴坐标，m。
图5中微元承受的挤压应力平行于x轴的分应力dF2x可表示为：
$$ \mathrm{d}{F}_{2{x}}=\sigma \mathrm{d}y\mathrm{d}z $$ (10) 式中：F2x为锥面CBO的煤渣压力沿x轴方向分力，N。
对式（10）进行积分，有：
$$ {F}_{2{x}}=\frac{\sigma }{2}（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\frac{r-h}{r}-\sqrt{2rh-{h}^{2}}（r-h））=\sigma {S }_{BDO} $$ (11) 式中：SBDO为割圆BDO的面积，m2。
割圆BDO为锥面CBO沿x轴的投影，因此，锥面CBO的煤渣压力沿x轴方向分力等于挤压应力σ与投影面积SBDO之积。
同理，锥面CBO的煤渣摩擦阻力在x轴、y轴2个方向的分力F3x、F3y等于挤压应力σ、摩擦阻力系数f和投影面积（投影平面与分力方向平行且包含z轴，投影面积分别为SCDO、SBDO）之积，可表示为：
$$ {F}_{3{x}}=\sigma f{S}_{CDO}=\mathrm{\sigma }f\frac{{h}^{2}}{2\mathrm{tan}\;\theta } $$ (12) $$ {F}_{3{y}}=\frac{\sigma }{2}f（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\frac{r-h}{r}-\sqrt{2rh-{h}^{2}}（r-h）） $$ (13) 式中：F3x为锥面CBO的煤渣摩擦阻力沿x轴方向分力，N；F3y为锥面CBO的煤渣摩擦阻力沿y轴方向分力，N；f为锥面端头与煤渣的最大静摩擦系数。
当封孔器推力F1达到一定数值（最小推力F）时，锥面端头才能突破煤渣阻力向前移动，x轴、y轴、z轴3个方向的受力均处于平衡状态。y轴方向受力呈对称分布，自然平衡；x轴、z轴方向受力平衡方程可表示为：
$$ F=2{F}_{2{x}}+2{F}_{3{x}} $$ (14) $$ G=2{F}_{2\textit{z}} $$ (15) 式中：F为锥面端头移动最小推力，N；G为锥面端头重力，N。
将式（8）、式（11）、式（12）代入式（14）和式（15），并解方程，有：
$$\begin{split} &\begin{gathered} F = G\left(\mathrm{t}\mathrm{a}\mathrm{n}\; \theta\left(\frac{r^2\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\; \frac{r-h}{r}}{h\sqrt{2rh-h^2}} - \frac{r-h}{h}\right) + f\sqrt{\frac{h}{2r-h}}\right) \\ 0 < h \leqslant r \end{gathered}\\[-16pt]& \end{split} $$ (16) 根据式（16）可知，锥面端头移动最小推力F与端头重力G、端头锥面角度的一半值θ、端头锥面与煤渣的最大静摩擦系数f属于正相关，与堆积煤渣厚度h和端头锥底半径r的相关性，无法直接判定。
假设端头重力G为50 N，端头锥面角度的一半值θ为30°，最大静摩擦系数f为0.1，通过数值计算并绘制得到的锥面端头移动最小推力F随堆积煤渣厚度h和端头锥底半径r的变化曲线如图6所示。
![]()
图 6 锥面端头移动最小推力与堆积煤渣厚度和端头锥底半径的关系曲线Figure 6. Relation curves between minimum thrust of cone end movement and height of accumulated coal dregs and cone base radius of end根据图6可知：锥面端头移动最小推力F与堆积煤渣厚度h和端头锥底半径r也属于正相关。
2.2 堆积煤渣超过锥顶时的受力分析
当堆积煤渣超过锥顶时，根据是否超过锥面端头锥底最高点，可将残留煤渣堆积情况进一步划分，本研究仅分析未超过锥面端头锥底最高点的情形。
锥面端头与煤渣发生力学作用的面积相当于整个锥面与非接触面的面积之差，并忽略锥面端头锥面上方煤渣重力。在式（16）中，令h=r，求得的F值乘以2，即为整个锥面与煤渣接触条件下的锥面端头移动最小推力；再令h=2r−h，求得的F值，即为非接触面与煤渣接触条件下的锥面端头移动最小推力；上述两者之差即为堆积煤渣超过锥顶时的锥面端头移动最小推力F，如式（17）：
$$\begin{split} &\begin{gathered} F=G\mathrm{t}\mathrm{a}\mathrm{n}\;\theta \left(\pi +\frac{f}{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta}\left(1-\sqrt{\frac{2r}{h}-1}\right)-\right.\\ \left.\left. \frac{{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\;\frac{h-r}{r}}{\left(2r-h\right)\sqrt{2rh-{h}^{2}}}+\frac{h-r}{2r-h}\right)\right. \qquad r < h \leqslant 2r \end{gathered}\\[-18pt]& \end{split} $$ (17) 3. 锥面端头移动最小推力试验研究
锥面端头移动最小推力F的影响因素较多，因篇幅受限，仅试验研究锥面端头移动最小推力随端头锥面角度的变化规律，并验证理论分析结果的可靠性。
试验方案如下：采用弹簧秤牵拉带有锥面端头的封孔器并让其在含煤渣的管道内匀速移动，弹簧秤显示读数即为锥面端头移动最小推力；锥面角度分别是60°、45°、30°、20°、18°、15°，共6组试验。试验装置设计图如图7所示，锥面端头移动最小推力与锥面角度的关系曲线如图8所示。
![]()
图 8 锥面端头移动最小推力与锥面角度的关系曲线Figure 8. Relation curves between minimum thrust of cone end movement and conical surface angle由图8可知：锥面角度增加，锥面端头移动最小推力升高，并且上升速率加快；移动距离越远，上升速率加快趋势越明显。根据式（16），锥面端头移动最小推力F与tan θ成正比，tan θ为正切函数，随θ增大而升高并且上升速度逐渐加快。由此可见，试验研究结果与理论分析结果基本一致。
4. 锥面端头减阻机理
与平面端头相比，锥面端头可以通过接触面对堆积煤渣施加平行于y轴的侧向力；因此，可通过分析侧向力来探索锥面端头减阻机理。因篇幅受限，以堆积煤渣未超过锥顶的情形为例，进行分析。
堆积煤渣承受的侧向力Fy是锥面端头y轴方向受力的反作用力，方向由接触面指向堆积煤渣，并平行于y轴，有：
$$ {F}_{{y}}={F}_{2{y}}-{F}_{3{y}} $$ (18) 将式（5）、式（13）、式（15）代入式（18），有式（19）：
$$\begin{split} &\begin{gathered} {F}_{{y}}=G\frac{h}{2\sqrt{2rh-{h}^{2}}}-\frac{G\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }{h\sqrt{2rh-{h}^{2}}}f（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\;\frac{r-h}{r}-\\ \sqrt{2rh-{h}^{2}}（r-h）） \end{gathered}\\[-19pt]& \end{split} $$ (19) 式中：Fy为堆积煤渣承受的侧向力，N。
根据式（19）可知，θ越小，Fy越大，即锥面端头越尖，Fy越大。
与锥面端头相接触的煤渣，在侧向力Fy的推动下，向钻孔孔壁侧移动（或被压缩）；导致锥面端头通过后，锥面端头后方残留煤渣增加，锥面端头前方堆积煤渣减少，堆积煤渣厚度h降低；锥面端头越尖，侧向力Fy越大，堆积煤渣厚度h越低。
锥面端头移动的突破点通常位于接触面顶端（堆积煤渣未超过锥顶时，在C点；堆积煤渣超过锥顶，在锥面端头锥顶），如果接触面顶端承受的应力超过该点的堆积煤渣剪切应力，则堆积煤渣遭到破坏、锥面端头可向前移动。
由于堆积煤渣相邻颗粒之间具有一定的咬合作用，图9（b）的某片煤渣薄板受侧向力挤压后，受压效应向前传递（距离越远，效应越弱），导致该煤渣薄板至钻孔孔壁一侧应力升高，而至钻孔中心一侧应力下降（相当于1根受压弹簧，从某个位置施加新力后，受压一侧弹簧绷紧程度加强，而另一侧弹簧绷紧程度减弱）；虽然煤渣属于散体，单片煤渣薄板微观上的受压效应传递能力较低，但是，当众多煤渣薄板的微观受压效应传递累积起来时，宏观上可以在一定程度上影响应力分布。
在全体煤渣薄板受压效应传递综合作用下，从C点到钻孔孔壁（E1、E2）的堆积煤渣应力变化如下：靠近C点一侧应力下降，靠近钻孔孔壁（E1、E2）一侧应力上升，C点（位于钻孔中心，并且距侧向力最近）应力下降大。根据散体动力学理论[17]，摩尔库伦破坏准则适用于散体，堆积煤渣剪切应力与正应力成正比。因此，受侧向力挤压效应传递的影响，堆积煤渣C点及其邻近区域的剪切应力降低，C点下降值最大。
综上所述，锥面端头减阻机理如下：锥面端头越尖（θ越小），侧向力Fy越大，堆积煤渣厚度h越低、接触面顶端及其临近区域堆积煤渣剪切应力越小，堆积煤渣对锥面端头移动的阻碍能力越弱。
注：C点及其邻近区域堆积煤渣受力破坏后，堆积煤渣整体失去稳定性，靠近钻孔孔壁（E1、E2）一侧堆积煤渣应力自然得到释放。
5. 封孔器端头优化改进与现场试验
根据锥面端头移动最小推力与影响因素的关系，可从以下几个方面对锥面端头开展优化设计：①采用低密度材料制作锥面端头，降低端头重力；②适当减小端头锥面角度；③采用喷漆等措施提高端头锥面的光滑性，降低端头锥面与煤渣的最大静摩擦系数；④采用更细的封孔器，从而减小端头锥底半径。
在文献[15]中，将1个空心外罩（锥面角度约20°）安设于钻孔瓦斯涌出初速度测定装置的封孔器端头，端头由台阶式转变为锥面式；在此基础上，根据锥面端头优化设计思路，对封孔器端头作进一步改进。
根据图8，当锥面角度较低时，降低锥面角度改变锥面端头移动最小推力的作用不大；由于端头锥底半径与封孔器主体部分一致，而主体部分需要考虑膨胀性能、耐压性能和膨胀液输送、内置管道气体流动阻力等因素，截面尺寸（目前为38 mm）很难进一步缩小。因此，从如下2个方面对封孔器端头进行改进：①将空心外罩的制作材料由不锈钢改合金，减轻端头质量；②对封孔器锥面端头打蜡，降低端头锥面与煤渣的最大静摩擦系数。
在开滦矿区钱家营矿1355工作面进行了钻孔瓦斯涌出初速度测定现场试验，现场试验数据见表1。
表 1 现场试验数据Table 1. Experimental data in field test封孔器型号 最高送入孔深/m 孔深/m 时间/min 原始型 16 8 1.2 10 1.3 12 1.5 14 1.7 16 2.0 改进型 >20 8 0.9 10 1.0 12 1.2 14 1.3 16 1.5 18 1.7 20 1.9 由表1可知：改进型封孔器送入阻力小、推进速度快、送入钻孔深度大。
6. 结 语
1）锥面端头与煤渣接触面所承受的煤渣压力沿x轴、y轴、z轴的分力均等于挤压应力与投影面积之积；锥面端头与煤渣接触面所承受的摩擦阻力沿x轴、y轴的分力也均等于挤压应力、摩擦阻力系数和投影面积之积。
2）从理论上推导了锥面端头移动最小推力计算公式，它与端头质量、端头锥面角度、端头锥面与煤渣的最大静摩擦系数、堆积煤渣厚度和端头锥底半径均具有正相关性，为锥面端头优化设计指明了方向。
3）锥面端头移动最小推力与锥面角度呈正切函数关系，随锥面角度增大而升高并且上升速度逐渐加快。
4）通过侧向力分析，揭示了锥面端头减阻机理；锥面端头越尖，侧向力越大，堆积煤渣厚度越低、接触面顶端及其临近区域堆积煤渣剪切应力越小，堆积煤渣对锥面端头移动的阻碍能力越弱。
5）根据锥面端头优化设计思路改进后的封孔器端头，可以降低送入阻力、提高送入速度与孔深。
-
![]()
图 2 不同初始状态下裂隙闭合处水锤压力波动幅值
Figure 2. Amplitudes of water hammer pressure fluctuation at fracture closure under different initial conditions
![]()
图 3 不同初始状态下裂隙中点处水锤压力波动幅值
Figure 3. Amplitude of water hammer pressure fluctuation at the fracture midpoint under different initial conditions
![]()
图 4 裂隙内初始水压与最大瞬时压力对比
Figure 4. Initial water pressure in the fracture versus maximum instantaneous pressure
![]()
图 5 水锤发生不同瞬时裂隙内水锤压力波传播规律
Figure 5. Water hammer pressure wave propagation laws within the fissure at different transients of water hammer occurrence
表 1 模型参数表
Table 1 Model parameters table
参数名称 方案1 方案2 方案3 方案4 煤体弹性模量/GPa 15 15 15 15 裂隙长度/m 10 10 10 10 裂隙宽度/m 0.1 0.1 0.1 0.1 初始压力/MPa 1.0 2.0 3.0 4.0 初始流量/（mL·min−1） 22.0 36.0 42.0 55.0 温度/K 293.15 293.15 293.15 293.15 流体密度/（kg·m−3） 1 000 1 000 1 000 1 000 泊松比 0.33 0.33 0.33 0.33 流体体积模量/GPa 2.10 2.10 2.10 2.10 
下载: 导出CSV
-
[1] 柴虎强. 高压脉动水力锤击对煤体介质损伤致裂特性的研究[D]. 徐州:中国矿业大学,2015. [2] 李全贵,林柏泉,翟成,等. 煤层脉动水力压裂中脉动参量作用特性的实验研究[J]. 煤炭学报,2013,38(7):1185−1190. LI Quangui, LIN Baiquan, ZHAI Cheng, et al. Experimental study on action characteristic of pulsating parameters in coal seam pulse hydraulic fracturing[J]. Journal of China Coal Society, 2013, 38(7): 1185−1190.
[3] 刘沛清,李忠义,冬俊瑞. 用二维瞬变流方程分析缝面层中脉动压力传播规律[J]. 水利学报,1996,27(4):27−33. LIU Peiqing, LI Zhongyi, DONG Junrui. On the mechanism of fluctuating pressure transmission in cracks of the bed rock[J]. Journal of Hydraulic Engineering, 1996, 27(4): 27−33.
[4] 李爱华,刘沛清. 脉动压力在板块缝隙中传播衰变机理研究[J]. 水利水电技术,2006,37(9):33−37. LI Aihua, LIU Peiqing. Study on mechanism of pulsating pressure propagation disintegration in joints of bottom slab[J]. Water Resources and Hydropower Engineering, 2006, 37(9): 33−37.
[5] 张建民,杨永全,戴光清,等. 水垫塘底板缝隙中脉动压力传播特性[J]. 四川大学学报(工程科学版),2000,32(3):5−8. ZHANG Jianmin, YANG Yongquan, DAI Guangqing, et al. Study on fluctuating pressure transmission in joint on the bottom of plunge pool[J]. Journal of Sichuan University (Engineering Science Edition), 2000, 32(3): 5−8.
[6] 贺培,熊继有,陆朝晖,等. 基于高压脉动水力压裂的脉动波传播机理研究[J]. 水动力学研究与进展(A辑),2017,32(1):117−123. HE Pei, XIONG Jiyou, LU Zhaohui, et al. Study of pulse wave transmission mechanism based on high pulse hydraulic fracturing[J]. Chinese Journal of Hydrodynamics, 2017, 32(1): 117−123.
[7] 赵耀南,梁兴蓉. 水流脉动压力沿缝隙的传播规律[J]. 天津大学学报,1988,21(3):55−65. ZHAO Yaonan, LIANG Xingrong. Investigation on fluctuating pressure transmission along joints[J]. Journal of Tianjin University, 1988, 21(3): 55−65.
[8] 姜文超,梁兴蓉. 应用紊流理论探讨脉动压力沿缝隙的传播规律[J]. 水利学报,1983,14(9):53−59. [9] 李爱华,朱江. 岩石裂隙中水流脉动压力传播的数值模拟[J]. 水利水电科技进展,2015,35(4):24−27. LI Aihua, ZHU Jiang. Numerical simulation on propagation of fluctuating pressure of flow within rock fractures[J]. Advances in Science and Technology of Water Resources, 2015, 35(4): 24−27.
[10] 满春雷,张景松,张飞. 煤层高压脉动水锤器的脉冲特性研究[J]. 矿山机械,2012,40(5):20−23. MAN Chunlei, ZHANG Jingsong, ZHANG Fei. Study on pulsation characteristics of high pressure pulsation water hammer for coal seams[J]. Mining & Processing Equipment, 2012, 40(5): 20−23.
[11] 杨春敏,张景松. 水锤及水锤压力波在水系中传播的分析与试验[J]. 流体机械,2014,42(9):10−13. YANG Chunmin, ZHANG Jingsong. Apparatus inducing water-hammer and the analysis and experiment of the water-hammer pressure wave propagation in the water system[J]. Fluid Machinery, 2014, 42(9): 10−13.
[12] 庄利,张景松,满春雷,等. 脉动水锤增透技术对煤体有效孔隙率影响实验研究[J]. 煤矿安全,2014,45(6):18−21. ZHUANG Li, ZHANG Jingsong, MAN Chunlei, et al. Experimental research of the impact of enhanced permeability technology by pulsation water-hammer on effective porosity of coal body[J]. Safety in Coal Mines, 2014, 45(6): 18−21.
[13] 胡斌,张景松,高飞. 高压脉动水锤装置压力冲击与柱塞结构的关系[J]. 煤矿安全,2013,44(7):220−222. HU Bin, ZHANG Jingsong, GAO Fei. The relation between pressure impact and plunger structure of high pressure pulsating water hammer device[J]. Safety in Coal Mines, 2013, 44(7): 220−222.
[14] FIOROTTO V, RINALDO A. Fluctuating uplift and lining design in spillway stilling basins[J]. Journal of Hydraulic Engineering, 1992, 118(4): 578−596. doi: 10.1061/(ASCE)0733-9429(1992)118:4(578)
[15] 张志军. 长距离供水管道水锤分析优化与应用研究[D]. 西安:西安建筑科技大学,2016. [16] E BENJAMIN Wylie, VICTOR L Streeter. Fluid transients[M]. New York: Mc Graw-Hill International Book Co, 1978.
[17] 余旭. 煤层脉动水力压裂压力波沿裂隙传播规律和解堵机理研究[D]. 徐州:中国矿业大学,2016.
计量
- 文章访问数: 56
- HTML全文浏览量: 3
- PDF下载量: 11
