黄、渤海盐度准三维数值预报模式

黄、渤海盐度的垂直结构具有典型的自模性,而其水平分布又受平流、水平扩散效应及径流等因素的影响。本文根据黄、渤海实测资料拟合了盐度垂直剖面的自模函数,并结合描述表、底层盐度及上均匀层厚度这三个特征量水平分布的方程,形成给出盐度三维结构的准三维模式。在模式中,综合考虑了海面风和热输入的强迫作用以及流场的平流、侧向混和及底层混和影响因素,同时还考虑了径流、蒸发及降水的作用,较客观地反映了盐度的三维分布及其变化的物理过程。试报结果分析表明,模式的功能较好,结果令人满意。     

层状饱和场地中圆柱形基础的阻抗函数

本文利用间接边界元法,建立层状饱和场地中三维土-基础动力相互作用模型,研究了孔隙中流体的存在、场地动力特性和圆柱形基础的高度对基础阻抗函数的影响。数值结果表明,孔隙中流体的存在对基础阻抗函数有明显的影响,尤其是水平和竖向分量;场地动力特性也对基础阻抗函数有明显的影响,层状饱和场地中阻抗函数围绕相应均匀半空间结果振荡,且基岩与土层刚度比越大或土层厚度越小,振荡的幅度越大,土层厚度越小,振荡的周期越小;另外,基础高度也对阻抗函数有明显的影响,随基础高度增大,阻抗函数逐渐增大。     

基于区域生长法和变差函数的PolSAR影像山脊线提取

山脊线通常是地震引发的山地灾害的起始部位。在监测此类山地灾害时,需要明确山脊线所处位置,从而为制定出相关预防措施提供支持。根据实际经验,由于合成孔径雷达(Synthetic Aperture Radar,SAR)的成像特殊性,在检测含建筑物的全极化SAR(Polarimetric Synthetic Aperture Radar,PolSAR)影像山脊线时,往往受到高散射强度的建筑物干扰,从而出现误判。针对PolSAR影像中山脊线受到建筑物干扰出现误识的问题,提出一种结合区域生长法和变差函数的识别方法。在识别过程中,首先,对采集的PolSAR影像通过区域生长法进行分割处理,再通过模糊C均值(Fuzzy C-Means,FCM)分类方法对变差纹理等特征进行聚类分析,并进行相似度判断,获取建筑物区域,在此基础上,对比建筑物、山脊线的混分成分,进而提取山脊线。通过对该方法进行实验对比分析,发现其相较于阈值分割方法提取精度有明显提高,可以为PolSAR影像中山脊线识别提供一种新思路。     

顾及有色噪声的GNSS时间序列时域信号提取

随着大型地震的发生,GNSS时间序列中除线性趋势和周期信号外,还存在大量震后瞬态,准确地提取各项时域信号是运用GNSS时间序列进行地学研究的关键.为此,本文提出了顾及有色噪声的GNSS时间序列时域信号提取法.该方法首先基于白噪声（White Noise,WN）+闪烁噪声（Flicker Noise,FN）模型,使用最大似然估计（Maximum Likelihood Estimation,MLE）对震前GNSS时间序列进行参数估计,并根据参数估值来去除震后时间序列中的震前信号,以此获取残差序列;然后将残差序列作为求解特征时间尺度的观测量,WN+FN模型作为观测量的随机模型,并采取非线性最小二乘法（Non-linear Least Squares,NLS）法估计特征时间尺度;最后利用估计的特征时间尺度构建GNSS时间序列函数模型,并采用MLE估计其未知参数,进而实现时域信号的提取.经模拟数据分析,考虑有色噪声时,特征时间尺度估计算法的收敛性提高了25%,各项未知参数的标准差（Standard Deviation,STD）显著下降.最后,将该算法应用于日本区域实测数据,并与传统方法进行了对比分析.

     

基于改进Harris函数的岩石统计损伤软化模型

为了反映岩石变形全过程,通过改进Harris函数,假设岩石微元强度服从改进的Harris概率密度分布,结合连续损伤理论,考虑损伤变量修正,建立岩石统计损伤软化模型。采用YSJ-01-00岩石三轴流变试验机进行红层泥岩常规三轴压缩试验,基于红层泥岩试验数据对模型参数进行敏感性分析,阐述模型参数的物理意义,分析累积损伤的扩展过程。引用相关文献中二长花岗岩和砂岩的试验数据,利用本文模型对3种岩石材料的试验数据进行辨识验证,给出初始损伤系数的约束范围,对比3种岩石材料的试验曲线和理论曲线,证明本文模型的合理性和适用性。     

远震体波波形反演获得的震源时间函数

地中海东部大地震成核问题的一个有趣的现象，是它们的震源时间函数的急始特性，使用同一反演方法确定了发生在世界其他一些地区破坏了地表的浅源大地震的震源时间函数都呈现出类似的特性，这些地震的断层滑动类型包括：正断层、走滑断层和冲断层型。这些急始型震源时间函数可能反映了破坏障碍体的困难性，即在破裂成核的动力断裂过程完成后，断层的局部地区仍未破坏，也可能反映了通常发生在大地震破裂期间的某些特殊的物理过程（如蠕变不稳定性）。发生在类似地质构造环境中，且震源机制也类似的地震，反映出共同的模式。     

广东省东部地区背景噪声面波群速度层析成像

选取广东省东部及福建和江西测震台网的26个固定台站2011年3～9月的地震背景噪声的垂直分量数据,用互相关方法获得瑞利面波的格林函数,再用频时分析方法提取瑞利面波的群速度频散曲线,最后用面波层析成像方法反演得到研究区周期为5～15s的瑞利面波群速度图像.反演结果揭示了广东省东部地区浅部地壳群速度结构存在明显横向不均匀性.短周期群速度异常与研究区内山脉、盆地的分布有较好的对应关系,沉积厚度较大的盆地地区表现为低速异常,而基底埋深较浅的山脉则显示了高速异常;地热分布对面波群速度产生重要影响,地热值较高的区域面波群速度一般表现为低速异常;粤东地区特别是粤东沿海地区地壳中部可能存在低速层.     

气候变化影响下极端水文事件的多变量统计模型研究

以黄河流域太原气象站和淮河流域鲁台子水文站为研究对象,利用Copula函数构建气候要素（降水）同极端水文事件（干旱和洪水）之间的多元统计模型,分析不同降水条件下不同等级干旱和洪水的发生概率变化。结果表明,Gumbel Copula能够较好地描述太原站7月份的前期累加降水量和帕默尔干旱指数（PDSI）的相关结构。随着降水量的增加,极端干旱的发生概率逐渐减小,重旱、中旱和轻旱的发生概率则先增加后减小。Clayton Copula能够较好地描述鲁台子水文站前期累加降水量和洪峰流量之间的相关结构。当前期累加降水量大于等于某一定值时,随着年最大洪峰x的增大,发生洪峰≥x的极端洪水事件的概率逐渐减小。在同一个极端洪水发生概率下,前期累加降水量越大,洪峰流量出现大值的可能性越大。     

Copula函数在多变量水文分析计算中的应用及研究进展

综述了多变量水文分析计算方法及存在的问题.系统地介绍了Copula函数的基本理论与方法,分析探讨Copula函数在多变量水文计算领域的适用性及优越性.总结了Copula函数应用的关键技术,包括参数估计、拟合检验及类型选择等,并展望了Copula函数的研究应用前景.     

利用接收函数研究郯庐断裂带鲁苏皖段及邻区地壳结构特征

利用郯庐断裂带鲁苏皖段及邻区的261个宽频带地震台站5年记录的763个远震波形数据,计算并筛选得到了10846条远震P波接收函数.采用P波接收函数H-κ法得到了该区的地壳厚度和壳内平均泊松比值,并采用共转换点叠加法进一步揭示台站下方Moho界面的起伏形态.研究发现：⑴研究区Moho界面埋深在27~40 km范围内变化,平均深度在~34 km,总体上以郯庐断裂带为界呈现出东薄西厚的特征.地壳厚度在不同块体之间或者是块体内部存在着明显差异,表明不同的地质构造单元经历了不同的构造演化过程.⑵研究区地壳泊松比在0.15~0.32之间变化,平均泊松比为0.24,略低于全球陆壳和中国陆壳平均泊松比值;然而,较大的泊松比浮动范围却意味着研究区内地壳物质具有强烈的横向非均匀性及物质组成的复杂性.沿郯庐断裂带展布着一条NNE-SSW方向的泊松比高值异常带,推测是镁铁质基性岩沿郯庐断裂带上涌至地壳所致,亦或是高温高压的幔源热物质底侵至下地壳所致.⑶研究区的地壳厚度和壳内平均泊松比存在着反相关的关系;地壳厚度和地表地貌特征呈镜向关系,即造山隆起区Moho界面埋藏较深,而平原盆地区Moho界面埋藏较浅.Moho埋深等值线分布图和研究区布格重力异常特征对应良好.⑷共转换点（Common Conversion Point,CCP）叠加法对Moho界面的刻画与H-κ叠加法求得的地壳厚度结果具有较好的一致性.CCP剖面表明郯庐断裂带不仅是扬子断块区和华北块体的分界断裂,更是一条切割Moho面、深抵上地幔的深大超壳断裂带,错距达4~7 km.⑸研究区内部分台站下方存在壳内分界面,仍能通过改变壳内速度后采用H-κ法获得其埋藏深度.