基于蒙特卡洛生存分析探究东北森林物候的影响因素

植被是生态环境变化的指示器,分析植被物候的影响因素不仅有助于气候变化分析,提高区域气候模式的模拟精度,而且对于准确评估植被生长趋势、生产力以及全球碳收支均具有重要意义。基于遥感的植物物候监测已取得了长足的发展和进步,但当前利用大范围、长时间序列的遥感数据分析植被物候影响因素的研究尚不多,采用线性回归模型对非线性的植被物候影响因素进行分析可能存在偏误。因此,本文提出一种基于蒙特卡洛模拟的生存分析方法,对东北森林物候的影响因素进行量化分析。首先利用东北森林地区1982-2009年间AVHRR GIMMS NDVI数据,应用双Logistic曲线拟合方法对植被春季返青期（SOS）、秋季落叶期（EOS）及植被生长期（GSL）进行提取;然后基于蒙特卡洛模拟和生存分析构建植被物候影响因素分析模型;最后运用所构建模型探讨了东北森林区春季返青期、秋季落叶期的可能影响因素。结果发现：温度、降水和风力对中国东北森林关键物候期有一定影响,其中温度是春季返青期和秋季落叶期的最主要驱动因素,长期平均温度比短期内的温度突变对物候影响更显著,落叶期前的风速增加有可能使落叶时间提前;除了环境因素,春季返青早的年间秋季落叶倾向于更晚。研究表明,结合蒙特卡洛方法的生存分析可以较好地对物候期的影响因素进行定量分析,可为物候现象的归因分析提供一种新的方法。     

大兴安岭诺敏河火山群远震P波衰减研究

板内火山活动是板内下方深部物质组成的重要视窗,而中国东北是研究板内火山的热点场所.中国东北火山活动主要沿大兴安岭及松辽盆地周缘分布,隐匿在大兴安岭内的各火山群落鲜有地球物理观测.本研究聚焦于大兴安岭北部的诺敏河火山群,利用围绕诺敏河火山布设的43个流动台站记录到的17个远震事件数据,通过时间域波形匹配法计算了300条直达P波震相的Δt~*,并进一步采用贝叶斯蒙特卡洛方法反演得到了研究区的二维Δt~*模型,为约束诺敏河火山地区的深部结构提供了新的观测证据.结果显示研究区高衰减区域与火山活动在空间位置上有很好的一致性.诺敏河火山与研究区北端临近小古里河火山区域均观测到高Δt~*值,可能由区域下方热地幔物质上涌导致.诺敏河火山东部、科洛河火山西部区域观测到低Δt~*值,与高衰减区域最大Δt~*差值达约0.3±0.05 s,可能由两部分原因导致:该区域下方残留未被侵蚀的岩石圈,以及局部区域或存在的部分熔融导致的残留地幔矿物脱水.     

FSRQ 0208-512的光学波段长周期光变分析和色指数变化研究

从SMARTS (The Small and Moderate Aperture Research Telescope System)数据库,收集到了平谱射电类星体(FRSQ) 0208-512天体光学波段B、 V、 R的观测数据,数据从2008年2月–2017年8月,时间跨度近10 yr.分别获得了光学B波段数据627组、V波段数据634组、 R波段数据639组.通过加权小波Z变换法(WWZ)和蒙特卡洛模拟的LombScargle periodogram (LSP)方法研究得出了如下结论:(1)在两种方法中都得到了光学B、V、R 3个波段有可能存在396 d左右约为1.08 yr的周期性光变;(2)得到了B波段星等B与光学波段色指数(B-V)呈正相关关系,当B星等越小时色指数(B-V)越小,即越亮越蓝;(3)通过光学3个波段的1.08 yr的长周期光变得到天体中心黑洞质量M约为0.12×106M⊙,辐射区半径RM=1.052×108km,其中M⊙为太阳质量.     

海原断裂带老虎山段晚第四纪滑动速率精确厘定与讨论

如何准确测定断裂滑动速率是近年来活动构造研究的前沿与热点.随着高精度地形数据获取手段与第四纪测年方法的不断进步,位错量和地貌面年龄的精度均得到大大提高.在进行滑动速率计算时还要考虑地质过程是否合理,蒙特卡洛方法为获取更加符合地质过程的滑动速率提供了重要工具.本文以滑动速率研究程度较低的海原断裂带老虎山段为例,基于LiDAR高精度地形数据,测得T1—T4阶地面年龄分别为1~3 ka,9~11 ka,15~17 ka,40~45 ka,陡坎前缘的位错分别为7~14 m,28~36 m,59~66 m,180~190 m.综合多地点的左旋走滑位错量及不同时代的地貌面年龄数据,并考虑滑动历史,利用蒙特卡洛模拟方法,将位错-时间两个参数的不确定性定量化,限定老虎山断裂45 ka以来平均滑动速率为4.3±0.16 mm·a^-1,17ka以来的平均滑动速率为4.0±0.15 mm·a^-1,与前人研究得到的狭义海原断裂滑动速率4.5±1.0 mm·a^-1基本一致.综合整个海原断裂带滑动速率,本文结果更支持低滑动速率变化趋势,即海原断裂带整体滑动速率趋于稳定,向东至六盘山断裂,滑动速率开始降低,推测海原断裂带的左旋走滑在尾端主要为马东山—六盘山隆起所吸收.结合老虎山断裂历史地震资料和深部锁闭浅部蠕滑的动力学特征,推测老虎山断裂具备与相邻断裂一起触发强震的能力.     

地震巨灾模型中的随机地震事件集模拟

地震巨灾保险是降低地震灾害风险的有效手段之一,而地震危险性分析是地震巨灾模型的主要分析模块之一。传统的概率地震危险性分析主要是基于潜在震源模型、地震活动性模型和地震动衰减模型等并采用概率方法得到场点的地震危险性值,该危险性表示的是未来所有地震对场点的综合影响。然而,在使用地震巨灾模型进行地震风险分析时需要用到单个地震事件对场点的影响,这就需要根据潜在震源区生成一系列单个地震事件,并计算每个事件对场点的影响。本研究采用蒙特卡洛方法,基于第五代中国地震动参数区划图中所采用的地震活动性模型(潜在震源区及其地震活动性参数),模拟符合我国地震活动时间、空间和强度分布特征的地震事件集。模拟时遵从的基本理论为:地震发生在时间上符合泊松分布,震级分布可用古登堡-里克特定律来描述,空间分布特征则用潜在震源区及其地震发生率来描述。模拟得到的地震事件包含以下参数:时间(年、月、日)、地点(经度、纬度)、深度、震级、断层走向以及衰减特征等。该模拟地震事件集可满足地震巨灾模型中地震风险分析的需求,已应用于我国地震巨灾模型中。     

同震崩塌滑坡的光学遥感影像多特征融合解译方法

同震崩塌滑坡的解译及定位是震区灾后恢复工作中需要重点解决的问题,如何在灾害快速、自动解译的基础上,不断提高解译精度,是目前同震崩塌滑坡解译的研究热点之一,也是促使地质灾害早期识别向智能化、科学化发展的必要前提。文章在团队前期所提出的遥感影像局部阈值二值化方法的基础上,针对同震崩塌滑坡解译结果假阳率偏高的问题,分析了假阳性地物的光学和几何特点,提出了融合目标区域光学影像灰度特征、区域坡度信息、NDVI指数特征及解译地物主轴特征的同震崩塌滑坡多特征融合解译方法。为验证所提出模型的准确性,以2014年云南鲁甸地震龙头山镇为研究区,利用震后获取的高分一号（GF-1）卫星影像数据及数字高程模型对该同震崩塌滑坡进行了解译识别,结果表明,文中提出的方法准确解译出了同震崩塌滑坡区域,并有效去除了假阳性地物干扰,提高了解译精度。     

广西合浦盆地石油勘探决策分析

将现代决策理论应用于合浦盆地石油勘探决策分析,提出了一个我国南方小盆地勘探中的决策分析模式。首先,分析盆地的资源量,在资源量可观的情况下,先对圈闭进行地质评价,然后,利用蒙特卡洛法模拟计算优选圈闭的勘探风险,通过给出各参数的概率分布,计算出各圈闭在不同置信水平下的储量和净现值,最后根据实际情况做出决策。     

提问蔚蓝

<正>~~     

基于蒙特卡洛模拟的概率断层位错危险性分析

同震位错对川藏铁路等跨断层工程的安全造成严重威胁,合理评价活断层的位错参数具有重要的应用价值.由于传统的确定性评价方法存在无法区分工程场点的重要程度和其在断层上的相对位置等缺陷,越来越多的学者推荐采用概率断层位错危险性分析（PFDHA）.然而基于经典的概率性方法开展PFDHA原理复杂且实现困难,不利于吸收断层地震活动性研究的最新成果,也有碍于PFDHA的普及和推广.相比经典的概率性方法,蒙特卡洛模拟具备逻辑清晰易懂、程序易于实现且兼容性和扩展性好的优点.本研究基于蒙特卡洛模拟实现了概率断层位错危险性分析的一般性算法,并将该方法应用于鲜水河断裂带的炉霍段.结果显示,PFDHA的结果随着超越概率水准或工程场点在断层上的相对位置的不同而显著变化.适当考虑最大同震位错和地表破裂长度的不确定性得到的位错参数更加合理.超越概率大于等于100年2%时,PFDHA的结果显著小于确定性方法的结果.然而随着断层活动性的提高,100年超越概率1%的结果可能会大于确定性方法的结果.按照不同类型工程的抗震设防水准选择相应的PFDHA评价结果,既有利于工程的安全,也有助于大多数工程节约成本.PFDHA相比确定性方法具备多种优势,有望为川藏铁路等重大工程的抗断参数评估提供技术支撑.      

宇生核素深度剖面暴露测年的原理与应用

宇生核素暴露测年技术是目前第四纪地貌年代测定最常用的方法之一,由于难以定量化继承性核素以及侵蚀对样品暴露年代的影响(通常假设继承性和侵蚀的影响为0),其测年结果难以反映地貌的真实暴露年代。通过测定地貌剖面不同深度样品的宇生核素浓度,建立“宇生核素浓度—深度”剖面曲线,可同时定量化研究地貌体的暴露年龄、继承性宇生核素浓度以及侵蚀速率,有效地弥补了宇生核素暴露测年中继承性核素以及侵蚀速率不确定性的缺陷。基于蒙特卡洛方法的宇生核素深度剖面暴露测年模型(简称蒙特卡洛深度剖面模型)是最常用的计算模式之一,然而,国内关于该方法的原理及其应用研究相对较少,大大限制了该方法的广泛使用。笔者等详细阐述了蒙特卡洛深度剖面模型的原理及其在冰川地貌、河流阶地、冲洪积扇等地貌测年中的应用。此外,介绍了深度剖面模型的复现、非稳态模型、线性反演模型等其他几种计算模式并进一步探讨了蒙特卡洛深度剖面模型暴露测年技术的野外采样方法,以使其能够广泛应用于第四纪地貌学和年代学研究中。