• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Deepwater gas kick simulation with consideration of the gas hydrate phase transition*

    2014-06-01 12:30:00WANGZhiyuan王志遠(yuǎn)SUNBaojiang孫寶江
    關(guān)鍵詞:志遠(yuǎn)

    WANG Zhi-yuan (王志遠(yuǎn)), SUN Bao-jiang (孫寶江)

    School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China, E-mail: wangzy1209@126.com

    Deepwater gas kick simulation with consideration of the gas hydrate phase transition*

    WANG Zhi-yuan (王志遠(yuǎn)), SUN Bao-jiang (孫寶江)

    School of Petroleum Engineering, China University of Petroleum (East China), Qingdao 266580, China, E-mail: wangzy1209@126.com

    (Received January 8, 2013, Revised March 26, 2013)

    The hydrate phase transition presents new problems and challenges for the deepwater well control in the drilling processes. A simulation model is built for deepwater gas kicks with consideration of the hydrate phase transition. The model is based on the multiphase flow governing equations and the hydrate phase transition calculation equations. The influence of the hydrate phase transition on the gas kick hydraulics is investigated through numerical simulations. It is shown that the diameter of the gas bubbles can significantly influence the hydrate phase transition effect. The influence of the hydrate phase transition on the gas kick hydraulics increases with the decrease of the average gas bubble diameter. The hydrate phase transition adds a “hidden” nature for the well kick in deepwater and hinders the early detection of the gas kick. The influence of the hydrate phase transition on the gas kick hydraulics is also studied in the case when the hydrate inhibitor is added to the drilling fluid.

    deepwater, gas kick, simulation, gas hydrate, phase transition

    Introduction

    The hydrocarbon resource in the offshore deepwater environment is one of important energy resources. In the deepwater drilling process, when the hydraulic pressure in the wellbore is below the drilling formation pressure, the formation fluids would flow into the wellbore. When the formation fluid is the vapour phase hydrocarbon, it is termed the gas kick in the industry, which can potentially cause a catastrophic well control event. Furthermore, the gas in the well annulus forms solid hydrates in the interval of low temperature and high pressure. The hydrates travel upward with the carrying fluids (drilling mud) along the well annulus. Once the hydrates are out of the hydrate formation region, they are disassociated into the vapour phase. The hydrate phase transition presents new challenges for deepwater drilling[1]. If the hydrates are not handled properly, we might face the well blowout and offshore pollution problems. The deepwater gas kick simulation is a possible way to determine the multiphase flow behaviour with consideration of the hydrate phase transition.

    Nunes et al.[2]established a kick analytical model and iteratively calculated the pressure distribution and the gas-liquid percentages in the well annulus and the choke line at each time step at a constant flow rate. Avelar et al.[3]proposed a set of deepwater kick simulation methods and analyzed the effect of the depth on the well control parameters. Johnny et al.[4]studied the hydrate formation potential and obtained the overcooled pressure of the hydrate formation in the well annulus using an advanced dynamic kick simulator. Wang et al.[5]presented a method for predicting the gas hydrate formation region in the deep-water well annulus. Yapa et al.[6]modelled underwater oil/gas jets and plumes with consideration of the gas hydrate formation and decomposition. Wang and Sun[7]established a deepwater gas kick model and conducted preliminary studies on the influence of the hydrate formation on kick flow parameters.

    Although some progresses were made in deepwater gas kick simulations, few previous studies con-sidered the influence of the hydrate phase transition on the gas kick hydraulics, which is a major problem in deepwater drilling and well control. In this study, a simulation model is established for the gas kick in the deepwater well annulus with consideration of the hydrate phase transition. The influences of the hydrate phase transition on the deepwater gas kick hydraulics are analysed using numerical calculations.

    1. Gas kick model with the gas hydrate phase transition in well annulus

    1.1Multiphase flow governing equations considering the gas hydrate phase transition

    Three components in the well annulus are usually considered, including the gas, the drilling fluid and the cuttings, where the mass and momentum conservation equations are used. However, the gases in the well annulus can form hydrate particles when they are transported into low-temperature and high-pressure regions, particularly near the seafloor, in the deepwater well annulus. The hydrate particles are lighter than the drilling fluids and will travel upward with the drilling fluids and the cuttings, and then disassociate after they leave the hydrate formation region. The hydrate phase transitions make the behaviour of the multiphase flow in the deepwater well differ from that in the well onshore or in the shallow water. In this study, we establish an improved annular multiphase flow governing equations by coupling the hydrate phase transition by adding the hydrate mass conservation equation, the hydrate term in the momentum conservation equation and the enthalpy term for the gas hydrate phase transition in the heat conservation equations.

    The following assumptions are made when the governing equations are derived.

    (1) The seawater temperature field is continuous and the formation temperature field varies linearly with the depth at a constant geothermal gradient from the mudline.

    (2) The water-based drilling fluid is used, and no mass transfer occurs between the natural gas and the drilling fluid.

    (3) The hydrates in the form of solid particles travel upward with the drilling fluid, neglecting the hydrate attachment on the well annulus or riser wall.

    (4) The hydrate particles will not be split into two or more parts during the hydrate formation or disassociation process.

    Based on the above assumptions, the annular multiphase flow governing equations with consideration of the hydrate phase transition, for the mass conservation, the momentum conservation, and the heat conservation, are obtained. The mass conservation equations of gases, hydrates, drilling fluids, and cuttings are as follows: whereρg,ρm,ρHandρcare the densities of the gas, the drilling fluids, the hydrates, and the cuttings, respectively,vg,vm,vHandvcare the velocities of the gas, the drilling fluids, the hydrates and the cuttings in the well annulus, respectively,Eg,Em,EHandEcare the volume fraction of the produced gas, the drilling fluids, the hydrates and the cuttings,Ais the annular cross-sectional area,qgis the gas production rate of unit length of formation,xgis the mass fraction of the gas in the hydrate,rHis the hydrate formation/decomposition rate of unit length in the well annulus, with its positive or negative value representing the formation or disassociation of the hydrates. andqcis the cutting production rate per unit length.

    Equation (2) is the hydrate mass conservation equation. The final terms on the right-hand side of Eqs.(1) and (3) represent the gas and water mass gain/ loss terms due to the hydrate disassociation/formation, respectively.

    The momentum conversation equations of the gases, the hydrates, the drilling fluids, and the cuttings are as follows whereαrepresents the well deviation angle,pfrepresents the annular friction, andprepresents the annular pressure.

    The hydrate phase transition is often accompanied by a heat exchange. The hydrate formation is anexothermic process[8], whereas the hydrate disassociation is an endothermic process[9]. Therefore, the enthalpy term for the gas hydrate phase transition, the third term in the left-hand side of Eq.(6), should be considered in the heat conservation equations. The heat conservation equations in the deepwater well annulus can be expressed as follows

    whereρlis the density of the liquid,wg,wlare the mass flow rates of the gas and liquid phases, respectively,Cg,Clare the specific heats of the gas and liquid phases, respectively,Tei,Ta,Ttare the temperatures of the formation, the well annulus and the drill string, respectively, and ΔHHis the enthalpy of the gas hydrate phase transition, which can be calculated using the Clausius-Clapeyron equation.A',B' are the intermediate coefficients. The expression ofB' is

    The deepwater well annulus section above the seabed is a riser surrounded by seawater and that below the seabed is surrounded by the formation, and the expression for these two sessions are defined as

    wherekeis the coefficient of heat transfer of the hydrate formation,rcois the casing external diameter,Uais the overall heat transfer coefficient of the annular fluid and the formation, andTDis a function of the transient heat transfer[11],his the well depth,Hsis the water depth.

    Because only the drilling fluid flows in the drill string, the temperature field equation can be expressed in a relatively simple manner as follows

    wheretArepresents the inner cross-sectional area of the drilling pipe.

    1.2Auxiliary equations

    Many previous studies considered the auxiliary equations for the multiphase flow in the well annulus, such as the flow pattern transition equations[12], the gas slip velocity equations, the fluid phase behaviour equations, and the flow friction calculation equations. In this study, we focus on the auxiliary equations related to the hydrate phase transition. These equations involve the hydrate phase equilibrium condition, the hydrate formation rate, and the hydrate disassociation rate.

    (1) The hydrate phase equilibrium conditions. The first issue is to identify the conditions under which the hydrates form/decompose. In this paper, the hydrate phase equilibrium conditions refer to the temperature and pressure conditions in the well annulus. A chemical potential balance is established among the aqueous phase, the gas phase, and the hydrate lattice system in the crystal structure of the hydrates. The equations for the hydrate phase equilibrium conditions are based on the thermodynamic equilibrium theory[13].

    where Δμ0is the difference between chemical potentials in the unoccupied lattice and the pure water in the reference state,Ris the gas constant,T0is the reference temperature,THis the hydrate-formation temperature, ΔH0is the difference between enthalpies in the unoccupied lattice and the pure water, ΔCKis the difference between the heat capacities in the unoccupied lattice and the pure water,pHis the hydrateformation pressure, ΔVis the molar volume differe-nce between the unoccupied lattice and the pure water,fwis the fugacity of water in the solution, Pascal,fwris the fugacity of water in the reference state (TH,pH), Pascal,lis the total number of hydrate species,Miis the ratio of the number of typeicavities to the number of water molecules in the hydrate phase,Lis the total number of gas types,θijis the fraction of thei-type cavities occupied by aj-type gas molecule,xwis the mole fraction of water; andywis the water activity coefficient in the solution.

    (2) Calculation of the hydrate formation rate.

    Englezos et al.[14]suggested two gas hydrate formation processes. The first process is the diffusion of the natural gas dissolved in the liquid into the hydrateliquid interface. The second process is the structural framework combination between gas and water molecules at the interface. Englezos proposed a method to calculate the hydrate growth rate from the aforementioned processes. The main formulas for the calculation of the hydrate formation rate are as follows:

    wherenis the hydrate formation rate,μnis thenth moment and the initialnth moment of the crystal distribution,Ag-lis the gas-liquid interface area,ais the gas-liquid interface area per unit liquid volume,Cwois the initial concentration of water molecules,Dis the mass transfer coefficient for gas in the dispersion film,fgis the gas fugacity,fbis the fugacity of the dissolved gas in water,feqis the three-phase equilibrium fugacity,yLis the thickness of the film at the gas-liquid interface,His the Henry’s law constant,γis the Hatta number, andGis an average growth rate for the computation of the first and second moments,Kis the combined rate parameter of hydrate formation.

    The hydrate formation region is usually near the seabed in the deepwater well annulus. The diameter of the gas bubbles is usually small because of the high pressure in the hydrate formation region. The number of gas bubbles is an important factor for calculating the hydrate formation rate per unit length in an annulus. The gas bubbles are assumed as the volume-equivalent spheres. The number of gas bubbles per unit length in the annulus can be described as

    whereNgis the number of gas bubbles of unit length in the annulus andDgis the initial mean diameter of the gas bubbles as they reach the hydrate formation region.

    Hinze[15]and Winterton[16]calculated the maximum and mean bubble diameters using the following Eqs.(19) and (20) to determine the diameter range of the gas bubbles when they reach the hydrate formation region in the deepwater well annulus.

    wheremaxDis the maximum bubble diameter,σis the surface tension,εis the energy dissipation rate per unit mass[15], andDis the diameter of the annulus.

    The hydrate growth rate for one gas bubble can be calculated using Eqs. (13)-(17). The equation of the hydrate formation rate of unit length in an annulus is as follows

    whereHMis the molar weight of the hydrate.

    (3) Calculation of the hydrate decomposition rate. Kim et al.[17]proposed a hydrate decomposition dynamic model including the following two parts: (1) the deconstruction process of the cage lattice structure of the hydrate particle surface and (2) the desorption process of guest molecules from the surface. The calculation formulas for the hydrate decomposition rate are as follows:

    wherenHis the total moles of CH4contained in hydrate particles,Asis the surface area of one hydrate particle, ΔEis the activation energy, 78.3 kJ/mol for CH4,kdis the hydrate decomposition rate constant.

    (4) The slip velocity of hydrate particles

    The hydrate particles are buoyant due to their low density (e.g., 900 kg/m3 for methane hydrate). The slip velocity of the gas bubbles and the hydrate particles is, according to Zheng and Yapa[18],

    wherevTis the slip velocity of the gas bubble or the hydrate particle,μis the viscosity of drilling fluid,ρmis the drilling fluid density,dis the diameter of the gas bubble or the hy drate particle,Ris the Reynolds number and Cliftprovided a correlation forRas

    whereρdis the density of the gas bubbles or the hydrate particles.

    Equations (24)-(30) can also be applied to calculate the slip velocity between the gas bubble and the hydrate shell whenρdin Eq.(29) is considered as the composite density of the bubbles and the hydrate shell (ρcom). The composite density of the bubbles with the hydrate shell is then

    whereρcomis the composite density of the bubble with hydrate shell,Rhis the initial radius of the gas bubble, andRbis the radius of the gas bubble with hydrate shell, which can be calculated easily from the hydrate growth rate.

    1.3Conditions of temperature field simulation

    (1) Initial conditions. In the drilling process, when the fluid flow state in the annulus changes because of the gas invasion, the initial condition of the transient temperature field is equivalent to the temperature calculated under the steady-state conditions before the gas enters the annulus.

    (2) Boundary conditions. The inlet temperature of the drill string can be measured directly as follows

    The temperature of the liquid in the drill string is similar to the liquid in the annulus at the well bottom, and can be expressed as

    The seawater temperature field can be expressed as

    wherehsis some depth in the seawater,Tseais the seawater temperature at the depthhs,T0is the seawater temperature at the surface, andfTseais a function of the seawater temperature.

    The formation temperature field can be expressed as

    whereTGis the formation temperature at some depth,Kgardis the geothermal gradient.

    Table 1 Basic data

    1.4Conditions of pressure field simulation

    (1) Initial conditions. Before the gas enters into the annulus, there are only the drilling fluids in the annulus. Then,

    whereVscis the apparent velocity of the cuttings,Ccis the velocity distribution factor,Vmmis the velocity of the mixture,Vcris the slippage velocity of the cuttings,Qmis the drilling fluid circulation rate, andQcis the volumetric flow rate of the cuttings.

    (2) Boundary conditions

    wherespis the casing pressure.

    2. Results of the deepwater gas kick simulation with hydrate phase transition

    The HG-X well in the South China Sea is used as an example to analyse the influence of the hydrate phase transition on the annulus flow parameters when a deepwater gas kick occurs. The well data are summarized in Table 1. The kick occurs during the drilling process and the invaded gas is methane.

    2.1Sensitivity study for bubble diameters

    The hydrate formation/disassociation rate is affected by the diameter of bubbles in the gas-liquid mixture because the bubble diameter is related to the gasliquid interface area. The maximum and mean bubble diameters are approximately 0.030 m and 0.009 m in the hydrate formation region, respectively, according to Eqs.(19) and (20). Thus, the gas bubble diameters are assumed in the range of 0.005 m-0.030 m to study the influence of the gas bubble diameter on the gas kick hydraulics, including the gas volume fraction, the pit gain, the bottom hole pressure and the shut-in casing pressure.

    Fig.1 Hydrate formation region in well annulus during drilling

    (1) The hydrate formation region in the annulus during the gas kick period. The hydrate formation region in the annulus can be determined from the hydrate phase equilibrium conditions and using the annulus temperature and pressure field calculation method. Figure 1 shows the hydrate formation region in the annulus during drilling for the case with no hydrate inhibitor added into the drilling fluid. The depth range enclosed by the hydrate phase equilibrium and the annulus temperature curves mark the hydrate formation region. In this case, the hydrate formation region is in the range of 600 m to 1 470 m. The temperature difference between the two curves in the closed region at the same depth is equivalent to the overcooled temperature for the hydrate formation. The hydrate formation region becomes larger when the closed region is longer in the vertical direction. On the other hand, the hydrate formation rate increases when the closed region is wider in the lateral direction.

    Fig.2(a) Gas volume fraction profiles at 95 min during gas kick

    Fig.2(b) Gas fraction profiles at 95 min during gas kick

    (2) The influence of the hydrate phase transition on the gas volume fraction during the gas kick period. Figures 2(a) and 2(b) show the gas fraction profiles in the annulus at the gas production rates of 0.15 Nm3/(MPa·s) and 1.0 Nm3/(MPa·s), respectively. The figures show that the gas volume fractions are different in cases with and without consideration of the hydrate phase transition. The influence of the hydrate phase transition on the gas volume fraction increases with the decrease of the average gas bubble diameter (Dg), because the smaller the gas bubbles diameters are at certain gas volume fractions, the larger the gasliquid interface areas will be, consequently, the hydrate formation rate will be increased. As seen in Fig.2(b), the ratios of the gas volume fractions with and without the hydrate phase transition are 0.26, 0.50, 0.62, 0.79, and 0.87 at a depth of 1 000 m and with bubble diameters of 0.005 m, 0.0075 m, 0.01 m, 0.02 m, and 0.03 m, respectively. The gas volume fraction curve has an inflection at points A and B, which are the lower and upper boundaries of the hydrate formation region, respectively. The gas volume fraction decreases when the gas enters the hydrate formation region because of the hydrate formation. Moreover, point A is near the mud line. The cross-sectional area of the riser section above the mud line is larger than the well annulus area. The sudden increase of the cross-sectional area reduces the gas volume fraction. The gas volume fraction increases when the gas leaves the hydrate formation region because of the hydrate disassociation. The hydrate formation would consume large amounts of gas and reduce the annular gas volume fraction, giving the deepwater kicks a “hidden” nature.

    Fig.3(a) Curve of pit gain versus kick time with different gas bubble diameters

    Fig.3(b) Curve of pit gain versus kick time with different gas bubble diameters

    (3) The influence of the hydrate phase transition on the pit gain during the gas kick period. Figures 3(a) and 3(b) show the pit gain versus the kick time at gas production rates of 0.15 Nm3/(MPa·s) and1.0 Nm3/(MPa·s), respectively. Each pit gain curve can be divided into three parts by points A and Bi(i=0,1,…5). Prior to point A, the invaded gas has not reached the hydrate formation region in the annulus, and thus, the pit gain curves with different bubble diameters are overlapped. When the gas enters the hydrate formation region (between A and Bi), the pit gain curves begin deviate from each other, and the deviation becomes larger over time because the hydrate formation consumes different amounts of gas. After point Bi, the front of the multiphase flow reaches the surface and some gase escapes from the annulus, which leads to the slow pit gain growth. The influence of the hydrate phase transition on the pit gain increases with the decreasing bubble diameter. The hydrate phase transition reduces the pit gain and delays the detection of the well ki3ck. As observed in Fig.3(a), the pit gain reaches 1 m after 98 min when the hydrate phase transition is not considered. On3the other hand, the pit gain takes 124 min to reach 1 m when the average bubble diameter is 0.005 m and the hydrate phase transition is considered. Therefore, the detection of the well kick is delayed by 26 min. The delays in the well kick detection greatly increase the well control risks. Moreover, the hydrate is more likely to be present in an annulus at a low gas production rate than at the high gas production rate when the gas kick is detected by using the pit gain monitoring equipment with certain preci3sions. For example, when the pit gain reaches 1 m, the gas would reach the hydrate formation region as shown in Fig.3(a) but not reach the hydrate formation region as shown in Fig.3(b).

    Fig.4 Curve of bottom hole pressure versus kick time with different gas bubble diameters

    (4) The influence of the hydrate phase transition on the bottom hole and shut-in casing pressures during the gas kick period. Figure 4 shows the curve of the bottom hole pressure versus the kick time. Figure 5 shows the curve of the shut-in casing pressure versus the kick time for different gas bubble diameters at the gas production rate of 1.0 Nm3/(MPa·s). As observed in these figures, the values of the bottom hole pressure and the shut-in casing pressure are dependent on the hydrate phase transition. The differences between the simulated results with and without the consideration of the hydrate phase transition increase with the decrease of the bubble diameter. The slopes of the curves decrease after point A (approximately 53 min) because the gas front enters the riser after 53 min and the section area of the riser is much larger than that of the annulus below the seafloor. Some wrong decision can be made when the degree of the gas kick or the type of invaded fluid is judged based on the pressure difference between the shut-in casing pressure and the shutin stand pressure because the hydrate phase transition reduces the shut-in pressure.

    Fig.5 Curve of shut-in casing pressure versus kick time with different gas diameters

    Fig.6 Hydrate formation regions with different methanol concentrations (Cm)

    2.2Result analysis of the deepwater gas kick simulations with different hydrate inhibitor concentrations

    Methanol is a common hydrate inhibitor. Figure 6 shows the hydrate formation region in the annulus with different methanol weight concentrations (Cm). The hydrate phase curves gradually move downward and the hydrate formation regions become smaller when the methanol concentration increases. Moreover, no hydrate is formed in the annulus when the methanol inhibitor concentration reaches 8% weight concentration. The hydrate formation regions and the hydrate overcooled temperatures decrease with the addition of inhibitors, consequently the hydrate formation ratesare decreased. Figure 7 shows the formation rate of the single bubble (diameter of 0.01 m) with different methanol concentrations (Cm) in the hydrate formation region. The hydrate formation rate decreases with the increase of the methanol concentration. Moreover, the maximum value of the hydrate formation rate with 5% methanol weight concentration is 1/7 of that without inhibitors being added.

    Fig.7 The single bubble’s formation rate with different methanol concentrations (Cm), when the bubble diameter is 0.01 m

    Fig.8 Gas fraction profiles with different methanol concentrations at 95 min during gas kick

    Fig.9 Pit gain versus kick time with different methanol concentrations

    The influence of the hydrate phase transition on the gas kick hydraulics, such as the gas volume fraction, the pit gain, is studied when the gas invasion rate is 0.10 Nm3/(MPa·s) and the average bubble diameter is 0.01 m. Figures 8 and 9 show the simulated results. The influence of the hydrate phase on the aforementioned parameters decreases with the increase of the methanol concentration. However, the influence of the hydrate phase transition on the flow parameters disappears when the methanol inhibitor weight concentration reaches 8%. As shown in Fig.8, the maximum overcooled temperature for the hydrate formation, the temperature difference between the annulus temperature curve and the hydrate phase equilibrium state curve (Cm=0) at depth of 1 057 m, is 6.2oC. Using the method proposed by Yamanlar et al.[19], the methanol weight concentration can be calculated as 7.6%, which is very close to the value 8% found in this study.

    The effect of the inhibitor concentration on the gas kick hydraulics should be considered in the deepwater gas kick simulations to better understand the deepwater gas kick process and provide guidances for future well control implementations.

    3. Conclusions

    A simulation model for deepwater gas kicks that accounts for the hydrate phase transition is established. The model is based on the improved annular multiphase flow governing equations and the hydrate phase transition calculation equations. The improved annular multiphase flow governing equations are derived by coupling the hydrate phase transition via adding the hydrate mass conservation equation, the hydrate term in the momentum conservation equation and the enthalpy term for the gas hydrate phase transition in the heat conservation equations. The model can be used quantitatively determine the relationship between the gas kick hydraulic parameters, such as the pit gain, the gas volume fraction, the bottom hole pressure and the shut-in casing pressure, versus the kick time and the impact of some hydrate inhibitor added into the drilling fluid.

    In the gas kick process, the gas bubble diameter can significantly influence the hydrate phase transition effect. Both the gas volume fraction in the annulus and the pit gain would be reduced because of the hydrate phase transition, which gives the well kicks in deepwater a “hidden” nature and increases the difficulty of the early detection of gas kicks. The hydrate phase curves gradually move downwards, and the hydrate formation region decreases when the hydrate inhibitor concentration increases. The hydrate overcooled temperature decreases and the hydrate formation rate and the influence of the hydrate phase transition on gas kick hydraulics are decreased. The influence of the hydrate phase transition on the kick parameters disappears when the methanol inhibitor weight concentration reaches 8% in this study.

    References

    [1] FOSSIL B., SANGESLAND S.. Managed pressure drilling for subsea applications: Well control challenges in deep waters[C]. SPE/IADC Underbalanced Technology Conference and Exhibition, 91633. Houston, Texas, USA, 2004.

    [2] NUNES J. O. L., BANNWART A. C. and RIBEIRO P. R. Mathematical model of a gas kick in deep water scenario[C]. IADC/SPE Asia Pacific Drilling Technology. Jakarta, Indonesia, 2002.

    [3] AVELAR C., RIBEIRO P. R. and SEPEHMOORI K. Deep-water gas kick simulation[J]. Journal of Petroleum Science and Engineering, 2009, 67(1-2): 13-22.

    [4] JOHNNY P., KNUT S. B. and KNUT L. Computing the danger of hydrate formation using a modified dynamic kick simulator[C]. SPE/IADC Drilling Conference. Amsterdam, The Netherlands, 2001.

    [5] WANG Z., SUN B. and CHENG H. Prediction of gas hydrate formation region in the well-bore of deepwater drilling[J]. Petroleum Exploration Development, 2008, 35(6): 731-735.

    [6] YAPA P. D., ZHENG L. and CHEN F. A model for deepwater oil/gas blowouts[J]. Marine Pollution Bulletin, 2001, 43(7): 234-241.

    [7] WANG Z., SUN B. Multiphase flow behavior in annulus with solid gas hydrate considering nature gas hydrate phase transition[J]. Petroleum Science, 2009, 6(1): 57-63.

    [8] JIANG G., WU Q. and ZHAN J. Experimental studies of the formation and dissociation of methane hydrate in loess[J]. Journal of Natural Gas Chemistry, 2010, 19(3): 217-223.

    [9] FAN S., ZHANG Y. and TIAN G. et al. Natural gas hydrate dissociation by presence of ethylene glycol[J]. Energy Fuels, 2006, 20(1): 324-326.

    [10] HAGOORT J. Prediction of wellbore temperatures in gas production wells[J]. Journal of Petroleum Science and Engineering, 2005, 49(1-2): 22-36.

    [11] HAGOORT J. Ramey’s wellbore heat transmission revisited. SPE Journal, 2004, 9(4): 465-474.

    [12] ZHANG H. Q., WANG Q. and SARICA C. et al. Unified model for gas-liquid pipe flow via slug dynamics-Part 1: Model development[J]. Journal of Energy Resources Technology, 2003, 125(4): 266-273.

    [13] JAVANMARDI J., MOSHFEGHIAN M. A new approach for prediction of gas hydrate formation conditions in aqueous electrolyte solutions[J]. Fluid Phase Equilibria, 2000, 168(2): 135-148.

    [14] ENGLEZOS P., KALOGERAKIS N. and DHOLABHAI P. D. et al. Kinetics of formation of methane and ethane gas hydrates[J]. Chemical Engineering Science, 1987, 42(11): 2647-2658.

    [15] HINZE J. O. Fundamentals of the hydrodynamic mechanism of splitting in dispersion processes[J]. AIChE Journal, 1955, 1(3): 289-295.

    [16] WINTERTON R. H. S. A simple method of predicting bubble size in bubble columns[J]. Chemical Engineering and Processing: Process Intensification, 1994, 33(1): 1-5.

    [17] KIM H. C., BISHNOI P. R. and HEIDEMANN R. A. et al. Kinetics of methane hydrate decomposition[J]. Chemical Engineering Science, 1987, 42(7): 1645-1653.

    [18] ZHENG L., YAPA P. D. Buoyant velocity of spherical and non-spherical bubbles/droplets[J]. Journal of Hydraic Engineering ASCE, 2000, 126(11): 825-854.

    [19] YAMANLAR S., POETTMANN F. H. and SLOAN E. D. Control hydrate formation[J]. Hydrocarbon Processing, 1991, 70(9): 155-161.

    10.1016/S1001-6058(14)60011-1

    * Project supported by the National Natural Science Foundation of China (Grant No. 51104172), the National Science and Technology Major Project of China (Grant No. 2011ZX05026-001) and the National High Technology Research and Development Program of China (863 Program, Grant No. 2013AA09A215).

    Biography: WANG Zhi-yuan (1981-), Ph. D.,

    Associate Professor

    猜你喜歡
    志遠(yuǎn)
    Corrigendum to“Atomic-scale electromagnetic theory bridging optics in microscopic world and macroscopic world”
    Atomic-scale electromagnetic theory bridging optics in microscopic world and macroscopic world
    禹志遠(yuǎn)作品
    大眾文藝(2022年24期)2023-01-09 09:27:16
    Topological photonic states in gyromagnetic photonic crystals:Physics,properties,and applications
    Quantum mechanical solution to spectral lineshape in strongly-coupled atom-nanocavity system
    呼志遠(yuǎn)美術(shù)作品
    我最喜愛(ài)的玩具①
    Atom interferometers with weak-measurement path detectors and their quantum mechanical analysis?
    香噴噴的年喲
    Functional Equivalence Theory and Its Limitations in Translation
    學(xué)周刊(2015年1期)2015-07-09 22:04:00
    成人特级av手机在线观看| 亚洲av男天堂| 亚洲乱码一区二区免费版| 精品久久久久久久久亚洲| 久久精品久久久久久噜噜老黄| 亚洲欧洲国产日韩| 国产 亚洲一区二区三区 | 久久精品国产自在天天线| 日韩欧美国产在线观看| 亚洲美女搞黄在线观看| 国产日韩欧美在线精品| 人人妻人人看人人澡| 18+在线观看网站| 国产成人精品福利久久| 插逼视频在线观看| 亚洲真实伦在线观看| 久久精品久久久久久噜噜老黄| 男人舔女人下体高潮全视频| 爱豆传媒免费全集在线观看| 精品一区在线观看国产| 在线观看免费高清a一片| 色播亚洲综合网| 99久国产av精品国产电影| 国产精品一二三区在线看| videos熟女内射| 极品少妇高潮喷水抽搐| 久久久精品免费免费高清| 欧美人与善性xxx| av天堂中文字幕网| 亚洲精品aⅴ在线观看| 91久久精品电影网| 国产 一区精品| 国产男人的电影天堂91| 日本av手机在线免费观看| 最近手机中文字幕大全| 最近的中文字幕免费完整| 99热这里只有是精品在线观看| 中文字幕制服av| 大又大粗又爽又黄少妇毛片口| 国产伦一二天堂av在线观看| av福利片在线观看| 欧美精品一区二区大全| 国产av在哪里看| 久久久久久久久大av| 国产精品爽爽va在线观看网站| 性插视频无遮挡在线免费观看| 国产精品久久久久久精品电影小说 | 97精品久久久久久久久久精品| 99热6这里只有精品| 最近最新中文字幕大全电影3| 只有这里有精品99| 精品人妻偷拍中文字幕| 婷婷色麻豆天堂久久| 亚洲国产精品专区欧美| 舔av片在线| 最后的刺客免费高清国语| 日韩一区二区三区影片| 乱码一卡2卡4卡精品| 色5月婷婷丁香| 日本-黄色视频高清免费观看| 国产黄色免费在线视频| 草草在线视频免费看| 日本猛色少妇xxxxx猛交久久| 午夜福利视频精品| 51国产日韩欧美| 直男gayav资源| 免费av观看视频| 亚洲最大成人中文| 永久免费av网站大全| 青春草视频在线免费观看| 色综合色国产| 精品人妻偷拍中文字幕| 97人妻精品一区二区三区麻豆| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产熟女欧美一区二区| 国产黄频视频在线观看| 我的女老师完整版在线观看| 男女国产视频网站| 超碰av人人做人人爽久久| 禁无遮挡网站| 国产精品伦人一区二区| 日韩欧美国产在线观看| 日韩在线高清观看一区二区三区| 免费观看在线日韩| 日本午夜av视频| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 两个人视频免费观看高清| 高清欧美精品videossex| 如何舔出高潮| 亚洲三级黄色毛片| 国产精品人妻久久久久久| 亚洲色图av天堂| 午夜日本视频在线| 又大又黄又爽视频免费| 床上黄色一级片| 一个人观看的视频www高清免费观看| 久99久视频精品免费| 国产精品久久久久久久久免| 久久精品熟女亚洲av麻豆精品 | 国产伦精品一区二区三区视频9| 舔av片在线| 国产精品久久视频播放| 秋霞在线观看毛片| 精品久久久久久久久久久久久| 欧美+日韩+精品| 国产有黄有色有爽视频| 亚洲精品成人久久久久久| 精品久久久久久电影网| 国产精品精品国产色婷婷| 一本久久精品| 亚洲成人一二三区av| 免费观看精品视频网站| 最近中文字幕2019免费版| 久久久午夜欧美精品| 一级毛片 在线播放| 美女xxoo啪啪120秒动态图| 麻豆成人午夜福利视频| 干丝袜人妻中文字幕| 成人美女网站在线观看视频| 国产成人午夜福利电影在线观看| 一个人看的www免费观看视频| 午夜福利高清视频| 久久这里只有精品中国| 久久久精品欧美日韩精品| 高清欧美精品videossex| 亚洲精品中文字幕在线视频 | 在线观看人妻少妇| 熟妇人妻久久中文字幕3abv| 成人毛片60女人毛片免费| 亚洲怡红院男人天堂| 欧美一区二区亚洲| 毛片女人毛片| av线在线观看网站| 国产在线男女| 99视频精品全部免费 在线| 国产 一区精品| av又黄又爽大尺度在线免费看| 最近视频中文字幕2019在线8| 少妇猛男粗大的猛烈进出视频 | 成人无遮挡网站| 高清欧美精品videossex| 亚洲成人一二三区av| 欧美一级a爱片免费观看看| 不卡视频在线观看欧美| 亚洲欧美一区二区三区黑人 | 精品少妇黑人巨大在线播放| 精品熟女少妇av免费看| 中文欧美无线码| 国产精品一区二区三区四区久久| 综合色丁香网| 2021少妇久久久久久久久久久| 18禁动态无遮挡网站| 成人综合一区亚洲| 亚洲av福利一区| 中文字幕av成人在线电影| 99久久中文字幕三级久久日本| 中国国产av一级| 久久午夜福利片| 一级av片app| a级一级毛片免费在线观看| 美女脱内裤让男人舔精品视频| 亚洲自拍偷在线| 777米奇影视久久| 午夜福利视频精品| 免费看日本二区| 亚洲欧美成人精品一区二区| 成年版毛片免费区| 国产真实伦视频高清在线观看| 国产视频首页在线观看| 日本一二三区视频观看| 久久99热6这里只有精品| 国产成人a区在线观看| 久久久久免费精品人妻一区二区| 日韩欧美一区视频在线观看 | 全区人妻精品视频| 国产 一区精品| 国产亚洲最大av| 69av精品久久久久久| 好男人在线观看高清免费视频| 久久精品久久久久久久性| 日韩电影二区| 三级经典国产精品| 91av网一区二区| 一边亲一边摸免费视频| 免费看美女性在线毛片视频| 免费观看的影片在线观看| 韩国高清视频一区二区三区| 免费少妇av软件| 国产精品一区二区在线观看99 | 夫妻午夜视频| 国产精品嫩草影院av在线观看| 别揉我奶头 嗯啊视频| 特级一级黄色大片| 一区二区三区高清视频在线| 少妇丰满av| 男女边吃奶边做爰视频| 最近手机中文字幕大全| 国产精品蜜桃在线观看| 秋霞伦理黄片| 熟女电影av网| 国产亚洲5aaaaa淫片| 有码 亚洲区| 精品欧美国产一区二区三| 熟妇人妻不卡中文字幕| 国产伦一二天堂av在线观看| 纵有疾风起免费观看全集完整版 | 亚洲高清免费不卡视频| 韩国高清视频一区二区三区| av专区在线播放| 一区二区三区高清视频在线| 日韩欧美三级三区| 亚洲国产最新在线播放| 欧美3d第一页| 亚洲人成网站在线观看播放| 日本免费在线观看一区| 99热网站在线观看| av在线播放精品| 日韩一区二区三区影片| a级一级毛片免费在线观看| 大话2 男鬼变身卡| 爱豆传媒免费全集在线观看| 日韩欧美精品免费久久| 国产在线一区二区三区精| 伊人久久精品亚洲午夜| 亚洲自拍偷在线| 欧美性猛交╳xxx乱大交人| 免费在线观看成人毛片| 99久久精品一区二区三区| 国产高潮美女av| 亚洲国产欧美人成| 水蜜桃什么品种好| 亚洲av成人精品一区久久| 精品熟女少妇av免费看| 亚洲av男天堂| 亚洲精品国产av成人精品| 丝瓜视频免费看黄片| 白带黄色成豆腐渣| 草草在线视频免费看| 国产免费视频播放在线视频 | 男人舔女人下体高潮全视频| 校园人妻丝袜中文字幕| 超碰97精品在线观看| 国产伦一二天堂av在线观看| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 在线免费观看不下载黄p国产| 亚洲人成网站在线观看播放| 日韩一本色道免费dvd| 少妇的逼好多水| 国产美女午夜福利| 精品一区二区三卡| 成人毛片60女人毛片免费| 麻豆国产97在线/欧美| 精品一区二区三区视频在线| 国产探花在线观看一区二区| 1000部很黄的大片| 精品久久久久久久久av| 能在线免费观看的黄片| 成年女人在线观看亚洲视频 | 国产av在哪里看| 日本爱情动作片www.在线观看| 中文字幕亚洲精品专区| 国产精品.久久久| www.色视频.com| 欧美成人一区二区免费高清观看| 男女那种视频在线观看| 国产女主播在线喷水免费视频网站 | 国产精品嫩草影院av在线观看| 欧美潮喷喷水| 日产精品乱码卡一卡2卡三| 日本av手机在线免费观看| 国产69精品久久久久777片| 亚洲国产精品国产精品| 日本免费在线观看一区| 777米奇影视久久| 久久久久精品性色| 久久亚洲国产成人精品v| 91av网一区二区| 欧美zozozo另类| 久久精品国产鲁丝片午夜精品| 老司机影院毛片| 国产精品人妻久久久影院| 熟女人妻精品中文字幕| 欧美高清成人免费视频www| 91aial.com中文字幕在线观看| 男女视频在线观看网站免费| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 青春草亚洲视频在线观看| 女人被狂操c到高潮| 白带黄色成豆腐渣| 国产男人的电影天堂91| 女人被狂操c到高潮| 免费观看av网站的网址| 五月玫瑰六月丁香| 亚洲高清免费不卡视频| 99热这里只有精品一区| 你懂的网址亚洲精品在线观看| 欧美精品国产亚洲| 麻豆精品久久久久久蜜桃| 看免费成人av毛片| 精品一区二区免费观看| 好男人在线观看高清免费视频| 国产伦理片在线播放av一区| 看非洲黑人一级黄片| 欧美日韩国产mv在线观看视频 | 毛片一级片免费看久久久久| 亚洲三级黄色毛片| 亚洲成人精品中文字幕电影| 一区二区三区乱码不卡18| 在线播放无遮挡| 亚洲一级一片aⅴ在线观看| 99热6这里只有精品| 九草在线视频观看| 久久久国产一区二区| 日本黄色片子视频| 欧美不卡视频在线免费观看| 精品午夜福利在线看| 亚洲欧美一区二区三区国产| 久久久久久久久久久免费av| 夫妻午夜视频| 亚洲精品中文字幕在线视频 | 美女国产视频在线观看| 亚洲va在线va天堂va国产| 在线播放无遮挡| 国产不卡一卡二| 简卡轻食公司| 欧美不卡视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 在线免费十八禁| 亚洲,欧美,日韩| 日韩不卡一区二区三区视频在线| 精品人妻熟女av久视频| 精品国产一区二区三区久久久樱花 | 97人妻精品一区二区三区麻豆| 性色avwww在线观看| 成人午夜精彩视频在线观看| 草草在线视频免费看| 国产美女午夜福利| 亚洲色图av天堂| 欧美另类一区| 成人二区视频| 免费观看在线日韩| 97在线视频观看| 尤物成人国产欧美一区二区三区| 看黄色毛片网站| 亚洲18禁久久av| 日韩一区二区三区影片| 一个人看的www免费观看视频| 国产精品国产三级国产av玫瑰| 久99久视频精品免费| 亚洲最大成人中文| 少妇的逼好多水| 久久精品国产亚洲网站| 你懂的网址亚洲精品在线观看| 免费大片18禁| 色尼玛亚洲综合影院| 亚洲精品乱码久久久久久按摩| 午夜福利在线观看吧| 亚洲精华国产精华液的使用体验| 欧美日韩在线观看h| 97热精品久久久久久| 中文字幕制服av| 日韩人妻高清精品专区| 免费观看性生交大片5| 国产日韩欧美在线精品| 热99在线观看视频| 超碰97精品在线观看| 在线观看av片永久免费下载| 国产视频首页在线观看| 亚洲精品色激情综合| a级毛色黄片| 在线天堂最新版资源| 国产精品爽爽va在线观看网站| 男女那种视频在线观看| 自拍偷自拍亚洲精品老妇| 高清欧美精品videossex| 不卡视频在线观看欧美| 麻豆国产97在线/欧美| av黄色大香蕉| 久久99蜜桃精品久久| 在线 av 中文字幕| 国产精品99久久久久久久久| 久久久久久伊人网av| 国产单亲对白刺激| 国产成人精品久久久久久| 免费观看性生交大片5| 成人一区二区视频在线观看| 黄片wwwwww| 人妻系列 视频| 国产精品久久久久久久电影| 超碰av人人做人人爽久久| 日韩电影二区| 精品一区二区免费观看| 建设人人有责人人尽责人人享有的 | 大陆偷拍与自拍| 欧美一区二区亚洲| 观看免费一级毛片| 久久99精品国语久久久| 亚洲四区av| 国产毛片a区久久久久| 国产成人福利小说| 日本三级黄在线观看| 日韩精品青青久久久久久| 最后的刺客免费高清国语| 内射极品少妇av片p| 麻豆国产97在线/欧美| 国产精品蜜桃在线观看| 中文字幕av成人在线电影| 国模一区二区三区四区视频| 亚洲欧美成人精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 美女高潮的动态| 久久久久网色| 精品国内亚洲2022精品成人| 国产一区二区在线观看日韩| 亚洲欧美中文字幕日韩二区| 色综合色国产| 国产成年人精品一区二区| 国产精品蜜桃在线观看| 身体一侧抽搐| 国产亚洲午夜精品一区二区久久 | 久久久色成人| 免费大片黄手机在线观看| 精品久久国产蜜桃| 欧美丝袜亚洲另类| 国产午夜福利久久久久久| 日本爱情动作片www.在线观看| 精品亚洲乱码少妇综合久久| 中文字幕av成人在线电影| 日本色播在线视频| 大片免费播放器 马上看| 亚洲性久久影院| 国产成年人精品一区二区| 亚洲精品456在线播放app| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| kizo精华| 97超视频在线观看视频| 欧美 日韩 精品 国产| 久久久久久九九精品二区国产| 久久韩国三级中文字幕| 免费av观看视频| 国产一区亚洲一区在线观看| 中文字幕制服av| 久久99蜜桃精品久久| 中文乱码字字幕精品一区二区三区 | 精品国产露脸久久av麻豆 | 亚洲高清免费不卡视频| 国产在视频线在精品| 午夜福利网站1000一区二区三区| av在线播放精品| 亚洲精华国产精华液的使用体验| 一级毛片 在线播放| 晚上一个人看的免费电影| 国产成人a∨麻豆精品| 亚洲在线观看片| 久久久久精品久久久久真实原创| 久久久久性生活片| 国产高潮美女av| 国语对白做爰xxxⅹ性视频网站| 国产伦精品一区二区三区四那| 精品国产三级普通话版| 成年免费大片在线观看| 国产女主播在线喷水免费视频网站 | 麻豆乱淫一区二区| 精品一区二区三区人妻视频| 好男人视频免费观看在线| 日韩,欧美,国产一区二区三区| or卡值多少钱| 午夜免费观看性视频| 99热网站在线观看| 成人亚洲精品av一区二区| 亚洲精品视频女| 成人av在线播放网站| 你懂的网址亚洲精品在线观看| 一个人观看的视频www高清免费观看| av在线亚洲专区| a级毛片免费高清观看在线播放| 国产免费又黄又爽又色| 久久久成人免费电影| 特大巨黑吊av在线直播| 汤姆久久久久久久影院中文字幕 | 在线观看人妻少妇| 男人舔女人下体高潮全视频| 成人鲁丝片一二三区免费| 男女边摸边吃奶| 精品欧美国产一区二区三| 肉色欧美久久久久久久蜜桃 | 六月丁香七月| 波多野结衣巨乳人妻| 亚洲精品日韩av片在线观看| 日日摸夜夜添夜夜爱| 欧美成人午夜免费资源| 看十八女毛片水多多多| 99久久精品热视频| 国产在视频线精品| 国产淫语在线视频| 只有这里有精品99| 啦啦啦中文免费视频观看日本| 欧美成人精品欧美一级黄| 男人舔奶头视频| 99热6这里只有精品| 男女边摸边吃奶| av女优亚洲男人天堂| 高清午夜精品一区二区三区| 欧美成人a在线观看| 伦理电影大哥的女人| 少妇裸体淫交视频免费看高清| 97超碰精品成人国产| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看| 色哟哟·www| 久久这里有精品视频免费| 日韩大片免费观看网站| 伊人久久国产一区二区| 99久久精品热视频| 尤物成人国产欧美一区二区三区| 久久99热这里只有精品18| 亚洲不卡免费看| 男女下面进入的视频免费午夜| 欧美成人精品欧美一级黄| 日韩中字成人| 欧美丝袜亚洲另类| 国产精品一区二区性色av| 好男人在线观看高清免费视频| 全区人妻精品视频| 午夜福利视频1000在线观看| 国产有黄有色有爽视频| 午夜福利视频精品| 久久久久久久久久人人人人人人| 久久韩国三级中文字幕| 网址你懂的国产日韩在线| 91久久精品国产一区二区成人| 免费看光身美女| 国精品久久久久久国模美| 黄色配什么色好看| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美日韩东京热| 欧美日韩视频高清一区二区三区二| 人妻一区二区av| 最近中文字幕高清免费大全6| 国产伦一二天堂av在线观看| 在线观看美女被高潮喷水网站| 国产精品99久久久久久久久| 丰满乱子伦码专区| 欧美日本视频| 国产精品女同一区二区软件| 国产伦一二天堂av在线观看| 69人妻影院| 国产毛片a区久久久久| 国产精品久久视频播放| 性插视频无遮挡在线免费观看| 女的被弄到高潮叫床怎么办| 日韩欧美精品v在线| 国产精品久久久久久久久免| 国产一区二区三区av在线| 欧美丝袜亚洲另类| 国产精品福利在线免费观看| 欧美日韩亚洲高清精品| 在线观看一区二区三区| 99热这里只有精品一区| 精品久久久久久电影网| 高清在线视频一区二区三区| 我的女老师完整版在线观看| 国产精品麻豆人妻色哟哟久久 | 国产亚洲av片在线观看秒播厂 | 2021少妇久久久久久久久久久| av网站免费在线观看视频 | 91狼人影院| 亚洲精品日韩在线中文字幕| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 99热这里只有是精品在线观看| 亚洲成人精品中文字幕电影| 亚洲av不卡在线观看| 成人国产麻豆网| 伊人久久国产一区二区| 日本wwww免费看| 日日啪夜夜爽| 国产男人的电影天堂91| 小蜜桃在线观看免费完整版高清| 国产视频首页在线观看| 亚洲在久久综合| 色播亚洲综合网| 肉色欧美久久久久久久蜜桃 | 欧美丝袜亚洲另类| 亚洲成人中文字幕在线播放| 伊人久久国产一区二区| 日韩大片免费观看网站| 免费人成在线观看视频色| 色综合站精品国产| 国产成人一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 国产精品1区2区在线观看.| 99久久精品国产国产毛片| 日韩一本色道免费dvd| 国产69精品久久久久777片| eeuss影院久久| 极品少妇高潮喷水抽搐| 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 国产黄色小视频在线观看| 成人亚洲精品一区在线观看 | 欧美最新免费一区二区三区| 男人和女人高潮做爰伦理| 男的添女的下面高潮视频| a级毛色黄片| 久久久久精品久久久久真实原创| 国产大屁股一区二区在线视频| 麻豆国产97在线/欧美| 99re6热这里在线精品视频|