首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 421 毫秒
1.
There are two main sources of non-orthogonality in multicomponent shear-wave seismics: inherent non-orthogonal split shear waves arising from substantial ray deviation in off-symmetry planes due to strong anisotropy or complex overburden, and apparent non-orthogonal split shear waves in the horizontal plane due to variation of the angle of incidence even if the two shear waves along the raypath are orthogonal. Many techniques for processing shear-wave splitting in VSP data ignore these kinds of non-orthogonality of the split shear waves. Assuming inherent non-orthogonality in zero-offset VSPs, and apparent non-orthogonality in offset VSPs, we derive equations for the four-component data matrix. These can be solved by extending the linear-transform technique (LTT) to determine the shear-wave polarizations in zero-offset and offset VSPs. Both full-wave synthetic and field data are used to evaluate the technique and to examine the effects of non-orthogonal polarized split shear waves. If orthogonality is incorrectly assumed, errors in polarization measurements increase with the degree of non-orthogonality, which introduces a consistent decreasing trend in the polarization measurements. However, the effect of non-orthogonality on the estimation of geophone orientation and time delays of the two split shear waves is small and negligible in most realistic cases. Furthermore, for most cases of weak anisotropy (less than 5% shear-wave anisotropy) apparent non-orthogonality is more significant than inherent non-orthogonality. Nevertheless, for strong anisotropy (more than 10% shear-wave anisotropy) with complicated structure (tilted or inclined symmetry axis), inherent non-orthogonality may no longer be negligible. Applications to both synthetic and real data show that the extended linear-transform techniques permit accurate recovery of polarization measurements in the presence of both significant inherent and apparent non-orthogonality where orthogonal techniques often fail.  相似文献   

2.
This is the final paper in a series on the 3D multicomponent seismic experiment in Oman. In this experiment a 3D data set was acquired using three-component geophones and with three source orientations. The data set will subsequently be referred to as the Natih 9C3D data set. We present, for the first time, evidence demonstrating that shear waves are sensitive to fluid type in fractured media. Two observations are examined from the Natih 9C3D data where regions of gas are characterized by slow shear-wave velocities. One is that the shear-wave splitting map of the Natih reservoir exhibits much larger splitting values over the gas cap on the reservoir. This increase in splitting results from a decrease in the slow shear-wave velocity which senses both the fractures and the fracture-filling fluid. Using a new effective-medium model, it was possible to generate a splitting map for the reservoir that is corrected for this fluid effect. Secondly, an anomaly was encountered on the shear-wave data directly above the reservoir. The thick Fiqa shale overburden exhibits a low shear-wave velocity anomaly that is accompanied by higher shear reflectivity and lower frequency content. No such effects are evident in the conventional P-wave data. This feature is interpreted as a gas chimney above the reservoir, a conclusion supported by both effective-medium modelling and the geology.
With this new effective-medium model, we show that introduction of gas into vertically fractured rock appears to decrease the velocity of shear waves (S2), polarized perpendicular to the fracture orientation, whilst leaving the vertical compressional-wave velocity largely unaffected. This conclusion has direct implications for seismic methods in exploration, appraisal and development of fractured reservoirs and suggests that here we should be utilizing S-wave data, as well as the conventional P-wave data, as a direct hydrocarbon indicator.  相似文献   

3.
Observations of shear-wave splitting at seismic stations above a swarm of small earthquakes on Hainan Island, China, and other examples world-wide, suggest that the time-delays of split shear-waves monitor the build up of stress before earthquakes and the stress release as earthquakes occur. Rock physics experiments on marble specimens also show variations of shear-wave time-delays with uniaxial pressure analogous to the field observations. The rock experiments show an abrupt decrease in time-delays immediately before fracturing occurs. Similar precursory behaviour has been observed before earthquakes elsewhere, and is believed to be important for two reasons. Precursory changes in shear-wave splitting could be used for short-term forecasting, but of greater importance may be the information such behaviour provides about the source processes in earthquake preparation zones.  相似文献   

4.
多分量地震波波场分解研究   总被引：17，自引：2，他引：15       下载免费PDF全文
提出了一种有效的多分量地震波波场分解方法, 该方法假设在相邻道集波场的振幅、相位变化不大的情况下, 根据纵波和转换波在传播中相对于不同的极化矢量和慢度得到的垂向和水平向分量, 在最小方差约束下从多分量输入数据中获取纵波和转换波的偏振方向和振幅并重构纵波和转换波, 从而分离纵波和转换波. 这种方法精度较高, 不仅能有效分解一次反射波的多分量波场, 多次反射的纵波和转换波也能被恰当地分解. 该方法不需要输入地层速度等参数信息, 减少了人为误差的影响. 通过对合成资料的处理, 验证了方法的可行性和正确性.  相似文献   

5.
The most difficult part of multicomponent processing is the estimation of the shear-wave velocity map for migration. We used refracted shear waves and a simple iterative method called wavefield continuation (WFC) to evaluate the shallow shear-wave velocity profile on a real data example. The WFC was developed in 1981 by Clayton and McMechan to determine compressional-wave velocity profiles from refracted compressional waves. The application to refracted shear waves is straightforward. The real data example shows that shear structure can be easily determined independently of the compressional structure.  相似文献   

6.
The attenuation of seismic waves propagating in reservoirs can be obtained accurately from the data analysis of vertical seismic profile in terms of the quality-factor Q. The common methods usually use the downgoing wavefields in vertical seismic profile data. However, the downgoing wavefields consist of more than 90% energy of the spectrum of the vertical seismic profile data, making it difficult to estimate the viscoacoustic parameters accurately. Thus, a joint viscoacoustic waveform inversion of velocity and quality-factor is proposed based on the multi-objective functions and analysis of the difference between the results inverted from the separated upgoing and downgoing wavefields. A simple separating step is accomplished by the reflectivity method to obtain the individual wavefields in vertical seismic profile data, and then a joint inversion is carried out to make full use of the information of the individual wavefields and improve the convergence of viscoacoustic full-waveform inversion. The sensitivity analysis of the different wavefields to the velocity and quality-factor shows that the upgoing and downgoing wavefields contribute differently to the viscoacoustic parameters. A numerical example validates our method can improve the accuracy of viscoacoustic parameters compared with the direct inversion using full wavefield and the separate inversion using upgoing or downgoing wavefield. The application on real field data indicates our method can recover a reliable viscoacoustic model, which helps reservoir appraisal.  相似文献   

7.
利用二维非均匀介质地震波传播的伪谱和有限差分混合方法,通过数值计算,讨论了松软覆盖层对隐伏断层带围陷波特征的影响.在没有覆盖层的情况下,围陷波振幅和围岩上相比明显增加,持续时间变长.覆盖层造成围岩上地面运动振幅增大,围陷波的部分能量传播到覆盖层中,使得围陷波的能量变小.随着覆盖层厚度增加,围陷波的振幅越来越小,和围岩上...  相似文献   

8.
Seismic anisotropy which is common in shale and fractured rocks will cause travel-time and amplitude discrepancy in different propagation directions. For microseismic monitoring which is often implemented in shale or fractured rocks, seismic anisotropy needs to be carefully accounted for in source location and mechanism determination. We have developed an efficient finite-difference full waveform modeling tool with an arbitrary moment tensor source. The modeling tool is suitable for simulating wave propagation in anisotropic media for microseismic monitoring. As both dislocation and non-double-couple source are often observed in microseismic monitoring, an arbitrary moment tensor source is implemented in our forward modeling tool. The increments of shear stress are equally distributed on the staggered grid to implement an accurate and symmetric moment tensor source. Our modeling tool provides an efficient way to obtain the Green’s function in anisotropic media, which is the key of anisotropic moment tensor inversion and source mechanism characterization in microseismic monitoring. In our research, wavefields in anisotropic media have been carefully simulated and analyzed in both surface array and downhole array. The variation characteristics of travel-time and amplitude of direct P- and S-wave in vertical transverse isotropic media and horizontal transverse isotropic media are distinct, thus providing a feasible way to distinguish and identify the anisotropic type of the subsurface. Analyzing the travel-times and amplitudes of the microseismic data is a feasible way to estimate the orientation and density of the induced cracks in hydraulic fracturing. Our anisotropic modeling tool can be used to generate and analyze microseismic full wavefield with full moment tensor source in anisotropic media, which can help promote the anisotropic interpretation and inversion of field data.  相似文献   

9.
The scattering of elastic energy by random clusters of fractures and/or cavities in a massive rock is studied. The interpretation of the scattered seismic response reveals crucial information about the clusters of inhomogeneities (fractures/cavities), which may correspond to reservoirs. The study is based on a new two‐dimensional numerical‐modelling method that relaxes the constraints on the location and orientation of the inhomogeneities, accounts for inhomogeneities that have almost no volume but a finite surface area (fractures) and improves the accuracy of the calculation when the size of the inhomogeneities is comparable to the mesh size. It is shown that the nature of the seismic response of zones of diffuse fracturing and/or cavities is associated with the non‐uniformity of micro‐inhomogeneities in such zones; accumulations of these micro‐inhomogeneities are known as clusters. The relationship between the non‐uniformity of micro‐inhomogeneities and the strength of the seismic response has been established and measured. Considerable differences in the structure of the seismic response of zones of diffuse fracturing and diffuse cavities have been identified. Converted PS‐waves dominate in the scattered wavefield associated with fractures. This is explained, as the modelling results show, by a greater transparency of fluid‐filled fractures, which reduces the reflected energy of compressional waves. The wavefield associated with cavities is characterized by the predominance (in terms of strength) of compressional PP‐waves. The strength of converted PS‐waves in the scattered wavefields for both media is approximately the same. On the whole, according to the results of the modelling, the energy of the scattered response of fractured reservoirs is considerably less (about two times) than that of cavernous reservoirs.  相似文献   

10.
In this paper, we present results from the analysis of a multicomponent VSP from a fractured gas reservoir in the Bluebell-Altamont Field, Utah. Our analysis is focused on frequency-dependent anisotropy. The four-component shear-wave data are first band-pass filtered into different frequency bands and then rotated to the natural coordinates so that the fast and slow shear-waves are effectively separated. We find that the polarisations of the fast shear-waves are almost constant over the whole depth interval, and show no apparent variation with frequency. In contrast, the time delays between the split shear-waves decrease as the frequency increases. A linear regression is then applied to fit the time-delay variations in the target and we find that the gradients of linear fits to time delays show a decrease as frequency increases. Finally, we apply a time-frequency analysis method based on the wavelet transform with a Morlet wavelet to the data. The variation of shear-wave time delays with frequency is highlighted in the time-delay and frequency spectra. We also discuss two mechanisms giving rise to dispersion and frequency-dependent anisotropy, which are likely to explain the observation. These are scattering of seismic waves by preferentially aligned inhomogeneneities, such as fractures or fine layers, and fluid flow in porous rocks with micro-cracks and macro-fractures.  相似文献   

11.
各向异性介质qP波传播描述II:分离纯模式标量波   总被引：1，自引：1，他引：0       下载免费PDF全文
在各向异性地震波场中,qP波与qS波常常是耦合在一起的.多分量地震数据处理中一个关键环节就是波型分离(即模式解耦),以纵波成分为主的常规单分量地震数据的成像则需要合理描述标量qP波的传播算子.本文作者曾构建了在运动学上同弹性波动方程等价,动力学上突出标量qP波的伪纯模式波动方程.为了彻底消除qS波残余,本文根据波矢量与qP波偏振矢量之间的偏差,提出从伪纯模式波场提取纯模式标量qP波的方法.数值分析展示了投影偏差算子在波数域和空间域的特征.基于不同复杂程度理论模型的试验结果表明,联合"伪纯模式传播算子"与"投影偏差校正"可为各向异性介质分离模式波场传播过程提供一种简便的描述工具.  相似文献   

12.
裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.  相似文献   

13.
A vector convolutional model for multicomponent data acquired in an anisotropic earth is used as a basis for developing algebraic solutions to interpret near-offset VSP data. This interpretation of the cumulative or interval medium response (Green's tensor) for shear waves, determines a polarization azimuth for the leading shear wave and the time-delay between the fast and slow split waves. The algebraic solutions effectively implement least-squares eigenanalysis or singular value decomposition. Although the methodology for shear-wave analysis is strictly relevant to a transmission response, it can be adapted to surface data for a uniform anisotropic overburden. The techniques perform well when calibrated and tested using synthetic seismograms from various anisotropic models. Noise tests demonstrate the sensitivity of the interval measurements to local interferences, particularly if the shear waves are generated by one source. Although the algorithms are faster than numerical search routines, this is not seen as their major advantage. The solutions may have potential in near real-time interpretation of shear-wave data in well logging, where they may be coded on a microchip to provide a direct stream of separated shear waves, or polarization and birefringence information. There may also be some benefit for large prestack multicomponent surface data sets, where the solutions provide a direct transformation to the split-shear-wave components, reducing the storage space for further processing.  相似文献   

14.
随着我国勘探开发难度逐步增大,勘探目标开始向裂缝油气藏、岩性油气藏等复杂探区转移,研究高精度、适应性强的多波多分量深度偏移算法在后续的地震解释、属性分析及储层预测中具有重要意义.针对多波多分量地震数据,本文提出了一种二维弹性波时空域高斯束偏移方法.时空域高斯束沿中心射线传播时能够面向成像目标描述局部波场,且对振幅和频率可调制的Gabor基函数有天然的适应性,因而将基于Gabor分解的子波重构方法应用于震源波场构建,从而得到任意点源函数产生的时空域高斯束波场.该方法由于直接在时间域进行计算,可以避开频率域中出现的假频和边缘截断效应等问题.基于各向同性弹性波动方程的Kirchhoff-Helmholtz积分解,利用矢量时空域高斯束传播算子构建格林函数和格林位移张量,并结合上行射线追踪策略,实现了检波点波场的反向延拓.针对矢量波成像问题,本文借鉴弹性波逆时偏移方法从矢量延拓波场中分离出纯纵波分量和纯横波分量,进而采用修改后的内积成像条件产生具有明确物理意义的PP、PS成像结果,避免了转换波成像的极性反转问题.最后利用简单两层模型和不含盐体构造的部分Sigsbee2a模型的成像结果,并将其与应用近似纵横波成像条件、标量和矢量势成像条件的偏移剖面进行对比,验证了本文方法的正确性和有效性.  相似文献   

15.
The most diagnostic effect of anisotropy on shear waves is shear-wave splitting. This phenomenon creates a distinctive signature in the 3D particle motion. Methods to extract the effects of anisotropy from shear-wave data attempt to measure details of this motion. Many techniques have been published recently which process the shear waves in the time or frequency domain. Here we examine the way in which information on the interference effects between the split shear waves is contained within the frequency domain, and suggest some criteria which may be used in future processing algorithms. The time-delay between the split shear waves, and the polarization direction of the leading shear wave can be converted into easily measured features from analysis of the Fourier spectrum of the shear-wave signal on each component of motion. These features arise in the spectral interference patterns which are formed by the interaction between the two closely-spaced and similar waveforms. The interference patterns are interpreted for synthetic and observed seismogram data.  相似文献   

16.
A marine VSP is designed to estimate the orientation and density of fracturing within a gas-producing dolomite layer in the southern North Sea. The overburden anisotropy is firstly estimated by analysing shear waves converted at or just below the sea-bed, from airgun sources at four fixed offset azimuths. Full-wave modelling helps confirm that the background has no more than 3% vertical birefringence, originating from TIH anisotropy with a symmetry axis orientated perpendicular to the maximum horizontal compressive stress of NW–SE. This finding concurs with current hypotheses regarding the background rock matrix in the upper crust. More detailed anisotropy estimates reveal two thin zones with possible polarization reversals and a stronger anisotropy. The seismic anisotropy of the dolomite is then determined from the behaviour of locally converted shear waves, providing a direct link with the physical properties of its fractures. It is possible to utilize this phenomenon due to the large seismic velocity contrast between the dolomite and the surrounding evaporites. Two walkaway VSPs at different azimuths, recorded on three-component receivers placed inside the target zone, provide the appropriate acquisition design to monitor this behaviour. Anisotropy in the dolomite generates a transverse component energy which scales in proportion to the degree of anisotropy. The relative amplitudes, for this component, between the different walkaway azimuths relate principally to the orientation of the anisotropy. Full-wave modelling confirms that a 50% vertical birefringence from TIH anisotropy with a similar orientation to the overburden is required to simulate the field observations. This amount of anisotropy is not entirely unexpected for a fine-grained brittle dolomite with a potentially high fracture intensity, particularly if the fractures contain fluid which renders them compliant to the shear-wave motion.  相似文献   

17.
The study of seismic body waves is an integral aspect in global, exploration and engineering scale seismology, where the forward modeling of waves is an essential component in seismic interpretation. Forward modeling represents the kernel of both migration and inversion algorithms as the Green’s function for wavefield propagation and is also an important diagnostic tool that provides insight into the physics of wave propagation and a means of testing hypotheses inferred from observational data. This paper introduces the one-way wave equation method for modeling seismic wave phenomena and specifically focuses on the so-called operator-root one-way wave equations. To provide some motivation for this approach, this review first summarizes the various approaches in deriving one-way approximations and subsequently discusses several alternative matrix narrow-angle and wide-angle formulations. To demonstrate the key strengths of the one-way approach, results from waveform simulation for global scale shear-wave splitting modeling, reservoir-scale frequency-dependent shear-wave splitting modeling and acoustic waveform modeling in random heterogeneous media are shown. These results highlight the main feature of the one-way wave equation approach in terms of its ability to model gradual vector (for the elastic case) and scalar (for the acoustic case) waveform evolution along the underlying wavefront. Although not strictly an exact solution, the one-way wave equation shows significant advantages (e.g., computational efficiency) for a range of transmitted wave three-dimensional global, exploration and engineering scale applications.  相似文献   

18.
Results of a passive microtremor survey at an oil and gas field in Voitsdorf, Austria, are presented. The survey consists in six parallel profiles approximately 9 km long over two hydrocarbon reservoirs. For each profile the seismic wavefield was recorded synchronously at 11 in-line stations. The measurements were conducted with broadband seismometers and lasted, for each profile, at least 12 hours overnight. Data interpretation is based on a comprehensive data set and on the analysis of four different spectral attributes. These attributes quantify the characteristic features of the wavefield's Fourier spectra in the low-frequency range (<10Hz). One attribute quantifies the spectral energy in the vertical wavefield component, another attribute quantifies the maxima in vertical-to-horizontal spectral ratios and two attributes describe the frequency shifts of peaks within the spectra of vertical and horizontal wavefield components. Due to temporal variations of the signals we combine the long-term measurements (several hours of continuous records) of multiple profiles. This procedure considerably enhances the consistency of each spectral attribute and makes them suitable to quantify lateral variations of the wavefield. The results show that using a combination of several attributes significantly increases the reliability of detecting anomalies in the microtremor wavefield that are presumably caused by hydrocarbon reservoirs. A numerical study of two-dimensional seismic wave propagation is applied to investigate the peak frequency shift attributes. The results of the study indicate that the attributes may contain information on the depth of hydrocarbon reservoirs, assuming that the reservoir acts as a (secondary) source of low-frequency seismic waves.  相似文献   

19.
A global optimization method incorporating a ray-tracing scheme is used to invert observations of shear-wave splitting from two near-offset VSPs recorded at the Conoco Borehole Test Facility, Kay County, Oklahoma. Inversion results suggest that the seismic anisotropy is due to a non-vertical fracture system. This interpretation is constrained by the VSP acquisition geometry for which two sources are employed along near diametrically opposite azimuths about the well heads. A correlation is noted between the time-delay variations between the fast and slow split shear waves and the sandstone formations.  相似文献   

20.
Although the shear-wave birefringence phenomenon affects the imaging of converted shear waves, it also provides a considerable amount of information on subsurface fracture development. Therefore, it is significant to separate split shear waves before seismic interpretation and reservoir prediction. In this paper, we propose a new method of split shear waves separation based on the polarization directions derived from hodogram analysis. Through the hodogram analysis, we find that the split shear-wave particle motions are within the range of a specific and fixed rectangle, which have relations with the fracture azimuth in strata. In addition, we found that a couple of split shear waves can only be fitted to the unique trajectory rectangle through the theoretical derivation. Based on this, we establish the trajectory rectangle through the wave vector calculation and calculate the fracture azimuth according to the fact that the one edge of the trajectory rectangle is along or perpendicular to the fracture azimuth. Synthetic data analysis shows that the calculation accuracy of fracture azimuth under the constraint of trajectory rectangle is less affected by the time delay between split shear waves than using the method of eigenvector–eigenvalue decomposition (EED). Therefore, we can obtain better results for separation of split shear waves using our method than using EED. Eventually, we propose an approach of layer stripping to deal with the problem that shear wave split several times due to the situation that different strata have different fracture azimuths. Synthetic data test indicates that our method can achieve higher calculation efficiency and faster convergence speed than the conventional eigenvector–eigenvalue decomposition method, even though the data are of a low signal-to-noise ratio. Moreover, field data applications show the effectiveness and potential of our method.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号