An integral equation solution for three‐dimensional heat extraction from planar fracture in hot dry rock

In the numerical simulation of heat extraction by circulating water in a fracture embedded in geothermal reservoir, the heat conduction in the reservoir is typically assumed to be one‐dimensional and perpendicular to the fracture in order to avoid the discretization of the three‐dimensional reservoir geometry. In this paper we demonstrate that by utilizing the integral equation formulation with a Green's function, the three‐dimensional heat flow in the reservoir can be modelled without the need of discretizing the reservoir. Numerical results show that the three‐dimensional heat conduction effect can significantly alter the prediction of heat extraction temperature and the reservoir life as compared to its one‐dimensional simplification. Copyright © 2003 John Wiley & Sons, Ltd.     

Integral equation solution of heat extraction‐induced thermal stress in enhanced geothermal reservoirs

During fluid injection in enhanced geothermal systems, thermo‐mechanical processes can play an important role. In fact, the phenomena of reservoir seismicity and the variation of injectivity with respect to injection water temperature can be attributed to the induced thermal stresses. In this paper, a three‐dimensional integral equation formulation is presented for calculating thermally induced stresses associated with the cooling of a fracture in a geothermal reservoir. By utilizing Green's function in the integral equation, the three‐dimensional heat flow and stresses in the reservoir are modelled without discretizing the reservoir. The formulation is implemented in a computer program for the solution of injection into an infinite fracture as well as for the injection/extraction in an arbitrarily shaped fracture. Copyright © 2005 John Wiley & Sons, Ltd.     

A three‐dimensional integral equation model for calculating poro‐ and thermoelastic stresses induced by cold water injection into a geothermal reservoir

Poro‐mechanical and thermo‐mechanical processes change the fracture aperture and thus affect the water flow pattern in the fracture during the cold water injection into enhanced geothermal systems (EGS). In addition, the stresses generated by these processes contribute to the phenomenon of reservoir seismicity. In this paper, we present a three‐dimensional (3D) partially coupled poro‐thermoelastic model to investigate the poroelastic and thermoelastic effects of cold water injection in EGS. In the model, the lubrication fluid flow and the convective heat transfer in the fracture are modeled by the finite element method, while the pore fluid diffusion and heat conductive transfer in the reservoir matrix are assumed to be 3D and modeled by the boundary integral equation method without the need to discretize the reservoir. The stresses at the fracture surface and in the reservoir matrix are obtained from the numerical model and can be used to assess the variation of in situ stress and induced seismicty with injection/extraction. Application of the model shows that rock cooling induces large tensile stresses and increases fracture conductivity, whereas the rock dilation caused by fluid leakoff decreases fracture aperture and increases compressive total stresses around the injection zone. However, increases in pore pressure reduce the effective stresses and can contribute to rock failure, fracture slip, and microseismic activity. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical Modeling and Spatial Moment Analysis of Thermal Fronts in a Coupled Fracture-Skin-Matrix System

In this study, the behavior of thermal fronts along the fracture is studied in the presence of fracture-skin in a coupled fracture-matrix system. Cold water is injected into the fracture, which advances gradually towards production well, while extracting heat from the surrounding reservoir matrix. The heat conduction into the fracture-skin and the rock-matrix from the high permeability fracture is assumed to be one dimensional perpendicular to the axis of fluid flow along the fracture. Constant temperature cold water is injected through an injection well at the fracture inlet. The fluid flow takes place along the horizontal fracture which ensures connectivity between the injection and production wells. Since the rock-matrix is assumed to be tight, the permeability of fracture-skin as well as the rock-matrix is neglected. The present study focuses on the heat flux transfer at the fracture-skin interface as against the earlier studies on fracture-matrix interface, and the sensitivity of additional heterogeneity in the form of fracture skin in a conventional fracture-matrix coupled system is studied. The behavior of thermal fronts for various thermal conductivity values of the fracture-skin and rock-matrix is analyzed. Spatial moment analysis is performed on the thermal distribution profiles resulting from numerical studies in order to investigate the impact on mobility and dispersion behavior of the fluid in the presence of fracture-skin. The presence of fracture skin affects the heat transfer significantly in the coupled fracture-matrix system. The lower order spatial moments indicate that the effective thermal velocity increases with increase in skin thermal conductivity and a significant thermal dispersion is observed at the inlet of the fracture owing to the high thermal conductivity of the fracture-skin at the early stages. Furthermore the higher spatial moments indicate that the asymmetricity increases with decrease in skin thermal conductivity unlike the case with half fracture aperture and fluid velocity and the kurtosis is maximum with higher skin thermal conductivity which implies enhanced heat extraction from the fracture-skin into the fracture. Results suggest that the amount of heat extraction by the circulating fluid within the fracture from the reservoir not only depends on the rock-matrix module of the reservoir characteristics but also the fracture-skin characteristics of the system and subsequently influence the reservoir efficiency.     

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.     

Integral equation solution for three-dimensional heat transfer in multiple-fracture enhanced geothermal reservoirs

For the prediction of energy production from multiple-fractured geothermal reservoirs, previous models basically focused on the one-dimensional conduction in the rock containing evenly distributed fractures of equal scale. Here, a novel model is described to numerically investigate the three-dimensional heat transfer in geothermal reservoirs with unevenly spaced disc fractures of various sizes including the aperture and radius. In terms of the water flow through each fracture, an approximate analytical solution is obtained on the assumption that the water pressure disturbances, induced by the fracture margin and extraction (injection) operation, at the injection (extraction) well center and at different locations within the injection (extraction) well range were approximately equal. By the integral equation scheme for two-dimensional planar fractures, the three-dimensional problem of heat exchange is simulated without the reservoir discretization. The singular integral is analytically calculated in polar coordinates whereas the nonsingular integrand is numerically estimated by the Gaussian quadrature method in Cartesian coordinates. Compared with the one-dimensional simplification, the three-dimensional heat conduction remarkably alters the prediction of extraction temperature. In addition, the reservoir temperature field is also significantly influenced by the spacings and dimensions of fractures. The present model may be used for the estimation, design, and optimization of a geothermal reservoir.     

Seepage-heat transfer coupling process of low temperature return water injected into geothermal reservoir in carbonate rocks in Xian County,China

Fracture seepage and heat transfer in the geothermal reservoir of carbonate rocks after the reinjection of low temperature geothermal return water is a complex coupling process,which is also the frontier of geothermal production and reinjection research. Based on the research of cascade comprehensive development of geothermal resources in Beijing-Tianjin-Hebei (Xian County),the carbonate geothermal reservoir of Wumishan formation in the geothermal field in Xian County is investigated. With the development of the discrete fracture network model and the coupling model of seepage and heat transfer,the numerical solution of seepage field and temperature field with known fracture network is reached using the finite element software COMSOL,and the coupling process of seepage flow and heat in carbonate rocks is revealed. The results show that the distribution of temperature field of fractured rocks in geothermal reservoir of carbonate rocks has strong non-uniformity and anisotropy. The fracture network is interpenetrated,which constitutes the dominant channel of water conduction,and along which the fissure water moves rapidly. Under the influence of convective heat transfer and conductive heat transfer,one of the main factors to be considered in the study of thermal breakthrough is to make the cold front move forward rapidly. When the reinjection and production process continues for a long time and the temperature of the geothermal reservoir on the pumping side drops to a low level,the temperature of bedrocks is still relatively high and continues to supply heat to the fissure water,so that the temperature of the thermal reservoir on the pumping side will not decrease rapidly to the water temperature at the inlet of reinjection,but will gradually decrease after a long period of time,showing an obvious long tail effect. The distribution of fractures will affect the process of seepage and heat transfer in carbonate reservoirs,which should be considered in the study of fluid thermal coupling in carbonate reservoirs.     

Late-Stage Reservoir Formation Effect and Its Dynamic Mechanisms in Complex Superimposed Basins

<正>Complex superimposed basins exhibit multi-stage tectonic events and multi-stage reservoir formation;hydrocarbon reservoirs formed in the early stage have generally late-stage genesis characteristics after undergoing adjustment,reconstruction and destruction of later-stage multiple tectonic events.In this paper,this phenomenon is called the late-stage reservoir formation effect.The late-stage reservoir formation effect is a basic feature of oil and gas-forming reservoirs in complex superimposed basins,revealing not only multi-stage character,relevance and complexity of oil and gas-forming reservoirs in superimposed basins but also the importance of late-stage reservoir formation. Late-stage reservoir formation is not a basic feature of oil and gas forming reservoir in superimposed basins.Multi-stage reservoir formation only characterizes one aspect of oil and gas-forming reservoir in superimposed basins and does not represent fully the complexity of oil and gas-forming reservoir in superimposed basins.We suggest using "late-stage reservoir formation effect" to replace the "late-stage reservoir formation" concept to guide the exploration of complex reservoirs in superimposed basins. Under current geologic conditions,the late-stage reservoir formation effect is represented mainly by four basic forms:phase transformation,scale reconstruction,component variation and trap adjustment.The late-stage reservoir formation effect is produced by two kinds of geologic processes: first,the oil and gas retention function of various geologic thresholds(hydrocarbon expulsion threshold,hydrocarbon migration threshold,and hydrocarbon accumulating threshold) causes the actual time of oil and gas reservoir formation to be later than the time of generation of large amounts of hydrocarbon in a conventional sense,producing the late-stage reservoir formation effect;second, multiple types of tectonic events(continuously strong reconstruction,early-stage strong reconstruction, middle-stage strong reconstruction,late-stage strong reconstruction and long-term stable sedimentation) after oil and gas reservoir formation lead to adjustment,reconstruction and destruction of reservoirs formed earlier,and form new secondary hydrocarbon reservoirs due to the late-stage reservoir formation effect.     

Early Palaeozoic carbonate reservoirs from the Yingshan Formation of Well block ZG‐43 in Tazhong Low Rise,Central Uplift,Tarim Basin,NW China: geological features and controlling factors

Marine carbonate reservoirs, as a focus of petroleum exploration and development in China, are involved with high exploration risk and prediction difficulty owing to high heterogeneity and diversity of reservoir beds. In the Tarim Basin, NW China, carbonate reservoirs host about 38% of the whole basin's hydrocarbon resources in a large prospecting area mainly distributed in the Cambrian and Ordovician rocks in central (Tazhong) and northern (Tabei) Tarim. Recently, a better understanding has been made of the karsted weathering crust at the top of the Lower Ordovician Yingshan Formation in the northern slope area of the Tazhong Low Rise, Central Uplift, Tarim Basin. As a new frontier of exploration, oil/gas distribution and controlling factors of carbonate reservoirs in the Yingshan Formation are not clearly understood. In this work, we investigated the reservoir beds and oil/gas properties in 13 wells in Well block ZG‐43 on the No. 10 structural belt in the Tazhong Low Rise, and studied hydrocarbon accumulation characteristics with seismic and geochemical data. The Yingshan Formation in Well block ZG‐43 is mainly composed of calcarenite, dolomitic limestone, dolomite, cryptite, as low porosity and low permeability reservoir beds, with fracture‐void porosity constituting the main reservoir pore space. Oil/gas is quasi‐layer distributed beneath the unconformity between the Yingshan and Lianglitag formations to a depth of 140 m. The oil in Well block ZG‐43 is condensate with low density, low viscosity, low sulphur, low resin, low asphaltene, and high wax. The gas is 87.3% methane, generally containing H₂S. The oil/gas distribution pattern is oil in the east and gas in the west, and H₂S content in the west is lower than that in the east. The controlling factors for hydrocarbon are multi‐source supply and multi‐phase charging, interstratal karstification, hydrothermal activity, structural location and sealing condition. A structural–lithological trap is the main type of oil/gas accumulation. Oil/gas distribution was clearly affected by strike–slip faults. Oil/gas with multi‐source supply and multi‐phase charging was controlled by favourable local palaeo‐highs, and affected by later karsting and hydrothermal activity, as well as gas invasion in the Himalayan (Cenozoic) period. Under the caprock of compact limestone in the third to fifth members of the Lianglitag Formation, oil/gas migrated up along the strike–slip fault planes, and moved laterally to both sides in a ‘T’ shape, and formed large‐scale quasi‐layer condensate gas reservoirs controlled by reservoir bed quality. Copyright © 2013 John Wiley & Sons, Ltd.     

米尺度裂隙岩体模型水流-传热试验的数值模拟分析

为了研究高放射性核废物地下处置库近场的水流-传热耦合问题,采用国内高放废物地下处置库预选场址--甘肃北山地区的花岗岩石块体,加工组合成米尺度的规则裂隙岩体模型,设置边界热源和裂隙水流,试验模拟裂隙水水流与传热之间的相互作用。作为该室内模型试验的前期理论研究,采用等效孔隙介质数值模型,着重分析了裂隙开度、裂隙流量和热源功率对流场和温度场的影响。在设定条件下,计算分析表明：热传导和裂隙水水流由热源作用初期的不耦合很快转化为耦合;不流动的裂隙水主要表现为热存储和热传导,而流动的裂隙水还引起流动传热和水与岩石之间的对流换热,使岩体温度场明显不同于单纯热传导的情况;如果保持裂隙水流量不变,则裂隙开度的变化对水流-传热影响不大;如果保持裂隙水流速不变,则裂隙开度的变化对水流-传热影响显著;热源功率越大,通过裂隙水的热流量越大,裂隙水压强越大,而当温度超过100 ℃时,裂隙水会因汽化而压强显著增大;加热7 d时,热量的输入和输出几乎相等,裂隙水流带走的热量接近热源供给的热量,模型系统基本达到了热平衡。     

准噶尔盆地西北缘二叠系火山岩储层裂缝发育特征及分布预测——以金龙2井区佳木河组为例

准噶尔盆地西北缘金龙2井区二叠系佳木河组裂缝是该区火山岩储层油气主要的渗流通道。综合岩心、岩石薄片及成像测井等资料,识别出该区主要发育的裂缝类型为半充填或未充填高角度缝,其次为半充填低角度斜交缝与网状缝。成像测井解释裂缝方位近东西向,与岩心古地磁解释现今地应力最大主应力方向近似平行,有效性开启较好。火山岩储层裂缝发育主要受构造与岩性两种因素影响。距离断层越近,由于构造曲率增大,裂缝越发育,裂缝多沿断裂呈条带状分布。不同的火山岩类型,裂缝发育程度也不同。通过成像测井资料分析,认为研究区中-酸性火山熔岩及火山碎屑熔岩裂缝较发育,并进一步定量计算出单井裂缝密度、裂缝倾角、裂缝孔隙度等,确定单井裂缝发育特征。结合叠前地震预测方法,即叠前方位各向异性法(AVAZ),优选衰减起始频率属性,预测了佳木河组火山岩储层裂缝分布特征。     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermal effects on fluid flow and hydraulic fracturing from wellbores and cavities in low-permeability formations

The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to low permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimize a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. Copyright © 1999 John Wiley & Sons, Ltd.     

Smooth particle hydrodynamics simulations of low Reynolds number flows through porous media

In this paper, a three‐dimensional smooth particle hydrodynamics (SPH) simulator for modeling grain scale fluid flow in porous media is presented. The versatility of the SPH method has driven its use in increasingly complex areas of flow analysis, including the characterization of flow through permeable rock for both groundwater and petroleum reservoir research. SPH provides the means to model complex multi‐phase flows through such media; however, acceptance of the methodology has been hampered by the apparent lack of actual verification within the literature, particulary in the three‐dimensional case. In this paper, the accuracy of SPH is addressed via a comparison to the previously recognized benchmarks of authors such as Sangani and Acrivos (Int. J. Multiphase Flow 1982; 8 (4): 343–360), Zick and Homsy (J. Fluid Mech. 1982; 115 :13–26) and Larson and Higdon (Phys. Fluids A 1989; 1 (1):38–46) for the well‐defined classical problems of flow through idealized two‐ and three‐dimensional porous media. The accuracy of results for such low Reynolds number flows is highly dependent on the implementation of no‐slip boundary conditions. A new, robust and numerically efficient, method for implementing such boundaries in SPH is presented. Simulation results for friction coefficient and permeability are shown to agree well with the available benchmarks. Copyright © 2010 John Wiley & Sons, Ltd.     

稀疏不规则裂隙岩体模型三维水流-传热对温度和应力影响的试验研究

选取高放射核废物处置库重要预选场区甘肃北山地区的花岗岩,制作750 mm（宽）×300 mm（厚）×1 000 mm（高）的稀疏不规则裂隙岩体模型,该模型由18块花岗岩和竖向与斜向各两条裂隙组成,在裂隙及岩石内部埋置温度传感器、水压计、直角应变花,并在模型一侧设置局部热源,研究热源温度和裂隙水流速对岩石温度和应力的影响。结果表明,竖裂隙水主要从顶部进水口流向底部出水口,斜裂隙水主要从侧部进水口流向侧部出水口,竖裂隙与斜裂隙在交汇处存在微小流量交换;由于热源处在两条斜裂隙进水口之间,并且斜裂隙长度小于竖裂隙,岩石热传导与斜裂隙水流对岩石温度分布起控制作用,竖裂隙水流对岩石横向热传导起阻滞作用;由于热传导和水流传热的不规则性,上层岩石形成从左向右为主的传热路径,中层和下层岩石形成从上向下为主的传热路径;由于上、下层岩石温度梯度较小,岩石收缩受热拉应力,而中层岩石温度梯度较大,岩石膨胀受热压应力,大主应力的方向大致垂直于斜裂隙面与竖裂隙面的交线,岩石应力增量随斜平面方向的温度梯度增大而增大;热源温度越高,裂隙水流速越低,岩石温度越高、岩石应力越大,系统达到稳态需要的时间越长。     

塔河油田12区块奥陶系裂缝分布规律研究

通过对露头、铸体薄片、岩性分析和成像测井资料观察和分析发现,塔河油田12区块奥陶系储层内缝洞发育,并且裂缝在一间房和鹰山组的储层中起了至关重要的作用,它既是地层流体的存储空间,又是连通孔洞缝的主要通道。裂缝的发育也增加了储层的非均质性。所以,精细研究裂缝的发育和分布规律,并定量地计算裂缝的相关参数对于识别奥陶系储层和对...     

饱和稀疏裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法

针对高放射核废深地质处置库近场环境,建立分布热源作用下饱和裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法：采用Goodier热弹性位移势和Laplace变换计算由温度梯度产生的温梯位移和应力;考虑单一裂隙的情况,利用经典弹性力学的Boussinesq解和Cerruti解计算为满足边界条件的约束位移和应力,与温梯位移和应力叠加,可得总体热位移和应力;把裂隙面离散为矩形单元集合,采用极坐标系下的解析法计算包含奇点的单元积分,采用数值法计算与分布热源有关和不含奇点的单元积分。与基于裂隙面法向一维热传导假设的一种解析解对比,结果表明,半解析法与解析法的计算结果基本一致,但由于半解析法考虑岩石的三维热传导,因温度时空分布和演变的不同而导致不同的温梯应力。针对一个假想单裂隙岩体三维水流-传热过程,计算温梯位移和应力、约束位移和应力、总体位移和应力;结果表明,裂隙水流-传热可能对位移和应力的分布和演变有显著影响,距离分布热源较近的岩石因升温膨胀受到约束而出现压应力,而距离分布热源较远的岩石则可能因协调收缩受到约束而出现拉应力。     

分布热源作用下裂隙岩体渗流-传热的拉氏变换-格林函数半解析计算方法

以裂隙岩体高放射性核废物地下处置库性能评估为目标,提出了分布热源作用下单裂隙岩体渗流-传热的简化概念模型、控制微分方程和拉氏变换-格林函数半解析法,为进一步采用半解析法计算分布热源作用下多裂隙岩体的渗流-传热问题奠定了基础。针对单裂隙岩体的渗流-传热问题,建立考虑岩石内热源和二维热传导的控制微分方程,利用拉氏变换域微分方程的基本解建立格林函数积分方程,采用解析法处理其中的奇点,通过数值积分和拉氏数值逆变换求解,计算任意时刻裂隙水和岩石的温度分布。通过算例,与基于岩石一维热传导假定的解析解进行了对比,并计算分析了分布热源作用下单裂隙岩体的渗流-传热特征及其对裂隙开度、岩石热传导系数和热流集度的敏感度。算例表明,（1）就裂隙水温度而言,由于考虑了岩石的二维热传导,拉氏变换-格林函数半解析解小于基于岩石一维热传导假定的解析解;（2）裂隙水温度和岩石温度对裂隙开度和热流集度的敏感度较大,对岩石热传导系数的敏感度较小。     

分布热源作用下单裂隙岩体三维水流-传热的半解析计算方法

针对高放射性核废物地下处置库近场饱和裂隙岩体环境,提出一种由分布热源、饱和单裂隙和两侧无限大岩石构成的三维水流-传热简化模型,建立了控制微分方程和基于拉氏变换域格林函数的积分方程;采用矩形单元把裂隙面域离散化,利用极坐标下的解析方法计算包含奇点的单元积分,利用数值方法计算分布热源和不包含奇点的单元积分,建立拉氏变换域的线性代数方程组,求解后,利用拉氏数值逆变换,计算任意时刻裂隙水和岩石的温度分布。对两个无内热源、流场确定的计算模型进行了计算,与仅考虑岩石沿裂隙面法向一维热传导的解析解进行了对比。计算分析了分布热源作用下饱和单裂隙岩体的三维水流-传热特征及其对裂隙水流速、岩石热传导系数和热源热流集度的敏感度。计算结果表明：与直接采用高斯数值积分相比,提出的解析法奇异积分精度较高;就裂隙水温度而言,单裂隙岩体三维水流-传热半解析计算方法与解析法得到的结果基本一致,但由于半解析计算方法考虑了岩石的三维热传导,使得裂隙水的上游温度较低,而下游温度较高;无分布热源作用时,岩石热传导系数越大,裂隙水温度越低;裂隙水流速越大,裂隙进水温度对裂隙水和岩石温度分布的影响越明显;由于受到裂隙水流动传热的作用,分布热源对裂隙水温度和岩石温度的影响在裂隙水流的下游区域比较显著。     

共和盆地干热岩体人工裂隙带结构的控热机理与产能优化

人工压裂是获取干热岩型地热资源的关键环节,压裂后的人工裂隙带结构对开采条件下水热传递过程具有重要控制作用。结合我国共和盆地干热岩储层地质条件,采用数值模拟方法着重分析干热岩不同产状人工裂隙带的渗透率与宽度对热储中水热传递过程的影响机理,明确不同人工裂隙结构条件下水热产出能力,进而优化井间距。结果表明:当人工裂隙带渗透率较小时（小于5 D）,裂隙带规模越大,开采井温度越高;当渗透率较大时（大于10 D）,在水平裂隙带中,随着裂隙带规模的增加,由于注入冷水的快速扩散导致整体低温区域增加,开采井温度反而降低。在水平裂隙带中注入冷水主要为水平向流动,随着渗透率的增加,开采井温度更易受注入冷水的影响而降低;但在垂直裂隙带及倾斜裂隙带中,随着渗透率的增加,垂向自由对流增强,注入冷水更易于向储层底部高温区域流动,经加热后到达开采井,使得开采温度提升。综合比较,同一井间距条件下,低渗水平裂隙带以及高渗垂直裂隙带的产热能力较其他裂隙带更强。