含瓦斯煤卸压拉伸响应：释放时序、扩散时间尺度与力学约束的影响

《Results in Engineering》：Tensile response of gas-bearing coal to pressure release: Effects of release timing, diffusion timescale and mechanical restraint

　　

匹配的压降和持续时间既可以增强也可以抑制拉伸响应；仅凭平均释放速率本身无法决定拉伸响应；第二时间矩（second temporal moment）使36个测试案例的均方根误差（RMSE）降低79.5%；横向约束（transverse restraint）改变了增强与抑制之间的边界；具有有限上升时间（finite rise time）和Biot耦合（Biot coupling）的选定测试保持了响应符号。

含瓦斯煤卸压过程中的拉伸响应受到释放时序、扩散时间尺度与力学约束的共同影响。煤与瓦斯突出涉及气体饱和多孔固体的快速破坏，已有能量分析、多物理场耦合及物理实验表明，卸压速率与拉伸损伤有关。然而，在多数实验中，压力、水分、装置体积、阀门响应、试件强度和释放持续时间等因素同时变化，无法判定具有相同总压降、相同持续时间和相同名义平均释放速率的不同压力历史是否产生等效的力学响应。为此，研究人员建立了线性孔隙弹性模型，在固定压降和释放持续时间条件下，比较不同释放波形对内部拉伸应力增量的影响，并考察扩散时间尺度、横向约束及Biot耦合的作用。研究结果表明，匹配的压降和持续时间既能增强也能抑制拉伸响应；平均释放速率单独不能决定响应；第二时间矩可显著改善预测；横向约束改变增强与抑制的边界；有限上升时间和Biot耦合不改变所选历史的响应方向。这一研究为评价卸压历史对煤体稳定性的影响提供了定量框架，强调应记录完整的压力–时间曲线和约束条件。

研究人员采用数值模拟方法（无样本队列），在单位正方形区域（归一化坐标X=x/L，Y=y/L）内设置一个位于暴露边界中央、占边界三分之一长度的卸压窗口，其余压力边界不渗透。孔隙压力扩散采用常扩散系数的线性扩散方程，用伽辽金有限元（线性三角形）离散，时间积分采用带行和集中质量矩阵的向后欧拉法。力学求解为准静态平面应变弹性，本构关系为Δσ=C:ε-αδpI，右边界法向位移增量为零，上、下边界用无量纲横向弹簧刚度κ_s=k_sL/E模拟不同约束程度。研究还实现了完全耦合的线性Biot公式与单向扩散模型的对比，耦合数C_B=α²M/K_d。主要无量纲参数为E/ΔP=20，ν=0.25，α=0.80；释放傅里叶数F_ou=cτ_u/L²取0.02、0.10和0.50，代表快速、中等和慢速释放相对扩散时间尺度。

研究结果部分按论文标题分述：

4 扩散时间尺度对响应转变的影响。4.1 压力记忆与应力重分布：有限卸压窗口在材料内部产生二维扩散场；相同终点压力下，不同历史在释放结束时保留不同的内部压力场，并通过全局平衡转化为应力差异。4.2 直接计算的响应图：对49种跳–保–斜坡（jump–hold–ramp, B3）组合直接计算，在平面应变和横向辊约束下，P₉₅最大主应力相对效应在F_ou=0.02时从-43.76%到+21.17%，在F_ou=0.10时从-22.20%到+11.08%，在F_ou=0.50时从-14.70%到+8.49%；较大的初始跳跃趋向增强，较长的保持趋向抑制，零效应曲线随F_ou移动；有限上升时间ε=0.025在F_ou=0.10下使零效应曲线平均移动0.00238，拓扑不变。4.3 匹配平均速率的不足：因为所有历史具有相同的ΔP/τ_u，若平均速率足以描述输入，相对效应应为零；实际计算范围跨越正负符号且比时间步进误差大两个数量级，表明需要更高阶的时序信息。

5 时间矩跨波形族评估。5.1 模型定义与拟合系列：对B3系列中9个案例拟合无截距的矩模型ê_P95=a1(F_ou)Δm₁+a2(F_ou)Δm₂，训练RMSE在0.199到0.674个百分点之间，留一法误差在0.241到0.814个百分点之间；二次曲面留一法误差更大，条件数57.8远大于矩模型的5.015。5.2 测试波形结果：对12种不同单调历史（4种指数、4种幂律、4种两段线性）在3个F_ou下共36个测试案例进行无重拟合预测，两矩模型得到整体RMSE为0.871个百分点，最大绝对误差2.934个百分点，36个案例符号分类全部正确；仅用一阶矩的模型RMSE为4.252个百分点，第二矩使RMSE降低79.5%。5.3 降维的诠释与局限：压力历史经扩散核卷积后，应力指标汇总空间滤波结果；两矩模型在所指定线性模型和波形范围内有效，但不适用于振荡释放或非线性扩散等情形。

6 连续约束与孔隙弹性耦合的稳健性。6.1 平面应变下最大主应力随横向弹簧刚度变化：对于增强、近转变和抑制历史，P₉₅效应分别从+6.54%到+9.38%、-3.78%到-3.94%、-9.95%到-12.12%，符号不变，更强的约束增大分离；零效应边界从牵引自由极限到辊约束移动，例如在F_ou=0.10、r_h=0.80时χ₀从0.69151变为0.65146。6.2 应力指标选择造成的表观反转：仅使用面内最大主应力时，增强和抑制历史会出现跨零反转；而重建平面应变下的完整张量最大主应力后，所选历史符号保持稳定，表明必须指定张量定义。6.3 对ROI、百分位数和应力分量的稳健性：使用4种ROI、4种统计量（均值、P90、P95、P99）和两种应力定义共32种组合，所选增强和抑制历史在所有组合中方向一致；大小范围变化较大，方向比幅值更稳健。6.4 完全耦合线性Biot响应：在F_ou=0.10下，C_B从0.05增至0.50使增强效应从+9.34%到+10.00%，近转变效应从-4.09%到-4.71%，抑制效应从-12.30%到-13.66%，符号不变；单向模型在强耦合下略微低估幅值。

7 工程诠释与局限。7.1 边界历史作为工程输入：卸压边界应指定完整的时间历史，而不仅仅是压降和持续时间；需要记录压力–时间曲线，并用F_ou、Δm₁、Δm₂作为紧凑摘要；横向约束也应一并说明；实验验证需同时测量变形、声发射或数字图像相关等，并固定ΔP、τ_u、几何、材料批次、水分和约束。7.2 范围与局限：计算量是相对B1的增量拉伸应力；绝对破裂条件需要叠加初始有效应力并使用强度准则；模型采用常扩散系数、均质材料、平面应变、中心卸压窗口，未包含塑性、损伤、渗透率演化、惯性等；测试历史均为单调单次释放。

讨论部分总结：研究人员指出，名义平均释放速率不能唯一描述边界输入，相同平均速率的历史可产生相反的增量应力；两时间矩在指定波形族内给出良好预测；横向约束和应力张量定义影响临界分类；单向模型在有限Biot耦合范围内保持方向。研究结论翻译如下：1. 在有限卸压窗口模型中，具有相同总压降、持续时间和名义平均释放速率的压力历史产生不同的力学响应；直接计算49种B3历史，平面应变下释放结束时的P₉₅效应包含增强和抑制，全局极值为-43.76%和+21.17%；在F_ou=0.10时，有限上升分数ε=0.025保持零效应曲线拓扑，共同轮廓节点上跳跃分数的平均绝对移动为0.00238。2. 扩散时间尺度控制响应幅值和转变位置；相对扩散更快的释放保留更大的边界历史差异；零效应曲线在所有三个F_ou下均存在。3. 释放的两个时间矩在测试范围内推广到拟合边界族之外；36个测试案例（12种不同历史×3个F_ou）无重拟合下符号分类准确率100%，RMSE为0.871个百分点，最大绝对误差2.934个百分点；第二矩比仅用第一矩的模型降低RMSE 79.5%。4. 连续横向约束改变响应幅值并移动零效应边界；在F_ou=0.10、r_h=0.80时，χ₀从牵引自由极限的0.69151变为辊极限的0.65146；当从平面应变完整张量重建最大主应力时，所选增强和抑制效应保持符号；仅使用面内应力会出现表观反转。5. 有限范围的完全耦合线性Biot计算（C_B=0.05–0.50）产生中等幅值变化，所选增强和抑制效应保持符号。6. 这些结论适用于指定线性模型中的增量拉伸响应；煤的剥落或突出评估还需初始应力、强度、非均质性和损伤/裂纹扩展物理。