我国西北部分地区地壳-上地幔电性横向变化特征与地震活动的关系

由50余个测点的大地电磁测深资料,讨论了该区的地壳-上地幔电性横向变化特征。按照上地幔第一低阻层顶面埋深,将测区划分为四类地区:浅埋深(55-90公里)、中浅埋深(90-110公里)、中深埋深(110-160公里)和深埋深(160-250公里)。讨论了本区六次大震例的深部电性背景。指出了上地幔顶部的梯度带地区、地壳内电性横向变化剧烈的地区和地壳内存在局部特殊增厚的低阻层地区将可以作为潜在震源区的深部电性判据之一     

波动方程的高阶广义屏叠前深度偏移

不同于常规广义屏传播算子的推导中使用散射理论,本文利用单平方根算子的渐近展开,推导出了单程波方程广义屏传播算子的高阶表达式.高阶广义屏传播算子不仅可提高常规广义屏传播算子的计算精度,而且还能改善广义屏传播算子对速度强横向变化介质的适应性.把高阶广义屏传播算子应用于波动方程叠前深度偏移,可得到比常规广义屏传播算子更好的效果.高阶广义屏传播算子的阶数越高,计算精度越高,但计算量也越多.以SEG EAGE二维盐丘模型数据的波动方程叠前深度偏移为例,二阶广义屏传播算子相对于常规（一阶）广义屏传播算子增加了30%的计算量.高阶广义屏传播算子是常规广义屏传播算子理论的发展和完善.     

Comparison between the Fourier finite-difference method and the generalized-screen method

The Fourier finite‐difference propagator and the generalized‐screen propagator are two general high‐order forms of one‐way dual‐domain methods. We compare these two propagators mainly on phase accuracy, computational efficiency and 3D extension. A comparison of phase accuracy shows that the high‐order generalized‐screen propagator is preferable to the Fourier finite‐difference propagator for heterogeneous media with a weak velocity contrast and wide dip angle. With increasing velocity contrast, the accuracy improvement gained by the high‐order generalized‐screen propagator declines rapidly. The Fourier finite‐difference propagator is more robust and flexible to lateral velocity variations than the generalized‐screen propagator. The 2D Fourier finite‐difference propagator is superior to the 2D generalized‐screen propagator when the velocity contrast is stronger than 23%. Despite the two‐way splitting error, the 3D Fourier finite‐difference propagator is more accurate than the second‐order generalized‐screen propagator when the velocity contrast is stronger than 20% and is more accurate than the fourth‐order generalized‐screen propagator when the velocity contrast is stronger than 40%. Numerical experiments on the SEG/EAGE salt model demonstrate that the Fourier finite‐difference propagator behaves better than the generalized‐screen propagator when imaging steep salt boundary and faults beneath the salt body. Under the same hardware and software conditions, the computational cost of the Fourier finite‐difference propagator in our implementation is greater than that of the second‐order generalized‐screen propagator but smaller than that of the third‐order generalized‐screen propagator. Compared with the Fourier finite‐difference propagator, the generalized‐screen propagator requires fewer grid points per wavelength and has more potential to improve running speed in the presence of a much faster Fourier transform. These analyses are applicable for both forward modelling and depth migration.     

闾阳—梨树断面深,浅部构造特征及其与地震活动的关系

本文利用壳－幔波速,电性等深部地球物理信息,结合区域地质及实验岩石学研究结果等资料,对闾阳－梨树断面进行了地质推断解释,并简略讨论了壳－幔结构对浅部构造运动的影响及与地震活动的关系。     

Review of the Anisotropic Interface Ray Propagator: Symplecticity,Eigenvalues, Invariants and Applications

A review of the 6 × 6 anisotropic interface ray propagator matrix in Cartesian coordinates and within the framework of the Hamiltonian formalism shows that there is one unique propagator satisfying the symplectic property. This is essential, since the symplecticity furnishes an exact inverse, while an eigenvalue analysis indicates that the propagator may be arbitrarily ill-conditioned. As such, the symplectic interface propagator naturally connects to symplectic ray integration algorithms for smooth media, designed to maintain accuracy. Moreover, several ray invariants for smooth media remain invariant across interfaces. It is straightforward to derive expressions for the interface propagator, both explicit and implicit. Symplecticity is equivalent to the condition that the propagator preserves the eikonal constraint across the interface. The symplectic interface propagator complies with phase matching of the incident and reflected/transmitted ray field, and is therefore in accordance with the earlier derived 4 × 4 matrix in ray-centred coordinates. The symplectic property is related to the symmetry of the second derivative matrix of the reflected/transmitted traveltime field. Thanks to the analytic expression of the symplectic interface propagator, relating interface curvature directly to second derivatives of traveltimes observed at a datum level, numerous applications are available in the area of processing and inversion.     

面向目标的小束源照明和成像

结合小束源和Fourier传播子,应用Fourier传播子进行面向目标的小束源照明和成像.小束源的合成通过小波束变换中的小束函数获得,而Fourier传播子进行波场外推,完成小束源的波传播.小束源具有空间位置和方向的双局部特性,对于照明和成像,它更加灵活,并有更多的控制方法.通过照明分析,选择面向目标的有效照明小束源,进行部分小束源偏移,可以提供更好的成像质量和计算效率.作为数值试验,我们选用Fourier有限差分传播子,对 Marmousi模型和二维SEG-EAGE盐丘模型的数据,试验了小束源对目标结构的方向照明和成像.获得的结果表明,结合小束源和Fourier传播子进行面向目标的照明和成像是切实有效的.     

局部角度域波传播步进算法研究

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.     

三维非均匀地幔中深源电场和磁场的谱理论(Ⅱ)——横向缓变摄动和极型场快速扩散

It is one of the classical problems in Earth′s electromagnetism that continuation of the observatory data of magnetic field into the conducting region. This paper built the coupling vector equations governing the poloidal, toroidal and potential fields in 3-D inhomogeneous conducting mantle. Considering the limitation on variable scale of inhomogeneities in global mantle from study of the Earth′s deep interior, a perturbation　theory of gradually lateral variation was presented. It is unnecessary that the 1-D spherical symmetric distribution of electric conductivity as zero-degree approximation. Serving as a example of solvability of the zero-degree approximation, it was demonstrated how the gradually lateral variation of electric conductivity effect on anti-diffusive problem of poloidal field in Earth′s mantel.     

三维非均匀地幔中深源电场和磁场的谱理论（Ⅱ）──-横向缓变摄动和极型场快速扩散

把磁场观测数据向导电区域延拓,是地球电磁学的经典问题之一．从准静态近似和非零矢势规范的电磁场方程出发,依据矢量的球面分解唯一性定理,本文建立支配三维非均匀电寻率分布全球地慢中的环型场、极型场和电位势场的耦合方程组烤虑地球深部研究对认识全球地幔非均匀性横向变化尺度的限制,提出横向缓变意义下三维非均匀地幔中电场和磁场的摄动理论,其零级近似不要求电导率分布一维球对称．作为零级近似可解的例证,研究了地幔深源极型场的反扩散问题,在利用地面磁场观测反演到核幔边界时,可以考虑地幄横向非均匀性影响．     

Thermal field and electric conductors in the crust and upper mantle of Ukraine

The paper presents data on the distribution of the deep heat flow in Ukraine, thermal models of the tectonosphere, and electric conductors in the crust and mantle that are mostly related to partial melting and dehydration of rocks.     

The state of the lithosphere in the junction zone of Tarim and Tien Shan according to the petrological interpretation of the magnetotelluric data

The geological-geophysical and petrological-geochemical studies of the Earth’s crust and upper mantle are combined to estimate the state of the lithosphere at the junction zone of Tarim and Tien Shan. The laboratory measurements of electric conductivity in the rocks sampled from the upper mantle and lower crust considered against the geoelectrical and thermal models revealed lherzolite, granulite, and eclogite massifs in the deep section of the Tarim and Tien Shan junction zone. The experimental results suggest that the crustal thickness in the southern Tien Shan attained 35–40 km 70 Ma ago.     

层状介质中地震面波频散函数和体波广义反射系数的计算

本文在B、P、C坐标中给出了弹性波在横向均匀介质中的传播矩阵,并将其表示为五个形式简单的矩阵乘积,其中有四个矩阵是与频率无关的。我们用传播矩阵的分离方式,将Abo-Zena（1979）算法大量简化。用本文的算法计算综合地震图时,计算量要比文献[2,7,8]的算法少一半以上。 文中还给出了地震面波频散函数和体波广义反射系数的快速计算步骤。     

中国南北地震带北段地壳和上地幔的电性特征

本文介绍了在我国南北地震带北段地区进行大地电磁测深所得到的一些非常初步的结果,主要是这一地区内的地壳和上地幔的电性特征以及震区和非震区的电性差异。     

Lateral velocity variation related correction in asymptotic true‐amplitude one‐way propagators

It was mathematically proved that the asymptotic true‐amplitude one‐way wave equation could provide the same amplitude as the full‐wave equation in heterogeneous lossless media in the sense of high‐frequency asymptotics. Much work has been done on the vertical velocity variation related amplitude correction term but the lateral velocity variation related term has not received much attention, even being excluded in some asymptotic true‐amplitude one‐way propagator formulations. Here we analyse the effects of different amplitude correction terms in the asymptotic true‐amplitude one‐way propagator, especially the effect related to the lateral velocity variation, by comparing the wavefield amplitude from the one‐way propagator with that from full‐wave modelling. We derive a dual‐domain wide‐angle screen type asymptotic true‐amplitude one‐way propagator and evaluate two implementations of the amplitude correction. Numerical examples show that the lateral velocity variation related correction term can play a significant role in the asymptotic true‐amplitude one‐way propagator. Optimization of the expansion coefficients in the asymptotic true‐amplitude one‐way propagator can improve both the amplitude and phase accuracy for wide‐angle waves.     

基于耦合反射/透射系数单程波传播算子的地震波模拟研究

单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合.     

长江中下游成矿带中段岩石圈电性结构研究

长江中下游成矿带位于大别造山带、长江中下游凹陷、江南隆起带等大地构造单元结合部位,通过在研究区内布设两条首尾相接共计150km长的大地电磁剖面,获得了50km以浅岩石圈尺度的电性分布.长江中下游地区中段地下电性结构显示出在地下10km和30km处分别存在明显的圈层结构,以此认为现今横向稳定的"电莫霍"反映了研究区经历燕山期陆内构造-岩浆活动后已基本上完成壳幔重新平衡;而分隔大地构造单元的郯庐断裂带、长江断裂带以及江南断裂带在电性上具有特征的梯度显现,在印支造山期后的引张背景下,断裂带成为强伸展活动带与控制了燕山期大范围的陆内岩浆活动;高导地幔的局域性存在以及从北向南地幔导电性的变化反映了在经受深部动力学过程中处于不同大地构造部位的地幔所遭受的不同类型的改造以及地幔深部的构造极性.     

Synthetic radargrams from electrical conductivity and magnetic permeability variations

The objective of this work is to assess the importance of electrical conductivity and magnetic permeability variations in ground penetrating radar (GPR) reflections commonly interpreted only in terms of permittivity variations. We use the matrix propagator approach to obtain the surface electric field associated with a horizontally layered model whose three electromagnetic properties vary from layer to layer. The solution is based on the plane wave boundary value problem using inverse Fourier transformation to accommodate particular GPR pulses. Our results indicate that while magnetic permeability is unimportant, reflections from electrical conductivity variations can be of the same order as those associated with electrical permittivity boundaries. In particular, we show that a realistic ground model composed of thin conductive layers can produce radargrams similar to those caused by a lossless permittivity contrast.     

滇西地区壳幔解耦与腾冲火山区岩浆活动的深部构造研究

根据青藏东部边缘的深部地球物理资料,分析了滇西地区壳幔耦合和腾冲火山区岩浆活动的深部构造特征,确认了地幔各向异性与上地幔速度结构(包括P波速度和S波速度)的内在联系,指出产生这一结果的原因与以腾冲火山区为中心的地幔热物质上涌有关:上地幔顶部平均温度升高导致介质强度降低,在印支块体的侧向挤压或印缅块体的向东俯冲作用下发生韧性变形,造成滇西地区地幔各向异性的快波方向与青藏东部地壳块体的旋转方向不一致.此外,鉴于中下地壳低速层的横向非均匀性,估计韧性流动并非贯通青藏高原的东部边缘,而是被不同的构造块体和边界断裂限定在局部地区.总体而言,滇西地区下地壳的地震波速度和电阻率偏低,具备发生韧性变形的构造条件.作为地壳和上地幔之间的解耦层,它使得青藏东部地壳块体旋转产生的构造应力未能传输至上地幔.腾冲火山区的地壳结构与不同时期的岩浆活动有关,火山区东侧的高速结构代表了上新世时期火山通道内冷凝固结的岩浆侵入体或难以挥发的高密度残留物质,火山区西侧的低速结构反映了更新世以来持续至今的岩浆活动,壳内岩浆源主要分布在10～20km的深度范围内,横向尺度约为15～20km,有可能通过地壳深部的断裂与上地幔岩浆源区相连,估计腾冲火山区下方的岩浆活动将持续进行.     

Electromagnetic excitation of asthenospheric conducting zones

The relations between galvanic and inductive electromagnetic excitation of 3-D asthenospheric conducting zones are considered. Mathematical modeling with the application of integral equations shows that the rather intense magnetotelluric anomalies reflecting the topography of the conducting asthenosphere arise via the galvanic channel, which has two branches: (1) the lithospheric branch (an electric current flows from the sedimentary cover into the lithosphere and concentrates in uplifts of the conducting asthenosphere), which is inactive in the case of a very high resistivity of the lithosphere (the screening effect), and (2) the mantle branch (the deep mantle current flows from below into uplifts of the conducting asthenosphere), which is free from the screening effect of a high resistivity lithosphere. In previous works, the mantle branch of the galvanic channel was erroneously associated with local induction in uplifts of the asthenosphere. According to modern estimates, local induction in the asthenosphere generates weak magnetotelluric anomalies that cannot be reliably measured.     

Electromagnetic study of lithospheric structure in the marginal zone of East European Craton in NW Poland

The marginal zone of the East European Platform, an area of key importance for our understanding of the geotectonic history of Europe, has been a challenge for geophysicists for many years. The basic research method is seismic survey, but many important data on physical properties and structure of the lithosphere may also be provided by the electromagnetic methods. In this paper, results of deep basement study by electromagnetic methods performed in Poland since the mid-1960s are presented. Over this time, several hundred long-period soundings have been executed providing an assessment of the electric conductivity distribution in the crust and upper mantle. Numerous 1D, 2D, and pseudo-3D electric conductivity models were constructed, and a new interpretation method based on Horizontal Magnetic Tensor analysis has been applied recently. The results show that the contact zone is of lithospheric discontinuity character and there are distinct differences in geoelectric structures between the Precambrian Platform, transitional zone (TESZ), and the Paleozoic Platform. The wide-spread conducting complexes in the crust with integral conductivity values reaching 10 000 S at 20–30 km depths are most spectacular. They are most likely consequences of geological processes related to Caledonian and Variscan orogenesis. The upper mantle conductivity is also variable, the thickness of high-resistive lithospheric plates ranging from 120–140 km under the Paleozoic Platform to 220–240 km under the East European Platform.