时间域航空电磁各向异性大地三维自适应有限元正演研究   总被引：1，自引：0，他引：1       下载免费PDF全文

齐彦福  李貅  殷长春  戚志鹏  周建美  孙乃泉  刘云鹤  张博 《地球物理学报》2020,63(6):2434-2448

各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.  相似文献   

回线源瞬变电磁法有限体积三维任意各向异性正演及分析          下载免费PDF全文

刘亚军  胡祥云  彭荣华  Pritam Yogeshwar 《地球物理学报》2019,62(5):1954-1968

目前,瞬变电磁法(TEM)数据基本都是基于各向同性模型进行反演解释,这对于存在明显电性各向异性的勘探区域会产生较大的反演解释误差.为分析电各向异性对回线源瞬变电磁信号的影响方式与程度,本文通过求解离散化的全张量电导率时间域Helmholtz方程,实现了基于有限体积法的TEM任意各向异性的三维正演算法.该算法采用基于交错网格的拟态有限体积法(MFV)对时域Maxwell方程组进行空间域离散,并利用后退欧拉算法(Backward Euler Method)进行时间域离散.为提高时域电磁场的求解精度与效率,该算法将时间分段等步长算法与方程直接求解法相结合.通过对一维各向异性模型以及三维复杂各向同性模型进行测试,验证了本算法对于回线源瞬变电磁响应计算的正确性及有效性.最后,通过对几类典型电各向异性介质中大回线源瞬变电磁信号响应的分析,总结了不同电各向异性类型对TEM电磁信号的影响模式,结果表明,主轴各向异性情况下TEM信号主要受水平方向电导率的影响,倾斜各向异性对TEM信号的影响程度远大于水平各向异性,而通过水平各向异性信号能较清晰判断出各向异性主轴方向.  相似文献   

用于声波正演模拟的时空域高精度交错网格有限差分方法     

李斌  温明明  牟泽霖 《CT理论与应用研究》2017,26(3):259-266

地震正演模拟是逆时偏移和全波形反演中的核心问题之一,因为它们都需要高效、高精度地模拟波场正向和反向传播。为了提高数值模拟的精度,人们广泛采用高阶有限差分方法,但是大多数方法仅在空间上具有更高的精度,在时间上只有二阶精度。首先系统介绍时空域高精度交错网格有限差分方法的基本原理,然后利用模型验证方法的有效性,结果表明:时空域高精度交错网格有限差分方法拥有比常规交错网格有限差分方法更低的数值频散。   相似文献   

多频可控源电磁法三维有理函数Krylov子空间模型降阶正演算法研究          下载免费PDF全文

周建美  刘文韬  刘航  李貅  戚志鹏 《地球物理学报》2018,61(6):2525-2536

本文采用有理函数Krylov子空间模型降阶算法实现了同时求解多频可控源电磁法三维正演响应的快速计算.首先采用基于Yee氏交错网格的拟态有限体积法实现控制方程的空间离散,将任意频率的电场响应表示为关于频率参数的传递函数.采用有理函数Krylov子空间算法求解该传递函数.针对构建m维有理函数Krylov子空间需要求解m次（几十到上百）关于有理函数极点和离散控制方程系数矩阵的线性方程组的问题,本文提出采用单个重复极点的有理函数Krylov子空间模型降阶算法,结合直接法求解器PARDISO,采用Gram-Schmidt方法,只需要1次系数矩阵分解和m次矩阵回代即可实现有理函数Krylov子空间的构建,极大地减少了计算量.针对最优化有理函数极点选取问题,本文根据传递函数的有理函数Krylov子空间投影算法的误差分析理论,引入关于单个重复极点的收敛率函数,通过求解有理函数的最大收敛率直接给出最优化的单个重复极点公式.最终实现了不同发射频率的可控源电磁法三维正演响应的快速计算.分别计算了典型层状模型多发射频率的CSAMT和海洋CSEM的正演响应,通过与解析解的对比验证了本文算法在多发射频率正演的计算精度和计算效率;并通过一个三维海洋CSEM勘探设计最优化发射频率和接收区域选取的例子进一步说明本文算法的优点.  相似文献   

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL          下载免费PDF全文

TANG Xin-gong  ZHANG Rui-feng  WAN Wei  CHEN Qing-li 《地震地质》2017,39(4):712-720

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.  相似文献   

断层识别技术及其在MB油气田的应用   总被引：1，自引：2，他引：1       下载免费PDF全文

刘彦  孟小红  胡金民  吴江南 《地球物理学进展》2008,23(2):515-521

结合MB凸起油气勘查项目,总结断层解释方法:既要依据地震时间剖面特征来分析判断,更要了解盆地范围内应力场、变形机理以及工区范围内的构造样式.针对微小断层识别困难这一问题,采取多种技术手段,包括:(1)叠后目标处理提高地震资料对小断层的分辨能力;(2)正演模型分析指导小断层的解释;(3)相干数据体水平切片和近反射道几何特征分析识别微小断层;(4)三维数据体任意线剖面检验微小断层,识别了常规地震剖面难以确定的小断层,达到了断层精细解释的目的.  相似文献   

Linear dynamic modeling of a uni‐axial servo‐hydraulic shaking table system     

J. P. Conte  T. L. Trombetti 《地震工程与结构动力学》2000,29(9):1375-1404

This paper focuses on the development of a linear analytical model (even though servo‐hydraulic actuation systems are inherently non‐linear, especially for large amplitude simulations — near the performance capacity of the system — linearized models proved experimentally to be quite effective overall in capturing the salient features of shaking table dynamics) of a uni‐axial, servo‐hydraulic, stroke controlled shaking table system by using jointly structural dynamics and linear control theory. This model incorporates the proportional, integral, derivative, feed‐forward, and differential pressure gains of the control system. Furthermore, it accounts for the following physical characteristics of the system: time delay in the servovalve response, compressibility of the actuator fluid, oil leakage through the actuator seals and the dynamic properties of both the actuator reaction mass and test structure or payload. The proposed model, in the form of the total shaking table transfer function (i.e. between commanded and actual table motions), is developed to account for the specific characteristics of the Rice University shaking table. An in‐depth sensitivity study is then performed to determine the effects of the table control parameters, payload characteristics, and servovalve time delay upon the total shaking table transfer function. The sensitivity results reveal: (a) a potential strong dynamic interaction between the oil column in the actuator and the payload, and (b) the very important effect of the servovalve time delay upon the total shaking table transfer function. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

前散式和透射式能见度仪的综合观测试验及对比分析北大核心          下载免费PDF全文

王惠颖  祖繁  袁成松  包云轩  吴泓  褚进华 《气象科学》2022,42(2):225-234

通过在野外观测基地和能见度计量检测实验室开展的前散式和透射式能见度仪对比观测试验,获取了能见度和其他气象与环境要素的同步观测数据,分析了不同大气环境条件下能见度测量值差异及其产生原因。结果表明:在能见度环境模拟舱内,两种能见度仪的测量值具有较好的一致性,但在自然环境中,两种能见度仪的测量值差异随能见度上升而快速增大。在不同天气现象出现时,前向散射仪与大气透射仪的测量值比值(VIS_(FS)/VIS_(T))规律性变化明显,总体呈现随能见度上升而增大的趋势。无降水时,低湿情况下大气颗粒物吸收作用明显,透射能见度仪测量值远小于前散能见度仪测量值,两者均随PM_(2.5)质量浓度的增加而下降,且其差异也随之减小;高湿情况下大气颗粒物的散射作用在大气消光中占绝对优势,两种能见度仪测量值趋近。  相似文献