国家数字地震台数据处理研究

为了突显被淹没在噪声背景下的有用信息,对国家数字地震台的数据处理方法进行了研究,并用Matlab程序给出了详细的设计过程。     

航空重力数据谐波小波滤波中卷积运算的处理

数据的滤波提取是航空重力测量数据处理的一个重要环节。谐波小波滤波器在频域具有良好的盒形谱特征,可以较好地实现信号的滤波。但谐波小波滤波器是一个连续函数,航空重力测量采集的信号是离散的,无法直接进行卷积运算。讨论了谐波小波用于航空重力数据滤波处理时的卷积运算处理方法,给出了谐波小波频域卷积的算法实现过程和仿真算例。     

Automated threshold selection methods for extreme wave analysis

The study of the extreme values of a variable such as wave height is very important in flood risk assessment and coastal design. Often values above a sufficiently large threshold can be modelled using the Generalized Pareto Distribution, the parameters of which are estimated using maximum likelihood. There are several popular empirical techniques for choosing a suitable threshold, but these require the subjective interpretation of plots by the user.In this paper we present a pragmatic automated, simple and computationally inexpensive threshold selection method based on the distribution of the difference of parameter estimates when the threshold is changed, and apply it to a published rainfall and a new wave height data set. We assess the effect of the uncertainty associated with our threshold selection technique on return level estimation by using the bootstrap procedure. We illustrate the effectiveness of our methodology by a simulation study and compare it with the approach used in the JOINSEA software. In addition, we present an extension that allows the threshold selected to depend on the value of a covariate such as the cosine of wave direction.     

改进的DisTrad模型在地形起伏区地表温度空间分辨率提升的应用

DisTrad（Disaggregation procedure for radiometric surface temperature）模型是常用于遥感地表温度空间分辨率提升的主要模型之一。DisTrad模型常面向空间范围有限、地形相对平坦的研究区域,且常选用植被参数（如植被指数或植被覆盖度等）作为关键参数。然而在空间范围较大、地形起伏地区,地表温度的空间变异可能无法完全通过植被参数解释。本研究选取四川盆地及毗邻地区为研究区,通过模拟数据研究DisTrad模型在地形起伏区地表温度空间分辨率提升中的适用性。数字高程模型（Digital Elevation Model,DEM）、归一化差值植被指数（Normalized Difference Vegetation Index,NDVI）等参数,采用滑动窗口逐步回归,将空间分辨率为6km的地表温度提升至空间分辨率为1km。研究结果表明,改进的模型在平原及海拔较低的高原地区提升获得的地表温度空间分辨率具有较高精度,均方根误差（Root Mean Square Error,RMSE）为0.5K左右;在地形起伏较大的地区,RMSE为4K,验证了改进的模型提升的可行性。     

Cavity expansion analyses of crushable granular materials with state‐dependent dilatancy

Crushability is one of the important behaviors of granular materials particularly under high stress states, and affects both the deformability and strength of the materials that are in essence associated with state‐dependent dilatancy. In this presentation, first, a new critical state model is proposed to take into account the three different modes of compressive deformation of crushable granular materials, i.e. particle rearrangement, particle crushing and pseudo‐elastic deformation. Second, the governing equations for cavity expansion in crushable granulates are introduced, in which the state‐dependent dilatancy as well as the bounding surface plasticity model are used. Then, the procedure to obtain semi‐analytical solutions to cavity expansion in the material is described in detail, in which a commercial differential equation solver is employed. Finally, cavity expansion analyses are carried out on Toyoura sand, a well‐documented granular material, to demonstrate the effects of crushability and state‐dependent dilatancy. The study shows that particle crushing does occur at both high stress and critical states and affects the stress fields and the deformation behavior of the material surrounding the cavity in association with state‐dependent dilatancy. This leads to conclusion that particle crushing and state‐dependent dilatancy have to be taken into account when cavity expansion theory is used to interpret cone penetration tests and pressuremeter tests. Copyright © 2011 John Wiley & Sons, Ltd.     

Stress‐dependent hardening and failure surfaces of dry sand

During several triaxial compression experiments on plastic hardening, softening, and failure properties of dense sand specimens, it was found on various stress paths that the size of the failure surface was not constant. Instead, it changed depending on the current state of hydrostatic pressure. This finding is in contrast to the standard opinion consisting of the fact that the failure surface remains constant, once it has been reached during an experiment or in situ. In general, the behaviour of cohesionless granular‐material‐like sand is somehow characterised in between fluid and solid, where the solid behaviour results from the angle of internal friction and the confining pressure. Although the friction angle is an intrinsic material property, the confining pressure varies with the boundary conditions, thus defining different solid properties like plastic hardening, softening, and also failure. Based on our findings, it was the goal of the present contribution to introduce an improved setting for the plastic strain hardening and softening behaviour including the newly found yield properties at the limit state. For the identification of the material parameters, a complete triaxial experimental analysis of the tested sand is given. The overall elasto‐plasticity concept is validated by numerical computations of several laboratory foundation‐ and slope‐failure experiments. The performance of the proposed approach is compared with the standard concept of a constant failure surface, where the corresponding yield surfaces are understood as contours of equivalent plastic work or plastic strain. Copyright © 2012 John Wiley & Sons, Ltd.     

漳州-华安“地震窗”及其机制

漳州-华安是一个中强地震活动区,现今小震活动频繁。本文分析该区历史和现今地震活动特点,发现该区地震活动性涨落与闽台强震活动存在一定的相关性,因而得出该区是一个“地震窗”的认识。进而对“地震窗”的形成条件进行讨论,认为该区“地震窗”的形成是与其深部热作用所产生的向上垂直作用力紧密相关     

Archaeomagnetic secular variation in the UK during the past 4000 years and its application to archaeomagnetic dating

This paper examines the limitations and deficiencies of the current British archaeomagnetic calibration curve and applies several mathematical approaches in an attempt to produce an improved secular variation curve for the UK for use in archaeomagnetic dating. The dataset compiled is the most complete available in the UK, incorporating published results, PhD theses and unpublished laboratory reports. It comprises 620 archaeomagnetic (directional) data and 238 direct observations of the geomagnetic field, and includes all relevant information available about the site, the archaeomagnetic direction and the archaeological age. A thorough examination of the data was performed to assess their quality and reliability. Various techniques were employed in order to use the data to construct a secular variation (SV) record: moving window with averaging and median, as well as Bayesian statistical modelling. The SV reference curve obtained for the past 4000 years is very similar to that from France, most differences occurring during the early medieval period (or Dark Ages). Two examples of dating of archaeological structures, medieval and pre-Roman, are presented based on the new SV curve for the UK and the implications for archaeomagnetic dating are discussed.     

A Least-squares Window Curves Method for Interpretation of Magnetic Anomalies Caused by Dipping Dikes

We have developed a least-squares method to determine simultaneously the depth and the width of a buried thick dipping dike from residualized magnetic data using filters of successive window lengths. The method involves using a relationship between the depth and the half-width of the source and a combination of windowed observations. The relationship represents a family of curves (window curves). For a fixed window length, the depth is determined for each half-width value by solving one nonlinear equation of the form f (z) = 0 using the least-squares method. The computed depths are plotted against the width values representing a continuous curve. The solution for the depth and the width of the buried dike is read at the common intersection of the window curves. The method involves using a dike model convolved with the same moving average filter as applied to the observed data. As a result, this method can be applied to residuals as well as to measured magnetic data. Procedures are also formulated to estimate the amplitude coefficient and the index parameter. The method is applied to theoretical data with and without random errors. The validity of the method is tested on airborne magnetic data from Canada and on a vertical component magnetic anomaly from Turkey. In all cases examined, the model parameters obtained are in good agreement with the actual ones and with those given in the published literature.