结合序贯指示模拟的贝叶斯高斯混合线性反演

高斯混合模型(Gaussian Mixture Model, GMM)可以用来描述储层性质的多峰分布特性,多峰特性主要是由于它们在不同离散变量内的变化而引起的.在高斯混合模型中,高斯分量的权值代表离散变量的概率.然而,基于高斯混合模型的贝叶斯线性反演可能会对某些点的离散变量错误地分类,进而影响连续变量的反演结果,尤其存在强噪声的时候.在本文中,我们考虑了离散变量的空间变化性,并将高斯混合模型与序贯指示模拟(Sequential Indicator Simulation, SIS)相结合来确定离散变量的后验条件权值,形成了结合序贯指示模拟的贝叶斯高斯混合线性反演方法.该方法能够准确地对离散变量进行归类,且具有良好的抗噪性.通过模型试算,我们证明了这种方法的可行性,并在实际资料中取得了较好的结果.     

青藏高原东缘及邻区强震构造:专辑序言

青藏高原东缘由多个次级构造单元组成了独特的“多层次挤出-旋转活动构造体系”,是调节高原物质向东挤出的构造转换带,发育了高密度的活动断裂,强震频度高且强度大,而且地质地貌复杂,成为我国地震灾害问题最为突出的地区之一。近年来,伴随全国地震灾害风险普查、城市活断层探测、地震科学实验场建设、工程场地的地震安全性评价以及重大工程与城镇区的活断层鉴定与地壳稳定性评价等工作的深入开展,高精度遥感、构造地貌、古地震和第四纪年代学等方法的广泛应用,显著提升了青藏高原东缘及邻区强震构造的调查研究程度。为了及时交流这方面的最新研究成果,支撑区域防震减灾及重要工程和基础设施建设的地质安全评价等工作,《地震科学进展》编辑部组织了“青藏高原东缘及邻区强震构造”成果专辑,征集了活断层与地震、甘肃积石山地震以及相关领域的研究综述等代表性学术论文20余篇,本期《青藏高原东缘及邻区强震构造专辑Ⅰ》优选了10篇论文,后续还将推出《青藏高原东缘及邻区强震构造专辑Ⅱ》,希望这些成果可提升对区域强震活动特征与孕震构造机制的理解,并为区域强震危险性分析和有效防范强震灾害风险提供科学依据。     

基于次序统计量理论的洪水频率分析方法研究

针对目前参数估计方法存在稀遇洪水频率设计值偏大的问题,提出采用基于数值次序统计量期望值的适线法计算洪水频率设计初始值;针对目前参数估计方法存在常遇洪水频率设计值偏小的问题,提出采用周文德公式修正目前频率设计偏小的情形,计算洪水频率设计修正值;针对目前传统方法计算洪水频率设计值抽样误差大的问题,提出采用基于数值次序统计量标准差的方法计算洪水频率设计保证值。定量分析表明,提出的基于数值次序统计量期望值的优化适线法具有良好的无偏性和有效性;寸滩站洪水频率设计修正值大于频率设计初始值,且两者偏差随着频率减小而减小;寸滩站洪水频率设计保证值明显大于修正值。     

贵州山区降水集中度和降水集中期的时空变化特征

利用1961~2015年贵州省81个气象观测站的逐日降雨资料,计算了降水集中度(PCD)和集中期(PCP),并采用线性趋势、空间分析、M-K突变检验和R/S分析等方法 ,对贵州山区复杂地形下的降水集中度和降水集中期进行了分析。结果表明:贵州1961~2015年PCD波动幅度介于0.30~0.59之间,空间分布上由东向西逐渐增加的趋势,PCP的波动范围在16.5~21.0旬之间,空间呈现出由东南向西北逐渐增加的趋势;PCD与降水量的相关性高于PCP与降水量的相关性,表明贵州地区的降水量越多,降水越趋于集中;R/S分析表明贵州地区PCD可能隐含着周期规律,PCP会保持提前的趋势。     

哀牢山-红河剪切带渐新世的构造体制转换与剥露历史:来自哀牢山南段磷灰石裂变径迹的证据

哀牢山-红河剪切带是东南亚重要的构造边界,其记录了青藏高原东南缘新生代以来的陆内变形和地貌演化。本次研究对该剪切带哀牢山南段开展了基于LA-ICPMS法测试的磷灰石裂变径迹低温年代学分析。磷灰石裂变径迹年龄数据和热史反演模拟揭示哀牢山段存在晚始新世-早中新世(40～20Ma)的快速剥露事件,而早中新世(大约20Ma)之后处于稳定的慢速剥露过程。磷灰石裂变径迹年龄-海拔分布曲线特征暗示:快速剥露机制存在差异,早期阶段(40～26Ma)的剥露过程受控于伸展为主的左旋走滑体制影响;晚阶段(26～20Ma)的快速剥露归因于简单剪切为主的左旋走滑剪切体制,上述结果暗示哀牢山-红河构造带在晚渐新世发生了一次重要的构造体制转换,即从走滑伸展变形转换为简单剪切变形。哀牢山杂岩带北段、中段、南段冷却路径对比,表明北-中段可能存在两阶段快速冷却作用,而南段只发生单一快速冷却作用;结合青藏高原东南缘低温热年代学数据,暗示自中-晚中新世,青藏高原中、下地壳物质可能向东南缘扩展,并已到达哀牢山中段,同时诱发哀牢山杂岩带以北广大地区的抬升和快速冷却。     

X射线荧光光谱法同曲线测定海洋沉积物和陆地地化样品中的29个主次痕量元素

本法采用粉末压片,使用Axios X射线荧光光谱仪同曲线测定海洋沉积物和陆地地球化学样品中钠、镁、铝、硅、磷、氯、钾、钙、钛、钒、铬、锰、铁、钴、镍、铜、锌、镓、砷、溴、钇、铌、锆、锶、铷、铅、钍、钡、铪共29个主次痕量元素。针对地质样品的多样性,通过多种校准样品拟合同曲线测量,降低了检出限,并有效解决了Cl、As等元素的测量范围局限性。用经验系数、理论alpha系数与经验系数相结合和康普顿散射线内标法校正元素间的基体效应。分析结果的精密度和准确度能满足《地质矿产实验室测试质量管理规范》中的相关要求。     

渤、黄、东海高营养层次重要生物资源种类的营养级研究

利用2000年和2001年2次大面调查所收集的11970个胃含物样品分析结果,计算了黄海和东海生态系统高营养层次35个重要生物资源种类的营养级,同时,结合对渤海和黄海39个种类营养级历史数据的修正,讨论研究了我国海洋高营养层次生物资源种类营养级的研究策略和计算方法。主要研究结果为：(1)渤海重要生物资源种类营养级的变化范围为3．12～4．9,黄海为3．2～4．9,东海为3．29～4．55。近年来各海域高营养层次的营养级呈下降趋势,如渤海从1959年的4．1下降到1998～1999年的3．4,黄海从1985～1986年的3．7下降到2000～2001年的3．4;(2)高营养层次营养级波动主要是由于群落种类组成变化及单种类营养级年间波动引起的,而单种类营养级年间波动又直接与群体个体变小以及摄食食物的低营养层次化有关。因此,高营养层次的营养级变化是认识海洋生态系统生物生产动态的重要指标,需要对其进行长期和系统的监测;(3)建议在今后的研究中,根据简化食物网的概念,对占生物量绝对多数的重要生物资源种类的营养级进行重点研究并采用国际通用的标准划分计算营养级。     

高碳铬铁EDXRF分析中的制样条件优化

利用能量色散X射线荧光光谱分析技术(EDXRF)检测高碳铬铁,能实现对主次组分的同时测量,且无需对铬铁样品进行化学处理.实验以荧光强度和测试精密度为指标,采用单因素试验考察了粉末压片制样过程中称样量、粘结剂种类与用量、研磨时间、压片压力、保压时间等制样条件对测量结果的影响趋势,并将优化后的条件用于对高碳铬铁的检测.结果...     

VLBI数据处理中的剩余钟行为与剩余大气效应建模

抽样选取了27次上海天台佘山站参与的天测与测地VLBI实验，分别进行了单次解算。通过分析解算参数随剩余钟行为和剩余大气效应分段拟合长度的变化，得到以下初步结论：（1）选取不同分段拟合长度时，站坐标解算结果和时延残差加权均方根存在差异。最大分别至百米级和数十皮秒，因而分段拟合长度不能随意选取。（2）分段拟合时段长度存在某一合理取值范围。它不宜过长，否则钟和大气的剩余效应短周期变化不能很好地模型化。为了保证待估参数解算时有足够的自由度，拟合时段不宜过短，融将导致法方程近于或出现奇异，达不到较好控制噪声的效果。（3）由于各次实验，同一次实验中的不同台站相应的钟和大气条件存在的差异。有必要对每次实验以及每次实验中的各观测台站分别分析。寻找合适的分段拟合长度，这在实际操作中显然相当烦琐。（4）一般而言，在剩余钟行为拟合长度缺省值60min情况下，剩余大气效应分段拟合长度以介于10min至40min为宜；在剩余大气效应分段拟合长度缺省值20min情况下。剩余钟行为分段拟合长度以介于20min至100min为宜。