Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

Non-point source pollution estimation in the Pingqiao River Basin,China, using a spatial hydrograph-separation approach

For efficient and targeted management, this study demonstrates a recently developed non-point source (NPS) pollution model for a year-long estimation in the Pingqiao River Basin (22.3 km²) in China. This simple but physically reasonable model estimates NPS export in terms of land use by reflecting spatial hydrological features and source runoff measurements under different land-use types. The NPS export was separately analysed by a distributed hydrological model, a spatial hydrograph-separation technique, and an empirical water quality sub-model. Simulation results suggest that 57 890 kg of total nitrogen (TN) and 1148 kg of total phosphorus (TP) were delivered. The results, validated with observed stream concentrations, show relative errors of 23.3% for TN and 47.4% for TP. Countermeasures for urban areas (5.3% of total area) were prioritized because of the high contribution rate to TN (14.1%) and TP (26.2%) which is caused by the high degree of runoff (8.5%) and pollution source.     

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到时差计算震源深度的结果与双差定位法结果接近,计算快速、简单,所需数据量小。因此,将该方法应用到实际工作中,可以快速、准确地测定近震震源深度。     

地磁低点时间与经度的相关性分析

地磁低点时间在空间分布上的主要特征是随着经度变化而变化。为了求出这种变化关系的解析式,首先使用国家地磁台网的数据计算甘肃省近三年13个地磁台站的低点时间序列和低点时间的期望,并绘制了低点时间序列的频次分布图和概率密度曲线;然后对低点时间与台站经度的关系做线性回归,并指出这个线性回归方程应无限趋近于“格林尼治时间-经度”关系式。文章第一次提出台站经度与地磁Z分量低点时间的关系式,并提出它的趋近方程。各台站低点时间的期望就是地磁低点时间的正常背景值,它是判定地磁低点时间异常与否的参考标准。本文的工作完善和发展了丁鉴海先生的地磁低点位移法。     

固结时间对软土动力特性影响的试验研究

固结试验是研究软土动力特性的重要环节,固结时间是影响软土动力特征参数的重要因素。围绕岩土工程及科研工作中对动三轴试验中固结时间影响规律的认识较为模糊的问题,拟对天津滨海软土开展动三轴平行试验,分别采用8h、12h、24h、48h进行固结,针对软土的最大动剪切模量G_max、最大阻尼比λ_max及动剪切模量比G_d/G_max和阻尼比λ随动剪切应变γ_d的变化规律,分析固结时间对软土动力特性参数的影响。试验结果表明:G_max、λ_max、G_d/G_max-γ_d曲线、λ-γ_d曲线受固结时间的影响较为显著。G_d/G_max-γ_d曲线随着固结时间的增加而减小。阻尼比试验结果表明动剪切应变存在明显的分界点,分界点以下,中长期固结时间所得阻尼比较大,而分界点以上,短期固结时间所得阻尼比较大。     

逆时偏移在声反射成像测井中的应用方法研究

声反射成像测井中常用的基于射线理论的Kirchhoff积分偏移算法和基于单程波理论的F-K偏移算法均可实现井旁缝洞反射体的快速偏移成像,但其仅适用于地层垂向变化较弱的速度场和高陡角裂缝的偏移成像,无法实现低角度反射体的准确偏移归位,产生偏移假象误导测井解释.逆时偏移基于全波动方程,可适应强垂向变化速度场,实现近似水平反射体的偏移成像.本文详细分析了将逆时偏移应用于声反射成像测井时存在的数据准备、时间采样间隔匹配和成像条件改进等若干问题,通过设置多组理论模型来说明算法对井旁不同反射体的识别能力.模拟资料和实际资料处理结果证实,较F-K偏移算法,逆时偏移算法成像精度更高、收敛性更好,可有效实现近似水平构造偏移归位.改进的归一化互相关成像条件可解决深部地层的远井壁成像衰减问题,降低测井解释的多解性.逆时偏移将成为声反射成像测井高精度偏移技术的发展方向.     

主动震源在中国区域尺度及水库地震监测中的应用

综述精密可控震源与陆上气枪2种主动震源特性及在国内区域尺度上的应用,估算观测介质波速随时间变化的精度,探讨主动震源在构建水库地震观测介质波速随时间变化中的系统应用。     

关于土层地震反应分析中阻尼矩阵建模方式的讨论

阻尼矩阵的选择对土层时域内地震反应计算的准确性有重要意义。本文构造了质量比例阻尼、刚度比例阻尼、Rayleigh阻尼、Caughey阻尼和Clough阻尼矩阵,研究各阻尼对土层地震反应分析精度的影响。以苏通大桥5号桥塔基础处的工程场地为例,分别以人工合成基岩波,Northridge波、Parkfield波和汶川波为输入,得到不同阻尼模型下土层的地震反应,结果表明Clough阻尼在计算精度和误差稳定性上要优于其他阻尼模型,针对Clough阻尼的不足之处提出了改进方法,算例表明改进后的Clough阻尼能在保证计算精度的前提下减少计算机的存储量和计算量。