高速屏蔽层下广角地震波场分析及成像方法

在高速玄武岩屏蔽区进行勘探存在一定困难,本文以我国南 方典型靶区为例,提出了解决问题的方法. 首先从理论上探讨了广角地震信号存在的条件以 及在应用过程中所面临的问题. 利用正演理论模型分析了广角地震信号的运动学特征以及噪 声对广角反射的影响. 利用物理模型对广角地震的动力学特征、波场特征及动校拉伸现象进 行了分析. 在上述基础上,对广角地震数据作了精细处理,完成了玄武岩下的地震成像.     

Sub-basalt decoupled walkaway VSP imaging

Walkaway vertical seismic profile (VSP) data acquired in basalt‐covered areas can be used to improve knowledge of the sub‐basalt structure. A synthetic example and a case study from the North Atlantic (UK) show that elastic two‐way downward‐continuation migration combined with the stationary‐phase principle is well suited to the processing of VSP data. Vector data are processed using decoupled elastic migration algorithms in both isotropic and anisotropic media. To illustrate the value of decoupled imaging equations, conventional PP imaging is carried out on the enhanced VSP data and compared with the decoupled scheme. Decoupled vector migration operates directly on the displacement vector, and uses various wave modes. Downgoing waves are migrated to image basalt lava flows and measure their anisotropy. Upgoing waves are used for high‐resolution sub‐basalt imaging.     

Seeing below the basalt – offshore Faroes

Exploration in the basalt covered areas of the Faroes offshore has always suffered from poor seismic imaging below the basalt. Long offset 2D and 3D seismic data were acquired and a significant improvement in the seismic image below top basalt has been achieved. Deep towing of the source and receiver cables helped by extending the seismic bandwidth towards lower frequencies. Bubble‐tuned rather than conventional peak‐tuned source arrays gave little, if any, incremental benefit. The improvement in the imaging comes primarily from the approach to processing the data. High frequencies (dominantly noise) are filtered out of the data early in the processing to concentrate on the low frequency data. Careful multiple removal is important with several passes of demultiple being applied to the data using both Surface‐Related Multiple Elimination (SRME) and Radon techniques. Velocity analysis is performed as an iterative process taking into account the geological model. Reprocessing legacy 2D surveys, acquired with wide‐ranging parameters, using these processing techniques improved these datasets significantly, indicating that sub‐basalt imaging seems to be more sensitive to processing than to the choice of acquisition parameters.     

伽师强震群区上地壳三维速度层析成像——人工爆破和天然地震的联合反演

伽师震区位于天山褶皱、帕米尔构造弧与塔里木块体三个构造单元的交接地带,近年来该区发生了一系列的强震活动.为进一步获得该震区详细的地壳速度结构,本文利用人工爆破和天然地震资料联合反演的方法,对1997年新疆伽师震区布设的三维人工地震透射台阵和流动地震台网的资料进行处理,重建了台阵下方上地壳三维速度扰动图像,并结合地震活动分布,对伽师强震群的地震成因作出进一步分析.结果表明研究区上地壳速度结构在纵向和横向上具有明显的非均匀性,随着深度的逐渐加深,震区下方以萨如锡为中心的低速异常体逐步被高速异常体所替代.自12 km深度开始,在与强震群震中相应的位置上,明显出现沿北北西向的高P波速度异常体,在其周围为相对低速分布,呈现出低速条带环绕高速条带的分布格局,V_P/V_S在相同的位置上也表现为高值分布.这种结构上的差异可能与伽师强震群发生有密切关系.16 km深度的P波速度层析图表明,伽师强震群发生在地壳相对高速扰动区内或是高速扰动向低速扰动过渡的边缘,壳内高速体的存在为强震的孕育和发生提供了重要基础.     

应用低频信号提高高速玄武岩下的成像质量

针对高速玄武岩屏蔽层下深层成像困难的实际问题,采用波动方程波场数值模拟技术,根据玄武岩地层的特点,设计三个相应的简单高速玄武岩模型,通过对深层反射地震信号能量的分析,说明了低频地震信号既具有较强的穿透薄高速玄武岩屏蔽层的能力,也具有减弱因粗糙表面所产生的绕射噪音的能力。一个完整的2D玄武岩模型的模拟试验证明了利用低频信号可以提高高速玄武岩屏蔽层下深层成像的质量,实际资料的低通滤波处理也取得了预期的效果。     

Sub-basalt imaging problems and the application of Artificial Neural Networks

In many areas of the world, the presence of shallow high velocity, highly heterogeneous layers complicate seismic imaging of deeper reflectors. Of particular economic interest are areas where potentially hydrocarbon-bearing strata are obscured by layers of basalt. Basalt layers are highly reflective and heterogeneous. Using reflection seismic, top basalt is typified by a high-amplitude, coherent reflector with poor resolution of reflectors below the basalt, and even bottom basalt. Here, we present a new approach to the imaging problem using the pattern recognition abilities of a back-propagation Artificial Neural Network (ANN). ANNs are computational systems that attempt to mimic natural biological neural networks. They have the ability to recognize patterns and develop their own generalizations about a given data set. Back-propagation neural networks are trained on data sets for which the solution is known and tested on the data that are not previously presented to the ANN in order to validate the network result. We show that Artificial Neural Networks, due to their pattern recognition capabilities, can invert the medium statistics based on the seismic character. We produce statistically defined models involving a basalt analogous layer, and calculate full wavefield finite difference synthetic seismograms. We vary basalt layer thickness and source frequency to generate a synthetic model that produces seismic that is similar to real sub-basalt seismic, i.e. high amplitude top basalt reflector and the absence of base basalt and sub-basalt events. Using synthetic shot gathers, generated in a synthetic representation of the sub-basalt case, we can invert the velocity medium standard deviation by using an ANN. By inverting the velocity medium standard deviation, we successfully identified the transition from basalt to sub-basalt on the synthetic shot gathers. We also show that ANNs are capable of identifying the basalt to sub-basalt transition in the presence of incoherent noise. This is important for any future applications of this technique to the real-world seismic data, as this data is never completely noise-free. There is always a certain level of residual (noise remaining after initial noise filtering) environmental/ambient noise present on the recorded seismics, hence, neural network training with noise-free synthetic seismic is less than optimal.     

基于无人机影像的九寨沟地震建筑物震害定量评估

利用2017年8月8日九寨沟7.0级地震震后获取的无人机影像,结合地面震害调查资料,分析各类建筑物震害特征,建立建筑物震害无人机遥感解译标志;选取地震灾区漳扎镇(部分区域)和荷叶寨2个区域作为研究区,进行了无人机遥感建筑物震害提取,基于遥感震害指数进行了震害定量评估,并与现场建筑物震害调查统计结果进行了比较验证。结果显示,遥感解译建筑物震害与实际震害程度相吻合,表明利用震后快速获取的高分辨率无人机影像,可以较为准确地识别建筑物震害,进而为地震灾害定量评估和应急救援辅助决策提供重要参考。     

玄武岩内非均质性对地震资料的影响:印度德干圈闭的一个研究实例(英文)

油气勘探需要重新开发玄武岩下目的层的成像技术。我们今天遇到的最重要的问题之一是玄武岩下地震成像。近年来,由于利用长炮检距这个问题部分已经得到解决。然而,由于地表以及玄武岩的内部非均质性引起的复杂波形,近炮检距的资料仍然不能充分地被利用。影响玄武岩下成像的近垂直入射资料对于了解玄武岩层内部结构是十分有用的。对比较均匀的玄武岩目标层可选用转换波。这里我们利用几个合成模型例子重点描述了更接近现实的非均质玄武岩流引起的实际困难。模拟计算了全波地震记录以有助于了解玄武岩内沉积物对地震资料的影响。本文介绍了印度德干圈闭的一个研究实例。首先探讨了夹层沉积物对整个地震成像的影响。其次利用该区声测井资料以反射系数法计算全波场响应与实际地震资料对比, 通过一系列速度-深度剖面探讨了利用模式转换波(顶底玄武岩界面上的P波转化为S波或相反)对印度库奇地区玄武岩下成像的可行性。通过与野外资料相比较证明玄武岩中多个薄层的效应严重恶化我们所要解释和利用的图象的质量。     

First seismic imaging results of tectonically complex structures at shallow depths beneath the northwest Canterbury Plains,New Zealand

Seismogenic structures underlie many regions of the vast Canterbury Plains on the South Island of New Zealand. Most of these structures are hidden beneath a layer of rapidly deposited Late Pleistocene sediments, the youth and thickness of which make the general application of conventional paleoseismological studies impractical. In an attempt to improve our understanding of potentially active structures in this region, we have acquired, processed and interpreted shallow seismic reflection data across the northwest Canterbury Plains. To separate the useful reflected signals from unusually high amplitude ambient and source-generated noise, we subjected the data to a specially tailored processing scheme that included time- and space-variant spectral balancing, custom static corrections and mutes, F-X deconvolution, DMO corrections and finite-difference migration. The final stacked and migrated seismic sections supply high-resolution images of the basement and overlying layered Cretaceous- to Quaternary-age supracrustal rocks that have been complexly faulted and folded. At one location, the uppermost Late Pleistocene layers appear to have been gently buckled.     

Quantitative evaluation of crosshole seismic reflection images using physical model data1

Physical models give us a known geometry with which to compare our processed reflection images and therefore our imaging techniques. We show how this comparison may be quantified in order to evaluate processed images properly. A crosshole data set was acquired through a model interrogated at ultrasonic frequencies using Durham University's physical modelling system. Various reflectivity images were obtained using processing sequences which include deconvolution, wavefield separation and migration. An error-energy scheme was used to assess the quality of these images, by comparing them against a best-fit depth model. A synthetic data set was also used to evaluate the imaging capability of the crosshole geometry and the effectiveness of the different processing schemes.     

玄武岩内非均质性对地震资料的影响:印度德干圈闭的一个研究实例

Hydrocarbon exploration interests have renewed the need for developing new sub basalt imaging techniques. One of the most important problems encountered today is seismic imaging below basalt. In recent years, this problem appears to have been overcome partly by using long offset seismic data. However near offset data are yet to be fully utilised due to the complex waveform caused by the surface as well as internal heterogeneity of the basalts. The near normal incidence data, which influence the sub-basalt imaging, are highly useful to understand the internal structure within a basalt layer. The use of converted waves for such targets has been proposed as an alternative in a rather homogeneous basalt layer. With a few synthetic modelling exercises here we highlight the practical difficulties in dealing with more realistic and heterogeneous basalt flow. Full waveform seismograms are computed to understand the effects of intra-trappean sediments on the seismic data. A case study from the Deccan Traps of India is presented in this paper. First, we discuss the effects of intercalated sediments on the overall seismic image. Later, the sonic log data from the field are used to compute the full wave-field response using the reflectivity method and compared with the field data. The feasibility of using mode converted waves （P to S and vice-versa at the top and bottom basalt interfaces） for sub-basalt imaging in Kutch region is discussed through a series of velocity-depth profiles. By comparing with the field data we demonstrate that the effects of multiple thin layering within the basalt can strongly deteriorate the image we seek to interpret and exploit.     

P‐wave velocity distribution in basalt flows of the Enni Formation in the Faroe Islands from refraction seismic analysis

The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. V_P, V_S and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data.     

Seismic attenuation in Faroe Islands basalts

We analysed vertical seismic profiling (VSP) data from two boreholes at Glyvursnes and Vestmanna on the island of Streymoy, Faroe Islands, to determine the magnitude and causes of seismic attenuation in sequences of basalt flows. The work is part of SeiFaBa, a major project integrating data from vertical and offset VSP, surface seismic surveys, core samples and wireline log data from the two boreholes. Values of effective seismic quality factor (Q) obtained at Glyvursnes and Vestmanna are sufficiently low to significantly degrade the quality of a surface reflection seismic image. This observation is consistent with results from other VSP experiments in the North Atlantic region. We demonstrate that the most likely cause of the low values of effective Q at Glyvursnes and Vestmanna is a combination of 1D scattering and intrinsic attenuation due to seismic wave‐induced fluid flow within pores and micro‐cracks. Tests involving 3D elastic wave numerical modelling with a hypothetical basalt model based on field observations, indicate that little scattering attenuation is caused by lateral variations in basalt structure.     

Enhancing hardrock seismic images: Reprocessing of high resolution seismic reflection data from Vihanti,Finland

Seismic reflection data were acquired across volcanic hosted massive sulfides (VHMS) of Vihanti in order to improve the understanding of the regional geological setting. Commercially processed seismic data from Vihanti are of good quality, but reprocessing can be used to extract additional information about geological structures. Especially, careful velocity analysis influences the quality of seismic images. Differentiating reflections caused by fractures from those caused by lithological contacts is very important for exploration and geological modeling. Reflections from fracture zones known from drilling stack with lower velocity (~ 5100 m/s) compared to typical stacking velocities of the Vihanti area (> 5500 m/s). The reprocessing also indicated that fracture zones are better imaged with low frequencies due to the better overall continuity of the fault zones at scales of hundreds of meters rather than at shorter seismic wavelengths.In full stacks, long offset data can mask structure close to the surface. More detailed seismic images of the shallow subsurface emerged by preferentially stacking short offset data wherever acquisition and processing lines lay close together and were nearly straight. Long offset data remains valuable for imaging deeper structures as well as dipping reflectors. Cross-dip-analysis revealed a bright diffractor located near the base of the Vihanti volcanic basin at 1.5 km depth. The seismic data allow a geological interpretation in which the Vihanti structure has developed through significant thrust faulting and displacement of the lithological contacts. Gentle folds that were formed prior to faulting are visible as undulating reflectivity in seismic sections. The reprocessed seismic section indicates a potential deep extension of the ore-hosting altered volcanic and calc-silicate rocks previously unexplored.     

雷琼拗陷地壳结构特征

根据1984～1985年在雷琼地区完成、长约210 km的白延—江洪深地震测深资料,应用有限差分法、时间项法反演和二维非均匀动力学射线追踪等方法,得到白延—江洪地壳剖面的二维非均匀构造图像.结果显示：琼北及雷州半岛地壳厚约25～26 km,地壳平均速度在琼州海峡南段琼北凸起约616 km·s-1,海峡北段雷琼拗陷区约570 km·s-1;地壳速度偏低,下地壳尤为明显,仅约63～65 km·s-1;结晶基底显示为厚约2～4 km、速度55～60 km·s-1的速度梯度带,Moho界面上、下界面速度分别是65 km·s-1和805 km·s-1的一级速度间断面,没有发现明显的壳幔过渡带.白延—江洪剖面位于雷琼拗陷更新世玄武岩覆盖区,推测地壳厚度减薄以及低速度构造主要反映了该区域自新生代以来一直延续至第四系晚期的多期火山岩浆喷发热地幔对地壳底部的底浸、拆层与地壳仍处于相对温热状态相关.     

2013年四川芦山7.0级地震烈度遥感评估

2013年4月20日四川芦山MS7.0级地震发生后,在灾区应急获取了多种高分辨率航空和无人机遥感影像,并快速解译提取了灾区建筑物震害信息.采用地震烈度遥感定量评估方法,利用2008年汶川8.0级地震等震后震害遥感解译和现场调查研究确定的经验震害遥感定量评估模型,获得了芦山地震灾区126个主要居民点的地震烈度遥感评估结果,并据此圈画了地震烈度分布遥感评估图.结果显示,本次地震Ⅸ度区面积约150km2,Ⅷ度区面积约900km2.该结果在第一时间(4月21日晚)提供给了中国地震局地震现场应急指挥部.对比分析显示,地震烈度遥感快速评估结果与中国地震局4月25日公布的地震烈度图,以及与笔者在现场实地进行的建筑物震害详细调查结果基础上评定的地震烈度具有较高的一致性.表明强烈地震发生后,借助于快速获取的灾区高分辨率遥感影像,可以快速估计地震烈度分布,对地震灾区灾情估计和抗震救灾工作具有十分重要的参考意义.     

The application of combined gravity and seismic data formation separation for revealing deep structure

In some oilfields with 3D seismic data, the deeper structure cannot be observed due to poor quality deep seismic data. Layer stripping using both seismic and gravity data is a solution for this problem but it cannot get satisfactory results because the horizontal variations in formation density are ignored. We present a variable-density formation separation technique to address this problem. Based on 3D seismic depth data and laterallyvariable density derived from 3D seismic velocity data, the upper formation gravity effect is calculated by forward modeling and removed from the Bouguer gravity. The formation-separated gravity anomaly with variable density is obtained, which mainly reflects the deeper geological structure. In block XX of North Africa, the shallow formations seismic data is excellent but the data at the top of basement is poor. The formation-separated gravity anomaly processed under the control of 3D seismic data fits well with the known seismic interpretation and wells. It makes the geological interpretation more reliable.     

A real-time AI-assisted seismic monitoring system based on new nodal stations with 4G telemetry and its application in the Yangbi MS 6.4 aftershock monitoring in southwest China

A rapidly deployable dense seismic monitoring system which is capable of transmitting acquired data in real time and analyzing data automatically is crucial in seismic hazard mitigation after a major earthquake. However, it is rather difficult for current seismic nodal stations to transmit data in real time for an extended period of time, and it usually takes a great amount of time to process the acquired data manually. To monitor earthquakes in real time flexibly, we develop a mobile integrated seismic monitoring system consisting of newly developed nodal units with 4G telemetry and a real-time AI-assisted automatic data processing workflow. The integrated system is convenient for deployment and has been successfully applied in monitoring the aftershocks of the Yangbi M_S 6.4 earthquake occurred on May 21, 2021 in Yangbi County, Dali City, Yunnan in southwest China. The acquired seismic data are transmitted almost in real time through the 4G cellular network, and then processed automatically for event detection, positioning, magnitude calculation and source mechanism inversion. From tens of seconds to a couple of minutes at most, the final seismic attributes can be presented remotely to the end users through the integrated system. From May 27 to June 17, the real-time system has detected and located 7 905 aftershocks in the Yangbi area before the internal batteries exhausted, far more than the catalog provided by China Earthquake Networks Center using the regional permanent stations. The initial application of this integrated real-time monitoring system is promising, and we anticipate the advent of a new era for Real-time Intelligent Array Seismology (RIAS), for better monitoring and understanding the subsurface dynamic processes caused by Earth's internal forces as well as anthropogenic activities.     

Deep seismic reflection profiling of sedimentary basins offshore Brazil: Geological objectives and preliminary results in the Sergipe Basin

The first deep seismic reflection profiles offshore Brazil were acquired in Campos Basin and processed to 10 s TWT in 1984. Starting in 1989, Petrobrás acquired an extensive data set of deep seismic profiles using special acquisition equipment capable of effectively penetrating through the sedimentary layers and imaging the whole crustal architecture. These deep (18 s TWT) seismic reflection profiles extend across the Atlantic-type marginal basins, from the platform to the deepwater province, presently considered frontier regions for petroleum exploration. This work addresses the geological objectives of a deep seismic profile in the Sergipe Basin and discusses the results obtained by integrating regional seismic, gravity and magnetic data. When combined, these data provide evidence that deep seismic reflectors observed in the Sergipe Basin are related to intracrustal-upper mantle structures rather than sedimentary features. The deep seismic reflection profile in the Sergipe Basin also suggests that, rather than a non-volcanic passive margin, the deepwater extension of this basin is marked by several magmatic structures, including thick wedges of seaward-dipping reflectors and volcanic plugs. These magmatic features are associated with basinforming processes resulting from lithospheric extension during the breakup of Gondwana in the Early Cretaceous and subsequent emplacement of oceanic crust. These results are compared to the crustal scale structures observed in the Campos Basin, in the southeastern margin of Brazil. The interpretation of the deep structure of these basins indicates that final separation between the South American and African plates formed passive margins characterized by different patterns of crustal attenuation underlying the rift blocks.     

重力-地震联合变密度剥层技术研究深层构造

在某些已经完成三维地震勘探工作的油气田区,因深层地震资料不好,深层构造难以得到解决。虽然结合地震资料与重力资料的重力剥皮技术得到开发和使用,但由于它没有考虑地层密度的横向变化,因此未能获得满意的处理效果。为了改善该类地区的重力勘探效果,本文提出了变密度剥层方法,它利用各地层的三维地震深度资料和三维地震层速度转换的具横向变化的密度资料,通过正演其重力效应,并从布格重力异常中去除,从而得到主要反映深部地质构造信息的变密度剥层重力异常。在北非某区块,除基底顶面外,各层地震资料良好,通过重力-地震联合变密度剥层处理得到的剥层重力异常,与已知地震和钻井资料吻合良好,为深部基底研究奠定了良好的基础。