一项关于元素硫在硫化氢中的溶解度的分子模拟研究，该硫化氢被限制在方解石纳米孔内

《Fluid Phase Equilibria》：A molecular simulation study on solubility of elemental sulfur in hydrogen sulfide confined in calcite nanopores

硫在受限碳酸盐纳米孔隙中的溶解度演变及分子机制研究，通过全原子分子动力学模拟发现孔隙越小硫溶解度越高，表面吸附效应和分子聚类结构变化是主因，为高硫气藏开发提供理论依据。

　　
该研究聚焦于高硫气田开发中的关键科学问题——受限纳米孔隙内硫溶度演变规律及其作用机制。通过分子动力学模拟方法，系统揭示了方解石纳米孔隙中硫溶度受孔隙尺寸、温度与压力三重调控的规律，并首次阐明了纳米受限环境下硫溶度异常增大的微观机理。研究结果为优化高硫气田开发技术提供了重要理论支撑。

在气田开发过程中，硫的沉积会引发孔隙 throat 堵塞，造成气体渗透率下降和安全生产风险。现有研究表明，纳米孔隙的受限效应会显著改变溶质在介质中的溶解行为，但针对硫分子在碳酸盐纳米孔隙中的溶解特性研究存在空白。该团队突破传统模拟方法的局限性，采用改进的直接关联法（DC）模拟技术，有效解决了硫分子结晶过程的长周期模拟难题。

研究构建了具有真实晶格结构的方解石纳米孔隙模型，其晶面参数与实验测定值完全吻合。通过调控纳米孔隙的开口尺寸（1-5 nm），在373.2 K（100℃）和30 MPa典型储层条件下，发现硫溶度呈现显著的空间受限效应。当孔隙尺寸降至1 nm时，硫溶度较无约束的 bulk 系统提升达3倍。这种反常增溶现象源于纳米孔隙特有的三维限域效应：在孔隙表面形成的硫吸附层通过氢键和范德华力产生协同作用，促使更多硫分子进入溶剂相。

分子动力学模拟揭示了多级作用机制：首先，方解石（1014）晶面的负电荷表面通过静电吸引作用富集硫分子；其次，孔隙内壁的羟基基团与硫分子形成氢键网络，降低硫分子间的聚集势垒；再者，纳米孔隙的几何约束抑制了硫簇的晶格重构，使单个S8分子保持高溶解活性。这种"吸附驱动-空间限制"的复合效应，使得纳米孔隙内的硫溶解度突破传统热力学平衡预期。

温度与压力的协同调控作用呈现出非线性特征。当储层埋深增加导致压力升高时，孔隙内H2S溶剂分子密度增大，形成更稳定的硫-溶剂复合体。温度升高则通过增强分子动能促进硫分子解离，但这一过程在纳米尺度下因空间受限而呈现延迟效应。模拟数据显示，在3 nm孔隙中，50℃时的硫溶解度较25℃提升37%，而相同温度下压力每增加10 MPa，溶解度增幅可达15%。

研究创新性地提出"硫溶度放大倍数"概念，量化了纳米孔隙对硫溶解行为的增强效应。通过追踪硫分子在孔隙中的动态行为，发现其存在两种主要存在形式：一种是分散于溶剂相的亚稳态单分子硫；另一种是受孔隙几何约束难以结晶的微小硫簇（平均直径<0.5 nm）。当孔隙宽度减小时，单分子硫比例从 bulk 系统的12%提升至5 nm孔隙的38%，而硫簇尺寸的分布标准差降低至0.21 nm，表明孔隙限制有效抑制了硫分子聚集。

该成果在气田开发领域具有重要指导价值：其一，通过孔隙尺寸调控储层渗透性，选择性的纳米孔隙发育区域可作为硫沉积的"安全带"；其二，在气井完井作业中，采用纳米孔隙特征材料进行选择性封堵，可有效阻止硫沉积向主渗流通道扩散；其三，揭示的硫分子解离动力学参数，为优化驱油剂配方提供了新依据。研究建立的纳米孔隙硫溶度预测模型，可拓展应用于其他硫化合物在碳酸盐岩储层中的相态行为预测。

当前研究仍存在待深化方向：首先，对非立方晶系碳酸盐岩的孔隙特性研究不足；其次，复杂储层条件下多组分耦合作用机制尚未完全明晰；再者，实验验证中如何准确表征纳米孔隙内的局部相态仍是技术难点。未来研究可结合原位表征技术，在微观尺度建立硫溶解动力学模型，同时拓展至其他硫化物在复杂储层环境中的相态行为研究。

该研究不仅完善了纳米限域效应下多相化学反应的理论体系，更直接回应了"如何量化纳米孔隙硫溶度提升倍数"这一工程界亟待解决的核心问题。通过建立孔隙尺寸-温度-压力与硫溶度的量化关系模型，为制定储层改造方案提供了关键参数。例如，在3-5 nm的临界孔隙尺寸范围内，硫溶解度增幅达300%，这为确定储层改造的纳米孔隙阈值提供了理论支撑。研究揭示的硫分子吸附解离路径，对开发基于分子识别原理的硫阻垢剂具有重要启示。