青海玛多MS7.4地震对周边活动断裂的库仑应力加载及发震概率增量的计算

2021年5月22日2时4分位于青藏高原巴颜喀拉地块内部青海果洛州玛多县发生了M_S7.4地震(玛多M_S7.4地震),震中位于98.34°E,34.59°N,震源深度17 km(中国地震台网中心测定).震后区域应力调整直接导致了周边断层的库仑应力变化和主震断层外余震的发生,进而影响到区域地震活动性并改变了即发地震概率.本文利用有限断层滑移模型计算了青海玛多地区主要断层的库仑应力变化,结合Dieterich地震发生率模型,给出了库仑应力扰动下超越某一特定震级地震的发震概率公式,得到了受青海玛多M_S7.4地震所产生的库仑应力场的影响周边8条断裂(段)分别发生M_S≥7.0和M_S≥6.0地震的概率变化.结果表明,8条断裂(段)的发震概率均有不同程度的增加.对于发震概率增幅较为明显的甘德南缘断裂、玛多—甘德断裂和西藏大沟—昌马河断裂在震后短时间内(约10年内)发震概率迅速提升,之后趋于平稳,不排除有潜在发生破坏性地震的可能,特别是发生M_S≥6.0地震的可能性.针对达日断裂,由于库仑应力增加造成的影响不明显,短期内发生M_S≥7.0或M_S≥6.0地震的可能性不大,但随着时间的推移,不排除几十年后有潜在发生破坏性地震的可能,尤其要关注发生M_S≥6.0地震的可能性;东昆仑断裂特别是玛沁—玛曲段仍然是未来强震发生的可能区段,需重点关注和防范M_S≥6.0地震乃至M_S≥7.0地震的发生;玉树—甘孜断裂发生M_S≥6.0地震,尤其是发生M_S≥7.0地震的危险性不高,而乌兰乌拉湖—玉树南断裂有潜在发生破坏性地震的风险,对于M_S≥6.0以上地震需要加强防范.     

Observation and modeling of NPP for Pinus elliottii plantation in subtropical China

Based on the stem analysis of 59 individuals of Pinus elliottii in combination with tree biomass models, we calculated annual biomass increment of forest plots at Qianyanzhou Ecological Station, Chinese Academy of Sciences in subtropical China. In addition, canopy layer and community NPP were calculated based on 12 years’ litter fall data. NPP of the 21-year-old forest was estimated by using the BIOME BGC model; and both measured NPP and estimated NPP were compared with flux data. Community biomass was 10574 g · m⁻²; its distribution patterns in tree layer, shrub layer, herbaceous layer, tree root, herbaceous and shrub roots and fine roots were 7542, 480, 239, 1810, 230, 274 and 239 g · m⁻², respectively. From 1999 to 2004, the average annual growth rate and litter fall were 741 g · m⁻² · a⁻¹ (381.31 gC · m⁻² · a⁻¹) and 849 g · m⁻² · a⁻¹ (463 gC · m⁻² · a⁻¹), respectively. There was a significant correlation between annual litter fall and annual biomass increment; and the litter fall was 1.19 times the biomass increment of living trees. From 1985 to 2005, average NPP and GPP values based on BGC modeling were 630.88 (343.31–906.42 gC · m⁻² · a⁻¹) and 1 800 gC · m⁻² · a⁻¹ (1351.62–2318.26 gC · m⁻² · a⁻¹). Regression analysis showed a linear relationship (R ²=0.48) between the measured and simulated tree layer NPP values. NPP accounted for 30.2% (25.6%–32.9%) of GPP, while NEP accounted for 57.5% (48.1%–66.5%) of tree-layer NPP and 41.74% (37%–52%) of stand NPP. Soil respiration accounted for 77.0% of measured tree NPP and 55.9% of the measured stand NPP. NEE based on eddy covariance method was 12.97% higher than the observed NEP. Supported by the National Key Basic Research Special Foundation of China (Grant No. 2002CB4125), International Joint Research Project under Ministry of Science and Technology of China (Grant No. 2006DFB91920)     

一种以低秩矩阵重建的图像混合噪声去除算法

近年来,低秩矩阵重建在机器学习、图像处理、计算机视觉与生物信息学等众多科学与工程应用领域中,迅速发展为一个新的研究热点,其主要涉及矩阵填充与稳健主成分分析2大问题,即分别从精确且不完全的采样矩阵元与从大误差矩阵元的分布较为稀疏的观测矩阵中恢复出原始低秩矩阵。鉴此,本文定义了稳健矩阵填充,即从非完全且存在稀疏误差的采样矩阵元中精确恢复出原始低秩矩阵,通过最小化核范数与l₁-范数的组合构建了相应的凸优化模型,并提出了一种新颖的增广分部拉格朗日乘数法来求解此类最优化问题。通过将其应用于混合高斯与椒盐噪声去除的问题中表明,此算法对具有规则纹理及相似结构内容等低秩特征的影像中混合噪声的去除效果较好,其能同时去除影像中的椒盐噪声与高斯噪声,且有效保留影像中的纹理细节等信息;当影像中椒盐噪声密度较高而高斯噪声相对较小时,其去噪性能更佳。     

结构性土压缩变形本构描述

在孔隙比与应力半对数坐标中,重塑土等向压缩线通常从正常固结线(NCL)下方逐渐趋近于NCL,整个过程中压缩线切线斜率大小单调增大。而对结构性土进行等向压缩时,因结构性的存在,压缩线会首先越至NCL之上,而后再从上方逐渐趋近于NCL,整个过程中压缩线切线斜率先增大后减小。重塑土与结构性土压缩线切线斜率演化在单调性方面的差异,可通过数学上增量函数构造的方法加以实现。如果在平面直角坐标系中以一水平线作为参考线,描述一过原点的函数曲线,使其上任意一处的切线斜率由该处函数值与参考线函数值的差决定,则参考线是否移动决定了所描述函数切线斜率演化的单调性。受此启发,在孔隙比与应力半对数坐标中提出移动正常固结线(MNCL)的概念。用MNCL代替NCL作为参考线,结合统一硬化模型(UH模型)的理论框架,建立了结构性土压缩变形的本构描述。     

不同方式处理后软土地基侧向变形规律

高速公路建设期侧向位移的控制是评价路基变形的重要依据。依据四川省遂－资（遂宁至资阳）高速公路软基沉降变形观测的数据,通过对现场监测资料的分析,分别探讨西南地区软土路基经塑料排水板（PVD）、碎石桩处理后,路基侧向变形y与深度z,最大侧向变形增量△ym与地表沉降增量△SD,平均侧向变形量Sy与地表沉降量Sf的变化规律。 研究表明,（1）根据Boussinesq解答和实测数据得到路堤荷载作用面积下变形y随着深度z大致呈“S”形变化,近地表变化剧烈,随荷载增加最大侧向位移深度基本保持不变,土体最大侧向变形深度与路堤填筑高度无关,与土性有关;（2）经PVD处理后的软基在填筑期、预压期均有大量侧向位移发生,各占总位移量的50%,最大侧向位移深度在距离地表1.0～1.5 m处;（3）经碎石桩处理后软基侧向位移在填筑期已完成总位移量的75%～80%,总位移量明显小于PVD处理后软基侧向位移量,最大侧向位移深度在距离地表2.5～3.5 m处;（4）PVD处理软基在固结阶段△ym/△SD、Sy/Sf分别为0.15～0.30和0.10～0.20,大于填筑前期2～4倍,而碎石桩处理软基固结阶段△ym/△SD、Sy/Sf分别为0.10～0.15和0.03～0.05,各个阶段变化相当。     

在大地电磁二维Occam反演中求取拉格朗日乘子方法改进

Occam反演具有可靠性和加上各种先验信息可行性的优点,但是因为求取拉格朗日乘子需要很长时间,所以前人提出了几个直接选取方法。笔者改进了二维大地电磁测深Occam反演中求取拉格朗日乘子的方法。在搜索过程中加上观测数据拟合差与其期望值的比较。这样,在Occam反演的拟合差下降阶段的最后迭代和模型光滑阶段,能够排除多余的正演计算。这种改进对反演结果没有影响,但却使得反演速度加快并且更符合于Occam思想。根据模型实验和野外数据处理,可以减少正演次数20%~50%。     

基于年际增量方法的云南省普洱雨季开始期气候预测

基于普洱雨季开始期年际增量变化规律和影响雨季开始期的环流形势及物理过程,采用年际增量方法和多元线性回归分析方法,选取5个具有物理意义的预测因子(包括前期1月南半球绕极环流、前期2月南太平洋高压、前期4月孟加拉湾至南海海平面气压、前期冬季加拿大北部海冰和前期冬季伊朗高原积雪深度),建立了普洱雨季开始期的预测模型,并对预测模型进行1967—2017年的交叉检验和1998—2017年的逐年独立样本检验。交叉检验中,雨季开始期预测值和观测值年际增量的相关系数为0.84,相对均方根误差为24%;独立样本检验中,雨季开始期年际增量的相对均方根误差为15%,模型对雨季开始期异常年份的预测误差小于7 d,表明该预测模型能很好再现1967—2017年雨季开始期的变化趋势。     

加锚节理面的抗剪试验研究

通过一系列剪切模型试验,研究在不同法向力、预应力和锚固情况下锚索对平滑节理面抗剪性能的影响,利用对比试验将销钉力、预应力和轴力增量3种力的抗剪作用分离出来,着重探讨了在剪切过程中3种力对节理面的抗剪贡献,并分析了锚索模拟体的破坏形式。试验结果表明,预应力锚固明显提高了节理面剪切力-位移曲线初始段的抗剪力和剪切刚度,并使得曲线末段表现出塑性强化的特征;预应力大小对加锚节理面抗剪强度的影响与试验条件、节理面状况和锚固角等有关。试验和理论分析表明,剪切位移开始阶段预应力和销钉作用发挥了主要作用,杆体形成塑性铰后轴力增量则起了重要的作用,锚索的破坏一般发生在节理面附近,属于拉剪复合破坏模式。     

太阳系真实力模型下Halo轨道的寻找

在总结现阶段共线平动点附近周期和拟周期轨道寻找方法的基础上,在JPL星历提供的真实力模型下利用目标轨道法(Target Shooting Method)实现了Halo轨道的寻找。即利用Lagrange乘数法作为迭代工具,利用圆型限制性三体问题下的分析解结果作为迭代初始值,迭代出真实力模型下的拟Halo轨道。重点构建了真实力模型下的虚拟平动点准惯性坐标系,极大地简化了轨道坐标从CRTBP模型向真实力模型的转换,加快了迭代收敛速度。