《Plants》：SCEA-YOLO: A General-Purpose Maturity Grading Model of Multi-Crop Greenhouse Robots
Tianyuan Li,
Ping Liu,
Dongfang Song,
Xingtian Zhao,
Xiangyu Lyu and
Kun Zhang
编辑推荐：
本刊推荐一项关于大型调水工程（SNWDP）调控下地面变形动态响应的创新研究。为揭示水资源政策转变如何影响地面变形，研究者综合运用InSAR（干涉合成孔径雷达）时序分析与ICA（独立成分分析），系统剖析了北京平原2015-2023年的沉降-回弹时空模式。研究发现，南水北调工程实施后，区域变形呈现准线性沉降、潮白河回弹及“沉降-回弹”转换三主导模式，并首次量化了调水启动与区域回弹响应间存在的5.5年时滞。该成果为类似缺水型大都市的地下水管理与地面沉降控制提供了科学范式。
想象一下，一座拥有两千多万人口的超级大都市，正因长期“口渴”而缓慢下沉。北京，作为中国政治文化中心，长期以来面临严重的水资源短缺。为满足城市发展需求，地下水被大量开采，导致地下水位持续下降，并引发了显著的地面沉降。这种沉降不仅威胁城市建筑与基础设施安全，还可能加剧沿海地区的海平面上升风险。然而，自2014年底南水北调中线工程（South-to-North Water Diversion Project-Middle Route, SNWDP-MR）通水以来，大量长江水北送，深刻改变了北京的用水结构，部分替代了地下水开采。一个令人关注的新现象随之出现：一些地区的地面沉降得到了遏制，甚至出现了回弹的迹象。在这种新的水文背景下，北京平原的地面变形机制发生了怎样的转变？从沉降到回弹的动态过程是如何演化的？其背后的驱动因素又是什么？回答这些问题，对于科学评估调水工程效益、优化区域水资源管理和防控地质灾害具有至关重要的意义。近期，一项题为“SCEA-YOLO: A General-Purpose Maturity Grading Model of Multi-Crop Greenhouse Robots”的研究（注：此处用户提供的文档标题与后续实际内容明显不符，实际内容为关于北京地面变形的研究）在《Plants》上（注：此期刊信息与文档内容不符，文档内容更可能出自遥感或地球科学类期刊）发表，系统揭示了这一转变过程的时空特征与内在机制。
为了深入探究北京平原地面变形（沉降-回弹）的时空特征、主导模式及驱动机制，研究人员开展了一项综合性研究。他们首先利用C波段的哨兵1号A星（Sentinel-1A）2015年至2023年的合成孔径雷达（Synthetic Aperture Radar, SAR）数据，采用改进的小基线集干涉合成孔径雷达（SBAS-InSAR）技术，获取了北京平原长时间序列、高精度的地面形变信息。在此基础上，研究运用K-means聚类方法对变形时空特征进行分析。更为关键的是，为了从复杂的混合信号中剥离出独立的驱动源，研究引入了独立成分分析（Independent Component Analysis, ICA）这一信号分解方法，成功提取了不同物理机制主导的变形时空成分。整个研究框架整合了InSAR时序分析、ICA以及区域水文地质数据，旨在系统解析由用水变化驱动的地面变形转变机制。
主要关键技术方法
研究采用的核心技术方法包括：1. 改进的SBAS-InSAR处理：基于172景Sentinel-1A降轨SAR数据（2015-2023），利用ISCE和MintPy软件进行干涉处理与时序反演，获取垂直向形变速率和时间序列，并通过GPS（BJFS站）数据验证了结果的可靠性。2. K-means聚类分析：对InSAR形变时间序列进行无监督聚类，依据“肘部法则”确定最优聚类数（K=4），从现象学角度识别了变形严重程度不同的空间区域。3. ICA信号分解：首先对InSAR时序数据进行主成分分析（PCA）降维（保留前4个主成分，涵盖96.62%方差），随后利用FastICA算法进行盲源分离，提取了统计独立的时空成分，从而分离出不同物理过程驱动的变形信号。
研究结果
3.1. InSAR结果分析
2015-2023年间，北京平原的垂直变形速率分布呈现高度空间异质性。沉降区主要集中在南部、西北部和中部偏东区域，形成了明显的条带状格局和多个沉降中心。最大沉降速率（-108 mm/yr）出现在朝阳-通州走廊地带。选取潮白河上游（C1）、朝阳-通州沉降带内（Z1）和永定河沿岸（Y1）三个代表点的位移时间序列分析显示，Z1点沉降显著但2019年后速率明显减缓，Y1点无明显沉降，而C1点则呈现短期回弹趋势。与BJFS站GPS数据的对比验证表明，InSAR反演的垂直位移具有高可靠性。
3.2. K-means结果分析
K-means聚类将研究区划分为4类。第1类（Cluster 1）代表沉降最严重的区域（累计位移26-50 cm），空间上集中于朝阳-通州走廊；第2、3类为过渡区域，沉降程度依次减轻；第4类为相对稳定区。该分析从现象上揭示了北京平原地表变形的空间分异格局，识别了严重沉降带、过渡区和稳定区。
3.3. ICA结果分析
ICA成功分解出四个独立成分（IC）。IC1和IC3被识别为具有明确物理意义的变形信号，而IC2和IC4表现为随机波动，被解释为噪声。IC1的空间贡献图显示，沉降（蓝色）主要集中于朝阳区东部和通州区北部，其时间序列呈近乎线性的单调下降趋势，代表了“准线性沉降”模式。IC3的空间贡献图以蓝色（负贡献）为主，结合其时间序列（先增后减），揭示了一种广泛的“沉降-回弹”转换模式。
4. 讨论
4.1. 准线性沉降模式的驱动因素
IC1所代表的准线性沉降模式，其驱动机制是深层地下水（第三、四含水组）的长期超量开采。该区域深层含水组（由中粗砂和砾石构成）上覆巨厚（>150米）的可压缩粘土地层。由于埋深大（180-300米），自然补给能力极弱。长期超采导致孔隙水压力下降，有效应力增加，引发粘土地层的不可逆（塑性）压实，从而表现为持续的地面沉降。该模式揭示了与短期气候波动或地表水管理政策相对“脱钩”的、根植于深层地下水开采的沉降驱动力。
4.2. 潮白河地区回弹模式的驱动因素
自2015年以来，在潮白河上游砂卵石含水层系统实施了多次人工地下水回灌。ICA结果显示该区域叠加了IC1（回弹）和IC3（沉降-回弹转换）模式。C1点显示最大垂直回弹速率可达20 mm/yr，证实了人工回灌措施的有效性。该区域回弹与沉降速率的不对称性（最大沉降108 mm/yr vs. 峰值回弹20 mm/yr）归因于上游潮白河地区以砂卵石地层为主，其压缩变形以弹性为主，不同于富含厚层塑性粘土的区域。
4.3. 沉降-回弹转换模式的驱动因素
IC3分量捕捉了北京平原地表对南水北调工程最广泛和显著的水文响应，阐明了其驱动机制、量化了时滞并强调了管理意义。
4.3.1. 浅层地下水调控与降水下渗的耦合
“沉降-回弹”转换主要受浅层地下水恢复和增强的降水下渗共同驱动。南水北调中线通水后，地下水开采减少和人工回灌（通过未衬砌河道约4.5×108m3）导致浅层含水层水位持续上升。同时，2014年以来年降水量频繁超过500毫米，为渗透性高的砂砾石地层提供了充足补给。这种双重机制促进了含水层骨架的弹性膨胀，逆转了长期沉降趋势。IC3的空间分布与受益于这种“政策+气候”效应的区域高度相关。
4.3.2. 南水北调工程诱发的滞后回弹
一个关键发现是量化了水文干预与区域地质力学响应之间的显著时滞。对IC3时间序列的分析揭示，南水北调供水开始（2014年12月）与区域回弹广泛开始（2019年5月）之间存在约5.5年的滞后期。这一时滞体现了两个关键水文地质过程：1) 调水入渗通过粘性弱透水层并补给浅层含水层所需的时间；2) 含水层骨架在有效应力减小下从压缩状态逐渐调整至弹性膨胀的过程。分段线性回归显示，IC3的变形速率在2019年前后存在不对称性（沉降趋势斜率绝对值|0.0179| > 回弹趋势斜率绝对值|0.0097|），这直接体现了可压缩粘土层的非弹性（塑性）变形，符合太沙基（Terzaghi）有效应力原理。塑性变形导致颗粒排列稳定，即使地下水位回升也难以完全恢复，造成回弹速率慢于沉降速率。
4.4. 基于变形模式的管理建议
ICA解码的变形模式时空异质性，为制定有针对性的地下水管理策略提供了科学依据。对于受IC1信号主导的区域（如朝阳-通州走廊），应严格控制第三、四含水组的开采，优先使用替代水源（如南水北调水），以降低地下岩层的有效应力。对于受IC3信号主导的区域（如潮白河、永定河流域），应加强工程化回灌实践，并利用未衬砌河道下渗，放大“降水+人工回灌”的耦合效应，加速浅层含水层回弹。需要建立分区管理框架，实现水资源安全与地质灾害防控的科学平衡。
研究结论与意义
本研究系统揭示了新水文背景下北京平原地表变形（沉降-回弹）的时空特征、主导模式及驱动机制。主要结论包括：识别了准线性沉降、潮白河回弹及广泛“沉降-回弹”转换三主导模式；阐明其驱动机制分别为深层地下水超采引发的塑性压实、人工回灌诱发的弹性回弹、以及政策调水与降水增加共同驱动的浅层含水层弹性膨胀；首次量化了南水北调启动与区域回弹间约5.5年的响应时滞，并揭示了变形不对称性（沉降斜率大于回弹斜率）与含水层岩性（粘土 vs. 砂砾石）的内在联系。
这项研究的意义重大。它从机理上解析了大型调水工程如何通过改变区域水循环，进而影响地质环境过程，为评估此类巨型工程的综合环境效应提供了新颖的视角和定量的方法。研究所采用的InSAR-ICA集成分析框架，能够有效分离复杂变形信号中的独立驱动源，不仅适用于北京，也可为全球面临水资源短缺、地下水超采和地面沉降威胁的类似大城市提供可转移的科学参考。最终，研究强调必须转向分区的地下水管理策略：在东部沉降走廊严格限制深层开采，在西北冲积扇区主动增强人工回灌。这标志着在地面变形动力学领域，从被动监测到基于机理理解进行主动调控的范式转变。
