2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method          下载免费PDF全文

Yuan Yuan  Jianke Qiang  Jingtian Tang  Zhengyong Ren  Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.  相似文献   

上地壳纵横波速度结构相关反演成像方法   总被引：3，自引：0，他引：3       下载免费PDF全文

王晓  周小鹏  张新彦  白志明  滕吉文 《地球物理学报》2015,58(10):3553-3570

基于纵横波初至走时数据的层析成像方法越来越广泛地被应用于揭示不同构造域壳幔速度结构特征.我们从同一地质体的纵横波速度属性相关这一基本思想出发,提出一种相关反演成像的方法:纵横波速度反演交替进行,在迭代反演过程中每通过一次反演获得相应的纵波速度(或横波速度)结构后,更新相应的纵横波速度比模型以及相应的横波(或纵波)速度反演的初始模型,然后继续开展后续横波(或纵波)速度反演工作.在反演过程中依据纵横波速度的相关性信息和射线路径长度将走时残差以不同权重分配到射线路径经过的单元,依据网格节点周围平均的慢度扰动更新速度模型.正反演过程分别基于有限差分走时计算方法和反投影成像方法.两种典型模型试验表明,该技术应用于上地壳速度结构反演成像过程,可有效提高反演结果的可靠性,在很大程度上避免了常规单独反演纵波和横波速度过程容易带来的畸变和失真.该方法应用于重建青藏高原西部札达—泉水沟深地震测深(DSS)剖面下方的上地壳速度结构,揭示出与青藏高原西缘板块碰撞相关的上地壳速度结构特征.  相似文献   

长江中下游成矿带中段岩石圈电性结构研究   总被引：2，自引：2，他引：0       下载免费PDF全文

王显莹  汤井田  张林成  薛帅  潘伟  周聪  肖晓  任政勇  吕庆田 《地球物理学报》2015,58(12):4403-4414

长江中下游成矿带位于大别造山带、长江中下游凹陷、江南隆起带等大地构造单元结合部位,通过在研究区内布设两条首尾相接共计150km长的大地电磁剖面,获得了50km以浅岩石圈尺度的电性分布.长江中下游地区中段地下电性结构显示出在地下10km和30km处分别存在明显的圈层结构,以此认为现今横向稳定的"电莫霍"反映了研究区经历燕山期陆内构造-岩浆活动后已基本上完成壳幔重新平衡;而分隔大地构造单元的郯庐断裂带、长江断裂带以及江南断裂带在电性上具有特征的梯度显现,在印支造山期后的引张背景下,断裂带成为强伸展活动带与控制了燕山期大范围的陆内岩浆活动;高导地幔的局域性存在以及从北向南地幔导电性的变化反映了在经受深部动力学过程中处于不同大地构造部位的地幔所遭受的不同类型的改造以及地幔深部的构造极性.  相似文献   

浅部频率域电磁勘探方法综述   总被引：11，自引：4，他引：7       下载免费PDF全文

汤井田  任政勇  周聪  张林成  原源  肖晓 《地球物理学报》2015,58(8):2681-2705

适用于近地表(2000m以内)勘探的频率域电磁法主要有音频大地电磁法(audio-frequency magnetotellurics,AMT),无线电大地电磁法(radio-magnetotellurics,RMT),可控源音频大地电磁法(controlled source audiofrequency magnetotellurics,CSAMT),广域电磁法(Wide Field Electromagnetic Method,WFEM).本文拟从最新的数据采集技术、数据处理技术、正反演算法、实例等四个方面,论述适用于浅部勘探的AMT,RMT,CSAMT和WFEM方法的国内外最新进展,总结目前AMT,RMT,CSAMT和WFEM方法遇到的困难,对潜在的发展方向提出建议.综述表明:(1)张量测量、多站阵列、多站叠加可提高AMT、RMT和CSAMT数据的质量.利用近区数据WFEM法可获得良好的效果.国产与国外仪器在质量方面的差距正在逐步缩小.(2)数学形态滤波技术、Hilbert-Huang变换等可有效分离出有用的数据,局部畸变仍然是亟待解决的难题,需要更为深入的研究.(3)矢量有限元与非结构网格的出现大幅度提高了有限元处理复杂电磁问题模拟的精度与应用范围,成为目前电磁正演的首选工具.完全非线性反演算法仍然局限于1D、2D问题,共轭梯度法和高斯牛顿算法等为解决3D问题的发展趋势.地质约束的引入和多数据联合反演可以减小反演的非唯一性.各向异性的反演为目前反演研究的热点之一.(4)野外数据解释的正确性严重依赖于对地下结构先期的维性判别,在2D特性不明显、3D特性明显时,需要采用3D进行反演解释.  相似文献   

Time–Space Decoupled Explicit Method for Fast Numerical Simulation of Tsunami Propagation     

Anxin Guo  Shengchao Xiao  Hui Li 《Pure and Applied Geophysics》2015,172(2):569-587

  相似文献   

塔里木盆地南缘2000-2013年人工绿洲扩张特点     

任晓  穆桂金  徐立帅  林永崇  赵雪 《干旱区地理》2015,38(5):1022-1030

利用2000、2005、2011、2013年4期MODIS13Q1、Landsat TM/ETM+/OLI遥感影像,基于RS、GIS技术,分别提取土地利用/覆被类型信息,探讨塔里木盆地南缘人工绿洲扩张特点。结果表明:塔里木盆地南缘人工绿洲在2000-2013年不断扩张,各时段扩增来源中,天然绿洲比重不断减小,绿洲-沙漠过渡带比重逐步上升,荒漠比重快速增加。各流域人工绿洲扩张的阶段性与区域性差异明显,且扩增来源结构由于扩张的相对区位不同而异。人口增长、政策引导使人工绿洲不断扩张,并具有阶段性与区域差异性特点。人工绿洲扩张过程中,天然绿洲和绿洲-沙漠过渡带有明显退化现象,使其生态屏障功能下降,威胁到绿洲安全和经济社会的可持续发展。  相似文献   

TRMM降水数据的空间降尺度方法研究   总被引：1，自引：1，他引：0  

李净  张晓 《地理科学》2015,35(9):1164-1169

以天山中段为研究区,以降水与高程（DEM）、植被指数（NDVI）、坡向（Aspect）、经纬度与之间存在的相关关系为基础,构建了TRMM卫星3B43降水数据与NDVI、DEM和坡向等相关因子的回归模型,对TsHARP统计降尺度算法引入DEM和局部Moran指数进行改进,得到2001~2006年研究区的250 m高分辨率降水空间分布数据。最后利用研究区气象站点降水数据对降尺度结果进行验证,得出降尺度的结果和实测值的误差明显小于原始降水数据和实测值的误差,年均降水最大改善程度是70 mm,因此构建的降尺度方法是合理可行的,可用于山区降水数据的降尺度研究。  相似文献   

单、双向振动条件下饱和粉砂变形特性对比试验研究          下载免费PDF全文

刘潇  张学  赵俭斌  侯世伟  郝哲 《地震工程学报》2015,37(3):748-753

由于仪器水平的限制,关于土体的动力特性的研究目前多采用单向激振方式替代双向激振方式。而这两种激振方式在试验原理方面存在明显不同。很多研究表明,应力路径对土的响应有着不可忽视的影响,这两种方式相替代的合理性尚待证明。本文利用DSC2000多功能动三轴试验系统,选取尾矿砂分别在单向及双向振动条件下进行一系列饱和土的动力特性试验,对比两种激振方式作用下土体的滞回曲线以及骨干曲线的差异,分析土体在这两种振动方式下产生差异的原因,并应用双曲线模型对骨干曲线进行拟合,最后给出在进行土体变形特性试验时以单向激振方式替代双向激振方式的适用条件。  相似文献   

地震速报微博发布系统的研发   总被引：1，自引：0，他引：1  

肖健  侯建民 《中国地震》2015,31(2):456-460

微博已成为当今数据发布和传播的最佳自媒体平台,将地震速报信息及时发布到微博上有利于信息的传播,使民众及时获取地震消息,提高地震系统的社会服务能力。地震速报微博发布系统能实时读取EQIM地震台网速报信息并将共享服务系统中的地震数据发送到指定微博账号,同时发布具有中国地震台网中心LOGO标识的震中位置图;在读取EQIM和发送微博时充分考虑到震级、台网代码、震中中文地名等多种过滤条件。该系统已运行近3年时间,运行稳定、可靠,达到了预期的效果,特别是在朝鲜核爆、芦山地震、鲁甸地震等地震事件中及时向社会公众提供了地震速报信息,取得了较好的社会效益。  相似文献   

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应     

杨骁  周磊  张敏 《岩土力学》2015,36(7):2013-2020

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。  相似文献