用振幅变化长短时均值比实现地震预警中P波自动拾取

地震P波到时快速、精确的自动拾取是实现地震预警的第一步。目前所有P波到时自动拾取方法都需要进一步考虑准确度问题,据此提出用"振幅变化"长短时均值比方法拾取P波到时。研究表明引进"振幅变化"作为特征函数比Allen的"振幅变化平方"更能突显地震波初至时刻的"突变"特征。对"振幅变化长短时均值比方法"和Allen的"振幅变化平方长短时均值比方法"进行抗噪音分析表明,在较低信噪比条件下,Allen方法不能准确地拾取P波到时,而新方法依然能较准确地拾取P波到时。在较高信噪比条件下,用两种方法拾取的160个M_S4.0地震P波到时,其误差大都在允许范围内,但新方法拾取的准确度略高。新方法可以看做是对长短时均值比方法的拓展。     

CHARACTERISTICS OF TECTONIC STRESS FIELD OF THE XIAOWAN RESERVOIR BEFORE AND AFTER THE IMPOUNDMENT

Using the seismic waveform data of Xiaowan seismic network and Yunnan seismic network, we determined the focal mechanisms of 36 earthquakes(M_L ≥ 3.0)from Jun. 2005 to Dec. 2008 and 51 earthquakes(M_L ≥ 2.5)from Jan. 2009 to Dec. 2015 by generalized polarity and amplitude technique. We inverted tectonic stress field of the Xiaowan reservoir before impounding, using the focal mechanisms of 36 earthquakes(M_L ≥ 3.0)from Jun. 2005 to Dec. 2008 and CAP solutions of 58 earthquakes(M_L ≥ 4.0)collected and the solutions in the Global Centroid Moment Tensor(GCMT)catalog; We inverted local stress field of the reservoir-triggered earthquake clustering area, using 51 earthquakes(M_L ≥ 2.5)from Jan. 2009 to Dec. 2015. Focal mechanisms statistics show that, the Weixi-Qiaohou Fault is the seismic fault. Focal mechanisms were strike-slip type in initial stage, but normal fault type in later stage. Focal depths statistics of 51 earthquakes(M_L ≥ 2.5)show that, the average value of focal depths in period Ⅰ, period Ⅱ and period Ⅲ are 8.2km, 7.3km and 7.8km respectively and the standard deviations are 4.3km, 3.5km and 6.0km respectively. The average value of focal depths is basically stable in different period, only the standard deviation is slightly different. Therefore, there is not positive connection between focal depth and deviation of focal mechanisms. What's more, there are 2 earthquakes(number 46 and number 47 in Fig.5 and Table 3)with almost the same magnitude, epicenter and focal depth, but they have different faulting types as normal and strike-slip. The focal mechanism of event No.46 is strike:302°, dip:40° and rake:-97° for plane Ⅰ, however, the focal mechanism of event No.47 is strike:292°, dip:82° and rake:140° for plane Ⅰ. Likewise, earthquake of number 3 and number 18 have similar characteristic. Therefore, the obvious focal mechanism difference of similar earthquake pair indicates the complexity of Weixi-Qiaohou Fault. Considering the quiet-active character of reservoir-triggered earthquakes, we discussed the change of local stress field in different period. The σ₁ of tectonic stress field was in the near-south direction, with a dip angle of 14° before the impoundment, however, the direction of σ₁ of local stress field changed continuously, with the dip angle getting larger after the impoundment. The direction of σ₁ of local stress field of reservoir-triggered earthquake clustering area is close to the strike of Weixi-Qiaohou Fault, and reservoir impoundment increased the shear stress in the fault, so the weakening of fault was beneficial to trigger earthquakes. Comprehensive analysis suggests that fluid permeation and pore pressure diffusion caused by the water impounding, and the weakening of fault caused by local stress field are the key factors to trigger earthquake in the Xiaowan reservoir.     

Predictive optimal control for seismic analysis of non-linear and hysteretic structures

In this paper, an effective active predictive control algorithm is developed for the vibration control of non-linear hysteretic structural systems subjected to earthquake excitation. The non-linear characteristics of the structural behaviour and the effects of time delay in both the measurements and control action are included throughout the entire analysis (design and validation). This is very important since, in current design practice, structures are assumed to behave non-linearly, and time delays induced by sensors and actuator devices are not avoidable. The proposed algorithm focuses on the instantaneous optimal control approach for the development of a control methodology where the non-linearities are brought into the analysis through a non-linear state vector and a non-linear open-loop term. An autoregressive (AR) model is used to predict the earthquake excitation to be considered in the prediction of the structural response. A performance index that is quadratic in the control force and in the predicted non-linear states, with two additional energy related terms, and that is subjected to a non-linear constraint equation, is minimized at every time step. The effectiveness of the proposed closed-open loop non-linear instantaneous optimal prediction control (CONIOPC) strategy is presented by the results of numerical simulations. Since non-linearity and time-delay effects are incorporated in the mathematical model throughout the derivation of the control methodology, good performance and stability of the controlled structural system are guaranteed. Copyright © 1999 John Wiley & Sons, Ltd.     

季风指数及其年际变化Ⅱ.南海夏季风分量动能强度指数及其年际变化

计算了各年南海夏季风建立前后流场的场相似度、场比幅、季风分量动能强度指数和突变度.指出按变差度最大或相似度绝对值最小及其变化最陡以及比幅最小,可客观定量地定出季风来临的预兆日期,在大多年份该日期比用天气气候学方法得到的季风来临(爆发)日期要早些,且两者有较好正相关.绝大多数年份季风建立时有环流突变发生,但也有少数年份呈调和变化或二次突变.季风分量动能强度指数能够反映各年南海夏季风建立后的强度.最后分析指出,南海850hPa夏季风的前兆日期,突变度和强度指数都有明显的年际和年代际变化.     

A semi‐distributed parallel‐type linear reservoir rainfall‐runoff model and its application in Taiwan

The main purpose of this paper is to introduce a semi‐distributed parallel surface rainfall‐runoff conceptual model. In this paper, a general solution of the instantaneous unit hydrograph (IUH) has been derived successfully for N linearly connected reservoirs, each having a different storage constant. The solution is a function of geomorphologic parameters, meteorologic factors and roughness coefficients. The model also takes into account the hydrologic response which is influenced by outflow downstream of a reservoir. For calibration, the shuffled complex evolution (SCE) algorithm is used to search for the global optimal parameters of the model. Because of the parallel structure, the mean roughness parameter of the channel becomes a “conceptual” parameter without a real physical meaning. To evaluate the adaptability of the model adopted, three watersheds around the city of Taipei in Taiwan were chosen to test the effectiveness of the model. The study provides an appropriate rainfall‐runoff model for planning flood mitigation in Taiwan. Copyright © 1999 John Wiley & Sons, Ltd.     

引入非线性瞬时弹性模量的软土流变模型

软土是一种非线性流变物质,但现行流变模型却忽略了加载瞬时的非线性变形;软土流变试验参数是数学模型参数而非物理模型参数,一个土样一般只对应一种材料参数,但现行流变模型的模型参数常随着荷载变化,且变化无规律,离散性大。基于以有机质为流变相物质的人工土试样室内一维压缩蠕变试验,建立了各级荷载作用下具有相同模型参数的7组件流变模型,用非线性弹簧H描述软土加载瞬时弹塑性变形、并联的三元件组件模拟软土随时间变化的黏弹塑变形。结果表明：引入三次非线性瞬时弹性模量的7组件流变模型能够较准确地反映软土的流变特性。     

Seafloor and buried mounds on the western slope of the Niger Delta

Seafloor mounds are potential geohazards to offshore rig emplacement and drilling operations and may contain evidence of underlying petroleum systems. Therefore, identifying and mapping them is crucial in de-risking exploration and production activities in offshore domains.A 738 km² high resolution three-dimensional seismic dataset was used to investigate the occurrence, seismic characteristics and distribution of features interpreted as seafloor and buried sediment mounds, at water-depths of 800–1600 m, on the western Niger Delta slope. Fifteen seafloor mounds and eighteen shallowly buried mounds were identified. The seafloor mounds are characterised by lower seismic amplitude anomalies than the surrounding seabed sediments, and overlie vertical zones of acoustic blanking. The buried mounds in contrast are characterised by high amplitude anomalies; they also directly overlie sub-vertical zones of acoustic blanking. Seismic evidences from the features, their distribution patterns and tectono-stratigraphic associations suggest that their formation was controlled by the juxtaposition of buried channels and structural highs and their formation caused by focused fluid flow and expulsion of entrained sediments at the seabed.Considering the acoustic and geometrical characteristics of the mounds and comparing them with mound-shaped features from around the world, we conclude that the mounds most likely comprise heterolithic seafloor extrusions of muds and sands from the Agbada Formation with gas and possibly oil in some of the pore space giving rise to the acoustic characteristics.     

一种地方与区域地震震源机制反演技术:广义极性振幅技术(二)——对实际震例的应用

为了检验广义极性振幅技术(GPAT)的实用性,我们利用GPAT反演确定了49次实际地震的震源机制、矩震级和震源深度.为了检验GPAT对地方地震、区域地震以及远震的实用性,我们选用了震级范围约为ML0.2~MS7.0之间的地震与震中距范围约在5~8000km之间的观测资料.对反演结果的分析表明,利用GPAT获取的震源机制结果在合理误差范围内是正确的,利用GPAT获取的矩震级结果是可靠的,而利用GPAT获取的震源深度总体上似乎比常规定位深度深约0.6km.总体而言,GPAT在震源机制、矩震级和震源深度的反演方面表现出良好的实用性.     

数据空间磁异常模量三维反演

强剩磁的存在通常导致了总磁化强度方向未知,进而影响了磁异常的反演和解释.磁异常模量是一种受磁化方向影响小的转换量,可以在强剩磁条件下通过反演三维磁化强度大小分布来推测场源分布状态.我们提出了一种数据空间磁异常模量反演算法来减少剩磁的影响.与标准的模型空间L2范数正则化反演方法相比,我们的方法有两个优点:一是无需搜索正则化参数(需要反复求解非线性反演问题),因而可以减少计算时间;二是反演结果更加聚焦,深度分辨率更高,我们对此进行了原因分析.通过模型和实测数据测试证明了该算法的有效性和更好的反演效果.     

基于强震和测震数据对鲁甸6.5级地震发震断层方向的研究

2014年8月3日16时30分,云南省昭通市鲁甸县发生M6.5地震,宏观震中位于鲁甸县龙头山镇(27.1°N,103.3°E),深度约12km。发震断层为NW向包谷垴-小河断裂,是NE向昭通-鲁甸断裂系的NW向次级走滑断层。文中利用鲁甸地震50个强震台站观测记录与135个测震台站观测记录,分别进行峰值加速度计算和S波辐射图案匹配,对此次地震发震断层方向进行讨论。结果表明,"NW向"的断层方向更符合强震、测震仪器观测值,同时与地震灾害现场调查结果一致。