短周期地震仪接收函数的可行性分析——以新疆和田地震台阵为例

利用新疆和田地震台阵3 km孔径范围内架设的9个子台(包括1套宽频带和9套短周期地震仪)记录的3年远震波形数据,对比研究和分析了短周期地震仪接收函数的稳定性和可靠性.通过比较和分析短周期和长周期地震仪获取的接收函数波形,结果发现:(1)短周期地震仪记录与宽频带地震仪记录得到的接收函数有很好的一致性,且具有较高的稳定性,无论α取值为1.5还是2.5,短周期地震仪接收函数与宽频带地震仪接收函数都具有较强的线性相关性(相关系数0.9),但Ps震相均存在小幅的振幅差(约20%);(2)采用接收函数震相到时信息的方法(如H-κ叠加搜索),短周期地震仪可以代替宽频带地震仪;(3)由于短周期地震仪缺乏0.155 Hz以下的低频信号和在1 Hz以下频段非线性的振幅响应,仅仅采用短周期地震仪接收函数波形反演台站下方S波速度结构,获得下地壳到上地幔顶部的速度偏差较大(约0.3 km/s),可能会造成错误解释(如下地壳低速层),因此需要和其它对波速值敏感的数据(如面波频散)进行联合分析.     

利用中国地震科学台阵研究青藏高原东南缘地壳各向异性:第一期观测资料的剪切波分裂特征

中国地震科学台阵第一期(2011-2013年)布设在南北地震带南段,本研究利用中国地震科学台阵布设在云南及相邻地区的部分流动台站记录到的2011年6月至2013年3月的数字地震波形资料,开展地壳各向异性分析.本文使用剪切波分裂系统分析方法(SAM方法),获得了研究区内67个台站的剪切波分裂参数.研究结果表明,受到云南及周边地区复杂的构造、应力环境和纵横交错的断裂分布的影响,该地区快剪切波偏振方向(PAZ)整体上显示出NNE向和NE向的优势取向,但在空间分布上比较复杂,虽然大部分台站的PAZ与构造应力场方向一致,但部分断裂附近台站的PAZ受到断裂的影响.结果显示,本研究区内不同区域的PAZ有一定差异性.本研究划分了5个子区,西部3个不同区域的PAZ从北到南分别为NNW向、近N-S向和NE向,有顺时针旋转的趋势,而东部的2个区域PAZ分别为NEE向和NNW向.研究证实,青藏东南缘地区的地壳各向异性空间分布虽然非常复杂,但大体上与区域内的主压应力的方向和断裂分布相关.     

Recent progress in imaging polarimetry

Recent instrumental developments in imaging polarimetry allow array detectors to reach a polarimetric sensitivity of 1 × 10^–4 of the intensity. New instrumental effects appear at these levels of sensitivity and generate spurious polarization signals with amplitudes of up to 5 × 10^–4. Here I discuss these effects and present methods to avoid them. Polarized spectra with an rms noise of 6 × 10^–6 may then be obtained. Furthermore a method is brought to the reader's attention that allows polarization measurements at the 1 × 10^–4 level with regular array detectors, e.g. in the near-infrared.     

基于贝叶斯理论的接收函数与环境噪声联合反演

基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27.     

轴对称地层中高分辨率阵列侧向测井信赖域反演法

本文研究轴对称地层中高分辨率阵列侧向测井（HRLA）的多参数信赖域反演方法.首先对前期HRLA的有限元正演方法进行改进,提出基于叠加原理和并行直接稀疏求解器PARDISO的快速正演方案,更适合于反演计算.将HRLA反演问题转化为非线性最小二乘问题,利用信赖域算法求解.为提高反演速度,推导了目标函数对优化参数偏导数的具体计算公式.对典型地层模型,与已有文献中Jacobi预条件共轭梯度法（JCG）计算结果比较,发现PARDIDO比JCG快10倍以上,验证了本文正演程序的正确性与高效性.利用信赖域算法求解了电阻率四参数反演和传统的三参数反演.研究结果表明：并行直接稀疏求解器PARDISO能有效求解此类HRLA正演问题,对6次不同探测深度的测井模拟,所形成的有限元刚度矩阵完全相同,只须进行一次矩阵分解,大大加快了正反演的速度.信赖域算法收敛速度快,且具有全局收敛性. HRLA的信赖域反演结果几乎不依赖于初值的选取,从较差初值出发仍能得到满意的反演结果.另外信赖域算法抗噪能力比较强,即使对测井数据添加信噪比为10dB（甚至5dB）的高斯白噪声,仍能通过反演得到较为准确的地层参数.     

Site characterization by seismic noise in Istanbul,Turkey

Single station seismic noise measurements were carried out at 192 sites in the western part of Istanbul, Turkey. This extensive survey allowed the fundamental resonance frequency of the sedimentary cover to be mapped, and identify areas prone to site amplification. The results are in good agreement with the geological distribution of sedimentary units, indicating a progressive decrease of the fundamental resonance frequencies from the northeastern part, where the bedrock outcrops, towards the southwestern side, where a thickness of some hundreds meters for the sedimentary cover is estimated. The particular distribution of fundamental resonance frequencies indicates that local amplification of the ground motion might play a significative role in explaining the anomalous damage distribution after the 17 August 1999 Kocaeli Earthquake.     

Optimisation of arrays of flap-type oscillating wave surge converters

In this paper a finite array of hinged flap-type wave energy converters are modelled using a mathematical approach. These are illustrative of the Oyster device of Aquamarine Power Ltd.¹ A novel semi-analytic solution method is presented for a set of boundary-value problems involving the scattering and radiation of waves by thin barriers used to model the device hydrodynamics. The approach makes use of the geometry to apply Fourier transforms, deriving non-singular integral equations in terms of the jumps in pressure across the flaps. These are then solved numerically using a highly efficient Galerkin expansion method. The focus of the results is on optimisation. We suggest optimal parameters for a single device, identifying flap length as crucial to device performance. This optimisation is then carried through to arrays with optimal arrangements and spacings being determined for a model sea state. Here, the lateral displacement of the devices emerges as a critical factor in optimal array configuration.     

斜井泥浆侵入仿真及其阵列侧向测井响应数值模拟

为了研究斜井泥浆侵入储层的井周岩石物理特征,采用广义曲线坐标系统模拟斜井泥浆侵入特性,进而研究泥浆侵入过程的阵列侧向测井响应.数值仿真表明,采用广义曲线坐标系消除了笛卡儿坐标系统在斜井边界处网格划分的锯齿现象;储层非均质造成泥浆侵入不均匀,储层渗透性越好,泥浆滤液侵入越深;泥浆滤液侵入使得井周地层饱和度、孔隙水矿化度的分布发生变化,进而造成斜井井周岩石电阻率分布剖面复杂变化;泥浆侵入过程中阵列侧向视电阻率变化以及深浅视电阻率的幅度差异,有效反映了斜井泥浆侵入特性;文中采用了基于Marquardt方法的阵列侧向测井四参数反演,有效恢复了斜井储层侵入剖面,得到了原状地层电阻率.斜井泥浆侵入特征及其阵列侧向测井响应模拟分析,有助于正确认识斜井井周岩石物理特征和利用阵列侧向测井进行斜井储层评价.     

Software Development For The Square Kilometre Array

Software development costs for the Square Kilometre Array are likely to be very large – in the range of 1000–2000 person-year a total. This level of software effort is unprecedented in radio astronomy. Consequently the risk associated with software development is very large. This is common to many large science projects and so we can learn from such projects how to best mitigate against the risk. We present a shopping list of suggestions drawn from the experience in other projects.     

利用阵列感应测井进行储层渗透率评价

钻井过程中储层受到泥浆侵入影响的程度与储层岩性有着密切关系,其中储层渗透率对侵入深度有着较大影响,因此若可以获知泥浆侵入深度,则有望对储层渗透率进行评估.本文首先建立含泥饼增长的泥浆侵入数值模型,然后建立阵列感应测井数值模型,两者的联合正演模拟显示泥浆侵入对地层的影响可以反映在阵列感应测井响应上,利用阻尼最小二乘法对阵列感应测井响应进行反演可以得到侵入深度.对侵入深度和储层渗透率的关系进行分析发现：在渗透率为1~100 mD（1 mD=0.987×10-3μm2）数量级的储层中,渗透率的变化可以在侵入深度上得到反映.以储层和井数据进行二维数值模拟发现：利用阵列感应测井响应反演出来的侵入深度曲线反映了渗透率在地层上的变化趋势,采用解释图版的方法可以对储层各层段的渗透率进行粗略估算.