Research Article  Open Access
Investigation on the Transition Criterion of Smooth Stratified Flow to Other Flow Patterns for GasHydrate Slurry Flow
Abstract
A stability criterion for gashydrate slurry stratified flow was developed. The model was based on onedimensional gasliquid twofluid model and perturbation method, considering unstable factors including shear stress, gravity, and surface tension. In addition, mass transfer between gas and liquid phase caused by hydrate formation was taken into account by implementing an inward and outward natural gas hydrates growth shell model for waterinoil emulsion. A series of gashydrate slurry flow experiments were carried out in a highpressure (>10 MPa) horizontal flow loop. The transition criterion of smooth stratified flow to other flow patterns for gashydrate slurry flow was established and validated and combined with experimental data at different water cuts. Meanwhile, parameters of this stability criterion were defined. This stability criterion was proved to be efficient for predicting the transition from smooth to nonsmooth stratified flow for gashydrate slurry.
1. Introduction
Gas hydrates are icelike crystals formed by inclosing gas molecules (guests) in clathrates of water molecules (host) under high pressure and low temperature [1]. Gas hydrates are known to cause blockage in oil and gas transportation pipelines [2]. With the tendency of offshore oil industry towards deeper water, gathering systems are confronting unprecedented flow assurance challenges on hydrate formation and plugging. Safety has become an unparalleled problem in offshore oil and gas facilities under undulating terrain condition, especially for highpressure and low temperature operation conditions owing to the risk of hydrate blockage.
Conventionally, chemicalbased injection and insulation are two major techniques to prevent hydrate plugging in offshore production pipelines. However, these two techniques are of high capital expenditure and technical limitations [1, 3]. Recently, hydrate slurry technique has turned out to be a new option to avoid hydrate blocking. In this method, hydrates are not prevented from formation but are made transportable suspending in liquid phase under complex multiphase flow condition [4–9]. Obviously, multiphase flow analysis is an important factor for implementing hydrate slurry technique in deep water petroleum explorations and developments.
Flow pattern is the key issue in characterizing multiphase flow. There are two generally recognized methods for flow pattern determination: plotting flow pattern map according to the experimental data [10–13] and establishing flow pattern transition criterion based on flow mechanism [14–16]. These two methods have been widely implemented in gasliquid twophase flow. However, the flow would be much too complicated, when hydrate forms in the liquid phase.
Though lots of investigations have been made on thermodynamic and kinetic of gashydrate formation and decomposition [17–24], study on multiphase flow of gashydrate slurry is far from sufficient. Tajima et al. [25] observed CO_{2} hydrate slurry flow patterns in CO_{2}water twophase flow using different types of static mixers. Zerpa et al. [26] and Joshi et al. [27] proposed a conceptual model for hydrate plug formation mechanism in waterdominated system, where gasliquid and liquidsolid flow pattern were defined and characteristics of flow pattern transition with increasing of hydrates volume fraction were studied. Zhao [28] and Gong et al. [7, 8] numerically simulated gashydrate slurry stratified flow, and Kwon et al. [29] simulated the oilgaswaterhydrate annularmist flow. Ding et al. [30] reported the difference of the flow patterns with and without hydrates. And, it was concluded that hydrates formation resulted in the decrease of stratified smooth flow area on flow map and less stability of stratified smooth flow. However, little attention has been paid to the flow pattern transition criterion of gashydrate slurry flow, especially the stability criterion of stratified flow.
In this paper, gashydrate slurry interface stability was analyzed by deriving dispersion function of the interface wave, and then the transition criterion from smooth stratified flow to other flow patterns was obtained. Thereafter, a set of gashydrate slurry flow experiments in highpressure (>10 MPa) horizontal flow loop were carried out to investigate the flow pattern characteristics of gashydrate slurry flow. The transition criterion of smooth stratified flow to other flow patterns for gashydrate slurry flow was established and validated, combined with experimental data at different water cuts.
2. Model Development
A onedimensional twofluid model describing gashydrate slurry stratified flow was developed. The schematic description is shown in Figure 1. In this work, the slurry phase was simplified as incompressible. Moreover, work exchange between the fluid and the environment as well as temperature change caused by hydrates formation was ignored.
In Figure 1, and represent the height of gas phase and hydrate slurry phase, and represent the cross section area of gas phase and hydrate slurry phase, and represent velocity of gas phase and hydrate slurry phase, , , and are shear friction at interface of gas phasepipe wall, hydrate slurrypipe wall, and gashydrate slurry phase, , , and represent wetted perimeter of gas phase, hydrate slurry phase, and gashydrate slurry interface, , , and represent pressure of gas phase, hydrate slurry phase, and pressure at the interface, and represent relative pressure at the interface of gas phase and hydrate slurry phase, represents angle of inclination, represents distance to the infinitesimal segment, and and represent axial and radial direction of the pipe. In addition, the cross section geometry can be found in the supporting information document in Supplementary Material available online at https://doi.org/10.1155/2017/9846507 (cross section geometry).
2.1. Flow Analysis
2.1.1. Continuity Equation
For the gashydrate slurry stratified flow shown in Figure 1, continuity equation for each phase can be derived as follows:where and are densities of gas and hydrate slurry phase, represents overall cross section area, and represents mass transfer rate, which is sectional average weight of gas transferred from gas to liquid per second for hydrate formation. The mass transfer rate, , can be estimated using the inward and outward gashydrate growth shell model [23], since this model was developed to describe the hydrates formation in waterinoil emulsion system. The process of hydrates nucleating, forming, and growing as shell covering the water droplets suspended in the continuous liquid phase can be simulated. And, this inward and outward gashydrate growth shell model is proved suitable and preferred to be applied in the hydrates slurry flow for the calculation of hydrates fraction and gas consumption. And, the application of this hydrate growth shell model is available in the supporting information (Hydrate Growth Model Implementation).
By introducing slurry phase height , which is a key parameter for stratified flow, (1a) and (1b) can be transformed into (2a) and (2b) as function of .where .
2.1.2. Momentum Balance Equation
Momentum balance for each phase can be described as follows:Assuming the interface to be flat in radial direction [31] and ignoring the influence of shorter waves on flow stability [32], the pressure difference at the interface can be approximated aswhere represents interfacial tension and is considered constant in the control volume.
Interfacial shear stress consists of smooth interface shear stress and instability contribution. The following method, (5a) and (5b), is used in this work [33] to determine interfacial shear stress , introducing a correlated coefficient where , , and refer to corresponding density, velocity, and viscosity of the phase with higher velocity and subscript represents G or SL according to the specific situation. Interfacial frication factor can be calculated after obtaining empirical parameters and through experiments.
Substituting (4) and (5a) and (5b) into (3a) and (3b), the key equation of momentum balance relationship can be obtained as follows:where and refer to liquid holdup and void fraction, respectively. For steady state, , (6a), (6b), and (6c) reduce to steadystate momentum equation.
Gas and hydrate slurry shear stresses at the pipe wall can be expressed as (7a) and (7b) by introducing fanning friction factor. For gas phase, the following methods as (8a) and (8b) are applied for friction coefficient .For hydrate slurry phase, a friction coefficient method, (9) [34], is used, considering the influence of liquid flow as well as hydrate formation and agglomeration.where and are hydrate and liquid densities, and are hydrate volume fraction and its maximum value, is hydrate particle diameter, is wetted perimeter diameter of hydrate slurry phase, and is the absolute roughness of the pipe wall.
Substituting (7a) and (7b) into (6c), (10) can be obtained.where the first term is positive when and third term is negative when , and vice versa [35].
2.2. Stability Analysis
2.2.1. Flow Stability Equation
Referring to the method used in interface stability analysis, where the steadystate contribution and perturbation contribution are separated, , , and can be rewritten as combinations of steadystate part and perturbed part shown inBased on (11a), (11b), and (11c), continuous equations (2a) and (2b) and momentum equations (6a), (6b), and (6c) can be linearized and simplified. The following fifthorder linear equation can be obtained to describe onedimensional steadystate gashydrate slurry stratified flow. Characterize the perturbation of liquid holdup in exponential form as in the following [36]: where , , and are parameters in wave equation.
Substituting (13) into (12), (14a), (14b), (14c), (14d), (14e), (14f), (14g), (14h), (14i), and (14j) can be obtained to describe the stratified smooth flow stability of gashydrate slurry stratified flow.where following parameters are defined to simplify expression (14a):
2.2.2. Stability Criterion
Analytically solve quadratic equations (14a), (14b), (14c), (14d), (14e), (14f), (14g), (14h), (14i), and (14j), and roots can be expressed as follows:Substituting , (14a) can be split into imaginary part and real part as follows: Imaginary part: Real part:It could be found in (16a) and (16b) that if , perturbation contribution of liquid holdup increases exponentially with time, resulting in flow instability. With approaching zero, stability criterion for gashydrate slurry stratified smooth flow can be derived as follows:Analyzing (17), two critical wave velocities at the inception of instability are introduced as (18a) and (18b). is the critical wave velocity considering the viscous effect, while is the one without the influence of viscosity [33, 37, 38].Substituting (14a), (14b), (14c), (14d), (14e), (14f), (14g), (14h), (14i), and (14j) and (18a) and (18b) into (17), the stability criterion for gashydrate slurry stratified smooth flow can be obtained as follows:In (19), the first term corresponds to the viscous shear stress that amplifies any disturbance on the surface; the second and third terms are stability criterion disregarding viscous effects, also known as the onedimensional interface IKH instability criterion; the fourth term stands for the surface tension contribution and depends mainly on wave variable ; the last term amplifies instability brought by other factors and is independent of the direction of relative velocity as well as shear stress between the two phases.
Neglecting interfacial tension, viscous shear stress, and interface shear stress, (20) can be simplified.Considering the large density difference between gas and hydrate slurry phases, (20) can be further simplified as follows, shown in (21). This equation matches the criterion presented by Taitel and Dukler [14] in 1976.Calculation procedure based on the criterion is presented in Scheme 1. For stratified flow boundary determination, superficial velocity of one phase can be calculated when that of the other phase is given. In implanting this model, bisection method is applied in the iteration calculation. Scheme 1 shows the flowchart to calculate liquid superficial velocity at a given gas superficial velocity , while at given can be obtained through similar procedure.
3. Experiments
To analyze stability transition criterion for stratified smooth flow using (19), parameters including , , , and need to be determined. In this paper, these parameters were obtained based on the experiments at different water cuts using highpressure hydrates slurry flow loop.
3.1. Material and Facility
A doublepass highpressure horizontal flow loop shown in Figure 2 was designed and established in China University of Petroleum, Beijing, supported by CNOOC [24, 39–42]. This flow loop is 30 m long, with inner diameter of 1 inch and pressure exceeding 10 MPa. Test section temperature can be maintained from −20°C to 80°C by circulating waterglycol blend with the countercurrent flowing surrounding it in a 2inch jacket. Fluid inside the loop can be observed through highpressure windows at two different locations. Gas and liquid can be injected separately by a plunge compressor (2200 m^{3}/h) and a magnetic centrifugal pump (12.0 m^{3}/h). Gas is injected at the inlet of the test section. At the outlet of the test section, gas and liquid flow into an insulated separator and are redirected towards the test section after pressurization.
Thermocouples are placed along the pipe, inside the separator, inside the waterglycol system, and on different gas utilities. A Coriolis flow meter is stalled to measure liquid mixture density and flow rate. Two FM1000 gamma ray densitometers are available to measure the mean density of the multiphase fluid. Differential pressure sensors are installed to follow the evolution of pressure. Rapid data acquisition system permits the detection of quickly occurring events. A Focused Beam Reflectance Measurements (FBRM) probe and Particle Video Microscope (PVM) are installed to capture the evolution of the droplets, bubbles, or solid particles in fluid.
Deionized water, civil natural gas, and −20# diesel (compositions listed in Table 1) were used for the experiments. The dosage of added antiagglomerate was 3 wt% (measuring error ± 0.01 g) to prevent the hydrate crystals from aggregating [43].

3.2. Experiment Procedure
Four groups of experiments at different water cuts (15%, 20%, 25%, and 30%) were carried out with gas, −20# diesel oil, and deionized water. Experimental conditions are listed in Table 2. Specific experimental procedure is detailed as follows:(1)Vacuum the experimental system one hour before injecting the required amount of −20# diesel oil, deionized water, and antiagglomerates.(2)Start and set the temperature controller at 18°C, the magnetic centrifugal pump at 40 Hz, and the control valve at 100%, get the liquid phase well circulated for least 5 hours to form waterinoil emulation, and then open FBRM to monitor changes of the particles in the fluid.(3)Open the gas injection valve and pressure the system up to experimental set point, after the temperature of the system stays constant at 18°C and particles size observed using FBRM becomes stable.(4)Set the system temperature to experimental set point and cool the flow loop; hydrate would form as the temperature drops below the hydrate equilibrium temperate at the system pressure.(5)Maintain temperature and pressure at experimental set point for at least 5 hours to ensure the fully hydrate formation, adjust compressor inlet valve and pump speed to attain a gas flow rate within 20~195 kg/h and liquid flow rate within 75~860 kg/h, and then observe the flow pattern and record the data; for each data point presented in this work, flow rate of both gas and slurry phase were controlled and flow pattern was observed when the flow rates, pressure, and temperature were stabilized.(6)Increase the system temperature, stop the compressor and pump, evacuate the residual gas, discharge the experiment liquid, clean and flush the flow loop with compressed air, and replace it with nitrogen.(7)Repeat the procedure at different water cuts.
 
: number of experimental data points. 
4. Results and Discussions
4.1. Experiment Results
In the experiments, 293 smooth and nonsmooth stratified flow pattern data points at different water cuts were obtained under steadystate conditions (Tables S.1–S.4), which was visually observed through the sight glass of the loop. Four types of flow patterns were found, including two typical flow patterns (stratified smooth flow and slug flow) and two transitional flow patterns (stratified wavy flow and short slug flow). Distinct gasliquid interface can be observed, and the interface was flat for the stratified smooth flow (Figure 3(a)) and was wavy for stratified wave flow (Figure 3(b)). A liquid film zone and a liquid slug were found intermittently for slug flow (Figure 4), while a relatively longer liquid film and shorter slug length were found for short slug flow.
(a)
(b)
(a)
(b)
(c)
The morphologies, sizes, and distributions of fully formed hydrate particles in the slurry were recorded by PVM as shown in Figure 5. The chord length distribution before and after hydrates formation is displayed in Figure 6. According to experimental observations, it is can be concluded that hydrates particles form as shells covering water droplets in the waterinoil emulsion, and the size of hydrates particles will be enlarged compared with water droplet. In addition, hydrate particle morphologies, sizes, and distributions showed no distinguishable difference at different flow patterns, implying that their influences on flow pattern transition were not signification.
4.2. Model Implementation
Based on experimental conditions listed in Table 2, hydrate growth parameters and gas consumption rate at different water cuts were calculated using the inward and outward hydrate growth shell model [23]. Results presented in Table 3 showed that gas consumption rate increased with water cut. Detailed calculation procedure including water consumption and hydrate fraction is presented in the supporting information (Hydrate Growth Model Implementation). Gasslurry stratified smooth flow region transition boundary at different water cuts was obtained by implementing the calculation procedure presented in Scheme 1. Parameters in the model are obtained, where is 0.01, is −0.2, is 0.01, and is 10.

4.3. Model Evaluation
To compare the transition criterion of smooth stratified flow to other flow patterns for gashydrate slurry flow developed in this work with the classical flow pattern distribution models, flow region boundaries calculated using this model as well as that calculated using TaitelDukler [14] and Mandhane [13] methods are presented in Figure 7.
(a)
(b)
(c)
(d)
As is shown in Figure 7, Mandhane’s model can roughly divide stratified flow from nonstratified flow at four water cuts, yet the results lead to a larger smooth stratified flow region. TaitelDukler model cannot give a satisfying stratified flow transition boundary, due to not considering the effect of hydrate formation on flow stability. Transition boundaries obtained based on the model developed in this work show good agreement with experimental data.
All data points are summarized in Figure 8. It can be seen that no distinctive effect of water cut on flow pattern division is observed, though gas consumption rate differs evidently according to Table 3. This may be due to the similar water consumption rates and hydrate fractions as shown in Figure 9. As calculated flow pattern boundary is primarily affected by fluid viscosity, and hydrate slurry viscosity is further determined by hydrate volume fraction, proximal transition lines are obtained under similar hydrate fraction though at different water cuts.
The pattern experiment data of gashydrate slurry multiphase flow pattern in this work are not effectively enough to obtain more precise correlation parameters or investigate the flow pattern transition mechanism deeply. More experimental and theoretical research should be carried out in the future. Nevertheless, the good agreement with experimental data proved this work to be both feasible and significant.
5. Conclusion
In this work, a onedimension twofluid model for gashydrate slurry stratified flow was developed based on perturbation method, and a stability criterion for smooth stratified flow was proposed. In establishing the criterion, mass transfer between gas and slurry phase caused by hydrate formation was considered, and various mathematical techniques were applied in linearizing the equation sets. As the formation of hydrate would evidently complicate the flow, the influences of shear stress, gravity, surface tension, hydrate formation, and other unstable factors were considered. Compared with the classical twophase gasliquid stratified flow stability criteria proposed by TaitelDukler [14] and Mandhane [13], transition boundary based on this study showed better agreement with experiment results.
Groups of gashydrate slurry multiphase flow experiments were carried out on a most advanced highpressure (>10 MPa) hydrate slurry flow loop in China, and 293 experimental data points of both smooth and nonsmooth stratified flow were obtained. Model parameters were generated in combining the experiment data with the gashydrate slurry stratified flow stability creation developed in this work, and the model was proved applicable in gashydrate slurry multiphase flow numerical simulation and characteristic study.
Nomenclature
:  Cross area of the pipeline, m^{2} 
:  Cross area of gas phase, m^{2} 
:  Cross area of hydrate slurry phase, m^{2} 
:  Viscosity critical wave velocity at the inception of instability, m s^{−1} 
:  Nonviscosity critical wave velocity at the inception of instability, m s^{−1} 
:  Correlated coefficient 
:  Diameter of hydrates particle, m 
:  Diameter of the wetted perimeter of hydrate slurry phase, m 
:  Friction coefficient of interface gasslurry phase 
:  Friction coefficient of interface gas phase at pipe wall 
:  Friction coefficient of interface slurry at pipe wall 
:  Height of gas phase, m 
:  Height of hydrate slurry phase, m 
:  Parameter in wave equation 
:  Pressure of gas phase 
:  Pressure of hydrate slurry phase 
:  Pressure at the interface 
:  Relative pressure at the interface of gas phase 
:  Relative pressure at the interface of hydrate slurry phase 
:  Gas phase fraction of pipe cross area 
:  Hydrate slurry phase fraction of pipe cross area 
:  Reynolds number of gas phase 
:  Wetted perimeter of gas phase, m 
:  Wetted perimeter of interface gasslurry phase, m 
:  Wetted perimeter of hydrate slurry phase, m 
:  Time, s 
:  Velocity of gas phase, m·s^{−1} 
:  Velocity of hydrate slurry phase, m·s^{−1} 
:  Parameter in wave equation 
:  Axial direction of the pipe, m 
:  Distance to the infinitesimal segment, m 
:  Radial direction of the pipe, m 
:  Density of hydrates, kg/m^{3} 
:  Density of oil phase, kg/m^{3} 
:  Hydrate volume fraction 
:  Maximum hydrate volume fraction 
:  Mass transfer rate between gas phase and hydrate slurry phase, kg·s^{−1}·m^{−3} 
:  Angle of inclination, rad 
:  Interfacial tension, N/m 
:  Density of gas phase, kg·m^{−3} 
:  Density of hydrate slurry phase, kg·m^{−3} 
:  Density of faster phase: subscript f represents the faster phase, kg·m^{−3} 
:  Shear friction of interface gas phase at pipe wall, N 
:  Shear friction of interface of gasslurry phase, N 
:  Shear friction of interface of slurry at pipe wall, N 
:  Absolute roughness of the pipe wall, m 
:  Friction calculated empirical parameter 
:  Friction calculated empirical parameter 
:  Parameter in wave equation. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (51306208, 51534007, and 51274218), National Science and Technology Major Project of China (2016ZX05028004001), National Key Research and Development Plan of China (SQ2016YFSF010222), and Science Foundation of China University of PetroleumBeijing (2462015YQ0404 and 201602), which are gratefully acknowledged.
Supplementary Materials
Based on the experiments carried out in the highpressure hydrate flow loop, 293 smooth and nonsmooth stratified flow pattern data points at water cuts of 15%, 20%, 25%, and 30% were obtained under steadystate conditions listed in Table S.1–Table S.4.
References
 E. Dendy Sloan and C. A. Koh, Clathrate Hydrates of Natural Gas, Crc Press, 2008.
 E. G. Hammerschmidt, “Formation of gas hydrates in natural gas transmission lines,” Industrial & Engineering Chemistry, vol. 26, no. 8, pp. 851–855, 1934. View at: Publisher Site  Google Scholar
 D. Sloan et al., Natural Gas Hydrates in Flow Assurance, 2010.
 Y. Peysson, “Collision process between particles in the transport of dispersed hydrates in production lines,” in Proceedings of the 5th International Conference on Gas Hydrate, Houston, Tex, USA, 2005. View at: Google Scholar
 O. C. Hernandez, Investigation of Hydrate Slurry Flow in Horizontal Pipelines, 2006.
 E. Dellecase, G. Geraci, L. Barrios, D. Estanga, R. Domingues, and M. Volk, Hydrate plugging or slurry flow: Effect of key variables, 2008.
 J. Gong, B. Shi, and J. Zhao, “Natural gas hydrate shell model in gasslurry pipeline flow,” Journal of Natural Gas Chemistry, vol. 19, no. 3, pp. 261–266, 2010. View at: Publisher Site  Google Scholar
 J. Gong, B. Shi, and W. Wang, “Numerical study the flow characteristics of gashydrate slurry two phase stratified flow,” in Proceeding of the 7th International Conference on Multiphase Flow, Tampa, Fla, USA, 2010. View at: Google Scholar
 M. Gainville, A. Sinquin, and M. Darbouret, Hydrate Slurry Characterisation for Laminar and Turbulent Flows in Pipelines, 2011.
 O. Baker, “Simultaneous flow of oil and gas,” Oil & Gas Journal, vol. 53, p. 185, 1954. View at: Google Scholar
 H. Duns and N. C. Ros, Vertical Flow of Gas and Liquid Mixtures in Wells, 1963.
 G. W. Govier, K. Aziz, and W. R. Schowalter, The Flow of Complex Mixtures in Pipes, Van Nostrand Reinhold Co., 1972.
 J. M. Mandhane, G. A. Gregory, and K. Aziz, “A flow pattern map for gas—liquid flow in horizontal pipes,” International Journal of Multiphase Flow, vol. 1, no. 4, pp. 537–553, 1974. View at: Publisher Site  Google Scholar
 Y. Taitel and A. E. Dukler, “A model for predicting flow regime transitions in horizontal and near horizontal gas‐liquid flow,” AIChE Journal, vol. 22, no. 1, pp. 47–55, 1976. View at: Publisher Site  Google Scholar
 Y. Taitel, D. Bornea, and A. E. Dukler, “Modelling flow pattern transitions for steady upward gas‐liquid flow in vertical tubes,” AIChE Journal, vol. 26, no. 3, pp. 345–354, 1980. View at: Publisher Site  Google Scholar
 D. Barnea, “A unified model for predicting flowpattern transitions for the whole range of pipe inclinations,” International Journal of Multiphase Flow, vol. 13, no. 1, pp. 1–12, 1987. View at: Publisher Site  Google Scholar
 A. L. Ballard and E. D. Sloan Jr., “The next generation of hydrate prediction: part III. Gibbs energy minimization formalism,” Fluid Phase Equilibria, vol. 218, no. 1, pp. 15–31, 2004. View at: Publisher Site  Google Scholar
 L. Ballard and E. D. Sloan Jr., “The next generation of hydrate prediction IV: a comparison of available hydrate prediction programs,” Fluid Phase Equilibria, vol. 216, no. 2, pp. 257–270, 2004. View at: Publisher Site  Google Scholar
 G.J. Chen and T.M. Guo, “A new approach to gas hydrate modelling,” Chemical Engineering Journal, vol. 71, no. 2, pp. 145–151, 1998. View at: Publisher Site  Google Scholar
 P. D. Dholabhai, P. Englezos, and N. Kalogerakis, Kinetics of formation of methane and ethane gas hydrates, 1987.
 P. Englezos et al., “Kinetics of gas hydrate formation from mixtures of methane and ethane,” Chemical Engineering Science, vol. 42, no. 11, pp. 2659–2666, 1987. View at: Publisher Site  Google Scholar
 H. C. Kim, P. R. Bishnoi, R. A. Heidemann, and S. S. H. Rizvi, “Kinetics of methane hydrate decomposition,” Chemical Engineering Science, vol. 42, no. 7, pp. 1645–1653, 1987. View at: Publisher Site  Google Scholar
 B.H. Shi, J. Gong, C.Y. Sun, J.K. Zhao, Y. Ding, and G.J. Chen, “An inward and outward natural gas hydrates growth shell model considering intrinsic kinetics, mass and heat transfer,” Chemical Engineering Journal, vol. 171, no. 3, pp. 1308–1316, 2011. View at: Publisher Site  Google Scholar
 X.S. Li and Y. Zhang, “Study on dissociation behaviors of methane hydrate in porous media based on experiments and fractional dimension shrinkingcore model,” Industrial and Engineering Chemistry Research, vol. 50, no. 13, pp. 13–17, 2011. View at: Google Scholar
 H. Tajima, A. Yamasaki, and F. Kiyono, “Effects of mixing functions of static mixers on the formation of CO2 hydrate from the twophase flow of liquid CO2 and water,” Energy and Fuels, vol. 19, no. 6, pp. 2364–2370, 2005. View at: Publisher Site  Google Scholar
 L. E. Zerpa et al., “Overview of Csmhyk: a Transient Hydrate Formation Model Developed Over a Decade of Hydrate Research,” 2011. View at: Google Scholar
 S. Joshi, I. Rao, L. Zerpa et al., “Understanding hydrate plug formation from high water cut systems with a four inch flowloop,” in Proceedings of the 7th International Conference on Gas Hydrates (ICGH '11), Edinburgh, UK, July 2011. View at: Google Scholar
 J. Zhao, “Study on flow properties of hydrate slurry in multiphase pipeline,” in Oil and Gas Storage and Transportation, China University of Petroleum, Beijing, China, 2009. View at: Google Scholar
 O. Kwon, S. Ryou, and W. Sung, “Numerical modeling study for the analysis of transient flow characteristics of gas, oil, water, and hydrate flow through a pipeline,” Korean Journal of Chemical Engineering, vol. 18, no. 1, pp. 88–93, 2001. View at: Publisher Site  Google Scholar
 L. Ding, B. Shi, X. Lv et al., “Investigation of natural gas hydrate slurry flow properties and flow patterns using a high pressure flow loop,” Chemical Engineering Science, vol. 146, pp. 199–206, 2016. View at: Publisher Site  Google Scholar
 P. Becher, Principles of Colloid and Surface Chemistry, Dekker, 1977.
 D. Barnea and Y. Taitel, “Nonlinear interfacial instability of separated flow,” Chemical Engineering Science, vol. 49, no. 14, pp. 2341–2349, 1994. View at: Publisher Site  Google Scholar
 N. Brauner and D. M. Maron, “Stability of twophase stratified flow as controlled by laminar/turbulent transition,” International Communications in Heat and Mass Transfer, vol. 21, no. 1, pp. 65–74, 1994. View at: Publisher Site  Google Scholar
 P. Yannick et al., Flow of Hydrates Dispersed in Production Lines, 2003.
 D. Deng, “Simulation study on natural gasgas condensate flow characteristics in multiphase transportation pipeline,” in Oil and Gas Storage and Transportation, China University of Petroleum, Beijing, China, 2005. View at: Google Scholar
 P. Drazin and W. Reid, “Hydrodynamic stability,” Journal of Fluid Mechanics, vol. 124, 1981. View at: Google Scholar
 D. Barnea, “On the effect of viscosity on stability of stratified gas—liquid flow—application to flow pattern transition at various pipe inclinations,” Chemical Engineering Science, vol. 46, no. 8, pp. 2123–2131, 1991. View at: Publisher Site  Google Scholar
 D. Barnea and Y. Taitel, “KelvinHelmholtz stability criteria for stratified flow: viscous versus nonviscous (inviscid) approaches,” International Journal of Multiphase Flow, vol. 19, no. 4, pp. 639–649, 1993. View at: Publisher Site  Google Scholar
 W. Li, “Experimental study on multiphase flow mechanism of hydrate slurry in pipes,” in Oil and Gas Storage and Transportation Engineering, China University of Petroleum, Beijing, China, 2012. View at: Google Scholar
 X. Lv et al., “Experimental study on natural gas hydrate slurry flow,” Spe Journal, vol. 2, no. 2, pp. 199–216, 2012. View at: Google Scholar
 X. Lv, D. Yu, W. Li, B. Shi, and J. Gong, “Experimental study on blockage of gas hydrate slurry in a flow loop,” in Proceeding of the 9th International Pipeline Conference (IPC '12), pp. 37–43, 2012. View at: Publisher Site  Google Scholar
 X. Lv, J. Gong, W. Li, Y. Tang, and J. Zhao, “Focusedbeam reflectance method aids hydrate blockage prediction,” Oil and Gas Journal, vol. 111, no. 1, pp. 99–107, 2013. View at: Google Scholar
 C. Y. Sun et al., “The experimental evaluation for controlling hydrate plug in oilgaswater multiphase pipeline,” in Proceeding of the 2nd Chinese National Chemical and Biochemical Engineering Annual Meeting, Beijing, China, 2005. View at: Google Scholar
Copyright
Copyright © 2017 Bohui Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.