Quilghini变换法在求解密度跃变的水合物热分解Stefan模型精确解中的应用

doi:10.11949/0438-1157.20240509

化工学报 ›› 2024, Vol. 75 ›› Issue (10): 3763-3774.DOI: 10.11949/0438-1157.20240509

Quilghini变换法在求解密度跃变的水合物热分解Stefan模型精确解中的应用

李明川¹^,²(), 樊栓狮³, 徐赋海⁴, 卢惠东⁴, 李晓军⁴

^1.中国石油大学（华东）石油工程学院，山东青岛 266580
^2.非常规油气开发教育部重点实验室（中国石油大学（华东）），山东青岛 266580
^3.华南理工大学化学与化工学院，广东广州 510641
^4.中国石化胜利油田分公司东辛采油厂，山东东营 257094

收稿日期:2024-05-09 修回日期:2024-06-17 出版日期:2024-10-25 发布日期:2024-11-04
通讯作者: 李明川
作者简介:李明川（1976—），男，博士，副教授，iceswpi@126.com
基金资助:
中国科学院知识创新工程基金项目(KGCX2-SW-309)

Application of Quilghini transformation method to exact solutions of Stefan's models for thermal dissociation hydrate with density jump

Mingchuan LI¹^,²(), Shuanshi FAN³, Fuhai XU⁴, Huidong LU⁴, Xiaojun LI⁴

^1.School of Petroleum Engineering，China University of Petroleum, Qingdao 266580, Shandong, China
^2.Key Laboratory of Unconventional Oil & Gas Development (China University of Petroleum), Ministry of Education, Qingdao 266580, Shandong, China
^3.School of Chemistry and Chemical Engineering, South China University of Technology, Guangzhou 510641, Guangdong, China
^4.Dongxin Production Plant, Shengli Oil Field, SINOPEC, Dongying 257094, Shandong, China

Received:2024-05-09 Revised:2024-06-17 Online:2024-10-25 Published:2024-11-04
Contact: Mingchuan LI

摘要/Abstract

摘要：

利用质量守恒原理建立了考虑密度跃变的天然气水合物热分解Stefan模型，通过对有弱解的均匀散度方程分部积分详细推导一般形式的Rankine-Hugoniot跳跃关系，获得了关于界面速度项三次方形式的Stefan耦合条件，但在工程需要误差范围内三次方可以忽略。在Stefan模型中引入速度项使得模型的精确解难以得到。Quilghini变换将质量转化为空间变量，Stefan模型的速度项消失，成为经典的类Stefan模型，再经Quilghini逆变换得到了密度跃变Stefan模型的精确解析解。通过Matlab编程，结合实例研究了超越方程解、温度分布和分解界面的规律性，探索注入温度和密度跃变对超越方程解、穿透深度、穿透时间和水合物位移的影响，密度跃变引起的误差在工程允许范围内。

关键词: Quilghini变换, 密度跃变, Stefan模型, Rankine-Hugoniot跳跃关系, 天然气水合物

Abstract:

Stefan's models for thermal decomposition of natural gas hydrate with density jump is established by using the principle of mass conservation. The general form of Rankine-Hugoniot jump relation was derive in detail, and the Stefan coupling condition in cubic form with respect to the velocity of interface was obtained, but the cubic term could be ignored within the errors of engineering requirements. It is the introduction of velocity terms in the Stefan model that makes it difficult to obtain an exact solution to the model. The mass as a space variable was introduced by the Quilghini transformation, the velocity term was vanished, and become the classical Stefan-like's models. Followed by the inverse Quilghini transformation, the accurate analytical solutions of the Stefan's models with density jump could be successfully gained. Through Matlab programming, a concrete example was taken, laws on the solutions of transcendental equation, temperature distribution and dissociation interface were studied, the influences of injection temperature and density jump on solutions of transcendental equation, penetration depth, penetration time and hydrate displacement were probed. The errors caused by density jump were within the permissible ranges of engineering.

Key words: Quilghini transformation, density jump, Stefan's model, Rankine-Hugoniot jump relation, natural gas hydrate

中图分类号:

TQ 511.1

李明川, 樊栓狮, 徐赋海, 卢惠东, 李晓军. Quilghini变换法在求解密度跃变的水合物热分解Stefan模型精确解中的应用[J]. 化工学报, 2024, 75(10): 3763-3774.

Mingchuan LI, Shuanshi FAN, Fuhai XU, Huidong LU, Xiaojun LI. Application of Quilghini transformation method to exact solutions of Stefan's models for thermal dissociation hydrate with density jump[J]. CIESC Journal, 2024, 75(10): 3763-3774.

图/表 12

图1 分解界面S和穿越曲线Σ

Fig.1 Dissociation interface S and an across curve Σ

图2 计算流程

Fig.2 Calculation flowchart

表1 一维平面水合物藏热物性参数

Table 1 Thermo-physical parameter of slab one-dimensional hydrate reservoir

初始压力/MPa

初始温度/℃

密度/(kg·m^-3)

液相水合物分解区

热扩散系数/(m²·s^-1)

固相水合物区

热扩散系数/(m²·s^-1)

液相水合物分解区

热传导率/(W·m^-1·K^-1)

固相水合物区热传导率/(W·m^-1·K^-1)

图3 超越方程f(λ) vsλ

Fig.3 Transcendental equation f(λ) vsλ

图4 考虑Tinj、ρH影响的超越方程f(λ) vsλ

Fig.4 Transcendental equation f(λ) vsλ with effects of Tinj and ρH

图5 温度vs距离，δxs(t)、xd vs时间

Fig.5 Temperature vs distance, δxs(t) and xdvs time

图6 考虑Tinj、ρH影响的温度vs距离

Fig.6 Temperature vs distancewith effects of Tinj and ρH

图7 温度vs时间，δxs (t)、tdvss(t)

Fig.7 Temperature vs time,δxs(t) and tdvss(t)

图8 考虑Tinj、ρH影响的温度vs时间

Fig.8 Temperature vs time with effects of Tinj and ρH

图9 考虑Tinj、ρH影响的s(t) vs时间

Fig.9 s(t) vs time with effects of Tinj and ρH

表A1 受Tinj、ρH影响的变量拟合方程

Table 1 Variable fitting equations with effects of Tinj，ρH

序号	变量	拟合函数(影响因素)		ρ_H引起的误差/%	其他	来源
序号	变量	T_inj	ρ_H	ρ_H引起的误差/%	其他	来源
1	λ	λ=-2.713 $T i n j - 0.1068$ +4.469	λ=0.2857exp(-0.002816ρ_H) + 2.902exp(-4.352×10^-5ρ_H)	0.532999463		图4
2	x_d	x_d=0.4952 $T i n j - 0.1079$ +0.9248	x_d=0.6255exp(0.0007971ρ_H)-39.31exp(-0.009085ρ_H)	1.897220537	x_d=0.3797t^0.4847+0.1226	图5, 图6
3	t_d	t_d=-2.399×10^-7 $T i n j 3$ +0.0001528 $T i n j 2$ -0.03428T_inj+5.262	t_d=0.03634 $ρ H 0.682$	2.977667494	t_d=-0.01667[s(t)]³+0.8089[s(t)]²-9.323s(t)+39.67	图7,图8
4	s(t)	s(t)=-8.573 $T i n j - 0.1068$ +14.12	s(t)= -5.773 $ρ H 0.6051$ +17.61	0.776644191		图9
5	t	t=21.07 $T i n j - 0.2943$ +7.232	t=-8.963×10^-6 $ρ H 2$ +0.003279ρ_H+10.42	0.535378675		图9

表A1 受Tinj、ρH影响的变量拟合方程

Table 1 Variable fitting equations with effects of Tinj，ρH

序号	变量	拟合函数(影响因素)		ρ_H引起的误差/%	其他	来源
序号	变量	T_inj	ρ_H	ρ_H引起的误差/%	其他	来源
1	λ	λ=-2.713 $T i n j - 0.1068$ +4.469	λ=0.2857exp(-0.002816ρ_H) + 2.902exp(-4.352×10^-5ρ_H)	0.532999463		图4
2	x_d	x_d=0.4952 $T i n j - 0.1079$ +0.9248	x_d=0.6255exp(0.0007971ρ_H)-39.31exp(-0.009085ρ_H)	1.897220537	x_d=0.3797t^0.4847+0.1226	图5, 图6
3	t_d	t_d=-2.399×10^-7 $T i n j 3$ +0.0001528 $T i n j 2$ -0.03428T_inj+5.262	t_d=0.03634 $ρ H 0.682$	2.977667494	t_d=-0.01667[s(t)]³+0.8089[s(t)]²-9.323s(t)+39.67	图7,图8
4	s(t)	s(t)=-8.573 $T i n j - 0.1068$ +14.12	s(t)= -5.773 $ρ H 0.6051$ +17.61	0.776644191		图9
5	t	t=21.07 $T i n j - 0.2943$ +7.232	t=-8.963×10^-6 $ρ H 2$ +0.003279ρ_H+10.42	0.535378675		图9

表A2 移动位移δxs（t）拟合方程

Table 2 Fitting equation of δxs（t）

序号	变量	其他拟合函数		来源
1	δx_s₍_t₎	δx_s₍_t₎=-0.07689exp(-0.1695Δt)-0.03523exp(-0.01509Δt)	δx_s₍_t₎=0.0009647[s(t)]³-0.02562[s(t)]²+0.1798s(t)-0.4583	图5,图7

参考文献 33

1	Sloan E D, Koh C A. Clathrate Hydrates of Natural Gases[M]. 3rd ed. Boca Raton, FL: CRC Press, 2008:124-132.
2	Makogon Y F. Hydrates of Natural Gas[M]. Moscow: Nedra, 1974:73-90.
3	李明川, 樊栓狮, 徐赋海, 等. 天然气水合物热分解Stefan相变数学模拟研究[J]. 化工学报, 2021, 72(6): 3252-3260.
	Li M C, Fan S S, Xu F H, et al. Mathematical modeling of Stefan phase change for thermal dissociation of natural gas hydrate[J]. CIESC Journal, 2021, 72(6): 3252-3260.
4	Chelliah H K, Williams F A. Asymptotic analysis of two-reactant flames with variable properties and Stefan-Maxwell transport[J]. Combustion Science and Technology, 1987, 51(4/5/6): 129-144.
5	Mazzeo D, Oliveti G. Parametric study and approximation of the exact analytical solution of the Stefan problem in a finite PCM layer in a steady periodic regime[J]. International Communications in Heat and Mass Transfer, 2017, 84: 49-65.
6	Blanchard D, Porretta A. Stefan problems with nonlinear diffusion and convection[J]. Journal of Differential Equations, 2005, 210(2): 383-428.
7	Briozzo A C, Natale M F, Tarzia D A. Existence of an exact solution for a one-phase Stefan problem with nonlinear thermal coefficients from Tirskii's method[J]. Nonlinear Analysis: Theory, Methods & Applications, 2007, 67(7): 1989-1998.
8	Cherniha R, Kovalenko S. Exact solutions of nonlinear boundary value problems of the Stefan type[J]. Journal of Physics A: Mathematical and Theoretical, 2009, 42(35): 355202.
9	Myers T G, Hennessy M G, Calvo-Schwarzwälder M. The Stefan problem with variable thermophysical properties and phase change temperature[J]. International Journal of Heat and Mass Transfer, 2020, 149: 118975.
10	Li M C, Fan S S, Su Y L, et al. Mathematical models of the heat-water dissociation of natural gas hydrates considering a moving Stefan boundary[J]. Energy, 2015, 90: 202-207.
11	Kanel G I. Shock Waves in Solid State Physics[M]. Boca Raton, FL: CRC Press, 2019: 73-84.
12	Rubinshteĭn L I. The Stefan Problem[M]. Providence, RI: American Mathematical Society, 1971: 24-28.
13	Tritscher P, Broadbridge P. A similarity solution of a multiphase Stefan problem incorporating general non-linear heat conduction[J]. International Journal of Heat and Mass Transfer, 1994, 37(14): 2113-2121.
14	Quilghini D. Su di un nuovo problema del tipo di Stefan[J]. Annali Di Matematica Pura Ed Applicata, 1963, 61(1): 59-97.
15	Quilghini D. Sul problema della conduzione del calore nel muro[J]. Bollettino Della Unione Matematica Italiana, 1965, 20(1): 96-105.
16	Li M C, Fan S S, Xu F H. Interface coupling model of Stefan phase change during thermal dissociation of natural gas hydrate[J]. International Journal of Heat and Mass Transfer, 2021, 175: 121403.
17	Mastroianni G. Sull'approssimazione di funzioni continue mediante operatori lineari[J]. CALCOLO, 1977, 14(4): 343-368.
18	Stefan J. Über die theorie der eisbildung, insbesondere über die eisbildung impolarmeere[J]. Annalen der Physik, 1889, 278(2): 269-286.
19	李明川, 樊栓狮. 天然气水合物注热水分解径向数学模型[J]. 高校化学工程学报, 2013, 27(5): 761-766.
	Li M C, Fan S S. Radial mathematical model for hot water dissociation frontal brim of natural gas hydrates[J]. Journal of Chemical Engineering of Chinese Universities, 2013, 27(5): 761-766.
20	Li M C, Fan S S, Su Y L, et al. The Stefan moving boundary models for the heat-dissociation hydrate with a density difference[J]. Energy, 2018, 160: 1124-1132.
21	Shaikh J, Sharma A, Bhardwaj R. On sharp-interface level-set method for heat and/or mass transfer induced Stefan problem[J]. International Journal of Heat and Mass Transfer, 2016, 96: 458-473.
22	Zhu X J, Dai X R, Liu L, et al. A simplified coupling model of carbonation and chloride ingress based on Stefan-like condition[J]. European Journal of Environmental and Civil Engineering, 2022, 26(7): 2654-2670.
23	Dokoza T, Plümacher D, Smuda M, et al. Solution to the 1D Stefan problem using the unified transform method[J]. Journal of Physics A: Mathematical and Theoretical, 2021, 54(37): 375203.
24	Jeffreys H, Jeffreys B. Methods of Mathematical Physics[M]. Cambridge, UK: Cambridge University Press, 1999:333-374.
25	Needham C E. The Rankine-Hugoniot relations[M]//Blast Waves. Cham: Springer, 2018: 9-17.
26	Chao J K, Zhang X X, Song P. Derivation of temperature anisotropy from shock jump relations: theory and observations[J]. Geophysical Research Letters, 1995, 22(17): 2409-2412.
27	Mitchell S L, Myers T G. Application of standard and refined heat balance integral methods to one-dimensional Stefan problems[J]. SIAM Review, 2010, 52(1): 57-86.
28	Alexiades V, Solomon A D. Mathematical Modeling of Melting and Freezing Processes[M]. London: Routledge, 2018:74-80.
29	Quilghini D. Una analisi Fisico-Matematica del processo del cambiamento di fase[J]. Annali Di Matematica Pura Ed Applicata, 1965, 67(1): 33-74.
30	李明川, 樊栓狮, 徐赋海, 等. 水合物热分解Stefan相变模型解的存在性及Laplace变换求解[J]. 化工学报, 2023, 74(4): 1746-1754.
	Li M C, Fan S S, Xu F H, et al. Existence and Laplace transform of the solution to Stefan phase change model in thermal dissociation hydrate[J]. CIESC Journal, 2023, 74(4): 1746-1754.
31	Kamath V A, Godbole S P. Evaluation of hot-brine stimulation technique for gas production from natural gas hydrates[J]. Journal of Petroleum Technology, 1987, 39(11): 1379-1388.
32	Deaton W M, Frost E M. Gas Hydrates and Their Relation to the Operation of Natural-gas Pipe Lines[M]. Westin Copley: American Gas Association, 1949: 25-28.
33	Reech F. Théorie générale des effets dynamiques de la chaleur[J]. Journal de Mathématiques pures et appliquées, 1853, 18: 357-568.

[1]	李明川, 樊栓狮, 徐赋海, 卢惠东, 李晓军. 水合物热分解Stefan相变模型解的存在性及Laplace变换求解[J]. 化工学报, 2023, 74(4): 1746-1754.
[2]	唐翠萍, 张雅楠, 梁德青, 李祥. 聚乙烯己内酰胺链端改性及其对甲烷水合物形成的抑制作用研究[J]. 化工学报, 2022, 73(5): 2130-2139.
[3]	李明川, 樊栓狮, 徐赋海, 严科, 黄爱先. 天然气水合物热分解Stefan相变数学模拟研究[J]. 化工学报, 2021, 72(6): 3252-3260.
[4]	唐翠萍, 周雪冰, 梁德青. X射线衍射分析聚乙烯吡咯烷酮对水合物分解过程的影响[J]. 化工学报, 2021, 72(2): 1125-1131.
[5]	辛亚男, 张建文, 张淑珍, 姜爱国. 螺旋内槽管内天然气-水-表面活性剂体系的水合物生成动力学计算[J]. 化工学报, 2018, 69(6): 2463-2473.
[6]	阮徐可, 李小森, 徐纯刚, 张郁, 颜克凤. 天然气水合物降压联合井壁加热开采的数值模拟[J]. 化工学报, 2015, 66(4): 1544-1550.
[7]	徐纯刚, 李小森, 蔡晶, 陈朝阳, 陈超. 二氧化碳置换法模拟开采天然气水合物的研究进展[J]. 化工学报, 2013, 64(7): 2309-2315.
[8]	李刚, 李小森. 单井热吞吐开采南海神狐海域天然气水合物数值模拟 [J]. 化工学报, 2011, 62(2): 458-468.
[9]	林曦，刘振宇. 天然气水合物的有效能量及CO2排放强度分析 [J]. CIESC Journal, 2010, 29(12): 2265-.
[10]	杜燕;何世辉;黄冲;冯自平. 多孔介质中水合物生成与分解二维实验研究 [J]. CIESC Journal, 2008, 59(3): 673-680.
[11]	冯自平, 沈志远, 唐良广, 李小森, 樊栓狮, 李清平. 水合物降压分解的实验及数值模拟 [J]. 化工学报, 2007, 58(6): 1548-1553.
[12]	刘有智，邢银全，崔磊军. 超重力旋转填料床中天然气水合物含气量研究 [J]. CIESC Journal, 2007, 26(6): 853-.
[13]	李刚;唐良广;黄冲;冯自平;樊栓狮 . 热盐水开采天然气水合物的热力学评价 [J]. CIESC Journal, 2006, 57(9): 2033-2038.
[14]	周锡堂，樊栓狮，梁德青. CO₂置换开采天然气水合物研究进展 [J]. CIESC Journal, 2006, 25(5): 524-.
[15]	孙始财, 刘玉峰, 吕爱钟, 樊栓狮, 徐海燕. 超声波与表面活性剂协同影响水合物诱导期 [J]. 化工学报, 2006, 57(1): 160-162.

Quilghini变换法在求解密度跃变的水合物热分解Stefan模型精确解中的应用

Application of Quilghini transformation method to exact solutions of Stefan's models for thermal dissociation hydrate with density jump

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 12

参考文献 33

相关文章 15

编辑推荐

Metrics

本文评价