起伏地表组合震源地震波场定向方法

多个震源组合激发通过改变激发延时,可以得到沿某方向传播的地震波场,即定向地震波场激发技术,它可用于特定目标体的照明与探测之中.本文以水平地表的组合震源定向激发原理为基础,推导了倾斜地表情况下的组合震源定向激发公式,并计算绘制了其理论方向图.另外本文将组合震源波场定向方法推广至任意起伏地表,根据惠更斯菲涅尔原理,提出了旋转坐标方法,即将水平坐标旋转至定向波场传播方向法向的倾斜坐标,可方便地计算震源传播至定向波场波前面的走时,作为组合震源的激发延时.根据本文提出的方法,我们分别计算了倾斜地表条件与复杂地表条件定向地震波场的组合震源延时参数,通过波动方程数值模拟技术得到的波场快照验证了本文方法的有效性.     

基于起伏地表的层析速度反演方法研究(英文)

叠前深度偏移速度分析是地震数据处理方法研究的重点,是影响地震成像效果的关键技术之一。对于地表起伏、地下构造复杂的双复杂地区,常规的叠前偏移速度分析方法是将起伏面校正到固定基准面上或进行表层建模再进行偏移与速度分析。本文提出的基于起伏地表的层析速度反演方法从起伏地表直接进行速度场的更新,可以提高层析速度反演的精度与效率。首先介绍基于起伏地表的角度域共成像点道集的提取方法,以此为基础阐述了起伏地表层析速度反演方法。起伏地表模型和实际资料试算验证了本文方法的有效性。     

起伏地表下地震波传播数值模拟方法研究进展

起伏地表是地震数据的采集、处理和解释中普遍遇到的难题.起伏地表下的地震波传播数值模拟,对起伏地表观测的地震资料处理解释有重要意义.地震波场模拟和地震波走时场分别描述地震波的动力学和运动学信息,为研究地震波传播理论的两种重要途径.本文从地震波场和地震波走时场两方面回顾和总结了起伏地表下的地震波传播数值模拟方法的研究进展,并展示了该领域的一些最新研究成果,为使读者能从中找到突破点,为起伏地表这一勘探领域的经典难题做出贡献.     

起伏地表采集数据的三维直接叠前时间偏移方法

提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性.     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.     

一种起伏地表条件下的照明补偿方法

起伏的地表条件限制了采集孔径范围并造成深层地震照明不足,为改善该类地区的成像质量,本文提出了一种起伏地表条件下的照明补偿方法.首先,基于小波束波场延拓算子和逐步累加的外推方法在波场延拓过程中解决起伏地表面的影响,并引入空间滤波函数压制虚拟层内的偏移噪音;其次,利用局部指数标架对上、下行波场分解,得到局部角度域成像和照明补偿因子.再次,利用计算出的成像值和照明补偿因子,在局部倾角域完成照明补偿.SEG起伏地表模型测试证明了本方法的有效性,深层构造照明度明显加强,不同角度成像振幅更加均衡,该技术为提高起伏地表地区的成像品质提供了新的手段.

     

三维起伏地表条件下地震波走时计算的不等距迎风差分法

三维起伏地表条件下的地震波走时计算技术是研究三维起伏地表地区很多地震数据处理技术的基础性工具.为了获得适应于任意三维起伏地表且计算精度高的走时算法,提出三维不等距迎风差分法.该方法采用不等距网格剖分三维起伏地表模型,通过在迎风差分格式中引入不等距差分格式、Huygens原理及Fermat原理来建立地表附近的局部走时计算公式,并通过在窄带技术中设定新的网格节点类型来获得三维起伏地表条件下算法的整体实现步骤.精度及算例分析表明:三维不等距迎风差分法具有很高的计算精度且能够适应于任意三维起伏地表模型.     

起伏地表复杂介质波动方程有限元数值模拟方法

波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时,求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场,本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面;把单元内的场和波速均看作单元上的线性函数,以适应复杂介质压制边角散射;采用吸收边界条件去除来自截断边界上的反射;采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆,提高了计算效率.对模型的计算表明该方法正确有效.     

起伏地表下的直接叠前时间偏移

提出了一种新的叠前时间偏移方法和流程,可不必应用野外静校正,直接对起伏地表采集的地震数据进行叠前时间偏移.本文采用输入道成像方式,通过基于稳相点原理给出单道数据的走时和振幅计算方法,发展了一个表驱动的叠前时间偏移算法.偏移方法可依据同相轴是否被拉平确定叠加速度和修正近地表速度模型,也可依据拟成像的构造倾角,自适应地确定偏移孔径;后者既减少了偏移计算量,也压制了偏移噪声.文中用二维起伏地表的断陷盆地模型的理论数据验证了所发展方法的成像效果.     

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

可控震源定向照明方法的仿真研究

当野外噪声很强,即使使用组合震源地震也无法获得满意信噪比的地震数据时,本文提出了一种基于可控震源阵列的定向照明控制方法,采用该方法可形成定向地震波.通过仿真研究合成了8激震器可控震源阵列分别采用简单组合及定向照明技术得到的单炮地震记录,可以看出采用合适的延时参数,定向照明单炮地震记录的反射波信噪比高于组合地震情况.定量的计算结果表明,实验条件下采用0.89 ms延时参数,各反射波信噪比分别提高了10.19 dB,3.23 dB和1.02 dB.由此可见,可控震源定向照明地震技术是一种提高原始地震资料信噪比的有效方法.     

时空阵列混场源电磁法理论及模拟研究

针对电磁勘探中的混场源输入条件,提出了一种时空阵列混场源电磁勘探方法.假设输入端同时包含天然电磁场源和可控人工场源,输出端进行多站同步观测;基于多输入-多输出系统的分析方法,推导了混场源条件下的时空阵列方程组.提出了系统响应的求解策略,一次性获得所有同步测道对各个场源独立激励时的响应,并分别求取了大地电磁法及可控源电磁法的各种解释参数.在此理论基础上,设计了针对性的施工方案;可控人工场源的选择形式多样,其类型、数量、激励方式及布设位置等不受限制,所有场源的激励可同时进行,无需分步;测区内各测站的布设方式灵活,可根据需要布设单分量测道或多分量测道;同时在测区外布设张量远参考站.利用数值模拟方法,论证了方法的有效性,并与常规方法进行了对比.结果表明,本方法不仅可提高人工场响应的处理精度,还具有一次野外数据采集,可同时获得天然场及人工场电磁响应结果的突出特点,进而提高采集效率,压制噪声影响;本方法集合了天然场源电磁法和人工场源电磁法的优点,为实现不同类型频率域电磁测深方法数据的统一处理提供了思路.

     

Monitoring of seismic events from a specific source region using a single regional array: A case study

In the monitoring of earthquakes and nuclear explosions using a sparse worldwide network of seismic stations, it is frequently necessary to make reliable location estimates using a single seismic array. It is also desirable to screen out routine industrial explosions automatically in order that analyst resources are not wasted upon detections which can, with a high level of confidence, be associated with such a source. The Kovdor mine on the Kola Peninsula of NW Russia is the site of frequent industrial blasts which are well recorded by the ARCES regional seismic array at a distance of approximately 300 km. We describe here an automatic procedure for identifying signals which are likely to result from blasts at the Kovdor mine and, wherever possible, for obtaining single array locations for such events. Carefully calibrated processing parameters were chosen using measurements from confirmed events at the mine over a one-year period for which the operators supplied Ground Truth information. Phase arrival times are estimated using an autoregressive method and slowness and azimuth are estimated using broadband f{-}k analysis in fixed frequency bands and time-windows fixed relative to the initial P-onset time. We demonstrate the improvement to slowness estimates resulting from the use of fixed frequency bands. Events can be located using a single array if, in addition to the P-phase, at least one secondary phase is found with both an acceptable slowness estimate and valid onset-time estimate. We evaluate the on-line system over a twelve month period; every event known to have occured at the mine is detected by the process and 32 out of 53 confirmed events were located automatically. The remaining events were classified as “very likely” Kovdor events and were subsequently located by an analyst. The false alarm rate is low; only 84 very likely Kovdor events were identified during the whole of 2003 and none of these were subsequently located at a large distance from the mine. The location accuracy achieved automatically by the single-array process is remarkably good, and is comparable to that obtained interactively by an experienced analyst using two-array observations. The greatest problem encountered in the single array location procedure is the difficulty in determining arrival times for secondary phases, given the weak Sn phase and the complexity of the P-coda. The method described here could be applied to a wide range of locations and sources for which the monitoring of seismic activity is desirable. The effectiveness will depend upon the distance between source and receiver, the nature of the seismic sources and the level of regional seismicity.     

Identification of contaminant point source in surface waters based on backward location probability density function method

A backward location probability density function (BL-PDF) method capable of identifying location of point sources in surface waters is presented in this paper. The relation of forward location probability density function (FL-PDF) and backward location probability density, based on adjoint analysis, is validated using depth-averaged free-surface flow and mass transport models and several surface water test cases. The solutions of the backward location PDF transport equation agreed well to the forward location PDF computed using the pollutant concentration at the monitoring points. Using this relation and the distribution of the concentration detected at the monitoring points, an effective point source identification method is established. The numerical error of the backward location PDF simulation is found to be sensitive to the irregularity of the computational meshes, diffusivity, and velocity gradients. The performance of identification method is evaluated regarding the random error and number of observed values. In addition to hypothetical cases, a real case was studied to identify the source location where a dye tracer was instantaneously injected into a stream. The study indicated the proposed source identification method is effective, robust, and quite efficient in surface waters; the number of advection–diffusion equations needed to solve is equal to the number of observations.     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.

     

基于伴随方法的线性台阵背景噪声面波和远震体波联合成像研究

伴随层析成像（Adjoint Tomography）通过求解全波方程来准确模拟地震波在复杂介质中的传播,并利用波形信息来反演地下结构,是新一代的高分辨率成像方法.其中3-D伴随层析成像需要庞大的计算资源,而2-D反演相对更具计算效率.面波和远震体波是研究地壳上地幔速度结构的重要方法,它们对S波速度及Moho面的敏感度不同,通过联合反演,可以得到更为准确的S波速度结构及Moho面.通过两种数据的高度互补性,本文提出基于伴随方法的线性台阵背景噪声面波和远震体波联合成像方法,同时约束台阵下方S波速度结构及Moho面形态.我们将该方法应用到符合华北克拉通岩石圈典型结构特征的理论模型上,测试结果表明联合反演方法优势明显,相比于面波伴随层析成像,能获得更高分辨率的S波速度结构,同时能精准约束Moho面形态.相比于体波伴随层析成像,联合反演能有效压制高频假象,降低波形反演过程中的非线性化程度.本研究有望提供一种更为高效精准的线性台阵成像方法,搭建联合伴随层析成像理论框架,提升岩石圈成像分辨率,并为后续其他类型波形数据的引入提供思路和方法.     

基于重磁异常相关分析的场源位置及属性识别方法

针对传统的重磁对应分析在场源区域外趋于高相关度的缺点,本文提出了一种基于相关系数和垂向导数的重磁场源位置及属性识别方法,给出了一个新的重磁场源平面位置及属性判别参数（GMCP）,该参数能够有效地缩小重磁场源的识别范围.判别参数GMCP值非零分布范围反映了场源的规模.GMCP值的正负反映了场源属性,当GMCP大于零时为正相关,代表了该区域存在高磁高密度或低磁低密度的重磁同源地质体;当GMCP小于零时为负相关,代表了该区域存在高磁低密度或低磁高密度的重磁同源地质体;GMCP趋于零,则表示无重磁同源地质体.两组不同干扰水平的同源复杂组合模型实验、两组不同源组合模型试验以及南海中央海盆实际重磁数据处理结果均验证了该重磁场源位置及属性识别方法的正确性和有效性.

     

Body-wave waveform constraints on the source parameters of the Yangjiang,China, earthquake of July 25, 1969: A devastating earthquake in a stable continental region

On July 25, 1969, anM _s 5.9 earthquake shook the Kwangtung Province near Yangjiang, China. Casualties and extensive damage were reported in the epicentral area. The Yangjiang earthquake occurred within the passive margin located along China's southeast coast. This stable continental setting is seismically one of the most quiet regions of China; historic records indicate this earthquake to have been the first devastating one in the area. A remarkable feature of this earthquake sequence is that its foreshocks and aftershocks are relatively small in terms of number and magnitude despite the relatively large main shock.Waveforms of long-and short-period teleseismicP and long-period teleseismicSH waves have been modeled to estimate the source parameters for this earthquake. The focal mechanism determined is predominately strike-slip with a small normal component (strike=263°, dip=79°, rake=194°) and is in good agreement with observations in the epicentral area regarding dislocation direction of fissures, distribution of aftershocks, and isoseismals. The average seismic moment is 5.15×10²⁴ dyne-cm and the focal depth is estimated to be 9 km. The earthquake is characterized by a rather short source-time function and a high stress drop of approximately 380 bars. The observed pattern of the foreshock-main shock-aftershock sequence is interpreted in terms of the large stress drop associated with the main shock and the material properties at and near the source. TheP-axis orientation lies roughly NW-SE and is consistent with the maximum compressive stress observed along the coast of southeast China. This stress orientation is consistent with the NW subducting Philippine Sea plate and the SE ward push from central China due to the convergence between Indian and Eurasian subcontinents.     

Enhancement of real‐time hybrid simulation on a shaking table configuration with implementation of an advanced control method

This paper presents the implementation of Three Variable Control (TVC), an advanced control method, to the existing hybrid simulation (HS) system at the University of California, Berkeley. Motivation, background, and implementation of the TVC are explained together with modifications in the existing HS system. An application, which consists of the real‐time HS of electrical disconnect switches on a shaking table configuration, demonstrates successful implementation of the TVC. The presented application also covers other HS‐related features, namely employment of a three‐dimensional analytical substructure, real‐time HS‐compatible operator‐splitting integration method, and an efficient equation solver for faster computations. Copyright © 2014 John Wiley & Sons, Ltd.