CN111781651A - 联合三种物探方法和两种地质方法的岩溶探测方法 - Google Patents
联合三种物探方法和两种地质方法的岩溶探测方法 Download PDFInfo
- Publication number
- CN111781651A CN111781651A CN202010530172.5A CN202010530172A CN111781651A CN 111781651 A CN111781651 A CN 111781651A CN 202010530172 A CN202010530172 A CN 202010530172A CN 111781651 A CN111781651 A CN 111781651A
- Authority
- CN
- China
- Prior art keywords
- karst
- apparent resistivity
- burial depth
- depth interval
- resistivity profile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 103
- 238000001514 detection method Methods 0.000 title claims abstract description 38
- 238000011161 development Methods 0.000 claims abstract description 93
- 238000009933 burial Methods 0.000 claims abstract description 66
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000010291 electrical method Methods 0.000 claims abstract description 19
- 230000002159 abnormal effect Effects 0.000 claims description 34
- 230000001052 transient effect Effects 0.000 claims description 33
- 238000009826 distribution Methods 0.000 claims description 20
- 230000002547 anomalous effect Effects 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 abstract description 17
- 230000005540 biological transmission Effects 0.000 abstract 1
- 230000009347 mechanical transmission Effects 0.000 abstract 1
- 230000018109 developmental process Effects 0.000 description 69
- 238000005259 measurement Methods 0.000 description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 14
- 230000005856 abnormality Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 238000011835 investigation Methods 0.000 description 7
- 239000011435 rock Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000013049 sediment Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 3
- 239000003673 groundwater Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241000508725 Elymus repens Species 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011981 development test Methods 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 201000008827 tuberculosis Diseases 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及页岩气钻前勘探领域，尤其是一种联合三种物探方法和两种地质方法的岩溶探测方法：S1，开展工程地质和水文地质踏勘获取综合地质信息；S2，将地下空间由浅至深划分为第一、第二和第三埋深区间，设计物探数据集成观测系统；S3，根据物探数据集成观测系统分别进行：高密度电法、TEM、AMT数据采集和数据处理反演，获取第一、第二和第三视电阻率剖面；步骤四，确定第一、第二和第三埋深区间内的岩溶发育范围；步骤五，综合第一、第二和第三埋深区间内的岩溶发育范围以及综合地质信息确定第三埋深区间以浅的岩溶发育范围。本发明拓展了采用单一物探方法探测深度局限问题，同时实现了多种物探方法的优势互补，提高了风险识别精度。
Description
技术领域
本发明涉及页岩气钻前勘探领域，尤其是一种联合三种物探方法和两种地质方法的岩溶探测方法。
背景技术
中国的页岩气主要在四川盆地，产量占全国的60％，目前页岩气产地主要集中在四川南部地区，2014年开始页岩气规模性开发，2016年已建成30亿立方米页岩气生产能力。“四川省川南地区页岩气勘查开发试验区”和“四川长宁—威远国家级页岩气示范区”已成为国家级试验区和示范区。
四川省南部属于中-高山地区，地形切割深、高差大，页岩气开发目的层一般埋深2000-3500米，钻井平台从自流井组到茅口组地层都有开孔，不同的开孔层位，平台下伏地层岩性的组合特征不同，钻进风险不同。碳酸盐岩地区岩溶极发育，具有多期性、复杂性，有充水岩溶、干枯岩溶，容易引起掉钻、卡钻事故。岩溶管道(地下暗河)发育极具隐蔽性，发生井液漏失导通性极强，对水源及环境影响巨大。
页岩气开发进程中为了规避以上风险，急需对井口到1000米深度的地层岩溶发育形态进行精细解释及预测钻井风险。目前，国内外现有技术中多采用单一的地质调查方法、地球物理勘探方法或者简单的物探+地质方法进行页岩气钻前岩溶勘探，仍然存在探测深度、识别分辨率和解释准确度方面的缺陷。
发明内容
针对现有技术方案存在的上述不足，本发明提供一种有效识别页岩气钻前勘探岩溶发育范围的联合三种物探方法和两种地质方法的岩溶探测方法。
为了实现上述目的，本发明采用的技术方案为：
联合三种物探方法和两种地质方法的岩溶探测方法，包括：
步骤一，以探测目标为中心，规划踏勘区域，在所述踏勘区域内开展工程地质踏勘和水文地质踏勘获取研究区综合地质信息；
步骤二，将所述探测目标的地下空间由浅至深划分为第一埋深区间、第二埋深区间和第三埋深区间，根据所述工程地质信息和水文地质信息设计物探数据集成观测系统；所述物探数据集成观测系统包括用于探测所述第一埋深区间的高密度电法观测系统，用于探测所述第二埋深区间的瞬变电磁法观测系统，以及用于探测所述第三埋深区间的音频大地电磁法观测系统；
步骤三，根据所述物探数据集成观测系统分别进行：
高密度电法数据采集和数据处理反演，获取第一视电阻率剖面；
瞬变电磁法数据采集和数据处理反演，获取第二视电阻率剖面；
音频大地电磁法数据采集和数据处理反演，获取第三视电阻率剖面；
步骤四，根据所述第一视电阻率剖面结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围；
步骤五，综合所述第一埋深区间内的岩溶发育范围、第二埋深区间内的岩溶发育范围和第三埋深区间内的岩溶发育范围，以及所述综合地质信息确定所述第三埋深区间以浅的岩溶发育范围。
所述探测目标为井平台。所述踏勘区域为以平台为中心1～2km2的范围。
进一步地，所述步骤四为：根据所述第一视电阻率剖面以及所述第二视电阻率剖面和所述第三视电阻率剖面中的至少一个，再结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面以及所述第三视电阻率剖面，再结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围。
进一步地，所述第一埋深区间为地表以下0～100m，所述第二埋深区间为地表以下100～400m，所述第三埋深区间为地表以下400～100m。选择100m和400m作为埋深区间划分的界限，一方面考虑了三种物探方法的探测有效深度区间，另一方面考虑了山区井平台对于三种物探方法观测系统的空间限制。因此，100m和400m的区间划分是适宜且具有施工便捷性的。
进一步地，所述综合地质信息包括地表岩溶发育地层及其产状、岩溶埋深、岩溶类型和岩溶发育主方向。所述岩溶发育主方向为三种物探方法观测系统的有效布置提供了依据，为高精度解释岩溶电性异常区提供了基础，岩溶类型对于根据三种物探方法的组合识别电性异常区提供了有力依据，有助于降低地球物理反演多解性。此外，三种物探方法组合识别地下电性异常区需要结合地下地质构造和岩性特征，综合地质信息提供的地层发育及其产状和研究区岩层信息同样为解释地下地质构造和识别电性异常区提供了依据。
进一步地，所述高密度电法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，在所述布设主方向布设3～6条测线，测线间距为20m；在所述岩溶发育主方向布设1条测线；每条测线道间距6m，至少同时布置120道。高密度电法观测系统呈十字型布设，以便于有效确认岩溶发育的实际主方向，并有效圈定岩溶发育在长轴方向和短轴方向的范围。
进一步地，所述瞬变电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设方向，布设9条测线，测线长度400m，测线间距10m，测点间距5m；所述瞬变电磁法的发射边框边长为400m，发射频率为25Hz，发射电流不低于15A。沿岩溶发育的长轴方向布设较为密集的瞬变电磁法测线，能够控制岩溶沿长轴方向的若干剖面图，以形成岩溶发育拟三维效果预测图。
进一步地，所述音频大地电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，布设3条测线：中间主测线长度为2000m，两侧辅助测线长度为1000m，点距40米，线距40米；在所述岩溶发育主方向布设1条测线，测线长度2000m，点距40m；所述音频大地电磁法的采集最低频率为0.1Hz，采集时间大于30min。音频大地电磁法观测系统也是呈十字型布置，以获取岩溶沿长轴方向和短轴方向的较深深度范围的电性剖面信息。
进一步地，步骤五还包括：确认地表岩溶类型，在所述第一视电阻率剖面、第二视电阻率剖面和第三视电阻率剖面中分别识别对应该地表岩溶的第一异常区、第二异常区和第三异常区，根据所述第一异常区、第二异常区和第三异常区确认该地表岩溶的发育范围。三种物探方法结合两种地质方法综合确认地下岩溶异常时，对于有地表岩溶出露的目标，首先从地表地质岩溶信息出发，在三种物探方法获取的视电阻率剖面中分别识别出对应该地表岩溶的地下电性异常区，如此能够在保证电性剖面解释以工区地质概况为依据的基础上，合理且全深度范围地识别出该地表岩溶对应的地下岩溶发育情况。
进一步地，步骤五还包括：首先通过第一视电阻率剖面识别所述第一异常区，再通过第二视电阻率剖面识别所述第二异常区，最后通过所述第三视电阻率剖面识别所述第三异常区。由于高密度电法具有有效探测范围较浅但反映异常精度较高的特点，以地表出露岩溶为基础，由浅至深分别利用高密度电法、瞬变电磁法和音频大地电磁法获取地表岩溶在地下的分布范围，能够进一步提升地下岩溶分布的识别准确度，有效预测钻井风险。
进一步地，根据所述工程地质信息和/或通过视电阻率剖面确认该地表岩溶在剖面的发育范围之后，通过输出沿岩溶发育主方向的连续瞬变电磁视电阻率剖面图，得到该地表岩溶的立体分布范围。由于瞬变电磁法采用均匀网度观测系统，能够快速进行平面网度测量，能够有效地控制岩溶异常在平面和深度上的发育情况。
进一步地，若岩溶类型为干枯岩溶，则根据所述第一视电阻率剖面和所述第三电阻率剖面重合的高阻异常区识别为干枯岩溶发育范围。瞬变电磁法能够有效探测高阻异常中的低阻异常，但对于低阻异常中的高阻异常识别效果低于高密度电法。因此，对于干枯岩溶等高阻异常，采用高密度电法结合音频大地电磁法重合剖面进行对比分析识别，能够达到更好的识别准确度。
本发明的有益效果：
本发明提供联合三种物探方法和两种地质方法的岩溶探测方法实现了采用三种物探方法在由浅至深的三段埋深区间分段识别地下岩溶异常区，拓展了采用单一物探方法进行岩溶勘探的探测深度局限问题，实现了满足页岩气钻井钻前风险预测的风险识别范围，降低了相关钻探风险。同时，三种物探方法之间在识别准确度上形成互补，相互印证，进一步提高了对于岩溶探测的准确度。
此外，综合工程地质踏勘成果和水文地质踏勘成果一方面对三种物探方法的观测系统设计提供了理论支撑，也对三种物探方法的成果解译提供了依据，降低了三种物探方法的反演多解性，最终实现了页岩气钻井预测下伏1000m以浅的地层钻井风险识别。
附图说明
图1为本发明实施例1的物探数据集成观测系统示意图；
图2为本发明实施例1高密度电法勘探流程图；
图3为本发明实施例1高密度电法数据处理流程图；
图4为本发明实施例1瞬变电磁法大定源回线装置示意图；
图5为本发明实施例1瞬变电磁法数据处理流程图；
图6为本发明实施例1音频大地电磁野外工作示意图；
图7为本发明实施例1音频大地电磁法数据处理流程图；
图8为本发明实施例1的物探解释成果一；
图9为本发明实施例1的物探解释成果二；
图10为本发明实施例1的物探解释成果三。
具体实施方式
下面结合附图，对本发明作详细的说明。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
实施例1
实施例1提供一种联合三种物探方法和两种地质方法的岩溶探测方法，包括：
步骤一，以探测目标为中心，规划踏勘区域，在所述踏勘区域内开展工程地质踏勘和水文地质踏勘获取研究区综合地质信息；
步骤二，将所述探测目标的地下空间由浅至深划分为第一埋深区间、第二埋深区间和第三埋深区间，根据所述工程地质信息和水文地质信息设计物探数据集成观测系统；所述物探数据集成观测系统包括用于探测所述第一埋深区间的高密度电法观测系统，用于探测所述第二埋深区间的瞬变电磁法观测系统，以及用于探测所述第三埋深区间的音频大地电磁法观测系统；
步骤三，根据所述物探数据集成观测系统分别进行：
高密度电法数据采集和数据处理反演，获取第一视电阻率剖面；
瞬变电磁法数据采集和数据处理反演，获取第二视电阻率剖面；
音频大地电磁法数据采集和数据处理反演，获取第三视电阻率剖面；
步骤四，根据所述第一视电阻率剖面结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围；
步骤五，综合所述第一埋深区间内的岩溶发育范围、第二埋深区间内的岩溶发育范围和第三埋深区间内的岩溶发育范围，以及所述综合地质信息确定所述第三埋深区间以浅的岩溶发育范围。
所述探测目标为井平台。所述踏勘区域为以平台为中心1～2km2的范围。
所述第一埋深区间为地表以下0～100m，所述第二埋深区间为地表以下100～400m，所述第三埋深区间为地表以下400～100m。选择100m和400m作为埋深区间划分的界限，一方面考虑了三种物探方法的探测有效深度区间，另一方面考虑了山区井平台对于三种物探方法观测系统的空间限制。因此，100m和400m的区间划分是适宜且具有施工便捷性的。
所述综合地质信息包括地表岩溶发育地层及其产状、岩溶埋深、岩溶类型和岩溶发育主方向。所述岩溶发育主方向为三种物探方法观测系统的有效布置提供了依据，为高精度解释岩溶电性异常区提供了基础，岩溶类型对于根据三种物探方法的组合识别电性异常区提供了有力依据，有助于降低地球物理反演多解性。此外，三种物探方法组合识别地下电性异常区需要结合地下地质构造和岩性特征，综合地质信息提供的地层发育及其产状和研究区岩层信息同样为解释地下地质构造和识别电性异常区提供了依据。
所述工程地质踏勘具体包括以下内容：
①确定探测目标周围地层划分，开始地质测量；②在野外对岩石进行观测，研究岩层顺序，成因、结构、产状、厚度、岩性和褶皱类型，断层性质等。且对每一岩层特征详细描述，如：颜色、层次、组成成分，颗粒形状与大小，分选程度，岩相变化，化石的种属和分布情况，岩石的胶结程度。结核的性质和分布，上下地层的接触关系，岩石节理，裂隙的形状和发育程度，矿化程度及风化现象等；③进行地质构造测绘；④通过地面观察，研究岩溶发育与岩性、层理、褶皱、断裂破碎带，以及新的构造关系等；⑤调查溶洞大小、形状、高程、育填物以及其它各种溶蚀现象的特征；⑥对每种典型的地貌单元都应进行形态描述和测量，如：阶地宽度、阶地台面倾斜方向及角度，阶地前缘高度；冲沟宽度(顶宽、底宽)及深度，山坡倾斜方向及坡度，山脊走向，相对高差等；⑦确定第四纪沉积物的分布规律及成因类型，查明不同成因类型第四纪沉积物的地层层序及时代，了解第四纪沉积物的成分，结构及厚度变化，确定岩相变化的规律及第四纪地质发展史，从水文地质角度出发，查明第四纪沉积物的水文地质特征。
所述水文地质踏勘具体包括以下内容：
①对探测目标周围地下水类型、含水层、隔水层调查主要包括调查地下水类型、调查城镇及工矿地下水水源地的位置和用途、机井及民井调查、泉的调查；②对探测目标周围地下水补给、径流、排泄调查主要包括调查地下水的补给来源、补给方式或途径，补给区分布和补给量；地下水的径流条件、径流分带规律和流向；地下水的排泄形式、排泄途径和排泄区(带)分布；不同含水层之间、地下水和地表水之间水力联系。
如图1所示为物探数据集成观测系统，其中：所述高密度电法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，在所述布设主方向布设3～6条测线(图1中G起头的测线)，测线间距为20m；在所述岩溶发育主方向布设1条测线；每条测线道间距5m，至少同时布置120道。高密度电法观测系统呈十字型布设，以便于有效确认岩溶发育的实际主方向，并有效圈定岩溶发育在长轴方向和短轴方向的范围。
高密度电法基本原理为对称四极测深，为达到控制深度100米，道间距5米，至少同时布置120道，即最多能测到39层数据，每层探测深度为道间距的0.6～0.7倍，因此测线长度一般布置600-800米，线距20米。
如图2所示，高密度电阻率法现场数据采集时，为达到探测100米深度，同时布设120根及以上电极，并用多心电缆将所有电极与多路电极转换器连接。测量时，通过测量主机按照预先设定好道间距5米、测量层数39层及以上、连续滚动测量、采样周期200毫秒的参数，将指令传递给多路电极转换器，所述多路电极转换器根据指令控制各电极的通断，从而实现测量装置、极距以及测点的自动转换；测量主机通过接入回路的电极对地下供电并进行数据采集，将采集到的数据以一定的顺序存储在高密度主机仪器内部；仪器内部存储的数据，再转存到计算机中，在计算机中利用完成对数据的处理，生成反演断面图；解释人员根据反演断面图结合所述综合地质信息进行定性或半定量解释，得出地下0-100米深度岩溶发育位置、影响范围相关地质问题。
如图3所示为计算机中进行的高密度电法数据处理流程图，包括以下步骤：对原始数据剔除坏数据后进行插值处理；校正插值后数据的畸变数据；建立地形模型；根据地形模型得到正演理论值；采用有限差分法、有限单元法和边界元法进行正演模拟，并对比分析选择效果最佳的正演方法完成正演；反演迭代输出高密度电法视电阻率剖面图。
所述瞬变电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设方向，布设9条测线(图1中S起头的测线)，测线长度400m，测线间距10m，测点间距5m；所述瞬变电磁法的发射边框边长为400m，发射频率为25Hz，发射电流不低于15A。沿岩溶发育的长轴方向布设较为密集的瞬变电磁法测线，能够控制岩溶沿长轴方向的若干剖面图，以形成岩溶发育拟三维效果预测图。
所述瞬变电磁法采用不接地回线向地下发送脉冲电磁场(一次场)，用接收线圈来观测由一次场感应而产生的二次电磁场，通过观测二次场的变化规律来研究地下介质的电性分布特征，从而探测地下100-400米深度岩溶发育位置、影响范围；
瞬变电磁法探测深度由发射边框长度、发射频率、发射电流决定，为探测地下400米深度，发射边框边长400米，发射频率25H，发射电流不低于15安培；
接收框边长1米，接收框匝数为100匝，等效接收线圈接收面积为100m2；
采集方式采用大定源回线装置，如图4所示，在探测区首先布置400为边长的发射框，接收线圈在发射框内1/3区域进行逐点测量，采集地下地质体因发射框一次场感应而产生的二次电磁场电信号，数据处理可分为四个步骤：(1)原始数据的录入；(2)一次数据处理：包括数据编辑，数据滤波、均滑处理；(3)二次数据处理：包括建立反演文件，数据经Surfer处理成图视电阻率断面图，解译地下100-400米深度岩溶发育位置、影响范围相关地质问题。
如图5所示为瞬变电磁法数据处理流程图，包括以下步骤：根据工程测量资料和归一化感应电压对原始数据进行整理；数据空间属性建立；选择测道范围、畸变数据剔除和滤波；绘制多测道曲线，反演视电阻率拟断面图，并判定所述视电阻率拟断面图是否满足成果解释要求；若是，绘制解释用多测道曲线和视电阻率拟断面图。
所述音频大地电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，布设3条测线(图1中L起头的测线)：中间主测线长度为2000m，两侧辅助测线长度为1000m，点距40米，线距20米；在所述岩溶发育主方向布设1条测线，测线长度2000m，点距40m；所述音频大地电磁法的采集最低频率为0.1Hz，采集时间大于30min。音频大地电磁法观测系统也是呈十字型布置，以获取岩溶沿长轴方向和短轴方向的较深深度范围的电性剖面信息。所述音频大地电磁法勘探野外工作示意图如图6所示。图7为音频大地电磁法数据处理流程图，包括以下步骤：对原始数据进行静态校正和预处理；对预处理后的数据根据钻孔资料和地质资料进行定性分析和质量检查；对通过质量检查的数据进行反演输出视电阻率剖面图，对所述视电阻率剖面图进行电性参数定量分析，划分研究区不同地质体对应的电性参数区间范围；结合高密度电法视电阻率剖面、瞬变电磁法视电阻率剖面、综合地质信息中的任意一项进行综合分析。
步骤五还包括：确认地表岩溶类型，在所述第一视电阻率剖面、第二视电阻率剖面和第三视电阻率剖面中分别识别对应该地表岩溶的第一异常区、第二异常区和第三异常区，根据所述第一异常区、第二异常区和第三异常区确认该地表岩溶的发育范围。三种物探方法结合两种地质方法综合确认地下岩溶异常时，对于有地表岩溶出露的目标，首先从地表地质岩溶信息出发，在三种物探方法获取的视电阻率剖面中分别识别出对应该地表岩溶的地下电性异常区，如此能够在保证电性剖面解释以工区地质概况为依据的基础上，合理且全深度范围地识别出该地表岩溶对应的地下岩溶发育情况。
步骤五还包括：首先通过第一视电阻率剖面识别所述第一异常区，再通过第二视电阻率剖面识别所述第二异常区，最后通过所述第三视电阻率剖面识别所述第三异常区。由于高密度电法具有有效探测范围较浅但反映异常精度较高的特点，以地表出露岩溶为基础，由浅至深分别利用高密度电法、瞬变电磁法和音频大地电磁法获取地表岩溶在地下的分布范围，能够进一步提升地下岩溶分布的识别准确度，有效预测钻井风险。
优选的实施方式中，所述步骤四为：根据所述第一视电阻率剖面以及所述第二视电阻率剖面和所述第三视电阻率剖面中的至少一个，再结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面以及所述第三视电阻率剖面，再结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围。
如图8所示为工程地质及水文地质+高密度电法+瞬变电磁法(TEM)+音频大地电磁法(AMT)探测含水裂隙的成果图。此裂隙在地表的可见宽度约0.5米，其间被松散物质充填。该剖面从上至下分别为属于同一地表裂隙上的高密度电法、瞬变电磁法及音频大地电磁测深法视电阻率反演剖面图，图中虚线圈出区域即为含水裂隙分布区。
通过高密度电法反演成果，发现在地表裂隙下方存在一处视电阻率低阻异常，判定为裂隙深部含水引起，该裂隙影响深度约地下120米，通过高密度电法、瞬变电磁法与音频大地电磁法相重合剖面对比分析确认出地下裂隙含水带在深部沿地层走向发育。
根据所述工程地质信息和/或通过视电阻率剖面确认该地表岩溶在剖面的发育范围之后，通过输出沿岩溶发育主方向的连续瞬变电磁视电阻率剖面图，得到该地表岩溶的立体分布范围。由于瞬变电磁法采用均匀网度观测系统，能够快速进行平面网度测量，能够有效地控制岩溶异常在平面和深度上的发育情况。如图9所示为含水裂隙瞬变电磁发反演成果连续视电阻率剖面图，图中虚线圈出区域即为含水裂隙分布区。
若岩溶类型为干枯岩溶，则根据所述第一视电阻率剖面和所述第三电阻率剖面重合的高阻异常区识别为干枯岩溶发育范围。瞬变电磁法能够有效探测高阻异常中的低阻异常，但对于低阻异常中的高阻异常识别效果低于高密度电法。因此，对于干枯岩溶等高阻异常，采用高密度电法结合音频大地电磁法重合剖面进行对比分析识别，能够达到更好的识别准确度。如图10所示，通过高密度电法与音频大地电磁法相重合剖面(上两幅剖面)对比分析可以得出井台下方干枯岩溶较为发育，两种方法吻合性较好。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.联合三种物探方法和两种地质方法的岩溶探测方法，其特征在于，包括：
步骤一，以探测目标为中心，规划踏勘区域，在所述踏勘区域内开展工程地质踏勘和水文地质踏勘获取研究区综合地质信息；
步骤二，将所述探测目标的地下空间由浅至深划分为第一埋深区间、第二埋深区间和第三埋深区间，根据所述工程地质信息和水文地质信息设计物探数据集成观测系统；所述物探数据集成观测系统包括用于探测所述第一埋深区间的高密度电法观测系统，用于探测所述第二埋深区间的瞬变电磁法观测系统，以及用于探测所述第三埋深区间的音频大地电磁法观测系统；
步骤三，根据所述物探数据集成观测系统分别进行：
高密度电法数据采集和数据处理反演，获取第一视电阻率剖面；
瞬变电磁法数据采集和数据处理反演，获取第二视电阻率剖面；
音频大地电磁法数据采集和数据处理反演，获取第三视电阻率剖面；
步骤四，根据所述第一视电阻率剖面结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围；
步骤五，综合所述第一埋深区间内的岩溶发育范围、第二埋深区间内的岩溶发育范围和第三埋深区间内的岩溶发育范围，以及所述综合地质信息确定所述第三埋深区间以浅的岩溶发育范围。
2.根据权利要求1所述的方法，其特征在于：所述第一埋深区间为地表以下0～100m，所述第二埋深区间为地表以下100～400m，所述第三埋深区间为地表以下400～100m；
所述步骤四为：根据所述第一视电阻率剖面以及所述第二视电阻率剖面和所述第三视电阻率剖面中的至少一个，再结合所述综合地质信息确定所述第一埋深区间内的岩溶发育范围；
根据所述第二视电阻率剖面以及所述第三视电阻率剖面，再结合所述综合地质信息确定所述第二埋深区间内的岩溶发育范围；
根据所述第三视电阻率剖面结合所述综合地质信息确定所述第三埋深区间内的岩溶发育范围。
3.根据权利要求1或2所述的方法，其特征在于：所述综合地质信息包括地表岩溶发育地层及其产状、岩溶埋深、岩溶类型和岩溶发育主方向。
4.根据权利要求3所述的方法，其特征在于：所述高密度电法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，在所述布设主方向布设3～6条测线，测线间距为20m；在所述岩溶发育主方向布设1条测线；每条测线道间距6m，至少同时布置120道。
5.根据权利要求3所述的方法，其特征在于：所述瞬变电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设方向，布设9条测线，测线长度400m，测线间距10m，测点间距5m；所述瞬变电磁法的发射边框边长为400m，发射频率为25Hz，发射电流不低于15A。
6.根据权利要求3所述的方法，其特征在于：所述音频大地电磁法观测系统的布设方式为：以探测目标为中心，以所述岩溶发育主方向的垂直方向为布设主方向，布设3条测线：中间主测线长度为2000m，两侧辅助测线长度为1000m，点距40米，线距40米；在所述岩溶发育主方向布设1条测线，测线长度2000m，点距40m；所述音频大地电磁法的采集最低频率为0.1Hz，采集时间大于30min。
7.根据权利要求3所述的方法，其特征在于，步骤五还包括：确认地表岩溶类型，在所述第一视电阻率剖面、第二视电阻率剖面和第三视电阻率剖面中分别识别对应该地表岩溶的第一异常区、第二异常区和第三异常区，根据所述第一异常区、第二异常区和第三异常区确认该地表岩溶的发育范围。
8.根据权利要求7所述的方法，其特征在于，步骤五还包括：首先通过第一视电阻率剖面识别所述第一异常区，再通过第二视电阻率剖面识别所述第二异常区，最后通过所述第三视电阻率剖面识别所述第三异常区。
9.根据权利要求7所述的方法，其特征在于：步骤五还包括：根据所述工程地质信息和/或通过视电阻率剖面确认该地表岩溶在剖面的发育范围之后，通过输出沿岩溶发育主方向的连续瞬变电磁视电阻率剖面图，得到该地表岩溶的立体分布范围。
10.根据权利要求7所述的方法，其特征在于：若岩溶类型为干枯岩溶，则根据所述第一视电阻率剖面和所述第三电阻率剖面重合的高阻异常区识别为干枯岩溶发育范围。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010530172.5A CN111781651B (zh) | 2020-06-11 | 2020-06-11 | 联合三种物探方法和两种地质方法的岩溶探测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010530172.5A CN111781651B (zh) | 2020-06-11 | 2020-06-11 | 联合三种物探方法和两种地质方法的岩溶探测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111781651A true CN111781651A (zh) | 2020-10-16 |
| CN111781651B CN111781651B (zh) | 2022-08-05 |
Family
ID=72756191
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010530172.5A Active CN111781651B (zh) | 2020-06-11 | 2020-06-11 | 联合三种物探方法和两种地质方法的岩溶探测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111781651B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112305624A (zh) * | 2020-10-30 | 2021-02-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种非接触式跨障碍的高密度测量方法 |
| CN112363245A (zh) * | 2020-10-26 | 2021-02-12 | 王昌益 | 一种地下水物探测量方法及系统 |
| CN112817057A (zh) * | 2020-12-31 | 2021-05-18 | 中国地质调查局天津地质调查中心 | 经济快速精准探测填埋坑塘的地下空间展布特征的方法 |
| CN113049783A (zh) * | 2021-03-23 | 2021-06-29 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | 一种地质数据获取方法 |
| CN113420456A (zh) * | 2021-07-07 | 2021-09-21 | 核工业航测遥感中心 | 基于反演电阻率断面的物探地质数据库合并方法 |
| CN114386154A (zh) * | 2022-01-14 | 2022-04-22 | 中水珠江规划勘测设计有限公司 | 坝基现代岩溶层影响深度确定的方法、装置及电子设备 |
| CN115016017A (zh) * | 2022-07-15 | 2022-09-06 | 中国地质调查局油气资源调查中心 | 一种泥页岩层系识别方法及装置 |
| CN117872494A (zh) * | 2024-03-12 | 2024-04-12 | 北京科技大学 | 深部金属矿花岗岩岩石内部节理信息获取方法、预测方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103064120A (zh) * | 2012-12-29 | 2013-04-24 | 福州华虹智能科技开发有限公司 | 煤矿井下磁电综合探测仪及磁电综合探测方法 |
| CN107305257A (zh) * | 2016-04-21 | 2017-10-31 | 新疆维吾尔自治区煤炭科学研究所 | 高密度电阻率法与瞬变电磁法联合反演技术 |
| CN108345049A (zh) * | 2018-02-12 | 2018-07-31 | 山东大学 | 地下工程不良地质探测多方法约束反演与联合解释方法 |
| CN109613616A (zh) * | 2018-12-20 | 2019-04-12 | 中交第三航务工程局有限公司 | 一种浅层岩溶隧道地质综合立体探测方法 |
| CN110823962A (zh) * | 2019-11-14 | 2020-02-21 | 山东大学 | 一种滑坡体的三维成像方法及系统 |
| CN111007578A (zh) * | 2019-12-06 | 2020-04-14 | 山东滨莱高速公路有限公司 | 一种路基岩溶及采空区的综合探测方法及系统 |
-
2020
- 2020-06-11 CN CN202010530172.5A patent/CN111781651B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103064120A (zh) * | 2012-12-29 | 2013-04-24 | 福州华虹智能科技开发有限公司 | 煤矿井下磁电综合探测仪及磁电综合探测方法 |
| CN107305257A (zh) * | 2016-04-21 | 2017-10-31 | 新疆维吾尔自治区煤炭科学研究所 | 高密度电阻率法与瞬变电磁法联合反演技术 |
| CN108345049A (zh) * | 2018-02-12 | 2018-07-31 | 山东大学 | 地下工程不良地质探测多方法约束反演与联合解释方法 |
| CN109613616A (zh) * | 2018-12-20 | 2019-04-12 | 中交第三航务工程局有限公司 | 一种浅层岩溶隧道地质综合立体探测方法 |
| CN110823962A (zh) * | 2019-11-14 | 2020-02-21 | 山东大学 | 一种滑坡体的三维成像方法及系统 |
| CN111007578A (zh) * | 2019-12-06 | 2020-04-14 | 山东滨莱高速公路有限公司 | 一种路基岩溶及采空区的综合探测方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 朱自强: "综合物探在隧道岩溶水害勘察中的应用", 《建筑技术开发》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112363245A (zh) * | 2020-10-26 | 2021-02-12 | 王昌益 | 一种地下水物探测量方法及系统 |
| CN112305624A (zh) * | 2020-10-30 | 2021-02-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种非接触式跨障碍的高密度测量方法 |
| CN112305624B (zh) * | 2020-10-30 | 2024-02-06 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种非接触式跨障碍的高密度测量方法 |
| CN112817057A (zh) * | 2020-12-31 | 2021-05-18 | 中国地质调查局天津地质调查中心 | 经济快速精准探测填埋坑塘的地下空间展布特征的方法 |
| CN113049783A (zh) * | 2021-03-23 | 2021-06-29 | 中国地质调查局南京地质调查中心(华东地质科技创新中心) | 一种地质数据获取方法 |
| CN113420456A (zh) * | 2021-07-07 | 2021-09-21 | 核工业航测遥感中心 | 基于反演电阻率断面的物探地质数据库合并方法 |
| CN113420456B (zh) * | 2021-07-07 | 2022-05-03 | 核工业航测遥感中心 | 基于反演电阻率断面的物探地质数据库合并方法 |
| CN114386154A (zh) * | 2022-01-14 | 2022-04-22 | 中水珠江规划勘测设计有限公司 | 坝基现代岩溶层影响深度确定的方法、装置及电子设备 |
| CN114386154B (zh) * | 2022-01-14 | 2023-12-15 | 中水珠江规划勘测设计有限公司 | 坝基现代岩溶层影响深度确定的方法、装置及电子设备 |
| CN115016017A (zh) * | 2022-07-15 | 2022-09-06 | 中国地质调查局油气资源调查中心 | 一种泥页岩层系识别方法及装置 |
| CN117872494A (zh) * | 2024-03-12 | 2024-04-12 | 北京科技大学 | 深部金属矿花岗岩岩石内部节理信息获取方法、预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111781651B (zh) | 2022-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111781651B (zh) | 联合三种物探方法和两种地质方法的岩溶探测方法 | |
| Dahlin | The development of DC resistivity imaging techniques | |
| Sharlov et al. | Transient electromagnetic surveys for highresolution near-surface exploration: basics and case studies | |
| CN111045114B (zh) | 一种玄武岩覆盖区砂岩型铀矿成矿有利砂体识别定位方法 | |
| CN110018528B (zh) | 幔源油气综合地质分析方法 | |
| Bernstone et al. | Assessment of two automated electrical resistivity data acquisition systems for landfill location surveys: two case studies | |
| Dahlin | The development of electrical imaging techniques | |
| Gao et al. | Water detection within the working face of an underground coal mine using 3D electric resistivity tomography (ERT) | |
| Raji | Review of electrical and gravity methods of near-surface exploration for groundwater | |
| Su et al. | Detection method of pile foundation on subway lines based on cross-hole resistivity computed tomography | |
| Negri et al. | INTEGRATED ANALYSIS OF GEOLOGICAL AND GEOPHYSICAL DATA FOR THE DETECTION OF UNDERGROUND MAN-MADE CAVES IN AN AREA IN SOUTHERN ITALY. | |
| Longo et al. | Prospecting for clay minerals within volcanic successions: application of electrical resistivity tomography to characterise bentonite deposits in northern Sardinia (Italy) | |
| Asry et al. | Groundwater exploration using 2-D geoelectrical resistivity imaging technique at Sungai. Udang, Melaka | |
| Kingdon et al. | UKGEOS Cheshire Energy Research Field Site: science infrastructure: version 2 | |
| Ramalho et al. | Using geophysical methods to characterize an abandoned uranium mining site, Portugal | |
| Kelly et al. | Continuous electrical imaging for mapping aquifer recharge along reaches of the Namoi River in Australia | |
| Zhou et al. | Multi-turn small-loop transient electromagnetic data processing using constraints from borehole and electrical resistivity tomography data | |
| Kiflu | Improved 2D and 3D resistivity surveys using buried electrodes and optimized arrays: The multi-electrode resistivity implant technique (MERIT) | |
| Feofilaktov et al. | The deep structure of the Pauzhetka hydrothermal system area, Southern Kamchatka | |
| Amini et al. | Anomaly enhancement in 2D electrical resistivity imaging method using a residual resistivity technique | |
| Majzoub | Characterization and delineation of karst geohazards along RM652 using electrical resistivity tomography, Culberson County, Texas | |
| Al-Hadithi et al. | Using source parameter imaging technique to the aeromagnetic data to estimate the basement depth of Tharthar Lake and surrounding area in Central Iraq | |
| Zhang et al. | A novel three-dimensional mine area hydrogeological model based on groundwater modeling systems | |
| Al Amoush et al. | Using geophysical methods to image near-surface cylindrical pipeline: A case study on engineering applications, Jordan | |
| Robison | Delineating Subsurface Karst Potential Using Electrical Resistivity in the Shell Mountain Province, Fort Hood Military Installation, Texas |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: No.39 Qinghua Road, Qingyang District, Chengdu, Sichuan 610072 Patentee after: Sichuan Natural Resources Investment Group Geophysical Exploration Institute Co.,Ltd. Address before: 610000 No. 39, Qinghua Road, Qingyang District, Chengdu, Sichuan Patentee before: SICHUAN ZHONGCHENG COALFIELD GEOPHYSICAL ENGINEERING INSTITUTE CO.,LTD. |
|
| CP03 | Change of name, title or address |