Time-domain electromagnetic migration in the solution of inverse problems

Time-domain electromagnetic (TDEM) migration is based on downward extrapolation of the observed field in reverse time. In fact, the migrated EM field is the solution of the boundary-value problem for the adjoint Maxwell's equations. The important question is how this imaging technique can be related to the solution of the geoelectrical inverse problem. In this paper we introduce a new formulation of the inverse problem, based on the minimization of the residual-field energy flow through the surface or profile of observations. We demonstrate that TDEM migration can be interpreted as the first step in the solution of this specially formulated TDEM inverse problem. However, in many practical situations this first step produces a very efficient approximation to the geoelectrical model, which makes electromagnetic migration so attractive for practical applications. We demonstrate the effectiveness of this approach in inverting synthetic and practical TDEM data.     

Comment on 'Estimation of resolution and covariance for large matrix inversions' by J. Zhang and G. A. McMechan

When inverting large matrices, iterative techniques are necessary because of their speed and low memory requirements, as opposed to singular value decomposition (SVD). Recently, there have been attempts to obtain information on the quality of the solutions calculated using conjugate gradient (CG) methods such as LSQR. The purpose of this note is to comment on the paper titled "Estimation of resolution and covariance for large matrix inversions' by Zhang & McMechan (1995), who extend Paige and Saunders' LSQR algorithm to obtain an orthonormal basis used to approximate resolution and covariance. We show that for larger problems, where the number of orthogonal vectors is several orders of magnitude smaller than the number of model parameters, the vectors obtained do not adequately span the range of the model space. We use a synthetic borehole experiment to illustrate the differences between the singular value spectrum obtained through the more complete method of SVD and the Ritz value spectrum that results from a simple extension of LSQR, We also present a trivial numerical example to illustrate the differences between Zhang & McMechan's approximate resolution matrix and the true resolution.     

黄汲清院士与第四纪冰川研究

黄汲清院士是在多方面作出重大贡献的地质学大师，也是第四纪冰川研究的先驱者。他于１９３２￣１９３５年在瑞士完成博士论文期间，受到良好的冰川地质学训练。于１９４３￣１９４４年与合作者对天山南麓台兰河谷中游地质调查中，对冰碛与非冰碛沉积物进行仔细的正确的调查和精美的地质制图，撰写出高水平的论文，由此对中国第四纪冰川研究作出了重要贡献。他关注中国东部第四纪冰川研究的进展，指出忽视古气候研究是工作中的重要缺     

三维溶质运移问题的分步广义迎风解法

对对流占优的三维溶质运移问题提出了分步广义迎风解法，首先利用Ｎ．Ｎ．Ｙａ－ｎｅｎｋｏ对水动力弥散方程分步求解的思想，将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题，对对流定解问题采用广义迎风对偶单元均衡法求解，对扩散定解问题采用一般的Ｇａｌｅｒｋｉｎ有限元法求解，不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题，而且避免了求节点速度这一步，节省运算步骤，对井点的浓度变化给出了更合适的求解方法。     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解，在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式，从而为研究物理大地测量中边值问题的求解提供了新的依据     

3-D rheologic model of earthquake preparation (II): Strain field and its applications

Introduction Based on the elastic theory of the hard inclusion (Dobrovolskii, 1991), we developed an inclusion theory of rheologic medium, and applied the results of bulk-strain field of a rheologic inclusion model to explain the spatial-temporal evolution process of earthquake precursors (SONG, et al, 2000). In the former paper (SONG, et al, 2003), we derived the viscoelastic displacement field of the rheologic inclusion model on the basis of the analytic expression of displacement field o…     

花岗岩研究进展 ———ISMA花岗岩类分类综述

本文在论述花岗岩主要分类的基础上论述了ISMA花岗岩类的基本含义,讨论了人们在以后应用和推广ISMA这类花岗岩分类方案时存在的一些问题,指出了ISMA花岗岩类本身在分类过程中存在的一些基本问题。     

地热井常见主要问题分析与研究

针对目前地热井施工技术现状及地热井使用中常见的问题，进行了较为详细的分析和研究，并就今后发展方向和建议进行了初步探讨。