[发明专利]用于煤层气的交互型解吸曲率变化的表征方法

摘要

本发明提供一种用于煤层气的交互型解吸曲率变化的表征方法，包括：建立原始实验数据库；将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积VL和所述兰氏压力PL分别赋值给a、b；求解煤层气解吸量V；求解煤层气解吸量V的一阶导数V'和二阶导数V"；求得煤层气解吸曲率KL和KL的一阶导数KL'与二阶导数KL"，分别计算所述KL、KL'以及KL"表达式的值；对KL'和KL"进行多次多项式拟合，求得KL'和KL"的根并用find()命令去除复数项；输出计算结果；绘制煤层气解吸特征系列曲线并对其中p‑V曲线进行阶段划分。与相关技术相比，本发明提供的用于煤层气的交互型解吸曲率变化的表征方法易学习、易使用、高效且准确性高。

主权项

1.一种用于煤层气的交互型解吸曲率变化的表征方法，该方法包括如下步骤：步骤S1、收集langmuir等温吸附数据，建立原始实验数据库，所述langmuir等温吸附数据包括兰氏体积V_L和兰氏压力P_L；步骤S2、将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积V_L和所述兰氏压力P_L分别赋值给a、b；步骤S3、将a、b、p代入到langmuir吸附等温式中求解煤层气解吸量V，其中，压力p范围为0MPa≤p≤15MPa]，将计算结果存入V数据库；步骤S4、用MATLAB中diff()求导命令求解煤层气解吸量V的一阶导数V'和二阶导数V"，将结果分别存入V'数据库和V"数据库中；步骤S5、将所述一阶导数V'和所述一阶导数V"代入曲率公式中求得煤层气解吸曲率KL，并用diff()命令求解所述煤层气解吸曲率KL的一阶导数KL'和二阶导数KL"，用eval()命令分别计算所述KL、KL'以及KL"表达式的值，将计算结果分别存入KL数据库、KL'数据库和KL"数据库中；步骤S6、用polyfit()命令对KL'和KL"进行多次多项式拟合，用roots()命令求得KL'和KL"的根并用find()命令去除复数项；步骤S7、输出计算结果；将所述步骤S6中输出的KL'的根赋值给ptu，命名为转折压力；KL"的实根中小值赋值给pse，大值赋值给pst，分别命名为敏感压力和启动压力；步骤S8、用subplot()和plot()命令绘制煤层气解吸特征系列曲线并对其中p‑V曲线进行阶段划分，以获取用于表征煤层气的交互型解吸曲率变化的系列曲线。