Theory and Method for Urban–Rural Integration, Evaluation and Cooperative Detection of Ground and Underground Space: Example from the Urban Geological Survey of Guanzhong Plain
-
摘要:
城乡融合、地上与地下空间协同探测、评价、规划是节约集约利用空间与自然资源,增强城市韧性,促进人与自然和谐共生现代化建设的重要途径。针对中国新型城镇化与乡村振兴协同发展面临着“多规合一”的国土空间规划、地上与地下空间协调集约利用等新的重大科技问题,以关中平原城市群为研究对象,基于地球系统科学和地球关键带理论,汲取系统论、协同论、信息论等理论,提出了城乡融合、地上与地下空间协同探测、综合评价、统筹规划的新时期城市地质理论,以及城市强干扰环境下地下空间精细化探测和建模关键技术、基于地下水位控制的地上与地下空间风险防控技术、“取热不取水”的地下空间地热能绿色开发新技术,并在多个城市得到成功应用。
Abstract:The collaborative detection and evaluation planning of ground and underground space is an important way to enhance the comprehensive carrying capacity of the city and promote the high–quality development of the city. In view of the new major scientific and technological problems such as "multi–plan integration" territorial spatial planning, coordinated and intensive utilization of ground and underground space aimed at the coordinated development of new urbanization and rural revitalization in China, taking Guanzhong plain urban agglomeration as the research area, based on the theories of the earth system science, the earth key belt, system, synergy and information, this paper puts forward the urban geological theory in the new period of urban–rural integration, coordinated detection, comprehensive evaluation and overall planning of ground and underground space; breaks through the key technologies of fine detection and modeling of underground space under the strong urban interference environment; proposes the risk prevention and control technology of ground and underground space based on groundwater; reveal the heat conduction law of geothermal fluid–exchanger in underground space, and raise the new green development technology of geothermal energy in underground space of "taking heat but not taking water". Which will provide the technique support for the urban planning, construction and operation management, land space planning, natural resource management, etc. And the theroies and the technologies have been successfully applied in many cities.
-
Keywords:
- urban geology /
- underground space /
- new urbanization /
- rural revitalization /
- Guanzhong plain
-
多年冻土是影响寒区生态环境最活跃的因素。在全球变暖的背景下,世界范围内的多年冻土发生了不同程度的退化,主要表现为多年冻土面积减少、厚度减薄、活动层厚度增加、多年冻土温度升高等。青藏高原是世界上中低纬度面积最大、海拔最高的多年冻土区(Zhao et al., 2020;吴青柏等, 2021),对气候变化更敏感(Peng et al., 2018),变暖速率是全球平均水平的两倍(Zhang et al., 2019; You et al., 2021),多年冻土面积减少速率约为1.15×104 km2/a(Gao et al., 2023)。多年冻土会直接或间接影响区域生态水文过程、碳循环乃至工程建设的稳定性,活动层热融和冻胀的特性也会改变地表形态,形成冻胀丘、热融滑塌、冻土沼泽等不良地质现象。为此许多国际科学机构将多年冻土作为重要监测和研究内容,包括极地系统科学计划(ARCSS)、环极地活动层监测计划(CALM)、多年冻土热状态计划(TSP)等。
传统的监测手段依赖大量的野外作业,然而受限于恶劣的自然环境,野外作业具有成本高、覆盖范围小等缺点,往往集中在青藏公路沿线、黑河流域等重点区域。卫星遥感和大地测量技术的迅速发展,为大范围多年冻土观测提供了新手段。合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar, InSAR)是一种基于微波观测地表变形的技术,具有全天候、精度高、范围广等优势(冯旻譞等, 2023),可以很好地观测冻融过程由水冰相变引起的地表周期性变形(张廷军, 2012; Wang et al., 2020),在区域尺度多年冻土调查中的应用日趋成熟。Daout等(2017)利用InSAR技术反演了2003至2011年青藏高原区域大范围的季节形变量,并对季节性形变延迟进行了探讨;Zhang等(2019)、王志伟等(2021)探讨了不同下垫面以及土壤孔隙度等因素对地表变形分布特征的影响;范鹏等(2021)、逯中香等(2022)通过InSAR技术研究了地表冻胀融沉作用对青藏工程设施稳定性的影响;Liu等(2012)、Wang等(2018)基于水质量守恒建立了InSAR地表变形与活动层厚度的函数关系来反演活动层厚度。目前采用InSAR技术对多年冻土的研究主要聚焦多年冻土地表变形空间分布特征、对工程稳定性的影响、多年冻土参数反演等方面,针对地表冻胀融沉与活动层水热过程定量关系的研究较少。为此,笔者采用时序InSAR技术对祁连山地区野牛沟地表变形进行5年的连续观测,分析野牛沟地表变形的成因;结合气象站观测数据,通过相关分析等方法研究地表变形与土壤水热动态变化的关系,对研究长时序、大范围多年冻土变化具有参考意义。
1. 研究区概况
野牛沟位于青海省祁连县,地处青藏高原东北部祁连山地区(图1),地形地貌主要由新生代构造活动形成,野牛沟南北两侧山脉呈NW—SE走向,地貌类型复杂,第四季冰川和冰碛堆积地貌大量发育(郭怀军, 2017)。河谷内部是黑河上游西支冲积平原,地势较为平坦。野牛沟受青藏高原和大陆性气候影响,为高原(亚)高寒带半干旱气候,年降雨量约400 mm,多年平均气温约为−4 ℃(王庆锋等, 2016)。研究区土壤质地粗糙,有机物含量少,主要为砂质土、黏土、亚粘土(陈玉兴等, 2019),植被类型为高寒草甸、高寒沼泽草甸(Cao et al., 2018; 马素辉等, 2018)。野牛沟属阿尔金–祁连山高寒带高山冻土区,多年冻土十分发育,多年冻土下界约为
3650 ~3700 m,活动层厚度为1~4.9 m,平均厚度约为2.7 m(Cao et al., 2017; Han et al., 2022)。2. 数据与方法
2.1 数据
(1)SAR数据
SAR数据采用Sentinel-1A卫星降轨单视复数(single look complex, SLC)影像,时间跨度为2017年10月至2022年10月,共152期,设置最大时间基线为60 d,数据时空基线分布及干涉像对组合见图2。
(2)土壤温湿度数据
通过国家青藏高原科学数据中心(https://doi.org/10.11888/Atmos.tpdc.300628)黑河流域地表过程综合观测网,收集大沙龙气象站2017~2022年不同深度的土壤温度、土壤含水量数据(Liu et al., 2018; Che et al., 2019; Liu et al., 2021),研究活动层水热过程。
2.2 小基线子集
小基线子集(small baseline subsets, SBAS)-InSAR方法采用多主影像组合,通过设置较短的时间和空间基线生成干涉像对,保证干涉图具有较高的相干性。与永久散射体(persistent scatters, PS)-InSAR方法相比,SBAS方法更关注分布式特征,在缺少高相干点的区域依然可以获取较可靠的结果,更适用于自然特征(张林梵, 2023)。野牛沟地表主要为高寒草甸等自然地表,人工建筑等高相干点十分稀少,因此笔者选择SBAS-InSAR技术作为地表变形监测的技术手段。
假设有N+1幅SAR影像,通过设置合适的时空基线得到M幅差分干涉图,M满足:
$$ (N+1)/2\leqslant M\leqslant N(N+1)/2 $$ (1) 用$ {v}_{i} $表示相邻某一时段的平均相位速率,第i幅和第j幅干涉对差分相位$ \Delta {\varnothing }_{i,j} $可以表示为该时段平均相位速率在时间间隔内的和:
$$ \left[\begin{array}{c}\Delta {\varnothing }_{12}\\ \Delta {\varnothing }_{13}\\ \Delta {\varnothing }_{23}\\ \cdots \\ \cdots \\ \Delta {\varnothing }_{mj}\end{array}\right]=\left[\begin{array}{ccccc}{t}_{2}-{t}_{1}& 0& 0& \cdots & 0\\ {t}_{2}-{t}_{1}& {t}_{3}-{t}_{2}& 0& \dots & 0\\ 0& {t}_{3}-{t}_{2}& 0& \cdots & 0\\ \cdots & \cdots & \cdots & \cdots & \cdots \\ \cdots & \cdots & \cdots & \cdots & \cdots \\ 0& 0& 0& \cdots & {t}_{j}-{t}_{j-1}\end{array}\right]\left[\begin{array}{c}{v}_{1}\\ {v}_{2}\\ {v}_{3}\\ \cdots \\ \cdots \\ {v}_{M-1}\end{array}\right] $$ (2) 式(2)系数矩阵可用B表示。在解算地表形变时,对于小基线子集内部,N≥M可以使用最小二乘法求解,当各子集间N<M,矩阵B容易产生秩亏,因此采用奇异值分解方法(SVD)得到速率的最小范数解,从而获得某一时间段的形变相位和地形等噪声相位。
3. 野牛沟地表变形特征与成因
本研究提及的InSAR地表变形方向均为降轨LOS(line of sight)向,InSAR解算野牛沟地表变形见图3,受地形地貌、积雪等因素影响,南北两侧山脉失相干较严重。河谷内部地势较平坦,植被类型主要为高寒草甸,地表变形主要受冻融作用影响,年均形变速率较小,为−10~10 mm/a;南北两侧山脉部分沟谷、边坡受融水侵蚀作用,年均形变速率较大,最大可达−162 mm/a。
3.1 冻融作用
为探究野牛沟地区冻融过程,采用季节–趋势分解法(Seasonal-Trend decomposition using LOESS, STL)对A区域2021年时序变形进行季节性分解见图4。A区域位于河谷内部山前平原(图3),植被类型为高寒草甸(Wang et al., 2017),西、北区域海拔位于多年冻土下界以上,东南区域位于多年冻土下界以下。该区域地表变形表现出明显的季节性特征,4~5月随着温度升高地表开始出现融沉,季节变形量逐渐增大并持续至9~10月,东南部季节冻土区融沉量更大。分别在多年冻土区与季节冻土区选择典型点PT1和PT2进行长时间序列地表变形分析如图5所示,季节冻土区与多年冻土区地表变形呈现相同的季节性特征,受局地因素影响,地表冻融变形量级有较大差异,年际地表冻胀、融沉幅度约为10~20 mm。
3.2 水力侵蚀作用
水力侵蚀主要表现为冰川积雪融水对坡体的侵蚀和河流对边坡的侵蚀作用。由于靠近冰川积雪区域失相干较为严重,笔者仅对河流边坡侵蚀变形进行分析。图6所示为图3中B、C两区域河流侵蚀边坡变形空间特征。B区域两条河流边坡均出现明显变形且只发生在一侧,形变速率约为50 mm/a(图6c);C区域两侧坡体均有变形分布,呈现漏斗状特征,漏斗中心变形速率可达80~100 mm/a,两侧边坡形变信号一正一负是坡向与雷达视角原因所致(王燕燕等, 2019)(图6d)。受青藏高原和大陆性气候影响,野牛沟属高原(亚)高寒带半干旱气候,年降水量少且集中,河流补给主要来源于冰川和积雪融水(Cheng et al., 2014)。夏季温度升高时,冰雪融水补给河流,河流流量增加对地表侵蚀作用增强导致边坡形变速率较大;而冬季冰雪以固态形式存在,河流流量减少,对边坡的侵蚀作用减弱,地表基本处于平稳状态。因此河流边坡地表变形曲线呈现明显的季节性特征(图7)。其中C区域地表变形在冬季出现了小幅度的反弹,可能与地表水位较浅、坡体局部区域含水量较高有关,降温时地表首先冻结阻碍了地下水的流动,使得地下水在此区域进一步富集冻结,导致地表出现冻胀抬升(吴吉春等, 2007)。
4. 地表变形与土壤水热过程
4.1 土壤冻融过程
大沙龙站位于祁连县西侧沙龙滩地区,海拔
3739 m,下垫面为沼泽化高寒草甸。大沙龙站在0 cm、−4 cm、−10 cm、−20 cm、−40 cm、−80 cm、−120 cm、−160 cm处记录土壤温度和土壤含水量,观测间隔10 min,数据周期为2017年1月1日至2022年12月31日。为方便对比,笔者根据SAR成像时间对气象观测数据进行了筛选。由于浅层土壤温度的日变化十分活跃,难以反映季节特征,仅对−20 cm深度及以下的土壤温度进行分析。图8为2017~2021年大沙龙站不同深度土壤温度随时间的变化过程(2022年数据缺失较为严重)。以2017年为例,从季节冻融过程来看,5月6日深度−20 cm处开始融化,随着时间的推移融化锋面逐渐下移,该过程土体处于吸热状态,融化过程是自上而下的,6月25日融化深度超过−160 cm,至此观测深度的土体全部处于负温状态,历时50 d。11月1日深度−20 cm处开始冻结,在−160 cm深度冻结始日则为11月16日,冻结过程是自上而下和自下而上的双向过程,土体处于上下温度低、中部温度高的状态(王庆锋等, 2016; 马晶晶等, 2022),冻结锋面交汇于深度约−120 cm处。值得注意的是,由于温度向下传输需要一定的过程(Li et al., 2015),因此地表冻融始日较−20 cm处应有所提前。
其他年份冻融过程与2017年类似,差异主要体现在冻融过程的时间节点和不同深度的温度分布。从历年数据对比来看,以−4 ℃等温线为例,其下探深度逐年增加、历时变长。从冻结锋面来看,每年交汇深度基本不变,而交汇时间逐渐提前,分别为2017年11月30日、2018年11月20日、2019年11月15日、2020年11月15日、2021年11月15日。以上数据表明,大沙龙站具有较小的降温趋势,与Sun等(2022)在野牛沟的研究结论一致。
4.2 地表变形与土壤水热过程
大沙龙站InSAR地表变形与土壤温湿度分层数据随时间变化的关系见图9,地表融沉变形始日大约在4月底至5月初,冻胀变形始日在10月初,较深度−20 cm处0 ℃等温线出现始日有所提前,原因已在前文述及;各层土壤温度在2017和2021年略高,2018年与2019年略低,土壤含水量则呈现明显的降低趋势。地表冻胀、融沉幅度与土壤含水量变化过程一致,呈逐年降低的趋势。从表1地表变形与土壤水热过程的相关性来看,地表变形与不同深度的土壤温度均表现出较高的负相关性,在−40 cm处相关性最显著(R=−0.513, P<0.001);与不同深度的土壤含水量相关性差异较大,在−10~−40 cm处相关性最显著(R=−0.415, P<0.001)。
表 1 地表形变与土壤温湿度相关系数Table 1. Correlation coefficient between surface deformation with soil temperature and humidity, respectively深度(cm) 0 4 10 20 40 80 120 160 土壤温度 −0.479 −0.455 −0.493 −0.513 −0.522 −0.517 −0.492 −0.471 土壤含水量 / −0.293 −0.415 −0.396 −0.409 −0.240 −0.209 −0.186 土壤水分在热量传输过程中发挥着重要作用(杜冉等, 2022),季节变形幅度与土壤含水量有密不可分的关系。假设地表冻胀、融沉主要是水冰相变体积差异导致的,理想状态下季节融沉量$ {d}_{season} $与单位体积水$ {z}_{water} $的关系可以表示为:
$$ \Delta {d}_{season}=\frac{{\rho }_{w}-{\rho }_{i}}{{\rho }_{i}}{z}_{water}\approx \text{0.09}{z}_{water}dh $$ (3) 其中:$ {\rho }_{w} $和$ {\rho }_{i} $分别为水和冰的密度,即100 mm的水柱冻结为冰将产生约9 mm的体积变化(Chen et al., 2020)。考虑到土壤孔隙度等因素,并非所有水分的冻融过程均会对地表冻胀融沉产生贡献,因此采用极差表示年内水分完全冻结和完全融化的量以及由此产生的冻胀融沉量并研究两者的关系(Chen et al., 2022a)。2018~2021年4个完整观测年的地表变形极差与−40 cm土壤含水量极差具有良好的线性关系(图10),地表季节变形幅度与活动层土壤含水量满足式(3)正比关系,即土壤含水量越高的年份,水和冰转换所产生的体积变化更大,进而导致地表冻胀、融沉幅度更大(陈玉兴等, 2019)。Chen等(2022b)在对青藏公路沿线InSAR地表变形监测中发现,土壤含水量少的区域季节融沉幅度较小,而在土壤含水量明显饱和的区域更大。深度−10~−40 cm(尤其是−40 cm)是该区域土壤水分富集区(图11),是地表变形与土壤温湿度在−40 cm处相关性更显著的主要原因。同时,图11融化期深度−40 cm土壤含水量呈现较大幅度的逐年降低趋势,这也是图9地表冻胀、融沉幅度逐年减小的主要成因。
5. 结论
(1)SBAS-InSAR技术较好地揭露了野牛沟的地表变形特征。结合遥感影像与地表时序变形,发现由冻融过程和水力侵蚀过程引起的地表变形最为显著,地表变形呈现明显的季节性特征,冻融作用引起的年际地表冻胀、融沉幅度为10~20 mm;水力侵蚀引起的年际地表变形通常超过50 mm。
(2)野外观测数据表明大沙龙站具有降温趋势,2017~2021年活动层负温等温线下探深度逐渐增加、历时加长,地表冻结始日和活动层冻结锋面交汇日期逐渐提前。
(3)地表变形与土壤温度、土壤水分具有较好的相关性,在土壤水分富集区(−40 cm处)更显著,相关系数分别为−0.522、−0.415(P<0.001)。土壤水分的变化对地表冻胀、融沉幅度变化的影响更显著,2018~2021年4个完整观测年的地表冻胀、融沉幅度与−40 cm土壤含水量极差具有良好的线性关系。
致谢:大沙龙站观测数据数据来源于“国家青藏高原科学数据中心”(http://data.tpdc.ac.cn),在此一并致谢。
-
-
董英,刘洁,王化齐. 西北地区城市地质调查与城市规划建设[J]. 西北地质,2022,55(03):200–209. DONG Ying, LIU Jie, WANG Huaqi. Urban Geological Survey and Urban Planning and Construction in Northwest China[J]. Northwestern Geology, 2022, 55(03): 200-209.
董英,张茂省,李宁,等. 城市地下空间开发利用的地质安全评价内容与方法[J]. 水文地质工程地质,2020,47(05):161–168. DONG Ying, ZHANG Maosheng, LI Ning, et al. . Methods and contents of geological safety evaluation for urbanunderground space development and utilization[J]. Hydrogeology & Engineering Geology, 2020, 47(05): 161-168.
黄强兵,彭建兵,范文,等. 西安地铁二号线沿线地裂缝未来位错量估算及工程分级[J]. 工程地质学报,2007,15(04):469–474. HUANG Qiangbing, PENG Jianbing, FANwen, et al. Estimation of the maximum displacement of ground fissuresalong Xi'an metro line 2anditsengineeringclassification. Journal of Engineering Geology, 2007(04): 469-474.
马琰,连皓,雷振东,等. 西咸城乡融合发展试验区融合发展路径与策略[J]. 规划师,2021,37(09):61–67. Ma Yan, Lian Hao, Lei Zhendong, et al. . Development Path and Strategy of Xi’an-Xianyang Urban-rural Integration Pilot Area[J]. Planners, 2021, 37(09): 61-67.
彭建兵, 苏生瑞, 张骏. 渭河盆地活动断裂与地质灾害[M]. 西安: 西北大学出版社, 1992 彭建兵,范文,李喜安,等. 汾渭盆地地裂缝成因研究中的若干关键问题[J]. 工程地质学报,2007,(04):433–440. PENG Jianbing, FAN Wen, LI Xi'an, et al. Some Key Questionsin the Formation of Ground Fissuresin the Fen-Wei Basin[J]. Journal of Engineering Geology, 2007(04): 433-440.
王亚辉,张茂省,师云超,等. 基于综合物探的城市地下空间探测与建模[J]. 西北地质,2019,52(02):83–94. WANG Yahui, ZHANG Maosheng, SHI Yunchao, et al. . Precise Detection and Modeling of Urban Underground Space Based onIntegrated Geophysical Exploration[J]. Northwestern Geology, 2019, 52(02): 83-94.
张家明. 西安地裂缝研究[M]. 西安: 西北大学出版社, 1990, 68−88 张茂省, 董英, 刘江, 等. 关中平原城市地质[M]. 北京: 科学出版社, 2021 张茂省,董英,刘洁. 论新型城镇化中的城市地质工作[J]. 兰州大学学报(自然科学版),2014,50(05):581–587. ZHANG Maosheng, DONG Ying, LIU Jie. Discussion of urban geological work in new urbanlization[J]. Journal of Lanzhou University (Natural Sciences), 2014, 50(05): 581-587.
张茂省,董英,张新社,等. 地面沉降预测及其风险防控对策—以大西安西咸新区为例[J]. 中国地质灾害与防治学报,2013,24(04):115–118+126. ZHANG Maosheng, DONG Ying, ZHANG Xinshe, et al. Prediction of land subsidence and Its mitigation methods-a case studyin the new urban district of Xi'an-Xianyang. The Chinese Journal of Geological Hazard and Control, 2013, 24(04): 115-118+126.
张茂省,王化齐,王尧,等. 中国城市地质调查进展与展望[J]. 西北地质,2018,51(04):1–9. ZHANG Maosheng, WANG Huaqi, DONG Ying, et al. . Progress and Prospect of Urban Geological Survey in China[J]. Northwestern Geology, 2018, 51(04): 1-9.