2017年九寨沟MS 7.0地震震源深度测定

为了更好地确定2017年8月8日九寨沟M_(S )7.0地震震源深度其发震机理,利用四川、甘肃和青海区域地震台网的观测波形数据,采用多种方法研究了此次地震的震源深度。首先,采用gCAP方法反演了九寨沟M_(S )7.0地震的震源机制解和矩心深度,结果显示,节面Ⅰ走向243°/倾角87°/滑动角-158°,节面Ⅱ走向151°/倾角68°/滑动角-3°,矩震级为M_(W )6.5,矩心深度为8 km;然后,采用ISOLA近震全波形方法反演了此次地震的震源机制解,反演结果与gCAP方法结果相差不大,矩心深度为7 km;最后,通过sPn震相与Pn震相之间的走时差测定此次地震初始破裂震源深度,结果显示深度约为12 km。研究表明,九寨沟M_(S )7.0地震的矩心深度为7—8 km,初始破裂深度约为12 km。     

SOURCE PARAMETERS OF THE CANGWU MS5.4 EARTHQUAKE, 31 JULY, 2016

On July 31^st, 2016, an earthquake of M_S5.4 occurred in Cangwu County, Guangxi Zhuang Autonomous Region, which is the first M_S ≥ 5.0 earthquake in coastal areas of southern China in the past 17a. The moderate earthquake activities have come into a comparatively quiet period in coastal areas of southern China for decades, so the study about the Cangwu M_S5.4 earthquake is very important. However, differernt research institutions and scholars have got different results for the focal depth of the Cangwu M_S5.4 earthquake. For this reason, we further measured the focal depth by using CAP method and sPL phase method. sPL phase was first put forward by Chong in 2010. It is often observed between P and S wave of continental earthquakes with epicentral distance of about 30km to 50km. The energy of sPL phase is mainly concentrated on the radial component. Arrival time difference between sPL phase and direct P wave is insensitive to epicentral distancs, but increases almost linearly with the increase of focal depth. Based on these characteristics and advantages, sPL phase method is chosen to measure the focal depth of Cangwu M_S5.4 earthquake in the paper. First of all, we selected the broadband waveform data through seismic stations distributed mainly in Guangxi and adjacent provinces from Data Management Centre of China National Seismic Network and Guangxi Earthquake Networks Center. And an appropriate velocity model of Cangwu area was constructed by the teleseismic receiver function method. Then, the focal mechanism and focal depth of Cangwu M_S5.4 earthquake were determined by using the CAP(Cut and Paste)method. Next, we compared the synthetic waveforms simulated by F-K forward method of different focal depth models with the actual observed waveforms. According to the difference of arrival times between sPL and Pg phases, we finally obtained the focal depth of Cangwu earthquake. The results show that the focal depth is 11km measured by CAP method and 9km by sPL phase method. Based on the focal mechanism solution, isoseismal shapes, aftershocks distributions and investigation on spot, we conclude that the Cangwu M_S5.4 earthquake is a left-lateral strike-slip earthquake which occurred in the upper crust. Our preliminary analysis considers that the seismogenic structure of Cangwu earthquake is a north-northwest branch fault, and the control fault of this earthquake is the Hejie-Xiaying Fault.     

以颗粒物理原理认识汶川M_W 7.9地震发生前临震预滑和震颤现象

发生在孕震区周边地块上的临震预滑和震颤现象,对破坏性地震预测有一定前兆意义,是值得地震学界关注的问题。选取2008年5月12日汶川M_W 7.9地震发生前,临夏和湟源地震台分量应变仪记录与临夏、恩施和西安地震台数字地震仪记录以及临夏和周至地震台深井水位仪记录,分析发现,在临震前数天至数小时,上述各地震台不同学科观测仪器均记录到一些"跃变"和"震颤"震相。文中试图以颗粒物理原理,来认识不同距离、不同台站、不同学科的观测仪器在临震前相近时间段内记录的低频和高频震相,可能是不同地块在临震前发生预滑错动后激发的预滑震相Xp和地下气体在裂隙内流动激发的震颤震相Tp。观测结果表明:2008年5月8日03时至主震发生,各地震台所处地块在相近时段内逐次发生次数不等的预滑错动,其中1—2次较大错动可在噪声背景中被识别;各地震台预滑错动方向指向或背向主震震中。据此认为:汶川M_W7.9地震前,上述各地震台所处地块在不同大小、不同方向的力链驱动下,发生指向或背向主震震中的临震预滑现象。     

Non-invasive lab test in the monitoring of vadose zone contaminated by light non-aqueous phase liquid

We discuss a strategy capable of a quantitative long-term monitoring of water saturation and volume variation of light non-aqueous phase liquids in the soil. The goal was reached monitoring a controlled sand cell contaminated with classical gasoline over 124 days, using geophysical methods such as electrical resistivity tomography, induced polarization and ground penetrating radar. We show that empirical relations, linking the water saturation to the physical parameters measured as resistivity from electrical resistivity tomography and travel time from georadar with advanced processing, are good tools for this purpose. The consistence of the proposed process is validated by both good overlap of results carried out from electrical resistivity tomography and georadar and theoretical models simulating the actual experiment.     

基于最小二乘的岩石频谱参数反演初始值选取的讨论

频谱激电法作为地球物理勘探的有效手段,在寻找激电异常体以及探寻深部矿产方面有重大贡献,岩石频谱参数的确定是实现该方法的重要步骤。在利用最小二乘法对复电阻率、相位参数进行反演时,各频谱参数初始值选取不同对反演效果影响较大,选取不好甚至会发散。在各频谱参数所属范围内采用随机选取初始值的方法进行反演,对复电阻率幅值曲线以及相位曲线拟合效果较好,不易陷入局部极小值的特点,能克服了较为盲目选取初始值而对结果带来的不良影响。     

Structure of residual oil as a function of wettability using pore-network modelling

In the water flooding of mixed-wet porous media, oil may drain down to relatively low residual oil saturations (S_or). Various studies have indicated that such low saturations can only be reached when oil layers in pore corners are included in the pore-scale modelling. These processes within a macroscopic porous medium can be modelled at the pore-scale by incorporating the fundamental physics of capillary dominated displacement within idealised pore network models. Recently, the authors have developed thermodynamic criteria for oil layer existence in pores with non-uniform wettability which takes as input geometrically and topologically representative networks, to calculate realistic S_or values for mixed-wet and oil-wet sandstones [16, 21]. This previous work is developed in this paper to include (i) the visualisation of the 3D structure of this residual oil, and (ii) a statistical analysis of this “residual/remaining” oil. Both the visualisation and the statistical analysis are done under a wide range of wettability conditions, which is reported for the first time in this paper.The structure of residual oil for strongly water wet systems is well known (where residual = remaining oil) and our model agrees with this but this structure changes radically for mixed wet systems (where residual ≠ remaining) and this has not yet been visualised experimentally. We find that for more water-wet systems high final residual oil saturations are reached at relatively small amounts of water injected and this oil is present in the pores as bulk oil. On the other hand, for more oil-wet systems we find a slow decrease of the amount of remaining oil with increasing amounts of injected water. During the process, the remaining connectivity of the oil phase is increasingly provided by oil layers only, hence the slow drainage. The final residual oil saturation, only reached in the theoretical limit of an infinite amount of injected water, is almost entirely contained in large number of (relatively low volume) oil layers, which are present in pores of most radius sizes.     

TTI介质的相速度研究

推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．     

Analysis of Short-Period Rayleigh Waves Recorded in the Bohemian Massif Area During Celebration 2000 Experiment

The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint inversion of phase and group velocity data and to provide γ_R data to be used for Q_β inversion. Validity of obtained V_S and Q_β models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian Unit. Quasi-2D model based on the individual 1D V_S models is well correlated with the surface geology. Lower V_S are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and Moldanubian Plutons, the near-surface V_S values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result of the weathering and cracks within the granitoid rocks.     

Kirchhoff型反偏移场稳相分析(英文)

在文献中,Kirchhoff型反偏移场的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像场信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像场的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移场进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。