基于多轨道卫星观测数据分析尼泊尔地震长波辐射特征

综合已经在红外异常提取中应用的涡度和RST（Robust Satellite Technique）算法优点,提出了红外异常指数算法.并基于长时间尺度的中国静止卫星FY-2D和美国极轨卫星NOAA长波辐射数据,应用RST和异常指数算法,分别对2015年4月25日尼泊尔M_S8.1和5月12日M_S7.5地震前后卫星长波辐射变化特征进行了分析,开展了多轨道、多时空分辨率长波辐射同步地震热红外特征研究.结果表明,运用RST算法,两次地震前后,未能在震中周围发现明显的长波辐射异常.运用异常指数算法：（1）对于NOAA卫星,4月15日在M_S8.1地震震中以西出现热红外异常,到4月24日震中以西约100 km处出现异常最大值,随后逐渐消失.5月10日在M_S7.5地震震中以东约200 km发现异常;（2）在NOAA卫星长波辐射异常发现最大值当日,采用FY-2D卫星每3 h的数据分析可发现红外异常的动态演化过程,弥补NOAA卫星分辨率不足.以上结果为利用多轨道卫星监测地震热辐射变化提供了依据.

     

Specific features of surface wave radiation by a shallow source

The uncertainty of the seismic moment tensor inversion from records of long period surface waves radiated by a shallow source is considered. The special cases of sources of pure thrust, normal, and strike-slip faults are examined in detail. Results of numerical modeling of Love and Rayleigh wave radiation patterns and frequency dependences of the radiation intensity are presented for these types of sources. It is shown that, in the general case, the focal mechanism and seismic moment of a shallow double couple can be uniquely determined from long period surface waves only if one of the nodal planes is subhorizontal. The threshold (maximum) value of the dip angle of this plane is determined by the source depth, the spectral range of observed surface waves, and the model structure in the vicinity of the source.     

经验方法研究下垫面变化对洪水的影响

1980s以来,人类活动对下垫面的影响加剧,研究下垫面变化对流域产汇流规律的影响具有重要意义.通过对海河流域中紫荆关、阜平两个流域分别用经验方法进行洪水模拟演算,并把洪水资料按照年代不同分时段,建立P+Pa～Rs相关关系图,分析流域下垫面变化对降雨径流相关关系的影响.建立产流量R与洪峰Qm之间的相关关系,通过统计比较不同年代的趋势线发现,在产流量相同情况下,1980s后的洪峰流量较1980s前有所减少,说明海河流域的调蓄作用有增强的趋势.同时建立流域下垫面条件改变后产流量与洪峰变化幅度的相关关系,发现流域产流量变化与洪峰变化呈正相关关系.     

An empirical method to estimate shear wave velocity of soils in the New Madrid seismic zone

In this study, a set of charts are developed to estimate shear wave velocity of soils in the New Madrid seismic zone (NMSZ), using the standard penetration test (SPT) N values and soil depths. Laboratory dynamic test results of soil samples collected from the NMSZ showed that the shear wave velocity of soils is related to the void ratio and the effective confining pressure applied to the soils. The void ratio of soils can be estimated from the SPT N values and the effective confining pressure depends on the depth of soils. Therefore, the shear wave velocity of soils can be estimated from the SPT N value and the soil depth. To make the methodology practical, two corrections should be made. One is that field SPT N values of soils must be adjusted to an unified SPT N′ value to account the effects of overburden pressure and equipment. The second is that the effect of water table to effective overburden pressure of soils must be considered. To verify the methodology, shear wave velocities of five sites in the NMSZ are estimated and compared with those obtained from field measurements. The comparison shows that our approach and the field tests are consistent with an error of less than of 15%. Thus, the method developed in this study is useful for dynamic study and practical designs in the NMSZ region.     

探地雷达探测地下目标时的波速估计

在利用探地雷达(GPR)进行工程检测过程中，电磁波在地下介质中的传播速度是影响埋深计算精度的关键参数之一。本文通过吸取几种波速求取方法的优点，给出一种精度高、操作性强的波速求取法——刚性界面反射系数法；并利用此法，由LTD探地雷达实测数据求得整条测线的速度曲线。经检验，借助于少量钻芯数据，利用此法可使公路面层厚度检测的相对误差小于5％。     

用波函数展开法求解界面圆孔的SH波散射问题

研究了ＳＨ波绕界面圆孔散射的波函数展开理论。由入射,反射和透射波组成的自由波场与孔的散射场叠加成总波场。将定义于两个半平面的散射场按一定方式延拓到全平面。通过Ｈａｎｋｅｌ－Ｆｏｕｒｉｅｒ型积分来满足界面与孔缘处的边界条件,得到了确定特定系数的封闭代数方程组。给出了具体算例,计算了孔缘处的动应力系数分布。     

内蒙古自治区东部两次Ms5．9地震前后长波辐射时序特征分析

以卫星对地遥测长波辐射(OLR)信息为基础,分析了2003年8月16日巴林左旗MS5.9地震和2004年3月24日东乌珠穆沁MS5.9地震前后震中附近的45°N,120°E点位长波辐射的时序特征。结果表明,在这两次MS5.9地震前后,位于45°N,120°E点位上长波辐射(OLR)月平均值、月距平值异常特征比较明显,这种基于不同点位上的时序特征分析对地震监视和预报具体孕震区具有一定的实际意义。     

Estimation of shear wave velocity profiles by the inversion of spatial autocorrelation coefficients

The subsurface shear-wave velocity (Vs) is considered to be a key parameter for site characterization and assessment of earthquake hazard because of its great influence on local ground-motion amplification. Array microtremor measurements are widely used for the estimation of shear-wave velocities. Compared to other methods such as frequency-wavenumber (f-k) methods, the spatial autocorrelation (SPAC) method requires fewer sensors and thus is relatively easier to implement and gives robust estimations of shear-wave velocity profiles for depths down to a few hundred meters. The quantity derived from observed data is the SPAC coefficient, which is a function of correlation distance, frequency and phase velocity. Generally, estimation of Vs profiles is a two stage process: Estimation of the dispersion data from the SPAC coefficients and inversion of the dispersion data for shear-wave velocity structure. In this study, instead of inverting dispersion curves, a more practical approach is used; that is, observed SPAC coefficients are directly inverted for the S-wave velocities. A synthetic case and a field data application are presented to test the potential of the inversion algorithm. We obtain an iterative damped least-squares solution with differential smoothing. The differential smoothing approach constrains the change in shear-wave velocities of the adjacent layers and thus stabilizes the inversion.     

River channel and bar patterns explained and predicted by an empirical and a physics‐based method

Our objective is to understand general causes of different river channel patterns. In this paper we compare an empirical stream power‐based classification and a physics‐based bar pattern predictor. We present a careful selection of data from the literature that contains rivers with discharge and median bed particle size ranging over several orders of magnitude with various channel patterns and bar types, but no obvious eroding or aggrading tendency. Empirically a continuum is found for increasing specific stream power, here calculated with pattern‐independent variables: mean annual flood, valley gradient and channel width predicted with a hydraulic geometry relation. ‘Thresholds’, above which certain patterns emerge, were identified as a function of bed sediment size. Bar theory predicts nature and presence of bars and bar mode, here converted to active braiding index (B_i). The most important variables are actual width–depth ratio and nonlinearity of bed sediment transport. Results agree reasonably well with data. Empirical predictions are somewhat better than bar theory predictions, because the bank strength is indirectly included in the empirical prediction. In combination, empirical and theoretical prediction provide partial explanations for bar and channel patterns. Increasing potential‐specific stream power implies more energy to erode banks and indeed correlates to channels with high width–depth ratio. Bar theory predicts that such rivers develop more bars across the width (higher B_i). At the transition from meandering to braiding, weakly braided rivers and meandering rivers with chutes are found. Rivers with extremely low stream power and width–depth ratios hardly develop bars or dynamic meandering and may be straight or sinuous or, in case of disequilibrium sediment feed, anastomosing. Copyright © 2010 John Wiley & Sons, Ltd.     

A non-contact mine pressure evaluation method by electromagnetic radiation

An electromagnetic radiation evaluation method for the relative stress state of coal bed under stress was proposed in this paper. The stress distribution of mine roadway or working face, as well as high stress zone or stress gradient zone, was analyzed by the method. The main advantages of the technique are its characteristics of non-contact, orientability, and regional monitoring. Correlation analysis of electromagnetic radiation with relative stress was carried out in coal mines and tunnels. The results indicate that the electromagnetic radiation technology has a wide application prospect in the evaluation of mine pressure.     

Estimation of coda wave attenuation in East Central Iran

The attenuation of coda waves, Q _c , has been estimated in Zarand, Jiroft, and Bam regions of east central Iran using a single back-scattering model of S-coda envelopes. For this purpose, the recordings of 97 earthquakes by three seismic networks and a local strong ground motion network have been used. In this research, the frequency-dependent Q _c values are estimated at central frequencies of 1.5, 3, 6, 8, 12, 16, and 24 Hz using different lapse time windows from 20 to 60 s. The frequency-dependent relationships obtained are for Zarand, for Jiroft, and for Bam region. From the strong ground motion data, we obtain the relation . The Q _c frequency-dependent relationship for the entire region of east central Iran from all data (both seismograms and accelerograms) is . The average Q _c values estimated and their frequency dependent relationships correlate well with a highly heterogeneous and highly tectonically active region. Results also show that the attenuation is higher in Bam region compared to Zarand and Jiroft regions.     

利用最小二乘自适应滤波实现绕射波分离

尺度目标体的绕射波能量较弱,常被能量较强的反射波淹没,与反射波一起成像难度较大.因而,实现绕射波成像的核心问题之一为绕射波与反射波分离.通常建立反射波模型的方法有Radon变换等域变换类方法,但该类方法在减去反射波能量上存在问题.为更彻底实现压制反射波目的,本文引入自适应滤波方法,该方法较简单的减去法能更好预测反射波能量.实际资料应用效果表明,最小二乘自适应滤波方法能更好突显单炮记录上绕射波特征,得到的叠加剖面中反射波去除更为彻底干净,偏移剖面中小断层、断点、尖灭点等小尺度地质体成像清晰.     

Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations.However, the computation cost generally increases linearly with increased order of accuracy.Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note th...     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

Estimation of source function and medium response function by autoregressive method

By modelling seismograms as “low” and “high” order autoregressive (AR) processes, the source function and the medium response function are separated from a single channel seismogram. Akaike's final prediction error is used as a statistic to select the appropriate “low” and “high” AR order of the process. Case studies of synthetic data show that the recovered source and reflectivity functions compare very well with the input functions. Using this method, arrivals of the surface reflected P phases of five explosions from the Soviet region and of two earthquakes from Kamchatka, recorded at Gauribidanur Seismic Array (GBA), India, are identified. Certain features of the source and source region of these events are also inferred.     

An empirical method for approximating canopy throughfall

Rainfall replenishes surface and subsurface water but is partially intercepted by a canopy. However, it is challenging to quantify the rainfall passing through the canopy (i.e. throughfall). This study derives simple‐to‐use empirical equations relating throughfall to canopy and rainfall characteristics. Monthly throughfall is calculated by applying a mass balance model on weather data from Singapore; Vancouver, Canada; and Stanford, USA. Regression analysis is then performed on the calculated throughfall with three dependent variables (i.e. maximum canopy storage, average rainfall depth and time interval between two consecutive rainfall in a month) to derive the empirical equations. One local equation is derived for each location using data from that particular location, and one global equation is derived using data from all three locations. The equations are further verified with calculated monthly throughfall from other weather data and actual throughfall field measurements, giving an accuracy of about 80–90%. The global equation is relatively less accurate but is applicable worldwide. Overall, this study provides a global equation through which one can quickly estimate throughfall with only information on the three variables. When additional weather data are available, one can follow the proposed methodology to derive their own equations for better estimates. Copyright © 2012 John Wiley & Sons, Ltd.     

Wavelet packet polarization method for an automatic detection of P and S wave arrivals

A new method is proposed for detecting P and S wave arrivals in the problem of passive seismic monitoring. This method is based on the wavelet packet decomposition of three-component seismic traces, their automatic quality control, and the calculation of the principal components of the traces in overlapping packet frequency bands and in scale-dependent moving time windows. The wave arrivals times are estimated by means of a robust iterative procedure of adjusting hyperbolic traveltime curves to the initial estimates of arrival times determined from the maximum values of the multilevel measure of the nonstationarity of each trace. The method has been developed for applications in the cases of a large number of bad single-component traces and a high noise level. To illustrate its application, the method is applied to the data set of the Cotton Valley experiment on the seismic monitoring of hydraulic fractures in a fractured gas reservoir.