海底电缆多次波压制方法研究

多次波一直是影响地震处理与解释的主要因素之一.本文基于表层相关多次波衰减(Surface Related Multiple Elimination,SRME)理论,对其进行一定的改造,使其可以预测出海底电缆(Ocean Bottom Cable,OBC)数据中所有与表层相关的多次波,而后将预测出的多次波自适应减去.本文对OBC波场,OBC多次波的压制方法给出了介绍,对SRME方法压制OBC多次波理论给出了详尽的阐述.利用SRME理论预测OBC多次波完全数据驱动,不需要地下介质的任何信息,理论数据验证了方法的有效性.     

多次波压制方法及优化组合压制技术研究

讨论了4种重要多次波压制方法,即预测反褶积、T-P域预测反褶积、高分辨率Radon滤波以及SRME方法。在Pluto模型数据处理基础上,多角度研究了每一方法压制多次波的适用条件及参数选择问题,以指导资料处理中有针对性组合压制方法。研究表明将SRME、高分辨率Radon滤波和近道内切相结合可有效压制各种类型的多次波。     

改进的多道预测算子压制浅水多次波方法

在浅水情况下,由于观测数据中缺少近偏移距信息,水层多次波的压制面临挑战.利用多道预测算子压制水层多次波是浅水环境下压制多次波的重要方法之一,这种方法先从输入数据中估计出多道预测算子,再将预测算子和输入数据做褶积预测出水层相关多次波.然而,估计的多道预测算子很容易受噪声污染,从而影响多次波模型的精度.所以,我们提出了改进的多道预测算子压制浅水多次波方法.该方法先从数据中估计出多道预测算子,并利用估计的算子构建出精确的水层模型;然后,通过计算算子的走时信息、估计振幅信息、合成新算子三个步骤来修正原始的多道预测算子.修正的算子不仅不受噪声影响,还含有精确的走时信息、可靠的振幅信息;最后,该方法用修正的算子来预测多次波,并结合自适应相减,将预测的多次波从输入数据中去除.通过合成数据和实际资料的验证表明,相比于原始的多道预测算子压制多次波方法,改进的方法能够取得更好的压制效果.

     

基于反馈环的多次波压制方法研究

基于波动方程理论压制多次波,综合考虑了多次波传播的运动学特性和动力学特性,避免了非波动法解决问题所带来的局限性,已成为多次波去噪的发展趋势.首先对多次波压制国内外研究现状做了简单的阐述.在此基础上,给出了反馈环波动方程预测减去法的方法原理.介绍了抛物线拉东(Radon)变换外推近道数据、谱模拟法提取震源子波、均衡伪多道匹配和多次波减去法效果评价等四大关键技术的思路及主要计算公式.用JAVA语言自行开发了相应的模块,并应用射线法多次波模型和SEG/EAGEPluto模型,对以上方法技术进行了检验,去噪效果良好.对海上某多次波特发育实际资料进行了处理,取得了满意的效果.     

多次波压制技术在南海北部陆缘深水区的应用

南海北部洋陆过渡带水深变化剧烈,海底地形崎岖,构造复杂多变,地震勘探中多次波十分发育,常规处理难以取得很好的成像效果.由于方法的局限性,任何一种多次波压制方法都很难彻底地压制全部多次波.本文通过对国家基金委共享航次在南海北部陆缘采集的二维反射地震资料的分析,总结了不同多次波压制技术的优缺点和适用范围,结合滤波类方法和波动理论类方法,综合应用了表面相关多次波衰减法(SRME,Surface Related Multiple Elimination)、τ-p域预测反褶积、高分辨率抛物线Radon变换以及叠后压制残余多次波等方法去除多次波,有效突出了一次反射波,证明了方法的有效性和实用性,同时也形成了针对南海深水资料多次波压制的一种有效的处理流程和方法组合.     

反馈迭代法压制表面多次波效果分析

地震数据处理通常假定反射数据仅由一次波组成,地震勘探资料中的表面多次波通常被视为相干噪音予以压制.基于波动方程预测的反馈迭代法为数据驱动的方法,可有效压制复杂地下介质的表面多次波,通常可由级数展开法和迭代法两种方法实现.文中分别从理论阐述以及应用效果分析方面对基于波动方程预测的反馈迭代法的两种实现途径进行研究,并采用GPU/CPU协同并行加速计算预测表面多次波,改善以往CPU计算效率低的状况.文中对含表面多次波的炮数据和复杂的SMAART模型进行了计算,并将利用级数展开法与迭代法压制表面多次波的效果进行对比分析,结果表明,迭代形式反馈迭代法压制表面多次波的方法效果更佳.由于实现迭代法和级数展开法需要完成的自适应匹配相减的次数不同,迭代法的计算成本略高于级数展开法.     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

平面波域反数据处理压制多次波方法研究

在地震勘探领域,尤其是海洋地震勘探中,多次波一直是影响地震处理与解释的主要因素之一.本文基于x-t域的反馈模型及多次波衰减的理论,详尽推导了x-t域反数据处理的方法;基于平面波域多次波的产生机制,借鉴x-t域反数据处理的方法,推导出在平面波域进行反数据处理的原理.文中给出一个有限差分的模拟数据进行测试,处理结果表明,本方法可以有效地衰减表层相关多次波,提高地震数据分析的精确性,在保护了一次波能量的同时,可以更加有效快捷地去除多次波.     

利用最小平方基准面延拓压制自由表面多次波II: 海洋拖缆观测数据

压制自由表面多次波一直是海洋拖缆地震数据处理中非常关键的环节.浅水等复杂条件下的多次波压制给现有方法构成了严峻的挑战.本文姊妹篇基于卷积型互易定理建立了地下基准面上的虚拟观测信号与实际观测地震数据之间的转化方程,进而借助约束优化反演理论提出最小平方基准面延拓(LSR)方法,用于消除海底观测地震数据中的表面多次波.本文遵循相同的原理,分别针对双检波器同时记录声压场和质点垂直速度信号、单检波器仅记录声压场信号两种拖缆采集方式,探讨基于LSR的表面多次波压制方法及其影响因素.理论模型合成数据实验展示了方法特点与有效性,南黄海浅水拖缆双检地震数据处理结果表明方法具有良好的应用前景.

     

多次波压制的研究现状与进展

目前地震资料处理技术主要利用一次波对地下地质体进行成像或反演,多次波的存在会严重影响地震成像、反演与解释结果.因此,多次波通常作为噪声在地震数据叠前处理阶段进行压制.为了合理选择不同的多次波压制方法,提高地震资料信噪比和保真度,本文系统综述了多次波压制的国内外研究现状,简要介绍了滤波法、预测相减法、稀疏反演法的基本原理,并围绕理论假设条件与实际数据要求,分析比较了各种方法的适用性与优缺点.基于波动方程的预测相减法克服了滤波法关于地震信号统计特征的假设,增强了对复杂地震数据的适应性,但最小二乘自适应相减会损害有效信号.稀疏反演法采用全波形反演的方法估计一次波与多次波,避免了自适应相减步骤,能够更加保真地恢复有效信号.本文通过实际算例重点阐述了稀疏反演法相对于反馈迭代法在多次波压制中的突破性进展,同时指明了稀疏反演法存在的问题并展望了未来的发展方向.     

控制照明与面向目标成像的观测系统设计

常规地震观测系统设计方法基于地下水平层 状介质的假设,通常不能适应复杂构造情况。我们 从控制照明的思想出发,提出了一种面向目标成像 的地震观测系统设计方法,该方法需要一个由初步 地震解释得到的速度模型。利用单程傅立叶有限差 分波场传播算子将目标层的平面源延拓到地表,通 过分析从目标层延拓到地表的波场能量的分布,可 以确定目标层成像所需要的炮点或者检波点的位置。 利用SEG-EAGE盐丘模型数值试算结果表明,该方法 用于设计面向目标成像的特定地震采集系统。     

Equilibrium bank geometry and the width of shallow sandbed streams

Cross-sections of 16 straight sandbed streams in Minnesota, Iowa, and Nebraska were surveyed. Two stratigraphic horizons were found in the banks at each site, an upper cohesive unit usually composed of silt and clay and a lower unit composed of sand. Bank erosion on these rivers occurs when the upper cohesive unit is undercut by scour at bends. The overhanging cohesive block fails by beam or cantilever failure. As upper bank failure is a direct result of undercutting, the stability and rate of retreat of the bank are largely determined by erosion of the sandy part of the bank. The cohesive layer has little influence on bank retreat and width adjustment on the rivers studied here. A quantitative lateral sediment transport model developed by Parker (1978a) is used to calculate the steady-state geometry of the sandy part of the bank. Results are obtained for the shape, length, and height of the sandy part of the bank. The model predicts the length of the bank fairly well, and the theoretical equation for the height of the bank is of the correct form. The model, however, overestimates the slope of the bank. The height of the sandy part of the bank (D_b) is approximately equal to the depth of the mean annual flow. Since D_b is determined by the lateral sediment transport model, the width (W) may be obtained from the equation of continuity (Q = WD_bV), published flow (Q) data, and a resistance equation for the mean velocity, V. The calculated widths are similar to those measured in the field.     

Robust f‐x projection filtering for simultaneous random and erratic seismic noise attenuation

Linear prediction filters are an effective tool for reducing random noise from seismic records. Unfortunately, the ability of prediction filters to enhance seismic records deteriorates when the data are contaminated by erratic noise. Erratic noise in this article designates non‐Gaussian noise that consists of large isolated events with known or unknown distribution. We propose a robust f‐x projection filtering scheme for simultaneous erratic noise and Gaussian random noise attenuation. Instead of adopting the ?₂‐norm, as commonly used in the conventional design of f‐x filters, we utilize the hybrid ‐norm to penalize the energy of the additive noise. The estimation of the prediction error filter and the additive noise sequence are performed in an alternating fashion. First, the additive noise sequence is fixed, and the prediction error filter is estimated via the least‐squares solution of a system of linear equations. Then, the prediction error filter is fixed, and the additive noise sequence is estimated through a cost function containing a hybrid ‐norm that prevents erratic noise to influence the final solution. In other words, we proposed and designed a robust M‐estimate of a special autoregressive moving‐average model in the f‐x domain. Synthetic and field data examples are used to evaluate the performance of the proposed algorithm.     

The attenuated Ricker wavelet basis for seismic trace decomposition and attenuation analysis

The widely used wavelets in the context of the matching pursuit are mostly focused on the time–frequency attributes of seismic traces. We propose a new type of wavelet basis based on the classic Ricker wavelet, where the quality factor Q is introduced. We develop a new scheme for seismic trace decomposition by applying the multi-channel orthogonal matching pursuit based on the proposed wavelet basis. Compared with the decomposition by the Ricker wavelets, the proposed method could use fewer wavelets to represent the seismic signal with fewer iterations. Besides, the quality factor of the subsurface media could be extracted from the decomposition results, and the seismic attenuation could be compensated expediently. We test the availability of the proposed methods on both synthetic seismic record and field post-stack data.     

DINI12电子水准仪在台站短水准测量中的应用

主要介绍了DINI12电子水准仪的特性、原理、检测方法及代替原有光学水准仪在台站短水准测量中的应用,并结合电子水准仪的实际观测数据,对该仪器的使用方法和测量精度进行了探讨。     

地层衰减定量模拟的地震物理模拟方法

Q值是描述地层衰减特征的重要地震物理参数,研究Q值对提高地震资料分辨率,提高地震成像精度具有重要意义.本文从地震物理模拟中模型材料的动力学特征出发,提出了一种地层衰减定量模拟的地震物理模拟方法.通过铝样验证了衍射校正方法的准确性,利用有机玻璃实验分析了不同测量方法求取Q值的测量精度问题,从而优选测量方法;然后基于复合材料的实验数据,建立复合材料配比与各物理参数之间的函数关系,利用该函数关系指导地层衰减定量模拟.结合单层地震物理模型和近地表地震物理模型,研究了模拟地震记录及其频谱特征与野外实际情况之间的相似性,基于建立的函数关系研究分析了地层衰减定量模拟的准确性.实验结果表明,基于地震衰减定量模拟方法构建近地表地震物理模型时,模型中模拟参数不仅与设计参数保持很好的一致性,同时具有很高的准确性,模拟地层的频谱衰减特征与野外地层保持一致.

     

小波尺度域含气储层地震波衰减特征

黏弹性衰减因子Q的可靠估计可通过Q反褶积来提高地震资料的分辨率并有助于振幅分析. 本文从小波理论出发,结合地震波在黏弹性介质中的传播方程,推导出小波尺度域地震波能量衰减公式. 能量衰减公式具有下列性质：（1）Q值越大,能量衰减得越慢;Q值越小,能量衰减越严重;（2）尺度越小,信号中保留的能量越少;（3）对于脉冲源来说在理想的无衰减介质(即Q趋近于∞)中传播时,信号在不同尺度内的能量相同. 利用尺度能量公式,可从反射地震资料中直接估计品质因子Q（即衰减因子）,也可以提取不同尺度的能量衰减剖面作为储层描述的属性参数,用来进行岩性识别和指示气藏,与经典的谱比法相比,避免了谱比法所面临的双时窗问题以及进行谱估计的窗选择问题. 理论模型试验表明了本文方法的正确性和有效性.     

浙江地区地震波衰减和场地响应研究

根据浙江数字地震台网记录的发生在浙江省及其附近海域的中小地震的波形资料，采用互相衔接的三段几何衰减模型，利用遗传算法联合计算，通过在频率域内的分析，测定和研究了浙江地区剪切波的几何衰减和非弹性衰减以及场地响应的特征。     

Faulting, shallow seismicity and seismic hazard analysis for the Costa Rican Central Valley

Faulting, shallow seismicity (0–30 km), and seismic hazard of the Costa Rican Central Valley were analyzed. Faults in the study area are oriented northwest or northeast. There is an active fault system in the south flank of the Central Volcanic Ridge and another in the north flank of the Talamanca Ridge. Faults of these systems have generated 15 destructive earthquakes in the area during the last 228 years all of them shallow and their locations show one cluster near the Poas Volcano and another southward the Central Valley. These earthquakes have damaged cities of the Central Valley, two of them destroyed Cartago city and almost 1000 people were killed. Regarding recent seismicity, there are three main seismic sources at the Central Volcanic Ridge: Irazu, Bajo de la Hondura and Poas and other three in the Talamanca Ridge: Puriscal, Los Santos and Pejibaye.A seismic hazard map for the Metropolitan Area of San José has been elaborated, based on local tectonic and seismic information. The area for the hazard computation covers an area of 20×15 km² and includes the zone where the most population and socioeconomic activities are concentrated. The computation analysis are based on areas zones and faults, each one characterized by recurrence parameters, geometry, minimum and maximum magnitude and source depth. A recent local spectral attenuation model, which includes relations for shallow crustal sources and subduction zone earthquakes, has been applied in this study. The seismic hazard results are presented in terms of contour plots of estimated peak ground acceleration (PGA) for bedrock conditions for return period of 50, 100 and 500 years. In the Central Park of San Jose City the following PGA values were found: 0.29g for 50 years, 0.36g for 100 years, and 0.53g for 500 years.     

井-地观测系统（3D+VSP）和基于模型的处理技术有关问题的讨论

对了解地下构造和描绘地下油气分布来说,用地面地震勘探技术可以分辩15～30 m的地层.低分辨率和缺乏详细的速度信息会降低地面地震勘探的效果,在研究更薄细的油气储层和研究已研究过的老区时,我们要采用一些更先进的勘探技术,对储层进行更精确的预测和描绘.本文分析了 VSP、测井和地面地震勘探各自的优势、不足和局限,采用VSP方法,我们只能得到井周围一定空问内的丰富速度信息,而详细的速度模型可以极大地提高地面地震勘探的效果.据此分析,组合VSP和地面地震方法优势的井-地观测系统(3D VSP),可有效地增加信息量,并可通过引入VSP数据处理时所用的基于模型的处理技术,来使井-地观测系统的优势进一步扩大,在此基础上,如果与精细测井相结合,可以更好的寻找和开发剩余油气,这将是我们今后重要的发展方向.