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.     

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

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

Development of A Depth-Integrated Longshore Current Model with Unstructured Grids

A depth-integrated model for simulating wave-induced longshore current was developed with unstructured grids. Effects of surface roller and horizontal mixing under combined waves and currents were incorporated in the numerical model. Recommended values of model coefficients were also proposed based on sensitivity analysis. Field observations and three series of laboratory measurements including two cases conducted on the plane beach and one implemented on the ideal inlet were employed to examine the predictive capability of this model. For the field case and laboratory cases conducted on the plane beach, numerical results were compared favorably with the measured data. For the case with an ideal inlet, simulated circulation pattern is supposed to be reasonable although some deviations between numerical results and measured data still can be detected.     

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

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

异形柱抗震性能分析

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

无结构网格二维河口海岸水动力数值模式的建立及其应用

为完全拟合河口近海复杂岸线和工程结构以及有效局部加密,设计并建立了一个无结构三角形网格二维河口海岸水动力数值模式。空间离散主要基于有限体积法以保证守恒性,时间积分采用预估修正法以提高精度。水位在三角形网格中心通过连续方程求解;水平x方向和y方向的流速U和V均在网格边中点上通过动量方程求解。流速平流项的求解中采用了TVD格式。TVD流速平流通量为一个一阶迎风格式通量和一个二阶格式通量的组合,一阶格式通量和二阶格式通量根据流速的局部分布情况得出配比,最终组合得到TVD通量。TVD格式具有低耗散和无频散的优点,提高了模式的稳定性。应用实测资料验证建立的模式,结果显示水位、流速和流向的计算值与实测值均符合良好。     

非结构化网格下近岸波生流数值模拟

波浪破碎产生的近岸流是近岸海域关键的水动力因素之一。基于近岸波浪的椭圆型缓坡方程和二维近岸波生流方程,建立了非结构化网格下近岸波浪破碎形成的近岸流数值模型。数值模型中,在空间上采用有限体积法进行数值离散,在时间上采用欧拉向前格式数值离散。数值计算结果表明,该数值模型可以有效地模拟近岸波浪破碎产生的近岸流。