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.     

库车坳陷新生代构造应力场恢复

库车坳陷的构造变形及演化与南天山造山带的发育密切相关。库车坳陷新近纪以来的喜马拉雅晚期构造运动最为强烈,形成了天山山前大型冲断带,并造就了现今的构造格局。通过对库车坳陷北部构造样式的识别及各种应力场指示标志的测量、统计和构造解析,对野外获得的应力场指示标志划分了期次,认为喜马拉雅晚期应力场标志为近南北向挤压。结合库车坳陷区域构造要素,如地质体几何形态、边界条件、岩石力学参数等,运用弹性力学有限元数值模拟的方法获得了喜马拉雅晚期库车坳陷的区域应力场。模拟结果表明,库车坳陷喜马拉雅晚期最大主压应力方向为近南北向,与古应力场标志拟合较高,可以为库车坳陷裂缝预测和评价提供依据,对勘探开发有应用价值。     

水平地震作用下无锚固储罐应力与应变响应分析

在考虑地基与储罐相互作用的情况下,采用有限元法对储罐在水平地震荷载作用下的应力及应变反应进行了数值计算。对3×104m3和2×103m3罐壁应力及应变的分析结果表明:环向、轴向应力及应变的分布形式呈现出明显的下部大上部小的特点,在偏底部的位置出现应力和应变的峰值;储液罐在水平地震作用下“象足”变形是由纵向压应力达到屈曲临界应力导致的屈曲破坏,不是强度破坏,即破坏属于失稳破坏而非强度破坏。     

水平地震激励下储罐液体晃动与提离分析

在考虑地基与储罐相互作用的情况下,采用有限元法对储罐在水平地震荷载作用下的液面晃动及储罐提离反应进行了分析。结果表明:罐内液体的晃动是长周期运动,其周期受地震动影响。无论是小体积罐还是大体积罐,在一定的地震烈度下均可以发生提离,而且储罐发生提离的时刻大多数是在地震动的峰值过后的一段时间内。大体积的储罐的提离明显小于小体积的罐。底板提离区域为月牙形。     

电缆地层测试最小测试时间的确定

电缆地层测试器是重要的地层评价仪器,可以测量地层压力和渗透率、识别流体类型、确定油水界面.然而在低孔低渗地层中,由于泥饼的封闭性差而造成的地层增压现象,致使无法测量到地层真实压力.同时由于长时间的测试,容易将仪器卡在地层中而造成事故.本文提出一种确定电缆地层测试最小测试时间的方法.利用有限元方法计算泥饼和原状地层中的压力分布情况,求取泥浆柱静压力的影响半径,继而确定储层的最小测试时间,并且针对渗透性不同的储层给出了最小测试时间的范围.电缆地层测试最小测试时间的确定,可以确保电缆地层测试器在低孔低渗储层中测量到地层的真实信息,减小由于测量时间过长而引起工程事故发生率.     

异形柱抗震性能分析

异形柱框架结构是1种新型节能住宅结构,本文利用有限元软件,采用反应谱方法,从异形柱与矩形柱框架结构在地震作用下的变形特性和位移着手,研究2种结构的抗震性能,得出了异形柱的抗震性能好的结论。     

A review of numerical experiments on seismic wave scattering

This paper reviews applications of the finite-difference and finite-element methods to the study of seismic wave scattering in both simple and complex velocity models. These numerical simulations have improved our understanding of seismic scattering in portions of the earth where there is significant lateral heterogeneity, such as the crust. The methods propagate complete seismic wavefields through highly complex media and include multiply scattered waves and converted phases (e.g.,P toSV, SV toP, body wave to surface wave). The numerical methods have been especially useful in cases of moderate and strong scattering in complex media where multiple scattering becomes important. Progress has been made with numerical methods in understanding how near-surface, low-velocity basin structures scatter surface waves and vertically-incident body waves. The numerical methods have proven useful in evaluating scattering of surface waves and body waves from topography of both the free surface and interfaces buried at depth. Numerical studies have demonstrated the importance of conversions from body waves to surface waves (andvice versa) when lateral heterogeneities and topographic relief are present in the uppermost crust. Recently, several investigations have applied numerical methods to study seismic wave propagation in velocity models which vary randomly in space. This stochastic approach seeks to understand the effects of small-scale complexity in the earth which cannot be resolved deterministically. These experiments have quantified the relationships between the statistical properties of the random heterogeneity and the measurable properties of high-frequency (1 Hz) seismograms. These simulations have been applied to the study of many features observed in actual high-frequency seismic waves, including: the amplitude and time decay of seismic coda, the apparent attenuation from scattering, the dispersion of waveforms, and the travel time and waveform variations across arrays of receivers.