物探方法在强风化区填图中的应用探索

陈松; 陈长敬; 吴俊; 刘恋; 卜建军

物探方法在强风化区填图中的应用探索

陈松^1,,
陈长敬¹,
吴俊¹,
刘恋²,
卜建军^1, ,

1.
中国地质调查局武汉地质调查中心, 湖北武汉 430205
2.
安徽省勘查技术院, 安徽合肥 230031

基金项目:

中国地质调查局地质调查项目"特殊地质地貌区填图试点" DD20160060

"广东1：5万筋竹圩、连滩镇、泗纶圩、罗定县幅强烈风化区填图试点" 12120114005901

"珠三角阳江-珠海地区海岸带1：5万填图试点" DD20160064

详细信息

作者简介:
陈松(1985-), 男, 安徽淮南人, 工程师, 硕士, 主要从事地球物理勘察与研究工作。E-mail:anhuisongchen@163.com

通讯作者:
卜建军(1972-), 男, 高级工程师, 主要从事区域地质和地层古生物方面的研究工作。E-mail:jianjunbu@cug.edu.cn

中图分类号: P631
计量
- 文章访问数: 223
- HTML全文浏览量: 162
- PDF下载量: 36
- 被引次数: 0
出版历程
- 收稿日期: 2016-09-25
- 刊出日期: 2017-04-01

APPLICATION AND EXPLORATION OF GEOPHYSICAL METHODS IN GEOLOGICAL MAPPING IN STRONGLY WEATHERED AREA

1.
Wuhan Center of Geological Survey, CGS, Wuhan 430205, Hubei, China
2.
Geological Exploration Technologies institute of Anhui Province, Hefei 230031, Anhui, China

摘要

摘要: 华南强风化区覆盖较厚，露头少，植被发育，通行困难，利用传统的地质填图方法很难有效查明风化覆盖层及其下伏地层、岩石、构造、矿产等地质特征。笔者利用浅层地震、高密度电法、EH4电磁成像等综合物探方法，对风化层厚度、基岩的岩性和产状、隐伏断裂等进行探测试验。通过对3种方法所得数据的分析和野外验证结果，了解了3种物探方法的有效性和适用性，初步总结了在强风化区填图运用多种物探手段填图的方法。
- 浅层地震 /
- 高密度电法 /
- EH4电磁成像 /
- 强风化区 /
- 填图
Abstract: Due to the unique geological feature of Southern China area, such as thick overburden layerfew outcrops, thick-growing vegetation and difficult access, it is difficult to effectively identify the geological characteristics of weathered layers and its underlying strata, rocks, structures and minerals by means of traditional geological mapping methods. Through comprehensive geophysical methods including shallow earthquake, high density resistivity and EH4 electromagnetic imaging, the geological characteristics like weathering layer thickness, bedrock lithology and occurrence, and buried fault are tested. By analyzing the data obtained from this 3 methods and the results of field investigation, the validity and applicability of this 3 geophysical methods are studied, and a series of geophysical methods for geological mapping in strongly weathered area are preliminarily summarized.
- shallow earthquake /
- high density resistivity /
- EH4 electromagnetic imaging /
- strongly weathered area /
- mapping

HTML全文

0. 引言

在煤炭开采过程中，随着工作面推移对覆岩体自然应力平衡状态的破坏，采空区上覆岩层的次生应力场重新分布，岩层内部的应力也重新分布以达到新的平衡。在此过程中上覆岩层要发生断裂、冒落、弯曲和下沉，这种运动变化由下向上逐步发展，出现不同步的下沉，上覆岩层由于岩性、距采空区距离、厚度等方面的不同，在垂直方向出现了不同的位移差，进而形成了离层空间。离层带中经常有大量的积水，这些离层水对于下伏的岩石形成巨大压力，当压力积累到一定程度时有可能破坏断裂带，形成离层水害，是巨大的煤矿安全隐患。因此离层的量化研究对煤矿安全生产和采区人民生活的影响具有积极的意义。

科研工作者对于覆岩离层问题进行了大量的研究。从最初的开采沉陷方式^[1]，到岩层运动的物理及力学模型，岩层运动过程，离层发育规律，覆岩的离层条件，以及离层所带来的危害都做了很多探讨。杨伦等^[2]提出的岩层二次压缩理论，将地表沉陷与岩层的物理力学性质联系起来。宋杨等^[3]曾利用有限元方法研究支承压力显现的过程，揭示其变化规律与上覆岩层运动间的关系。苏仲杰^[4]等将离层岩层假设为岩板建立了岩板挠度的力学模型。梁运培^[5]提出了覆岩层移动的组合岩梁理论。张健全^[6]分析了覆岩离层可能出现的位置，并提出了计算最大离层的方法。然而这些研究大多停留于定性阶段，在定量研究方面却刚刚起步，对于离层进行三维模型分析的研究尚属空白，其主要困难在于覆岩薄层众多，离层量小可至毫米级大可到米级，而三维空间体积巨大，现有的三维有限元接触算法无法解决。

针对上述问题，笔者进行了两个方面的研究，其一为层状界面弱接触类三维力学有限元算法，其二为应用该算法对三维力学有限元模拟的计算机程序和网格自动生成进行研究。在研究基础上进行的大规模并行计算结果表明，上述算法可行，计算成果对煤炭安全生产具有指导性意义。

1. 模型建立与接触算法

1.1 层内模型

假定岩层内部介质变形服从弹性力学基本规律，其力学变形服从弹性力学方程。

1.1.1 平衡方程

描述弹性体内部的力学平衡性质，根据牛顿第二定律可得出：

$\left\{ \begin{array}{l} \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{xx}}}}{{\partial x}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{xy}}}}{{\partial y}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{xz}}}}{{\partial z}} + {f_x} = 0\\ \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{yx}}}}{{\partial x}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{yy}}}}{{\partial y}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{yz}}}}{{\partial z}} + {f_y} = 0\\ \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{zx}}}}{{\partial x}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{zy}}}}{{\partial y}} + \frac{{\partial {\mathit{\boldsymbol{\sigma }}_{zz}}}}{{\partial z}} + {f_z} = 0 \end{array} \right.$

(1)

式中，σ_xx，σ_yy，σ_zz，σ_yz，σ_xz，σ_xy为应力矢量，Pa；f_x，f_y，f_z为内力密度（N/m³），这里为重力；{f_x f_y f_z}={0 0 -ρg}，ρ为密度（kg/m³），g为重力加速度（N/kg）。

1.1.2 几何方程

描述应变和变形（位移）的关系。介质变形的幅度远小于空间尺度可以忽略，因此本文采用了小变形的几何方程，进行准确的模拟。几何方程如下：

$\left\{ \begin{array}{l} {\varepsilon _{xx}} = \frac{{\partial u}}{{\partial x}}\\ {\varepsilon _{yy}} = \frac{{\partial u}}{{\partial y}}\\ {\varepsilon _{zz}} = \frac{{\partial w}}{{\partial z}}\\ {\varepsilon _{yz}} = \frac{{\partial w}}{{\partial x}} + \frac{{\partial v}}{{\partial z}}\\ {\varepsilon _{xz}} = \frac{{\partial w}}{{\partial x}} + \frac{{\partial u}}{{\partial z}}\\ {\varepsilon _{xy}} = \frac{{\partial u}}{{\partial y}} + \frac{{\partial v}}{{\partial x}} \end{array} \right.$

(2)

式中，ε_xx，ε_yy，ε_zz，ε_yz，ε_xz，ε_xy为应变；u，v，w分别为x, y, z方向的位移，m。

1.1.3 本构方程

描述物体应力与应变（ε_xx，ε_yy，ε_zz，ε_yz，ε_xz，ε_xy）的关系，这里假定地下介质为线弹性体，E为杨氏模量，ν为泊松比，物性参数的选择直接决定了变形的分布情况，相关物性参数由钻井数据，以及查阅资料得出。本构方程如下：

$\left( \begin{array}{l} {\mathit{\sigma }_{\mathit{xx}}}\\ {\mathit{\sigma }_{\mathit{yy}}}\\ {\mathit{\sigma }_{\mathit{zz}}}\\ {\mathit{\sigma }_{\mathit{yz}}}\\ {\mathit{\sigma }_{\mathit{xz}}}\\ {\mathit{\sigma }_{\mathit{xy}}} \end{array} \right) = \frac{E}{{\left( {1 + v} \right)\left( {1 - 2v} \right)}}\left( {\begin{array}{*{20}{c}} {1 - v}&v&v&{}&{}&{}\\ v&{1 - v}&v&{}&{}&{}\\ v&v&{1 - v}&{}&{}&{}\\ {}&{}&{}&{0.5 - v}&{}&{}\\ {}&{}&{}&{}&{0.5 - v}&{}\\ {}&{}&{}&{}&{}&{0.5 - v} \end{array}} \right)\left( \begin{array}{l} {\varepsilon _{xx}}\\ {\varepsilon _{yy}}\\ {\varepsilon _{zz}}\\ {\varepsilon _{yz}}\\ {\varepsilon _{xz}}\\ {\varepsilon _{xy}} \end{array} \right)$

(3)

1.1.4 边界条件

第一类边界条件适用于给定边界面位移的情况：

$\begin{array}{l} u = {u_0}\\ v = {v_0}\\ w = {w_0} \end{array}$

(4)

其中，u₀，v₀，w₀分别为初始位移在x, y, z三个方向上的分量，m。

第二类边界条件适用于给定边界面力密度的情形：

$\left\{ \begin{array}{l} {T_x} = {\sigma _{xx}}{n_x} + {\sigma _{xy}}{n_y} + {\sigma _{xz}}{n_z} = {g_x}\\ {T_y} = {\sigma _{yx}}{n_x} + {\sigma _{yy}}{n_y} + {\sigma _{yz}}{n_z} = {g_y}\\ {T_z} = {\sigma _{zx}}{n_x} + {\sigma _{zy}}{n_y} + {\sigma _{zz}}{n_z} = {g_z} \end{array} \right.$

(5)

其中，g_x，g_y，g_z为表面面力密度（Pa）在x, y, z方向的分量；n_x, n_y, n_z为表面外法向的x, y, z方向分量。

对于如上弹性问题，其变分形式为：

$\begin{array}{l} \int\limits_v {{\sigma _{xx}}\delta {\varepsilon _{xx}} + } {\sigma _{yy}}\delta {\varepsilon _{yy}} + {\sigma _{zz}}\delta {\varepsilon _{zz}} + {\sigma _{yz}}\delta {\varepsilon _{yz}} + {\sigma _{xz}}\delta {\varepsilon _{xz}} + {\sigma _{xy}}\delta {\varepsilon _{xy}}\\ = \int\limits_V {{f_x}\delta u + {f_y}\delta v + } {f_z}\delta w + \int\limits_F {{T_x}\delta u + {T_y}\delta v + } {T_z}\delta w \end{array}$

(6)

1.2 层间模型

假定介质为n层，则第i层和第i+1层之间需满足接触边界条件。以w_i，1(x，y)表示第i层下边界垂直方向位移，w_{i, 2}(x，y)表示第i层上边界垂直方向位移，则接触边界条件表示为：

$0 ＞ {w_{i,1}}\left( {x,\mathit{y}} \right) ＞ {w_{i + 1,2}}\left( {x,\mathit{y}} \right)$

(7)

若第n层下边界为采空区，则有：

${W_{n,2}}\left( {x,\mathit{y}} \right) ＞ - D$

(8)

其中，D为采空高度，m。

1.3 接触问题有限元计算方法

首先假定相邻岩层之间不会发生强制粘结，各层可以发生自由形变，因此各层可进行独立的有限元变形计算。

根据各个岩层发生的形变进行比较分析，在边界上满足接触面位移条件，即0＞w_{i, 1}(x，y)＞w_{i+1, 2}(x，y)。若该式成立，则两个岩层可能出现离层，继续视为两个不同表面，表面为自由表面边界条件，亦即表面力为0。若相邻两层接触面位移满足条件w_{i, 1}(x，y)≤w_{i+1, 2}(x，y)＜0，表明i层下表面发生的形变大于i+1层上表面，则粘合分裂的两个点，其坐标均为(x，y)，在下一次迭代过程中两个点则视为一个点。考虑最底层的下边界垂直方向位移w_n，2(x，y)，如果满足w_n，2(x，y)≤-D，则将该点的垂直方向位移w设为第一类边界条件，否则考虑为自由表面边界条件，如此反复迭代直至收敛，则接触边界条件得到正确处理。

2. 三维离层模拟算例

本文以济宁二号煤矿西部采区（十一采区）存在的覆岩离层问题为例进行定量化模拟研究。利用地震成像理论确定离层水平空间范围，确定计算区域和适当的边界条件（见图 1），利用有限元预处理软件建立二维模型（见图 2）。

图 1 矿区研究范围

Figure 1. Studied mining area

下载: 全尺寸图片幻灯片

图 2 二维网格示意图

Figure 2. A sketch map of 2D grid

下载: 全尺寸图片幻灯片

利用生成的二维网格沿第三方向延展得到三维网格，同时记录两个相邻地层的网格对应关系，以便计算时处理接触边界。本文采用网格粘合技术（见图 3）处理接触边界条件，首先将相互连接的两个地层的同一个点分裂为2个不同的点，进行有限元计算，然后再根据计算结果接触情况进行网格粘合，图 4为研究区域的三维网格示意图。

图 3 网格粘合技术示意图

Figure 3. A sketch map of Grid-bonding

下载: 全尺寸图片幻灯片

图 4 三维网格示意图

Figure 4. A sketch map of 3D grid

下载: 全尺寸图片幻灯片

根据采区部分钻孔的测井数据（杨氏模量、密度和泊松比等），建立岩层空间分布信息以及不同岩层的岩性，再查阅相关资料得到不同岩层的力学模拟参数（见表 1、图 5）。

表 1 采区覆岩层空间分布及力学参数

Table 1. Spatial distribution and mechanical parameters of each layer

下载: 导出CSV

| 显示表格

图 5 覆岩层空间分布及力学参数示意图

Figure 5. A sketch map of spatial distribution of each layer

下载: 全尺寸图片幻灯片

本开采工作面采掘3_下煤层，采高5 m左右。利用上述三维网格以及物理参数，在大规模计算集群上进行了数值模拟。图 6列出了多次迭代收敛过程。

图 6 迭代结果示意图

Figure 6. Sketch of iterative results

下载: 全尺寸图片幻灯片

图 6a—图 6e代表每次迭代过程，从收敛过程可以看出，尽管第1次迭代时计算结果有上层穿插下层的情况，但是随着迭代步数增加，在第5次迭代完成后，收敛到最终的变形情况，显示出3个不同离层位置（见图 6e），分别位于累积深度-657.5 m，-707.55 m，-766.7 m处，相应的最大离层高度分别为0.2144 m、0.5261 m和0.1809 m（见表 2）。

表 2 离层数据

Table 2. tatistics of separation layers

下载: 导出CSV

| 显示表格

应该注意，模拟过程中应用了弹性力学方程，没有能力模拟采煤过程中冒落带形成的各种物理过程和获得其几何参数，因此离层3参数具有较大的不确定性，它可能已经位于冒落带和裂隙带范围（见图 6e）。

计算结果进行可视化，得到反映离层形态和体积的直观图，直观显示了离层形状特性以及水平位置分布。本文给出了模拟出的离层1和离层2的变形俯视图和变形三维图（见图 7、图 8），从图中可以看到离层位于采掘面后方，先变大，后闭合的过程（为了绘图的需要将x，y坐标归一化到[0, 1]）。

图 7 离层1变形图

Figure 7. Deformation of the 1^st layer

下载: 全尺寸图片幻灯片

图 8 离层2变形图

Figure 8. Deformation of the 2^nd layer

下载: 全尺寸图片幻灯片

为了显示离层，将离层2高度放大20倍，可以清楚地看到其几何特征（见图 9）。分析表明，该形态符合离层结构力学特征。

图 9 离层2剖面图（高度方向扩大20倍）

Figure 9. Profile view of the 2^nd layer

下载: 全尺寸图片幻灯片

根据唐山矿覆岩离层井下视频观测，覆岩离层的最大高度一般为几厘米到几十厘米^[7]。济宁二号矿与唐山矿均位于华北地层区，开采层位均为P—C煤层，采掘深度唐山矿500 m左右，济宁二号矿700 m左右，地质和几何空间尺度等条件相似，因此上述模拟得到的覆岩离层几何参数具有可信度，在煤矿安全生产方面也具有重要参考价值。

3. 结论

在覆岩离层的接触问题中，垂直方向的位移较大，水平方向位移相对较小，利用这一特点，首先考虑相邻层只在重力作用下的变形，然后判断垂直位移的大小，如果满足接触条件，则采用网格粘合技术将其粘合视为一个点，很好地解决了这一类接触问题。通过二维网格的三维延展，实现了大规模接触问题的网格生成问题。首次使用真正三维模型进行模拟，得到了离层体积量、最大离层高度等参数。

图 1 六次覆盖观测系统示意图

Figure 1. The sketch of six covering observation

下载: 全尺寸图片幻灯片

图 2 温纳装置测量原理示意图

Figure 2. The sketch of measuring principle of wenner device

下载: 全尺寸图片幻灯片

图 3 1线地震反射波法叠加剖面图

Figure 3. The stacked profile of seismic reflection wave method in line 1

下载: 全尺寸图片幻灯片

图 4 1线高密度电法电阻率反演断面图

Figure 4. Inversion section of high density resistivity in line 1

下载: 全尺寸图片幻灯片

图 5 1线EH4反演断面图

Figure 5. Inversion section of EH4 in line 1

下载: 全尺寸图片幻灯片

图 6 2线地震反射波法叠加剖面图

Figure 6. The stacked profile of seismic reflection wave method in line 2

下载: 全尺寸图片幻灯片

图 7 2线高密度电法电阻率反演断面图

Figure 7. Inversion section of high density resistivity in line 2

下载: 全尺寸图片幻灯片

图 8 2线EH4反演断面图

Figure 8. Inversion section of EH4 in line 2

下载: 全尺寸图片幻灯片

参考文献(20)

[1]	刘菁华, 王祝文.地面综合物探方法在浅覆盖区地质填图单元的划分研究[J].中国地质, 2005, 32(1):162~167. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200501022.htm LIU Jing-hua, WANG Zhu-wen. Division scheme of the geological mapping units in shallow cover areas based on ground integrated geophysical methods[J]. Geology in China, 2005, 32(1):162~167. http://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200501022.htm
[2]	任杰, 董瀚, 田建荣.浅论高密度电法测量在划分不同岩性界面中的应用[J].甘肃地质, 2010, 19(3):85~88. http://www.cnki.com.cn/Article/CJFDTOTAL-GSDZ201003017.htm REN Jie, DONG Han, TIAN Jian-rong. Application of high density electric survey on different lithological interface[J]. Gansu Geology, 2010, 19(3):85~88. http://www.cnki.com.cn/Article/CJFDTOTAL-GSDZ201003017.htm
[3]	朱首峰, 盛君.第四系覆盖区地质调查中的物探方法研究[J].江苏科技信息, 2016, (3):70~75. http://www.cnki.com.cn/Article/CJFDTOTAL-KJXY201603027.htm ZHU Shou-feng, SHENG Jun. Research on geophysical prospecting methods used in geological survey of the quaternary coverage area[J]. Jiangsu Science & Technology Information, 2016, (3):70~75. http://www.cnki.com.cn/Article/CJFDTOTAL-KJXY201603027.htm
[4]	石战结, 田钢, 薛建, 等.地震与电法技术在浅覆盖区区域地质调查中的应用[J].物探与化探, 2003, 27(3):194~201. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH200303009.htm SHI Zhan-jie, TIAN Gang, XUE Jian, et al. The application of seismic and electric methods to the regional geological survey in the shallow fields[J]. Geophysical & Geochemical Exploration, 2003, 27(3):194~201. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH200303009.htm
[5]	赵玉岩. 区域地球物理地球化学资料在浅覆盖区地质调查中的应用方法研究[D]. 长春: 吉林大学, 2009. ZHAO Yu-yan. Application of regional geophysical and geochemical data to geological survey in shallow overlay area[D]. Changchun:Jilin University, 2009.
[6]	陈浩辉. 综合物探方法在深覆盖区地质构造划分中的应用[D]. 北京: 中国地质大学(北京), 2014. CHEN Hao-hui. Comprehensive geophysical methods applied in deep geological formations coverage area is divided in[D]. Beijing:China University of Geosciences (Beijing), 2014.
[7]	卓武, 陈长敬, 陈松, 等.浅层地震反射波法和高密度电法联合在地质灾害评价中的应用[J].工程地球物理学报, 2013, 10(5):625~630. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201305007.htm ZHUO Wu, CHEN Chang-jing, CHEN Song, et al. The application of joint method of shallow seismic and high-density resistivity method to evaluation of geological disasters[J]. Chinese Journal of Engineering Geophysics, 2013, 10(5):625~630. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201305007.htm
[8]	陈长敬, 黄理善, 罗士新.隐伏断裂勘察中的综合地球物理方法研究[J].工程地球物理学报, 2011, 8(3):348~353. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201103023.htm CHEN Chang-jing, HUANG Li-shan, LUO Shi-xin. Integrated geophysical methods in buried fault exploration[J]. Chinese Journal of Engineering Geophysics, 2011, 8(3):348~353. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201103023.htm
[9]	刘保全.浅层地震勘探在岩溶调查中的应用效果[J].物探与化探, 1996, 20(6):474~477. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH606.009.htm LIU Bao-quan. The application of shallow seismic exploration to Karst invetigation[J]. Geophysical & Geochemical Exploration, 1996, 20(6):474~477. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH606.009.htm
[10]	陈蜀俊, 罗登贵, 甘家思, 等.浅层地震勘探在工程场地地震安全性评价近场工作中的应用[J].物探与化探, 2004, 28(5):463~466. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH200405025.htm CHEN Shu-jun, LUO Deng-gui, GAN Jia-si, et al. The application of shallow seismic exploration to the near-site researches on the seismic safety evaluation of the construction sites[J]. Geophysical & Geochemical Exploration, 2004, 28(5):463~466. http://www.cnki.com.cn/Article/CJFDTOTAL-WTYH200405025.htm
[11]	张胜业, 潘玉玲.应用地球物理学原理[M].武汉:中国地质大学出版社, 2004, 299~300. ZHANG Sheng-ye, PAN Yu-ling. Application Geophysics[M]. Wuhan:China University of Geosciences Press, 2004, 299~300.
[12]	刘艳章, 王水林, 葛修润, 等.高密度电法在三峡库区滑坡勘察中的应用[J].岩土力学, 2006, 27(S):519~523. http://www.cnki.com.cn/Article/CJFDTOTAL-DZHB201602019.htm LIU Zhang-yan, WANG Shui-lin, GE Xiu-run, et al. Application of high density resistivity method to landslide investigation in the three gorges reservoir area[J]. Rock and Soil Mechanics, 2006, 27(S):519~523. http://www.cnki.com.cn/Article/CJFDTOTAL-DZHB201602019.htm
[13]	刘晓东, 张虎生, 邢应太.高密度电法在岩溶地质调查中的应用[J].地质与勘探, 2003, 39(S):65~68. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200311003015.htm LIU Xiao-dong, ZHANG Hu-sheng, XING Ying-tai. The application of resistivity imaging to geological survey of Karst[J]. Geology and Prospecting, 2003, 39(S):65~68. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200311003015.htm
[14]	董浩斌, 王传雷.高密度电法的发展与应用[J].地学前缘, 2003, 10(1):171~176. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200301031.htm DONG Hao-bin, WANG Chuan-lei. Development and Application of 2D resistivity imaging surveys[J]. Earth Science Frontiers, 2003, 10(1):171~176. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200301031.htm
[15]	王志豪, 周振广, 张贺飞. EH4电磁测深在探测断层中的应用[J].资源环境与工程, 2013, 27(4):532~534. http://www.cnki.com.cn/Article/CJFDTOTAL-HBDK201304043.htm WANG Zhi-hao, ZHOU Zhen-guang, ZHANG He-fei. Application of EH4 electromagnetic sounding in faults[J]. Resources Environment & Engineering, 2013, 27(4):532~534. http://www.cnki.com.cn/Article/CJFDTOTAL-HBDK201304043.htm
[16]	蒋鑫, 黄小喜. EH4在覆盖区寻找隐伏矿床中的应用[J].矿产勘查, 2011, 2(6):784~788. http://www.cnki.com.cn/Article/CJFDTOTAL-YSJS201106019.htm JIANG Xin, HUANG Xiao-Xi. Application of EH4 system for prospecting concealed ore deposit in the overburden area[J]. Mineral Exploration, 2011, 2(6):784~788. http://www.cnki.com.cn/Article/CJFDTOTAL-YSJS201106019.htm
[17]	王玉和, 王厚臣, 程久龙, 等. EH4电磁方法在断层定位及导水性探测中的应用[J].山东科技大学学报(自然科学版), 2010, 29(4):35~39. http://www.cnki.com.cn/Article/CJFDTOTAL-SDKY201004009.htm WANG Yu-he, WANG Hou-chen, CHENG Jiu-long, et al. Application of EH4 electromagnetic method in locating of faults and detection of its water transmissibility[J]. Journal of Shandong University of Science and Technology(Natural Science), 2010, 29(4):35~39. http://www.cnki.com.cn/Article/CJFDTOTAL-SDKY201004009.htm
[18]	谭红艳, 吕骏超, 刘桂香, 等. EH4音频大地电磁测深方法在鄂东南地区寻找隐伏矿体的应用[J].地质与勘探, 2011, 47(6):1133~1141. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201106023.htm TAN Hong-yan, LV Jun-chao, LIU Gui-xiang, et al. Application of the Audio-frequency magnetotelluric method to search for concealed ore bodies in southeastern Hubei Province[J]. Geology and Exploration, 2011, 47(6):1133~1141. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKT201106023.htm
[19]	张莹, 张胜业. EH-4资料处理解释系统的研究[J].工程地球物理学报, 2005, 2(4):311~315. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ200504010.htm ZHANG Ying, ZHANG Sheng-ye. A study of the EH-4 processing and interpretation system[J]. Chinese Journal of Engineering Geophysics, 2005, 2(4):311~315. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ200504010.htm
[20]	柳建新, 蒋玲霞, 严家斌, 等. EH4电磁成像系统在高速公路隧道工程勘察中的应用[J].工程地球物理学报, 2008, 5(6):652~656. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201003013.htm LIU Jian-xin, JIANG Ling-xia, YAN Jia-bin, et al. Application of EH4 electromagnetic image system in survey of tunnels engineering for expressways[J]. Chinese Journal of Engineering Geophysics, 2008, 5(6):652~656. http://www.cnki.com.cn/Article/CJFDTOTAL-GCDQ201003013.htm

被引

资源附件(0)

访问统计

图(8)

计量

文章访问数: 223
HTML全文浏览量: 162
PDF下载量: 36
被引次数: 0

0. 引言
1. 模型建立与接触算法
1.1 层内模型
1.2 层间模型
1.3 接触问题有限元计算方法
2. 三维离层模拟算例
3. 结论

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

物探方法在强风化区填图中的应用探索

作者简介:
陈松(1985-), 男, 安徽淮南人, 工程师, 硕士, 主要从事地球物理勘察与研究工作。E-mail:anhuisongchen@163.com

通讯作者:
卜建军(1972-), 男, 高级工程师, 主要从事区域地质和地层古生物方面的研究工作。E-mail:jianjunbu@cug.edu.cn

计量