Investigations in isochemical contact metamorphism

Isochemical contact metamorphism was observed in 3 localities based on 184 chemical analyses of rock samples and the significance tested using simple statistical techniques. The intrusions included gabbro, granosyenite, and granitic type rocks being intruded into respective schists, clays and shales, and schists and gneisses. Conductive heat transfer appeared to be the most important heat transfer mechanism in isochemical contact metamorphism. During isochemical metamorphism both H₂O and CO₂ were rather mobile and impoverished in the metamorphic rocks.     

壁面局部接触裂隙岩体水流-传热的一种理论模型及计算分析

核废物地质处置、地热开发、石油开采等工程领域都可能涉及稀疏裂隙岩体中的水流-传热过程。现有的裂隙岩体水流-传热理论模型和计算方法基本上都是以平行光滑壁面裂隙模型为基础的,没有考虑裂隙的壁面局部接触对水流、水-岩热交换以及岩体传热的影响。针对粗糙壁面裂隙水流过程,阐述了基于Stokes方程的Reynolds润滑方程及Hele-Shaw裂隙模型,采用MATLAB软件中的PDE工具求解,并与Walsh的等效水力开度公式进行对比;分析壁面局部接触裂隙水流-传热与填充裂隙水流-传热的相似性,提出了瞬时局部热平衡假设的适用条件,并在裂隙局部接触体传热满足Biot数条件的前提下,计算分析裂隙局部接触体与水流之间的局部热平衡时间及其影响因素;在裂隙局部接触体与水流之间满足瞬时热平衡假设的前提下,利用填充裂隙水流-传热的解析解,计算了壁面局部接触裂隙水及两侧岩石的温度分布,并分析了裂隙局部接触面积率、裂隙开度、裂隙水平均流速对岩石温度和裂隙水温度的影响特征,结果表明：（1）在设定条件下,由于裂隙局部接触体与裂隙水流之间的热交换,裂隙水流对其两侧岩石温度的影响范围随接触面积率的增大而减小,裂隙两侧岩石对裂隙水流温度的影响程度随接触面积率的增大而增大;（2）裂隙开度和裂隙水流速对岩石温度和裂隙水温度的影响方式的影响是一致的,即由于裂隙水流量随裂隙开度和裂隙水流速的增大而增大,裂隙水流对其两侧岩石温度的影响范围随裂隙开度和裂隙水流速的增大而增大,裂隙两侧岩石对裂隙水流温度的影响程度随裂隙开度和裂隙水流速的增大而减小。     

Dependence of thermal conductivity of snow on microstructure

A geometrical model, including different geometrical shapes influencing thermal conductivity of snow is proposed. The geometrical model has been assumed to comprise of unit cells having solid (ice) inclusion as an aggregation of spherical, cylindrical or cubical shapes with vertical connection, arranged in a cubic packing. From the geometrical model and one-dimensional heat transfer theory, the effective thermal conductivity has been computed. For this purpose, coupled one-dimensional heat transfer equations have been solved for steady-state condition to account for conduction in ice, conduction in air and latent heat transfer due to water vapour sublimation through air. The model demonstrates the dependency of thermal conductivity on density, grain-spacing, grain contact ratio and temperature. Spherical inclusions give highest conductivity while cubical inclusion estimates lowest value for the same density. Thermal conductivity has been found increasing sharply near to the packing density for all three shapes. Empirical model results and results obtained from existing microstructure based models have also been compared with the present model.     

花岗岩粗糙单裂隙对流换热特性的数值模拟

在干热岩储层中开采地热能,往往需要对储层进行人工水力压裂以形成贯穿的换热通道。然而,热储中的对流换热对干热岩的采热率有重要影响,经过人工刺激的储层会形成几何形态各异的裂隙面,而裂隙粗糙程度的不同则会引起换热性能的显著差异。因此,选取4条Barton的经典岩石裂隙粗糙度曲线,在试验室条件下建立一个单裂隙对流换热模型。详细分析了花岗岩粗糙裂隙中热工质的换热特性。结果表明：局部对流换热系数沿着裂隙长度方向逐渐降低;节理粗糙系数JRC值越大,平均对流换热系数就越大,表明换热性能越好;局部对流换热系数的分布与JRC曲线的几何轮廓形态有很好的相关性,波峰波谷的变化趋势相一致;相对于温度而言,高流速对局部对流换热系数具有放大效应,流速越大,局部对流换热系数波动越大。     

A three-dimensional heat transfer and thermal cracking model considering the effect of cracks on heat transfer

A simple three-dimensional heat transfer model is developed to consider the hindering effect of cracks on heat transfer. The 3D heat transfer model can also be applied to numerical methods such as the combined finite-discrete element method (FDEM), discrete element method (DEM), discontinuous deformation analysis (DDA), the numerical manifold method (NMM), and the finite element method (FEM) to construct thermo-mechanical coupling models that allow these methods to solve thermal cracking problems and dynamically consider the hindering effect of cracks on heat transfer. In the 3D heat transfer model, the continuous-discontinuous medium is discretized into independent tetrahedral elements, and joint elements are inserted between adjacent tetrahedral elements. Heat transfer calculations for continuous-discontinuous media are converted to heat conduction in tetrahedral elements and the heat exchange between the adjacent tetrahedral elements through the joint element. If the joint element between adjacent tetrahedral elements breaks (ie, a crack generates), the heat exchange coefficient of the joint element is reduced to account for the hindering effect of cracks on heat conduction. Then the model and the FDEM are combined to build a thermo-mechanical coupling model to simulate thermal cracking. The thermally induced deformation, stress, and cracking are investigated by the thermo-mechanical coupling model, and the numerical results are compared with analytical solutions or experimental results. The 3D heat transfer model and thermo-mechanical model can provide a powerful tool for simulating heat transfer and thermal cracking in a continuous-discontinuous medium.     

基于分布式光纤传感技术的渗流监测理论研究

从阐述多孔介质传热过程出发,逐一分析了在非渗流及渗流情况下,光纤与多孔介质之间的传热过程,详细分析了光纤与水流对流的传热过程,推导了多孔介质中渗流水和光纤之间的换热系数计算公式。在一定假设条件下,推导了导热系数（渗流流速）、加热功率和温升的非渗流情况与渗流情况监测理论方程,为分布式光纤传感技术监测渗流的应用提供了坚实的理论基础。     

透壁通风管路堤的对流换热和蒸发散热

针对透壁通风管路堤中透壁通风管管壁与空气之间的对流换热和土体水分通过管壁小孔的蒸发散热机制,分析了开孔率、风速及含水率等因素对透壁通风管管壁对流换热和水分蒸发散热的影响,并具体给出了管壁对流换热系数和蒸发散热系数的计算公式。冬季路堤由于冻土层未冻水含量较小而使管壁小孔的水分蒸发散热较弱,路堤总的降温效果主要由管壁对流换热效应控制,而暖季通风管内空气与管壁的对流换热效应可使路堤土体增温,同时,由于通风管周围融土的未冻水含量较大,而使得通过管壁小孔的水分蒸发散热较强,可部分或全部抵消对流换热引起的增温效应,而有利于路堤的稳定。     

中深层地热井下同轴换热器长期换热性能研究

中深层地热井下同轴换热器具有取热量大、出口水温高的特点,近年来受到供热行业的高度关注。目前的研究主要是对单个供热周期内的中深层地热井下同轴换热器的换热性能进行分析,缺乏对其长期换热性能的相关研究。本文根据中深层地热井下同轴换热器供暖季供热、非供暖季停歇运行的特点,基于能量守恒方程建立换热器的数值换热模型,并采用有限体积法对模型进行离散,通过Matlab平台数值分析其在长期运行过程中换热性能的变化规律。结果表明:换热器的换热性能随着运行年份逐渐下降,且下降程度逐年减小,最终达到准稳态。其中,次年平均换热量的下降比例最大,且换热器埋深越大,换热量的下降比例越小。埋深为2 000、2 500、3 000 m换热器的次年平均换热量的下降比例依次为4.00%,3.78%和3.56%,第30年的平均换热量较第1年分别下降13.7%、13.1%、12.4%。岩体温度逐年下降,其受干扰半径逐年增加。在30 a运行期间,埋深2 000 m的换热器在每年供热季结束时的岩体温度受干扰半径从13 m增加至105 m。此外,换热器深度越大,其周围岩体温度受干扰半径越大。本研究结果阐明了中深层地热井下同轴换热器在长期取热过程中换热性能的变化规律,对换热器长期取热的设计具有指导意义。      

花岗岩单裂隙渗流传热特性试验

基于岩石裂隙渗流传热试验系统,以蒸馏水为换热工质,针对预制平滑裂隙与粗糙裂隙的花岗岩岩样,开展不同试验温度与水流条件下的单裂隙对流换热特性试验研究。试验结果表明：①同一温度水平下,对流换热系数的量值与流量呈正相关关系,但对流换热系数的增幅大于流量的增幅,随岩石基质温度的升高,相同流量水平下的对流换热系数呈增大趋势;②与平滑裂隙试验结果相比,粗糙裂隙面的对流换热强度有所提高,但增幅不大;③沿渗流路径,裂隙面局部对流换热强度的演化呈现多波峰的非线性特征,其中粗糙裂隙的非线性更加强烈。但二者均存在由进水口至出水口,换热强度逐渐减弱的特征。④努塞尔数与普朗特数1/3次方的商与雷诺数较好地满足幂指数关系,随着试验温度的升高,这种关系逐渐向线性转变。     

Discrete fracture model for simulating waterflooding processes under fracturing conditions

In our study, we develop a model for simulating fracturing processes in a poroelastic medium. The proposed approach combines the discrete fracture model enriched with contact plane mechanics. The model captures mechanical interactions of fractures and a deformable medium, fluid, and heat transfer in fractures and in a porous medium. Both effects of poroelasticity and thermoelasticity are accounted in our model. Mass and heat conservation equations are approximated by the finite volume method, and mechanical equilibrium equations are discretized by means of the Galerkin finite element approach. Two‐dimensional grid facets between 3‐dimensional finite elements are considered as possible fracture surfaces. Most of these facets are inactive from the beginning and are activated throughout the simulation. A fracture propagation criterion, based on Irwin's approach, is verified on each nonlinear iteration. When the criterion is satisfied, additional contact elements are added into finite element and discrete fracture model formulations respectively. The proposed approach allows modeling of existing natural and artificially created fractures within one framework. The model is tested on single‐ and multiple‐phase fluid flow examples for both isothermal and thermal conditions and verified against existing semianalytical solutions. The applicability of the approach is demonstrated on an example of practical interests where a sector model of an oil reservoir is simulated with different injection and production regimes.     

Thermodynamic analysis of heat transfer in a wellbore combining compressed air energy storage

Compressed air energy storage (CAES) is a potential energy storage technology. The gas phase and short cycle period are two key factors affecting heat transfer loss in the wellbore of CAES. A semi-analytical solution was developed by using the convolution method considering gas movement in this study to describe the transient behavior of heat transfer with a short cycle period. The comparative analysis of the presented solution with two published solutions showed that the solution matched well with the previous solutions under steady state. Parametric studies were carried out to investigate the impact of injection rate, overall heat transfer coefficient and thermal diffusivity of the formation on heat loss in the wellbore. The results indicated that a low overall heat transfer coefficient and thermal diffusivity of the formation with an appropriate injection rate can efficiently reduce the heat loss. A hypothetical case study with a short cycle period of injection and production was conducted to demonstrate the applicability of the developed solution in CAES. The results suggest that the semi-analytical solution is applicable for heat transfer in the wellbore of CAES.     

岩石导热热阻对裂隙对流换热的影响机制

增强型地热系统(EGS)中高温岩石与流体之间的对流换热特征一直以来是干热岩(HDR)研究的重要基础内容。岩石导热热阻对裂隙对流换热特征具有重要影响。为研究其具体影响,综合运用理论解析与数值模拟2种研究方法,通过对解析解讨论以及建立数值模型,研究两平行光滑平板之间的换热规律。结果表明：流体速度、传热边界层充分发展时,局部努塞尔特准数Nu_x为定值,与其他因素无关;局部对流换热系数h_x仅与流体热导率k和裂隙开度e有关,与其他因素无关。上下平板壁面热流恒定时,Nu_x为8.235;温度恒定时,Nu_x为7.54。然后建立多组导热热阻不同的岩石裂隙对流换热数值模型,发现岩石导热热阻增大,温度场进口段延长,对流换热系数h增大。岩石长度显著影响进口段占比,进而影响h的大小。h随着长度增大而减小;当岩石长度足够长时,进口段占比足够小,此时除k与e之外的参数对h基本没影响。并且发现实验室常用岩石长度为100 mm,而典型EGS工程中裂隙长度是米级的,建议室内实验重视岩石长度对裂隙对流换热特征的影响。     

挤密条件下U型地埋管换热器换热效率的理论分析及数值模拟

土层的导热系数是影响地源热泵地埋管换热器换热效率的重要因素之一,与土样的密度密切相关。为了提高地源热泵换热器换热效率,本文基于柱热源理论,建立了挤密钻进条件下换热器的传热理论模型,对换热器的换热效率进行了计算,并利用Ansys有限元软件对土层挤密条件下换热器的换热效果进行了数值分析。结果表明,与传统方法相比,土层挤密可有效降低钻孔周围的热阻,U型管换热器的换热效率可提高17%~20%。     

地热井U 型管周围土壤温度场的ANSYS 模拟

埋管换热器和岩土的传热是个复杂过程,利用有限元仿真分析两者的传热过程是有效的办法。通过与单井二维模型的对比分析,建立二维群井传热模型,模型为9 口间距为4 m 的储热井,模型传热过程遵循经典传热理论。按照ANSYS 网格划分原则,将求解区域划分为最大边长0. 1 m,最小边长 0. 01 m 的单元。在地埋管内侧加载热流密度,外边界加载地温边界条件,经过模拟72 h 的储热过程,精确得出群井周围岩土体的温度场变化规律,求解区域最高温度为40. 263℃,最低温度为10℃,影响半径约1. 4 m,温度的升高幅度会减缓,并得到了温度及热量迁徙的规律。模拟结果与实际测量结果相符。     

The flow and heat transfer characteristics of compressed air in high-pressure air injection wells

High-pressure air injection (HPAI) is a significant enhanced oil recovery (EOR) technology of light oils especially in deep, thin, low-permeability reservoirs. The flow and heat transfer behaviors of compressed air in wellbore is essential to maximize performance of air in EOR. Due to strong compressibility of air and high injection pressure, wellbore temperature and pressure are greatly affected by friction and gas compression. However, the available models of wellbore flow and heat transfer are only accurate for thermal fluid, such as saturated steam and superheated steam, injected at relatively low pressure and high temperature. In this paper, a novel model is proposed to characterize wellbore pressure and temperature distribution for HPAI wells with consideration of dynamic behaviors of injected air. Flow and heat transfer in depth direction are coupled with air properties by iterative technique, and heat transfer in radial direction is treated as steady state in wellbore and transient state in formation. The mathematical model is solved by employing finite difference method and it is validated by field data. Then, integrated analyses of flowing pressure, heat transfer mechanism, and interaction between pressure and temperature are conducted. Results indicate that (1) as well depth increases, temperature difference between formation and air tends to become constant, and the radial heat transfer tends to reach an equilibrium state. The higher the flow rate is, the deeper the equilibrium depth is. (2) Air temperature is dominated by heat transmission from formation at low flow rates and dominated by frictional heat and gas compression effect at high flow rates. Fictional heat begins to affect air temperature at an injection rate beyond the critical value, while gas compression effect can increase air temperature in the whole calculated injection rate range. (3) Interaction between wellbore temperature and pressure is mainly achieved by altering air density. The effect of injection pressure on air temperature can be negligible, while the influence of injection temperature shows strong dependency on injection rate.     

Free convection along a vertical wall in a porous medium with periodic permeability variation

This communication investigates the effect of permeability variation on free convection flow and heat transfer in a porous medium bounded by a vertical porous wall. A transverse periodic variation in permeability is considered to study the effect on heat transfer rate and skin friction. The problem becomes three-dimensional due to the variation in permeability. Analytical expressions for velocity, temperature, skin friction and rate of heat transfer are obtained using the perturbation technique. Effects of permeability and suction parameters on skin friction and heat transfer is studied and it is found that the increase in permeability creates a higher heat transfer rate and greater skin friction. However, the increase in suction parameter causes the skin friction values to decrease. The results obtained will be useful in the design of steam displacement processes in oil recovery and various geothermal systems.     

水合物钻井液微观结构热传导分形模型研究

[摘 要] 在水合物勘探开发过程中,钻井液除常规作用外,还必须具有抑制水合物分解与再生成 的作用。钻井液的这些宏观作用由其内部微观结构所决定,而水合物的分解和生成与介质的传热过程 有关。由于钻井液微观结构与多孔介质相似,故其传热过程也相似。利用多孔介质的性质及其传热理 论,系统分析了水合物钻井液微观结构与多孔介质的异同,探讨了水合物在钻井液中的分解特性。接着 阐述了分形理论,对钻井液微观结构的传热分形模型进行了研究。最后得出钻井液中水合物分解的传 热过程,建立了钻井液微观结构的一维热传导分形模型,并得出有效导热系数的分形表达式。     

桩埋管参数对渗流下能量桩热-力耦合特性的影响

为获得地下水渗流作用下桩埋管参数对能量桩热-力耦合特性的影响,建立了不同埋管参数的能量桩数值模型,分析了桩埋管数量、埋管布置形式、埋管管径对单位桩深换热量、日换热量、桩截面平均温升、桩身位移增量及桩身附加温度荷载的影响。结果表明：增加埋管数量可以增大能量桩换热量,但也会加剧桩内不同埋管间的热干扰,导致换热性能下降及桩身...     

实验室尺度高地温梯度模拟地层的实现方法研究

地热井与周围热储层的传热过程对地热井产热性能研究有重要意义。由于实际工程中在地热井周边布置测点较难,无法获取地热井周围热储层的参数变化,进而为地热热储模拟结果提供验证,故以往大多地热热储模拟仅将地热井作为源项处理,未考虑地热流体和储层的耦合流动换热。实验室条件下的模拟试验方便布置测点,可为热储-井筒耦合流动传热模型提供试验验证,其中如何实现实验室尺度下有温度梯度的模拟地层是试验研究的关键,目前尚未有类似研究。基于传热学基本原理,研究了实验室条件下有较高温度梯度多孔地层的快速实现方法,通过确定模拟热储层和热储盖层几何尺寸、优选填充多孔介质和实现恒定温度的模拟热储层,设计了一套实验室尺度下有高温度梯度的模拟地层系统,通过分层加热与边界动态热补偿方法,较快实现了热储层温度分别为60,65,70°C下模拟地层的线性温度分布,采用有限体积法得到的数值模拟与试验结果的相对误差在±2.5%范围内,二者吻合较好。文章设计搭建的模拟地层系统可为开展地热井筒-热储耦合模拟试验提供条件,进而为开发的地热热储-井筒耦合传热数值软件提供试验验证。     

不同土壤类型与含水率对水平埋管换热性能影响数值分析

为揭示地源热泵系统水平埋管换热器在不同土壤类型中的换热性能,基于土壤毛管水理论知识,结合数值模拟的研究手段,探讨了蓄能不同类型土体内(砂土、壤土、黏土)三相组成的差异对水平埋管换热器换热特性的影响规律。结果表明,在通入308.15 K制冷工况下,水平管在壤土中的出水温度降低至303.3 K,进出口水温差为4.9 K,埋管单位延米换热量37.1 W/m,水平管在壤土中的制冷换热效益显著;不同土壤(砂土、壤土、黏土)在经历相同制冷周期下,水平管的换热过程对壤土的温度场分布影响最小,管体在壤土中运行时热堆积风险系数最低。研究表明,水平管与土壤的换热性能同时受土壤比热容与土壤导热系数的影响,提高土壤导热系数比提高土壤比热容获得的效益更加显著。可以通过压实回填、减少土壤孔隙率、提高固相回填材料导热系数、加大布管深度以提高回填材料含水率等方法来强化埋管的换热性能。