一种新型海底沉积物声学原位测量系统的研制及应用

为了能够精确地测量海底表层沉积物的声学参数,自主研制了一种新型海底沉积物声学原位测量系统,与国内外传统的声学原位测量系统相比,该系统能够实时显示声波波形,调整测量参数,其工作方式除了站位式测量之外,还实现了拖行式连续测量,极大地提高了工作效率.根据前期海试情况,对海底仪器结构进行了重新设计,使之可以同时测量海底沉积物及海底海水的声学参数,同时建立了双向数字信道,解决了测量过程中系统信号的干扰问题.该系统的结构分为两部分:甲板控制单元和水下测量单元,整套系统通过主机控制程序进行控制,采用GPS定位系统测定仪器的大地坐标.为了检验系统的稳定性及准确性,分别进行了实验室水槽实验和海上试验.利用水声测量设备对测量系统进行实验室水槽标定分析,实验结果表明系统测量值相对误差仅为0.04%,测量结果具有较高的精度.海上试验在青岛胶州湾和东海海域进行,获得了试验区域海底沉积物声速和声衰减系数的测量数据,将测量数据与他人的研究结果进行对比分析,结果表明测量数据与前人研究结果一致,较为准确.该原位测量系统在站位式测量和拖行式测量中都能够快速准确地测量出沉积物声速和声衰减系数,可以作为海底底质声学测量的调查设备.     

Retrieval of super‐virtual refraction by cross‐correlation

Topography and severe variations of near‐surface layers lead to travel‐time perturbations for the events in seismic exploration. Usually, these perturbations could be estimated and eliminated by refraction technology. The virtual refraction method is a relatively new technique for retrieval of refraction information from seismic records contaminated by noise. Based on the virtual refraction, this paper proposes super‐virtual refraction interferometry by cross‐correlation to retrieve refraction wavefields by summing the cross‐correlation of raw refraction wavefields and virtual refraction wavefields over all receivers located outside the retrieved source and receiver pair. This method can enhance refraction signal gradually as the source–receiver offset decreases. For further enhancement of refracted waves, a scheme of hybrid virtual refraction wavefields is applied by stacking of correlation‐type and convolution‐type super‐virtual refractions. Our new method does not need any information about the near‐surface velocity model, which can solve the problem of directly unmeasured virtual refraction energy from the virtual source at the surface, and extend the acquisition aperture to its maximum extent in raw seismic records. It can also reduce random noise influence in raw seismic records effectively and improve refracted waves’ signal‐to‐noise ratio by a factor proportional to the square root of the number of receivers positioned at stationary‐phase points, based on the improvement of virtual refraction's signal‐to‐noise ratio. Using results from synthetic and field data, we show that our new method is effective to retrieve refraction information from raw seismic records and improve the accuracy of first‐arrival picks.     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.     

非对称圆弧片状声源在井旁地层中产生的声场

提出了一种可用于随钻三维反射声波成像测井的圆弧片状声源,推导了该声源声学性质在波数-频率域内的数学描述,并利用实轴积分的方法对该声源在井旁地层中产生的声场进行求解,考察了声源的线度、频率等因素对该声源向地层中辐射的声场的影响.研究结果表明,特定尺寸和特定频率下的圆弧片状声源产生的纵波场的水平和垂直指向性图中仅具有一个明显的主瓣,且主瓣三分贝角宽较窄,方位分辨率较高,指向性良好,适用于随钻三维反射成像测井;SV波场和SH波场的指向性图中存在两个或者多个角瓣,SV波场在垂直于井轴方向上辐射的能量为0,另外在井中发射和接收的反射SH波场和SV波场会相互干扰,所以难以利用SV波场和SH波场进行三维反射声波成像;声源的线度和频率等因素对其辐射声场影响较大,考虑到激励效率和方位分辨率等因素,对于本文描述的井孔模型,选择圆周角在75°和90°之间, 主频为12 kHz左右的声源是合适的.     

标本非均匀性对岩石变形声发射时空分布的影响及其地震学意义

基于一组实验结果,讨论了岩石标本的非均匀性（预存微裂纹和宏观裂纹）对声发射 时空分布的影响. 对比分析表明,预存微裂纹的存在使得破裂成核前声发射率快速增加、 b值表现出一种短期异常现象,即在下降背景上出现起伏,从而增加了破裂时间的可预报 性. 宏观结构（节理、层面等）的存在对声发射率和b值的影响与微裂纹相同,而且宏 观结构对声发射的空间分布具有控制作用,声发射丛集的宏观构造通常控制着未来的主破裂 . 这意味着宏观构造的存在使得主破裂的时间和位置预报成为可能.      

Assessing the Validity of Seismo-Acoustic Predictor Equations for Obtaining Seabed and Sub-Surface Sediment Physical Properties

The interdependence between the seismo-acoustic properties of a marine sediment and its geotechnical/physical parameters has been known for many years, and it has been postulated that this should allow the extraction of geotechnical information from seismic data. Though in the literature many correlations have been published for the surficial layer, there is a lack of information for greater sediment depths. In this article, a desktop study on a synthetic seafloor model illustrates how the application of published near-surface prediction equations to subsurface sediments (up to several tens of meters burial depth) can lead to spurious predictions. To test this further, acoustic and geotechnical properties were measured on a number of sediment core samples, some of which were subjected to loading in acoustically-equipped consolidation cells (oedometers) to simulate greater burial depth conditions. For low effective pressures (representing small burial depths extending to around 10 meters subsurface), the general applicability of established relationships was confirmed: the prediction of porosity, bulk density, and mean grain size from acoustic velocity and impedance appears generally possible for the investigated sedimentary environments. As effective pressure increases through, the observed relationships deviate more and more from the established ones for the near-surface area. For the samples tested in this study, in some instances increasing pressure even resulted in decreasing velocities. There are several possible explanations for this abnormal behavior, including the presence of gas, overconsolidation, or bimodal grain size distribution. The results indicate that an appropriate depth correction must be introduced into the published prediction equations in order to obtain reliable estimates of physical sediment properties for greater subsurface depths.     

Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

Seismic characterization of switching platform geometries and dominant carbonate producers (Miocene,Las Negras,Spain)

The primary goals of seismic interpretation and quantification are to understand and define reservoir architecture and the distribution of petrophysical properties. Since seismic interpretation is associated with major uncertainties, outcrop analogues are used to support and improve the resulting conceptual models. In this study, the Miocene carbonates of Cerro de la Molata (Las Negras, south‐east Spain) have been selected as an outcrop analogue. The heterogeneous carbonate rocks of the Cerro de la Molata Platform were formed by a variety of carbonate‐producing factories, resulting in various platform morphologies and a wide range of physical properties. Based on textural (thin sections) and petrophysical (porosity, density, carbonate content and acoustic properties) analyses of the sediments, eleven individual facies types were determined. The data were used to produce synthetic seismic profiles of the outcrop. The profiles demonstrate that the spatial distribution of the facies and the linked petrophysical properties are of key importance in the appearance of the synthetic seismic sections. They reveal that carbonate factory and facies‐specific reflection patterns are determined by porosity contrasts, diagenetic modifications and the input of non‐carbonate sediment. The reflectors of the seismograms created with high‐frequency wavelets are coherent with the spatial distribution of the predefined facies within the depositional sequences. The synthetic seismograms resulting from convolution with lower frequency wavelets do not show these details – the major reflectors coincide with: (i) the boundary between the volcanic basement and the overlying carbonates; (ii) the platform geometries related to changes in carbonate factories, thus sequence boundaries; and (iii) diagenetic zones. Changes in seismic response related to diagenesis, switching carbonate producers and linked platform geometries are important findings that need to be considered when interpreting seismic data sets.     

Acoustic full waveform inversion of synthetic land and marine data in the Laplace domain

Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain.