Study on aggregate accumulation and growth mechanism in underground dynamic water cutting-off construction
-
摘要:
灌注骨料形成阻水屏障实现截流降速是井下治水工程的关键所在。为定量描述骨料堆积生长的时空演化机理,提出了将骨料沉降堆积过程切分为具有一定厚度的片状颗粒层堆叠过程的分析计算方法,以及动水环境中骨料灌注后从沉降到堆积全过程关键参量(颗粒水平移距xd,动水休止角ψ,未接顶区流速U,骨料留存临界流速Ucr)的理论计算公式,在此基础上构建了骨料堆积体生长预测模型。通过CFD-DEM双向耦合算法研究了骨料灌注期间流场分布演化特征和骨料堆积形态差异,并验证了预测模型合理性。研究表明:空间上,依据初始流速Us是否超过骨料起动流速Uc或骨料留存临界流速Ucr,可将截流初期骨料沉积主域的所在位置划分为3个区域,若沉积主域位于③区,则将其视为无效灌注;时间上,骨料堆积体生长过程可归结为3个阶段:高度快速增长阶段、高度长度同步生长阶段、仅水平向伸长阶段;初始流速决定了第1阶段的有无,而灌注条件(骨料粒径、灌注速度)主导着第1阶段所持续的时间;基于泥沙运动力学中推移质输沙率、起动流速等概念,得出了堵孔现象发生及骨料能否留存的临界判据;数值模拟试验结果表明,预测模型能够较好地刻画骨料灌注后的堆积生长规律,骨料堆积稳定后,理论计算值与模拟试验值相对误差小于10%。
Abstract:The key to inrush water control engineering in coal mine lies in forming a water-resistant barrier through the pouring of aggregate, which achieves flow interception and reduction. In order to quantitatively describe the spatiotemporal evolution mechanism of aggregate accumulation, an analytical calculation method is proposed to divide the process of aggregate settling and stacking into the process of particle layers superposition with a certain thickness and theoretical formulas for the entire process from settlement to pile after aggregate pouring in a dynamic water environment about critical parameters (horizontal displacement of particle xd, dynamic repose angle ψ, flow velocity in the topping zone U, critical flow velocity for aggregate retention Ucr) are proposed. Based on this, a prediction model about the growth of aggregate accumulation is constructed. The distribution characteristics of the flow field and the differences in aggregate accumulation morphology during the pouring period are studied with CFD-DEM, and the rationality of the prediction model is verified. The study reveals that spatially, the location of the main sedimentary domain during the initial interception phase can be divided into three regions, depending on whether the initial flow velocity Us exceeds the incipient flow velocity for particles Uc or the critical flow velocity for aggregate retention Ucr. That is considered ineffective pouring if the main sedimentary domain lies in region ③. Temporally, the growth process of aggregate accumulation can be summarized into three stages: a height rapid increase stage, a synchronous growth stage in height and length, and a stage with only horizontal elongation. The presence of the first stage is determined by the initial flow velocity, while the duration of the first stage is predominantly governed by the pouring conditions (particle size, pouring rate). A critical criterion for pore clogging and aggregate retention is derived based on sediment transport rate and incipient flow velocity from sediment dynamics. The numerical simulation experimental results indicate that the prediction model effectively characterizes the growth law of aggregate accumulation because the relative error between the theoretical calculation values and the simulated experiment values is less than 10% after the aggregate stacking stabilizes.
-
0. 引 言
煤矿巷道发生突水淹井事故后,通常先灌注骨料形成高阻弱渗堆积体,将管道流变为渗流,为后期注浆工作中水泥浆液的有效积聚和凝固创造条件[1]。
其中,骨料堆积体的建成是决定封堵成败的关键,众多学者在此方面进行了大量的理论研究和模拟试验工作。邵红旗等[2]提出在静水条件下骨料灌注效率低时采用双液浆法快速建造阻水墙骨架,再采用综合注浆法灌注水泥浆液加固墙体,达到了快速堵水抢险救援的目的;牟林[3]从骨料灌注期间的水力学状态出发,对截流工程技术综合优化等方面进行了研究,将截流过程分为了4个阶段,即骨料灌注前的施工准备阶段、铺底–充填阶段、骨料接顶阶段、注浆加固阶段;董书宁等[4]基于流体力学理论,考虑水头压差、渗流突破和浆液初凝时间等要素,建立了截流所需骨料堆积段最小临界长度估算模型;杨志斌等[5]针对突水巷道截流接顶困难这一技术难点,提出了截流、堵源同步进行的治理方案。为研究骨料灌注后在动水中的堆积过程,惠爽[6]设计研制了可视化的矿井淹没巷道动水骨料灌注试验平台,对骨料截流过程中的影响因素及堵水效果进行了分析;李军[7]采用FLUENT软件开展骨料灌注阶段的数值模拟试验,研究了骨料颗粒的运动轨迹、运动速度及管道内水压力变化情况。
以往研究成果以定性研究截流工作中骨料灌注后其运移堆积过程及封堵效果评价为主,未能提出能够描述骨料堆积生长时空演化过程的定量分析方法,导致其内在力学机理仍缺乏深入研究。鉴于此,通过理论分析满流状态下骨料灌注后的沉积主域位置与动水流速间的响应关系,建立骨料可留存堆积的临界判据;采用以单位时间灌注量对骨料颗粒分层的处理手段,将连续的骨料灌注堆积过程简化为颗粒层堆叠过程,并构建了能够反映随骨料持续灌注,堆积体高度及长度增长情况的骨料堆积生长预测模型,以期为煤矿水害治理工程施工提供参考。
1. 堆积形态时空演化规律
1.1 基本假定
1) 采用泥沙运动力学中输沙率(kg/s)表示水流携砂能力qs* (kg/s),且仅考虑水流携带骨料颗粒中推移质[8]的能力;
2) 堆积体形成稳定结构后,其内部骨料堆积密实,阻水能力强,同一横截面不同高度处流速相同;
3) 实测数据[9]表明,截流成功前的流量基本稳定,水流绝大部分从巷道顶部未接顶区域通过,堆积体接顶成功后流量快速下降。
1.2 灌注初期沉积主域空间分布
对灌注初期沉积主域的研究,能够帮助预判骨料灌注后能否顺利堆积形成具有一定规模的阻水体。事实上,粒径较大的均匀颗粒群在沉降时由于颗粒间的碰撞挤压,会出现明显分层凝聚[10]。据此将连续的骨料沉降及堆积过程拆解简化为无限个具有相同厚度δ的片状颗粒层堆叠及塌落过程,如图1所示。
取首个颗粒层为隔离体进行分析,如图2所示,轨迹线1、2、3为颗粒层在不同流场中的沉降路径。静水中(Us=0),颗粒沿轨迹线1沉降,后一颗粒层堆叠在前层之上,两颗粒层中心线重合,床面上部分颗粒(红色)沿水下休止角θ塌落(图3a);当初始流速Us大于零时,颗粒层沉降在轨迹线1右侧(①、②、③区),沉积主域与孔口中心水平距离记为xd。由于颗粒层铺底后主流区被压缩,床面瞬时水流速升高,导致后一颗粒层着床前水平位移略微增大,与首层在垂向上有微小 “错动”(图3b),两颗粒层水平移距差记为Δxd;轨迹线2为初始流速Us与颗粒临界起动流速Uc相同时颗粒层的沉降路径,将该情形下的xd记为xdc,若Us超过Uc,水流携砂能力qs不再为0,颗粒层在②区沉降后床面层部分骨料(红色颗粒)被二次起动带离至下游(图3c);轨迹线3为初始流速Us与颗粒层能够有效堆叠的临界流速Ucr(骨料留存临界流速)相同时的颗粒层沉降路径,若Us超过Ucr,则骨料在③区沉降触底后无法在原位置附近稳定沉积,而被高速水流加速搬运至下游(红色颗粒),导致后一层颗粒无法着床(图3d),该情形为无效灌注。当灌注速度Φ小于水流携砂能力qs时,骨料无法留存,因此本文将qs与Φ相同时的水流速定义为骨料留存临界流速Ucr,水流携砂能力为单位时间内被水流搬运带离至下游的骨料质量。
1.3 骨料堆积生长影响因素分析及预测模型建立
从分析静水条件下影响骨料堆积的关键因素出发,通过增加水流对骨料沉降堆积的作用,逐步过渡到动水中骨料堆积生长机制的研究中,并建立动水中骨料灌注后的堆积体生长预测模型。
将进入巷道的骨料按单位时间切分,每单位时间内进入巷道的骨料为同一颗粒层,首个颗粒层在巷底沉积的瞬间记为0时刻,首层颗粒沉积后,其高度ys和长度ls可通过灌注速度、粒径和孔隙率三者的关系得到:
$$ {y_{\text{s}}} = \frac{{{V_{\text{P}}}}}{{{A_{\text{H}}}}} = \frac{{{V_{\text{V}}}}}{{(1 - \varepsilon ){A_{\text{H}}}}} = \frac{\phi }{{(1 - \varepsilon ){\rho _{\text{s}}}{A_{\text{H}}}}} $$ (1) $$ {l_{\text{s}}} = \frac{{{V_{\text{P}}}}}{{{y_{\text{s}}}L}} = \frac{{{V_{\text{V}}}}}{{(1 - \varepsilon ){y_{\text{s}}}L}} = \frac{\phi }{{(1 - \varepsilon ){\rho _{\text{s}}}{y_{\text{s}}}L}} $$ (2) 式中:VP为颗粒占据空间体积,m3;VV为颗粒总体积,m3;ϕ为灌注速度,kg/s;L为巷道宽度,m;ε为孔隙率,%;AH为投料孔面积,m2;ρs为颗粒密度,kg/m3。
后续的骨料颗粒着床后,一部分稳定沉积在床面以增加堆积高度,另一部分在重力作用下沿休止角θ塌落(图3a),则1时刻的堆积高度增量$\Delta {y_{(1)}} $为
$$ \Delta {y_{(1)}} = \frac{{m_{\text{r}}}_{(1)}}{{(1 - \varepsilon ){\rho _{\text{s}}}L({l_{\text{s}}} - 2{y_{\text{s}}}{\text{/}}\tan\; \theta )}} $$ (3) 其中:
$$ {m_{\text{r}}}_{(1)} = \phi - m{'_{(1)}} $$ (4) $$ m{'_{(1)}} = \dfrac{{{l_{\text{H}}} - \left({l_{\mathrm{s}}} - \dfrac{{4{y_{\mathrm{s}}}}}{{\tan\; \theta }}\right)}}{{{l_{\text{H}}}}}\phi $$ (5) 式中:mr(t)为t时刻稳定着床的颗粒质量,kg;m′(t)为t时刻沿休止角塌落的颗粒质量,kg;θ为水下休止角,(°);lH为投料孔边长,m。
1时刻堆积高度为:
$$ {y_1} = {y_{\text{s}}} + \Delta {y_{(1)}} $$ (6) 骨料颗粒层沉降着床后,部分颗粒塌落并附着在堆积体两侧坡面使休止角大小保持不变(图4a、图4b)。为便于描述并计算这一过程的堆积高度及长度增量,可将这部分颗粒视为“颗粒桩”竖直嵌入堆积体内部,使堆积体向两侧伸展,如图4c所示,图中红色颗粒质量与图4b相同。
则堆积体长度增量计算公式可表示为:
$$ \Delta {l_{(1)}} = \frac{{m{'_{(1)}}}}{{(1 - \varepsilon ){y_1}{\rho _{\text{s}}}L}} $$ (7) 第二层颗粒沉降后的堆积体长度为:
$$ {l_1} = {l_{\text{s}}} + \Delta {l_{(1)}} $$ (8) t时刻堆积体高度增量为:
$$ \Delta {y_{(t)}} = \frac{{{m_{\text{r}}}_{(t)}}}{{(1 - \varepsilon ){\rho _{\text{s}}}L({l_{t - 1}} - 2{y_{t - 1}}{\text{/}}\tan\; \theta )}} $$ (9) t时刻堆积体长度增量为:
$$ \Delta {l_{(t)}} = \frac{{m{'_{(t)}}}}{{(1 - \varepsilon ){y_t}{\rho _{\text{s}}}L}} $$ (10) 其中:
$$ {m_r}_{(t)} = \phi - m{'_{(t)}} $$ (11) $$ m{'_{(t)}} = \frac{{{l_{\text{H}}} - \left[ {{l_{t - 1}} - \frac{2}{{\tan\; \theta }}({y_{t - 1}} + {y_s})} \right]}}{{{l_{\text{H}}}}}\phi $$ (12) t时刻的堆积高度和长度分别为:
$$ {y_t} = {y_{t - 1}} + \Delta {y_{(t)}} = {y_{\text{s}}} + \sum\limits_{x = 1}^t {\Delta {y_{(x)}}} $$ (13) $$ {l_t} = {l_{t - 1}} + \Delta {l_{(t)}} = {l_{\text{s}}} + \sum\limits_{x = 1}^t {\Delta {l_{(x)}}} $$ (14) 动水环境中,若初始流速较小,则可以认为首层高度ys及长度ls与静水中相同。考虑到因骨料铺底后主流区被压缩,流速升高使后续颗粒层在沉降时产生较大水平位移从而使部分颗粒无法着床,将这部分颗粒质量记为m′′。另外,由于迎水坡处塌落量较少(图3b),因此在计算时将其忽略以简化计算,得到高度增量计算公式:
$$ \Delta {y_{(t)}} = \frac{{{m_r}{{_{(t}}_)}}}{{(1 - \varepsilon ){\rho _{\text{s}}}L({l_{t - 1}} - {y_{t - 1}}{\text{/}}\tan\; \psi )}} $$ (15) 其中:
$$ {m_{\text{r}}}_{(t)} = \phi - m{'_{(t)}} - m'{'_{(t)}} $$ (16) $$ m'{'_{(t)}} = \frac{{\Delta {x_{\text{d}}}_{{\text{(}}t{\text{)}}}}}{{{l_{\text{H}}}}}\phi = \left[\frac{{{u_{(t)}}(H - {y_{t - 1}})}}{{{u_{(t - 1)}}(H - {y_{t - 2}})}} - 1\right]\frac{{{x_{\text{d}}}_{{\text{(}}t - 1{\text{)}}}{y_{t - 1}}}}{{{y_{\text{s}}}{l_{\text{H}}}}}\phi $$ (17) 式中:xd(t)为t时刻的颗粒沉底时的水平移距,m;H为巷道高度,m。
由于水流速未达到颗粒起动流速,水下休止角保持不变,则m′计算公式如下:
$$ m{'_{(t)}} = \frac{{{l_{\text{H}}} - \left[ {{l_{t - 1}} - \dfrac{2}{{\tan\; \theta }}({y_{t - 1}} + {y_{\mathrm{s}}})} \right]}}{{{l_{\text{H}}}}}(\phi - m'{'_{(t)}}) $$ (18) t时刻堆积体长度增量为:
$$ \Delta {l_{(t)}} = \frac{{m{'_{(t)}} + m'{'_{(t)}}}}{{(1 - \varepsilon ){y_t}{\rho _{\text{s}}}L}} $$ (19) 式(15)—式(19)成立前提为随堆积高度增加,顶区流速U始终未达到Uc。随着堆积高度增加,水流速进一步抬升,当U ≥ Uc时,堆积体的生长进入第2阶段;若堆积高度达到H-D时流速仍小于Uc,则将会引发堵孔,导致有效堆积段过短,即:
$$ \left\{ \begin{gathered} {y_t} \geqslant H - D \\ U \leqslant {U_{\text{c}}} \\ \end{gathered} \right. $$ (20) 式中:D为颗粒粒径,m。
当堆积体生长进入第2阶段时,由水流速度所主导的携砂能力qs > 0,因此在计算高度和长度增量时,需要将其考虑在内,高度增量公式同式(15)。
其中:
$$ {m_r}_{(t)} = \phi - m{'_{(t)}} - m'{'_{(t)}} - {q_{{\text{s}}(t)}} $$ (21) 由于初始沉积主域右侧长度远超过左侧长度,因此在计算长度增量时省略左侧部分,计算公式如下:
$$ \Delta {l_{(t)}} = \frac{{m{'_{(t)}} + m'{'_{(t)}} + {q_{{\text{s}}(t)}}}}{{(1 - \varepsilon ){y_t}{\rho _{\text{s}}}L}} $$ (22) 如图5所示,当床面长度超过灌注孔边长时,即式(23)成立,则可以消除颗粒层塌落和偏移所导致的增量变化。
$$ {l}_{t}-\frac{{y}_{t}}{\mathrm{tan}\;\psi }>{l}_{\text{H}} $$ (23) 高度和长度增量计算公式变为:
$$ \Delta {y_{(t)}} = \frac{{\phi - {q_{{\text{s}}(t)}}}}{{(1 - \varepsilon ){\rho _{\text{s}}}L({l_{t - 1}} - {y_{t - 1}}{\text{/}}\tan\; \psi )}} $$ (24) $$ \Delta {l_{(t)}} = \frac{{{q_{{\text{s}}(t)}}}}{{(1 - \varepsilon ){\rho _{\text{s}}}L{y_t}}} $$ (25) t时刻堆积高度和长度计算公式同式(13)、式(14)。当携砂能力qs*不断增大至与灌注速度ϕ相同时,即顶区流速达到Ucr,骨料灌注后无法在堆积体顶部留存,则预示着生长进入第3阶段,Δy减小为0,此阶段堆积体仅水平向长度增加。
若Us超过Uc,灌注初期沉积主域位于②区,则堆积体生长直接从第2阶段开始,高度Δy、长度Δl增量计算同式(15)、式(19),首层高度ys、长度ls为:
$$ {y_{\text{s}}} = \frac{{{V_{\text{P}}}}}{{{A_{\text{H}}}}} = \frac{{{V_{\text{V}}}}}{{(1 - \varepsilon ){A_{\text{H}}}}} = \frac{{\phi - {q_{\text{s}}}}}{{(1 - \varepsilon ){\rho _{\text{s}}}{A_{\text{H}}}}} $$ (26) $$ {l_{\text{s}}} = \frac{{{V_{\text{P}}}}}{{{y_{\text{s}}}L}} = \frac{{{V_{\text{V}}}}}{{(1 - \varepsilon ){y_{\text{s}}}L}} = \frac{{\phi + {q_{\text{s}}}}}{{(1 - \varepsilon ){\rho _{\text{s}}}{y_{\text{s}}}L}} $$ (27) 工程现场流速一般小于0.5 m/s,小于截流工程常用粒径(5~30 mm)骨料的起动流速[3],因此堆积生长过程多从第1阶段开始,若灌注方案制定不当,则极易造成堵孔,因此截流初期选取小粒径骨料灌注,当堆积体具有一定高度及长度后,再行更换大粒径骨料或增大灌注速度直至堆积体接顶。
2. 关键参量的确定
一个适用性强、准确度高的骨料堆积体生长预测模型,不仅取决于正确的理论基础、科学的假定处理及合理的公式结构,还依赖于所取参量是否合理。对于本文预测模型,主要涉及颗粒动水沉降时水平移距xd,散体颗粒堆积动水休止角ψ,巷道顶区未接顶空间流速U,以及水流携砂能力qs等关键参量的确定。
2.1 颗粒动水沉降时水平移距计算公式
基于恒定流条件计算颗粒在满流巷道中自由沉降时的水平移距。假设投料孔较高,孔内充满水,颗粒由投料孔进入巷道前已经达到颗粒沉降末速ug (m/s)。由于颗粒水平向作变加速运动,水流拖曳力FD相比于其他作用力具有主导作用[11],因此影响作用较小的力不予考虑。另外,需要考虑因颗粒作变加速运动而产生的与运动方向相反的惯性力FG,如图6所示。
作用于颗粒上的水流拖曳力:
$$ {F_{\text{D}}} = {C_{\text{D}}}\frac{{\pi {D^2}}}{4}\frac{{\rho {{(U - u)}^2}}}{2} $$ (28) 式中:ρ为流体密度,kg/m3;u为颗粒水平向移速,m/s;CD为拖曳力系数,其经验公式为[12]:
$$ {C}_{\text{D}}=\left\{\begin{array}{l}\dfrac{24}{Re}(1+0.15R{e}^{0.687}),Re<\text{1 000}\\ 0.44,\text{ }Re\geqslant \text{1 000}\end{array}\right. $$ (29) 式中:Re为雷诺数。
颗粒变加速运动产生的惯性力[13]:
$$ {F_G} = (M' + \lambda )\frac{{{\mathrm{d}}u}}{{{\mathrm{d}}t}} $$ (30) 式中:M'为颗粒质量,kg;λ为颗粒因不定常运动带来的虚质量[14],kg;分别表示为
$$ M' = \frac{{\pi {D^3}{\rho _{\text{s}}}}}{6} $$ (31) $$ \lambda = \frac{{\pi {D^3}\rho }}{{12}} $$ (32) 根据牛顿定律得到颗粒运动平衡方程:
$$ {F_{\text{D}}} = {F_G} $$ (33) 将式(28)、式(30)—(32)代入式(33)得到:
$$ k{(U - u)^2} = \frac{{{\mathrm{d}}u}}{{{\mathrm{d}}t}} $$ (34) 其中:
$$ k = \frac{{3{C_{\text{D}}}}}{{4d\left(\dfrac{{{\rho _{_{\mathrm{S}}}}}}{\rho } + \dfrac{1}{2}\right)}} $$ (35) 对式(34)积分,并代入边界条件 t = 0,u = 0;t = 0,x = 0,可得:
$$ x = \int_0^t {u{\mathrm{d}}t} = Ut - \frac{{\ln (kUt + 1)}}{k} $$ (36) 将t = H/ug代入式(36),并改用颗粒群沉降末速代替单颗粒沉降末速,得到颗粒层中心处颗粒沉底时的水平移距:
$$ \begin{array}{c} {x_{\text{d}}} = \dfrac{{UH}}{{{u_{\mathrm{g}}}}} - \\ \dfrac{{4\ln \{3{C_{\text{D}}}UH{\text{/}}\{[(4D{\rho _{\text{s}}}{\text{/}}\rho + 2D) + 1]{u_{\text{g}}}\}\}({\rho _{\text{s}}}{\text{/}}\rho + 1{\text{/}}2)}}{{3{C_{\text{D}}}{\text{/}}D}} \end{array}$$ (37) 颗粒沉底时的水平移距xd决定了骨料灌注后沉积主域的所在位置,若xd大于xdc,主域落于②区,则预示着堆积体生长不存在第1阶段(堆积高度迅速增长阶段),而直接从第2阶段(高度长度同步生长阶段)开始。
2.2 散体颗粒堆积动水休止角计算公式
动水中骨料散体堆积休止角ψ是水流和堆积体相互作用塑造的一种动态平衡形态,坡面上的颗粒处于临界起动状态,动水休止角是描述这一临界状态的关键物理参数[15]。本节从堆积体稳定坡面上颗粒受力分析出发,引入作用流速Ub建立动水休止角计算公式。
如图7所示,颗粒主要受到有效重力G′,水流拖曳力FD、上举力FL的作用。其中,拖曳力FD计算公式同式(28),式中U = Ub;u = 0。
$$ {F_{\text{L}}} = {C_{\text{L}}}\frac{{\pi {D^2}}}{4}\frac{{\rho U_{\mathrm{b}}^2}}{2} $$ (38) $$ G' = \frac{{\pi {D^3}\left( {{\gamma _{\text{s}}} - \gamma } \right)}}{6} $$ (39) 式中:FL为上举力,kN; CL为上举力系数,对暴露于床面的单球体颗粒一般取0.1[16];Ub为作用流速,m/s;G′为有效重力,kN;γs为颗粒容重,kN/m3;γ为流体容重,kN/m3。
突水稳定阶段,巷道内水流流动为管道流;骨料灌注阶段,堆积尚未接顶时,巷道内部具有床沙的二元结构,未接顶区流速较高为管道流,水流可将骨料携带至堆积体前端,堆积稳定的骨料内部为低速渗流[17]。考虑管道携砂的阻力形式与明渠流相似,管流与明渠流流速分布相似[18],因此作用流速采用指数流速分布公式[19-20]:
$$ {U_{\text{b}}} = (1 + \beta ){\alpha ^\beta }U{\left(\frac{D}{h}\right)^\beta } $$ (40) 式中:α为近底颗粒高度系数,取值1/4;β为指数,取值1/6;h为水深,m。
根据斜面上骨料颗粒起动临界条件:
$$ F = G'\sin\; \psi + {F_{\text{D}}} $$ (41) $$ N = G'\cos\; \psi - {F_{\text{L}}} $$ (42) $$ F/N \geqslant f $$ (43) 式(41)—式(43)联立得:
$$\begin{array}{c} \tan\; \dfrac{\psi }{2} = \\ \dfrac{{ - G' + \sqrt {{{\left( {G'} \right)}^2} - \left[ { - {F_{\text{D}}} + f(G' + {F_{\text{L}}})} \right]\left[ { - {F_{\text{D}}} - f(G' - {F_{\text{L}}})} \right]} }}{{ - {F_{\text{D}}} + f(G' + {F_{\text{L}}})}} \end{array} $$ (44) 上述公式中,若动水流速相对于骨料粒径过大,计算结果会小于0,这与事实不符。因此若计算值小于0,将其规定为0。
淹没状态下,骨料颗粒处于饱和状态,颗粒间不存在黏聚力,水下休止角与内摩擦因数满足式(45),此时的内摩擦因数取决于颗粒形态特征及成分[21]:
$$ f = \tan\; \theta $$ (45) 式中:θ为水下休止角,由张红武[21]公式计算得到:
$$ \theta = 29.5 + 4.5\lg D $$ (46) 2.3 巷道顶区未接顶空间流速估算公式
巷道内主流区流速由初始流速决定,而影响瞬时携砂能力的主要因素为未接顶区所能达到的最大流速。掌握流速的变化规律可进行骨料粒径的动态匹配,保证堆积体快速长距离接顶[17],若顶区流速相对骨料粒径过低,则会引发堵孔问题。
设定初始动水流速Us,则矩形断面巷道内满流状态下单位时间流量为:
$$ Q = {U_{\text{s}}} L H $$ (47) 随着堆积高度y增加,剩余过水断面高度减小为H-y,假设堆积体内部孔隙率ε为0,属于完全不透水情况,则顶部流速U升高至:
$$ U = \frac{Q}{{L\left( {H - y} \right)}} $$ (48) 流速增幅ΔU = U−Us,牟林[9]采用此方法估算了截流过程中顶部未接顶空间的流速。事实上水流在堆积体内部低速渗流,孔隙率不会为0,若堆积体内孔隙率为1,堆积体不存在,则ΔU = 0。因此,根据孔隙率不同,流速增幅处于(0~1)ΔU范围内。假设ΔU与孔隙率ε线性相关,随ε增大,ΔU减小,则骨料堆积体顶部流速U可由下式计算得到:
$$ U = \left[ {\frac{Q}{{L\left( {H - y} \right)}} - {U_{\text{s}}}} \right]\left( {1 - \varepsilon } \right) + {U_{\text{s}}} $$ (49) 2.4 水流携砂能力及骨料留存临界流速计算公式
骨料留存临界流速Ucr关系着骨料灌注后能否留存堆积,同时,水流携砂能力qs与堆积能力(由灌注速度和骨料粒径所决定)的相互作用主导着堆积体能够达到的极限高度。因此对上述参量的把控,可以帮助预判骨料的留存和堆积状态,减少无效灌注。
单位时间内单位床面上水流势能损失E表示为[22]:
$$ E = {\tau _1}U $$ (50) $$ {\tau _1} = \tau (1 - {h_{\text{s}}}/h) $$ (51) 式中:τ1为主流区与床层颗粒交界面处的流体剪切应力,kN/m2;hs为床面可动层厚度,m;τ为流体剪切应力,kN/m2。
其中一部分能量用于床面颗粒的起动Ec,另一部分用于床面推移质输运Es,分别表示为:
$$ {E_{\text{c}}} = {\tau _{\text{c}}}{U_{\text{c}}} $$ (52) $$ {E_{\text{s}}} = \frac{{{\gamma _{\text{s}}} - \gamma }}{\gamma }{q_{\text{s}}}^ * \tan\; \theta $$ (53) 式中:τc为颗粒起动剪切应力,kN/m2。
E′表示除Es,Ec以外的能量损耗[23],根据能量守恒原理列出如下方程:
$$ E = {E_{\text{c}}} + {E_{\text{s}}} + E' $$ (54) 由于推移质输沙问题的复杂性,E′的形式尚且难以给出,因此依照物理学常用处理方式将式表示为[24]:
$$ {E_{\text{s}}} = \eta (E - {E_{\text{c}}}) $$ (55) 其中,η为输沙效率系数,η=ebFr。根据Bagnold假定,取eb=0.5[25];弗汝德数Fr作为一定水深矫正的水流强度指标[26]:
$$ {F_{\text{r}}} = \frac{U}{{\sqrt {gh} }} $$ (56) 将式(50)—式(53)代入式(55)可得:
$$ \frac{{{\gamma _{\text{s}}} - \gamma }}{\gamma }{q_{\text{s}}}^ * \tan\;\theta = \eta (\tau U - {\tau _{\text{c}}}{U_{\text{c}}}) $$ (57) 其中:
$$ \tau = \gamma RJ = \gamma \frac{{{U^2}}}{{{C^2}}} $$ (58) $$ {\tau _{\text{c}}} = \gamma R{J_{\text{c}}} = \gamma \frac{{U_{^{\text{C}}}^2}}{{{C^2}}} $$ (59) 式中:C为谢才系数;R为水力半径。
实际上由于水流的冲刷作用,不同流速下骨料堆积体休止角变化范围较大,在计算携砂能力时采用水下休止角略显粗糙,因此将式(57)中水下休止角θ改用式(44)计算得到的动水休止角ψ代替,将式(56)、式(58)、式(59)代入式(57)得水流携砂能力:
$$ {q_{\text{s}}}^ * = {e_{\text{b}}}\frac{{{\gamma ^2}}}{{{\gamma _{\text{s}}} - \gamma }}{\mathrm{ctan}} \;\psi \frac{{{U^3}(1 - {h_{\text{s}}}/H) - U_{\text{c}}^{\text{3}}}}{{{C^2}}}\frac{U}{{\sqrt {gh} }} $$ (60) 其中,床面可动层厚度采用Bagnold[27]建议的公式计算:
$$ {h_{\text{s}}} = 1.4{({U_*}/{U_{*{\text{c}}}})^{0.6}}D $$ (61) 式中:U*为水流摩阻流速,U*=(τ/ρ)1/2,m/s;U*c为泥沙起动摩阻流速,m/s。
以曼宁—斯特里克勒公式为基础,参考钱宁、麦乔威等学者取D65作为粗糙特征长度的处理,谢才系数C可表示为[28]:
$$ C = {A_{65}}{(\frac{R}{{{D_{65}}}})^{1/6}} $$ (62) 式中: D65为小于该粒径的颗粒占总质量百分比65%的粒径,m; A65为摩阻参数。
根据大量床沙资料统计结果,摩阻参数A65同泥沙中值粒径D50的关系式如下[23]:
$$ {A_{65}} = 1.55h\left( {\frac{{{D_{50}}}}{{{D_0}}}} \right) + 32.1 $$ (63) 式中:D50为中值粒径,m;D0为避免量纲不和谐而引入的参考粒径,为保持定量关系不变,取D0=1 m。
颗粒被起动时,除考虑拖曳力、上举力和水下重力外,还应考虑由于颗粒间相互影响以及颗粒在床面的相对暴露度而产生的附加质量力,因此何文社[29]提出颗粒临界起动流速Uc计算公式,上下限分别为:
$$ {U_{\text{c}}} = \sqrt {2\frac{{{\gamma _{\text{s}}} - \gamma }}{\gamma }g{D_{50}} + 5.5 \times {{10}^{ - 7}}\frac{{10 + h}}{{{D_{50}}^{0.7}}}} {\left( {\frac{h}{{{D_{50}}}}} \right)^{1/6}} $$ (64) $$ {U_{\text{c}}} = \sqrt {1.134\frac{{{\gamma _{\text{s}}} - \gamma }}{\gamma }g{D_{50}} + 3.1 \times {{10}^{ - 7}}\frac{{10 + h}}{{{D_{50}}^{0.7}}}} {\left( {\frac{h}{{{D_{50}}}}} \right)^{1/6}} $$ (65) 由于床面可动层厚度hs较薄,根据牟林[30]的研究,颗粒层厚度为((1~2)Dmax)时,起动流速靠近拟合曲线的下限,因此在后续计算中取颗粒起动流速下限。随着堆积体高度的增加,主流场进一步被压缩,流速升高,导致水流携砂能力迅速增加,当其与灌注速度达到动态平衡状态时,得到骨料留存的临界流速:
$$ {U_{{\text{cr}}}} = \frac{{\sqrt r + \sqrt {r - 2r + 4k\sqrt r } }}{2} $$ (66) 其中:
$$ r = \sqrt[3]{{ - \frac{q}{2} + \sqrt {\frac{{{q^2}}}{4} + \frac{{{p^3}}}{{27}}} }} + \sqrt[3]{{ - \frac{q}{2} - \sqrt {\frac{{{q^2}}}{4} + \frac{{{p^3}}}{{27}}} }} $$ (67) $$ k = \frac{{U_{\text{c}}^{\text{3}}}}{{2r(1 - {h_{\text{s}}}{\text{/}}H)}} $$ (68) $$ q = - {\left( {\frac{{U_{\text{c}}^{\text{3}}}}{{1 - {h_{\text{s}}}/H}}} \right)^2} $$ (69) $$ p = \frac{{4\phi {C^2}\tan\;\psi \sqrt {gh} }}{{{e_{\text{b}}}(1 - {h_{\text{s}}}/H)}}\frac{{({\gamma _{\text{s}}} - \gamma )}}{{{\gamma ^2}}} $$ (70) 式(60)中携砂能力qs*表示一定水流条件下,单位宽度上单位时间内通过断面的骨料颗粒重量,量纲为kg/(m·s)。取巷道宽度为L,则 qs=Lqs*,量纲与骨料灌注速度(kg/s)量纲一致。
3. 骨料灌注堆积数值模拟试验
结合CFD-DEM双向瞬态耦合仿真方法[31],研究骨料堆积对流场分布的影响,以及不同流场中堆积体形态差异,并对预测模型进行验证。此方法中流体运动采用欧拉方法描述,固体颗粒运动采用拉格朗日方法描述,可以同时反映颗粒–颗粒以及颗粒–流体间的相互作用,能够准确获取各网格单元的流速和压强数据,以及颗粒的位置、速度以及受力情况[32]。
3.1 模型建立
对流体计算域进行非结构化立方体网格剖分,网格尺寸30 mm。将骨料颗粒简化为球体,颗粒间接触选用Hertz-Mindlin无滑移“软球”模型[33],拖曳力模型采用精度较高的Gidaspow曳力模型[34],材料基本参数均按表1设定[12]。
表 1 计算模型材料参数Table 1. Material parameters used in the model颗粒参数 流体参数 颗粒密度ρs /
(kg·m−3)弹性模量E /
(MN·m−2)泊松比v 恢复系数e 滑动摩擦因数μr 滚动摩擦因数μf 密度ρ/
(kg·m−3)动力黏度μf/
(kg·m−1·s−1)2 650 5 0.3 0.3 0.5 0.1 1 000 0.001 实际工程中煤岩体巷道多为断面尺寸4×4 m的矩形巷道,本试验按相似比1∶20设置计算域,如图8所示,灌注速度Φ同样依据现场数据按比例1∶20取值[30]。
3.2 巷道内流场时空演化分析
以初始流速0.5 m/s动水环境中灌注中值粒径D50为15 mm骨料(粒径范围13~17 mm,为模拟一般性工况粒径级配按线性分布[30])为例,探究巷道内流速分布的变化规律。
如图9所示,骨料灌注期间,巷道顶部高速水流区和底部低速区同步扩大,由初始流速所决定的中部主流区(动水流速位于0.26~0.52 m/s之间)不断被“压缩”。如图9a所示,t=8 s时,骨料堆积体已接近巷道顶部,但由于堆积长度较短,流场受影响范围较小,程度较弱。t=10 s时,过水断面高度不足原来的20%,顶部高速水流区范围(动水流速大于1.04 m/s)扩大,背水坡处流场受到的压缩程度加剧。t=15 s水流携砂能力超过骨料堆积能力,骨料颗粒已无法稳定沉积在堆积体顶部。其中,突出在床面上的和位于堆积体前沿的颗粒都是最容易起动的,从图9c中可以发现以上两位置的颗粒速度均已接近0.75 m/s,颗粒在堆积体顶部表面做接触质、跃移质运动。
为定量分析骨料堆积生长对断面流速分布的影响,创建监测线实时监测巷道x=1.2 m处底部到顶部(z=0.2~0.18 m)的流速变化情况。如图10所示,t=2 s时,骨料开始堆积,水流流经堆积体时受阻,进而从堆积体顶部“绕行”,导致堆积高度以上区域流速增大,堆积体内部流速降低。t=15 s时,骨料堆积稳定,高度达到极限,堆积体内部为低速渗流流态,且不同高度流速大致相同,稳定在0.09~0.11 m/s之间,由此证明了1.1基本假定2的合理性。
3.3 骨料堆积体形态差异分析
分别在初始流速0.3、0.5、0.8 m/s及静水条件下以相同灌注速度灌注中值粒径10、15、20 mm骨料颗粒,测试其堆积形态的差异。
如图11中虚线所示,根据张红武公式[23],静水中3种粒径骨料颗粒的水下休止角分别为:33.57°、34.56°、35.27°。理论公式值,模拟试验值及张红武公式值吻合程度较好,相对误差最大值仅为1.26%。
如图11中实线所示,动水休止角ψ随流速增大呈抛物线型趋势减小,根据詹义正[19]的研究,将动水休止角为0时所对应的流速定义为平床流速Up。3种粒径颗粒的Up分别为:1.03、1.19、1.33 m/s。流速与平床流速比值K为相对流速,当K∈(0, 0.5)时,动水休止角的大致范围为27°~35°(图12),颗粒自身内摩擦系数作为影响动水休止角的主导因素,内摩擦系数越大则动水休止角越大,此时流速对动水休止角影响较小。当K∈(0.5, 1)时,流速对动水休止角的影响作用不断增强,相同流速条件下,粗颗粒的动水休止角大于细颗粒,对流速变化的敏感度低于细颗粒(图11、12),该趋势与文献[19]中得到的结论一致,但流速增大对于堆积体内部孔隙率的影响较弱,如表2所示,最大变化值仅有2.02%。
表 2 骨料堆积体内部孔隙率Table 2. Porosity of aggregate accumulation粒径/mm 初始流速/(m·s−1) 孔隙率/% 10 0 41.17 0.3 42.48 0.5 42.74 0.8 43.06 15 0 41.53 0.3 43.05 0.5 43.27 0.8 43.55 20 0 42.25 0.3 42.74 0.5 43.02 0.8 43.69 3.4 骨料堆积生长预测模型验证分析
1) 试验方案设计。依据以往工程案例,截流初期巷道内流速一般为0~0.5 m/s,但骨料堆积临近接顶时,顶区流速可达初始流速的数倍以上[4]。因此,同时在高流速(流速大于 0.5 m/s)条件下开展骨料堆积过程及截流阻水机理研究对截流工程施工更具有指导意义[17]。截流堵水工作中,由于初期流速较小,先灌注5~10 mm (小石子)骨料进行铺底,抬升巷道内水流流速,然后灌注粒径10~30 mm(中石子)骨料对巷道进行充填[35],接顶阶段会混合灌注粒径超过30 mm(大石子)的骨料。重点研究内容为充填阶段的骨料堆积生长,因此选取中值粒径D50为15 mm 的骨料颗粒开展模拟试验。试验过程中保持进水口初始流速及灌注速度恒定不变,试验方案见表3。
表 3 模拟试验方案Table 3. Simulation test scheme编号 初始流速/(m·s−1) 粒径/mm 灌注速度/(kg·s−1) 方案一 0.5 15 0.5、1.0、2.0 方案二 0.3、0.5、0.8 15 1.0 2) 结果对比分析。预测模型中各参变量取值与模拟试验一致,关键参量取值由理论公式计算得到,孔隙率ε根据3.3节模拟试验结果取值。由图13可以看出预测模型能够较好的反映堆积体高度及长度随时间的演化规律。灌注速度越大,水流速度越低,堆积高度增长速度越快,第1阶段持续时间越短;相同灌注速度下,初始流速越高,堆积体极限高度越小,堆积长度增速越快。如图13a所示,灌注开始后,堆积高度呈线性迅速增长,Φ=2.0 kg/s条件下,t=8 s时,已达极值(168.33 mm)的93.5%,此时骨料堆积生长仍处于高度快速增长阶段;8 s后过渡至第2阶段,高度增速放缓呈抛物线型增长趋势,长度与高度同步生长;直至水流携砂能力qs与灌注速度Φ达到平衡,顶区流速U达到骨料留存临界流速Ucr,生长进入第3阶段,t=15 s时,堆积体垂向高度已达到极限。在水平长度方面,其初期增长呈“凹”型曲线(图13b、图13d);若堆积体生长顺利过渡至第3阶段,堆积长度将呈线性增长趋势。若初始流速与灌注条件不匹配,灌注过程将面临堵孔风险,如图13c,t=25 s时,骨料堵塞投料孔,而此时堆积段长度不足,增势趋于平缓(图13d中红色曲线),导致堆积体抵抗水力冲击的能力较弱。
灌注初期,堆积高度模拟试验值略低于理论公式计算值(图13a、图13c),长度方面则情况相反(图13b、图13d)。分析原因:由于前期堆积还未形成稳定结构,颗粒间存在相对移动,空隙较大,且实际的坡面角度小于理论公式动水休止角计算值,因此两结果的相对误差较大(表4)。但随着骨料的持续灌注,内部骨料受到压力不断增大,颗粒相互挤压,空隙随之减小达到稳定,堆积体坡角接近动水休止角理论计算值。10 s后,预测模型值与模拟试验值之间的“误差带”带宽逐渐缩小(图13),相对误差最大值(ey,max、el,max)均小于10%,由此可以认定截流初期的结果误差并未对后续计算产生较大影响,预测模型具有较高的计算精度。
表 4 预测模型与模拟试验误差统计Table 4. Statistical analysis of relative errors between prediction model and simulation testt /s ey,max / % el,max / % 方案一 方案二 方案一 方案二 5 15.2 18.0 16.9 17.7 10 6.6 6.6 8.8 9.9 15 3.0 3.2 9.7 8.6 20 2.7 1.9 9.3 8.1 25 3.0 3.8 8.5 7.7 30 1.0 1.0 6.6 3.4 35 1.6 1.6 7.4 5.5 40 1.8 0.8 4.4 5.4 45 1.6 0.7 2.6 4.5 50 1.2 1.1 6.0 3.3 注:ey,max为堆积高度相对误差最大值;el,max为堆积长度相对误差最大值 4. 结 论
1) 分析了灌注初期骨料沉积主域的所在位置与初始流速间响应关系,将其在空间上分为3个区域,若沉积主域位于区域③,则标志着骨料灌注后无法留存堆积,属于无效灌注。截流降速工程中应根据水流条件制定合理的骨料灌注方案,包括灌注速度、骨料粒径。
2) 研究了骨料灌注过程中由于堆积体床面颗粒沿休止角塌落,流速抬升导致颗粒沉降时水平移距增大,以及颗粒着床后被高速水流二次起动带离等因素在堆积体生长过程中的作用机理,据此提出了静水和动水条件下堆积高度及长度预测模型;分析了导致堵孔现象发生的原因,建立了动水中骨料灌注堵塞投料孔的临界判据。
3) 考虑了动水环境中骨料堆积体内孔隙率等因素,得出了巷道顶部未接顶区流速估算方法;对堆积体坡面颗粒受力分析,提出了骨料动水休止角计算公式;根据颗粒沉降时的受力情况,提出了颗粒水平移距计算公式;基于床面层能量平衡方程,提出了水流携砂能力和骨料留存临界流速计算公式。
4) 建立了考虑初始流速、灌注速度、骨料粒径、骨料摩擦系数和水流黏度系数等因素的数值模型,模拟了骨料在动水环境中的堆积过程,研究了流场随骨料灌注的演化规律,对比分析了不同流场中骨料堆积形态的差异,验证了预测模型计算结果的可靠性。
-
表 1 计算模型材料参数
Table 1 Material parameters used in the model
颗粒参数 流体参数 颗粒密度ρs /
(kg·m−3)弹性模量E /
(MN·m−2)泊松比v 恢复系数e 滑动摩擦因数μr 滚动摩擦因数μf 密度ρ/
(kg·m−3)动力黏度μf/
(kg·m−1·s−1)2 650 5 0.3 0.3 0.5 0.1 1 000 0.001 表 2 骨料堆积体内部孔隙率
Table 2 Porosity of aggregate accumulation
粒径/mm 初始流速/(m·s−1) 孔隙率/% 10 0 41.17 0.3 42.48 0.5 42.74 0.8 43.06 15 0 41.53 0.3 43.05 0.5 43.27 0.8 43.55 20 0 42.25 0.3 42.74 0.5 43.02 0.8 43.69 表 3 模拟试验方案
Table 3 Simulation test scheme
编号 初始流速/(m·s−1) 粒径/mm 灌注速度/(kg·s−1) 方案一 0.5 15 0.5、1.0、2.0 方案二 0.3、0.5、0.8 15 1.0 表 4 预测模型与模拟试验误差统计
Table 4 Statistical analysis of relative errors between prediction model and simulation test
t /s ey,max / % el,max / % 方案一 方案二 方案一 方案二 5 15.2 18.0 16.9 17.7 10 6.6 6.6 8.8 9.9 15 3.0 3.2 9.7 8.6 20 2.7 1.9 9.3 8.1 25 3.0 3.8 8.5 7.7 30 1.0 1.0 6.6 3.4 35 1.6 1.6 7.4 5.5 40 1.8 0.8 4.4 5.4 45 1.6 0.7 2.6 4.5 50 1.2 1.1 6.0 3.3 注:ey,max为堆积高度相对误差最大值;el,max为堆积长度相对误差最大值 -
[1] 姬中奎. 矿井大流量动水注浆细骨料截流技术[J]. 煤炭工程,2014,46(7):43−45. JI Zhongkui. Fine aggregate grouting closure technology of high flow running water in mine[J]. Coal Engineering,2014,46(7):43−45.
[2] 邵红旗,王 维. 双液注浆法快速建造阻水墙封堵突水巷道[J]. 煤矿安全,2011,42(11):40−43. SHAO Hongqi,WANG Wei. Rapid construction of water blocking wall by double liquid grouting method to block water inrush roadway[J]. Safety in Coal Mines,2011,42(11):40−43.
[3] 牟 林. 过水大断面截流堵巷工程若干关键技术问题的探讨[J]. 煤炭学报,2021,46(11):3525−3535. MOU Lin. Discussion on key technical problems of water cutting-off engineering in submerged roadway with large-cross section[J]. Journal of China Coal Society,2021,46(11):3525−3535.
[4] 董书宁,牟 林. 突水淹没矿井动水巷道截流阻水墙建造技术研究[J]. 煤炭科学技术,2021,49(1):294−303. DONG Shuning,MOU Lin. Study on construction technology of water blocking wall in hydrodynamic pathway of submerged mine due to water inrush[J]. Coal Science and Technology,2021,49(1):294−303.
[5] 杨志斌,董书宁. 动水大通道突水灾害治理关键技术[J]. 煤炭科学技术,2018,46(4):110−116. YANG Zhibin,DONG Shuning. Key technology of water inrush disaster control under hydrodynamic large channel condition[J]. Coal Science and Technology,2018,46(4):110−116.
[6] 惠 爽. 矿井淹没巷道多孔灌注骨料封堵模拟试验[D]. 徐州:中国矿业大学,2018:14-57. HUI Shuang. An experimental investigation on pouring aggregate to plug an inundated mine tunnel through boreholes[D]. Xuzhou:China University of Mining and Technology,2018:14-57.
[7] 李 军. 地铁突水抢险中粗骨料注浆快速封堵数值模拟研究[D]. 徐州:中国矿业大学,2021:21-48. LI Jun. Numerical simulation study on rapid plugging of coarse aggregate grouting in subway water inrush rescue[D]. Xuzhou:China University of Mining and Technology,2021:21-48.
[8] 沈 淇,顾峰峰,万远扬,等. 基于泥沙随机交换过程及不同推移形式的推移质输沙公式[J]. 应用基础与工程科学学报,2023,31(3):599−610. SHEN Qi,GU Fengfeng,WAN Yuanyang,et al. A formula for bedload transport based on the stochastic sediment interchanges and the different motion patterns[J]. Journal of Basic Science and Engineering,2023,31(3):599−610.
[9] 牟 林,董书宁. 截流巷道骨料堆积体中浆液运移规律与阻水机制[J]. 地下空间与工程学报,2020,16(6):1891−1900. MOU Lin,DONG Shuning. Migration rule and water blocking mechanism of cement slurry in aggregate accumulation of underground tunnel closure[J]. Chinese Journal of Underground Space and Engineering,2020,16(6):1891−1900.
[10] 陈光国,阳 宁,唐达生,等. 垂直管道颗粒及颗粒群沉降运动规律研究[J]. 泥沙研究,2010(4):16−21. CHEN Guangguo,YANG Ning,TANG Dasheng,et al. Study on the settling regularity of solid particles in vertical pipelines[J]. Journal of Sediment Research,2010(4):16−21.
[11] 苏东升. 基于CFD-DEM耦合模拟方法的水流泥沙运动研究[D]. 天津:天津大学,2015:8−21. SU Dongsheng. Study on flow and sediment movement based on CFD-DEM coupling simulation method[D]. Tianjin:Tianjin University,2015:8−21.
[12] 周宗青,李利平,石少帅,等. 隧道突涌水机制与渗透破坏灾变过程模拟研究[J]. 岩土力学,2020,41(11):3621−3631. ZHOU Zongqing,LI Liping,SHI Shaoshuai,et al. Study on tunnel water inrush mechanism and simulation of seepage failure process[J]. Rock and Soil Mechanics,2020,41(11):3621−3631.
[13] 詹义正,寇树萍. 球体的移距及稳定移距公式[J]. 武汉水利电力大学学报,1996,29(2):85−90. ZHAN Yizheng,KOU Shuping. Equations of sphere’s falling distance[J]. Journal of Wuhan University of Hydraulic and Electric Engineering,1996,29(2):85−90.
[14] 姚仕明,梁 兰,刘卫峰,等. 抛石移距规律初探[J]. 武汉水利电力大学学报,1997,30(6):24−27. YAO Shiming,LIANG Lan,LIU Weifeng,et al. Preliminary study of throwing distance of riprap[J]. Engineering Journal of Wuhan University,1997,30(6):24−27.
[15] 姚令侃,李仕雄,蒋良潍. 自组织临界性及其在散粒体研究中的应用[J]. 四川大学学报(工程科学版),2003,35(1):8−14. YAO Lingkan,LI Shixiong,JIANG Liangwei. Self-organized criticality and its application in granular mixtures[J]. Journal of Sichuan University (Engineering Science Edition),2003,35(1):8−14.
[16] 何文社,方 铎,曹叔尤,等. 泥沙起动判别标准探讨[J]. 水科学进展,2003,14(2):143−146. HE Wenshe,FANG Duo,CAO Shuyou,et al. Study on standards for incipient motion of sediment[J]. Advances In Water Science,2003,14(2):143−146.
[17] 牟 林. 动水条件巷道截流阻水墙建造机制与关键技术研究[D]. 北京:煤炭科学研究总院,2021:29−92. MOU Lin. Study on construction mechanism and key technology of water-blocking wall in hydrodynamic pathway[D]. Beijing:Chinese Institute of Coal Science,2021:29−92.
[18] 范 玉. 水沙两相流理论及引水工程管道输沙问题的研究[D]. 天津:天津大学,2014:88−110. FAN Yu. Study on the theory of water-sediment two-phase flow and the problem of sediment transport in the pipeline of water diversion project[D]. Tianjin:Tianjin University,2014:88−110.
[19] 陈 立,宋 涛,李东锋,等. 侧向水流作用下均匀沙休止角变化的试验研究[J]. 泥沙研究,2017,42(3):1−6. CHEN Li,SONG Tao,LI Dongfeng,et al. Experiment study on repose angle of uniform sand under influence of lateral flow[J]. Journal of Sediment Research,2017,42(3):1−6.
[20] 马子普,张宝森,邓 宇,等. 适用于明渠流及冰盖流的统一泥沙颗粒起动流速公式[J]. 水利学报,2021,52(8):969−978. MA Zipu,ZHANG Baosen,DENG Yu,et al. Unified incipient velocity formula of non-cohesive sediments for open channel flow and ice-cover flow[J]. Journal of Hydraulic Engineering,2021,52(8):969−978.
[21] 张红武,汪家寅. 沙石及模型沙水下休止角试验研究[J]. 泥沙研究,1989(3):90−96. ZHANG Hongwu,WANG Jiayin. Experimental study on underwater angle of repose of sand and model sand[J]. Journal of Sediment Research,1989(3):90−96.
[22] 夏华永,廖世智,肖志建. 基于床面层能量平衡关系的推移质输沙率计算式[J]. 水利水运工程学报,2006(4):1−9. XIA Huayong,LIAO Shizhi,XIAO Zhijian. Formula for bedload sediment transport rate based on energy balance relationship within bed-surface layer[J]. Hydro-Science and Engineering,2006(4):1−9.
[23] 张红武,张俊华,卜海磊,等. 试论推移质输沙率公式[J]. 南水北调与水利科技,2011,9(6):140−145. ZHANG Hongwu,ZHANG Junhua,BU Hailei,et al. Discussion of bed-load transport equations[J]. South-to-North Water Diversion and Water Science & Technology,2011,9(6):140−145.
[24] 张红武,张 清. 黄河水流挟沙力的计算公式[J]. 人民黄河,1992,14(11):7−9,61. ZHANG Hongwu,ZHANG Qing. Formula of sediment carrying capacity of the Yellow River[J]. Yellow River,1992,14(11):7−9,61.
[25] Bagnold. The flow of cohesionless grains in fluids[J]. Philosophical Transactions of the Royal Society of London Series A,Mathematical and Physical Sciences,1956,249(964):235−297.
[26] 贺梦杨,尚海鑫,张宽地,等. 仿机翼形便携式量水槽水力特性试验与数值模拟[J]. 农业工程学报,2021,37(12):117−124. HE Mengyang,SHANG Haixin,ZHANG Kuandi,et al. Hydraulic performance experiments and numerical simulation of portable water measuring flume of imitating airfoil shape[J]. Transactions of the Chinese Society of Agricultural Engineering,2021,37(12):117−124.
[27] BAGNOLD. The nature of saltation and of ‘bed-load’ transport in water[J]. Proceedings of the Royal Society of London A Mathematical and Physical Sciences,1973,332(1591):473−504. doi: 10.1098/rspa.1973.0038
[28] 张红武,张罗号,彭 昊,等. 冲积河流糙率由来与计算方法研究[J]. 水利学报,2020,51(7):774−787. ZHANG Hongwu,ZHANG Luohao,PENG Hao,et al. Research on cognition and calculation method of alluvial river roughness[J]. Journal of Hydraulic Engineering,2020,51(7):774−787.
[29] 何文社,方 铎,杨具瑞,等. 泥沙起动流速研究[J]. 水利学报,2002,33(10):51−56. HE Wenshe,FANG Duo,YANG Jurui,et al. Study on incipient velocity of sediment[J]. Journal of Hydraulic Engineering,2002,33(10):51−56.
[30] 牟 林,董书宁,郑士田,等. 基于CFD-DEM耦合模型的阻水墙建造过程数值模拟[J]. 岩土工程学报,2021,43(3):481-491,F0002,F0003. MOU Lin,DONG Shuning,ZHENG Shitian,et al. Numerical simulation of construction of water-blocking wall based on CFD-DEM coupling method[J]. Chinese Journal of Geotechnical Engineering,2021,43(3):481-491,F0002,F0003.
[31] 王 胤,艾 军,杨 庆. 考虑粒间滚动阻力的CFD-DEM流–固耦合数值模拟方法[J]. 岩土力学,2017,38(6):1771−1780. WANG Yin,AI Jun,YANG Qing. A CFD-DEM coupled method incorporating soil inter-particle rolling resistance[J]. Rock and Soil Mechanics,2017,38(6):1771−1780.
[32] 刘巨保,王 明,王雪飞,等. 颗粒群碰撞搜索及CFD-DEM耦合分域求解的推进算法研究[J]. 力学学报,2021,53(6):1569−1585. LIU Jubao,WANG Ming,WANG Xuefei,et al. Research on particle swarm collision search and advancement algorithm for cfd-Dem coupling domain solving[J]. Chinese Journal of Theoretical and Applied Mechanics,2021,53(6):1569−1585.
[33] TSUJI Y,TANAKA T,ISHIDA T. Lagrangian numerical simulation of plug flow of cohesionless particles in a horizontal pipe[J]. Powder Technology,1992,71(3):239−250. doi: 10.1016/0032-5910(92)88030-L
[34] GIDASPOW D. Hydrodynamics of fiuidizatlon and heat transfer:supercomputer modeling[J]. Applied Mechanics Reviews,1986,39(1):1−23. doi: 10.1115/1.3143702
[35] 蒋勤明. 大埋深突水巷道“阻水段”骨料灌注技术[C]//2017年钻探工程学术研讨会论文集. 鄂尔多斯:中国煤炭学会钻探工程专业委员会,2017:210−212. JIANG Qinming. Aggregate pouring technology for "water blocking section" of deep water inrush pathway[C]//Proceedings of the 2017 Symposium on Drilling Engineering. Ordos:Drilling Engineering Professional Committee of China Coal Society,2017:210−212.