SEISMOLOGY AND GEOLOGY ›› 2024, Vol. 46 ›› Issue (4): 972-992.DOI: 10.3969/j.issn.0253-4967.2024.04.012

• Research paper • Previous Articles    

2.5-D CSAMT FORWARD NUMERICAL SIMULATION BASED ON ANISOTROPIC AND POLARIZED MEDIA

YU Sheng-hong(), TANG Xin-gong*(), XIONG Zhi-tao   

  1. Key Laboratory Exploration Technologies for Oil and Gas Resource, Ministry of Education, Yangtze University, Wuhan 430100, China
  • Received:2023-06-20 Revised:2023-10-13 Online:2024-08-20 Published:2024-09-23
  • Contact: TANG Xin-gong

基于各向异性和极化介质的2.5维CSAMT正演数值模拟

余胜红(), 唐新功*(), 熊治涛   

  1. 长江大学, 油气资源与勘探技术教育部重点实验室, 武汉 430100
  • 通讯作者: 唐新功
  • 作者简介:

    余胜红, 男, 1997年生, 现为长江大学固体地球物理学专业在读硕士研究生, 研究方向为电磁法勘探, E-mail:

  • 基金资助:
    国家自然科学基金(42274087); 国家自然科学基金(42174083)

Abstract:

The Controlled Source Audio Magnetotelluric(CSAMT)method is a frequency-domain electromagnetic sounding technique that utilizes an artificial field source. It has been widely used in the exploration of various resources. Due to tectonic stress, temperature variations, material migration, and geological deposition within the Earth, the electrical parameters of underground media may vary with the direction of current, exhibiting electrical anisotropy. Additionally, rocks and ores exhibit complex electrochemical phenomena under the influence of an artificial field source, resulting in resistivity being a frequency-dependent complex number, known as the induced polarization(IP)effect. Electrical anisotropy and the IP effect are common phenomena in Earth's internal media, significantly impacting CSAMT observation data. Traditional CSAMT methods typically assume that the underground medium is electrically isotropic and lacks IP effects when performing forward and inversion numerical simulations. However, the actual underground medium often exhibits electrical anisotropy, the IP effect, or both. Adopting an isotropic model assumption during data processing and interpretation can therefore lead to serious errors in the results. Current research on CSAMT mainly considers the influence of either anisotropy or the IP effect individually. Since the anisotropy and polarization effect model involves multiple parameters, this complexity increases the difficulty of research. Consequently, there are relatively few studies on the combined influence of anisotropy and the IP effect on CSAMT responses. With increasing demands for exploration accuracy, it is imperative to study the response characteristics of the CSAMT electromagnetic field under the combined effects of anisotropy and the IP effect. In this paper, we employ a 2.5-dimensional model to derive the electromagnetic field partial differential equation for a 2.5-dimensional anisotropic medium based on Maxwell's equations with a source. The 2.5-dimensional CSAMT numerical simulation of anisotropy is achieved using the Galerkin finite element method. The real resistivity is replaced by complex resistivity using the Cole-Cole model. We discuss the effects of the Euler angle and principal axis resistivity in electrical anisotropy parameters, as well as three polarization parameters in the Cole-Cole model, on the CSAMT electromagnetic field response. Finally, we study the controllable source electromagnetic field response considering both anisotropy and IP effects.
The results show that electrical anisotropy significantly impacts the response of 2.5-D CSAMT, primarily depending on the Euler angle and principal axis resistivity values. Depending on these values, electrical anisotropy can either increase or decrease the CSAMT response value. Among the three polarization parameters, polarizability has the greatest influence on the calculation results, with the apparent resistivity value gradually decreasing as polarizability increases. The time constant and frequency correlation coefficient have relatively smaller influences on the calculation results, and their effects are consistent. The combined effect of anisotropy and the IP effect on the 2.5-D CSAMT response is more complex. Electrical anisotropy, influenced by the Euler angle and principal axis resistivity, can either increase or decrease the CSAMT response result, while the polarization effect generally decreases the CSAMT response result. Thus, when both factors act simultaneously, their effects can overlap or offset each other depending on the values of the Euler angle and principal axis resistivity. The findings of this paper have significant theoretical implications for improving the accuracy of CSAMT data processing and interpretation.

Key words: CSAMT, anisotropy, induced polarization, finite element, 2.5-D, forward

摘要:

传统的CSAMT在进行正、 反演数值模拟时通常假设地下介质为电各向同性且无IP效应的介质, 而实际的地下介质往往同时存在电各向异性和IP效应, 因此开展各向异性和IP效应共同作用下的CSAMT电磁场响应特征的研究势在必行。文中基于含源的Maxwell方程组, 推导了2.5维各向异性介质的电磁场偏微分方程, 采用伽辽金有限元法实现了各向异性的2.5维CSAMT数值模拟, 进而讨论了欧拉角和各向异性系数对2.5维CSAMT响应的影响。通过引入Cole-Cole模型分别研究了不同极化参数对CSAMT响应的影响, 最后讨论了各向异性和IP效应同时存在时的电磁场响应特征。研究结果表明, 各向异性和IP效应均会对CSAMT的响应结果产生较为严重的影响。当各向异性和IP效应同时存在时, 根据欧拉角和主轴电阻率的取值, 二者的效果将相互叠加或抵消。文中的研究结果对于提高CSAMT资料处理与解释精度具有重要的理论意义。

关键词: 可控源音频大地电磁法, 各向异性, 激发极化, 有限元, 2.5维, 正演