利用钻孔温度梯度重建过去地表温度变化研究

摘要

对过去气候变化的研究有助于我们理解现在气候变化和预测未来气候。作为地球“第三极”的青藏高原，具有独特的地理位置和热力系统，并对气候变化非常敏感。了解青藏高原过去气候变化历史、变化特征，对于准确评估目前气候变化与预测未来青藏高原气候变化提供了重要依据。然而青藏高原的气象观测资料观测时间短，观测站点少且覆盖不均，因此仅靠气象观测不足以让我们了解青藏高原长时间尺度的气候变化。过去地表温度变化的气候信息可通过分析近期观测的钻孔温度剖面进行重建。青藏高原作为多年冻土大区，多年冻土内部土壤冻结，地温热传递主要以热传导方式进行，适宜利用钻孔温度方法的开展。钻孔温度方法较其他代用指标方法是基于地温剖面和地表温度变化的物理联系之上的，具有更强的物理意义。本文利用青藏高原多个钻孔温度剖面和地热梯度，重建青藏高原不同地区过去地表温度变化历史。
　　地温垂直温度剖面受地中热流和地表温度变化影响。地球内部热流通过地热梯度影响地温形成稳态温度，在此基础上地表温度变化以热传导方式向地下传播使稳态地温发生扰动而产生偏离，偏离稳态温度的这部分地温作为瞬时温度记录了过去地表温度变化的信息。钻孔温度方法利用近期观测到的钻孔温度剖面，根据地温梯度分离出稳态温度，进而利用一维热传导模型分析瞬时温度剖面，重建过去地表温度变化。由于年际地表温度变化随深度明显衰减，而十年或更长时期的温度变化（“古气候信号”）向多年冻土更深处传播，多年冻土温度分布是气候变化和地表能量平衡长期变化的敏感指示。陆地作为天然的气候低通滤波器，使得多年冻土区的地温剖面可用来重建过去地表温度低频变化趋势。
　　在研究钻孔温度重建过去地表温度变化的反问题之前，首先需要理解正向问题，地温如何响应地表温度变化。对于简单一次气候变化情形，根据一维热传导模型，对于不同类型地表温度变化地下温度场具有显式的解析解。对于多年冻土区域热传导且考虑相变的复杂过程，仅能通过数值方法求解。本文利用控制体积方法数值模拟多年冻土地温相变问题。控制体积方法在离散程度上介于有限差分和有限元方法之间，具有更直接的物理解释。计算时考虑未冻水含量及相变潜热，并随时间重新计算各深度热物理参数进而准确计算地温变化。通过地表周期变化情形阐明控制体积方法数值模拟多年冻土相变问题的过程及未冻水含量对热参数、地温的影响。
　　基于对正问题的理解，本文对钻孔温度重建过去地表温度变化的反问题提出改进的Tikhonov方法。基于过去地表温度变化和近期钻孔地温剖面间的物理联系，在对问题参数化后我们利用Tikhonov正则化方法来重建过去地表气候变化。此方法是基于奇异值分解（SVD）方法之上的，具有相同的参数化，都将问题转化为求解不适定的矩阵方程。本文利用两个数值例子模拟地表升温和复杂气候事件来验证方法有效性以及与奇异值分解方法相比的改进效果。由于钻孔温度观测具有无法避免的观测误差，我们对模拟的地温剖面添加随机扰动误差来模拟观测误差。利用扰动的地温剖面重建地表温度变化并与假设的地温剖面比较从而验证方法有效性。两个实验例子中重建的地表温度及相关不确定性分析表明Tikhonov方法可较好地重建地表温度，并且改进方法可成功压制噪音导致的不稳定性得到更平滑的地表温度变化。通过比较地表温度误差可选出更适合 Tikhonov正则化的正则化参数选取方法。此外，本文利用 Tikhonov方法分析了钻孔温度方法的求解能力，可更好地理解重建的地表温度。
　　钻孔温度重建地表温度变化的结果是依赖于所选取的反演方法，我们基于热传导方程反边界值问题给出创新的基本解方法。基本解方法根据热传导方程的基本解将问题参数化之后转化为待定线性系统，由于反问题的不适定性，此方程组受钻孔温度观测误差影响无法直接求解。利用 Tikhonov正则化和广义交叉核实方法选取正则化参数求解待定参数，进而同时重建地表温度变化和地表温度热流变化。数值模拟例子表明基本解方法是可行且稳定的，并且对模拟钻孔温度剖面添加不同水平随机误差扰动后仍能有效降低误差扰动带来的不适定性，得到精确的地表温度变化。
　　利用钻孔温度重建地表温度变化与其他地球物理反问题相同，最大的求解难点在于观测误差导致的结果不稳定性。由于不同方法采用不同的参数化和优化方法，因此钻孔温度问题的结果依赖于所选取的方法。本文综合比较应用较广泛的泛函空间反演（FSI）、奇异值分解方法（SVD）、改进的Tikhonov方法以及创新的基本解方法（MFS）。通过五类不同类型地表温度变化的模拟例子来比较各方法数值结果：（1）阶梯变化；（2）线性升温；（3）光滑线性升温；（4）周期变化；（5）复杂周期变化，并在模拟例子中添加不同水平模拟观测误差。重建的地表温度结果比较表明，在钻孔温度剖面具有较小误差扰动下，所有方法均能给出较精确的地表温度变化重建结果。尽管四种方法具有不同的参数化方式和正则化方法选取，重建地表温度变化具有相似结果，仅在气候时间和温度幅度上有细微差别。方法的有效性是依赖于地表温度变化类型的。基本解方法更适用于重建周期变化和复杂周期变化信号。对于其他类型地表温度变化，Tikhonov方法在较小钻孔温度误差0.001℃和0.01℃情形下，结果最精确。泛函空间反演方法较在钻孔温度剖面误差较大时仍能重建地表温度变化趋势，但对初始地表温度重建较其他方法相比有较大误差。并且泛函空间反演在重建近期地表温度变化时具有更高分辨率，更精确。
　　基于钻孔温度方法研究，本文根据青藏高原不同地区钻孔温度剖面对各研究点进行过去地表温度变化的单点重建研究，利用钻孔温度方法反演得到不同时间区间古气候信息：1）黑河上游100米钻孔PT1钻孔由于气候变暖导致进入多年冻土的长期净热流约为0.014 Wm-2，深处稳态热流约0.0247 Wm-2。PT1钻孔1952年至2012年地表温度由-2.7℃线性升高约0.5至0.65℃；2）黑河上游150米钻孔PT9钻孔地热梯度为2.25℃/100m，1895年至2015年地表温度由-2.3℃升温至-1.5℃；3）奇异值分解方法和Tikhonov方法根据五道梁120米钻孔温度剖面重建地表温度结果表明在过去1930年至2013年间地表温度升温1.8(±0.2)℃，且剧烈升温过程开始于1980年代。五道梁气象观测站的气温观测结果验证了Tikhonov方法重建2008年至2012年间的地表温度波动，且在时间重合阶段气温和重建的地表温度具有相同趋势。4）根据昆仑山钻孔220米钻孔温度剖面，奇异值分解方法和 Tikhonov方法重建地表温度结果表明，1700年至2013年地表温度由-6.5(±0.8)℃升高至-2.8(±0.2)℃。两方法重建的地表温度变化具有相同趋势，具体升温时间和幅度略有差别。根据五道梁观测站的气温观测数据对比表明，Tikhonov方法重建的地表温度更可靠。5）柴达木盆地7个钻孔（最大深度220米至400米）温度反演表明此区域过去514年地表温度由6.1℃升高了1.2℃(-0.11~2.21℃)，并表现出1500年至1900年间的小冰期寒冷信号。最冷时期发生在1780至1790年间，当时的地表温度为5.4℃。在19世纪和20世纪间，重建的地表温度具有升温趋势，且在20世纪末达到最高值，随后开始降温。重建的地表温度变化幅度已由EdGCM模式模拟的地表平均气温所验证，细节温度特征得到代用指标结果验证。
　　基于钻孔温度方法研究，本文给出钻孔温度方法的改进方法和创新方法，并进行方法比较。本论文根据青藏高原中部不同地区钻孔温度剖面进行过去地表温度变化的单点重建研究，并利用地温梯度反演得到各钻孔位置的过去地表温度变化信息。

目录

封面

中文摘要

英文摘要

目录

第一章 绪论

1.1 研究背景及研究意义

1.2 研究内容与论文结构

第二章 钻孔温度方法原理及进展

2.1 模型描述

2.2 钻孔温度方法

2.3 影响因素

2.4 研究现状及发展趋势

2.5 方法局限性及存在的问题

第三章 野外钻孔资料

3.1 黑河钻孔

3.2 五道梁钻孔

3.3 昆仑山钻孔

3.4 柴达木盆地

第四章 钻孔温度方法改进及创新

4.1 方法改进——Tikhonov方法

4.2 方法创新——基本解方法

4.3 方法比较

第五章 应用钻孔温度梯度重建地表温度变化历史

5.1 黑河上游PT1 100米深孔

5.2 黑河上游PT9 150米钻孔

5.3 五道梁钻孔

5.4 昆仑山钻孔

5.5 柴达木盆地

5.6 地表温度重建小结

第六章 结论与展望

6.1 主要结论

6.2 存在的问题与研究展望

参考文献

在学期间的研究成果

致谢