基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移

近年来,面向实际应用的TI介质准P波正演模拟与逆时偏移成像技术受到空前的关注.基于常规耦合型传播方程的正演模拟方法不仅存在伪横波及频散假象干扰,而且还遭受模型参数限制(η0)和不稳定影响;而纯qP波方程的推导繁琐,且由于方程中包含拟微分算子造成求解难度大且精度有限.为此,本文首先构建了一种适用于任意TI介质的纯qP波传播算子,然后借助Low-rank分解求取该算子中的空间-波数域矩阵,同时引入Cerjan衰减边界条件来压制边界反射干扰,最终实现了一种间接的纯qP波波场外推方案,并将其成功应用于复杂TI介质正演模拟与逆时偏移成像中.通过开展数值模拟,并与其他方法对比表明:①该方法既避免了纯qP波方程的繁琐推导,又克服了耦合型方程对模型参数的限制;②还彻底消除了残余伪横波噪音及数值频散;③且能适应较大时间或空间步长及高频震源,是一种相对准确且稳定的各向异性纵波正演与成像策略.     

水平层状非均质TI地层中仪器偏心情况下三维感应测井响应高效数值模拟与响应特征分析

本文针对中国石油集团测井有限公司(CPL)三维阵列感应测井仪器(Three Dimensional Induction Tool,TDIT)的真实结构与参数,利用2.5维数值模式匹配技术研究建立水平层状非均质TI地层中仪器偏心情况下三维感应测井响应的高效数值模拟算法.首先,利用2.5维数值模式匹配技术给出水平层状非均质TI模型中偏心磁流源并矢Green函数的半解析解.在此基础上,结合三维阵列感应测井仪器参数以及测量过程中出现的仪器偏心与旋转等实际情况,通过叠加原理推导出该仪器响应的有效计算公式.最后,通过数值模拟结果分析考察仪器偏心、仪器自旋角等参数对测井响应的影响以及变化特征,为多分量感应资料处理与解释提供理论基础.     

大同地震窗P波频散衰减的异常变化特征

通过测量介质品质因子来开展大同窗附近测震台站记录P波衰减频散特征的研究,并分析大同窗近场区和外围区震例前后的变化。对山自皂台记录的2001年以来大同窗地震序列的研究结果表明,近场区3次和外围区6次地震前,介质品质因子显示超出1倍均方差的高值变化特征,震后逐渐恢复到背景起伏状态;与大同窗传统的频次指标相比,介质品质因子的映震效果更好;阳原台的结果表明,仅近场区1次和外围区2次地震前,介质品质因子显示超出1倍均方差的高值变化特征,且有虚报、漏报。所处地质构造位置不同、传播路径介质差异可能是引起上述两个台站附近介质品质因子的地震对应效果明显不同的原因。     

中、强震前川滇块体遥感红外辐射异常分析——以宁蒗M5.7地震和鲁甸M6.5地震为例

许多研究表明中强地震前出现了红外辐射异常。近年来,川滇块体中、强地震频发,利用遥感数据对其红外辐射异常变化开展深入分析十分必要。结合13年MODIS/Terra遥感红外数据,采用空间距平分析法、辐射增强阈值面积指数法及辐射亮温均值低频时间序列法,对2004年4月—2017年6月红外亮温数据进行深入分析,研究川滇块体及周边构造断裂M≥5.5地震的震前异常信息。结果表明:(1)M≥6.0地震在震前3～6个月川滇块体会出现较为明显的红外辐射增强异常现象;(2)震前红外异常并非随机分布,而可能是受到相关块体的控制;(3)以宁蒗M5.7地震和鲁甸M6.5地震为例研究发现,震前辐射增强较高的小范围区域与震中位置存在一定关联。     

A continuum model of layered rock masses with non-associative joint plasticity

Layered rock masses can be modelled either as standard, orthotropic continua if the layer bending can be neglected or as Cosserat continua if the influence of layer bending is essential. This paper presents a finite element smeared joint model based on the Cosserat theory. The layers are assumed to be elastic with equal thickness and equal mechanical properties. All the cosserat parameters are expressed through the elastic properties of layers, layer thickness and joint stiffness. Plastic-slip as well as tensile-opening of layer interface (joint) are accounted for in a manner similar to the conventional non-associative plasticity theory. As an application, the behaviour of an excavation in a layered rock mass is examined. The displacement and stress fields given by smeared joint models based on the Cosserat continuum and the conventional anisotropic continuum approaches are compared with those obtained from the discrete joint model. The conventional anisotropic continuum model is found to break-down completely when the effective shear modulus in the direction parallel to layering is low in comparison to the shear modulus of the intact layer, whereas the Cosserat model is found to be capable of accurately reproducing complex load–deflection patterns irrespective of the differences in shear moduli. © 1998 John Wiley & Sons, Ltd.     

基于Ruger近似式预测裂缝型储层的可行性研究

叠前纵波方位各向异性检测裂缝是目前应用最为广泛的方法,该方法主要是利用AVAZ地震资料椭圆拟合的长轴与短轴信息来评价裂缝,该方法在具体应用时,认为裂缝型储层AVAZ响应特征的变化仅由各向异性参数决定,而没有考虑非各向异性参数的影响,从而引起裂缝评价的多解性。为提高该技术的裂缝评价精度,从模型分析入手,将模型中各向异性参数及岩性参数设置为概率密度分布函数,然后采用Monte Carlo随机方法进行叠前AVAZ正演模拟。首先,对各向异性参数对椭圆拟合的影响进行简单分析,得出各向异性参数γ对椭圆扁率B/A与各向异性因子B影响最大、δ次之、ε最小。然后,进一步重点分析不同标准差的速度与密度模型的B/A与B响应特征,得出地层的速度与密度的变化在较小范围内时,该技术评价裂缝才具有一定的可行性,其中,地层纵波和横波速度的变化（尤其是纵波速度）对叠前AVAZ响应影响较大,地层密度几乎没有影响,故当横向岩性变化较大时,该方法在评价裂缝型储层时的结果具有一定的不确定性。      

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年8月冀鲁豫交界区流动地磁总强度数据进行处理并解析延拓,研究该区岩石圈磁场及其深部、浅表磁场异常分布特征与地下介质磁化率结构的性质。结果表明:冀鲁豫交界区岩石圈磁场异常分布是岩石圈浅表与深部磁异常综合叠加的结果,岩石圈地下介质的浅表与深部磁化率结构具有一定差异性,且受断裂控制呈分段性特征,可能与该区域复杂的岩石圈磁化率结构及地质构造有关。     

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

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

衬砌隧道中移动轴向激励作用下两相多孔介质动力响应

为了研究衬砌隧道中地铁振动在饱和地层中的传播情况,采用移动轴向激励模拟技术,建立了隧道-衬砌-两相介质的动力分析模型;利用波函数展开法、傅里叶变换法等,推导了频域内衬砌隧道移动轴向激励作用下两相多孔介质动力响应的解析解,并给出了两相多孔介质临界速度的经验公式;通过离散快速傅里叶逆变换得到了时-空域内两相介质的动力响应。结果表明：无衬砌隧道中移动的轴向常激励作用下,两相介质临界速度只与介质的剪切模量和密度有关,且数值接近其剪切波速的1.1倍;对于衬砌隧道,介质的临界速度随着衬砌剪切模量的增大而增大,随着衬砌密度的增大而减小,衬砌对振动向介质中的传播有一定的削弱作用,衬砌剪切模量与介质剪切模量相差越大,削弱越明显;动力响应频率越接近激振频率,其幅值越大,所对应的临界速度越小。