基于BSWT-DDTFA方法的地球天然脉冲电磁场震前信号时频分析研究

地球天然脉冲电磁场（ENPEMF）信号,可理解为地球天然变化磁场的瞬间扰动,携带了大量有用的地质构造及其动力学信息.研究ENPEMF信号所蕴含的时间-频率联合分布特点,有利于深入了解目标对象的地球物理现象及其地质动力学原理.本文针对ENPEMF信号的非平稳特点,在数据驱动时频分析方法（DDTFA）的基础上提出了基于二值化同步压缩小波变换的改进算法（BSWT-DDTFA）.该算法可以实现数据驱动初始相位自动赋值的功能,具有自适应性.实验仿真和实际数据均证明了该改进算法不仅能够得到较为精确的频率曲线和更加清晰的时频分布,而且具有较强的抗噪声能力.以2013年芦山M_S7.0地震为例,利用BSWT-DDTFA方法提取ENPEMF信号的时频特性,结果表明ENPEMF信号的时间-频率-幅度分布在震前有明显的异常特征.     

Investigation of ionospheric TEC changes related to the 2008 Wenchuan earthquake based on statistic analysis and signal detection

Ionospheric TEC (total electron content) time series are derived from GPS measurements at 13 stations around the epicenter of the 2008 Wenchuan earthquake. Defining anomaly bounds for a sliding window by quartile and 2-standard deviation of TEC values, this paper analyzed the characteristics of ionospheric changes before and after the destructive event. The Neyman-Pearson signal detection method is employed to compute the probabilities of TEC abnormalities. Result shows that one week before the Wenchuan earthquake, ionospheric TEC over the epicenter and its vicinities displays obvious abnormal disturbances, most of which are positive anomalies. The largest TEC abnormal changes appeared on May 9, three days prior to the seismic event. Signal detection shows that the largest possibility of TEC abnormity on May 9 is 50.74%, indicating that ionospheric abnormities three days before the main shock are likely related to the preparation process of the M_S8.0 Wenchuan earthquake.     

Influence of the soil–structure interaction on the fundamental period of buildings

This paper includes an investigation of the influence of the soil–structure interaction (SSI) on the fundamental period of buildings. The behaviour of both the soil and the structure is assumed to be elastic. The soil‐foundation system is modelled using translational and rotational discrete springs. Analysis is first conducted for one‐storey buildings. It shows that the influence of the SSI on the fundamental frequency of building depends on the soil–structure relative rigidity K_ss. Analysis is then extended for multi‐storey buildings. It allows the generalization of the soil–structure relative rigidity K_s to such complex structures. Charts are proposed for taking into account the influence of the SSI in the calculation of the fundamental frequency of a wide range of buildings. Copyright © 2007 John Wiley & Sons, Ltd.     

Ground motion directionality in the 2010–2011 Canterbury earthquakes

This paper examines the observed directionality of ground motions in the Christchurch urban area during the 2010–2011 Canterbury, New Zealand earthquakes. A dataset of ground motions recorded at 20 strong motion stations over 10 different earthquake events is utilized to examine the ratios of various response spectral directionality definitions and the orientation of the maximum direction. Because the majority of previous related studies have utilized overlapping ground motion datasets from the NGA database, the results of this study provide a largely independent assessment of these ground motion aspects. It is found that the directionality ratio between the maximum (100th percentile) and 50th percentile orientation‐independent spectral acceleration is similar to that obtained from recent studies. Ground motions from the 4 September 2010 Darfield earthquake are shown to exhibit strong directionality for source‐to‐site distances up to R_rup = 30 km, notably further than results from a previous study, which suggests that such effects are generally limited to R_rup < 5 km. The adopted dataset also offers the unique potential to consider site‐specific effects on directionality ratios and maximum direction orientations; however, in both cases, site‐specific effects are found not to be significant in the observed empirical results. Copyright © 2014 John Wiley & Sons, Ltd.     

青海玛多MS7.4地震前LURR与OLR短临异常关联性研究

加卸载响应比(LURR)是根据岩石应力-应变关系的动态演化规律提出的地震预测方法。以Benioff应变为响应量,分析青海玛多M_S7.4地震前震中及其周边400 km内的LURR异常的演化特征发现,在震前一个月LURR值达到最高值,半个月之后出现高值回落的现象,说明此时孕震区岩石介质已进入屈服阶段的末端。在此基础上,利用NOAA卫星地面长波辐射(OLR)数据,分析高值回落后(5月14—25日)地面长波辐射值的短临异常演化规律。结果显示:震前在空间上仅仅在震中北部出现OLR显著且连续的短临异常变化并呈现向震中扩展的趋势,在时间上整个演化过程经历了初始增温-异常扩展-达到峰值-异常减弱-地震发生-异常消失等6个阶段。LURR和OLR震前短临异常特征说明,回落前的LURR高值反映了岩石介质屈服阶段的结束,OLR在一定程度上可以表征处于临界状态的活动断层的构造应力累积状态,综合研究两种物理参量震前异常演化特征,有助于更准确评估地震活动趋势。     

Bridge–embankment interaction under transverse ground excitation

Seismic performance and dynamic response of bridge–embankments during strong or moderate ground excitations are investigated through finite element (FE) modelling and detailed dynamic analysis. Previous research studies have established that bridge–embankments exhibit increasingly flexible performance under high‐shear deformation levels and that soil displacements at bridge abutment supports may be significant particularly in the transverse direction. The 2D equation of motion is solved for the embankment, in order to evaluate the dynamic characteristics and to describe explicitly the seismic performance and dynamic response under transverse excitations accounting for soil nonlinearities, soil–structure interaction and imposed boundary conditions (BCs). Using the proposed model, equivalent elastic analysis was performed so as to evaluate the dynamic response of approach embankments while accounting for soil–structure interaction. The analytical procedures were applied in the case of a well‐documented bridge with monolithic supports (Painter Street Overcrossing, PSO) which had been instrumented and embankment participation was identified from its response records after the 1971 San Fernando earthquake. The dynamic characteristics and dynamic response of the PSO embankments were evaluated for alternative BCs accounting for soil–structure interaction. Explicit expressions for the evaluation of the critical embankment length L_c are provided in order to quantify soil contribution to the overall bridge system under strong intensity ground excitations. The dynamic response of the entire bridge system (deck–abutments–embankments) was also evaluated through simplified models that considered soil–structure interaction. Results obtained from this analysis are correlated with those of detailed 3D FE models and field data with good agreement. Copyright © 2007 John Wiley & Sons, Ltd.     

汶川地震频次场异常特征分析

地震频次场是描述地震发生频次时空特征的一种数学方法.将2008年5月12日汶川8.0级地震震中附近区域(30.0°～33.0°N、101.5°～105.5°E)作为研究对象,以自然正交函数展开方法分析频次场典型场时间因子的时间变化特征.当取前8个特征值对应的典型场时,拟合精度可达0.9368;其中6个典型场显示有异常变...     

Simple generalization approach for intensity–duration–frequency relationships

Rainfall intensity–duration–frequency (IDF) relationships describe rainfall intensity as a function of duration and return period, and they are significant for water resources planning, as well as for the design of hydraulic constructions. In this study, the two‐parameter lognormal (LN2) and Gumbel distributions are used as parent distribution functions. Derivation of the IDF relationship by this approach is quite simple, because it only requires an appropriate function of the mean of annual maximum rainfall intensity as a function of rainfall duration. It is shown that the monotonic temporal trend in the mean rainfall intensity can successfully be described by this parametric function which comprises a combination of the parameters of the quantile function a(T) and completely the duration function b(d) of the separable IDF relationship. In the case study of Aegean Region (Turkey), the IDF relationships derived through this simple generalization procedure (SGP) may produce IDF relationships as successfully as does the well‐known robust estimation procedure (REP), which is based on minimization of the nonparametric Kruskal–Wallis test statistic with respect to the parameters θ and η of the duration function. Because the approach proposed herein is based on lower‐order sample statistics, risks and uncertainties arising from sampling errors in higher‐order sample statistics were significantly reduced. The authors recommend to establish the separable IDF relationships by the SGP for a statistically favorable two‐parameter parent distribution, because it uses the same assumptions as the REP does, it maintains the observed temporal trend in the mean additionally, it is easy to handle analytically and requires considerably less computational effort. Copyright © 2012 John Wiley & Sons, Ltd.     

The design of visco‐acoustic weighted L1‐error frequency–space wavefield extrapolators

Seismic imaging is an important step for imaging the subsurface structures of the Earth. One of the attractive domains for seismic imaging is explicit frequency–space (f – x) prestack depth migration. So far, this domain focused on migrating seismic data in acoustic media, but very little work assumed visco‐acoustic media. In reality, seismic exploration data amplitudes suffer from attenuation. To tackle the problem of attenuation, new operators are required, which compensates for it. We propose the weighted L ₁ ‐error minimisation technique to design visco‐acoustic f – x wavefield extrapolators. The L ₁ ‐error wavenumber responses provide superior extrapolator designs as compared with the previously designed equiripple L ₄ ‐norm and L_∞‐norm extrapolation wavenumber responses. To verify the new compensating designs, prestack depth migration is performed on the challenging Marmousi model dataset. A reference migrated section is obtained using non‐compensating f – x extrapolators on an acoustic dataset. Then, both compensating and non‐compensating extrapolators are applied to a visco‐acoustic dataset, and both migrated sections are then compared. The final images show that the proposed weighted L ₁ ‐error method enhances the resolution and results in practically stable images.     

Shaking table tests of a reinforced concrete waffle–flat plate structure designed following modern codes: seismic performance and damage evaluation

Shaking table tests were conducted on a scaled reinforced concrete waffle–flat plate structure. It represented a conventional construction design under current building codes in the Mediterranean area. The test structure was subjected to a sequence of four seismic simulations of increasing magnitude. Each simulation was associated with a seismic hazard level characterized by the mean return period P_R. The test structure performed well for the simulations associated with P_R = 95, 475 and 975 years but collapsed under the maximum considered earthquake of P_R = 2475 years. Damage concentrated at column bases, where the maximum chord rotation reached 93% of the ultimate capacity, and at the transverse beams of the exterior plate‐to‐column connection that failed in torsion. It is shown that most (from 85% to 90%) of the energy input by the earthquake that contributes to damage is dissipated by the plate. The capacity curve of the tested structure estimated from the experimental base shear vs. top displacement relationship allowed us to compute the overstrength (1.4). It is close to the maximum established by European code EN 1998‐1 (1.5). Based on a detailed study of the test results, potential updates to current codes and design recommendations are suggested. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupling model of FE–BE–IE–IBE for non-linear layered soil–structure interactions

A coupling model of Finite Elements (FEs), Boundary Elements (BEs), Infinite Elements (IEs) and Infinite Boundary Elements (IBEs) is presented for analysis of soil–structure interaction (SSI). The radiation effects of the infinite layered soil are taken into account by FE–IE coupling, while the underlying bed rock half-space is discretized into BE–IBE coupling whereby the non-horizontal bed rock surface can be accounted for. Displacement compatabilities are satisfied for all types of aforementioned elements. The equivalent linear approach is employed for approximation of nonlinearity of the near field soil. This model has some advantages over the current SSI program in considering the bed rock half-space and non-vertical wave incidence from the far field. Examples of verification demonstrate the applicability and accuracy of the method when compared with the FLUSH program. Finally, the effects of the relative modulus ratio E_r/E_s of rock and soil and the incident angles of non-vertical waves on the responses of the structure and the soil are examined. Copyright © 1999 John Wiley & Sons, Ltd.     

Comparison of methods for applying the Priestley–Taylor equation at a regional scale

The scenario assumed for this study was that of a region with a complete or first‐order weather station surrounded by a network of second‐order stations, where only monthly air temperature data were available. The objective was to evaluate procedures to estimate the monthly α parameter of the Priestley–Taylor equation in the second‐order stations by adjusting and extrapolating α values determined at the first‐order station. These procedures were applied in two climatic zones of north‐east Spain with semi‐arid continental and semi‐arid Mediterranean climates, respectively. Procedure A assumed α to be constant over each zone for each month (direct extrapolation). Procedure B accounted for differences in vapour pressure deficit and available energy for evapotranspiration between the first‐ and second‐order stations. Procedure C was based on equating the Penman–Monteith (P–M) and Priestley–Taylor (P–T) equations on a monthly basis to solve for α. Methods to estimate monthly mean vapour pressure deficit, net radiation and wind speed were developed and evaluated. A total of 11 automated first‐order weather stations with a minimum period of record of 6 years (ranging from 6 to 10 years) were used for this study. Six of these stations were located in the continental zone and five in the Mediterranean zone. One station in each zone was assumed to be first‐order whereas the remainder were taken as second‐order stations. Monthly α parameters were calibrated using P–M reference crop evapotranspiration (ET₀) values, calculated hourly and integrated for monthly periods, which were taken as ‘true’ values of ET₀. For the extrapolation of monthly α parameters, procedure A was found to perform slightly better than procedure B in the Mediterranean zone. The opposite was true in the continental zone. Procedure C had the worst performance owing to the non‐linearity of the P–M equation and errors in the estimation of monthly available energy, vapour pressure deficit and wind speed. Procedures A and B are simpler and performed better. Overall, monthly P–T ET₀ estimates using extrapolated α parameters and R_n?G values were in a reasonable agreement with P–M ET₀ calculated on an hourly basis and integrated for monthly periods. The methods presented for the spatial extrapolation of monthly available energy, vapour pressure deficit and wind speed from first‐ to second‐order stations could be useful for other applications. Copyright © 2001 John Wiley & Sons, Ltd.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO₂ concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO₂ on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.     

2013年岷县—漳县6.6级地震前地磁日变化异常及机理分析

利用空间相关法对2013年岷县—漳县6.6级地震前后151个台站地磁垂直分量日变化进行分析,结果显示:(1)2013年7月3日、4日在震中附近分别出现16和17个台站的地磁Z分量日变化明显反相位异常,具有明显的分界线.(2)利用S变换方法对分界线两侧的天水和兰州台的Z分量日变化进行分频率相关系数计算,发现7月3日、4日...     

Geomagnetic diurnal-variation anomalies and their relation to strong earthquakes

The diurnal-variation anomalies of the vertical-component in geomagnetic field are mainly the changes of phase and amplitude before strong earthquakes. On the basis of data recorded by the network of geomagnetic observatories in China for many years, the anomalous features of appearance time of the minima of diurnal variations (i.e, low-point time) of the geomagnetic vertical components and the variation of their spatial distribution (i.e, phenomena of low-point displacement) have been studied before over 30 strong earthquakes with M _S≥6.6 such as Kunlunshan M _S=8.1 earthquake on November 14, 2001; Bachu-Jashi M _S=6.8 earthquake on February 24, 2003; Xiaojin M _S=6.6 earthquake on September 22, 1989, etc. There are good relations between such rare phenomena of geomagnetic anomalies and the occurrence of earthquakes. It has been found that most earthquakes occur in the vicinity of the boundary line of sudden change of the low-point displacement and generally within four days before and after the 27th or 41st day counting from the day of appearance of the anomaly. In addition, the anomalies of diurnal-variation amplitude near the epicentral area have been also studied before Kunlunshan M _S=8.1 earthquake and Bachu-Jiashi M _S=6.8 earthquake. Foundation item: National Science Technology Tackle Key Project during the Tenth Five-year Plan (2001BA601B01-05-04)     

2021年青海玛多MS7.4地震前地下流体异常特征分析

地震地下流体在地震预测研究与震情跟踪方面发挥着重要作用,本文对2021年5月22日青海玛多MS7.4地震前震中附近地下流体观测资料及异常特征进行系统分析,并结合震前预测过程进行了回顾性总结。结果表明,玛多MS7.4地震震中距500km范围内存在5项异常,其中4项为短期异常,主要出现在震前2～3个月。异常在空间上分布不均匀,主要位于玛多MS7.4地震震中的东北部和西南部,整体呈现出由外围向震中收缩的迁移特征,且玛多MS7.4地震和2010年青海玉树MS7.1地震在异常特征和应力加载作用方面具有较好的相似性。研究结果为地震监测能力较低地区积累了震例资料,对于提升强震地下流体前兆异常认识及未来震情趋势判定水平具有一定的参考意义。     

2017年四川九寨沟7.0级地震前地震应变场分析

以地震应变场作为地震活动的变量,通过自然正交函数展开方法,计算2017年8月8日四川九寨沟7.0级地震前的地震应变场,提取出震前时间因子的异常变化。计算结果发现应变场前4个时间因子在震前1~3年的中短期异常并不显著,只有第2个和第4个时间因子震前有小幅度的异常变化,分析其原因,时间因子可能受到2008年汶川8.0级地震和2013年岷县漳县6.6级地震前大幅度异常的影响。与时间因子对应的空间等值线形成局部应变高值异常危险区,可能是2013年岷县漳县6.6级地震和2017年九寨沟7.0级地震空间异常的中短期特征。再对比九寨沟7.0级地震前后与松潘-平武7.2级双震的空间异常随时间演变,分析异常发展变化模式的差异,最后应用地震资料的累积频次从物理角度解释异常形成的机制。     

2019年9月16日张掖5.0级地震前多尺度电磁异常现象

2019年9月16日甘肃张掖发生5.0级地震，地震发生前既有空间范围跨度较大的地磁低点位移异常，也有单台电磁扰动异常；既有震前2年开始出现的山丹台视电阻率的年变形态畸变，也有震前3个月之内出现的高台钻孔应变和高台大庄低频电磁扰动异常。本次地震震级不高却记录到丰富的电磁场的异常信号，表明该地区处于地震电磁异常的敏感区。高台钻孔应变在震前观测到同步的异常变化，揭示孕震区应力状态发生了显著的变化。     

Step–pool evolution in the Rio Cordon,northeastern Italy

The principle that formative events, punctuated by periods of evolution, recovery or temporary periods of steady‐state conditions, control the development of the step–pool morphology, has been applied to the evolution of the Rio Cordon stream bed. The Rio Cordon is a small catchment (5 km²) within the Dolomites wherein hydraulic parameters of floods and the coarse bedload are recorded. Detailed field surveys of the step–pool structures carried out before and after the September 1994 and October 1998 floods have served to illustrate the control on step–pool changes by these floods. Floods were grouped into two categories. The first includes ‘ordinary’ events which are characterized by peak discharges with a return time of one to five years (1·8–5·15 m³ s^?1) and by an hourly bedload rate not exceeding 20 m³ h^?1. The second refers to ‘exceptional’ events with a return time of 30–50 years. A flood of this latter type occurred on 14 September 1994, with a peak discharge of 10·4 m³ s^?1 and average hourly bedload rate of 324 m³ h^?1. Step–pool features were characterized primarily by a steepness parameter c = (H/L_s)/S. The evolution of the steepness parameter was measured in the field from 1992 to 1998. The results indicate that maximum resistance conditions are gradually reached at the end of a series of ordinary flood events. During this period, bed armouring dominate the sediment transport response. However, following an extraordinary flood and unlimited sediment supply conditions, the steepness factor can suddenly decrease as a result of sediment trapped in the pools and a lengthening of step spacing. The analogy of step spacing with antidune wavelength and the main destruction and transformation mechanism of the steps are also discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

河北尚义M4.0地震加速度记录时频特性分析

2016年6月23日,河北张家口尚义县发生M 4.0地震,发震构造为尚义-平泉断裂。采用基于聚类经验模态分解（EEMD）的HHT计算方法,对尚义M4.0地震强震动记录进行EEMD分解、Hilbert谱变换及谱分析,得到强震动台站的加速度记录分量的时程、FFT谱和边际谱、瞬时能量谱及时频谱,对比分析并认识各种谱的时频特性及规律,做好灾前预防工作,以减少人民生命财产损失。