A model for downhole fluid and rock temperature prediction during circulation

Abstract In this paper, we present a model to investigate the evolution of fluid and rock temperature during fluid circulation in a wellbore. The analysis considers circulation of a fluid down a centralized drill or tubing string with the returned fluid travelling up the annular space between the tubing and the wellbore. Under such conditions, which typically occur during drilling, the cooler injected fluid is heated as it travels down the tubing and cools the wellbore rock as it returns up the annulus. Based on the established governing equations for heat transfer between fluid and rock, a semi-analytic method is developed by applying Laplace transformation and numerical inversion to find the results in time and space. The heat transfer coefficients between rock and fluid are dependent on flow behaviour and material properties, characterizing advective heat transfer under complex flow. A dimensional analysis is conducted to identify the controlling dimensionless parameters. The solutions are validated by comparisons with theoretical predictions of heat diffusion inside the rock and with measured downhole temperature variations. The results show that the injection rate plays an important role in the downhole temperature evolution. The surface outlet temperature of the fluid from the annulus typically reaches a pseudo-steady state in a relatively rapid manner. Additionally, the strong cooling resulting from injection of cold fluid with circulation back up the annulus may cause significant thermoelastic changes in rock stress near the wellbore, potentially leading to tensile hydraulic fracturing initiation.

[1]  Robert G. Jeffrey,et al.  A semi-analytic solution of a wellbore in a non-isothermal low-permeability porous medium under non-hydrostatic stresses , 2012 .

[2]  F. Incropera,et al.  Fundamentals of Heat Transfer , 1981 .

[3]  H. J. Ramey.  Wellbore Heat Transmission , 1962 .

[4]  Hassan Hassanzadeh,et al.  Comparison of different numerical Laplace inversion methods for engineering applications , 2007, Appl. Math. Comput..

[5]  C. S. Kabir,et al.  Determining Circulating Fluid Temperature in Drilling, Workover, and Well Control Operations , 1996 .

[6]  H. R. Parkison,et al.  Calculation of Formation Temperature Disturbances Caused by Mud Circulation , 1962 .

[7]  W. M. Cobb,et al.  Static Formation Temperature From Well Logs - An Empirical Method , 1975 .

[8]  Edgar Santoyo,et al.  Convective heat-transfer coefficients of non-Newtonian geothermaldrilling fluids , 2003 .

[9]  S. Churchill.  A comprehensive correlating equation for laminar, assisting, forced and free convection , 1977 .

[10]  Lap Mou Tam,et al.  Transitional Heat Transfer in Plain Horizontal Tubes , 2006 .

[11]  L. R. Raymond.  Temperature Distribution in a Circulating Drilling Fluid , 1969 .

[12]  G. Willhite.  Over-all Heat Transfer Coefficients in Steam And Hot Water Injection Wells , 1967 .

[13]  A. Ghassemi,et al.  A three-dimensional thermo-poroelastic model for fracture response to injection/extraction in enhanc , 2011 .

[14]  A. Ghassemi,et al.  Propagation of a System of Cracks Under Thermal Stress , 2011 .

[15]  J. M. Coulson,et al.  Heat Transfer , 2018, Finite Element Method for Solids and Structures.

[16]  F. C. Arnold.  Temperature Variation in a Circulating Wellbore Fluid , 1990 .

[17]  B. B. Williams,et al.  Prediction of Wellbore Temperatures During Mud Circulation and Cementing Operations , 1973 .

[18]  Curtis M. Oldenburg,et al.  Process Modeling of CO2 Injection into Natural Gas Reservoirs for Carbon Sequestration and Enhanced Gas Recovery , 2001 .

[19]  Edgar Santoyo,et al.  Estimation of Temperatures in Geothermal Wells During Circulation and Shut-in in the Presence of Lost Circulation , 1998 .

[20]  V. Gnielinski.  New equations for heat and mass transfer in turbulent pipe and channel flow , 1976 .

[21]  Ahmad Ghassemi,et al.  Poro-thermoelastic borehole stress analysis for determination of the in situ stress and rock strength , 2010 .

[22]  T. Hashida,et al.  Heat flow rate at a bore-face and temperature in the multi-layer media surrounding a borehole , 2003 .

[23]  E. C. Bullard,et al.  The Time Necessary for a Bore Hole to Attain Temperature Equilibrium , 1947 .

[24]  S. C. Swift,et al.  Calculation of Circulating Mud Temperatures , 1970 .

[25]  Zhongming Chen,et al.  Accurate Prediction Wellbore Transient Temperature Profile Under Multiple Temperature Gradients: Finite Difference Approach and Case History , 2003 .

[26]  Sergei Fomin,et al.  A borehole temperature during drilling in a fractured rock formation , 2005 .

[27]  R. M. Beirute.  A circulating and shut-in well-temperature-profile simulator , 1991 .

[28]  E. N. Sieder,et al.  Heat Transfer and Pressure Drop of Liquids in Tubes , 1936 .

[29]  E. J. Couch,et al.  Temperature Distribution in Circulating Mud Columns , 1973 .

[30]  Y. Shibuya,et al.  Thermoelastic Analysis of a Cracklike Reservoir in a Hot Dry Rock During Extraction of Geothermal Energy , 1983 .

[31]  E. Sparrow,et al.  Internal-flow Nusselt numbers for the low-Reynolds-number end of the laminar-to-turbulent transition regime , 2011 .

[32]  Gilberto Espinosa-Paredes,et al.  TEMPLOP/V.2: a computer program for estimation of fully transient temperatures in geothermal wells during circulation and shut-in , 2001 .

[33]  G. R. Wooley,et al.  Computing Downhole Temperatures in Circulation, Injection, and Production Wells , 1980 .

[34]  Paul B. Crawford,et al.  A Method for Calculating Circulating Temperatures , 1967 .