首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 156 毫秒
1.
自然电场对多种地下污染物敏感,开展相应的场源反演成像研究可帮助有效监测污染源活动.把目标区域内的多种自然电场场源进行统一处理,利用有限单元构建二维地电模型,实现自然电位正演的高精度计算.采用最小二乘正则化反演迭代,实现自然电场场源的二维反演.在设计二维自然电场场源反演算法的基础上,利用加入高斯噪声的合成模拟数据,对反演算法进行测试.测试结果表明:基于有限单元二维模型的最小二乘正则化反演算法能有效地实现对自然电场场源的反演成像,并能准确地确定自然电位异常源的位置和深度,且算法收敛快速稳定,具有一定的抗噪声能力.  相似文献   

2.
研究表明,同一区域地电暴记录与地磁场变化率在形态、幅度及周期成分上呈现出较强的相关性。采用磁暴期间电磁场观测数据,基于一维水平层状电性介质模型对地电暴波形进行拟合,反演地下介质电性结构。结果表明,地电暴拟合波形与观测记录一致性较好,利用反演确定的大地电性结构参数和其他磁暴事件,计算得到的理论感应地电场与观测结果相符合。  相似文献   

3.
为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.  相似文献   

4.
华北地区地电暴时GIC及涡旋电流响应分析   总被引:4,自引:0,他引:4       下载免费PDF全文
章鑫  杜学彬  刘君 《地球物理学报》2017,60(5):1800-1810
通过地电场台址近地表介质电阻率和地电场值计算大地电流,从大地电流场中分离出涡旋电流;根据平面波理论和水平导电层模型,使用地磁暴观测数据在频率域计算地磁感应电场(GIE),由GIE计算地磁暴感应电流(GIC).计算结果与实测值对比分析表明:GIE计算结果与实测地电暴具有很好的相位一致性;GIC涡旋中心相对地电暴涡旋中心存在向SE漂移约3°的现象;磁暴时地磁场Z分量的幅值分布图中极大值区域与涡旋中心重合,可能是GIC涡旋中心偏移的原因.另外,根据电磁感应原理提出的等效环电流模型,在一定程度上解释了涡旋大地电流的形成机制.本项工作可应用于地磁观测与地电观测的相互校正,同时有助于认识地电暴对大地电流分布的影响.  相似文献   

5.
CSAMT单分量数据解释方法   总被引:3,自引:1,他引:2       下载免费PDF全文
可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.  相似文献   

6.
瞬变电磁反演存在高度的非线性特征,常用的最小二乘等线性反演方法往往对初始模型高度依赖,并且极易陷入局部最优解.本文基于观测数据与模拟数据的L1范数建立目标函数,采用模拟退火非线性全局最优化方法实现瞬变电磁一维反演.初始模型完全随机产生,通过指数函数退温机制模拟系统能量最小实现迭代,通过接收概率函数评价当前模型,实现局部最优解的跳出,最终实现全局最优化求解.通过数值算例发现,无论给定的反演层数等于还是大于设计模型,都可以获得较好的反演效果,因而可以在反演初始就设计较多的层数,实现反演模型的自动拟合;同时,利用含噪声数据反演进一步验证算法的稳定性.最后,对实测数据进行了反演测试,结果与钻孔编录基本一致,表明提出的基于L1范数的模拟退火反演可用于实测数据处理.  相似文献   

7.
Kalman滤波在沉降监测数据处理中的应用   总被引:1,自引:0,他引:1  
将卡尔曼滤波应用于建筑物变形监测数据分析,给出了离散卡尔曼滤波模型的建立思路及相应的精度评定公式,结合西安市某建筑物沉降观测数据,建立了相应的离散卡尔曼滤波数学模型,通过Matlab编程对观测数据进行处理,成果图像显示滤波值曲线与原始观测数据曲线的变化趋势基本一致。该模型较好地模拟了建筑物沉降的变化规律,对于改善沉降监测数据处理的精度效果非常理想。  相似文献   

8.
基于非均一场源的球坐标大地电磁模拟方法   总被引:1,自引:0,他引:1       下载免费PDF全文
大地电磁测深理论与数据处理解释均假定平面电磁波垂直入射地下空间,但随着研究尺度的逐渐扩大,使得因地球弧度产生的影响难以忽略.此时,传统笛卡尔坐标体系及平面波场源不再适用于大尺度的大地电磁数据正反演解释.本文提出并实现了一种基于球坐标系的大地电磁交错网格有限差分三维正演算法,并对电场进行极向-环向分解,结合球谐函数和贝塞尔函数构建了可取代平面波的场源模型.首先利用经纬度信息构建三维地电模型,将场源设置于模型空间正上方,然后通过直接求解球坐标系下麦克斯韦方程来获得大地电磁响应.在此基础上,本文设计了球坐标下具有不同分辨率的多个三维地电模型,阐述了由球体模型到笛卡尔模型的转换方法,详细对比了两种坐标体系在计算效率、所求得的电场和视电阻率方面的差异.结果表明二者差异度主要与电性横向分布和地图投影方法有关,与周期并不存在明显的单调递增关系.  相似文献   

9.
在地电观测仪器技术比测规范的编写中,为了衡量比测仪器的稳定性和观测数据的一致性,会设置具体的技术指标值来进行约束。为了验证技术指标设置值的合理性,利用地电场仪、直流地电阻率仪对比观测所得数据,按照规范编写中的方法计算技术指标值,经分析认为,地电观测仪器技术比测规范中的技术指标值设置合理,为地电观测仪器技术比测规范技术的编写和地电观测仪器台站比测提供依据与支撑,从而为规范地电观测仪器和产出用于地震监测、预报的高质量地电观测数据提供保障。  相似文献   

10.
多通道瞬变电磁m序列全时正演模拟与反演   总被引:8,自引:7,他引:1       下载免费PDF全文
传统瞬变电磁方法主要用于金属矿勘查,无法满足油气资源高阻目标体的勘探需要.多通道瞬变电磁(MTEM)系统的出现解决了这一问题.该方法采用伪随机序列发射波形和拟地震观测方式,测量同线电场分量,记录全时发射电流及多道观测数据,实现对高阻薄层的高精度探测.鉴于国内对此方法的研究还处于理论探索阶段,尚未进行相应的仿真模拟和数据处理工作,本文针对m序列发射波形多通道瞬变电磁法的全时正演模拟和反演解释进行研究,为国内正在进行的MTEM仪器系统研发及数据解释提供理论指导.我们利用方波响应移位叠加和电流导数与阶跃响应褶积两种方法实现理论m序列和实际发射波形的全时正演模拟;再通过相关辨识技术,削弱噪声影响,计算脉冲响应;最后对积分得到的阶跃响应进行共中心点道集数据联合反演,获取地下电性分布信息.  相似文献   

11.
The assessment of hydraulic conductivity of heterogeneous aquifers is a difficult task using traditional hydrogeological methods (e.g., steady state or transient pumping tests) due to their low spatial resolution. Geophysical measurements performed at the ground surface and in boreholes provide additional information for increasing the resolution and accuracy of the inverted hydraulic conductivity field. We used a stochastic joint inversion of Direct Current (DC) resistivity and self-potential (SP) data plus in situ measurement of the salinity in a downstream well during a synthetic salt tracer experiment to reconstruct the hydraulic conductivity field between two wells. The pilot point parameterization was used to avoid over-parameterization of the inverse problem. Bounds on the model parameters were used to promote a consistent Markov chain Monte Carlo sampling of the model parameters. To evaluate the effectiveness of the joint inversion process, we compared eight cases in which the geophysical data are coupled or not to the in situ sampling of the salinity to map the hydraulic conductivity. We first tested the effectiveness of the inversion of each type of data alone (concentration sampling, self-potential, and DC resistivity), and then we combined the data two by two. We finally combined all the data together to show the value of each type of geophysical data in the joint inversion process because of their different sensitivity map. We also investigated a case in which the data were contaminated with noise and the variogram unknown and inverted stochastically. The results of the inversion revealed that incorporating the self-potential data improves the estimate of hydraulic conductivity field especially when the self-potential data were combined to the salt concentration measurement in the second well or to the time-lapse cross-well electrical resistivity data. Various tests were also performed to quantify the uncertainty in the inverted hydraulic conductivity field.  相似文献   

12.
Time-lapse monitoring is a powerful tool for observing dynamic changes in the subsurface. In particular it offers the potential for achieving inversion results with increased fidelity through the inclusion of complementary information from multiple time-steps. This inclusion of complementary information can reduce the need for spatial smoothing, without adding inversion artifacts to the resulting images. Commonly used time-lapse inversion methods include the ratio method, cascaded time-lapse inversion, difference inversion and differencing independent inversions. We introduce two additional methods in which both time-lapse data sets are inverted simultaneously. In the first, called temporally constrained time-lapse inversion, inversion of both datasets is done under a single optimization procedure and constraints are added to the regularization to ensure that the changes from one time to another are smooth. In the second method, called simultaneous time-lapse inversion, the inversions at time 1 and time 2 are performed simultaneously and constraints of smoothness and closeness to a reference model are applied to the difference image produced at each iteration, and subsequently, the constraints are updated at each iteration. Through both a numerical and a field example we compare the results of common time-lapse inversion methods as well as the introduced approaches. We found that of the commonly used time-lapse inversion methods the difference inversion method produced the best resolution of time-lapse changes and was the most robust in the presence of noise. However, we found that the alternative approach of simultaneous time-lapse inversion produced the best reconstruction of modeled EC changes in the numerical example and easily interpretable high resolution difference images in the field example. Moreover, there was less tailoring of regularization parameters with our simultaneous time-lapse approach, suggesting that it will lend itself well to an automated inversion code.  相似文献   

13.
抗差自适应Kalman滤波算法中,抗差等价权矩阵和自适应因子的计算,要求观测信息具有多余观测量且准确可靠,但在动态变形监测应用中,通常滤波观测值仅为三维坐标且存在较强噪声和粗差的影响。为此,先对该算法中的自适应因子和抗差等价权矩阵的计算进行研究和改进,然后计算了某高速公路边坡的GPS动态监测数据。结果表明,抗差自适应Kalman滤波能够有效地抵制动态变形监测中观测值异常的影响。  相似文献   

14.
Time-lapse seismic data is useful for identifying fluid movement and pressure and saturation changes in a petroleum reservoir and for monitoring of CO2 injection. The focus of this paper is estimation of time-lapse changes with uncertainty quantification using full-waveform inversion. The purpose of also estimating the uncertainty in the inverted parameters is to be able to use the inverted seismic data quantitatively for updating reservoir models with ensemble-based methods. We perform Bayesian inversion of seismic waveform data in the frequency domain by combining an iterated extended Kalman filter with an explicit representation of the sensitivity matrix in terms of Green functions (acoustic approximation). Using this method, we test different strategies for inversion of the time-lapse seismic data with uncertainty. We compare the results from a sequential strategy (making a prior from the monitor survey using the inverted baseline survey) with a double difference strategy (inverting the difference between the monitor and baseline data). We apply the methods to a subset of the Marmousi2 P-velocity model. Both strategies performed well and relatively good estimates of the monitor velocities and the time-lapse differences were obtained. For the estimated time-lapse differences, the double difference strategy gave the lowest errors.  相似文献   

15.
应用基于EMD的小波阈值去噪方法,去除地电场观测资料中轨道交通干扰,并将小波阈值去噪法和EMD去噪法的效果相比较,结果表明:该方法能够滤除地电场信号中的地铁干扰,同时保留原始信号中微小的突变,突出有用信息,提高地电场台站观测数据的使用率,有较好的去噪效果。基于EMD的小波阈值去噪方法可推广到其他地球物理观测资料的去噪分析,甚至地电场与地电阻率同场地观测中人工供电干扰信号的剔除。  相似文献   

16.
For the exploration of near-surface structures, seismic and geoelectric methods are often applied. Usually, these two types of method give, independently of each other, a sufficiently exact model of the geological structure. However, sometimes the inversion of the seismic or geoelectric data fails. These failures can be avoided by combining various methods in one joint inversion which feads to much better parameter estimations of the model than the independent inversions. A suitable seismic method for exploring near-surface structures is the use of dispersive surface waves: the dispersive characteristics of Rayleigh and Love surface waves depend strongly on the structural and petrophysical (seismic velocities) features of the near-surface Underground. Geoelectric exploration of the structure Underground may be carried out with the well-known methods of DC resistivity sounding, such as the Schlumberger, the radial-dipole and the two-electrode arrays. The joint inversion algorithm is tested by means of synthetic data. It is demonstrated that the geoelectric joint inversion of Schlumberger, radial-dipole and two-electrode sounding data yields more reliable results than the single inversion of a single set of these data. The same holds for the seismic joint inversion of Love and Rayleigh group slowness data. The best inversion result is achieved by performing a joint inversion of both geoelectric and surface-wave data. The effect of noise on the accuracy of the solution for both Gaussian and non-Gaussian (sparsely distributed large) errors is analysed. After a comparison between least-square (LSQ) and least absolute deviation (LAD) inversion results, the LAD joint inversion is found to be an accurate and robust method.  相似文献   

17.
为了提高AVO(amplitude versus offset)反演结果的精度和横向连续性,本文提出了一种新的AVO反演约束方法,该方法结合贝叶斯原理和卡尔曼滤波算法实现了对反演参数纵向和横向的同时约束.文章首先结合反演参数的纵向贝叶斯先验概率约束和反演参数的横向连续性假设建立了与卡尔曼滤波算法对应的AVO反演系统的数学模型,然后将该数学模型代入卡尔曼滤波算法框架,利用卡尔曼滤波算法实现了双向约束AVO反演.二维模型测试和实际数据测试结果表明,相对于单纯的纵向贝叶斯先验概率约束,双向约束能更准确地刻画参数的横向变化,得到更准确、横向连续性更好的反演结果.  相似文献   

18.
林旭  罗志才 《地球物理学报》2016,59(5):1608-1615
多速率Kalman滤波方法可用于低采样率的位移和高采样率的加速度数据融合,而未知的噪声协方差信息则显著制约着多速率Kalman滤波精度.本文通过将多速率Kalman滤波转换为传统的单速率Kalman滤波,建立了Kalman滤波增益的自协方差矢量与未知的加速度谱密度和观测噪声参数间的线性函数模型,并采用最小二乘估计方法对未知的噪声协方差参数进行估计,进而有效地提高了多速率Kalman滤波精度.数值仿真和震动台实验结果验证了本文方法的正确性和有效性.  相似文献   

19.
崇明与南京台震前地电场变化异常分析   总被引:12,自引:6,他引:12       下载免费PDF全文
2001年11月3日和12月25日江苏南通地区分别发生了ML3.8和ML4.1的地震群事件. 本文以上海崇明地电场台和南京地电场台ZD9A大地电场仪的观测资料为研究对象,在阐述了多极距地电场观测去噪音原理的基础上, 研究了这两个台站地电场多极距观测系统的资料在震前的异常变化特征. 结果表明: ① 多极距观测去噪音方法是一种很好的识别地电场噪音的方法;② 在震前20天、2天及42天离震中较近的崇明台NS和NE方向上分别接收到了地震电信号, 其持续时间最长的为9天,而较远的南京台出现的异常信号都是噪音信号而非地震电信号;③ 地电场多极距观测系统中公用电极的布极方式是不可取的.   相似文献   

20.
Introduction Geoelectric field is one of the proper physical fields of the Earth. It is controlled by the vari-ous current systems originating outside of the Earth and the internal electrical structure of the Earth as well as the underground current origin. It carries a great amount of information coming from the Earths interior. At present, a lot of observations and studies about the geoelectric field have been made in the world. Many applications and studies of the geoelectric field in the…  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号