2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

位场各阶垂向导数换算的新正则化方法

位场垂向导数大量应用于位场数据处理与解释中.当前广泛采用的位场各阶垂向导数换算方法为基于Laplace方程并结合波数域和空间域方法的具有递推特性的ISVD(integrated second vertical derivative)算法.本文在位场垂向导数换算的正则化方法和径向平均功率谱的基础上,提出一种位场各阶垂向导数换算的新正则化方法.新正则化方法仅需通过分析位场径向平均功率谱来确定一个截止波数,即可稳定换算位场各阶垂向导数.理论模型和实测数据实验结果表明:(1)新正则化方法物理意义明确、计算简单,且各阶垂向导数换算的稳定性和精度明显优于ISVD算法;(2)在用新正则化方法求得各阶垂向导数的基础上,利用泰勒级数法可以获得大深度、高精度的位场向下延拓结果.     

地球形成和演化过程中的分异能计算方法研究

地球形成初期,构成地球的物质在组成上是大致均一的.目前地球的地核-地幔-地壳圈层结构,是由分异作用形成的.分异过程释放的能量称为分异能.Sorokhtin和Chilingarian等人从行星吸积的定义出发,导出了基于地球内部密度分布的势能计算公式,计算出的分异能大小为1.698×1031J.本文采用计算球体势能的思路,导出分异能计算的解析公式和数值计算公式,通过求取原始地球模型与均匀分层模型、PREM模型的势能差计算分异能.两种方法的计算结果分别为1.535×1031J和1.698×1031J.前者与Sorokhtin等的结果相近,后者与之相同.本文初步分析了方法间的异同以及造成结果偏差的主要原因.     

利用加强水平方向总水平导数对位场全张量数据进行边界识别

位场全张量梯度数据以其信息量大、含有更高频的信号成分,能更好地描述小的异常特征等优点在地球物理领域中得到广泛应用.边界检测是位场解释中不可缺少的任务,需要新的边界探测器来处理位场梯度张量数据.为了充分利用位场梯度张量数据的多信息成分,本文定义了方向总水平导数和加强方向总水平导数,并利用其定义新的边界检测器.为了能同时显示不同振幅大小异常的边界,本文对其进行了归一化处理.通过模型试验,证明了归一化方法能更加清晰准确地显示浅部和深部的地质体边界信息.最后将该边界检测方法用于加拿大圣乔治湾实际测得全张量重力梯度数据和中国朱日和地区的磁异常数据中,并得到了较好的边界检测结果.     

Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

龙门山和成都地震构造区的划分

新编制的地震动参数区划图采用了潜在震源区三级划分方案,以体现背景地震活动空间分布的不均匀性,并在地震构造区内归纳出统一的地震构造模型.本文根据西南地区潜在震源区三级划分的成果,分析了龙门山地震统计区内的龙门山和成都地震构造区的基本特征,历史地震活动强度及频度,主要活动构造的构造变形样式,建立了地震构造区的发震构造模型,确定了构造区的本底地震及划分构造源的地震构造标志.同时,提出了确定背景源空间分布函数的简单方法.     

敦煌莫高窟潜在地震变形数值模拟

在研究莫高窟地区地震地质条件及潜在地震危险性分析的基础上,基于地震的断层弹性位错理论,采用位移水平分布的抛物线模型,假定三危山断层分别发生5.5、6.0和6.5级地震情况下,用三维有限元的方法分析计算了莫高窟地区的潜在地震变形。给出了在震级5.5、6.0及6.5级三种情况断层逆冲和水平走滑两种活动方式下,莫高窟区的地震位移和应变值,供石窟的抗震设防加固设计参考。     

On-chip earthquake simulation model using potentials

A two-dimensional (2-D) Cellular Automata (CA) dynamic system constituted of cells-charges has been proposed for the simulation of the earthquake process. The CA model has been calibrated with the use of real data. The calibration incorporates major seismic characteristics of the area under test. The simulation results are found in good quantitative and qualitative agreement with the recorded Gutenberg–Richter (GR) scaling relations. The model is enriched with a powerful multi-parameter interface that enables the user to load real data from different regions. This paper examines the on-chip realisation of the model and its instrumentation. The CA model hardware implementation is based on Field Programmable Gate Array (FPGA) logic. It utilises an array of 32 × 32 cells. Parameters that construct the local CA rule constitute the input data. The initial seed, which to some extent corresponds to the seismic features of the area under test, is loaded in a semi-parallel way and the process is completed in a certain number of time steps. The automatic response of the processor provides the corresponding GR scaling law of the area under study. The hardware implementation of the CA-based earthquake simulation model is advantageous in terms of low-cost, high-speed, compactness and portability features. It can operate as a preliminary data-acquisition filter that accelerates the evaluation of recorded data as far as its origin time, spatial and magnitude completeness and quality are concerned. Software that performs reliable automatic phase picking, as well as data elaboration, can be assembled next to the earthquake recording instruments (the whole network) output to assure a quick and reliable iteration of the focal parameters of a recorded earthquake (epicentre coordinates, focal depth and magnitude). The dedicated processor can be accommodated right after this stage (before any manual elaboration) focusing on the near real-time development of a reliable qualitative dynamical seismic record and a mapping of the seismic characteristics of the area.     

海南省大广坝水电站集水区降水特征及人工增雨潜力分析

本文分析了海南省大广坝水电站集水区降水分布特征并研究了人工增雨的潜力。结果表明,大广坝水电站集水区降水年际差异较大,常年无法满负荷发电,影响了水电站效益的正常发挥。4～7月集水区实施人工增雨的潜力和成云致雨的条件较好,通过人工增雨的手段可以有效提高水电站的发电效益。     

海水ORP对3C钢腐蚀行为的影响

研究了氧化还原电位(Oxidation Reduction Potential, ORP)对海水中碳钢腐蚀电化学特征的影响,应用动电位极化曲线法探讨了 ORP 与碳钢腐蚀行为之间的关系。结果表明,海水的 ORP 参数对碳钢腐蚀行为具有较大的影响,其中 O2/OH?电对的作用占主导,其他电对如 Fe3+/Fe2+也有一定的影响。随着海水 ORP 的增大,其氧化能力越强,碳钢腐蚀电流增大,腐蚀速度加快。水质相对稳定的开放性大洋海水中,可以考虑海水氧化还原电位在碳钢腐蚀速度预测评估中的应用。