用于煤层气的交互型解吸曲率变化的表征方法技术

本发明专利技术提供一种用于煤层气的交互型解吸曲率变化的表征方法，包括：建立原始实验数据库；将压力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曲线进行阶段划分。与相关技术相比，本发明专利技术提供的用于煤层气的交互型解吸曲率变化的表征方法易学习、易使用、高效且准确性高。

Characterization of Interactive Desorption Curvature Change for Coalbed Methane

【技术实现步骤摘要】
用于煤层气的交互型解吸曲率变化的表征方法
本专利技术涉及交互型解吸曲率的表示方法领域，尤其涉及一种用于煤层气的交互型解吸曲率变化的表征方法。
技术介绍
煤层吸附气含气总量的80％以上，主要成分为甲烷，煤岩储层的吸附量和吸附特征是煤层气地质评价的重要依据。在煤层气研究中主要采用兰格缪尔(Langmuir)等温吸附理论表征其吸附行为，煤层气等温吸附与解吸理论上是完全可逆的。目前主要根据“排水—降压—解吸—扩散—渗流”的机理对煤层气进行开发利用。因此，煤层气的吸附与解吸行为特征对其开发有重大影响。大量煤层气研究主要关于等温吸附(解吸)实验及其主要的影响因素，对解吸过程理论研究较少，且没有很好的方法用于处理实验数据，很难实现对不同煤层气藏解吸过程进行对比分析及指导实际开发。根据等温吸附实验可知，不同煤层兰氏体积和兰氏压力不同，从而导致煤层产能不同，故在开发过程中不同煤层气藏的排采条件亦不同。煤层气井产能主要受渗透率和解吸效率的共同控制，若解吸效率小于煤层裂缝的渗流能力，即使有较高的渗透率，煤层气井产能依旧很低。因此，提高产能的主要途径是提高甲烷的解吸效率。影响解吸效率高低的关键因素有：兰氏体积、兰氏压力和煤储层压力，兰氏体积和兰氏压力越高，解吸效率越高；煤储层压力越高，解吸效率越低，反之亦然。故在埋深相似的条件下，高阶煤储层比低阶煤储层在形成高产能井更有优势；在相似的煤阶储层条件下，埋深不同的煤储层其产气能力截然不同。深入探究煤层气的解吸特征，提供一种可以表征任何煤阶储层煤层气的解吸特征至关重要，只有了解煤储层的解吸特征，结合压降漏斗才能实现对煤层气井产能的定量评价，为气井开发提供重要的依据。因此，有必要提供一种新的用于煤层气的交互型解吸曲率变化的表征方法以解决上述问题。
技术实现思路
本专利技术需要解决的技术问题是提供一种易学习、易使用、高效且准确性高的用于煤层气的交互型解吸曲率变化的表征方法。为解决上述问题，本专利技术提供了一种用于煤层气的交互型解吸曲率变化的表征方法，该方法包括如下步骤：步骤S1、收集langmuir等温吸附数据，建立原始实验数据库，所述langmuir等温吸附数据包括兰氏体积VL和兰氏压力PL；步骤S2、将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积VL和所述兰氏压力PL分别赋值给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曲线进行阶段划分，以获取用于表征煤层气的交互型解吸曲率变化的系列曲线。优选的，在所述步骤S6中，用polyfit()命令对KL'和KL"进行多项式拟合的次数分别为10次。优选的，在所述步骤S6中，用polyfit()命令对KL'和KL"进行多项式拟合的次数为分别5次和7次。优选的，在所述步骤S8中，具体包括：分别调用所述V数据库、所述V'数据库、所述KL数据库、所述KL'数据库以及所述KL"数据库，以绘制p-V曲线、p-V'曲线、p-KL曲线、p-KL'曲线及p-KL"曲线，其中p-V'曲线即为p-η曲线，η为煤层气解吸效率。与相关技术比较，本专利技术的用于煤层气的交互型解吸曲率变化的表征方法中，将交互设计用于表征煤层气的解吸曲率的变化，能快速得到高精度的计算结果，并将计算结果赋值给转折压力、敏感压力和启动压力实现自动对解吸曲线进行阶段划分，同时可以实现对任意储煤层是否存在转折压力、敏感压力和启动压力的判断，实现对任意煤阶储层煤层气吸附行为的表征，该方法简单、快速，且结果精确度高；对于零基础的人均可使用，操作简单，结果直观明了；结合降压漏斗可对储煤层的产气潜能和产气量进行定量分析；该方法为煤层气井排采工艺设计提供了便利。附图说明为了更清楚地说明本专利技术实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图，其中：图1为本专利技术的用于煤层气的交互型解吸曲率变化的表征方法的流程框图。图2为本专利技术的用于煤层气的交互型解吸曲率变化的表征方法实施例一的用于表征煤层气的交互型解吸曲率变化的系列曲线。其中，图2a为煤层气吸附曲线图；图2b为煤层气解吸效率曲线图；图2c为煤层气解吸曲率曲线图；图2d为煤层气解吸曲率一阶导数曲线图；图2e为煤层气解吸曲率二阶导数曲线图。图3为本专利技术的用于煤层气的交互型解吸曲率变化的表征方法实施例二的用于表征煤层气的交互型解吸曲率变化的系列曲线。其中，图3a为煤层气吸附曲线图；图3b为煤层气解吸效率曲线图；图3c为煤层气解吸曲率曲线图；图3d为煤层气解吸曲率一阶导数曲线图；图3e为煤层气解吸曲率二阶导数曲线图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本专利技术保护的范围。请参阅图1，本专利技术提供了一种用于煤层气的交互型解吸曲率变化的表征方法，以下结合具体实例进行说明：实施例一以山西毛则渠矿5号煤层瘦煤为例，该方法具体包括如下步骤：步骤S1、收集langmuir等温吸附数据，建立原始实验数据库，所述langmuir等温吸附数据包括兰氏体积VL和兰氏压力PL。本步骤中，调用所述原始实验数据库中30℃等温吸附数据的兰氏体积VL和兰氏压力PL，分别为VL＝20.60m3/t，PL＝3.46MPa。步骤S2、将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积VL和所述兰氏压力PL分别赋值给a、b。本步骤中的目的是为了实现高精度运算，即a＝20.60，b＝3.46。步骤S3、将a、b、p代入到langmuir吸附等温式中求解煤层气解吸量V，其中，压力p范围为0MPa≤p≤15MPa]，将计算结果存入V数据库。步骤S4、用MATLAB中diff()求导命令求解煤层气解吸量V的一阶导数V'和二阶导数V"，将结果分别存入V'数据库和V"数据库中。本步骤中，具体求导命令为diff(V)和di本文档来自技高网...

【技术保护点】
1.一种用于煤层气的交互型解吸曲率变化的表征方法，该方法包括如下步骤：步骤S1、收集langmuir等温吸附数据，建立原始实验数据库，所述langmuir等温吸附数据包括兰氏体积VL和兰氏压力PL；步骤S2、将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积VL和所述兰氏压力PL分别赋值给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曲线进行阶段划分，以获取用于表征煤层气的交互型解吸曲率变化的系列曲线。...

【技术特征摘要】
1.一种用于煤层气的交互型解吸曲率变化的表征方法，该方法包括如下步骤：步骤S1、收集langmuir等温吸附数据，建立原始实验数据库，所述langmuir等温吸附数据包括兰氏体积VL和兰氏压力PL；步骤S2、将压力p定义为符号变量，定义a、b为常数，并将所述兰氏体积VL和所述兰氏压力PL分别赋值给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()命令去...

【专利技术属性】
技术研发人员：张学梅，李东，郝静远，马青华，
申请(专利权)人：西安思源学院，
类型：发明
国别省市：陕西,61