基于局部采样MCMC方法的时移探地雷达反演

马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...     

平流及层结在阿尔卑湖泊风成扩散问题中的作用

A semi-implicit semi-spectral hydrodynamic primitive equation model is used to study tracer diffusion in homogeneous and stratified lakes. Impulsively applied spatially uniform wind is applied in the long direction of a rectangular basin with constant depth and Lake Constance. Tracer mass is released at various locations of the free surface and its spreading under the action of dijfusion and advection is studied. We show that for homogeneous basins the tracer quickly spreads over the entire water depth. Inertial waves can also be detected in the tracer concentration. On the other hand for stratified waters the tracer mass is largely confined to the epilimnion with occasional penetration into the hypolimnion where large downwelling occurs. Here Kelvinand Poincar ''e-type wave dynamics is discernable in the tracer-concentratwn-time series, the former more conspicuously at nearshore, the latter at off-shore positions.     

含流体致密砂岩的纵波频散及衰减:基于双重双重孔隙结构模型描述的特征分析

致密砂岩普遍具有低孔、低渗及微裂缝发育的地质特征,并且呈现出很强的非均匀性.致密砂岩储层与常规砂岩储层比较,具有明显的岩石物理性质、渗流力学性质方面的差异.致密砂岩内部的非均匀性对弹性波频散、衰减有显著影响,其中包括孔隙结构的非均匀性,即岩石内部孔隙参数的不均一性,以及孔隙内部不相混溶流体的非均匀分布;此外,非均匀性的尺度也决定了波出现显著频散与衰减的频段.综合考虑致密砂岩孔隙结构非均匀性及流体斑块状饱和的非均匀性,本文采用双双重孔隙介质结构模拟了致密砂岩的弹性波响应,分析了同时具备两类非均质性岩石中的波传播特征.调查分析了两组分别来自中国鄂尔多斯盆地苏里格气田及四川盆地广安气田的不同类型致密砂岩储层的岩芯超声波实验数据,给出了岩石样本的弹性波速度频散与衰减曲线.结果显示理论模型预测结果与完全饱和、部分饱和岩石的实验数据吻合良好.对两个地区致密砂岩岩芯数据进行对比分析,苏里格致密砂岩样本总体上比广安致密砂岩渗透率高,在各孔隙度范围内,特征模拟显示苏里格样本的裂隙尺寸明显大于广安样本.广安致密砂岩在低孔隙度范围内发育了更多、更小的颗粒裂隙/接触.致密砂岩的速度频散与衰减结果受流体黏度、晶体破裂及流体斑块状饱和的共同影响.此外,孔隙度越大,部分饱和岩石中斑块状饱和机制对总衰减的贡献越低,与之相对,结构非均质性所占的比重则有所增强.     

基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析

基于青藏高原东北缘密集宽频带野外流动观测台阵以及固定台站资料,利用双差层析成像方法对地震位置和研究区的地壳速度结构进行了反演.最终用于联合反演的地震事件合计9644个.结果显示青藏高原东北缘速度结构具有明显的横向不均匀性.从整体上看,青藏高原地区表现为低速异常,鄂尔多斯表现为高速异常,而扬子地块亦表现为高速异常.不同深度处速度结构表现不一致,同一深度处P波速度结构和S波速度结构也有明显差异.由西秦岭北缘断裂带、临潭-宕昌断裂以及礼县-罗家堡断裂围限的地震活动强烈的区域中,P波速度结构由深度0 km时呈现的低速异常,逐渐过渡到5 km时高低速相间分布的特征;而S波速度结构在此区域中,由近地表0 km时高低速相间分布的特征,逐渐过渡到30 km时几乎表现为低速异常.2017年8月8日九寨沟7级地震所在的塔藏断裂、岷江断裂和雪山断裂围限区域,在深度20 km处的P波速度结构和周围存在明显差异,九寨沟地震处于高速异常与低速异常的过渡带内.此外,2013年7月22日发生在青藏高原东北缘的岷漳县6.6级地震,震源区所在的临潭-宕昌断裂附近的P波速度结构在15 km深度处也有明显特征,震源位置所在区域也处于高低速过渡带.该区域这种地壳内部高低速过渡带可能是应力比较容易积累而发生中强地震的一个重要场所.     

3D P-WAVE VELOCITY STRUCTURE AT THE NORTHEASTERN MARGIN OF ORDOS BLOCK

Using the 7 100 absolute first arrivals of P waves and 91 513 relative P arrival times of 726 events at the northeastern margin of the Ordos block since 2009, the 3D fine structure of P wave velocity within the depth of 15km in the crust was inverted by the double difference seismic tomography method. The results show that there exist obvious high-speed continuous bodies in the northwest of the study area, and their lateral areas increase gradually with depth, while the velocity of east and south is relatively low. The velocity inhomogeneity exists and differs at different depths. The lateral differences of velocity are related to seismicity and faults. The 5~15km depth profile shows that earthquakes tend to occur in the area with relatively high velocity or high speed transition zones, which to some extent reflects the fragility of regional crustal media and the strong differential movement of faults in vertical and horizontal directions where the crust body is easy to absorb and store strain energy and generate major earthquakes. A "Y"-shape low-velocity channel is present in the lower crust around Liangcheng, corresponding to the NW-trending Heilaoyao-Shahukou fault set, which may reveal the migration path of the Late Tertiary-Quaternary basalt eruption. The Helingeer M6.2 earthquake in 1976 was related to the formation of the locking section of the thermal welding in this area. The three-dimensional fine structure of P wave velocity presented in this paper provides intuitive seismological evidence for physical and chemical properties of crustal media and the deep tectonic environment of earthquake preparation.     

双稀疏字典和FISTA的地震数据去噪

地震数据的随机噪声去除是地震数据处理中的一项重要步骤,双稀疏字典提供了两层稀疏模型,比单层稀疏模型可以更好地去除噪声.该方法首先利用contourlet变换对地震数据进行稀疏表示,然后在contourlet域中使用快速迭代收缩阈值算法（fast iterative shrinkage-thresholding algorithm,FISTA）对初始字典系数进行更新,接着采用数据驱动紧标架（data-driven tight frame,DDTF）在contourlet域中得到DDTF字典并通过FISTA得到更新后的字典系数,最后通过DDTF字典和更新后的字典系数获得新的contourlet系数,并对新的contourlet系数进行硬阈值和contourlet反变换得到去噪后的数据.通过模拟数据和实际数据的实验证明：与固定基变换去噪方法相比,该方法可以自适应地对地震数据进行稀疏表示,在地震数据较为复杂时得到更高的信噪比;与字典学习去噪方法相比,该方法不仅拥有较快的去噪速度,而且克服了字典学习因为缺少先验约束造成瑕疵的缺点.     

龙门山断裂带精细速度结构的双差层析成像研究

利用川西流动地震台阵、汶川地震震后应急台网记录到的P波到时资料,对2008年5月至2008年10月期间发生的汶川地震余震序列应用双差层析成像方法进行了地震震源和三维P波速度结构的联合反演.结果显示,联合反演获得的地震重定位结果与基于一维地壳参考模型的双差定位方法结果相近;研究区15 km以上速度结构与地表断裂分布密切相关,20 km以下深度呈现北东向和北西向交错结构.汶川地震破裂带南段龙门山断裂带之间上地壳呈现高速异常,速度结构的非均匀变化是控制余震分布和主震破裂传播的主要因素;联合反演结果给出了小鱼洞-理县方向存在隐伏断裂的速度结构证据,同时发现,破裂带北东段可能沿新发断裂扩展;结果确认了汶川地震起始段的高角度逆冲断裂特征,也确认了前山断裂和中央断裂在约20 km深度合并到脆韧转换带的特征.     

热扩散驱动的元素分异和同位素分馏:一种不容忽视的硅酸盐成分分异机制

热扩散在地质过程中是否发挥重要作用一直存有争议。本文回顾了热扩散的研究历史和现状,重点总结了热扩散驱动的元素和同位素行为规律,并探讨了温度、硅酸盐组分、压力和氧逸度等因素对热扩散行为的影响。已有的研究表明,稳定热梯度下的硅酸盐热扩散效应类似于结晶分异或AFC过程,可以造成轻、重同位素分别在高温端和低温端富集,而主、微量元素的扩散方向则取决于两端化学势的高低和熔体中的电价平衡。从基性岩浆到酸性岩浆,熔体聚合度增大,黏度增加,热扩散速率明显降低,成网元素的热扩散效应减弱,变网元素则反之;水、氟、氯和硫化氢等挥发组分能增加熔体的非桥氧比例,降低熔体聚合度,因而能显著增强硅酸盐熔体中元素和同位素的热扩散效应。在此基础上,本文提出了当前硅酸盐体系热扩散研究中存在的五个亟需解决的问题,即:1)对不同硅酸盐体系的热扩散规律的研究还不够全面;2)对微量元素的热扩散行为认识不足;3)硅酸盐体系热扩散作用的影响因素及尺度还不够明确;4)热扩散作用的地质与地球化学关键识别标志有待确立;5)硅酸盐体系热扩散作用的理论模型有待建立。尽管硅酸盐体系热扩散的研究还存在诸多不足,但越来越多的证据表明,热扩散是地质过程中不容忽视的一种成分分异机制。这种机制会造成岩浆房或岩浆通道中的元素分异和同位素分馏,可能对于一些成分分异的岩石和矿床的形成具有重要的意义。     

水动力条件对沉积物-水界面氧通量的影响

氧环境决定了水体沉积物中各种生命所需元素的最终归趋,沉积物-水界面是水相与沉积物相氧传递的重要场所,而水动力条件是影响沉积物-水界面氧传递的重要因素.选择三峡库区一级支流御临河为研究对象,根据长年监测数据建立实验室模型,采用声学多普勒流速测试仪及微电极测试系统构建了非侵入式涡度相关测试系统,探究了不同水动力条件对沉积物-水界面氧通量的影响.结果表明：水体静止状态下沉积物-水界面溶解氧浓度随时间的增加而减少,非静止状态下随时间的增加而增加;沉积物-水界面氧通量随水体流速的增加而增加.根据氧通量求解对应流速下垂直涡动扩散系数并进行线性拟合,当水体流速为0.01~0.14 m/s时,垂直涡动扩散系数与水体流速的相关性最好,此时沉积物-水界面氧通量的传递以涡动扩散为主导.     

2017年9月4日河北临城M_L4.4地震震源深度测定

震源深度是地震定位中的重要参数之一。对于一个地震,sPn与Pn震相到时差是一个常量,与震中距无关,只与震源深度和地壳模型有关。基于该理论,以河北数字台网数据为基础,识别sPn震相和Pn震相,计算2017年9月4日河北临城ML4.4地震震源深度,结果为6.86km。双差定位法是目前确定震源深度比较精确的方法,双差定位得到的震源深度为7.01km。双差定位法所用数据量比较大,计算过程繁杂。利用sPn与Pn到时差计算震源深度的结果与双差定位法结果接近,计算快速、简单,所需数据量小。因此,将该方法应用到实际工作中,可以快速、准确地测定近震震源深度。