[发明专利]二氧化硅表面官能团修饰任意可控且位置随机的建模方法

说明书

本发明涉及二氧化硅表面官能团修饰任意可控且位置随机的建模方法，包括：导入二氧化硅的原始晶胞结构模型；建立二氧化硅表面结构；通过原始晶胞表面结构XY二维方向上进行5×5超胞化得到表面为100个氧原子的超晶胞结构；定义表面氧原子组，编程得到表面各氧原子信息；设置变量N用于确定表面生成N个甲基官能团；提取乱序数组中的氧原子，得到被提取原子的XYZ坐标信息；在被提取原子上方生成一个甲基官能团的C原子，Z方向上与被提取原子距离为一个键长距离；甲基官能团的C原子与原表面氧原子成键，对其加氢处理为‑CH₃官能团，对其余氧原子加氢处理成为‑OH官能团。本发明通过准确控制CH₃官能团修饰表面结构，构建油气藏储层矿物不同润湿性微观模型。

技术领域

本发明涉及石油与天然气行业微观模拟技术领域，具体涉及矿物表面不同润湿性官能团修饰的建立方法。

背景技术

润湿性是指不相混的两相流体与岩石固相表面接触时，其中某一相流体沿着岩石表面铺开的现象，其结果使体系的表面自由能降低，沿着岩石表面铺开的一相称为润湿相[杨胜来, 魏俊之. 油层物理学[M], 石油工业出版社]。润湿性是油气藏储层岩石的重要基础参数之一，其特征直接影响流体在岩石孔道内的微观分布和宏观分布特征，也制约着油气开发效率。因此，在油气开发技术中，往往添加表面活性剂来改变矿物表面的润湿性。

表面活性剂是兼具疏水基和亲水基的两亲型化合物，通过让表面活性剂流入及扩散进入基质，可改变储集岩的润湿性，同时降低界面张力减小毛管压力。表面活性剂改变润湿性的主要机理为：离子对形成、表面活性剂吸附和胶束增溶作用，研究表明，阳离子表面活性剂、非离子表面活性剂、两性离子表面活性剂和阴离子表面活性剂具有较好效果，其中，非离子表面活性剂主要含甲基或羟基[Nguyen D, Wang D, Oladapo A, et al.Evaluation of Surfactants for Oil Recovery Potential in Shale Reservoirs[C]//SPE Improved Oil Recovery Symposium. Society of Petroleum Engineers, 2014]。此外，在一定温度压力条件下，二甲基二氯硅烷可用于纳米二氧化硅表面改性，随着二甲基二氯硅烷用量的增加，纳米二氧化硅表面羟基数下降[唐洪波,李萌,马冰洁.二氯二甲基硅烷改性纳米二氧化硅工艺研究[J].精细石油化工,2007(06):44-47]。因此，从微观角度而言，可通过不同比例的甲基与羟基来表征不同润湿性矿物。

分子模拟方法从微观状态分布来计算热力学参数，可体现微纳米尺度界面效应的影响。近年来，分子模拟被广泛运用于研究油气领域中流体的界面效应及微纳米孔隙内的界面特性，研究者们采用分子模拟建立矿物结构并进行了不同表面润湿性的模拟研究，考察了不同晶面、亲水和亲油官能团表征不同润湿性石英吸附甲烷的影响[汪周华,赵建飞,白银,郭平,刘煌.不同润湿性修饰石英吸附甲烷的模拟研究[J].西南石油大学学报(自然科学版),2019,41(06):28-34]，但其亲润性模型为特定官能团对表面全覆盖修饰。研究人员对不同润湿性进行了详细探究[焦红岩,董明哲,刘仲伟,靳彦欣,韩旭.水环境下甲烷在不同润湿性石英表面吸附行为的分子动力学模拟[J].中国石油大学学报(自然科学版),2014,38(05):178-183]，构建了几种特定比例不同润湿性的石英表面模型，石英表面修饰程度n(-CH₃): n(-OH)为0:100、25:75、50:50、75:25、100:0，但其所研究模型表面修饰比例特定，且官能团分布规则，未体现任意比例以及官能团出现位置的随机性。

到目前为止，在油气储层矿物研究背景下，不同润湿性二氧化硅界面及其亲润性的模拟计算研究过程中，研究人员虽然建立了不同二氧化硅表面结构模型，但由于二氧化硅表面氧原子结构的特殊性，这些模型针对不同润湿性修饰而言均有不足之处，主要表现为：特定官能团全覆盖修饰表面、CH₃/OH特定比例修饰表面、官能团分布规则无随机性，无法充分考虑和有效控制模型中的不同亲润性官能团比例。

发明内容