高速非均质体多尺度散射的实验研究（英文）

地下介质内部存在着复杂的非均质体,其非均质尺度对地震波传播具有显著影响。本文设计了一组11种直径的玻璃珠高速非均质模型,研究了这种三维非均质体对纵波传播的影响。每个模型中玻璃珠的体积分数（10%）相同,玻璃珠直径在0.18-11mm范围内变化。通过三种基质频率（0.83、0.61和0.34MHz）的超声脉冲测试,分析了纵波随非均质体尺度变化的波场特征及基质波数与非均质尺度乘积（ka）对纵波速度的影响。实验发现,非均质尺度影响波的传播,其波场变化较为复杂,记录的波形、振幅、速度等的变化与非均质尺度有明显的相关性;对于非均质体尺度较大的前向散射,在实验室中观察到明显的直达波与衍射波,这表明对于大尺度非均匀体,不能忽略直达波与衍射波的影响;速度与ka的关系表现出频率依赖性,这主要是由波数k、非均质尺度a引起的速度变化率不同引起的;不论是依据波形、振幅变化规律,还是不同频率速度与非均质尺度a之间的关系,确定的射线近似转换点在ka=10附近;当ka小于1时,速度缓慢变化并逐渐接近有效介质速度;对于速度扰动较大的非均匀介质,实验室测量的射线速度要明显小于微扰理论预测的速度。     

随机-均匀介质散射波特征相干分析研究

基于统计学方法建立了随机-均匀介质模型,利用有限差分法进行正演模拟,根据地震剖面中散射波同相轴随着模型横向自相关长度的增大其连续性逐渐变好的特点,通过相干技术研究了散射波同相轴的连续性.研究结果表明,当模型横向自相关长度与地震波长比值较小（a/λ<1）时,散射波同相轴连续性较差,其相干值变化不明显;当模型横向自相关长度与地震波长比值较大（a/λ>1）时,散射波同相轴连续性明显变好,其相干值明显增大.     

不同纵波速度介质中炸药震源激发数值模拟研究

震源激发是影响地震勘探效果的一个重要因素,在油气勘探中需要选择合适的震源参数以获取较好的地震记录。本文以数值模拟方法为基础,通过使用ANSYS/LS-DYNA软件,对单井条件下不同岩土介质中柱状炸药震源激发进行了对比研究,计算得到相应的地震子波。通过对比分析各组参数对应的地震子波振幅和频谱,选择合适的激发参数。研究结果表明,在纵波速度较低的介质中,中低爆速、细长药柱炸药激发产生的地震子波频率高、能量大,而在纵波速度较高的介质中,高爆速、短粗药柱炸药激发产生的地震子波频率高、能量大。该研究结果与实际勘探结果比较符合,可为实际勘探中激发参数的选择提供参考。      

井间地震方向性震源模拟方法

为了能够在水平方向上传播更远的距离,井间地震震源辐射能量具有水平方向强的特性,这是井间地震资料大入射角反射波能量强的主要原因.为了研究井间地震震源方向特性对地震信号的影响,本文提出了圆形震源阵列模拟方法.基于声波方程的一阶应力一速度方程的交错网格有限差分算法,进行了均匀各向同性无限介质的直达波和反射波正演数值模拟.数值模拟结果表明,井间方向性震源辐射能量的方向特性影响了地震记录的能量分布,为进一步研究方向性震源的能量补偿提供了理论依据.     

各向异性介质地震波场数值模拟及特征分析

由于各向异性广泛存在于地下岩石中,随着勘探精度的不断提高,对地下介质的各向同性假设越来越不能够满足于现状,因此对各向异性介质的数值模拟显得更为重要。本文推导了各向异性介质的弹性波动方程,总结了震源类型,通过PML方法处理了人工边界问题,通过快照分析验证了数值频散、稳定性条件。研究结果表明:① PML完全匹配层,可较好地解决人工边界问题;②减小空间采样间隔压制数值频散比减小时间采样间隔效果要好得多,盲目减小时间采样间隔会大大降低数值模拟的运算效率;③各向异性介质中弹性波场中除含有准纵波外,还含有速度较慢的准横波;④准纵波波前能量要比由各向异性引起的准横波能量强,准纵波和准横波的波前随着各向异性介质参数的变化而变化。      

井间高分辨率纵横波层析成像研究井间油藏

采用视速度-偏振井间波场分离法分离得到纵波和横波波场,使用跨孔地震走时层析成像级联方法计算井间纵波和横波速度结构以及泊松比空间分布.综合纵横波速度、泊松比和测井参数,分析两井间岩性和储层特性.数值模型试验表明该方法垂直分辨率达到主频波长的1/4,并能准确分辨倾斜断层和垂直断层.处理相同纵波资料结果表明本文方法分辨率远远高于以往的井间走时层析成像方法,比井间声波层析成像分辨率更高.使用该方法处理垦71区井间资料,得到同台同源纵波和横波速度结构和泊松比分布.泊松比的低/高变化与测井参数指示的砂/泥岩层基本吻合.根据泊松比分布,参考岩石地震物理学测试数据,可以区分井间介质中的泥岩层、砂岩层,以及砂岩层中的饱水和饱油区,确定储层的连通性,圈定井间的剩余油藏.     

地震散射波的高精度数值模拟与振幅分析

随着地震勘探要求和技术的不断提高,散射波地震勘探方法逐渐引起人们的重视,尤其在构造复杂地区的高精度勘探中,散射波理论具有广泛的应用前景.本文在讨论了现阶段几种散射波数值模拟方法的利弊后,设计了一种基于微扰论的FK域积分法,在散射场的二次震源和空间能量衰减处理两方面做了改进,其计算精度和效率均得到了提高.此外,采用该算法针对散射波动力学进行了正演模拟试算,从扰动介质空间分布和扰动量两个角度定量地分析了散射波振幅动力学方面的特征.在此基础上得到了五个经验常数,并给出了对应的计算公式和物理意义,通过这五个常数可以划分散射波不同形态的动力学特征.     

超声实验中谱比法衰减的散射与本征吸收特性

实际地球介质的弹性波衰减,无论是对于地球内部结构反演、含油气储层描述,还是对于地震岩石物理本身的基础研究而言,都十分重要.通过岩石物理实验测量岩石的弹性波衰减是相关研究的基本手段,而谱比法是一种最常用的测量方法.利用谱比法测量的弹性波衰减,简称谱比法衰减,由介质的本征吸收特性和非均质散射特性共同决定.然而,在谱比法理论中,没有将这两种衰减特性的作用区分对待.因此,实验室测量的谱比法衰减中是否含有散射特性,其散射特性是否可以忽略,常不明确.本文基于对岩石物理实验中样品波形特征的基本认识,通过能流理论描述了弹性波的本征吸收和非均质散射过程,以及对应的尾波形成机理.文章还应用能流理论测量岩石的弹性波衰减,即能流法衰减.无论是通过理论推导,还是通过实验验证,谱比法衰减与能流法衰减都显示出良好的一致性.这两种直达波衰减的一致性暗示了谱比法衰减中存在着两种衰减特性.文中还对比分析了谱比法衰减与喷射流本征衰减的实验结果,该结果同样展示了谱比法衰减中可能存在着散射特性.因此,通过能流理论,可以利用样品波形的尾波特征定性判断谱比法衰减中可能存在的散射特性.即异常明显的尾波振幅可能意味着强烈的直达波非均质散射,谱比法衰减中的散射特性因而不可忽略.     

运用分形方法模拟地震勘探中的三维随机介质

实际的地震勘探遇到的是复杂随机介质,常规层状介质地震波场理论已经不能适应于高精度地震工作,需要进一步发展非均匀介质地震波场理论.本文首次给出了运用分形方法建立三维随机介质的具体算法;其次引入了三维锥形函数来压制建模时产生的离散误差;最后,通过建立几种三维随机介质模型来讨论尺度系数和豪斯特指数对模拟结果的影响,并且简单构建了区域多尺度随机介质模型.模拟的结果表明,分形方法可以更方便、灵活地模拟三维随机介质.针对南海某地区生物礁储层非均质性强的特点,在建立生物礁速度模型时,对目标区充填了多尺度的分形随机介质,波动方程正演模拟结果表明,分形随机介质可以用来描述生物礁等非均质性储层.     

井间定向震源正演模拟与特性分析

为减少管波对采集数据的影响,井间地震采用了定向震源,其能量辐射特性影响接收波场的能量分布与后续成像效果.研究定向震源的正演模拟并分析能量辐射特性,为井间地震成像过程中能量补偿提供基础数据.采用圆形震源阵列法,在近源区内按照质点与震源点连线的方向设置振动函数强度,水平方向能量最强,垂直方向能量最弱,模拟的辐射能量具有随方向改变的特点.实际井间地震能量方向性非常强,因而同时采用吸收衰减层法进一步减弱纵向的地震波能量.吸收衰减层法在震源周围设置具有衰减因子的海绵层,通过不同方向设定不同的衰减因子来控制震源的方向特性.衰减因子采用了渐变函数,保证海绵层边界不产生人为反射干扰.将圆形震源阵列法与吸收衰减层法联合进行双层介质模型和实际速度场模拟,得到的炮记录方向性更加明显.两种方法的有机结合,能够更好地模拟井间定向震源.在圆形震源阵列法可行性验证基础上,分析了阵列半径对接收波场能量的影响,初步确定了模拟定向震源的最佳阵列半径.通过分析方向性参数和井间距对接收波场能量分布的影响,获得了定向震源激发下接收波场的能量分布,为在井间地震偏移成像中进行能量补偿建立了基础.     

EXPERIMENTAL STUDY ON THE CHANGES OF ULTRASONIC CODA WAVE AND ACOUSTIC EMISSION DURING ROCK LOADING AND DEFORMATION

The coda wave propagation path has received extensive attention as it is more sensitive to small changes in the medium than the direct wave. During the process of loading, the wave velocity, medium or source changes may cause the coda wave to change. The physical mechanism of change in the ultrasonic coda wave varies during different deformation stages. Meanwhile, there exist local damages in the rock sample during the deformation, and it will be accompanied by acoustic emission. Combining the ultrasonic coda wave and acoustic emission is beneficial to characterize the coda wave characteristics and damage degree of the sample at different deformation stages. In this paper, three kinds of rocks, including granodiorite, marble and sandstone with the sizes of 50mm×50mm×150mm, are used to carry out observations of ultrasonic coda wave and acoustic emission during the whole process of loading so as to study characteristics of the coda wave at different deformation stages. The major results are given below: 1)There is a good correspondence between the coda wave variation and the acoustic emission evolution process. When the acoustic emission frequency increases, the coda wave changes accordingly. In particular, the coda wave changes in the early stages of increased acoustic emission frequency, which indicates that the early damage information of rock can be obtained by analysis of the coda wave. 2)The physical mechanism of the coda wave change is different in different deformation stages. At the initial stage of loading, there are obvious scatterer changes in the coda wave change; then, in the linear elastic deformation stage, the wave velocity change is dominating; in the late-stage of loading, the scatterer change increases and coexists with the wave velocity change, the scatterer change effect is related with the rock micro-fracture degree, the rock will locally be damaged before rupturing, and the role of the scatterer will be enhanced. 3)With the increase of loading, the amplitude of increase of the wave velocity generally decreases gradually, which is basically consistent with the understanding obtained through the direct wave. The interference of acoustic emission can be eliminated because of the Kaiser effect when analyzing the coda wave. The consistency of the wave velocity change and stress loading and unloading is further verified. 4)The micro-fracture generated during rock deformation will change the physical mechanism of the coda wave change, and the scatterer effect will be significantly enhanced. At the same time, the acoustic emission waveform will cause interference to the ultrasonic coda wave. This means that attention needs to be paid when analyzing rock damage using only coda wave data. In short, the ultrasonic coda wave and acoustic emission can reflect the damage inside the rock, and the change mechanism of the coda wave in different deformation stages is different. The joint observation of the two can play a mutual verification role, which is conducive to improving the reliability of the observation results.     

Digital core based transmitted ultrasonic wave simulation and velocity accuracy analysis

Transmitted ultrasonic wave simulation (TUWS) in a digital core is one of the important elements of digital rock physics and is used to study wave propagation in porous cores and calculate equivalent velocity. When simulating wave propagates in a 3D digital core, two additional layers are attached to its two surfaces vertical to the wave-direction and one planar wave source and two receiver-arrays are properly installed. After source excitation, the two receivers then record incident and transmitted waves of the digital rock. Wave propagating velocity, which is the velocity of the digital core, is computed by the picked peak-time difference between the two recorded waves. To evaluate the accuracy of TUWS, a digital core is fully saturated with gas, oil, and water to calculate the corresponding velocities. The velocities increase with decreasing wave frequencies in the simulation frequency band, and this is considered to be the result of scattering. When the pore fluids are varied from gas to oil and finally to water, the velocity-variation characteristics between the different frequencies are similar, thereby approximately following the variation law of velocities obtained from linear elastic statics simulation (LESS), although their absolute values are different. However, LESS has been widely used. The results of this paper show that the transmission ultrasonic simulation has high relative precision.     

中层大气中重力波饱和机制的数值分析

通过对向上传播的重力波波包在中层大气中的非线性传播过程进行数值模拟,讨论中层大气中重力波的饱和机制.数值模拟结果表明,向上传播的重力波波包的扰动振幅在接近波包的本征水平相速度之前随高度单调增长,而当波振幅接近本征水平相速度时,在对流不稳定区域出现等位温面的翻转,同时波振幅的增长达到饱和（波振幅随高度不发生变化）.小尺度对流在等位温面的翻转和波饱和发生后产生,随后波包开始破碎,这些非线性过程的最终结果产生湍流.表明导致重力波饱和的关键因素是等位温面的翻转而不是诸如波破碎、湍流、波-流相互作用等其他的一些物理过程.     

Numerical simulation of attenuation and group velocity of guided ultrasonic wave in grouted rock bolts

In this paper, the guided ultrasonic wave propagating in grouted rock bolts was simulated with finite element method. An 800 mm partially grouted cylindrical rock bolt model was created. Dynamic input signals with frequency from 25 to 100 kHz were used to excite ultrasonic wave. The simulated waveform, group velocity and amplitude ratio matched well with the experimental results. This model made it possible to study the behaviour of the guided waves in the grouted bolt along its central axis. Analysis of the simulated results showed that the group velocity in grouted rock bolts is constant along the grouted length, and the boundary effect on the group velocity is negligible. This paper also presents methods to determine the attenuation coefficient from simulation and to determine the boundary effect on attenuation at the bolt ends. The analysis showed that the attenuation of the guided wave propagating inside the grouted bolts is similar to the theoretical solution in steel bar with infinite length. After correction for the boundary effects the grout length of a grouted rock bolt can be determined using the measured attenuation, with sufficient accuracy.     

Boundary-reflected waves and ultrasonic coda waves in rock physics experiments

Ultrasonic coda waves are widely used to study high-frequency scattering. However, ultrasonic coda waves are strongly affected by interference from by boundary-reflected waves. To understand the effect of boundary-reflected waves, we performed ultrasonic experiments using aluminum and shale samples, and the rotating staggered-mesh finite-difference method to simulate the wavefield. We analyzed the wavefield characteristics at the different receiving points and the interference characteristics of the boundary-reflected waves with the ultrasonic coda wave, and the effect of sample geometry on the ultrasonic coda waves. The increase in the aspect ratio of the samples delays the interference effect of the laterally reflected waves and reduces the effect on the ultrasonic coda waves. The main waves interfering with the ultrasonic coda waves are laterally reflected PP-, PS-, PPP-, and PPS-waves. The scattering and attenuation of the high-frequency energy in actual rocks can weaken the interference of laterally reflected waves with the ultrasonic coda waves.     

含混合裂隙、孔隙介质的纵波衰减规律研究

地下多孔介质中的孔隙类型复杂多样,既有硬孔又有扁平的软孔.针对复杂孔隙介质,假设多孔介质中同时含有球型硬孔和两种不同产状的裂隙(硬币型、尖灭型裂隙),当孔隙介质承载载荷时,考虑两种不同类型的裂隙对于孔隙流体压力的影响,建立起Biot理论框架下饱和流体情况含混合裂隙、孔隙介质的弹性波动方程,并进一步求取了饱和流体情况下仅由裂隙引起流体流动时的含混合裂隙、孔隙介质的体积模量和剪切模量,随后,在此基础上讨论了含混合裂隙、孔隙介质在封闭条件下地震波衰减和频散的高低频极限表达式.最后计算了给定模型的地震波衰减和频散,发现地震波衰减曲线呈现"多峰"现象,速度曲线为"多频段"频散.针对该模型分析讨论了渗透率参数、裂隙纵横比参数以及流体黏滞性参数对于地震波衰减和频散的影响,表明三个参数均为频率控制参数.     

基于含气泡液体声波方程的海底冷泉数值模拟

海底冷泉羽状流与海底天然气水合物的分布密切相关,对水合物稳定带的边界具有指示作用,是未来能源勘探的重要领域.研究海底冷泉羽状流的地震响应特征,对确定天然气水合物的储集区域及成藏环境等均有重要意义.当前获得海底冷泉羽状流的地震响应主要通过数值模拟进行,然而该过程所依据的含气泡介质声速模型及随机介质理论不能完整描述海底冷泉的物理性质,采用的声波方程也不适用于高频地震波数值模拟.为了准确地实现海底冷泉羽状流地震波数值模拟,精确分析其地震响应特征,本文提出利用Keller-Miksis气泡振动模型来描述气泡在声波作用下的运动状态,同时考虑气泡间的相互作用,建立海底冷泉气泡模型.在此基础上,本文创新性地采用含气泡液体声波方程进行海底冷泉高频地震波数值模拟.数值模拟结果表明,本文提出的方法能够实现海底冷泉羽状流地震响应的高精度数值模拟.     

部分饱和孔隙岩石中声波传播数值研究

利用基于Biot理论的孔隙弹性介质的高阶交错网格有限差分算法,模拟了具有随机分布特征的多种流体饱和岩石中声波在中心频率分别为25,50,75,100kHz时的声场特点. 对于一个由两种成分（气和水）饱和的岩石模型, 假设含不同流体的孔隙介质随机分布在不同的宏观区域,该区域尺度远小于计算的声波波长;组成模型的两种随机分布介质具有相同的固体骨架参数、渗透率和孔隙度,但分别被具有不同压缩性、密度和黏滞系数特性的水和气饱和. 计算和统计分析结果表明,在两种孔隙成分随机分布的部分饱和条件下纵波速度比较复杂,除骨架参数外,其变化主要依赖于中心频率、各种孔隙成分饱和度及饱和介质的速度. 比较该随机分布模型、Gassmann理论模型和White的“气包”模型,发现三种模型得到的纵波速度和衰减规律有较好的定性对应关系. 其次,按照这种随机计算模型的处理方法,本文还首次计算了一个三种流体成分充填饱和的例子,即岩石模型中的孔隙被水、油和气部分饱和,计算时保持模型含水饱和度不变而只改变含油和含气饱和度. 在这种计算条件下,纵波速度随中心频率呈增大的趋势但有起伏变化. 声场快照显示了各种转换波在多种孔隙成分充填（两种和三种孔隙成分）岩石中的声场特征,复杂的水-油-气界面的非均匀分布对声场有重要影响,纵波能量主要转换形成了较为复杂的多种慢纵波和横波.     

套管井偶极弯曲波频散向高频偏移的特性

对套管井偶极弯曲模式波的频散特性进行了系统的数值考察、实例对比和分析.发现套管井弯曲波频散曲线随地层横波速度的降低,特别是地层横波速度小于2000 m·s^-1以下,会迅速移向高频区,偶极弯曲波基础模式主频散区（或截止频率）可出现在13 kHz以上,以致超出了现行低频偶极子声波测井仪的激发与接收频带,这是一过去没有被研究者注意到的现象,并进一步被现场实例所证实.研究表明控制套管井弯曲波频散曲线主频散区位置的主要是钢套管的厚度和地层横波速度.对地层横波速度大于井孔流体声速的快速地层,在钢套管壁厚一定（8 mm）的情况下,频散曲线主频散区可移至11 kHz以上,可能出现的最大可能频域位置是同一井孔内径,井外全钢时的频散曲线上等于、小于地层横波速度那一段,这对各种地层和套管参数都是适用的.对地层横波速度小于等于井孔流体声速（1500 m·s^-1）的慢速地层,弯曲波频散曲线随地层横波速度的降低移向高频区的特点更为明显,可能移至16 kHz以上;而套管厚度的影响,也比快速地层大的多,对地层横波速度小于1380 m·s^-1的慢速地层,无论用多高的频率激发,都不能在现行使用的各类套管井（壁厚6~12 mm）中用偶极声波测井仪测到弯曲模式波.     

The influence of mesoscopic flow on the P-wave attenuation and dispersion in a porous media permeated by aligned fractures

In sedimentary rocks attenuation/dispersion is dominated by fluid-rock interactions. Wave-induced fluid flow in the pores causes energy loss through several mechanisms, and as a result attenuation is strongly frequency dependent. However, the fluid motion process governing the frequency dependent attenuation and velocity remains unclear. We propose a new approach to obtain the analytical expressions of pore pressure, relative fluxes distribution and frame displacement within the double-layer porous media based on quasi-static poroelastic theory. The dispersion equation for a P-wave propagating in a porous medium permeated by aligned fractures is given by considering fractures as thin and highly compliant layers. The influence of mesoscopic fluid flow on phase velocity dispersion and attenuation is discussed under the condition of varying fracture weakness. In this model conversion of the compression wave energy into Biot slow wave diffusion at the facture surface can result in apparent attenuation and dispersion within the usual seismic frequency band. The magnitude of velocity dispersion and attenuation of P-wave increases with increasing fracture weakness, and the relaxation peak and maximum attenuation shift towards lower frequency. Because of its periodic structure, the fractured porous media can be considered as a phononic crystal with several pass and stop bands in the high frequency band. Therefore, the velocity and attenuation of the P-wave show an oscillatory behavior with increasing frequency when resonance occurs. The evolutions of the pore pressure and the relative fluxes as a function of frequency are presented, giving more physical insight into the behavior of P-wave velocity dispersion and the attenuation of fractured porous medium due to the wave-induced mesoscopic flow. We show that the specific behavior of attenuation as function of frequency is mainly controlled by the energy dissipated per wave cycle in the background layer.