裂隙参数对衰减各向异性影响的数值模拟

理论、观测和实验均证实,在地壳和上地幔存在对应力变化敏感的直立裂隙,在整体上呈现方位各向异性.在Hudson裂隙理论基础上,系统全面地数值模拟了裂隙介质几何、物性和弹性波频率等参数对各向异性衰减的影响.结果显示,裂隙密度、裂隙纵横比、泊松比、裂缝填充物、弹性波频率和未破裂岩石母体纵横波速度等对各向异性衰减有着显著影响,...     

裂隙参数对衰减各向异性影响的数值模拟

Propagation through stress-aligned fluid-filled cracks and other inclusions have been claimed to be the cause of azimuthal anisotropy observed in the crust and upper mantle.This paper examines the behavior of seismic waves attenuation caused by the internal structure of rock mass,and in particular,the internal geometry of the distribution of fluid-filled openings Systematic research on the effect of crack parameters,such as crack density,crack aspect ratio(the ratio of crack thickness to crack diameter),pore fluid properties(particularly pore fluid velocity),VP/VS ratio of the matrix material and seismic wave frequency on attenuation anisotropy has been conducted based on Hudson’s crack theory.The result shows that the crack density,aspect ratio,material filler,seismic wave frequency,and P-wave and shear wave velocity in the background of rock mass,and especially frequency has great effect on attenuation curves.Numerical research can help us know the effect of crack parameters and is a good supplement for laboratory modeling.However,attenuation is less well understood because of the great sensitivity of attenuation to details of the internal geometry.Some small changes in the characteristics of pore fluid viscosity,pore fluids containing gas and liquid phases and pore fluids containing clay can each alter attenuation coefficients by orders of magnitude.Some parameters controlling attenuation are therefore necessary to make reasonable estimations,and anisotropic attenuation is worth studying further.     

MODELLING CHANNEL WAVES WITH SYNTHETIC SEISMOGRAMS IN AN ANISOTROPIC IN-SEAM SEISMIC SURVEY1

In-seam seismic surveys with channel waves have been widely used in the United Kingdom and elsewhere to map coal-seams and to detect anomalous features such as dirt bands, seam thinning and thickening, and particularly in-seam faulting. Although the presence of cleat-induced anisotropy has been recognized in the past, almost all previous analyses have assumed homogeneous isotropic or transversely isotropic coal-seams. Channel waves, however, exhibit properties which cannot be fully explained without introducing anisotropy into the coal-seam. In particular, Love-type channel waves are observed for recording geometries where, in a homogeneous isotropic or transversely isotropic structure, the source would not be expected to excite transverse motion. Similarly, modes of channel-wave propagation display the coupled three-component motion of generalized modes in anisotropic substrates, which would not be expected for Rayleigh and Love wave motion in isotropy or in transversely isotropic media with azimuthal isotropy. We model the observed in-seam seismic channel waves with synthetic seismograms to gain an understanding of the effects of cleat-induced anisotropy on the behaviour of channel waves. The results show a reasonable good match with the observations in traveltime, relative amplitudes, dispersion characteristics and particle motions. We demonstrate that anisotropy in the surrounding country rocks contributes significantly to the coupling of channel wave particle motion, although its effect is not as strong as the anisotropy in the coal-seam. We conclude that the effects of cleat- and stress-induced anisotropy are observed and can be modelled with synthetic seismograms, and that anisotropy must be taken into account for the detailed interpretation of channel waves.     

裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2-3kHz）和长源距（3-4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

PARAMETER ESTIMATION AND FAULT DETECTION BY THREE-COMPONENT SEISMIC AND GEOELECTRICAL SURVEYS IN A COAL MINE*

Three-component seismic and geoelectrical in-mine surveys were carried out in Lyukobanya colliery near Miskolc, Hungary to determine the in situ petrophysical parameter distributions and to detect inhomogeneities in the coal seam. The seismic measurements comprise an underground vertical seismic profile, using body waves, and an in-seam seismic amplitude-depth distribution and transmission survey, using channel waves. The geoelectrical measurements are based on the drift- and seam-sounding method. Interval traveltime-, amplitude-, multiple-filter- and polarization analysis methods are applied to the seismic data. They lead to a five-layer model for the strata including the coal seam. The coal seam and two underlying beds act as a seismic waveguide. The layer sequence supports the propagation of both normal and leaky mode channel waves of the Love- and Rayleigh type. A calculation of the total reflected energy for each interface using Knott's energy coefficients shows that the velocity ranges of high reflection energy and of normal and leaky mode wavegroups coincide. The excitation of wavegroups strongly depends on the seismic source. A simultaneous inversion of a geoelectrical drift- and seam-sounding survey prevents misinterpretations of the seismic data by clearly identifying the low-velocity coal seam as a high-resistivity bed. Calculations of dispersion and sounding curves improve the resolution of the slowness and resistivity in each layer. Both diminished amplitudes and distortions in the polarization of transmission seismo-grams and decreasing resistivities in a geoelectrical pseudosection of the coal seam are related to an inhomogeneity. A calculation of synthetic seismograms for Love and Rayleigh channel waves with the finite-difference and the Alekseev-Mikhailenko method agrees well with the field data for the main features, i.e., particular arrivals in the wave train, wavegroups, velocities and symmetries or asymmetries. This in-mine experiment demonstrates that the simultaneous acquisition, processing and interpretation of seismic and geoelectrical data improve the lithological interpretation of petrophysical parameter distributions. Coal seam inhomogeneities can also be detected more reliably by the two independent surveys than by one alone.     

煤巷小构造Rayleigh型槽波超前探测数值模拟

对煤巷小构造地震波场进行了数值模拟研究,分析了层状煤层中地震波的传播特征.研究表明:(1)在煤巷迎头前方煤层内以纵波震源激发的Rayleigh型槽波相对于体波能量较强,波列较长,波速较低.(2)沿煤层传播的Rayleigh型槽波在小构造面上产生Rayleigh型槽波反射波,反射Rayleigh型槽波垂直分量相对于水平分量能量较强.沿煤层反向传播的反射Rayleigh型槽波在煤巷迎头面上转换为沿煤巷底板传播的Rayleigh面波.沿煤巷底板可以接收到能量较强的反射Rayleigh型槽波产生的Rayleigh面波,其可以作为超前探测小构造面的特征波.在地震记录上反射Rayleigh型槽波产生的Rayleigh面波波至最迟,在时间域与其他波列时间间隔较大,其垂直分量能量相对于水平分量较强,在地震记录上容易识别.(3)在相同的地质条件下应用反射地震超前探测方法,标志煤巷迎头前方存在小构造面的反射地震波能量较弱,受煤巷顶、底板界面和采煤迎头面的强反射波干扰,在地震记录中难以识别.     

ALIGNMENT OF NEAR-SURFACE INCLUSIONS AND APPROPRIATE CRACK GEOMETRIES FOR GEOTHERMAL HOT-DRY-ROCK EXPERIMENTS1

Ubiquitous splitting of seismic shear-waves indicates that most rocks in the upper half of the crust are pervaded by stress-aligned fluid-filled inclusions, called EDA-cracks. These inclusions are expected to be aligned perpendicular to the minimum compressional stress by stress relationships similar to those aligning industrial hydraulic fractures. At depths where the overburden stress is sufficiently large (typically below a few hundred metres), this minimum stress is usually horizontal, so that the EDA-cracks and hydraulic fractures are typically aligned vertically, striking parallel, or subparallel, to the direction of maximum compression. This is confirmed by the polarizations of the split shear-waves along raypaths at depth in the crust. At the free surface, however, the vertical stress is zero (or approximately zero) and cracks (and hydraulic fractures) at shallow depths in intact rock tend to be horizontal. Thus, the directions of minimum stress, and the orientations of hydraulic fractures, are likely to swing through 90° near the surface of the Earth. Since the behaviour of cracks and stress is often crucial to drilling operations, the rotation of the crack- and stress-geometry near-surface has important implications, particularly for optimizing hydrocarbon production and geothermal reservoir management. Consequently, evidence gained from experiments, for example in hot-dry-rock geothermal heat extraction, in inappropriate crack geometries at shallow depths, may not be valid when applied to other crack- and stress-geometries at depth in hot rock.     

IMAGING CAPABILITY OF CROSS-HOLE SEISMIC REFLECTION SURVEYS1

In seismic migration, it is important to sample a range of dips around the local structural dip at each image point. Meaningful images are obtained only where this condition holds. For cross-hole seismic reflection surveys, the distribution of dips sampled at each image point is controlled principally by the survey geometry, including source and receiver array lengths and their element spacings. Using a real data set as an example, we show how survey geometry can limit imaging capability close to the boreholes and even in the middle of the section between the boreholes. At the processing stage, effective removal of direct waves and accurate estimation of the velocity field are essential for optimizing image quality. For migration, we propose a generalized Berryhill (GB) scheme which is based on the Kirchhoff integral and takes into account both the near-field and far-field terms. This should improve the ability to image close to source and receiver arrays, provided that the element spacing in the nearby array is small enough.     

Shallow geology characterization using Rayleigh and Love wave dispersion curves derived from seismic noise array measurements

The local geology and shallow S-wave velocity structure of a site are recognized to be key factors for the increase in the damaging potential of seismic waves. Indeed, seismic amplitudes may be amplified in frequency ranges unfavorable for building stock by the presence of soft sedimentary covers over lying hard bedrock. Hence, microzonation activities, which aim at assessing the site response as accurately as possible, have become a fundamental task for the seismic risk reduction of urbanized areas. Methods based on the measurement of seismic noise, which typically are fast, non-invasive, and low cost, have become a very attractive option in microzonation studies.Using observations derived from seismic noise recordings collected by two-dimensional arrays of seismic stations, we present a novel joint inversion scheme for surface wave curves. In particular, the Love wave, the Rayleigh wave dispersion and the HVSR curves are innovatively combined in a joint inversion procedure carried out following a global search approach (i.e., the Genetic Algorithm).The procedure is tested using a data set of seismic noise recordings collected at the Bevagna (Italy) test-site. The results of the novel inversion scheme are compared with the inversion scheme proposed by Parolai et al. (2005), where only Rayleigh wave dispersion and HVSR curves are used, and with a cross-hole survey.     

基于矢量波数变换法(VWTM)的多道Rayleigh波分析方法

在近二十年来,多道面波分析法(MASW)由于其便捷、高效等特性在浅层地震勘探领域得到了广泛的应用.本文基于多道面波勘探的采集方式,提出了一种新的面波多道分析方法——矢量波数变换法(VWTM).该方法通过对震源的近似,基于水平层状模型得到台站与震源间近似格林函数,然后进行矢量波数变换得到含有高阶模态Rayleigh波(频率-相速度)频散能量图.本研究首先利用合成地震数据到频散能量图与理论频散曲线进行叠加分析该方法的有效性和正确性;然后与相移法进行对比分析,我们发现在频散能量图中VWTM法对基阶、高阶模态成像均具有更高的分辨率和成像质量;最后我们将其应用于实际多道瞬态面波探测中,通过与相移法进行对比分析,发现VWTM法是一种方便、实用、有效的Rayleigh波频散提取方法.VWTM法提取多模态的Rayleigh波频散特征具有巨大潜力,可为基阶、高阶面波频散联合反演提供丰富的高阶模态频散信息.     

CROSS-BOREHOLE OBSERVATION OF MODE CONVERSION FROM BOREHOLE STONELEY WAVES TO CHANNEL WAVES AT A COAL LAYER1

Conversion of borehole Stoneley waves to channel waves was observed in data from a seismic cross-borehole experiment conducted between wellbores penetrating a thin coal layer at 2022 m depth, near Rifle, Colorado. Traveltime moveout observations show that borehole Stoneley waves underwent partial conversion to channel waves at the coal layer. The channel waves were detected directly in an adjacent borehole 35 m away at receiver positions within the coal. Stoneley waves, subsequently produced by partial conversion of channel waves, were also detected at receiver positions located up to 50 m above and below the coal layer in the adjacent borehole. We infer the channel wave to be the first-higher Rayleigh mode by comparing the observed group velocity with theoretically derived dispersion curves. Identifying the conversion between borehole and stratigraphically guided waves is significant because coal penetrated by multiple wells may be detected without placing a transmitter or receiver in the coal itself.     

RESEARCH ON TIME-DISTANCE INTERPRETATION OF MULTIPLES AND GHOSTS BASED ON MARINE SEISMIC DATA

Currently, the study on the active fault in the land areas is relatively mature, while there is still lack of detection and research on active faults in the sea areas. Marine exploration, which is different from land areas, has a prominent problem due to the existence of strong reflecting interfaces such as water surface and seafloor in the sea, thus the recording is often accompanied by interference of multiples on seafloor reflections. In addition, because of the characteristics of marine seismic exploration, the source exciting in the water and the geophone receiving in the water, ghost wave usually can be recorded simultaneously during the reflected wave propagation. This phenomenon makes it difficult to distinguish the effective waves and the noise, and has always plagued the data and seriously affects the quality of records. In the offshore and other regions of complex structures, such as inclined interfaces, it is difficult to eliminate the interference of multiples accurately by traditional multiples suppression methods, which are based on the horizontal interface assumption. This paper combines the sea area seismic data and its acquisition method, uses simplified model to simulate the multiples based on the time-distance analysis of multiples and their ghost wave in inclined interface. The time-distance characteristics of the multiples and their ghost waves from different interfaces(including the inclined interface)are obtained, and they are consistent with the actual records. The multiples time-distance simulation can help to distinguish the causes of reflected waves, summarize the multiple-wave time-distance characteristics from different interfaces(including inclined interfaces), and analyze the relationship between the characteristics of multiple waves and primary waves. In particular, this simulation has a significant effect on characterizing the internal multiples that are difficult to identify due to inconspicuous periodicity and the multiples of the inclined interface which present the phenomenon that the vertex of the time-distance curve is shifted. On this basis, relying on the time-distance analysis of ghost wave, we analyze the travel time difference characteristics between reflected waves and their accompanying ghost waves. The differences of the travel time characteristics of different orders ghost wave and reflected wave are summarized and the symmetry of the travel time difference between inclined interface and horizontal interface of ghost waves and reflected waves is analyzed. We simulate the distraction of the ghost wave event with the event of the reflected wave and analyze the influence of the ghost wave on the sea area seismic records. These results can improve the practical interpretation of seismic data. At last, the time-distance information is used to synthesize sea area seismic records, which can help us carry out the effective data processing and understand the characteristics of the time-distance and velocity of multiples in different interfaces and the layer artifact caused by multiples. This study combines the time-distance simulation of multiples and their ghost wave with conventional seismic data processing to analyze the pre-stack and post-stack features of multiple waves and their ghost waves in the seismic records of the sea area. The results of this study are conducive to the effective identification of multiples in seismic records in the sea, provide a theoretical basis for multi-wave suppression and prediction, and may facilitate the future study of sea-area seismic activity detection.     

Tomographic inversion of near-surface <Emphasis Type="Italic">Q</Emphasis> factor by combining surface and cross-hole seismic surveys

The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer. We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation. In this study, we drilled number of receiver holes around the source hole, each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys. We also propose a novel tomographic inversion of the Q factor without the effect of the source signature, and examine its stability and reliability using synthetic data. We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield. The results show that seismic absorption in the near-surface layer is much greater than that in the subsurface strata. Thus, it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption. In addition, we derive different Q factors from two frequency bands, which can be treated, to some extent, as evidence of a frequency-dependent Q.     

Seismic on floating ice: data acquisition versus flexural wave noise

Geophysical surveying of the Arctic will become increasingly important in future prospecting and monitoring of the terrestrial and adjacent areas in this hemisphere. Seismic data acquired on floating ice are hampered with extensive noise due to ice vibrations related to highly dispersive ice flexural waves generated by the seismic source. Several experiments have been conducted on floating ice in van Mijenfjorden in Svalbard in the Norwegian Arctic to specifically analyse the extent of flexural waves recorded with various seismic receivers and sources deployed both on top of ice and in the water below. The data show that flexural waves are severely damped at 5 m or deeper below the ice and hydrophone data suffer less from these vibrations compared with data recorded on the ice. Aliasing of single receiver hydrophone data can to some extent be suppressed using an in-line line source of detonating cord. Experiments on ice on shallow water show prominent guided wave modes often referred to as Scholte waves propagating along the seabed. In this case, both flexural and Scholte waves interfere and make a complicated pattern of coherent noise. On shallow water, the positioning and type of the seismic source must be evaluated with respect to the coherent noise generated by these waves. Geophone strings of 25 m effectively suppress both flexural and Scholte waves due to their relative short wavelengths. An airgun generates relative more low-frequency energy than a surface source of detonating cord. Accordingly, seismic mapping of deep seismic horizons seem to be best achieved using geophone strings of such length and an airgun source. For shallow targets, the use of hydrophones in combination with detonating cord is an appropriate solution. Seismic surveying in the Arctic always have to follow environmental restrictions of not disturbing or harming wildlife and not causing permanent footprints into the vulnerable tundra, which implies that the choice of seismic acquisition strategy might occur as a trade-off between optimum data quality and environmental constraints.     

迭部5.9级地震爆破地震波动力学特性的变化

本文利用遥测传输台网可见记录和磁带记录对厂坝炮爆破进行地震波动力学特性分析,发现各次爆破记录有一致性。对失真信号进行滤波处理,经凝聚函数分析证明了爆源的一致性。对1987年1月8日迭部Ms=5.9地震的分析,表明由于介质状态的改变,记录到的地震波优势频率向高频方向移动。在近震中有关台站上,地震波位移谱出现多峰与能量分散现象,拐角频率也向高频方向移动,振幅比呈现低—高—发震的形态。在地震前二十几天,主频处地震波衰减小,而震后明显变大。上述特征可能与该地震孕育过程中应力水平和介质状态变化有关。     

充液井孔中的相控线阵声波辐射器在地层中产生的声场

应用有限元法研究了充液井孔中的相控线阵声波辐射器在井外地层中产生的声场,探讨了增强声源向地层的某一特定方向辐射的声波能量的方法. 采用相控线阵声波辐射器并调节其参数可以控制由充液井孔向地层中辐射的声波能量的方向,使声源向地层中的某个特定方向辐射的纵波能量得到明显提高. 在井间地震勘探和反射声波成像测井等领域中应用相控线阵声波辐射器,将有利于增大探测深度、提高探测分辨率和接收信号的信噪比.     

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.     

Using an airgun array in a land reservoir as the seismic source for seismotectonic studies in northern China: experiments and preliminary results

This paper reports the field setup and preliminary results of experiments utilizing an airgun array in a reservoir in north China for a seismotectonic study. Commonly used in offshore petroleum resource exploration, the airgun source was found to be more useful than a traditional explosive source for large‐scale and long offset land seismic surveys. The airgun array, formed by four 1,500 in³ airguns (a total of 6,000 in³ in volume) was placed at a depth of 6–9 m into the reservoir to generate the pressure impulse. No direct evidence was found that the airgun source adversely affected the fish in the reservoir. The peak ground acceleration recorded on the top of the reservoir dam 100 m away was 17.8 gal in the horizontal direction; this is much less than the designed earthquake‐resistance threshold of 125 gal for this dam. The energy for one shot of this airgun array is about 6.68 MJ, equivalent to firing a 1.7 kg explosive. The seismic waves generated by the airgun source were recorded by receivers of the regional seismic networks and a temporary wide‐angle reflection and refraction profile formed by 100 short‐period seismometers with the maximum source‐receiver offset of 206 km. The seismic wave signature at these long‐offset stations is equivalent to that generated by a traditional blast source in a borehole with a 1,000–2,000 kg explosive. Preliminary results showed clear seismic phases from refractions from the multi‐layer crustal structures in the north China region. Forward modelling using numerical simulation confirms that the seismic arrivals are indeed from lower crustal interfaces. The airgun source is efficient, economical, environmentally friendly and suitable for being used in urbanized areas. It has many advantages over an explosive source for seismotectonic studies such as the high repeatability that is supreme for stacking to improve signal qualities. The disadvantage is that the source is limited to existing lakes or reservoirs, which may restrict experimental geometry.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.