Next Article in Journal
A Novel Acoustic Method for Cavitation Identification of Propeller
Next Article in Special Issue
High-Frequency Dependence of Acoustic Properties of Three Typical Sediments in the South China Sea
Previous Article in Journal
Effects of the Ionic Liquid [BMIM]Cl on the Baltic Microphytobenthic Communities
Previous Article in Special Issue
Normalized Stress–Strain Behavior of Deep-Sea Soft Soils in the Northern South China Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Rheology of CO2 Hydrate Slurry by Using the Capillary Method

1
State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, Xuzhou 221116, China
2
CNPC Research Institute of Exploration and Development (Langfang), Langfang 065000, China
3
School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(9), 1224; https://doi.org/10.3390/jmse10091224
Submission received: 7 July 2022 / Revised: 29 August 2022 / Accepted: 30 August 2022 / Published: 1 September 2022
(This article belongs to the Special Issue Advances in Marine Engineering: Geological Environment and Hazards)

Abstract

:
The rheological properties of the CO2 hydrate slurry are experimentally investigated by using a new flow loop with an inner diameter of 50 mm. The pressure drops of the CO2 hydrate slurry are measured experimentally under the volumetric hydrate fraction, ranging from 1.4 to 17.2 vol%, and the shear rate ranging from 40 to 590 s−1. Using the capillary method, the rheology of the CO2 hydrate slurry is analyzed based on pressure drop. The CO2 hydrate slurry is identified as a power-law fluid and exhibits strong shear-thickening behavior. The increase in the volumetric hydrate concentration not only ascends the apparent viscosity, but also makes the non-Newtonian behavior of the hydrate slurry become more obvious. The non-Newtonian index and the consistence factor of fluid increase exponentially with the volumetric hydrate concentration increasing. A Herschel–Bulkley-type rheological model of the CO2 hydrate slurry is correlated from experimental data. The developed model performs the average discrepancy of less than 16.3% within the range of the experiment.

1. Introduction

Initially, the CO2 hydrate is a secondary refrigerant used to reduce the usage of the primary refrigerants and to improve the efficiency of the secondary system [1]. Understanding the rheology of the CO2 hydrate slurry can predict the pressure drop of refrigerating pipes accurately and help to optimize the efficiency of the refrigerating system. In recent years, as the development of oil and gas in deep-water regions progresses [1], the reservoir gases become compounds, including methane, ethane, propane, butane, CO2 and H2S [2,3]. The acid gas (CO2 and H2S) can form gas hydrates in the wellbore much more easily than methane. Thus, the hydrate flow assurance problem in the wellbore and pipe becomes more complicated when natural gas contains CO2.
Currently, the rheology of the hydrate slurry is a fundamental issue in calculating the pressure drop along the pipe length and time, and plays an important role in predicting hydrate plugging in the wellbore. However, CO2 hydrate accelerates the process of gas hydrate formation and amplifies the pressure drop induced by hydrate slurry flowing. Thus, the current prevention method for gas hydrate underestimates the risk of gas hydrate containing acid gas. The current hydrate management plan becomes improper. Studies focusing on the rheology of the rCO2 hydrate slurry contributes to the development of a high-efficiency hydrate management plan. Moreover, the rheology of the CO2 hydrate slurry participates in a number of industrial applications as an important fundamental concept, such as CO2 secondary refrigeration and CO2 storage in deep-water stratum.
Generally, the rheology of the hydrate slurry is influenced by multiple factors, including mesoscopic structure between hydrate particles, hydrate concentration, rheology of the carry fluid, additive concentrations, etc. [4,5,6]. Since the hydrate slurry only keeps stable under the high-pressure and low-temperature environment, only two kinds of viscosimeters can satisfy the experimental requirement after a customizing upgrade, which are the Ostwald viscosimeter and rotatory viscosimeter [4,5,6,7,8,9,10,11]. Based on the kinds of carry fluid, the studies of hydrate slurry rheology can be divided into two parts [11]. The first part is the hydrate slurry formed from the water-in-oil emulsion system [5,7,8,9,10]. The other is the hydrate slurry formed from the aqueous solution [12,13,14,15,16,17,18,19,20,21]. In the water-in-oil emulsion system, Sinquin et al. (2004) pointed out that the increase in pressure drop by the hydrate slurry is dominated by the friction factor under turbulent flow and the apparent viscosity under laminar flow [7]. In the most of their experiments, the hydrate slurry exhibits Newtonian behavior, and the relative viscosity of the hydrate slurry increases with the volume fraction of hydrate [7]. In some specific cases, the hydrate slurry exhibits shear-thinning behavior due to high enough interparticle forces between hydrate particles [7]. Peng et al. (2012) and Yan et al. (2014) conducted a series of experiments to investigate rheological properties of natural gas hydrate slurries with varied anti-agglomerants [5,8]. They also used the flow loop to test flow performance of the hydrate slurry under the shut-down and the restarting condition [5,8]. The hydrate slurry in a water-in-diesel oil emulsion exhibits obvious shear-thinning behavior and the increase in hydrate volume fraction makes the hydrate slurry become more shear-thinning. An empirical Herschel–Bulkley-type equation is developed, considering the hydrate volume fraction to describe the rheological property of the hydrate slurry [5,8]. Webb et al. (2012, 2014) experimentally investigated the rheology of the hydrate slurry in the water-in-crude oil emulsion, which varied in time, shear rate (1~500 s−1), water cut (0~50%), temperature (0~6 °C) and initial pressure of methane (750~1500 psig) [22,23]. The yield stress of the hydrate slurry was found to increase with an annealing time of up to 8 h, and then remain constant. The hydrate slurry performed shear-thinning behavior and the increase in water cut resulted in an increase in the apparent viscosity and the yield stress [22,23]. Qin et al. (2017) used a control-stress rheometer to study the rheology of hydrate slurries in water-in-oil emulsion, which found that the hydrate formation induced the significant increase in apparent viscosity [9]. The yield stress of the hydrate slurry in oil varied between 3 and 25 Pa at the water cut, ranging from 0.05 to 0.3 [9]. Liu et al. (2020) experimentally analyzed the influence of the surfactant’s concentration on the rheology of the hydrate slurry. The surfactant had a positive contribution to the stability of the hydrate slurry [24]. When the surfactant concentration was beyond 5 wt%, the viscosity of the slurry decreased significantly [24]. The yield stress of the hydrate slurry was related to the shut-in time. In the first 10 min of shut-in, the yield stress was constant. When the shut-in time is beyond 20 min, the yield stress increased greatly. Liu et al. (2021) studied the rheology of the hydrate slurry containing wax and defined the three stages of apparent profile, which are the initial growth and aggregation, dynamic recovery and the balanced stage [10]. Wax can increase the viscosity of the hydrate slurry but weaken the thixotropy of the hydrate slurry [10]. The most proper reason that wax can induce an extra increment of viscosity for the hydrate slurry is due to the emergence of wax-hydrate coupling aggregates [10].
In the aqueous solution, Andersson and Gudmundsson (2000) firstly used the flow loop to text the rheology of hydrate-in-water and pointed out that the hydrate-in-water exhibits Bingham behavior [12]. Darbouret et al. (2005) experimentally investigated the rheology of TBAB hydrate slurries as the requirement from the secondary refrigeration field [25]. The Bingham behavior of TBAB hydrate slurries was found in experiments of the ranges of particle contents and shear rates, and a corresponding method was proposed to decide its apparent viscosity and yield stress [25]. Delahaye et al. (2008, 2011) conducted research on the rheology of the CO2 hydrate slurry with additives [1,15]. The rheology of the CO2 hydrate slurry is highly dependent on the volumetric hydrate fraction [1,15]. Without the effect of additives, the CO2 hydrate slurry exhibits shear-thickening behavior with a volumetric hydrate fraction ranging from 4 to 10 vol%, and shear-thinning behavior with volumetric hydrate fraction ranging from 11 to 20 vol% [1,15]. Hashimoto et al. (2011) used a stress-control-type rheometer to study the rheology of TBAB and TBAF hydrate slurries, and pointed out their rheology exhibits pseudoplastic [16]. Clain et al. (2012, 2015) experimentally verified the shear-thinning behavior of the TBPB hydrate slurry with a volumetric hydrate concentration ranging from 0 to 28.2 vol% [4,25]. Meanwhile, the particle size of the TBPB hydrate slurry was confirmed at range between 75 to 160 um by focused-beam reflectance measurement (FBRM) [25]. The hydrate particle growth and agglomeration were indicated to have a potential relation with the rheological properties of the hydrate slurry. Jerbi et al. (2013) used the capillary viscometer method to reveal the rheological properties of the CO2 hydrate slurry by using a flow loop with an inner diameter of 7.74 mm [17]. The CO2 hydrate slurry is a power fluid and the non-Newtonian index decreased from 1 to 0.6 as the volumetric hydrate concentration increased from 0 to 22 vol% [17]. Sun et al. (2020) experimentally revealed the rheology of the methane hydrate slurry in the bubbly flow and pointed out that the methane hydrate slurry with bubbly exhibited shear-thinning under low-shear-rate conditions and shear-thickening under high shear rate condition [18]. Fu et al. (2020a,b, 2022) used the same method as Sun et al., to investigate the influence of hydrate formation on water-based drilling fluid [6,19,21]. The hydrate has an obvious effect on the rheology of the drilling fluid. The increase in volumetric hydrate concentration enhances the non-Newtonian behavior of the drilling fluid. Under the high shear rate condition, the water-based drilling fluid with hydrate also tends to the shear-thickening behavior [6,19,21]. Muhlstedt et al. (2021) used a stress-control rotational rheometer to study the rheology of water-based drilling fluid with THF hydrate. The influence of hydrate formation on the shear-stress and shear rate were analyzed [26]. Based on Fu et al., and Muhlstedt et al.’s work, although the effect of the hydrate on the rheology of the drilling fluid is negligible, the rheology of water-based drilling fluid is highly dependent on the additives of drilling fluid. Sakurai et al. (2021) revealed the hydrate formation process in the water-dominated bubbly flow and verified Sun et al.’s result in which the hydrate formation in the bubbly flow can increase the shear stress of the bubbly flow [27].
Although previous researchers have conducted experiments to reveal the rheology of the CO2 hydrate slurry, a consistent conclusion has not been reached. In Delahaye et al.’s work, the rheology of the CO2 hydrate slurry could be shear-thinning and shear-thickening depending on volumetric hydrate concentration [1,15]. In Jerbi et al., and Lv et al.’s work, the CO2 hydrate slurry was found to be shear-thinning fluid [27,28]. In this work, a new flow loop with an inner diameter of 50 mm is built up to study the rheology of the CO2 hydrate slurry by using the capillary method, where the 50 mm inner diameter of the flow loop is quietly closed to the industrial conditions and the experiments are designed based on the Reynolds number similar criterion. Different from previous works, the CO2 hydrate slurry exhibits an apparent shear-thickening behavior in the whole experiment. An empirical Herschel–Bulkley-type equation is developed to predict the shear stress and the apparent viscosity of the CO2 hydrate slurry.

2. Experiment

In this work, a new flow loop was built up to study the rheology of the CO2 hydrate slurry. Compared to our previous flow loop, the new loop had a larger inner diameter and longer test section, which were 50 mm and 15 m, respectively. The large diameter of test section pipe provided a flow condition closed to the industrial condition, where the drainage line in deep-water well had a diameter of 80 to 100 mm. The centrifugal pump equipped in the flow loop could tolerate 10 vol% volumetric gas fraction which enabled the void fraction and hydrate fraction inside the flow loop to be up to 10 vol% and beyond 12 vol%, correspondingly. The accuracies of temperature, pressure, pressure drop and flow rate were 0.01 °C, 0.001 MPa, 0.01 Kpa and 0.01 m3/s. The tolerances of pressure and temperature were 5 Mpa and 0 °C. The interval of data collection was 1 s. The transparent PVC pipe enabled the observation of flow pattern and morphology of the hydrate slurry but limited the maximum pressure of the flow loop. The total volume of the flow loop was 43,000 mL. The image and schematic of the new flow loop are given in Figure 1a,b.
The experiments were designed based on the Reynolds number similar criterion, where the CO2 hydrate slurry had the same Reynolds number with fluid in the drainage line of a deep-water natural-gas hydrate well. The rheology of the CO2 hydrate slurry was tested by the capillary method. The capillary method uses the pressure drop of the hydrate slurry in the horizontal pipe to evaluate its rheology. Since the rheology is measured in the pipe, the flow condition in the experiment was the same as the deep-water well. Thus, the rheological model developed based on the capillary method is more reliable for the petroleum engineers. The experiments had void fractions ranging from 0 to 7.09 vol%, flow velocity ranging from 0 to 4 m/s, volumetric hydrate concentration ranging from 1.4 to 17.2 vol% and pressure ranging from 0 to 4 Mpa. The experimental procedure is illustrated as follows:
(1)
In the preparation stage, the flow loop is cleaned up by the deionized water and vacuumed by the vacuum pump. When the pressure of the flow loop drops to 0.1 MPa, the flow loop is vacuumed. Then, the plunger pump is used to pump water inside the flow loop until the expected volume of water and an electric weight can weigh the pumped water in real time. At this stage, the data collection system starts to record experimental data, such temperature, pressure, pressure drop, flow rate, etc.
(2)
Then, the valve of the gas cylinder opens and CO2 is injected inside the flow loop. At the same time, the cooling system starts to cool down the test section and the pressure of the flow loop is kept under 1.2 MPa. In this stage, the temperature of the flow loop decreases and CO2 dissolves into water. When the temperature drops to the target and the pressure is kept the constant, the water is saturated with CO2.
(3)
CO2 is injected into the flow loop until the pressure reaches the target and the valve of the gas cylinder closes. CO2 hydrate starts to crystalline in the flow loop, and the pressure of the flow loop decreases correspondingly.
(4)
When the pressure keeps the constant, the hydrate formation stops in the flow loop. Thus, the flow loop varies from 37 Hz to 7 Hz and records the flow velocity and pressure drop of the hydrate slurry at the same time.
(5)
CO2 is injected into the flow loop again and the second hydrate formation experiment is conducted to increase the volumetric hydrate concentration in the flow loop. When the void fraction decreases to be negligible, and the volumetric hydrate concentration becomes hard to increase, the experiment stops.
The rheological experiments of the CO2 hydrate slurry consist of a CO2 hydrate formation experiment and rheological examination experiment. The CO2 hydrate formation experiment increases the volumetric concentration in the flow loop. The rheological examination experiment tests the pressure drop of the hydrate slurry under a different flow velocity, where such pressure drops can evaluate the rheology of the hydrate slurry by the capillary method. The whole procedure of the rheological experiment is given in Figure 2a,b.
In Figure 2a, the experiment consists of 20 times of CO2 injections, which means 20 times of CO2 hydrate formation experiments. The volumetric hydrate concentration is increased after each hydrate formation experiment and gradually ascends from 1.4 to 17.2 vol%. After hydrate formation, the pump frequency (fpump) decreases from 37 Hz to 7 Hz and the flow velocity of the hydrate slurry also decreases correspondingly. In Figure 2a, as the volumetric hydrate fraction increases, the range of flow velocity shrinks at the same time. This is because the increased viscosity of the hydrate slurry induces higher friction force with the pipe wall, and the centrifugal pump needs more frequency to maintain the velocity of the hydrate slurry. At 37 Hz, the velocity of the hydrate slurry decreases from 4 to 2.16 m/s, which is about a half, as the volumetric hydrate concentration increases from 1.4 to 17.2 vol%. Figure 2b exhibits the effect of volumetric hydrate concentration on the pressure drop of the hydrate slurry, where first and third hydrate formation experiments are marked. The pressure drop of the hydrate slurry increases greatly as the volumetric hydrate concentration increases. For instance, under the hydrate slurry velocity of 15.5 m/s, the pressure drops are 6 and 65.8 kPa at the volumetric hydrate concentration of 1.4 and 17.2 vol%. The pressure drops increase by about 11 times.

3. Rheological Model

Previous researchers have performed numbers of experiments to study the rheology of hydrate slurries in aqueous solutions, as shown in Table 1 [1,4,6,12,13,14,15,18,19,21]. All those works use the capillary method (Ostwald viscosimeter) to investigate the rheology of the hydrate slurry. However, their rheology does not meet a unique conclusion which exhibits Bingham, shear-thinning and shear-thickening in different aqueous solutions. If the aqueous solution contained additives, the hydrate slurry would be highly related to the physical and chemical properties of additives. Compared with Delahaye et al.’s work, we used the flow loop with a larger inner diameter to research the rheology of the CO2 hydrate slurry and a wider range of apparent viscosities. The flow loop with a larger dinner diameter enables the flow condition of the hydrate slurry in experiments to be closed to the industrial condition and the experimental results become reliable.
Although no clear agreement for the CO2 hydrate slurry rheology is reached, the hydrate slurry is recognized as a power-law fluid which is a non-Newtonian fluid. In this model, the rheology of the CO2 hydrate slurry is assumed to be a laminar flow regime in the cylindrical pipe under most of the experimental conditions. Since the wall slip is widely considered in the multiphase flow in porous media or in the micro-tube, the flow behavior of the hydrate slurry in the pipeline is a macroscopic problem, and the other influencing factors are more important than the wall slip, such as hydrate concentrations and flow velocity [18]. Therefore, we ignore the wall slip in this work. The Metzner and Reed model (1955) is a theoretical and fundamental base to describe the rheology of the hydrate slurry which revealed the relationship between flow rate, shear rate and shear stress, as shown in Equation (1) [29]:
Q π R 3 = 1 π 3 0 τ w τ 2 γ · d τ
where Q is the flow rate of the hydrate slurry, m3/s, R is the inner diameter of the pipe, m, τ is the shear stress, Pa, γ is the shear rate, s−1.
The shear stress of the hydrate slurry is formulated by the pressure drop of the hydrate slurry, test section length and the diameter of the pipe [30], as shown in Equation (2):
τ H S = D 4 Δ P L
where τHS is the shear stress of the hydrate slurry, Pa, ΔP is the pressure drop of the hydrate slurry, Pa, D is the diameter of the pipe, m, L is the length of the test section pipe, m.
Moreover, the Herschel–Bulkley equation is a classical expression for the non-Newtonian fluid to illustrate the relationship between shear rate and shea stress [30], as shown in Equation (3):
τ H S = τ 0 + K γ H S n
where τ0 is the yield stress of the hydrate slurry, Pa, γHS is the shear rate of the hydrate slurry, s−1, n is the non-Newtonian fluid index, K is the consistency factor. Equation (4) defines the n index [30]:
n = n = d l n D 4 Δ P L d l n 8 v H S D
In our previous works, the apparent viscosity of the hydrate slurry was defined as the ratio of shear stress and shear rate of the hydrate slurry [6], as shown in Equation (5):
μ a p p = τ H S γ H S = τ 0 γ H S + K γ H S n 1
where μapp is the apparent viscosity of the hydrate slurry, mPa s.

4. Experimental Results and Discussion

The rheology of the CO2 hydrate slurry in current experiments are studies within volumetric hydrate concentrations from 1.4 to 17.2 vol%, shear rates from 40 to 590 s−1 and void fractions from 7.09 to 2.32 vol%. The rheology model of the CO2 hydrate slurry is developed based on the Herschel–Bulkley equation (power-law model). The volumetric hydrate concentration is taken into account when describing the non-Newtonian index (n factor) and the consistency factor (K factor).
The n index is obtained from the differential ratio of Ln(DΔP/(4L)) and Ln(8vHS/D) which is the slope of the curve in Ln(DΔP/(4L)) verse Ln(8vHS/D) as an expression of Equation (4). The pressure drop is directly measured in experiments and the flow velocity of the hydrate slurry is calculated based on the flow rate measured in experiments. As shown in Figure 3, under each volumetric hydrate concentration condition, Ln(DΔP/(4L)) and Ln(8vHS/D) have a linear relationship between them, indicating that they have an unique n index under a certain hydrate concentration condition. When CO2 hydrate starts to form in the flow and the volumetric hydrate concentration increases from 1.4 to 17.2 vol%, the n index varies from 1.53 to 1.79. The n index is always above 1, which indicates that the CO2 hydrate slurry performs shear-thickening behavior.
In Figure 4, the relationship between shear stress and shear rate is exhibited. The shear stress increases exponentially while the shear rate increases, which identifies that the CO2 hydrate slurry is a typical power-law model. Compared with Equation (3), the coefficient of the regression equation for shear stress and shear rate is the consistency factor of the CO2 hydrate slurry. As the volumetric hydrate concentration ranges from 1.4 to 17.2 vol%, the K factor varies from 0.0016 to 0.0095. Moreover, the regression equation has no intercept with shear stress. The yield stress of the CO2 hydrate slurry is zero, and Equation (5) is simplified as μapp = KγHSn−1.
Based on the simplified Equation (5), Figure 5 is plotted to present the relationship between shear rate and the apparent viscosity of the CO2 hydrate slurry. The apparent viscosity increases exponentially with the shear rate increasing, where the CO2 hydrate slurry is demonstrated by shear-thickening behavior. As the volumetric hydrate concentration increases from 1.4 to 17.2 vol%, the ascending trend of the apparent viscosities along the shear rate become more obvious. Generally, the increase in the volumetric hydrate concentration enhances the non-Newtonian behavior of the CO2 hydrate slurry. This result reaches an agreement with our previous works.
In this work, the formation of CO2 hydrate in the flow loop induces a CO2–water two-phase flow transits to become a CO2–water–CO2 hydrate three-phase flow. CO2 hydrates result in the non-Newtonian behavior of the hydrate slurry and the increase in the apparent viscosity. The CO2 hydrate slurry exhibits obvious shear-thickening behavior throughout both the low-shear-rate region and high-shear-rate region. This result shows an apparent difference to our previous work [18,19] and Delahaye et al.’s work [1,15]. In the previous work, the rheology of the hydrate slurry would transit between the shear-thinning behavior and the shear-thickening behavior according to the variations of the volumetric hydrate concentration or the shear rates [15,18]. Generally, the shear-thinning behavior of the hydrate slurry depends on the inner structure between hydrate particles. If hydrate particles have agglomerated structure between each other, the increase in flow velocity and shear rate can break up the agglomerate of hydrate particles and make the hydrate slurry flow more easily [5,8,11]. Thus, the hydrate slurry would exhibit shear-thinning behavior. The shear-thickening behavior of the hydrate slurry depends on the friction and collision between hydrate particles [31,32]. If hydrate particles are suspended in the hydrate slurry as individual particles or small agglomerates, the flow of the hydrate slurry would induce frictions and collisions between hydrate particles and cause additional energy dissipations and pressure drops [33,34]. Thus, the higher flow velocity of the hydrate slurry causes the higher intensity of frictions and collisions between hydrate particles. The increase in flow velocity results in more increments of pressure drop, and the hydrate slurry starts to exhibit shear-thickening behavior. In this work, since the hydrate slurry is always under high circulation condition, the hydrate particles suspended in the flow loop are fully dispersed. Thus, the CO2 hydrate slurry in the current experiment performs shear-thickening behavior.
The effect of volumetric hydrate concentration and shear rate on the apparent viscosity of the CO2 hydrate slurry is revealed in Figure 6. Because of the shear-thickening behavior of the hydrate slurry, the increase in shear rate results in the apparent viscosity increasing for all hydrate concentration conditions. Moreover, the ascent of volumetric hydrate concentration also induces the apparent increase in viscosity. The increase in volumetric hydrate concentration means the quantity of suspended hydrate particles in the hydrate slurry increases. The probability of frictions and collisions between hydrate particles is improved and the increment of pressure drops induced by the hydrate slurry flowing is enhanced. Thus, the increase in volumetric hydrate concentration drives the hydrate slurry performing more non-Newtonian.
From Equation (4) and Figure 3, the non-Newtonian fluid index (n index) under varied volumetric hydrate concentration condition can be obtained. Figure 7 exhibits variations of the n index of the CO2 hydrate slurry along the volumetric hydrate concentration. The n index increases exponentially with the increase in volumetric hydrate concentration. Different from our previous work, the trend in the n index is discontinuity along the volumetric hydrate concentration, increasing from 1.4 to 17.2 vol%. The regressed correlations of the n index are given in Equation (6):
n = 1.9239 α h y d r 0.048 α h y d r < 8.6 % v o l n = 2.0753 α h y d r 0.1224 8.6 % v o l α h y d r
The discontinuity of the n index of the CO2 hydrate slurry indicates the change of rheological property of the hydrate slurry. The critical volumetric hydrate concentration is 8.6 vol% in Figure 7. Although no solid evidence was found in the experiment to explain this discontinuity, the existing theory does. The main reason for this is that the CO2 hydrate slurry is not fully converted to the non-Newtonian fluid when the volumetric hydrate concentration is less than 8.6 vol%. A similar conclusion is received in the CH4 hydrate slurry [18]. Under 8.6 vol%, the CO2 hydrate slurry starts to exhibit non-Newtonian fluid behavior due to hydrate particles but is still not fully converted where the flow pattern is located in the transient zone. At this stage, the pressure drop of the hydrate slurry is mainly contributed by turbulent flow and its rheology. Above 8.6 vol%, the CO2 hydrate slurry is fully converted to the non-Newtonian fluid, where flow pattern is a normally laminar flow. At this stage, the pressure drop of hydrate slurry is dominated by the rheology. However, in this work, once the hydrate forms in the flow loop, the CO2 hydrate slurry is considered a non-Newtonian fluid, and the same assumption is made in previous work. This assumption takes an advantage for modeling the rheology of the hydrate slurry easily and has good accuracy within the experimental range. Local engineers can utilize this model easily and conveniently.
From Equation (3) and Figure 4, the consistency factor (K) factor under varied volumetric hydrate concentration condition can be plotted in Figure 8. Based on the critical volumetric hydrate concentration, the correlations of the K factor are also regressed differently as shown in Figure 7. The K factor increases exponentially as the volumetric hydrate concentration increases. Moreover, by substituting Equations (6) and (7) into simplified Equation (5), the expression of the apparent viscosity of the CO2 hydrate slurry is obtained, as shown in Equation (8).
K = 0.0034 α h y d r 0.1763 α h y d r < 8.6 % v o l K = 0.0144 α h y d r 0.4381 8.6 % v o l α h y d r
μ a p p = 0.0034 α h y d r 0.1763 γ H S 1.9239 α h y d r 0.048 1 α h y d r < 8.6 % v o l μ a p p = 0.0144 α h y d r 0.4381 γ H S 2.0753 α h y d r 0.1224 1 8.6 % v o l α h y d r

5. Validation

With Equation (8), the apparent viscosities of the CO2 hydrate slurry are well predicted. Comparisons between experimental data and predicted values are shown in Figure 9. Most of the predicted apparent viscosities exhibit good accuracies and have discrepancies within ±20%. The average discrepancies of the predicted values are within 16.3%. The developed correlation of the apparent viscosity satisfies the application in the industrial condition.
Generally, the developed model shows a good performance at the low apparent viscosity condition, but the absolute discrepancy between the predicted value and the experimental value is higher at the high apparent viscosity condition. The high apparent viscosity means a high shear rate and flow velocity of the hydrate slurry, where the flow regime tends to transit from the laminar to the turbulent. The reliability of the developed model decreases gradually.
Moreover, the developed correlation is empirical and correlated from experimental data. The accuracy of the correlation will have good performance within the range of the experiment. However, the experimental range in the current work covers most parts of industrial conditions in the aspects of flow velocity, void fraction and volumetric concentration. Additionally, Figure 10 gives the distribution of Reynolds numbers along the shear rate of the hydrate slurry under different hydrate concentration conditions. It shows that all Reynolds number are less than 4500 and the most of them are less than 2100, where the hydrate slurry is under laminar flow condition. The increase in hydrate concentrations results in the reduction of Reynolds numbers. In the future, the inner mechanism is studied necessarily, where the CO2 hydrate slurry exhibits shear-thickening behavior. The application scope of the rheology model needs to improve by developing a mechanistic rheological model.

6. Conclusions

Considering that CO2 is one of reservoir fluid components, CO2 can amplify the risk of hydrate plugging in the production line of deep-water oil and gas wellbore under a low-temperature and high-pressure environment. The rheology of the CO2 hydrate slurry is experimentally studied by a new flow loop with an inner diameter of 50 mm. The CO2 hydrate slurry is a power-law fluid and performs obvious shear-thickening behavior for the shear rate ranging from 40 to 590 s−1, and volumetric hydrate concentration ranging from 1.4 to 17.2 vol%. The increases in shear rate and volumetric hydrate concentration both result in the apparent viscosity of the hydrate slurry increasing. The n index and the K factor of the hydrate slurry increase exponentially as the volumetric hydrate concentration increases, where their empirical correlations are obtained from experimental data. An empirical correlation for the rheology of the CO2 hydrate slurry is developed based on the power-law model, with an average discrepancy of 16.3%. This work provides a fundamental theoretical base and a new sight for deep-water oil and gas production, CO2 storage in the deep-water environment and secondary refrigeration.

Author Contributions

Conceptualization, W.F. and W.W.; methodology, W.F. and Z.W.; formal analysis, W.F., B.H. and Z.W.; investigation, W.F. and W.W.; data curation, H.W.; writing—original draft preparation, W.F.; writing—review and editing, W.F., W.W. and H.W.; supervision, B.H.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) Youth Fund (52104047), CNPC Innovation Fund (2021DQ02—1005), the Natural Science Foundation of Jiangsu Province Youth Fund (BK20210507), and the Independent Research Project of State Key Laboratory of Coal Resources and Safe Mining, CUMT (SKLCRSM22X002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is from the experiments we conducted in this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Delahaye, A.; Fournaison, L.; Marinhas, S.; Martínez, M.C. Rheological study of CO2 hydrate slurry in a dynamic loop applied to secondary refrigeration. Chem. Eng. Sci. 2008, 63, 3551–3559. [Google Scholar]
  2. Sun, B.; Song, R.; Wang, Z. Overflow behaviors of natural gas kick well with high content of H2S gas. China Univ. Pet. 2012, 36, 73–79. [Google Scholar]
  3. Sun, B.; Sun, X.; Wang, Z.; Chen, Y. Effects of phase transition on gas kick migration in deepwater horizontal drilling. J. Nat. Gas Sci. Eng. 2017, 46, 710–729. [Google Scholar]
  4. Clain, P.; Delahaye, A.; Fournaison, L.; Mayoufi, N.; Dalmazzone, D.; Fürst, W. Rheological properties of tetra-n-butylphosphonium bromide hydrate slurry flow. Chem. Eng. J. 2012, 193–194, 112–122. [Google Scholar] [CrossRef]
  5. Peng, B.Z.; Chen, J.; Sun, C.Y.; Dandekar, A.; Guo, S.H.; Liu, B.; Mu, L.; Yang, L.A.; Li, W.Z.; Chen, G.J. Flow characteristics and morphology of hydrate slurry formed from (natural gas+diesel oil/condensate oil+water) system containing anti-agglomerant. Chem. Eng. Sci. 2012, 84, 333–344. [Google Scholar]
  6. Fu, W.; Wang, Z.; Zhang, J.; Cao, Y.; Sun, B. Investigation of rheological properties of methane hydrate slurry with carboxmethylcellulose. J. Pet. Sci. Eng. 2020, 184, 106504. [Google Scholar]
  7. Sinquin, A.; Palermo, T.; Peysson, Y. Rheological and Flow Properties of Gas Hydrate Suspensions. Oil Gas Sci. Technol. 2004, 59, 41–57. [Google Scholar] [CrossRef]
  8. Yan, K.-L.; Sun, C.-Y.; Chen, J.; Chen, L.-T.; Shen, D.-J.; Liu, B.; Jia, M.-L.; Niu, M.; Lv, Y.-N.; Li, N.; et al. Flow characteristics and rheological properties of natural gas hydrate slurry in the presence of anti-agglomerant in a flow loop apparatus. Chem. Eng. Sci. 2014, 106, 99–108. [Google Scholar]
  9. Qin, Y.; Zachary, M.A.; Paul, F.P.; Michael, L.J.; Eric, F.M. High Pressure Rheological Measurements of Gas Hydrate-in-Oil Slurries. J. Non-Newton. Fluid Mech. 2017, 248, 40–49. [Google Scholar]
  10. Liu, Y.; Lv, X.; Shi, B.; Zhou, S.; Lei, Y.; Yu, P.; Chen, Y.; Song, S.; Ma, Q.; Gong, J.; et al. Rheological study of low wax content hydrate slurries considering phase interactions. J. Nat. Gas Sci. Eng. 2021, 94, 104106. [Google Scholar]
  11. Bohui, S.; Shangfei, S.; Yuchuan, C.; Xu, D.; Qingyun, L.; Shunkang, F.; Lihao, L.; Jinhao, S.; Junpeng, J.; Haotian, L.; et al. Status of Natural Gas Hydrate Flow Assurance Research in China: A Review. Energy Fuels 2021, 35, 3611–3658. [Google Scholar]
  12. Andersson, V.; Gudmundsson, J.S. Flow Properties of Hydrate-in-Water Slurries. Ann. N. Y. Acad. Sci. 2000, 912, 322–329. [Google Scholar] [CrossRef]
  13. Darbouret, M.; Cournil, M.; Herri, J.-M. Rheological study of TBAB hydrate slurries as secondary two-phase refrigerants. Int. J. Refrig. 2005, 28, 663–671. [Google Scholar] [CrossRef]
  14. Wang, W.; Fan, S.; Liang, D.; Yang, X. Experimental study on flow characters of CH3CCl2F hydrate slurry. Int. J. Refrig. 2008, 31, 371–378. [Google Scholar] [CrossRef]
  15. Fournaison, L.; Jerbi, S.; Delahaye, A. Rheological Properties of CO2 Hydrate Slurry Flow in the Presence of Additives. Ind. Eng. Chem. Res. 2011, 50, 8344–8353. [Google Scholar]
  16. Hashimoto, S.; Kawamura, K.; Ito, H.; Nobeoka, M.; Ohgaki, K.; Inoue, Y. Rheological Study of TNBA Slaly Semi-Clathrate Hydrate Slurries. In Proceedings of the 7th International Conference on Gas Hydrates (ICGH 2011), Edinburgh, Scotland, UK, 17–21 July 2011. [Google Scholar]
  17. Jerbi, S.; Delahaye, A.; Oignet, J.; Fournaison, L.; Haberschill, P. Rheological properties of CO2 hydrate slurry produced in a stirred tank reactor and a secondary refrigeration loop. Int. J. Refrig. 2013, 36, 1294–1301. [Google Scholar] [CrossRef]
  18. Baojiang, S.; Weiqi, F.; Zhiyuan, W.; Jianchun, X.; Litao, C.; Jintang, W.; Jianbo, Z. Characterizing the Rheology of Methane Hydrate Slurry in a Horizontal Water-Continuous System. SPE J. 2020, 25, 1026–1041. [Google Scholar]
  19. Weiqi, F.; Zhiyuan, W.; Baojiang, S.; Jianchun, X.; Litao, C.; Xuerui, W. Rheological Properties of Methane Hydrate Slurry in the Presence of Xanthan Gum. SPE J. 2020, 25, 2341–2352. [Google Scholar]
  20. Cao, X.; Yang, K.; Bian, J. Investigation of CO2 hydrate slurry flow characteristics with particle dissociation for carbon storage and transportation. Process Saf. Environ. Prot. 2021, 152, 427–440. [Google Scholar]
  21. Fu, W.; Chen, B.; Zhang, K.; Liu, J.; Sun, X.; Huang, B.; Sun, B. Rheological Behavior of Hydrate Slurry with Xanthan Gum and Carboxmethylcellulose under High Shear Rate Conditions. Energy Fuels 2022, 36, 3169–3183. [Google Scholar] [CrossRef]
  22. Webb, E.B.; Rensing, P.J.; Koh, C.A.; Sloan, E.D.; Sum, A.K.; Liberatore, M.W. High-Pressure Rheology of Hydrate Slurries Formed from Water-inOil Emulsions. Energy Fuels 2012, 26, 3504–3509. [Google Scholar]
  23. Webb, E.B.; Koh, C.A.; Liberatore, M.W. High Pressure Rheology of Hydrate Slurries Formed from Water-in-Mineral Oil Emulsions. Ind. Eng. Chem. Res. 2014, 53, 6998–7007. [Google Scholar]
  24. Liu, Z.; Song, Y.; Liu, W.; Liu, R.; Lang, C.; Li, Y. Rheology of methane hydrate slurries formed from water-in-oil emulsion with different surfactants concentrations. Fuel 2020, 275, 117961. [Google Scholar]
  25. Clain, P.; Ndoye, F.T.; Delahaye, A.; Fournaison, L.; Lin, W.; Dalmazzone, D. Particle size distribution of TBPB hydrates by focused beam reflectance measurement (FBRM) for secondary refrigeration application. Int. J. Refrig. 2015, 50, 19–31. [Google Scholar] [CrossRef]
  26. Muhlstedt, G.; Galdino, J.F.; Andrade, D.E.V.; Negrao, C.O.R. Rheological Properties of Hydrate Slurry in the Water-Based Drilling Fluid. Energy Fuels 2021, 35, 10488–10497. [Google Scholar]
  27. Sakurai, S.; Hoskin, B.; Choi, J.; Nonoue, T.; May, E.F.; Kumar, A.; Norris, B.W.E.; Aman, Z.M. Investigating hydrate formation rate and the viscosity of hydrate slurries in water-dominant flow: Flowloop experiments and modelling. Fuel 2021, 292, 120193. [Google Scholar]
  28. Lv, X.; Zhang, J.; Zuo, J.; Zhao, D.; Liu, Y.; Zhou, S.; Du, H.; Song, S. Study on the Formation Characteristics of CO2 Hydrate and the Rheological Properties of Slurry in a Flow System Containing Surfactants. ACS Omega 2022, 7, 2444–2457. [Google Scholar] [CrossRef]
  29. Metzner, A.B.; Reed, J.C. Flow of non-newtonian fluids—correlation of the laminar, transition, and turbulent-flow regions. Aiche J. 1955, 1, 434–440. [Google Scholar] [CrossRef]
  30. Dodge, D.W.; Metzner, A.B. Turbulent flow of non-newtonian systems. Aiche J. 1959, 5, 189–204. [Google Scholar] [CrossRef]
  31. Metzner, A.B.; Whitlock, M. Flow Behavior of Concentrated (Dilatant) Suspensions. Trans. Soc. Rheol. 1958, 2, 239. [Google Scholar] [CrossRef]
  32. Griskey, R.G.; Green, R.G. Flow of dilatant (shear-thickening) fluids. Aiche J. 1971, 17, 725–728. [Google Scholar] [CrossRef]
  33. Fu, W.; Yu, J.; Xu, Y.; Wang, C.; Huang, B.; Sun, B. A Pressure Drop Prediction Model for Hydrate Slurry Based on Energy Dissipation under Laminar Flow Condition. SPE J. 2022, 27, 209586. [Google Scholar]
  34. Fu, W.; Yu, J.; Xiao, Y.; Wang, C.; Huang, B.; Sun, B. A pressure drop prediction model for hydrate slurry based on energy dissipation under turbulent flow condition. Fuel 2022, 311, 122188. [Google Scholar] [CrossRef]
Figure 1. Introduction of high-pressure and low-pressure multiphase flow loop. (a) Image of flow loop. (b) Schematic of flow loop.
Figure 1. Introduction of high-pressure and low-pressure multiphase flow loop. (a) Image of flow loop. (b) Schematic of flow loop.
Jmse 10 01224 g001
Figure 2. Variations of key flow parameters in the experiment. (a) Variations of pressure drop and hydrate slurry velocity in the experiment. (b) Variations of pressure drop and volumetric hydrate fraction in the experiment.
Figure 2. Variations of key flow parameters in the experiment. (a) Variations of pressure drop and hydrate slurry velocity in the experiment. (b) Variations of pressure drop and volumetric hydrate fraction in the experiment.
Jmse 10 01224 g002
Figure 3. Representation of Ln(DΔP/(4 L)) and Ln(8vHS/D) for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Figure 3. Representation of Ln(DΔP/(4 L)) and Ln(8vHS/D) for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Jmse 10 01224 g003
Figure 4. Representation of shear stress and shear stress for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Figure 4. Representation of shear stress and shear stress for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Jmse 10 01224 g004
Figure 5. Representation of shear rate and apparent viscosity for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Figure 5. Representation of shear rate and apparent viscosity for CO2 hydrate concentrations from 1.39 to 17.19 vol%.
Jmse 10 01224 g005
Figure 6. Evolution of apparent viscosity of CO2 hydrate slurry along with hydrate concentration increasing under shear rate at 100, 200 and 300 s−1.
Figure 6. Evolution of apparent viscosity of CO2 hydrate slurry along with hydrate concentration increasing under shear rate at 100, 200 and 300 s−1.
Jmse 10 01224 g006
Figure 7. Evolution of the non-Newtonian index n of a function of volumetric hydrate fraction for CO2 hydrate slurry.
Figure 7. Evolution of the non-Newtonian index n of a function of volumetric hydrate fraction for CO2 hydrate slurry.
Jmse 10 01224 g007
Figure 8. Evolution of the consistency factor K of a function of volumetric hydrate fraction for CO2 hydrate slurry.
Figure 8. Evolution of the consistency factor K of a function of volumetric hydrate fraction for CO2 hydrate slurry.
Jmse 10 01224 g008
Figure 9. Comparison between experimental data and predicted apparent viscosity.
Figure 9. Comparison between experimental data and predicted apparent viscosity.
Jmse 10 01224 g009
Figure 10. Reynold numbers along shear rate under different hydrate concentration conditions.
Figure 10. Reynold numbers along shear rate under different hydrate concentration conditions.
Jmse 10 01224 g010
Table 1. Rheology of hydrate slurries in aqueous phase [1,4,6,12,13,14,15,18,19,21].
Table 1. Rheology of hydrate slurries in aqueous phase [1,4,6,12,13,14,15,18,19,21].
ReferencesHydrate Slurry TypesExperimental MethodHydrate Fraction (vol%)Rheological BehaviorApparent Viscosity (mPa s)Shear Rate (s−1)
Andersson and Gudmundsson [12]NG hydrateOstwald viscosimeter
(Di = 8 mm and 20 mm)
1~11Bingham plastic1~3.50~800
Darbouret et al. [13]TBAB hydrateOstwald viscosimeter
(Di = 8 mm and 20 mm)
4~53Bingham plastic8~1705~270
Fukushima et al. [15]TBAB hydrateOstwald viscosimeter26–35Pseudoplastic50~148/
Delahaye et al. [1]CO2 hydrateOstwald viscosimeter
(Di = 16 mm)
4~20αhydr < 5 vol%: Dilatant
5 vol% ≤ αhydr < 10 vol%: HB-type (dilatant trend)
αhydr = 10 vol%: Bingham plastic
10 vol% < αhydr: HB-type (pseudoplastic trend)
3.8~42.2500~1000
Delahaye et al. [15]CO2 hydrate with additivesOstwald viscosimeter (Di = 5 mm)4~104 vol% ≤ αhydr < 10 vol%: Newtonian (Caflon)7~15 at 600 s−1500~750
Wang et al. [14]CH3CCl2F hydrateOstwald viscosimeter
(Di = 42 mm)
10~68Dilatant1.1~1.7/
Clain et al. [4]TBPB hydrateOstwald viscosimeter (Di = 8 mm)0~28.2Pseudoplastic4–41130~700
Sun et al. [18]CH4 hydrateOstwald viscosimeter (Di = 25.4 mm)6~1120 ≤ γHS < 180 s−1: Pseudoplastic
180 s−1 ≤ γHS < 300 s−1: HB-type (transient region)
300 s−1 ≤ γHS: Dilatant
3~5020~700
Fu et al. [6,19,21]CH4 hydrate with XG and CMCOstwald viscosimeter (Di = 25.4 mm)1~9Both pseudoplastic and dilatant depending on shear rate and additive concentration12~50010~750
Fu et al.
(present work)
CO2 hydrateOstwald viscosimeter (Di = 50 mm)1.4~17.2Dilatant12~52040~590
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, W.; Wei, W.; Wang, H.; Huang, B.; Wang, Z. Study on the Rheology of CO2 Hydrate Slurry by Using the Capillary Method. J. Mar. Sci. Eng. 2022, 10, 1224. https://doi.org/10.3390/jmse10091224

AMA Style

Fu W, Wei W, Wang H, Huang B, Wang Z. Study on the Rheology of CO2 Hydrate Slurry by Using the Capillary Method. Journal of Marine Science and Engineering. 2022; 10(9):1224. https://doi.org/10.3390/jmse10091224

Chicago/Turabian Style

Fu, Weiqi, Wei Wei, Hao Wang, Bingxiang Huang, and Zhiyuan Wang. 2022. "Study on the Rheology of CO2 Hydrate Slurry by Using the Capillary Method" Journal of Marine Science and Engineering 10, no. 9: 1224. https://doi.org/10.3390/jmse10091224

APA Style

Fu, W., Wei, W., Wang, H., Huang, B., & Wang, Z. (2022). Study on the Rheology of CO2 Hydrate Slurry by Using the Capillary Method. Journal of Marine Science and Engineering, 10(9), 1224. https://doi.org/10.3390/jmse10091224

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop