经验模态分解法在地震资料处理中的应用试验

通过对多种连续小波基的平稳序列时频分析,提出经验模态的非平稳序列的时频分析方法,并在地震资料处理解释应用中作探索性试验.试验结果表明：对地震记录道长为5 s和2 ms采样资料,经验模态分解后可获得6-7个固有模态函数单分量,分解的第一个分量C1（t）具有较宽的频率成分,可作为地震属性的主分量.     

基于宽范围岩石物性约束的大地电磁和地震联合反演

岩石物性约束是地球物理联合反演重要的实现方式.在以往联合反演研究中,这种约束通常以点对点的映射方式出现,要求经验性的岩石物性关联是相对精确的,无形中提高了岩石物性约束应用于联合反演的条件.基于此,本文结合全局寻优的模拟退火算法,提出了宽范围岩石物性约束技术.该技术将岩石物性约束与模拟退火模型扰动、模型接受准则相结合,实现了空间对空间的映射,将电阻率和速度耦合在一起.文章以MT和地震正则化同步联合反演为例,检验了不同的先验信息条件下的宽范围岩石物性约束技术的适用性.模型试验表明:新技术容易实现,容错性高,可降低"不精确"岩石物性关联信息引入带来的风险,提高岩石物性约束应用的灵活性.     

下垫面对WRF模式模拟黑河流域区域气候精度影响研究

利用高精度的土地覆盖、土壤质地类型和地形高度值替换了天气研究和预报模式Weather Research and Forecasting Model(WRF)中的相关数据,通过数值模式试验检验了下垫面数据对WRF模拟精度的影响。同时,通过与黑河综合遥感联合试验中7个测站观测数据的比较,以平均误差、均方根误差和相关系数为指标,分析了WRF模式下垫面数据改变对近地表气象要素的模拟精度的影响。结果表明:(1)WRF模式本身的地形高度信息在黑河流域上游地区有较大误差,造成了一定的模拟误差。而使用高精度的下垫面数据可以提高WRF模式在黑河流域上游复杂区域的模拟能力;(2)2m气温除了随地形高度递减外,还受土壤质地和土地覆盖小幅度影响,而且进行地形订正后的2m气温与2m湿度的模拟在下垫面为水体的区域对比强烈,因此为模式提供准确的水体分布信息也至关重要;(3)2m气温和湿度等要素的模拟差异值与地形高度资料的差异呈负相关,而降雨量的差异与地形高度差异呈微弱的正相关,与土壤质地差异和土地覆盖差异的相关性也比较弱。     

基于改进经验模态分解的地震弱信号增强处理方法

针对地震反射信号经验模态(EMD)分解过程中出现的混叠及分解重建失真等问题,本文提出一种改进的EMD分解方法.首先对地震信号进行小波转换,对转换分解得到的信号高频部分进行噪声压制处理,然后将转换后的结果再进行EMD分解;针对频带宽度、振幅大小变换,提出了地震信号EMD分解的分段控制停止准则,根据主频带宽,确定相应的时窗长度,并采用各个时窗段部分叠置衔接的方法,最后对每段结果再利用单调函数停止准则进行变换,该方法有效避免了计算结果的时窗效应.通过实际资料处理对比分析,改进的EMD分解方法较以往常规EMD分解方法,在分解结果的可靠性、分解精度及非线性弱信号时频分析能力等方面都具有显著的改善,为提高薄层地质体识别能力提供了有力的技术手段.     

基于希尔伯特-黄变换的九寨沟M7.0地震加速度记录时频分析

2017年8月8日四川九寨沟M7.0地震是继2008年汶川M8.0地震和2013年芦山M7.0地震之后,青藏高原东缘在不到10年的时间内发生的第3个震级M7.0以上的强震,震中位于青藏高原巴颜喀拉块体东缘东昆仑断裂带东端的塔藏断裂、岷江断裂和虎牙断裂交汇部位,四川省地震局的数字强震台网共有37个台站获取了主震的三分量强震加速度记录。由于傅里叶（Fourier）变换仅能提供强震记录的频域信息,故本文在对九寨沟M7.0地震的加速度记录进行时频分析时采用了一种基于聚类经验模态分解（EEMD）的希尔伯特-黄变换（HHT）方法提取信号时频特性,通过对震中附近台站的加速度记录进行EEMD分解和希尔伯特（Hilbert）变换及谱分析,最终有效获得了信号能量的时频分布特征,量化提取了中心频率、Hilbert能量、最大振幅对应的时间等特性,并与Fourier变换进行了对比研究。研究结果表明:对于非线性的强震记录采用EEMD能够有效抑制经验模态分解（EMD）中存在的模态混叠问题,FFT谱与Hilbert边际谱相比,它在低频处会低估地震动的幅值,随着频率的增加,FFT谱又会放大其幅值。      

基于正交设计的复杂坝基弹塑性力学参数反演

复杂坝基弹塑性力学参数反演的计算工作量直接取决于采用何种优化方法。根据大坝在运行期的实测资料，运用有限元正交数值试验、回归分析和优化相结合的方法，反演了坝基岩体的弹塑性力学参数。实例研究表明了该方法的有效性和实用性。     

格陵兰冰芯氧同位素显示近千年气候变化的多尺度分析

用经验模态分解(EMD)方法对格陵兰冰盖GISP2冰芯古气候代用指标δ~ (18) O 序列进行分解, 结果表明: 格陵兰近10 ka来气候变化总趋势出现历时约490 a的中世纪暖期和历时约570 a的小冰期, 其间还存在次级的冷暖期变化;气候波动具有3 a、 6.5 a、 12 a、 24 a、 49 a、 96 a、 213 a、 468 a准周期, 既有NESO的因素, 也受制于太阳活动周期的影响. 第7, 第8内在模函数(IMF7, IMF8)波动振幅以及总体趋势分量res在1350A.D.出现明显的转型, 表明1350A.D.为中世纪暖期和小冰期的分界.     

卫星遥感反演土壤水分研究综述

土壤水分是影响地表过程的核心变量之一。精准地测量土壤水分及其时空分布,长期以来是定量遥感研究领域的难点问题。简要回顾基于光学、被动微波、主动微波和多传感器联合反演等卫星遥感反演土壤水分的主要反演算法、存在的难点和前沿性研究问题,介绍了应用土壤水分反演算法所形成的3种主要全球土壤水分数据集,包括欧洲气象业务卫星（ERS/MetOp）数据集、高级微波扫描辐射计（AMSR-E）数据集、土壤湿度与海洋盐分卫星（SMOS）数据集,并结合目前存在的问题探讨卫星遥感反演土壤水分研究的发展趋势。     

矿柱稳定性影响因素敏感性分析及其应用研究

为研究矿柱稳定性影响因素的敏感程度,以招金大尹格庄金矿为研究背景,从矿柱载荷、强度、失稳形式及影响因素确定四方面推导出两种矿柱形式的安全系数计算公式,设计采用6因素5水平正交试验对矿柱稳定性影响因素敏感性进行分析,并研究主要影响因素与矿柱安全系数之间的函数关系。研究表明,矿柱宽度、矿体开采深度和矿房宽度3个主要因素对矿柱稳定性影响程度是最为剧烈的;矿柱安全系数随着矿柱宽度地增加呈指数形式递增且递增速率不断加大,随矿体开采深度地增加呈幂函数形式递减且递减速率逐渐变缓,随矿房宽度地增加呈负指数形式递减且降低速率逐渐变慢。利用DPS和Matlab建立的矿柱安全系数回归方程得出保证安全生产的合理矿柱尺寸,矿房宽度应不超过8 m,布置条形矿柱宽度应不小于3.6 m,布置方形矿柱宽度应不小于5.9 m。结合采场矿柱布置实际,应用效果较为理想,可为下中段矿体回采矿块布置提供较好的依据。     

水文双变量重现期分析及在干旱中应用

单变量水文统计中一些广为接受的概念在多变量环境下尚缺乏深入分析,也易被误解,如N年内重现期大于等于T的多变量事件发生的次数与N/T的关系。实践中,多变量联合重现期与其边缘分布变量重现期的一些经验关系被发现并通过了案例验证分析,但缺乏解释和推导。基于GH Copula推导了双变量联合重现期与边缘分布变量重现期的关系以及双变量事件发生次数与其重现期、变量相关程度间的定量关系。以昆明56年的逐月SPI(Standardized Precipitation Index)和SRI(Standardized Runoff Index)识别了干旱事件,采用GH Copula构建了干旱历时和烈度的联合分布函数,验证了双变量联合重现期与边缘分布变量重现期的关系以及多变量事件发生次数与其重现期的定量关系。表明不宜以“and”第1重现期是否接近于比该干旱事件的旱情更重的干旱发生的平均时间间隔来说明干旱特征值重现期分析的合理性。变量的相关性不强时,需谨慎采用边缘分布变量重现期的较大值近似代替“and”事件的第1重现期。