线源三维井地电法异常边界增强及地形校正技术

熊治涛 唐新功 张连群 李丹丹 余俊虎

熊治涛,唐新功,张连群,李丹丹,余俊虎. 2022. 线源三维井地电法异常边界增强及地形校正技术. 地震学报,44(0):1−15 doi: 10.11939/jass.20220074
引用本文: 熊治涛,唐新功,张连群,李丹丹,余俊虎. 2022. 线源三维井地电法异常边界增强及地形校正技术. 地震学报,44(0):1−15 doi: 10.11939/jass.20220074
Xiong Z T,Tang X G,Zhang L Q,Li D D,Yu J H. 2022. Anomaly edge enhancement and topographic correction technology of linear source 3D borehole-to-surface electrical method. Acta Seismologica Sinica,44(0):1−15 doi: 10.11939/jass.20220074
Citation: Xiong Z T,Tang X G,Zhang L Q,Li D D,Yu J H. 2022. Anomaly edge enhancement and topographic correction technology of linear source 3D borehole-to-surface electrical method. Acta Seismologica Sinica44(0):1−15 doi: 10.11939/jass.20220074

线源三维井地电法异常边界增强及地形校正技术

doi: 10.11939/jass.20220074
基金项目: 国家自然科学基金(42274087,41874119,42174083)联合资助
详细信息
    作者简介:

    熊治涛,博士后,主要从事电磁法勘探研究,e-mail:1228409354@qq.com

    通讯作者:

    唐新功,博士,教授,博士生导师,主要从事电磁法勘探与地球动力学研究,e-mail:tangxg@yangtzeu.edu.cn

  • 中图分类号: P139

Anomaly edge enhancement and topographic correction technology of linear source 3D borehole-to-surface electrical method

  • 摘要: 基于非结构网格有限元方法开展了三维复杂地电模型的线源井地电法的高效正演模拟研究,探讨了通过求取电场响应导数来刻画目标体边界范围、采用差异场地形校正技术消除地形影响等措施对井地电法成像的效果与精度的影响。并通过与解析解的对比,验证了本文数值解算法的有效性。模型计算结果表明:积水巷道的空间位置与走向均会引起视电阻率的显著变化,视电阻率变化率的极值准确且清晰地指示了巷道边界的位置;电位的归一化总水平导数极大地提高了井地电法对目标体复杂边界位置的识别能力;地形对井地电场分布的影响也是严重的,其视电阻率响应与地形形状近似呈对称关系,利用差异场技术能有效地削弱地形对井地电法高精度成像的影响。

     

  • 图  1  均匀半空间模型的电极布设及计算结果

    (a) 电极布设方式;(b) 数值解与解析解对比及其相对误差;(c) 电流密度分布

    Figure  1.  Electrode layout and calculation results of uniform half-space model

    (a) Electrode layout;(b) Comparison of numerical solution and analytical solution and their relative error; (c) Current density distribution

    图  2  模型剖面(a)和平面(b)以及测点分布(c)示意图

    Figure  2.  Schematic diagram of profile (a) and plane (b) of the model and distribution of measuring points (c)

    图  3  地表视电阻率等值线图

    Figure  3.  Contour map of apparent resistivity on surface

    图  4  模型剖面(a)和平面(b)示意图

    Figure  4.  Schematic diagram of profile (a) and plane (b) of the model

    图  5  巷道模型地表视电阻率平面等值线图(白色方框代表巷道在地表的投影)

    (a) 模型1;(b) 模型2;(c) 模型3;(d) 模型4

    Figure  5.  Contour maps of apparent resistivity on surface of roadway model where the white box represents the projection of the roadway on the surface

    (a) Model 1;(b) Model 2;(c) Model 3;(d) Model 4

    图  6  巷道模型视电阻率一阶导数曲线图

    Figure  6.  Curves of the first derivative of apparent resistivity of roadway models

    图  7  巷道模型3的视电阻率二阶导数曲线图

    Figure  7.  Curve of second derivative of apparent resistivity of the roadway model 3

    图  8  垮塌巷道模型平面示意图

    Figure  8.  Plan view of collapsed roadway model

    图  9  平面等值线图

    (a) 电位;(b) 视电阻率;(c) 异常电位;(d) 电位总水平导数负对数;(e) 电位的归一化总水平导数

    Figure  9.  Plane contour maps

    (a) Electric potential;(b) Apparent resistivity;(c) Abnormal electric potential;(d) Negative logarithm of the total horizontal derivative of the electric potential;(e) Normalized total horizontal derivative of the electric potential

    图  10  地形模型剖面(a)和平面(b)示意图

    Figure  10.  Schematic diagram of profile (a) and plane (b) of the topographic model

    图  11  当地垒(a,c)和地垫(b,d)深为20 m和40 m时纯地形模型地表视电阻率等值线图

    Figure  11.  Surface apparent resistivity contour maps of pure topography model with horst (a,c) and graben (b,d) being 20 m and 40 m deep

    图  12  纯地形模型视电阻率曲线图

    Figure  12.  Apparent resistivity curves of pure topography model

    图  13  带地形模型剖面(a)和平面(b)示意图

    Figure  13.  Schematic diagram of profile (a) and plane (b) of the topography model

    图  14  地垒(a)和地堑(b)地形未校正的地表视电阻率等值线图

    Figure  14.  Contour maps of surface apparent resistivity for uncorrected horst (a) and graben (b)

    图  15  地形校正之后以及水平地表的视电阻率等值线图

    (a) 地垒地形校正后结果;(b) 地堑地形校正后结果;(c) 水平地表模型结果

    Figure  15.  Apparent resistivity contour maps after topography correction and without topography

    (a) Results after topography correction for horst;(b) Results after topography correction for graben; (c) Model results without topography

    图  16  模型剖面(a)和平面(b)示意图

    Figure  16.  Schematic diagram of model profile (a) and plane (b)

    图  17  地表等值线图

    (a,b) 地形校正前后的视电阻率;(c,d) 地形校正前后的视电阻率归一化总水平导数

    Figure  17.  Surface contour map

    (a,b)Topographically corrected front Apparent resistivity before and after topographical correction;(c,d) Normalized total horizontal derivative of apparent resistivity before and after terrain correction

    表  1  巷道规模及位置

    Table  1.   Roadway scale and location

    模型巷道规模(按xyz顺序)巷道左侧距源/m
    1200m×20m×10m100
    2200m×20m×10m0
    320m×200m×10m100
    4500m×20m×10m-200
    下载: 导出CSV
  • 戴前伟,陈德鹏,陈勇雄,侯智超. 2013. 电法勘探中异常响应特征的增强算法及其实现[J]. 煤田地质与勘探,41(3):75–78. doi: 10.3969/j.issn.1001-1986.2013.03.018
    Dai Q W,Chen D P,Chen Y X,Hou Z C. 2013. The enhanced algorithms and its implementation for the abnormal response characteristics in electrical exploration[J]. Coal Geology &Exploration,41(3):75–78 (in Chinese).
    李长伟,阮百尧,吕玉增,段长生,杨庭伟,王建历. 2010. 点源场井-地电位测量三维有限元模拟[J]. 地球物理学报,53(3):729–736.
    Li C W,Ruan B Y,Lü Y Z,Duan C S,Yang T W,Wang J L. 2010. Three-dimensional FEM modeling of point source borehole-ground electrical potential measurements[J]. Chinese Journal of Geophysics,53(3):729–736 (in Chinese).
    谭河清,沈金松,周超,董辉,房锡业,张福莱. 2004. 井-地电位成像技术及其在孤东八区剩余油分布研究中的应用[J]. 石油大学学报(自然科学版),28(2):31–37.
    Tan H Q, Shen J S, Zhou C, Dong H, Fang X Y, Zhang F L. 2004. Borehole-to-surface electrical imaging technique and its application to residual oil distribution analysis of the eighth section in Gudong Oilfield[J]. Journal of the University of Petroleum, China, 28(2): 31-37 (in Chinese).
    汤井田,张继锋,冯兵,林家勇,刘长生. 2007. 井地电阻率法歧离率确定高阻油气藏边界[J]. 地球物理学报,50(3):926–931. doi: 10.3321/j.issn:0001-5733.2007.03.035
    Tang J T,Zhang J F,Feng B,Lin J Y,Liu C S. 2007. Determination of borders for resistive oil and gas reservoirs by deviation rate using the hole-to-surface resistivity method[J]. Chinese Journal of Geophysics,50(3):926–931 (in Chinese).
    王智,潘和平. 2014. 三维井地电阻率法异常响应特征增强算法模拟研究[J]. 石油物探,53(4):491–500. doi: 10.3969/j.issn.1000-1441.2014.04.016
    Wang Z,Pan H P. 2014. Research on the enhanced algorithms of the abnormal response characteristics for 3D borehole-to-surface resistivity method[J]. Geophysical Prospecting for Petroleum,53(4):491–500 (in Chinese).
    王智,吴爱平,李刚. 2018. 起伏地表条件下的井中激电井地观测正演模拟研究[J]. 石油物探,57(6):927–935. doi: 10.3969/j.issn.1000-1441.2018.06.015
    Wang Z,Wu A P,Li G. 2018. Forward modeling of borehole-ground induced polarization method under undulating topography[J]. Geophysical Prospecting for Petroleum,57(6):927–935 (in Chinese).
    薛国强,闫述,陈卫营. 2016. 电磁测深数据地形影响的快速校正[J]. 地球物理学报,59(12):4408–4413. doi: 10.6038/cjg20161202
    Xue G Q,Yan S,Chen W Y. 2016. A fast topographic correction method for electromagnetic data[J]. Chinese Journal of Geophysics,59(12):4408–4413 (in Chinese).
    杨沁润,谭茂金,李桂山,张福莱,白泽. 2020. 大斜度井和水平井井地三维电阻率数值模拟和联合反演[J]. 地球物理学报,63(12):4540–4552. doi: 10.6038/cjg2020O0137
    Yang Q R,Tan M J,Li G S,Zhang F L,Bai Z. 2020. Numerical simulation and joint inversion of three-dimensional borehole-to-surface resistivity of high deviated or horizontal wells[J]. Chinese Journal of Geophysics,63(12):4540–4552 (in Chinese).
    张天伦,张伯林. 1995. 消除直流电阻率三极梯度法中各种干扰的实验与研究[J]. 石油地球物理勘探,30(1):100–110. doi: 10.13810/j.cnki.issn.1000-7210.1995.01.013
    Zhang T L,Zhang B L. 1995. Research on removing noises in DC resistivity trielectrode gradient measurement[J]. Oil Geophysical Prospecting,30(1):100–110 (in Chinese).
    Kong F N,Johnstad S E,Røsten T,Westerdahl H. 2008. A 2.5D finite-element-modeling difference method for marine CSEM modeling in stratified anisotropic media[J]. Geophysics,73(1):F9–F19. doi: 10.1190/1.2819691
    Ku C C,Hsieh M S,Lim S H. 1973. The topographic effect in electromagnetic fields[J]. Can J Earth Sci,10(5):645–656. doi: 10.1139/e73-065
    Li L L,Han L G,Huang D N. 2014. Normalized edge detection,and the horizontal extent and depth of geophysical anomalies[J]. Appl Geophys,11(2):149–157. doi: 10.1007/s11770-014-0436-2
    Spitzer K. 1995. A 3-D finite-difference algorithm for DC resistivity modelling using conjugate gradient methods[J]. Geophys J Int,123(3):903–914. doi: 10.1111/j.1365-246X.1995.tb06897.x
    Xiong Z T,Tang X G,Li D D,Zhang L Q. 2019. Linear source CSAMT response simulation in the 2D anisotropic formation with topography[J]. J Appl Geophys,171:103861. doi: 10.1016/j.jappgeo.2019.103861
    Xiong Z T,Tang X G,Qiu W Z,Zhao C Y,Zhang L Q. 2020. New algorithm for three-dimensional borehole-to-surface apparent resistivity imaging based on unstructured mesh finite-element method[J]. IEEE Access,8:184–194. doi: 10.1109/ACCESS.2019.2961806
    Yang Q R,Tan M J,Zhang F L,Bai Z. 2021. Wireline logs constraint borehole-to-surface resistivity inversion method and water injection monitoring analysis[J]. Pure Appl Geophys,178(3):939–957. doi: 10.1007/s00024-021-02674-6
  • 加载中
图(17) / 表(1)
计量
  • 文章访问数:  23
  • HTML全文浏览量:  19
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-19
  • 修回日期:  2022-08-05
  • 网络出版日期:  2022-11-07

目录

    /

    返回文章
    返回