Molecular dynamics study on wettability difference between lignite dust and silica dust
-
摘要:
利用分子动力学模拟和量子化学方法，研究了褐煤煤尘与矽尘润湿机理的差异性；构建了水−表面活性剂−褐煤煤尘/矽尘模型，通过分子动力学模拟研究平衡后100 ps内各组分间相互作用能、径向分布函数（RDF）、均方位移（MSD）、自扩散系数（D）和相对数密度。结果表明：相对于矽尘，十六氨基三甲基溴化铵（CTAB）溶液与褐煤分子的相互作用能和自扩散系数较大、相对数密度峰值较高；基于褐煤煤尘与矽尘的静电势，SiO2的电势差变化幅度较大，可知活性剂吸附导致煤尘表面与水分子的电势差大小是影响表面润湿性差异的因素之一，但褐煤中官能团增加了分子间的相互作用，并在润湿过程中发挥着巨大作用。
Abstract:In this paper, the differences in the wetting mechanisms of lignite dust and silica dust were investigated using molecular dynamics simulations and quantum chemical methods. A water-surfactant-lignite coal dust/silica dust model was constructed and the interaction energy, radial distribution function (RDF), mean square displacement (MSD), self-diffusion coefficient (D) and relative number density among the components within 100 ps after equilibrium were investigated by molecular dynamics simulations. The results show that the interaction energy and self-diffusion coefficient between cetyl trimethyl ammonium bromide(CTAB) solution and lignite molecules are larger and the peak value of relative number density is higher compared to silica dust; by comparing the electrostatic potential of lignite dust with that of silica dust, the magnitude of the potential difference of SiO2 varies widely, which shows that the size of the potential difference between the dust surface and water molecules due to active agent adsorption is one of the factors affecting the difference in surface wettability, but the functional groups in lignite increase the intermolecular interactions and play a significant role in the wetting process.
-
Keywords:
- silica dust /
- coal dust /
- wettability /
- molecular dynamics /
- quantum chemistry /
- dust control
-
声波传输是实现矿用声波随钻测量的关键，钻柱作为信息传输通道，由不同尺寸的钻杆相连接而成，以钻杆中传播的声波为信息载体进行数据传输[1]。为了改善目前煤矿井下瓦斯抽采的钻孔施工中大部分处于“盲钻”的问题，有必要研究煤矿井下钻杆中声波的传输机制和特性。
国内外学者针对石油随钻信息传输技术领域的钻杆柱中声波传播特性开展了大量研究。针对钻杆中声波的传播特性和频谱特性的研究表明，物性参数和钻杆结构对声波的传输过程有显著影响，印证了钻杆用作煤矿井下无线数据传输信道的可行性[2]。为进一步深入了解钻杆内声波的传播规律，通过有限元法分析影响电磁随钻测量信号强度的因素，针对钻杆等周期性管柱中的信息传输，设计井下声波发生器以有效地产生声波和高效地传输声能，建立频谱特性的判断标准，通过评价函数对通带和频谱特性进行定量表征[3-5]，有助于选取最优载波频率，减少井下因素对钻柱信号的不良效应。在实际工况中，钻柱振动会引起钻柱故障、井眼轨迹恶化、钻头过度磨损等问题，学者们通过建立振动模型和实验钻柱系统，从轴向力、钻柱水平位移和摩擦扭矩等角度分析钻柱在不同振动形式下的特性，确定导致钻井过程不稳定的主要因素[6-9]，更好地降低水平井钻井时的钻柱振动。此外，考虑到钻柱波与井孔环境的相互作用，学者们通过建立物理模型和理论计算模型，研究了弯曲钻柱声波传播特性，得出了不同声源频率下声波信号的处理方法、衰减规律和有效传输距离[10-11]；为进一步深入了解钻杆内声波的传播规律和井下信息传输的载体信号频率选择提供了参考。
考虑到煤矿井下实际工况条件的差异、不同信号加载方式等因素影响，钻杆柱中声波信号衰减和色散等问题也不可避免地有所差别，因此，有必要研究煤矿井下钻杆柱中声波传输特性。为此，基于声波在周期性结构中的波动方程分析声波在钻杆柱中的传播特性；通过有限元分析软件对声波在钻杆中的传播过程进行仿真，分析输出波形幅值衰减特性以及频散特性；搭建试验平台得到了声波幅值衰减特性和波形变化规律，验证了理论分析和建模仿真结果；研究结果将为声波传输系统的设计和通信的模式选择、编码与波形识别提供理论依据。
1. 钻杆柱中声波传播特性
对实际结构进行调研发现，钻具之间通过接头连接，因此钻柱并不具有本质上的连续性。针对钻杆中的声波传播特性研究，通常将钻杆看作周期性管串结构，仅考虑由等长度的钻杆组成的钻柱[1-2,5]，周期性理想钻杆模型如图1所示。钻杆主要由管体和接头连接而成，管体部分为等截面的圆管，接头处截面积发生改变。结合周期性理想钻杆模型，建立等效周期性系统，周期性钻杆系统结构示意图如图2所示。图中：Li、Si、ρi分别为第i段的长度、横截面积、材料密度；uti为第i个单元内透射波的位移；uri为第i个单元内反射波的位移。
由图2可得：i为奇数时为钻杆段时，Li= L1，Si= S1，ρi=ρ1；i为偶数时为接头段时，Li= L2，Si= S2，ρi=ρ2。
1.1 声波在周期性结构中的波动方程
纵波常作为井下声波传输的载波形式，利用纵波在圆截面的直杆中的传播特性来近似表示声波在周期性结构钻柱中的传播特性。假设钻柱的阻抗为密度ρ和横截面积S的综合函数为F，针对纵波在钻柱结构中的传播，仅考虑质点的轴向位移$ u(x,t) $，钻杆中纵波波动方程为[12]：
$$ \rho s\frac{{\partial v}}{{\partial x}} = - \frac{{\partial F}}{{\partial x}} $$ (1) 式中：v为质点速度，$ v = \dfrac{{\partial u}}{{\partial t}} $；t为时间。
用x代替位置坐标，将油管柱的单位质量表示为：
$$ m = \int_0^{{x}} \rho (\xi )S(\xi ){\text{d}}\xi $$ (2) 式中：m为油管柱的质量坐标。
对式（2）求偏导可得：
$$ \frac{\partial F}{{\partial {{x}}}} = \frac{{\partial m}}{{\partial {{x}}}} \cdot \frac{\partial F}{{\partial m}} = \rho S\frac{\partial F}{{\partial m}} $$ (3) 式（1）在质量坐标下转换为：
$$ \frac{{\partial v}}{{\partial t}} = - \frac{{\partial F}}{{\partial m}} $$ (4) 式中：$ F = - {{\textit{z}}^2}\dfrac{{\partial u}}{{\partial m}} $，z为阻抗，$ {\textit{z}} = S\rho c $；c为声波波速。
由式（4）可得：
$$ \frac{{{\partial ^2}u}}{{\partial {t^2}}} = \frac{\partial }{{\partial m}}\left( {{{\textit{z}}^2}\frac{{\partial u}}{{\partial m}}} \right) $$ (5) $$ \frac{{{\partial ^2}F}}{{\partial {t^2}}} = {{\textit{z}}^2}\frac{{{\partial ^2}F}}{{\partial {m^2}}} $$ (6) 由上述公式可知，声波在钻杆中的衰减特征与其频率、阻带、带宽等因素密切相关。式（6）即为质量坐标系下的经典波动方程，是理论分析纵波周期性钻柱结构中的传播特性的基础。
1.2 钻柱中声波信道的频率传输
在零初始条件下，用Fourier变换求解质量坐标系下的经典波动方程可得[13]：
$$ \frac{{{{\mathrm{d}}^2}f}}{{{\mathrm{d}}{m^2}}} + {k^2}f = 0 $$ (7) 式中：$ f(m,\omega ) $为$ F(m,t) $的傅立叶变换；$ \omega $为角频率。
令$ K = \omega /z $，$ K(m,\omega ) $是$ \omega $的周期函数K，用$ {K_\xi } $表示为：
$$ {K_\xi } = \frac{\omega }{{{{\textit{z}}_\xi }}} = \frac{\omega }{{{\rho _\xi }{S_\xi }{c_\xi }}} $$ (8) 利用傅里叶定理，将式（7）的通解表示为：
$$ f = f(m,\omega ){{\mathrm{e}}^{ - {\mathrm{j}}Km}} $$ (9) 周期性钻杆系统由2个部分单元组成，将每一部分的基本特性视为是1组常数，用下标$ \xi $标记。将各部分长度记为$ {L_\xi } $，质量记为$ {m_\xi } $，则：
$$ {m_\xi } = {\rho _\xi }{S_\xi }{L_\xi } $$ (10) 在给定的长度L上，方程式（7）的解也可表示为：
$$ {f_L} = {A_L}{{\mathrm{e}}^{ - {\mathrm{j}}KSm}} + {B_L}{{\mathrm{e}}^{{\mathrm{j}}KSm}} $$ (11) 将式（11）代入到式（9），可得：
$$ {f_L} = {f_L}(m,\omega ){{\mathrm{e}}^{( - {\mathrm{j}}knL/m)}} $$ (12) 式中：$ {f_L}(m,\omega ) = {A_L}{{\mathrm{e}}^{{\mathrm{j}}( {kL/{m{ { - {K_\xi }} }}} )m}} + {B_L}{{\mathrm{e}}^{{\mathrm{j}}( {kL/{m{{ + K{K_\xi }} )m}}}}} $。
由式（12）计算可得s段质点的振速为：
$$ {v_L} = \left( {1/{z_\xi }} \right){g_L}(m,\omega ){{\mathrm{e}}^{( - {\mathrm{j}}kmL/m)}} $$ (13) 式中：$ {g_L}(m,\omega ) = {A_L}{{\mathrm{e}}^{{\mathrm{j}}\left( {kL//{K_s}} \right)m}} - {B_L}{{\mathrm{e}}^{{\mathrm{j}}\left( {kL//m + {K_\xi }} \right)m}} $。
在周期性钻杆系统x=0和x=Li 2个接头处，作用力和质点振速满足均是连续的边界条件[14]，结合式（9）和式（13），应用Floquet理论，可得以下方程组：
$$ {f_{i - 1}}(0,\omega ) = {f_i}(0,\omega ) $$ (14-a) $$ \left( {1/{{\textit{z}}_2}} \right){g_{i - 1}}(0,\omega ) = \left( {1/{{\textit{z}}_1}} \right){g_i}(0,\omega ) $$ (14-b) $$ {f_{i - 1}}\left( { - {m_2},\omega } \right) = {f_i}\left( {{m_1},\omega } \right) $$ (14-c) $$ \left( {1/{{\textit{z}}_2}} \right){g_{i - 1}}\left( { - {m_2},\omega } \right) = \left( {1/{{\textit{z}}_1}} \right){g_i}\left( {{m_1},\omega } \right) $$ (14-d) 再将上面的方程组写成矩阵形式，即：
$$ \left( {\begin{array}{*{20}{c}} {{{\textit{z}}_1}}&{{{\textit{z}}_1}}&{{{\textit{z}}_2}}&{{{\textit{z}}_2}} \\ 1&{ - 1}&1&{ - 1} \\ {{{\textit{z}}_1}{{\mathrm{e}}^{{\alpha _1}{y_1}}}}&{{{\textit{z}}_1}{{\mathrm{e}}^{\beta {\gamma _1}}}}&{{{\textit{z}}_2}{{\mathrm{e}}^{ - {\alpha _2}{\gamma _2}}}}&{{{\textit{z}}_2}{{\mathrm{e}}^{ - {\beta _2}{y_2}}}} \\ {{{\mathrm{e}}^{{\alpha _1}{\gamma _1}}}}&{ - {{\mathrm{e}}^{{\beta _1}{\gamma _1}}}}&{{{\mathrm{e}}^{ - {\alpha _2}{y_2}}}}&{ - {{\mathrm{e}}^{ - {\beta _y}{y_2}}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{A_i}/{{\textit{z}}_1}} \\ {{B_i}/{{\textit{z}}_1}} \\ { - {A_{i - 1}}/{{\textit{z}}_2}} \\ { - {B_{i - 1}}/{{\textit{z}}_2}} \end{array}} \right) = \left( {\begin{array}{*{20}{l}} 0 \\ 0 \\ 0 \\ 0 \end{array}} \right) $$ (15) 式中：$ {\alpha _\xi } = {\mathrm{j}}\left( {kL/m - {K_\xi }} \right) $，$ {\beta _\xi } = {\mathrm{j}}\left( {kL/m - {K_\xi }} \right) $。
令式（15）=0，可得：
$$ \begin{split} & \begin{aligned} \cos\; kL = \cos \left( {\frac{{\omega {L_1}}}{{{c_1}}}} \right)\cos \left( {\frac{{\omega {L_2}}}{{{c_2}}}} \right) -\\ \frac{1}{2}\left( {\frac{{{{\textit{z}}_1}}}{{{{\textit{z}}_2}}} + \frac{{{{\textit{z}}_2}}}{{{{\textit{z}}_1}}}} \right)\sin \left( {\frac{{\omega {L_1}}}{{{c_1}}}} \right)\sin \left( {\frac{{\omega {L_2}}}{{{c_2}}}} \right) \end{aligned}\\[-14pt]& \end{split} $$ (16) 式中：c1、c2为阻抗z1、z2对应的声波波速。
当${\textit{z}}_1 \neq {\textit{z}}_2 $时，式（16）的解为非常数，即纵波声速与频率有关，发生了色散，此时方程称为色散方程。当$ |\cos\; kL| \lt 1 $时，方程有解，波数k为实数；反之，方程无解，波数k为复数。在波动方程中，当波数k为实数时，则其相应的频率信号相对无衰减通过，形成1条通带；当波数k为复数时，相应的频率由于衰减较大，形成了阻带，证实了声波在钻杆柱中的传播过程中，频率呈现出通带与阻带相间隔交替的梳状滤波器特征。
1.3 声波传播的衰减特性
声波沿钻杆传播的过程中，其信号强度随时间而衰减，主要受到接箍、声波频率及环境介质的阻尼等因素的影响[15]。声波在钻杆中传播时，其信号幅值随距离增加呈指数衰减。声波传输相对应的波动方程为：
$$ A = {A_0}{{\mathrm{e}}^{ - \alpha x}}{{\mathrm{e}}^{{\mathrm{j}}(\omega t - kx)}} $$ (17) 式中：A为声波在传播距离x处的波振幅；$ {A_0}{{\mathrm{e}}^{ - \alpha x}} $为波动方程中的包络线函数；A0为声波初始波振幅；α为声波的衰减系数。
声波的衰减系数可表示为：
$$ \alpha = \frac{1}{x}\ln \frac{{{A_0}}}{A} = \frac{{2.302\;6}}{x}\lg \frac{{{A_0}}}{A} $$ (18) 设声波信号的传播距离为xi，不同距离处的波振幅设为Ai，则衰减系数的一般表达式为：
$$ \alpha = \frac{1}{{\left( {{x_{i + 1}} - {x_i}} \right)}}\ln \frac{{{A_i}}}{{{A_{i + 1}}}}\quad i = 1,2,3, \cdots n$$ (19) 2. 钻杆声波传输有限元分析
钻杆连接示意如图3所示，钻杆两端通过螺纹与信号源（敲击端）以及公接头堵头（接收端）紧密扣接，振动传感器作为声波信号接收器安装于公接头堵头的外侧端面中心位置。
2.1 求解条件设定
采用通用隐式静力学求解钻杆重力载荷加载和接触状态，基于模态叠加法的模态动力学求解声波传递过程。计算时假设所有连接螺纹配合足够紧密的，不存在松动的可能，主要考虑钻杆自然搁置在水平地面/试验台上的情况。提取频率前定义钻杆与地面在重力载荷下的接触关系并计算其接触状态，获得真实的接触刚度。
钻杆材质为钢材，密度为7 850 kg/m3，弹性模量为210 GPa，泊松比为0.3。由于钻杆属于小质量构件，选用对高频信号较为敏感的beta材料刚度阻尼，阻尼系数由各个钻杆装配体的最低阶频率值和临界阻尼系数得出，临界阻尼系数取0.05。
为避免沙漏现象对声波传递能量的影响，单元类型选择全积分单元进行求解。将地面简化为解析刚体，对其参考点施加完全约束，并约束钻杆顶部中间节点的侧向自由度，避免钻杆在侧向产生不必要的滚动。实际计算时，将实测的加速度信号作为基础激励加载到敲击点，计算钻杆结构的实际响应过程。
2.2 不同因素的影响
为研究声波在钻杆中传播的影响因素，以钻杆有限元模型为基础，研究不同钻杆直径、不同敲击位置、不同放置形式以及埋沙注水2种工况下对声波传播特性和衰减规律的影响。选用最大幅值（对接收信号的正负最大值取绝对值后再平均）和时均幅值（接收信号围绕时间轴的面积除以信号接收时间0.05 s）2个指标来评估信号接收效果，以便对比不同因素作用下接收到的信号强度变化规律。
1）不同敲击位置的影响。以原始有限元模型为基础，改变敲击位置（竖直方向的最高点（a点）、敲击端面几何中心、竖直方向的最低点，得到的接收信号情况如图4所示。由图4可知，敲击位置从上到下，钻杆在重力的作用下，截面刚度沿着竖直方向逐渐增加，即越靠近接地侧截面刚度越大，越有利于信号传输，接收信号的最大幅值和时均幅值均逐渐增大，信号衰减变小。
2）钻杆直径的影响。增大钻杆直径分别进行计算，钻杆直径每次增加10 mm，即73.5 mm和83.5 mm。得到的不同钻杆直径下接收信号强度数据见表1。由表1可以看出，钻杆直径越大，钻杆的接收信号的时均幅值明显增大，信号衰减变小，当直径为83.5 mm时，接收信号的时均幅值为0.080。当钻杆直径越大时，在竖向重力作用下，钻杆的接地面积和接地刚度均增大，且直径增大的同时保持钻杆厚度不变，必然造成钻杆整体材料体积增加，材料的材料阻尼耗能增加，当阻尼耗能的影响幅度小于接地刚度的贡献幅度时，信号的整体衰减降低，表现出钻杆直径增大有利于信号传输的特性。
表 1 不同钻杆直径的信号强度数据Table 1. Signal strength data for different drill pipe diameters钻杆直径/mm 最大幅值 时均幅值 63.5 0.89 0.071 73.5 0.94 0.076 83.5 0.94 0.080 3）钻杆不同放置形式的影响。将钻杆竖直放置，钻杆上部为悬挂端，下部为敲击端，重力方向竖向向下。对比钻杆水平放置和竖直放置条件下的声波接收信号可得，竖直放置条件下声波信号强度低于水平放置，信号时均幅值降低了32%。在重力作用下，钻杆属于大柔度结构且轴向主要承受预拉力，钻杆整体刚度降低，由于声波在刚度较大的介质中传播距离更远，因此竖直放置条件下不利于声波信号的传输。
4）外部埋沙对声波传输的影响。将埋沙环境简化为一定大小的竖向载荷，并以体力的方式叠加在钻杆重力载荷上，总体力载荷大小为10倍钻杆重力。计算得到的接收信号强度数据表明，在外部埋沙条件下，受外部埋沙竖向压力载荷的影响，钻杆与地面接触压力增大，接触刚度有所增加，信号的最大幅值略有增加，但幅度较小，时均幅值大小与无埋沙情况基本一致，接收信号整体与无埋沙情况一致，故外部埋沙条件对钻杆信号接收影响很小，其影响基本可以忽略。
5）内部注水条件对声波传输的影响。针对钻杆内部注水的工作条件，计算时将内部水域简化为仅抗压的材料，模拟水的流体特性，计算得到的钻杆内部有无注水条件下的声波接收信号对比情况如图5所示。由图5可以看出，在内部注水条件下，接收端信号明显有所衰减，说明钻杆内部注水条件不利于信号传输。内部注水首先会增大钻杆的整体质量，使得钻杆和地面的接触压力增大，接触刚度增大，有利于信号的传输，但是内部注水同时也会增大水的黏性耗能，部分信号进入水域后必然造成信号能量的衰减。刚度增加的贡献小于水黏性耗能的贡献时，黏性耗能占主导因素，最终造成信号的整体衰减。
![]()
图 5 内部是否注水条件下接收信号对比Figure 5. Comparison of received signals with and without internal water injection condition2.3 多根钻杆连接分析
简化钻柱有限元模型，侧重声波在钻柱远距离传输中的衰减特性研究。针对不同钻杆根数计算得到钻柱整体频率变化规律，随着连接钻杆数量的增加，整体钻柱的质量增大，钻杆与地面在连接头位置的接触点数增多，钻杆侧向刚度相对变小，故整体频率值降低，自振周期变长。接收信号强度随传输距离的衰减曲线如图6所示。
![]()
图 6 接收信号强度随传输距离的衰减曲线Figure 6. Attenuation curves of received signal strength with transmission distance由图6可以看出，在较短距离（30 m以内）传输时，接收信号的最大幅值衰减较快，在较远距离（60 m以外）传输时，接收信号的最大幅值衰减幅度变缓，信号强度保持较好。在最大计算距离（300 m）时，接收信号的强度仍不低于60%，衰减倍数约为1.6倍。
考虑到接收信号幅值存在正负，且存在一定的波动，将接收信号的正负最大值取绝对值并平均作为信号强度指标SI（$0 \leqslant {\mathrm{SI}} \leqslant 1$）[16]。对计算传输距离内的最大接收信号强度进行拟合，得到信号强度SI随传输距离d的衰减关系式为：
$$ {\text{SI}} = 0.988{{{d}}^{0.094}} $$ 由上式可以看出，衰减服从典型的幂函数衰减规律。与试验结果相比，仿真分析得到信号的衰减幅度略小，这是由于试验中可能存在接头扣紧力矩施加不足引起的局部松动，试验信号在这些位置的传输时衰减增大，而仿真模型中各钻杆的连接较试验连接更为理想，信号衰减较小。另外，试验条件下，信号总能量的一部分会以声波等形式传递到试验环境中，造成接收信号强度的进一步损失。
3. 钻杆声波传播试验
为研究声波在钻杆中的传播和衰减规律，搭建试验平台进行测试，试验平台结构示意图如图7所示。平台由声波信号发生设备、简易钻杆系统（不同数量的ϕ73 mm×1.5 m常规外平钻杆通过螺纹相连接）和信号接收模块3个部分组成。
声信号发生设备由驱动电源、信号发生器和换能器组成：驱动电源发出高压脉冲驱动信号；信号发生器按照固定频率发出控制信号；换能器将高压脉冲驱动信号转化为脉冲振动信号纵向敲击钻杆公接头端面。在较大的量程宽度内，选取了具有良好低频响应性能的压电式加速度传感器用于信号接收模块，利用示波器采集声波传感器信号进行输出显示和保存。
通过所搭建的测试平台，在不改变信号源敲击强度的条件下，以5 Hz的频率进行连续敲击，并逐渐增加钻杆柱的长度，记录接收到的原始波形。结果表明，4.5 m和289.5 m处所接收到的原始波形都较为清晰，信噪比相近，而信号的背景噪声稍大，毛刺较多。随后，为了降低信噪比、滤除背景噪声，截取峰值段波形进行傅里叶变换得到频域波形，频域波形如图8所示。
由图8可以看出，峰值段信号频域成分主要分布在1~3 kHz之间，因此对原始信号1~3 kHz以外的信号进行滤波，原始波形和滤波后波形叠加图如图9所示。
![]()
图 9 滤波前后波形叠加图（红色为滤波后波形）Figure 9. Superposition of waveforms before and after filtering (red is the waveform after filtering)从滤波后得到的信号波形可以看出，背景噪声得到大幅消除，信号峰值基本没有衰减，信噪比得到提高，289.5 m处接收信号滤波效果更加明显。随后试验过程中，增加钻杆柱连接长度测试接收波形，每增加4根钻杆，即6 m钻杆柱，测试1组接收波形，统计所得的钻杆柱接收信号幅值变化与仿真分析结果一致，呈指数衰减趋势。
4. 结 语
1）以质量坐标系下的经典波动方程作为理论分析纵波周期性钻柱结构中的传播特性的基础。声波的频率传输过程中，其频带呈现出通、阻带交替的梳状滤波器结构特性，声波信号幅值近似按照指数规律衰减。
2）通过有限元分析钻杆声波传输过程得到：钻杆外直径越大、信号敲击端的敲击点位置距离地面越近，越有利于信号的传输，信号衰减少；钻杆内部注水条件下，接收信号明显衰减；而外部埋沙环境对信号强度的影响可以忽略；当多根钻杆连接时，接头部位钻杆与地面的接触点增多，钻杆侧向刚度相对变小，导致整体频率下降，自振周期变长。
3）钻杆柱中纵向传播的声波波形比较稳定，适合进行声波传输。在基本没有外部噪声的环境中传播时，短距离（4.5 m）和长距离（289.5 m）接收到的波形轮廓非常相似；接收信号峰值部分的频域成分集中在1个较窄的频域范围内，对此频域范围以外的频域进行滤波，能有效提高信噪比，获得良好的通频滤波效果。
-
![]()
图 3 模拟平衡后水的密度场分布图
Figure 3. Density field distribution of water after simulated equilibrium
![]()
图 4 CTAB头基和H2O的径向分布曲线
Figure 4. Radial distribution curves of CTAB head groups and H2O
![]()
图 5 CTAB尾基和H2O的径向分布曲线
Figure 5. Radial distribution curves of CTAB tail groups and H2O
![]()
图 9 各个分子或晶胞静电势图
Figure 9. Diagrams of the electrostatic potential of individual molecules or crystals
表 1 盒子或超晶胞参数
Table 1 Box or supercell parameters
名称 长/nm 宽/nm 高/nm 分子数量 α、β、γ/（°） 密度/ （g·cm−3） 水盒子 2.946 0 3.241 2 1.566 4 500 90 1.000 CTAB盒子 2.946 0 3.241 2 0.239 7 5 90 1.322 煤分子盒子 2.946 0 3.241 2 1.342 6 8 90 1.350 SiO2超晶胞 2.946 0 3.241 2 2.288 3 — 90 2.200 注：α、β、γ为三维空间中与x、y、z轴的夹角。 
下载: 导出CSV
表 2 水−CTAB−煤尘/矽尘的相互作用能
Table 2 Interaction energies of water-CTAB-coal dust/silicon dust
体系 能量/（kcal·mol−1） 总能量 水和CTAB 粉尘 相互作用能 水−CTAB−煤尘 −2 739.640 2 −4 053.856 8 2 107.936 5 −793.720 4 水−CTAB−矽尘 −24 706.914 0 −4 065.075 7 −20 451.980 0 −135.598 0 注：1 kcal=4.184 kJ。
下载: 导出CSV
表 3 各个分子或晶胞的静电势对比表
Table 3 Comparison of the electrostatic potential of individual molecules or crystals
分子或晶胞种类 静电势/（Ha·e−1） 正 负 尾链 H2O 0.091 56 −0.065 44 — CTAB 0.203 40 — 0.021 89 WD 0.124 20 −0.070 96 — SiO2 0.568 80 −0.020 57 —
下载: 导出CSV
-
[1] CHEN Z R, QIU H X, HONG Z X, et al. Molecular dynamics simulation of a lignite structure simplified model absorbing water[J]. Molecular Simulation, 2020, 46(1): 71−81. doi: 10.1080/08927022.2019.1674851
[2] QIU H X, CHEN Z R, WANG G H. Effect of dodecyl trimethyl ammonium bromide on the migration of water molecules in the pores of lignite: an experimental and molecular simulation study[J]. ACS Omega, 2020, 5(39): 25456−25466. doi: 10.1021/acsomega.0c04012
[3] 高正阳,杨维结. 不同煤阶煤分子表面吸附水分子的机理[J]. 煤炭学报,2017,42(3):753−759. GAO Zhengyang, YANG Weijie. Adsorption mechanism of water molecule on different rank coals molecular surface[J]. Journal of China Coal Society, 2017, 42(3): 753−759.
[4] XU L, PEI Z H. Preparation and optimization of a novel dust suppressant for construction sites[J]. Journal of Materials in Civil Engineering, 2017, 29(8).
[5] HUANG Z A, ZHANG L H, YANG Z J, et al. Preparation and properties of a rock dust suppressant for a copper mine[J]. Atmospheric Pollution Research, 2019, 10(6): 2010−2017. doi: 10.1016/j.apr.2019.09.008
[6] CRAWFORD R J, MAINWARING D E. The influence of surfactant adsorption on the surface characterisation of Australian coals[J]. Fuel, 2001, 80(3): 313−320. doi: 10.1016/S0016-2361(00)00110-1
[7] ZHANG L, LI B, XIA Y C, et al. Wettability modification of Wender lignite by adsorption of dodecyl poly ethoxylated surfactants with different degree of ethoxylation: A molecular dynamics simulation study[J]. Journal of Molecular Graphics and Modelling, 2017, 76: 106−117. doi: 10.1016/j.jmgm.2017.06.028
[8] 陈璐瑶. 表面活性剂对二氧化硅粉尘润湿性影响的分子动力学研究[D]. 兰州:兰州大学,2021. [9] 孟筠青,王丽娟,王顾儒,等. 表面活性剂对硅尘润湿性影响的分子模拟[J]. 安全与环境学报,2023,23(7):2360−2365. MENG Junqing, WANG Lijuan, WANG Guru, et al. Molecular simulation of the influence of surfactants on the wettability of silica dust[J]. Journal of Safety and Environment, 2023, 23(7): 2360−2365.
[10] 邱玄,钱玉鹏,陈彬,等. Pluronic F-68在萤石和石英表面的吸附特性及分散作用研究[J]. 金属矿山,2021(12):56−61. QIU Xuan, QIAN Yupeng, CHEN Bin, et al. Study on adsorption characteristics and dispersion of pluronic F-68 on the surface of fluorite and quartz[J]. Metal Mine, 2021(12): 56−61.
[11] WANG Z Q, WANG P Y, SONG H, et al. Dynamic wetting behavior of nanofluid droplet on a vertically vibrating surface: A molecular dynamics study[J]. Journal of Molecular Liquids, 2022, 347: 118360. doi: 10.1016/j.molliq.2021.118360
[12] ZHOU G, XING M Y, WANG K L, et al. Study on wetting behavior between CTAC and BS-12 with gas coal based on molecular dynamics simulation[J]. Journal of Molecular Liquids, 2022, 357: 118996. doi: 10.1016/j.molliq.2022.118996
[13] DELLEY B. From molecules to solids with the DMol3 approach[J]. 2000, 113(18): 7756-7764.
[14] GRIMME S. Accurate description of van der Waals complexes by density functional theory including empirical corrections[J]. Journal of Computational Chemistry, 2004, 25(12): 1463−1473. doi: 10.1002/jcc.20078
[15] LIU L, JIN J, LIN Y Y, et al. The effect of calcium on nitric oxide heterogeneous adsorption on carbon: A first-principles study[J]. Energy, 2016, 106: 212−220. doi: 10.1016/j.energy.2016.02.148
[16] GRIMME S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction[J]. Journal of Computational Chemistry, 2006, 27(15): 1787−1799. doi: 10.1002/jcc.20495
[17] LI W F, WANG H N, LI X, et al. Effect of mixed cationic/anionic surfactants on the low-rank coal wettability by an experimental and molecular dynamics simulation[J]. Fuel, 2021, 289: 119886. doi: 10.1016/j.fuel.2020.119886
[18] XU J F, ZHANG Y, CHEN H X, et al. Effect of surfactant headgroups on the oil/water interface: an interfacial tension measurement and simulation study[J]. Journal of Molecular Structure, 2013, 1052: 50−56. doi: 10.1016/j.molstruc.2013.07.049
[19] 朱蓓蓓. 溴化锂水溶液微观特性的分子动力学研究[D]. 大连:大连海事大学,2010. [20] 韩方伟,胡福宏,奚志林,等. AEO3润湿褐煤性能及机理[J]. 煤炭学报,2023,48(10):3766−3775. HAN Fangwei, HU Fuhong, XI Zhilin, et al. Study on wettability and mechanism of AEO3 to lignite[J]. Journal of China Coal Society, 2023, 48(10): 3766−3775.
计量
- 文章访问数: 36
- HTML全文浏览量: 2
- PDF下载量: 1
