VTI介质中的地层弹性参数反演

在VTI介质中，用非双曲线可以很好地近似P－P波中长排列情况下的旅行时；而对P－SV波来说，应用非双曲线却不能很好地近似P－SV波的施行时。对P－P波、P－SV波的旅行时进行三项泰勒展开时，它们的系数却包含了地震波的各向异性信息及其垂直速度信息，这为地震波的弹性参数反演提供了基础。实验分析表明，在中等排列长度下可以精确地得到了P－P波，P－SV波旅行时的三项泰勒展式系数，从而可以精确地进行地层弹性参数反演。     

条形均布荷载下横观各向同性弹性地基中附加应力的分布规律

本文评述了条形均布荷载下横观各向同性弹性地基中附加应力的诸户［２］解析解和近似解，并提出了笔者的近似简化公式，该公式简明地显示出附加应力在分布规律上与各向同性地基之异同点，可供工程界方便地加以利用。     

Quasi-Rayleigh Waves in Transversely Isotropic Half-Space with Inclined Axis of Symmetry

A method for determination of characteristics of quasi-Rayleigh (qR) wave in a transversely isotropic homogeneous half-space with inclined axis of symmetry is outlined. The solution is obtained as a superposition of qP, qSV and qSH waves, and surface wave velocity is determined from the boundary conditions at the free surface and at infinity, as in case of Rayleigh wave in an isotropic half-space. Though the theory is simple enough, a numerical procedure for calculation of surface wave velocity presents some difficulties. The difficulty is caused by necessity to calculate complex roots of a non-linear equation, which in turn contains functions determined as roots of non-linear equations with complex coefficients. Numerical analysis shows that roots of the equation corresponding to the boundary conditions do not exist in the whole domain of azimuths and inclinations of the symmetry axis. The domain of existence of qR wave depends on the ratio of the elastic parameters: for some strongly anisotropic models the wave cannot exist at all. For some angles of inclination qR-wave velocities deviate from those calculated on the basis of the perturbation method valid for weak anisotropy, though they have the same tendency of variation with azimuth. The phase of qR wave varies with depth unlike Rayleigh wave in an isotropic half-space. Unlike Rayleigh wave in an isotropic half-space, qR wave has three components - vertical, radial and transverse. Particle motion in horizontal plane is elliptic. Direction of the major axis of the ellipsis coincides with the direction of propagation only in azimuths 0° (180°) and 90° (270°).     

The reflection and transmission of elastic waves at a interface between a transversely isotropic elastic soils and saturated soils

ｎｔｒｏｄｕｃｔｉｏｎＴｈｅｒｅｆｌｅｃｔｉｏｎａｎｄｔｒａｎｓｍｉｓｉｏｎｏｆｅｌａｓｔｉｃｗａｖｅａｔａｉｎｔｅｒｆａｃｅｂｅｔｗｅｅｎｔｗｏｄｉｆｅｒｅｎｔｍｅｄｉａａｒｅａｎｉｍｐｏｒｔａｎｔａｎｄｂａｓｉｃｑｕｅｓｔｉｏｎｉｎＥａｒｔｈｑ...     

均匀横向各向同性介质中P波及SV波的射线规律

常规的反射数据处理都假定介质为均匀且各向同位。研究表明，各向异性的存在会严重影响处理和解释的结果。在已有知识的基础上，研究了二维均匀横向各向同性介质中的射线传播规律并描述了速度各向异性波场的特征。对于准P波，各向异性参数ε，δ的绝对值决定波前面的形状，而ε与δ之差决定射线密度；对于准SV波，σ则是主要的控制因素。     

Amplitude variation with incident angle and azimuth inversion for Young's impedance,Poisson's ratio and fracture weaknesses in shale gas reservoirs

By analogy with P- and S-wave impedances, the product of Young's modulus and density can be termed as Young's impedance, which indicates the rock lithology and brittleness of unconventional hydrocarbon reservoirs. Poisson's ratio is also an effective indicator of rock brittleness and fluid property of unconventional reservoirs, and fracture weaknesses indicate the fracture properties (fracturing intensity and fracture fillings) in fracture-induced unconventional reservoirs. We aim to simultaneously estimate the Young's impedance, Poisson's ratio and fracture weaknesses from wide-azimuth surface seismic data in a fracture-induced shale gas reservoir, and use the horizontal transversely isotropic model to characterize the fractures. First, the linearized PP-wave reflection coefficient in terms of Young's impedance, Poisson's ratio, density and fracture weaknesses is derived for the case of a weak-contrast interface separating two weakly horizontal transversely isotropic media. In addition, an orthorhombic anisotropic case is also discussed in this paper. Then a Bayesian amplitude variation with incident angle and azimuth scheme with a model constraint is used to stably estimate Young's impedance, Poisson's ratio and fracture weaknesses with only PP-wave azimuthal seismic data. The proposed approach is finally demonstrated on both synthetic and real data sets with reasonable results.     

3D waveform inversion of downhole microseismic data for transversely isotropic media

3D anisotropic waveform inversion could provide high-resolution velocity models and improved event locations for microseismic surveys. Here we extend our previously developed 2D inversion methodology for microseismic borehole data to 3D transversely isotropic media with a vertical symmetry axis. This extension allows us to invert multicomponent data recorded in multiple boreholes and properly account for vertical and lateral heterogeneity. Synthetic examples illustrate the performance of the algorithm for layer-cake and ‘hydraulically fractured’ (i.e. containing anomalies that simulate hydraulic fractures) models. In both cases, waveform inversion is able to reconstruct the areas which are sufficiently illuminated for the employed source-receiver geometry. In addition, we evaluate the sensitivity of the algorithm to errors in the source locations and to band-limited noise in the input displacements. We also present initial inversion results for a microseismic data set acquired during hydraulic fracturing in a shale reservoir.