生产油井井下温度场数值模拟分析

An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature field model in cylindrical coordinates is developed. In the model, we considered general heat conduction as well as the heat convection due to fluid flow from porous formation to the borehole. We also take into account the fluid velocity variation in the wellbore due to multiple pay zones. We present coupled boundary conditions at the interfaces between the wellbore and adjacent formation, the wellbore and pay zone, and the pay zone and adjacent formation. Finally, an alternating direction implicit difference method （ADI） is used to solve the temperature model for the downhole temperature distribution. The comparison of modeled temperature curve with actual temperature log indicates that simulation result is in general quite similar to the actual temperature log. We found that the total production rate, production time, porosity, thickness of pay zones, and geothermal gradient, all have effects on the downhole temperature distribution.     

Estimation of vertical water and heat fluxes in the semi-confined aquifers in tokyo metropolitan area,japan

Groundwater circulation is known to be one of the agents responsible for the redistribution of geothermal energy by acting as a source or sink in the course of its movement through porous media. Heat transport in groundwater systems is considered to be a coupled process and the theory based on this was used to analyse temperature profiles of 30 thermally stable observation wells in a deep, semi-confined aquifer system in the Tokyo Metropolitan area. Vertical water fluxes in the semi-confined aquifers and the associated upward heat fluxes were estimated from a heat flux equation that describes convection and conduction processes of heat transport in one dimension. The vertical downward water fluxes in Shitamachi lowland, Musashino and Tachikawa terraces were 0.69.26.91 × 10^?9, 1.46-70.92 × 10^?9 and 2.61.2204 × 10^?9 m/s, respectively. A vertical upward water flux of 1.80-33.60 × 10^?9 m/s was estimated in Shitamachi lowland. The water flux generally decreased with increasing depth for observation wells which intercepted more than one semi-confining layer. The estimated upward heat fluxes for Shitamachi lowland, Musashino and Tachikawa terraces were 0.32-1.12, 0.49-1.21 and 1.00-11.62 W/m², respectively. The heat flux was highest in Tachikawa terrace where a major fault, the Tachikawa fault, is located. Generally, the estimated heat flux was higher in the semi-confining layers than in the aquifers. Areas with heat sources and sinks as well as groundwater flow patterns in the semi-confined aquifers were revealed by heat flux and temperature distributions in the study area.     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.     

Non-equilibrium interphase heat and mass transfer during two-phase flow in porous media—Theoretical considerations and modeling

Mass and heat transfer occurring across phase-interfaces in multi-phase flow in porous media are mostly approximated using equilibrium relationships or empirical kinetic models. However, when the characteristic time of flow is smaller than that of mass or heat transfer, non-equilibrium situations may arise. Commonly, empirical approaches are used in such cases. There are only few works in the literature that use physically-based models for these transfer terms. In fact, one would expect physical approaches to modeling kinetic interphase mass and heat transfer to contain the interfacial area between the phases as a variable. Recently, a two-phase flow and solute transport model was developed that included interfacial area as a state variable [36]. In that model, interphase mass transfer was modeled as a kinetic process.     

Numerical analysis for cooling effect of open boundary ripped-rock embankment on Qinghai-Tibetan railway

At present, the Qinghai-Tibetan railway is being built, and it will pass across more than 550-km perma-frost regions. Therefore, the key to the stability of therailway embankment lies in solving the permafrost problem. Because global warming and existence of railway tend to degrade the permafrost in these re-gions[1], more difficulties and problems are induced in the construction and maintenance of railway. In the area where the mean annual air temperature is higher than a certain value, the …     

Topographic forcing of supercritical convection in a porous medium such as the oceanic crust

We have performed laboratory experiments using a Hele-Shaw cell to model a saturated, porous layer with various sinusoidal upper boundaries. Our intent was to determine the range of conditions over which boundary topography can control the pattern of thermal convection within a porous layer, and thereby take the first step toward understanding why heat flow seems correlated with hypsography in many areas of the ocean floor.These experiments indicate that above the critical Rayleigh number, topography does not control the convection pattern, except when the topographic wavelength is comparable to the depth of water penetration. Scaled to the depth of the layer, the convective wavenumbers are restricted to values between 2.5 and 4.8—a range which brackets π, the natural wavenumber for convection in a porous slab with planar, isothermal, impermeable boundaries. Topographies within this range control the circulation pattern perfectly, with downwelling under valleys and upwelling aligned with topographic highs. Other topographies do not force the pattern, although in some cases, the convection wavenumber may be a harmonic of the topographic wavenumber. Unforced circulation cells wander and vary in size, because they are not locked to the topography.For these experiments we employed eight different topographies with non-dimensional wavenumbers between 1.43 and 8.17, and we studied the flow at Rayleigh numbers between zero and five times the critical Rayleigh number. The amplitude of each topography tapered linearly (over a factor of three to six) from one end of the apparatus to the other, and the mean topographic amplitude was 0.05 times the depth of the layer. Under these conditions, amplitude has only a minor effect on the structural form and vigor of supercritical convection.Our results may apply to submarine geothermal systems, sealed by a thin layer of impermeable sediment draped over the basement topography. In this case, the convection wavelength—as measured perhaps by the spatial periodicity of conductive heat flow—may be a good measure of the depth to which water penetrates the crust. Where the circulation correlates with the bottom topography, it may be because the topographic wavelength is comparable to the depth to which water penetrates the porous crust.     

Physical,radiative, and dynamical processes within a nighttime marine stratus cloud

Observations taken by aircraft and conventional platforms are used to investigate dynamical, physical, and radiative processes within a marine stratus cloud during the Canadian Atlantic Storms Program (CASP) II field project which took place over the east coast of Canada. Stratus which formed over the ocean on February 6, 1992 during the nighttime, is studied to analyze cloud top and base processes. The cloud was supercooled during the study period. Fluctuations and fluxes are calculated along constant flight altitude legs approximately 100 km long in space. The scales of structures larger than 5 km are removed from the analysis using a running average technique. Droplet spectra obtained by a forward scattering spectrometer probe (FSSP) were used in a 1-D radiative transfer model to calculate infrared (IR) fluxes and radiative heating rates. A heat conservation equation was used to estimate vertical air velocity (w _a ) within the cloud. The results showed that, because of a warmer ocean surface, significant moisture and heat were transferred from the ocean surface to the boundary layer. The cloud base was at about 400 m height and the top was at about 1.4 km.w _a at the cloud base was estimated about 5 cm s^–1. Strong IR cooling rate at the cloud top was calculated to be 75°C day^–1 for a 100 m thick layer. Negative skewness inw _a , suggesting narrow downdrafts, was likely due to radiative cooling at the cloud top. The entrainment velocity was found to be about 1.5 cm s^–1 at cloud top. Mean moisture and heat fluxes within the cloud were estimated to be comparable to those from the ocean surface. Vertical air velocity at the cloud top due to radiative cooling was found to be about –40 cm s^–1.     

Data collection methodology for dynamic temperature model testing and corroboration

This article describes a data collection approach for determining the significance of individual heat fluxes within streams with an emphasis on testing (i.e. identification of possible missing heat fluxes), development, calibration and corroboration of a dynamic temperature model. The basis for developing this approach was a preliminary temperature modelling effort on the Virgin River in southwestern Utah during a low‐flow period that suggested important components of the energy balance might be missing in the original standard surface‐flux temperature model. Possible missing heat fluxes were identified as bed conduction, hyporheic exchange, dead zone warming and exchange and poor representation of the amount of solar radiation entering the water column. To identify and estimate the relative importance of the missing components, a comprehensive data collection effort was developed and implemented. In particular, a method for measuring shortwave radiation behaviour in the water column and an in situ method for separating out bed conduction and hyporheic influences were established. The resulting data and subsequent modelling effort indicate that hyporheic and dead zone heat fluxes are important, whereas solar radiation reflection at the water surface was found to be insignificant. Although bed conduction can be significant in certain rivers, it was found to have little effect on the overall heat budget for this section of the Virgin River. Copyright © 2009 John Wiley & Sons, Ltd.     

Transport of nonsorbing solutes in a streambed with periodic bedforms

Previous studies of hyporheic zone focused largely on the net mass transfer of solutes between stream and streambed. Solute transport within the bed has attracted less attention. In this study, we combined flume experiments and numerical simulations to examine solute transport processes in a streambed with periodic bedforms. Solute originating from the stream was subjected to advective transport driven by pore water circulation due to current–bedform interactions as well as hydrodynamic dispersion in the porous bed. The experimental and numerical results showed that advection played a dominant role at the early stage of solute transport, which took place in the hyporheic zone. Downward solute transfer to the deep ambient flow zone was controlled by transverse dispersion at the later stage when the elapsed time exceeded the advective transport characteristic time t_c (= L/u_c with L being the bedform length and u_c the characteristic pore water velocity). The advection-based pumping exchange model was found to predict reasonably well solute transfer between the overlying water and streambed at the early stage but its performance deteriorated at the later stage. With dispersion neglected, the pumping exchange model underestimated the long-term rate and total mass of solute transfer from the overlying water to the bed. Therefore both advective and dispersive transport components are essential for quantification of hyporheic exchange processes.     

The effect of stratification and compressibility on anelastic convection in a rotating plane layer

We study the effect of stratification and compressibility on the threshold of convection and the heat transfer by developed convection in the nonlinear regime in the presence of strong background rotation. We consider fluids both with constant thermal conductivity and constant thermal diffusivity. The fluid is confined between two horizontal planes with both boundaries being impermeable and stress-free. An asymptotic analysis is performed in the limits of weak compressibility of the medium and rapid rotation (τ^?1/12???|θ|???1, where τ is the Taylor number and θ is the dimensionless temperature jump across the fluid layer). We find that the properties of compressible convection differ significantly in the two cases considered. Analytically, the correction to the characteristic Rayleigh number resulting from small compressibility of the medium is positive in the case of constant thermal conductivity of the fluid and negative for constant thermal diffusivity. These results are compared with numerical solutions for arbitrary stratification. Furthermore, by generalizing the nonlinear theory of Julien and Knobloch [Fully nonlinear three-dimensional convection in a rapidly rotating layer. Phys. Fluids 1999, 11, 1469–1483] to include the effects of compressibility, we study the Nusselt number in both cases. In the weakly nonlinear regime we report an increase of efficiency of the heat transfer with the compressibility for fluids with constant thermal diffusivity, whereas if the conductivity is constant, the heat transfer by a compressible medium is more efficient than in the Boussinesq case only if the specific heat ratio γ is larger than two.     

Soret and Dufour effects on thermohaline convection in rotating fluids

Using linear and weakly nonlinear stability theory, the effects of Soret and Dufour parameters are investigated on thermohaline convection in a horizontal layer of rotating fluid, specifically the ocean. Thermohaline circulation is important in mixing processes and contributes to heat and mass transports and hence the earth’s climate. A general conception is that due to the smallness of the Soret and Dufour parameters their effect is negligible. However, it is shown here that the Soret parameter, salinity and rotation stabilise the system, whereas temperature destabilises it and the Dufour parameter has minimal effect on stationary convection. For oscillatory convection, the analysis is difficult as it shows that the Rayleigh number depends on six parameters, the Soret and Dufour parameters, the salinity Rayleigh number, the Lewis number, the Prandtl number, and the Taylor number. We demonstrate the interplay between these parameters and their effects on oscillatory convection in a graphical manner. Furthermore, we find that the Soret parameter enhances oscillatory convection whereas the Dufour parameter, salinity Rayleigh number, the Lewis number, and rotation delay instability. We believe that these results have not been elucidated in this way before for large-scale fluids. Furthermore, we investigate weakly nonlinear stability and the effect of cross diffusive terms on heat and mass transports. We show the existence of new solution bifurcations not previously identified in literature.     

Simultaneous crystallization and melting at both the roof and floor of crustal magma chambers: Experimental study using NH4Cl–H2O binary eutectic system

Thermal and compositional evolution of magmas after emplacement of basalt into continental crust has been investigated by means of fluid dynamic experiments using a cold solid mixture with eutectic composition and a hot liquid with higher salinity in the NH₄Cl–H₂O binary eutectic system. The experiments were designed to simulate cases where crystallization of a basalt magma is accompanied by melting at both the roof and floor of a crustal magma chamber. The results show that thermal and compositional convection occur simultaneously in the solution; the thermal convection is driven by cooling at the roof and the compositional convection is driven by melting and crystallization at the floor. The roof was rapidly melted by the convective heat flux, which resulted in formation of a separate eutectic melt layer (the upper liquid layer) with negligible mixing of the underlying liquid (the lower liquid layer). On the other hand, a mushy layer formed at the floor. The compositional convection at the floor carried a low heat flux, so that the heat transfer at the floor was basically explained by simple heat conduction. The thermal boundary layer in the lower liquid layer at the interface with the upper liquid layer became thicker with time and subsequently temperature decreased upward throughout the lower liquid layer. Compositional gradient with NH₄Cl content decreasing upward formed by compositional convection in the lower liquid layer. The formation of these gradients resulted in formation of double-diffusive convecting layers in the lower liquid layer. The upward heat transfer was suppressed when compared with the case where the liquid region is homogenized by vigorous convection.These experimental results imply that, when a basalt magma is emplaced in continental crust, floor melting does not always enhance the cooling of the magma, but it may even reduce the total heat loss from the magma to the crusts due to suppression of convection by formation of a stabilizing compositional gradient.     

The effects of small water bodies on the atmospheric heat and water budgets over the MacKenzie River Basin

The effects of small water bodies or lakes on the surface sensible and latent heat fluxes and the transport of heat and water vapour in the atmospheric boundary layer (ABL) over the Mackenzie River Basin (MRB) are studied from two cases, which occurred on 2 and 8 June 1999 during the warm season. The synoptic condition for the cases is representative of about 33% of the synoptic situation over the MRB. The two events are simulated using the Canadian mesoscale compressible community (MC2) model. A one‐way nesting grid approach is employed with the highest resolution of 100 m over a domain of 100 km². Experiments were conducted with (LAKE) and without (NOLAKE) the presence of small water bodies, whose size distribution is obtained through an inversion algorithm using information of their linear dimension determined from aircraft measurement of surface temperature during MAGS (the Mackenzie GEWEX (Global Energy and Water Cycle Experiment) Study) in 1999. The water bodies are assumed to be distributed randomly in space with a fractional area coverage of 10% over the MRB. The results show that, in the presence of lakes, the domain‐averaged surface sensible heat flux on 2 June 1999 (8 June 1999) decreases by 9·3% (6·6%). The surface latent heat flux is enhanced by 18·2% (81·5%). Low‐level temperature advection and the lake surface temperature affect the air–land/lake temperature contrast, which in turn controls the sensible heat flux. In the absence of lakes the surface wind speed impacts the latent heat flux, but in the presence of lakes the moisture availability and the atmospheric surface layer stability control the latent heat flux. The enhancement is smaller on 2 June 1999 as a result of a stable surface layer caused by the presence of colder lake temperatures. The domain‐averaged apparent heat source and moisture sink due to turbulent transports were also computed. The results show that, when lakes are present, heating and drying occur in the lowest 100 m from the surface. Above 100 m and within the ABL, there was apparent cooling. However, the apparent moistening profiles reveal that lakes tend to moisten the ABL through transfer of moisture from the lowest 50–100 m layer. Copyright © 2004 John Wiley & Sons, Ltd.     

A land surface model incorporated with soil freeze/thaw and its application in GAME/Tibet

Land surface process is of great importance in global climate change, moisture and heat exchange in the interface of the earth and atmosphere, human impacts on the environment and eco- system, etc. Soil freeze/thaw plays an important role in cold land surface processes. In this work the diurnal freeze/thaw effects on energy partition in the context of GAME/Tibet are studied. A sophisti- cated land surface model is developed, the particular aspect of which is its physical consideration of soil freeze/thaw and vapor flux. The simultaneous water and heat transfer soil sub-model not only reflects the water flow from unfrozen zone to frozen fringe in freezing/thawing soil, but also demon- strates the change of moisture and temperature field induced by vapor flux from high temperature zone to low temperature zone, which makes the model applicable for various circumstances. The modified Picard numerical method is employed to help with the water balance and convergence of the numerical scheme. Finally, the model is applied to analyze the diurnal energy and water cycle char- acteristics over the Tibetan Plateau using the Game/Tibet datasets observed in May and July of 1998. Heat and energy transfer simulation shows that: (i) There exists a negative feedback mechanism between soil freeze/thaw and soil temperature/ground heat flux; (ii) during freezing period all three heat fluxes do not vary apparently, in spite of the fact that the negative soil temperature is higher than that not considering soil freeze; (iii) during thawing period, ground heat flux increases, and sensible heat flux decreases, but latent heat flux does not change much; and (iv) during freezing period, soil temperature decreases, though ground heat flux increases.     

A preliminary study on the heat storage fluxes of a tropical seasonal rain forest in Xishuangbanna

In order to discuss the values and daily variation characteristics of heat storage fluxes in a tropical seasonal rain forest in Xishuangbanna, the sensible and latent heat storage flux within air column, canopy heat storage flux, energy storage by photosynthesis and ground heat storage above the soil heat flux plate, as well as the ratios of these heat storage fluxes to the net radiation in the cool-dry, hot-dry and rainy season were compared and analyzed based on the observation data of carbon fluxes, meteorological factors and biomass within this tropical seasonal rain forest from January 2003 to December 2004. The findings showed that heat storage terms ranged significantly in the daytime and weakly in the nighttime, and the absolute values of sensible and latent heat storage fluxes were obviously greater than other heat storage terms in all seasons. In addition, the absolute values of total heat storage fluxes reached the peak in the hot-dry season, then were higher in the rainy season, and reached the minimum in the cool-dry season. The ratios of heat storage fluxes to net radiation generally decreased with time in the daytime, moreover, the sensible and latent heat storage dominated a considerable fraction of net radiation, while other heat storage contents occupied a smaller fraction of the net radiation and the peak value was not above 3.5%. In the daytime, the ratios of the total heat storage to net radiation were greater and differences in these ratios were distinct among seasons before 12:00, and then they became lower and differences were small among seasons after 12:00. The energy closure was improved when the storage terms were considered in the energy balance, which indicated that heat storage terms should not been neglected. The energy closure of tropical seasonal rain forest was not very well due to effects of many factors. The results would help us to further understand energy transfer and mass exchange between tropical forest and atmosphere. Moreover, they would supply a research basis for studying energy closure at other places.     

Combined Effect of Climatic Variations and Groundwater Movement on Observed Heat Flow

In most practical situations, the upper part of a geological section consists of loose sediments, in which heat transfer cannot be described as a purely conductive process. To investigate such situations a one-dimensional numerical model of terrestrial temperature field formation under the combined influence of vertical groundwater filtration and ground surface temperature changes has been developed. The model allows one to consider the perturbation of heat flow interval values resulting from short- and long-period temperature waves propagating into permeable rocks under conditions of advective heat transfer, caused by vertical groundwater filtration. The results show that temperature profiles and interval heat flow values are sensitive to both the paleoclimatic history and the rate of groundwater filtration. The latter plays the prevailing role in the variations of geothermal field parameters, especially within the uppermost part of the loose sediments in unconfined aquifers. The problem was solved for a permeable layer, underlaid by an impermeable layer. This schematisation of water exchange is the typically accepted for hydrogeological analysis. Even at very low rates of filtration the intensity of this effect is enhanced substantially for long-period variations. In the extreme case (for periods of temperature variations of the order of 100,000 years) at typical rates of filtration within the permeable layer, an almost gradient-free zone can be formed down to depths of a few hundred metres. For the case of upward filtration, on the contrary, the influence of climatic variations on the terrestrial temperature field becomes substantially attenuated.     

Sensitivity of model parameterizations for simulated latent heat flux at the snow surface for complex mountain sites

The snowcover energy balance is typically dominated by net radiation and sensible and latent heat fluxes. Validation of the two latter components is rare and often difficult to undertake at complex mountain sites. Latent heat flux, the focus of this paper, is the primary coupling mechanism between the snow surface and the atmosphere. It accounts for the critical exchange of mass (sublimation or condensation), along with the associated snowcover energy loss or gain. Measured and modelled latent heat fluxes at a wind‐exposed and wind‐sheltered site were compared to evaluate variability in model parameters. A well‐tested and well‐validated snowcover energy balance model, Snobal, was selected for this comparison because of previously successful applications of the model at these sites and because of the adjustability of the parameters specific to latent heat transfer within the model. Simulated latent heat flux and snow water equivalent (SWE) were not sensitive to different formulations of the stability profile functions associated with heat transfer calculations. The model parameters of snow surface roughness length and active snow layer thickness were used to improve latent heat flux simulations while retaining accuracy in the simulation of the SWE at an exposed and sheltered study site. Optimal parameters for simulated latent heat flux and SWE were found at the exposed site with a shorter roughness length and thicker active layer, and at the sheltered site with a longer roughness length and thinner active layer. These findings were linked to physical characteristics of the study sites and will allow for adoption into other snow models that use similar parameters. Physical characteristics of wind exposure and cover could also be used to distribute critical parameters in a spatially distributed modelling domain and aid in parameter selection for application to other watersheds where detailed information is not available. Copyright © 2012 John Wiley & Sons, Ltd.     

TOGA-COARE IOP期间的海气通量观测结果

在TOGA－COAREIOP期间用涡度相关法对海气热通量进行了船载直接观测．对垂直风速、温度和湿度湍流脉动观测数据的谱分析显示它们在高频区基本满足"－2／3次方律"．对船体简谐震荡影响的讨论从理论上证明该影响在热通量计算中可被忽略．根据以此方法得到的通量求出了中性层结条件下感热和潜热的整体输送系数分别为2．25×10^－3和1．26×10^－3．对海气边界层特性的分析表明该海域的近海层主要呈中性或弱不稳定层结．海气通量的变化与背景环流形势密切相关,潜热通量主要受海面风场强度的影响,而感热通量变化除了风场的影响外,层结变化也是一个重要因素．用整体输送法计算TOGA－COAREIOP期间以及TOGA期间8个航次的通量结果而得到的Bowen比约为0．1,显示潜热通量是暖池大气的主要热源．