塔里木盆地深反射地震资料处理

反射地震资料的属性与油气勘探的最终目标——油、气有密切联系。利用地震波场的属性中隐藏着识别油气藏的密码进行油气检测,是一项新的油气检测技术。有油气存在时其地震属性会在局部产生非常规变异,寻找这种变异的差异,就可以进行油气检测,这也就为油气圈闭的地球物理评价提供了基础。通过在塔里木油田的试验表明,上述检测油气的属性差异分析方法技术是有效的。建议在油气钻探之前,首先进行油气检测和圈闭评价,以减少风险,提高油气勘探的成功率。     

TSP原始数据文件的读取及与国产地震仪数据之间的转换

近年来,TSP(Tunnel Seismic Prediction)超前地质预报系统在我国得到了比较广泛的应用,但厂家未提供其数据文件存储格式。为此,采用UltraEdit软件对多个TSP数据文件进行分析,弄清了其原始数据的总体存储结构。同时还简单介绍了国际标准地震数据SEG-2文件的存储格式,并在Matlab语言环境下编制了读取SEG-2数据文件与TSP数据文件的程序,实现了将SEG-2文件转换成TSP数据文件的功能。因此,在今后的工作中,可直接应用TSP的数据处理软件包TSPwin来处理国产的超前地震数据。     

一种地震共反射点道集数据的叠前反演方法

叠前地震资料含有地层的纵波、横波速度和密度等信息。利用叠前反演获得隐藏在地震数据中的这些基本参数后,即可揭示大量岩性及孔隙流体性质的信息。这里推导了平面波在层状弹性介质中传播的正演算子,提出了一种基于Brent正交搜索方向组的叠前三参数反演方法．该方法不需要求解庞大而复杂的导数矩阵。通过自适应退火因子和罚函数来处理约束条件,提高了算法的稳定性。将K—L变换引入到方向置换过程,有效防止了搜索方向组的线性相关。经理论模型和油田实际数据的反演结果表明,该反演方法是一种利用叠前地震数据进行储层预测的有效手段。     

Identification of karst features with spectral analysis on the seismic reflection data

The use of high-resolution seismic tools to investigate karst terrain has had limited success due to the ambiguous seismic events caused by lateral heterogeneities of the sedimentary cover, as well as the superimposition of surface waves. To address these problems, a spectral analysis scheme is proposed on the seismic reflection data to characterize karst features in this paper. First, the S-transform is applied to suppress the ground roll. As a type of time–frequency transform, the S-transform provides frequency-dependent resolution while maintaining a direct relationship with the Fourier spectrum. Application of a filter based on the S-transform attenuates ground roll in a time–frequency domain. Then, the continuous wavelet transform is employed for the time–frequency mapping of potential karst features, which are featured as frequency differentiation anomalies. Synthetic model and real-data example validate the efficiency of the spectral analysis scheme for the detection of karst features.     

Heterogeneous structure of western Nankai seismogenic zone deduced by multichannel reflection data and wide-angle seismic data

The 1946 Nankai earthquake (Ms=8.2) at the forearc region of the western Nankai Trough showed slow slip deformation off Cape Muroto, which did not propagate until the western end of the Nankai seismogenic zone. New seismic investigations show a low-velocity layer (LVL) on the subducting oceanic crust in the coseismic area. Two prestack depth-migrated sections show reflectivity events in the clay-rich boundary layer on the oceanic crust. Narrowly spaced imbricated slices develop in the nonrupture area. The reflective boundary layer indicates probably that underplating develops in the nonrupture area rather than the coseismic area. It is suggested that the friction is larger in the nonrupture area than the coseismic area because of the lack of LVL on the oceanic crust, the well developed underplating and the narrowly spaced imbricated thrusts in the nonrupture area. The topographic high of the oceanic crust with about 50 km width and maximum 3 km height is also revealed and is related to bending and thickening of the oceanic crust, the well developed underplating and the narn spaced imbricated thrusts in the nonrupture area. These structural characters may be the reason why the slow slip deformation did not propagate until the western end of the Nankai seismogenic zone and toward the trough side.     

A method for avoiding artifacts in the migration of deep seismic reflection data

Conventional wave-equation-based migration of deep seismic reflection data can produce severe artifacts, which appear as broad circular arcs or “smiles”, due to the existence of apparent truncations of reflections on the stack section arising from poor signal penetration, changes in orientation of the acquisition profile, and the existence of strong overlying lateral velocity variations. These artifacts limit the interpretation of deep seismic profiles, because they obscure weak reflections and reflection truncations that may, e.g., indicate the presence of subsurface faults. Here I present a new migration algorithm, in which each sample of the stack is migrated to a short linear segment whose position and dip are determined by its original position on the stack, an estimate of the local apparent dip at that point, and a user-specified migration velocity. No subjective interpretation of reflections on the stack section is required, and the algorithm produces no arc-like migration artifacts. The degree of lateral smearing can be easily controlled, allowing reflection truncations to be revealed. In practice, the algorithm is most effectively applied to data that have been coherency-filtered to remove low amplitude noise, which would otherwise be preserved.     

Active tectonics of the Yizre'el valley, Israel, using high-resolution seismic reflection data

We present a series of high-resolution seismic reflection lines across the Yizre'el valley, which is the largest active depression in Israel, off the main trend of the Dead Sea rift. The new seismic reflection data is of excellent quality and shows that the valley is dissected into numerous small blocks, separated by active faults. The Yizre'el valley is found to consist of a series of half grabens, rather than a single half graben, or a symmetrical graben. The faults are generally vertical and appear to have a dominant strike-slip component, but some dip-slip is also evident. A marked zone of compression near Megido is associated with the intersection of the two largest faults in the valley, the Carmel fault and the Gideon fault. Variable trend of the faults reflects the complexity of the local geology along the boundary between the wide NW–SE trending Farah–Carmel fault zone and the E–W trending basins and ranges in the Lower Galilee. This tectonic complexity is likely to result from a highly variable stress pattern, modified by the structures inside it. Normal faulting in the valley occurred at an early stage of its development as a tectonic depression. However, strike-slip motion on the Carmel fault, and possibly also on some of the other faults, appears to have started together with the onset of normal faulting. Earthquake hazard in the area appears to be uniform as faults are distributed throughout the Yizre'el valley.     

利用XML格式解析原理批量进行成果资料格式检查

由于科学技术的进步,各种项目成果资料获取的周期缩短,数量成逐年上升趋势.所有地震勘探项目的开始和结束都是地质资料的开始获取和完成获取的过程,因而地震勘探项目的成果资料编辑、检查、存档管理工作一直处于最重要的位置.自从归档资料统一格式之后,各种归档的维护工作就相对变得容易,但在归档之前的文档检查校正工作却显得繁琐和不易.这里试图从常见的Word电子文档入手,利用其公开的Word XML格式解析原理,通过编程实现了常见排版和错误的批量检查工作,并能够极大提高在资料归档之前的检查工作效率和减少人工翻阅检查的难度.经实际工作验证,效果很好.     

New insights into the crustal structure of the North German Basin from reprocessing of seismic reflection data using the Common Reflection Surface stack

The influence of deep crustal processes on basin formation and evolution and its relation to current morphology is not well understood yet. A key feature to unravel these issues is a detailed seismic image of the crust. A part of the data recorded by the hydrocarbon industry in the late 1970s and 1980s in the North German Basin were released to the public recently. The seismic reflection data were recorded down to 15 s two-way travel time. The mean Common Midpoint fold of about 20 is relatively low compared to contemporary seismic acquisitions. The processing of the 1980s focussed on the sedimentary structures to explore the hydrocarbon potential of this area. We applied the Common Reflection Surface stack technique to the data sets, which is well suited for low-fold data. The reprocessing was focussed on the imaging of the subsedimentary crustal range. The reprocessed images show enhanced reflections, especially in the mid and lower crustal part. Also, the image of the salt structures in the graben area was improved. Furthermore, the reprocessed images indicate an almost flat Moho topography in the area of the Glückstadt Graben and an additional lower crustal structure, which can be correlated with a high-density body found in recent gravity modeling studies.     

Cordilleran front range structural features in northwest Montana interpreted from vintage seismic reflection data

Industry seismic reflection data spanning the Rocky Mountain Cordillera front ranges of northwest Montana were reprocessed and interpreted in this study. Five seismic profiles represent 160 km of deep reflection data collected in 1983 that span the eastern Purcell anticlinorium, Rocky Mountain Trench (RMT), Rocky Mountain Basal Décollement (RMBD), and Lewis thrust. The data were reprocessed using modern techniques including refraction statics, pre-stack time migration (PSTM), and pre- and post-stack depth migration. Results indicate the RMBD is 8–13 km below the Earth's surface and dip 3–10° west. Evidence for the autochthonous Mesoproterozoic Belt and basal Cambrian rocks beneath the RMBD is present in all of the profiles and appears to extend east of the RMT. The Lewis thrust was identified in the seismic profiles and appears to sole into the RMBD east of the RMT. The RMT fault system has a dip displacement of 3–4 km and forms a half graben filled with 1 km of unconsolidated Tertiary sedimentary deposits. The RMT and adjacent Flathead fault systems are interpreted to be structurally linked and may represent a synthetic, en echelon fault system.     

松辽盆地北部石炭—二叠纪地层厚度：来自深反射地震的证据

在松Ⅰ～松Ⅵ6条深反射地震剖面上依据震相特征并结合地质演化过程分析识别石炭—二叠系,分析其层位反射特征及同相轴特征,建立松辽盆地北部地区石炭—二叠系地震震相特征识别标志。主要震相特征为中—强反射能量,局部存在高值,整体同相轴连续—较连续,由南到北、由西向东规律变化。利用已有的钻遇石炭—二叠系的探井资料和地震波传播速度,得到研究区时间-深度转换关系。由研究区石炭—二叠系相位追踪对比后的反射时间分布,计算石炭—二叠纪地层顶界面与底界面在深反射地震时间剖面上的走时之差,通过时深转换得出研究区基底石炭—二叠纪地层厚度和分布。全区地层分布不很均匀,主要在中央隆起带以东地区及西部断陷区出现2个厚度高值区;地层整体由浅到深大致可分为上、中、下3层,且地震震相特征互不相同。松辽盆地北部石炭—二叠系分布对东北地区主力油层之下的深层油气勘探提供了有价值的依据。     

基于A- CGAN的深反射地震数据随机噪声压制方法研究

基于深度学习的地震数据噪声压制方法是当前地震数据去噪处理的重要方向。深度学习方法突破了传统滤波处理的局限,在对常规地震数据的噪声压制中表现出效率高、信噪分离效果好的特点。但针对深部弱有效反射数据,当前的深度学习方法特征提取能力有限,难以取得较好的去噪效果。笔者等结合深反射地震数据特点,针对当前深度学习噪声压制方法在特征提取及对数据集依赖上的局限,提出了基于注意力循环生成对抗网络(Attention Cycle- Consistent Generative Adversarial Networks,A- CGAN)的深反射地震数据随机噪声压制方法。借助循环一致生成对抗网络(Cycle- Consistent Generative Adversarial Networks,Cycle- GAN)的域映射思想,降低对数据集的要求。为了构建适用于深反射地震数据的去噪网络,从3个方面对Cycle- GAN进行改进：在Cycle- GAN的生成器（去噪器）中加入残差结构和注意力机制,用于加深网络深度和提高其特征提取能力;在Cycle- GAN的鉴别器中使用块判决,提高鉴别精度和准确度;在损失函数部分加入感知一致性损失函数,提升网络模型恢复纹理细节信息的能力。通过合成地震数据和实际深反射地震数据测试,验证了优化算法的有效性,体现了良好的应用价值。     

基于A- CGAN的深反射地震数据随机噪声压制方法研究

基于深度学习的地震数据噪声压制方法是当前地震数据去噪处理的重要方向。深度学习方法突破了传统滤波处理的局限,在对常规地震数据的噪声压制中表现出效率高、信噪分离效果好的特点。但针对深部弱有效反射数据,当前的深度学习方法特征提取能力有限,难以取得较好的去噪效果。笔者等结合深反射地震数据特点,针对当前深度学习噪声压制方法在特征提取及对数据集依赖上的局限,提出了基于注意力循环生成对抗网络(Attention Cycle- Consistent Generative Adversarial Networks,A- CGAN)的深反射地震数据随机噪声压制方法。借助循环一致生成对抗网络(Cycle- Consistent Generative Adversarial Networks,Cycle- GAN)的域映射思想,降低对数据集的要求。为了构建适用于深反射地震数据的去噪网络,从3个方面对Cycle- GAN进行改进：在Cycle- GAN的生成器（去噪器）中加入残差结构和注意力机制,用于加深网络深度和提高其特征提取能力;在Cycle- GAN的鉴别器中使用块判决,提高鉴别精度和准确度;在损失函数部分加入感知一致性损失函数,提升网络模型恢复纹理细节信息的能力。通过合成地震数据和实际深反射地震数据测试,验证了优化算法的有效性,体现了良好的应用价值。     

A possible lower crustal flow channel in the Middle Urals based on reflection seismic data

The Europrobe Seismic Reflection Profiling in the Urals Experiments (ESRU) reflection seismic data from the Middle Urals images c. 10‐km thick band of strong, subhorizontal lower crustal reflectivity and a thinning of the crust that is associated with the East Uralian Zone, a broad strike‐slip fault system containing high‐grade metamorphic rocks and syn‐orogenic to post‐orogenic granitoids. The lower crustal reflectivity consists of discontinuous to continuous, high‐amplitude reflections. Reflections are subparallel to slightly oblique and have a layered to oblate appearance. Geometrical relationships indicate that the reflectivity post‐dates fault activity, suggesting that late‐orogenic processes modified the lower crust. The surface geology indicates that the conditions for lower crustal flow were met in the East Uralian Zone. We suggest that the lower crustal reflectivity imaged by the ESRU data is related to a flow channel that developed at the base of the crust in the interior of the orogen.     

Some interpretational aspects of time cross-sections based on seismic reflection data from the Sea of Okhotsk

The results of investigations in the Sea of Okhotsk by the seismic reflection method are considered with special attention paid to the acoustic anomalies. Two types of such anomalies are definable in the basin: local (“gas channel”) and areal (“turbid sequences”). It is shown that the acoustic anomalies distinguishable in the single channel seismic profiling (SCSP) records are characteristic of sediments with partial gas filling. Many local anomalies documented on the Sakhalin slope are spatially connected with zones of faulting. The simultaneous presence of acoustic anomalies and BSR in the upper part of the section serves as a reliable diagnostic feature of gas hydrates. The dispersion of the velocities in the wave packet characterizing the acoustic anomaly associated with the lower part of the sedimentary cover in the Kurile Basin points to the significant thickness of the gas-bearing sediments in the latter.     

Low velocity layer characterization in the Niger Delta: Implications for seismic reflection data quality

The problems of poor data quality and statics in seismic surveys have been attributed to lack of proper understanding of the low velocity layer characteristics of the area from which such data was acquired. Downhole seismic refraction survey was therefore conducted at twenty (20) borehole locations within parts of Niger Delta, Nigeria to determine the low velocity layer characteristics of the area using the Geometrics Stratavisor NZ11 instrument. The data was processed using Udysys software with spatial representation of the results presented. Static corrections were carried out on reflection seismic data acquired from the study area using Geoscribe II software to determine the effects of the low velocity layer on reflection seismic data in the area. Results of the study revealed that the velocity of the low velocity layer ranged between 144 and 996m/s with a regional average of 407m/s. The thickness of the low velocity layer varied between 3.0 and 9.6m with a mean value of 5.0m. Similarly, the velocity of the consolidated layer ranged between 1449 and 1812m/s with a mean value of 1738m/s. Results of the static correction carried out on the seismic reflection data revealed a substantial improvement in the resolution of the data after static correction. Based on these findings, it is therefore, recommended that shots for reflection seismic survey should be located at a minimum depth of 9.6m in the area to eliminate the effects associated with the low velocity layer.     

地震反射数据Q值估计及其在油气勘探中的应用

油气地震勘探中,地震波在实际介质中传播时会产生明显的吸收衰减现象,导致地震信号中的主频向低频移动,频带宽度变窄,相位发生畸变,制约了地震勘探识别薄层的分辨能力。为了获得更高分辨率的地震成像,文章介绍了一种在地震反射数据成像域进行Q值估计建模的方法与流程,并利用估计出的Q值通过偏移成像对数据进行衰减补偿,实现高分辨率成像。该方法在时间域引入等效Q值的概念,首先在初步黏弹性时间偏移成像域的时窗内,通过数据在时窗内的补偿效果来确定时间域的等效Q值参数,接着通过在深度域层速度上计算成像射线获得时深转换关系,进而对转换到时间域层Q值的等效参数进行时深转换,完成最终深度域Q值建模,最后将该关键参数作为黏弹性叠前深度偏移的输入,进行复杂构造的黏弹性补偿成像。同时使用中国东部某实际地震数据来验证方法的有效性,验证结果表明发展的流程和方法可较好实现深度域负责构造Q值建模和成像域补偿,实现复杂构造高分辨成像。     

Study of subsurface features using seismic reflection data in Erbil area,Kurdistan Region,Iraq

The present study seeks to interpret seismic reflection data for an area measuring approximately about 1,450 km² within the Erbil plain in the Kurdistan Region of Iraq to produce subsurface geological pictures. Thirteen seismic lines were identified by the Iraqi Oil Exploration Company using dynamite sources and then processed; then, a seismic line, BH-15, of around 432.5 km total length was used in this study. Three horizons (reflectors) were selected and identified by tying them directly to Demir Dagh well no. 1, these are H1, a reflector that represents the near top of the Lower Fars Formation; H2, a reflector that represents the top of the Pila Spi Formation; and H3, a reflector that represents the top of the Shiranish Formation. The qualities of the reflectors range between fair and good. All the isochrone and depth maps and 3D view pictures of the reflectors show the same main subsurface structural features trending in a NW–SE direction, namely, Erbil trough, Sherawa trough, and the southeastern plunge of the Demir Dagh anticline. Ninety-five normal and reverse-type faults trending in a NW–SE direction were detected. Also, the maps show that the reflectors tend to increase in dipping and deepening towards the southeastern and eastern parts.     

MapGIS数据向ArcGIS数据格式转换方法的实践与探讨

由于在地图编辑和数据分析方面各具优势,MapGIS和ArcGIS逐渐发展成为目前国内被广泛使用的地理信息系统平台,实现二者间数据安全转换将有利于对两种系统功能的充分利用,大大提高空间数据库的生产效率和使用率。在分析比较两种系统空间建模思想及专有数据格式的基础上,提出了将空间数据从MapGIS平台转换到ArcGIS平台上的基本思路,并通过工作实例进行了数据转换验证,对过程中存在的问题进行了总结,形成了一套较为成熟的技术方法。     

Comparison of P- and S-wave velocity profiles obtained from surface seismic refraction/reflection and downhole data

High-resolution seismic-reflection/refraction data were acquired on the ground surface at six locations to compare with near-surface seismic-velocity downhole measurements. Measurement sites were in Seattle, WA, the San Francisco Bay Area, CA, and the San Fernando Valley, CA. We quantitatively compared the data in terms of the average shear-wave velocity to 30-m depth (V_s30), and by the ratio of the relative site amplification produced by the velocity profiles of each data type over a specified set of quarter-wavelength frequencies. In terms of V_s30, similar values were determined from the two methods. There is <15% difference at four of the six sites. The V_s30 values at the other two sites differ by 21% and 48%. The relative site amplification factors differ generally by less than 10% for both P- and S-wave velocities. We also found that S-wave reflections and first-arrival phase delays are essential for identifying velocity inversions. The results suggest that seismic reflection/refraction data are a fast, non-invasive, and less expensive alternative to downhole data for determining V_s30. In addition, we emphasize that some P- and S-wave reflection travel times can directly indicate the frequencies of potentially damaging earthquake site resonances. A strong correlation between the simple S-wave first-arrival travel time/apparent velocity on the ground surface at 100 m offset from the seismic source and the V_s30 value for that site is an additional unique feature of the reflection/refraction data that could greatly simplify V_s30 determinations.