Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

频率域位场快速正则化反演界面和物性分布

从理论与模型算例两方面对一种以Ｐａｒｋｅｒ公式为基础的频率域位场反演方法做了进一步深入研究。表明它是一种能够更有效地引入低通滤波因子的迭代方法。与改进的正则化稳定因子相结合，提出新的快速正则化反演方法。除了将其用于反演求解均质界面以外，它也被推广至变物性的界面反演以及物性反演     

含油气盆地构造样式研究中几个问题的讨论

对我国含油气盆地构造样式研究中存在的一些主要问题进行了讨论。包括：关于构造样式的确切含义及类型的划分，由于对构造样式概念理解的不一致，导致了构造样式类型划分的混乱；构造样式的鉴别，往往由于过分强调了某些单一“鉴定标志”而忽略了总体变形特征，从而导致对构造变形机理的错误判断；中国含油气盆地的构造反转问题。     

The effects of pre-existing normal faults on thrust ramp development: An example from the northern Apennines,Italy

The Umbria-Marche foreland fold-and-thrust belt in the northern Apennines of Italy provides excellent evidence to test the hypothesis of synsedimentary-structural control on thrust ramp development. This orogenic belt consists of platform and pelagic carbonates, Late Triassic to Miocene in age, whose deposition was controlled by significant synsedimentary extension. Normal faulting, mainly active from Jurassic through Late Cretaceous-Paleogene time, resulted in significant lateral thickness variability within the related stratigraphic sequences. By Late Miocene time the sedimentary cover was detached from the underlying basement and was deformed by east-verging folds and west-dipping thrusts. Two restored balanced cross sections through the southernmost part of the belt show a coincidence between the early synsedimentary normal faults and the late thrust fault ramps. These evidences suggest that synsedimentary tectonic structures, such as faults and the related lithological lateral changes, can be regarded as mechanically important controlling factors in the process of thrust ramp development during positive tectonic inversion processes.     

The transition from passive to active margin tectonics: a case study from the zone of Samedan (eastern Switzerland)

The Zone of Samedan is part of a fossil, early Mesozoic rift system originally situated in the distal, Lower Austro-Alpine domain of the Adriatic passive continental margin. An early Mesozoic configuration of asymmetrical rift basins bounded by relative structural highs compartmentalized Late Cretaceous active margin tectonics; Jurassic half-grabens were folded into arcuate synclines, whereas relative structural highs engendered thin, imbricated thrust sheets. West-directed thrusting and folding initiated at the surface and continued to depths favoring mylonitization under lower greenschist-facies conditions. At this time Liguria-Piemontese ophiolites were accreted to Lower Austro-Alpine units directly underlying the Zone of Samedan. Late Cretaceous orogenic collapse of the Adriatic active margin involved the reactivation of west-directed thrusts as low-angle, top-to-the-east, normal faults. These faults accommodated extensional uplift of Liguria-Piemontese ophiolites and Lower Austro-Alpine units beneath and within the Zone of Samedan. During Paleogene collision, some Late Cretaceous faults in the Zone of Samedan were reactivated under lower anchizonal conditions as north-directed thrusts. The latter stages of this early Tertiary thickening were transitional to brittle, high-angle normal faulting associated with top-to-the-east extension and spreading above the warm, uplifting Lepontine dome.     

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

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

模拟退火剩余静校正的一步算法

当地震资料的剩余静校正量大、信噪比低时，利用传统的线性反演方法不能有效地拾取静校正量。本文介绍一种非线性反演方法（一步模拟退火法）来估算炮点和检波点的剩余静校正量。通过对合成地震记录的试算表明，该算法可获得近似于真实的剩余静校正量。     

层析成像中的波形反演

在波动方程有限差分数值模拟基础上，用波形记录残差平方并建立目标函数。对目标函数进行梯度优化，我们得出了介质速度的迭代公式。通过计算目标函数梯度我们可以对介质进行速度校正，从而获得反演解。经数值试验证明这是一种很有效的反演方法。     

井间电阻率层析成象的平滑度约束反演

井间电阻率层析成象是一种探测地下浅部精细结构的物探方法，它主要用于解决工程地质问题。本文采用２．５－Ｄ有限元法，针对井间电阻率层析成象中的ＡＭ观测系统进行了反演成象计算。文中首先计算了灵敏度矩阵，然后给出了一种电阻率层析成象反演算法—平滑度约束反演，理论模型的计算和实际资料的处理，都证明了该算法的有效性     

地球物理反问题中的随机反演理论

本文从最大后验概率密度观点出发，在数据噪音向量和待求模型向量为具有零均值的独立高斯随机过程的假设前提下，建立起了随机反演的非线性系统方程；给出了模型方差估计的函数表达式，并在文章最后，证明了反演解的稀疏性，即解释了随机反演的输出解的高分辨率特征。文章在最小二乘反演方法的基础上，发展并完善了随机反演方法的理论基础；揭示了随机反演方法与最小二乘反演方法之间的本质区别；阐述了随机反演方法的优越性，并指出了其广阔的应用前景。