CN107152276A - 用于随钻电磁波电阻率测量的信号处理方法与系统 - Google Patents
用于随钻电磁波电阻率测量的信号处理方法与系统 Download PDFInfo
- Publication number
- CN107152276A CN107152276A CN201610118358.3A CN201610118358A CN107152276A CN 107152276 A CN107152276 A CN 107152276A CN 201610118358 A CN201610118358 A CN 201610118358A CN 107152276 A CN107152276 A CN 107152276A
- Authority
- CN
- China
- Prior art keywords
- mrow
- signal
- frequency
- msub
- component
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 21
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 11
- 230000003321 amplification Effects 0.000 claims description 10
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 10
- 239000013078 crystal Substances 0.000 claims description 7
- 238000005265 energy consumption Methods 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 description 16
- 238000006243 chemical reaction Methods 0.000 description 8
- 241001074085 Scophthalmus aquosus Species 0.000 description 6
- 238000005562 fading Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 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
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/13—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling by electromagnetic energy, e.g. radio frequency
Landscapes
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种用于随钻电磁波电阻率测量的信号处理方法与系统，方法包括：将预设频率的发射信号发射到待测地层；获取经待测地层反射回来的第一电磁波信号和第二电磁波信号；对两电磁波信号进行放大，得到第一放大信号和第二放大信号；利用数字开关模块来分离两放大信号中的高、低频分量，数字开关模块由数字开关控制信号控制；过滤掉两放大信号中的高频分量；计算两放大信号的低频分量的相位差和幅度比；根据相位差和幅度比，得到待测地层的电阻率。本发明采用的数字开关控制模块与滤波模块所需的计算量小，使得系统的功耗低，电路简单，可靠性高。
Description
技术领域
本发明涉及油气开发与勘探技术领域，尤其涉及一种用于随钻电磁波电阻率测量的信号处理方法，还涉及一种用于随钻电磁波电阻率测量的信号处理系统。
背景技术
地层电阻率是反映地层特性的主要参数，随钻地层电阻率的测量是随钻测井的关键组成部分。通过随钻测量地层的电阻率判断地层特性变化，再结合其它地层资料，经过综合分析，可以及时地判断钻头是否钻出储层，从而及时调整井眼轨迹，确保井眼最大程度穿行于储层中有利于油气开采的最佳位置。目前使用的随钻电阻率测量工具以随钻电磁波电阻率测量应用最为广泛。由于电磁波在穿越地层时产生幅度衰减和相位偏移，并且产生幅度衰减和相位偏移由地层的电阻率和介电常数的特性所决定，电磁波在穿越不同的地层介质时产生的幅度衰减和相位偏移不同。不同的频率的电磁波对地层介质的响应特征也不同，当10MHz以上频率的电磁波，其产生的幅度衰减和相位偏移主要与地层的介电常数相关；而10MHz以下的电磁波，其产生的电磁波幅度衰减和相位偏移主要与地层的电阻率相关。地层电阻率对于地质导向和地层评价是重要地质参数。用于随钻测量电阻率的方法也很多，随钻电磁波测量方式可以应用在油基钻井液，这是随钻侧向电阻率测量方式所不能实现的。
现有的随钻电磁波电阻率测量方法的不足在于：现有的随钻电磁波电阻率测量信号发射和接收电路为了实现一个线圈及其匹配电路能够在两个不同的频率(2000KHz和400KHz)下实现谐振，结构较为复杂，功耗较大。由于2000KHz和400KHz的信号直接采集需要复杂的电路系统。为了便于采样，现有技术都是使用混频器来实现高频向低频的转换，将2000KHz或400KHz转换为几KHz的低频信号。但是混频器的功耗较大，而且需要较高的供电电压，与前后的模拟放大和滤波电路不兼容。在实际操作过程中，通常使用高温FPGA来实现数字混频。因为高温FPGA属于非常复杂的芯片，其功耗非常大、芯片引脚非常多，对电路板设计要求极高。这些都会增加电路系统的复杂程度，降低电路系统的可靠性。
因此，期望提供一种简单有效的随钻电磁波电阻率仪器的信号处理方法及系统。
发明内容
本发明所要解决的技术问题是：现有的随钻电磁波电阻率测量方法涉及数字混频，计算复杂，需要功耗大、芯片引脚多的高级处理器来实现，从而在增加了电路的复杂程度的同时，降低了电路的可靠性。
为了解决上述技术问题，本发明提供了一种简单有效的随钻电磁波电阻率仪器的信号处理方法及系统。
根据本发明的一个方面，提供了一种用于随钻电磁波电阻率测量的信号处理方法，其包括：
将预设频率的发射信号发射到待测地层；
获取经待测地层反射回来的第一电磁波信号和第二电磁波信号；
对所述第一电磁波信号和所述第二电磁波信号进行放大，得到第一放大信号和第二放大信号；
利用数字开关模块来分离所述第一放大信号和所述第二放大信号中的高、低频分量，得到对应所述第一放大信号的第一高频分量和第一低频分量，以及对应所述第二放大信号的第二高频分量和第二低频分量，所述数字开关模块由数字开关控制信号控制；
过滤掉所述第一高频分量和所述第二高频分量，以提取出所述第一低频分量和所述第二低频分量；
计算所述第一低频分量和所述第二低频分量的相位差和幅度比；
根据所述相位差和幅度比，得到所述待测地层的电阻率。
优选的是，所述数字开关模块的输入为第一/第二放大信号，所述数字开关控制信号F(t)满足：
其中：
并且，t表示时间，f0表示第一频率，f1表示第二频率，所述第一频率大于所述第二频率，并且第一/第二放大信号中频率大于所述第一频率的分量为所述高频分量，频率小于或者等于所述第一频率的分量为所述低频分量。
优选的是，将预设频率的发射信号发射到待测地层，包括：
生成预设频率的发射信号；
对所述发射信号进行放大；
对放大后的发射信号进行阻抗匹配；
由发射线圈将经阻抗匹配的信号发射到所述待测地层。
优选的是，所述预设频率为2000KHz或者400KHz。
根据本发明的另一个方面，提供了一种用于随钻电磁波电阻率测量的信号处理系统，其包括：
发射器，设置为将预设频率的发射信号发射到待测地层；
第一接收天线，设置为获取经待测地层反射回来的第一电磁波信号；
第二接收天线，设置为获取经待测地层反射回来的第二电磁波信号；
接收放大器，设置为对所述第一电磁波信号和所述第二电磁波信号进行放大，得到第一放大信号和第二放大信号；
数字开关模块，设置为分离所述第一放大信号和所述第二放大信号中的高、低频分量，得到对应所述第一放大信号的第一高频分量和第一低频分量，以及对应所述第二放大信号的第二高频分量和第二低频分量，所述数字开关由数字开关控制信号控制；
数字控制振荡器，设置为产生所述数字开关控制信号；
滤波模块，设置为过滤掉所述第一高频分量和所述第二高频分量，以提取出所述第一低频分量和所述第二低频分量；
计算模块，设置为计算所述第一低频分量和所述第二低频分量的相位差和幅度比；并根据所述相位差和幅度比，得到所述待测地层的电阻率。
优选的是，所述数字开关模块的输入为第一/第二放大信号，所述数字开关控制信号F(t)满足：
其中：
并且，t表示时间，f0表示第一频率，f1表示第二频率，所述第一频率大于所述第二频率，并且第一/第二放大信号中频率大于所述第一频率的分量为所述高频分量，频率小于或者等于所述第一频率的分量为所述低频分量。
优选的是，所述发射器包括：
晶体振荡器，设置为生成预设频率的发射信号；
发射放大器，设置为对所述发射信号进行放大；
阻抗匹配电路，设置为对放大后的发射信号进行阻抗匹配；
发射线圈，设置为将经阻抗匹配的信号发射到所述待测地层。
优选的是，所述预设频率为2000KHz或者400KHz。
与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
本发明采用数字开关控制通断的方式实现信号中高频分量和低频分量的有效分离，结合滤波器，可实现低频分量的有效提取(即高频分量向低频分量的转换)。与现有技术中采用混频器来实现高频分量向低频分量的转换的技术方案相比，本发明采用的数字开关控制模块与滤波模块所需的计算量小，因此采用普通的处理器即可满足数字开关控制模块和滤波模块的计算需求。可以看出，采用本本发明的系统的功耗低，电路简单，可靠性高。
本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中：
图1示出了本发明实施例用于随钻电磁波电阻率测量的信号处理方法的流程示意图；
图2示出了本发明实施例中将预设频率的发射信号发射到待测地层的方法的流程示意图；
图3示出了本发明实施例用于随钻电磁波电阻率测量的信号处理系统的发射部分的结构示意图；
图4示出了本发明实施例用于随钻电磁波电阻率测量的信号处理系统的接收部分的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
本发明所要解决的技术问题是：现有的随钻电磁波电阻率测量方法涉及数字混频，计算复杂，需要功耗大、芯片引脚多的高级处理器来实现，从而在增加了电路的复杂程度的同时，降低了电路的可靠性。为解决现有技术中，本发明实施例提供了一种用于随钻电磁波电阻率测量的信号处理方法。
图1示出了本发明实施例用于随钻电磁波电阻率测量的信号处理方法的流程示意图。参照图1，本实施例用于随钻电磁波电阻率测量的信号处理方法主要包括步骤10至步骤70。
在步骤10中，将预设频率的发射信号发射到待测地层。
参照图2，本步骤进一步包括步骤11至步骤14。
在步骤11中，生成预设频率的发射信号。
具体地，晶体振荡器经过分频生成预设频率的发射信号。特别地，预设频率优选为2000KHz或400KHz，此种频率的电磁波的幅度衰减和相位偏移主要与地层的电阻率相关，从而有助于进行地质导向和地层评价。
在步骤12中，对发射信号进行放大。发射信号经过功率放大器进行放大。
在步骤13中，对放大后的发射信号进行阻抗匹配。将放大后的发射信号输入至阻抗匹配电路，以使阻抗匹配电路对该信号进行阻抗匹配，从而使信号更有效地耦合至发射线圈。
在步骤14中，由发射线圈将经阻抗匹配的信号发射到待测地层。
在步骤20中，获取经待测地层反射回来的第一电磁波信号和第二电磁波信号。
具体地，通过两个接收天线与匹配电路接收到携带地层电阻率信息的微弱电磁波信号，第一接收天线接收的信号为第一电磁波信号，第二接收天线接收的信号为第二电磁波信号。
在步骤30中，对第一电磁波信号和第二电磁波信号进行放大，得到第一放大信号和第二放大信号。由前级低噪声放大器对接收天线接收到的电磁波信号进行放大，得到对应第一接收天线的第一放大信号和对应第二接收天线的第二放大信号。
在步骤40中，利用数字开关模块来分离第一放大信号和第二放大信号中的高、低频分量，得到对应第一放大信号的第一高频分量和第一低频分量，以及对应第二放大信号的第二高频分量和第二低频分量；数字开关模块由数字开关控制信号控制。
在本步骤中，高频分量和低频分量是相对于某个预设的频率而言。一般来讲，信号中大于预设的频率(例如2000KHz)的分量统称为高频分量，信号中小于或者等于预设的频率的分量统称为低频分量。可以通过合理地设计数字开关模块的函数模型，来使某一信号中的高频分量和低频分量分离开来。
在步骤50中，过滤掉第一高频分量和第二高频分量，以提取出第一低频分量和第二低频分量。利用带通滤波器将分离开来的高频分量过滤掉，只剩下低频分量，通过此种方式完成信号中低频分量的有效提取。
在步骤60中，计算第一低频分量和第二低频分量的相位差和幅度比。具体地，对得到的两路低频信号(即第一低频信号和第二低频信号)进行ADC采样后，计算两者的相位差和幅度比。
在步骤70中，根据相位差和幅度比，得到待测地层的电阻率。具体地，利用计算得到的相位差和幅度比，再对照地层电阻率图板，反演得到待测地层的地层电阻率。由于利用两组信号的相位差的幅度比，结合地层电阻率图板来反演得到地层电阻率是本领域技术人员获得地层电阻率所采用的惯常技术手段，因此在本文中不再对本步骤涉及的方法进行展开说明。
本实施例所述的信号处理方法，采用数字开关控制通断的方式实现信号中高频分量和低频分量的有效分离，结合滤波器，可实现低频分量的有效提取(即高频分量向低频分量的转换)。与现有技术中采用混频器来实现高频分量向低频分量的转换的技术方案相比，本实施例的数字开关控制模块与滤波模块所需的计算量小，因此采用普通的处理器即可满足数字开关控制模块和滤波模块的计算需求。可以看出，采用本实施例所述方法的系统的功耗低，电路简单，可靠性高。
此外，本实施例还保留幅度比和相位差信息的随钻电磁波电阻率信号处理方法与电路系统，这样的模拟信号处理方式可以有效降低后续模拟到数字转换，以及数字信号处理的性能要求。这样可以提高采样精度，进一步提高地层电阻率的测量精度，降低随钻电磁波电阻率测量仪器的功耗，延长工作时间。
由于本发明使用了数字开关模块对高频电磁波信号进行处理，因此需要介绍这种开关通断处理实现高频信号到低频信号转换的原理。在此之前，首先介绍传统混频器的工作原理，以作用本实施例的对比。
传统的采用混频器的方式实现高频到低频转换，首先需要使用一个锁相频率综合器(PLL,Phase Locked Loop)或者直接数字频率综合器(DDS,Direct DigitalSynthesizer)来生成一个与发射信号频率具有较小频率偏差的高频模拟信号。将该信号与接收的信号通过混频器进行混频，混频后生成一个低频信号和一个高频信号。
具体计算方法如下：
假设发射的信号频率为f0，f0通常取2000KHz或者400KHz。接收的信号与发射信号频率一致，其函数表达式表示为F1(t)＝A1sin(2πf0t+ω1)，其中A1是该信号的幅度，ω1是该信号的初始相位。由锁相频率综合器PLL或者直接数字频率综合器DDS生成的信号频率为f0+f1，f1通常取2KHz，其函数表达式表示为F2(t)＝A2sin(2π(f0+f1)t+ω2)，其中A2是该信号的幅度，ω2是该信号的初始相位。
接收的信号与锁相频率综合器PLL或者直接数字频率综合器DDS生成的信号通过混频器进行混频得到一个新的信号F3(t)，其函数表达式如下：
F3(t)＝F1(t)*F2(t)
＝A1sin(2πf0t+ω1)*A2sin[2π(f0+f1)t+ω2]
＝-1/2A1A2{cos[2πf0t+ω1+2π(f0+f1)t+ω2]-cos[2πf0t+ω1-2π(f0+f1)t-ω2]}
＝-1/2A1A2{cos[2π(2f0+f1)t+ω1+ω2]-cos(2πf1t+ω1-ω2)}
＝-1/2A1A2cos[2π(2f0+f1)t+ω1+ω2]+1/2A1A2cos(2πf1t+ω1-ω2)
从生成的信号F3(t)的表达式可以看出，该信号分成两个部分，前半部分信号的频率为2f0+f1，后半部分的信号的频率为f1，再将该信号输入一个低通滤波器，可以将前半部分滤出，只留下后半部分。该后半部分仍包含接收信号的幅度和初始相位信息。因此通过采样两路该信号可以进行一步计算出这两路信号的幅度比和相位差。
上述传统方法需要采用锁相频率综合器PLL或者直接数字频率综合器DDS等模拟电路，生成的频率为f0+f1的模拟信号，而且需要使用混频器进行混频得到一个新的信号F3(t)，该信号频率为远低于f0的f1。
本发明实施例采用数字开关控制信号对接收的信号进行开关控制，数字开关控制信号为频率为f0+f1的方波信号，该信号使用数字控制振荡器(NCO,Numerically Controlled Oscillator)产生，该电路为数字电路，比传统方法的模拟电路(例如PLL或DDS)简单。由数字开关控制信号信号的函数表达式为：
其中：
并且，t表示时间，f0表示第一频率，f1表示第二频率，第一频率大于第二频率，并且第一/第二放大信号中频率大于第一频率的分量为高频分量，频率小于或者等于第一频率的分量为低频分量。
本发明使用该数字开关控制信号控制一个开关电路(即数字开关模块)，开关电路的输入为接收的信号(即第一放大信号或第二放大信号)。该接收信号的表达式与传统方法的接收信号表达式一致：F1(t)＝A1sin(2πf0t+ω1)。该接收信号作为输入，数字开关控制信号作为控制信号，连接到开关电路的控制端口，由此，
开关电路的输出信号表达式如下：
F4(t)＝F1(t)*F3(t)
其中F3(t)经过复立叶变换展开表达式如下：
从上述表达式可以看出，该信号就是由一个直流信号与频率为f0+f1的1、3、5、…各次谐波之和。
因此有：
由数字控制开关信号的函数表达式，可知上述表达式仅有第二项：
经过三角函数的积化和差运算得到：
从上述表达式可以看出：表达式的后半部分的频率为f1(低频分量)，其它的组成频率都超过f0(高频分量)，因此经过后续连接的带通滤波器可以将这些频率都超过f0的高频分量全部滤除。经过滤波，上述表达式只剩下：
可以看出，放大的电磁波信号(例如，第一放大信号、第二放大信号)被频率为f0+f1(例如2002KHz或402KHz)的数字开关控制信号进行通断处理后，得到彼此分离的高频分量和低频分量。
并且，上述表达式中包含了接收信号的幅度和相位信息，因此通过采样两路该信号可以进行一步计算出这两路信号的幅度比和相位差，进而得到待测地层的地层电阻率。
相应地，本发明实施例还提供了用于随钻电磁波电阻率测量的信号处理系统。图3和图4分别示出了该信号处理系统的发射部分和接收部分的结构示意图。
参照图3，晶振(即晶体振荡器)110，用于产生原始的32MHz频率信号。16倍分频器121，对晶振产生的32MHz信号进行16分频得到2000KHz的信号。80倍分频器122，对晶振产生的32MHz信号进行80分频得到400KHz的信号。多路选择器130，用于2000KHz信号和400KHz信号的选择，输出为2000KHz信号和400KHz信号中的一个。功率放大器(即发射放大器)140，用于对多路选择器130的输出信号进行功率放大。发射匹配电容(即阻抗匹配电路)150，用于对功率放大器140的输出进行信号匹配，使得信号有效耦合到发射线圈160。发射线圈160，用于对发射匹配电容150的输出信号发射到地层当中。本实施例中；功率放大器140和发射匹配电容150等器件或芯片的温度特性非常重要，其温度特性直接影响最终发射信号的在不同温度下的发射效率。
参照图4，接收线圈对(包括第一接收天线和第二接收天线)230，用于接收从地层中反射过来的电磁波信号。接收匹配电容240，用于提高接收线圈对的接收效率。低噪声放大器(即接收放大器)250，用于对接收的信号进行低噪声放大。数字开关模块260，在数字控制振荡器220产生数字开关控制信号的控制下对低噪声放大器250的输出进行开关控制。带通滤波器(即滤波模块)270，用于对数字开关模块260的输出进行滤波，滤除不需要的信号。ADC采样器280，用于对带通滤波器270的输出进行数字化采样。微处理器(即计算模块)290，用于对ADC采样器280输出的数字信号进行计算处理，计算出两路接收信号的相位差和幅度比，再对照地层电阻率图板，反演得到地层电阻率。
本实施例还中包括：数字控制振荡器220，用于产生2002KHz和402KHz的数字控制信号，其输出连接到数字开关模块260。
值得说明的是，数字开关模块260、数字控制振荡器220、带通滤波器270可集成在同一个微处理器290中，也可以分别由不同的处理器进行处理。
在本实施例中，数字开关模块的输入为第一/第二放大信号，数字开关控制信号F(t)满足：
其中：
并且，t表示时间，f0表示第一频率，f1表示第二频率，第一频率大于第二频率，并且第一/第二放大信号中频率大于第一频率的分量为高频分量，频率小于或者等于第一频率的分量为低频分量。
上述各模块中的操作的具体细化，可参见上面结合图1和图2对本发明方法的说明，在此不再详细赘述。
应用本发明实施例提供的信号处理系统，采用数字开关控制通断的方式实现信号中高频分量和低频分量的有效分离，结合滤波器，可实现低频分量的有效提取(即高频分量向低频分量的转换)。与现有技术中采用混频器来实现高频分量向低频分量的转换的技术方案相比，本实施例的数字开关控制模块与滤波模块所需的计算量小，因此采用普通的处理器即可满足数字开关控制模块和滤波模块的计算需求。可以看出，采用本实施例所述方法的系统的功耗低，电路简单，可靠性高。
本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (8)
1.一种用于随钻电磁波电阻率测量的信号处理方法，其特征在于，包括：
将预设频率的发射信号发射到待测地层；
获取经待测地层反射回来的第一电磁波信号和第二电磁波信号；
对所述第一电磁波信号和所述第二电磁波信号进行放大，得到第一放大信号和第二放大信号；
利用数字开关模块来分离所述第一放大信号和所述第二放大信号中的高、低频分量，得到对应所述第一放大信号的第一高频分量和第一低频分量，以及对应所述第二放大信号的第二高频分量和第二低频分量，所述数字开关模块由数字开关控制信号控制；
过滤掉所述第一高频分量和所述第二高频分量，以提取出所述第一低频分量和所述第二低频分量；
计算所述第一低频分量和所述第二低频分量的相位差和幅度比；
根据所述相位差和幅度比，得到所述待测地层的电阻率。
2.根据权利要求1所述的方法，其特征在于，所述数字开关模块的输入为第一/第二放大信号，所述数字开关控制信号F(t)满足：
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中：
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
并且，t表示时间，f0表示第一频率，f1表示第二频率，所述第一频率大于所述第二频率，并且第一/第二放大信号中频率大于所述第一频率的分量为所述高频分量，频率小于或者等于所述第一频率的分量为所述低频分量。
3.根据权利要求1或2所述的方法，其特征在于，将预设频率的发射信号发射到待测地层，包括：
生成预设频率的发射信号；
对所述发射信号进行放大；
对放大后的发射信号进行阻抗匹配；
由发射线圈将经阻抗匹配的信号发射到所述待测地层。
4.根据权利要求1至3中任一项所述的方法，其特征在于，所述预设频率为2000KHz或者400KHz。
5.一种用于随钻电磁波电阻率测量的信号处理系统，其特征在于，包括：
发射器，设置为将预设频率的发射信号发射到待测地层；
第一接收天线，设置为获取经待测地层反射回来的第一电磁波信号；
第二接收天线，设置为获取经待测地层反射回来的第二电磁波信号；
接收放大器，设置为对所述第一电磁波信号和所述第二电磁波信号进行放大，得到第一放大信号和第二放大信号；
数字开关模块，设置为分离所述第一放大信号和所述第二放大信号中的高、低频分量，得到对应所述第一放大信号的第一高频分量和第一低频分量，以及对应所述第二放大信号的第二高频分量和第二低频分量，所述数字开关由数字开关控制信号控制；
数字控制振荡器，设置为产生所述数字开关控制信号；
滤波模块，设置为过滤掉所述第一高频分量和所述第二高频分量，以提取出所述第一低频分量和所述第二低频分量；
计算模块，设置为计算所述第一低频分量和所述第二低频分量的相位差和幅度比；并根据所述相位差和幅度比，得到所述待测地层的电阻率。
6.根据权利要求5所述的系统，其特征在于，所述数字开关模块的输入为第一/第二放大信号，所述数字开关控制信号F(t)满足：
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>n</mi>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中：
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
并且，t表示时间，f0表示第一频率，f1表示第二频率，所述第一频率大于所述第二频率，并且第一/第二放大信号中频率大于所述第一频率的分量为所述高频分量，频率小于或者等于所述第一频率的分量为所述低频分量。
7.根据权利要求5或6所述的系统，其特征在于，所述发射器包括：
晶体振荡器，设置为生成预设频率的发射信号；
发射放大器，设置为对所述发射信号进行放大；
阻抗匹配电路，设置为对放大后的发射信号进行阻抗匹配；
发射线圈，设置为将经阻抗匹配的信号发射到所述待测地层。
8.根据权利要求5至7中任一项所述的系统，其特征在于，所述预设频率为2000KHz或者400KHz。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610118358.3A CN107152276B (zh) | 2016-03-02 | 2016-03-02 | 用于随钻电磁波电阻率测量的信号处理方法与系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610118358.3A CN107152276B (zh) | 2016-03-02 | 2016-03-02 | 用于随钻电磁波电阻率测量的信号处理方法与系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107152276A true CN107152276A (zh) | 2017-09-12 |
| CN107152276B CN107152276B (zh) | 2020-06-02 |
Family
ID=59792020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610118358.3A Active CN107152276B (zh) | 2016-03-02 | 2016-03-02 | 用于随钻电磁波电阻率测量的信号处理方法与系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107152276B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109838229A (zh) * | 2017-11-27 | 2019-06-04 | 中石化石油工程技术服务有限公司 | 一种电磁波电阻率数据处理方法 |
| CN113109878A (zh) * | 2021-04-09 | 2021-07-13 | 电子科技大学 | 一种用于井间电磁远探测收发天线的优化系统及方法 |
| CN113138425A (zh) * | 2020-01-19 | 2021-07-20 | 中国石油化工股份有限公司 | 一种随钻测井电磁波数据采集方法及装置 |
| CN113593213A (zh) * | 2021-09-05 | 2021-11-02 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种适用于钻孔无线遥测系统的地面发射装置及方法 |
| CN115685355A (zh) * | 2022-11-01 | 2023-02-03 | 中国科学院地质与地球物理研究所 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
| CN115903588A (zh) * | 2022-11-01 | 2023-04-04 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器的信号采集方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4916616A (en) * | 1986-12-08 | 1990-04-10 | Bp Exploration, Inc. | Self-consistent log interpretation method |
| CN1589521A (zh) * | 2001-10-01 | 2005-03-02 | 飞思卡尔半导体公司 | 双控制频率合成器 |
| CN1687816A (zh) * | 2005-06-17 | 2005-10-26 | 哈尔滨工业大学 | 超精密回转轴与激光直写机直写光轴空间对准方法 |
| CN101049013A (zh) * | 2004-10-22 | 2007-10-03 | 三洋电机株式会社 | 噪声去除装置 |
| CN201004168Y (zh) * | 2007-01-29 | 2008-01-09 | 大庆油田装备制造集团 | 测井车智能控制系统 |
| CN102062797A (zh) * | 2009-11-17 | 2011-05-18 | 北京普源精电科技有限公司 | 一种具有高频低频路径分离电路的示波器 |
| CN103701434A (zh) * | 2013-12-18 | 2014-04-02 | 贝兹维仪器(苏州)有限公司 | 一种用于随钻电阻率测量的调频装置及方法 |
-
2016
- 2016-03-02 CN CN201610118358.3A patent/CN107152276B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4916616A (en) * | 1986-12-08 | 1990-04-10 | Bp Exploration, Inc. | Self-consistent log interpretation method |
| CN1589521A (zh) * | 2001-10-01 | 2005-03-02 | 飞思卡尔半导体公司 | 双控制频率合成器 |
| CN101049013A (zh) * | 2004-10-22 | 2007-10-03 | 三洋电机株式会社 | 噪声去除装置 |
| CN1687816A (zh) * | 2005-06-17 | 2005-10-26 | 哈尔滨工业大学 | 超精密回转轴与激光直写机直写光轴空间对准方法 |
| CN201004168Y (zh) * | 2007-01-29 | 2008-01-09 | 大庆油田装备制造集团 | 测井车智能控制系统 |
| CN102062797A (zh) * | 2009-11-17 | 2011-05-18 | 北京普源精电科技有限公司 | 一种具有高频低频路径分离电路的示波器 |
| CN103701434A (zh) * | 2013-12-18 | 2014-04-02 | 贝兹维仪器(苏州)有限公司 | 一种用于随钻电阻率测量的调频装置及方法 |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109838229A (zh) * | 2017-11-27 | 2019-06-04 | 中石化石油工程技术服务有限公司 | 一种电磁波电阻率数据处理方法 |
| CN113138425A (zh) * | 2020-01-19 | 2021-07-20 | 中国石油化工股份有限公司 | 一种随钻测井电磁波数据采集方法及装置 |
| CN113138425B (zh) * | 2020-01-19 | 2024-02-27 | 中国石油化工股份有限公司 | 一种随钻测井电磁波数据采集方法及装置 |
| CN113109878A (zh) * | 2021-04-09 | 2021-07-13 | 电子科技大学 | 一种用于井间电磁远探测收发天线的优化系统及方法 |
| CN113593213A (zh) * | 2021-09-05 | 2021-11-02 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | 一种适用于钻孔无线遥测系统的地面发射装置及方法 |
| CN115685355A (zh) * | 2022-11-01 | 2023-02-03 | 中国科学院地质与地球物理研究所 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
| CN115903588A (zh) * | 2022-11-01 | 2023-04-04 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器的信号采集方法及装置 |
| CN115903588B (zh) * | 2022-11-01 | 2023-09-15 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器的信号采集方法及装置 |
| CN115685355B (zh) * | 2022-11-01 | 2024-02-27 | 中国科学院地质与地球物理研究所 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107152276B (zh) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107152276A (zh) | 用于随钻电磁波电阻率测量的信号处理方法与系统 | |
| CN102704921B (zh) | 一种随钻电磁波电阻率的测量装置和测量方法 | |
| CN100386645C (zh) | 用雷达电波探测江河湖泊表面流速的方法及其雷达系统 | |
| CN106154341B (zh) | 一种核磁共振与瞬变电磁一体化探测仪器及工作方法 | |
| CN107765270A (zh) | 一种基于卡尔曼滤波的卫星导航接收机跟踪环 | |
| CN107134996A (zh) | 一种多参考模式的数字锁相放大器 | |
| CN103217578B (zh) | 基于锁相环技术测量信号间相位差的数字相位计及其方法 | |
| CN108055058A (zh) | 载波多普勒及其变化率的高精度测量方法 | |
| CN103840795A (zh) | 一种基于dds芯片移相的正交检波电路 | |
| Wan et al. | A compact magnetic resonance system with inside-out sensor for soil moisture measurement | |
| CN117328862B (zh) | 应用于随钻方位测井仪的全息测井方法和发射接收电路 | |
| CN102928887A (zh) | 一种无线电坑道透视仪一体化信号接收装置 | |
| CN106990372A (zh) | 一种核磁共振射频天线电路及其阻抗匹配方法 | |
| CN104316969B (zh) | 一种智能型数字化氦光泵磁力仪 | |
| Li et al. | Development of a distributed hybrid seismic–electrical data acquisition system based on the Narrowband Internet of Things (NB-IoT) technology | |
| CN203759264U (zh) | 基于工频整周期触发的核磁共振双极性叠加消噪装置 | |
| CN203399086U (zh) | 基于dds谐波提取技术的频率源 | |
| CN205561748U (zh) | 一种电流测距装置 | |
| CN104678388A (zh) | 一种用于提高测高精度的gnss-r正交化交错处理装置 | |
| CN202679345U (zh) | 一种Ku波段超低相噪频率源组件 | |
| CN206710579U (zh) | 拉莫尔磁共振信号的高精度检测装置 | |
| CN205377837U (zh) | 一种锁相频率合成器 | |
| CN104090163A (zh) | 一种高稳定度、高精度的相位幅度检测装置 | |
| CN204731426U (zh) | 一种甚低频探水仪 | |
| CN107121586A (zh) | 一种双锁相技术20Hz~20kHz多频信号幅度相位实时检测分布式系统 |
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 |