[發明專利]用于煤層氣的交互型解吸曲率變化的表征方法

有效

• 申請號: 201910161837.7
• 申請日: 2019-03-04
• 公開（公告）號: CN110059343B
• 公開（公告）日: 2022-11-25
• 發明（設計）人: 張學梅;李東;郝靜遠;馬青華
• 申請（專利權）人: 西安思源學院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F16/21
• 代理公司: 暫無信息
• 代理人: 暫無信息

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曲線進行階段劃分，以獲取用于表征煤層氣的交互型解吸曲率變化的系列曲線。

2.根據權利要求1所述的用于煤層氣的交互型解吸曲率變化的表征方法，其特征在于，在所述步驟S6中，用polyfit()命令對KL'和KL進行多項式擬合的次數為分別為10次。

3.根據權利要求1所述的用于煤層氣的交互型解吸曲率變化的表征方法，其特征在于，在所述步驟S6中，用polyfit()命令對KL'和KL進行多項式擬合的次數為分別5次和7次。

4.根據權利要求1所述的用于煤層氣的交互型解吸曲率變化的表征方法，其特征在于，在所述步驟S8中，具體包括：分別調用所述V數據庫、所述V'數據庫、所述KL數據庫、所述KL'數據庫以及所述KL數據庫，以繪制p-V曲線、p-V'曲線、p-KL曲線、p-KL'曲線及p-KL曲線，其中p-V'曲線即為p-η曲線，η為煤層氣解吸效率。