地球物理学报  2016, Vol. 59 Issue (3): 1003-1016   PDF    
引用本文
常健, 邱楠生, 赵贤正等. 2016. 渤海湾盆地冀中坳陷现今地热特征. 地球物理学报, 59(3): 1003-1016, doi: 10.6038/cjg20160322   
Chang J, Qiu N S, Zhao X Z, et al. 2016. Present-day geothermal regime of the Jizhong depression in Bohai Bay basin, East China. Chinese J. Geophys. (in Chinese), 59(3): 1003-1016, doi: 10.6038/cjg20160322   
渤海湾盆地冀中坳陷现今地热特征
常健1,2, 邱楠生1,2, 赵贤正3, 许威1,2, 徐秋晨1,2, 金凤鸣3, 韩春元3, 马学峰3, 董雄英3, 梁小娟3    
1. 中国石油大学(北京)油气资源与探测国家重点实验室, 北京 102249;
2. 中国石油大学(北京)盆地与油藏研究中心, 北京 102249;
3. 中国石油华北油田公司, 河北 任丘 062552
收稿日期: 2015-03-31, 2016-01-27 收修定稿
基金项目: 国家自然科学基金项目(41125010,91114202,41402112)和中国石油大学(北京)科研基金(2462013YJRC034)资助.
作者简介: 常健,男,1982年生,河北人,讲师,主要从事构造地质学和低温热年代学方面的研究.E-mail:changjian@cup.edu.cn
摘要: 渤海湾盆地冀中坳陷是我国最典型的潜山油气藏富集区.本文借助117口钻井地层测温资料和45块实测岩石热导率数据系统研究了冀中坳陷现今地温梯度、大地热流、热岩石圈厚度、岩石圈热结构等地热特征参数.研究表明,冀中坳陷0~3000 m统一深度现今地温梯度为20.8~41.0℃·km-1,平均值为31.6℃·km-1,比未校正值减小1~3℃·km-1;现今大地热流介于48.7~79.7 mW·m-2,平均值为59.2 mW·m-2.平面上,冀中坳陷现今地温梯度和热流由西向东(从盆地边缘向内部)逐渐增大,并且凸起区地温梯度和热流相对较高,而凹陷区则偏低,与基底地形起伏具有很好的对应关系.同时,冀中坳陷腹部高热流凸起区广泛分布地热田.冀中坳陷现今热岩石圈厚度为98~109 km,其岩石圈热结构为一典型的"冷壳热幔"型.本研究不仅对冀中坳陷油气勘探与地热能开发具有重要的指导意义,而且为深部岩石圈研究(华北克拉通破坏科学问题)提供了新依据.
关键词: 冀中坳陷     地温梯度     大地热流     岩石圈热结构     渤海湾盆地    
Present-day geothermal regime of the Jizhong depression in Bohai Bay basin, East China
CHANG Jian1,2, QIU Nan-Sheng1,2, ZHAO Xian-Zheng3, XU Wei1,2, XU Qiu-Chen1,2, JIN Feng-Ming3, HAN Chun-Yuan3, MA Xue-Feng3, DONG Xiong-Ying3, LIANG Xiao-Juan3     
1. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum, Beijing 102249, China;
2. Basin and Reservoir Research Center, China University of Petroleum, Beijing 102249, China;
3. PetroChina Huabei Oilfield Company, Hebei Renqiu 062552, China
Abstract: As a second major tectonic unit of the Bohai Bay basin,the Jizhong depression is a very typical area in China with enrichment of buried-hill hydrocarbon reservoirs.Based on 117 borehole temperature data and 45 measured thermal conductivity data,we studied present-day geothermal characteristics of the Jizhong depression systematically such as geothermal gradient,terrestrial heat flow,thermal lithosphere thickness and lithospheric thermal structure. The results show that the present-day geothermal gradients at depths 0~3000 m in this region range from 20.8℃·km-1 to 41.0℃·km-1 with a mean of 31.6℃·km-1, lower 1~3℃·km-1 than uncorrected value. The heat flow varies between 48.7~79.7 mW·m-2 with an average of 59.2 mW·m-2. Present-day geothermal gradients and heat flow of this depression gradually increase from west(basin margin) to east(basin interior) horizontally. These values in sub-uplift areas of the Jizhong depression are usually higher than those of the sub-depression regions,and their variations correspond to the basement relief very well. Meanwhile,the uplift areas with high heat flow in the central Jizhong depression host amounts of geothermal fields. In addition,we calculated the thickness of present-day thermal lithosphere of study area to be 98~109 km,and its lithospheric thermal structure is a "cold crust and hot mantle" type now. This study has an important guiding significance for hydrocarbon exploration and geothermal energy exploitation in the Jizhong depression,and provides new insights for deep lithosphere research(especially for the scientific problem-North China Craton destruction).
Key words: Jizhong depression     Geothermal gradient     Heat flow     Lithospheric thermal structure     Bohai Bay basin    
1 引言

随着1975年任丘古潜山油田的发现，冀中坳陷成为我国最重要的潜山油气勘探区. 冀中坳陷主要发育Es1，Es3，Es4-Ek和C-P四套烃源岩，而中元古界长城系的高于庄组、蓟县系的雾迷山组、下古生界寒武系的府君山组和奥陶系是主要的油气储集层(杜金虎等，2002; 赵利杰等，2012; 尘福艳等，2013). 近年来，华北油田加大了深潜山和潜山内幕油气藏的勘探力度，陆续发现了牛东1、长3、宁古8等多个隐蔽性高产潜山油气藏(赵贤正等，2014; Zhao et al.，2015)，这表明冀中坳陷还具有良好的勘探前景. 温度是制约烃源岩成熟演化的最关键因素，因此正确认识一个地区的现今地热特征对油气成熟演化研究具有重要意义. 一些学者曾对冀中坳陷现今地热特征进行过研究(陈墨香等，1982; 冯石和郝石生，1983; 周瑞良，1987; 龚育龄，2003; 林世辉和龚育龄，2005). 陈墨香等(1982)根据钻孔实测温度曾研究过冀中坳陷牛驼镇凸起地区现今地温梯度在纵向和平面上的分布特征，并认为岩石热导率垂向-侧向变化是造成这一地区地热异常形成的关键因素; 而周瑞良(1987)通过对华北平原北部深层地下水活动的细致研究，认为水动力条件是造成冀中坳陷不同地区地温场差异分布的主要因素. 结合系统测温数据和岩石热导率数据，冯石和郝石生(1983)曾探讨过冀中坳陷饶阳凹陷和廊固凹陷现今热流分布特征和新生界地层古地温演化史; 而龚育龄(2003)及林世辉和龚育龄(2005)对整个冀中坳陷现今地温梯度或大地热流分布特征进行过研究，认为冀中坳陷现今地温梯度为25.8~48.8 ℃·km-1，大地热流为43.9~90.0 mW·m-2，并认为其平面分布与基底地形起伏有关.

本文结合近10年来新增的一大批测温数据和实测岩石热导率数据，不仅重新剖析了冀中坳陷现今地温梯度和大地热流分布特征和影响因素，而且研究了冀中坳陷0～3000 m统一深度的地温梯度平面分布特征，从而更真实地反映冀中坳陷现今热状态.同时，根据生热率、热导率及分层数据，我们分析了冀中坳陷现今热岩石圈厚度和岩石圈热结构.本研究对冀中坳陷油气和地热能勘探与开发具有重要的指导意义.

2 地质背景

冀中坳陷是发育在华北克拉通之上的一个中、新生代断陷-坳陷盆地. 作为渤海湾盆地的一个次级构造单元，它北依燕山隆起，南抵邢衡隆起，西临太行山隆起，东至沧县隆起，呈北东-南西走向，面积约3.2×104 km2(图 1). 区域上，冀中坳陷被两条近东西向-北西西向变换带(即无极-衡水变换带和徐水-安新-文安变换带)又分割成南、中、北等三个区，其中，南区包括石家庄凹陷、晋县凹陷、束鹿凹陷、无极-藁城低凸起和宁晋凸起; 中区分布有饶阳凹陷、保定凹陷、深县凹陷、高阳凸起、深泽低凸起及蠡县斜坡; 北区包括北京凹陷、大厂凹陷、廊固凹陷、徐水凹陷、武清凹陷、霸县凹陷、大兴凸起、牛驼镇凸起、容城凸起、牛北斜坡、杨村斜坡及文安斜坡. 另外，研究区还包括沧县隆起的大城凸起和兴起凸起(图 1).

图 1 (a)渤海湾盆地构造单元图(据徐守余和严科，2005): ①辽河坳陷; ②渤中坳陷; ③黄骅坳陷; ④沧县隆起; ⑤冀中坳陷; ⑥埕宁隆起; ⑦济阳坳陷; ⑧临清坳陷; ⑨内黄隆起; ⑩东濮坳陷.(b)冀中坳陷构造分区及测温井和地热田分布(修改自杜金虎等(2002)，赵贤正等(2010)，梁宏斌等(2010)) Fig. 1 (a)Tectonic units of Bohai Bay basin(Xu and Yan，2005): ①Liaohe depression; ② Bozhong depression; ③Huanghua depression; ④Cangxian uplift; ⑤Jizhong depression; ⑥ Chengning uplift; ⑦Jiyang depression; ⑧ Linqing depression; ⑨ Neihuang uplift; ⑩Dongpu depression.(b)Tectonic units of Jizhong depression and temperature-measurement wells and geothermal fields(modified from Du et al.(2002)，Zhao et al.(2010)，Liang et al.(2010))

冀中坳陷基底主要发育一套太古界及下元古界变质岩. 中、晚元古代-中奥陶世，受海侵作用影响，冀中坳陷地区主要发育浅海相碳酸盐岩，沉积厚度达~5000 m(赵贤正等，2010; 陆诗阔等，2011). 晚奥陶世-晚石炭世时期的加里东、海西运动致使冀中坳陷整体抬升而遭受长期剥蚀(杜金虎等，2002; 陆诗阔等，2011). 二叠纪时期，冀中坳陷沉积了一套海陆交互相碎屑岩夹碳酸盐岩，厚约700~800 m(杜金虎等，2002). 早中生代的强烈挤压作用(印支运动引起的)造成冀中坳陷在三叠纪-侏罗纪形成一系列褶皱，并遭受了强烈剥蚀，仅在其南部沉积有少量这一时期地层(杨明慧等，2001; 吴智平等，2007). 晚白垩世，冀中坳陷地区的构造应力场由挤压作用转变为伸展作用，岩浆活动强烈，而发育一套陆源碎屑岩夹煤层和火山岩地层(孙冬胜等，2004a，2004b). 始新世晚期，受太平洋板块俯冲诱发深部地幔物质对流的影响，整个华北地区处于弧后伸展状态，而冀中坳陷也进入强烈伸展断陷演化阶段，受断块掀斜的控制，发育大量NE-NNE及NW向潜山构造，这一时期主要沉积陆相湖盆砂岩和泥岩，夹部分泥灰岩和石膏(杨明慧等，2002; 崔周旗，2005; 张文朝等，2008).新近纪，华北地区的弧后拉张作用减弱，冀中坳陷地区进入坳陷演化阶段，以整体沉降为主，主要发育河流相的砂泥岩(梁苏娟，2001; 董大伟等，2013).

3 热导率测试与分析

本研究共选取了冀中坳陷苏50、西59、馆深1、高深1、葛2、大参1、京101、丘3、王9、磁1、榆31、西55、胜1等13口井45块样品进行热导率测试，层位涉及到新生界、中生界、石炭-二叠系、下古生界、中-上元古界和太古宙; 岩性主要有砂岩、泥岩、火山岩、白云岩、灰岩等. 热导率测试工作是在中国科学院地质与地球物理研究所进行的，使用仪器是德国生产的TCS(Thermal Conductivity Scanning)热导率自动扫描仪，其测量范围为0.2~25 W/(m·K)，精度为±3%.

图 2a显示冀中坳陷45块样品的岩石热导率集中于1.5~3.5 W/(m·K)，主要对应于中-新生界碎屑岩. 图 2b表明冀中坳陷热导率与样品深度无明显相关性. 对于浅部碎屑岩，热导率分布散乱; 而不同深度的碳酸盐岩样品具有相似的热导率. 结合龚育龄(2003)已发表的热导率数据，我们构建了冀中坳陷热导率柱(表 1)，为研究冀中坳陷现今地温梯度和大地热流提供了可能.

图 2 冀中坳陷实测岩石热导率直方图(a)及不同岩性热导率-深度关系图(b) Fig. 2 (a) Histogram of measured thermal conductivity in Jizhong depression and (b) relationship between thermal conductivity and depths by different lithology

表 1 冀中坳陷岩石热导率柱 Table 1 Thermal conductivity data of different formations in Jizhong depression
4 冀中坳陷现今地温梯度分布特征

沉积盆地现今地温梯度主要借助各种钻孔测试温度数据(包括系统稳态测温、试油温度、准稳态测温、瞬态测温等四类)进行研究.其中，系统稳态测温资料和试油温度数据最可靠.系统稳态测温指钻孔热平衡时间以后的测温，一般每隔10 m或20 m测量一个温度点.它反映了研究区最真实的地温状况，测温资料最可靠，精度最高.本次研究在华北油田勘探开发研究院收集到了冀中坳陷28口井的系统稳态测温数据，通过温度-深度线性拟合得到了每口井的现今地温梯度(图 3a-3e; 表 2).其中牛驼镇凸起地区的容1、雄101、雄104、淀6、淀6-3、霸8、牛8-2等钻井的现今地温梯度介于48.1~72.8 ℃·km-1，属于高异常值.周瑞良(1987)研究认为这与区域水循环系统有关.由于这些钻井地区受水动力条件影响严重，且热传递方式以热对流为主，不能代表真实的热状态，在研究冀中坳陷地温梯度分布规律时应忽略掉.其余钻孔的地温梯度介于27.4~39.7 ℃·km-1之间，平均值为34.7 ℃·km-1.

图 3 冀中坳陷钻孔温度-深度关系图
(a—e) 系统测温曲线; (f) 试油温度数据. Fig. 3 Relationship between borehole temperatures and depths in Jizhong depression
(a—e) Steady-state temperature curves of boreholes; (f) Oil-testing temperature data.

表 2 冀中坳陷28口系统连续测温钻孔地温梯度和大地热流计算结果 Table 2 Calculated geothermal gradients and heat flow values of 28 systematically-continuously measured temperatures of boreholes in Jizhong depression

试油温度是在开发井或生产井中的含油层段进行试油测压时获得的低温数据，它也能够较真实地反映盆地的温度状况.本文除收集到林世辉和龚育龄(2005)已发表的400个试油温度数据外，还搜集了华北油田近10年来新钻的200多口探井的试油温度数据(图 3f)，测温深度跨度为1000~6000 m，而地层温度介于40~200 ℃之间.另外，利用公式(1)计算了其中89口代表性钻孔的地温梯度(图 1; 表 3).

其中，G为地温梯度(℃·km-1); Ts为试油温度; T0为近地表恒温带温度(设为14 ℃); Z为测温深度(km); Z0为恒温带深度(设为20 m).由于许多井在不同深度上有几个温度值，因此在同一口井可计算出多个梯度值，为此还需将同一口井不同深度上的多个地温梯度值加以算术平均，从而得到一口井的地温梯度值.由试油温度计算得到的地温梯度介于28.5~43.6 ℃·km-1(不考虑牛28井，其地温梯度为高异常值55.0 ℃·km-1)，平均值为33.2 ℃·km-1.综合两种测温数据得到的地温梯度，冀中坳陷现今地温梯度为27.4~43.6 ℃·km-1，平均值为33.5 ℃·km-1.
表 3 冀中坳陷89口典型试油温度钻孔地温梯度和大地热流计算结果 Table 3 Calculated geothermal gradients and heat flow values of 89 typical oil-testing boreholes in Jizhong depression

图 4是计算获得的117口钻孔地温梯度平面分布图.从平面上整体看，冀中坳陷现今地温梯度分布与下伏基底构造基本一致.由西到东，现今地温梯度总体变化趋势为低-高-低-高，与凹-凸-凹-凸的地貌特征具有很好的对应关系.冀中坳陷西部(太行山东麓)北京-徐水-保定-石家庄凹陷一带，现今地温梯度介于28~32 ℃·km-1; 其东部紧邻的牛驼镇-容城-高阳凸起区地温梯度升高，达36~38 ℃·km-1; 而牛驼镇-容城-高阳凸起区东部的廊固、霸县和饶阳三大含油气坳陷区现今地温梯度又降至32~34 ℃·km-1; 之后向东地温梯度又增大，到沧县隆起区已增大到38~40 ℃·km-1.

图 4 冀中坳陷现今地温梯度分布平面图 Fig. 4 Distribution of present-day geothermal gradients in Jizhong depression

另外，沉积盆地现今地温梯度还与深度有关，一般浅部地温梯度大于深部.因此不同钻孔测温深度的差异会导致直接利用钻井测温数据计算的地温梯度在空间上不具有可对比性，应当统一深度范围之后再分析地温梯度平面分布规律.对某一口钻井，其统一深度0~3000 m地温梯度计算方法如下: 首先，结合实测岩石热导率和地层分层数据，利用一维热稳态传导方程(5)计算获得3000 m深度处的温度; 然后，用计算获得的温度减去地表温度后除以深度值(3 km)即为0~3000 m地温梯度.据此，我们计算得到了冀中坳陷117口代表性钻孔统一深度0~3000 m的地温梯度(表 2和3).校正后的地温梯度分布范围为20.8~41.0 ℃·km-1，平均值为31.6 ℃·km-1.大部分系统测温钻孔校正后的地温梯度比未校正前普遍减小1~3 ℃·km-1.图 5显示冀中坳陷校正后地温梯度平面分布规律与基底构造起伏呈现出更好的对应关系.与校正前相比，冀中坳陷现今地温梯度降低了~2 ℃·km-1.冀中坳陷西缘地温梯度为26~30 ℃·km-1; 牛驼镇-容城-高阳凸起一带地温梯度增大为32~36 ℃·km-1; 廊固、霸县和饶阳地区的地温梯度为26~30 ℃·km-1，而东部沧县隆起一带地温梯度增至34~40 ℃·km-1.

图 5 冀中坳陷0~3000 m统一深度现今地温梯度平面分布图 Fig. 5 Distribution of present-day geothermal gradients in Jizhong depression at depths 0~3000 m
5 冀中坳陷现今大地热流分布特征

在大地热流计算过程中，针对不同类型的测温数据采用了不同的计算方法.对于系统测温数据，大地热流计算公式如下:

其中，qs为大地热流(mW·m-2); Kt为调和岩石热导率(W/(m·K)); dT/dZ为地温梯度(℃·km-1).

而对于试油温度，采用的是Chapman等(1984)提出的热阻法，其计算公式为:

式中，qs同上; ΔZi为相邻构造层间的深度差; Ki为Zi深度处地层热导率(W/(m·K)); Ts为试油温度; T0为近地表恒温带温度(14 ℃).

根据这两种计算方法，获得了冀中坳陷117口钻井大地热流值(表 2和3).冀中坳陷大地热流值介于48.7~79.7 mW·m-2之间(未考虑容1、雄101、雄104、淀6、淀6-3、霸8、牛8-2、牛28井)，平均值为59.2 mW·m-2.在此基础上，绘制了冀中坳陷大地热流平面分布图(图 6).冀中坳陷西部(即太行山东麓)大地热流偏低，仅为54~57 mW·m-2，向东逐渐增大，到沧县隆起区，大地热流增大到72~75 mW·m-2.在冀中坳陷内部，廊固、霸县、饶阳三大凹陷区及周边斜坡带大地热流为51~57 mW·m-2; 而凸起区(牛驼镇、容城、高阳低凸起)大地热流偏高，为60~69 mW·m-2，容城凸起南部高达72 mW·m-2以上.与地温梯度分布规律相似，冀中坳陷大地热流展布也与基底构造起伏具有很好的对应关系.另外，冀中坳陷腹部高热流区通常分布有地热田(图 6).

图 6 冀中坳陷现今大地热流平面分布图 Fig. 6 Distribution of present-day heat flow values in Jizhong depression
6 冀中坳陷岩石圈热结构和热岩石圈厚度

沉积盆地大地热流主要由两部分组成: 一个是地壳热流，其热量来自地壳浅部U、Th、K等放射性元素的衰变; 另一个是地幔热流，其热量来自深部地幔.岩石圈热结构就是指一个地区壳、幔两部分热流的配分比例及其组构关系(Blackwell，1971).理论上，地幔热流可以通过实测热流与岩石生热率之间的线性关系获取(Birch et al.，1968).但在沉积盆地区这种线性关系不一定存在.汪缉安等(1992)提出了一种基于地表热流、生热率和地壳分层厚度来计算沉积盆地区地幔热流的“回剥”方法，其计算公式如下:

式中，qc，qm和qs分别为地壳、地幔和地表热流(mW·m-2); Ai为第i层的放射性生热率(μW·m-3); Di为第i层的厚度(km).在计算冀中坳陷现今岩石圈热结构时，构造层分层数据参考的是龚育龄(2003)和徐杰等(2001); 浅部地层热导率由本次实测值获得，而深部热导率参考龚育龄(2003); 生热率数据来自邱楠生等(2015)，详细数据见表 4.经计算，冀中坳陷地幔热流为32.43 mW·m-2(图 7)，从而得到地壳热流为26.77 mW·m-2.Wang(1996)根据地幔和地壳热流比例曾提出了“冷壳热幔”和“热壳冷幔”的热结构类型.本次研究表明冀中坳陷现今地壳热流小于地幔热流，比值为0.83，为一典型的“冷壳热幔”型热结构.

表 4 冀中坳陷主要构造层厚度、岩石生热率及热导率 Table 4 Thicknesses,heat generations and conductivity of main tectonic layers in Jizhong depression

图 7 冀中坳陷各构造层热流计算结果
qs为冀中坳陷平均大地热流. Fig. 7 Heat flow calculation results of tectonic layers in Jizhong depression
qs is mean heat flow of Jizhong depression.
热岩石圈是指以热传导方式进行热传递的岩石圈层，其下为以热对流为主的软流圈(Morgan，1984).利用一维稳态热传导方程(5)得到冀中坳陷岩石圈温度随深度变化的曲线，它们与地幔绝热线或干玄武岩固相线的交点所处的深度即为“热”岩石圈厚度.

式中，i为构造层数; qi上为第i构造层上界面的热流值(图 7); Zi为第i层的厚度(km)(表 4); Ti上、Ti下分别为第i构造层上、下界面的温度(℃); Ai为第i层的岩石生热率(μW/m3)(表 4);Ki为第i层岩石热导率(W/(m·K))(表 4).臧绍先等(2002)将热传导地温线与地幔绝热曲线相交的深度定义为岩石圈底界面，同时采用两条绝热线分别作为热岩石圈底面温度上限T1和下限T2的判据，其计算公式如下:

本文采用这种方法计算获得，冀中坳陷现今热岩石圈厚度平均值为98~109 km(图 8)，比地震岩石圈(厚度为110~120 km，Zhu等(2011))薄10 km左右.

图 8 冀中坳陷热岩石圈厚度计算结果 Fig. 8 Calculated thermal lithospheric thickness in the Jizhong depression
7 讨论与结论

我们通过对钻孔测温数据的校正计算得到，渤海湾盆地西部冀中坳陷0～3000 m统一深度现今地温梯度为20.8~41.0 ℃·km-1，平均值为31.6 ℃·km-1，比未校正值减小~2 ℃·km-1.邱楠生等(2009)在研究渤海海域现今地温梯度分布特征也曾利用相同的方法并得到相似的结果.这说明研究沉积盆地现今地温梯度时进行统一深度温度校正是必要的.冀中坳陷现今大地热流介于48.7~79.7 mW·m-2，平均值为59.2 mW·m-2.另外，前人还曾对渤海湾盆地其他坳陷现今地热特征进行过研究.Guo等(2011)通过线性拟合得到黄骅坳陷板桥凹陷地区的现今地温梯度为~32.6 ℃·km-1; 龚育龄等(2003)计算得到济阳坳陷的现今地温梯度为~35.5 ℃·km-1，现今大地热流为65.8 mW·m-2.区域上，渤海湾盆地现今地温梯度和大地热流具有从西向东逐渐增大的趋势.这除了与沉积层生热率有关外，还可能与华北克拉通岩石圈厚度变化有关.Zhu等(2011)通过地震资料得到，渤海湾盆地从西到东岩石圈厚度从~120 km 减小到80 km.渤海湾盆地东西向现今地热特征与岩石圈厚度具有明显的负相关性.然而，渤海湾盆地最东部渤海海域(岩石圈厚度仅65~80 km，Zhu等(2011))现今地温梯度仅为30.8 ℃·km-1(邱楠生等，2009)，比其他二级构造单元都小.这可能与后期快速沉降作用造成这一地区现今还未达到热平衡有关(邱楠生等，2007; 李文正等，2014).与冀中坳陷南部相邻的临清坳陷现今地温梯度为~31.1 ℃·km-1(龚育龄，2003);同时，Zuo等(2014)研究认为渤海湾盆地最南端东濮坳陷现今地温梯度为~32 ℃·km-1.它们与冀中坳陷现今地温梯度相近，这表明渤海湾盆地现今地温梯度并未受到南北向区域构造的影响.另外，渤海湾盆地冀中坳陷的现今热岩石圈厚度为98~109 km，而其西部的鄂尔多斯盆地热岩石圈厚度为130~140 km(汪洋和程素华，2011)，其东部的济阳坳陷热岩石圈厚度仅为71~90 km(刘绍文等，2005; Qiu et al.，2014).这表明，与地震岩石圈相似，华北克拉通热岩石圈也具有西厚东薄的特征(何丽娟，2014; He，2015)，从而为华北克拉通破坏研究提供了新依据.

综合前人研究成果认为，渤海湾盆地现今地温梯度为30~36 ℃·km-1，大地热流为58~69 mW·m-2，与西部塔里木盆地和准噶尔盆地相比(冯昌格等，2009; 饶松等，2013)，属于高地温梯度和高热流区，这主要是与区域构造演化背景有关.早期研究证实古近纪时期，渤海湾盆地处于裂陷发育期(漆家福等，1995; 侯贵廷等，2001; 丁增勇等，2008)，此时地幔上拱和岩浆活动频繁，不仅造成深部岩石圈减薄，而且导致这一地区地表热流增大到80 mW·m-2以上(左银辉，2011).新近纪以来，随着太平洋板块向欧亚大陆俯冲带后退，区域应力场由拉张转变为挤压，渤海湾盆地进入坳陷发育阶段.同时伴随着岩石圈厚度增厚和火山活动减弱，渤海湾盆地从新近纪至今一直处于热衰减状态.但由于冷却时间相对短，这一地区现今仍具有较高的地温梯度和大地热流(王良书等，2002; 龚育龄等，2003; 邱楠生等，2009).

冀中坳陷现今地温梯度和大地热流平面分布与基底地形起伏具有很好的对应关系.凸起区地温梯度和大地热流偏大，而凹陷区偏小.这主要是因为凹陷区沉积了相对较厚的沉积物，多以砂、泥岩为主，这些岩性热导率通常较低，而凸起区基岩热导率相对较大，沉积盖层和基岩热导率的差异产生了“热折射”效应和热流的再分配，于是造成热流向隆起部位聚集，结果坳陷部位形成低温区，凸起部分形成高温区(熊亮萍和高维安，1982; Lin et al.，2000; 冯昌格等，2009; 饶松等，2013).图 6表明冀中坳陷腹部凸起区广泛分布的地热田与高热流具有很好的对应关系.这除了与凸起区高热导率有关外，还与水动力条件有关(周瑞良，1987).在凸起区，深部地下水顺着开启性较好的断裂溶洞等通道垂向流动到浅部新生界地层时，会同时把深部热量带到浅部，从而形成高温异常区，即地热田(周瑞良，1987).

致谢 感谢中国科学院地质与地球物理研究所胡圣标研究员和博士生高鹏在岩石热导率测试方面提供的帮助; 感谢华北油田研究院提供的大量基础地质数据和测温数据！
参考文献
[1] Birch F, Roy R F, Decker E R. 1968. Heat flow and thermal history in New England and New York.//Zen E, White W S, Hadley J B, et al. Studies of Appalachian Geology:Northern and Maritime. New York:Interscience, 437-451.
[2] Blackwell D D. 1971. The thermal structure of the continental crust.//Heaoock J G. The Structure and Physical Properties of the Earth's Crust. AGU, 169-184.
[3] Chapman D S, Keho T H, Bauer M S, et al. 1984. Heat flow in the Uinta Basin determined from bottom hole temperature(BHT) data. Geophysics, 49(4):453-466.
[4] Chen F Y, Ding W L, Qian Z, et al. 2013. Control factors of burial-hydrocarbon generation history of Carboniferous-Permian source rocks in Jizhong Depression. Fault-Block Oil & Gas Field(in Chinese), 20(2):145-149.
[5] Chen M X, Huang G S, Zhang W R, et al. 1982. The temperature distribution pattern and the utilization of geothermal water at Niutuozhen Basement protrusion of Central Hebei Province. Scientia Geologica Sinica(in Chinese),(3):239-252.
[6] Cui Z Q. 2005. Depositional systems and exploration of subtle oil-gas pools of Paleogene in Jizhong Depression, Bohai Bay Basin[Ph. D. thesis](in Chinese). Xi'an:Northwest University.
[7] Ding Z Y, Wang L S, Hu W X, et al. 2008. Reconstruction of Cenozoic thermal history of Bohai Bay Basin with a transient heat flow model. Acta Petrolei Sinica(in Chinese), 29(5):650-656.
[8] Dong D W, Li L, Liu J, et al. 2013. Cenozoic tectonic evolution in the north-central Jizhong Depression. Oil & Gas Geology(in Chinese), 34(6):771-780.
[9] Du J H, Zou W H, Fei B S, et al. 2002. Buried Hill Oil and Gas Complex Zone in Jizhong Depression(in Chinese). Beijing:Science Press.
[10] Feng C G, Liu S W, Wang L S, et al. 2009. Present-day geothermal regime in Tarim Basin, northwest China. Chinese J. Geophys.(in Chinese), 52(11):2752-2762, doi:10.3969/j.issn.0001-5733.2009.11.010.
[11] Feng S, Hao S S. 1983. Paleotemperature evolution and petroleum generation potential of Raoyang and Langgu subsidences of North China Basin. Journal of China University of Petroleum (in Chinese),(1):1-10.
[12] Gong Y L, Wang L S, Liu S W, et al. 2003. Distribution characteristics of geotemperature field in Jiyang Depression, Shandong, North China. Chinese J. Geophys. (in Chinese), 46(5):652-658.
[13] Gong Y L. 2003. The thermal structure and thermal evolution of Bohai Bay Basin in East China[Ph. D. thesis](in Chinese). Nanjing:Nanjing University.
[14] Guo X W, He S, Liu K Y, et al. 2011. Modelling the petroleum generation and migration of the third member of the Shahejie Formation(Es3) in the Banqiao Depression of Bohai Bay Basin, Eastern China. Journal of Asian Earth Sciences,40(1):287-302.
[15] He L J. 2014. The Rheological boundary layer and its implications for the difference between the thermal and seismic lithospheric bases of the North China Craton. Chinese J. Geophys. (in Chinese), 57(1):53-61, doi:10.6038/cjg20140106.
[16] He L J. 2015. Thermal regime of the North China Craton:Implications for craton destruction. Earth-Science Reviews, 140:14-26.
[17] Hou G T, Qian X L, Cai D S. 2001. The tectonic evolution of Bohai Basin in Mesozoic and Cenozoic Time. Acta Scientiarum Naturalium Universitatis Pekinensis(in Chinese), 37(6):845-851.
[18] Li W Z, Jiao Y X, Zuo Y H, et al. 2014. Effect of deposition rate on geothermal field in Bozhong depression, Bohai Bay Basin. Chinese J. Geophys. (in Chinese), 57(5):1568-1577, doi:10.6038/cjg20140520.
[19] Liang H B, Qian Z, Xin S L, et al. 2010. Assessment and development of geothermal resources in Jizhong Depression. New Energy(in Chinese),(5):63-68.
[20] Liang S J. 2001. The characteristics of tectonics and hydrocarbon accumulation of Jizhong Depression in Late Cenozoic Era[Master's thesis](in Chinese). Xi'an:Northwest University.
[21] Lin G C, Nunn J A, Deming D. 2000. Thermal buffering of sedimentary basins by basement rocks:implications arising from numerical simulations. Petroleum Geoscience, 6(4):299-307.
[22] Lin S H, Gong Y L. 2005. Distribution Characteristics of Geotemperature field in Jizhong Depression, North China. Journal of East China Institute of Technology (in Chinese), 28(4):359-364.
[23] Liu S W, Wang L S, Gong Y L, et al. 2005. Thermal-rheological structure of the lithosphere beneath Jiyang depression:its implications for geodynamics. Science China Earth Sciences, 48(10):1569-1584.
[24] Lu S K, Li J Y, Wu K Y, et al. 2011. Tectonic evolution of buried hill in Jizhong Depression and the petroleum geological significance. Journal of Oil and Gas Technology(in Chinese), 33(11):35-40.
[25] Morgan P. 1984. The thermal structure and thermal evolution of the continental lithosphere. Physics and chemistry of the Earth, 15:107-193.
[26] Qi J F, Zhang Y W, Lu K Z, et al. 1995. Cenozoic tectonic evolution in Bohai Bay Basin province. Journal of the University of Petroleum, China(in Chinese), 19(S):1-6.
[27] Qiu N S, Su X G, Li Z Y, et al. 2007. The Cenozoic tectono-thermal evolution of depressions along both sides of mid-segment of Tancheng-Lujiang Fault Zone, East China. Chinese J. Geophys.(in Chinese), 50(5):1497-1507.
[28] Qiu N S, Wei G, Li C C, et al. 2009. Distribution features of current geothermal field in the Bohai Sea waters. Oil & Gas Geology(in Chinese), 30(4):412-419.
[29] Qiu N S, Zuo Y H, Chang J, et al. 2014. Geothermal evidence of Meso-Cenozoic lithosphere thinning in the Jiyang sub-basin, Bohai Bay Basin, eastern North China Craton. Gondwana Research, 26(3-4):1079-1092.
[30] Qiu N S, Zuo Y H, Chang J, et al. 2015. Characteristics of Meso-Cenozoic thermal regimes in typical eastern and western sedimentary basins of China. Earth Science Frontiers (in Chinese), 22(1):157-168.
[31] Rao S, Hu S B, Zhu C Q, et al. 2013. The characteristics of heat flow and lithospheric thermal structure in Junggar Basin, northwest China. Chinese J. Geophys. (in Chinese), 56(8):2760-2770, doi:10.6038/cjg20130824.
[32] Sun D S, Liu C Y, Yang M H, et al. 2004a. Study on complex extensional structures in the middle Jizhong Depression in the Bohai Bay Basin. Geological Review(in Chinese), 50(5):484-491.
[33] Sun D S, Liu C Y, Yang M H, et al. 2004b. Evidences for the large-scale detachment structure in middle area of Jizhong Depression in Mid-Late Mesozoic period. Geotectonica et Metallogenia(in Chinese), 28(2):126-133.
[34] Wang J A, Wang J Y, Xiong L P, et al. 1992. Geothermal characteristics and petroleum resources of Liaohe Graben.//Bulletin(5) of Institute of Geology, Chinese Academy of Sciences(in Chinese). Beijing:Science Press, 1-77.
[35] Wang J Y. 1996. Geothermics in China. Beijing:Seismological Press.
[36] Wang L S, Liu S W, Xiao W Y, et al. 2002. Distribution characteristics of terrestrial heat flow in Bohai Basin. Chinese sciences Bulletin (in Chinese), 47(2):151-155.
[37] Wang Y, Cheng S H. 2011. Thermal state and rheological strength of the lithosphere beneath the Eastern China. Geotectonica et Metallogenia (in Chinese), 35(1):12-23.
[38] Wu Z P, Hou X B, Li W. 2007. Discussion on Mesozoic basin patterns and evolution in the eastern North China Block. Geotectonica et Metallogenia (in Chinese), 31(4):385-399.
[39] Xiong L P, Gao W A. 1982. Characteristics of geotherm in uplift and depression. Chinese J. Geophys.(Acta Geophysica Sinica)(in Chinese), 25(5):448-456.
[40] Xu J, Gao Z W, Sun J B, et al. 2001. A preliminary study of the coupling relationship between basin and mountain in extensional environments-a case study of the Bohai Bay Basin and Taihang Mountain. Acta Geologica Sinica(in Chinese), 75(2):165-174.
[41] Xu S Y, Yan K. 2005. Structural system and hydrocarbon distribution in the Bohai Gulf Basin. Journal of Geomechanics(in Chinese), 11(3):259-265.
[42] Yang M H, Liu C Y, Yang B Y. 2001. Tectonic transform of the Mesozoic and their relations to hydrocarbon traps in the Jizhong(Central Hebei Province) Basin, China. Geotectonica et Metallogenia(in Chinese), 25(2):113-119.
[43] Yang M H, Liu C Y, Yang B Y, et al. 2002. Extensional structures of the Paleogene in the Central Hebei Basin, China. Geological Review(in Chinese), 48(1):58-67.
[44] Zang S X, Liu Y G, Ning J Y. 2002. Thermal structure of the lithosphere in North China. Chinese J. Geophys.(in Chinese), 45(1):56-66.
[45] Zhang W C, Yang D X, Chen Y J, et al. 2008. Sedimentary structural characteristics and hydrocarbon distributed rules of Jizhong depression. Acta Geologica Sinica(in Chinese), 82(8):1103-1112.
[46] Zhao L J, Jiang Y L, Liu H, et al. 2012. Thermal evolution of Paleogene source rocks and relationship with reservoir distribution in Raoyang sag, Bohai Bay Basin. Petroleum Geology and Recovery Efficiency(in Chinese), 19(4):1-4.
[47] Zhao X Z, Jin F M, Wang Q, et al. 2014. Exploration and discovery of subtle deep buried hill and buried hill inner-curtain hydrocarbon reservoir in Jizhong Sag. China Petroleum Exploration (in Chinese), 19(1):10-21.
[48] Zhao X Z, Jin F M, Zhang Y M, et al. 2010. Fine Exploration of Unconventional Buried Hill Reservoirs in Oil-Enriched Sags(in Chinese). Beijing:Petroleum Industry Press.
[49] Zhao X Z, Jin F M, Wang Q, et al. 2015. Buried-hill play, Jizhong subbasin, Bohai Bay basin:A review and future prospectivity. AAPG Bulletin, 99(1):1-26.
[50] Zhou R L. 1987. The activity of deep underground water in the northern part of the North China Plain and its effect on the geothermal field. Bulletin of the 562 Comprehensive Geological Brigade Chinese Academy of Geological Sciences(in Chinese),(6):17-35.
[51] Zhu R X, Chen L, Wu F Y, et al. 2011. Timing, scale and mechanism of the destruction of the North China Craton. Science China Earth Sciences, 54(6):789-797.
[52] Zuo Y H. 2011. Geothermal evidence of the Mesozoic and Cenozoic lithospheric evolution in the Bohai Bay Basin, North China[Ph. D. thesis](in Chinese). Beijing:China University of Petroleum.
[53] Zuo Y H, Qiu N S, Hao Q Q, et al. 2014. Present geothermal fields of the Dongpu Sag in the Bohai Bay Basin. Acta Geologica Sinica-English Edition, 88(3):915-930.
[54] 尘福艳, 丁文龙, 钱铮等. 2013. 冀中坳陷石炭-二叠系烃源岩埋藏-生烃史控制因素. 断块油气田, 20(2):145-149.
[55] 陈墨香, 黄歌山, 张文仁等. 1982. 冀中牛驼镇凸起地温场的特点及地下热水的开发利用. 地质科学,(3):239-252.
[56] 崔周旗. 2005. 渤海湾盆地冀中坳陷古近系沉积体系与隐蔽油气藏勘探[博士论文]. 西安:西北大学.
[57] 丁增勇, 王良书, 胡文瑄等. 2008. 利用岩石圈瞬态热流模型重建渤海湾盆地新生代热历史. 石油学报, 29(5):650-656.
[58] 董大伟, 李理, 刘建等. 2013. 冀中坳陷中北部新生代构造演化特征. 石油与天然气地质, 34(6):771-780.
[59] 杜金虎, 邹伟宏, 费宝生等. 2002. 冀中拗陷古潜山复式油气聚集区. 北京:科学出版社.
[60] 冯昌格, 刘绍文, 王良书等. 2009. 塔里木盆地现今地热特征. 地球物理学报, 52(11):2752-2762, doi:10.3969/j.issn.0001-5733.2009.11.010.
[61] 冯石, 郝石生. 1983. 饶阳凹陷和廊固凹陷的古地温演化与生油气潜力. 华东石油学院学报,(1):1-10.
[62] 龚育龄, 王良书, 刘绍文等. 2003. 济阳坳陷地温场分布特征. 地球物理学报, 46(5):652-658.
[63] 龚育龄. 2003. 中国东部渤海湾盆地热结构和热演化[博士论文]. 南京:南京大学.
[64] 何丽娟. 2014. 流变边界层及其对华北克拉通热/地震岩石圈底界差异的意义. 地球物理学报, 57(1):53-61, doi:10.6038/cjg20140106.
[65] 侯贵廷, 钱祥麟, 蔡东升. 2001. 渤海湾盆地中、新生代构造演化研究. 北京大学学报(自然科学版), 37(6):845-851.
[66] 李文正, 焦亚先, 左银辉等. 2014. 沉积速率对渤中坳陷大地热流的影响. 地球物理学报, 57(5):1568-1577, doi:10.6038/cjg20140520.
[67] 梁宏斌, 钱铮, 辛守良等. 2010. 冀中坳陷地热资源评价及开发利用. 新能源,(5):63-68.
[68] 梁苏娟. 2001. 冀中坳陷晚新生代地质构造特征及其油气赋存[硕士论文]. 西安:西北大学.
[69] 林世辉, 龚育龄. 2005. 冀中坳陷现今地温场分布特征. 东华理工学院学报, 28(4):359-364.
[70] 刘绍文, 王良书, 龚育龄等. 2005. 济阳坳陷岩石圈热-流变学结构及其地球动力学意义. 中国科学D辑——地球科学, 35(3):203-214.
[71] 陆诗阔, 李继岩, 吴孔友等. 2011. 冀中坳陷潜山构造演化特征及其石油地质意义. 石油天然气学报, 33(11):35-40.
[72] 漆家福, 张一伟, 陆克政等. 1995. 渤海湾盆地新生代构造演化. 石油大学学报(自然科学版), 19(S):1-6.
[73] 邱楠生, 苏向光, 李兆影等. 2007. 郯庐断裂中段两侧坳陷的新生代构造-热演化特征. 地球物理学报, 50(5):1497-1507.
[74] 邱楠生, 魏刚, 李翠翠等. 2009. 渤海海域现今地温场分布特征. 石油与天然气地质, 30(4):412-419.
[75] 邱楠生, 左银辉, 常健等. 2015. 中国东西部典型盆地中-新生代热体制对比. 地学前缘, 22(1):157-168.
[76] 饶松, 胡圣标, 朱传庆等. 2013. 准噶尔盆地大地热流特征与岩石圈热结构. 地球物理学报, 56(8):2760-2770, doi:10.6038/cjg20130824.
[77] 孙冬胜, 刘池阳, 杨明慧等. 2004a. 渤海湾盆地冀中坳陷中区中新生代复合伸展构造. 地质论评, 50(5):484-491.
[78] 孙冬胜, 刘池阳, 杨明慧等. 2004b. 冀中坳陷中区中生代中晚期大型拆离滑覆构造的确定. 大地构造与成矿学, 28(2):126-133.
[79] 汪缉安, 汪集旸, 熊亮萍等. 1992. 辽河断陷地热特征与油气资源.//中国科学院地质研究所集刊(5). 北京:科学出版社, 1-77.
[80] 王良书, 刘绍文, 肖卫勇等. 2002. 渤海盆地大地热流分布特征. 科学通报, 47(2):151-155.
[81] 汪洋, 程素华. 2011. 中国东部岩石圈热状态与流变学强度特征. 大地构造与成矿学, 35(1):12-23.
[82] 吴智平, 侯旭波, 李伟. 2007. 华北东部地区中生代盆地格局及演化过程探讨. 大地构造与成矿学, 31(4):385-399.
[83] 熊亮萍, 高维安. 1982. 隆起与坳陷地区地温场的特点. 地球物理学报, 25(5):448-456.
[84] 徐杰, 高战武, 孙建宝等. 2001. 区域伸展体制下盆-山构造耦合关系的探讨——以渤海湾盆地和太行山为例. 地质学报, 75(2):165-174.
[85] 徐守余, 严科. 2005. 渤海湾盆地构造体系与油气分布. 地质力学学报, 11(3):259-265.
[86] 杨明慧, 刘池阳, 杨斌谊. 2001. 冀中坳陷中生代构造变形的转换及油气. 大地构造与成矿学, 25(2):113-119.
[87] 杨明慧, 刘池阳, 杨斌谊等. 2002. 冀中坳陷古近纪的伸展构造. 地质论评, 48(1):58-67.
[88] 臧绍先, 刘永刚, 宁杰远. 2002. 华北地区岩石圈热结构的研究. 地球物理学报, 45(1):56-66.
[89] 张文朝, 杨德相, 陈彦均等. 2008. 冀中坳陷古近系沉积构造特征与油气分布规律. 地质学报, 82(8):1103-1112.
[90] 赵利杰, 蒋有录, 刘华等. 2012. 饶阳凹陷烃源岩热演化特征及其与油藏分布的关系. 油气地质与采收率, 19(4):1-4.
[91] 赵贤正, 金凤鸣, 王权等. 2014. 冀中坳陷隐蔽深潜山及潜山内幕油气藏的勘探发现与认识. 中国石油勘探, 19(1):10-21.
[92] 赵贤正, 金凤鸣, 张以明等. 2010. 富油凹陷隐蔽型潜山油气藏精细勘探. 北京:石油工业出版社.
[93] 周瑞良. 1987. 华北平原北部深层地下水活动及其对地温场的影响. 中国地质科学院562综合大队集刊,(6):17-35.
[94] 左银辉. 2011. 渤海湾盆地中、新生代岩石圈演化的地热学证据[博士论文]. 北京:中国石油大学.