共查询到20条相似文献，搜索用时 140 毫秒
1.
本文针对Groom-Bailey分解法在消除地表电性不均匀体对MT数据的畸变效应时存在的问题引入遗传算法,对基于传统线性最优化方法的GB分解算法进行改进,提出基于遗传算法的大地电磁阻抗张量分解方法.通过对理论合成数据以及三维/二维模型正演数据的分解试验,并且对青藏高原北缘阿尔金断裂地区实际MT数据进行分解处理,证明了基于遗传算法的GB分解能够更加有效地校正三维/二维情况下近地表三维电性不均匀体所造成的畸变影响.最后,在已有算法基础上,研究了基于遗传算法的多频率GB分解算法和MT数据静校正算法,并通过实际MT数据的处理证明了这些算法的有效性. 相似文献
2.
地壳和上地幔普遍存在的电性各向异性是大地电磁测深方法在实践应用中的关键和难点问题,而由近地表局部地电异常引起的电流型畸变可进一步加剧各向异性反演的多解性和不确定性.为压制电流型畸变影响,本文实现了一种利用大地电磁相位张量和倾子数据获取地下介质主轴电性各向异性结构的三维反演算法.该算法正演问题的求解采用交错采样有限差分法,反演问题的求解采用非线性共轭梯度法.理论模型测试表明,相位张量和倾子的三维各向异性反演不仅恢复了真实的各向异性异常分布特征,且压制了为拟合电流型畸变出现的假异常.由于相位张量和倾子数据不含幅值信息,先验模型选取对反演模型电阻率值的恢复效果影响较大,但初始模型的改变基本不影响反演结果.为进一步验证该算法的有效性,将其应用于USArray项目在美国西北部采集的109个长周期大地电磁实测数据.反演结果表明,相对于阻抗数据,相位张量和倾子数据的各向异性反演结果对浅层结构的刻画更为准确,与地表地质特征更为吻合. 相似文献
3.
由于浅部、局部不均匀三维异常体的影响,大地电磁测深(Magnetotelluric,MT)观测数据往往需要进行畸变校正.本文针对相位张量不受局部电场畸变影响的特点,研究相位张量分析约束下的大地电磁测深(MT)阻抗张量GB(Groom-Bailey)分解方法.分析由相位张量得到的相关参数,结合初步电性结构维性判断及二维或近似二维电性结构条件下电性主轴方向等信息,设计较合理的GB阻抗张量分解初始模型,提高GB算法的稳定性及可靠性.同时,在传统GB分解求解模型的基础上,将各向异性参数引入目标函数,并采用共轭梯度法进行模型的最优化求解.理论数据计算结果表明,MT阻抗张量GB分解方法的计算结果,依赖于初始模型中电性主轴角、阻抗张量值等参数的选择,相位张量分析约束下的GB分解方法能有效改善算法的稳定性和计算结果的可靠性.此外,研究表明,在合理的初始模型选择下,较之传统的GB分解方法,尽管各向异性参数仍旧无法求出,但带各向异性参数的GB分解方法在计算效率和结果可靠性方面均具有一定的优势.选择实际大地电磁测深点数据进行方法验证,取得了理想的电性分析结果. 相似文献
4.
通过引入模糊均值聚类(FCM)模型约束函数对电阻率与速度进行约束,开展二维直流电阻率法与地震初至波走时成像法联合反演研究.在地下浅层结构勘探中,通常低电阻率的地质体具有低速特征,较高电阻率的地质体表现为较高的地震波速度.直流电阻率法因为低电阻率区域吸引电流而对其敏感,地震走时成像法因为射线集中在高波速区而对高速体敏感,因此,两者联合成像能够大幅度提高反演效果.合成数据反演表明,直流电阻率法和地震初至波走时联合反演对于两类地质体的分辨能力均有提升,能够优势互补.尤其是引入FCM模型约束进行联合反演,根据已知物性进行监督学习,进一步提高了反演质量,改善了成像模型的分辨率. 相似文献
5.
在地质情况复杂化以及地质勘探精细化的要求下,有必要开展张量CSAMT正反演研究.考虑到各向异性对正反演结果影响较大,本文开展基于主轴各向异性介质的张量CSAMT正演研究.根据含源的电场控制方程,引入主轴各向异性介质的张量电阻率,将电场总场分解为一次场和二次场,一次场由Key的算法得到,二次场利用有限差分法计算,即实现了主轴各向异性介质的张量CSAMT三维正演.通过与二维有限元算法结果比较来验证算法的正确性;通过设计几组主轴各向异性三维模型,分析卡尼亚电阻率和相位的响应特征.结果表明卡尼亚电阻率ρxy和ρyx分别对异常体X方向和Y方向的电阻率变化敏感,说明张量CSAMT方法对识别各向异性具有优势.本研究为开展任意各向异性介质的张量CSAMT正反演研究提供基础. 相似文献
6.
电阻率法和背景噪声法是通过获得地下介质的电阻率和速度参数的分布来探究地球内部物质分布的非均匀性特征,联合反演可以有效地发挥两种方法的优势,减小反演多解性.本文采用有限内存拟牛顿（L-BFGS）算法实现了电阻率法和背景噪声法的单方法三维反演,然后基于电阻率和速度模型结构耦合约束,采用交替迭代方式实现了电阻率法和背景噪声法的三维联合反演.通过设计规则体组合模型和不规则体组合模型进行理论模型合成数据三维反演试算,结果表明：联合反演可以获得同时满足数据拟合差和结构更为相似的速度-电阻率模型;联合反演可以压制背景噪声单方法反演出现的假异常,获得更精细的速度结构;联合反演获得的电阻率模型对倾斜异常体、高阻覆层下方异常体、圈闭的高/低阻体等边界信息有明显的提升,有效克服电阻率法单方法反演的局限,提高深部电阻率的分辨率. 相似文献
7.
直流电阻率法以其对含水构造响应敏感且适应性强的特点,是隧道施工期超前地质预报的常用方法.针对目前定点源三极观测方式易受到测线附近旁侧异常干扰与聚焦观测方式探测距离过短无法三维成像的问题,提出了一种基于多同性源阵列电阻率的隧道超前探测新方法,利用环形布置的多同性源供电压制隧道掌子面后方的异常体的干扰,同时增强掌子面前方远处异常体的敏感度,阵列电极测量可为含水构造的三维反演成像提供三维数据.首先,开展了含水构造超前探测的数值正演模拟,发现多同性源方法的视电阻率纵向微分曲线极小值以及与横轴的交点位置可对异常体位置进行定性判断,探讨了多同性源方法对后方及侧向异常干扰的响应特征及识别方法.其次,将电阻率光滑约束与轨迹光滑策略引入蚁群反演算法,提出了一种最小二乘方法与改进蚁群算法相结合的混合反演算法,反演数值算例表明混合反演算法能够实现含水构造的三维成像,较好地反演出不同距离的含水构造位置,并使其边界较为清晰.为了验证多同性源探测方法的探测效果开展了物理模型试验,多同性源阵列电阻率法能较好的反映含水构造的存在,与实际含水构造的位置较为相符. 相似文献
8.
Z轴倾子电磁法(ZTEM)是一种极具前景的新型天然源频率域航空电磁法,具有勘探深度较大,工作效率高的特点.本文基于大地电磁法(MT)三维数据空间OCCAM反演算法,实现了三维ZTEM倾子资料的数据空间反演算法.将该算法应用到理论合成算例中进行反演,反映出ZTEM倾子反演对模型横向边界的约束优于垂向的特性,并且与MT阻抗反演结果进行对比,验证了所实现的倾子资料三维数据空间反演算法的有效性,表明ZTEM倾子反演在约束电阻率模型横向位置能力上更具优越性.该反演算法可用于ZTEM所整理得到的倾子数据进行三维定量反演,实现大区域的地质构造三维成像,获得地下地质模型的电阻率结构信息. 相似文献
9.
局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。 相似文献
10.
时间域激发极化法(Time-domain induced polarization method,简称为TDIP)已有的反演算法采用的是分步反演的思路,即先由视电阻率资料反演电阻率,固定电阻率再由视极化率资料反演极化率,这样就存在极化率结果严重依赖于电阻率反演结果的问题.为了有效解决这一问题,本文实现了TDIP二维数据空间分步反演算法,提出了基于交叉梯度约束的TDIP二维同步反演策略,实现了交叉梯度约束的电阻率和极化率二维同步反演算法.分别用电阻率和极化率结构一致和不一致的二维模型合成数据进行了分步和同步反演试算,对不同模型试算结果进行了对比分析.结果表明:对于电阻率和极化率结构一致和不一致模型,同步反演结果比分步反演结果能更好地确定异常体的空间分布范围,反演得到的电阻率和极化率值更接近真值.理论模型算例表明本文提出的同步反演算法有效解决了分步反演的问题,优于分步反演算法,具有更好的实用性. 相似文献
11.
Since the true Earth is 3-D in nature, a three-dimensional (3-D) inversion has clear advantages over lower dimensional inversions. We utilized a 3-D magnetotelluric (MT) inversion code, WSINV3DMT, to obtain a realistic resistivity model using a long period MT data set collected in the Northwest Sabalan geothermal field in Ardabil, Iran. The apparent resistivity and phase curves, the magnetic induction vectors, the impedance polar diagrams and the rotational invariant of impedance tensor, indicate a complex 3-D conductivity structure. After setting up the model parameters and designing the appropriate block discretization, we performed the 3-D inversions for two sets of observed data; one set includes the full MT impedance tensor and another set contains only off-diagonal elements of the MT impedance. The final model was selected according to the relative magnitude of the data misfit and the model norm with respect to various Lagrangian multipliers. The results of this study illustrate the 3-D inversion of the off-diagonal elements of MT impedance tensor is precisely enough to explain the structures related to the geothermal source. The obtained results were compared with the results of available 2-D models and they are then interpreted using all of the geological and drilling data of the area. The main outcome of this study is the precise delineation of the geometry of geothermal source that is located at the center of the study area with a surface coverage of about 7 km2. 相似文献
12.
George R. Jiracek 《Surveys in Geophysics》1990,11(2-3):163-203
13.
14.
目前大地电磁(MT)测深资料反演主要基于各向同性介质,但随着MT实际应用的需要,各向异性研究已逐渐引起关注.我们采用广泛应用的广义逆法对一维MT水平层状各向异性介质模型反演进行了探索性研究,并实现了MT全张量响应(即所有的阻抗张量的视电阻率和相位)的一维各向异性反演.理论模型试验表明,无论理论观测值中是否含有噪声,这种方法都能够较好地恢复真实模型,验证了其正确性和有效性.将此方法用于MT实测资料时,能够同时拟合4对视电阻率和阻抗相位曲线,说明本方法可以用于实测资料的处理解释,具有一定的实用价值. 相似文献
15.
16.
针对三维／二维地质结构开发了大地电磁阻抗张量分解技术,利用该技术仅需对各个独立测点的大地电磁资料进行分析即可分离和消除由三维局部异常体引起的其附近测点的ＭＴ响应函数的形态扭曲畸变。 相似文献
17.
18.
A method for obtaining the galvanic distortion matrix is presented so that the regional impedance tensor (free of distortion) is recovered. The method is a constrained stochastic heuristic method, which consists in randomly exploring the space of the distortion parameters. Constraints are imposed on the shortest periods of the regional impedance tensor that, at these short periods, tends to be two dimensional (or one dimensional). Depending on the constraints used, two different methods to recover the regional impedance tensor in this 2D/3D case are presented. Method 1 needs to find the strike of the short periods and Method 2 applies to the measurement directions. Twist, shear and anisotropy parameters are obtained. Thus, the regional impedance tensor is recovered with the only exception being the vertical shift due to the gain, which is equal for all the components of the tensor. Examples with synthetic impedance tensors from 2D/3D models perturbed with galvanic distortion are presented to illustrate how the algorithm works. The presence of noise in data is considered and rules for proceeding are provided. The same examples perturbed by Gaussian noise together with experimental data illustrate the capabilities of the algorithm. 相似文献
19.
本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术. 相似文献
20.
APPLYING 3D INVERSION OF SINGLE-PROFILE MAGNETOTELLURIC DATA TO IDENTIFY THE SHADE AND YUNONGXI FAULTS 
下载免费PDF全文
下载免费PDF全文 JIANG Feng CHEN Xiao-bin DONG Ze-yi CUI Teng-fa LIU Zhong-yin WANG Pei-jie 《地震地质》2019,41(6):1444-1463
Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures.
In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults. 相似文献