基于不同地形的GPS对流层延迟插值方法研究

通过利用北京、武汉两地的GPS数据,采用反距离加权法、克里金法研究不同地形下GPS对流层延迟插值效果.经比较,反距离加权法与克里金法的插值效果基本一致,在对武汉地区进行GPS对流层延迟插值时插值效果均在毫米级,在北京地区的插值效果则有不同:平坦地区插值效果达到毫米级;山区插值效果为厘米级,无法满足应用需求.     

Finite-Difference Modeling and Dispersion Analysis of High-Frequency Love Waves for Near-Surface Applications

Love-wave propagation has been a topic of interest to crustal, earthquake, and engineering seismologists for many years because it is independent of Poisson’s ratio and more sensitive to shear (S)-wave velocity changes and layer thickness changes than are Rayleigh waves. It is well known that Love-wave generation requires the existence of a low S-wave velocity layer in a multilayered earth model. In order to study numerically the propagation of Love waves in a layered earth model and dispersion characteristics for near-surface applications, we simulate high-frequency (>5 Hz) Love waves by the staggered-grid finite-difference (FD) method. The air–earth boundary (the shear stress above the free surface) is treated using the stress-imaging technique. We use a two-layer model to demonstrate the accuracy of the staggered-grid modeling scheme. We also simulate four-layer models including a low-velocity layer (LVL) or a high-velocity layer (HVL) to analyze dispersive energy characteristics for near-surface applications. Results demonstrate that: (1) the staggered-grid FD code and stress-imaging technique are suitable for treating the free-surface boundary conditions for Love-wave modeling, (2) Love-wave inversion should be treated with extra care when a LVL exists because of a lack of LVL information in dispersions aggravating uncertainties in the inversion procedure, and (3) energy of high modes in a low-frequency range is very weak, so that it is difficult to estimate the cutoff frequency accurately, and “mode-crossing” occurs between the second higher and third higher modes when a HVL exists.     

东北三省耕地物候期对热量资源变化的响应

选取1998-2010 年期间的SPOT/VGT NDVI逐旬时间序列数据、物候观测数据和气候数据,利用函数拟合方法、GIS空间分析和统计方法,分析了东北三省不同耕地物候期对热量资源变化的响应关系。首先,基于拟合数据提取了耕地物候期特征(包括生长季开始期、峰值期、结束期与生长季长度等),验证结果显示与不同作物的观测物候期存在明显对应关系,能较好地反映作物物候期特征的宏观时空分异。其次,分析了1991-2009 年期间区内主要农业热量资源(≥ 10 oC初日,初霜日,温度生长期天数和≥ 10 oC积温) 的时空变化特征。最后,重点分析了1998-2009年期间耕地物候期特征对热量资源变化的响应。结果表明,生长季开始期与≥10 oC初日在松嫩平原北部、三江平原东部和吉林省中东部的正相关系显著(P < 0.01);生长季长度与温度生长期天数亦在上述区域内表现为显著正相关关系(P < 0.05);然而,生长季结束期与初霜日的相关关系整体上较不明显;此外,生长季长度与≥ 10 oC积温的相关关系存在两种相反的情形,在松嫩平原西北部、三江平原、吉林省中部和辽宁省中部的负相关关系显著(P < 0.05),而黑龙江省中南部和三江平原部分区域,两者相关系数达到0.55 (P < 0.05),分别反映了作物生长发育周期的缩短态势和中早熟作物被中晚熟作物替代的趋势。     

2013年吉林前郭M_S5.8震群序列研究

对前郭M_S5.8地震后6个月(2013年10月31日~2014年4月30日)内的地震事件进行了重新定位,并利用吉林区域地震台网提供的地震目录研究了序列活动特征。精定位结果显示,余震区呈NW走向分布,长轴约11km,短轴约6km。震源深度为5~15km。结合区域构造特征认为,该序列可能与NW走向的通榆-长春隐伏断裂活动有关,序列具有b值较低、5级地震频次高、发震时间集中的特点。序列主要起伏活动可以分为3个时段:主震后10天,余震强度衰减不明显;此后的2次起伏活动,地震活动分别表现出平静-增强(震级爬升)-发生强震和增强(震级爬升)-平静-再增强(震级爬升)-发生强震的特点;在序列早期阶段,震级-频度无法拟合成1条直线,G-R关系在低震级段和高震级段呈现2个线性段;在M_S5.8地震前h值有变小的趋势,且接近1;而M_S5.8地震后h值明显增大;M_S5.5地震后15天序列衰减比较慢,p值为0.76;M_S5.8地震后15天序列衰减较快,p值为1.17。     

背景场项目中iDirect卫星通信系统设计

介绍背景场项目中i Direct卫星通信系统,主要从背景场项目对卫星通信的需求、系统设计、数据接入、系统应用等方面进行阐述,具体介绍卫星信系统为行业网备份的数据接入方案以及数据应用案例,可独立为国家测震台提供数据传输信道,保障测震速报业务的正常开展。     

Experimental analysis on the impact force of viscous debris flow

A miniaturized flume experiment was carried out to measure impact forces of viscous debris flow. The flow depth (7.2–11.2 cm), velocity (2.4–5.2 m/s) and impact force were recorded during the experiment. The impact process of debris flow can be divided into three phases by analyzing the variation of impact signals and flow regime. The three phases are the sudden strong impact of the debris flow head, continuous dynamic pressure of the body and slight static pressure of the tail. The variation of impact process is consistent with the change in the flow regime. The head has strong–rapid impact pressure, which is shown as a turbulent‐type flow; the body approximates to steady laminar flow. Accordingly, the process of debris flows hitting structures was simplified to a triangle shape, ignoring the pressure of the tail. In order to study the distribution of the debris flow impact force at different depths and variation of the impact process over time, the impact signals of slurry and coarse particles were separated from the original signals using wavelet analysis. The slurry's dynamic pressure signal appears to be a smooth curve, and the peak pressure is 12–34 kPa when the debris flow head hits the sensors, which is about 1.54 ± 0.36 times the continuous dynamic pressure of the debris flow body. The limit of application of the empirical parameter α in the hydraulic formula was also noted. We introduced the power function relationship of α and the Froude number of debris flows, and proposed a universal model for calculating dynamic pressure. The impact pressure of large particles has the characteristic of randomness. The mean frequency of large particles impacting the sensor is 210 ± 50–287 ± 29 times per second, and it is 336 ± 114–490 ± 69 times per second for the debris flow head, which is greater than that in the debris flow body. Peak impact pressure of particles at different flow depths is 40–160 kPa, which is 3.2 ± 1.5 times the impact pressure of the slurry at the bottom of the flow, 3.1 ± 0.9 times the flow in the middle, and 3.3 ± 0.9 times the flow at the surface. The differences in impact frequency indicate that most of the large particles concentrate in the debris flow head, and the number of particles in the debris flow head increases with height. This research supports the study of solid–liquid two phase flow mechanisms, and helps engineering design and risk assessment in debris flow prone areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Recursive impedance inversion of ground-penetrating radar data in stochastic media

The travel time and amplitude of ground-penetrating radar (GPR) waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data.     

频率域航空电磁法地形影响和校正方法

为了研究地形对频率域航空电磁法的影响程度,以及如何有效提取由于地形产生的航空电磁场,笔者利用频率域有限差分法开展了复杂地形条件下的频率域航空电磁响应计算和分析研究。计算结果表明,坡度较陡的地形对频率域航空电磁响应影响较大(尤其在沟谷位置地形影响最大),以至于增加了电磁异常的复杂程度,影响了电磁异常的解释结果。在综合研究的基础上,采用二维和三维的地形校正方法对典型地形地电模型进行校正,有效去除和压制了地形引起的电磁干扰异常,突出了异常体引起的电磁异常,提高了频率域航空电磁的解释效果。     

Comparison of ground truth location of earthquake from InSAR and from ambient seismic noise:A case study of the 1998 Zhangbei earthquake

Because ambient seismic noise provides estimated Greenos function (EGF) between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB (a station about 40 km away from the epicenter) and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Greenos functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Greenos functions with EGFs in period band of 10-30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 (ground truth location with accuracy 5 km) with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.