基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

下扬子及周边地区深部泊松比结构及深部动力过程约束研究

下扬子及周边地区存在着丰富的地质构造和多金属矿产资源,其深部结构和动力学过程已成为地学界研究的热点.为了更好地讨论下扬子及周边地区的深部动力学过程和岩浆活动机制,本文基于均匀网格层析成像方法提出了非均匀网格远震层析成像方法,利用大量的天然地震相对走时残差数据反演获得了下扬子及周边地区深至700 km范围内的三维S波和P波速度结构,并根据纵横波的比值关系计算出泊松比异常.由于S波速度比P波对流体的反应更加敏感,所以泊松比异常反映了物质是否包含流体或者物质的软硬、冷热程度.本研究结果显示：（1）长江中下游成矿带下方的上地幔内存在明显的高泊松比异常,而地幔过渡带内则存在明显的低泊松比异常;（2）大别造山带及其南侧的中扬子地块的上地幔中下部及地幔过渡带内都存在明显低泊松比异常,且呈现东深西浅的空间分布特征.结合已有的地质、地球物理和地球化学等资料,我们认为长江中下游成矿带下方的地幔过渡带内滞留着古太平洋俯冲板块,其上地幔内则赋存着软的上地幔热物质,为深部成矿提供了热量或幔源物质.因此,古太平洋板块的俯冲对长江中下游成矿带的形成发挥了至关重要的作用.     

多通道瞬变电磁法频率域比值算法及视电阻率计算研究

视电阻率在电磁法勘探中起重要作用,可以作为初步解释结果,多通道瞬变电磁法通过大地脉冲响应的峰值时刻定义视电阻率.鉴于传统方法计算峰值时刻视电阻率需要恢复完整大地脉冲响应,本文通过研究大地脉冲响应及其频谱特征,提出一种新的频率域比值法.利用参考谱及其与一般谱的伸缩因子计算一般谱的峰值时刻,进而得到视电阻率.该方法只需大地脉冲响应单个频点幅度谱,从而不用估计整个大地脉冲响应.实际计算时,可以用大地脉冲响应多个频点幅度谱得到平均视电阻率.模拟数据与实测数据结果表明,在一维介质中,频率域比值法可以得到合理的视电阻率.     

Selection between the generalized Pareto and kappa distributions in peaks-over-threshold hydrological frequency modelling

Hydrologists use the generalized Pareto (GP) distribution in peaks-over-threshold (POT) modelling of extremes. A model with similar uses is the two-parameter kappa (KAP) distribution. KAP has had fewer hydrological applications than GP, but some studies have shown it to merit wider use. The problem of choosing between GP and KAP arises quite often in frequency analyses. This study, by comparing some discrimination methods between these two models, aims to show which method(s) is (are) recommended. Three specific methods are considered: one uses the Anderson-Darling goodness-of-fit (GoF) statistic, another uses the ratio of maximized likelihood (closely related to the Akaike information criterion and the Bayesian information criterion), and the third employs a normality transformation followed by application of the Shapiro-Wilk statistic. We show this last method to be the most recommendable, due to its advantages with sizes typically encountered in hydrology. We apply the simulation results to some flood POT datasets.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR E. Volpi     

宾川地区的地壳厚度和泊松比分布特征研究

This paper presents the changes of crust thickness and Poissons ratios distribution in the Binchuan region,where the first air-gun transmitting station and its a small dense array were deployed.From September 2011 to January 2014,more than 239 teleseismic events of M≥6.0 were recorded in 16 stations in the Binchuan region.Their P-wave receiver functions were analyzed respectively.The high spatial resolution result shows that the average crust thickness of Binchuan region is 45.3km,it follows the rule of"deeper in the north and east part,shallower in the south and west part."The deepest region is in Xiaoyindian Station;the crust thickness is 47.9km;the shallowest region is in Paiying Station,it has the thickness of 42.1km.It shows that the deeper Moho surface nearby the Chenghai fault and shallower nearby the Honghe fault;the isoline distribution of thickness changes greatly nearby the Chenghai fault and slowly nearby the Honghe fault.From the distribution of Poissons ratios,it is unevenly in the study area with a great difference from the north part to the south part,which shows a characteristic of"lower in the north,higher in the south".The Poissons ratio nearby the Honghe fault is medium too high(0.26≤σ≤0.29);lower nearby the Chenghai fault(σ≤0.26).This paper concludes the possible reason of different characteristic between Poissons ratio and crust thickness is thicker in the upper crust in the Binchuan region.     

新疆天山地区的波速比异常分析

Based on the seismic observation report data provided by the Xinjiang Digital Seismic Network from 2009 to 2014,we calculate the wave velocity ratio and its background value for medium and small earthquakes by using the multi-station method in Tianshan,Xinjiang.This paper analyzes the variation of the wave velocity ratio disturbance value to highlight the abnormal,and also back-traces 7 moderate earthquakes at the research area.The results show that:(1)the background value of the wave velocity ratio is almost 1.70,the wave velocity ratio obviously decreases in the middle-eastern part of Tianshan and the region near the Puchang fault;(2)the wave velocity ratio disturbance value is mostly low in the epicenter before four earthquakes of M≥5.0 from 2011 to 2013 in the study area;(3)before 7 moderate strong earthquakes,the earthquake events with low value of the wave velocity ratio account for over 60% of corresponding total events near the epicenters,and the low value of the wave velocity ratio is relatively obvious before moderate earthquakes.     

月球辐射状撞击坑经度方向分布特征与同步旋转轨道状态模拟分析

月球目前的同步旋转轨道状态使得形成的撞击坑分布满足一定的不对称性.本文利用最新的LRO影像和地形数据,结合早期的Clementine影响数据,分析了月球辐射状撞击坑经度方向成坑率分布,结果表明所识别的辐射状撞击坑的.年龄为O.9 Ga内,西一东半球方向存在明显不对称性,比值约为1.35~1.53;同时利用辐射状撞击坑模拟了0.9 Ga内的同步旋转轨道的平均状态,得到该时间内月球的向点-背点为70°W—110°E附近;最后利用撞击坑数据库资料对月球形成以来各地质世纪时间尺度内的撞击坑进行了成坑率分布分析,结论表明月球在大爆炸中后期间可能处于过近似的同步轨道旋转状态,但是与现在的轨道运行状态相反,月球在之后的地质时期内经历了翻转,其诱因可能是雨海和东方海盆地遭受的撞击.     

青藏高原东缘的地壳流及动力过程

黏滞性地壳流对地壳及上地幔变形作用及动力机制,是大陆新生代造山带的一个重要研究内容.青藏高原中下地壳存在部分熔融或含水物质的黏滞性流体,已为一系列地球物理及岩石学研究所证实.为研究青藏高原东缘地壳流的动力作用,本文用密集的被动源宽频带地震台的观测数据,反演了地壳上地幔精细速度结构和泊松比.研究表明,川西及滇西北高原的中地壳内普遍存在低速层,而高泊松比的地壳只分布在川西北地区.位于中地壳的黏滞性地壳流从青藏高原腹地羌塘高原流出,自北西向南东流入青藏高原东缘.这些黏滞性地壳流带动了上地壳块体水平移动,当它们受到刚强的四川盆地及华南地块阻挡时将发生分层作用,地壳流将分为二或更多分支不同方向的分流,向上的一支地壳流将对上地壳产生挤压,引起地面隆升,向下的一支地壳流将使莫霍面下沉加厚下地壳·黏滞性地壳流的运动在地壳中产生应变破裂发生强烈地震活动,地震的空间分布与震源机制也受到地壳流动力作用控制.     

基于动态图像变形的PP与PS波层位直接匹配

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分:首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明:该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.     

中国东北地区高分辨率地壳结构:远震接收函数

利用分布在东北地区的国家地震局台网、NECESSArray台网、吉林大学在长白山及其周边地区布设的26个临时台站总计259个台站接收到的16,070条高质量的P波接收函数,采用H-k和CCP(Common Conversion Point,共转换点)叠加成像方法,获得该区高分辨率的地壳结构.观测结果显示,东北地区莫霍界面深度和地表高程总体呈镜像关系;西部大兴安岭—太行山重力梯级带附近存在莫霍界面深度陡变带;中部的松辽盆地地区受晚中生代的地壳拉伸作用影响,地壳厚度较薄,北部的小兴安岭地区和南部的华北北缘造山带可能同样受拉伸运动影响,具有较小的地壳厚度;松辽盆地莫霍界面深度由西向东逐渐减小,推测这与太平洋板块俯冲作用有关;东部地区莫霍界面结构比较复杂,依兰—伊通断裂与敦化—密山断裂之间出现复杂震相,可能与该区存在地幔物质的底侵作用有关;长白山火山地区地壳厚度较大,对应较高的波速比,推测在该区地壳内存在岩浆囊.