Anisotropic viscoelastic models with singular memory

The physical background of singular memory models and in particular the Cole–Cole model is discussed. Three models of anisotropic linear viscoelasticity with frequency-dependent stiffness coefficients are considered. The models are constructed in such a way that anisotropic properties are separated from anelastic effects. Two of these models represent finite-speed wave propagation with singularities at the wavefronts (the exponential relaxation model) and without singularities at the wavefronts (the Cole–Cole model), while a third model called the fractional model is related to the constant Q with unbounded propagation speed. The Cole–Cole and fractional models belong to the class of singular memory models studied earlier because of their applications in polymer rheology, poroelasticity, poroacoustics, seismic wave propagation and other applications. Well-posedness of initial boundary value problems with mixed Dirichlet–Neumann boundary conditions is established for the three models. Regularity properties of the three models are examined.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

粘弹各向异性介质中地震波场模拟与特征

通过引入记忆变量,可以避免粘弹性应力-应变关系中的褶积运算,使波场数值模拟易于实现.通过伪谱法对粘弹各向异性介质中的qP波、qS波数值模拟,结合理论分析,研究了粘弹各向异性介质中速度各向异性和衰减各向异性.衰减各向异性要比速度各向异性更为显著,并且qS波比qP波的衰减各向异性明显.粘弹各向异性介质中,粘弹性对波的影响主要在于波的衰减,各向异性主要影响波前面形状.     

Source-independent waveform inversion for attenuation estimation in anisotropic media

In previous publications, we presented a waveform-inversion algorithm for attenuation analysis in heterogeneous anisotropic media. However, waveform inversion requires an accurate estimate of the source wavelet, which is often difficult to obtain from field data. To address this problem, here we adopt a source-independent waveform-inversion algorithm that obviates the need for joint estimation of the source signal and attenuation coefficients. The key operations in that algorithm are the convolutions (1) of the observed wavefield with a reference trace from the modelled data and (2) of the modelled wavefield with a reference trace from the observed data. The influence of the source signature on attenuation estimation is mitigated by defining the objective function as the ℓ₂-norm of the difference between the two convolved data sets. The inversion gradients for the medium parameters are similar to those for conventional waveform-inversion techniques, with the exception of the adjoint sources computed by convolution and cross-correlation operations. To make the source-independent inversion methodology more stable in the presence of velocity errors, we combine it with the local-similarity technique. The proposed algorithm is validated using transmission tests for a homogeneous transversely isotropic model with a vertical symmetry axis that contains a Gaussian anomaly in the shear-wave vertical attenuation coefficient. Then the method is applied to the inversion of reflection data for a modified transversely isotropic model from Hess. It should be noted that due to the increased nonlinearity of the inverse problem, the source-independent algorithm requires a more accurate initial model to obtain inversion results comparable to those produced by conventional waveform inversion with the actual wavelet.     

非弹性层状介质地震波频变AVO响应模拟及分析

以非弹性层状介质为模型,基于广义传播矩阵理论计算地震波频变反射系数,算法中同时考虑了与频散和衰减有关的地层岩性因素,以及与薄互层有关的地层结构因素.实现了岩石物理模型、反射系数这两个计算过程的"无缝"连接,精确考虑了由复数弹性模量表示的地层非弹性因素,也为在同一反射模型中考虑源于不同物理机制的频散与衰减提供了方法.数值模拟结果验证了算法的有效性和稳定性,计算结果表明,非弹性薄层的反射振幅随频率先增加后减小,不一定表现常规"低频亮点"异常;同时,薄互层条件下的频散与衰减使得地震反射波的频谱以及AVO特征呈现与频率相关的复杂变化.本文完善了频变AVO算法,为含油气储层频变AVO响应的模拟和分析提供了方法.     

线性粘弹介质中地震波场数值模拟

本文将以往两种粘弹介质中地震波模拟方法的优点结合起来，以模型理论和积分本构方程为基础，从理论上分析了模型对地震波场的影响；采用交错网格有限差分法对粘弹介质中的地震波进行数值模拟.数值计算结果表明该方法不仅便于计算，同时也便于从力学的角度来分析地震波的传播.数值计算结果与理论分析一致，说明这种方法可以更为有效地模拟粘弹介质中地震波的传播.     

Effective bulk modulus of the pore fluid at patchy saturation

We explore a package of parallel porous layers, each filled with a different fluid. Assume that this package is sampled by an elastic wave with the wavelength much larger than the thickness of an individual layer. Also assume that the layers are hydraulically isolated from each other, meaning that the diffusion length is smaller than that of the individual layer. This assumption is relevant to a patchy saturation scenario. Suppose that we wish to conduct the fluid substitution operation on this package treated as a single porous elastic body. What is the effective bulk modulus of the pore fluid to be used in this operation that will result in the same elastic modulus as computed by Backus averaging the individual moduli of the layers? We address this question analytically by assuming that the porosity, dry frame, and the mineral matrix properties of the individual layers are the same for all layers. The only difference between the layers is the pore fluid. We find that the resulting effective bulk modulus of the fluid thus derived falls between the arithmetic and harmonic averages of the fluid bulk moduli in the layers. It can be approximated by a linear combination of these two bounds where the weights are 0.50 and 0.50 or 0.75 for the arithmetic average and 0.25 for the harmonic average, depending on the elastic moduli of the dry frame, the mineral, and the pore fluids. This solution also provides a relation between the effective bulk modulus of the pore fluid in the system under examination and water saturation to be used in the fluid substitution operation at a coarse spatial scale.     

Seismic wave attenuation in fluid-saturated porous media

Contrary to the traditional view, seismic attenuation in Biot's theory of fluid-saturated porous media is due to viscous damping of local (not global) pore-fluid motion. Since substantial inhomogeneities in fluid permeability of porous geological materials are to be expected, the regions of highest local permeability contribute most to the wave energy dissipation while those of lowest permeability dominate the fluid flow rate if they are uniformly distributed. This dichotomy can explain some of the observed discrepancies between computed and measured attenuation of compressional and shear waves in porous earth. One unfortunate consequence of this result is the fact that measured seismic wave attenuation in fluid-filled geological materials cannot be used directly as a diagnostic of the global fluid-flow permeability.     

黏弹性介质瑞雷波有限差分模拟与特性分析

本文通过数值模拟研究了介质黏弹性对瑞雷波传播的影响.模拟采用结合了交错Adams-Bashforth时间积分法、应力镜像法和多轴完美匹配层的标准交错网格高阶有限差分方案.通过模拟结果和理论结果对比,测试了方法的精度,验证了结果的正确性.在均匀半空间模型中,分别从波场快照、波形曲线及频散能量图三个角度,对黏弹性介质瑞雷波衰减和频散特性进行了详细分析.两层速度递增模型被用于进一步分析瑞雷波在黏弹性层状介质中的特性.结果表明：由于介质的黏弹性,瑞雷波振幅发生衰减,高频成分比低频成分衰减更剧烈,衰减程度随偏移距增大而增强;瑞雷波相速度发生频散,且随频率增大而增大,频散能量的分辨率有所降低;黏弹性波动方程中的参考频率,不会影响瑞雷波振幅衰减和相速度频散的程度,但决定了黏弹性和弹性介质瑞雷波相速度相等的频率位置.本研究有助于人们更好地理解地球介质中瑞雷波的行为,并为瑞雷波勘探的应用和研究提供了科学和有价值的参考.

     

基于地震正交各向异性的裂缝介质岩石力学建模研究

岩石力学模型是描述地层原地应力状态的基础,而常规各向同性模型难以刻画地层本征横观各向同性(VTI)和天然高角度裂缝的耦合作用,建立更为准确的正交各向异性模型变得尤为重要.本文利用VTI介质中发育单组垂直缝的Schoenberg裂缝等效模型,简化了一般正交各向异性介质的表征参数,推导了由背景VTI介质弹性参数和裂缝参数的各向异性杨氏模量、泊松比和水平地应力的精确方程.借鉴Thomsen弱各向异性近似思路,舍去裂缝弱度参数高阶扰动项,推导了正交各向异性介质岩石力学近似方程,由裂缝弱度表示的近似方程更为直观地反映了裂缝对背景介质力学性质的影响,杨氏模量和泊松比值随裂缝法向弱度的增大而显著减小.选取页岩实验数据进行数值实验,结果表明本文方程与实际物理规律相吻合,随裂缝弱度的增加,形成相同应变所需的水平地应力值降低,且近似方程与模拟结果相对误差小于5%,有助于提高岩石力学参数估算和地应力预测精度.     

Generalized anisotropy parameters and approximations of attenuations and velocities in viscoelastic media of arbitrary anisotropy type – theoretical and experimental aspects*

Seismic anisotropy in geological media is now widely accepted. Parametrizations and explicit approximations for the velocities in such media, considered as purely elastic and moderately anisotropic, are now standards and have even been extended to arbitrary types of anisotropy. In the case of attenuating media, some authors have also recently published different parametrizations and velocity and attenuation approximations in viscoelastic anisotropic media of particular symmetry type (e.g., transversely isotropic or orthorhombic). This paper extends such work to media of arbitrary anisotropy type, that is to say to triclinic media. In the case of homogeneous waves and using the so‐called ‘correspondence principle’, it is shown that the viscoelastic equations (for the phase velocities, phase slownesses, moduli, wavenumbers, etc.) are formally identical to the corresponding purely elastic equations available in the literature provided that all the corresponding quantities are complex (except the unit vector in the propagation direction that remains real). In contrast to previous work, the new parametrization uses complex anisotropy parameters and constitutes a simple extension to viscoelastic media of previous work dealing with non‐attenuating elastic media of arbitrary anisotropy type. We make the link between these new complex anisotropy parameters and measurable parameters, as well as with previously published anisotropy parameters, demonstrating the usefulness of the new parametrization. We compute the explicit complete directional dependence of the exact and of the approximate (first and higher‐order perturbation) complex phase velocities of the three body waves (qP, qS1 and qS2). The exact equations are successfully compared with the ultrasonic phase velocities and phase attenuations of the three body waves measured in a strongly attenuating water‐saturated sample of Vosges sandstone exhibiting moderate velocity anisotropy but very strong attenuation anisotropy. The approximate formulas are checked on experimental data. Compared to the exact solutions, the errors observed on the first‐order approximate velocities are small (<1%) for qP‐waves and moderate (<10%) for qS‐waves. The corresponding errors on the quality factor Q are moderate (<6%) for qP‐waves but critically large (up to 160%) for the qS‐waves. The use of higher‐order approximations substantially improves the accuracy, for instance typical maximum relative errors do not exceed 0.06% on all the velocities and 0.6% on all the quality factors Q, for third‐order approximations. All the results obtained on other rock samples confirm the results obtained on this rock. The simplicity of the derivations and the generality of the results are striking and particularly convenient for practical applications.     

Scattering and attenuation of elastic waves in random media

In seismic exploration, elastic waves are sent to investigate subsurface geology. However, the transmission and interpretation of the elastic wave propagation is complicated by various factors. One major reason is that the earth can be a very complex medium. Nevertheless, in this paper, we model some terrestrial material as an elastic medium consisting of randomly distributed inclusions with a considerable concentration. The waves incident on such an inhomogeneous medium undergo multiple scattering due to the presence of inclusions. Consequently, the wave energy is redistributed thereby reducing the amplitude of the coherent wave.The coherent or average wave is assumed to be propagating in a homogeneous continuum characterized by a bulk complex wavenumber. This wavenumber depends on the frequency of the probing waves; and on the physical properties and the concentration of discrete scatterers, causing the effective medium to be dispersive. With the help of multiple scattering theory, we are able to analytically predict the attenuation of the transmitted wave intensity as well as the dispersion of the phase velocity. These two sets of data are valuable to the study of the inverse scattering problems in seismology. Some numerical results are presented and also compared, if possible, with experimental measurements.     

Physical modeling observations of wave velocity and apparent attenuation in isotropic and anisotropic two-phase media

We study wave propagation through isotropic and anisotropic scatterer distributions in order to observe azimuthal variations in velocity and apparent attenuation. Using thin aluminum plates as physical models, we obtained seismograms for compressional and shear wave propagation through heterogeneous media. Three random distributions of scatterers are studied: circular scatterers in isotropic distributions (modeling circular scatterers), elongated scatterers in isotropic distributions (modeling randomly oriented elliptical scatterers), and elongated scatterers in anisotropic distributions (modeling aligned elliptical scatterers). All scatterers had approximately the same cross-sectional area and were filled with epoxy in order to reduce the impedance contrast. In addition to seismograms recorded for no scatterers, seismograms were recorded for several scatterer volume fractions. Azimuths were measured relative to the long axis of the aligned elongated scatterers. Velocities were calculated using travel times and phase shifts at low frequencies. The velocities measured from the data were compared to simple low-frequency average-velocity theories based on thin lamellae or on distributions of penny-shaped cracks. The apparent attenuation for different scatterer distributions was computed using spectral ratios.Comparisons of the results for circular and randomly oriented elongated scatterers were made to determine the effects of scatterer shape. As expected, the circular and randomly oriented elongated scatterers showed no systematic azimuthal variation in velocity. The velocity anomalies were systematically larger for the randomly oriented elongated scatterers than for the circular scatterers. Both methods of theoretical estimation for the isotropic velocities produced velocities significantly larger than those measured. The spectral ratios showed more apparent attenuation for the randomly oriented elongated scatterers than for the circular scatterers.Comparisons of the results for the randomly oriented and aligned elongated scatterers were made to determine the effects of anisotropy in the scatterer distribution. Compressional waves for the aligned elongated scatterers with wave propagation parallel to the scatterers had larger velocities than for the aligned elongated scatterers with wave propagation perpendicular to the scatterers for all velocity calculations. Shear wave velocities were complicated by an anomalous phase change in the shear wave seismograms for azimuths less than 40° and were not as conclusive. The general trend of the theoretical velocities is similar to the velocities calculated from the data. There are, however, what appear to be significant differences. The spectral ratios showed more apparent attenuation for the randomly oriented elongated scatterers than for the aligned elongated scatterers with wave propagation parallel to the scatterers, and less attenuation than for the aligned elongated scatterers with wave propagation perpendicular to the scatterers.     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明：本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.

     

黏弹性双相介质参数反演的同伦方法

本文以介质的动态时域响应（位移、速度、加速度）作为反演的依据,对黏弹性双相介质的材料参数进行了反分析研究.根据附加条件,并考虑到由于测量误差和测量数据的不充分造成的反演问题的不适定性,将黏弹性双相介质的参数反演问题转化为非线性泛函的极小值问题,然后应用大范围收敛的同伦方法求出使泛函极小的根作为反问题的解.最后,以二维半空间黏弹性双相介质模型为例进行了数值反演分析,数值结果表明本文方法在参数的反演过程中目标函数稳定收敛,且具有一定的抗“噪声”干扰性能;与利用速度和加速度信息反演的结果相比,利用位移信息反演的结果具有较高的精度和抗噪声能力.     

P‐ and S‐wave anisotropy to characterise and quantify damage in media: laboratory experiment using synthetic sample with aligned microcracks

Damage characterisation in solid media is studied in this work through ultrasonic measurements. A synthetic three‐dimensional printed sample including a system of horizontally aligned microcracks is used. In contrast to other manual fabrication methods presented in the literature, the construction process considered here ensures a better control and accuracy of size, shape, and spatial distribution of the microcrack network in the synthetic sample. The acoustic measurements were conducted through a specific device using triple acoustic sensors, which allows capturing at each incident direction three wave modes. The evolution of the ultrasonic velocities with respect to incident angle accounted for the damage‐induced anisotropy. The experimental results are then compared with some well‐known effective media theories in order to discuss their potential use for the following studies. Finally, we highlighted and compared the accuracy of these theories used for inversion procedure to quantify damage in the medium.     

Modelling of viscoacoustic wave propagation in transversely isotropic media using decoupled fractional Laplacians

A new wave equation is derived for modelling viscoacoustic wave propagation in transversely isotropic media under acoustic transverse isotropy approximation. The formulas expressed by fractional Laplacian operators can well model the constant-Q (i.e. frequency-independent quality factor) attenuation, anisotropic attenuation, decoupled amplitude loss and velocity dispersion behaviours. The proposed viscoacoustic anisotropic equation can keep consistent velocity and attenuation anisotropy effects with that of qP-wave in the constant-Q viscoelastic anisotropic theory. For numerical simulations, the staggered-grid pseudo-spectral method is implemented to solve the velocity–stress formulation of wave equation in the time domain. The constant fractional-order Laplacian approximation method is used to cope with spatial variable-order fractional Laplacians for efficient modelling in heterogeneous velocity and Q media. Simulation results for a homogeneous model show the decoupling of velocity dispersion and amplitude loss effects of the constant-Q equation, and illustrate the influence of anisotropic attenuation on seismic wavefields. The modelling example of a layered model illustrates the accuracy of the constant fractional-order Laplacian approximation method. Finally, the Hess vertical transversely isotropic model is used to validate the applicability of the formulation and algorithm for heterogeneous media.     

层状黏弹性介质中Rayleigh波频散曲线“交叉”现象分析

Rayleigh波勘探方法在探测近地表横波速度、动力学特征等环境与工程地球物理领域获得了广泛应用.这种方法以弹性层状介质理论为基础,然而实际介质具有黏弹性,研究面波在层状黏弹性介质中的传播特征,将为近地表面波勘探提供有益帮助.在某些弹性层状介质模型中,例如存在低速夹层和强波阻抗差异地层模型,Rayleigh波相邻两条频散曲线彼此会非常靠近,产生看似彼此"交叉"的现象,即"osculation"现象,但对于黏弹性介质中的这种现象并没有进行相关的研究.本文利用Muller法计算层状黏弹性介质Rayleigh波频散方程,基于层状介质模型中Rayleigh波频散和衰减曲线连续的性质,结合本征位移曲线特征,分析二层黏弹性介质模型中Rayleigh波频散曲线"交叉"现象以及"交叉"点附近的波动特性.结果表明:与弹性介质相比,黏弹性介质中Rayleigh波的波动特性存在明显差异,随着介质对地震波的损耗越来越强,将导致Rayleigh波频散曲线发生"交叉"现象.     

On occurrence of point singularities in orthorhombic media

Degeneracies of the slowness surfaces of shear (and compressional) waves in low-symmetry anisotropic media (such as orthorhombic), known as point singularities, pose difficulties during modelling and inversion, but can be potentially used in the latter as model parameter constraints. I analyse the quantity and spatial arrangement of point singularities in orthorhombic media, as well as their relation to the overall strength of velocity anisotropy. A classification scheme based on the number and spatial distribution of singularity directions is proposed. In normal orthorhombic models (where the principal shear moduli are smaller than the principal compressional moduli), point singularities can only be arranged in three distinct patterns, and media with the theoretical minimum (0) and maximum (16) number of singularities are not possible. In orthorhombic models resulting from embedding vertical fractures in transversely isotropic background, only two singularity distributions are possible, in contrast to what was previously thought. Although the total number of singularities is independent of the overall anisotropy strength, for general (non-normal) orthorhombic models, different spatial distributions of singularities become more probable with increasing magnitude of anisotropy.     

黏弹性VTI介质微地震正演模拟与震源机制全波形反演

微地震矩张量能够描述储层岩石破坏的细节,有助于了解水力压裂储层的地质力学特征,震源机制矩张量反演对于非常规油气开发具有重要的作用.水力压裂通常在页岩油气储层中进行,震源区介质表现出速度和衰减各向异性,需要研究衰减各向异性介质中微地震信号的传播规律.基于观测到的微地震记录,通过波形反演能够获得矩张量,忽视震源区介质的衰减和速度各向异性会导致震源机制反演结果的误差.本文提出了一种考虑地层各向异性吸收衰减作用的震源机制全波形反演方法,对于黏弹性VTI介质,采用GPU并行的伪谱法计算合成地震记录,通过全波形反演方法匹配合成波形和记录波形,借助优化算法获得震源机制矩张量反演结果.本文通过实验测试验证了该方法的有效性,分析了震源区各向异性速度和衰减参数、噪声以及观测系统等对于微地震震源机制各分量及成分反演结果的影响,研究表明考虑模型中的黏弹性各向异性对于获得可靠的震源机制反演结果非常重要.