Low-frequency symmetric waves in fluid-filled boreholes and pipes with radial layering

Many tasks in geophysics and acoustics require estimation of mode velocities in cylindrically layered media. For example, acoustic logging or monitoring in open and cased boreholes need to account for radial inhomogeneity caused by layers inside the borehole (sand screen, gravel pack, casing) as well as layers outside (cement, altered and unaltered formation layers). For these purposes it is convenient to study a general model of cylindrically layered media with inner fluid layer and free surface on the outside. Unbounded surrounding media can be described as a limiting case of this general model when thickness of the outer layer is infinite. At low frequencies such composite media support two symmetric modes called Stoneley (tube) and plate (extensional) wave. Simple expressions are obtained for these two mode velocities valid at zero frequency. They are written in a general form using elements of a propagator matrix describing axisymmetric waves in the entire layered composite. This allows one to apply the same formalism and compute velocities for n -layered composites as well as anisotropic pipes. It is demonstrated that the model of periodical cylindrical layers is equivalent to a homogeneous radially transversely isotropic media when the number of periods increases to infinity, whereas their thickness goes to zero. Numerical examples confirm good validity of obtained expressions and suggest that even small number of periods may already be well described by equivalent homogeneous anisotropic media.     

不同模式自组织介质中声波传播特性的比较研究

中型和大型尺度结构构造分区性, 使得地球内部介质地震波速度表征出大尺度的确定性模型, 而小尺度上微孔和裂隙的存在会导致地震波速度分布的随机扰动, 即自组织模型. 目前地震学中常用的6种模式自组织模型为: 高斯型、指数型、自相似型、白噪声、Flicker噪声和Brown 噪声. 我们从其滤波因子入手, 比较了不同模式的自组织模型特征, 结合有限差分法地震记录模拟及地震波特征分析, 对比了上述6种模式自组织介质中声波传播特性. 结果表明：地震波对不同模式的自组织介质的地震响应不同, 如波形、能量、振幅等. 由此我们可以通过分析实际地震记录的复杂性特征来鉴别其自组织类型, 从而更好地解释实际地震资料.     

偶应力理论框架下的弹性波数值模拟与分析

在经典地震学理论框架下,先人发展了数不胜数的地震技术,为社会发展做出了巨大贡献;可是,当前的技术仍有难以逾越的障碍,亟需破局.高铁地震学联合研究组,在河北定兴采集到大量数据;其中可见含有大量的旋转运动分量.由于基于经典连续介质力学推导弹性波动方程时,从理论出发点上就去除了旋转项,且在其理论框架内,介质被视为一个连续的质量体.而事实上,不论是人造还是天然的介质都存在着复杂的内部结构,广义连续介质力学更适合描述具有更加复杂内部结构的情况.于是,我们尝试启用广义连续介质力学理论,推导偶应力理论框架下的弹性波动方程;将其数学表达形式以及数值模拟结果与传统弹性波动方程进行对比.研究探索推导的具非对称性的波动方程所具有的理论意义和应用价值.     

基于弹性波全波形反演的主被动源多分量混采地震数据速度建模

在常规地震采集中,被动源地震波场往往被视为噪声而去除,这就造成了部分有用信息的丢失.在目标区进行主动源和被动源弹性波地震数据的多分量混合采集,并对两种数据进行联合应用,使其在照明和频带上优势互补,能显著提高成像和反演的质量.本文针对两种不同类型的主被动源混采地震数据,分别提出了相应的联合全波形反演方法.首先,针对主动源与瞬态被动源弹性波混采地震数据,为充分利用被动源对深部照明的优势,同时有效压制被动震源点附近的成像异常值,提出了基于动态随机组合的弹性波被动源照明补偿反演策略.然后,针对低频缺失主动源与背景噪声型被动源弹性波混采地震数据,为充分利用被动源波场携带的低频信息,并避免对被动源的定位和子波估计,提出了基于地震干涉与不依赖子波算法的弹性波主被动源串联反演策略.最后,分别将两种方法在Marmousi模型上进行反演测试.结果说明,综合利用主动源和被动源弹性波混采地震数据,不仅能增强深部弹性参数反演效果,还能更好地构建弹性参数模型的宏观结构,并有助于缓解常规弹性波全波形反演的跳周问题.     

小湾库区水库诱发地震的地质环境分析

在文献[1]的基础上,通过对小湾库区岩石介质、地质构造、构造应力场、库区及周围地震活动特点的分析,认为小湾库区水库诱发地震的地质环境比较复杂,存在发生岩溶型、裂隙型、构造型等三种成因类型的水库诱发地震的基本条件,而沿通过库区的澜沧江断裂带及其次级断层,和无量山断裂带发生构造型水库诱发地震,是小湾库区水库地震活动的主要诱发因素：     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

MSH法直接将地球勘探测线转化为深度,密度剖面

根据运动学中的动量守恒定律、分析地震勘探过程中地震波在地下介质中的传播过程,利用反射能量、透射能量与人射能量的计算公式和检波器转化的电信号与地下密度界面反射地震波波动能量的等效原理,推导出直接求层速度和层密度公式。利用某一反射界面直接测得的波速或依据钻孔资料求得的平均速度,在计算机上就可直接计算出各地震测线物性层的埋深和密度（速度、波阻抗）,由此地质工作人员依据已知地质柱状剖面（钻孔资料）可进一步     

库车坳陷高陡构造地震勘探复杂性定量分析

波与非均匀介质的相互作用是多年来地球物理研究的理论问题之一.本文基于地震成像算子,定量分析了库车坳陷盐相关逆冲推覆构造的复杂性.首先,根据地震资料计算了大北构造、博孜构造、却勒构造、西秋10构造和西秋4构造的速度横向变化非均质谱和地层倾角变化非均质谱,这些量化表征高陡构造地质复杂性的地质非均质谱定最描述了速度横向变化和...     

Non‐iterative equivalent linearization of inelastic SDOF systems for earthquakes in Japan and California

The seismic performance of existing structures can be assessed based on nonlinear static procedures, such as the Capacity Spectrum Method. This method essentially approximates peak responses of an inelastic single‐degree‐of‐freedom (SDOF) system using peak responses of an equivalent linear SDOF model. In this study, the equivalent linear models of inelastic SDOF systems are developed based on the constant strength approach, which does not require iteration for assessing the seismic performance of existing structures. To investigate the effects of earthquake type and seismic region on the equivalent linear models, four ground‐motion data sets—Japanese crustal/interface/inslab records and California crustal records—are compiled and used for nonlinear dynamic analysis. The analysis results indicate that: (1) the optimal equivalent linear model parameters (i.e. equivalent vibration period ratio and damping ratio) decrease with the natural vibration period, whereas they increase with the strength reduction factor; (2) the impacts of earthquake type and seismic region on the equivalent linear model parameters are not significant except for short vibration periods; and (3) the degradation and pinching effects affect the equivalent linear model parameters. We develop prediction equations for the optimal equivalent linear model parameters based on nonlinear least‐squares fitting, which improve and extend the current nonlinear static procedure for existing structures with degradation and pinching behavior. Copyright © 2010 John Wiley & Sons, Ltd.     

High-order generalized screen propagator migration based on particle swarm optimization

Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unresolved. To describe the seismic wave propagation in media with lateral velocity variations and to image high-angle structures, we propose the generalized screen propagator based on particle swarm optimization (PSO-GSP), for the precise fitting of the single-square-root operator. We use the 2D SEG/EAGE salt model to test the proposed PSO-GSP migration method to image the faults beneath the salt dome and compare the results to those of the conventional high-order generalized screen propagator (GSP) migration and split-step Fourier (SSF) migration. Moreover, we use 2D marine data from the South China Sea to show that the PSO-GSP migration can better image strong reflectors than conventional imaging methods.     

地震波在孔隙介质中低频衰减现象的粘弹性特征分析及近似（英文）

由于介观尺度的孔隙流体流动,弹性波传播过孔隙岩层时在地震频段表现出较强的频散和衰减。Johnson理论给出了在任意孔隙形状的条件下,部分气水饱和孔隙介质的理论相速度和品质因子的解析解。本文在Johnson模型的基础上,通过对Q值曲线的低频和高频近似,推导了Q值曲线的近似公式,以及基于孔隙介质基本地球物理参数和孔隙斑块几何形态参数T和比表面积S/V的最大衰减Qmin近似公式。通过与理论值的对比,对Qmin近似公式存在的线性误差进行改正,进一步提高了精度。复杂的斑块形态对最大衰减Qmin和过渡频率ftr的都产生一定影响,且对ftr影响更大。因为数值模拟直接求解介观尺度的Biot孔隙介质方程需要极大的计算量,我们使用Zener模型建立了等效粘弹模型,有效地模拟了地震频带内的衰减和频散现象。     

三维复杂介质的块状建模和试射射线追踪

为了解决三维复杂介质的射线追踪,本文改变了传统的三维层状地层的建模描述方式,提出了块状结构的建模描述方法,结合三角形面片来描述地质界面,可以构造非常复杂的三维地质模型．为了满足射线追踪的需要,本文对模型界面内的法向量进行光滑处理,光滑后的法向量在界面内是连续变化的．在块状模型的基础上,本文运用三角形的面积坐标,提出了几种试射角度的修正方法：细分三角形法、分割三角形法和子三角形法,计算表明子三角形法最好．文中给出了三维块状模型和射线追踪实例．     

Study on equivalent velocity pulse of nearfault ground motions

Introduction By analyzing earthquake motions, we could find that earthquake motions near the causativefault have two characteristics. One is the remarkable directivity effect. The amplitude of thefault-normal component is larger than that of the fault-parallel one; the other is obvious pulse mo-tions. Bertero, et al (1977) studied the earthquake records of the 1971 San Fernando earthquake.They first pointed out that some ground motions recorded near the causative fault is characterizedb…     

Seismic dynamics of offshore breakwater on liquefiable seabed foundation

Offshore structures, such as composite breakwaters, are generally vulnerable to strong seismic wave propagating through loose or medium-dense seabed foundation. However, the seismically induced failure process of offshore structures is not well understood. In this study, seismic dynamics of a composite breakwater on liquefiable seabed foundation is investigated using a fully coupled numerical model FSSI-CAS 2D. The computation results show that the numerical model is capable of capturing a variety of nonlinear interaction phenomena between the composite breakwater and its seabed foundation. The numerical investigation demonstrates a three-stage failure process of the breakwater under seismic loading. In this process, the far-field seabed can become fully liquefied first, inducing excessive settlement of the structure, followed by significant lateral movement and tilting of the structure when the near-field soil progressively liquefies. The study demonstrates great promise of using advanced numerical analysis in geotechnical earthquake design of offshore structures.     

TI介质局部角度域射线追踪与叠前深度偏移成像

研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations

In the context of wide-angle seismic profiling, the determination of the physical properties of the Earth crust, such as the elastic layer depth and seismic velocity, is often performed by inversion of P- and/or S-phases propagation data supplying the geometry of the medium (reflector depths) or any other structural parameter (P- or S-wave velocity, density...). Moreover, the inversion for velocity structure and interfaces is commonly performed using only seismic reflection travel times and/or crustal phase amplitudes in isotropic media. But it is very important to utilize more available information to constrain the non-uniqueness of the solution. In this paper, we present a simultaneous inversion method of seismic reflection travel times and polarizations data of transient elastic waves in stratified media to reconstruct not only layer depth and vertical P-wave velocity but also the anisotropy feature of the crust based on the estimation of the Thomsen’s parameters. We carry out a checking with synthetic data, comparing the inversion results obtained by anisotropic travel-time inversion to the results derived by joint inversion of seismic reflection travel times and polarizations data. The comparison proves that the first procedure leads to biased anisotropic models, while the second one fits nearly the real model. This makes the joint inversion method feasible. Finally, we investigate the geometry, P-wave velocity structure and anisotropy of the crust beneath Southeastern China by applying the proposed inversion method to previously acquired wide-angle seismic data. In this case, the anisotropy signature provides clear evidence that the Jiangshan-Shaoxing fault is the natural boundary between the Yangtze and Cathaysia blocks.     

薄覆盖层地区隐伏断层及其上断点探测的地震方法技术——以废黄河断层为例

通过在同一条测线上应用三种不同地震勘探手段(共偏移距地震反射法、横波反射法与高分辨率折射法)联合反演的方法,获得了测线控制地段内废黄河断层的确切位置、上断点埋深以及速度分布图像.探测结果表明:在薄覆盖层地区的断裂调查中,上述三种技术手段的联合反演要比单独使用其中任何一种手段更加可靠,并能从不同角度查明断层的位置、性质及其特征,为钻孔联合剖面位置的布设和钻孔深度的设计提供地震学依据.经高精度钻孔联合地质剖面证实,三种地震勘探方法反演得到的主要地层界面和构造特征都与钻孔联合地质剖面吻合较好.试验表明了上述三种地震勘探方法在基岩面埋深较浅地区联合反演的可行性以及地震勘探与钻孔联合地质剖面相结合的工作方法的有效性.     

基于牛顿迭代法求解群速度的地震波射线追踪模拟

地下介质中普遍存在着各向异性,当前基于各向异性的地震波射线追踪多是在弱各向异性介质中进行且采用群速度近似表示方法,这些近似方法在强各项异性介质中会导致很大误差而无法真正模拟地震波的传播规律。根据地下普遍存在各向异性的事实和地震波基本传播规律,提出利用牛顿迭代法高效求解群速度,基于Paraview平台自动化构建三维地质模型,采用最短路径法进行地震波射线追踪模拟及可视化,实现对复杂三维地质的速度不均匀性和各向异性的表达,为三维地质模型的构建和地震波射线追踪模拟及可视化提供一种新思路,并以华北克拉通山西断陷带北部局部区域为例进行研究。结果表明,该方法能够减少由各向异性对地震波传播模拟造成的影响,清晰表达了研究区地质结构和各向异性特点,在对复杂三维地质结构的解读中能够较好应用。     

基于地震成像数据稀疏反演的不连续及非均质地质体检测方法

本文针对油气藏储层预测中的不连续及非均质地质信息识别问题,研究基于地震成像数据的稀疏反演方法.由于该类地质体的地震响应特征为弱信号,因此利用平面波破坏滤波器由地震成像数据中去除强反射同相轴.在此基础上,对剩余的地震数据进行非线性加强滤波,并构建L1稀疏反演模型.为有效求解L1模型,采用非光滑泛函L1范数逼近和拟牛顿求解算法.该方法考虑稀疏先验信息,能够提高反演结果信噪比.缝洞模型测试验证该方法在检测断点、微断裂、散射点等小尺度地质体上的有效性,塔北缝洞型碳酸盐岩储层预测的应用效果进一步证实该方法的实用性.     

Relationship between the transition frequency of local fluid flow and the peak frequency of attenuation

Local fluid flow (LFF) at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band. LFF is easily influenced by the structure and boundary conditions of the porous media, which leads to different behaviors of the peak frequency of attenuation. The associated transition frequency can provide detailed information about the trend of LFF; therefore, research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation (i.e., inverse of quality factor) facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media. In this study, we firstly obtain the transition frequency of fluid flux based on Biot’s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell (RUC). We then analyze changes of these two frequencies in porous media with different porous properties. Finally, we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs. This setup can facilitate the understanding of the effect from the undrained boundary condition. Results demonstrate that these two frequencies have the same trend at low water saturation, but amplitude variations differ between the frequencies as the amount of saturation increases. However, for cases of high water saturation, both the trend and the amplitude variation of these two frequencies fit well with each other.