# Hydrate Dissociation Model with Time Fractional Derivative

Authors: Xinyu Fang; Hairong Lian; Wanjing Luo; Mingzhu Liu; Changfu Chen; Qian Wang
Journal: Geofluids
 (2022)

---

## Abstract

In this paper, we shall investigate fractional partial differential equations with fractional moving boundary condition to study the dissociation of natural gas hydrate under heat injection. The moving boundary separates the hydrate reservoir into the dissociated zone and the hydrate one. By using the self-similar transformation and Wright function, we obtain the explicit solutions for two zones. We present simulations with steam and hot water injection and show the dissociation temperature in graphical mode from injection temperature to reservoir temperature with respect to the time, distance, and fractional order. Our analysis of fractional model turns out to be a successful generalization of the classical one; i.e., it can well describe the dissociation of natural gas hydrate and is theoretically consistent with the existing integer hydrate dissociation model. When the factional order tends to 1, the “limit solution” becomes the classical one.


---

## Body

## 1. Introduction

Amounts of natural gas hydrates have been found in permafrost and ocean bottom sediments. Natural gas hydrates are an unconventional source of energy and belong to the clean energy. Because of the characteristics of large storage scale and high energy density, natural gas hydrate become one of the potential new energy sources. There are many researchers focusing their study on the exploitation and mathematical model of natural gas hydrate.Thermal models are investigated for the dissociation of natural gas hydrate. Holder et al. [1] established a zero-dimensional hydrate model. Ignoring energy loss, heat, and mass transfer, the authors revealed the feasibility of hydrate dissociation by thermal method. Selim and Sloan [2] viewed the hydrate decomposition process as a moving boundary ablation process and established one-dimensional hydrate decomposition models for temperature change under heat injection method. The classical parabolic differential equations based on the mass, momentum, and energy balances were used to describe the dissociation process. More precisely, let T be temperature and x=χt be a moving boundary, and the moving boundary condition was written as
(1)khII∂TII∂x−khI∂TI∂x=ϕρHQHddχdt,x=χt,t>0,where some physical quantitieskhI,khII,ϕ,ρH,andQHd are of given values. An analytical solution was derived, and results were used to evaluate gas production. Later, the authors in [3–8] generalized the analytical model to radial symmetric decomposition model, well structure, etc. For hydrate decomposition model solved numerically, we refer to references [9, 10] and their related references.Fractional calculus has been proposed to scientists and engineers for decades. And very soon it turned out to be a very powerful tool both in theoretical analysis and wide applications, as one can see its influence in the study of solute transport, proteins diffusion, viscoelasticity, and other problems [11–13]. For example, in [11], fractional differential equations were established for fractal mobile/immobile solute transport. The solutions were gained by performing an integral transform. Results captured the anomalous behavior of tracer plumes in heterogeneous aquifers, which showed that a fractional model can describe the power law breakthrough curve more accurately than the classical model. Roscani et al. [14] investigated two kinds of fractional differential equations with Stefan-like boundary condition. One model was the Caputo one, written as
(2)0cDtαu1x,t=λ12∂2u1x,t∂x2,0<x<χt,0<t<T,0cDtαu2x,t=λ22∂2u2x,t∂x2,x>χt,0<t<T,whereα∈0,1, D0ctαux,t is the α‐order Caputo fractional derivative of ux,t with respect to t. The author presented explicit solutions of fractional two-phase Stefan-like problems in terms of Wright functions.Heat equations have many similarities to diffusion ones; therefore, fractional models are naturally introduced to study heat transfer problem [15]. Apart from the aforementioned fields, we found that very few study of fractional order models have been carried out on hydrate dissociation.Motivated by the literatures listed above, we will discuss a time fractional heat conduction model in this paper to describe the hydrate decomposition process for heat injection case. Because the moving boundary is usually linked to melting (change from solid to liquid) or freezing (change from liquid to solid). We also consider the fractional moving boundary. Results from our model are used to calculate temperature distributions from hydrate reservoirs by steam or hot water thermal stimulation.There are also many literatures discussing the depressurization model of hydrate decomposition [16–20]. As an attempt to focus on the discussion of fractional methods, in this paper, we shall not divert our attention away from the methods to cumbersome calculation. That is, we shall not consider the depressurization and the mass transfer model.


## 2. Hydrate Decomposition Model

Consider the hydrates that are uniformly distributed in a semi-infinite porous medium (0<x<∞). It has an initial reservoir temperature Ti. At time t=0, hot water is injected into the wellbore and reservoir temperature near the wellbore is immediately raised to injection temperature T0T0>Ti. Then, the balance is disrupted and the hydrates begin to dissociate. The hydrate decomposition can be considered a moving boundary value problem. It has a decomposition front, which divides the reservoir into two parts: the decomposition zone (zone I, 0<x<χt) and the hydrate zone (zone II, χt<x<∞), which are shown in Figure 1.Figure 1
Hydrate decomposition diagram.Before establishing the mathematical model, we give some assumptions here.(1)
The hydrate reservoir is a one-dimensional infinite region. It is saturated(2)
Only water and ideal state natural gas are dissociated from the hydrate, without carrying ice(3)
The thermal conductivity and diffusivity are constant. The thermal convection is ignored(4)
There is no energy loss for water dissociated diffusion, heat transfer from the down-hole to the dissociation and from hydrate zones into the top or bottom caps(5)
The gravity effect, the viscous dissipation, and reversible work of compression are all neglectedWith the restrictions listed above, the energy conservation can be written as(3)ρH∂u∂t+divq=0.where u is the total internal energy and q is heat conduction flux per unit volume. The relationship between total internal energy and temperature is du=CdT, where C is the specific heat capacity. Noticing q=−k∇T and we have the one-dimensional heat conduction equation as follows:
(4)∂Tx,t∂t=α∂2Tx,t∂x2,whereα=k/ρC. Furthermore, we consider a memory effect of power law decay on temperature variation with respect to time. Let ft=tα/Γ1−α,α∈0,1 be the ‘memory function.’ Then, the convolution of temperature derivative with respect to time and power function turns out to the Caputo fractional derivative, that is,
(5)∂Tx,t∂t∗ft=0cDtαTx,t,where ∗ denotes the convolution. Here, we claim that the α-order Caputo fractional derivative of Tx,t on t for α∈0,1 is defined by
(6)D0ctαTx,t=1Γ1−α∫0tt−τ∂Tx,τ∂τdτ.For more details, we refer the author to see the reference [21].So the Caputo fractional temperature distribution equations are established for the dissociation zone and hydrate zone, respectively.(7)0cDtαTIx,t=αI∂2TIx,t∂x2,0<x<χt,0<t<T,(8)0cDtαTIIx,t=αII∂2TIIx,t∂x2,x>χt,0<t<T.The boundary and initial conditions associated with the above equations are(9)TIIx,0=Ti,(10)TI0,t=T0,TII+∞,t=Ti,(11)TIχt,t=TIIχt,t=Td,as well as the fractional-order moving boundary condition:
(12)ϕρHQHd0cDtαχt=kII∂∂xTIIx,tχt+,t−kI∂∂xTIx,tχt−,t.The unknown functions of our model are the temperaturesTIx,t,TIIx,t and moving boundary χt.


## 3. A Similarity Solution

For everyx∈ℝ, the Wright function is defined as
(13)Wx,ρ,β=∑k=0∞xkk!Γρk+β,where the parameters satisfyρ>−1 andβ∈ℝ. Mainardi function is a special Wright function, written as
(14)Mρx=W−x,−ρ,1−ρ=∑k=0∞−xkk!Γ−ρk+1−ρ,0<ρ<1.By direct calculations, we can obtain the following results. The proof is omitted here.Proposition 1.
Letβ≥0,0<ρ<1. The Wright function has the following properties:
(1)
W−x,−ρ,β is positive and strictly decreasing on x∈0,+∞:(15)ρxW−x,−ρ,β−ρ=W−x,−ρ,β−1+1−βW−x,−ρ,β,x≥0,(16)limx⟶+∞W−x,−ρ,β=0(2)
limα⟶11−W−2x,−α/2,1=erfx, where erfx=2/π∫0xe−u2du(17)∂∂xWx,ρ,β=Wx,ρ,ρ+β,x∈RNext, we will solve our fractional model (7)–(12) and present explicit solutions in terms of the Wright functions. Adopting the self-similar transformation:
(18)x=2ηαIItα/2,and we assume that form of the moving boundary is χt=2ξαIItα/2, where ξ>0 is undetermined. In fact, we can solve equation (12) to get ξ later. The dissociation front of hydrate can be determined by ξ, which is related to the physical properties of hydrate reservoir.Suppose that the self-similar solution to problem (7)–(12) can be written by
(19)TIη=A+B1−W−2λη,−α2,1,TIIη=C+D1−W−2η,−α2,1,whereA,B,C,andD are some unknown constants, which just makes the solution a little bit more elegant.Letλ=αII/αI and it is related to the diffusivity of decomposition zone and hydrate zone. From equations (4) and (5), noticing that when x=0 or x=χt, η=0 or η=ξ, respectively, it has
(20)A+B1−W0,−α2,1=T0,A+B1−W−2λξ,−α2,1=Td,which arrives atA=T0,B=−T0−Td/1−W−2λξ,−α/2,1. Similar, from
(21)TIIx,0=C+D1−W−∞,−α2,1=C+D=Ti,(22)TIIχt,t=C+D1−W−2ξ,−α2,1=Td.We can obtain thatC=Ti+Td−Ti/W−2ξ,−α/2,1,D=−Td−Ti/W−2ξ,−α/2,1.So the solutions of problems (7)–(12) are given by
(23)TIη=T0−T0−Td1−W−2λξ,−α/2,11−W−2λη,−α/2,1,TIIη=Ti+Td−TiW−2ξ,−α/2,1W−2η,−α2,1.The constantξ>0 can be calculated by the transcendental equation, which can be obtained by equation (23). The transcendental equation is expressed by
(24)khIT0−TdΓ1−α/2αIαIIF12λξ−khIITd−TiΓ1−α/2αIIF22ξ=2Γ1+α2ϕρHQHdξ,where(25)F1ξ=Mα/2ξ1−W−ξ,−α/2,1,F2ξ=Mα/2ξW−ξ,−α/2,1.Whenα=1 in equations (7)–(12), the Caputo fractional temperature distribution model becomes the classical one [2].


## 4. Result and Discussion

Two cases are simulated in this section to show the temperature distribution when hot steam and hot water are injected, respectively. The relationship of time, decomposition distance, and temperature are investigated for both cases under given fractional order. Initial heat injection temperature and the decomposition front are also discussed. The influence of fractional order is also shown in the last for the hot-water-injecting case.ChooseαI=2.89×106μm2/s, αII=0.697×106μm2/s, ρH=913kg/m3, ϕ=0.3,.kI=5.57W/m⋅K, kII=2.73W/m⋅K, and QHd=446120+132.638×0.5T0+Ti. The injection temperature, boundary temperature, and reservoir temperature are 563.5 K, 282.5 K, and 275 K for steam case, respectively. And for hot water case, they are 450 K, 285 K, and 280 K, respectively.Letα=0.9. When hot steam and hot water are injected, by Newton iterative computation, the values of ξ are 1.97 and 1.675, respectively. Furthermore, when decomposition distance and time are given, the temperature can be calculated from (23) and are shown in Figure 2. Figure 2(a) shows the temperature with distance for steam case and the right for hot water case.Figure 2
The temperature of steam and hot water varies with time.Easily, the dotted line, in both figures of Figure2, divides the temperature distribution into the above and the below area. The above one corresponds to the temperature of the hydrate decomposition area, and the below is the hydrate one. The temperature decreases rapidly in the decomposition zone and very slowly in the hydrate zone when the axial distance is increasing. It does not change and keeps the initial temperature all the time in most hydrate area. This is because of the thermal diffusion coefficient is very small in the hydrate zone and the heat here cannot be conducted out. Comparing the left and right Figure 1, it can be obtained that the higher the initial heat injection temperature is, the slower the decomposition zone temperature decreases, and the longer the axial distance of the decomposition leading edge reaches, the larger the decomposition zone is. Obviously, increasing the initial heat injection temperature can make hydrates decompose more quickly.Figure3 shows that the temperature distribution at x=0.2,x=0.5,andx=0.8 meter distances. The small figure in the box is the local enlarged one at the beginning.Figure 3
The temperature of steam and hot water varies with time.With the increase of time, the reservoir temperature in some fixed location slowly rises to the hydrate decomposition temperature and then rises rapidly to tend to the initial heat injection temperature. Meanwhile, the closer the location is to the wellbore, the faster the temperature rises and the shorter the time needed to approach toT0. On the contrary, the further away from the wellbore, the longer the time required for the temperature rise to approach the initial heat injection temperature.When the fractional order is increasing, the value ofξ is increasing too. In fact, for the hot-water-injecting case, when α=0.8,0.9,and1, the values of ξ are 1.905, 1.97, and 2.01, respectively. The temperature distribution for dissociation, calculated by equations (23) and (24), are presented in Figures 4 and 5.Figure 4
The temperature distribution at different distances.Figure 5
The temperature distribution at different times.Figure4 shows the change of temperature when t=105s and α=0.8,0.9,and1, respectively. When the fractional order α increases, the temperature decreases more slowly, and the dissociation interface χt is farther from the wellbore. Figure 5 shows the change of temperature when x=0.8m and α=0.8,0.9,and1, respectively. And we also know from Figure 5 that when α increases, the temperature will reach the initial heat injection temperature in a shorter time.Figure5 shows the change of temperature of steam case when x=0.8 and α=0.8,0.9,and1, respectively. Obviously, with the increase of time, the temperature of hydrate decomposition increases slowly and then increases rapidly close to the initial heat injection temperature. The higher the time fractional order, the faster the temperature increases.The constantξ affects the move speed of the dissociation front. Figure 6 shows that, with the increase of initial heat injection temperature, ξ will increase and the decomposition rate will also accelerate. Meanwhile, the value of ξ is also increasing when the fractional order α is increasing. When α tends to 1, the value of ξ is the same as that of the classical hydrate decomposition.Figure 6
The variation of the decomposition front with the initial temperature and fractional order.Whenα=0.9 and 1, the ξ value of the corresponding hot steam injection is 1.97 and 2.01, respectively. Hydrate saturation can be calculated from moving boundary, and both are shown in Figure 7 for both fractional and integral hydrate models. They decrease gradually with the increase of time, while the fractional one has a tailing phenomenon.Figure 7
The variation of hydrate saturation with time.


## 5. Conclusions

By using time fractional derivative, a new mathematical model is investigated for the temperature variation in decomposed and undecomposed (hydrate) zones of the semi-infinite hydrate reservoir. The study shows that the fractional-order model can well describe the change rule of temperature on hydrate decomposition. When the fractional order and the time are given, the temperature drops sharply to the dissociation temperature and then gradually to the reservoir temperature as the distance from the decomposition front increases. And when the decomposition distance is fixed, the temperature slowly rises to the decomposition temperature, and then rapidly to the initial heat injection temperature with the passage of time. With the increase of the fractional order at different times, the temperature gradient becomes smaller, the temperature drops more slowly, and the decomposition distance can be reached further. With the increase of the fractional order at different distances, the temperature will reach the initial heat injection temperature in a shorter time.The variation of temperature with time during the classical hydrate dissociation model was obtained in [2, 4]. The fractional-order model is consistent with the classical dissociation model in terms of temperature variation and is an extension of the classical hydrate model. Furthermore, there is a need to extend the present analysis to both heat and mass transfer, which will be directly applicable to evaluate gas production by thermal stimulation from hydrate reservoirs.



---

*Source: 5598287-2022-04-12.xml*