Terrestrial heat flow and geothermal gradients in relation to hydrodynamics in the Alberta basin, Canada

Geothermal gradients in the Alberta part of the Western Canadian sedimentary basin have been studied on the basis of 55,244 bottom-hole temperature values from 28,260 petroleum exploration wells. Gradient estimates for different depth and stratigraphic intervals together with a study of the heat conductivity distribution indicate both regional heat flow variations and variations with depth. The regional hydrodynamics of the basin strongl influences both grad ifT gradient and heat flow increase with depth in water recharge areas to the west and decrease with depth in discharge areas to the north and east. The results indicate that heat flow in the central part of the basin should be approximately equal to the deep crustal heat flow.     

Geothermal gradients and heat flow in the Beaufort-Mackenzie Basin,Arctic Canada

Thermal gradients have been calculated and heat flow estimates made for 34 petroleum exploration wells along four regional profiles crossing the Mesozoic-Cenozoic Beaufort-Mackenzie Basin of northern Canada. The geothermal gradients vary from 22 mKm^–1 to 44 mKm^–1. Four sets of possible thermal conductivity values were used to calculate a range of heat flow values for each well. Generally low heat flow is observed in wells within the deeper portions of the basin and higher heat flow values occur in wells along the Aklavik Arch Complex which forms the southeastern margin of the basin.The contribution to heat flow by heat generation below the Mesozoic-Cenozoic basin fill sediments has been considered. The heat flow contribution from sub-Mesozoic sedimentary strata and underlying basement is highest along the basin-bounding Aklavik Arch Complex. The decrease in heat flow from below the basin fill sediments toward the basin depocenter may be related to basinward crustal thinning and corresponding reductions in intra-crustal radiogenic heat production.     

A geothermal study of the Jharia Gondwana basin (India): Heat flow results from several holes and heat production of basement rocks

Temperatures have been measured in nine boreholes (ranging from 400 to 900 m in depth) in the Jharia Gondwana sedimentary basin of the Indian shield. About two hundred thermal conductivity determinations have been made on core samples from these holes. Temperature profiles, Bullard plots and heat flow profiles of these holes indicate different types of disturbances in the shallow geothermal regime, attributable principally to groundwater movement. Heat flow in the region of the “anticlinal high” is about 0.4 HFU higher than the heat flow in the main synclinal region of the basin. The possible sources for this variation are regional groundwater movement and upwelling of thermal waters through a deep-seated fault/fracture system. The heat flow of 1.9 HFU characterizing the main synclinal region, taken as the regional value for the basin as a whole, has been related to the heat generation of the Precambrian basement rocks. A plot of heat flow vs. heat generation falls in line with three plots for the Precambrian complexes of the Indian shield, indicating the absence of a thermal anomaly due to deeper crustal conditions underneath this basin.     

The geothermal regime of the northern Yukon and Mackenzie delta regions of northwest Canada — studies of two regional profiles

Temperature data from deep petroleum exploration wells and thermal conductivity estimates based on net rock analysis data have been used to make terrestrial heat flow estimates along two profiles across the sedimentary strata of the Mackenzie Delta, northern Yukon, and offshore Beaufort Sea regions.Both profiles exhibit low heat flow values that range from 34 mWm^–2 to 58 mWm^–2, and little change occurs over large distances in the continental part of the area. Low heat flow values (<40 mWm^–2) occur in the Beaufort-Mackenzie Basin and Rapid Depression, both of which are areas of thick successions of Cretacecus and Tertiary clastic sedimentary strata. High heat flow values of almost 80 mWm^–2 occur to the south in the Taiga Nahoni Foldbelt and values as high as 60 mWm^–2 are indicated along the Aklavik Arch Complex, northeast of Aklavik.The regional variations of effective thermal conductivity are insufficient to account for the heat flow variations along the profiles, and so these may indicate deep radiogenic or other heat sources.     

辽河裂谷盆地地幔热流

辽河盆地是一个在前中生代基底上发展起来的裂谷盆地,实测大地热流平均值为65mW/m2,变动于44—83mW/m2之间。在给出地壳结构模型并确定各岩层放射性生热率的基础上,采用“剥层”法从地表开始,自上而下,由浅及深地扣除各岩层所提供的热量,从而得出地幔热流值。结果表明,辽河裂谷盆地地幔热流为41mW/m2,占整个地表总热流量的63％。可见,本区热量大部分来自地幔。与世界上其它地质构造单元相比,辽河裂谷盆地无论地幔热流绝对值或其与地表热流之比值,都具有介于稳定地区和构造活动区之间的特点。作者认为,辽河裂谷盆地地幔热流的上述特点,乃是中、新生代以来本区长期地质历史发展的产物。     

Three-dimensional structure of the Laguna Salada Basin and its thermal regime

A comprehensive reinterpretation of the available gravity, magnetic, geothermal, geological and borehole information has been made of the Laguna Salada Basin to establish a 3D model of the basement and sedimentary infill. According to statistical spectral analysis, the residual gravity anomaly is due to sources with a mean regional depth of 2.8 km. The topography of the basement was obtained from a three‐dimensional inversion carried out in the wavenumber domain using an iterative scheme. The maximum density contrast of ?300 kg/m³ estimated from previous studies and the mean depth of 2.5 km finally constrained this inversion. The resulting model indicated that the sedimentary infill is up to 4.2 km thick at its deepest point. According to the gravity‐derived basement topography, the basin presents an asymmetry (i.e. it is of the half‐graben type). It is deeper to the east, where it is delimited from the Sierra Cucapah by a step fault. By contrast, the limit with the Sierra de Juarez is a gently sloping fault (i.e. a listric fault). The basement is not even, but it comprises a series of structural highs and lows. N–S to NW–SE and E–W to NE–SW faults delimit these structural units. The magnetic modelling was constrained by (i) the gravity‐derived basement topography; (ii) a Curie isotherm assumed to be between 7 km and 10 km; (iii) assuming induced magnetization only; (iv) the available geological and borehole information. The magnetic anomalies were interpreted successfully using the gravity‐derived basement/sedimentary interface as the top of the magnetic bodies (i.e. the magnetic modelling supports the gravity basement topography). An elongated N–S to NW–SE trending highly magnetized body running from south to north along the basin is observed to the west of the basin. This magnetic anomaly has no gravity signature. Such a feature can be interpreted as an intrusive body emplaced along a fault running through the Laguna Salada Basin. Treatment of the gravity and magnetic information (and of their horizontal gradients) with satellite image processing techniques highlighted lineaments on the basement gravity topography correlating with mapped faults. Based on all this information, we derived detailed geological models along four selected profiles to simulate numerically the heat and fluid flow in the basin. We used a finite‐difference scheme to solve the coupled Darcy and Fourier differential equations. According to our results, we have fluid flow in the sedimentary layers and a redistribution of heat flow from the basin axis toward its rims (Sierra de Juárez and Sierra Cucapah). Our model temperatures agree within an error of 4% with the observed temperature profiles measured at boreholes. Our heat‐flow determinations agree within an error of ±15% with extrapolated observations. The numerical and chemical analyses support the hypothesis of fluid circulation between the clay–lutite layer and the fractured granitic basement. Thermal modelling shows low heat‐flow values along the Laguna Salada Basin. Deep fluid circulation patterns were observed that redistribute such flow at depth. Two patterns were distinguished. One displays the heat flow increasing from the basin axis towards its borders (temperature increase of 20°C). The second pattern shows an increasing heat flow from south to north of the basin. Such behaviour is confirmed by the temperature measurements in the thermometric boreholes.     

磷灰石裂变径迹分析在恢复盆地沉降抬升史中的应用──以鄂尔多斯盆地为例

对鄂尔多斯盆地磷灰石裂变径迹资料深入分析表明．最迟２３Ｍａ以来盆地发生了一期由于快速抬升剥蚀引起的冷却事件．盆地东部以９５ｍ／Ｍａ的速率抬升，造成约２０００ｍ的剥蚀量；而盆地西部则以５６ｍ／Ｍａ的速率抬升，导致了约１０００ｍ的剥蚀量．盆地东、西部的差异抬升剥蚀导致了盆地现今微微西倾的构造面貌．这一抬升剥蚀事件是印度板块与欧亚板块碰撞引起亚洲构造运动形式以挤压为主，转换为中新世以来以地壳增厚为主的结果．Ｋ－Ａｒ年龄和镜质体反射率资料分析表明，盆地在１７０－１６０Ｍａ（中侏罗末）曾发生一期热事件，使古地温梯度达５７℃／ｋｍ，古热流值达９６－１０９ｍｗ／ｍ^２．     

渤海唐山地区的地下热结构及地表热流分布的数值模拟

本文采用空间有限元和时间有限差分法模拟了渤中盆地区域的演化历史和地下热结构,并分别计算了渤中上地幔的热隆起与浅部的沉积层对渤中盆地内部温度场的影响。文中得到的计算结果较好地解释了渤中盆地地区热流正异常的观测事实,据此,对渤中盆地在形成过程中产生的热应力与地震孕震力源的关系也进行了一些探讨     

珠江口盆地大地热流特征及其与热岩石圈厚度的关系

沉积盆地现今热流特征是岩石圈构造-热演化过程的综合反映和盆地热史恢复的必要约束条件,其总体变化趋势与热岩石圈厚度密切相关.本文根据新收集的珠江口盆地19口钻井温度数据,新增计算了19个大地热流数据,其中12个数据位于深水区（水深大于300 m）,丰富了该盆地深水区钻井地热数据.结合前人研究成果,绘制了该盆地的大地热流图,并分析了其热流分布特征.在此基础上通过求解一维热传导方程,计算得到36口井位处的热岩石圈厚度,量化了盆地大地热流与热岩石圈厚度间的关系.结果显示,珠江口盆地大地热流值介于24.2~121.0 mW·m-2,平均71.8±13.6 mW·m-2,新增盆地深水区钻井平均热流值高达84.5±4.4 mW·m-2.大地热流分布整体上从陆架区到陆坡区升高,而热岩石圈厚度整体分布趋势与大地热流相反.大地热流与热岩石圈厚度间存在良好的指数相关性.     

张裂型盆地地热参数的垂向变化与琼东南盆地热流分布特征

海洋热流数据是开展海洋地球动力学研究和油气资源评价的基础数据.为深入认识琼东南盆地的地热特征,本文首先利用耦合沉积作用与岩石圈张裂过程的数值模型分析了张裂型盆地主要地热参数的垂向变化特征;并通过钻孔资料的详细分析,获得了琼东南盆地44口钻孔的热流数据;结合海底地热探针获取的热流数据,对琼东南盆地地热特征及其主要影响因素进行了简要分析.结果表明:沉积作用的热披覆效应对表层热流有较明显的抑制作用,由于沉积物生热效应与披覆效应的共同作用,同一钻孔处海底表层热流与钻孔深度3000~4000 m处热流或与海底间的平均热流差异很小,可以一起用于分析琼东南盆地的热流分布特征;莺歌海组、乐东组热导率随深度变化小于黄流组及其下地层热导率的变化,钻孔沉积层平均热导率约为1.7 W·(m·K)-1,钻孔地层生热率一般低于2.5 μW·m-3,平均生热率为1.34 μW·m-3,平均地温梯度主要介于30~45 ℃/km,热流介于50~99 mW·m-2,陆架区热流主要集中于60~70 mW·m-2,深水区钻孔具有较高的地温梯度和热流值;从北部陆架与上陆坡区往中央坳陷带,热流值从50~70 mW·m-2,增高为65~85 mW·m-2,并且往东有升高趋势,在盆地东部宝岛凹陷、长昌凹陷与西沙海槽北部斜坡带构成一条热流值高于85 mW·m-2的高热流带.进一步分析认为,琼东南盆地现今热流分布特征是深部热异常、强烈减薄岩石圈的裂后冷却作用、晚期岩浆热事件、地壳与沉积层的生热贡献以及沉积作用的热披覆效应等多种主要因素综合作用的结果.     

济阳坳陷沾化东区块现今地温场及热历史

根据区内钻井测温和镜质体反射率RO资料,在该区现今地温场分析的基础上,采用岩石圈和盆地尺度相结合的多期热演化模型对该区的热史进行了恢复.结果表明：(1)现今地温梯度为35.8℃/km.孤岛及垦东地区地温梯度较高,大于3℃/km.(2)古新世早期的大地热流值为83.6mW/m2,相当于现代活动裂谷的热流值.自古新世以来,盆地有逐步变冷的趋势,中间有二次回升,但回升的幅度逐渐变小.现今大地热流值为63mW/m2,接近全球大地热流的平均值.(3)该区主力烃源岩经历了持续的受热过程,现今仍处于“生油窗口”内,在深度上具有较大的油气赋存空间,热演化背景对油气生成有利.上述模拟结果可以用渤海湾盆地的裂谷演化模式进行较好的解释,并对胜利油田沾化凹陷的油气勘探工作具有一定的实际意义.     

Environmental isotopic and hydrochemical characteristics of groundwater systems in Daying and Qicun geothermal fields,Xinzhou Basin,Shanxi, China

The conceptual hydrogeological model of the low to medium temperature Daying and Qicun geothermal fields has been proposed, based on hydrochemical characteristics and isotopic compositions. The two geothermal fields are located in the Xinzhou basin of Shanxi, China and exhibit similarities in their broad‐scale flow patterns. Geothermal water is derived from the regional groundwater flow system of the basin and is characterized by Cl·SO₄‐Na type. Thermal water is hydrochemically distinct from cold groundwater having higher total dissolved solids (TDS) (>0·8 g/l) and Sr contents, but relatively low Ca, Mg and HCO₃ contents. Most shallow groundwater belongs to local flow systems which are subject to evaporation and mixing with irrigation returns. The groundwater residence times estimated by tritium and ¹⁴C activities indicate that deep non‐thermal groundwater (130–160 m) in the Daying region range from modern (post‐1950s) in the piedmont area to more than 9·4 ka BP (Before Present) in the downriver area and imply that this water belong to an intermediate flow system. Thermal water in the two geothermal fields contains no detectable active ¹⁴C, indicating long residence times (>50 ka), consistent with this water being part of a large regional flow system. The mean recharge elevation estimated by using the obtained relationship Altitude (m) = ? 23·8 × δ²H (‰ ) ? 121·3, is 1980 and 1880 m for the Daying and Qicun geothermal fields, respectively. The annual infiltration rates in the Daying and Qicun geothermal fields can be estimated to be 9029 × 10³ and 4107 × 10³ m³/a, respectively. The variable ⁸⁶Sr/⁸⁷Sr values in the thermal and non‐thermal groundwater in the two fields reflect different lithologies encountered along the flow path(s) and possibly different extents of water‐rock interaction. Based on the analysis of groundwater flow systems in the two geothermal fields, hydrogeochemical inverse modelling was performed to indicate the possible water‐rock interaction processes that occur under different scenarios. Copyright © 2010 John Wiley & Sons, Ltd.     

济阳坳陷沾化东区块现今地温场及热历史

根据区内钻井测温和镜质体反射率RO资料，在该区现今地温场分析的基础上，采用岩石圈和盆地尺度相结合的多期热演化模型对该区的热史进行了恢复.结果表明：(1)现今地温梯度为35.8℃/km.孤岛及垦东地区地温梯度较高，大于3℃/km.(2)古新世早期的大地热流值为83.6mW/m²，相当于现代活动裂谷的热流值.自古新世以来，盆地有逐步变冷的趋势，中间有二次回升，但回升的幅度逐渐变小.现今大地热流值为63mW/m²，接近全球大地热流的平均值.(3)该区主力烃源岩经历了持续的受热过程，现今仍处于“生油窗口”内，在深度上具有较大的油气赋存空间，热演化背景对油气生成有利.上述模拟结果可以用渤海湾盆地的裂谷演化模式进行较好的解释，并对胜利油田沾化凹陷的油气勘探工作具有一定的实际意义.     

CHANGES IN MORPHOLOGY AT THE MARGIN OF THE GREENLAND ICE SHEET (LEVERETT GLACIER), IN THE PERIOD 1943–1992: A QUANTITATIVE ANALYSIS

Changes in ice-marginal morphology near Leverett glacier, a small outlet glacier at the western margin of the Greenland ice sheet, are determined from a photogrammetrical analysis. To be able to compare two datasets from subsequent years with measurements at different coordinates, kriging was used for interpolation. In this study the kriging standard error is used to evaluate the relative accuracy of the resulting maps. Aerial photographs of 1943, 1968 and 1985 were compared. In the period 1943–1968 an area of 0.2 × 10⁶ m² was deglaciated. Approximately 1.1 × 10 m³ of material is deposited in this area. The southern part of the deglaciated area is characterized by ice-cored moraines, while moraines without ice core were formed in the north. Differences in depositional products reflect differences in meltwater activity and probably ice-marginal thermal regime. During deglaciation a small proglacial sandur decreased in altitude by 3.2 ± 0.1 m. From the early 1970s Leverett glacier advanced over a previously deglaciated area. During this advance, small ice-marginal accumulations were incorporated and eroded by the advancing glacier. Erosion products were for a substantial part stored in the proglacial sandur. About 1.2 ×10⁵ m² of the northern part of an ice-cored moraine complex decreased in altitude by −3.6 ± 0.1 m from 1943 to 1968 and over 2.7 × 10⁴ m² by −2.7 ± 0.1 m during 1968–1985. The spatial patterns of altitude change were analysed in relation to topomorphological parameters as exposition and slope angle and areas occupied by lakes. The estimated energy used to melt the subsurface ice of the ice-cored moraine is 1.4–2.2 W m² (1943–1968) and 1.0–1.6 W m² (1968–1985). These values are 30–50 times larger than the geothermal heat flux. For the expected average debris concentration of the ice core (< 10 per cent by volume) the deviation of the surface energy balance forced by climate change will be small and encompass an insignificant part of the total estimated energy used for melting.     

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.     

Thermal structure and energy of the Hakone volcano,Japan

The thermal energy balance and the temperature profile of the Hakone volcano are considered quantitatively. Across the Hakone volcano and its surroundings the heat flow values vary from 10^–1 to 10³ mW/m², due to thermal conduction and mass flow involving volcanic steam and hot spring discharge. An area with extremely low heat flow is observed in the western side of the caldera showing the presence of percolating meteoric water. The hydrothermal activity is intense in the eastern half of the caldera.The total heat discharge from the high temperature zone (discharge area) of the Hakone volcano amounts to 11.0×10⁷ W. The magmatic steam energy discharge is 95.0×10⁶ W. The thermal energy by redistribution of the terrestrial heat flow by the lateral deep ground water flow is calculated to be 9.00×10⁶ W. For the model having the vertical vent in the volcano's central part up to 1 km depth below the ground surface from a magma reservoir the computed temperature distribution is consistent with the observed values. The depth of the magma reservoir is 7 km below the ground surface and the diameter is 5 km.     

Geoelectromagnetic measurements across the southern Senegal basin (West Africa)

Magnetotelluric and differential geomagnetic sounding surveys, consisting of nine soundings, were performed in 1984 along a 200-km profile across the southern Senegal basin. They were intended to obtain information concerning the resistivity structure of the crust and upper mantle and the distribution of the induced electric currents. Magnetotelluric data indicate that two-dimensional resistivity models are appropriate for the region. The zone above the basement is highly inhomogeneous in geoelectrical structure. Very conductive sediments (0.6-3 ohm m) appear in the Mesozoic-Cenozoic Senegal basin. These sediments lie at depths of up to 4500 m on the west end of the profile. Below this, a modest resistivity material (10–30 ohm m) extends to a maximum depth of about 3000 m. The material at depth on the cast part of the traverse line is thought to be Palaeozoic sediments of the Bove basin. The depth of the magnetotelluric basement lies between about 250 m (in the east) and 4800 m (in the west). The crust is characterised by a drop in electrical resistivity at a depth of 15 km below the east part of the profile. Considering the total section, we observe a general trend towards lower resistivities at depths in excess of 100 km, the transition from 2000 ohm m to about 2 ohm m occurs in the depth range 100 to 175 km. An analysis of the geomagnetic variation field has identified a concentration of telluric current flow beneath the deep basin. It appears that the additional currents flowing in the striking direction of the Senegal basin are largely controlled by sedimentary rocks of high conductivity lying at depths less than 5 km. Model studies show that the local conductivity distribution is able to explain the currents circulating in the thick well-conducting sediments.     

辽河裂谷盆地地幔热流

辽河盆地是一个在前中生代基底上发展起来的裂谷盆地,实测大地热流平均值为65mW/m²,变动于44-83mW/m²之间。在给出地壳结构模型并确定各岩层放射性生热率的基础上,采用“剥层”法从地表开始,自上而下,由浅及深地扣除各岩层所提供的热量,从而得出地幔热流值。结果表明,辽河裂谷盆地地幔热流为41mW/m²,占整个地表总热流量的63％。可见,本区热量大部分来自地幔。与世界上其它地质构造单元相比,辽河裂谷盆地无论地幔热流绝对值或其与地表热流之比值,都具有介于稳定地区和构造活动区之间的特点。作者认为,辽河裂谷盆地地幔热流的上述特点,乃是中、新生代以来本区长期地质历史发展的产物。     

Deep circulation driven by strong vertical mixing in the Timor Basin

The importance of deep mixing in driving the deep part of the overturning circulation has been a long debated question at the global scale. Our observations provide an illustration of this process at the Timor Basin scale of ～1000 km. Long-term averaged moored velocity data at the Timor western sill suggest that a deep circulation is present in the Timor Basin. An inflow transport of ～0.15 Sv is observed between 1600 m and the bottom at 1890 m. Since the basin is closed on its eastern side below 1250 m depth, a return flow must be generated above 1600 m with a ～0.15 Sv outflow. The vertical turbulent diffusivity is inferred from a heat and transport balance at the basin scale and from Thorpe scale analysis. Basin averaged vertical diffusivity is as large as 1 × 10^?3 m² s^?1. Observations are compared with regional low-resolution numerical simulations, and the deep observed circulation is only recovered when a strong vertical diffusivity resulting from the parameterization of internal tidal mixing is considered. Furthermore, the deep vertical mixing appears to be strongly dependent on the choice of the internal tide mixing parameterization and also on the prescribed value of the mixing efficiency.     

塔里木盆地航磁场分析与磁源体结构

本文研究塔里木盆地区域磁异常图的反演及磁源体结构.由于众多异常的叠加和反演固有的多解性,区域磁异常图的准确解释是非常困难的.三维欧拉反褶积是一种确定地质体位置和埋藏深度的自动定量反演方法,比较适用于计算区域磁异常源的埋藏深度.由于大型克拉通沉积盆地地层具有上新下老的规律性,将磁异常源分解为三个深度层次,圈定它们各自的分布区域,便可将它们与形成的地质作用及时代联系起来,为准确解释区域磁异常图提供可靠的依据.本文应用三维欧拉反褶积反演方法,计算出的塔里木盆地深度为2~5 km、5~10 km、10~20 km三个等级的磁异常源,它们与形成的地质作用及时代分别为: 中生代构造运动,海西期玄武岩侵位和太古代结晶基底的变质作用;圈定了它们各自的分布区域.