共查询到19条相似文献，搜索用时 144 毫秒
1.
岩石圈各向异性主要由上地幔矿物晶格优势排列方向和上地壳内裂缝、裂隙的定向分布造成.在各向异性特征显著区域,利用SKS震相剪切波分裂获得的延迟时间高达1.5s以上.本文根据方位各向异性,利用广义反射透射系数矩阵方法正演S波接收函数,研究各向异性对不同反方位角接收函数转换震相走时的影响.我们发展了基于HTI模型各向异性走时校正的方法,成功的在单层和多层(快轴方向相同或不同)的各向异性介质中对齐不同反方位角接收函数的Moho面和LAB的转换震相走时.我们将该方法应用于在青藏高原东北缘的流动台站,试图在各向异性强度较大区域对实测数据转换震相走时的校正效果进行测试.结果表明:各向异性走时校正能够加强单台接收函数转换震相的可追踪性,能量增强的叠后转换震相在时深转换后更利于对界面深度的识别与判断;在考虑三维成像的情况下,我们的各向异性校正方法对提高成像结果的准确性有重要意义. 相似文献
2.
利用横波响应进行裂缝性各向异性介质的检测在实际应用中取得了很好的效果,但技术复杂、成本较高使该方法的广泛使用受到限制,而纵波资料采集和处理技术的精细有效保持了纵波的各种属性,这为直接利用纵波资料进行裂缝检测创造了条件.地下垂直定向裂缝通常用HTI介质模型来描述,为此,本文利用射线追踪和反射率法计算了层状各向同性介质背景下的HTI介质顶、底界面的反射纵波旅行时和反射系数,并分析了这些属性随观测方位的变化规律.研究表明,HTI介质底界面反射纵波旅行时和HTI介质顶界面反射系数表现出了明显的方位各向异性;旅行时、振幅和AVO梯度属性均在0°观测方位和90°观测方位上存在最大差异,可以用多种属性联合来精确判定裂缝的发育方向. 相似文献
3.
对Christoffel公式进行Bond变换得到EDA介质的Christoffel方程,并由其非零解推导出EDA介质中视横波(qSV)、横波(SH)、视纵波(qP)的相速度、群速度、偏振向量(质点的振动方向)的三维计算公式.通过模型计算分析了具有水平对称轴的各向异性(HTI)介质和EDA介质中介质对称轴的极角和方位角对相速度、群速度及偏振向量的影响,对其随极角、方位角的变化特征进行了分析,并采用Matlab进行了数值计算,对其特征采用三维显示.通过取极角或方位角为零简化得到HTI介质和具有垂直对称轴的各向异性(VTI)介质中地震波的相速度、群速度,对EDA介质中的三维计算结果进行退化验证. 通过数值计算进一步验证了地震波相速度与EDA介质对称轴的相互关系. 结果表明,通过广角地震勘探可探明地下介质的裂隙走向及密度,从而确定灾害体产状. 相似文献
4.
正交各向异性介质P波走时分析及Thomsen参数反演 总被引：10，自引：3，他引：10
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层，正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时，时距曲线表现出强方位各向异性.采用小生境遗传算法，对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明，此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数. 相似文献
5.
转换波方位各向异性裂缝检测技术研究及应用 总被引：1，自引：0，他引：1
HTI裂缝各向异性介质中,转换波随方位角的变化比较复杂,目前还没有解析公式可以表达其变化特征,只能通过物理实验或数值模拟来分析其应用的可能性.数值模拟结果表明,转换波在裂缝各向异性介质中传播时,其R分量和T分量的振幅属性都具有明显的方位各向异性特征,R分量振幅方位各向异性拟合椭圆的长轴方向指示裂缝方位,这与纵波方位各向异性特征相似;根据P波AVAZ方位各向异性分析原理,对转换波R分量振幅方位各向异性曲线进行方位椭圆拟合,寻找椭圆的长轴方向,即裂缝主方位,再由P波AVAZ技术中振幅响应与炮检方向和裂缝走向之间的夹角关系式得到裂缝的发育密度,从而构建转换波方位各向异性AVAZ裂缝检测技术.该技术已用于川西新场气田某区块的裂缝储层预测,取得了较好的应用效果. 相似文献
6.
在裂缝诱导各向异性理论研究中通常使用等效HTI介质来近似多组裂缝所引起的综合效应.由于构造运动的复杂性,多组裂缝普遍存在于地壳与油气储层中.为了研究多组裂缝的地震属性特征,分析常用的等效HTI模型对于多组裂缝近似精度及附加裂缝对介质属性特征的影响,本文利用线性滑移模型进行了多组垂直裂缝的单斜各向异性等效介质理论计算,并利用空间搜索方法求取与其最为接近的HTI介质各向异性弹性参数.重点研究了在两种各向异性介质中纵波速度、快慢横波速度和极化特征及其差异,量化分析附加裂缝对于地震属性如速度、极化方向和走时等的影响,研究对附加裂缝敏感的地震属性.此研究结果和方法为进一步研究多组裂缝的反演及识别方法提供基础,同时对于将高阶对称性各向异性介质中已存在的计算方法应用于低阶对称性时的适用程度、精度分析及相关方法研究具有重要作用. 相似文献
7.
常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。 相似文献
8.
利用地震数据较为准确地预测地下介质的裂缝发育信息,是裂缝型储层预测的关键手段,也是非常规页岩油气储层压裂改造的关键参数.方位AVO(AVAZ)在描述HTI介质的裂缝分布及发育方向等方面有重要应用,本文提供一种基于全方位地震数据的方位AVO(AVAZ)反演方法,为裂缝预测提供参考资料.首先利用Ruger近似方程正演计算不同方位角及入射角下的反射系数,继而与方位角度道集中地震反射振幅建立目标函数,然后基于正演方程及目标函数采用最小二乘法进行AVAZ反演.实现了从全方位地震资料中同步反演纵波反射系数、各向同性梯度、各向异性梯度及裂缝发育方向.模型数据和实际地震资料应用表明该方法能够有效的预测HTI介质的裂缝分布及方向. 相似文献
9.
基于AVO梯度、走时及速度的纵波方位角裂缝检测已有较多研究和应用.由于裂缝及所含流体的综合作用,碳酸岩裂缝性储层的地震响应具有高频衰减、频带低移及方位各向异性的特点,因此应用衰减属性进行裂缝预测具有重要价值.基于小波变换计算的瞬时频域振幅谱及以此为基础提取的衰减属性,由于在时域和频域都有较好的分辨率能力,能够准确地刻画出碳酸岩裂缝性储层的非均质性.结合三维地震纵波数据的分方位角处理和分析技术,可以由方位衰减属性估算裂缝相对密度.对Z区下古生界潜山碳酸岩储层进行的裂缝预测表明,由衰减梯度和85%能量比频率两种属性导出的裂缝发育密度具有很好的相关性,且与实际情况符合. 相似文献
10.
煤层的含气性及渗透性可以利用纵波AVO技术及基于各向异性理论的方位AVO技术进行研究.煤层气属于吸附气,储层厚度较薄,地震异常响应较弱,对地震资料的品质及其携带的信息量提出了更高的要求.本文对AVO与方位AVO反演理论的近似条件及模型试算结果进行了分析,提出了对地震纵波数据信噪比、入射角、偏移距分布及方位角分布等参数的要求:1.观测系统设计及野外采集时应保证地震波的有效入射角在0.～30.或更大范围内,偏移距应大中小均匀分布,避免过于集中或缺失.2.地震数据应有足够宽分布的方位角,煤层气地震勘探中应选择宽方位观测系统.3、地震资料信噪比应足够高. 相似文献
11.
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP- and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities. 相似文献
12.
提出了利用人工爆破P波走时反演地壳介质方位各向异性参数的方法. 在假定介质是弱各向异性介质的情况下, 使用扰动理论得到了线性化的反演公式, 其中待反演的弱各向异性参数是P波走时的线性函数. 如果在反演公式中参考走时取相同震中距接收点的P波平均走时, 那么所获得的弱各向异性参数与参考介质速度的选取无关. 反演得到的弱各向异性参数可以看作是不同震中距和不同深度范围内介质的等效弱各向异性参数. 等效弱各向异性参数在一定程度上反映了不同深度范围内水平方向相速度随方位的变化. 这种变化可能是不同时期构造应力作用的结果. 2007年中国地震局在首都圈怀来地区实施了一次大吨位人工爆破实验, 以爆破点为中心, 布设了高密度的地震观测台网和台阵. 台站相对于爆破点具有360°的全方位覆盖, 所得到的地震记录数据为研究怀来、 延庆地区地壳介质P波方位各向异性提供了必要条件. 我们通过走时反演获得了与水平方位相关的弱各向异性参数, 并对弱各向异性参数进行坐标变换, 得到了能够直观描述岩石弱各向异性的具有水平对称轴的横向各向同性介质, 给出了对应的3个独立弱各向异性参数及其对称轴方位, 讨论了介质各向异性与构造应力场的关系. 结果表明该地区地壳介质存在明显的方位各向异性, 其最大值约为4.6%. 相似文献
13.
Bayesian Markov chain Monte Carlo inversion for anisotropy of PP- and PS-wave in weakly anisotropic and heterogeneous media 
下载免费PDF全文
下载免费PDF全文 A single set of vertically aligned cracks embedded in a purely isotropic background may be considered as a long-wavelength effective transversely isotropy (HTI) medium with a horizontal symmetry axis. The crack-induced HTI anisotropy can be characterized by the weakly anisotropic parameters introduced by Thomsen. The seismic scattering theory can be utilized for the inversion for the anisotropic parameters in weakly anisotropic and heterogeneous HTI media. Based on the seismic scattering theory, we first derived the linearized PP- and PS-wave reflection coefficients in terms of P- and S-wave impedances, density as well as three anisotropic parameters in HTI media. Then, we proposed a novel Bayesian Markov chain Monte Carlo inversion method of PP- and PS-wave for six elastic and anisotropic parameters directly. Tests on synthetic azimuthal seismic data contaminated by random errors demonstrated that this method appears more accurate, anti-noise and stable owing to the usage of the constrained PS-wave compared with the standards inversion scheme taking only the PP-wave into account. 相似文献
14.
We use residual moveouts measured along continuous full azimuth reflection angle gathers, in order to obtain effective horizontal transversely isotropic model parameters. The angle gathers are generated through a special angle domain imaging system, for a wide range of reflection angles and full range of phase velocity azimuths. The estimation of the effective model parameters is performed in two stages. First, the background horizontal transversely isotropic (HTI)/vertical transversely isotropic (VTI) layered model is used, along with the values of reflection angles, for converting the measured residual moveouts (or traveltime errors) into azimuthally dependent normal moveout (NMO) velocities. Then we apply a digital Fourier transform to convert the NMO velocities into azimuthal wavenumber domain, in order to obtain the effective HTI model parameters: vertical time, vertical compression velocity, Thomsen parameter delta and the azimuth of the medium axis of symmetry. The method also provides a reliability criterion of the HTI assumption. The criterion shows whether the medium possesses the HTI type of symmetry, or whether the azimuthal dependence of the residual traveltime indicates to a more complex azimuthal anisotropy. The effective model used in this approach is defined for a 1D structure with a set of HTI, VTI and isotropic layers (with at least one HTI layer). We describe and analyse the reduction of a multi‐layer structure into an equivalent effective HTI model. The equivalent model yields the same NMO velocity and the same offset azimuth on the Earth's surface as the original layered structure, for any azimuth of the phase velocity. The effective model approximates the kinematics of an HTI/VTI layered structure using only a few parameters. Under the hyperbolic approximation, the proposed effective model is exact. 相似文献
15.
Constraining the anisotropy structure of the crust by joint inversion of seismic reflection travel times and wave polarizations 总被引：2，自引：0，他引：2
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. 相似文献
16.
Three‐dimensional residual moveout analysis is the basic step in velocity model refinement. The analysis is generally carried out using horizontal and/or vertical semblances defined on a sparse set of in‐lines or cross‐lines with densely sampled source–receiver offsets. An alternative approach, which we call dense residual moveout analysis (DRMA), is to use all the bins of a three‐dimensional survey but sparsely sampled offsets. The proposed technique is very fast and provides unbiased and statistically efficient estimates of the residual moveout. Indeed, for the sparsest possible offset distribution, when only near‐ and far‐angle stacks are used, the variance of the residual moveout estimate is only 1.4 times larger than the variance of the least‐squares estimate obtained using all offsets. The high performance of DRMA makes it a useful tool for many applications, of which azimuthal velocity analysis is considered here. For a horizontal transverse isotropy (HTI) model, a deterministic procedure is proposed to define, at every point of residual moveout estimation, the azimuthal angle of the HTI axis of symmetry, the Thomsen anisotropy coefficients, and the interval (or root‐mean‐square) velocities in both the HTI isotropy and symmetry planes. The procedure is not restricted by DRMA assumptions; for example, it is also applicable to semblance‐based residual moveout estimates. The high resolution of the technique is illustrated by azimuthal velocity analysis over an oilfield in West Siberia. 相似文献
17.
Effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves 下载免费PDF全文
下载免费PDF全文 With the development of the dense array, the surface wave velocity and azimuthal anisotropy under the array can be directly obtained by beamforming the noise cross-correlation functions (NCFs). However, the retrieval of the Green's function by cross-correlating the seismic noise requires that the noise source has a uniform distribution. For the case with uneven noise source, the azimuthal dependence on the sources in the expression for the spatial coherence function, which corresponds to the NCF in the time domain, has the same form as the azimuthal dependence of the surface wave velocity in weakly anisotropic media. Therefore, the uneven noise source will affect the surface wave anisotropy extraction. In this study, three passive seismic methods, i.e., beamforming, SPAC (spatial autocorrelation), and NCF, are compared to demonstrate that an uneven source distribution and uneven station distribution have equivalent effects on the outcome from each method. A beamforming method is proposed to directly extract the velocity and azimuthal anisotropy of surface waves. The effect of uneven noise source and/or station distribution on estimating the azimuth anisotropy of surface waves was investigated using data from the ChinArray Phase II. A method for correcting the apparent anisotropy in beamforming results caused by an uneven station distribution is suggested. 相似文献
18.
19.
Estimating azimuthal stress‐induced P‐wave anisotropy from S‐wave anisotropy using sonic log or vertical seismic profile data 下载免费PDF全文
下载免费PDF全文 Most sedimentary rocks are anisotropic, yet it is often difficult to accurately incorporate anisotropy into seismic workflows because analysis of anisotropy requires knowledge of a number of parameters that are difficult to estimate from standard seismic data. In this study, we provide a methodology to infer azimuthal P‐wave anisotropy from S‐wave anisotropy calculated from log or vertical seismic profile data. This methodology involves a number of steps. First, we compute the azimuthal P‐wave anisotropy in the dry medium as a function of the azimuthal S‐wave anisotropy using a rock physics model, which accounts for the stress dependency of seismic wave velocities in dry isotropic elastic media subjected to triaxial compression. Once the P‐wave anisotropy in the dry medium is known, we use the anisotropic Gassmann equations to estimate the anisotropy of the saturated medium. We test this workflow on the log data acquired in the North West Shelf of Australia, where azimuthal anisotropy is likely caused by large differences between minimum and maximum horizontal stresses. The obtained results are compared to azimuthal P‐wave anisotropy obtained via orthorhombic tomography in the same area. In the clean sandstone layers, anisotropy parameters obtained by both methods are fairly consistent. In the shale and shaly sandstone layers, however, there is a significant discrepancy between results since the stress‐induced anisotropy model we use is not applicable to rocks exhibiting intrinsic anisotropy. This methodology could be useful for building the initial anisotropic velocity model for imaging, which is to be refined through migration velocity analysis. 相似文献