凤凰山地区蒲江-新津-德阳隐伏断裂反射波法地震勘探与断裂活动性研究

利用反射波法地震勘探确定了在凤凰山地区蒲江-新津-德阳隐伏断裂通过的具体位置以及断裂的规模和性质,证明了该断裂在向东俯冲的过程中错切了晚更新统沉积物,具有明显的第四纪活动性,对凤凰山台地的形成起到了一定的控制作用。地面地质调查及历史地震研究也证明蒲江-新津-德阳隐伏断裂具有晚更新世活动性,且断裂南段的活动性明显强于北段,表明该断裂的北段是成都平原内部的一条弱地震活动带。同时,文中还对在城市近郊人烟繁杂地带反射波法地震勘探中野外观测系统的布置及地震参数的选取进行了探索,小道间距、小偏移距、多接收道、短排列、多次覆盖、高频检波器接收的工作方法和野外地震参数的合理选取对实现主要干扰波和反射波信号的有效分离,提高地震勘探分辨率起到了至关重要的作用     

喀斯特地区煤田多孔联采3D-RVSP勘探实践:以中国五轮山为例（英文）

喀斯特地貌等复杂地表地区的常规地震勘探一直面临着激发和接收条件差的难题。RVSP (Reverse Vertical Seismic Profile)地震勘探采用井中激发、地面接收的观测方式,能够一定程度上减少复杂地表对地震波传播的影响。本文通过数值模拟和实际资料分析指出:在岩溶发育的喀斯特地貌等复杂地表地区,浅层激发的地震波受各种多次波、面波和岩溶地质体的散射干扰作用,反射地震数据品质极低,且难以通过去噪方法提高信噪比。当加大激发深度,避开浅表层复杂地质体时,能够提高反射波品质。RVSP方法采用井中激发,正好具有上述潜在优势。根据其观测系统得特殊性,利永有效反射波与大部分干扰波存在的视速度差异进行滤波可进一步提高反射波信噪比。本研究在典型的喀斯特地貌地区开展了8孔联采3D-RVSP勘探研究,结果表明:在岩溶发育的复杂地表条件地区,3D-RVSP观测能够获得较高质量的反射地震资料;采用多孔联采的3D-RVSP方式,具有数据采集效率高、改善地下覆盖均匀性的特点;获得的成像结果分辨率高,能够对细小断裂构造进行精细刻画。     

地震干涉技术被动源地震成像

被动源地震成像是基于地下反射波响应和透射波响应之间的关系,通过在地表接收由地下非人工震源发出的透射波响应,利用互相关运算来合成反射波响应的方法.合成的反射波响应,又称虚炮集记录等价于地表地震剖面记录.本文研究了由地下随机分布的被动源发出的透射波记录来合成虚炮集记录的基本原理,推导了反射波记录和透射波记录的关系公式;并通过数值模拟方法,分析了被动源成像的可行性.结果表明:利用被动源地震信息进行成像是可行的.     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.     

大庆电场地震实验资料的面波特征

通过对电场作用下的地震面波特征的研究，发现电场对地震面波的影响比对地震反射波的影响要大，这对进一步研究电场地震的机制有重要作用。     

基于形态成分分析地震信号二维域面波分离方法研究（英文）

地震数据中的面波是严重降低地震资料信噪比的干扰波,它的存在影响了后续地震资料的处理与解释。本文根据地震记录中面波与反射波信号形态结构的差异,采用基于二维字典形态成分分析方法对面波噪声与反射波进行分离。根据面波信号的低频、低视速度和频散的特性,选择二维非抽样离散小波变换作为面波的稀疏表示字典,根据反射波局部相关性较强的特点,选择二维局部离散余弦变换作为反射波的稀疏表示字典,构建地震记录在联合二维字典下的稀疏表示模型并采用块协调松弛算法进行求解,将地震记录分解为反射波部分和面波部分。对合成地震信号以及实际地震资料的处理结果表明本文方法不仅能有效压制强能量的面波干扰,而且还能很好保护反射波信号的波形。     

南海深水多次波成像

南海海域水深可达5 km,是我国未来油气勘探最有前景的区域.南海海域地震资料含有各种类型多次波,在地震成像前需将多次波去除.然而与传统的一次反射波偏移比较,多次波携带更多的地下结构信息,多次波偏移可提供更好的地下成像.本文提出了一种修改传统逆时偏移的方法,利用多次波进行成像.该方法用包含一次波和多次波的记录代替震源子波,用预测多次波代替一次反射波作为输入数据.在逆时偏移处理流程中,将地表多次波记录沿时间延拓到各个深度层,通过互相关成像条件将包含一次波和多次波的地震记录沿时间延拓到相同深度层.将多次波偏移结果与一次波偏移结果的差应用目标函数进行匹配,弥补一次波偏移的不足.对Sigsbee2B数据做数值测试的结果表明与传统一次反射波偏移相比,该方法可以得到范围更大、均衡性更好的地下照明度.将该方法应用到南海深水实际资料成像,获得了南海深水地震资料多次波成像剖面.偏移结果表明,多次波成像可以弥补反射波照明的不足,由此说明多次波偏移是未来的一个重要研究方向.     

地震数据的反射波动方程最小二乘偏移

基于反射波动方程,本文提出了一种估计地下反射率分布的地震数据最小二乘偏移方法.高频近似下,非齐次的一次反射波动方程的源项是由反射率与入射波场的时间一阶导数相互作用产生的.根据反射波动方程,利用线性最小二乘反演方法由地震反射数据重建出地下产生反射波的反射源,再结合波场正演计算出的地下入射波场,得到地下反射率分布的估计.在地下反射源的线性最小二乘反演重建中,我们采用迭代求解方法,并以地震波的检波器单向地下照明强度作为最小二乘优化问题中Hessian矩阵的近似.     

高密度面波勘探和小型台阵地脉动观测在西部地区地震防灾中的应用

鉴于我国西部地区地震防灾需要一种简便快速、低成本的地下构造调查方法,本文利用对甘肃岷县漳县地震灾害科考机会,对高密度面波勘探和小型台阵地脉动观测在西部山区地下构造调查中的有效性进行了试验.试验选取房屋倒塌和滑坡等不同灾害类型以及不同地质条件的场地进行;2种方法的数据采集都使用普通小型便携式面波仪和4 Hz地震检波器,其中,小型台阵采用底边长16m的正三角形台阵,观测时间约30分钟.试验结果显示,高密度面波勘探能够探测地下约30～40m以内的剪切波速度构造,小型台阵探测深度稍深,约为40～50m,但是浅层分辩率较低.试验还显示,两种勘探方法使用同一台普通面波仪和地震检波器,设备小、成本低、便于携带,非常适合西部地区地震防灾之地下构造调查.     

地震反射波成像新方法及其理论基础

提出了一个地震反射波的成像公式及其理论基础. 波的传播现象应该由瞬时谱来刻画, 在数学上表示成三元(即时间、空间和频率)复函数. 介质的物理参数, 如反射系数波和阻抗等均是空间和频率的复函数. 给出了反射波谱和介质反射系数间的关系, 反演公式中还考虑到多次反射和结构滤波作用. 同时也指出传统的卷积模型和波动方程反演方法中存在的问题. 各种反演公式均可利用小波变换和三基色方法作出地下介质参数分布的剖面图, 不同颜色表示不同介质.     

基于匹配滤波的多次波压制方法研究

常规滤波方法对浅海多次波压制不能取得理想效果.通过改变模型表面边界条件,求解波动方程实现多情景地震数值模拟.利用自由边界条件模拟时获得含海平面强反射多次波地震记录,而利用吸收边界条件时得到理想条件下无多次波地震记录.对比分析多次波对有效地震信号的影响得到海平面强反射多次波改变了有效地震信号的振幅、频率和相位.利用归一化方法计算得到压制多次波滤波器算子,并作用于含海平面强反射多次波地震记录.叠加地震剖面对比分析表明,归一化滤波器算子较好的压制了地震记录中的多次波,提高地震资料分辨率.最后,根据多次波产生条件,给出了利用实际区域简化地质模型的地震数值模拟结果计算归一化滤波器算子的方法,从而利用该方法压制实际地震数据中多次波.     

海上地震资料子波提取鬼波压制技术

海洋地震勘探由于受到海面强反射界面的影响,记录数据中存在鬼波,致使数据频谱中产生陷波效应,限制数据的频带宽度,影响数据的分辨率,干扰地层反射的识别.近年来发展的鬼波压制方法,能够较好的压制地震数据中的鬼波,然而现有鬼波压制技术容易产生延续相位,且计算效率较低.为进一步提高鬼波压制效果和工作效率,本文对海上地震资料特性进行分析,研究认为,海上地震资料中的海底反射信号真实记录了海底有效反射、震源鬼波和检波点鬼波的形态,可以利用提取的海底反射子波对数据进行确定性子波反褶积处理,能够有效的压制震源鬼波和检波点鬼波.模型海上实际数据测试结果表明,该方法具有鬼波压制效果好、计算效率高的优势,处理之后的数据频带得到有效拓宽,低频和高频信息均得到加强.     

THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA

In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research.     

浅层地震勘探方法在活动断层中的应用——以邯郸市活动断层探测为例

浅层地震勘探方法在城市活动断层探测中发挥着重要作用,实践表明,浅层地震反射波法可以获取深度十米至几十米的浅层地层反射信号,且大部分反射剖面均可较清楚地揭示出浅部断层位置和断层特征,地震数据处理是准确识别近地表活动断层的重要环节,也是提高资料分辨率和信噪比的有效途径。本文应用邯郸市活动断层浅层地震资料,运用多途径、分步骤的去噪技术和方法,压制干扰,极大地提高了资料的分辨率和信噪比,并对活动断层的典型剖面进行重点研究和验证。     

浅层地震反射资料的多阶振型面波反演

浅层地震反射波法和面波方法是两种相互独立发展的地震勘探方法,在各自的数据采集和处理中,对方都是作为干扰信号而存在. 本文利用浅层地震反射资料中被视为干扰的面波信号,通过成熟的多道面波勘探技术处理浅层地震反射资料,在频率-波数域中提取多阶振型面波的频散曲线,并基于该曲线反演浅地表S波速度结构. 这种方法充分开发利用了已有数据,无需单独的面波数据采集系统,同时为解释浅层地震反射资料提供了额外的信息约束. 结果表明:浅层地震反射资料中可提取出可靠的多阶振型面波频散曲线,并能给出稳定的反演结果,同时,面波反演的多解性可以通过高阶振型反演得以进一步约束;低速层的存在是观测频散曲线出现振型跳跃或呈“之”字形回折的必要条件而非充分条件.      

非零偏VSP技术在油田复杂断块开发中的应用

非零VSP是一种地面激发,井中接收的地震观测技术.其信号能量强,频带宽,波场信息丰富.能够获得“三高”的地震反射信息.针对复杂断块油田开发后期的精细描述,非零VSP是一项有效的开发地震技术手段.多年来,在三维地震资料品质较差情况下,VSP是对三维地震信息的一种补充,解决了辽河油田复杂断块开发中存在的小断层、薄砂体认识不清的许多地质问题.为油田复杂断块开发提供了一种有效的技术手段.     

分段光滑曲线边界波动方程数值模拟研究

矩形网格有限差分法在地震波传播数值模拟方面具有计算速度快的显著优势,但该方法在处理复杂边界问题上存在着效率低的严重缺陷.本文针对分段光滑曲线边界定义了尖点处的一种正则导数,给出了矩形网格情形分段光滑曲线网格边界点法向导数的一种插值计算方法.采用矩形网格有限差分法对复杂边界地球介质模型进行地震波场数值模拟,并采用波场系列快照技术揭示地震波在起伏地表和复杂介质中的传播规律.模拟结果表明：法向导数插值计算方法为矩形网格有限差分法处理复杂边界提供了有效途径,采用波场系列快照技术可以清晰地展现地震波在反射界面的反射和透射规律、在尖点的绕射规律以及在自由表面的直达波和多次反射规律.     

基于波动方程的上下缆地震数据鬼波压制方法研究

本文发展基于波动方程的上下缆鬼波压制方法,推导了上下缆地震波场频率波数域波动方程延拓合并公式.基于Fourier变换的波场解析延拓确保上下缆资料振幅相位的一致性,消除了长拖缆远偏移距信号的计算误差,同时具有较高的计算效率;上下缆地震波场的波动方程法合并有效解偶鬼波干涉,实现综合利用上下缆地震数据压制鬼波.理论模型数据和实际采集地震数据的测试表明了方法的有效性.     

Radon变换在非规则观测系统中波场分离的应用

Radon变换是数据处理中广泛应用的一种方法技术.本文介绍了Radon变换在规则及其非规则观测系统中波场分离中应用,讨论了Radon变换的公式、实现方法和有关计算参数的选择.通过理论模型的试算,阐明了Radon变换在缺失道、道距变化情况下依然可以高精度的分离波场,并对线性、抛物Radon变换进行了比较,对隧道超前预报反射波资料进行了处理,给出了隧道超前预报处理流程中提取反射波的算例.     

再论地震数据偏移成像

利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高.