Seismicity rate and wave-velocity variations as consequences of rainfall: the case of the catastrophic storm of September 2002 in the Nîmes Fault region (Gard, France)   总被引：1，自引：0，他引：1  

Alexis Rigo  Nicole Béthoux  Frédéric Masson  Jean-François Ritz 《Geophysical Journal International》2008,173(2):473-482

  相似文献   

Post-seismic rebound of a spherical Earth: new insights from the application of the Post–Widder inversion formula     

D. Melini  V. Cannelli  A. Piersanti  G. Spada 《Geophysical Journal International》2008,174(2):672-695

The post-seismic response of a viscoelastic Earth to a seismic dislocation can be computed analytically within the framework of normal-modes, based on the application of propagator methods. This technique, widely documented in the literature, suffers from several shortcomings; the main drawback is related to the numerical solution of the secular equation, whose degree increases linearly with the number of viscoelastic layers so that only coarse-layered models are practically solvable. Recently, a viable alternative to the standard normal-mode approach, based on the Post–Widder Laplace inversion formula, has been proposed in the realm of postglacial rebound models. The main advantage of this method is to bypass the explicit solution of the secular equation, while retaining the analytical structure of the propagator formalism. At the same time, the numerical computation is much simplified so that additional features such as linear non-Maxwell rheologies can be simply implemented. In this work, for the first time, we apply the Post–Widder Laplace inversion formula to a post-seismic rebound model. We test the method against the standard normal-mode solution and we perform various benchmarks aimed to tune the algorithm and to optimize computation performance while ensuring the stability of the solution. As an application, we address the issue of finding the minimum number of layers with distinct elastic properties needed to accurately describe the post-seismic relaxation of a realistic Earth model. Finally, we demonstrate the potentialities of our code by modelling the post-seismic relaxation after the 2004 Sumatra–Andaman earthquake comparing results based upon Maxwell and Burgers rheologies.  相似文献   

一种用于TEM高灵敏度感应式磁场传感器设计   总被引：2，自引：1，他引：1       下载免费PDF全文

刘凯  米晓利  朱万华  闫彬  刘雷松  方广有 《地球物理学报》2014,57(10):3485-3492

瞬变电磁法(TEM)作为地球物理非地震探勘主要电磁方法之一,广泛应用于油气、矿产等地下资源探测.本文分析TEM磁场传感器频率响应特点,采用欠阻尼匹配拓展传感器线性响应频率范围;分析磁传感器各类噪声源分布,提出不同频段上影响磁传感器性能的主要因素.最终研制一款高磁导率磁芯的TEM感应式磁场传感器,工作频段为10 mHz~1 kHz,等效输入噪声水平为4pT/√Hz@1Hz, 15fT/√Hz@200 Hz,传感器直径为50 mm,长度800 mm,重量3.2 kg,相比现有TEM接收磁场传感器,体积小、重量轻、工作频段低.实验测试与生产作业表明:最新研制的高灵敏度TEM磁场传感器更加适合大深度、3D阵列电磁勘探.  相似文献   

Coil response inversion for very early time modelling of helicopter‐borne time‐domain electromagnetic data and mapping of near‐surface geological layers     

Cyril Schamper  Esben Auken  Kurt Sørensen 《Geophysical Prospecting》2014,62(3):658-674

Very early times in the order of 2–3 μs from the end of the turn‐off ramp for time‐domain electromagnetic systems are crucial for obtaining a detailed resolution of the near‐surface geology in the depth interval 0–20 m. For transient electromagnetic systems working in the off time, an electric current is abruptly turned off in a large transmitter loop causing a secondary electromagnetic field to be generated by the eddy currents induced in the ground. Often, however, there is still a residual primary field generated by remaining slowly decaying currents in the transmitter loop. The decay disturbs or biases the earth response data at the very early times. These biased data must be culled, or some specific processing must be applied in order to compensate or remove the residual primary field. As the bias response can be attributed to decaying currents with its time constantly controlled by the geometry of the transmitter loop, we denote it the ‘Coil Response’. The modelling of a helicopter‐borne time‐domain system by an equivalent electronic circuit shows that the time decay of the coil response remains identical whatever the position of the receiver loop, which is confirmed by field measurements. The modelling also shows that the coil response has a theoretical zero location and positioning the receiver coil at the zero location eliminates the coil response completely. However, spatial variations of the coil response around the zero location are not insignificant and even a few cm deformation of the carrier frame will introduce a small coil response. Here we present an approach for subtracting the coil response from the data by measuring it at high altitudes and then including an extra shift factor into the inversion scheme. The scheme is successfully applied to data from the SkyTEM system and enables the use of very early time gates, as early as 2–3 μs from the end of the ramp, or 5–6 μs from the beginning of the ramp. Applied to a large‐scale airborne electromagnetic survey, the coil response compensation provides airborne electromagnetic methods with a hitherto unseen good resolution of shallow geological layers in the depth interval 0–20 m. This is proved by comparing results from the airborne electromagnetic survey to more than 100 km of Electrical Resistivity Tomography measured with 5 m electrode spacing.  相似文献   

Transient process and optimal design of receiver coil for small‐loop transient electromagnetics     

Cigong Yu  Zhihong Fu  Huaiqing Zhang  Heng‐Ming Tai  Xuegui Zhu 《Geophysical Prospecting》2014,62(2):377-384

Distributed parameters of the receiver coils greatly affect transient electromagnetic signals over short time periods, causing a delay in the signal's effective sampling time and the loss of shallow exploration information. This paper investigates the influence of transient process on apparent resistivity calculation and analyses the relations between the error of apparent resistivity and receiver coil design. We find that, under the same effective area, different radii of the receiver coils lead to different levels of impact on the estimation of the apparent resistivity. An optimization model is proposed to determine the optimal receiver coil size that gives rise to the smallest estimation error of the apparent resistivity. The relationship between the optimal radius and the effective areas is developed, which serves as a guideline for the optimal receiver coil design. The results may provide a useful means for improving the accuracy of the small loop transient electromagnetic instrumentation for shallow‐depth mapping.  相似文献   

电性源瞬变电磁地空逆合成孔径成像   总被引：2，自引：1，他引：1       下载免费PDF全文

李貅  张莹莹  卢绪山  姚伟华 《地球物理学报》2015,58(1):277-288

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.  相似文献   

基于等值反磁通原理的浅层瞬变电磁法   总被引：7，自引：0，他引：7       下载免费PDF全文

席振铢  龙霞  周胜  黄龙  宋刚  侯海涛  王亮 《地球物理学报》2016,59(9):3428-3435

基于等值反磁通原理的瞬变电磁法是一种新的探测地下纯二次场的方法.该方法采用上下平行共轴的两个相同线圈通以反向电流作为发射源,且在该双线圈源合成的一次场零磁通平面上,测量对地中心耦合的纯二次场.理论计算和物理实验论证了该方法能够有效消除接收线圈本身的感应电动势,从而获得地下纯二次场的响应.理论推导和数值计算证明了该方法采用的双线圈源比传统瞬变电磁法采用的单线圈源对地中心耦合场能量更集中,因而有利于减少旁侧影响、提高探测的横向分辨率.实测试验表明该方法是浅层探测的一种有效方法.  相似文献   

横观各向同性饱和地基中埋置荷载的非轴对称瞬态响应   总被引：2，自引：0，他引：2       下载免费PDF全文

王小岗 《地球物理学报》2009,52(8):2084-2092

基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础.  相似文献   

Study of transient poroviscoelastic soil motions by semi-analytical and numerical approaches     

A. Mesgouez  G. Lefeuve-Mesgouez 《Soil Dynamics and Earthquake Engineering》2009

In this paper, the authors compare results obtained by semi-analytical and numerical approaches for the dynamic response of a poroviscoelastic soil under transient loads. The behaviour of the medium is governed by complete Biot formalism. The semi-analytical approach is based on Helmholtz decompositions and Fourier transforms, and yields exact solid and fluid displacements in the transformed domain. The numerical approach uses a C++$\mathrm{C}++$ object oriented programming finite element–finite difference code. Both methods give concurring results. Moreover, influence of viscous coupling on the response of the ground and visualization of the compressional wave of the second kind are discussed.  相似文献   

An efficient space-fractional dispersion approximation for stream solute transport modeling     

Chaopeng Shen  Mantha S. Phanikumar 《Advances in water resources》2009

The success of transient storage (TS) modeling for natural streams depends, in part, on the ability to describe the dispersion process accurately. Evidence based on stream tracer data shows that solute transport processes often do not follow the classical second-order dispersion model (e.g., early breakthrough and faster than Fickian travel times were observed). While models based on space-fractional dispersion are a promising alternative, different definitions of fractional derivatives exist in the literature. Unlike integer-order derivatives, fractional derivatives represent convolutions of concentration with long-range spatial correlation and numerical approximations can produce dense matrices. Therefore issues of both accuracy and computational efficiency need to be examined to successfully identify model parameters for natural streams. In this paper, we first compare the performance of several numerical approaches for solving the space-fractional dispersion equation. We examine three different numerical approaches to approximate the space-fractional derivatives including: (a) a fully-implicit scheme based on the shifted Grünwald–Letnikov (GL) approximation (b) a three-point implicit representation based on the GL formula and (c) a three-point implicit scheme based on mass conservation and the Caputo definition of the fractional derivative. We then use an operator-splitting technique to evaluate a TS model based on space-fractional dispersion (the FSTS model) and test the model against analytical solutions and stream tracer data. A sequence acceleration method (Richardson extrapolation) significantly improves the performance of all schemes examined. Results indicate that the fully-implicit GL method with Richardson extrapolation produces the most accurate solutions while the three-point implicit GL scheme has a stringent time-step restriction to produce acceptable solutions. The three-point implicit scheme based on the Caputo derivative produces accurate solutions in a fraction of the time taken by the fully-implicit GL method and represents the best trade-off between accuracy and computational efficiency for practical applications. The scheme is suitable for parameter estimation and is used to successfully describe tracer data in a natural stream.  相似文献