复杂地表边界元-体积元波动方程数值模拟

复杂近地表引起来自深部构造的地震反射信号振幅和相位的异常变化,是影响复杂近地表地区地震资料品质的主要原因.本文采用边界元-体积元方法,通过求解含复杂地表的波动积分方程,来模拟地震波在复杂近地表构造中的传播.其中,边界元法模拟地形起伏和表层地质结构对地震波传播的影响;体积元法模拟起伏地表下非均质低降速层的影响.与其他数值...     

自相关法单频干扰识别与消除方法

常规单频干扰压制方法需要估算频率、振幅和相位,由于频率和相位与单频干扰呈非线性函数关系,因此任何参数估算方法都会非常费时,且计算效率低下.本文提出的自相关法单频干扰识别与消除方法,利用地震数据自相关函数计算单频干扰自相关函数,采用余弦函数自适应减算法快速估算单频干扰.本方法不需要估算单频干扰的频率、相位和振幅参数,因此...     

Refractor velocity analysis: a signal processing procedure

This paper presents a signal processing procedure to perform refractor velocity analysis. The procedure enables one to obtain the seismic velocity from the refracted wavefield without the picking of refracted arrival times. Two processing procedures are derived, one starting from a seismic interferometric approach and another, from the conventional reciprocal method and generalized reciprocal method approaches. The theoretical equivalence of the two approaches is also demonstrated. The proposed processing procedure is applied to synthetic data in order to test the influence of some procedural parameters and its capability to reconstruct a known velocity model starting from refracted signals, without and with perturbations, in arrival times and noise; finally, it is applied to a field data set.     

Seismic interferometry experiment in a shallow cased borehole using a seismic vibrator source‡

We present the results of a seismic interferometry experiment in a shallow cased borehole. The experiment is an initial study for subsequent borehole seismic surveys in an instrumented well site, where we plan to test other surface/borehole seismic techniques. The purpose of this application is to improve the knowledge of the reflectivity sequence and to verify the potential of the seismic interferometry approach to retrieve high‐frequency signals in the single well geometry, overcoming the loss and attenuation effects introduced by the overburden. We used a walkaway vertical seismic profile (VSP) geometry with a seismic vibrator to generate polarized vertical and horizontal components along a surface seismic line and an array of 3C geophones cemented outside the casing. The recorded traces are processed to obtain virtual sources in the borehole and to simulate single‐well gathers with a variable source‐receiver offset in the vertical array. We compare the results obtained by processing the field data with synthetic signals calculated by numerical simulation and analyse the signal bandwidth and amplitude versus offset to evaluate near‐field effects in the virtual signals. The application provides direct and reflected signals with improved bandwidth after vibrator signal deconvolution. Clear reflections are detected in the virtual seismic sections in agreement with the geology and other surface and borehole seismic data recorded with conventional seismic exploration techniques.     

Repeatability analysis of land time‐lapse seismic data: CO2CRC Otway pilot project case study

Time‐lapse seismics is the methodology of choice for remotely monitoring changes in oil/gas reservoir depletion, reservoir stimulation or CO₂ sequestration, due to good sensitivity and resolving power at depths up to several kilometres. This method is now routinely applied offshore, however, the use of time‐lapse methodology onshore is relatively rare. The main reason for this is the relatively high cost of commercial seismic acquisition on land. A widespread belief of a relatively poor repeatability of land seismic data prevents rapid growth in the number of land time‐lapse surveys. Considering that CO₂ sequestration on land is becoming a necessity, there is a great need to evaluate the feasibility of time‐lapse seismics for monitoring. Therefore, an understanding of the factors influencing repeatability of land seismics and evaluating limitations of the method is crucially important for its application in many CO₂ sequestration projects. We analyse several repeated 2D and 3D surveys acquired within the Otway CO₂ sequestration pilot project (operated by the Cooperative Research Centre for Greenhouse Technologies, CO2CRC) in Australia, in order to determine the principal limitations of land time‐lapse seismic repeatability and investigate the influence of the main factors affecting it. Our findings are that the intrinsic signal‐to‐noise ratio (S/N, signal to coherent and background noise levels) and the normalized‐root‐mean‐square (NRMS) difference are controlled by the source strength and source type. However, the post‐stack S/N ratio and corresponding NRMS residuals are controlled mainly by the data fold. For very high‐fold data, the source strength and source type are less critical.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

多重散射问题的改进算法以及双排非连续弹性屏障对水平向剪切波的隔离研究

与Aviles等学者的多重散射隔振理论相比较,采用更为完整的Fourier-Bessel函数系的级数展开表达式来研究弹性波入射下的多重散射问题。引入Graf加法定理,并结合位移和应力边界条件,推导了无限均匀弹性介质中任意分布的多个弹性圆柱体对于平面水平向剪切波(SH波)入射时散射系数的理论解析解。随后,取圆柱体分布于两条相互平行的直线上,则问题演变为双排弹性桩屏障对于弹性波的隔离问题。引入无量纲位移(屏障后土体总的位移与仅由入射平面水平向剪切波引起的位移的比值)和透射系数(屏障宽度范围内无量纲位移的几何平均值)的概念来研究屏障的隔振效果。重点讨论了双排桩间距h以及桩土的剪切模量比ξ对于屏障整体隔振效果的影响,阐述了双排桩屏障的隔振性状。特别地,当h=0时,该问题退化为常见的单排非连续屏障的隔振问题。     

龙门山中段彭灌断裂带汶川Ms8.0地震同震破裂断层识别及其浅、表构造特征

2008年汶川Ms8.0地震在龙门山中段的彭灌断裂带产生的地表破裂,是该地震产生的第二大地表破裂带.综合应用地质、钻井以及二维、三维地震数据,利用横贯前山带的多条人工地震反射剖面,对彭灌断裂带产生同震破裂的断层进行准确识别和解释.研究表明,龙门山中段的彭灌断裂带是一套由3条主要断层和次级广泛发育的断裂组合构成,浅层表现...     

海洋微生物偏振光散射的双组分模型

光学散射方法被广泛用于获取海洋微生物的形态和浓度信息.通过测量散射光的空间分布可以测量海洋微生物的形态和浓度;而最近发现,散射光的偏振特征却对细胞内部亚微米级散射颗粒更加敏感.把海洋微生物等效为粒径远大于入射波长的“米氏”粒子和远小于入射波长的“瑞利”粒子的混合体,以两种不同粒径的实心球作为微生物偏振光散射模型.利用蒙...     

琼东南盆地长昌凹陷三维地震区中新统微观构造成因探讨

基于高精度三维地震资料的精细解释和方差体顺层切片技术,在琼东南盆地长昌凹陷中新统首次发现高密度小断裂,它们近SN走向,断距6～28m,长0.5～4km,倾角约60°.这些小断裂既不是渐新统继承性断裂,也不是近期热点讨论的多边形断层,而是在中新世特定地质条件下由局部引张形成的微观构造.其成因为:1)中新世琼东南盆地处于坳...