VARIATION OF THE TRAVEL TIME DIFFERENCE BETWEEN TELESEISMIC PS CONVERTED WAVE AND PP TRANSMITTED WAVE IN FOCAL REGION BEFORE AND AFTER WENCHUAN MS8.0 EARTHQUAKE

When P waves from distant earthquakes meet a velocity discontinuity in the earth's crust and upper mantle, they give rise to a series of converted PS waves besides PP refracted waves. It is possible to monitor the variation of the physical properties of the medium in the limited formation space above the transition zone of the seismogenic zone by measuring the time difference between the teleseismic PS converted wave and the first arrival P wave, that is, time-variation Δt_PS=t_PS-t_P. The advantage of this method is that the transition point of the teleseismic source with similar source is relatively stable at the transition interface, and the accuracy of the measured relative time is high, and the change of the medium in a small range of the seismogenic zone above the conversion interface can be monitored. This paper studies the variation of the travel time difference Δt_PS in focal region before and after Wenchuan M_S8.0 earthquake. We select 2001 to 2012 as the research period, use teleseismic waveforms which occurred in the southern region of Sumatra and Hindu Kush area recorded by Sichuan station YZP and JJS. These teleseisms satisfy 5.0 ≤ M ≤ 6.5, and their waveform signal-to-noise ratio is high with clear initial P-wave motion. The epicentral distance of teleseisms is less than 3 degrees. Then we obtain the variation of the travel time difference Δt_PS between teleseismic PS converted wave and PP transmitted wave recorded during the study period in the two stations. The results show that there is a slow increasing trend of Δt_PS before 2006, and an obvious low value process of Δt_PS appeared in the period about 2 years before the Wenchuan earthquake. The maximum decline was about 0.2~0.3s, more than 4~5 times the measurement error. The low value has a certain degree of return about 2~3 months before the earthquake. The change of arrival time difference indicates that the medium is in different states in different periods of seismogenic process. The sharp decrease in Δt_PS from 2006 to January 2008 may be due to the strong disturbance caused by the stress accumulation of the medium. At this stage, the velocity of P wave and S wave increases with the increase of stress, and the increase of S wave velocity will result in the decrease of Δt_PS. The change of Δt_PS is greatly affected by S wave velocity, so Δt_PS appears to decrease rapidly. Regarding the low value that has a certain degree of return about 2~3 months before the Wenchuan earthquake, the possible reason is that the release of stress is much higher than the accumulation of stress in meta-instability stage. At this stage, the velocity of S wave decreases and the decrease of S wave causes Δt_PS to increase. Then, the Wenchuan earthquake of magnitude 8.0 occurred. It is shown that the teleseismic converted wave method in this paper can monitor the variation of medium's wave velocity before large earthquakes, and it has a good prospect in seismic monitoring and worth further experimental study.     

2017年9月4日河北临城M_L4.4地震震源深度测定

震源深度是地震定位中的重要参数之一。对于一个地震,sPn与Pn震相到时差是一个常量,与震中距无关,只与震源深度和地壳模型有关。基于该理论,以河北数字台网数据为基础,识别sPn震相和Pn震相,计算2017年9月4日河北临城ML4.4地震震源深度,结果为6.86km。双差定位法是目前确定震源深度比较精确的方法,双差定位得到的震源深度为7.01km。双差定位法所用数据量比较大,计算过程繁杂。利用sPn与Pn到时差计算震源深度的结果与双差定位法结果接近,计算快速、简单,所需数据量小。因此,将该方法应用到实际工作中,可以快速、准确地测定近震震源深度。     

断裂带首波研究进展

断裂带首波是沿着存在物性差异界面传播的一种地震折射波,在传播过程中携带了断裂带的重要信息,对分析和研究断裂带以及附近区域的精细结构提供了一种新的分析方法。本文主要阐述断裂带首波的产生原理、波形特征、识别及分析方法等,介绍目前国际上识别及利用断裂带首波开展断裂带特征方面的研究现状,并针对地震危险区域存在物性差异的断裂带,提出可结合密集台阵观测技术,利用断裂带首波进行断裂带精细结构探测及其变化监测研究,提高潜在孕震环境及发生机理的认识水平。     

京西北地区地震重定位及构造特征

采用双差定位法对京西北地区（39.5°—41.5° N,113°—117° E）2013年1月1日至2017年12月31日6 223次有效地震进行精确定位,得到该区震源分布的精细图像和震源深度剖面图。结果显示,重新定位后地震的水平分布更集中,沿断裂带分布特征更加明显,震中在断裂带呈更明显的条带状、簇状分布,地震与线状的深浅断裂构造的关系密切;大部分地震发生在中上地壳,震源分布为典型震源密集区的空间形态,呈纵深约15 km、直径20—40 km的近圆形“厚饼状”。     

湖泊三维风生流隐式差分模型的研究

建立了湖泊三维风生流的隐格式差分模型,通过对该模型计算稳定性的分析发现,它的时间步长的取值具有较大的灵活性．该模型波用于武汉市墨水湖的民生流模拟,实例计算表明该模型计算结果合理,能较好地反映湖泊风生流的流态特征．     

Pathlines of pollutant migration in groundwater of the Visakhapatnam urban area,India

The groundwater regime in the north‐western part of the Visakhapatnam urban area was polluted as early as 1981 by discharge of untreated industrial effluent from a Hindustan Polymers Limited (HPL) plant. A total dissolved solids (TDS) concentration of surficial effluent up to 6500 mg/l and of groundwater in the range 3000–4200 mg/l has been reported in the environs of the HPL plant during May 1992. Groundwater occurs under water table conditions. The groundwater flow model was simulated as a single layer aquifer, using MODFLOW and FLOWPATH computer models. Aquifer parameters were estimated and the value of effective porosity assumed. The water table configuration was mapped in May 1981 using 33 observation wells. The calibrated hydraulic head distribution was used to compute the velocity field using the effective porosity values. The flow paths of groundwater migration from the source locations up to 2002 AD were predicted. Pathlines of particles in the groundwater indicated a predominant north‐east and south‐west migration of groundwater pollution in the area. Copyright © 1999 John Wiley & Sons, Ltd.     

两个30年气候平均值的差异及其对气候业务的影响

对陕西短期气候预测业务使用的39个气象站气候要素1971～2000年平均值与1961～1990年平均值进行比较,发现大部分地区年平均气温升高,晴天日数增加;年降水量、阴天日数、降水日数、年日照时数减少;年大风日数、沙尘暴日数、冰雹日数、雷暴日数减少等。把两个平均值放在长序列中分析,发现某些气象要素最近的气候平均值的差值在20世纪80～90年代期间的5～6个气候平均值每两个相邻平均值的差值中是最大的,说明陕西20世纪90年代气候较60年代干燥、温暖,且90年代气温高、降水少的特点十分突出。还对某些要素平均值改变对气候影响评价、气候预测业务产生的影响进行简要分析。     

Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.     

HY-1A卫星COCTS数据条带消除的两种定量化方法比较

HY-1A卫星是中国第一颗海洋遥感卫星,水色水温扫描仪(COCTS)是其主要探测器。由于各探测单元的信号响应特性不同,导致COCTS数据含有大量条带噪声。这严重影响了图像的质量和可视性,对资料的定量化应用也是十分不利的。为了解决这个问题,本研究采用矩匹配及直方图匹配法,对COCTS数据进行条带消除。结果表明,两种方法都是有效的,能够在消除条带的同时保持数据原有的物理意义,并避免细节特征的丢失。但是从整体效果看,直方图匹配法更好一些。本研究的结果同样适用于其他以并扫方式获取数据的遥感器的条带消除。     

京津风沙源区2000—2015年植被覆盖状况的区域差异研究

植被覆盖状况是监测与表征区域生态建设成效的重要指标。以往研究注重区域植被覆盖状况的整体变化分析,对长时期区域内部植被覆盖变化及其差异研究较少。文中基于京津风沙源区遥感影像数据,采用GIS空间分析技术,重点评估了2000—2015年植被覆盖度的年际变化及其区域差异。研究结果表明,2000—2015年京津风沙源区植被覆盖度变化为35%~45%,且随年份变化呈波动增加趋势,年均增速为0.4%(P<0.05),生态治理取得明显植被恢复成效。8个治理分区植被覆盖均有所增加,但区域差异明显,晋北山地丘陵亚区和燕山丘陵山地水源保护亚区植被覆盖度年均增速超过0.6%,浑善达克沙地亚区与荒漠草原亚区植被覆盖度年均增速不及0.2%,这与区域地表组成和气候背景有关。从地市来看,北京、天津和承德植被覆盖度较高,但朔州和张家口植被覆盖度增速明显,而乌兰察布和包头植被覆盖度年均增速低于0.1%。相比2000年,2015年京津风沙源区有51%的区域植被覆盖度增加,49%区域植被覆盖度未变或降低,主要集中在京津风沙源区的中部和西部县市(旗),未来生态治理过程中应加以重点关注。