Amplitude-compensated Laplacian filtering of reverse time migration and its application

Reverse time migration is an advanced seismic migration imaging method. When the source wavefield and the receiver wavefield are cross-correlated, the cross-correlations of direct arrivals, backscattered waves and overturned waves will produce a lot of low-frequency noise, which will mask the final imaging results. Laplacian filtering, as a common method to suppress low-frequency noise, can adapt to any complex media, just adding a little computational cost. However, simple direct Laplacian filtering will destroy the characteristics of the useful signals. Therefore, the amplitude needs to be compensated before filtering when using the Laplacian filtering method. Zhang and Sun proposed an improved Laplacian filtering method and gave a simple calculation formula and explanation. This method can effectively suppress the low-frequency noise in reverse time migration while retaining the useful signal characteristics, but lacks detailed and strict mathematical derivation. Therefore, this paper gives a detailed and rigorous mathematical derivation of the amplitude-compensated Laplace filtering method from the point of view of amplitude-preserved filtering. The source wavelet is used instead of the source wavefield to compensate amplitude, just adding a little calculation cost. Finally, the amplitude-compensated Laplace filtering method is verified by two theoretical models and compared with the direct Laplacian filtering method.     

长江入海口打击非法捕捞政策系统的协调性研究

长江流域禁捕为全局计、为子孙谋,但在禁捕后却存在非法捕捞的问题,尤其在渔业资源较为丰富的长江入海口处尤为严重,阻碍政策目标的实现。基于现实层面的矛盾,文章选取24份国家层面、30份地方层面有关长江入海口“打击非法捕捞”的政策文本,对政策系统的协调性进行分析与评价。结果表明：长江入海口打击非法捕捞政策系统协调性较好,但存在部分不协调问题。文章创新性提出政策系统协调性整合框架,从公共政策视角为长江入海口非法捕捞问题的治理提供优化方向。     

STUDY ON THE DISTRIBUTION PATTERN OF EARTHQUAKE- TRIGGERED LANDSLIDES BASED ON SEISMIC LANDSLIDE SUSCEPTIBILITY ANALYSIS: A CASE STUDY OF LANDSLIDES TRIGGERED BY THE MS6.5 LUDIAN EARTHQUAKE IN 2014

On August 3, 2014, an M_W6.5 earthquake occurred in Ludian County, Yunnan Province, which triggered significant landslides and caused serious ground damages and casualties. Compared with the existing events of earthquake-triggered landslides, the spatial distribution of co-seismic landslides during the Ludian earthquake showed a special pattern. The relationship between the co-seismic landslides and the epicenter or the known faults is not obvious, and the maximum landslide density doesn't appear in the area near the epicenter. Peak ground acceleration (PGA), which usually is used to judge the limit boundary of co-seismic landslide distribution, cannot explain this distribution pattern. Instead of correlating geological and topographic factors with the co-seismic landslide distribution pattern, this study focuses on analyzing the influence of seismic landslide susceptibility on the co-seismic distribution. Seismic landslide susceptibility comes from a calculation of critical acceleration values using a simplified Newmark block model analysis and represents slope stability under seismic loading. Both DEM (SRTM 90m)and geological map (1 ︰ 200^￣￣000)are used as inputs to calculate critical acceleration values. Results show that the most susceptible slopes with the smallest critical accelerations are generally concentrated along the banks of rivers. The stable slopes, which have the larger critical accelerations and are comparably stable, are in the places adjacent to the epicenter. Comparison of the distribution of slope stability and the real landslides triggered by the 2014 M_W6.1 Ludian earthquake shows a good spatial correlation, meaning seismic landslide susceptibility controls the co-seismic landslide distributions to a certain degree. Moreover, our study provides a plausible explanation on the special distribution pattern of Ludian earthquake triggered landslides. Also the paper discusses the advantages of using the seismic landslide susceptibility as a basic map, which will offer an additional tool that can be used to assist in post-disaster response activities as well as seismic landslides hazards zonation.     

2016年8月唐山震群视应力变化特征

2016年8月河北唐山发生震群活动,具有频次高、持续时间长等特点。利用首都圈数字化地震台网监测的2016年1月—2017年5月河北唐山地区M_L≥2.5地震,采用遗传算法反演震源参数,计算相应地震视应力值。分析结果如下：①此次地震序列为前震—主震—余震型;②视应力值在震群发生前明显增高,发生过程中出现高低相间现象,在M_L 4.3主震出现峰值,震源区应力状态不稳定,震后视应力恢复正常;③震群视应力与震级呈现一定正相关。     

我国和东亚地区硫化物跨边界输送态势研究

本文利用污染物三维欧拉长距离输送实用模式,较详细地模拟了我国和东亚地区硫化物跨地区、跨边界输送态势,分析了东亚地区硫化物输送通量随高度、季节变化的分布特点,讨论了不同地区边界上硫化物的跨边界输送通量及流的收支平衡,并给出了我国和周边国家和地区之间硫化物的相互输送量。结果表明,我国向外输送的硫化物占周边国家和地区总硫沉降的比例并不大,各地区硫沉降的主要来源是自身排放;硫化物长距离输送中硫酸根的输送占重要地位,各地区硫酸根沉降中的外来比例较大。我们初步总结出一个东亚地区硫化物输送的概念模式：在低层,夏季主要自西南向东北输送,冬季长江以北自北向南输送,长江以南由西南向北输送,并在长江中下游地区形成一个较强的硫化物辐合带,在高层,无论冬夏均由西向东输送,并随高空带变化,中层为二者的过渡,春秋季的情形界于冬夏之间,偏向冬季。     

太湖流域河流鱼类群落的时空分布

确定河流鱼类群落的时空分布格局及其形成机制是开展鱼类物种多样性保护与管理的科学基础.基于2013年10月和2014年5月共2次对太湖流域57个河道样点的调查数据,初步研究太湖流域河流鱼类群落结构及其多样性的季节动态和空间分布特点.共采集鱼类5051尾,计46种,其中鲤科鱼类26种,占全部物种数的57%.10月份的鱼类多样性显著高于5月份,且2个季度的鱼类群落结构存在显著性差异.5个主要水系间的鱼类多样性差异显著,总体上,沿江水系和洮滆水系鱼类多样性较低,黄浦江水系居中,而南河水系和苕溪水系较高;鱼类群落结构也随水系而显著变化,主要表现为黄浦江水系与洮滆、苕溪和沿江水系呈显著差异.在2个一级生态分区之间,鱼类多样性无显著差异但群落结构显著不同,主要因、鲫、似鳊等优势种及宽鳍鱲、尖头鱥、中华青鳉、食蚊鱼等偶见种的空间分布差异所引起;在4个二级生态分区之间,鱼类多样性和群落结构均存在显著的空间变化.     

Upper crustal P-wave velocity structure beneath two volcanic areas in northern Iran

Detailed information about the crustal structure is essential for better understanding the occurrence and mechanisms of earthquakes and volcanoes. Here we present a study of the upper crustal P-wave velocity structure of two seismically and volcanically active areas in northern Iran using the two-dimensional Pg travel time tomography method. The imaging results suggest low velocities in the upper crust beneath the Damavand and Sahand-Sabalan volcanic areas in the central and western parts of northern Iran, respectively. The upper crustal low velocities in these two areas roughly coincide with previously imaged low Pn velocity anomalies, suggesting that the Late Cenozoic volcanic activity was probably caused by the upwelling of hot materials from the mantle. The image feature of the Pg velocity structure beneath the Sahand-Sabalan volcanic area further indicates that the hot materials stored in the upper crust beneath Sahand may be larger in size than those stored beneath Sabalan. Comparison of the Pg velocity images with the earthquake distribution in north Iran suggests that earthquakes mainly occur at moderately low velocity or low to high velocity boundary areas instead of significantly low or high velocity regions. The anisotropy results show that the Pg wave fast direction is consistent with the GPS direction at high Pg velocity areas and the fast direction is inconsistent with the GPS direction but consistent with the strike direction of faults at low velocity areas. Our new upper crustal structural images provide the basic observation for better understanding of the regional seismicity and volcanism, and link the surface geological phenomena to deep crustal and mantle processes associated with the active tectonics in northern Iran.     

利用分段时频峰值滤波法抑制磁共振全波信号随机噪声

磁共振信号极其微弱,容易受到周围环境中各种电磁噪声干扰.其中随机噪声,由于频带宽、不规则、无规律、与有效信号混叠,难以抑制.近年来,采用数量级为10⁴~10⁵ Hz采样频率收录的全波磁共振信号,以其携带丰富全面的信息量,为数据处理及解释提供了更多的潜能.然而,只要随机噪声的幅度大于信号幅度,拟合得到的信号特征参数准确度就会降低.目前普遍采用的数据叠加方法仅能抑制部分随机噪声,且需要多次采集信号,探测效率低.本文针对全波磁共振信号采样点数多和信号非线性强的特点,提出采用分段时频峰值滤波（STFPF）法消噪,将全波磁共振信号分成若干段,编码为解析信号的瞬时频率,采用短窗长PWVD计算解析信号的时频分布,利用时频分布沿瞬时频率集中的特性,通过提取时频分布的峰值获得信号的无偏估计,达到抑制全波磁共振信号随机噪声的目的.为了验证消噪效果,与传统叠加法进行对比分析,仿真结果表明,对于单次采集信号,信噪比低至-5 dB时,STFPF方法依然能有效抑制信号中的随机噪声,消除随机噪声后信噪比提高23.19 dB,信号的初始振幅拟合误差为3.03%,平均横向弛豫时间拟合误差为2.7%,消噪效果优于传统叠加法,且由于无需多次采集磁共振信号,可有效提高探测效率.模型数据的反演解释进一步验证了STFPF方法的有效性,本文研究结果为实际数据处理奠定了良好的基础.     

STUDY ON THE MRF-BASED METHOD FOR DAMAGED BUILDINGS EXTRACTION FROM THE SINGLE-PHASE SEISMIC IMAGE

Earthquake events are one of the most extraordinarily serious natural calamities, which not only cause heavy casualties and economic losses, but also various secondary disasters. Such events are devastating, and have far-reaching influences. As the main disaster bearing body in earthquake, buildings are often seriously damaged, thus it can be used as an important reference for earthquake damage assessment. Identifying damaged buildings from post-earthquake images quickly and accurately is of real importance, which has guidance meaning to rescue and emergency response. At present, the assessment of earthquake damage is mainly through artificial field investigation, which is time-consuming and cannot meet the urgent requirements of rapid emergency response. Markov Random Field(MRF)combines the neighborhood system of pixels with the prior distribution model to effectively describe the dependence between spatial pixels and pixels, so as to obtain more accurate segmentation results. The support vector machine(SVM)model is a simple and clear mathematical model which has a solid theoretical basis; in addition, it also has unique advantages in solving small sample, nonlinear and high-dimensional pattern recognition problems. Thus, in this paper, a Markov random field-based method for damaged buildings extraction from the single-phase seismic image is proposed. The framework of the proposed method has three components. Firstly, Markov Random Field was used to segment the image; then, the spectral and texture features of the post-earthquake damaged building area are extracted. After that, Support Vector Machine was used to extract the damaged buildings according to the extracted features. In order to evaluate the proposed method, 5 areas in ADS40 earthquake remote sensing image were selected as experimental data, this image covers parts of Wenchuan City, Sichuan Province, where an earthquake had struck in 2008. And in order to verify the applicability of this method to different resolution images, an experimental area was selected from different resolution images obtained by the same equipment. The experimental results show that the proposed method has good performance and could effectively identify the damaged buildings after the earthquake. The average overall accuracy of the selected experimental areas is 93.02%. Compared with the result extracted by the widely used eCognition software, the proposed method is simpler in operation and can improve the extraction accuracy and running time significantly. Therefore, it has significant meaning for both emergency rescue work and accurate disaster information providing after earthquake.     

Crustal structure study based on principal component analysis of receiver functions

The receiver function(RF) technique is an effective method for studying crustal structure. For a single station, the average 1-D crustal structure is usually derived by stacking the radial RFs from all back-azimuths, whereas structural variations(such as dipping discontinuities or anisotropy) can be constrained through analysis of waveform dependence on the backazimuth of both the radial and tangential RFs. However, it is often difficult to directly extract information about structural variations from the waveform of RF, due to the common presence of noise in real data. In this study, we proposed a new method to derive structural variation information for individual stations by applying principal component analysis(PCA) to RFs sorted by back-azimuth. In this method(termed as RF-PCA), a set of principal components(PCs), which are uncorrelated with each other and reflect different characteristics of the RF data, were extracted and utilized separately to reconstruct new RFs. Synthetic tests show that the first PC of the radial RFs contains the average structural information of the crust beneath the corresponding station, and the second PC of the radial RFs and the first PC of the tangential RFs both reflect the variations of the crustal structure. Our synthetic modeling results indicate that the new RF-PCA method is valid for a variety of synthetic models with intra-crustal dipping discontinuities and/or anisotropy. We applied this method to the real data from a broadband temporary seismic station(s233) in the central part of the Sichuan Basin. The results suggest that the RF data can be best explained by the presence of two nearly parallel dipping discontinuities within the crust. Combining with previous logging data, seismic exploration and deep sounding observations, we interpret the shallow dipping discontinuity as the top boundary of the Precambrian crystalline basement of the Sichuan Basin and the deep one corresponding to the Conrad interface between the upper and lower crust, consistent with the geological feature of the study area. In this work, both synthetic tests and real data application results demonstrate the effectiveness of the RF-PCA method for studying crustal structures.