提高任意广角波动方程计算效率及压制倏逝波干扰方法

任意广角波动方程（Arbitrarily Wide-angle Wave Equation,AWWE）是一种高精度的时空域单程波动方程,与其他时空域单程波方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文改进了现有的AWWE有限差分计算方案,避免了原计算方案中对矩阵反复求逆的过程,同时减少了矩阵相乘的次数,显著地提高了AWWE数值计算的效率.数值实验证实,采用改进之后的计算方案,其计算结果与原方案完全一致.此外,本文采用了一种方法,在有倏逝波干扰的情况下,能在频率-波数（f-k）域区分出倏逝波与非倏逝波,为设计f-k滤波器提供合理的视速度门槛信息.偏移实验证实,在深度延拓时仅仅进行数次f-k滤波就能很好地抑制倏逝波干扰,提高偏移成像质量.     

基于选权拟合法的电离层电子密度层析重构

附加约束的电离层层析算法是解决电离层电子密度反演中不适定问题的主要方法,为避免此类方法中约束权阵的选取不当对电子密度分布重构产生的不良影响,本文将选权拟合法应用到电离层层析成像技术中.该方法特别设计了依据电子密度空间分布特性构造参数权矩阵的方案.新方法有明确的物理意义,挖掘了隐含的信息量,为解决电离层电子密度反演中由于观测数据的不足等因素引起的不适定问题提供了一种新途径,可以得到符合客观实际的结果.数值模拟实验和实测数据的反演结果证实了该算法有效性、可靠性和优越性.     

中国西部地区地壳结构特征与强震活动相关性研究

中国西部地区是地震活动十分强烈的地区,天山、阿尔泰、帕米尔和西昆仑都是著名的地震构造带,在这些地震构造带和周边地区发生了多次震级大于5级的强震.本文通过分析西部地区的重力场特征,根据重力数据结合地震剖面、应用Parker-Oldenburg方法反演得到了研究区莫霍面深度,通过对比地震层析成像的反演结果,分析了研究区的地壳结构特征.计算结果表明,研究区地壳结构不均匀特征明显,在造山带地区一般是莫霍面坳陷区,盆地则是莫霍面隆起区,主要造山带地壳速度结构表现为高速区,盆地和主要凹陷区为低速区.根据计算结果和以往强震震中位置分析了地壳构造与强震活动的相关性,西部地区的地震活动与地壳结构的横向不均匀密切相关,强震主要发生在地壳速度变化带附近和地壳速度结构差异较大的地区,在构造应力作用下,这些地壳介质非均匀地区易发生强震,这是中国西部造山带和盆-山边界附近频发强震的构造原因之一.     

岩体间超低摩擦发生机理的摆型波理论

针对深部冲击地压发生时煤岩体有明显的错动现象,基于块系岩体中的摆型波传播理论,研究摆型波在块系岩体中传播时诱发岩块间的超低摩擦发生机理。通过分析动力作用下岩块间的相对位移最大值,研究块系岩体的超低摩擦发生规律。当岩块间的软弱介质处于拉伸最大值状态时,在侧向扰动作用下岩体极易发生超低摩擦滑动,由此给出了块系岩体的超低摩擦发生判据。最后,通过计算分析了块系岩体自身特性及外界扰动特性对岩块间超低摩擦的影响。基于摆型波理论的超低摩擦发生机理研究,对认识深部岩体错动型冲击地压启动规律具有重要意义。     

New evidence of delamination in the Western Alboran Sea: Geodynamic evolution of the Alboran domain and its margins

The presence of continuous upper crustal blocks between the Iberian Betics and Moroccan Rif in the western and middle Alboran Sea, detected with tomography, can add new information about the lithosphere structure and geodynamic evolution in this region. A large volume of seismic data (P and S wave arrival times) has been collected for the period between 1 December 1988 and 31 December 2008 by 57 stations located in northern Morocco (National Institute of Geophysics, CNRST, Rabat), southern Portugal (Instituto de Meteorologia, Lisbon) and Spain (Instituto Geografico National, Madrid) and used to investigate the lithosphere in the western Alboran Sea region. We use a linearized inversion procedure comprising two steps: (1) finding the minimal 1-D model and simultaneous relocation of hypocenters and (2) determination of local velocity structure using linearized inversion. The model parameterization in this method assumes a continuous velocity field. The resolution tests indicate that the calculated images give near true structure imaged at 5 km depth for the Tanger peninsula, the Alhoceima region and southern Spain. At 15, 30 and 45 km depth we observe a near true structure imaged in northern Morocco, and southern Spain. At 60 and 100 km, southern Spain and the SW region of the Alboran Sea give a near true structure. The resulting tomographic image shows the presence of two upper crustal bodies (velocity 6.5 km/s) at 5–10 km depth between the Betics, Rif, western and central Alboran Sea. Low velocities at the base of these two bodies favor the presence of melt. This new evidence proves that the Tethysian ocean upper crust was not totally collapsed or broken down during the late Oligocene–early Miocene. These two blocks of upper crust were initially one block. The geodynamic process in the eastern of the Mediterranean is driven by slab rollback. The delamination process of the lithospheric mantle terminates with the proposed slab rollback in the western part of the Mediterranean. This can be explained by the removal of the major part of the lithosphere beneath the area, except in the SW part of the Alboran Sea where a small part of the lithospheric mantle is still attached and is extends and dips to SE beneath the Rif, slowly peeled back to the west. A second detached lithospheric mantle is located and extends to eastern part of the Rif and dips to the SE. The removal of lithosphere mantle from the base of the crust was replaced and heated by extrusion of asthenospheric material coming from depth to replace the part of crust detached. A combination of isostatic surface/topographic uplift and erosion induced a rapid exhumation and cooling of deep crustal rocks.     

Detailed crustal structure in the area of the southern Apennines–Calabrian Arc border from local earthquake tomography

We present a new seismic velocity model for the southern Apennines–Calabrian Arc border region with the aim to better define the crustal structures at the northern edge of the Ionian subduction zone. This sector also includes the Pollino Mts. area, where a seismic sequence of thousands of small to moderate earthquakes has been recorded between spring 2010 and 2013. In this sector a seismic gap was previously hypothesized by paleoseismological evidences associated with the lack of major earthquakes in historical catalogs.To perform the tomographic inversion we selected ca. 3600 earthquakes that have occurred in the last thirty years and recorded by permanent and temporary networks managed by INGV and Calabria University. Using for the first time the Local Tomography Software for passive tomography inversion (LOTOS hereinafter) to crustal analysis in southern Italy, we have computed the distribution of Vp, Vs, and the Vp/Vs ratio. The obtained velocity model, jointly evaluated with results of synthetic modeling, as well as with the hypocenter distribution and geological information, gives us new constraints on the geodynamical and structural knowledge of the study area.The comparison between the shallow tomography sections and surface geology shows good correlation between velocity patterns and the main geological features of the study area. In the upper crust a low-velocity anomaly of P- and S-waves is detectable beneath the Pollino Mts. area and seems to separate the Calabrian and southern Apennines domains, characterized by higher velocities. The distributions of high Vp/Vs ratio, representing strongly fractured rocks with likely high fluid content, clearly correlate with areas of significant seismicity.In the lower crust we detect a clear transition from high to low seismic velocities in correspondence with the Tyrrhenian coast of the study area, which may represent the transition from the thinner Tyrrhenian crust to the thicker one beneath Calabria. In this area, also characterized by a progressive detachment of a retreating lithospheric slab, the generation of a Subduction-Transform Edge Propagator (STEP) fault zone, that laterally decouples subducting lithosphere from non-subducting lithosphere in a scissor type of fashion, may have taken place. These conditions imply the existence of a kinematic decoupling which allows for differential movement between the Calabrian Arc and the southern Apennine chain. The low velocity anomaly separating the southern Apennines and the Calabrian Arc domain may be related to fluid upwelling occurring in correspondence with the northern edge of the Calabrian subducting slab.     

Isolation of plane shear wave using water saturated trench barrier

The article presents a theoretical investigation of the propagation of normally and obliquely incident plane shear waves past a rectangular trench filled up with water saturated soil sandwiched between anisotropic elastic media. The motivation for this work is due to the effectiveness of the isolation of waves by the in-filled rectangular trench. Unlike the most of the previous researchers, this model considers the soil–structure interaction effects and directly determines the influence of barrier in the form of in-filled trench on the mode of wave propagation. It is of interest to determine the reflection and transmission coefficients, and the energy partition distribution of shear waves in the in-filled rectangular trench showing the influence of barrier on the propagation of waves. An extensive parametric study through numerical computation is carried out to investigate the influence of the material properties of the in-filled trench and the amplitude ratios on shear waves. The in-filled trench barrier directly declines the intensity of waves significantly in such a way that the waves do not create any hazards to the nearby structures, if exists at all.     

断层带震相及其应用研究现状与展望

Deep structure and material properties of faults can be understood by observing and simulating the particular phase in a fault fracture zone. This paper reviews the development of fault-zone seismic waves in the seismological domain. The present research status of fault-zone head wave and trapped wave are summarized systematically. Based on recent progress in this field, the paper discusses the prospect on the utilization of seismic wave in fault structure research.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

3D Modeling of sediment movement by ships-generated wakes in confined shipping channel

Ship-generated waves and return currents are capable of re-suspending significant quantities of bottom and bank sediments.However,most of the previous studies done on the subject do not show how and where sediment is re-suspended by the wakes and the directions of net transport.In this paper,a 3D numerical model based on hydro-sedimentary coupling is presented to search the relationship between the sediment movement,and the pattern of ship-generated waves around and far away from the vessel and the return currents around the ships.The hydrodynamic model is based on 3D Navier-Stokes equations including the standard k-ε model for turbulence processes,and the sediment transport model is based on a 3D equation for the re-suspended sediment transport.The computation results show that the areas of sediment concentration and transport(whether by resuspension or by the bedload) depend mainly on the position,the speed of the ship in the waterways,the kinematics of ship-generated waves and on the return flows.Thus,a map of sediment distribution and the modes of sediment transport generated by the passage of the ship are presented.