煤矿三维地震勘探技术发展趋势

我国煤田三维地震勘探经过20余年的发展,已成为当前煤田地质勘探的主流技术,随着煤矿采掘现代化及安全保障系统对勘探精度要求的提高,三维地震勘探技术也须与之相适应,不断地探索与发展。参考近年来相关研究成果与勘探实例,对可在煤田地质勘探中大面积推广应用的三维地震勘探技术进行了介绍,其中包括：高密度空间采样、数字检波器、高精度静校正技术、叠前时间（深度）偏移技术、分频解释技术、波阻抗反演技术、相干体技术,同时也简要介绍了在煤田勘探方面有发展前景的多分量及各向异性技术。     

Gas hydrate saturation from acoustic impedance and resistivity logs in the Shenhu area, South China Sea

During the China’s first gas hydrate drilling expedition -1 (GMGS-1), gas hydrate was discovered in layers ranging from 10 to 25 m above the base of gas hydrate stability zone in the Shenhu area, South China Sea. Water chemistry, electrical resistivity logs, and acoustic impedance were used to estimate gas hydrate saturations. Gas hydrate saturations estimated from the chloride concentrations range from 0 to 43% of the pore space. The higher gas hydrate saturations were present in the depth from 152 to 177 m at site SH7 and from 190 to 225 m at site SH2, respectively. Gas hydrate saturations estimated from the resistivity using Archie equation have similar trends to those from chloride concentrations. To examine the variability of gas hydrate saturations away from the wells, acoustic impedances calculated from the 3 D seismic data using constrained sparse inversion method were used. Well logs acquired at site SH7 were incorporated into the inversion by establishing a relation between the water-filled porosity, calculated using gas hydrate saturations estimated from the resistivity logs, and the acoustic impedance, calculated from density and velocity logs. Gas hydrate saturations estimated from acoustic impedance of seismic data are ∼10-23% of the pore space and are comparable to those estimated from the well logs. The uncertainties in estimated gas hydrate saturations from seismic acoustic impedances were mainly from uncertainties associated with inverted acoustic impedance, the empirical relation between the water-filled porosities and acoustic impedances, and assumed background resistivity.     

Application of constrained sparse spike inversion in reservoir predication:a case study of Oriente Basin in South America

In the conditions of low Signal-to-Noise Ratio(SNR) of seismic data and a small quality of log information,the consequences of seismic interpretation through the impedance inversion of seismic data could be more precise. Constrained sparse spike inversion(CSSI) has advantage in oil and gas reservoir predication because it does not rely on the original model. By analyzing the specific algorithm of CSSI,the accuracy of inversion is controlled. Oriente Basin in South America has the low amplitude in geological...     

三维电磁阵列剖面法的基本原理及应用

本文论述了三维电磁阵列剖面法的基本原理，并从理论上论述了大地电磁测深中的静态效应是由电磁响应中与频率无关的静态项引起的，基于这一原理，本文提出了三维ＥＭＡＰ野外采集方式和数据处理方程。三维ＥＭＡＰ法不仅能消除ＭＴ中的静态效应，而且由于采集信息量大大增加，较常规ＭＴ有更好的地质效果，文中给出的实例证明了这一点。     

Dipping-interface Mapping Using Mode-separated Rayleigh Waves

Multichannel analysis of surface waves (MASW) method is a non-invasive geophysical technique that uses the dispersive characteristic of Rayleigh waves to estimate a vertical shear (S)-wave velocity profile. A pseudo-2D S-wave velocity section is constructed by aligning 1D S-wave velocity profiles at the midpoint of each receiver spread that are contoured using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. Based on the assumption that a dipping-layer model can be regarded as stepped flat layers, high-resolution linear Radon transform (LRT) has been proposed to image Rayleigh-wave dispersive energy and separate modes of Rayleigh waves from a multichannel record. With the mode-separation technique, therefore, a dispersion curve that possesses satisfactory accuracy can be calculated using a pair of consecutive traces within a mode-separated shot gather. In this study, using synthetic models containing a dipping layer with a slope of 5, 10, 15, 20, or 30 degrees and a real-world example, we assess the ability of using high-resolution LRT to image and separate fundamental-mode Rayleigh waves from raw surface-wave data and accuracy of dispersion curves generated by a pair of consecutive traces within a mode-separated shot gather. Results of synthetic and real-world examples demonstrate that a dipping interface with a slope smaller than 15 degrees can be successfully mapped by separated fundamental waves using high-resolution LRT.     

Simple formulas for the dynamic stiffness of pile groups

Simple formulas are derived for the dynamic stiffness of pile group foundations subjected to horizontal and rocking dynamic loads. The formulations are based on the construction of a general model of impedance matrices as the condensation of matrices of mass, damping, and stiffness, and on the identification of the values of these matrices on an extensive database of numerical experiments computed using coupled finite element–boundary element models. The formulations obtained can be readily used for the design of both floating piles on homogeneous half‐space and end‐bearing piles and are applicable for a wide range of mechanical and geometrical parameters of the soil and piles, in particular for large pile groups. For the seismic design of a building, the use of the simple formulas rather than a full computational model is shown to induce little error on the evaluation of the response spectra and time histories. Copyright © 2009 John Wiley & Sons, Ltd.     

川西地区壳幔结构与汶川MS8.0级地震的孕震背景

收集了四川数字地震台网记录的57个远震事件,并从宽频带数字化三分量地震记录中计算出了马尔康（MEK）、都江堰（YZP）、中江（JJS）、江油（ZJG）、广元（YTS）、康定（GDS）、汉源（XJP）、雅安（MDS）、峨眉山（EMS）、沐川（WMP）、仁寿（YGD）、荣县（HMS）等12个台站下方的远震P波径向接收函数.另外,引入地震勘探中的动校正技术,将各台的接收函数校正到67°的参考震中距处,然后对接收函数进行叠加以增强信号,并把叠加接收函数作为台站下方的平均接收函数.最后,利用台站下的平均接收函数反演得到S波速度结构.反演结果表明：以锦屏山—龙门山断裂为界,其西侧地壳厚达70 km,而东侧仅为50 km左右,Moho面在断裂下方形成了一个陡坎;川西地区的地壳速度结构与川中地区差异较大,主要表现在川西地区的中地壳存在厚度为8~22 km的低速层.在都江堰、雅安一带,低速层的厚度最大,其厚度在20~22 km之间,其上地壳为一个坚硬固体,在区域构造应力场作用下,形成了孕育大地震的构造环境.     

Probabilistic impedance of foundation: Impact of the seismic design on uncertain soils

In the linear seismic design of buildings, the (deterministic) substructure method is a customary and efficient approach. However, the existence of spatial variability in the parameters of the mechanical model of the soil, as well as parametric errors, calls for the use of probabilistic approaches in order to provide a reliable design of the structure. The construction of probabilistic models of the soil impedance matrix provides a natural path to such approaches within the context of the substructure method. Two main techniques are described in this paper: a parametric one, typically using the stochastic finite element method, and a nonparametric one, which was introduced more recently. The latter is explored more specifically, and the possibilities it offers in terms of seismic design are presented. In particular, it is shown that it allows for the estimation of quantiles of the quantities of interest, rather than confidence intervals, which lead to highly conservative design. Copyright © 2008 John Wiley & Sons, Ltd.     

Identifying causes of ground-penetrating radar reflections using time-domain reflectometry and sedimentological analyses

Ground-penetrating radar (GPR) is a geophysical technique widely used to study the shallow subsurface and identify various sediment features that reflect electromagnetic waves. However, little is known about the exact cause of GPR reflections because few studies have coupled wave theory to petrophysical data. In this study, a 100- and 200-MHz GPR survey was conducted on aeolian deposits in a quarry. Time-domain reflectometry (TDR) was used to obtain detailed information on the product of relative permittivity (ɛ_r) and relative magnetic permeability (μ_r), which mainly controls the GPR contrast parameter in the subsurface. Combining TDR data and lacquer peels from the quarry wall allowed the identification of various relationships between sediment characteristics and ɛ_rμ_r. Synthetic radar traces, constructed using the TDR logs and sedimentological data from the lacquer peels, were compared with the actual GPR sections. Numerous peaks in ɛ_rμ_r, which are superimposed on a baseline value of 4 for dry sand, are caused by potential GPR reflectors. These increases in ɛ_rμ_r coincide with the presence of either organic material, having a higher water content and relative permittivity than the surrounding sediment, or iron oxide bands, enhancing relative magnetic permeability and causing water to stagnate on top of them. Sedimentary structures, as reflected in textural change, only result in possible GPR reflections when the volumetric water content exceeds 0·055. The synthetic radar traces provide an improved insight into the behaviour of radar waves and show that GPR results may be ambiguous because of multiples and interference.