CN108071385B - 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 - Google Patents
一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 Download PDFInfo
- Publication number
- CN108071385B CN108071385B CN201611019419.7A CN201611019419A CN108071385B CN 108071385 B CN108071385 B CN 108071385B CN 201611019419 A CN201611019419 A CN 201611019419A CN 108071385 B CN108071385 B CN 108071385B
- Authority
- CN
- China
- Prior art keywords
- sampling
- analog
- digital
- signal
- digital conversion
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005553 drilling Methods 0.000 title claims abstract description 26
- 238000005070 sampling Methods 0.000 claims abstract description 172
- 238000006243 chemical reaction Methods 0.000 claims abstract description 113
- 238000012545 processing Methods 0.000 claims abstract description 57
- 238000009825 accumulation Methods 0.000 claims abstract description 53
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 239000000126 substance Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 6
- 239000003990 capacitor Substances 0.000 description 78
- 230000015572 biosynthetic process Effects 0.000 description 12
- 239000013078 crystal Substances 0.000 description 11
- 239000000919 ceramic Substances 0.000 description 9
- 238000003786 synthesis reaction Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 230000010363 phase shift Effects 0.000 description 6
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 5
- 238000007667 floating Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
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
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analogue/Digital Conversion (AREA)
Abstract
本发明公开了一种随钻电磁波电阻率幅值比和相位差的测量方法和电路。测量方法包括:将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频得到每一路电磁波信号中的高频信号和中频信号;对每一路电磁波信号中的高频信号和中频信号进行中频带通滤波,保留中频信号;对每一路电磁波信号中的中频信号进行模数转换和采集,每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理;根据两路采样数字信号累加运算处理的结果,生成两路采样数字信号的幅值比和相位差。采用本发明可以降低随钻电磁波电阻率测量***中的幅值比和相位差测量与计算处理的复杂程度。
Description
技术领域
本发明涉及油气开发与勘探技术领域,尤其涉及一种随钻电磁波电阻率幅值比和相位差的测量方法和电路。
背景技术
随钻电阻率仪器可以通过测定地层电阻率的方法,指导随钻井眼轨迹尽可能地在存储内钻进,从而有效增大储层的泄油面积,提高采收率。目前使用的随钻电阻率测量工具以随钻电磁波电阻率测量应用最为广泛。由于电磁波在穿越地层时产生幅度衰减和相位偏移,并且产生幅度衰减和相位偏移由地层的电阻率和介电常数的特性所决定。电磁波在穿越不同的地层介质时产生的幅度衰减和相位偏移不同,不同的频率的电磁波对地层介质的响应特征也不同。10MHz以上频率的电磁波,其产生的幅度衰减和相位偏移主要与地层的介电常数相关;而10MHz以下的电磁波,其产生的电磁波幅度衰减和相位偏移主要与地层的电阻率相关。
随钻电磁波电阻率测量技术需要通过数字信号处理器、现场可编程门阵列/复杂可编程逻辑器件,以及大容量的随机存储器,进行电磁波在穿越地层时产生幅度衰减和相位偏移的测量和计算。该技术方案首先对接收的电磁波信号进行高速采样,再利用数字信号处理器的高速浮点运算、现场可编程门阵列/复杂可编程逻辑器件的高速逻辑运算能力,以及大容量随机存储器的存储能力进行大数据的存储与运算。具体的,首先对两个接收天线线圈接收到的两路电磁波信号进行大量地采样,然后将采样数据存储到大容量随机存储器中。在采集完大量的数据后将它们提供给数字信号处理器、现场可编程门阵列/复杂可编程逻辑器件进行复杂的浮点运算。最后计算出两个接收天线线圈所接收的两路电磁波信号的幅度衰减和相位偏移。
现有技术的不足在于:由于数字信号处理器、现场可编程门阵列/复杂可编程逻辑器件,以及大容量的随机存储器这些芯片的功耗很高,电路设计起来也极为复杂,对随钻测量***这样需要电池供电,电路板安装空间有限,工作在高温、高振动的环境是非常不利的。
发明内容
针对上述技术问题,本发明提供了一种随钻电磁波电阻率幅值比和相位差的测量方法,包括以下步骤:
将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
对所述每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
对每一路电磁波信号中的中频信号进行模数转换和采集,其中每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理;
根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
一个实施例中,在将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理之前,还包括以下步骤:
对所述接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,生成两路待进行混频处理的电磁波信号。
一个实施例中,对每一路电磁波信号中的中频信号进行模数转换和采集,是将每一路电磁波信号中多个周期的中频信号进行模数转换和采集。
一个实施例中,对每一路电磁波信号中的中频信号进行模数转换和采集,包括:
每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
一个实施例中,每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,包括:
每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
一个实施例中,根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差,包括:
确定第一路采样数字信号和第二路采样数字信号的幅值,分别定义为AMPLITUDE1和AMPLITUDE2,其中,
确定第一路采样数字信号和第二路采样数字信号的相位,分别定义为PHASE1和PHASE2,其中,
根据AMPLITUDE1、AMPLITUDE2、PHASE1和PHASE2确定两路采样数字信号的幅值比和相位差,其中,
两路采样数字信号的幅值比=AMPLITUDE1/AMPLITUDE2,两路采样数字信号的相位差为PHASE1-PHASE2;
其中,表示第一路采样数字信号的N个采样数据累加和的同相分量进行累加处理得到的总和,表示第一路采样数字信号的N个采样数据累加和的正交分量进行累加处理得到的总和,表示第二路采样数字信号的N个采样数据累加和的同相分量进行累加处理得到的总和,表示第二路采样数字信号的N个采样数据累加和的正交分量进行累加处理得到的总和。
根据本发明的另一方面,还提供了一种随钻电磁波电阻率幅值比和相位差的测量电路,包括:
混频器,其用于将两个接收天线接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
滤波器,其与所述混频器电性连接,用于对所述混频器输出的每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
模数转换器,其与所述滤波器电性连接,用于对所述滤波器输出的每一路电磁波信号中的中频信号进行模数转换和采集;
微处理器,其与所述模数转换器电性连接,用于每进行一次模数转换利用该次模数转换与下次模数转换的间隙,对所述模数转换器采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,并根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
一个实施例中,还包括:信号放大器,其用于对两个接收天线接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,并将两路放大后的信号提供给所述混频器。
一个实施例中,所述模数转换器具体用于:每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
一个实施例中,所述微处理器具体用于:
每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
第一,本发明实施例中利用两次模拟数字转换的间隙分别对采集到的两路采样数字信号进行运算,每次采集到的数据就不需要存储,微处理器内部的随机存储器的存储能力就可以满足存储需求了,也即不需要大容量随机存储器存储连续采样的模拟数字转换器的输出信号。这样,可以降低随钻电磁波电阻率幅值比和相位差的测量电路的复杂度,进而降低电路***功耗,提高电路***的寿命和可靠性。
第二,本发明实施例中利用两次模拟数字转换的间隙进行运算,无需进行大数据的计算处理,也即无需增加现场可编程门阵列/复杂可编程逻辑器件这类功耗较高的芯片,优化了电路结构,降低电路***功耗,提高电路***的寿命和可靠性。
第三,本发明实施例中将每一路含有被测地层信息的电磁波信号与设定频率的基准信号进行混频与滤波得到含有被测地层信息的中频信号,该中频信号既保留了被测地层信息,又降低了电磁波信号的频率,可以有效降低数据处理的复杂度。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例共同用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明第一实施例的随钻电磁波电阻率幅值比和相位差的测量方法的流程图;
图2是本发明第一实施例的中频信号Fi的采样方法的流程图;
图3是本发明第二实施例的随钻电磁波电阻率幅值比和相位差的测量电路结构示意图;
图4是本发明第二实施例的模数转换器的电路示意图;
图5是本发明第二实施例的直接数字频率合成器的电路示意图;
图6是本发明第二实施例的微处理器的电路示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下结合附图对本发明作进一步地详细说明。
第一实施例
图1是本发明第一实施例的随钻电磁波电阻率幅值比和相位差的测量方法的流程图。下面结合图1对本方法进行说明。
步骤S110和步骤S120主要介绍如何对接收到的两路含有被测地层信息的电磁波信号进行预处理,具体的预处理过程可以包括混频处理和滤波处理,详细说明如下:
步骤S110,将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
采用远接收天线线圈和近接收天线线圈分别接收含有被测地层信息的电磁波信号,这里含有被测地层信息的电磁波信号为产生幅度衰减和相位偏移的电磁波信号。
在进行混频处理前需要生成设定频率的基准信号,可以通过直接数字频率合成器(DDS)产生,设定频率可以为2.001953MHz,近似于2MHz。
假设接收到的电磁波信号为频率为Ft=Fx+Fi的高频信号,与频率接近2MHz的基准信号进行混频处理后,得到频率为2Fx+Fi的高频信号与频率为Fi的中频信号。
优选的,在将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理之前,还包括以下步骤:
对所述接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,生成两路待进行混频处理的电磁波信号。
具体的,接收到的电磁波信号经过信号放大器放大后得到频率为Ft=Fx+Fi的高频信号,然后将其与频率接近2MHz的基准信号进行混频变换后,分别得到频率为2Fx+Fi的高频信号与频率为Fi的中频信号。
步骤S120,对所述每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
在本步骤中,将每一路频率为2Fx+Fi的高频信号与频率为Fi的中频信号通过中频带通滤波器进行滤波后仅输出频率为Fi=1.953kHz的中频信号。该中频信号既保留了被测地层信息,又降低了电磁波信号的频率,可以有效降低数据处理的复杂度。
对接收到的两路含有被测地层信息的电磁波信号进行预处理后,进行数据采集与运算,详细过程如下:
步骤S130,对每一路电磁波信号中的中频信号进行模数转换和采集,其中每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理。
具体的,采用远接收和近接收的两路模数转换器分别对两路电磁波信号中的中频信号进行模数转换,两路模数转换器采用同时转换启动设计,使得两个模数转换器的采集完全同步。
优选的,对每一路电磁波信号中的中频信号进行模数转换和采集,是对每一路电磁波信号中多个周期的中频信号进行模数转换和采集。
在对中频信号进行模数转换和采样时,为避免外界干扰,需要对多个周期进行模数转换和采样。
优选的,对每一路电磁波信号中的中频信号进行模数转换和采集,包括:
每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
模数转换器的采样频率可以为生成的中频信号的128倍。在实际应用中,也可以选择64倍,256倍等。该处倍数的选择主要考虑以下两个因素:
(1)必须是4的整数倍,这样采样一个周期可以分成4等份。从角度方面理解,相当于一个周期360度,可以等分成4个90度,这样就可以通过一个数组表示正弦和余弦,即sin(90°+α)=cos(α)。
(2)倍数不能小于64,否则会引入太大的采样误差。也不能大于4096,这样会超出模数转换器的采样能力和微控制器的计算能力。
图2为本发明第一实施例的中频信号Fi的采样方法的流程图。参照图2,可以包括以下步骤:
步骤S210,输入Fi中频信号到16位模数转换器;
步骤S220,模数转换器以250KHz的采样频率对Fi中频信号进行采样;
步骤S230,判断是否已采集完成256个周期,若是,转入步骤S240,若否,转入步骤S220;
S240,结束。
由于信号Fi=1.953kHz,模数转换器的采样频率为250kHz,因此对于Fi信号进行采样时,Fi信号的每个周期内的采样点数为:250kHz/1.953kHz=128.0082,接近128。模数转换器每完成一次对Fi信号的采样时共得到32768个采样数据,因此需对Fi信号的256(32768/128)个周期进行采样。
选取的是采样数是128,本领域技术人员还可以在不付出创造性劳动的前提下设计出采样数是64、256、512等的计算***。
下面对于步骤S130中每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理(边转换采集边运算处理)的具体实施过程进行说明。
首先,利用两次模拟数字转换的间隙分别对采集到的两路采样数字信号进行运算,每次采集到的数据就不需要存储,微处理器内部的随机存储器的存储能力就可以满足存储需求了,也即不需要大容量随机存储器存储连续采样的模拟数字转换器的输出信号。这样,可以降低随钻电磁波电阻率幅值比和相位差的测量电路的复杂度,进而降低电路***功耗,提高电路***的寿命和可靠性。
其次,利用两次模拟数字转换的间隙进行运算,每采集一个数据进行一次运算,无需进行大量数据的计算处理,也即无需增加现场可编程门阵列/复杂可编程逻辑器件这类功耗较高的芯片,优化了电路结构,降低电路***功耗,提高电路***的寿命和可靠性。
优选的,每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,包括:
(1)每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
在(1)中,由于Fi信号的每个周期内采样点数都为128,因此为降低噪声的影响,将每个Fi信号周期(共256个周期)内第i(i=1,2,3…128)个采样值分别进行累加。将32768个采样点分成128组,每组分别定义为一个“bin”,并将每组的数据进行累加,得到128个采样数据累加和。
其中,bini表示信号Fi每个周期内第i个采样点数据的累加和。由于每个周期采样128个点,因此每个“bin”的相位相差360/128=2.8°。
(2)每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
由于将所有采样点分成128组,因此每32组“bin”相差π/2。根据正交采样方法分别定义“X”信号和“R”信号。“X”信号和“R”信号分别为采样数字信号的同相分量和正交分量。
下面对正交采样的原理进行说明。
假设信号S(t)=Acos(ωt+φ),分别求得t=0和t=T/4时的信号。
基于上述分析对(1)和(2)的具体实施过程进行说明。
“X”信号为:R_SIN[n]*bini,n=0,1,2...127;
“R”信号为:R_SIN[(n+32)\128]*bini,n=0,1,2...127。
R_SIN[n]=INT(32768*sin(n*2π/128)),n=0,1,2...127。通过乘以32768,然后取整数部分,将正弦或余弦的浮点数转换成整型数。取整可以减少后续的计算量。
1)初始化1个整型数组R_SIN[128]。初始化4个长整型数,XF0=0,RF0=0,XN0=0,RN0=0。
2)同步启动远接收和近接收的两路模数转换器,分别采集第1次模数转换后的远接收和近接收的两路采样数字信号第1个周期内的第1个采样点。每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点。进行128次模数转换后,分别采集到远接收和近接收的两路采样数字信号第1个周期内的128个数据,并将这128个数据依次存放在微处理器内部的内存中。
在分别采集两路采样数字信号的第2个周期的采样点时,每采集一个采样点将该采样点与第1个周期的该位置的采样点进行累加,这样,采集到的数据就不用接着存储,累加后的数据继续存储到微处理器内部的内存中。累加运算是在模数转换的间隙进行的。例如,采集第2个周期的第3个采样点,利用该次模数转换与下次模数转换的间隙将第2个周期的第3个采样点与第1个周期的第3个采样点的数据进行累加。
在进行第3个周期采样时,每采集一个采样点将该采样点的数据和内存相应位置的数据进行累加,并将累加后的结果存放到内存相应位置中。例如,采集第3个周期的第3个采样点,利用该次模数转换与下次模数转换的间隙将第3个周期的第3个采样点的数据与内存存放的第3个数据进行累加,再将结果存储到该位置。
依次重复至第255个周期。
3)在进行第256个周期采样时,将第256个周期的第1个采样点与内存存放的第1个数据相加得到AF1、AN1,分别乘以R_SIN[0]和R_SIN[32],计算得到:
XF1=XF0+AF1*R_SIN[0],RF1=RF0+AF1*R_SIN[32],
XN1=XN0+AN1*R_SIN[0],RN1=RN0+AF1*R_SIN[32]。
将第256个周期的第2个采样点与内存存放的第2个数据相加得到AF2、AN2,分别乘以R_SIN[1]和R_SIN[33],计算得到:
XF2=XF1+AF2*R_SIN[1],RF2=RF1+AF2*R_SIN[33],
XN2=XN1+AN2*R_SIN[1],RN2=RN1+AN2*R_SIN[33]。
以上相乘和累加的运算都是在两次模数转换的间隙完成的。
依次重复,将第256个周期的第95个采样点与内存存放的第95个数据相加得到AF95,AN95,分别乘以R_SIN[95]和R_SIN[0],再依次重复到128次。
最后:计算得到:
XF128=XF127+AF128*R_SIN[127],RF128=RF127+AF128*R_SIN[31],
XN128=XN127+AN128*R_SIN[127],RN128=RN127+AN128*R_SIN[31]。
步骤S140,根据所述两路采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
步骤S140的具体实施过程如下:
确定第一路采样数字信号和第二路采样数字信号的幅值,分别定义为AMPLITUDE1和AMPLITUDE2,其中,
确定第一路采样数字信号和第二路采样数字信号的相位,分别定义为PHASE1和PHASE2,其中,
根据AMPLITUDE1、AMPLITUDE2、PHASE1和PHASE2确定两路采样数字信号的幅值比和相位差,其中,
两路采样数字信号的幅值比=AMPLITUDE1/AMPLITUDE2,两路采样数字信号的相位差为PHASE1-PHASE2;
其中,表示第一路采样数字信号的N个采样数据累加和的同相分量进行累加处理得到的总和,表示第一路采样数字信号的N个采样数据累加和的正交分量进行累加处理得到的总和,表示第二路采样数字信号的N个采样数据累加和的同相分量进行累加处理得到的总和,表示第二路采样数字信号的N个采样数据累加和的正交分量进行累加处理得到的总和。
具体的,根据上述步骤步骤S130中得到的XF128、RF128、XN128、RN128确定两路采样数字信号的幅值比和相位差。
首先,计算远接收的幅值和相位:
AMPLITUDEF=SQRT(XF128 2+RF128 2),PHASEF=ARCTAN2(XF128,RF128);
再次,计算近接收的幅值和相位:
AMPLITUDEN=SQRT(XN128 2+RN128 2),PHASEN=ARCTAN2(XN128,RN128);
最后,计算出两个接收天线的幅值比和相位差:
幅值比=AMPLITUDEF/AMPLITUDEN,相位差=PHASEF-PHASEN。
以上步骤S130和S140中的计算过程既可以在井下的随钻电磁波电阻率仪器的微处理器中实现,也可以把XF、RF、XN、RN通过随钻的传输***传输到地面,在地面的计算机上计算获得。
由此可知,本发明实施例提供的随钻电磁波电阻率幅值比和相位差的测量方法,可以降低随钻电磁波电阻率测量***中的幅值比与相位差测量与计算处理的复杂程度,最终降低整个电磁波电阻率***的复杂程度,降低***功耗,提高***的寿命和可靠性。
综上所述,本实施例的随钻电磁波电阻率幅值比和相位差的测量方法,在油气开发与勘探中具有实际的指导意义。
第二实施例
图3为本发明第二实施例的随钻电磁波电阻率幅值比和相位差的测量电路结构示意图。参照图3,该电路可以包括混频器310、滤波器320、模数转换器330、微处理器340,其中,
混频器,其用于将两个接收天线接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
滤波器,其与所述混频器电性连接,用于对所述混频器输出的每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
模数转换器,其与所述滤波器电性连接,用于对所述滤波器输出的每一路电磁波信号中的中频信号进行模数转换和采集;
微处理器,其与所述模数转换器电性连接,用于每进行一次模数转换利用该次模数转换与下次模数转换的间隙,对所述模数转换器采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,并根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
优选的,还包括:信号放大器350,其用于对两个接收天线接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,并将两路放大后的信号提供给所述混频器。
优选的,所述模数转换器具体用于:每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
优选的,所述微处理器具体用于:
每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
图3中的各个电路的具体设计可以如下所示:
混频器310和信号放大器350的电路设计是现有技术,此处不再赘述。滤波器320采用中频带通滤波器。
图4为本发明第二实施例的模数转换器的电路示意图。参照图4,模数转换器包括第一电容C14、第二电容C15、第三电容C16、第四电容C17、第五电容C18、第六电容C19、第七电容C20、第八电容C21、第一模数转换芯片U3、第二模数转换芯片U4、第一基准信号产生芯片U5,其中,
第一电容、第二电容、第四电容至第六电容中任一电容一端接+5V电源,另一端接地;第三电容、第八电容中任一电容一端与第一基准信号产生芯片第6脚连接,另一端接地;第七电容一端与第二模数转换芯片第37脚连接,另一端接地;第一模数转换芯片和/或第二模数转换芯片第1、4至6、8、17、20、30、31、33、36、38、39脚接地,第2、7、18、19脚接+5V电源,第9至16、21至29、32、34、35脚分别与主控制芯片第15至18、32、43、44、46、48、49、60、61、62、63、2、3、30、29、53、54脚连接,第4、41、42脚与输入口连接,第37、40脚与第一基准信号产生芯片第6脚连接;第一基准信号产生芯片第2脚接+5V电源,第4脚接地;
具体的,第一电容至第八电容采用0.1μF电容、第一模数转换芯片和第二模数转换芯片采用AD7665模数转换芯片,第一基准信号产生芯片采用LT1019ACN8-2.5基准信号产生芯片。
模数转换电路的作用是将输入的模拟信号转换为数字信号,其中:第一电容至第八电容的作用是对信号滤波,第一模数转换芯片和第二模数转换芯片的作用是将模拟信号转换为数字信号,第一基准信号产生芯片作用是产生2.5V的基准信号。
在图3的电路结构中还可以增加一个直接数字频率合成器,具体的电路设计如下所示:
图5为本发明第二实施例的直接数字频率合成器的电路示意图。参照图5,直接数字频率合成器包括第一电阻R5、第二电阻R6、第三电阻R7、第八电容C10、第九电容C11、第十电容C12、第十一电容C13、第一晶振Y2、第一数字频率合成芯片U2,其中,
第一电阻至第三电阻的一端分别与第一数字频率合成芯片第12、21、20脚连接,另一端接地;第八电容、第十一电容一端分别接+5V电源,另一端接地;第九电容、第十电容一端分别与第一数字频率合成芯片第21、20脚连接,另一端接地;第一晶振第1脚接+5V电源,第3脚接地,第4脚与第一数字频率合成芯片第9脚连接;第一数字频率合成芯片第5、10、19、24脚接地,第6、11、18脚接+5V电源,第7、8、22脚分别与主控制芯片第7、13、14脚连接;
具体的,直接数字频率合成电路第一电阻采用11KΩ电阻,第二电阻采用100Ω电阻,第三电阻采用100Ω电阻,第八电容采用0.1μF电容,第九电容采用220pF电容,第十电容采用220pF电容,第十一电容采用0.1μF电容,第一晶振采用32MHz晶振,第一数字频率合成芯片采用AD9850数字频率合成芯片。
直接数字频率合成电路作用是产生2.001953MHz的频率信号,其中:第一电阻作用是限制电流,第二电阻作用是将电流信号转换为电压信号,第三电阻作用是将电流信号转换为电压信号,第八电容作用是对信号进行滤波,第九电容作用是对信号进行滤波,第十电容作用是对信号进行滤波,第十一电容作用是对信号进行滤波,第一晶振作用是产生32MHz频率的信号,第一数字频率合成芯片作用是合成2.001953MHz频率的数字信号。
图6为本发明第二实施例的微处理器的电路示意图。参照图6,微处理器包括第一限流电阻R1、第二限流电阻R2、第一电源滤波电容C1、第二电源滤波电容C2、第三电源滤波电容C3、第四电源滤波电容C4、第五电源滤波电容C5、第六电源滤波电容C6、第七电源滤波电容C7、第八电源滤波电容C8、第九电源滤波电容C9、第一晶振Y1、第一插件ICSP、微控制器U1,其中,
第一限流电阻一端和微控制器第7脚连接,另一端与与第一电源滤波电容一端相连接;第二限流电阻一端与第一电源滤波电容一端相连接,另一端接5V电源,第一电源滤波电容另一端接地;第二电源滤波电容至第四电源滤波电容一端分别与微控制器第56、39、40脚连接,另一端接地;第五电源滤波电容至第九电源滤波电容一端与5V电源连接,另一端接地;第一晶振一端与微控制器第39脚连接,另一端与微控制器第40脚连接;第一插件第1、4、5脚分别与微控制器第7、15、16脚连接,第2脚与5V电源连接,第3脚接地;微控制器第10、19、26、38、57脚与5V电源连接,第9、20、25、41脚接地。
具体的,第一限流电阻R1采用620Ω电阻,第二限流电阻R2采用4.7KΩ电阻,第一电源滤波电容采用0.1μF瓷片电容,第二电源滤波电容采用10μF瓷片电容,第三电源滤波电容采用10pF瓷片电容,第四电源滤波电容采用10pF瓷片电容,第五电源滤波电容采用0.1μF瓷片电容,第六电源滤波电容采用0.1μF瓷片电容,第七电源滤波电容采用0.1μF瓷片电容,第八电源滤波电容采用0.1μF瓷片电容,第九电源滤波电容采用0.1μF瓷片电容、第一晶振采用8MHz晶振、第一插件采用单排6针的插件,微控制器采用dsPIC33EP512GM706处理器。
微处理器用来控制电路板其他电路的工作情况,其中的第一限流电阻、第二限流电阻用来限制电源电流,第一电源滤波电容至第九电源滤波电容用来对微处理器的电源滤波,第一晶振用来为微处理器提供时钟信号,第一插件用来作为微处理器的在线串行编程接口,微控制器作为微处理器电路的总处理器。
以上所述,仅为本发明的具体实施案例,本发明的保护范围并不局限于此,任何熟悉本技术的技术人员在本发明所述的技术规范内,对本发明的修改或替换,都应在本发明的保护范围之内。
Claims (10)
1.一种随钻电磁波电阻率幅值比和相位差的测量方法,包括以下步骤:
将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
对所述每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
对每一路电磁波信号中的中频信号进行模数转换和采集,其中每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理;
根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
2.根据权利要求1所述的方法,其特征在于,在将接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理之前,还包括以下步骤:
对所述接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,生成两路待进行混频处理的电磁波信号。
3.根据权利要求1所述的方法,其特征在于,对每一路电磁波信号中的中频信号进行模数转换和采集,是对每一路电磁波信号中多个周期的中频信号进行模数转换和采集。
4.根据权利要求1所述的方法,其特征在于,对每一路电磁波信号中的中频信号进行模数转换和采集,包括:
每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
5.根据权利要求1所述的方法,其特征在于,每进行一次模数转换利用该次模数转换与下次模数转换的间隙对采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,包括:
每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
6.根据权利要求1所述的方法,其特征在于,根据所述两路采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差,包括:
确定第一路采样数字信号和第二路采样数字信号的幅值,分别定义为AMPLITUDE1和AMPLITUDE2,其中,
确定第一路采样数字信号和第二路采样数字信号的相位,分别定义为PHASE1和PHASE2,其中,
根据AMPLITUDE1、AMPLITUDE2、PHASE1和PHASE2确定两路采样数字信号的幅值比和相位差,其中,
两路采样数字信号的幅值比=AMPLITUDE1/AMPLITUDE2,两路采样数字信号的相位差为PHASE1-PHASE2;
7.一种随钻电磁波电阻率幅值比和相位差的测量电路,其特征在于,包括:
混频器,其用于将两个接收天线接收到的两路含有被测地层信息的电磁波信号分别与一设定频率的基准信号进行混频处理,得到每一路电磁波信号中的高频信号和中频信号;
滤波器,其与所述混频器电性连接,用于对所述混频器输出的每一路电磁波信号中的高频信号和中频信号进行中频带通滤波处理,保留所述中频信号;
模数转换器,其与所述滤波器电性连接,用于对所述滤波器输出的每一路电磁波信号中的中频信号进行模数转换和采集;
微处理器,其与所述模数转换器电性连接,用于每进行一次模数转换利用该次模数转换与下次模数转换的间隙,对所述模数转换器采集到的该次模数转换后的每一路采样数字信号进行累加运算处理,并根据两路所述采样数字信号累加运算处理的结果,生成所述两路采样数字信号的幅值比和相位差。
8.根据权利要求7所述的测量电路,其特征在于,还包括:信号放大器,其用于对两个接收天线接收到的两路含有被测地层信息的电磁波信号分别进行放大处理,并将两路放大后的信号提供给所述混频器。
9.根据权利要求7所述的测量电路,其特征在于,所述模数转换器具体用于:每进行一次模数转换,采集该次模数转换后的每一路采样数字信号一个周期内的一个采样点,每个周期内有N个采样点。
10.根据权利要求7所述的测量电路,其特征在于,所述微处理器具体用于:
每采集一个采样点利用该次模数转换与下次模数转换的间隙,将当前周期内的该采样点与前各周期内对应位置的采样点的数据进行累加处理,累加处理完成后得到N个采样数据累加和;
每完成一个采样数据累加和的累加处理,将每一路采样数字信号的该采样数据累加和与之前各个采样数据累加和的同相分量和正交分量分别进行累加处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611019419.7A CN108071385B (zh) | 2016-11-17 | 2016-11-17 | 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611019419.7A CN108071385B (zh) | 2016-11-17 | 2016-11-17 | 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108071385A CN108071385A (zh) | 2018-05-25 |
CN108071385B true CN108071385B (zh) | 2021-03-16 |
Family
ID=62160480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611019419.7A Active CN108071385B (zh) | 2016-11-17 | 2016-11-17 | 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108071385B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111155982B (zh) * | 2020-01-03 | 2023-01-20 | 电子科技大学 | 一种多频电阻率测量方法 |
CN113138425B (zh) * | 2020-01-19 | 2024-02-27 | 中国石油化工股份有限公司 | 一种随钻测井电磁波数据采集方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN2858950Y (zh) * | 2005-05-27 | 2007-01-17 | 中国石化集团胜利石油管理局钻井工艺研究院 | 随钻双感应电阻率测量仪 |
WO2010083388A2 (en) * | 2009-01-19 | 2010-07-22 | Services Petroliers Schlumberger | Estimating petrophysical parameters and invasion profile using joint induction and pressure data inversion approach |
CN102590869A (zh) * | 2012-03-19 | 2012-07-18 | 中国科学院电工研究所 | 一种人工场源频率域电法勘探方法及勘探*** |
CN102704921A (zh) * | 2012-05-28 | 2012-10-03 | 中国石油天然气集团公司 | 一种随钻电磁波电阻率的测量装置和测量方法 |
CN103297372A (zh) * | 2013-05-22 | 2013-09-11 | 西南石油大学 | 一种井下数据高速无线传输***及其降低峰均比方法 |
CN105089646A (zh) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | 一种集成有数据传输功能的随钻电阻率测量装置及方法 |
-
2016
- 2016-11-17 CN CN201611019419.7A patent/CN108071385B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN2858950Y (zh) * | 2005-05-27 | 2007-01-17 | 中国石化集团胜利石油管理局钻井工艺研究院 | 随钻双感应电阻率测量仪 |
WO2010083388A2 (en) * | 2009-01-19 | 2010-07-22 | Services Petroliers Schlumberger | Estimating petrophysical parameters and invasion profile using joint induction and pressure data inversion approach |
CN102590869A (zh) * | 2012-03-19 | 2012-07-18 | 中国科学院电工研究所 | 一种人工场源频率域电法勘探方法及勘探*** |
CN102704921A (zh) * | 2012-05-28 | 2012-10-03 | 中国石油天然气集团公司 | 一种随钻电磁波电阻率的测量装置和测量方法 |
CN103297372A (zh) * | 2013-05-22 | 2013-09-11 | 西南石油大学 | 一种井下数据高速无线传输***及其降低峰均比方法 |
CN105089646A (zh) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | 一种集成有数据传输功能的随钻电阻率测量装置及方法 |
Non-Patent Citations (1)
Title |
---|
随钻方位电磁波电阻率测量***发展进展;张晓彬等;《测井技术》;20160229;第40卷(第1期);第12-17页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108071385A (zh) | 2018-05-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1130844C (zh) | 用于扩频接收机的方法和结构 | |
CN108071385B (zh) | 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 | |
CN1143448C (zh) | 高功效接收机 | |
CN100353169C (zh) | 测试电子元件的方法及其仪器 | |
CN104898172B (zh) | 一种基于互相关的核磁共振全波信号噪声滤除方法 | |
CN202041640U (zh) | 一种基于gpu的卫星导航软件接收机 | |
CN115685355B (zh) | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 | |
CN86105364A (zh) | 数字球定位接收机中多普勒搜索的方法及装置 | |
CN1545786A (zh) | 监控上行流频带 | |
CN102364420A (zh) | 高精度ect智能传感器 | |
CN105136276A (zh) | 一种非接触式的球磨机旋转筒体振动信号采集*** | |
CN203133273U (zh) | 基于cpci总线的高频地波雷达数据采集和处理装置 | |
CN102801459A (zh) | 多模卫星信号接收装置及其方法 | |
CN106323447A (zh) | 一种基于手机的便携式激光测振仪及其方法 | |
CN103888167A (zh) | 信号捕获方法、装置及带有该装置的设备 | |
EP0648357A1 (en) | HIGH SPEED PROCESSOR AND METHOD, SIGNAL ANALYSIS SYSTEM AND MEASURING DEVICE AND METHOD. | |
CN106324343B (zh) | 基于频移集合经验模态分解的谐波检测方法及检测*** | |
CN111796305A (zh) | 一种高动态导航卫星信号接收方法及接收机 | |
CN101789833B (zh) | 抖动传输特性测量装置和方法 | |
CN114576568B (zh) | 基于次声波的管道泄漏检测方法及装置 | |
CN115963185A (zh) | 基于fpga控制的超声信号采集处理***和方法 | |
CN115903588A (zh) | 一种随钻方位电磁波电阻率仪器的信号采集方法及装置 | |
CN108521388B (zh) | 一种基于tc-ofdm的频率捕获方法、装置、电子设备及存储介质 | |
CN109633707A (zh) | 一种基于预平均处理的变系数匹配滤波的快速捕获方法 | |
Krishnan et al. | Hilbert–Huang transform and its variants in engineering data analytics: State of the art and research challenges |
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 |