CN111963152A - 基于水平孔分层测温数据的地温梯度计算方法 - Google Patents

基于水平孔分层测温数据的地温梯度计算方法 Download PDF

Info

Publication number
CN111963152A
CN111963152A CN202010581211.4A CN202010581211A CN111963152A CN 111963152 A CN111963152 A CN 111963152A CN 202010581211 A CN202010581211 A CN 202010581211A CN 111963152 A CN111963152 A CN 111963152A
Authority
CN
China
Prior art keywords
horizontal
hole
temperature
temperature measurement
geothermal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010581211.4A
Other languages
English (en)
Other versions
CN111963152B (zh
Inventor
张晓宇
孟祥连
杜世回
黄勇
常帅鹏
马涛
李俊青
苗晓岐
王生仁
张文忠
王杜江
贾琦
李其帅
霍欣
杜庆者
张玉玺
周泉
梁树文
陈兴强
周福军
王栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway First Survey and Design Institute Group Ltd
Original Assignee
China Railway First Survey and Design Institute Group Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Railway First Survey and Design Institute Group Ltd filed Critical China Railway First Survey and Design Institute Group Ltd
Priority to CN202010581211.4A priority Critical patent/CN111963152B/zh
Publication of CN111963152A publication Critical patent/CN111963152A/zh
Application granted granted Critical
Publication of CN111963152B publication Critical patent/CN111963152B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/16Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements
    • G01K7/22Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements the element being a non-linear resistance, e.g. thermistor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/10Geothermal energy

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Nonlinear Science (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Agronomy & Crop Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明为一种基于水平孔分层测温数据的地温梯度计算方法,解决了现有技术中复杂艰险山区地热异常隧道工程水平及垂向地温梯度计算困难的问题,能够实现水平孔水平地温梯度和垂直地温梯度的计算,提高地温预测精度。本发明包括以下步骤:(1)由于极复杂艰险山区地形地面陡峻,深埋隧道工程部分垂直孔难以实施,为了查明地热异常区的水平向地温分布特征,选择合适的位置布置水平深孔勘探,(2)设计水平孔的孔深、方向、角度、隧道洞身位置、钻探测试要求;(3)开展水平孔的分层测温工作;(4)分析水平孔的分层测温数据变化特征,根据变化特征计算垂向地温梯度以及水平向地温梯度。

Description

基于水平孔分层测温数据的地温梯度计算方法
技术领域:
本发明属于深埋隧道高温热害评估领域,涉及一种地温梯度计算方法,尤其是涉及一种基于水平孔分层测温数据的地温梯度计算方法。
背景技术:
目前,地温梯度计算方法主要存在以下问题:
1、高原复杂艰险山区地热异常深埋隧道工程由于位于山岭地区,地形陡峻,垂直钻孔难以实施,开展水平孔的勘探工作能有效解决垂直钻孔难以实施的问题。
2、以往的地温计算中多是根据垂直孔地温曲线来拟合垂直孔的地温梯度,难以获取水平向地温梯度,对地温预测精度有一定影响。
3、目前还未开展过水平孔分层测温,对于水平孔分层测温获取的数据分析也是一个难点,尤其是受地热影响的水平孔,同时如何利用水平孔测温数据计算水平向及垂向地温梯度,也是行业领域空白。
发明内容:
本发明的目的在于提供一种基于水平孔分层测温数据的地温梯度计算方法,解决了现有技术中复杂艰险山区地热异常隧道工程水平及垂向地温梯度计算困难的问题,能够实现水平孔水平地温梯度和垂直地温梯度的计算,提高地温预测精度。
为实现上述目的,本发明采用的技术方案为:
一种基于水平孔分层测温数据的地温梯度计算方法,其特征在于:包括如下步骤:
S1,布置水平深孔勘探,勘探位置达到隧道设计标高附近或以下;
S2,设计水平孔的孔深、方向、角度以及钻探测试要求;
S3,进行水平孔的分层测温工作;
S4,分析水平孔的分层测温数据变化特征,根据变化特征计算垂向地温梯度以及水平向地温梯度。
S3步骤中,对水平深孔,采用无缆随钻测试方法,测温设备为热敏电阻,测温步骤为:将绳索钻杆下到预定位置,将探管连接好以后,向钻杆内注入高压水推送,一直将探头推出钻孔,然后提钻,将所有钻杆提出,再将改装电缆与绞车连接,上提测试,上提至孔口位置时进行深度归零;施钻过程中每100~200m进行一次测温,获取水平方向的地温数据及变化特征。
S4步骤中,绘制水平孔地温随深度变化曲线,判断地热异常区范围,提取单次水平孔测温的孔底地温,依据多次测试结果计算单孔的水平地温梯度和垂向地温梯度。
S4步骤中,整理出水平孔每次测温最深处的孔深、埋深、孔底标高、地温基础数据,根据式1计算水平向地温梯度;利用水平孔测温结果及埋深反算垂向地温梯度,根据式2确定测温点对应垂直高程处地表温度,根据式3结合测温点埋深反算水平孔的垂向地温梯度;
G水平=(Ti-T1)/(Hi-H1) 式1
式中:G水平:水平向地温梯度(℃/m);
Ti:第i次测温的孔底温度(℃);
Hi:第i次测温的孔底深度(m);
H1:根据地温曲线取200m;
T1:第一次分层测温200m的温度(℃);
TH=-0.0039Z0+21.331 式2
Z0:海拔高程(m)
G垂直=(Ti-Th)/(Hi-H0) 式3
G垂直:垂向地温梯度(℃/m)
H0:恒温层温度(℃)。
与现有技术相比,本发明具有的优点和效果如下:
1、本发明提出了水平孔分层测温,避免了受地热影响,钻孔成孔后一次测温,地温数据失真难以利用的难题。
2、本发明利用不同时间不同孔数获取的地温数据,从而能够实现水平孔水平地温梯度和垂直地温梯度的计算,提高地温预测精度。
附图说明:
图1是本发明提出的一种拉月隧道超深定向钻设计成果图(比例尺:水平1:1000;垂直1:200);
图2是本发明提出的一种水平深孔地温分层测试及监测结果。
具体实施方式:
下面结合具体实施方式对本发明进行详细的说明。这些实施例是用于说明本发明而不限于本发明的范围。
本发明为一种基于水平孔分层测温数据的地温梯度计算方法,具体包括以下步骤:
S1,由于极复杂艰险山区地形地面陡峻,深埋隧道工程部分垂直孔难以实施,为了查明地热异常区的水平向地温分布特征,选择合适的位置布置水平深孔勘探,勘探位置达到隧道设计标高附近或以下;
S2,设计水平孔的孔深、方向、角度以及钻探测试要求;
S3,开展水平孔的分层测温工作;
S4,分析水平孔的分层测温数据变化特征,根据变化特征计算垂向地温梯度以及水平向地温梯度。
S4步骤中,绘制水平孔地温随深度变化曲线,判断地热异常区范围,提取单次水平孔测温的孔底地温,依据多次测试结果计算单孔的水平地温梯度和垂向地温梯度。
S4步骤中,整理出水平孔每次测温最深处的孔深、埋深、孔底标高、地温基础数据,根据式1计算水平向地温梯度;利用水平孔测温结果及埋深反算垂向地温梯度,根据式2确定测温点对应垂直高程处地表温度,根据式3结合测温点埋深反算水平孔的垂向地温梯度;式1是本发明提出的一种利用水平孔分层测温计算水平地温梯度;式2是本发明提出的一种某高原艰险复杂山区地表恒温温度与海拔的关系式;式3是本发明提出的一种利用垂直孔分层测温计算垂直地温梯度。
G水平=(Ti-T1)/(Hi-H1) 式1
式中:G水平:水平向地温梯度(℃/m);
Ti:第i次测温的孔底温度(℃);
Hi:第i次测温的孔底深度(m);
H1:根据地温曲线取200m;
T1:第一次分层测温200m的温度(℃);
TH=-0.0039Z0+21.331 式2
Z0:海拔高程(m)
G垂直=(Ti-Th)/(Hi-H0) 式3
G垂直:垂向地温梯度(℃/m)
H0:恒温层温度(℃)。
实施例:
参见图1,本发明包括如下步骤:
S1,高原复杂艰险山区地热异常深埋隧道工程由于位于山岭地区,由于部分地区地形陡峻,垂直钻孔难以实施,开展了水平孔的勘探工作。
S2,参见图1,拉月隧道CSDXZ-1钻孔,位于林芝市巴宜区鲁朗镇拉月存北侧拉月2沟沟口。DK1192+615左1387m,水平角N48°W,下倾0.1°,孔口高程2606m,设计深度1400m。钻探目的查明钻孔所在位置地层岩性特征、层厚、完整程度、地质构造、富水性等工程地质条件和水文地质条件。
S3,对水平深孔,采用无缆随钻测试技术,测温设备为热敏电阻,主要测温步骤为:将绳索钻杆下到预定位置,将探管连接好以后,向钻杆内注入高压水推送,一直将探头推出钻孔,然后提钻,将所有钻杆提出,在将改装电缆与绞车连接,上提测试,上提至孔口位置时进行深度归零。施钻过程中每100~200m进行一次测温,获取水平方向的地温数据及变化特征,水平孔的分层测温数据变化巨大,存在同一位置后一次测温高于前一次的特点(图2)。
S4,绘制水平孔地温随深度变化曲线,判断地热异常区范围,提取单次水平孔测温的孔底地温,依据多次测试结果计算单孔的水平地温梯度和垂向地温梯度。
绘制水平孔地温随深度变化曲线,根据受地热影响区水平孔的测温结果,埋深越大温度明显增加,同时每次测出的地温均比上次的高(图2),通过拟合恒温层深度以下的曲线发现获取的水平地温梯度一次比一次高,分析认为这样的结论是不合理的,不符合自然规律。
这主要是由于水平孔内热交换强烈,深部的高温引起整个孔内地温发生变化,可以确定的是,每次分层测定的孔底温度是准确的,其余段落可能受到深部温度的影响,因此应整理处水平孔每次测温最深处的孔深、埋深、孔底标高、地温等基础数据,利用终孔温度减去地温线性变化起点温度除以水平距离来计算水平向地温梯度(式1),地温变化深度取200m,温度取第一次测得的200m深度温度。计算成果见表1,从计算结果看,米林断裂带(F51)范围为1100~950m,计算的水平地温梯度为4.44~6.20℃/100m,垂直地温梯度最大为14.7℃/100m,大于两侧完整基岩的水平地温梯度和垂直地温梯度,反映出米林断裂带为地热异常热源。
G水平=(Ti-T1)/(Hi-H1) 式1
公式中:G水平:水平向地温梯度(℃/m);
Ti:第i次测温的孔底温度(℃);
Hi:第i次测温的孔底深度(m);
H1:根据地温曲线取200m;
T1:第一次分层测温200m的温度(℃)。
利用水平孔测温结果及埋深反算垂向地温梯度,根据式2确定测温点对应垂直高程处地表温度,恒温层深度按照地形地貌及海拔确定为20~80m,结合测温点埋深反算水平孔的垂向地温梯度(式3)。
TH=-0.0039Z0+21.331 式2
Z0:海拔高程(m)
G垂直=(Ti-Th)/(Hi-H0) 式3
G垂直:垂向地温梯度(℃/m)
H0:恒温层温度(℃)。
表1利用水平孔测温数据计算水平及垂向地温梯度
Figure BDA0002552381420000071
本发明提出了水平孔分层测温,解决了受地热影响的钻孔成孔后一次测温,地温数据失真难以利用的难题。利用不同时间不同孔数获取的地温数据,从而能够实现水平孔水平地温梯度和垂直地温梯度的计算,提高地温预测精度。本发明方法成功实验应用于拉月隧道CSDXZ-1水平孔水平和垂直地温梯度计算、地热成因分析以及拉月隧道地温预测研究。
本发明的内容不限于实施例所列举,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。

Claims (4)

1.一种基于水平孔分层测温数据的地温梯度计算方法,其特征在于:包括如下步骤:
S1,布置水平深孔勘探,勘探位置达到隧道设计标高附近或以下;
S2,设计水平孔的孔深、方向、角度以及钻探测试要求;
S3,进行水平孔的分层测温工作;
S4,分析水平孔的分层测温数据变化特征,根据变化特征计算垂向地温梯度以及水平向地温梯度。
2.根据权利要求1所述的基于水平孔分层测温数据的地温梯度计算方法,其特征在于:
S3步骤中,对水平深孔,采用无缆随钻测试方法,测温设备为热敏电阻,测温步骤为:将绳索钻杆下到预定位置,将探管连接好以后,向钻杆内注入高压水推送,一直将探头推出钻孔,然后提钻,将所有钻杆提出,再将改装电缆与绞车连接,上提测试,上提至孔口位置时进行深度归零;施钻过程中每100~200m进行一次测温,获取水平方向的地温数据及变化特征。
3.根据权利要求1所述的基于水平孔分层测温数据的地温梯度计算方法,其特征在于:
S4步骤中,绘制水平孔地温随深度变化曲线,判断地热异常区范围,提取单次水平孔测温的孔底地温,依据多次测试结果计算单孔的水平地温梯度和垂向地温梯度。
4.根据权利要求3所述的基于水平孔分层测温数据的地温梯度计算方法,其特征在于:
S4步骤中,整理出水平孔每次测温最深处的孔深、埋深、孔底标高、地温基础数据,根据式1计算水平向地温梯度;利用水平孔测温结果及埋深反算垂向地温梯度,根据式2确定测温点对应垂直高程处地表温度,根据式3结合测温点埋深反算水平孔的垂向地温梯度;
G水平=(Ti-T1)/(Hi-H1) 式1
式中:G水平:水平向地温梯度(℃/m);
Ti:第i次测温的孔底温度(℃);
Hi:第i次测温的孔底深度(m);
H1:根据地温曲线取200m;
T1:第一次分层测温200m的温度(℃);
TH= -0.0039Z0 + 21.331 式2
Z0:海拔高程(m)
G垂直=(Ti-Th)/(Hi-H0) 式3
G垂直:垂向地温梯度(℃/m)
H0:恒温层温度(℃)。
CN202010581211.4A 2020-06-23 2020-06-23 基于水平孔分层测温数据的地温梯度计算方法 Active CN111963152B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010581211.4A CN111963152B (zh) 2020-06-23 2020-06-23 基于水平孔分层测温数据的地温梯度计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010581211.4A CN111963152B (zh) 2020-06-23 2020-06-23 基于水平孔分层测温数据的地温梯度计算方法

Publications (2)

Publication Number Publication Date
CN111963152A true CN111963152A (zh) 2020-11-20
CN111963152B CN111963152B (zh) 2023-12-08

Family

ID=73360745

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010581211.4A Active CN111963152B (zh) 2020-06-23 2020-06-23 基于水平孔分层测温数据的地温梯度计算方法

Country Status (1)

Country Link
CN (1) CN111963152B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114382540A (zh) * 2021-12-22 2022-04-22 中铁二院工程集团有限责任公司 一种确定高地温区铁路空间线位的方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103132990A (zh) * 2011-11-25 2013-06-05 普拉德研究及开发股份有限公司 井下温度分布的动态预测
CN104500045A (zh) * 2014-12-15 2015-04-08 中铁第四勘察设计院集团有限公司 一种深埋洞室内钻孔注浆的地温测试方法
CN105652342A (zh) * 2016-01-13 2016-06-08 中国石油化工股份有限公司 一种基于地层不整合面的分段式地温梯度拟合方法
CN108875995A (zh) * 2018-03-19 2018-11-23 中铁隧道局集团有限公司 一种未开挖条件下高地温隧道原岩温度的预测方法
CN110487443A (zh) * 2019-09-26 2019-11-22 中铁隧道局集团有限公司 一种高地温状态下隧道开挖检测方法
CN209724420U (zh) * 2019-04-12 2019-12-03 中国矿业大学 高地温矿井地温测试系统
CN210685972U (zh) * 2019-08-21 2020-06-05 中煤地华盛水文地质勘察有限公司 一种用于l型地热井的测温装置
US20200190968A1 (en) * 2018-12-12 2020-06-18 Exxonmobil Upstream Research Company Method for Zonal Injection Profiling and Extraction of Hydrocarbons in Reservoirs

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103132990A (zh) * 2011-11-25 2013-06-05 普拉德研究及开发股份有限公司 井下温度分布的动态预测
CN104500045A (zh) * 2014-12-15 2015-04-08 中铁第四勘察设计院集团有限公司 一种深埋洞室内钻孔注浆的地温测试方法
CN105652342A (zh) * 2016-01-13 2016-06-08 中国石油化工股份有限公司 一种基于地层不整合面的分段式地温梯度拟合方法
CN108875995A (zh) * 2018-03-19 2018-11-23 中铁隧道局集团有限公司 一种未开挖条件下高地温隧道原岩温度的预测方法
US20200190968A1 (en) * 2018-12-12 2020-06-18 Exxonmobil Upstream Research Company Method for Zonal Injection Profiling and Extraction of Hydrocarbons in Reservoirs
CN209724420U (zh) * 2019-04-12 2019-12-03 中国矿业大学 高地温矿井地温测试系统
CN210685972U (zh) * 2019-08-21 2020-06-05 中煤地华盛水文地质勘察有限公司 一种用于l型地热井的测温装置
CN110487443A (zh) * 2019-09-26 2019-11-22 中铁隧道局集团有限公司 一种高地温状态下隧道开挖检测方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114382540A (zh) * 2021-12-22 2022-04-22 中铁二院工程集团有限责任公司 一种确定高地温区铁路空间线位的方法

Also Published As

Publication number Publication date
CN111963152B (zh) 2023-12-08

Similar Documents

Publication Publication Date Title
Kang et al. In-situ stress measurements and stress distribution characteristics in underground coal mines in China
Zhang et al. Estimation of in situ stress along deep tunnels buried in complex geological conditions
CN109828316B (zh) 一种钙结岩型铀矿找矿勘查方法
CN106401574B (zh) 一种钻前高温地热井地层压力的预测方法
CN102175832A (zh) 一种确定典型储层最佳饱和度计算模型的方法
CN110595887B (zh) 一种浅埋硬岩大型地下洞室成洞条件量化评价方法
CN103410516A (zh) 小直径隧洞全断面开挖围岩变形预警工法
Shuai et al. Key oil accumulation periods of ultra-deep fault-controlled oil reservoir in northern Tarim Basin, NW China
CN108169456A (zh) 一种湿陷性黄土预估大压力评价方法
CN104459823A (zh) 一种综合物探测试海水侵入界面方法
CN111963152B (zh) 基于水平孔分层测温数据的地温梯度计算方法
CN111912540B (zh) 高原复杂艰险山区地热异常隧道工程地温测试及预测方法
Duchkov Review of Siberian heat flow data
CN111089662A (zh) 一种用于浅层地温测量的方法
CN108343431B (zh) 一种基岩裂隙地下水系统流网探究方法
Xing et al. Evaluation of the modification potential for carbonate geothermal rocks by considering multiple mechanics properties and acid dissolution in the Jizhong Depression, China
CN111768104B (zh) 单斜砂泥岩互层型坝基岩性变化程度定量评价方法
Wei et al. Application of geological model dynamic updating method in Guanjiapu Oilfield, Dagang
Obara et al. Induced stress measurement around a tunnel by the compact conical-ended borehole overcoring (CCBO) technique
CN114510808A (zh) 一种循缝找洞酸压数值模拟缝洞属性参数精细表征方法
Zheng et al. Normalization processing and research of multi-field data in reservoir seepage investigation
CN111929742B (zh) 高原复杂山区地热异常深埋隧道工程高温热水判别方法
Ai et al. Study on the Method for Measuring Stress on Rock Mass Excavation Surface under Extremely High Stress Conditions
CN104615900A (zh) 一种地下水排入湖泊营养盐的定量计算方法
CN104881536B (zh) 一种岩体结构面吻合度系数测量方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant