CN1252492C - 大地电磁波电阻率测量方法及其仪器 - Google Patents
大地电磁波电阻率测量方法及其仪器 Download PDFInfo
- Publication number
- CN1252492C CN1252492C CNB2003101130358A CN200310113035A CN1252492C CN 1252492 C CN1252492 C CN 1252492C CN B2003101130358 A CNB2003101130358 A CN B2003101130358A CN 200310113035 A CN200310113035 A CN 200310113035A CN 1252492 C CN1252492 C CN 1252492C
- Authority
- CN
- China
- Prior art keywords
- depth
- degree
- resistivity
- curve
- data
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000005684 electric field Effects 0.000 claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 9
- 230000008859 change Effects 0.000 claims abstract description 7
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 238000005553 drilling Methods 0.000 claims description 31
- 238000005070 sampling Methods 0.000 claims description 22
- 238000013500 data storage Methods 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 17
- 241000209094 Oryza Species 0.000 claims description 16
- 235000007164 Oryza sativa Nutrition 0.000 claims description 16
- 235000009566 rice Nutrition 0.000 claims description 16
- 230000003321 amplification Effects 0.000 claims description 13
- 238000011835 investigation Methods 0.000 claims description 13
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000000691 measurement method Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 7
- 239000002184 metal Substances 0.000 claims description 5
- 229910052751 metal Inorganic materials 0.000 claims description 5
- 239000003245 coal Substances 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 4
- 238000007689 inspection Methods 0.000 claims description 4
- 230000000644 propagated effect Effects 0.000 claims description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 229910052802 copper Inorganic materials 0.000 claims description 3
- 239000010949 copper Substances 0.000 claims description 3
- 230000003412 degenerative effect Effects 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 239000003102 growth factor Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 3
- 239000011707 mineral Substances 0.000 abstract description 3
- 230000008901 benefit Effects 0.000 abstract description 2
- 239000002344 surface layer Substances 0.000 abstract 1
- 230000001276 controlling effect Effects 0.000 description 13
- 239000010410 layer Substances 0.000 description 7
- 230000001105 regulatory effect Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 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/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是大地电磁波电阻率测量方法和仪器。方法是建立线性坐标观察***，预置基本采集参数；进行深度系数的标定、表层深度系数的校正、确定采集开始/结束深度、进行通道增益设置和低通滤波器设置；仪器包括接收电场强度信号的电场传感器，接收磁场强度信号的磁场传感器，两个分别与电场传感器和磁场传感器的输出端连接的前置放大与滤波器、与前置放大与滤波器的输出端连接的数据采集***和采集控制、数据存储及处理***。该方法由实际测量数据出发建立线性坐标观察***，使大地电磁测深测量精度大大提高，直接应用予矿藏勘探，减少部分钻井，提高勘探整体效益；仪器实时测量随深度变化的电场强度、磁场强度曲线；利用仪器对采集资料进行现场处理，获得随深度变化的电磁波电阻率曲线。
Description
技术领域
本发明涉及一种地球物理勘探的方法和仪器，具体地说，涉及一种直接将频率转变为深度的大地电磁波电阻率测量方法及其仪器。
背景技术
传统的大地电磁测深法始于上世纪五十年代，距今已有五十多年历史。它是基于电磁感应原理，用于研究地球电性的一种地球物理方法。在研究地球电性的多种电法和电磁法中，大地电磁测深法具有探测深度大，不受高阻层屏蔽，分辨能力较强，等值范围较窄，工作成本低廉和野外装备较轻便等特点。在很大程度上，这是由于大地电磁测深法利用频带非常宽阔的，能量巨大的天然交变电磁作场源的结果。大地电磁测深法建立在电磁波波动方程的基础上，它的标量阻抗理论的假设前提是平面电磁波垂直入射到地球表面，地层为水平层状介质。当地层为水平非均匀构造及存在各向异性时，地层阻抗随时间及测量方向而变化，这时采用大地电磁测深张量阻抗理论。根据大地电磁测深基本理论，苏、美、法、德、加等国家研制了各种大地电磁测深仪器，我国也研制了大地电磁测深仪器。根据大地测深理论，提出了对野外观测资料的解释方法。其中定量解释主要根据视电阻率曲线反演地下的电性分层。大地电磁测深实测曲线的电子计算机自动反演，通常根据给定的初始参数(Hi ρi)计算出理论曲线，然后与一条实测曲线进行对比，按最小方差原理不断修改参数，找出一条与实测曲线拟合最好的理论曲线，这条理论曲线所对应的参数就是反演的解。后来又提出了一维连续模型的大地电磁测深反演方法。把地球介质的电阻率视为各向同性随深度连续变化的函数。它仍然是从初始模型出发，计算理论曲线，与实测曲线对比，修改初始模型，逐次逼近实测曲线，最后获得经过唯一性评价的电阻率随深度连续变化的曲线。大地电磁测深法应用于地壳上地幔研究，普查勘探沉积盆地和含油气构造，寻找和勘探热田，进行地震预报和寻找金属矿产，地下水源等方面。
虽然大地电磁测深法已经得到较广泛应用。但是由于它的勘探精度不高，其分层厚度往往是以百米计，甚致数千米，很少见到数十米的情况。因此它在矿藏勘探中难于发挥作用。此外其勘探效率有待进一步提高。
发明内容
本发明的目的是克服上述缺陷，提供一种探测深度大，不受高阻层屏蔽，高精度、高效率的大地电磁波电阻率测量方法加强仪器。
本发明所述的大地电磁波电阻率测量方法，其步骤是：
第一步，建立线性坐标观察***，预置初始采集参数：
进入一个勘探区，首先建立线性坐标观察***，在这个***中频率与深度的关系表示为：
F＝a+bK′/H (1)
式中：a-为表层系数，是一个较小的值，在第一步工作中可忽略
b-为区域性地层电阻率随深度的增长系数，其值设定为0.01
H-为深度(单位米)
K′-为电磁波在地层中传播定义的探测深度值
在上式中H、H1、H2、S单位为米，F单位为赫兹(HZ)
(1)式可以进-步表示为：
F＝a+B/H
B定义为深度系数，在第一步中，B＝bK′
测量时，分别输入设计要求的起始测量深度值(H1)、终止测量深度值(H2)、采样间隔值(S)，将由关系式确定的F参数预置到电阻率测量仪器各参数界面，使仪器正常运行。
第二步，进行深度系数B的标定：
进入采集后，首先应当在探区或相邻地区，选择一口已知钻井，利用第一步基本关系式确定的参数，在选择的已知钻井旁探测，测取一条电磁波电阻率曲线，用这条曲线与已知钻井的测井电阻率曲线进行对比，如此反复，直到找到一条相关性最好的曲线为止，按下列步骤进行B参数标定：
选择已知钻井的条件是：(1)地层有代表性；(2)有电阻率测井资料或钻井取芯等资料；(3)地层倾角≤15度；(4)地面无强工业电网等干扰。
将测得的随深度变化的电阻率曲线，截取一段与已知钻井中的电测电阻率曲线或者具有电性特征的其他地层分层资料进行比较，求得采样间隔S′(单位米)，该S′用下式表示：
S′＝(Hp2-Hp1)/(L2-L1) (2)
式中：Hp1-为已钻探井的第1个测井曲线特征点的深度，单位米
Hp2-为已钻探井的第2个测井曲线特征点的深度，单位米
L1-为新测探的与第1个测井特征点相似特征的电磁波电阻率曲线的采集点号数
L2-为新测探的与第2个测井特征点相似特征的电磁波电阻率曲线的采集点号数
B＝S′/S (3)
B为实际确定的深度系数，在以后的工作中用B取代(1)式中的bK值进行频率深度关系变换；在没有深度标定井的新区，可以利用以前的与这个区块类似的其他区域的深度系数，对基本关系式进行标定后进行测量。
第三步，进行表层深度系数Ha的校正：
利用经过第二步参数校正的曲线与已知钻井电测井电阻率曲线进行全测量段比较，各特征深度点之间存在的***误差值用下式表示：
Ha＝Hd-HC (4)
Hd--为已知钻井电测井曲线特征点深度
Hc--为探测电磁波电阻率曲线同一特征点深度
式中Ha、Hd、Hc单位为米；
利用Ha对探测电磁波电阻率曲线测量深度进行表层深度校正；
第四步，其他采集参数的选择：
(1)采集开始/结束深度根据设计要求确定；
(2)采样间隔：
a、进行区域地层对比，追踪大地层电性界面使用：如5m、10m、20m。
b、进行矿层(油层、煤层、金属矿层)、裂缝带追踪与探测0.2m、0.5m、1m。
c、采样点数的设置 结束深度≤1000m取8点，或者
结束深度≤2000m取16点，或者
结束深度≤4000m取32点，或者
结束深度≤8000m取64点
d、采样周期数对应值为：128 64 32 16
(3)、通道增益设置：在保证目的层段信号采集基础上，选择较低增益可以抑制交流工业电源干扰；
(4)、低通滤波器设置：在工业电源干扰较强地区及浅层测量时使用自动跟踪滤波，在工业电源干扰弱地区及深层测量时可以不使用自动跟踪滤波；
第五步，资料的处理：
对野外采集资料经过质量验收之后，进行资料的处理：将采集记录文件调入解释处理程序中，输入起始深度和经过校正的采样间隔，进行适当数值处理，调节纵横向比例，获得采集点最终处理成果图。
本发明所述的大地电磁波电阻率测量仪器，包括接收电场强度信号的电场传感器，接收磁场强度信号的磁场传感器，两个分别与电场传感器和磁场传感器的输出端连接的前置放大与滤波器、与前置放大与滤波器的输出端连接的数据采集***和采集控制、数据存储及处理***；所述两个前置放大与滤波器的输出端分别连接到数据采集***的两个输入端，数据采集***的并行端口与采集控制、数据存储及处理***的打印机端口连接，或者数据采集***的并行端口与采集控制、数据存储及处理***的并行端口连接。
所述电场传感器(1)包括一对铜电极或者铅电极，或者不极化电极；磁场传感器采用多匝环形感应线圈，或者具有磁通负反馈的多匝环形感应线圈。
所述前置放大与滤波***包括放大器、多级陷波器、低通滤波器及输出级。
所述数据采集***用通用的高速数据采集器。
所述采集控制、数据存储与处理***是高笔记本电脑。
在所述采集控制、数据存储与处理***中收录有深度与传播频率之间的基本关系式，电磁波电阻率与深度的关系，深度系数以及参数选择所覆盖的各项界面参数。
本发明所述的方法与现有大地电磁测量法相比，探测深度相同时，分层精度大大提高。使得大地电磁测深由只能进行区域普查进入勘探阶段，甚至可以少打一部分钻井，大大提高勘探整体效益；测取同样的深度，对仪器灵敏度要求低得多。有利于仪器的制造；测量时间缩短，效率极大提高，外界干扰、仪器漂移、电极极化等影响减少；解释变得简单，在一维解释中，当深度确定以后，只需要读取深度点上的电场强度和磁场强度便可容易的获得视电阻率。可以方便地实现实时处理解释。由于不进行钻井便可以测得一条随深度变化的电阻率曲线，大大拓宽了利用地面地球物理勘探解决地质问题的能力。例如，油田水淹层监测，寻找裂缝带的分布，层间对比，煤层、金属矿层的追踪等等。大地电磁波电阻率测量方法和仪器成为一种运用于各个勘探阶段的有效方法和手段。从而提高总体勘探效益，同时为深部勘探提供一种有效手段。此外本发明也可以运用到有源的电磁测深中。
本发明所述的大地电磁波电阻率测量仪器经过大量实测资料的验证，达到下列指标：
深度误差：有井标定的地区，其误差不大于5％；
分层精度：高电阻率地层≥1米
低电阻率地层≥0.5米
探测深度40-8000m
300-8000m(在地面电磁干扰严重地区)
探测时间：每个测量点次，(测量深度≤4000m)工作时间≤1h；
仪器器重量：20kg±，对施工场地要求简单。
附图说明
图1是大地电磁波电阻率测量仪器示意图；
图2是前置放大滤波器原理框图；
图3数据解释处理界面图。
具体实施方式
本发明所述的大地电磁波电阻率测量方法和仪器技术思路，是利用实际数据直接标定频率与深度的关系。改变传统利用理论公式求取地层深度厚度的方法，使地层电阻率成为唯一的变量，只要获得每个采集深度点上的电场强度、磁场强度，便可以得到随深度连续变化的电阻率曲线。
电磁波传播理论指出，场强在地层中传播呈指数规律变化，或者说电场的振幅随地层深度呈指数规律衰减，传统电磁波方法，利用趋肤深度估计仪器的探测深度。本发明所指的探测深度取决于仪器的灵敏度和入射场源的强度。入射场源强度一定时，仪器灵敏度越高探测深度越深。例如定义振幅衰减到1/e3时(5％)为仪器的探测深度，则探测深度K′值是三倍趋肤深度K值。这个探测深度是现在的仪器灵敏度能够达到的。在一定的条件下，地层深度与电磁波的传播频率呈线性关系。
本发明所述方法的步骤是：
第一步，建立线性坐标观察***，预置初始采集参数：
进入一个勘探区，考虑到地层电阻率随深度的变化趋势，建立一个线性坐标观察***，在这个***中频率(F)与深度(H)的关系表示为：
F＝a+bK′/H (1)
式中：a-为表层系数，是一个较小的值，在第一步工作中可忽略
b-为区域性地层电阻率随深度的增长系数，其值设定为0.01
H-为深度，第一步测量时分别输入设计要求的起始测量深度值(H1)、终止测量深度值(H2)、采样间隔值(S)
K′-为电磁波在地层中传播定义的探测深度值。
在上式中H、H1、H2、S单位为米，F单位为赫兹(HZ)；(1)式可以进一步表示为：
F＝a+B/H
B定义为深度系数，在第一步中，B＝bk′
测量前，将由关系式确定的F参数预置到大地电磁波电阻率测量仪器各参数界面，分别输入设计要求的起始测量深度值(H1)、终止测量深度值(H2)、采样间隔值(S)，使仪器正常运行。
第二步，进行深度系数B的标定
为了使深度标定更加准确，需要进行深度系数B的标定。进入采集后，首先应当在探区或相邻地区，选择一口已知钻并。该已知钻井的条件应当符合以下条件：1、地层有代表性；2、有电阻率测井资料或钻井取芯等资料；3、地层倾角≤15度；4、地面无强工业电网等干扰。
利用第一步基本关系式确定的参数，在选择的已知钻井旁，测取一条电磁波电阻率曲线，用这条曲线与已知钻井的测井电阻率曲线进行对比，两条曲线是否在各特征点上是相似的。如果他们的相关性不好，按一定角度旋转电场传感器和磁场传感器的方向，重新进行采集，直到找到一条相关性最好的曲线为止。此时所测曲线纵坐标是电磁波电阻率、电场强度、磁场强度，横坐标是采集点序号。按下列步骤进行B参数标定：
将测得的随深度变化的电阻率曲线，截取一段与已知钻井中的电测电阻率曲线或者具有电性特征的其他地层分层资料进行比较，求得采样间隔S′(单位米)，该S′用下式表示：
S′＝(Hp2-Hp1)/(L2-L1) (2)
式中：Hp1-为已钻探井的第1个测井曲线特征点的深度(单位米)
Hp2-为已钻探井的第2个测井曲线特征点的深度(单位米)
L1-为新测探的与第1个测井特征点相似特征的电磁波电阻率曲线的采集点号数
L2-为新测探的与第2个测井特征点相似特征的电磁波电阻率曲线的采集点号数深度系数：
B＝S′/S (3)
B为实际确定的深度系数，在以后的工作中用B取代(1)中的bK′值进行频深关系变换。在没有深度标定井的新区，可以利用以前的与这个区块类似的其他区域的深度系数，对基本关系式进行标定后进行测量，也能获得较为满意的效果。仪器预置了一套经过实践检验的深度系数参数供使用。
第三步，进行表层深度系数Ha的校正：利用经过第二步参数校正的电磁波电阻率曲线与已知钻井电测井电阻率曲线进行全测量段比较，各特征深度点之间往往存在一个***误差，求取该误差值Ha：
Ha＝Hd-HC (4)
Hd--为已知钻井电测井曲线特征点深度
Hc--为探测电磁波电阻率曲线同一特征点深度
式中Ha、Hd、Hc单位为米。
利用表层深度系数Ha对探测电磁波电阻率曲线测量深度进行表层深度校正；
第四步，其他采集参数的选择：
1、采集开始/结束深度根据设计要求确定：
2、采样间隔：
(1)进行区域地层对比，追踪大地层电性界面使用：如5m、10m、20m。
(2)进行矿层(油层、煤层、金属矿层)、裂缝带追踪与探测0.2m、0.5m、1m。
(3)采样点数的设置 结束深度≤1000m取8点；
≤2000m取16点
≤4000m取32点
≤8000m取64点
(4)采样周期数对应值为：128 64 32 16
3、通道增益设置：在使用双传感器时，选择通道1(CH1)、通道2(CH2)，在使用多传感器时，选择通道1-8(CH1-8)；在保证目的层段信号采集基础上，选择较低增益可以抑制交流工业电源干扰。
4、低通滤波器设置：在工业电源干扰较强地区及浅层测量时使用自动跟踪滤波，在工业电源干扰弱地区及深层测量时可以不使用自动跟踪滤波。
第五步，资料的处理：
对野外采集资料经过质量验收之后，进行资料的处理：将采集记录文件调入解释处理程序中，输入起始深度和经过校正的采样间隔，进行适当数值处理，调节纵横向比例，便可以获得采集点最终处理成果图。见图3，其中横坐标为线性坐标，表示深度，深度比例可以根据要求进行调节；纵坐标：在上曲线中表示电磁波电阻率，具有线性和对数两种坐标供选择，比例可以调节，在中曲线中表示电场强度，比例可以调节，在下曲线中表示磁场强度，比例可以调节。
为完成本发明所述的方法，设计了大地电磁波电阻率测量仪器。
参见图1，本发明所述的大地电磁波电阻率测量仪器由五部分组成：接收电场强度信号的电场传感器1，接收磁场强度信号的磁场传感器2，两个前置放大与滤波器3，数据采集***4和采集控制、数据存储及处理***5。
所述电场传感器1的输出端连接到前置放大与滤波器3的输入端，磁场传感器2的输出端连接到另一个前置放大和滤波器3的输入端，两个前置放大与滤波器3的输出端分别连接到数据采集***的两个输入端CH1CH2，数据采集***4与采集控制及数据存储与处理***5的连接有两种方式：方式1，数据采集***4的并行端口(PARALLEL PORT)与采集控制及数据存储与处理***5的打印机端口连接；方式2，数据采集***4的并行端口(PARALLEL PORT)与采集控制和数据存储与处理***5的并行端口连接。数据采集***的数据传输给采集控制及数据存储与处理***，同时采集控制及数据存储与处理***向数据采集***输送控制、采集指令。
所述电场传感器1由一对铜电极或者铅电极组成，或者采用其他不极化电极。磁场传感器一般情况用多匝环形感应线圈，或者使用具有磁通负反馈的多匝环形感应线圈。两个传感器分别接收电场强度信号(Ex)和磁场强度信号(Hy)。
参见图2，前置放大与滤波器3由放大器、多级陷波器、低通滤波器及输出级组成，对由电场、磁场传感器接收的信号进行放大，并滤掉由工业电源及其谐波造成的干扰，防止混迭干扰的产生。
所述数据采集***4用通用的高速数据采集器。其主要性能：采样率：1MHz并可受程序控制调节，通道数：≥6，A/D转换：12或16位，可调增益：1-100，波形同步，程控抗混迭滤波器。数据采集***在程序控制下，对前置放大滤波器输入的模拟信号数字化，把按频率分组的数字信号输出到采集控制及数据存储与处理***5。
所述采集控制、数据存储及处理***5为高性能笔记本电脑，其内存储有根据方法中频率深度基本方程式F＝＝a+bK/H的程序，以及预置各项测量参数的程序，其主要任务是对数据采集***进行程序控制，按照要求完成预定的数据采集任务，并将成果数据纪录在硬盘中。电脑内还存储有数据处理的程序，根据大地电磁测深中的基本关系式：
ρ＝|Ex/Hy|2/ωμ，将采集的电场强度Ex和磁场强度Hy转换为电磁波电阻率。并可将最终成果打印输出。
仪器控制软件的设计充分体现技术思路，为了采集到不同深度的信息，不断改变采集频率，可以采集到相应深度点上的电场强度、磁场强度。最高可以按0.2米的深度等间隔采样；可以从技术指标范围内任意深度点对任意深度段采样；实时显示随深度线性变化的电场强度、磁场强度曲线；利用数据处理软件对采集资料进行现场处理，获得随深度H变化的电磁波电阻率曲线ρ。
Claims (10)
1、大地电磁波电阻率测量方法，其步骤是：
第一步，建立线性坐标观察***，预置初始采集参数：
进入一个勘探区，首先建立线性坐标观察***，在这个***中频率(F)与深度(H)的关系表示为：
F＝a+bK′/H (1)
式中：a-为表层系数，是一个较小的值，在第一步工作中可忽略
b-为区域性地层电阻率随深度的增长系数，其值设定为0.01
H-为深度(单位米)
K′-为电磁波在地层中传播定义的探测深度值
在上式中H、H1、H2、S单位为米，F单位为赫兹(HZ)
(1)式可以进一步表示为：
F＝a+B/H
B定义为深度系数，在第一步中，B＝bK′
测量时，分别输入设计要求的起始测量深度值(H1)、终止测量深度值(H2)、采样间隔值(S)，将由关系式确定的F参数预置到电阻率测量仪器各参数界面，使仪器正常运行。
第二步，进行深度系数B的标定：
进入采集后，首先应当在探区或相邻地区，选择一口已知钻井，利用第一步基本关系式确定的参数，在选择的已知钻井旁探测，测取一条电磁波电阻率曲线，用这条曲线与已知钻井的测井电阻率曲线进行对比，如此反复，直到找到一条相关性最好的曲线为止，按下列步骤进行B参数标定：
选择已知钻井的条件是：(1)地层有代表性；(2)有电阻率测井资料或钻井取芯等资料；(3)地层倾角≤15度：(4)地面无强工业电网等干扰。
将测得的随深度变化的电阻率曲线，截取一段与已知钻井中的电测电阻率曲线或者具有电性特征的其他地层分层资料进行比较，求得采样间隔S′(单位米)，该S′用下式表示：
S′＝(Hp2-Hp1)/(L2-L1) (2)
式中：Hp1-为已钻探井的第1个测井曲线特征点的深度，单位米
Hp2-为已钻探井的第2个测井曲线特征点的深度，单位米
L1-为新测探的与第1个测井特征点相似特征的电磁波电阻率曲线的采集点号数
L2-为新测探的与第2个测井特征点相似特征的电磁波电阻率曲线的采集点号数
B＝S′/S (3)
B为实际确定的深度系数，在以后的工作中用B取代(1)式中的bK值进行频率深度关系变换；在没有深度标定井的新区，可以利用以前的与这个区块类似的其他区域的深度系数，对基本关系式进行标定后进行测量。
第三步，进行表层深度系数Ha的校正：
利用经过第二步参数校正的曲线与已知钻井电测井电阻率曲线进行全测量段比较，各特征深度点之间存在的***误差值用下式表示：
Ha＝Hd-HC (4)
Hd--为已知钻井电测井曲线特征点深度
Hc--为探测电磁波电阻率曲线同一特征点深度
式中Ha、Hd、Hc单位为米；
利用Ha对探测电磁波电阻率曲线测量深度进行表层深度校正；
第四步，其他采集参数的选择：
(1)采集开始/结束深度根据设计要求确定；
(2)采样间隔：
a、进行区域地层对比，追踪大地层电性界面使用：如5m、10m、20m。
b、进行矿层(油层、煤层、金属矿层)、裂缝带追踪与探测0.2m、0.5m、1m。
c、采样点数的设置 结束深度≤1000m取8点，或者
结束深度≤2000m取16点，或者
结束深度≤4000m取32点，或者
结束深度≤8000m取64点
d、采样周期数对应值为：128 64 32 16
(3)、通道增益设置：在保证目的层段信号采集基础上，选择较低增益可以抑制交流工业电源干扰；
(4)、低通滤波器设置：在工业电源干扰较强地区及浅层测量时使用自动跟踪滤波，在工业电源干扰弱地区及深层测量时可以不使用自动跟踪滤波；
第五步，资料的处理：
对野外采集资料经过质量验收之后，进行资料的处理：将采集记录文件调入解释处理程序中，输入起始深度和经过校正的采样间隔，进行适当数值处理，调节纵横向比例，获得采集点最终处理成果图。
2、大地电磁波电阻率测量仪器，其特征是包括接收电场强度信号的电场传感器(1)，接收磁场强度信号的磁场传感器(2)，两个分别与电场传感器(1)和磁场传感器(2)的输出端连接的前置放大与滤波器(3)、与前置放大与滤波器(3)的输出端连接的数据采集***(4)和采集控制、数据存储及处理***(5)；所述两个前置放大与滤波器(3)的输出端分别连接到数据采集***的两个输入端(CH1 CH2)，数据采集***(4)的并行端口与采集控制、数据存储及处理***(5)的打印机端口连接，或者数据采集***(4)的并行端口与采集控制、数据存储及处理***(5)的并行端口连接。
3、根据权利要求2所述的测量仪器，其特征是所述电场传感器(1)包括一对铜电极或者铅电极，或者不极化电极；磁场传感器采用多匝环形感应线圈，或者具有磁通负反馈的多匝环形感应线圈。
4、根据权利要求2或3所述的测量仪器，其特征是所述前置放大与滤波***包括放大器、多级陷波器、低通滤波器及输出级。
5、根据权利要求2或3所述的测量仪器，其特征是所述数据采集***用通用的高速数据采集器。
6、根据权利要求4所述的测量仪器，其特征是所述数据采集***用通用的高速数据采集器。
7、根据权利要求5所述的测量仪器，其特征是所述采集控制、数据存储与处理***是高笔记本电脑。
8、根据权利要求6所述的测量仪器，其特征是所述采集控制、数据存储与处理***是高笔记本电脑。
9、根据权利要求7所述的测量仪器，其特征是在所述采集控制、数据存储与处理***中收录有深度(H)与传播频率(F)之间的基本关系式，电磁波电阻率(ρ)与深度(H)的关系，深度系数(B)以及参数选择所覆盖的各项界面参数。
10、根据权利要求8所述的测量仪器，其特征是在所述采集控制、数据存储与处理***中收录有深度(H)与传播频率(F)之间的基本关系式，电磁波电阻率(ρ)与深度(H)的关系，深度系数(B)以及参数选择所覆盖的各项界面参数。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2003101130358A CN1252492C (zh) | 2003-12-25 | 2003-12-25 | 大地电磁波电阻率测量方法及其仪器 |
| EP04802419A EP1703304A1 (en) | 2003-12-25 | 2004-12-02 | The method and apparatus for measuring resistivity of earth by elelctromagnetic waves |
| PCT/CN2004/001403 WO2005083468A1 (fr) | 2003-12-25 | 2004-12-02 | Procede et appareil destines a mesurer la resistivite de la terre par ondes electromagnetiques |
| CA002551588A CA2551588A1 (en) | 2003-12-25 | 2004-12-02 | Method and apparatus for measuring the resistivity of electromagnetic waves of the earth |
| US10/583,807 US7519474B2 (en) | 2003-12-25 | 2004-12-02 | Method and apparatus for measuring the resistivity of electromagnetic waves of the earth |
| JP2006545885A JP2007517201A (ja) | 2003-12-25 | 2004-12-02 | 大地電磁波抵抗率測定方法及びその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2003101130358A CN1252492C (zh) | 2003-12-25 | 2003-12-25 | 大地电磁波电阻率测量方法及其仪器 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1632617A CN1632617A (zh) | 2005-06-29 |
| CN1252492C true CN1252492C (zh) | 2006-04-19 |
Family
ID=34843400
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2003101130358A Expired - Fee Related CN1252492C (zh) | 2003-12-25 | 2003-12-25 | 大地电磁波电阻率测量方法及其仪器 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US7519474B2 (zh) |
| EP (1) | EP1703304A1 (zh) |
| JP (1) | JP2007517201A (zh) |
| CN (1) | CN1252492C (zh) |
| CA (1) | CA2551588A1 (zh) |
| WO (1) | WO2005083468A1 (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101482013B (zh) * | 2009-02-23 | 2010-05-12 | 中国石化集团胜利石油管理局钻井工艺研究院 | 一种随钻井眼补偿电磁波电阻率测量装置 |
| CN102498418A (zh) * | 2009-09-10 | 2012-06-13 | 布西鲁斯欧洲有限公司 | 用于原始矿床地电勘探的传感器设备和方法 |
Families Citing this family (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1252492C (zh) | 2003-12-25 | 2006-04-19 | 周仁安 | 大地电磁波电阻率测量方法及其仪器 |
| CN100480734C (zh) * | 2007-03-08 | 2009-04-22 | 刘俊昌 | 一种高分辨率去静态频率域大地电磁法 |
| CN101050700B (zh) * | 2007-05-09 | 2010-11-17 | 中国海洋石油总公司 | 一种用高导磁材料的电磁波测井仪校验方法 |
| JP5662804B2 (ja) | 2007-12-18 | 2015-02-04 | シュルンベルジェ ホールディングス リミテッドＳｃｈｌｎｍｂｅｒｇｅｒ Ｈｏｌｄｉｎｇｓ Ｌｉｍｉｔｅｄ | 表面電磁探査を改善するシステム及び方法 |
| CN101556340B (zh) * | 2008-04-10 | 2011-08-03 | 中国石油集团东方地球物理勘探有限责任公司 | 三维小面元大地电磁连续阵列数据采集方法 |
| US9250352B2 (en) * | 2008-04-17 | 2016-02-02 | Richard H. Hardman | Methods for producing a log of material properties |
| US8242781B2 (en) * | 2008-08-20 | 2012-08-14 | Lockheed Martin Corporation | System and method for determining sub surface geological features at an existing oil well site |
| US8521484B2 (en) * | 2010-06-02 | 2013-08-27 | Livermore Software Technology Corp. | Curve matching for parameter identification |
| BR112013032514A2 (pt) | 2011-06-21 | 2020-12-08 | Groundmetrics, Inc. | Sistema, método e aparelho para medição ou geração de campo elétrico em fundo de poço |
| CN102419455B (zh) * | 2011-08-23 | 2013-08-14 | 安徽理工大学 | 井间并行电阻率ct测试方法 |
| CN102998709B (zh) * | 2012-10-31 | 2015-09-16 | 江西理工大学 | 硐室***盲炮检测地面接收*** |
| CN103015995B (zh) * | 2012-12-14 | 2015-06-24 | 中国电子科技集团公司第二十二研究所 | 通过带倾斜线圈的电阻率设备的地层倾角测量方法及装置 |
| ES2499915B1 (es) * | 2013-03-27 | 2015-04-29 | Fundación Attico | Medidor de espesor de estratos de materiales |
| US10444404B2 (en) * | 2013-07-26 | 2019-10-15 | Halliburton Energy Services Inc. | System, method and computer-program product for in-situ calibration of a wellbore resistivity logging tool |
| US10270314B2 (en) * | 2013-10-02 | 2019-04-23 | Kawasaki Jukogyo Kabushiki Kaisha | Cooling for drive motor and transmission of electric vehicle |
| US9293843B2 (en) * | 2014-03-21 | 2016-03-22 | Yi Lu | Non-polarized geophysical electrode |
| CN104375191B (zh) * | 2014-12-08 | 2017-09-26 | 中国石油天然气集团公司 | 一种大地电磁数据采集***及其采集方法 |
| CA2968501C (en) * | 2015-01-07 | 2022-05-24 | Halliburton Energy Services, Inc. | Functional earth model parameterization for resistivity inversion |
| CN105607130A (zh) * | 2015-09-08 | 2016-05-25 | 浙江华东工程安全技术有限公司 | 一种用钻孔高密度电法检测孔周岩溶或空洞的方法 |
| CN107085240B (zh) * | 2017-03-30 | 2020-01-10 | 湖南科技大学 | 一种边坡磁流体探测***及方法 |
| CN108873083B (zh) * | 2018-05-06 | 2019-11-15 | 东华理工大学 | 一种人工场源频率域电磁视电阻率测量方法 |
| US11435497B2 (en) * | 2018-09-22 | 2022-09-06 | Halliburton Energy Services, Inc. | Three dimensional visualization from point-by-point one dimensional inversion with bed azimuth |
| CN111983674A (zh) * | 2019-05-22 | 2020-11-24 | 中国石油天然气股份有限公司 | 地层倾斜区的地层真实厚度确定方法和装置 |
| CN112539056B (zh) | 2019-09-05 | 2022-11-25 | 中国石油化工股份有限公司 | 用于多维井下的成像特征提取方法及成像装置 |
| CN110703342B (zh) * | 2019-11-20 | 2020-06-30 | 中国科学院地质与地球物理研究所 | 一种全覆盖抗干扰人工源电磁探测方法 |
| CN111208572A (zh) * | 2020-02-28 | 2020-05-29 | 重庆慧勘工业技术研究院有限公司 | 一种基于矿用瞬变电磁仪的多剖面动态成像方法及装置 |
| CN111736225A (zh) * | 2020-06-11 | 2020-10-02 | 中化地质矿山总局地质研究院 | 一种中-低温热液填充型萤石矿勘探方法及勘探装置 |
| CN112114367B (zh) * | 2020-09-25 | 2023-05-26 | 中国人民解放军海军工程大学 | 一种利用主动电场探测海域中电学分界面的方法 |
| CN112363245A (zh) * | 2020-10-26 | 2021-02-12 | 王昌益 | 一种地下水物探测量方法及*** |
| CN112946762B (zh) * | 2021-03-17 | 2022-03-25 | 中南大学 | 基于水平电场和磁场获取电阻率的频率域电磁勘探方法 |
| CN113238290B (zh) * | 2021-05-11 | 2023-10-27 | 中煤科工集团西安研究院有限公司 | 一种煤矸石场内部高温区双参数联合无损探测方法 |
| CN113325486A (zh) * | 2021-06-03 | 2021-08-31 | 青海省第一地质勘查院 | 覆盖区下的构造蚀变岩型矿物勘查方法、***及装置 |
| CN114089429B (zh) * | 2021-10-27 | 2024-02-09 | 中煤科工集团西安研究院有限公司 | 一种利用钻孔雷达构建煤矿透明工作面的方法 |
| CN114137036B (zh) * | 2021-11-29 | 2024-03-29 | 安徽理工大学 | 一种基于邻源电位电阻率的注浆范围快速检测方法 |
| CN114779352B (zh) * | 2022-05-31 | 2022-11-22 | 成都理工大学 | 一种基于空-地-井多维电磁法的地质体勘测*** |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4107598A (en) * | 1976-02-02 | 1978-08-15 | Texaco Inc. | Electromagnetic wave logging system for determining resistivity and dielectric constant of earth formations |
| US4107597A (en) * | 1976-12-13 | 1978-08-15 | Texaco Inc. | Electromagnetic wave propagation well logging utilizing multiple phase shift measurement |
| US4663593A (en) * | 1985-03-15 | 1987-05-05 | Standard Oil Company | Method of electromagnetic exploration |
| US4899112A (en) * | 1987-10-30 | 1990-02-06 | Schlumberger Technology Corporation | Well logging apparatus and method for determining formation resistivity at a shallow and a deep depth |
| WO1992022833A1 (en) * | 1991-06-14 | 1992-12-23 | Baroid Technology, Inc. | Method and apparatus for detecting boundary stratum |
| US5469062A (en) * | 1994-03-11 | 1995-11-21 | Baker Hughes, Inc. | Multiple depths and frequencies for simultaneous inversion of electromagnetic borehole measurements |
| US5777478A (en) * | 1996-08-27 | 1998-07-07 | John R. Jackson | Passive geophysical prospecting apparatus and method based upon detection of discontinuities associated with extremely low frequency electromagnetic fields |
| DE19652835C1 (de) | 1996-12-18 | 1998-03-26 | Bauer Spezialtiefbau | Schlitzwandfräse |
| US6114972A (en) * | 1998-01-20 | 2000-09-05 | Halliburton Energy Services, Inc. | Electromagnetic resistivity tool and method for use of same |
| US6476609B1 (en) * | 1999-01-28 | 2002-11-05 | Dresser Industries, Inc. | Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone |
| US6163155A (en) | 1999-01-28 | 2000-12-19 | Dresser Industries, Inc. | Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations |
| CN1185508C (zh) * | 2000-07-07 | 2005-01-19 | 石油大学(北京) | 无趋肤效应电磁波测井方法 |
| US6892137B2 (en) * | 2003-04-29 | 2005-05-10 | Pathfinder Energy Services, Inc. | Adjustment for frequency dispersion effects in electromagnetic logging data |
| CN1252492C (zh) | 2003-12-25 | 2006-04-19 | 周仁安 | 大地电磁波电阻率测量方法及其仪器 |
| GB2411728B (en) * | 2004-03-01 | 2008-05-28 | Pathfinder Energy Services Inc | Azimuthally focussed electromagnetic measurement tool |
| US7385400B2 (en) * | 2004-03-01 | 2008-06-10 | Pathfinder Energy Services, Inc. | Azimuthally sensitive receiver array for an electromagnetic measurement tool |
| GB2413851B (en) * | 2004-05-06 | 2006-08-09 | Ohm Ltd | Electromagnetic surveying for hydrocarbon reservoirs |
-
2003
- 2003-12-25 CN CNB2003101130358A patent/CN1252492C/zh not_active Expired - Fee Related
-
2004
- 2004-12-02 WO PCT/CN2004/001403 patent/WO2005083468A1/zh active Application Filing
- 2004-12-02 CA CA002551588A patent/CA2551588A1/en not_active Abandoned
- 2004-12-02 EP EP04802419A patent/EP1703304A1/en not_active Withdrawn
- 2004-12-02 US US10/583,807 patent/US7519474B2/en active Active
- 2004-12-02 JP JP2006545885A patent/JP2007517201A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101482013B (zh) * | 2009-02-23 | 2010-05-12 | 中国石化集团胜利石油管理局钻井工艺研究院 | 一种随钻井眼补偿电磁波电阻率测量装置 |
| CN102498418A (zh) * | 2009-09-10 | 2012-06-13 | 布西鲁斯欧洲有限公司 | 用于原始矿床地电勘探的传感器设备和方法 |
| CN102498418B (zh) * | 2009-09-10 | 2016-02-03 | 卡特彼勒全球矿场欧洲有限公司 | 用于原始矿床地电勘探的传感器设备和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2007517201A (ja) | 2007-06-28 |
| US20080228401A1 (en) | 2008-09-18 |
| CA2551588A1 (en) | 2005-09-09 |
| CN1632617A (zh) | 2005-06-29 |
| WO2005083468A1 (fr) | 2005-09-09 |
| EP1703304A1 (en) | 2006-09-20 |
| US7519474B2 (en) | 2009-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1252492C (zh) | 大地电磁波电阻率测量方法及其仪器 | |
| CN1276377C (zh) | 使用虚拟源阵列对地下岩层成像的方法 | |
| CN105221133B (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| CN107578343B (zh) | 一种碳酸盐岩缝洞型储层含水饱和度计算方法及其装置 | |
| CN1017916B (zh) | 测井方法及装置 | |
| CN100351650C (zh) | 一种利用叠前地震波形反演构建虚拟井数据的方法 | |
| CN201334906Y (zh) | 三相位补偿电阻率测井仪器 | |
| CN1950721A (zh) | 关于碳氢化合物储藏的电磁勘探 | |
| CN1740517A (zh) | 高频阵列电磁波传播测井的方法及其测井仪 | |
| CN109782360A (zh) | 一种低阻覆盖区深部砂体探测方法 | |
| CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
| Cui et al. | Application of the Ground Penetrating Radar ARMA power spectrum estimation method to detect moisture content and compactness values in sandy loam | |
| CN206378448U (zh) | 基于综合物探方法的污染土检测测线布置结构 | |
| WO2020106287A1 (en) | Enhanced anisotropy analysis with multi-component dipole sonic data | |
| CN106802432A (zh) | 土层结构的探查方法和装置 | |
| CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
| CN1209499A (zh) | 阵列电极系成象测井方法 | |
| CN111927443A (zh) | 一种基于电场测量地层真电阻率的测井***与方法 | |
| CN110133727A (zh) | 超深油气储层高精度电磁频谱探测方法 | |
| CN209264613U (zh) | 一种基于电阻率ct成像的垃圾填埋场污染物扩散监测装置 | |
| CN108594301B (zh) | 一种具有差异特征的地震数据融合的方法及处理终端 | |
| CN1245637C (zh) | 管波探测法 | |
| CN113216945A (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
| CN108761540A (zh) | 一种频率域天然电场三维勘探方法 | |
| CN1877366A (zh) | 重磁延拓回返垂直导数目标优化处理技术 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20060419 |