E-mail Alert Rss
 

物探与化探, 2024, 48(4): 934-944 doi: 10.11720/wtyht.2024.1316

地质调查·资源勘查

反射地震剖面揭示容城凸起基岩面和断裂结构及其对地热成因的启示

刘宏凯,1, 高磊1, 张杰2, 侯贺晟1, 谢民英3, 李洪强,1

1.中国地质科学院,北京 100037

2.中国地质科学院 地球物理地球化学勘查研究所,河北 廊坊 065000

3.中国自然资源航空物探遥感中心,北京 100083

Bedrock surface and fault structures in the Rongcheng uplift revealed from reflection seismic profiles and their implications for the geothermal origin

LIU Hong-Kai,1, GAO Lei1, ZHANG Jie2, HOU He-Sheng1, XIE Min-Ying3, LI Hong-Qiang,1

1. Chinese Academy of Geological Sciences, Beijing 100037,China

2. Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000, China

3. China Aero Geophysical Survey and Remote Sensing Center for Natural Resources, Beijing 100083, China

通讯作者: 李洪强(1984-),男,博士,副研究员,主要从事深部地质结构探测工作。Email:hongqiangli@126.com

第一作者: 刘宏凯(2000-),男,中国地质科学院硕士研究生,主要研究方向为主动源地震探测。Email:CAGSLiuHongkai@163.com

责任编辑: 叶佩, 沈效群

收稿日期: 2023-07-25   修回日期: 2024-02-25  

基金资助: 国家重点研发计划项目(SQ2023YFF0800060)
国家自然科学基金项目(42074097)
国家自然科学基金项目(42230303)
国家自然科学基金项目(41904083)
国家自然科学基金项目(41704089)
国家自然科学基金项目(42004079)
中国地质调查局项目(DD20230008)
中国地质调查局项目(DD20189132)

Received: 2023-07-25   Revised: 2024-02-25  

摘要

容城凸起是华北地区地热资源富集区域之一,探测容城凸起基岩面的精细结构、埋深和断裂,对理解该区地热资源的分布和富集以及地热资源的勘查与开发具有重要意义。本文基于中国地质科学院2018年在容城凸起采集的SN向反射地震剖面数据,通过对反射数据精细处理,获得容城凸起4 km以浅的高精度几何结构,经地热孔资料标定后开展解释。主要获得以下4点认识:容城凸起基岩面以上新生代沉积地层呈近水平层状分布,这些近水平地层对容城地热田起到保温盖层作用;容城凸起基岩面埋深在700~3 000 m之间,其凸起中心部位变化相对平缓,向周边迅速加深到3 000 m左右;牛南断裂和容东断裂在深部延伸交汇,与其他中、小型断裂组成的断裂系统为水和热的传导提供了通道和空间;基岩面凸起的几何结构有利于热流在凸起下方汇聚。

关键词: 容城凸起; 反射地震剖面; 基岩面; 地热; 断裂系统

Abstract

The Rongcheng uplift in North China boasts abundant geothermal resources. Research indicates that the Rongcheng uplift exhibits significantly different physical properties between the bedrock surface and the overlying Cenozoic strata. Moreover, the bedrock surface serves as the primary top boundary of the geothermal reservoir in the Wumishan Formation. Investigating the fine-scale structures, burial depths, and faults of the bedrock surface in the Rongcheng uplift holds critical significance for understanding the distribution and enrichment of geothermal resources in the area and guiding their exploration and production. Through elaborative processing of the north-south reflection seismic profile data of the Rongcheng uplift, collected by the Chinese Academy of Geological Sciences in 2018, this study obtained the high-precision geometric structure of the Rongcheng uplift within a depth of 4 km. The geometric structure was calibrated using geothermal borehole data before interpretation. Key findings are as follows: (1) The Cenozoic sedimentary strata overlying the bedrock surface of the Rongcheng uplift exhibit a nearly horizontal layered distribution, serving as cap rocks of the Rongcheng geothermal field; (2) The bedrock surface of the Rongcheng uplift manifests burial depths ranging from 700 to 3 000 m, with gentle changes in the central portion, and rapidly deepening to around 3 000 m towards the periphery; (3) The Niunan and Rongdong faults converge in the deep part, constituting a fault system along with other medium and small faults, thus facilitating the conduction of water and heat; (4) The geometric structure of the Rongcheng uplift on the bedrock surface contributes to the convergence of heat flow beneath the uplift.

Keywords: Rongcheng uplift; reflection seismic profile; bedrock surface; geothermal; fault system

PDF (9585KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘宏凯, 高磊, 张杰, 侯贺晟, 谢民英, 李洪强. 反射地震剖面揭示容城凸起基岩面和断裂结构及其对地热成因的启示[J]. 物探与化探, 2024, 48(4): 934-944 doi:10.11720/wtyht.2024.1316

LIU Hong-Kai, GAO Lei, ZHANG Jie, HOU He-Sheng, XIE Min-Ying, LI Hong-Qiang. Bedrock surface and fault structures in the Rongcheng uplift revealed from reflection seismic profiles and their implications for the geothermal origin[J]. Geophysical and Geochemical Exploration, 2024, 48(4): 934-944 doi:10.11720/wtyht.2024.1316

0 引言

充分利用地热资源对绿色低碳城市建设具有重要意义。雄安新区地热资源分布广泛,具有埋藏浅、储量大、水质优的特点,是我国中东部地热资源开发利用条件较好的区域[1-2]。容城地热田是雄安新区3个大中型地热田之一,该地热田东部以容东断裂为界,南部以牛南断裂为界,西部以容西断裂为界[3-4],位于热背景值较高的华北东部热盆地区[1,5-6]。容城地热田主要地热储层是元古代蓟县系雾迷山组和高于庄组[4]。已有研究表明容城凸起的基岩面与上部新生代地层存在显著的物性差异,该基岩面同时也是地热储层雾迷山组的主要顶界面。目前对容城地热田的研究主要是从开发潜力评估、大范围地热特征等方向进行,很少聚焦地下构造来对容城地热田进行讨论[7-8]。反射地震剖面可以提供清晰的地下构造几何形态,如断裂分布、地层形态等[9],而且基岩面结构、断裂等均与地热系统有着紧密的关系,所以,探测容城凸起基岩面的精细结构、埋深和断裂,对理解该区地下结构对地热资源分布、富集的影响有重要意义。本文基于中国地质科学院于2018年采集的一条NS向过容城凸起的反射地震剖面,通过对深地震反射数据进行系统处理,获得了容城凸起3 s(约4 km)以浅地层的高精度几何结构,结合钻孔资料标定,研究容城凸起地热田的结构,探讨基岩面结构与断裂对地热田的影响。

1 地质背景

1.1 构造背景

容城凸起构造位置上属于华北断坳(Ⅱ级)、冀中坳陷(Ⅲ级)的中区。冀中坳陷为渤海湾盆地的次级构造单元,是发育在华北克拉通基底之上近SN向的中新生代断陷—坳陷盆地。容城凸起属于冀中坳陷中的一个次级正向构造,由于处于华北克拉通破坏的中心部位置,该凸起经历了克拉通演化的主要过程。自太古宙以来,渤海湾盆地经历了多次构造运动[10-12] 。太古宙至古元古代为克拉通结晶基底的形成阶段 [13];新元古代至古生代,构造活动以整体性的升降运动为主[14]。受印支运动的影响,渤海湾盆地受到近SN向挤压而转变为近EW向伸展 [15],并且在早白垩世出现大规模的伸展构造[16]。中生代晚期,受太平洋板块俯冲的影响,克拉通东部岩石圈发生巨量减薄[17];新生代早期,冀中坳陷进入断陷阶段,处于冀中坳陷的雄安新区在古近纪早期也开始发生断陷活动;古新世至早始新世,冀中坳陷内容城、徐水等处开始一系列正断层活动,并发育徐水凹陷、容城凸起等次级构造单元[18];新近纪至晚第四纪,由于华北地区的弧后拉张作用减弱,冀中坳陷地区进入坳陷演化阶段,以整体沉降为主 [19],冀中坳陷的活动明显减弱;晚更新世以来,冀中坳陷地质构造相对稳定 [20]

容东断裂和牛南断裂是研究区内2条较大断裂(图1中F1和F2)。容东断裂位于容城东侧,其上断点位于容城县以东西里村北部附近;在地震测线附近其走向为NE向,延伸约15 km,倾向SE;在大河镇东部附近,该断裂走向变化为近NS向,倾向近EW向,延伸约7.8 km[21]。该断层切穿古近系和中—新元古界地层,上盘为白洋淀洼陷内古近系地层,下盘为容西凸起内中—新元古界地层;断层形态上陡下缓,与太行山山前断裂形态一致,推测为同时期在拉张应力环境下形成的伸展构造,主要控制白洋淀洼陷的发育。牛南断裂位于安新县至雄县赵北口南一带,是控制牛驼镇凸起西南边界的正断裂,长约23 km,断层走向近EW,断面南倾,倾角45°左右,垂直段距1 200~3 200 m,水平段距1 000~2 500 m。此断裂控制古近系的厚度和构造,断开了结晶基底,为长期活动的大断裂[22]

图1

图1   容城构凸起及周缘主要构造单元和断裂分布(底图据文献[23]修改)

Fig.1   Rongcheng tectonic unit and main fault distribution map(revised according to reference[23])


当前,对容东断裂的性质和分布认识相对清楚,但受限于高精度探测资料的不足,对牛南断裂与容东断裂的深部延伸尚不清晰。

根据容城凸起内的4口地热井岩性、地层和物性资料分析(图2),容城凸起区地层由浅至深为第四系,新近系明化镇组、馆陶组,古近系东营组、沙河街组,蓟县系雾迷山组、杨庄组,长城系高于庄组、大红峪组,以及太古界变质岩和中生界、古生界。第四系主要为灰黄、黄灰色砂泥互层,新近系明化镇组厚度约340~640 m,棕红色泥岩与杂色砂岩互层,测井曲线整体表现为中低伽马、低阻、低波速特征,从第四系至新近系明化镇组的电阻率整体呈降低趋势,声波速度整体呈增大趋势。古近系厚度约为934 m,主要发育泥岩与砂岩,测井曲线整体表现为低伽马、低阻、低波速。元古界蓟县系雾迷山组厚度约100~845 m,以厚层灰色、灰白色白云岩为主,以极低的自然伽马测井响应为典型特征,电阻率、声波速度整体呈中高特征,多见指状低阻、低波速异常,显示出雾迷山组良好的裂隙发育特征;杨庄组地层较薄,厚度常小于50 m,岩性为棕红色(白云质)泥岩,测井曲线整体表现为中低伽马、低阻、中低速特征,与上部雾迷山组和下部长城系地层物性差异明显,为典型的隔水层;长城系高于庄组在该区分布不均,部分钻孔缺失,钻遇最厚约1 226 m,以灰色白云岩夹白云质泥岩为主,同样以极低的自然伽马响应为典型特征,电阻率、声波速度整体呈中高特征,多见指状伽马高异常,偶见指状低阻、低波速异常;长城系大红峪组上部岩性以玄武岩为主,下部岩性以白云岩为主,该区钻遇最厚366 m,测井曲线上部表现为高伽马、低阻、高波速特征,下部表现为低伽马、高阻、高波速特征。其中,D13揭露了太古界变质岩地层,以片麻岩为主,整体为高伽马、高阻、高波速特征。这些为后续反射地震处理、标定和解释奠定了基础。

图2

图2   容城凸起综合钻井柱状图(据文献[13,24]修改)

Fig.2   Comprehensive drilling histogram of Rongcheng uplift(revised according to references [13,24])


1.2 地热储层特征

容城地热田热储层包括新近系明化镇组孔隙型砂岩热储以及寒武系、蓟县系基岩热储。其中,明化镇组热储由于离地表近,难以回灌而全面禁止开采[4];寒武系热储在研究区域分布小,热储裂隙不发育,不属于地热田的主要热储层;蓟县系雾迷山组与高于庄组是容城地热田主要开发的深部基岩热储层,地层岩性主要为灰白色含燧石条带白云岩,是热传导和热对流共同作用下形成的裂隙型深部热储层。从岩心CT扫描结果可知,蓟县系的裂缝长度和宽度较大,和孔洞相连形成良好的缝洞型储存空间[25] 。该地层由于长期受到风化、剥蚀、淋滤的作用,裂隙较为发育,为地下热水储存和运移提供空间[26]

容城凸起与牛驼镇凸起区域的雾迷山组顶界埋深一般在2 500 m以浅,而其他区域埋深较深,最深超过3 000 m;雾迷山组厚度约为100~845 m。雾迷山组的孔隙度主要分布于6%之内,占比85.9%,最大孔隙度可达22.4%;渗透率主要分布于(0.01~100)×10-3 μm2,占比87.3%;地温梯度为0.2~1.0 ℃/100 m。深度在1 100~2 100 m之间的井口水温为52~61 ℃,热水出水量为80~120 m3/h,矿化度为1 900~3 100 mg/L [22]。高于庄组顶界埋深为2 078 m,揭露厚度925.24 m,孔隙度在1.82%~3.76%之间,平均渗透率为1.43×10-3μm2,井口温度71℃;该地层是容城地热田的另一主力热储[27]

2 数据采集及处理

2.1 反射地震数据采集

反射地震测线南起保定凹陷端村镇东南部,经容城凸起,北至梁家营镇北部(图1)。测线满覆盖长度为22.91 km,采用炸药震源激发、多次覆盖方式进行数据采集。震源井深3 045 m,药量10~18 kg,采用中间激发、两边对称接收的观测系统,检波点采用10个10 Hz检波器串联组合,检波点距40 m,720道接收,采样率间隔2 ms,记录长度30 s,满覆盖次数180次,详细采集参数见表1

表1   容城凸起反射地震剖面数据采集参数

Table 1  Parameters of Rongcheng uplift reflection seismic profile data

采集系统检波器
主频/Hz
检波点
间距/m
炮点
深度/m
接收
道数
记录长
度/s
最小偏
移距/m
最大偏
移距/km
炮间距/m覆盖次数药量/kg采样间
隔/ms
Sercel 428XL104030~45720304014.388018010~182

新窗口打开| 下载CSV


由于测线经过区域村庄众多、交通发达,人类活动范围广,工业干扰比较严重,导致数据的线性干扰、50 Hz工业干扰和随机噪声发育。但采集的单炮数据浅层资料连续性较好,信噪比高,可清晰看到基岩面反射信号。由于此次研究针对容城凸起的浅层地质结构,所以选取了3 s以上、具有代表性的原始地震单炮数据(图3)。

图3

图3   原始单炮记录(3 s内)

Fig.3   Original single-shot record map(within 3 s)


2.2 数据处理

围绕研究区的基岩面结构和其物性特征,对本区地震地质条件和数据进行分析,而数据处理的关键是静校正和噪声衰减。运用初至波层析进行静校正,并通过噪声衰减与低频面波去除的组合去噪方法来衰减噪声、提高信噪比。反射数据经静校正、噪声压制、速度分析、动校正、叠加、偏移等处理,获得了高信噪比成果剖面。处理流程见图4

图4

图4   研究区采用的反射地震数据处理流程

Fig.4   Reflection seismic data processing flow


3 结果分析

经数据处理后获得的高信噪比成果剖面——反射地震剖面(图5),清晰反映了容城凸起地层的结构特征,剖面整体信噪比高,由浅至深反射界面清晰,反射能量强、同相轴连续性好,有利于进行地层划分、构造解释与钻井资料对比。结合D17钻孔数据对地震资料进行标定,可以获得地层、断裂的空间展布与岩性特征。

图5

图5   容城凸起反射地震偏移剖面

Fig.5   The reflection seismic migration profile of Rongcheng uplift


3.1 地层结构

地震剖面解释如图6所示。地层由浅至深主要有第四系、新近系明化镇组(上下段)、新近系馆陶组、古近系东营组、古近系沙河街组(上下段)、蓟县系雾迷山组。通过强反射界面(明化镇组底界)可将剖面分为上层和下层两个部分,剖面上层部分主要为近水平层状地层,整体表现为从北向南倾斜;下层部分CMP10~1 000之间的深部地层产状较为复杂,并且有多条断裂穿切地层,CMP1 000~1 800深部地层没有明显的反射同相轴,且内部反射信号振幅频率均匀,推测为岩体。

图6

图6   反射地震剖面及其解释结果

Fig.6   seismic reflection profile interpretation map


第四系(Q)在研究区内全区分布,岩性主要为灰黄、黄灰色的泥砂岩互层,表现为上粗下细沉积韵律,结构较为松散。厚度分布比较均匀,在凸起区域受剥蚀有一定减薄,与下方新近系明化镇组呈不整合接触,整体呈现为近平行的地层结构。

新近系明化镇组(Nm)全区沉积,岩性主要发育棕红色泥岩与灰黄色砂岩互层,沉积地层整体由厚转薄,变化较为缓慢。由于泥岩和砂岩互层的影响,明化镇组层内可见多组较强的反射同相轴。剖面显示明化镇组的底界是一条强能量的连续同相轴,该界面同时也是容城凸起下方的基岩面,下上地层之间巨大的物性差异使其成为剖面上反射能量最强的一条反射界面,清晰地呈现了容城凸起基岩面的构造形态。D17钻井资料显示该界面埋深为1 520 m,从南向北埋深由1 800 m缓慢变浅,在CMP1 580左右逐渐过渡到1 500 m之后,构造形态开始急剧变化,基岩面发育成背斜构造,背斜顶部埋深大约在1 200 m。隆起的基岩面下方地层振幅和频率分布均较为均匀,可推测下方为基岩地层,并且由于基岩面的隆起,导致上方的明化镇组地层出现了强烈的剥蚀减薄作用,最薄处厚度仅为300 m。

新近系馆陶组(Ng)在研究区广泛分布,该地层主要发育有杂色砂岩、泥岩,底部发育有砾岩,与下伏地层呈角度不整合接触。由于地质作用和断裂的影响,馆陶组底界形态起伏较大,沿地震剖面由南向北埋深从2 400 m左右逐渐过渡到1 800 m,至F2断裂后地层沉积逐渐增厚,在基岩面隆起处厚度达到400 m左右。

古近系东营组(Ed)在研究区内主要分布于保定凹陷处,在容城凸起上方已经全部剥蚀,主要发育浅灰色砂岩,紫红色、褐色泥岩。钻井揭示该地层顶部界面埋深在1 624 m,底界埋深沿测线从2 550 m先向下增加再变浅,CMP 900后已经被完全剥蚀。东营组平均厚度较薄,整体呈由北向南倾斜的层状地层。

古近系沙河街组(Es)上段主要分布于保定凹陷,下段分布较为广泛。该地层岩性以砂岩、泥岩、泥质砂岩为主。沙河街组沉积厚度较大,钻井揭示该地层顶部界面埋深为1 948 m,并且从剖面可以观察到该地层于蓟县系雾迷山组呈不整合接触。层内存在多套反射能量强,同相轴清晰,发育多条中小断裂,推测是由于受到当时强烈地质作用的影响所导致。

蓟县系雾迷山组(Jxw)分布广泛,主要发育灰白色、含燧石条带的裂隙白云岩。该地层表现为巨厚沉积,形成了优良的热储条件。钻井揭示该地层顶部界面埋深在2 558 m,且变化较为平缓,顶界埋深没有明显起伏,地层内未见明显的反射同相轴。

3.2 断裂结构

研究区断裂发育,基岩面以浅发育许多中小断裂,多为产状相似的正断裂,均切穿新近系明化镇组,垂直断距最大约2 000 m,最小约500 m。基岩面以深的中小断裂基本在古近系地层内发育,该地层断裂产状各异,倾角较大,断裂之间相互连通,组成了一个较大范围的断裂系统。

研究区内大断裂沿地震剖面主要分布有F1、F2、F3和F4(图6)。F1为牛南断裂,走向近EW,倾向NE,上断点埋深500 m左右,垂直断距2 000 m,延伸到蓟县系雾迷山组,控制研究区内第四系至古近系的沉积,表现为典型的上陡下缓产状构造。F2为容东断裂,走向近SN,倾向SW,倾角50°,上断点埋深600 m左右,垂直断距超过3 000 m,切穿剖面上所有的反射界面。地震剖面显示F1与F2在蓟县系雾迷山组交汇,二者共同组成较大断裂带。F3上断点埋深在1 500 m左右,垂直断距超过2 000 m,倾向NE,倾角约55°。该断裂切穿新近系馆陶组、古近系沙河街组和蓟县系雾迷山组。F4上断点埋深约为1 300 m,垂直断距2 000 m,倾向NE,倾角约60°。该断裂顶部伴生有小型断裂,呈斜列状展布,共同组成“Y”字型断裂构造。

4 讨论

容城地热田是中国华北地区具有代表性的低温热水型地热田之一[4],有较好的开发前景。通过地震剖面和钻井资料的结合,获得容城凸起南部精细的基岩面、断裂构造结构,研究结果表明容城地热田是多个因素共同作用下形成的结果,主要影响因素有盖层、基岩面和断裂结构。其中盖层具有良好的保温能力,能够起到良好的热阻隔作用;热储结构为地热系统聚热、储热等提供重要支撑;牛南断裂与容东断裂组成的断裂带具有导水、导热的通道作用,能够为地下热水提供储存和对流的空间,并且大断裂与研究区中其他中小型断裂形成的断裂系统更加有利于大规模地热田的形成。

4.1 盖层

第四系—新近系构造层与古近系构造层的岩性主要为致密的泥岩与粉砂岩,两构造层直接覆盖在中元古界构造层发育的裂隙白云岩之上,形成盖层。盖层岩性致密、孔隙较少,且热导率低,对下方的热量具有良好的保存作用。热导率低、致密性高的盖层覆于热导率高、裂隙发育的热储层之上,形成了容城地热田的基本地质条件。除了岩性的影响以外,盖层的厚度对隔热效果也存在影响,若盖层太厚则下部热量难以向上发散。由于容城凸起基岩面向上隆起,上方沉积遭受剥蚀,沉积的盖层较薄,使得容城地热田的热量在良好保存的同时也能有良好的地温效果。

4.2 基岩面结构

基岩面的结构对地热田的影响是多方面的,主要包括:基岩面起伏与构造形态对于地温场分布的影响、埋深对于地温梯度的影响、稳定结构对地热系统的支撑和保护作用,以及高孔隙度岩层对热流优良的储存作用。

基岩的起伏与构造形态对地温场分布的影响是区域性的,对地壳浅部地温起着主导作用。地震波走时携带有传播过程中的地下结构信息,利用反射剖面大量的初至波信息进行层析反演,得到容城凸起下方的纵波速度结构(图7,其中主要层位和断裂参考反射地震剖面)。初至波层析成像方法是利用接收地震波的到时数据反演出地下的速度结构,最终得到的P波速度结构剖面与地震剖面可以较好地结合起来,展示了容城凸起下方纵波速度分布特征和基底的展布。

图7

图7   容城凸起纵波速度剖面

Fig.7   P-wave velocity profile of Rongcheng uplift


图7可以看出基底的展布与地形地貌较为匹配:保定凹陷处基底埋深较深,断裂F2处基底开始迅速抬升,在容城凸起区域有明显背斜构造,这验证了F2是控制保定凹陷和容城凸起边界的重要构造。研究区浅部是一层较薄的松散堆积物,速度一般小于1 500 m/s,最上层的速度结构与第四系地层一致,全区均匀分布,变化较小,但在凸起上方遭受剥蚀而减薄。明化镇组平均速度在2 200 m/s左右,速度层与明化镇组地层形状基本一致,沉积厚度从南向北逐渐减薄,在凸起区域剥蚀减薄作用严重。随着深度的下降,由于重力压实等作用,压力的增加导致岩性发生明显的变化,基岩面(明化镇组下界)下方速度梯度明显增加,地震剖面上基岩面的隆起区和纵波速度结构凸起的高速异常基本重合,且凸起下方的速度梯度达到最大。

已有研究表明,凸起区的热流值一般要高于凹陷区[28]。容城凸起区域地温异常带展布的主体方向也反映了构造的基本轮廓。容城县的地温梯度,在西南部约为<4 ℃/100 m,由西到东升高到>7 ℃/100 m,相对高的地温梯度区的分布范围与下伏的基底构造基本一致。地温异常高低相间分布与地质构造的凸凹分布相一致,最大值与构造高点基本吻合。造成地温展布与构造轮廓相匹配的原因可能与热流的折射有关,表现为在同一深度上热流会从热导率低的盖层凹陷区流向热导率高的基岩凸起区[29]。从图6中可以观察到CMP1430~1800处基岩面向上隆起,馆陶组的砾岩热导率显著高于在同一深度上的明化镇组下段泥岩,热流可以由热导率低的明化镇组下段流向热导率高的馆陶组,在容城凸起下方汇聚。

基岩面结构对地热的影响还体现在其埋深对地温梯度的影响。新生界盖层地温梯度与基岩埋深的关系十分密切,表现为基岩面的埋深变浅使上方的盖层受到剥蚀而发生减薄。盖层的减薄令下方热量得以更好地散发,而地温梯度逐渐增大,这与前文中盖层对地温梯度的影响亦可相互印证。

由于基岩面具有坚硬稳固的特点,对地层能够起到支撑和保护作用,可以很好地维持整个地热系统的稳定,并且高热容的基岩比盖层的热容更大,这为地热系统能够长期稳定提供了基础。

基岩的岩性是影响地热的另一个重要因素。容城凸起下方的蓟县系雾迷山组发育的裂隙白云岩具有极好的储热效果。白云岩本身具有良好的渗透性,加之裂隙的发育使雾迷山组的裂隙白云岩具有极高的热导率,使得热量可以迅速转移给地下水层,基底的热量可以向上运输并且保存,对顶部的盖层通过热传导及热对流的方式进行热量传递。

4.3 断裂

断裂可为热能运移、流体运移富集提供通道,由于断裂带中的裂隙和空隙可以充当流体运移的通道,热水和热能可以沿着断裂带流动和传导。研究区内断裂分析表明,牛南断裂(F1)和容东断裂(F2)是该地区主要的断裂构造。牛南断裂是容城—牛驼镇凸起与保定凹陷边界的断裂构造,而容东断裂则是容城凸起与牛驼镇凸起边界的断裂构造,这两条断裂在蓟县系雾迷山组处交汇在一起,二者共同形成了一组大型断裂带。该断裂带对地热田起影响主要表现在以下3个方面:断裂带由浅至深穿切了整条剖面,下至结晶基底,沟通地下水层与深部热源,提供热流上升通道;该断裂带作为容城凸起与保定凹陷的边界,在控制着容城地热田的空间分布的同时,也防止热流向侧向流失;大型断裂带宽度大,能够为热水对流提供良好的环境,提供更大的储存空间,有利于形成更大规模的地热系统。研究区内相对发育的断裂带对容城地热田起到了重要的影响,它们提供了热流通道、控制了地热田的空间分布,并为热水对流和储存提供了有利条件。但值得注意的是,F1与F2上断点都在盖层内,已错断盖层,这可能是影响容城当地形成中低温型地热田的因素之一。由于断裂切穿盖层,热流在向上运移的过程中对盖层进行加热的同时,也存在一定程度上的热量流失。在此背景下,地下热水的温度也会相对降低。经过长时间演化,当热流值平衡之后,最终形成了现如今中低温型的容城地热田。

5 结论

容城地热田受基岩面结构、盖层、断裂等多因素相互作用的影响。通过对容城凸起反射地震资料的处理解释,获得其下方基岩面和主要断裂的几何形态,对容城凸起基岩面的结构和下方地层的岩性以及断裂发育和分布情况进行了分析,获得如下认识:①容城凸起地区的基岩面上方覆盖近水平层状泥、砂岩地层,这些地层直接覆盖在具有高孔隙度的裂隙白云岩地层上,这种盖层结构对地热资源具有良好的保温效果;②容城凸起断裂发育,多条断裂切穿基岩面,其中由牛南断裂与容东断裂组成的大断裂带为地下热水的储存和热对流提供空间和通道;③由于热流的折射效应,容城凸起下方热流易从热导率低的地层向热导率高的地层流动,使得热流更好地向凸起区汇聚。

本文研究成果对加强地热资源的开发利用、深入理解地热能的形成过程以及研究地热富集规律有指导意义。但文章讨论未涉及热源以及深部地质结构,有待今后从这些方面对容城地热田热异常成因机制继续研究。

参考文献

邱楠生, 许威, 左银辉, .

渤海湾盆地中—新生代岩石圈热结构与热—流变学演化

[J]. 地学前缘, 2017, 24(3):13-26.

DOI:10.13745/j.esf.2017.03.002      [本文引用: 2]

文章主要利用中&mdash;新生代热史、地壳分层结构以及流变学参数,模拟计算渤海湾盆地中&mdash;新生代岩石圈热结构和热流变结构演化特征。结果表明,盆地由三叠纪&mdash;侏罗纪时期的&ldquo;冷幔热壳&rdquo;型岩石圈热结构转变为白垩纪至今的&ldquo;热幔冷壳&rdquo;型岩石圈热结构。从济阳坳陷岩石圈热流变结构演化特征来看,中生代早期上地壳上部、中地壳上部及上地幔顶部表现为厚的脆性层;早白垩世初期中地壳上部及上地幔顶部的脆性层完全转变为韧性层;晚白垩世开始,中地壳上部出现薄层的脆性层;古近纪早期中地壳上部脆性层变薄变浅;现今则除了发育上地壳上部、中地壳上部脆性层外,上地幔顶部开始在浅部发育薄的脆性层。中&mdash;新生代岩石圈总强度演化表明在早白垩世晚期和古近纪早期经历了两期减弱,中生代早期岩石圈总强度远大于中侏罗世之后的岩石圈总强度。岩石圈热流变结构和强度演化与华北克拉通破坏过程中岩石圈厚度的变化具有良好的对应关系,从侧面反映太平洋板块俯冲和回撤导致华北克拉通东部破坏的地球动力学过程。因此,岩石圈热流变结构可以为盆地形成、大陆边缘和造山带等的动力学演化过程研究提供科学依据。

Qiu N S, Xu W, Zuo Y H, et al.

Evolution of Meso-Cenozoic thermal structure and thermal-rheological structure of the lithosphere in the Bohai Bay Basin,eastern North China Craton

[J]. Earth Science Frontiers, 2017, 24(3):13-26.

[本文引用: 2]

王贵玲, 高俊, 张保建, .

雄安新区高阳低凸起区雾迷山组热储特征与高产能地热井参数研究

[J]. 地质学报, 2020, 94(7):1970-1980.

[本文引用: 1]

Wang G L, Gao J, Zhang B J, et al.

Study on the thermal storage characteristics of the Wumishan Formation and huge capacity geothermal well parameters in the Gaoyang low uplift area of Xiong’an New Area

[J]. Acta Geologica Sinica, 2020, 94(7):1970-1980.

[本文引用: 1]

陈墨香, 黄歌山, 张文仁, .

冀中牛驼镇凸起地温场的特点及地下热水的开发利用

[J]. 地质科学, 1982, 17(3):239-252.

[本文引用: 1]

Chen M X, Huang G S, Zhang W R, et al.

The temperature distribution pattern and the utilization of geothermal water at Niutuozhen basement protrusion of central Hebei Province

[J]. Chinese Journal of Geology, 1982, 17(3):239-252.

[本文引用: 1]

马峰, 王贵玲, 张薇, .

雄安新区容城地热田热储空间结构及资源潜力

[J]. 地质学报, 2020, 94(7):1981-1990.

[本文引用: 4]

Ma F, Wang G L, Zhang W, et al.

Structure of geothermal reservoirs and resource potential in the Rongcheng geothermal field in Xiong’an New Area

[J]. Acta Geologica Sinica, 2020, 94(7):1981-1990.

[本文引用: 4]

朱日祥, 徐义刚, 朱光, .

华北克拉通破坏

[J]. 中国科学:地球科学, 2012, 42(8):1135-1159.

[本文引用: 1]

Zhu R X, Xu Y G, Zhu G, et al.

Destruction of the North China Craton

[J]. Sci China Earth Sci, 2012, 42(8):1135-1159.

[本文引用: 1]

常健, 邱楠生, 赵贤正, .

渤海湾盆地冀中坳陷现今地热特征

[J]. 地球物理学报, 2016, 59(3):1003-1016.

DOI:10.6038/cjg20160322      [本文引用: 1]

渤海湾盆地冀中坳陷是我国最典型的潜山油气藏富集区.本文借助117口钻井地层测温资料和45块实测岩石热导率数据系统研究了冀中坳陷现今地温梯度、大地热流、热岩石圈厚度、岩石圈热结构等地热特征参数.研究表明,冀中坳陷0~3000 m统一深度现今地温梯度为20.8~41.0℃&#183;km-1,平均值为31.6℃&#183;km-1,比未校正值减小1~3℃&#183;km-1;现今大地热流介于48.7~79.7 mW&#183;m-2,平均值为59.2 mW&#183;m-2.平面上,冀中坳陷现今地温梯度和热流由西向东(从盆地边缘向内部)逐渐增大,并且凸起区地温梯度和热流相对较高,而凹陷区则偏低,与基底地形起伏具有很好的对应关系.同时,冀中坳陷腹部高热流凸起区广泛分布地热田.冀中坳陷现今热岩石圈厚度为98~109 km,其岩石圈热结构为一典型的"冷壳热幔"型.本研究不仅对冀中坳陷油气勘探与地热能开发具有重要的指导意义,而且为深部岩石圈研究(华北克拉通破坏科学问题)提供了新依据.

Chang J, Qiu N S, Zhao X Z, et al.

Present-day geothermal regime of the Jizhong depression in Bohai Bay Basin,East China

[J]. Chinese Journal of Geophysics, 2016, 59(3):1003-1016.

[本文引用: 1]

王凯, 张杰, 白大为, .

雄安新区地热地质模型探究:来自地球物理的证据

[J]. 中国地质, 2021, 48(5):1453-1468.

[本文引用: 1]

Wang K, Zhang J, Bai D W, et al.

Geothermal-geological model of Xiongan New Area:Evidence from geophysics

[J]. Geology in China, 2021, 48(5):1453-1468.

[本文引用: 1]

郭世炎, 李小军.

河北保定容城凸起地热田储层属性与资源潜力

[J]. 地质科学, 2013, 48(3):922-931.

[本文引用: 1]

Guo S Y, Li X J.

Reservoir stratum characterstics and geothermal resources potential of Rongcheng uplift geothermal field in Baoding,Hebei

[J]. Chinese Journal of Geology:Scientia Geologica Sinica, 2013, 48(3):922-931.

[本文引用: 1]

高锐, 周卉, 卢占武, .

深地震反射剖面揭露青藏高原陆—陆碰撞与地壳生长的深部过程

[J]. 地学前缘, 2022, 29(2):14-27.

DOI:10.13745/j.esf.sf.2021.7.13      [本文引用: 1]

印度板块与亚洲板块的碰撞使喜马拉雅-青藏高原隆升,地壳增厚并生长扩展。探测青藏高原深部结构,揭露两个大陆如何碰撞以及碰撞如何使大陆变形的过程,是对全球关切的科学奥秘的探索。深地震反射剖面探测是打开这个科学奥秘的最有效途径之一。二十多年来,运用这项高技术探测到青藏高原巨厚地壳的精细结构,攻克了难以得到下地壳和Moho面信息的技术瓶颈,揭露了陆-陆碰撞过程。本文在探测研究成果的基础上,从青藏高原南北-东西对比,再到高原腹地,系统地综述了青藏高原之下印度板块与亚洲板块碰撞-俯冲的深部行为。印度地壳在高原南缘俯冲在喜马拉雅造山带之下,亚洲板块的阿拉善地块岩石圈在北缘向祁连山下俯冲,祁连山地壳向外扩展,塔里木地块与高原西缘的西昆仑发生面对面的碰撞,在高原东缘发现龙日坝断裂(而不是龙门山断裂)是扬子板块的西缘边界,高原腹地Moho面厚度薄而平坦,岩石圈伸展垮塌。多条深反射剖面揭露了在雅鲁藏布江缝合带下印度板块与亚洲板块碰撞的行为,不仅沿雅鲁藏布江缝合带走向印度地壳俯冲行为存在东西变化,而且印度地壳向北行进到拉萨地体内部的位置也不同。在缝合带中部,研究显示印度地壳上地壳与下地壳拆离,上地壳向北仰冲,下地壳向北俯冲,并在俯冲过程中发生物质的回返与构造叠置,这导致印度地壳减薄,喜马拉雅地壳加厚。俯冲印度地壳前缘与亚洲地壳碰撞后沉入地幔,处于亚洲板块前缘的冈底斯岩基与特提斯喜马拉雅近于直立碰撞,冈底斯下地壳呈部分熔融状态,近乎透明的弱反射和局部出现的亮点反射以及近于平的Moho面都反映出亚洲板块南缘处于伸展构造环境。

Gao R, Zhou H, Lu Z W, et al.

Deep seismic reflection profile reveals the deep process of continent-continent collision on the Tibetan Plateau

[J]. Earth Science Frontiers, 2022, 29(2):14-27.

DOI:10.13745/j.esf.sf.2021.7.13      [本文引用: 1]

The collision between the Indian and Asian plates uplifted the Himalayan-Tibetan Plateau, thickening and expanding the crust. The deep structure of the plateau-a scientific mystery of global concern-can reveal how the two continents collide and how the collision deforms the continent. One of the most effective ways to unlock this scientific mystery is deep seismic reflection profiling. For more than 20 years using this technology we have detected the fine structure of the thick crust of the Tibetan Plateau, overcoming the technical bottleneck of accessing the lower crust and the Moho and revealing the intracontinental collision process. Based on the research results, this paper systematically summarizes the deep behavior of the collision and subduction between the Indian and Asian plates under the Tibetan Plateau, on topics ranging from the NS-EW correlation to the plateau’s hinterland. The deep behavior includes the underthrusting of the Indian crust beneath the Himalayan orogenic belt on the southern margin of the plateau, the subduction of the lithosphere of the Alashan block in the Asian plate under the Qilian mountain in the northern rim of the plateau, the outward crustal extension in the Qilian mountain range, and the face-to-face collision between the Tarim block and the West Kunlun in the northwestern margin of the plateau; in the eastern margin of the plateau it was discovered that the Longriba fault, not the Longmenshan fault, is the western border fault of the Yangtze plate. In the plateau’s hinterland the Moho lines are thin and flat and lithosphere extension collapses. Multiple deep reflection profiles reveal the collisional behavior of the Indian and Asian plates under the Yalungzangbo suture zone. The subduction of the Indian continental crust varies from east to west, and the crust travels northward to different locations inside the Lhasa terrain. In the central suture zone, the upper and lower parts of the Indian crust become detached as the upper crust overthrusts while the lower crust subducts northward. During the subduction, material return and structural superposition cause the thinning and thickening of the Indian and Himalayan crust, respectively. The leading edge of the subducting Indian crust collides with the Asian crust and sinks into the mantle, while a near-vertical collision occurs between the Gangdise at the leading edge of the Asian plate, and the Tethys. Together, the partial melting of the lower crust of the Gangdise, the seismic data that show near transparent weak reflection and localized bright spot reflection, and the near-flat Moho, all reflect the extensional tectonic environment of the southern margin of the Asian plate.

陈墨香, 汪集旸, 汪缉安, .

华北断陷盆地热场特征及其形成机制

[J]. 地质学报, 1990, 64(1):80-91.

[本文引用: 1]

Chen M X, Wang J Y, Wang J A, et al.

The characteristics of the geothermal field and its formation mechanism in the North China down-faulted basin

[J]. Acta Geological Sinica, 1990, 64(1):80-91.

[本文引用: 1]

唐博宁, 朱传庆, 邱楠生, .

雄安新区雾迷山组岩溶裂隙发育特征

[J]. 地质学报, 2020, 94(7):2002-2012.

[本文引用: 1]

Tang B N, Zhu C Q, Qiu N S, et al.

Characteristics of the Karst thermal reservoir in the Wumishan Formation in the Xiong’an New Area

[J]. Acta Geologica Sinica, 2020, 94(7):2002-2012.

[本文引用: 1]

孙冬胜, 刘池阳, 杨明慧, .

渤海湾盆地冀中坳陷中区中新生代复合伸展构造

[J]. 地质论评, 2004, 50(5):484-491.

[本文引用: 1]

Sun D S, Liu C Y, Yang M H, et al.

Study on complex extensional structures in the middle Jizhong depressionin the Bohai Bay basin

[J]. Geological Review, 2004, 50(5):484-491.

[本文引用: 1]

何登发, 单帅强, 张煜颖, .

雄安新区的三维地质结构:来自反射地震资料的约束

[J]. 中国科学:地球科学, 2018, 48(9):1207-1222.

[本文引用: 3]

He D F, Shan S Q, Zhang Y Y, et al.

3D geologic architecture of Xiongan New Area:Constraints from seismic reflection data

[J]. Scientia Sinica:Terrae, 2018, 48(9):1207-1222.

[本文引用: 3]

翟明国.

华北克拉通的形成演化与成矿作用

[J]. 矿床地质, 2010, 29(1):24-36.

[本文引用: 1]

Zhai M G.

Tectonic evolution and metallogenesis of North China Craton

[J]. Mineral Deposits, 2010, 29(1):24-36.

[本文引用: 1]

于福生, 漆家福, 王春英.

华北东部印支期构造变形研究

[J]. 中国矿业大学学报, 2002, 31(4):402-406.

[本文引用: 1]

Yu F S, Qi J F, Wang C Y.

Tectonic deformation of indosinian period in eastern part of North China

[J]. Journal of China University of Mining & Technology, 2002, 31(4):402-406.

[本文引用: 1]

安美建, 赵越, 冯梅, .

什么控制了华北克拉通东部在新近纪的构造活动?

[J]. 地学前缘, 2011, 18(3):121-140.

[本文引用: 1]

An M J, Zhao Y, Feng M, et al.

What resulted in new tectonic activities in the eastern North China Craton in the Neogene?

[J]. Earth Science Frontiers, 2011, 18(3):121-140.

[本文引用: 1]

朱日祥, 陈凌, 吴福元, .

华北克拉通破坏的时间、范围与机制

[J]. 中国科学:地球科学, 2011, 41(5):583-592.

[本文引用: 1]

Zhu R X, Chen L, Wu F Y, et al.

Timing,scale and mechanism of the destruction of the North China Craton

[J]. Sci China Earth Sci, 2011, 41(5):583-592.

[本文引用: 1]

商世杰, 丰成君, 谭成轩, .

雄安新区附近主要隐伏断裂第四纪活动性研究

[J]. 地球学报, 2019, 40(6):836-846.

[本文引用: 1]

Shang S J, Feng C J, Tan C X, et al.

Quaternary activity study of major buried faults near Xiongan new area

[J]. Acta Geoscientica Sinica, 2019, 40(6):836-846.

[本文引用: 1]

梁苏娟. 冀中坳陷晚新生代地质构造特征及其油气赋存[D]. 西安: 西北大学, 2001.

[本文引用: 1]

Liang S J. The Characteristics of Tectonics and Hydrocarbon Accumulation of Jizhong Depression in late Cenozoic era[D]. Xi’an: Northwest University, 2001.

[本文引用: 1]

索艳慧, 李三忠, 曹现志, .

中国东部中新生代反转构造及其记录的大洋板块俯冲过程

[J]. 地学前缘, 2017, 24(4):249-267.

DOI:10.13745/j.esf.yx.2017-3-17      [本文引用: 1]

中新生代期间,中国东部可以划分西部沉降带(包括东北、华北和中下扬子盆地群)和东部沉降带(包括东海和南海盆地群)。各盆地内部反转构造表现为东强西弱、自东向西迁移的规律,反映了构造反转的动力来自东部,证实了反转构造是古太平洋或太平洋板块的俯冲效应。但反转构造总体又具有自西部沉降带向东部沉降带迁移的特征:西部沉降带以中生代末期的挤压反转为主,而东部沉降带以渐新世末和中新世末的新生代反转最为典型。可见,虽然在同一俯冲板块作用下,仰冲板块的反转时间和迁移规律却存在差异性,这可能与反转构造所隶属的仰冲板内背景或俯冲板块自身行为在不同地区存在时空差异有关。为确定哪一种反转机制,本文搜集了古地磁条带和最新的板块重建方案、太平洋海山KAr年龄等数据,初步厘定出中新生代期间太平洋板块运动学行为;也对中新生代期间中国东部的东北、华北和华南的地质特征等开展对比研究。本文认为,伊泽奈崎板块早期为平板式俯冲,160 Ma左右板块拆沉并折返及岩石圈减薄,是西部沉降带的成盆机制;80 Ma左右伊泽奈崎板块向西的运动速度达到极值,对中国东部造成统一向西的挤压应力,是西部沉降带中生代的反转机制;大洋板块上发育的NW向破碎带的俯冲,可能加剧了这一期反转强度。60 Ma左右伊泽奈崎板块俯冲潜没于欧亚大陆之下及太平洋板块俯冲启动并后撤,这是东部沉降带的成盆机制;中新世中晚期,菲律宾海板块向欧亚板块的俯冲挤压导致了中国东部大陆左旋压扭的应力场,这是东部沉降带新生代构造的反转机制。

Suo Y H, Li S Z, Cao X Z, et al.

Mesozoic-Cenozoic inversion tectonics of East China and its implications for the subduction process of the oceanic plate

[J]. Earth Science Frontiers, 2017, 24(4):249-267.

[本文引用: 1]

胡秋韵, 高俊, 马峰, .

雄安新区容城凸起区地热可采资源量动态预测

[J]. 地质学报, 2020, 94(7):2013-2025.

[本文引用: 1]

Hu Q Y, Gao J, Ma F, et al.

Dynamic prediction of geothermal recoverable resources in the Rongcheng uplift area of the Xiong’an New Area

[J]. Acta Geologica Sinica, 2020, 94(7):2013-2025.

[本文引用: 1]

戴明刚, 雷海飞, 胡甲国, .

雄安新区顶面埋深在3500m以浅的中元古界热储可采地热资源量和开发参数评估

[J]. 地质学报, 2019, 93(11):2874-2888.

[本文引用: 2]

Dai M G, Lei H F, Hu J G, et al.

Evaluation of recoverable geothermal resources and development parameters of Mesoproterozoic thermal reservoir with the top surface depth of 3500 m and shallow in Xiongan New Area

[J]. Acta Geologica Sinica, 2019, 93(11):2874-2888.

[本文引用: 2]

黄元溢, 孟凡厚, 张占恩, .

雄安新区深部结构三维探测——雄安新区深地震反射剖面数据采集与常规处理项目采集施工总结

[R]. 中国石油集团东方地球物理勘探有限责任公司, 2018.

[本文引用: 2]

Huang Y Y, Meng F H, Zhang Z E, et al.

3 D detection of deep structure in Xiongan New Area:The collection and construction summary of deep seismic reflection profile data collection and routine processing project in Xiongan New Area

[R]. China National Petroleum Corporation. 2018.

[本文引用: 2]

单帅强. 太行山山前断层的构造几何学、运动学及其对渤海湾盆地发育的控制作用[D]. 北京: 中国地质大学(北京), 2018.

[本文引用: 2]

Shan S Q. Structural geometry and kinematics of the Taihang Mountain piedmont fault and its controlling on the development of the Bohai Bay basin[D]. Beijing: China University of Geosciences, 2018.

[本文引用: 2]

马岩, 张保建, 闫金凯, .

雄安新区深部储热构造探测研究与地热井优选技术

[J]. 地球学报, 2022, 43(5):699-710.

[本文引用: 1]

Ma Y, Zhang B J, Yan J K, et al.

Deep geothermal reservoir structure detection and geothermal well optimization technology in Xiongan new area

[J]. Acta Geoscientica Sinica, 2022, 43(5):699-710.

[本文引用: 1]

鲁锴, 鲍志东, 季汉成, .

雄安新区蓟县系雾迷山组岩溶热储特征、主控因素及有利区预测

[J]. 古地理学报, 2019, 21(6):885-900.

[本文引用: 1]

Lu K, Bao Z D, Ji H C, et al.

Characteristics,main controlling factors and favorable area prediction of karstic geothermal reservoirs of the Jixianian Wumishan Formation in Xiongan New Area

[J]. Journal of Palaeogeography, 2019, 21(6):885-900.

[本文引用: 1]

王贵玲, 李郡, 吴爱民, .

河北容城凸起区热储层新层系——高于庄组热储特征研究

[J]. 地球学报, 2018, 39(5):533-541.

[本文引用: 1]

Wang G L, Li J, Wu A M, et al.

A study of the thermal storage characteristics of Gaoyuzhuang Formation,A new layer system of thermal reservoir in Rongcheng uplift area,Hebei Province

[J]. Acta Geoscientica Sinica, 2018, 39(5):533-541.

[本文引用: 1]

熊亮萍, 高维安.

隆起与拗陷地区地温场的特点

[J]. 地球物理学报, 1982, 25(5):448-456.

[本文引用: 1]

Xiong L P, Gao W A.

Characteristics of geotherm in uplift and depression

[J]. Chinese Journal of Geophysics, 1982, 25(5):448-456.

[本文引用: 1]

熊亮平, 张菊明.

热流的折射和再分配的数学模拟

[J]. 地质科学, 1984, 19(4):445-454.

[本文引用: 1]

Xiong L P, Zhang J M.

Mathematical simulation of refract and redistribution of heat flow

[J]. Chinese Journal of Geology, 1984, 19(4):445-454.

[本文引用: 1]

/

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com , whtbjb@163.com