长江三角洲地区TZK3孔碎屑锆石U-Pb年龄及其物源意义

以TZK3孔的磁性地层学为基础,结合沉积物的岩性、结构构造及锆石年龄谱系特征,探讨沉积物的物质来源及长江贯通的时限。古地磁结果表明:96.7 m、263.3 m、603.75 m分别对应B/M界线、M/G界线、晚上新世/早上新世的界线（3.58 Ma）。锆石年龄谱系数据显示,TZK3孔的U-Pb锆石年龄主要分为5组:100~300 Ma,400~500 Ma,700~850 Ma,1800~2000 Ma,2400~2600 Ma。其中3.7 Ma的锆石年龄谱相对简单,以白垩纪（100~150 Ma）为主,物质主要来自于长江中下游的火山盆地,为近源沉积。TZK3孔3.04 Ma以来,锆石年龄谱变得复杂且主峰相对较多,表明物源区更广且加入了远源的成分。3.04 Ma的锆石年龄谱中开始出现峨眉山玄武岩年龄段（251~260 Ma）的锆石,表明在此时期长江上游的物质就已到达了长江三角洲地区,即长江贯通的时限为3.04~3.7 Ma。      

冬季大风影响下的渤黄海水交换特征

利用ROMS海洋数值模式对2006年冬季渤黄海的海洋动力环境进行模拟,基于温度、盐度模拟结果,使用谱混合模型进行水团分析,定义了渤海海峡地区的水交换区。并进一步讨论了冬季大风事件对水交换区的影响,给出了冬季大风影响下的渤黄海水交换特征。研究得出,冬季的黄海水团以“舌”形分布于渤海海峡地区,水交换区则表现为沿“舌”形边缘呈带状分布,具有西北——东南的走向趋势,并且在“舌”尖处的水交换面积最大。通过缩小研究范围,发现位于黄海最北部的沿岸海域并不参与渤黄海之间的水体交换。最后研究发现,冬季大风事件对渤海水交换具有促进作用,具体表现为：大风过程使黄海暖流对渤海的入侵更加深入,水交换区向渤海方向伸展,南部的水交换带变宽,河流径流进入渤海后与渤海水的混合区加大,并发生北移。     

地球自由振荡的保结构模拟

大地震等诸多激励均能激发全球自由振荡现象,通常表现为驻波形式的全球整体振荡.现有的地震波数值模拟方法多为非保结构方法,无法压制长时程计算中的积累误差.本文采用优化的三阶辛格式谱元法,对地球自由振荡及全球尺度的地震波传播进行了长时程模拟.通过与传统的基于Newmark算法的谱元方法结果对比分析,明确验证了本文所得优化的三阶辛格式谱元法在模拟地球自由振荡等大规模长时程问题上的优越性和准确性.上述进展在方法论层面为今后探测、刻画全球尺度地球非均匀结构的驻波数值方法奠定了部分基础,并为相关研究领域提供了新的选择.     

位场各阶垂向导数换算的新正则化方法

位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.     

大功率电波加热电离层中热自聚焦不稳定性的理论研究和数值模拟

本文首先从电子密度及电子温度的输运方程和考虑自作用时的电磁波波动方程出发,利用简正模展开的方法推导出泵波在反射区域激发出热自聚焦不稳定性(thermal self-focusing instabilities,TSFI)所需电场阈值以及其增长率的完整数学表达式,并估算了TSFI激发阈值及所对应的有效辐射功率(ERP)的量级.随后利用三维垂直加热的理论模型,结合国际参考电离层(IRI-2012)和中性大气模型(MSIS-E-00)给出的背景参数,数值模拟了大功率高频泵波加热电离层时泵波反射区域电子密度及电子温度因TSFI而产生的变化及发展的过程,并对比分析了不同背景参数对较热效果的影响.结果表明:当高频泵波的加热阈值达到或超过百毫伏每米的量级时,即可激发TSFI,发展出大尺度电子密度及温度不均匀体,这些不均匀体内的密度耗空约为4%~10%,而电子温度剧烈增长,到达背景温度值的1.6~2.1倍;且在相当的加热条件下,背景电子温度越低、电子密度越小,加热效果越显著;电子密度及电子温度的扰动幅度随着加热时间的推移而逐渐减小,即扰动逐渐趋于饱和,且电子温度要快于电子密度达到饱和状态.本文还对泵波反射高度处的电子密度及电子温度变化率进行采样并求得其功率谱密度,分析结果表明:TSFI发展出的大尺度不均匀体满足幂律谱结构,谱指数随着加热的进行逐渐趋于稳定,白天与夜间的幂律谱指数区别不大,但电子密度与电子温度的幂律谱有所区别.     

山脊线与坡度和峰值速度放大系数的相关性研究

地形效应是地震工程和地震学的重要研究内容.数值方法是研究地形效应的重要工具.然而以往关于地形特征对地形效应的影响的研究大多是基于二维简单模型.对三维真实地形特征和地震动放大的关系缺乏详细的研究.为了更好地描述三维地形特征,本文将地形用不同的地形特征因子来描述,并用算法将它们从地形数据中提取出来.从而使地形和放大系数的相关性的研究转换为不同地形特征因子和放大系数的相关性的研究上.本文中,我们选择芦山地区作为研究对象,用算法提取出特征因子山脊线和坡度.它们分别表征了山脊形态和山体的陡峭程度.我们又用谱元法分别模拟三种不同主频Ricker子波的地震波在起伏地表和水平地表模型中的传播,并得到了它们各自的峰值速度(PGV)放大系数.通过分析三个主频的PGV放大系数在山脊线上分布,我们发现放大系数幅值在山脊线上分布不均匀.放大系数较高的地方位于山脊线弯曲或宽度发生变化的附近.我们又统计了三个放大系数在滑动平均坡度上分布,并得到了放大系数均值和坡度的正相关曲线.结果表明坡度和放大系数两者在幅值分布上具有正相关性.     

基于FFT-MA谱模拟的快速随机反演方法研究

虽然基于地质统计学的随机反演方法能够有效融合测井资料中的高频信息,但计算效率低,占用内存大,限制了它在实际资料中的应用.本文在保留传统随机反演方法优点的基础上,创造性地引入傅里叶滑动平均(Fast Fourier Transform-Moving Average,FFT-MA)谱模拟进行频率域的地质统计模拟,并利用逐步变形算法(Gradual Deformation Method,GDM)确保模拟结果与实际地震数据的匹配,构建了基于FFT-MA谱模拟的新的快速随机反演方法.与常规随机反演相比,新方法不仅分辨率高,而且能够使反演解得到快速收敛,有效提高计算效率,减少内存占用.模型试算获得了与理论模型吻合度较好的高分辨率反演结果.实际资料分析也表明新方法所得到的高分辨率反演结果能够对薄互储层进行良好的展示,为薄储层的识别提供高效可靠的技术支持.     

GNMF小波谱分离在地震勘探噪声压制中的应用

地震勘探资料噪声压制及信噪比提高是整个地震勘探信号处理过程中的重要任务,随着地震勘探深度的增加及其复杂性,人们对地震数据质量的要求越来越高.勘探环境的复杂化使得采集到的地震资料中有效信号被大量噪声淹没,无法清晰辨识,严重影响后续的数据处理与解释.小波去噪是地震勘探中常用且发展较成熟的一种方法,但是其涉及到的阈值函数选取问题一直令人困扰,虽然已有多种阈值函数被提出,但仍存在各自的缺陷.本文利用小波分解在时域及频域良好的信号细节体现特性,引入模式识别中的非负矩阵分解(NMF)谱分离思想,针对小波系数阈值优化问题,提出了一种小波域图非负矩阵分解(GNMF)消噪算法.该方法首先在小波分解基础上,利用GNMF算法实现小波分解系数谱中信号分量与噪声分量的谱分离,然后通过反变换重构各分离子谱对应的子信号,最后利用K均值聚类算法将得到的多个子信号划分为信号类及噪声类,最终得到重构信号及分离噪声.合成记录和实际地震资料的消噪结果验证了新方法在提高信号与噪声分离准确性和精度方面的有效性,同时新方法避免了阈值选取造成的噪声压制不理想或有效成分损失问题.与小波消噪结果的对比及数值分析也说明了新方法在噪声压制及有效成分保持方面的优势.     

不同仪器烈度算法在四川地区 历次地震中的比较应用

强震观测资料包括强震动时程记录资料和宏观震害调查资料两大类, 前者是定量的微观数据, 后者则是定性的宏观指标. 本文明晰给出了微观仪器烈度与宏观仪器烈度、 微观考察烈度与宏观考察烈度的概念, 并在此基础上提出了只有在宏观仪器烈度与宏观考察烈度之间进行比较才具有实质性的意义;同时对国内外常见的仪器烈度算法进行了比较性研究, 得出了袁一凡仪器烈度算法可靠性更高的结论;最后以四川地区历年来重要震例的强震动记录为依据, 对修正的袁一凡仪器烈度算法的可靠性进行了比较应用. 结果表明, 对袁一凡仪器烈度算法修正与扩展的应用是可行的, 同时也验证了将微观仪器烈度与宏观考察烈度直接进行比较会存在较大的差距.      

基于强震和测震数据对鲁甸6.5级地震发震断层方向的研究

2014年8月3日16时30分,云南省昭通市鲁甸县发生M6.5地震,宏观震中位于鲁甸县龙头山镇(27.1°N,103.3°E),深度约12km。发震断层为NW向包谷垴-小河断裂,是NE向昭通-鲁甸断裂系的NW向次级走滑断层。文中利用鲁甸地震50个强震台站观测记录与135个测震台站观测记录,分别进行峰值加速度计算和S波辐射图案匹配,对此次地震发震断层方向进行讨论。结果表明,"NW向"的断层方向更符合强震、测震仪器观测值,同时与地震灾害现场调查结果一致。