基于最大特征值估计的C3算法及应用

第三代相干体算法(C3算法)具有分辨率高、压制噪声能力强的优点,但是该算法需要计算协方差矩阵的特征值,所以耗时较多。为提高大矩阵特征值分解的运算速度,提出了一种估计最大特征值的快速收敛算法,通过运用A2k矩阵的迹来构造一种收敛速度更快的特征值估计方法,并给出了这种方法的误差估计表达式,实现计算精度的有效控制。文中还讨论了地震数据体边缘补地震道方法,解决了将递推算法应用到平面方向时空间窗越过边界的问题。实例表明,本算法的计算效率明显高于常规的C3算法。     

中国县域生态资产损益的影响机制及优化提升路径

生态资源资产化、生态资产资本化是实现“两山”转化的一个重要途径,开展生态资产及其损益核算有助于生态系统的有效保护、科学管理和可持续利用。本文在借鉴国内外生态资产相关研究的基础上,构建了县域尺度生态资产及其损益核算技术体系,分析了1990—2018年中国县域生态资产的时空变化特征,明晰了造成生态资产损益的影响因素及作用机制,并针对具有不同生态资源类型、损益状况、驱动要素特征的县域提出了差异化的优化提升路径。结果表明：① 以农田、森林生态资源为主的县域数分别约占县域总数的45%和37%;1990—2018年中国县域生态存量资产的质量整体呈波动中上升趋势,而生态流量资产的变化趋势差异显著,近70%县域水源涵养量减少。② 生态流量资产价值呈现损与益态势的县域数分别约占44%和37%,空间分布呈现与“胡焕庸线”一致的分割特征,即西北广大地区县域明显增益,而东部和南部县域有所减损。③ 1990—2018年中国超过70%县域生态资产变化受气候变化与人类活动共同驱动,人类活动对生态资产增益县域的平均作用率约80%,而气候变化对生态资产减损县域的平均作用率约为60%。④ 综合分析,可将全国县域分为适应、减缓气候变化和生态资源修复、保护、管理5种类型,采用差异化的生态资产优化提升路径。     

基于主体功能区的京津冀城市群碳收支时空分异与碳补偿分区

从主体功能区划视角研究碳收支和碳补偿分区,对于制定适应各主体功能区低碳发展策略,推动区域生态环境协同治理,实现高质量发展具有重要意义,也是地理学思想对于实现“双碳”目标的重要贡献。本文首先构建主体功能区视角下碳收支与碳补偿的理论框架,然后以京津冀城市群157个县级单元为研究区,引入集中化指数、标准显性比较优势指数、SOM-K-means聚类等方法研究功能区视角下京津冀城市群碳收支时空分异与碳补偿分区,并提出以低碳发展为导向的碳减排空间优化方案。结果表明：① 2000—2017年京津冀城市群碳收支量呈现波动上升态势,其集中化指数均高于0.4的“警戒线”,碳收支地区差异整体偏大。② 碳收支时空分异显著,碳排放高值区呈现以京津唐为中心,向外逐步降低的“核心—外围”空间格局;而碳吸收空间格局趋于稳定,总体呈现东、北、西部高,而中、南部低的倒“U”型格局。③ 京津冀城市群碳收支与主体功能区战略定位较为吻合,优化开发区和重点开发区是碳排放的主要承压区,而重点生态功能区是碳吸收的优势主导区,各功能区碳吸收集中化指数差别较碳排放集中化指数小。④ 京津冀城市群共有53个支付区、64个平衡区和40个获补区,结合主体功能区规划战略目标,最终形成9类碳补偿空间优化区,并提出每一类型区低碳发展方向及策略。⑤ 未来要加强更微观尺度的碳收支及碳补偿研究,丰富和完善碳补偿理论框架,将碳补偿融入到碳交易市场,探索实现“双碳”目标的多元化路径。     

基于贝叶斯网络的生态系统服务空间格局优化——以泾河流域为例

城市高速发展引起了区域土地利用格局的改变,不仅影响生态环境质量,而且还对生态系统服务的空间格局产生影响,在此基础上优化生态系统服务显得至关重要。在泾河流域2000—2020年净初级生产力（Net primary productivity, NPP）、农业生产力、土壤保持和产水服务空间评估的基础上,将贝叶斯网络和生态系统服务相结合,在关键变量子集和可视化的最优状态子集的基础上,评估了4种生态系统服务需要优化的区域,为区域经济和生态和谐发展提供参考。结果表明：（1） 水文评价模型（Soil and water assessment tool, SWAT）模型能较准确地模拟区域的径流量。通过模拟值和观测值对比分析,该模型具有较高的决定性系数（R²>0.6）和纳什效率系数（NSE>0.5）,可为进一步评估产水服务提供保障。（2） 2000—2020年泾河流域4种生态系统服务的时空差异性较为显著。在时间尺度上,4种生态系统服务均呈现波动中上升的趋势,在空间尺度上呈现较为稳定的变化趋势。（3） 通过对4种生态系统服务优化区域进行叠加分析,发现综合优化区域集中在彭阳县的中部和西南部以及环县的零星区域。研究结果对指导优化区域进行生态系统的可持续管理以及改善生态系统的退化状况具有重要意义。     

国内土地评价与土地优化利用研究综述

土地利用优化是当前土地资源管理和调控过程中必须解决的难题,本文在对国内土地评价和土地利用优化文献回顾的基础上,总结了土地评价研究内容和技术方法以及土地利用优化研究内容和技术方法的发展。     

Multi-objective autocalibration of SWAT model for improved low flow performance for a small snowfed catchment

ABSTRACT

A reliable modelling framework needs to ensure that the model is simulating reality with limited uncertainty, thus enhancing its predictive ability. In the literature, hydrological model assessment using one or more metrics is reported to be inadequate when the river flow regime is required to be reproduced comprehensively. This research is aimed to: (a) calibrate the Soil and Water Assessment Tool (SWAT) based on the concept of multi-objective optimization by applying the Borg multi-objective evolutionary algorithm (MOEA); (b) apply hydrological signatures as objective functions; and (c) adopt a multi-metric approach for model evaluation. The SWAT model was coupled with a relatively newer and powerful Borg MOEA. The inclusion of hydrological signatures as objective functions along with the conventional statistical functions assisted in improving the performance for low flows by 135% in terms of volume efficiency and 65% for flow time series simulation.     

How the Method of Minimization of Action Avoids Singularities

The method of minimization of action is a powerful technique of proving the existence of particular and interesting solutions of the n-body problem, but it suffers from the possible interference of singularities. The minimization of action is an optimization and, after a short presentation of a few optimization theories, our analysis of interference of singularities will show that:(A) An n-body solution minimizing the action between given boundary conditions has no discontinuity: all n-bodies have a continuous and bounded motion and thus all eventual singularities are collisions;(B) A beautiful extension of Lambert's theorem shows that, for these minimizing solutions, no double collision can occur at an intermediate time;(C) The proof can be extended to triple and to multiple collisions. Thus, the method of minimization of action leads to pure n-body motions without singularity at any intermediate time, even if one or several collisions are imposed at initial and/or final times.This method is suitable for non-infinitesimal masses only. Fortunately, a similar method, with the same general property with respect to the singularities, can be extended to n-body problems including infinitesimal masses.     

基于流溪河模型的新丰江水库入库洪水预报模型研究

基于流溪河模型构建了不同空间分辨率的新丰江水库入库洪水预报模型,采用粒子群算法优选模型参数进行洪水模拟,对比实测资料验证了模型性能。结果表明,建立的新丰江水库入库洪水预报流溪河模型模拟效果优良,符合新丰江水库入库洪水预报的精度要求,可用于新丰江水库入库洪水预报：90 m空间分辨率的模型模拟效果最好,200 m空间分辨率的模型也具有良好的模拟性能,500 m空间分辨率的模型模拟性能有明显降低。模型计算工作量随分辨率的提高呈指数增加,90 m空间分辨率的参数优选时间是200 m空间分辨率参数优选时间的7.04倍,而200 m空间分辨率的参数优选时间是500 m空间分辨率参数优选时间的3.53倍。     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.