加权最小范数平差的伪观测解法

首先证明了加权最小范数平差一般解的一个重要性质，由此提出了利用附加的观测值来求解加权最小范数平差问题的方法，并证明了该法与利用一般解法等价性的充要条件。     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.     

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.     

用非线性遗传反演方法进行自动静校正

静校正处理是消除近地层影响的一个重要地震资料数字处理方法，在近地表层纵，横向变化复杂，静校正量变化较大的地区，常规静校正方法和线性反演方法均不能解决问题，遗传算法是近年发展起来的一种模拟生物生化过程的较先进的非线性反演方法，对用遗传算法进行静校正的若干问题进行了讨论，说明它是可以解决表层情况复杂地区静校正问题的一个有效方法。     

混合岩中斜长石的交代净边结构和倒转双晶研究——以大别罗田黄土岭长英片麻岩为例

湖北罗田黄土岭长英片麻岩是与紫苏石榴黑云片麻岩共生的浅色片麻岩，成分上近于Ａ型花岗岩．受到浅混合岩化，钾交代现象显著．在斜长石与钾长石新成体接触处有清楚的交代净边结构，背散射电子探针分析表明，净边成分为钠长石（Ａｂ８９～９７）．在净边结构上，发现叠加倒转双晶．净边结构产生于长英片麻岩部分熔融的高温条件下，与整个大别地体隆升过程和Ａ型花岗岩的生成有成因联系     

论广义逆矩阵阶数增大的递推算法

本文针对加拿大地球物理学家ＥｄｏＮｇｌａｎｄ教授提出的问题，对Ｍ—Ｐ广义逆矩阵的递推算法进行了探分论述探讨的初步结果和对Ｇｒｅｖｉｌｌｅ递推算法的扩充。文中对算法的扩充给出严格的证明。     

电法勘探的发展和展望

电法勘探在经历了近一个世纪的发展后，其方法理论、仪器设备、野外数据采集、处理和解释等方面都经历了一系列重大变化．本文以方法理论的进展为主线，回顾、展望了目前电法勘探中几个重要而令人关注的研究焦点．这些问题的研究进展将会对21世纪的电法勘探产生深远的影响．